4. 差分定位原理
RTK(Real-Time Kinematic,实时动态载波相位差分)技术之所以能够实现厘米级定位,其核心并不在于接收机硬件本身有多么特殊,而在于巧妙地利用了空间相关性来消除误差。卫星信号从距地面约2万公里的轨道传播至地面接收机,途中会穿过电离层和对流层,这些大气层对信号传播速度的影响是缓慢变化的。如果两台接收机相距不远(例如20公里以内),它们接收同一颗卫星信号时,所经历的大气延迟几乎相同。差分定位的基本思想,就是利用基准站(已知精确坐标)的观测值,计算出这些公共误差,并在流动站(待定位)的观测值中将其扣除,从而实现高精度定位。本章将从基本思想出发,逐步推导单差、双差和三差观测方程,并深入分析每一级差分对误差的消除效果。
4.1 差分定位的基本思想:从“测量”到“相对测量”
GNSS定位的本质是测量接收机到卫星的伪距(或载波相位),然后通过空间距离交会解算位置。然而,伪距观测值中包含了卫星钟差、接收机钟差、电离层延迟、对流层延迟、多路径效应等多种误差。对于单点定位而言,这些误差只能通过模型或参数估计来削弱,精度通常在米级。差分定位的核心思路是:既然两台接收机距离较近时,卫星钟差、大气延迟等误差具有很强的空间相关性,那么就可以通过作差的方式将这些公共误差直接消除,而不是费力去估计它们。
差分定位的基本架构包含一个坐标精确已知的基准站和一个或多个流动站。基准站连续观测所有可见卫星,并利用其已知坐标计算出每颗卫星的“误差改正数”或“观测值残差”。这些改正数通过数据链(如UHF电台、4G/5G网络)实时播发给流动站。流动站接收改正数后,对其自身的观测值进行修正,再进行定位解算。根据作差方式的不同,差分定位可分为单差、双差和三差,其误差消除能力和数学模型复杂度逐级递增。
4.2 单差观测方程:消除卫星钟差
单差(Single Difference, SD)是指在同一历元,两台接收机(基准站b和流动站r)对同一颗卫星i的观测值作差。以载波相位观测为例,基准站和流动站对卫星i的载波相位观测方程分别为:
φ_b^i = (1/λ)·ρ_b^i + f·δt_b - f·δt^i - (1/λ)·I_b^i + (1/λ)·T_b^i + N_b^i + ε_b^i
φ_r^i = (1/λ)·ρ_r^i + f·δt_r - f·δt^i - (1/λ)·I_r^i + (1/λ)·T_r^i + N_r^i + ε_r^i
其中,λ为载波波长,ρ为站星几何距离,f为载波频率,δt_b和δt_r分别为基准站和流动站接收机钟差,δt^i为卫星钟差,I为电离层延迟,T为对流层延迟,N为整周模糊度,ε为观测噪声。将两式相减得到单差观测方程:
Δφ_b,r^i = (1/λ)·(ρ_b^i - ρ_r^i) + f·(δt_b - δt_r) - (1/λ)·(I_b^i - I_r^i) + (1/λ)·(T_b^i - T_r^i) + (N_b^i - N_r^i) + (ε_b^i - ε_r^i)
观察单差方程可以发现:卫星钟差δt^i被完全消除。这是因为在同一历元,同一颗卫星的钟差对两台接收机的影响完全相同,作差后自然抵消。同时,如果两台接收机距离较近(<20km),电离层延迟I和对流层延迟T的差值(I_b^i - I_r^i)和(T_b^i - T_r^i)会非常小,可以近似忽略或作为残余误差处理。单差方程中仍然保留了接收机钟差之差(δt_b - δt_r),这是一个未知参数,需要在解算时一并估计。
| 误差源 | 单差消除效果 | 残余影响 |
|---|---|---|
| 卫星钟差 | 完全消除 | 无 |
| 接收机钟差 | 转化为单差钟差参数 | 需作为未知数估计 |
| 电离层延迟 | 大幅削弱(短基线) | 残余延迟与基线长度成正比 |
| 对流层延迟 | 大幅削弱(短基线) | 残余延迟与基线高度差相关 |
| 整周模糊度 | 转化为单差模糊度 | 仍为整数,但需重新估计 |
4.3 双差观测方程:消除接收机钟差
单差虽然消除了卫星钟差,但引入了接收机钟差之差这一新参数。为了进一步简化模型,可以在单差的基础上,对两颗不同卫星(卫星i和卫星j)的单差观测值再次作差,得到双差(Double Difference, DD)观测方程:
∇Δφ_b,r^i,j = Δφ_b,r^i - Δφ_b,r^j
= (1/λ)·(ρ_b^i - ρ_r^i - ρ_b^j + ρ_r^j) + (N_b^i - N_r^i - N_b^j + N_r^j) + ε_dd
双差方程中,接收机钟差之差(δt_b - δt_r)被完全消除,因为它在两颗卫星的单差方程中是相同的,作差后抵消。同时,卫星钟差在单差阶段已被消除,因此双差观测方程中不再包含任何钟差项。这是双差模型最突出的优点:未知参数大幅减少,仅包含双差模糊度(N_b^i - N_r^i - N_b^j + N_r^j)和流动站坐标。
双差模糊度仍然保持整数特性,这是RTK能够实现快速固定解的关键。在实际解算中,通常先通过最小二乘或卡尔曼滤波获得浮点解,然后利用整数搜索算法(如LAMBDA方法)将双差模糊度固定为整数,最后代入方程求解高精度坐标。双差模型是RTK定位中最常用的数学模型,绝大多数商用RTK接收机都采用双差解算引擎。
4.4 三差观测方程:消除整周模糊度
双差方程中仍然包含整周模糊度这一未知参数。如果接收机在观测过程中发生周跳(信号失锁后重新锁定,导致整周计数跳变),模糊度参数就需要重新初始化。为了规避这一问题,可以对相邻历元(t和t+1)的双差观测值再次作差,得到三差(Triple Difference, TD)观测方程:
δ∇Δφ_b,r^i,j(t,t+1) = ∇Δφ_b,r^i,j(t+1) - ∇Δφ_b,r^i,j(t)
= (1/λ)·[ρ_b^i(t+1) - ρ_r^i(t+1) - ρ_b^j(t+1) + ρ_r^j(t+1)
- ρ_b^i(t) + ρ_r^i(t) + ρ_b^j(t) - ρ_r^j(t)] + ε_td
由于在连续观测且未发生周跳的情况下,整周模糊度保持不变,因此三差方程中模糊度项被完全消除。三差方程仅包含流动站坐标(隐含在几何距离ρ中)和观测噪声,是纯粹的几何观测方程。三差解的精度低于双差固定解(因为模糊度未被固定为整数),但它的优势在于:对周跳不敏感,不需要模糊度初始化,非常适合用于动态定位的初始阶段或周跳频繁发生的恶劣环境。
在实际RTK处理流程中,三差通常不用于最终定位,而是用于:① 探测和修复周跳;② 在模糊度固定失败时提供降级解;③ 作为模糊度浮点解的初始值。三差解虽然精度有限(通常为分米级),但它的鲁棒性使其成为RTK算法中不可或缺的“安”。
4.5 案例拆解:城市峡谷环境下的差分策略选择
某测绘单位在市中心进行管线测量,作业环境为高楼密集区,GNSS信号遮挡严重,且多路径效应明显。作业要求平面精度优于3cm,但观测条件恶劣,频繁发生信号失锁和周跳。项目组最初采用标准双差RTK模式,结果发现:由于周跳频繁,模糊度反复初始化,固定成功率仅约60%,作业效率极低。
项目组分析后调整策略:首先利用三差观测值进行周跳探测与修复,在数据预处理阶段将周跳标记并修复,确保进入双差解算器的数据是连续的。同时,在模糊度固定阶段,采用“部分固定”策略——只固定连续观测时间较长、质量较高的卫星对的模糊度,对于频繁失锁的卫星对保持浮点解。经过调整后,固定成功率提升至92%,平面精度达到2.5cm,作业效率提升约50%。
这一案例说明:差分层级的选择不是越高级越好,而是要根据实际观测条件灵活组合。三差虽然精度不高,但在数据质量控制和周跳处理方面具有不可替代的价值;双差是精度核心,但需要高质量的数据支撑。
4.6 差分定位的误差消除效果对比与总结
为了更直观地理解各级差分的误差消除能力,下表汇总了单差、双差和三差对主要误差源的消除效果:
| 误差源 | 单差(站间) | 双差(站间+星间) | 三差(站间+星间+历元间) |
|---|---|---|---|
| 卫星钟差 | 消除 | 消除 | 消除 |
| 接收机钟差 | 保留(需估计) | 消除 | 消除 |
| 电离层延迟 | 大幅削弱 | 大幅削弱 | 大幅削弱 |
| 对流层延迟 | 大幅削弱 | 大幅削弱 | 大幅削弱 |
| 整周模糊度 | 保留(整数) | 保留(整数) | 消除 |
| 周跳影响 | 需重新初始化 | 需重新初始化 | 不敏感 |
| 观测噪声 | 放大√2倍 | 放大2倍 | 放大2√2倍 |
| 典型定位精度 | 分米级 | 厘米级(固定解) | 分米级 |
从表中可以看出,双差在精度和模型复杂度之间取得了最佳平衡,因此成为RTK定位的标准模型。三差虽然消除了模糊度,但噪声放大更严重,精度反而不如双差固定解。在实际RTK算法设计中,通常采用“三差预处理 + 双差解算”的组合策略:先用三差进行周跳探测和数据清洗,再用双差进行高精度定位解算。
此外,需要特别强调的是,差分定位的误差消除效果高度依赖于基线长度。当基线超过20-30公里时,电离层和对流层的空间相关性显著下降,差分残余误差迅速增大,RTK的固定成功率会明显降低。对于长基线场景,需要引入电离层加权模型或非差非组合模型,这将在后续章节中详细讨论。