一种基于成像测井的孔隙度刻度方法与流程

本发明涉及油气勘探，更具体的说，它涉及一种基于成像测井的孔隙度刻度方法。

背景技术：

1、对缝洞型碳酸盐岩储层、裂缝性致密砂岩储层、裂缝性火山岩储层等次生裂缝发育的致密岩石储层，裂缝发育情况是度量此类储层好坏的关键因素。

2、在裂缝性储层中，定量计算裂缝参数是测井精细解释的基础，其中最重要的是定量计算裂缝孔隙参数。裂缝孔隙度和裂缝面孔率是最常用的裂缝孔隙参数,裂缝孔隙度为裂缝孔隙体积与岩石体积之比，描述的是裂缝在三维立体空间发育情况；裂缝面孔率为单位井段上裂缝在井壁上所占面积与微电阻率成像测井覆盖井壁的面积之比，描述的是裂缝在二维平面发育的情况。

3、中国专利cn106014397a，提出了一种成像测井缝洞储层参数提取方法，包括以下步骤：准确提取孔洞、裂缝图像；输入经浅电阻率刻度后的成像测井钮扣电极曲线段；对成像测井钮扣电极曲线段进行一维二进小波变换分解；选取对孔洞、裂缝敏感的一维二进小波变换谱阶；寻找已选择的阶小波变换谱所有的峰谷点，并根据峰谷点之间小波变换谱值的大小对这些峰谷点进行剔除；保留筛选下来的像素值；并绘制成成像图；除掉成像图中背景岩石的测量值；运行单目标边缘跟踪算法，并在其执行的过程中分别计算出单目标的周长、面积、长度及宽度数据。

4、上述专利是从电成像电导率图像上分割缝洞子图像、提取缝洞参数的成像图像分割和缝洞参数，应用相应的技术从电成像资料上提取的是缝洞面孔率。但不论是在数量级上还是具体的物理含义上，电成像资料上提取的缝洞面孔率与岩石缝洞孔隙度均有较大差异，缝洞面孔率无法准确描述裂缝在三维立体空间发育情况。

5、为此，本发明旨在于提供一种基于成像测井的孔隙度刻度方法，利用岩心ct扫描数字图像提取的缝洞孔隙度将成像测井提取的缝洞面孔率刻度成缝洞孔隙度，为更好地综合应用面孔率和孔隙度评价次生裂缝发育的致密岩石储层提供依据。

技术实现思路

1、针对现有技术存在的不足，本发明的目的在于提供一种基于成像测井的孔隙度刻度方法，其根据岩心ct扫描数字图像提取的缝洞孔隙度和通过成像测井图像提取的缝洞面孔率拟合出缝洞孔隙度与缝洞面孔率之间的统计关系式，并根据获得的统计关系式将成像测井提取的缝洞面孔率刻度成缝洞孔隙度，为更好地综合应用面孔率和孔隙度评价次生裂缝发育的致密岩石储层提供了依据。

2、为实现上述目的，本发明提供了如下技术方案：一种基于成像测井的孔隙度刻度方法，具体步骤如下：

3、s1、根据岩心ct扫描数字图像，获取岩心ct缝洞孔隙度；

4、s2、根据成像测井图像，获取成像缝洞面孔率；步骤s1、s2的顺序可调换；

5、s3、利用岩心ct缝洞孔隙度将成像缝洞面孔率刻度成缝洞孔隙度；

6、步骤s3具体如下：

7、s31、岩心归位：获取岩心ct扫描数字图像对应测井深度的成像缝洞面孔率；

8、s32、以缝洞孔隙度和缝洞面孔率为横纵坐标建立坐标系，将步骤s1获取的岩心ct缝洞孔隙度和步骤s31获取的成像缝洞面孔率按照对应关系填入坐标系中构成数据点；

9、s33、根据步骤s32填入坐标系中的数据点拟合出缝洞孔隙度与缝洞面孔率之间的统计关系式；

10、s34、根据步骤s33获得的统计关系式，将步骤s2获得的成像缝洞面孔率刻度转换成缝洞孔隙度。

11、本发明进一步设置为：步骤s33中的统计关系式为三次多项式。

12、本发明进一步设置为：所述步骤s1具体步骤如下：

13、s11、输入选定的岩心ct扫描数字图像，对图像进行降噪处理；

14、s12、对降噪处理后的图像进行图像分割处理，获得缝洞二值化图像；

15、s13、统计缝洞二值化图像中像素灰度值为0的个数即为该岩心切片中的缝洞孔隙面积s缝洞，灰度值为255的个数即为该岩心切片中的岩石骨架的面积s骨架；

16、s14、判断处理的岩心ct扫描数字图像是否是最后一个岩心切片，如果不是最后一个岩心切片，则移动至下一个岩心切片重复s11-s13；如果是最后一个岩心切片，则按公式计算岩心ct缝洞孔隙度。

17、本发明进一步设置为：步骤s11的降噪处理采用中值滤波法。

18、本发明进一步设置为：步骤s12的图像分离处理采用区域生长法。

19、本发明进一步设置为：所述步骤s2具体步骤如下：

20、s21、输入一个图像框的电导率图像数据，对输入的电导率图像数据进行图像分割，得到裂缝子图像和孔洞子图像；

21、s22、给定步长，按步长移动图像框重复步骤s21，直至测井终点，并输出裂缝子图像和孔洞子图像；

22、s23、输入一个图像框的裂缝子图像数据和孔洞子图像数据，按如下公式计算裂缝面孔率和孔洞面孔率：

23、

24、

25、s24、给定步长，按步长移动图像框重复步骤s23，直至测井终点深度，输出裂缝面孔率曲线和孔洞面孔率曲线，所述裂缝面孔率曲线和孔洞面孔率曲线即为成像缝洞面孔率。

26、本发明进一步设置为：步骤s21中图像分割采用二维小波变换法。

27、本发明进一步设置为：步骤s21中图像框的窗长为256个采样点；步骤s22中的步长为256个采样点。

28、本发明进一步设置为：步骤s23中图像框的窗长为96个采样点，步骤s24中的步长为48个采样点。

29、综上所述，本发明相比于现有技术具有以下有益效果：本发明根据岩心ct扫描数字图像提取的缝洞孔隙度和通过成像测井图像提取的缝洞面孔率拟合出缝洞孔隙度与缝洞面孔率之间的统计关系式，并根据获得的统计关系式将成像测井提取的缝洞面孔率刻度成缝洞孔隙度，为更好地综合应用面孔率和孔隙度评价次生裂缝发育的致密岩石储层提供了依据。

技术特征：

1.一种基于成像测井的孔隙度刻度方法，其特征在于，具体步骤如下：

2.根据权利要求1所述的一种基于成像测井的孔隙度刻度方法，其特征在于：步骤s33中的统计关系式为三次多项式。

3.根据权利要求1-2任一所述的一种基于成像测井的孔隙度刻度方法，其特征在于，所述步骤s1具体步骤如下：

4.根据权利要求3所述的一种基于成像测井的孔隙度刻度方法，其特征在于：步骤s11的降噪处理采用中值滤波法。

5.根据权利要求4所述的一种基于成像测井的孔隙度刻度方法，其特征在于：步骤s12的图像分离处理采用区域生长法。

6.根据权利要求1-2任一所述的一种基于成像测井的孔隙度刻度方法，其特征在于，所述步骤s2具体步骤如下：

7.根据权利要求6所述的一种基于成像测井的孔隙度刻度方法，其特征在于：步骤s21中图像分割采用二维小波变换法。

8.根据权利要求7所述的一种基于成像测井的孔隙度刻度方法，其特征在于：步骤s21中图像框的窗长为256个采样点；步骤s22中的步长为256个采样点。

9.根据权利要求7所述的一种基于成像测井的孔隙度刻度方法，其特征在于：步骤s23中图像框的窗长为96个采样点，步骤s24中的步长为48个采样点。

技术总结本发明公开了一种基于成像测井的孔隙度刻度方法，其技术方案要点包括S1、获取岩心CT缝洞孔隙度；S2、获取成像缝洞面孔率；S3、获取岩心CT扫描数字图像对应测井深度的成像缝洞面孔率；根据步骤S1获取的岩心CT缝洞孔隙度和成像缝洞面孔率拟合统计关系式；根据统计关系式，将成像缝洞面孔率刻度转换成缝洞孔隙度。本发明根据岩心CT扫描数字图像提取的缝洞孔隙度和通过成像测井图像提取的缝洞面孔率拟合出缝洞孔隙度与缝洞面孔率之间的统计关系式，并根据获得的统计关系式将成像测井提取的缝洞面孔率刻度成缝洞孔隙度，为更好地综合应用面孔率和孔隙度评价次生裂缝发育的致密岩石储层提供了依据。技术研发人员：蔺学旻,孙兆辉,王谦,周威,秦伟强,赵冉受保护的技术使用者：中国石油化工股份有限公司技术研发日：技术公布日：2024/7/23