暗河裂缝发育预测方法、系统、存储介质以及电子设备与流程

1.本发明属于地质勘探技术领域，尤其涉及一种暗河裂缝发育预测方法、系统、存储介质以及电子设备。


背景技术：

2.古暗河的形成与断裂破碎带相关，可溶性大气水经断裂进入碳酸盐岩地层，对形成的断裂系统进行溶蚀改造，最终形成管道状水流通道。断裂及其伴生裂缝增加了地表水及地下水与碳酸盐岩的接触面积和溶蚀范围，从而改善了碳酸盐岩的渗流作用。断裂及其伴生断裂带是裂缝与岩溶发育的密集带，其中，断裂作为主要的运移通道，一方面可以使地下水发生侧向运移，另一方面为地表淡水的垂向运移提供了条件，使得断裂带附近多发育岩溶。
3.可见，断裂及裂缝对岩溶发育具有重要的控制作用，不同尺度的裂缝发育对古暗河储层具有不同程度的影响和改善作用。因此，建立一种能够准确分析不同尺度裂缝对古暗河的影响的方法，对古暗河有利储集体发育区的预测具有重要的意义。


技术实现要素：

4.本发明正是基于不同尺度的裂缝发育对古暗河储层具有不同程度的影响和改善作用的技术问题，提出了一种暗河裂缝发育预测方法、系统、存储介质以及电子设备。
5.第一方面，本发明实施例提供了一种暗河裂缝发育预测方法，包括：
6.获取暗河区域的叠后地震数据；
7.对所述叠后地震数据分别进行多个预设方向的曲波变换，以获得多个预设方向下的地震数据体；
8.分别对每个预设方向下的地震数据体进行地震相干处理，以获得每个预设方向下的相干体；
9.根据各个预设方向下的相干体，获得该暗河区域的裂缝发育密度数据体以及裂缝发育方向数据体；
10.根据所述裂缝发育密度数据体，获取目标岩层对应的裂缝发育密度平面属性；以及
11.根据所述裂缝发育方向数据体，获取目标岩层对应的裂缝发育方向平面属性；
12.将所述裂缝发育密度平面属性和所述裂缝发育方向平面属性与该暗河区域对应的暗河展布图进行叠合显示，获得第一叠合展布图，以根据所述第一叠合展布图分析小尺度裂缝在暗河区域的发育情况。
13.可选地，所述方法还包括：
14.从所述裂缝发育密度数据体中提取目标岩层的蚂蚁体属性，获得该目标岩层的裂缝密度蚂蚁体；
15.将所述裂缝密度蚂蚁体与该暗河区域对应的暗河展布图进行叠合显示，获得第二
叠合展布图，以根据所述第二叠合展布图分析大尺度裂缝在暗河区域的发育情况。
16.可选地，从所述裂缝发育密度数据体中提取目标岩层的蚂蚁体属性，获得该目标岩层的裂缝密度蚂蚁体之后，所述方法还包括：
17.将所述裂缝密度蚂蚁体与所述第一叠合展布图进行叠合显示，获得第三叠合展布图，以根据所述第三叠合展布图分析小尺度裂缝以及大尺度裂缝在暗河区域的发育情况。
18.可选地，所述方法还包括：
19.获取所述暗河区域的溶洞分布图；
20.将所述溶洞分布图与所述第三叠合展布图进行叠合显示，获得第四叠合展布图，以根据所述第四叠合展布图分析小尺度裂缝以及大尺度裂缝在暗河区域的溶洞中的发育情况。
21.可选地，根据所述裂缝发育密度数据体，获取目标岩层对应的裂缝发育密度平面属性；以及
22.根据所述裂缝发育方向数据体，获取目标岩层对应的裂缝发育方向平面属性，包括：
23.从所述裂缝发育密度数据体中提取目标岩层上与裂缝相关的沿层地震属性，以根据所述沿层地震属性获得所述目标岩层对应的裂缝发育密度平面属性；以及
24.从所述裂缝发育方向数据体中提取目标岩层上与裂缝相关的沿层地震属性，根据所述沿层地震属性获得所述目标岩层对应的裂缝发育方向平面属性。
25.可选地，所述预设方向包括0
°
、45
°
、90
°
以及135
°
中的至少一个方向。
26.第二方面，本发明实施例还提供了一种暗河裂缝发育预测系统，包括：
27.获取单元，用于获取暗河区域的叠后地震数据；
28.曲波变换模块，用于对所述叠后地震数据分别进行多个预设方向的曲波变换，以获得多个预设方向下的地震数据体；
29.相干处理模块，用于分别对每个预设方向下的地震数据体进行地震相干处理，以获得每个预设方向下的相干体；
30.裂缝数据获取模块，用于根据所述裂缝发育密度数据体，获取目标岩层对应的裂缝发育密度平面属性；以及
31.用于根据所述裂缝发育方向数据体，获取目标岩层对应的裂缝发育方向平面属性；
32.叠合显示模块，用于将所述裂缝发育密度平面属性和所述裂缝发育方向平面属性与该暗河区域对应的暗河展布图进行叠合显示，获得第一叠合展布图，以根据所述第一叠合展布图分析小尺度裂缝在暗河区域的发育情况。
33.可选地，所述系统还包括：
34.蚂蚁体提取模块，用于从所述裂缝发育密度数据体中提取目标岩层的蚂蚁体属性，获得该目标岩层的裂缝密度蚂蚁体；
35.所述叠合显示模块还用于将所述裂缝密度蚂蚁体与该暗河区域对应的暗河展布图进行叠合显示，获得第二叠合展布图，以根据所述第二叠合展布图分析大尺度裂缝在暗河区域的发育情况。
36.可选地，所述叠合显示模块还用于将所述裂缝密度蚂蚁体与所述第一叠合展布图
进行叠合显示，获得第三叠合展布图，以根据所述第三叠合展布图分析小尺度裂缝以及大尺度裂缝在暗河区域的发育情况。
37.可选地，所述获取模块还用于获取所述暗河区域的溶洞分布图；
38.所述叠合显示模块还用于将所述溶洞分布图与所述第三叠合展布图进行叠合显示，获得第四叠合展布图，以根据所述第四叠合展布图分析小尺度裂缝以及大尺度裂缝在暗河区域的溶洞中的发育情况。
39.可选地，所述裂缝数据获取模块具体用于从所述裂缝发育密度数据体中提取目标岩层上与裂缝相关的沿层地震属性，根据所述沿层地震属性获得所述目标岩层对应的裂缝发育密度平面属性；以及
40.从所述裂缝发育方向数据体中提取目标岩层上与裂缝相关的沿层地震属性，根据所述沿层地震属性获得所述目标岩层对应的裂缝发育方向平面属性。
41.可选地，所述预设方向包括0
°
、45
°
、90
°
以及135
°
中的至少一个方向。
42.第三方面，本发明实施例还提供了一种存储介质，所述存储介质上存储有程序代码，所述程序代码被处理器执行时，实现如上述实施例中任一项所述的暗河裂缝发育预测方法。
43.第四方面，本发明实施例还提供了一种电子设备，所述电子设备包括存储器、处理器，所述存储器上存储有可在所述处理器上运行的程序代码，所述程序代码被所述处理器执行时，实现如上述实施例中任一项所述的暗河裂缝发育预测方法。
44.在本发明实施例提供的一种暗河裂缝发育预测方法，通过对暗河区域的叠后地震数据进行多个预设方向的曲波变换，从而获得多个预设方向下的地震数据体；进而分别对多个预设方向下的地震数据体进行地震相干处理，从而获得该暗河区域的裂缝发育密度数据体以及裂缝发育方向数据体，并且从裂缝发育密度数据体以及裂缝发育方向数据体中分别提取目标岩层对应的裂缝发育密度平面属性以及目标岩层对应的裂缝发育方向平面属性，从而利用目标岩层对应的裂缝发育密度平面属性以及裂缝发育方向平面属性来分析小尺度裂缝在暗河区域的分布状况。通过将目标岩层对应的裂缝发育密度平面属性以及裂缝发育方向平面属性与暗河区域的暗河展布图进行叠合显示，可以直观地了解小尺度裂缝在暗河各个区域的发育情况，从而为后续钻井设置提供数据依据。
附图说明
45.通过结合附图阅读下文示例性实施例的详细描述可更好地理解本公开的范围。其中所包括的附图是：
46.图1示出了本发明实施例一提出的一种暗河裂缝发育预测方法的流程示意图；
47.图2示出了本发明实施例一提出的第一叠合展布图的示意图；
48.图3示出了本发明实施例二提出的一种暗河裂缝发育预测方法的流程示意图；
49.图4示出了本发明实施例二提出的第二叠合展布图的示意图；
50.图5示出了本发明实施例三提出的一种暗河裂缝发育预测方法的流程示意图；
51.图6示出了本发明实施例三提出裂缝发育密度平面属性、裂缝发育方向平面属性以及裂缝密度蚂蚁体叠合显示的叠合图的示意图；
52.图7示出了本发明实施例三提出的第四叠合展布图的示意图；
53.图8示出了某工区中测井的裂缝地震预测示意图。
具体实施方式
54.为使本发明的目的、技术方案和优点更加清楚，以下将结合附图及实施例来详细说明本发明的实施方法，借此对本发明如何应用技术手段来解决技术问题，并达成技术效果的实现过程能充分理解并据以实施。
55.在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是，本发明还可以采用其他不同于在此描述的其他方式来实施，因此，本发明的保护范围并不受下面公开的具体实施例的限制。
56.实施例一
57.根据本发明的实施例，提供了一种暗河裂缝发育预测方法，图1示出了本发明实施例一提出的一种暗河裂缝发育预测方法的流程示意图，如图1所示，该暗河裂缝发育预测方法可以包括：步骤110至步骤170。
58.在步骤110中，获取暗河区域的叠后地震数据。
59.这里，利用采集设备对地震波信号进行采集得到地震数据体，然后对该地震数据体进行叠加处理得到叠后地震数据。通过对叠后地震数据进行处理，能够从中提取出反映地震属性的数据，该地震属性是有关地震波的几何形态、运动学特征、动力学特征和统计学特性的特殊度量。
60.其中，地震波信号在地下地层传播过程中受地层岩性、物性等因素影响会产生相应变化，它是对地下储层综合特征的一种复杂反映。地下地层岩石物理等性质的空间变化，必然导致地震反射波特征的变化，进而影响到地震属性。特别是当储层含油气时，其地震响应特征会发生相应变化，相对应的地震属性也会有所体现。地震属性技术预测油气的理论基础是：地震属性携带有地下地层信息，同时地震属性和储层的含油气性之间也必然存在某种形式的内在联系。
61.在步骤120中，对所述叠后地震数据分别进行多个预设方向的曲波变换，以获得多个预设方向下的地震数据体。
62.这里，由于裂缝的发育方向以及发育密度具有随机性，同一个岩层的裂缝方向发育也会存在随机性。通过对所述叠后地震数据分别进行多个预设方向的曲波变换，可以获得全方位描述裂缝发育密度以及发育方向的数据。
63.其中，曲波变换是继小波变换、脊波变换等发展起来的一种能够对高维信号提供近乎最优稀疏表示的数学变换方法，具有多分辨率特性，时频局部性、多方向性和各向异性，而且克服了小波变换表示边缘、轮廓等高维奇异时存在的局限性。曲波变换作为一种多尺度多方向分析方法，依据曲波变换角度、方向、位置等特征，对信号进行分析，能够有效保留地震信号中的有效信号。
64.这里，所述预设方向包括0
°
、45
°
、90
°
以及135
°
中的至少一个方向。
65.通过对暗河区域的叠后地震数据进行0
°
、45
°
、90
°
以及135
°
方向下的曲波变换，能够获取到不同角度方向下的地震数据体，从而通过不同角度方向下的地震数据体反映岩层中裂缝的发育密度以及发育方向。
66.在步骤130中，分别对每个预设方向下的地震数据体进行地震相干处理，以获得每
个预设方向下的相干体。
67.这里，相干体是指由常规地震数据体经过相干处理而得到的一个新的数据体，通过对各个方向下的地震数据体进行相干处理，获得各个方向下的地震数据体对应的相干体。
68.其中，相干体技术是一项重要的地震属性技术，它通过计算相邻地震道波形的相似性将三维地震数据体转换为相干数据体，突出了波形的不连续性特征。因此，相干体能够度量由于构造、地层、岩性、油气等因素的变化所引起的地震响应的横向变化，从而有效揭示断层、裂缝、岩性体边缘和不整合等地质现象，反映地质异常特征的平面展布。
69.优选地，可以利用第三代相干算法对每个预设方向下的地震数据体进行地震相干处理，以分别对经过不同预设方向曲波变换得到的地震数据进行地震相干处理，从而获得每个预设方向下的相干体。
70.在步骤140中，根据各个预设方向下的相干体，获得该暗河区域的裂缝发育密度数据体以及裂缝发育方向数据体。
71.这里，该暗河区域的裂缝发育密度数据体指的是在裂缝在该暗河区域中的发育密度数据，该裂缝发育方向数据体则指的是在裂缝在该暗河区域中的发育方向数据。
72.其中，在获得各个方向下的相干体后，对各个相干体进行进一步的分析比较，寻找出变化最强的数据子体及其对应的方向，得到裂缝发育带及其走向的分布图。
73.获得该暗河区域的裂缝发育密度数据体以及裂缝发育方向数据体的整个具体流程是：对叠后地震数据体进行曲波变换，得到多个方向的地震数据体，然后利用相干体算法对多个方向的地震数据体进行相干处理获得各个方向下的相干体，然后计算各个相干体中的裂缝发育带走向以及强度参数，分别得到所述裂缝发育密度数据体以及所述裂缝发育方向数据体。
74.在步骤150中，根据所述裂缝发育密度数据体，获取目标岩层对应的裂缝发育密度平面属性。
75.这里，目标岩层是指将要研究的储层，即在采油或采气时，油气可能会蕴藏的岩层。
76.其中，裂缝发育密度数据体中包括多个不同深度的岩层的裂缝发育密度，从裂缝发育密度数据体中获取目标岩层对应的裂缝发育密度平面属性，从而根据该裂缝发育密度平面属性对目标岩层所在平面的裂缝密度进行分析。该裂缝发育密度平面属性可以是目标岩层所在平面的裂缝密度展布图，用于反映裂缝在目标岩层的发育密度。
77.值得说明的是，该裂缝发育密度平面属性反映了裂缝在目标岩层的发育密度，即目标岩层的裂缝密度。裂缝密度是表示裂缝发育程度的一个概念值。由于计算时使用的资料不同,又有线密度、面密度、体积密度之分,但在实际使用中并不那么严密，各种计算方法都是为了反映裂缝发育程度。根据这个概念值，可将研究区划分出裂缝发育带、次发育带及不发育带。
78.在一个可选的实施方式中，步骤150中，根据所述裂缝发育密度数据体，获取目标岩层对应的裂缝发育密度平面属性，包括：
79.从所述裂缝发育密度数据体中提取目标岩层上与裂缝相关的沿层地震属性，以根据所述沿层地震属性获得所述目标岩层对应的裂缝发育密度平面属性。
80.这里，沿层地震属性是以解释层位为基础，沿目标层段开一时窗，对时窗内的记录作自相关、自回归等统计特征分析而提取的地震属性，其主要包括振幅、倾角分析、方位分析、边缘检测、差异检测等沿层构造属性和频率、波阻抗等沿层岩性。该沿层地震属性是一种要用来对时窗内包含的不止一个蜂或谷的地震指令进行定量描述的属性。
81.其中，通过从所述裂缝发育密度数据体中提取目标岩层上与裂缝相关的沿层地震属性，然后将该沿层地震属性用来描述目标岩层的裂缝发育密度，从而获得该目标岩层的裂缝发育密度平面属性，即得到裂缝在该目标岩层的发育密度情况。
82.值得说明的是，在这里仅研究暗河区域裂缝的发育情况，因此，提取目标岩层上与裂缝相关的沿层地震属性，从而减少计算量。
83.在步骤160中，根据所述裂缝发育方向数据体，获取目标岩层对应的裂缝发育方向平面属性。
84.这里，裂缝发育密度数据体中包括多个不同深度的岩层的裂缝发育方向平面属性，从裂缝发育密度数据体中获取目标岩层对应的裂缝发育密度平面属性，从而根据该裂缝发育方向平面属性对目标岩层所在平面的裂缝方向进行分析。该裂缝发育密度平面属性可以是目标岩层所在平面的裂缝方向展布图，用于反映目标岩层的裂缝的发育方向。
85.应当理解的是，裂缝性气藏裂缝是天然气主要的渗流通道，裂缝的发育方向和发育密度决定了气井的产能。通过对目标岩层的裂缝发育密度以及裂缝发育方向进行分析，可以得到目标岩层的裂缝发育情况，从而为后续探井的设置提供依据。
86.其中，该裂缝发育方向平面属性反映了裂缝在目标岩层的发育方向，即目标岩层的裂缝发育方向。
87.在一个可选的实施方式中，步骤160中，根据所述裂缝发育方向数据体，获取目标岩层对应的裂缝发育方向平面属性，包括：
88.从所述裂缝发育方向数据体中提取目标岩层上与裂缝相关的沿层地震属性，根据所述沿层地震属性获得所述目标岩层对应的裂缝发育方向平面属性。
89.这里，沿层地震属性的一般性解释已在上述实施方式中进行了详细说明，在此不再赘述。
90.其中，通过从所述裂缝发育方向数据体中提取目标岩层上与裂缝相关的沿层地震属性，然后将该沿层地震属性用来描述目标岩层的裂缝发育方向，从而获得该目标岩层的裂缝发育方向平面属性，即得到裂缝在该目标岩层的发育方向情况。
91.值得说明的是，在这里仅研究暗河区域裂缝的发育情况，因此，提取目标岩层上与裂缝相关的沿层地震属性，从而减少计算量。
92.值得说明的是，步骤150与步骤160的执行可以没有先后顺序，可以同时执行步骤150与步骤160。
93.在步骤170中，将所述裂缝发育密度平面属性和所述裂缝发育方向平面属性与该暗河区域对应的暗河展布图进行叠合显示，获得第一叠合展布图，以根据所述第一叠合展布图分析小尺度裂缝在暗河区域的发育情况。
94.这里，裂缝裂缝按尺度进行划分，可分为微裂缝、小尺度裂缝以及大尺度裂缝。
95.其中，微裂缝是指在岩心薄片上能观察到的裂缝，一般是微米级的裂缝，微裂缝在储层中数量庞大，通常是与基质相伴生的，因此微裂缝实际是基质的一部分。
96.小尺度裂缝一般是指在岩心和成像测井上能观察到的裂缝，长度一般在几十厘米到几十米。成像测井和岩心上看到的是小尺度裂缝的高度。小尺度裂缝在三维空间交织成网，形成油气的渗流通道。裂缝密度、裂缝开度、裂缝长度和裂缝方位变化程度是决定裂缝渗透率的关键因素。
97.大尺度裂缝一般是指在地震级别的裂缝，长度一般在几十米到上千米。大尺度裂缝横向延伸长度大，垂向切割层位深，开度大的大尺度裂缝有很高的渗透率率，在油藏中小尺度裂缝提供了油藏的渗透率，而大尺度裂缝决定了油藏的非均质性。
98.其中，该该暗河区域对应的暗河展布图是指暗河在该区域的展布图，其可以通过地震资料获得，为现有技术，在此不多赘述。
99.在步骤170中，对于小尺度裂缝的地震预测可以通过裂缝发育密度平面属性和裂缝发育方向平面属性来进行。通过将所述裂缝发育密度平面属性和所述裂缝发育方向平面属性与该暗河区域对应的暗河展布图进行叠合显示，得到第一叠合展布图，可以根据该第一叠合展布图来分析小尺度裂缝。
100.图2示出了本发明实施例一提出的第一叠合展布图的示意图，如图2所示，可以看出在暗河发育区域周边的小尺度裂缝发育情况。例如，在暗河区域的拐弯处，小尺度裂缝发育密度大，对于暗河的改向起到一定的作用。
101.在本实施例中，通过对暗河区域的叠后地震数据进行多个预设方向的曲波变换，从而获得多个预设方向下的地震数据体；进而分别对多个预设方向下的地震数据体进行地震相干处理，从而获得该暗河区域的裂缝发育密度数据体以及裂缝发育方向数据体，并且从裂缝发育密度数据体以及裂缝发育方向数据体中分别提取目标岩层对应的裂缝发育密度平面属性以及目标岩层对应的裂缝发育方向平面属性，从而利用目标岩层对应的裂缝发育密度平面属性以及裂缝发育方向平面属性来分析小尺度裂缝在暗河区域的分布状况。通过将目标岩层对应的裂缝发育密度平面属性以及裂缝发育方向平面属性与暗河区域的暗河展布图进行叠合显示，可以直观地了解小尺度裂缝在暗河各个区域的发育情况，从而为后续钻井设置提供数据依据。
102.实施例二
103.在上述实施例的基础上，本发明的实施例二还可以提供一种暗河裂缝发育预测方法。图3示出了本发明实施例二提出的一种暗河裂缝发育预测方法的流程示意图，如图3所示，该暗河裂缝发育预测方法可以包括：步骤210至步骤260。
104.在步骤210中，获取暗河区域的叠后地震数据。
105.这里，利用采集设备对地震波信号进行采集得到地震数据体，然后对该地震数据体进行叠加处理得到叠后地震数据。通过对叠后地震数据进行处理，能够从中提取出反映地震属性的数据，该地震属性是有关地震波的几何形态、运动学特征、动力学特征和统计学特性的特殊度量。
106.其中，地震波信号在地下地层传播过程中受地层岩性、物性等因素影响会产生相应变化，它是对地下储层综合特征的一种复杂反映。地下地层岩石物理等性质的空间变化，必然导致地震反射波特征的变化，进而影响到地震属性。特别是当储层含油气时，其地震响应特征会发生相应变化，相对应的地震属性也会有所体现。地震属性技术预测油气的理论基础是：地震属性携带有地下地层信息，同时地震属性和储层的含油气性之间也必然存在
某种形式的内在联系。
107.在步骤220中，对所述叠后地震数据分别进行多个预设方向的曲波变换，以获得多个预设方向下的地震数据体。
108.这里，由于裂缝的发育方向以及发育密度具有随机性，同一个岩层的裂缝方向发育也会存在随机性。通过对所述叠后地震数据分别进行多个预设方向的曲波变换，可以获得全方位描述裂缝发育密度以及发育方向的数据。
109.其中，曲波变换是继小波变换、脊波变换等发展起来的一种能够对高维信号提供近乎最优稀疏表示的数学变换方法，具有多分辨率特性，时频局部性、多方向性和各向异性，而且克服了小波变换表示边缘、轮廓等高维奇异时存在的局限性。曲波变换作为一种多尺度多方向分析方法，依据曲波变换角度、方向、位置等特征，对信号进行分析，能够有效保留地震信号中的有效信号。
110.在步骤230中，分别对每个预设方向下的地震数据体进行地震相干处理，以获得每个预设方向下的相干体。
111.这里，相干体是指由常规地震数据体经过相干处理而得到的一个新的数据体，通过对各个方向下的地震数据体进行相干处理，获得各个方向下的地震数据体对应的相干体。
112.其中，相干体技术是一项重要的地震属性技术，它通过计算相邻地震道波形的相似性将三维地震数据体转换为相干数据体，突出了波形的不连续性特征。因此，相干体能够度量由于构造、地层、岩性、油气等因素的变化所引起的地震响应的横向变化，从而有效揭示断层、裂缝、岩性体边缘和不整合等地质现象，反映地质异常特征的平面展布。
113.在本实施方式中，可以利用第三代相干算法对每个预设方向下的地震数据体进行地震相干处理，以分别对经过不同预设方向曲波变换得到的地震数据进行地震相干处理，从而获得每个预设方向下的相干体。
114.在步骤240中，根据各个预设方向下的相干体，获得该暗河区域的裂缝发育密度数据体以及裂缝发育方向数据体。
115.这里，该暗河区域的裂缝发育密度数据体指的是在裂缝在该暗河区域中的发育密度数据，该裂缝发育方向数据体则指的是在裂缝在该暗河区域中的发育方向数据。
116.在步骤250中，从所述裂缝发育密度数据体中提取目标岩层的蚂蚁体属性，获得该目标岩层的裂缝密度蚂蚁体。
117.这里，蚂蚁追踪技术是一种模仿蚁群寻找食物的仿生算法，其能够突出地震数据的不连续性，是一种强化断裂特征的新属性，能够快速了解区域内断层的发育和展布情况，能够提高断层的解释进度，使地质构造的细节尽可能全面地展示出来。
118.其基本原理是在预处理地震数据体的基础上产生地震蚂蚁属性体，自动提取一套详细地反映不连续断片的地震属性体。在地震属性体上撒播蚂蚁种子点、定义蚂蚁觅食线路偏移度、蚂蚁搜索步长、法定和非法定范围、搜索终止门槛值，最终得到蚂蚁属性体
119.其中，基于多方向曲波变换获得的裂缝发育密度平面属性以及裂缝发育方向平面属性反映了暗河小尺度裂缝的发育情况。而通过对裂缝发育密度数据体进行线性增强，提取蚂蚁体属性，并调整不同方向的参数，获得裂缝密度蚂蚁体，该裂缝密度蚂蚁体能够反映暗河区域中大尺度裂缝的发育情况，从而为开发人员在暗河储层设置钻井提供数据支持。
120.在步骤260中，将所述裂缝密度蚂蚁体与该暗河区域对应的暗河展布图进行叠合显示，获得第二叠合展布图，以根据所述第二叠合展布图分析大尺度裂缝在暗河区域的发育情况。
121.这里，通过将能够反映暗河中大尺度裂缝的发育情况，裂缝密度蚂蚁体与该暗河区域对应的暗河展布图进行叠合显示，获得第二叠合展布图，使得开发人员能够根据该第二叠合展布图直观地了解大尺度裂缝在暗河区域中的发育情况，从而为开发人员如何设置钻井提供数据支持。
122.图4示出了本发明实施例二提出的第二叠合展布图的示意图，如图4所示，可以看出暗河的走向与大尺度裂缝的展布方向具有相关性。在图4中，大尺度裂缝的发育方向主要有两个方向，第一个发育方向是北西方向，第二个发育方向是北东方向。从图4中可以看出沿两个裂缝发育方向的暗河展布方向，说明大尺度裂缝对暗河的改向具有显著影响。
123.实施例三
124.在上述实施例的基础上，本发明的实施例二还可以提供一种暗河裂缝发育预测方法。图5示出了本发明实施例三提出的一种暗河裂缝发育预测方法的流程示意图，如图5所示，该暗河裂缝发育预测方法可以包括：步骤310至步骤390。
125.在步骤310中，获取暗河区域的叠后地震数据。
126.这里，利用采集设备对地震波信号进行采集得到地震数据体，然后对该地震数据体进行叠加处理得到叠后地震数据。通过对叠后地震数据进行处理，能够从中提取出反映地震属性的数据，该地震属性是有关地震波的几何形态、运动学特征、动力学特征和统计学特性的特殊度量。
127.其中，地震波信号在地下地层传播过程中受地层岩性、物性等因素影响会产生相应变化，它是对地下储层综合特征的一种复杂反映。地下地层岩石物理等性质的空间变化，必然导致地震反射波特征的变化，进而影响到地震属性。特别是当储层含油气时，其地震响应特征会发生相应变化，相对应的地震属性也会有所体现。地震属性技术预测油气的理论基础是：地震属性携带有地下地层信息，同时地震属性和储层的含油气性之间也必然存在某种形式的内在联系。
128.在步骤320中，对所述叠后地震数据分别进行多个预设方向的曲波变换，以获得多个预设方向下的地震数据体。
129.这里，由于裂缝的发育方向以及发育密度具有随机性，同一个岩层的裂缝方向发育也会存在随机性。通过对所述叠后地震数据分别进行多个预设方向的曲波变换，可以获得全方位描述裂缝发育密度以及发育方向的数据。
130.其中，曲波变换是继小波变换、脊波变换等发展起来的一种能够对高维信号提供近乎最优稀疏表示的数学变换方法，具有多分辨率特性，时频局部性、多方向性和各向异性，而且克服了小波变换表示边缘、轮廓等高维奇异时存在的局限性。曲波变换作为一种多尺度多方向分析方法，依据曲波变换角度、方向、位置等特征，对信号进行分析，能够有效保留地震信号中的有效信号。
131.值得说明的是，关于如何对叠后地震数据进行曲波变换的操作已在上述实施例中进行了说明，在此不再赘述。
132.在步骤330中，分别对每个预设方向下的地震数据体进行地震相干处理，以获得每
个预设方向下的相干体。
133.这里，相干体是指由常规地震数据体经过相干处理而得到的一个新的数据体，通过对各个方向下的地震数据体进行相干处理，获得各个方向下的地震数据体对应的相干体。
134.其中，相干体技术是一项重要的地震属性技术，它通过计算相邻地震道波形的相似性将三维地震数据体转换为相干数据体，突出了波形的不连续性特征。因此，相干体能够度量由于构造、地层、岩性、油气等因素的变化所引起的地震响应的横向变化，从而有效揭示断层、裂缝、岩性体边缘和不整合等地质现象，反映地质异常特征的平面展布。
135.优选地，可以利用第三代相干算法对每个预设方向下的地震数据体进行地震相干处理，以分别对经过不同预设方向曲波变换得到的地震数据进行地震相干处理，从而获得每个预设方向下的相干体。
136.在步骤340中，根据各个预设方向下的相干体，获得该暗河区域的裂缝发育密度数据体以及裂缝发育方向数据体。
137.这里，该暗河区域的裂缝发育密度数据体指的是在裂缝在该暗河区域中的发育密度数据，该裂缝发育方向数据体则指的是在裂缝在该暗河区域中的发育方向数据。
138.值得说明的是，关于如何根据各个预设方向下的相干体，获得该暗河区域的裂缝发育密度数据体以及裂缝发育方向数据体的详细步骤已在上述实施例中进行了详细说明，在此不再赘述。
139.在步骤350中，根据所述裂缝发育密度数据体，获取目标岩层对应的裂缝发育密度平面属性。
140.这里，目标岩层是指将要研究的储层，即在采油或采气时，油气可能会蕴藏的岩层。
141.其中，裂缝发育密度数据体中包括多个不同深度的岩层的裂缝发育密度，从裂缝发育密度数据体中获取目标岩层对应的裂缝发育密度平面属性，从而根据该裂缝发育密度平面属性对目标岩层所在平面的裂缝密度进行分析。该裂缝发育密度平面属性可以是目标岩层所在平面的裂缝密度展布图，用于反映裂缝在目标岩层的发育密度。
142.值得说明的是，该裂缝发育密度平面属性反映了裂缝在目标岩层的发育密度，即目标岩层的裂缝密度。裂缝密度是表示裂缝发育程度的一个概念值。由于计算时使用的资料不同,又有线密度、面密度、体积密度之分,但在实际使用中并不那么严密,各种计算方法都是为了反映裂缝发育程度。据这个概念值,可将研究区划分出裂缝发育带、次发育带及不发育带。
143.在步骤360中，根据所述裂缝发育方向数据体，获取目标岩层对应的裂缝发育方向平面属性。
144.这里，裂缝发育密度数据体中包括多个不同深度的岩层的裂缝发育方向平面属性，从裂缝发育密度数据体中获取目标岩层对应的裂缝发育密度平面属性，从而根据该裂缝发育方向平面属性对目标岩层所在平面的裂缝方向进行分析。该裂缝发育密度平面属性可以是目标岩层所在平面的裂缝方向展布图，用于反映目标岩层的裂缝的发育方向。
145.应当理解的是，裂缝性气藏裂缝是天然气主要的渗流通道，裂缝的发育方向和发育密度决定了气井的产能。通过对目标岩层的裂缝发育密度以及裂缝发育方向进行分析，
可以得到目标岩层的裂缝发育情况，从而为后续探井的设置提供依据。
146.其中，该裂缝发育方向平面属性反映了裂缝在目标岩层的发育方向，即目标岩层的裂缝发育方向。
147.值得说明的是，步骤350与步骤360的执行可以没有先后顺序，可以同时执行步骤350与步骤360。
148.在步骤370中，将所述裂缝发育密度平面属性和所述裂缝发育方向平面属性与该暗河区域对应的暗河展布图进行叠合显示，获得第一叠合展布图，以根据所述第一叠合展布图分析小尺度裂缝在暗河区域的发育情况。
149.这里，裂缝裂缝按尺度进行划分，可分为微裂缝、小尺度裂缝以及大尺度裂缝。
150.其中，微裂缝是指在岩心薄片上能观察到的裂缝，一般是微米级的裂缝，微裂缝在储层中数量庞大，通常是与基质相伴生的，因此微裂缝实际是基质的一部分。
151.小尺度裂缝一般是指在岩心和成像测井上能观察到的裂缝，长度一般在几十厘米到几十米。成像测井和岩心上看到的是小尺度裂缝的高度。小尺度裂缝在三维空间交织成网，形成油气的渗流通道。裂缝密度、裂缝开度、裂缝长度和裂缝方位变化程度是决定裂缝渗透率的关键因素。
152.大尺度裂缝一般是指在地震级别的裂缝，长度一般在几十米到上千米。大尺度裂缝横向延伸长度大，垂向切割层位深，开度大的大尺度裂缝有很高的渗透率率，在油藏中小尺度裂缝提供了油藏的渗透率，而大尺度裂缝决定了油藏的非均质性。
153.其中，在步骤370中，对于小尺度裂缝的地震预测可以通过裂缝发育密度平面属性和裂缝发育方向平面属性来进行。通过将所述裂缝发育密度平面属性和所述裂缝发育方向平面属性与该暗河区域对应的暗河展布图进行叠合显示，得到第一叠合展布图，可以根据该第一叠合展布图来分析小尺度裂缝。
154.在步骤380中，从所述裂缝发育密度数据体中提取目标岩层的蚂蚁体属性，获得该目标岩层的裂缝密度蚂蚁体。
155.这里，这里，蚂蚁追踪技术是一种模仿蚁群寻找食物的仿生算法，其能够突出地震数据的不连续性，是一种强化断裂特征的新属性，能够快速了解区域内断层的发育和展布情况，能够提高断层的解释进度，使地质构造的细节尽可能全面地展示出来。
156.其基本原理是在预处理地震数据体的基础上产生地震蚂蚁属性体，自动提取一套详细地反映不连续断片的地震属性体。在地震属性体上撒播蚂蚁种子点、定义蚂蚁觅食线路偏移度、蚂蚁搜索步长、法定和非法定范围、搜索终止门槛值，最终得到蚂蚁属性体
157.其中，基于多方向曲波变换获得的裂缝发育密度平面属性反映了暗河小尺度裂缝的发育情况。而通过对裂缝发育密度数据体进行线性增强，提取蚂蚁体属性，并调整不同方向的参数，获得裂缝密度蚂蚁体，该裂缝密度蚂蚁体能够反映暗河区域中大尺度裂缝的发育情况，从而为开发人员在暗河储层设置钻井提供数据支持。
158.在步骤390中，将所述裂缝密度蚂蚁体与所述第一叠合展布图进行叠合显示，获得第三叠合展布图，以根据所述第三叠合展布图分析小尺度裂缝以及大尺度裂缝在暗河区域的发育情况。
159.这里，通过将所述裂缝密度蚂蚁体与所述第一叠合展布图进行叠合显示，获得第三叠合展布图，以根据所述第三叠合展布图分析小尺度裂缝以及大尺度裂缝在暗河区域的
发育情况，可以根据该第三叠合展布图分析不同尺度的裂缝在暗河区域的发育情况。其中，基于多方向曲波变换获得的裂缝发育密度平面属性以及裂缝发育方向平面属性反映了暗河小尺度裂缝的发育情况。而通过对裂缝发育密度数据体进行线性增强，提取蚂蚁体属性，并调整不同方向的参数，获得裂缝密度蚂蚁体，该裂缝密度蚂蚁体能够反映暗河区域中大尺度裂缝的发育情况。因此，通过该第三叠合展布图可以展示出不同尺度裂缝与暗河发育的关系。
160.图6示出了本发明实施例三提出裂缝发育密度平面属性、裂缝发育方向平面属性以及裂缝密度蚂蚁体叠合显示的叠合图的示意图，如图6所示，可以看出，在暗河发育的边部大尺度裂缝发育情况，暗河的走向与大尺度裂缝的展布方向关系密切。这与暗河发育的边部应力释放有关，地表水及地下水沿裂缝发育带进入碳酸盐岩储层，形成管道状地下暗河，从而影响暗河的发育情况。由此，通过将该第三叠合展布图，可以直观地分析出小尺度裂缝、大尺度裂缝在暗河区域中的发育情况。
161.在一个可选的实施方式中，所述方法还可以包括：
162.获取所述暗河区域的溶洞分布图；
163.将所述溶洞分布图与所述第三叠合展布图进行叠合显示，获得第四叠合展布图，以根据所述第四叠合展布图分析小尺度裂缝以及大尺度裂缝在暗河区域的溶洞中的发育情况。
164.这里，所述暗河区域的溶洞分布图是指该暗河区域内的溶洞展布图。其中，该溶洞分布图可以通过地震资料获得，其为现有技术，在此不多赘述。
165.其中，通过将所述溶洞分布图与所述第三叠合展布图进行叠合显示，获得第四叠合展布图，以根据所述第四叠合展布图分析小尺度裂缝以及大尺度裂缝在暗河区域的溶洞中的发育情况，是要在一个图像中综合展示出裂缝、暗河以及溶洞的分布状况，从而综合分析裂缝、溶洞以及暗河三者之间的影响。
166.图7示出了本发明实施例三提出的第四叠合展布图的示意图，如图7所示，可以看出，在暗河的周边落水洞的发育情况，由于暗河的周边地应力释放，导致溶洞周边小尺度裂缝发育密度大，在图中往往呈现圆圈状。因此，在该第四叠合展布图中不仅可以看出小尺度裂缝、大尺度裂缝对暗河发育的影响，而且可以直观地看出溶洞周边不同尺度的裂缝发育情况。
167.图8示出了某工区中测井的裂缝地震预测示意图，如图8所示，其中，图a为某工区暗河区域的第四叠合展布图，从图a中可以看出，小尺度裂缝发育方向为nne向。图b为某工区中测井获得的测井成像结果，从该井的测井解释及成像测井结果可以看到，其展示了该井区水平最大主应力方位为nne-ssw向，约22-37
°
或202-217
°
，井段顶部5380-5428m各向异性响应明显。可见，第四叠合展布图的地震预测结果与测井成像解释结果相吻合，说明了基于地震进行裂缝预测方法的可靠性。
168.实施例四
169.根据本发明的实施例四，还可以提供一种暗河裂缝发育预测系统，该系统包括：
170.获取单元，用于获取暗河区域的叠后地震数据；
171.曲波变换模块，用于对所述叠后地震数据分别进行多个预设方向的曲波变换，以获得多个预设方向下的地震数据体；
172.相干处理模块，用于分别对每个预设方向下的地震数据体进行地震相干处理，以获得每个预设方向下的相干体；
173.裂缝数据获取模块，用于根据所述裂缝发育密度数据体，获取目标岩层对应的裂缝发育密度平面属性；以及
174.用于根据所述裂缝发育方向数据体，获取目标岩层对应的裂缝发育方向平面属性；
175.叠合显示模块，用于将所述裂缝发育密度平面属性和所述裂缝发育方向平面属性与该暗河区域对应的暗河展布图进行叠合显示，获得第一叠合展布图，以根据所述第一叠合展布图分析小尺度裂缝在暗河区域的发育情况。
176.可选地，所述系统还包括：
177.蚂蚁体提取模块，用于从所述裂缝发育密度数据体中提取目标岩层的蚂蚁体属性，获得该目标岩层的裂缝密度蚂蚁体；
178.所述叠合显示模块还用于将所述裂缝密度蚂蚁体与该暗河区域对应的暗河展布图进行叠合显示，获得第二叠合展布图，以根据所述第二叠合展布图分析大尺度裂缝在暗河区域的发育情况。
179.可选地，所述叠合显示模块还用于将所述裂缝密度蚂蚁体与所述第一叠合展布图进行叠合显示，获得第三叠合展布图，以根据所述第三叠合展布图分析小尺度裂缝以及大尺度裂缝在暗河区域的发育情况。
180.可选地，所述获取模块还用于获取所述暗河区域的溶洞分布图；
181.所述叠合显示模块还用于将所述溶洞分布图与所述第三叠合展布图进行叠合显示，获得第四叠合展布图，以根据所述第四叠合展布图分析小尺度裂缝以及大尺度裂缝在暗河区域的溶洞中的发育情况。
182.可选地，所述裂缝数据获取模块具体用于从所述裂缝发育密度数据体中提取目标岩层上与裂缝相关的沿层地震属性，以根据所述沿层地震属性获得所述目标岩层对应的裂缝发育密度平面属性；以及
183.从所述裂缝发育方向数据体中提取目标岩层上与裂缝相关的沿层地震属性，根据所述沿层地震属性获得所述目标岩层对应的裂缝发育方向平面属性。
184.可选地，所述预设方向包括0
°
、45
°
、90
°
以及135
°
中的至少一个方向。
185.实施例五
186.根据本发明的实施例，还提供了一种存储介质，所述存储介质上存储有程序代码，所述程序代码被处理器执行时，实现如上述实施例任一项所述的暗河裂缝发育预测方法。
187.实施例六
188.根据本发明的实施例，还提供了一种电子设备，所述电子设备包括存储器、处理器，所述存储器上存储有可在所述处理器上运行的程序代码，所述程序代码被所述处理器执行时，实现如上述实施例任一项所述的暗河裂缝发育预测方法。
189.以上结合附图详细说明了本发明的技术方案，考虑到相关技术中，基于不同尺度的裂缝发育对古暗河储层具有不同程度的影响和改善作用，而现有技术中，基本通过测井来测定区域的裂缝发育情况，而通过测井的方式只能反映部分区域的裂缝发育情况，无法对整个暗河区域的裂缝发育情况进行分析。本发明提供一种暗河裂缝发育预测方法、系统、
存储介质以及电子设备，通过对暗河区域的叠后地震数据进行多个预设方向的曲波变换，从而获得多个预设方向下的地震数据体；进而分别对多个预设方向下的地震数据体进行地震相干处理，从而获得该暗河区域的裂缝发育密度数据体以及裂缝发育方向数据体，并且从裂缝发育密度数据体以及裂缝发育方向数据体中分别提取目标岩层对应的裂缝发育密度平面属性以及目标岩层对应的裂缝发育方向平面属性，从而利用目标岩层对应的裂缝发育密度平面属性以及裂缝发育方向平面属性来分析小尺度裂缝在暗河区域的分布状况。通过将目标岩层对应的裂缝发育密度平面属性以及裂缝发育方向平面属性与暗河区域的暗河展布图进行叠合显示，可以直观地了解小尺度裂缝在暗河各个区域的发育情况，从而为后续钻井设置提供数据依据。
190.在本技术所提供的几个实施例中，应该理解到，所揭露的装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。
191.作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本发明实施例方案的目的。
192.另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以是两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
193.集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分，或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台电子设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read-only memory)、随机存取存储器(ram，random access memory)、磁碟或者光盘等各种可以存储程序代码的介质。
194.虽然本发明所公开的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所公开的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的保护范围，仍须以所附的权利要求书所界定的范围为准。