配电网重构与辐射状拓扑约束：断线解环建模及Matlab实现

配电网重构、网架规划、分布式电源接入优化只要你在做这几个方向一定被“辐射状拓扑约束”折磨过。需求听起来很简单配电网正常运行必须是树状结构不能有环所有负荷还得连通到电源。可要把这句大白话翻译成优化模型里可求解的数学约束很多人当场卡住——树约束天然是离散的、组合的直接枚举生成树在工程规模下完全不可行。早些年我复现某EI期刊里的配电网重构方法时就被一个反直觉的退化拓扑坑过明明加了“闭合支路数等于节点数减一”的经典约束求解器却给出一堆孤岛加一个环的“非法树”。后来换成基于断线解环思想的建模方式整个问题才清爽起来。这篇文章把该方法的建模原理、线性化技巧、Matlab实现和调试经验完整梳理一遍给正在做配电网优化方向的朋友一个可以直接照抄的模板。1. 先搞清楚辐射状约束到底难在哪1.1 一个反直觉的翻车案例很多做配电网重构的人第一次写的约束是这样的闭合支路数等于节点数减一再加一个“每个节点只能有一个父节点”的有向化处理。理论上这像是树的定义实际求解时却经常跑出一个环加几个孤立点的拓扑。原因很简单只约束“边数 N-1”并不等价于树。一个环加上若干孤立点边数照样等于 N-1。而“单父节点”约束在没有强制根节点可达性时孤岛照样出现因为岛内每个节点也可以满足单父节点约束。这个坑我印象非常深。当时用的是某个商业求解器最优解的目标函数值看起来合理拓扑图一画出来却完全不能用于实际调度。后来仔细检查才发现问题不在求解器而是约束组本身不完整。树的定义是两件事同时成立连通且边数为节点数减一。缺了连通性边数约束就形同虚设。这看起来像基础图论但在优化模型里把“连通性”用线性约束完整表达出来才是真正麻烦的开始。1.2 树约束的数学本质从图论角度看一个 N 个节点、E 条支路的连通图是一棵树当且仅当 E N-1。注意是“连通图”加“边数条件”两个条件同时满足。经典的生成树约束通常写成割集形式对任意非空节点子集 S从 S 连向外部的闭合支路数量至少为 1。这类约束能完美刻画连通性加无环但问题在于所有非空子集 S 的数量是 2^N 级别不可能全部写进模型。常见的做法是让求解器在分支定界过程中动态生成约束即割平面。这种方案实现复杂度高对新手很不友好。断线解环走的完全是另一条路线先承认原始网络是闭环设计存在 L 个独立回路然后通过“每个回路至少断一条边”加上“断边总数等于回路数”再配合虚拟潮流保证连通性最后剩下的子图自动是一棵树。这个角度和配电网工程师的实际操作习惯完全吻合——本来就是闭环设计、开环运行你打开联络开关的过程就是解环。1.3 断线解环与工程直觉的契合配电网有一种与输电网完全不同的运行特性闭环设计开环运行。变电站到用户常常有多条备用通道正常运行时只保留一条树形通路。这意味着辐射状结构不是一个随意的数学假设而是物理运行方式的必然结果。断线解环模型复用的正是这个工程直觉给定整个闭环网络你需要在每个基本回路上“开一刀”打开的总开关数正好等于独立回路数同时保证所有负荷仍然带电。这套语言在跟现场工程师沟通时极其好用。那次我尝试用割集约束解释模型对方听得一头雾水换成“每个环必须断一刀总断刀数等于原有环数最后还要保证全网络带电”的说法对方立刻理解了。这也是很多EI论文倾向于用断线解环形式建模的原因之一——物理意义清晰便于工程人员理解和接受。2. 断线解环建模三组约束就够2.1 变量定义与网络描述设配电网有 N 个节点M 条支路节点 1 为根节点变电站母线或馈线首端。这里每条支路可以是分段开关、联络开关也可以是从未安装开关的固定线路统一用一个二进制变量 z_i 表示状态z_i 1 表示支路 i 闭合z_i 0 表示支路 i 断开。为了表达连通性再引入一组连续辅助变量 f_i表示支路 i 上的虚拟潮流。这组变量不是真实功率它唯一的作用是在抽象的“流”层面保证网络的连通性。真实的潮流约束后面会另加这里先把拓扑关系做干净。从图论角度看原始闭环网络如果是连通的那么它的独立回路数 L M - N 1。这个公式来自经典欧拉公式也可以在生成树视角下理解一棵树只需要 N-1 条边多余的 M-(N-1) 条边每增加一条就多出一个独立回路。2.2 连通性用虚拟潮流替代割集虚拟潮流法的核心思想非常朴素给每个非根节点发一个单位需求从根节点流出 N-1 个单位流量流量只能通过闭合支路传输。如果某个节点和根之间没有闭合路径它就无法收到这 1 个单位的需求约束方程必然矛盾。数学上写成矩阵形式很简洁。定义节点-支路关联矩阵 A行是节点列是支路。对每条支路人为规定一个正方向起点节点处 A 取 1终点节点处取 -1。那么对所有非根节点有A(2:N,:) * f ones(N-1,1)根节点则有A(1,:) * f -(N-1)再把支路状态耦合进去-M * z_i f_i M * z_i这里 M 直接取 N-1 就够因为整个虚拟流的总量不会超过 N-1。当 z_i 0 时f_i 被强制为 0虚拟流不能穿过断开支路z_i 1 时f_i 可以在上下限之间自由取值。为什么这组约束能证明网络连通我给一个直观证明假设闭合子图中存在一个不包含根节点的连通分量 S。S 内所有节点都要满足净流出等于 1那么对 S 内所有节点求和内部支路的正负贡献相互抵消边界支路成为唯一净输出通道。为了凑出至少 |S| 的净流出量S 整体必须向外部输出大量虚拟流。但外部没有对应的注入源所有流量最终只能来自根节点。如果 S 与根不连通这个净流出无法被满足方程组无解。因此所有非根节点必然处于包含根节点的同一个连通分量中。2.3 断线数量解环的核心等式原始闭环网络的独立回路数是 L M - N 1要把闭环变成树必须恰好断掉 L 条支路。写成约束就是sum(1 - z) L等价于闭合支路数等于 N-1sum(z) N - 1这一步直接对应工程上的操作你要恰好打开 L 个开关。注意这里务必用等号不能写成大于等于或小于等于。用小于等于可能出现断线不够、拓扑仍有环的情况用大于等于则可能断线过多形成多个孤岛。这个我在调试中吃过亏后面会专门讲。2.4 可选增强基本回路约束断线数量加连通性在整数可行解层面已经严格保证树形。但混合整数规划的线性松弛往往很弱求解器要花费大量节点做分支。一个工程上常用的增强方法是按图论中的基本回路加约束对每个独立回路要求回路内至少有一条断开支路写成sum_{i in loop}(1 - z_i) 1这组约束在整数解上其实是冗余的因为断线总数已经保证每个回路平均被断一次且连通性会排除异常结构。但它对线性松弛的收紧效果明显。求解器在松弛解里最容易出现的非法结构就是“闭合边数满足要求但拓扑是一个大环”这条约束能直接切掉这类不可行解大幅减少分支数量。2.5 与其他建模方法的横向对比我整理了四种常见辐射状约束方案的对比方便大家选型建模方案约束数量变量规模松弛质量实现难度割集约束动态生成指数级需割平面O(M)强高单商品虚拟流O(NM)O(M)中等偏弱低断线解环 断线总数O(NML)O(M)中低MTZ 有向流约束O(M) 逻辑约束O(NM)弱中从我的实际使用感受来说动态割集方案在小规模系统上表现很好但代码复杂需要自己实现求解器回调。MTZ 方法在某些路线规划问题里常见但在配电网重构中松弛质量差大系统上收敛慢。断线解环加虚拟流的组合在规模和实现难度之间取得了最好的平衡——约束条数线性增长不需要回调YALMIP 配合商业求解器直接就能跑是复现 EI 论文时最高效的选择。3. Matlab代码实现从基本回路到求解器对接3.1 整体框架设计我习惯把代码拆成四个模块。数据读入模块负责处理节点坐标、支路端点、电阻电抗、负荷等基础信息拓扑预处理模块负责构建关联矩阵、生成 BFS 生成树、提取基本回路建模模块用 YALMIP 把辐射状约束和潮流约束拼接成优化问题求解与校验模块调用求解器并对结果做连通性、无环性校验。模块之间用函数隔离调试时能快速定位问题。这里有一个重要的工程习惯建模之前先写好拓扑校验函数而不是等求解完再检查。因为校验函数不仅用于最终结果验证调试约束时也能派上大用场。我先手动构造几个已知的合法树和非法拓扑让校验函数跑一遍确认校验函数本身没毛病再拿它去验证求解结果。3.2 基本回路生成基本回路约束可以增强求解但前提是能正确找到独立回路。最简单的方案是用 BFS 生成一棵树每条非树边对应一个基本回路。代码如下function loops fundamental_loops(nb, branch) % nb: 节点数量 % branch: Mx2 矩阵每行是一条支路的两个端点 nr size(branch, 1); % 构建邻接表 adj cell(nb, 1); for i 1:nr adj{branch(i,1)}(end1) i; adj{branch(i,2)}(end1) i; end % BFS生成树 visited false(nb, 1); parent_edge zeros(nb, 1); queue 1; visited(1) true; tree_flag false(nr, 1); while ~isempty(queue) u queue(1); queue(1) []; for e adj{u} v branch(e, 1) branch(e, 2) - u; if ~visited(v) visited(v) true; parent_edge(v) u; tree_flag(e) true; queue(end1) v; end end end % 非树边对应基本回路 loops {}; for e 1:nr if ~tree_flag(e) u branch(e, 1); v branch(e, 2); % 回溯v到根的路径再回溯u到根的路径取对称差得到回路 path_v trace_path(parent_edge, v, nb); path_u trace_path(parent_edge, u, nb); loop_edges build_loop(branch, e, path_v, path_u); loops{end1} loop_edges; end end endtrace_path和build_loop是两个辅助函数前者利用父节点指针回溯到根记录路径经过的支路编号后者结合去找这两条路径的公共前缀把非树边加进去得到一个闭合回路的支路索引集合。实现细节有点繁琐但思路本身不复杂核心就是“树 一条非树边 一个基本回路”。这段代码有一个值得注意的点BFS 树的构造依赖于根节点可达所有节点。对此我默认原始闭环网络是连通的。如果你的网架本身有多个馈线组需要先把多个根节点合并成一个超级根节点或者分别对每个连通分量建模否则 BFS 生成树会漏掉非树边。3.3 辐射状约束的YALMIP实现核心建模代码不长我直接贴出可用的 YALMIP 片段z binvar(nr, 1); % 支路开关状态1闭合 f sdpvar(nr, 1); % 虚拟潮流 Constraints []; % 构造节点-支路关联矩阵 A zeros(nb, nr); for i 1:nr A(branch(i,1), i) 1; A(branch(i,2), i) -1; end % 连通性非根节点净流出1 Constraints [Constraints, A(2:end,:) * f ones(nb-1, 1)]; % 根节点净流入 N-1 Constraints [Constraints, A(1,:) * f -(nb-1)]; % 支路状态与虚拟流量耦合 M nb - 1; Constraints [Constraints, -M * z f M * z]; % 断线总数等于独立回路数 L nr - nb 1; Constraints [Constraints, sum(1 - z) L]; % 基本回路增强约束 for k 1:numel(loops) Constraints [Constraints, sum(1 - z(loops{k})) 1]; end这段代码可以直接嵌入任意配电网优化模型。实际使用时我会把辐射状约束封装成一个函数输入节点数、支路矩阵和回路列表返回约束集合和 z, f 变量。这样在做网架规划、配电网重构、分布式电源优化时一行调用就能复现同样的拓扑建模。3.4 与潮流约束联动的完整模型辐射状约束本身不包含电气量它只负责拓扑结构。实际优化问题中需要把虚拟流和真实潮流一起放进模型。配电网重构里最常用的是线性化 DistFlow 模型。支路 i 从节点 k 流向节点 j有功和无功功率分别为 P_i、Q_i节点电压幅值近似为 U根节点电压为 U0。线性化 DistFlow 的约束如下-U(j) U(k) (r_i * P_i x_i * Q_i) / U0 P_inj(k) sum(P_out(k)) load_P(k) Q_inj(k) sum(Q_out(k)) load_Q(k)如果支路 i 断开真实功率 P_i 和 Q_i 必须为 0这里同样用大 M 耦合P sdpvar(nr, 1); Q sdpvar(nr, 1); U sdpvar(nb, 1); Mp 2 * sum(load_P); % 有功功率上限取负荷总和的2倍 Mq 2 * sum(load_Q); Constraints [Constraints, -Mp*z P Mp*z]; Constraints [Constraints, -Mq*z Q Mq*z];大 M 取值不能拍脑袋。取太小会切掉可行解取太大会造成数值病态、求解器收敛慢甚至报错。我个人经验是有功功率上限取所有负荷有功总和的两倍无功同理。这样留给网络损耗和传输裕度足够的空间又不会让数值问题失控。虚拟流的大 M 直接取 N-1因为虚拟流总量有严格的上界不需要额外余量。4. 算例验证与性能实测4.1 5节点小型系统测试先用一个最小规模的系统验证约束正确性。5 个节点6 条支路根节点 1。原始网络闭环连接独立回路数 L 6 - 5 1 2。在这个系统上辐射状解表示要断开 2 条支路剩下的 4 条支路必须连通全部 5 个节点。我把约束单独跑一遍目标函数设为最小化开关操作次数。求解器给出的断开支路组合经校验函数验证闭合子图完全连通边数刚好等于 4无环满足树形要求。这套小型系统的意义在于快速暴露建模 bug。如果约束代写错误5 节点系统规模下很容易手工追踪验证。我在确认这个最小系统全绿之后才把同一套约束直接移植到更大系统大大减少了调试时间。4.2 33节点标准测试系统接下来是经典的 33 节点标准测试系统支路数 37独立回路数 L 37 - 33 1 5也就是需要打开 5 个联络开关。我同时测试了两种约束形式一组只含虚拟流加断线总数另一组额外加入基本回路增强约束。实测下来不加基本回路约束时求解器在分支定界过程中会反复探索环状松弛解中等负载水平下大约要探索上百个分支节点才收敛。加入基本回路约束后由于松弛解中的环状结构被提前剪枝分支数明显下降求解时间大概快 3 到 5 倍。这个加速效果在不同负载场景下都稳定出现所以我后来的所有模型都默认带基本回路约束。4.3 结果校验函数怎么写不管求解器给出的最优解目标多漂亮我都要跑一遍独立的校验程序确认闭合子图是合法的辐射状拓扑。校验用最朴素的 DFSfunction [conn_flag, is_tree] verify_topology(z_value, branch, nb) closed_edges find(z_value 0.5); % 闭合支路 % 检查闭合支路数 n_closed length(closed_edges); % DFS连通性 adj cell(nb, 1); for e closed_edges k branch(e,1); j branch(e,2); adj{k}(end1) j; adj{j}(end1) k; end visited false(nb, 1); dfs(1, adj, visited); conn_flag all(visited); is_tree conn_flag (n_closed nb - 1); end这函数只有二三十行却是我每一版模型都跑一遍的安全网。尤其是当你把辐射状约束和其他模型混合在一起时目标函数权重稍微一变求解器可能就走到了某个数值边缘结果违反约束却不自知。校验函数能第一时间发现这种问题。5. 常见问题与调试速查表5.1 无解或者根节点方向报错最容易踩的坑是关联矩阵符号不一致。假设支路 i 方向从节点 1 指向节点 2你却把节点 2 填成 1那虚拟潮流平衡方程就变成负的需求量模型大概率无解。调试方法单独检查关联矩阵每一行是否满足行和等于 0因为每条支路在一个节点是起点、另一个节点是终点所有节点行和必须为 0。如果不为 0说明方向填错了。另一个常见问题是将根节点约束方向写反。非根节点净流出 1根节点净流入 N-1二者的符号恰好相反。每次换根节点时都要重新确认符号定义。5.2 出现孤岛或者仍含环如果求解结果里闭合子图不连通或者有环优先检查断线总数约束是否写成了不等式。用sum(1-z) L的后果是求解器可能少开几个开关拓扑上保留环。用sum(1-z) L又可能多断开形成孤岛。这个约束必须严格等于。我在改模型时有过一次把等号误写成大于等于结果同一个系统跑出了两个孤岛后来靠校验函数才发现。5.3 求解太慢数值病态求解慢往往和大 M 取值有关。虚拟流的大 M 我固定用 N-1这是有理论依据的不要盲目放大。功率潮流的大 M 要与总负荷匹配取过大会让线性松弛失去意义分支定界效率暴跌。我见过有人把功率大 M 取到 1e6结果求解器在 33 节点系统上跑了半小时都不收敛改成总负荷的两倍后几十秒出解。5.4 多馈线系统怎么处理实际配电网可能有多个变电站电源点也就是多个根节点。如果直接把每个电源点都作为“根”虚拟潮流模型就崩溃了因为每个节点需求 1多个根节点同时注入会导致流量无法平衡。我采用的方法是把多个变电站母线合并成一个超级根节点所有到变电站母线的支路都连接到这个超级节点虚拟流模型照常使用。这个方法在网架规划中非常好用因为规划阶段本身就要考虑未来可能接入的电源位置。5.5 问题排查速查表现象可能原因排查方法模型无解关联矩阵方向错检查 A 每行求和是否为0模型无解根节点约束符号反检查 A(1,:)*f 是否为 -(N-1)结果有孤岛断线总数用了 改为严格等于结果仍有环断线总数用了 改为严格等于求解慢大M过大虚拟流取 N-1功率取总负荷2倍求解慢未加基本回路约束生成基本回路并追加增强约束BFS生成回路失败网络不连通先合并多个根为超级根节点最后再分享一个实际操作中的体会。复现这类建模方法时最大的风险往往不是理论看不懂而是符号、方向、大M这些细节在集成到大模型后突然爆发问题。我在每一个项目里都会保留一个独立的最小测试算例任何一次模型改动都先跑小系统再上大系统。这个习惯帮我省下了大量排错时间。断线解环的方法论本身并不复杂它最值钱的地方在于把“树的拓扑约束”变成了“三组线性约束”而且物理意义还特别清晰。如果你正在做配电网方向建议直接把它沉淀成自己的通用工具函数以后无论遇到重构、规划还是分布式电源优化都能用得上。