地震图像的振幅共生矩阵分析及其应用.pdfVIP

学兔兔 石油天然气学报 (江汉石油学院学报) 2011年6月 第33卷 第6期 · 64 · JournalofOilandGasTechnology (J．JPI) Jun．2011 Vo1．33 No．6 地震 图像 的振幅共生矩阵分析及其应用 雷英成 (成都理工大学地球物理学院，四川成都610059) 曹 芝兴 (成都理工大学油气藏地质及开发工程国家重点试验室，四川成都610059) [摘要]介绍 了地震 图像的振幅共生矩阵分析方法及其在地震属性提取等方面的应用。详细讨论 了地震 图 像振幅共生矩 阵分析 的实现方法和 时一空窗等参数的选择。基于实际资料处理 的试验结果表 明，振幅共生 矩 阵分析获得 的均一性剖面能在一定程度上提高地震 图像 的纵 向分辨率 ，增强同相轴的连续性。 [关键词]地震 图像 ；振幅共生矩阵；均一性；分辨率 [中图分类号]P631．44 [文献标识码]A [文章编号]1000—9752 (2011)06—0064—05 地震图像纹理是指在地震数据的某个体元 (或像元)处，反射振幅与波形在水平方向和垂直方向的 变化情况L】]。振幅共生矩阵分析是地震图像纹理特征提取的一种常用方法，GaoⅢ、West等 把地震 图像纹理属性应用于地震相的划分，Chopraf。把地震图像纹理属性应用于油气识别等。高分辨率的地 震资料对薄层的解释 ，层位的标定 ，储层结构的识别等有其重要的意义。研究应用表明，振幅共生矩阵 分析获得的均一性剖面能在一定程度上提高地震图像的纵向分辨率 ，增强同相轴的连续性。 l 二维地震图像振幅共生矩阵分析原理 L 1 构建地震子剖面 经过偏移处理后的二维地震数据剖面，每个结点 (z．， )数据可用A(z．， )来表示 ，其中，X方向为测 线方向，y方向为时间方 向。每一个结点则刚好可以看 作为一个像元，对应于图中的每一个小正方形 (图 1)。 这样 ，一幅地震剖面即被划分为许多个像元，每个像 元的值为各地震道在每个采样点处的振幅值。 所谓地震子剖 面，是指 由有 限个像元组成的一个 二维数组 ，相对于像元足够 大，但相对地震剖面又足 够小的子剖面。首先在地震剖面上选 出一个像元 (图 1)，进而可 以以该像元 为中心 ，建立一个 (2m4-1)× (2n+1)的地震子剖面 (图1)。该地震子剖面将作为步 图 1 地震剖面与地震子剖面 骤 2的输入 ，用于建立振幅共生矩阵。 1．2 建立振幅共生矩阵 振幅共生矩阵是指从地震子剖面上振幅值为A 的像元 ( ，)出发，统计与其距离为 d、方 向为 a、振 幅值为 A 的像元 (p，q)，同时出现的频度 E(A ，A ，d，)，其数学表达式Ⅲ 为 ： ‘ E(A，A。，，a)一∑ {[(，)，(，q)==)(．)lA(m，”)=A，A(p，q)一A]} (1) 式 中，d取 1；a取0、45、90、135。。由此在每一个像元处便产生 4个振幅共生矩阵 ，这4个振幅共生矩阵便作 [收稿日期]2010—11—18 [基金项目]国家自然科学基金项 目 。 [作者简介]雷英成 (1986一)，男，2009年大学毕业 ，硕士生，现主要从事油气地球物理勘探与储层信息提取研究工作 。 学兔兔 第 33卷第 6期 雷英成等 ：地震图像的振 幅共生矩阵分析及其应用 为步骤 3的输入 。 1．3 计算二维地震图像均一性属性值 南振幅共生矩阵去评估该地震子剖面的均一性纹理属性，进而建立地震图像均一性属性剖面。均一性 (Homogeneity，记为 H)的表达式 为 ： HH：：∑∑I[ ]I。