基于小波包变换与近似熵的储集层识别.pdfVIP

第 15卷 第6期 重庆科技学院学报 自然科学版 2013年 l2月 基于小波包变换与近似熵的储集层识别 裴小刚 陈学华 侯 磊 邓吉刚 i．成都理工大学地球物理学院，成都 610059；2．川庆地球物理勘探公司物探研究中心，成都 610213 摘 要：通常不同地层特征对应的地震波复杂程度均有差异。近似熵是一种反映信号复杂程度的特征量，对储层特 征有一定反映。但是，地震波是有一定带宽的信号，在此带宽内，可能某一部分频段对油气层较为敏感，而其他频段 的信号则对油气层的反应不够敏感，如果直接进行近似熵的计算，其结果往往不能很好地反映储层特征。鉴于这种 情况，对小波包变换和近似熵结合进行储集层预测的方法进行研究，即对小波包分解的结果进行近似熵计算。实际 资料计算结果表明，结合小波包分解和近似熵的储层识别方法，能够精确刻画储层形态。 关键词：小波包变换；近似熵；高频；储集层识别 中图分类号：TEl9：P628 文献标识码：A 文章编号：1673—1980 2013 06—0062—04 地面接收到的地震波波形包含丰富的地下地质 可积函数构成的函数空间L2 R 的标准正交基库， 信息，通常储集层 内部的非均匀性越强时，地震波 从此库中可以选出L R 的许多组标准正交基 ，常 就会越复杂。因此，近年来有地球物理研究者 用的正交小波基是其中的一组。小波函数是小波 通过分析地震波波形的复杂性 以达到储集层检测的 包函数族 中的一个 ，所 以，小波包是小波函数的 目的，并且取得了良好的效果。然而地震波是频带 推广 J。 有限的信号 ，在这个有限的频带内，部分频段信 用 h 代表正交尺度函数 t 对应的正 号对储层特征反映较为敏感，其他频段的信号则对 交低通实系数滤波器， g 代表正交小波函数 储层特征的反应不够敏感 “』，甚至会造成干扰。 t 对应的高通滤波器，其中g ～1 h一，则 因此直接用地震道数据进行近似熵的计算，其结果 它们满足 以下两尺度方程和小波方程 ¨ ： 可能达不到理想的储层检测效果。对此本文对利用 f咖 ∑hkch 2t一 小波包分频的方法对地震道数据进行分解 J，然 J【 1 后利用分解后的分量的近似熵 开展储集层特性 砂 ∑gk 2一 KEZ 检测的方法进行了研究。考虑到小波变换只分解地 其中： t 为正交尺度函数， t 为正交小波函 震信号的低频部分，对地震信号的高频成分即细节 数。为了便于表示小波包函数，本文引入以下符号： 部分不做进一步分解，因此小波变换对 以低频信息 为主要成分的信号能较好地描述，而不能较好地描 f 2 述信号的细节信息 如较小的边缘、纹理等 。小波 t 包变换可以对高频信号进行无冗余、无疏漏分解，使 于是，式 1 可表示为： 信号能有更好的时间和频率分辨率 ，从而达到 f。 f ∑ 。 2一 3 精细刻画储层形态和边界的效果。 。 ∑g ￡ 2一 1 小波包变换原理 通过 。 、h、g在固定尺度的条件下可定义一 1．1 小波包变换的思想 组称为小波包的函数。 正交小波包是一个函数族 ，由它们构成平方 利用小波包对小波空间 进行更精细的分解。 收稿13期：2013—09—02 基金项目：国家 自然科学基金项 目 座机电话号码 ；国家科技重大专项项 目 2011ZX05023—005—010 。 作者简介：裴小刚 1989一 ，男，河南省濮阳人，成都理工大学在读硕士研究生，研究方向为储层预测。 · 62 · 裴小刚，等：基于小波包变换与近似熵的储集层识别 对任意非负整数 ∈Z，令 表示由小波包 。的二 复杂度越低。 进伸缩和平移2 。 2Jt—k ，k∈Z的线形组合生成 地震信号的近似熵计算步骤如下 ， ： 的 R 的闭子空间，则 第一步：用地震道数据构造一个相空间，设地 f vj，∈z … 震道信号离散采样长度为Ⅳ，时间序列为 ， ，…， I ： ，∈z ，重构长度为 Ⅳm、维数为m的相空间，重构的伪 相空间轨道为： 其中， 是由尺度函数tz。 生成的L R 的多 X Ix，X十【￡，…， 1，2，…，Ⅳm 9 分辨率分析，而 是由小波函数tz 咖生成的正 上一一—— ] ．二一 。—～一一树一～ 交小波子空间序列。利用新的符号，可将正交分解 式 9 中，为时间延迟，m为嵌人到欧氏空间 vjo ，J∈Z改写： ① ，∈Z，其 R 的维数，其中N N一 m一1 ‘。 + 第二步：对每一个 i，计算矢量与其余矢量之间 中 是由小波包 的二进伸缩和平移的线性组 的距离 合生成的L R 的闭子空间，①表示正交和，并且， d X i ，X j nmax f