留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

SBAS星历改正数及UDRE参数生成算法分析

金彪 魏巍 陈姗姗 李东俊

金彪, 魏巍, 陈姗姗, 李东俊. SBAS星历改正数及UDRE参数生成算法分析[J]. 武汉大学学报 ( 信息科学版), 2021, 46(1): 111-117. doi: 10.13203/j.whugis20190033
引用本文: 金彪, 魏巍, 陈姗姗, 李东俊. SBAS星历改正数及UDRE参数生成算法分析[J]. 武汉大学学报 ( 信息科学版), 2021, 46(1): 111-117. doi: 10.13203/j.whugis20190033
JIN Biao, WEI Wei, CHEN Shanshan, LI Dongjun. Analysis of SBAS Ephemeris Correction and UDRE Generation Algorithm[J]. Geomatics and Information Science of Wuhan University, 2021, 46(1): 111-117. doi: 10.13203/j.whugis20190033
Citation: JIN Biao, WEI Wei, CHEN Shanshan, LI Dongjun. Analysis of SBAS Ephemeris Correction and UDRE Generation Algorithm[J]. Geomatics and Information Science of Wuhan University, 2021, 46(1): 111-117. doi: 10.13203/j.whugis20190033

SBAS星历改正数及UDRE参数生成算法分析

doi: 10.13203/j.whugis20190033
基金项目: 

国家自然科学基金 41974041

装备预研基金 61405180103

详细信息
    作者简介:

    金彪,博士生,高级工程师,主要从事卫星导航完好性和卫星定轨研究。jinbiao366788@126.com

  • 中图分类号: P228

Analysis of SBAS Ephemeris Correction and UDRE Generation Algorithm

Funds: 

The National Natural Science Foundation of China 41974041

Equipment Preresearch Foundation of China 61405180103

More Information
    Author Bio:

    JIN Biao, PhD candidate, senior engineer, specializes in the integrity of satellite navigation and satellite orbit determination. E-mail: jinbiao366788@126.com

  • 摘要: 星基增强系统(satellite based augmentation system,SBAS)通过地球同步轨道卫星实时播发导航卫星星历改正数和完好性参数,以提升用户定位精度和完好性。采用最小方差法解算GPS星历改正数,利用卡方统计进行改正数完好性检核,并依据星历改正数方差-协方差信息计算SBAS用户差分距离误差(user differential range error,UDRE)和信息类型28(message type 28, MT28)等完好性参数。利用中国区域27个监测站的实测数据,首先以国际GNSS服务组织的精密轨道和钟差产品为参考解算星历改正数,结果表明,钟差改正精度优于0.1 m,轨道改正精度优于0.4 m;然后解算广播星历改正数,并生成UDRE和MT28参数,广播星历残余误差卡方检验值均小于告警门限,保证了改正数的完好性;最后利用生成的改正数进行SBAS定位解算,得到定位结果的水平精度优于0.7 m,垂直精度优于1.0 m,对比GPS单点定位,所提算法的水平和垂直方向精度分别提升了30%和40%。
  • 图  1  GPS PRN01卫星慢变钟差改正数

    Figure  1.  Long-Term Clock Corrections for GPS PRN01 Satellite

    图  2  GPS PRN01卫星慢变轨道改正数

    Figure  2.  Long-Term Orbit Corrections for GPS PRN01 Satellite

    图  3  GPS PRN01卫星钟差快变改正数

    Figure  3.  Fast Clock Corrections for GPS PRN01 Satellite

    图  4  GPS PRN21卫星钟差慢变改正数

    Figure  4.  Long-Term Clock Corrections for GPS PRN21 Satellite

    图  5  GPS PRN21卫星轨道慢变改正数

    Figure  5.  Long-Term Orbit Corrections for GPS PRN21 Satellite

    图  6  GPS PRN21卫星钟差快变改正数

    Figure  6.  Fast Clock Corrections for GPS PRN21 Satellite

    图  7  星历改正数完好性检测结果

    Figure  7.  Integrity Test Results of Ephemeris Corrections

    图  8  卫星UDRE参数随高度角的变化

    Figure  8.  Satellite UDRE Change with Elevation

    图  9  模拟卫星钟故障下GPS PRN21卫星钟差慢变改正数

    Figure  9.  Long-Term Corrections for GPS PRN21 with Satellite Clock Failure

    图  10  模拟卫星机动时GPS PRN21卫星慢变改正数

    Figure  10.  Long-Term Corrections for GPS PRN21 with Satellite Maneuvering

    图  11  GPS单点定位和SBAS定位误差

    Figure  11.  GPS Single Point Positioning and SBAS Positioning Errors

  • [1] Burns J.Wide Area Augmentation System (WAAS) Update[C].SBAS IWG 30, Changsha, China, 2016
    [2] Schempp T.WAAS Development Changes Since Commissioning[C].ICAO GBAS/SBAS Implementation Workshop, Seoul, Korea, 2019
    [3] 吴云, 杨鑫春, 陈慧.多系统多频率的EGNOS系统完备性模拟分析[J].武汉大学学报·信息科学版, 2012, 37(3):269-273 http://ch.whu.edu.cn/article/id/135

    Wu Yun, Yang Xinchun, Chen Hui.Simulation and Analysis of EGNOS System's Integrity Under Multi-system with Multi-frequency[J].Geomatics and Information Science of Wuhan University, 2012, 37(3):269-273 http://ch.whu.edu.cn/article/id/135
    [4] Seynat C, Flament D.European Geostationary Navigation Overlay Service EGNOS Status Update[C].The 23rd International Technical Meeting of the Institute of Navigation, Portland, Oregon, USA, 2010
    [5] Rao K.GAGAN-The Indian Satellite Based Augmentation System[J].Indian Journal of Radio & Space Physics, 2007, 36:293-302
    [6] Saito S.MSAS System Development[C].ICAO GBAS/SBAS Implementation Workshop, Seoul, Korea, 2019
    [7] 楼益栋, 郑福, 龚晓鹏.QZSS系统在中国区域增强服务性能评估与分析[J].武汉大学学报·信息科学版, 2016, 41(3):298-303 doi:  10.13203/j.whugis20140273

    Lou Yidong, Zheng Fu, Gong Xiaopeng.Evaluation of QZSS System Augmentation Service Performance in China Region[J].Geomatics and Information Science of Wuhan University, 2016, 41(3):298-303 doi:  10.13203/j.whugis20140273
    [8] Cao Yueling, Hu Xiaogong, Wu bin.The Wide Area Difference System for the Regional Satellite Navigation System of Compass[J].Science China Physics, Mechanics and Astronomy, 2012, 55(7):1 307-1 315 doi:  10.1007/s11433-012-4746-1
    [9] RTCA SC-159.Minimum Operational Performance Standards for Global Positioning System/Wide Area Augmentation System Airborne Equipment[S/OL].[2013-06-25].https://global.ihs.com/doc_detail.cfm?item_s_key=00233409
    [10] Tsai Y.Wide Area Differential Operation of the Global Positioning System: Ephemeris and Clock Algorithms[D].Stanford: Stanford University, 1999
    [11] 陈俊平, 杨赛男, 周建华, 等.综合伪距相位观测的北斗导航系统广域差分模型[J].测绘学报, 2017, 46(5):537-546 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201705002.htm

    Chen Junping, Yang Sainan, Zhou Jianhua, et al.A Pseudo Range and Phase Combined SBAS Differential Correction Model[J].Acta Geodaetica et Cartographica Sinica, 2017, 46(5):537-546 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201705002.htm
    [12] Juan C.Hughes Aircraft's Architectural Design of the Federal Aviation Administration Wide Area Augmentation System:An International System[J].Acta Astronautica, 1997, 41(4):335-345
    [13] Rouanet L A, Lembachar R, Authié T, et al.New Orbit Determination and Clock Synchronisation Modules for EGNOS[C].The 29th International Technical Meeting of the Institute of Navigation, Portland, Oregon, USA, 2016
    [14] Walter T, Hansen A, Enge P.Message Type 28[C].The 2001 National Technical Meeting of the Institute of Navigation, Long Beach, California, USA, 2001
    [15] Authié T, Trilles S, Fort J C, et al.Integrity Based on MT28 for EGNOS: New Algorithm Formulation & Results[C].The 30th International Technical Meeting of the Institute of Navigation, Portland, Oregon, USA, 2017
    [16] Blanch J, Walter T, Enge P.Evaluation of a Covariance Based Clock and Ephemeris Error Bounding Algorithm for SBAS[C].The 27th International Technical Meeting of the Institute of Navigation, Tampa, Florida, USA, 2014
    [17] Wu J T, Peck S.An Analysis of Satellite Integrity Monitoring Improvement for WAAS[C].The 15th International Technical Meeting of the Institute of Navigation, Portland, Oregon, USA, 2002
    [18] Jefferson D, Bar-Sever Y.Accuracy and Consistency of Broadcast GPS Ephemeris Data[C].The 13th International Technical Meeting of the Institute of Navigation, Salt Lake City, Utah, USA, 2000
    [19] Schempp T R, Peck S R, Fries R M.WAAS Algorithm Contribution to Hazardously Misleading Information[C].The 14th International Technical Meeting of the Institute of Navigation, Salt Lake City, Utah, USA, 2001
    [20] Tang H, Pullen S, Enge P.Ephemeris Type a Fault Analysis and Mitigation for LAAS[C].IEEE/ION Position, Location and Navigation Symposium, Indian Wells, California, USA, 2010
    [21] Hansen A, Walter T, Enge P.GPS Satellite Clock Event on SVN 27 and Its Impact on Augmented Navigation Systems[C].The 11th International Technical Meeting of the Institute of Navigation, Nashville, Tennessee, USA, 1998
  • [1] 刘路, 郭金运, 周茂盛, 鄢建国, 纪兵, 赵春梅.  GNSS广播星历轨道和钟差精度分析 . 武汉大学学报 ( 信息科学版), 2022, 47(7): 1122-1132. doi: 10.13203/j.whugis20200166
    [2] 王兵浩, 周建华, 陈俊平, 张益泽, 王阿昊.  北斗WADS分区综合改正数拟合模型及定位分析 . 武汉大学学报 ( 信息科学版), 2020, 45(4): 565-573. doi: 10.13203/j.whugis20180260
    [3] 赵昂, 杨元喜, 许扬胤, 景一帆, 杨宇飞, 马越原, 徐君毅.  GNSS单系统及多系统组合完好性分析 . 武汉大学学报 ( 信息科学版), 2020, 45(1): 72-80. doi: 10.13203/j.whugis20180425
    [4] 张勇, 王庆, 黄永江.  非差增强信息与差分虚拟观测数据的等价变换理论 . 武汉大学学报 ( 信息科学版), 2019, 44(2): 221-227. doi: 10.13203/j.whugis20170087
    [5] 张益泽, 陈俊平, 杨赛男, 陈倩.  北斗广域差分分区综合改正数定位性能分析 . 武汉大学学报 ( 信息科学版), 2019, 44(2): 159-165. doi: 10.13203/j.whugis20170044
    [6] 宫晓春, 王宇谱, 王宁, 翟树峰, 吕志平.  GPS卫星钟差改正数实时预报算法 . 武汉大学学报 ( 信息科学版), 2018, 43(6): 867-873. doi: 10.13203/j.whugis20160154
    [7] 王甫红, 龚学文, 刘万科.  顾及广播星历误差分离吸收的分米级星载gps实时定轨新方法 . 武汉大学学报 ( 信息科学版), 2015, 40(9): 1230-1236. doi: 10.13203/j .whu g is20130682
    [8] 龚航, 刘增军, 彭竞, 王飞雪.  利用平滑广播星历的GNSS星载钟短稳评估方法 . 武汉大学学报 ( 信息科学版), 2013, 38(7): 837-841.
    [9] 牛飞, 陈金平, 高为广, 张义生.  非正态分布假设下的完好性阈值计算与分析 . 武汉大学学报 ( 信息科学版), 2012, 37(12): 1429-1433.
    [10] 孟领坡, 吴杰, 袁义双.  载波相位差分接收机自主完好性监测研究 . 武汉大学学报 ( 信息科学版), 2011, 36(3): 271-275.
    [11] 宋小勇, 毛悦, 贾小林, 吴显兵.  基于星间测距的分布式自主星历更新算法 . 武汉大学学报 ( 信息科学版), 2010, 35(10): 1161-1164.
    [12] 高书亮, 李锐, 黄智刚.  利用卡尔曼滤波估计导航信号电离层差分改正 . 武汉大学学报 ( 信息科学版), 2010, 35(9): 1021-1023.
    [13] 李作虎, 郝金明, 李建文, 张成军.  QZSS对GPS区域性能增强分析 . 武汉大学学报 ( 信息科学版), 2010, 35(1): 17-20.
    [14] 刘峻宁, 翟传润, 宋嫡儿, 胡华.  基于改进时星差分的室内伪卫星精密单点定位研究 . 武汉大学学报 ( 信息科学版), 2009, 34(1): 105-108.
    [15] 牛飞, 韩春好, 张义生, 常守峰.  导航卫星星载原子钟异常监测分析 . 武汉大学学报 ( 信息科学版), 2009, 34(5): 585-588.
    [16] 牛飞, 韩春好, 张义生.  局部电离层异常对局域增强系统的影响及其监测方法 . 武汉大学学报 ( 信息科学版), 2009, 34(4): 444-448.
    [17] 蔡艳辉, 程鹏飞.  差分GPS水下定位系统的解析法网形分析 . 武汉大学学报 ( 信息科学版), 2008, 33(8): 824-827.
    [18] 范磊, 符养, 杜晓勇, 李黄.  雾灵山山基掩星观测反演误差分析 . 武汉大学学报 ( 信息科学版), 2008, 33(1): 89-92.
    [19] 刘晖, 刘经南, 祁芳.  广域差分GPS定位算法在通用DSP系统上的实现 . 武汉大学学报 ( 信息科学版), 2005, 30(7): 597-600.
    [20] 金彪, 陈姗姗, 李祝莲, 李语强, 李子潇.  SBAS GEO卫星URE精度及定位增强研究 . 武汉大学学报 ( 信息科学版), 0, 0(0): 0-0. doi: 10.13203/j.whugis20210091
  • 加载中
图(11)
计量
  • 文章访问数:  1257
  • HTML全文浏览量:  361
  • PDF下载量:  70
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-01-03
  • 刊出日期:  2021-01-05

SBAS星历改正数及UDRE参数生成算法分析

doi: 10.13203/j.whugis20190033
    基金项目:

    国家自然科学基金 41974041

    装备预研基金 61405180103

    作者简介:

    金彪,博士生,高级工程师,主要从事卫星导航完好性和卫星定轨研究。jinbiao366788@126.com

  • 中图分类号: P228

摘要: 星基增强系统(satellite based augmentation system,SBAS)通过地球同步轨道卫星实时播发导航卫星星历改正数和完好性参数,以提升用户定位精度和完好性。采用最小方差法解算GPS星历改正数,利用卡方统计进行改正数完好性检核,并依据星历改正数方差-协方差信息计算SBAS用户差分距离误差(user differential range error,UDRE)和信息类型28(message type 28, MT28)等完好性参数。利用中国区域27个监测站的实测数据,首先以国际GNSS服务组织的精密轨道和钟差产品为参考解算星历改正数,结果表明,钟差改正精度优于0.1 m,轨道改正精度优于0.4 m;然后解算广播星历改正数,并生成UDRE和MT28参数,广播星历残余误差卡方检验值均小于告警门限,保证了改正数的完好性;最后利用生成的改正数进行SBAS定位解算,得到定位结果的水平精度优于0.7 m,垂直精度优于1.0 m,对比GPS单点定位,所提算法的水平和垂直方向精度分别提升了30%和40%。

English Abstract

金彪, 魏巍, 陈姗姗, 李东俊. SBAS星历改正数及UDRE参数生成算法分析[J]. 武汉大学学报 ( 信息科学版), 2021, 46(1): 111-117. doi: 10.13203/j.whugis20190033
引用本文: 金彪, 魏巍, 陈姗姗, 李东俊. SBAS星历改正数及UDRE参数生成算法分析[J]. 武汉大学学报 ( 信息科学版), 2021, 46(1): 111-117. doi: 10.13203/j.whugis20190033
JIN Biao, WEI Wei, CHEN Shanshan, LI Dongjun. Analysis of SBAS Ephemeris Correction and UDRE Generation Algorithm[J]. Geomatics and Information Science of Wuhan University, 2021, 46(1): 111-117. doi: 10.13203/j.whugis20190033
Citation: JIN Biao, WEI Wei, CHEN Shanshan, LI Dongjun. Analysis of SBAS Ephemeris Correction and UDRE Generation Algorithm[J]. Geomatics and Information Science of Wuhan University, 2021, 46(1): 111-117. doi: 10.13203/j.whugis20190033
  • 为满足民航用户对全球导航卫星系统(global navigation satellite system,GNSS)导航定位精度、完好性、连续性及可用性的需求,很多国家相继建立了星基增强系统(satellite based augmentation system,SBAS),用于提升GNSS基本导航的性能。目前已建成的SBAS包括美国的广域增强系统(wide area augmentation system,WAAS)[1-2]、欧洲地球静止轨道导航重叠服务(European geostationary navigation overlay service,EGNOS)[3-4]、印度GPS辅助地球静止轨道卫星增强导航(GPS aided geo augmented navigation,GAGAN)[5]和日本多功能卫星增强系统(multi-functional satellite augmentation system,MSAS)[6-7]等,中国北斗星基导航增强系统(BeiDou satellite based augmentation system,BDSBAS)、俄罗斯差分改正和监测系统(system of differential correction and monitoring,SDCM)以及韩国增强卫星系统(Korea augmentation satellite system,KASS)也在建设之中。SBAS利用区域监测站数据生成导航卫星星历改正数、电离层改正数及相应的完好性参数,终端利用SBAS改正数及完好性参数进行定位解算,并生成保护级以保证用户定位的完好性。

    SBAS通过播发卫星钟差和轨道改正数,以及用户差分距离误差(user differential range error,UDRE)和改正数协方差信息保证用户定位的完好性[8-9]。受限于SBAS信息播发带宽,星历改正数分为快变和慢变参数。慢变参数用于描述卫星轨道及钟差的缓慢变化,快变参数主要是卫星钟差的快变部分。快变参数播发周期为6 s,慢变参数播发周期为120 s[9]

    利用监测站数据解算星历改正数主要有快照法[10-11]和卫星动力学[12-13]方法。基于卫星动力学模型解算星历改正数,精度较高,但算法较为复杂;快照法的实现相对简单,同时能够快速反应卫星机动。本文采用基于最小方差估计的快照法估计星历改正数,首先利用GPS观测数据解算卫星钟差改正数,并利用经钟差改正后的伪距残差解算卫星轨道改正数,将此时得到的轨道和钟差改正数作为慢变参数;然后利用慢变改正后的伪距残差再次进行卫星钟差解算,将得到的改正数作为快变改正数。

    SBAS向用户播发完好性参数UDRE表示改正数的精度。UDRE是标量[9],若仅向用户提供UDRE参数,则需要播发服务区域内UDRE的最大值,将影响系统的可用性。信息类型28(message type 28,MT28)可对UDRE参数进行调整,用户利用MT28信息恢复改正数方差协方差信息,并根据自身位置对UDRE参数进行调整[14-15],避免了采用区域内UDRE最大值导致可用性降低的问题[16]

    本文利用经过慢变和快变参数改正后的伪距残差解算星历残余误差,并进行完好性检验;通过检验后,基于改正数方差-协方差矩阵生成该卫星的UDRE参数和MT28信息;并利用中国区域实测数据对本文算法进行了验证。

    • GPS双频无电离层伪距观测方程可表示为[11]

      $$\left\{ {\begin{array}{*{20}{l}} {\rho _i^j = {\rho ^0} - \mathit{\boldsymbol{G}}{\rm{\Delta }}\mathit{\boldsymbol{E}} + {\rm{cd}}{t_i} - {\rm{cd}}{t^j} + \varepsilon }\\ {\mathit{\boldsymbol{G}} = \left[ {\frac{{{x^{{\rm{sat}}}} - x}}{{{\mathit{\boldsymbol{r}}^{{\rm{sat}}}} - \mathit{\boldsymbol{r}}}}, \frac{{{y^{{\rm{sat}}}} - y}}{{{\mathit{\boldsymbol{r}}^{{\rm{sat}}}} - \mathit{\boldsymbol{r}}}}, \frac{{{z^{{\rm{sat}}}} - z}}{{{\mathit{\boldsymbol{r}}^{{\rm{sat}}}} - \mathit{\boldsymbol{r}}}}} \right]}\\ {{\rm{\Delta }}\mathit{\boldsymbol{E}} = {{[{\rm{d}}{x^{{\rm{sat}}}}, {\rm{d}}{y^{{\rm{sat}}}}, {\rm{d}}{z^{{\rm{sat}}}}]}^{\rm{T}}}} \end{array}} \right. $$ (1)

      式中,ρij为观测值;ρ0为卫星和测站间的几何距离;G为卫星至测站的视线向量;ΔE为卫星轨道误差;dti为接收机钟差;dtj为卫星钟差;c为光速;ε包括了观测噪声、多路径及模型残余误差等;(xsat, ysat, zsat)为卫星位置;(x, y, z)为测站位置;rsatr分别为卫星、测站的位置向量;(dxsat, dysat, dzsat)为卫星轨道改正数。在解算卫星及测站钟差时,忽略卫星轨道误差,则式(1)可表示为:

      $$\rho _i^j = {\rho ^0} + c{\rm{d}}{t_i} - c{\rm{d}}{t^j} + \varepsilon $$ (2)

      此时卫星钟差吸收了轨道误差在站星视线方向的投影。解算时选定某一测站作为主控站估计接收机钟差,再以主控站为基准解算卫星及其他测站钟差。解算完成后,对式(1)中的伪距残差进行卫星和测站钟差改正,此时认为残余误差仅包括卫星轨道误差,利用最小方差法估计轨道误差,计算如下:

      $${\mathit{\boldsymbol{\hat X}}_{{\rm{MV}}}} = {({\mathit{\boldsymbol{ \boldsymbol{\varLambda} }}^{ - 1}} + {\mathit{\boldsymbol{G}}^{\rm{T}}}{\mathit{\boldsymbol{W}}^{ - 1}}\mathit{\boldsymbol{G}})^{ - 1}}{\mathit{\boldsymbol{G}}^{\rm{T}}}{\mathit{\boldsymbol{W}}^{ - 1}}\mathit{\boldsymbol{z}} $$ (3)
      $${\mathit{\boldsymbol{\hat P}}_{{\rm{MV}}}} = {({\mathit{\boldsymbol{ \boldsymbol{\varLambda} }}^{ - 1}} + {\mathit{\boldsymbol{G}}^{\rm{T}}}{\mathit{\boldsymbol{W}}^{ - 1}}\mathit{\boldsymbol{G}})^{ - 1}} $$ (4)

      式中,${\mathit{\boldsymbol{\hat X}}_{{\rm{MV}}}}$为轨道误差估计值;Λ为广播星历先验轨道精度,取值为径向2.61 m、切向13.25 m、法向5.45 m[17-18],实际使用时需将该误差转到XYZ方向;W为观测值权矩阵;z为消除卫星和测站钟差后的伪距残差。

      上述轨道和钟差改正数每秒解算一次,实际电文播发时需进行120 s拟合得到慢变改正数。完成慢变改正数解算之后,对伪距残差进行慢变误差改正,根据式(2)再次解算卫星钟差,并将其作为快变改正数,每6 s播发一次。本文在解算过程中未对慢变和快变改正数进行拟合,直接使用解算结果作为慢变和快变改正数。

    • 完成SBAS星历改正数解算后,对观测值进行修正得到消除星历误差的伪距残差,并利用此残差进行卫星星历残余误差的解算,计算如下:

      $$\mathit{\boldsymbol{e}} = {({\tilde P^{ - 1}} + {\mathit{\boldsymbol{H}}^{\rm{T}}}\mathit{\boldsymbol{WH}})^{ - 1}}{\mathit{\boldsymbol{H}}^{\rm{T}}}\mathit{\boldsymbol{W}}{\mathit{\boldsymbol{z}}_c} $$ (5)
      $$\mathit{\boldsymbol{P}} = {({\tilde P^{ - 1}} + {\mathit{\boldsymbol{H}}^{\rm{T}}}\mathit{\boldsymbol{WH}})^{ - 1}} $$ (6)
      $$\begin{array}{*{20}{c}} {\tilde P = {\rm{diag}}(\mathit{\Lambda }_X^2 + L_X^2, \mathit{\Lambda }_Y^2 + L_Y^2, \mathit{\Lambda }_Z^2 + L_Z^2, }\\ {{{2.61}^2} + L_{{\rm{clock}}}^2)} \end{array} $$ (7)

      式中,e为星历残余误差;P为星历残余误差的后验方差矩阵;$\tilde P$为广播星历先验方差;ΛXΛYΛZ为广播星历误差在XYZ方向的投影,广播星历钟差先验精度取值2.61 m;(LX, LY, LZ)、Lclock分别为§1.1中解算的轨道、钟差改正数;Hm×4观测矩阵,m为观测值数量,H矩阵中前3列为各测站至卫星的方向余弦,第4列向量元素均为1;W为观测值权矩阵;zcm×1矩阵,是经卫星轨道、钟差以及接收机钟差改正之后的伪距残差。

      利用星历残余误差构造统计量,计算如下:

      $$\left\{ {\begin{array}{*{20}{l}} {{k_{{\rm{state}}}} = \sqrt {{\mathit{\boldsymbol{e}}^{\rm{T}}}{\mathit{\boldsymbol{P}}^{ - 1}}\mathit{\boldsymbol{e}}} }\\ {{k_{{\rm{user}}}} = \frac{{{\mathit{\boldsymbol{u}}^{\rm{T}}}\mathit{\boldsymbol{e}}}}{{\sqrt {{\mathit{\boldsymbol{u}}^{\rm{T}}}\mathit{\boldsymbol{Pu}}} }}} \end{array}} \right. $$ (8)

      式中,u的前3列为用户至卫星的视线向量,第4列为1;当残余误差符合高斯分布时,eT P-1e符合卡方分布。当kstatekFA时,根据柯西-施瓦茨定理有kuserkFA,即星历残余误差投影至用户视线方向引起的测距误差小于kFA倍方差在视线方向的投影。kFA为检测门限,当误警率为1×10-3,自由度为4时,kFA为4.3[17]。因此,当kstatekFA时认为该卫星正常,否则将该卫星标识为未监测状态。

    • SBAS播发星历改正数的同时向用户提供相应的完好性信息,即UDRE参数和MT28信息。用户可利用MT28信息恢复改正数的方差-协方差矩阵,同时结合自身位置对UDRE参数进行调整,提升系统的可用性[9],计算如下:

      $$\left\{ {\begin{array}{*{20}{l}} {{U_{{\rm{UDRE}}}} = {\sigma _{{\rm{UDRE}}}} \cdot {\delta _{{\rm{UDRE}}}}}\\ {{\delta _{{\rm{UDRE}}}} = \sqrt {{\mathit{\boldsymbol{u}}^{\rm{T}}}\mathit{\boldsymbol{Cu}}} }\\ {\mathit{\boldsymbol{C}} = \mathit{\boldsymbol{M}}_{28}^{\rm{T}}{\mathit{\boldsymbol{M}}_{28}}} \end{array}} \right. $$ (9)

      式中,UUDRE为用户根据自身位置确定的UDRE参数;σUDRE为SBAS信息类型2~6中播发的卫星UDRE参数;σUDRE为调整因子;C为改正数方差-协方差矩阵;M28为矩阵C的上三角分解矩阵,即信息类型28的内容[9]

      星历改正数方差-协方差信息Pbroadcast的计算如下:

      $${\mathit{\boldsymbol{P}}_{{\rm{broadcast}}}} = \sigma _{{\rm{UDRE}}}^2\mathit{\boldsymbol{C}} $$ (10)
      $${\mathit{\boldsymbol{P}}_{{\rm{broadcast}}}} = {\left( {\frac{1}{{5.33}}} \right)^2}{({k_{md}} + {k_{FA}})^2}\mathit{\boldsymbol{P}} $$ (11)

      式中,系数5.33是考虑用户在后续定位解算时为满足1.0×10-7/150 s的完好性要求,取5.33倍定位误差作为用户保护级[9]kmd为高斯分布漏警门限值,漏警率取值为4.5×10-10时,kmd为6.13[19]。为保障系统服务的完好性,在计算矩阵P时需假设某一监测站i发生故障,此时解算得到的星历改正数及协方差信息记为eiPi,协方差矩阵Pbroadcast需满足:

      $$\sqrt {{\mathit{\boldsymbol{u}}^{\rm{T}}}{\mathit{\boldsymbol{P}}_{{\rm{broadcast}}}}\mathit{\boldsymbol{u}}} \ge \sqrt {\mathit{\boldsymbol{e}}_i^{\rm{T}}{\mathit{\boldsymbol{P}}_i}{\mathit{\boldsymbol{e}}_i}} $$ (12)

      使用数学方法计算Pbroadcast较为复杂,本文采用文献[17]提出的方法计算Pbroadcast矩阵。假设某一监测站i故障时,SBAS服务区域内用户j处改正数误差σj最大,即:

      $${\sigma _j} = {\rm{max}}\left( {\sqrt {\mathit{\boldsymbol{u}}_j^{\rm{T}}{\mathit{\boldsymbol{P}}_i}{\mathit{\boldsymbol{u}}_j}} } \right) $$ (13)

      利用σj以及所有监测站均无故障时解算得到的方差矩阵P计算扩展系数F0

      $${F_0} = {\rm{max}}\left( {\frac{{{\sigma _j}}}{{\sqrt {\mathit{\boldsymbol{u}}_j^{\rm{T}}\mathit{\boldsymbol{P}}{\mathit{\boldsymbol{u}}_j}} }}} \right) $$ (14)

      此时,系统播发的协方差矩阵计算如下:

      $${\mathit{\boldsymbol{P}}_{{\rm{broadcast}}}} = {\left( {\frac{1}{{5.33}}} \right)^2}{({k_{md}} + {k_{FA}})^2}F_0^2\mathit{\boldsymbol{P}} $$ (15)

      Pbroadcast进行分解得到上三角分解矩阵U以及系统播发的UDRE参数σUDRE和信息M28

      $${\mathit{\boldsymbol{P}}_{{\rm{broadcast}}}} = {\mathit{\boldsymbol{\bar U}}^{\rm{T}}}\mathit{\boldsymbol{\bar U}} $$ (16)
      $$\left\{ {\begin{array}{*{20}{l}} {{\sigma _{{\rm{UDRE}}}} = {{\mathit{\boldsymbol{\bar U}}}_{4, 4}}}\\ {{\mathit{\boldsymbol{M}}_{28}} = \mathit{\boldsymbol{\bar U}}/{{\mathit{\boldsymbol{\bar U}}}_{4, 4}}} \end{array}} \right. $$ (17)
    • 为了验证本文算法的有效性和正确性,使用2018-06-01中国区域27个监测站的实测数据进行GPS改正数的解算并生成完好性参数。

    • 利用IGS精密产品计算卫星位置和钟差,得到测站和卫星间的几何距离及伪距残差并解算改正数。图 1为GPS PRN01号卫星的慢变钟差改正数,慢变钟差改正数均值为-0.02 m,标准差为0.07 m。图 2为GPS PRN01号卫星慢变轨道改正数,XYZ方向的轨道改正数均值分别为-0.19 m、0.39 m及0.17 m,方差分别为0.28 m、0.33 m和0.23 m,三轴平均精度优于0.4 m。在卫星观测开始和结束弧段改正数精度较低,这主要是由于卫星高度角较低导致观测值噪声大以及可视监测站数量少。

      图  1  GPS PRN01卫星慢变钟差改正数

      Figure 1.  Long-Term Clock Corrections for GPS PRN01 Satellite

      图  2  GPS PRN01卫星慢变轨道改正数

      Figure 2.  Long-Term Orbit Corrections for GPS PRN01 Satellite

      利用经过慢变改正后的伪距残差重新计算卫星钟差,并将其作为快变改正数,由于采取每历元解算一次慢变改正数,因此此处快变钟差改正数可视为星历残余误差导致的测距误差。图 3为GPS PRN01号卫星的钟差快变改正数,均值为-0.02 m,方差为0.04 m,其综合影响小于0.1 m。

      图  3  GPS PRN01卫星钟差快变改正数

      Figure 3.  Fast Clock Corrections for GPS PRN01 Satellite

    • 利用GPS广播星历计算卫星慢变和快变改正数及星历残余误差,对星历残余误差进行完好性检验并生成UDRE参数。图 4为利用广播星历解算的GPS PRN21号卫星慢变钟差改正数,改正数小于1 m,与广播星历卫星时钟精度较为一致,且能够反映出卫星钟差的变化趋势。图 5为PRN21号卫星的慢变轨道改正数,XYZ方向改正数统计值分别为(0.27±0.43) m、(0.23±0.72) m和(-0.46±0.64) m,该改正数中包含了广播星历的真实轨道误差。图 6为该卫星的快变改正数,均值为0,标准差为0.07 m,其余卫星情况类似。

      图  4  GPS PRN21卫星钟差慢变改正数

      Figure 4.  Long-Term Clock Corrections for GPS PRN21 Satellite

      图  5  GPS PRN21卫星轨道慢变改正数

      Figure 5.  Long-Term Orbit Corrections for GPS PRN21 Satellite

      图  6  GPS PRN21卫星钟差快变改正数

      Figure 6.  Fast Clock Corrections for GPS PRN21 Satellite

      通过精密星历和广播星历改正数的解算结果可以看出,本文解算的卫星钟差和轨道改正数精度较高,残余误差小于0.1 m,能够保证SBAS系统服务的精度和可用性。

      采用§1.2的方法对星历改正数进行完好性检验并计算UDRE参数。图 7为所有卫星24 h完好性检核结果,红色线为检测门限,绿色点为卡方统计值。由图 7可以看出,统计值均小于检测门限,表明该时段内各卫星改正数均能够满足完好性要求。

      图  7  星历改正数完好性检测结果

      Figure 7.  Integrity Test Results of Ephemeris Corrections

      图 8给出了卫星UDRE值与高度角的关系,红色点表示各卫星的UDRE值,绿色点为该高度角下卫星UDRE的均值,评估时段内UDRE最大值为20.96 m,最小值为1.27 m,对应的UDRE指数分别为12和3[9]。从图 8中可以看出,卫星UDRE值随高度角的增大而减小。

      图  8  卫星UDRE参数随高度角的变化

      Figure 8.  Satellite UDRE Change with Elevation

    • 为了验证本文算法应对卫星故障[20-21]的能力,对GPS PRN21号卫星钟和轨道分别加入误差。首先对卫星钟加入0.5 m/s的钟漂误差以模拟星钟故障,然后在卫星Y坐标上加入0.5 m/s的轨道变化模拟卫星机动。两次误差模拟起始时刻均为475 200 s,持续时间为100 s。

      图 9为加入卫星钟误差后解算的慢变钟差改正数。由图 9可以看出,钟差改正数在475 200 s时能够准确反映出模拟误差并在475 300 s处达到最大值。图 10为加入轨道误差后解算得到的改正数。由图 10可以看出,轨道径向误差被钟差改正数吸收,残余误差体现在轨道改正数上。通过上述分析,本文算法能够有效针对卫星钟故障或轨道机动导致的测距误差进行改正。

      图  9  模拟卫星钟故障下GPS PRN21卫星钟差慢变改正数

      Figure 9.  Long-Term Corrections for GPS PRN21 with Satellite Clock Failure

      图  10  模拟卫星机动时GPS PRN21卫星慢变改正数

      Figure 10.  Long-Term Corrections for GPS PRN21 with Satellite Maneuvering

    • 利用广播星历及§2.2解算的星历改正数进行SBAS定位解算,并与GPS单点定位结果进行比较,电离层误差使用Klobuchar模型进行改正。图 11给出了27个监测站的定位结果。从图 11中可以看出,SBAS水平定位精度为0.7 m,高程定位精度为1.0 m,GPS水平定位精度为1.1 m,高程定位精度为1.8 m。相对于GPS单点定位,利用本文解算的改正数进行SBAS定位时,水平和高程精度分别提升了30%和40%。

      图  11  GPS单点定位和SBAS定位误差

      Figure 11.  GPS Single Point Positioning and SBAS Positioning Errors

    • 本文利用中国区域27个监测站数据进行GPS星历慢变、快变改正数解算,首先利用精密星历验证了算法的正确性,利用广播星历解算改正数进行完好性检测,并生成UDRE参数,结果显示改正后的星历残余误差小于0.1 m,卫星改正数能够满足完好性要求,UDRE参数随卫星高度角增大而减小;在广播星历的基础上加入误差以仿真卫星钟故障或轨道机动,结果表明,本文算法能够有效改正由卫星异常引起的测距误差。利用GPS广播星历及改正数进行SBAS定位解算,与GPS单点定位相比,SBAS定位在保证完好性的同时能够有效提升用户定位精度,水平和高程方向分别提升了30%和40%。

参考文献 (21)

目录

    /

    返回文章
    返回