一种基于分层融合的地震低频分量恢复方法技术

本发明专利技术提供了一种基于分层融合的地震低频分量恢复方法，通过正演方式获得全频段地震记录；将全频段地震记录滤除低频成分，得到对应的带限地震记录；将带限地震记录按照能量相近性做分层处理；将分层处理后的带限地震记录划分为训练集、验证集和测试集；利用训练集、验证集训练设计好的卷积神经网络，获得训练好的分层网络；将测试集带限数据输入训练好的分层网络，得到预测的分层的全频地震记录；将预测的分层的全频地震记录融合为完整的预测的地震记录；对完整的预测的地震记录做全波形反演，得到预测速度场。可有效避免由地震数据能量不均衡导致的卷积神经网络参数学习受干扰的现象，且恢复较深层低频分量的效果更优。且恢复较深层低频分量的效果更优。且恢复较深层低频分量的效果更优。

【技术实现步骤摘要】
一种基于分层融合的地震低频分量恢复方法


[0001]本专利技术属于全波形反演领域中地震低频分量恢复的方法，具体涉及到一种基于分层融合的地震低频分量恢复方法。

技术介绍

[0002]随着石油勘探的不断深入以及对高精度成像要求的不断提高，如何获取一个精细的地下介质模型成为当今地震勘探学家们思考的问题。地震勘探是通过人工方法在地下介质中激发地震波，地震波经过地下介质的反射、透射、折射等传播到布设在地面的检波器并被记录下来，形成地震记录。依据地震波在介质中的传播规律，设计相应的数值算法对该地震记录进行处理解释，得到地下介质的物性信息。作为现代地震学的杰出代表，全波形反演技术(Full waveform inversion)直接利用观测记录的全波信息重建地下介质的速度和密度等模型，可得地下介质的物性信息，是目前分辨率最高的一种精细建模技术。由于野外地震信息采集受限制，采集到的观测记录频率通常在3HZ以上，低频信息的缺失会严重影响全波形反演的精度。地震子波在地层向下传播时，会产生衰减，深度越大，采集到的地震记录能量越少。
[0003]卷积神经网络可以自动挖掘高频地震记录与相应全频地震记录之间的关系，被广泛应用于全波形反演领域的低频拓展中。但是由于地震记录的能量不均衡性，在训练神经网络的过程中，上层能量强的数据会对权重参数的学习贡献更大，对恢复下层能量低的数据起到一个反作用。因此，怎样训练卷积神经网络模型，能恢复出低频地震信息是目前亟待解决的问题。

技术实现思路

[0004]本专利技术提供了一种基于分层融合的地震低频分量恢复方法，克服训练卷积神经网络恢复全频地震记录存在的低能量数据被高能量数据抑制，深层数据难以恢复的问题。
[0005]本专利技术提供的一种基于分层融合的地震低频分量恢复方法，依次包括以下步骤：
[0006]S1、通过正演方式获得全频段地震记录；
[0007]S2、将全频段地震记录滤除低频成分，得到对应的带限地震记录；
[0008]S3、将带限地震记录按照能量相近性做分层处理，为减弱边缘效应，可使相连的两层设置重叠部分；
[0009]S4、将分层处理后的带限地震记录划分为训练集、验证集和测试集；
[0010]S5、利用训练集、验证集训练设计好的卷积神经网络，获得训练好的分层网络；
[0011]S6、将测试集带限数据输入训练好的分层网络，得到预测的分层的全频地震记录；
[0012]S7、将预测的分层的全频地震记录融合为完整的预测的地震记录，重叠部分在融合时给与不同权重；
[0013]S8、对完整的预测的地震记录做全波形反演，得到预测速度场。
[0014]优选的，步骤1正演设置包括地震子波设置为主频为10hz的雷克子波；所用地质模
型选择像素点数为155*461的Marmousi模型，设置炮数为112，检波点数为461；所用地震波数值模拟方法为有限差分法，使用时域2阶、空间10阶的PML有限差分建模算子，采样时间间隔为0.8ms，采样点数为4001。
[0015]优选的，步骤3所述分层处理将带限地震记录按照能量相近性共划分了五层。
[0016]优选的，步骤5所设计的卷积神经网络的工作流程为：
[0017]S1、输入分层的带限地震记录；
[0018]S2、将分层的带限地震记录经过一层二维卷积，卷积核大小为640*1，滤波器个数为128，使用PReLu激活函数，然后将得到的数据进行批处理归一化；
[0019]S3、步骤2结果经过一层Dropout层，以概率50％舍弃部分神经元；
[0020]S4、步骤3结果经过一层二维卷积，卷积核大小为640*1，滤波器个数为64，使用PReLu激活函数，然后将得到的数据进行批处理归一化；
[0021]S5、步骤4结果经过一层二维卷积，卷积核大小为640*1，滤波器个数为128，使用PReLu激活函数，然后将得到的数据进行批处理归一化；
[0022]S6、步骤5结果经过一层二维卷积，卷积核大小为640*1，滤波器个数为64，使用PReLu激活函数，然后将得到的数据进行批处理归一化；
[0023]S7、步骤6结果经过一层二维卷积，卷积核大小为640*1，滤波器个数为32，使用PReLu激活函数，然后将得到的数据进行批处理归一化；
[0024]S8、步骤7结果经过一层全连接层，神经元个数等于每层对应的采样点数；
[0025]S9、得到卷积神经网络预测出的该层对应的全频地震记录。
[0026]本专利技术的有益效果：本专利技术所提出的基于分层融合恢复低频信息的方法可有效避免由地震数据能量不均衡导致的卷积神经网络参数学习受干扰的现象，且恢复较深层低频分量的效果更优。
附图说明
[0027]图1所示为基于分层融合的地震低频分量恢复方法的流程图，
[0028]图2所示为正演方式获得全频段地震记录所需设置的参数示意图，
[0029]图3所示为本专利技术将全频段地震记录滤除低频成分的流程图，
[0030]图4所示为本专利技术实施例中验证集分层处理示意图，
[0031]图5所示为本专利技术设计的卷积神经网络的工作流程图，
[0032]图6所示为本专利技术实施例步骤6的流程图，
[0033]图7所示为本专利技术实施例预测的完整的地震记录的全频图，
[0034]图8所示为未采用本专利技术实施方法预测的完整的地震记录的全频图，
[0035]图9所示为本专利技术实施例方法所得结果所对应的全频地震记录的真实标签图。
具体实施方式
[0036]为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用于解释本专利技术，不能理解为对本专利技术具体保护范围的限定。
[0037]实施例
[0038]本实施例提供了一种基于分层融合的地震低频分量恢复方法，如图1所示，依次包括如下步骤：
[0039]S1、通过正演方式获得全频段地震记录；
[0040]S2、将全频段地震记录滤除低频成分，得到对应的带限地震记录；
[0041]S3、将带限地震记录按照能量相近性做分层处理，为减弱边缘效应，可使相连的两层设置重叠部分；
[0042]S4、将分层处理后的带限地震记录划分为训练集、验证集和测试集；
[0043]S5、利用训练集、验证集训练设计好的卷积神经网络，获得训练好的分层网络；
[0044]S6、将测试集带限数据输入训练好的分层网络，得到预测的分层的全频地震记录；
[0045]S7、将预测的分层的全频地震记录融合为完整的预测的地震记录，重叠部分在融合时给与不同权重；
[0046]S8、对完整的预测的地震记录做全波形反演，得到预测速度场。
[0047]步骤1采用现有C语言编程的正演程序去获得全频段地震记录，具体需要设置以下参数，如图2所示，所用地震子波为主频为10hz本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于分层融合的地震低频分量恢复方法，其特征在于，依次包括以下步骤：S1、通过正演方式获得全频段地震记录；S2、将全频段地震记录滤除低频成分，得到对应的带限地震记录；S3、将带限地震记录按照能量相近性做分层处理，为减弱边缘效应，可使相连的两层设置重叠部分；S4、将分层处理后的带限地震记录划分为训练集、验证集和测试集；S5、利用训练集、验证集训练设计好的卷积神经网络，获得训练好的分层网络；S6、将测试集带限数据输入训练好的分层网络，得到预测的分层的全频地震记录；S7、将预测的分层的全频地震记录融合为完整的预测的地震记录，重叠部分在融合时给与不同权重；S8、对完整的预测的地震记录做全波形反演，得到预测速度场。2.根据权利要求1所述的基于分层融合的地震低频分量恢复方法，其特征在于，步骤1正演设置包括地震子波设置为主频为10hz的雷克子波；所用地质模型选择像素点数为155*461的Marmousi模型，设置炮数为112，检波点数为461；所用地震波数值模拟方法为有限差分法，使用时域2阶、空间10阶的PML有限差分建模算子，采样时间间隔为0.8ms，采样点数为4001。3.根据权利要求1所述的基于分层融合的地震低频分量恢复方法，其特征在于：步骤3所述分层处理将...

【专利技术属性】
技术研发人员：刘瑶，刘伟锋，黄建平，刘宝弟，王珺，齐玉娟，
申请(专利权)人：中国石油大学华东，
类型：发明
国别省市：