做GNSS定位的朋友应该都遇到过这种情况明明是晴天、卫星数也足够固定解却突然“飘”出去几米或者单点定位精度肉眼可见地变差。如果恰好又赶上太阳活动活跃的年份这种情况会更频繁甚至连续几天都没法达到厘米级。十有八九问题出在电离层延迟上。电离层延迟估计是GNSS高精度定位绕不开的核心环节。无论你是做测量测绘、无人机自主作业还是搞自动驾驶、地基增强系统只要遇到单频接收机或长基线差分就会被它“上一课”。这篇内容我想从一个干了多年GNSS数据处理的从业者角度把电离层延迟估计这个事讲透从物理机制、数学模型到实际软件配置和踩坑经验尽量让刚入行的朋友能看懂也让有经验的朋友能对照查漏补缺。1. 先搞清楚电离层延迟是怎么“拖累”定位的1.1 信号的“减速”与“超前”折射效应到底是什么电离层是地球大气层中距离地面约60到1000公里的一段区域这里的空气分子被太阳辐射电离形成了大量自由电子和正离子。GNSS卫星信号从约2万公里高空传到地面必须穿过这片区域而电磁波在含有自由电子的介质中传播时会发生折射速度发生变化这就是电离层延迟产生的根本原因。这里有个很有意思的细节对不同的信号成分电离层的影响方向是相反的。GNSS信号包含载波相位和伪码两种测量值载波相位在电离层中传播速度比真空中快结果就是载波相位测距结果“变短”而伪码的群速度比真空中慢伪距测距结果会“变长”。二者符号相反大小基本相同这个特性在后面很多数据处理策略中非常关键。从数学上可以表达为这样在GNSS使用的L波段1到2GHz电离层对伪距和载波相位的一阶延迟量与信号频率的平方成反比比例系数约为40.3。也就是说频率越低延迟越大。GPS的L1频率是1575.42MHzL2频率是1227.60MHz同样一段路径上L2的延迟大约是L1的1.65倍。正是这个频率相关性给了我们消除电离层延迟的“抓手”。1.2 量级有多大什么时候最严重很多人对电离层延迟量级没有直观概念我举个例子在太阳活动平静期的夜间天顶方向卫星在头顶的电离层延迟大约只有1到3米而在白天、特别是太阳活动高峰期的下午天顶方向延迟可以达到10到30米。对于低仰角卫星信号斜穿电离层的路径更长延迟会被放大数倍极端情况下超过50米也不稀奇。影响电离层延迟大小的因素有几个一是太阳活动强度太阳黑子数、耀斑爆发都会显著改变电离层电子密度二是地方时白天太阳辐射强电子密度大延迟明显高于夜间三是地理纬度低纬度和赤道附近电离层最活跃高纬度经常受到极区扰动影响四是季节春秋分前后电离层电子密度往往更高。还有一个容易被忽视的点电离层延迟的空间梯度。通过两个接收机做差分定位时如果不能完全消除电离层残差基线越长、梯度越明显残差就越大。这也是为什么长基线RTK比短基线RTK难做很多时候不是模糊度解算的问题而是电离层残差把整周模糊度“带偏”了。1.3 为什么说单频和双频用户“待遇不同”单频接收机只有L1观测值没有第二个频率去“做差值”所以只能依赖外部模型或先验信息来估计电离层延迟。经典的广播模型Klobuchar模型就是干这个的GPS广播电文里带了8个参数用户用这8个参数去计算当前电离层延迟但模型通常只能消除50%到60%的电离层误差剩下的残余量在太阳活动高峰期可能仍有5到15米。双频接收机就不一样了。因为电离层延迟和频率平方成反比用L1和L2两个频率的观测值做线性组合就能把一阶项几乎完全消掉——消掉99.9%以上。这就是双频消电离层组合的基本思想也是过去几十年高精度静态定位、RTK和PPP最主流的处理方式。随着多系统多频率的普及现在的趋势更灵活了可以用L1/L5、L2/L5等不同频率组合还可以用三频甚至更多频率做无电离层组合或者干脆不消电离层而是把延迟当作参数直接估计出来。这几种思路没有绝对优劣关键在于使用场景和算力约束。2. 电离层延迟估计的主流技术路线2.1 经验模型不依赖实测也能“估个大概”所谓经验模型就是用长期观测数据拟合出的电离层变化规律不需要当前实测数据就能给出一个延迟估计值。这类模型在实时单频定位里使用最广因为接收机没有第二个频率也没有外部改正服务只能用广播电文里的参数。最典型的就是Klobuchar模型。它的设计思路非常巧妙把夜间电离层延迟假设为一个常数大约5纳秒约1.5米白天用一个余弦函数来描述延迟随地方时的变化余弦函数的振幅和周期则根据广播参数计算。之所以这么设计是因为GPS系统设计时希望用尽量少的比特位传输参数8个参数就能让全球用户在单频接收机上大致改正电离层延迟在当时是很经济实用的方案。Galileo系统使用的是NeQuick模型这是一种三维电子密度模型比Klobuchar复杂一些可以描述电子密度随高度、纬度、经度和地方时的变化。使用NeQuick模型用户需要知道大致的地理位置和日期时间模型会给出总电子含量再换算成延迟。实测下来NeQuick在中纬度的精度略好于Klobuchar但这两种模型的本质都只是“统计意义上的估算”无法反映当天电离层的真实异常结构比如赤道异常、电离层暴等。2.2 双频消电离层组合把延迟项“物理消除”既然一阶电离层延迟与频率平方成反比那用两个频率的观测值做一个加权相减就可以把共同的延迟项消去。以GPS L1/L2为例消电离层伪距组合可以写成P_if (f1² × P1 - f2² × P2) / (f1² - f2²)对载波相位观测值也做同样的组合注意符号相反的处理。组合之后电离层一项就被消掉了但同时也会带来一个副作用组合观测值的噪声会被放大。具体来说伪距IF组合噪声大约是原始L1伪距噪声的3倍左右相位IF组合噪声也类似。这是因为加权系数都大于1把噪声也按比例放大了。因为这个噪声放大效应实际处理中并不是“有了双频就必须消电离层”。在短基线RTK里由于基线两端电离层延迟强相关差分后残差本来就很小用不用IF组合影响不大反而因为噪声放大可能降低固定成功率。所以很多RTK算法在短基线下用无电离层组合之外的原始频率组合而在长基线或PPP场景下才用IF组合。2.3 非差非组合PPP把电离层当参数直接估近年来随着PPP技术从“双频消电离层”向“非差非组合”演进越来越多研究人员把电离层延迟当成一个未知参数直接估计。这种做法有几个明显好处第一能够保留原始频率观测值和模糊度的整数特性更有利于模糊度固定第二可以更细致地刻画电离层的变化特征第三单频PPP也能在无外部产品的情况下通过估计电离层参数来改善定位精度。但“直接估计”并不是每个历元都把斜路径电离层延迟当作独立参数那样参数太多、无法解算。实际做法通常是把斜延迟投影成垂直方向总电子含量VTEC用随机游走过程或白噪声模型来约束它的时间变化或者引入外部GIM全球电离层地图作为先验约束。这样既能让电离层参数有一定的自由度去吸收真实变化又不会因为参数过密而使法方程奇异。在软件实现上RTKLIB较新的版本在单频PPP里就支持这种电离层估计模式默认用随机游走模型并可以输入外部IONEX文件作为约束。多频多系统PPP中也有类似设计比如用L1、L2、L5多个频率联立观测方程同时估计位置、钟差、模糊度、对流层湿延迟和电离层延迟。2.4 格网电离层产品区域增强的工程化方案除了接收机端用模型和组合消除还有一个工程化思路是“由外部系统提供电离层改正数”格网电离层产品就是这类方案的代表。IGS国际GNSS服务每天发布全球电离层地图GIM以2小时为时间间隔、2.5°经度×5°纬度的格网形式提供垂直总电子含量用户插值得到当前区域的VTEC再结合投影函数转换为斜路径延迟。地基增强系统也广泛使用这类格网。比如WAAS、欧洲EGNOS以及国内一些区域增强系统都通过监测站估算电离层延迟生成区域格网改正数通过卫星或互联网播发给用户。单频用户接收这些改正数后定位精度可以从模型改正的米级提升到亚米级甚至分米级。这种“区域格网”方案对监测站密度和分布要求很高尤其是在电离层梯度大的区域。如果监测站间距过大格网插值就会失真。我见过某区域增强系统的测试数据低纬度地区电离层梯度大时格网改正效果还不如Klobuchar——这不是系统不行而是格网分辨率不够电离层结构变化太快插值根本追不上。3. 关键实现细节与数学基础3.1 从观测方程看电离层延迟的位置要想把电离层延迟估计做好先得把它在观测方程里的位置搞清楚。GNSS伪距和载波相位的观测方程可以简化为P ρ c·(dt - dT) T I B ε_P L ρ c·(dt - dT) T - I N·λ b ε_L其中ρ是卫星到接收机的几何距离dt和dT分别是接收机和卫星钟差T是对流层延迟I是电离层延迟B和b分别是伪距和相位的硬件延迟N是整周模糊度λ是波长。注意I在伪距和相位方程里的符号正好相反。这个符号关系非常重要如果程序里把符号弄错了双频IF组合非但消不掉电离层反而会引入2倍的电离层延迟误差。我在代码审查时就发现过有人把相位方程中的负号漏掉导致高精度定位结果完全发散。另外还要注意这里I表示的是给定频率上的斜路径电离层延迟它与垂直方向总电子含量之间有一个投影关系。常见的投影函数是卫星高度角的三角函数关系不过在低高度角时投影函数的误差会显著增大所以很多精密数据处理中会把高度角截止角设到10°甚至15°以上。3.2 DCB藏在TEC里的“捣乱分子”如果要做电离层TEC监测或者向用户提供格网电离层产品一定会遇到一个绕不开的参数DCB差分码偏差。不同频率的伪距观测值在卫星端和接收机端的硬件延迟并不同这个差值就叫差分码偏差。它们会混进电离层延迟估计里如果不扣除算出来的TEC和真实值之间就会有一个系统性偏差。问题在于DCB和电离层延迟在数学上是强耦合的。在同一个测站、同一颗卫星的数据里只能观测到“电离层延迟 卫星DCB 接收机DCB”的组合无法直接分开。解决这个问题通常需要引入基准约束比如假设全球所有卫星DCB之和为零或者通过多站联合解算把卫星DCB和接收机DCB分离出来。对于普通定位用户来说DCB的影响往往体现在“电离层估计值有偏”进而可能影响模糊度固定和定位结果。对于做电离层研究和产品的开发者DCB处理就更加关键了。IGS发布的GIM产品一般已经扣除了DCB用户用的是“纯净”的电离层产品而自己解算时一定要检查DCB是否处理正确否则反演出的VTEC可能出现负值或不合理跳变。3.3 几种估计策略的精度对比我把几种策略在实际测试中的典型表现做一个粗略对比。注意这是基于中等纬度、中等太阳活动条件下的静态观测数据具体数值会随环境变化但趋势具有参考价值。策略典型水平精度典型垂直精度收敛/固定时间电离层残差量级单频 Klobuchar2~5 m4~8 m即时可见的模型残余双频 IF 静态PPP0.5~2 cm1~3 cm20~40分钟基本消除非差非组合 PPP估计电离层0.5~2 cm1~3 cm15~30分钟有约束参数吸收短基线 RTK差分后残差小1~2 cm2~4 cm数秒残余很小单频 PPP 外部GIM约束0.3~1 m0.5~2 m10~30分钟取决于GIM精度从表里可以看出双频IF组合和估计电离层两种方式在“最终精度”上差别不大真正的差别在收敛时间和模糊度固定成功率上。非差非组合PPP因为保留了原始频率和模糊度整数特性通常能够更快固定这也是它成为现在PPP技术主流的原因之一。4. 实操在RTKLIB里怎么设置与评估4.1 单频、双频、PPP分别怎么选RTKLIB是开源GNSS数据处理软件里用得最多的一套很多人第一时间接触电离层处理就是在它的选项里。打开RTKPOST在Options的“Settings1”页签里有一项“Ionosphere Correction”下拉选项分别是None、Broadcast、SBAS、IFLC消电离层组合、Estimate、IONEX。单频定位时如果做的是标准单点定位SPP选择Broadcast就是使用Klobuchar模型RTKLIB会读取GPS广播星历里的电离层参数自动计算。如果手里有IONEX格式的外部产品也可以选择IONEX通过插值得到每个观测值的斜路径电离层延迟。实测下来用IGS的最终GIM产品比Klobuchar明显好很多只是需要提前下载文件。双频定位时RTK/PPP模式下选择IFLC就会用L1/L2的消电离层组合。这个选项下RTKLIB会直接使用双频组合观测值解算电离层一项被组合公式天然消掉不需要外部参数。静态测量、长基线RTK和双频PPP都可以这么设置。如果是单频PPPRTKLIB还支持“Estimate”选项把电离层垂直延迟作为参数采用随机游走过程估计。使用这个模式时我建议同时加载GIM文件作为初始约束不然收敛时间会非常长而且高度角比较低的时候参数可能吸收很大的噪声。实际测试中有GIM约束的单频PPP能在30分钟内达到亚米级没有约束可能要1小时以上还晃来晃去。4.2 解算结果里怎么看电离层改正效果RTKLIB输出文件很多多数人只关心位置解忽略了电离层相关输出。如果要评估电离层改正有没有生效有几种直观方法。一是看解算结果中的残差序列。如果伪距残差里还有明显的系统性趋势尤其是随高度角变化大概率是电离层改正不彻底。二是用RTKLIB调试输出看“实时电离层参数”。在Estimate模式下RTKLIB能够把每个历元估计的电离层延迟输出到文件里单位为米。你可以把这个值画成时间序列白天高、夜间低、太阳活动强时数值大这样的曲线基本就是合理的。如果曲线跳动剧烈或者数值明显超出正常范围说明先验约束或随机游走噪声设置有问题。三是做“自检”实验用同一组观测数据分别选择Broadcast和IFLC跑一遍对比定位结果。如果两者水平坐标差很小说明电离层影响不大如果差出去好几米说明当前时段电离层活跃必须依赖双频消电离层或外部产品。这种对比测试是判断“到底需不需要担心电离层”最直接的方法。5. 常见问题与排查经验5.1 太阳活动高峰期单频定位精度骤降2024年开始太阳活动进入第25个周期的极大期我在处理野外单频采集数据时明显感觉到下午时段单点定位误差会异常增大部分日子甚至无法达到标称精度。排查这类问题的思路比较简单先看当天太阳活动指标如F10.7、太阳黑子数再看电离层TEC地图是否有明显异常最后对照双频接收机或IGS测站数据看是否同步出现异常。如果你手头只有单频接收机又必须在太阳活动高峰期做测量我有几个经验建议尽量把观测时段安排在上午和夜间避开当地下午14到18点多测几个时段、延长观测时间通过多时段取平均抵消部分电离层残差有条件的就把单频数据和多频参考站做差分利用差分后残留电离层较小的特性改善精度。5.2 双频周跳导致IF组合模糊度错位双频消电离层组合有一个隐蔽的坑L1和L2任何一个频率发生周跳IF组合的模糊度都会变化而这个变化不是整数很难通过常规整周模糊度搜索直接修复。处理这种问题通常会在做IF组合之前先做周跳探测最常用的是GF几何无关组合即L1和L2相位之差。GF组合已经消掉了几何距离和钟差只受电离层变化影响电离层在短时间内的变化是平缓的所以GF组合的跳变很容易识别出周跳。实际操作中如果发现RTK或PPP在某个时刻后突然解算失败、残差变大优先检查前后历元的GF组合是否有跳变。如果是数据文件里已经有周跳标记RTKLIB会按新弧段处理这个没问题但如果周跳没被探测出来IF组合就会带着错误的模糊度“硬算”结果自然是发散的。5.3 低纬度地区电离层闪烁与空间梯度低纬度地区是电离层处理的“重灾区”。这里存在所谓的赤道异常区电子密度高、梯度大而且经常发生电离层闪烁——信号幅度和相位快速随机起伏严重时接收机直接失锁。我处理过一组低纬度地区的静态观测数据晚上8点到10点之间信号中断频繁周跳密度是白天的几十倍定位精度完全没法保障。对这个问题的应对思路有几个层次。接收机层面选择支持抗闪烁算法的设备或者使用能够锁定弱信号的机型数据层面缩短采样间隔、加强周跳探测如果采用差分定位尽量缩短基线长度并选择靠近测站的参考站后处理层面可以使用载波相位平滑伪距或者用卡尔曼滤波对闪烁噪声做自适应调整虽然不能完全消除闪烁影响但能显著提升数据利用率。5.4 电离层参数与DCB的耦合陷阱做电离层TEC反演或者单频PPP电离层估计时最常见的问题就是忽略DCB导致结果异常。有一次我帮朋友排查一个自定义电离层监测程序的异常某颗卫星的斜TEC持续偏低甚至出现负值。检查半天后发现卫星DCB改正表过期了同一颗卫星的DCB值换了新的还没导入导致TEC计算结果系统性偏了约2TECU。这里分享一个排查技巧对一个稳定测站把连续多天的TEC时间序列画出来正常情况下每天的形状应该非常相似只是大小随太阳活动略有浮动。如果某天曲线突然出现整体垂直偏移或者某颗卫星的TEC曲线与相邻卫星不一致很可能是DCB参数出了问题。在做单频PPP时也是这样如果估计出的电离层延迟序列有个明显的系统性偏置而在其他测站没有同样现象要先检查接收机DCB是否混入了估计结果。5.5 多系统组合和未来趋势现在的GNSS接收机基本都能收到GPS、BDS、Galileo、GLONASS多个系统多个频率的信号。多系统组合对电离层处理也带来了好处卫星数多几何结构好即使某个系统受到电离层闪烁影响其他系统还能维持定位。另外随着三频接收机普及可以构造多个频率上的无电离层组合或者在不同组合之间选择噪声和残留最小的策略处理方式的灵活性大大增加。从趋势上看电离层延迟估计正在从“当作误差消除”走向“当作信息利用”。比如PPP-RTK技术就把区域电离层改正作为一项增强信息通过参考站网络生成高分辨率的斜路径电离层延迟并实时播发给用户用户在本地用这些改正值能大幅缩短PPP收敛时间。这类技术在电离层活跃地区和低纬度地区尤其有优势是未来精密定位服务的重要方向。我再分享一个自己习惯的做法只要条件允许我就把IGS的GIM产品下载下来存成本地文件作为电离层约束的“保底方案”。不管用不用得上有这个先验信息在手处理各种异常数据时心里都会踏实很多。电离层延迟估计说到底就是在“观测值”“模型”“约束”三者之间找平衡理解了这个逻辑遇到再复杂的情况也不会慌了。
