PSCAD双指数雷电波建模：1.2/50与8/20波形参数详解与避坑指南

简介：面向电力系统工程师、电气工程专业学生及相关从业者的PSCAD雷电波模拟建模资料，围绕雷击过电压导致电气设备损坏的问题，聚焦浪涌保护器分类试验。内容覆盖I级50kA 1.2/50冲击电流与II级20kA 8/20冲击波形的数学表达式、建模参数设置，并配有仿真图和逐步分析，梳理了从波形公式到模型设定的完整流程，便于理解雷电波传播过程，评估系统耐雷能力并设计防雷措施。资源包仅含1个doc文档，约120KB，内容集中，适合快速查阅与对照实操。目前已有2005人学习下载，可作为课程设计、论文参考或工程方案比选时的速查资料。文档还针对PSCAD中时间常数以秒为单位的特点，说明了1.2/50波形中波前时间采用负值833333.33倍、半峰值时间取负值1/20000等易错设置，并结合指数衰减公式展示波形拟合思路，可帮助读者避开建模误区、提升冲击电流仿真准确度，尤其适合电磁暂态仿真的初学者快速上手。

1. 雷击波建模的起点：先从 1.2/50 和 8/20 两个波形说起

做 PSCAD 建模的雷电波模拟，第一步往往不是搭线路、挂避雷器，而是先把“雷电流波形”本身复现对。前阵子帮一个朋友复盘某配电房避雷器解体事故，翻记录时发现当初选型用的冲击波形是 8/20，而实际雷击波的波前只有不到 2 微秒，电流上升率差了好几倍。这件事让我重新翻了一遍 GB18802 里 I 级、II 级分类试验的冲击波形定义，也在 PSCAD 里把 50kA 1.2/50 和 20kA 8/20 两条双指数雷电波完整地建了出来。这篇笔记就围绕这两个波形，讲清双指数数学模型、PSCAD 参数换算、受控源搭建，以及拟合时最容易踩的坑，适合刚开始做防雷仿真、又想确认自己模型参数没填错的人。

2. 双指数雷电波模型：1.2/50 与 8/20 的数学拆解与参数换算

2.1 雷电波为什么能用双指数函数描述

真实雷电流最明显的特征就是“快上慢下”：波前阶段电流在微秒级冲到峰值，波尾阶段再缓慢衰减到零。工程上为了统一描述这种形态，引入了视在波前时间 T1 和半峰值时间 T2 两个指标。所谓 1.2/50，指的就是视在波前时间为 1.2 微秒、从起始点到波形下降到半峰值的时间为 50 微秒；8/20 则是波前 8 微秒、半峰值时间 20 微秒。雷电流的现场实测波形千变万化，但做设备型式试验和仿真校验时，必须用这种标准化的波形代替，否则不同厂家、不同项目之间的数据完全没法对比。

双指数函数 I(t) = A × [exp(-t/τ_tail) - exp(-t/τ_rise)] 恰好能同时表达“快速上升”和“缓慢衰减”两个特征。t = 0 时两项相消，电流为零；t 稍微增大，衰减更快的 τ_rise 项迅速退去，波形出现峰值；之后主要由 τ_tail 项主导，缓慢下降。这个数学形态和雷电流的物理过程高度相似，所以 PSCAD、ATP-EMTP 这类电磁暂态软件普遍采用双指数模型作为标准冲击源的基础表达式。

需要特别注意的是：控制波前和控制波尾的两个时间常数并不相同，τ_rise 决定了上升段的陡峭程度，τ_tail 决定了波尾衰减快慢。模型里哪个填错，波形就会变成“慢上快下”甚至完全发散。双指数函数还有一个容易被忽略的性质——它算出来的峰值并不等于幅值系数 A，因为上升项退去时衰减项已经开始起作用，两者叠加后峰值会出现明显跌落。这个跌落幅度在不同波形参数下差异很大，是后续建模调参最需要关注的问题。

2.2 GB18802 试验波形的关键参数：从微秒到秒的换算

GB18802 中 I 级分类试验对应 50kA、1.2/50 冲击电流，II 级分类试验对应 20kA、8/20 冲击电流。PSCAD 里所有时间相关参数都按“秒”处理，所以第一步是把微秒换算成秒，再换算成双指数函数里的时间常数。

对于 1.2/50 波形，τ_rise = 1.2e-6 秒，τ_tail = 50e-6 秒。如果写成 I = A × [exp(-α·t) - exp(-β·t)] 这种系数形式，那么 α = 1/τ_tail = 1/50e-6 = 20000，β = 1/τ_rise = 1/1.2e-6 ≈ 833333.33。这正好对应了资料里出现的两个关键数字：833333.33 来自 1/1.2e-6，1/20000 就是半峰值时间 50 微秒的另一种写法。

同样地，8/20 波形对应 τ_rise = 8e-6 秒，τ_tail = 20e-6 秒。按系数形式可以写成 I = A × [exp(-50000·t) - exp(-125000·t)]，其中 50000 = 1/20e-6，125000 = 1/8e-6。

这里要特别提醒一个常见的混淆点：同样一个波形，有的资料用时间常数 τ 填写，有的资料用系数 1/τ 填写，两者是倒数关系，差着数量级。比如 1.2/50 的波前，如果 PSCAD 组件要求填“时间常数”，就填 1.2e-6；如果要求填“指数系数”，就填 833333.33。填错后波形要么快得看不见，要么慢得像工频。资料里将 833333.33 和 1/20000 混在一起写，实际就是因为他们一个用的是系数形式，一个用的是时间常数形式。

2.3 用 NumPy 把两条波形先画出来

在进 PSCAD 之前，我习惯先用 Python 把两条理论波形算一遍，确认参数方向和幅值跌落情况。下面这段脚本用 NumPy 实现双指数波形，并自动计算峰值与实际时间参数。

import numpy as np def lightning_current(amp, tau_rise, tau_tail, t): # amp: 幅值系数（kA） # tau_rise: 波前时间常数（s） # tau_tail: 半峰值时间常数（s） # t: 时间数组（s） return amp * (np.exp(-t / tau_tail) - np.exp(-t / tau_rise)) t = np.linspace(0, 200e-6, 40000) # 1.2/50 波形，幅值系数先给 50kA，观察实际峰值 i_12_50 = lightning_current(50, 1.2e-6, 50e-6, t) # 8/20 波形，幅值系数先给 20kA，观察实际峰值 i_8_20 = lightning_current(20, 8e-6, 20e-6, t) def analyze_waveform(t, i): peak_idx = int(np.argmax(i)) i_peak = i[peak_idx] t_peak = t[peak_idx] # 10% 和 90% 峰值点 i_10 = 0.1 * i_peak i_90 = 0.9 * i_peak idx_10 = int(np.argmin(np.abs(i[:peak_idx] - i_10))) idx_90 = int(np.argmin(np.abs(i[:peak_idx] - i_90))) t_10 = t[idx_10] t_90 = t[idx_90] # 视在波前时间按 10%~90% 外推 t_front = (t_90 - t_10) / 0.8 # 半峰值时间：从 0 到下降到 50% 峰值的时间 half_tail = 0.5 * i_peak idx_half = int(np.argmin(np.abs(i[peak_idx:] - half_tail))) + peak_idx t_tail = t[idx_half] return i_peak, t_front, t_tail res_12_50 = analyze_waveform(t, i_12_50) res_8_20 = analyze_waveform(t, i_8_20) print(f"1.2/50: 实际峰值={res_12_50[0]:.2f}kA, 波前时间={res_12_50[1]*1e6:.2f}us, 半峰值时间={res_12_50[2]*1e6:.2f}us") print(f"8/20: 实际峰值={res_8_20[0]:.2f}kA, 波前时间={res_8_20[1]*1e6:.2f}us, 半峰值时间={res_8_20[2]*1e6:.2f}us")

这段脚本的逻辑很简单：先用lightning_current生成波形，再用analyze_waveform从数值上提取三个关键指标。argmax找到峰值索引，峰值之前找 10% 和 90% 点，外推得到视在波前时间；峰值之后找 50% 点，得到半峰值时间。参数说明放在函数注释里，实际运行时注意 t 的单位必须是秒，tau_rise和tau_tail也要用科学计数法。

跑完你会发现两组结果都对不上目标值：1.2/50 用 50kA 幅值系数算出来峰值只有约 44.5kA，8/20 用 20kA 幅值系数算出来峰值只有约 6.5kA。这就是我在 2.1 里提到的“峰值跌落”。8/20 跌落尤其严重，因为它的两个时间常数相差只有 2.5 倍，上升项还没消失衰减项就追上来了。所以进 PSCAD 之前，脑海里要先有这样一个概念：幅值系数 A 不等于目标峰值，需要单独做峰值校正。

3. PSCAD 搭建 I 级试验雷电波：50kA 1.2/50 的受控源建模

3.1 PSCAD 里实现双指数雷电波的三种路径

在 PSCAD 里实现双指数雷电波，常见做法有三种。第一种是直接用软件自带的冲击源组件，好处是参数卡直观，不用写代码；坏处是不同版本对时间常数的填写规则不统一，有的要求填时间常数，有的要求填系数，遇到版本升级容易翻车。第二种是用受控电流源加 Expression 表达式，直接在信号框里写数学公式，灵活但公式长了不便于复用，改动参数要进表达式逐项改。第三种是我最常用的方式：用 CSMF 自定义 Fortran 模型，把双指数公式封装成一个带参数输入的组件，项目里多个测点、多条线路需要重复使用时，拖出来改参数就行。

三种方式各有适用场景：快速验证用内置源，临时算一版用表达式，正式项目中要做多工况对比就老老实实建 CSMF。下面这张表可以帮你快速决策：

3.2 用 CSMF 写双指数雷电波模型

CSMF 在 PSCAD 里可以直接挂到 Master Library 中，右键新建组件后选择 CSMF，把下面的 Fortran 逻辑写进 Script 窗口。这个模型接收三个输入参数：幅值系数、波前时间常数、半峰值时间常数，输出一路电流信号。实际使用时，把这路信号接到受控电流源上，再串联到被测回路即可。

SUBROUTINE R_LIGHTNING(IAMP, TAU_RISE, TAU_TAIL, TIME, OUT) C C 双指数雷电波源模型 C IAMP: 幅值系数（kA） C TAU_RISE: 波前时间常数（s），对应 1.2e-6 或 8e-6 C TAU_TAIL: 半峰值时间常数（s），对应 50e-6 或 20e-6 C TIME: PSCAD 仿真时间变量（s） C OUT: 输出电流瞬时值（kA） C REAL IAMP, TAU_RISE, TAU_TAIL REAL TIME, OUT OUT = IAMP * (EXP(-TIME / TAU_TAIL) - EXP(-TIME / TAU_RISE)) RETURN END

这段代码的核心只有一行：OUT = IAMP * (EXP(-TIME / TAU_TAIL) - EXP(-TIME / TAU_RISE))。TIME 是 PSCAD 仿真时间变量，单位是秒，不需要自己定义。逻辑上先算衰减项EXP(-TIME / TAU_TAIL)，再算上升项EXP(-TIME / TAU_RISE)，两项相减得到完整的双指数波形。参数方面，IAMP 不要一上来就填目标峰值，1.2/50 波形的跌落系数约 0.89，想得到 50kA 峰值，IAMP 可以先预填 56.2kA，具体值建议跑一版后再校准。

3.3 仿真配置与波形观测：步长、时长与标定

模型建好之后，仿真配置直接影响波形质量。针对 1.2/50 这种波前只有 1.2 微秒的冲击波，仿真步长必须足够小，否则上升沿会变成一截折线。我一般把仿真总时长设为 120 微秒到 200 微秒，既能完整看到半峰值点又不浪费计算时间；步长设 0.01 微秒到 0.05 微秒，至少比波前时间小一个数量级。时间步长太大会出现两个问题：一是峰值点可能恰好落在采样间隔之间，导致仿真峰值偏小；二是波形上升沿出现台阶，后续测量波前时间时误差放大。

接好模型后，在输出通道上加一个电压/电流测量元件，运行仿真后用游标读取峰值、10% 和 90% 点。第一次跑出来的峰值大概率不是 50kA，这时不需要重新调 IAMP，在输出侧串联一个 Gain 模块即可。Gain 设置为 50 除以实测峰值，比如实测 44.5kA，Gain 就填 1.124。这样做的优势是不动源模型，保持双指数函数的数学形态不变，只做输出幅值标定。标定完再检查波形上 10% 到 90% 的上升时间是否接近 1.08 微秒（1.2 微秒的 0.9 倍），半峰值点是否落在 50 微秒附近，两个指标都满足说明模型搭建正确。

4. II 级试验 20kA 8/20 波形：拟合难点与调参路径

4.1 8/20 为什么比 1.2/50 难拟合

8/20 波形试验在防雷工程中很常见，但它的拟合难度明显高于 1.2/50。原因在于两个时间常数太接近：波前 8 微秒，半峰 20 微秒，比值只有 2.5，而 1.2/50 的比值大于 40。时间常数接近意味着上升项还没来得及完全衰减，衰减项已经明显拉低整体幅值，结果就是“峰值跌落”极其严重。

用实际数字说明：把 8/20 代入双指数公式，幅值系数 A = 20kA 时，实际峰值只有约 6.5kA，连目标值的三分之一都不到。如果像 1.2/50 那样简单放大到 31kA 左右，虽然峰值能怼到 20kA，但波形整体形态会变得很平缓，波前时间可能已经漂到 10 微秒开外。这就是资料里说“拟合这个曲线有点难度”的根本原因——三个参数互相牵制，动一个另外两个跟着偏。

4.2 两条可行路线：双指数修正与分段线性插值

面对 8/20，我通常先试双指数修正路线：保持双指数函数形式，但用预偏置幅值加参数微调来逼近目标波形。幅值系数先按 61kA 起跳，这是基于跌落系数 0.326 估算的初值。跑完一版后看三个指标：峰值不够就增大 IAMP；半峰时间偏长就把 TAU_TAIL 调小；波前时间偏慢就把 TAU_RISE 调小。每次只动一个参数，观察对应指标的变化方向。

如果双指数怎么调都差一口气，就换第二条路线：分段线性插值。用 Python 把目标 8/20 波形离散成 20 到 30 个点，导出成 PSCAD 可以读取的数据文件，再用线性插值信号源加载。这种方法不依赖双指数函数的数学拟合，只要离散点足够密，波形精度几乎不受限制。

import numpy as np # 用修正幅值生成离散点，目标是让峰值接近 20kA t_points = np.linspace(0, 100e-6, 30) i_points = 61.3 * (np.exp(-t_points / 20e-6) - np.exp(-t_points / 8e-6)) # 归一化，确保最大值为目标峰值 20kA i_points = i_points / i_points.max() * 20.0 # 打印成 PSCAD 数据文件点表格式 for t_val, i_val in zip(t_points, i_points): print(f"{t_val:.8e}, {i_val:.6f}")

这段脚本的核心是生成一组坐标点，供 PSCAD 分段线性源使用。t_points在 0 到 100 微秒内均匀取 30 个点，i_points先按双指数公式计算，再归一化到 20kA 峰值。打印结果可以复制保存为文本文件。参数上注意：点太密文件会大且可能出现抖动，点太疏波前不够精确，30 个点对 8/20 波形是经验上比较稳的选择。

4.3 调参顺序：先半峰、再波前、最后校准幅值

不论走哪条路线，调参顺序都必须固定下来，否则会陷入“动 A 影响 B，改 B 又带偏 C”的循环。我的习惯是先调半峰值时间，再调波前时间，最后校幅值。

半峰值时间主要由 TAU_TAIL 控制。目标 20 微秒，初始 TAU_TAIL = 20e-6，如果仿真测出来半峰时间是 24 微秒，就往小调 10% 左右。波前时间主要由 TAU_RISE 控制，目标 8 微秒，注意波前时间不等于 TAU_RISE，两者之间有约 1.1 到 1.2 倍的关系，所以填 7 微秒左右通常能测出接近 8 微秒的波前。幅值校准放在最后，因为在前面两步调整过程中波形峰值本身也在变化，过早校准等于白做。

这套顺序适用于固定波形形态的逼近过程，做完一轮如果误差还在 20% 以上，就回到 4.2 的离散点方案，不要再在双指数参数上继续死磕。

5. 避坑指南：雷电波建模里五个高频踩坑点

5.1 指数系数忘加负号，波形一路暴涨

有一次在参数卡里填时间常数时把负号落下了，仿真结果完全发散，电流一路涨到几千千安还不收敛。

原因是双指数函数要求指数项为exp(-t/τ)，如果写成了exp(t/τ)，衰减项变成增长项，波形自然爆炸。

解决方法是统一按“先负号后数值”的习惯填写：在 CSMF 里写-TIME / TAU_TAIL，在参数卡里填负的时间常数或负的系数。PSCAD 版本不同，有的组件要求参数本身就是负值，填-1.2e-6而不是1.2e-6，动手前先看组件帮助里的示例。

5.2 把微秒当成秒填，波形快到看不见

有次同事说模型“没有输出”，调了半天发现波形其实存在，只是压缩在 0.1 微秒内结束了。

原因是把 1.2 直接当成了 1.2e-6 填成“1.2 秒”，实际含义差了百万倍；另一种情况是反过来，把秒当微秒，波形慢到半天不衰减。

解决方法是所有时间参数先在草稿纸上换算成科学计数法，填进 PSCAD 后再检查参数显示是否有科学计数格式。我自己的习惯是在模型旁边挂一张参数说明表，把目标波形、时间常数、系数三列写清楚，填参数时逐项对照。

5.3 “乘 2”参数：不同版本里的经验修正值

资料里提到 833333.33 在实际设置时要再乘 2 并取负值，这个描述初看很费解。

调一个老版本的 PSCAD 内置组合源时发现，它的内部公式把时间扫描步长做了二分，导致同一个波形在参数填写上要和理论值相差 2 倍来补偿，所以文档里才要求“833333.33 又乘了 2”。

但这并不是通用规则。新版组件或 CSMF 自定义模型里直接填1 / 1.2e-6就是对的，不需要乘 2。如果是在老版本里复现，建议先用一个最简单的测试网络跑一下，实测波形峰值和波前时间都对得上，再确认是否需要乘 2。

提示：遇到“乘 2”“加系数”这类经验修正值，先记录它适用的软件版本，换版本后重新验证，不要直接沿用。

5.4 仿真步长太大，波前变钝、峰值丢失

雷电波波前只有几微秒，如果仿真步长取 0.5 微秒，上升沿只有两三个采样点，波形看起来不是平滑曲线而是锯齿状，峰值也可能偏低。

原因是离散采样间隔太大，峰值点落在采样点之间的概率很高。

解决方法是把步长压到波前时间的 1/50 以下：1.2/50 建议 0.01 微秒到 0.02 微秒，8/20 可以放松到 0.05 微秒。仿真总时长不变，计算时间增加不多，但波形质量提升明显。

5.5 半峰值时间测错参考点

给定 1.2/50 波形，有人用游标量出来的半峰时间是 48 微秒，有人量出来是 51 微秒，排除模型差异后大概率是参考点取错。

原因是半峰值时间的定义是从波形起点到下降到 50% 峰值的时间，有人习惯从峰值开始算，这样量出来会多出波前那一段。

解决方法是测量前先把游标零点锁定在波形起始点，再量到下降沿 50% 处。PSCAD 的波形窗口支持设置两个游标并显示差值，核对一下差值是否等于起始点到半峰值点的时差，避免视觉误差。

6. 波形验证：用三个指标量化仿真结果与目标偏差

6.1 三个指标定义

雷电波建完不等于结束，必须用量化指标验证波形是否正确。我固定看三个值：峰值电流、视在波前时间、半峰值时间。峰值直接决定 SPD 选型和应力分析；视在波前时间对应电流上升率；半峰值时间反映波尾能量。两个时间指标在数值提取上要按“波形定义”来做，峰值容易找，波前时间按 10% 到 90% 上升段外推，半峰值时间从波形起点计。

6.2 快速验证脚本

从 PSCAD 里导出波形 CSV 后，我习惯直接跑下面这段脚本，自动算三个指标和误差。

import pandas as pd import numpy as np def verify_waveform(path, target_peak, target_front, target_tail): df = pd.read_csv(path, skiprows=2) t = df.iloc[:, 0].values i = df.iloc[:, 1].values i_peak = np.max(i) peak_idx = int(np.argmax(i)) i_10 = 0.1 * i_peak i_90 = 0.9 * i_peak idx_10 = int(np.argmin(np.abs(i[:peak_idx] - i_10))) idx_90 = int(np.argmin(np.abs(i[:peak_idx] - i_90))) t_front = (t[idx_90] - t[idx_10]) / 0.8 half = 0.5 * i_peak idx_half = int(np.argmin(np.abs(i[peak_idx:] - half))) + peak_idx t_tail = t[idx_half] print(f"峰值: {i_peak:.2f} kA, 误差 {100 * (i_peak - target_peak) / target_peak:+.2f}%") print(f"波前: {t_front * 1e6:.2f} us, 误差 {100 * (t_front - target_front) / target_front:+.2f}%") print(f"半峰: {t_tail * 1e6:.2f} us, 误差 {100 * (t_tail - target_tail) / target_tail:+.2f}%") verify_waveform("output_12_50.csv", 50, 1.2e-6, 50e-6)

脚本核心是verify_waveform函数：读 CSV，跳过 PSCAD 导出的前两行表头，取时间列和电流列，依次计算三个指标和误差。skiprows=2是 PSCAD 导出文件的通用格式，如果你用的版本表头行数不同需要调整。参数target_front和target_tail都按秒传，浮点计算保持量纲一致，输出时再转成微秒方便阅读。

6.3 校准习惯

从那以后，我每次在 PSCAD 里搭完雷电波模型，都会把波形导出来跑一遍这个验证脚本，三个指标全部落在容差内，才敢把波形往系统模型里丢。峰值对不上就先查幅值系数，半峰对不上就查 TAU_TAIL，波前对不上就查 TAU_RISE，按这个顺序排查基本没有例外。这套验证流程看起来多花两分钟，实际上省掉了后面整板仿真的所有返工，已经成了我的默认习惯。希望帮到你。

本文还有配套的精品资源，点击获取