基于多任务的断层智能识别方法

1.本发明属于地质勘探与解释领域，特别涉及一种断层的识别技术。


背景技术：

2.断层的识别解释是地质勘探与解释中重要的组成部分。长期以来，断层解释通过应用图像分析，图像增强的方法对相应的断层属性进行增强。按照这种思路发展而来的断层解释方法，有着计算步骤冗长，流程复杂，参数众多难以调整，导致断层解释的自动化，智能化程度不高的缺点。随着深度学习的兴起，其凭借着强大的高维复杂特征提取能力，以及端到端智能化的特点被快速引入到断层识别的领域中，并取得了不错的成绩。
3.在断层解释领域，最开始是由地质专家对地震数据中的断层结构进行人工解释。人工解释将三维地震数据进行二维切片，依靠专家相关地学经验对地震数据中反射层的不连续特征进行分析，最后得到断层解释的结果。使用人工解释的方法大量消耗人力资源，费时费力，且断层解释结果依赖于解释人员的主观判断。基于人工进行断层识别在断层交错等复杂结构下难以实现准确的识别。因此人工识别的方法越来越难以适应规模越来越大的断层识别任务的需求，学术界也在一直探究使用自动方法高效的对断层进行识别。其中相干体属性，方差体属性，基于地震数据的倾角体和方位角体属性，以及利用蚂蚁追踪得到的蚂蚁体属性等多种的属性体增强数据能够有效的对断层进行识别。这些技术中以相干体数据与蚂蚁体数据影响较大。随着卷积神经网络的提出，深度学习网络被广泛应用于计算机视觉领域，并在自然图像，医学图像的分类，目标检测，语义分割，实例分割等领域取得较大突破。各个学科与图像处理相关交叉领域中也开始应用深度学习技术，并取得巨大成功。断层识别依靠识别地震图像中由于断层构造所产生的不连续信息，非常适合深度学习中计算机视觉相关技术，对断层进行识别类似于深度学习中的语义分割任务。
4.k.m.tingdahl,m.de rooij等在2005年提出了一种基于多属性数据的断层半自动算法。该网络将人为挑选的属性作为输入输入到多层感知机中，最终输出断层识别结果。该算法需要人工选择属性无法做到自动的断层识别，并且断层的识别结果依赖于人为经验。2017年，l.huang,x.dong将理论正演数据作为训练数据，使用卷积神经网络对分块后的多种断层属性数据进行逻辑回归运算，得到断层识别结果。2017年w.xiong等在地震属性数据中选取一个待分类点，并以该点为中心inline，xline，timeline三个方向进行切片，将切片数据作为神经网络的输入数据利用分类网络对中心点是否为断层进行分类。该方法使用相干体数据作为输入，在实际工区数据中也取得了一定的效果。上述的两种方法都需要将原始数据进行分块或切片，且每次预测仅预测原始数据中一个点的类别，会有大量的重复计算以及存储的额外开销，时间成本较高。2019年，x.m.wu将3d u-net应用于断层识别任务并对原始的网络进行简化以及改进提出了faultseg3d模型，该模型充分利用了地震数据中的三维空间信息，在识别的效果以及精度上达到了较高水平，并在实际的工区数据中得到了较好的效果，该网络以一定大小的三维地震数据作为输入，输出断层的语义分割结果。该算法使用原始的地震数据作为输入，做到了对断层数据端到端的识别。
5.虽然目前基于faultseg3d的网络取得了较好的效果但是其仍然存在一定的问题：首先地震数据随着工区的不同，以及地质结构的不同具有较大的差异，在分布上具有较大的不同，而现阶段使用的训练数据大多数是基于正演模型的理论数据，使用faultseg3d容易造成模型的过拟合，降低模型的泛化性。在实际不同工区中识别结果表现出较大差异。


技术实现要素：

6.为解决上述技术问题，本发明提出一种基于多任务的断层智能识别方法，过引入倾角方位角检测这两个辅助任务，利用主任务以及辅助任务任务间的相关性以及公共特征的提取，提高了模型的泛化性，使得模型在不同分布，不同特征的地震数据中都有较好的表现。
7.本发明采用的技术方案为：基于多任务的断层智能识别方法，包括：
8.s1、构建理论数据集：通过正演模型构建理论数据集，所述理论数据集中包括若干组理论数据，每一组理论数据中包括四个数据样本，具体为：地震振幅数据、断层标签数据、倾角体数据、方位角体数据；
9.s2、对理论数据集中的理论数据进行增强处理；
10.s3、构建网络模型：对3d u-net网络进行改进，将其扩展为多任务网络，得到多任务的u-net网络；
11.s4、对网络进行训练：将步骤s2增强后的理论数据输入多任务的u-net网络中进行训练；
12.s5、断层解释阶段：将实际工区数据经过数据预处理之后送入到训练好的多任务的u-net网络中进行预测，得到断层属性体、断层倾角体、断层方位角体。
13.步骤s2所述的增强处理包括：
14.s21、对三维数据体进行旋转变换；
15.s22、采取对应角度间隔分别对倾角体数据、方位角体数据进行分类取整。
16.步骤s3所述多任务的u-net网络结构为：包括编码器部分与解码器部分；所述编码器部分包括四个编码器层，所述解码器部分包括四个解码器层；
17.编码部分相邻的2个编码器层，通过个2
×2×
2大小、步长为2的池化层将上一层编码器层的输出进行下采样后输入下一层的编码器层，下一层的编码器层中卷积层的通道数是上一层编码器层中卷积层通道数的2倍；
18.解码部分相邻的2个解码器层，通过一个2
×2×
2的上采样层将上一层解码器层的输出进行上采样后输入下一层的解码器层，下一层的解码器层中卷积层的通道数是上一层解码器层中卷积层通道数的一半；
19.每个编码器层包括：两个3
×3×
3的卷积层，每个卷积层使用relu作为激活函数；
20.每个解码器层包括：一个与对应编码层相连的拼接层，以及两个3
×3×
3的卷积层，每个卷积层均使用relu作为激活函数。
21.在断层解释任务中将第四个解码器层输出的特征图进行两次残差模块的卷积后将特征图的通道数量减少到16个，得到用于断层识别任务的特征。
22.在倾角检测任务中，将第三个解码器层输出的特征图和第四个解码器层输出的特征图进行拼接，再将拼接后的特征图进行三次残差卷积块的卷积后得到最终的倾角检测结
果。
23.在方位角体检测任务中，将第三个解码器层输出的特征图和第四个解码器层输出的特征图进行拼接，再将拼接后的特征图进行三次残差卷积块的卷积后得到最终的方位角体检测结果。
24.本发明的有益效果：本发明基于三维地震振幅，提出了一种基于3d u-net多任务网络的断层识别方法。该网络以三维的理论地震数据作为输入，以断层识别，倾角检测，方位角检测作为网络输出的任务。该算法能够对断层进行准确的分割，并有效的检测断层所对应的倾角和方位角。在该算法下，引入了基于滑动窗口的断层预测拼接方法对传统的断层预测进行优化，解决了断层预测中不连续过渡不自然的缺点。通过多任务学习的方式提高了模型的鲁棒性，使其在不同特征分布的实际工区数据的表现上有较大的进步。
附图说明
25.图1为本发明的方法流程图；
26.图2为地震数据及相关属性示意图；
27.其中，(a)为地震振幅数据，(b)为断层标签数据，(c)为方位角体，(d)为倾角体；
28.图3为基于多任务的u-net网络示意图；
29.图4为基于滑动窗口的断层预测流程图；
30.图5为基于多任务的断层识别网络分割结果；
31.其中，(a)为原始地震振幅数据，(b)为断层标签数据，(c)为倾角体标签数据，(d)为倾角体预测数据，(e)为方位角标签数据，(f)为方位角预测数据；
32.图6为断层建模结果；
33.图7为断层预测结果对比；
34.其中，(a)为原始地震数据，(b)为断层标签数据，(c)为3d u-net预测结果，(d)为多任务网络预测结果；
35.图8为在不同频率下的预测结果对比；
36.其中，(a)为地震振幅数据，(b)为断层标签数据，(c)为3d u-net预测数据，(d)为多任务网络预测数据；
37.图9为f3数据断层分割对比；
38.其中，(a)为原始地震数据，(b)为3d u-net预测结果，(c)为多任务网络识别结果；
39.图10为西南某地区地震数据断层分割对比；
40.其中，(a)为原始地震数据，(b)为3d u-net预测结果，(c)为多任务网络识别结果；
41.图11为三重阈值预测示意图；
42.图12为地震振幅数据及信噪比测试结果；
43.其中，(a)为地震振幅数据，(b)为信噪比测试结果；
44.图13为断层分割结果对比图；
45.其中，(a)为直接预测结果，(b)为三重阈值过滤结果。
具体实施方式
46.为便于本领域技术人员理解本发明的技术内容，下面结合附图对本发明内容进一
步阐释。
47.实施例1
48.本发明使用基于3d u-net网络的多任务学习网络同时对断层属性体以及倾角体和方位角体进行预测。针对使用传统的扫描法对断层倾角方位角体的计算时间复杂度大，计算时间长且依赖于断层属性体提取结果的缺点，本发明提出了一种基于3d u-net网络的多任务学习网络。该网络以三维地震振幅数据作为输入，输出断层识别结果以及倾角体和方位角体。针对训练的地震数据不同工区之间的数据分布差异，使用多任务学习网络，提取不同相关任务中的公共特征，减少模型过拟合的风险，提高模型在不同工区数据之间的泛化性。
49.如图1所示，本发明的方法包括以下步骤：
50.1.构建理论数据集：通过正演模型构建理论数据集的模型，共包含四个部分，地震振幅数据，断层标签数据，理论倾角体数据，理论方位角体数据。
51.2.构建网络模型：对3d u-net网络进行改进，并将其扩展为多任务网络。
52.3.对网络进行训练：对理论数据进行数据增强之后送入到多任务学习网络中进行训练。
53.4.断层解释阶段：将实际工区数据送入到多任务学习网络中进行预测，得到断层属性体，断层倾角体，断层方位角体，三个属性体。
54.步骤1还包括对理论数据集进行处理，在断层智能识别任务中使用人工标注的方法受标注者自身主观影响较大，很难对数据中的断层标签进行准确的标注。在真实地震振幅数据中求解倾角体和方位角工作量大，计算精度较差，无法获得足够且准确的训练数据对网络进行训练。因此使用理论地震数据是一种更合理的选择。本发明使用正演模型构建了220个128
×
128
×
128的理论地震数据。
55.使用其中的前200个数据作为训练数据，后20个数据作为测试数据。每一组训练数据中包含四个数据样本，地震振幅数据，断层标签数据，倾角体数据，方位角体数据。在本发明中断层的倾角范围为[30
°
,90
°
]，方位角范围为[0
°
,360
°
]，其训练样本如图2所示。
[0056]
数据集具体处理过程为：
[0057]
(1)数据旋转处理
[0058]
对一个三维数据体进行旋转变换即选定一个旋转点，并将数据体中所有的样本点围绕旋转点旋转一个角度。本步骤应用于地震振幅数据、断层标签数据、倾角体以及方位角体数据的数据旋转处理。
[0059]
其中c(x0，y0)是旋转点，θ是旋转角，样本数据中的任意点为p1(x，y)，其旋转后的数据为p2(x
′
，y
′
)，x0、y0是旋转点的横纵坐标，x、y是样本数据的横纵坐标，其中x是地震振幅数据，y是断层标签数据，倾角体旋转公式为：
[0060]
p2＝m(p
1-c) c
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0061]
其中，m为旋转矩阵，旋转矩阵如下式所示：
[0062][0063]
由于在本实施例中使用的数据集均为正方体数据，为了避免旋转任意角度所带来的插值问题，以及数据空缺问题，本实施例仅仅使用90
°
，180
°
，270
°
，三个角度对数据体进
行旋转。并且由于地震数据在x轴与y轴进行旋转之后其分布不具有实际的物理意义，所以旋转的方向仅限于z轴方向。
[0064]
在沿z轴方向旋转之后倾角方向没有变化，方位角会相应的增加旋转的角度。方位角体经过旋转变换后的表达式如下所示：
[0065]
p2＝m(p
1-c) c θ
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)
[0066]
(2)对倾角体和方位角体数据的处理
[0067]
对三维地震数据进行断层识别可以看作是一个分割问题，同时也可以看作是在原始数据大小的结果图中对每一个像素点的分类问题。而对倾角体，方位角体的检测问题可以看作是在原始数据大小的结果图上的回归问题。使用回归问题对断层的倾角方位角体进行训练波动大，并且由于回归问题与分类问题相似度不高，使用同一主干网络(backbone)对数据进行公共特征提取会导致后续不同任务间相互干扰，导致网络训练不收敛，无法达到预期的效果。因此本发明将倾角体与方位角体的识别转换为分类问题。
[0068]
在本发明构建的地震数据中断层的倾角范围为[30
°
,90
°
]，本实施例采取2
°
为间隔，对断层倾角数据进行分类取整，最终将断层的倾角体分为32类(无断层数据的为一类)。同理，本发明所构建地震数据中方位角的范围是[0
°
,360
°
]，本实施例采取5
°
为间隔将其进行分类取整，将方位角体分为74类(无断层数据的单独为一类)。
[0069]
将原有的倾角体方位角体的有回归任务转换为分类任务，增加了不同任务之间的相似性，有利于公共特征的提取同时加速了模型的训练和收敛。
[0070]
步骤2所述的基于多任务的u-net网络如图3所示，该网络基于u-net网络演化而来，包括两个部分，第一个部分是网络的编码器部分，网络输入为128
×
128
×
128的三维地震数据，编码器的第一层包括两个3
×3×
3的卷积层，每个卷积层均使用relu作为激活函数，然后采用2
×2×
2大小，步长为2的池化层对网络进行下采样。在编码器部分每向下一层就会将该层卷积层的通道数翻倍。该网络一共包括四层这样的编码器层。另一个部分为解码器部分，每个解码层都包括一个2
×2×
2的上采样层，一个与左边的编码层相连的拼接层，以及两个3
×3×
3的卷积层，每个卷积层均使用relu作为激活函数。在解码器部分每向下一层就会将该层卷积层的通道数减少一半。该网络一共包含四层这样的解码器层。
[0071]
本发明采取这样简化的3d u-net网络作为断层识别的主干网络对多个任务间公共特征进行提取。使用这样的3d u-net网络在断层识别任务中取得过良好的效果。在特征提取中通过对上采样和下采样的使用可以使得特征提取器能够提取到不同尺度上的断层图像特征。在断层解释任务中将骨干网络的特征图进行两次残差模块的卷积后将特征图的通道数量减少到16个，得到用于断层识别任务的特征。在倾角检测任务中，将骨干网络的特征图和断层识别的特征图进行拼接，在将拼接后的特征图进行三次残差卷积块的卷积后得到最终的倾角检测结果。同理使用类似的网络结构得到方位角体的检测结果。
[0072]
相较于倾角体和方位角体的检测任务，断层的识别任务是相对更简单的任务。因此在对断层进行识别的时候使用了相对较少的卷积层数。而对于倾角体和方位角体的识别使用相对较多的卷积网络。为了避免过深的网络导致网络退化，本发明使用残差块代替传统卷积网络缓解网络过深带来的模型退化问题。在传统的倾角方位角检测算法中，为了提高检测效果，往往在断层属性体上进行计算。因此在本发明的断层网络结构构建中将断层识别的特征作为输入提供到倾角与方位角的检测任务中。断层识别特征与骨干网络的公共
特征共同作为输入，为倾角与方位角的检测任务提供了不同尺度下的特征以及部分先验信息，减小了后续网络学习的难度。
[0073]
步骤3训练过程中所涉及的损失函数具体为：
[0074]
断层识别任务
[0075]
三维的断层识别任务与二维断层识别任务具有相似的特点。在损失函数的选择上需要考虑到断层与非断层间类别不均衡问题以及从整体图像相似性上描述断层识别结果。因此在本发明网络训练中选取bce dice损失函数。
[0076]
方位角体的检测任务
[0077]
断层的方位角检测问题转换为了对像素点的分类问题。由于断层在方位角上的分布较为均匀，没有明显的分布不均匀的情况，并且在数据增强中使用了旋转的增强方式进一步缓解了类别不均衡问题，因此在方位角检测任务中使用交叉熵函数作为损失函数，其数学表达式如下所示：
[0078][0079]
其中s示方位角数据总类别数，y
ic
表示样本i的类别是否是类别c，p
ic
表示方位角体样本i属于类别c的概率值，n代表类别总数。
[0080]
倾角体检测任务
[0081]
同方位角检测任务不同，在理论数据中倾角范围被限制在了30
°
～90
°
。其中大部分断层位于60
°
～90
°
，有明显的类别不均衡问题。且不同于断层分割这类的二分类问题，方位角检测问题中一共有32个类别，较多的分类类别使得平衡交叉熵损失函数效果欠佳，因此本发明使用focal loss作为倾角检测的损失函数，其数学表达式如下：
[0082][0083]
其中α代表该类别的权重，γ代表该样本分类的难易程度，y代表标签的类别。
[0084]
总体上来讲，focal loss能够有效的缓解分类问题中类别问题不均衡，不同样本间难易程度不均衡的情况。
[0085]
实施例2
[0086]
针对传统预测方法在对实际数据进行预测时，模型的边界会出现不连续，过渡不自然的情况，本发明中提出了一种基于滑动窗口的断层预测方法，最终有效提高了断层识别的连续性。
[0087]
具体包括：
[0088]
1、测试时增强
[0089]
对预测数据进行测试时增强(tta)。本发明选取的测试时增强算法如下：
[0090]
11.将数据在水平方向上随机旋转90
°
，180
°
，270
°
；
[0091]
12.将原始地震振幅数据分别在xline,inline,timeline方向上随机平移5像素以内[-5,5]。
[0092]
2、基于滑动窗口的断层预测方法
[0093]
由于u-net网络的在进行卷积，下采样以及上采样操作时图像边缘信息相较于中心位置信息会有较大程度上的损失，因此u-net网络的预测结果边缘的效果相较于中心的
效果较差。使用传统方法将地震数据进行简单分块后再拼接容易在拼接边缘产生较多的噪声，且不同分块之间的断层识别不连续，断层属性数据间的过渡不自然的缺点。
[0094]
本发明采用如图4所示的基于滑动窗口的方法，将窗口的大小设置为基本单位的大小，同时将滑动窗口的步长设置为窗口大小的一半，在对地震体数据进行预测时不断移动窗口，并对窗口内的地震数据进行预测，并将预测结果最中心的1/2部分保留赋权值为1，边缘部分线性降低其权值，将每个点的值累加，待整个数据体完成预测之后除以每个点的权值得到最后的预测结果。
[0095]
实施例3
[0096]
在使用基于滑动窗口的断层智能预测方法后位于预测图图像居中位置的断层分割结果具有了良好的连续性和自然的过渡，但是在分割结果的最外层效果依然不理想。本发明在对地震数据进行预测时首先需要将预测数据的外层宽度为模型宽度1/2的数据进行镜像的扩展。在将扩展之后的数据作为预测数据进行预测，最后将扩展后的地震数据进行切割，获得原始地震数据的断层识别结果。
[0097]
实施例4
[0098]
以下结合理论数据对本发明的效果进行说明
[0099]
本发明将综合基于3d u-net的多任务网络在理论数据上的断层预测效果和其在实际工区数据上的断层识别效果进行多方面的验证。为了验证本发明提出的基于3d u-net的多任务网络在断层分割任务中的对模型泛化性能的改进，本发明选取两组不同特征，不同频率的数据集对网络对断层分割结果进行测试。其中一组数据集的频率与训练数据一致，而另一组数据在地震波频率，以及断层的分布规律上与训练数据不同。
[0100]
首先分析模型在三维理论数据的效果。
[0101]
首先选择一组理论数据，该理论数据与训练数据的频率以及数据分布特征相同，如图5为其对应的辅助任务，倾角体检测结果，以及方位角检测结果。
[0102]
理论地震数据的网络输出结果如图5所示。图中(a)表示原始理论地震振幅数据；(b)是理论地震振幅数据对应的断层标签数据；(c)是倾角体的标签数据。(d)是倾角体预测数据。(e)方位角的标签数据。(f)方位角的预测数据。
[0103]
从该数据可以看出使用基于多任务的断层分割网络能够有效的提取出断层相关任务的公共特征。利用该公共特征，后续的倾角体检测网络以及方位角体的检测网络能够对该地震数据的倾角体以及方位角体做出有效的预测。从该理论数据来看本发明提出的基于多任务的断层分割网络在断层的倾角识别以及方位角识别上能够做到较为准确的识别，识别效果能够有效的表征断层的走向以及趋势。
[0104]
倾角方位角属性以及断层的预测结果是地质建模中重要的属性数据。利用这些属性数据能够完成断层的细化，追踪，重建等地质建模工作。网络中的倾角，以及方位角解释结果能否用于地质建模并得出正确的建模结果是判断网络识别效果的重要标准之一。因此，此处选取上述理论数据标签结果对数据进行断层追踪建模，其三维断层建模结果如图6所示，可以看出根据网络输出的断层建模结果能够准确追踪出不同的断层以及断层间的位置关系。
[0105]
多任务学习中更重要的是利用辅助任务和主任务之间的相关性提取公共特征进一步改善主任务效果。在本发明网络中倾角方位角检测任务的公共特征有利于促进断层分
u-net网络在横向断层的分割结果中断层的连续性不足，在xline[130,140]，inline[100,250]范围内未能够将断层进行有效分割，出现了遗漏的现象。同时在断层间的交叉位置，以及断层数据的边缘位置上断层分割结果的连续性较差。而如图9中(c)所示的本发明提出的多任务网络的识别结果中，则能够对断层进行准确的识别以及分割，在xline-inline方向上以及inline-timeline方向上断层分割的连续性都优于3d u-net网络的识别结果。
[0116]
如图10所示为西南某地区的原始地震数据以及不同模型的断层识别结果对比。西南某地区工区实际数据大小为521
×
581
×
1000，具有良好的断层发育以及复杂的地质构造，其中包含了正断层，逆断层，溶洞以及褶皱等地质构造，其中溶洞与断层具有相似的特征对断层的识别造成较大的难度。在该地震数据中，地质结构较为复杂产生较大的相干噪声对断层的识别也会造成较大的影响。且该工区数据的数据频率与网络训练数据集中的数据频率不同对断层的预测的影响也较大。
[0117]
从断层预测的结果上来看基于多任务网络的断层识别结果明显优于3d u-net，3du-net网络丢失了大量的断层识别结果，仅仅完整的识别出了一条主要断层。而本发明提出的网络则能够对断层进行准确的识别，保留了较多的断层细节。3d u-net网络的断层识别结果边缘较为模糊，断层识别宽度较大，在大跨度断层的连续性上表现不佳，总体上较为粗造。本发明提出的网络分割结果较为细致，边缘相对清晰，连续性较好。在timeline[600,800]的范围内，地质结构较为复杂，存在较多细碎的地质结构，地震数据噪声偏大断层识别难度较高。3d u-net网络在该地区的识别效果欠佳，出现大量的识别错误以及遗漏。本发明网络则能较为准确的对主要断层进行分割，识别错误以及噪声相对较少。
[0118]
综上可以看出多任务学习网络相较于传统的3d u-net网络在断层的识别结果以及模型的泛化性上都优于传统的3d u-net网络。
[0119]
同时在西南某地区数据集上可以看出深度学习网络在地质结构较为复杂的区域(timeline[600,800]的范围内)预测结果不佳。复杂地质结构产生较大的相干噪声，影响断层识别效果，预测结果产生较多的噪声以及细碎的结构，这样的断层识别结果会对后续的地质建模工作产生不利影响。因此本发明在下一节中提出了一种基于地震数据信噪比估算的预测改良方法，该方法对复杂区域内的地质结构进行过滤，过滤细碎结构，保留主要断层，缓解噪声对断层建模的影响。
[0120]
实施例6
[0121]
针对网络预测结果在地质结构复杂区域内产生较多细碎结果，对断层建模产生不利影响的情况，本发明提出一种基于地震数据信噪比估算的预测改良，能够从细碎结构提取出主要断层结构。具体包括：
[0122]
1、相关平移法估算地震数据信噪比
[0123]
在对地震数据进行信噪比估算时，首先选取一定大小的估算时窗，时窗范围内地震记录的总道数为n道，其中每道共有m个采样点。其中任意地震道记录xi(t)由有效信号si(t)和数据噪声ni(t)组成。
[0124]
xi(t)＝si(t) ni(t)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0125]
将道间的最大的互相关值q
i,i 1
(t
i,i 1
)作为功率谱：
[0126]qi,i 1
(t
i,i 1
)＝max{q
i,i 1
(τ)}
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(7)
[0127]
其中i为第i道地震道，τ为最大互相关时的时偏移量，t
为一道时间序列。
[0128]
以相关值作为地震信号的总功率谱。则平均信号能量为：
[0129][0130]
噪音能量为：
[0131][0132]
则信噪比为：
[0133][0134]
2、复杂区域预测方法
[0135]
按照上述理论在对复杂区域的预测结果进行处理时应该适度降低该区域内的置信度。传统的做法是在该区域内设定阈值，将高于阈值水平的点进行保留，低于该阈值的数据进行去除。该方法具有较大的局限性，受噪声，复杂结构等多种因素的影响，在同一窗口内可以出现预测值整体较高，或者整体偏低的情况，简单的使用单一阈值无法满足在面对不同预测窗口时将主要断层结构进行有效的提取。例如，使用较高的阈值，会导致整体预测值偏低的窗口的主要断层被过滤，使用较低的阈值会导致在整体预测值偏高的窗口保留过多的细碎结构。
[0136]
针对单一阈值的上述问题本发明提出了基于三重阈值的预测方法，在该方法中共包含三个阈值：高阈值h，低阈值l，总和阈值s。当预测结果中值低于低阈值的部分直接进行去除，高于高阈值部分的结果直接予以保留。计算出当前保留的值的总和，如果总和大于总和阈值s则返回预测结果，否则按预测值大小由高到低从高阈值和低阈值间的数据进行保留直至保留数据总和大于总和阈值s，或者没有数据满足条件。
[0137]
具体算法如下所示：
[0138][0139]
如图11为该算法的简单示例，选择大小为3*3，l＝0.3，h＝0.6，s＝2对该算法进行演示。
[0140]
最后将复杂区域预测结果中的离散点进行过滤得到最终的结果图。经过三重阈值预测方法的过滤后，能够将该区域内的地质结构进行有效的筛选，得到主要的断层结构。
[0141]
3、实际数据仿真结果
[0142]
如图12所示为某真实工区的地震数据,该数据在inline范围[100,200]间存在两条跨度较大的断层，在该地震振幅数据的上半部分的其他区域地质结构较为简单，除主要断层外无其他复杂地质结构的影响。在timeline范围[200,600]存在着一系列的地质破碎带，在该范围内存在这一系列复杂，细小的地质结构。这些地质结构产生较多的相干噪声，不利于模型对断层的预测。对该数据从timeline进行划分，每间隔100单位进行划分，对不同的划分区域使用相关时移法进行信噪比测试。
[0143]
从信噪比测试结果来看，在数据timeline[600,1000]范围内的信噪比较高，地质结构较为复杂的timeline[200,600]范围内的信噪比较低，在信噪比较低的区域地质结构也较为复杂与上文的分析结果相符合。如图13所示，直接使用模型对该地震数据进行预测，在结构复杂的范围内产生大量的细碎的断层分割结果。在主要断层之间的也出现了若干细
碎的断层影响，这些细碎断层会对后续的断层建模造成干扰，影响主要断层的走向以及趋势。在该区域内使用三重阈值对断层数据的预测结果进行过滤之后在该区域范围内的细碎断层预测结果显著减少，在该范围内主要断层的周围的干扰显著减少分割结果也较为清晰，相较于直接预测的结果更有利于断层的建模。
[0144]
图12、13中的inline表示主测线，timeline表示时间，fault probability表示断层概率.
[0145]
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的权利要求范围之内。