Global Generalized R/T Recursive Propagation Matrix Method For 2-D Heterogeneous Media Sh Seismic Wave Simulation

Liu Bin,Rong Mianshui,Sun Weijia,Huang Jianping,Ba Zhenning

DOI: https://doi.org/10.6038/cjg2021N0459

2021-01-01

Abstract:The deep earth spheres and sedimentary basins are partitioned heterogeneous medium systems, in which irregular stratigraphic boundaries ( including undulating surface) have a significant influence on the main characteristics of seismic waves, while the random inhomogeneity of the strata mainly affects the scattering and attenuation of seismic waves feature. In order to characterize the seismic wave reflection and transmission affected by irregular stratum boundaries and the seismic wave attenuation caused by heterogeneous body scattering accurately, the global generalized R/T recursive propagation matrix (GGRTM) method was proposed and gradually developed into another semi-analytical solution method for high-precision seismic propagation in complex media after the finite element and finite difference methods. Among the existing methods, the advantage of the GGRTM method for homogeneous medium with irregular boundaries situation lies in the accurate simulation of the reflection and transmission effects of the irregular boundaries, while slicing heterogeneous media into a series of horizontal and heterogeneous slabs and cooperate with GGRTM method can accurately describe the scattering of heterogeneous bodies. Here we absorb the advantages of these two existing methods, and propose a global generalized R/T recursive propagation matrix hybrid method that suits for inhomogeneous media and irregular boundaries, and then use it to simulate the 2-D SH wavefield of a complex model with irregular boundaries and inhomogeneous layers. Subsequently, the simulation accuracy of the method is discussed based on the comparative study of the method in this paper and the boundary element method. The results show that the hybrid method proposed in this paper is an effective way to solve the high-precision seismic simulation of complex models.