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

说明书

技术领域

本发明属于地质图像处理技术领域，具体涉及一种生物礁储层的几何形态识别方法。

背景技术

在漫长的地质历史中，从太古代的微生物藻礁到现代的珊瑚和珊瑚藻礁，可以说生物礁是地球上最早的生命作用的产物。生物礁常形成隆凸的构造，再加上原生和成岩后生作用产生的孔隙度和渗透率，有利于石油天然气聚集，并且在被其它沉积盖层覆盖之后常形成良好的储油构造，因此生物礁的研究具有非常重要的意义。

近年来在地震勘探方面针对生物礁储层的研究方法层出不穷，其研究主要集中在如下方面：（1）储层地震响应特征的数值模拟；（2）古地貌分析，如地层厚度、层拉平等；（3）基于地震波形分类的地震相分析；（4）地震属性分析；（5）叠后地震反演分析；（6）时频分析方法；（7）地震波吸收衰减分析；（8）基于Gassman方程、岩石物理模型的流体替换分析；（9）叠前弹性参数反演；（10）流体识别因子分析等。但生物礁储层的钻井成功率依然很低，这是因为其地球物理预测难度十分大：

（1）由于在形成的过程中会受到地质、沉积、气候等诸多因素的影响，能够成形生物礁储层的空间区域很苛刻；

（2）生物礁储层空间分布很不规则；

（3）生物礁储层埋藏的深度大，在获得的地震资料上的分辨率和信噪比较低；

（4）生物礁的类型多，内部结构复杂，地球物理响应特征变化大。

因此，现有的侧重于地层地质方面的研究并不能很好的识别生物礁储层。

发明内容

本发明的目的在于克服现有技术中的上述问题，提供一种基于纹理属性、傅里叶频谱特征分析及混沌纹理属性等原理的一种识别生物礁储层几何形态的方法。

为解决上述技术问题，本发明采用以下技术方案：

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

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

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

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

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

进一步地，S2中采用基于灰度共生矩阵导出的纹理属性特征提取法进行特征提取。

进一步地，所述基于灰度共生矩阵导出的纹理属性的特征提取法中作为提取特征的为灰度级为16、元素为两个间距d为1的像素的联合概率密度的灰度共生矩阵的0°、45°、90°、135°四个方向导出的四组特征参数的均值。

进一步地，所述特征参数包括：角二阶矩（ASM）、对比度、相关、熵、逆差距和均一性度量。

进一步地，所述S1中对原始地震剖面数据进行预处理前先提取原始地震剖面数据的混沌纹理属性。

进一步地，S2中采用基于傅里叶频谱特征的特征提取法进行特征提取。

进一步地，所述基于傅里叶频谱特征的特征提取法包括以下步骤：

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