[发明专利]基于方位角的灰度共生矩阵预测碳酸盐岩裂缝储层的方法

权利要求书

1.一种基于方位角的灰度共生矩阵预测碳酸盐岩裂缝储层的方法，其特征在于：包括以下步骤：

(1)读取一时长为T的3D地震数据体；

(2)求取3D地震数据体中所有像素点的方位角，得到方位角体，并将3D地震数据体灰度化，得到3D灰色图像体；

(3)抽取t＝0时，方位角体和3D灰色图像体的水平切片，t为采样时间；

(4)将方位角体的水平切片标记为Qt，再将3D灰色图像体V的水平切片压缩到N阶，得到N阶灰度图像，标记为Vt，且N阶灰度图像中每个像素点对应一灰度值；

(5)计算Vt上所有像素点带方位角的灰度共生矩阵Lt，并计算Lt对应的特征参数Pt；

(51)在Vt上设定一窗口W，其高度和宽度为N的倍数，其中心像素点对应在Qt中的方位角为Az；

(52)构建窗口W的灰度共生矩阵Lw，Lw为一N×N的二维数组，其元素为Lw(i，j)，i＝1，2，…，N，j＝1，2，…，N；

对于每个元素Lw(i，j)，在W窗口中，统计方向为Az、距离为d的两个像素值为(i，j)的像素对个数，并将该个数作为该元素的元素值；

(53)以Lw为基础，计算W中心像素点对应的特征参数Pw，所述Pw包括对比度CON、能量ASM、熵ENT、相关值COR；

(54)在Vt上移动窗口W，重复步骤(52)(53)，得到Vt中所有像素点的特征参数Pw、按行列排列得到t＝0时，Vt特征参数Pt；

(6)依次抽取t＝1，2，3…，T时，方位角体Q和3D灰色图像体V的水平切片，每抽取一次，重复步骤(4)(5)，最终得到t＝1，2，3…，T时，Vt的特征参数Pt；

(7)将所有特征参数构成3D特征参数属性体P_D。

2.根据权利要求1所述的基于方位角的灰度共生矩阵预测碳酸盐岩裂缝储层的方法，其特征在于：所述步骤(2)中，方位角Az采用下式得到：

其中，v₁(x)、v₁(y)是通过梯度结构张量方法计算出x,y方向上的特征向量。

3.根据权利要求1所述的基于方位角的灰度共生矩阵预测碳酸盐岩裂缝储层的方法，其特征在于：还包括以下步骤：

(8)选择一待测时间层作为目标层，从P_D中抽取目标层的特征参数；

(9)将特征参数中，对比度CON、能量ASM、熵ENT、相关值COR四种参数均一化后等权加融合得到综合参数P，利用综合参数P对碳酸盐岩裂缝储层发育特征进行预测。

4.根据权利要求1所述的基于方位角的灰度共生矩阵预测碳酸盐岩裂缝储层的方法，其特征在于：所述N＝8或16。