馬劍飛，顏冰，林春生

(海軍工程大學(xué) 兵器工程學(xué)院，湖北 武漢 430033)

0 引言

磁異常探測技術(shù)具有不受空氣、海水、泥沙、土壤等介質(zhì)影響的優(yōu)勢，在勘探礦物、航空反潛以及未爆彈探測等方面具有廣泛而深入的應(yīng)用，其作為一種被動探測技術(shù)具有隱蔽性好的特點(diǎn)[1-3]。

目前主要的標(biāo)量磁異常探測技術(shù)大致分為兩類:第一類是目標(biāo)磁模型正交基函數(shù)(OBF)與測量信號乘積的積分方法。文獻(xiàn)[4]提出以磁偶極子為模型的標(biāo)量OBF檢測器，該方法基于磁偶極子模型并結(jié)合目標(biāo)與測量系統(tǒng)的相對運(yùn)動信息對信號進(jìn)行正交基分解，適用于隱蔽磁性目標(biāo)的檢測；文獻(xiàn)[5]提出了以線列陣偶極子為模型的主成分檢測器，該方法的OBF是在特定距離條件下仿真計(jì)算得到，不適用于對未知距離條件下的目標(biāo)檢測。第二類是基于信號和噪聲統(tǒng)計(jì)特性不同的檢測方法，包括最小熵檢測器(MED)[6]以及高階過零檢測器[7]，這類方式相對于OBF檢測器具有計(jì)算量較小的優(yōu)點(diǎn)，但需要對環(huán)境噪聲進(jìn)行長時(shí)間的統(tǒng)計(jì)特性分析，而且檢測能力比較有限。在實(shí)際的地磁環(huán)境中，噪聲功率譜一般為1/fα結(jié)構(gòu)(f表示信號頻率，α表示功率譜衰減常數(shù))，文獻(xiàn)[8]以自回歸(AR)模型白化濾波器對環(huán)境噪聲以及OBF進(jìn)行了校正，該方法在1/fα環(huán)境噪聲條件下相比于OBF檢測器取得了更優(yōu)的檢測結(jié)果，但是求解AR模型白化濾波器需要對噪聲進(jìn)行長時(shí)間的統(tǒng)計(jì)特性分析，也要求噪聲是平穩(wěn)的，但實(shí)際的測量難以滿足這一條件。

目前主要的磁性目標(biāo)遠(yuǎn)程探測方法主要是基于磁標(biāo)量或者標(biāo)量磁梯度的測量，標(biāo)量測量所獲取的目標(biāo)磁場信息較為有限，不能從磁場信息的角度充分解釋目標(biāo)；矢量磁測量雖然可以提高對目標(biāo)磁場信息的解釋能力，但地磁場相對于磁性目標(biāo)產(chǎn)生的磁異常為強(qiáng)磁場，測量平臺輕微的晃動就會導(dǎo)致矢量磁傳感器輸出產(chǎn)生很大的變化。

隨著超導(dǎo)磁力儀技術(shù)的發(fā)展，基于矢量磁梯度的磁異常探測方式愈發(fā)受到重視[9-10]。矢量磁梯度相對磁標(biāo)量探測以及磁矢量探測的優(yōu)勢在于：1)測量結(jié)果當(dāng)?shù)卮艃A角和地磁偏角無關(guān)；2)具有更豐富的磁場信息，可以利用幾個(gè)測量點(diǎn)的磁梯度張量快速反演目標(biāo)的位置和磁矩[11-13]；3)具備很強(qiáng)的共模抑制能力，特別是當(dāng)測量平臺晃動時(shí)，測量值可以有效避免地磁場和背景磁場噪聲影響[14]；4)日變效應(yīng)會導(dǎo)致假的磁異常，需要使用額外的基站進(jìn)行校正，磁梯度張量不受日變效應(yīng)影響，自然環(huán)境的磁梯度一般小于0.02 nT/m，而且在一個(gè)區(qū)域內(nèi)的值相對穩(wěn)定[15]。

在磁梯度OBF探測算法層面，文獻(xiàn)[4]提出了一種基于5個(gè)正交基的標(biāo)量磁梯度OBF檢測方法，該方法能夠有效消除環(huán)境磁噪聲的影響，可以顯著提升檢測的信噪比。文獻(xiàn)[16]在此基礎(chǔ)提出了一種基于4個(gè)正交基的標(biāo)量磁梯度OBF檢測方法，該方法消耗更小的計(jì)算資源，而且具有更好的魯棒性，限制條件是目標(biāo)與系統(tǒng)之間距離大于磁梯度陣列基線長度的3倍(一般情況下此條件是很容易滿足的)。文獻(xiàn)[17]進(jìn)一步研究了消除觀測信號趨勢和偏差的方法，并評估了觀測窗長度對算法性能的影響，發(fā)現(xiàn)當(dāng)觀測窗長為最近點(diǎn)(CPA)距離的1～4倍時(shí)，算法的檢測性能相對于最優(yōu)窗長(2～2.5倍CPA距離)的檢測性能不會發(fā)生明顯下降。

本文旨在研究一種基于矢量磁梯度測量的磁性目標(biāo)遠(yuǎn)程探測技術(shù)，可以在低信噪比條件下能夠提高提升探測的穩(wěn)定性(降低檢測虛警率)以及準(zhǔn)確估計(jì)目標(biāo)的位置參數(shù)和磁矩參數(shù)。

1 OBF檢測器

1.1 傳統(tǒng)的磁標(biāo)量OBF檢測器

當(dāng)探測距離大于2.5倍磁性目標(biāo)長度時(shí)，磁性目標(biāo)可以視為一個(gè)磁偶極子[18]，此時(shí)距離磁性目標(biāo)r處的磁場矢量B可以表示為

(1)

式中：μ0為空間磁導(dǎo)率，μ0=4π×10-7H/m；M表示目標(biāo)磁矩矢量；r表示磁性目標(biāo)到傳感器的位移矢量r的模值。當(dāng)B?T時(shí)(T表示地磁場矢量)，標(biāo)量磁場S的測量可視為磁場矢量在地磁場方向上的投影，即

(2)

圖1所示為運(yùn)動載體沿直線運(yùn)動時(shí)與磁性目標(biāo)(目標(biāo)視為偶極子)的相對幾何關(guān)系。圖1中，坐標(biāo)原點(diǎn)O表示目標(biāo)的位置(偶極子)，X軸表示運(yùn)動載體沿直線運(yùn)動的方向，Y軸表示運(yùn)動載體的左舷方向，Z軸表示豎直向上方向，R0表示傳感器和磁性目標(biāo)之間的CPA距離，D0表示傳感器當(dāng)前位置和CPA點(diǎn)之間的距離，θM表示目標(biāo)磁矩方向與Z軸的夾角，θT表示地磁場方向與Z軸的夾角，φM表示目標(biāo)磁矩矢量在OXY面的投影與X軸的夾角，φT表示地磁場方向在OXY平面的投影與X軸的夾角。

(1)式代入(2)式，標(biāo)量磁場可表示為基函數(shù)φ(ω)的加權(quán)和形式：

(3)

對基函數(shù)φ(ω)進(jìn)行施密特正交化，可得

(4)

圖2所示為磁偶極子模型對應(yīng)的正交基函數(shù)f1(ω)、f2(ω)與f3(ω)。

圖2 正交基函數(shù)

由(4)式可知，正交化系數(shù)aj*可通過OBF與測量信號乘積的積分進(jìn)行解算：

(5)

對(5)式離散化處理，可得能量系數(shù)

(6)

式中：Δω表示空間采樣長度，Δω=ωi+1-ωi，ωi表示i時(shí)刻X軸方向距離與CPA距離的比值，i的取值范圍為[-k,k]。在[ω-k,ωk]觀測窗區(qū)間內(nèi)，信號的能量Eb可以表示為能量系數(shù)的平方和

(7)

1.2 磁矢量OBF檢測器及性能分析

當(dāng)目標(biāo)磁場遠(yuǎn)小于地磁場時(shí)，磁標(biāo)量的測量可等效于目標(biāo)磁場矢量在地磁場方向上的投影，同理三分量磁場的測量SX、SY和SZ也可以看作是目標(biāo)磁場在X軸、Y軸和Z軸方向上的投影。以磁矢量測量SX為例，結(jié)合(3)式可知

(8)

式中：bX,j為分量SX對應(yīng)的基函數(shù)系數(shù)，bX,1=3sinθM·cosφMsinθXcosφX，bX,2=3(cosθMsinθXcosφX+sinθMcosφMcosθX)，bX,3=3cosθMcosθX，bX,4=-sinθM·sinθXcos(φM-φX)-cosθMcosθX，θX表示矢量磁力儀X軸方向與坐標(biāo)系Z軸的夾角(θX=π/2 rad)，φX表示矢量磁力儀X軸方向與坐標(biāo)系X軸的夾角(φX=0 rad)。在[ω-k,ωk]的觀測窗區(qū)間內(nèi)，X軸方向的OBF能量檢測器為輸入信號的OBF加權(quán)能量：

(9)

式中：αX,j*為分量SX對應(yīng)的能量系數(shù)。

同理，也可以構(gòu)造Y軸方向的OBF檢測器EY以及Z軸方向的OBF檢測器EZ，則矢量OBF能量檢測器的輸出可以構(gòu)造為

Ev=EX+EY+EZ.

(10)

下面以X軸的分量SX為例(無晃動條件下與磁標(biāo)量測量等效)，推導(dǎo)OBF檢測結(jié)果所服從的概率分布并對檢測性能的影響因素進(jìn)行分析。

假設(shè)X軸測量的環(huán)境噪聲服從高斯獨(dú)立同分布SX(ωi)～N(0,σ2)，σ為高斯噪聲的標(biāo)準(zhǔn)差，代入(6)式可得

(11)

由于

(12)

故

αX,j*～N(0,σ2Δω).

(13)

進(jìn)一步有

(14)

下面證明αX,1、αX,2和αX,3之間相互獨(dú)立，以αX,1和αX,2為例，αX,1和αX,2的期望分別為

E(αX,1)=0，E(αX,2)=0，

(15)

αX,1和αX,2乘積的期望為

(16)

又因?yàn)?/p>

(17)

故

E(αX,1αX,2)=0，

(18)

E(αX,1αX,2)-E(αX,1)E(αX,2)=0.

(19)

由高斯分布相互獨(dú)立的充要條件可知，αX,1、αX,2和αX,3之間相互獨(dú)立，聯(lián)立(14)式，可得

(20)

由卡方分布的性質(zhì)可知，X軸測量噪聲經(jīng)過OBF加權(quán)后其能量的期望EX,n和方差DX，n分別為

(21)

(22)

進(jìn)一步，綜合(9)式和(21)式可知，噪聲OBF加權(quán)能量與信號OBF加權(quán)能量的比值可表示為

(23)

式中：EX為X軸測量的信號能量；η越小，表示信號與噪聲的區(qū)分度越高。

聯(lián)立(9)式、(21)式可知，信號OBF檢測能量僅與信號本身的能量相關(guān)，而高斯噪聲的OBF檢測能量卻與空間采樣長度有關(guān)，因此空間采樣長度Δω越小則OBF檢測器對噪聲和信號的區(qū)分度越高。另外由(22)式可知，空間采樣長度越小則同分布噪聲經(jīng)過OBF檢測器后概率分布的方差越小。因此不難得出結(jié)論：1)在同樣測量環(huán)境條件下，空間采樣長度越小則OBF檢測器對環(huán)境噪聲和目標(biāo)信號的區(qū)分越明顯；2)空間采樣長度越小則檢測器的穩(wěn)定性越好。

假設(shè)三軸測量的噪聲服從獨(dú)立同分布且三軸測量的信號能量相同，對于矢量檢測器而言，

(24)

(25)

式中：En和Dn分別為矢量檢測器噪聲經(jīng)OBF加權(quán)后其能量的期望和方差。則矢量檢測器的噪聲OBF加權(quán)能量En和信號OBF加權(quán)能量Ev的比值為

(26)

由(26)式可以看出，在同等空間采樣長度條件下，矢量測量相對于標(biāo)量檢測并不能改變區(qū)分度，但聯(lián)立(22)式和(25)式可得

3D

En=DX,n/EX,n.

(27)

1.3 磁梯度OBF檢測器

鑒于矢量磁梯度測量兼具磁標(biāo)量檢測的抗晃動性以及磁矢量檢測信息豐富的優(yōu)勢，據(jù)此在傳統(tǒng)標(biāo)量磁梯度OBF檢測器的基礎(chǔ)之上提出了一種矢量磁梯度OBF檢測器，該算法可在低信噪比以及平臺晃動條件下可靠檢測磁性目標(biāo)，并能夠反演得到目標(biāo)的矢量磁梯度信號。文獻(xiàn)[16]提出了一種基于4個(gè)OBF的標(biāo)量磁場梯度的OBF檢測器，標(biāo)量磁梯度測量等效于磁場矢量在地磁場方向上投影的差分G，即

(28)

式中：Br+ΔZ為距離磁性目標(biāo)r+ΔZ處的磁場矢量；Br為距離磁性目標(biāo)r處的磁場矢量。

磁梯度傳感器由兩個(gè)相距ΔZ的磁傳感器構(gòu)成，圖3所示為磁梯度傳感器沿直線運(yùn)動時(shí)與磁性目標(biāo)(目標(biāo)視為偶極子)的相對幾何關(guān)系，其中X軸為運(yùn)動載體沿直線運(yùn)動的方向，Y軸為運(yùn)動載體的左舷方向，Z軸為豎直向上方向，R0為傳感器和磁性目標(biāo)之間的CPA距離。

圖3 磁梯度傳感器與目標(biāo)的相對運(yùn)動關(guān)系

標(biāo)量磁梯度的正交基形式為

(29)

圖4 磁梯度正交基函數(shù)

(30)

標(biāo)量磁梯度檢測器的輸出Es可表示為能量系數(shù)的平方和，

(31)

記矢量磁梯度的三分量分別為GZ,X、GZ,Y和GZ,Z，結(jié)合矢量磁場OBF檢測器的推導(dǎo)過程可知，矢量磁梯度OBF檢測器的輸出可以構(gòu)造為矢量磁梯度三分量對應(yīng)能量系數(shù)的平方和形式

(32)

由(27)式可知，矢量磁梯度OBF檢測器相比于標(biāo)量磁梯度OBF檢測器改善了檢測的穩(wěn)定性，矢量OBF檢測器反演的磁梯度信號可以表示為

(33)

2 磁性目標(biāo)參數(shù)估計(jì)

反演的磁梯度信號雖然相比于原始的測量信號在很大程度上提高了信噪比，但受檢測器窗長以及噪聲非嚴(yán)格高斯等因素的影響，反演信號相較于真實(shí)信號產(chǎn)生了一定的畸變。注意到由目標(biāo)產(chǎn)生的信號是可通過偶極子模型進(jìn)行擬合的，而信號畸變的部分則不可以，因此可以通過偶極子模型對反演信號進(jìn)行擬合從而實(shí)現(xiàn)對目標(biāo)磁矩參數(shù)和位置參數(shù)的估計(jì)。

以磁性目標(biāo)中心為原點(diǎn)建立坐標(biāo)系，X為運(yùn)動載體沿直線運(yùn)動的方向，Y為運(yùn)動載體的左舷方向，Z為豎直向上方向，磁性目標(biāo)等效于單偶極子模型，設(shè)磁傳感器l的坐標(biāo)為(xl，yl，zl)，則i時(shí)刻磁傳感器l的測量為

(34)

式中：MX、MY和MZ分別表示目標(biāo)在X軸方向、Y軸方向和Z軸方向的磁矩；磁偶極子系數(shù)的計(jì)算系數(shù)分別為

bX,i,l=aY,i,l，cX,i,l=aZ,i,l，cY,i,l=bZ,i,l，

Fi·M=ΔBi=Bi,1-Bi,2，

(35)

式中：

進(jìn)一步，聯(lián)立N個(gè)測量點(diǎn)的磁梯度信號接收模型，可得

F·M=ΔB，

(36)

式中：F=[F1;F2; …;FN]；ΔB=[B1,1-B1,2;B2,1-B2,2; …;BN,1-BN,2].

顯然，(36)式中的方程個(gè)數(shù)為3N，G的列數(shù)為3.磁性目標(biāo)參數(shù)估計(jì)問題可表示為

(37)

式中：U為估計(jì)向量。

由于載體做直線運(yùn)動且各測量點(diǎn)在X軸方向的相對坐標(biāo)為已知參數(shù)，故U=[YCPA,ZCPA,MX,MY,MZ]。當(dāng)N>1時(shí)測量模型為超定方程組，利用最小二乘法求解(36)式可得

M=F+ΔB，

(38)

式中：F+為F的M-P廣義逆矩陣，由于F是列滿秩矩陣，F(xiàn)+存在且唯一，

F+=(FTF)-1FT.

(39)

(39)式代入(37)式，磁性目標(biāo)參數(shù)估計(jì)問題可轉(zhuǎn)換為

(40)

式中：I為單位矩陣；U=[YCPA,ZCPA].(40)式是一個(gè)典型的多變量優(yōu)化問題，相比于(37)式表述的優(yōu)化模型，其估計(jì)的向量維數(shù)由5維降到了2維，極大地降低了算法的運(yùn)算量。

粒子群優(yōu)化[19](PSO)算法是一種群體智能優(yōu)化算法，算法中的每個(gè)粒子都代表一個(gè)問題的潛在解，粒子速度決定了粒子移動的方向和距離，速度隨自身和其他粒子的移動經(jīng)驗(yàn)進(jìn)行動態(tài)調(diào)整，從而實(shí)現(xiàn)個(gè)體在可解空間內(nèi)的動態(tài)尋優(yōu)，故采用PSO算法能夠解算(40)式表述的多變量優(yōu)化問題。估計(jì)向量U代入(38)式可以得到目標(biāo)的磁矩參數(shù)，進(jìn)一步將位置參數(shù)和磁矩參數(shù)代入偶極子模型，可得到磁性目標(biāo)的擬合信號，而且磁擬合信號的能量就可以作為區(qū)分目標(biāo)信號和噪聲的檢測量。

3 仿真實(shí)驗(yàn)

設(shè)磁性目標(biāo)的中心為空間直角坐標(biāo)系的原點(diǎn)，X軸與平臺的軸線重合并指向平臺的艏部，Y軸指向左舷，Z軸豎直向上。仿真條件設(shè)置為：初始時(shí)刻目標(biāo)的位置為(-200 m，-50 m，5 m)；目標(biāo)的磁矩為(-3×103A·m2,1×103A·m2,5×103A·m2)，磁梯度張量系統(tǒng)測量基線的尺寸為0.5 m；空間采樣長度為0.1 m；觀測窗區(qū)間為[-2R0,2R0]。下面分5個(gè)層面對本文算法的性能進(jìn)行仿真實(shí)驗(yàn)驗(yàn)證。

3.1 與傳統(tǒng)檢測方法性能的對比

圖5(a)和圖5(b)所示分別為信噪比為6 dB和-6 dB時(shí)，磁標(biāo)量OBF檢測器、傳統(tǒng)的特征值檢測方法和MED檢測方法的檢測性能對比，此處的信噪比定義為

(41)

式中：Amax表示信號的峰值。OBF檢測器和特征值檢測器的凸出峰值越突出，表示檢測的區(qū)分度越大；MED檢測器向下的峰值越明顯，表示此時(shí)的信息熵越小，即檢測的區(qū)分度越大。在圖5(a)和圖5(b)中，可以看出在高信噪比條件下，3種檢測算法都能準(zhǔn)確檢測目標(biāo)，但OBF檢測性能更加穩(wěn)定而且具有更高的檢測信噪比；在低信噪比條件下，特征值檢測方法與MED檢測方法皆已失效，此時(shí)OBF檢測器可以準(zhǔn)確檢測目標(biāo)，而且仍然具有較高的檢測信噪比。

圖5 標(biāo)量OBF檢測器與傳統(tǒng)檢測方法性能的對比

3.2 空間采樣長度對OBF檢測器的影響

圖6(a)和圖6(b)所示分別為信噪比為-6 dB和-12 dB時(shí)，空間采樣長度為0.1 m、0.2 m以及0.5 m時(shí)所對應(yīng)的磁標(biāo)量OBF檢測器性能。從圖6中可以看出，在同等信噪比條件下，空間采樣長度越小，則檢測的信噪比越高而且檢測的波動性越小，此結(jié)果與1.2節(jié)中“空間采樣長度越小，則信號與噪聲區(qū)分越明顯以及檢測穩(wěn)定性越好”的理論證明是一致的。

圖6 空間采樣長度對標(biāo)量OBF檢測器性能的影響

3.3 矢量與標(biāo)量OBF檢測器對比

圖7(a)和圖7(b)所示為空間采樣長度為0.1 m、信噪比分別為-6 dB和-12 dB時(shí)的標(biāo)量OBF和矢量OBF檢測器性能比較，從中可以看出在沒有目標(biāo)時(shí)矢量OBF的檢測曲線的波動性較小，具有更穩(wěn)定的檢測性能，此結(jié)果與1.2節(jié)中“矢量OBF檢測器相對于標(biāo)量OBF檢測器的優(yōu)勢是改善了標(biāo)量檢測的穩(wěn)定性”的理論證明是一致的。

圖7 標(biāo)量和矢量檢測器性能對比(空間采樣長度0.1 m)

然而受地磁場的影響，當(dāng)矢量磁傳感器發(fā)生輕微晃動時(shí)，其輸出就會產(chǎn)生較大的變化，晃動幅值為0°、0.01°和0.1°時(shí)的矢量OBF檢測器性能如圖8所示，輕微的晃動都會導(dǎo)致檢測器性能發(fā)生惡化，當(dāng)晃動幅值達(dá)到0.1°時(shí)，矢量OBF檢測算法已經(jīng)失效。

圖8 矢量磁力儀晃動對矢量OBF檢測器性能的影響

3.4 磁梯度OBF檢測器的抗晃動性

圖9(a)、圖9(b)、圖9(c)、圖9(d)分別為信噪比為-6 dB和-12 dB時(shí)，矢量磁梯度OBF檢測器與標(biāo)量磁梯度OBF檢測器的性能比較，從中可以看出在同等信噪比和空間采樣長度條件下，當(dāng)沒有目標(biāo)時(shí)，矢量梯度檢測相對于標(biāo)量梯度檢測的波動要小得多，即矢量梯度檢測具有更低的虛警率，此結(jié)果符合1.2節(jié)中“矢量OBF檢測器相對于標(biāo)量OBF檢測器的優(yōu)勢是改善了標(biāo)量檢測的穩(wěn)定性”的結(jié)論。

圖9 標(biāo)量和矢量磁梯度OBF檢測器性能對比

圖10所示為-6 dB信噪比時(shí)不同晃動角度對應(yīng)的磁梯度矢量OBF檢測曲線，從中可以看出在10°晃動角以下，算法的性能沒有發(fā)生明顯的變化，仍然保持較高的檢測信噪比。通過與磁標(biāo)量OBF檢測器、磁矢量OBF檢測器以及磁梯度標(biāo)量檢測器進(jìn)行對比，可以發(fā)現(xiàn)矢量磁梯度OBF檢測器既可以探測到更豐富的目標(biāo)磁信息，同樣也具有較好的抗晃動性能。

3.5 目標(biāo)磁矩和位置參數(shù)估計(jì)

在信噪比為-12 dB、晃動幅度為1°時(shí)，PSO算法優(yōu)化前后的檢測性能對比如圖11所示。由圖11可以看出：PSO算法優(yōu)化后的檢測具有更高的信噪比；PSO優(yōu)化前的峰值持續(xù)約為50 m，而PSO優(yōu)化后其檢測峰值更為尖銳，與真實(shí)CPA點(diǎn)在X軸方向的距離僅為3 m，意味著PSO優(yōu)化后的檢測器能夠給出更為準(zhǔn)確的目標(biāo)范圍。

圖11 PSO優(yōu)化前后檢測性能對比(SNR=-12 dB，晃動幅度1°)

表1所示為基于矢量磁梯度OBF檢測器的反演信號，利用PSO算法估計(jì)得到的目標(biāo)位置參數(shù)與磁矩參數(shù)，可以看出該算法在-12 dB信噪比條件下，對目標(biāo)的位置參數(shù)以及磁矩參數(shù)的估計(jì)誤差皆小于20%.

表1 位置和磁矩參數(shù)估計(jì)結(jié)果

綜上所述，磁梯度OBF檢測方法可以克服平臺晃動的影響以及提供豐富的目標(biāo)磁場信息，與PSO優(yōu)化算法相結(jié)合，不但能夠提升進(jìn)一步檢測的信噪比，而且可以在低信噪比條件下準(zhǔn)確估計(jì)目標(biāo)的位置參數(shù)和磁矩參數(shù)，可應(yīng)用于遠(yuǎn)程探測磁性目標(biāo)與確定目標(biāo)的搜索范圍。

4 實(shí)測試驗(yàn)

圖12所示為午夜時(shí)分在某消磁實(shí)驗(yàn)室進(jìn)行的探測試驗(yàn)，磁梯度測量系統(tǒng)為4個(gè)傳感器呈十字型布置的Mag-03磁通門傳感器，采集設(shè)備為HIOKI高精度24位AD采集系統(tǒng)，受試驗(yàn)條件限制，探測的目標(biāo)為一個(gè)運(yùn)動的制式滅火器(幾何尺寸與一般的未爆彈相近)，磁梯度測量的噪聲標(biāo)準(zhǔn)差約為0.2 nT/m.

圖12 高精度磁梯度測量

圖13所示為HIOKI高精度采集系統(tǒng)采集的磁性目標(biāo)磁梯度信號，其中GY,X、GY,Y和GY,Z分別表示矢量磁信號在Y軸方向上的梯度，磁梯度系統(tǒng)與磁力儀目標(biāo)的CPA距離為1.2 m。圖13中可以明顯觀測到目標(biāo)的通過特性曲線，不適合檢驗(yàn)算法在低信噪比條件下的檢測性能。

圖13 HIOKI高精度采集系統(tǒng)測量的磁梯度信號

為進(jìn)一步檢驗(yàn)算法的性能，如圖14所示，中午時(shí)分在某院的小樹林進(jìn)行磁性目標(biāo)磁梯度探測試驗(yàn)，AD采集設(shè)備為16位NI USB-6216采集卡，共進(jìn)行了6航次測量試驗(yàn)，磁梯度測量的噪聲標(biāo)準(zhǔn)差約為65 nT/m.

如圖15(a)、圖15(b)、圖15(c)所示，分別是采樣率為5 Hz、50 Hz以及500 Hz時(shí)，6個(gè)航次的磁梯度OBF檢測結(jié)果，從中可以看出采樣頻率越高則檢測器的檢測能力越強(qiáng)，此結(jié)果驗(yàn)證了1.2節(jié)中“空間采樣長度越小，則信號與噪聲區(qū)分越明顯以及檢測穩(wěn)定性越好”的理論證明。

表2所示為利用PSO磁性目標(biāo)參數(shù)估計(jì)算法對目標(biāo)參數(shù)的估計(jì)結(jié)果，其中位置參數(shù)的相對估計(jì)誤差小于30%，受感應(yīng)磁場的影響，橫向目標(biāo)的磁矩模值在2 A·m2左右，垂向目標(biāo)的磁矩在0.3 A·m2左右。

表2 目標(biāo)參數(shù)估計(jì)結(jié)果

這里通過磁梯度傳播衰減理論來推算矢量磁梯度OBF檢測器對水雷目標(biāo)的探測距離(目標(biāo)磁矩在102A·m2量級)，對比兩次試驗(yàn)的磁梯度測量噪聲可知，第1次磁梯度測量噪聲為第2次測量噪聲的1/320，其探測距離同等條件下會增加至之前的4.2倍。由圖15可知，當(dāng)目標(biāo)磁矩在0.3 A·m2且測量噪聲標(biāo)準(zhǔn)差在65 nT/m左右時(shí)探測的距離大于2 m，因此當(dāng)目標(biāo)磁矩大于100 A·m2且測量噪聲標(biāo)準(zhǔn)差在0.2 nT/m左右時(shí)，磁梯度OBF檢測器的探測距離大于35 m.

圖15 不同采樣率下矢量磁梯度OBF檢測

5 結(jié)論

本文針對標(biāo)量OBF檢測器磁場解釋性弱的問題以及磁矢量OBF檢測器抗晃動性能差的缺點(diǎn)，提出了一種矢量磁梯度的OBF檢測器，并分析指出空間采樣長度是影響OBF檢測器性能的關(guān)鍵因素且矢量OBF檢測器相較于標(biāo)量OBF檢測器可以降低檢測的虛警概率。得出主要結(jié)論如下：

1)OBF檢測器相對于傳統(tǒng)的檢測算法具有更優(yōu)的檢測性能。

2)相對于傳統(tǒng)的OBF檢測器，矢量磁梯度OBF檢測器的優(yōu)勢在于兼具抗晃動性和磁場解釋性。

3)矢量磁梯度OBF檢測器與PSO目標(biāo)參數(shù)估計(jì)算法相結(jié)合，可以在低信噪比條件下比較準(zhǔn)確地估計(jì)磁性目標(biāo)的位置參數(shù)以及磁矩參數(shù)。

參考文獻(xiàn)(References)

[1] LIN J，WANG M C，ZHAO J.Review: progress in SQUID-based geophysical precision measurement technology[J].Journal of Harbin Institute of Technology(New Series)，2020，27(3):101-115.

[2] ALLEN G I，MATTHEWS R，WYNN M.Mitigation of platform generated magnetic noise impressed on a magnetic sensor mounted in an autonomous underwater vehicle[C]∥Proceedings of MTS/IEEE Oceans 2001.Honolulu，HI，US: IEEE，2001.

[3] ALLEN G I，SULZBERGER G，BONO J T，et al.Initial evaluation of the new real-time tracking gradiometer designed for small unmanned underwater vehicles[C]∥Proceedings of MTS/IEEE Oceans 2005.Washington，DC，US：IEEE，2005.

[4] GINZBURG B，F(xiàn)RUMKIS L，KAPLAN B Z.Processing of magnetic scalar gradiometer signals using orthonormalized functions[J].Sensors and Actuators A: Physical，2002，102(1/2):67-75.

[5] SHEINKER A，MOLDWIN M B.Magnetic anomaly detection(MAD)of ferromagnetic pipelines using principal component analysis(PCA)[J].Measurement Science & Technology，2016，27(4): 045104.

[6] SHEINKER A，SALOMONSKI N，GINZBURG B，et al.Magnetic anomaly detection using entropy filter[J].Measurement Science and Technology，2008，19(4): 045205.

[7] SHEINKER A，GINZBURG B，SALOMONSKI N，et al.Magnetic anomaly detection using high-order crossing method[J].IEEE Transactions on Geoscience & Remote Sensing，2012，50(4):1095-1103.

[8] SHEINKER A，SHKALIM A，SALOMONSKI N，et al.Processing of a scalar magnetometer signal contaminated by 1/fαnoise[J].Sensors and Actuators A: Physical，2007，138(1): 105-111.

[9] SCHMELZ M，STOLZ R，ZAKOSARENKO V，et al.Field-stable SQUID magnetometer with sub-fT Hz-1/2resolution based on sub-micrometer cross-type Josephson tunnel junctions[J].Superconductor Science & Technology，2011，24(6):065009.

[10] SCHMELZ M，ZAKOSARENKO V，CHWALA A，et al.Thin-film based ultralow noise SQUID magnetometer[J].IEEE Transactions on Applied Superconductivity，2016，26(5): 1600804.

[11] 張光，張英堂，李志寧，等.載體平動條件下的磁梯度張量定位方法[J].華中科技大學(xué)學(xué)報(bào)(自然科學(xué)版)，2013，41(1): 21-24.

ZHANG G，ZHANG Y T，LI Z N，et al.Localizing method of magnetic field gradient tensor under carriers moving parallelly[J].Journal of Huazhong University of Science&Technology(Natural Science Edition)，2013，41(1):21-24.(in Chinese)

[12] 于振濤，呂俊偉，許素琴，等.運(yùn)動平臺的磁性目標(biāo)實(shí)時(shí)定位方法[J].哈爾濱工程大學(xué)學(xué)報(bào)，2015，36(5): 606-610.

YU Z T，Lü J W，XU S Q，et al.A real-time localization method of a magnetic target based on moving flat[J].Journal of Harbin Engineering University，2015，36(5):606-610.(in Chinese)

[13] 遲鋮，任建存，呂俊偉，等.基于磁梯度張量的目標(biāo)多測量點(diǎn)線性定位方法[J].探測與控制學(xué)報(bào)，2017，39(5): 58-62，70.

CHI C，REN J C，Lü J W，et al.Linear localization method based on magnetic gradient tensor of muliti-points[J].Journal of Detection & Control，2017，39(5):58-62，70.(in Chinese)

[14] CHWALA A，SCHMELZ M，ZAKOSARENKO V，et al.Underwater operation of a full tensor SQUID gradiometer system[J].Superconductor Science and Technology，2019，32(2): 024003.

[15] 鐘煬，管彥武，石甲強(qiáng)，等.利用IGRF模型計(jì)算全張量地磁梯度[J].物探與化探，2020，44(3):582-590.

ZHONG Y，GUAN Y W，SHI J Q，et al.The calculation method of full tensor geomagnetic gradient based on IGRF model[J].Geophysical and Geochemical Exploration，2020，44(3):582-590.(in Chinese)

[16] GINZBURG B，F(xiàn)RUMKIS L，KAPLAN B Z.An efficient method for processing scalar magnetic gradiometer signals[J].Sensors and Actuators A: Physical，2004，114(1):73-79.

[17] FRUMKIS L，GINZBURG B，SALOMONSKI N，et al.Optimization of scalar magnetic gradiometer signal processing[J].Sensors and Actuators A: Physical，2005，121(1):88-94.

[18] 張朝陽，肖昌漢，高俊吉，等.磁性物體磁偶極子模型適用性的試驗(yàn)研究[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào)，2010，18(5):862-868.

ZHANG C Y，XIAO C H，GAO J J，et al.Experiment research of magnetic dipole model applicability for a magnetic object[J].Journal of Basic Science and Engineering，2010，18(5): 862-868.(in Chinese)

[19] KENNEDY J，EBERHART R.Particle swarm optimization[C]∥Proceedings of the 4th International Conference on Neural Networks.Perth，WA，Australia：IEEE，1995：1942-1948.