艦載反魚雷魚雷對抗行為分析與建模

陳衛偉，閔紹榮，李明輝，謝紅勝

中國艦船研究設計中心，湖北武漢430064

0 引 言

魚雷技術的發展使得魚雷的智能化程度不斷提高，使其具有很好的識別和跟蹤目標的能力。在這種情況下，魚雷對抗聲誘導、噪聲干擾等軟殺傷的能力大大加強。采用以反魚雷魚雷（Anti-Torpedo Torpedo，ATT）為主的硬殺傷裝備對抗來襲魚雷日益成為艦艇反魚雷技術采用的重要手段之一［1］。

現有的關于艦載ATT 對抗魚雷的研究通常不考慮對抗誤差，在二維平面上建立艦載ATT 對抗數學模型，對其對抗行為進行建模與研究［2-4］。這些研究取得了一些成果，但也存在著不足：首先，艦艇反魚雷作戰的戰場環境涉及水下三維立體空間，而在二維平面上建立的對抗數學模型卻忽略了魚雷、反魚雷魚雷等各作戰實體的深度變化帶來的影響，因而難以準確地對其作戰行為進行描述；其次，忽略了對抗誤差，可能會帶來仿真結果的失真。

在上述研究的基礎上，本文擬基于蒙特卡洛法建立三維態勢下的艦載ATT 對抗數學模型，在考慮魚雷報警距離誤差、報警舷角誤差、速度誤差和航向誤差等誤差的基礎上，對ATT 直接攔截與迎面逆彈道攔截兩種攔截方案進行分析與建模，從而為進一步深入研究反魚雷對抗方案的效能評估問題提供參考。

1 ATT 攔截方案描述

由于來襲魚雷與水面艦艇之間對抗態勢的多樣性，ATT 的使用存在多種方案。根據來襲魚雷的攻擊彈道，以及魚雷和水面艦艇的運動參數，ATT 的攔截彈道主要包括以下2 種方案［5］：

1）直接攔截方案。

設艦艇發射ATT 時，來襲魚雷已進入自導導引狀態，并以固定提前角彈道攻擊艦艇。如圖1所示，以水面艦艇發射ATT 時的位置為坐標原點，艦艇航向Cw為X0軸建立右手直角坐標系。

圖1 ATT 直接攔截方案示意圖Fig.1 Schematic plan of ATT direct interception

艦艇以βatt的發射角發射ATT，沿直線搜索目標，在A 點附近區域與來襲魚雷相遇。

2）迎面逆彈道攔截方案。

按照上述方法建立相同的坐標系，如圖2 所示。ATT 發射后，以半徑Ratt和航速Vatt做旋回運動，轉至來襲魚雷的逆彈道，在A 點附近區域與來襲魚雷相遇。

圖2 ATT 迎面攔截方案示意圖Fig.2 Schematic plan of ATT head-on interception

2 對抗數學模型

2.1 ATT 運動學模型

ATT 運動模型分為3 段：入水尋深段、程序搜索段和自導追蹤段。分段依次對其運動過程進行建模，具體模型如下。

1）ATT 入水尋深段彈道模型。

式中：Vatr為ATT 入水速度；Az為入水后的阻尼系數；Catt為ATT 發射航向。

2）ATT 程序搜索段彈道模型。

根據不同的攔截方案，ATT 程序搜索段采用不同的搜索模式。若ATT 采用直接攔截方案，則其程序搜索彈道為直航搜索，如下式所示：

式中，Vatt為ATT 航速。在直航搜索段，ATT 的航向保持不變。

若ATT 采用迎面逆彈道攔截方案，則其程序搜索彈道為旋回搜索。設ATT 在水平面內進行旋回，旋回角速度為ωatt，旋回半徑為ratt，旋回角度為αatt。以ATT 開始旋回的點為原點，以該時刻（ts）ATT 的水平面航向為y0軸建立直角坐標系，則其旋回搜索彈道模型為：

旋回搜索過程中，ATT 的航向為

3）ATT 自導追蹤彈道模型。

在近距離魚雷報警情況下，來襲魚雷的速度較快，與ATT 呈相向運動態勢，視角變化大，因此ATT 需在較短的時間和航程內攔截魚雷。常規的尾追彈道由于需要消耗大量的尾追時間和航程，因此不再適用于ATT 追蹤模型。本文選用固定提前角彈道作為ATT 的追蹤彈道。

根據魚雷與ATT 之間的相對運動關系，ATT在固定提前角式導引下航向Catt為：

2.2 來襲魚雷運動學模型

目前，采用定深航直航方式搜索目標，采用尾追式的追蹤彈道跟蹤目標直至命中目標是魚雷最常見的作戰方式，為簡化模型并不失一般性，本文將魚雷的彈道模型分為定向直航階段與尾追階段。

魚雷的直航搜索彈道模型與ATT 相同，見式（2）。

尾追法導引關系是指魚雷速度向量CT與目標視線（魚雷與目標的連線）相重合［6］。

設t 時刻魚雷捕獲到的目標坐標為(XM，YM，ZM)，魚雷坐標為(XT，YT，ZT)，則在t+1時刻，魚雷尾追跟蹤彈道的軌跡可表示為：

尾追跟蹤段的魚雷航向CT為：

式中：(XM，YM，ZM)為追蹤目標t時刻的空間坐標；DMT為兩者之間的距離。

2.3 ATT/魚雷自導模型

根據目前的自導信號處理水平［7］，無論是艦載ATT 還是來襲魚雷，主動自導與被動自導都是最主要的自導方式。

1）被動自導模型。

被動自導模型由被動聲吶方程描述：

式中：SLT為目標的輻射噪聲聲源級；TL(Rta)為魚雷至目標距離Rta的聲傳播損失；DI 為魚雷的指向性因子；NL 為魚雷的自噪聲級；DT 為魚雷的被動檢測閾。

2）主動自導模型。

魚雷聲自導主動檢測模型是以主動聲吶方程為基礎進行建模，本文的主動自導方程采用噪聲掩蔽下的主動聲吶方程：

式中：SL 為魚雷自導發射聲源級；TS 為目標反射強度。

2.4 對抗誤差模型

艦艇在水下對抗魚雷過程中的誤差模型主要包括魚雷報警距離誤差ΔDM、魚雷報警舷角誤差ΔQM、魚雷速度誤差ΔVMT和魚雷航向誤差ΔCMT。設ΔDM，ΔVMT和 ΔCMT分別服從的正態分布，則艦艇聲吶測量得到的魚雷報警距離、報警舷角、魚雷速度和魚雷航向分別為：

式中：DM，QM，VMT和CMT分別為真實情況下的相關參數；σi為標準正態分布隨機數。

3 仿真算例及結果分析

根據ATT 攔截策略的不同，本文分別以迎面攔截和直接攔截制定2 個推演算例。

上述模型中，假設水面艦艇的初始速度為18 kn，最高規避航速為32 kn，旋回半徑為300 m，旋回角速度為1.5（°）/s；艦艇拖曳陣列聲吶線纜長500 m，深40～60 m；來襲魚雷速度45 kn，自導扇面角90°，旋回角速度6（°）/s，直航搜索深度40～80 m，主動聲源級180 dB，自噪聲60 dB，主、被動檢測閾20 dB，魚雷自導頻率30 kHz，方向性指數30 dB，魚雷總航程20 000 m；ATT 速度50 kn，自導扇面角60°～150°，旋回角速度6（°）/s，自噪聲60 dB，主、被動檢測閾20 dB。

為了保證實驗數據的精度，本文取蒙特卡洛實驗次數為8 000 次，仿真步長0.5。

圖3～圖4 所示為第89 次仿真中，ATT 采用直接攔截方案成功攔截來襲魚雷的三維與二維態勢圖。圖5～圖6 所示為第173 次仿真中，ATT 采用直接攔截方案未成功規避來襲魚雷的三維與二維態勢圖。

圖3 ATT 采用直接攔截方案成功攔截三維態勢圖Fig.3 The 3D situation map of ATT's successful intercepting under the direct interception scheme

圖4 ATT 采用直接攔截方案成功攔截二維態勢圖Fig.4 The 2D situation map of ATT's successful intercepting under the direct interception scheme

圖5 ATT 采用直接攔截方案未成功攔截三維態勢圖Fig.5 The 3D situation map of ATT's unsuccessful intercepting under the direct interception scheme

圖6 ATT 采用直接攔截方案未成功攔截二維態勢圖Fig.6 The 2D situation map of ATT's unsuccessful intercepting under the direct interception scheme

由圖3～圖6 可知，在直接攔截方案下，ATT 彈道以直線為主，在該方案下ATT 能夠在較短的時間內有效攔截來襲魚雷。但直接攔截對魚雷運動要素的解算要求較高，且追蹤彈道很容易形成尾追態勢而導致攔截失敗，甚至還會對本方艦艇造成傷害（圖5）。

圖7～圖8 所示為第51 次仿真中，ATT 采用迎面逆彈道攔截方案成功攔截來襲魚雷的三維與二維態勢圖。圖9～圖10 所示為第78 次仿真中，ATT采用迎面逆彈道攔截方案未成功規避來襲魚雷的三維與二維態勢圖。

由圖7～圖10 可知，與直接攔截方案相比，迎面攔截方案對魚雷運動要素解算的精確性要求不高，能在一定程度上避免采用解相遇時，因對來襲魚雷態勢估算誤差而導致攔截失敗的問題。但在迎面攔截方案下，ATT 的彈道較為彎曲，需要ATT有更好的法向過載能力，否則，會使ATT 在搜索過程中丟失目標（圖7）。

圖7 ATT 采用迎面攔截方案成功攔截三維態勢圖Fig.7 The 3D situation map of ATT's successful intercepting under the head-on interception scheme

圖8 ATT 采用迎面攔截方案成功攔截二維態勢圖Fig.8 The 2D situation map of ATT's successful intercepting under the head-on interception scheme

圖9 ATT 采用迎面攔截方案未成功攔截三維態勢圖Fig.9 The 3D situation map of ATT's unsuccessful intercepting under the head-on interception scheme

圖10 ATT 采用迎面攔截方案未成功攔截二維態勢圖Fig.10 The 2D situation map of ATT's unsuccessful intercepting under the head-on interception scheme

由圖3～圖10 可知，目前常用的二維艦艇—魚雷對抗模型為簡化的魚雷、反魚雷魚雷運動方程，未考慮其在空間垂直方向上的運動特性，其實質上是以實體運動參數的水平面分量來替代其三維參數，因此帶來了速度、旋回半徑、旋回角速度及航程等方面的誤差。

同時，魚雷在垂直平面上合理的機動體現在二維平面，則表現出轉彎半徑過小，甚至是小于其最小轉彎半徑（圖8）。可見在二維態勢下建立的數學模型無法準確描述三維立體戰場環境下的各實體作戰行為。

基于此，本文建立的三維艦載ATT 對抗魚雷數學模型充分考慮了艦艇、魚雷、反魚雷魚雷等作戰實體在三維的運動特性，能夠更為直觀、逼真地顯示各作戰實體的空間位置關系，且以此建立的對抗誤差模型還增加了對抗模型的精度，彌補了傳統二維模型的不足。

圖11 所示為魚雷報警舷角為120°、ATT 搜索扇面為120°時，不同攔截方案下艦艇單次攔截概率隨ATT 發射距離變化的趨勢圖。

比較迎面攔截與直接攔截這2 種攔截方案，發現當發射距離較近（1 500～3 000 m）時，ATT 采用直接攔截方案要優于迎面攔截方案，這是因為ATT 采用迎面攔截需旋回至來襲魚雷的逆彈道，花費時間較長。但當發射距離較近且ATT 的速度較低時，ATT 很容易在搜索扇面內錯失來襲魚雷，從而導致攔截概率下降。而當發射距離較遠（3 000～6 000 m）時，ATT 采用迎面攔截方案則要明顯優于直接攔截。因此，可根據發射距離的不同來選擇合適的攔截方案。

圖11 不同攔截方案下ATT 發射距離對艦艇單次攔截概率的影響Fig.11 The effects of ATT's transmitting distance to single interception probability under the different interception scheme

4 結 語

本文對傳統的二維艦艇—魚雷對抗模型進行了創新性的改進和提高，提出并建立了三維態勢下的艦載ATT 對抗魚雷數學模型，包括ATT 運動學模型、魚雷運動模型、ATT/魚雷自導模型及誤差模型。同時，基于蒙特卡洛法對ATT 采用直接攔截與迎面攔截這2 種攔截方案進行了仿真分析，定性與定量對比了2 種方案對ATT 攔截概率的影響，得到了有一定價值的結論。該模型可用于ATT 對抗效能的定量分析與評價，從而進一步為ATT 的設計與使用提供參考。

［1］錢東，張起. 歐洲反魚雷魚雷研發展望［J］. 魚雷技術，2006，14（5）：1-5，11.QIAN Dong，ZHANG Qi. Development of anti-torpedo torpedo in Europe［J］. Torpedo Technology，2006，14（5）：1-5，11.

［2］高學強，楊日杰，張林琳. 反魚雷方案對抗效能定量評價模型［J］.探測與控制學報，2010，32（2）：88-91.GAO Xueqiang，YANG Rijie，ZHANG Linlin.Quantita?tive evaluation models for countermeasure effective?ness of anti-torpedo solution［J］. Journal of Detection and Control，2010，32（2）：88-91.

［3］由大德，徐德民.反魚雷魚雷攔截概率影響因素仿真分析［J］.魚雷技術，2010，18（4）：312-315.YOU Dade，XU Demin. Analysis on parameters influ?encing intercept probability of anti-torpedo torpedo［J］.Torpedo Technology，2010，18（4）：312-315.

［4］姜凱峰，周明，葛津華，等.艦載反魚雷魚雷齊射作戰能力研究［J］.艦船科學技術，2010，32（11）：64-68.JIANG Kaifeng，ZHOU Ming，GE Jinhua，et al. Re?search on operational capability of shipboard anti-tor?pedo torpedo［J］. Ship Science and Technology，2010，32（11）：64-68.

［5］王新宏，劉慶，劉建國. 基于相遇區域的反魚雷魚雷攔截彈道［J］.魚雷技術，2013，21（4）：241-245.WANG Xinhong，LIU Qing，LIU Jianguo. Intercept tra?jectory of anti-torpedo torpedo based on meeting area［J］.Torpedo Technology，2013，21（4）：241-245.

［6］周明，徐德民.火箭助飛魚雷彈道的仿真實現與應用研究［J］. 彈箭與制導學報，2006，27（2）：235-238，267.ZHOU Ming，XU Demin. Infection of impact dispersion to the rocket assisted torpedo detection probability［J］.Journal of Projecties，Rockets，Missiles and Guidance，2006，27（2）：235-238，267.

［7］康鳳舉，楊惠珍，高立娥，等. 現代仿真技術與應用［M］.北京：國防工業出社，2006.