碳酸盐岩储层渗透率的计算方法及计算装置、设备和介质与流程

本发明涉及测井，具体涉及一种碳酸盐岩储层渗透率的计算方法及计算装置、设备和介质，尤其涉及一种基于电成像测井孔隙度谱的碳酸盐岩储层渗透率的计算方法及计算装置、设备和介质。

背景技术：

1、世界上比较大型的、高产的油气田多来自于碳酸盐岩储集层，其在地壳中的分布范围仅次于碎屑岩，其中中东地区碳酸盐岩油气藏的储量又占全球总储量的70％。受到沉积过程、成岩作用及构造运动等方面的影响，致使碳酸盐岩储层具有岩性复杂、孔隙类型多样、各向异性强以及非均质性强等特点，而这些特点使得储层的测井响应比较复杂，给测井评价带来较大的困难，尤其使碳酸盐岩地层的渗透率评价研究成为测井评价的难题之一。

2、目前，评价碳酸盐岩储层渗透率方法主要分为基于常规测井资料和基于特殊测井资料两种思路。最直接的方法是分析岩心渗透率与岩心孔隙度或者测井响应参数之间的相关性，进而建立岩心刻度统计模型；或应用经典的岩石物理统计模型建立渗透率与储层特征参数之间的关系。在碳酸盐岩储集层中，分析发现，岩心物性实验数据的岩心渗透率与孔隙度之间的相关性较差，且应用常规测井模型预测的地层渗透率误差大、精度低，说明利用常规测井方法计算的储层渗透率多为受基质孔隙影响较大的基质渗透率，这类方法预测的渗透率不能表征次生孔隙对地层渗透率的贡献，所以计算的渗透率不能满足地层渗透率评价的精度。

3、近年来，对基于特殊测井资料的储层渗透率评价方法开展了相应的研究，主要有基于核磁共振测井、电成像测井和阵列声波测井资料的储层渗透率评价方法。电成像测井是近十几年发展了起来的新兴储层测井技术，它是利用井眼测得的大量井壁高分辨率图像和地层电阻率的参数测量值，做成一幅覆盖整个井壁的二维图像。电成像测井资料具有高分辨率、能够直观地观察地层构造特征和孔隙、孔洞空间特征等特点，因此，可以通过电成像测井资料获取地层中有效孔隙的信息，进而利用这些信息直接观察和分析渗透率的影响因素，从而定量识别碳酸盐岩储层中孔隙空间的分布状况及其内部的几何结构，对碳酸盐岩储层测井评价具有重要的意义，对储层渗透率评价有非常重要的应用价值。

4、目前，本领域已出现多种使用成像测井资料进行碳酸盐岩储层渗透率计算的方法，如秦瑞宝等(潜山油田裂缝孔隙度和渗透率测井评价新方法[j].中国海上油气,2015.doi:cnki:sun:zhsd.0.2015-03-005.)通过电成像测井孔隙度谱类似于核磁测井构建伪毛管压力曲线，从而利用毛管压力曲线提取swanson参数，建立了利用电成像测井预测渗透率的新方法；cn108828687a公开了通过成像测井得到孔洞面孔率，基于地层孔隙结构以及孔洞面孔率，通过地层有效渗透率公式计算渗透率；cn112098291a公开了一种基于电成像孔隙度谱的渗透率计算方法，针对已知岩心的电成像测井资料处理获得孔隙度谱，进而将孔隙度谱进行反向累加，构建伪毛管压力曲线，从而根据伪毛管压力曲线获得孔隙结构参数，通过将孔隙结构参数与已知岩心的渗透率进行拟合，构建孔隙结构-渗透率模型。

5、cn117345220a公开了一种根据孔隙直径谱分别计算碳酸盐岩储层中微孔、细孔、中孔、粗孔、洞五种不同孔隙组分的孔隙度，进而将根据各孔隙组分的孔隙度与渗透率之间的关系计算的渗透率求和，获得储层的最终渗透率。

6、综上可知，针对电成像测井渗透率解释方面从不同角度开展了大量的研究和工作，但由于不同地区储层岩石背景的差异和电成像测井测量仪器的不同，几种方法在渗透率解释精度方面虽有一定的提升，但从应用的推广程度来看存在一定的局限性。同时，碳酸盐岩储层渗透率与储层孔隙类型以及孔隙结构有着密切关系。由此可见，如何提供一种基于成像测井资料并综合考虑储层孔隙类型与孔隙结构的方法，提高碳酸盐岩储层渗透率计算精度，成为了目前本领域技术人员迫切需要解决的问题。

技术实现思路

1、鉴于现有技术中存在的问题，本发明的目的在于提供一种碳酸盐岩储层渗透率的计算方法及计算装置、设备和介质，以解决现有碳酸盐岩储层渗透率计算精度差，适应性差的缺陷。

2、为达此目的，本发明采用以下技术方案：

3、第一方面，本发明提供了一种碳酸盐岩储层渗透率的计算方法，所述计算方法包括：

4、构建电成像的绝对刻度孔隙度谱和浮动刻度孔隙度谱，之后结合孔分布占比特征划分孔隙结构类型；

5、提取各深度点对应的绝对孔隙度谱的特征参数，之后结合基于成像测井图所得面孔率构建不同孔隙结构类型的渗透率模型；

6、依据实际工区岩心渗透率标定所得渗透率模型的待定系数，得到标定渗透率模型，进而计算非取芯井渗透率。

7、本发明利用电成像测井资料具有较高储层纵向分辨率、能直观地反映储层中溶蚀孔洞信息的优势，通过构建电成像绝对与浮动刻度孔隙度谱，并根据提取的谱特征参数划分不同孔隙结构类型。考虑碳酸盐岩储层的强非均质性，以及渗透率与储层孔隙类型和孔隙结构之间的密切关系，利用提取的孔隙度谱特征参数并结合面孔率，构建了一种基于电成像孔隙度谱以及孔隙结构类型的多参数渗透率模型，该方法能够突显碳酸盐岩储层中不同孔隙大小及孔隙结构类型对渗透率的影响。

8、作为本发明优选的技术方案，所述构建电成像的绝对刻度孔隙度谱和浮动刻度孔隙度谱包括构建绝对刻度孔隙度谱和构建浮动刻度孔隙度谱。

9、作为本发明优选的技术方案，所述构建绝对刻度孔隙度谱包括利用计算各深度段内绝对孔隙度谱的最大值和最小值作为原始孔隙度谱的左右刻度值。

10、其中，所述构建浮动刻度孔隙度谱包括利用各深度点绝对孔隙度谱的最大值和最小值作为各深度点处原始孔隙度谱的左右刻度值，重新进行数据排列，放大谱峰形态。

11、作为本发明优选的技术方案，所述孔分布占比特征包括依据孔径大小划分大孔、中孔和小孔；

12、其中，所述大孔为孔直径为dl≥500μm；

13、其中，所述中孔的孔直径为100μm≤dm＜500μm；

14、其中，所述小孔的孔直径为ds＜100μm。

15、作为本发明优选的技术方案，所述特征参数包括对应深度点处电成像最大孔隙度对应深度点处电成像最小孔隙度对应深度点处电成像孔隙度平均值对应深度点电成像孔隙度谱几何均值pigm。

16、作为本发明优选的技术方案，所述结合基于成像测井图所得面孔率构建不同孔隙结构类型的渗透率模型中所得渗透率模型的表达式为：

17、

18、式中，ai为待定系数；bi+ci＝1且bi≥ci且ci≠0；若大孔占主导，则若大孔和中孔占主导，则和若中孔占主导，则若小孔和中孔占主导，则和若小孔占主导，则为电成像某一深度点绝对孔隙度最大值；为电成像某一深度点绝对孔隙度最小值；为电成像某一深度点绝对孔隙度均值；pigm为电成像某一深度点绝对孔隙度谱几何均值，数值范围为0-200；fiv为某一深度点面孔率；i取值为1,2,3,4……n，n为正整数；j取值为1,2,3,4……n，n为正整数。

19、作为本发明优选的技术方案，所述依据实际工区岩心渗透率标定所得渗透率模型的待定系数包括ai；

20、其中，bi和ci根据计算中渗透率和岩心渗透率的误差进行调整。

21、第二方面，本发明提供了碳酸盐岩储层渗透率的计算装置，所述计算装置包括：

22、孔隙结构类型获取模块，用于构建电成像的绝对刻度孔隙度谱和浮动刻度孔隙度谱，之后结合孔分布占比特征划分孔隙结构类型；

23、渗透率模型获取模块，用于提取各深度点对应的绝对孔隙度谱的特征参数，之后结合基于成像测井图所得面孔率构建不同孔隙结构类型的渗透率模型；

24、非取芯井渗透率获取模块，用于依据实际工区岩心渗透率标定所得渗透率模型的待定系数，得到标定渗透率模型，进而计算非取芯井渗透率。

25、第三方面，本发明提供了一种电子设备，所述电子设备包括：

26、至少一个处理器；以及与所述至少一个处理器通信连接的存储器；

27、其中，所述存储器存储有可被所述至少一个处理器执行的计算机程序，所述计算机程序被所述至少一个处理器执行，以使所述至少一个处理器能够执行第一方面所述碳酸盐岩储层渗透率的计算方法。

28、第四方面，本发明提供了一种计算机存储介质，所述计算机存储介质中存储有计算机可执行指令，所述计算机可执行指令被处理器执行时实现第一方面所述碳酸盐岩储层渗透率的计算方法。

29、与现有技术方案相比，本发明具有以下有益效果：

30、本发明提供的计算方法，综合考虑了影响碳酸盐岩储层渗透率参数的储层孔隙类型与孔隙结构等因素，具有渗透率解释精度高的优点，在实际应用中取得了良好的效果，对后续储层评价及三维空间表征具有指导意义。