匹配逼近法频谱分析效果对比分析.pdfVIP

学兔兔学兔兔 石油天然气学报 (江汉石油学院学报) 2009年8月 第31卷 第4期 Journal of Oil and Gas Technology(J．JPI) Aug．2009 Vo1．3 1 No．4 匹配逼近法频谱分析效果对比分析 张铁强 (中国地质大学(北京)地下信息探测技术与仪器教育部重点实验室，北京1 00083) 张 崴 (中国石化石油勘探开发研究院，北京1 00083) 牛滨华 (中国地质大学(北京)地下信息探测技术与仪器教育部重点实验室，北京1 00083) [摘要]在地震数据处理中，频谱分解是一个重要的信号处理工具。简要介绍了窗口Fourier变换 (WFT)，连续小波变换 (CWT)和匹配逼近 (Matching pursuit，MP)等3种频谱分析方法，比较了它 们的特点，并提出一种实现匹配逼近分解的计算方法。用合成信号对这3种方法进行试算，结果表明了 匹配逼近方法具有良好的时频特性，验证了本研究提出的实现匹配逼近算法的有效性。 [关键词]匹配逼近；频谱分析；子波 [中图分类号]P631．44 文献标识码]A [文章编号]1000—9752(2009)O4—0247—04 时间域信号序列的频率域分布通常表现出很多时间域难以观察的特征。时间序列转换到频率域后， 可以从其频率特征获得大量的新信息。Fourier变换是地震数据处理方法的基础，如频率滤波L1]。但是， 这些特征是标量参数，不能描述每一个时间点的频谱属性。 Fourier变换的振幅谱反映了不同频率的存在，但是却不能表示这些频率的时间分布。如果假定在 一 个小时窗内的信号是平稳的，那么它的Fourier变换展示在这个时间周期内信号的频率成分。wFT (窗口Fourier变换)方法[2 是在时间域用移动的时窗来计算Fourier谱，那么时间和频率的二维分布就 产生了，这个二维分布就是WFT。一旦为WFT选定了窗口函数，时频分辨率也就固定了，因此WFT 计算的时频谱局部特征较差，且受时间带宽分辨率的限制[3]。与WFT固定的时窗不同，CWT(连续 小波变换)的时窗是变化的。小波分析所使用的一系列函数在时间域是紧密支架，在频率域是带限 的[4 ]。这些窗口函数像细小的波，在短时间内发生和衰减，因此被称为小波。小波分析如同数学显微 镜。当一个窗口函数在时间上的宽度很大，则能得到信号的粗略特征。随着尺度变化，分析窗I：1变窄， 信号细微特征的分辨得到增强[6]。CWT的频率参数与尺度参数成反比例，因而不能提供精确的频率特 征 尤其在高频情况下[7卅]。MP(匹配逼近)的出现，目的在于克服WFT和CWT的缺点L1。。。 MP是一种新算法，它能将任何信号分解成一系列小波形函数的线性展开，这些小波函数是从大量 的函数字典中挑选出来的[1o]。对小波函数挑选的依据是它们能最佳地匹配信号特征。MP是一个迭代分 解的过程，传统的实现方法计算量很大，研究中用两步法实现MP算法。笔者简要介绍了MP方法的原 理，用3种方法对理论数据进行了试算和比较。 1 匹配逼近 (MP)方法 1．1 用Morlet函数作子波基的MP方法 MP将信号展开为一系列子波，这些子波是从函数字典中挑选出来的 ]。借助Gabor子波 ]，MP 对信号厂(￡)进行分解的公式如下： N—l 厂( )一∑a g ( )+R 厂 (1) 一 O [收稿日期]2009一O3一l8 [作者简介]张铁强 (1976一)，男，1999年大学毕业，博士生，现主要从事应用地球物理方法研究工作。 学兔兔学兔兔 石油天然气学报 (江汉石油学院学报) 2009年8月 式