双相各向异性随机介质伪谱法地震波场特征分析
- 格式:pdf
- 大小:412.68 KB
- 文档页数:7
模拟地震波场的伪谱和高阶有限差分混合方法
魏星;王彦宾;陈晓非
【期刊名称】《地震学报》
【年(卷),期】2010(032)004
【摘要】伪谱法是一种高效、高精度计算非均匀介质地震波传播的数值方法, 但是由于它的微分算子的全局性, 使得该方法不适用于分散内存的并行计算. 本文将有限差分算子的局部性和伪谱法算子的高效、高精度相结合, 发展基于两种方法的伪谱/有限差分混合方法. 该方法在一个空间坐标方向上利用交错网格高阶有限差分算子, 在另外的空间坐标方向上利用交错网格伪谱法算子, 既保留了后者的高效、高精度优势, 又便于在PC集群上实现并行计算. 对二维模型的计算显示, 混合方法能有效处理介质不连续面, 在保证伪谱法计算精度的情况下, 提供了一种并行计算的可能途径.
【总页数】9页(P392-400)
【作者】魏星;王彦宾;陈晓非
【作者单位】中国北京,100045,中国地震台网中心;中国北京,100871,北京大学地球与空间科学学院地球物理学系;中国合肥,230026,中国科学技术大学地球和空间科学学院
【正文语种】中文
【中图分类】P315.3+1
【相关文献】
1.基于图形处理器的伪谱和高阶有限差分混合方法地震波数值模拟 [J], 崔丛越;张献兵;王彦宾
2.地震波传播的三维伪谱和高阶有限差分混合方法并行模拟 [J], 秦艳芳;王彦宾
3.二维非均匀介质地震波传播的伪谱和有限差分混合方法的应用研究 [J], 李少华;王彦宾;吴志坚
4.基于伪谱和有限差分混合方法的兰州盆地强地面运动二维数值模拟 [J], 严武建;王彦宾;石玉成
5.基于伪谱和有限差分混合方法的兰州盆地强地面运动二维数值模拟 [J], 严武建;王彦宾;石玉成;
因版权原因,仅展示原文概要,查看原文内容请购买。
第57卷第2期 2018年3月石油物探G E O P H Y S I C A L P R O S P E C T I N G F O R P E T R O L E U MVol.57,No.2Mar. ,2018宋利伟,石颖,柯璇,等.基于检查点方法的各向异性介质逆时偏移[J].石油物探,2018,57(2): 274-282S O N G L i w e l S H I Y i n g^K E X u a n=t al.^Reverse time migration of anisotropic media using the checkpoint method[J#(jeophysical Prospecting for Petroleum,2018,57(2) :274-282基于检查点方法的各向异性介质逆时偏移宋利伟,石颖,柯璇,王维红,李松龄(东北石油大学,黑龙江大庆163318)摘要:基于准P(quas-P,qP)波方程研究了 TTI介质逆时偏移方法,通过引人横波速度,改善数值模拟的稳定性。
提出了联合有 效内边界存储和检查点(heckpont)方法降低存储需求的TTI介质逆时偏移成像方法,实现了震源波场的精确重建,为高精度 逆时偏移成像奠定了基础,也使方法更好地适应大数据计算。
借助图形处理器(Graphc P r〇C e S S mgUmt,GPU)的高效并行计算 能力,大幅度提升震源波场重建以及TTI介质逆时偏移成像的计算效率。
应用BP2007各向异性理论模型数据进行测试的结果 表明,TTI介质逆时偏移成像方法成像精度高、适用性强、稳定性好。
关键词:各向异性介质;逆时偏移;检查点;存储策略;图形处理器;波场重建中图分类号:P631文献标识码:A文章编号:1000-1441 (2018)02-0274-09D01:10.3969/j.iss m1000-1441.2018.02.013Reverse time migration of anisotropic media using the checkpoint method S O N G L i w e s S H I Y i n g,K E X u a n,W A N G W e ih o n g,L I S o n g lin g(.Northeast Petroleum University ^Daqing 163318,—$.)Abstract: A reverse time migration method for tilted transversely isotropic (TTI) media that combines the efficient internal boundary storage and checkpoint method to reduce storage i s proposed. In this method, shear-wave velocity i s introduced into t^h e qP-wave equation to improve the stability of the numerical smulatiomTo reduce the required storage,a combined efficient mternal boundary storage and checkpoint method i s proposed,thus realizing accurate reconstruction of source wave fields and improving i t s applicability to large data computation Parallel computing of graphic processing units s used to enhance the efficiency of the source wave field recalculation and reverse time migration for T T I media. A test using B P 2007 Anisotropic Velocity Benchmark data revealed the feasibility and effectiveness of the proposed method.Keywords: anisotropic media, reverse time migration, checkpoint ^storage strategy, graphic processing unit, wave field reconstruction基于双程波动方程的逆时偏移方法,无地层倾角限制,适应于各种波现象,是复杂地下构造成像的有效技术[13]。
基于双相介质的地震波衰减特性油气检测方法
汪佳蓓;黄捍东;郭飞
【期刊名称】《东北石油大学学报》
【年(卷),期】2016(040)002
【摘要】油气储层是典型的多相介质,目前油气检测中依据的单相介质理论容易造成预测结果的多解性.为了更准确地检测油气,以Biot双相介质理论为基础,推导吸收系数表达式,建立能量密度衰减方程,利用广义S变换对地震资料进行频谱分解,分析介质含不同流体后地震波各频率成分的衰减特征.结果表明:当双相介质中所含流体性质不同时,地震波的频率衰减也不尽相同,当储集层中含油气时,地震波频率衰减量显著增大,低频能量衰减小,高频能量衰减大,在频谱上表现为低频能量相对增强,高频能量相对减弱,整体能量向低频移动.对比实际工区的低、高频能量剖面和衰减剖面,将地震波的衰减规律运用到油气检测中取得较好的效果,理论模型和实际资料的运用表明该方法可行.
【总页数】8页(P19-26)
【作者】汪佳蓓;黄捍东;郭飞
【作者单位】中国地质大学(北京)能源学院,北京100083;中国石油大学(北京)油气资源与探测国家重点实验室,北京102249;中国石油大学(北京)提高采收率研究院,北京 102249;中海石油深圳分公司研究院,广东广州510240
【正文语种】中文
【中图分类】P631.4;TE122.2
【相关文献】
1.基于储层预测的双相介质理论油气检测方法应用研究
2.基于双相介质的地震油气检测方法——以Y925地区为例
3.利用地震波在双相介质中的衰减特性检测油气
4.基于双相介质的地震波衰减特性油气检测方法
5.基于双相介质油气检测方法在塔东GC地区的应用
因版权原因,仅展示原文概要,查看原文内容请购买。
三维各向异性介质中的地震波模拟
谢桂生
【期刊名称】《石油物探》
【年(卷),期】1994(033)003
【摘要】本文给出了一种在三维各向异性介质中模拟地震波传播的方案,该方案是利和Fourier方法计算空间导数和时间导数,用蛙跳技术计算来完成模拟的。
允许弹性常数和密度在横向上和纵向上任意变化,由于采用了Fourier方法计算空间导数,可获得较高的模拟精度,并可有效地减少工作空间,这在三维模拟中是十分重要的。
模拟方案是在Micro-VAX3500机上实现的,尽管本方案可用于计算任意弹性各向异性介质,文中给出了
【总页数】5页(P54-58)
【作者】谢桂生
【作者单位】无
【正文语种】中文
【中图分类】P631.41
【相关文献】
1.复杂各向异性介质中的地震波场有限差分模拟 [J], 张会星;何兵寿;张晶;苑春方
2.二维各向异性介质中地震波场的高阶同位网格有限差分模拟 [J], 祝贺君;张伟;陈晓非
3.改进的FCT方法在粘滞各向异性介质地震波模拟中的研究 [J], 李勤武;周怀来
4.基于GPU和CPU协同并行的三维各向异性介质地震波场正演模拟 [J], 刘春成;
顾汉明;陈宝书;焦振华;马凯;蔡志成;张立
5.粘弹各向异性介质中地震波场模拟与特征 [J], 郭智奇;刘财;杨宝俊;刘洋;王典因版权原因,仅展示原文概要,查看原文内容请购买。
第45卷 第1期2023年1月物探化探计算技术COMPUTINGTECHNIQUESFORGEOPHYSICALANDGEOCHEMICALEXPLORATIONVol.45 No.1Jan.2023收稿日期:2021 10 23基金项目:国家自然科学基金创业创新发展联合基金(U19B603)第一作者:王鸣谦(1995-),男,硕士,主要研究方向为地震波正演及偏移成像,E mail:285639863@qq.com。
文章编号:1001 1749(2023)01 0001 07基于双相介质弹性波波场分离的逆时偏移方法王鸣谦1,张胜寒1,张济东1,付继有2(1.中国石油化工股份有限公司 石油物探技术研究院,南京 2111032.中国石油 长庆油田页油岩开发分公司,甘肃 745000)摘 要:随着地震勘探的发展,以单相介质模型为基础的声波逆时偏移和弹性波逆时偏移很难比较精确地描述地下介质和满足复杂的工程要求,尤其是在处理地下含流体多孔介质时,所以需要引入双相介质模型。
首先基于Biot理论双相介质一阶速度-应力弹性波方程,利用散度场和旋度场的波场分离理论将固相和流相中的纵横波场解耦分离,并由波场快照图得出固相与流相中的纵波波场与横波波场均完成分离,而快纵波与慢纵波相互耦合均存在于纵波波场中。
采用高阶交错网格有限差分方法,结合PML吸收边界条件,推导了具有波场分离的双相介质弹性波正传与反传的偏移延拓算子,结合互相关成像条件,可以得到双相介质弹性波固相与流相中的各分量偏移成像剖面。
最后对模型使用该双相介质弹性波波场分离的逆时偏移方法进行偏移成像试算并取得成功,效果较常规双相介质弹性波逆时偏移成像效果有了明显提升。
关键词:Biot理论;双相介质;波场分离;弹性波逆时偏移中图分类号:P631.4 文献标志码:A 犇犗犐:10.3969/j.issn.1001 1749.2023.01.010 前言由于实际情况中含有油气的储层是具有固体与流体两种状态的双相介质,因此建立在均匀各向同性的理论基础上研究波动传播规律并不是十分准确的,只是将储层简化为单相弹性介质,从而忽略了固相和流相之间的耦合作用,因此对地震波的传播描述有所偏差,所以建立更加符合实际地层情况的双相介质模型具有重大意义。
非均匀弹性介质地震波场伪谱法正演模拟
逄硕;刘财
【期刊名称】《世界地质》
【年(卷),期】2014(33)4
【摘要】基于一阶速度—应力弹性波波动方程,利用傅里叶变换将空间导数变换到频率域求解,对介质波场进行正演模拟.本文结合实际勘探中可能存在的地质问题,设计了多种地质模型,包括背斜模型、向斜模型、透镜体模型和断层-侵入体模型等,进行了伪谱法地震波场正演模拟、合成地震记录的实现,并分析了其波场动力学变化特征.结果表明,应用伪谱法对复杂模型的模拟有较高的精度和可靠性,可直接用于辅助地震野外工作设计、反演解释和资料解释等工作.
【总页数】6页(P934-939)
【作者】逄硕;刘财
【作者单位】吉林大学地球探测科学与技术学院,长春130026;吉林大学地球探测科学与技术学院,长春130026
【正文语种】中文
【中图分类】P631.443
【相关文献】
1.各向异性粘弹性孔隙介质伪谱法波场正演模拟 [J], 杜启振;刘莲莲;孙晶波
2.一种基于Kirchhoff积分的快速三维非均匀弹性介质地震波场模拟方法 [J], 张钋;Mathias Alerini;Gilles Lambaré
3.各向异性弹性介质方向行波波场分离正演数值模拟 [J], 陈可洋
4.含直立裂缝粘弹性介质地震波场正演模拟 [J], 韩令贺;何兵寿;张会星
5.伪谱法运算效率的提高及井间波场的正演模拟 [J], 田仁飞;曹俊兴
因版权原因,仅展示原文概要,查看原文内容请购买。
二维随机介质模型正演模拟及其波场分析殷学鑫;刘洋【期刊名称】《石油地球物理勘探》【年(卷),期】2011(046)006【摘要】本文通过弹性波方程交错网格有限差分正演,模拟了地震波在二维弹性随机介质中的传播,得到了球面波传播的单炮记录和波场快照以及平面波自激自收记录.为了考察二维弹性随机介质的波场特征,将平面波自激自收记录分成直达波区、散射波区、反射波区三个区域,并在三个区域中分别计算了能量相对值、主频、有效频带宽度三个特征量.通过研究这三个特征量与随机介质模型的自相关长度、粗糙度因子、扰动标准差等因素的关系,得到以下结论:直达波、散射波、反射波的能量相对值随随机介质模型自相关长度、粗糙度因子、扰动标准差以及震源主频的变化有明显的变化.直达波的主频和频带宽度受这些因素的影响很小;散射波的主频主要受自相关长度的影响,散射波的有效频带宽度很窄,受各因素影响也很小;反射波的主频随自相关长度、粗糙度因子、扰动标准差的变化有明显的变化,而反射波的频带宽度主要受扰动标准差的影响.【总页数】11页(P862-872)【作者】殷学鑫;刘洋【作者单位】中国石油大学(北京)油气资源与探测国家重点实验室,北京昌平,102249;中国石油大学(北京)CNPC物探重点实验室,北京昌平,102249;中国石油大学(北京)油气资源与探测国家重点实验室,北京昌平,102249;中国石油大学(北京)CNPC物探重点实验室,北京昌平,102249【正文语种】中文【中图分类】P631【相关文献】1.基于随机介质模型的近地表波场模拟及分析 [J], 肖云飞;殷厚成;孙成禹;陈林;常鉴;徐国庆;彭代平2.随机介质模型正演模拟及其地震波场分析 [J], 姚姚;奚先3.随机介质模型中瑞雷波正演模拟及波场特性分析 [J], 赵宏欣;熊章强;张大洲4.随机介质背散射二维Mueller矩阵的数值模拟与分析 [J], 王清华;李振华;来建成;贺安之5.二维粘弹性随机介质中的波场特征分析 [J], 奚先;姚姚因版权原因,仅展示原文概要,查看原文内容请购买。
VTI介质多波各向异性参数分析李勤*① 李庆春② 张林③【摘要】摘要在地震数据处理中,通过各向同性参数分析无法校平大炮检距地震记录,会引起动校拉伸、叠加结果失真,而通过各向异性参数分析,能实现对大炮检距地震反射资料的有效处理。
本文针对VTI介质中PP波和PSV波各向异性参数分析过程,提出充分利用多波数据信息丰富的优势,采用预估—分析—提取的策略限定各向异性参数的范围,简化各向异性参数分析的流程,达到提高多波各向异性参数提取精度和分析效率的目的;并通过优化参数谱,增强各向异性参数谱的分辨率。
通过理论模型和实际资料验证了方法的正确性和可行性。
【期刊名称】石油地球物理勘探【年(卷),期】2014(049)003【总页数】5【关键词】关键词 VTI介质参数分析多波1 引言VTI介质是指具有垂向对称轴的横向各向同性介质,是当前研究的热点。
在地面多波多分量地震勘探中,PP波和PSV波联合勘探成为一种趋势。
关于VTI 介质PP波和PSV波各向异性参数分析研究,早期主要集中在对各向异性介质中的时距曲线方程进行推导和简化[1~4],其中,Alkhalifah[3]通过推导得到的PP波时距曲线方程沿用至今;Li等[4]通过引入χeff参数对PSV波时距曲线进行简化,为VTI介质中进行转换波参数分析提供了新的思路。
近几年尽管许多学者都对各向异性技术进行了系统的研究,也取得了丰硕的成果[5~17],但是对于越来越精细的地震勘探,人们仍然期待更高精度和效率的各向异性参数分析技术。
本文改进了VTI介质中多波各向异性参数分析的方法,在PP波和PSV波参数分析时,充分发挥多波信息量丰富的优势,利用PP波信息,采用预估—分析—提取的策略,简化PSV波参数分析过程,增强其对地层的分辨能力,从而达到提高多波各向异性参数分析效率的目的。
2 参数分析方法原理及改进措施在VTI各向异性介质中,无论是纵波还是转换波,其时距曲线都是非双曲线,要比各向同性介质的时距曲线复杂得多。
2020年4月第55卷 第2期 辽宁省沈阳市黄河北大街280号中国地质调查局沈阳地质调查中心,163712。
Email:moonliyan@sohu.com本文于2019年5月10日收到,最终修改稿于2020年1月2日收到。
本项研究受国家自然科学青年基金项目“油气勘探中宽方位OVT处理方法及其应用研究”(41804138)和中国地质调查局项目“松辽盆地北部页岩油战略调查与评价”(DD20190114)联合资助。
·处理技术·文章编号:1000 7210(2020)02 0292 10宽方位地震资料方位各向异性处理———以大庆长垣三分量地震数据为例张丽艳 ① 李 昂① 刘建颖② 杨建国① 陈志德②(①中国地质调查局沈阳地质调查中心,辽宁沈阳110000;②中国石油大庆油田勘探开发研究院,黑龙江大庆163712)张丽艳,李昂,刘建颖,杨建国,陈志德.宽方位地震资料方位各向异性处理———以大庆长垣三分量地震数据为例.石油地球物理勘探,2020,55(2):292 301,310.摘要 地震纵波在HTI介质中传播时,沿垂直于裂缝方向的速度小于沿平行于裂缝方向,且反射强度降低和频率变低;横波在HTI介质中传播时分裂成快、慢横波。
以大庆长垣喇嘛甸地区的宽方位三分量地震资料为基础,分析了方位各向异性的表现特征、影响及不同方位角范围内属性的变化,并根据这些性质的变化利用椭圆拟合对该地区的各向异性的方位和强度进行定量预测,与利用横波分裂预测的结果一致。
为了消除方位各向异性对地震成像分辨率的影响,利用椭圆拟合准确求出了HTI介质的方位速度,进行方位速度动校正,消除了方位各向异性引起的快慢波时差,提高了资料的成像质量。
关键词 宽方位 各向异性方向 各向异性强度 快慢波时差 方位速度 椭圆拟合中图分类号:P631 文献标识码:A 犱狅犻:10.13810/j.cnki.issn.1000 7210.2020.02.0080 引言随着常规油气产量的不断下降,勘探理念的不断转变和创新,页岩油、页岩气勘探、开发逐渐成为油田增储的热点领域[1]。
正交各向异性介质中的剪张源震源机制与矩张量特征唐杰;温雷;王浩;张文征【摘要】水力压裂区在整体上呈现各向异性特征,深入研究和分析水力压裂区的各向异性参数对震源机制和矩张量的影响成为必要.本文研究了干燥和完全饱水垂直正交裂缝介质的地震波响应特性,通过各向异性介质参数和源张量获得各向异性介质中的微地震矩张量,分析了震源区各向异性对双力偶分量(DC)、补偿线性偶极子分量(CLVD)和各向同性分量(ISO)的影响,并对比了裂缝及其所合流体的特征参数对地震波响应特性的影响效果.研究结果表明:各向异性介质中由于剪张破裂产生的非DC分量依赖于破裂方式、各向异性介质的类型、裂缝参数及断面方向等,干燥裂纹介质中的ISO分量较高,水饱和裂纹介质中产生较高的CLVD分量、而ISO较小;各向同性与各向异性介质中剪张源的远场P波辐射花样存在差异,震源特征和介质各向异性会对地震波走时、极性和振幅产生明显影响.探究水力压裂区各向异性参数对震源机制的影响效果对微地震震源描述具有重要意义.【期刊名称】《石油地球物理勘探》【年(卷),期】2018(053)006【总页数】9页(P1247-1255)【关键词】正交各向异性;微地震;震源机制;地震矩张量;剪张源【作者】唐杰;温雷;王浩;张文征【作者单位】中国石油大学(华东)地球科学与技术学院,山东青岛266580;中国石油大学(华东)地球科学与技术学院,山东青岛266580;中国石油大学(华东)地球科学与技术学院,山东青岛266580;中国石油大学(华东)地球科学与技术学院,山东青岛266580【正文语种】中文【中图分类】P6311 引言地壳和上地幔广泛存在着方位各向异性,产生原因包括沉积物分层、微裂纹定向排列以及岩石矿物的结构特征等[1-3]。
各向异性会对震源信号及地震波的传播特征产生明显影响,以往的研究主要针对各向异性介质中的波传播特征分析,而各向异性对于震源信号的影响效果需做进一步分析。
利用Hartley变换模拟各向异性介质地震波场
周辉;何樵登
【期刊名称】《石油地球物理勘探》
【年(卷),期】1995(030)005
【摘要】伪谱法在地震波场的模拟中具有很重要的作用。
然而,基于傅氏变换的伪谱法运算速度慢,占用计算机内存较多。
Hartley变换是一种快速的变换方法,能克服傅氏变换的上述缺陷。
【总页数】10页(P593-601,614)
【作者】周辉;何樵登
【作者单位】不详;不详
【正文语种】中文
【中图分类】P631.443
【相关文献】
1.二维各向异性介质中地震波场的高阶同位网格有限差分模拟 [J], 祝贺君;张伟;陈晓非
2.运用Hartley变换模拟地震波在正交各向异性介质中的传播 [J], 徐天吉;李显贵;范菊芬;程冰洁
3.基于双孔模型的含流体各向异性介质地震波场数值模拟 [J], 阮韵淇;李世凯;李天;晓涛
4.基于GPU和CPU协同并行的三维各向异性介质地震波场正演模拟 [J], 刘春成;顾汉明;陈宝书;焦振华;马凯;蔡志成;张立
5.各向异性介质地震波场的优化褶积微分算子法数值模拟 [J], 朱童;李小凡;张美根;龙桂华
因版权原因,仅展示原文概要,查看原文内容请购买。
利用地震波在双相介质中的衰减特性检测油气张会星;何兵寿;姜效典;张锦昌【期刊名称】《石油地球物理勘探》【年(卷),期】2010(045)003【摘要】在平面波假设条件下,本文通过求解波动方程推导了双相介质中地震波衰减系数计算公式,并分析了衰减系数随频率的变化关系.结果表明,在常规地震勘探的频率范围内,地震波在双相介质中传播时,存在介质吸收效应,且衰减系数近似与圆频率成正比.衰减系数同时受地震波频率和耗散系数的影响,表现为频率越高,衰减系数越大,耗散系数越大.衰减系数随频率的变化越大,表明地震波在穿过双相介质时其不同频率成分对应的能量要发生变化.低频能量衰减小,高频能量衰减大,在频谱上表现为低频能量相对增强,高频能量相对减弱.据此提出了一种利用地震资料的低、高频段信息反演地下介质双相特性或含油气性的方法.理论模型与实际资料的反演结果均验证了文中方法是可行的.【总页数】7页(P343-349)【作者】张会星;何兵寿;姜效典;张锦昌【作者单位】中国海洋大学海底科学与探测技术教育部重点实验室,山东青岛266100;中国矿业大学煤炭资源与安全开采国家重点实验室,北京100083;中国海洋大学海底科学与探测技术教育部重点实验室,山东青岛266100;中国海洋大学海底科学与探测技术教育部重点实验室,山东青岛266100;中国海洋大学海底科学与探测技术教育部重点实验室,山东青岛266100【正文语种】中文【相关文献】1.双相介质油气检测技术在龙西—杏西地区含油性预测中的应用 [J], 李永义2.基于双相介质的地震波衰减特性油气检测方法 [J], 汪佳蓓;黄捍东;郭飞;;3.基于双相介质的地震波衰减特性油气检测方法 [J], 汪佳蓓;黄捍东;郭飞4.利用宁夏大峰矿特大爆破资料研究地震波的衰减特性 [J], 金春华;赵卫明;盛菊琴;马禾青5.利用时间分裂的错格伪谱法模拟地震波在基于改进BISQ模型的双相介质中的传播 [J], 兰慧田;刘财;郭智奇因版权原因,仅展示原文概要,查看原文内容请购买。
地震各向异性研究进展一、引言2011年第81届SEG 年会“各向异性”专题涉及文章32篇(含张贴报告),讨论的内容主要包括各向异性介质波动方程解耦、地震波在各向异性介质中的传播、各向异性介质偏移和成像方法研究、各向异性介质时差速度反演与各向异性参数建模方法、各向异性弹性参数反演与裂缝预测、近地表各向异性分析等诸多方面。
在本文中,我们对专题论文涉及的主要内容进行了归类,总结了各向异性介质地震波传播、成像与反演方面的研究进展。
二、各向异性介质伪P 波正演和逆时偏移地震波在各向异性介质中传播是一种弹性现象,由于弹性各向异性介质波动方程十分复杂,所以大多数各向异性正演和偏移方法通常使用伪纵波方程。
Chu 等(2011)给出了使用伪差分算子的泰勒级数展开近似解耦得到纯P 波方程,级数展开产生了一个新的纯纵波各向异性方程如式(1)所示,对于一定范围的波数和各向异性参数都是收敛的。
4222222222222ˆˆˆˆˆˆ()()()px x y pz z pz px px x y z v k k v k v v v k k k ωω⎡⎤-=-++--+⎣⎦(1) 使用有限差分法求解导出的方程(1),并用纯P 波方程进行数值模拟,其结果与伪纵波方程的结果十分接近,精度较高,并且能保证完全不受横波的影响。
图1为BP 倾斜各向异性(TTI )模型的一部分,(a)为速度模型,(b )、(c )、(d )分别为各向异性参数模型,其中各向异性参数θ变化很大。
模拟结果如图2所示,其中(a )为伪声波波动方程(0sz v =)的模拟结果,(b )为伪声波波动方程(42()3v v sz pz εδ=-)的模拟结果;(c )为纯声波波动方程模拟结果;(d )为图(c )与图(b )叠加结果。
从图2中可以看出当横波速度为零时,用伪声波波动方程产生的结果不稳定,引入非零S 波速度后的结果稳定,但是用该方程得到的结果包含了震源产生的S 波和P-SV 转换波。