利用SLR對GLONASS—M衛(wèi)星系統(tǒng)精密定軌研究

樊東昊　姜家慶

【摘 要】利用人衛(wèi)激光測距（SLR）技術(shù)，構(gòu)建了GLONASS-M衛(wèi)星的精密定軌模型，并以CODE GNSS數(shù)據(jù)分析中心發(fā)布的GLONASS精密星歷作為標準軌道，得到SLR定軌誤差在測距方向精度可以達到 10cm，對于連續(xù)的21天時間序列SLR定軌精度徑向可以達到30cm，三維方向可以達到1m以內(nèi)，以CODE的軌道產(chǎn)品為標準軌道為SLR定軌結(jié)果進行精度評定，其定軌精度相當。對比分析了不同弧段、不同軌道面的GLONASS-M衛(wèi)星的定軌精度，采用7天的定軌弧長對GLONASS導(dǎo)航衛(wèi)星進行SLR定軌可得到最優(yōu)解。

【關(guān)鍵詞】GLONASS-M；人衛(wèi)激光測距（SLR）；精密定軌

【Abstract】The precision orbit determination model of GLONASS-M satellite is established by using the SLR technology.The GLONASS precision ephemeris published by the CODE GNSS data analysis center is used as the standard orbit to obtain the SLR orbit determination error.The precision of direction can be up to 10cm. For the continuous 21-day time series，the accuracy of SLR orbit determination can reach 30cm in radial direction and less than 1m in three-dimensional direction.The rail product of CODE is the standard orbit for accuracy evaluation of SLR orbit determination.Accuracy.The orbit determination accuracy of the GLONASS-M satellites with different arcs and orbits is analyzed and compared.The optimal solution of the GLONASS navigation satellites is obtained by using the 7-day orbit length.

【Key words】GLONASS-M；Satellite laser ranging（SLR）；Precise orbit

0 引言

人衛(wèi)激光測距（SLR）由于測距精度高等特點，廣泛應(yīng)用于衛(wèi)星精密軌道確定[1-2]、地球自轉(zhuǎn)參數(shù)解算[3]、地球重力場低階球諧系數(shù)的確定[4]、測定臺站坐標[5]、參考框架的建立[6]及軌道檢驗[7-8]等衛(wèi)星大地測量領(lǐng)域。雖然SLR觀測覆蓋區(qū)域有限，且受天氣影響嚴重，但隨著觀測密度的加大，測站數(shù)的增加及測站位置幾何分布合理化，它將繼續(xù)在地學科學研究中發(fā)揮不可替代的重要作用。

由于GNSS全球衛(wèi)星導(dǎo)航系統(tǒng)的微波信號定軌存在一個嚴重的問題，就是軌道的徑向誤差和星載的原子鐘中差高度相關(guān)性，不易分離，相位中心誤差也很難模擬。SLR觀測沒有這種相關(guān)性，搭載激光角反射器的GNSS衛(wèi)星不僅僅能檢核GNSS微波定軌結(jié)果，而且能作為一種完全獨立的精密定軌手段。

目前GPS衛(wèi)星只有兩顆衛(wèi)星裝有激光反射器，GPS-35/36，數(shù)據(jù)量少，對其進行SLR精密定軌很難得到較為理想的定軌結(jié)果.對比而言GLONASS導(dǎo)航衛(wèi)星均裝有激光反射器，并且目前GLONASS上裝載的激光反射器的反射信號的反射率較高，所以相對來說GLONASS的激光測距數(shù)據(jù)相對較多，雖然GLONASS衛(wèi)星屬于中高軌衛(wèi)星，觀測數(shù)據(jù)較低軌衛(wèi)星少，也可以得到精度較為理想的定軌結(jié)果。不僅可以利用SLR進行單獨定軌還可以與微波數(shù)據(jù)進行聯(lián)合定軌。對同一顆導(dǎo)航衛(wèi)星用2種不同的觀測手段，可以獲得2種不同類型的觀測數(shù)據(jù)。利用這2種觀測數(shù)據(jù)，可以分別解算同一顆衛(wèi)星在2個不同系統(tǒng)中的軌道，用此來直接地比較這2個系統(tǒng)間的差異，或用一個系統(tǒng)的解算值去評估另一系統(tǒng)解算值的精度。還可對這2種數(shù)據(jù)聯(lián)合求解，去解算一個公共軌道，估算一些公共參數(shù)[9]。所以針對GLONASS導(dǎo)航衛(wèi)星的SLR精密定軌的研究是很有意義的。

1 GLONASS-M衛(wèi)星及星座特征

24顆衛(wèi)星分布在三個軌道面上，1-8號，9-16號，17-24號分別分布于三個軌道面上，目前GLONASS在軌運行的衛(wèi)星，除了正在測試的GLONASS-k（125）外，均為GLONASS-M星，而GLONASS-M星不同于早期的GLONASS衛(wèi)星和GLONASS-K衛(wèi)星。GLONASS-M衛(wèi)星的基本參數(shù)：帆板面積為30*50cm，反射器個數(shù)112。GLONASS星固系定義為：X軸遠離天底，Z軸沿太陽帆板方向，Y軸反向太陽成右手系。（區(qū)別于GPS35，36 X軸朝向太陽，Y軸沿帆板方向，Z軸朝向天底，成右手系），反射陣列在星固系中的坐標是（-1.874，-0.137，+0.003）.這樣我們就可以根據(jù)反射陣列在星固系中的位置建立GLONASS-M衛(wèi)星的質(zhì)心改正。

2 定軌策略

基于衛(wèi)星動力學原理，采用全球SLR觀測數(shù)據(jù)，計算出GLONASS-M衛(wèi)星的軌道根數(shù)及部分力學參數(shù)，實現(xiàn)其精密定軌。定軌過程中涉及：1）對SLR數(shù)據(jù)做觀測改正，包括地球潮汐改正、板塊運動影響、對流層延遲改正、相對論效應(yīng)、衛(wèi)星激光偏心改正等。2）建立攝動力模型，包括N體攝動、太陽光壓攝動、潮汐攝動、地球形狀攝動和廣義相對論攝動，經(jīng)驗力攝動等。因為GLONASS-M衛(wèi)星屬于中高軌衛(wèi)星，所以可以忽略大氣阻力。為提高定軌精度，實際計算中把部分動力學參數(shù)作為待估計量，與衛(wèi)星軌道參數(shù)一起解。

采用的力學模型：地球非球形攝動（GGM03C模型，取15×15階）；固體潮攝動（IERS2010規(guī)范）；海潮攝動（FES2004模型）；太陽、月亮等引力攝動（采用JPL DE/LE 200大行星歷表中給出的位置）；相對論攝動（IERS2010規(guī)范）；經(jīng)驗RTN攝動；太陽光壓攝動（Box-Wing模型）。

在處理太陽光壓攝動時，我們采用的是Box-Wing 攝動模型，其中也嘗試用類比GPS35，36的光壓攝動模型，建立了GLONASS-M的光壓攝動模型[13]，經(jīng)試驗，結(jié)果沒有采用Box-Wing攝動模型的結(jié)果理想?？赡苁怯捎贕LONASS-M衛(wèi)星與GPS衛(wèi)星存在著一定的差異，所以結(jié)果不理想。GLONASS-M光壓模型的完善也是以后研究應(yīng)該考慮的[13]。

解算參數(shù)：衛(wèi)星初始位置和速度（7天估算一次初軌）（本文也對比了其它定軌弧長）；T、N方向經(jīng)驗力（7天估算一組）；光壓參數(shù)（7天估算一組）；Y方向偏差攝動系數(shù)每7天估算一個系數(shù)。由于SLR觀測數(shù)據(jù)量不多，所以不適宜解算過多的參數(shù)。

3 定軌精度分析及比較

定軌殘差的均方根（即內(nèi)符精度）是評定SLR資料確定的衛(wèi)星軌道精度的常用手段之一，由于內(nèi)符精度所體現(xiàn)的是觀測數(shù)據(jù)的擬合程度，不能定量評估衛(wèi)星軌道的真實精度，因此不能作為評估定軌結(jié)果的唯一和絕對手段，必須通過和其它軌道精度評定方法結(jié)合的方式確保軌道誤差分析的可信度。由于GLONASS衛(wèi)星為導(dǎo)航衛(wèi)星，IGS分析中心會產(chǎn)生相應(yīng)的精密星歷等產(chǎn)品，能提供高精度、連續(xù)的觀測數(shù)據(jù)，實現(xiàn)高精度定軌，其結(jié)果可為評價SLR定軌精度的依據(jù)之一（軌道互比即獨立軌道比較）。另外，重疊弧段比較也是檢驗軌道精度的常用手段之一，因此，本文通過內(nèi)符精度和獨立軌道比較兩種方法可以較準確的分析GLONASS衛(wèi)星SLR定軌精度。

3.1 內(nèi)符精度

衛(wèi)星精密定軌通常是通過定軌殘差的差異判斷其收斂性，用定軌殘差的均方根表示內(nèi)符精度。如用2013年5月19日至6月9日21天的SLR資料分7天弧段對GLONASS123進行定軌，具體每個弧段的內(nèi)符精度.情況見表1。

分析表1，可得出：GLONASS衛(wèi)星的測站數(shù)和觀測值數(shù)目都是比較稀少的，數(shù)據(jù)量較少，然而內(nèi)符精度較高，這并不能說明定軌精度較高，需要結(jié)合獨立軌道檢驗才能全面的評定軌道的精度。

3.2 獨立軌道檢驗

3.2.1 不同弧段對比分析

過短的弧長會因觀測數(shù)據(jù)缺乏導(dǎo)致參數(shù)解算的困難，過長的弧長會因力學模型不精確導(dǎo)致軌道偏離，二者均會導(dǎo)致定軌精度明顯下降。所以選擇一個合適的解算弧段對保證軌道精度是十分重要的。本文以CODE的微波軌道作為精確值，進行精度分析。具體見表2。

分析表2，不難看出：分別采用5天、7天、10天、14天、21天定軌弧長，對GLONASS123衛(wèi)星進行SLR精密定軌徑向精度分別為36cm，20cm，31cm，22cm，28cm。三維方向精度分別為109cm，68cm，106cm，72cm，88cm。對比該21天的定軌數(shù)據(jù)，7天的定軌弧長，定軌的精度最好。徑向精度可以達到20cm，沿軌方向精度可達到40cm，法向可達到50cm，三維方向可達到60cm。考慮到連續(xù)的21天觀測數(shù)據(jù)的質(zhì)量不一，可能有些弧段的定軌精度會相對較低，影響整體的定軌精度。但量級上不會有太大的差異。

3.2.2 不同軌道面衛(wèi)星的精度統(tǒng)計

因為GLONASS衛(wèi)星星座為24顆衛(wèi)星分布在三個軌道面上，為了考慮數(shù)據(jù)的全面性，分別選取了位于三個不同軌道面上的衛(wèi)星數(shù)據(jù)進行軌道計算。其中R07位于第一軌道面，R12位于第二軌道面，R17位于第三軌道面。位于同一軌道面上的衛(wèi)星徑向的系統(tǒng)差相近。（秦顯平，SLR資料精密測定GLONASS衛(wèi)星軌道）。本文也通過對不同軌道面衛(wèi)星的計算驗證了這個結(jié)論。本文選取7天弧段的三顆衛(wèi)星的定軌精度進行了比較分析，具體見表3。經(jīng)對比分析GLONASS的徑向精度最優(yōu)可以達到10cm。

4 結(jié)束語

1）本文針對GLONASS-M衛(wèi)星進行了基于SLR的精密軌道的確定，并且以SHA的GNSS數(shù)據(jù)分析中心的軌道作為標準軌道進行了計算，得到從2013年5月20日至6月9日21天連續(xù)時間序列的定軌結(jié)果，分析在徑向、沿軌方向、法向誤差平均分別可達到37cm和40cm、50cm。徑向最好，法向精度最差，切向次之，是因為 SLR 觀測資料主要反映的是徑向方向的測距值。從徑向的較差圖中可以看出存在明顯的徑向系統(tǒng)差。說明SLR軌道和微波觀測軌道徑向存在系統(tǒng)誤差。

2）針對同一衛(wèi)星采用了不同的定軌弧長進行軌道確定，通過對定軌精度的分析，可以得到針對GLONASS-M衛(wèi)星7天弧長精度略優(yōu)于5、10、14、21天弧長的定軌精度。綜合考慮采用7天的定軌弧長對GLONASS導(dǎo)航衛(wèi)星進行SLR定軌最理想。

3）分別選取了位于三個不同軌道面上的衛(wèi)星數(shù)據(jù)進行軌道計算，得到7天弧段的軌道精度徑向誤差可達到10cm。

4）分別以SHA和CODE的軌道產(chǎn)品為標準軌道為SLR定軌結(jié)果進行精度評定，定軌精度相當。

高精度激光測距資料是檢驗衛(wèi)星軌道精度的重要技術(shù)手段。人衛(wèi)激光測距技術(shù)不僅可以作為軌道檢驗手段，而且還可以作為衛(wèi)星精密定軌的主要手段之一，值得廣泛推廣與應(yīng)用。而且SLR數(shù)據(jù)應(yīng)用于導(dǎo)航衛(wèi)星的精密定軌不僅僅可以檢核導(dǎo)航衛(wèi)星軌道，而且可以將兩種觀測資料聯(lián)合求解導(dǎo)航衛(wèi)星的軌道，可以使定軌精度得到提高[12]。

【參考文獻】

[1]李濟生.人造衛(wèi)星精密軌道確定[M].北京：解放軍出版社，1995.

[2]秦顯平，楊元喜，焦文海，等.利用SLR與偽距資料綜合定軌[J].武漢大學學報，信息科學版，2003，28（6）：745-748.

[3]朱文耀，馮初剛，張華，等.中國SLR網(wǎng)獨立測定ERP和定軌的精度估計[J].測繪學報，1989，18（4）：305-312.

[4]吳斌，彭碧波，許厚澤.地心變化的測定[J].科學通報，1999，44（10）：1106-1108.

[5]朱元蘭，馮初剛，楊福民，等.中國SLR流動站7343，7355地心站坐標的精密測定[J].中國科學院上海天文臺年刊，2003，23：5-11.

[6]何麗娜，王解先，宋淑麗，等.基于SLR時間序列構(gòu)建STRF及其與ITRF2008轉(zhuǎn)換關(guān)系的研究[J].大地測量與地球動力學，2011，31（6）：104-108.

[7]谷德峰，涂先勤，易東云.利用SLR數(shù)據(jù)校準GPS精密定軌系統(tǒng)誤差[J].國防科技大學學報，2008，30（6）：14-18.

[8]孫明國，劉承志，范存波，等.基于SLR精密軌道的天文定位精度分析[J].天文學報，2012，53（2）：153-160.

[9]瞿鋒，王譚強，陳現(xiàn)軍，劉乃苓，程伯輝，等.用SLR資料精密確定GPS35衛(wèi)星軌道[J].測繪學報，1001-1595（2003）03-0224-05.

[10]Claudia Urschl.Gerhard Beutler etal.Contribution of SLR tracking data to GNSS orbit determination[J].Advances in Space Research，2007（39）：1515-1523.

[11]F.Dilssner，T.Springer，G.Gienger，J.Dow etal.The GLONASS-M satellite yaw-attitude model[J].Advances in Space Research，2011（47）：160-171.

[12]秦顯平，楊元喜，王剛，焦文海.SLR資料精密測定GLONASS衛(wèi)星軌道[J].武漢大學學報，1671-8860（2003）04-0440-04.

[13]劉玉瓊，張耀文，吳吉賢.利用SLR檢核GPS35衛(wèi)星精密軌道[J].測繪工程， 2007，16（2）：36-38.

[14]徐克紅，王赫，王永富.基于SLR的GRACE衛(wèi)星定軌中重力場模型對軌道精度的影響[J].測繪工程，2011，20（3）：15-20.

[責任編輯：田吉捷]