2. 电磁暂态仿真基础

电磁暂态仿真的核心任务，是在微秒级时间尺度上求解电网中由操作、故障或雷击引起的宽频暂态过程。与机电暂态关注发电机功角、转速等秒级慢变量不同，电磁暂态必须保留线路分布参数、开关电弧、变压器饱和等快速动态。EMTP（Electromagnetic Transients Program）类程序正是为此设计的求解框架，其算法基础、积分方法、元件模型和步长选择共同决定了仿真的精度与稳定性。

一、EMTP类算法的基本思想

EMTP类算法采用节点法建立网络方程，将电感、电容等储能元件通过数值积分转化为等效电阻与历史电流源并联的诺顿等效支路，再对每个时间步求解节点导纳方程 Y·v = i。由于导纳矩阵 Y 在步长不变时保持恒定，只需在仿真开始时进行一次三角分解，后续每步仅做前代回代，计算效率高。这一“固定导纳矩阵+历史项更新”的结构，是EMTP类程序能够处理数千节点系统的关键。

与状态变量法相比，节点法不需要显式列出状态方程，元件建模更直观；代价是当网络含理想电压源回路或纯电容割集时，导纳矩阵可能奇异，需要引入支路或小电阻处理。实际使用中，开关用可变电阻或二值电阻表示，变压器用耦合电感或理想变压器加漏抗表示，线路用分布参数模型或集中参数π型等值。

二、数值积分方法：梯形法与后向欧拉法

EMTP最常用的积分方法是梯形法（Trapezoidal Rule）。对电感 L，其伏安关系 v = L·di/dt 经梯形积分离散后得到：

i(t) = (Δt / 2L)·[v(t) + v(t-Δt)] + i(t-Δt)
     = v(t)/R_eq + I_hist
其中 R_eq = 2L/Δt，I_hist = i(t-Δt) + v(t-Δt)/R_eq

对电容 C 同理可得 R_eq = Δt/(2C)，历史电流源由上一时刻电压和电流决定。梯形法具有二阶精度，且对纯电感、纯电容回路不产生数值阻尼，适合分析振荡型暂态。但它对刚性系统或含突变开关的网络可能产生数值振荡，此时可改用后向欧拉法（Backward Euler）：

i(t) = v(t)·Δt/L + i(t-Δt)   （电感）
R_eq = L/Δt，I_hist = i(t-Δt)

后向欧拉法是一阶精度，但具有强阻尼特性，能有效抑制开关操作后的非物理振荡。实际EMTP程序常采用临界阻尼调整（CDA）或梯形法与后向欧拉法组合：在开关动作后的几个步长内切换到后向欧拉，随后恢复梯形法，兼顾精度与稳定。

三、元件建模要点

RLC元件：电阻直接作为导纳 1/R；电感、电容按上述诺顿等效处理。串联RLC支路可先合并为等效阻抗再离散，减少节点数。

开关：理想开关用二值电阻模拟，闭合时取 R_on（如1 mΩ），断开时取 R_off（如1 MΩ）。开关状态变化时导纳矩阵改变，需重新分解或采用补偿法。补偿法将开关作为诺顿电流源补偿，避免每步重新分解矩阵，是EMTP中处理开关的常用技术。

变压器：需考虑漏抗、励磁支路和饱和特性。饱和用分段线性或非线性电感表示，通过磁链-电流曲线查表。三相变压器还需处理接线组别和零序通路。

线路：短线路可用π型集中参数；长线路必须用分布参数模型（Bergeron模型），考虑波过程、波阻抗和传播延时。Bergeron模型将线路视为波阻抗 Z_c 与历史电流源并联，两端通过延时 τ = l/v 耦合，能准确模拟行波传播。

四、步长选择与数值稳定性

步长 Δt 是电磁暂态仿真最关键的参数之一。它受三方面约束：

• 行波传播精度：分布参数线路的传播延时为 τ = l/v，通常要求 Δt ≤ τ/10，否则行波波头会被严重平滑。例如300 km架空线，波速约 3×10^5 km/s，τ ≈ 1 ms，则 Δt 宜取 50~100 μs；若线路仅10 km，τ ≈ 33 μs，Δt 需取 3 μs 左右。
• 最高关注频率：若需观察2 kHz的暂态分量，采样率至少10 kHz，即 Δt ≤ 100 μs；若关注雷电冲击（MHz级），则需 Δt ≤ 0.1 μs。
• 数值稳定性：梯形法对无阻尼LC回路临界稳定，可能产生持续振荡；后向欧拉法无条件稳定但精度低。当 Δt 过大时，离散化引入的虚拟阻尼或负阻尼会导致结果失真。实际中可通过比较不同步长结果是否收敛来判断步长是否合适。

五、一个简单示例：合闸空载线路

设一条300 km、波阻抗 Z_c = 300 Ω、传播时间 τ = 1 ms 的架空线，末端开路，电源经开关在 t=0 时合闸，电源电压为 e(t)=E·sin(ωt)，E=500 kV，f=50 Hz。采用Bergeron模型，步长取 Δt = 50 μs。合闸后，电压波以波速向末端传播，1 ms 后到达末端并反射，末端电压逐渐升高，理论上稳态末端电压可达电源电压的2倍左右（空载线路电容效应）。仿真中可观察线路首端电流和末端电压的波形：首端电流在合闸瞬间出现高频振荡，随后趋于稳态；末端电压在 2τ 后达到第一个峰值。若步长取 200 μs，波头会被明显拉平，峰值偏低，说明步长过大。

该示例中，开关用 R_on=0.01 Ω、R_off=1 MΩ 模拟，合闸时刻设为 t=0.001 s 以避开零时刻的数值启动问题。仿真总时长取 0.02 s，足以观察多次反射过程。通过对比不同步长下的末端电压峰值，可以直观理解步长对行波仿真精度的影响。

六、小结与注意事项

EMTP类算法以节点法为骨架，以梯形法/后向欧拉法为积分核心，通过诺顿等效将储能元件转化为电阻与历史电流源，实现高效求解。元件建模需根据关注频段选择模型复杂度：短线路可用π型，长线路必须用分布参数；开关需处理状态突变；变压器需考虑饱和与接线。步长选择需兼顾行波传播、最高频率和数值稳定性，建议通过步长减半对比来验证结果收敛性。掌握这些基础，才能为后续电磁-机电暂态协同仿真中的接口设计打下可靠基础。