基于经验模态分解的地震前兆资料处理与分析.pdfVIP

基于经验模态分解的地震前兆资料处理与分析 武安绪。 吴培稚” 宋治平2’ 1)北京市地震局北京 1000802)上海市地震局上海200062 地震前兆观测中干扰因素较多，如果观测的目标是固体潮汐，则温度、气压、降水、地下水、地震波等信号 都是干扰。前兆数字化设备大量投入观测之后，由于采样率与灵敏度的大幅度提高，干扰因素也会随之增多。在 数据处理中，排除这些干扰因素的传统方法有光滑、拟合、滤波、差分等，期望达到分离干扰、压制噪声、提取 异常等目的。在实际应用中，为了排除非线性动态变化的干扰信息，一些学者先后引入了可处理非线性资料的多 道维纳滤波、卡尔曼滤波、小波分析等一系列方法，取得了一定效果。 et Mode 近年出现的经验模态分解(Huangal，1998)(Empirical Mode 列信号分解为一组不同尺度的内在的、客观的固有模态函数(IntrinsicFunction，IMF)。EMD是一种新的非 线性信号分解方法，适用于非稳态、非线性的动态数据。目前，EMD分解已用于多项研究。 在EMD分解基础上，把IMF进行适当组合，便可构成低通、带通与高通自适应滤波器。本文利用EMD处理地 震前兆资料。首先利用前兆数字化观测资料，探索消除前兆资料中的干扰成份，如对同震信号剥离、单点突跳消 除、高频干扰分离、低频或超低频背景趋势消除。最后对前兆观测数据的异常信息提取作了初步尝试。 经以上详细计算与分析认为，EMD分解与IMF组合滤波方法在处理地震前兆资料时具有以下主要特点： (1)EMD分解中越是早分解出来的IMF频率越高，第一个分解出来的IMF函数代表原信号的最高频率成份，分 解到最低频为止(为一单调函数或很小的量)。 (2)每阶IMF为一调频调幅稳态信号，不但是希尔伯特变换的前提，也是其它稳态方法进一步分析与处理的基础。 (3)在刻画非线性与动态性资料方面，EMD方法与小波多尺度分解相比较丽言具有较大的优势与潜在的应用价 值(朱光明等，1993)。EMD分解过程是内在的、客观的和自适应的，反应序列的内在固有特征，因此具有较高的 分解、组合滤波精度与频率分辨率。每阶IMF序列与原数据闾的有机组合可取得不同的滤波效果，满足不同的后 续处理与分析研究的需要。 最后通过综合分析，得以下结论： (1)本文介绍TEMD分解方法的原理及其相应算法，给出了已知信号的仿真处理结果，然后把基于EMD的滤波 技术引入到前兆观测时间序列信号处理领域之中； (2)EMD方法与小波多尺度分解具有一定类似之处，但由于EMD方法没有采用特定的核函数，分解的思想不一 etal，1998)。 样，具有较高的分解性能。EMD分解出来的IMF信息不会重叠，不存在信息沉余的问题(Huang (3)对前兆观测资料中存在的典型干扰实例进行EMD分解、IMF组合滤波的分析与处理，初步展示了基于EMD 分解的地震前兆信号滤波方法与技术思路，通过滤波前后信号的实际处理与应用效果比较，认为取得了～定成效； (4)EMD分解获得的IMF序列本身就是多通带滤波的结果，再包括一定的IMF组合，可用来展示前兆异常信息 的频带分布与时回进程演化特征，揭示异常特征、信息量、时间因子与频率之间的特定关系，对地震前兆信息的 解释、异常特征的分析与提取带来了方便。通过对IMF序列及其组合滤波的初步研究，认为EMD方法可作为地震 前兆分析中异常信息提取的一种得力工具。 以上研究结果表明，EMD分解与IMF组合滤波方法将是一重要的数字地震前兆信号处理方法与异常信息提取 的得力工具。 本研究由北京市自然科学基金项目(8041001)与上海市财政专项(030204)资助。 参考文献 朱光明，高静怀，王玉贵．1993．小波变换及其在一维滤波中的应用IJl．石油物探，32． N R a1．1998．The mode andtheHilbert for E。SheenZ，StevenL，et empiricaldecomposition spectrum Huang