电磁暂态与机电暂态

很多读者第一次听到“电磁暂态仿真”时，会觉得这个词有点大。电磁、暂态、仿真，三个词叠在一起，好像马上就要进入一套很重的理论体系。

电磁暂态仿真常写作 EMT，全称是 Electromagnetic Transient simulation；机电暂态仿真的英文全称是 Electromechanical Transient simulation。

关于电磁暂态

笔者的经验定义

下面这句话不是某一本书或某个标准里的权威定义，而是笔者为了讨论方便，把课程资料、EMT 工具说明、行业资料和自己的学习经验合在一起形成的经验定义：

电磁暂态仿真，简称 EMT 仿真，英文全称是 Electromagnetic Transient simulation，是用时域瞬时值方程，逐时间步计算电力系统在开关、故障、储能交换、线路传播、非线性元件和快速控制作用下的电压电流波形。

这个定义里有几个关键词。

• 时域：直接沿时间轴一步一步计算，而不是只在某个频率上求稳态。
• 瞬时值：保留 $u_a(t)$、$u_b(t)$、$u_c(t)$、$i_a(t)$、$i_b(t)$、$i_c(t)$ 这样的波形。
• 电路方程：每一步都要求解由元件、节点、支路和控制器共同形成的方程。
• 暂态事件：开关、故障、投切、保护动作、控制限幅、非线性饱和等都会触发暂态。

这里说的“方程”不是单指某一种形式。对电阻、电感、电容、控制器这类集总模型来说，底层常常是常微分方程 ODE。最常见的例子就是电感和电容：

整理成状态方程，就是：

但完整电路不只有元件自己的微分关系，还要加上节点电压法、支路约束、控制器限幅等代数关系；装配到完整网络以后，更接近一组微分-代数方程 DAE。对线路、电缆这类分布参数模型来说，原始物理模型还可能来自偏微分方程 PDE 或传输线方程，但在 EMT 程序里通常会被离散成特征线模型、历史源或等效网络。最后落到每一个时间步，求解器面对的往往是一组代数方程，例如后面文章反复出现的节点方程 $Gv_k=b_k$。

这并不是说普通 EMT 算例里经常直接求 PDE。对入门级 RLC、节点电压法、Dommel 伴随模型来说，ODE 和 DAE 才是主线。PDE 更常出现在它的“前一层物理来源”里：比如架空线、电缆、母线和 GIS 中的行波传播，雷电冲击和操作过电压中的波过程，高频暂态下变压器绕组、接地网、电缆屏蔽层的分布效应，以及少数需要把有限元电磁场模型同外部电路联立的设备级仿真。进入系统级 EMT 工具以后，这些模型往往已经被整理成可装配的端口模型，而不是让节点电压法在全网层面直接求 Maxwell 方程。

业内资料如何描述 EMT

不同资料给 EMT 下定义时，侧重点不完全一样。它们不一定互相矛盾，只是站在不同位置说同一件事。

MathWorks 对 EMT simulation 的公开解释是： 它用于建模和分析电力系统中快速、短时的 electromagnetic transients。这些暂态常由 IBR，也就是 inverter-based resources，逆变器型资源的开关动作或设备投入、退出引起，并可能造成过电压、过电流或系统稳定问题。¹

PSCAD/EMTDC 本身是 EMTDC 求解器和图形界面 PSCAD 的组合。 它的核心是对电力系统电磁暂态过程进行时域仿真。换句话说，它关心的不是只求一个基波相量结果，而是在给定网络、元件模型和控制模型后，沿时间轴一步步求出瞬时电压、电流和内部状态。

浙江大学徐政教授在相关课程中介绍了机电暂态和电磁暂态的区别。 机电暂态从大系统稳定性出发，研究发电机转子之间的相对摇摆、机电能量交换和同步稳定；电磁暂态则更关注局部、详细、短时间尺度的电磁过程，网络通常采用三相模型或全相模型，物理量为瞬时值而不是相量。这个角度特别适合初学者理解：EMT 并不只是“更精细”，而是研究对象、主导物理过程和变量选择都发生了变化。²

《电力系统电磁暂态仿真》第 1 章第 1-2 页讨论了电力系统暂态过程的时间范围。 书中把偏短时间尺度上、主要由电感磁场能量和电容电场能量相互作用主导的过程称为电磁暂态；把偏长时间尺度上、主要受旋转电机中机械能量和网络中电能量相互作用影响的过程称为机电暂态；中间还有一个受两类作用共同影响的暂态稳定区域。³

从系统分析角度看，电磁暂态问题通常可以理解为：求解一组满足基尔霍夫定律的一阶微分方程组，用来描述 RLC 电路或更一般电力网络在特定激励之后的动态行为。这里的“特定激励”可以是开关操作、故障、雷击、断路器重燃、绝缘配合、操作过电压、故障恢复或控制器动作。不同激励落在不同时间尺度上，所需的元件模型、仿真步长和等效方式也会不同。

两个小例子：EMT 和机电暂态到底差在哪

先看一个 EMT 味道很强的小例子。假设一个直流电源通过电阻 $R$ 和电感 $L$ 给支路供电，开关在 $t_s$ 时刻断开。为了让电感电流有去处，开关断开后电感通过一个吸收电阻 $R_d$ 续流。断开前，如果支路已经接近稳态，电感电流近似为：

断开后，电感中的磁场能量通过 $R_d$ 释放，满足：

于是：

吸收电阻上的电压为：

这个例子里，我们关心的是开关断开之后每一个时刻的电流衰减、电压尖峰和能量释放。也就是说，问题的输出不是“系统还稳不稳定”这种宏观判断，而是 $i_L(t)$、$u_d(t)$ 这样的瞬时波形。这就是典型的电磁暂态视角。

再看一个最经典的机电暂态例子：暂态功角稳定。假设一台同步发电机接在一个很大的电网后面，系统在 $t_f$ 时刻发生三相短路，保护在 $t_c$ 时刻切除故障线路。故障期间，发电机的电磁输出功率会明显下降，但汽轮机或水轮机给进来的机械功率 $P_m$ 来不及立刻改变，于是转子会被加速；故障切除后，电磁功率恢复，转子又开始减速和摆动。

这里所谓“简化”，不是说不解微分方程了，而是说不再详细追踪发电机定子绕组、电网三相电压电流和开关触头附近的电磁波形。我们把外部电网等效成一个无限大母线，把发电机内部电磁过程压缩成“电气输出功率和转子功角之间的关系”，于是问题的主要状态量变成转子功角 $\delta$ 和转速偏差。

简化到单机无穷大系统时，常见的摆动方程可以写成：

其中 $\delta$ 是发电机转子功角，$M$ 是惯性系数，$D$ 是阻尼系数，$P_m$ 是机械输入功率，$P_e(\delta)$ 是电气输出功率。常见简化下：

为了数值求解，通常会把这个二阶微分方程改写成两个一阶方程。令：

则：

如果采用最简单的显式欧拉法，时间步长为 $\Delta t$，那么第 $k$ 步可以写成：

故障发生、故障切除、线路退出运行，都会改变网络等效电抗，因而改变 $P_{\max}$。程序会按不同阶段使用不同的 $P_{\max}$ 继续推进 $\delta_k$ 和 $\omega_k$。于是问题变成：转子角 $\delta(t)$ 会不会越摆越大？发电机还能不能保持同步？这个问题关心的是几十毫秒到数秒内的机电能量交换，而不是断路器触头附近微秒级的电压波头。

还有一些问题会落在两类暂态的交界处，比如次同步谐振 SSR，英文是 subsynchronous resonance。它一边涉及发电机转子轴系的机械扭振，另一边又涉及串补线路、电网电气谐振和电磁能量交换。只看机电暂态，可能看不清电气网络在次同步频率附近的细节；只看局部 EMT，又可能需要额外保留轴系模态和机械阻尼。工程上常见的处理，是根据研究目标选择小信号分析、机电暂态、EMT，或者机电-电磁混合仿真。

所以，同样是“开关在 $t_s$ 时刻动作”，两类仿真问的问题不一样：

从特征时间尺度理解两类暂态

从数学上看，这两类问题都可以落到微分方程或微分-代数方程上。差别不在于“谁有微分方程”，而在于微分方程里保留了哪些状态量，以及这些状态量对应的时间常数有多大。

先看数学上的时间尺度

把一个非线性系统在某个工作点附近线性化，常常可以写成：

如果 $A$ 有一个特征值 $\lambda$，那么对应模态大致长成：

当 $\lambda=-\alpha+j\omega$ 时，$\alpha$ 决定衰减速度，$\omega$ 决定振荡频率。粗略地说：

这就是时间尺度的数学来源。系统里如果存在很大的 $\alpha$ 或很大的 $\omega$，就说明有很快的衰减或振荡。为了把这个快模态算出来，仿真步长 $\Delta t$ 必须比它的时间常数或周期小很多。

电磁暂态：电感、电容和线路带来的快模态

对电磁暂态来说，状态量通常来自电感、电容、线路和开关。一个最简单的 $RL$ 回路有：

它的时间常数是：

一个 $RC$ 回路有：

它的时间常数是：

如果是 $LC$ 振荡，理想情况下有：

这些时间常数可能很小。电力电子开关、杂散电感电容、线路行波、断路器重燃、雷电冲击和操作过电压，都会把问题推到微秒甚至更短的时间尺度。于是 EMT 的步长常常要取到微秒级，甚至更小；如果研究对象只是较低频的电磁振荡，也可能用几十微秒。

机电暂态：转子惯性带来的慢模态

机电暂态的主角则通常是转子惯性。把前面的摆动方程在工作点 $\delta_0$ 附近线性化。设：

由于：

在 $\delta_0$ 附近有：

其中：

$K_s$ 通常可以理解成工作点附近的功角同步刚度。它描述的是：当转子功角 $\delta$ 偏离原来的平衡点一点点时，电气输出功率 $P_e$ 会增加得有多快。如果 $K_s$ 很大，说明功角稍微变化，电气功率就会明显变化，电磁转矩对转子的“拉回作用”比较强；如果 $K_s$ 很小，说明功角变化以后，电气功率给出的恢复作用很弱，系统就更容易出现低频、大幅度摆动。这里的“刚度”不是机械弹簧的刚度，而是功角和电磁功率之间的局部斜率。

代回摆动方程，可以得到一个类似二阶振子的形式：

这个式子和机械振动里的“质量-阻尼-弹簧”模型很像。$M$ 对应惯性，$D$ 对应阻尼，$K_s$ 对应恢复刚度。区别只是这里振动的变量不是位移，而是功角偏差 $\Delta\delta$；恢复力也不是弹簧力，而是由电磁输出功率随功角变化产生的同步作用。

如果暂时忽略阻尼，它的自然振荡角频率近似为：

这里的 $M$ 来自转子惯性。更直观地说，$M$ 越大，同样的不平衡功率 $P_m-P_e$ 造成的角加速度就越小，转子越“不愿意被加速或减速”。所以 $\omega_n$ 的表达式里，$K_s$ 在分子，$M$ 在分母：同步刚度越大，摆动越快；惯性越大，摆动越慢。

转子很重，惯性很大，所以 $\omega_n$ 通常不会特别高。很多功角摆动、低频振荡和区域间振荡落在零点几赫兹到几赫兹的量级，对应周期是几百毫秒到几秒。因此机电暂态仿真可以使用毫秒级、甚至更大的步长，而不必像 EMT 那样追踪微秒级电压波头。

用一组数量级把差别拉开

为了让这个数量级更具体，取几组并不极端的参数做个估算。

先看一个电磁暂态里的 $LC$ 振荡。若：

则：

对应频率为：

周期大约是：

如果希望一个周期内至少有几十个采样点，步长自然会落到几微秒到十几微秒的量级。若再考虑更小的杂散电感和杂散电容，例如：

则：

对应周期只有约 $6.3\ \mu s$。这类模态如果对工程结论重要，仿真步长就必须进一步缩小。

再看一个机电暂态里的功角摆动。采用常见的标幺化摆动方程写法，可以把惯性项近似写成：

这里的 $H$ 是同步机惯性常数，单位通常写成秒。它表示额定转速下转子储存的动能和额定容量的比值。$H$ 越大，说明同样容量的机组转子储能越多，受到功率不平衡时速度变化越慢。除以同步电角速度 $\omega_s$ 后，就得到前面摆动方程里使用的惯性系数 $M$。

若系统频率为 $50\ \mathrm{Hz}$，则：

取同步机惯性常数：

则：

如果在某个运行点附近，功角同步刚度取一个量级值：

则机电模态的自然角频率约为：

对应频率为：

周期大约是：

这个数量级和前面的 $LC$ 振荡相比，已经相差约 $5.6\times 10^3$ 倍；如果和高频杂散参数模态相比，则可以相差 $10^5$ 倍以上。

模型取舍要跟着时间尺度走

所以，“快”和“慢”不是凭感觉分出来的，而是由模型保留的状态变量决定的。EMT 保留电磁储能和开关波形细节，特征时间常数小；机电暂态保留转子机械运动和功角摆动，惯性大、特征频率低。真正复杂的工程问题，比如次同步谐振、HVDC 控制和新能源并网振荡，往往会把这两类时间尺度耦合在一起，这也是为什么有时需要机电-电磁混合仿真。

小结

电磁暂态仿真不是一个神秘软件类别，而是一种看待电力系统动态的方式。

当问题可以被工频相量充分描述时，没有必要强行使用 EMT。相量法更轻、更快，也更适合长时间大系统分析。

但当波形细节本身会改变结论时，例如开关点波、暂态过电压、谐振、故障电流、保护测量、功率电子快控制和非线性饱和，电磁暂态仿真就变得必要。

从算法角度看，EMT 的核心工作可以概括成一句话：

把电路中的连续时间微分关系，按固定步长翻译成每一拍都能求解的网络代数方程，并输出电压电流的瞬时波形。

理解了这一点，再往后看 Dommel 形式、节点电压法、$G$ 矩阵、历史项和实时仿真，就不再是零散概念，而是一条连续的工程路线。

参考资料与术语来源

• resource_forAI/deep-research-report.md：电磁暂态实时仿真的定义、范围、数学基础与术语边界。
• H. W. Dommel, “Digital Computer Solution of Electromagnetic Transients in Single- and Multiphase Networks,” IEEE Transactions on Power Apparatus and Systems, 1969.
• IEEE/CIGRE Joint Task Force, “Definition and Classification of Power System Stability,” IEEE Transactions on Power Systems, 2004.
• N. Hatziargyriou et al., “Definition and Classification of Power System Stability Revisited & Extended,” IEEE Transactions on Power Systems, 2020.
• CIGRE WG C4.56, Technical Brochure 881, Electromagnetic transient simulation models for large-scale system impact studies in power systems having a high penetration of inverter-connected generation, 2022.
• NERC, Electromagnetic Transient Modeling for BPS-Connected Inverter-Based Resources, 2023.
• NERC, EMT Analysis in Operations Planning, 2025.
• PSCAD/EMTDC、EMTP、RTDS、OPAL-RT HYPERSIM/eHS 等厂商公开技术资料。