全波形反演方法技术研究

论文摘要

地震波全波形反演是一种基于波动方程所描述的关于地震波在地下介质中传播过程,借助于地震波形所包含的丰富信息得到地下介质的相关参数的方法。其核心思想是利用观测到的数据和模拟数据的最优匹配进行地下介质的模型的建立。事实上,全波形反演的实现目前还存在着诸多困难,比如反演问题解的不稳定性;反演问题的强非线性;所处理数据的海量性和时间空间的高复杂度等。因此,提高地震波全波形反演的效率和精度是非常有理论意义和实用价值的工作。正演数值模拟是全波形反演的重要理论基础,其结果的准确性影响着地震数据处理的后期工作。由于计算机发展水平的提高,波动方程正演模拟方法发展较快。由于此方法不限制模型结构,一方面保持了地震波的运动学特征,另一方面保持了地震波的动力学特征,因此在正演模拟中有着广泛应用。波动方程正演可以在时间域进行,也可以在频率域进行。在时间域弹性波正演时,本文利用二维各向同性介质中的一阶速度-应力弹性波方程,推导了不同差分精度条件下的有理格式差分系数和相应的差分公式,分析了在不同差分精度条件下关于高阶交错网格的稳定性条件。利用高阶交错网格有限差分方法,在不同的差分精度条件下对波场进行了数值模拟,从波场快照、合成地震记录和CPU时间等方面对数值模拟结果进行了分析和对比,指出时间四阶和空间十二阶的差分精度具有较高计算精度和计算效率。在频率弹性波正演时,作为一种改进,本文在频率域25点差分格式的基础上分别提出了17点差分格式和15点差分格式的正演方法,重新计算了差分系数和频散条件,并对比了二者的阻抗矩阵及数值解。为提升效率,还给出了弹性波频率正演的多炮多核并行算法。最后,通过复杂模型的实验表明,在精度基本不变的前提下,采用并行策略的频率域15点格式计算效率更高。在反演部分,首先将声波方程看作弹性波方程的简化,给出了基于声波方程在时间域和频率域的反演算法和算例。在时间域声波全波形反演部分,采用基于L-BFGS优化算法的多尺度反演方法,其核心思想是使用低通滤波器将将观测数据,分解到不同的频带上,从低频带开始反演,将反演结果作为高频带的初始模型,然后用高频段反演出模型构造的细节部分,逐步提高成像的分辨率,从而提高反演效率。之后用Marmousi模型进行了反演实验,通过反演结果和单道速度信息的对比,证明该反演算法具有较高的反演精度。在声波频率全波形反演部分,给出了频率域多尺度的反演算法。核心思想是将全部频率分为几个频段分别反演,该方法使用较少的频点就可以完成反演过程,且能达到比使用全部频点一次反演还要好的效果。用简单模型进行的反演实验验证了算法的有效性,用Marmousi模型反演的结果表明该算法对复杂模型也能够达到较好的反演效果。在弹性波全波形反演时,参考声波多尺度反演算法的思想,分别设计了基于L-BFGS方法的弹性波时间域和频率域的多尺度反演算法。之后利用简单模型和修改后的Marmousi2模型进行反演实验,证明两种算法对于简单模型和复杂模型都能达到不错的反演效果。文章最后还对弹性波频率域和时间域的正演结果和反演初始模型、成像精度、反演效率等进行了对比和分析。

论文目录

摘要

Abstract

第一章绪论

1.1 论文的研究背景

1.2 国内外研究现状

1.2.1 时间域正演研究进展

1.2.2 频率域正演研究进展

1.2.3 时间域弹性波全波形反演发展

1.2.4 频率域弹性波全波形反演发展

1.3 本文主要研究内容及创新点

1.4 论文的框架结构

第二章弹性波时间域正演?

2.1 速度-应力方程高阶差分公式推导

2.1.1 时间上的4阶差分近似

2.1.2 2M-1阶导数的2N阶精度差分系数计算

2.1.3 2M阶导数的2N阶精度差分系数计算

2.1.4 混合偏导数的2N阶精度系数计算

2.2 稳定性条件

2.3 算例

2.3.1 模型参数设置

2.3.2 数值模拟结果

2.4 小结

第三章弹性波频率域正演

3.1 两种差分格式

3.1.1 二阶时间域弹性波方程

3.1.2 频率域17点有限差分格式

3.1.3 频率域15点有限差分格式

3.2 加权系数、网格频散、吸收边界

3.2.1 17点格式加权系数

3.2.2 15点格式加权系数

3.2.3 频散分析

3.2.4 PML吸收边界

3.3 数值模拟

3.3.1 17点格式数值模拟

3.3.2 15点格式数值模拟

3.3.3 复杂模型数值模拟

3.4 弹性波频率域并行计算

3.4.1 弹性波频率域多炮并行

3.4.2 弹性波频率域多炮多核并行

3.5 小结

第四章声波方程时间域多尺度全波形反演

4.1 声波方程的有限差分格式正演

4.1.1 二阶时间导数的2M阶差分展开

4.1.2 二阶空间导数的2N阶差分展开

4.1.3 时间域二阶声波方程PML边界构造

4.2 反演优化算法

4.2.1 共轭梯度法

4.2.2 步长的搜索条件

4.2.3 拟牛顿算法及L-BFGS算法

4.3 声波时间域全波形反演

4.3.1 反演基本流程

4.3.2 时间域方程梯度计算

4.3.3 时间域多尺度的反演策略

4.3.4 简单高速体模型时间域多尺度反演

4.3.5 Marmousi模型时间域多尺度反演

4.4 小结

第五章声波方程频率域多尺度全波形反演

5.1 频率域声波方程的有限差分法

5.1.1 常规9点差分格式

5.1.2 扩展的25点差分格式

5.1.3 四阶精度17点差分格式

5.1.4 优化15点差分格式

5.1.5 频率域声波方程PML边界构造

5.2 频率域声波方程全波形反演

5.2.1 反演一般流程

5.2.2 反演过程的梯度计算

5.2.3 频率域多尺度反演策略

5.2.4 双异常体模型反演

5.2.5 凹陷模型反演

5.2.6 正断层模型反演

5.2.7 Marmousi模型反演

5.3 小结

第六章弹性波方程时间域多尺度反演

6.1 弹性波时间域全波形反演理论推导

6.1.1 梯度的计算

6.2 Marmousi2模型时间域的多尺度反演

6.2.1 Marmousi2模型的预处理

6.2.2 Marmousi2模型的反演结果

6.3 小结

第七章弹性波方程频率域多尺度反演

7.1 弹性波频率域多炮并行正演

7.2 梯度的计算

7.3 反演频段的选择

7.4 简单异常体模型多尺度反演

7.5 Marmousi2模型多尺度反演

7.5.1 小尺寸复杂模型多尺度反演

7.5.2 Marmousi2模型多尺度反演

7.6 小结

第八章时间域和频率域正演、反演对比

8.1 弹性波时间域和频率域正演对比分析

8.2 反演的初始模型要求

8.3 反演精度和效率对比

8.3.1 声波时间域和频率域反演精度和效率对比

8.3.2 弹性波时间域和频率域反演精度和效率对比

第九章结论和展望

9.1 主要结论

9.1.1 关于正演

9.1.2 关于反演

9.2 对未来工作展望

参考文献

攻读学位期间的学术成果

致谢

文章来源

类型: 博士论文

作者: 岳晓鹏

导师: 白超英

关键词: 弹性波方程,全波形反演,时间域,频率域,多尺度反演

来源: 长安大学

年度: 2019

分类: 基础科学,工程科技Ⅰ辑

专业: 地质学,地球物理学,矿业工程

单位: 长安大学

分类号: P631.4

总页数: 178

文件大小: 16925K

下载量: 471

全波形反演方法技术研究

论文摘要

论文目录

文章来源

相关论文文献

猜你喜欢