劉宇宸，劉 利，周善石，徐君毅，齊 巍，胡小工

(1.中國科學院 上海天文臺，上海200030；2.北京衛(wèi)星導航中心，北京100094；3.中國科學院大學，北京100049；4.北京跟蹤與通信技術研究所，北京100094)

1 引 言

國際地球參考框架(international terrestrial reference frame,ITRF)是國際地球參考系(international terrestrial reference system,ITRS)的實現(xiàn)，通過全球導航衛(wèi)星系統(tǒng)(global navigation satellite system,GNSS)、甚長基線干涉測量(very long baseline interferometry,VLBI)、SLR以及星載多普勒定軌和無線電定位(doppler orbitography and radio-positioning integrated by satellite,DORIS)四種空間大地測量技術建立和維持。2016年1月ITRF發(fā)布最新版國際地球參考框架的ITRF2014[1?3]。

不同的導航系統(tǒng)采用不同的空間基準。如GPS采用WGS84坐標系，GLONASS采用PZ-90坐標系，Galileo采用GTRF坐標系。為保證各衛(wèi)星導航系統(tǒng)間的兼容和互操作，目前各導航系統(tǒng)參考框架采用向ITRF對齊的方式實現(xiàn)導航系統(tǒng)之間空間基準的一致性。對導航系統(tǒng)參考框架的精度可以采用兩種方法進行估算：(1)利用地面站在兩個參考框架的精密坐標之間直接計算轉換參數(shù)；(2)通過包含坐標框架信息的衛(wèi)星軌道產(chǎn)品和定位結果間接監(jiān)測估算框架對齊精度。對于導航系統(tǒng)供應商，一般采用第一種方法獲得與ITRF對齊的精確估算。該方法需要參考框架實現(xiàn)過程中選用的監(jiān)測站數(shù)據(jù)，而該數(shù)據(jù)一般用戶不可獲得。為實現(xiàn)對導航系統(tǒng)參考框架精度的監(jiān)測，一般采用第二種方式，即通過廣播星歷等產(chǎn)品間接估算。兩種方法的監(jiān)測精度和時效性有所區(qū)別：前一種方法可實現(xiàn)毫米級精度估算，一般用于框架實現(xiàn)時的估算；第二種方法精度較第一種低，但更方便用戶對參考框架精度進行監(jiān)測和估算。本文重點分析第二種方法的估算精度，并利用實測數(shù)據(jù)對BDCS開展評估[4]。

北斗衛(wèi)星導航系統(tǒng)的空間基準是北斗坐標系。北斗坐標系采用120余個分布于中國境內及全球的GNSS監(jiān)測站觀測數(shù)據(jù)并實現(xiàn)與ITRF2014的對齊。最新版北斗坐標系通過IGS監(jiān)測站實現(xiàn)與ITRF2014對齊，對齊精度約為毫米級。

除導航系統(tǒng)公布的導航系統(tǒng)參考框架的精度，IGS從用戶角度也提供GPS空間基準的精度監(jiān)測結果。IGS通過比較GPS廣播星歷與IGS快速星歷產(chǎn)品比從而計算7參數(shù)，實現(xiàn)對空間基準對齊精度監(jiān)測，在一定程度上代表了GPS衛(wèi)星對RNSS用戶的空間基準傳遞精度[5?7]。

本文對上述方法的估算精度開展分析，同時結合BDS-3系統(tǒng)導航電文軌道精度較高以及所有衛(wèi)星搭載激光角反射器的特點，給出廣播星歷定位方法以及SLR坐標方法開展對BDS空間基準傳遞精度進行估算。結合中俄空間基準互操作應用背景，利用SLR方法估算了BDS與GLONASS系統(tǒng)空間基準對齊精度，為導航系統(tǒng)兼容互操作應用提供支撐。

2 參考框架實現(xiàn)精度監(jiān)測方法

2.1 廣播星歷定位法

廣播星歷定位法的處理采用固定廣播星歷衛(wèi)星軌道和鐘差，估算GNSS監(jiān)測站每天一組位置參數(shù)、接收機各歷元鐘差、每6 h大氣折射誤差和相位模糊度。式(1)為觀測方程：

式中，PCi為偽距無電離層組合值，LCi為相位無電離層組合值；ρ為星地間幾何距離；tr和ts分別為接收機鐘差和衛(wèi)星鐘差；dtrop為對流層延遲；dother為其他誤差項，如相對論效應，固體潮海潮改正以及多徑噪聲等；λi為無電離層組合頻率載波波長，Ni為組合模糊度未知數(shù)；c為光速，c=299 792 458 m·s?1。為達到對導航系統(tǒng)空間基準傳遞實時監(jiān)測估算的目的，廣播星歷定位法采用廣播星歷計算各歷元衛(wèi)星軌道和鐘差。

利用IGS臺站一天弧長觀測數(shù)據(jù)估算坐標參數(shù)，解算出監(jiān)測站在導航系統(tǒng)參考框架下的坐標。利用上述坐標及該監(jiān)測站在ITRF2014框架下的坐標估算轉換參數(shù)，即得到該參考框架相對于ITRF2014的平移參數(shù)、旋轉參數(shù)和尺度因子。7參數(shù)轉換公式如下：

式中，Xtrans為平移參數(shù)；Xrot為旋轉參數(shù)；Xscale為尺度因子；XGNSS為廣播星歷定位法得到的監(jiān)測站坐標估計值；XITRF為監(jiān)測站在ITRF2014框架下的坐標。

2.2 星歷比對法

廣播星歷空間基準為各導航系統(tǒng)參考框架，精密星歷空間基準主要由IGS維持。目前IGS空間基準與ITRF2014一致，因此可以通過廣播星歷計算導航衛(wèi)星的位置與精密星歷進行比較，獲得導航系統(tǒng)參考框架與ITRF的轉換關系。由于IGS的精密軌道是衛(wèi)星質心坐標，而廣播星歷為衛(wèi)星天線相位中心位置，因此需要進行衛(wèi)星天線相位中心到質心的改正?？紤]到北斗不同類型衛(wèi)星姿態(tài)控制模式不同，衛(wèi)星天線相位中心修正具體如下。

(1)IGSO/MEO衛(wèi)星天線相位中心改正模型：

其中，r,rsun為慣性系中衛(wèi)星位置向量和太陽位置向量；ex,ey,ez為星固系到地固系的轉換矩陣。

(2)GEO衛(wèi)星天線相位中心改正模型：

其中，r,v為慣性系中衛(wèi)星位置向量和速度向量。相位中心改正模型為：

其中，Rciscts為慣性系與地固系轉換旋轉矩陣，xphs,yphs,zphs為星固坐標系下衛(wèi)星天線相位中心偏差，X,Y,Z為衛(wèi)星天線相位中心改正值。

歸算到同一中心后，計算互差并進行7參數(shù)轉換，即得到該參考框架對齊ITRF的平移參數(shù)、旋轉參數(shù)和尺度因子。7參數(shù)轉換公式同公式(2)。

2.3 SLR坐標估算法

SLR坐標估算法與廣播星歷定位法類似，通過將SLR臺站在導航系統(tǒng)參考框架中的坐標與ITRF坐標進行比較，獲得兩個參考框架的轉換參數(shù)。激光數(shù)據(jù)為地面站到衛(wèi)星角反射器之間的測距，具有受電離層折射誤差影響小、測量精度高、且不受鐘差影響等特點。北斗三號衛(wèi)星基于星間鏈路技術提供實時全球高精度廣播星歷產(chǎn)品，軌道精度高[8]。同時，北斗導航系統(tǒng)各衛(wèi)星均配備了角反射器，因此可以采用廣播星歷軌道產(chǎn)品和SLR高精度測量數(shù)據(jù)獲得SLR臺站在BDCS框架下的坐標。將該坐標與SLR臺站在ITRF框架下的坐標進行比較，可獲得兩參考框架的轉換參數(shù)[9,10]。

將廣播星歷計算的軌道由相心坐標改正至質心坐標，方法同上，并進行衛(wèi)星質心改正。采用一個月觀測弧長對SLR臺站坐標進行估算，并將估算坐標與ITRF2014給出的當前時刻坐標進行比較，最后計算7參數(shù)，公式同(2)。

基于中俄合作的背景，對偽距用戶兼容互操作提出了要求。采用SLR坐標估計法對BDSC與PZ-90之間的轉換7參數(shù)進行了比較。

下面給出不同參考框架之間的轉換公式，以BDCS和PZ-90為例。由公式(2)可得到PZ-90與ITRF的轉換公式：

將公式(6)代入到公式(2)中，可得到：

由此可得到BDCS和PZ-90的轉換參數(shù)：

3 BDCS實現(xiàn)精度監(jiān)測

3.1 試驗條件

本文試驗需求的監(jiān)測站需要具備穩(wěn)定接收北斗三號B1I/B3I信號，且歷元可用的北斗三號衛(wèi)星至少4顆。經(jīng)過篩選后選用15個臺站數(shù)據(jù)，分別為境內站wuh2,urum；境外站tuva,zim2,ulab,arht,met3,savo,bor1,pove,ganp,moiu,sthl,pots,sgoc。表1為各站歷元可用北斗三號衛(wèi)星數(shù)量統(tǒng)計。圖1為監(jiān)測站站址分布圖。

表1 歷元可視北斗三號衛(wèi)星數(shù)量統(tǒng)計表

圖1 監(jiān)測站站址分布

3.2 廣播星歷定位方法結果

采用2019年1月到12月全球監(jiān)測站BDS衛(wèi)星B1/B3I雙頻偽距相位30 s采樣數(shù)據(jù)及廣播星歷，計算單天監(jiān)測站坐標。統(tǒng)計了全年的坐標估計值與ITRF坐標互差，并每天計算一組7參數(shù)。

因為北斗三號衛(wèi)星與北斗二號衛(wèi)星空間信號精度存在差異，所以我們僅對北斗三號衛(wèi)星與北斗二號、三號聯(lián)合兩種處理策略相比較。第一種全部使用北斗三號衛(wèi)星，為C19-C37，采用B1I/B3I雙頻組合；第二種采用北斗二號/北斗三號混合星座，為C01-C37，其中北斗二號采用B1I/B3I和B1I/B2I雙頻組合，北斗三號為B1I/B3I雙頻組合。

北斗三號處理各站2019年全年的坐標估計值與ITRF互差序列的統(tǒng)計結果，見表2。

表2 北斗三號廣播星歷定位坐標互差

我們利用2019年2月北斗三號數(shù)據(jù)，通過廣播星歷定位法扣除衛(wèi)星鐘差、接收機鐘差、大氣延遲以及相位模糊度等誤差，并得到O－C時間序列，結果如圖2所示，其平均殘差在?0.2～0.2 m之間。

圖2 北斗三號廣播星歷定位法殘差序列

表3 、圖3給出了北斗三號廣播星歷定位法的轉換7參數(shù)結果。

表3 北斗三號廣播星歷定位法7參數(shù)結果

圖3 北斗三號廣播星歷定位法得到的平移參數(shù)、旋轉參數(shù)和尺度因子隨時間的變化

北斗二號/三號聯(lián)合處理各站2019年全年的坐標值與ITRF互差序列見表4。

表4 北斗二號/三號聯(lián)合廣播星歷定位坐標互差

我們利用2019年2月北斗二號/三號數(shù)據(jù)，通過廣播星歷定位法扣除衛(wèi)星鐘差、接收機鐘差、大氣延遲以及相位模糊度等誤差，并得到O－C時間序列，結果如圖4所示，北斗三號平均殘差在?0.5～0.5 m之間，北斗二號平均殘差在?0.7～0.7 m之間。

圖4 北斗二號/三號衛(wèi)星定位法北二/三衛(wèi)星殘差序列

由于北斗二號是區(qū)域服務系統(tǒng)，沒有境外觀測和注入，因此進行全球評估時采用的是北斗二號長期預報產(chǎn)品。而北斗三號使用星間鏈路觀測，地面站結合星地星間數(shù)據(jù)聯(lián)合定軌，其軌道精度較高。因此相對于北斗二號，北斗三號衛(wèi)星空間信號精度較高，見表5及圖5。

表5 北斗二號/三號聯(lián)合廣播星歷定位法7參數(shù)結果

圖5 聯(lián)合廣播星歷定位法得到的平移參數(shù)、旋轉參數(shù)和尺度因子隨時間的變化

3.3 星歷比對法結果

采用2019年1月3日到1月27日總計25 d的精密星歷及廣播星歷，精密星歷為中國科學院上海天文臺事后處理精密星歷。將廣播星歷與精密星歷相擬合可計算出衛(wèi)星在BDCS框架下和ITRF框架下的坐標，并統(tǒng)計25 d的坐標互差，每天計算一組7參數(shù)。

將上海天文臺北斗事后精密星歷與廣播星歷進行比較，計算單天7參數(shù)，并計算一個月平均值，結果如表6及圖6所示。

表6 BDS廣播星歷與精密星歷7參數(shù)結果

3.4 SLR坐標估算法結果

采用SLR測量數(shù)據(jù)以及廣播星歷，分別計算BDS和GLONASS系統(tǒng)的SLR臺站坐標，結果見表7,8。對于BDS結果，4顆北斗三號衛(wèi)星，共3 330個標準點數(shù)據(jù)；對于GLONASS結果，21顆在軌衛(wèi)星，共4 317個標準點數(shù)據(jù)。

我們采用SLR估算法解得2019年全年的殘差時間序列(見圖7)，表明該方法下BDS的空間信號精度的殘差在?0.15～0.15 m范圍內。

圖6 星歷比對法得到的平移參數(shù)、旋轉參數(shù)和尺度因子隨時間的變化

表7 GLONASS SLR坐標估計互差

表8 2019年BDS SLR坐標估計互差

圖7 SLR坐標估算法殘差序列

采用SLR測量數(shù)據(jù)以及廣播星歷，計算SLR臺站坐標。數(shù)據(jù)為2019年BDS全年ILRS激光聯(lián)測數(shù)據(jù)，由于SLR數(shù)據(jù)采樣較少，一個月估計一組臺站坐標，計算了BDCS與ITRF之間的轉換7參數(shù)。此外，采用GLONASS 2019年1月ILRS激光聯(lián)測數(shù)據(jù)，計算了PZ-90與ITRF，以及BDCS與PZ-90之間的轉換7參數(shù)，結果見表9及圖8。

表9 SLR坐標估算法7參數(shù)結果

圖8 SLR坐標估算法得到的平移參數(shù)、旋轉參數(shù)和尺度因子隨時間的變化

3.5 結果分析

本文采用廣播星歷定位法、廣播星歷與精密星歷比較以及SLR臺站坐標估算三種方法對BDCS實現(xiàn)精度進行計算。第一種方法分別用北斗二號、北斗三號混合星座和北斗三號進行精密坐標估算。北斗三號7參數(shù)轉換結果平移參數(shù)為3.3 cm,?0.7 cm,3.9 cm；旋轉參數(shù)為?2.19 mas,0.08 mas,?0.43 mas；尺度因子為7.6×10?13。北斗二號和北斗三號混合星座7參數(shù)轉換結果平移參數(shù)為7.5 cm,7.3 cm,8.5 cm；旋轉參數(shù)為?2.65 mas,?2.14 mas,0.76 mas；尺度因子為?2.3×10?10。

第二種方法將廣播星歷與精密星歷進行比較。轉換結果平移參數(shù)為?3.4 cm,0.5 cm,3.5 cm；旋轉參數(shù)為0.01 mas,0.4 mas,0.04 mas；尺度因子為?2.6×10?11。

第三種方法采用精密星歷和SLR數(shù)據(jù)計算臺站坐標。BDS轉換結果平移參數(shù)為?4.8 cm,3.4 cm,2.2 cm；旋轉參數(shù)為2.5 mas,?2.6 mas,0.8 mas；尺度因子為4×10?12。GLONASS轉換結果平移參數(shù)為?0.045 cm,1.6 cm,?1.2 cm；旋轉參數(shù)為?1.31 mas,?0.99 mas,?0.94 mas；尺度因子為?1.802×10?9。

BDS與GLONASS轉換結果平移參數(shù)為?0.3 cm,1.7 cm,?6.5 cm；旋轉參數(shù)為?4.51 mas,0.01 mas,?0.85 mas；尺度因子為?2.786×10?9。

通過上述計算結果可以看出，三種方法計算的7參數(shù)的形式誤差均在厘米量級，因此BDCS與ITRF對齊精度在厘米量級。

對比上述三種方法，由于星歷比對法和SLR坐標估算法引入了事后精密產(chǎn)品，監(jiān)測精度比廣播星歷定位法高。北斗衛(wèi)星廣播星歷包含軌道誤差、PCO誤差等，因此所得精度不高，會對參考框架實現(xiàn)精度監(jiān)測造成一定的影響，但隨著廣播星歷精度的提高，其監(jiān)測精度也會逐漸提高。

4 總 結

本文分別給出了三種方法來對BDCS的實現(xiàn)精度進行計算，并采用了2019年全年的數(shù)據(jù)進行試驗，結論如下。

(1)通過上述三種方法對實測數(shù)據(jù)進行計算，均可滿足RNSS米級用戶兼容互操作的需求。廣播星歷定位法可以進行參考框架的實時監(jiān)測，但廣播星歷精度決定了監(jiān)測精度；軌道比對法通過引入事后精密產(chǎn)品來提高監(jiān)測精度，但時效性不強；SLR坐標估算法所得的精度是三種方法中最高的，但目前激光聯(lián)測數(shù)據(jù)量較少，因此只能進行長期的監(jiān)測結果分析。

(2)本文所選取的IGS站均需能夠接收到北斗三號衛(wèi)星B1I/B3I頻點，因此未選取太多站參與本次試驗。北斗二號和北斗三號聯(lián)合處理的方法與僅北斗三號處理的結果在厘米量級一致，因此在后續(xù)的工作中，可以加入更多以IGS核心站為主的地面監(jiān)測站進行全球網(wǎng)解計算，以進行參考框架監(jiān)測，有助于監(jiān)測結果的可靠性。