Compared with VTI qP-wave modeling,TTI qP-wave numerical simulation method which considers dip angle factor could be more accurate in describing the characteristics of wave-field propagation in real anisotropic media.Because of the rapid change of dip angle in TTI media,using the conventional acoustic approximated method may cause numerical instability.In this paper,the authors first derived the stable TTI second-order qP-wave equations by introducing an anisotropic control parameter to ensure the stability of qP-wave propagation.Then,the equivalent first-order stress-velocity forms were deduced through introducing the pseudo-velocity components of the wave-fields,and the corresponding PML boundary condition was given.Finally,the authors implemented accurate numerical simulation using the optimal rotated stagger-grid finite-difference (RSGFD) method.Numerical results demonstrate that TTI first-order qP-wave equations can describe the kinematics features of qP-wave effectively and stably,and the application of optimal RSGFD method can acquire accurate synthetic seismic recordings and enhance computational efficiency relatively.
[1] 李振春.地震叠前成像理论与方法[M].东营:中国石油大学出版社,2011:24-39.[2] 裴正林.三维各向异性介质中弹性波方程交错网格高阶有限差分法数值模拟[J].石油大学学报:自然科学版,2004,28(5):23-29.[3] 何兵寿,张会星.VTI介质中准P波方程的数值解法[J].煤炭学报,2006,31(4):446-450.[4] Alkhalifah T.An acoustic wave equation for anisotropic media[J].Geophysics,2000,65(4):1239-1250.[5] Zhang L B,Rector J W,Hoversten G M.An acoustic wave equation for modeling in tilted TI media[J].SEG,2003.[6] Zhou H B,Zhang G Q,Bloor R.An anisotropic acoustic wave equation for modeling and migration in 2D TTI media[C]//Expanded Abstracts of 76th SEG Annual Internet Meeting,2006:194-198.[7] Fletcher R,Du X,Fowler P J.A new pseudo-acoustic wave equation for TI media[C]//Expanded Abstracts of 78th SEG Annual Internet Meeting,2008:2082-3086.[8] Fowler P J,Du X,Fletcher R P.Coupled equations for reverse time migration in transversely isotropic media[J].Geophysics,2010,75(1):S11-S22.[9] Grechka V,Zhang L,Rector J W.Shear waves in acoustic anisotropic media[J].Geophysics,2004,69(2):576-582.[10] Fletcher R P,Du X,Fowler P J.Reverse time migration in tilted transversely isotropic(TTI) media[J].Geophysics,2009,74(6):WCA179-WCA187.[11] Yoon K,Sang S,Jean J,et al.Stability and speedup issues in TTI RTM implementation[C]//Expanded Abstracts of 80th SEG Annual Internet Meeting,2010:3221-3225.[12] Zhang Y,Zhang H Z,Zhang G Q.A Stable TTI Reverse time migration and its implementation[J].Geophysics,2011,76:WA3-WA11.[13] Zhan G,Pestana R C,Stoffa P L.An acoustic wave equation for pure P wave in 2D TTI media[C]//Expanded Abstracts of 81st SEG Annual Internet Meeting,2011:168-173.[14] Chu C L,Macy B K,Anno P D.An accurate and stable wave equation for pure acoustic TTI modeling[C]//Expanded Abstracts of 81st SEG Annual Internet Meeting,2011:179-184.[15] 程玖兵,康玮,王腾飞.各向异性介质qP波传播描述I:伪纯模式波动方程[J].地球物理学报,2013,56(10):3474-3486.[16] Madariage R.Dynamics of an expanding circular fault[J].Bulletin of the Seismological Society of American,1976,66:639-666.[17] Gold N,Shapiro S A,Burr E,et al.Modeling of high contrasts in elastic media using a modified finite difference scheme[C]//Expanded Abstracts of 67th SEG Annual Internet Meeting,1997:1850-1853.[18] Liu B,Yin X Y,Yang F Y.Optimization of rotated staggered-grid scheme for seismic forward modeling[J].SEG Technical Program Expanded Abstracts,2014:3308-3312.[19] Tsvankin I.Seismic Signatures and Analysis of Reflection Data in Anisotropic Media[J].Elsevier Science,2001:1-60.[20] Thomsen L.Weak elastic anisotropy[J].Geophysics,1986,51(10):1954-1966.[21] Berenger J P.A perfectly matched layer for the absorption of electromagnetic waves[J].Journal of Computational Physics,1994,114(2):185-220.