Abstract: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.

The Method Of Finite Difference Of High Precision Elastic Wave Equations In The Frequency Domain And Wave-Field Simulation

An Improved 25‐point Finite‐difference Scheme for Frequency‐domain Elastic Wave Modelling

Numerical simulation of elastic wave based on the staggered grid finite difference method

Quasi-P wave forward modeling in viscoelastic VTI media in frequency-space domain

Analysis of Numerical Dispersion in Wave-Field Simulation

Seismic modeling of scalar seismic wave propagation with finitedifference scheme in frequency-space domain.

2-D Elastic Wave Modeling with Frequency-Space 25-Point Finite-Difference Operators

High Precision Finite Difference Operators for Qp Wave Equation in VTI Media

Optimal finite-difference schemes for elastic wave based on improved cosine-combined window function

A 21-Point Finite Difference Scheme for 2d Frequency-Domain Elastic Wave Modelling

Optimized difference coefficient method seismic wave field numerical simulation in vertical transverse isotropic medium

The Precise Finite Difference Method for Seismic Modeling

High-order Temporal and Implicit Spatial Staggered-Grid Finite-Difference Operators for Modelling Seismic Wave Propagation

A high order finite difference method for the elastic wave equation in bounded domains with nonconforming interfaces

Finite difference method for acoustic-elastic coupled equations of seismic waves in the frequency domain

Simulation of Seismoelectric Waves Using Finite-Difference Frequency-Domain Method: 2-D SHTE Mode

Forward modeling of seismic wave fields in complex anisotropic media using finite difference method

Optimized staggered grid finite-difference numerical simulation of well seismic and analysis of numerical dispersion

Simulate the elastic wavefields in media with an irregular surface topography based on staggered grid finite difference

Multidomain Finite-Difference and Chebyshev Pseudospectral Hybrid Method for Elastic Wave Simulation: Two-Dimensional Case

A Staggered-Grid High-Order Difference Method Of One-Order Elastic Wave Equation