地震纹理属性应用研究.docVIP

地震纹理属性在JJD工区断层识别中的应用 张军华1 王伟2 谭明友3 崔世凌3 陈海云3 （1, 中国石油大学（华东）资源信息学院, 东营, 257061；2，东方地球物理公司研究院, 涿州, 072751； 3, 中国石化胜利油田物探研究院, 东营, 257022） 摘要：纹理分析是图像处理中的一种常用技术，通过构建地震纹理基元和灰度共生矩阵Bahorich和Farmer【6】的第一代的互相关算法、发展到Marfurt【7,8】的第二代多道相似算法和第三代特征值相干算法，近几年以Gao Dengliang【9-11】为首的研究小组利用纹理属性进行断层识别研究并取得了一系列成果。本文借鉴其纹理分析的相关方法原理，研究了三维纹理体的提取算法，获得了纹理能量、熵、对比度和相关性等纹理属性体，讨论了纹理基元大小选择、灰度级确定、RGB数据融合等应用要素，并在JJD实际工区中取得了较好的应用效果。 2．基于灰度共生矩阵的纹理分析方法 灰度共生矩阵是一种用来分析图像纹理特征的重要方法，由Haralick 于1973提出通过计算图像中一定距离和一定方向的两个像素之间的灰度相关性，对图像的所有像素进行统计，反映图像在方向、相邻间隔、变化幅度及快慢上的综合信息【】。对于一个图像来说，统计距离为δ的像素出现的概率(i, j)： （1） 式中i, j=0, 1, 2, …, L- 1，L为图像的灰度级数x，y是图像中的像素坐标 是处的灰度值；；Nx, Ny 分别为图像的行列数。灰度共生矩阵 （2） 容易推知，R是一对称阵，位置偏移由x方向和y方向的两个量组成。一般来说，较小则反应图像的整体纹理分布，而较大的则反应小区域的细微变化。 用灰度共生矩阵可提取多种特征值，最常用的是以下四个特征： （1）角二阶矩或能量 它是纹理灰度变化均一的度量，反映了灰度分布均匀程度和纹理粗细度。 （3） （2）熵 熵可以度量纹理的随机性。当共生矩阵中所有值均相等时，它取得最大值；相反，如果共生矩阵中的值非常不均匀时，其值较小。 （4） （3）对比度 对比度是灰度共生矩阵主对角线附近的惯性矩，它度量矩阵值的分布和局部变化，可反映了清晰度和纹理的沟纹深浅。 （5） （4）相关性 它度量空间灰度共生矩阵元素在行或列方向上的相似程度，因此，相关值大小反映了局部灰度相关性。 （6） 其中： ； ； ； 。 3．三维地震纹理体的提取原理 3.1 纹理基元的引入 地震数据体的振幅用灰度来表示，并按一定的窗去截取数据，可构建地震纹理基元（Texture element）Figure 1 the chart of seismic texture element Figure 2 the different directions of seismic texture element volume 3.2 地震资料灰度共生矩阵的构建 地震资料灰度共生矩阵可在地震纹理元基础上构建，它反映了地震纹理体在方向、间隔、变化幅度方面的综合信息。对于图1所示的纹理基元体，其灰度共生矩阵的矩阵元素值可通过下式来计算，即： （7） 式中，为纹理基元体中点的灰度值；和表示距离为δ的两个像素点；∑表示纹理体中沿着某一方向距离为δ的任意两个像素点满足以上条件的概率统计；,表示方向，特殊地，=90°，=0°对应X轴方向；=0°，=0°对应Y轴方向；=90°对应Z轴方向；L表示灰度级数。 4．实际资料应用及分析 选取胜利油田JJD工区作为应用实例，该区域断裂比较发育，断层和裂缝对油气储集和疏导起主要的作用，因此有效地识别断层和裂缝具有重要的研究价值。计算时采用滑动纹理窗，提取的属性包括纹理能量、对比、熵和纹理相关性等四种属性。此外，还讨论了不同的计算方向、时窗大小和灰度级别等应用要素问题。 4.1 不同纹理属性的对比分析 图3(a)为某工区一沿层振幅切片，为了检测断层，用纹理属性来进行分析。通过试验，采用16级灰度级别和9×9×9的三维滑动窗提取纹理体属性，然后抽取沿层切片进行对比分析。这里提取了纹理能量属性（图3(b)），纹理对比度属性（图3(c)），纹理熵属性（图3(d)）和纹理相关性属性（图3(e)）。对比图3种纹理属性我们发现，以上四种属性基本上都能较好地识别南部大的断层，中部多条断层层次感也比原始切片要好，纹理对比度和相关属性还较好地指示了北部的向斜构造。 （a）原始切片