利用SLR檢核GOCE衛星精密軌道

衷路萍　鄒賢才　吳林沖　于　男

1　武漢大學測繪學院，武漢市珞喻路129號，430079

?

利用SLR檢核GOCE衛星精密軌道

衷路萍1鄒賢才1吳林沖1于男1

1武漢大學測繪學院，武漢市珞喻路129號，430079

摘要：采用2009-11-01～2010-01-31三個月的SLR觀測數據，對GOCE衛星運動學軌道 PKI(precise kinematic orbit)進行檢核?；跉埐罘治霭l現，SLR觀測存在測站時距系統性偏差。消除測站時距偏差后，GOCE衛星PKI精密軌道的外符合精度達到1.5 cm。

關鍵詞：精密軌道；衛星激光測距；外部檢核；時距偏差

GOCE衛星[1]精密軌道產品包括簡化動力學軌道PRD(precise reduced-dynamic orbit)和運動學軌道PKI，均是基于GPS衛星跟蹤，其軌道內符合精度為2 cm[2]。GOCE衛星搭載了激光后向反射器，可實施人衛激光測距。目前SLR測距精度為亞cm級，對于一些觀測條件好的測站甚至可以達到幾mm，成為單次測距精度最高的空間大地測量技術之一[3]。由于SLR的距離觀測獨立于GPS衛星跟蹤，故可利用SLR技術對GOCE衛星精密軌道作外部檢核。利用SLR檢核衛星軌道實際上就是比較SLR得到的星地距離觀測值和通過衛星精密軌道獲得的星地距離計算值之間的差異，兩個距離較差的統計結果就是衛星軌道精度。

1　方法與模型

1.1SLR檢核衛星軌道的方法

SLR技術的原理是觀測激光脈沖信號在地面激光發射站與衛星反射器間的傳播時間，然后乘以光速，從而轉換為地面測站至衛星的距離。其單程距離為:

(1)

式中，ρ為星地距離，c為光速，Δt為激光脈沖信號的雙程傳播時間。

根據地面測站位置和衛星軌道信息，可以計算地面測站至衛星的距離:

(2)

式中，ρ為星地距離，(xs,ys,zs)為衛星(質心)位置矢量，(xi,yi,zi)為第i個測站(系統參考中心SRP)位置矢量。

利用SLR檢核衛星軌道精度的原理就是比較上述星地距離的觀測值和計算值之間的差。為區分兩個距離，將星地距離觀測值記為ρo，計算值記為ρc，兩者之差表示為:

(3)

1.2改正模型

由于SLR觀測在傳播路徑和衛星上有各種誤差，因此需要對SLR原始距離觀測值作數據預處理，包括對流層改正[4]、衛星質心改正[5-6]和廣義相對論改正[4]:

(4)

式中，ρatm為對流層延遲改正，ρcom為衛星質心改正，ρrel為廣義相對論改正。

衛星位置矢量(xs,ys,zs)根據衛星軌道信息內插得到，而測站位置矢量(xi,yi,zi)則由地面測站系統的ITRF坐標歸算至系統參考點SRP(system reference point)，需作的改正包括地殼運動改正、潮汐位移改正[4]和測站偏心改正。各項改正模型如表1所示。

此外，地面測站由于系統設備物理性質、觀測條件等各種因素，在觀測過程中存在時間和距離上的系統誤差[3]，分別為時間偏差Tb和距離偏差Rb。時間偏差指計時器在事件歷元的記錄時刻提前或延遲了時間Tb，距離偏差指測距誤差Rb。在數據后處理分析時，采用時距偏差模型[7-9]估計Tb和Rb兩個參數:

(5)

表1　SLR數據處理中各項改正模型

注：IPIE(The Institute of Precision Instruments Engineering)。

2　計算與分析

GOCE衛星軌道數據采用AIUB(Astronomical Institute of the University of Bern)基于GPS跟蹤的PKI精密軌道。本文利用2009-11-01～2010-01-31三個月內所有GOCE衛星SLR跟蹤的標準點數據，對其軌道作外部檢核。由于GOCE衛星軌道較低、運行速度快，SLR跟蹤較難實現，每次過境只有幾分鐘的觀測，總觀測數據少。在計算時間段內所有的SLR觀測統計如圖1，全球有20個測站觀測，共5 897個標準點數據。由于南北半球測站分布不均勻，南半球只有4個觀測站，幾何結構較差。且各測站觀測能力差異較大，其中南半球Yarragadee(ID: 7090)站觀測標準點數近2 000個，占總觀測量的1/3，少數幾個測站觀測標準點數近900個，而大部分測站觀測標準點數不足300個。

圖1　2009-10-01～2010-01-31期間GOCE衛星SLR觀測數據統計Fig.1　SLR observations of GOCE satellite from 2009-10-01 to 2010-01-31

選取具有代表性的兩個測站——位于北半球的GRAZ站(7839)和位于南半球的Yarragadee站(7090)，各選擇一段降弧和升弧，分析單個觀測弧段的殘差特征?；谏鲜鰯祿A處理模型，取15°截止高度角，測站時距偏差消除前和消除后殘差特性如圖2所示。設置15°截止高度角是因為低高度角的觀測有較大的對流層延遲，帶來較大誤差?？梢钥吹剑?)在大于15°高度角范圍內，SLR觀測存在明顯的系統性誤差，呈線性變化趨勢。這與上述時距偏差模型(式(4))中殘差在一個弧段內線性變化的理論吻合，因此可采用該時距偏差模型擬合殘差，消除測站時距偏差的影響。

圖2　單弧段殘差序列Fig.2　Single pass residuals time series

2)利用式(4)擬合殘差后，殘差量級有明顯減小，基本小于5 cm，線性項趨勢消失，但還存在某些較小的系統性誤差，這可能是由于采用的對流層延遲模型未考慮對流層水平梯度所致[4,10]。3)同測站的升弧和降弧的殘差線性變化趨勢相反(圖2(a)和圖2(b)，圖2(c)和圖2(d))，不同半球同降弧或同升弧的殘差線性變化趨勢相反(圖2(a)和圖2(c)，圖2(b)和圖2(d))，這樣解出的時間偏差參數對于同一測站的升弧和降弧符號也相反。

采用上述數據處理模型，每站每弧段估計一組時距偏差參數。消除測站時距偏差影響后，軌道檢核結果的殘差序列如圖3所示。圖中殘差均小于1 dm，且絕大部分小于5 cm，集中在理論值0上下分布。

圖3　軌道檢核殘差序列Fig.3　Time series of orbit validation residuals

軌道精度檢核的統計結果如表2所示。截止高度角15°，剔除殘差值較大的標準點和觀測數過少的弧段后，總標準點數減少1/4。這主要是因為GOCE衛星軌道過低，較難跟蹤。表2檢核結果顯示，GOCE衛星的PKI軌道無系統性偏差，外符合精度為1.53 cm。與Bock等[2]的檢核結果(平均值0.88 cm，RMS為2.23 cm)相比，本文的結果進一步削弱了系統偏差，RMS提高約0.7 cm。這主要是因為本文相對文獻[2]設置了15°截止高度角，采用了測站時距偏差模型消除測站系統偏差。

表2　GOCE衛星軌道檢核統計結果

為進一步分析截止高度角和時距偏差模型對檢核結果的影響，設計4種實驗方案分別計算各方案中SLR檢核結果，如表3所示。結果顯示，設置一定的截止高度角，檢核精度有小幅度提高；采用時距偏差模型能有效消除系統偏差，且檢核精度有較大提高。表3定量分析了圖2中采用時距偏差模型的影響，也揭示了本文相對于Bock等[2]的檢核結果在系統偏差和精度上提高的主要原因。

表3　截止高度角和時距偏差模型的影響

時距偏差是測站系統性能的一個重要指標。圖4展示了Yarragadee(7090)站時間偏差和距離偏差參數的解算結果，該結果剔除了時間偏差>40 μs的弧段?？梢钥吹?，時間偏差除個別外均穩定在±10 μs內，距離偏差穩定在±5 cm內，其絕對值的均值分別為7.60 μs和1.11 cm，該結果與文獻[11]中的時距偏差結果類似?？傮w來說，測站時距偏差參數變化比較均勻，起伏不大，說明該測站系統性能較穩定。

圖4　Yarragadee站時距偏差變化Fig.4　Variation of station Yarragadee’s time bias and range bias

3　結　語

1)本文利用SLR對GOCE衛星的PKI精密軌道作外部檢核，數據處理過程中發現，檢核殘差存在系統誤差，原因是SLR測距系統存在測站時距偏差。引入時距偏差模型擬合后，殘差的系統誤差明顯減弱。

2)檢核結果表明，GOCE衛星PKI軌道外符合精度為1.53 cm，無系統性偏差。與Bock等[2]的檢核結果相比，進一步削弱了系統性偏差，且精度提高約0.7 cm，主要原因是本文引入了時距偏差模型并消除了測站時距偏差。

3)估算的測站時距偏差參數可作為測站系統質量和性能的參考指標。圖4顯示，Yarragadee測站系統性能整體穩定。

致謝：感謝ILRS數據中心提供SLR觀測數據！

參考文獻

[1]Gruber T, Rummel R, Abrikosov O. GOCE L2 Product Data Handbook[Z]. EGG-C, 2009

[2]Bock H, J?ggi A. GPS-derived Orbits for the GOCE Satellite[J]. Journal of Geodesy, 2011, 85(11): 807-818

[3]Xu G C. Science of Geodesy-I: Advances and Future Directions[M]. Berlin: Springer, 2010

[4]Petit G, Luzum B. IERS Conventions(2010)[R]. Bureau International des Poids et Mesures Sevres, France, 2010

[5]Montenbruck O, Neubert R. Range Correction for the CryoSat and GOCE Laser Retroreflector Arrays[Z].DLR/GSOC, 2011

[6]Bigazzi A, Frommknecht B. Note on GOCE Instruments Positioning-Memo_3.1_[XGCE-GSEG-EOPG-TN-09-0007 v3.1][Z]. ESA, 2010

[7]馮初剛, 朱元蘭. 1992年LAGEOS衛星全球SLR標準點資料的分析[J]. 中國科學院上海天文臺年刊, 1995, 16: 107-112(Feng Chugang, Zhu Yuanlan. Analysis for LAGEOS SLR Global Normal Point Data of 1992[J]. Annals of Shanghai Observatory Academia Sinica, 1995, 16: 107-112)

[8]馮初剛, 朱元蘭, 張飛鵬. LAGEOS衛星精密定軌及殘差分析[J]. 天文學報, 2003, 44(1): 55-64(Feng Chugang, Zhu Yuanlan, Zhang Feipeng. Determination of LAGEOS Satellite’s Precise Orbits and Residual Analysis[J]. Acta Astronomica Sinica, 2003, 44(1): 55-64)

[9]祝芙英, 吳云, 杜瑞林, 等. 流動SLR烏魯木齊站觀測精度評估[J]. 大地測量與地球動力學, 2008, 28(2): 127-131(Zhu Fuying, Wu Yun, Du Ruilin, et al. Accuracy of Mobile SLR at Urumchi Station[J]. Journal of Geodesy and Geodynamics, 2008, 28(2): 127-131)

[10]Hulley G C, Pavlis E C. A Ray-tracing Technique for Improving Satellite Laser Ranging Atmospheric Delay Corrections, Including the Effects of Horizontal Refractivity Gradients[J].Journal of Geophysical Research: Solid Earth, 2007, 112(B6)

[11]趙春梅, 瞿鋒, 程鵬飛, 等. 阿根廷圣胡安激光測距系統的SLR數據質量分析[J]. 測繪學報, 2008, 37(3): 338-341(Zhao Chunmei, Qu Feng, Cheng Pengfei, et al. Data Quality Analysis of Argentina San Juan Laser Ranging System[J]. Acta Geodaetica et Cartographica Sinica, 2008, 37(3): 338-341)

Foundation support:National Key Basic Research Program of China, No. 2013CB733301; National Natural Science Foundation of China, No. 41004007, 41274033; Special Fund for Public Welfare Industry (Surveying,Mapping and Geoinformation), No. 201512001; Fundamental Research Funds for the Central Universities, No. 2015214020202.

About the first author:ZHONG Luping, postgraduate, majors in satellite laser ranging, E-mail: lpzhong@whu.edu.cn.

收稿日期：2015-08-14

第一作者簡介：衷路萍，碩士生，主要研究方向為人衛激光測距，E-mail：lpzhong@whu.edu.cn。 通訊作者：鄒賢才，博士，副教授，主要研究方向為衛星大地測量和物理大地測量，E-mail：xczou@whu.edu.cn。

DOI：10.14075/j.jgg.2016.08.014

文章編號：1671-5942(2016)08-0719-04

中圖分類號：P228

文獻標識碼：A

Corresponding author:ZOU Xiancai, PhD, associate professor, majors in satellite geodesy and physical geodesy, Email: xczou@whu.edu.cn.

Evaluation of GOCE Precise Orbit Using SLR

ZHONGLuping1ZOUXiancai1WULinchong1YUNan1

1School of Geodesy and Geomatics, Wuhan University, 129 Luoyu Road, Wuhan 430079, China

Abstract:In this paper we present the method of satellite orbit validation with SLR measurements, and use SLR data observed from 1 November 2009 until 31 January 2010 (3 months) to validate the GOCE PKI (precise kinematic) orbit. Based on residual analysis, we find that system error due to station time bias and range bias exits and has not been removed in the SLR observation. After modeling and fitting the station time bias and range bias, we achieve the GOCE PKI orbit precision at the level of 1.5 cm.

Key words:precise orbit; SLR; validation; station time bias and range bias

項目來源：國家973計劃(2013CB733301)；國家自然科學基金(41004007, 41274033)；測繪地理信息公益性行業科研專項(201512001)；中央高?；究蒲袠I務費(2015214020202)。