何佳敏，齊紅麗，雷華明，吉小軍

(1.上海交通大學(xué)電子信息與電氣工程學(xué)院，上海 200240；2.中國船舶重工集團(tuán)公司第704研究所，上海 200031)

0 引言

在磁場測量領(lǐng)域，磁通門磁強(qiáng)計(jì)憑借其較高的靈敏度、精度和低成本、安全可靠等特點(diǎn)被廣泛應(yīng)用。三軸磁通門磁強(qiáng)計(jì)更是可以直接測得空間任意點(diǎn)的3個正交磁場分量，以獲得磁場的完整信息[1]。在使用磁強(qiáng)計(jì)對空間中的磁異常目標(biāo)進(jìn)行探測時，測得的磁場信號通常為地磁場信號與磁異常信號的疊加。一般地磁場大小范圍在40～60 μT，而磁異常信號大小根據(jù)磁強(qiáng)計(jì)與目標(biāo)之間距離的不同常小于1 μT，地磁場信號遠(yuǎn)大于磁異常信號[2]。此外，受制造工藝等因素影響，三軸磁強(qiáng)計(jì)還存在三軸非正交、靈敏度和零點(diǎn)等誤差，這些誤差使得磁強(qiáng)計(jì)在恒定磁場下以任意姿態(tài)旋轉(zhuǎn)時輸出的總磁場不恒定，引起所謂的“轉(zhuǎn)向誤差”。測量過程中磁強(qiáng)計(jì)自身的晃動和旋轉(zhuǎn)導(dǎo)致的轉(zhuǎn)向誤差也會疊加在測量輸出信號中，嚴(yán)重制約探測距離和探測精度。

為了減少地磁場等背景磁場對微弱磁異常信號探測造成的干擾，章志濤等提出了基于三端式磁通門的磁力梯度儀[3]，該梯度儀采用雙磁通門探頭的結(jié)構(gòu)形式，實(shí)現(xiàn)了對背景磁場和磁場梯度的閉環(huán)測量，但沒有校正磁通門本身存在的轉(zhuǎn)向差；高翔等提出了基于線性神經(jīng)網(wǎng)絡(luò)的轉(zhuǎn)向差校正算法[4-5]，建立了單個三軸磁強(qiáng)計(jì)轉(zhuǎn)向差的校正模型，并對構(gòu)成磁力梯度儀的兩個磁強(qiáng)計(jì)的非同軸誤差進(jìn)行了校正，但沒有考慮到兩個磁強(qiáng)計(jì)之間的靈敏度差異；黃玉等提出了基于函數(shù)鏈接型神經(jīng)網(wǎng)絡(luò)的校正算法[6-7]，構(gòu)建了上述幾種誤差的校正模型，具有一定的理論研究價值，但校正過程需要在上位機(jī)采集到傳感器的輸出數(shù)據(jù)后利用軟件進(jìn)行校正，實(shí)時性較差，且部分校正參數(shù)的獲得需要嚴(yán)格的實(shí)驗(yàn)環(huán)境，難以應(yīng)用于工程實(shí)際測量。

基于實(shí)際的工程應(yīng)用需求，本文提出了一種由間隔一定距離的兩個三軸磁通門磁強(qiáng)計(jì)(分別作為測量端和補(bǔ)償端)構(gòu)成的實(shí)時動態(tài)地磁補(bǔ)償系統(tǒng)[13]。工作過程中，系統(tǒng)將兩磁強(qiáng)計(jì)對應(yīng)敏感軸的輸出通過模擬電路進(jìn)行實(shí)時差分處理以減小地磁場影響，提高對磁異常信號處理的增益范圍，同時針對兩磁強(qiáng)計(jì)本身存在的轉(zhuǎn)向誤差和靈敏度差異，本文采用神經(jīng)網(wǎng)絡(luò)校正算法估計(jì)誤差參數(shù)，并通過模擬電路融合在地磁動態(tài)補(bǔ)償過程中完成，從而實(shí)現(xiàn)對微小磁異常目標(biāo)的遠(yuǎn)距離、實(shí)時、高精度檢測。

1 三軸磁強(qiáng)計(jì)轉(zhuǎn)向差校正模型建立

設(shè)理想三軸磁強(qiáng)計(jì)的測量坐標(biāo)系由3個相互正交的敏感軸OX0、OY0、OZ0構(gòu)成，受制造工藝影響，實(shí)際三軸磁強(qiáng)計(jì)的測量坐標(biāo)系由3個非正交敏感軸OX、OY、OZ構(gòu)成，二者關(guān)系如圖1所示。為便于建立校正模型，以O(shè)Z軸為基準(zhǔn)軸，YOZ平面為基準(zhǔn)平面，對實(shí)際三軸磁強(qiáng)計(jì)的測量坐標(biāo)系進(jìn)行校正。假設(shè)OZ0軸與OZ軸重合，OY軸在平面Y0OZ0內(nèi)，且與OY0軸夾角為α。OX軸與平面X0OZ0的夾角為β，OX軸在平面X0OZ0上的投影與OX0軸的夾角為γ。設(shè)理想磁通門磁強(qiáng)計(jì)的三軸輸出分別為BX0、BY0、BZ0，實(shí)際磁通門磁強(qiáng)計(jì)的三軸輸出分別BX、BY、BZ。

圖1 理想坐標(biāo)系與實(shí)際坐標(biāo)系的關(guān)系

由圖1可得：

(1)

考慮到實(shí)際磁強(qiáng)計(jì)的三軸敏感元件靈敏度以及激勵和信號提取電路特性參數(shù)不可能完全相同，導(dǎo)致三軸測量的靈敏度不同,輸出信號也存在零點(diǎn)誤差。為了便于分析，設(shè)磁強(qiáng)計(jì)OZ軸的靈敏度為1，OX軸和OY軸相對OZ軸的相對靈敏度為KX和KY，三軸的零點(diǎn)誤差分別為bX0、bY0、bZ0。磁強(qiáng)計(jì)實(shí)際輸出與理想輸出之間的關(guān)系可表示為[8-12]

(2)

即

B=K·S·B0+b0

(3)

(4)

由式(4)可知，要想對實(shí)測值B進(jìn)行校正得到校正值B′，需要求得K11、K12、K13、K22、K23、bX0、bY0、bZ08個誤差參數(shù)，它們在測試系統(tǒng)建立完成后表現(xiàn)為固定的系統(tǒng)誤差參數(shù)，可以利用實(shí)際測得的磁場數(shù)據(jù)，建立神經(jīng)網(wǎng)絡(luò)校正模型來得到其估計(jì)值。將式(4)展開可得：

(5)

(6)

(7)

可以認(rèn)為經(jīng)校正后的單個磁強(qiáng)計(jì)三軸輸出互相正交，靈敏度相同(同校正前的Z軸)且沒有零點(diǎn)誤差。

傳統(tǒng)的轉(zhuǎn)向差校正多采用軟件校正的方法，需要先將磁強(qiáng)計(jì)的實(shí)測信號B采集到上位機(jī)后再用軟件算法進(jìn)行校正，存在事后處理等問題，無法滿足現(xiàn)場實(shí)時動態(tài)檢測的技術(shù)要求。此外，在某些場合下，若上位機(jī)采集的信號不是磁強(qiáng)計(jì)的實(shí)測信號(如兩個磁強(qiáng)計(jì)對應(yīng)敏感軸輸出的差分信號)，則無法使用軟件對單個磁強(qiáng)計(jì)的轉(zhuǎn)向差進(jìn)行校正。為此，本文提出采用硬件電路方式實(shí)現(xiàn)轉(zhuǎn)向差校正，單個磁強(qiáng)計(jì)轉(zhuǎn)向差硬件校正電路結(jié)構(gòu)框圖如圖2所示，其中加法電路和放大電路均由運(yùn)算放大器構(gòu)成，三軸零點(diǎn)誤差bX0、bY0、bZ0為恒定值，在電路中對應(yīng)的恒定電壓由基準(zhǔn)電壓通過電阻網(wǎng)絡(luò)分壓得到，5個放大電路的放大倍數(shù)根據(jù)5個系統(tǒng)誤差參數(shù)K11、K12、K13、K22、K23的值進(jìn)行調(diào)節(jié)。

圖2 單磁強(qiáng)計(jì)轉(zhuǎn)向差硬件校正電路框圖

2 融合轉(zhuǎn)向差校正的地磁補(bǔ)償硬件系統(tǒng)構(gòu)建

地磁是影響磁強(qiáng)計(jì)測量精度的另一重要因素，因此需要對地磁進(jìn)行動態(tài)補(bǔ)償。本文采用的地磁補(bǔ)償系統(tǒng)由相距一定距離安裝在剛性件上的2個三軸磁通門磁強(qiáng)計(jì)組成，分別作為測量端和補(bǔ)償端。由于地磁場在一定區(qū)域內(nèi)是均勻的，故可以認(rèn)為系統(tǒng)在工作過程中測量端和補(bǔ)償端磁強(qiáng)計(jì)所處位置的地磁場大小相同而目標(biāo)磁異常信號不同(磁場梯度)。因此將測量端與補(bǔ)償端磁強(qiáng)計(jì)對應(yīng)敏感軸的輸出在模擬電路部分進(jìn)行差分，即可消除地磁信號，保留磁異常信號，在后續(xù)測量電路中對有效的磁異常信號進(jìn)行增益調(diào)整至模數(shù)轉(zhuǎn)換器的參考電壓范圍，這樣可以極大地減小地磁信號和噪聲對磁異常信號的干擾，實(shí)現(xiàn)對微小磁異常目標(biāo)遠(yuǎn)距離和實(shí)時高精度檢測。由于上位機(jī)最終采集到的是2個磁強(qiáng)計(jì)對應(yīng)敏感軸輸出的差分信號，故需要在信號差分之前采用第一節(jié)所設(shè)計(jì)的硬件轉(zhuǎn)向差校正電路分別對兩個磁強(qiáng)計(jì)的三軸輸出進(jìn)行校正。

在地磁補(bǔ)償系統(tǒng)中，兩個磁強(qiáng)計(jì)對應(yīng)敏感軸的不同向或靈敏度不同也會影響地磁補(bǔ)償?shù)男ЧＴ趦蓚€磁強(qiáng)計(jì)分別進(jìn)行轉(zhuǎn)向差校正后，可認(rèn)為單個磁強(qiáng)計(jì)的三軸輸出互相垂直且靈敏度相同，但兩個測量坐標(biāo)系之間存在不同軸和靈敏度的差異，因此實(shí)際應(yīng)用中還應(yīng)做進(jìn)一步的聯(lián)合補(bǔ)償。

(8)

(9)

(10)

地磁補(bǔ)償系統(tǒng)的硬件實(shí)現(xiàn)結(jié)構(gòu)如圖3所示。首先2個磁強(qiáng)計(jì)的三軸原始輸出均需要經(jīng)過硬件校正電路進(jìn)行轉(zhuǎn)向差校正，同時測量端磁強(qiáng)計(jì)的三軸校正輸出還需要經(jīng)過由運(yùn)算放大器構(gòu)成的放大電路以統(tǒng)一兩個磁強(qiáng)計(jì)的靈敏度，放大電路的放大倍數(shù)為相對靈敏度K2。最后將處理后的6路磁強(qiáng)計(jì)信號輸入到差分放大電路，使得信號采集電路采集到的是差分信號ΔBX、ΔBY和ΔBZ，從而實(shí)現(xiàn)式(6)所示的補(bǔ)償原理。

圖3 地磁補(bǔ)償系統(tǒng)硬件實(shí)現(xiàn)結(jié)構(gòu)框圖

3 實(shí)驗(yàn)驗(yàn)證

根據(jù)上述技術(shù)路線和系統(tǒng)實(shí)現(xiàn)的結(jié)構(gòu)原理框圖，本節(jié)采用硬件方式實(shí)現(xiàn)了對單個磁強(qiáng)計(jì)的轉(zhuǎn)向差自校正和地磁補(bǔ)償系統(tǒng)中2個磁強(qiáng)計(jì)的靈敏度互校正以及動態(tài)地磁補(bǔ)償，并通過實(shí)驗(yàn)對比了軟件和硬件2種轉(zhuǎn)向差校正方法的校正效果以及地磁補(bǔ)償系統(tǒng)在有無校正電路接入情況下的地磁補(bǔ)償效果，最后比較了單磁強(qiáng)計(jì)和地磁補(bǔ)償系統(tǒng)對磁異常信號的探測效果。

3.1 單個磁強(qiáng)計(jì)軟件和硬件轉(zhuǎn)向差校正效果對比

本實(shí)驗(yàn)的實(shí)驗(yàn)場所選擇在作者所在單位某空曠場地，遠(yuǎn)離道路、電線、建筑等可能會使磁場產(chǎn)生劇烈波動的干擾源，可認(rèn)為在測試過程中，磁強(qiáng)計(jì)所處位置的磁場是恒定的。將三軸磁通門磁強(qiáng)計(jì)通過電纜連接到接口電路，接口電路驅(qū)動磁通門磁強(qiáng)計(jì)，同時將磁強(qiáng)計(jì)的輸出轉(zhuǎn)化為和待測磁場大小成正比的直流電壓，最后通過串口將測量值傳輸?shù)缴衔粰C(jī)進(jìn)行保存[14]。將磁強(qiáng)計(jì)放置在某一特定位置以任意姿態(tài)旋轉(zhuǎn)，同時上位機(jī)記錄輸出，記錄時間約30 s，記錄300組數(shù)據(jù)。

若在磁場恒定的條件下以任意姿態(tài)旋轉(zhuǎn)三軸磁強(qiáng)計(jì)，理想三軸磁強(qiáng)計(jì)輸出的總磁場值是不變的，等于理論總磁場值。但由于誤差的存在，三軸磁強(qiáng)計(jì)實(shí)際輸出的總磁場值隨姿態(tài)的變化在理論總磁場值上下波動。因此，可以利用實(shí)際測得的總磁場數(shù)據(jù)，建立神經(jīng)網(wǎng)絡(luò)校正模型，得到第一節(jié)所述8個誤差參數(shù)的估計(jì)值。

將用于構(gòu)建地磁補(bǔ)償系統(tǒng)的2個磁強(qiáng)計(jì)分別進(jìn)行上述實(shí)驗(yàn)，得到的誤差參數(shù)估計(jì)值如表1所示，并通過軟件對測得數(shù)據(jù)進(jìn)行轉(zhuǎn)向差校正。得到校正前后總磁場波動幅值對比如表2、圖4和圖5所示。可以看出，經(jīng)過軟件轉(zhuǎn)向差校正后，每個磁強(qiáng)計(jì)測得的總磁場最大波動幅值大幅度減小，證明了該校正算法的有效性，其誤差參數(shù)的估計(jì)值可以作為硬件校正電路參數(shù)的設(shè)計(jì)依據(jù)。

表1 轉(zhuǎn)向差校正參數(shù)估計(jì)值

表2 單磁強(qiáng)計(jì)轉(zhuǎn)向差軟件校正和硬件校正效果對比

圖4 測量端磁強(qiáng)計(jì)軟件轉(zhuǎn)向差校正前后總磁場波動

圖5 補(bǔ)償端磁強(qiáng)計(jì)軟件轉(zhuǎn)向差校正前后總磁場波動

圖6 測量端磁強(qiáng)計(jì)硬件轉(zhuǎn)向差校正前后總磁場波動

圖7 補(bǔ)償端磁強(qiáng)計(jì)硬件轉(zhuǎn)向差校正前后總磁場波動

(1)8個誤差參數(shù)的估計(jì)值與真實(shí)值之間存在差異；

(2)磁強(qiáng)計(jì)在以任意姿態(tài)旋轉(zhuǎn)時，所處空間中的磁場大小存在微小變化；

(3)受實(shí)驗(yàn)條件的限制，對硬件校正電路參數(shù)的調(diào)節(jié)存在誤差。

3.2 融合轉(zhuǎn)向差校正的地磁補(bǔ)償硬件實(shí)現(xiàn)

分別將測量端和補(bǔ)償端磁強(qiáng)計(jì)靜止放在相同的均勻恒定磁場下，在有硬件轉(zhuǎn)向差校正的情況下計(jì)算2個磁強(qiáng)計(jì)輸出的總磁場大小比值，認(rèn)為該比值是兩磁強(qiáng)計(jì)的相對靈敏度K2，實(shí)驗(yàn)測得K2=1.164?；谑?6)和圖3調(diào)節(jié)硬件電路參數(shù)，令差分放大電路的放大倍數(shù)A=1。將測量端和補(bǔ)償端磁強(qiáng)計(jì)相距30 cm固定在無磁鋁制安裝構(gòu)件上，如圖8所示，構(gòu)成地磁補(bǔ)償系統(tǒng)，并在戶外場地進(jìn)行測試，認(rèn)為該測試場所無明顯干擾源且磁場梯度為零。

圖8 實(shí)時動態(tài)地磁補(bǔ)償系統(tǒng)

然后實(shí)驗(yàn)測試在采用和不采用硬件校正電路2種情況下進(jìn)行，記錄系統(tǒng)在2種不同姿態(tài)下的三軸差分輸出，結(jié)果如表3所示?？梢钥闯龅卮叛a(bǔ)償系統(tǒng)在每個磁強(qiáng)計(jì)轉(zhuǎn)向差硬件校正和靈敏度互校正后，其每一軸的差分輸出較校正之前有明顯減小，校正后的地磁補(bǔ)償效果更為明顯。導(dǎo)致補(bǔ)償后輸出不為零的原因可能有：

表3 地磁補(bǔ)償系統(tǒng)硬件校正前后輸出對比 nT

(1)受實(shí)驗(yàn)條件的限制，校正電路參數(shù)的調(diào)節(jié)存在誤差；

(2)實(shí)驗(yàn)所處的場地磁場并非完全均勻；

(3)地磁補(bǔ)償系統(tǒng)兩磁強(qiáng)計(jì)安裝存在一定的同軸度誤差。

3.3 微小磁異常探測實(shí)驗(yàn)

將地磁補(bǔ)償系統(tǒng)靜止放在3.2節(jié)所述實(shí)驗(yàn)環(huán)境下,實(shí)驗(yàn)者攜帶磁性物體從測量端磁強(qiáng)計(jì)前方約3 m處經(jīng)過2次，行走路線垂直于無磁鋁制安裝構(gòu)件。圖9(a)為單個磁強(qiáng)計(jì)在未進(jìn)行地磁補(bǔ)償時的測試結(jié)果，可以看出三軸輸出均無明顯變化，無法識別由磁性物體所引起的磁異常信號。

采用地磁補(bǔ)償系統(tǒng)后，調(diào)節(jié)差分放大環(huán)節(jié)的增益系數(shù)，并進(jìn)行上述相同實(shí)驗(yàn)，得到三軸輸出結(jié)果如圖9(b)所示?？梢钥闯龅卮叛a(bǔ)償系統(tǒng)的輸出可以明確地反映出磁性物體2次經(jīng)過時產(chǎn)生的磁異常信號。

(a)未做地磁補(bǔ)償?shù)拇女惓Ｌ綔y結(jié)果

4 結(jié)論

本文提出了一種融合轉(zhuǎn)向差校正的實(shí)時動態(tài)地磁補(bǔ)償系統(tǒng)，該系統(tǒng)由2個三軸磁通門磁強(qiáng)計(jì)和硬件校正電路組成。通過實(shí)驗(yàn)可以得知，硬件轉(zhuǎn)向差校正效果略差于軟件轉(zhuǎn)向差校正，但依然有較好的校正效果，且上位機(jī)在采集到磁強(qiáng)計(jì)的輸出信號之后不再需要進(jìn)行后續(xù)處理，可以滿足工程上實(shí)時、在線、動態(tài)的檢測需求。融合轉(zhuǎn)向差校正和靈敏度互校正的地磁補(bǔ)償系統(tǒng)有更明顯的地磁補(bǔ)償效果，且該系統(tǒng)在進(jìn)行磁異常探測和定位的過程中能顯著減小地磁場對探測結(jié)果的影響，同時可以為微小的磁異常信號提供更大的增益范圍，能夠?qū)崟r直觀地從系統(tǒng)的輸出判斷磁異常目標(biāo)的有無，提高磁異常目標(biāo)探測的準(zhǔn)確度和成功率。