摘" 要：

為改進(jìn)水下航行器磁異常探測技術(shù)中磁異常實(shí)測數(shù)據(jù)難獲得、現(xiàn)有磁體模型計(jì)算精度不高等問題，基于有限元仿真探討建立無需區(qū)分近、遠(yuǎn)場的高精度混合磁體模型方法。應(yīng)用有限元數(shù)值方法仿真復(fù)雜結(jié)構(gòu)水下航行器的空間磁異常分布，以磁場數(shù)值解為收斂目標(biāo)，建立以偶極子陣元數(shù)目、位置及磁矩為參數(shù)的均勻磁化橢球體與偶極子陣列混合的航行器磁異常解析模型，采用非線性最小二乘算法求解模型系數(shù)。仿真結(jié)果表明，基于該模型得到磁異常計(jì)算值與全空間內(nèi)數(shù)值解的擬合度較高，測試平面平均誤差為3%。該模型在磁場延拓、高精度建模等方面可以進(jìn)一步應(yīng)用。關(guān)鍵詞：

航空磁異常探測; 有限元仿真; 磁體模擬; 數(shù)值擬合

中圖分類號：

O 411.3

文獻(xiàn)標(biāo)志碼： A""" DOI：10.12305/j.issn.1001-506X.2024.07.03

Simulation of geomagnetic anomaly of underwater vehicle based on

finite element method

ZHAO Gaoyang1， LIU Yong2， ZHU Pingjie2， XIANG Bing2， ZHOU Hongjuan1，*

（1. School of Information Science and Engineering， Harbin Institute of Technology （Weihai）， Weihai 264209， China;

2. 722 Research Institute， China State Shipbuilding Corporation， Wuhan 430205， China）

Abstract：

To overcome the difficulties of obtaining the magnetic anomaly data of underwater vehicles and the low calculation precision of existing magnet simulation models， a high-precision hybrid analytical model suitable in both near-field and far-field based on finite element method is introduced. The spatial signature of magnetic anomaly of an underwater vehicle with complex structure is analyzed by the finite element method. Then， an analytical model combined with uniformly-magnetized spheroid and a dipole array is proposed with the numerical results as the convergence goal. This analytical model consists of several parameters， including the number of array elements， the positions of each element， and magnetic moments， which are solved by nonlinear least-square algorithm. The simulation results show that the magnetic anomaly calculated from this model has good consistencies with the numerical solution in full space， giving an average error of 3% at the measurement plane. This model could be potentially extended to applications like magnetic field continuation and high-precision modeling.

Keywords：

airborne magnetic anomaly detection; finite element method （FEM）; magnet simulation; numerical fitting

＼=

0" 引" 言

水下航行器在日益重要的海洋開發(fā)中承擔(dān)著水下目標(biāo)、資源探測、氣象、海洋數(shù)據(jù)搜集等任務(wù)［1］，這使得對水下航行器的探測具有深刻的民用和軍事意義。對水下大型航行器的探測來說，隨著非磁性或弱磁性材料的應(yīng)用與降噪技術(shù)的發(fā)展，傳統(tǒng)聲學(xué)探測受到前所未有的挑戰(zhàn)，特別是在對聲信號具有巨大影響的淺海區(qū)域［23］。航空磁異常探測技術(shù)利用機(jī)載磁探儀對鐵磁性目標(biāo)進(jìn)行探測與定位，具有全天時和全地域的特點(diǎn)。作為一種非接觸式的被動探測方法［45］，航空磁異常探測具有識別能力強(qiáng)、定位精度高、隱蔽性良好等優(yōu)點(diǎn)［6］，在一定程度上能彌補(bǔ)水聲探測的應(yīng)用限制。

常用的水下航行器磁異常估算方法有數(shù)值法和磁體模擬解析法。數(shù)值法主要包括邊界元法、有限元法等［78］，可以較好地解決不規(guī)則物體的磁場模擬問題，缺點(diǎn)是計(jì)算量大且實(shí)時性差，無法適應(yīng)實(shí)際應(yīng)用。磁體模擬解析法通過特定磁矩磁源產(chǎn)生的磁場來模擬目標(biāo)磁場［910］。常用的模型有單磁偶極子、磁偶極子陣列、均勻磁化的旋轉(zhuǎn)橢球體等［11］。目前，普遍采用的模擬方法是：在靠近航行器一倍艇長的距離內(nèi)，航行器近似為單個旋轉(zhuǎn)橢球，其精度為75%～85%［12］;當(dāng)距離遠(yuǎn)至可忽略航行器物理尺寸時，航行器可簡化為單個磁偶極子或磁偶極子陣列模型，其中對磁偶極子陣列模型還涉及到多個磁偶極子的磁矩估計(jì)問題?？紤]近、遠(yuǎn)場模型的切換問題，出現(xiàn)了磁偶極子和旋轉(zhuǎn)橢球體的組合磁體模型，常用的組合磁體擬合方法需要求解多維超定方程，存在穩(wěn)定性與方程系數(shù)矩陣復(fù)雜的問題，同時建模過程中偶極子數(shù)量及位置改變也會影響方程求解。

在確定磁體模型磁矩值方面，目前多數(shù)研究基于實(shí)測數(shù)據(jù)或采用經(jīng)驗(yàn)值來確定，如一艘中型航行器（艇長為70～80 m）磁矩值為2×108 mA·m2［13］。采用單一經(jīng)驗(yàn)值的做法在理論研究上可以得到具有指導(dǎo)意義的定性結(jié)論，但在實(shí)際應(yīng)用中具有精度不高、靈活性不足的缺點(diǎn)。為獲得精確磁矩值，已有針對遺傳算法［14］、模擬退火算法［15］等智能算法的相關(guān)研究，但在討論中以固定的磁偶極子個數(shù)為前提并基于固定數(shù)目數(shù)據(jù)點(diǎn)討論建模誤差［16］，不能達(dá)到全面優(yōu)化磁矩值的目的，模擬方法仍然缺乏靈活性。

本文在已有的航行器磁異常磁體模擬的研究基礎(chǔ)上［17］，基于COMSOL Multiphysics多物理場軟件（簡稱為COMSOL軟件）仿真航行器目標(biāo)磁異常分布，通過有限元仿真方法獲得數(shù)值解以解決難以獲取實(shí)測數(shù)據(jù)的困難，進(jìn)一步建立參數(shù)包括可變陣元數(shù)目的偶極子陣列、磁偶極子位置以及磁矩的旋轉(zhuǎn)橢球體與偶極子陣列混合磁體模型，以仿真數(shù)值解為收斂目標(biāo)，實(shí)現(xiàn)具有良好穩(wěn)定性并且不區(qū)分近、遠(yuǎn)場探測的高精度磁體模擬方法的數(shù)值擬合研究。

1" 混合磁體模型

由磁疇磁化理論可知，組成航行器外殼體及內(nèi)部設(shè)備儀器的鐵磁性材料內(nèi)部自發(fā)磁化形成磁疇結(jié)構(gòu)［1819］，在地磁場的作用下各磁疇的排列方向近似一致［20］，形成感應(yīng)磁場，該感應(yīng)場由內(nèi)部磁矩方向相同的磁疇結(jié)構(gòu)的感應(yīng)場疊加而成［21］，因此在磁體模擬法等效模型中，對各個磁體磁矩可選取相同的磁化方向。

為提高航行器磁異常的估算精度，建立偶極子陣列、均勻磁化旋轉(zhuǎn)橢球體的混合磁體模型（以下簡稱混合模型），如圖1所示。

在圖1中，以旋轉(zhuǎn)橢球體的中心為原點(diǎn)，n元偶極子陣列等間距排布在旋轉(zhuǎn)橢球體長軸，橢球體長半軸為a，短半軸為b，與目標(biāo)的幾何尺寸相當(dāng)。

混合磁體在探測點(diǎn)產(chǎn)生的磁感應(yīng)強(qiáng)度B表示為以下矩陣形式：

B（r，n，M）=∑nj=0Bj=∑nj=0ajxajyajz

bjxbjybjz

cjxcjycjzMj（1）

式中：r為空間探測點(diǎn)相對于原點(diǎn)或偶極子中心的位置矢量;n為偶極子陣列陣元個數(shù);Mj表示磁矩，其中j=0代表旋轉(zhuǎn)橢球體的磁矩矢量，1≤j≤n時表示第j個磁偶極子的磁矩矢量;Mjx、Mjy、Mjz分別表示沿x、y、z方向的磁矩分量。由旋轉(zhuǎn)橢球體和磁偶極子磁場的推導(dǎo)，可得式（1）中相關(guān)系數(shù)［22］表示如下。

（1） j=0

ajx=-3μ04π12g2lnan+gan-g-ang2t，ajy=3μ0xy4πanb2nt，ajz=3μ0xz4πanb2nt

bjx=ajy，bjy=-3μ08πanb2ng2-12g2lnan+gan-g-2any2b4nt，bjz=3μ0anxz4πanb2nt

cjx=ajz，cjy=bjz，cjz=-3μ08πanb2ng2-12g2lnan+gan-g-2anz2b4nt

式中：an=12（x2+y2+z2+g2+t）;bn=an－g2t=（x2+y2+z2+g2）2－4g2x2;g=a2－b2。

（2） 1≤j≤n

ajx=14π3（x-xj）2r5j-1r3j，ajy=3（x-xj）y4πr5j，ajz=3（x-xj）z4πr5j

bjx=ajy，bjy=14π3y2r5j-1r3j，bjz=3yz4πr5j

cjx=ajz，cjy=bjz，cjz=14π3z2r5j-1r3j

式中：rj=（x－xj）2+（y－yj）2+（z－zj）2;（xj，yj，zj）為第j個偶極子的位置坐標(biāo)。

為提高擬合精度，由混合模型的磁感應(yīng)場分量Bx、By、Bz組成含待求系數(shù)的方程組：

Bx=∑nj=0（ajxMjx+ajyMjy+ajzMjz）（2）

By=∑nj=0（bjxMjx+bjyMjy+bjzMjz）（3）

Bz=∑nj＝0（cjxMjx+cjyMjy+cjzMjz）（4）

由電磁場的唯一性定理可知，求解電磁場的定解問題時，只要滿足給定的控制方程和邊界條件，其解必定唯一。因此，在求解中將探測點(diǎn)（x，y，z）坐標(biāo)值及對應(yīng)磁場值代入式（2）～式（4），即可建立含偶極子陣元數(shù)目、對應(yīng)磁矩參數(shù)、陣元位置的非線性方程組。

2" 航行器磁異常的有限元仿真

磁體模擬法需要部分磁場數(shù)據(jù)，目前難以獲得實(shí)際海域中目標(biāo)磁異常實(shí)測數(shù)據(jù)，且實(shí)測數(shù)據(jù)往往存在外界干擾、儀器噪聲甚至較多無效數(shù)據(jù)，因此利用數(shù)值計(jì)算法獲得磁異常的仿真應(yīng)用具有較強(qiáng)的理論研究價值。文獻(xiàn)［23］中討論了以球體驗(yàn)證有限元軟件COMSOL仿真結(jié)果與理論值的一致性，為簡化運(yùn)算，將通過該軟件對目標(biāo)磁特性進(jìn)行仿真。有限元計(jì)算原理是求解如下偏微分方程。

在無電流區(qū)域內(nèi)，靜磁場的基本方程如下：

Δ×H=0（5）

定義標(biāo)量磁勢Vm的表達(dá)式為

H=－ΔVm（6）

設(shè)目標(biāo)殼體材料的相對磁導(dǎo)率為μr，利用磁通密度與磁場的本構(gòu)關(guān)系，可得到：

Δ·B=0（7）

B=μ0μrH（8）

由此，得到標(biāo)量磁勢方程：

－Δ·（μ0μrΔVm）=0（9）

在本研究背景下，水下航行器外加磁場為地磁場Bt，待求解的微分方程為

－Δ·（μ0μrΔVm+Bt）=0（10）

按照艦長為100 m、艇直徑為10 m的中型航行器建立幾何模型，對殼體參數(shù)、材料相對磁導(dǎo)率、背景磁場等參數(shù)盡量按照真實(shí)參數(shù)設(shè)置。如圖2所示，構(gòu)建水下航行器模型并進(jìn)行網(wǎng)格劃分，完成計(jì)算。

圖3以探測高度為500 m的平面目標(biāo)磁感應(yīng)場分布為例。受到目標(biāo)物理形狀的影響，無論沿橫向還是縱向，其感應(yīng)磁場的分布并不滿足嚴(yán)格對稱分布。在傳統(tǒng)磁體模擬方法中，通常采用單個均勻磁化橢球體或偶極子模型來進(jìn)行近似模擬，模擬的磁感應(yīng)場呈對稱分布，與實(shí)際磁異常分布不符，勢必導(dǎo)致在某些探測位置上磁場強(qiáng)度失真較大。傳統(tǒng)的單磁體模擬法對理論研究具有重要的意義，但不適用于對磁異常測量精度要求較高的實(shí)際應(yīng)用場合，因此構(gòu)建合適的混合磁體模型實(shí)現(xiàn)對目標(biāo)感應(yīng)磁場的高精度擬合非常必要。

3" 混合磁體模型系數(shù)求解及驗(yàn)證

3.1" 系數(shù)求解原理及過程

本文采用非線性最小二乘算法列文伯格馬夸爾特（Levenberg-Marquarat， L-M）算法對待定系數(shù)進(jìn)行求解［24］。非線性最小二乘算法的原理是從某個初值開始按照搜索方向持續(xù)迭代計(jì)算，再通過新方向的調(diào)整，直至目標(biāo)函數(shù)收斂到極小值點(diǎn)，較常用的算法有最速下降法［25］、牛頓法、高斯牛頓法等。上述算法在本文中均存在一定的使用缺陷：最速下降法在接近于最優(yōu)解時收斂速度較慢，因此需要較多迭代次數(shù);牛頓法需計(jì)算目標(biāo)函數(shù)的二階導(dǎo)矩陣，在數(shù)據(jù)規(guī)模較大時計(jì)算困難［26］;高斯牛頓法因近似二階泰勒展開，只在展開點(diǎn)附近有較好結(jié)果，在迭代中可能出現(xiàn)奇異或病態(tài)矩陣，或者因局部近似無法實(shí)現(xiàn)收斂［2728］。作為改進(jìn)，L-M算法為增量添加信任域，該區(qū)域范圍依近似模型與實(shí)際函數(shù)的差異自適應(yīng)調(diào)整［29］。

目標(biāo)函數(shù)定義為磁體模型與實(shí)際函數(shù)間誤差平方和R：

R=∑Ni=1［f2（ri）－f1（ri，n，M）］2（11）

式中：i=0，1，2，…，N表示用以求解的探測點(diǎn);ri為探測距離;n為偶極子陣元數(shù)目;M為對應(yīng)磁矩系數(shù);f1表示含待求參數(shù)的磁體模型;f2表示實(shí)際值函數(shù)，在本文中為COMSOL軟件磁場仿真值。該算法通過多次迭代計(jì)算，控制目標(biāo)函數(shù)最小以獲得待求參數(shù)，主要包括以下步驟。

步驟 1" f1待求參數(shù)p，初始點(diǎn)p0，終止控制常數(shù)ε0=f2－f1（pk），參數(shù)k=0，λ0=10-3，v可以是任何大于1的數(shù)。

步驟 2" Jacobi矩陣為Jk，構(gòu)造增量正規(guī)方程N(yùn)-kδk=JTkεk，其中N-k=JTkJk+λkI。

步驟 3" 求解增量正規(guī)方程得到δk。

步驟 3.1" 若f2－f1（pk+δk）lt;εk，令pk+1=pk+εk，當(dāng)δklt;ε時輸出結(jié)果;否則迭代繼續(xù)，并令λk+1=λk/v。

步驟 3.2" 如果f2－f1（pk+δk）≥εk，令λk+1=λk·v，重新解正規(guī)方程得到δk，返回步驟1。

如圖4所示，水下航行器外殼形狀可能不完全規(guī)則，參照其形狀特征建立適當(dāng)尺寸橢球體并取其長軸方向上的中心點(diǎn)為原點(diǎn)?？臻g探測點(diǎn)Qi對應(yīng)磁感應(yīng)場Bi=（Bix， Biy， Biz）。

圖4中，rij=（rijx， rijy， rijz，）為探測點(diǎn)Qi與磁偶極子陣列第j個陣元相對距離。

建立混合模型含參方程組如下，利用最小二乘算法求解待求磁矩參數(shù)，并建立橢球體與磁偶極子陣列的混合磁體模型。

Bix=∑nj＝0（ajxMjx+ajyMjy+ajzMjz）

Biy=∑nj＝0（bjxMjx+bjyMjy+bjzMjz）

Biz=∑nj＝0（cjxMjx+cjyMjy+cjzMjz）（12）

3.2" 擬合實(shí)驗(yàn)結(jié)果與分析

下面將通過實(shí)驗(yàn)驗(yàn)證擬合模型。通過L-M算法求解混合模型中的待定參數(shù)，在擬合過程中可根據(jù)應(yīng)用條件需求

設(shè)定陣元數(shù)目范圍、對最佳陣元數(shù)進(jìn)行取整求解，選擇與探測點(diǎn)磁場值誤差最小的L-M算法結(jié)果，算法初值參考經(jīng)驗(yàn)值。后文進(jìn)行了多組實(shí)驗(yàn)討論取點(diǎn)方式、陣元數(shù)目等參數(shù)影響，增大模型驗(yàn)證范圍。假定水下航行器艇長為100 m，以建立橢球體與五元偶極子陣列混合模型為例，基于數(shù)值計(jì)算結(jié)果隨機(jī)選取探測點(diǎn)，分別以完全規(guī)律性、一定規(guī)律性、完全隨機(jī)性3種方式選取探測點(diǎn)，討論選取探測點(diǎn)的規(guī)律性對最終模型擬合精度的影響。

（1） 擬合使用符合完全規(guī)律性探測點(diǎn)

在固定探測高度平面上選取規(guī)律分布探測點(diǎn)如圖5所示，在500 m探測高度平面上，以均勻間隔在與航行器運(yùn)行方向夾角為45°的航線上選取探測點(diǎn)Qi=（x，y，500），x，y∈（-500，500），i=0，1，2，…，100。將對應(yīng)位置坐標(biāo)與磁場值代入式（12），并利用L-M算法對磁矩參數(shù)進(jìn)行擬合。

在500 m探測高度平面上，對探測平面磁場分布的數(shù)值仿真值與混合模型計(jì)算值等效效果進(jìn)行評估，分別繪制磁感應(yīng)場與3個分量，如圖6所示。計(jì)算磁場計(jì)算值與仿真值間峰值差：x分量差值為0.050 nT，y分量差值為0.039 nT，z分量差值為0.066 nT，磁場總量差值為0.046 nT，通過繪制場強(qiáng)等值線分布圖可知磁場分布擬合結(jié)果較接近。

（2） 擬合使用符合一定規(guī)律性探測點(diǎn)

在固定探測高度平面上隨機(jī)選取相同數(shù)目探測點(diǎn)如圖7所示，探測點(diǎn)Qi=（x， y，500），x，y∈（-500，500），i=0，1，2，…，100;將對應(yīng)位置坐標(biāo)與磁場值代入式（11）進(jìn)行擬合。

進(jìn)行后續(xù)的磁場分布計(jì)算，同樣選取500 m探測高度測平面，對該取點(diǎn)方式計(jì)算得到混合模型磁場分布計(jì)算值，與數(shù)值仿真值等效效果進(jìn)行評估，分別繪制磁感應(yīng)場與3個分量，如圖8所示。計(jì)算磁場計(jì)算值與仿真值間峰值差：x分量差值為0.089 nT，y分量差值為0.029 nT，z分量差值為0.052 nT，磁場總量差值為0.004 nT，雖然場值分布有差異但差值較小。

（3） 擬合使用符合完全隨機(jī)性探測點(diǎn)

更一般地，在選取探測點(diǎn)時不限制探測高度，在探測空間中完全隨機(jī)選擇相同數(shù)目探測點(diǎn)，Qi=（x，y，z），x，y∈（-500，500），z∈（400，1 200），i=0，1，2，…，100。在本文的討論中，探測空間由磁探儀作用距離進(jìn)行定義：國外領(lǐng)先技術(shù)將有效作用范圍提升到1 200 m左右;有磁探儀據(jù)稱有效作用距離大約為400～600 m，因此將有效探測距離設(shè)定為航行器上方400～1 200 m距離。在該空間內(nèi)隨機(jī)選取探測點(diǎn)，如圖9所示，利用探測點(diǎn)坐標(biāo)信息與磁場數(shù)值進(jìn)行擬合。

在模型計(jì)算值與數(shù)值仿真值的對比中，同樣選取500 m探測高度測平面，分別繪制磁感應(yīng)場與3個分量等效效果，如圖10所示。計(jì)算磁場計(jì)算值與仿真值間峰值差：x分量差值為0.038 nT，y分量差值為0.033 nT，z分量差值為0.057 nT，磁場總量差值為0.016 nT。

在混合模型建模時分別用3種選取探測點(diǎn)方法計(jì)算待求參數(shù)，3種不同模型對目標(biāo)空間磁場分布計(jì)算值均較接近仿真值。圖11以平均絕對百分比誤差（mean absolute percentage error， MAPE）量化目標(biāo)磁場分布計(jì)算值與仿真值間誤差，在計(jì)算誤差的過程中擴(kuò)大探測平面數(shù)據(jù)密度，驗(yàn)證點(diǎn)坐標(biāo)范圍為x，y∈（-1 000， 1 000），z∈（400，1 200）。

分析不同取點(diǎn)方式，建立模型計(jì)算誤差，3種模型計(jì)算誤差均較穩(wěn)定，應(yīng)用在不同模型使用條件下可適當(dāng)調(diào)整擬合取點(diǎn)方式：采用完全規(guī)律與一定規(guī)律取點(diǎn)方法，擬合所用探測點(diǎn)來自同一高度平面，在該平面及附近高度平面內(nèi)誤差較小。采用完全隨機(jī)取點(diǎn)方法擬合所用探測點(diǎn)數(shù)據(jù)范圍更大，增強(qiáng)了模型的魯棒性，計(jì)算誤差比較穩(wěn)定，特別是應(yīng)用在較遠(yuǎn)探測距離時使模型誤差降低，對后續(xù)數(shù)值模型的應(yīng)用與研究有一定的啟發(fā)作用。

3.3" 擬合結(jié)果應(yīng)用

聯(lián)合有限元仿真與磁體模擬方法獲得高精度的數(shù)值擬合模型，理論上可以通過航行器目標(biāo)引起地磁異常分布，判斷大致位置。

以圖6為例進(jìn)行分析，航行器形狀因?yàn)榫哂袑ΨQ、規(guī)則結(jié)構(gòu)，磁感應(yīng)場及3個分量具有明顯對稱性;x分量具有很強(qiáng)的區(qū)域特征，在航行器位置處出現(xiàn)1個負(fù)峰值;y分量有4個較突出磁場峰值，其中2個正峰值分別出現(xiàn)在目標(biāo)的左前側(cè)和右后側(cè)、2個負(fù)峰值出現(xiàn)在左后側(cè)和右前側(cè)，且右側(cè)峰值較左邊更大，右側(cè)峰值位置以航行器位置為中心沿y軸對稱;磁異常z方向分量峰值最大，對目標(biāo)地磁異?？倛鲇绊懽畲?，在探測平面上各有1個磁場正峰值與負(fù)峰值，正峰值位置出現(xiàn)在航行器運(yùn)行方向的前側(cè)，與磁異?？偭康?/p>

峰值位置大致相同，負(fù)峰出現(xiàn)在后側(cè)，目標(biāo)大概位置在磁場負(fù)峰值中心與正峰值中心中間距離處。按照以上的分布特征，由磁場分布及磁場分量特征可確定航行器位于x∈（-80， 50）、y=0，與實(shí)際仿真條件大致相同，驗(yàn)證了該方法在確定目標(biāo)位置上的可行性。

4" 混合模型建模精度分析

4.1" 驗(yàn)證模型精度改善

本文在建模方法上對單個磁體模擬方法進(jìn)行改進(jìn)，經(jīng)前文對比與討論，實(shí)驗(yàn)中建立的五元偶極子陣列混合模型實(shí)現(xiàn)較好擬合。為進(jìn)一步對建模精度進(jìn)行分析，在相同的探測條件下對以下磁體模型等效性進(jìn)行評價。

Ⅰ模型：均勻磁化的單個旋轉(zhuǎn)橢球體模型。

Ⅱ模型：橢球體與五元偶極子陣列混合模型。

Ⅲ模型：橢球體與九元偶極子陣列混合模型。

求解磁體模型中的待求系數(shù)。Ⅰ模型中未知磁矩值參考常用經(jīng)驗(yàn)公式，Ⅱ、Ⅲ模型通過L-M算法求解參數(shù)，在建模時保證模型結(jié)構(gòu)尺寸等相同。建模完成后，在500 m探測高度測平面上分別計(jì)算3個模型磁感應(yīng)場及分量，比較與數(shù)值仿真值間誤差。通過決定系數(shù)R2、均方誤差（root mean square error， RMSE）、殘差平方和（residual sum of squares， RSS）量化誤差，擬合模型最理想誤差指標(biāo)應(yīng)為R2接近1、RMSE和RSS接近0，說明該磁體模型精度越高，擬合模型空間實(shí)際磁場分布特性越吻合。

表1中，B是數(shù)值仿真結(jié)果，平均值為B-，不同模型計(jì)算值為B′。3個模型磁感應(yīng)場在探測平面上的MAPE值分別為21.55%、3.59%、3.94%。本文提出的混合模型是基于單個磁體模擬方法的改進(jìn)，Ⅱ、Ⅲ模型的磁感應(yīng)場及分量的R2更接近1，誤差指標(biāo)RMSE、RSS較Ⅰ模型有較明顯優(yōu)化，完全滿足工程的高精度需要。作為對比，組成Ⅱ、Ⅲ模型的偶極子陣列陣元數(shù)目不同，Ⅱ模型在磁感應(yīng)場及y、z兩個方向分量上較Ⅲ模型誤差更小，因此并非一味靠增加磁偶極子數(shù)目來提高擬合精度。綜上，由橢球體與偶極子陣列組成的混合模型從精度與模型靈活度兩個方面實(shí)現(xiàn)了模型優(yōu)化：通過L-M算法確定磁矩參數(shù)，較磁矩經(jīng)驗(yàn)值方法進(jìn)一步提高擬合精度，通過調(diào)整偶極子陣列陣元個數(shù)，在滿足精度需要的前提下提高模型的靈活度。

4.2" 驗(yàn)證相關(guān)參數(shù)影響

在驗(yàn)證本文提出的混合模型提高磁場分布計(jì)算精度的基礎(chǔ)上，為避免單一尺寸模型結(jié)果的偶然性，增加不同尺寸航行器建模實(shí)驗(yàn)、不同偶極子陣元數(shù)目建模實(shí)驗(yàn)，討論磁體模型中偶極子陣列陣元個數(shù)與探測距離等相關(guān)參數(shù)與建模精度的相關(guān)關(guān)系。以探測平面上磁場分布計(jì)算值與仿真值間的平均絕對百分比誤差（后文簡稱平均誤差）作為主要評價指標(biāo)，分別進(jìn)行以下兩個實(shí)驗(yàn)。

實(shí)驗(yàn) 1" 研究不同尺寸水下航行器建模影響。參考文獻(xiàn)［30］中的實(shí)驗(yàn)選取3種尺寸的航行器，按照長寬尺寸分為小、中、大3種艇型，小型艇模型參數(shù)：艇長57.2 m，艇直徑8.6 m;中型艇模型參數(shù)：艇長100 m，艇直徑10 m;大型艇模型參數(shù)：艇長153 m，艇直徑17.3 m。為使結(jié)果具有可對比性，雖然目標(biāo)航行器尺寸改變，但保持探測高度與艇長的相對探測尺度不變，分析3種艇型磁場計(jì)算結(jié)果的探測平面高度分別為286 m、500 m、765 m。

在對不同艇型航行器目標(biāo)的混合模型建模中，分別調(diào)整混合模型中磁偶極子陣列的陣元個數(shù)n=2，3，…，12，在對應(yīng)探測平面上由模型計(jì)算得到的磁場值與有限元仿真值間平均誤差如圖12所示。通過對不同陣元個數(shù)模型誤差進(jìn)行分析可知，混合模型的磁體模擬方法受偶極子陣列陣元數(shù)目影響，3種艇型航行器混合模型對磁場分布的計(jì)算誤差均可控制在11%以內(nèi);當(dāng)偶極子陣列陣元數(shù)目達(dá)到一定程度時，建模誤差趨于穩(wěn)定，理論上可以通過調(diào)整陣元數(shù)目獲得最佳擬合模型。

實(shí)驗(yàn) 2" 研究不同偶極子陣元數(shù)目建模影響。建立某中型艇（艇長100 m，艇直徑10 m）的磁場分布混合模型，偶極子陣列陣元數(shù)目為n=2，3，…，12，在500 m探測高度探測平面上選取相同數(shù)目探測點(diǎn)。計(jì)算不同陣元數(shù)目的混合模型磁場值，與數(shù)值仿真值間誤差與探測高度間關(guān)系如圖13所示。雖然混合模型對磁場分布的計(jì)算受探測高度，即探測距離影響較大，但在400～700 m探測范圍內(nèi)磁場分布計(jì)算值平均誤差控制在10%內(nèi)，模型具有一定的精度穩(wěn)定性。保持相同探測高度，混合模型中偶極子陣列陣元數(shù)目對模型精度影響不明顯，但較多陣元數(shù)目的混合模型在探測距離較遠(yuǎn)表現(xiàn)更加穩(wěn)定。

混合模型磁體模擬法有較高的建模實(shí)用價值。在本探測條件下建立的四元偶極子陣列混合模型，在探測點(diǎn)所在探測平面上磁感應(yīng)場計(jì)算平均誤差為3%、在探測空間最遠(yuǎn)距離處為32%。作為對比，均勻磁化單個橢球體模型在相同探測距離時磁場計(jì)算誤差分別為22%、47%，在本實(shí)驗(yàn)條件下可認(rèn)為四元偶極子陣列混合模型為最佳模型。通過多組實(shí)驗(yàn)得到以下分析結(jié)論：

（1） 磁偶極子陣列陣元數(shù)目是影響混合模型計(jì)算精度的可變量，相關(guān)待求參數(shù)如磁矩值、陣元分布位置，同時受到陣元數(shù)目的影響。

（2） 探測距離對混合模型空間磁場分布影響較大，分布隨著探測距離增大，空間磁場對稱特性出現(xiàn)畸變或峰值轉(zhuǎn)移。

（3） 采用不同磁體模擬方法對磁場計(jì)算結(jié)果有較大影響。本文提出的通過L-M算法調(diào)整混合模型參數(shù)方法，比應(yīng)用經(jīng)驗(yàn)公式的單個磁體模擬方法計(jì)算更簡單、在探測距離增加時結(jié)果更有效。

通過以上在不同尺寸航行器、不同偶極子陣列陣元數(shù)目、不同探測距離等相關(guān)參數(shù)下的實(shí)驗(yàn)理論分析，證明在實(shí)際磁異常探測技術(shù)的條件下如獲得任意探測點(diǎn)的位置及磁異常探測數(shù)據(jù)，不僅可實(shí)現(xiàn)與水下航行器相對應(yīng)尺寸的混合模型空間磁場分布計(jì)算，也可通過對相關(guān)參數(shù)的調(diào)整獲得最佳擬合模型，滿足高精度應(yīng)用要求。

5" 結(jié)" 論

伴隨水下航行器目標(biāo)環(huán)境日益復(fù)雜，針對復(fù)雜磁源目標(biāo)的定位識別需建立更高精度磁場擬合模型。本文綜合有限元法與磁體模擬法的應(yīng)用研究，充分利用有限元法獲得復(fù)雜結(jié)構(gòu)航行器目標(biāo)的空間磁感應(yīng)強(qiáng)度分布，建立旋轉(zhuǎn)橢球體與偶極子陣列的高精度混合磁體數(shù)值模型與待求模型系數(shù)方法，討論不同航行器目標(biāo)混合陣列建模影響因素與最佳擬合模型，驗(yàn)證了該模型較好的擬合優(yōu)度與探測空間泛化性。本文討論了磁體模擬方法混合模型在全探測空間建模、模型中精確磁矩值、提升計(jì)算精度等多個方面的改進(jìn)，拓展實(shí)現(xiàn)磁場數(shù)據(jù)延拓、高精度磁體模擬方法的理論應(yīng)用。

參考文獻(xiàn)

［1］ 宋保維， 潘光， 張立川， 等. 自主水下航行器發(fā)展趨勢及關(guān)鍵技術(shù)［J］. 中國艦船研究， 2022， 17（5）： 2744.

SONG B W， PAN G， ZHANG L C， et al. Development trend and key technologies of autonomous underwater vehicles［J］. Chinese Journal of Ship Research， 2022， 17（5）： 2744.

［2］ WELLER R A， LEVIN L A， GATES A R. Global observing needs in the deep ocean［J］. Frontiers in Marine Science， 2019， 6（6）： 167218016.

［3］ WANG C， CUI Y， SONG X， et al. A novel underwater target detection method based on low-frequency magnetic signal generated by portable transmitter［J］. IEEE Sensors Journal， 2023， 23（8）： 84598465.

［4］ SANCHEZ P， PAPAELIAS M， MARQUEZ F. Autonomous underwater vehicles： instrumentation and measurements［J］. IEEE Instrumentation and Measurement Magazine， 2020， 23（2）： 105114.

［5］ 鞠建波， 郁紅波， 范趙鵬， 等. 吊放聲吶擴(kuò)展螺旋陣搜潛效能評估改進(jìn)方法［J］. 水下無人系統(tǒng)學(xué)報， 2020， 28（4）： 434439.

JU J B， YU H B， FAN Z P， et al. An improved method for evaluation the submarine searching effectiveness with extension spiral array of dipping sonar［J］. Journal of Unmanned Undersea Systems， 2020， 28（4）： 434439.

［6］ 胡平， 顧雪峰， 劉凱， 等. 國外航空反潛裝備現(xiàn)狀及發(fā)展趨勢［J］. 艦船電子工程， 2022， 42（9）： 1012， 86.

HU P， GU X F， LIU K， et al. Current situation and development trend of foreign aviation anti-submarine equipment［J］. Ship Electronic Engineering， 2022， 42（9）： 1012， 86.

［7］ BRUNOTTE X， MEUNIER G， BONGIRAUD J. Ship magnetizations modelling by finite element method［C］∥Proc.of the IEEE Conference on Electromagnetic Field Computation， 1992.

［8］ WANG J X， REN Q. A 3-D hybrid maxwell’s equations finite-difference time-domain （ME-FDTD）/wave equation finite-element time-domain （WE-FETD） method［J］. IEEE Trans.on Antennas and Propagation， 2023， 71（6）： 52125220.

［9］ CHADEBEC O， COULOMB J L， BONGIRAUD J P， et al. Recent improvements for solving inverse magnetostatic problem applied to thin shells［J］. IEEE Trans.on Magnetics， 2002， 38（2）： 10051008.

［10］ IOANNIDIS G. Identification of a ship or submarine from its magnetic signature［J］. IEEE Trans.on Aerospaceamp;Electronic Systems， 1977， 13： 327329.

［11］ 曲曉慧， 楊日杰， 單志超. 潛艇磁場建模方法的分析與比較［J］. 艦船科學(xué)技術(shù)， 2011， 33（3）： 711.

QU X H， YANG R J， SHAN Z C. Analysis and comparison on magnetic field modeling method of submarine［J］. Ship Science and Technology， 2011， 33（3）： 711.

［12］ 林春生. 艦船磁場深度換算中的單列橢球體和磁偶極子混合模型［J］. 水雷戰(zhàn)與艦船防護(hù)， 1996， 1996（3）： 5458.

LIN C S. Hybrid model of single-row ellipsoid and magnetic dipole in ship magnetic field depth conversion［J］. Mine Warfare ＆ Ship Self-defence， 1996， 1996（3）： 5458.

［13］ HOFFMANN A. Spin hall effects in metals［J］. IEEE Trans.on Magnetics， 2013， 49（10）： 51725193.

［14］ 徐杰. 基于遺傳算法的潛艇三維模型及其磁場延拓方法［J］. 艦船電子工程， 2009， 29（7）： 188191.

XU J. Application of genetic algorithm in continuation of magnetic field of submarines in high areas［J］. Ship Electronic Engineering， 2009， 29（7）： 188191.

［15］ 招妙妍. 基于模擬退火算法的航向?qū)崪y磁場數(shù)據(jù)建模分析［J］. 艦船電子工程， 2019， 39（9）： 188190.

ZHAO M Y. Modeling and analysis of directional field data based on simulated annealing algorithm［J］. Ship Electronic Engineering， 2019， 39（9）： 188190.

［16］ NARA T， SUZUKI S， ANDO S. A closed form formula for magnetic dipole localizati-on by measurement of its magnetic field and spatial gradients［J］. IEEE Trans.on Magnetics， 2006， 42（10）： 32913293.

［17］ QIN Y J， LI K Y， YAO C， et al. Magnetic anomaly detection using full magnetic gradient orthonormal basis function［J］. IEEE Sensors Journal， 2020， 20（21）： 1292812940.

［18］ 張志東. 磁性材料的磁結(jié)構(gòu)、磁疇結(jié)構(gòu)和拓?fù)浯沤Y(jié)構(gòu)［J］. 物理學(xué)報， 2015， 64（6）： 521.

ZHANG Z D. Magnetic structures， magnetic domains and topological magnetic textures of magnetic materials［J］. Acta Physica Sinica， 2015， 64（6）： 521.

［19］ JILES D， KIARIE W. An integrated model of magnetic hysteresis， the magnetomechanical effect， and the barkhausen effect［J］. IEEE Trans.on Magnetics， 2023， 57（2）： 111.

［20］ CHUNG H J， YANG C S， JUNG W J. A magnetic field separation technique for a scaled model ship through an earth’s magnetic field simulator［J］. Journal of Magnetics， 2015， 20（1）： 6268.

［21］ HALGEDAH F. The dependence of magnetic domain structure upon magnetization state with emphasis upon nucleation as a mechanism for pseudo-single-domain behavior［J］. Journal of Geophysical Research， 1983： 65056522.

［22］ XIN Z， XU Q， LI Q， et al. An orthonormalized basis function based narrowband filtering algorithm for magnetic anomaly detection［C］∥Proc.of the International Congress on Image amp; Signal Processing， 2016.

［23］ 孫曉永， 張琦， 潘孟春， 等. 基于COMSOL Multiphysics的目標(biāo)磁特性仿真分析［J］. 中國測試， 2017， 43（1）： 122126.

SUN X Y， ZHANG Q， PAN M C， et al. Simulation analysis of the magnetic properties based on COMSOL multiphysics［J］. China Measurementamp;Test， 2017， 43（1）： 122126.

［24］ RUBIO J D J. Stability analysis of the modified levenberg-marquardt algorithm for the artificial neural network training［J］. IEEE Trans.on Neural Networks and Learning Systems， 2020， 32（8）： 35103524.

［25］ ANDREI N. Steepest descent methods［J］. Springer Optimization and its Applications， 2022， 195（1）： 81107.

［26］ HANSEN E R， GREENBERG R I. An interval Newton method［J］. Applied Mathematics amp; Computation， 1983， 12（2/3）： 8998.

［27］ VASIN V V. Solving nonlinear inverse problems based on the regularized modified Gauss-Newton method［J］. Doklady Mathe-matics， 2022， 105（3）： 175177.

［28］ GE J， WANG S Q， DONG H B， et al. Real-time detection of moving magnetic target using distributed scalar sensor based on hybrid algorithm of particle swarm optimization and Gauss-Newton method［J］. IEEE Sensors Journal， 2020， 20（18）： 1071710723.

［29］ YANG J， ZOU J. Parameter estimation of a horizontally multilayered soil with a fast evaluation of the apparent resistivity and its derivatives［J］. IEEE Access， 2020， 8： 5265252662.

［30］ 戴忠華， 張曉兵， 周穗華. 艦船磁場磁偶極子陣列模型的適用性研究［J］. 武漢理工大學(xué)學(xué)報， 2020， 42（4）： 99104.

DAI Z H， ZHANG X B， ZHOU S H. Applicability of magnetic dipole array model for ship magnetic field［J］. Journal of Wuhan University of Technology， 2020， 42（4）： 99104.

作者簡介

趙高陽（1997—），女，碩士研究生，主要研究方向?yàn)榇女惓Ｌ綔y。

劉" 勇（1968—），男，研究員，碩士，主要研究方向?yàn)榈皖l電磁通信。

朱平杰（1978—），男，高級工程師，主要研究方向?yàn)榈皖l電磁天線。

向" 冰（1991—），女，高級工程師，碩士，主要研究方向?yàn)榈皖l電磁信號處理。

周洪娟（1980—），女，副教授，博士，主要研究方向?yàn)楹Ｑ蟠盘綔y、電波傳播。