2D three-component seismic forward modeling in TTI media based on the Lebedev grid

被引:0
作者
Liu D.-Y. [1 ]
Peng S. [1 ]
Shi S. [1 ]
Zhao T. [1 ]
机构
[1] State Key Laboratory of Coal Resources and Safe Mining, China University of Mining and Technology (Beijing), Beijing
来源
Liu, Dong-Yang (18910219317@163.com) | 2018年 / Science Press卷 / 53期
关键词
Lebedev grid; Lens model; Multi-axis perfectly matched layer (MPML); TTI media;
D O I
10.13810/j.cnki.issn.1000-7210.2018.02.009
中图分类号
学科分类号
摘要
The conventional staggered grid finite difference method is a numerical simulation method which is often used to study seismic wave propagation.However,different components of each variable in the staggered grid is staggered defined.For undefined points,their variables need to be interpolated,which reduces the accuracy of simulation.Therefore,the 2D 3C velocity-stress elastic wave equation of the TTI media is deduced on the basis of previous studies and the Lebedev grid is used to process the high accuracy differential separation.This scheme can avoid errors caused by the wave field interpolation in the conventional staggered grid when dealing with anisotropic media and it improves the modulated precision.Meanwhile,the multi-axis perfectly matched layer (MPML) method is introduced into this method.The forward modeling of the single layer TTI medium and the complex TTI medium with the lens is carried out.The results indicate as follows: A.The Lebedev grid does not need to interpolate the wavefield when the elastic wave equation of the anisotropic medium is discrete.Compared with the conventional staggered grid method,the simulation accuracy is higher.The test shows that P-wave,fast S-wave and slow S-wave can be clearly observed when the TTI medium is simulated and the polarization direction of fast S-wave and slow S-wave is opposite.The velocity characteristics of the three kinds of waves observed in single gun record also conform to the propagation regulation of seismic wavefield; B.After introducing MPML into this method,the phenomenon of boundary reflection is effectively suppressed without affecting the simulation effect. © 2018, Editorial Department OIL GEOPHYSICAL PROSPECTING. All right reserved.
引用
收藏
页码:288 / 296
页数:8
相关论文
共 26 条
[1]  
Tsvankin I., Thomsen L., Inversion of reflection traveltimes for transverse isotropy, Geophysics, 60, 4, pp. 1095-1107, (1995)
[2]  
Alterman Z., Karal F.C., Propagation of elastic waves in layered media by finite difference methods, Bulletin of the Seismological Society of America, 58, 1, pp. 367-398, (1968)
[3]  
Virieux J., SH wave propagation in heterogeneous media: velocity-stress finite-difference method, Geophysics, 49, 11, pp. 1933-1942, (1984)
[4]  
Virieux J., P-SV wave propagation in heterogeneous media: velocity-stress finite-difference method, Geophysics, 51, 4, pp. 889-901, (1986)
[5]  
Dablain M.A., The application of high-order differencing to scalar wave equation, Geophysics, 51, 1, pp. 54-66, (1986)
[6]  
Graves R.W., Simulating seismic wave propagation in 3D elastic media using staggered-grid finite differences, Bulletin of the Seismological Society of America, 86, 4, pp. 1091-1106, (1996)
[7]  
Dong L., Ma Z., Et al., Study on stability of staggered grid high order difference method for first order elastic wave equation, Chinese Journal of Geophysics, 43, 3, pp. 411-419, (2000)
[8]  
Pei Z., Numerical simulation of elastic wave equation with staggered grid high-order finite difference method in three dimensional anisotropic medium, Journal of China University of Petroleum(Edition of Natural Science), 28, 5, pp. 23-29, (2004)
[9]  
Pei Z., Dong Y., Peng S., Numerical simulation of azimuthal anisotropy of elastic wave field in fractured seams, OGP, 42, 6, pp. 665-672, (2007)
[10]  
Pei Z., Numerical simulation of splitting and splitting of transverse wave in layered anisotropic media, OGP, 41, 1, pp. 17-25, (2006)