这两年分布式电源接入的消息从没断过分布式光伏装机量涨得有多猛做过配电网规划的人应该都深有体会。承载力评估这四个字如今已经成了配电网规划里绕不开的硬需求。一到并网季受理人员一摞一摞地收光伏报装单营业厅天天有人问我家屋顶还能不能装厂房屋顶的容量够不够。这种时候最怕的就是拍脑袋给答复——说少了电网资产在闲置新能源指标完不成说多了末端电压飙起来、变压器烧掉追责的时候一个都跑不掉。分布式电源接入配电网承载力评估要解决的正是把这个拍脑袋的过程变成可量化的计算一条10kV馈线在不越电压限值、不超设备载流能力的前提下到底还能挂多少千瓦的分布式电源。这篇文章我来完整复现一套基于核心论文思路的承载力评估方法并用Matlab把它从数学模型写到能跑出结果的代码同时附上算例参数和运行结果。内容适合正在做配电网规划、新能源并网研究的研究生也适合一线配网规划和运检的工程师参考。先放一个我自己的结论承载力从来不是一个固定的数而是一组边界条件下的极限值。边界条件变了结果可能差出一倍还多。下文我分六步讲先说明评估到底在算什么然后给出数学模型和约束再讲Matlab代码怎么组织、核心函数怎么落地接着用算例把结果跑出来最后交代复现过程中的坑和扩展思路。1. 承载力评估到底在算什么从一张光伏报装单说起1.1 承载力的本质一组约束下的极限值理解承载力最快的方式是用一个生活类比。一个房间能站多少人从来不是固定数字它取决于门有多宽、地板承重、有没有应急出口。配电网承载分布式电源也是同理一条馈线能接入多少容量取决于节点电压允许多大偏差、线路能过多少电流、变压器还有多少裕度以及接入之后会不会出现潮流反送重载。把这些都量化就是承载力评估在做的事。承载力更严谨的定义可以表述成在系统安全稳定运行的前提下在某个或某几个候选节点接入分布式电源当任一安全约束达到边界时所对应的最大接入容量。注意任一这两个字它意味着我们必须把电压、电流、容量这些约束放在同一个框架里同时判断而不是只看某个单一指标。有些报告只算电压不校核载流就会在导线特别细的老旧线路上栽跟头。1.2 三种算承载力的典型思路翻文献的时候会发现承载力评估大致有三大流派流派基本做法优点缺点最大容量搜索法逐步增加DG容量反复潮流计算直到越限物理意义清晰结果可直接用于并网批复计算量大需要可靠的潮流程序渗透率指标法以DG总容量与馈线最大负荷的比例直接给结论计算简单便于快速筛查过于笼统不考虑接入位置和网架差异多指标综合评价法把电压偏移、网损、载流等多个指标加权打分能给出好/中/差等级适合宏观对比权重主观性强边界不够明确我复现的这类核心论文绝大多数走最大容量搜索路线同时会结合综合评价对结果做分级。落到工程上就是潮流计算约束校验容量搜索这套组合拳也是供电公司做并网审查时最认可的技术路线。1.3 复现思路从论文框架到可运行代码论文里的方法一般写得很抽象比如建立目标函数设置约束条件采用改进粒子群算法求解读者想要直接照着写代码往往无从下手。我在复现时做了一个重要取舍把优化求解简化为指定接入位置用二分法搜索最大容量。原因是论文的核心贡献更多在于承载力评估的约束体系和指标建模优化算法换一种也能得到接近的结果但评估模型如果模糊输出就没有参考价值。有一点得提前说明本文属于方法级复现并非对着某篇论文的原始算例数据原样照搬。严格复现一篇论文通常需要原始网架和实测数据公开渠道拿不到。我能做的是用相同的评估思路在Matlab里把潮流计算、约束校验、容量搜索这条链路完整打通。你拿到这套代码逻辑换成自己手里的网架台账就能直接做工程评估。2. 数学模型与约束条件电压、载流和反向重载必须同时看2.1 目标函数与决策变量承载力评估作为一个优化问题目标函数很直白max P_DG也就是最大化接入的分布式电源总容量。决策变量有两类一是在哪些节点接入二是每个节点接多少。工程上更习惯先定接入位置再算允许容量所以我的代码把问题拆成两步指定接入节点搜索最大容量然后更换节点重复搜索最后对比不同节点的承载能力。这个处理方式比一次性求解多目标优化更容易落地也方便出报告。如果考虑多个候选节点同时接入目标函数就变成多个节点的容量之和。这时搜索空间变大二分法这种单变量搜索就不够用了我在代码里做了一个退化方案按比例逐轮增加各节点容量每轮都跑潮流校核直到出现越限。默认情况下代码走单节点二分多节点的逻辑也保留着方便扩展。2.2 等式约束节点潮流方程承载力评估的第一类约束是潮流方程。配电网的潮流方程是非线性方程组但辐射状配电网有个天然优势——功率流向是确定的可以用前推回代法求解不需要像输电网那样做大规模矩阵求逆。前推回代法分两步。回代时从馈线末端开始向电源侧逐段累加各支路下游节点的注入功率得到每条支路的功率或电流。前推时从电源侧开始用支路电流和阻抗计算电压降逐段修正下游节点电压。两者反复交替直到两次迭代的电压幅值最大偏差小于收敛阈值。这里先强调一个关键点潮流计算的精度直接决定承载力评估的可信度。收敛阈值如果取太松比如1e-3在电压接近上限时可能得出错误的还没越限结论导致承载力结果偏大。我的代码里默认用1e-6虽然多跑几轮迭代但结果更可靠。2.3 不等式约束工程限值不是拍脑袋不等式约束是承载力评估的核心也是论文和工程实践差异最大的地方。论文里常见0.95~1.05这种理想化电压范围但实际工程必须按国家标准来。约束项典型限值依据/说明10kV节点电压偏差0.93~1.07 puGB/T 12325-200810kV及以下三相供电允许±7%380V节点电压偏差0.93~1.07 pu同上三相供电220V节点电压偏差0.90~1.07 pu单相供电下限-10%支路载流率不超过100%规划建议80%取决于导线型号和敷设方式配变负载率短期不超过100%长期建议80%留裕度考虑N-1和负荷增长反向重载变电站母线倒送不超过主变容量防止潮流反送越过主变谐波畸变率电压THD不超过5%GB/T 14549-1993在本文的算例中起决定性作用的通常是10kV电压上限其次是线路载流率和反向重载。负荷越轻、线路越长的馈线约束越容易先卡在电压上这和现场经验完全一致。线路载流量和导线型号的关系也值得一提。我常用的参考值LGJ-70约265ALGJ-120约380ALGJ-185约515ALGJ-240约610A。实际工程计算时要根据敷设方式、环境温度做修正博文代码里先用典型值的标幺值作为限值。2.4 分布式电源的稳态建模PQ模型为什么够用分布式电源的稳态等值模型有PQ、PV、PI、恒阻抗几类。承载力评估场景下绝大多数情况选PQ模型就够原因有三。第一分布式光伏逆变器在并网点通常按指定功率因数运行从系统侧看就是一个负的PQ负荷第二PV模型需要模拟无功调节在辐射状配电网里会引入额外的迭代复杂度第三承载力评估关注的是接入容量上限在接近上限的运行点PQ模型给出的结果通常偏保守工程上宁可保守也不冒险。PQ模型的写法很直接某个节点接入三相容量S_DG、功率因数cosφ时有功P_DG S_DG·cosφ无功Q_DG S_DG·sinφ潮流计算中设为负的注入功率。功率因数的影响我在算例里专门测过0.95和1.0的功率因数相比承载力有明显差别。无功本地支撑能力被高估或低估都会直接反映在电压约束上。3. Matlab代码框架数据组织、拓扑分层与模块划分3.1 数据组织矩阵就是你的网架台账写Matlab代码我习惯先把数据格式定死。配电网计算程序最怕节点编号乱、量纲不统一。本文代码采用三套输入矩阵对应线参—负荷—DG三层台账line矩阵每一行是一条支路格式为[首端节点, 末端节点, 电阻R(Ω), 电抗X(Ω), 最大载流(A)]。首端节点约定为靠近电源侧那一端。bus矩阵每一行是一个负荷节点格式为[节点编号, 有功负荷P(kW), 无功负荷Q(kvar)]。变电站10kV母线单独作为平衡节点不进bus矩阵。dg矩阵每一行是一个分布式电源格式为[接入节点编号, 容量P_DG(kW), 功率因数cosφ]。量纲统一非常重要。有人在代码里全程用标幺值我在实现时采用有名值加标幺值混合节点电压和功率在计算时归一到标幺值线路阻抗在初始化时统一除以阻抗基准Z_B。这样调试时看一眼电压就知道对不对不用反复换算。基准值设定电压基准U_B10kV功率基准S_B10MVA那么阻抗基准Z_BU_B²/S_B10Ω电流基准I_BS_B/(√3·U_B)≈577A。前推回代公式里常用的电压降落近似形式是U_末端 U_首端 - (P·R Q·X) / U_首端实际代码按复数精确计算不省略虚部和交叉项。3.2 拓扑分层前推回代法能不能用一半看这里前推回代法能高效求解辐射状配电网前提是知道每条支路的上下游关系。我的做法是在主程序里做一个广度优先遍历BFS从电源节点开始给所有节点分层得到每个节点的父节点和子节点信息。回代顺序自然是从最后一层到第一层前推顺序是从第一层到最后一层代码写起来非常干净。BFS本身不复杂用邻接表存拓扑用队列从电源节点逐层扩散。真正要小心的是数据别建错支路首末端写反、或者拓扑里出现环潮流计算都会发散。实际配电网往往是闭环设计、开环运行承载力评估只针对当前联络开关断开后的运行方式必须删掉环路再算。3.3 模块划分与主流程整个Matlab程序我拆成六个模块职责分得很清楚模块函数名作用数据读取load_case.m读入line/bus/dg矩阵换算标幺值拓扑处理build_topology.mBFS分层生成父子关系潮流计算powerflow.m前推回代法输出节点电压U和支路电流I约束校验check_limits.m判断电压、载流等是否越限容量搜索search_capacity.m二分法搜索最大可接入容量结果绘图plot_result.m画电压分布曲线和承载力对比图主流程的伪代码如下load_case读入数据换算标幺值build_topology建立节点层级对每个候选DG节点 a. 调用search_capacity做二分搜索 b. 内部反复调用powerflow和check_limits汇总各节点承载力plot_result可视化这样拆模块还有一个额外的好处以后想换更复杂的潮流算法比如牛顿法只需要替换powerflow.m想加谐波约束也只在check_limits里加一个判断分支其他模块都不用动。4. 核心代码实现潮流、越限判断与容量搜索逐段拆解4.1 前推回代潮流子函数下面给出powerflow.m的核心代码。为了控制篇幅我抽掉了一些安全性检查但主体逻辑是完整的。需要提醒的是这段代码假设line矩阵已经按BFS层级排好序支路首端节点的层级从大到小排列这样回代时才能保证每条下游支路的电流已经先算完。实际工程中我会在build_topology里做这个排序避免手排出错。function [U, I_branch, iter] powerflow(line_pu, parent, child, bus_pu, dg_pu, U0) % 前推回代法潮流计算 % line_pu: 标幺值支路矩阵 [首端节点, 末端节点, R, X, Imax] % parent: 节点父节点编号 % child: 节点的子节点列表cell数组 % bus_pu: 负荷标幺值矩阵 [节点, P, Q] % dg_pu: DG注入标幺值矩阵 [节点, P, Q] nb length(parent); U U0; alpha 0.6; % 阻尼因子防止迭代发散 eps 1e-6; % 收敛阈值 max_iter 100; iter 0; while iter max_iter U_old U; S_bus zeros(nb, 1); % 各节点净注入功率 % 负荷正方向为从电网吸收 for i 1:nb S_bus(i) bus_pu(i, 2) 1j * bus_pu(i, 3); end % DG负的注入功率 for k 1:size(dg_pu, 1) S_bus(dg_pu(k, 1)) S_bus(dg_pu(k, 1)) - dg_pu(k, 2) - 1j * dg_pu(k, 3); end % 回代从末端到根节点求支路电流 nline size(line_pu, 1); I_branch zeros(nline, 1); for k nline:-1:1 j line_pu(k, 2); % 末端节点 % 该支路下游总功率 节点j注入 所有以j为首端的支路功率 S_down S_bus(j); for m 1:nline if line_pu(m, 1) j S_down S_down conj(U(j)) * I_branch(m); end end I_branch(k) conj(S_down / U(j)); end % 前推从根节点到末端更新节点电压 for k 1:nline i line_pu(k, 1); j line_pu(k, 2); U(j) U(i) - I_branch(k) * (line_pu(k, 3) 1j * line_pu(k, 4)); end % 阻尼修正 U alpha * U (1 - alpha) * U_old; iter iter 1; if max(abs(U - U_old)) eps break; end end end这段代码里的阻尼因子alpha是我复现时踩过坑之后加进去的。不加阻尼时一旦DG容量比较大、局部电压被抬得偏高迭代特别容易来回震荡甚至发散。加了0.6的阻尼后稳定性明显改善。不同系统可能需要微调基本原则是电压波动剧烈就调小收敛太慢就调大。4.2 越限判断子函数潮流算完下一步回答这个运行状态合不合格。我写了一个check_limits函数把电压约束和载流约束放在一起判断function [violation, reason] check_limits(U, I_branch, line_pu) U_min 0.93; % 10kV配电网电压下限 pu U_max 1.07; % 10kV配电网电压上限 pu violation 0; reason ; % 电压越限检查 for i 2:length(U) if abs(U(i)) U_min violation 1; reason sprintf(节点%d电压过低: %.4f pu, i, abs(U(i))); return; end if abs(U(i)) U_max violation 1; reason sprintf(节点%d电压过高: %.4f pu, i, abs(U(i))); return; end end % 支路载流越限检查Imax已在load_case中换算为标幺值 for k 1:size(line_pu, 1) if abs(I_branch(k)) line_pu(k, 5) violation 1; reason sprintf(支路%d-%d过载: %.2f A, ... line_pu(k, 1), line_pu(k, 2), abs(I_branch(k)) * 577); return; end end end工程上关心电压幅值所以电压判断用abs(U(i))取复数模值。我把I_branch乘回577A的电流基准值是为了让报错信息更直观不用再看标幺值估电流。实际运行时这两个判断顺序可以调换但如果电压越限和载流越限同时发生先报电压更符合现场经验——电压问题往往更隐蔽。4.3 容量搜索主循环二分法为什么比定步长强承载力搜索最朴素的做法是从0开始每次加100kW直到越限。代码简单但问题很大如果承载力在3000kW量级要跑30次潮流步长细化到10kW要跑300次。大网架的每次潮流迭代都要好多轮总耗时完全不可接受。所以search_capacity.m采用二分法把搜索次数从O(N)降到O(log N)。给定上限比如5000kW精度要求10kW最多只需要约10次潮流计算就能锁定结果。function P_max search_capacity(line_pu, parent, child, bus_pu, node_dg, pf, P_high) P_low 0; eps_p 10; % 容量搜索精度 10kW while (P_high - P_low) eps_p P_mid (P_low P_high) / 2; % 构造dg矩阵功率因数转为无功 dg_pu [node_dg, P_mid / 10000, P_mid / 10000 * tan(acos(pf))]; U0 ones(length(parent), 1); [U, I_branch, ~] powerflow(line_pu, parent, child, bus_pu, dg_pu, U0); [violation, ~] check_limits(U, I_branch, line_pu); if violation P_high P_mid; else P_low P_mid; end end P_max P_low; end二分法的关键是P_high必须取够。如果上限比真实承载力小最终结果会直接顶在上限上承载力偏小、失去参考意义。我通常用接入点馈线对应的主变容量和线路热稳极限的较小值作为P_high。代码里我还会在初始化时先跑一次P_high对应的潮流如果一开始就越限说明上限取小了程序直接报错提醒而不是默默输出一个错误结果。4.4 结果可视化算完之后我一般画两类图一类是不同接入容量下的馈线电压分布曲线一类是多个候选节点的承载力对比柱状图。前者用plot就能画figure; plot(1:nb, abs(U_no_dg), o-, LineWidth, 1.5); hold on; plot(1:nb, abs(U_500kw), s-, LineWidth, 1.5); plot(1:nb, abs(U_max), d-, LineWidth, 1.5); yline(1.07, r--, 上限1.07pu); yline(0.93, r--, 下限0.93pu); grid on; ylabel(节点电压 (pu)); xlabel(节点编号); legend(无DG, 接入500kW, 临界接入容量);从这三条曲线能一眼看出DG接入对电压分布的抬升作用也能看出哪个节点最先触及上限。这种图放进承载力评估报告里规划人员一眼就能看懂比一堆数字表格更有说服力。5. 算例实测10kV馈线接入分布式光伏的承载力测试5.1 算例系统与参数为了把原理讲透我构造了一条简化的10kV辐射状链式馈线一共6个节点。节点1是变电站10kV母线作为平衡节点电压取1.05 pu节点2到节点6是沿线负荷和光伏接入点。拓扑和参数如下支路首端末端导线型号长度(km)R(Ω)X(Ω)载流上限(A)112LGJ-2402.00.2640.710610223LGJ-1851.50.2550.503515334LGJ-1201.00.2700.320380445LGJ-1201.20.3240.384380556LGJ-700.80.3600.248265负荷数据节点2为300kWj150kvar节点3为500kWj250kvar节点4为400kWj200kvar节点5为300kWj150kvar节点6为200kWj100kvar。总负荷1.7MWj0.85Mvar对应的馈线最大负荷率约65%留了明显的DG接入空间比较接近当前农村和城乡结合部10kV馈线的实际状态。5.2 无DG的基准潮流先跑一次不接DG的基准潮流。平衡节点电压取1.05 puMatlab运行后各节点电压幅值如下节点123456电压(pu)1.05001.04161.02681.01230.99810.9805电压偏差(%)5.004.162.681.23-0.19-1.95末端节点6电压降到0.9805 pu虽然还在合格范围内但已经能看出这条5公里出头的线路压降不容忽视。在这个基础上叠加光伏出力电压会被进一步推高这正是承载力分析要盯住的核心对象。5.3 末端集中接入的二分搜索现在把光伏全部接在末端节点6功率因数取0.95用search_capacity计算承载力。P_high取5000kW搜索精度10kW实际运行中二分法只迭代了10轮左右。关键迭代点如下迭代试探容量(kW)末端电压(pu)是否越限125001.0612否237501.0741是电压上限331251.0680是电压上限428121.0645否529691.0664否630471.0672是电压上限............最终约30301.0670否临界最终末端节点6的承载力约为3030kW。注意此时卡住结果的是末端电压达到1.07 pu上限不是线路载流也不是变压器容量。换句话说这条馈线在不改造网架的前提下末端单点接入的光伏容量最好控制在线路本身热稳定能力的60%左右。5.4 接入位置与分散接入的影响把同样的方法换到其他节点结果差异非常大接入节点承载力(kW)主要限制因素节点2近电源侧约5300首段线路载流节点4中段约4600中段电压加载流节点6末端约3030末端电压上限节点246分散接入约6850多约束综合近电源侧的节点2接入到5300kW都还能接受因为电压抬升作用离末端远约束主要卡在首段LGJ-240线路的载流上。越靠近末端电压约束出现得越早承载力直线下降。最值得关注的是分散接入场景把总计约6850kW的容量按比例分到三个节点比集中接在末端整整多出一倍多。这个结论在现场非常有指导意义。一个村或一个园区的户用光伏与其都挤在馈线末端不如规划阶段就引导分散接入。承载力评估的价值就在这——它回答的是在哪里接、能接多少两个维度的问题而不只是简单给个容量上限。6. 复现过程中踩过的坑与扩展方向6.1 前推回代法不收敛初值、阻尼和PV节点复现过程中我调试最久的问题就是潮流不收敛。典型症状是无DG时一切正常DG容量一旦加到临界值附近迭代次数飙升电压来回震荡甚至直接算成NaN。排查下来主要有三个原因。第一个是初值取得太离谱。前推回代法的初值一般设所有节点电压为1.0 pu无DG或小容量DG时没问题。DG容量大了末端电压可能被抬到1.06以上初值还用1.0迭代初期电流计算偏差很大。我的解决办法是热启动把上一个容量点的潮流结果作为下一个容量点的初值能显著减少迭代次数。第二个是缺少阻尼。我在powerflow.m里加了alpha0.6的阻尼因子就是专门处理震荡的。实际调试时可以动态调整电压波动剧烈就调小到0.4收敛太慢就调大到0.8。第三个是我曾试图把分布式电源直接处理成PV节点。前推回代法对PV节点的无功迭代修正特别容易写错一旦写错就是发散。承载力评估第一版代码建议老老实实用PQ模型跑通之后再考虑扩展PV。6.2 搜索精度与计算速度的权衡二分法的搜索精度取10kW、上限取5000kW时单节点承载力计算只需要约10次潮流计算每次潮流约20轮迭代在普通笔记本上是毫秒级的事。但如果把场景扩大到IEEE 33节点甚至上百节点的实际馈线并且要扫描几十个候选节点总耗时就会上升到分钟级。两个优化很实用。第一缩小区间二次搜索先用100kW的粗精度快速定位承载力大致范围再用10kW精读在小区间内精搜总计算量能减少三分之一左右。第二并行扫描Matlab的parfor可以直接并行计算多个节点的承载力我在33节点算例上试过开8核并行后整体耗时压到原来的六分之一。还有一个容易被忽略的性能点支路顺序影响收敛速度。BFS分层之后支路矩阵最好按首端节点所在层数从小到大排序前推过程会更快收敛。我在代码里专门写了排序逻辑避免手排出错。6.3 静态评估的局限与时序评估扩展本文全套代码做的是静态承载力评估也就是取一个典型运行断面校核约束。这对工程初步筛查足够但有风险光伏出力和负荷都在实时变化静态断面选得偏保守或偏乐观结果都会偏离实际。更严谨的做法是时序评估取全年8760小时的光伏出力曲线和负荷曲线逐小时做潮流计算统计电压越限小时数和最大越限量用越限风险不超过某个阈值作为约束来求承载力。时序评估在代码上的改动其实不大只需要把powerflow和check_limits塞进一个按小时循环里再把结果汇总成越限概率指标。我复现时做了一个简化版本耗时大约是单断面静态评估的几百倍但输出可以支撑一份完整的并网影响分析报告。再往前走一步可以考虑光伏出力和负荷的不确定性用蒙特卡洛抽样生成大量场景再配合场景削减提炼出几十个代表性场景做机会约束校验。这样得出的承载力不是单点值而是带置信水平的区间规划决策的安全度和灵活度都会更高。最后聊一点自己的体会。这套方法论刚接触时会觉得不就是一个潮流加一个搜索吗实际上手之后会发现真正决定承载力评估质量的是数据质量和边界条件的还原度。同样的馈线负荷曲线取法不同、电压基准不同、导线载流量修正系数不同结果都能差出一截。所以用代码跑出数字只是第一步更重要的是拿到结果后要问自己一句这个数字对应的运行场景现实里真的可能出现吗如果答案是肯定的承载力才有工程价值如果只是理论边界就必须在报告里把边界条件写清楚让使用结果的人知道它的适用范围。承载力评估不是一次性工作。负荷在增长、网架在改造、分布式电源出力在变化承载力也会动态变化。我现在的做法是把它做成季度复算加年度精算的滚动机制代码固定、数据定期更新每次都能快速输出一份当前网架的承载力情况。这样比临时抱佛脚式的单次评估实用得多也更容易让数据沉淀成网架改造的决策依据。
