吳 芳，王學(xué)敏，劉首善，王潁航，張欽程

（1.海軍航空大學(xué)，山東 煙臺(tái) 264001；2.92697 部隊(duì)，海南 陵水 572400；3.91445 部隊(duì)，遼寧 大連 116043）

邊界元法是一種在半空間區(qū)域基于格林函數(shù)，通過標(biāo)量磁位或矢量磁位分布求拉普拉斯的邊值問題[1]。根據(jù)電磁場(chǎng)的基本原理，可將計(jì)算潛艇高空磁場(chǎng)分布的問題轉(zhuǎn)化為求解空間磁位分布的問題。通過計(jì)算標(biāo)量磁位的負(fù)梯度或是計(jì)算矢量磁位的旋度[2-3]，可以得到潛艇目標(biāo)磁源在空間中的磁場(chǎng)分布情況，該方法具有較高的精度。

反潛機(jī)使用磁探儀對(duì)水下目標(biāo)進(jìn)行探測(cè)時(shí)，由于距目標(biāo)較遠(yuǎn)，磁探儀接收到的磁異常信號(hào)能量將直接影響設(shè)備的探測(cè)效能。文獻(xiàn)[4]根據(jù)邊界元法的基本原理，建立了潛艇磁場(chǎng)的高空預(yù)測(cè)模型，并對(duì)該模型進(jìn)行了理論驗(yàn)證。本文基于潛艇磁場(chǎng)的高空預(yù)測(cè)模型，結(jié)合航空磁探儀實(shí)際搜潛過程，仿真分析潛艇不同運(yùn)動(dòng)航向，以及反潛機(jī)不同搜索高度下，潛艇磁場(chǎng)的空間特征。

1 潛艇磁場(chǎng)組成

地球是一個(gè)巨大的磁場(chǎng)，潛艇上的大部分設(shè)備與材料都是良好的磁化材料，所以潛艇在磁場(chǎng)中極易被磁化，并在其周圍空間產(chǎn)生方向不同、強(qiáng)度各異的磁場(chǎng)。潛艇磁化后的磁場(chǎng)可以分為永久磁場(chǎng)、感應(yīng)磁場(chǎng)及變化磁場(chǎng)3 個(gè)部分[3]。

1.1 永久磁場(chǎng)

潛艇一般都是由高強(qiáng)度的合金鋼制成，屬于鐵磁性材料。在建造潛艇的時(shí)候，要經(jīng)過電焊和鉚接等工藝，會(huì)產(chǎn)生磁場(chǎng)。鐵合金含有“磁疇”，每一個(gè)磁疇都是一個(gè)小的磁體，有南北極。當(dāng)磁疇無規(guī)則排列的時(shí)候，會(huì)產(chǎn)生一個(gè)很小的磁力線圈圖。當(dāng)把這種含有磁疇的合金放在穩(wěn)定的磁場(chǎng)中，它的分子會(huì)被激勵(lì)，在被錘擊和受熱后，磁疇自己趨于定向，形成南北極，使所有磁疇的磁場(chǎng)都得到加強(qiáng)。地球的磁場(chǎng)雖然不是很強(qiáng)，但是潛艇是一個(gè)很大的鋼結(jié)構(gòu)體，不可避免地要在地磁場(chǎng)的作用下加工和建造，在這個(gè)過程中，材料內(nèi)部應(yīng)力的反復(fù)變化、溫度的升降變遷及局部磁場(chǎng)的影響，都會(huì)引起鐵磁材料內(nèi)無磁滯磁化的形成，在潛艇建成形成潛艇的永久磁場(chǎng)。由于潛艇的結(jié)構(gòu)很復(fù)雜及外形曲面的不規(guī)則性，現(xiàn)代消磁技術(shù)并不能完全消除掉潛艇的永久磁場(chǎng)，特別是潛艇的縱向永久磁場(chǎng)和垂直方向永久磁場(chǎng)還無法分離。潛艇的磁場(chǎng)有3 個(gè)主要的分量：垂直分量、縱向分量、橫向分量。這3 個(gè)部分合起來構(gòu)成了潛艇完整的永久磁場(chǎng)。

1.2 感應(yīng)磁場(chǎng)

潛艇服役后在地磁場(chǎng)的作用下工作，潛艇各部分的鐵磁材料內(nèi)部會(huì)引起可逆的磁化過程，從而產(chǎn)生潛艇的感應(yīng)磁特性。潛艇在水下航行時(shí)，受波浪和爆炸的沖擊，或者受高速水擊作用等振動(dòng)后，也會(huì)產(chǎn)生感應(yīng)磁場(chǎng)。由于感應(yīng)磁場(chǎng)的存在，潛艇在某一緯度海域長(zhǎng)期活動(dòng)，潛艇的固定磁場(chǎng)會(huì)慢慢接近某一固定值。但是，當(dāng)潛艇轉(zhuǎn)移到另一個(gè)磁緯度區(qū)航行一段時(shí)間后，潛艇的固定磁場(chǎng)就會(huì)再次變化，并逐漸地穩(wěn)定到另一個(gè)相應(yīng)的穩(wěn)定值上。

感應(yīng)磁場(chǎng)的大小主要與潛艇所在磁緯度的地磁場(chǎng)的大小、潛艇的航向、潛艇本身搖擺的變化、潛艇所用鋼材的磁性能和潛艇的形狀等因素有關(guān)。目前的常規(guī)潛艇受內(nèi)部空間及耗能的限制，大多數(shù)未裝消磁系統(tǒng)，因此感應(yīng)磁性無法及時(shí)消除。潛艇的感應(yīng)磁性與永久磁性共同作用，引起潛艇所在位置周圍地磁場(chǎng)的異常。

1.3 變化磁場(chǎng)

潛艇的感應(yīng)磁場(chǎng)和固定磁場(chǎng)由于地磁場(chǎng)磁化而產(chǎn)生，稱為磁化靜磁場(chǎng)，潛艇的渦流磁場(chǎng)、雜散磁場(chǎng)、低頻磁場(chǎng)是電流引起的磁場(chǎng)，稱為變化磁場(chǎng)。

潛艇磁場(chǎng)就是由于鐵磁性潛艇受地磁場(chǎng)磁化而產(chǎn)生的附加磁場(chǎng)。它將引起周圍空間磁場(chǎng)分布的改變。航空磁探儀通過探測(cè)潛艇目標(biāo)的磁異常信號(hào)來發(fā)現(xiàn)目標(biāo)。因此，潛艇磁場(chǎng)是航空磁探儀搜索目標(biāo)的信息源。由于磁場(chǎng)受海水、空氣、泥沙等介質(zhì)的影響小，使得基于磁異常信號(hào)的航空磁探測(cè)具有可靠性強(qiáng)、隱蔽性好、定位精度高、搜索連續(xù)和反應(yīng)迅速等特點(diǎn)，成為水下磁性目標(biāo)探測(cè)的重要手段。為提高航空磁探儀的檢測(cè)能力，進(jìn)一步對(duì)潛艇目標(biāo)進(jìn)行定位和識(shí)別，需要分析潛艇目標(biāo)的高空磁場(chǎng)分布特征，為航空磁探儀提供目標(biāo)信號(hào)參考。

2 基于邊界元法的潛艇高空磁場(chǎng)預(yù)測(cè)模型

由邊界元法原理，可得到P 點(diǎn)處潛艇的磁感應(yīng)強(qiáng)度為

式中：n 為包絡(luò)觀測(cè)面S 上的單位外法向量；BQ為包絡(luò)觀測(cè)面S 上Q（x0，y0，z0）處的磁感應(yīng)強(qiáng)度。

為得到實(shí)際條件下適用的潛艇磁場(chǎng)預(yù)測(cè)模型，對(duì)式（1）進(jìn)行離散化，則

式中：Mi為第i 個(gè)平面剖分的邊界單元數(shù)，i=1，2，…，6；nij為第ij 個(gè)平面的單位法向量；BQij為單元ij 的測(cè)量磁感應(yīng)強(qiáng)度數(shù)據(jù)，j=1，2，…，Mi；Sij為測(cè)量單元ij的面積。

潛艇磁場(chǎng)預(yù)測(cè)模型的包絡(luò)觀測(cè)面示意圖如圖1 所示。假設(shè)采用長(zhǎng)方體包絡(luò)面包圍潛艇磁源進(jìn)行三分量磁場(chǎng)測(cè)量，邊界面的劃分及測(cè)量點(diǎn)的布放規(guī)劃如圖1所示。測(cè)量得到面積單元ij 的磁感應(yīng)強(qiáng)度數(shù)據(jù)BQij，并記錄各觀測(cè)點(diǎn)的坐標(biāo)信息及該面積元的單位法向量，根據(jù)式（2）可以進(jìn)行潛艇空間磁場(chǎng)預(yù)測(cè)。

圖1 預(yù)測(cè)模型的觀測(cè)面

3 潛艇高空磁場(chǎng)信號(hào)仿真分析

本節(jié)主要根據(jù)潛艇高空磁場(chǎng)預(yù)測(cè)模型，即式（2），仿真分析潛艇高空磁異常信號(hào)的空間特征。

仿真數(shù)據(jù)：①潛艇目標(biāo)，長(zhǎng)半軸為2.8 m，短半軸為0.2 m 的鐵磁性旋轉(zhuǎn)橢球體；②長(zhǎng)方體測(cè)量區(qū)域，長(zhǎng)為7.50 m，寬與高均為1.50 m 的長(zhǎng)方體表面；③邊界面劃分，將上下兩表面和左右兩側(cè)面各分成125 個(gè)長(zhǎng)方形面積元，將前后兩側(cè)面各分成25 個(gè)長(zhǎng)方形面積元，則整個(gè)邊界面共分成550 個(gè)長(zhǎng)方形面積元；④地磁總場(chǎng)約為50 000 nT，磁傾角為52°，磁偏角為-3°。

3.1 潛艇一定航深、不同運(yùn)動(dòng)航向條件下的高空磁場(chǎng)分布

通過測(cè)量潛艇縮比模型在東、西、南和北4 個(gè)典型航向上的包絡(luò)面磁場(chǎng)數(shù)據(jù)，根據(jù)預(yù)測(cè)模型得到潛艇目標(biāo)在4 個(gè)航向上的空間磁場(chǎng)分布，并投影到地磁場(chǎng)方向得到總場(chǎng)值，如圖2—圖5 所示。

圖3 東航向

圖4 南航向

圖5 西航向

由圖2—圖5 對(duì)比潛艇的4 航向高空300 m 的總場(chǎng)分布可以發(fā)現(xiàn)，潛艇4 航向總場(chǎng)空間分布形狀基本相同，均具有1 個(gè)正峰和1 個(gè)負(fù)峰。4 航向之間的區(qū)別在于正峰與負(fù)峰之間的位置及幅值大小。就峰值而言，北航向空間磁場(chǎng)最大，東航向次之，南航向空間磁場(chǎng)的峰值最小。北航向、東航向、西航向的正峰幅值均小于負(fù)峰幅值，南航向反之。

3.2 潛艇一定航深、不同高度條件下的磁場(chǎng)分布

潛艇上方100、300、500 和700 m 的磁場(chǎng)分布如圖6 所示。

圖6 總場(chǎng)分布隨距離的變化

由圖6 可知，對(duì)比潛艇高空100、300、500 和700 m總場(chǎng)分布不難發(fā)現(xiàn)，隨著空間探測(cè)點(diǎn)與潛艇之間距離的不斷增大，信號(hào)衰減快。隨著高度不斷增加，潛艇磁場(chǎng)的空間波形基本保持一致。當(dāng)磁探儀以同等速度及采樣頻率按相同的采樣軌跡從潛艇上方經(jīng)過時(shí)，磁探儀采樣獲得的信號(hào)寬度不斷增大。

4 結(jié)論

由邊界元法建立的潛艇高空磁場(chǎng)預(yù)測(cè)模型可延拓得到潛艇磁場(chǎng)的空間分布，基于潛艇高空預(yù)測(cè)磁場(chǎng)模型，本文仿真分析了潛艇在不同運(yùn)動(dòng)航向、不同高度的空間磁場(chǎng)分布特征。由仿真結(jié)果可知，不同運(yùn)動(dòng)航向條件下，潛艇總場(chǎng)空間分布形狀基本相同，均具有1 個(gè)正峰和1 個(gè)負(fù)峰；隨著空間探測(cè)點(diǎn)與潛艇之間距離的增大，信號(hào)衰減，但潛艇磁場(chǎng)的空間波形基本保持一致。各航向之間高空磁場(chǎng)分布區(qū)域性明顯，具有一定的分布特征，能夠?yàn)楹娇沾盘絻x的目標(biāo)探測(cè)和識(shí)別提供支持。通過仿真分析潛艇目標(biāo)的高空磁場(chǎng)分布特征，為提高航空磁探儀的檢測(cè)能力，以及進(jìn)一步對(duì)潛艇目標(biāo)進(jìn)行定位和識(shí)別提供理論支持。