裂缝储层确定方法、装置、存储介质及电子设备与流程

1.本技术涉及油气勘探技术领域，特别地涉及一种裂缝储层确定方法、装置、存储介质及电子设备。


背景技术：

2.裂缝作为地球表面最小，分布最广，同时也是最复杂的构造。裂缝在油气勘探中，可以作为良好的储层空间和流体运移通道，对储层的连通性、产能和储量有着重要的影响。近年来，针对裂缝型油气藏的识别和描述形成了多种方法技术，包括利用边缘检测、本征值、曲率、蚂蚁体等叠后属性技术以及方位各向异性反演等叠前技术等，但是，由于裂缝成因类型复杂，针对裂缝的各种预测方法和描述参数也有侧重，因此无论是叠后的属性计算还是叠前的各向异性反演，单一的预测方法常常只能对某一类断裂-裂缝或者裂缝发育带进行预测，并且预测结果常常是不同级别、不同尺度的裂缝型储层的综合响应。利用叠后地震数据形成描述裂缝的技术，是定性刻画裂缝的几何学特征。利用叠前方位各向异性反演方法最终可以刻画裂缝的方位和密度，但利用叠后地震资料采用常规的反演技术只能刻画储层的横向特征，不能刻画类似裂缝型储层的纵向特征。


技术实现要素：

3.针对上述问题，本技术提供一种裂缝储层确定方法、装置、存储介质及电子设备。
4.第一方面，本技术提供了一种裂缝储层确定方法，所述方法包括：
5.根据获取的目标区域的地震资料及测井资料建立初始低频模型；
6.根据所述地震资料建立裂缝阻抗模型；
7.根据所述初始低频模型和所述裂缝阻抗模型，得到混合阻抗低频模型；
8.根据所述混合阻抗低频模型进行反演，确定所述目标区域的纵波阻抗数据体，以基于所述纵波阻抗数据体确定所述目标区域的空间展布；
9.获取表示纵波阻抗数据体与孔隙度之间关系的关系函数，以根据所述关系函数基于所述空间展布确定所述裂缝储层。
10.上述实施方式中，根据目标区域的地震资料及测井资料建立初始低频模型，同时根据地震资料建立裂缝阻抗模型，然后根据初始低频模型和裂缝阻抗模型得到的混合阻抗低频模型，从而能够根据混合阻抗低频模型在纵向和横向的立体空间上刻画目标区域的空间展布。对该混合阻抗低频模型进行反演得到目标区域的纵波阻抗数据体，纵波阻抗数据体能够反映目标区域的岩性、物性等信息。因此，根据纵波阻抗数据体及纵波阻抗数据体与孔隙度之间关系的关系函数则可以准确的对目标区域内部缝洞体的空间展布进行刻画，从而能够更直观的刻画储层的缝洞连通关系。
11.根据本技术的实施例，可选的，上述裂缝储层确定方法中，所述获取的目标区域中地震资料及测井资料建立初始低频模型，包括：
12.基于获取的所述测井资料和所述地震资料进行井震标定建立时-深关系；
13.根据所述时-深关系建立所述目标区域中碎屑岩地层的框架模型；
14.根据所述测井资料确定所述目标区域的深度方向的纵波阻抗数据和背景纵波阻抗数据；
15.利用所述纵波阻抗数据对所述框架模型进行横向内插，得到所述碎屑岩地层的低频阻抗模型；
16.基于所述背景纵波阻抗值，确定所述目标区域中碳酸岩地层的背景低频纵波阻抗模型；
17.基于所述低频阻抗模型与所述背景低频纵波阻抗模型确定初始低频模型。
18.根据本技术的实施例，可选的，上述裂缝储层确定方法中，所述测井资料包括：单井的声波数据和密度数据，所述地震资料包括：地震子波数据，所述基于获取的所述测井资料和所述地震资料进行井震标定建立时-深关系，包括：
19.基于单井的声波数据和密度数据确定反射系数；
20.对反射系数和所述地震子波进行褶积，得到合成地震记录；
21.基于所述合成地震记录建立所述时-深关系。
22.根据本技术的实施例，可选的，上述裂缝储层确定方法中，根据所述地震资料建立裂缝阻抗模型，包括：
23.对所述地震资料进行构造导向滤波处理，得到滤波后的地震数据；
24.基于所述滤波后的地震数据进行断层自动检测afe计算，得到所述目标区域的中间裂缝阻抗模型；所述中间裂缝阻抗模型包括afe值；
25.获取所述目标区域对应的afe阈值，并基于所述afe阈值以及所述中间裂缝阻抗模型中的afe值确定所述裂缝阻抗模型。
26.根据本技术的实施例，可选的，上述裂缝储层确定方法中，基于所述afe阈值以及所述中间裂缝阻抗模型中的afe值确定所述裂缝阻抗模型，包括：
27.获取所述目标区域对应的常数阻抗值；
28.将所述中间裂缝阻抗模型中大于所述afe阈值的afe值确定为所述常数阻抗值，并将所述中间裂缝阻抗模型中小于所述afe阈值的afe值确定为空值，以得到所述裂缝阻抗模型。
29.根据本技术的实施例，可选的，上述裂缝储层确定方法中，根据所述混合阻抗低频模型进行反演，确定所述目标区域的纵波阻抗数据体，以基于所述纵波阻抗数据体确定所述目标区域的空间展布包括：
30.基于所述地震资料采用所述混合阻抗低频模型进行约束稀疏脉冲反演，得到所述目标区域的纵波阻抗数据；
31.基于所述纵波阻抗数据和纵波阻抗阈值，确定目标纵波阻抗数据；
32.将所述目标纵波阻抗数据对应的空间位置确定为所述目标区域的空间展布。
33.根据本技术的实施例，可选的，上述裂缝储层确定方法中，获取表示纵波阻抗数据体与孔隙度之间关系的关系函数，以根据所述关系函数基于所述空间展布确定所述裂缝储层，包括：
34.根据所述关系函数以及所述目标区域的纵波阻抗数据体计算所述目标区域的孔隙度；
35.根据所述孔隙度确定所述裂缝储层。
36.第二方面，本技术提供了一种裂缝储层确定装置，所述装置包括：
37.初始低频模型建立模块，用于根据获取的目标区域的地震资料及测井资料建立初始低频模型；
38.裂缝阻抗模型建立模块，用于根据所述地震资料建立裂缝阻抗模型；
39.混合阻抗低频模型确定模块，用于根据所述初始低频模型和所述裂缝阻抗模型，得到混合阻抗低频模型；
40.纵波阻抗数据体确定模块，用于根据所述混合阻抗低频模型进行反演，确定所述目标区域的纵波阻抗数据体，以基于所述纵波阻抗数据体确定所述目标区域的空间展布；
41.裂缝确定模块，用于获取表示纵波阻抗数据体与孔隙度之间关系的关系函数，以根据所述关系函数基于所述空间展布确定所述裂缝储层。
42.根据本技术的实施例，可选的，上述裂缝储层确定装置中，所述初始低频模型建立模块包括：
43.井震标定单元，用于基于获取的所述测井资料和所述地震资料进行井震标定建立时-深关系；
44.框架模型建立单元，用于根据所述时-深关系建立所述目标区域中碎屑岩地层的框架模型；
45.阻抗数据确定单元，用于根据所述测井资料确定所述目标区域的深度方向的纵波阻抗数据和背景纵波阻抗数据；
46.低频阻抗模型获取单元，用于利用所述纵波阻抗数据对所述框架模型进行横向内插，得到所述碎屑岩地层的低频阻抗模型；
47.背景低频纵波阻抗模型确定单元，用于基于所述背景纵波阻抗值，确定所述目标区域中碳酸岩地层的背景低频纵波阻抗模型；
48.初始低频模型确定单元，用于基于所述低频阻抗模型与所述背景低频纵波阻抗模型确定初始低频模型。
49.根据本技术的实施例，可选的，上述裂缝储层确定装置中，所述测井资料包括：单井的声波数据和密度数据，所述地震资料包括：地震子波数据，所述井震标定单元包括：
50.反射系数确定子单元，用于基于单井的声波数据和密度数据确定反射系数；
51.合成地震记录获取子单元，用于对反射系数和所述地震子波进行褶积，得到合成地震记录；
52.时-深关系建立子单元，用于基于所述合成地震记录建立所述时-深关系。
53.根据本技术的实施例，可选的，上述裂缝储层确定装置中，裂缝阻抗模型建立模块包括：
54.构造导向滤波单元，用于对所述地震资料进行构造导向滤波处理，得到滤波后的地震数据；
55.断层自动检测单元，用于基于所述滤波后的地震数据进行断层自动检测afe计算，得到所述目标区域的中间裂缝阻抗模型；所述中间裂缝阻抗模型包括afe值；
56.裂缝阻抗模型确定单元，用于获取所述目标区域对应的afe阈值，并基于所述afe阈值以及所述中间裂缝阻抗模型中的afe值确定所述裂缝阻抗模型。
57.根据本技术的实施例，可选的，上述裂缝储层确定装置中，裂缝阻抗模型确定单元包括：
58.常数阻抗值获取子单元，用于获取所述目标区域对应的常数阻抗值；
59.裂缝阻抗模型子单元，用于将所述中间裂缝阻抗模型中大于所述afe阈值的afe值确定为所述常数阻抗值，并将所述中间裂缝阻抗模型中小于所述afe阈值的afe值确定为空值，以得到所述裂缝阻抗模型。
60.根据本技术的实施例，可选的，上述裂缝储层确定装置中，纵波阻抗数据体确定模块包括：
61.纵波阻抗数据获取单元，用于基于所述地震资料采用所述混合阻抗低频模型进行约束稀疏脉冲反演，得到所述目标区域的纵波阻抗数据；
62.目标纵波阻抗数据确定单元，用于基于所述纵波阻抗数据和纵波阻抗阈值，确定目标纵波阻抗数据；
63.空间确定单元，用于将所述目标纵波阻抗数据对应的空间位置确定为所述目标区域的空间展布。
64.根据本技术的实施例，可选的，上述裂缝储层确定装置中，裂缝确定模块包括：
65.孔隙度计算单元，用于根据所述关系函数以及所述目标区域的纵波阻抗数据体计算所述目标区域的孔隙度；
66.裂缝确定单元，用于根据所述孔隙度确定所述裂缝储层。
67.第三方面，本技术提供了一种存储介质，该存储介质存储的计算机程序，可被一个或多个处理器执行，可用来实现如上述的裂缝储层确定方法。
68.第四方面，本技术提供了一种电子设备，包括存储器和处理器，所述存储器上存储有计算机程序，该计算机程序被所述处理器执行时，执行上述的裂缝储层确定方法。
69.综上所述，本技术提供的一种裂缝储层确定方法、装置、存储介质及电子设备，所述方法包括：根据获取的目标区域的地震资料及测井资料建立初始低频模型；根据所述地震资料建立裂缝阻抗模型；根据所述初始低频模型和所述裂缝阻抗模型，得到混合阻抗低频模型；根据所述混合阻抗低频模型进行反演，确定所述目标区域的纵波阻抗数据体，以基于所述纵波阻抗数据体确定所述目标区域的空间展布；获取表示纵波阻抗数据体与孔隙度之间关系的关系函数，以根据所述关系函数基于所述空间展布确定所述裂缝储层。根据目标区域的地震资料及测井资料建立初始低频模型，同时根据地震资料建立裂缝阻抗模型，然后根据初始低频模型和裂缝阻抗模型得到的混合阻抗低频模型，从而能够根据混合阻抗低频模型在纵向和横向的立体空间上刻画储层展布。对该混合阻抗低频模型进行反演得到目标区域的纵波阻抗数据体，纵波阻抗数据体能够反映目标区域的岩性、物性等信息。因此，根据纵波阻抗数据体及纵波阻抗数据体与孔隙度之间关系的关系函数则可以准确的对目标区域内部缝洞体的空间展布进行刻画，从而能够更直观的刻画储层的缝洞连通关系。进一步则可以准确的计算油气的资源规模，为井位部署提供可靠依据，从而能够提高勘探开发钻井的成功率，继而提升勘探开发效益，且该方法不受地域限制适用范围广。
附图说明
70.在下文中将基于实施例并参考附图来对本技术进行更详细的描述。
71.图1为本技术实施例一提供的一种裂缝储层确定方法的流程示意图。
72.图2为本技术实施例四提供的一种裂缝储层确定装置的结构框图。
73.图3为本技术实施例六提供的一种电子设备的连接框图。
74.在附图中，相同的部件使用相同的附图标记，附图并未按照实际的比例绘制。
具体实施方式
75.以下将结合附图及实施例来详细说明本技术的实施方式，借此对本技术如何应用技术手段来解决技术问题，并达到相应技术效果的实现过程能充分理解并据以实施。本技术实施例以及实施例中的各个特征，在不相冲突前提下可以相互结合，所形成的技术方案均在本技术的保护范围之内。
76.实施例一
77.本发明提供一种裂缝储层确定方法，请参阅图1，该方法包括如下步骤：
78.步骤s110：根据获取的目标区域的地震资料及测井资料建立初始低频模型。
79.根据地震资料和测井资料可以建立基本反映沉积体地质特征的初始低频模型。低频模型可以是以地震资料解释层位和沉积规律为约束，将测井资料在整个数据体范围内进行内插和外推得到的，此外，还可以将地震资料中的地震速度谱信息与测井资料相结合来建立初始低频模型，可以在一定程度上补偿了地震资料中缺失的低频信息。在根据地震资料以及测井资料建立初始低频模型时，可以选取测井资料中质量较好的井曲线，在层位、断层等约束下按照一定算法进行插值，如权重法、克里金、反距离加权等。除了上述提到的借助地震速度转换建立低频模型，还可以使用单一井插值建立初始低频模型，通过正演对比合成记录与原始地震之间的差异，在差异较大的地方加入伪井，改变低频模型信息，通过不断更新得到相对真实的低频模型。为了考虑到地震反射特征，还可以通过地震振幅、频率等属性划定不同岩性、流体的边界，对各个相带定义不同的弹性性质，以得到初始低频模型。
80.步骤s120：根据所述地震资料建立裂缝阻抗模型。
81.构造形变是裂缝性储层形成的主要原因之一，构造形变和裂缝具有同源性，构造性变对裂缝空间展布特征影响较大，断层是构造形变的一种重要表现形式。为断层在勘探阶段受重视成都较小，但其对油田的开发影响较大，裂缝既能改善储层的渗透性，又能作为烃类储层。地震资料的几何学特征是有效描述裂缝储层的重要属性，常见的地震资料的几何学地震属性有地震资料的相干属性、曲率、倾角和方位角等，因此，在根据地震资料建立裂缝阻抗模型时，可以主要基于上述的几何学地震属性。具体地，确定地震资料中所需要的用到的几何学地震属性数据，然后进行本征值相干体制作，在相干体基础之上，进行断层自动检测计算(automatic fault extraction，afe)，继而根据断层自动检测计算结果建立缝隙阻抗模型。
82.步骤s130：根据所述初始低频模型和所述裂缝阻抗模型，得到混合阻抗低频模型。
83.具体地，可以利用所述低频模型替代所述裂缝阻抗模型中对应的空值，以得到混合阻抗低频模型。裂缝阻抗模型中有空值部分，把裂缝阻抗模型中的空值部分用步骤s110中的初始低频模型对应的值来代替，就可以得到混合阻抗低频模型。
84.根据初始低频模型可以刻画目标区域横向特征，而利用裂缝阻抗模型则可以定性刻画裂缝的纵向几何学特征，根据初始低频模型和裂缝阻抗模型得到的混合阻抗低频模型
则可以在纵向和横向的立体空间上刻画目标区域的空间展布。
85.步骤s140：根据所述混合阻抗低频模型进行反演，确定所述目标区域的纵波阻抗数据体，以基于所述纵波阻抗数据体确定所述目标区域的空间展布。
86.根据本技术的实施例，步骤s140包括以下步骤：
87.步骤s141：基于所述地震资料采用所述混合阻抗低频模型进行约束稀疏脉冲反演，得到所述目标区域的纵波阻抗数据。
88.根据地震资料对混合阻抗低频模型进行反演可以充分利用地震资料和测井资料提供的构造、层位、岩性等信息将常规的地震反射振幅的变化，转换成目标区域的纵波阻抗数据，以此来反映目标区域的岩性、物性等信息。目标区域的纵波阻抗数据可以准确的分析基岩矿物成分及孔隙度。在暴露期，基岩由浅至深受风化淋滤作用的强度越来越弱，孔隙度越来越小，纵波阻抗越来越大，当基岩完全不受地表水风化淋滤作用时，孔隙度接近于零，此时纵波阻抗值最大且只与组成岩石的矿物成分相关。
89.具体地，可以先对地震资料进行细致处理，以提取地震体、地震子波等参数，然后再对初始低频模型进行约束稀疏脉冲反演。约束稀疏脉冲反演是一种递推地震波阻抗反演方法，约束稀疏脉冲反演通过调整反射系数序列的稀疏性来拓宽输入地震数据的有效频带宽度，并得到一个弹性参数模型、稀疏性约束因子，地震信噪比、合并频率、子波刻度因子为该算法中的关键参数，地震信噪比用于约束反演结果与地震数据的相似性，信噪比设置的越高，就表示从反演结果中转换的合成记录与地震越相关，反之亦然、稀疏性约束因子是反射系数序列的稀疏性，稀疏性约束因子的值越小，反射系数序列越稀疏。
90.步骤s142：基于所述纵波阻抗数据和纵波阻抗阈值，确定目标纵波阻抗数据。
91.在通过约束稀疏脉冲反演得到反演结果后，可以获取纵波阻抗阈值，在获取纵波阻抗阈值时，可以根据实际钻探钻到的洞穴层储层，通过井震标定，确定“串珠状”反射对应洞穴层储层在反演得到的纵波阻抗剖面上的位置所对应的纵波阻抗值，结合实钻井得到洞穴层储层对应的纵波阻抗值，就确定洞穴层储层的纵波阻抗阈值。然后根据纵波阻抗数据和纵波阻抗阈值，从所述纵波阻抗数据中确定出目标纵波阻抗数据。
92.具体地，在确定目标纵波阻抗数据时，可以基于以下过程进行确定。首先比较所述纵波阻抗数据和所述纵波阻抗阈值的大小；然后将纵波阻抗数据中小于所述纵波阻抗阈值的部分确定为目标纵波阻抗数据。例如，在碳酸盐岩地层中，将纵波阻抗数据小于这个纵波阻抗阈值确定为洞穴层储层的目标纵波阻抗数据。根据纵波阻抗阈值对目标低频模型进行约束稀疏脉冲反演，得到所述目标区域的纵波阻抗数据进行刻画。纵波阻抗数据中小于该纵波阻抗阈值的保留下来并确定为目标纵波阻抗数据，纵波阻抗数据中大于纵波阻抗阈值的则舍弃掉，从而根据保留下的目标纵波阻抗数据部分就可以对目标区域的空间进行刻画。
93.步骤s143：将所述目标纵波阻抗数据对应的空间位置确定为所述目标区域的空间展布。
94.得到目标纵波阻抗数据之后，利用得到的全频带的目标纵波阻抗体能够更准确的对洞穴储层进行刻画。
95.步骤s150：获取表示纵波阻抗数据体与孔隙度之间关系的关系函数，以根据所述关系函数基于所述空间展布确定所述裂缝储层。
96.纵波阻抗数据体是与孔隙度之间存在关系的，因此，可以预先根据实际钻探获取资料中的纵波阻抗数据体及孔隙度，建立两者之间相对应的函数关系，由此可以根据目标纵波阻抗数据以及函数关系计算目标区域在空间上的孔隙分布，从而实现对目标区域中裂缝储层的确定。
97.具体的，在获取表示纵波阻抗数据体与孔隙度之间关系的关系函数，以根据所述关系函数基于所述空间展布确定所述裂缝储层时，可以先根据所述关系函数以及所述目标区域的纵波阻抗数据体计算所述目标区域的孔隙度；然后根据所述孔隙度确定所述裂缝储层。
98.综上所述，本技术提供一种裂缝储层确定方法，包括：根据获取的目标区域的地震资料及测井资料建立初始低频模型；根据所述地震资料建立裂缝阻抗模型；根据所述初始低频模型和所述裂缝阻抗模型，得到混合阻抗低频模型；根据所述混合阻抗低频模型进行反演，确定所述目标区域的纵波阻抗数据体，以基于所述纵波阻抗数据体确定所述目标区域的空间展布；获取表示纵波阻抗数据体与孔隙度之间关系的关系函数，以根据所述关系函数基于所述空间展布确定所述裂缝储层。根据目标区域的地震资料及测井资料建立初始低频模型，同时根据地震资料建立裂缝阻抗模型，然后根据初始低频模型和裂缝阻抗模型得到的混合阻抗低频模型，从而能够根据混合阻抗低频模型在纵向和横向的立体空间上刻画储层展布。对该混合阻抗低频模型进行反演得到目标区域的纵波阻抗数据体，纵波阻抗数据体能够反映目标区域的岩性、物性等信息。因此，根据纵波阻抗数据体及纵波阻抗数据体与孔隙度之间关系的关系函数则可以准确的对目标区域内部缝洞体的空间展布进行刻画，从而能直观刻画储层的缝洞连通关系。
99.实施例二
100.在实施例一的基础上，本实施例通过具体实施案例对实施例一中的方法进行说明。
101.本技术提供的裂缝储层确定方法包括：
102.步骤s110：根据获取的目标区域的地震资料及测井资料建立初始低频模型。
103.步骤s120：根据所述地震资料建立裂缝阻抗模型。
104.步骤s130：根据所述初始低频模型和所述裂缝阻抗模型，得到混合阻抗低频模型。
105.步骤s140：根据所述混合阻抗低频模型进行反演，确定所述目标区域的纵波阻抗数据体，以基于所述纵波阻抗数据体确定所述目标区域的空间展布。
106.步骤s150：获取表示纵波阻抗数据体与孔隙度之间关系的关系函数，以根据所述关系函数基于所述空间展布确定所述裂缝储层。
107.其中，上述裂缝储层确定方法中，所述步骤s110包括以下步骤：
108.步骤s1110：基于获取的所述测井资料和所述地震资料进行井震标定建立时-深关系。
109.利用地震资料和测井资料可以对目标区域的构造进行解释，并进行储层的预测，以对油藏进行精细的描述。但是测井资料描述的纵向刻度为深度，地震资料所描述的剖面为时间刻度，两种资料描述方式不同，因此无法直接联合应用，需要先根据测井资料与地震资料间的井震标定。在进行井震标定时，可以基于地震资料中的合成地震记录与井旁地震道波形相关性来标定。此外，还可以通过声波、密度测井资料计算反射系数，利用发射系数
与子波的褶积构建与地震道相似的合成地震记录，再通过将合成地震记录与井旁地震道进行对比的方式，对标定结果进行调整来实现井震标定。例如，可以先确定目标区域样本钻井，并依据该样本钻井获取单晶深度域合成地震记录，然后根据单井深度域合成地震记录，建立单井处时间和深度的对应关系。接着将单井的地质层位由深度域转为时间域，建立时间域地质层位模型，最后将研究区剩余的钻井层位自动匹配到时间域地质层位模型上，从而实现全部钻井的快速井震标定工作。
110.步骤s1120：根据所述时-深关系建立所述目标区域中碎屑岩地层的框架模型。
111.井震标定得到的时-深关系，可以将单井分析的不同地质界面对应到地震剖面上，也就是从深度域对应到时间域，因此，可以利用地震资料对不同地质界面进行空间追踪，得到地质界面在空间上的分布，利用这些不同的地质界面建立碎屑岩地层的框架模型。
112.步骤s1130：根据所述测井资料确定所述目标区域的深度方向的纵波阻抗数据和背景纵波阻抗数据；
113.步骤s1140：利用所述纵波阻抗数据对所述框架模型进行横向内插，得到所述碎屑岩地层的低频阻抗模型；
114.步骤s1150：基于所述背景纵波阻抗值，确定所述目标区域中碳酸岩地层的背景低频纵波阻抗模型；
115.步骤s1160：基于所述低频阻抗模型与所述背景低频纵波阻抗模型确定初始低频模型。
116.在框架模型的约束下，把不同单井纵向上的纵波阻抗值进行横向内插，得到井内插的低频阻抗模型。由于碳酸盐岩在横向上非均质性强，单井横向内插的模型不能反映碳酸盐岩地层横向分布的特征。因此，可以对碳酸盐岩顶面(t74地震反射界面)以上的碎屑岩地层用井内插的低频模型；t74以下是碳酸盐岩地层，用碳酸盐岩地层的背景纵波阻抗值，来替代低频模型，就是碳酸盐岩地层部分的低频模型是一个常数的直板低频模型。其中，背景纵波阻抗值在这里可以用一个常数代表碳酸盐岩背景值。通过以上方式，可以在碳酸盐岩地层中的地震资料中的异常值能反映储层的变化，消除了由于井内插的低频模型异常带来的假象，并真正反映地下的实际情况。
117.其中，上述裂缝储层确定方法中，所述测井资料包括：单井的声波数据和密度数据，所述地震资料包括：地震子波数据，所述步骤s1110包括以下步骤：
118.步骤s1111：基于单井的声波数据和密度数据确定反射系数。
119.将时间域子波转换到深度域后，由于速度的影响，波形随深度发生压缩或拉伸两种变化，因此可以基于声波数据和密度数据计算反射系数。在计算反射系数时，要在深度偏移剖面上选取品质较好的井旁地震道计算子波的主频，要选用零相位子波，在计算反射系数之前，要对声波时差曲线和密度曲线进行校正、去野值等处理。
120.步骤s1112：对反射系数和所述地震子波进行褶积，得到合成地震记录。
121.褶积，又名卷积，是一种积分变换的数学方法，是通过两个函数生成第三个函数的一种数学算子，表征两个函数之间经过翻转和平移的重叠部分函数值乘积对重叠长度的积分。通过褶积获取的合成地震记录是用声波测井或垂直地震剖面资料经过人工合成转换成的地震记录，即地震道。是地震模型技术中应用非常广泛的一种，也是层位标定、油藏描述等工作的基础，是把地质模型转化为地震信息的中间媒介。合成地震记录是联合高分辨率
的测井信息与区域性的地震信息的桥梁，其精度直接影响到地质层位的准确标定。
122.步骤s1113：基于所述合成地震记录建立所述时-深关系。
123.井资料是深度域，地震资料是时间域的。通过单井的声波和密度资料得到反射系数，地震子波与反射系数褶积得到合成地震记录，合成记录剖面与过井的地震剖面进行标定，把深度域的资料转换到时间域。通过对反射系数和地震子波进行褶积得到的合成地震记录，能够准确的建立出时-深关系，从而保证准确的对裂缝储层进行确定。
124.实施例三
125.在实施例一的基础上，本实施例通过具体实施案例对实施例一中的方法进行说明。
126.在根据所述地震资料建立裂缝阻抗模型时，可以采用以下过程。首先对所述地震资料进行构造导向滤波处理，得到滤波后的地震数据；然后基于所述滤波后的地震数据进行断层自动检测afe计算，得到所述目标区域的中间裂缝阻抗模型；所述中间裂缝阻抗模型包括afe值；最后获取所述目标区域对应的afe阈值，并基于所述afe阈值以及所述中间裂缝阻抗模型中的afe值确定所述裂缝阻抗模型。
127.其中，基于所述afe阈值以及所述中间裂缝阻抗模型中的afe值确定所述裂缝阻抗模型时，可以先获取所述目标区域对应的常数阻抗值；然后将所述中间裂缝阻抗模型中大于所述afe阈值的afe值确定为所述常数阻抗值，并将所述中间裂缝阻抗模型中小于所述afe阈值的afe值确定为空值，以得到所述裂缝阻抗模型。
128.具体地，对原始的地震资料进行构造导向滤波处理，通过对地震数据的构造导向滤波，可以提高地震数据的信噪比，使地震同相轴的连续特征或间断特征更明显，方便地震数据的层位和断层解释，同时也为后续相关性计算与断层识别提供资料基础。在进行滤波后的地震数据的基础上进行afe计算，以得到afe值。此外，在获取目标区域对应的afe阈值时，可以先获取已钻井钻到的裂缝相关资料，然后根据获取的相关资料进行井震标定，确定裂缝所在位置对应的afe的值，该值可以作为裂缝的afe阈值。基于所述afe阈值以及所述中间裂缝阻抗模型中的afe值确定所述裂缝阻抗模型时，将afe值中大于afe阈值的afe值用一个常数碳酸盐岩背景的阻抗值代替，即目标区域中的裂缝；afe值中小于等于afe阈值的afe值设置为空值，也就是目标区域中的非裂缝，从而建立裂缝的阻抗模型。
129.裂缝阻抗模型中有空值部分，把裂缝阻抗模型中的空值部分用步骤s110中的初始低频模型对应的值来代替，就可以得到混合阻抗低频模型。
130.实施例四
131.请参看图2，本技术提供了一种裂缝储层确定装置200，该装置包括：
132.初始低频模型建立模块210，用于根据获取的目标区域的地震资料及测井资料建立初始低频模型；
133.裂缝阻抗模型建立模块220，用于根据所述地震资料建立裂缝阻抗模型；
134.混合阻抗低频模型确定模块230，用于根据所述初始低频模型和所述裂缝阻抗模型，得到混合阻抗低频模型；
135.纵波阻抗数据体确定模块240，用于根据所述混合阻抗低频模型进行反演，确定所述目标区域的纵波阻抗数据体，以基于所述纵波阻抗数据体确定所述目标区域的空间展布；
136.裂缝确定模块250，用于获取表示纵波阻抗数据体与孔隙度之间关系的关系函数，以根据所述关系函数基于所述空间展布确定所述裂缝储层。
137.根据本技术的实施例，可选的，上述裂缝储层确定装置中，所述初始低频模型建立模块包括：
138.井震标定单元，用于基于获取的所述测井资料和所述地震资料进行井震标定建立时-深关系；
139.框架模型建立单元，用于根据所述时-深关系建立所述目标区域中碎屑岩地层的框架模型；
140.阻抗数据确定单元，用于根据所述测井资料确定所述目标区域的深度方向的纵波阻抗数据和背景纵波阻抗数据；
141.低频阻抗模型获取单元，用于利用所述纵波阻抗数据对所述框架模型进行横向内插，得到所述碎屑岩地层的低频阻抗模型；
142.背景低频纵波阻抗模型确定单元，用于基于所述背景纵波阻抗值，确定所述目标区域中碳酸岩地层的背景低频纵波阻抗模型；
143.初始低频模型确定单元，用于基于所述低频阻抗模型与所述背景低频纵波阻抗模型确定初始低频模型。
144.根据本技术的实施例，可选的，上述裂缝储层确定装置中，所述测井资料包括：单井的声波数据和密度数据，所述地震资料包括：地震子波数据，所述井震标定单元包括：
145.反射系数确定子单元，用于基于单井的声波数据和密度数据确定反射系数；
146.合成地震记录获取子单元，用于对反射系数和所述地震子波进行褶积，得到合成地震记录；
147.时-深关系建立子单元，用于基于所述合成地震记录建立所述时-深关系。
148.根据本技术的实施例，可选的，上述裂缝储层确定装置中，裂缝阻抗模型建立模块包括：
149.构造导向滤波单元，用于对所述地震资料进行构造导向滤波处理，得到滤波后的地震数据；
150.断层自动检测单元，用于基于所述滤波后的地震数据进行断层自动检测afe计算，得到所述目标区域的中间裂缝阻抗模型；所述中间裂缝阻抗模型包括afe值；
151.裂缝阻抗模型确定单元，用于获取所述目标区域对应的afe阈值，并基于所述afe阈值以及所述中间裂缝阻抗模型中的afe值确定所述裂缝阻抗模型。
152.根据本技术的实施例，可选的，上述裂缝储层确定装置中，裂缝阻抗模型确定单元包括：
153.常数阻抗值获取子单元，用于获取所述目标区域对应的常数阻抗值；
154.裂缝阻抗模型子单元，用于将所述中间裂缝阻抗模型中大于所述afe阈值的afe值确定为所述常数阻抗值，并将所述中间裂缝阻抗模型中小于所述afe阈值的afe值确定为空值，以得到所述裂缝阻抗模型。
155.根据本技术的实施例，可选的，上述裂缝储层确定装置中，纵波阻抗数据体确定模块包括：
156.纵波阻抗数据获取单元，用于基于所述地震资料采用所述混合阻抗低频模型进行
约束稀疏脉冲反演，得到所述目标区域的纵波阻抗数据；
157.目标纵波阻抗数据确定单元，用于基于所述纵波阻抗数据和纵波阻抗阈值，确定目标纵波阻抗数据；
158.空间确定单元，用于将所述目标纵波阻抗数据对应的空间位置确定为所述目标区域的空间展布。
159.根据本技术的实施例，可选的，上述裂缝储层确定装置中，裂缝确定模块包括：
160.孔隙度计算单元，用于根据所述关系函数以及所述目标区域的纵波阻抗数据体计算所述目标区域的孔隙度；
161.裂缝确定单元，用于根据所述孔隙度确定所述裂缝储层。
162.综上所述，本技术提供一种裂缝储层确定装置，包括：初始低频模型建立模块，用于根据获取的目标区域的地震资料及测井资料建立初始低频模型；裂缝阻抗模型建立模块，用于根据所述地震资料建立裂缝阻抗模型；混合阻抗低频模型确定模块，用于根据所述初始低频模型和所述裂缝阻抗模型，得到混合阻抗低频模型；纵波阻抗数据体确定模块，用于根据所述混合阻抗低频模型进行反演，确定所述目标区域的纵波阻抗数据体，以基于所述纵波阻抗数据体确定所述目标区域的空间展布；裂缝确定模块，用于获取表示纵波阻抗数据体与孔隙度之间关系的关系函数，以根据所述关系函数基于所述空间展布确定所述裂缝储层。根据目标区域的地震资料及测井资料建立初始低频模型，同时根据地震资料建立裂缝阻抗模型，然后根据初始低频模型和裂缝阻抗模型得到的混合阻抗低频模型，从而能够根据混合阻抗低频模型在纵向和横向的立体空间上刻画储层展布。对该混合阻抗低频模型进行反演得到目标区域的纵波阻抗数据体，纵波阻抗数据体能够反映目标区域的岩性、物性等信息。因此，根据纵波阻抗数据体及纵波阻抗数据体与孔隙度之间关系的关系函数则可以准确的对目标区域内部缝洞体的空间展布进行刻画，从而能直观刻画储层的缝洞连通关系。
163.实施例五
164.本实施例还提供一种计算机可读存储介质，如闪存、硬盘、多媒体卡、卡型存储器(例如，sd或dx存储器等)、随机访问存储器(ram)、静态随机访问存储器(sram)、只读存储器(rom)、电可擦除可编程只读存储器(eeprom)、可编程只读存储器(prom)、磁性存储器、磁盘、光盘、服务器、app应用商城等等，其上存储有计算机程序，所述计算机程序被处理器执行时可以实现如下方法步骤：
165.步骤s110：根据获取的目标区域的地震资料及测井资料建立初始低频模型。
166.步骤s120：根据所述地震资料建立裂缝阻抗模型。
167.步骤s130：根据所述初始低频模型和所述裂缝阻抗模型，得到混合阻抗低频模型。
168.步骤s140：根据所述混合阻抗低频模型进行反演，确定所述目标区域的纵波阻抗数据体，以基于所述纵波阻抗数据体确定所述目标区域的空间展布。
169.步骤s150：获取表示纵波阻抗数据体与孔隙度之间关系的关系函数，以根据所述关系函数基于所述空间展布确定所述裂缝储层。
170.可选的，上述裂缝储层确定方法中，步骤s110，包括以下步骤：
171.基于获取的所述测井资料和所述地震资料进行井震标定建立时-深关系；
172.根据所述时-深关系建立所述目标区域中碎屑岩地层的框架模型；
173.根据所述测井资料确定所述目标区域的深度方向的纵波阻抗数据和背景纵波阻抗数据；
174.利用所述纵波阻抗数据对所述框架模型进行横向内插，得到所述碎屑岩地层的低频阻抗模型；
175.基于所述背景纵波阻抗值，确定所述目标区域中碳酸岩地层的背景低频纵波阻抗模型；
176.基于所述低频阻抗模型与所述背景低频纵波阻抗模型确定初始低频模型。
177.可选的，上述裂缝储层确定方法中，所述测井资料包括：单井的声波数据和密度数据，所述地震资料包括：地震子波数据，所述基于获取的所述测井资料和所述地震资料进行井震标定建立时-深关系，包括：
178.基于单井的声波数据和密度数据确定反射系数；
179.对反射系数和所述地震子波进行褶积，得到合成地震记录；
180.基于所述合成地震记录建立所述时-深关系。
181.可选的，上述裂缝储层确定方法中，根据所述地震资料建立裂缝阻抗模型，包括：
182.对所述地震资料进行构造导向滤波处理，得到滤波后的地震数据；
183.基于所述滤波后的地震数据进行断层自动检测afe计算，得到所述目标区域的中间裂缝阻抗模型；所述中间裂缝阻抗模型包括afe值；
184.获取所述目标区域对应的afe阈值，并基于所述afe阈值以及所述中间裂缝阻抗模型中的afe值确定所述裂缝阻抗模型。
185.可选的，上述裂缝储层确定方法中，基于所述afe阈值以及所述中间裂缝阻抗模型中的afe值确定所述裂缝阻抗模型，包括：
186.获取所述目标区域对应的常数阻抗值；
187.将所述中间裂缝阻抗模型中大于所述afe阈值的afe值确定为所述常数阻抗值，并将所述中间裂缝阻抗模型中小于所述afe阈值的afe值确定为空值，以得到所述裂缝阻抗模型。
188.可选的，上述裂缝储层确定方法中，根据所述混合阻抗低频模型进行反演，确定所述目标区域的纵波阻抗数据体，以基于所述纵波阻抗数据体确定所述目标区域的空间展布包括：
189.基于所述地震资料采用所述混合阻抗低频模型进行约束稀疏脉冲反演，得到所述目标区域的纵波阻抗数据；
190.基于所述纵波阻抗数据和纵波阻抗阈值，确定目标纵波阻抗数据；
191.将所述目标纵波阻抗数据对应的空间位置确定为所述目标区域的空间展布。
192.可选的，上述裂缝储层确定方法中，获取表示纵波阻抗数据体与孔隙度之间关系的关系函数，以根据所述关系函数基于所述空间展布确定所述裂缝储层，包括：
193.根据所述关系函数以及所述目标区域的纵波阻抗数据体计算所述目标区域的孔隙度；
194.根据所述孔隙度确定所述裂缝储层。
195.上述方法步骤的具体实施例过程可参见实施例一，本实施例在此不再重复赘述。
196.实施例六
197.本技术实施例提供了一种电子设备，该电子设备可以是手机、电脑或平板电脑等，包括存储器和处理器，所述存储器上存储有计算器程序，该计算机程序被处理器执行时实现如实施例一中所述的裂缝储层确定方法。可以理解，如图3所示，该电子设备300还可以包括：处理器301，存储器302，多媒体组件303，输入/输出(i/o)接口304，以及通信组件305。
198.其中，处理器301用于执行如实施例一中的裂缝储层确定方法中的全部或部分步骤。存储器302用于存储各种类型的数据，这些数据例如可以包括电子设备中的任何应用程序或方法的指令，以及应用程序相关的数据。
199.处理器301可以是专用集成电路(application specific integrated circuit，简称asic)、数字信号处理器(digital signal processor，简称dsp)、数字信号处理设备(digital signal processing device，简称dspd)、可编程逻辑器件(programmable logic device，简称pld)、现场可编程门阵列(field programmable gate array，简称fpga)、控制器、微控制器、微处理器或其他电子元件实现，用于执行上述实施例一中的裂缝储层确定方法。
200.存储器302可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，例如静态随机存取存储器(static random access memory，简称sram)，电可擦除可编程只读存储器(electrically erasable programmable read-only memory，简称eeprom)，可擦除可编程只读存储器(erasable programmable read-only memory，简称eprom)，可编程只读存储器(programmable read-only memory，简称prom)，只读存储器(read-only memory，简称rom)，磁存储器，快闪存储器，磁盘或光盘。
201.多媒体组件303可以包括屏幕和音频组件，该屏幕可以是触摸屏，音频组件用于输出和/或输入音频信号。例如，音频组件可以包括一个麦克风，麦克风用于接收外部音频信号。所接收的音频信号可以被进一步存储在存储器或通过通信组件发送。音频组件还包括至少一个扬声器，用于输出音频信号。
202.i/o接口304为处理器301和其他接口模块之间提供接口，上述其他接口模块可以是键盘，鼠标，按钮等。这些按钮可以是虚拟按钮或者实体按钮。
203.通信组件305用于该电子设备300与其他设备之间进行有线或无线通信。无线通信，例如wi-fi，蓝牙，近场通信(near field communication，简称nfc)，2g、3g或4g，或它们中的一种或几种的组合，因此相应的该通信组件305可以包括：wi-fi模块，蓝牙模块，nfc模块。
204.综上，本技术提供的一种裂缝储层确定方法、装置、存储介质及电子设备，该方法包括：根据获取的目标区域的地震资料及测井资料建立初始低频模型；根据所述地震资料建立裂缝阻抗模型；根据所述初始低频模型和所述裂缝阻抗模型，得到混合阻抗低频模型；根据所述混合阻抗低频模型进行反演，确定所述目标区域的纵波阻抗数据体，以基于所述纵波阻抗数据体确定所述目标区域的空间展布；获取表示纵波阻抗数据体与孔隙度之间关系的关系函数，以根据所述关系函数基于所述空间展布确定所述裂缝储层。根据目标区域的地震资料及测井资料建立初始低频模型，同时根据地震资料建立裂缝阻抗模型，然后根据初始低频模型和裂缝阻抗模型得到的混合阻抗低频模型，从而能够根据混合阻抗低频模型在纵向和横向的立体空间上刻画储层展布。对该混合阻抗低频模型进行反演得到目标区域的纵波阻抗数据体，纵波阻抗数据体能够反映目标区域的岩性、物性等信息。因此，根据
纵波阻抗数据体及纵波阻抗数据体与孔隙度之间关系的关系函数则可以准确的对目标区域内部缝洞体的空间展布进行刻画，从而能够更直观的刻画储层的缝洞连通关系。
205.在本技术实施例所提供的几个实施例中，应该理解到，所揭露的系统和方法，也可以通过其它的方式实现。以上所描述的系统和方法实施例仅仅是示意性的。
206.需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个
……”
限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
207.虽然本技术所揭露的实施方式如上，但所述的内容只是为了便于理解本技术而采用的实施方式，并非用以限定本技术。任何本技术所属技术领域内的技术人员，在不脱离本技术所揭露的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本技术的专利保护范围，仍须以所附的权利要求书所界定的范围为准。