简介面向电力系统电磁暂态仿真初学者与防雷设计人员这份资源聚焦PSCAD平台下雷电波的建模与验证帮助读者理解冲击电流波形的数学表达与参数设置逻辑。包内为1个doc文档约120KB以图文结合方式呈现建模过程与仿真结果分析。内容围绕GB18802标准中浪涌保护器I级与II级分类试验展开分别给出50kA、1.2/50波形与20kA、8/20波形的双指数函数表达式并说明PSCAD中波前时间与半峰值时间需取负值输入的原因及具体数值换算同时附有仿真图与曲线拟合难点的讨论。已有2006人学习适合需要对照标准搭建雷电波模型、排查参数设置问题或补充防雷仿真思路的读者参考。1. 雷电波在 PSCAD 里到底怎么建从 GB18802 的 1.2/50 与 8/20 说起做低压配电防雷仿真的人绕不开一个具体动作在 PSCAD/EMTDC 里把一条符合标准要求的冲击电流波形搭出来然后拿它去考核浪涌保护器SPD的残压和通流。这份资源就是围绕这个动作展开的——一份 PSCAD 建模文档主题是雷电波模拟附带仿真图和分析。它解决的不是“防雷理论是什么”而是“标准里写的 50kA、1.2/50 和 20kA、8/20在 PSCAD 里到底填哪几个数、为什么填负数、波形对不上该改谁”。适合正在做 SPD 选型校核、低压系统过电压仿真、或者课程设计里要交一张像样冲击波形图的从业者和学生。下面按我拆这份文档的顺序讲先立住波形数学再落到 PSCAD 元件参数最后说拟合和踩坑。2. 双指数波形先立住1.2/50 与 8/20 的数学表达2.1 双指数公式的物理含义与两个时间常数标准里描述的冲击电流波形工程上普遍用双指数函数来近似形式是 I(t) A·[exp(-αt) − exp(-βt)]。前一项管衰减后一项管上升两个指数相减就得到一条先快速爬升、到峰值后缓慢拖尾的曲线。资源里给出的两条式子很直接I 级分类试验 50kA、1.2/50I 50 × [exp(-50e-6·t) − exp(-1.2e-6·t)]II 级分类试验 20kA、8/20I 20 × [exp(-20e-6·t) − exp(-8e-6·t)]这里有个容易看岔的地方1.2/50 里的 1.2 是视在波前时间微秒50 是半峰值时间微秒但公式里指数上的系数并不是 1/1.2 和 1/50 这么直白。资源正文里 50kA 那条式子指数系数写的是 50e-6 和 1.2e-6这其实是把“时间常数”和“波形标称时间”混在一处表述了。真正决定波形形状的是两个指数的时间常数标称的 1.2/50 是拟合后从曲线上量出来的特征量不是直接等于时间常数。这一点如果不先想清楚后面在 PSCAD 里填参数一定翻车。我一般会先把公式在纸上或脚本里画一遍确认峰值出现在哪个时刻、峰值大概多少再进 PSCAD。因为双指数波形对参数极其敏感指数系数差一个量级波形就从“冲击”变成“阶跃”或者“慢爬坡”。2.2 用 Python 先把两条曲线画出来核对进 PSCAD 之前先用几行脚本把波形跑出来是省时间的做法。下面这段只做一件事按资源给的式子生成两条曲线打印峰值和峰值时刻方便和标准波形对照。import numpy as np def double_exp(A, alpha, beta, t): # A 为幅值系数alpha 控制衰减beta 控制上升 return A * (np.exp(-alpha * t) - np.exp(-beta * t)) t np.linspace(0, 100e-6, 20000) # 100 微秒窗口足够看 1.2/50 和 8/20 # I 级50kA按资源式子 alpha50e-6, beta1.2e-6 i1 double_exp(50.0, 50e-6, 1.2e-6, t) # II 级20kA按资源式子 alpha20e-6, beta8e-6 i2 double_exp(20.0, 20e-6, 8e-6, t) for name, w in ((1.2/50, i1), (8/20, i2)): idx np.argmax(w) print(name, 峰值, round(w[idx], 3), kA, 峰值时刻, round(t[idx]*1e6, 3), us)跑完你会发现一个关键现象按资源正文这组系数直接代进去峰值时刻和标称的波前时间对不上峰值也不一定正好落在 50kA 或 20kA 上。这不是脚本写错而是双指数公式里“幅值系数 A”和“标称峰值”本来就不是一回事指数系数和标称时间也不是一一对应。资源正文里那句“拟和这个曲线有点难度”说的就是这件事。所以正确的心态是公式给的是形状骨架A 和两个指数系数都要按目标波形反调不能照抄。提示先用脚本把目标波形的峰值时刻、半峰值时刻量出来再进 PSCAD 调参比在仿真界面里盲试快得多。3. 落到 PSCAD 元件冲击电流源参数怎么填3.1 用哪个元件、时间常数为什么是负值PSCAD 里搭冲击电流源常见做法是用受控电流源配合自定义函数或者直接用库里带指数形式的源元件。资源正文提到的设置方式很有代表性波前时间相关参数填 833333.33 再乘 2、且取负值半峰值时间相关参数填 1/20000、也取负值。这两个数看着玄学其实来源是“把时间常数换算成 PSCAD 内部按秒计、且以负指数形式表达”的结果。PSCAD 的指数源一般写成 exp(-t/τ) 这种形式τ 是时间常数。要让波形按 exp(-αt) 衰减就得让 1/τ α也就是 τ 1/α。当 α 取 50e-6 时τ 1/(50e-6) 20000所以半峰值那一路出现 1/20000 这个数当 β 取 1.2e-6 时τ 1/(1.2e-6) ≈ 833333.33所以波前那一路出现 833333.33。取负值是因为元件内部按 exp(-t/τ) 计算而双指数里上升项是减号符号不处理就会得到一条方向相反的曲线。至于“又乘了 2”是部分源元件对上升沿时间常数的定义和标准时间常数差一个因子属于具体元件口径问题换元件就要重新核对不能当通用常数背。3.2 参数填写对照表与操作步骤把上面这层关系整理成一张表填参数时对着看比记数字可靠目标波形衰减项系数 α对应时间常数 1/α上升项系数 β对应时间常数 1/βPSCAD 填写负值50kA, 1.2/5050e-6200001.2e-6833333.33半峰值 -1/20000波前 -833333.33×220kA, 8/2020e-6500008e-6125000半峰值 -1/50000波前 -125000×2操作步骤按这个顺序走新建工程放置冲击电流源元件确认它的指数表达式是 exp(-t/τ) 形式不是 exp(-τt) 形式。这一步决定后面填倒数还是填系数。按上表把两个时间常数填进去注意负号。波前那一路如果元件文档要求乘 2就乘 2并在工程里留注释说明依据。幅值系数先按标称峰值填 50 或 20跑一次短时仿真看峰值实际落在多少。用测量元件抓峰值和峰值时刻和 1.2/50、8/20 的标称值比对偏差大就回到第 2 步微调。仿真设置建议PSCAD 里对应字段 - 仿真步长1e-7 s 量级冲击波形前沿只有微秒级步长太粗会削峰 - 仿真总时长200e-6 s 起步够看完整拖尾 - 输出通道电流瞬时值 峰值保持方便量半峰值时刻步长这条是血泪经验。1.2 微秒的波前如果步长给到 1e-5 秒一个上升沿只采到一两个点画出来的波形又矮又胖你会误以为是参数填错其实是采样不够。我一般把步长压到 1e-7 秒总时长 200 微秒单次仿真几秒钟出结果调参节奏才起得来。3.3 幅值系数与标称峰值的换算很多人卡在“我填了 50为什么峰值只有三十几”。原因是双指数在峰值时刻 t_p 处两个指数项都没到 1括号里的值小于 1所以 A50 得到的峰值必然小于 50。正确做法是先算出括号项在峰值时刻的值再用标称峰值除以它得到应该填的 A。峰值时刻可以对 I(t) 求导令其为零解出来也可以直接用脚本扫点找最大值。资源正文里没有展开这一步但实际建模绕不过去属于“合格从业者最可能用的补全”。import numpy as np from scipy.optimize import minimize_scalar def shape(t, alpha, beta): return np.exp(-alpha * t) - np.exp(-beta * t) alpha, beta 50e-6, 1.2e-6 # 找括号项最大值对应的时刻 res minimize_scalar(lambda x: -shape(x, alpha, beta), bounds(0, 100e-6), methodbounded) tp res.x peak_shape shape(tp, alpha, beta) A_needed 50.0 / peak_shape # 目标峰值 50kA 对应的幅值系数 print(峰值时刻 us:, round(tp*1e6, 3), 括号项峰值:, round(peak_shape, 4), 应填 A:, round(A_needed, 3))这段的逻辑是波形形状由 α、β 定死峰值时刻和括号项峰值也就定了剩下的只是用一个比例系数把高度拉到目标值。参数说明alpha、beta 换成 20e-6、8e-6 就得到 8/20 那组的 A。跑出来的 A 通常明显大于标称峰值这是正常的不是填错。4. 波形拟合与验证把“像不像”变成可量化的比对4.1 用特征量比对代替肉眼判断资源正文说“拟和这个曲线有点难度”难点不在公式而在“像不像”没有量化标准。我的做法是固定四个特征量来比对峰值、峰值时刻近似视在波前时间、半峰值时刻、以及 t0 附近的上升斜率。前三个直接量第四个看波形起步是否过冲。把这四个量和标准标称值列在一起偏差超过 10% 就回去调 α、β。特征量1.2/50 目标8/20 目标量取方法峰值50 kA20 kA电流通道最大值视在波前时间1.2 us8 us10%到90%上升时间×1.67 近似半峰值时间50 us20 us峰值后降到一半的时刻起步过冲无无看 t 很小时是否反向视在波前时间严格定义是按 10% 和 90% 两点连线的斜率外推到零点和峰值点的时间工程上常用 1.67 倍上升时间近似够用。半峰值时间就是峰值之后电流降到峰值一半的那个时刻直接从波形上读。4.2 拟合不上的三种典型原因第一种元件表达式口径没核对。有的源元件内部是 exp(-τt)你按 exp(-t/τ) 填倒数波形会完全走样。第二种幅值系数没换算峰值永远差一截越调越乱。第三种步长太粗波形被采样削平特征量全偏。这三种在资源正文里没有逐条列但都是实际建模时最先撞上的。排查顺序建议先看元件表达式再算幅值系数最后压步长别一上来就怀疑公式。注意不要用“波形看起来差不多”作为验收标准。把四个特征量写进工程注释下次改参数时一眼能看出动了哪个量。5. 避坑与排查五个真实会翻车的点5.1 现象波形方向反了电流先负后正原因双指数里上升项是减号PSCAD 源元件按 exp(-t/τ) 计算符号没处理减号变成了加号方向。解决把上升项对应的时间常数取负值或在自定义函数里显式写成相减形式跑一次短时仿真确认 t 很小时电流从零往正方向走。5.2 现象峰值怎么调都到不了 50kA原因把标称峰值直接当幅值系数填了忽略了括号项峰值小于 1。解决按 3.3 的脚本算出应填的 A通常比 50 大一截填进去再微调。5.3 现象波前时间明显偏大波形又矮又胖原因仿真步长太粗微秒级前沿被采样削平。解决步长压到 1e-7 秒量级总时长 200 微秒重新跑特征量会立刻回到合理区间。5.4 现象换了源元件后原参数全失效原因不同元件的指数表达式口径不同有的按时间常数有的按系数有的自带 2 倍因子。解决换元件先读它的表达式定义把 1/α、1/β 重新换算一遍别把上一套数字直接搬过去。5.5 现象半峰值时刻量出来和 50 微秒差很多原因把“半峰值时间”和“衰减时间常数”混为一谈直接拿 1/α 当半峰值时间。解决半峰值时间是从波形上量出来的特征量和时间常数不是一回事以波形实测为准反调 α 直到量出来的值对上。6. 进阶技巧把调参过程固化成可复用的参数扫描调双指数波形最耗时的不是建一次模而是反复试参数。我的习惯是把 α、β、A 做成三个可扫的变量用脚本批量跑特征量把满足峰值和半峰值时间容差的那组参数挑出来再回填 PSCAD。下面这段是思路骨架跑的是纯数学层不依赖 PSCAD 接口import numpy as np def features(A, alpha, beta): t np.linspace(0, 200e-6, 200000) w A * (np.exp(-alpha*t) - np.exp(-beta*t)) idx np.argmax(w) peak, tp w[idx], t[idx] # 半峰值时刻峰值后首次降到一半 after np.where(w[idx:] peak/2)[0] th t[idx after[0]] if len(after) else np.nan return peak, tp, th best None for alpha in np.linspace(40e-6, 60e-6, 21): for beta in np.linspace(0.8e-6, 1.6e-6, 21): A 50.0 / (np.exp(-alpha*1.2e-6) - np.exp(-beta*1.2e-6)) # 粗估 peak, tp, th features(A, alpha, beta) err abs(tp*1e6 - 1.2) abs(th*1e6 - 50)/50 if best is None or err best[0]: best (err, alpha, beta, A, peak, tp, th) print(best)逻辑说明外层扫 α、β内层用目标峰值反算 A再用特征量误差挑最优。参数说明容差权重可以按项目要求调波前时间卡得紧就把第一项权重加大。跑出来的最优组不一定物理上最“标准”但一定最贴近你设定的验收口径这就是把玄学调参变成可复现流程的价值。从那以后我每次建冲击波形都强制先跑一遍特征量脚本再进 PSCAD不再靠肉眼在图上比划。这份资源的仿真图和分析可以当起点但真正让模型站得住的是你自己那套量化比对习惯。希望帮到你。本文还有配套的精品资源点击获取