基于粒子群的非規(guī)則區(qū)域連續(xù)覆蓋星座設(shè)計

劉 明 楊 明 高 興 楊 丁 葛亞杰

1.國防科技大學(xué), 長沙410073 2.空間物理重點實驗室, 北京100076

0 引言

采用對地靜止軌道衛(wèi)星進行網(wǎng)絡(luò)通信傳輸服務(wù)具有發(fā)射成本高、空間資源有限等缺點，近年來，隨著各類微小衛(wèi)星和低成本、快速響應(yīng)運載器技術(shù)的發(fā)展，利用近地軌道衛(wèi)星組網(wǎng)來執(zhí)行網(wǎng)絡(luò)通信傳輸服務(wù)已成為各國的研究熱點。美國Space-x公司提出的“星鏈計劃”以及中國的“虹云工程”即屬于此類低軌星座，且都已成功部署相關(guān)技術(shù)驗證衛(wèi)星。近地軌道衛(wèi)星組網(wǎng)實現(xiàn)區(qū)域覆蓋時，要求單顆衛(wèi)星的星下點能夠?qū)崿F(xiàn)周期性的重復(fù)，從而保證在某天的固定時段或者全天候?qū)μ囟▍^(qū)域進行連續(xù)覆蓋，圓回歸軌道的軌道高度均勻，適合進行對地觀測，此類星座中的單顆衛(wèi)星軌道通常基于圓回歸軌道進行設(shè)計[1]。共地面軌跡星座是一類特殊的星座形式，星座中的所有衛(wèi)星具有一致的星下點運行軌跡，適用于特定區(qū)域的覆蓋系統(tǒng)[2]。近地軌道衛(wèi)星區(qū)域周期性連續(xù)覆蓋星座的設(shè)計方法主要分為解析法和優(yōu)化法。朱仁璋討論了衛(wèi)星循環(huán)軌道問題并給出了一種覆蓋軌道的設(shè)計方法[3]，吳廷勇等針對高軌通信衛(wèi)星對規(guī)則區(qū)域連續(xù)覆蓋問題，采用遺傳算法求解了共星下點軌跡星座的設(shè)計問題[2,4-5]。Gao等考慮傳感器模型，針對特定緯度范圍內(nèi)具有一定時間間隔的覆蓋問題，采用解析法設(shè)計了星座參數(shù)[6]。文獻[7-8]將包含特定目標城市的規(guī)則區(qū)域連續(xù)覆蓋星座設(shè)計問題轉(zhuǎn)化為多目標優(yōu)化問題進行了求解。張文敦基于解析法設(shè)計了針對特定目標點覆蓋任務(wù)的共星下點軌跡星座，并采用遺傳算法對星座組網(wǎng)形式進行了優(yōu)化[9]。馬劍等考慮重訪周期約束，以星座內(nèi)衛(wèi)星數(shù)量最少為目標，優(yōu)化設(shè)計了共星下點軌跡星座實現(xiàn)了對規(guī)則區(qū)域的覆蓋[10]。范麗提出一種搭載不同有效載荷的區(qū)域覆蓋混合星座設(shè)計方法[11]。

針對非規(guī)則區(qū)域連續(xù)覆蓋星座設(shè)計問題，受推導(dǎo)條件限制，以上研究成果的解析法難以直接解析求解星座設(shè)計參數(shù)[12]；而采用數(shù)值尋優(yōu)的方式，較多優(yōu)化變量的星座設(shè)計參數(shù)將會導(dǎo)致計算效率較差[13-15]。針對該問題，本文采用了一種具有回歸軌道屬性的共星下點軌跡星座方案，提出了一種優(yōu)化法與解析法相結(jié)合的星座優(yōu)化設(shè)計方法，基于衛(wèi)星對地觀測的傳感器視場覆蓋模型，采用粒子群算法(PSO)優(yōu)化算法，以非規(guī)則區(qū)域覆蓋范圍最大為性能指標，選擇軌道傾角和升交點赤經(jīng)為優(yōu)化變量，獲得單顆衛(wèi)星的最優(yōu)軌道參數(shù)，在此基礎(chǔ)上，根據(jù)共星下點軌跡星座約束關(guān)系，解析求解星座軌道參數(shù)，可實現(xiàn)星座軌道的快速求解，同時滿足非規(guī)則區(qū)域連續(xù)覆蓋的性能要求。

1 回歸軌道計算模型

針對圓軌道衛(wèi)星，建立地心慣性坐標系E-xyz，原點位于地心，其中x軸在赤道平面內(nèi)指向春分點方向，z軸由地心指向北極，y軸構(gòu)成右手坐標系。如圖1所示的軌道半徑a、軌道傾角i、升交點赤經(jīng)Ω、緯度幅角u(升交點到衛(wèi)星所在位置的地心夾角)4個參數(shù)即可表示衛(wèi)星在空間中的位置。

圖1 圓軌道描述參數(shù)示意圖

地球引力非球形攝動(地球并不是一個正球體，其內(nèi)部的密度分布不均勻而產(chǎn)生的攝動稱為非球形攝動)對近地軌道具有較大的影響，衛(wèi)星升交點赤經(jīng)Ω與緯度幅角u隨在軌時間的變化而變化，其變化趨勢由軌道半徑a和軌道傾角i確定。J2項是非球形攝動的主要攝動項，考慮J2項攝動可滿足近地軌道設(shè)計時的精度要求，相應(yīng)升交點赤經(jīng)與緯度幅角平均攝動為

(1)

(2)

其中，RE為軌道半徑，ω為圓軌道角速度。

軌道周期T定義為衛(wèi)星由升交點出發(fā)繞行一周再次回到升交點的時長，可根據(jù)下式進行計算

(3)

因此，一個軌道周期內(nèi)升交點赤經(jīng)的增量ΔΩ為

(4)

式中，ωe為地球自轉(zhuǎn)角速度。

回歸軌道或準回歸軌道是星下點按照一定周期重復(fù)的軌道，回歸軌道衛(wèi)星的星下點軌跡重復(fù)周期為1天，準回歸軌道的重復(fù)周期為多天。當滿足式(5)時，星下點軌跡在d天后恰好轉(zhuǎn)過n圈回到原位。

n·ΔΩ=2π·d

(5)

當d為1時，式(5)表示回歸軌道條件，當d大于1時，表示準回歸軌道條件。以下討論的所有問題均為回歸軌道問題，因此，式(5)可以改寫為

n·ΔΩ=2π

(6)

即衛(wèi)星在一天時間內(nèi)運行n圈后，星下點軌跡回歸，聯(lián)立式(3)和式(4)得到

(7)

將式(1)與式(2)代入式(7)，得到考慮J2項攝動的回歸軌道條件

(8)

通常而言，軌道高度H小于500km時，氣動力攝動項對軌道的影響較大，將引起衛(wèi)星能量的耗散，軌道高度難以長時間保持，因此近地軌道高度H一般在500km至1000km之間。由于J2項攝動在式(8)中的影響較小，可忽略其影響，估算H∈[500km,1000km]時，所對應(yīng)的回歸周期T內(nèi)的運行圈數(shù)n∈[13.66,15.18]，取整后n=14，15，其相應(yīng)的軌道高度約為800km以及500km。

2 共星下點軌跡星座約束條件

近地回歸軌道衛(wèi)星只能以回歸周期一天為時間間隔對覆蓋區(qū)域進行重訪，無法在較長時間內(nèi)實現(xiàn)持續(xù)性的覆蓋。針對此類任務(wù)，具有回歸軌道屬性的共星下點軌跡星座可實現(xiàn)對特定區(qū)域的周期性連續(xù)覆蓋。此衛(wèi)星組網(wǎng)星座的優(yōu)點是滿足單顆衛(wèi)星的覆蓋效果最優(yōu)時，組網(wǎng)形成的星座覆蓋效果最好。共星下點軌跡星座設(shè)計一系列軌道高度與軌道傾角相同的近地回歸軌道，每個軌道面內(nèi)一顆衛(wèi)星，考慮地球自轉(zhuǎn)特性、回歸軌道參數(shù)、傳感器模型的條件下，通過設(shè)計各軌道面間的升交點經(jīng)度差ΔΩ以及緯度幅角差Δu，實現(xiàn)對特定區(qū)域的周期性連續(xù)覆蓋。

如圖2所示為相鄰軌道面內(nèi)兩顆衛(wèi)星的空間幾何關(guān)系，衛(wèi)星Si的升交點赤經(jīng)為Ωi，衛(wèi)星Sj的升交點赤經(jīng)為Ωj，兩軌道面的升交點赤經(jīng)差值ΔΩ=|Ωi-Ωj|，兩軌道同為順行軌道，Δu是衛(wèi)星Sj滯后于衛(wèi)星Si的緯度幅角差值。

圖2 相鄰軌道面衛(wèi)星空間幾何關(guān)系

由圖2可知，若衛(wèi)星Sj由當前位置運行到其升交點λj的時間和地球自轉(zhuǎn)ΔΩ的時間相同，則衛(wèi)星Si和Sj具有相同的星下點軌跡，由此可得共星下點軌跡星座相鄰軌道面衛(wèi)星應(yīng)滿足的條件為

(9)

考慮地球非球形引力攝動J2項，將式(1)和式(2)代入式(9)得到

(10)

依據(jù)式(10)即可在已知相鄰軌道緯度幅角相位差Δu的條件下，快速解析計算升交點赤經(jīng)差值ΔΩ。

3 非規(guī)則區(qū)域連續(xù)覆蓋星座設(shè)計方法

3.1 基于傳感器特性的星座分布設(shè)計

衛(wèi)星通常攜帶各類傳感器作為有效載荷，實現(xiàn)對地面目標的通信中繼、對地觀測等任務(wù)，其覆蓋范圍通常由軌道高度、傳感器視場角、傳感器視場形狀等因素決定，典型的對地觀測為圓錐形視場傳感器，圓錐形傳感器的視場如圖3所示，衛(wèi)星所在位置為S，星下點為S′，其軌道高度為SS′長度H，視場半錐角(Cone Half Angle)定義為θ，每一時刻傳感器所覆蓋的圓形區(qū)域半徑AS′長度為d，掃過的地面幅寬為2d，視場投影半徑d可根據(jù)集合關(guān)系計算

d=H·tan(θ)

(11)

圖3 圓錐形傳感器示意圖

(12)

圖4 圓錐形傳感器地面覆蓋區(qū)域示意圖

在確定兩星下點軌跡距離2Δs的條件下，可求出相應(yīng)的地心夾角β

β=2Δs/R

(13)

其中，R為地球平均半徑。地心夾角β即為相鄰軌道面緯度幅角差Δu。

一條完整的共地面軌跡衛(wèi)星環(huán)所分布的衛(wèi)星數(shù)量由式(14)計算

m=2π/β

(14)

若m不為整數(shù)，應(yīng)對其向上取整。

3.2 單顆衛(wèi)星軌道參數(shù)優(yōu)化數(shù)學(xué)模型

由于星座內(nèi)衛(wèi)星星下點軌跡相同，只要單顆衛(wèi)星的覆蓋效果最佳，即可保證星座對于特定區(qū)域的覆蓋效果最佳。分析可知，共星下點軌跡星座的設(shè)計參數(shù)主要包含：軌道半徑a；軌道傾角i；衛(wèi)星總數(shù)；星座內(nèi)第1顆衛(wèi)星的升交點赤經(jīng)Ωo；相鄰軌道面間的升交點赤經(jīng)差ΔΩ。

已知傳感器特性，對覆蓋效果產(chǎn)生影響的參數(shù)主要包括a、i和Ωo。定義傳感器1個軌道周期內(nèi)所掃過的面積在總面積中的占比為特定區(qū)域周期累計覆蓋率η。在軌道半徑以及傳感器視場角選定的情況下，累計覆蓋率僅與i及Ωo相關(guān)，因此在合理選取軌道半徑以及傳感器視場角的情況下，只需對i和Ωo進行優(yōu)化設(shè)計。建立如下的優(yōu)化問題數(shù)學(xué)模型，性能指標為特定區(qū)域周期累計覆蓋率η最大。

(15)

其中，[imin,imax]為優(yōu)化i的取值范圍，[Ωmin,Ωmax]為優(yōu)化Ωo的取值范圍。以上優(yōu)化變量的約束根據(jù)任務(wù)區(qū)域的形狀及位置確定。

3.3 基于粒子群算法(PSO)求解優(yōu)化問題

選用粒子群算法(PSO)求解式(15)所示的參數(shù)優(yōu)化問題，PSO算法起源于對鳥群捕食行為的觀察和研究，其基本設(shè)計思想來源于2個方面：人工生命和進化計算。在PSO算法中，用沒有質(zhì)量和體積的粒子來模擬捕食的鳥，每個粒子只有簡單的行為規(guī)則，通過粒子間的集體協(xié)作使群體達到最優(yōu)[15]。

PSO算法中，每個粒子根據(jù)它自身和群體的經(jīng)驗，向更好的位置“飛行”。粒子本身找到的最優(yōu)解稱為個體最優(yōu)位置，就是每個粒子在飛行過程所經(jīng)歷過的最好位置。整個群體目前找到的最優(yōu)解稱為全局最優(yōu)值，就是整個群體所經(jīng)歷過的最好位置。分別用D維速度矢量Vq=(v1,v2,…,vqd,…,vqD)和位置矢量Xq=(x1,x2,…,xqd,…,xqD)表示第q個粒子的狀態(tài)，則每個粒子根據(jù)以下公式來更新自己的速度和位置，從而產(chǎn)生新一代群體。

(16)

PSO算法的主要運算過程如下：

1)初始化。在問題搜索空間允許范圍內(nèi)隨機設(shè)置粒子的初始位置和速度，設(shè)定PSO算法參數(shù)和算法終止條件。

2)適應(yīng)度評價。根據(jù)目標函數(shù)計算每個粒子的合適值，根據(jù)個體極值和全局極值評選條件，計算每個粒子的最優(yōu)和群體的最優(yōu)。

3)更新粒子。根據(jù)粒子速度和位置更新公式計算，并根據(jù)粒子速度和位置的限定范圍，調(diào)整粒子位置和速度。

4)算法結(jié)束判斷。檢查終止條件，如果滿足條件則終止算法，輸出最優(yōu)解，否則轉(zhuǎn)到步驟2，繼續(xù)迭代。

(17)

3.4 非規(guī)則區(qū)域連續(xù)覆蓋星座設(shè)計流程

采用回歸軌道的共地面軌跡星座優(yōu)化設(shè)計流程為：

1)依據(jù)覆蓋區(qū)域的位置形狀，確定軌道傾角i以及首顆衛(wèi)星的升交點赤經(jīng)Ωo取值范圍，基于3.3節(jié)所述的參數(shù)優(yōu)化方法求取一個軌道周期內(nèi)累計覆蓋率最大所對應(yīng)的軌道傾角i以及Ωo；

2)已知軌道傾角i，選取回歸周期內(nèi)運行圈數(shù)n，迭代求解式(8)等于0時所對應(yīng)的軌道半徑a，進而求得軌道高度H；

3)已知軌道高度H，依據(jù)傳感器類型以及覆蓋區(qū)域幅寬約束，設(shè)計合理的傳感器視場角，基于3.1節(jié)方法計算共地面軌跡衛(wèi)星環(huán)內(nèi)所包含的衛(wèi)星總數(shù)m及緯度幅角差值Δu；

4)根據(jù)持續(xù)覆蓋時間長度，計算共地面軌跡衛(wèi)星環(huán)數(shù)量C，從而得到衛(wèi)星總數(shù)N；

5)基于式(10)計算得到各衛(wèi)星升交點赤經(jīng)Ωq(q=1,…,N)以及緯度幅角uq(q=1,…,N)，完成共地面軌跡星座設(shè)計。

圖5 星座設(shè)計流程

4 仿真驗證

仿真中的覆蓋區(qū)域為以Place1(18.45°N,108.97°E)，Place2(36.71°N ,126.46°E)，Place3(40.24°N ,115.46°E)，Place4(22.23°N ,97.36°E)四地位置坐標連線構(gòu)成的非規(guī)則封閉區(qū)域，覆蓋時間為每天的00：00(UTCG)至03：00(UTCG)的不間斷覆蓋，衛(wèi)星傳感器類型為圓錐形傳感器，視場半錐角為50°。

按照3.4節(jié)所述的設(shè)計步驟，對首顆衛(wèi)星的升交點赤經(jīng)Ωo以及軌道傾角i進行優(yōu)化設(shè)計，相應(yīng)的取值范圍為Ωo∈[105,150]，i∈[30,70]。在衛(wèi)星軌道高度為800km，選擇傳感器視場半錐角(Cone Half Angle)θ為50°的條件下，一個周期的累計覆蓋率η隨Ωo及i變化的趨勢如圖6所示。

圖6 單周期累計覆蓋率隨升交點赤經(jīng)及軌道傾角變化趨勢

采用粒子群算法(PSO)優(yōu)化得到首顆衛(wèi)星升交點赤經(jīng)Ωo為127.5°，軌道傾角i為55°。選取回歸周期內(nèi)運行圈數(shù)n為14。采用迭代方法求解等式約束式(8)，獲得對應(yīng)的軌道半徑a為7211.14km，相應(yīng)的軌道高度H為833km。

傳感器視場半錐角(Cone Half Angle)θ為50°的條件下，基于3.1節(jié)方法計算共地面軌跡星座內(nèi)相鄰衛(wèi)星的緯度幅角差值Δu應(yīng)為18°，相鄰衛(wèi)星的升交點赤經(jīng)差值ΔΩ為1.286°。

相鄰衛(wèi)星空間相對位置及覆蓋特性如圖7所示，一條完整的共地面軌跡衛(wèi)星環(huán)所需衛(wèi)星數(shù)量m為20顆，為實現(xiàn)3小時連續(xù)覆蓋，共地面軌跡衛(wèi)星環(huán)數(shù)量C=1.5，從而得到衛(wèi)星總數(shù)N=30，基于STK的仿真表明，所設(shè)計的共地面軌跡星座可以實現(xiàn)所要求的非規(guī)則區(qū)域每日固定時間連續(xù)覆蓋3小時的需求。

圖7 相鄰衛(wèi)星空間位置及覆蓋特性

表1 星座組網(wǎng)參數(shù)

上述對非規(guī)則區(qū)域3小時的連續(xù)覆蓋仿真的問題中，30顆衛(wèi)星可實現(xiàn)100%非規(guī)則區(qū)域的連續(xù)覆蓋要求，對應(yīng)星座軌道參數(shù)的計算時間在5s以內(nèi)，表明該方法求解效率高，星座覆蓋效果好。

5 結(jié)論

針對非規(guī)則區(qū)域連續(xù)覆蓋的問題，采用回歸軌道的共星下點軌跡星座方案，設(shè)計了一種優(yōu)化法與解析法相結(jié)合的星座優(yōu)化設(shè)計方法。該方法基于粒子群算法(PSO)優(yōu)化求解單顆衛(wèi)星軌道參數(shù)，共星下點軌跡星座約束關(guān)系解析求解星座軌道參數(shù)，具有較高的求解效率，為衛(wèi)星組網(wǎng)的工程實現(xiàn)提供了新的設(shè)計思路和新途徑。