Hilbert_Huang变换与地震信号的时频分析.docVIP

[文章编号]100124683(2005022207209 [收稿日期]2004206208;[修定日期]2005203230。 [项目类别]北京市自然科学基金项目(8041001、地震科学联合基金项目(604022、中国地震局三结合项目。 [第一作者简介]武安绪,男,生于1967年,副研究员,研究方向为地震预报、地震波形处理与应用。 Hilbert 2Huang 变换与地震信号的时频分析 武安绪1,2　吴培稚1　兰从欣1　徐　平1,3　林向东1 1北京市地震局,北京市苏州街28号　100080 2中国地震局地球物理研究所,北京　1000813中国科学院地质与地球物理研究所,北京　100029 摘要　本文介绍了HHT 时频分析方法及瞬时频率的概念,给出了已知信号的经验模态分 解及其时频分布,并对实际地震波形信号进行了HHT 时频处理与剖析。结果表明,HHT 方法能准确描述地震波形信号的非线性时变特征,是地震信号时频分析的有效工具。 关键词:　H ilbert 2H u ang 变换　瞬时频率　地震波形　时频分析 [中图分类号]P315　　　[文献标识码]A 0　引言 随着数字化地震台网的不断建设,采用新方法对高精度、高采样率地震数据的分析研究具有重要的现实意义(吴书贵等,2003。地震波形是具有时变特性(或称非稳态性质的典型信号(沈萍等,1999;刘希强等,2000,对于这类信号,不仅需要从总体上了解它的频率成分,而且还需要了解每一时刻信号中所包含的频率成分。 目前对地震信号进行分析的主要工具是傅里叶变换(胡广书,1997;郑治真,1998、现代谱估计(皇甫堪等,2003、G abor 变换(郑治真等,1996;科恩,1988、Wigner 2Ville 分布(沈萍等,1999;科恩,1988;郑治真等,1993、小波变换(刘希强等,2000;章珂等,1996;李宪优等,1999;Mallat ,1989;Daubechies ,1988;C oifman ,1990等。傅里叶变换和现代谱估计是建立在稳态信号处理基础之上,它仅能给出信号总体所包含的各种频率成分,不能解决何时出现何频率的问题;G abor 变换是在傅里叶变换基础上发展起来的一种时频分析方法,它在时域和频域采用固定的分辨率,或对不同的频率成分,在时域上的取样步长不变,有时会出现G ibbs 现象;Wigner 2Ville 分布也是一种局部化的时频分析方法,它因交叉项的引入会导致时频平面上的伪影现象,使分辨率受到一定影响;小波理论与以往其他研究方法相比,在频谱分辨率和时间定位精度方面都得到了显著提高,但小波变换的有效性依赖于小波函数的选取,有时还会存在随着尺度增大相应正交基函数的频谱局部性变差的缺陷,使其对信号的更精细分解受到了一定的限制(秦前清等,1995;李世雄等,1995。 第21卷　第2期(207~215 2005年6月中国地震E ARTH QUAKE RESE ARCH I N CHI NA V ol.21　N o.2Jun.2005 近年发展起来的希尔伯特-黄变换(Hilbert 2Huang T rans form ,HHT 是一种新的非线性信号处理方法(Huang et al.,1998;Y ue et al.,2001;Deng Y ongjun et al.,2001。HHT 方法将信号作经验模态分解(Em pirical M ode Decom position ,E MD ,能有效地将信号的各种频率成分以固有模态函数(Intrinsic M ode Function ,I MF 形式从时间曲线中分离出来。不同的I MF 分量是平稳信号或简单的非线性信号,具有简单的非线性特征。其缓变波包特征意味着不同特征尺度波动的波幅随时间变化,因而也具有时间上的局域化特征,而I MF 则属于窄带信号,正好满足Hilbert 变换(胡广书,1997的要求。对I MF 序列进行Hilbert 变换,得到包含时间、频率、振幅的三维离散骨架谱,可提供更加清晰的局部细节时频特征。HHT 方法具有较好的客观性、内在性与自适应性,对信号的非线性反映能力较好,适合于对具有非线性和非平稳动态变化特征的地震信号的描述与刻划。 本文试图在地震波形分析与研究中引入HHT 时频分析方法。首先通过对已知信号的仿真试验,探讨该方法的有效性;然后对地震波形信号进行多尺度经验模态分解和希尔伯特变换,剖析地震波形在不同I MF 尺度上的变化特征,研究地震波形在时频谱上的演化特性;最后讨论HHT 方法的适用性及其对非平稳动态变化的地震波形信号的时频刻划能力。1　基于HHT 方法的时频分析 HHT 方法由Huang 变换和Hilbert 变换两部分组成。 111　H uang 变换 Huang