崔紅正唐歌實宋柏延劉薈萃葛茂榮

(1北京航天飛行控制中心航天飛行動力學(xué)技術(shù)重點實驗室,北京100094)

(2北京航空航天大學(xué)宇航學(xué)院,北京100191) (3德國地學(xué)研究中心,波茨坦14473)

北斗衛(wèi)星導(dǎo)航系統(tǒng)實時定軌與鐘差處理策略

崔紅正1,2唐歌實1宋柏延1劉薈萃1葛茂榮3

(1北京航天飛行控制中心航天飛行動力學(xué)技術(shù)重點實驗室,北京100094)

(2北京航空航天大學(xué)宇航學(xué)院,北京100191) (3德國地學(xué)研究中心,波茨坦14473)

目前鮮有對北斗衛(wèi)星導(dǎo)航系統(tǒng)(BeiDou Navigation Satellite System,BDS)實時精密定軌與鐘差確定的研究,文章提出了BDS實時軌道與實時鐘差處理策略,包括了觀測與動力學(xué)模型、實時軌道與實時鐘差處理流程與評估方法。尤其對于實時鐘差,為了提高計算效率,聯(lián)合使用兩個獨立并行的線程估計非差絕對鐘差和歷元間相對鐘差。利用多模全球衛(wèi)星導(dǎo)航系統(tǒng)試驗(MGEX)與全球連續(xù)檢測評估系統(tǒng)(iGMAS)實測數(shù)據(jù)進行了北斗實時軌道與鐘差解算,BDS實時軌道徑向平均精度對于GEO衛(wèi)星優(yōu)于20 cm,對于IGSO與MEO一般優(yōu)于10 cm;鐘差精度對于GEO衛(wèi)星為0.5～4.5 ns,對于IGSO/MEO為0.2～2.0 ns?；谀壳暗能壍琅c鐘差結(jié)果,實時精密單點定位(Precise Point Positioning,PPP)結(jié)果可以達到分米量級。

北斗衛(wèi)星導(dǎo)航系統(tǒng);實時確定;衛(wèi)星鐘差;精密單點定位;全球連續(xù)檢測評估系統(tǒng);衛(wèi)星軌道

1 引言

截止到2012年底,北斗衛(wèi)星導(dǎo)航系統(tǒng)(BDS,本文簡稱“北斗”系統(tǒng))共有14顆衛(wèi)星,并已經(jīng)開始提供亞太區(qū)域服務(wù),計劃到2020年全球衛(wèi)星星座部署完畢,開始提供全球?qū)Ш蕉ㄎ环?wù)。為開展全球范圍全球衛(wèi)星導(dǎo)航系統(tǒng)(GNSS)開放服務(wù)的監(jiān)測評估,2011年,中國在ICG(International Committee on GNSS)大會上提出了建立iGMAS系統(tǒng)的建議[1]。經(jīng)初步論證,全球連續(xù)監(jiān)測評估系統(tǒng)包括30個全球跟蹤站(國外21個、國內(nèi)9個,大部分支持四大系統(tǒng))、3個數(shù)據(jù)中心、8個分析中心、1個監(jiān)測評估中心、1個產(chǎn)品綜合與服務(wù)中心、1個運行控制管理中心和通信鏈路等。國際GNSS服務(wù)組織(IGS)也一直關(guān)注中國“北斗”系統(tǒng)的建設(shè),IGS發(fā)起了MGEX(Multi GNSS Experiment)實驗,并在2012年召開了以Multi-GNSS為主體的工作組,目前有近60個跟蹤站具“北斗”系統(tǒng)數(shù)據(jù)跟蹤能力,對本文的研究提供了重要的數(shù)據(jù)源。隨著“北斗”系統(tǒng)空間段與地面段建設(shè)的日趨完善,為了改善“北斗”系統(tǒng)在亞太區(qū)域的定位精度、可用性、連續(xù)性以及完好性,我國也開始籌建“北斗”系統(tǒng)廣域增強系統(tǒng),其核心是生成高精度實時衛(wèi)星軌道、鐘差、電離層產(chǎn)品信息,本文正是基于此開展“北斗”系統(tǒng)實時衛(wèi)星軌道與衛(wèi)星鐘差確定技術(shù)的研究。

2000年JPL研制的實時GPS衛(wèi)星定軌與鐘差確定軟件RTG(Real Time GIPSY)[2],基于美國NASA全球GPS網(wǎng)絡(luò)(GGN)60個左右觀測站的實時數(shù)據(jù),可以提供GPS衛(wèi)星準實時與實時精密軌道與鐘差產(chǎn)品。2002年IGS成立了實時產(chǎn)品工作組(RTWG),為促進GNSS實時產(chǎn)品的應(yīng)用,開始了原型實時產(chǎn)品生成系統(tǒng)的研究。在IGS組織2002—2007戰(zhàn)略計劃中,實時產(chǎn)品處理被確定為GNSS高精度數(shù)據(jù)處理技術(shù)發(fā)展趨勢之一。文獻[3]在Bernese軟件的基礎(chǔ)上,對軌道、鐘差、地球自轉(zhuǎn)參數(shù)、測站坐標以及大氣參數(shù)進行了準實時處理分析。文獻[4]采用擴展卡爾曼濾波開發(fā)了實時衛(wèi)星軌道和鐘差生成的Auto-BAHN系統(tǒng),實時衛(wèi)星軌道鐘差更新率為5 min,實時軌道精度15 cm左右,鐘差精度為0.3 ns。文獻[5]開發(fā)了實時精密軌道和鐘差確定系統(tǒng)(RTOD),快速預(yù)報軌道作為實時軌道,精度為20 cm。德國地學(xué)研究中心(GFZ)開發(fā)了EPOS-(RT)軟件,實時定軌精度為5 cm以下,鐘差精度為0.1 ns[6]。國內(nèi)方面,導(dǎo)航衛(wèi)星定軌技術(shù)研究開始于20世紀80年代,伴隨GPS進入中國,許多學(xué)者開始對導(dǎo)航衛(wèi)星精密定軌展開研究。雖然近年來GNSS數(shù)據(jù)分析工作在國內(nèi)取得了長足的進度,如武漢大學(xué)基于實時軌道與鐘差產(chǎn)品[7-8],開展了中國及周邊區(qū)域的廣域?qū)崟r精密定位原型系統(tǒng)的研究與建設(shè),但實時GNSS數(shù)據(jù)處理工作尚處于起步階段,尤其對于具有自主知識產(chǎn)權(quán)的“北斗”系統(tǒng)。

本文給出了實時軌道與鐘差確定策略,利用iGMAS、MGEX實時數(shù)據(jù)流,對全球?qū)崟r“北斗”系統(tǒng)衛(wèi)星軌道與鐘差處理結(jié)果進行了分析,并給出了實時精密單點定位(Precise Point Positioning, PPP)結(jié)果。

2 實時軌道與鐘差確定策略

導(dǎo)航衛(wèi)星經(jīng)過長期的模型精化,軌道動力學(xué)模型已經(jīng)非常成熟,短期軌道預(yù)報達到很高的精度,本文實時軌道每次處理繼承之前70 h的觀測弧段,結(jié)合2 h的最新觀測弧段,合并成72 h觀測弧段(3天弧段)利用最小二乘估計器進行軌道確定,并預(yù)報2 h弧段作為實時軌道產(chǎn)品。

實時“北斗”系統(tǒng)衛(wèi)星鐘差處理流程如圖1所示,為了提高計算效率,尤其是在實時應(yīng)用中,精密鐘差實時估計采用固定衛(wèi)星軌道和站坐標參數(shù),僅估計衛(wèi)星鐘差、接收機鐘差、對流層和模糊度參數(shù)的方式。為了進一步提高計算效率采用非差模型估計與歷元差分模型估計兩個線程相結(jié)合的模式,對于非差模式,采用無電離層的非差相位相位(LC)和偽距組合(PC)觀測值,同時估計的參數(shù)有衛(wèi)星鐘差、接收機鐘差、實數(shù)無電離層組合模糊度參數(shù)、對流層參數(shù):

對于非差模式的處理模式,大量待估參數(shù)影響了實時鐘差的處理速度,此處采用參數(shù)消除方法來降低法方程中參數(shù)的維數(shù)[9],由于鐘差逐歷元變化,在進入法方程之前將其預(yù)先消除,待模糊度參數(shù)、對流層參數(shù)求解完成后,然后恢復(fù)出消除的鐘差參數(shù)。為了進一步提高處理效率,采用如下歷元差分估計模式:

此處消去了模糊度參數(shù),進而有效提高處理效率。聯(lián)合使用兩個獨立并行的線程估計非差和歷元相對鐘差,非差線程估計5 s采樣率鐘差,而在歷元差分線程中1 s估計一組相對衛(wèi)星鐘差,即得到高效率生成高采樣率實時鐘差產(chǎn)品。

圖1 “北斗”系統(tǒng)衛(wèi)星實時鐘差處理流程Fig.1 Diagram of BDS satellite real time clock determination

固定利用后處理PPP解算或從SINEX(Solution INdependent EXchange Format)文件中讀取的測站坐標,并固定預(yù)報實時軌道,利用“北斗”系統(tǒng)廣播星歷作為待估鐘差的初值,采用均方根信息濾波器(SRIF)估計器[10]對實時數(shù)據(jù)流進行處理,得出“北斗”系統(tǒng)實時鐘差產(chǎn)品,具體估計策略如表1所示。由于網(wǎng)絡(luò)延遲與數(shù)據(jù)處理延遲,實時鐘差產(chǎn)品有幾秒至幾十秒的延遲,需要對其進行超短期預(yù)報。此處,由于鐘差估計時間基準的不穩(wěn)定性,對歷元間鐘差差分序列進行二次多項式擬合,求出歷元間差分預(yù)報值,之后再還原成預(yù)報鐘差。

表1 實時 “北斗”系統(tǒng)衛(wèi)星軌道與衛(wèi)星鐘差估計策略Tab.1 Processing strategies for BDS real time orbit determination and clock determination

續(xù)表1

對于實時軌道與鐘差,利用GFZ解算的“北斗”系統(tǒng)最終軌道與鐘差產(chǎn)品進行外符合精度對比,并利用仿真實時PPP結(jié)果驗證“北斗”系統(tǒng)實時軌道與鐘差的解算效果。用1D RMS來評價定軌結(jié)果,1D RMS=((RM+RM+RM)/3)1/2,其中,RMSR為徑向RMS,RMSC為法向RMS, RMSA為跡向RMS。

3 算例測試與結(jié)果評估

算例測試結(jié)果如圖2～圖7所示?！氨倍贰毕到y(tǒng)實時軌道處理采用MGEX與iGMAS一共55個站的“北斗”系統(tǒng)觀測數(shù)據(jù),如圖2所示,計算時間為“北斗”系統(tǒng)時間(BDT):2015年第060天到066天。2 h預(yù)報軌道與GFZ的最終“北斗”系統(tǒng)軌道產(chǎn)品[12]對比。軌道精度對于GEO衛(wèi)星(C01C02C04C05)為100～400 cm(1D RMS),對于IGSO(C06C07C08C09C10)為80 cm(1D RMS)以下,MEO(C11C12C14)為30 cm(1D RMS)以下,如圖4所示。對于定位精度,主要考察軌道徑向精度,BDS實時軌道徑向平均精度如圖5所示,GEO衛(wèi)星為20 cm以下,對于IGSO與MEO為10 cm以下。

“北斗”系統(tǒng)實時鐘差處理采用如圖1所示的策略,MGEX與iGMAS一共19個站具有“北斗”系統(tǒng)實時流觀測數(shù)據(jù),如圖3所示,計算時間為“北斗”系統(tǒng)時間:2015年第060天到066天,實時“北斗”系統(tǒng)鐘差與GFZ的最終“北斗”系統(tǒng)鐘差產(chǎn)品[12]對比。鐘差精度對于GEO衛(wèi)星為0.5～4.5 ns,對于IGSO與MEO為0.2～2.0 ns以下,如圖6所示。

實時PPP結(jié)果用來驗證“北斗”系統(tǒng)實時軌道與鐘差的解算效果,利用北京航天飛行控制中心自己維護的多模GNSS跟蹤站BJSH站的實時數(shù)據(jù),進行仿真實時PPP解算。計算時間為“北斗”系統(tǒng)時間:2015年第060天到066天,利用GPS靜態(tài)PPP結(jié)果作為測站參考坐標,每天重新啟動,24 h后的定位結(jié)果如圖7所示:水平東方向(East)低于5 cm,水平北方向(North)低于10 cm,垂直向上方向(Up)低于10 cm,但由于063天“北斗”系統(tǒng)鐘差缺少部分衛(wèi)星,導(dǎo)致定位精度稍微偏差。

圖2 用于BDS實時軌道解算的跟蹤站分布Fig.2 Ground tracking stations distribution for real time orbit determination

圖3 用于BDS實時鐘差解算的跟蹤站分布Fig.3 Ground tracking stations distribution for real time clock determination

圖4 BDS實時軌道與GFZ的最終北斗軌道產(chǎn)品對比(1D RMS)Fig.4 1D RMS of BDS real time orbit compared with GFZ final orbit

圖5 BDS實時軌道徑向平均精度統(tǒng)計Fig.5 Mean radial RMS of BDS real time orbit compared with GFZ final orbit

圖6 BDS實時鐘差精度統(tǒng)計Fig.6 STD of BDS real time clock compared with GFZ final clock

圖7 基于BDS實時軌道與實時鐘差的實時PPP精度Fig.7 Real time PPP accuracy for BJSH station

4 結(jié)束語

本文針對“北斗”系統(tǒng)廣域增強系統(tǒng)建設(shè)的需要,提出了“北斗”系統(tǒng)實時軌道與實時鐘差處理流程,尤其對于實時鐘差,為了提高計算效率聯(lián)合使用兩個獨立并行的線程估計非差絕對鐘差和歷元相對鐘差。利用MGEX與iGMAS實測數(shù)據(jù)進行了“北斗”系統(tǒng)實時軌道與鐘差解算,BDS實時軌道徑向平均精度:GEO衛(wèi)星為20 cm以下,對于IGSO與MEO為10 cm以下;鐘差精度對于GEO衛(wèi)星為0.5～4.5 ns,對于IGSO與MEO為0.2～2.0ns以下?；谀壳暗能壍琅c鐘差結(jié)果,實時PPP結(jié)果可以達到分米量級。結(jié)果表明本文提出的處理流程與策略可很好完成“北斗”系統(tǒng)廣域增強系統(tǒng)對核心實時軌道與鐘差的需求,可以用于北斗二代重大專項中國廣域增強系統(tǒng)關(guān)鍵技術(shù)攻關(guān)與演示驗證。后續(xù),隨著“北斗”系統(tǒng)空間段與地面段加速建設(shè),實時軌道與鐘差精度能夠進一步提高。為了滿足不同導(dǎo)航系統(tǒng)的互操作,需要進一步對多系統(tǒng)的實時產(chǎn)品聯(lián)合解算效果進行分析。

致謝 作者感謝國際GNSS監(jiān)測評估系統(tǒng)(iGMAS)提供的觀測數(shù)據(jù)(www.igmas.org)。

[1] YANG.Working progress of iGMAS[C].IGS Workshop,Olsztyn,Poland,2012.

[2] SHARP T,HATCH R,NELSON F.Star Fire and real-time GIPSY:A high accuracy Global Differential GPS System[C].Proceedings of the 5th International Symposium on Satellite Navigation Technology& Application,Canberra,Australia,July 24-27,2001.

[3] DACH R,HUGENTOBLER U,FRIDEZ P,et al.Bernese GPS software version 5.0 user manual[R]. Bern:University of Bern,2007.

[4] ZHANG Q,MOORE P,HANLEY J,et al.Auto-BAHN:Software for near real-time GPS orbit and clock computations[J].Advanced in Space Research,2007,39(10):1531-1538.

[5] SCOTT L.Real time precise orbit&clock determination[R].UK:Thales Research and Technology Limited,2005.

[6] GE M,CHEN J,GENDT G.EPOS-RT:Software for real-time GNSS data processing[C].EGU General Assembly,Vienna,Austria,2009.

[7] 施闖,樓益棟,宋偉偉,等.廣域?qū)崟r精密定位原型系統(tǒng)及初步結(jié)果[J].武漢大學(xué)學(xué)報 (信息科學(xué)版), 2009,30(11):1271-1274.

SHI CHUANG,LOU YIDONG,SONG WEIWEI,et al.A wide area real-time differential GPS prototype system and the initial results[J].Geomatics and Information Science of Wuhan University,2009,30(11): 1271-1274.

[8] 樓益棟,施闖,周小青,等.GPS精密衛(wèi)星鐘差估計與分析[J].武漢大學(xué)學(xué)報(信息科學(xué)版),2009,34(1): 88-91.

LOU YIDONG,SHI CHUANG,ZHOU XIAOQING,et al.Realization and analysis of GPS precise clock products[J].Geomatics and Information Science of Wuhan University,2009,34(1):88-91.

[9] 葛茂榮.GPS衛(wèi)星精密定軌理論及軟件研究[D].武漢:武漢測繪科技大學(xué),1995.

GE MAORONG.Research on GPS precise orbit determination theory and software[D].Wuhan:Wuhan Technological University of Survey and Mapping,1995.

[10] BIERMAN G J.Factorization methods for discrete sequential estimation[M].New York:Academic Press,1997.

[11] RALF S,PETER S,GERD G,et al.Generation of a consistent absolute phase-center correction model for GPS receiver and satellite antennas[J].Journal of Geodesy,2007,81(12):781-798.

[12] DENG ZHIGUO.Precise orbit determination of BeiDou satellites at GFZ[C].IGS Workshop,Pasadena, USA,2014.

作者簡介

崔紅正 1983年生,2012年獲北京航空航天大學(xué)飛行器設(shè)計專業(yè)博士學(xué)位,現(xiàn)為北京航天飛行控制中心博士后。研究方向為GNSS高精度數(shù)據(jù)處理、衛(wèi)星星座的優(yōu)化設(shè)計、航天器精密軌道確定、航天器姿態(tài)軌道動力學(xué)與控制、航天任務(wù)分析與設(shè)計。

BDS Satellite Real Time Orbit and Clock Determination and Initial Results Analysis

CUI Hongzheng1,2TANG Geshi1SONG Baiyan1LIU Huicui1GE Maorong3
(1 Flight Dynamics Laboratory,Beijing Aerospace Control Center,Beijing 100094)
(2 School of Astronautics,Beijing University of Aeronautics and Astronautics,Beijing 100191)
(3 German Research Center for Geosciences,Potsdam 14473)

Navigation satellite real time orbit determination and clock determination are research highlights in navigation areas.The processing strategies and processes for Beidou Navigation Satellite System(BDS)real time orbit determination and clock determination were provided,especially for the clock processing with zero-differenced absolute clock and epoch-differenced relative clock.Using the data from the Multi-GNSS Experiment(MGEX)and international GNSS monitoring&Assessment Service(iGMAS)tracking networks,the average orbit accuracy for real time BDS in radial direction is better than 20 cm,and 10 cm for GEO,IGSO and MEO,respectively.And the clock accuracy for real time BDS can reach 0.5～4.5 ns,and 0.2～2.0 ns for GEO,IGSO and MEO,respectively. Based on the this result,the real time precise pointing position(PPP)performance following BDS can achieve decimeter level.

BDS;Real time determination;Satellite clock;Precise point positioning; iGMAS;Satellite orbit

10.3780/j.issn.1000-758X.2015.05.001

(編輯:車曉玲)

中國博士后科學(xué)基金第52批面上(2012M521883),全球連續(xù)監(jiān)測評估系統(tǒng)分析中心建設(shè)與運行維護(GFZX0301040308-05)資助項目

2015-01-07。收修改稿日期:2015-05-20