基于MATLAB/Simulink的单机无穷大系统暂态稳定仿真分析

搞电力系统稳定性的，第一个正经模型基本都是从单机无穷大系统开始的。这模型看着简单——一台发电机、一段线路、一个无穷大母线，但真要在 MATLAB/Simulink 里把它搭起来，再把不同故障下的暂态过程仿真明白，涉及的细节一点都不少。我自己第一次做这个题目的时候，光是在故障模块的切除时间设置上就卡了大半天，后来又把等面积定则和功角曲线反复对照，才算真正把暂态稳定这回事串通了。这篇文章就以“暂态稳定性仿真分析·基于 MATLAB/Simulink 的单机无穷大系统”为核心，把模型搭建、故障设置、结果分析这三条线完整展开。适合正在做电力系统课程设计、毕业设计，或者刚开始接触 Simulink 电力系统仿真模块的同学参考。

1. 暂态稳定的本质：先搞懂它到底在分析什么

1.1 大扰动下的稳定问题

暂态稳定研究的是系统在遭受严重扰动之后的动态行为，比如线路发生短路、发电机突然跳闸、大容量负荷突然投切等。这类扰动跟小扰动完全不是一回事，小扰动关心的是系统在某个平衡点附近能不能通过线性化分析判断收敛，而暂态稳定要回答的问题更加直接：故障清除之后，发电机还能不能继续跟系统同步运行。

你可以这么理解：发电机转子上同时作用着原动机的机械转矩和发电机的电磁制动转矩。正常运行时，这两个转矩平衡，转速稳定在同步转速。故障一发生，由于电压跌落或短路导致功率传输受阻，电磁转矩会急剧下降，但机械转矩基本不变，转子就获得了一个净加速转矩，转速开始上升。故障切除后，电磁功率恢复，但转子已经带着多出来的动能往前摆，如果这个摆动的角度超过了某个极限，系统就没法回到同步状态了。

在数值上，这个“角度”就是发电机的功角，即发电机转子位置与无穷大母线参考位置之间的夹角。暂态稳定分析的大多数工作，最终都归结为观察功角随时间的运动曲线，也就是所谓摇摆曲线。功角能越过峰值后回落并逐步衰减振荡，系统稳定；功角一旦超过极限还继续单调增大，就说明发电机跟系统失去同步了。这就是为什么做这个仿真，第一件事不是急着搭模型，而是先把“功角”这个概念吃透。

1.2 转子运动方程与等面积定则

转子运动方程，也叫摇摆方程，是整个暂态稳定分析的核心算法基础。它描述的是发电机转子在机械转矩和电磁转矩不平衡时的运动规律，最常见的标幺值形式可以写成：

T_J · d²δ/dt² = P_m - P_e

其中 T_J 是发电机组的惯性时间常数，δ 是功角，P_m 是原动机机械功率，P_e 是发电机电磁功率。这个方程的物理含义很直白：只要机械功率和电磁功率不相等，转子就会加速或减速，功角就会变化。要是 P_m 一直大于 P_e，功角就会一直增大，系统最终失步。

H 这个惯性时间常数值得多聊几句。它代表了发电机转子储存动能的能力，H 越大，转子越“沉”，同样的功率差下加速越慢，系统越不容易失稳。这也是为什么在做实际工程分析时，机组的转动惯量是一个十分重要的参数。课程设计里如果导师给了你不同容量的机组数据，你可以对比看看 H 值对临界切除时间的影响，那结果非常直观。

等面积定则是从转子运动方程推导出来的一种图解分析方法。它的思路是：加速过程中转子的动能增量对应于加速面积，减速过程中动能的减少量对应于减速面积。如果故障切除后系统能够提供的减速面积大于或等于故障期间积累的加速面积，系统就能暂态稳定。临界情况就是加速面积刚好等于减速面积，这时对应的功角就是极限切除角。

用这个定则可以在不跑仿真的情况下，快速估算出系统能否在给定的切除角下保持稳定，也能帮助你理解为什么故障切除越快越稳定。它是在手工计算时代最核心的暂态稳定判据，现在虽然有了计算机仿真，但，这个物理图像依然很重要——你在 Simulink 里看到功角曲线回摆或者发散，本质上都可以回到等面积定则来解释。建议在跑仿真之前，先用等面积定则在纸上算一遍某个故障的极限切除角，然后跟仿真结果对比，你会对这套理论有完全不同的体会。

2. 单机无穷大系统在 Simulink 里怎么搭

2.1 为什么选单机无穷大模型

单机无穷大系统（Single Machine Infinite Bus，SMIB）是暂态稳定分析最经典的简化模型。它把复杂的电力系统等效成一台发电机通过线路接入一个容量无限大、电压和频率恒定的母线。说白了，无穷大母线就是一个“永远不会垮”的系统，它能吸收或提供任何功率，所以我们的关注点只需要集中在发电机这一台设备上。

现实中几乎没有真正的无穷大母线，但当一个局部系统接入一个大电网时，这个大电网相对于单台机组来说，确实有类似无穷大的性质。用这个模型做暂态稳定分析，好处是能把问题压缩到最小但仍不失物理本质的维度上。功角、电磁功率、故障清除时间这些因素之间的因果关系一目了然。如果一上来就搭一个多机系统，那光潮流计算、发电机之间的摇摆耦合就够你折腾的了，反而看不清核心机制。

另外，单机无穷大系统也是教材上推导等面积定则时的标准场景，软件算出来的结果跟理论推导能严格对上。课程设计或者入门阶段，用它作为仿真对象是最合理的选型。等以后需要研究机组间的暂态振荡或者互联电网稳定性，再扩展成多机模型也不迟。

2.2 模型分步搭建与参数设置

在 Simulink 里搭 SMIB 系统，推荐用 Specialized Power Systems 库下的电力系统专业模块。模型整体结构从左到右是：同步发电机、升压变压器、输电线路、无穷大母线。故障模块放在输电线路上。

需要用到的模块主要就是这些：

搭建步骤按这个顺序走：

第一步，把同步发电机模块拖进来。双击打开，机器类型选“Synchronous Machine（pu units）”，额定功率可以设为 100 MVA，额定电压 13.8 kV，频率 50 Hz。在“Parameters”页签里设置 Xd'、Xq'、H 等电感和惯性参数。注意初始条件标签里的功率和电压要跟后面潮流状态匹配，否则初始化会报错，这个后面我会单独说。

第二步，接变压器。发电机出口电压一般比较低，通过升压变压器升到输电电压，比如 230 kV。打开变压器模块，设容量、频率，原副边电压填 13.8/230 kV，漏抗选 0.1 pu 左右的典型值。接线时把发电机输出接到变压器原边，副边接输电线路。

第三步，接输电线路。线路用 Distributed Parameters Line，长度可以设 200 km，正序电阻 0.02 Ω/km，正序电感 0.8 mH/km，正序电容 0.012 μF/km，零序参数可以按比例设成 3 倍正序。线路长度和单位阻抗直接决定了发电机能向无穷大母线输送多少功率，也就直接影响暂态稳定裕度。

第四步，把无穷大母线模块接在线路末端。Three-Phase Programmable Voltage Source 内部的基波电压设为 230 kV，相角 0 度作为参考轴，频率 50 Hz。这里有一个关键细节：这个电压源的内部阻抗要设得比较小，比如 0.001 Ω 左右，否则它就不是“无穷大”了，会对仿真结果产生额外的影响。

第五步，在输电线路的适当位置插入 Three-Phase Fault 模块。建议把线路分成两段，故障点放中间。比如把 200 km 线路拆成两段 100 km，中间用故障模块连接，这样可以方便地改变故障位置的比例。

第六步，加入测量模块。发电机模块的 m 输出端口接到一个 Bus Selector，可以提取转子角、转速这些内部状态量。线路末端用 Three-Phase V-I Measurement 测量电压电流，方便后续算功率。

模型搭好之后，建议先把故障模块的切换时间设成无穷大（也就是不触发），先做一次稳态仿真，确认发电机能正常输出设定功率。这一步能帮你把建模和初始化问题跟故障仿真问题隔离开来，排查起来会快很多。

2.3 功角测量与换算

功角是判断暂态稳定的核心量，所以在模型中一定要准确测出来。Synchronous Machine 模块的 m 端口输出的信号里，有转子角位移、转速、电磁功率等。具体哪个信号取决于你的机器类型，一般有一个“Rotor angle deviation”（标幺或弧度格式）的状态量。

但这里有个容易踩的坑：发电机转子角是相对于它自己初始位置的角度变化值，而暂态稳定分析里的功角是发电机 q 轴与系统参考电压矢量之间的夹角。直接拿转子角输出去画摇摆曲线，严格来说是不对的，需要做一个偏移换算。

实际操作中，简单可靠的做法是在模型里加一个常数模块，作为初始功角偏移量。你把稳态时发电机输出的功角值（可以从潮流计算结果或 Powergui 的 Load Flow 工具里读到）记下来，然后在测量通道里把转子角偏差加上这个初始功角，就得到了对无穷大母线的实际功角。这个换算看起来简单，但我不止一次看到有同学直接把偏差值当成真实功角，结果做出的曲线跟理论分析对不上。

做换算时建议直接在模型里用加法模块和增益模块实现，输出到总线后记录成角度制单位。方便归一化处理还是做成 rad，但显示的时候我喜欢转成度（deg），因为功角超过 180 度、90 度这些概念用度来说更直观。

3. 故障仿真：从故障类型到切除时间

3.1 Three-Phase Fault 模块配置

单机无穷大系统的暂态过程是由故障触发和切除来推动的，Three-Phase Fault 模块就是用来模拟这个过程的。这个模块支持多种故障类型，通过勾选“Fault A/G”、“Fault B/G”、“Fault C/G”三个复选框来组合。三个相都勾上就是三相短路，只勾 A 相就是 A 相接地短路，勾 A 相和 B 相就形成两相相间短路。

关键参数有两个。一个是 Fault resistance，一般设 0.001~0.01 欧姆，代表金属性短路的接地阻抗。如果你特意模拟带弧光电阻的短路，可以设个较大的值，比如 1 欧姆，这时故障电流和功率传输情况会有区别，可以作为扩展研究的方向。

另一个是 Transition Times，这是个数组，表示断路器的动作时间序列。典型设置是两个元素，比如 [0.1 0.25]，意思是 0.1 秒时刻合闸（故障发生），0.25 秒时刻跳闸（故障切除）。如果要在一次仿真里设置多个动作事件，比如重合闸过程，就需要更复杂的时间序列，不过那是后面的扩展话题。

这里有个细节很多人容易忽略：Transition Times 的第一项是故障发生时间，不是仿真开始时间。如果你想看正常运行一小段再故障，就必须把第一项设成大于零的数值，比如 0.1。同时，仿真总时长要覆盖故障切除后的振荡过程，一般是 3~5 秒，否则功角曲线的末段趋势看不清，没法判断稳定还是失稳。

3.2 不同故障类型的严重程度对比

不同故障类型对暂态稳定性的冲击差异非常明显，这也是这题目最有教学价值的部分之一。在同一个单机无穷大模型上，只改变故障模块的勾选组合，其他条件完全不变，就可以对比各种故障下的功角响应。

从故障期间电磁功率跌落的角度看，严重程度排序是这样的：

这里面最值得注意的是三相短路和单相接地之间的差距。实际运行中输电线路故障里单相接地占绝大多数，因为雷击、树枝碰线大多只影响一相。三相短路很少见，但一旦发生，基本是最恶劣的工况，工程上的校核往往都按三相短路来做。做仿真的时候可以每次只改故障类型，切换时间固定，观察功角最大摆幅的变化，你会直观地看到系统在单相故障下稳定裕度大得多。

另外，故障位置对严重程度的影响也很大。故障点离发电机越近，故障期间发电机端的等效电势到故障点的转移电抗越小，电磁功率跌落越严重，暂态稳定问题越严峻。在教学模型里，可以在不同位置各放一个故障模块，分别仿真，对比功角曲线。试验下来，故障发生在变压器高压侧出口处的效果，通常比线路末端严重得多。

3.3 切除时间扫描与临界切除时间

故障切除时间是暂态稳定的关键控制变量，也是继电保护动作速度的依据。同一故障，切除越快，系统越稳；切除慢了，功角摆过去之后就救不回来了。系统能保持稳定的最大故障切除时间，就是临界切除时间（Critical Clearing Time，CCT），这是暂态稳定分析最常用的量化指标。

手工用等面积定则求 CCT 要做不少计算，但有了 Simulink，可以很方便地做参数扫描。我习惯写一个 MATLAB 脚本，循环修改故障模块的 Transition Times，然后运行仿真，读取功角曲线的最大值，判断系统是否失稳。

下面这个脚本可以直接用，前提是模型中已经用日志方式把功角信号记录到了 logsout 里：

% 暂态稳定性参数扫描：估算临界切除时间 % 模型名：SMIB_Stability clear; clc; model = 'SMIB_Stability'; load_system(model); % 切除时间扫描范围 clear_times = 0.10:0.02:0.40; results = []; for i = 1:length(clear_times) t_clear = clear_times(i); % 故障0.1s发生，t_clear时刻切除 set_param([model '/Three-Phase Fault'], ... 'TransitionTimes', mat2str([0.1, t_clear])); % 运行仿真 simOut = sim(model); % 提取功角数据（记录信号名为 delta_deg） delta = simOut.logsout.get('delta_deg').Values.Data; delta_max = max(delta); results(i, :) = [t_clear, delta_max]; if delta_max < 180 fprintf('切除时间 %.2f s：稳定，最大功角 %.1f deg\n', ... t_clear, delta_max); else fprintf('切除时间 %.2f s：失稳，功角越过180度\n', t_clear); end end % 输出临界切除时间范围 stable_idx = find(results(:, 2) < 180); if ~isempty(stable_idx) last_stable = results(stable_idx(end), 1); fprintf('临界切除时间大约在 %.2f s 到 %.2f s 之间\n', ... last_stable, clear_times(stable_idx(end) + 1)); end

判断稳定与失稳的阈值，这里用了功角 180 度。严格来说，功角超过 180 度后系统还能通过减速恢复到同步，但常规的单机无穷大系统暂态稳定判据中，功角一旦超过 180 度，继续发展下去基本就是失步，所以用 180 度作为快速判据是工程上常用的做法。你要是想严谨一点，可以看功角曲线是否在仿真结束时已经回摆并收敛，脚本里也可以加一个判断，看最后 0.5 秒内功角的变化幅度是否小于某个值。

实测典型结果是这样的：一个 H=3.2s、双回线路 200km、满发的单机无穷大系统，三相短路发生在线路中点，切除时间 0.15 秒时功角最大约 120 度，稳定；切除时间 0.25 秒时功角已经超过 180 度，系统失稳。也就是说，这个工况的 CCT 差不多在 0.18~0.22 秒之间。你可以对比不同 H 值、不同故障位置下的结果，会发现 CCT 的变化规律完全符合等面积定则的预期。

4. 仿真结果怎么看：从曲线形态到稳定判据

4.1 功角摇摆曲线的判稳方法

仿真跑完，第一件事就是画功角随时间变化的摇摆曲线。稳定和失稳两种情况在曲线上表现得非常清晰，见多了之后一眼就能分辨。

稳定工况下，功角在故障切除后继续往前摆动一段时间，到达一个最大值后回摆，然后围绕一个新的稳态功角做衰减振荡。比如前面提到的切除时间 0.15 秒工况，功角从初始稳态值（比如 30 度）开始上升，故障期间慢慢加速，切除后继续冲高到 120 度左右，然后回弹，随后在 35 度附近振荡并逐渐衰减，这就是典型的暂态稳定响应。

失稳工况下，功角曲线则是单调上涨的形态，超过 90 度后还在加速，越过 180 度后继续增大。这时候代表发电机跟无穷大母线已经“脱钩”了，转速越拉越高，电磁功率振荡剧烈但无法建立同步。如果仿真时间再拉长，功角会继续一圈一圈地转动，就像两个齿轮脱开了一样。

判稳还有一个非常实用的技巧：看功角曲线回摆之后是不是形成了清晰的“包络衰减”。稳定系统在头 1~2 次振荡中衰减明显，如果功角振荡幅度长时间不衰减，甚至越来越大，那也可能是出现了振荡失稳或者弱阻尼问题。单机无穷大系统的机械阻尼相对较小，但你可以在发电机模块里设置阻尼系数，看看它对曲线的影响，这也是很好的扩展实验。

4.2 辅助判据：转速、电磁功率和电压

功角是主判据，但只看功角有时候不够，我还习惯同时看发电机转速偏差和电磁功率的波形。这三个量放在同一个时间轴上，相互印证，整个动态过程就会清楚很多。

转速偏差（Speed deviation）在故障瞬间会突然上跳，因为电磁转矩丢了，转子加速。切除故障后，转速偏差继续增大一段时间，直到功角摆到最大位置附近才开始回落。稳定工况下转速偏差会围绕零值振荡并逐渐收敛；失稳工况下转速偏差会持续为正，说明转子一直在超速，同步已经不可能维持。

电磁功率波形也能说明问题。故障发生时电磁功率突然跌落，故障切除时电磁功率又跳升到较高水平。这个“跳升”的幅度跟切除时的功角位置有关，功角摆得越高，切除后的电磁功率峰值反而可能越低，因为功率跟功角的正弦值有关。你结合等面积定则来看，就明白为什么切除慢了恢复能力变差。

电压波形主要是辅助确认故障模块是否按预期工作，以及断路器的开断是否正确。比如三相短路期间，故障点附近电压显著跌落甚至趋于零；切除后电压逐渐恢复。如果电压波形出现不正常的尖峰或者持续振荡，往往是数值计算问题或者参数设置不当，需要回到求解器设置里检查。

4.3 影响暂态稳定性的主要因素

既然已经能跑仿真，那接下来最值得做的事情就是做几个“单变量扫描”，把影响暂态稳定性的各因素挨个验证一遍，这也是这题目的核心价值之一。以下这几个因素对仿真结果影响最显著：

故障切除时间。这个不用多说，切除越快，加速面积越小，系统越稳定。CCT 是暂态稳定分析最有代表性的指标。

发电机惯性时间常数 H。H 越大，转子动能越大，同样的加速转矩下角加速度越小，CCT 越大，系统越稳。换个通俗说法，就是“重车惯性大，不容易被推倒”。你可以把 H 从 3 秒调到 6 秒，看看 CCT 能提升多少，非常直观。

发电机暂态电抗 Xd'。Xd' 越小，故障期间和故障切除后的功率传输能力越强，稳定裕度越大。这可能有点反直觉——电抗不是越小越容易被短路吗？但暂态稳定关心的不是短路电流大小，而是故障期间和故障后发电机与系统之间能传递多少功率，所以 Xd' 小反而有优势。

系统电抗（变压器漏抗、线路电抗）。线路越长、电抗越大，发电机与无穷大母线之间的传输极限越低，初始功角越大，稳定裕度也越小。你把线路从 100 公里加到 400 公里，CCT 会明显下降。

出力水平 P_m。发电机带得越重，初始功角越大，同样故障下越容易失稳。这也是为什么发电机在接近满发时，对继电保护动作时间的要求更高。

做这类扫描实验时，记得一次只改一个参数，否则几个因素混在一起，你就无法判断到底是哪个变量引起的效果了。记录数据时，把 CCT 和最大功角整理成表格，你会发现规律非常清晰。

5. 常见问题与调试经验

5.1 求解器与仿真步长：仿真跑得快慢全是它决定的

电力系统 Simulink 模型里包含同步电机、变压器、分布参数线路，还有一个开关性质的故障模块，这个系统在数值上属于刚性系统。刚性系统用普通的 ode45 求解器会跑得极其慢，甚至出现数值发散。我建议直接用 ode23tb 或 ode15s，这两种是专门处理刚性问题的隐式求解器，速度和解的可靠性都靠谱得多。

最大步长参数也值得单独说。默认的 auto 可能会让故障时刻附近的事件检测不够精细，导致波形出现高频毛刺。实践中我通常把最大步长设成 1e-3 或 5e-4 秒，仿真时间 3 秒也就几千步，完全跑得动。不建议设成 1e-5 以下，那会让仿真慢到你怀疑人生。

还有一个小坑：故障模块的动作时间设定值如果在 ode23tb 的步长网格里没被精确检测到，波形会出现跳变或者丢失事件。Simulink 对这类开关事件有零交叉检测机制，但前提是你没有关闭它。保持零交叉检测默认开启，然后通过最大步长保证时间精度，故障切除时刻的判断就会很准确。

5.2 初始条件是仿真能跑起来的大前提

关于同步发电机模块的初始条件，这可以说是新手翻车最多的地方。Synchronous Machine 模块内部自带初始功率和电压设置，如果跟外部线路和无穷大母线算出的潮流状态不一致，初始化的时候就报错，提示“Initialization commands cannot be evaluated”之类。

解决办法是先把模型里的 Powergui 模块找出来，双击打开，选择“Load Flow”工具，设置发电机的有功、无功输出以及母线电压幅值，然后让它自动计算整个系统的初始潮流。计算完成后再点击“Machines”按钮，把发电机模块的初始状态量写回模型。之后运行仿真，发电机就能直接在稳态平衡点启动，而不需要经历一个过渡过程或者直接初始化失败。

我在搭模型时有个习惯：先把故障模块的 Transition Times 设成很大的值，比如 [100 101]，然后做一次纯稳态仿真。此时功角应该是一条平坦的直线，转速偏差接近 0，电磁功率等于设定值。确认稳态没问题之后，再设置真实的故障和切除时间。这一步能帮助你把“模型本身有问题”和“故障动态不符合预期”两种情况区分清楚。

5.3 常见报错与排查技巧速查表

以下是这模型在实际调试中我最常遇到的一批问题，整理成速查表，碰到了直接对照处理：

排查问题有一个思路上的建议：先从稳态部分查起。把故障功能临时关掉，看系统能不能稳定运行。这一步没问题，再进入故障仿真。故障仿真出了问题，先看故障发生前的那一段波形是否正常，如果故障前就已经在振荡，说明还是稳态初始化的问题。如果故障前正常、故障中异常，再怀疑故障模块和数值设置。这样一层一层剥开，定位问题很快。

还有一个经验，做批量参数扫描时不要用sim命令连续运行同一个模型而每次都加载一遍模型，那样又慢又容易因为内存累积出错。正确做法是先用load_system加载一次模型，然后在循环里用sim配合set_param修改变量。仿真结束后用bdclose统一关闭模型。脚本跑批量扫描并自动输出报告，可以省下大量手动操作。

关于把功角数据拿出来画图，我习惯用simOut.logsout.get('delta_deg').Values.Data这个接口，前提是模型中了启用 Signal Logging 并把功角信号标记为记录信号。你也可以用 To Workspace 模块，简单直接，但批量扫描时 logsout 更规范，而且不会在多次仿真间产生变量覆盖的问题。

最后再分享一点个人体会。这类仿真分析真正让你受益的，往往不是最后那张漂亮的功角曲线，而是你在反复调参、对比、排查的过程中建立起来的物理直觉。当我第一次跑出稳定和失稳两条曲线放在一起对比时，等面积定则里那些抽象的“加速面积”和“减速面积”一下子就具体化了。后续你做多机系统或者含新能源的稳定性分析，很可能仍然要回到单机无穷大模型来找感觉。所以这题值得沉下心把它吃透，现在多花点时间搞清楚“故障—切除—稳定”这个链条的每个环节，后面看任何暂态分析文献都会顺手很多。