简介：在PSCAD/EMTDC环境下开展雷电波模拟的专题资料，面向电力系统设计、过电压防护与电磁暂态分析方向的工程师、科研人员及高校师生。文档以GB18802标准中浪涌保护器（SPD）的I级（50kA冲击电流，1.2/50波形）与II级（20kA冲击电流，8/20波形）分类试验为背景，系统梳理雷电波双指数函数表达式，并结合PSCAD模型设置详细解释波前时间与半峰值时间的负值换算原理，帮助读者规避常见参数配置错误。针对II级波形拟合难度较大的情况，文档还给出反复仿真、逐步调整参数以获得高吻合度结果的实操思路，并延伸到直击雷、感应雷等典型雷电波的构建方法。整包仅1个doc文件，大小约120KB，内含仿真图和配套分析，可对照文档说明逐项复现模型并验证波形输出。该资源已有2006人学习，适合希望快速掌握雷电波建模方法、深入理解SPD试验波形机理并优化仿真准确度的读者。
1. 雷电波模拟在PSCAD里难的不是建模，而是让波形和边界条件都对
不少新接触PSCAD建模仿真的人，拿到“雷电波的模拟”这个题目后，第一反应是找元件库里有没有现成的雷电源，拖出来接上线就运行。真正跑起来才发现，仿真图画是画出来了，可讲不出每个折点对应什么物理过程，换一组参数波形就面目全非。雷电波模拟和工频短路仿真最大的区别在于时间尺度：工频仿真看的是毫秒级的变化，雷电波看的却是微秒级的波头与反射过程，步长、线路模型、接地阻抗任何一个环节选错，结果都会偏离工程实际。这篇笔记会把雷电波模拟从“波形怎么选”讲到“仿真图怎么读”，再落到PSCAD里的建模步骤、参数设置和几次血泪踩坑，适合正在做输电线路防雷设计、设备绝缘配合，以及用PSCAD做毕业设计高压研究的人直接参考。
2. 雷电流波形与等效电路：先把几个关键物理量定下来
2.1 雷电流波形“长什么样”才值得仿真
雷电波的工程等效波形有很多种，但在PSCAD这类电磁暂态仿真里最常用的还是双指数波。双指数波的表达式可以写成：
I(t) = I0 × ( e^(-αt) - e^(-βt) )
其中 I0 是雷电流峰值，α 控制波尾衰减快慢，β 控制波头上升快慢。常见的参数组合有 2.6/50 μs 和 1.2/50 μs 两套：斜杠前的 2.6 μs 是波头时间，也就是电流从零上升到峰值的时间；斜杠后的 50 μs 是波尾时间，描述峰值衰减到一半的时间。需要提醒的是，2.6/50 μs 通常用于雷电流波形，1.2/50 μs 更多用于雷电过电压波形，做线路防雷分析时不要混用。
在PSCAD里把双指数波变成实际电路，通常不是直接写这个公式，而是用一个受控电流源，把波形信号作为控制输入。受控电流源本质上是个理想的电流注入源，雷电击中导线时表现为向线路注入电流，所以用电流源比电压源更贴近物理过程。如果你用了电压源，波形会带上电源内阻的约束，实测的雷电流陡波头特征很容易被削平。
模拟雷电波的电源在PSCAD里有三种常见实现方式：
| 实现方式 | 精度 | 操作复杂度 | 适用场景 |
|---|---|---|---|
| 元件库自带的雷电源模块 | 中 | 低 | 快速摸底、论文示意、初学练习 |
| 受控电流源 + 波形控制信号 | 较高 | 中 | 参数可调范围大，工程分析常用 |
| 外部实测波形数据文件导入 | 高 | 较高 | 需要复现现场实测数据的专项研究 |
如果你所在版本的元件库里找不到雷电源模块，不用纠结，用受控电流源加波形控制信号的方式会更可靠，参数全在自己手里。
2.2 回击通道与接地回路：先决定等效电路
雷电波仿真里最容易忽略的是回击通道的等效阻抗。雷电回击通道在电磁暂态分析中通常被简化成一个波阻抗，常见取值为 300Ω 到 1500Ω 之间，建模时有人干脆忽略它。忽略通道阻抗意味着雷电流源直接并联在导线上，所有电流全部注入线路，这会让线路上的过电压幅值偏大，偏保守，但不符合实际分流关系。做设备绝缘配合时这类保守可以接受，做雷击跳闸率评估时就会明显失真。
另一个关键参数是接地电阻。雷电流流过接地装置时，由于土壤放电和冲击特性的影响，实际呈现的冲击接地电阻往往小于工频接地电阻。常见做法是直接用线路杆塔的工频接地电阻 10Ω 或 15Ω 参与仿真，这往往会抬高接地点的电位，属于偏保守的估算。如果想更贴近实际，可以考虑在PSCAD里用非线性电阻近似冲击接地电阻，或者至少把工频接地电阻和冲击接地电阻的差异写进分析报告的假设说明里。
2.3 三种雷电流源实现：按仿真目标选型
先说明选择原则：如果你的目的是课程设计或方案论证，元件库自带的雷电模块足够用，省时间；如果你要改变波头时间、波尾时间、峰值，还要研究幅值变化对过电压的影响，就选受控电流源加控制信号；如果你有实测的雷电流示波图数据，或者想做更严谨的复现，用数据文件导入是最佳路线。
受控电流源加波形信号的具体做法是：在PSCAD里放置一个受控电流源，控制端口的输入信号来自控制系统的波形生成模块。波形信号可以用分段线性描述来近似双指数波，例如在 0 到 2.6 μs 之间用一条快速上升的折线，之后按指数衰减分成若干线段。分段越细，波形越接近真实的双指数波，但仿真步长也要相应缩小。这种方式的好处是每个时间点的电流值都清楚，波形哪里不对可以直接看到控制信号的曲线。
外部数据文件导入方式适合对波形真实性有要求的场景。PSCAD支持从文本文件读取数据序列作为信号源，文件里两列数据，第一列是时间，第二列是电流值。实测数据往往带噪声，导入前需要做平滑和去毛刺处理，否则高频噪声会干扰仿真收敛性。
3. PSCAD里搭一个雷电波模型：元件选择、参数设置、运行出图
3.1 从新建工程到第一个模型：关键参数与顺序
先把模型的结构讲清楚。一套最简单的雷电波模拟模型包含四个部分：雷电流源、被击中的输电线路、线路末端的设备等效阻抗、测量与输出通道。我一般会在PSCAD里按这个顺序建：
- 新建工程，先改仿真时间参数。
- 放置雷电流源，接在雷击点的节点上，电流源另一端接地。
- 放置输电线路模型，一端接雷击点，另一端接负载。
- 放置电压表和电流表，设置输出通道。
- 运行仿真，打开波形窗口。
仿真参数设置是关键的第一步。雷电波波头只有几微秒，主仿真步长必须足够小才抓得住陡波头。一般来说主步长不要超过 0.1 μs，我习惯先设 0.05 μs，仿真总时长用 40 μs 到 60 μs。60 μs 可以看到完整的波尾衰减和至少一两次的末端反射。如果把主步长设成工频仿真习惯的 10 μs 或 20 μs，波形直接变成一条圆滑的曲线，波头细节全部丢失。
具体参数可以参考下面这张表：
| 参数项 | 建议值 | 说明 |
|---|---|---|
| 仿真总时长 | 40 ~ 60 μs | 覆盖波头和波尾完整过程 |
| 主仿真步长 | 0.05 μs | 至少保证波头有 50 个采样点 |
| 绘图步长 | 0.01 ~ 0.02 μs | 画图精度更高，波形更平滑 |
| 雷电流峰值 | 10 ~ 20 kA | 典型线路雷击电流幅值 |
| 波头时间 | 2.6 μs | 按 2.6/50 μs 标准波形 |
| 波尾时间 | 50 μs | 同上 |
PSCAD里主仿真步长是整个电磁暂态求解的推进步长，绘图步长只影响波形输出和文件大小。绘图步长如果比主步长大，波形会丢失细节；比主步长小，软件会做插值输出。雷电波这种陡波前信号，两个步长设置得越接近越好。
3.2 线路模型：Bergeron 与频率相关模型怎么选
输电线路在PSCAD里有多种模型，雷电波模拟中常用的有Bergeron模型和频率相关模型。Bergeron模型把线路等效成波阻抗和两端反射网络，适合分析行波传播和反射过程，计算速度快，参数直观，适合初版仿真。频率相关模型能反映雷电波高频分量沿线路衰减更快的特性，精度更高，但参数多、计算慢，而且需要准确的线路几何参数。
我的选择建议是：第一版模型用Bergeron，先把波形趋势跑出来，如果分析重点是过电压幅值的准确性，再换成频率相关模型对比一次。两个模型在末端开路时的反射电压会有明显差异，这正是线路模型对高频分量处理不同造成的。
Bergeron模型的参数主要有线路长度、波阻抗、波速。架空线路的波速可以近似取 300 m/μs，线路长度按实际档距设置。比如线路长度取 1 km，雷电波从雷击点传播到末端大约需要 3.3 μs，这个数值在后面的仿真图分析里可以用来验证模型是否正确。
线路末端设备一般用等效阻抗表示。如果是变压器入口，可以简化成一个小电容并联高阻，电容值在几百 pF 到几千 pF 之间；如果是避雷器保护，需要单独放置避雷器模型，避雷器的参考电压和残压按所选型号参数填。注意末端等效元件设置不当会造成非物理的电压振荡，这一点在避坑章节细说。
3.3 测量与输出通道：做出一张能放进报告里的仿真图
测量环节直接影响仿真图的质量。电压表要并联在被测节点和地之间，电流表串联在支路里。PSCAD里默认的输出通道单位可能是 p.u.，这一点必须改掉。雷电过电压分析要用绝对单位 kV，雷电流用 kA，否则仿真图上的数值没有物理意义，写报告时自己都容易算错。
输出通道命名也要注意。多个通道混在一起，波形窗口里靠颜色根本分不清，最好在通道名称里写清楚节点和量纲，例如“母线入口过电压_kV”“雷电流源输出_kA”。波形输出后，横轴时间单位显示为 s，需要手动调整成 μs，仿真图更直观，也符合雷电波分析的习惯。
PSCAD波形窗口支持直接导出图片，导出时勾选“包含坐标轴名称和刻度”，这样图片放进Word后不需要再补坐标说明。波形比较密集时，可以用缩放工具局部分析波头段，报告中可以放一张全过程的波形图，再加一张波头局部放大图，这是评审最容易认可的形式。
4. 仿真图怎么读：从波形反推模型对不对
4.1 源头波形检查：先确认雷电流源输出“像不像”
仿真出来的第一张图，应该先看雷电流源本身的输出波形。正常的双指数波是迅速上升到峰值，然后按指数规律衰减，整个过程平滑无毛刺。如果源波形上升沿不是一次到位，而是出现明显台阶，通常是波形信号的分段线性近似点太少。解决办法是增加分段数，或者在上升段单独加密时间点。
还要看波形起始时刻。雷电流源如果带有初始延时设置，波形会晚于 0 时刻出现。延时本身不影响分析，但读传播时间时必须减去这段延时，否则算出来的波速对不上。把源波形和线路末端波形放在同一个窗口对比，能看到明显的传播延迟，这个延迟有没有落在预期范围内，是判断整条线路参数是否正确的第一把尺子。
雷电流源输出波形的幅值也需要核对。受控电流源的输入信号如果单位不匹配，波形可能差好几个数量级。PSCAD的控制信号通常按 p.u. 处理，受控电流源内部有比例系数，检查这个系数是否设成了 1 kA / p.u.，如果没有，记得换算。
4.2 沿线传播与末端反射：读图的边界条件
雷电波沿线路传播到末端时，会发生折射和反射，这是仿真图里最值得分析的部分。末端开路时，反射系数接近 1，末端电压会接近入射波的两倍；末端接有低阻负载时，反射系数为负，电压可能不升反降。很多初学PSCAD的人看到末端电压超过源电压就以为模型发散，其实是边界条件没有被理解。
我在分析仿真图时通常会做一张表，把不同末端条件下的理论预期写下来，再和仿真结果对比：
| 末端类型 | 理论预期 | 仿真图上对应特征 |
|---|---|---|
| 开路 | 电压约翻倍，电流为零 | 末端电压峰值约为入射波两倍，电流波形归零 |
| 阻性负载匹配 | 无反射，电压幅值约等于入射波 | 末端电压和源电压波形几乎一致 |
| 变压器入口电容 | 电压上升后有小幅振荡 | 波头后出现高频衰减振荡 |
| 避雷器动作 | 电压被限制在残压附近 | 电压波形顶部出现明显削平 |
读了这张表再回到仿真图上，你就能判断末端等效阻抗设得对不对。比如末端接了电容，波形没有任何振荡，往往说明电容值太小或者被并联的大电阻完全旁路了。
4.3 FFT 和传播时间验证：把“看起来对”变成“数值对”
仿真图“看起来差不多”是不够的，还需要做两个数值层面的验证。第一个是传播时间验证：线路长度除以波速得到理论传播时间，在仿真图上量出入射波到达末端的时刻，两者误差应该在仿真步长的几倍以内。如果误差大，先检查线路长度单位，PSCAD里长度单位有 km 和 m 之分，填错会导致传播时间差 1000 倍。
第二个是频谱验证。雷电波的等效频率范围和波头时间有关，波头越陡，高频分量越丰富。可以把末端电压波形做 FFT，观察主要频率分量是否在几十 kHz 到几百 kHz 之间。如果 FFT 结果显示主要能量集中在上 MHz 甚至更高，基本可以断定仿真步长过大，波形里混入了数值振荡。这种振荡在时域图上是锯齿状毛刺，在频域图上就是高频尾巴。
有些波形看起来平滑但数值偏高，用 FFT 能看出原因：高频衰减不足。频率相关模型和Bergeron模型在这个问题上的表现差异很大，建议在分析阶段至少要跑一次频率相关模型做对比，否则你无法判断当前结果是因为模型简化还是真实物理过程。波形里的振荡到底是什么，靠肉眼常分不清，在PSCAD里用FFT工具一验就露馅，这不是玄学，是边界条件没理清。
5. PSCAD雷电波仿真的避坑记录：常见问题与排查方法
5.1 仿真波形变成光滑曲线，看不到波头尖峰
现象：雷电流源设置的是 2.6 μs 波头，但仿真图上的电流上升沿被拉得很平，波头细节完全消失，整个波形像被低通滤波过一样。
原因：最典型的原因是主步长或绘图步长设置过大。我见过把主步长设成 10 μs 跑雷电波仿真的，2.6 μs 的波头连一个采样点都分不到，波形自然被抹平。绘图步长过大即使主步长合理，也会出现同样问题。
解决：把主步长调到 0.05 μs，绘图步长调到 0.01 μs，重新运行。如果波形仍然平滑，再去检查受控电流源的控制信号本身，看波形生成模块里的时间轴单位是不是毫秒，单位写错会让整个波形时间尺度偏大三个数量级。
5.2 雷电源接成了串联元件，波形到处是毛刺
现象：仿真电流波形出现连续的尖锐振荡，波形整体形态完全不像双指数波，电压波形在某些时刻出现非物理的突变。
原因：雷电流源被串联进了主回路。雷电波模拟的电流源是并联使用的，雷电流源一端接导线节点，另一端接地，相当于一个向节点注入电流的并联支路。串联进去后，电源内阻和线路阻抗直接分压，波形彻底失真。
解决：把电流源改回并联接法，确认电流源的一端在雷击节点上，另一端直接接地。查看受控电流源元件方向，PSCAD里带方向的电流源如果接反了，波形会反向，这也是毛刺来源之一。
5.3 末端开路电压翻倍，被误判为模型发散
现象：线路末端电压峰值是入射波的两倍多，波形出现明显的台阶状上升，有人担心是数值发散。
原因：这是传输线末端开路时的正常反射现象。行波到达开路末端时反射系数接近 1，反射波和入射波叠加，电压接近翻倍。如果完全没有翻倍，反倒说明末端等效负载设置得有问题。
解决：先不要急着改模型，用折射反射公式手算一遍理论峰值。确认末端确实是开路，确认波阻抗参数和线路长度与实际一致，确认幅值落在理论值范围内，然后继续分析。只有峰值超过理论值两倍以上才需要怀疑数值稳定。
5.4 忽略冲击接地电阻，过电压结果整体偏高
现象：仿真得到的杆塔地电位升很高，绝缘子串两端电压明显超过实际运行经验值，换不同接地电阻值时结果波动剧烈。
原因：模型中使用了工频接地电阻作为接地阻抗，而雷电流作用下接地体的冲击阻抗往往低于工频电阻，高频电流下还有电感效应和火花效应，差距可能达到 2 倍以上。
解决：如果手里有现场接地装置的冲击接地电阻实测值，直接填入模型；没有实测值就在报告里明确说明用的是工频电阻的保守假设。参数敏感性分析可以给接地电阻设几组数据，10Ω、5Ω、2Ω 各跑一次，看过电压对接地电阻的依赖程度。
5.5 改变线路长度，波形传播时间却没有变化
现象：把线路长度从 500 m 改成 2 km，末端电压波形到达时刻没有明显变化，波速验证对不上。
原因：线路长度参数没有真正传给线路模型。常见情况是线路元件内部还有一个独立的总长度参数，外部输入的长度没有关联过去，或者线路模型被设置成了集中参数模式，不再按分布参数计算传播延时。
解决：在线路元件内部参数中确认使用分布参数模型，并把长度参数关联到外部输入变量。改完长度后重新运行，用波速 300 m/μs 验证传播时间是否线性变化，这一步能筛出多数线路参数设置错误的场景。
6. 进阶：从标准波形走向实测波形，做一个能说服人的对比仿真
如果只停留在双指数波的阶段，仿真结果只能说明“雷电波大概长什么样”。要做更扎实的结论，我建议在模型稳定后做两件事：一是导入实测雷电流波形做对比仿真，二是叠加避雷器和绝缘子闪络判据，让仿真图里反映出真正的保护动作过程。
导入实测波形的方式不复杂：把雷电流示波数据整理成两列文本文件，第一列时间、第二列电流，在PSCAD的数据读取模块中加载。实测波形往往有噪声和测量毛刺，入模之前先做平滑，不然这些毛刺经过传输线模型后会被放大成高频振荡。同一套模型跑两次，一次用标准双指数波，一次用实测波形，两张仿真图放在一起对比过电压峰值和波形的差异，这种对比在报告里非常有说服力，也能帮助你判断采用标准波形估算是否偏危险。
避雷器和绝缘子闪络的加入会让仿真图的分析层次更丰富。避雷器模型设好参考电压和残压参数后，可以在图上看到原本持续上升的过电压波形被削平在残压附近；绝缘子闪络判据可以按过电压达到闪络电压时触发导通来处理。有了这两个环节，仿真图就不只是“雷电流注入后的电压响应”，而是一整套“雷击→波传播→保护动作→设备承压”的完整过程。
我自己的习惯是保留第一版最简单的模型作为备案，每次加复杂元件之前复制一个工程副本。这个过程我翻过不少次车，有加避雷器后不收敛的，有闪络判据触发条件写错导致波形乱跳的，但都不亏——因为每一次翻车都把模型参数理解得更清楚。做仿真的经验感和直觉就是这样堆出来的，模型越复杂，留的退路越多越好。这套从双指数波到实测波形的对比方法，已经帮我处理过好几轮雷击过电压分析的评审意见，希望帮到你。
本文还有配套的精品资源，点击获取
