鄭鵬飛, 陳宏宇, 郭崇濱
(1.中國科學院微小衛星創新研究院, 上海 201203; 2.中國科學院大學, 北京 100049;3.上海微小衛星工程中心, 上海 201203)
目前,低軌巨型星座發展迅速,如中國航天科技集團“鴻雁”星座(360余顆)、中國航天科工集團“虹云”星座(156顆)、美國SpaceX公司“Starlink”星座(4.2萬顆)等。該類星座一般采用大量的低軌微小衛星組網實現近實時覆蓋,經典的Walker-δ、共地面軌跡等星座衛星分布均勻性較差且難以利用解析法快速設計。
應用于低重訪時間覆蓋的經典星座類型主要有2種:①典型的Walker星座[1],包括星形構型、δ構型、σ構型、玫瑰構型等;Peng等[2]提出了Brush Line構型星座,采用多星沿經度掃描的思想實現全球覆蓋,但其連續覆蓋性較差;He等[3]針對通信定位問題,綜合考慮軌道高度、傾角等參數設計了多傾角組合的低軌δ星座實現連續覆蓋,但復合星座構型不穩定且維持較復雜,難以建立星間鏈路。②共地面軌跡(repeat ground track,RGT)星座,該星座所有衛星的地面軌跡重復,但每個軌道面上僅有一顆衛星;Lee[4]針對特定目標的連續覆蓋對比分析了RGT星座和δ星座的數量差異,在重訪時間要求較大(20 min以上)時,RGT星座規模顯著小于δ星座,在重訪時間要求較小(10 min以內)時二者差別不大;SpaceX公司的專利[5]顯示Starlink的Ⅲ期星座(7 518顆)欲采用多個共地面軌跡的星簇組成。
針對以上星座的常規設計方法主要有2種:①利用解析法設計星座參考碼并進行數值仿真、人工選優,一般能夠快速得到滿足任務需求的星座;王茂才等[6]根據任務需求人工選取幾組參數,仿真選優設計了雙層對地觀測星座;Ulybyshev利用解析法分析覆蓋性[7],設計了近極軌和赤道圓軌道弱對稱的非均勻衛星星座。②利用優化算法搜索軌道的最優解或星座參考碼的最優解;馬劍、Lake、賀波勇、沈欣等[8-11]用遺傳算法、多目標進化算法、改進粒子群優化算法等設計了不同任務需求的星座,該方法存在算法求解時間較長、算法復雜、優化目標單一等問題。
Walker[1]提出了覆蓋性和經濟性均較好的Walker-σ星座,該星座的星下點軌跡重合且形成一條類似正弦曲線的不自相交封閉曲線,覆蓋特性均勻,但該構型依賴于回歸軌道,維持該軌道所消耗的燃料較大,不利于低軌星座的長期維護。因此,參考σ星座的設計思想,本文提出了一種低軌鏈形星座,并基于覆蓋要求和傳感器視場角推導了解析設計方案,以對地偵察任務為例設計了一個包含132顆衛星的鏈形星座,從覆蓋性和星間鏈路角度與常規δ星座和RGT星座對比分析,驗證了鏈形星座效能。
Walker-σ星座軌道一般選取中、高回歸地球軌道,從而實現所有衛星的星下點軌跡重合,而低軌非回歸軌道Walker星座無法實現。借鑒σ星座封閉曲線覆蓋特性的思想,設計合理的相位參數,使得低軌Walker星座的星間鏈路連線在二維地圖的投影類似于一條或多條自封閉的正弦曲線,衛星在該曲線上等時間間隔分布,該鏈狀構型隨時間動態移動實現對所有經度范圍的覆蓋,稱之為Walker-鏈形星座,如圖1所示。

圖1 鏈形星座示意圖
為描述鏈形星座構型,參考σ星座描述方法[12],引入相位配置參數M(M為正整數)、衛星數量N,星座中所有衛星的軌道傾角均為i,軌道平均高度均為h,偏心率均為e,近地點幅角均為ω,升交點赤經為Ω,相位角為u,則順序相隨的2顆衛星之間的升交點赤經差和相位差為

(1)
通過選取不同的N和M的值,可得到不同覆蓋性能的鏈形星座。
1.2.1 軌道設計
1) 軌道高度
低地球軌道衛星一般處于Van-Allen內帶以下避免高強度輻射,同時過低的軌道受到的大氣阻力較大,因此可根據傳感器探測距離選擇平均高度為300~1 500 km內的軌道。
2) 軌道傾角
Walker星座一般應用于緯度帶覆蓋,傾角的選擇取決于緯度帶的覆蓋要求和傳感器的視場角,在覆蓋要求為南北緯θN,S角度之間時,設傳感器的半視場角為α,對應的地心角為β,如圖2所示。

圖2 單星地面覆蓋性示意圖
由幾何關系可得

(2)
式中,地球平均半徑R?=6 378.137 km。
則軌道傾角值的選取可參考
i=θN,S-β+ζ
(3)
式中,ζ的一般取值為1°~3°,以完全覆蓋緯度帶的邊緣部分。
3) 軌道偏心率和近地點幅角
通常低地球軌道選取圓軌道,此時偏心率e=0,近地點幅角此時無意義,令ω=0°。
1.2.2 升交點赤經和相位設計
在鏈形星座的相位配置中,為實現對緯度帶的低重訪覆蓋,需使鏈路中順序相隨的衛星覆蓋范圍、鏈與鏈之間的衛星覆蓋范圍在空間上鄰接。
令任意順序相隨的2顆衛星鄰接,可使整個鏈形星座順序鄰接,順序相隨的2顆衛星距離遠近決定了重訪時間的大小。在地心第三赤道坐標系(衛星Sat1位于升交點時)中表示出順序相隨的2顆衛星位置:以赤道面為XOY面,原點O為地心,X軸指向衛星Sat1的升交點,Z軸指向北極。
由圖3可以得到衛星Sat1和衛星Sat2的坐標(φ為赤經,λ為赤緯)為

圖3 地心第三赤道坐標系下的衛星示意圖
(4)

表1 圖3中的主要符號及其含義
將順序相隨的衛星鏈路連線可形成一條閉合的曲線,隨著M取值的不同,鏈路的起伏程度不同,對地面的均勻覆蓋性也不同,用赤經差和赤緯差的比率γ來描述鏈路的起伏程度

(5)
由圖4可知,γ越大,則鏈路的瞬時覆蓋間隙越大,重訪時間越高;γ越小,則鏈路起伏越陡峭,瞬時覆蓋均勻性越好,重訪時間越小。

圖4 M取不同值時鏈形星座起伏程度示意圖
為保證順序相隨的衛星覆蓋范圍鄰接,要求順序相隨的衛星赤緯差符合(6)式要求
Δλ=|λ2-λ1|=2β·ε
(6)
ε定義為連續覆蓋控制因子,ε越小,則星座的連續覆蓋性越好,重訪時間越小,但所需的衛星數量N增加。
由(4)式和(6)式可得Δu的取值范圍

(7)
結合(4)式和(5)式可得M值的選取依據為

(8)
(8)式為超越方程,無法直接求解,可繪制γ關于M和Δu的圖像,根據圖像選取γ取最小值時對應的M值,然后根據Δu的范圍結合(1)式得到N值的范圍,從而得到一組或多組鏈形星座參考碼。
1.2.3 單鏈形星座與多鏈形星座的轉換
鏈形星座是一種特殊的δ星座,可用δ星座參考碼表示。包含單條鏈路的鏈形星座瞬時覆蓋范圍可能無法覆蓋所有經度,可將多個鏈形星座組合實現緯度帶的瞬時全覆蓋,此時仍然為一個特殊的δ星座,因此利用δ星座參考碼進行轉換。
1) 單鏈形星座參考碼轉換為δ星座參考碼
δ星座參考碼為N/P/F,N為衛星數量,P為軌道面數量,F為相位因子,單鏈形星座與δ星座的參考碼對應為[12]

(9)
(9)式中相位因子F是[0,P-1]內的整數,由此整數系數K可唯一確定。
2) 多鏈形星座參考碼轉換為δ星座參考碼
單鏈形星座對應的δ星座參考碼為N1/P1/F1,則包含n條鏈(n為奇數)的多鏈形星座對應的δ星座參考碼N2/P2/F2為
(10)
式中,當n為偶數時,多個單鏈形星座在空間上重疊,無法增加瞬時覆蓋率,因此n僅能取奇數,此時多條鏈形在空間上錯開,瞬時覆蓋率增加。
相比于單鏈形星座,多鏈形星座每個軌道面上的衛星數量擴大n倍,但軌道面數量不變,有利于一箭多星技術的利用,極大地降低發射成本。
1.2.4 設計流程
綜上所述可得鏈形星座設計流程為:
step1 確定傳感器視場角和覆蓋緯度需求,選取合適的軌道;
step2 根據重訪時間的需求選取ε取值范圍,利用圖解法得到一組或多組鏈形星座參考碼;
step3 根據瞬時覆蓋率選取是否轉換為多鏈形星座,并從多組結果中仿真選優,一般選取衛星數量最少的方案作為設計結果。
空間對地偵察星座具有覆蓋范圍廣、探測距離遠、隱蔽性強等優點,已成為重要的情報來源。第16任北約盟軍最高司令斯塔夫里迪斯上將[13]闡述了未來幾十年美軍在太平洋對華作戰的設想,美軍的總體海上戰略態勢是建立全球海上聯盟,并有針對性地增強美方在西太平洋的前沿存在和作戰能力;彭耿等[14]針對衛星偵察時效需求建立了偵察衛星支援反艦導彈作戰的基本模型,計算結果表明當衛星偵察定位精度為4 km時,目指信息要求的最大時延約為4~6 min;郭福成等[15]在分析偵察衛星瞬時覆蓋范圍時提到對于高度在500~1 000 km的衛星,其偵察范圍可達2 500~3 000 km;高世杰等[16]提到對于大規模星座,需具備星間鏈路以降低地面測控壓力并提高數據實時回傳能力。
綜上所述,可得到偵察任務的需求為:
1) 面向全球熱點地區和海洋地區,即在南緯60°和北緯60°之間實現低重訪時間覆蓋;
2) 對南北緯60°之間的任意地點平均重訪時間在4 min以內、最大重訪時間在6 min以內;
3) 星載傳感器的半視場角為50°,可簡化為對地圓錐視場,探測距離小于1 500 km;
4) 星座需在當前星間激光通信載荷能力下建立永久性星間鏈路;
根據任務需求,設計低軌鏈形星座,并按照經典方法設計常規δ星座和RGT星座,對比分析效能。
2.1.1 軌道設計
傳感器視場角固定時,探測距離小于1 500 km范圍內,軌道高度越高,瞬時覆蓋范圍越大,設計軌道高度為800 km,偏心率為0。
結合(2)式和(3)式可得軌道傾角為
(11)
將數值代入(11)式可得軌道傾角為50.44°,令i=53°以完全覆蓋緯度帶的邊緣部分。
2.1.2 升交點赤經和相位設計
由(8)式可知γ的取值主要取決于M,因此連續覆蓋控制因子ε可選取較大的范圍:[0.8,1.2],進而由(7)式可得Δu的范圍為[19.62°,29.73°]。在Δu的取值范圍內,根據(8)式繪制γ關于M和Δu的圖像如圖5所示。

圖5 鏈路起伏程度γ在Δu范圍內關于M的圖像
由圖5可知,γ取最小時對應M值為2,則根據(1)式可得N的取值范圍為(37,55)之間的偶數,可得如圖6所示的單鏈形星座。

圖6 單鏈形星座設計結果M=2
從圖6可以看出,單鏈無法滿足緯度帶的瞬時覆蓋,需設計多鏈形星座,根據(9)式和(10)式將單鏈形星座轉換為多鏈形星座,得到初步設計結果。
依據表2由衛星數量最少的方案開始在STK(美國Analytical Graphics公司衛星工具包)中仿真計算重訪時間,直至符合任務需求,可得M=2,N=44對應δ星座參考碼為132/22/13的多鏈形星座符合要求,根據(1)式可得星座的升交點赤經差和相位差為

表2 鏈形星座初步設計結果

(12)
星座在二維地圖上的投影如圖7所示:

圖7 多鏈形星座設計結果(132顆)
2.2.1 Walker-δ星座常規設計
δ星座的參考碼為N/P/F,P為軌道面數量,F為相位因子,確定δ星座升交點赤經差和同軌道面相位差的方程為[12]

(13)
設計δ星座軌道與鏈形星座相同,采用遺傳算法設計得到的δ星座結果與鏈形星座完全相同,但算法復雜、優化計算時間超過2天;采用經驗法設計時[12],相位因子通常取0或1,通過枚舉法仿真不同衛星數量的常規δ星座,結果如圖8所示。

圖8 平均/最大重訪時間隨衛星數量變化圖
由圖8綜合權衡可得參考碼為170/17/0的常規δ星座基本滿足任務需求且數量最優,如圖9所示。

圖9 常規δ星座設計結果(170顆)
2.2.2 共地面軌跡星座設計
共地面軌跡星座中所有衛星具有相同的地面軌跡,順序相隨的衛星升交點赤經差和相位差需滿足[17]

(14)
式中:ω?=7.292 115×10-5rad/s為地球平均自轉角速度;ωs為衛星的在軌角速度,由軌道高度h決定。
根據(14)式,按照衛星數量由少至多順序,即相位差Δu值由大至小枚舉仿真,可得滿足任務要求的共地面軌跡星座衛星數量N=160,每個軌道面僅有1顆衛星,相位差Δu=32.11°,如圖10所示。

圖10 RGT星座設計結果(160顆)
星座效能主要體現為經濟性、覆蓋性和星間鏈路特性,在均滿足對地偵察任務需求的條件下,對比分析3種星座效能。
由表3可知,鏈形星座所需的衛星數量最少,僅為132顆,而常規δ星座數量是其1.3倍,RGT星座數量是其1.2倍;鏈形星座和常規δ星座每個軌道面均包含多顆衛星,可利用一箭多星技術降低發射成本,而RGT星座每個軌道面僅有1顆衛星,星座部署成本大大增加。

表3 對地偵察星座設計結果
針對緯度帶的低重訪覆蓋,主要的覆蓋性指標有重訪時間和覆蓋時間比,設置仿真時長為一天,通過STK仿真得到3種星座的覆蓋時間比、重訪時間隨緯度變化如圖11~12所示。
由圖11可知,鏈形星座覆蓋時間比在73%以上,而常規δ星座為63%以上,RGT星座僅為56%以上。

圖11 3種星座覆蓋時間比對比圖
由圖12可知,除RGT星座的重復地面軌跡特性使其在極少數地區重訪時間過大,3種星座均滿足最大重訪時間小于6 min、平均重訪時間小于4 min,鏈形星座的平均重訪時間在2 min以內,且在不同緯度上較均勻,而常規δ星座和RGT星座在赤道地區的覆蓋性起伏較大、劣于中高緯度地區。

圖12 3種星座平均和最大重訪時間對比圖
重訪時間越小、覆蓋時間比越大則對地偵察時間越長,因此鏈形星座的覆蓋性能優于常規δ星座和RGT星座。
衛星可以與前后左右4顆衛星建立星間鏈路,攜帶多個激光通信載荷分別對準不同的衛星,因此,通過衛星姿態機動僅能實現2顆衛星的對準,無法同時完成所有激光鏈路的持續對準,這就要求激光通信載荷自身需具備一定的轉動能力[16]。通信距離、轉向能力是限制建立星間激光鏈路的主要指標,如德國Mynaric公司[18]研制的星載激光通信終端通信距離小于4 500 km,方向角轉動范圍為±175°,俯仰角轉動范圍為+5°~-25°。激光通信終端的俯仰改變靠自身轉動實現,大范圍的方向角改變需依賴衛星的偏航姿態機動,而星載傳感器為對地指向,偏航機動對星座覆蓋性能無影響。
依據可見性約束設計永久性鏈路。
1) 鏈形星座星間鏈路設計
鏈形星座的分布特性使其可以按照順序相隨的衛星建立1條閉合的星間鏈路;則3條鏈可分別建立3條閉合的星間鏈路。另外可在鏈與鏈之間建立鏈路,共可建4條閉合鏈路;則鏈形星座的全連通星間鏈路如圖15所示。

圖13 鏈形星座同鏈星間鏈路(1條) 圖14 鏈形星座異鏈星間鏈路(1條) 圖15 鏈形星座星間鏈路設計結果
2) 常規δ星座星間鏈路設計
針對170/17/0的常規δ星座,其不同軌道面對應序號的衛星始終處在同一赤緯上,因此不同軌道面上同序號的衛星共計可組成10條異軌道面星間鏈路;該星座同軌道面上的衛星共計可組成17條閉合的同軌道面星間鏈路;則常規δ星座的全連通星間鏈路如圖18所示。

圖16 常規δ星座異軌道面星間鏈路(10條) 圖17 常規δ星座同軌道面星間鏈路(17條) 圖18 常規δ星座星間鏈路設計結果
3) RGT星座星間鏈路設計
RGT星座具有重復地面軌跡的特性,因此可將RGT星座的所有衛星按軌道面順序相連,得到1條完整的星間鏈路;同時,RGT星座可按照間隔地面軌跡的順序形成1條閉合的星間鏈路;則RGT星座的全連通星間鏈路如圖21所示。

圖19 RGT星座同軌跡鏈路(1條) 圖20 RGT星座異軌跡鏈路(1條) 圖21 RGT星座星間鏈路設計結果
4) 對比分析
對以上3種星座建立的6種鏈路進行分析,結果如表4所示。由表可知,3種星座均能建立全連通的永久性星間鏈路,鏈路指標均在正常范圍內,鏈形星座的同鏈星間鏈路觀測方位角、俯仰角、角速率的變化范圍均較小,能夠較容易地建立3條完整的鏈路實現全部連接,而常規δ星座至少需要10條異軌道鏈路,RGT星座的鏈路距離較大,且僅為1條閉合鏈路、穩健性不高;鏈形星座的異鏈星間鏈路的距離較大,可作為同鏈鏈路的補充,實現星座全連通。

表4 3種星座星間鏈路對比
本文面向低重訪時間覆蓋需求提供了一種低軌巨型鏈形星座及解析設計方法,該方法依據覆蓋要求和傳感器視場角,在超越方程的多組解中仿真選優,能快速設計出符合重訪時間需求的低軌巨型星座。以對地偵察任務為例設計了一個132顆衛星的低軌鏈形星座,相較于常規δ星座、RGT星座所需衛星數量更少、覆蓋更均勻、重訪時間更小、容易建立永久性星間鏈路。未來可針對重訪時間要求θN,S與鏈形星座相位差Δu、相位配置參數M的關系作進一步研究,以縮小連續覆蓋控制因子ε的取值范圍,從而減小超越方程解的數量。