地震资料的Hilbert—Huang变换及其应用.pdfVIP

石油天然气学报 2012年 12月 第34卷 第 12期 JournalofOilandGasTechnology Dec．2012 Vo1．34 No．12 地震资料的Hilbert—Huang变换及其应用 刘浩杰 (q]石化胜利油田分公司物探研究院，山东东营257022) [摘要]地震信号能量和频率随着时间的推移而改变，仅仅从时间域或频率域来分析信号已不能充分地描 述地震信号的特征 。基于 Hilbert—Huang变换 (HHT)，把非稳态 的地震信号分解为不 同分辨尺度 的基本 模态分量，利用不 同模态分量 的瞬时时频特征精细刻画地震信号的时间一频率一能量之间的关系，探讨 了 HHT在面波压制方面的应用。应用结果表明，HHT在时频域可以较好地达到衰减面波的 目的。 [关键词]地震信号；Hilbert-Huang变换 (HHT)；时频分析；时间分辨率；频率分辨率 [中图分类号]P631．44 [文献标志码]A [文章编号]1000—9752(2012)l2—0076一O5 地震信号具有短时 突变等特点，其能量和频率随着时间的推移而变化 ，仅仅从时间域或频率域来 分析信号已不能充分地描述信号的特征。时频分析技术将地震数据变换到二维时间一频率域，得到信号 的频谱 、能量特征随时间的变化 ，在时间一频率域 中研究地震资料的各种特征。时频分析方法有很多， 主要有短时傅里叶变换 、小波变换、S变换以及广义s变换等方法 ]。传统的傅里叶变换建立在稳态 信号处理基础上 ，分析信号总体所包含 的各种频率成分 ，但不能给出不 同频率成分的时间局部化信息； 短时傅里fl-r变换 以傅里叶变换为基础 ，通过窗函数作用，分析固定分辨率条件下信号的时频特征 ；小波 变换在时域和频域都具有很好的局部化性质，提高了频谱分辨率和时间定位精度，但本质上仍是一种窗 口可调的傅里叶变换 ；S变换和广义 S变换改进了小波的形式，分辨率与频率 (尺度)有关 ，但本质上 仍是窗 口可调的傅里叶变换 ，分辨率受窗口长度的限制。常规时频分析方法在分辨率和信号平稳性方面 的限制影响了对地震信号时频特征刻画的精确性。 Hilbert—Huang变换 (HHT)是一种新兴的时频分析方法 ，被认为是近年来对 以傅里叶变换为基 础进行线性和稳态谱分析 的一个重大突破 。它最初是 由美 国国家宇航局的 Huang_l4于 1998年提出的。 该方法不受傅里叶分析的全局性限制，将任意信号依据信号本身的时问尺度特征，进行多尺度时频特征 的刻画，得到的时频谱有较高的时间分辨率与频率分辨率，也使得信号瞬时频率具有实际的物理意义。 HHT已应用到地震力学、地球物理学等领域，并取得 了较好的效果[ 。该次研究首先介绍了HHT 的原理 ，研究了HHT频谱和时频谱 ，探讨 了HHT变换在面波压制等方面的应用。 1 Hilbert—Huang变换 (HHT) HHT包括两部分l_4：第 】部分是经验模态分解 (empiricalmodedecomposition，EMD)；第 2部分 是 Hilbert谱分析 (Hilbertspectrumanalysis，HSA)。 为了利用 HHT来分析非线性和非平稳信号以获得有明确物理意义的瞬时频率，必须保证信号处理 方法 的完备性 、正交性、局部性和 自适应性 。基于此 ，Huang提 出了本征模态 函数 (intrinsicmode function，IMF)的概念 。本征模态函数必须满足 2个条件 ：一是在信号的有限时间内，其零交点的数 目与极值点的数 目必须相等或者至多相差一个 ；二是在任意时刻，由极大值得到的上包络线和由极小值 得到的下包络线的平均值为零 。非平稳信号并不满足 IMF分量 的两个条件 ，因此必须先将信号分解成 符合上述条件的 IMF分量，然后再对各分量进行 HHT。基于 IMF的定义 ，Huang提出了一种将任意 [收稿日期]2012—08—11 [基金项目]国家 8“63 计划项目 (2011AA060302)。 [作者简介]刘浩杰 (1976一)，男，1998年江汉石油学院毕业，博士，研究员，现主要从事油藏综合地球物理研究。 第 34卷第 12期 刘浩杰 ：地震资料的 Hilbert—Huang变换及其应用