EMD新技术在数字波形预处理中的初步应用.pdfVIP

第26卷第1期 华 南 地 震 V01．26，No．1 2006年3月 Mar．，2006 SOUTHCHINA OFSELSMOLOGY JOURNAL EMD新技术在数字波形预处理中的初步应用 武安绪，林向东，穆会泳，赵桂儒，柴金翼，吴培稚 (北京市地震局，北京100080) 摘要：介绍了EMD技术，给出由该技术产生的固有模态函数重构数字地震波形信号处理中 所用滤波器算法。然后利用该方法对实际波形信号进行了初步处理，处理结果表明，基于 EMD的滤波技术具有许多其它分析手段所不具备的特点，是一种新的滤波方法，可用于数 字地震波形信号的预处理。 关键词：滤波；数字地震波形；EMD技术 中图分类号：P315．61文献标识码：A 日U吾 由于数字化设备采样率与灵敏度的大幅度提高，在数字地震波形记录中，干扰因素也 随之增多。在一些数字地震波形分析与研究工作中，如初动检测、震相识别、Q值计算等， 都需要对波形进行一定的滤波、分频或降噪等预处理【”，消除噪声或滤去不需要的频率成 分，分离干扰，突出信号，压制噪声，提高信噪比。可采用的方法较多，如滤波、基线校 正、灭尖、光滑、拟合等【2】，但一般线性处理方法不能很好地排除这些非规则干扰，只能 消除某单一性质的干扰，且可能出现相位偏移现象，为此一些学者先后引入多道维纳滤波、 卡尔曼滤波[31、高斯滤波【21、小波分析【4咽等非线性资料处理方法，取得了一定效果。本 Mode 文则基于经验模态分解(Empirical 波方法，进行数字地震波形初步处理试验。 EMD技术是一种新的非线性多尺度分解处理方法，适用于处理非稳态、高度非线性的 动态变化数据。该技术首先把一时间序列分解成一组不同尺度的内在的客观的固有模态函 Mode 数(Intrinsic 成高通、低通、带通(包括多通带，可以不相邻)自适应滤波器。本文基于该技术，通过对 数字波形的分解，得到各阶IMF函数，进行低通、高通、带通为目标的组合与重构，实现 数字地震波形的滤波、频带分离和降噪。 1 EMD技术简介与IMF滤波器构成原理 收稿日期：2005—12—29 基金项目：上海地区震源参数及介质研究(上海市财政专项030204)；北京市tl然科学基金项目(8041001) 作者简介：武安绪(1967一)，男，副研究员．主要从事地震预报研究和软件研制。 华南地震 26卷 1．1 EMD技术简介[6--s] EMD技术的本质是通过特征时间尺度获得本征震荡模式，然后由本征震荡模式来分解 时间序列资料。设地震波形信号时间序列为X(力，则： (1)找出X(t)的所有极大值点和极小值点，将其用三次样条函数分别拟合为原数据序 列的上、下包络线。上、下包络线的均值为平均包络线m，。，将原序列减去m，。便可得到一 个去掉低频的新序列hm即 hlo=XO)-mlo (1) 一般h，。不一定是一个平稳序列，为此需要对它重复上述过程。如果h，。的平均包络线 为mn，则去除该包络线所代表的低频成份后的序列为矗n，即 危11=^10一，n11 (2) 重复上述过程，经k次循环后，使得到的平均包络m。t趋向于零，此时的h。t为第一阶 IMF序列，定义为分量C，，它表示在信号数据序列XO)中最高频的成份。 解，便得到第二阶IMF分量C：。如此重复直到最后一个序列r。不可再被分解时为止。这时 的k代表序列X(t)的残余项，即通常为X0)的趋势项或均值。 上述过程可以表述为 rl=x