1.本发明涉及地震速度分析技术领域，尤其是涉及一种基于卷积神经网络深度学习的反射波形反演方法和系统。
背景技术：2.全波形反演(full waveform inversion，fwi)可以看作是走时层析和偏移技术的组合。然而，传统的fwi无法准确估计更深地层的背景速度模型。通过分离层析成像和偏移分量，反射波形反演(reflection waveform inversion，rwi)能够有效地利用走时层析分量来更新更深地层的背景速度模型。现有技术中，一般利用卷积神经网络(convolutional neural networks，cnn)对背景速度模型进行预测，但是传统用于预测速度模型的cnn仅仅是经过一次训练和一次预测过程，因而得到的预测结果与真实速度模型误差依然较大。
技术实现要素：3.有鉴于此，本发明的目的在于提供一种基于卷积神经网络深度学习的反射波形反演方法和系统，以缓解现有技术中存在的预测的速度模型不准确的技术问题。
4.第一方面，本发明实施例提供了一种基于卷积神经网络深度学习的反射波形反演方法，包括：对原始地震数据进行预处理，得到初始速度模型；区块化步骤：通过将所述初始速度模型区块化为空间相关且连接不相交的特征集，创建速度模型基；基于所述速度模型基，生成训练集；所述训练集包括训练速度模型集和逆时偏移图像集；所述逆时偏移图像集为对所述训练速度模型集进行逆时偏移成像得到的图像集；基于所述训练集对预设卷积神经网络进行训练，得到训练之后的卷积神经网络；以所述初始速度模型和初始逆时偏移图像作为输入，利用所述训练之后的卷积神经网络进行预测，得到预测速度模型；所述初始逆时偏移图像为对所述初始速度模型进行逆时偏移成像得到的图像；判断基于所述预测速度模型正演模拟的地震数据，与所述原始地震数据之间的误差是否小于预设阈值；如果否，则将预测速度模型作为新的初始速度模型，并返回执行所述区块化步骤。
5.进一步地，所述区块化步骤包括：基于k均值聚类，将所述初始速度模型的网格点划分为两个不相交的聚类簇；分别确定每个聚类簇中网格点的最大空间相连区域，得到两个相连区域；基于网格点与所述两个相连区域的曼哈顿距离，将所述初始速度模型区块化为空间相关且连接不相交的两个子模型；对所述两个子模型分别进行递归区块化处理，得到所述特征集；基于所述特征集，创建所述速度模型基。
6.进一步地，基于网格点与所述两个相连区域的曼哈顿距离，将所述初始速度模型区块化为空间相关且连接不相交的两个子模型，包括：将所述初始速度模型的网格点划分为三个区域，分别为所述两个相连区域和第三区域；所述第三区域为所述初始速度模型的网格点中除所述两个相连区域内的网格点之外的所有网格点组成的区域；分别计算所述第三区域内的每个网格点与所述两个相连区域的曼哈顿距离；将所述第三区域内的网格点合并到曼哈顿距离最小的相连区域，得到所述两个子模型。
7.进一步地，基于所述速度模型基，生成训练集，包括：基于所述速度模型基的每个分量中速度的平均值和方差，创建所述训练速度模型集；对所述训练速度模型集中的每个训练速度模型进行逆时偏移成像，得到所述逆时偏移图像集；将所述训练速度模型集和所述逆时偏移图像集作为所述训练集。
8.进一步地，基于所述训练集对预设卷积神经网络进行训练，包括：基于所述训练集并通过迭代最小化目标损失函数的方式，对所述预设卷积神经网络进行训练；其中，所述目标损失函数包括：g是所述预设卷积神经网络，w为所述预设卷积神经网络的参数化的权重，w
*
为所述训练之后的卷积神经网络的参数化的权重，m为所述训练集中训练速度模型的数量，为第k个初始训练速度模型，为第k个训练速度模型，为第k个逆时偏移图像。
9.第二方面，本发明实施例还提供了一种基于卷积神经网络深度学习的反射波形反演系统，包括：预处理模块，区块化模块，生成模块，训练模块，预测模块和迭代模块；其中，所述预处理模块，用于对原始地震数据进行预处理，得到初始速度模型；所述区块化模块，用于通过将所述初始速度模型区块化为空间相关且连接不相交的特征集，创建速度模型基；所述生成模块，用于基于所述速度模型基，生成训练集；所述训练集包括训练速度模型集和逆时偏移图像集；所述逆时偏移图像集为对所述训练速度模型集进行逆时偏移成像得到的图像集；所述训练模块，用于基于所述训练集对预设卷积神经网络进行训练，得到训练之后的卷积神经网络；所述预测模块，用于以所述初始速度模型和初始逆时偏移图像作为输入，利用所述训练之后的卷积神经网络进行预测，得到预测速度模型；所述初始逆时偏移图像为对所述初始速度模型进行逆时偏移成像得到的图像；所述迭代模块，用于判断基于所述预测速度模型正演模拟的地震数据，与所述原始地震数据之间的误差是否小于预设阈值；如果否，则将预测速度模型作为新的初始速度模型，并返回执行所述区块化模块。
10.进一步地，所述区块化模块，还用于：基于k均值聚类，将所述初始速度模型的网格点划分为两个不相交的聚类簇；分别确定每个聚类簇中网格点的最大空间相连区域，得到两个相连区域；基于网格点与所述两个相连区域的曼哈顿距离，将所述初始速度模型区块化为空间相关且连接不相交的两个子模型；对所述两个子模型分别进行递归区块化处理，得到所述特征集；基于所述特征集，创建所述速度模型基。
11.进一步地，所述生成模块，还用于：基于所述速度模型基的每个分量中速度的平均值和方差，创建所述训练速度模型集；对所述训练速度模型集中的每个训练速度模型进行逆时偏移成像，得到所述逆时偏移图像集；将所述训练速度模型集和所述逆时偏移图像集作为所述训练集。
12.第三方面，本发明实施例还提供了一种电子设备，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述第一方面所述的方法的步骤。
13.第四方面，本发明实施例还提供了一种具有处理器可执行的非易失的程序代码的计算机可读介质，所述程序代码使所述处理器执行上述第一方面所述方法。
14.本发明的目的在于提供一种基于卷积神经网络深度学习的反射波形反演方法和
系统，通过应用所提出的区块化方法对上一次迭代中获得的cnn预测速度模型进行分区，以生成随机训练速度模型，从而制备训练速度模型；然后，在对最新生成的随机训练速度模型进行cnn训练后，应用cnn从初始速度模型和相应的初始逆时偏移图像中再次预测未知的真实速度模型。本发明通过上述方式，能够动态调整选择偏置训练速度模型，以从原始起始速度模型和相应的原始偏移图像逐渐预测更准确的速度模型，从而缓解了现有技术中存在的预测的速度模型不准确的技术问题。
附图说明
15.为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
16.图1为本发明实施例提供的一种基于卷积神经网络深度学习的反射波形反演方法的流程图；
17.图2为本发明实施例提供的一种cnn-rwi的工作流的示意图；
18.图3为本发明实施例提供的一种基于卷积神经网络深度学习的反射波形反演系统的示意图。
具体实施方式
19.下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
20.实施例一：
21.图1是根据本发明实施例提供的一种基于卷积神经网络深度学习的反射波形反演方法的流程图。如图1所示，该方法具体包括如下步骤：
22.步骤s102，对原始地震数据进行预处理，得到初始速度模型。
23.可选地，通过对原始地震数据进行反射走时层析或偏移速度分析(migration velocity analysis，mva)处理，得到初始速度模型。
24.区块化步骤s104，通过将初始速度模型区块化为空间相关且连接不相交的特征集，创建速度模型基。
25.步骤s106，基于速度模型基，生成训练集；训练集包括训练速度模型集和逆时偏移图像集；逆时偏移图像集为对训练速度模型集进行逆时偏移成像得到的图像集。
26.步骤s108，基于训练集对预设卷积神经网络进行训练，得到训练之后的卷积神经网络。
27.步骤s110，以初始速度模型和初始逆时偏移图像作为输入，利用训练之后的卷积神经网络进行预测，得到预测速度模型；初始逆时偏移图像为对初始速度模型进行逆时偏移成像得到的图像。
28.步骤s112，判断基于预测速度模型正演模拟的地震数据，与原始地震数据之间的误差是否小于预设阈值；如果否，则将预测速度模型作为新的初始速度模型，并返回执行区
块化步骤s104；如果是，则将预测速度模型作为输出结果。
29.本发明提供了一种基于卷积神经网络深度学习的反射波形反演方法，通过应用所提出的区块化方法对上一次迭代中获得的cnn预测速度模型进行分区，以生成随机训练速度模型，从而制备训练速度模型；然后，在对最新生成的随机训练速度模型进行cnn训练后，应用cnn从初始速度模型和相应的初始逆时偏移图像中再次预测未知的真实速度模型。本发明实施例提供的方法通过上述方式，能够动态调整选择偏置训练速度模型，以从原始起始速度模型和相应的原始偏移图像逐渐预测更准确的速度模型，从而缓解了现有技术中存在的预测的速度模型不准确的技术问题。
30.本发明实施例提供的基于卷积神经网络深度学习的反射波形反演(cnn-rwi)是一种迭代反演方法，用于应用cnn从初始速度模型和相应的初始逆时偏移图像重复预测速度模型。假设原始的初始速度模型是通过反射走时层析或基于偏移速度分析(mva)从观测数据中获得的，相应的原始偏移成像可以通过逆时偏移(reverse time migration，rtm)获得，即：
[0031][0032]
其中us(x,t)和ur(x,t)分别是传播到原始初始速度模型中的源波场和接收器波场，分别位于位置x和时间t的位置。t是最长旅行时间。
[0033]
图2为根据本发明实施例提供的一种cnn-rwi的工作流的示意图。如图2所示，该工作流包含一个外部循环和一个内部循环。其中，外部循环在每次迭代时包含四个主要部分：训练速度模型的准备、cnn训练、cnn预测和误差评估。内部循环是cnn的迭代训练循环，用于迭代更新cnn以拟合训练速度模型。
[0034]
cnn-rwi的一个关键新颖之处在于从先验的初始速度模型中准备训练速度模型。对于cnn-rwi的第一次迭代，最佳先验速度模型是原始初始速度模型，该模型通常通过从观测到的原始地震数据中经过预处理(例如，走时层析或mva)来获得，以估计未知的真实速度模型。对于cnn-rwi的迭代过程，最佳先验速度模型将成为最新的cnn预测速度模型。
[0035]
训练速度模型的准备包括两个主要过程：从先验模型生成模型基(model basis)，以及从模型基生成训练速度模型。通过将先验的速度模型区块化(parcellating)为空间相关和连接不相交的特征集来创建模型基。然后，使用模型基创建训练速度模型。
[0036]
具体的，区块化步骤s104包括如下步骤：
[0037]
步骤s1041，基于k均值聚类，将初始速度模型的网格点划分为两个不相交的聚类簇。
[0038]
步骤s1042，分别确定每个聚类簇中网格点的最大空间相连区域，得到两个相连区域。
[0039]
步骤s1043，基于网格点与两个相连区域的曼哈顿距离，将初始速度模型区块化为空间相关且连接不相交的两个子模型。具体的，该步骤包括：
[0040]
将初始速度模型的网格点划分为三个区域，分别为两个相连区域和第三区域；第三区域为初始速度模型的网格点中除两个相连区域内的网格点之外的所有网格点组成的区域；
[0041]
分别计算第三区域内的每个网格点与两个相连区域的曼哈顿距离；
[0042]
将第三区域内的网格点合并到曼哈顿距离最小的相连区域，得到两个子模型。
[0043]
步骤s1044，对两个子模型分别进行递归区块化处理，得到特征集。
[0044]
步骤s1045，基于特征集，创建速度模型基。
[0045]
本发明实施例提供了一种空间约束、可分割的分层k均值区块化方法(parcellation method)，将速度模型作为模型基，将速度模型区块化为空间相关和连接不相交的特征集，此方法区块化的特征数量取决于图像的渐变。本发明实施例提供的区块化方法将k均值与区域增长相结合，形成分裂(自上而下)的层次结构。
[0046]
最小特征大小的阈值为关键超参数，这是判断是否进一步细分特征的条件，此阈值在权衡每个特征的均匀性和最小尺寸方面起着关键作用。本发明实施例将阈值设置为最大化特征数量，这样就不会创建无意义的小散射区域。本发明实施例提供的区块化步骤的数学解释分为以下6个步骤：
[0047]
步骤1、实施k均值聚类。
[0048]
应用k-均值聚类(k＝2)将先验的初始速度模型v划分为两个不相交的聚类簇c1和c2，其中，v＝c1∪c2，
[0049]
通过最小化下述损失函数获得两个不相交的聚类簇c1和c2：
[0050][0051]
其中，v(p)表示初始速度模型网格点p的速度，μi表示速度的平均值，单位为ci。
[0052]
k-均值聚类方法根据速度分布将初始速度模型v划分为聚类簇c1和c2，而不考虑空间关系和连接。因此，聚类簇c1和c2可能包括许多空间上不相关且互不相连的小区域，应避免在相应的地震数据中产生不需要的散射伪影。
[0053]
步骤2、查找最大的区域。
[0054]
通过独立计算每个聚类簇c1和c2中每个空间不相交区域的网格点数，将最大的空间相连区域分别选择为s1和s2。
[0055]
步骤3、区域增长。
[0056]
初始速度模型v被重新划分为三个不相交的簇v＝{s0,s1,s2}，其中，和而然后，每个模型网格点p∈s0，合并为s1或s2，基于曼哈顿距离dm，使得：
[0057][0058]
曼哈顿距离dm(p,si)，i∈{1,2}定义为(对于2-d模型)：
[0059][0060]
其中，x和y是水平和垂直坐标，s表示群集si中的网格点(对于i∈{1,2})。
[0061]
区域增长步长是通过数值上迭代合并分别与每个聚类簇s1和s2毗邻的模型网格点(p∈s0)实现的。
[0062]
步骤4、将初始速度模型拆分为两个子模型。
[0063]
在区域增长完成后，初始速度模型v最终被区块化成两个空间相关且连接不相交的子模型(s1和s2)：v＝s1∪s2，
[0064]
步骤5和步骤6、递归区块化(parcellate)子模型。
[0065]
对空间相关和连接的不相交子模型s1和s2分别递归包化(区块化)，直到子模型大小(要分区的子模型中的网格点数)小于阈值(允许的最小特征大小)。
[0066]
步骤7、将特征集另存为模型基的组件。
[0067]
如果子模型大小(要分区的子模型中的格网点数)小于阈值(允许的最小特征大小)，则不会递归地对子模型进行区块化分配。相反，它将保存为模型基的组件(即它成为分层树中的叶节点)。
[0068]
通过步骤1-6中的递归区块化，初始速度模型v将被区块化为许多空间相关且相互连接的不相交特征集(子模型)作为模型基：v＝{s1,s2,
…
,sn}，其中，n是特征的数量，并且：如果i≠j。
[0069]
上述步骤1-6的算法描述(空间约束的分裂层次结构k-均值算法)：
[0070]
开始：当模型的大小(网格点数)小于阈值时，将阈值设置为停止条件。
[0071]
运行过程：空间约束的k-均值(模型)。
[0072]
当模型大小≥阈值：
[0073]
①
基于速度模型分布，应用k＝2的k均值方法，将模型划分为两个聚类a和b。
[0074]
②
分别查找聚类a和b中的最大区域x和y。
[0075]
③
让区域x和y竞争以合并模型中的其他区域。
[0076]
④
根据区域x和y将模型拆分为子模型a和b。
[0077]
⑤
递归区块化，得到空间约束k均值(子模型a)。
[0078]
⑥
递归区块化，得到空间约束k均值(子模型b)。
[0079]
结束。
[0080]
⑦
输出：将模型另存为叶节点(模型基的组件)。
[0081]
结束运行过程。
[0082]
具体的，步骤s106还包括如下步骤：
[0083]
步骤s1061，基于速度模型基的每个分量中速度的平均值和方差，创建训练速度模型集。
[0084]
步骤s1062，对训练速度模型集中的每个训练速度模型进行逆时偏移成像，得到逆时偏移图像集。
[0085]
步骤s1063，将训练速度模型集和逆时偏移图像集作为训练集。
[0086]
本发明实施例提供的方法在获得模型基{s1,s2,
…
,sn}后，可以基于模型基的每个分量中速度的平均值e(si)和方差var(si)创建训练速度模型：
[0087][0088]
和
[0089][0090]
其中，ni是每个组件(parcel)si中的网格点数。然后，随机训练速度模型集(m是随机训练速度模型的数量)是通过分配随机值来获得的，该值是从模型基中组件的速度分布中得出的：
[0091][0092]
此处：
[0093][0094]
其中n表示正态分布。
[0095]
为每个特征分配一个随机值，不仅可以增加训练速度模型的多样性，还可以创建具有许多不同同质层的训练速度模型(每个同质层具有均匀速度)。任意两个相邻同质层之间的速度对比在训练速度模型中创建具有不同反射率的反射体，这确保了相应的训练地震数据包含反射事件，就像观测到的数据一样。
[0096]
在前面准备阶段创建训练速度模型后，训练地震数据通过相应的训练速度模型集计算。接着，初始训练速度模型集计算。接着，初始训练速度模型集可以通过一致的预处理(例如，反射走时层析或mva)获得，rtm图像集可以通过逆时偏移(rtm)获得。因此，训练用的初始速度模型集和训练用的rtm图像集是从合成数据获得的，这与从观测数据获得原始初始速度模型v0和rtm图像i0的方式相同。
[0097]
仿照原始处理方式以获得训练初始速度模型(例如，反射走时层析或mva)有利于减少选择偏差，因为训练初始速度模型的获得方式与原始初始速度模型相同。但是，对于所有训练速度模型而言，这种模仿成本很高。本发明实施例中，应用高斯滤波器来模仿昂贵的反射走时层析或mva，通过平滑训练速度模型来获得训练初始速度模型。
[0098]
在训练数据集准备就绪后，训练初始速度模型v和训练rtm图像就形成输入数据，训练cnn以准确再现训练速度模型方法是最小化损失函数。
[0099]
具体的，基于训练集并通过迭代最小化目标损失函数的方式，对预设卷积神经网络进行训练；其中，目标损失函数包括：
[0100][0101]
g是所述预设卷积神经网络，w为所述预设卷积神经网络的参数化的权重，w
*
为所述训练之后的卷积神经网络的参数化的权重，m为所述训练集中训练速度模型的数量，为第k个初始训练速度模型，为第k个训练速度模型，为第k个逆时偏移图像。通过迭代最小化损失函数，cnn可以认为是在训练数据集中精确地逼近这个从和到的映射。
[0102]
在cnn预测中，原始初始速度模型v0和相应的rtm图像i0的组合被输入到训练得到的cnn来预测未知的真实速度模型：
[0103]vp
＝g(v0,i0；w
*
)
[0104]
其中v
p
是cnn的预测速度模型。
[0105]
cnn预测速度模型v
p
后，需要进行正演建模以获得相应的模拟地震数据来进行误差评估(例如，计算最小二乘数据残差)。与通过最小化每次迭代时的数据误差来更新速度模型的传统fwi不同，cnn-rwi计算数据误差用以间接评估cnn预测的速度模型v
p
的精度。如果从cnn预测的速度模型v
p
获得的数据误差不符合收敛条件(例如，最小数据误差的阈值)，则cnn预测的速度模型v
p
将成为最新生成的先验模型，并进行打包以生成训练速度模型集，以供下一次迭代使用。
[0106]
可选地，本发明实施例提供的方法中，cnn-rwi的另一个关键新颖之处在于，在外循环中的每次迭代中，基于最新的先验速度模型对训练数据集进行动态调整。最新的先验速度模型定义为第一次迭代时的原始初始速度模型v0，或在每次迭代中获得的cnn的预测速度模型v
p
。
[0107]
由于原始的初始速度模型包含不准确的先验信息内容(例如，不准确的背景速度)，因此在第一次迭代时，相应的训练速度模型集是有选择倾向的，这是因为是原始初始速度模型v0的良好表示，而非未知的真实速度模型v
t
的表示。相反，在理想情况下，从训练速度模型集获取训练用的初始速度模型集和训练用的偏移图像集i的预处理过程可与从未知真实速度模型v
t
获取原始初始速度模型v0和偏移图像i0的预处理过程高度一致。换言之，空间关系不像训练数据集中的信息内容那样有选择倾向，因为前者与实际应用程序共享类似的预处理，但后者与实际数据集{v0,i0,v
t
}中的信息内容不同。因此，通过最小化训练数据集中的损失函数以精确地逼近从和到的映射，cnn就预测了未知速度模型v
p
，该模型通常比原始初始速度模型v0更准确，但仍与未知的真实速度模型v
t
相去甚远。然后，cnn-rwi动态调整训练数据集方法是将上一次迭代中获得的最新cnn预测速度模型v
p
作为最新的先验速度模型，为cnn训练创建新的模型基和相应的训练速度模型集
[0108]
cnn-rwi的收敛依赖于这样的假设：cnn的预测速度模型v
p
比训练速度模型集更准确，后者表示的是在上一次迭代(或第一次迭代时的原始初始速度模型)获得的cnn的预测速度模型。该假设是合理的，因为原始速度模型v0和原始偏移图像i0输入到cnn以预测速度模型v
p
，在每次迭代时，都是从观测到的数据获得的，而这些数据是从未知的真实速度模型v
t
获得的。相比之下，cnn训练用的cnn输入数据(v和)是从合成数据中获得的，这些数据是从训练速度模型集计算得出的。因此，cnn的预测速度模型v
p
将比训练速度模型集更接近未知的真实速度模型v
t
。
[0109]
因此，在每次迭代时，这种训练速度模型的动态调整隐式地减少了训练速度模型集的选择倾向，从而使cnn更能代表未知的真实速度模型。从而cnn通过隐式减少选择倾向，而不是像fwi那样最小化数据误差，从而更准确地逐步预测未知的真实速度模型v
p
。
[0110]
由以上描述可知，本发明实施例提供了一种基于卷积神经网络深度学习的反射波形反演方法，与常规的fwi、rwi、以及普通的基于cnn的用于预测速度模型的技术相比，本发明实施例提出的cnn-rwi具有以下优势：
[0111]
首先，对于深层速度模型的反演，fwi和rwi算法主要是分别更新界面和背景速度，而cnn-rwi可以有效地更新二者。
[0112]
其次，fwi和rwi都依赖于低频数据，使得旅行时差局限于半个周期内，以避免在数据拟合过程中周期跳跃；而cnn-rwi更倾向于高频数据以获得更高分辨率的偏移图像作为cnn输入数据，因为cnn-rwi不需要数据拟合程序。
[0113]
第三，用于预测速度模型的常规cnn技术通常无法表示目标区未知的真实速度模型(因此表现出选择偏差)，相比之下，基于cnn的rwi结合了初始速度模型和偏移图像中的先验信息，以准确预测所有深度的速度模型。
[0114]
实施例二：
[0115]
图3为根据本发明实施例提供的一种基于卷积神经网络深度学习的反射波形反演系统的示意图。如图3所示，该系统包括：预处理模块10，区块化模块20，生成模块30，训练模块40，预测模块50和迭代模块60。
[0116]
具体的，预处理模块10，用于对原始地震数据进行预处理，得到初始速度模型。可选地，通过对原始地震数据进行反射走时层析或偏移速度分析(migration velocity analysis，mva)处理，得到初始速度模型。
[0117]
区块化模块20，用于通过将初始速度模型区块化为空间相关且连接不相交的特征集，创建速度模型基。
[0118]
生成模块30，用于基于速度模型基，生成训练集；训练集包括训练速度模型集和逆时偏移图像集；逆时偏移图像集为对训练速度模型集进行逆时偏移成像得到的图像集。
[0119]
训练模块40，用于基于训练集对预设卷积神经网络进行训练，得到训练之后的卷积神经网络。
[0120]
预测模块50，用于以初始速度模型和初始逆时偏移图像作为输入，利用训练之后的卷积神经网络进行预测，得到预测速度模型；初始逆时偏移图像为对初始速度模型进行逆时偏移成像得到的图像。
[0121]
迭代模块60，用于判断基于预测速度模型正演模拟的地震数据，与原始地震数据之间的误差是否小于预设阈值；如果否，则将预测速度模型作为新的初始速度模型，并返回执行区块化模块20。
[0122]
本发明提供了一种基于卷积神经网络深度学习的反射波形反演系统，通过应用所提出的区块化方法对上一次迭代中获得的cnn预测速度模型进行分区，以生成随机训练速度模型，从而制备训练速度模型；然后，在对最新生成的随机训练速度模型进行cnn训练后，应用cnn从初始速度模型和相应的初始逆时偏移图像中再次预测未知的真实速度模型。本发明通过上述方式，能够动态调整选择偏置训练速度模型，以从原始起始速度模型和相应的原始偏移图像逐渐预测更准确的速度模型，从而缓解了现有技术中存在的预测的速度模型不准确的技术问题。
[0123]
可选地，区块化模块20，还用于：
[0124]
基于k均值聚类，将初始速度模型的网格点划分为两个不相交的聚类簇；
[0125]
分别确定每个聚类簇中网格点的最大空间相连区域，得到两个相连区域；
[0126]
基于网格点与两个相连区域的曼哈顿距离，将初始速度模型区块化为空间相关且连接不相交的两个子模型；具体的，包括：将所述初始速度模型的网格点划分为三个区域，分别为所述两个相连区域和第三区域；所述第三区域为所述初始速度模型的网格点中除所述两个相连区域内的网格点之外的所有网格点组成的区域；分别计算所述第三区域内的每个网格点与所述两个相连区域的曼哈顿距离；将所述第三区域内的网格点合并到曼哈顿距离最小的相连区域，得到所述两个子模型；
[0127]
对两个子模型分别进行递归区块化处理，得到特征集；
[0128]
基于特征集，创建速度模型基。
[0129]
具体的，生成模块30，还用于：
[0130]
基于速度模型基的每个分量中速度的平均值和方差，创建训练速度模型集；
[0131]
对训练速度模型集中的每个训练速度模型进行逆时偏移成像，得到逆时偏移图像集；
[0132]
将训练速度模型集和逆时偏移图像集作为训练集。
[0133]
具体的，训练模块40还用于：
[0134]
基于训练集并通过迭代最小化目标损失函数的方式，对预设卷积神经网络进行训练；其中，目标损失函数包括：
[0135][0136]
g是预设卷积神经网络，w为预设卷积神经网络的参数化的权重，w
*
为训练之后的卷积神经网络的参数化的权重，m为训练集中训练速度模型的数量，为第k个初始训练速度模型，为第k个训练速度模型，为第k个逆时偏移图像。
[0137]
本发明实施例还提供了一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现上述实施例一中的方法的步骤。
[0138]
本发明实施例还提供了一种具有处理器可执行的非易失的程序代码的计算机可读介质，程序代码使处理器执行上述实施例一中的方法。
[0139]
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
技术特征：1.一种基于卷积神经网络深度学习的反射波形反演方法，其特征在于，包括：对原始地震数据进行预处理，得到初始速度模型；区块化步骤：通过将所述初始速度模型区块化为空间相关且连接不相交的特征集，创建速度模型基；基于所述速度模型基，生成训练集；所述训练集包括训练速度模型集和逆时偏移图像集；所述逆时偏移图像集为对所述训练速度模型集进行逆时偏移成像得到的图像集；基于所述训练集对预设卷积神经网络进行训练，得到训练之后的卷积神经网络；以所述初始速度模型和初始逆时偏移图像作为输入，利用所述训练之后的卷积神经网络进行预测，得到预测速度模型；所述初始逆时偏移图像为对所述初始速度模型进行逆时偏移成像得到的图像；判断基于所述预测速度模型正演模拟的地震数据，与所述原始地震数据之间的误差是否小于预设阈值；如果否，则将预测速度模型作为新的初始速度模型，并返回执行所述区块化步骤。2.根据权利要求1所述的方法，其特征在于，所述区块化步骤包括：基于k均值聚类，将所述初始速度模型的网格点划分为两个不相交的聚类簇；分别确定每个聚类簇中网格点的最大空间相连区域，得到两个相连区域；基于网格点与所述两个相连区域的曼哈顿距离，将所述初始速度模型区块化为空间相关且连接不相交的两个子模型；对所述两个子模型分别进行递归区块化处理，得到所述特征集；基于所述特征集，创建所述速度模型基。3.根据权利要求2所述的方法，其特征在于，基于网格点与所述两个相连区域的曼哈顿距离，将所述初始速度模型区块化为空间相关且连接不相交的两个子模型，包括：将所述初始速度模型的网格点划分为三个区域，分别为所述两个相连区域和第三区域；所述第三区域为所述初始速度模型的网格点中除所述两个相连区域内的网格点之外的所有网格点组成的区域；分别计算所述第三区域内的每个网格点与所述两个相连区域的曼哈顿距离；将所述第三区域内的网格点合并到曼哈顿距离最小的相连区域，得到所述两个子模型。4.根据权利要求1所述的方法，其特征在于，基于所述速度模型基，生成训练集，包括：基于所述速度模型基的每个分量中速度的平均值和方差，创建所述训练速度模型集；对所述训练速度模型集中的每个训练速度模型进行逆时偏移成像，得到所述逆时偏移图像集；将所述训练速度模型集和所述逆时偏移图像集作为所述训练集。5.根据权利要求1所述的方法，其特征在于，基于所述训练集对预设卷积神经网络进行训练，包括：基于所述训练集并通过迭代最小化目标损失函数的方式，对所述预设卷积神经网络进行训练；其中，所述目标损失函数包括：
g是所述预设卷积神经网络，w为所述预设卷积神经网络的参数化的权重，w
*
为所述训练之后的卷积神经网络的参数化的权重，m为所述训练集中训练速度模型的数量，为第k个初始训练速度模型，为第k个训练速度模型，为第k个逆时偏移图像。6.一种基于卷积神经网络深度学习的反射波形反演系统，其特征在于，包括：预处理模块，区块化模块，生成模块，训练模块，预测模块和迭代模块；其中，所述预处理模块，用于对原始地震数据进行预处理，得到初始速度模型；所述区块化模块，用于通过将所述初始速度模型区块化为空间相关且连接不相交的特征集，创建速度模型基；所述生成模块，用于基于所述速度模型基，生成训练集；所述训练集包括训练速度模型集和逆时偏移图像集；所述逆时偏移图像集为对所述训练速度模型集进行逆时偏移成像得到的图像集；所述训练模块，用于基于所述训练集对预设卷积神经网络进行训练，得到训练之后的卷积神经网络；所述预测模块，用于以所述初始速度模型和初始逆时偏移图像作为输入，利用所述训练之后的卷积神经网络进行预测，得到预测速度模型；所述初始逆时偏移图像为对所述初始速度模型进行逆时偏移成像得到的图像；所述迭代模块，用于判断基于所述预测速度模型正演模拟的地震数据，与所述原始地震数据之间的误差是否小于预设阈值；如果否，则将预测速度模型作为新的初始速度模型，并返回执行所述区块化模块。7.根据权利要求6所述的系统，其特征在于，所述区块化模块，还用于：基于k均值聚类，将所述初始速度模型的网格点划分为两个不相交的聚类簇；分别确定每个聚类簇中网格点的最大空间相连区域，得到两个相连区域；基于网格点与所述两个相连区域的曼哈顿距离，将所述初始速度模型区块化为空间相关且连接不相交的两个子模型；对所述两个子模型分别进行递归区块化处理，得到所述特征集；基于所述特征集，创建所述速度模型基。8.根据权利要求6所述的系统，其特征在于，所述生成模块，还用于：基于所述速度模型基的每个分量中速度的平均值和方差，创建所述训练速度模型集；对所述训练速度模型集中的每个训练速度模型进行逆时偏移成像，得到所述逆时偏移图像集；将所述训练速度模型集和所述逆时偏移图像集作为所述训练集。9.一种电子设备，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现上述权利要求1至5任一项所述的方法的步骤。10.一种具有处理器可执行的非易失的程序代码的计算机可读介质，其特征在于，所述
程序代码使所述处理器执行所述权利要求1-5任一项所述方法。
技术总结本发明提供了一种基于卷积神经网络深度学习的反射波形反演方法和系统，包括：对原始地震数据进行预处理，得到初始速度模型；区块化步骤：通过将初始速度模型区块化为空间相关且连接不相交的特征集，创建速度模型基；基于速度模型基，生成训练集；基于训练集对预设卷积神经网络进行训练，得到训练之后的卷积神经网络；以初始速度模型和初始逆时偏移图像作为输入，利用训练之后的卷积神经网络进行预测，得到预测速度模型；判断基于预测速度模型正演模拟的地震数据，与原始地震数据之间的误差是否小于预设阈值；如果否，则将预测速度模型作为新的初始速度模型，并返回执行区块化步骤。本发明缓解了现有技术中存在的预测的速度模型不准确的技术问题。型不准确的技术问题。型不准确的技术问题。
技术研发人员：王彦飞 吴毓朗
受保护的技术使用者：中国科学院地质与地球物理研究所
技术研发日：2022.04.06
技术公布日：2022/7/5
