[发明专利]基于GST的碳酸盐岩非均质储层中溶洞识别标定方法

权利要求书

1.一种基于GST的碳酸盐岩非均质储层中溶洞识别标定方法，其特征在于，包括以下步骤：

(1)根据三维叠后地震数据进行结构梯度张量的构建，并计算其特征值，所述特征值为λ1、λ2、λ3，且λ1≥λ2≥λ3；

(2)根据λ1、λ2，构建溶洞识别算子值根据C_cave值，构建C_cave属性数据体；

(3)从C_cave属性数据体中抽取数个溶洞识别属性剖面，将其处理为二值化溶洞剖面，溶洞区域表现为1值，其他区域表现为0值；

(4)用二值图像连通区域标记法对二值化溶洞剖面进行连通域标记，得到连通域标记后的剖面，标记后的连通域为C，该二值化溶洞剖面为目标空间D，所述二值图像连通区域标记法为一步法；

(5)从连通域标记后的剖面中进行溶洞定位及几何参数的提取；

(51)连通域质心确定，质心(x_m，y_m)

式中，x_m为连通域质心的横坐标，y_m为质心的纵坐标；p_i为连通域C内点的像素值，在二值图像中，图像的像素值为1，即p_i＝1；

(52)连通域宽度确定：连通域宽度指在一个连通域内其最左边的像素点到最右边的像素点的X坐标之差的的绝对值；

w_c＝|x_max-x_min|

式中，x_min为在连通域C中最小横坐标值，即最左边的横坐标值，x_max为在连通域C最大横坐标值，即最右边的横坐标值；w_c为连通域C的宽度；

(53)连通域高度确定：连通域高度指最顶端的纵坐标到最低端的纵坐标的差的绝对值；

h_c＝|y_max-y_min|

式中，y_min为在连通域C中最小横坐标值，即最顶端的纵坐标值，y_max为在连通域C最底端纵坐标值，即最右边的横坐标值；

(54)连通域的外切圆半径R_外及内切圆半径R_内：连通域的最左、最右和最顶、最低的四组坐标，构成四个顶点绘制一个矩形，R_内为该矩形短边长的1/2，R_外为该矩形长对角线的1/2。

2.根据权利要求1所述的基于GST的碳酸盐岩非均质储层中溶洞识别标定方法，其特征在于：还包括以下步骤：

(6)定义溶洞的发育程度评估参数

(61)扁平度e：

当连通域形状接近圆形，其标记矩形为正方形，e＝1，当连通域形状越扁平，扁平度e越大；

(62)溶洞体面孔率CVPA：

在连通域标记后的剖面上开设窗口，窗口长度为窗长，A_i为窗长内第i个孔洞的面积；A_w为窗长面积；N为窗长内的孔洞个数，CVPA为面孔率，根据下式计算面孔率

(63)溶洞的空间聚集度I_c：

式中，A_i为窗长内第i个孔洞的面积，x_i为连通域的质心的x坐标，y_i为各连通域在质心的y坐标，为该区域的连通域的面积加权中心的横坐标、纵坐标，σ_x、σ_y为各连通域质心位置的方差，I_c为连通域的空间聚集度；

(64)所选窗口内溶洞的分选度S_C：指储层中溶洞尺度的均匀程度、或溶洞尺度围绕集中趋势的偏离程度；

其中，以目标空间D中，所有连通域内切圆半径或者外接圆半径做粒度累计频率曲线，r_25％和r_75％为粒度累计频率曲线中，累计频率在25％和75％所对应的溶洞内切圆半径或者外接圆半径。