Seismic wavefield simulation in 2D elastic and viscoelastic tilted transversely isotropic media: comparisons between four different kinds of finite-difference grid schemes

2013 ◽  
Vol 10 (4) ◽  
pp. 045005 ◽  
Author(s):  
Zhong-sheng Li ◽  
Chao-ying Bai ◽  
Yao-chong Sun
Geophysics ◽  
2010 ◽  
Vol 75 (2) ◽  
pp. S51-S59 ◽  
Author(s):  
Daniela Amazonas ◽  
Rafael Aleixo ◽  
Jörg Schleicher ◽  
Jessé C. Costa

Standard real-valued finite-difference (FD) and Fourier finite-difference (FFD) migrations cannot handle evanescent waves correctly, which can lead to numerical instabilities in the presence of strong velocity variations. A possible solution to these problems is the complex Padé approximation, which avoids problems with evanescent waves by rotating the branch cut of the complex square root. We have applied this approximation to the acoustic wave equation for vertical transversely isotropic media to derive more stable FD and hybrid FD/FFD migrations for such media. Our analysis of the dispersion relation of the new method indicates that it should provide more stable migration results with fewer artifacts and higher accuracy at steep dips. Our studies lead to the conclusion that the rotation angle of the branch cut that should yield the most stable image is 60° for FD migration, as confirmed by numerical impulse responses and work with synthetic data.


2012 ◽  
Vol 433-440 ◽  
pp. 4656-4661
Author(s):  
Qiang Zhang ◽  
Qi Zhen Du ◽  
Xu Fei Gong

We present a staggered-grid finite difference scheme for velocity-stress equations to simulate the elastic wave propagating in transversely isotropic media. Instead of the widely used temporally second-order difference scheme, a temporally fourth-order scheme is obtained in this paper. We approximate the third-order spatial derivatives with 2N-order difference rather than second-order or other fixed order difference as before. Thus, it could be possible to make a balanced accuracy of O (Δt4+Δx2N) with arbitrary N. Related issues such as stability criterion, numerical dispersion, source loading and boundary condition are also discussed in this paper. The numerical modeling result indicates that the scheme is reliable.


Geophysics ◽  
2009 ◽  
Vol 74 (5) ◽  
pp. WB3-WB17 ◽  
Author(s):  
Peter M. Bakker

For the purpose of one-way wave-equation imaging, a pseudoscreen propagator is developed for transversely isotropic media with vertical axes of symmetry. The phase shift for propagation through a depth slice is decomposed into three terms: a Gazdag phase shift for propagation in a laterally homogeneous reference medium, a correction for lateral variability of vertical propagation, and a remaining wide-angle term for oblique directions of propagation. Based on rational function approximation for this remaining wide-angle term, a Fourier finite-difference (FFD) approach with four-way splitting is applied. Fourth-order Padé approximation is unsatisfactory in anelliptic media for large propagation angles with respect to the vertical direction. Therefore, a method of coefficient optimization is developed in conjunction with a method of choosing an adequate homogeneous reference medium in a depth slice. By symmetrizing the finite-difference operators, and because of the choice of the optimized coefficients, the propagator is stable in the sense that the least-squares norm of the wavefield, measured for a frequency-depth slice, does not grow with increasing depth of propagation. A small amount of artificial damping is applied to suppress artifacts that appear at the critical angle defined by the velocities in the reference medium and the actual medium. Synthetic examples confirm that good kinematic accuracy is achieved for a wide range of propagation angles (typically up to 60°).


Sign in / Sign up

Export Citation Format

Share Document