No.19 基于电力系统仿真软件ETAP的潮流计算，短路计算，稳定计算仿真，包含以下仿真文件: 单机无穷大系统的潮流计算、短路计算、稳定计算仿真。 多机IEEE39节点系统的潮流计算、短路计算、稳定计算仿真。 短路计算包含不同的短路类型。 稳定性分析，包含临界切除时间分析，提高暂态稳定性的措施验证

搞电力系统仿真的人都知道，ETAP这玩意儿就像电工的瑞士军刀。今天咱们直接上硬菜，拿单机无穷大系统和IEEE39多机系统当案例，手把手盘一盘潮流、短路、稳定这三个核心模块怎么玩。

先说单机无穷大系统的潮流计算。别被名字唬住，其实就是发电机接个无穷大母线。在ETAP里新建项目时，发电机参数设置有个坑要注意：

Generator.Voltage = 22.0;  // kV
Generator.Xd = 0.3;       // 同步电抗别手滑输成0.03

这里X_d的单位是标幺值，新手容易搞错量纲。跑完潮流记得看母线电压是否在1.0±0.05pu范围内，超出就得调变压器分接头或者加无功补偿。

短路计算最刺激的是故障类型切换。三相短路虽然凶残但计算简单，ETAP的脚本控制可以这么写：

Fault.Type = LLG;  // 两相接地故障
Fault.Location = Bus23;
TimeDelay = 0.1;   // 100ms后动作

接地电阻参数千万别漏，实际系统中这个值直接影响短路电流幅值。对比不同故障类型的计算结果会发现，两相短路电流可能比三相还大——这反常识的现象跟系统零序阻抗有关。

No.19 基于电力系统仿真软件ETAP的潮流计算，短路计算，稳定计算仿真，包含以下仿真文件: 单机无穷大系统的潮流计算、短路计算、稳定计算仿真。 多机IEEE39节点系统的潮流计算、短路计算、稳定计算仿真。 短路计算包含不同的短路类型。 稳定性分析，包含临界切除时间分析，提高暂态稳定性的措施验证

稳定计算才是重头戏。拿临界切除时间（CCT）来说，用ETAP的时域仿真功能，逐步增加故障持续时间：

for (t=0.1; t<1.0; t+=0.05) {
   ApplyFault(t);
   if (RotorAngle > 120deg) break;
}

当功角曲线发散时，前一个t值就是CCT。想提升稳定性？试试在发电机出口加串联电容补偿，动态调整参数：

SeriesCapacitor.Enable = True;
SeriesCapacitor.Xc = GetOptimalXc(RotorSpeed);

这代码里的GetOptimalXc需要自己写算法，实时计算最佳容抗值。实测发现补偿量超过临界值反而会引起次同步振荡，把握度很重要。

切换到IEEE39节点系统时，数据导入建议用Python预处理：

import pandapower as pp
net = pp.convert.from_etap("IEEE39.ewd")
pp.drop_slack_transformers(net)  # 清理冗余数据

多机系统震荡模式分析得看特征值，ETAP的脚本输出要配合MATLAB做模态分析：

[V,D] = eig(JacobianMatrix);
damp_ratio = -real(D)./abs(D);

重点跟踪阻尼比小于5%的振荡模式。某次仿真发现0.8Hz的低频振荡，最后定位到是区域间联络线传输功率超限。

玩仿真最忌只跑标准案例。尝试修改IEEE39的负荷增长曲线，观察电压稳定裕度的变化：

Load.GrowthRate = 2.5%;  // 年增长率
TimeHorizon = 10;         // 十年规划
RunContinuationPowerFlow();

当PV曲线出现拐点时，系统就到了崩溃临界点。这时候加装STATCOM比扩建线路更划算，不信你试试把补偿容量从50MVar调到80MVar，电压曲线立马支棱起来。

说到底，仿真软件只是工具，真正的功夫在参数设置背后的物理直觉。下次遇到计算结果诡异时，先别怀疑软件bug，检查接地方式设置或者变压器分接头位置，往往会有惊喜（或者惊吓）。