标量声波波动方程高阶交错网格有限差分法 被引量：58

1

作者 陈可洋 《中国海上油气》 CAS 北大核心 2009年第4期232-236,共5页

为进一步提高标量声波波动方程正演数值模拟精度,减小数值频散和计算量,改善边界吸收效果,提出将高阶交错网格有限差分法直接用于求解标量声波波动方程的方法,并推导出了时间和空间导数的高阶交错网格有限差分格式,给出了数值模拟计算... 为进一步提高标量声波波动方程正演数值模拟精度,减小数值频散和计算量,改善边界吸收效果,提出将高阶交错网格有限差分法直接用于求解标量声波波动方程的方法,并推导出了时间和空间导数的高阶交错网格有限差分格式,给出了数值模拟计算所需的稳定性条件和吸收边界条件。模型试算结果表明,与传统的中心有限差分法相比,本文方法不仅可以快速提高时空导数的差分精度,而且具有频散小、计算效率高等特点,可推广应用于三维标量声波高精度正演数值模拟和叠后逆时深度偏移。 展开更多

关键词 标量声波波动方程 高阶交错网格 有限差分格式 正演数值模拟 叠后逆时深度偏移

在线阅读 下载PDF 职称材料

利用高阶交错网格有限差分法数值模拟VTI介质井孔声场 被引量：1

2

作者 岳崇旺 王飞 《煤田地质与勘探》 CAS CSCD 北大核心 2016年第4期125-131,共7页

横向各向同性(TI)介质是岩石地球物理中常见的一种现象,研究其井孔声场传播特征对声波测井理论以及为声波测井解释提供依据具有重要意义。针对具有垂直对称轴的横向各向同性(VTI)介质,根据柱坐标系条件下的弹性波波动方程,推导了速度—... 横向各向同性(TI)介质是岩石地球物理中常见的一种现象,研究其井孔声场传播特征对声波测井理论以及为声波测井解释提供依据具有重要意义。针对具有垂直对称轴的横向各向同性(VTI)介质,根据柱坐标系条件下的弹性波波动方程,推导了速度—应力交错有限差分公式,采用时间二阶、空间十阶的交错有限差分算法对VTI介质中的井孔声场进行数值模拟。给出了在均匀介质中井孔声场不同时刻的波场快照,以及不同各向异性系数的VTI介质中的波场快照,计算了井轴上声源激发出的声波全波列波形。结果表明,在其他条件不变的条件下,VTI地层的各向异性系数的增大对横波的传播影响不大,但会使得纵波在纵向上的传播速度相对变小,径向上变化不大。各向异性系数的增大会使声波测井全波列首波信号时差变大,声波幅度略变小。 展开更多

关键词 交错网格有限差分 高阶 VTI 介质 井孔声场

在线阅读 下载PDF 职称材料

高阶交错网格有限差分法的地震波场数值模拟 被引量：3

3

作者 杨圳 李蓉艳 +1 位作者 岳继光 萧蕴诗 《北京交通大学学报》 CAS CSCD 北大核心 2012年第5期68-72,83,共6页

基于地震波波动方程和有限差分法,针对多种不同的复杂地质模型,进行高阶交错网格地震波场模拟仿真研究.首先对地震波波动方程原始形式进行推演,得到一阶速度应力弹性波方程,然后用时空二元高阶交错差分网格的方法对其进行离散化,最后对... 基于地震波波动方程和有限差分法,针对多种不同的复杂地质模型,进行高阶交错网格地震波场模拟仿真研究.首先对地震波波动方程原始形式进行推演,得到一阶速度应力弹性波方程,然后用时空二元高阶交错差分网格的方法对其进行离散化,最后对三层背斜模型、断层模型和含油模型进行了数值模拟实验,获取相应的波场快照和地震道图像进行对比分析,并讨论了差分阶次对于仿真精度的影响.仿真结果表明:算法计算精度准确,计算效率高,适合于模拟非均匀介质、复杂构造和复杂地质体的地震波场. 展开更多

关键词 地震波 波动方程 高阶 交错网格 有限差分法

在线阅读 下载PDF 职称材料

基于高阶交错网格有限差分法的海洋声场计算 被引量：4

4

作者 孟路稳 程广利 +1 位作者 张明敏 尚建华 《海军工程大学学报》 CAS 北大核心 2017年第6期55-59,共5页

为更好地认识和解释海洋声场的传播规律,基于数值计算的方法对声场的传播过程进行了模拟,给出了二维理想流体介质中的一阶声波方程,推导出相应的高阶交错网格有限差分格式,以及海水自由界面和完全匹配层吸收边界条件(PML)公式,并采用通... 为更好地认识和解释海洋声场的传播规律,基于数值计算的方法对声场的传播过程进行了模拟,给出了二维理想流体介质中的一阶声波方程,推导出相应的高阶交错网格有限差分格式,以及海水自由界面和完全匹配层吸收边界条件(PML)公式,并采用通量校正传输方法(FCT)压制数值色散。通过对海洋声场计算,验证了所采用的自由界面条件、吸收边界条件、FCT,以及声场模拟的正确性。 展开更多

关键词 声场计算 交错网格 高阶有限差分 PML FCT

在线阅读 下载PDF 职称材料

弹性波交错网格高阶有限差分法波场分离数值模拟 被引量：59

5

作者 李振春 张华 +1 位作者 刘庆敏 韩文功 《石油地球物理勘探》 EI CSCD 北大核心 2007年第5期510-515,共6页

地震波中的弹性波传播会产生纵波和横波,所以用完全弹性波波动方程进行弹性波波场数值模拟时,只能得到纵横波耦合的混合波场。本文从P波波场为无旋场,而S波波场为无散场的思路出发,推出满足此条件的一阶速度—应力弹性波波场分离方程,... 地震波中的弹性波传播会产生纵波和横波,所以用完全弹性波波动方程进行弹性波波场数值模拟时,只能得到纵横波耦合的混合波场。本文从P波波场为无旋场,而S波波场为无散场的思路出发,推出满足此条件的一阶速度—应力弹性波波场分离方程,并利用交错网格高阶有限差分法对波场分离方程进行数值模拟。模拟实例表明,此方法不但成功地将P波波场和S波波场从混合波场中分离出来,而且这种方法的稳定性好、模拟精度高,可用于弹性波传播规律研究及地震资料处理。 展开更多

关键词 弹性波 交错网格 高阶有限差分 波场分离

在线阅读 下载PDF 职称材料

任意倾斜各向异性介质中弹性波波场交错网格高阶有限差分法模拟 被引量：18

6

作者 裴正林 王尚旭 《地震学报》 CSCD 北大核心 2005年第4期441-451,共11页

给出了任意倾斜各向异性介质中二维三分量一阶应力速度弹性波方程交错网格任意偶数阶精度有限差分格式及其稳定性条件,并推导出了二维任意倾斜各向异性介质完全匹配吸收层法边界条件公式和相应的交错网格任意偶数阶精度差分格式.数值模... 给出了任意倾斜各向异性介质中二维三分量一阶应力速度弹性波方程交错网格任意偶数阶精度有限差分格式及其稳定性条件,并推导出了二维任意倾斜各向异性介质完全匹配吸收层法边界条件公式和相应的交错网格任意偶数阶精度差分格式.数值模拟结果表明,该方法模拟精度高,计算效率高,边界吸收效果好.各向异性介质中弹性波波前面形态复杂,且qP波波速不总是比qS波波速快.qS波波前面和同相轴的三分叉现象普遍,且其同相轴一般不是双曲线型.当TI介质倾斜时,3个分量上均能够观测到横波分裂现象,而且各波形的同相轴变得不对称. 展开更多

关键词 各向异性介质 交错网格 高阶有限差分 完全匹配吸收层边界条件

在线阅读 下载PDF 职称材料

非均匀介质地震波传播交错网格高阶有限差分法模拟 被引量：26

7

作者 裴正林 牟永光 《石油大学学报（自然科学版）》 CSCD 北大核心 2003年第6期17-21,共5页

采用常规的二阶声波方程有限差分方法对于非均匀介质进行了数值模拟时 ,其数值模拟精度较低。而采用一阶双曲型标量波动方程 ,则无须对介质的弹性常数进行空间求导。根据Taylor级数展开式 ,推导出了交错网格一阶空间导数的任意偶数阶精... 采用常规的二阶声波方程有限差分方法对于非均匀介质进行了数值模拟时 ,其数值模拟精度较低。而采用一阶双曲型标量波动方程 ,则无须对介质的弹性常数进行空间求导。根据Taylor级数展开式 ,推导出了交错网格一阶空间导数的任意偶数阶精度展开式和相应差分系数计算式以及一阶双曲型标量波动方程交错网格任意偶数阶精度差分格式 ,并给出了该差分算法的稳定性条件。用该差分算法对均匀介质模型、非均匀介质模型和Marmousi模型进行了数值模拟试验 ,并与伪谱法进行了对比。结果表明 ,一阶双曲型标量波动方程交错网格高阶差分法的模拟精度与伪谱法的精度非常接近 ,计算效率高 ,且适合于模拟非均匀介质。 展开更多

关键词 非均匀介质 地震波 交错网格 高阶有限差分法 数值模拟 一阶双曲型标量波动方程 油气勘探 地震勘探

在线阅读 下载PDF 职称材料

双相各向异性介质弹性波传播交错网格高阶有限差分法模拟 被引量：24

8

作者 裴正林 《石油地球物理勘探》 EI CSCD 北大核心 2006年第2期137-143,共7页

含流体孔隙各向异性介质中的弹性波传播问题是目前石油勘探和地震学研究的难点之一。基于Biot理论,本文提出了二维双相任意倾斜各向异性介质三分量弹性波方程交错网格任意偶阶精度有限差分解法,并对均匀及两层双相VTI介质和TTI介质中的... 含流体孔隙各向异性介质中的弹性波传播问题是目前石油勘探和地震学研究的难点之一。基于Biot理论,本文提出了二维双相任意倾斜各向异性介质三分量弹性波方程交错网格任意偶阶精度有限差分解法,并对均匀及两层双相VTI介质和TTI介质中的弹性波场进行了模拟。波场快照和合成记录表明:当横波源的偏振方向与TI介质的对称轴存在一定夹角、且耗散系数较小时,在固相与流相中存在快横波、慢横波、快纵波和慢纵波;双相TI介质中横波分裂现象及波前面尖角和三分叉现象与TI介质中的基本相同,只是黏滞相界中弹性波速度有所减小,且具有频散和衰减特性;从直达波转换为反射波后,该波型具有反射波的性质。 展开更多

关键词 双相各向异性介质 弹性波场 交错网格 高阶有限差分

在线阅读 下载PDF 职称材料

三维双相各向异性介质弹性波方程交错网格高阶有限差分法模拟 被引量：26

9

作者 裴正林 《中国石油大学学报（自然科学版）》 EI CAS CSCD 北大核心 2006年第2期16-20,共5页

基于Biot理论,给出了三维双相各向异性介质应力速度弹性波方程交错网格任意偶阶精度有限差分解法,并对三维双相各向异性(TIH)介质中弹性波场进行了模拟。结果表明,弹性波在三维双相各向异性介质中传播时存在快纵波qP1、快横波qS1、慢横... 基于Biot理论,给出了三维双相各向异性介质应力速度弹性波方程交错网格任意偶阶精度有限差分解法,并对三维双相各向异性(TIH)介质中弹性波场进行了模拟。结果表明,弹性波在三维双相各向异性介质中传播时存在快纵波qP1、快横波qS1、慢横波qS2和慢纵波qP2,并清楚地观测到了横波分裂、横波分裂盲点、波面三分叉等特殊现象。快纵波在固相和流相中的相位相同;而慢纵波在固相和流相中的相位相反,即慢纵波在流相中振幅大,而在固相中的振幅较小。另外,三维双相各向异性介质中弹性波场的快纵波和快横波的耦合关系、波的类型、能量分布和相位等都是在三维空间中变化的。 展开更多

关键词 三维双相各向异性介质 弹性波场 交错网格 高阶有限差分

在线阅读 下载PDF 职称材料

三维各向同性介质弹性波方程交错网格高阶有限差分法模拟 被引量：43

10

作者 裴正林 《石油物探》 EI CSCD 2005年第4期308-315,共8页

给出了三维各向同性介质中一阶应力-速度弹性波方程交错网格任意偶数阶精度有限差分格式,推导出了三维各向同性弹性介质完全匹配层吸收边界条件公式和相应的交错网格高阶有限差分格式。对三维French 模型进行了弹性波模拟,结果表明,该... 给出了三维各向同性介质中一阶应力-速度弹性波方程交错网格任意偶数阶精度有限差分格式,推导出了三维各向同性弹性介质完全匹配层吸收边界条件公式和相应的交错网格高阶有限差分格式。对三维French 模型进行了弹性波模拟,结果表明,该方法模拟精度高,边界吸收效果好。在三维French模型的xoz平面和yoz 平面的弹性波场快照中可以见到断面反射波、侧面反射波和散射波等波场特征。多次侧面反射波和多次散射波说明不同三维构造所形成的波场是相互影响的。 展开更多

关键词 三维各向同性介质 弹性波 交错网格 高阶有限差分 完全匹配层吸收边界条件

在线阅读 下载PDF 职称材料

基于递归积分的复频移PML弹性波方程交错网格高阶差分法 被引量：23

11

作者 张显文 韩立国 +2 位作者 黄玲 李翔 刘前坤 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2009年第7期1800-1807,共8页

常规完全匹配吸收边界(PML)对以近掠射角入射到界面上的波以及低频波、损耗波都会产生虚假边界反射.基于递归积分的不分裂复频移PML算法,利用复频移拉伸函数,极大地改善了 PML边界条件的性能,我们进一步推导出基于递归积分的不分裂复频... 常规完全匹配吸收边界(PML)对以近掠射角入射到界面上的波以及低频波、损耗波都会产生虚假边界反射.基于递归积分的不分裂复频移PML算法,利用复频移拉伸函数,极大地改善了 PML边界条件的性能,我们进一步推导出基于递归积分的不分裂复频移PML弹性波方程交错网格高阶差分法,对长条形介质模型进行数值模拟,与常规PML算法进行比较说明该算法对以掠射角入射到PML界面的波以及PML层内损耗波的吸收效果. 展开更多

关键词 交错网格 高阶有限差分 完全匹配层 递归积分

在线阅读 下载PDF 职称材料

三维感应测井响应计算的交错网格有限差分法 被引量：48

12

作者 汪功礼 张庚骥 +1 位作者 崔锋修 高峰 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2003年第4期561-567,共7页

应用交错网格有限差分法计算三维复杂环境中的感应测井响应 .其中 ,利用Krylov子空间不变性求解离散得到的大型稀疏复对称线性方程组 .在构造Krylov子空间时使用其系数矩阵的伪逆以改善迭代的收敛性 .迭代中 ,使用不完全Cholesky分解共... 应用交错网格有限差分法计算三维复杂环境中的感应测井响应 .其中 ,利用Krylov子空间不变性求解离散得到的大型稀疏复对称线性方程组 .在构造Krylov子空间时使用其系数矩阵的伪逆以改善迭代的收敛性 .迭代中 ,使用不完全Cholesky分解共轭梯度法求解 4个三维Poisson方程以得到新的Lanczos向量 .通常迭代不超过 2 0次可得到理想结果 .另外 ,提出一种新的物质平均公式以计算电导率平均值 ,可保证电流守恒 . 展开更多

关键词 三维感应测井响应 交错网格有限差分法 计算 KRYLOV子空间 电导率平均值 地球物理勘探

在线阅读 下载PDF 职称材料

交错网格有限差分法模拟微动信号 被引量：6

13

作者 周晓华 陈祖斌 +1 位作者 曾晓献 焦健 《吉林大学学报（地球科学版）》 EI CAS CSCD 北大核心 2012年第3期852-857,共6页

为了较好地模拟给定模型在随机源激发下产生的稳定波场信息,为微动信号的检验和数据特征分析等提供理论依据,采用交错网格有限差分方法模拟微动信号,并根据瑞雷波和勒夫波的产生机理,通过测区内震源的分布情况估测微动信号中面波的组成... 为了较好地模拟给定模型在随机源激发下产生的稳定波场信息,为微动信号的检验和数据特征分析等提供理论依据,采用交错网格有限差分方法模拟微动信号,并根据瑞雷波和勒夫波的产生机理,通过测区内震源的分布情况估测微动信号中面波的组成成分。首先建立交错网格有限差分法模拟微动信号的数值模型,然后给出估测微动信号中面波组成成分的算式,最后采用典型方法数值模拟微动信号,应用空间自相关方法提取模拟信号中面波的频散曲线进行分析对比。结果表明,交错网格有限差分方法模拟得到的微动信号更符合实测结果。 展开更多

关键词 交错网格 有限差分法 微动 面波 频散曲线

在线阅读 下载PDF 职称材料

多分量感应测井响应的交错网格有限差分法模拟 被引量：8

14

作者 王昌学 杨韡 +3 位作者 储昭坦 陶果 沈金松 朱益华 《石油大学学报（自然科学版）》 EI CAS CSCD 北大核心 2005年第3期35-40,共6页

从电磁响应方程出发,应用交错网格有限差分法推导了三维各向异性地层多分量感应测井响应的差分计算格式,并利用所形成的算法考察了多种环境条件下多分量感应测井的响应规律。计算结果表明,多分量感应测井各磁场分量不仅与地层水平电阻... 从电磁响应方程出发,应用交错网格有限差分法推导了三维各向异性地层多分量感应测井响应的差分计算格式,并利用所形成的算法考察了多种环境条件下多分量感应测井的响应规律。计算结果表明,多分量感应测井各磁场分量不仅与地层水平电阻率和垂直电阻率有关,还与仪器方位角、井斜角、井眼轴与地层走向的相对位置有关。根据仪器方位角与井斜角间的关系和各测量量,可反演出地层水平电阻率和垂直电阻率,进而改善电性各向异性储层的含油饱和度评价。计算结果表明该算法是准确可靠的。 展开更多

关键词 多分量感应测井响应 电磁响应方程 交错网格有限差分法 电阻率 各向异性储层 含油饱和度 电测井

在线阅读 下载PDF 职称材料

基于横向各向同性BISQ模型的弹性波高阶交错网格有限差分数值模拟 被引量：9

15

作者 李红星 刘财 陶春辉 《石油地球物理勘探》 EI CSCD 北大核心 2007年第6期686-693,共8页

孔隙介质模型主要有Biot模型与BISQ模型。BISQ模型同时考虑了Biot流动与喷射流动,能够获得比Biot理论更为准确的结果。弹性波数值模拟有利于认识波的传播特性与过程,而采用高阶交错网格有限差分法数值模拟的空间频散性比低阶差分法小。... 孔隙介质模型主要有Biot模型与BISQ模型。BISQ模型同时考虑了Biot流动与喷射流动,能够获得比Biot理论更为准确的结果。弹性波数值模拟有利于认识波的传播特性与过程,而采用高阶交错网格有限差分法数值模拟的空间频散性比低阶差分法小。本文从BISQ模型弹性波的本构方程和运动方程出发,推导出了基于横向各向同性BISQ模型的弹性波高阶交错网格有限差分算法,还进行了数值模拟。通过调整差分阶数、频率以及黏滞系数,研究了兼顾两种流动机制时弹性波的传播特性,并分析了BISQ模型下弹性波的透射与反射。 展开更多

关键词 数值模拟 高阶有限差分 BISQ模型 Biot模型 交错网格 孔隙介质 横向各向同性

在线阅读 下载PDF 职称材料

高阶有限差分法管波传播数值模拟 被引量：5

16

作者 孟庆生 樊玉清 +1 位作者 张珂 张盟 《吉林大学学报（地球科学版）》 EI CAS CSCD 北大核心 2011年第1期292-298,共7页

为正确认识流体充填钻孔中管波的传播规律,首先给出均匀各向同性介质中包含流体充填钻孔条件下的弹性波动力学方程和二阶速度-应力公式及差分格式的稳定性条件和边界条件;然后利用高阶交错网格有限差分法,结合典型模型,模拟了不同时刻... 为正确认识流体充填钻孔中管波的传播规律,首先给出均匀各向同性介质中包含流体充填钻孔条件下的弹性波动力学方程和二阶速度-应力公式及差分格式的稳定性条件和边界条件;然后利用高阶交错网格有限差分法,结合典型模型,模拟了不同时刻管波的传播快照和复杂层状介质中管波的合成地震记录,并与实际管波地震记录进行了对比。结果表明:在地层分界面、溶洞的边界等波阻抗差异界面处,管波都会产生反射;反射管波的振幅、频率特征受界面两侧的岩石物理性质差异的影响;在理论合成记录中,反射管波与直达管波同相轴的交点与理论模型中地层分界面或异常体的边界一致。 展开更多

关键词 勘查地球物理 交错网格 有限差分法 管波 数值模拟

在线阅读 下载PDF 职称材料

高阶交错网格有限差分弹性波场模拟的精度分析 被引量：2

17

作者 岳晓鹏 白超英 岳崇旺 《煤田地质与勘探》 CAS CSCD 北大核心 2017年第1期125-130,共6页

交错网格波场数值模拟是目前地震正演中广泛使用的方法,为对比分析不同阶数的差分格式下产生的计算效率和精度差异,重新推导了弹性波方程的4种时间4阶、空间2N阶的差分公式及系数,并计算了他们的稳定性条件。利用这4种差分格式进行弹性... 交错网格波场数值模拟是目前地震正演中广泛使用的方法,为对比分析不同阶数的差分格式下产生的计算效率和精度差异,重新推导了弹性波方程的4种时间4阶、空间2N阶的差分公式及系数,并计算了他们的稳定性条件。利用这4种差分格式进行弹性波场数值模拟,对比分析了波场快照、合成地震记录及CPU时间。结果表明:时间4阶、空间6+6阶精度的交错网格有限差分方法在进行地震波场数值模拟时具有较高的计算精度和计算效率。 展开更多

关键词 精度分析 弹性波场 交错网格 高阶有限差分

在线阅读 下载PDF 职称材料

非均匀TI介质P-SV波传播交错网格高阶有限差分数值模拟 被引量：4

18

作者 黄翼坚 《西北地震学报》 CSCD 北大核心 2008年第1期11-16,共6页

给出了在非均匀横向各向同性(TI)介质情况下,四阶时间精度、高阶空间精度的一阶速度-应力P-SV波的波动方程交错网格有限差分解法。首先根据一阶速度(应力)波动方程把速度(应力)对时间的一阶和三阶导数转换为应力(速度)对空间的导数,从... 给出了在非均匀横向各向同性(TI)介质情况下,四阶时间精度、高阶空间精度的一阶速度-应力P-SV波的波动方程交错网格有限差分解法。首先根据一阶速度(应力)波动方程把速度(应力)对时间的一阶和三阶导数转换为应力(速度)对空间的导数,从而在使用四阶时间精度有限差分格式计算某一时刻的波场时只需要前面两个时间步的波场值;然后在空间上采用高阶有限差分格式以提高数值模拟的精度。数值模拟结果和实测垂直地震剖面(VSP)记录符合得很好,说明该方法是可行的。 展开更多

关键词 非均匀TI介质 P-SV波 交错网格 高阶有限差分 数值模拟

在线阅读 下载PDF 职称材料

一阶弹性波方程的变网格高阶有限差分数值模拟 被引量：30

19

作者 李振春 张慧 张华 《石油地球物理勘探》 EI CSCD 北大核心 2008年第6期711-716,611+742,共8页

使用可变网格的有限差分法进行地震模拟有许多独特的优点,主要表现为对地质模型的离散化更为合理,在低速带和复杂构造区域,可将局部网格划分得相对精细些,不仅提高了模拟精度,消除了因采样不足导致的频散现象,而且可以减少计算机内存需... 使用可变网格的有限差分法进行地震模拟有许多独特的优点,主要表现为对地质模型的离散化更为合理,在低速带和复杂构造区域,可将局部网格划分得相对精细些,不仅提高了模拟精度,消除了因采样不足导致的频散现象,而且可以减少计算机内存需求,保持模型计算的灵活性。本文提出一种新的基于高阶交错网格技术的弹性波数值模拟方法,通过改变网格的空间步长实现了局部网格加密技术,弥补了常规网格的缺陷和不足。试算结果表明,本文提出的算法稳定性较好,且能够提高模拟精度,减少计算时间,提高计算效率。 展开更多

关键词 变网格 交错网格 高阶有限差分 弹性波

在线阅读 下载PDF 职称材料

求解二阶解耦弹性波方程的低秩分解法和低秩有限差分法 被引量：8

20

作者 袁雨欣 胡婷 +2 位作者 王之洋 郭鹏 刘洪 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2018年第8期3324-3333,共10页

时间域的波场延拓方法在本质上都可以归结为对一个空间-波数域算子的近似.本文基于一阶波数-空间混合域象征,提出一种新的方法求解解耦的二阶位移弹性波方程.该方法采用交错网格,连续使用两次一阶前向和后向拟微分算子,推导得到了解耦... 时间域的波场延拓方法在本质上都可以归结为对一个空间-波数域算子的近似.本文基于一阶波数-空间混合域象征,提出一种新的方法求解解耦的二阶位移弹性波方程.该方法采用交错网格,连续使用两次一阶前向和后向拟微分算子,推导得到了解耦的二阶位移弹性波方程的波场延拓算子.由于该混合域象征在伪谱算子的基础上增加了一个依赖于速度模型的补偿项,可以补偿由于采用二阶中心差分计算时间微分项带来的误差,有效地减少模拟结果的数值频散,提高模拟精度.然而,在非均匀介质中,直接计算该二阶的波场延拓算子,每一个时间步上需要做N次快速傅里叶逆变换,其中N是总的网格点数.为了减少计算量,提出了交错网格低秩分解方法;针对常规有限差分数值频散问题,本文将交错网格低秩方法与有限差分法结合,提出了交错网格低秩有限差分法.数值结果表明,交错网格低秩方法和交错网格低秩有限差分法具有较高的精度,对于复杂介质的地震波数值模拟和偏移成像具有重要的价值. 展开更多

关键词 波场延拓 波数-空间混合域象征 交错网格低秩分解方法 交错网格低秩有限差分法

在线阅读 下载PDF 职称材料