电成像缝洞参数的提取方法、系统及计算机可读存储介质与流程

1.本技术实施例涉及测井技术，特别涉及电成像缝洞参数的提取方法、系统及计算机可读存储介质。


背景技术：

2.碳酸盐岩缝洞储层具有基块孔隙度低，裂缝、溶蚀孔洞等次生储渗空间发育的特点。对于这类裂缝、孔洞发育的碳酸盐岩缝洞储层，裂缝、溶蚀孔洞的发育和组合情况是控制这类储层有效性的主要因素。因此，如何从储层测井资料中提取裂缝、溶蚀孔洞信息成为评价这类储层有效性即储层能否产出工业油流或经储层改造后能否产出工业油流的关键问题。电成像测井通过多纽扣电极测量井壁附近地层的电导率曲线，经过预处理，结合伴测的井斜、井径及方位测量数据将电导率曲线以图像的形式直观地表达出来。由于导电泥浆对裂缝、溶蚀孔洞的侵入，裂缝、溶蚀孔洞处的导电性与基块岩石的导电性有差别，因此，可以从电成像测井图像中提取裂缝、溶蚀孔洞信息来评价碳酸盐岩缝洞储层的有效性。
3.相关技术中，对电成像缝洞参数的提取以小波变换模极大值图像分割技术为基础。具体地，以纽扣电池电极电导率曲线为对象，先消除井壁凹凸不平导致的地层背景噪声的影响，利用小波变换模极大值图像分割方法得到包含裂缝和溶蚀孔洞目标的子图像，根据子图像提取缝洞参数信息。
4.然而，相关技术对电成像缝洞参数的提取方法存在的问题是：1、电成像测井仪的纽扣电极极板并未对井眼形成全覆盖，因而在极板间存在空白带，没有测量信息。由于这种极板间空白带的存在，相关技术对电成像资料的处理并不清晰完整，影响了电成像缝洞参数提取结果的准确性。2、图像分割出来的高电导目标中除包含有裂缝和溶蚀孔洞之外，还包含有水平状高电导泥质纹层、缝合线及诱导缝等无效信息。相关技术将所有分割出来的高电导目标都当作裂缝、溶蚀孔洞对待，因而在提取的电成像缝洞参数中包含有这些无效信息。


技术实现要素：

5.本技术实施例提供了电成像缝洞参数的提取方法、系统及计算机可读存储介质，可用于完整地获得电成像缝洞参数提取结果，并提高电成像缝洞参数提取结果的准确性。
6.一方面，本技术实施例提供了电成像缝洞参数的提取方法，该方法包括：获取电成像测井全井眼覆盖图像；对电成像测井全井眼覆盖图像进行初分割，获取全井眼电导率图像初分割缝洞子图像；提取全井眼电导率图像初分割缝洞子图像的目标形状参数；基于全井眼电导率图像初分割缝洞子图像的目标形状参数，对全井眼电导率图像初分割缝洞子图像进行二次分割，获取裂缝、孔洞子图像；基于裂缝、孔洞子图像，提取电成像缝洞参数。
7.可选地，获取电成像测井全井眼覆盖图像，包括：获取电成像测井全井眼电导率数据；基于电成像测井全井眼电导率数据，获取电成像测井全井眼覆盖图像。
8.可选地，获取电成像测井全井眼电导率数据，包括：输入原始电导率数据；基于原
始电导率数据，绘制原始电成像电导率图像；基于原始电成像电导率图像，应用电成像极板间空白充填的方法，获取电成像测井全井眼电导率数据。
9.可选地，对电成像测井全井眼覆盖图像进行初分割，获取全井眼电导率图像初分割缝洞子图像，包括：基于电成像测井全井眼覆盖图像，获取一个图像框所在区域的全井眼电导率图像初分割缝洞子图像；移动图像框，获取其他区域的全井眼电导率图像初分割缝洞子图像；基于一个图像框所在区域的全井眼电导率图像初分割缝洞子图像以及其他区域的全井眼电导率图像初分割缝洞子图像，输出全井眼电导率图像初分割缝洞子图像。
10.可选地，基于电成像测井全井眼覆盖图像，获取一个图像框所在区域的全井眼电导率图像初分割缝洞子图像，包括：输入一个图像框的电成像测井全井眼覆盖图像数据；对一个图像框的电成像测井全井眼覆盖图像数据进行二维二进小波变换，对于阶数为2的小波变换计算每个电导率数据点小波变换系数的模值，得到2阶小波变换系数模值；基于2阶小波变换系数模值，在深度方向上利用三点寻峰获取模极大值点；基于模极大值点，进行图像初分割，获取一个图像框所在区域的全井眼电导率图像分割缝洞子图像。
11.可选地，提取全井眼电导率图像初分割缝洞子图像的目标形状参数，包括：基于全井眼电导率图像初分割缝洞子图像，提取高电导目标的角度参数，作为全井眼电导率图像初分割缝洞子图像的目标形状参数。
12.可选地，基于全井眼电导率图像初分割缝洞子图像的目标形状参数，对全井眼电导率图像初分割缝洞子图像进行二次分割，获取裂缝、孔洞子图像，包括：在全井眼电导率图像初分割缝洞子图像上，去掉目标形状参数小于有效阈值的非有效目标，得到裂缝、孔洞子图像。
13.可选地，基于裂缝、孔洞子图像，提取电成像缝洞参数，包括：基于裂缝、孔洞子图像，对一个图像框所在区域提取电成像缝洞参数；移动图像框，对其他区域提取电成像缝洞参数。
14.一方面，本技术实施例提供了电成像缝洞参数的提取系统，该系统包括：全井眼覆盖图像生成模块，用于获取电成像测井全井眼覆盖图像；初分割模块，用于对电成像测井全井眼覆盖图像进行初分割，获取全井眼电导率图像初分割缝洞子图像；目标形状参数提取模块，用于提取所述全井眼电导率图像初分割缝洞子图像的目标形状参数；二次分割模块，用于基于全井眼电导率图像初分割缝洞子图像的目标形状参数，对全井眼电导率图像初分割缝洞子图像进行二次分割，获取裂缝、孔洞子图像；缝洞参数提取模块，用于基于裂缝、孔洞子图像，提取电成像缝洞参数。
15.本技术实施例还提供了一种计算机可读存储介质，该计算机可读存储介质中有至少一条程序指令或代码，所述程序指令或代码由处理器加载并执行时以使计算机实现如上述的电成像缝洞参数的提取方法。
16.本技术实施例提供的技术方案至少包括以下有益效果：1、以完整的二维方式(全井眼电导率数据)进行图像分割，并以目标形状参数为基础进行二次分割，分割出来的缝洞子图像更为清晰完整。2、分割出来的高电导目标剔除了水平状高电导泥质纹层、缝合线及诱导缝等无效信息，使电成像缝洞参数提取结果更为准确。
附图说明
17.为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
18.图1是本发明实施例提供的一种电成像缝洞参数提取方法的流程图；
19.图2是本发明实施例提供的一种电成像缝洞参数提取方法所获取的原始电成像电导率图像以及电成像测井全井眼覆盖图像；
20.图3是本发明实施例提供的一种电成像缝洞参数提取方法所获取的原始电成像电导率图像、电成像测井全井眼覆盖图像以及全井眼电导率图像初分割缝洞子图像；
21.图4是本发明实施例提供的一种电成像缝洞参数提取方法所提取的目标形状参数的示意图；
22.图5是本发明实施例提供的一种电成像缝洞参数提取方法所获取的原始电成像电导率图像、电成像测井全井眼覆盖图像、全井眼电导率图像初分割缝洞子图像以及裂缝、孔洞子图像；
23.图6是本发明实施例提供的一种电成像缝洞参数提取方法所获取的电成像测井全井眼覆盖图像、全井眼电导率图像初分割缝洞子图像、裂缝、孔洞子图像以及电成像缝洞参数的提取结果；
24.图7是本发明实施例提供的一种电成像缝洞参数提取方法所获取的电成像缝洞参数的提取结果与常规测井资料对比结果；
25.图8是本发明实施例提供的一种电成像缝洞参数提取方法所获取的电成像缝洞参数的提取结果与相关技术所获取的缝洞参数的提取结果的对比；
26.图9是本发明实施例提供的一种电成像缝洞参数提取方法在techlog环境中实现的流程图；
27.图10是本发明实施例提供的一种电成像缝洞参数的提取系统的示意图。
具体实施方式
28.为使本技术的目的、技术方案和优点更加清楚，下面将结合附图对本技术实施方式作进一步地详细描述。
29.请参考图1，其示出了本技术一个实施例提供的一种电成像缝洞参数提取方法的流程图，该方法包括如下几个步骤：
30.步骤101，获取电成像测井全井眼覆盖图像。
31.相关技术中，电成像测井仪的纽扣电极极板并未对井眼形成100％的覆盖，因此在极板间存在空白带，没有测量信息。由于这种极板间空白带的存在，沿井壁正北方向剖开的电成像资料就不能以完整的二维方式进行处理。基于此，本技术实施例首先获取电成像测井全井眼覆盖图像。
32.获取电成像测井全井眼覆盖图像，包括但不限于如下几个子步骤：
33.1011、获取电成像测井全井眼电导率数据。
34.获取电成像测井全井眼电导率数据，包括但不限于如下子步骤：
35.1011
‑①
输入原始电导率数据。
36.其中，输入的原始电导率数据为经浅电阻率刻度后的极板电导率数据。由于电成像测井仪的纽扣电极系是非聚集电极系，测量值仅与井壁附近地质体的电导率成比例变化，因此要用浅电阻率进行刻度，得到经浅电阻率刻度后的极板电导率数据。
37.1011
‑②
基于原始电导率数据，绘制原始电成像电导率图像。
38.基于原始电导率数据，根据参考极板中间纽扣电极电导率曲线的方位，沿井筒正北剖开，绘制电成像电导率图像。
39.参考极板是在测井过程中为其他极板提供方位参考的极板。可选地，参考极板可以是1号极板，1号极板是电成像测井仪器中编号为1的极板，测井过程中测量并记录1号极板的方位，即1号极板中点相对于正北方向的方位，将1号极板的方位为其他极板提供方位参考。在本领域的绘图处理中，通常要根据1号极板的方位确定其它极板相对于正北方向的方位。在根据1号极板中间纽扣电极电导率曲线的方位绘制电成像测井的电导率图像时，需要沿井筒的正北方向剖开进行绘制。
40.可选地，不同电成像测井仪不同参考极板中间纽扣电极不同，如增强型微电阻率成像仪(extended range micro-resistivity imager，xrmi)将1号极板作为参考极板，xrmi中1号极板第13个纽扣电极为中间纽扣电极；地层微电阻率扫描成像仪(formation microscanner imager，fmi)将1号极板作为参考极板，fmi中1号极板第12个纽扣电极为中间纽扣电极。
41.1011
‑③
基于原始电成像电导率图像，应用电成像极板间空白充填的方法，获取电成像测井全井眼电导率数据。
42.可选地，电成像极板间空白充填的方法指的是基于原始电成像电导率图像，应用极板边缘向空隙中间生长的方法，生成电成像测井仪极板间空白带数据，从而得到电成像测井全井眼覆盖的电导率数据，即电成像测井全井眼电导率数据。
43.1012、基于电成像测井全井眼电导率数据，获取电成像测井全井眼覆盖图像。
44.在一种可能的实现方法中，根据电成像测井全井眼电导率数据，利用参考极板方位数据计算每个极板相对正北方向的方位，绘制得到极板间空白带充填后的全井眼覆盖的电导率图像，即电成像测井全井眼覆盖图像。
45.请参考图2，其示出了本发明的一个实施例提供的一种电成像缝洞参数提取方法所获取的原始电成像电导率图像以及电成像测井全井眼覆盖图像。在一种可能的实现方法中，根据fmi电成像测井仪1号极板中间纽扣电极电导率曲线的方位，沿井筒正北剖开，绘制出了如图2所示的沿井眼正北剖开的电导率图像，即原始电成像电导率图像。
46.可选地，根据电成像测井全井眼电导率数据，利用fmi电成像测井仪1号极板方位数据计算每个极板相对正北方向的方位，绘制得到如图2所示的全井眼覆盖的电导率图像，即电成像测井全井眼覆盖图像。
47.其中，图2第一道为沿井眼正北剖开的电导率图像，即原始电成像电导率图像；第二道为深度道，表示测井深度；第三道为全井眼覆盖的电导率图像，即电成像测井全井眼覆盖图像。图2中，代号p1az表示1号极板方位，deg表示度，为角度单位，heated表示图头，fmi_256表示沿井眼正北剖开的电导率图像，即原始电成像电导率图像，fmi_fill_liu3_032_256表示电成像测井全井眼覆盖图像，md代表测量深度。
48.步骤102，对电成像测井全井眼覆盖图像进行初分割，获取全井眼电导率图像初分割缝洞子图像。
49.对电成像测井全井眼覆盖图像进行初分割，获取全井眼电导率图像初分割缝洞子图像，包括但不限于如下几个子步骤：
50.1021、基于电成像测井全井眼覆盖图像，获取一个图像框所在区域的全井眼电导率图像初分割缝洞子图像。
51.其中，一个图像框相当于常规测井资料的采样间隔，一个图像框所在区域是指一个图像框在电成像测井全井眼覆盖图像上所围成的区域。本技术实施例首先对一个图像框所在区域的电成像测井全井眼覆盖图像进行处理。
52.可选地，基于电成像测井全井眼覆盖图像，获取一个图像框所在区域的全井眼电导率图像初分割缝洞子图像包括但不限于如下子步骤：
53.1021
‑①
输入一个图像框的电成像测井全井眼覆盖图像数据。
54.1021
‑②
对一个图像框的电成像测井全井眼覆盖图像数据进行二维二进小波变换，对于阶数为2的小波变换计算每个电导率数据点小波变换系数的模值，得到2阶小波变换系数模值。
55.在一个示例中，将一个图像框的电成像测井全井眼覆盖图像中的电导率数据设为:
[0056][0057]
进行二次小波变换，得到该电导率数据对应的小波系数:
[0058][0059]
其中，j为小波变换阶数，取电导率数据2阶小波变换的小波系数：
[0060][0061]
计算2阶小波变换系数的模值：
[0062][0063]
1021
‑③
基于2阶小波变换系数模值，在深度方向上利用三点寻峰法获取模极大值点。
[0064]
三点寻峰法是一种寻找极大值点的方法，利用深度方向上曲线的三个点拟合一个二次多项式，再求出二次多项式的极大值点及对应的坐标，该点即为模极大值点。
[0065]
1021
‑④
基于模极大值点，获取模极大值点的坐标像素值，根据模极大值点的坐标像素值确定分割阈值进行图像初分割，获取一个图像框所在区域的全井眼电导率图像分割缝洞子图像。
[0066]
模极大值点是全井眼电导率图像电导率值突变点，即分割点。可选地，根据分割点位置确定分割阈值，可选地，根据分割阈值进行图像初分割。
[0067]
1022、移动图像框，获取其他区域的全井眼电导率图像初分割缝洞子图像。
[0068]
在一个可能的实现方式中，移动图像框，获取其他区域的全井眼电导率图像初分割缝洞子图像包括：移动图像框，重复上述子步骤1021的方法，直至测井终点深度，得到其他区域的全井眼电导率图像初分割缝洞子图像。其中，测井终点深度从图像中获得。
[0069]
1023、基于一个图像框所在区域的全井眼电导率图像初分割缝洞子图像以及其他
区域的全井眼电导率图像初分割缝洞子图像，输出全井眼电导率图像初分割缝洞子图像。
[0070]
请参考图3，其示出了本发明的一个实施例提供的一种电成像缝洞参数提取方法所获取的原始电成像电导率图像、电成像测井全井眼覆盖图像以及全井眼电导率图像初分割缝洞子图像。该实施例的电成像测井仪为fmi测井仪，参考极板为1号极板。
[0071]
其中，图3第一道为沿井眼正北剖开的电导率图像，即原始电成像电导率图像；第二、四为深度道，表示测井深度；第三道为全井眼覆盖的电导率图像，即电成像测井全井眼覆盖图像；第五道为二维二进小波变换图像初分割结果图，即全井眼电导率图像初分割缝洞子图像。图3中，代号p1az表示1号极板方位，deg表示度，为角度单位，heated表示图头，fmi_256表示沿井眼正北方向剖开的电导率图像，fmi_fill_liu3_032_256表示电成像测井全井眼覆盖图像，md代表测量深度，fmi_fill_liu3_032_256_seg_fmi_v表示二维二进小波变换初分割缝洞子图像，即全井眼电导率图像初分割缝洞子图像。
[0072]
步骤103，提取全井眼电导率图像初分割缝洞子图像的目标形状参数。
[0073]
在一个可能的实现方式中，提取全井眼电导率图像初分割缝洞子图像的目标形状参数包括：基于全井眼电导率图像初分割缝洞子图像，提取高电导目标的角度参数，将高电导目标的角度参数作为全井眼电导率图像初分割缝洞子图像的目标形状参数。
[0074]
在剖开的全井眼覆盖的二维电导率图像上，水平状高电导泥质纹层、缝合线等非有效高电导目标与裂缝、孔洞目标的角度不同。为了应用初步分割目标形状参数区分非有效高电导目标与裂缝、孔洞目标，对初分割出来的缝洞子图像进行进一步处理，提取高电导目标的角度参数。
[0075]
请参考图4，其示出了本发明的一个实施例提供的一种电成像缝洞参数提取方法所提取的目标形状参数的示意图。
[0076]
在一种可能的实现方法中，提取全井眼电导率图像初分割缝洞子图像的目标形状参数包括：基于全井眼电导率图像初分割缝洞子图像，提取高电导目标的角度参数。
[0077]
可选地，高电导目标的角度参数定义为高电导目标长轴，即目标边缘上距离最大的两点之间的连线，与水平方向的夹角。进一步可选地，按照如下公式定义角度参数：
[0078][0079]
其中，水平方向为x轴，高电导目标长轴为y轴，(x2，y2)为高电导目标长轴的终点，(x1，y1)为高电导目标长轴的起点。
[0080]
步骤104，基于全井眼电导率图像初分割缝洞子图像的目标形状参数，对全井眼电导率图像初分割缝洞子图像进行二次分割，获取裂缝、孔洞子图像。
[0081]
在一种可能的实现方式中，基于全井眼电导率图像初分割缝洞子图像的目标形状参数，对全井眼电导率图像初分割缝洞子图像进行二次分割，获取裂缝、孔洞子图像，包括：在全井眼电导率图像初分割缝洞子图像上，去掉目标形状参数小于有效阈值的非有效目标，得到裂缝、孔洞子图像。
[0082]
在直井中，例如水平状或倾角较低的高电导泥质纹层、缝合线等非有效目标的角度参数很小。示例性地，以10
°
作为目标形状参数的有效阈值。进一步地，将目标形状参数小于10
°
，即角度参数小于10
°
的高电导目标作为非有效目标，对全井眼电导率图像初分割缝洞子图像进行二次分割，去掉非有效目标，得到裂缝、孔洞子图像。其中，有效阈值可基于经
验设置，本技术实施例对此不进行限定。
[0083]
请参考图5，其示出了本发明的一个实施例提供的一种电成像缝洞参数提取方法所获取的原始电成像电导率图像、电成像测井全井眼覆盖图像、全井眼电导率图像初分割缝洞子图像以及裂缝、孔洞子图像。该实施例的电成像测井仪为fmi测井仪，参考极板为1号极板。
[0084]
其中，图5第一道为沿井眼正北剖开的电导率图像，即原始电成像电导率图像；第二、四、六为深度道，表示测井深度；第三道为全井眼覆盖的电导率图像，即电成像测井全井眼覆盖图像；第五道为二维二进小波变换图像初分割结果图，即全井眼电导率图像初分割缝洞子图像；第七道为以目标形状参数为基础的二次分割结果图，即裂缝、孔洞子图像。可见，全井眼电导率图像初分割缝洞子图像上既包含裂缝、溶蚀孔洞，也包含水平状高电导纹层；以目标形状参数为基础的二次图像分割子图像即裂缝、孔洞子图像则去掉了水平状高电导纹层，仅包含裂缝、溶蚀孔洞。图5中，代号p1az表示1号极板方位，deg表示度，为角度单位，heated表示图头，fmi_256表示沿井眼正北方向剖开的电导率图像，fmi_fill_liu3_032_256表示电成像测井全井眼覆盖图像，md代表测量深度，fmi_fill_liu3_032_256_seg_fmi_v表示二维二进小波变换初分割缝洞子图像，即全井眼电导率图像初分割缝洞子图像，fmi_fill_liu3_032_256_seg_fmi_v_sd表示以目标形状参数为基础的二次分割子图像，即裂缝、孔洞子图像。
[0085]
步骤105，基于裂缝、孔洞子图像，提取电成像缝洞参数。
[0086]
在一种可能的实现方式中，基于裂缝、孔洞子图像，提取电成像缝洞参数包括：对于二次分割出来的裂缝、孔洞子图像，在单目标参数计算的基础上，对一个图像框所在区域提取电成像缝洞参数；移动图像框，对其他区域提取电成像缝洞参数。
[0087]
可选地，电成像缝洞参数包括但不限于：图像框内多目标的平均圆度、图像框内多目标的平均长度、图像框内多目标的平均宽度、图像框内孔洞-裂缝面孔率、图像框内孔洞面孔率及图像框内裂缝面孔率参数。进一步可选地，上述电成像缝洞参数用于表征缝洞储层的好坏。
[0088]
本技术实施例提供的技术方案至少包括以下有益效果：1、以完整的二维方式(全井眼电导率数据)进行图像分割，并以目标形状参数为基础进行二次分割，分割出来的缝洞子图像更为清晰完整。2、分割出来的高电导目标剔除了水平状高电导泥质纹层、缝合线及诱导缝等无效信息，使电成像缝洞参数提取结果更为准确。
[0089]
请参考图6，其示出了本发明的一个实施例提供的一种电成像缝洞参数提取方法所获取的电成像测井全井眼覆盖图像、全井眼电导率图像初分割缝洞子图像、裂缝、孔洞子图像以及电成像缝洞参数的提取结果。该实施例的电成像测井仪为fmi测井仪，参考极板为1号极板。
[0090]
其中，图6第一道为全井眼覆盖的电导率图像，即电成像测井全井眼覆盖图像；第二、四、六为深度道，表示测井深度；第三道为二维二进小波变换图像初分割结果图，即全井眼电导率图像初分割缝洞子图像；第七、八道为电成像缝洞参数的提取结果。进一步地，第七、八道中的代号aphit为缝洞面孔率，代号vvphit为溶蚀孔洞面孔率，代号ffphit为裂缝面孔率，代号xphit为缝洞管流模型产量指数。代号p1az表示1号极板方位，deg表示度，为角度单位，heated表示图头，fmi_fill_liu3_032_256表示电成像测井全井眼覆盖图像，md表
示测量深度，fmi_fill_liu3_032_256_seg_fmi_v表示二维二进小波变换初分割缝洞子图像，即全井眼电导率图像初分割缝洞子图像，fmi_fill_liu3_032_256_seg_fmi_v_sd表示以目标形状参数为基础的二次分割子图像，即裂缝、孔洞子图像。
[0091]
请参考图7，其示出了本发明的一个实施例提供的一种电成像缝洞参数提取方法所获取的电成像缝洞参数的提取结果与常规测井资料对比结果。该实施例的电成像测井仪为fmi测井仪，参考极板为1号极板。
[0092]
其中，图7第一道为岩性曲线段，其中gr(gamma ray)为自然伽马曲线，kth(kalium thorium)为去铀自然伽玛曲线，cal(caliper)为井径曲线；第二道为深度道；第三道为电阻率曲线道，其中代号rt为深侧向电阻率测井曲线，rxo为浅侧向电阻率测井曲线；第四道为孔隙度曲线道，其中ac(acoustictime)为声波时差曲线，den(density)为密度曲线，cnl(compensated neutron log)为中子孔隙度曲线；第五道为计算的孔隙度曲线；第六道为深度道；第七道为电成像测井全井眼覆盖图像；第八道为全井眼电导率图像初分割缝洞子图像；第九道为裂缝、孔洞子图像；第十道为电成像缝洞参数提取结果图，其中aphit为缝洞面孔率，vvphit为溶蚀孔洞面孔率，ffphit为裂缝面孔率。代号p1az表示1号极板方位，deg表示度，为角度单位，heated表示图头，fmi_fill_liu3_032_256表示电成像测井全井眼覆盖图像，fmi_fill_liu3_032_256_seg_fmi_v表示二维二进小波变换初分割缝洞子图像，即全井眼电导率图像初分割缝洞子图像，fmi_fill_liu3_032_256_seg_fmi_v_sd表示以目标形状参数为基础的二次分割子图像，即裂缝、孔洞子图像，reference表示深度，gapi表示自然伽马测井值的单位，in表示英寸，ohmm为欧姆.米，表示电阻率单位，us/ft为微妙每英尺，表示声波时差单位，v/v为体积/体积，表示中子孔隙度测井单位，。
[0093]
从电成像测井资料上来看，7064.4至7069.5米井段溶蚀裂缝和溶蚀孔洞发育，渗透性好，钻井过程中泥浆侵入溶蚀裂缝和溶蚀孔洞，测量的双侧向曲线值降低。另外，溶蚀孔洞和溶蚀裂缝发育在孔隙度测井曲线上也有响应。溶蚀裂缝和溶蚀孔洞发育处对应密度曲线(代号为：den)值减小、声波曲线(代号为：ac)值增大、中子孔隙度(代号为：cnl)曲线值略有增大，对应计算的孔隙度值(代号为：pige_qepp)增大。本发明实施例所获取的电成像缝洞参数的提取结果与常规常规双侧向电阻率和孔隙度曲线是对应的。7064.4至7069.5米井段二次分割缝洞子图像上见分割出来的裂缝和溶蚀孔洞目标，提取的的裂缝、孔洞面孔隙度大，与双侧向曲线有很好的反相关关系，与孔隙度曲线有很好的正相关关系。
[0094]
请参考图8，其示出了本发明的一个实施例提供的一种电成像缝洞参数提取方法所获取的电成像缝洞参数的提取结果与相关技术所获取的缝洞参数的提取结果的对比。该实施例的电成像测井仪为fmi测井仪，参考极板为1号极板。
[0095]
其中，图8第二道为电成像电导率静态图像，第四道为相关技术处理得到的电成像缝洞参数提取结果，第六道为本发明实施例得到的电成像缝洞参数提取结果。与相关技术相比，本发明实施例分割出来的缝洞子图像更为清晰完整，并且不包含水平状高电导纹层等非有效目标。图8中，array_wbi_scale2表示浅电阻率刻度后的电导率图像，即电成像电导率静态图像，array_abi_seg表示相关技术处理得到的电成像缝洞参数提取结果，image orientation表示图像方位，代号p1az表示1号极板方位，deg表示度，为角度单位，heated表示图头，fmi_fill_liu3_032_256_seg_fmi_v_sd表示以目标形状参数为基础的二次分割子图像，即裂缝、孔洞子图像，reference表示深度。
[0096]
请参考图9，其示出了本发明的一个实施例提供的一种电成像缝洞参数提取方法在斯伦贝谢techlog环境中实现的流程图，该方法包括如下几个步骤：
[0097]
步骤901，在techlog软件中加载资料。
[0098]
可选地，测井资料包括：常规资料以及电成像测井资料。
[0099]
步骤902，基于所加载的资料，用techlog软件进行预处理。
[0100]
可选地，通过techlog软件中的电成像测井资料预处理工作流对电成像测井资料进行预处理。
[0101]
步骤903，基于预处理后的资料，用techlog软件获取电成像测井全井眼覆盖图像。
[0102]
基于预处理后的资料，用techlog软件获取电成像测井全井眼覆盖图像，包括但不限于如下几个子步骤：
[0103]
9031、基于预处理后的资料，用techlog软件获取电成像测井全井眼电导率数据。
[0104]
基于预处理后的资料，用techlog软件获取电成像测井全井眼电导率数据，包括但不限于如下子步骤：
[0105]
9031
‑①
从预处理后的资料中，用techlog软件获取原始电导率数据。
[0106]
其中，原始电导率数据为经浅电阻率刻度后的极板电导率数据。由于电成像测井仪的纽扣电极系是非聚集电极系，测量值仅与井壁附近地质体的电导率成比例变化，因此要用浅电阻率进行刻度，得到经浅电阻率刻度后的极板电导率数据。
[0107]
9031
‑②
基于原始电导率数据，用techlog软件绘制原始电成像电导率图像。
[0108]
基于原始电导率数据，根据参考极板中间纽扣电极电导率曲线的方位，沿井筒正北剖开，用techlog软件绘制电成像电导率图像。
[0109]
可选地，不同仪器不同参考极板中间纽扣电极不同，如xrmi电成像测井仪将1号极板作为参考极板，xrmi中1号极板第13个纽扣电极为中间纽扣电极；fmi电成像测井仪将1号极板作为参考极板，fmi中1号极板第12个纽扣电极为中间纽扣电极。
[0110]
9031
‑③
基于原始电成像电导率图像，应用电成像极板间空白充填的方法，在techlog软件中获取电成像测井全井眼电导率数据。
[0111]
可选地，电成像极板间空白充填的方法指的是基于原始电成像电导率图像，应用极板边缘向空隙中间生长的方法，生成电成像测井仪极板间空白带数据，从而得到电成像测井全井眼覆盖的电导率数据，即电成像测井全井眼电导率数据。
[0112]
9032、基于电成像测井全井眼电导率数据，用techlog软件获取电成像测井全井眼覆盖图像。
[0113]
根据电成像测井全井眼电导率数据，利用1号极板方位数据计算每个极板相对正北方向的方位，在techlog软件中绘制得到极板间空白带充填后的全井眼覆盖的电导率图像，即电成像测井全井眼覆盖图像。
[0114]
步骤904，用techlog软件对电成像测井全井眼覆盖图像进行初分割，获取全井眼电导率图像初分割缝洞子图像。
[0115]
用techlog软件对电成像测井全井眼覆盖图像进行初分割，获取全井眼电导率图像初分割缝洞子图像，包括但不限于如下几个子步骤：
[0116]
9041、基于电成像测井全井眼覆盖图像，用techlog软件获取一个图像框所在区域的全井眼电导率图像初分割缝洞子图像。
[0117]
可选地，基于电成像测井全井眼覆盖图像，用techlog软件获取一个图像框所在区域的全井眼电导率图像初分割缝洞子图像包括但不限于如下子步骤：
[0118]
9041
‑①
在techlog软件中输入一个图像框的电成像测井全井眼覆盖图像数据；
[0119]
9041
‑②
用techlog软件对一个图像框的电成像测井全井眼覆盖图像数据进行二维二进小波变换，对于阶数为2的小波变换计算每个电导率数据点小波变换系数的模值，得到2阶小波变换系数模值。
[0120]
在一个示例中，将一个图像框的电成像测井全井眼覆盖图像中的电导率数据设为:
[0121][0122]
进行二次小波变换，得到该电导率数据对应的小波系数:
[0123][0124]
其中，j为小波变换阶数，取电导率数据2阶小波变换的小波系数：
[0125][0126]
计算2阶小波变换系数的模值：
[0127][0128]
9041
‑③
基于2阶小波变换系数模值，用techlog软件在深度方向上利用三点寻峰获取模极大值点。
[0129]
9041
‑④
基于模极大值点，用techlog软件获取模极大值点的坐标像素值，并根据模极大值点的坐标像素值确定分割阈值进行图像初分割，获取一个图像框所在区域的全井眼电导率图像分割缝洞子图像。
[0130]
模极大值点，即全井眼电导率图像电导率值突变点，进一步地，即分割点。可选地，根据分割点位置确定分割阈值，进一步可选地，根据分割阈值进行图像初分割。
[0131]
9042、移动图像框，用techlog软件获取其他区域的全井眼电导率图像初分割缝洞子图像。
[0132]
可选地，移动图像框，用techlog软件获取其他区域的全井眼电导率图像初分割缝洞子图像包括：移动图像框，在techlog软件中重复上述子步骤9041的方法，直至测井终点深度，得到其他区域的全井眼电导率图像初分割缝洞子图像。其中，测井终点深度从图像中获得。
[0133]
9043、基于一个图像框所在区域的全井眼电导率图像初分割缝洞子图像以及其他区域的全井眼电导率图像初分割缝洞子图像，在techlog软件中输出全井眼电导率图像初分割缝洞子图像。
[0134]
步骤905，用techlog软件提取全井眼电导率图像初分割缝洞子图像的目标形状参数。
[0135]
可选地，用techlog软件提取全井眼电导率图像初分割缝洞子图像的目标形状参数包括：基于全井眼电导率图像初分割缝洞子图像，用techlog软件提取高电导目标的角度参数，将高电导目标的角度参数作为全井眼电导率图像初分割缝洞子图像的目标形状参数。
[0136]
在剖开的全井眼覆盖的二维电导率图像上，水平状高电导泥质纹层、缝合线等非
有效高电导目标与裂缝、孔洞目标的角度不同。为了应用初步分割目标形状参数区分非有效高电导目标与裂缝、孔洞目标，对初分割出来的缝洞子图像进行进一步处理，提取高电导目标的角度参数。
[0137]
进一步可选地，高电导目标的角度参数定义为高电导目标长轴(目标边缘上距离最大的两点之间的连线)与水平方向的夹角。进一步可选地，按照如下公式定义角度参数：
[0138][0139]
其中，水平方向为x轴，高电导目标长轴为y轴，(x2，y2)为高电导目标长轴的终点，(x1，y1)为高电导目标长轴的起点。
[0140]
步骤906，基于全井眼电导率图像初分割缝洞子图像的目标形状参数，用techlog软件对全井眼电导率图像初分割缝洞子图像进行二次分割，获取裂缝、孔洞子图像。
[0141]
可选地，基于全井眼电导率图像初分割缝洞子图像的目标形状参数，用techlog软件对全井眼电导率图像初分割缝洞子图像进行二次分割，获取裂缝、孔洞子图像，包括：在全井眼电导率图像初分割缝洞子图像上，用techlog软件去掉目标形状参数小于有效阈值的非有效目标，得到裂缝、孔洞子图像。
[0142]
在直井中，例如水平状或倾角较低的高电导泥质纹层、缝合线等非有效目标的角度参数很小。在一个示例中，以10
°
作为目标形状参数的有效阈值。进一步地，将目标形状参数小于10
°
，即角度参数小于10
°
的高电导目标作为非有效目标，在techlog软件中对全井眼电导率图像初分割缝洞子图像进行二次分割，去掉非有效目标，得到裂缝、孔洞子图像。其中，有效阈值可基于经验设置，本技术实施例对此不进行限定。
[0143]
步骤907，基于裂缝、孔洞子图像，用techlog软件提取电成像缝洞参数。
[0144]
可选地，基于裂缝、孔洞子图像，用techlog软件提取电成像缝洞参数包括：对于二次分割出来的裂缝、孔洞子图像，在单目标参数计算的基础上，用techlog软件对一个图像框所在区域提取电成像缝洞参数，其中，一个图相框相当于常规测井资料的采样间隔；移动图像框，用techlog软件对其他区域提取电成像缝洞参数。
[0145]
可选地，电成像缝洞参数包括：图像框内多目标的平均圆度、图像框内多目标的平均长度、图像框内多目标的平均宽度、图像框内孔洞-裂缝面孔率、图像框内孔洞面孔率及图像框内裂缝面孔率参数。进一步可选地，上述电成像缝洞参数用于表征缝洞储层的好坏。
[0146]
步骤908，用techlog软件对所提取的电成像缝洞参数，进行绘图。
[0147]
请参考图10，其示出了一种电成像缝洞参数的提取系统100，该系统包括如下几个模块：全井眼覆盖图像生成模块1001、初分割模块1002、目标形状参数提取模块1003、二次分割模块1004以及缝洞参数提取模块1005。
[0148]
全井眼覆盖图像生成模块1001，用于获取电成像测井全井眼覆盖图像。可选地，全井眼覆盖图像生成模块1001用于但不限于：
[0149]
10011、获取电成像测井全井眼电导率数据。
[0150]
进一步可选地，全井眼覆盖图像生成模块1001用于但不限于：
[0151]
10011
‑①
输入原始电导率数据。
[0152]
在一种可能实现的方式中，输入的原始电导率数据为经浅电阻率刻度后的极板电导率数据。由于电成像测井仪的纽扣电极系是非聚集电极系，测量值仅与井壁附近地质体
的电导率成比例变化，因此要用浅电阻率测井资料进行刻度，得到经浅电阻率刻度后的极板电导率数据。
[0153]
10011
‑②
基于原始电导率数据，绘制原始电成像电导率图像。
[0154]
基于原始电导率数据，根据参考极板中间纽扣电极电导率曲线的方位，沿井筒正北剖开，绘制电成像电导率图像。
[0155]
可选地，不同电成像测井仪不同参考极板中间纽扣电极不同，如xrmi电成像测井仪将1号极板作为参考极板，xrmi中1号极板第13个纽扣电极为中间纽扣电极，fmi电成像测井仪将1号极板作为参考极板，fmi中1号极板第12个纽扣电极为中间纽扣电极。
[0156]
10011
‑③
基于原始电成像电导率图像，应用电成像极板间空白充填的方法，获取电成像测井全井眼电导率数据。
[0157]
可选地，电成像极板间空白充填的方法指的是基于原始电成像电导率图像，应用极板边缘向空隙中间生长的方法，生成电成像测井仪极板间空白带数据，从而得到电成像测井全井眼覆盖的电导率数据，即电成像测井全井眼电导率数据。
[0158]
10012、基于电成像测井全井眼电导率数据，获取电成像测井全井眼覆盖图像。
[0159]
在一种可能实现的方式中，根据电成像测井全井眼电导率数据，利用参考极板方位数据计算每个极板相对正北方向的方位，绘制得到极板间空白带充填后的全井眼覆盖的电导率图像，即电成像测井全井眼覆盖图像。
[0160]
初分割模块1002，用于对电成像测井全井眼覆盖图像进行初分割，获取全井眼电导率图像初分割缝洞子图像。可选地，初分割模块1002用于但不限于：
[0161]
10021、基于电成像测井全井眼覆盖图像，获取一个图像框所在区域的全井眼电导率图像初分割缝洞子图像。
[0162]
进一步可选地，初分割模块1002用于但不限于：
[0163]
10021
‑①
输入一个图像框的电成像测井全井眼覆盖图像数据；
[0164]
10021
‑②
对一个图像框的电成像测井全井眼覆盖图像数据进行二维二进小波变换，对于阶数为2的小波变换计算每个电导率数据点小波变换系数的模值，得到2阶小波变换系数模值。
[0165]
在一个示例中，将一个图像框的电成像测井全井眼覆盖图像中的电导率数据设为:
[0166][0167]
进行二次小波变换，得到该电导率数据对应的小波系数:
[0168][0169]
其中，j为小波变换阶数，取电导率数据2阶小波变换的小波系数：
[0170][0171]
计算2阶小波变换系数的模值：
[0172][0173]
10021
‑③
基于2阶小波变换系数模值，在深度方向上利用三点寻峰获取模极大值点。
[0174]
10021
‑④
基于模极大值点，获取模极大值点的坐标像素值，根据模极大值点的坐
标像素值确定分割阈值进行图像初分割，获取一个图像框所在区域的全井眼电导率图像分割缝洞子图像。
[0175]
模极大值点，即全井眼电导率图像电导率值突变点，进一步地，即分割点。可选地，根据分割点位置确定分割阈值，进一步可选地，根据分割阈值进行图像初分割。
[0176]
10022、移动图像框，获取其他区域的全井眼电导率图像初分割缝洞子图像。
[0177]
在一种可能实现的方式中，初分割模块1002用于：移动图像框，重复上述子步骤10021的方法，直至测井终点深度，得到其他区域的全井眼电导率图像初分割缝洞子图像。其中，测井终点深度从图像中获得。
[0178]
10023、基于一个图像框所在区域的全井眼电导率图像初分割缝洞子图像以及其他区域的全井眼电导率图像初分割缝洞子图像，输出全井眼电导率图像初分割缝洞子图像。
[0179]
目标形状参数提取模块1003，用于提取全井眼电导率图像初分割缝洞子图像的目标形状参数。
[0180]
可选地，目标形状参数提取模块1003用于：基于全井眼电导率图像初分割缝洞子图像，提取高电导目标的角度参数。在剖开的全井眼覆盖的二维电导率图像上，水平状高电导泥质纹层、缝合线等非有效高电导目标与裂缝、孔洞目标的角度不同。为了应用初步分割目标形状参数区分非有效高电导目标与裂缝、孔洞目标，对初分割出来的缝洞子图像进行进一步处理，提取高电导目标的角度参数。
[0181]
进一步可选地，高电导目标的角度参数定义为高电导目标长轴(目标边缘上距离最大的两点之间的连线)与水平方向的夹角。
[0182]
二次分割模块1004，用于基于全井眼电导率图像初分割缝洞子图像的目标形状参数，对全井眼电导率图像初分割缝洞子图像进行二次分割，获取裂缝、孔洞子图像。
[0183]
可选地，二次分割模块1004用于：在全井眼电导率图像初分割缝洞子图像上，去掉目标形状参数小于有效阈值的非有效目标，得到裂缝、孔洞子图像。在直井中，例如水平状或倾角较低的高电导泥质纹层、缝合线等非有效目标的角度参数很小。在一个示例中，以10
°
作为目标形状参数的有效阈值。进一步可选地，二次分割模块1004用于将目标形状参数小于10
°
，即角度参数小于10
°
的高电导目标作为非有效目标，对全井眼电导率图像初分割缝洞子图像进行二次分割，去掉非有效目标，得到裂缝、孔洞子图像。
[0184]
缝洞参数提取模块1005，用于基于裂缝、孔洞子图像，提取电成像缝洞参数。可选地，缝洞参数提取模块1005用于：对于二次分割出来的裂缝、孔洞子图像，在单目标参数计算的基础上，对一个图像框所在区域提取电成像缝洞参数，其中，一个图相框相当于常规测井资料的采样间隔；移动图像框，对其他区域提取电成像缝洞参数。
[0185]
可选地，电成像缝洞参数包括：图像框内多目标的平均圆度、图像框内多目标的平均长度、图像框内多目标的平均宽度、图像框内孔洞-裂缝面孔率、图像框内孔洞面孔率及图像框内裂缝面孔率参数。进一步可选地，上述电成像缝洞参数用于表征缝洞储层的好坏。
[0186]
应理解的是，上述图10提供的系统在实现其功能时，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将设备的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。另外，上述实施例提供的系统与方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不
再赘述。
[0187]
在示例性实施例中，还提供了一种计算机可读存储介质，该存储介质中存储有至少一条程序代码，该至少一条程序代码由处理器加载并执行，以使计算机实现上述任一种电成像缝洞参数的提取方法。
[0188]
可选地，上述计算机可读存储介质可以是只读存储器(read-only memory，rom)、随机存取存储器(random access memory，ram)、只读光盘(compact disc read-only memory，cd-rom)、磁带、软盘和光数据存储设备等。
[0189]
在示例性实施例中，还提供了一种计算机程序或计算机程序产品，该计算机程序或计算机程序产品中存储有至少一条计算机指令，该至少一条计算机指令由处理器加载并执行，以使计算机实现上述任一种电成像缝洞参数的提取方法。
[0190]
在本技术所提供的几个实施例中，应该理解到，所揭露的系统和方法，可以通过其它的方式实现。例如，以上所描述的系统实施例仅仅是示意性的，例如，该模块的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个模块或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另外，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口、设备或模块的间接耦合或通信连接，也可以是电的，机械的或其它的形式连接。
[0191]
该作为分离部件说明的模块可以是或者也可以不是物理上分开的，作为模块显示的部件可以是或者也可以不是物理模块，即可以位于一个地方，或者也可以分布到多个网络模块上。可以根据实际的需要选择其中的部分或者全部模块来实现本技术实施例方案的目的。
[0192]
另外，在本技术各个实施例中的各功能模块可以集成在一个处理模块中，也可以是各个模块单独物理存在，也可以是两个或两个以上模块集成在一个模块中。上述集成的模块既可以采用硬件的形式实现，也可以采用软件功能模块的形式实现。
[0193]
还应理解，在本技术的各个实施例中，各个过程的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本技术实施例的实施过程构成任何限定。
[0194]
本技术中术语“至少一个”的含义是指一个或多个，本技术中术语“多个”的含义是指两个或两个以上，例如，多个数据是指两个或两个以上的数据。
[0195]
应理解，在本文中对各种所述示例的描述中所使用的术语只是为了描述特定示例，而并非旨在进行限制。如在对各种所述示例的描述和所附权利要求书中所使用的那样，单数形式“一个(“a”，“an”)”和“该”旨在也包括复数形式，除非上下文另外明确地指示。
[0196]
以上仅为本技术的示例性实施例，并不用以限制本技术，凡在本技术的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本技术的保护范围之内。