[发明专利]一种微震信号降噪方法

权利要求书

1.一种微震信号降噪方法，采用NLMS自适应滤波器，其特征在于，包括如下步骤：

步骤1：设微震信号的时序序列为x(t)，对微震信号的时序序列x(t)进行NLMS自适应滤波预处理，得到信号x'(t)，其中，t＝1，2，…，n，n为微震信号的采样点个数；

步骤2：给信号x'(t)叠加一组高斯白噪声信号m(t)，得到一个总体信号X(t)，即X(t)＝x′(t)+m(t)；

步骤3：根据公式(1)对总体信号X(t)进行EEMD分解，得到m阶IMF分量；

$X\left(t\right)=\underset{j=1}{\overset{n}{\Σ}}{c}_j+r_n---\left(1\right);$

其中，c_j表示第j个IMF分量；r_n为分解残余项；

步骤4：根据公式(3)对步骤3中的m阶IMF分量进行阈值处理；

$\stackrel{\‾}{{\mathrm{imf}}_k}\left(t\right)=\left\{\begin{array}{cc}{\mathrm{imf}}_k\left(t\right),& |{\mathrm{imf}}_k\left(t\right)\≥T_k|\\ 0,& |{\mathrm{imf}}_k\left(t\right)\<T_k|\end{array}\right.,(k=2,...7)---\left(3\right);$

其中，表示经阈值处理降噪后的IMF；T_k为通用阈值，其表达式如公式(4)所示；

$T_k=C{\widehat{\mathrm{\σ}}}_i\sqrt{(2\×\ln (n\left)\right)}---\left(4\right);$

其中，C为阈值系数；为第i个IMF分量的估计能量；

步骤5：重构出降噪后的微震信号y'(t)，即