一种利用时间相干性抑噪的磁异常检测方法技术

本发明专利技术公开了一种利用时间相干性抑制背景噪声的磁异常检测方法。该方法利用地磁背景噪声在间隔比较近的时间段内具有较强的的相干性的特点，针对同一传感器采集的不同时间段的信号，在频域自适应地计算地磁背景噪声在时间上的传递函数，利用传递函数消除目标信号中的地磁背景噪声，提高磁异信号的信噪比，以便检测到磁异常目标。本发明专利技术可以提高磁异常信号的信噪比，有助于找到磁性目标。磁性目标。磁性目标。

【技术实现步骤摘要】
一种利用时间相干性抑噪的磁异常检测方法


[0001]本专利技术涉及电子行业雷达
，涉及信号去噪、目标检测技术，属于微弱信号去噪技术。

技术介绍

[0002]人们称地球产生的磁场为地磁场，铁磁物体在地球上会影响自身周围的磁场分布。因为这个特性，磁异常探测技术被开发并不断发展，已有几十年的历史。由于磁场测量受的影响因素少，且磁场能穿透大部分固体和液体，磁异常探测技术被广泛应用于各种领域。如水下潜艇探测、目标车辆检测、地质资源勘探等。
[0003]自适应相干抑噪技术是前人利用地磁背景噪声的空间相干性而提出的一种抑噪技术。这种抑噪技术需要两组传感器分别放在两处相隔较远的位置。由于地磁背景噪声和磁异常目标信号的相干性极低，于是利用两处传感器采集的背景信号之间的相干性，可以从采集信号中估算出磁异常信号。由于不同采集卡、传感器之间存在不一致性问题，抑制噪声的工作会因此受到影响，偏离预期效果。
[0004]磁异常探测技术一般可以分为两大类，其中一类揭示了磁背景的性质，前面提到的自适应相干抑噪就属于这一类。另一类技术是利用目标信号的典型特征进行分析，这类方法中较为典型的方法是标准正交基函数(OBF)分解算法。标准正交基函数分解算法常常能在磁异常检测中提供重要的帮助，但是当信噪比过低的时候，标准正交基函数分解算法检测磁异常的能力受到限制。

技术实现思路

[0005]本专利技术提供了一种利用时间相干性抑制背景噪声并探测磁异常的方法
[0006]根据本专利技术的一个方面，提供了一种利用时间相干性抑制背景噪声的磁异常探测方法。该探测方法包括：
[0007]步骤1：利用磁传感器与采集卡采集得到磁信号数据y(n)；
[0008]步骤2：分别设置两个长度均为L的移动窗Wt、Wb，使两个窗口分别在信号上移动，得到Wb、Wt窗口截取的信号；
[0009]步骤3：将Wt窗内的信号作为目标信号，Wb窗内的信号作为背景信号，分别令为x(n,k)和x(n,k+Δk)，对目标信号x(n,k)和背景信号x(n,k+Δk)分别去均值，然后进行傅里叶变换，得到对应的频域信号X(ω,k)和X(ω,k+Δk)，其中，k表示窗口序数，ω表示频率；
[0010]步骤4：通过对两窗内信号的运算，即求信号x(n,k)和信号x(n,k+Δk)的互谱以及窗Wb内信号x(n,k+Δk)的自谱可以得到目标信号窗Wt与背景信号窗Wb中的背景噪声在时间轴上的传递函数H
k+Δk,k
(ω,k)；
[0011]步骤5：在频域中，利用背景信号以及步骤4得到的传递函数估计目标信号中的背景噪声，即通过H
k+Δk,k
(ω,k)
×
X(ω,k+Δk)估计窗Wt中信号的背景噪声，并用目标信号减去估计得到的背景噪声，得到抑噪后的目标频域信号，即用公式T(ω,k)＝X(ω,k)
‑
H
k+Δk,k
(ω,k)
×
X(ω,k+Δk)计算得到抑噪后的目标频域信号；
[0012]步骤6：对步骤5得到的目标频域信号进行傅里叶逆变换，得到抑噪后的目标时域信号t(n,k)；
[0013]步骤7：对步骤6得到的目标信号t(n,k)进行标准正交基函数分解，得到这段信号的能量E；
[0014]步骤8：保持目标信号窗口Wt位置固定不变，将背景信号窗Wb向右平移mb长度；
[0015]步骤9：重复步骤3～步骤8，直到背景信号窗Wb平移至信号的末尾，将信号从头到尾完整地遍历一遍；
[0016]步骤10：累加每次循环计算的能量E，然后求平均值，得到Wt窗口所截取的信号段的E_aver；
[0017]步骤11：将目标信号窗口Wt位置向右平移mt，然后重复步骤步骤3～步骤10，直到目标信号窗Wt平移至信号的末尾，将信号从头到尾完整地遍历一遍；
[0018]步骤12：累加并整理每次循环计算的能量E_aver，然后求平均值，得到整个信号段的全局平均值E_aver_all，对E_aver_all进行析并设定阈值，判断是否存在磁异常以及磁异常出现的时刻。
[0019]本专利技术的优点是：1、本专利技术通过利用背景噪声的时间相干性抑噪，只需要使用一组采集卡和传感器设备，相比传统的抑噪方法，实验所需的设备更少，节约了实验的成本。
[0020]2、本专利技术由于只使用一组采集卡和传感器，引入的干扰来源更少，相比传统的抑噪的方法，不会受多组采集卡传感器设备带来的信号不一致性的影响。
[0021]3、本专利技术利用时间相干性抑噪，由于数据由同一组传感器采集卡获取，所以在去噪过程中，不仅能对地磁背景噪声进行抑制，还对设备产生的噪声有一定的抑制效果。
[0022]4、本专利技术通过将目标信号与多组背景信号分别使用算法运算，最后取平均值，降低了过程中可能出现的误判概率。
附图说明
[0023]图1为利用时间相干性自适应抑噪检测磁异常方法的流程图；
[0024]图2为设备摆放示意图；
[0025]图3为窗口移动示意图；
具体实施方式
[0026]为使本专利技术的目的，技术方案和优点更加清晰明了，下面将结合具体实例，并参照附图,对本专利技术进一步详细说明。需要说明的是,在附图或说明书描述中,相似或相同的部分都使用相同的图号。附图中未绘出或描述的实现方式,为所属
中普通技术人员所知的形式。另外,虽然本文可提供包含特定值的参数的示范,但应了解,参数无需确切等于相应的值,而是可在可接受的误差容限或设计约束内近似于相应的值。实施例中提到的方向用语,例如“上”“下”“前”、“后”、“左”“右”等,仅是参考附图的方向。因此,使用的方向用语是用来说明并非用来限制本专利技术的保护范围。
[0027]本专利技术利用地磁背景噪声具有较高的时间相干性，自适应地滤除背景噪声。
[0028]在本专利技术的一个示例性实施例中，提供了一种利用时间相干性抑制地磁背景噪声的磁异常探测方法。图1为根据本专利技术实例利用时间相干性自适应抑噪探测磁异常方法的流程图。如图1所示，本实例利用时间相干性自适应抑噪探测磁异常方法包括：
[0029]步骤1：设置磁传感器，使传感器的探测范围覆盖磁性目标所在区域。
[0030]图2为传感器等设备的摆放示意图。按照图2中的方式摆放设备，确保在磁性目标经过传感器正前方的时候，传感器的探测范围能够覆盖磁性目标。
[0031]利用磁传感器与采集卡采集得到磁信号数据y(n)；
[0032]步骤2：分别设置两个长度均为L的移动窗Wt、Wb，使两个窗口分别在信号y(n)上移动，得到Wb、Wt窗口截取的信号；
[0033]步骤3：将Wt窗内的信号作为目标信号，Wb窗内的信号作为背景信号，并将其分别表示为x(n,k)和x(n,k+Δk)，两组信号具体可以写成如下形式：
[0034]x(n,k)＝t(n,k)+e(n,k)(1)
[0035]x(n,k+Δk)＝t(n,k+Δ本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种利用时间相干性抑噪的磁异常检测方法，其特征在于实现步骤如下：步骤1：利用磁传感器与采集卡采集得到磁信号数据y(n)；步骤2：分别设置两个长度均为L的移动窗Wt、Wb，使两个窗口分别在信号y(n)上移动，得到Wb、Wt窗口截取的信号；步骤3：将Wt窗内的信号作为目标信号，Wb窗内的信号作为背景信号，对目标信号和背景信号分别去均值，然后进行傅里叶变换，得到对应的频域信号；步骤4：分别计算目标信号与背景信号的互谱，计算背景信号的自谱，然后求目标信号和背景信号中的背景噪声在时间轴上的传递函数；步骤5：在频域中，利用背景信号以及步骤4得到的传递函数估计目标信号中的背景噪声，并用目标信号减去估计得到的背景噪声，得到抑噪后的目标频域信号；步骤6：对步骤5得到的目标频域抑噪信号进行傅里叶逆变换，得到抑噪后的目标时域信号t(n,k)；步骤7：对步骤6得到的目标信号t(n,k)进行标准正交基函数分解，得到这段信号的能量E；步骤8：保持目标信号窗Wt位置固定不变，将背景信号窗Wb向右平移mb长度；步骤...

【专利技术属性】
技术研发人员：王成栋，彭根斋，雷廷宇，闫欢，陈勇，杨冬辰，周家齐，
申请(专利权)人：中国电子科技集团公司第九研究所，
类型：发明
国别省市：