1. 先把问题拆清楚分布式电源为什么会让配电网“不舒服”1.1 核心矛盾配电网从“单向被动”变成“双向主动”传统配电网的设计思路特别清晰变电站是唯一电源电能顺着馈线往用户侧单向流动电压沿着线路一路递减继电保护、无功补偿、线路选型全部基于这个“单向”前提。但分布式光伏、风电、储能大量接入后配电网的角色从“一根管子只往一个方向送水”变成了“中间有人往回倒水”潮流方向不再固定电压分布也不一定是单调递减了。我做项目时喜欢用一个比喻原来整条馈线就是一条下坡路水从高处往低处流你站在路上看越往下游水位越低现在半山腰突然出现一个泵站往上打水下游的水位可能比中游还高甚至会倒灌到上游。分布式电源在配电网里的作用就有点像这个“半山腰的泵站”——它能抬升局部电压也可能让整段线路的潮流方向反过来。这就是为什么要专门研究“分布式电源接入对配电网的影响”。不是简单地把DG看成“某个节点少用了一点电”就完事了而是要从电压质量、网络损耗、线路负载、保护配合等多个维度重新审视配电网的运行状态。本文用Matlab实现了一套完整的分析流程能帮你定量回答三个问题DG接入后电压到底偏了多少网损是降了还是升了这个DG容量和接入位置到底合不合理。1.2 影响分析的三个主要维度做这类研究我习惯先把影响面拆开这样后面仿真才有明确的目标。配电网领域最关心的指标一般集中在三个方面电压质量各节点电压幅值是否越限电压偏移率多少。中压配电网国标一般要求电压偏差不超过标称电压的±7%实际工程里还要留裕度。网络损耗系统总有功损耗、线损率。这直接关系到DG接入的经济账也是规划人员最敏感的指标之一。设备负载与潮流方向重点线路的载流量、是否有功率倒送。倒送功率会影响线路保护定值、变压器分接头甚至计量装置。如果只算一个潮流就结束那论文或报告根本没有说服力。完整的做法是先在无DG的基准场景下跑一遍再改变DG的接入位置、接入容量、功率因数三个变量分别对比电压分布曲线和网损曲线最后总结规律。我们这套Matlab实现就是沿着这个思路搭的。1.3 这套Matlab实现能做什么代码实现部分完成了几件具体的事搭建了IEEE 33节点典型辐射状配电网模型实现了基于前推回代法的三相平衡潮流计算考虑分布式电源作为有功无功注入源支持任意指定DG接入节点、容量和功率因数自动输出各节点电压幅值、支路潮流和全网有功损耗并提供了画电压分布曲线、网损随容量变化曲线的一键脚本。这套代码适合正在做课程设计、本科毕业设计或者刚接触配电网仿真的人。如果你在电力公司、设计院或者新能源单位做接入方案校核也可以把节点数据和负荷数据替换成自己负责的实际馈线结构前推回代法在辐射状网里照样能跑只是数据整理的工作量多了一些。2. 配电网模型与分布式电源建模2.1 网络建模为什么选IEEE 33节点仿真分析的第一步是选一个能代表实际配电网特性的测试系统。我选的是IEEE 33节点标准算例这也是配电网研究里被引用最多的馈线模型原因有三个第一它是典型的辐射状结构33个节点、32条支路包含多条分支线很贴近实际的中压馈线形态第二基准数据公开电压等级12.66kV基准容量10MVA总负荷约3715kW加2300kvar各大论文里都能查到第三它的末端电压在无DG时本身就偏低通常在0.93左右正好方便观察DG接入后电压抬升的效果。搭建模型时输入数据用每行一条支路的方式组织格式是[首端节点号末端节点号电阻(Ω)电抗(Ω)]节点负荷则用[节点号有功(kW)无功(kvar)]表示。IEEE 33节点系统的支路阻抗数据和负荷数据在网上能搜到这里只摘几条说明格式% 支路数据片段1-2, 2-3, 3-4 ... branch [ 1 2 0.0922 0.0470; 2 3 0.4930 0.2511; 3 4 0.3660 0.1864; ... ];这里有个很关键的细节容易踩坑阻抗单位是欧姆还是标幺值不同资料写法不一样。我在代码里统一采用实际值输入然后在程序里做标幺化处理基准阻抗由电压和容量推出当U_base取12.66kV、S_base取10MVA时Z_base 12.66² / 10 ≈ 16.02Ω。这么处理的好处是以后换成实际馈线数据时只需要修改支路阻抗和负荷表主程序完全不用动。2.2 DG在潮流计算里的等效模型分布式电源在电力系统分析里怎么建模是个很讲究的问题。光伏、风电、储能、燃气轮机物理特性完全不同但在稳态潮流计算这个层面我们关心的只是它们向电网注入的有功和无功所以一般按节点类型来处理恒功率控制的分布式光伏逆变器和储能变流器在潮流里最常处理成PQ节点。有功出力固定比如光照稳定时就是最大功率跟踪出来的值无功按功率因数或运行策略给定。如果逆变器参与了电压控制比如恒电压控制或无功下垂控制那就不能简单当成PQ节点要按PV节点处理潮流算法需要增加无功修正环节。大型风电场如果涉及异步发电机组还需要考虑无功吸收特性但配电网层面的小型分布式风电一般也用恒功率因数模型。这套Matlab实现默认把分布式电源处理为PQ节点节点注入功率为有功加无功的复数形式并且按“电源注入为负、负荷消耗为正”的统一约定参与潮流计算。这样处理的好处是简洁、稳定、好理解也覆盖了大部分配电网稳态分析场景。如果你确实需要模拟DG的电压控制特性可以在潮流迭代循环里加一个无功修正函数这就是后话了。2.3 负荷与基础场景设置负荷大小和分布直接决定潮流结果这部分不能随便填。我在这套代码里保留了IEEE 33节点标准负荷数据分布在节点2到33每个节点的负荷大概是几十到几百千瓦不等。实际做研究时你可以按项目需要把负荷整体放大或缩小也可以把某条分支线上的负荷集中到一个节点只要保证总负荷量级合理就行。DG接入场景我建议按“对照实验”的方式设计否则结果容易乱。以这个项目为例我设计了这样几组场景无DG的基准场景DG接在馈线首端比如节点2、中段比如节点18、末端比如节点33保持容量不变对比位置影响DG接在同一个位置容量从200kW逐步增加到1200kW观察容量变化的影响。每组场景跑完记录最大电压偏差、系统总网损、最小节点电压三个指标最后整理成表格或曲线这种“控制变量”的方法写进报告里逻辑清楚评审一看就懂。3. 基于前推回代法的潮流计算与Matlab实现3.1 为什么选前推回代法配电网潮流算法有很多种牛顿-拉夫逊法、PQ分解法、前推回代法都能算。我写这个项目时选了前推回代法原因是它和配电网的辐射状结构天然匹配实现简单、收敛性好、内存占用低。牛顿类方法虽然通用性强但需要形成雅可比矩阵对刚接触配电网计算的人不太友好而且配电网R/X比值比较大PQ分解法还可能遇到收敛性问题。前推回代法本质上就是反复利用两个基本关系一是节点功率平衡KCL二是线路电压降落公式。它不涉及矩阵求逆每次迭代只做标量复数的加减乘除计算效率非常高。对IEEE 33节点这种规模收敛通常只需要4到8次迭代速度是毫秒级的。如果你后面要做蒙特卡洛抽样或者大量场景遍历用这个方法能把计算时间压缩很多。3.2 算法流程和数学原理前推回代法每一步分两段前推过程从线路末端向电源端推进。已知各节点负荷和DG注入功率再结合上一次迭代得到的节点电压从最末端的支路开始逐条计算流经支路的复功率。某条支路的末端功率等于其末端节点的净注入功率加上从这个节点继续往下游流出的所有子支路功率之和然后再叠加支路本身的功率损耗。回代过程从电源端向末端推进。已知始端节点电压变电站母线电压设为恒定结合前推算出的支路功率用电压降落公式逐条支路计算下一节点电压。这里要特别留意支路功率、支路阻抗都用复数表示电压降落包含纵分量和横分量。收敛判断比较本轮与前一轮各节点电压幅值差取最大值小于给定阈值比如1e-6则认为收敛否则带着新电压进入下一轮迭代。用公式写支路损耗近似为S_loss (|S_end| / |V_end|)² × (R jX)电压降落近似为V_end V_start - (P·R Q·X) / |V_start| - j(P·X - Q·R) / |V_start|这里的标幺值计算很容易出错我在代码里专门做了单位统一阻抗除以基准阻抗功率除以基准功率电压本身就是标幺值根节点电压恒为1.0∠0°。3.3 核心代码实现整个程序拆成两个文件思路更清晰。一个是潮流计算函数负责给定网络数据和DG配置后返回电压、支路功率、网损另一个是主脚本负责定义系统数据、循环遍历不同DG场景、画图。潮流计算函数的核心部分长这样function [V, S_branch, P_loss_total] solve_powerflow(branch, R_pu, X_pu, load_data, dg, S_base) % 前推回代法求解配电网潮流 % 输入: % branch: n_branch x 2, 每行[首端节点, 末端节点] % R_pu, X_pu: 支路电阻、电抗标幺值 % load_data: [节点, P(kW), Q(kvar)] % dg: [节点, P(kW), Q(kvar)] % S_base: 基准容量(kVA) % 输出: % V: 节点电压复数标幺值列向量 % S_branch: 各支路首端功率标幺值列向量 % P_loss_total: 全网有功损耗(kW) n_node max(branch(:)); n_branch size(branch, 1); % 建立父子关系默认每条支路第一个节点是靠近电源侧的父节点 parent zeros(n_node, 1); child_of cell(n_node, 1); for k 1:n_branch s branch(k,1); e branch(k,2); parent(e) s; child_of{s}(end1) e; end % 初始化净注入功率负荷为正DG注入为负 S_net zeros(n_node, 1); for i 1:size(load_data,1) bus load_data(i,1); S_net(bus) S_net(bus) (load_data(i,2) 1j*load_data(i,3)) / S_base; end for i 1:size(dg,1) bus dg(i,1); S_net(bus) S_net(bus) - (dg(i,2) 1j*dg(i,3)) / S_base; end % 平启动初始化 V ones(n_node, 1); V(1) 1.0; % 平衡节点电压 max_iter 100; tol 1e-6; for iter 1:max_iter V_old V; % ---------- 前推: 从末端向首端计算支路功率 ---------- S_branch zeros(n_branch, 1); for k n_branch:-1:1 e branch(k,2); % 末端节点总流出功率 节点净注入 所有子支路首端功率 S_end S_net(e); for ch child_of{e} idx find(branch(:,2) ch, 1); S_end S_end S_branch(idx); end % 叠加支路损耗(近似用上一轮电压幅值) Vm abs(V(e)); Zk R_pu(k) 1j*X_pu(k); S_branch(k) S_end (abs(S_end)/Vm)^2 * Zk; end % ---------- 回代: 从首端向末端更新节点电压 ---------- for k 1:n_branch s branch(k,1); e branch(k,2); Zk R_pu(k) 1j*X_pu(k); V(e) V(s) - conj(S_branch(k)) * Zk / conj(V(s)); end V(1) 1.0; % 保持平衡节点电压 % ---------- 收敛判断 ---------- if max(abs(abs(V) - abs(V_old))) tol break; end end % 计算全网有功损耗(kW) P_loss_total 0; for k 1:n_branch e branch(k,2); Vm abs(V(e)); Ploss_pu (abs(S_branch(k)) / Vm)^2 * R_pu(k); P_loss_total P_loss_total Ploss_pu * S_base; end end这段代码有几个地方值得解释。第一前推时用abs(S_end)/Vm算电流幅值进而求支路损耗这是前推回代法最常见的简化处理迭代收敛后误差很小。第二回代更新电压用了共轭形式这是严谨的复功率电压降落公式如果偷懒写成V(e) V(s) - S_branch(k)*Zk/V(s)在轻载时差别不大但在重载或DG渗透率高的场景可能会偏差几个百分点。第三收敛判据我取的是电压幅值差不比较相角因为配电网分析更关心幅值如果要做严格的潮流比较可以考虑同时比较实部和虚部。主脚本负责数据定义和场景遍历%% 主脚本DG接入影响分析 clc; clear; close all; % 系统基准值 U_base 12.66; % kV S_base 10000; % kVA Z_base U_base^2 / (S_base/1000); % 欧姆 % 这里填完整的IEEE33支路数据 branch [...]; % 每行[首端, 末端, R(ohm), X(ohm)] Z_pu (branch(:,3) 1j*branch(:,4)) / Z_base; R_pu real(Z_pu); X_pu imag(Z_pu); % 完整负荷数据 load_data [...]; % 每行[节点, P(kW), Q(kvar)] % 场景1: 无DG基准 V0 solve_powerflow(branch, R_pu, X_pu, load_data, [0 0 0], S_base); % 场景2: DG接在节点18容量从200kW到1200kW dg_bus 18; cap_list 200:200:1200; V_dg zeros(length(branch)1, length(cap_list)); loss_dg zeros(size(cap_list)); for i 1:length(cap_list) dg [dg_bus, cap_list(i), 0]; % 单位kW功率因数1 [V_dg(:,i), ~, loss_dg(i)] solve_powerflow(branch, R_pu, X_pu, ... load_data, dg, S_base); end % 画电压分布曲线 figure; plot(0:33, abs(V0), k-o, LineWidth, 1.5); hold on; plot(0:33, abs(V_dg(:,1)), b--s); plot(0:33, abs(V_dg(:,end)), r--^); xlabel(节点编号); ylabel(电压幅值(标幺值)); legend(无DG, DG 200kW, DG 1200kW); grid on;这段主脚本只是骨架实际使用时你还需要把IEEE33的完整支路和负荷数据填进去。我强烈建议第一次运行时先用无DG场景和已发表论文的基准结果做对比确认自己的潮流结果和文献一致比如末端节点电压大约在0.90到0.93之间、总网损大约200kW左右。基准对不上后面一切分析都没有意义。3.4 结果输出与画图Matlab仿真的好习惯是把所有结果组织成结构体或者表格方便后续统计分析。我一般会在主脚本结束位置加一段汇总代码把每组场景的最大电压、最小电压、网损率、倒送功率等关键指标存到一个表里然后用writetable导出成Excel这样论文写作时可以直接引用数据。画图方面最常用的三张图是不同DG接入方案下的电压沿馈线分布曲线系统总有功损耗随DG容量变化的曲线不同接入位置对末端节点电压提升效果的对比柱状图。这三张图基本就能覆盖一篇报告里需要展示的大部分信息。画图时注意用不同的线型或颜色区分场景图例、坐标标签、单位都要给全否则过两天自己回来看都容易忘。4. 影响评估仿真结果怎么看才有意义4.1 无DG基准潮流结果算完无DG场景后你手里就有了一把“比较的尺子”。IEEE 33节点系统在标准负荷下电压沿馈线从根节点1.0开始递减最末端的节点电压通常会跌到0.90到0.93之间系统总网损大概在180到220kW的区间。这个末端电压其实是偏低的接近运行下限恰好说明配电网有提高电压质量的需求DG接入可以从某个角度缓解这个问题。我习惯先画基准场景的电压分布曲线再把后续每个方案的曲线叠上去一张图就能看出DG对电压的整体影响。如果你算出来的无DG末端电压偏离这个区间太多先别急着分析DG回头检查负荷数据和支路数据是不是填错了或者看回代公式里的共轭有没有处理好。4.2 DG接入对电压分布的影响规律DG接入后最直观的变化就是电压被“抬”起来。用前推回代法多次实验后你能看到非常明显的规律接入点离根节点越远对电压抬升的拉动作用越明显接入点附近节点电压抬升最大远离接入点的上游节点影响逐渐减弱。原因是线路压降和输送功率成正比当DG就地供给了一部分负荷线路上流过的功率减小压降自然变小越靠近DG接入点的支路这种“减载”效果越显著。值得一提的是DG容量不能无限加。当某个节点DG注入功率超过该节点以下区域总负荷时线路上会出现功率倒送潮流反向导致DG接入点附近的电压可能被抬高到接近甚至超过1.05标幺值。我在仿真里把DG接到节点33线路末端并逐步加大容量时电压几乎呈线性上涨1200kW时电压已经顶到1.08以上直接越限。这种“过电压”现象在实际分布式光伏高渗透率的台区里是真实存在的尤其是节假日负荷低、光伏出力又最高的时段。电压分布曲线还有一个很有意思的现象DG接入后整条馈线的电压最低点不一定在末端而是可能转移到DG接入点与下游某个负荷集中区域之间。分析时不能只看末端电压要把全线路的电压曲线拉出来看哪个节点最低、哪个节点最高都要心里有数。4.3 DG接入对网络损耗的影响规律网损是分布式电源接入经济性最核心的指标。从仿真结果看DG容量从零开始增大时系统总网损先下降、后上升存在一个最优容量点。这个现象背后的物理机制不难理解初始阶段DG就地平衡负荷缩小了线路输送功率电流减小焦耳损耗下降但当DG出力超过本地负荷时多余功率被迫倒送到上级电网电流不降反升损耗又开始增加。不同接入位置的“最优容量”差别很大。DG接在馈线末端附近对网损的降低效果特别明显但最优容量窗口很窄容量一超出就会迅速恶化DG接在靠近变电站的首端节点对网损的影响幅度较小基本是平稳变化。这个结论对规划特别有指导意义如果业主想在末端台区装光伏必须严格校核接入容量不能简单按“屋顶面积能装多少就装多少”来定。计算网损随容量变化的曲线时我建议每间隔50到100kW取一个点把曲线画平滑一点。尤其是寻找网损最低点时间隔太粗容易漏掉峰值区间影响结论精度。如果你的研究涉及经济性分析还可以把网损降低带来的收益量化成电费和DG投资做对比这就超出了纯潮流仿真的范畴但方向上很有价值。4.4 接入位置与功率因数的敏感性对比做敏感性分析时我会把“位置”和“容量”两个维度做成热点图或对比表这样规律一眼就能看出来。比如把DG分别接到节点8、节点18、节点33分别计算三个位置下的网损-容量曲线放在同一张图里对比可以直观看到末端节点DG的“双刃剑”效应最明显。功率因数的影响也不容忽视。在同样的有功出力下如果DG能发出容性无功功率因数超前对电压支撑的效果会进一步放大反之如果DG以滞后功率因数运行吸收无功电压提升效果会打折扣。这套Matlab代码里可以直接修改dg输入的第三个分量来模拟不同功率因数运行状态做方案对比时多跑几组会得到更全面的结论。5. 常见问题与排查技巧实录5.1 潮流不收敛怎么办前推回代法在辐射状配电网里一般收敛性很好但如果仿真中迭代发散八成是下面几个原因DG注入功率设置过大导致局部节点电压严重越限数学上无法收敛支路首末端方向填反父子关系混乱负荷数据里出现孤立节点或回路。排查时我会先把DG容量设为0跑基准潮流如果能收敛就说明网络建模没问题问题出在DG参数上逐步降容量找到临界值如果基准场景就不收敛优先检查支路方向。还有一个小警告迭代初始值用平启动全部电压设为1.0通常没问题但如果配电网负荷太重平启动可能导致初期迭代振荡这时可以改用上一轮场景的潮流结果作为初始值特别是在连续扫描DG容量时能显著减少迭代次数。5.2 单位不统一导致结果离谱这是新手最容易踩的坑。支路阻抗用欧姆、基准容量用MVA、负荷用好几百kW混在一起计算数值会错到离谱。我的做法是统一按下面这套规则导线电阻电抗用欧姆负荷有功用kW、无功用kvar基准容量取10000kVA基准电压取12.66kV程序开头就计算基准阻抗后续全部转标幺值参与迭代只在最后输出时乘回基准值。代码里加一行注释写明“全程序统一使用标幺值仅在输出端恢复实际单位”能省掉大量排查时间。有时电压结果看起来正常但网损异常偏大我会怀疑是不是负荷数据里某些节点写错了或者支路阻抗多了一段。用前推回代结果与理论功率平衡对比一下最方便所有负荷功率之和加上全网有功损耗应该等于根节点注入的总有功功率误差大于0.1%就说明哪里错了。5.3 DG模型选择导致的偏差把分布式电源一律处理成PQ节点在稳态分析里是最稳妥的做法但也要知道它的局限。实际光伏逆变器在电网电压异常时有低电压穿越逻辑可能瞬间切换出力模式恒功率控制在大扰动场景下也不成立。如果你的项目需要评估暂态响应或者电压控制策略本文这套前推回代代码就不够了需要转Matlab/Simulink搭详细变流器模型。即便是稳态分析也要注意PQ节点模型的适用边界在DG容量远小于本地负荷时模型和实际基本一致但当DG导致局部电压越限时真实逆变器的无功限幅、降额运行等特性会让结果比模型保守。所以仿真结论里出现电压越限不要直接判定“实际一定会越限”只能说“如果不加控制措施存在过电压风险”。5.4 画图与数据导出的细节Matlab画电压分布图时建议横坐标从0号根节点开始这样能看到从变电站母线到线路末端的完整电压剖面。如果只想展示某个局部线路段可以切片处理但对比图里最好保留全段。多条曲线叠在同一个figure里注意不同方案用不同颜色和线型避免全是一色实线没法区分。数据导出方面我强烈推荐用表格保存每轮仿真指标而不要只存图。图会丢信息表格才是可以二次分析的基础。代码里建一个矩阵每一行对应一种方案列包括DG接入节点、容量、功率因数、电压最大值、电压最小值、系统总网损、网损率等7到8个字段最后用writetable导出成Excel写报告时直接从表里取数高效且不容易错。6. 几个提升研究深度的扩展方向代码跑通、基础结果拿到之后如果想让研究更上一个层次可以在当前基础上扩展几个方向。一是时序仿真把光伏出力曲线和负荷曲线按24小时分段接入观察电压和网损在一天内的波动这对评估高渗透率光伏台区的电压越限风险特别有用。二是随机性分析用蒙特卡洛方法对DG出力和负荷进行大量抽样统计电压越限概率这个方向在可靠性评估和规划风险分析里比较受认可。三是选址定容优化把本文的潮流计算封装成目标函数外层用遗传算法或粒子群算法搜索最优的DG接入位置和容量很多硕士论文就是在这套逻辑上展开的。我个人在实际操作中的体会是仿真程序的架构远比一次两次跑通重要。前推回代法本身不难难的是把网络数据、潮流核心、后处理分析分层拆开让代码能够灵活支持不同场景。这套Matlab实现把这三层都做了简化拆解你在自己项目里可以直接套用框架。最后再分享一个小技巧每次修改参数后顺手把关键结果和迭代次数打印到命令行不要只看最后的曲线迭代次数突然增多往往是运行状态逼近极限的信号提前发现能省很多排查时间。