简介：大气波导高度对气温、气压等气候参数的响应是电磁波超视距传播研究中的关键问题。这份 MATLAB 脚本资源面向通信工程、气象及雷达探测领域的学习者和工程师，提供从折射率计算、大气不均匀性分析到折射角求解的完整流程，并通过主脚本直接绘制高度随多参数变化的趋势图。压缩包为 rar 格式，体积仅 2KB，内部包含 4 个 m 文件，分别承担参数计算与结果整合功能，结构精简，适合快速仿真和二次开发。目前已有 919 人浏览学习。利用该资源，读者可自行调整温压输入，观察大气波导高度与传播条件的变化关系，进而预测无线信号的路径损耗和覆盖距离；所采用的 CCIR 估算思路也便于教学演示与算法验证，对理解波导形成机理、优化通信系统设计具有直接帮助。 搞通信和雷达的朋友，应该都遇到过这种诡异情况：明明设备指标正常、链路余量充足，信号却突然消失了，或者平时收不到的远方信号，忽然强得不像话。这种“神出鬼没”的现象，十有八九是大气波导在作怪。我在做微波链路评估和雷达覆盖分析时，跟大气波导打了多年交道，最头疼的还不是“有没有波导”，而是“波导高度到底在哪”。波导高度直接决定了超视距传播的通道位置、雷达盲区的范围，以及干扰信号的来向。这篇文章就专门聊聊，大气波导高度到底受哪些参数影响、怎么变化、又该怎么算，全程干货，没有废话。
坦率说，很多资料把大气波导讲得玄而又玄，其实核心就那么几个参数：温度、湿度、气压，以及它们随高度的梯度。谁在主导、谁在辅助、谁在捣乱，搞清楚之后，你再看探空数据、再翻数值预报产品，思路会完全不一样。这篇内容适合做微波通信、雷达探测、频谱管理、无线电环境评估的工程师，也适合气象专业想跨界了解电波传播的朋友参考。
1. 先搞懂波导高度到底在说什么
1.1 一个公式看懂折射率剖面
大气波导的本质，是大气折射率随高度异常变化，导致电波射线向下弯曲的程度超过地球曲率，电波被“困”在某一层大气里来回折射传播。要量化这个现象，绕不开两个物理量：大气折射率 N 和修正折射率 M。
大气折射率的经典公式是：
N = 77.6 × (p / T) + 3.73 × 10^5 × (e / T²)
其中 p 是气压（hPa），T 是绝对温度（K），e 是水汽压（hPa）。前一项叫干项，主要反映气压和温度的影响；后一项叫湿项，反映湿度的影响。湿项虽然数值上占比不算最大，但它对高度变化的响应极快，是形成波导的关键推手。
但用 N 判断波导并不直观，因为 N 本身随高度总是减小的，只是减小得快慢不同。为了方便判断，引入修正折射率 M：
M = N + (z / a) × 10^6 ≈ N + 0.157 × z
这里的 z 是海拔高度（m），a 是地球半径。加上这一项后，标准大气下 M 随高度单调增加，电波不会被陷获。当 dM/dz < 0 时，即修正折射率随高度下降，就具备了形成波导的条件。波导高度，说白了就是 M 剖面上出现负梯度的那一层的底高和顶高。
提示：在实际工程里，我从来不看 N 剖面判断波导，只看 M 剖面。原因很简单，N 在标准大气里也是递减的，容易误判；M 剖面一旦出现“负斜率段”，波导层一目了然。
1.2 三种波导类型的高度定义
大气波导按形成机制和位置，通常分三类：表面波导、抬升波导、蒸发波导。它们的高度定义略有差别，工程上很容易混淆。
表面波导的底高贴近地面或海面，顶高一般在几十米到几百米。M 剖面上表现为从地面开始的一段负梯度层，波导高度就是这段负梯度层的顶高。抬升波导底高在地面以上，往往在几百米甚至上千米，M 剖面上表现为中间一段负梯度夹在两个正梯度段之间，此时需要同时关注波导底高和顶高。蒸发波导特指海洋近海面层内由湿度急剧递减形成的波导，底高就是海面，顶高通常在 5 到 30 米之间。
我见过不少刚入行的同事把“波导高度”默认当成波导顶高，结果做雷达覆盖评估时把盲区范围算错了。实际业务里，如果是表面波导，重点关注顶高；如果是抬升波导，波导底高决定了电磁波从哪个高度“漏”进去，顶高决定了被束缚的大致厚度，两个都要交代清楚。
2. 哪些参数在真正决定波导高度
2.1 温度梯度：逆温层是波导的“骨架”
温度对波导高度的影响，核心不在温度本身，而在温度随高度的变化率，也就是温度梯度。标准大气里气温随高度递减，大约每上升 100 米降 0.65 摄氏度。一旦出现逆温，即气温随高度增加而升高，N 的干项就会随高度快速减小，为波导形成创造条件。
逆温层的底部，基本就是波导的下边界；逆温层的顶部，往往对应波导的上边界。常见的辐射逆温在夜间晴朗风小的条件下发展最旺盛，早晨日出前后逆温最强，波导高度也最清晰。我在做陆地微波链路评估时，最怕的就是清晨时段的探测数据，因为此时的辐射逆温会让表面波导顶高冲到 200 到 400 米，而白天混合层充分发展后，波导又消失得无影无踪。
下沉逆温则是抬升波导的主要推手。副热带高压控制区域，空气大范围下沉，在几百米到一两千米高度形成稳定的下沉逆温层，这个高度区间内极易出现抬升波导。这也解释了为什么海上远距离通信异常增强的案例，大多发生在夏季副高控制的时段。
2.2 湿度梯度：对波导高度影响最敏感的参数
如果说温度梯度是波导的骨架，湿度梯度就是波导的灵魂，尤其对低层波导而言。水汽压随高度递减得越剧烈，N 的湿项下降得越快，dM/dz 越容易转负。大量实测表明，许多表面波导的形成，湿度骤减的贡献甚至超过温度逆温。
这里有一个经常被忽略的点：湿度梯度对波导高度的“塑形”能力很强。近地面层内，如果湿度从近饱和迅速下降到干燥状态，波导顶高往往被压在几十米以内，形成典型的薄波导；如果湿度的递减是缓慢持续的，则波导层可能被拉得很厚，顶高相应抬升。海上蒸发波导就是最典型的例子：海面温度决定饱和水汽压，近海面几米内湿度快速下降，波导顶高基本由这个湿度梯度决定。
我自己的经验是，做波导高度估算时，优先看湿度廓线，再看温度廓线。很多数值预报产品对低层湿度的模拟误差偏大，直接导致波导高度算偏，这也是后面要讲的常见坑之一。
2.3 气压与风：常被忽略的调节器
气压对 N 的直接影响是气压越高，干项越大。但气压随高度的递减基本稳定，不像温度和湿度那样容易出现异常梯度，所以气压本身不是产生波导的决定性参数。它更多是背景环境，影响的是波导出现的海拔区间，而不是决定负梯度是否出现。
风的作用相对微妙。一方面，适中的风可以通过动力混合把逆温层和逆湿层维持在特定高度，帮助波导结构保持完整；另一方面，风速过大时湍流混合加强，会把温湿梯度“抹平”，波导被破坏。我见过不少海上观测报告，风速从 3 米/秒升到 8 米/秒后，蒸发波导高度明显下降甚至消失，就是这个道理。
在实际工程判断中，风速是一个很好的“否决指标”：当低层风速普遍超过 10 米/秒时，即便温湿条件有利于波导形成，实际出现强波导的概率也非常低。反之，微风或静风条件下，一旦探空显示有逆温和湿度骤减，就该警惕波导的超视距传播效应了。
3. 不同场景下波导高度怎么变
3.1 海上蒸发波导的高度规律
海上蒸发波导是通信和雷达领域最关心的波导类型，因为它常年存在于全球大部分海域，高度虽低，但对海面低空目标探测和舰载通信影响极大。它的形成机制是：海面水汽饱和，湿度随高度急剧下降，使得低层 M 剖面出现负梯度。
蒸发波导高度主要受海气温差、风速和湿度差控制。海气温差越大，波导高度越高；风速适当时，蒸发波导可以维持在 10 到 20 米；风速过大会破坏近海面层的湿度梯度，波导高度被打压到 5 米以下甚至消失。还有一个容易忽略的因素是海流，暖流区域海气温差大，蒸发波导普遍偏高；冷流区域则相反。
在工作实践中，我很少用探空数据去判断蒸发波导高度，因为探空气球上升速度快，近海面层的细节根本测不准。更实用的办法是用海面气象观测数据结合蒸发波导诊断模型来计算，常用的有 PJ 模型、NPS 模型等。这些模型输入海面温度、气温、风速、湿度，输出就是蒸发波导高度，做面覆盖预测图时非常方便。
3.2 陆地表面波导与抬升波导
陆地表面波导多在夜间和清晨出现，典型场景是晴朗夜间地面辐射冷却，近地面形成逆温层，加上植被和土壤释放的水汽在近地面堆积，湿度随高度迅速下降，两个条件叠加，M 剖面负梯度段出现。波导顶高通常在 50 到 300 米之间，日出后地温升高、湍流增强，波导很快消散。
抬升波导的情况更复杂一些。它不贴地，下方还有一段正常折射区域，因此电波不会在地面反复反射，而是在波导层内“跳跃式”传播。抬升波导的高度区间受下沉逆温控制和锋面系统影响，变化幅度很大，从几百米到两千米都可能出现。这类波导对直升机、无人机和低轨卫星通信的影响比对地面微波链路更大。
我做过一个典型的案例：某年春季一次冷锋过境后，探空数据显示在 800 到 1200 米高度存在明显的抬升波导，结果当天无人机测控链路出现间歇性中断，地面站却收到了 300 公里外某调频广播的异常强信号。这就是抬升波导把电波“锁”在特定高度层里的典型表现。
3.3 季节、昼夜与下垫面的影响
波导高度有非常明显的季节和昼夜节律。夏季和初秋是波导的高发季节，尤其是副高边缘的暖湿气流区域，抬升波导高度动辄冲破 1000 米。冬季则以近地面辐射逆温形成的表面波导为主，高度普遍偏低。
昼夜变化则与太阳辐射和边界层演变高度相关。夜间到清晨，表面波导发展成熟，高度和强度都在早晨达到峰值；午后对流旺盛，边界层内混合均匀，波导基本消失。这个规律对监测窗口的选择很重要——你白天去测波导，大概率什么都测不到，所有关键数据都集中在夜间和清晨。
下垫面性质也不能忽视。海洋、沙漠、平原、城市下垫面的热力性质和湿度条件差异巨大，波导高度差异可以达到一个数量级。沙漠地区白天地表极端高温，夜间辐射冷却也极端，表面波导反而更强烈；城市下垫面由于人为热和建筑遮挡，波导结构更破碎，高度判断也更困难。
4. 实际工作中如何获取和计算波导高度
4.1 探空数据计算 M 剖面
探空数据是计算波导高度的最经典来源。常规气象探空每天释放两次气球，时间为世界时 00 时和 12 时，对应北京时间 08 时和 20 时。数据里包含各高度的气压、温度、露点，换算成水汽压后，就可以逐层计算 N 和 M。
计算步骤并不复杂：先把露点温度转换成水汽压，常用的经验公式是马格努斯公式，然后逐层代入 N 的计算公式，再加上高度贡献项得到 M。关键在插值处理——探空数据的分辨率在高空比较粗，尤其近地面层，可能只有几百米一层，直接算 M 梯度会漏掉薄波导。我的做法是先对温湿度剖面做平滑插值，按 10 米间隔重采样，再计算 M 剖面。
判断波导时，我习惯直接用数值方法扫描 M 剖面：找到 dM/dz ≤ 0 的连续段，记录起始高度和结束高度，就是波导的底高和顶高。如果这段贴地，就是表面波导；如果上方和下方都有正梯度段，就是抬升波导。这一套逻辑写成脚本处理历史探空数据，几分钟就能统计出一个站点多年的波导高度气候特征。
注意：探空数据的时间代表性有限，一天两次的采样很可能错过波导的峰值时刻。如果做精细评估，建议叠加数值再分析资料和地面自动站的连续观测数据一起判断。
4.2 再分析资料与蒸发波导诊断模型
只用探空数据做波导高度分析，时空覆盖都太有限。现在大家更常用再分析资料，比如 ERA5，它能提供三维的温湿度风场，时间分辨率到小时级，可以弥补探空的稀疏采样问题。但从再分析资料算波导高度有一个陷阱：它的低层湿度在海洋上误差偏大，直接算出来的蒸发波导高度经常系统性偏低或偏高。
我的处理办法是，用再分析资料判断波导的大尺度分布和演变趋势，用实测探空和气象站数据进行校准。对于海洋蒸发波导，则用诊断模型单独计算。PJ 模型相对保守，适合中性层结；NPS 模型对不稳定层结的响应更灵敏，在中低纬度海域表现更好。无论用哪个模型，输入数据的质量直接影响结果，海面温度和湿度务必用高质量的浮标或船舶观测数据。
在实际项目中，我通常把波导高度预测做成一张随时间变化的廓线图，叠加不同高度上的反射路径示意。这张图在电磁兼容分析、干扰排查和雷达覆盖评估中非常实用，比单纯报告“今天存在波导”有价值得多。
5. 常见问题与排查技巧实录
5.1 几个典型问题
这些年处理了不少波导相关的工程问题，挑几个高频问题整理成表，供大家排查参考。
| 问题现象 | 可能原因 | 排查方向 |
|---|---|---|
| 雷达突然出现超视距探测回波 | 表面波导或抬升波导形成 | 查看探空 M 剖面、海杂波特征 |
| 微波链路清晨误码率飙升 | 辐射逆温导致表面波导顶高变化 | 对比 08 时探空与链路余量数据 |
| 海上通信距离异常增大 | 蒸发波导高度抬升 | 用海面气象数据跑蒸发波导模型 |
| 无人机测控信号间歇中断 | 抬升波导把信号困在特定层 | 关注 500 米以上逆温和逆湿层 |
| 干扰信号来向无法定位 | 波导带来多路径和超视距干扰 | 结合波导高度做射线追踪分析 |
5.2 我的几条实操建议
第一，永远不要只看单一高度的气象数据来判断波导。无论是地面站还是船舶观测，只能反映一个高度上的状态，波导判断需要完整的垂直廓线。没有探空数据的时候，至少结合再分析资料看温湿度剖面。
第二，关注探空仪的湿度响应延迟。湿度传感器在低温环境下响应慢，上升过程中测到的湿度廓线会偏“平滑”，计算出的波导高度往往偏高。遇到明显不合理的波导层，先怀疑湿度数据的质量。
第三，多和频谱监测人员联动。波导来临时，频谱仪上通常有大量远距离广播和通信信号同时出现，这是最直接的旁证。有一次我们怀疑某个波导高度判断不准确，后来通过频谱监测到的异常信号频率反推传播路径，验证了波导顶高确实比模型算出的高了几十米。
第四，做波导高度预测时，务必区分“能形成波导”和“实际形成强波导”。前者只需 dM/dz 转负，后者还需要负梯度段足够厚、强度足够大，才能有效陷获电波。薄得可怜的负梯度层，对实际传播的影响可以忽略。
最后，如果你需要长期做波导高度分析，建议自己积累一个本地化修正系数库。每个地区的下垫面、气候背景不同，再分析资料和模型的误差特征也不一样，把本地实测数据与模型结果对比积累几个月，后续预测精度会明显提升。这个活没有捷径，但做完了真的香。
本文还有配套的精品资源，点击获取
