配电网重构、网架规划、分布式电源接入优化，只要你在做这几个方向，一定被“辐射状拓扑约束”折磨过。需求听起来很简单：配电网正常运行必须是树状结构，不能有环，所有负荷还得连通到电源。可要把这句大白话翻译成优化模型里可求解的数学约束，很多人当场卡住——树约束天然是离散的、组合的，直接枚举生成树在工程规模下完全不可行。早些年我复现某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-1：
sum(z) == N - 1这一步直接对应工程上的操作：你要恰好打开 L 个开关。注意这里务必用等号，不能写成大于等于或小于等于。用小于等于可能出现断线不够、拓扑仍有环的情况；用大于等于则可能断线过多，形成多个孤岛。这个我在调试中吃过亏，后面会专门讲。
2.4 可选增强：基本回路约束
断线数量加连通性，在整数可行解层面已经严格保证树形。但混合整数规划的线性松弛往往很弱，求解器要花费大量节点做分支。一个工程上常用的增强方法是按图论中的基本回路加约束：对每个独立回路，要求回路内至少有一条断开支路，写成：
sum_{i in loop}(1 - z_i) >= 1这组约束在整数解上其实是冗余的，因为断线总数已经保证每个回路平均被断一次，且连通性会排除异常结构。但它对线性松弛的收紧效果明显。求解器在松弛解里最容易出现的非法结构就是“闭合边数满足要求但拓扑是一个大环”，这条约束能直接切掉这类不可行解，大幅减少分支数量。
2.5 与其他建模方法的横向对比
我整理了四种常见辐射状约束方案的对比，方便大家选型：
| 建模方案 | 约束数量 | 变量规模 | 松弛质量 | 实现难度 |
|---|---|---|---|---|
| 割集约束（动态生成） | 指数级，需割平面 | O(M) | 强 | 高 |
| 单商品虚拟流 | O(N+M) | O(M) | 中等偏弱 | 低 |
| 断线解环 + 断线总数 | O(N+M+L) | O(M) | 中 | 低 |
| MTZ 有向流约束 | O(M) + 逻辑约束 | O(N+M) | 弱 | 中 |
从我的实际使用感受来说，动态割集方案在小规模系统上表现很好，但代码复杂，需要自己实现求解器回调。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)}(end+1) = i; adj{branch(i,2)}(end+1) = 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(end+1) = 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{end+1} = 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 结果校验函数怎么写
不管求解器给出的最优解目标多漂亮，我都要跑一遍独立的校验程序，确认闭合子图是合法的辐射状拓扑。校验用最朴素的 DFS：
function [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}(end+1) = j; adj{j}(end+1) = 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这些细节在集成到大模型后突然爆发问题。我在每一个项目里都会保留一个独立的最小测试算例，任何一次模型改动都先跑小系统再上大系统。这个习惯帮我省下了大量排错时间。断线解环的方法论本身并不复杂，它最值钱的地方在于把“树的拓扑约束”变成了“三组线性约束”，而且物理意义还特别清晰。如果你正在做配电网方向，建议直接把它沉淀成自己的通用工具函数，以后无论遇到重构、规划还是分布式电源优化，都能用得上。
