IEEE 33节点配电网Simulink建模：数据清单、搭建步骤与调试经验

做配电网仿真这几年，我电脑里最常被翻出来的一个模型就是IEEE 33节点系统。不管是写论文、做课程设计，还是给工程项目验证算法，它都像是配电领域的一个“标准实验台”。这篇内容就围绕IEEE33节点的配电网Simulink模型，把节点数据、支路数据、数据来源和建模过程完整梳理一遍，并附上我在实际搭建中踩过的坑和解决办法。无论你是刚接触配电网仿真的学生，还是想快速搭一个基准算例做分布式电源、重构或故障分析的工程师，这份内容应该都能帮你少走不少弯路。

1. IEEE33节点系统：为什么它几乎成了配电仿真的行业默认算例

1.1 一个从论文附录里走出来的标准

IEEE 33节点系统最早出现在1989年Baran和Wu发表的论文中，原意是拿来做配电网重构和降损研究。当时作者给出了一个完整的33节点、32条支路的辐射状网架参数，包括每条线路的电阻、电抗，以及每个节点的有功、无功负荷。这个算例最难得的地方在于：规模不大不小，既不会让你看花眼，又能把电压越限、网损偏大、三相不平衡这些配电网典型问题都体现出来。

由于数据是公开的，后续做配电网分布式电源接入、故障定位、无功优化、网络重构的人，几乎都拿它当基准。我自己接触过的几十篇论文里，光是“基于IEEE 33节点”做仿真的就占了相当比例。说白了，这就是配电领域的“MNIST数据集”，大家用它不是为了它多复杂，而是为了结果可对比、可复现。

1.2 拓扑结构其实非常好记

整个系统是一个辐射状网络，电源从变电站节点1出发，通过一条主馈线逐级延伸到各个末端。节点1通常被当作平衡节点或松弛节点，电压保持在12.66 kV（标幺值就是1.0）。主馈线从节点1一直延伸到节点18，中间在节点2、3、6等处引出多条分支线，分别带到节点19到33。

除了32条常规闭合分段支路，系统还定义了5条联络开关支路。正常运行时这些联络开关是断开的，分别连接8-21、9-15、12-22、18-33、25-29。这个设计非常巧妙，因为网络重构研究的核心就是“打开哪些分段开关、闭合哪些联络开关”能让网损最小。所以你在Simulink里建模时，联络开关一定别漏掉，后续做重构会直接用到。

1.3 为什么推荐用Simulink而不是纯程序算潮流

纯MATLAB脚本算潮流当然可以，比如用Newton-Raphson法写个几十行的程序就能算完。但Simulink的最大价值在于：你可以把控制策略、电力电子变流器、保护逻辑都拖进同一个模型里跑。比如做分布式光伏接入，你可以先搭好33节点的被动网络，然后在节点18挂一个光伏逆变器模型，观察电压抬升和反向潮流。这种“网络平台+控制电路”的玩法，纯程序很难做到这么直观。

另外，Simulink里的Simscape Electrical（早期叫SimPowerSystems）提供了现成的三相电源、线路、负荷、断路器模块，搭建速度比手写微分方程快得多，而且可视化程度高，答辩、汇报、写文档都方便。

2. 节点与支路数据完整清单：来源、取值与单位

2.1 基准值与系统总览

在动手建模之前，必须先把基准值定死。IEEE 33节点系统的常用基准值是：

你可能会困惑为什么既有50 Hz又有60 Hz。原因很简单：原始Baran-Wu文献是北美系统，默认60 Hz；国内很多论文和教材在引用时会直接改成50 Hz。如果只在Phasor模式下做稳态潮流仿真，R和X是直接给定的欧姆值，频率影响不大；但如果你要算线路电感L = X / (2πf) 搭电磁暂态模型，就必须先确认自己用哪个频率。

2.2 32条常规支路的阻抗数据

以下是我整理的标准支路数据，R和X单位均为欧姆。首端节点到末端节点的方向不影响仿真结果，但建议保持一致，方便对照文献。

注意支路18是从节点2接到节点19，支路22是从节点3接到节点23，这两条是主馈线上的“T接点”，建模时千万别接错位置。很多人一不留神就把节点19接到了节点1或者节点3，导致后面的潮流结果完全对不上。

2.3 联络开关支路与节点负荷数据

5条联络支路的阻抗如下，默认断开：

节点负荷数据是另一个关键表。这里P单位是kW，Q单位是kVar。注意节点1是变电站节点，本身不带负荷。

这里有个容易忽略的点：节点30的无功是600 kVar，在有功只有200 kW的情况下，功率因数只有0.32左右，明显偏低。这是原始数据故意设置的“重无功节点”，用来考察配电系统的无功补偿能力。仿真时如果你发现节点30附近电压掉得特别厉害，不要怀疑数据错了，这是正常的。

2.4 数据来源与版本差异

这套数据最权威的来源是Baran和Wu在1989年发表的论文，文章附录直接给出了完整参数。后来MATPOWER里收录了case33bw，和原始数据基本一致，只是单位做了统一处理。很多开源仓库、中文论文里引用的也是同源数据，但需要注意三个常见差异：

第一，部分文献会把频率从60 Hz改成50 Hz，这会影响你按物理H/m去折算电感时的数值；第二，部分版本会把节点24和节点25的负荷写成别的大小，尤其是当作者为了研究分布式电源场景而主动增大负荷时；第三，个别版本对节点30的600 kVar做了调整，比如改成200 kVar。所以我建议以原始论文和MATPOWER的case33bw为准，当你从CSDN或其他地方下载模型时，最好先核对一下这两个关键点。

3. Simulink建模过程：从数据表到可运行模型

3.1 建模路线选型

在Simulink里建配电网模型，主流路线是使用Simscape Electrical库里的Specialized Power Systems。这套库提供三相电源、三相串联RLC支路、三相并联RLC负荷、三相变压器等模块，模型结构和你看到的单线图几乎一一对应，排查问题非常方便。

对于IEEE 33节点这种规模的网络，我强烈建议做“单相拓扑的三相等效”：就是用三相互感、三相线路模块搭建整个系统，所有模块都选Three-Phase类型，而不是把A/B/C三相分开建模。分开建模会让模型膨胀三倍，而且很容易因为某一相接错导致三相不平衡，纯属给自己找麻烦。

另外要提前想清楚仿真模式。如果只是算稳态潮流，用Powergui的Phasor模式就够了，速度快、结果准；如果要看故障暂态或逆变器控制动态，就必须用Discrete模式或Continuous模式。我的做法是先建一个Phasor模式的版本，把潮流结果验证通过后，再复制一份改成Discrete模式做动态分析。

3.2 电源与负荷建模细节

电源选用Three-Phase Source模块，参数设置：相电压幅值按线电压12.66 kV折算，频率按你的系统设定，内部阻抗可以设一个很小的值代表变电站变压器阻抗。节点1作为平衡节点，电压保持在1.0 p.u.。

负荷模块有两种选择。一种是Three-Phase Parallel RLC Load，直接填有功P和无功Q；另一种是Three-Phase Dynamic Load，可以设置负荷随电压变化的特性。稳态分析用第一种足够，动态分析推荐第二种。填写时单位要看清，模块默认是kW和kVar，如果是标幺值还要先除以基准值。

特别提醒：Three-Phase Parallel RLC Load模块里的Q是一个净无功，正值表示感性无功。你把节点30的600 kVar填进去之后，它会自动折算成等效电感值。如果后面要做无功补偿，可以在同一母线上再并联一个容性负荷，等效于投切电容器。

3.3 线路建模：R、X与L的换算

线路模块建议用Three-Phase Series RLC Branch，每个模块代表一条支路。R就是表里的电阻值，X对应的电感L要用公式 L = X / (2πf) 换算得到。比如支路1的X是0.0470 Ω，在50 Hz下L = 0.0470 / 314.16 ≈ 0.00015 H。这个值很小，但确实不能省。

电容C一般直接设0。原因是IEEE 33节点数据的原始模型默认忽略对地电容，尤其是中压架空线路，容性效应本来就弱。你硬填一个电容值进去，反而会让潮流结果出现偏差。

还有一个操作技巧：不要在模型里手动建32个RLC Branch然后逐个填参数，太容易错。我一般先在MATLAB工作区定义一个结构体数组，把32条支路的R、X、L算好，然后用循环语句自动生成并连接模块。虽然前期要写一点代码，但后期改数据、换场景特别方便。

3.4 测量模块与数据采集

验证模型必须测电压和功率。在每个节点接一个Three-Phase V-I Measurement，然后通过Bus Selector把三相电压RMS值取出来。如果你想在MATLAB脚本里处理结果，就把测量信号接到To Workspace模块里，仿真结束后直接在命令行分析。

功率测量稍微复杂一点，因为三相瞬时功率是波动的。Phasor模式下可以直接读模块输出的P、Q标幺值；Discrete模式下需要用Three-Phase Power模块做dq变换再取均值。我个人经验是：验证网损时直接比较电源侧输出的三相有功总和与所有负荷有功总和的差值，比逐个支路累加功率要稳定得多。

3.5 仿真参数与Powergui设置

模型里必须放一个Powergui模块，这是Simscape Electrical的心脏。Phasor模式下，你不需要设置步长细节，Powergui会把它当作稳态相量问题来处理。Discrete模式下，采样步长一般设为50微秒到100微秒，求解器选ode23tb或ode15s，这两种对电力电子和刚性系统都比较友好。

仿真时长也分场景：只算潮流，仿真时间设0.1秒都够，因为Phasor模式直接收敛到稳态；动态仿真，至少跑0.5到1秒才能看到过渡过程。别把仿真时间设太长，否则波形数据会非常大，浪费内存。

4. 模型验证：拿什么判断自己建对了

4.1 先对潮流结果，再谈其他

模型搭完之后，第一件事不是接分布式电源，而是验证基础潮流。IEEE 33节点基准场景有公认的参考值：全网有功损耗约202.7 kW，无功损耗约135.1 kVar，最低电压出现在节点18附近，大约0.913 p.u.。这三个数字就是你的“体检指标”。

我在仿真后一般这样操作：从Scope里导出各节点电压RMS值，再用MATLAB命令行计算总损耗。如果损耗偏差在1%以内，说明线路参数和连接关系基本正确；如果偏差超过5%，大概率是某条支路接错了节点，或者某个负荷填成了标幺值而没转换。先把基础场景调到接近参考值，后面所有扩展工作才有意义。

4.2 电压分布曲线的读法

把33个节点的电压画成曲线，你会看到一个典型的“下垂+锯齿”形状。主馈线从节点1到节点18，电压一路缓慢下降，在节点18处达到最低；支线上的节点33、22、25也偏低。这是因为越靠近电源，线路电流越大，压降越明显；而末端节点由于支线阻抗大、负荷重，压降容易叠加。

如果你发现电压曲线整体偏高或偏低，先检查变电站电压设置是不是1.0 p.u.；如果只是某一个支线明显异常，就检查那条支线的负荷和线路参数是否填反了。我遇到过一位同学，把节点24和25的420 kW写成了42 kW，结果整条支线电压几乎不降，和标准曲线差得很明显。

4.3 功率平衡检查

还有一个非常直观的验证方法：看电源输出的总有功功率是否等于全网负荷有功加网损。33个节点的负荷有功总和是3715 kW，如果仿真结果显示电源输出了约3917 kW，两者之差就是网损，应该和前面说的202 kW基本吻合。这本质上就是能量守恒，也是最容易定位问题的方法。

如果电源输出功率远大于负荷和损耗之和，说明模型里存在额外接地的通路，或者某个负荷设成了负极；反过来，如果电源输出功率过小，说明有支路断路或者负荷没有接入。这套查错逻辑我百试不爽。

5. 我在搭建和调试中踩过的坑

5.1 数据索引错位是最隐蔽的坑

Simulink和MATLAB的数组索引默认从1开始，但很多开源数据文件在导出时会把节点0作为电源点。如果你直接从某个CSV文件导入数据，编号-1的错位会让你所有的连接关系全部错乱。我当时就因为这事排查了两个小时：所有负荷值看起来都对，网损却奇高，后来才发现是节点号整体往前挪了一位。

解决方法是建模前先画一张节点连接图，把33个节点的上下游关系列出来，然后做一个“支路连接矩阵”。核对连接矩阵里第1行是1-2、第2行是2-3、第18行是2-19，确认无误后再生成模型。这一步看似繁琐，但省下的排查时间远超投入。

5.2 标幺值与欧姆值的混乱

Simscape电气模块虽然允许直接填欧姆，但很多参考值的潮流结果是标幺值。如果你在负荷模块里不小心用了标幺值而基准容量设成了1 MVA，那节点24的420 kW就会变成420 kW/1 MVA = 0.42 p.u.，而正确做法是420 kW/100 MVA = 0.0042 p.u.。整整差了100倍，所有电压都会严重越限。

我的习惯是在模型里统一使用物理单位：线路阻抗填欧姆，负荷填kW和kVar，测量结果再统一折算到标幺值。这样虽然界面显得“不专业”，但对接数据表最直观，也不容易出错。

5.3 联络开关必须保持断开

搭建初始模型时，5条联络开关如果不处理，默认状态下可能被当成普通支路接入系统，形成环网。环路会让潮流计算变得完全不一样，网损暴增不说，某些电流甚至会反向。特别是支路18-33和25-29的阻抗比较小，一旦闭合，环流非常明显。

如果你用的是Three-Phase Breaker模块来做联络开关，记住初始状态设为open，控制信号保持0。如果你只是想快速验证，最省事的做法是直接用线路断开，也就是不连接这5个支路。后面做重构研究时再换成Breaker，通过控制信号模拟开关动作。

5.4 求解器选错导致仿真崩溃

Phasor模式基本不会出现数值问题，但换成Discrete模式之后，如果你的步长设置太大，比如超过200微秒，很容易出现振荡甚至仿真直接报错。这时候不是模型有问题，而是数值稳定性问题。把采样步长降到50微秒，求解器选ode23tb，绝大多数情况都能解决。

另外，Zero-Order Hold和离散积分器混用也容易出问题。建议保持整个模型统一使用同一个采样时间，别在一个模块里用连续时间、另一个模块里用离散时间，否则你会在“蓝色实线”和“红色虚线”的报错提示里浪费大量时间。

5.5 负荷接地点别乱加

Simscape Electrical的三相负荷模块一般不需要额外接地，因为负荷内部已经规定了中性点处理方式。但如果你在某一节点额外添加了接地故障模块，或者把某根导线直接接地，就相当于人为制造了一个单相短路，模型会瞬间出现巨大的故障电流。做稳态仿真时，确保所有负荷都通过三相线路模块接入，母线不要随便加接地符号。

6. 从基础模型延伸：分布式电源、故障与重构

6.1 接入分布式光伏或储能

IEEE 33节点系统最经典的扩展就是在节点18、22、25、33这些电压偏低的位置接入分布式电源。你可以用Three-Phase Source模块模拟一个恒压源，也可以用受控电流源模拟光伏逆变器的并网电流。接入后你会看到两个典型现象：局部电压抬升，甚至可能超过1.0 p.u.；主馈线反向潮流增加，靠近变电站的支路电流可能减小，形成新的网损分布。

我这几年做分布式电源接入研究时，最喜欢用这个模型来对比不同接入位置的电压支撑效果。比如把同样容量的光伏分别放在节点18和节点24，电压改善幅度完全不同。这种对比在真实电网数据上反而不容易做，因为变量太多。

6.2 故障仿真与保护整定

在Specialized Power Systems里可以直接用Three-Phase Fault模块做三相短路、两相短路、单相接地故障。你把故障点设在节点8或者节点25，跑一次Discrete仿真，就能看到短路电流的暂态波形。这些波形可以用来做保护定值校验、馈线自动化策略验证。

需要注意的是，故障仿真必须放在Discrete模式下，Phasor模式看不到暂态过程。故障起始时间一般设在0.2秒，持续0.1秒，这样波形既有故障前的稳态，也有故障中的暂态和故障切除后的恢复过程。

6.3 网络重构与降损对比

有了带联络开关的完整模型，重构实验就顺理成章了。思路很简单：初始场景下打开5条联络开关，记录网损；然后按重构算法给出的方案，闭合某几条联络开关并打开对应的分段开关，再跑一次，比较网损变化。比如经典结论里，通过把联络开关8-21、9-15、12-22、18-33、25-29中的部分闭合，配合断开相应分段，网损可以从202 kW降到150 kW左右，降幅非常可观。

你甚至可以在此基础上做“故障后恢复供电”的仿真：在某条支路设置故障，然后通过联络开关倒负荷，观察非故障区域是否因为重构而恢复供电。这也是配电网自愈研究中最常被演示的场景。

最后分享一个经验

我在实际使用中发现，最顺手的工作流是这样的：先在MATLAB脚本里维护一份IEEE 33节点的数据表，包括支路矩阵和负荷矩阵，然后用脚本自动生成Simulink模型，而不是在模型界面里手动拖拽32条支路。手动搭一次模型至少要半天，脚本生成只要几分钟，而且改起参数来特别方便，换一套数据重新生成即可。

如果你只是急着交作业或者做毕设，直接下载一个别人搭好的模型然后核对数据也是可行的。但不管用哪种方式，都一定要自己动手跑一遍潮流，亲眼看一看节点18的电压有没有到0.913附近，亲眼看一看损耗是不是200多千瓦。只有自己验证过，后面接分布式电源、做故障分析、写小论文，你心里才有底。IEEE 33节点这套数据我用了很多年，最大的感受就是：资料不难找，难的是把每个细节都做对。希望这篇分享能帮你把这一步走稳。