跳至主要內容
Accelerating Full Waveform Inversion By Transfer Learning

Accelerating Full Waveform Inversion By Transfer Learning

  • DOI:10.1007/s00466-025-02600-w

  • 全波形反演(FWI)是一种强大的工具,能够基于波传播获得的稀疏测量数据来重建物质场。对于特定问题,使用神经网络(NN)对物质场进行离散化可以提高相应优化问题的鲁棒性和重建质量。本文将这种方法称为基于神经网络的FWI。从初始预测开始,迭代更新神经网络的权重,以使模拟的波形信号与稀疏测量数据集相拟合。对于基于梯度的优化,选择合适的初始预测(即合适的神经网络权重初始化)对于实现快速且稳健的收敛至关重要。

    本文提出了一种新颖的迁移学习方法,以进一步改进基于神经网络的FWI。该方法利用监督预训练来提供更好的神经网络权重初始化,从而加快后续优化问题的收敛速度。此外,反演结果还能得到物理上更有意义的局部最小值。网络被预训练以使用传统FWI第一次迭代中的梯度信息来预测未知的物质场。在二维域计算实验中,训练数据集由具有任意位置、不同形状和方向的椭圆形空洞的参考模拟组成。

    本文将所提出的迁移学习基于神经网络的FWI与另外三种方法进行了性能比较:传统FWI、未进行预训练的基于神经网络的FWI,以及使用预训练神经网络预测的初始猜测的传统FWI。结果表明,迁移学习基于神经网络的FWI在收敛速度和重建质量方面均优于其他方法。


Xenny原创大约 7 分钟FWIFWI迁移学习
Deep Learning for Low-Frequency Extrapolation of Multicomponent Data in Elastic FWI

Deep Learning for Low-Frequency Extrapolation of Multicomponent Data in Elastic FWI

  • DOI: 10.1109/TGRS.2021.3135790

  • FWI的成功强依赖于一个准确的初始模型,在弹性体制下尤为重要——由于S-波波长较短,弹性FWI中的周期跳跃现象比声波FWI更严重。本文中,作者通过提出合成多分量弹性低频地震记录,并将这些“人工” 低频数据作为弹性FWI的频率扫描种子,进一步扩展了外推弹性FWI的工作。作者的方案包含深度学习——可以在两个训练数据集上训练同一个卷积神经网络(CNN),一个是粒子速度的垂直分量,另一个是水平分量,或将两个分量放在一起训练,以推断2D 弹性 FWI的低频弹性数据。CNN架构通过空洞卷积获得大感受野。在Marmousi2上的实验表明从44 Hz 以上的带限(band-limited)数据推断出的 242-4 Hz 低频数据为P波和S波速度的弹性 FWI 提供了良好的初始模型。 此外,我们还研究了网络从声波数据到弹性数据的泛化能力,在弹性测试数据上, 通过弹性模拟收集训练数据集比声学模拟显示出更好的外推精度,即更小的泛化差距。


Xenny原创大约 19 分钟FWIEFWIFWI深度学习
CAUC: Combining Channel Attention U-Net and Convolution for Seismic Data Resolution Improvement

CAUC: Combining Channel Attention U-Net and Convolution for Seismic Data Resolution Improvement

  • DOI:10.1109/LGRS.2023.3322263

  • 地震分析和解释对数据分辨率很敏感。由于复杂的自然环境和有限的采集技术,原始地震数据分辨率通常较低。但地震子波的频率带宽是可调的,所以根据这些特征来提高数据的分辨率是可行的。本文提出了一种通道注意U-Net和物理卷积组合(CAUC)算法来增强地震数据的分辨率。


Xenny原创大约 4 分钟FWIFWIU-Net
地震反演基本知识

地震反演基本知识

基本概念

  1. 地震数据

    波向地下传播时,由于地下结构不同导致地震波会反射回地表,通过在地表设置检波器收集到的波数据即为地震数据。

  2. 速度模型

    在地震反演中,我们用速度模型来描述地下地质结构,不同的介质波的传播速度不一样,故而可以由速度模型反推介质成分。

  3. 波阻抗(Impedance)

    波阻抗等于速度乘密度,波阻抗越大代表要产生单位振动速度所需的应力越大。

  4. 反射系数

    反射系数代表地震波在地下不同介质界面上发生折射、反射的能量损失程度。计算公式为

    ri=AIi+1AIiAIi+1+AIi,(1) r_i = \frac{AI_{i+1} - AI_{i}}{AI_{i+1} + AI_i},\tag{1}

    其中AIAI代表波阻抗。

  5. 地震褶积

    其实就是子波和反射系数进行卷积运算,一般是在时间域进行运算。本质是把反射系数序列作为一个函数,是地震波传播的数字信号模拟。可以理解为地震波在地下垂直传播过程,在经过不同界面时和反射系数进行叠加。


Xenny原创大约 1 分钟FWIFWI
地震波阻抗反演实验

地震波阻抗反演实验

准备工作

数据集

  • 我选用Convolutional neural network for seismic impedance inversion中的数据,其中包含2020道一维地震波以及波阻抗数据。

    图1. 数据展示

    正演参数为频率为30Hz的雷克子波,Δt=4.3875e4\Delta t = 4.3875e-4。图1中为数据集中的4道数据展示,其中波阻抗数据已进行缩放。

  • 随后将这些数据集划分为500个验证集、500个测试集和1020个训练集。


Xenny大约 4 分钟FWIFWICNN深度学习
SG-PML吸收边界条件(边界反射问题)

SG-PML吸收边界条件(边界反射问题)

  • 本文将对博文 波动方程的有限差分解(地震数据正演) 中提到的边界反射问题进行研究。解决地震正演中的边界反射带来的干扰。

    首先我们要知道波为什么会反射,当从一种介质进入另一种介质时,在介质边界由于阻抗的不连续性,导致一部分能量继续进入新的介质称为透射波,一部分能量返回原介质称为反射波。而我们将需要了解反射部分具体是多少,以及如何消除边界反射问题。

吸收边界条件

  • 自然界中地震波传播到介质边界(中断层)时,界面的波阻抗出现差异导致波发生反射和透射,地球介质可以看作“无限大”的区域,但在模拟计算中受到计算量的限制,模拟的传播区域大小是由界限的,这个人为截取的区域便是人工边界。

    在人工边界处,边界外的波速可以看为0,即边界反射系数绝对值为1,此时表现在数值模拟中便产生了全反射现象,对地震正演中的波场计算将产生影响。此时我们需要引入吸收边界的概念来吸收掉这些传播至边界的波。

  • 吸收边界条件即人为的在模拟计算区域的边界处设置一定厚度的吸收层,让能量在吸收层呈指数曲线衰减,最终在吸收层边界处近似衰减到0而避免边界反射。

    注意这里我们的衰减函数也是有讲究的,因为如果衰减不够连续的话吸收层也将看成断层产生反射,虽然这里使用指数曲线进行衰减,但是实际模拟时还是会产生一定量的反射波。


Xenny原创大约 14 分钟FWIPMLFWI