• 中国科学引文数据库(CSCD)
  • 中文科技期刊数据库
  • 中国核心期刊(遴选)数据库
  • 日本科学技术振兴机构数据库(JST)
  • 中国学术期刊(网络版)(CNKI)
  • 中国学术期刊综合评价数据库(CAJCED)
  • 中国超星期刊域出版平台

GPS/BDS-3系统时差稳定性分析

刘明玥, 涂锐, 洪菊, 苗亚格, 李芳馨

刘明玥, 涂锐, 洪菊, 苗亚格, 李芳馨. GPS/BDS-3系统时差稳定性分析[J]. 全球定位系统, 2022, 47(4): 55-63. DOI: 10.12265/j.gnss.2022055
引用本文: 刘明玥, 涂锐, 洪菊, 苗亚格, 李芳馨. GPS/BDS-3系统时差稳定性分析[J]. 全球定位系统, 2022, 47(4): 55-63. DOI: 10.12265/j.gnss.2022055
LIU Mingyue, TU Rui, HONG Ju, MIAO Yage, LI Fangxin. Stability analysis of GPS / BDS-3 system time offset[J]. GNSS World of China, 2022, 47(4): 55-63. DOI: 10.12265/j.gnss.2022055
Citation: LIU Mingyue, TU Rui, HONG Ju, MIAO Yage, LI Fangxin. Stability analysis of GPS / BDS-3 system time offset[J]. GNSS World of China, 2022, 47(4): 55-63. DOI: 10.12265/j.gnss.2022055

GPS/BDS-3系统时差稳定性分析

基金项目: 国家自然科学基金项目(41674034,41974032)
详细信息
    作者简介:

    刘明玥: (1995—),女,硕士,研究方向为卫星导航定位方法与技术

    涂锐: (1985—),男,博士,研究员,研究方向为GNSS精密定位、测速、时间传递、灾害监测等

    洪菊: (1994—),女,博士,研究方向GNSS精密单点定位增强技术

    通信作者:

    涂锐 E-mail:turui-2004@126.com

  • 中图分类号: P228.4

Stability analysis of GPS / BDS-3 system time offset

  • 摘要: 不同卫星导航系统之间时差的稳定性分析对时差的建模预报及应用具有重要价值, 因此基于由国际GNSS服务(IGS)中心提供的数据和产品从不同测站、不同观测量、不同时间长度三个方面对GPS和北斗三号全球卫星导航系统(BDS-3)的系统时差进行稳定性分析. 结果表明:不同测站因接收机、天线、时钟三种设备型号不同,提取的时差值存在差异,但稳定性相当;基于载波相位观测量解算的时差结果稳定性优于伪距观测量解算的时差结果,伪距得到的时差结果万秒频率稳定度在10−12~10−13量级,相位得到的时差结果万秒频率稳定度在10−14量级;通过单天与多天数据的时差解算结果比较,得到不同时长的时差稳定性一致.
    Abstract: The stability analysis of time offset between different satellite navigation systems is of great value in modeling prediction and application of time offset. Therefore, based on the data and products provided by the International GNSS Service Center (IGS), the stability analysis of system time offset between the GPS and BeiDou-3 Navigation Satellite System (BDS-3) from three aspects: different stations, different observations, and different time lengths. The results show: due to the different models of receivers, antennas and clocks, the time offset values extracted from different stations are different, but the stability is similar; The stability of the time offset result calculated based on the carrier phase observation is better than that of the pseudo range observation, the time offset result obtained by pseudo range has a frequency stability of the order of 10−12 to 10−13 in 10 000 seconds, and the time difference obtained by phase has a frequency stability of 10−14 seconds. By comparing the time offset calculation results of single-day and multi-day data, it is obtained that the time difference stability of different durations is consistent.
  • 随着全球卫星导航系统(GNSS)的不断发展和完善,GNSS所提供的定位、导航和授时(PNT)等服务已成为当今社会不可或缺的一部分[1]. 目前在轨运行的GNSS主要有美国的GPS、俄罗斯的GLONASS、欧洲的Galileo和中国的北斗卫星导航系统(BDS). 多GNSS兼容互操作是当前卫星导航的发展态势也是研究热点,而系统间的时差监测分析是实现GNSS兼容互操作的基础前提[2]. 系统间时差是由于不同的卫星导航系统都有各自的系统时,但不同导航系统时间的产生机制、时间尺度算法和溯源方式不相同,从而导致不同导航系统的系统时存在着几十纳秒的差异,该差异称为系统间时差[3].

    目前国内外学者对GNSS间时差监测的研究大多针对时差监测方法、时差监测精度评估以及时差应用等方面. 陈俊平等[4]分析了目前GNSS系统间时差监测主要的方法,并进行GPS/GLONASS系统时差解算,最后将GNSS的系统时差应用到多模导航定位,详细讨论了GPS/GLONASS时差以及测站硬件延迟. 于碧云[5]研究了Galileo系统时差监测与数据处理方法,同时验证GPS/Galileo系统时差在组合导航中的应用. 章洁君等[6]比较基于GNSS空间信号的三种时差监测方法,在距离不同的时间实验室之间进行时差试验,分析不同方法时差监测结果的差异. 同时,RICARDO等[7]分析了对GIOVE-A和GIOVE-B两颗实验卫星的组合定位精度以及时间偏差. MOUDRAK等[8]提出系统时差在组合导航定位中的应用的必要性,并对比利用广播星历进行组合导航定位的精度.

    时差监测结果对于多系统兼容互操作是必不可少的一个环节[9-10],同时,时差的稳定性在一定程度上影响着多系统融合的应用情况[11]. 在多系统GNSS时差结果实时预报应用中,时差稳定性是一个十分重要的影响因素[12-13]. 然而,当前在时差稳定性分析方面的相关研究较少.

    此外,北斗三号全球卫星导航系统(BDS-3)于2020年7月全面建成并开始提供全球服务,GPS与BDS-3的时差分析对实现GPS与BDS-3之间的兼容互操作,为全球用户提供更加优质的服务提供支撑作用. 因此,本文将对GPS与BDS-3系统时差稳定性进行分析.

    由于不同测站数据质量存在差异,导致某些天数据不能正常解算,但整体变化趋势分析不受个别结果的影响,因此论文选用13个国际GNSS服务(IGS)中心的测站从年积日第183—220天其中30天的数据和产品,分别采用伪距法和相位法获得钟差结果,再通过单站空间信号法得到系统时差监测结果,从设备型号、观测量、时间长度三个方面进行时差的稳定性分析.

    要分析GPS和BDS-3系统间时差稳定性,首先要得到两系统间的时差结果. 本文采用单站空间信号法以获得GPS和BDS-3的时差值.

    单站空间信号法是指单一测站的多模接收机通过接收卫星导航系统信号,对该信号经一系列处理后得到原始观测数据,然后再解算得到卫星导航系统时间与接收机本地参考时间之间的钟差,利用不同导航系统数据得到的钟差值作差即可得到系统时差[14-15]. 空间信号法的原理如图1所示,解算方式如式(1)所示. 本文分别解算GPS和BDS-3单系统时间与接收机本地参考时间之间的钟差,再将钟差值作差后得到GPS与BDS-3两个导航系统之间的时差结果.

    图  1  单站空间信号法原理图
    $$ \begin{split} {T}_{\text{GPS}}-{T}_{\text{BDS-3}}&= ({T}_{\mathrm{Re}\text{f}}-{T}_{\text{BDS-3}})-({T}_{\mathrm{Re}\text{f}}-{T}_{\text{GPS}})\\&= {\delta }_{{T}}{}_{{}_{\text{R,BDS-3}}}-{\delta }_{{T}}{}_{{}_{\text{R,GPS}}}. \end{split}$$ (1)

    式中:$ {T_{{\text{GPS}}}} $为GPS系统时间;${T_{{\text{BDS-3}}}}$为BDS-3时间;$ {T_{{\text{Ref}}}} $为接收机本地参考时间;${\delta _{{T}}}_{_{{\text{R,BDS-3}}}}$为接收机本地参考时间与BDS-3时间之间的时差值;${\delta _{{T}}}_{_{{\text{R,GPS}}}}$为接收机本地参考时间与GPS时间之间的时差值[16].

    为得到接收机本地参考时间与GNSS系统时间之间的钟差值,本文分别采用伪距法和相位法进行解算[17].

    伪距法以测距码为观测量进行测距. 该方法是卫星导航定位中最基本的方法,具有解算速度快、解算结果唯一等优势,但测量精度有限[18]. 观测方程为

    $$ \widetilde \rho = \rho _{\text{R}}^{\text{S}} + c \cdot ({\delta _{{{{t}}_{\text{R}}}}} - {\delta _{{{{t}}^{\text{S}}}}}) - {\tau _{{\text{ion}}}} + {\tau _{{\text{trop}}}} + \varepsilon . $$ (2)

    式中:$\;\widetilde\rho$为伪距观测值;$ \;\rho _{\text{R}}^{\text{S}} $为卫星S至接收机R的几何距离;c为光速;${\delta _{{{{T}}_{\text{R}}}}}$为接收机钟相对于标准卫星导航系统时间的偏差,${\delta _{{{{T}}^{\text{S}}}}}$为卫星钟相对于标准卫星导航系统时间的偏差;$ {\tau _{{\text{ion}}}} $为电离层延迟;$ {\tau _{{\text{trop}}}} $为对流层延迟;$ \varepsilon $为观测误差.

    相位法以载波相位为观测量进行测距,载波相位观测量是卫星导航系统接收机接收的卫星载波信号与接收机本振信号的相位差,再乘以波长,可以得到卫星至接收机的粗略距离[6],观测方程为

    $$ \lambda \varphi = \rho _{\rm{R}}^{\rm{S}} + c \cdot ({\delta _{{{{t}}_{\rm{R}}}}} - {\delta _{{{{t}}^{\rm{S}}}}}) + \lambda N_{\rm{R}}^{\rm{S}} - {\tau _{{\rm{ion}}}} + {\tau _{{\rm{trop}}}} + \varepsilon . $$ (3)

    式中:$ \lambda $为载波观测量的波长;$ \varphi $表示相位差;$ N_{\text{R}}^{\text{S}} $表示整周模糊度.

    本文采用IGS提供的精密产品和观测数据,利用双频无电离层组合观测消除电离层误差,将对流层干分量延迟误差和部分湿分量延迟误差、相对论效应、地球自转误差、天线相位缠绕误差、地球潮汐误差等进行模型改正,卫星端和接收机端的天线相位中心使用IGS提供的igs14.atx文件进行改正[19],对流层残差采用随机游走过程估计,并强约束测站坐标,通过标准单点定位求解得到伪距法的时差结果,通过精密单点定位(PPP)求解得到相位法的时差结果[18]. 同时,采用正反二次解算来减小载波相位模糊度收敛对结果的影响.

    每个GNSS测站所采用的设备存在差异,为研究该差异是否对系统时差稳定性产生影响,选取13个IGS测站年积日从第183—220天其中30天的数据进行分析,数据采样间隔为30 s,软件采用自研的GNSS精密时差分析软件. 根据设备型号将所选测站分为四组,分别为型号完全相同组、接收机型号不同组、天线型号不同组、时钟型号不同组,测站和设备详细信息如表1所示.

    表  1  测站信息表
    参数测站名接收机型号时钟类型天线型号
    型号
    相同
    METGSEPT POLARX5INTERNALTRM59800.00
    NKLGSEPT POLARX5INTERNALTRM59800.00
    PTGGSEPT POLARX5INTERNALTRM59800.00
    SEYGSEPT POLARX5INTERNALTRM59800.00
    接收机不同NKLGSEPT POLARX5INTERNALTRM59800.00
    PTGGSEPT POLARX5INTERNALTRM59800.00
    ASCGTRIMBLE ALLOYINTERNALTRM59800.00
    CHPGTRIMBLE ALLOYINTERNALTRM59800.00
    天线
    不同
    ARUCSEPT POLARX5INTERNALASH701945C_M
    FAIRSEPT POLARX5INTERNALASH701945G_M
    DGARSEPT POLARX5INTERNALASH701945E_M
    MDO1SEPT POLARX5INTERNALJAVRINGANT_DM
    PTGGSEPT POLARX5INTERNALTRM59800.00
    时钟
    不同
    GAMGSEPT POLARX5TRINTERNALLEIAR25.R4
    MGUESEPT POLARX5TREXTERNAL H-MASERLEIAR25.R4
    YEL2SEPT POLARX5TREXTERNAL VCH-1008 MASERLEIAR25.R4
    下载: 导出CSV 
    | 显示表格

    对四组测站分别采用伪距法和相位法两种方式解算GPS/BDS-3系统间时差,图2~3分别为各测站第183天伪距法和相位法解算的时差结果,图4~5分别为各测站30天伪距法和相位法解算的时差结果. 由图可知,不管是伪距法还是相位法,不同的设备型号对时差结果稳定性的影响较小. 然而,由图4~5可知,第23天的时差结果发生一个跳跃性波动,发生该变化的测站采用接收机型号均为SEPT POLARX5型接收机,该现象可能与接收机固件升级有着一定的关系,这对长期时差稳定性的监测带来一定影响. 虽然不同型号设备对时差稳定性影响较小,但因不同型号设备时延的差异导致求解的系统间时差值有着较大偏差,由图2(d)图4(d)可知,伪距法的时差结果相差近40 ns;由图3(d)图5(d)可知,相位法的时差结果相差近百纳秒.

    图  2  基于伪距法解算的第183天时差结果
    图  3  基于相位法解算的第183天时差结果
    图  4  基于伪距法解算的30天时差结果
    图  5  基于相位法解算的30天时差结果

    为研究伪距和相位观测值解算的钟差结果对时差稳定性的影响,对所选测站分别使用伪距法和相位法进行单天和多天两种情况下系统时差稳定性分析,由图2~5可知,使用载波相位观测量所解算的时差结果稳定性优于使用伪距测量得到的时差结果. 单天时差监测结果中,使用伪距观测量解算的时差结果呈现较大噪声的波动,波动幅度可达25 ns,而使用载波相位观测量得到的结果中,时差值在图中近似一条直线,波动较小. 多天监测结果有着同样的结论,采用载波相位观测量解算的时差结果在30天内的稳定性明显优于采用伪距测量得到的时差结果. 其原因是在GNSS信号观测中,伪距测量因码元宽度较大测量精度不高所致,而载波相位观测量精度通常是伪距的100倍.

    为进一步分析不同观测量对时差结果的影响,对第183天各测站时差结果采用修正阿伦方差计算频率稳定度,并对其进行分析. 图6为伪距法得到的修正阿伦方差结果,图7为相位法得到的修正阿伦方差结果. 由图6~7可知,伪距法得到的时差结果万秒频率稳定度在10−12~10−13量级,而采用相位法得到的时差结果万秒频率稳定度在10−14量级.

    图  6  伪距法第183天修正阿伦方差结果
    图  7  相位法第183天修正阿伦方差结果

    为研究不同时间长度的时差结果的稳定性,将单天结果和30天结果进行对比分析. 从表2~3的时差标准差可得,采用相位法解算的单天时差标准差为亚纳秒级,30天时差标准差约在5 ns;伪距法得到的单天和30天时差标准差相差不大,约在5~10 ns. 结合图5可知,天与天之间相连的地方存在不同大小的时差值波动. 这种波动是由于采用事后单天独立解算策略而非实时连续多天解算使得相邻两天的整周模糊度求解值存在差异导致的,通过对时差结果进行归算处理后,其30天时差值变化曲线与单天相似,近似一条直线,同时归算后时差标准差与单天时差标准差接近,为亚纳秒级.

    表  2  各测站各天时差结果标准差统计
    测站ARUCASCGCHPGDGARFAIRGAMGMDO1METGMGUENKLGPTGGSEYGYEL2
    解算
    方法
    PPPSPPPPPSPPPPPSPPPPPSPPPPPSPPPPPSPPPPPSPPPPPSPPPPPSPPPPPSPPPPPSPPPPPSPPPPPSPP
    1 0.30 4.45 0.31 4.22 0.38 5.84 0.27 2.47 0.31 6.11 0.50 4.27 0.31 5.06 0.35 3.49 0.50 5.71 0.18 4.82 0.39 3.98 0.21 3.34 0.21 4.29
    2 0.40 3.98 0.24 5.53 0.22 4.21 0.14 3.11 0.39 5.29 0.31 3.15 0.52 5.06 0.35 3.77 0.23 4.19 0.42 4.10 0.31 3.31 0.25 3.39 0.68 4.58
    3 0.46 5.43 0.48 5.25 0.35 6.41 0.13 2.90 0.51 5.17 0.24 3.92 0.59 5.56 0.37 4.90 0.32 4.70 0.34 5.34 0.16 3.44 0.59 3.41 0.49 5.24
    4 0.52 5.65 0.48 5.91 0.51 5.78 0.19 3.10 1.63 6.48 0.23 3.33 0.56 5.57 0.52 5.65 0.45 5.66 0.25 4.90 0.18 3.68 0.19 3.66 0.62 5.57
    5 0.38 4.96 0.45 5.23 0.36 4.42 0.15 3.43 0.54 6.62 0.37 3.17 0.33 6.57 0.40 5.12 0.45 5.98 0.26 4.67 0.21 3.25 0.19 3.69 1.20 4.73
    6 0.44 5.27 0.32 6.00 0.27 6.98 0.19 3.62 0.36 5.81 0.46 3.74 0.18 5.89 0.39 4.51 0.41 5.21 0.18 4.80 0.31 3.52 0.13 3.75 0.85 4.45
    7 0.31 4.68 0.32 4.10 0.38 5.89 0.26 2.43 0.28 5.77 0.48 4.22 0.34 5.19 0.35 4.20 0.38 5.86 0.17 4.69 0.36 3.82 0.23 3.36 0.35 4.28
    8 0.42 3.93 0.19 5.47 0.24 4.26 0.25 3.12 0.37 5.42 0.36 3.07 0.46 4.79 0.33 3.68 0.23 4.16 0.34 3.98 0.27 3.13 0.18 3.40 0.40 4.47
    9 0.51 5.35 0.42 5.22 0.49 6.30 0.12 2.94 0.52 5.29 0.25 3.95 0.55 6.02 0.40 4.99 0.36 4.63 0.38 5.27 0.48 3.46 0.19 3.34 1.30 5.34
    10 0.40 5.63 0.50 5.90 0.47 5.67 0.19 3.06 0.56 6.34 0.26 3.30 0.54 5.37 0.54 5.63 0.44 5.63 0.31 4.83 0.16 3.45 0.21 3.71 0.51 5.67
    11 0.45 4.79 0.26 5.20 0.39 4.29 0.16 3.08 0.51 6.45 0.42 3.35 0.35 6.49 0.47 4.96 0.53 6.14 0.27 4.54 0.21 3.26 0.22 3.71 0.44 4.69
    12 0.35 5.16 0.24 5.98 0.28 7.12 0.16 3.57 0.41 5.99 0.44 3.77 0.12 5.92 0.45 4.58 0.35 5.14 0.19 4.85 0.25 3.38 0.17 3.59 1.29 4.45
    13 0.35 4.48 0.32 4.46 0.29 5.91 0.39 2.36 0.33 5.78 0.63 4.87 0.29 5.31 0.36 4.22 0.29 5.72 0.18 4.60 0.51 4.32 0.33 3.89 0.18 4.27
    14 0.36 3.89 0.20 5.49 0.26 4.47 2.79 3.35 0.40 5.37 0.36 3.21 0.48 4.91 0.35 3.64 0.23 4.06 0.32 4.27 0.31 3.14 0.17 3.70 0.42 4.51
    15 0.55 5.62 0.44 5.20 0.48 6.47 0.14 4.03 0.50 5.33 0.24 4.59 0.53 5.67 0.40 5.32 0.37 4.62 0.35 5.97 0.20 3.83 0.15 4.13 0.58 5.30
    16 0.54 5.37 0.43 6.09 0.50 5.64 0.30 3.19 0.54 6.28 0.32 3.35 0.53 5.46 0.43 5.41 0.48 5.67 0.25 4.91 0.15 3.45 0.26 3.98 0.92 5.41
    17 0.57 4.71 0.34 5.07 0.37 4.25 0.40 3.18 0.58 6.53 0.72 3.45 0.38 6.53 0.47 4.84 0.48 6.27 0.36 4.42 0.75 3.07 0.47 3.73 0.33 4.76
    18 0.31 5.21 0.25 6.11 0.25 7.24 0.20 3.63 1.61 6.04 0.46 3.76 0.11 6.60 0.43 4.66 0.41 4.90 0.23 4.85 0.29 3.41 0.16 3.61 0.47 4.42
    19 0.26 4.33 0.31 4.62 0.19 6.13 0.25 2.28 0.33 5.80 0.50 4.26 0.30 5.51 0.38 4.60 0.25 5.71 0.19 4.42 0.37 4.14 0.23 3.48 0.85 4.14
    20 0.35 3.81 0.23 5.65 0.27 4.41 1.66 3.18 0.43 5.53 0.37 3.19 0.53 4.90 0.40 3.70 0.24 3.94 0.33 4.12 0.28 3.12 0.19 3.56 0.49 4.83
    21 0.46 5.38 0.44 5.10 0.37 6.50 1.26 2.86 0.49 5.49 0.27 4.06 0.56 5.48 0.42 4.95 0.36 4.70 0.34 5.33 0.17 3.46 0.10 3.38 0.48 5.02
    22 0.44 4.61 0.30 4.76 0.44 4.30 0.20 3.18 0.52 6.54 0.32 3.47 0.38 7.38 0.40 4.78 0.56 6.36 0.31 4.31 0.20 3.15 0.61 3.75 1.37 4.61
    23 0.80 5.17 0.24 6.02 0.24 7.17 0.32 3.68 0.57 5.79 0.69 3.78 0.15 6.61 0.77 4.66 0.39 4.84 0.27 4.94 0.68 4.14 0.35 3.80 0.29 4.50
    24 0.68 71.90 0.34 21 119.66 0.38 11.07 1.11 13.17 0.45 18.89 1.51 90 064.09 0.34 33.23 0.69 10.95 0.48 11 176.48 0.27 9.20 3.27 41.64 0.52 9.77 0.46 8.75
    25 0.44 3.79 0.24 5.71 0.51 4.41 0.15 3.19 0.40 5.82 0.41 3.23 0.52 5.14 0.38 3.91 0.25 3.94 0.34 4.35 0.59 3.16 0.16 3.53 0.39 5.18
    26 0.53 5.80 0.44 5.09 0.34 6.55 0.16 3.12 0.46 5.37 0.28 4.36 0.59 5.80 0.43 5.08 0.37 4.65 0.33 5.33 0.18 3.72 0.16 3.92 0.46 4.80
    27 0.56 5.63 0.47 6.28 0.57 5.61 0.21 3.16 0.57 5.93 0.25 3.67 0.55 5.52 0.46 5.61 0.50 5.53 0.27 4.91 0.07 3.93 0.22 4.20 1.18 5.89
    28 0.52 5.56 0.27 4.62 0.43 4.32 0.22 3.38 0.62 6.49 0.38 3.49 0.45 6.32 0.48 5.08 0.59 6.34 0.28 4.17 0.24 2.71 0.26 4.19 0.42 4.64
    29 0.50 5.50 0.26 5.87 0.71 7.17 0.20 3.60 0.48 5.78 0.47 3.79 0.15 6.05 0.54 4.64 0.41 4.86 0.29 4.91 0.49 3.62 0.20 3.81 0.22 4.38
    30 0.61 5.62 0.43 6.02 0.52 5.67 0.18 4.23 0.50 5.92 0.27 3.67 0.52 5.47 0.45 5.97 0.52 5.52 0.29 5.50 0.10 3.50 0.22 5.19 2.26 5.73
    下载: 导出CSV 
    | 显示表格
    表  3  各测站30天时差结果标准差统计
    测站 ARUCASCGCHPGDGARFAIRGAMGMDO1METGMGUENKLGPTGGSEYGYEL2
    30天相位法标准差4.7784.7484.6114.6734.8004.6945.1114.3244.6494.6884.6104.7524.755
    归算后30天相位法标准差0.4780.5540.5860.7040.6280.4900.4400.4530.4120.2990.6890.2840.809
    30天伪距法标准差14.0983856.0346.1454.4856.89816440.6808.5165.2562040.5685.1998.4614.6805.284
    下载: 导出CSV 
    | 显示表格

    图8~9分别为伪距法和相位法得到的多天时差结果进行相邻两天的时差一阶差分的结果,图10是归算后相位法相邻两天的时差一阶差分结果,图11~12分别为伪距法和相位法单天时差结果进行相邻历元一阶差分结果. 从图可以得出,伪距法得到的差分结果中,单天与30天的差分波动范围相近,在±4 ns内波动. 相位法未进行归算处理的差值结果较为离散,但经归算后得到的结果与单天一致,都在零值附近.

    图  8  伪距法相邻两天的时差一阶差分结果
    图  9  相位法相邻两天的时差一阶差分结果
    图  10  归算后相位法相邻两天的时差一阶差分结果
    图  11  伪距法单天相邻历元时差一阶差分结果
    图  12  相位法单天相邻历元时差一阶差分结果

    本文使用IGS提供的多测站多天的数据和产品,分别采用伪距法和相位法两种方法解算钟差结果,再通过单站空间信号法得到GPS/BDS-3系统间时差. 分别从不同设备型号、不同观测量、不同时间长度三个方面对时差稳定性进行了分析,得到如下结论:

    1)不同测站在接收机、天线、钟三种设备有所差异的情况下GPS/BDS-3时差稳定性相当,即不同设备型号对其稳定性影响较小. 但不同设备型号时延的差异会造成时差值的不同.

    2)通过比较伪距法和相位法对时差稳定性的影响,结果表明,采用载波相位观测量解算的时差结果稳定性优于伪距测量方法,伪距法得到的时差结果万秒频率稳定度在10−12~10−13,相位法得到的时差结果万秒频率稳定度约在10−14.

    3)对比单天与30天时差解算结果,结果表明,除相位法因数据处理策略导致的天与天之间的跳跃外,不同时间长度监测的时差稳定性基本相同,伪距法约为5~10 ns,相位法为亚纳秒级.

    本文研究的GPS和BDS-3系统间时差的稳定性,不仅为GNSS时间兼容互操作提供了技术支撑,同时为下一步的时差建模预报以及应用奠定了基础.

    致谢:本文的研究得到了国家自然科学基金项目(41674034,41974032)的支持.

  • 图  1   单站空间信号法原理图

    图  2   基于伪距法解算的第183天时差结果

    图  3   基于相位法解算的第183天时差结果

    图  4   基于伪距法解算的30天时差结果

    图  5   基于相位法解算的30天时差结果

    图  6   伪距法第183天修正阿伦方差结果

    图  7   相位法第183天修正阿伦方差结果

    图  8   伪距法相邻两天的时差一阶差分结果

    图  9   相位法相邻两天的时差一阶差分结果

    图  10   归算后相位法相邻两天的时差一阶差分结果

    图  11   伪距法单天相邻历元时差一阶差分结果

    图  12   相位法单天相邻历元时差一阶差分结果

    表  1   测站信息表

    参数测站名接收机型号时钟类型天线型号
    型号
    相同
    METGSEPT POLARX5INTERNALTRM59800.00
    NKLGSEPT POLARX5INTERNALTRM59800.00
    PTGGSEPT POLARX5INTERNALTRM59800.00
    SEYGSEPT POLARX5INTERNALTRM59800.00
    接收机不同NKLGSEPT POLARX5INTERNALTRM59800.00
    PTGGSEPT POLARX5INTERNALTRM59800.00
    ASCGTRIMBLE ALLOYINTERNALTRM59800.00
    CHPGTRIMBLE ALLOYINTERNALTRM59800.00
    天线
    不同
    ARUCSEPT POLARX5INTERNALASH701945C_M
    FAIRSEPT POLARX5INTERNALASH701945G_M
    DGARSEPT POLARX5INTERNALASH701945E_M
    MDO1SEPT POLARX5INTERNALJAVRINGANT_DM
    PTGGSEPT POLARX5INTERNALTRM59800.00
    时钟
    不同
    GAMGSEPT POLARX5TRINTERNALLEIAR25.R4
    MGUESEPT POLARX5TREXTERNAL H-MASERLEIAR25.R4
    YEL2SEPT POLARX5TREXTERNAL VCH-1008 MASERLEIAR25.R4
    下载: 导出CSV

    表  2   各测站各天时差结果标准差统计

    测站ARUCASCGCHPGDGARFAIRGAMGMDO1METGMGUENKLGPTGGSEYGYEL2
    解算
    方法
    PPPSPPPPPSPPPPPSPPPPPSPPPPPSPPPPPSPPPPPSPPPPPSPPPPPSPPPPPSPPPPPSPPPPPSPPPPPSPP
    1 0.30 4.45 0.31 4.22 0.38 5.84 0.27 2.47 0.31 6.11 0.50 4.27 0.31 5.06 0.35 3.49 0.50 5.71 0.18 4.82 0.39 3.98 0.21 3.34 0.21 4.29
    2 0.40 3.98 0.24 5.53 0.22 4.21 0.14 3.11 0.39 5.29 0.31 3.15 0.52 5.06 0.35 3.77 0.23 4.19 0.42 4.10 0.31 3.31 0.25 3.39 0.68 4.58
    3 0.46 5.43 0.48 5.25 0.35 6.41 0.13 2.90 0.51 5.17 0.24 3.92 0.59 5.56 0.37 4.90 0.32 4.70 0.34 5.34 0.16 3.44 0.59 3.41 0.49 5.24
    4 0.52 5.65 0.48 5.91 0.51 5.78 0.19 3.10 1.63 6.48 0.23 3.33 0.56 5.57 0.52 5.65 0.45 5.66 0.25 4.90 0.18 3.68 0.19 3.66 0.62 5.57
    5 0.38 4.96 0.45 5.23 0.36 4.42 0.15 3.43 0.54 6.62 0.37 3.17 0.33 6.57 0.40 5.12 0.45 5.98 0.26 4.67 0.21 3.25 0.19 3.69 1.20 4.73
    6 0.44 5.27 0.32 6.00 0.27 6.98 0.19 3.62 0.36 5.81 0.46 3.74 0.18 5.89 0.39 4.51 0.41 5.21 0.18 4.80 0.31 3.52 0.13 3.75 0.85 4.45
    7 0.31 4.68 0.32 4.10 0.38 5.89 0.26 2.43 0.28 5.77 0.48 4.22 0.34 5.19 0.35 4.20 0.38 5.86 0.17 4.69 0.36 3.82 0.23 3.36 0.35 4.28
    8 0.42 3.93 0.19 5.47 0.24 4.26 0.25 3.12 0.37 5.42 0.36 3.07 0.46 4.79 0.33 3.68 0.23 4.16 0.34 3.98 0.27 3.13 0.18 3.40 0.40 4.47
    9 0.51 5.35 0.42 5.22 0.49 6.30 0.12 2.94 0.52 5.29 0.25 3.95 0.55 6.02 0.40 4.99 0.36 4.63 0.38 5.27 0.48 3.46 0.19 3.34 1.30 5.34
    10 0.40 5.63 0.50 5.90 0.47 5.67 0.19 3.06 0.56 6.34 0.26 3.30 0.54 5.37 0.54 5.63 0.44 5.63 0.31 4.83 0.16 3.45 0.21 3.71 0.51 5.67
    11 0.45 4.79 0.26 5.20 0.39 4.29 0.16 3.08 0.51 6.45 0.42 3.35 0.35 6.49 0.47 4.96 0.53 6.14 0.27 4.54 0.21 3.26 0.22 3.71 0.44 4.69
    12 0.35 5.16 0.24 5.98 0.28 7.12 0.16 3.57 0.41 5.99 0.44 3.77 0.12 5.92 0.45 4.58 0.35 5.14 0.19 4.85 0.25 3.38 0.17 3.59 1.29 4.45
    13 0.35 4.48 0.32 4.46 0.29 5.91 0.39 2.36 0.33 5.78 0.63 4.87 0.29 5.31 0.36 4.22 0.29 5.72 0.18 4.60 0.51 4.32 0.33 3.89 0.18 4.27
    14 0.36 3.89 0.20 5.49 0.26 4.47 2.79 3.35 0.40 5.37 0.36 3.21 0.48 4.91 0.35 3.64 0.23 4.06 0.32 4.27 0.31 3.14 0.17 3.70 0.42 4.51
    15 0.55 5.62 0.44 5.20 0.48 6.47 0.14 4.03 0.50 5.33 0.24 4.59 0.53 5.67 0.40 5.32 0.37 4.62 0.35 5.97 0.20 3.83 0.15 4.13 0.58 5.30
    16 0.54 5.37 0.43 6.09 0.50 5.64 0.30 3.19 0.54 6.28 0.32 3.35 0.53 5.46 0.43 5.41 0.48 5.67 0.25 4.91 0.15 3.45 0.26 3.98 0.92 5.41
    17 0.57 4.71 0.34 5.07 0.37 4.25 0.40 3.18 0.58 6.53 0.72 3.45 0.38 6.53 0.47 4.84 0.48 6.27 0.36 4.42 0.75 3.07 0.47 3.73 0.33 4.76
    18 0.31 5.21 0.25 6.11 0.25 7.24 0.20 3.63 1.61 6.04 0.46 3.76 0.11 6.60 0.43 4.66 0.41 4.90 0.23 4.85 0.29 3.41 0.16 3.61 0.47 4.42
    19 0.26 4.33 0.31 4.62 0.19 6.13 0.25 2.28 0.33 5.80 0.50 4.26 0.30 5.51 0.38 4.60 0.25 5.71 0.19 4.42 0.37 4.14 0.23 3.48 0.85 4.14
    20 0.35 3.81 0.23 5.65 0.27 4.41 1.66 3.18 0.43 5.53 0.37 3.19 0.53 4.90 0.40 3.70 0.24 3.94 0.33 4.12 0.28 3.12 0.19 3.56 0.49 4.83
    21 0.46 5.38 0.44 5.10 0.37 6.50 1.26 2.86 0.49 5.49 0.27 4.06 0.56 5.48 0.42 4.95 0.36 4.70 0.34 5.33 0.17 3.46 0.10 3.38 0.48 5.02
    22 0.44 4.61 0.30 4.76 0.44 4.30 0.20 3.18 0.52 6.54 0.32 3.47 0.38 7.38 0.40 4.78 0.56 6.36 0.31 4.31 0.20 3.15 0.61 3.75 1.37 4.61
    23 0.80 5.17 0.24 6.02 0.24 7.17 0.32 3.68 0.57 5.79 0.69 3.78 0.15 6.61 0.77 4.66 0.39 4.84 0.27 4.94 0.68 4.14 0.35 3.80 0.29 4.50
    24 0.68 71.90 0.34 21 119.66 0.38 11.07 1.11 13.17 0.45 18.89 1.51 90 064.09 0.34 33.23 0.69 10.95 0.48 11 176.48 0.27 9.20 3.27 41.64 0.52 9.77 0.46 8.75
    25 0.44 3.79 0.24 5.71 0.51 4.41 0.15 3.19 0.40 5.82 0.41 3.23 0.52 5.14 0.38 3.91 0.25 3.94 0.34 4.35 0.59 3.16 0.16 3.53 0.39 5.18
    26 0.53 5.80 0.44 5.09 0.34 6.55 0.16 3.12 0.46 5.37 0.28 4.36 0.59 5.80 0.43 5.08 0.37 4.65 0.33 5.33 0.18 3.72 0.16 3.92 0.46 4.80
    27 0.56 5.63 0.47 6.28 0.57 5.61 0.21 3.16 0.57 5.93 0.25 3.67 0.55 5.52 0.46 5.61 0.50 5.53 0.27 4.91 0.07 3.93 0.22 4.20 1.18 5.89
    28 0.52 5.56 0.27 4.62 0.43 4.32 0.22 3.38 0.62 6.49 0.38 3.49 0.45 6.32 0.48 5.08 0.59 6.34 0.28 4.17 0.24 2.71 0.26 4.19 0.42 4.64
    29 0.50 5.50 0.26 5.87 0.71 7.17 0.20 3.60 0.48 5.78 0.47 3.79 0.15 6.05 0.54 4.64 0.41 4.86 0.29 4.91 0.49 3.62 0.20 3.81 0.22 4.38
    30 0.61 5.62 0.43 6.02 0.52 5.67 0.18 4.23 0.50 5.92 0.27 3.67 0.52 5.47 0.45 5.97 0.52 5.52 0.29 5.50 0.10 3.50 0.22 5.19 2.26 5.73
    下载: 导出CSV

    表  3   各测站30天时差结果标准差统计

    测站 ARUCASCGCHPGDGARFAIRGAMGMDO1METGMGUENKLGPTGGSEYGYEL2
    30天相位法标准差4.7784.7484.6114.6734.8004.6945.1114.3244.6494.6884.6104.7524.755
    归算后30天相位法标准差0.4780.5540.5860.7040.6280.4900.4400.4530.4120.2990.6890.2840.809
    30天伪距法标准差14.0983856.0346.1454.4856.89816440.6808.5165.2562040.5685.1998.4614.6805.284
    下载: 导出CSV
  • [1] 袁海波, 张继海, 广伟. 不同GNSS时差监测方法研究[C]//第九届中国卫星导航学术年会, 2018.
    [2] 朱琳, 张慧君, 李孝辉, 等. 全球卫星导航系统时差监测与修正方法研究[C]//第四届中国卫星导航学术年会, 2013.
    [3] 尹海博, 迅郭杭, 罗孝文. BDS/GNSS融合精密单点定位性能分析[J]. 全球定位系统, 2021, 46(3): 66-71. DOI: 10.12265/j.gnss.2020121501
    [4] 陈俊平, 吴斌, 胡小工, 等. GPS/GLONSS时差监测及其在多模定位中的应用[C]//第三届中国卫星导航学术年会, 2012.
    [5] 于碧云. Galileo系统时差监测方法研究与实现[D]. 西安: 中国科学大学, 2016.
    [6] 章洁君, 陈俊平, 张益泽. GNSS空间信号法时差监测方法与结果分析[J]. 测绘通报, 2017(1): 26-29,92. DOI: 10.13474/j.cnki.11-2246.2017.0006
    [7]

    PÍRIZ R, GARCÍA Á M, TOBÍAS G, et al. GNSS interoperability: offset between reference time scales and timing biases[J]. Metrologia, 2008, 45(6): S87. DOI: 10.1088/0026-1394/45/6/s14

    [8]

    MOUDRAK A, KONOVALTSEV A, FURTHNER J, et al. Timing aspects of GPS Galileo interoperability: challenges and solutions[C]// The 36th Annual Precise Time and Time Interval (PTTI) Systems and Applications Meeting, 2004: 279-292. DOI: 10.1008/0026-1394/45/6/s14

    [9] 张慧君, 李孝辉. GNSS系统时差监测与预报研究进展[C]//2013中国天文学会学术年会, 2013.
    [10] 广伟. GNSS时间互操作关键技术研究[D]. 西安: 中国科学院大学(中国科学院国家授时中心), 2019.
    [11]

    HAHN J H , POWERS E D . Implementation of the GPS to Galileo Time Offset (GGTO)[C]// Proceeding of the 2005 International Frequency Control Symposium and Exposition, 2005. DOI: 10.1109/FREQ.2005.1573899

    [12]

    HEO Y J, CHO J, HEO M B. Improving prediction accuracy of GPS satellite clocks with periodic variation behaviour[J]. Measurement science and technology, 2010, 21(7): 073001. DOI: 10.1088/0957-0233/21/7/073001

    [13]

    HUANG G W, ZHANG Q, XU G C. Real-time clock offset prediction with an improved model[J]. GPS solutions, 2014, 18(1): 95-104. DOI: 10.1007/s10291-013-0313-0

    [14] 李红涛. 基于GPS和GLONASS的单站授时和时差监测研究[D]. 西安: 长安大学, 2012.
    [15] 孙广, 陆华, 沈利荣, 等. 基于不同GNSS系统的单站时差监测评估[C]//新型导航技术及应用研讨会, 2015.
    [16] 孔思嘉, 谢维华. 基于GLONASS频间偏差的GNSS时差监测方法研究[J]. 全球定位系统, 2020, 45(4): 21-28. DOI: 10.13442/j.gnss.1008-9268.2020.04.004
    [17]

    HUANG G W, ZHANG Q, FU W J, et al. GPS/GLONASS time offset monitoring based on combined precise point positioning (PPP) approach[J]. Advances in space research, 2015, 55(12): 2950-2960. DOI: 10.1016/j.asr.2015.03.003

    [18] 杨松. 多模多频GNSS环境下精密单点定位技术研究[D]. 贵阳: 贵州大学, 2017.
    [19] 洪菊, 涂锐, 王星星, 等. 系统偏差常数弧段估计策略对多系统静态精密单点定位的影响[J]. 大地测量与地球动力学, 2019, 39(9): 947-951. DOI: 10.14075/j.jgg.2019.09.014
  • 期刊类型引用(1)

    1. 戴群雄,解晋,丁格,左兆辉,尹继凯. 一种变温环境下高精度铷钟温度特性自动标定及补偿方法. 无线电工程. 2025(05): 1021-1029 . 百度学术

    其他类型引用(0)

图(12)  /  表(3)
计量
  • 文章访问数:  750
  • HTML全文浏览量:  119
  • PDF下载量:  59
  • 被引次数: 1
出版历程
  • 收稿日期:  2022-04-05
  • 录用日期:  2022-07-04
  • 网络出版日期:  2022-07-27
  • 刊出日期:  2022-09-04

目录

/

返回文章
返回
x 关闭 永久关闭