基于改进的EEMD实现地震信号去噪.pdfVIP

刘伟等：基于改进的EEMD 实现地震信号去噪 对于信号x(t )(t 0,1, ,N ) ，定义 为循环平 p m m S n b ik r x(t)tr (8)    移n 次的平移算子，即 k 0 pr t m r m S (x(t)) x(t n)mod N (1) 因此，给定需要拟合的点数m ，多项式阶数p ， n 则相应的逆平移算子为 求解线性方程组(8) ，就可以求出系数序列 (S )1 S (2) n n (b ,b ,b b ) ，从而多项式F (t ) 就可以确定 p 0 p 1 p 2 p p 通过平移消除伪吉布斯现象的过程可表示为 x S (T (S (x))) (3) 了。 n n 式中： 为分析算子； 为原始信号 基于平移 T x x (t ) 2.4 改进的EEMD 去噪 不变阈值去噪后的结果。 改进的EEMD 具体算法步骤为 基于平移不变小波阈值去噪的步骤如下： （1）对待分析信号作EEMD 分解，得到n 个 （1）确定平移范围，对含噪信号作循环平移； IMF 分量； （2 ）对每次循环平移得到的信号进行小波阈 （2 ）对IMF 分量进行分类，对高频IMF 分量 值去噪处理； 作平移不变小波阈值处理，对 IMF 低频分量作 （3 ）对去噪处理后的信号作相应的逆循环平 Savitzky-Golay 滤波处理； 移，使其与原始信号具有相同的相位； （3 ）利用处理后的IMF 分量进行信号重构得 （4 ）对所得到的结果求平均值，即可得到去 到去噪后的结果。 噪后的结果。 3 仿真试验 2.3 Savitzky-Golay 滤波 Savitzky-Golay 滤波的基本原理为：设有一组 现采用同一个含噪合成地震记录分别使用传 数据为x (t ) ，然后选取该数据上任意一点附近固定 统EEMD 和改进EEMD 进行去噪处理，结果如图1 个数的点拟合一个多项式，该多项式在任意一点的 所示。由图可见，传统EEMD 的去噪结果中还残留 值就是平滑处理后的数据。取某点附近2m 1个连 部分高频噪声，而改进EEMD 较好地压制了随机噪 续的点进行p 阶多项式(p 2m 1) 拟合，记该多项 声，且与原始信号的吻合度更高，去噪效果更为理 式为 想。表1 显示了两种方法去噪前、后信噪比和标准 p