本發(fā)明屬于磁探測技術領域，更具體地，涉及一種基于矢量磁梯度計的運動式磁性目標定位方法。

背景技術：

磁探測法常用于地質勘探、生物醫(yī)療、導航定位、未爆物探測等領域；在這些領域應用高分辨率的磁場總強度探測方法可方便地進行磁性物體的移動探測，但實現(xiàn)磁性物體的精確定位較為困難?；诖盘荻扔嫷奶綔y方法能夠有效實現(xiàn)目標的定位，但矢量磁梯度計的安裝制作需要保證構成矢量磁梯度計的兩組矢量磁強計的軸線嚴格平行、電氣性能嚴格對稱；現(xiàn)有技術通過加工工藝來保證軸線平行度以及電氣對稱性，存在結構誤差，難以實現(xiàn)高分辨率的探測，實際的探測分辨率在10nT～100nT的量級，難以在未爆物的探測領域應用。

技術實現(xiàn)要素：

針對現(xiàn)有技術的以上缺陷或改進需求，本發(fā)明提供了一種基于矢量磁梯度計的運動式磁性目標定位方法，其目的在于解決矢量磁梯度計的結構誤差補償問題。

為實現(xiàn)上述目的，按照本發(fā)明的一個方面，提供了一種基于矢量磁梯度計的運動式磁性目標定位方法，利用目標磁場全張量信息進行目標定位，包括如下步驟：

(1)對構成矢量磁梯度計的兩組矢量磁強計的結構誤差進行補償，使兩組矢量磁強計的信號輸出保持一致；

(2)使矢量磁梯度計在磁異?？臻g繞豎直軸線旋轉，采集磁梯度計在旋轉方位角的磁場方向梯度數(shù)據(jù)，并根據(jù)方向梯度與方位角的關系獲取磁場全張量；

(3)利用磁場全張量，根據(jù)GR＝-3B進行磁性目標定位；其中，B是磁異常量，R是探測點相對于磁性目標的位置；G是指磁場全張量。

優(yōu)選地，上述基于矢量磁梯度計的運動式磁性目標定位方法，其步驟(1)包括如下子步驟：

(1.1)使矢量磁梯度計在均勻磁場的空間做三自由度旋轉，根據(jù)構成矢量磁梯度計的兩組磁強計在旋轉方位角處對應軸的數(shù)據(jù)相等關系，獲取矢量磁梯度計的三軸補償系數(shù)；

(1.2)根據(jù)補償系數(shù)對構成磁梯度計的兩組矢量磁強計的結構誤差進行補償，使得兩組矢量磁強計的輸出一致。

優(yōu)選地，上述基于矢量磁梯度計的運動式磁性目標定位方法，其步驟(1.1)包括如下子步驟：

(1.1.1)獲取三分量的磁梯度計輸出

其中，B_ax、B_ay、B_az分別是第一矢量磁強計a的x、y、z三軸數(shù)據(jù)；

B_bx、B_by、B_bz分別是第二矢量磁強計b的x、y、z三軸數(shù)據(jù)，ΔB_x是指磁梯度計x軸的數(shù)據(jù)，ΔB_y是指磁梯度計y軸的數(shù)據(jù)，ΔB_z是指磁梯度計z軸的數(shù)據(jù)；

(1.1.2)獲取空間兩點的實際磁場梯度與磁梯度計輸出之間的誤差

其中，ΔB^*_x是指補償后的磁梯度計x軸的數(shù)據(jù)，ΔB^*_y是指補償后的磁梯度計y軸的數(shù)據(jù)，ΔB^*_z是指補償后的磁梯度計z軸的數(shù)據(jù)，k_mn(m,n＝1、2、3)為磁梯度計的三軸補償系數(shù)；

(1.1.3)使矢量磁梯度計在均勻磁場的空間做三自由度旋轉，采集三軸磁強計和矢量磁梯度計的輸出樣本，根據(jù)輸出樣本建立三軸補償系數(shù)方程；

其中，j＝1、2、3……M為樣本序號；ΔB_xj、ΔB_yj、ΔB_zj是指樣本j中未經(jīng)過補償?shù)氖噶看盘荻扔嫷膞、y、z軸數(shù)據(jù)；

(1.1.4)求解上述三軸補償系數(shù)方程，獲得磁梯度計的三軸補償系數(shù)k_mn(m,n＝1、2、3)。

優(yōu)選地，上述基于矢量磁梯度計的運動式磁性目標定位方法，其步驟(2)包括如下子步驟：

(2.1)使矢量磁梯度計在磁異?？臻g繞豎直軸線旋轉，采集矢量磁梯度計在各旋轉方位角的磁場方向梯度數(shù)據(jù)；

(2.2)將矢量磁梯度計軸線與豎直軸線的夾角θ、矢量磁梯度計轉動的方位角與矢量磁梯度計x方向的梯度之間的關系作為參量，建立方程：

其中，L是指第一矢量磁強計a與第二矢量磁強計b的距離，G_xx指x軸方向的磁場分量在x方向上的梯度，G_xy指x軸方向的磁場分量在y方向上的梯度，G_xz指x軸方向的磁場分量在z方向上的梯度；

(2.3)根據(jù)矢量磁梯度計的旋轉方位角與第一矢量磁強計a和第二矢量磁強計b的x軸輸出的對應關系，建立方程組：

其中，B_axj是指旋轉方位角為時第一矢量磁強計a的x軸數(shù)據(jù)；B_bxj是指旋轉方位角為時第二矢量磁強計b的x軸數(shù)據(jù)；

(2.4)求解所述方程組，獲得G_xx、G_xy、G_xz；

(2.5)以矢量磁強計旋轉方位角與第一矢量磁強計a和第二矢量磁強計b的y軸數(shù)據(jù)對應的關系作為參量，重復步驟(2.2)～(2.3)，獲得G_yx、G_yy、G_yz；

(2.6)以矢量磁強計旋轉方位角與磁強計a和磁強計b的z軸數(shù)據(jù)對應的關系作為參量，重復步驟(2.2)～(2.3)，獲得G_zx、G_zy、G_zz；

(2.7)獲取磁場全張量

優(yōu)選地，上述基于矢量磁梯度計的運動式磁性目標定位方法，探測點相對于磁性目標的位置

其中，B_x是指磁異常量B的x軸方向分量、B_y是指磁異常量B的y軸方向分量、B_z是指磁異常量B的z軸方向分量。

總體而言，通過本發(fā)明所構思的以上技術方案與現(xiàn)有技術相比，能夠取得下列有益效果：

(1)本發(fā)明提供的基于矢量磁梯度計的運動式磁性目標定位方法，由于采用三軸系數(shù)補償?shù)姆椒ǎM行矢量磁梯度計的結構誤差補償；有效降低了空間梯度的測量誤差，能夠更加準確地獲得目標磁場全張量信息；

(2)本發(fā)明提供的基于矢量磁梯度計的運動式磁性目標定位方法，由于對矢量磁梯度計進行了補償校正，能有效提高目標磁場梯度的測量精度，因此具有更加準確實現(xiàn)目標探測定位的效果。

附圖說明

圖1是實施例中矢量磁梯度計結構誤差補償?shù)氖疽鈭D；

圖2是實施例中轉動矢量磁梯度計進行探測定位的示意圖；

在所有附圖中，相同的附圖標記用來表示相同的元件或結構，其中：a-第一矢量磁強計，b-第二矢量磁強計，c-探棒。

具體實施方式

為了使本發(fā)明的目的、技術方案及優(yōu)點更加清楚明白，以下結合附圖及實施例，對本發(fā)明進行進一步詳細說明。應當理解，此處所描述的具體實施例僅僅用以解釋本發(fā)明，并不用于限定本發(fā)明。此外，下面所描述的本發(fā)明各個實施方式中所涉及到的技術特征只要彼此之間未構成沖突就可以相互組合。

實施例中，如圖1所示的，矢量磁梯度計由平行放置的矢量磁強計構成；其中，a、b兩個模塊是指第一矢量磁強計a、第二矢量磁強計b；B_ax、B_ay、B_az是磁強計a的x、y、z三軸的輸出數(shù)據(jù)；B_bx、B_by、B_bz是磁強計b的x、y、z三軸的輸出數(shù)據(jù)；α、β、γ是第一矢量磁強計a與第二矢量磁強計b的三軸夾角。

實施例提供的基于矢量磁梯度計的運動式磁性目標定位方法，包括如下步驟：

步驟(1)：矢量磁梯度計的結構誤差補償：

三分量的磁梯度計輸出為：

由于存在安裝結構誤差和電氣一致性的誤差，空間兩點的實際磁場梯度與磁梯度計輸出之間的誤差為：

其中，ΔB^*_x、ΔB^*_y、ΔB^*_z為磁強計a和b的實際磁場x、y、z三分量的差值，k_mn(m,n＝1、2、3)為磁場梯度的補償系數(shù)；

使矢量磁梯度計在均勻磁場的空間做三自由度旋轉，采集三軸磁強計與矢量磁梯度計的輸出樣本，以求解磁場梯度補償系數(shù)；

由于磁場均勻空間中真實的磁梯度為零，三分量磁場差值與三軸補償系數(shù)的關系具體如下：

其中，j＝1、2、3……M為樣本序號，求解上述方程，獲得補償系數(shù)。

步驟(2)：根據(jù)方位角與磁場的方向倒數(shù)之間的關系計算獲取磁場全張量：

實施例中，將矢量磁梯度計傾斜懸掛在運動探測平臺(如旋翼無人機)上，使磁梯度計旋轉，采集磁梯度計在處于方位角時的磁場方向梯度數(shù)據(jù)，根據(jù)方向梯度與方位角的關系，建立方程解算目標的磁場全張量，具體如下：

如圖2所示，當磁梯度計軸線與豎直軸線的夾角為θ、磁梯度計轉動的方位角為時，根據(jù)全微分的概念有：

聯(lián)立方位角對應的關系式，建立以下方程組：

求解上述方程組，獲得G_xx、G_xy、G_xz；

根據(jù)矢量磁梯度計的旋轉方位角與矢量磁梯度計的y軸梯度對應的關系方程組：

獲得G_yx、G_yy、G_yz；

根據(jù)矢量磁梯度計的旋轉方位角與矢量磁梯度計的z軸梯度對應的關系方程組：

獲得G_zx、G_zy、G_zz；

由此獲得磁場全張量

步驟(3)：磁性目標定位：

依據(jù)GR＝-3B，利用磁場全張量進行磁性目標定位；其中，B是磁異常量，R是探測點相對于磁性目標的位置；

獲得相對位置R，實現(xiàn)目標的定位。

本領域的技術人員容易理解，以上所述僅為本發(fā)明的較佳實施例而已，并不用以限制本發(fā)明，凡在本發(fā)明的精神和原則之內所作的任何修改、等同替換和改進等，均應包含在本發(fā)明的保護范圍之內。