企业官网建设流程全解析

刚才把69节点系统的数据跑通时，我其实已经踩了一路的坑。从case文件加载失败，到快速解耦算法怎么调都不收敛，再到最后发现是负荷单位出了问题，每一步都值得记录。这篇不算教程，更像一份从零到一、把IEEE 69节点配电系统用Matpower完整玩转的实操笔记，把计算原理、数据转换、算法选型、结果解读这些事一次性讲透，适合电力专业学生、刚接触配电网分析的工程师，以及想用分布式电源做扩展研究的人。

1. 为什么要玩69节点系统：教科书案例与真实配网的差距

1.1 潮流计算到底在算什么

潮流计算本质上是在解一个大型非线性方程组。电网的节点电压和注入功率不是线性关系，我们已知的是线路参数、变压器参数、负荷大小和发电机出力，未知的是每个节点的电压幅值与相角，然后求支路功率和网损。因为节点注入功率方程里全是电压乘积与三角函数的组合，没法直接求解析解，只能靠迭代逼近数值解。

理解潮流计算的价值，只要想清楚一件事：电网规划和运行中的绝大多数决策，都要建立在"知道稳态电压和潮流分布"之上。线路会不会过载、电压会不会越限、无功补偿装在哪里、分布式电源接到哪个位置最合适，这些问题的第一步都是先做一次可靠的潮流计算。没有这个基础，后续任何优化算法都像在沙子上盖楼。

1.2 IEEE 69节点系统是个什么样的"小白鼠"

IEEE 69节点测试系统是配电网领域非常经典的一套公开算例，早期主要出现在配电网重构的研究文献中，是从真实配电馈线抽象出来的标准测试网络。整套系统共69个节点、68条支路，基准电压12.66 kV，总负荷规模大约在3.8 MW + 2.7 Mvar量级，属于典型的中压辐射状配电网络。

它还有5个平时断开的联络开关，这些开关一旦闭合，就能重新组织网络拓扑。正因如此，IEEE 69节点系统成了配电网重构、故障恢复、电压治理、分布式电源接入研究最常用的通用实验平台之一。跟IEEE 14、30、118节点这些输电系统算例完全不同：输电系统是环形网络、电压等级高、R/X比值很小；69节点系统是单电源辐射状供电、R/X比值明显偏大、末端电压偏低、网损占比高，典型的"配网性格"。

1.3 为什么选它作为Matpower入门对象

做潮流计算，先选输电系统的算例上手当然也行，但配电网的很多现象在输电算例上看不到。比如电压沿馈线逐渐降低，末端节点电压可能逼近0.9 p.u.；再比如线路损耗占比明显高于输电网。69节点系统规模不大，数据量刚刚好，出问题能快速定位，同时又能体现配电网分析中的核心矛盾。

我在实际带项目的过程中发现，用这套系统做Matpower练习，最大的好处在于"能看见问题"。比如网络重构前后网损变化了多少、接入DG后末端电压抬升了多少，这些数字在69节点系统上非常直观，做课程设计、写论文、做项目预研都拿得出手。如果你已经在IEEE 14节点上跑过常规潮流，那么69节点会把你从"会按按钮"拉到"理解配网特性"这一层。

2. Matpower的数据准备：从原始数据到可计算的案例文件

2.1 认识Matpower的算例结构

Matpower的算例文件本质上是一个MATLAB结构体，主要在四个字段中存放数据：

• baseMVA：基准功率，通常单位是MVA；
• bus：母线数据矩阵，每一行代表一个节点；
• branch：支路数据矩阵，每一行代表一条线路或变压器；
• gen：发电机数据矩阵，每一行代表一台发电机或无功源。

bus矩阵的列顺序非常固定，我最常用到的列是这些：

branch矩阵的列含义里，重点是前五列：首端节点、末端节点、支路电阻R（p.u.）、支路电抗X（p.u.）、对地电纳B（p.u.）。R和X都是标幺值，必须和baseMVA对应。

2.2 单位问题是最容易翻车的坑

我在指导别人用Matpower做配网潮流时，遇到最多的错误不是算法选错，而是单位没统一。Matpower内部所有功率单位都按MW、Mvar处理，电压按p.u.处理。而IEEE 69节点系统的原始数据往往以kVA、kVAR、欧姆为单位，必须转换成标幺值。

标幺值转换的逻辑很简单，但要自己动手算一遍才记得住：选定基准功率S_base和基准电压U_base之后，基准阻抗就是：

Z_base = (U_base)^2 / S_base

以69节点系统为例，如果取U_base = 12.66 kV，S_base = 1 MVA，则基准阻抗约为160.28 欧姆。一根实际电阻为0.5欧姆的线路，它的标幺电阻就是0.5 / 160.28 ≈ 0.00312 p.u.。一个40 kW + 30 kvar的负荷，在1 MVA基准下就是0.04 + j0.03 p.u.。这些数字看着小，但放在一起才能形成正确的潮流方程。

有个常见的错误是把baseMVA设成100，还沿用从别处抄来的标幺值阻抗数据，结果负荷也按100倍放大，潮流直接不收敛或者得到一堆离谱的负电压、超大相角。49节点的数据不是输电网络，基准功率通常选1 MVA，最大也就10 MVA。拿到数据先问自己三个问题：阻抗数据是在什么基准下换算的？负荷数据是三相总功率还是单相功率？节点电压是用线电压还是相电压？这三个问题不确认清楚，后面全白跑。

2.3 转换数据的一次完整操作

我分享一个自己常用的转换套路。假设你手上有一份Excel格式的IEEE 69节点数据，其中包括节点负荷表和线路阻抗表。

第一步，在MATLAB中把Excel读入为表格，拆成bus矩阵和branch矩阵的雏形。 第二步，设置基准值：baseMVA = 1；线路电压等级统一填12.66 kV。 第三步，把所有负荷从kW、kvar换算成MW、Mvar，写入bus矩阵第3、4列。 第四步，把线路欧姆值换算成标幺值：

R_pu = R_ohm / 160.28; X_pu = X_ohm / 160.28;

第五步，在bus矩阵第2列里，把变电站出口节点设为3（平衡节点），其余节点全部设为1（PQ节点）。69节点系统正常情况下只有一个电源点，不存在PV节点。 第六步，在branch矩阵最后几列设置线路投运状态和容量限值，状态列填1表示投运，联络开关所在支路填0表示断开。 第七步，检查一遍：是否有节点编号缺失、是否有孤岛、是否有节点电压初始值全为1。

检查有没有孤岛，最直接的办法是把case文件加载进Matpower后，调用连通性检查函数，或者干脆先跑一次潮流，如果在某个节点附近报出"island"相关警告，那就说明数据拓扑一定有问题。

2.4 保存成case文件的小技巧

保存成m文件可以用Matpower自带的savecase函数。我自己更习惯直接写一个简单的脚本，把bus和branch矩阵手工拼好，然后存成case69_xxx.m。文件名只要不是Matpower内置的case名称，就不会覆盖自带算例。最后用loadcase验证一下：

mpc = loadcase('case69_xxx'); mpc.baseMVA mpc.bus(1:5, 1:10) mpc.branch(1:5, 1:5)

能正常加载，且数据格式符合预期，数据准备这一关就算过了。这个过程看起来繁琐，但在后续所有分析中每分每秒都在回报你。

3. 跑通第一次潮流：runpf的核心参数与输出解读

3.1 最精简的三行代码

数据准备好了，跑一次潮流其实只需要三行代码：

mpc = loadcase('case69_xxx'); results = runpf(mpc);

如果顺利，命令行会输出迭代信息，最后显示Converged in X.X X seconds之类的内容。results结构体里最重要的一组字段是：

• results.success：1表示收敛，0表示不收敛；
• results.iterations：迭代次数；
• results.et：计算耗时；
• results.bus、results.branch、results.gen：收敛后的节点电压、支路潮流和发电机出力。

我建议第一次跑的时候不要关输出，把完整的迭代过程看一遍。观察每次迭代的最大功率失配量是怎么下降的，这对理解牛顿-拉夫逊法的收敛过程非常有帮助。

3.2 算法选型：牛顿-拉夫逊最稳

Matpower的runpf默认算法就是牛顿-拉夫逊法。它利用功率方程对电压幅值和相角的偏导数构成雅可比矩阵，在每个迭代点线性化，然后反复修正电压。对IEEE 69节点这样的配电网，默认的牛顿-拉夫逊算法通常都能在几次迭代内收敛，不需要额外调整。

当然也可以用mpoption显式指定算法：

mpopt = mpoption('PF_ALG', 1); % 1为牛顿-拉夫逊 results_nr = runpf(mpc, mpopt); mpopt2 = mpoption('PF_ALG', 2); % 2为快速解耦 results_fd = runpf(mpc, mpopt2);

这两个算法的差异，放到下一章详细讲，因为它是配电网潮流中最容易踩中的"隐藏地雷"。

3.3 结果怎么读

读结果的核心是把results.bus和results.branch两个矩阵的列搞清楚。

results.bus的第8列是电压幅值，第9列是电压相角。把它们提取出来就可以分析系统整体的电压水平：

Vm = results.bus(:, 8); Va = results.bus(:, 9); [minVm, minIdx] = min(Vm); fprintf('最低电压: %.4f p.u.，出现在节点%d\n', minVm, results.bus(minIdx, 1));

results.branch比输入数据多了几列，其中第14到17列是首末端的有功、无功潮流，第18、19列是支路有功损耗和无功损耗。支路损耗直接累加，就能得到全网的有功损耗：

loss = results.branch(:, 18); total_loss = sum(loss);

也可以用平衡节点的注入功率减去总负荷来验证网损。69节点系统在基础状态下，网损通常在0.2 MW量级，占总负荷的比例能到5%左右，这个数字在输电系统里属于偏大的，但在配电网里很常见，也是网络重构想优化的目标。

3.4 第一个让人惊喜的发现

跑通第一次潮流后，我建议你立刻做一件事：把电压分布画出来。

plot(results.bus(:, 1), Vm, 'o-'); xlabel('节点编号'); ylabel('电压幅值 (p.u.)');

画完你会发现，69节点系统的电压呈现出明显的"沿馈线降低"趋势，某些末端节点的电压甚至会降到0.91 p.u.左右。这个现象在输电系统算例里很少见，但在真实配电网里非常普遍。看到这条曲线往下掉的时候，你会真正理解什么叫"无功不足导致电压偏低"、什么叫"末端电压支撑困难"。这就是69节点系统对初学者最大的教育意义。

4. 配电网的倔脾气：R/X比高导致的收敛问题与调试思路

4.1 快速解耦法为什么在配网里容易翻车

快速解耦潮流算法是牛顿-拉夫逊法的一种简化，它基于两个简化假设：线路电抗远大于电阻，也就是R/X约等于0；节点电压相角差很小。由此可以把有功和无功解耦，雅可比矩阵变成两个常数矩阵，迭代一次的计算量大大降低。

这套假设在输电网里基本成立，但在配电网上经常失效。69节点系统的很多馈线段电阻和电抗在同一数量级，甚至R比X还大。这时候快速解耦法的常数矩阵和真实系统偏差太大，迭代容易振荡、收敛慢，甚至直接发散。我见过不少同学用PF_ALG = 2去跑69节点，出来的结果要么报错，要么电压曲线完全失真。这不是Matpower有问题，而是算法本身的适用边界被突破了。

我在实际项目中测试过一个场景：同一份69节点数据，牛顿-拉夫逊法迭代4到5次收敛，快速解耦法则在10次迭代后仍然无法满足收敛精度。所以做配电网潮流，首选就是牛顿-拉夫逊法，或者直接用默认设置，别为了一点速度去冒不收敛的风险。

不过遇到一个特殊场景我会主动推荐快速解耦法：当网络规模特别大，单次牛顿迭代的矩阵分解开销已经难以忍受，并且你已经验证过该网络属于典型输电网结构时，快速解耦法才有实用价值。69节点根本不属于这个场景。

4.2 不收敛时的完整排查链路

如果results.success = 0，不要慌，按下面这条链路走一遍，绝大多数问题都能定位。

第一步，先区分是"数据问题"还是"算法问题"。把牛顿-拉夫逊法的迭代输出打开，观察最大功率失配量。如果迭代次数已经很多、失配量先降后升，可能是初始值离真实解太远；如果失配量一直不下降甚至出现NaN，大概率是数据错误。

第二步，检查负荷单位和基准功率。这是我最先怀疑的地方。比如原本20 kW + 10 kvar的负荷，如果填成20 MW + 10 Mvar，总负荷达到数千兆瓦，几乎所有节点电压都会被拉到很低的水平，潮流解不存在或者不收敛。检查方式很简单，把mpc.bus(:, 3:4)加起来，和已知总负荷对比，数量级不对就是单位错了。

第三步，检查平衡节点的设置。配电网只有一个变电站电源点，如果把平衡节点漏掉，或者把多个节点都设为平衡节点，网络失去参考或者产生矛盾约束，也会造成不收敛。69节点系统通常以变电站出口节点作为平衡节点，其余全是PQ。

第四步，检查支路状态。联络开关断开的支路，其BR_STATUS列要填0。如果一个开关支路忘记断开，网络拓扑变成了环，部分线路参数又不符合环网条件，也可能出现数值问题。

第五步，考虑给一个更好的初始值。Matpower默认采用平启动，也就是所有节点电压幅值取1 p.u.、相角取0度。配电系统在重负荷下，末端电压可能只有0.9 p.u.，平启动仍然是可以的，毕竟0.9离1不远。但如果网络包含很重的负荷，比如某些节点负荷标幺值超过1，平启动可能失败。可以先用一个更小的负荷运行，把结果作为下一次迭代的初始值。

% 把已有结果作为初始值 mpc.bus(:, 8) = results.bus(:, 8); mpc.bus(:, 9) = results.bus(:, 9); results2 = runpf(mpc, mpopt);

这种"先轻载后满载"的续用方法，在配电网很多重载场景下非常管用。

4.3 一个真实的调试故事

我印象比较深的一次，是有个做课程项目的同学拿了一份改动过的69节点数据，跑出来电压曲线全乱，有几处节点电压超过1.2 p.u.，还有几处直接为负。排查了半天，最终发现他把负荷功率填进bus矩阵第5、6列，也就是把负荷当成发电功率填了。Matpower的bus矩阵第3、4列是负荷，第5、6列才是发电注入，填反了，等于每个节点都凭空多出发电机，潮流自然完全失真。

这种错误其实很好理解，因为在一些教程里，bus数据表被简化成"有功、无功"两列，没有明确说是负荷还是注入。所以我的经验是：每一次从外部数据源转换到Matpower格式，都先打印出矩阵的前几行，对着列名逐列核对。宁可多花两分钟，也不要跑出一堆无法解释的结果后再返工。

4.4 收敛容差的设置经验

Matpower默认的收敛判据是功率失配量小于一个阈值，通常已经够用。如果你需要更精确的结果，可以调小容差：

mpopt = mpoption('PF_ALG', 1, 'PF_TOL', 1e-10);

调小容差会多几次迭代，但对结果影响往往只在网损的第五位小数以后。我的建议是，计算时用默认容差，发布结果时用较小的容差复核一遍，避免审稿人或导师追问时还要重新算。对于配电网分析，1e-8已经足够，再小对实际工程没有意义。

5. 潮流结果能干什么：电压分布、网损定位与DG接入评估

5.1 从电压曲线找薄弱点

跑通潮流只是起点，真正的价值在于解读结果。绘制电压分布曲线是最直观的第一步。69节点系统的电压曲线通常会呈现几个明显的"洼地"，对应供电半径较大的末端节点。这些节点就是电网的薄弱点，做无功补偿规划时应该优先关注。

如果把电压曲线和节点负荷叠加起来看，还能判断电压偏低究竟是因为下游负荷太重，还是因为线路阻抗太大。比如某个末端节点本身负荷不大，但电压偏低，说明问题出在上面几段线路，可能是线路过细、线路过长。这种判断在配电网规划中很常用，用Matpower跑一次潮流就能形成初步结论。

5.2 支路损耗定位与网络重构动机

逐个检查results.branch第18列的有功损耗，把损耗最大的几条支路找出来，也能讲述一个完整的故事。69节点系统的损耗通常集中在主干线的前几段，因为越靠近电源的线路承担着越大的传输功率，电流大，损耗自然高。

损耗最大的支路对应着网络结构的脆弱环节。如果把这个结论延伸到网络重构场景：通过联络开关改变供电路径，让一部分末端负荷改由其他馈线供电，就可能降低主干线电流、减少网损。这就是很多重构论文为什么拿69节点系统做验证的原因，因为它的损耗分布规律足够明显，任何重构策略效果都能直接通过网损下降幅度来评估。

[~, idx] = sort(results.branch(:, 18), 'descend'); % 显示损耗最大的5条支路 results.branch(idx(1:5), [1 2 14 18]);

5.3 接入分布式电源后会发生什么

这是69节点系统最吸引人的扩展玩法。把某个负荷节点换成分布式电源节点，最稳妥的方法是先把DG当成"负的负荷"处理：在原来负荷有功功率基础上减去DG出力，无功功率同理。这种方式改动最小，不需要动发电机数据矩阵，也完全符合恒功率逆变器的实际运行特点。

% 假设DG接在节点50，有功0.5 MW，无功0 mpc2 = mpc; bus_sel = 50; mpc2.bus(bus_sel, 3) = mpc2.bus(bus_sel, 3) - 0.5; % 节点负荷有功减少0.5 MW results_dg = runpf(mpc2);

跑完之后对比results和results_dg的电压曲线，你会发现DG接入点附近的电压被明显抬高，全网总网损也可能下降。这就是"分布式电源就地支撑电压、减少送端电流"的直观证据。

如果你想把DG建模为PV节点，也就是恒压源，那就需要动gen矩阵：把接入DG的节点类型改为2，设置PG为DG出力、VG为目标电压幅值，同时设置合理的无功上下限。但这种建模有个坑：DG的无功能力是有限的，如果系统为了维持目标电压需要大量无功，而gen矩阵里QMAX设得不够，潮流又会不收敛或者无功越限。所以我给初学者的建议是：先用负负荷模型体会DG接入效果，再尝试PV节点模型并观察无功越限问题，后者会加深你对"电压支撑是需要代价的"这句话的理解。

5.4 进一步可以做的分析

跑通基础潮流后，延展方向非常多。比如固定DG总容量、遍历不同接入位置，观察网损和最低电压的变化，这就成了分布式电源选址定容的简化版预实验。再比如修改联络开关组合、重新运行潮流，比较不同拓扑的网损，这就是网络重构思想的最初体验。还可以给末端节点增加并联电容器，观察无功补偿对电压的支撑作用。

做这些分析时，我推荐把runpf封装成一个函数，输入网络结构、DG配置，输出网损和最低电压，然后用一个循环去批量计算。这时候记得设置mpoption('OUT_ALL', 0, 'VERBOSE', 0)关掉全部输出，否则几百次循环下来，命令行会被刷爆，计算速度也会被拖慢。

6. 几个容易踩坑的实际经验

6.1 版本差异与函数变化

Matpower的版本迭代中，mpoption的字段名和runpf的部分行为在不同版本间有过调整。比如在某些旧版本中，选项的写法更依赖于mpoption的字符串参数；新版则统一用PF_ALG、PF_TOL这类字段。我的建议是，写代码时尽量用官方文档中当前版本对应的写法，并且把matpower的路径和版本信息记录在项目说明里。这样隔几个月回来看结果，至少能复现当时的计算环境。

6.2 别把配电系统数据当成输电系统处理

69节点系统只有12.66 kV，很多线路参数在标幺值下很小，而有功和无功负荷却可能达到同样的数量级。这导致潮流计算中无功功率对电压的影响非常显著，不能像输电网那样简单认为"无功就地平衡就行"。在分析结果时，一定要同时看有功和无功，特别是无功的流向。有些节点虽然有功送入，但无功却有功地流出，说明它附近存在较大的充电功率或容性补偿，这种细节在输电系统算例中容易被忽略。

6.3 数据来源的复核

网上流传的69节点系统数据版本不止一个，有的总负荷略有差异，有的联络开关编号不同。我在实际使用中发现，不同版本的数据跑出来的网损数字可能差百分之几，所以写报告或论文时一定要在附录里注明自己使用的是哪个版本的数据，最好附上总负荷量级。不要小看这一步，很多论文被质疑复现不了，根源就在数据版本不一致。

6.4 批量计算时的性能习惯

Matpower对69节点系统来说非常轻量，单次潮流基本是毫秒级。但在做选址遍历这类大批量计算时，有几个坏习惯会让计算变得很慢：一是在循环里反复调用loadcase加载case文件，二是不关输出，三是每次都构建新的mpoption对象。把case加载一次、选项设置一次，循环体内只修改需要修改的字段，然后复用，这是最基本的性能习惯。

mpc = loadcase('case69_xxx'); mpopt = mpoption('OUT_ALL', 0, 'VERBOSE', 0); for loc = 1:69 mpc_temp = mpc; % 修改数据... res = runpf(mpc_temp, mpopt); % 记录结果... end

这个小习惯在69节点上体现不出多大差距，但换到上千节点的配网上，差别就是几分钟和几个小时的差距。趁早养成，后面不吃亏。

7. 最后再说一点

Matpower跑通69节点系统，技术上的门槛其实很低，难点全在数据、单位和算法选择这些看似琐碎的地方。我个人在实际操作中最大的体会是：潮流计算不是"按个按钮出个结果"就结束了，每一份结果都要能回答一个问题——这个电压为什么这么低？这条支路为什么损耗最大？如果你能从电网结构和运行方式的角度解释这些现象，那么无论是做课程设计、写论文还是处理真实配电系统问题，都已经完成了最关键的一步。下一步，就可以放心地去碰重构算法、DG优化这些真正有挑战的方向了。