一种基于变分模态分解的瞬变电磁探测信号的降噪方法技术

本发明专利技术涉及一种基于变分模态分解的瞬变电磁探测信号的降噪算法，将采空区强干扰瞬变电磁信号的分解过程转换至变分问题的求解。通过对强干扰瞬变电磁进行变分构造，引入二次惩罚因子和拉格朗日惩罚算子将约束性变分问题转化为非约束性问题,通过乘法算子交替方向算法获得有限个模态分量，实现信号的有效分离，将噪声和有效信号分离开来。通过本发明专利技术，能够自适应地对信号进行分解，不需要先验知识，实际操作更为简单方便；变分模态分解算法算法简单，计算机运行用时短，节省储存空间和计算时间；有效解决了经验模态分解中的模态混叠问题，提取的信号更加精确，表现出更好的噪声鲁棒性。

【技术实现步骤摘要】
一种基于变分模态分解的瞬变电磁探测信号的降噪方法
本专利技术涉及信号处理
，具体是一种基于变分模态分解的瞬变电磁探测信号的降噪方法。
技术介绍
瞬变电磁法作为一种新兴的地球物理勘探方法，因随机干扰小、对低阻体反应灵敏、穿透高阻能力强、勘探深度大和施工方便等优点被世界各国广泛应用。国内主要利用瞬变电磁方法精确探测含水采空区，防止突水事故的发生，保证煤矿工作的安全运行。然而，瞬变电磁探测在现场实施中总会受到周围噪声的干扰，导致TEM信号发生畸变或失真，影响数据的处理和解释，从而得到与实际不符的结论，降低了物探解释的准确性，甚至可能会导致事故发生，因此很有必要对瞬变电磁信号进行降噪处理。然而，目前针对瞬变电磁信号频带范围宽、晚期信号弱、非平稳性、非线性的特点提出的主要降噪方法都存在着一些自身的缺陷。小波变换降噪根据噪声与信号的小波分解系数的不同，去除各频带上的噪声对应的小波系数，保留原始信号的小波分解系数，然后对处理后的系数进行小波重构，得到纯净信号，相比起传统傅里叶变换的纯频域分析，小波变换能同时在时域和频域中对信号进行分析，有效区分信号中的突变部分和噪声，但在降噪时需预先选定小波基和分解层数，缺乏自适应性。经验模态分解方法把瞬变电磁信号分解成若干个从高频到低频的本征模态函数，再对每一个本征模态分量进行频谱分析，根据条件筛选、组合，得到重构结果，这种方法具有自适应性、不需要先验信息、多分辨、多分解特性，十分适于非线性非平稳信号的分解，然而分解过程中存在的模态混叠和端点效应正成为经验模态分解方法亟待解决的问题。数字滤波将瞬变电磁信号及其噪声看成随机信号，利用其统计算法和傅里叶变换估计出信号本身，从而提高信噪比，虽然可靠性高，但需要噪声和有效信号的先验知识，不适应非平稳过程。独立成分分析降噪通过将接收到的瞬变电磁信号进行某种线性分解，使其分离为独立统计成分，去除噪声基向量，保留有效信号，达到分离源信号的目的，具有自适应性好，收敛速度快的优点，但由于其算法的限制导致分解的幅值和次序具有不确定性。卡尔曼滤波将瞬变电磁信号及其噪声看成随机信号，利用其统计算法和傅里叶变换估计出信号本身，从而提高信噪比，虽然可靠性高，但需要噪声和有效信号的先验知识，数值稳定性差，计算效率低。
技术实现思路
本专利技术要解决的技术问题在于，针对现有技术的上述缺陷，提供一种基于变分模态分解的瞬变电磁探测信号的降噪方法。本专利技术解决其技术问题所采用的技术方案是：构造一种基于变分模态分解的瞬变电磁探测信号的降噪方法，包括步骤：获取采空区强干扰瞬变电磁信号，预设瞬变电磁信号包括K个单分量信号，对强干扰瞬变电磁信号进行变分模态分解，得到强干扰瞬变电磁信号的变分约束函数；将拉格朗日惩罚算子和二次惩罚因子引入强干扰瞬变电磁信号的变分约束函数，得到增广拉格朗日函数，使变分约束函数不受约束；采用乘法算子交替方法(ADMM)，按预设方式交替迭代更新增广拉格朗日函数的参数，求解增广拉格朗日函数的鞍点，直至满足预设停止条件后停止迭代更新；其中，增广拉格朗日函数的参数至少包括瞬变电磁信号单分量信号、中心角频率和拉格朗日惩罚算子；提取迭代至增广拉格朗日函数满足预设停止条件时的全部单分量信号，剔除其中的噪声信号，将剩余单分量信号重组，得到去除噪声后的瞬变电磁信号。其中，在对强干扰瞬变电磁信号进行变分模态分解，得到强干扰瞬变电磁信号的变分约束函数的步骤中，设瞬变电磁信号为x(t)，由K个单分量信号ui(t)组成，利用希尔伯特变换计算每个单分量信号的相关分析信号yi(t)，得到单分量信号的单侧频谱：其中，δ(t)为狄克拉冲击函数，*为卷积符号；根据调制和调解原理将相关分析信号yi(t)与指数混合进行修正，把ui(t)频率中心调谐到零频位置，得到相应的解调信号y′i(t)：其中，ωi(t)为单分量信号ui(t)的中心角频率；通过解调信号y′i(t)进行高斯平滑度计算，得到的约束变分函数为：其中，{uk}＝{u1,u2,…,uk}为K个单分量信号ui(t)的集合，{ωk}＝{ω1,ω2,…,ωk}为K个单分量信号ui(t)中心角频率的集合，变分函数约束条件是K个单分量信号之和等于瞬变电磁信号。其中，增广拉格朗日函数的公式表示为：其中，λ为拉格朗日惩罚算子，α为二次惩罚因子；利用拉格朗日惩罚算子λ和二次惩罚因子α使变分函数不受约束，将求解变分函数的解变为求解增广拉格朗日函数的鞍点。其中，在采用乘法算子交替方法通过交替更新单个分量信号ui(t)、中心角频率ωi(t)和拉格朗日惩罚算子λ求解增广拉格朗日函数鞍点时，uk根据式(5)进行迭代更新，ωk根据式(6)进行迭代更新，λ根据式(7)进行迭代更新，直到迭代满足停止条件(8)；公式(5)-(8)的表示形式如下：其中，当对分解结果保真度要求低时，将式(7)中的τ值置为0。区别于现有技术，本专利技术提出一种基于变分模态分解的瞬变电磁探测信号的降噪算法，将采空区强干扰瞬变电磁信号的分解过程转换至变分问题的求解。通过对强干扰瞬变电磁进行变分构造，引入二次惩罚因子和拉格朗日惩罚算子将约束性变分问题转化为非约束性问题,通过乘法算子交替方向算法获得有限个模态分量，实现信号的有效分离，将噪声和有效信号分离开来。通过本专利技术，能够自适应地对信号进行分解，不需要先验知识，实际操作更为简单方便；变分模态分解算法算法简单，计算机运行用时短，节省储存空间和计算时间；有效解决了经验模态分解中的模态混叠问题，提取的信号更加精确，表现出更好的噪声鲁棒性。附图说明下面将结合附图及实施例对本专利技术作进一步说明，附图中：图1是本专利技术提供的一种基于变分模态分解的瞬变电磁探测信号的降噪方法的流程示意图。图2是本专利技术提供的一种基于变分模态分解的瞬变电磁探测信号的降噪方法的流程示意图。图3是本专利技术提供的一种基于变分模态分解的瞬变电磁探测信号的降噪方法中采用采用乘法算子交替方法更新增广拉格朗日函数的流程示意图。图4是本专利技术提供的一种基于变分模态分解的瞬变电磁探测信号的降噪方法中瞬变电磁探测信号含随机噪声干扰曲线图。图5是本专利技术提供的一种基于变分模态分解的瞬变电磁探测信号的降噪方法中瞬变电磁探测信号的理想曲线图。图6是本专利技术提供的一种基于变分模态分解的瞬变电磁探测信号的降噪方法中采用变分模态分解方法(VMD)对含随机噪声干扰的瞬变电磁信号的降噪结果图。图7是本专利技术提供的一种基于变分模态分解的瞬变电磁探测信号的降噪方法中采用经验模态分解(EMD)对含随机噪声干扰的瞬变电磁信号的降噪结果图。图8是本专利技术提供的一种基于变分模态分解的瞬变电磁探测信号的降噪方法中采用小波阈值降噪(WTD)对含随机噪声干扰的瞬变电磁信号的降噪结果图。图9是本专利技术提供的一种基于变分模态分解的瞬变电磁探测信号的降噪方法本文档来自技高网...

【技术保护点】
1.一种基于变分模态分解的瞬变电磁探测信号的降噪方法，其特征在于，包括：/n获取采空区强干扰瞬变电磁信号，预设瞬变电磁信号包括K个单分量信号，对强干扰瞬变电磁信号进行变分模态分解，得到强干扰瞬变电磁信号的变分约束函数；/n将拉格朗日惩罚算子和二次惩罚因子引入强干扰瞬变电磁信号的变分约束函数，得到增广拉格朗日函数，使变分约束函数不受约束；/n采用乘法算子交替方法(ADMM)，按预设方式交替迭代更新增广拉格朗日函数的参数，求解增广拉格朗日函数的鞍点，直至满足预设停止条件后停止迭代更新；其中，增广拉格朗日函数的参数至少包括瞬变电磁信号单分量信号、中心角频率和拉格朗日惩罚算子；/n提取迭代至增广拉格朗日函数满足预设停止条件时的全部单分量信号，剔除其中的噪声信号，将剩余单分量信号重组，得到去除噪声后的瞬变电磁信号。/n

【技术特征摘要】
1.一种基于变分模态分解的瞬变电磁探测信号的降噪方法，其特征在于，包括：
获取采空区强干扰瞬变电磁信号，预设瞬变电磁信号包括K个单分量信号，对强干扰瞬变电磁信号进行变分模态分解，得到强干扰瞬变电磁信号的变分约束函数；
将拉格朗日惩罚算子和二次惩罚因子引入强干扰瞬变电磁信号的变分约束函数，得到增广拉格朗日函数，使变分约束函数不受约束；
采用乘法算子交替方法(ADMM)，按预设方式交替迭代更新增广拉格朗日函数的参数，求解增广拉格朗日函数的鞍点，直至满足预设停止条件后停止迭代更新；其中，增广拉格朗日函数的参数至少包括瞬变电磁信号单分量信号、中心角频率和拉格朗日惩罚算子；
提取迭代至增广拉格朗日函数满足预设停止条件时的全部单分量信号，剔除其中的噪声信号，将剩余单分量信号重组，得到去除噪声后的瞬变电磁信号。


2.根据权利要求1所述的基于变分模态分解的瞬变电磁探测信号的降噪方法，其特征在于，在对强干扰瞬变电磁信号进行变分模态分解，得到强干扰瞬变电磁信号的变分约束函数的步骤中，
设瞬变电磁信号为x(t)，由K个单分量信号ui(t)组成，利用希尔伯特变换计算每个单分量信号的相关分析信号yi(t)，得到单分量信号的单侧频谱：



其中，δ(t)为狄克拉冲击函数，*为卷积符号；
根据调制和调解原理将相关分析信号yi(t)与指数混合进行修正，把ui(t)频率中心调谐到零频位置，得到相应的解调信号y′i(t)：



其中，ω...

【专利技术属性】
技术研发人员：戚庭野，卫会汝，冯国瑞，白锦文，赵德康，廉旭刚，张建忠，王朋飞，
申请(专利权)人：太原理工大学，
类型：发明
国别省市：山西;14