Please wait a minute...
E-mail Alert Rss
 
物探与化探  2018, Vol. 42 Issue (1): 134-143    DOI: 10.11720/wtyht.2018.1.16
  本期目录 | 过刊浏览 | 高级检索 |
TTI介质拟声波方程数值模拟
黄杰(), 杨国权, 李振春, 谷丙洛
中国石油大学(华东) 地球科学与技术学院,山东 青岛 266580
Numerical simulation of pseudo acoustic wave equation in TTI medium
Jie HUANG(), Guo-Quan YANG, Zhen-Chun LI, Bing-Luo GU
School of Geosciences,China University of Petroleum (East China),Qingdao 266580,China
全文: PDF(2200 KB)   HTML
输出: BibTeX | EndNote (RIS)      
摘要 

TTI介质qP波数值模拟方法因为考虑了倾角因素,可以比VTI介质qP波数值模拟方法更加准确地描述各向异性介质中地震波场的传播规律。文中用拟声波方程对TTI介质中的地震波场进行了高阶有限差分数值模拟,在改进衰减函数分布方式后,通过坐标变换,利用改进的完全匹配层(perfectly matched layer,PML)边界控制方程对波场边界进行吸收处理,取得了良好的效果;然后分析了拟声波方程数值模拟中的稳定性问题,并对波场中的伪横波进行压制。通过对不同模型的数值模拟,验证了文中使用的TTI介质拟声波波动方程的稳定性以及所采用的PML边界控制方程的可靠性和适用性。

服务
把本文推荐给朋友
加入引用管理器
E-mail Alert
RSS
作者相关文章
黄杰
杨国权
李振春
谷丙洛
关键词 TTI介质拟声波方程PML边界控制方程稳定性伪横波    
Abstract

TTI qP wave numerical simulation method is more accurate than VTI qP wave numerical simulation method in describing the propagation law of seismic wavefield in anisotropic medium because of considering the dip angle factor.Using the pseudo acoustic wave equation,the authors carried out high order finite difference numerical simulation of the seismic wavefield in the TTI medium in this paper.After improving the distribution of the attenuation function,the authors used the improved perfectly matched layer (PML) boundary control equation to deal with the wavefield boundary and achieved good results.Then the pseudo shear wave in the numerical simulation of pseudo acoustic wave equation was suppressed and the stability problem was analyzed.The numerical simulation of different models prove that the TTI medium pseudo acoustic wave equation is stable and the PML boundary control equations have high reliability and applicability.

Key wordsTTI medium    pseudo acoustic wave equation    PML boundary control equation    stability    pseudo shear wave
收稿日期: 2016-11-18      出版日期: 2018-02-20
:  P631.4  
基金资助:国家重点基础研究发展计划(“973”计划)课题(2014CB239006),国家自然科学基金项目(41504100)
作者简介:

作者简介: 黄杰(1991-),男,硕士研究生,主要从事地震波正演模拟与逆时偏移工作。Email:18205420120@163.com

引用本文:   
黄杰, 杨国权, 李振春, 谷丙洛. TTI介质拟声波方程数值模拟[J]. 物探与化探, 2018, 42(1): 134-143.
Jie HUANG, Guo-Quan YANG, Zhen-Chun LI, Bing-Luo GU. Numerical simulation of pseudo acoustic wave equation in TTI medium. Geophysical and Geochemical Exploration, 2018, 42(1): 134-143.
链接本文:  
https://www.wutanyuhuatan.com/CN/10.11720/wtyht.2018.1.16      或      https://www.wutanyuhuatan.com/CN/Y2018/V42/I1/134
  PML边界处理示意
  本文采用的PML边界衰减函数分布方式示意

a—x方向衰减函数dx;b—z方向衰减函数dz

  稳定性对比

a—ε<δ;b—εδ

  二维TTI介质中传播时间为400 ms时p的波场快照

a—θ=0°;b—θ=30°;c—θ=45°;d—θ=90°

  二维TTI介质中传播时间为400 ms时q的波场快照

a—θ=0°;b—θ=30°;c—θ=45°;d—θ=90°

  θ=30°时800 ms的波场快照对比

a—未加边界时;b—应用改进的PML边界

  θ=30°时的单炮地震记录对比

a—未加边界时;b—应用改进的PML边界

  θ=0°时应用本文的PML边界条件的波场快照

a—800 ms;b—2000 ms

  压制伪横波后400 ms的波场快照

a—θ=0°;b—θ=30°

  BP_partI模型参数

a—δ;b—ε;c—θ;d—vp

  BP_partI模型波场快照

a—1200 ms;b—1400 ms;c—1600 ms;d—1800 ms

  BP_partI模型单炮地震记录
  BP_partII模型参数

a—δ;b—ε;c—θ;d—vp

  BP_partII模型波场快照

a—800 ms;b—1000 ms;c—1400 ms;d—1600 ms

  BP_partII模型单炮地震记录
  BP模型不同介质情况下1 200 ms的波场快照

a—各向同性介质;b—VTI介质;c—TTI介质

[1] 牟永光,裴正林.三维复杂介质地震数值模拟[M].北京:石油工业出版社,2005:1-13.
[2] Alkhalifah T.Acoustic approximations for processing in transversely isotropic media[J].Geophysics,1998,63(2):623-631.
[3] Alkhalifah T.An acoustic wave equation for anisotropic media[J].Geophysics,2000,65(4):1239-1250.
[4] Zhou H,Zhang G,Bloor R.An anisotropic acoustic wave equation for VTI media[C]//68th EAGE Conference and Exhibition Incorporating SPE EUROPEC 2006.Vienna,Austria:EAGE,2006:H033.
[5] Du X,Fletcher R,Fowler P J.A new pseudo-acoustic wave equation for VTI media[C]//70th EAGE Conference and Exhibition Incorporating SPE EUROPEC 2008.Rome:European Association of Geoscientists and Engineers,2008:1-5.
[6] Hestholm S.Acoustic VTI modeling using high-order finite differences[J].Geophysics,2009,74(5):T67-T73.
[7] 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 International Meeting.Tulsa,OK:SEG,2006:194-198.
[8] 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.Tulsa,OK:SEG,2008:2082-3086.
[9] Duveneck E,Milcik P,Bakker P M,et al.Acoustic VTI wave equations and their application for anisotropic reverse-time migration[C]//Expanded Abstracts of 78th SEG Annual International Meeting.Tulsa,OK:SEG,2008:2186-2190.
[10] Duveneck E,Bakker P M.Stable P-wave modeling for reverse-time migration in tilted TI media[J].Geophysics,2011,76(2):S65-S75.
[11] 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.
[12] Liu F Q,Morton S A,Jiang S S,et al.Decoupled wave equations for P and SV waves in an acoustic VTI media[C]//SEG Technical Program Expanded Abstracts 2009.Tulsa,OK:SEG,2009:2844-2848.
[13] Pestana R C,Ursin B,Stoffa P L.Separate P- and SV-wave equations for VTI media[C]//SEG Technical Program Expanded Abstracts 2011.Tulsa,OK:SEG,2011:163-167.
[14] Chu C L,Macy B K,Anno P D.An accurate and stable wave equation for pure acoustic TTI modeling[C]//SEG Technical Program Expanded Abstracts 2011.Tulsa,OK:SEG,2011:179-184.
[15] Zhan G,Pestana R C,Stoffa P L.An acoustic wave equation for pure P wave in 2D TTI media[C]//SEG Technical Program Expanded Abstracts 2011.Tulsa,OK:SEG,2011:168-173.
[16] Barrera Pacheco D F,Pestana R,Vivas F.New pseudo-acoustic wave equations for modeling and reverse time migration in TTI media[C]//75th EAGE Conference & Exhibition Incorporating SPE EUROPEC 2013.London:SPE,2013.
[17] Thomsen L.Weak elastic anisotropy[J].Geophysics,1986,51(10):1954-1966.
[18] 张岩,吴国忱.TTI介质qP波逆时偏移中伪横波噪声压制方法[J].地球物理学报,2013,56(6):2065-2076.
[19] 刘洋,李承楚,牟永光.任意偶数阶精度有限差分法数值模拟[J].石油地球物理勘探,1998,33(1):1-10.
[20] Berenger J P.A perfectly matched layer for the absorption of electromagnetic waves[J].Journal of Computational Physics,1994,114(2):185-200.
[21] Collino F,Tsogka C.Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media[J].Geophysics,2001,66(1):294-307.
[22] Komatitsch D,Tromp J.A Perfectly Matched Layer absorbing boundary condition for the second-order seismic wave equation[J].Geophysical Journal International,2003,154(1):146-153.
[23] Hastings F D,Schneider J B,Broschat S L.Application of the perfectly matched layer (PML) absorbing boundary condition to elastic wave propagation[J].The Journal of the Acoustical Society of America,1996,100(5):3061-3069.
[24] Groby J P,Tsogka C.A time domain method for modeling viscoacoustic wave propagation[J].Journal of Computational Acoustics,2004,14(2):201-236.
[25] Tsvankin I.Seismic signatures and analysis of reflection data in anisotropic media[M].London:Pergamon,2001.
[26] Grechka V,Zhang L B,Rector J W.Shear waves in acoustic anisotropic media[J].Geophysics,2004,69(2):576-582.
[27] Xu S,Zhou H B.Accurate simulations of pure quasi-P-waves in complex anisotropic media[J].Geophysics,2014,79(6):T341-T348.
[1] 单希鹏, 谢汝宽, 梁盛军, 余学中. 直升机TEM测量影响因素分析[J]. 物探与化探, 2021, 45(1): 178-185.
[2] 陈晓晶, 虎新军, 李宁生, 安百州, 白亚东. 银川平原基于地球物理资料三维建模的深部地质构造研究[J]. 物探与化探, 2020, 44(2): 245-253.
[3] 刘芬, 王万银, 纪晓琳. 空间域和频率域平面位场延拓影响因素和稳定性分析[J]. 物探与化探, 2019, 43(2): 320-328.
[4] 何涛, 王万银, 黄金明, 张明华, 杨敏. 正则化方法在比值类位场边缘识别方法中的研究[J]. 物探与化探, 2019, 43(2): 308-319.
[5] 王彦国, 张瑾. 位场高阶导数的波数域迭代法[J]. 物探与化探, 2016, 40(1): 143-147.
[6] 王小丹, 印兴耀, 杨富森. 稳定的二维TTI介质一阶qP波方程RSGFD数值模拟[J]. 物探与化探, 2016, 40(1): 135-142.
[7] 李长江, 李庆春, 王亚妮. 旋转交错网格TTI介质波场模拟与波场分解[J]. 物探与化探, 2015, 39(3): 553-557.
[8] 张伟, 甘伏平, 刘伟, 郑智杰. 双相介质瑞雷面波有限差分正演模拟[J]. 物探与化探, 2014, 38(6): 1275-1283.
[9] 郭伟立, 薛国强, 周楠楠, 肖宏跃. 利用瞬变电磁法监测煤矿含水采空区[J]. 物探与化探, 2012, 36(S1): 114-118.
[10] 刘百红, 周巍, 马中高. 含倾斜定向排列裂隙介质纵波特征[J]. 物探与化探, 2012, 36(3): 467-472.
[11] 刘妹, 顾铁新, 史长义, 鄢卫东, 鄢明才. 我国主要土壤类型元素地球化学形态成分标准物质研制[J]. 物探与化探, 2008, 32(5): 492-496,508.
[12] 解华明, 丁华. PLS法与隧道围岩稳定性分类[J]. 物探与化探, 2003, 27(4): 320-322.
[13] 孙秀容, 刘志高, 许东儿, 戴纯俊. 上海市区区域工程场地剪切波地震稳定性评价[J]. 物探与化探, 2002, 26(5): 398-402.
[14] 王卫平. 辽宁西部地区地壳结构及其稳定性的初步研究[J]. 物探与化探, 1997, 21(5): 354-362.
[15] 熊章强. 218Po法测氡对甘肃疏勒河断裂带中段的稳定性研究[J]. 物探与化探, 1997, 21(3): 219-224.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
京ICP备05055290号-3
版权所有 © 2021《物探与化探》编辑部
通讯地址:北京市学院路29号航遥中心 邮编:100083
电话:010-62060192;62060193 E-mail:whtbjb@sina.com