一种基于自编码器和深度神经网络的全波形反演方法技术

本发明专利技术公开了一种基于自编码器和深度神经网络的全波形反演方法，包括以下步骤：构建输入矩阵；利用去噪编码模块数据去噪；利用卷积自编码器模块进行特征提取；使用卷积自编码器模块的特征作为输入，利用深度神经网络正演模型进行波场值预测；利用预测的模拟值与真实地震数据记录，基于目标函数进行两者之间的偏差计算；利用优化算法进行梯度反向传播，进行神经网络的权重优化；将经过训练得到的最佳神经网络模型来提取反演效果。本发明专利技术可广泛应用于各类波形反演场景中。于各类波形反演场景中。

【技术实现步骤摘要】
一种基于自编码器和深度神经网络的全波形反演方法


[0001]本专利技术涉及人工智能与地球物理领域，具体涉及一种基于自编码器和深度神经网络的全波形反演算法。

技术介绍

[0002]全波形反演(Full Waveform Inversion，FWI)方法和理论是地球物理领域的研究热点之一，具有很高的实用价值。全波形反演是利用地震资料中的全部的动力学和运动学信息来进行反演地层构成模型的方法。在一些实际的地质资源勘探问题中，如油气储藏勘探、金属矿脉勘探、地质构成勘察等等，由于不同尺度的非均匀地质体共生，会引起各种反射波、绕射波及其对应的多次波、转换波等的相互叠加、干涉，可形成极为复杂的地震波场。FWI首先通过正演模拟获得模拟地震预测结果，然后使用目标函数度量模拟地震结果与实际观测地震结果之间的差异，再通过特定的优化算法来更新模型参数。
[0003]传统的FWI使用模拟值与真实值之间的偏差作为目标函数，并使用伴随状态法来进行梯度的计算，沿着负梯度方向不断更新正演模型的相关参数。其中，以L2范数作为目标函数的FWI方法被广泛地使用，但是其是一个复杂的非线性优化的问题，且在反演中为了减小运算量通常使用一些局部优化算法，若初始模型的选择不恰当，那么反演过程就容易遇到周期跳变(cycle skipping)问题，从而导致无法取得精确的反演结果。另外，地震记录中蕴含全部的运动学信息和动力学信息，而不同的信息之间往往相互耦合、相互影响并且存在不同程度的噪声的，如果不能对这些信息进行合理的特征提取和利用则会导致反演结果的精度进一步降低。
[0004]目前在地球物理领域，已有相当一部分的工作尝试使用各种深度学习的手段来进行各类传统问题的反演。现有的成熟的深度学习框架例如Pytorch针对GPU等底层设备进行了专有的运算加速，可以方便地进行自动的梯度计算和反向传播，为基于神经网络的反演模型提供了便利。并且这种自动求导和反演的过程与常用的伴随状态矩阵的求解结果具有等价性。深度学习领域完善的优化方法也相较于传统模型的局部搜索方法显示出了优越性。传统的全局优化方法由于运算复杂度过高而难以应用，而局部搜索方法又由于模型的非凸性和非线性而容易陷入局部最优。而深度学习领域的Adam等优化器可以基于动量等角度来进行梯度估计，从而可以在较小的运算成本下达到全局最优。

技术实现思路

[0005]针对现有技术的不足，本专利技术旨在提供一种基于自编码器和深度神经网络的全波形反演方法。
[0006]为了实现上述目的，本专利技术采用如下技术方案：
[0007]一种基于自编码器和深度神经网络的全波形反演方法，包括以下步骤：
[0008]步骤1、对原始地震数据进行处理，构建输入矩阵；
[0009]步骤2、利用去噪编码模块进行退化实现数据去噪，去除原始地震数据中的干扰信
息和噪声；
[0010]步骤3、利用卷积自编码器模块进行特征提取；
[0011]步骤4、使用卷积自编码器模块提取的特征作为输入，利用深度神经网络正演模型进行波场值预测；
[0012]步骤5、利用步骤4中预测得到的模拟值与真实地震数据记录的波场值进行比较，利用目标函数进行两者之间的偏差计算；
[0013]步骤6、根据步骤5所得的偏差值，利用优化算法进行梯度反向传播，进行神经网络的权重优化；
[0014]步骤7、将测试数据输入步骤6中经过训练得到的最佳神经网络模型，以提取反演效果；通过查询测试阶段的输出值和真实值之间的偏差，即可得到全波形预测的效果。
[0015]进一步地，步骤2的具体过程如下：利用去噪编码模块对步骤1所得的输入矩阵进行初步编码嵌入，并利用退化策略进行数据去噪；所述退化策略为利用随机dropout策略对神经元进行随机失活。
[0016]进一步地，步骤3的具体过程如下：
[0017]所述卷积自编码器模块利用卷积层和池化层作为编码器，利用上采样层和卷积层作为解码器进行输入重构，以提取准确的输入特征；假设X表示输入样本，W
k
与b
k
分别表示第k个卷积核的权值与偏置，*表示卷积运算，g(
·
)表示池化函数，h
k
表示第k个卷积核所提取的抽象特征，则编码过程为:
[0018]h
k
＝g(W
k
*X+b
k
)
ꢀꢀꢀ
#(1)
[0019]假设W
′
k
与b
′
k
分别表示解码器中第k个卷积核的权值与偏置，g
′
(
·
)表示上采样函数，H表示抽象特征集合，则解码过程为:
[0020][0021]表示重构的样本；解码器将各个卷积核提取的抽象特征进行解码重构，并将其合并为最终的重构样本；重构损失函数为:
[0022][0023]J(
·
)表示损失函数，CoAE表示自编码器；为权值L2范数正则化项，W表示所有卷积核的参数；λ表示衰减系数，用于控制权值的衰减程度。
[0024]进一步地，步骤4的具体过程如下：在步骤3训练好的卷积自编码器模块的基础上，重新将输入样本X输入卷积自编码器模块以提取特征编码矩阵h，将所获得的编码矩阵输入深度神经网络正演模型中进行波场值预测，得到模拟地震记录分布f。
[0025]进一步地，步骤5的具体过程如下：
[0026]通过考虑模拟值的累积分布F(t)和真实值的累积分布G(t)之间的最小二乘差异，使用积分波场目标函数对地震数据做线性正变换与积分归一化，计算公式为：
[0027][0028]分别表示F(t)和G(t)的密度函数，T0表示地震数据中的积分上界，故基于积分归一化度量的全波形反演目标函数为：
[0029][0030]其中t表示时间，θ表示模型参数，x
r
表示自变量，和分别表示关于s的密度函数，s和r分别表示两个自变量，N
s
和N
r
分别表示s和r的最大值。
[0031]进一步地，步骤6的具体过程为：利用Adam算法从矩估计的角度来自适应地调整学习率；利用梯度的一阶矩和二阶矩估计动态来调整学习率，给予每个参数的不同的学习率；
[0032]最终形成的Adam算法的公式为：
[0033][0034]t表示时间，w
t
表示为t时刻神经网络的权重，和是常数，α为控制指数衰减系数，∈为一极小值常数，其作用为避免分母为0；m
t
为t时刻梯度的一阶指数移动的值，v
t
为t时刻梯度的二阶指数移动的值，g
t
为为t时刻的梯度值，则有
[0035][0036]本专利技术的有益效果在于：本专利技术通过去噪自编码器模块中的随机退化思想来对原始地震数据中的蕴含的随机噪声进行去除，并利用卷积自编码器进行高级特征的提取和编码。本专利技术利用深度神经网络来构建波形正本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于自编码器和深度神经网络的全波形反演方法，其特征在于，包括以下步骤：步骤1、对原始地震数据进行处理，构建输入矩阵；步骤2、利用去噪编码模块进行退化实现数据去噪，去除原始地震数据中的干扰信息和噪声；步骤3、利用卷积自编码器模块进行特征提取；步骤4、使用卷积自编码器模块提取的特征作为输入，利用深度神经网络正演模型进行波场值预测；步骤5、利用步骤4中预测得到的模拟值与真实地震数据记录的波场值进行比较，利用目标函数进行两者之间的偏差计算；步骤6、根据步骤5所得的偏差值，利用优化算法进行梯度反向传播，进行神经网络的权重优化；步骤7、将测试数据输入步骤6中经过训练得到的最佳神经网络模型，以提取反演效果；通过查询测试阶段的输出值和真实值之间的偏差，即可得到全波形预测的效果。2.根据权利要求1所述的方法，其特征在于，步骤2的具体过程如下：利用去噪编码模块对步骤1所得的输入矩阵进行初步编码嵌入，并利用退化策略进行数据去噪；所述退化策略为利用随机dropout策略对神经元进行随机失活。3.根据权利要求1所述的方法，其特征在于，步骤3的具体过程如下：所述卷积自编码器模块利用卷积层和池化层作为编码器，利用上采样层和卷积层作为解码器进行输入重构，以提取准确的输入特征；假设X表示输入样本，W
k
与b
k
分别表示第k个卷积核的权值与偏置，*表示卷积运算，g(
·
)表示池化函数，h
k
表示第k个卷积核所提取的抽象特征，则编码过程为:h
k
＝g(W
k
*X+b
k
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
#(1)假设W
′
k
与b
′
k
分别表示解码器中第k个卷积核的权值与偏置，g
′
(
·
)表示上采样函数，H表示抽象特征集合...

【专利技术属性】
技术研发人员：马飞，朱红达，史浩，张海阳，郝骞，李斌，陈钢，
申请(专利权)人：北京华电力拓能源科技有限公司，
类型：发明
国别省市：