印度反衛事件分析

王偉林，何海峰，楊永嘉，周士豪

（航天工程大學，北京101400）

0 引言

世界協調時(UTC)2019年3月27日，印度反衛星試驗成為世人關注的焦點。根據國外報道以及衛星軌道數據的計算結果可以確認，撞擊時間為2019年3月 27日 5時 42分（UTC），撞擊地點為（17.68°N，87.65°E），高度約為 283 km，相 對 撞 擊速度約為9.8 km/s，目標衛星是印度國防研究與發展組織（DRDO）生產研制的“Microsat-R”衛星。

1 試驗中目標衛星與攔截彈分析

1.1 目標衛星分析

Microsat-R衛星于2019年1月24日晚11點37分在斯里赫戈達島蒂什·達萬航天中心搭載PSLV C-44進入高度為277.2 km的軌道，該衛星是由印度國防研究與發展組織（DRDO）實驗室開發，在發射時印度聲稱此衛星是一顆對地觀測衛星，但是歐洲咨詢公司在2015年發布的《全球政府航天計劃預算分析》表明，印度對地觀測衛星絕大多數由印度空間研究組織（ISRO）研制（見表1），而Microsat-R并不在此計劃之列。結合人造衛星造價昂貴且Microsat-R衛星僅在太空中存在61天的事實，可以推測出Microsat-R衛星可能本來就是作為反衛試驗的靶星發射，印度對于此次反衛星試驗計劃已久。

1.2 攔截彈分析

印度公開的紀錄片顯示，此次試驗中的攔截彈的正式名稱為：PDV-MK-II。它很可能是2014年4月27日DRDO首次進行大氣層外反導試驗的PDV（“大地防御飛行器”）攔截彈的改進型，其主要性能指標如表2所示。PDV的攔截高度小于150 km，如果不加改進，它的攔截高度達不到300 km，因此認為PDV-MKII導彈可能使用了新的助推段。根據網站公布的信息，印度的反衛星導彈配備了兩級固體火箭發動機和一個攔截彈頭，這種導彈能夠捕獲高達1000 km的目標。通過芯級側面向下延伸的線槽，可以得知在PDV衍生的第三級之下有兩級發動機，線槽中導線將控制命令從最頂級傳送到第一級的推力控制系統。線槽的中斷標志著第一級和第二級之間的級間段。通過此次反衛星導彈的照片和之前PDV的照片比對(見圖1)可以推測，在彈徑1.1 m的情況下，印度反衛星導彈的長度大約在11 m左右。如果印度反衛導彈的質量真的如印度媒體所言達到18 t，那么相比于美國彈徑1 m、彈長 10.6 m、質量不到 7.5 t的“潘興”-Ⅱ，印度的反衛導彈顯得過重，說明印度在彈體結構和材質上還有很大差距，因此，可推測：印度此次反衛星導彈的攔截高度不會高出300 km太多，達到印度聲稱1000 km可能性更小。

表2 印度PDV攔截彈主要性能指標

印度DRDO公開的紀錄片顯示，該彈頭采用了末端紅外制導的方式，Microsat-R衛星作為近地軌道的對地觀測衛星，太陽能電池板通常指向太陽，在攔截彈末制導過程中，太陽能電池板不應該是完全可見的。利用STK對此次碰撞進行了模擬，若如印度所言的采用頂點攔截，則碰撞前目標衛星與攔截彈太陽電池板的相對位置關系應如圖2所示。

而在DRDO公開的視頻資料中，衛星的太陽帆板面向彈頭，且面積張開幾乎最大，如圖3所示。

這與理論情況不符。所以推測這可能是由以下兩種情況造成，一是在碰撞發生前，該衛星太陽帆板進行了有目的性的姿態調整以增強紅外特性，如圖4所示，側面證明目標衛星在此次打擊過程中可能是合作靶星。

表1 印度后續對地觀測衛星發展計劃

圖1 印度反衛攔截彈與PDV彈等比比對圖

圖2 理論碰撞時衛星姿態示意

圖3 DRDO視頻顯示導彈上攝像頭拍攝目標衛星太陽能電池板的最后圖像

圖4 預測情況一碰撞時衛星姿態示意

二是此次反衛實驗采用的是上升段攔截方式，而與印度所宣稱的頂點攔截方式不符，如圖5所示。

圖5 預測情況二碰撞時衛星姿態示意

2 模型校對與驗證

2.1 模型來源及實驗目的

利用本課題組的航天器解體模型碎片生成軟件、碎片演化分析軟件對此次碰撞事件進行情況再現和碎片演化分析，仿真模擬出印度反衛星試驗可能的打擊方式，對本課題組模型的準確性進行校對、驗證，進而預測印度此次試驗產生的碎片在軌數量及其變化趨勢，分析印度反衛星試驗對空間環境造成的威脅。

2.2 航天器解體模型碎片生成軟件、碎片演化分析軟件簡介

2.2.1 碎片生成軟件的功能

1）解體程度分析：根據碰撞時刻兩物體的運動參數，評估由于碰撞造成航天器解體的程度，是完全解體還是部分解體。2）解體碎片生成：根據航天器解體模型，抽樣計算碰撞后所產生的每個碎片的特性參數，包括質量、尺寸、面質比和飛散速度。3）碎片分布統計：對所有生成的碎片分別統計質量、尺寸、面質比和飛散速度等特性的分布。4）結果輸出與顯示：將產生的碎片數據按照一定的格式進行存儲，對于碎片統計分布則可以曲線形式進行顯示。

2.2.2 碎片演化軟件的功能

1）碎片在軌分析：根據輸入的碎片數據文件，利用空間物體受力模型，對在軌碎片進行演化，計算出其留軌數量、質量的變化。2）結果顯示：根據動態演化時長，以圖形界面對空間碎片的位置、運動狀態進行動態顯示，并實時顯示留軌數量、留軌質量。

2.3 參數定義

在仿真模擬中，主要關心碰撞發生時衛星和攔截彈的相對位置、相對速度對碰撞碎片生成、演化的影響。對此首先對坐標參數進行定義說明，以目標衛星為原點，建立目標衛星軌道坐標系，如圖6所示，原點O1為目標衛星與攔截彈碰撞發生位置，其中X軸表示為目標運動速度矢量方向，Z軸表示為地球矢徑方向，Y軸根據右手定則確定。碰撞時刻攔截彈速度方向可在該坐標系下進行分解，在平面XO1Y投影，與X軸夾角定義為方位角α，在平面XO1Z投影，與X軸夾角定義為航跡角β。當β=0時，即可表示頂點攔截；當0＜β＜90°時，定義該攔截為上升段攔截；當-90°＜β＜0°時，定義為下降段攔截。當 α=180°，β=0°時，即表示頂點逆軌攔截，如圖 7所示；當 α=0°，β=0°時，即表示頂點順軌攔截，如圖8所示。

2.4 試驗初始參數設計

由于缺乏印度反衛星試驗中Microsat-R衛星與攔截彈的具體詳細信息，僅能根據有關推測及碰撞后的初步探測結果對撞擊場景做如下想定：

對于目標衛星軌道信息，為盡可能接近碰撞時刻衛星運動狀態實際，采用美國www.spacetrack.com網站在公布的世界協調時2019年3月26日（碰撞發生前一天）更新的Microsat-R兩行TLE數據，由于人造地球衛星與空間探測器的結構材料大多采用鋁合金，假設衛星材質為鋁合金，根據衛星的材質和質量、有關圖像資料（見圖9），假設衛星為標準立方體，估計出其尺寸。

圖6 目標衛星軌道坐標系

圖 7 頂點逆軌攔截（α=180°,β=0°）

圖 8 頂點順軌攔截（α=0°,β=0°）

圖9 Microsat-R衛星

1）目標衛星TLE軌道根數

0 MICROSAT-R

1 43947U 19006A 19084.98720666.00152043 11367-4 27860-3 0 9992

2 43947 96.6362 358.3958 0015832 220.6056 155.1371 16.01110841 9636

2）材質：鋁合金；

3）質量：740 kg；

4）尺寸：72 cm×72 cm×65 cm；

5）撞擊時刻：2019年3月27日05∶42(UTC)；

6）笛卡爾J2000坐標系下撞擊時刻的軌道參數，如表3所示。

表3 撞擊時刻的軌道參數

如果印度此次使用的反衛戰斗部為動能攔截方式，那么在行動中使用的反衛星武器類型可能與美國洲際彈道導彈防御實驗地面攔截導彈（GBI）攔截器相似。由于缺乏對印度反衛星試驗所使用攔截彈的確切情報，在此假設此次試驗使用的攔截彈尺寸、質量與其均大致相似。

相關資料表明，美國洲際彈道導彈防御實驗地面攔截導彈（GBI）攔截器為EKV。在此暫用維基百科對該地面攔截導彈（GBI）攔截彈的描述進行仿真實驗。

EKV參數：質量約為64 kg；長度為1.4 m；直徑為0.6 m。

2.5 實驗過程仿真模擬

對印度反衛星試驗的具體細節進行研究，碎片產生數量、在軌時間可能與攔截彈與目標衛星的相對速度、角度等相關因素有關，由于缺乏對碰撞時刻攔截彈的位置、運動狀態的相關具體情報，在已知目標衛星速度、角度的前提下，在仿真模擬過程中利用筆者課題組航天器解體模型、演化模型，假定目標衛星、攔截彈的相對方向（α、β）、速度大小，并設置多組攔截彈速度、角度等參數，對攔截彈分別在不同的的撞擊速度、角度進行了仿真模擬，分析其碎片生成數量、質量、留軌時間，得出有關結論。

2.5.1 實驗過程

根據以上假定設置碰撞的初始參數數據，結合速度、角度的改變，計算出模型應用的瞬時速度、位置參數，輸入到碎片生成軟件中得出，統計生成碎片數量、質量、面質比等參數，統計結果如表，而后利用碎片演化軟件，對攔截彈速度為2.5 km/s、生成尺寸大于5 cm的碎片的演化進行分析，仿真其在軌數量、質量隨演化時長變化情況，統計結果如表4所示。

2019年5月19日美國SpaceTrack網站公布印度反衛星試驗碎片數量為92個(其中在軌63個，墜落29個)，利用本課題組航天器解體模型進行模擬仿真，在10 cm、5 cm碎片尺寸以上的數量與美國公布數據相比相對較多，主要原因可能是由于碰撞軌道較低，碰撞發生初始階段碎片數量下降較快，一些速度較低的碎片在被NASA探測到之前就已經墜毀，可以認為此模型較為準確，本次仿真結果與印度反衛星試驗實際相差不大。

假定美國低軌空間碎片觀測能力在5 cm左右，為了進一步驗證模型的準確性，在仿真實驗中對5 cm以上碎片的演化進行了模擬，表5為攔截彈V=2.5 km/s時，尺寸大于5 cm的碎片的演化結果。

表4 碎片的生成

表5 攔截彈V=2.5 km/s時，尺寸大于5 cm的碎片的演化

如表5所示，在碰撞發生24 h內，碎片大量墜落，數量、質量迅速下降，大于5 cm的碎片數量基本穩定在30～110之間，與美國公布數據數量級基本一致，進一步驗證解體模型對碎片數量預報的準確性。截至2019年05月15日編目的在軌碎片數量為63，遠大于表5中1～4組的數量，因此傾向認為印度反衛星試驗更可能采取了上升段攔截的方式，其攔截角度可能在40°～60°之間。

2.5.2 結論分析

1）表第 1～3組數據中，第一組 α=180°，β=0°時，即頂點逆軌攔截時，此時目標衛星和攔截彈的相對速度最大，產生的碎片數量最多，第4組當α=60°，β=0°時，此時相對速度較小，產生的碎片數量最少，可知頂點攔截時（即β=0°），其他變量不變的情況下，當α越大時，相對速度越大，產生的碎片越多，即順軌攔截相對于逆軌攔截生成的碎片數量更少。在表5演化實驗結果中，第1、2、3組表明，逆軌攔截雖然產生碎片數量較多，但其碎片留軌時間較短，對空間環境的影響更小。

2）對比表 4第 1、11、12組數據，當 α、β一定時，改變碰撞發生時刻攔截彈的速度，對于統計大于0.1 cm以上碎片，攔截彈速度越小，即攔截彈與目標衛星的相對速度越小，產生的碎片數量越少.

3）對比表4第5～10組數據，碎片在較大尺寸（10 cm以上）時，β對碎片的產生數量影響不明顯，而在小尺寸（0.1 cm以上）時，無論是上升段攔截還是下降段攔截相對于頂點攔截都有利于減少碎片產生數量，且角度越大，對減少碎片數量的影響越強烈；而β＞0（上升段攔截）相對于β＜0（下降段攔截）時，碎片數量減小更加明顯。在表5演化實驗結果中第1、4、5、6、7組表明，頂點攔截產生的碎片留軌時間明顯相對更短；下降段次之，且角度越大，留軌時間越短；上升段產生的碎片留軌時間最久，且角度越大，留軌時間更久。總體而言，逆軌頂點攔截方式對空間環境的影響更小。

3 碰撞碎片對空間安全的影響

印度反衛試驗產生的碎片分布在200～2500 km（遠地點）的范圍內，其中大部分碎片集中在300～800 km的軌道高度上。在這個范圍的軌道上，運行著大量的低軌道衛星，碎片對衛星的運行產生了巨大的威脅。例如，在遠地點高度為409 km左右的近地軌道上，有著人類寶貴的太空資產——國際空間站，而在403 km（遠地點）和415 km（遠地點）的軌道上存在著2個碰撞碎片，給國際空間站的運行帶來了巨大的風險。

根據SpaceTrack網站(2019-04-18)更新的數據，如圖10、11所示，印度反衛試驗所產生的碎片被SpaceT-rack編目的有59個，其軌道傾角主要集中在94°～97°范圍內，軌道高度在300～800 km范圍內的有37個，在300 km以下的有3個，最高的軌道高度為2226.8 km，對該范圍內在軌的航天器構成了巨大的威脅。

圖10 碎片遠地點分布（2019-04-18）

圖11 碎片軌道傾角分布（2019-04-18）

根據SpaceTrack網站(2019-12-23)更新的數據，如圖12、13所示，時隔9個月后，印度反衛試驗所產生的碎片被SpaceTrack編目的仍然有23個，軌道高度在300 km以上的有16個，最高的軌道高度為1766.7 km，碎片在軌時間遠遠超過其宣稱的“45天內消失”，雖然目前沒有任何消息披露在軌衛星受到了印反衛試驗的影響，但碎片對該范圍內在軌的航天器依然具有潛在威脅。

圖12 碎片遠地點分布（2019-12-23）

圖13 碎片軌道傾角分布（2019-12-23）

4 結束語

印度反衛試驗產生了大量空間碎片，對低地球軌道運行的航天器造成了長期的安全威脅，對人類航天活動將產生深遠的影響，引起了國際宇航界和國際社會的高度關注和普遍擔憂，印度借發展反導發展反衛星能力，嚴重威脅太空安全態勢。近一年以來，對于印度此次反衛試驗，美國等國僅僅是稍加指責，其他國家也沒有更大的反應，印度幾乎沒有遭到處罰，這很可能加劇各國在反衛星武器研發領域的競賽，未來很可能有更多的國家擁有同樣的武器，更是在國際上埋下了一個重大隱患。

本文通過模擬實驗對該次碰撞事件進行情況再現、碎片演化分析，進一步驗證了本課題組模型的準確性。不同撞擊方式對空間環境的威脅程度不同。逆軌攔截相對于順軌攔截產生的碎片留軌時間短，能夠在較短時間內大量墜落。頂點攔截、下降段攔截和上升段攔截在碎片產生數量上差別不大，但頂點攔截碎片最易墜落。總體而言，頂點逆軌攔截對環境的影響最小。在模擬仿真中，對于較大尺寸的碎片生成、演化數量統計結果，本課題組模型模擬結果與NASA觀測實際基本一致，具有較高的準確性；但對于更小尺寸的碎片，由于缺乏有效實際數據支撐，暫時無法進行對比驗證。通過碎片數量的分析以及紅外制導圖像的分析，傾向認為印度反衛試驗采用了上升段攔截，碎片在軌時間長，嚴重影響空間環境安全。但仿真實驗基于開源資料以及假定分析，缺乏印度反衛星試驗的確切詳細情報，還需進一步跟蹤研究。■