基于广义S变换的地震资料谱分解技术研究.pdfVIP

维普资讯 石油天然气学报 (江汉石油学院学报) 2008年6月 第 30卷 第3期 JournalofOilandGasFechnology(J．JP1) Jun．2008 Vo1．30 No．3 基于广义 S变换的地震资料谱分解技术研究 郭爱华 (冀东油田分公司勘探开发研究院，河北唐山063004) 路鹏飞 ，刘 鹏 (中国科学院地质与地球物理研究所，北京100029) 张会来 (中国石油大学资源与信息学院，北京102249) [摘要]谱分解技术一般采用短时Fourier变换和连续小波变换。短时Fourier变换的分析时窗是恒定的， 不随频率的变化而 自动调节，连续小波变换 的尺度与频率 的关系并不确定 。广义 S变换具有优 良的时频 聚集性和 时频分辨率，它综合 了短 时Fourier变抉和连续小波变换的一些优 点，而克服 了它们 的不足 ，同 时通过调节 日寸窗宽度 因子可 以获得更好 的时频分辨率。研究将其应用于地震资料 的谱分解 中，给 出了基 于广义 S变换的谱分解 业务 流程 ，并对实 际地 震资料进行 了分析 。结果表明该方 法是可行 的。 [关键词]谱分解；时频分析；连续小波变换；广义 S变换 [中图分类号]P631．44 [文献标识码]A [文章编号]1000—9752 (2008)03—0083一O4 谱分解技术是近年来发展起来 的一项基于频率谱分解的储层特色解释技术 ，是地震属性分析 中的重 要组成部分L1]。1999年，BP．AI]IOCO石油公司 的GregPartyka等人利用短 时 Fourier变换 (Short TimeFourierFransform，STF~I、)，研究在短时窗内通过频谱分解来研究薄层变化 ，并提出了谱分解概 念L5]。2001年 ，Landmark在此基础上推 出了 SpectraDecom工具 ；随后 ，Schlumberger公司与 Land～ mark公司合作在 GeoFrame软件 中加入 了 SpectraI)ecom 工具；2004年 ，Paradigm 公司利用 Gabor Morlet小波变换推出了SpectralDecomposition工具 。地震波在地层介质 中传播 ，由于散射和吸收衰减 等均与频率有关，严格意义上属非平稳信号。由于Fourier分析使用的是一种全局变换，无法表达信号 的时频局部性质，而这种性质恰恰又是非平稳信号最根本和关键 的性质。由此，研究人员分别应用短时 Fourier变换、Gabor变换 ]、连续小波变换 (ContinuousWaveletTransform，CWT)对谱分解 中的 频谱分析方法进行改进L7一。1996年，R．G．Stockwell综合短时Fourier变换和小波变换的特点，提出 了一种非平稳信号分析和处理的方法——广义 S变换，它是以Morlet小波为基本小波的连续小波变换 的延伸[8]。同小波变换、短时Fourier变换等方法相比，广义 s变换有其独特的优点，如信号的广义 S 变换的时频谱分辨率与频率相关 ，不但有多尺度聚焦性 ，而且与其 Fourier谱直接相联系，保持频率的 绝对相位 ；其基本变换 函数不必满足容许性条件等。同Wigner分布 、Cohen等双线性变换相 比，广义 S变换是线性的，不存在交叉项 ，更有利于分离信号 ，特别是存在噪声 的情形_g 。谱分解技术 的核心 模块为地震信号的时频分析，笔者将广义 s变换应用于地震资料的谱分解 中，米提高刻画和识别特殊地 质体的水平 。 1 广义S变换基本原理 设 ()∈I。(R)，(L。(R)表示实数域上的平方可积函数空间)，^()相对于分析小波 硼()的广义 s 变换定义为 ： s(r，厂)一f。。(￡)×gf(￡一r)×exp(一2nifofi)d￡ (1) J … 式 中，S(r，_厂)是 ^()的广义 S变换 ；gr(￡一r)一A {．厂!exp[一 f(t—r)-p)]是高斯窗函数 ；t是时问； [收稿 日期]2007—12一O7 [作者简介]郭爱华 (1981一)，女，2003年大学毕业，工程师 ，现主要从事地震地质综合解释方面的研究工作。 维普资讯