For simplicity, I use the acoustic wave equation (11), which is second order in t, as opposed to equation (7) [fourth order in derivates of t]. Adding a force function, f(x,y,z,t), to equation (11) yields
In addition, we must solve
for F along with solving equation (20).
Using the second-order finite-difference approach, we can approximate partial derivatives with difference equations, as follows
where i and j are indexes corresponding to steps in position or time. The later equation is derived from the 2-D convolution of the second-order approximation for each of x and z derivatives, and it is conveniently of a fourth-order accuracy.
Therefore, F in equation (21) is calculated using the following recursive formula:
and P is calculated using
where is the finite-difference approximation of the new Laplacian equation (20).
Higher-order finite-difference approximations (4th order) can be used with equation (7) directly. The set of finite-difference equations (20-21) are subjected to the same constraints and rules used in the isotropic case (such as the CFL condition) to avoid numerical dispersion and instability.
Figure 2 shows the wavefield at time 1 second caused by an impulse force excited at time 0. On the left side, the medium is homogeneous and VTI with vv=2 km/s, v=2 km/s, and =0.4. On the right side, the medium is isotropic with vv=2 km/s, v=2 km/s, and =0. Both wavefields are calculated using the second-order finite difference applied to the new acoustic wave equation. In the case of VTI, an additional wave type appears in the section and travels at a speed that is lower than the P-wave velocity. This artifact is the additional solution, mentioned earlier, that behaves like a wave for positive and exponentially decays or grows for negative .Such an artifact does not appear in the solution for an isotropic medium. Therefore, we may want to place the source in an isotropic layer and take advantage of the evanescent nature of this wave in isotropic media. The black curves in Figure 2 correspond to solutions of the eikonal equation. The solid curves correspond to the solutions that use the acoustic assumption, in equation (16), in which the shear wave velocity equals zero, while the dashed curves correspond to a shear wave velocity equal to half the P-wave velocity. As expected, in the isotropic case both curves exactly coincide and are therefore indistinguishable. In the VTI case, differences between the two curves exist, but are hardly noticeable. The independence of the eikonal equation on the shear wave velocity in VTI medium is in agreement with the results I obtained in an earlier study 1997.
Figure 3 The z-component of the elastic wavefield at 1 s caused by a source at time 0 located at the center. The VTI model is the same as in Figure 2, with v=1000 m/s and =0.4. The solid curve corresponds to the solution of the acoustic eikonal equation for the same medium.
Figure 3 shows the z-component of the elastic wavefield (computed using the elastic wave equation) for the same model used in Figure 2. The solid curve is the solution of the acoustic eikonal equation. Kinematically, for P-waves, the acoustic and elastic wavefields are similar. Dynamically, they differ considerably; the elastic wavefield includes S-waves (the slower wave), here with triplication, a phenomenon common to S-waves in strongly anisotropic media.