部分6-全波形反演.pptVIP

二、频率-空间域全波形反演 Overthrust model test 右下图所示为原始3D 逆掩断层模型的一个侧面，目标区为右侧河道、中间断层及深部薄层。 为了测试近地表速度对反演的影响，在其余参数相同的情况下，分别采用表层速度与真速度相同的初始模型与表层速度与真速度不同的速度模型进行速度反演。 通过反演结果的差异可以看到表层速度的重要性，同时通过阻尼波场的引入可以看到反演效果的改善。 二、频率-空间域全波形反演 表层速度正确时的反演结果 利用如左下图所示的初始模型，利用频率空间域单频迭代波形反演，所得结果如右下图所示，可以看到右侧河道、中部断层及深部薄层都能基本得到。 二、频率-空间域全波形反演 表层速度不准确时的反演结果 利用如左下图所示的初始模型，表层速度与真实速度场并不相同，右下图为反演结果，可以看到若表层速度不正确时，不做任何处理，常规的反演算法并不能得到理想的结果。 基于阻尼波场的层剥离方法 左下图为阻尼系数为5的反演结果，可以看到模型的左侧和右侧已相对准确的反演出来，而中部断层则趋于局部极小值；右下图为以左下图为初始模型进行的阻尼系数为1.25的波形反演，可以看到左下图中表层速度反演正确部分的下部正确反演，而对于中部错误区域，下部的反演结果也受到影响。 二、频率-空间域全波形反演 二、频率-空间域全波形反演 4、小 结 1、频率-空间域全波形的理论推导基于线性代数的相关知识，推导较为简单，但可证明其与基于泛函分析的求导方法是等价的； 2、频率-空间域波形反演在处理小尺度问题时，利用线性方程组的直接求解方法可以方便的处理多炮的问题；在处理大尺度问题时，该优势不再显著； 3、频率-空间域波形反演方法可以只对频带中的几个频率进行反演，可以看做是一种压缩域的反演； 4、当从低频到高频逐次反演时，频率-空间域波形反演可看做是一种多尺度的反演方法； 5、Laplace-Fourier域反演方法可以看做是一种层剥离的反演方法。 提 纲 一、前言 二、频率-空间域全波形反演 三、时间域全波形反演 1、方法原理与数值实现 2、模型试算 3、小结 四、结论与认识 三、时间域全波形反演 1、方法原理与数值实现 首先定义模型空间 及数据空间 ，则模型参数属于模型空间，观测数据属于数据空间，即： 通过正演算子可以将模型参数映射到观测数据，即： 分别定义模型空间和数据空间的标量乘算子为： 三、时间域全波形反演 最小二乘意义下的最优模型就是寻求如下L2范数的最小值： 基于最速下降法可以得到其模型更新过程为： 其中： 为更新步长; 为预处理算子; 为Frechet 求导算子的转置算子; 如何构建Frechet求导算子及其转置算子？ 三、时间域全波形反演 Frechet求导算子及其转置算子 正演算子 进行Talyor展开得到： Frechet求导算子即为联系模型扰动与波场扰动的线性算子。 对于任意的 和 其转置算子满足下式： 通过上式可以进一步求出转置算子。 三、时间域全波形反演 对于声波介质全波形反演，其正问题由波动方程控制： 对于声波全波形反演问题的Frechet求导算子及其转置算子推导 基于Born近似，可求得介质参数扰动与波场扰动之间的关系： 可以通过格林函数显示写出： 三、时间域全波形反演 将上式代入算子与转置算子之间的关系： 将标量乘算子显示写开，并进行进一步整理可得： 物理解释： 正向波场在地下某点散射之后的传播至接收点的波场。 物理解释： 扰动波场，在梯度计算中为剩余波场。 对于转置算子的物理解释： 将地下任一点视为绕射点得到的波场与扰动波场的零延迟互相关。 这种显示的计算方式，对于m个介质参数需要m+1次正演过程。 物理解释： 正向传播波场二阶时间导数。 物理解释： 剩余波场的反向传播。 转置算子的物理意义： 剩余波场的反向传播与正向波场二阶时间导数的零延迟互相关； 此时计算梯度向量只需要一次正向传播和一次反向传播。 三、时间域全波形反演 借助波场的互易性，并作相关推导可得： 三、时间域全波形反演 基于相位编码的全波形反演 与频率域全波形反演相比，时间域全波形反演的效率与炮点的个数成正比，特别是在炮数较多的情况下，时域反演算法计算量非常大。通过引入相位编码技术可在一定基础上降低反演计算效率，特别是3D情况下效果更明显。 通过相位编码将多炮记录组成一超道集，在提高计算效率的同时，也引入了不同炮之间的串扰噪声。通过在迭代的过程中改变编码函数，即应用动态编码技术，可以在一定程度上降低串扰噪声的干扰。 下面两个示意图，展示了传统全波形反演与基于相位编码的全波形反演的异同。 三、时间域全波形反演 shot1 shot1 shot1 shot1 sh