衛(wèi)星導(dǎo)航系統(tǒng)與射電天文兼容性分析

車戰(zhàn)武，劉志儉，申杰敏

(北京衛(wèi)星導(dǎo)航中心，北京100094)

0 引言

射電天文業(yè)務(wù)(RAS)［1］本身并不需要主動向外發(fā)射無線電信號，不會對其他無線電業(yè)務(wù)造成有害干擾;同時，由于射電天文業(yè)務(wù)所接收的是來自宇宙的無線電波，這些無線電電波與人為發(fā)射的無線電電波相比，信號強度極其微弱，大約在10～30 Wm2Hz量級，使得射電天文觀測非常容易受到其他無線電業(yè)務(wù)的干擾［2］。因此，世界無線電通信大會在《無線電規(guī)則》中給出了明確的射電天文業(yè)務(wù)頻率劃分［3］，并形成相關(guān)建議書，對射電天文業(yè)務(wù)進行特殊保護。

北斗衛(wèi)星導(dǎo)航區(qū)域系統(tǒng)已開始正式提供服務(wù)，其使用的1 559～1 610 MHz頻段［4］與射電天文業(yè)務(wù)的1 610．6 ～1 618．8 MHz頻段為相鄰頻段，北斗衛(wèi)星導(dǎo)航系統(tǒng)在確保自身系統(tǒng)穩(wěn)定運行的條件下，必須遵守《無線電規(guī)則》對射電天文業(yè)務(wù)的特殊保護要求。

國際電聯(lián)739決議僅對RAS鄰頻內(nèi)的衛(wèi)星導(dǎo)航業(yè)務(wù)(RNSS)帶外發(fā)射的門限提出限值要求，并在1583建議書(ITU-R M．1583)中給出了非GSO系統(tǒng)在射電天文站址處產(chǎn)生的無用發(fā)射電平的分析方法，基于上述決議和建議書給出的計算方法，對北斗衛(wèi)星導(dǎo)航系統(tǒng)(GSO、MEO、IGSO［5］)為保護 RAS 所需進行的帶外抑制進行了仿真計算。

1 射電天文對RNSS的帶外發(fā)射限制

為保護射電天文業(yè)務(wù)，要求與射電天文業(yè)務(wù)鄰近頻段內(nèi)的RNSS業(yè)務(wù)帶外發(fā)射不能超過下面門限［6］:

對于使用1 559～1 610 MHz頻段的GSO衛(wèi)星而言，其泄露到射電天文頻段 1 610．6～1 613．8 MHz的功率必須小于:

①對于單反射面譜線觀測而言，在基準(zhǔn)帶寬為20 kHz、累積時間為2 000 s的條件下，pfd在98%的概率條件下不得超過－194 dB(W/m2);

② 對于VLBI而言，在基準(zhǔn)帶寬為20 kHz、累積時間為2 000 s的條件下，pfd在98%的概率條件下不得超過－166 dB(W/m2)。

(2)對于使用1 559～1 610 MHz頻段的NGSO衛(wèi)星而言，其泄露到射電天文頻段 1 610．6～1 613.8 MHz的功率必須小于:

①對于單反射面譜線觀測而言，在基準(zhǔn)帶寬為20 kHz、累積時間為2 000 s的條件下，epfd不得超過－258 dB(W/m2);

② 對于VLBI而言，在基準(zhǔn)帶寬為20 kHz、累積時間為 2 000s的條件下，epfd不得超過－230 dB(W/m2)。

此外，不論何時收到俄羅斯 GLONASS/GLONASS-M衛(wèi)星無線電導(dǎo)航系統(tǒng)的相關(guān)協(xié)調(diào)或通知資料，該決議均不適用于1 559～1 610 MHz頻段的GLONASS/GLONASS-M衛(wèi)星無線電導(dǎo)航系統(tǒng)目前和未來的指配。由于在1 610．6～1 613．8 MHz頻段內(nèi)GLONASS/GLONASS-M衛(wèi)星無線電導(dǎo)航系統(tǒng)已經(jīng)給予了射電天文業(yè)務(wù)足夠的保護，且這種保護將繼續(xù)沿用到俄羅斯政府以及GLONASS/GLONASS-M系統(tǒng)的通知主管部門與IUCAF(射電天文學(xué)和空間科學(xué)頻率分配聯(lián)合委員會)之間、以及之后與其他主管部門之間的雙邊協(xié)議中。

2 分析方法及過程

ITU-R M．1583建議書給出了評估由某個非GSO系統(tǒng)在射電天文站址處產(chǎn)生的無用發(fā)射電平的方法［7］。該方法的基礎(chǔ)是把天空分成尺寸幾乎相同的小區(qū)并進行統(tǒng)計分析，在統(tǒng)計分析中RAS天線的指向和衛(wèi)星星群的啟動時間都是隨機變量。對于每次評估，無用發(fā)射電平(以 epfd表示)都在2 000 s的時段內(nèi)平均。

關(guān)于對GPS、GLONASS和GALILEO的仿真，所選的特性相當(dāng)于德國Effelsberg射電望遠鏡，采用100 m直徑和約64 dBi的最大增益。天線方向性圖和最大天線增益取自ITU-R RA．1631建議書［8］。

2．1 分析方法(ITU-R M．1583)

考慮衛(wèi)星系統(tǒng)和射電天文望遠鏡天線特性參數(shù)后，采用等效功率通量密度(epfd)的概念來計算NGSO衛(wèi)星系統(tǒng)無用發(fā)射對射電天文望遠鏡的干擾。epfd是所有衛(wèi)星發(fā)射等效源(主瓣峰值功率)的集合對射電天文望遠鏡的貢獻。

2．1．1 輸入條件

在計算對射電天文望遠鏡干擾時，其累積統(tǒng)計時間要在2 000 s。2種業(yè)務(wù)需要考慮的輸入條件有:

①NGSO衛(wèi)星系統(tǒng):射電天文望遠鏡觀測站的可見衛(wèi)星數(shù);衛(wèi)星的軌道參數(shù);每顆衛(wèi)星的pfd值。

②射電天文望遠鏡:接收天線站址;接收天線旁瓣和天線增益;天線實際指向范圍;接收天線視軸指向;射電天文望遠鏡天線視軸與發(fā)射衛(wèi)星的夾角;累積時間(2 000 s)。

2．1．2 射電天文望遠鏡接收端的epfd值計算

射電天文望遠鏡指向NGSO衛(wèi)星方向的天線接收增益隨著衛(wèi)星的運動發(fā)生變化，但在分析干擾時只考慮望遠鏡指向衛(wèi)星方向，且接收增益大于0 dBi時的影響。計算epfd的數(shù)學(xué)公式和模型來自無線電規(guī)則的第22.5C款［9］。

式中，Na為接收端可觀測到的衛(wèi)星總數(shù);i為用于計算的發(fā)射衛(wèi)星序號;Pi為在參考帶寬內(nèi)，發(fā)射衛(wèi)星天線輸入端的射頻功率譜，dB(W/MHz);θi為發(fā)射衛(wèi)星視軸與接收指向視軸間的離軸角(°);Gt(i)為指向接收機方向的衛(wèi)星發(fā)射天線離軸增益(dB);di為發(fā)射衛(wèi)星與接收機間距離(m);φi為接收機指向軸與發(fā)射衛(wèi)星視軸間的離軸角(°);Gr(φi)為指向發(fā)射衛(wèi)星的接收機接收天線旁瓣增益(dB);Gr，max為接收機的最大增益(dB);epfd為某一時刻的接收機等效功率通量密度(dB(W/(m2·MHz)))。

但式(1)的假設(shè)前提條件是干擾產(chǎn)生在接收天線的視軸方向。然而就射電天文這一保護標(biāo)準(zhǔn)來說，則要求考慮到射電天文天線的0 dBi圈。因此式(1)的干擾源在0 dBi接收增益方向可以表示為:

式中，epfdGr=0dBi是2 000 s累積時間內(nèi)的平均計算結(jié)果。

2．2 分析過程

首先假定天線將位于某一個特定區(qū)域，針對射電天文業(yè)務(wù)天線的指向進行隨機選擇，并隨機選擇觀測星座的起始時間，根據(jù)2 000 s的積分時間對每一時間抽樣的epfd進行評估。然后再針對所選擇的指向和觀測星座的起始時間，進行與此試驗相對應(yīng)的平均epfd的計算。

重復(fù)進行此操作，以便獲得所考慮區(qū)域的epfd分布統(tǒng)計。此方法涉及到多次試驗，其中每個試驗均根據(jù)2 000 s的積分間隔計算平均epfd值。試驗的次數(shù)越多，此分布情況越精確。需要進行足夠次數(shù)的試驗才能獲得所需的可信值。尤其是，試驗數(shù)目與2 000 s的積分時間的乘積應(yīng)顯著高于星座的階段。而且，還有必要確保星座整個階段有適當(dāng)?shù)慕y(tǒng)計抽樣。一當(dāng)發(fā)現(xiàn)在分布中沒有進一步的顯著變化，即可得出結(jié)論，試驗的次數(shù)已經(jīng)足夠。此檢查可作為模擬工作的一項不可分割的內(nèi)容自動進行，或通過定期停止模擬工作人工進行。

為了全面的彌補隨機選擇的可能帶來的分析不足，對任意一個天線位置，采用遍歷算法，對天線指向以空間5°間隔建立格網(wǎng)，在一個衛(wèi)星系統(tǒng)的回歸周期以100 s為間隔進行積分周期2 000 s的試驗并統(tǒng)計。

3 仿真計算

3．1 假設(shè)條件

所有系統(tǒng)均為假設(shè)特性，不代表最終北斗全球系統(tǒng)特性。北斗全球系統(tǒng)B1信號發(fā)射特性如表1所示。

表1 B1信號發(fā)射特性(未濾波)

3．2 北斗全球系統(tǒng)MEO軌道特性

BD全球MEO星座由27顆衛(wèi)星組成，軌道高度27 878 km，3 個軌道面，間隔120°，軌道傾角 55°，同一軌道面衛(wèi)星間隔為45°，偏心率為0。軌道回歸周期為8天。

3．3 北斗全球系統(tǒng)IGSO軌道特性

BD全球IGSO星座由3顆衛(wèi)星組成，軌道高度42 164.2 km，3 個軌道面，間隔120°，軌道傾角 55°，偏心率為0，與地球赤道交點為118°E。軌道回歸周期為1天。

3．4 分析計算

3．4．1 GEO 衛(wèi)星

GEO衛(wèi)星遵循表1規(guī)則，直接采用鏈路計算方法計算，可知在假設(shè)條件下，GEO衛(wèi)星在帶外頻譜自然降低的條件下抑制16 dB后可滿足RAS的保護需求。

3．4．2 MEO 衛(wèi)星

采用第3節(jié)計算方法，對于某一位置，需進行9 455 616次2 000 s的積分計算，計算量巨大，極為耗時，因此僅僅選取高、中和低緯3個位置進行了計算，結(jié)果如表2所示。

表2 MEO衛(wèi)星計算結(jié)果

3．4．3 IGSO 衛(wèi)星

采用第3節(jié)計算方法，對于某一位置，需進行9 455 616次2 000 s的積分計算，計算量巨大，極為耗時。IGSO衛(wèi)星僅僅覆蓋有限區(qū)域，因此僅對覆蓋區(qū)域內(nèi)中緯度的一個位置進行了計算，結(jié)果如表3所示。

表3 IGSO衛(wèi)星計算結(jié)果

3．5 結(jié)論

根據(jù)計算結(jié)果可見，在假定的信號調(diào)制方式和星座構(gòu)型條件下，MEO在帶外頻譜自然降低的條件下再抑制39 dB即可滿足RAS保護要求;IGSO衛(wèi)星在帶外頻譜自然降低的條件下再抑制28 dB即可滿足RAS保護要求，GEO衛(wèi)星在帶外頻譜自然降低的條件下抑制16 dB即可滿足RAS保護要求。

綜合考慮以下4方面因素，建議將MEO帶外需求設(shè)定為40 dB，IGSO帶外需求設(shè)定為30 dB，GEO帶外需求設(shè)定為16 dB:

①由于計算時間的關(guān)系，沒有對RAS天線位置進行遍歷仿真，有可能存在比39 dB更為嚴格的抑制需求;

②多GNSS系統(tǒng)的存在，導(dǎo)致多系統(tǒng)同時使用時RAS可能受到干擾;

③在計算中使用了最大EIRP，未考慮發(fā)射天線離軸角，使計算條件過于嚴苛，根據(jù)RAD文件，北斗發(fā)射天線增益最大最小差異為2 dB;

④在MEO和IGSO共同覆蓋區(qū)域，帶外抑制要求將大于以上計算結(jié)果。

4 結(jié)束語

隨著全球衛(wèi)星導(dǎo)航系統(tǒng)的廣泛使用，人們在活動的各個領(lǐng)域?qū)πl(wèi)星定位導(dǎo)航服務(wù)的依賴越來越強。國際電聯(lián)739決議要求RAS業(yè)務(wù)鄰頻的RNSS業(yè)務(wù)對其進行保護，這不僅是對北斗導(dǎo)航信號發(fā)射功率的約束，也是實現(xiàn)北斗衛(wèi)系統(tǒng)服務(wù)可用性的重要因素，采用1583建議書的計算方法進行計算，并綜合考慮相關(guān)因素，正是為北斗導(dǎo)航信號帶外抑制設(shè)計提供了合理、可行的仿真數(shù)據(jù)，保證北斗全球系統(tǒng)合法有效的使用1 559～1 610 MHz導(dǎo)航頻段，為用戶提供可靠、穩(wěn)定的服務(wù)。

［1］ International Telecommunication Union．Radio Regulations Articles［R］．Geneva，2012．

［2］ 劉 斌，方 正．射電天文業(yè)務(wù)的無線電管理淺析［J］．中國無線電，2009(3):51 －56．

［3］ Intemationd Telecommunication union．Radio Regulations Articles［R］．Geneva，2012．

［4］ 譚述森．衛(wèi)星導(dǎo)航定位工程(第2版)［M］．北京:國防工業(yè)出版社，2010．

［5］ 朱立東，吳廷勇．衛(wèi)星通信導(dǎo)論［M］．北京:電子工業(yè)出版社，2010．

［6］ International Telecommunication Union．Radio Regulations Resolutions and Recommendations［R］．Geneva，2012:367－372．

［7］ International Telecommunication Union．Radio Regulations ITU-R Recommendations［R］．Geneva，2012:401 －407．

［8］ International Telecommunication Union．Radio Regulations ITU-R Recommendations［R］．Geneva，2012:433 －436．

［9］ International Telecommunication Union．Radio Regulations Articles［R］．Geneva，2012:272．