一种基于正交基函数处理标量磁异常梯度信号的方法技术

本发明专利技术涉及一种基于正交基函数处理标量磁异常梯度信号的方法，通过分析目标产生的磁异常梯度，设计了磁异常梯度的4个基函数。由于4个基函数为线性无关，通过格拉姆‑施密特正交化和归一化，获得标量磁异常梯度信号的正交基函数。根据正交基函数的性质和测量磁异常梯度信号，计算出对应正交基函数的系数。利用该系数设计了磁性目标的探测指数及目标方位参数。当探测指数大于设置阈值时，实现磁异常的探测。然后，根据正交基函数的系数获得目标的方位信息。本发明专利技术能够有效抑制磁噪声对目标探测的影响，提高输出信噪比，实现对弱磁异常的探测与定位。

【技术实现步骤摘要】
一种基于正交基函数处理标量磁异常梯度信号的方法
本专利技术所涉及到一种基于正交基函数处理标量磁异常梯度信号的方法。该方法可广泛用于能源矿藏勘测、水下各种管线探测、水下目标探测、水下考古、沉船勘测、扫雷反潜等领域。
技术介绍
地磁场是地球的一个天然的物理场，它有各种不同的起源，由不同变化规律的磁场成分叠加而成。按照场源位置划分，地磁场可以分为内源场和外源场。如果考虑地磁场随时间的变化特征，将随时间变化较快的地磁场成为地球的变化磁场，随时间变化较慢或者基本不变的地磁场成为地球的稳定磁场。同时地磁场是反映宇宙演变、地球演变、地质结构演变以及地震活动等过程的重要物理量之一。地磁场研究成果在航海、航空、航天、能源矿产、安全、考古等领域中有着广泛而重要的应用。磁性目标产生的磁场叠加到地磁场上，使得其周围空间的磁场分布发生变化，形成磁异常。通过对磁异常的探测与反演，可以实现目标的识别、定位与追踪。由于磁异常随距离的三次方快速衰减，这就导致远距离目标产生的磁异常通过掩埋在磁噪声中。传统的噪声抑制方法在一定程度上能够提高磁异常的探测能力。但当信噪比小于1dB，基于噪声抑制的探测方法对目标的探测能力就变得非常有限。同时，传统的探测方法不能给出目标的方位信息，不能满足对目标定位的需求。因此，需要有效的弱磁异常探测方法提高低信噪比下的磁性目标的探测与定位的能力。
技术实现思路
要解决的技术问题针对现有技术在低信噪比下的目标探测能力不足，本专利技术提供一种基于正交基函数处理标量磁异常梯度信号的方法，提高低信噪比下目标探测能力。技术方案一种基于正交基函数处理标量磁异常梯度信号的方法，其特征在于步骤如下：步骤1：采集磁场数据，获得磁性目标在X、Y和Z方向磁异常梯度Gx,Gy,Gz；步骤2：构建4个正交基函数，如下：其中，D＝|D|表示位于测线上的梯度传感器距离最近点的距离，R0＝|R0|表示目标到测线的最近点的距离；步骤3：利用正交基函数和磁场梯度信号，获得对应正交基函数的系数，如下：其中，m表示窗口的中心值，l表示数据处理窗口的长度；步骤4：利用磁异常梯度的正交基函数对应的系数，构建磁性目标的探测指数，如下：将E与阈值T进行比较来判断是否探测到磁异常：当E≥T，磁异常存在，即存在磁性目标；当E＜T，无磁异常，即不存在磁性目标；步骤5：当E≥T，根据正交基函数的系数的关系，构建目标的方位指数：当Dθ＞0时，目标位于磁梯度仪运动测线的北侧；当Dθ＜0时，目标位于磁梯度仪运动测线的南侧。有益效果为了解决低信噪比下弱磁异常探测能力不足，同时无法确定目标方位的问题，本专利技术提出了一种基于正交基函数处理标量磁异常梯度信号的方法，该方法分析目标产生的磁异常梯度，设计了磁异常梯度的4个正交基函数，根据正交基函数的性质和测量磁异常梯度信号，可以计算出对应正交基函数的系数，利用正交基函数的系数，构建磁性目标的探测指数及目标方位指数。该方法提高低信噪比下磁性目标的探测性能，同时能过确定目标的方位信息，提升了目标的探测与定位能力。由于目标产生的磁异常随距离的三次方快速衰减，远距离目标产生的磁异常信号通常掩埋在磁噪声中。本方法中，将采集的磁场梯度信号，通过构建的一组正交基函数进行处理，获得对应正交基函数的系数，根据该系数的性质，构建了弱磁异常的目标探测指数和目标方位指数，提高了弱磁异常的探测与定位能力。通过分解磁异常梯度信号，获得4个线性无关的基函数，通过格拉姆-施密特正交化和归一化，获得对应的正交基函数，利用正交基函数处理磁异常梯度信号，实现对弱磁异常的探测与定位。附图说明图1基于正交基函数的探测方法示意图图2磁性目标探测模型示意图图3磁性目标产生的磁异常梯度信号图4磁异常的探测结果具体实施方式现结合实施例、附图对本专利技术作进一步描述：一种基于正交基函数处理标量磁异常梯度信号的方法，包括以下步骤：步骤一：构建标量磁异常梯度的探测模型和探测坐标系，磁性目标位于坐标系的原点，其磁矩的单位适量为m；磁梯度传感器沿着平行于X方向的测线运动，目标到测线的最近点(CPA)的距离为R0＝|R0|；s表示R0在Y轴的投影大小，即最近点到目标的水平距离；h表示R0在Z轴的投影大小，即测线到目标的垂直距离。步骤二：采集磁场数据，获得磁性目标在X、Y和Z方向磁异常梯度Gx,Gy,Gz，可以消除地磁场空间分布的影响及随时间变化的影响。步骤三：获得标量磁异常梯度信号的正交基函数，其过程如下：Gx＝k(m·r)[5(u·r)D-uxr2]-kr2(u·m)D-kr2(u·r)mx(1)其中：μ0表示真空磁导率，r＝R0+D表示磁性目标到梯度传感器的位置矢量，r＝||r||，D＝|D|表示位于测线上的梯度传感器距离最近点的距离，M表示磁性目标的磁矩大小。m＝[mx,my,mz]表示目标的磁矩单位矢量，u＝[ux,uy,uz]表示地磁场的单位矢量。令其中：hx1＝p(5a2b2-uxa2-c-mxb2)，hx2＝p(5(a1b2+a2b1)-uxa1-mxb1)，hx3＝p(5a1b1-uxa2-c-mxb2)，hx4＝-p(uxa1+mxb1)。由于φ1(ω),φ2(ω),φ3(ω)和φ4(ω)构成的朗斯基行列式恒不等于0，即：φ1(ω),φ2(ω),φ3(ω)和φ4(ω)为线性无关的函数。因此，φ1(ω),φ2(ω),φ3(ω)和φ4(ω)是磁异常梯度的基函数。对基函数进行格拉姆-施密正交化和归一化，获得磁异常梯度的正交基函数，如下：因此，X方向上的磁异常梯度可以表示为：其中：同理，Y和Z方向上的磁异常梯度可以表示为：其中：hy1＝p(-uya2-myb2)hy2＝p(5θa2b2-a1uy-cθ-myb1)hy3＝p(5θ(a1b2+a2b1)-a2uy-myb2)hy4＝p(5θa1b1-a1uy-cθ-myb1)步骤四：根据正交基函数的性质，构建磁性目标探测指数E：由正交基函数的性质，可知：因此，根据公式(4)和(6)，可以获得系数ηxi：由于传感器的测量到的信号通常是离散信号，因此，公式(7)可以表示为：同理：磁异常梯度在Y和Z方向的系数可以表述为：利用磁异常梯度的正交基函数对应的系数，可以构建磁性目标的探测指数，如下：步骤五：通过阈值T，判断是否探测到磁异常。当E≥T，磁异常存在，即存在磁性目标；当E＜T，无磁异常，即不存在磁性目标。步骤六：当E≥T(即磁异常存在时)，根据正交基函数的系数的关系，构建目标的方位指数：当Dθ＞0时，目标位本文档来自技高网...

【技术保护点】
1.一种基于正交基函数处理标量磁异常梯度信号的方法，其特征在于步骤如下：/n步骤1：采集磁场数据，获得磁性目标在X、Y和Z方向磁异常梯度G

【技术特征摘要】
1.一种基于正交基函数处理标量磁异常梯度信号的方法，其特征在于步骤如下：
步骤1：采集磁场数据，获得磁性目标在X、Y和Z方向磁异常梯度Gx,Gy,Gz；
步骤2：构建4个正交基函数，如下：



其中，D＝|D|表示位于测线上的梯度传感器距离最近点的距离，R0＝|R0|表示目标到测线的最近点的距离；
步骤3：利用正交基函数和磁场梯度信号，获得对应正交基函数的系数，如下：



其中，m表示窗口的中...

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