叠前地震数据多频段相干值获取方法及断裂预测方法与流程

1.本发明属于石油地震勘探技术领域，特别涉及一种叠前地震数据多频段相干值获取方法及断裂预测方法。


背景技术：

2.相干处理在地震资料解释中具有广泛的应用价值，可以用来做滤波处理、速度谱分析和地质体边缘检测等。自上世纪九十年代末，相干体属性作为最重要的地震属性之一被用来预测河道、断层、重力流、差异压实、盐丘等与油田勘探息息相关的地质现象，是地震定性和定量解释的关键技术。传统的相干计算方法包括相邻地震道互相关法(bahorich and farmer,1995,1996)、特征结构分析法、相似性计算法(kirlin,1992；gersztenkorn andmarfurt,1996a,b；marfurt et al.,1998)、梯度矢量场计算法(randen et al.,2000)等。这些方法的共同缺点是针对于叠加处理后的地震数据计算的。由于地震叠加处理把不同炮检距、不同方位、不同频带的信号累加到了一起，造成了很多有用弱信号对应的地质信息被消除掉。所以常规叠后相干结果并不能有效捕捉到对地震资料解释有用的微弱信号，而叠前地震相干技术通过对不同炮检距不同频带的地震信号形成的协方差矩阵进行分析，得出的相干结果可以很好捕捉这些微弱信号产生的地质异常。
3.常规叠后相干的另外一个问题是预测得到的各种地质信息混杂在一起，不能有效的将断裂单独区分开来，而且预测的断裂仅在水平切片或层切片上可解释性强，断裂在纵向剖面和三维空间的可解释性不够。


技术实现要素：

4.针对现有技术的不足，本发明的目的在于提供一种叠前地震数据多频段相干值获取方法；该相干值获取方法能够更好的保留各频段、各炮检距范围的有效地震弱信号，使用该相干值获取方法获取得到的相干值数据体进行断裂预测，可以很好的提高断裂的预测精度，同时压制与断裂无关的地质特征。
5.为了实现上述目的，本发明提供了一种叠前地震数据多频段相干值获取方法，其中，该方法包括：
6.步骤1：获取叠前地震数据体，并将叠前地震数据体按照炮检距、频带的不同进行分解得到不同炮检距、不同频带的叠前地震分数据体；
7.步骤2：选取所述叠前地震数据体中的某一个地震样点，选取分析时窗，由各个不同炮检距、不同频带的叠前地震分数据体分别构建协方差矩阵；
8.步骤3：对不同炮检距、不同频带的叠前地震分数据体对应的协方差矩阵进行求和得到基准协方差矩阵；
9.步骤4：根据基准协方差矩阵确定该地震样点的相干值；
10.步骤5：重复步骤2-4，历遍所述叠前地震数据体中的每一个地震样点，确定并记录所述每一个地震样点的相干值，形成相干值数据体。
11.在上述叠前地震数据多频段相干值获取方法中，步骤2中构建协方差矩阵采用常规方法进行即可；在一优选实施方式中，构建得到的协方差矩阵为：
[0012][0013]
其中，c
o，f，i，j
为第o个炮检距、第f个频带的叠前地震分数据体对应的协方差矩阵，即矩阵d
o,f
第j行与第i列的相关；d
o，f，k，i
为矩阵d
o,f
第k行第i列的值；d
o，f，k，j
为矩阵d
o,f
第k行第j列的值；为矩阵d
o,f
的转置矩阵第k行第i列的值；为矩阵d
o,f
的转置矩阵第k行第j列的值；其中，矩阵d
o,f
为由第o个炮检距、第f个频带的叠前地震分数据体形成的矩阵，记为：
[0014][0015]
其中，分析时窗内有横向有m道，纵向2k 1个样点。
[0016]
在上述叠前地震数据多频段相干值获取方法中，步骤4根据基准协方差矩阵确定该地震样点的相干值采用常规方法进行即可；优选地，步骤4包括：
[0017]
对基准协方差矩阵进行特征值分解，确定基准协方差矩阵的特征向量和特征值；
[0018]
由基准协方差矩阵的第一特征向量确定相干能量；
[0019]
由基准协方差矩阵确定总能量；
[0020]
基于所述相干能量以及所述总能量确定所述相干值。
[0021]
更优选地，所述相干能量通过下述公式确定：
[0022][0023][0024][0025]
其中，为矩阵csum的第一特征向量；αk为中间变量；e
pre_coh
为相干能量；为基准协方差矩阵csum的第一特征向量形成的矩阵第k行第m列的值；为基准协方差矩阵csum的第一特征向量形成的矩阵的转置在第k行第m列的值；csum
km
为基准协方差矩阵csum第k行第m列的值；分析时窗内有横向有m道，纵向2k 1个样点。
[0026]
更优选地，所述总能量通过下述公式确定：
[0027][0028]
其中，csum
km
为基准协方差矩阵csum第k行第m列的值；为基准协方差矩阵csum的转置矩阵第k行第m列的值；分析时窗内有横向有m道，纵向2k 1个样点。
[0029]
更优选地，所述相干值通过下述公式确定：
[0030][0031]
其中，e
pre_coh
为相干能量；e
pre_tot
为总能量；ε为一个等于或接近于0的值。其中，优选地，ε的绝对值小于等于0.00005。
[0032]
本发明还提供了一种基于叠前地震数据多频段相干的断裂预测方法，其中，该方法包括：
[0033]
基于上述叠前地震数据多频段相干值获取方法获取工区的相干值数据体；
[0034]
对获取的工区的相干值数据体沿着不同倾角和方位角进行扫描，根据扫描结果确定断裂面的位置。
[0035]
在上述基于叠前地震数据多频段相干的断裂预测方法中，优选地，所述对获取的工区的相干值数据体沿着不同倾角和方位角进行扫描，根据扫描结果确定断裂面的倾位置包括：
[0036]
1)选定获取的工区的相干值数据体的一个样点，沿着不同倾角和方位角的窗口进行扫描分别获取各倾角、各方位角的扫描窗口内相干值之和；
[0037]
2)确定步骤1)获取的所有扫描窗口对应的相干值之和的最小值；如果该最小值小于特定门槛值，则判断该样点为断面位置；
[0038]
3)重复步骤1)-2)直至历遍获取的工区的相干值数据体中每一个样点，完成断裂面的位置预测。
[0039]
在上述基于叠前地震数据多频段相干的断裂预测方法中，优选地，所述特定门槛值为0.25-0.3。
[0040]
在上述基于叠前地震数据多频段相干的断裂预测方法中，优选地，所述对获取的工区的相干值数据体沿着不同倾角和方位角进行扫描时，扫描倾角扫描增量不超过为5
°
，方位角扫描增量不超过5
°
。
[0041]
本发明还提供了一种叠前地震数据多频段相干值获取系统，其中，该系统包括：
[0042]
数据获取模块：用于获取叠前地震数据体，并将叠前地震数据体按照炮检距、频带的不同进行分解得到不同炮检距、不同频带的叠前地震分数据体；
[0043]
协方差矩阵构建模块：用于选取所述叠前地震数据体中的某一个地震样点，选取分析时窗，由各个不同炮检距、不同频带的叠前地震分数据体分别构建协方差矩阵；
[0044]
基准协方差矩阵构建模块：用于对不同炮检距、不同频带的叠前地震分数据体对应的协方差矩阵进行求和得到基准协方差矩阵；
[0045]
相干值计算模块：用于根据基准协方差矩阵确定该地震样点的相干值；
[0046]
相干值数据体获取模块：用于重复运行所述协方差矩阵构建模块、基准协方差矩阵构建模块和相干值计算模块，历遍所述叠前地震数据体中的每一个地震样点，确定并记录所述每一个地震样点的相干值，形成相干值数据体。
[0047]
在上述叠前地震数据多频段相干值获取系统中，协方差矩阵构建模块在构建协方差矩阵时采用常规方法进行即可；在一优选实施方式中，构建得到的协方差矩阵为：
[0048][0049]
其中，c
o，f，i，j
为第o个炮检距、第f个频带的叠前地震分数据体对应的协方差矩阵，
即矩阵d
o,f
第j行与第i列的相关；d
o，f，k，i
为矩阵d
o,f
第k行第i列的值；d
o，f，k，j
为矩阵d
o,f
第k行第j列的值；为矩阵d
o,f
的转置矩阵第k行第i列的值；为矩阵d
o,f
的转置矩阵第k行第j列的值；其中，矩阵d
o,f
为由第o个炮检距、第f个频带的叠前地震分数据体形成的矩阵，记为：
[0050][0051]
其中，分析时窗内有横向有m道，纵向2k 1个样点。
[0052]
在上述叠前地震数据多频段相干值获取系统中，相干值计算模块在根据基准协方差矩阵确定该地震样点的相干值时采用常规方法进行即可；优选地，相干值计算模块包括：
[0053]
特征值分解子模块：用于对基准协方差矩阵进行特征值分解，确定基准协方差矩阵的特征向量和特征值；
[0054]
相干能量计算子模块：用于由基准协方差矩阵的第一特征向量确定相干能量；
[0055]
总能量计算子模块：用于由基准协方差矩阵确定总能量；
[0056]
相干值计算子模块：用于基于所述相干能量以及所述总能量确定所述相干值。
[0057]
其中，所述相干能量优选通过下述公式确定：
[0058][0059][0060][0061]
其中，为矩阵csum的第一特征向量；αk为中间变量；e
pre_coh
为相干能量；为基准协方差矩阵csum的第一特征向量形成的矩阵第k行第m列的值；为基准协方差矩阵csum的第一特征向量形成的矩阵的转置在第k行第m列的值；csum
km
为基准协方差矩阵csum第k行第m列的值；分析时窗内有横向有m道，纵向2k 1个样点。
[0062]
其中，所述总能量优选通过下述公式确定：
[0063][0064]
其中，csum
km
为基准协方差矩阵csum第k行第m列的值；为基准协方差矩阵csum的转置矩阵第k行第m列的值；分析时窗内有横向有m道，纵向2k 1个样点。
[0065]
其中，所述相干值优选通过下述公式确定：
[0066][0067]
其中，e
pre_coh
为相干能量；e
pre_tot
为总能量；ε为一个等于或接近于0的值。其中，优选地，ε的绝对值小于等于0.00005。
[0068]
本发明还提供了一种基于叠前地震数据多频段相干的断裂预测系统，其中，该系统包括：
[0069]
上述叠前地震数据多频段相干值获取系统：用于获取工区的相干值数据体；
[0070]
断裂面预测模块：用于对获取的工区的相干值数据体沿着不同倾角和方位角进行扫描，根据扫描结果确定断裂面的位置。
[0071]
在上述基于叠前地震数据多频段相干的断裂预测系统中，优选地，所述断裂面预测模块包括：
[0072]
样点相干值扫描子模块：用于选定获取的工区的相干值数据体的一个样点，沿着不同倾角和方位角的窗口进行扫描分别获取各倾角、各方位角的扫描窗口内相干值之和；
[0073]
样点断面判定子模块：用于确定样点相干值扫描子模块获取的所述样点所有扫描窗口对应的相干值之和的最小值；如果该最小值小于特定门槛值，则判断该样点为断面位置；
[0074]
断裂面位置确定子模块：用于重复运行样点相干值扫描子模块、样点断面判定子模块，直至历遍获取的工区的相干值数据体中每一个样点，完成断裂面的位置预测。
[0075]
在上述基于叠前地震数据多频段相干的断裂预测系统中，优选地，对获取的工区的相干值数据体沿着不同倾角和方位角进行扫描时，扫描倾角扫描增量不超过为5
°
，方位角扫描增量不超过5
°
。
[0076]
本发明还提供了一种叠前地震数据多频段相干值获取装置，包括处理器及存储器；其中，
[0077]
存储器，用于存放计算机程序；
[0078]
处理器，用于执行存储器上所存放的程序时，实现上述叠前地震数据多频段相干值获取方法的步骤。
[0079]
本发明还提供了一种基于叠前地震数据多频段相干的断裂预测装置，包括处理器及存储器；其中，
[0080]
存储器，用于存放计算机程序；
[0081]
处理器，用于执行存储器上所存放的程序时，实现上述基于叠前地震数据多频段相干的断裂预测方法的步骤。
[0082]
本发明还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有一个或者多个程序，所述一个或者多个程序可被一个或者多个处理器执行，以实现上述叠前地震数据多频段相干值获取方法的步骤。
[0083]
本发明还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有一个或者多个程序，所述一个或者多个程序可被一个或者多个处理器执行，以实现上述基于叠前地震数据多频段相干的断裂预测方法的步骤。
[0084]
本发明提供的技术方案，以叠前地震数据为基础，对不同频带、不同炮检距的地震振幅求取协方差矩阵，进而求取叠前相干结果。本发明提供的技术方法通过叠前地震数据的多频段、不同炮检距来建立协方差矩阵从而求取了更高精度相干数据体。进一步地，本发明提供的断裂预测方法在获取得到的上述叠前相干数据基础上，通过不同倾角、不同方位扫描寻找断裂面的位置。本发明提供个的的断裂预测方法对河道等层序地质特征不敏感，可用来识别不同尺度断裂在三维空间的展布规律，尤其是对通源断裂识别具有独到之处。
该方法可作为油气储层预测的核心技术手段之一，用来指导油田勘探开发。与传统的叠后相干相比，该新方法的优势主要体现在：
[0085]
1)从多频带、不同炮检距地震数据出发计算相干，保留了所有各频段、各炮检距范围的有效地震弱信号；
[0086]
2)假设断裂在局部空间为面的特征，从不同倾角、方位角分别进行搜索来判断断裂面的空间展布，可以很好的提高断裂的预测精度，同时压制与断裂无关的地质特征。
附图说明
[0087]
图1为本发明一实施例提供的叠前地震数据多频段相干值获取方法的流程示意图。
[0088]
图2为本发明一实施例提供的叠前地震数据多频段相干值获取方法中步骤s14的优化示意图。
[0089]
图3为本发明一实施例提供的基于叠前地震数据多频段相干的断裂预测方法的流程示意图。
[0090]
图4为本发明一实施例提供的叠前地震数据多频段相干值获取系统的结构示意图。
[0091]
图5为本发明一实施例提供的叠前地震数据多频段相干值获取系统中相干值计算模块的优化示意图。
[0092]
图6为本发明一实施例提供的含基于叠前地震数据多频段相干的断裂预测系统的结构示意图。
[0093]
图7a为本发明一实施例中地震数据体剖面图。
[0094]
图7b为本发明一实施例中断裂预测剖面图。
[0095]
图7c为本发明一对比例中断裂预测剖面图。
[0096]
图8a为本发明一实施例中地震数据体切片图。
[0097]
图8b为本发明一实施例中断裂预测切片图。
[0098]
图8c为本发明一对比例中断裂预测切片图。
[0099]
图9为本发明一实施例提供的叠前地震数据多频段相干值获取装置的结构示意图。
具体实施方式
[0100]
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚完整的描述。显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明的保护范围。
[0101]
下面参考本发明的若干代表性实施方式，详细阐述本发明的原理和精神。
[0102]
参见图1，本发明一实施例提供了一种叠前地震数据多频段相干值获取方法，其中，该方法包括：
[0103]
步骤s11：获取叠前地震数据体，并将叠前地震数据体按照炮检距、频带的不同进行分解得到不同炮检距、不同频带的叠前地震分数据体；
[0104]
步骤s12：选取所述叠前地震数据体中的某一个地震样点，选取分析时窗，由各个不同炮检距、不同频带的叠前地震分数据体分别构建协方差矩阵；
[0105]
步骤s13：对不同炮检距、不同频带的叠前地震分数据体对应的协方差矩阵进行求和得到基准协方差矩阵；
[0106]
步骤s14：根据基准协方差矩阵确定该地震样点的相干值；
[0107]
步骤s15：重复步骤s12-s14，历遍所述叠前地震数据体中的每一个地震样点，确定并记录所述每一个地震样点的相干值，形成相干值数据体。
[0108]
其中，对不同炮检距、不同频带的叠前地震分数据体对应的协方差矩阵进行求和得到基准协方差矩阵可以记为：
[0109][0110]
其中，csum为基准协方差矩阵；c
o,f
为第o个炮检距、第f个频带的叠前地震分数据体对应的协方差矩阵；不同炮检距、不同频带的叠前地震分数据体包括jo个炮检距，nf个频带的叠前地震分数据体。
[0111]
其中，步骤s12中构建协方差矩阵可以采用常规方法进行。在一个实施方式中，构建得到的协方差矩阵为：
[0112][0113]
其中，d
o，f，i，j
为第o个炮检距、第f个频带的叠前地震分数据体对应的协方差矩阵，即矩阵d
o,f
第j行与第i列的相关；d
o，f，k，i
为矩阵d
o,f
第k行第i列的值；d
o，f，k，j
为矩阵d
o,f
第k行第j列的值；为矩阵d
o,f
的转置矩阵第k行第i列的值；为矩阵d
o,f
的转置矩阵第k行第j列的值；其中，矩阵d
o,f
为由第o个炮检距、第f个频带的叠前地震分数据体形成的矩阵，记为：
[0114][0115]
其中，分析时窗内有横向有m道，纵向2k 1个样点。
[0116]
其中，步骤s14中根据基准协方差矩阵确定该地震样点的相干值可以采用常规方法进行。在一实施方式中，参照图2，步骤s14包括：
[0117]
步骤s141：对基准协方差矩阵进行特征值分解，确定基准协方差矩阵的特征向量和特征值；
[0118]
步骤s142：由基准协方差矩阵的第一特征向量确定相干能量；
[0119]
步骤s143：由基准协方差矩阵确定总能量；
[0120]
步骤s144：基于所述相干能量以及所述总能量确定所述相干值。
[0121]
在一实施方式中，相干能量通过下述公式确定：
[0122][0123][0124][0125]
其中，为矩阵csum的第一特征向量；αk为中间变量；e
pre_coh
为相干能量；为基准协方差矩阵csum的第一特征向量形成的矩阵第k行第m列的值；为基准协方差矩阵csum的第一特征向量形成的矩阵的转置在第k行第m列的值；csum
km
为基准协方差矩阵csum第k行第m列的值；分析时窗内有横向有m道，纵向2k 1个样点。
[0126]
在一实施方式中，总能量通过下述公式确定：
[0127][0128]
其中，csum
km
为基准协方差矩阵csum第k行第m列的值；为基准协方差矩阵csum的转置矩阵第k行第m列的值；分析时窗内有横向有m道，纵向2k 1个样点。
[0129]
在一实施方式中，相干值通过下述公式确定：
[0130][0131]
其中，e
pre_coh
为相干能量；e
pre_tot
为总能量；ε为一个等于或接近于0的值。其中，ε可以选用绝对值小于等于0.00005的数值。
[0132]
参见图3，本发明一实施例提供了一种基于叠前地震数据多频段相干的断裂预测方法，其中，该方法包括：
[0133]
步骤s31：获取工区叠前地震数据体，并将叠前地震数据体按照炮检距、频带的不同进行分解得到不同炮检距、不同频带的叠前地震分数据体；
[0134]
步骤s32：选取所述叠前地震数据体中的某一个地震样点，选取分析时窗，由各个不同炮检距、不同频带的叠前地震分数据体分别构建协方差矩阵；
[0135]
步骤s33：对不同炮检距、不同频带的叠前地震分数据体对应的协方差矩阵进行求和得到基准协方差矩阵；
[0136]
步骤s34：根据基准协方差矩阵确定该地震样点的相干值；
[0137]
步骤s35：重复步骤s32-s34，历遍所述叠前地震数据体中的每一个地震样点，确定并记录所述每一个地震样点的相干值，形成相干值数据体；
[0138]
步骤s36：对相干值数据体沿着不同倾角和方位角进行扫描，根据扫描结果确定断裂面的位置。
[0139]
其中，对不同炮检距、不同频带的叠前地震分数据体对应的协方差矩阵进行求和得到基准协方差矩阵可以记为：
[0140]
[0141]
其中，csum为基准协方差矩阵；c
o,f
为第o个炮检距、第f个频带的叠前地震分数据体对应的协方差矩阵；不同炮检距、不同频带的叠前地震分数据体包括jo个炮检距，nf个频带的叠前地震分数据体。
[0142]
其中，步骤s32中构建协方差矩阵可以采用常规方法进行。在一个实施方式中，构建得到的协方差矩阵为：
[0143][0144]
其中，c
o，f，i，j
为第o个炮检距、第f个频带的叠前地震分数据体对应的协方差矩阵，即矩阵d
o,f
第j行与第i列的相关；d
o，f，k，i
为矩阵d
o,f
第k行第i列的值；d
o，f，k，j
为矩阵d
o,f
第k行第j列的值；为矩阵d
o,f
的转置矩阵第k行第i列的值；为矩阵d
o,f
的转置矩阵第k行第j列的值；其中，矩阵d
o,f
为由第o个炮检距、第f个频带的叠前地震分数据体形成的矩阵，记为：
[0145][0146]
其中，分析时窗内有横向有m道，纵向2k 1个样点。
[0147]
其中，步骤s34中根据基准协方差矩阵确定该地震样点的相干值可以采用常规方法进行。
[0148]
在一实施方式中，步骤s34包括：
[0149]
对基准协方差矩阵进行特征值分解，确定基准协方差矩阵的特征向量和特征值；
[0150]
由基准协方差矩阵的第一特征向量确定相干能量；
[0151]
由基准协方差矩阵确定总能量；
[0152]
基于所述相干能量以及所述总能量确定所述相干值。
[0153]
在一实施方式中，相干能量通过下述公式确定：
[0154][0155][0156][0157]
其中，为矩阵csum的第一特征向量；αk为中间变量；e
pre_coh
为相干能量；为基准协方差矩阵csum的第一特征向量形成的矩阵第k行第m列的值；为基准协方差矩阵csum的第一特征向量形成的矩阵的转置在第k行第m列的值；csum
km
为基准协方差矩阵csum第k行第m列的值；分析时窗内有横向有m道，纵向2k 1个样点。
[0158]
在一实施方式中，总能量通过下述公式确定：
[0159][0160]
其中，csum
km
为基准协方差矩阵csum第k行第m列的值；为基准协方差矩阵csum的转置矩阵第k行第m列的值；分析时窗内有横向有m道，纵向2k 1个样点。
[0161]
在一实施方式中，相干值通过下述公式确定：
[0162][0163]
其中，e
pre_coh
为相干能量；e
pre_tot
为总能量；ε为一个等于或接近于0的值。其中，ε可以选用绝对值小于等于0.00005的数值。
[0164]
在一实施方式中，步骤s36对相干值数据体沿着不同倾角和方位角进行扫描，根据扫描结果确定断裂面的位置包括：
[0165]
1)选定获取的工区的相干值数据体的一个样点，沿着不同倾角和方位角的窗口进行扫描分别获取各倾角、各方位角的扫描窗口内相干值之和；
[0166]
2)确定步骤1)获取的所有扫描窗口对应的相干值之和的最小值；如果该最小值小于特定门槛值，则判断该样点为断面位置；
[0167]
3)重复步骤1)-2)直至历遍获取的工区的相干值数据体中每一个样点，完成断裂面的位置预测。
[0168]
其中，分别计算各倾角、各方位角的扫描窗口内相干值的和可以记为：
[0169][0170]
其中，θi为倾角；为方位角；为扫描窗口内相干值的和；为θi倾角、方位角的扫描窗口中iw位置处的相干值；扫描窗口的长度为nw。
[0171]
其中，断裂强度数学上可以定义为：
[0172][0173]
其中，为扫描窗口内相干值的和。
[0174]
在一实施方式中，对相干值数据体沿着不同倾角和方位角进行扫描时，扫描倾角扫描增量不超过为5
°
，方位角扫描增量不超过5
°
。
[0175]
在一实施方式中，特定门槛值为0.25-0.3。
[0176]
本发明的又一实施例提供了一种基于叠前地震数据多频段相干的断裂预测方法，用于对工区a中的断裂进行预测，该方法包括：
[0177]
1)：获取工区三维叠前地震数据体(该叠前地震数据体如若进行叠加处理得到的叠后地震数据体剖面图如图7a所示、切片图如图8a所示)，并将叠前地震数据体按照炮检距、频带的不同进行分解得到jo个炮检距，nf个频带的叠前地震分数据体；
[0178]
2)：选取所述叠前地震数据体中的某一个地震样点，选取分析时窗，由各个不同炮检距、不同频带的叠前地震分数据体分别构建协方差矩阵；具体而言：
[0179]
假设分析窗口内横向有m道，纵向2k 1个样点；这样第o个炮检距、第f个频带的地震数据形成的矩阵为：
[0180][0181]
其中，d
o,f,i,,j
(-k≤i≤k；1≤j≤m)为第o个炮检距、第f个频带的地震数据在第i个样点、第j道的地震振幅；
[0182]
由矩阵d
o,f
组建协方差矩阵c
o,f
：
[0183][0184]
其中，c
o，f，i，j
为第o个炮检距、第f个频带的叠前地震分数据体对应的协方差矩阵，即矩阵d
o,f
第j行与第i列的相关；d
o，f，k，i
为矩阵d
o,f
第k行第i列的值；d
o，f，k，j
为矩阵d
o,f
第k行第j列的值；为矩阵d
o,f
的转置矩阵第k行第i列的值；为矩阵d
o,f
的转置矩阵第k行第j列的值；其中，矩阵d
o,f
为由第o个炮检距、第f个频带的叠前地震分数据体形成的矩阵；
[0185]
3)：对不同炮检距、不同频带的叠前地震分数据体对应的协方差矩阵进行求和得到基准协方差矩阵；具体而言：
[0186][0187]
其中，csum为基准协方差矩阵；c
o,f
为第o个炮检距、第f个频带的叠前地震分数据体对应的协方差矩阵；
[0188]
4)：根据基准协方差矩阵确定该地震样点的相干值；具体而言：
[0189]
4.1：对基准协方差矩阵进行特征值分解，确定基准协方差矩阵的特征向量和特征值；
[0190]
4.2：由基准协方差矩阵的第一特征向量确定相干能量；其中，相干能量通过下述公式确定：
[0191][0192][0193][0194]
其中，为矩阵csum的第一特征向量；αk为中间变量；e
pre_coh
为相干能量；为基准协方差矩阵csum的第一特征向量形成的矩阵第k行第m列的值；为基准协方差矩阵csum的第一特征向量形成的矩阵的转置在第k行第m列的值；csum
km
为基准协方差矩阵csum第k行第m列的值；
[0195]
4.3：由基准协方差矩阵确定总能量；其中，总能量通过下述公式确定：
[0196][0197]
其中，csum
km
为基准协方差矩阵csum第k行第m列的值；为基准协方差矩阵csum的转置矩阵第k行第m列的值；
[0198]
4.4：基于所述相干能量以及所述总能量确定所述相干值。其中，相干值通过下述公式确定：
[0199][0200]
其中，e
pre_coh
为相干能量；e
pre_tot
为总能量；ε为0.00005；
[0201]
5)：重复步骤2)-4)，历遍叠前地震数据体中的每一个地震样点，确定并记录所述每一个地震样点的相干值，形成相干值数据体；
[0202]
6)：对相干值数据体沿着不同倾角和方位角进行扫描，根据扫描结果确定断裂面的位置；具体而言：
[0203]
6.1：选定获取的工区的相干值数据体的一个样点，沿着不同倾角和方位角的窗口进行扫描分别获取各倾角、各方位角的扫描窗口内相干值之和；其中，扫描倾角扫描增量为5
°
、倾角θi的取值范围为0
°-
360
°
，方位角扫描增量为5
°
、方位角的取值范围为0
°-
360
°
；
[0204]
各倾角、各方位角的扫描窗口内相干值的和通过下述公式计算：
[0205][0206]
其中，为倾角；为方位角；为扫描窗口内相干值的和；为θi倾角、方位角的扫描窗口中iw位置处的相干值；扫描窗口的长度为nw。
[0207]
6.2：确定步骤1)获取的所有扫描窗口对应的相干值之和的最小值；如果该最小值小于特定门槛值0.25，则判断该样点为断面位置；断裂强度数学上定义为：
[0208][0209]
6.3：重复步骤6.1-6.2直至历遍获取的工区的相干值数据体中每一个样点，完成断裂面的位置预测。
[0210]
断裂预测结果的剖面图如图7b所示，断裂预测的切片图如图8b所示。
[0211]
本发明的一对比例提供了一种基于常规叠后地震相干的断裂预测方法，用于对工区a中的断裂进行预测，该方法包括：
[0212]
1)：获取工区叠前地震数据体，对工区叠前地震数据体进行处理得到工区叠后地震数据体(剖面图如图7a所示，切片图如图8a所示)；
[0213]
2)：选取所述叠后地震数据体中的某一个地震样点，选取分析时窗，构建协方差矩阵；具体而言：
[0214]
假设分析窗口内横向有m道，纵向2k 1个样点；这些点的地震振幅可组成矩阵d：
[0215][0216]
其中，d
ij
(-k≤i≤k；1≤j≤m)为第i个样点、第j道的地震振幅；
[0217]
由d
ij
可组建矩阵c：
[0218][0219]
其中，c
ij
为协方差矩阵，即矩阵d第j行与第i列的相关；d
ki
为矩阵d第k行第i列的值；d
kj
为矩阵d第k行第j列的值；为矩阵d的转置矩阵第k行第i列的值；为矩阵d的转置矩阵第k行第j列的值；
[0220]
3)：根据协方差矩阵确定该地震样点的相干值；具体而言：
[0221]
3.1：对协方差矩阵进行特征值分解，确定基准协方差矩阵的特征向量和特征值；
[0222]
3.2：由基准协方差矩阵的第一特征向量确定相干能量；其中，相干能量通过下述公式确定：
[0223][0224][0225][0226]
其中，为矩阵c的第一特征向量；βk为中间变量；e
coh
为相干能量；为协方差矩阵c的第一特征向量形成的矩阵第k行、第m列的值；为协方差矩阵c的第一特征向量形成的矩阵的转置在第k行，第m列的值；c
km
为协方差矩阵c第k行第m列的值；
[0227]
3.3：由协方差矩阵确定总能量；其中，总能量通过下述公式确定：
[0228][0229]
其中，e
tot
为总能量；c
km
为协方差矩阵c第k行第m列的值；为协方差矩阵c的转置矩阵第k行第m列的值；
[0230]
3.4：基于所述相干能量以及所述总能量确定所述相干值。其中，相干值通过下述公式确定：
[0231][0232]
其中，e
coh
为相干能量；e
tot
为总能量；ε为0.00005；
[0233]
4)：重复步骤2)-3)，历遍叠后地震数据体中的每一个地震样点，确定并记录所述每一个地震样点的相干值，形成相干值数据体；
[0234]
5)：对相干值数据体进行扫描，根据扫描结果确定断裂面；
[0235]
断裂预测结果的剖面图如图7c所示，断裂预测的切片图如图8c所示。
[0236]
由图7a-图7c的垂直剖面可见，浅层的走滑断裂断距大，能够产生地震同相轴明显错断，常规相干能识别断点位置，但是在深层常规相干识别断裂效果较差。另外，常规相干在剖面上识别的断裂整体解释性不强。相比之下，本发明提供的断裂预测方法能够较好克服常规相干的问题，在剖面上断裂的延展性和可解释性较强，可以较好的预测深大断裂的发育。
[0237]
图8a-图8c为地震数据、本发明提供的预测方法和基于常规叠后相干法进行的断裂预测在沿层切片上的展示。由于地震数据在叠加时消除了好多与小断裂有关的微弱信号，常规叠后相干不能较好的识别尺度较小的断裂，反之由于叠前相干算法兼顾了不同炮检距、不同频段的地震信号，较好保留了小断裂对应的微弱信号，所以本发明提供的断裂预测方法可以很好的捕捉到小尺度断裂信息，同时由于采用了相干体的空间搜索，断裂预测结果的整体性、延展性、可靠性更好。
[0238]
本发明实施例还提供了一种叠前地震数据多频段相干值获取系统，优选地，该系统用于实现上述的方法实施例。
[0239]
图4是根据本发明实施例的叠前地震数据多频段相干值获取系统的结构框图，如图4所示，该系统包括：
[0240]
数据获取模块41：用于获取叠前地震数据体，并将叠前地震数据体按照炮检距、频带的不同进行分解得到不同炮检距、不同频带的叠前地震分数据体；
[0241]
协方差矩阵构建模块42：用于选取所述叠前地震数据体中的某一个地震样点，选取分析时窗，由各个不同炮检距、不同频带的叠前地震分数据体分别构建协方差矩阵；
[0242]
基准协方差矩阵构建模块43：用于对不同炮检距、不同频带的叠前地震分数据体对应的协方差矩阵进行求和得到基准协方差矩阵；
[0243]
相干值计算模块44：用于根据基准协方差矩阵确定该地震样点的相干值；
[0244]
相干值数据体获取模块45：用于重复运行协方差矩阵构建模块42、基准协方差矩阵构建模块43和相干值计算模块44，历遍所述叠前地震数据体中的每一个地震样点，确定并记录所述每一个地震样点的相干值，形成相干值数据体。
[0245]
其中，协方差矩阵构建模块43在构建协方差矩阵时可以采用常规方法进行。
[0246]
在一实施方式中，构建得到的协方差矩阵为：
[0247][0248]
其中，c
o，f，i，j
为第o个炮检距、第f个频带的叠前地震分数据体对应的协方差矩阵，即矩阵d
o,f
第j行与第i列的相关；d
o，f，k，i
为矩阵d
o,f
第k行第i列的值；d
o，f，k，j
为矩阵d
o,f
第k行第j列的值；为矩阵d
o,f
的转置矩阵第k行第i列的值；为矩阵d
o,f
的转置矩阵第k行第j列的值；其中，矩阵d
o,f
为由第o个炮检距、第f个频带的叠前地震分数据体形成的矩阵，记为：
[0249][0250]
其中，分析时窗内有横向有m道，纵向2k 1个样点。
[0251]
其中，相干值计算模块44在根据基准协方差矩阵确定该地震样点的相干值时可以采用常规方法进行。
[0252]
在一实施方式中，参照图5，相干值计算模块44包括：
[0253]
特征值分解子模441：用于对基准协方差矩阵进行特征值分解，确定基准协方差矩阵的特征向量和特征值；
[0254]
相干能量计算子模块442：用于由基准协方差矩阵的第一特征向量确定相干能量；
[0255]
总能量计算子模块443：用于由基准协方差矩阵确定总能量；
[0256]
相干值计算子模块444：用于基于所述相干能量以及所述总能量确定所述相干值。
[0257]
其中，所述相干能量可以通过下述公式确定：
[0258][0259][0260][0261]
其中，为矩阵csum的第一特征向量；αk为中间变量；e
pre_coh
为相干能量；为基准协方差矩阵csum的第一特征向量形成的矩阵第k行第m列的值；为基准协方差矩阵csum的第一特征向量形成的矩阵的转置在第k行第m列的值；csum
km
为基准协方差矩阵csum第k行第m列的值；分析时窗内有横向有m道，纵向2k 1个样点。
[0262]
其中，所述总能量可以通过下述公式确定：
[0263][0264]
其中，csum
km
为基准协方差矩阵csum第k行第m列的值；为基准协方差矩阵csum的转置矩阵第k行第m列的值；分析时窗内有横向有m道，纵向2k 1个样点。
[0265]
其中，所述相干值可以通过下述公式确定：
[0266][0267]
其中，e
pre_coh
为相干能量；e
pre_tot
为总能量；ε为一个等于或接近于0的值。其中，ε优选绝对值小于等于0.00005的数值。
[0268]
本发明实施例还提供了一种基于叠前地震数据多频段相干的断裂预测系统，优选地，该系统用于实现上述的方法实施例。
[0269]
图6是根据本发明实施例的基于叠前地震数据多频段相干的断裂预测系统的结构框图，如图6所示，该系统包括：
[0270]
数据获取模块61：用于获取叠前地震数据体，并将叠前地震数据体按照炮检距、频带的不同进行分解得到不同炮检距、不同频带的叠前地震分数据体；
[0271]
协方差矩阵构建模块62：用于选取所述叠前地震数据体中的某一个地震样点，选取分析时窗，由各个不同炮检距、不同频带的叠前地震分数据体分别构建协方差矩阵；
[0272]
基准协方差矩阵构建模块63：用于对不同炮检距、不同频带的叠前地震分数据体对应的协方差矩阵进行求和得到基准协方差矩阵；
[0273]
相干值计算模块64：用于根据基准协方差矩阵确定该地震样点的相干值；
[0274]
相干值数据体获取模块65：用于重复运行协方差矩阵构建模块62、基准协方差矩阵构建模块63和相干值计算模块64，历遍所述叠前地震数据体中的每一个地震样点，确定并记录所述每一个地震样点的相干值，形成相干值数据体；
[0275]
断裂面预测模块66：用于对相干值数据体沿着不同倾角和方位角进行扫描，根据扫描结果确定断裂面的位置。
[0276]
其中，协方差矩阵构建模块62在构建协方差矩阵时可以采用常规方法进行。
[0277]
在一实施方式中，构建得到的协方差矩阵为：
[0278][0279]
其中，c
o，f，i，j
为第o个炮检距、第f个频带的叠前地震分数据体对应的协方差矩阵，即矩阵d
o,f
第j行与第i列的相关；d
o，f，k，i
为矩阵d
o,f
第k行第i列的值；d
o，f，k，i
为矩阵d
o,f
第k行第j列的值；为矩阵d
o,f
的转置矩阵第k行第i列的值；为矩阵d
o,f
的转置矩阵第k行第j列的值；其中，矩阵d
o,f
为由第o个炮检距、第f个频带的叠前地震分数据体形成的矩阵，记为：
[0280][0281]
其中，分析时窗内有横向有m道，纵向2k 1个样点。
[0282]
其中，相干值计算模块64在根据基准协方差矩阵确定该地震样点的相干值时可以采用常规方法进行。
[0283]
在一实施方式中，相干值计算模块64包括：
[0284]
特征值分解子模：用于对基准协方差矩阵进行特征值分解，确定基准协方差矩阵的特征向量和特征值；
[0285]
相干能量计算子模块：用于由基准协方差矩阵的第一特征向量确定相干能量；
[0286]
总能量计算子模块：用于由基准协方差矩阵确定总能量；
[0287]
相干值计算子模块：用于基于所述相干能量以及所述总能量确定所述相干值。
[0288]
其中，所述相干能量可以通过下述公式确定：
[0289][0290][0291][0292]
其中，为矩阵csum的第一特征向量；αk为中间变量；e
pre_coh
为相干能量；为基准协方差矩阵csum的第一特征向量形成的矩阵第k行第m列的值；为基准协方差矩阵csum的第一特征向量形成的矩阵的转置在第k行第m列的值；csum
km
为基准协方差矩阵csum第k行第m列的值；分析时窗内有横向有m道，纵向2k 1个样点。
[0293]
其中，所述总能量可以通过下述公式确定：
[0294][0295]
其中，csum
km
为基准协方差矩阵csum第k行第m列的值；为基准协方差矩阵csum的转置矩阵第k行第m列的值；分析时窗内有横向有m道，纵向2k 1个样点。
[0296]
其中，所述相干值可以通过下述公式确定：
[0297][0298]
其中，e
pre_coh
为相干能量；e
pre_tot
为总能量；ε为一个等于或接近于0的值。其中，ε优选绝对值小于等于0.00005的数值。
[0299]
在一实施方式中，断裂面预测模块66包括：
[0300]
样点相干值扫描子模块：用于选定获取的工区的相干值数据体的一个样点，沿着不同倾角和方位角的窗口进行扫描分别获取各倾角、各方位角的扫描窗口内相干值之和；
[0301]
样点断面判定子模块：用于确定样点相干值扫描子模块获取的所述样点所有扫描窗口对应的相干值之和的最小值；如果该最小值小于特定门槛值，则判断该样点为断面位置；
[0302]
断裂面位置确定子模块：用于重复运行样点相干值扫描子模块、样点断面判定子模块，直至历遍获取的工区的相干值数据体中每一个样点，完成断裂面的位置预测。
[0303]
在一实施方式中，对获取的工区的相干值数据体沿着不同倾角和方位角进行扫描时，扫描倾角扫描增量不超过为5
°
，方位角扫描增量不超过5
°
。
[0304]
在一实施方式中，所述特定门槛值为0.25-0.3。
[0305]
图9是根据本发明实施例的叠前地震数据多频段相干值获取装置的示意图。图9所示的叠前地震数据多频段相干值获取装置为通用数据处理装置，其包含通用的计算机硬件结构，其至少包含处理器1000、存储器1111；所述处理器1000用于执行所述存储器中存储的叠前地震数据多频段相干值获取程序，以实现各方法实施例所述的叠前地震数据多频段相干值获取方法(具体方法参见上述方法实施例的描述，在此不再赘述)。
[0306]
本发明一实施例提供了一种基于叠前地震数据多频段相干的断裂预测装置，其包
含通用的计算机硬件结构，其至少包含处理器、存储器；所述处理器用于执行所述存储器中存储的基于叠前地震数据多频段相干的断裂预测程序，以实现各方法实施例所述的基于叠前地震数据多频段相干的断裂预测方法(具体方法参见上述方法实施例的描述，在此不再赘述)。
[0307]
本发明实施例还提供了一种计算机可读存储介质，所述存储介质存储有一个或者多个程序，所述一个或者多个程序可被一个或者多个处理器执行，以实现各方法实施例所述的叠前地震数据多频段相干值获取方法(具体方法参见上述方法实施例的描述，在此不再赘述)。
[0308]
本发明实施例还提供了一种计算机可读存储介质，所述存储介质存储有一个或者多个程序，所述一个或者多个程序可被一个或者多个处理器执行，以实现各方法实施例所述的基于叠前地震数据多频段相干的断裂预测方法(具体方法参见上述方法实施例的描述，在此不再赘述)。
[0309]
以上参照附图描述了本发明的优选实施方式。这些实施方式的许多特征和优点根据该详细的说明书是清楚的，因此权利要求旨在覆盖这些实施方式的落入其真实精神和范围内的所有这些特征和优点。此外，由于本领域的技术人员容易想到很多修改和改变，因此不是要将本发明的实施方式限于所例示和描述的精确结构和操作，而是可以涵盖落入其范围内的所有合适修改和等同物。
[0310]
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0311]
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0312]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0313]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0314]
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内
容不应理解为对本发明的限制。