[发明专利]地球物理信号去噪方法

说明书

技术领域

本发明涉及工程地球物理勘探技术领域，具体地指一种地球物理信号去噪方法。

背景技术

工程地球物理勘探，简称工程物探，它的探测目标是地下岩土层或者建筑结构，对探测目标加载自然或人工物理场后，通过仪器观测其变化，以确定隐蔽目标的空间范围(大小、形状、分布)，并测定目标体的物理参数，达到解决地质问题的一种物理勘探方法。现阶段工程物探的技术手段主要包括弹性波方法和电磁波方法两大类。

由于工程物探中待勘探目标通常埋深较浅、缺陷较小，所以不管选用弹性波还是电磁波，工程物探采用的发射源能量较小，频率较高，这也导致了工程物探所使用更容易受外界各种因素干扰。所以，在工程物探信号处理中，对采集到的信号中的噪音进行压制显得尤为重要。但现有的各类去噪方法不能很好的满足工程应用的需要。

发明内容

本发明的目的就是要提供一种地球物理信号去噪方法，该方法能简便高效的去除地球物理信号中的噪声。

为实现此目的，本发明所设计的一种地球物理信号去噪方法，其特征在于，它包括如下步骤：

步骤1：对原始地球物理信号进行stockwell变换，得到原始地球物理信号的时频谱；

步骤2：将步骤1所得的原始地球物理信号的时频谱进行时频滤波；

步骤3：将步骤2的时频滤波后信号进行阈值滤波；

步骤4：将步骤3所得的阈值滤波的结果进行stockwell变换逆变换，得到初步去噪信号；

步骤5：将步骤4所得的初步去噪信号进行经验模态分解，得到初步去噪信号的若干个本征模态函数分量；

步骤6：将步骤5所得的初步去噪信号的所有本征模态函数分量进行希尔伯特变换，对各个本征模态函数的希尔伯特谱进行分析，分析方法为将各个本征模态函数的希尔伯特谱中的畸变点分别与被检测对象的可能存在的畸变点进行对比，留下畸变点与可能存在的畸变点一致的本征模态函数；

步骤7：将步骤6留下的本征模态函数相加，即可实现信号重构，得到所需的去噪信号。

本发明针对不同的观测系统和相应的工程物探采集信号，在步骤2、步骤3中选取合适的时频滤波和阈值滤波参数，在步骤6中选择相应的本征模态函数，可以获得较好的信号去噪效果，有利于下一步工作的进行，使工程物探探测结果更加准确、更加符合实际情况。

附图说明

图1为本发明的信号去噪方法流程图

图2为采集信号波形图；

图3为信号stockwell变换后的时频谱；

图4为经过时频滤波和阈值滤波后得到初步去噪信号的波形图；

图5为初步去噪信号经验模态分解的结果；

图6为重构后去噪信号的波形图。

具体实施方式

以下结合附图和具体实施例对本发明作进一步的详细说明：

如图1所示的地球物理信号去噪方法，它包括如下步骤：

步骤1：对原始地球物理信号(地震信号、锚杆检测信号、基桩检测信号或雷达信号，这里以地震信号为例，如图2所示，其横坐标为时间，纵坐标为传感器位移量)进行stockwell变换，得到原始地球物理信号的时频谱(如图3所示，某地震采集信号S变换谱，其横坐标表示时间，纵坐标表示频率，亮度表示该点的能量大小)；

步骤2：将步骤1所得的原始地球物理信号的时频谱进行时频滤波；

步骤3：将步骤2的时频滤波后信号进行阈值滤波；

步骤4：将步骤3所得的阈值滤波的结果进行stockwell变换逆变换，得到初步去噪信号(如图4所示，某地震采集信号初步去噪后的结果，其横坐标表示时间，纵坐标表示传感器位移量)；

步骤5：将步骤4所得的初步去噪信号进行经验模态分解(Empirical Mode Decomposition,EMD)，得到初步去噪信号的若干个本征模态函数分量(如图5所示，某地震采集信号初步去噪后经验模态分解的结果，其横坐标表示时间，纵坐标表示传感器位移量，第一条曲线表示初步去噪的信号，第二条曲线表示模态1，第三条曲线表示模态2，由此类推)；

步骤6：将步骤5所得的初步去噪信号的所有本征模态函数分量进行希尔伯特变换，对各个本征模态函数的希尔伯特谱进行分析，分析方法为将各个本征模态函数的希尔伯特谱中的畸变点分别与被检测对象(原始地球物理信号是对被检测对象的检测信号)的可能存在的畸变点进行对比，留下畸变点与可能存在的畸变点一致的本征模态函数；