Finite‐difference modeling of elastic wave propagation: A nonsplitting perfectly matched layer approach

Geophysics ◽  
2003 ◽  
Vol 68 (5) ◽  
pp. 1749-1755 ◽  
Author(s):  
Tsili Wang ◽  
Xiaoming Tang

In this paper, we present a nonsplitting perfectly matched layer (NPML) method for the finite‐difference simulation of elastic wave propagation. Compared to the conventional split‐field approach, the new formulation solves the same set of equations for the boundary and interior regions. The nonsplitting formulation simplifies the perfectly matched layer (PML) algorithm without sacrificing the accuracy of the PML. In addition, the NPML requires nearly the same amount of computer storage as does the split‐field approach. Using the NPML, we calculate dipole and quadrupole waveforms in a logging‐while‐drilling environment. We show that a dipole source produces a strong pipe flexural wave that distorts the formation arrivals of interest. A quadrupole source, however, produces clean formation arrivals. This result indicates that a quadrupole source is more advantageous over a dipole source for shear velocity measurement while drilling.

1964 ◽  
Vol 31 (2) ◽  
pp. 213-222 ◽  
Author(s):  
J. P. Jones

Elastic wave propagation in a medium consisting of two finite layers is considered. Two types of solutions are treated. The first is a Rayleigh train of waves. It is seen that for this case, when the wavelength becomes short, the waves approach two Rayleigh waves plus a possible Stoneley wave. When the wavelength becomes large, there are two waves; i.e., a flexural wave and an axial wave. Calculations are presented for this case. The propagation of SH waves is treated, but no calculations are presented.


Geophysics ◽  
2017 ◽  
Vol 82 (5) ◽  
pp. T207-T224 ◽  
Author(s):  
Zhiming Ren ◽  
Zhen Chun Li

The traditional high-order finite-difference (FD) methods approximate the spatial derivatives to arbitrary even-order accuracy, whereas the time discretization is still of second-order accuracy. Temporal high-order FD methods can improve the accuracy in time greatly. However, the present methods are designed mainly based on the acoustic wave equation instead of elastic approximation. We have developed two temporal high-order staggered-grid FD (SFD) schemes for modeling elastic wave propagation. A new stencil containing the points on the axis and a few off-axial points is introduced to approximate the spatial derivatives. We derive the dispersion relations of the elastic wave equation based on the new stencil, and we estimate FD coefficients by the Taylor series expansion (TE). The TE-based scheme can achieve ([Formula: see text])th-order spatial and ([Formula: see text])th-order temporal accuracy ([Formula: see text]). We further optimize the coefficients of FD operators using a combination of TE and least squares (LS). The FD coefficients at the off-axial and axial points are computed by TE and LS, respectively. To obtain accurate P-, S-, and converted waves, we extend the wavefield decomposition into the temporal high-order SFD schemes. In our modeling, P- and S-wave separation is implemented and P- and S-wavefields are propagated by P- and S-wave dispersion-relation-based FD operators, respectively. We compare our schemes with the conventional SFD method. Numerical examples demonstrate that our TE-based and TE + LS-based schemes have greater accuracy in time and better stability than the conventional method. Moreover, the TE + LS-based scheme is superior to the TE-based scheme in suppressing the spatial dispersion. Owing to the high accuracy in the time and space domains, our new SFD schemes allow for larger time steps and shorter operator lengths, which can improve the computational efficiency.


Sign in / Sign up

Export Citation Format

Share Document