基于标量磁力仪阵列的弱磁异常自适应实时探测方法技术

本发明专利技术公开了一种基于标量磁力仪阵列的弱磁异常自适应实时探测方法，用于解决现有磁异常探测方法探测概率低的技术问题。技术方案是通过经验模式分解，将采集的磁场信号分解成多个固有模态函数(IMF)，根据IMF的性质，重构为两部分：磁噪声和磁异常信号。由于经验模式分解的模态混叠问题，重构的磁异常信号包含部分噪声。利用重构的磁噪声实时更新噪声的概率密度函数。基于更新的概率密度函数，利用一个移动的窗口计算对应的地磁熵。当地磁熵小于设置阈值时，实现磁异常的探测。本发明专利技术提高了磁异常探测的概率，提升了低信噪比下目标探测能力。

Adaptive real-time detection method of weak magnetic anomaly based on scalar magnetometer array

【技术实现步骤摘要】
基于标量磁力仪阵列的弱磁异常自适应实时探测方法
本专利技术涉及一种磁异常探测方法，特别涉及一种基于标量磁力仪阵列的弱磁异常自适应实时探测方法。
技术介绍
文献“Magneticanomalydetectionusingentropyfilter，Measurementscienceandtechnology,2008,Vol19(4),045205”公开了一种基于地磁熵的磁异常探测方法。该方法是基于磁性目标产生的磁异常引起磁噪声模式变化的机理构建了一种磁异常探测方法。该方法将磁场背景噪声视为平稳随机过程，以概率密度函数作为先验条件，以地磁熵作为磁异常探测的准则。该方法在探测过程中，算法复杂度低，计算过程简单，能够满足实时探测的需求。文献所述方法以背景磁噪声的概率密度函数作为先验条件当环境条件发生变化时或者基于运动平台的探测时，作为先验条件的背景磁噪声的概率密度函数的获取将变得困难，其检测性能会因概率密度函数的不准确而受到影响，导致该方法的探测概率降低。
技术实现思路
为了克服现有磁异常探测方法探测概率低的不足，本专利技术提供一种基于标量磁力仪阵列的弱磁异常自适应实时探测方法。该方法通过经验模式分解，将采集的磁场信号分解成多个固有模态函数(IMF)，根据IMF的性质，重构为两部分：磁噪声和磁异常信号。由于经验模式分解的模态混叠问题，重构的磁异常信号包含部分噪声。利用重构的磁噪声实时更新噪声的概率密度函数。基于更新的概率密度函数，利用一个移动的窗口计算对应的地磁熵。当地磁熵小于设置阈值时，实现磁异常的探测。本专利技术提高了磁异常探测的概率，提升了低信噪比下目标探测能力。本专利技术解决其技术问题所采用的技术方案：一种基于标量磁力仪阵列的弱磁异常自适应实时探测方法，其特点是包括以下步骤：(a)利用经验模式分解对采集的磁场差值进行自适应分解，获得多个固有模态函数和一个残差项式中，ci(t)表示第i个固有模态函数分量，r(t)表示残差分量。d表示经验模式分解的层数；(b)根据固有模态函数的性质，分别重构磁噪声和磁异常信号：式中，N(t)表示重构的磁噪声，S(t)表示重构的磁异常信号；(c)利用重构的磁噪声N(t)，统计分析并更新该噪声的概率密度函数p(·)；(d)设置滑动窗口的长度L，根据更新的噪声概率密度函数，计算重构磁异常信号S(t)的地磁熵：式中，i表示第i时刻。本专利技术的有益效果是：该方法通过经验模式分解，将采集的磁场信号分解成多个固有模态函数(IMF)，根据IMF的性质，重构为两部分：磁噪声和磁异常信号。由于经验模式分解的模态混叠问题，重构的磁异常信号包含部分噪声。利用重构的磁噪声实时更新噪声的概率密度函数。基于更新的概率密度函数，利用一个移动的窗口计算对应的地磁熵。当地磁熵小于设置阈值时，实现磁异常的探测。本专利技术由于采用经验模式分解，将采集的磁场差值进行自适应分解，获得多个固有模态函数和一个残差项，利用分解得到的固有模态函数重新构建了磁噪声部分，通过对构建的磁噪声进行统计分析，获得其近似的概率密度函数，解决了由外部环境条件发生变化或者基于运动平台的磁异常探测过程中磁噪声概率密度函数获取困难的关键问题，提高了磁异常探测的概率，提升了低信噪比下目标探测能力。下面结合附图和具体实施方式对本专利技术作详细说明。附图说明图1是本专利技术基于标量磁力仪阵列的弱磁异常自适应实时探测方法的流程图。图2是本专利技术方法一维磁力仪阵列的目标探测示意图。图3是远距离目标产生的磁异常。图4是本专利技术方法弱磁异常探测的结果示意图。具体实施方式参照图1-4。本专利技术基于标量磁力仪阵列的弱磁异常自适应实时探测方法具体步骤如下：1、磁场差值的采集。采用两个磁力仪构成一维标量磁传感器阵列。阵列中的传感器间距Bx＝1m，磁传感器采用灵敏度为0.6pT的CS-L光泵磁力仪。磁性目标在水平面内移动，其磁矩为(0,0,25)A·m2。目标从起始位置(-180,30,-1)m，以匀速运动的方式到达终点位置(180,30,-1)m。2、磁噪声与磁异常信号的重构。利用经验模式分解对采集的磁场差值进行自适应分解，分解层数d＝8，获得8固有模态函数和1个残差项：式中，ci(t)表示第i个固有模态函数分量，r(t)表示残差分量。d表示经验模式分解的层数。根据固有模态函数的性质，固有模态函数的前4项用于重构磁噪声，后4项用于重构磁异常信号：式中，N(t)表示重构的磁噪声，S(t)表示重构的磁异常信号；3、磁噪声的概率密度函数更新。利用重构的磁噪声N(t)，统计分析并更新该噪声的概率密度函数p(·)；4、地磁熵的计算。设置滑动窗口的长度L，根据更新的噪声概率密度函数p(·)，计算重构磁异常信号S(t)的地磁熵：式中，i表示第i时刻；5、磁异常探测。通过阈值T，判断是否探测到磁异常。当I(Si)＜T，磁异常存在，即存在磁性目标；当I(Si)≥T，无磁异常，即不存在磁性目标。本文档来自技高网...

【技术保护点】
1.一种基于标量磁力仪阵列的弱磁异常自适应实时探测方法，其特征在于包括以下步骤：/n(a)利用经验模式分解对采集的磁场差值进行自适应分解，获得多个固有模态函数和一个残差项/n

【技术特征摘要】
1.一种基于标量磁力仪阵列的弱磁异常自适应实时探测方法，其特征在于包括以下步骤：
(a)利用经验模式分解对采集的磁场差值进行自适应分解，获得多个固有模态函数和一个残差项



式中，ci(t)表示第i个固有模态函数分量，r(t)表示残差分量；d表示经验模式分解的层数；
(b)根据固有模态函数的性质，分别重构磁噪...

【专利技术属性】
技术研发人员：樊黎明，赵维娜，王惠刚，刘星，刘建国，
申请(专利权)人：西北工业大学青岛研究院，西北工业大学，
类型：发明
国别省市：山东;37