简介本资源是一套基于MATPOWER工具箱的电力系统潮流计算实践代码包面向电气工程专业本科生、研究生及电力系统仿真初学者用于掌握多机多节点系统的建模、导纳矩阵构建与稳态潮流分析全流程。压缩包共6个文件含3个核心MATLAB脚本用于构建节点导纳矩阵、数据校验及主案例定义、1个MATPOWER 7.1安装压缩包、1份详细使用说明文档.docx和1个文本说明文件.txt总大小32.09MB。已有4151人学习下载反映出其在教学实验与课程设计中的广泛适用性。用户可直接运行Case28_08.m调用runpf完成8机28节点系统的潮流计算代码中完整实现了发电机出力约束设置、含变压器支路建模、节点无功补偿配置及三类基础数据矩阵节点、发电机、支路的规范构造具备良好的可读性与复现性是理解MATPOWER案例文件结构与电力系统建模逻辑的实用参考。 Matpower自带案例库里有9节点、14节点、30节点、39节点这些熟面孔但提到“8机28节点”老手第一反应通常是这不是标准型号。这意味着你得自己搭一个可运行的电力系统模型。这不算难但里面琐碎的坑非常多尤其是数据结构不熟的时候一个列填错就能让潮流迭代半天还不收敛。我最近在一个项目里正好用Matpower完整搭建了一套8机28节点系统从母线数据到发电机组、支路参数、最优潮流再到结果可视化整个过程记录下来给要做论文仿真、毕业设计或者工程测试算例的朋友做个参考。1. 为什么是8机28节点自建测试系统的现实需求1.1 Matpower到底是什么Matpower是一个基于MATLAB的开源电力系统分析工具箱核心功能包括潮流计算power flow、最优潮流OPF、连续潮流CPF、状态估计等。它最大的优势是数据接口极其简洁把整个电力系统定义成一个结构体mpc包含母线、发电机、支路、成本函数四张核心矩阵然后调用一句runpf(mpc)就能完成交流潮流计算runopf(mpc)完成最优潮流。对研究者来说不需要从零写牛顿-拉夫逊迭代也不需要自己处理稀疏矩阵的编号问题省下来的时间可以全部花在算法和策略上。1.2 28节点为什么值得“手搓”IEEE标准测试系统里确实没有28节点实际项目中遇到8机28节点一般有三种来源一是从30节点系统删改而来二是来自某篇论文的算例数据三是仿真者根据实际电网抽象出来的简化模型。无论哪种来源核心都指向一个需求需要一个规模适中、机组数量足够多、能体现多机协调特性的测试系统。8台发电机意味着系统可以在潮流层面讨论机组出力分配、经济调度在动态层面可以研究机组间的低频振荡、暂态稳定甚至故障后的功率转移。28个节点又保证了拓扑结构有一定的复杂度——可以设置不同的电压等级、多条联络线、多个负荷中心而不像9节点或14节点那样“一眼看穿”。1.3 这类自建系统最常用的场景我的实际体验是这类系统最集中的需求集中在三个方向算法验证比如粒子群优化求解经济调度、神经网络做负荷预测后的潮流校验、强化学习做电压控制都需要一个“比小系统有说服力又比大系统跑得快”的算例。教学实验本科生或者研究生做课程设计需要一个能手动改写参数、观察不同运行方式影响的中等规模系统。工程预研在真实电网数据不方便公开时用一个结构类似的抽象系统做方案验证。所以8机28节点不是凭空想出来的数字它恰好是“教学与研究之间”的一个平衡点。2. bus、gen、branch三张表Matpower数据结构的底层逻辑所有Matpower程序设计都绕不开三张表mpc.bus、mpc.gen、mpc.branch。这三张矩阵就是整个电力系统的“骨架”搞不清楚它们的列含义后面所有工作都会像在迷宫里打转。2.1 全局参数和数据结构先看一个典型的Matpower数据文件开头function mpc case28_8m mpc.version 2; mpc.baseMVA 100;baseMVA是系统功率基准值通常取100 MVA。潮流计算中电压是标幺值功率是有名值MW/Mvar支路阻抗是标幺值。很多人第一次接触时容易把Pd、Qd写成标幺值这是非常常见的错误——Matpower要求母线负荷直接填有名值MW/Mvar阻抗才是标幺值。2.2 bus矩阵电力系统的“节点档案”mpc.bus一共13列每一行代表一个节点。列含义如下列号含义说明1bus_i节点编号必须正整数且互不重复2type节点类型1-PQ节点2-PV节点3-平衡节点4-隔离节点3Pd有功负荷MW4Qd无功负荷Mvar5Gs并联电导MW1pu6Bs并联电纳Mvar1pu7area分区编号默认18Vm电压幅值初始值pu9Va电压相角初始值度10baseKV基准电压kV11zone损耗区编号默认112Vmax电压上限pu13Vmin电压下限pu写文件时最容易犯的错误是type选错。功率平衡节点只能有一个通常选在装机容量最大的母线上负责平衡全网功率差额PV节点接发电机给定有功出力和电压幅值PQ节点是普通负荷节点。如果平衡节点缺失潮流计算会直接报“slack bus not found”之类的错误。2.3 gen矩阵发电机的“运行档案”mpc.gen矩阵每一行对应一台发电机核心前10列分别是母线号、有功出力Pg、无功出力Qg、无功上限Qmax、无功下限Qmin、电压设定值Vg、机组容量基准mBase、运行状态、有功上限Pmax、有功下限Pmin。这里有一个非常关键的细节mpc.bus里的Vm和mpc.gen里的Vg都要给一个合理的电压初值。PV节点的电压设定值在bus矩阵第8列和gen矩阵第6列都有体现我个人的习惯是两处保持一致避免排查问题时被自己搞晕。实际项目里很多人只改了bus表里的电压忘了gen表里的Vg最后结果和预期对不上。2.4 branch矩阵电力系统的“连接关系”mpc.branch每一行代表一条输电线路或者变压器支路前13列的含义是列号含义说明1fbus首端节点编号2tbus末端节点编号3r电阻pu4x电抗pu5b充电电纳pu6rateA长期载流量限制MVA7rateB短期载流量限制MVA8rateC紧急载流量限制MVA9ratio变压器变比10angle变压器移相角度11status投运状态1-投运0-停运12angmin相角下限度13angmax相角上限度纯线路的ratio填0有载调压变压器则填变比值。计算支路负载率时Matpower调用的就是rateA这一列如果全填0运行结果里支路潮流依然能算出来但负载率显示为0这一细节在做安全分析时特别容易踩。2.5 gencost最优潮流的经济参数mpc.gencost是成本函数矩阵runopf求解最优潮流的依据就在它里面。一般用多项式模型指定机组有功出力的二次成本函数形如c2*P^2 c1*P c0。由于不同版本列数定义略有差异我建议直接用loadcase(case30)打开自带案例仿照它的gencost格式去改永远比对着文档手搓要稳。3. 从零到可运行8机28节点系统的搭建全流程在搞清楚三张表之后剩下的工作就是“填格子”。但怎么填、填多少、怎么验证是有讲究的。3.1 最稳的起点从case30改造我强烈不建议从空白文件开始盲写28节点。实际项目里最快的路径是以case30为模板改造成8机28节点。理由很简单case30的参数是经过验证的潮流能收敛负荷水平、阻抗数量级、电压范围都是合理的。你只需要做三件事删掉两个节点、增加两台发电机、重新整理编号。在MATLAB里手动改case30.m比较麻烦更推荐用代码的方式处理% 加载自带30节点系统 mpc loadcase(case30); % 删除第29号和第30号节点以及连接到它们的支路 bus_del [29; 30]; branch_del []; for i 1:size(mpc.branch, 1) if ismember(mpc.branch(i, 1), bus_del) || ismember(mpc.branch(i, 2), bus_del) branch_del [branch_del; i]; end end mpc.branch(branch_del, :) []; mpc.bus(bus_del, :) []; % 删除与这两个节点关联的发电机组如果有 gen_del []; for i 1:size(mpc.gen, 1) if ismember(mpc.gen(i, 1), bus_del) gen_del [gen_del; i]; end end mpc.gen(gen_del, :) []; % 修改节点编号删除中间编号后会产生跳号 % 使用 ext2int 重新编号 mpc ext2int(mpc);ext2int这个函数是Matpower内部的编号统一工具可以把不连续的外部节点编号映射成1到N的连续内部编号。很多人在删除节点后不重新编号结果bus矩阵第1列出现了29、30这种空洞虽然Matpower多数情况下能处理但后续可视化、结果导出到Excel时非常容易出错。3.2 增加两台发电机删除节点后原case30系统是6台发电机离8机还差两台。加机组的时候要考虑节点类型的变化% 在第5号和第12号母线上分别加一台发电机 newgen1 [5, 40, 0, 60, -60, 1.0, 100, 1, 120, 0, ... 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]; newgen2 [12, 50, 0, 80, -40, 1.0, 100, 1, 150, 0, ... 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]; mpc.gen [mpc.gen; newgen1; newgen2];注意mpc.gen一个矩阵有21列版本不同列数可能略有差异补齐时必须把后面不用的列都填0。同时要把mpc.bus中对应节点的type从1改为2表示它变成了PV节点mpc.bus(5, 2) 2; mpc.bus(12, 2) 2;这一步很多人会忘结果新增的发电机在潮流计算里根本没被当成PV节点系统里的机组数量看着是8台实际参与电压调节的只有原来的6台。3.3 从零开始写28节点数据文件如果手里已经有文献里的系统参数不愿意从case30改也可以直接写数据文件。我的建议是先把单线图拓扑画清楚标好节点号、支路阻抗、负荷大小、发电机接入点然后按顺序填表。写之前先规划好3个问题系统的基准电压用什么等级全部用同一个电压等级最省事比如230kV否则变压器支路的变比参数要仔细算。平衡节点放在哪选一台容量大、有调节能力的机组所在节点通常放节点1。系统总负荷量级是多少自建系统建议控制在300~500MW之间太大会导致机组出力超过Pmax太小则电压普遍偏高。实际编写时数据文件的格式就是把我前面讲的三张表依次填好保存成case28_8m.m即可。写完之后用loadcase加载再调用check_mpc做基础校验这个函数会检查矩阵尺寸、节点编号、支路端点是否越界等常见问题。不同版本check_mpc调用方式略有差异但载入数据时如果数据有致命错误命令窗口通常会直接报错或者给出warning。3.4 跑通第一个潮流计算数据文件准备好了运行潮流只需要几行代码mpc loadcase(case28_8m); opt mpoption(verbose, 2); res runpf(mpc, opt);verbose设置为2可以在命令窗口看到牛顿法每次迭代的最大功率失配量。正常情况下3~5次迭代就收敛了res.success返回1。如果verbose0则只显示收敛结果不显示迭代过程。跑通之后先别急着高兴先做一次基础自检查看发电机总出力、负荷总量、网损、各节点电压幅值确认数据没有“表面收敛但结果离谱”的情况。4. 潮流结果怎么读从res结构体到可视化runpf返回的结果是一个结构体res里面包含了完整的潮流计算结果。很多新手跑完只看一个res.success就结束了这太浪费了。4.1 结果结构体的关键字段res里的核心字段和mpc是同构的也包含bus、gen、branch三张表但各表内容已经被刷新成了潮流计算后的运行状态res.bus(:, 8)各节点电压幅值pures.bus(:, 9)各节点电压相角度res.gen(:, 2)各发电机实际有功出力MWres.gen(:, 3)各发电机实际无功出力Mvarres.branch(:, 14)和res.branch(:, 15)支路首末端有功潮流MWres.branch(:, 16)和res.branch(:, 17)支路首末端无功潮流Mvar注意不同版本Matpower里branch矩阵在结果中的实际潮流列位置可能有差异。最稳妥的方式是在命令窗口输入res.branch(1, :)看一下前几列确认哪些列是有功、无功再看结果。4.2 怎么判断结果是否合理一个可用的潮流结果必须满足几个基本条件电压幅值基本在0.95~1.05 pu之间如果出现0.8或1.2这种数值说明系统无功分布有问题。相角是逐步变化的相邻节点相角差通常在5度以内出现几十度的相角差说明支路参数异常。发电机无功出力在Qmin和Qmax范围内。如果某台机组无功出力打到上限说明该节点附近无功不足系统可能处于临界运行状态。总有功出力减总负荷等于网损网损一般在总负荷的2%~5%之间。如果网损达到负荷的30%以上支路阻抗的数据几乎可以肯定填错了。判断网损在Matpower里可以直接算P_gen_total sum(res.gen(:, 2)); P_load_total sum(res.bus(:, 3)); P_loss P_gen_total - P_load_total; fprintf(系统总负荷: %.2f MW\n, P_load_total); fprintf(发电机总出力: %.2f MW\n, P_gen_total); fprintf(网损: %.2f MW (%.2f%%)\n, P_loss, P_loss / P_load_total * 100);4.3 可视化与报告导出Matpower本身不自带特别花哨的画图工具但用MATLAB基础绘图功能就足够完成大部分展示需求。我常用的三个可视化手段节点电压幅值柱状图快速判断全网电压水平是否正常。figure; bar(res.bus(:, 1), res.bus(:, 8)); xlabel(节点编号); ylabel(电压幅值(pu)); title(28节点系统各母线电压幅值); ylim([0.9 1.1]);支路负载率排序图找到系统的重载线路。loading res.branch(:, 14) ./ res.branch(:, 6); figure; bar(loading * 100); xlabel(支路编号); ylabel(负载率(%)); title(各支路负载率);注意如果rateA填了0这里的负载率全是Inf或NaN所以我前面才强调rateA一定要填。相角分布曲线能看到以平衡节点为参考的相角变化趋势。在做课程报告或者论文的时候我通常还会把res.bus和res.branch导出为Excel或者CSV。writetable(array2table(res.bus), result_bus.xlsx); writetable(array2table(res.branch), result_branch.xlsx);这样后续在Word或者LaTeX里引用数据很方便不用每次重新跑一遍再手动抄写。5. 排错实战搭建与仿真中最容易踩的坑这部分是我最想写的因为8机28节点系统本身不算复杂真正的门槛全在那些“看似不起眼的错误”上。5.1 潮流不收敛先查数据再调迭代好几次项目里我搭建的自定义系统第一次跑runpf都会碰到不收敛。Matpower的提示一般是“Newtons method power flow did not converge in 10 iterations”。碰到这个提示很多人的第一反应是增加迭代次数这完全是本末倒置。我的排查顺序是固定的检查bus表里有没有type4的隔离节点如果负荷节点被错误标为隔离节点系统会出现功率缺口。检查每个节点是否至少连了一条支路孤立节点会让雅可比矩阵奇异。检查支路x是否为0或者极小值如果某条支路电抗填了0.0001相当于一条短路支路矩阵条件数急剧恶化。检查负荷总量和发电机总出力是否匹配如果负荷远大于发电机出力上限潮流很难收敛到可行域。检查是否有节点电压初值设置得过低比如0.8pu导致迭代起点太差。有一次我排查了很久最后发现是某条支路我复制粘贴时把x0.06写成了x0.006数量级差了一个数量级潮流就怎么都算不拢。所以数据文件的校验一定要做。5.2 PV节点无功越限一个最容易忽略的细节Matpower的牛顿法在迭代过程中会考虑发电机的无功上下限。当某台发电机需要发出的无功超过Qmax时它会自动把该节点从PQ节点处理给定无功在Qmax处电压由潮流自然决定。这个机制本身是合理的但它带来的坑是你的8机系统里某台发电机可能名义上是PV节点实际上因为Qmax设置太小在结果里已经变成了无功封顶的状态电压远远低于初始设定值。很多人做最优潮流或者静态安全分析时看不到这个现象是因为默认情况下Matpower输出里只在“Generator Reactive Power Limits”部分显示Limits信息不仔细看很容易忽略。我的经验是新增发电机时Qmax/Qmin不要随意填要参考原有case30里的数值量级。比如100MVA基准下Qmax给60~80Mvar、Qmin给-60~-40Mvar是比较常见的。如果填了±9999这种极端值潮流是很容收敛但结果没有任何工程意义。5.3 节点编号与行号的混淆bus矩阵的第1列是节点编号而矩阵的行号是顺序索引。删除节点后如果不重新编号节点编号可能变成[1,2,3,5,6,7,...]这种不连续序列当你在后续程序里用res.bus(find(res.bus(:,1)5), 8)这种方式取数据时容易漏项。更麻烦的是gen矩阵和branch矩阵里引用节点编号时必须和bus表第1列保持一致。如果节点编号不连续只要有任何一个地方写错就会引用到错误的节点。所以我强烈建议在删改节点后立刻调用ext2int重新编号或者在bus表里把编号重新写成连续的1到N。5.4 支路阻抗的单位和基准问题branch表里的r、x、b都是标幺值不是欧姆也不是西门子。很多从实际工程数据换算过来的人容易把线路的欧姆值直接填进表里。标幺值和有名值的换算公式是Zpu Z_ohm * baseMVA / (baseKV^2)一个230kV、100MVA基准下的系统如果一条线路阻抗是10欧姆对应的标幺值是10 * 100 / (230^2) ≈ 0.0189 pu。如果直接把10填进去潮流结果会显示电压崩溃因为系统的电抗被放大了500多倍。5.5 建议用自带工具做一次完整检查Matpower有check_mpc这样的数据校验函数虽然不同版本调用方法不一致但它能检查矩阵尺寸、支路端点引用、重复节点编号等明显问题。我自己一般在写完数据文件后会先跑一下这个函数确认没有结构性问题后再跑潮流。另外一个笨办法是把系统节点数、支路数、发电机数打印出来和设计图对照一遍fprintf(节点数: %d\n, size(mpc.bus, 1)); fprintf(支路数: %d\n, size(mpc.branch, 1)); fprintf(发电机数: %d\n, size(mpc.gen, 1));如果数量都对不上那大概率是自己写数据文件时少了行或者多了行先在这里排查比在潮流结果里排查效率高得多。6. 从潮流到更多玩法OPF与多机系统扩展8机28节点系统的价值不止于潮流计算。把runpf换成runopf就能做经济调度意义上的最优潮流。6.1 最优潮流的快速上手调用方式非常简单res_opf runopf(mpc);但要做好OPFmpc.gencost必须填写完整否则Matpower会按默认成本计算结果可能与你预期的不符。gencost里定义的是发电机有功出力的成本曲线单位一般换算成$/h。用二次多项式表示% 三台发电机为例c2, c1, c0 分别为二次项、一次项、常数项 mpc.gencost [ 2 0 0 3 0.005 20 50; 2 0 0 3 0.010 15 40; 2 0 0 3 0.002 25 60; ];这里的MODEL2表示多项式成本模型NCOST3表示有三个系数从最高次数往常数项排列。不同版本对gencost列数的定义有细微差别我的习惯是打开自带case30.m把它的gencost直接复制过来改成自己的机组参数这样永远不会因为列数不一致报错。6.2 静态安全分析和N-1校验当系统搭好并跑通潮流后我通常还会做一轮N-1静态安全校验遍历每一条支路让它退出运行重新跑一次潮流看是否出现电压越限或支路过载。实现方式很简单n_branch size(mpc.branch, 1); for i 1:n_branch mpc_i mpc; mpc_i.branch(i, 11) 0; % 把第i条支路状态设为停运 res_i runpf(mpc_i); if res_i.success 0 fprintf(支路%d退出运行后潮流不收敛\n, i); end end这个循环虽然朴素但在8机28节点这种小系统上跑得很快是做“系统薄弱环节分析”时非常实用的工具。6.3 多机系统动态仿真的延伸如果下一步要做暂态稳定或者小干扰稳定分析8机28节点这套数据同样可以作为基础。常见做法是把构建好的mpc对象导出成PSS/E的raw文件或者直接导入到PSAT、PST等动态仿真工具里。因为在潮流层面已经验证过系统是收敛、可运行的导入到动态仿真时只需要关心动态参数惯性时间常数、励磁系统模型、调速器参数的填写不需要再回头调潮流初始值。6.4 代码包的组织习惯最后说说代码组织。我习惯把“数据文件”和“计算脚本”分开目录结构大致是这样的case28_8m.m % 系统数据文件只定义mpc结构 main_powerflow.m % 潮流计算主程序 main_opf.m % 最优潮流主程序 static_security.m % N-1校验脚本 plot_results.m % 可视化脚本数据文件里只有mpc的定义不放任何计算逻辑。主程序通过loadcase加载数据再调用Matpower的函数。这样做的好处是换一个系统数据文件主程序一行都不用改直接切换loadcase的参数就行。在实际使用中我最深的体会是Matpower项目里最花时间的永远不是runpf这一行命令而是把数据文件整理清楚。8机28节点这种自建系统关键就是先把bus、gen、branch这三张表的结构吃透后面的程序调用都是锦上添花。希望这篇分享能帮你少走点弯路。本文还有配套的精品资源点击获取
Matpower搭建8机28节点系统：从数据结构到潮流计算实战
简介本资源是一套基于MATPOWER工具箱的电力系统潮流计算实践代码包面向电气工程专业本科生、研究生及电力系统仿真初学者用于掌握多机多节点系统的建模、导纳矩阵构建与稳态潮流分析全流程。压缩包共6个文件含3个核心MATLAB脚本用于构建节点导纳矩阵、数据校验及主案例定义、1个MATPOWER 7.1安装压缩包、1份详细使用说明文档.docx和1个文本说明文件.txt总大小32.09MB。已有4151人学习下载反映出其在教学实验与课程设计中的广泛适用性。用户可直接运行Case28_08.m调用runpf完成8机28节点系统的潮流计算代码中完整实现了发电机出力约束设置、含变压器支路建模、节点无功补偿配置及三类基础数据矩阵节点、发电机、支路的规范构造具备良好的可读性与复现性是理解MATPOWER案例文件结构与电力系统建模逻辑的实用参考。 Matpower自带案例库里有9节点、14节点、30节点、39节点这些熟面孔但提到“8机28节点”老手第一反应通常是这不是标准型号。这意味着你得自己搭一个可运行的电力系统模型。这不算难但里面琐碎的坑非常多尤其是数据结构不熟的时候一个列填错就能让潮流迭代半天还不收敛。我最近在一个项目里正好用Matpower完整搭建了一套8机28节点系统从母线数据到发电机组、支路参数、最优潮流再到结果可视化整个过程记录下来给要做论文仿真、毕业设计或者工程测试算例的朋友做个参考。1. 为什么是8机28节点自建测试系统的现实需求1.1 Matpower到底是什么Matpower是一个基于MATLAB的开源电力系统分析工具箱核心功能包括潮流计算power flow、最优潮流OPF、连续潮流CPF、状态估计等。它最大的优势是数据接口极其简洁把整个电力系统定义成一个结构体mpc包含母线、发电机、支路、成本函数四张核心矩阵然后调用一句runpf(mpc)就能完成交流潮流计算runopf(mpc)完成最优潮流。对研究者来说不需要从零写牛顿-拉夫逊迭代也不需要自己处理稀疏矩阵的编号问题省下来的时间可以全部花在算法和策略上。1.2 28节点为什么值得“手搓”IEEE标准测试系统里确实没有28节点实际项目中遇到8机28节点一般有三种来源一是从30节点系统删改而来二是来自某篇论文的算例数据三是仿真者根据实际电网抽象出来的简化模型。无论哪种来源核心都指向一个需求需要一个规模适中、机组数量足够多、能体现多机协调特性的测试系统。8台发电机意味着系统可以在潮流层面讨论机组出力分配、经济调度在动态层面可以研究机组间的低频振荡、暂态稳定甚至故障后的功率转移。28个节点又保证了拓扑结构有一定的复杂度——可以设置不同的电压等级、多条联络线、多个负荷中心而不像9节点或14节点那样“一眼看穿”。1.3 这类自建系统最常用的场景我的实际体验是这类系统最集中的需求集中在三个方向算法验证比如粒子群优化求解经济调度、神经网络做负荷预测后的潮流校验、强化学习做电压控制都需要一个“比小系统有说服力又比大系统跑得快”的算例。教学实验本科生或者研究生做课程设计需要一个能手动改写参数、观察不同运行方式影响的中等规模系统。工程预研在真实电网数据不方便公开时用一个结构类似的抽象系统做方案验证。所以8机28节点不是凭空想出来的数字它恰好是“教学与研究之间”的一个平衡点。2. bus、gen、branch三张表Matpower数据结构的底层逻辑所有Matpower程序设计都绕不开三张表mpc.bus、mpc.gen、mpc.branch。这三张矩阵就是整个电力系统的“骨架”搞不清楚它们的列含义后面所有工作都会像在迷宫里打转。2.1 全局参数和数据结构先看一个典型的Matpower数据文件开头function mpc case28_8m mpc.version 2; mpc.baseMVA 100;baseMVA是系统功率基准值通常取100 MVA。潮流计算中电压是标幺值功率是有名值MW/Mvar支路阻抗是标幺值。很多人第一次接触时容易把Pd、Qd写成标幺值这是非常常见的错误——Matpower要求母线负荷直接填有名值MW/Mvar阻抗才是标幺值。2.2 bus矩阵电力系统的“节点档案”mpc.bus一共13列每一行代表一个节点。列含义如下列号含义说明1bus_i节点编号必须正整数且互不重复2type节点类型1-PQ节点2-PV节点3-平衡节点4-隔离节点3Pd有功负荷MW4Qd无功负荷Mvar5Gs并联电导MW1pu6Bs并联电纳Mvar1pu7area分区编号默认18Vm电压幅值初始值pu9Va电压相角初始值度10baseKV基准电压kV11zone损耗区编号默认112Vmax电压上限pu13Vmin电压下限pu写文件时最容易犯的错误是type选错。功率平衡节点只能有一个通常选在装机容量最大的母线上负责平衡全网功率差额PV节点接发电机给定有功出力和电压幅值PQ节点是普通负荷节点。如果平衡节点缺失潮流计算会直接报“slack bus not found”之类的错误。2.3 gen矩阵发电机的“运行档案”mpc.gen矩阵每一行对应一台发电机核心前10列分别是母线号、有功出力Pg、无功出力Qg、无功上限Qmax、无功下限Qmin、电压设定值Vg、机组容量基准mBase、运行状态、有功上限Pmax、有功下限Pmin。这里有一个非常关键的细节mpc.bus里的Vm和mpc.gen里的Vg都要给一个合理的电压初值。PV节点的电压设定值在bus矩阵第8列和gen矩阵第6列都有体现我个人的习惯是两处保持一致避免排查问题时被自己搞晕。实际项目里很多人只改了bus表里的电压忘了gen表里的Vg最后结果和预期对不上。2.4 branch矩阵电力系统的“连接关系”mpc.branch每一行代表一条输电线路或者变压器支路前13列的含义是列号含义说明1fbus首端节点编号2tbus末端节点编号3r电阻pu4x电抗pu5b充电电纳pu6rateA长期载流量限制MVA7rateB短期载流量限制MVA8rateC紧急载流量限制MVA9ratio变压器变比10angle变压器移相角度11status投运状态1-投运0-停运12angmin相角下限度13angmax相角上限度纯线路的ratio填0有载调压变压器则填变比值。计算支路负载率时Matpower调用的就是rateA这一列如果全填0运行结果里支路潮流依然能算出来但负载率显示为0这一细节在做安全分析时特别容易踩。2.5 gencost最优潮流的经济参数mpc.gencost是成本函数矩阵runopf求解最优潮流的依据就在它里面。一般用多项式模型指定机组有功出力的二次成本函数形如c2*P^2 c1*P c0。由于不同版本列数定义略有差异我建议直接用loadcase(case30)打开自带案例仿照它的gencost格式去改永远比对着文档手搓要稳。3. 从零到可运行8机28节点系统的搭建全流程在搞清楚三张表之后剩下的工作就是“填格子”。但怎么填、填多少、怎么验证是有讲究的。3.1 最稳的起点从case30改造我强烈不建议从空白文件开始盲写28节点。实际项目里最快的路径是以case30为模板改造成8机28节点。理由很简单case30的参数是经过验证的潮流能收敛负荷水平、阻抗数量级、电压范围都是合理的。你只需要做三件事删掉两个节点、增加两台发电机、重新整理编号。在MATLAB里手动改case30.m比较麻烦更推荐用代码的方式处理% 加载自带30节点系统 mpc loadcase(case30); % 删除第29号和第30号节点以及连接到它们的支路 bus_del [29; 30]; branch_del []; for i 1:size(mpc.branch, 1) if ismember(mpc.branch(i, 1), bus_del) || ismember(mpc.branch(i, 2), bus_del) branch_del [branch_del; i]; end end mpc.branch(branch_del, :) []; mpc.bus(bus_del, :) []; % 删除与这两个节点关联的发电机组如果有 gen_del []; for i 1:size(mpc.gen, 1) if ismember(mpc.gen(i, 1), bus_del) gen_del [gen_del; i]; end end mpc.gen(gen_del, :) []; % 修改节点编号删除中间编号后会产生跳号 % 使用 ext2int 重新编号 mpc ext2int(mpc);ext2int这个函数是Matpower内部的编号统一工具可以把不连续的外部节点编号映射成1到N的连续内部编号。很多人在删除节点后不重新编号结果bus矩阵第1列出现了29、30这种空洞虽然Matpower多数情况下能处理但后续可视化、结果导出到Excel时非常容易出错。3.2 增加两台发电机删除节点后原case30系统是6台发电机离8机还差两台。加机组的时候要考虑节点类型的变化% 在第5号和第12号母线上分别加一台发电机 newgen1 [5, 40, 0, 60, -60, 1.0, 100, 1, 120, 0, ... 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]; newgen2 [12, 50, 0, 80, -40, 1.0, 100, 1, 150, 0, ... 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]; mpc.gen [mpc.gen; newgen1; newgen2];注意mpc.gen一个矩阵有21列版本不同列数可能略有差异补齐时必须把后面不用的列都填0。同时要把mpc.bus中对应节点的type从1改为2表示它变成了PV节点mpc.bus(5, 2) 2; mpc.bus(12, 2) 2;这一步很多人会忘结果新增的发电机在潮流计算里根本没被当成PV节点系统里的机组数量看着是8台实际参与电压调节的只有原来的6台。3.3 从零开始写28节点数据文件如果手里已经有文献里的系统参数不愿意从case30改也可以直接写数据文件。我的建议是先把单线图拓扑画清楚标好节点号、支路阻抗、负荷大小、发电机接入点然后按顺序填表。写之前先规划好3个问题系统的基准电压用什么等级全部用同一个电压等级最省事比如230kV否则变压器支路的变比参数要仔细算。平衡节点放在哪选一台容量大、有调节能力的机组所在节点通常放节点1。系统总负荷量级是多少自建系统建议控制在300~500MW之间太大会导致机组出力超过Pmax太小则电压普遍偏高。实际编写时数据文件的格式就是把我前面讲的三张表依次填好保存成case28_8m.m即可。写完之后用loadcase加载再调用check_mpc做基础校验这个函数会检查矩阵尺寸、节点编号、支路端点是否越界等常见问题。不同版本check_mpc调用方式略有差异但载入数据时如果数据有致命错误命令窗口通常会直接报错或者给出warning。3.4 跑通第一个潮流计算数据文件准备好了运行潮流只需要几行代码mpc loadcase(case28_8m); opt mpoption(verbose, 2); res runpf(mpc, opt);verbose设置为2可以在命令窗口看到牛顿法每次迭代的最大功率失配量。正常情况下3~5次迭代就收敛了res.success返回1。如果verbose0则只显示收敛结果不显示迭代过程。跑通之后先别急着高兴先做一次基础自检查看发电机总出力、负荷总量、网损、各节点电压幅值确认数据没有“表面收敛但结果离谱”的情况。4. 潮流结果怎么读从res结构体到可视化runpf返回的结果是一个结构体res里面包含了完整的潮流计算结果。很多新手跑完只看一个res.success就结束了这太浪费了。4.1 结果结构体的关键字段res里的核心字段和mpc是同构的也包含bus、gen、branch三张表但各表内容已经被刷新成了潮流计算后的运行状态res.bus(:, 8)各节点电压幅值pures.bus(:, 9)各节点电压相角度res.gen(:, 2)各发电机实际有功出力MWres.gen(:, 3)各发电机实际无功出力Mvarres.branch(:, 14)和res.branch(:, 15)支路首末端有功潮流MWres.branch(:, 16)和res.branch(:, 17)支路首末端无功潮流Mvar注意不同版本Matpower里branch矩阵在结果中的实际潮流列位置可能有差异。最稳妥的方式是在命令窗口输入res.branch(1, :)看一下前几列确认哪些列是有功、无功再看结果。4.2 怎么判断结果是否合理一个可用的潮流结果必须满足几个基本条件电压幅值基本在0.95~1.05 pu之间如果出现0.8或1.2这种数值说明系统无功分布有问题。相角是逐步变化的相邻节点相角差通常在5度以内出现几十度的相角差说明支路参数异常。发电机无功出力在Qmin和Qmax范围内。如果某台机组无功出力打到上限说明该节点附近无功不足系统可能处于临界运行状态。总有功出力减总负荷等于网损网损一般在总负荷的2%~5%之间。如果网损达到负荷的30%以上支路阻抗的数据几乎可以肯定填错了。判断网损在Matpower里可以直接算P_gen_total sum(res.gen(:, 2)); P_load_total sum(res.bus(:, 3)); P_loss P_gen_total - P_load_total; fprintf(系统总负荷: %.2f MW\n, P_load_total); fprintf(发电机总出力: %.2f MW\n, P_gen_total); fprintf(网损: %.2f MW (%.2f%%)\n, P_loss, P_loss / P_load_total * 100);4.3 可视化与报告导出Matpower本身不自带特别花哨的画图工具但用MATLAB基础绘图功能就足够完成大部分展示需求。我常用的三个可视化手段节点电压幅值柱状图快速判断全网电压水平是否正常。figure; bar(res.bus(:, 1), res.bus(:, 8)); xlabel(节点编号); ylabel(电压幅值(pu)); title(28节点系统各母线电压幅值); ylim([0.9 1.1]);支路负载率排序图找到系统的重载线路。loading res.branch(:, 14) ./ res.branch(:, 6); figure; bar(loading * 100); xlabel(支路编号); ylabel(负载率(%)); title(各支路负载率);注意如果rateA填了0这里的负载率全是Inf或NaN所以我前面才强调rateA一定要填。相角分布曲线能看到以平衡节点为参考的相角变化趋势。在做课程报告或者论文的时候我通常还会把res.bus和res.branch导出为Excel或者CSV。writetable(array2table(res.bus), result_bus.xlsx); writetable(array2table(res.branch), result_branch.xlsx);这样后续在Word或者LaTeX里引用数据很方便不用每次重新跑一遍再手动抄写。5. 排错实战搭建与仿真中最容易踩的坑这部分是我最想写的因为8机28节点系统本身不算复杂真正的门槛全在那些“看似不起眼的错误”上。5.1 潮流不收敛先查数据再调迭代好几次项目里我搭建的自定义系统第一次跑runpf都会碰到不收敛。Matpower的提示一般是“Newtons method power flow did not converge in 10 iterations”。碰到这个提示很多人的第一反应是增加迭代次数这完全是本末倒置。我的排查顺序是固定的检查bus表里有没有type4的隔离节点如果负荷节点被错误标为隔离节点系统会出现功率缺口。检查每个节点是否至少连了一条支路孤立节点会让雅可比矩阵奇异。检查支路x是否为0或者极小值如果某条支路电抗填了0.0001相当于一条短路支路矩阵条件数急剧恶化。检查负荷总量和发电机总出力是否匹配如果负荷远大于发电机出力上限潮流很难收敛到可行域。检查是否有节点电压初值设置得过低比如0.8pu导致迭代起点太差。有一次我排查了很久最后发现是某条支路我复制粘贴时把x0.06写成了x0.006数量级差了一个数量级潮流就怎么都算不拢。所以数据文件的校验一定要做。5.2 PV节点无功越限一个最容易忽略的细节Matpower的牛顿法在迭代过程中会考虑发电机的无功上下限。当某台发电机需要发出的无功超过Qmax时它会自动把该节点从PQ节点处理给定无功在Qmax处电压由潮流自然决定。这个机制本身是合理的但它带来的坑是你的8机系统里某台发电机可能名义上是PV节点实际上因为Qmax设置太小在结果里已经变成了无功封顶的状态电压远远低于初始设定值。很多人做最优潮流或者静态安全分析时看不到这个现象是因为默认情况下Matpower输出里只在“Generator Reactive Power Limits”部分显示Limits信息不仔细看很容易忽略。我的经验是新增发电机时Qmax/Qmin不要随意填要参考原有case30里的数值量级。比如100MVA基准下Qmax给60~80Mvar、Qmin给-60~-40Mvar是比较常见的。如果填了±9999这种极端值潮流是很容收敛但结果没有任何工程意义。5.3 节点编号与行号的混淆bus矩阵的第1列是节点编号而矩阵的行号是顺序索引。删除节点后如果不重新编号节点编号可能变成[1,2,3,5,6,7,...]这种不连续序列当你在后续程序里用res.bus(find(res.bus(:,1)5), 8)这种方式取数据时容易漏项。更麻烦的是gen矩阵和branch矩阵里引用节点编号时必须和bus表第1列保持一致。如果节点编号不连续只要有任何一个地方写错就会引用到错误的节点。所以我强烈建议在删改节点后立刻调用ext2int重新编号或者在bus表里把编号重新写成连续的1到N。5.4 支路阻抗的单位和基准问题branch表里的r、x、b都是标幺值不是欧姆也不是西门子。很多从实际工程数据换算过来的人容易把线路的欧姆值直接填进表里。标幺值和有名值的换算公式是Zpu Z_ohm * baseMVA / (baseKV^2)一个230kV、100MVA基准下的系统如果一条线路阻抗是10欧姆对应的标幺值是10 * 100 / (230^2) ≈ 0.0189 pu。如果直接把10填进去潮流结果会显示电压崩溃因为系统的电抗被放大了500多倍。5.5 建议用自带工具做一次完整检查Matpower有check_mpc这样的数据校验函数虽然不同版本调用方法不一致但它能检查矩阵尺寸、支路端点引用、重复节点编号等明显问题。我自己一般在写完数据文件后会先跑一下这个函数确认没有结构性问题后再跑潮流。另外一个笨办法是把系统节点数、支路数、发电机数打印出来和设计图对照一遍fprintf(节点数: %d\n, size(mpc.bus, 1)); fprintf(支路数: %d\n, size(mpc.branch, 1)); fprintf(发电机数: %d\n, size(mpc.gen, 1));如果数量都对不上那大概率是自己写数据文件时少了行或者多了行先在这里排查比在潮流结果里排查效率高得多。6. 从潮流到更多玩法OPF与多机系统扩展8机28节点系统的价值不止于潮流计算。把runpf换成runopf就能做经济调度意义上的最优潮流。6.1 最优潮流的快速上手调用方式非常简单res_opf runopf(mpc);但要做好OPFmpc.gencost必须填写完整否则Matpower会按默认成本计算结果可能与你预期的不符。gencost里定义的是发电机有功出力的成本曲线单位一般换算成$/h。用二次多项式表示% 三台发电机为例c2, c1, c0 分别为二次项、一次项、常数项 mpc.gencost [ 2 0 0 3 0.005 20 50; 2 0 0 3 0.010 15 40; 2 0 0 3 0.002 25 60; ];这里的MODEL2表示多项式成本模型NCOST3表示有三个系数从最高次数往常数项排列。不同版本对gencost列数的定义有细微差别我的习惯是打开自带case30.m把它的gencost直接复制过来改成自己的机组参数这样永远不会因为列数不一致报错。6.2 静态安全分析和N-1校验当系统搭好并跑通潮流后我通常还会做一轮N-1静态安全校验遍历每一条支路让它退出运行重新跑一次潮流看是否出现电压越限或支路过载。实现方式很简单n_branch size(mpc.branch, 1); for i 1:n_branch mpc_i mpc; mpc_i.branch(i, 11) 0; % 把第i条支路状态设为停运 res_i runpf(mpc_i); if res_i.success 0 fprintf(支路%d退出运行后潮流不收敛\n, i); end end这个循环虽然朴素但在8机28节点这种小系统上跑得很快是做“系统薄弱环节分析”时非常实用的工具。6.3 多机系统动态仿真的延伸如果下一步要做暂态稳定或者小干扰稳定分析8机28节点这套数据同样可以作为基础。常见做法是把构建好的mpc对象导出成PSS/E的raw文件或者直接导入到PSAT、PST等动态仿真工具里。因为在潮流层面已经验证过系统是收敛、可运行的导入到动态仿真时只需要关心动态参数惯性时间常数、励磁系统模型、调速器参数的填写不需要再回头调潮流初始值。6.4 代码包的组织习惯最后说说代码组织。我习惯把“数据文件”和“计算脚本”分开目录结构大致是这样的case28_8m.m % 系统数据文件只定义mpc结构 main_powerflow.m % 潮流计算主程序 main_opf.m % 最优潮流主程序 static_security.m % N-1校验脚本 plot_results.m % 可视化脚本数据文件里只有mpc的定义不放任何计算逻辑。主程序通过loadcase加载数据再调用Matpower的函数。这样做的好处是换一个系统数据文件主程序一行都不用改直接切换loadcase的参数就行。在实际使用中我最深的体会是Matpower项目里最花时间的永远不是runpf这一行命令而是把数据文件整理清楚。8机28节点这种自建系统关键就是先把bus、gen、branch这三张表的结构吃透后面的程序调用都是锦上添花。希望这篇分享能帮你少走点弯路。本文还有配套的精品资源点击获取
