用單頻GPS數據實現低軌衛星動力學法定軌研究

張德成，鄭作亞，王霞迎，李 偉，3

(1.山東科技大學，山東 青島 266590；2.中國測繪科學研究院，北京 100830；3.蘭州交通大學，蘭州 730070)

1 引言

目前，低軌衛星精密定軌技術已經發展的比較成熟，定軌精度可以達到厘米級水平[1-4]。在利用星載全球定位系統(global positioning system，GPS)觀測數據進行定軌時，許多都是基于雙頻GPS觀測數據解算的，基于單頻觀測數據的定軌研究相對較少[5]。文獻[6]利用HY-2和ZY-3單頻觀測數據，得到的簡化動力學軌道三維精度可以達到1～2 dm[6]，其中用到了半合組合觀測量(group and phase ionospheric correction， GRAPHIC)消去電離層一階項的影響，對采用單頻星載GPS數據進行低軌衛星定軌具有一定參考價值。

本文通過GRACE A衛星的星載GPS觀測數據，分別采用單頻GRAPHIC組合觀測量和雙頻消電離層組合觀測量，運用動力學低軌衛星定軌方法進行解算，證明了利用單頻觀測數據(C/A和L1)進行定軌的可行性，并對定軌結果進行了精度評估。

2 定軌方法

2.1 動力學法低軌衛星定軌

動力學法是比較傳統的定軌方法?；贕PS觀測數據的動力學低軌衛星定軌方法原理為：首先利用加載在低軌衛星上的力學模型及描述低軌衛星的物理模型等，計算出低軌衛星的參考軌道；然后通過最小二乘方法，使動力法獲得的軌道解和GPS跟蹤數據達到最佳擬合，從而得到較為準確的低軌衛星軌道。

在動力學法低軌衛星定軌中，力學模型精度的高低是制約定軌精度的關鍵問題[7]。因此，對于高速運行的低軌衛星來說，必須分析衛星所受到的各種攝動力、精化各種力學模型，才能保證定軌的精度和可靠性。一般將低軌衛星所受攝動力分為保守力和非保守力，保守力包括地球引力攝動、潮汐攝動、N體攝動和相對論效應攝動，非保守力包括太陽光壓攝動、地球反照與紅外輻射壓攝動和大氣阻力攝動等。

衛星在地心慣性系下的運動方程和初始狀態可表示為[8]

(1)

2.2 GRAPHIC觀測模型

在利用GPS雙頻觀測數據進行低軌衛星定軌解算的時候，通常采用無電離層延遲線性組合(PC/LC)的方式消除電離層誤差的影響。而當利用單頻數據(C/A和L1)進行解算時，通常有兩種方法。一種是通過電離層模型加以改正，像全球電離層模型(global ionospheric maps，GIM)對電離層延遲改正效果可以達到90%以上[9]；另一種方法就是通過C/A碼偽距和L1載波相位線性組合的方式，這一方法是由Yunck于1993年首先提出的，稱為GRAPHIC觀測量[10-12]。

對于碼偽距和載波相位觀測量，電離層延遲一階項的影響占到整個電離層延遲的99.9%，并且大小相等，符號相反。因此，通過求碼觀測量和載波相位觀測量的平均值，即采用GRAPHIC線性組合觀測量，就可以消除大部分電離層延遲的影響?？紤]到一般低軌衛星的軌道高度，不考慮對流層延遲的影響。同時，對于多路徑效應、碼偏差和相對論效應也不加考慮。C/A碼偽距和載波相位觀測方程可表示為

C=ρ+ctr-cts+Vion+εC

(2)

Φ=ρ+δtr-δts-Vion+λ1N+εΦ

(3)

式中，C和Φ分別為C/A碼偽距和載波相位觀測量，ρ為GPS衛星到低軌衛星接收機的幾何距離，tr為接收機鐘差，ts為衛星鐘差，Vion為電離層延遲一階項，λ1為L1載波相位波長，N為L1載波整周模糊度，εC和εΦ分別為C/A碼偽距和L1載波相位的觀測噪聲,c為光速。

由此可得到GRAPHIC組合觀測量方程為

(4)

3 實驗算例分析

選取2013-01-11-17GRACE A衛星的GPS跟蹤數據，分別利用單頻和雙頻觀測數據，采用動力學低軌衛星定軌方法，解算了7 d的衛星軌道，并采用兩種方法對定軌結果進行精度評定。

3.1 與RSO軌道進行比較

從德國地學研究中心下載了對應定軌時間段內的快速科學軌道(rapid science orbit， RSO)，以此作為參考軌道分別與兩種策略定軌得到的結果進行比較。得到如下統計結果，圖1和圖2是分別通過GRAPHIC單頻組合觀測量定軌的結果和雙頻消電離層組合觀測量定軌的結果與參考軌道比較得到的殘差曲線圖；圖3和圖4分別是對應的均方根誤差和三維位置誤差統計圖。

圖1 GRAPHIC觀測量解算軌道與RSO軌道比較殘差圖

圖2 雙頻消電離層組合觀測量解算軌道與RSO軌道比較殘差圖

圖3 GRAPHIC觀測量解算軌道與RSO軌道比較誤差統計

圖4 雙頻組合觀測量解算軌道與RSO軌道比較誤差統計

由圖1和圖2可以看出，兩種定軌策略得到的結果在R、T、和N方向上的殘差變化均勻，無明顯系統偏差，但采用GRAPHIC觀測量解算的軌道與參考軌道相比較，殘差較大。由圖3和圖4可以看出，相比雙頻數據，采用單頻數據的定軌結果較差，這主要是因為在GRAPHIC觀測量中雖然消去了電離層延遲的影響，但是碼和載波相位觀測噪聲卻依然存在。由于載波相位觀測噪聲相對較小，因此GRAPHIC觀測量的觀測噪聲約為獨立碼觀測噪聲的一半，這與實際情況是相符的。與RSO軌道相比較，由單頻數據和雙頻數據分別解算得到的定軌結果在R、T、N方向上均方根(root mean square，RMS)的平均值分別為7.9 cm、20.1 cm、5.5 cm和3.5 cm、5.8 cm、3.3 cm，三維定軌精度平均值分別為22.8 cm和9.0 cm。

需要注意的是，在進行動力學定軌時，沒有考慮地球反照與紅外輻射壓和大氣阻力這些難以模型化的攝動力的影響。在優化各種力學模型后，定軌的三維精度可以優于20 cm，能夠滿足一般低軌衛星定軌精度的要求。

3.2 通過衛星激光測距觀測數據進行檢核

通過比較衛星激光測距(satellite laser ranging，SLR)直接測得的站星距與GRACE衛星反算的站星距得到GRACE衛星的軌道檢核結果，可以進一步評價衛星軌道的精度[13]。表1給出了采用SLR數據校核GRACE A衛星軌道的統計結果。圖5為對應的殘差分布圖，其中圖5(a)和圖5(b)分別為雙頻定軌結果和單頻定軌結果的殘差圖。

本算例采用定軌時間段內11個SLR觀測站的402個地殼動態數據信息系統(crustal dynamics data information system，CDDIS)標準點數據，得到軌道檢核結果。對單頻GRAPGIC組合觀測量定軌結果進行檢核，殘差平均值為-1.8 cm，均方根誤差RMS為8.6 cm，標準差為13.7 cm；對采用雙頻消電離層組合觀測量定軌結果進行檢核，殘差平均值為-0.7 cm，均方根誤差RMS為3.3 cm，標準差為5.3 cm?？梢钥闯?，各個觀測站臺不存在較大系統偏差，檢核結果與采用RSO比較的結果一致。

表1 SLR數據檢核GRACE A衛星軌道的比較結果/cm

圖5 GRACE A衛星站星距與SLR測距信息的殘差圖

4 結束語

本文基于單頻GRAPHIC組合觀測量和雙頻消電離層組合觀測量，通過動力學低軌衛星定軌的方法分別進行解算，證明了運用GRAPHIC組合觀測量進行動力學定軌的可行性，并且精度滿足一般低軌衛星定軌的要求。這對只裝載有單頻接收機的低軌衛星定軌，或者對于雙頻接收機在非完好條件下(不能獲得完好的雙頻觀測數據)的低軌衛星定軌，都具有一定的參考價值。

[1] 趙春梅，程鵬飛，益鵬舉.基于偽隨機脈沖估計的簡化動力學衛星定軌方法[J].宇航學報，2011(4)：762-766.

[2] 盛傳貞，甘衛軍，趙春梅，等.利用星載GPS數據確定厘米級Jason-1衛星精密軌道[J].測繪通報，2012(4)：11-14.

[3] 朱俊，王家松，陳建榮，等.HY-2衛星DORIS厘米級精密定軌[J].宇航學報，2013，34(2)：163-169.

[4] 韓保民.基于雙頻P碼的星載GPS低軌衛星相對定軌研究[J].測繪通報，2005(9)：6-9.

[5] 胡國榮，歐吉坤，崔偉宏.星載單頻GPS接收機低軌衛星幾何法定軌研究[J].天文學報，2000，41(3)：278-288.

[6] 郭向，張強，趙齊樂，等.基于單頻星載GPS數據的低軌衛星精密定軌[J].中國空間科學技術，2013(2)：41-46.

[7] 羅志才，鐘波，寧津生，等.GOCE衛星軌道攝動的數值模擬與分析[J].武漢大學學報：信息科學版，2009，34(7)：757-760.

[8] 王威，于志堅.航天器軌道確定：模型與算法[M].北京：國防工業出版社，2007.

[9] 彭冬菊，吳斌.GIM在LEO衛星單頻GPS定軌中的應用[J].天文學報，2012，53(1)：36-50.

[10] MONTENVRUCK O.Kinematic GPS Positioning of LEO Satellites Using Ionosphere-free Single Frequency Measurements[J].Areospace Science and Technology，2003，7(5)：396-405.

[11] BOCK H，JAGGI A，DACH R，et al.GPS Single-frequency Orbit Determination for Low Earth Orbiting Satellites[J].Advances in Space Research，2009，43(5)：783-791.

[12] 宋偉偉，施闖，姚宜斌，等.單頻精密單點定位電離層改正方法和定位精度研究[J].武漢大學學報：信息科學版，2009，34(7)：778-781.

[13] 秦顯平，焦文海，程蘆穎，等.利用SLR檢核CHAMP衛星軌道[J].武漢大學學報：信息科學版，2005，30(1)：92-95.