基于时频谱熵的储层分布检测.pdfVIP

石油天然气学报 (江汉石油学院学报) 2010年4月 第32卷 第2期 ·66 · JournalofOilandGasTechnology(J．jPI) Apr．2010 Vo1．32 No．2 基于时频谱熵的储层分布检测 蔡涵鹏 ，贺振华 黄德济 ，李 瑞 (成都理工大学油气藏地质及开发工程国家重点实验室，四¨I成都。10059) [摘要]地震波在不同介质中传播表现出不同的频谱特征。目前，时频分析技术作为地震波频谱特征分析 的一种必要手段 ，在储层分布检测 中起着重要 的作用，但是地震数据体经谱分解后生成若干共频率谱数 据体 占用大量储存资源 ，而且使得解释成为一项繁重的工作。谱熵作为信息 的期望值 ，从熵的角度看来 ， 有序信号具有小的熵值 ，而无序信号具有大的熵值。因此为 了克服时频分析带来 的不便 ，提 出利用 s变 换完成时频谱分解 ，然后在时频谱的基础上进行谱熵运算 的时频谱熵方法。通过计算谱分解后频谱 的谱 熵 ，将若干共频率数据体合成一个 时频谱熵数据体，减少 了数据 的存储量 ，提高 了解释 的效率。在实 际 资料处理 中用于分析岩性、砂岩储层 的空间展布和地质构造等，取得 了良好 的效果。在实际资料应用 中 发现 ，砂岩储层的时频谱熵值最高’。‘ [关键词]时频分析；S变换；储层分布；谱熵；时频谱熵 、 [中圈分类号]P631．44 [文献标识码]A [文章编号]1000—9752(2010J0‘2—0066—03 传统的地震方法主要用于构造解释，然而地震波除了含有反射层构造信息外，还包含有许多与储层 性质相关的信息、。储集体中细微的构造、岩性及油气水状况的变化等 ，都不同程度地影响地震反射波的 能量、波形、频率以及相位等特征。地震波的频谱有时比波形本身更能反映出波动的·些特征，而这些 特征又与地层的岩性有着一定的关系。并且表征岩石物性的地层吸收衰减的品质因子 Q值越大，能量 衰减得越慢 ；Q值越小 ，能量衰减越快 。由于地层含流体 ，其品质因子较小 ，地震波在传播过程中，地 层对地震波的吸收也加强。高频地震波在地层传播的过程 中震动较快 ，使得 当地震波穿过含油气地层 时，高频能量衰减快 】『]，从而使能量相对集中到低频部分 ，导致地震波在不同的介质 中传播表现出不 同的频谱特征。时频分析技术给信号局部频谱分析指明了方向。然而对于大规模的地震数据 ，尤其是三 维地震数据 ，通过时频分析技术进行谱分解后得到的共频率瞬时能量谱数据体不仅 占用大量 的储存资 源，而且对于每一个频率的时间切片和沿层切片进行资料详细解释也是一件繁重的工作 ，不利于大规模 地震资料的处理和应用 。而谱熵特征作为信息的期望值 ，从熵的角度看来，有序信号具有小的熵值 ，而 无序信号具有大的熵值_3“]。不 同的地层 、岩性处的频谱特征具有不 同的有序度 ，因此通过计算该处频 谱 的谱熵来度量频谱的总体特征 ，可 以避免时频分析技术带来 的不利 。笔者在时频分析技术的基础上， 提出利用时频谱熵进行储层分布检测的方法 。应用该方法对实际地震数据进行处理后生成时频谱熵数据 体 ，将其用于砂岩储层的空间展布、地质构造和岩性分析，取得 良好效果。 1 算法原理 1．1 时频分析技术的原理 频谱分解技术是一种频率域的解释方法 ，是地震属性分析中的重要组成部分。到 目前为止 ，用频率 信息来研究储层仍然是地质 、地球物理工作者关注的焦点 引。时频分析最早是由短时傅里叶变换来实 现 ，然后发展到由小波变换来实现．，到 目前的S变换，乃至广义 S变换来实现。S变换和广义 S变换的窗函 [收稿日期]2009—04—13 [基金项目]国家 863”计划项目 (2o06AAo9A1o2—12) [fie者简介]蔡涵鹏 (1983一)，男，2006年大学毕业，博士生，现主要从事石油地球物理勘探方面的研 工作。 第 32卷第 2期 蔡涵鹏等：基于日寸频谱熵的储层分布检测 数随着频率的增加而 自适应地调整分析时窗 ，时频窗 口具有一定 的自适应性 ，不需要满足小波容许性条件 以及与传统傅里叶变换 的直接关系 。由于与短时傅里叶变换和小波变换等相 比较 ，S变换和广义 S变换 有其独特的优越性，因此该研究采用 S变换计算时频谱。一维连续 S正变换和反变换的表