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

1.本发明属于全波形反演领域中地震低频分量恢复的方法，具体涉及到一种基于分层融合的地震低频分量恢复方法。


背景技术：

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


技术实现要素：

4.本发明提供了一种基于分层融合的地震低频分量恢复方法，克服训练卷积神经网络恢复全频地震记录存在的低能量数据被高能量数据抑制，深层数据难以恢复的问题。
5.本发明提供的一种基于分层融合的地震低频分量恢复方法，依次包括以下步骤：
6.s1、通过正演方式获得全频段地震记录；
7.s2、将全频段地震记录滤除低频成分，得到对应的带限地震记录；
8.s3、将带限地震记录按照能量相近性做分层处理，为减弱边缘效应，可使相连的两层设置重叠部分；
9.s4、将分层处理后的带限地震记录划分为训练集、验证集和测试集；
10.s5、利用训练集、验证集训练设计好的卷积神经网络，获得训练好的分层网络；
11.s6、将测试集带限数据输入训练好的分层网络，得到预测的分层的全频地震记录；
12.s7、将预测的分层的全频地震记录融合为完整的预测的地震记录，重叠部分在融合时给与不同权重；
13.s8、对完整的预测的地震记录做全波形反演，得到预测速度场。
14.优选的，步骤1正演设置包括地震子波设置为主频为10hz的雷克子波；所用地质模
型选择像素点数为155*461的marmousi模型，设置炮数为112，检波点数为461；所用地震波数值模拟方法为有限差分法，使用时域2阶、空间10阶的pml有限差分建模算子，采样时间间隔为0.8ms，采样点数为4001。
15.优选的，步骤3所述分层处理将带限地震记录按照能量相近性共划分了五层。
16.优选的，步骤5所设计的卷积神经网络的工作流程为：
17.s1、输入分层的带限地震记录；
18.s2、将分层的带限地震记录经过一层二维卷积，卷积核大小为640*1，滤波器个数为128，使用prelu激活函数，然后将得到的数据进行批处理归一化；
19.s3、步骤2结果经过一层dropout层，以概率50％舍弃部分神经元；
20.s4、步骤3结果经过一层二维卷积，卷积核大小为640*1，滤波器个数为64，使用prelu激活函数，然后将得到的数据进行批处理归一化；
21.s5、步骤4结果经过一层二维卷积，卷积核大小为640*1，滤波器个数为128，使用prelu激活函数，然后将得到的数据进行批处理归一化；
22.s6、步骤5结果经过一层二维卷积，卷积核大小为640*1，滤波器个数为64，使用prelu激活函数，然后将得到的数据进行批处理归一化；
23.s7、步骤6结果经过一层二维卷积，卷积核大小为640*1，滤波器个数为32，使用prelu激活函数，然后将得到的数据进行批处理归一化；
24.s8、步骤7结果经过一层全连接层，神经元个数等于每层对应的采样点数；
25.s9、得到卷积神经网络预测出的该层对应的全频地震记录。
26.本发明的有益效果：本发明所提出的基于分层融合恢复低频信息的方法可有效避免由地震数据能量不均衡导致的卷积神经网络参数学习受干扰的现象，且恢复较深层低频分量的效果更优。
附图说明
27.图1所示为基于分层融合的地震低频分量恢复方法的流程图，
28.图2所示为正演方式获得全频段地震记录所需设置的参数示意图，
29.图3所示为本发明将全频段地震记录滤除低频成分的流程图，
30.图4所示为本发明实施例中验证集分层处理示意图，
31.图5所示为本发明设计的卷积神经网络的工作流程图，
32.图6所示为本发明实施例步骤6的流程图，
33.图7所示为本发明实施例预测的完整的地震记录的全频图，
34.图8所示为未采用本发明实施方法预测的完整的地震记录的全频图，
35.图9所示为本发明实施例方法所得结果所对应的全频地震记录的真实标签图。
具体实施方式
36.为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用于解释本发明，不能理解为对本发明具体保护范围的限定。
37.实施例
38.本实施例提供了一种基于分层融合的地震低频分量恢复方法，如图1所示，依次包括如下步骤：
39.s1、通过正演方式获得全频段地震记录；
40.s2、将全频段地震记录滤除低频成分，得到对应的带限地震记录；
41.s3、将带限地震记录按照能量相近性做分层处理，为减弱边缘效应，可使相连的两层设置重叠部分；
42.s4、将分层处理后的带限地震记录划分为训练集、验证集和测试集；
43.s5、利用训练集、验证集训练设计好的卷积神经网络，获得训练好的分层网络；
44.s6、将测试集带限数据输入训练好的分层网络，得到预测的分层的全频地震记录；
45.s7、将预测的分层的全频地震记录融合为完整的预测的地震记录，重叠部分在融合时给与不同权重；
46.s8、对完整的预测的地震记录做全波形反演，得到预测速度场。
47.步骤1采用现有c语言编程的正演程序去获得全频段地震记录，具体需要设置以下参数，如图2所示，所用地震子波为主频为10hz的雷克子波；所用地质模型(速度场)为像素点数为155*461的marmousi模型，设置炮数为112，检波点数为461；所用地震波数值模拟方法为有限差分法，使用时域2阶、空间10阶的pml有限差分建模算子，采样时间间隔为0.8ms，采样点数为4001。
48.步骤2将全频段地震记录滤除低频成分的流程如图3所示，设计一个频域滤波器，滤除全频段地震记录中频率在5hz以下的数据(即低频成分)，得到对应的带限地震记录。
49.步骤3所述的将带限地震记录和全频地震记录按照能量相近性做分层处理，如图4所示的是4炮地震记录时域图，从左到右依次为第一炮的地震记录，第一炮与中间炮之间的中间炮地震记录，中间炮地震记录，最后一炮地震记录，此4炮地震记录大致可体现全部地震记录的特征。地震记录会随地层深度的增大产生不同程度的衰减，根据地震记录的能量相近性，按图4中横线所示将完整的地震记录分为5层或者其他层数，相连的两层可设置重叠部分，以此减弱边缘效应。本实施例通过5层分层处理，从而获取分层处理后的带限、全频地震记录。
50.本实施例的全频地震记录共有112炮，我们将奇数炮数据作为训练集，偶数炮数据作为测试集，从偶数炮里抽取4炮作为验证集，图4所示的即为验证集数据。
51.利用训练集、验证集训练设计好的卷积神经网络，获得训练好的分层网络的流程如图5所示，步骤1501是起始状态；步骤1502输入分层的带限地震记录；步骤1503将带限地震记录经过一层二维卷积，卷积核大小为640*1，滤波器个数为128，使用prelu激活函数，然后将得到的数据进行批处理归一化，prelu激活函数的公式为：
[0052][0053]
其中，y是非线性激活函数g(α,y)的输入，α是一个可以学习的参数，当输入大于等于零时，prelu就相当于relu，输出等于输入，当输入小于零时，要将输入乘以一个参数α得到输出；步骤1504将上一步结果经过一层dropout层，以概率50％舍弃部分神经元；步骤1505将上一步结果经过一层二维卷积，卷积核大小为640*1，滤波器个数为64，使用prelu激活函数，然后将得到的数据进行批处理归一化；步骤1506将上一步结果经过一层二维卷积，
卷积核大小为640*1，滤波器个数为128，使用prelu激活函数，然后将得到的数据进行批处理归一化；步骤1507将上一步结果经过一层二维卷积，卷积核大小为640*1，滤波器个数为64，使用prelu激活函数，然后将得到的数据进行批处理归一化；步骤1508将上一步结果经过一层二维卷积，卷积核大小为640*1，滤波器个数为32，使用prelu激活函数，然后将得到的数据进行批处理归一化；步骤1509将上一步结果经过一层全连接层，神经元个数等于每层对应的采样点数；步骤1510得到卷积神经网络预测出的该层对应的全频地震记录；步骤1511是结束状态。
[0054]
如图6所示，本实施例分别将五层带限地震记录的测试数据放入对应的分层网络，得到这5层的预测的全频地震记录，然后将这5层预测的全频地震记录融合为完整的预测的地震记录，重叠部分在融合时给与不同权重。对完整的预测的地震记录做全波形反演，得到预测速度场，将预测速度场与步骤1正演时所用的速度场进行比较，看最终的预测结果的准确性。
[0055]
图7给出了使用本实施例基于分层融合的地震低频分量恢复方法预测的完整的地震记录的全频图。图7所示结果未设置重叠区域，使用fimage软件画图，便于与现有方法的结果直观对比。
[0056]
图8给出了对应的未使用本实施例基于分层融合的地震低频分量恢复方法预测的全频图。具体步骤如下：
[0057]
s1、通过正演方式获得全频段地震记录；
[0058]
s2、将全频段地震记录划分为训练集、验证集和测试集；
[0059]
s3、利用训练集、验证集训练设计好的卷积神经网络；
[0060]
s4、将测试集带限数据放入训练好的卷积神经网络，得到预测的全频地震记录；
[0061]
s5、对预测的全频地震记录做全波形反演，得到预测速度场。
[0062]
图8使用的实验数据与图7所采用的实验数据保持一致，训练集、验证集、测试集的划分方式也保持一致，只是数据不曾经过图1中步骤s3所示的分层处理，用完整的训练集训练相同的卷积神经网络，验证集的预测结果为完整的全频地震记录，将其用fimage画图，与图7对比。
[0063]
图9给出了本实施例方法所得结果所对应的全频地震记录的真实标签图。将验证集对应的标签用fimage画图，与图7、图8对比。结果显示，图7的的下层恢复效果较好，上层有少量噪声，恢复效果较下层稍差；图8的下层恢复结果差，上层恢复效果也不如图7，且图8整体噪声明显。实验结果表明，本发明所提出的基于分层融合恢复低频信息的方法可有效避免由地震数据能量不均衡导致的卷积神经网络参数学习受干扰的现象，且恢复较深层低频分量的效果更优。