Pure quasi-P wave stable simulation in complex TTI media
Hu Shuhua1, Wang Yuchao1, Liu Wenqing1, Zhou Yang2
1. Northwest Branch, Research Institute of Petroleum Exploration & Development, PetroChina, Lanzhou, Gansu 730020, China; 2. Wave Phenomena and Inversion Imaging Group, School of Ocean and Earth Science, Tongji University, Shanghai 210092, China
Abstract:Facing to more complicated targets, seismic exploration drops on an increasing demand for accurate subsurface imaging.We analyse the dispersion features of the decomposed scalar operator and elliptic differential operator produced by elliptic pseudo differential decomposition method proposed by Xu et al (2015).Then the corresponding first-order differential equation of pure quasi-P wave for tilted transverse isotropic (TTI) media is derived by introducing self-adjoint operators under rotated coordinates to preserve stability.Under the Lebedev staggered grid framework, we deduce the high-order finite difference scheme of the equation and give out the strategy of 2 times calculation of wavefield gradient to preserve the simulating accuracy.We test our algorithm in homogeneous and complex TTI media.The results show that the quasi-P wavefiled is not affected by pseudo-S-waves.Compared with other quasi-P wave simulating methods, our solution produces more stable wavefield in complex anisotropic media which shows that it has a high tolerance to the spatial variation of anisotropic parameters, especially the abrupt changes of the TTI symmetric axis.Our solution can be used for accurate pure quasi-P wave reverse time migration (RTM) to improve the imaging accuracy of subsurface with complex anisotropy.
Baysal E,Kosloff D D,Sherwood J W C.Reverse time migration.Geophysics,1983,48(11):1514-1524.
[2]
刘文卿,王宇超,雍学善等.基于GPU/CPU叠前逆时偏移研究及应用.石油地球物理勘探,2012,47(5):712-716.Liu Wenqing,Wang Yuchao,Yong Xueshan et al.Prestack reverse time migration on GPU/CPU co-parallel computation.OGP,2012,47(5):712-716.
[3]
刘文卿,王西文,刘洪等.盐下构造速度建模与逆时偏移成像研究及应用.地球物理学报,2013,56(2):612-625.Liu Wenqing,Wang Xiwen,Liu Hong et al.Application of velocity modeling and reverse time migration to subsalt structure.Chinese Journal of Geophysics,2013,56(2):612-625.
[4]
刘定进,蒋波,李博等.起伏地表逆时偏移在复杂山前带地震成像中的应用.石油地球物理勘探,2016,51(2):319-324.Liu Dingjin,Jiang Bo,Li Bo et al.Rugged-topography reverse time migration in complex piedmont zone.OGP,2016,51(2):219-324.
[5]
姜国博,曾庆芹,李虹等.逆时偏移技术在土库曼斯坦地区的应用.石油地球物理勘探,2017,52(增刊2):81-85,103.Jiang Guobo,Zeng Qingqin,Li Hong et al.Application of reverse time migration in Turkmenistan area.OGP,2017,52(S2):81-85,103.
[6]
韩令贺,何兵寿.VTI介质一阶准P波方程正演模拟及边界条件.石油地球物理勘探,2010,45(6):819-825.Han Linghe,He Bingshou.First-order Quasi-P wave equation forward modeling and boundary conditions in VTI medium.OGP,2010,45(6):819-825.
[7]
Alkhalifah T.Acoustic approximation for processing in transversely isotropic media.Geophysics,1998,63(2):623-631.
[8]
Alkhalifah T.An acoustic wave equation for aniso-tropic media.Geophysics,2000,65(4):1239-1250.
[9]
Zhou H,Zhang G,Bloor R.An anisotropic acoustic wave equation for modeling in 2D TTI Media.SEG Technical Program Expanded Abstracts,2006,25:194-198.
[10]
Du X,Fletcher R P,Fowler P J.A new pseudo-acoustic wave equation for VTI media.70th EAGE Conference & Exhibition,Extended Abstracts,2008:H033.
[11]
Fowler P J,Du X,Fletcher R P.Coupled equations for reverse time migration in transversely isotropic media.Geophysics,2010,75(1):S11-S22.
[12]
Fletcher R P,Du X,Fowler P J.Reverse time migration in tilted transversely isotropic (TTI) media.Geo-physics,2009,74(6):WCA179-WCA187.
[13]
Duveneck E,Bakker P M.Stable P-wave modeling for reverse time migration in tilted TI media.Geophy-sics,2011,76(2):S65-S75.
[14]
Zhang Y,Zhang H Z,Zhang G Q.A stable TTI re-verse time migration and its implementation.Geophysics,2011,76(3):WA3-WA11.
[15]
Alkhalifah T,Ma X,Waheed U B et al.Efficient anisotropic wavefield extrapolation using effective isotropic models.75th EAGE Conference & Exhibition,2013:Tu-01-16.
[16]
Xu S,Zhou H.Accurate simulations of pure quasi-P-waves in complex anisotropic media.Geophysics,2014,79(6):1-8.
[17]
Xu S,Tang B,Mu J et al.Elliptic decomposition of quasi-P wave equation.77th EAGE Conference & Exhibition,2015:We P612.
[18]
Bube K P,Nemeth T,Stefani J P et al.On the instability in second-order systems for acoustic VTI and TTI media.Geophysics,2012,77(5):T171-T186.
[19]
Hu S H,Wang X W,Sun J Q et al.Stable simulating algorithm of pure quasi-P wavefield in complex anisotropic media.78th EAGE Conference & Exhibition,2016:We SRS304.
[20]
李娜,李振春,黄建平等.Lebedev网格与标准交错网格耦合机制下的复杂各向异性正演模拟.石油地球物理勘探,2014,49(1):126-131.Li Na,Li Zhenchun,Huang Jianping et al.Numerical simulation with coupling Lebedev and standard staggered grid schemes for complex anisotropic media.OGP,2014,49(1):126-131.
[21]
董良国,马在田,曹景忠等.一阶弹性波方程交错网格高阶差分解法.地球物理学报,2000,43(3):411-419.Dong Liangguo,Ma Zaitian,Cao Jingzhong et al.A staggered-grid high-order difference method of one-order elastic wave equation.Chinese Journal of Geophysics,2000,43(3):411-419.