3. 机电暂态仿真基础
机电暂态仿真关注的是电力系统在扰动后数秒到数十秒时间尺度上的动态行为，其核心特征是：发电机转子运动方程与网络代数方程联立求解，网络以工频相量形式描述，电磁暂态过程被平均化处理。理解机电暂态仿真，需要从三个层面入手：模型简化假设、同步发电机与励磁系统的建模方式、以及稳定分析的组织流程。
3.1 机电暂态模型的基本假设与相量域特点
机电暂态仿真建立在准稳态假设之上。其核心简化包括：网络元件用工频下的相量阻抗或导纳表示，忽略线路分布参数引起的波过程；发电机定子绕组中的直流分量和负序分量被忽略或近似处理，只保留与转子磁链相关的缓慢变化量；网络方程用节点导纳矩阵描述，形如 I = YV，其中 I、V 为节点注入电流和节点电压相量。
相量域仿真的关键特点是：电压、电流均以幅值和相角表示，时间步长通常取 1～10 ms，远大于电磁暂态仿真中微秒级的步长。这使得机电暂态仿真可以覆盖更长的物理时间，但代价是无法准确反映行波、开关暂态和高频振荡。因此，机电暂态仿真适用于功角稳定、电压稳定、频率稳定和低频振荡分析，而不适用于雷电冲击、操作过电压等场景。
在数值求解上，机电暂态仿真通常采用交替求解法：每个时步先由网络方程求出各发电机端电压和电流，再代入发电机微分方程积分求转子角和转速，然后更新注入电流，进入下一时步。这种交替求解结构清晰，但需注意步长与迭代收敛性的配合。
3.2 同步发电机与励磁系统建模
同步发电机是机电暂态仿真的核心动态元件。在机电暂态时间尺度下，定子暂态被忽略，发电机通常采用三阶或五阶模型。三阶模型包含转子运动方程和转子磁链方程，状态量为转子角 δ、转速偏差 Δω 和暂态电势 Eq′。其方程可写为：
dδ/dt = ω − ω₀
M dω/dt = Pm − Pe − D·Δω
Td0′ dEq′/dt = Efd − Eq′ − (Xd − Xd′)·Id其中 M 为惯性时间常数，D 为阻尼系数，Td0′ 为 d 轴开路暂态时间常数，Xd 和 Xd′ 分别为 d 轴同步电抗和暂态电抗，Id 为 d 轴电流分量。电磁功率 Pe 可近似表示为 Eq′ 与端电压、功角相关的函数。若需更精确地反映 q 轴动态和阻尼绕组效应，可扩展为五阶模型，增加 Ed′、Eq″、Ed″ 等状态量。
励磁系统建模通常采用 IEEE 标准模型，如 IEEE AC1A、DC1A、ST1A 等。以常见的 AC1A 为例，其结构包括：电压测量环节（一阶惯性，时间常数 Tr）、电压调节器（超前滞后或 PID 结构）、励磁机（一阶惯性，时间常数 Te）、以及励磁限幅环节。在相量域仿真中，励磁系统以微分方程形式与发电机方程联立，输出励磁电压 Efd 作用于发电机转子磁链方程。
实际建模时，需要根据设备铭牌和试验数据填写参数。例如某 600 MW 汽轮发电机，惯性时间常数 M 取 10 s，d 轴暂态电抗 Xd′ 取 0.28 p.u.，Td0′ 取 6.5 s；励磁系统 AC1A 的 Tr 取 0.02 s，Te 取 0.5 s，放大倍数 Ka 取 200。这些参数需与现场试验或厂家数据核对，不可随意假设。
3.3 电力系统稳定分析流程
机电暂态稳定分析通常按以下流程组织：
- 确定分析目标与故障场景：明确是功角稳定、电压稳定还是频率稳定，选择典型故障（三相短路、单相重合闸、切机、切负荷等），并确定故障位置、持续时间和保护动作时间。
- 建立系统数据模型：包括潮流数据（节点电压、注入功率）、网络参数（线路阻抗、变压器变比）、发电机及其控制系统参数、负荷模型（恒功率、恒阻抗或电动机模型）。
- 潮流计算获取初始状态：通过潮流计算得到各节点电压相量和发电机初始功角、电磁功率，作为动态仿真的初值。
- 设置仿真参数：选择积分方法（如隐式梯形法）、步长（通常 1～10 ms）、仿真时长（通常 10～30 s），并设定故障时序。
- 执行时域仿真：按交替求解结构逐步积分，记录发电机转子角、转速、母线电压、线路功率等关键变量。
- 结果判读与稳定判断：观察转子角差是否随时间发散，若最大功角差超过 180° 或持续增大，则判断为失稳；若振荡衰减并趋于稳定值，则判断为稳定。同时检查电压是否恢复至 0.9 p.u. 以上、频率是否在允许范围内。
3.4 相量域仿真的适用边界与注意事项
相量域仿真以工频相量为基本变量，其优势在于计算效率高、可覆盖长过程，但存在明确边界。当系统中含有大量电力电子设备（如 HVDC、FACTS、新能源变流器）时，相量域模型可能无法准确反映其快速控制动态，此时需与电磁暂态仿真协同。此外，相量域仿真对不对称故障的处理能力有限，通常需采用序分量或三相相量扩展模型。
在实际工程中，机电暂态仿真常用于：大扰动功角稳定校验、低频振荡分析、AGC 与调频策略验证、以及新能源并网后的系统稳定性评估。使用时应注意：模型参数需与现场一致；故障时序需与保护定值配合；仿真步长需与系统最小时间常数匹配，避免数值失稳。