基于iGMAS北斗新一代卫星精密定轨初步结果与分析

1、基于iGMAS观测网络数据的北斗新一代卫星（BeiDou-3）精密定轨初步结果与分析潭冰峰，袁运斌第二届大地第二届大地测量与地球物理学学量与地球物理学学术大会大会-IAG-G02江江苏 南京南京中国科学院测量与地球物理研究所2016年9月24日Institute of Geodesy and Geophysics (IGG) Slide 2报告目录报告目录p研究背景研究背景pBeiDou-3卫星精密定轨方法与流程卫星精密定轨方法与流程pBeiDou-3卫星精密定轨初步结果展示卫星精密定轨初步结果展示p总结与展望总结与展望Institute of Geodesy and Geophysics (

2、IGG) Slide 3研究背景研究背景(1/5)北斗卫星导航系统北斗卫星导航系统-BeiDou Navigation Satellite System (BDS)BDS系统建设系统建设区域系统 全球系统BDS区域系统情况区域系统情况5 GEO + 5 IGSO + 4 MEO；2012年12月27日起正式向亚太区域提供定位、导航、授时服务；BDS全球系统情况全球系统情况2015年3月30日，第一颗新一代北斗全球导航卫星（BeiDou-3）- BeiDou I1S-C31发射成功；标志着北斗卫星导航系统正式从区域走向全球化建设；2015年2月10日起，国际海事组织（IMO）正式认可BDS并将其

3、纳入全球无线电导航系统。目前在轨运行BeiDou-3 Satellites： C31, C32, C33, C34, C35；Institute of Geodesy and Geophysics (IGG) Slide 4研究背景研究背景(2/5)北斗卫星导航系统星座现状（北斗卫星导航系统星座现状（2016-9）Beidou-2SatellitePRNSVNTypeNotesLaunch DateG01C01C003GEO140.0 EJan 17, 2010G02C02C002GEONot in operationApr 15, 2009G03C03C004GEO84.0 EJan 2,

4、2010G04C04C006GEO160.0 ENov 1, 2011G05C05C011GEO58.75 EFeb 25, 2012G06C02C016GEO110.5 EOct 25, 2012G07C17C017GEOIn orbit testJun 12, 2016I01C06C005IGSO122 EAug 1, 2010I02C07C007IGSO119 EDec 18, 2010I03C08C008IGSO120 EApr 10, 2011I04C09C009IGSO96.5 EJul 27, 2011I05C10C010IGSO92.5 EDec 2, 2011I06C15C0

5、17IGSO95 EMar 29, 2016M01C30C001MEONot in operationApr 4, 2007M03C11C012MEO Apr 30, 2012M04C12C013MEO Apr 30, 2012M05C13C014MEONot in operationSep 18, 2012M06C14C015MEO Sep 18, 2012Beidou-3SatellitePRNSVNTypeNotesLaunch DateI1-SC31C101IGSO Mar 30, 2015M1-SC33C102MEO Jul 25, 2015M2-SC34C103MEO Jul 25

6、, 2015I2-SC32C104IGSO Sep 29, 2015M3-SC35C105MEOIn orbit testFeb 1, 2016* /web/eoportal/satellite-missions/c-missions/cnssInstitute of Geodesy and Geophysics (IGG) Slide 5研究背景研究背景(3/5)BeiDou-3 Satellites新特性新特性保留保留BeiDou-2 信号传输体系信号传输体系三个频点：1561.098 MHz (B1), 1207.140 MHz

7、(B2) and 1268.520 MHz (B3)；新的信号体制新的信号体制B1变更为1575.42 MHz，与GPS L1及Galileo E1频率一致；B1信号调制方式从QPSK变为MBOC；B2信号调制方式采用新的四种调制方式；地面跟踪站地面跟踪站iGMAS：目前有9个跟踪站能接收BeiDou-3 Satellites 旧的B1及B3频点数据；MGEX： 目前有4个跟踪站能接收BeiDou-3 C31 B1频点数据；B2a及B2b两种备选频率；两种B3c备选频率；Institute of Geodesy and Geophysics (IGG) Slide 6研究背景研究背景(4/5)

8、iGMAS BDS跟踪站跟踪站-9个个Institute of Geodesy and Geophysics (IGG) Slide 7研究背景研究背景(4/5)iGMAS BDS跟踪站跟踪站-9个个Station Abb.LocationCountryReceiverAntennaBJF1BeijingChinaCETC-54 GMR-4011LEIAR25.R4 LEITBRCHBraunschweigGermanyCETC-54 GMR-4011LEIAR25.R4 LEITCLGYCalgaryCanadaCETC-54 GMR-4011LEIAR25.R4 LEITGUA1Urumq

9、iChinaGNSS_GGRRINT-8CH CETDLHA1LhasaChinaCETC-54 GMR-4011NOV750.R4 NOVSTAHTTahitiFrance GNSS_GGRRINT-8CH CETDWUH1WuhanChinaCETC-54 GMR-4011LEIAR25.R4 LEITXIA1XianChinaGNSS_GGRRINT-8CH CETDXHONAntarcticaUnited NationsCETC-54 GMR-4011LEIAR25.R4 LEITInstitute of Geodesy and Geophysics (IGG) Slide 8研究背景

10、研究背景(5/5)MGEX BDS跟踪站跟踪站-52个，其中有个，其中有4个能接收个能接收C31 B1频点数据频点数据Institute of Geodesy and Geophysics (IGG) Slide 9BeiDou-3卫星精密定轨方法与流程卫星精密定轨方法与流程(1/2)BeiDou-3卫星精密定轨卫星精密定轨-两步法间接联合两步法间接联合GPS单系统精密定轨单系统精密定轨 BeiDou-3 Satellites地面跟踪站目前数量和分布均有限； 将所选取的iGMAS 9个站及MGEX 52个站用于GPS单系统精密定轨；BDS单系统精密定轨单系统精密定轨 将得到的测站坐标等参数引入

11、，在北斗单系统定轨中固定或作为先验值； 将得到的测站坐标、测站天顶对流层延迟参数、地球自转参数输出。 观测数据：B1、B3消电离层双差相位及伪距数据；定轨弧段3天； 动力学参数：选用五参数ECOM光压模型以及伪随机脉冲参数。Institute of Geodesy and Geophysics (IGG) Slide 10BeiDou-3卫星精密定轨方法与流程卫星精密定轨方法与流程(2/2)BeiDou-3卫星精密定轨策略卫星精密定轨策略StrategyDescribeObservationsDouble-differenced ionosphere-free code and phase c

12、ombination of B1 and B3 in deployment phase 1Elevation angle cut-off7Sampling rate30 sData span3 daysStation coordinates and ERPsFixed using estimates from the GPS daily POD processWeighting4 mm for phase observations, 2 m for code observations and elevation-dependent weightingSatellite antenna PCO

13、and PCVPCO values are given by the Operation Control Department (apply for authorization), no PCV appliedSLR retro-reflector offsetsGiven by the Operation Control Department (apply for authorization)Attitude modelYaw-steering attitude modeGround antenna PCO and PCVNot appliedTroposphereDry and wet G

14、MF mapping functions; ZTDs estimated for each station at intervals of 2 hours; GPS daily ZTDs used as a priori valuesPrecession and nutationIAU 2010 precession and IAU 2010 nutation modelGeopotentialEGM2008 12*12Solid Earth tides, ocean tides and solid Earth pole tidesIERS Conventions 2010N bodyJPL

15、DE405 ephemerisOrbital parameters6 orbit elements and 5 ECOM SRP parameters were estimated:constants in the D, Y and X directions; periodic terms in the X directionPseudo-stochastic orbit parametersEvery 12 hours; constrained to 1E-6 m/sec in the radial direction, 1E-5 m/sec in the along-track direc

16、tion, and 1E-8 m/sec in the cross-track directionAmbiguityReal constant value for each ambiguity arcInstitute of Geodesy and Geophysics (IGG) Slide 11BeiDou-3卫星精密定轨初步结果展示卫星精密定轨初步结果展示 (1/13)定轨数据选取定轨数据选取 iGMAS跟踪站网 9个跟踪站2016.07.16-08.14观测数据；BeiDou-3轨道检核轨道检核 内符合（重叠弧段） 外符合（SLR检核） MGEX跟踪站网 52个跟踪站2016.07.1

17、6-08.14观测数据；BeiDou-2Institute of Geodesy and Geophysics (IGG) Slide 12BeiDou-3卫星精密定轨初步结果展示卫星精密定轨初步结果展示 (2/13)内符合内符合-重叠弧段互差重叠弧段互差Institute of Geodesy and Geophysics (IGG) Slide 13BeiDou-3卫星精密定轨初步结果展示卫星精密定轨初步结果展示 (3/13)外符合外符合-SLR检核检核本次处理期间，C31, C32, C33, C34, C01, C08, C10及C11卫星分别有12, 29, 45，10，114，36

18、，125及104个SLR观测值normal points (NPs)，这些观测值来自ILRS观测网络的Yarragadee, Monument, Papeete, Changchun, Shanghai, Mount Stro, Herstmonce and Matera测站 (见下图).Institute of Geodesy and Geophysics (IGG) Slide 14BeiDou-3卫星精密定轨初步结果展示卫星精密定轨初步结果展示 (4/13)内符合内符合-重叠弧段互差结果重叠弧段互差结果-BeiDou 3C31C32C33C34Institute of Geodesy a

19、nd Geophysics (IGG) Slide 15BeiDou-3卫星精密定轨初步结果展示卫星精密定轨初步结果展示 (5/13)内符合内符合-24h重叠弧段互差结果重叠弧段互差结果-BeiDou 3PRN of BeiDou SatellitesRadial (cm)Along(cm)Cross (cm)C3112.5124.9523.82C3211.8825.8124.32C3340.7558.6751.82C3441.1659.8859.11Institute of Geodesy and Geophysics (IGG) Slide 16BeiDou-3卫星精密定轨初步结果展示卫星

20、精密定轨初步结果展示 (6/13)外符合外符合-SLR结果结果-BeiDou 3C31C32C33C34Institute of Geodesy and Geophysics (IGG) Slide 17BeiDou-3卫星精密定轨初步结果展示卫星精密定轨初步结果展示 (7/13)外符合外符合-SLR结果结果-BeiDou 3PRN of BeiDou SatellitesMean (cm)RMS (cm)C31-3.8211.25C32-3.5111.28C3315.6438.61C3415.4440.26Institute of Geodesy and Geophysics (IGG) S

21、lide 18BeiDou-3卫星精密定轨初步结果展示卫星精密定轨初步结果展示 (10/13)内符合内符合-重叠弧段重叠弧段-BeiDou 2123456789101511121400.7PRN of BeiDou 2 SatellitesRMS/m RadialAlongCross1DInstitute of Geodesy and Geophysics (IGG) Slide 19BeiDou-3卫星精密定轨初步结果展示卫星精密定轨初步结果展示 (11/13)内符合内符合-重叠弧段重叠弧段-BeiDou 2PRN of BeiDou 2 Satellit

22、esRadial (cm)Along(cm)Cross (cm)C010.27390.58720.3649C020.25720.38790.2909C030.27320.30870.3678C040.41550.69190.4712C050.33580.41250.348C060.12350.20050.1077C070.11460.16120.1057C080.08490.2070.1503C090.11120.16670.1317C100.11230.20650.1465C150.09260.15630.122C110.05690.12450.0637C120.04460.13460.0731C140.05220.14410.1219Institute of Geodesy and Geophysics (IGG) Slide 20BeiDou-3卫星精密定轨初步结果展示卫星精密定轨初步结果展示 (12/13)C01C08C10C11外符合外符合-SLR结果结果-BeiDou 2Institute of Geodesy and Geoph