GNSS电离层延迟估计全解析：从物理机制到RTKLIB工程实践

做GNSS定位的朋友应该都遇到过这种情况：明明是晴天、卫星数也足够，固定解却突然“飘”出去几米，或者单点定位精度肉眼可见地变差。如果恰好又赶上太阳活动活跃的年份，这种情况会更频繁，甚至连续几天都没法达到厘米级。十有八九，问题出在电离层延迟上。

电离层延迟估计是GNSS高精度定位绕不开的核心环节。无论你是做测量测绘、无人机自主作业，还是搞自动驾驶、地基增强系统，只要遇到单频接收机或长基线差分，就会被它“上一课”。这篇内容我想从一个干了多年GNSS数据处理的从业者角度，把电离层延迟估计这个事讲透：从物理机制、数学模型，到实际软件配置和踩坑经验，尽量让刚入行的朋友能看懂，也让有经验的朋友能对照查漏补缺。

1. 先搞清楚：电离层延迟是怎么“拖累”定位的

1.1 信号的“减速”与“超前”：折射效应到底是什么

电离层是地球大气层中距离地面约60到1000公里的一段区域，这里的空气分子被太阳辐射电离，形成了大量自由电子和正离子。GNSS卫星信号从约2万公里高空传到地面，必须穿过这片区域，而电磁波在含有自由电子的介质中传播时会发生折射，速度发生变化，这就是电离层延迟产生的根本原因。

这里有个很有意思的细节：对不同的信号成分，电离层的影响方向是相反的。GNSS信号包含载波相位和伪码两种测量值，载波相位在电离层中传播速度比真空中快，结果就是载波相位测距结果“变短”；而伪码的群速度比真空中慢，伪距测距结果会“变长”。二者符号相反，大小基本相同，这个特性在后面很多数据处理策略中非常关键。

从数学上可以表达为这样：在GNSS使用的L波段（1到2GHz），电离层对伪距和载波相位的一阶延迟量与信号频率的平方成反比，比例系数约为40.3。也就是说，频率越低，延迟越大。GPS的L1频率是1575.42MHz，L2频率是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 几种估计策略的精度对比

我把几种策略在实际测试中的典型表现做一个粗略对比。注意这是基于中等纬度、中等太阳活动条件下的静态观测数据，具体数值会随环境变化，但趋势具有参考价值。

从表里可以看出，双频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都可以这么设置。

如果是单频PPP，RTKLIB还支持“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产品下载下来存成本地文件，作为电离层约束的“保底方案”。不管用不用得上，有这个先验信息在手，处理各种异常数据时心里都会踏实很多。电离层延迟估计说到底就是在“观测值”“模型”“约束”三者之间找平衡，理解了这个逻辑，遇到再复杂的情况也不会慌了。