VTI介质qP波方程高精度有限差分算子
高阶旋转交错网格有限差分方法模拟TTI介质中横波分裂
高阶旋转交错网格有限差分方法模拟TTI介质中横波分裂李敏;刘洋【期刊名称】《物探与化探》【年(卷),期】2012(36)6【摘要】笔者给出了一种能够模拟弹性波在任意各向异性介质中传播的二维三分量高阶有限差分算法.相对于常规交错网格有限差分方法,旋转交错网格有限差分方法在介质具有强差异性时能更精确地模拟地震波的传播,避免常规交错网格中因对弹性系数进行插值而带来的误差.采用高阶旋转交错网格有限差分方法模拟并分析了零偏移距横波分裂现象随裂缝介质方位角和倾角变化的响应特征.结果表明:结合完全匹配层(PML)吸收边界条件的高阶旋转交错网格有限差分方法能获得高精度的地震波场模拟数据,并且在边界具有良好的吸收效果；横波分裂现象主要受裂缝走向与波的极化方向之间的夹角影响,受裂缝倾角影响较小,且快慢横波的能量也跟裂缝走向与波极化方向间的夹角有关.具有倾斜对称轴的横向各向同性(TTI)介质倾角的变化可能会导致记录中波到达时的变化,影响快慢横波的时差.利用横波分裂的能量分布和方位各向异性特征,可以帮助检测裂缝的方位角和倾角.横波在多层TTI介质中传播时会发生多次分裂的现象.%A high-order rotated staggered grid scheme ( RSG) has been implemented to simulate the shear-wave splitting in tilted transversely isotropic (TTI) media. The high-order RSG can simulate wave propagation in media that contain high-contrast discontinuities like cracks more precisely than the standard staggered grid scheme (SSG) by avoiding the unstableness of the staggered grid scheme (SSG). The authors conducted a study of zero-offset S-wave splitting withthe high-order RSG. The S-wave splitting study was mainly focused on fractured media which, on the scale of seismic wavelength, could be regarded as transversely isotropic (TI) media. The results of numerical modeling show that the high-order RSG scheme can be used to simulate waves' propagation in general anisotropic media. The perfect matched layer (PML) absorbing boundary condition combined with the high order RSG scheme can well attenuate reflections from the artificial boundary. The S-wave splitting is mainly affected by the angle between polarization direction of incoming wave and strike of the TTI media, and the energy of fast and slow shear waves is also associated with this angle. The dipping angle of ITI media may affect time lag between the fast and slow waves, which may result in variation of arrival time of waves from the same interface. Thus, the analysis of energy distribution of the fast and slow waves and the variation of arrival time may help detect the strike and dipping angle of the fracture. Besides, when propagating in the media that contain more than one layer of TTI media, the S-wave splitting will occur more than once.【总页数】7页(P934-940)【作者】李敏;刘洋【作者单位】中国石油大学油气资源与探测国家重点实验室,北京102249;中国石油大学CNPC物探重点实验室,北京102249;中国石油大学油气资源与探测国家重点实验室,北京102249;中国石油大学CNPC物探重点实验室,北京102249【正文语种】中文【中图分类】P631.4【相关文献】1.黏弹TTI介质旋转交错网格微地震波场模拟 [J], 姚振岸;孙成禹;谢俊法;唐杰2.黏弹TTI介质中旋转交错网格高阶有限差分数值模拟 [J], 严红勇;刘洋3.旋转交错网格在横波分裂和再分裂模拟与分析中的应用 [J], 张建利;刘志斌;周超;田小波;李维新;王赟4.基于卷积完全匹配层的旋转交错网格高阶差分法模拟弹性波传播 [J], 冯德山;王向宇5.基于CPML-RML组合边界条件粘弹TTI介质旋转交错网格有限差分正演模拟[J], 张奎涛;顾汉明;刘少勇;刘春成;陈宝书;张立;肖逸飞因版权原因，仅展示原文概要，查看原文内容请购买。
TTI介质Lebedev网格高阶有限差分正演模拟及波型分离
TTI介质Lebedev网格高阶有限差分正演模拟及波型分离黄金强;李振春;黄建平;张金淼;孙文博【期刊名称】《石油地球物理勘探》【年(卷),期】2017(052)005【摘要】针对TTI介质交错网格高阶有限差分正演模拟精度低及矢量波型分离不彻底的问题,发展了一种TTI介质Lebedev网格高阶有限差分正演模拟方法及矢量波型分离一体化处理流程.首先从一阶偏导数弹性波速度-应力方程出发,构建了TTI 介质Lebedev网格高阶差分波场递推格式;然后借助Low rank分解策略处理由极化矢量构成的空间-波数域分离算子,进而对速度分量实施纵、横波型分离;最终实现了Lebedev网格高阶差分方案与纵、横波分离模拟思路的有机结合,形成了复杂TTI介质高精度正演模拟与波型分离一体化处理流程.在实现方法的基础上,对均匀介质、层状介质以及复杂BP2007模型进行了正演及分离测试.计算结果表明:选取Lebedev网格高阶有限差分能够减小交错网格剖分方式带来的插值误差,从而获得高精度矢量波场;其次,通过应用Low-rank分解波型分离方法能够实现矢量波彻底分离,得到完全解耦的纵、横波场;文中方法对非均匀复杂模型具有良好的适应能力.【总页数】13页(P915-927)【作者】黄金强;李振春;黄建平;张金淼;孙文博【作者单位】中国石油大学(华东)地球科学与技术学院地球物理系,山东青岛266580;中国石油大学(华东)地球科学与技术学院地球物理系,山东青岛266580;中国石油大学(华东)地球科学与技术学院地球物理系,山东青岛266580;中海油研究总院,北京100028;中海油研究总院,北京100028【正文语种】中文【中图分类】P631【相关文献】1.基于Lebedev网格的TTI介质二维三分量正演模拟 [J], 刘东洋;彭苏萍;师素珍;赵太郎2.二维三分量TTI介质Lebedev网格正演模拟 [J], 李娜;黄建平;李振春;李庆洋;郭振波;田坤3.Lebedev网格改进差分系数TTI介质正演模拟方法研究 [J], 李娜;黄建平;李振春;田坤;李庆洋4.变网格波动方程高阶有限差分正演模拟研究 [J], 段宏凯;张晓丹;佘翼翀;谢宝林;杨白雪5.基于CPML-RML组合边界条件粘弹TTI介质旋转交错网格有限差分正演模拟[J], 张奎涛;顾汉明;刘少勇;刘春成;陈宝书;张立;肖逸飞因版权原因，仅展示原文概要，查看原文内容请购买。
VTI介质中多波多分量数值模拟及PML吸收边界条件
VTI介质中多波多分量数值模拟及PML吸收边界条件张为彪;洪宇;郑洁;夏弋峻;郭飞;邹清文【摘要】为了充分认识多波多分量地震勘探技术中各种波及其分量的传播特点,以便深入挖掘来自地下复杂介质信号中的有效信息,解决实际勘探中遇到的复杂问题,选择了一种特殊的各向异性介质——具有垂直对称轴的横向各向同性介质(VTI)中传播的波及其分量作为研究对象,对其波场进行了数值模拟,分析总结了其传播特点.为了解决边界反射的干扰问题,构建并应用了相应的PML吸收边界条件,模拟结果表明,吸收效果较好.【期刊名称】《石家庄经济学院学报》【年(卷),期】2017(040)005【总页数】6页(P7-12)【关键词】多波;多分量;各向异性;吸收边界条件;数值模拟【作者】张为彪;洪宇;郑洁;夏弋峻;郭飞;邹清文【作者单位】中海石油(中国)有限公司深圳分公司,广东深圳 518054;中海石油(中国)有限公司深圳分公司,广东深圳 518054;中海石油(中国)有限公司深圳分公司,广东深圳 518054;中海油能源发展股份有限公司工程技术分公司,广东深圳518054;中海石油(中国)有限公司深圳分公司,广东深圳 518054;中海石油(中国)有限公司深圳分公司,广东深圳 518054【正文语种】中文【中图分类】TE132随着油气勘探程度的增加，勘探难度不断加大，出现了很多具有挑战性的地质勘探问题，目前最为常用、理论和技术发展也最为成熟的地震勘探方法—纵波勘探已无法完全解决这些问题，而多波多分量地震勘探技术在解决这些问题上有其独特的优势，如：消除模糊带、识别真假亮点、探测裂缝、获得与横波速度有关的各种岩性参数等。
波场数值模拟是人们理解波在不同介质中传播规律，认识和正确解释采集波信号携带的重要信息的重要手段[1-2]。
为了理解多波多分量在地下介质中的传播特性，需对其波场进行数值模拟，模拟时，为了避免人工边界产生的反射干扰介质内的波场，必须引入吸收边界条件来处理边界。
含煤地层各向异性介质有限差分数值模拟
含煤地层各向异性介质有限差分数值模拟钱进;崔若飞;陈同俊【期刊名称】《煤田地质与勘探》【年(卷),期】2010(038)002【摘要】煤储层是一种双重孔隙介质,其孔隙由基质孔隙和裂隙组成,并具有自身独特的割理系统.传统的煤田地震勘探将煤层假设为单一的各向同性介质,没有考虑其中裂隙对煤层的影响.利用交错网格高阶有限差分算法,分别对不同厚度VTI和HTI 型煤层的弹性波传播特征进行数值模拟,获得波场快照和含煤层状介质的地震模拟记录.从波场快照中可以看出,弹性波在VTI和HTI型煤层中传播存在明显差异;煤层厚度变化对波场的敏感性要强于煤层类型的变化.【总页数】5页(P63-67)【作者】钱进;崔若飞;陈同俊【作者单位】中国矿业大学资源与地球科学学院,江苏,徐州,221116;煤层气资源与成藏过程教育部重点实验室,江苏,徐州,221116;中国矿业大学资源与地球科学学院,江苏,徐州,221116;煤层气资源与成藏过程教育部重点实验室,江苏,徐州,221116;中国矿业大学资源与地球科学学院,江苏,徐州,221116;煤层气资源与成藏过程教育部重点实验室,江苏,徐州,221116【正文语种】中文【中图分类】P631【相关文献】1.各向异性介质弹性波高阶交错网格有限差分模拟 [J], 霍凤斌;李振鹏;徐发;张涛;2.三维各向异性介质中弹性波方程交错网格高阶有限差分法数值模拟 [J], 裴正林3.正交各向异性介质中地震记录的有限差分数值模拟 [J], 邵治龙;贺振华4.双相各向异性介质中偶数阶精度有限差分数值模拟 [J], 刘洋;魏修成5.各向异性介质Low-rank有限差分法纯qP波叠前平面波最小二乘逆时偏移 [J], 黄金强;李振春因版权原因，仅展示原文概要，查看原文内容请购买。
二维TTI介质的纯P波波动方程数值模拟
二维TTI介质的纯P波波动方程数值模拟张千祥;王德利;周进举【摘要】声波各向异性数值模拟对地震数据处理和解释起着重要的作用.基于Tsvankin提出的精确色散关系,通过平方根近似,在时间-波数域中推导出二维TTI 介质纯P波声波波动方程,并利用快速展开法(Rapid Expansion Method,REM)进行了数值模拟.与传统的有限差分法求解二维TTI介质耦合方程和傅里叶有限差分法在时间上进行波场外推相比,该方法的模拟结果精度更高,计算速度更快,并且成功去除横波分量.【期刊名称】《石油物探》【年(卷),期】2015(054)005【总页数】8页(P485-492)【关键词】声波各向异性数值模拟;纯P波声波方程;快速展开法;有限差分法;傅里叶有限差分法【作者】张千祥;王德利;周进举【作者单位】吉林大学地球探测科学与技术学院,吉林长春130026;吉林大学地球探测科学与技术学院,吉林长春130026;吉林大学地球探测科学与技术学院,吉林长春130026【正文语种】中文【中图分类】P631地震数值模拟是地震勘探方法研究的前提和基础,在地震勘探和地震学的各项研究及生产工作中都扮演着重要的角色[1]。
常用的地震波场数值模拟方法主要有几何射线法、波动方程法和积分方程法[2]。
波动方程模拟方法中的有限差分法由于计算速度快、占用计算机内存小而被广泛应用。
很早时候各国地球物理学家就对各向同性介质和各向异性介质弹性波地震数值模拟进行了深入研究。
近年来,周进举等[3]利用高阶旋转网格有限差分法研究了复杂介质下弹性波数值模拟。
在对地下的各向异性介质进行弹性波数值模拟时,由于弹性波方程复杂,各向异性参数多,导致模拟计算量大,耗时长,增加了弹性波偏移和干涉的难度。
为了解决这些问题,我们采用声波各向同性近似理论,通过设定弹性波中的横波速度为零来简化计算参量,在保证模拟精度的条件下提高计算效率。
然而,由于地下介质的不均匀性,这种各向同性的声学假设常常是不恰当的。
高精度有限差分瑞雷面波模拟及频散特征提取
文章标题：深度解析高精度有限差分瑞雷面波模拟及频散特征提取一、介绍高精度有限差分瑞雷面波模拟及频散特征提取是当今地球物理领域的热门研究课题之一。
通过对地球内部的雷达波传播特性进行模拟和分析，可以更好地理解地球的结构和性质，为地质勘探、地震监测等领域提供重要支持。
本文将从深度和广度的角度，对高精度有限差分瑞雷面波模拟及频散特征提取进行全面解析，并结合个人观点进行探讨。
二、原理及模拟技术1. 高精度有限差分瑞雷面波模拟简介在地球物理探测中，雷达波在地下介质中的传播特性对地下结构的识别和勘探具有重要意义。
高精度有限差分瑞雷面波模拟是一种常用的地球物理探测技术，通过数值模拟地下介质对雷达波的响应，可以获得地下结构的相关信息，包括速度、密度等。
2. 频散特征提取方法频散是指不同频率的波在介质中传播速度不同，导致波形畸变的现象。
在雷达波传播过程中，频散特征的提取对于识别地下结构、解释波形数据具有重要作用。
通过合适的算法和分析方法，可以有效提取频散特征，从而更准确地理解地下介质的特性。
三、从简到繁的探讨1. 我们可以从传统的有限差分模拟方法入手，介绍基本的数值模拟原理和计算技术，以及其在地球物理勘探中的应用。
2. 深入讨论高精度有限差分模拟技术的优化和改进，包括并行计算、高性能计算等方法，以提高模拟结果的精度和准确度。
3. 结合实际案例，探讨高精度有限差分模拟在地质勘探、地下水资源调查等领域的应用和效果，以及对相关领域的意义和挑战。
四、频散特征提取分析1. 介绍常见的频散特征提取方法，包括谱分析、小波变换等技术，分析其应用范围和局限性。
2. 深入探讨频散特征在地球物理学中的重要性，以及在雷达波数据解释和地下结构识别中的作用。
3. 结合实际案例，剖析频散特征提取在地震监测、地下勘探等领域的应用效果和前景。
五、总结与展望通过对高精度有限差分瑞雷面波模拟及频散特征提取的全面分析，我们可以更好地理解这一研究领域的理论基础和应用技术，为相关领域的研究和实践提供重要参考。
基于VTI介质的多波速度分析
基于VTI介质的多波速度分析摘要：在各向异性介质中，特别是大偏移距道集，纵波的时距方程不再满足双曲规律，因此对纵波数据的处理应采用非双曲时差速度分析方法；而对转换波来说，即使在水平各向同性介质中转换波的射线路径也是不对称的。
本文采用横向各向同性（VTI）介质，利用非双曲时距方程实现了纵波和转换横波的速度分析与参数反演，并通过模型试算和实际资料处理说明了所用方法的必要性和正确性。
关键词：各向异性；速度分析；时差；纵波；转换波1 引言速度分析技术是地震数据处理的关键技术之一。
常规速度分析是基于各向同性介质的假设，然而在实际勘探中发现地下介质广泛存在各向异性[1]。
各向异性的存在，使得地震波的传播变得更加复杂，最突出的特点是速度的各向异性，此时的速度是传播方向的函数，而不再是一个标量。
研究表明，如果忽略速度各向异性会造成速度提取的不准确，进而影响后续的处理流程。
因而，对于存在各向异性的地震资料需要用各向异性的速度分析方法进行处理。
2 纵波非双曲时差速度分析5 结论（1）对于纵波，在双曲速度分析的基础上引用非双曲时差方程进行速度分析，充分利用短、中、长排列数据提高了速度分析的精度。
（2）对于转换波，引用了高精度的泰勒级数展开方程，在短排列分析的基础上，充分利用中、长排列的信息，设计合理的速度分析流程并通过模型及实际资料试处理取得较好的效果。
（3）通过对比分析，证实了当介质存在各向异性时，各向异性非双曲时差速度分析相对于双曲时差速度分析而言能够提供更加准确的叠加速度，并且能够反演出较准确的各向异性参数。
参考文献[1] Thomsen L. Weak elastic anisotropy[J]. Geophysics，1986，51（10）：1954～1966.[2] Alkhalifah T and Tsvankin I. Velocity analysis for transversely isotropic media[J]. Geophysics，1995，60（5）：1550～1566.[3] Li X Y. Converted-wave moveout analysis revisited：the search for a standard approach[J].73nd annual international SEG meeting，2002，6～11.。
三维双相各向异性介质弹性波方程交错网格高阶有限差分法模拟
三维双相各向异性介质弹性波方程交错网格高阶有限差分法模
拟
裴正林
【期刊名称】《中国石油大学学报（自然科学版）》
【年(卷),期】2006(030)002
【摘要】基于Biot理论,给出了三维双相各向异性介质应力-速度弹性波方程交错网格任意偶阶精度有限差分解法,并对三维双相各向异性(TIH)介质中弹性波场进行了模拟.结果表明,弹性波在三维双相各向异性介质中传播时存在快纵波qP1、快横波qS1、慢横波qS2和慢纵波qP2,并清楚地观测到了横波分裂、横波分裂盲点、波面三分叉等特殊现象.快纵波在固相和流相中的相位相同;而慢纵波在固相和流相中的相位相反,即慢纵波在流相中振幅大,而在固相中的振幅较小.另外,三维双相各向异性介质中弹性波场的快纵波和快横波的耦合关系、波的类型、能量分布和相位等都是在三维空间中变化的.
【总页数】5页(P16-20)
【作者】裴正林
【作者单位】中国石油大学,中国石油天然气集团公司物探重点实验室,北
京,100083
【正文语种】中文
【中图分类】P631.443
【相关文献】
1.三维各向异性介质中弹性波方程交错网格高阶有限差分法数值模拟 [J], 裴正林
2.双相各向异性介质弹性波传播交错网格高阶有限差分法模拟 [J], 裴正林
3.任意起伏地表弹性波方程交错网格高阶有限差分法数值模拟 [J], 裴正林
4.三维各向同性介质弹性波方程交错网格高阶有限差分法模拟 [J], 裴正林
5.任意倾斜各向异性介质中弹性波波场交错网格高阶有限差分法模拟 [J], 裴正林;王尚旭
因版权原因，仅展示原文概要，查看原文内容请购买。
二维波动方程的一种高精度紧致差分方法
（一 ） —２） ＝ ） Ａ （ ， （ ； Ｔ２ ＋ ×
ｅ
，
Ｊ
＋Ｒ
０≤ ， ≤ 一１， Ｏ≤ ｎ≤ Ⅳ 一 １
用 ｒ表示时 间步 长 ， ｈ表 示空 间步 长 ， １给 出 了 当 ７＝ 表 －
（） ８
ｅ ＝ｏ
ｈ ［ （ Ｌ 曰
Ａ
（
３ 数值 验 证
对于式 （ ）一（ ） 令 （ ｙ ￣ ， ｉ（ ）ｉ（ｒ） （ １ ４， ， ）＝ ２ ｒｎ ｓ Ｔ ， ， ｒ ｓ ｎ ｙ
Ｙ ）＝０ｇ ， ，） ， 问题 的精确解 为 ／ ， ，）＝ ｉ（ 订 ） ，（ Ｙ ｔ ＝０ 则 ２ Ｙ ｔ ｓ ｔ （ ｎ ｓ （Ｔ）ｉ（ｒ） ｉ ＇ ｓ 叮 。数值 实验计算是用 Ｆ ￣ａ ７语言进行编程 ｎＩ ｎ ｙ Ｘ ｏ ｒｎ７
ｈ） 。 阶精度的数值 解。孙 志 忠 提 出 了求 解二 维波 动方 程 的 高精度交替方 向隐式 方法 ， 并且是无条件稳定 的。有关这方 面 最新 的一 些 工作 可参 见 文献 ［ ６—８ 。本文 在 此工 作基 础 之 ］
上， 利用 Ｒｃ ａ ｓｎ外推法进 一步 提 高计 算精 度 ， ｉ ｒｏ ｈｄ 最终 可得 到
ｏ≤ √≤ 一１
４ ｔ ｈ ，＝１时刻 ， 本文格 式在不 同网格步 长下误 差 的 、 、
范数 ， 以及与四阶 Ａ Ｉ Ｄ 格式 计 算结 果 的 比较 。Ｌ 范 数定 义 ２
厂］ｉ 面＝广——一 『
ｄ］ ＝ ｏ 。 ｅ
＼ ｅ ＝０
。 ≤ Ｍ
一
１
层的。即每一次时间推进都需要知道前 两个 时间步 的值 ， ０ 第
求解该 问题 精度为 Ｏ（ ＋ 。 的数值解 。 ｈ）
三维波动方程有限差分正演方法
三维波动方程有限差分正演方法
首先，将三维波动方程进行离散化处理，将连续的空间和时间域离散
为网格点的集合。
在空间域中，将地下介质划分为多个均匀网格，每个网
格点上对应一个地震波场的值。
在时间轴上，将时间域离散为多个时间步长，每个时间步长表示一个时刻。
然后，利用有限差分算子将三维波动方程转化为离散形式。
常用的有
限差分格式有常系数差分格式和变系数差分格式。
常系数差分格式适用于
各向同性介质，而变系数差分格式适用于各向异性介质。
在三维波动方程有限差分正演方法中，需要考虑各向异性介质的模型。
各向异性介质与各向同性介质不同，其波速和波阻抗在不同方向上具有不
同的值。
因此，在离散化过程中，需要引入特殊的差分算子来考虑各向异
性的影响。
最后，利用迭代求解的方法，按时间步长依次求解离散化的波动方程。
利用差分算子和初始条件，根据时空的变化逐步更新波场的数值，并计算
出每个网格点的波动值。
通过不断迭代求解，最终可以得到地震波在地下
的传播和反射结果。
三维波动方程有限差分正演方法在地震勘探中具有重要的应用价值。
它可以模拟地震波在地下介质中的传播和反射过程，帮助研究人员了解地
下结构、油气储层等。
同时，它也可以用来研究地震波在各向异性介质中
的传播规律，提高地震勘探的准确性和效率。
总之，三维波动方程有限差分正演方法是一种常用的数值模拟方法，
可以用来模拟地震波在地下介质中的传播和反射。
它具有广泛的应用价值，对地震勘探、岩石物理、地质测井等领域具有重要的意义。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
题， 会 降低 地 震 波场 模 拟 的精 度 与 分 辨 率 ． 为 了克服 常 规 有 限 差 分 算 子 的 数 值 频 散 ， 本 文 针 时 ＶＴＩ 介 质 地 震 波 数 值 模 拟 问题 ， 构 造 了频 率 空 间 域 ｑ Ｐ波 波 动 方程 高精 度 有 限 差 分 优 化 算 子 ， 根 据 最 优 化 理 论 中 高 斯一 牛 顿 法 确 定 了 高 精 度 有 限 差 分 算 子 的 优 化 系数 ． 利 用 常 规 差 分 算 子 和 高精 度 优 化 差 分 算 子 时 归 一 化 相 速 度 的 频 散 关 系精 度进 行 了 时 比 分 析， 并 对 均 匀 各 向 同性 介 质 和 均 匀 ＶＴ Ｉ 介 质 中的 ｑ Ｐ 波地 震 波 场进 行 了有 限 差 分 数 值 模 拟 ， 通 过 频 散 关 系精 度 分 析 和 波场 数 值 模 拟 结 果表 明 ： 有限差分优化算子具有较 高的波场 数值模 拟精 度， 有 效 压 制 了传 统 有 限 差 分 算 子 数 值 模 拟 中 的数 值 频 散 现 象 ， 提 高 了有 限 差 分 算 子 精 度 ， 为 ＶＴ Ｉ 介 质 频 率一 空 间域 ｑ Ｐ 波 正 演模 拟 奠 定 了基 础 、 关键词 ＶＴ Ｉ 介 质， 有 限 差 分算 子 ， 优 化 系数 ， 数 值 频散 Ｐ ６ ３ １ 文献 标 识 码 Ａ 文 章 编 号 １ ０ ０ ４ — ２ ９ ０ ３ （ ２ ０ ０ ７ ） ０ ３ — ０ ８ ９ ６ — ０ ９
ＶＴＩ 介质 ｑ Ｐ波 方 程 高精 度 有 限差分 算 子
吴国忱 ， 梁 锴
（ 中国石油大学（ 华 东） 地 球 资 源 与 信 息学 院 ， 东营 ２ ５ ７ ０ ６ １ ）
摘
要
波 动 方 程 有 限 差 分 法 是 一 种 使 用 广 泛 的 地 震 波 数 值 模 拟 方 法． 但 是 有 限 差 分 法 本 身 固 有 存 在 着 数 值 频散 问
维普资讯
第 ２ ２卷
第 ３期
地
球
物
理
学
进
展
Ｖｏ １ ． ２ ２ Ｎｏ ． ３
２ ０ ０ ７年 ６月 （ 页码 ： ８ ９ ６ ～９ ０ ４ ）
ＰＲ０Ｇ ＲＥＳＳ Ｉ Ｎ ＧＥ０Ｐ Ｈ ＹＳＩ ＣＳ
Ｊ ｕ ｎ ｅ ２ ０ ０ ７
中 图分 类 号
Ｈｉ ｇ ｈ ｐｒ ｅ ｃ ｉ ｓ ｉ ｏ ｎ ｆ ｉ ｎｉ ｔ ｅ ｄｉ ｆ ｆ ｅ ｒ ｅ ｎ ｃ ｅ ｏ ｐ ｅ ｒ ａ ｔ ｏ ｒ ｓ ｆ ｏ ｒ
ｑＶＴＩ ｍｅ ｄｉ ａ
Ｗ Ｕ ＧＵ Ｏ — ｔ ｈ ｅ ｎ。 ＬＩ Ａ Ｎ Ｇ Ｋａ ｉ
Ａｂ ｓ ｔ ｒ ａ ｃ ｔ Ｆｉ ｎ ｉ ｔ ｅ － ｄｉ ｆｅ ｒ ｅ ｎ ｃ ｅ ｍｅ ｔ ｈ ｏｄ ｆ ｏｒ ｗａ ｖｅ ｅ ｑｕ ａ ｔ ｉ ｏｎ ｉ ｓ ｗｉ ｄｅ ｌ ｙ ｉ ｍ ｐｌ ｅ ｍｅ ｎｔ ｅ ｄ ｉ ｎ ｓ ｅ ｉ ｓ ｍｉ ｃ ｗａｖ ｅ ｆ ｉ ｅ ｌ ｄ ｎｕｍ ｅ ｒ ｉ ｃ ａｌ ｓ ｉ ｍ ｕｌ ａ — ｔ ｉ ｏ ｎ． Ｆｉ ｎｉ ｔ ｅ － ｄ ｉ ｆ ｆ ｅ ｒ ｅ ｎ ｃ ｅ ｍｅ ｔ ｈｏ ｄ ｈ ａ ｓ ｉ ｎｈ ｅ ｒ ｅ ｎ ｔ ｎｕｍ ｅ ｒ ｉ ｃ ａ ｌ ｄｉ ｓ ｐｅ ｒ ｓ ｉ ｏｎ。 ｗｈｉ ｃ ｈ ｒ ｅ ｄｕ ｃ ｅ ｓ ｔ ｈｅ ａ ｃ ｃ ｕｒ ａ ｃ ｙ ａ ｎ ｄ ｒ ｅ ｓ ｏｌ ｕ ｔ ｉ ｏｎ ｏｆ ｓ ｅ ｉ ｓ ｍｉ ｃ ｗａ ｖ ｅ ｆ ｉ ｅ ｌ ｄ ｓ ｉ ｍ ｕｌ ａ ｔ ｉ ｏ ｎ． Ｉ ｎ ｏ ｒ ｄｅ ｒ ｔ ｏ ｄｅ ｃ ｒ ｅ ａ ｓ ｅ ｔ ｈｅ ｎｕ ｍｅ ｒ ｉ ｃ ａ ｌ ｄｉ ｓｐ ｅ ｒ ｓ ｉ ｏｎ ｏｆ ｃ ｏｎ ｖｅ ｎ ｔ ｉ ｏｎａ ｌ ｆ ｉ ｎｉ ｔ ｅ － ｄ ｉ ｆｅ ｒ ｅ ｎｃ ｅ ｏｐ ｅ ｒ ａｔ ｏｒ ｓ ， ｗｅ ｐ ｒ ｅ ｓ ｅ ｎｔ ｈ ｉ ｇ ｈ ｐｒ ｅ ｃ ｉ ｓ ｉ ｏｎ ｆ ｉ ｎ ｉ ｔ ｅ － ｄｉ ｆ ｆ ｅ ｒ ｅ ｎ ｃ ｅ ｏｐ ｅ ｒ ａｔ ｏｒ ｓ ｏｆ ｑ Ｐ ｗａ ｖ ｅ ｅ ｑｕ ａ ｔ ｉ ｏ ｎ ｉ ｎ ｆ ｒ ｅ ｑｕ ｅ ｎ ｃ ｙ ’ 。ｓ ｐ ａ ｃ ｅ ｄｏ ｍａ ｉ ｎ ｆ ｏ ｒ ｓ ｅ ｉ ｓ ｍｉ ｃ ｗａ ｖ ｅ — — ｆ ｉ ｅ ｌ ｄ ｎｕ ｍｅ ｒ ｉ ｃ ａ ｌ ｓ ｉ ｍｕ ｌ ａ ｔ ｉ ｏ ｎ ｉ ｎ Ｖ ＴＩ ｍｅ ｄｉ ａ，ａｎ ｄ ｏ ｂｔ ａ ｉ ｎ ｔ ｈｅ ｏ ｐｔ ｉ ｍａ ｌ ｃ ｏ ｅ ｆ ｆ ｉ ｃ ｉ ｅｎ ｔ ｓ ｏ ｆ ｈｉ ｇｈ ｐｒ ｅ ｃ ｉ ｓ ｉ ｏｎ ｆ ｉ ｎ ｉ ｔ ｅ －ｄ ｉ ｆ ｆ ｅ ｒ ｅ ｎｃ ｅ ｏｐｅ ｒ ａ — ｔ ｏｒ ｓ ａ ｃ ｃ ｏｒ ｄｉ ｎ ｇ ｔ ｏ ｔ ｈｅ Ｇａｕ ｓ ｓ — Ｎｅ ｗｔ ｏ ｎ ｍｅ ｔ ｈ ｏ ｄ ｉ ｎ ｏｐｔ ｉ ｍｉ ｚ ａ ｔ ｉ ｏ ｎ ｔ ｈｅ ｏｒ ｙ ．Ｗ ｅ ａ ｎ ａｌ ｙ ｚ ｅ ｄ ｉ ｓ ｐ ｅ ｒ ｓ ｉ ｏ ｎ ｒ ｅ ｌ ａ ｔ ｉ ｏ ｎ ａ ｃ ｃ ｕｒ ａ ｃ ｙ ｏ ｆ ｎｏ ｒ ｍａ ｌ —
（ Ｃ ｏ ｌ ｌ ｅ ｇ ｅ ｏ ｆ Ｇ ｅ ｏ — Ｒ ｅ ｓ ｏ ｕ ｒ ｃ ｅ ｓ ａ ｎ ｄ ｆ ｎ ｆ ｏ ｒ ｍａ ｔ ｉ ｏ ｎ， Ｃｈ ｉ ｎ ａ Ｕ ｎ ｉ ｖ ｅ ｒ ｓ ｉ ｔ ｙ ｏ ｆ Ｐ ｅ ｔ ｒ ｏ ｌ ｅ ｕ ｍ（ Ｅ ａ ｓ ｔ Ｃ ｈ ｉ ｎ ａ ） ，Ｄｏ ｎ ｇ ｙ ｉ ｎ ｇ ２ ５ ７ ０ ６ １ ，Ｃ ｈ ｉ ｎ ａ ）