基于Gabor变换的时频分析在气层识别中的应用.pdfVIP

石油天然气学报 2012年3月 第34卷 第3期 JournalofOilandGasTechnology Mar．2012 Vo1．34 No．3 基于 Gabor变换的时频分析 在气层识别 中的应用 唐湘蓉 (地球探测与信息技术 教育部重点实验室；成都理工大学地球物理学院，四JIl成都610050) 彭真明 (电子科技大学光电信息学院，四JII成都610054) 蔡涵鹏 (成都理工大学地球物理学院，四川成都610059) [摘要]地震波在地层传播过程中能量发生衰减，且高频能量衰减强，低频能量衰减弱，尤其在含气储层 中高频能量急剧衰减特征更为突出。借助于Gabor变换优 良的时间域与频率域局部化特性，将基于Gabor 变换的时频分析方法引入到地震波吸收衰减分析 中，对南海莺歌海盆地某 区块碎屑岩储层进行 了油气预 测，预测结果与实际钻井资料符合 ，证 明了该方法的有效性和在油气预测 中的应用前景 。 [关键词]Gabor变换；时频分析；能量衰减；碎屑岩储层；油气预测 [中图分类号]P631．44 [文献标识码]A [文章编号]1000—9752 (2012)03—0065一O6 地震波在地层传播过程 中，能量发生不 同程度地衰减 ，其波形和瞬时频谱均发生变化 。引起能量衰 减的原因分为两类 ：一类是与传播特性有关的各种散射、反射和透射 ；另一类是由于介质 自身非完全弹 性所引起的地震波振幅的吸收衰减。影响地震波吸收衰减的因素较多，但在相同地质 、沉积环境下，地 层结构相对较为稳定，纵、横 向岩性变化不大的条件下，地震波能量衰减主要由流体性质引起 。大量 的 正演模型和理论研究表 明[】 ，地震波在聚集 了天然气的储层 中传播时，地震波的高频成分多被储层 所吸收，往往伴随着高频能量的剧烈衰减 ，且频率越高的能量衰减越快。因此 ，如果能很好地检测含气 储层对地震波高频成分吸收衰减特征，就能预测油气储层分布。笔者提出了利用基于Gabor变换的时 频分析方法进行能量衰减分析与谱分解 ，预测结果跟实际钻井资料非常符合，取得了良好的预测效果。 1 基本原理 传统的傅里叶变换是以时间为变量的信号转换到频率域 ，虽然得到 了能量在各频率成分 的分布情 况，但缺乏时间信息；时间信号展示了信号随时间的变化情况，却缺乏能量在各频率成分的分布信息。 时频分析方法则兼具二者的优点，它建立了一个以时间和频率为变量的窗函数，通过窗函数的时移和频 移，在时间域和频率域对信号进行局部分析，从而得到信息在不同时刻的局部频谱特性。通过对信号进 行时频分析，就能够获得信号能量在时间域和频率域的分布情况，以达到对信号细微特征的解剖和更全 面 的认识 。 Gabor变换是 Heisenberg不确定准则下的最优短时傅里叶变换 ，是具有高斯窗函数的短时傅里叶 变换 ，高斯窗函数是短时傅里叶变换时间分辨率和频率分辨率最优的窗函数。与传统的时频分析方法 (如短时傅里叶变换)相 比，Gabor变换具有更好的时间一频率分辨率 ，在诸如地震信号的非平稳信号处 理中，Gabor变换发挥着积极的作用 。 对于样点个数为N 的离散信号z()一 {。， ，z。，…，z~-。)，其 Gabor展开式为： [收稿 日期]2011—10—27 [基金项 目]国家自然科学基金项 目 40839905)。 [作者简介]唐湘蓉 (1967一)，女，1989年大学毕业，副教授，博士生，现从事油气地震勘探信号处理及综合解释工作。 石油天然气学报 石油物探与测井 2012年 3月 z()一∑∑G (n) (1) = 0 = 0 ，女 ()一 fi(n—mAM)e~k 一 2~kAK／N (2) 式 中，h()是综合窗函数^(n)以N为周期的周期延拓 ； 、k分别为时域 、频域的采样点序号 ；M、K分别为 时域 、频域的采样点个数 ；G ，为 ()的Gabor变换 ；AM、AK分别为时域 、频