希尔伯特_黄变换谱及其在地震信号分析中的应用
- 格式:pdf
- 大小:170.25 KB
- 文档页数:5
希尔伯特黄变换及其应用希尔伯特黄变换及其应用希尔伯特黄变换(Hilbert-Huang Transform,HHT)是一种用于分析非线性和非平稳信号的方法,它由黄其森(Norden E. Huang)和希尔伯特(Hilbert)共同提出。
该方法通过将信号分解为一组固有模态函数(Intrinsic Mode Functions,IMF)来提取信号中的模式和趋势。
本文将介绍希尔伯特黄变换的应用,并详细讲解其中的几个应用领域。
应用一:信号处理•希尔伯特黄变换可以用于音频信号处理,通过提取信号的固有模态函数,可以分离出音频信号中的主要频率成分,从而实现去噪、降噪等处理。
•在图像处理中,希尔伯特黄变换可以用于边缘检测和纹理分析。
通过提取图像的固有模态函数,可以分离出图像中的纹理信息和边缘信息,从而实现图像增强和分割等操作。
应用二:地震学•地震学中的信号分析是一项重要的任务,希尔伯特黄变换可以用于地震信号的分析和处理。
通过将地震信号分解为固有模态函数,可以提取出地震信号中的地震波的时频特征,从而实现地震信号的分类和识别。
•希尔伯特黄变换还可以用于地震信号的时频谱分析,通过将地震信号分解为固有模态函数,并对每个分量进行傅里叶变换,可以得到地震信号的时频谱图,从而更好地理解地震信号的时频特性。
应用三:医学工程•在医学工程中,希尔伯特黄变换可以用于生物信号的分析和处理,如心电图(ECG)和脑电图(EEG)等。
通过将生物信号分解为固有模态函数,可以提取出信号中的重要特征,如心跳频率、脑电波的频率等,从而实现疾病的诊断和监测。
•希尔伯特黄变换还可以用于生物信号的时频谱分析,通过将生物信号分解为固有模态函数,并对每个分量进行傅里叶变换,可以得到信号的时频谱图,从而更好地分析信号的时频特性。
应用四:金融市场•在金融市场中,希尔伯特黄变换可以用于股票价格的分析和预测。
通过将股票价格分解为固有模态函数,可以提取出股票价格的趋势和周期成分,从而更好地预测股票价格的走势。
Hilbert_Huang变换与地震信号的时频分析[⽂章编号]100124683(2005)022207209[收稿⽇期]2004206208;[修定⽇期]2005203230。
[项⽬类别]北京市⾃然科学基⾦项⽬(8041001)、地震科学联合基⾦项⽬(604022)、中国地震局三结合项⽬。
[第⼀作者简介]武安绪,男,⽣于1967年,副研究员,研究⽅向为地震预报、地震波形处理与应⽤。
Hilbert 2Huang 变换与地震信号的时频分析武安绪1),2) 吴培稚1) 兰从欣1) 徐 平1),3) 林向东1)1)北京市地震局,北京市苏州街28号 1000802)中国地震局地球物理研究所,北京 1000813)中国科学院地质与地球物理研究所,北京 100029摘要 本⽂介绍了HHT 时频分析⽅法及瞬时频率的概念,给出了已知信号的经验模态分解及其时频分布,并对实际地震波形信号进⾏了HHT 时频处理与剖析。
结果表明,HHT ⽅法能准确描述地震波形信号的⾮线性时变特征,是地震信号时频分析的有效⼯具。
关键词: H ilbert 2H u ang 变换 瞬时频率 地震波形 时频分析[中图分类号]P315 [⽂献标识码]A0 引⾔随着数字化地震台⽹的不断建设,采⽤新⽅法对⾼精度、⾼采样率地震数据的分析研究具有重要的现实意义(吴书贵等,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)等。
基于希尔伯特-黄变换的地震信号时频谱分析
基于希尔伯特-黄变换的地震信号时频谱分析
希尔伯特-黄变换(Hilbert-Huang Transform,HHT)是分析非线性、非稳定信号的一种新方法,能清晰地刻画地震信号的时频能量分布.首先将信号分解为有限数量的固有模态函数IMF,再对这些IMF求解瞬时频率,进而获得信号的时频谱.应用理论模型和实际地震道数据进行了试算,并与S变换谱进行了对比,证明该方法比S变换具有更好的`时频域刻画能力.对实际二维地震剖面做HHT变换求得希尔伯特谱,提取分频剖面分析认为,HHT瞬时谱具有一定的油气检测能力.
作者:侯斌桂志先胡敏王鹏陈小军 Hou Bin Gui Zhixian Hu Min Wang Peng Chen Xiaojun 作者单位:长江大学油气资源与勘探技术教育部重点实验室,湖北,荆州,434023 刊名:勘探地球物理进展英文刊名: PROGRESS IN EXPLORATION GEOPHYSICS 年,卷(期): 2009 32(4) 分类号: P631.4 关键词:希尔伯特-黄变换固有模态信号经验模态分解时频谱。
油气地球物理2015年10月PETROLEUM GEOPHYSICS第13卷第4期收稿日期:2015-07-20;改回日期:2015-08-25作者简介:冯玉苹,女,工程师,地球探测与信息技术专业硕士研究生,主要从事地球物理方法研究及软件开发工作。
基金项目:中国石化科技攻关项目(P12078)资助。
时频分析以联合时频分布的形式表示信号的特性,克服了一维傅里叶分析时域和频域完全分离的缺陷,可以准确定位某一时刻出现哪些频率分量,以及某一频率分量分布在哪些时刻上[1]。
目前,比较典型的时频分析方法有短时傅里叶变换法、小波变换法、S 变换法、改进S 变换法、Wigner-Ville 时频分布法、匹配追踪和Cohen 类时频分析法等。
时频分析方法在地震资料高分辨率处理、含油气性检测及油气藏识别、低幅度构造识别、地层层序划分、沉积旋回划分、随钻地震等方面取得了良好的应用效果[2-6]。
Huang 等[7]创立了Hilbert-Huang 变换(Hil-bert-Huang Transforms,HHT ),提出了经验模态分解(Empirical Mode Decomposition,EMD )和固有模态函地震数据的Hilbert-Huang 变换与瞬时属性提取冯玉苹(中国石化石油工程地球物理公司胜利分公司,山东东营257086)摘要:常规Hilbert-Huang 变换法中固有模态函数分量包含太宽的频率范围,且在低频区域会产生错误的频谱特征,不能分离出低能量信号成分,为克服常规方法存在的问题,须将小波分解法与Hilbert-Huang 变换法相结合。
首先利用小波分解方法将地震数据分解成窄带信号,然后进行经验模态分解,基于小波分解的Hilbert-Huang 变换方法充分发挥两种时频分析方法的优势,有利于提高时频分析和瞬时属性提取的精度。
将改进的Hil-bert-Huang 变换方法应用到模拟信号和实际地震数据处理中,时频分析和地震属性提取结果表明,二维时频谱具有较好的时频分辨率,瞬时地震属性具有明确的物理意义,在断层预测和储层精细描述中具有较高的应用价值。
希尔伯特黄变换及其在地震资料分析处理中的应用1.信号分解和频谱分析:希尔伯特黄变换能够将地震信号分解成不同的本征模态函数,这些函数表示了信号在不同频率上的振动成分。
通过分析每个本征模态函数的频谱,可以获得地震信号在不同频率段上的能量分布和频率特征,进而揭示地震事件的产生机制和地下构造。
2.成分提取和分析:地震数据通常包含多个成分,如地震波到达时间、振幅、频率等。
通过希尔伯特黄变换,可以将地震信号中的每个成分分离出来,并进行独立的分析。
例如,可以将地震波到达时间作为一个本征模态函数提取出来,用于地震定位和速度模型构建。
3. 时频特性分析:由于地震信号往往具有非线性和非稳定特性,传统的Fourier变换无法很好地揭示其时频特性。
而希尔伯特黄变换能够在时间和频率领域上同时揭示信号的特征,可以更准确地描述地震信号的时频演化过程。
4.脉冲识别和信号提取:地震信号中常常存在短时脉冲和异常事件,如地震的震相到达和异常地震事件。
希尔伯特黄变换可以对这些脉冲进行精确的提取和识别,从而帮助地震学家定位和研究这些地震事件。
5.地震预测和监测:希尔伯特黄变换可以对地震信号的非线性特性进行全面的分析,探索地震事件的前兆信号和变化规律。
通过对地震信号的长期监测和分析,可以预测地震活动的趋势和可能的突发事件,为地震预警和防灾减灾提供支持。
总结来说,希尔伯特黄变换在地震资料的分析处理中能够提供更准确和全面的信号分析结果,从而揭示地震事件的时频特性、构造演化和异常变化。
它在地震领域的应用有助于地震研究、地震监测和地震预测等方面的工作。
标题:hilber-huang 变换应用实例: 从信号处理到气候预测一、引言hilber-huang 变换(Hilbert-Huang Transform,简称HHT)是一种新兴的信号处理和数据分析方法,它结合了经验模态分解(Empirical Mode Dposition,简称EMD)和希尔伯特变换(Hilbert Transform)两种技术,被广泛应用于不同领域的数据处理和分析中。
本文将从信号处理到气候预测的多个实际应用实例,介绍hilber-huang 变换在各个领域的重要性和价值。
二、hilber-huang 变换在地震信号处理中的应用1. 地震信号特点分析地震信号是一种典型的非线性和非平稳信号,传统的频域分析方法往往难以准确捕捉其特征。
hilber-huang 变换通过EMD将地震信号分解成若干个本征模态函数(Intrinsic Mode Functions,简称IMF),再通过希尔伯特变换进行精确的瞬时频率分析,可以更好地分析地震信号的频率特性和瞬时变化规律。
2. 地震预测模型改进基于hilber-huang 变换的地震信号处理技术,可以为地震预测模型提供更精确、更全面的地震监测数据。
利用该方法分析地震监测数据,可以提高地震预测的准确性和可靠性,为地震防灾减灾工作提供重要支持。
三、hilber-huang 变换在金融领域的应用实例1. 股市交易信号提取股市交易数据具有复杂的非线性和非平稳特性,传统的技术分析方法难以全面把握市场的变化。
hilber-huang 变换的EMD分解可以将股市交易数据分解成不同的频率成分,希尔伯特变换可以提取出股市交易信号的瞬时特征,为投资者提供更准确的交易决策依据。
2. 金融时间序列预测利用hilber-huang 变换的信号处理技术,可以更准确地分析金融时间序列数据的周期性和趋势性,提高金融预测模型的准确度和鲁棒性,为金融市场参与者提供更可靠的预测和决策支持。
基于希尔伯特—黄变换的汶川地震记录分析裴强;胡波【摘要】利用HHT方法对6个台站的汶川8.0级地震加速度记录进行EMD分解,进而对地震波输入结构的总能量进行分解,得到地震动能量在IMF分量之间的分布关系.同时通过对各分量边际能量谱的分析,可以从另一方面得到地震波的频谱特性,边际能量谱峰值频率与地震波傅立叶谱的卓越频率相同,边际能量谱也可以作为一种估算震动卓越周期的方法.%Using the Hilbert-Huang Transform, acceleration time histories recorded by 6 stations in Wenchuan M08. 0 earthquake are decomposed into several Intrinsic Mode Function components by Empirical Mode Decomposition. The total energy of structure inputted by seismic wave is decomposed to get the distribution relation of the ground motion energy in the different IMF component By analysis of marginal energy spectrum of each IMF component, we can obtain the spectral characteristics of seismic wave. The frequency of marginal energy spectrum peak is the same as the predominant period of the Fourier spectrum of seismic wave, so it can also be used as the way to estimate the predominant period of vibration.【期刊名称】《地震研究》【年(卷),期】2012(035)002【总页数】5页(P246-250)【关键词】Hilbert-Huang Transform;经验模态分解;能量分布;卓越频率;汶川8.0级地震【作者】裴强;胡波【作者单位】大连大学土木工程研究与技术开发中心,辽宁大连116622;大连大学土木工程研究与技术开发中心,辽宁大连116622【正文语种】中文【中图分类】P315.9采用恰当的数据分析方法,从地震记录中提取地面运动特性是非常重要的:作为模拟地震动的输出和地震现象的解释,有助于对地震震源机制、在非线性土壤中传播等问题的理解;作为岩土和结构工程系统的输入,可以用来计算结构的动力非线性反应,从而评估结构的抗震性能。
地震信号时频分析中的希尔伯特黄变换研究
周竹生;罗勇涛
【期刊名称】《物探化探计算技术》
【年(卷),期】2016(038)001
【摘要】随着时频分析方法的发展,生产研究上对复杂信号的时频分析有了更高的要求.这里简要介绍了希尔伯特—黄变换的原理和实现步骤,然后对合成信号进行经验模态分解(EMD)和希尔伯特谱分析,进而将HHT方法运用到实际地震信号的时频分析中.在此基础上,运用经验模态分解对合成地震信号进行了阀值去噪,证明了该去噪方法的有效性.
【总页数】8页(P59-66)
【作者】周竹生;罗勇涛
【作者单位】中南大学地球科学与信息物理学院,长沙 410083;中南大学地球科学与信息物理学院,长沙 410083
【正文语种】中文
【中图分类】P631.4
【相关文献】
1.基于小波变换算法的时频分析技术在地震信号处理中的应用 [J], 杨立强
2.基于希尔伯特-黄变换的地震信号时频谱分析 [J], 侯斌;桂志先;胡敏;王鹏;陈小军
3.希尔伯特—黄变换在地震信号时频分析中的应用研究 [J], 冯红武;王建昌
4.旋转机械振动信号基于EMD的希尔伯特变换和小波变换时频分析比较 [J], 杨
世锡;胡劲松;吴昭同;严拱标
5.希尔伯特-黄变换地震信号时频分析与属性提取 [J], 杨培杰;印兴耀;张广智因版权原因,仅展示原文概要,查看原文内容请购买。
基于Hilbert-Huang变换的地震瞬时属性提取方法及应用黄亚平;耿建华;钟广法;郭彤楼;蒲勇;丁孔芸;麻纪强【期刊名称】《应用地球物理(英文版)》【年(卷),期】2011(008)002【摘要】希尔伯特—黄变换(Hilbert-Huang transform,HHT)是一种新的适合非平稳和非线性信号的分析方法,由于地震信号一般呈现出非平稳与非线性特性,因此HHT非常适合地震信号的分析。
本文首先介绍了HHT中关于经验模态分解(Empirical Mode Decomposition,EMD)的实现过程,在此基础上分析了几种基于EMD获得本征模态函数(Intrinsic Mode Functions,IMF)来计算瞬时频率的算法,其中利用了两个采样间隔瞬时频率的平均来计算瞬时频率,较好地反映了地震信号频率成分随时间变化的特征。
将该方法应用于四川东北部某地区海相碳酸盐岩地层三维地震叠后偏移数据处理,提取"三瞬"地震属性,与传统的希尔伯特变换提取的"三瞬"地震属性进行对比,结果表明基于HHT的"三瞬"地震属性结果具有更高的分辨率,IMF2的瞬时相位能够较好地刻画台地边缘生物礁相,IMF2的瞬时频率亦具有较好的分带性。
将IMF2的"三瞬"地震属性与钻井等资料结合分析,能够更好地识别沉积相的分布。
【总页数】9页(P125-133)【作者】黄亚平;耿建华;钟广法;郭彤楼;蒲勇;丁孔芸;麻纪强【作者单位】同济大学海洋地质国家重点实验室海洋与地球科学学院,上海200092;同济大学海洋地质国家重点实验室海洋与地球科学学院,上海200092;同济大学海洋地质国家重点实验室海洋与地球科学学院,上海200092;中国石化勘探南方分公司,成都610041;中国石化勘探南方分公司,成都610041;同济大学海洋地质国家重点实验室海洋与地球科学学院,上海200092;同济大学海洋地质国家重点实验室海洋与地球科学学院,上海200092【正文语种】中文因版权原因,仅展示原文概要,查看原文内容请购买。
希尔伯特黄变换及其应用1. 应用背景希尔伯特黄变换(Hilbert-Huang Transform,HHT)是一种用于非平稳和非线性信号分析的方法,由中国科学家黄钧提出。
传统的傅里叶变换等线性方法仅适用于平稳信号,而在实际应用中,许多信号都是非平稳的,因此需要一种更加灵活和准确的分析方法。
希尔伯特黄变换结合了经验模态分解(Empirical Mode Decomposition,EMD)和希尔伯特变换(Hilbert Transform),能够有效地分解非平稳信号,并提取出其局部特征。
2. 应用过程希尔伯特黄变换的应用过程主要包括以下几个步骤:2.1 数据采集与预处理首先需要采集到待分析的非平稳信号,并进行预处理。
预处理包括去除噪声、滤波等操作,以提高信号的质量和准确性。
2.2 经验模态分解(Empirical Mode Decomposition,EMD)经验模态分解是希尔伯特黄变换的核心步骤,用于将非平稳信号分解成一系列固有模态函数(Intrinsic Mode Functions,IMF)。
IMF是一组具有局部特征的函数,它们能够准确地描述信号的本质。
经验模态分解的具体步骤如下: - 将信号的极大值点和极小值点连接起来,得到信号的上包络线和下包络线; - 计算信号的局部平均值(上包络线加下包络线的平均值),得到信号的均值函数; - 用原始信号减去均值函数,得到第一次分解得到的第一固有模态函数(IMF1); - 对IMF1进行局部极值点的连接和平均值的计算,得到IMF1的上包络线和下包络线; - 用IMF1减去上包络线和下包络线的平均值,得到第二次分解得到的第二个固有模态函数(IMF2); - 重复以上步骤,直到最后得到的IMF满足一定的停止准则。
2.3 希尔伯特变换(Hilbert Transform)希尔伯特变换是一种用于计算信号的分析信号的方法,可以将实数信号转换为复数信号,并提取出信号的相位信息。
汶川地震基岩强震记录的Hilbert-Huang变换分析
裴强;胡波
【期刊名称】《煤炭学报》
【年(卷),期】2011()S2
【摘要】详细解释了Hilbert-Huang变换(HHT)的思想和基本理论,利用HHT方法对汶川8.0级的地震记录进行了经验模态分解(EMD),并对其Hilbert-Huang三维时间-频率-幅值(能量)图谱、边际能量谱和瞬时能量谱进行了分析讨论。
研究表明,基于HHT的Hilbert谱能够精准和清晰地揭示地震记录的时-频-能量分布。
【总页数】6页(P268-273)
【关键词】Hilbert-Huang变换;经验模态分解;边际能量谱
【作者】裴强;胡波
【作者单位】大连大学土木工程技术研究与开发中心
【正文语种】中文
【中图分类】TD82;TQ53
【相关文献】
1.基于汶川地震远场强震动记录的厚覆盖土层对长周期地震动影响分析 [J], 徐扬;赵晋泉;李小军;马秀芳;赵向佳
2.汶川地震泉州台强震仪与地震仪的记录分析 [J], 陈先兰
3.汶川Ms 8.0级地震强震记录所揭示的地震断层特征分析 [J], 张冬丽;徐锡伟;周正华;于贵华
4.基于基岩衰减公式的汶川·芦山地震强震动评价 [J], 司宏俊;纐缬一起;三宅弘惠;李小军
5.汶川地震在江苏的强震动记录及地震影响场 [J], 朱升初;杨伟林;唐忠良;朱祥健;孔建国
因版权原因,仅展示原文概要,查看原文内容请购买。
命、财产造成危害。
对于土木工程结构物,一旦建成投入使用后,除了材料自身性能会不断退化、老化外,还会受到风、地震、疲劳、超限、地震、人为因素等作用,从而导致结构或构件有不同程度的自然累积损伤和突发性损伤。
在科技如此发达的今天,人类仍无法完全避免爆发性的天灾人祸的袭击,往往造成结构突发性地严重损伤,甚至完全毁坏,对人类的生命、财产安全造成严重危害。
一系列自然界和人为的灾祸表明对结构物随时间变化的破坏程度作出具有工程意义的损伤指标并进行识别的重要性,尤其是对地震过程中以及震后结构破坏情况识别。
因为它不仅涉及到是否在震后对建筑物进行拆除或采取相应的加固措施,当地震荷载施加到重要性等级高的建筑时,我们更需要实时识别破坏的发生时刻,严重程度和大致位置。
=十世纪末的台湾集集大地震令中国抗震界人士记忆犹新【3】。
1999年9月21日凌晨1时47分12.6秒我国台湾省发生了震级ML7.3级的集集地震,震中位于北纬23.87度、东径120.75度的台湾南投县集集镇、车笼埔断层上,震源深约为7.0公里,主震持续102秒。
此次地震导致2412人死亡,11305人受伤;26831栋房屋全倒,24495栋房屋半倒。
图卜一1显示了震后房屋不同程度破坏情况。
蹦1—1台湾集集地震房屋破坏情况理想的损伤检测目标是对图l—l中2-#建筑在地震时的破坏过程作出实时识别,对1#建筑震后损伤情况作出评估。
对于一些担任着重要社会和政治功能的生命线工程,在设计时采取了较高的抗震等级措施。
一次地震灾害的发生一般伴随着一次或数次强震和若干次余震,建筑物可能不会在首次强震时完全破坏,但是由于损伤积累,不能保证它在后续地震时不随时倒塌。
1994年美国加州Northridge地震、1995年日本神户地震和1999年“9.21”台湾大地震中,一些0·80·70.6量0.5兰0.4△墨0·30·20.10.0MarginalSpectrum(C)FourierSpectrum5lO1520250Frequegcy(D)InstaneousEnergy510152025Frequency30405060Time(s)(E)续图3—4Elccntro分析结果65432lO枷枷啪啪啪枷姗。
希尔伯特-黄变换在地震资料处理中的应用
刘庆敏;杨午阳;李书平;何欣
【期刊名称】《地震工程学报》
【年(卷),期】2009(031)003
【摘要】希尔伯特-黄变换是最新发展起来的处理非线性非平稳信号的时频分析方法.本文介绍了该方法的基本原理和几个关键技术,以及其优势所在,并将它应用于含随机噪音的地震信号和实际地震资料的分解.分解后的信号能够重构原始信号,重构后在信噪比和分辨率方面有了一定的提高,对地震资料的去噪有重要意义.
【总页数】6页(P211-216)
【作者】刘庆敏;杨午阳;李书平;何欣
【作者单位】中国石油勘探开发研究院西北分院,甘肃,兰州,730020;中国石油勘探开发研究院西北分院,甘肃,兰州,730020;中国石油玉门油田分公司勘探开发研究院,甘肃,酒泉,735019;中国石油勘探开发研究院西北分院,甘肃,兰州,730020
【正文语种】中文
【中图分类】P681.4
【相关文献】
1.基于希尔伯特变换的信号解调算法及其在飞机供电特性参数测试系统中的应用[J], 李小舟; 金海彬
2.空间域加窗二维希尔伯特变换在三维地震资料体边缘检测中的应用 [J], 吕丙南;陈学华;徐赫;刘芸菲;罗鑫;周晨
3.改进希尔伯特-黄变换方法在钢轨波磨检测中的应用 [J], 吴传奇;柴晓冬;李立明;
郑树彬
4.希尔伯特—黄变换在地震资料去噪中的应用 [J], 王强;刘金辉;叶恒
5.基于CEEMDAN的希尔伯特变换在海底天然气水合物地震探测中的应用 [J], 夏秋萍;刘怀山
因版权原因,仅展示原文概要,查看原文内容请购买。
第34卷第2期福州大学学报(自然科学版)Vol.34No.2 2006年4月Journal of Fuzhou University(Natural Science)Apr.2006文章编号:1000-2243(2006)02-0260-05希尔伯特-黄变换谱及其在地震信号分析中的应用陈子雄,吴琛,周瑞忠(福州大学土木建筑工程学院,福建福州350002)摘要:介绍了希尔伯特-黄变换(HHT)这一非线性、非平稳信号处理方法,并利用HHT处理了地震工程中常用的El Centro地震波,得到了该信号的Hilbert谱、边际谱和能量谱,提取了该信号的主要动力特性,并与该信号的Fourier分析结果进行了对比,显示出HHT这一方法的优越性.关键词:希尔伯特-黄变换;经验模态分解;固有模态函数;地震信号中图分类号:TU311.3文献标识码:AHilbert-Huang transform spectru m and its application in seismic signal analysisCHEN Zi-xiong,W U Chen,ZHOU Rui-zhong(College of Civil Engineering and Architecture,Fuzhou University,Fuzhou,Fujian350002,China) Abstract:HHT is a ne w method to deal with non-linear and non-stationary data.El Centro earth-quake wave is analyzed by HHT.Through the way,Hilbert spectrum,marginal spectrum and energyspec trum are got and dynamic property is extrac ted.The comparison between HHT spectrum and Fourierspec trum is made and the superiority of HHT is demonstrated.Keyw ords:Hilbert-Huang transform;empirical mode decomposition;intrinsic mode function;seismic signal地震信号具有短时、突变等特点,是一种典型的非平稳随机信号,必须对其进行分析与处理,才可以提取信号的主要特征.传统的Fourie r变换能够表述信号的频率特性,但不提供任何时域信息[1],而小波分析虽然在时域和频域都具有很好的局部化性质,但本质上仍是一种窗口可调的Fourier变换,在小波窗内的信号必须是平稳的,因而没有根本摆脱Fourier分析的局限[2].小波基的选择也是信号分析中的一个重要问题,另外,小波基的有限长会造成信号能量的泄漏,使信号的能量-频率-时间分布很难定量表述.Hilbert-Huang变换(HH T)的信号处理方法被认为是近年来对以Fourier变换为基础对线性和稳态谱分析的一个重大突破[2].它由经验模态分解(E mpirical Mode Decomposition,E MD)方法和Hilbert变换(H T)两部分组成,其核心是E MD分解.该方法采用了固有模态函数(Intrinsic Mode Function,I MF)概念以及将任意信号分解为I MF组成的思想,即E MD法,使得瞬时频率具有实际的物理意义[3].它不受Fourier分析的局限,可依据数据本身的时间尺度特征进行模态分解,分解过程中保留了数据本身的特性,再对各I MF分量进行Hilbert变换,得到信号能量在时间尺度上的分布规律,实现地震动力特性的提取.1Hilbert-Huang变换1.1经验模态分解和固有模态函数经验模态分解(EMD)的目的是通过对非线性非平稳信号的分解获得一系列表征信号特征时间尺度的固有模态函数(I MF),使得各个I MF是窄带信号,可以进行Hilbert分析.首先设定两个条件:¹整个时间序列的极大极小值数目与过零点数目相等或最多相差一个;º时间序列的任意点上,由极大值确收稿日期:2005-07-27作者简介:陈子雄(1981-),男,硕士研究生;通讯联系人:周瑞忠,教授.基金项目:教育部博士点专项科研基金资助项目(20040386004)定的包络与由极小值确定的包络的均值始终为零.能满足以上两个条件的信号称为I MF 信号,用Hilbert 变换求I MF 信号的瞬时参数,其结果是准确的.求取I MF 的具体步骤如下:1)找出原始序列x (t)的各个局部极大值.这里,为更好地保留原序列的特性,局部极大值定义为时间序列中的某个时刻的值,其前一时刻的值不比它大,后一时刻的值也不比它大.然后用三阶样条函数进行插值,得到原序列的上包络序列值x max (t).同理,可以得到下包络序列值x min (t ).2)对每个时刻的x ma x (t)和x min (t)取平均,得到瞬时平均值m(t):m (t)=[x max (t)+x min (t)]/2(1) 3)用原序列x (t)减去瞬时平均值m(t),得到一个去掉低频的新序列h(t ):h(t)=x (t)-m(t)(2)此时,如果h(t)满足I MF 的两个条件,那它就是固有模态函数.否则,把h (t)当作原序列,重复以上步骤,直至满足I MF 的定义,得到第一个固有模态函数,并以c 1(t)记之.一般来说,c 1(t)代表了原始序列中的高频部分.然后,用原序列减去c 1(t),得到剩余值序列r 1(t):r 1(t)=x (t)-c 1(t)(3)至此,提取第1个固有模态函数的过程全部完成.然后,把r 1(t)作为一个新的原序列,按照以上步骤,依次提取第2,第3,,,,直至第n 个固有模态函数c n (t).之后,r n (t)变成一个单调序列,再也没有内在模函数能被提取出来.如果把分解后的各分量合并起来,就得到原序列x (t):x (t )=E nj =1c j (t)+r n (t )(4)1.2 Hilbert 变换和Hilbert 谱通过E MD 分解得到的I MF 分量非常适合作Hilbert 变换,进而求出瞬时频率,得到HHT 谱.Hilbert 变换是一种线性变换,如果输入信号是平稳的,那么输出信号也应该是平稳的;Hilbert 变换强调局部属性,这避免了Fourier 变换时为拟合原序列而产生的许多多余的、事实上并不存在的高、低频成分.对于任一固有模态函数c(t),其Hilbert 变换^c (t)定义为[4]:^c (t )=1P P V Qc(S )t -Sd S (5)式中:P V 代表柯西主值,则对于c(t)的解析信号z (t)为:z (t)=c(t)+i^c (t )=a(t )ei H (t )(6)式中:a(t)和H (t)分别为信号x (t)的瞬时振幅和瞬时相位,按下式计算:a (t)=c 2(t)+^c 2(t)(7)H (t)=arctan (^c (t)/c(t))(8)由瞬时相位可得到信号的瞬时频率:X (t)=d H (t)/d t (9)可见,由Hilbert 变换得到的振幅和频率都是时间的函数,如果把振幅显示在频率-时间平面上,就可以得到Hilbert 幅值谱H (X ,t),简称Hilbert 谱,记作:H (X ,t)=Re E nj=1a j (t)e i Q X j(t)d t(10)H (X ,t)精确地描述了信号的幅值随时间和频率的变化规律.将H (X ,t)对时间积分,就得到Hilbert边际谱h (X ):h(X )=Q TH (X ,t)d t (11)边际谱提供了对每个频率的振幅量测,表达了在整个时间长度内振幅的累积.将振幅的平方对频率积#261#第2期陈子雄,等:希尔伯特-黄变换谱及其在地震信号分析中的应用分,可定义Hilbert 瞬时能量I E (t),如下式所示:I E (t)=QXH 2(X ,t)d X (12)Hilbert 瞬时能量提供了信号能量随时间的变化情况.将振幅的平方对时间积分,可得到Hilbert 能量谱E S (X ),如下式所示:E S (X )=Q TH 2(X ,t)d t (13)Hilbert 能量谱提供了对于每个频率的能量量测,表达了每个频率在整个时间长度内所累积的能量.2 EMD 前的准备工作2.1 端点效应的抑制在运用EMD 方法对非线性资料进行分解时,必须进行端点抑制.否则,要么会因在端点处弃值而严重影响资料的完整性;要么会因在端点处的发散而使运算溢出;如果抑制的不好,又会因/污染0度过大而使分解严重失真.本文采用了文献[5]介绍的在端点处依据端点附近数据变换插入样条插值起点的方法.在每一次I MF 分解前,先判断待分解系列起始段的变化趋势,若为上升,则将第一个局部极大值作为上包络样条插值的起点,而直接把端点作为下包络样条插值的起点.若为下降,直接把端点作为上包络样条插值的起点,将第一个局部极小值作为下包络样条插值的起点,信号另一端按同样方法处理.该法虽具有一定强制性,也缺乏明确的理论依据,但确能有效抑制端点效应.2.2 确定固有模态函数的判据要分解出固有模态函数,必须确定I MF 的判据,否则,无限重复的结果只可能得到定常振幅的调频波,这就失去了应有的物理意义.习惯上,用前后h(t)的标准差S 的大小作为I MF 的判据:S =ETt=0h k-1(t)-h k (t)2ETt=0h k-1(t)2(14)一般说来,S 的值越小,所得的固有模态函数的线性和稳定性就越好,能够分解出的固有模态函数的个数就越多.实践表明,当S 的值介于0.2~0.3时,既能保证固有模态函数的线性和稳定性,又能使所得的固有模态函数具有相应的物理意义[6].3 HHT 在地震信号分析中的应用作为Hilbert-Huang 变换在地震信号分析中的一个简单实例,对地震工程中最常用的一条地震记录El Centro 地震波(图1)的南北向分量进行分析.El Centro 波经E MD 分解,形成8个I MF 分量和一个残量,如图2所示.将分解的I FM 分量实行Hilbert 变换,得到地震波的时间-频率-振幅分布图(图3),#262#福州大学学报(自然科学版)第34卷利用公式(11)、(12)、(13)分别得到边际谱(图4)、瞬时能量谱(图5)、Hilbert 能量谱(图6),为比较还作了地震信号的Fourier 频谱(图7)、功率谱密度(图8).表1、表2分别列出了信号能量集中区的一些能量值.表1 信号能量集中区各时段所对应的能量Tab.1 Energy of the every period of centralized district of signal energyt /s 0~11~22~33~44~55~66~77~88~99~10能量412802990900325290074387012072006382807491294053482610584550t /s 10~1111~1212~1313~1414~1515~1616~1717~1818~1919~20能量22562013448001919202142802321901368201524001444108677299118t /s 20~2121~2222~2323~2424~2525~2626~2727~2828~2929~30能量174330954289803460215257140460500294810551894921530926#263#第2期陈子雄,等:希尔伯特-黄变换谱及其在地震信号分析中的应用表2 信号能量集中区各频段所对应的能量Tab.2 Energy of the every frequency band of centralized district of signal energyf /Hz 0~11~22~33~44~55~66~77~88~99~10能量78036004567400121920040473027728025035011326015446023282368814 分析与讨论1)由图2可知,E MD 首先将地震波中的高频分量分离出来,然后,所分解的IMF 分量的频率依次降低,这说明EMD 分解的效率是非常高的.同时,由于它不依赖如Fourier 、小波那样所需的基函数,自适应性好,因此更为有效地反映信号的内部特征,避免了信号能量的扩散与泄漏.2)由图2可见,E MD 方法清晰地说明了地震波中无周期性的正弦分量,如果有的话,必然满足I MF 条件分离出来,事实上并没有.这一事实证明:Fourier 分析对地震波这种非平稳信号的分析是不适用的,对非平稳信号的解释是不合理的.3)从所分解出来的各I MF 分量可以看到,前几个IMF 分量体现了原始信号的主要特征,该频率是地震信号的主要贡献频率.4)由图3可见,Hilbert 谱是随着频率分布而瞬时变化的,图中显示了该信号的大部分能量集中在10Hz 以下,30s 以内,这在边际谱、瞬时能量谱和Hilbert 能量谱中也得到体现.这在以往的信号处理方法中是不能得到的.5)由图7的Fouier 频谱和图4的边际谱对比可见,Fourier 谱给出的结果往往会低估低频部分的幅值,而高频部分又往往高于实际的幅值,这应该引起地震工程界的重视.6)图5所示的瞬时能量谱表达了信号的能量随时间的变化过程,其峰值高达155860,而在大部分能量集中的前30s 内的波能平均值只有9676,比值超过16倍,如此巨大的瞬时峰值能量势必对结构安全产生不利影响.由图6也可见,信号的能量基本上集中在10Hz 以下,这与信号的功率谱的表现情况基本一致.7)对图5所示的瞬时能量谱每时刻能量求和,可得其结果为1.4787@107;同理,对图6所示的Hilbert 能量谱的求和,其结果仍然相同,相对误差不到十万分之一,其中的误差是由于计算截断引起的,这也证实了Paraseval 公式的正确性.参考文献:[1]大崎顺彦(日).地震动的谱分析入门[M].吕敏申,谢礼立,译.北京:地震出版社,1980.[2]Huang N E,Shen Z,Long S R,et al.The empirical mode decomposi tion and Hilbert spectrum for nonlinear and non-stationaryti me series analysis[J].Proc R Soc Lond,1998,A454:903-995.[3]钟佑明,秦树人,汤宝平.一种振动信号新变换法的研究[J].振动工程学报,2002,15(2):233-238.[4]郑君里,应启珩,杨为理.信号与系统[M].北京:高等教育出版社,2000.[5]杜寿昌.基于Hilbert-Huang 变换的心音信号时频分析研究[D].昆明:云南大学,2003.[6]熊学军,郭炳火,胡筱敏,等.E MD 方法和Hilbert 谱分析法的应用与探讨[J].黄渤海海洋,2002,20(2):12-21.#264#福州大学学报(自然科学版)第34卷。