史俊波 董新莹 欧阳晨皓 彭文杰 姚宜斌
1 武汉大学测绘学院,武汉市珞喻路129号,430079
北斗增强系统是北斗卫星导航系统BDS的重要组成部分,包括地基增强系统和星基增强系统,这2类增强系统对高精度“北斗+”应用的推广至关重要[1]。北斗地基增强系统以基于网络RTK技术的连续运行参考站系统为主,通过地面通信网络,为用户提供实时cm级、事后mm级的定位增强服务[2]。截至2019年底,国内连续运行参考站申请备案数量已达7 223座[3],同时出现了多家覆盖全国的商用北斗地基增强系统[4]。与地基增强系统不同,北斗星基增强系统通过卫星通信网络为用户提供定位增强服务[5],包括BDSBAS服务[6]和PPP服务[7],这2种服务均通过3颗北斗地球静止轨道GEO卫星(C59、C60、C61)向我国及周边地区广播增强信息[8]。2种星基增强服务的区别在于:BDSBAS服务通过BDSBAS-B1C和BDSBAS-B2a增强信号播发伪距差分增强信息,而PPP服务通过PPP-B2b增强信号播发载波相位差分增强信息[9]。
目前已有学者对北斗PPP服务性能展开了评估[10-16],这些研究主要是使用2个独立接收机分别接收导航信号和增强信号,评估场景主要为长时段跟踪站静态场景和模拟动态场景。本文利用一台独立接收机同时接收GNSS导航信号和PPP-B2b增强信号,并设计108个快速静态实验和1个低动态导轨实验,进一步丰富北斗星基PPP服务定位性能评估的流程和手段。
基于一台能同时接收GNSS导航信号和PPP-B2b增强信号的接收机开展BDS-3精密单点定位服务的流程如图1所示,具体步骤包括:
图1 基于PPP-B2b信号的PPP服务使用流程
1)PPP-B2b增强信号和导航信号的获取:通过串口命令设置板卡,接收原始报文,解码得到BDS-3+GPS观测信息、广播星历和PPP-B2b增强信息;
2)PPP-B2b增强信息的使用:使用步骤1)得到的PPP-B2b SSR改正数,与广播星历进行匹配,恢复得到精密的轨道/钟差/DCB产品;
3)精密单点定位的实现:使用步骤1)得到的BDS-3+GPS观测信息,和步骤2)恢复出的PPP-B2b精密轨道/钟差/DCB产品,进行精密单点定位解算。
图2为信号采集与解码平台的实现流程。首先使用单台接收机同时接收PPP-B2b增强信号和导航信号;然后通过串口将数据流传输给电脑,电脑运行自编软件实现PPP-B2b增强信号和导航信号的实时解码,得到GNSS观测值、北斗CNAV广播星历、GPS LNAV广播星历和PPP-B2b增强信息。
图2 增强信号和导航信号同时采集与解码示意图
根据空间信号接口控制文件,PPP-B2b增强信号提供了卫星的轨道、钟差和码间偏差改正数[7]。其中,卫星轨道改正信息δO用于计算卫星位置改正向量δX,以修正广播星历卫星位置Xbroadcast,从而得到精密卫星位置Xorbit:
(1)
式中,ei为方向余弦矢量,i={radial,along,cross}分别对应径向、切向和法向。
PPP-B2b增强信号提供的卫星钟差改正信息C0,使用方法如下:
(2)
式中,tsatellite为精密卫星钟差,tbroadcast为广播星历卫星钟差,c为真空中的光速。
PPP-B2b增强信号还提供卫星信号的码间偏差改正数DCBsig,使用方法如下:
(3)
目前北斗PPP-B2b增强信号仅提供BDS-3和GPS的改正数,表1为本文采用的PPP解算策略。由于PPP-B2b提供的BDS-3和GPS卫星钟差基于B3I和L1/L2消电离层(ionosphere-free,IF)组合观测值,而本文实验采用的观测值分别为BDS-3 B1C/B2a和GPS L1/L2消电离层组合,因此需要改正BDS-3卫星码间偏差,不需要改正GPS卫星码间偏差。
表1 实时星基PPP解算策略
在武汉大学校园内开展实验数据采集工作,天线和导轨放置在测绘学院楼顶,四周环境开阔。本文基于司南K803板卡和华信GPS 1000导航型天线开展实验。
快速静态数据采集时间为2022-05-28~06-05(doy148~156),采样频率为1 Hz,采集现场和观测到的卫星数如图3(a)和3(b)所示,天线在三脚架上保持静止。从UTC 00:00开始,每2 h为1个时段,每天划分出12个时段,在每个时段的开始时刻重启接收机。需要指出的是,本文设计的快速静态测量潜在场景是1 h左右观测时长、cm~dm级定位精度的像控点测量,与现有研究中24 h观测时长、mm~cm级定位精度的静态测量场景侧重点不同。
图3 快速静态/低动态星基PPP数据采集情况
低动态数据采集时间为2022-06-13(doy 164)UTC 04:30~05:31,采样频率为1 Hz,采集现场和观测到的卫星数如图3(c)和3(d)所示。天线在长为2.2 m的电控滑动导轨上按照设定的程序进行运动。首先在起点静止34 min,然后打开导轨开关,使天线以稳定的速度在起点和终点之间运动,完成5次往返后回到起点,动态时段为UTC 05:05~05:31。
由图3可见,2个实验场景都有约20颗BDS-3+GPS卫星能接收到增强信号。另外,相比于快速静态场景,低动态场景多次出现增强卫星数低至5颗的突变。
将连续9 d的快速静态数据每天划分为12个时段,共计108个2 h样本。以连续24 h静态事后PPP坐标为参考,图4为doy148在E、N、U方向的误差序列,每个时段开始都存在因接收机重启导致的重收敛和误差趋于平稳的现象。
图4 基于PPP-B26b信号的快速静态星基PPP误差(doy 148)
将收敛条件定义为平面精度优于10 cm、高程精度优于20 cm、持续10 min以上。除去3个(占比2.8%)不满足收敛条件的样本,针对105个(占比97.2%)可收敛样本的快速静态星基PPP解算结果,统计收敛时间和收敛后的定位精度。图5和6分别为105个样本收敛时间的分时段汇总折线图和直方图,平均收敛时间为34.6 min,59.3%的样本可以在30 min内完成收敛。
图5 105个快速静态星基PPP样本的收敛时间折线
图6 105个快速静态星基PPP样本的收敛时间直方图
图7为105个样本的快速静态星基PPP收敛后的定位精度,其中E、N、U方向的RMS分别为5.2 cm、2.9 cm、9.7 cm,水平及高程方向的RMS分别为5.9 cm、9.7 cm。
图7 105个快速静态星基PPP样本的收敛后的定位精度
为进一步验证北斗星基PPP快速静态定位的性能,本文统计了不同重启时间后全部108个样本水平及高程方向的定位精度及分布情况,如图8和9所示。其中,接收机重启10 min、20 min、30 min后,水平方向RMS分别为9.6 cm、8.0 cm、7.3 cm,优于10 cm的比例分别为78%、83%、86%;高程方向RMS分别为15.5 cm、13.1 cm、11.6 cm,优于15 cm的比例分别为75%、84%、86%。
图8 108个快速静态星基PPP样本的定位精度
图9 全部108个样本的接收机重启后水平及高程方向的定位精度直方图
动态数据共计61 min,划分为2个阶段:1)收敛阶段:34 min(UTC 04:30~05:04),天线在导轨起点保持静止;2)运动阶段:26 min(UTC 05:05~05:31),天线在导轨起点和终点之间进行5次匀速往返运动,最后一次回到起点。
以RTK结果为参考,图10为动态星基PPP解算结果相对于起点E、N、U方向的位移量序列,可以看出位移量在2个阶段与实际采集过程相吻合。
图10 基于PPP-B2b信号的动态星基PPP轨迹(doy164)
针对运动阶段进行动态星基PPP精度分析,图11为动态PPP解算结果和参考结果的平面轨迹对比。可以看出,解算的轨迹存在整体偏移,即星基PPP与地基RTK之间存在系统性偏差。
图11 基于PPP-B2b信号的动态星基PPP平面轨迹
以RTK结果为参考,逐历元统计动态PPP在E、N、U方向的误差序列,如图12所示。由图可见,经过34 min的收敛后,运动阶段的动态星基PPP误差波动较为平稳,E、N、U方向的RMS分别为24.3 cm、1.1 cm、12.2 cm。
图12 基于PPP-B2b信号的动态星基PPP误差
图13为运动阶段参与定位的卫星数和PDOP值变化情况,由图可见,误差序列在UTC 05:06和05:14出现的2处跳动是由于对应时刻卫星数减少和PDOP值发生了变化所致。
图13 运动阶段的卫星数和PDOP值变化情况
1)在快速静态实验中,97%的样本平均收敛时间为34.6 min,收敛后的水平及高程方向RMS可达5.9 cm、9.7 cm;接收机重启10 min、20 min、30 min后,水平方向优于10 cm的比例分别为78%、83%、86%,高程方向优于15 cm的比例分别为75%、84%、86%。
2)在低动态导轨实验中,以RTK结果为参考,30 min的收敛时间后E、N、U方向的RMS分别为24.3 min、1.1 min、12.2 cm。
相比于现有基于改正数地面播发的PPP模式,北斗三号PPP服务提供了一种公开的、无需任何网络成本的定位模式,为北斗高精度定位应用提供了良好的基础设施,后续将研究北斗三号PPP服务在车载导航中的应用。