劉錦昌,黃樹彩,龐 策,黃 達
(空軍工程大學防空反導學院,陜西西安 710051)
天基信息引導下反導武器系統攔截窗口分析
劉錦昌,黃樹彩,龐 策,黃 達
(空軍工程大學防空反導學院,陜西西安 710051)
基于分析天基信息對反導武器系統的引導方式和引導精度,提出了能夠評估攔截效能的反導武器系統攔截窗口的計算方法;建立了不同預警方式下的預警信息發出時刻和天基信息引導下的武器系統雷達發現目標時間的計算模型。通過仿真表明天基預警信息能夠為反導武器系統提供充足的作戰準備時間,在天基引導信息引導下,武器系統跟蹤制導雷達能夠更準確快速地跟蹤制導雷達發現目標,改善引導信息的精度,增大攔截時間窗口和空間窗口并提高攔截彈攔截更遠射程目標的能力。
天基信息;引導精度;反導武器系統;攔截窗口
及時并準確地感知目標運動信息是反導作戰的先決條件,但是,彈道導彈目標相比于一般的空氣動力學目標具有射程遠,威脅大,速度快,體積小的特點,更加難于探測、識別和攔截[1-3]。當武器系統制導雷達的探測能力無法滿足攔截彈的攔截能力需求時,需要外部信息提供信息支持[4-5]。
天基預警系統具有探測范圍廣,監視空間大,檢測識別準確快速的特點,將天基信息引入作戰系統,為跟蹤制導雷達提供及時準確的引導可有效提高其探測發現目標的能力。在天基信息支持反導作戰的研究方面,文獻[6]采用系統效能分析(SEA)方法分析了天基信息對反導武器系統作戰的支持作用,文獻[7]采用多智能體(MAS)仿真了天基信息支持下的反導作戰系統的組成及流程,文獻[8]基于OODA方法建立了天基信息支持下的反導作戰過程模型,但是這些文獻都是從定性的角度考慮天基信息對反導攔截武器系統的作用,沒有從定量的角度分析天基信息的引導精度對反導攔截武器攔截作戰的影響。
天基預警信息引導對反導攔截系統作戰效能的影響表現在反導攔截武器的攔截窗口[9]大小,天基預警信息提供的引導精度和預警時間會影響攔截彈攔截窗口的大小。本文以反導攔截的時間和空間窗口為天基信息引導作用的指標,以天基預警信息的引導精度和預警時間為輸入,建立天基預警信息引導反導攔截作戰效能評價的方法,分析了天基信息對于反導攔截武器系統的攔截窗口的支持作用。
(1)主動段運動模型
在彈道導彈主動段,采用零攻角模型建立彈道導彈的運動模型如下:
(1)
式中,rm和vm分別表示導彈目標的位置矢量和速度矢量,μ是地球的引力系數,F為導彈推力,m為導彈的質量,β為彈道導彈阻力系數,ρ為空氣密度,Sm為導彈參考面積。在彈道導彈發射參數和初始狀態x0已知的情況下,由運動模型(1)可通過數值積分方法法計算出主動段任意時刻的彈道導彈運動狀態x(t),即
x(t)=f(x0,t)
(2)
(2)預警系統測量模型

(3)
衛星測量方程可表示為
z=[α,e]T+ω=h(xture)+ω
(4)

天基信息對后續武器系統的引導主要采用預測引導的方式進行,目標預測信息精度為天基信息的引導精度。首先,分析天基預警系統對主動段關機點狀態的估計及誤差。在關機點狀態和誤差的基礎上建立預測狀態及誤差的計算方法。
1.2.1目標關機點狀態及誤差
關機點狀態誤差來源有兩方面:一方面是初始狀態估計誤差,另一方面是關機點時間估計誤差。首先分析這兩方面的誤差,然后分析這兩方面的誤差對關機點狀態誤差的影響進而得到關機點狀態的誤差。
(1)導彈初始狀態估計及誤差
假定彈道導彈發射時刻為t0,主動段關機點時刻為tb,預警衛星在導彈穿出云層后首次檢測到目標的時刻為t1,測量周期為Ts,預警衛星在主動段共進行M次測量,測量時刻為ti,i=1,2,…,M,并且tM (5) (6) 彈道導彈初始狀態的最大似然估計為 (7) 式(7)是非線性最小二乘問題,可用Levenberg-Marquart算法求解[11],這里不再贅述。預警衛星初始狀態的極大似然估計誤差為 (8) (2)關機點時間估計誤差 關機點時刻tb無法由衛星準確觀測到,通過估計得到,預警衛星最后一次觀測到目標的時刻為tM,且tM (9) (3)關機點狀態及誤差 (10) 1.2.2天基信息導引誤差 (11) 目標在RSW坐標系下跟據關機點時刻信息預測t時刻的狀態誤差協方差為PRSW(t),令J4是xRSW(tb)到xRSW(t)的狀態轉移矩陣,則 (12) (13) 要計算攔截彈發射時間窗口和攔截空間窗口,首先要對反導攔截的作戰時序進行分析,圖1 給出了反導攔截作戰的作戰時序。 圖1中,twarn是預警信息發出時刻,反導武器系統收到預警信息后開始從戰備狀態轉為戰斗狀態,從戰備狀態轉為作戰狀態的時間稱為戰斗準備時間Tzh。 2.1.1反導攔截彈發射時間窗口 反導攔截彈發射窗口根據來襲彈道目標彈頭運動軌跡和反導攔截武器系統的作戰性能計算出的攔截彈攔截來襲目標的發射時間區間,攔截彈發射時間窗口為[tlaumin,tlaumax],其計算過程如下。 (1)計算攔截彈最早發射時刻tlaumin。攔截彈在最遠殺傷區殺傷目標的發射時刻為tfar,反導武器系統達到戰斗狀態的最早時刻為tzh=twarn+Tzh,攔截彈發射必須同時滿足跟蹤制導雷達發現目標、彈目遭遇點在殺傷區內、武器系統處于戰斗狀態三個條件,因此,攔截彈最早發射時刻tlaumin可由下式計算 tlaumin=max{tzh,tfar,tdis} (14) (2)計算攔截彈最晚發射時刻tlaumax。攔截彈在最近殺傷區殺傷目標的發射時刻為tnea,當tdis>tnea或tzh>tnea時,說明跟蹤制導雷達發現目標過晚或者提供的預警時間太晚不足以準備攔截,因此攔截彈不具備發射設時間窗口,無法攔截目標。只有當tdis tlaumax=tnea (15) (3)二次攔截發射條件分析。攔截彈第一次發射時刻為tlau1=tlaumin,第一次遭遇時刻為tint1,攔截后效能評估時間為tpg,則第二次最早發射時刻為tlau2=tint1+tpg,如果tlau2 2.1.2反導攔截彈攔截空間窗口 空間窗口是指攔截彈攔截來襲目標的空間遭遇點的高度區間與遠程區間。攔截彈在tlaumin時刻發射,彈目遭遇點位置的高度為hmax,與發射點的距離為dmax。攔截彈在tlaumax時刻發射,彈目遭遇點高度為hmax,與發射點的距離為dmax,反導攔截作戰的高度窗口為[hmin,hmax],距離窗口為[dmin,dmax]。 2.2.1跟蹤制導雷達發現目標時間計算模型 跟蹤制導雷達發現目標的過程分三步,首先確定雷達的截獲窗口,接著雷達向截獲窗口內進行波束掃描,最后發射若干脈沖對目標進行檢測。記雷達對目標的發現概率為Pa,Pb表示目標處于截獲窗口的概率,Pc表示目標被掃描波束截獲的概率。截獲窗口表示雷達搜索的空域,Pb定義為目標落入雷達截獲窗口Ω1中的概率[5],Ω為屬于Ω1中的任意區域,用f(Ω)表示在區域Ω內目標存在的概率密度,用sΩ表示區域Ω的面積大小,當Ω足夠小時,Pb的表達式為 (16) 1-P[(tc (17) 根據德摩根定律可得 (18) 由于在[tmin,t]內沒有發現目標和在[t,t+Δt]內沒有發現目標這兩事件相互獨立,所以有 (19) 令單位時間目標被掃描波束截獲概率為Ps,則有 (20) (21) 整理可得 (22) 當Δt→0時,因為雷達探測具有連續性,則有 (23) (24) 積分可得 (25) 波束駐留時間為Tdwell,雷達每次波束搜索的面積為S,則雷達單位時間內對雷達截獲窗口的掃描次數為 (26) 目標被雷達單個波束檢測出來的概率為Pdect,因此目標單位時間內被截獲的概率為 Ps=1-(1-Pdect)ns (27) Pdect與雷達距離目標的距離有關,如果雷達最遠探測距離處的單波束檢測概率為Pdect0,最大探測距離為rdect0,則雷達在距離目標為r處的單波束檢測概率為 (28) 式中:Pfa為雷達探測虛警概率,目標與雷達的距離r只是時間的函數,因此Pdect也只是時間的函數。 雷達發現目標的概率為 Pa(t)=PbPc(t) (29) 發現目標的概率的閾值設為PT,目標被制導雷達發現的時刻為tdis,則Pa(tdis)需滿足Pa(tdis)≥PT,即tdis為滿足制導雷達發現目標概率閾值的最早時刻。如果雷達開始探測目標的時刻為tsta,那么雷達發現目標所用時間為Ta=tdis-tsta。下面分別分析無天基信息引導和有天基信息引導情況下計算跟蹤制導雷達發現目標時間的所需參數。 (1)無天基引導信息 跟蹤制導雷達沒有外部信息進行引導提示時,只能采用全域掃描的方式進行目標搜索,截獲窗口為雷達的整個可觀察區域,目標在可觀察區域內服從均勻分布,Pb可由下式計算 (30) S總為目標可能出現的區域面積,S可視是雷達可觀察總面積,當目標進入雷達探測范圍時,雷達開始探測目標,tsta=tmin。 (2)有天基引導信息 (31) (32) 慣性坐標系下目標位置預報誤差協方差為PECI,r,雷達極坐標系下的預報位置誤差為PR,r,則 (33) (34) 引導信息提供給制導雷達的時刻為tyin,tyin 2.2.2預警信息發出時刻計算模型 反導系統的預警信息可由天基預警衛星提供,也可由遠程預警雷達提供,下面分別分析由天基預警系統和遠程預警雷達提供預警信息的預警信息發出時刻計算模型。 (1)天基預警衛星提供預警信息 天基預警系統在第一次觀測到導彈尾焰后,經過四次連續觀測確認目標,再通過估計可初步得到導彈的射程、射向、類型等信息,并提示反導攔截武器系統開始進行戰斗準備,信息處理時間為tps,由前文可知預警衛星首次探測到目標時刻為t1,天基預警系統的預警信息發出時刻為 twarn=t1+4Ts+tps (35) (2)遠程預警雷達提供預警信息 (36) 當目標從搜索遠界進入雷達探測范圍時,雷達的預警信息發出時刻為 (37) 預警衛星運行軌道為地球同步軌道,衛星觀測周期Ts=2 s,衛星觀測視線角誤差為σLOS=40 μrad,云層高度設為7 km,預警信息處理時間為5 s。預警雷達的波束寬度為1°,波束駐留時間為16 ms,本文設自由飛行段的彈頭的RCS為0.25 m2,預警雷達對RCS為0.25 m2的目標的最遠搜索距離為2617 km,雷達虛警概率為10-7,最遠搜索距離處單波束檢測概率為0.8,觀察區域方位角為-60°~+60°,俯仰角為3°~85°,預警信息處理時間為10 s。THAAD[1]系統是較為典型的反導攔截武器系統,其性能參數如表1所示。預警雷達、跟蹤制導雷達部都署在目標落點位置。 表1 THAAD系統性能參數Table 1 The performance parameters of THAAD system 隨著彈道導彈的射程不同,其彈道總飛行時間、導彈出云時刻會有很大的差別,以目標發射時刻為零時刻,不同射程的目標進入預警雷達探測范圍的時刻及其與預警雷達的距離如表2所示。 從表2可以看出當到導彈的射程小于3500 km時,預警雷達可以在搜索低界探測到目標,隨著射程的增加,目標進入預警雷達探測范圍的時間會增大,當導彈射程為4000 km時,預警雷達無法在低界搜索到目標,只能在遠界最早搜索到目標,導致目標進入預警雷達探測范圍的時間大大增加,目標進入預警雷達探測范圍后被預警雷達發現的概率隨時間變化的曲線如圖3所示。 表2 不同射程導彈飛行時間參數Table 2 The flight time parameters of different range missle 圖3中橫坐標以目標進入預警雷達探測范圍為零時刻,目標被雷達發現的概率閾值設為0.95,從圖上可以看出,目標被雷達發現的概率隨時間的增加而增加,隨著目標射程的增加,預警雷達發現目標所用時間會增大,這是由于目標距離雷達越遠,單波束檢測概率越低所導致的。天基衛星和遠程預警雷達對不同射程彈道導彈預警信息發出時刻如圖4所示。圖中縱坐標以彈道導彈發時刻為零時刻,可以看出,天基衛星的預警信息發出時刻明顯早于預警雷達,這是因為天基衛星的軌道高,探測范圍大,能夠在導彈出云時刻附近迅速捕獲到目標,尤其是當導彈的射程增加時,天基衛星信息預警時刻幾乎不受影響,而預警雷達受發現目標時間及地球曲率的影響,導彈射程增加后,預警信息發出時刻大大增加。 無信息引導時,跟蹤制導雷達目標發現概率隨目標進入雷達探測范圍時間的變化曲線如圖5所示。 圖中5橫坐標以目標進入跟蹤制導雷達探測范圍為零時刻,不同射程的彈道導彈在跟蹤制導雷達探測范圍內飛行的時間不同,因此每條曲線的結束時間不同,本文以目標落地時刻為仿真結束時刻,從圖中可看出在沒有天基信息對制導雷達進行引導時,跟蹤制導雷達在目標落地前無法發現目標,攔截系統將不具備攔截能力,因此,必須在外界信息的引導下反導攔截系統才能發揮攔截效能。采用天基信息引導下對不同射程的目標引導位置誤差大小、跟蹤制導雷達截獲窗口大小及跟蹤制導雷達發現目標所用時間如表3所示。 從表中數據分析可得,天基信息引導的位置誤差隨著目標射程的增加而增加,跟蹤制導雷達截獲窗口大小隨著引導誤差的增大而增大,進而導致目標發現時間的增大。 攔截系統對不同射程目標的攔截發射時間窗口如圖6所示。 表3 天基信息引導下制導雷達發現目標時間Table 3 The time of guidance radar using to discover target under the guidance of space-based information 圖中縱坐標以目標發射時刻為零時刻,從圖6可以看出,采用遠程預警雷達進行預警,當目標的射程較短時,反導攔截武器系統達到戰斗狀態的最早時刻tzh晚于tfar,從而導致攔截系統的最早發射時刻變晚,時間窗口變小。引入天基信息進行預警時,反導攔截武器系統達到戰斗狀態的最早時刻tzh遠遠小于tfar,可為攔截系統提供更充分的作戰準備時間。天基信息發出預警后,當目標的射程小于3000 km時,制導雷達發現目標時刻tdis小于tfar,可以在攔截彈遠界攔截目標;當射程大于等于3000 km時,由于引導信息位置誤差的增大,目標發現時間的增長而導致制導雷達發現目標時刻tdis大于tfar,攔截時間窗口變小。通過以上分析可知,天基信息引導下的制導雷達發現時刻與引導信息的精度密切相關,此時的引導信息精度無法滿足攔截系統攔截射程大于3000 km的目標的需求,需要提高天基信息引導精度,而預測信息的精度主要與天基預警系統的觀測視線角誤差和觀測周期有關,下面仿真不同觀測條件下的引導效果,仿真結果如圖7所示。 表4 觀測條件參數設置Table 4 The observation condition parameters setting 觀測條件3000km射程目標3500km射程目標4000km射程目標距離窗口高度窗口距離窗口高度窗口距離窗口高度窗口條件1[90,177][40,139][90,108][40,51]——————條件2[90,200][40,150][90,167][40,131]——————條件3[90,200][40,150][90,200][40,150][90,152][40,117]條件4[90,200][40,150][90,200][40,150][90,121][40,103]條件5[90,200][40,150][90,200][40,150][90,200][40,150] 仿真結果表明,通過減小觀測周期和觀測視線角誤差可以改善天基信息的引導效果,縮短跟蹤制導雷達發現目標時間,使攔截系統攔截空間窗口變大,能夠更早發射攔截彈并在更遠的地方攔截目標,并且能夠攔截射程更遠的目標。 本文建立了信息引導下反導武器系統攔截窗口的計算模型,分析了天基信息對反導攔截作戰支持作用。結果表明:1)天基預警系統相對于遠程預警雷達,能夠更早發現目標,為反導武器系統提供充足的作戰準備時間;2)反導武器系統獨立作戰時,沒有攔截發射窗口,不具備攔截作戰能力,當有天基引導信息提供引導時,跟蹤制導雷達能夠準確快速地發現目標,具備攔截一定射程內的彈道導彈的能力;3)改善天基預警系統的觀測條件可以提高天基信息的引導效果,增大攔截時間窗口和空間窗口,提高攔截彈攔截更遠射程目標的能力。 [1] Gilmore J M. 2014 assessment of the ballistic missile defense system [R]. Washington, DC: National Defense Authorization Act. March 2015. [2] 許登榮, 程水英, 包守亮. 一種多級助推段彈道導彈跟蹤算法[J]. 宇航學報, 2016, 37(9): 1114-1121. [Xu Deng-rong, Cheng Shui-ying, Bao Shou-liang. A tracking algorithm of ballistic missile with multi-stage boosters [J]. Journal of Astronautics, 2016, 37(9):1114-1121.] [3] 申鎮, 強勝, 張寅生.單星無源探測彈道導彈射向估計新方法[J]. 宇航學報, 2011, 32(7): 1451-1456. [Shen Zhen, Qiang Sheng, Zhang Yin-sheng. New method on course heading estimation with single satellites passive detection [J]. Journal of Astronautics, 2011, 32(7): 1451-1456.] [4] 李志匯, 劉昌云, 倪鵬, 等. 反導多傳感器協同任務規劃綜述[J]. 宇航學報, 2016, 37(1): 29-38. [Li Zhi-hui, Liu Chang-yun, Ni Peng, et al. Review on multisensor cooperative mission planning in anti-TBM system [J]. Journal of Astronautics, 2016, 37(1): 29-38. ] [5] 唐書娟, 許蘊山, 畢篤彥. 信息引導條件下雷達搜索空域及策略[J]. 西安電子科技大學學報(自然科學版), 2016, 43(1): 173-179. [Tang Shu-juan, Xu Yun-shan, Bi Du-yan. Airspace and searching strategy of radar under the information guiding condition [J]. Journal of Xidian University, 2016, 43(1): 173-179.] [6] 黃樹彩, 周延延, 韋剛. 基于天基信息的反導作戰效能分析[J]. 系統工程與電子技術, 2009, 31(10): 2414-2417. [Huang Shu-cai, Zhou Yan-yan, Wei Gang. Operation effectiveness analysis of anti-TBM system with space-based information support [J]. Systems Engineering and Electronics, 2009, 31(10): 2414-2417.] [7] 楊娟, 池建軍. 基于 MAS 的空間信息支援反導作戰效能評估[J]. 兵工自動化, 2013, 32(7): 9-12.[Yang Juan, Chi Jianjun. Evaluation of space information support anti-TBM operation efficiency based on MAS[J]. Ordnance Industry Automation, 2013, 32(7): 9-12.] [8] 吳忠望, 張漢峰, 姚宏林. 基于OODA 環的天基信息支援防空反導作戰建模[J]. 四川兵工學報, 2013, 34(9): 8-16.[Wu Zhong-wang, Zhang Han-feng, Yao Hong-lin. Research on the operation model of air defense and anti-missile under space-based information support with OODA loop [J]. Journal of Sichuan Ordnance, 2013, 34(9): 8-16.] [9] 邵建兆, 畢義明, 錢大慶. 雙層火力協同反導作戰中的火力調度[J]. 系統工程與電子技術, 2016, 38(8): 1840-1846. [Shao Jian-zhao, Bi Yi-ming, Qian Da-qing. Firepower scheduling research on double-layer firepower coordination anti-missile combat [J]. Systems Engineering and Electronics, 2016, 38(8): 1840-1846.] [10] 盛衛東, 徐洋, 周一宇. 天基光學傳感器的視線測量誤差分析[J]. 宇航學報, 2011, 32(1): 129-135. [Sheng Wei-dong, Xu Yang, Zhou Yi-yu. Analysis of LOS measurement error for space-based optical sensor[J]. Journal of Astronautics, 2011, 32(1): 129-135.] [11] 張峰, 田康生. 預警衛星對彈道導彈主動段狀態估計算法[J]. 紅外與激光工程, 2013, 42(7): 1912-1922. [Zhang Feng, Tian Kang-sheng. Estimation algorithm of boost-phase ballistic missile state from satellite observations [J]. Infrared and Laser Engineering, 2013, 42(7): 1912-1922. ] [12] 白顯宗, 陳磊, 張翼. 空間目標碰撞預警技術研究綜述[J]. 宇航學報, 2013. 34(8): 1027-1039. [Bai Xian-zong, Chen Lei, Zhang Yi. Survey on collision assessment and warning techniques for space object [J]. Journal of Astronautics, 2013, 34(8): 1027-1039.] [13] 白顯宗, 陳磊. 天基空間監視交會計算和可觀測時段預報的誤差分析[J]. 宇航學報, 2015, 36(5): 574-582. [Bai Xian-zong, Chen Lei. Error Analysis of conjunction calculation and visible time duration prediction for space-based space surveillance[J]. Journal of Astronautics, 2015, 36(5): 574-582.] [14] 白顯宗, 陳磊, 唐國金. 碰撞概率對軌道誤差和接近幾何條件的靈敏度分析[J]. 宇航學報, 2016, 37(4): 427-434. [Bai Xian-zong, Chen Lei, Tang Guo-jin. Sensitivity of collision probability to orbit error and close approach geometry [J]. Journal of Astronautics, 2015, 37(4): 427-434.] AnalysisofInterceptWindowofAnti-MissileSystemundertheCueingofSpace-BasedInformation LIU Jin-chang, HUANG Shu-cai,PANG Ce, Huang Da (Air and Missile Defense College, Air Force Engineering University, Xi’an 710051, China) Firstly,the cueing method and cueing accuracy of the space-based information for anti-missile weapon system are analyzed. Then,the method of calculating the intercept window of anti-missile system which can evaluate the interception efficiency is put forward. The calculation model of the early warning information issuing time with different warning mode and the time that guidance radar need to find target under the cueing of space-based information is established. Lastly, the simulation results indicate that the space-based warning information can provide sufficient preparation time for the anti-missile system and the guidance radar cued by the space-based information can be faster and more accurate to find the target. Increased accuracy of cueing information can increase the intercept time window and space window, and raise the ability of the interceptor to intercept farther range target. space-based information;cueing information accuracy;anti-missile system;intercept window 2017- 07- 12; 2017- 09- 20 國家自然科學基金(61573374);青年國家自然科學基金(61503408) TJ761.7 A 1000-1328(2017)12- 1339- 09 10.3873/j.issn.1000- 1328.2017.12.011 劉錦昌(1993-),男,碩士,主要從事空天目標協同跟蹤與攔截引導研究。 通信地址:陜西省西安市長樂東路甲字一號(710051) 電話:15129650516 E-mail:15129650516@163.com



2 天基信息引導下攔截窗口模型
2.1 攔截彈發射時間窗口與空間窗口模型
2.2 時間參數計算模型














3 仿真校驗
3.1 參數設定

3.2 預警信息發出時刻

3.3 跟蹤制導雷達發現目標所用時間
3.4 天基信息引導對攔截窗口分析



4 結 論