Two-dimensional seismic wave simulation in anisotropic media by non-staggered finite difference method
被引:21
|
作者:
Zhu He-Jun
论文数: 0引用数: 0
h-index: 0
机构:
Peking Univ, Sch Earth & Space Sci, Lab Computat Geodynam, Beijing 100871, Peoples R ChinaPeking Univ, Sch Earth & Space Sci, Lab Computat Geodynam, Beijing 100871, Peoples R China
Zhu He-Jun
[1
]
Zhang Wei
论文数: 0引用数: 0
h-index: 0
机构:
Peking Univ, Sch Earth & Space Sci, Lab Computat Geodynam, Beijing 100871, Peoples R China
Univ Rhode Isl, Grad Sch Oceanog, Kingston, RI 02881 USAPeking Univ, Sch Earth & Space Sci, Lab Computat Geodynam, Beijing 100871, Peoples R China
Zhang Wei
[1
,3
]
Chen Xiao-Fei
论文数: 0引用数: 0
h-index: 0
机构:
Univ Sci & Technol China, Sch Earth & Space Sci, Mengcheng Natl Geophys Observ, Hefei 230026, Peoples R ChinaPeking Univ, Sch Earth & Space Sci, Lab Computat Geodynam, Beijing 100871, Peoples R China
Chen Xiao-Fei
[2
]
机构:
[1] Peking Univ, Sch Earth & Space Sci, Lab Computat Geodynam, Beijing 100871, Peoples R China
[2] Univ Sci & Technol China, Sch Earth & Space Sci, Mengcheng Natl Geophys Observ, Hefei 230026, Peoples R China
[3] Univ Rhode Isl, Grad Sch Oceanog, Kingston, RI 02881 USA
来源:
CHINESE JOURNAL OF GEOPHYSICS-CHINESE EDITION
|
2009年
/
52卷
/
06期
In this paper, a DRP/opt MacCormack finite difference scheme is used to simulate seismic wave propagation in two-dimensional anisotropic media. DRP/opt MacCormack is a non-staggered finite difference scheme which avoids the stress/strain interpolation in traditional staggered finite difference method. Compared to low-order non-staggered finite different scheme, it has low dispersion and dissipation. This scheme divides center difference to forward and backward one-side difference operators which are combined in 4-6 steps Runge-Kutta time integration. In a transverse isotropic media with vertical axis (VTI), we validate the high accuracy and stability of DRP/opt MacCormack scheme by comparing with spectral element method. In the realistic geological situations, the symmetric axes of transverse isotropic media are tiled (TTI), we simulate three components seismic wavefields in two-dimensional TTI media. The results show shear wave splitting and decoupling of in-plane/anti-plane movements. The numerical simulations show that DRP/opt MacCormack scheme is an efficient tool to study the wavefields in arbitrary anisotropic media. This scheme is suitable for media with arbitrary topographic variations and can be easily implemented in three-dimensional simulations.
机构:
Univ Calif Santa Barbara, Inst Crustal Studies, Santa Barbara, CA 93106 USAUniv Calif Santa Barbara, Inst Crustal Studies, Santa Barbara, CA 93106 USA
机构:
Univ Calif Santa Barbara, Inst Crustal Studies, Santa Barbara, CA 93106 USAUniv Calif Santa Barbara, Inst Crustal Studies, Santa Barbara, CA 93106 USA