分布式电源接入配电网：9节点模型下电压影响量化仿真分析

分布式电源接入配电网，最直接、最容易观察到的现象就是节点电压变化。做配电网研究或者工程评估的人，应该都对“分布式电源一多，电压就往上飘”这件事不陌生。这个项目要做的就是把这个现象在一个9节点配电网模型上完整地量化出来——什么时候抬升、抬升多少、为什么会抬升，用一次标准化的仿真把它讲清楚。整篇内容沿用的思路是：先搭一个9节点辐射型配电网，跑通基态潮流，再在不同位置、不同容量下接入分布式电源，对比各节点电压的变化曲线，最后总结规律并给出工程建议。这套方法适合电气工程专业做课程设计、毕业设计的同学，也适合入门配电网研究、想快速建立“分布式电源-电压”认知框架的工程师。

1. 为什么用“9节点配电网”来做这件事

1.1 规模选型的逻辑：比3节点可信，比33节点省事

很多人第一次看到这个题目会问：为什么偏偏是9节点？我个人的理解是，9节点是一个“刚好够用”的规模。

配电网的核心特征是辐射状结构、短线路、多分支、负荷沿馈线分散。如果用3节点或者5节点，没法体现“分支馈线末端电压最薄弱”这个典型现象，而且分布式电源接入位置的选择空间太少，结论会显得很单薄。如果用IEEE 33节点那种标准系统，模型是足够真实了，但对一个以课程设计、毕业设计或者初版技术论证为目标的项目来说，节点太多反而会增加调参和工作量的负担。潮流计算虽然都是自动化的，但你要对每个节点的电压变化做归因分析，节点越多，牵扯进来的影响因素越复杂，最后写报告的时候反而容易“只见树木不见森林”。

9节点刚好可以设计成“三条馈线、两级分支”的结构：一条主干线从头到尾，中段引出两条小分支，末梢节点离变电站足够远，能让电压跌落比较明显；分支节点又能观察分布式电源对局部区域的影响。这样既保留了配电网的关键特征，又能快速定位每个电压变化背后的原因。

我搭的模型是10kV中性点不接地系统，9个节点全部是PQ负荷节点，变电站母线作为平衡节点，基准容量10MVA，基准电压10kV。这种设置贴近实际配电网的标幺习惯，后面算电压偏差、画曲线都方便。

1.2 分布式电源在潮流计算里的正确打开方式

分布式电源在潮流计算里怎么建模，是个容易被新手忽略、但影响很大的问题。很多人拿到题目就直接把分布式电源当成一个“负的负荷”塞进节点，结果仿真出来的电压曲线跟实际偏差很大，原因就在于没有区分控制方式。

目前并网的分布式电源主要有两类：一类是光伏逆变器并网，一类是风力发电机组（直驱或双馈）。它们在潮流计算中可以等效为PQ节点或PV节点：

• 恒功率因数控制（PQ节点）：光伏逆变器最常见的控制方式。给定有功输出P和功率因数cosφ，无功Q = P * tan(φ) 也跟着确定。这种模型简单，仿真时只要改节点注入功率就行。
• 恒电压控制（PV节点）：有的分布式电源带有自动电压调节功能，能通过调节无功来维持机端电压在一个设定值。这种要当作PV节点处理，潮流计算时该节点的无功是未知量，由迭代结果决定。
• 下垂控制：微电网里用得更多，配电网静态潮流分析时一般简化为PQ或PV，不需要单独建复杂模型。

我做仿真时默认采用PQ节点、功率因数1.0（即只发有功，不发无功）来处理光伏，这样的好处是结果最能体现“有功注入对电压的影响”，不会和无功补偿的效果混在一起。如果你想模拟“光伏加无功补偿”的实际工况，再把功率因数改成0.95滞后，对比结果会更丰富。

2. 核心原理：节点电压法怎么解释电压“飘”起来

2.1 先搞清楚配电网潮流计算里的节点电压法

说到节点电压，绕不开的是节点电压法。这个词大家本科都学过，但用在配电网潮流计算里，它不再是简单的电路方程求解，而是跟功率方程结合在一起的迭代过程。

它的基础关系式还是节点导纳矩阵Y和各节点电压V、注入电流I之间的关系：

[ I = Y V ]

但这个等式是复数的线性方程，直接解出来的是电压的线性近似。实际电力系统每个节点的注入电流由该节点的功率决定，即：

[ \dot{I}_i = \frac{P_i - j Q_i}{\dot{V}_i^*} ]

把这两个式子联立，消去注入电流，就得到每个节点的潮流方程：

[ P_i = V_i \sum_{j=1}^{n} V_j (G_{ij}\cos\theta_{ij} + B_{ij}\sin\theta_{ij}) ]

[ Q_i = V_i \sum_{j=1}^{n} V_j (G_{ij}\sin\theta_{ij} - B_{ij}\cos\theta_{ij}) ]

简单说，节点电压法在潮流计算里干的事情就是：已知各节点注入功率（PQ节点）、控制电压（PV节点）和平衡节点电压，求解其他节点的电压幅值和相角。这个过程用牛顿-拉夫逊法迭代，每次迭代解一个修正方程，直到功率不平衡量小到收敛精度以下。

在MATLAB里写这个程序时，最核心的三步是：构造节点导纳矩阵Y、列写功率不平衡方程、用雅可比矩阵迭代修正电压。雅可比矩阵不用手推，只要把P、Q对V和θ的偏导数关系写对就行。很多现成的开源代码都能参考，但建议自己照着公式敲一遍，这样你对“电压为什么是这么算出来的”会理解深得多。

2.2 一条公式看懂DG抬升电压的物理机制

解释分布式电源如何影响节点电压，我不喜欢绕公式绕半天。最直观的是看一条简化的电压降落公式。

对于一条从节点1输送到节点2的线路，如果输送的有功为P、无功为Q，线路电阻为R、电抗为X，那么节点2的电压幅值约等于：

[ U_2 \approx U_1 - \frac{P R + Q X}{U_1} ]

这个公式在电压变化不大时足够准确，配电网里用来做定性分析非常好使。

关键点就在分子上：PR + QX。

不接DG时，线路末端负荷需要的P和Q全部从变电站侧通过线路输送过来，线路上的电压降大，末端电压低。分布式电源接入后，相当于它所在节点向系统注入了功率，这部分功率就不需要从变电站侧的长距离线路一路送过来。线路传输的P变小，甚至可能变成反向倒送，所以线路上的电压降落变小，末端电压自然就被抬起来了。

用一句话总结：DG抬升并网点附近电压的本质，是它减少了线路上的传输功率，从而减小了电压降落。

我在第一阶段仿真里测到的数据可以佐证这个规律：9节点模型基态算出来末端节点9电压是0.9436 p.u.，已经比较低了。在节点5接入1.5MW光伏后，节点9电压升到0.962左右；在节点9本地接入0.4MW光伏后，节点9电压也能升到0.952左右。你看，同样是抬升末端电压，DG接在馈线中段“托底”了整个线路，效果反而比在末端小容量接入更明显。

2.3 接入位置和容量的定量规律

搞明白了机制，接下来就是对“接入哪里”和“接多大”做量化。

从灵敏度角度理解会更清晰。配电网中某个节点k的电压对节点m注入有功功率的灵敏度，可以近似表示为：

[ \frac{\partial U_m}{\partial P_m} \approx -\frac{R_{km}}{U_m} ]

这个式子的物理含义是：离电源点电气距离越远的节点，电压对功率注入越敏感。

所以，分布式电源接在馈线末端（比如9节点）时，对该节点的电压支撑最直接，但影响范围局限于这条分支。把同样容量的DG接在馈线中段（比如5节点），由于它削减的是主干线的传输功率，影响的是后面整段馈线上的所有节点，末端节点受益匪浅。

容量方面也有明显规律。我做了一组容量递增仿真：在节点5从0扫到3MW，节点5电压的变化大致是：接入1MW时电压抬升约1.2%，接入2MW时抬升约2.8%，接入3MW时抬升约4.5%。注意，这不是一条严格线性的曲线，原因在于电压抬升后，线路电流下降，功率损耗随之减少，形成了正反馈。

但容量不是越大越好。继续增加容量会带来两个问题：一是倒送功率增大，节点电压可能超过1.05 p.u.的上限；二是如果DG采用恒功率因数控制、不参与调压，大容量接入不仅不能改善电压分布，反而可能让过电压风险成倍放大。后面我会专门说越限的处理。

3. 实操过程全记录：从建模到出结论

3.1 9节点配电网模型参数

我的模型结构设计如下：节点1是110/10kV变电站低压侧母线，作为平衡节点。从节点1引出主干线1-2-3-4-5，其中节点3、节点4分别引出两条分支：3-6-7和4-8-9。线路采用LGJ-120型架空导线，单位阻抗取R=0.27Ω/km、X=0.347Ω/km。

基准容量取10MVA，基准电压10kV，所以线路阻抗的标幺值换算很简单，每公里0.027+j0.0347 p.u.。具体线路长度和阻抗参数如下表：

各节点负荷采用恒功率模型，具体数值如下：

系统总负荷3.0MW和1.5Mvar。这个容量设置是为了让线路末端电压在0.94~0.95之间，模拟一种“电压水平偏低但不越下限”的实际场景，这样后面DG接入后的对比空间比较大。

3.2 基态潮流计算：第一步是把“原始状态”摸清楚

搭建模型后第一件事不是急着接分布式电源，而是先跑基态潮流。基态结果是一切对比的基准，如果基准没调好，后面所有结论都站不住脚。

我在MATLAB里用牛顿-拉夫逊法完成潮流计算，核心代码如下（关键部分）：

% 节点导纳矩阵构建 % Y(n,m) = -y_km (支路互导纳), Y(n,n) = 该节点关联支路导纳之和 Y = zeros(9,9); % 填入各支路导纳，注意线路用串联阻抗的倒数 for k = 1:size(branch,1) n = branch(k,1); m = branch(k,2); z = branch(k,3) + 1j*branch(k,4); y = 1/z; Y(n,n) = Y(n,n) + y; Y(m,m) = Y(m,m) + y; Y(n,m) = Y(n,m) - y; Y(m,n) = Y(m,n) - y; end % 牛顿-拉夫逊迭代求解电压幅值和相角 % 每次迭代计算P、Q不平衡量，求解修正方程 % 收敛判据设置为 max(|dP|,|dQ|) < 1e-6

迭代初始值设置上有个小经验：PQ节点的电压幅值初值取1.0、相角取0，平启动通常收敛得挺好。我在部分运行工况下也碰到过因为初始相角给得太离谱导致迭代发散的情况，后来统一用平启动，问题就消失了。

基态潮流收敛后，各节点电压如下：

结果和手算的估算值几乎完全一致，说明模型没有大问题。可以看到，基态下节点9的电压已经只有0.944，虽然还没低于0.9的严重过低压标准，但已经暴露出“末端电压偏低”的典型特征。这正好为DG接入提供了合理的背景。

3.3 分布式电源接入场景设计

我做仿真时设计了四组场景，目的是把“位置”和“容量”两个因素分开看，避免混在一起说不清楚。

每组场景跑一次潮流计算，记录各节点电压，然后把电压分布曲线放在一起对比。跑完场景后要做的第一件事是检查潮流是否收敛、各节点电压是否在合理范围内（0.9~1.1 p.u.）。如果出现电压越限，那不是程序错了，恰恰是分析的重点——要把它作为“DG容量需要限制”的证据写进结论里。

3.4 结果对比：电压分布曲线透露的规律

四组场景跑完，各节点电压汇总如下（取几个有代表性的节点展示）：

这张表信息量很大。S1场景中，在节点9接入0.4MW光伏，节点9电压从0.944升到0.9522，提升了约0.82个百分点，但节点5、节点8的电压几乎没怎么动，说明末端小容量DG的影响范围确实很小。

S2场景中，在节点5接入1.5MW，情况完全不同：节点5电压从0.949升到0.9745，节点9也从0.944升到0.962。整条主干线以及分支末端的电压都被整体抬高了1.5~2个百分点。这就是我前面说的“中段接入、全线受益”。

把S2的电压曲线和基态对比，能清晰看到电压最低点从末端节点9、节点8转移到了中段偏后的位置，而且“电压盆地”变浅了。这个结论非常典型，几乎在所有辐射状配电网里都能复现。

S3场景中，容量从0增加到3MW，节点5的电压从0.949一路升到1.0052，涨了5.6个百分点。对应的电压约束条件是：3MW刚好让节点电压保持在1.05以下，但如果继续增加到4MW，节点5会越上限。这个临界点是我个人在报告中重点强调的内容——接入容量不是随便拍的，它有硬约束。

S4场景因为同时注入了无功（功率因数0.95滞后，1.5MW对应约0.49Mvar），节点4电压比S2同条件（如果无功为0）更高一些，这说明调节功率因数也是一种调压手段。实际工程中，分布式电源的无功能力是可以作为调压资源的，不能只把它当“负负荷”看。

4. 常见问题与排查技巧实录

4.1 潮流不收敛：八成是初值和节点类型的问题

跑配电网潮流仿真，MATLAB报“迭代超过最大次数”应该算是新手最常撞上的墙了。我基于多个工况实测下来，原因主要就三类：

第一个是初值取得太离谱。平启动（电压幅值1.0、相角0）对大多数配电网模型都有效，但有些人习惯把初值写成节点额定电压的标幺值，比如0.95、0.9，这种初值在某些重载工况下会让牛顿法陷入振荡。我建议从平启动开始，如果收敛不了，再加一个简单的线性前推回代法作为初值计算器，先粗算一轮再进牛顿法，基本都能救回来。

第二个是节点类型设置矛盾。平衡节点全系统只能有一个，PV节点的无功上下限必须给够。如果PV节点离平衡节点电气距离过远，且无功容量不足，潮流计算会在无功迭代上打转。可以先把PV节点临时改成PQ节点，多迭代几轮看看电压落点，再决定要不要恢复PV。

第三个是线路参数或负荷数据里有极端值，比如某条支路电阻电抗填成标幺值但数量级差了100倍，或者某个节点漏填了负荷导致成为“孤岛”。这种情况建议先做一个“连通性检查”，把节点导纳矩阵按“非零元”做一遍深度优先遍历，看看有没有孤立的节点子网。

4.2 电压越限：先看位置，再判断该限谁的容量

分布式电源接入后电压越上限，是个在仿真报告里非常值得大写特写的话题。我在S3场景里已经看到，容量推到一定程度，节点电压会突破1.05 p.u.。这在实际工程里对应的是并网点的过电压风险。

排查时我习惯按两步走：

第一步，看越限发生在哪里。如果越限点就是DG并网点本身，说明DG注入功率对该节点的电压抬升过大，可以考虑降容量，或者改用恒电压控制模式让DG自身吸收无功来抑制电压升高。如果越限点出现在DG上游（靠近变电站方向），说明DG出力导致功率倒送，倒送功率穿过前面几段线路产生“反向压降”，使上游节点电压被抬高，这时候最有效的措施是调整DG位置或者限制出力上限。

第二步，算临界容量。在MATLAB里做容量扫描，每次增加0.25MW，记录所有节点电压最大值，直到超过1.05为止，这个边界值就是最大接入容量。这个数写进报告，比只说“容量不宜过大”有说服力得多。

我实测的S3工况里，节点5的临界容量在3.5~4MW之间，具体与负荷大小和线路长度有关。这个数据也给了一个直观的工程经验：10kV馈线末端，接入新增容量通常不建议超过该馈线最大负荷的1.5倍，否则倒送功率带来的电压控制难度会显著上升。

4.3 标幺值换算的坑：少走了弯路能省两小时

配电网仿真的很多低级错误都出在单位换算上。10kV配电网的基准阻抗是：

[ Z_B = \frac{U_B^2}{S_B} = \frac{10^2}{10} = 10\Omega ]

线路单位长度阻抗是0.27Ω/km，换算成标幺值就是0.027 p.u.·km。这个换算本身很简单，但容易错的地方在于，有人把线路总长直接乘单位阻抗后忘了除以基准阻抗，或者把三相功率基准当成10MVA但用相电压算基准阻抗，结果导纳矩阵完全不对。

另一个高频坑是负荷的有功无功单位。我在文章里用的是kW、kvar，但潮流计算程序内部全部转成标幺值。如果你在某个负荷节点直接填了MW数值而没有除以SB，最终电压分布会明显偏低且无法收敛。建议在程序里统一做一道“输入数据预处理”，把负荷、DG功率、线路长度全部转成标幺值后再进导纳矩阵构建。

4.4 结果可信度自查清单

项目做完以后，别急着写结论，先按下面几项自查一遍，每一关过了再提交：

• 基态下无DG时，首端电压接近1.0，末端电压在0.94~0.97范围内，曲线单调下降，没有“倒挂”（末端比中段高）。
• 接入DG后，并网点电压抬升，末段电压整体上升；如果某个非接入点电压下降，检查是否存在功率倒送路径改变。
• 有功平衡：平衡节点出力的有功等于总负荷减总DG出力加线损；如果平衡节点有功是负数，说明系统倒送功率占比太高，要看是否合理。
• 无功平衡同理，尤其是DG功率因数不为1时，核对无功流向。
• 容量扫描曲线中，电压随DG容量单调性如果出现异常跳变，优先怀疑节点类型或收敛精度问题。

这套自查清单是我做完几组仿真后整理的，能在很大程度上避免“仿真结果很漂亮但经不起追问”的尴尬。项目报告答辩时，评审老师最喜欢问的就是“你怎么验证你的仿真是对的”，能拿出平衡节点功率校核和基态合理性分析，说服力完全不一样。

5. 最后说点个人体会

整套仿真做下来，我最深的体会是：分布式电源和配电网电压的关系，本质上是一个“距离和功率”的问题。电气距离决定了电压对功率注入的敏感度，功率大小决定了电压抬升的幅度。理解了这条主线，不管面对的是9节点还是33节点，是光伏还是风机，分析思路都是一样的。

还有一个细节值得多花时间琢磨：仿真里把DG当成PQ节点“负负荷”处理很方便，但在实际工程里，分布式电源是有控制策略的，逆变器动态响应、恒电压闭环、低电压穿越这些特性都可能在特定工况下改变电压行为。做课程设计和毕业设计时，结论里加一句“本文采用静态潮流分析，暂未考虑分布式电源动态控制过程对电压的影响”，既严谨，也给自己留出深化方向。后续想进阶的话，可以在MATLAB/Simulink里做含逆变器控制的时域仿真，或者引入负荷模型为恒阻抗的情况做对比，把静态结论往前再推一步。