陸軼材,高成發(fā),郭 奇
(東南大學(xué)交通學(xué)院,江蘇 南京 210096)
IGS軌道和鐘差產(chǎn)品是精密定位等重要的基礎(chǔ)數(shù)據(jù),也是維持衛(wèi)星導(dǎo)航系統(tǒng)時空基準(zhǔn)的重要產(chǎn)品。目前IGS運行的數(shù)據(jù)處理分析中心不少于12家,主要包括歐洲定軌中心(CODE)、加拿大自然資源部(EMR)、歐洲空間局(ESA)、德國波茲坦地學(xué)研究中心(GFZ)、法國空間研究中心(GRG)、美國噴氣動力實驗室(JPL)、美國麻省理工學(xué)院(MIT)、美國國家大地測量局(NGS)、美國斯克里普斯海洋學(xué)研究所(SIO)、美國海軍天文臺(USN)、捷克大地實驗室(GOP)和武漢大學(xué)(WHU)均能獨立解算,生成精度優(yōu)于5 cm的GNSS精密軌道[1-3]。由于各家分析中心使用的基準(zhǔn)站坐標(biāo)、數(shù)據(jù)處理軟件及定軌策略等不盡相同,定軌結(jié)果會存在明顯的差異[4]。IGS對多個分析中心產(chǎn)品進(jìn)行質(zhì)量分析,然后通過綜合解算生成參考框架和時間基準(zhǔn)穩(wěn)定、統(tǒng)一且1 d精度優(yōu)于2.5 cm的綜合產(chǎn)品[5-7]。
隨著IGS分析中心的發(fā)展,IGS產(chǎn)品綜合的方法不斷得到改進(jìn)和完善。文獻(xiàn)[4]表明分析中心的軌道坐標(biāo)存在系統(tǒng)偏差,且該偏差與極移參數(shù)之差相關(guān);文獻(xiàn)[8]提出了一種利用L1G范數(shù)法求取轉(zhuǎn)換參數(shù)進(jìn)行軌道產(chǎn)品加權(quán)平均的綜合算法。文獻(xiàn)[9—11]探討和分析了IGS產(chǎn)品的一致性問題,提出利用旋轉(zhuǎn)變換和附加改正等方法維持IGS最終軌道/鐘差和周解SINEX等產(chǎn)品的一致性。文獻(xiàn)[12]總結(jié)了IGS軌道產(chǎn)品綜合的3種常用的加權(quán)平均定權(quán)方法:等權(quán)法、基于Helmert轉(zhuǎn)換的驗后方差定權(quán)法及方差因子定權(quán)法。文獻(xiàn)[13]通過選擇最優(yōu)的分析中心組合和框架差異消除方法進(jìn)行精密軌道綜合,極大提升了合成軌道的精度?,F(xiàn)有文獻(xiàn)中常用的定權(quán)方法均是以各分析中心等權(quán)平均軌道為參考軌道,進(jìn)行參數(shù)估計,很少基于分析中心框架、軌道RMS趨于穩(wěn)定的方向?qū)壍谰C合的研究。IGS最終軌道產(chǎn)品在每周四發(fā)布且具有12~18 d的時延性,相比之下分析中心的軌道產(chǎn)品則提前2~3 d,即IGS的最終綜合產(chǎn)品通常需要在分析中心產(chǎn)品全部提交后的2~3 d后才能得到,因此對分析中心的軌道產(chǎn)品進(jìn)行綜合生成精度更優(yōu)的產(chǎn)品是極具工程意義的。故本文提出了基于IGS分析中心產(chǎn)品的軌道綜合算法,并自行編寫程序?qū)PS衛(wèi)星軌道進(jìn)行綜合,為北斗系統(tǒng)中不同類型衛(wèi)星(如GEO/IGSO/MEO)軌道綜合方法提供參考,為iGMAS中心的建設(shè)提供建議[14-16]。
目前IGS僅發(fā)布GPS、GLONASS綜合軌道產(chǎn)品,產(chǎn)品綜合方法有兩種:動力學(xué)法和加權(quán)平均法。第一種方法是將各分析中心產(chǎn)品作為偽觀測數(shù)據(jù)進(jìn)行軌道確定,估計衛(wèi)星初始狀態(tài)和其他動力模型參數(shù),保持了軌道動力學(xué)特性,可靠性較高,但增加了算法復(fù)雜度,且軌道和鐘差綜合無法采用同一算法。第二種方法是將分析中心產(chǎn)品首先轉(zhuǎn)換到同一參考框架,然后進(jìn)行加權(quán)平均得到綜合產(chǎn)品,是IGS目前采用的方法。該方法簡單靈活,軌道和鐘差綜合可采用類似算法,軌道綜合產(chǎn)品在特定條件下滿足動力學(xué)特性,但喪失了軌道動力學(xué)的高階項特性。本文將使用加權(quán)平均法進(jìn)行軌道綜合研究。
雖然分析中心定軌策略與解算軟件不同,導(dǎo)致精度存在明顯差異,但從GPS 1506周后,分析中心軌道精度趨于穩(wěn)定,WRMS曲線趨于平緩(如圖1所示),分析中心的權(quán)值也在整體上趨于穩(wěn)定。
圖1 分析中心軌道相對于IGS綜合軌道的WRMS
基于文獻(xiàn)[2,8,17—20],歸納IGS軌道產(chǎn)品綜合的主要步驟如下:①統(tǒng)一分析中心軌道坐標(biāo)的參考框架,通常選取IERS的極移解(Bulletin A/B)或IGS分析中心的ERP均值作為“真值”,分析中心的軌道坐標(biāo)根據(jù)其極移解與“真值”的差值分別繞X軸和Y軸作相應(yīng)旋轉(zhuǎn);②進(jìn)行Helmert轉(zhuǎn)換,進(jìn)一步統(tǒng)一分析中心的參考框架并確定分析中心的初始權(quán);③基于加權(quán)平均得到組合軌道。
1.1.1 參考框架統(tǒng)一
由于不同分析中心的軌道坐標(biāo)參考框架之間存在旋轉(zhuǎn)、平移等系統(tǒng)性偏差,在軌道綜合之前需先通過旋轉(zhuǎn)變換消除坐標(biāo)參考框架不一致的影響。本文最終軌道綜合選取IGS發(fā)布的周極移解為“真值”,以當(dāng)天中間時刻ERP參數(shù)為參考,計算分析中心極移解與“真值”的偏差,然后對軌道坐標(biāo)作旋轉(zhuǎn)變換。由于極移偏差均為小角度,則旋轉(zhuǎn)變換的公式可按下式簡化
(1)
轉(zhuǎn)之后的分析中心衛(wèi)星坐標(biāo);dXp、dYp分別為各家分析中心極移解相較于IGS發(fā)布的周極移在X、Y方向分量的差值。
1.1.2 Helmert轉(zhuǎn)換定初始權(quán)
在進(jìn)行極移旋轉(zhuǎn)變換后,分析中心衛(wèi)星坐標(biāo)系統(tǒng)間仍存在系統(tǒng)偏差。以IGS發(fā)布的最終軌道為“真值”,將分析中心軌道坐標(biāo)通過Helmert轉(zhuǎn)換統(tǒng)一到相同的坐標(biāo)參考框架消除偏差,并依據(jù)最小二乘準(zhǔn)則,用間接平差法求解參數(shù),誤差方程如下
V3×N,1=B3×N,7·X7,1-L3×N,1
(2)
(3)
基于Helmert旋轉(zhuǎn)后的驗后方差對各分析中心進(jìn)行定權(quán)
(4)
(5)
1.1.3 綜合軌道合成
(6)
(7)
根據(jù)各分析中心變換處理后的坐標(biāo)與綜合軌道間X、Y、Z3個方向坐標(biāo)分量的關(guān)系,可得各分析中心與參考軌道間的RMS及各顆衛(wèi)星的RMS,公式如下
(8)
(9)
式中,mAC為分析中心AC的衛(wèi)星個數(shù)。
按照圖2流程,利用自編的軌道綜合算法,利用2017年2月26日—2017年4月8日(GPS 1938—1940周)的9家分析中心事后精密軌道為試驗數(shù)據(jù),進(jìn)行精密軌道綜合試驗。其中分別以第K周的數(shù)據(jù)作為先驗觀測數(shù)據(jù),進(jìn)行轉(zhuǎn)換參數(shù)的求解與定權(quán)處理,對第K+1周的數(shù)據(jù)進(jìn)行軌道合成。
定權(quán)結(jié)果見表1—表2,為了易于比較和分析,將分析中心ESA的權(quán)值定為單位權(quán)(1.000)。經(jīng)分析得到,分析中心NGS、CODE、GFZ權(quán)值基本大于1,即相對精度較好,剩余分析中心的相對精度則稍差,其中分析中心NGS在所有分析中心的相對權(quán)值最大,即相對精度最高。該結(jié)果與IGS發(fā)布的Weighted Average圖具有一致性,驗證了此算法的可行性(DOY:年積日)。
對GPS 1939周各個分析中心的軌道產(chǎn)品進(jìn)行軌道綜合,圖3為PRN08號衛(wèi)星7天內(nèi)所有歷元的綜合軌道與IGS最終軌道在X、Y、Z方向的差值,其余衛(wèi)星的效果與圖3一致。圖4為綜合軌道、分析中心與IGS最終軌道的單天三維RMS(GPS周1939)。圖5為分析中心綜合軌道相對于IGS最終軌道各顆衛(wèi)星單天平均三維RMS(GPS周1939),其中PRN 04為非健康狀態(tài)。表3給出了合成軌道的單天平均三維RMS(GPS周1939—1943)。
圖2 軌道合成程序流程
ACDOY057058059060061062063064065066067068069070CODE1.0631.0941.5281.5681.4501.4561.7151.8151.5841.1931.4991.5120.9951.231EMR0.3770.7310.6960.6210.6110.4820.7930.6940.5330.3120.5280.5540.4070.529ESA1.0001.0001.0001.0001.0001.0001.0001.0001.0001.0001.0001.0001.0001.000GFZ1.2611.3541.3451.3481.0151.1611.3471.2201.5341.3020.9881.1950.9561.218GRG0.5530.6090.8560.6840.6660.3860.7230.6220.3180.7090.6250.7870.4830.723JPL0.4061.0450.8050.9620.9140.8670.9250.7990.4390.6100.8290.8220.7650.950MIT0.5000.5550.5650.5100.5300.6400.5800.5510.6690.4820.6040.6840.4860.647NGS1.5401.5021.4221.3281.4331.7631.6821.2590.4601.4131.7531.5681.1861.466SIO0.5360.5420.5670.6420.5690.5490.7270.5560.6170.6050.6480.8780.6940.907
表2 分析中心軌道產(chǎn)品定權(quán)(1940—1941周)
圖3 綜合軌道與IGS最終軌道的差值
圖4 綜合軌道、分析中心與IGS最終軌道的單天3D RMS
分析以上數(shù)據(jù)處理結(jié)果,圖4、圖5表明該算法得到的綜合軌道相比于各分析中心的軌道精度有顯著的提高。從圖4可以發(fā)現(xiàn),分析中心軌道與IGS最終軌道比較的平均三維RMS均大于15 mm,即表明分析中心與IGS發(fā)布的綜合軌道間存在框架差異問題。圖5中表明綜合軌道與IGS最終軌道相比,各GPS衛(wèi)星的軌道精度在3~7 mm之間,絕大部分優(yōu)于5 mm。
圖5 綜合軌道各顆衛(wèi)星單天平均3D RMS
表3以周為單位,對綜合軌道的整體精度進(jìn)行統(tǒng)計計算,結(jié)果表明綜合軌道單周精度優(yōu)于4.5 mm,證明了對分析中心軌道產(chǎn)品進(jìn)行加權(quán)綜合的有效性。
表3 合成軌道精度統(tǒng)計(1939—1944)
本文研究了IGS分析中心軌道產(chǎn)品的綜合方法,提出了基于IGS分析中心產(chǎn)品的軌道綜合算法。利用自編算法對IGS各個分析中心2017年2月26日—2017年4月8日精密軌道進(jìn)行綜合,獲得GPS綜合軌道。同時以IGS最終軌道為參考,比較分析了綜合軌道的精度。結(jié)果表明,IGS進(jìn)行軌道綜合的9個分析中心中,NGS、GFZ、CODE和ESA 4個分析中心的軌道產(chǎn)品精度相對較好,剩余5個分析中心的軌道產(chǎn)品相對較差;分析中心軌道與IGS發(fā)布的綜合軌道間存在框架差異性;使用本文算法計算得到的合成軌道結(jié)果與IGS綜合軌道作比較,二者三維差異小于5 mm,優(yōu)于快速星歷,證明了本文算法的可行性。目前,我國將建成多個北斗數(shù)據(jù)分析中心并提供軌道等各類產(chǎn)品,本文研究將為分析中心軌道產(chǎn)品的綜合提供有益的參考[21]。
參考文獻(xiàn):
[1] DOW J M,NEILAN R E,GENDT G.The International GPS Service:Celebrating the 10th Anniversary and Looking to the Next Decade[J].Advances in Space Research,2005,36(3):320-326.
[2] BEUTLER G,KOUBA J,SPRINGER T A.Combining the Orbits of the IGS Analysis Centers[J].Bulletin Géodésique,1995,69(4):200-222.
[3] LI X,GE M,DAI X,et al.Accuracy and Reliability of Multi-GNSS Real-time Precise Positioning:GPS,GLONASS,BeiDou,and Galileo[J].Journal of Geodesy,2015,89(6):607-635.
[4] SPRINGER T A,BEUTLER G.Towards an Official IGS Orbit by Combining the Results of All IGS Processing Centers[C]∥Proceedings of the 1993 IGS Workshop.[S.l.]:IGS,1993:24-26.
[5] GRIFFITHS J.Status of IGS Core Products [C]∥2013 AGU Fall Meeting.[S.l.]:IGS,2013.
[6] 李征航,黃勁松.GPS測量與數(shù)據(jù)處理[M].3版.武漢:武漢大學(xué)出版社,2016.
[7] 高成發(fā),胡伍生.衛(wèi)星導(dǎo)航定位原理與應(yīng)用[M].北京:人民交通出版社,2011.
[8] KOUBA J,MIREAULT Y,LAHAYE F.IGS Orbit/Clock Combination and Evaluation:Appendix I of the Analysis Coordinator Report,International GPS Service for Geo-dynamics (IGS) 1994 Annual Report[M].[S.l.]:Jet Propulsion Laboratory Publication,1994:70-94.
[9] SPRINGER T A,ZUMBERGE J F,KOUBA J.The IGS Analysis Products and the Consistency of the Combined Solutions[C]∥Proceedings of the 1998 IGS Analysis Center Workshop.[S.l.]:IGS,1998.
[10] 姚宜斌.GPS精密定位定軌后處理算法與實現(xiàn)[D].武漢:武漢大學(xué),2004.
[11] 施闖,鄒蓉,姚宜斌,等.基于SINEX解的數(shù)據(jù)組合及系統(tǒng)誤差分析[J].武漢大學(xué)學(xué)報(信息科學(xué)版),2008,33(6):608-611.
[12] 陳康慷.IGS分析中心軌道/鐘差產(chǎn)品綜合及網(wǎng)解模式動態(tài)精密單點定位研究[D].西安:長安大學(xué),2014.
[13] 耿濤,徐夏炎.IGS 分析中心軌道綜合算法實現(xiàn)及精度分析[J].大地測量與地球動力學(xué),2017,37(4):369-373.
[14] STEIGENBERGER P,HUGENTOBLER U,HAUSCHILD A,et al.Orbit and Clock Analysis of Compass GEO and IGSO Satellites[J].Journal of Geodesy,2013,87(6):515-525.
[15] 譚暢,陳國,張強,等.iGMAS軌道產(chǎn)品綜合及精度初步分析[J].武漢大學(xué)學(xué)報(信息科學(xué)版),2016,41(11):1469-1475.
[16] CHEN Kangkang,XU Tianhe,Chen G,et al.The Orbit and Clock Combination of iGMAS Analysis Centers and the Analysis of Their Precision[C]∥中國衛(wèi)星導(dǎo)航學(xué)術(shù)年會.西安:[s.n.],2015.
[17] GRIFFITHS J.A Review of the IGS Final Orbit Combination Procedures (2012),Report from the Analysis Coordinator[EB/OL].[2017-08-07].http:∥acc.igs.org.
[18] WEBER R,F(xiàn)RANGER E.Combination GLONASS Orbits [C]∥IGS International GLONASS Experiment IGEX-98 Workshop.Tennessee:IGS,1999.
[19] 劉偉平,郝金明,于合理,等.導(dǎo)航衛(wèi)星精密軌道與鐘差確定方法研究及精度分析[J].測繪通報,2014(5):5-7.
[20] FERLAND R,PIRASZEWSKI M.The IGS-combined Station Coordinates,Earth Rotation Parameters and Apparent Geo-center[J].Journal of Geodesy,2009,83(3):385-392.
[21] 于素梅.iGMAS產(chǎn)品精度評定和一致性分析[D].西安:長安大學(xué),2016.