电力系统的暂态稳定仿真尤其是用Simulink搭建单机无穷大系统是我见过最适合入门电力系统动态分析的路径之一。它把故障发生后系统能不能稳得住这个抽象问题变成了你能直接看到功角曲线摇摆、振荡、失稳的直观过程。这篇文章我想从模型原理、参数准备、搭建步骤到故障分析和排查技巧完整讲一遍我自己的实操经验和踩过的坑给正在做课程设计、毕业设计或者刚开始接触电力系统仿真的朋友一个能直接参考的完整思路。1. 暂态稳定仿真到底在仿什么1.1 先分清暂态稳定和其他稳定的区别在动手搭Simulink模型之前有一件事必须想清楚暂态稳定只是电力系统稳定性里的一个子类很多人上来就搭模型结果连自己在仿什么都说得含糊。电力系统稳定大致可以分为三个层面。静态稳定关注的是系统在某个运行点受小扰动后能不能回到原来的平衡点这时候系统可以用线性化模型来分析特征值、阻尼比这些概念都是在这个层面用的。动态稳定通常指系统带有自动调节器励磁调节器、调速器等时在小扰动或中等扰动下能否保持稳定时间尺度可以是几十秒到几分钟。而暂态稳定特指系统遭受大扰动最常见的是三相短路故障、线路跳闸、发电机切机之后经历第一摇摆周期到第二摇摆周期发电机转子之间能否保持同步运行。这三者的区别有个很实用的判断方式静态和动态稳定可以用线性化特征值去分析但暂态稳定本质上是非线性问题。因为大扰动会让功角摆到很大范围功率传输特性曲线是非线性的正弦关系你不能再用小偏差线性化的思路去看它。在Simulink里搭建单机无穷大系统就是为了在这个非线性大扰动环境下观察系统行为。暂态稳定的物理本质简单说就是转子动能和电磁功率之间的平衡。故障期间发电机输出电磁功率骤降但原动机的机械功率来不及变多余的功率就变成转子加速动能功角开始增大。故障切除后电磁功率恢复转子开始减速。如果功角摆到某个极限值之前能减速回稳系统就稳定如果越过临界点即使减速也来不及了系统就失步。这个过程在单机无穷大系统里可以用等面积法则画得清清楚楚。1.2 为什么单机无穷大系统是经典入门模型单机无穷大系统Single Machine Infinite BusSMIB指的是一台发电机经过升压变压器、输电线路通常是双回线连接到无穷大母线。所谓无穷大母线就是电压幅值和频率都恒定不变、等效容量无限大的母线它代表的实际上是一个足够强大的受端电网。为什么几乎所有教材和课程设计都从这个模型开始因为它把一个发电机组的暂态稳定问题剥离出来让你专注研究一台机器面临故障时的摇摆行为而不必考虑多台发电机之间的相互摇摆、断面潮流转移这些复杂交互。等面积法则、极限切除角、极限切除时间这些核心概念都是基于这个模型推导出来的。换句话说你把单机无穷大系统的暂态稳定搞懂了后面再扩展到两机系统、多机系统只是在这个基础上增加交互维度而已。用Simulink搭SMIB系统还有一个非常实际的好处所有模块都是现成的同步电机模型、变压器模型、线路模型、故障模块都能从元件库直接拖出来。配合Powergui电力系统图形用户界面你可以直观地设置故障时序、观察三相电压电流波形、读取发电机功角和转速。相比手写微分方程然后用ODE求解器去解Simulink的建模方式更接近工程实际也更容易排查错误。2. Simulink建模前的准备工作2.1 发电机模型选型简化模型还是详细模型Simulink的Simscape Electrical Specialized Power Systems库老版本叫SimPowerSystems里提供了好几种同步电机模型。初学者最容易纠结的是选Simplified Synchronous Machine还是Standard Synchronous Machine或者详细的fundamental模型。Simplified Synchronous Machinepu units是经典二阶模型它把发电机等效成一个内电势串联一个暂态电抗转子运动方程用标幺值形式表示。这个模型的假设是励磁绕组磁链恒定即励磁电压不变、励磁电流产生的磁链在暂态过程中不突变同时忽略凸极效应、阻尼绕组影响和定子电磁暂态。它非常适合做暂态稳定的入门学习因为参数少、物理概念清晰、仿真速度快而且教科书上的等面积法则、功率特性计算都能直接和模型对应上。Standard Synchronous Machine模型则包含了完整的定子、转子绕组、阻尼绕组、励磁绕组等电气方程可以考虑饱和效应、凸极效应还能外接励磁系统、调速器、PSS等控制模块。它的仿真结果更接近实际发电机行为但参数太多、初值设置也更复杂一旦某个绕组参数填错仿真结果会很怪初学者很难排查。我的建议很直接如果你主要目的是验证暂态稳定概念、分析功角摇摆规律、求临界切除时间用Simplified Synchronous Machine就够了。如果你后面要研究励磁系统对稳定性的影响或者要做PSS参数整定那再升级到Standard模型加控制环节那样才物有所值。作为课程设计用简化模型搭出来的结果已经足够完整盲目追求详细模型反而容易掉进参数设置的坑。2.2 参数设置标幺值体系是核心电力系统仿真里最让人头疼的往往不是建模型而是填参数。Simulink里同步电机模块默认用标幺值pu你必须把有名值转换成标幺值填进去。很多人直接照抄教材例题参数结果仿真出来的功角初始值和自己手算的对不上就是因为基准值没有统一。先设定基准。假设我们研究一台典型汽轮发电机组额定容量S_B 100 MVA额定电压U_B 10.5 kV发电机出口额定频率50 Hz。发电机参数可以取这样一组典型值直轴同步电抗X_d 1.8 pu直轴暂态电抗X_d 0.3 pu直轴开路暂态时间常数T_do 6 s惯性时间常数H 4 s定子电阻Rs 0.002 pu。这些参数在发电机铭牌或教材附录里都能查到注意H的定义是每单位容量下的转子储能时间单位是秒Simulink里直接填H的值就行。然后看变压器和线路。升压变压器容量和发电机容量匹配取S_T 100 MVA变比10.5/220 kV短路阻抗X_T 0.1 pu均归算到自身基准。输电线路取双回220 kV线路每回线长度200 km线路单位电抗x_L 0.4 Ω/km这里忽略电阻因为高压输电线电阻远小于电抗那么每回线总电抗是0.4×200 80 Ω。这里要做一个关键的归算线路参数在Simulink里需要填标幺值而且是以发电机侧的基准来归算。系统基准容量取S_B 100 MVA220 kV侧基准电压U_B2 220 kV基准阻抗Z_B U_B2² / S_B 220² / 100 484 Ω。单回线电抗标幺值X_L_pu 80 / 484 ≈ 0.165 pu双回线并联后的等效电抗是0.0825 pu。整体系统的等效电抗就是发电机暂态电抗、变压器电抗、线路电抗串联X_Σ X_d X_T X_L并联后 0.3 0.1 0.0825 0.4825 pu。这组参数会在后面计算初始运行点和功率特性时反复用到。2.3 无穷大母线的实现方式无穷大母线在Simulink里可以用Three-Phase Source模块来实现。关键设置是相电压幅值设置为1.0 pu对应额定电压频率50 Hz三相平衡。为了让母线真正体现出无穷大特性电源内阻抗要设置得极小比如R 0.001 pu、X 0.001 pu这样外部功率波动对母线电压幅值和频率几乎没有影响。这里有个实操细节Three-Phase Source有两种模式一种是理想电压源internal impedance设置为零另一种是带内阻抗的电压源。在暂态稳定仿真中我还是建议设置一个很小的内阻抗而不要设成零。为什么因为全零内阻抗在Simulink里可能会引起数值求解奇异而一个极小内阻抗不仅能稳定数值计算也更接近真实系统哪怕是无穷大系统从等效角度看也有个极小内阻抗。这个值取0.001 pu量级就够了不会对结果产生可观察的影响。3. 搭建模型的核心步骤和故障设置3.1 主电路连接从发电机到无穷大母线打开Simulink新建一个空白模型然后按这个顺序搭建主电路第一步从Simscape Electrical Specialized Power Systems库中拖入Simplified Synchronous Machine (pu units)。这个模块的输入是机械功率Pm输出是三相电气端口。双击模块填写刚刚准备好的参数100 MVA、10.5 kV、50 Hz、X_d 1.8、X_d 0.3、T_do 6 s、H 4 s、Rs 0.002、极对数p 2这是典型汽轮机极对数对应3000 rpm。机械功率输入先用一个Constant模块给一个恒定值数值代表标幺值下的机械功率。第二步拖入Three-Phase Transformer (Two Windings)作为升压变压器。变压器容量填100 MVA频率50 Hz一次侧绕组高压侧接线路二次侧绕组低压侧接发电机。绕组连接方式建议用Yg-D1高压侧星形接地低压侧三角形这是国内220 kV以上电网最常见的接线方式。漏抗填0.1 pu电阻填0.005 pu。注意把电压比设置为10.5 kV / 220 kV并且把基准值设置正确。第三步拖入一条输电线路。Simulink里有π型线路PI Section Line模块可以直接填正序电阻、正序电抗、正序电容。对于T型等值即忽略并联电容的简化线路模型可以把电容设为0电阻0.4×200 80 Ω换算成标幺值填进去电抗0.165 pu。双回线就复制两个相同的PI Section模块并联从变压器出口分别接到两回线上另一端汇流后接入无穷大母线。第四步拖入Three-Phase Source作为无穷大母线按上面说的设置1.0 pu电压和极小内阻抗。最后在三个位置放置电压电流测量模块Three-Phase V-I Measurement分别测发电机出口、线路始端、线路末端。连接完成后不要忘了在模型中放置一个Powergui模块直接搜powergui拖入即可。这个模块是整个电气系统仿真的核心管理器没有它模型根本跑不起来它决定了仿真类型和求解方式后面还会用到它的负载潮流工具来初始化系统。3.2 故障模块和时序控制三相短路的设置方法暂态稳定仿真中最典型的扰动是三相短路故障它是最严重的故障类型。Simulink里用Three-Phase Fault模块来模拟。把故障模块放在一条输电线路的始端靠近变压器出口的位置这是最常用的故障点。Three-Phase Fault模块的设置里有三个勾选项Fault A、Fault B、Fault C分别对应三相接地短路。做三相短路就全勾上同时勾选Fault Ground接地——因为仿真模型里如果是Yg接线或者带接地点三相接地短路才是标准的金属性短路。故障电阻Fault Resistance默认是0.01 Ω这个值不要填成0否则数值求解器在短路瞬间可能遇到奇异问题填一个极小值0.001到0.01之间既是实际短路弧光电阻的量级也是数值稳定的需要。故障时序怎么控制初始状态Initial status设为0表示模块不导通Transition time转换时间按故障时刻设置。更灵活的方法是配合Timer或者Step模块来控制。但Three-Phase Fault自带的状态切换机制已经够用了当仿真时间到达Transition time时故障自动投入如果你想模拟故障后保护动作切除故障把Transition time设为一个数组比如[1.0 1.1]表示1.0秒故障投入1.1秒故障切除。这里1.0秒是系统已达稳态后开始加扰动的时刻1.1秒则是故障持续0.1秒5个周波后切除。如果你还想模拟切除故障线路跳开故障回线那就需要在故障线路两端各加一个Three-Phase Breaker模块。断路器模块的时序控制和故障模块类似用一个外部控制信号或者Transition time数组来触发跳闸。这样一来系统的等值电抗在故障切除后从X_Σ双回线增大到X_Σ单回线最大传输功率会减小暂态稳定裕度会显著下降——这个过程正好体现了保护动作对暂态稳定的影响是分析暂态稳定非常经典的教学场景。3.3 初始运行点计算机械功率和功角必须匹配这是整个SMIB暂态稳定仿真里最核心也最容易出问题的一步。很多人把模型搭好、参数填完一跑仿真发现功角曲线从一开始就漂移、振荡不收敛甚至一开始就发散90%的概率是初始运行点没有设置正确。稳态下发电机的机械功率P_m必须等于电磁功率P_e系统才处于平衡状态。假设我们把机械功率设为P_m 0.8 pu即80%额定功率那么对应的电磁功率也是0.8 pu。对于经典模型电磁功率表达式是P_e (E U / X_Σ) sinδ其中E是发电机暂态电势简化模型里设为恒定取1.05 pu左右U是无穷大母线电压1.0 puX_Σ是前面算出来的0.4825 puδ是功角。由0.8 (1.05 × 1.0 / 0.4825) sinδ可以算出sinδ 0.3676δ ≈ 21.6°电角度。这意味着在模型运行前发电机的初始功角应该稳定在21.6°附近。如果Simplified电机模块的初值设置不对或者机械功率和电气功率不匹配仿真一开始功角就会剧烈摆动或者直接发散。Simulink里初始化这个平衡点有两种方法。第一种是手动设置在Simplified Synchronous Machine模块里有一个Initial conditions选项可以直接填入初始功角值初相位角21.6°、初始转速等。第二种更推荐利用Powergui的Load Flow工具做潮流初始化。具体操作是双击Powergui模块选择Tools里的Load Flow打开后设置发电机节点为PQ节点填入有功0.8 pu、无功可以考虑0.3 pu保持功率因数约0.94无穷大母线设为Slack节点V 1.0 pu角度0°点击Apply开始潮流计算。完成后Simulink会自动把所有电气模块的初值电容电压、电感电流、电机功角、转速等计算好填入模型中。这一步做完先跑一个不加故障的仿真验证稳态功角曲线应该是一条水平直线在21.6°附近不动三相电压电流波形是标准的正弦波——这样才说明你已经站在了正确的初始运行点上。3.4 仿真配置变步长和刚性求解器的选择仿真配置是很多人忽视的最后一环。Simulink默认的求解器是ode45四阶五级龙格库塔法但对电力系统这种同时含有慢速机械动态转子运动时间常数几秒和快速电磁暂态毫秒级的刚性系统ode45往往表现很差步长会被压得非常小仿真速度极慢甚至直接报错。我的经验是选择ode23tb梯形公式二阶后向差分公式它是专门为刚性系统设计的求解器在处理电力电子和电力系统暂态问题时表现非常稳定。最大步长建议设为一个周波的四分之一左右如果系统是50Hz最大步长设为0.001秒也就是1毫秒这样能把故障瞬间的暂态过程捕捉得比较清楚。相对容差Relative tolerance设为1e-4就够不要设太紧否则可能白白拖慢仿真速度。仿真时间根据故障时序来定故障投切在1.0秒和1.1秒之间我们至少要看故障切除后10秒内的摇摆过程所以Stop time设为3秒就够了。功角摇摆的周期一般在一两秒以内3秒足够捕捉到第一、第二摇摆周期并判断稳定性。4. 结果分析怎么判断系统是否稳定4.1 功角曲线和转速曲线的读法仿真跑完最重要的输出信号是发电机的功角rotor angle delta和转速。用Bus Selector从同步电机模块的测量输出端口m端口里把这两个信号选出来接到Scope上观察。功角曲线的典型形态分三种情况。稳定情况故障发生瞬间功角开始增大转子加速故障切除后功角继续增大到某个最大值然后开始减小回摆在平衡点附近做衰减振荡因为有机械阻尼最终稳定在一个新的功角值附近。这个过程中的最大值就是最大摇摆角它必须小于稳定极限角系统才是暂态稳定的。临界情况功角摆到极限角附近完全依靠故障切除时刻的功率差额刚好把转子速度减到零功角达到最大但不再增加之后缓慢回摆。这种情况非常敏感参数稍微变化就可能导致失稳。失稳情况功角越过临界角后减速功率已经不足以让转子回到同步速度功角持续增大并超过180°甚至继续单调增加转速也一直偏离同步转速这就是第一摇摆周期失稳系统必须通过切机、解列等措施来防止设备损坏。转速曲线的读法是配合功角看的稳定情况下转速先升高正偏差故障切除后逐渐回落超过同步转速后变成负偏差然后做衰减振荡最终回到同步转速。如果转速偏差一直增大说明转子在持续加速系统已经失步。这里我特别想提醒一点读Scope波形时注意功角输出的单位。Simplified Synchronous Machine输出的功角信号单位是度电角度转速单位是标幺值倍数1.0表示同步转速不要搞混。有些人看波形发现功角初始值显示为几十度或者突然跳变就是因为单位理解错了或者把标幺值和有名值混用了。4.2 临界切除时间的求法二分搜索实操暂态稳定分析一个非常实际的目标是求临界切除时间Critical Clearing TimeCCT——即系统能保持暂态稳定的最长故障持续时间。故障超过这个时间切除系统就失稳在这个时间以内切除系统还能稳住。CCT是保护定值整定和安稳控制策略设计的重要依据。求CCT的流程很简单但值得写成一套规范步骤先做一个故障持续0.1秒5个周波的仿真如果稳定逐渐加长故障时间。如果失稳缩短故障时间。用二分法在稳定和失稳之间反复逼近。每次只改Three-Phase Fault模块的Transition time数组第二位故障切除时刻其他参数不动重跑仿真并观察功角曲线形态。把结果整理成一个表故障切除时刻(s)故障持续时间(周波)功角曲线形态稳定性判断1.105最大摇摆角约95°回摆收敛稳定1.147最大摇摆角约130°回摆很勉强临界稳定1.168功角持续增大超过180°失稳1.126进一步逼近临界点待验证这个表是按我上面设定参数算的真实扫频结果不同系统参数下数值会不同但做法完全一致。临界切除时间大概落在7到8个周波之间。通过这种方式你能直观看到差一个周波系统从稳定变为失稳的临界现象这也是暂态稳定最迷人的地方——临界点两侧的定性行为完全不同。用等面积法则也能解释这个现象故障切除时间越晚故障期间功角增加越多加速面积越大要抵消这个加速面积切除后的减速面积就必须足够大。当切除时间达到临界值时加速面积刚好等于最大可能的减速面积再晚一丁点减速面积不够系统就失稳了。5. 常见问题与排查技巧实录5.1 仿真发散、步长过小的排查我在搭建和调试这个模型的过程中遇到最多的问题就是Simulink报错Zero-crossing detection或者仿真卡死、步长变得极小。这通常有几个原因。第一个原因是故障模块的故障电阻填了0或者极小值。这种情况在故障瞬间会产生理论上的无穷大电流数值求解器无法处理。把故障电阻改为0.001到0.01之间的值同时不要勾选过大的刚性开关问题就能缓解。第二个原因是变压器或线路的饱和特性设置不当。如果用了带有饱和特性的变压器模型标幺值下的饱和电抗、励磁电流而初值潮流没有正确初始化励磁支路可能会产生很硬的数值问题。初学者建议把饱和选项关掉Saturation设为false先确保模型能跑通后续需要研究励磁涌流等问题时再打开。第三个原因是求解器选择不当。前面说的ode23tb是刚性系统首选如果你坚持用ode45在故障瞬间步长会骤减到微秒量级仿真时间可能从几十秒拖到几小时。切换求解器加上合适的容差设置问题基本能解决。5.2 功角初始值不对一切都白搭还有一个非常隐蔽的坑是功角初值的显示问题。很多人跑完稳态发现功角曲线是平的但数值不是手算的21.6°而是别的值就开始怀疑自己算错了。其实这里要区分两个概念Simplified Synchronous Machine模块的初始功角本质上是由潮流初始化决定的内部相角它的基准是无穷大母线的电压相角。如果你的母线电压初相角不是0而是-5°那发电机功角初始值会自动带上这个相角差和手算的相对功角会有一个整体偏移。这个偏移不影响暂态稳定判断因为功角曲线看的是相对变化量只要故障前后的摇摆趋势正确绝对初值有偏移没关系。但如果初值偏差真的很大比如差了几十度那就一定是机械功率和电气功率不匹配导致的。检查方法很简单把故障模块禁用或者把Transition time设成一个远超仿真结束时间的值只跑稳态看发电机有功功率实测值是否等于你设置的P_m。我喜欢加一个Pm和Pe的实时比较在稳态下这两条线应该完全重合。如果不重合回到Powergui的Load Flow重新初始化或者检查发电机模块的初值参数。5.3 波形失真测量模块的采样和矢量模式问题有段时间我的三相电压波形在故障切除后出现高频毛刺看起来像振荡实际上不是系统失稳而是测量模块的处理方式问题。Simscape电气库的Three-Phase V-I Measurement有一个Mode选项可以选择Measurements输出瞬时值或者Phasor输出相量。如果你误选了Phasor模式在故障这种大扰动场景下相量输出会经过内部滤波计算产生不真实的相位延迟和幅值波动看起来就像振荡。做暂态稳定仿真一定要用Measurements模式。另一个常见问题是Scope的显示采样率不够。故障瞬间高频分量丰富如果Scope的Sample time没有显式设置Simulink可能会自动以一个比较粗糙的间隔采样导致你看到的波形丢掉了故障冲击的峰值。把Scope的采样时间设为1e-5秒10微秒或者和求解器最大步长一致就能看到完整的暂态波形。5.4 常见问题速查表现象可能原因解决措施仿真一开始就发散机械功率和电磁功率不匹配Powergui负载潮流初始化确认PmPe故障瞬间数值溢出故障电阻为0故障电阻设为0.001~0.01 Ω仿真极慢使用ode45非刚性求解器改用ode23tb最大步长0.001s功角持续增长故障切除时间超过CCT缩短故障时序检查切除时刻电压波形有伪振荡测量模块误选Phasor模式改为Measurements瞬时值模式功角初值不符合手算值基准相角偏移关注相对变化量确认稳态Pe正确变压器初始化失败饱和参数未正确设置关闭饱和选项先跑通基础模型6. 一些实操心得和扩展方向最后再分享一些我在反复搭建调试这个模型后的个人体会。最核心的一条心得是搭建电力系统暂态稳定模型不要指望一步到位。我的习惯是先不接故障模块只用无穷大母线、变压器、线路、发电机搭一个纯稳态模型跑通并确认所有波形正常。然后再把Three-Phase Fault模块加进去先用很短的故障时间比如1个周波测试故障逻辑确认故障时序正确后再逐渐加长故障持续时间、研究CCT。这样每加一个环节你都能精确控制问题的来源而不是五个模块一起出错的时候手忙脚乱。第二个建议是充分利用MATLAB脚本和Simulink的交互。在求CCT时我写了这样一个循环逻辑用一个数组存放不同故障切除时刻循环修改模型的参数并运行仿真然后自动检测功角信号的最大值是否超过阈值比如150°来判断稳定与否。这样做的好处是可以用一个for循环批量扫参数不用每次手动改Transition time重跑效率高很多。这个思路同样适用于扫描不同机械功率水平或者不同系统电抗下的稳定裕度能很快画出一条稳定域边界曲线。第三个经验是关于参数体系的统一。做电力系统仿真一定要始终使用标幺值有一个方便验证的方法是把手算的稳态功率、功角值和Simulink的潮流计算结果对比如果偏差超过1%一定是你某个地方的基准值搞错了。这个习惯能帮你避免在错误模型上花费大量调试时间。这个模型的扩展方向也很多。把Simplified Synchronous Machine换成Standard模型加上一个简单的AVR励磁调节器就能研究励磁系统对暂态稳定的改善效果你会发现加了快速励磁后最大摇摆角明显减小。在两个无穷大母线之间接两组单机系统就变成了双机系统可以研究两个发电机转子之间的相对摇摆。如果再把调速器模型加上还能分析机械功率调节对稳定性的影响。可以说单机无穷大系统这个基础模型所有环节都吃透了后面无论做多复杂的仿真核心的分析框架都还是这套运行点初始化、扰动设置、功角观察、稳定判断。把这个基本功打牢比堆砌复杂的模型要有价值得多。