基于最佳平方逼近的TTI介质解耦qP波与qSV波逆时偏移

doi:10.13810/j.cnki.issn.1000-7210.2019.06.011

摘要
图/表
参考文献
相关文章 (15)

全文: PDF (6886 KB) HTML []
输出: BibTeX | EndNote (RIS)

摘要在实际地震资料处理中，若忽略各向异性影响、TTI介质对称轴倾角变化，会造成目标体成像位置偏差，甚至由焦散问题造成振幅不均衡。传统TI介质耦合拟声波方程正演模拟与逆时偏移存在伪横波干扰，且在各向异性参数ε<δ时数值模拟不稳定。为此，借鉴前人经验，首先使用最佳平方逼近（OQA）近似VTI介质精确相速度公式，分别得到较高精度的近似qP波、qSV波相速度公式。然后从近似相速度公式出发，推导了二维TTI介质解耦qP波、qSV波波动方程。为实现高效、稳定的数值模拟，使用有限差分-伪谱混合法求解波动方程。最终，建立了高精度、高效率、稳定的TTI介质正演模拟与逆时偏移成像算法。数值模拟实验表明，所提方法可实现复杂TTI介质偏移成像。

	服务

	把本文推荐给朋友
	加入我的书架
	加入引用管理器
	E-mail Alert
	RSS
	作者相关文章
	慕鑫茹
	黄建平
	李振春
	郭旭
	雍鹏
	袁双齐

关键词 ：逆时偏移, TTI介质, qP波, qSV波, 最佳平方逼近, 有限差分—伪谱混合法

Abstract：In the seismic data processing,if the influence of anisotropy and the symmetry inclination of TTI medium changes are ignored,the imaging position of targets will be deviated,and even the amplitude imbalance may be caused by the caustic problem.The pseudo-shear-wave interference exists in the forward modeling and RTM when using the conventional TI coupled pseudo-acoustic wave equations,and the numerical simulation is unstable when the anisotropy parameter ε<δ.To solve the problem,a new approach is proposed.By referring to the experience of predecessors,the idea of the optimal quadratic approximation (OQA) is used to approximate the accurate phase velocity formula in VTI medium,and approximate phase velocity formulas of qP-wave and qSV-wave with relatively high accuracy are obtained respectively.Then,based on the approximate phase velocity formula,the decoupled qP-wave and qSV-wave equations in TTI medium are derived.In order to realize efficient and stable numerical simulation,a hybrid finite-difference/pseudo-spectral method is used to solve the wave equations.Finally,a high-precision,high-efficiency,and stable TTI media forward modeling and RTM is established.Numerical simulation results show that the proposed method can achieve migration imaging in complex TTI media.

Key words： reverse time migration (RTM) TTI medium pure qP-wave pure qSV-wave optimal quadratic approximation (OQA) hybrid finite-difference/pseudo-spectral method

收稿日期: 2019-03-20

PACS:

P631

基金资助:本项研究受泰山学者青年人才工程专项基金项目"泰山学者青年专家人才计划"（Tsqn20161016））、国家油气重大专项"基于多次散射理论的散射波地震成像技术"（2016ZX05014-008-005）和"薄互层全波形反演和最小二乘偏移联合成像"（2016ZX05002-005-002））联合资助。

通讯作者: 黄建平,山东省青岛市经济技术开发区长江西路66号中国石油大学(华东)地球科学与技术学院,266580。Email:jphuang@upc.edu.cn E-mail: jphuang@upc.edu.cn

作者简介: 慕鑫茹,博士研究生,1995年生;2017年获中国石油大学(华东)勘查技术与工程专业学士学位;现在该校攻读地质资源与地质工程专业(物探方向)博士学位,主要研究方向为复杂各向异性介质地震波传播与成像。

引用本文:

慕鑫茹, 黄建平, 李振春, 郭旭, 雍鹏, 袁双齐. 基于最佳平方逼近的TTI介质解耦qP波与qSV波逆时偏移[J]. 石油地球物理勘探, 2019, 54(6): 1280-1292. MU Xinru, HUANG Jianping, LI Zhenchun, GUO Xu, YONG Peng, YUAN Shuangqi. Reverse time migration of decoupled qP- and qSV-waves in TTI media with the optimal quadratic approximation. Oil Geophysical Prospecting, 2019, 54(6): 1280-1292.

链接本文:

http://www.ogp-cn.com/CN/10.13810/j.cnki.issn.1000-7210.2019.06.011 或 http://www.ogp-cn.com/CN/Y2019/V54/I6/1280

[1]	Thomsen L.Weak elastic anisotropy[J].Geophysics,1986,51(10):1954-1966.
[2]	吴国忱.各向异性介质地震波传播与成像[M].山东东营:石油大学出版社,2006.
[3]	Du X,Fletcher R P,and Fowler P J.A new pseudo-acoustic wave equation for VTI media[C].Extended Abstracts of 70th EAGE Conference & Exhibition,2008,H033.
[4]	段新意,李振春,黄建平,等.各向异性介质共炮域高斯束叠前深度偏移[J].石油物探,2014,53(5):579-586.DUAN Xinyi,LI Zhenchun,HUANG Jianping,et al.A prestack Gaussian beam depth migration in common-shot domain for anisotropic media[J].Geophysical Prospecting for Petroleum,2014,53(5):579-586.
[5]	黄金强,李振春.基于Low-rank分解的复杂TI介质纯qP波正演模拟与逆时偏移[J].地球物理学报,2017,60(2):704-721.HUANG Jinqiang,LI Zhenchun.Modeling and reverse time migration of pure quasi-P-waves in complex TI media with a low-rank decomposition[J].Chinese Journal of Geophysics,2017,60(2):704-721.
[6]	Crampin S.Anisotropy in exploration seismics[J].First Break,1984,doi:10.3997/1365-2397.1984006.
[7]	Winterstein D F.Velocity anisotropy terminology for geophysicists[J].Geophysics,1990,55(8):1070-1088.
[8]	杜启振,秦童.横向各向同性介质弹性波多分量叠前逆时偏移[J].地球物理学报,2009,52(3):801-807.DU Qizhen,QIN Tong.Multicomponent prestack reverse-time migration of elastic waves in transverse isotropic medium[J].Chinese Journal of Geophysics,2009,52(3):801-807.
[9]	Yong P,Huang J P,Li Z C,et al.Elastic-wave reverse-time migration based on decoupled elastic-wave equations and inner-product imaging condition[J].Journal of Geophysics & Engineering,2016,13(6):953-963.
[10]	Zhou H,Zhang G,and Bloor R.An anisotropic acoustic wave equation for VTI media[C].Extended Abstracts of 68th EAGE Conference & Exhibition,2006,H033.
[11]	Du X,Bancroft J C,and Lines L R.Anisotropic reverse-time migration for tilted TI media[J].Geophysical Prospecting,2007,55(6):853-869.
[12]	Duveneck E,Bakker P M.Stable P-wave modeling for reverse-time migration in tilted TI media[J].Geophysics,2011,76(2):S65-S75.
[13]	王伟国,熊水金,徐华宁,等.TTI介质各向异性伪谱法逆时偏移[J].石油地球物理勘探,2012,47(4):566-572.WANG Weiguo,XIONG Shuijin,XU Huaning,et al.Reverse-time migration using pseudo-spectral method for tilted TI media[J].Oil Geophysical Prospecting,2012,47(4):566-572.
[14]	程玖兵,康玮,王腾飞.各向异性介质qP波传播描述Ⅰ:伪纯模式波动方程[J].地球物理学报,2013,56(10):3474-3486.CHENG Jiubing,KANG Wei,WANG Tengfei.Description of qP-wave propagation in anisotropic media,Part Ⅰ:Pseudo-pure-mode wave equations[J].Chinese Journal of Geophysics,2013,56(10):3474-3486.
[15]	刘文卿,王西文,王宇超,等.带正则化校正的TTI介质qP波方程及其逆时偏移方法和应用[J].地球物理学报,2016,59(3):1059-1069.LIU Wenqing,WANG Xiwen,WANG Yuchao,et al.A regularized qP-wave equation for TTI media and its application to reverse time migration[J].Chinese Journal of Geophysics,2016,59(3):1059-1069.
[16]	李振春,杨富森,王小丹.基于LS-RSGFD方法优化的横向各向同性(TI)介质一阶qP波高精度数值模拟[J].地球物理学报,2016,59(4):1477-1490.LI Zhenchun,YANG Fusen,WANG Xiaodan.High-precision numerical simulation of first-order qP-waves in the transversely isotropic (TI) medium optimized by the LS-RSGFD method[J].Chinese Journal of Geo-physics,2016,59(4):1477-1490.
[17]	胡书华,王宇超,刘文卿,等.复杂TTI介质稳定的纯qP波波场模拟方法[J].石油地球物理勘探,2018,53(2):280-287.HU Shuhua,WANG Yuchao,LIU Wenqing,et al.Pure quasi-P wave stable simulation in complex TTI media[J].Oil Geophysical Prospecting,2018,53(2):280-287.
[18]	Tsvankin I.P-wave signatures and notation for transversely isotropic media:An overview[J].Geophysics,1996,61(2):467-483.
[19]	Alkhalifah T.Acoustic approximations for processing in transversely isotropic media[J].Geophysics,1998,63(2):623-631.
[20]	Alkhalifah T.An acoustic wave equation for aniso-tropic media[J].Geophysics,2000,65(4):1239-1250.
[21]	Klíe H and W Toro.A new acoustic wave equation for modeling in anisotropic media[C].SEG Technical Program Expanded Abstracts,2001,20:1171-1174.
[22]	Hestholm S.Acoustic VTI modeling using high-order finite differences[C].SEG Technical Program Expanded Abstracts,2007,26:139-143.
[23]	Duveneck E,Milcik P,Bakker P M,et al.Acoustic VTI wave equations and their application for anisotropic reverse-time migration[C].SEG Technical Program Expanded Abstracts,2008,27:2186-2190.
[24]	Zhang Y and Zhang H.A stable TTI reverse time migration and its implementation[C].SEG Technical Program Expanded Abstracts,2009,28:2794-2798.
[25]	杨富森,李振春,王小丹.TTI介质一阶qP波稳定方程波场数值模拟及逆时偏移[J].石油地球物理勘探,2016,51(3):497-505.YANG Fusen,LI Zhenchun,WANG Xiaodan.Wave-field numerical simulation and reverse-time migration based on the stable TTI first-order qP-wave equations[J].Oil Geophysical Prospecting,2016,51(3):497-505.
[26]	张岩,吴国忱.TTI介质qP波逆时偏移中伪横波噪声压制方法[J].地球物理学报,2013,56(6):2065-2076.ZHANG Yan,WU Guochen.Methods of removing pseudo SV-wave artifacts in TTI media qP-wave reverse-time migration[J].Chinese Journal of Geophy-sics,2013,56(6):2065-2076.
[27]	Zhang H,Zhang G,and Zhang Y.Removing S-wave noise in TTI reverse time migration[C].SEG Technical Program Expanded Abstracts,2009,28:2849-2853.
[28]	李博,李敏,刘红伟,等.TTI介质有限差分逆时偏移的稳定性探讨[J].地球物理学报,2012,55(4):1366-1375.LI Bo,LI Min,LIU Hongwei,et al.Stability of reverse time migration in TTI media[J].Chinese Journal of Geophysics,2012,55(4):1366-1375.
[29]	Chu C,Macy B K,and Anno P D.Approximation of pure acoustic seismic wave propagation in TTI media[J].Geophysics,2011,76(5):WB97-WB107.
[30]	Zhan G,Pestana R C,and Stoffa P L.Decoupled equations for reverse time migration in tilted transversely isotropic media[J].Geophysics,2012,77(2):37-45.
[31]	Song G,Huang R,Tian J,et al.A new qP-wave equation for 2D VTI media[C].SEG Technical Program Expanded Abstracts,2016,35:3977-3981.
[32]	Mu X R,Huang J P,Huang J Q,et al.A precise si-mulation of pure qP-wave in VTI media[C].Extended Abstracts of 80th EAGE Conference & Exhibition,2018,doi:10.3997/2214-4609.201800970.
[33]	Li X and Zhu H.A finite-difference approach for solving pure quasi-P wave equations in transversely isotropic and orthorhombic media[J].Geophysics,2018,83(4):C161-C172.
[34]	Zhan G,Pestana R C,and Stoffa P L.An efficient hybrid pseudo spectral/finite difference scheme for solving TTI pure P wave equation[C].13rd International Congress of the Brazilian Geophysics Society & EXPOGEF,Rio de Janeior,Brazil,2013,1322-1327.
[35]	Kosloff D D,Baysal E.Forward modeling by a Fourier method[J].Geophysics,1982,47(10):1402-1412.
[36]	Yang J,Zhang S,Zhu H.Isotropic elastic wavefields decomposition using fast Poisson solvers[C].SEG Technical Program Expanded Abstracts,2017,36:4716-4720.
[37]	Baysal E,Dan D K,Sherwood J W C.Reverse time migration[J].Geophysics,1983,48(11):1514-1524.
[38]	刘红伟,李博,刘洪,等.地震叠前逆时偏移高阶有限差分算法及GPU实现[J].地球物理学报,2010,53(7):1725-1733.LIU Hongwei,LI Bo,LIU Hong,et al.The algorithmof high order finite difference pre-stack reverse time migration and GPU implementation[J].Chinese Journal of Geophysics,2010,53(7):1725-1733.
[39]	李振春,李庆洋,黄建平,等.一种稳定的高精度双变网格正演模拟与逆时偏移方法[J].石油物探,2014,53(2):127-136.LI Zhenchun,LI Qingyang,HUANG Jianping,et al.A stable and high-precision dual-variable grid forward modeling and reverse time migration method[J].Geophysical Prospecting for Petroleum,2014,53(2):127-136.
[40]	Feng B,Wang H Z.Reverse time migration with source wavefield reconstruction strategy[J].Journal of Geophysics & Engineering,2011,9(1):69-74.
[41]	王保利,高静怀,陈文超,等.地震叠前逆时偏移的有效边界存储策略[J].地球物理学报,2012,55(7):2412-2421.WANG Baoli,GAO Jinghuai,CHEN Wenchao,et al.Efficient boundary storage strategies for seismic reverse time migration[J].Chinese Journal of Geophy-sics,2012,55(7):2412-2421.