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

王尔申，刘慧超，雷 虹，韩 琳，宋 建，徐 嵩

（1. 沈阳航空航天大学 a.电子信息工程学院，b.辽宁通用航空研究院，沈阳 110136；2. 沈阳飞机设计研究所 电磁环境效应航空科技重点实验室， 沈阳 110035；3. 中国人民解放军131423部队 电磁环境效应航空科技重点实验室，沈阳 110163）

为满足广大实时用户导航定位需求，国际全球导航卫星系统服务于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 的完好性支持消息（In‐tegrity 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/SSRA00CNE 0）在CNES 中的采样间隔为5 s。因此，理论上每颗卫星每天将收到17 280 次修正［11］。实时轨道校正信息包含径向、切向和法向校正参数，实时用户根据SSR 信息中的数据龄期来选择相应的卫星轨道参数［12］。SSR轨道电文给出的是卫星星固坐标系下的径向、切向和法向坐标改正分量。卫星在t 时刻的轨道改正数可以通过式（1）计算

式中：δO 与δȮ为参考时刻t0的3 个方向改正分量和速度分量，可以在SSR修正信息1 259消息类型中获得。卫星修正位置计算如式（2）所示

式中：e 为卫星轨道坐标系的旋转矩阵；Xbroadcast为广播星历计算出的卫星位置。t 时刻的卫星钟差改正信息δC计算为

式中：t0为SSR 轨道修正信息获得的参考时间；Ci为从SSR 时钟校正消息的多项式系数。将SSR电文计算的钟差信息用于改正广播星历卫星钟差即可恢复精密卫星钟差

式中：tbroadcast为卫星时间由广播时钟参数；tsatellite为由SSR时钟修正信息修正的卫星时间；tbroadcast为广播星历卫星钟差；δC为从SSR时钟校正信息中获得的时钟校正；c为光速。

1.3 高精度SISRE异常排除

SISRE是卫星信号传输过程中空间段的误差源，它可以描述在统计中由星历误差引起空间段误差的不确定性［13-14］。位于卫星覆盖区内的每个用户拥有不同的视线（Line-of-sight，LOS）矢量，因此将经历不同的瞬时用户测距误差（Instantaneous User Range Error， IURE）［15］。GPS、SPS、PS 将SISRE 定义为位于卫星可见范围内用户的所有IURE 的平均值。最坏情况下的URE（Worst User Range Error， WURE）代表特定时间卫星覆盖区内用户的最大IURE［16］。实时数据流中会播发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 卫星发生数据缺失，这部分是由于剔除了未满足上述条件的广播、精密星历所导致。

图1 数据可用性信息

图2 给出了研究所需BDS-2 和BDS-3 卫星数据0～100%区间范围内的可用性统计信息。结果表明，BDS-3 MEO 的可用率最高，为93.39%。BDS-2 IGSO 和BDS-2 MEO 卫星的北斗数据可用性分别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 BDS-2和BDS-3数据可用性统计信息

2.2 长期高精度SISRE特征统计分析

精密数据的准确性直接决定了估计高精度URA值的大小，并且不同分析中心提供的数据精度也各不相同。为了准确地估计高精度URA 值的大小，首先对精密数据精度进行评估。利用CNES分析中心的实时轨道和时钟数据， 对2021 年1 月～2022 年1 月的18 颗BDS-3 MEO、3 颗BDS-2 MEO、7 颗BDS-2 IGSO 和5 颗BDS-2 GEO 卫星进行了分析。图3～6 提供了北斗各轨道类型1年的RMS统计值。

图3 BDS-2 GEO轨道RMS统计值

图4 BDS-2 IGSO轨道RMS统计值

图5 BDS-2 MEO轨道RMS统计值

图6 BDS-3 MEO轨道RMS统计值

在排除空间信号异常后，卫星的轨道、时钟误差在相对稳定的范围内变化，并呈现时间相关性。其中，BDS-2 定轨精度与轨道类型有关。由于GEO 卫星静地几何特性和光压模型精度较低，导致定轨精度较其他类型卫星差。虽然BDS-2/BDS-3 MEO 的广播轨道差异较为明显，但修正之后误差差异值明显缩小。尽管某些卫星在某时刻的误差分量数值较大，例如C03、C16、C24，但由于比例和数量较小，这些数值不会显著改变统计结果［18］。BDS-2的7颗IGSO 卫星轨道精度在3 个方向上的误差均低于BDS-2的MEO卫星。

表1 提供了BDS-2 和BDS-3 轨道误差分量的平均值。与BDS-2 相比，由于BDS-3 卫星搭载了星间链路和稳定性更高的卫星时钟，其SISRE 值也得到了改进。但由于引入实时轨道、时钟修正后，这种差异明显缩小。

表1 BDS-2和BDS-3轨道误差分量的平均值

2.3 高斯包络模型效果分析

图7 列出BDS-2 和BDS-3 卫星完好性信息的估计结果。估计结果以用户格网URE 与最坏情况URE 的68%的误差包络统计值。最坏情况URE具有最大的包络值，它也代表着用户可以受到最安全的保护，但该方法是以牺牲卫星故障概率为前提的。其中，由于BDS-2 GEO 卫星的静地特性导致轨道、时钟误差过大，所以其体现出最大的68%的误差包络统计值。在统计的所有卫星高精度URA 中，C01 具有最大的URA值，两种情况下的包络值分别为0.732 m 与0.961 m。C29 卫星具有最小的URA 值，两种情况下的包络值分别为0.146 m与0.162 m。不同星座的URA 值有较大的差异，相同星座的不同卫星URA值也有较小的差异。

图7 BDS-2和BDS-3卫星完好性信息的估计结果

传统事后北斗广播星历中大多数卫星保持着2.4～3.2 m 的URA［11，17］，通过本研究所介绍的完好性性能评估方法评估高精度URA 大多数保持在0.15～0.96 m。实验结果表明经过精密单点定位（PPP）的轨道、钟差改正数修正后的高精度URA 在理论上较广播星历URA值会有一个比较大的提升，总体精度在0.15～0.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.15～0.96 m。