High-Order Finite-Difference Numerical Modeling of Wave Propagation in Viscoelastic TTI Media Using Rotated Staggered Grid

2012 ◽  
Vol 55 (2) ◽  
pp. 252-265 ◽  
Author(s):  
Hong-Yong YAN ◽  
Yang LIU
Keyword(s):  
Numerical Modeling ◽  
Wave Propagation ◽  
Finite Difference ◽  
High Order ◽  
Staggered Grid ◽  
Tti Media ◽  
High Order Finite Difference

On the construction and efficiency of staggered numerical differentiators for the wave equation

Geophysics ◽  
1990 ◽  
Vol 55 (1) ◽  
pp. 107-110 ◽  
Author(s):  
M. Kindelan ◽  
A. Kamel ◽  
P. Sguazzero
Keyword(s):  
Numerical Modeling ◽  
Wave Propagation ◽  
Wave Equation ◽  
Finite Difference ◽  
Computational Efficiency ◽  
Differential Operators ◽  
High Order

Finite‐difference (FD) techniques have established themselves as viable tools for the numerical modeling of wave propagation. The accuracy and the computational efficiency of numerical modeling can be enhanced by using high‐order spatial differential operators (Dablain,1986).

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.

Analyses of the Dispersion Overshoot and Inverse Dissipation of the High-Order Finite Difference Scheme

2013 ◽  
Vol 5 (06) ◽  
pp. 809-824 ◽  
Author(s):  
Qin Li ◽  
Qilong Guo ◽  
Hanxin Zhang
Keyword(s):  
Finite Difference ◽  
Difference Scheme ◽  
Finite Difference Scheme ◽  
High Order ◽  
Staggered Grid ◽  
Third Order ◽  
Order Scheme ◽  
Relationship Of ◽  
High Order Finite Difference ◽  
The Relationship

AbstractAnalyses were performed on the dispersion overshoot and inverse dissipation of the high-order finite difference scheme using Fourier and precision analysis. Schemes under discussion included the pointwise- and staggered-grid type, and were presented in weighted form using candidate schemes with third-order accuracy and three-point stencil. All of these were commonly used in the construction of difference schemes. Criteria for the dispersion overshoot were presented and their critical states were discussed. Two kinds of instabilities were studied due to inverse dissipation, especially those that occur at lower wave numbers. Criteria for the occurrence were presented and the relationship of the two instabilities was discussed. Comparisons were made between the analytical results and the dispersion/dissipation relations by Fourier transformation of typical schemes. As an example, an application of the criteria was given for the remedy of inverse dissipation in Weirs & Martín’s third-order scheme.

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.

Sign in / Sign up

Export Citation Format

Share Document