复杂固-液介质的速度应力方程三角形单元间断伽辽金地震波模拟算法

1

作者 曹文忠 张伟 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2024年第2期684-695,共12页

间断伽辽金法可用于存在复杂边界条件的模型,同时具有高阶精度和易于并行计算的优点,因此近年来在地震波传播模拟研究中得到了快速发展.数值通量是间断伽辽金法的关键组成部分之一.相比于其他通量,基于Rankine-Hugoniot跳跃条件的通量(R... 间断伽辽金法可用于存在复杂边界条件的模型,同时具有高阶精度和易于并行计算的优点,因此近年来在地震波传播模拟研究中得到了快速发展.数值通量是间断伽辽金法的关键组成部分之一.相比于其他通量,基于Rankine-Hugoniot跳跃条件的通量(RH-condition通量)在固体和液体介质边界具有更宽的稳定性,可以使用更大的时间步长,尤其是液体和固体之间的波阻抗差异较大时.它已被用于基于四边形网格和速度-应变方程的间断伽辽金地震波固液介质模拟.本文为了模拟陆地自然水体等存在复杂固-液界面形状的地震波传播,并在未来与其他数值方法耦合,发展了基于三角形网格的RH-condition通量间断伽辽金方法,使用一阶速度-压力声波和一阶速度-应力弹性波方程模拟复杂固-液介质中地震波的传播.通过水平层状和sin型起伏固-液模型,验证了该方法模拟结果的准确性.通过数值模拟展示了所提出的间断伽辽金法在不同网格大小和阶数下的准确性和效率.最后通过一个复杂模型的例子表明存在复杂固-液界面时该方法可以准确施加固-液边界条件. 展开更多

关键词 间断伽辽金 地震波模拟 固-液介质 数值通量

在线阅读 下载PDF 职称材料

三维复杂构造中地震波模拟的单程波方法 被引量：25

2

作者 刘礼农 崔凤林 张剑锋 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2004年第3期514-520,共7页

复杂构造中单程波与双程波方法模拟结果的比较表明 ,就地震勘探中主要关心的一次反射波而言 ,单程波算法已具有足够的精度 .使用单程波方程将极大地减少数值计算的计算量 ,同时对介质的几何和物理参数建模也降低了要求 .单程波算法可视... 复杂构造中单程波与双程波方法模拟结果的比较表明 ,就地震勘探中主要关心的一次反射波而言 ,单程波算法已具有足够的精度 .使用单程波方程将极大地减少数值计算的计算量 ,同时对介质的几何和物理参数建模也降低了要求 .单程波算法可视为深度偏移的“逆运算” ,这样可以很好地借用已知的深度偏移方法及其程序系统 .基于计算效率和计算精度的双重考虑 ,本文在介质速度结构较复杂时采用显式短算子波场延拓方法 ,而在介质速度结构相对简单时采用分裂步相移法 . 展开更多

关键词 地震波模拟 单程波方法 三维地震模拟 显式短算子 分裂步相移

在线阅读 下载PDF 职称材料

界面起伏条件下反射/透射算子+单程波方程的地震波模拟方法 被引量：15

3

作者 谢桂生 刘洪 +1 位作者 李幼铭 胡润苗 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2005年第5期1172-1178,共7页

本文导出了一种由单程波方程利用反射/透射算子的可分表示方法模拟复杂介质中一次反射地震波的数值算法.文中利用算子可分表示理论将反射/透射算子分解成适合于双域(空间域和波数域)运算的表达形式,使得本文得到的地震波数值模拟算法可... 本文导出了一种由单程波方程利用反射/透射算子的可分表示方法模拟复杂介质中一次反射地震波的数值算法.文中利用算子可分表示理论将反射/透射算子分解成适合于双域(空间域和波数域)运算的表达形式,使得本文得到的地震波数值模拟算法可适应于一定程度横向非均匀介质和界面起伏情况,在入射角小于45°时能够准确模拟振幅随入射角(AVA)的关系.就模拟一次反射地震波而言,与前人研究的双程波动方程伪谱法地震波模拟相比,本文算法具有足够高的模拟精度,且计算效率成倍地提高. 展开更多

关键词 地震波模拟 单程波方程 反射/透射算子 可分表示

在线阅读 下载PDF 职称材料

细胞自动机地震波模拟的并行化算法 被引量：15

4

作者 王真理 李幼铭 《地球物理学报》 SCIE EI CAS CSCD 北大核心 1999年第3期410-415,共6页

采用区域分解的办法对格子气自动机地震波传播模拟的并行化算法进行了研究，提出了两种区域分解算法，并在SGIPowerChallenge并行计算机上进行了实际地质模型中平面波反射地震响应的计算，获得了较高的并行计算效率．

关键词 并行计算 地震波传播 地震波模拟 细胞自动化

在线阅读 下载PDF 职称材料

地震波模拟中多轴复频移近似完全匹配层的吸收效果及稳定性 被引量：4

5

作者 罗玉钦 刘财 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2020年第8期3078-3090,共13页

在进行地震波模拟计算的过程中用有限的计算区域模拟地下无限空间,需要进行边界截断.为了在边界处不产生虚假反射影响模拟结果,需要引入吸收边界条件.本文采用的近似完全匹配层是一种新型非分裂完全匹配层,计算效率较高.同时相比于其他... 在进行地震波模拟计算的过程中用有限的计算区域模拟地下无限空间,需要进行边界截断.为了在边界处不产生虚假反射影响模拟结果,需要引入吸收边界条件.本文采用的近似完全匹配层是一种新型非分裂完全匹配层,计算效率较高.同时相比于其他非分裂完全匹配层,其还具有不改变方程的形式、易于实现等优势.但是当入射波角度较大,边界吸收效果变弱,且残留在边界中的能量使近似完全匹配层变得极其不稳定.多轴复频移近似完全匹配层的提出就是为了改善对大角度入射波的吸收并且提高边界的稳定性.通过实验模拟和矩阵特征值灵敏度来研究多轴复频移近似完全匹配层的吸收效果及稳定性.结果表明该方法不仅能够吸收掠入波,而且对常规入射波的吸收也得到提升,同时拥有更好的稳定性. 展开更多

关键词 完全匹配层 地震波模拟 多轴 复频移

在线阅读 下载PDF 职称材料

反转构造模型的地震波模拟及成像 被引量：1

6

作者 董良国 周祖翼 《高校地质学报》 CAS CSCD 2001年第2期189-195,共7页

为了检验由地震剖面解释得出的反转构造的可靠性 ,探讨反转构造地区在地震数据采集、处理中应采用的有效方法 ,设计了一个典型的反转构造模型———突起构造。以该模型为基础 ,针对复杂构造地震波模拟中常见的数值频散严重、人为边界反... 为了检验由地震剖面解释得出的反转构造的可靠性 ,探讨反转构造地区在地震数据采集、处理中应采用的有效方法 ,设计了一个典型的反转构造模型———突起构造。以该模型为基础 ,针对复杂构造地震波模拟中常见的数值频散严重、人为边界反射吸收不理想两个问题 ,应用高阶差分法以及一种新的边界吸收条件 ,对反转模型进行了地震波场数值模拟和深度偏移成像处理。采用中心放炮两边接收的观测方式 ,共模拟了 185炮地震记录。研究表明 ,波动方程高阶差分法是模拟复杂构造模型中地震波传播的有效方法 ,除了近于直立的断层因难以接收到其反射波而较难成像外 ,其余的反转断层及其组合在理论上是可以通过合适方法 (如Fourier有限差分深度偏移 )较好地成像的。 展开更多

关键词 反转构造 地震波场模拟 深度偏移成像 突起构造 地震剖面

在线阅读 下载PDF 职称材料

有限元地震波模拟中的点状震源 被引量：2

7

作者 范书蕊 《石油地球物理勘探》 EI CSCD 北大核心 1992年第2期223-229,共7页

用有限元模拟地震波时,常采用点状震源。震源模拟的质量直接影响到数值模拟的结果。本文采用在震源点周围附加几个虚设震源的办法,使震源附近介质受力均匀,取得了满意的结果。

关键词 有限元法 地震波模拟 点状震源

在线阅读 下载PDF 职称材料

基于多辛结构谱元法的保结构地震波场模拟 被引量：19

8

作者 汪文帅 李小凡 +1 位作者 鲁明文 张美根 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2012年第10期3427-3439,共13页

近年来构造高精度、高效且具有长时程跟踪能力的保结构算法已逐渐成为地震波模拟算法发展的重要方向之一.本文基于谱元法(SEM)进行空间域离散结合新推导的三阶辛算法(NTSTO)进行时间域离散,构造了一种具有时-空保结构特性的新算法.本文... 近年来构造高精度、高效且具有长时程跟踪能力的保结构算法已逐渐成为地震波模拟算法发展的重要方向之一.本文基于谱元法(SEM)进行空间域离散结合新推导的三阶辛算法(NTSTO)进行时间域离散,构造了一种具有时-空保结构特性的新算法.本文给出的多组数值试验对比结果表明,本算法无论在内存消耗、稳定性及计算耗时,还是长时程跟踪能力方面都有上佳的表现;另外,本文给出的起伏地表多层介质模型的数值算例验证了该算法处理复杂几何形状和复杂介质时的有效性.该多辛结构谱元法的发展将为长时程地震波传播的计算及模拟提供更为广泛而有效的选择. 展开更多

关键词 谱元法 辛算法 地震波模拟

在线阅读 下载PDF 职称材料

基于Shannon奇异核理论的褶积微分算子在地震波场模拟中的应用 被引量：12

9

作者 龙桂华 李小凡 张美根 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2009年第4期1014-1024,共11页

本文在前人工作的基础上,建立了一种基于shannon奇异核的交错网格褶积微分算子方法.文中不仅详细讨论了影响算子精度的各种因素,同时也着重分析了其在弹性波模拟中的频散关系和稳定性条件.通过和交错网格有限差分算子比较,发现该算子即... 本文在前人工作的基础上,建立了一种基于shannon奇异核的交错网格褶积微分算子方法.文中不仅详细讨论了影响算子精度的各种因素,同时也着重分析了其在弹性波模拟中的频散关系和稳定性条件.通过和交错网格有限差分算子比较,发现该算子即使在高波数域也具有较高的精度.均匀介质中的数值试验也表明,该方法9点格式就基本上达到了解析解精度.而分层均匀介质和复杂介质中的地震波数值模拟也同时证实了该方法精度高,稳定性好,是一种研究复杂介质中地震波传播的有效数值方法. 展开更多

关键词 离散奇异核 褶积微分算子 地震波模拟

在线阅读 下载PDF 职称材料

粘弹介质中可变网格地震波传播数值模拟 被引量：7

10

作者 李晓波 董良国 《石油物探》 EI CSCD 北大核心 2012年第1期1-10,111,共10页

从Kelvin介质的本构方程出发推导了粘弹介质中的弹性波方程和声波方程,建立了粘滞弹性波方程的交错网格高阶差分解和粘滞声波方程的高阶差分解,并应用可变空间网格与局部可变时间步长数值计算技术,形成了适用于粘弹介质中弹性(声)波传... 从Kelvin介质的本构方程出发推导了粘弹介质中的弹性波方程和声波方程,建立了粘滞弹性波方程的交错网格高阶差分解和粘滞声波方程的高阶差分解,并应用可变空间网格与局部可变时间步长数值计算技术,形成了适用于粘弹介质中弹性(声)波传播的高效、高精度的数值模拟方法。数值实验结果表明,该方法可以有效模拟地震波在粘弹介质中的传播,并且可以比较准确地刻画精细介质结构,在保证模拟精度的同时极大地提高了计算效率。 展开更多

关键词 粘弹性 Kelvin介质 可变空间网格 局部可变时间步长 地震波模拟 高阶差分 交错网格

在线阅读 下载PDF 职称材料

地震波数值模拟中的最优Shannon奇异核褶积微分算子 被引量：3

11

作者 龙桂华 赵宇波 赵家福 《地震学报》 CSCD 北大核心 2011年第5期650-662,共13页

如何有效协调地震波数值模拟过程中的计算精度和计算效率,一直是许多地球物理工作者密切关注且努力寻求解决的问题.这一问题求解的好坏,将直接影响到人们正确认识地下构造的地球物理响应.本文从离散Shannon奇异核理论出发,推导了基于Sha... 如何有效协调地震波数值模拟过程中的计算精度和计算效率,一直是许多地球物理工作者密切关注且努力寻求解决的问题.这一问题求解的好坏,将直接影响到人们正确认识地下构造的地球物理响应.本文从离散Shannon奇异核理论出发,推导了基于Shannon奇异核的交错网格褶积微分算子,分析并求出了影响算子精度窗函数的最优化参数.通过将最优褶积微分算子与不同长度的交错网格有限差分算子进行比较,发现该算子8点格式的精度相当于16阶交错网格有限差分.对于同长度的算子,交错网格褶积微分算子可以在一个波长内具有较少的采样点数.将本文算子运用于Corner-Edge固液介质模型,其数值试验结果同样表明,该方法计算精度高,稳定性好,可作为研究复杂介质中地震波传播的又一有效数值方法. 展开更多

关键词 交错网格 Shannon奇异核 褶积微分算子 地震波模拟

在线阅读 下载PDF 职称材料

球坐标系下谱元法三维地震波场模拟 被引量：7

12

作者 董兴朋 杨顶辉 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2017年第12期4671-4680,共10页

谱元法已成为区域性乃至大陆性尺度地震波场模拟的重要工具.对于区域或大陆尺度层析成像而言,地球曲率不可忽略,此时模拟地震波传播采用球坐标系更为合适.本文从球坐标系下弹性波动方程弱形式出发,基于球坐标系变分原理给出了球坐标系... 谱元法已成为区域性乃至大陆性尺度地震波场模拟的重要工具.对于区域或大陆尺度层析成像而言,地球曲率不可忽略,此时模拟地震波传播采用球坐标系更为合适.本文从球坐标系下弹性波动方程弱形式出发,基于球坐标系变分原理给出了球坐标系下求解三维地震波方程的谱元法.另一方面,计算Fréchet敏感核是进行全波形反演的关键,本文借助伴随原理,推导了全波走时层析成像三维Fréchet敏感核表达式.为了验证球坐标系下谱元法的精度,我们将数值模拟结果与normal mode方法得到的解析解在1-D PREM模型下进行了对比.同时,我们将此方法应用到华北克拉通区域,以期获得地球内部结构精确成像.基于3-D全球径向各向异性地幔模型S362ANI和3-D地壳模型Crust1.0,我们建立了华北克拉通初始3-D背景模型,并将数值模拟结果与实际观测台站记录波形资料进行对比分析,利用互相关方法提取走时残差,最后给出了Fréchet敏感核在3-D空间中的分布,这些工作为下一步开展球坐标系下三维大尺度全波形反演奠定了基础. 展开更多

关键词 球坐标系 谱元法 Fréchet敏感核 地震波模拟

在线阅读 下载PDF 职称材料

梯形网格伪谱法地震波场模拟 被引量：6

13

作者 谭文卓 吴帮玉 +1 位作者 李博 雷军 《石油地球物理勘探》 EI CSCD 北大核心 2020年第6期1282-1291,1163-1164,共12页

波动方程的数值求解是地震勘探高精度成像和反演方法的核心。为保证精度,常规方法往往根据模型最小速度设置固定的空间采样间隔,易在高速区形成"过采样",计算存在一定冗余性。由于重力引起的岩石压实效应,波的传播速度由浅入... 波动方程的数值求解是地震勘探高精度成像和反演方法的核心。为保证精度,常规方法往往根据模型最小速度设置固定的空间采样间隔,易在高速区形成"过采样",计算存在一定冗余性。由于重力引起的岩石压实效应,波的传播速度由浅入深整体增大,梯形坐标变换更契合这种速度变化趋势。为此,提出梯形网格伪谱法地震波模拟方法。利用基于梯形坐标变换的梯形网格法对介质进行剖分,在浅层低速区采用细网格、深部高速区采用粗网格,可有效减少网格点数;同时,通过伪谱法求解坐标变换后的变系数声波方程,并引入完全匹配层消除人工边界虚假反射,以兼顾波场模拟的效率与精度。Marmousi模型的测试结果表明,相较于常规固定网格剖分,梯形网格剖分的网格数可减少约69%;相较于常规网格伪谱法和高阶有限差分法,梯形网格伪谱法的计算时间分别减少约58%和60%,同时具有较小的数值频散。可见,所提方法是一种高效、高精度地震波模拟方法。 展开更多

关键词 梯形坐标变换 地震波模拟 伪谱法 声波方程

在线阅读 下载PDF 职称材料

三维粘弹介质地震波场有限差分并行模拟 被引量：25

14

作者 王德利 雍运动 +1 位作者 韩立国 廉玉广 《西北地震学报》 CSCD 北大核心 2007年第1期30-34,共5页

有限差分法在三维粘弹性复杂介质正演模拟地震波的传播中对计算机内存和计算速度要求比较高,单个PC机或工作站只能计算较少网格内短时间的波场。本文介绍一种基于MPI的并行有限差分法,可在PCCluster上模拟较大规模三维粘弹性复杂介质中... 有限差分法在三维粘弹性复杂介质正演模拟地震波的传播中对计算机内存和计算速度要求比较高,单个PC机或工作站只能计算较少网格内短时间的波场。本文介绍一种基于MPI的并行有限差分法,可在PCCluster上模拟较大规模三维粘弹性复杂介质中地震波传播时的波场;可预测地震波在此类条件下传播时的运动学和动力学性质。对于更好地理解波动传播现象,解释实际地震资料及反问题的解决等均具有重要的理论与实际意义。 展开更多

关键词 三维粘弹性介质 地震波场模拟 并行计算 MPI 有限差分

在线阅读 下载PDF 职称材料

模拟地震波场的伪谱和高阶有限差分混合方法 被引量：17

15

作者 魏星 王彦宾 陈晓非 《地震学报》 CSCD 北大核心 2010年第4期392-400,共9页

伪谱法是一种高效、高精度计算非均匀介质地震波传播的数值方法,但是由于它的微分算子的全局性,使得该方法不适用于分散内存的并行计算.本文将有限差分算子的局部性和伪谱法算子的高效、高精度相结合,发展基于两种方法的伪谱/有限差分... 伪谱法是一种高效、高精度计算非均匀介质地震波传播的数值方法,但是由于它的微分算子的全局性,使得该方法不适用于分散内存的并行计算.本文将有限差分算子的局部性和伪谱法算子的高效、高精度相结合,发展基于两种方法的伪谱/有限差分混合方法.该方法在一个空间坐标方向上利用交错网格高阶有限差分算子,在另外的空间坐标方向上利用交错网格伪谱法算子,既保留了后者的高效、高精度优势,又便于在PC集群上实现并行计算.对二维模型的计算显示,混合方法能有效处理介质不连续面,在保证伪谱法计算精度的情况下,提供了一种并行计算的可能途径. 展开更多

关键词 交错网格 伪谱法 高阶有限差分法 混合方法 地震波场模拟

在线阅读 下载PDF 职称材料

地震波数值模拟的非规则网格PML吸收边界 被引量：32

16

作者 徐义 张剑锋 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2008年第5期1520-1526,共7页

以格子法为基础,以声波方程为例研究非规则网格PML(Perfectly Matched Layer)方法.本方法的核心是建立局部坐标系下的分裂方程和基于积分近似的微分方程弱形式.该非规则网格模拟方法允许在计算域内设置任意形状的人工边界.对于二维半空... 以格子法为基础,以声波方程为例研究非规则网格PML(Perfectly Matched Layer)方法.本方法的核心是建立局部坐标系下的分裂方程和基于积分近似的微分方程弱形式.该非规则网格模拟方法允许在计算域内设置任意形状的人工边界.对于二维半空间问题,与采用矩形人工边界相比,采用半圆形人工边界可减少计算量20%以上.采用光滑的曲边界,不仅可减少计算区域,还可避免常规的PML吸收边界在吸收带角点区域的特殊处理.本方法事先计算和存储边界单元的局部几何参数,在计算的每一时间步查表调用这些参数,与常规的直边界PML方法相比,不增加任何计算量. 展开更多

关键词 格子法 非规则网格 PML 吸收边界 声波方程 地震波场模拟

在线阅读 下载PDF 职称材料

两组任意夹角倾斜裂缝介质地震波场模拟 被引量：6

17

作者 王月英 孙赞东 白海军 《石油地球物理勘探》 EI CSCD 北大核心 2010年第1期47-54,共8页

两组任意夹角的倾斜裂缝在自然界中具有普遍性,本文针对两组任意夹角的倾斜裂缝介质模型,采用交错网格高阶差分法进行地震波场二维三分量数值模拟。根据Hudson理论和Bond矩阵变换得到裂缝型储层在观测坐标下的弹性系数,对比了在本构坐... 两组任意夹角的倾斜裂缝在自然界中具有普遍性,本文针对两组任意夹角的倾斜裂缝介质模型,采用交错网格高阶差分法进行地震波场二维三分量数值模拟。根据Hudson理论和Bond矩阵变换得到裂缝型储层在观测坐标下的弹性系数,对比了在本构坐标和观测坐标下地震波场之间的差异,分析了不同方位角下炮记录波场特征、不同炮检距方位角道集上的波场特征,以及改变裂缝夹角后在方位角道集记录上的响应特征。研究表明,裂缝型储层的各向异性特征是裂缝的综合反映,各向异性的强弱不仅与裂缝的发育强度有关,还与裂缝发育的复杂程度有关。通过对方位各向异性介质数值模拟和波场特征的分析,可为方位各向异性介质的裂缝反演提供理论基础。 展开更多

关键词 地震波场模拟 裂缝介质 夹角 方位各向异性介质 裂缝型储层 波场特征 观测坐标 数值模拟

在线阅读 下载PDF 职称材料

用于地震波场模拟的变网格边长声格固体模型 被引量：8

18

作者 刘劲松 许云 乌达巴拉 《地球物理学报》 SCIE EI CAS CSCD 北大核心 1999年第4期536-542,共7页

提出一种新的用于地震波场模拟的变网格边长声格固体模型（Phononiclatticesolidwithvariousgridlength，简称PLSVL模型）．该模型通过改变不同介质中网格的边长来体现介质速度的变化，粒子在所有网格中的运动均保持每一步运行一个网格... 提出一种新的用于地震波场模拟的变网格边长声格固体模型（Phononiclatticesolidwithvariousgridlength，简称PLSVL模型）．该模型通过改变不同介质中网格的边长来体现介质速度的变化，粒子在所有网格中的运动均保持每一步运行一个网格边长．详细推导了该模型的Boltzmann方程，证明该方程及其散射项的表达式与PLS模型相同，因而从该方程出发，可寻出宏观变量所满足的波动方程．新模型在同种网格内部的传输过程不存在误差，在界面产生的误差不会随时间而积累．本文同时给出了一个理论模型的波场模拟结果． 展开更多

关键词 地震波场模拟 声格固体模型 变网格边长

在线阅读 下载PDF 职称材料

声格固体模型地震波场模拟 被引量：6

19

作者 刘劲松 许云 乌达巴拉 《石油地球物理勘探》 EI CSCD 北大核心 1997年第3期370-375,共6页

PeterMora为模拟复杂介质中的纵波，提出了声格固体模型，即在格子Boltzman模型基础上引入变速粒子而建立的细胞自动机模型，从而为利用这种方法模拟地震波场奠定了理论基础。文中简要介绍了PeterMora的声格固体模型的基本理论框架及差... PeterMora为模拟复杂介质中的纵波，提出了声格固体模型，即在格子Boltzman模型基础上引入变速粒子而建立的细胞自动机模型，从而为利用这种方法模拟地震波场奠定了理论基础。文中简要介绍了PeterMora的声格固体模型的基本理论框架及差分算法，同时还介绍了黄联捷的插值算法。笔者的工作在于提出吸收边界的处理方法，并给出一个非均匀介质模型的数值模拟结果。由于此法允许在每个网格点上的状态可同时更新，因而适合在大型并行机上实现。对于波阻抗和形状变化剧烈的介质，不必作渐变和平滑处理，表明该法是有潜力的。 展开更多

关键词 声格固体模型 地震波场模拟 吸收边界 地震勘探

在线阅读 下载PDF 职称材料

标量地震波频率-空间域有限差分法数值模拟 被引量：17

20

作者 任浩然 王华忠 龚婷 《石油物探》 EI CSCD 北大核心 2009年第1期20-26,共7页

讨论了频率-空间域优化有限差分算子地震波场模拟方法。推导了频率-空间域双程波高阶有限差分算子,计算了9点有限差分系数;提出了一种节省内存的存储策略,对稀疏矩阵和稀疏矩阵LU分解后的矩阵采用了节省内存的存储策略,使内存空间的需... 讨论了频率-空间域优化有限差分算子地震波场模拟方法。推导了频率-空间域双程波高阶有限差分算子,计算了9点有限差分系数;提出了一种节省内存的存储策略,对稀疏矩阵和稀疏矩阵LU分解后的矩阵采用了节省内存的存储策略,使内存空间的需求量大幅度降低,可以利用本方法进行具有生产规模的地震波正演模拟;推导了依赖于频率的最佳匹配层吸收边界条件,获得了很好的吸收效果。洼陷模型和Marmousi模型数值模拟试验表明,方法可行且有效,在节约计算机资源的同时,计算精度也能得到保证。 展开更多

关键词 频率-空间域 有限差分法 标量波方程 最佳匹配层吸收边界条件 地震波正演模拟

在线阅读 下载PDF 职称材料