[发明专利]生物礁储层的几何形态识别方法

权利要求书

1.一种生物礁储层的几何形态识别方法，其特征在于：包括以下步骤：

S1：对原始地震剖面数据进行预处理，包括根据地震剖面大小选取滑动窗口和滑动窗口的移动步长；

S2：对S1中预处理后的地震剖面数据进行特征提取；

S3：对S2所提取的特征进行生物礁储层和非生物礁储层的聚类；

S4：对聚类后的生物礁储层和非生物礁储层贴以不同颜色的聚类标签以显示在地震剖面上。

2.根据权利要求1所述的生物礁储层的几何形态识别方法，其特征在于：S2中采用基于灰度共生矩阵导出的纹理属性特征提取法进行特征提取。

3.根据权利要求2所述的生物礁储层的几何形态识别方法，其特征在于：所述基于灰度共生矩阵导出的纹理属性的特征提取法中作为提取特征的为灰度级为16、元素为两个间距d为1的像素的联合概率密度的灰度共生矩阵的0、°45、°90、°135°四个方向导出的四组特征参数的均值。

4.根据权利要求2或3所述的生物礁储层的几何形态识别方法，其特征在于：所述特征参数包括：角二阶矩（ASM）、对比度、相关、熵、逆差距和均一性度量。

5.根据权利要求2所述的生物礁储层的几何形态识别方法，其特征在于：所述S1中对原始地震剖面数据进行预处理前先提取原始地震剖面数据的混沌纹理属性。

6.根据权利要求1所述的生物礁储层的几何形态识别方法，其特征在于：S2中采用基于傅里叶频谱特征的特征提取法进行特征提取。

7.根据权利要求6所述的生物礁储层的几何形态识别方法，其特征在于：所述基于傅里叶频谱特征的特征提取法包括以下步骤：

S21：将S1中预处理后的地震剖面图像的傅里叶功率谱分成若干个等间距同心矩形环，将若干个标记为M个；设图像f(x,y)的大小为L×W，将图像置于笛卡尔坐标系中，使图像中心的坐标为（L/2，W/2），其中，L表示图像的长度，W表示图像的宽度，则图像的离散傅里叶变换为：F(u,v)=Σx=0L-1Σy=0W-1f(x,y)exp[-j2π(uxL+vyW)]=R(u,v)+jI(u,v)---(1)]]>

其中，u＝0,...,L-1,v＝0,...,W-1；x，y是图像的横、纵坐标轴；F(u,v)是频域图像谱，通常F(u,v)是两个实频率变量u和v的复值函数，频率u是对应于x轴，频率v是对应于y轴；R(u,v)是F(u,v)的实部，I(u,v)是F(u,v)的虚部；

点(u,v)的能量谱为：

p(u,v)＝R²(u,v)+I²(u,v)     （2）

S22：求出每个矩形环功率谱的能量；

每个矩形环功率谱的能量为：

E_i＝∑p(u,v)        （3）

其求和范围为：L/2M*i≤|u-L/2|＜L/2M*(i+1)，

W/2M*i≤|v-W/2|＜W/2M*(i+1),i＝0,1,2,3M...-1；

S23：求每个矩形环功率谱能量占所有矩形环功率谱总能量的百分比；

百分比的计算公式如式（4）所示：

fhi=Ei/Σu=0LΣv=0Wp(u,v)---(4)]]>

将所求得百分比fh_i作为提取特征进行特征提取。