裂缝动态识别方法、装置、存储介质和计算机设备与流程

1.本发明涉及油藏开发技术领域，尤其涉及一种裂缝动态识别方法、装置、存储介质和计算机设备。


背景技术：

2.在缝洞型油气藏中，裂缝不仅是碳酸盐岩油气藏的流动通道，也是最主要的储集空间类型之一。因此，尽可能的识别储层中的裂缝分布对碳酸盐岩油气资源的勘探开发至关重要。
3.与其他类型油藏一样，缝洞油藏的裂缝识别成果可以用于油藏描述和油藏开发的方方面面。其裂缝识别成果不仅可以直接影响油藏地质模型的准确性和油藏数值模拟计算的可靠性，同时也是为油井制定合理的工作制度及开发调整对策的主要依据。但是，定量化裂缝识别目前仍是国际上公认的难题，也是油藏描述的热点问题之一。
4.常规的裂缝识别和描述手段很多，例如，包括岩心观察法、成像测井法、地震反演法、示踪剂监测法、吸水剖面分析法等。其中，岩心观察法和成像测井法仅限于观察井周的裂缝情况；地震反演法不仅识别精度低，并且难以识别判断裂缝导致的油井见水是底水还是注入水。这些方法都属于静态识别，无法识别不同开发阶段之间的裂缝的变化情况，而且，这些方法的内容和结论多数属于定性识别，定量程度低，这直接影响了裂缝识别成果的应用效果。
5.塔河油田和顺北地区断溶体油藏储层的埋藏深度多数达到5000米以上，最深可达8000-9000米左右。例如顺北蓬1井设计井深8593米。断溶体油藏研究结果表明，在油藏深部发育着大量的垂直裂缝，这些垂直裂缝与底水沟通，由于底水能量高体积大，沿着高角度裂缝水窜严重，造成油井快速见水，直接影响了断溶体油藏开发效果。在断溶体油藏开发的中后期，为了保持地层能量，需要将部分开发效果差、高含水油井转注，这样就产生了注水沿横向裂缝向油井指进的情况。
6.底水和注入水沿着裂缝窜流会导致油井快速水淹，如何识别注水和底水造成的水淹情况，并针对不同来水分别进行治理，以提高水淹治理方法的针对性和有效性是一项迫切的工作。大量的油田开发实践经验证明，明确的裂缝的动态发育特征可以为油井水淹治理的方法提供直接的依据。因此，裂缝发育的动态识别具有很高的应用价值。从生产过程看，裂缝对油水运动规律、含水上升规律、水驱采收率等都有很大影响，甚至可以是制约性的因素。不同方位的裂缝发育程度不同，含水变化过程也不同，相应的治理措施也不同。只有确立了不同方位的裂缝发育情况，并对生产过程中的裂缝动态变化规律进行分析，才能指导地质建模、方案编制和治理措施等。
7.从生产和科研分析的需求出发，裂缝识别是缝洞油藏描述的核心工作之一。然而，现有的裂缝识别和描述的方法是一个世界级的难题，特别是底水发育的断溶体油藏，开发时间短，经验少，裂缝识别存在大量关键性问题，包括裂缝动态识别的方法、裂缝发育程度的差异性分析、识别结果的验证、识别方法的量化表征等等，这些问题在不同程度上都给缝
洞型油藏储层的识别和开发带来了很大困扰。因此，亟需一种新的动态识别方法来识别裂缝。


技术实现要素：

8.本发明的主要目的是提供一种裂缝动态识别方法、装置、存储介质和计算机设备，以实现对油藏缝洞发育程度的量化识别。
9.第一方面，本技术提供一种裂缝动态识别方法，包括以下步骤：将目标井在不同生产阶段的生产动态数据从时域变换到频域，以得到目标井在不同生产阶段的生产谱图，其中，所述生产动态数据为表征目标井生产能力的一种参数的动态数据；对于每个生产阶段的生产谱图，对生产谱图进行量化分析，根据分析结果确定该生产阶段目标井下方的缝洞类型和各个缝洞类型之间的体积比；对于每个生产阶段，根据目标井在该生产阶段的总产量和目标井下方的缝洞类型及各个缝洞类型之间的体积比，确定在该生产阶段目标井下方各个缝洞类型的体积，从而得到目标井在各个生产阶段的各个缝洞类型的体积。
10.在一个实施例中，所述生产动态数据包括产量的动态数据、压力的动态数据或含水的动态数据；所述对生产谱图进行量化分析包括对生产谱图的谱宽、密度、振幅、功率或相位进行量化分析。
11.在一个实施例中，所述生产动态数据包括含水的动态数据，所述生产谱图为含水率谱图；当所述生产谱图为含水率谱图时，对生产谱图的振幅进行量化分析，根据分析结果确定该生产阶段目标井下方的缝洞类型和各个缝洞类型之间的体积比，包括：分析目标井的含水率谱图中各个不同频率段的平均振幅，对于所述含水率谱图的每个频率段，根据其平均振幅确定该频率段所对应的缝洞类型；分析所述含水率谱图中各个不同频率段的平均振幅之间的比例关系，根据所述比例关系，基于所述含水率谱图中每个频率段所对应的缝洞类型，确定目标井下方各个缝洞类型之间的体积比。
12.在一个实施例中，当目标井为多井缝洞单元的生产井时，将目标井在不同生产阶段的生产动态数据从时域变换到频率域，以得到目标井在各个生产阶段的生产谱图，包括：获取目标井在各个生产阶段注水前后的生产动态数据；将目标井在各个生产阶段注水前后的生产动态数据从时域变换到频域，以得到目标井在各个生产阶段注水前后的生产谱图；对于每个生产阶段，确定目标井在注水后的生产谱图减其在注水前的生产谱图所得的差值谱图，将所述差值谱图作为目标井在该生产阶段的生产谱图，从而得到目标井在各个生产阶段的生产谱图。
13.在一个实施例中，在获取目标井注水前的生产动态数据之后且在获取目标井注水后的生产动态数据之前，所述方法还包括：对多井缝洞单元进行连通性分析，确定与目标井连通的注水井，并向所述注水井注水；在目标井因所述注水井注水而受效后，获取目标井在该生产阶段的生产动态数据。
14.在一个实施例中，目标井的含水率谱图中各个不同频率段包括频率为0～0.02、0.02～0.1以及频率大于0.1的三个频率段；频率范围为0～0.02时所对应的缝洞类型为大型溶洞，频率范围为0.02～0.1时所对应的缝洞类型为大型裂缝，频率范围为大于0.1时所对应的缝洞类型为小型裂缝。
15.在一个实施例中，目标井下方的缝洞类型包括大型溶洞体、大尺度裂缝和中小尺
度裂缝，其中，所述大型溶洞体为直径大于500mm的溶洞，所述大尺度裂缝为裂缝宽度大于1mm的裂缝，所述中小尺度裂缝为裂缝宽度小于或等于1mm的裂缝。
16.第二方面，本技术提供一种裂缝动态识别装置，包括：数据变换模块，用于将目标井在不同生产阶段的生产动态数据从时域变换到频域，以得到目标井在不同生产阶段的生产谱图，其中，所述生产动态数据为表征目标井生产能力的一种参数的动态数据；数据分析模块，用于对于每个生产阶段的生产谱图，对生产谱图进行量化分析，根据分析结果确定该生产阶段目标井下方的缝洞类型和各个缝洞类型之间的体积比；体积计算模块，用于对于每个生产阶段，根据目标井在该生产阶段的总产量和目标井下方的缝洞类型及各个缝洞类型之间的体积比，确定在该生产阶段目标井下方各个缝洞类型的体积，从而得到目标井在各个生产阶段的各个缝洞类型的体积。
17.第三方面，本技术提供一种存储介质，存储有计算机程序，所述计算机程序被处理器执行时，实现如上文所述的裂缝动态识别方法的步骤。
18.第四方面，本技术提供一种计算机设备，包括处理器和存储有程序代码的存储介质，所述程序代码被所述处理器执行时，实现如上文所述的裂缝动态识别方法的步骤。
19.本发明提供的裂缝动态识别方法，将油田生产动态数据进行时频域变换处理得到频率域的生产谱图，对所得到的生产谱图进行分析，来实现对缝洞型油藏中的缝洞发育程度的量化分析，实现缝洞油藏动态差异化、定量化的识别，为缝洞油藏描述、地质建模和油井水淹治理提供依据。本发明的方法简单、可靠，且精确度高，所需资料均基于现场生产，来源广泛，获取资料的成本较低，该方法不仅适用于碳酸盐岩储层，尤其适用于断溶体油藏，而且适用于其他底水型油气藏描述。
附图说明
20.构成本技术的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定，在附图中：
21.图1为根据本技术一示例性实施方式的裂缝动态识别方法的流程图；
22.图2为根据本技术一具体实施例的裂缝动态识别方法的流程图；
23.图3为根据本技术一具体实施例的th10433h井开采曲线；
24.图4为根据本技术一具体实施例的th10433h井含水率变化曲线；
25.图5为根据本技术一具体实施例的th10433井的含水率变化频谱；
26.图6为根据本技术一具体实施例的th10433h井只包含裂缝时的含水率变化频谱；
27.图7为根据本技术一具体实施例的th10433h井地震处理结果；
28.图8为根据本技术一具体实施例的s86井开采曲线；
29.图9为根据本技术一具体实施例的s86井注水受效前的含水率频谱；
30.图10为根据本技术一具体实施例的s86井只包含裂缝时的注水受效前的含水率频谱；
31.图11为根据本技术一具体实施例的s86井注水受效后的含水率频谱；
32.图12为根据本技术一具体实施例的s86井地震处理结果。
具体实施方式
33.需要说明的是，在不冲突的情况下，本技术中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
34.实施例一
35.本实施例提供一种裂缝动态识别方法，图1为根据本技术一示例性实施方式的裂缝动态识别方法的流程图。如图1所示，裂缝动态识别方法可以包括以下步骤：
36.s100：将目标井在不同生产阶段的生产动态数据从时域变换到频域，以得到目标井在不同生产阶段的生产谱图，其中，所述生产动态数据为表征目标井生产能力的一种参数的动态数据；
37.s200：对于每个生产阶段的生产谱图，对生产谱图进行量化分析，根据分析结果确定该生产阶段目标井下方的缝洞类型和各个缝洞类型之间的体积比；
38.s300：对于每个生产阶段，根据目标井在该生产阶段的总产量和目标井下方的缝洞类型及各个缝洞类型之间的体积比，确定在该生产阶段目标井下方各个缝洞类型的体积，从而得到目标井在各个生产阶段的各个缝洞类型的体积。
39.实施例二
40.本实施例提供一种裂缝动态识别方法，包括以下步骤：
41.第一步，将目标井在不同生产阶段的生产动态数据从时域变换到频域，以得到目标井在不同生产阶段的生产谱图，其中，所述生产动态数据为表征目标井生产能力的一种参数的动态数据，例如可以包括产量的动态数据、压力的动态数据或含水的动态数据。
42.当目标井为单井缝洞单元的生产井时，将目标井在不同生产阶段的生产动态数据从时域变换到频率域，以得到目标井在各个生产阶段的生产谱图，包括：获取目标井在各个生产阶段注水前的生产动态数据；将目标井在各个生产阶段注水前的生产动态数据从时域变换到频域，以得到目标井在各个生产阶段注水前的生产谱图；对于每个生产阶段，将目标井在注水前的生产谱图作为目标井在该生产阶段的生产谱图，从而得到目标井在各个生产阶段的生产谱图。
43.当目标井为多井缝洞单元的生产井时，将目标井在不同生产阶段的生产动态数据从时域变换到频率域，以得到目标井在各个生产阶段的生产谱图，包括：获取目标井在各个生产阶段注水前后的生产动态数据；将目标井在各个生产阶段注水前后的生产动态数据从时域变换到频域，以得到目标井在各个生产阶段注水前后的生产谱图；对于每个生产阶段，确定目标井在注水后的生产谱图减其在注水前的生产谱图所得的差值谱图，将所述差值谱图作为目标井在该生产阶段的生产谱图，从而得到目标井在各个生产阶段的生产谱图。
44.其中，当目标井为多井缝洞单元的生产井时，在获取目标井注水前的生产动态数据之后且在获取目标井注水后的生产动态数据之前，所述方法还包括：对多井缝洞单元进行连通性分析，确定与目标井连通的注水井，并向所述注水井注水；在目标井因所述注水井注水而受效后，获取目标井在该生产阶段的生产动态数据。
45.可以通过多种方法将目标井在不同生产阶段的生产动态数据从时域变换到频域。例如可以通过傅里叶变换的方法。
46.作为一种数学处理方法和手段，傅里叶变换在数字信号处理领域是一种非常有效的分析方法。在复杂信号处理中，傅里叶变换的典型作用是将信号分解成频率谱——显示
与频率对应的幅值大小。它可以将时间域的信号转换到频率域进行分析(正变换)，也可以将频率域的信号转换为时间域的信号进行分析(逆变换)。
47.其中，傅立叶正变换的表达式为：
[0048][0049]
傅里叶逆变换的表达式为：
[0050][0051]
将钻遇缝洞储集体油井的生产动态数据看成是油藏生产信号的反映，对其进行傅里叶变换，可以得到不同频率的分布，分析不同开发阶段动态指标的频率变化，可以得到裂缝和溶洞的动态指标。
[0052]
第二步，对于每个生产阶段的生产谱图，对生产谱图进行量化分析，根据分析结果确定该生产阶段目标井下方的缝洞类型和各个缝洞类型之间的体积比。
[0053]
其中，对生产谱图进行量化分析包括对生产谱图的谱宽、密度、振幅、功率或相位进行量化分析。当生产动态数据包括含水的动态数据时，生产谱图为含水率谱图。当生产谱图为含水率谱图时，对生产谱图的振幅进行量化分析，根据分析结果确定该生产阶段目标井下方的缝洞类型和各个缝洞类型之间的体积比，包括：分析目标井的含水率谱图中各个不同频率段的平均振幅，对于含水率谱图的每个频率段，根据其平均振幅确定该频率段所对应的缝洞类型；分析所述含水率谱图中各个不同频率段的平均振幅之间的比例关系，根据所述比例关系，基于所述含水率谱图中每个频率段所对应的缝洞类型，确定目标井下方各个缝洞类型之间的体积比。
[0054]
目标井下方的缝洞类型可以包括大型溶洞体、大尺度裂缝和中小尺度裂缝，其中，所述大型溶洞体为直径大于500mm的溶洞，所述大尺度裂缝为裂缝宽度大于1mm的裂缝，所述中小尺度裂缝为裂缝宽度小于或等于1mm的裂缝。
[0055]
目标井的含水率谱图中各个不同频率段可以包括频率为0～0.02、0.02～0.1以及频率大于0.1的三个频率段，其中，频率范围为0～0.02时所对应的缝洞类型为大型溶洞，频率范围为0.02～0.1时所对应的缝洞类型为大型裂缝，频率范围为大于0.1时所对应的缝洞类型为小型裂缝。
[0056]
第三步，对于每个生产阶段，根据目标井在该生产阶段的总产量和目标井下方的缝洞类型及各个缝洞类型之间的体积比，确定在该生产阶段目标井下方各个缝洞类型的体积，从而得到目标井在各个生产阶段的各个缝洞类型的体积。
[0057]
其中，总产量可以包括总产油量、总产气量等。
[0058]
本发明提供的裂缝动态识别方法，将油田生产动态数据进行时频域变换处理得到频率域的生产谱图，对所得到的生产谱图进行分析，来实现对缝洞型油藏中的缝洞发育程度的量化分析，实现缝洞油藏动态差异化、定量化的识别，为缝洞油藏描述、地质建模和油井水淹治理提供依据。本发明的方法简单、可靠，且精确度高，所需资料均基于现场生产，来源广泛，获取资料的成本较低，该方法不仅适用于碳酸盐岩储层，尤其适用于断溶体油藏，而且适用于其他底水型油气藏描述。
[0059]
实施例三
[0060]
本实施例通过一个具体例子来对本技术的裂缝动态识别方法进行介绍。
[0061]
通过数值模拟计算发现，底水发育的缝洞油藏其开发特征与砂岩油藏完全不同。在底水发育的缝洞油藏中，底水和注入水沿裂缝快速推进，造成油井严重水淹，而且，裂缝开度随见水时间和见水强度不同也差异很大，裂缝一旦见水，裂缝将很少产油。
[0062]
通过生产实践发现，油井通过裂缝与底水和注水井相互沟通，在压力差的作用下，不同裂缝产水发生时间和强度变化呈不规则的谱图分布，裂缝发育均匀，谱宽越大，反之裂缝发育程度越集中，谱宽越窄，谱峰越高。谱图高度即振幅大小反映了裂缝的体积大小。因此，可以说裂缝见水分布频谱图直接反映了断溶体油藏缝洞储集体的发育程度。利用统计分析，借助频谱分析方法，对大量的谱图进行量化分析，可以得到不同储集体类型的截止值，以此为依据可以建立缝洞发育程度定量识别和描述。
[0063]
按照地球物理基础理论和后期缝洞油藏地震属性的研究成果，地震波的分频信号与裂缝大小有直接的关系，裂缝越大，地震反射波的频率越低，反之亦然。在油藏地震反演过程中，傅立叶变换被引进处理地震反射波信号，得到不同的反射波频率分布，根据不同反射波频率分布得到的地震成像可以用于识别不同尺度的缝洞结构，如孔洞型储层的“串珠状”地震反射特征。因此，地震波频率的大小反应了裂缝的大小分布。基于此原理可以将油井生产数据视为油藏系统输出的生产信号，对生产信号进行数学处理，得到生产频谱，为后期的裂缝量化描述提供条件。
[0064]
生产数据的动态属性，决定了基于生产数据频谱分析开展的裂缝识别是一种动态的识别方法。在该实施例中，我们从数学处理方法入手，以快速傅里叶变换为手段，对油井含水率的变化特征进行计算，实现裂缝的动态识别，为断溶体油藏地质建模和油井水淹治理提供依据。
[0065]
本实施例的裂缝动态识别方法包括：
[0066]
(1)生产谱图计算：对不同开发阶段的生产动态数据进行傅里叶变换，将生产曲线由时间域转换为频率域分布，形成生产谱图，一般包括产量、压力、含水三种谱图类型。
[0067]
(2)裂缝发育程度定量分析：利用频谱分析技术，对频谱进行定量分析，得到谱宽、密度、振幅、功率、相位等量化指标，实现量化计算。
[0068]
(3)不同裂缝方位发育程度分析：对单井缝洞单元生产频谱得到的是与底水沟通的裂缝信息，这类裂缝多数属于高角度缝；对多井缝洞单元，在注水受效期间，利用频谱反映与底水和注入水沟通的缝洞，得到中低角度发育的裂缝体系。
[0069]
(4)不同类型储集空间类型发育程度定量分析：按照频率差异，区分不同储集空间类型。其中，频率0-0.02的部分对应是大型溶洞体，频率0.02-0.1的部分对应的是大型裂缝，频率大于0.1的部分对应的是中小尺度裂缝。
[0070]
下面以断溶体油藏的生产井的含水率为变化为例对本技术的技术方案进行说明，不同缝洞单元，进行裂缝动态识别的流程包括以下步骤，(其中，对单井缝洞单元执行第(1)、(2)和(6)步即可，对多井缝洞单元则需要执行步骤(1)至(6)(如图2所示))：
[0071]
(1)选取断溶体油藏地质背景的一口油井，收集该井的生产动态数据，对生产动态数据进行筛选以删除工作制度和措施变化引起的异常点；
[0072]
(2)对生产井注水进行分析，取水井注水前的含水率变化部分进行傅里叶变换处理，得到与底水沟通的高角度裂缝的频谱曲线；
[0073]
(3)进行连通性分析，得到油水井间的连通情况，确定与生产井连通的注水井进行注水，对油井受效前后的含水曲线分别进行分析，确定含水来源是底水还是注入水。
[0074]
(4)对生产井因注水受效后的含水率变化部分进行傅里叶变换处理，得到与底水及注水井沟通的裂缝分布的频谱曲线；
[0075]
(5)将注水前后生产井的两个频谱曲线进行处理，分别得到与底水沟通的高角度裂缝及与注水井沟通的中低角度裂缝的动态分布结果。
[0076]
在注水前，生产井与底水沟通，频谱曲线反映生产井下方高角度裂缝的分布状况，在注水之后，生产井与底水和注入水沟通，频谱曲线反映生产井下方高角度裂缝和中低角度裂缝的总体情况。将注水之后的频谱曲线减去注水之前的频谱曲线，即可得到反映生产井下方中低角度裂缝的频谱曲线。
[0077]
(6)断溶体油藏缝洞发育程度分类评价，按照变换后的频率分布进行细分，频率为0-0.02的部分对应的振幅值高，说明含水变化很小，这部分对应的是大型溶洞体，频率增大到0.1后振幅值整体降低，但变化加快，对应频率为0.02-0.1的是大型裂缝，频率大于0.1的则是中小尺度裂缝。
[0078]
根据各个缝洞类型及其相互之间的比例关系，结合生产井的总产油量或总产气量等总产量信息，可以确定生产井下方各个缝洞类型的体积。
[0079]
傅里叶变换作为一种重要手段，为油藏工程师识别裂缝发育情况带来新思路和手段。通过傅里叶变换，可以将复杂无序的生产数据进行逻辑整合，得到频谱、振幅、功率、密度等频谱分析指标。再将通过本技术的裂缝动态识别结果与地震反射信号的分频结果进行对比，更进一步验证了本技术技术方案的有效性和可靠性。
[0080]
实施例四
[0081]
对于单井缝洞单元，以断溶体油藏th10433h单井缝洞单元为例对本技术的裂缝动态识别方法进行说明。
[0082]
th10433h井处于塔河油田十区，属于断溶体地质背景的一口水平井，th10433h井2014年1月3日完钻，完钻井深6375.00(斜)/6083.00(垂)，层位o2yj。7-5/8
″
套管井口直下，149.2mm裸眼油管测试完井，该井先后经过自喷、转抽、酸化、换井口、定向喷砂射孔完井等措施。
[0083]
该井连通性研究表明该井与周围井不连通，属于孤立的缝洞单井单元，反映的是该区域高角度裂缝的发育情况。考虑到该井2015年到2017年工作制度和措施变化很多，因此需要舍弃2015年到2017年间的数据，专门研究2015年10月之前的生产数据变化，再对生产数据变化进行频谱处理和分析，结果见图3-图6。
[0084]
从该井的含水频谱图分布结果看，频率为0的部分对应的是溶洞信号，根据分析结果，频率为0-0.02的频段对应的储集体空间类型为大型溶洞，其振幅平均值为26.28，频率0.02-0.1之间的频段对应的储集体空间类型为断裂带即大裂缝，其振幅平均值为1.64，剩余的裂缝为中小尺度裂缝，其振幅值平均为0.47。可以看出溶洞的体积是裂缝体积的55.9倍左右，是断裂带部分的16倍左右。再结合该井的总产量信息，能够确定该生产井下方大型溶洞、大裂缝和中小尺度裂缝的体积分别是多少。
[0085]
与如图7所示的地震资料进行对比，可以看出，地震反射图像显示，在th10433h井底附近有串珠状发射特征，而且，钻井过程中有1355方的泥浆漏失，这些信息均说明该井附
近发育一定数量的溶洞和裂缝系统。
[0086]
对于多井缝洞单元，以断溶体油藏s86多井缝洞单元为例对本技术的裂缝动态识别方法进行说明。
[0087]
基本情况：s86井于2001年8月21日完钻，完钻井深5881.50m，层位o1-2y。之后进行裸眼酸压完井，该井先后经过自喷、转抽、酸化、换管柱等措施。
[0088]
地质背景分析：该井处于塔河油田八区，地质背景为断溶体，是多井缝洞单元s86单元的主力产油井。
[0089]
产水分析：该井2010年10月份开始受效，因此之前的含水来源于底水，之后含水变化主要来源于注水井。
[0090]
连通性分析：利用注采关系曲线方法分析连通性，结果表明该井与tk836ch井连通，受到底水和tk836ch注水井影响，属于多井连通的缝洞单元，反映的该区域井间中低角度裂缝发育和分布的情况。
[0091]
注水受效分析：s86井在注水井tk836ch作用下受效，选取2010年10月2日至2011年8月9日含水率变化数据进行处理，结果见图8-图11。
[0092]
从该井的含水频谱图分布结果看，频率为0附近的部分对应的是溶洞信号，根据分析结果，频率为0-0.02的频率段对应的储集体空间类型为大型溶洞，其振幅平均值为31.88，频率为0.02-0.1的频率段对应的储集体空间类型为断裂带即大裂缝，其振幅均值为2.93，剩余的频率段对应的裂缝为中小尺度裂缝，其振幅值平均为0.64，可以看出溶洞的体积是裂缝体积的49.8倍左右，是断裂带的10.88倍左右。结合该井的总产量数据，可以得出，该井下方大型溶洞、大裂缝和中小尺度裂缝各自的体积。
[0093]
与如图12所示的地震资料进行对比，可以看出，地震反射图像在s86井底附近有串珠状发射特征，而且，钻井过程中有3024方的泥浆漏失，这些信息均说明该井附近发育一定数量的溶洞和裂缝系统。
[0094]
本实施例将基于生产动态得到的缝洞发育情况与地震解释中的分频解释结果进行相互对比，实现两种方法的相互验证。在大量的研究基础上，最终形成本技术的定量化的裂缝和溶洞动态识别和描述方法。
[0095]
本技术的裂缝动态识别方法在塔河油田断溶体油藏39个缝洞单元的208口油井展开了应用，以油井含水率分布为应用对象，得到含水率频谱和缝洞动态识别结果，与地震资料处理结果进行了对比，符合率高达92％以上，具有较高的一致性，该方法较好地反映了裂缝的动态变化，为断溶体油藏裂缝动态识别及油井水淹治理奠定了坚实的基础。
[0096]
实施例五
[0097]
本实施例提供一种裂缝动态识别装置，包括：数据变换模块，用于将目标井在不同生产阶段的生产动态数据从时域变换到频域，以得到目标井在不同生产阶段的生产谱图，其中，所述生产动态数据为表征目标井生产能力的一种参数的动态数据；数据分析模块，用于对于每个生产阶段的生产谱图，对生产谱图进行量化分析，根据分析结果确定该生产阶段目标井下方的缝洞类型和各个缝洞类型之间的体积比；体积计算模块，用于对于每个生产阶段，根据目标井在该生产阶段的总产量和目标井下方的缝洞类型及各个缝洞类型之间的体积比，确定在该生产阶段目标井下方各个缝洞类型的体积，从而得到目标井在各个生产阶段的各个缝洞类型的体积。
[0098]
在本实施例的另一个示例中，裂缝动态识别装置还包括存储器和处理器，其中，处理器执行存储在存储器中的以下程序模块：数据变换模块、数据分析模块和体积计算模块，以实现对生产井下方的缝洞储集体类型和体积的动态识别。
[0099]
实施例六
[0100]
本实施例提供一种存储介质，存储有计算机程序，其特征在于，所述计算机程序被处理器执行时，实现如上文所述的裂缝动态识别方法的步骤。
[0101]
存储介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。存储介质的例子包括但不限于相变内存(pram)、静态随机存取存储器(sram)、动态随机存取存储器(dram)、其他类型的随机存取存储器(ram)、只读存储器(rom)、电可擦除可编程只读存储器(eeprom)、快闪记忆体或其他内存技术、只读光盘只读存储器(cd-rom)、数字多功能光盘(dvd)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。
[0102]
实施例七
[0103]
本实施例提供一种计算机设备，包括处理器和存储有程序代码的存储介质，所述程序代码被所述处理器执行时，实现如上文所述的裂缝动态识别方法的步骤。
[0104]
在一个示例中，计算机设备包括一个或多个处理器(cpu)、输入/输出接口、网络接口和内存。
[0105]
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(ram)和/或非易失性内存等形式，如只读存储器(rom)或闪存(flash flash ram)。内存是计算机可读介质的示例。
[0106]
需要注意的是，这里所使用的的术语仅是为了描述具体实施方式，而非意图限制根据本技术的示例性实施方式，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
[0107]
应当理解的是，本说明书中的示例性实施方式可以由多种不同的形式来实施，并且不应当被解释为只限于这里所阐述的实施方式。提供这些实施方式是为了使得本技术的公开彻底且完整，并且将这些示例性实施方式的构思充分传达给本领域普通技术人员，而不应当理解为对本发明的限制。
[0108]
凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。
[0109]
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在如上所述的一个可读存储介质(如rom/ram、磁碟、光盘)中，包括若干指令用以使得一台系统设备(可以是手机，计算机，服务器，空调器，或者网络设备等)执行本发明各个实施例所述的方法。