scholarly journals Acoustic wave propagation with new spatial implicit and temporal high-order staggered-grid finite-difference schemes

2021 ◽  
Vol 18 (5) ◽  
pp. 808-823
Author(s):  
Jing Wang ◽  
Yang Liu ◽  
Hongyu Zhou
Keyword(s):  
Finite Difference ◽  
High Order ◽  
Staggered Grid ◽  
Finite Difference Schemes ◽  
Spatial Accuracy ◽  
Time Step ◽  
Dispersion Formula ◽  
Time Space ◽  
L2 Norm ◽  
Spatial Derivatives

Abstract The implicit staggered-grid (SG) finite-difference (FD) method can obtain significant improvement in spatial accuracy for performing numerical simulations of wave equations. Normally, the second-order central grid FD formulas are used to approximate the temporal derivatives, and a relatively fine time step has to be used to reduce the temporal dispersion. To obtain high accuracy both in space and time, we propose a new spatial implicit and temporal high-order SG FD stencil in the time–space domain by incorporating some additional grid points to the conventional implicit FD one. Instead of attaining the implicit FD coefficients by approximating spatial derivatives only, we calculate the coefficients by approximating the temporal and spatial derivatives simultaneously through matching the dispersion formula of the seismic wave equation and compute the FD coefficients of our new stencil by two schemes. The first one is adopting a variable substitution-based Taylor-series expansion (TE) to derive the FD coefficients, which can attain (2M + 2)th-order spatial accuracy and (2N)th-order temporal accuracy. Note that the dispersion formula of our new stencil is non-linear with respect to the axial and off-axial FD coefficients, it is complicated to obtain the optimal spatial and temporal FD coefficients simultaneously. To tackle the issue, we further develop a linear optimisation strategy by minimising the L2-norm errors of the dispersion formula to further improve the accuracy. Dispersion analysis, stability analysis and modelling examples demonstrate the accuracy, stability and efficiency advantages of our two new schemes.

Acoustic and elastic modeling by optimal time-space-domain staggered-grid finite-difference schemes

Geophysics ◽  
2015 ◽  
Vol 80 (1) ◽  
pp. T17-T40 ◽  
Author(s):  
Zhiming Ren ◽  
Yang Liu
Keyword(s):  
Finite Difference ◽  
Stability Condition ◽  
Computational Cost ◽  
Numerical Dispersion ◽  
Staggered Grid ◽  
Equation Modeling ◽  
Optimal Time ◽  
Finite Difference Schemes ◽  
Space Domain ◽  
Time Space

Staggered-grid finite-difference (SFD) methods are widely used in modeling seismic-wave propagation, and the coefficients of finite-difference (FD) operators can be estimated by minimizing dispersion errors using Taylor-series expansion (TE) or optimization. We developed novel optimal time-space-domain SFD schemes for acoustic- and elastic-wave-equation modeling. In our schemes, a fourth-order multiextreme value objective function with respect to FD coefficients was involved. To yield the globally optimal solution with low computational cost, we first used variable substitution to turn our optimization problem into a quadratic convex one and then used least-squares (LS) to derive the optimal SFD coefficients by minimizing the relative error of time-space-domain dispersion relations over a given frequency range. To ensure the robustness of our schemes, a constraint condition was imposed that the dispersion error at each frequency point did not exceed a given threshold. Moreover, the hybrid absorbing boundary condition was applied to remove artificial boundary reflections. We compared our optimal SFD with the conventional, TE-based time-space-domain, and LS-based SFD schemes. Dispersion analysis and numerical simulation results suggested that the new SFD schemes had a smaller numerical dispersion than the other three schemes when the same operator lengths were adopted. In addition, our LS-based time-space-domain SFD can obtain the same modeling accuracy with shorter spatial operator lengths. We also derived the stability condition of our schemes. The experiment results revealed that our new LS-based SFD schemes needed a slightly stricter stability condition.

Temporal high-order staggered-grid finite-difference schemes for elastic wave propagation

Geophysics ◽  
2017 ◽  
Vol 82 (5) ◽  
pp. T207-T224 ◽  
Author(s):  
Zhiming Ren ◽  
Zhen Chun Li
Keyword(s):  
Wave Propagation ◽  
Wave Equation ◽  
Finite Difference ◽  
Elastic Wave ◽  
High Order ◽  
Elastic Wave Propagation ◽  
Staggered Grid ◽  
S Wave ◽  
Order Accuracy ◽  
Spatial Derivatives

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.

Accurate and efficient seismic modeling in random media

Geophysics ◽  
1993 ◽  
Vol 58 (4) ◽  
pp. 576-588 ◽  
Author(s):  
Guido Kneib ◽  
Claudia Kerner
Keyword(s):  
Wave Propagation ◽  
Finite Difference ◽  
Random Media ◽  
Polynomial Interpolation ◽  
Random Medium ◽  
High Order ◽  
Numerical Dispersion ◽  
Staggered Grid ◽  
Finite Difference Schemes ◽  
Seismic Modeling

The optimum method for seismic modeling in random media must (1) be highly accurate to be sensitive to subtle effects of wave propagation, (2) allow coarse sampling to model media that are large compared to the scale lengths and wave propagation distances which are long compared to the wavelengths. This is necessary to obtain statistically meaningful overall attributes of wavefields. High order staggered grid finite‐difference algorithms and the pseudospectral method combine high accuracy in time and space with coarse sampling. Investigations for random media reveal that both methods lead to nearly identical wavefields. The small differences can be attributed mainly to differences in the numerical dispersion. This result is important because it shows that errors of the numerical differentiation which are caused by poor polynomial interpolation near discontinuities do not accumulate but cancel in a random medium where discontinuities are numerous. The differentiator can be longer than the medium scale length. High order staggered grid finite‐difference schemes are more efficient than pseudospectral methods in two‐dimensional (2-D) elastic random media.

Time-space-domain temporal high-order staggered-grid finite-difference schemes by combining orthogonality and pyramid stencils for 3D elastic-wave propagation

Geophysics ◽  
2019 ◽  
Vol 84 (4) ◽  
pp. T259-T282 ◽  
Author(s):  
Shigang Xu ◽  
Yang Liu ◽  
Zhiming Ren ◽  
Hongyu Zhou
Keyword(s):  
Dispersion Relation ◽  
Finite Difference ◽  
Elastic Wave ◽  
Large Time ◽  
High Order ◽  
Staggered Grid ◽  
Elastic Wave Equation ◽  
Space Domain ◽  
First Order ◽  
Time Space

The presently available staggered-grid finite-difference (SGFD) schemes for the 3D first-order elastic-wave equation can only achieve high-order spatial accuracy, but they exhibit second-order temporal accuracy. Therefore, the commonly used SGFD methods may suffer from visible temporal dispersion and even instability when relatively large time steps are involved. To increase the temporal accuracy and stability, we have developed a novel time-space-domain high-order SGFD stencil, characterized by ([Formula: see text])th-order spatial and ([Formula: see text])th-order temporal accuracies ([Formula: see text]), to numerically solve the 3D first-order elastic-wave equation. The core idea of this new stencil is to use a double-pyramid stencil with an operator length parameter [Formula: see text] together with the conventional second-order SGFD to approximate the temporal derivatives. At the same time, the spatial derivatives are discretized by the orthogonality stencil with an operator length parameter [Formula: see text]. We derive the time-space-domain dispersion relation of this new stencil and determine finite-difference (FD) coefficients using the Taylor-series expansion. In addition, we further optimize the spatial FD coefficients by using a least-squares (LS) algorithm to minimize the time-space-domain dispersion relation. To create accurate and reasonable P-, S-, and converted wavefields, we introduce the 3D wavefield-separation technique into our temporal high-order SGFD schemes. The decoupled P- and S-wavefields are extrapolated by using the P- and S-wave dispersion-relation-based FD coefficients, respectively. Moreover, we design an adaptive variable-length operator scheme, including operators [Formula: see text] and [Formula: see text], to reduce the extra computational cost arising from adopting this new stencil. Dispersion and stability analyses indicate that our new methods have higher accuracy and better stability than the conventional ones. Using several 3D modeling examples, we demonstrate that our SGFD schemes can yield greater temporal accuracy on the premise of guaranteeing high-order spatial accuracy. Through effectively combining our new stencil, LS-based optimization, large time step, variable-length operator, and graphic processing unit, the computational efficiency can be significantly improved for the 3D case.

Removing the stability limit of the time–space domain explicit finite difference schemes for acoustic modeling with stability condition-based spatial operators

Geophysics ◽  
2021 ◽  
pp. 1-82
Author(s):  
Yang Liu
Keyword(s):  
Finite Difference ◽  
Stability Condition ◽  
Finite Difference Schemes ◽  
Grid Spacing ◽  
Time Step ◽  
Space Domain ◽  
Wavenumber Range ◽  
Time Space ◽  
Explicit Finite Difference ◽  
The Stability

The time step and grid spacing in explicit finite-difference (FD) modeling are constrained by the Courant-Friedrichs-Lewy (CFL) condition. Recently, it has been found that spatial FD coefficients may be designed through simultaneously minimizing the spatial dispersion error and maximizing the CFL number. This allows one to stably use a larger time step or a smaller grid spacing than usually possible. However, when using such a method, only second-order temporal accuracy is achieved. To address this issue, I propose a method to determine the spatial FD coefficients, which simultaneously satisfy the stability condition of the whole wavenumber range and the time–space domain dispersion relation of a given wavenumber range. Therefore, stable modeling can be performed with high-order spatial and temporal accuracy. The coefficients can adapt to the variation of velocity in heterogeneous models. Additionally, based on the hybrid absorbing boundary condition, I develop a strategy to stably and effectively suppress artificial reflections from the model boundaries for large CFL numbers. Stability analysis, accuracy analysis and numerical modeling demonstrate the accuracy and effectiveness of the proposed method.

Sign in / Sign up

Export Citation Format

Share Document