高精度BDS空间信号完好性评估方法.pdf

第40卷 第3期2 0 2 3 年 6 月沈 阳 航 空 航 天 大 学 学 报Journal of Shenyang Aerospace UniversityVol.40 No.3Jun.2 0 2 3高精度BDS空间信号完好性评估方法王尔申1a，1b，刘慧超1a，雷虹2，韩琳3，宋建1b，徐嵩1a（1.沈阳航空航天大学 a.电子信息工程学院，b.辽宁通用航空研究院，沈阳 110136；2.沈阳飞机设计研究所 电磁环境效应航空科技重点实验室，沈阳 110035；3.中国人民解放军131423部队 电磁环境效应航空科技重点实验室，沈阳 110163）摘要：针对北斗卫星导航系统精密单点定位空间信号完好性评估方法研究中存在的问题，构建了一种高精度BDS空间信号完好性性能评估方法。通过轨道时钟误差恢复、空间信号异常排除以及建立误差包络模型，分析了无故障情况下的BDS-2和BDS-3高精度空间信号测距误差的完好性。为准确地剔除有故障的空间信号，提出了一种故障趋势与阈值组合判断的异常空间信号剔除方法，能比经验阈值法更准确剔除空间信号异常，提升了空间信号完好性性能评估精度。通过对2021年1月至2022年1月期间法国空间研究中心的实际数据进行分析，给出了高精度用户测距精度的评估结果，并分析了其值在轨道间显著差异性和星间差异性。关键词：北斗卫星导航系统；空间信号；轨道误差；完好性信息；用户测距精度中图分类号：TN967.1 文献标志码：Adoi：10.3969/j.issn.2095-1248.2023.03.007Assessment method of high-precision BDS signal-in-space integrity WANG Er-shen1a，1b，LIU Hui-chao1a，LEI Hong2，HAN Lin3，SONG Jian1b，XU Song1a（1a College of Electronic and Information Engineering，1b Liaoning General Aviation Academy，Shenyang Aerospace University，Shenyang 110136，China；2.Key Laboratory of Aeronautical Science and Technology for Electromagnetic environmental Effects，Shenyang Aircraft Design and Research Institute，Shenyang 110035，China；3.Computer Networks and Applications，Unit 31423 of the Chinese People s Liberation Army，Shenyang 110163，China）Abstract：To address the problem of whether the current BDS satellite navigation system high-precision spatial signal ranging accuracy can meet users requirements，A method was constructed to evaluate the signal-in-space integrity performance of high-precision BDS by means of orbital and clock recovery，signal in space anomaly rejection and error envelope modeling.An abnormal signal in space rejection method with combined fault trend and threshold judgment was proposed，which could accurately reject the start and end moments of signal in space abnormalities and improve the accuracy of signal in space integrity performance assessment compared with the empirical threshold method.By processing the actual data of CNES from January 2021 to January 2022，the evaluation results of the high-pre收稿日期:2022-07-01基金项目:国家自然科学基金(项目编号:62173237);民机专项子课题(项目编号:01020220627066-3);辽宁省应用基础研究计划项目(项目编号:2022020502-JH2/1013);辽宁省教育厅科研项目(项目编号:LJKZZ20220031);沈阳市科技计划项目(项目编号:22-322-3-34)作者简介:王尔申(1980-),男,辽宁辽阳人,教授,博士,主要研究方向:北斗/GNSS卫星导航定位理论与应用,E-mail:wanges_。文章编号：2095-1248（2023）03-0056-07王尔申，等：高精度BDS空间信号完好性评估方法第 3 期cision user ranging accuracy were given，and the values were analyzed for significant inter-orbit variability and weak inter-satellite variability.Key words：BeiDou navigation satellite system（BDS）；signal in space（SIS）；orbit error；integrity information；user ranging accuracy（URA）为满足广大实时用户导航定位需求，国际全球导航卫星系统服务于 2007年 6月正式启动 实 时 计 划 项 目（Real-Time Pilot Project，RTPP），并于2013年4月正式提供实时精密数据1。服务端提供的是基于状态空间表述（State Space Representation，SSR）的改正数，如卫星钟差改正、卫星轨道改正等这些误差共同构成了卫星导航定位的状态空间。为保障空域内飞行安全和无人机集群协同对抗的应用2，实时精密单点定位服务端在播发实时轨道和时钟数据的同时需提供相应的完好性信息来保证空间信号性能。为描述高精度空间信号测距误差（Signal In Space Range Error，SISRE）的大小，用户测距精度（User Ranging Accuracy，URA）提供了未知大小的高精度 SISRE 的保守估计。该参数虽然在导航信息中可用，但在提供给航空用户之前需经过空中导航服务提供商的验证或调 整3。已 有 一 些 针 对 GPS、GLONASS、GALILEO 和 BDS 的 SISRE 特性的研究工作。文献 4-5 验证了GPS的完好性支持消息（Integrity Support Message，ISM）参 数。Wang等6通过对比2013-2017年广播星历与精密星历，得到相应的SISRE大小并分析其误差包络情况。结果表明，2 m的URA不能完全满足在评估期间内的系统性能，但以2.4 m为阈值的URA可以满足且更适合用户。Chen等7分析了2016年1月以来大约5年的北斗广播星历和武汉大学的精密数据，结果表明，与BDS-2卫星相比，BDS-3卫星的广播轨道误差性能显著提高。Wang等8分析了 GPS卫星的 SISRE的包络参数，结果表明由于 LEO 卫星 GRACE FO-1的几何结构因素影响，计算其平均情况下的包络标准差在分米级。由于IGS不提供高精度URA信息，并且目前在高精度定位中，对于高精度URA的评估也比较少。为此，本文的研究重点是对实时轨道和时钟数据的长期特征进行统计，计算并分析高精度 SISRE 的大小并确定其可能分布的合适阈值，提供 BDS-2和 BDS-3的高精度 URA值的初步分析与估计。1评估方法1.1预处理策略预处理策略每年有135 000条导航信息因数据记录错误而破坏，11个空间信号的假异常对应一个真异常9。所以在计算SISRE之前，从原始数据中删除数据记录错误是一项重要的预处理步骤。由于任意一种数据错误或丢失都会导致此时间段内无法对其进行评估10，所以满足以下任意一个条件的数据将剔除同时间段其他数据：未接收到实时轨道、时钟改正数；精密星历/时钟丢失或设置为无效值；不在2小时内的广播星历；完好性状态标识不为零；导航消息不合理。1.2根据轨道改正数信息恢复精密轨道根据轨道改正数信息恢复精密轨道目前，BDS实时数据（CLK93/SSRA00CNE0）在CNES中的采样间隔为5 s。因此，理论上每颗卫星每天将收到 17 280次修正11。实时轨道校正信息包含径向、切向和法向校正参数，实时用户根据SSR信息中的数据龄期来选择相应的卫星轨道参数12。SSR轨道电文给出的是卫星星固坐标系下的径向、切向和法向坐标改正分量。卫星在t时刻的轨道改正数可以57沈 阳 航 空 航 天 大 学 学 报第 40 卷通过式（1）计算O=OradialOalongOcross+OradialOalongOcross(t-t0)（1）式中：O与O为参考时刻t0的 3个方向改正分量和速度分量，可以在SSR修正信息1 259消息类型中获得。卫星修正位置计算如式（2）所示Xorbit=Xbroadcast-eO（2）式中：e为卫星轨道坐标系的旋转矩阵；Xbroadcast为广播星历计算出的卫星位置。t时刻的卫星钟差改正信息C计算为C=C0+C1(t-t0)+C2(t-t0)2（3）式中：t0为SSR轨道修正信息获得的参考时间；Ci为从 SSR时钟校正消息的多项式系数。将SSR电文计算的钟差信息用于改正广播星历卫星钟差即可恢复精密卫星钟差tsatellite=tbroadcast-Cc（4）式中：tbroadcast为卫星时间由广播时钟参数；tsatellite为由SSR时钟修正信息修正的卫星时间；tbroadcast为广播星历卫星钟差；C为从SSR时钟校正信息中获得的时钟校正；c为光速。1.3高精度高精度SISRESISRE异常排除异常排除SISRE是卫星信号传输过程中空间段的误差源，它可以描述在统计中由星历误差引起空间段误差的不确定性13-14。位于卫星覆盖区内的每个用户拥有不同的视线（Line-of-sight，LOS）矢量，因此将经历不同的瞬时用户测距误差（Instantaneous User Range Error，IURE）15。GPS、SPS、PS将SISRE定义为位于卫星可见范围内用户的所有IURE的平均值。最坏情况下的 URE（Worst User Range Error，WURE）代表特定时间卫星覆盖区内用户的最大 IURE16。实时数据流中会播发URA信息，用来保守地估计SISRE的不确定性17。用于实时 GNSS应用的 RTCM SC-104消息类型定义中说明了消息1261为播发北斗的SSR URA值。由于评估期间空间信号会存在异常情况，空间信号异常部分参与评估时会对完好性信息的计算产生较大影响，从而导致完好性信息URA值异常大而影响评估结果。因此评估期间需对异常值进行剔除。为此，本文提出一种设置故障趋势与阈值组合判断方法。具体设置如下：（1）BDS-3/BDS-2 MEO、BDS-2 IGSO 轨道径向、切向和法向误差阈值分别定为 3、10、6 m，BDS-2 GEO轨道径向、切向和法向误差阈值分别定为15、45、15 m。所有卫星的时钟误差阈值设置为50 ns。评估期间内超过该阈值，即认为发生SIS故障。（2）连续时间内误差序列形成斜坡、正弦等非随机形状的曲线，截除非随机形状曲线起止点后，若残余中任何分量的绝对值超过P1-也认为其发生 SIS故障。其中，P1-是总样本1-的分位数，在本文中，由于残余粗差较小且少，故使用=0.006来增强统计结果的可靠性。2评估方法2.1数据质量分析数据质量分析根据数据预处理策略，通过滤除异常值，获得长期稳定的数据。图 1显示了 BDS-2和BDS-3卫星数据的可用性信息。有289 086个星历、组合数据点，有 28 503个异常数据被滤除，约占9.86%。从图1中可以看出，缺失的数据可以分为两种类型：一种是 BDS-3/BDS-2卫星在相同时间段缺失数据。这部分是由于NTRIP Caster的操作问题、传输网络、地面跟踪站广播的不稳定实时数据流、BNC软件的稳定性等因素导致的本时间段卫星数据不可用问题；另一种是不同时间段的 BDS-3/BDS-2 卫星发生数据缺失，这部分是由于剔除了未满足上述条件的广播、精密星历所导致。图 2给出了研究所需 BDS-2和 BDS-3卫星数据 0100%区间范围内的可用性统计信息。结果表明，BDS-3 MEO的可用率最高，为93.39%。BDS-2 IGSO和BDS-2 MEO卫星的58王尔申，等：高精度BDS空间信号完好性评估方法第 3 期北斗数据可用性分别89.42%和93.39%，其卫星 的 可 用 性 分 别 在 85.43%91.61%与92.16%93.87%。与 MEO 和 IGSO 卫星相比，GEO卫星缺失数据的频率更高，其中 C03卫星可用率仅为 80.43%。与 BDS-2 MEO 卫星相比，BDS-3 MEO 卫星有更多的数据缺失。数据缺失的主要原因是丢失实时轨道、时钟数据。2.2长期高精度长期高精度SISRESISRE特征统计分析特征统计分析精密数据的准确性直接决定了估计高精度URA值的大小，并且不同分析中心提供的数据精度也各不相同。为了准确地估计高精度URA 值的大小，首先对精密数据精度进行评估。利用CNES分析中心的实时轨道和时钟数据，对2021年1月2022年1月的18颗BDS-3 MEO、3颗BDS-2 MEO、7颗BDS-2 IGSO和5 颗 BDS-2 GEO 卫星进行了分析。图 36 提供了北斗各轨道类型1年的RMS统计值。在排除空间信号异常后，卫星的轨道、时钟误差在相对稳定的范围内变化，并呈现时间相关性。其中，BDS-2定轨精度与轨道类型有关。由于GEO卫星静地几何特性和光压模型精度较低，导致定轨精度较其他类型卫星差。虽然 BDS-2/BDS-3 MEO 的广播轨道差异较为明显，但修正之后误差差异值明显缩小。尽图3BDS-2 GEO轨道RMS统计值图1数据可用性信息图2BDS-2和BDS-3数据可用性统计信息59沈 阳 航 空 航 天 大 学 学 报第 40 卷管某些卫星在某时刻的误差分量数值较大，例如 C03、C16、C24，但由于比例和数量较小，这些数值不会显著改变统计结果18。BDS-2的7颗IGSO卫星轨道精度在3个方向上的误差均低于BDS-2的MEO卫星。表 1提供了 BDS-2和 BDS-3轨道误差分量的平均值。与 BDS-2 相比，由于 BDS-3 卫星搭载了星间链路和稳定性更高的卫星时钟，其 SISRE 值也得到了改进。但由于引入实时轨道、时钟修正后，这种差异明显缩小。2.3高斯包络模型效果分析高斯包络模型效果分析图 7列出 BDS-2和 BDS-3卫星完好性信息的估计结果。估计结果以用户格网URE与最坏情况URE的68%的误差包络统计值。最图4BDS-2 IGSO轨道RMS统计值图5BDS-2 MEO轨道RMS统计值图6BDS-3 MEO轨道RMS统计值表1BDS-2和BDS-3轨道误差分量的平均值系统轨道BDS-2 GEOBDS-2 IGSOBDS-2 MEOBDS-3 MEO径向 RMS/m0.5450.1130.0340.052法向 RMS/m1.1660.1660.0760.093切向 RMS/m1.6450.1780.1030.128钟差 RMS/ns1.230.460.290.27SISRE RMS/m0.7130.3490.2570.22160王尔申，等：高精度BDS空间信号完好性评估方法第 3 期坏情况URE具有最大的包络值，它也代表着用户可以受到最安全的保护，但该方法是以牺牲卫星故障概率为前提的。其中，由于 BDS-2 GEO 卫星的静地特性导致轨道、时钟误差过大，所以其体现出最大的68%的误差包络统计值。在统计的所有卫星高精度URA中，C01具有最大的URA值，两种情况下的包络值分别为0.732 m 与 0.961 m。C29 卫星具有最小的URA值，两种情况下的包络值分别为0.146 m与 0.162 m。不同星座的 URA 值有较大的差异，相同星座的不同卫星URA值也有较小的差异。传统事后北斗广播星历中大多数卫星保持着2.43.2 m的URA11，17，通过本研究所介绍的完好性性能评估方法评估高精度URA大多数保持在 0.150.96 m。实验结果表明经过精密单点定位（PPP）的轨道、钟差改正数修正后的高精度URA在理论上较广播星历URA值会有一个比较大的提升，总体精度在 0.150.96 m。3结论完好性监测（Integrity Monitoring，IM）是保障定位服务安全性和可靠性的一种重要方法。本文通过构建一种高精度BDS空间信号完好性性能评估方法，从完好性的角度分析了2021年 1月至 2022年 1月的 BDS-2和 BDS-3的高精度空间信号测距误差的性能。（1）对评估期间内 BDS-2 和 BDS-3 的数据质量进行分析。结果表明，BDS-2 MEO相对 于 其 他 轨 道 类 型 数 据 缺 失 率 较 高，为14.93%。除 C35外，BDS-3 MEO卫星数据可用 性 超 过 85%，平 均 值 为 92.68%。BDS-2 MEO 和 BDS-2 IGSO 卫星的实时数据平均可用率分别为93.39%、89.42%。（2）使用CNES的RTS数据研究了BDS-2和BDS-3轨道、时钟的性能。结果表明，BDS-2 GEO卫星轨道精度最差，径向轨道误差平均值为 0.545 m，法向与切向轨道误差分别为1.166、1.645 m。BDS-2 IGSO径向、法向和切向轨道误差分别为 0.113、0.166、0.178 m；BDS-2 MEO径向、法向和切向轨道误差平均值分别为0.056、0.103、0.113 m；BDS-3 MEO径向、法向和切向轨道误差平均值分别为0.052、0.093、0.128 m。（3）建立了误差包络模型，初步评估高精度空间信号测距精度值。结果表明，不同星座的 URA值有较大的差异，相同星座的不同卫星的URA值也有较小差异性。经过PPP轨道钟恢复后，得到其 URA 精度大多数保持在0.150.96 m。图7BDS-2和BDS-3卫星完好性信息的估计结果61沈 阳 航 空 航 天 大 学 学 报第 40 卷参考文献（References）：1舒宝，王利，张勤，等.SSR延迟下的轨道钟差外推误差及其对多 GNSS实时精密单点定位的影响评估J.测绘学报，2021，50(12):1738-1750.2王尔申，郭靖，宏晨，等.改进目标收益函数的无人机集群空地对抗模型J.南京航空航天大学学报，2021，53(6):888-897.3王尔申，孙彩苗，佟刚，等.基于 PSO 的多星座GNSS垂直保护级优化方法J.北京航空航天大学学报,2021,47(11):2175-2180.4Perea S，Meurer M，Rippl M，et al.URA/SISA analysis for GPS and Galileo to support ARAIMJ.Navigation，2017，64(2)：237-254.5Mistrapau F，Bija B，Cueto-Felgueroso G，et al.GPS SISRE/URA integrity analysis for ARAIMC/Proceedings of the 29th International Technical Meeting of The Satellite Division of the Institute of Navigation（ION GNSS+2016），ION GNSS+，The International Technical Meeting of the Satellite Division of The Institute of Navigation.Portland，USA：Institute of Navigation，2016：1793-1803.6Wang Z P，Shao W，Li R，et al.Characteristics of BDS signal-in-space user ranging errors and their effect on advanced receiver autonomous integrity monitoring performanceJ.Sensors，2018，18(12)：4475.7Chen G，Zhou R Y，Hu Z G，et al.Statistical characterization of the signal-in-space errors of the BDS：a comparison between BDS-2 and BDS-3J.GPS Solutions，2021，25(3)：112.8Wang K，El-Mowafy A，Rizos C.Integrity monitoring for precise orbit determination of LEO satellitesJ.GPS Solutions，2022，26(1):32.9Zhao Y X，Cheng C，Li L，et al.BDS signal-in-space anomaly probability analysis over the last 6 yearsJ.GPS Solutions，2021，25(2)：49.10张熙，刘长建，章繁.四大GNSS广播星历精度评估与对比分析J.武汉大学学报（信息科学版），2022，47(2):208-218.11王朝辉，马下平，严丽，等.北斗三号全球卫星导航系统的广播星历精度评估J.测绘通报，2021(1):59-65，98.12李广源，花向红，贺小星.北斗卫星导航系统空间信号测距误差评估J.测绘科学,2020,45(5):1-6.13Perea S，Meurer M，Pervan B.Impact of sample correlation on SISRE overbound for ARAIMJ.Navigation,2020,67(1):197-212.14田秋凝，曹月玲，胡小工，等.北斗三号系统广域差分服务精度评估J.北京航空航天大学学报,2021,47(5):967-976.15Xue B，Wang H T，Yuan Y B.Performance of BeiDou-3 signal-in-space ranging errors：accuracy and distributionJ.GPS Solutions，2021，25(1)：23.16Walter T，Gunning K，Eric Phelts R，et al.Validation of the unfaulted error bounds for ARAIMJ.Navigation，2018，65(1)：117-133.17Wang S Z，Zhai Y W，Zhan X Q.Characterizing BDS signal-in-space performance from integrity perspectiveJ.Navigation，2021，68(1)：157-183.18WANG Z，HU W，et al.Enhancement Dataset for Low Altitude Unmanned Aerial Vehicle Detection J.Transactions of Nanjing University of Aeronautics and Astronautics,2021,38(6):914-926.（责任编辑：刘划 英文审校：杜文友）62