接手这个题目我也比较有感触因为很多学习者刚拿到39节点相关模型的时候第一反应往往是双击打开然后对着满屏的三相线路和一堆彩色信号线发呆甚至直接卡在“模型初始化不收敛”这一步就放弃了。其实10机39节点系统并没有那么神秘它就是一个中等规模、结构清晰的交流电网测试算例。这篇内容我会尽可能把模型背后的系统结构、搭建路径、稳态初始化、暂态故障仿真以及我实际踩过的坑一次讲清楚目标是让一个能看懂单机无穷大系统仿真的读者也能比较从容地跑起来多机系统仿真。1. 拨开39节点的外衣这套模型到底在模拟什么电网1.1 经典“新英格兰”模型的来历与用途10机39节点系统通常也被称为新英格兰系统New England 39-bus System它最早是上世纪中叶为了研究大规模电力系统动态行为而简化的一个区域电网模型并不是凭空设计的理论算例而是从实际电网中提取了关键断面、主要电源点、负荷中心之后得来的。正因为如此它比纯理论的测试系统多了一层“真实感”。后来这套模型被广泛用于电力系统暂态稳定分析、低频振荡研究、继电保护动作逻辑验证、广域测量与控制等方向。你在IEEE期刊或学位论文里看到的“39节点系统”大量算例大部分本质都是同一套公开数据只是搭建平台不同有的用BPA、有的用PSS/E、有的用Matlab/Simulink还有的用Python调用或自编算法。这一点非常重要数据同源结论才可对比。对我来说这套系统的最大意义是提供了一个“不算太大但足够复杂”的电网样本。单台发电机带无穷大系统只能看单机行为4机2区域又太简单而39节点系统的机群之间、区域之间确实存在可观察的摇摆和振荡现象跑到暂态稳定阶段时曲线里能看出很多名堂。1.2 一张拓扑认知地图39个节点、10台机、46条支路动手仿真之前必须先搞清楚系统骨架。39节点系统共有39条母线其中10条是发电机母线通常编号为30到39系统基准容量是100 MVA输电网络电压等级约为345 kV。这10台发电机并不都挂在同一电压等级下它们通过升压变压器接入高压输电网发电机母线之外的1到29号母线则是负荷节点和中间联络节点。支路数量方面常见的公开数据版本里有46条支路其中大部分是交流输电线路还有一组双绕组变压器支路。系统的总负荷量级大约在6000 MW左右不同数据版本会有一点差异但整体规模就是这个级别。传统上这套系统还经常被划分为三个区域区域之间通过几条联络线相连用来研究区域间低频振荡。发电机在其中的角色也不一样31号母线发电机通常作为参考机Slack其余机组按给定有功出力运行有的参与电压调整有的设置成恒定无功或电压控制。理解这个区分很重要因为后面在Simulink里做潮流初始化时如果参考机设错了整个初始化很可能失败。记住一点在网络拓扑中所有母线的连接关系、线路电抗、变压器变比、负荷大小都是公开标定好的你不需要“发明”这些数据你需要做的是把它们正确映射到所选仿真平台的参数表里。1.3 为什么研究多机电力系统几乎都绕不开这一档规模单机无穷大系统主要用来理解基本原理比如功角稳定性的含义、等面积法则。但实际电网的问题从来不是单机行为而是多机之间的相互影响。39节点系统的节点规模刚好能支撑“不同机组群之间有相对摇摆”的场景故障发生在一个区域远处机组的功角摆动和近处机组的功角摆动幅度不一样你会看到有些机组在首摆阶段就冲得很远有些机组反而向反方向摆。这就是研究多机电力系统最核心的兴趣点之一失步模式。相比之下节点数太少体现不出这种模式差异节点数太多比如上百节点又会让初次接触者陷入参数海洋。39节点算是一个非常好的中间档位既能跑出差别的动态行为又不至于让普通笔记本电脑的仿真时间难以接受。另外值得一提的是很多经典控制策略的验证都在这套系统上做过比如电力系统稳定器PSS参数设计、切机切负荷策略、广域阻尼控制等。你后续想把研究往哪个方向深入39节点都能作为起步平台。2. 两条搭建路线怎么选自己拖元件还是先用数据文件打底2.1 路线A从Simscape Electrical库手搓全模型第一条路线是完全按物理结构在Simulink里搭建。打开Simulink左侧库浏览器里找到Simscape Electrical再找Specialized Power Systems模块集这里有三相电压源、三相变压器、三相π型线路、三相负荷、同步电机、断路器、故障模块等。早期版本的模块名可能叫SimPowerSystems功能基本一致。手搓模型的核心动作就是按拓扑结构把39条母线对应的电气节点用线路模块连接起来把10台同步电机接到对应升压变压器的高压侧或指定母线把负荷模块挂在负荷节点上。纯手工搭建的工程量不小尤其容易在接线阶段出现“漏了一条支路”或者“两条母线电气连接点看错了”的低级错误。但这条路线的优点也很明显所有元件都在眼前物理概念清晰想改哪条线路的参数、想换哪台发电机的模型都是直接双击模块改参数就行。对于学习电力系统建模的人来说这个过程本身就是在加深理解。不过纯手工接线之前建议先把母线编号和关联支路列成一张表再对照着表去连线不要凭记忆拖线。哪怕是少连一根联络线后面潮流初始化时都会表现为部分节点电压异常。2.2 路线B先用Matpower数据文件打好底子第二条路线是先借助静态分析工具箱把数据吃透。Matpower是一个高度成熟的开源潮流与最优潮流工具箱它自带一个名为case39的标准算例文件里面完整定义了39条母线的基准电压、负荷功率、发电机出力上限下限、线路阻抗、变压器变比等全部数据。在Matlab命令行里执行mpc case39; res runpf(mpc);就能得到整套系统的潮流解并打印出各母线电压幅值、相角和发电机出力。这些结果就是你后面做动态仿真初始化时需要输入的稳态工作点。这里要澄清一个容易混淆的地方Matpower本身不做电磁暂态或机电暂态仿真它给出的是静态潮流。Simulink才是做动态仿真的地方。所以Matpower的角色更像是“数据仓库”和“稳态核对工具”先用它确认正确的运行点长什么样再进入Simulink动态模型。实际上很多Simulink专家也是这么干的先在Matpower里跑通case39把母线和发电机数据导出成表格再用这份表格指导Simulink里的模块参数填写。这样能大大减少“一边建模一边迷茫”的时间。2.3 我的选型建议混合使用最稳如果只做潮流分析和最优潮流研究直接用Matpower就够了根本不需要Simulink。但如果目标是暂态稳定、故障仿真、控制策略设计那就绕不开Simulink动态模型。我的建议是先花半天时间把Matpower的case39数据完整跑一遍然后再进入Simulink手搓模型这两条路线并不是二选一而是前后接力。当然网上也能找到很多现成的Simulink版39节点模型下载之后不要急着运行先做三件事查看Powergui里选择的仿真类型查看同步电机模块是否带励磁和调速器查看负荷是否恒定PQ。这三个点决定模型能否代表实际物理系统也决定你后面的分析结论靠不靠谱。下表是我个人对两种路线的定位路线适合场景优点主要代价Simulink手搓暂态稳定、故障仿真、控制器设计物理结构直观、改参数灵活接线工作量大、排查错误耗时Matpower数据文件潮流分析、最优潮流、稳态校核数据标准、跑起来快不能做动态时域仿真3. 仿真前必须过的一道坎让39个母线先“停在”正确的稳态工作点3.1 为什么电力系统仿真不能上来就做短路不少人第一次做故障仿真时习惯性地搭好模型直接放一个短路模块点击运行结果看到发电机功率和转速波形像过山车一样乱冲甚至直接发散。其实这往往不是故障本身造成的而是因为模型根本没有处于正确的稳态初始状态。电力系统仿真和汽车启动的道理类似你得让发动机先稳定怠速一段时间才能去踩油门测试加速表现。如果发动机本身在剧烈抖动你记录下来的数据就分不清是故障响应还是初始扰动。在Simulink环境里同步电机模块内部有一套完整的电磁暂态方程它的初始状态必须与电网潮流解一致。也就是说在短路发生之前所有母线的电压幅值、相角所有发电机的有功输出、无功输出、励磁电压都要恰好落在潮流解上。这个“落在潮流解上”的过程就是初始化。3.2 Powergui的Load Flow初始化与参数填写原则如果模型里有专门处理电力系统的Powergui模块双击它就能看到Load Flow工具。不同Matlab版本界面文字稍有差异但核心操作逻辑差不多先设定系统频率通常50 Hz或60 Hz取决于你的数据版本39节点原版数据常按60 Hz处理然后在Bus页面填入各母线电压幅值与相角的初值估计再在Machine页面设置各同步电机的运行方式是作为Slack参考机还是按PV节点给定有功和机端电压或是按PQ节点给定有功无功。填写初值最稳妥的方法就是把Matpower跑出来的潮流结果直接抄进来。比如母线16的电压幅值如果潮流结果是0.98 pu角度是-5度那初始值就填接近这个的数。填好之后点击执行潮流若收敛再把潮流结果推送到同步电机的初始化状态中。这一步等于告诉Simulink所有发电机的转子位置、励磁、功率输出已经就位。实际操作时还要注意一个隐含条件模型中所有同步电机的额定容量可能不一样。case39数据里各台机的容量并不是完全统一的因此填发电机有功出力时要注意单位到底是MW还是标幺值是在系统基准100 MVA上标幺还是在发电机自身容量上标幺。这往往是初始化反复失败的隐藏原因之一。3.3 潮流初始化失败时我通常会按什么顺序排查初始化失败时的报错信息往往很笼统但我的排查顺序比较固定。先看参考机系统里必须有一个且仅有一个Slack母线参考机的有功和无功是自动平衡的不能把它也设置成定功率。再看发电机有功之和所有非参考机的有功出力加上网损要等于总负荷有功如果差太多潮流没有可行解。接着看无功很多初始化失败其实是无功功率不平衡尤其是负荷侧无功和发电机无功上限设置不当。最后看线路参数如果某条线路的电抗值填错了变成0.01和1.0的差别潮流会严重畸变或直接奇异。还有一个非常实用的习惯初始化完成后不要立刻设置故障先仿真前0.5秒或1秒只观察各母线电压和发电机电磁功率。如果波形在初始阶段就是平的说明初始化成功如果一开始就出现大幅振荡说明初始工作点没对准这时候去翻潮流解通常能发现问题。我在多台不同版本Matlab上都试过类似的流程结论是初始化这一步的耗时往往占整个建模周期的一半以上它不是一个可以跳过的“格式性操作”而是决定模型可信度的基础环节。4. 制造一次三相短路从故障设置到功角曲线判稳4.1 最直观的教学场景母线近端三相短路做到这一步说明前面的稳态工作点已经稳定住了可以开始给系统施加扰动。最典型、也最适合入门的故障类型是三相对称短路Simulink模块库里对应的是Three-Phase Fault。在仿真中这个模块可以串联或并联在指定位置最常见的接法是并联在某条母线或其出线上模拟该点发生三相接地短路。以一条靠近负荷中心的母线为例比如母线16附近发生三相短路我习惯的把故障时序设为0.1秒故障合上0.3秒故障断开。也就是说短路持续时间为0.2秒。之所以不从0秒开始故障是为了先给一段干净的稳态时段让初始化后的系统稳定运行也有利于后期对照故障前后的波形。Three-Phase Fault模块里有故障相选择这里应选三相短路故障阻抗很小比如0.001Ω。Transitional Time一般设为零表示理想瞬时切换。若想模拟保护动作后线路切除可以再放一个Three-Phase Breaker模块在0.3秒同时断开故障所在线路让系统在故障切除的同时改变拓扑。第一次做时建议先只加故障不加跳闸观察最纯粹的短路响应。4.2 解算器与仿真时长设置是另一条命脉很多初次跑通39节点模型的人会遇到仿真异常不是模型物理上错了而是解算器选得不对。Simscape Electrical的连续系统默认情况下建议使用变步长解算器我个人的起手提法是使用ode23tb相对容差设置1e-4最大步长根据故障暂态需要设置为1e-3秒量级。如果选显式ode45在故障瞬间很容易产生数值振荡或步长过小表现为仿真进度条几乎不动或波形出现锯齿状毛刺。那是因为系统在短路瞬间包含快速电磁暂态分量显式算法对这类刚性系统不太友好。ode23tb和ode15s这类隐式算法则更适合处理此类刚性问题。仿真时长方面观察暂态稳定问题通常至少要跑20秒以上。因为故障清除后的暂态过程不是马上结束机组之间的功率摇摆会持续好几秒甚至十几秒。只跑5秒常常看不到完整的首摆、反摆和后续衰减过程。如果你只关心故障瞬间电压跌落那短时间仿真也可以但若目标是功角稳定判断时间窗口必须给足。4.3 判断稳定与失稳不能只看电压波形故障发生后母线电压跌落是必然现象电压恢复曲线主要反映暂态电压特性。但电力系统暂态稳定最核心的判断是发电机转子之间是否保持同步也就是要看功角或转速相对变化。Simulink里同步电机模块通常提供测量信号包含转速、转子角、电磁功率等通道。通过这些通道可以看到各台发电机在故障冲击下的转速波动。稳定情况下故障清除后各机组转速会在同步转速附近摆动振幅逐渐衰减并趋于一致失稳情况下会出现某台或某几台机组转速持续偏离相对其他机组的角度越拉越大最终失去同步。我通常会把多台发电机的转速曲线画在同一张图里观察是否出现明显的“分群”趋势。若一部分机组摆向一侧、另一部分摆向另一侧那就是典型的机群间失步模式。临界切除时间CCT也是一个很值得做的延伸实验。它的含义是在故障持续时间的临界值以内清除故障系统还能保持稳定超过这个时间再切除系统就失稳了。操作上就是用0.15秒、0.2秒、0.25秒……逐渐增加故障持续时间找到那个“稳定”变成“失稳”的分界点。这个指标对理解系统稳定裕度非常有帮助。5. 结果不是看个热闹电压、功角、频率曲线怎么读才算数5.1 信号抓取和保存不要只在Scope里看个亮刚开始做仿真时我直接把Scope拖出来看波形觉得挺直观。但问题马上来了Scope里的波形放大不方便数据也不好做二次处理。后来我养成一个习惯关键信号都额外用To Workspace模块送进Matlab工作空间或者使用模型数据记录器把仿真数据存成结构体或数组再用plot命令自己画图。以母线电压为例可以在测量模块之后接一个求有效值或幅值的环节把故障前后的电压幅值变化存下来。发电机转速、电磁功率、转子角等信号同样处理。这样后期做论文插图、对比不同参数方案时不需要反复运行仿真直接从工作空间里调数据即可。5.2 母线电压曲线里藏着哪几层信息母线电压曲线是“最外行也能看到变化、内行能读出信息”的波形。我一般按三个层次来读。第一层是故障瞬间跌得多深近端三相短路时电压可能跌到0.2 pu甚至更低这反映了系统电气距离和短路电流水平第二层是故障切除后的恢复过程是立刻弹回来、小幅振荡几次后恢复还是一路下滑第三层是恢复后是否真正回到稳态电压持续偏低便可能涉及电压稳定问题。下面这个表是我常用的快速判读思路电压曲线现象可能的系统状态故障瞬间跌至0.2 pu以下近端短路、电气距离近切除后小幅振荡并逐渐平稳系统暂态稳定电压持续振荡衰减很慢系统阻尼偏弱可考虑PSS切除后电压爬不起来甚至继续下降可能趋向电压失稳5.3 从波形到论文级结论频率与有功响应的分析框架母线电压之外频率响应和有功功率响应也是解读重点。当系统受到扰动后发电机转速先快速调整进而影响系统频率随后调速器动作改变有功出力。频率波动的变化率通常反映系统惯性大小而波动的衰减速度反映阻尼强弱。如果你在波形里看到振荡衰减慢持续好几个周期都消不掉这就是一种阻尼偏弱的迹象。此时可以尝试在励磁系统旁边加电力系统稳定器PSS通过引入附加励磁信号来增加阻尼。把加PSS前后的波形放一起对比振荡衰减速度差异会非常直观。把仿真结果写成结论时我通常按照这样组织先写明工况故障位置、故障类型、持续时间、切除方式再给关键波形发电机功角、母线电压、系统频率然后对照稳定判据或性能指标做说明最后得出结论。比如“0.2秒切除可保持暂态稳定0.25秒切除失步该工况下临界切除时间约为0.22秒”这种写法在工程报告里很有用。6. 反复出现的仿真翻车现场与我最后留下的一版参数6.1 发电机模型“看起来一样”结果却天差地别多机系统仿真里最容易翻车的一个坎是发电机模型阶数不一致。同一个39节点数据用二阶经典模型和用六阶次暂态模型短路故障后的首摆可能差不多但后续的振荡频率、幅值、衰减特性都可能明显不同。这不是说谁对谁错而是你必须明确说明自己用的是哪一类模型。尤其在做论文时如果只写“采用同步电机模型”而不说明阶数和励磁系统类型结果基本不可复现。我的建议是核心算例用六阶模型加励磁系统至少要包含IEEE I型励磁或简化励磁快速试算阶段可以用二阶模型但要在报告中注明。6.2 基准值、单位、接线错误三座大山这三个坑我几乎每次带人做39节点仿真都会遇到。基准值的坑在于case39数据的很多电抗标幺值是基于系统基准100 MVA的而同步电机的额定容量可能是几十MVA到几百MVA不等。你把case39数据里的d轴电抗标幺值直接填进同步电机模块时模块内部默认的基准可能不是100 MVA若不换算短路电流和暂态响应都会失真。单位的坑更隐蔽Simulink的线路模块可选国际单位或标幺值单位比如π型线路模块用欧姆、亨利、法拉每公里而case39数据里线路阻抗常用标幺值。换算的时候基准电压345 kV、基准容量100 MVA阻抗基准大约是1190欧数值上接近把标幺电抗乘上这个基准得到欧姆值再除以线路长度得到每公里值。这个换算过程只要错一位数潮流就完全偏离。接线错误则是手搓模型的通病某条联络线少画一段、变压器极性接反、负荷节点挂错母线这些都不太容易一眼发现但会导致个别母线电压异常或潮流不收敛。我后来用过一个很土但有效的方法先用彩色笔在纸质拓扑图上把已连接支路逐个划掉确保所有支路都落了地再去谈仿真。6.3 仿真太慢和波形毛刺的处理方法39节点系统虽然规模不算大但20秒时域仿真在普通笔记本上可能需要跑很久特别是在故障瞬间变步长解算器会把步长压得很小。如果仿真时间实在难以接受我一般按这个顺序做优化先把发电机模型从六阶降为二阶试跑通再把相对容差从1e-4放宽到1e-3接着限制最大步长让它在非故障时段稍微大一点。三步做完仿真耗时往往能下降一个量级。波形毛刺的问题要分两类看一类是数值毛刺通常表现为极高频的锯齿波动这大概率是步长过大或解算器不适合减小最大步长能缓解另一类是真实的电磁暂态高频分量故障瞬间出现短暂高频波动是正常的不必刻意滤掉反而应该保留以便观察暂态细节。判断标准是看毛刺是否在故障后几十毫秒内消失。6.4 最后留下的一版“稳当起手参数”经过反复折腾我自己最终沉淀了一套比较稳妥的起手配置不敢说是最优但很适合作为第一次跑39节点的底稿。直接照这个配置跑通以后再根据你的具体问题调整才有意义。配置项我使用的参数解算器ode23tb变步长相对容差1e-4最大步长1e-3仿真时长20 s故障设置母线16三相短路0.1 s合0.3 s断故障阻抗0.001 Ω发电机模型六阶同步电机模型带简化励磁参考机设为31号母线机组负荷模型恒功率PQ负荷按case39数据设置信号记录各发电机转速/转子角关键母线电压幅值这套配置不是金标准但它跑通的概率很高。等你把波形跑出来以后再针对具体问题调整比如把负荷换成恒阻抗模型、给励磁系统加PSS、延长或缩短故障持续时间这时你已经有了一条可以参照的基线。最开始我在Simulink里搭39节点系统时最耗时间的其实不是画图而是初始化反复不收敛。后来我学乖了每做一套新模型都会先借助工具把潮流解跑出来再往Simulink里填数据整个过程顺畅很多。建议你也把“先稳态后暂态”这件事当成铁律而不是一上来就想着加短路故障看热闹。后面如果想向IEEE 118节点这种更大规模延伸思路是一样的只是数据准备和模型验证的工程量更大。
