企业官网建设流程全解析

前阵子有个电气专业的学生找我，说课程设计卡在“单机无穷大系统暂态稳定性仿真”上，老师给的题目偏偏不是最常规的三相短路，而是要分析两相接地性短路故障时发电机的转速变化。他试了两天，要么模型初始化报错，要么仿真出来转速飙升停不下来，连最基本的摇摆曲线都看不明白。我帮他从头把建模、故障设置、结果判读捋了一遍，顺便把这些年做暂态稳定仿真积累的细节也整理了一遍。这篇文章就把这套完整做法写出来，供正在做单机无穷大系统暂态稳定性仿真、需要观察两相接地故障下发电机转速摇摆曲线的读者参考。

先说结论：Simulink做这类题目，真正难的不是拖模型，而是你心里有没有一套“先物理、再模型、后参数”的思路。转速变化的背后是转子运动方程、电磁功率与机械功率的平衡关系；两相接地故障又比三相短路多了一层不对称的复杂性。把这些搞透了，建模和调参都只是顺势而为。

1. 单机无穷大系统为什么是暂态稳定研究的第一课——物理模型与等面积定则

1.1 单机无穷大模型里隐藏了什么假设

单机无穷大系统，英文常写成 SMIB（Single-Machine Infinite-Bus），结构上用一句话说就是：一台同步发电机，经升压变压器和双回输电线路，接到一个容量无限大、电压和频率恒定的理想母线上。

这个“无穷大母线”是个非常强的假设，意思是无论你这边怎么扰动，对侧的电压幅值、频率、相角都纹丝不动。实际电网当然不存在这种东西，但正因为这个假设，系统的动态特性完全由“一台机 + 它的外接电抗”决定，方便把暂态稳定的本质讲清楚。做研究和教学，第一步几乎都是它。

模型里另一个关键假设是机械功率输入恒定。原动机（汽轮机、水轮机）的调速器响应是秒级甚至更慢的，而暂态稳定的关心窗口是故障后第一秒到几秒的机电暂态过程，所以在这段时间里把机械功率当成常数，是工程上公认的简化。

发电机本身也做了简化：定子、转子绕组的电磁暂态可以忽略，只保留转子运动方程和与转子运动强相关的电磁功率表达式。这就是教科书里最经典的“实用暂态稳定模型”，核心状态量只有转子的相对角度和转速偏差。

1.2 两相接地短路为什么值得单独立项

电力系统中最常见的短路类型排序，按发生概率来算是：单相接地最多，两相短路、两相接地居中，三相短路最少。但按严重程度排则是反过来，三相短路最严重，两相接地排第二，然后是两相短路，单相接地最轻。

很多课程作业为了省事直接做三相短路，因为三相短路时电磁功率可以近似看成零，转子在故障期间纯粹加速，物理过程最简单。但三相短路在真实电网中占比很小，研究价值反而不如不对称故障。

两相接地短路（比如 A 相、B 相同时接地，记作 AB-g）的特殊之处在于它是不对称故障，故障点会同时出现负序和零序分量。负序分量在发电机转子上产生二倍频脉动转矩，零序分量又和中性点接地方式直接相关。对于单机无穷大系统，两相接地时的正序等效网络需要把负序、零序网络并联后串入，故障期间电磁功率的下降程度介于三相短路和两相短路之间。

用 Simulink 做这个题目的好处是：你不需要手推这些序分量，瞬时值模型直接就能把不对称故障的电磁过程算出来。但你如果完全不理解序分量的含义，后面看结果时会非常困惑——为什么转速的振荡曲线里带着高频抖动？为什么非故障相电压会抬高？这些都不是模型错了，而是不对称故障的真实特征。

1.3 转子运动方程与转速变化的物理直觉

发电机转速变化的源头，是转子上的力矩不平衡。转子运动方程写成惯用形式是：

[ \frac{2H}{\omega_s} \frac{d^2\delta}{dt^2} = P_m - P_e - D\frac{d\delta}{dt} ]

其中 H 是惯性时间常数，单位是秒，代表转子的转动惯量在额定转速下储存的动能与机组额定容量之比；(\omega_s) 是同步电角速度；(\delta) 是转子相对角度；(P_m) 是机械功率；(P_e) 是电磁功率；D 是阻尼系数。

把生活经验拿来类比：你推一辆手推车，持续往前推（机械功率），地面阻止它后退，但中途有人突然在后面拽它一下（电磁功率骤降），手推车就会往前冲一段。发电机的转速波动就是这么来的——正常运行时 (P_m = P_e)，转速稳定在同步速；故障发生瞬间，电磁功率大幅下降，而机械功率还是原来那么大，多出来的能量全部转化为转子动能，转速必然上升。

故障切除后，电磁功率恢复，转子受到制动，转速开始回摆。如果回摆后能稳定在同步速附近，系统暂态稳定；如果第一摆都没收住，转速持续偏离，就是失稳。Simulink 里观察的“发电机转速变化”，本质就是把这条加速、减速、再振荡的曲线完整记录下来。

2. Simulink建模实操：从同步发电机到故障模块的参数选择口径

2.1 元件级建模套路与推荐参数

Simulink 里做电力系统机电暂态仿真，走的是 Simscape > Electrical > Specialized Power Systems 这套电气专业库。模型链条非常清晰：原动机给机械功率，发电机把机械功率转成电功率，变压器升压，双回线送电，最后接无穷大母线。

我建议初学者直接采用一组可复制的典型参数，而不是一上来就纠结于精确复刻某台真实机组。下面这组参数我实际用过，跑两相接地短路非常顺：

为什么发电机用标幺值（pu）模型？因为电力系统分析的惯例就是标幺制，你看文献、对比结果、算等面积定则时全是 pu 值，统一口径后你可以把自己仿真里读到的功率、转速、功角直接代进公式去验证。如果选 SI 单位模型，电机内部一堆物理单位换算，很容易把励磁电压、转动惯量的量纲搞混。

原动机那边有个细节容易踩坑：同步电机模块的 Pm 输入端口期望的是 pu 值（标幺值模型下），你把常数模块设为 0.85 就行。但很多人不知道 0.85 这个值意味着发电机在额定容量下输出约 85% 的有功，这样发电机和无穷大系统之间有明确的功率传输，功角才会稳定在某一个非零角度上。如果你设成 0，那整个系统空载，后面故障时转速变化的表现和带载时完全不同，题目就失真了。

2.2 故障模块的时序编排与故障相选择

两相接地短路在 Simulink 里怎么设置，是很多人第一次卡住的地方。我以 Three-Phase Fault 模块为例讲透。

这个模块本质上是三组可控开关加一个接地开关的组合。老版本里你在参数面板上勾选 Phase A Fault、Phase B Fault，再勾选 Ground Fault，等于把 A、B 两相分别与地短接，这就是标准的两相接地短路。新版本里通常有 Fault Type 下拉选择，你直接选 Two-Phase-to-Ground（或者 Phase A-B to Ground），效果相同。

故障时序靠两个 Transition Time 控制：Transition Time 1 是故障开始时刻，Transition Time 2 是故障结束时刻。假设仿真从 0 秒开始，想让系统先在稳态下运行 10 秒，然后故障持续 0.2 秒，就设 Transition Time 1 = 10，Transition Time 2 = 10.2。

这个“先跑一段稳态再故障”的做法极其重要。因为电力系统仿真需要初始状态收敛，你一上来就故障，发电机内部的初始磁场、功角关系都还没建立，仿真结果基本是乱跳的，转速初始偏差可能都是几个 pu，完全不可用。

如果想做得更精细，可以用外部断路器模块组合代替 Three-Phase Fault 的简单时序。比如用三个单相断路器分别控制 A 相对地、B 相对地、以及中性点接地，再配合 Step 信号手动编排动作顺序。这样你能模拟出“A 相先接地，一个周波后 B 相也接地”这种更真实的故障发展过程。但对于课程作业级别，直接在 Three-Phase Fault 里设两个时间点就够了。

2.3 测量与数据导出设置

转速信号怎么取？同步电机模块的测量端口通过 Machines Measurement Demux 或者 Bus Selector 拉出来就行。不同版本、不同电机模型输出的信号名略有差异，但重点看两个：一个是机械转速 wm，另一个是电磁功率 Pe，有条件的话再把电机的 rotor angle 拉出来。

这里有个单位陷阱要特别提醒：在 pu Fundamental 模型下，wm 输出是标幺值，同步转速就是 1；在 SI Fundamental 模型下，wm 输出是 rad/s，同步转速取决于极对数和系统频率（2 极对 60 Hz 电机是 376.99 rad/s 电角速度，机械角速度是 188.5 rad/s，不同极对数还要再折）。很多人看到转速曲线在 376 附近跳，以为发电机转速是同步速的 376 倍，其实就是没注意单位。我自己的习惯是统一用 pu 标幺值，判稳直接看是否回到 1.0，省去换算。

Scope 适合人眼观察，但要做定量分析一定得用 To Workspace 模块把 wm、Pe、delta 这三个信号导出成时间序列。后面算最大转速偏差、临界切除时间、画摇摆曲线包络，都靠导出来的数据。别偷懒，这一步后面一定会用到。

3. 仿真结果怎么读：发电机转速摇摆曲线的物理含义与判稳逻辑

3.1 典型的两相接地故障仿真结果长什么样

模型搭好、参数照上面设置后，运行仿真，观察 0 到 20 秒的转速波形，你会看到三个阶段：

第一阶段是 0 到 10 秒的稳态段，转速恒定为 1 pu，一条水平直线。这说明初始潮流正确，发电机机械功率与电磁功率平衡，系统平稳运行。

第二阶段是 10 秒故障瞬间：转速曲线突然以近乎线性的斜率开始上升。时间很短可能看不出斜率变化，但你把波形放大，发现上升速率和 (P_m - P_e) 成正比，这就是转子加速的直接证据。两相接地故障下电磁功率并没有完全变为零，所以转速上升速度会比三相短路慢一些，这条曲线的斜率本身就能反映故障严重程度。

第三阶段是 10.2 秒故障切除后：转速先是继续冲高一段，然后掉头向下，越过 1 pu 继续下探，形成典型的“第一摆”振荡。如果系统有足够的减速能力，后续振荡幅度逐步减小，转速最终回到 1 pu；如果减速能力不足，转速在第一摆就突破某个极限，之后不再回头，单调上升直到失步保护动作。

我把这个过程的四个关键时间点总结一下：

3.2 转速变化背后的电磁-机械能量博弈

为什么转速不是故障切除后立刻回落，而是还会惯性冲高一段？这是很多初学者看曲线时最容易疑惑的地方。

故障切除后，系统的电磁功率曲线并不是简单回到故障前的状态。因为保护动作通常意味着跳开一回故障线路，只剩单回线运行，等效电抗增大，功角-功率特性曲线的幅值会比原先低。如果此时功角已经摆到比较大的角度，而该角度对应的电磁功率仍然小于机械功率，转子就还是继续加速，只不过加速度变小了。直到功角继续增大到电磁功率大于机械功率，转子才开始减速。所以转速曲线的拐点，和功角位置是有严格对应关系的，不是随便哪里就掉头。

用等面积定则的语言说，故障期间转子积累的动能对应“加速面积”，故障切除后转子释放动能对应“减速面积”。系统稳定的充要条件，就是最大减速面积大于等于加速面积。Simulink 算出来的转速曲线，就是这块能量差额的动态展示——转速升得多高，取决于加速面积有多大；转速能不能拉回来，取决于减速面积够不够。

两相接地故障和三相短路在等面积定则上的差异一目了然：三相短路时 (P_e) 约等于 0，加速面积几乎是最大；两相接地时 (P_e) 虽然减少但不为 0，加速面积更小。所以相同的切除时间下，两相接地的临界切除时间更长，稳定裕度更大。这也是为什么工程上对不对称故障的稳定问题相对不那么保守。

3.3 转速之外还要看什么

题目只让你看转速变化，但真正做研究时，只盯着转速一个量是不够的。最少再拉两条曲线：电磁功率 (P_e) 和功角 (\delta)。

电磁功率曲线能帮你确认故障起止时刻的物理过程是否如预期：故障瞬间 (P_e) 有没有骤降，切除瞬间有没有跳变，振荡过程是否与转速相位对应。功角曲线则是暂态稳定最经典的判据——功角单调越过 180 度，基本就宣告失稳。

功角的获取比转速稍微麻烦一点，需要把电机的转子角信号与无穷大母线正序电压相位做差，一般可以用 PLL 模块锁相无穷大母线电压再和电机功角信号相减。如果你只是做课程设计，不追求这个细节，直接用转速曲线判稳也可以：如果故障切除后转速第一摆的最大值不超过某个阈值且能回摆，就认为稳定。但写论文时建议把功角曲线也补上，审阅老师更认可这个。

4. 关键参数影响分析：故障切除时间、故障位置与系统惯量的扰动实验

4.1 临界切除时间的寻优方法

暂态稳定研究里最有工程价值的参数之一就是临界切除时间（Critical Clearing Time，CCT）。它的定义是：系统能保持暂态稳定的最晚故障切除时刻，超过这个时间再切除故障，系统就失稳。

仿真中找 CCT 的方法非常机械但有效——反复改 Transition Time 2，跑一系列仿真，观察转速曲线是否收敛。比如先用粗扫描，步长取 0.05 秒：10.05 秒切除稳定，10.10 秒切除稳定，10.15 秒切除失稳，那 CCT 就落在 10.10 到 10.15 之间。再在这个区间里用 0.005 秒的细扫描，逼近精确值。

做这个扫描时有个实用技巧：不要每次都打开 Scope 肉眼判断，写一个简单的 MATLAB 脚本批量跑仿真，然后读取 To Workspace 导出的 wm 序列，用程序判断失稳条件。判据可以定为：仿真结束时刻转速相对同步速的偏差超过某个阈值（比如 0.05 pu），或者转速在第一个长窗口内持续单调增长超过 0.3 pu。这样几轮下来，CCT 就能精确到小数点后两位。

我在自己的仿真里试过，上面那组参数下两相接地短路的 CCT 大概在 0.12 到 0.16 秒区间，具体取决于你设置的线路阻抗和发电机出力。你如果自己搭的模型参数跟我不同，CCT 差异大很正常，关键是掌握这个“粗扫定位区间、细扫逼近精度”的做法。

4.2 故障位置与过渡电阻的敏感性分析

同一个故障类型，故障点离发电机端越近，对稳定性的破坏越大。这个直觉你完全可以通过仿真验证：把 Three-Phase Fault 模块从线路末端挪到升压变压器出口，你会看到同样 0.2 秒的故障持续时间，机端故障时转速第一摆的峰值明显更高。原因很简单，故障点越靠近机端，机端电压被压得越低，电磁功率下降越狠，加速面积越大。

工程上更关注的是故障过渡电阻的影响。真实的短路都是通过电弧接地的，电弧电阻不是零，会限制故障电流。你在故障模块里可以设置 Fault Resistance，从 0 欧姆到几十欧姆扫一遍，转速上升幅度和 CCT 都会相应变化。这个结果导向一个工程结论：永久性故障的过渡电阻比较大时，系统承受的能力反而强一些，这也是继电保护整定中需要考虑的因素。

我建议同学们做这个敏感性分析时，一次性跑多组仿真，把不同故障位置、不同过渡电阻下的转速第一摆峰值汇总成一张表，这种对比在报告里非常有说服力。

4.3 惯性时间常数 H 对转速摆幅的影响

发电机组的惯性时间常数 H 是暂态稳定的另一个关键参数。物理上 H 越大，转子储存的动能越多，同样功率差额下角速度变化越慢。仿真实验很简单：把同步电机模块的 H 从 4.2 改成 2.5 和 6.0，分别跑一次相同故障，对比转速曲线。

结果你会发现：H 小的机组转速上升更快，第一摆峰值更大，看起来更容易失稳；H 大的机组转速摆动更平缓，振荡周期也更长。这也对应一个工程直觉——大惯量机组对暂态冲击更“抗造”。但注意，H 越大振荡周期越长，如果负阻尼或者弱阻尼存在，长时间慢速摇摆反而可能激发低频振荡，这就要靠 PSS 这类附加控制来处理了。

做这类参数扫描时，我强烈建议把每次仿真的机械功率、故障时序完全固定，只动 H 一个变量，保证是单变量对比。很多人喜欢同时改好几个参数，结果曲线对不上，自己也说不清是哪个参数起的作用，这是仿真分析的大忌。

5. 仿真中容易翻车的细节：求解器、标幺制、初始化与数值陷阱

5.1 求解器选择与容差设置

电力系统故障仿真属于典型的刚性系统：故障瞬间理想开关动作，电流电压波形里出现几乎瞬态的跳变，而机电暂态部分又是几百毫秒到几秒的慢过程。两种时间尺度耦合在一起，对数值积分器要求很高。

我的经验是首选变步长 ode23tb，它对这种刚性切换问题很稳，ode15s 在部分版本里也表现不错。定步长求解器在这种场景下不建议用，除非你步长取得非常小，否则故障瞬间的尖峰会被严重平滑，转速曲线的第一摆细节失真。

相对容差默认值是 1e-3，粗扫参数时够用；但当你需要精确读 CCT 或者做定量对比时，把相对容差调到 1e-4 到 1e-5，仿真结果会更平滑。代价是仿真速度明显变慢，故障瞬间步长会缩到微秒级。我的习惯是粗扫用 1e-3 快速定位，细扫用 1e-5 精确确认，两边兼顾。

5.2 初始化失败与功率潮流校验

Simulink 电力系统仿真最常见的报错，就是初始化不收敛。这个报错背后九成是初始潮流没建立起来，而非模型本身连接错误。

Standard 做法是双击 Powergui 模块，打开 Machines Initialization 或者 Load Flow，让软件根据你设定的发电机出力、机端电压和线路潮流自动计算初始功角和发电机内部状态。如果你在 Powergui 里看到潮流求解失败，优先检查三个地方：一是无穷大电源的内阻和接地方式是否正确，很多标准无穷大模型要求内部接线为 Yg；二是变压器接线组别与线路连接是否一致，高压侧 230 kV 是线电压，别用相电压去算标幺值基准；三是功率基准是否全系统统一，发电机 100 MVA、变压器 100 MVA 要一致，如果变压器是 50 MVA 而发电机 100 MVA，标幺值折算会乱套。

另一个更简单的办法是“空转校验”：先把故障模块禁用，让系统在无扰动下仿真 20 秒，确认转速始终为 1 pu、功率表读数与你设定的 0.85 pu 一致，再开故障。如果空转都稳不住，说明模型里某个参数或连接有问题，先排查掉再加故障，问题定位会快很多。

5.3 转速振荡不衰减是不是模型错了

很多同学做完仿真，发现故障切除后转速振荡完全不衰减，像一条永远不安静的水平弹簧，于是怀疑自己搭错了。其实这很可能是正确的物理结果，而不是模型错误。

同步电机模块的机械阻尼系数 D 默认是 0，转子运动方程里没有任何阻尼项，电磁暂态在实用模型中又被忽略，系统呈现的是无阻尼等幅振荡。真实机组有阻尼绕组、励磁系统、调速器、电网负荷特性提供的各种正阻尼，振荡才会逐渐衰减。你如果想看到更真实的衰减效果，可以在电机参数里给 D 设一个合理的正值（比如 1 到 2 pu 秒），或者给励磁系统加一个简化模型。这也是为什么进阶版的课程设计会要求你加入励磁系统甚至 PSS，不单纯是摆着好看，而是没有它们，摇摆曲线就不具备真实的收敛性。

另外，故障模块里理想开关的导通电阻默认为毫欧级，很小的值在数值上容易引起步长骤降。如果仿真速度非常慢，试着把故障电阻从 1e-4 欧提到 1e-3 欧，对转速这种机电量几乎没有影响，但数值稳定性会好很多。

6. 从两相接地仿真走向更实际的工程问题——扩展方向

6.1 把控制措施加进模型：强励与 PSS 的效果对比

做完基础的两相接地短路仿真，这个模型的使命还没结束。你可以在同步电机模块的励磁端口接上简化的励磁系统模型，观察故障期间快速强励对角速度摆动的抑制作用。强励的本质是在机端电压跌落时迅速提高励磁电压，维持发电机内电势，让电磁功率不要跌太狠，从而缩小加速面积。

仿真对比思路清晰：同一组故障参数下，分别跑“恒定励磁”和“带强励”两组，把两条转速曲线画在同一张图里。你会看到带强励的机组转速第一摆峰值明显更低，回摆更快。更进一步，在励磁系统输出后叠加一个 PSS 稳定器信号，还能观察到后续振荡被进一步压制。这三个层次的对比，放在研究报告里就是一个完整的“从无控制到有控制”提升暂态稳定的故事线。

6.2 从单机到多机的扩展思路

单机无穷大系统的天花板在于，它只有一台发电机，无法研究机群之间的相对摇摆。把无穷大母线替换成一台容量更大的等值发电机，或者干脆构建两台容量相近的发电机并联送电，就能看到两机之间的功角相对摆动问题。这是一个自然的升级路径。

Simulink 里多机建模比单机复杂的地方在于初始潮流分配：Powergui 的潮流计算需要你给每台机设定出力，并且要保证所有发电机的初始相角满足网络潮流方程。初学者很容易在这里卡住。我的建议是先从“两台机 + 公共母线”这种最简结构起步，确认空转稳定后再加双回线和故障模块。多机系统的故障分析不再只有一个转速曲线，而是每个电机的转速和相对功角，判稳标准也变成“任意两个发电机之间的相对功角是否发散”。

6.3 仿真与理论计算的交叉验证习惯

最后聊一个长期受益的习惯：每次仿真跑完，拿结果和等面积定则的理论估算对一下。比如用功角特性曲线的幅值算出加速面积和减速面积，比较是否临界；用转子运动方程估算故障期间的转速上升斜率，再和仿真曲线对比。这种交叉验证不需要多精确，对得上数量级和趋势，就说明模型设置没有大偏差。

我见过太多同学仿真跑得出来，但问他“为什么转速上升斜率是这么大”就答不上来。能做出结果和能解释结果是两码事，后者才是做仿真研究的根本目的。两相接地故障这个题目特别好，因为故障期间的电磁功率不是简单的零，你正好可以用序分量原理估算出故障时的功率特性，再和仿真曲线比对——如果理论预判和仿真差距过大，要么参数设置有问题，要么对故障类型的理解有偏差，排查的过程本身就是很大的收获。

我在实际做这套仿真时最深的体会是：不要急着把故障时间设得过于苛刻，先把无故障稳态跑平整，把转速单位、功率基准这些基础问题确认清楚，再一步步加复杂度。每一次改动只动一个变量，记录曲线变化，你对系统的理解会比刷十次仿真都深。这样一套流程走下来，两相接地短路转速变化这个题目，你不仅完成了仿真作业，实际也把电力系统暂态稳定分析的看家本领练了一遍。