葉楠，王帆，馮軍華，李海玥

(宇航動力學(xué)國家重點實驗室，西安710043)

1 引言

太空數(shù)以億計的空間碎片，嚴重威脅到在軌航天器的安全。為減小空間碎片的危害，針對厘米級以下空間碎片，可采用結(jié)構(gòu)防護方式等；10cm以上的空間碎片，航天器可通過主動變軌來進行規(guī)避；厘米級空間碎片目前由于缺乏有效防護手段，危害最大。有代表性的非激光空間碎片清除方法有太陽能帆、電動力纜繩、阻力增加裝置、網(wǎng)捕裝置和軌道拖船等，這些技術(shù)要么每次清除碎片數(shù)目有限，成本太高；要么仍然處于概念研究階段，技術(shù)尚未成熟。地基激光主動清除碎片技術(shù)由于其優(yōu)良的實用性和部署運營成本低廉，最具應(yīng)用前景[1]，受到廣泛關(guān)注。

1996年，在美國國家航空航天局 (NASA)和美國空軍的資助下，美國洛斯阿拉莫斯實驗室Phipps等人提出利用地基高能脈沖激光清除近地軌道空間碎片的獵戶座 (ORION)計劃[2]，如圖1所示。該計劃對地基激光清除空間碎片機理和環(huán)節(jié)進行了系統(tǒng)的分析論證，形成了包含系統(tǒng)設(shè)計、目標獲取、方案風險等內(nèi)容的系統(tǒng)報告。

圖1 ORION計劃示意圖Fig.1 Schematic diagram of ORION

美國國防支援計劃 (DSP)一期由于工作波長選擇不當，導(dǎo)致衛(wèi)星在軌部署之后，虛警概率過高而無法使用，說明在工程論證階段，參數(shù)設(shè)計和復(fù)核的重要性。本文擬就地基激光清除碎片的原理和參數(shù)設(shè)計依據(jù)，對ORION計劃公布的設(shè)計參數(shù)進行復(fù)核，檢驗其合理性。

2 關(guān)鍵技術(shù)

地基激光清除碎片原理為:地基高能激光輻照空間碎片，空間碎片表面形成等離子體反噴，產(chǎn)生沖量，降低碎片速度和軌道高度，使其近地點墜入大氣層最終燒毀。由設(shè)計概念可知，地基激光清除空間碎片，必須有足夠強的高能激光輻照空間碎片以形成等離子體反噴，這對遠場激光功率密度和光束質(zhì)量提出了要求。既有足夠高的激光功率密度突破碎片材料闕值，形成等離子體反噴；同時有足夠高的平均功率，能夠獲得誘導(dǎo)變軌所需的速度增量。高能激光遠距離輻照空間碎片，需要考慮大氣傳輸效應(yīng)對激光能量衰減和光束質(zhì)量的影響。在保證產(chǎn)生速度增量的遠場激光功率密度和平均功率前提下，進一步考慮大氣傳輸效應(yīng)的影響，可確定地基激光器功率密度和光束質(zhì)量等要求[3]。

2.1 激光燒蝕推進

2.1.1 最佳沖量耦合系數(shù)

高功率激光輻照碎片表面后，在短時間內(nèi)，光斑區(qū)的溫度急劇升高，使材料熔化和氣化，在激光作用下產(chǎn)生高溫高壓等離子體。等離子體向外膨脹噴射，形成羽流。羽流作用使碎片受到一個方向相反的動量作用而獲得速度增量，從而實現(xiàn)對碎片的驅(qū)動[4]，如圖2所示。

圖2 激光燒蝕驅(qū)動碎片示意圖Fig.2 Schematic diagram of laser ablation driving debris

在激光燒蝕碎片的過程中，激光能量與碎片獲得的速度增量由沖量耦合系數(shù)Cm表征，定義如下:

式中，m為目標碎片的質(zhì)量；Δv為速度增量；EL為輻照到碎片上的單脈沖激光能量；P為碎片表面的燒蝕壓力；I為入射激光的功率密度。沖量耦合系數(shù)表征激光能量轉(zhuǎn)化為碎片沖量的能力大小，單位為N·s/J。隨著激光能量密度的增加，碎片將逐漸熔化、氣化，在燒蝕氣化產(chǎn)物作用下，碎片獲得的沖量逐漸增加，使得沖量耦合系數(shù)增大。氣化產(chǎn)物在向外膨脹的同時將吸收入射激光能量發(fā)生電離而產(chǎn)生等離子體。然而當激光能量密度超過一定值時，產(chǎn)生的等離子體將屏蔽部分入射激光，導(dǎo)致碎片表面的吸收降低，使沖量耦合系數(shù)發(fā)生下降，效率降低。Phipps等人[5]試驗給出沖量耦合系數(shù)隨激光能量密度的關(guān)系如圖3所示。

2.1.2 最佳遠場功率密度

圖3 沖量耦合系數(shù)與激光能量密度的關(guān)系Fig.3 Relationship between impulse coupling coefficient and laser energy density

2006年，美國和歐洲16個不同研究小組開展了46組不同激光波長和靶材材料的真空燒蝕實驗，給出了達到最佳沖量耦合時，激光脈寬對應(yīng)的激光能量通量。Phipps等人[6]經(jīng)過分析發(fā)現(xiàn)，激光與碎片相互作用過程中普遍遵守耦合規(guī)律，即達到最佳沖量耦合時，激光與碎片相互作用的耦合時間隨著激光能量通量的增大而延長，在對數(shù)坐標中呈線性關(guān)系，如圖4所示。

圖4 激光能量通量與脈寬的關(guān)系 (最佳沖量耦合)Fig.4 Relationship between laser energy flux and pulse width(optimal impulse coupling)

圖4中線條為實驗數(shù)據(jù)的最佳擬合曲線，即最佳沖量耦合系數(shù)對應(yīng)的激光能量通量與脈寬的關(guān)系?？梢钥闯?，當激光能量通量小于等離子體的點火閾值時，激光與碎片作用主要取決于材料的性質(zhì)，即定態(tài)氣化過程；而當激光能量通量大于等離子體的點火閾值時，激光與碎片的作用演變?yōu)榧す馀c等離子體的相互作用，需要采用等離子體狀態(tài)定標率來確定耦合特性，與材料性質(zhì)關(guān)系不大。

2.2 激光大氣傳輸

高能激光經(jīng)過大氣傳輸時，由于大氣的非線性效應(yīng)，會出現(xiàn)能量衰減、光束彎曲等現(xiàn)象[4]。大氣的非線性效應(yīng)主要分為三類:第一類大氣衰減，包括氣體分子和氣溶膠的吸收和散射；第二類大氣端流，包括大氣閃爍、光束漂移和擴展等；第三類非線性效應(yīng)限制，包括光束熱暈、大氣擊穿和受激拉曼散射等。Phipps等人[2]通過對上述非線性因素進行試驗計算分析，給出了激光波長分別為1.06m和0.53m時，各種非線性因素的近場激光功率密度的限制曲線。

圖5 非線性因素對近場功率密度和脈寬的限制Fig.5 Restriction of nonlinear factors on near-field power density and pulse width

從圖5可以看出，對于波長1.06μm的激光，當脈寬在100ps～10ns之間時，受激拉曼散射是各種非線性因素中的主要限制條件，近場激光功率密度不超過3×106W/cm2時，各種非線性因素的影響可忽略不計。并且脈寬在100ps～10ns的條件下，脈寬越小對近場激光功率密度限制越有所放寬，近場激光功率密度可放寬到5×107W/cm2。對于波長0.53μm的激光有類似結(jié)論。

2.3 激光器指標

地基發(fā)射的高能激光，需克服大氣傳輸效應(yīng)的影響，以足夠高的功率密度輻照空間碎片，這對近場功率密度和光束質(zhì)量提出了要求。目前可以利用的激光器包括脈沖固體激光器、脈沖化學(xué)和氣體激光器以及連續(xù)波氣體激光器。其中脈沖化學(xué)和氣體激光波長較長，達到同樣的遠場效益，需要的發(fā)射鏡口徑為固體激光器的10倍以上。而連續(xù)波氣體激光器由于熱暈造成的傳輸功率密度較低，無法達到最佳沖量耦合功率密度。目前主要可使用的是釹玻璃激光器，激光器指標的選擇包括平均功率、重復(fù)頻率、單脈沖能量、波長、脈寬、透射率、光束直徑等。

2.4 目標捕獲、跟蹤與瞄準

空間碎片瞄準和高能激光發(fā)射，需要具有探測、照明、捕獲、跟蹤、瞄準和發(fā)射功能的自適應(yīng)光學(xué)系統(tǒng)。空間碎片的探測、捕獲、跟蹤，一般采用雷達、被動光學(xué)系統(tǒng)和主動光學(xué)激光雷達系統(tǒng)等。主動光學(xué)激光雷達尚處于技術(shù)論證階段，而美國的Haystack雷達和STARFIRE被動光學(xué)系統(tǒng)在技術(shù)上都有富余，它們的探測精度達到1～2cm，還具有一定的降軌效果評估功能。

3 方案設(shè)計

3.1 系統(tǒng)組成

地基激光空間碎片清除系統(tǒng)一般由5大分系統(tǒng)組成，如圖6所示。

圖6 地基激光空間碎片清除系統(tǒng)結(jié)構(gòu)Fig.6 Composition of space debris removal system using ground-based laser

控制通信系統(tǒng)負責收集和處理其他系統(tǒng)所提供的信息，協(xié)調(diào)各系統(tǒng)的運行，有計劃地組織清除工作。高能激光器是系統(tǒng)的核心，其主要部件是激光發(fā)生器，負責產(chǎn)生高能激光束。捕獲跟蹤瞄準系統(tǒng)負責目標指向和維持。光學(xué)發(fā)射系統(tǒng)負責把激光束發(fā)射到遠場，匯聚到空間碎片上，形成功率密度足夠高的光斑。自適應(yīng)光學(xué)系統(tǒng)負責校正從發(fā)射望遠鏡到空間碎片這一傳輸通道中大氣對激光束產(chǎn)生的畸變。

3.2 工作流程

地基激光空間碎片清除系統(tǒng)的一般工作流程[7]:首先由雷達或光學(xué)探測設(shè)備發(fā)現(xiàn)碎片目標，并將碎片信息數(shù)據(jù)傳送給控制通信系統(tǒng)，控制通信系統(tǒng)經(jīng)過目標確認，引導(dǎo)精密跟蹤瞄準系統(tǒng)捕獲并鎖定目標，精密跟蹤瞄準系統(tǒng)再引導(dǎo)光學(xué)發(fā)射系統(tǒng)對準碎片。當碎片處于適當位置時，控制通信系統(tǒng)發(fā)出清除指令，啟動激光器，激光器發(fā)出光束，對碎片目標進行清除，如圖7所示。

圖7 地基激光空間碎片清除系統(tǒng)工作流程Fig.7 Workflow of space debris removal system using ground-based laser

3.3 參數(shù)設(shè)計

根據(jù)最佳沖量耦合、能量關(guān)系、集合尺寸關(guān)系、速度增量和大氣傳輸?shù)确矫?，通過綜合協(xié)調(diào)Is/Ib，進行推進激光參數(shù)的設(shè)計，如圖8所示。

圖8 要素設(shè)計原則與約束Fig.8 Element design principle and constraints

圖8中I表示功率密度；E表示能量；Cm表示功率耦合密度；τ表示脈寬；A表示面積；d(D)表示直徑；z表示探測距離；m表示質(zhì)量；η表示大氣傳輸效率；α表示光束常數(shù)；S=1/(極限衍射倍數(shù))2表示斯特列爾比，越大光束質(zhì)量越好；下標s表示遠場，b表示近場。I0為最佳沖量耦合時的功率密度，c1、c2為擬合系數(shù)，一般情況下遠場功率工作在Is≥I0下。

3.3.1 波長選擇

根據(jù)大氣窗口，激光波長一般可選為1μm和3～5μm，為了減小遠場光斑尺寸，典型波長可選取為典型值1.06μm。

3.3.2 遠場功率密度

最佳沖量耦合時的遠場激光功率密度:I0=c1/τ1-c2。對于30種金屬材質(zhì)，試驗可得c1≈8.01×105，c2≈0.648；對于15種非金屬材質(zhì)，c1≈5.97×103，c2≈0.408，考慮到碎片材料的概率和分布，總的參數(shù)選取為c1≈2.30×104，c2≈0.446?？梢酝ㄟ^增大脈寬，來減小遠場功率密度。

3.3.3 近場功率密度

近場功率密度主要受非線性因素約束。對于波長1.06μm的激光，當脈寬在100ps～10ns之間時，受激拉曼散射是各種非線性因素中的主要限制條件，非線性因素要求了近場激光功率密度不超過3×106W/cm2?？梢酝ㄟ^減小脈寬，來增加近場功率密度的限制。很顯然，這與遠場功率的需求相互矛盾，需要協(xié)調(diào)考慮。

由于發(fā)射口徑一般固定，而激光波長確定后，可以根據(jù)探測距離z，來調(diào)整遠近場功率以滿足上式要求。

3.3.4 降軌速度增量

激光燒蝕所能提供的速度增量為:

該速度增量不僅和最佳沖量耦合系數(shù)，近場能量和碎片面質(zhì)比有關(guān)，還與大氣傳輸效率η和有效覆蓋率u有關(guān)。根據(jù)碎片的軌道運動特征、姿態(tài)和旋轉(zhuǎn)、質(zhì)量形狀分布特征等，計算出降軌控制所需的速度增量。在口徑、波長和探測距離確定的基礎(chǔ)上，可以通過調(diào)節(jié)發(fā)射功率和重頻，獲得滿足速度增量需求的近場能量。

空間碎片大多帶有自旋，對小尺度的碎片，自旋角速度可以非常大。為研究碎片自旋對激光碎片清除的影響，假設(shè)自旋為剛體轉(zhuǎn)動，激光光斑對碎片進行全覆蓋連續(xù)照射，且燒蝕帶來的碎片材質(zhì)變化可忽略。激光某一瞬間的有效作用面積為:

這里n為激光的方向矢量；nA為碎片平面的瞬時法向；A為碎片平面面積。當n與nA平行時，有效面積最大；當n與nA垂直時，有效面積為0。假設(shè)碎片的自旋軸在慣性空間的指向ns和激光入射方向矢量n短時間內(nèi)不變，則當n與ns平行時，等于ns·nAA，此時有效面積與自旋軸和平面法向的點乘有關(guān)，與自旋角度無關(guān)；當n與ns垂直時，有效面積隨自旋角度呈周期變化，在一個自旋周期內(nèi)，平均有效面積為:

當激光的入射方向n為其他方向時，可將其按這兩個方向分解，并利用按式 (1)計算對應(yīng)的Δv‖與Δv⊥，最終的速度增量為:

對Δv‖，一個周期內(nèi)平均速度增量，指向自旋的方向；對Δv⊥由于正反面的緣故，一個周期內(nèi)的平均速度增量指向激光入射的方向。注意這里的Δv‖與Δv⊥都是標量，實際得到的Δv方向始終指向碎片平面法向nA，呈周期變化。

需要注意的是，當碎片自旋頻率為激光重頻的整數(shù)倍時，碎片每次獲得的速度增量為一固定值。從式 (6)可以看出，當碎片平面法向與激光入射方向垂直時，獲得的瞬時速度增量為0，此時發(fā)生合拍現(xiàn)象，激光沖量的作用效果將大打折扣，在實際應(yīng)用中應(yīng)盡力避免。

4 參數(shù)復(fù)核

1996年，在NASA和美國空軍的聯(lián)合資助下，洛斯阿拉莫斯實驗室的Phipps等人提出利用地基高能脈沖激光清除近地軌道空間碎片的ORION計劃。該計劃設(shè)計配備一臺波長1.06m，發(fā)射口徑6m的釹玻璃激光器，遠場光斑直徑可達到60cm，配有自適應(yīng)光學(xué)系統(tǒng)，在一顆導(dǎo)向星的引導(dǎo)下能夠清除軌道高度1500km范圍內(nèi)的空間碎片[8]。下表給出ORION計劃公布的設(shè)計參數(shù)。

表1 ORION計劃激光器參數(shù)Table 1 Laser parameters of ORION

下面以作用距離z=600km的近期目標設(shè)計參數(shù)為例，對其合理性進行復(fù)核。作用距離z=3000km的遠期目標也有類似結(jié)論。

(1)激光波長1.06μs，處于合適的大氣窗口。

(2)波長1.06μs時，受激熱瑞利散色和熱暈等非線性效應(yīng)限制脈寬τ≤10ns。此時受激拉曼散色和非線性折射限制了近場功率密度Ib＜3～50MW/cm2，且隨著脈寬的減小，功率限制放寬。

(3)激光器近場功率密度:

滿足近場功率密度Ib＜3～50MW/cm2的限制。

(4) 由遠近場功率密度公式Is/Ib=(ηS/α2)，取斯特列爾比S=0.5，大氣傳輸效率η=0.85，高斯光束常數(shù)α=4/π，可以得到:

(5)由τ=0.1ns，取c1=2.30×104，c2=0.446，最佳沖量耦合時的遠場功率密度I0=c1/τ1-c2=7.28×109(W/cm2)， 與上面計算出的遠場功率密度吻合，即該參數(shù)設(shè)計滿足最佳沖量耦合，效率最高。

(6)取碎片為鋁材碎片，密度為ρ=2.7g/cm3，厚度為1cm，沖量耦合系數(shù)為Cm=4×10-5Ns/J，遠場光斑完全覆蓋的情況下，碎片單脈沖獲得的速度增量為:

由于重頻f=50Hz，激光脈沖作用時間為200s時，可以作用104個脈沖，獲得的速度增量累積到100m/s，對于600km的軌道高度，可降軌約400km，達到使碎片墜入大氣層燒毀的目的。

5 結(jié)論

本文首先梳理了地基激光清除空間碎片中涉及的物理事件及關(guān)鍵技術(shù)；然后設(shè)計了清理操作的實施方案，對要素參數(shù)的選取進行了細化；最后針對ORION計劃給出的參數(shù)，進行了逐一復(fù)核，驗證了其合理性。具體包括:波長滿足大氣窗口；近場功率密度滿足非線性效應(yīng)限制；遠場功率密度符合最佳沖量耦合；輻照時間可滿足預(yù)期速度增量等。地基激光除了清除空間碎片外，還有其他多種用途[9]，如對空間目標進行軌道修正以避免碰撞在軌航天器，精確控制大型空間目標再入大氣層墜毀的位置，將地球同步軌道上的空間碎片推入廢棄軌道等，在本文的基礎(chǔ)上可以進一步開展研究。