【EI复现】基于断线解环思想的配电网辐射状拓扑约束建模方法（Matlab代码实现）

 💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码、数据

---

💥1 概述

文章来源：

近年来，随着凸优化理论的不断发展和应用，很多学者采用数学规划方法解决配电网中的优化问

题[1-5]。使用数学规划方法的优势在于，只要将问题建模为凸优化模型或某些特定形式的混合整数凸优化模型，则可利用成熟的商业优化求解器求解并得到全局最优解。在数学规划模型中，约束条件都应是解析表达式，因此对于需要辐射状拓扑的优化问题，需要建立辐射状拓扑约束的数学模型，以确保最终得到的网络为辐射状。文献[6]提出配电系统重构问题的凸优化模型，其中表示辐射状拓扑约束的生成树约束(spanning tree constraints，ST)被广泛应用于其他研究中。文献[7]提出含 ST 约束的交流最优潮流问题的混合整数线性规划模型。文献[8-9]均提出配电网故障恢复问题的混合整数二阶锥规划模型，且其中辐射状约束为 ST 约束。文献[10]利用ST 约束作为辐射状拓扑约束，建立配电网恢复操作次序确定问题的混合整数线性规划模型；文献[11]建立了含 ST 约束的配电网重构和无功优化协调优化问题的混合整数半定规划模型。上述文献均采用ST 约束表示配电网辐射状拓扑约束。然而在大量测试的基础上，笔者发现 ST 约束在某些场景下无法保证得到的拓扑为辐射状。此外文献[12]亦给出了 ST 约束可能会得到含环的非连通图的算例，但该文献并未详细分析 ST 约束无法确保辐射状的原因及相关条件。文献[13]中提到，由于优化问题中潮流约束的存在，使得 ST 约束在某些条件下成立，但未明确给出在何种条件下潮流约束能够使其成立。对于辐射状拓扑建模，除 ST 约束外，文献[14]指出了满足辐射状约束的充分必要条件，即需要同时满足连通性和节点–边的数目关系，并从虚拟潮流的角度提出了确保图的连通性的单商品流约束 (single-commodity flow constraints，SCF)，与节点–边的数目关系共同组成辐射状约束(文中记作 SCF0)，被应用于相关研究中[13,15]。然而该组约束中仅用 1 个等式描述图中各个边的关系，包含表征线路状态变量的相关信息较少，可能导致混合整数规划求解过程较慢。

📚2 运行结果

部分代码：

duration_time = [2000,5760,1000];           % 三种负荷水平各自的持续时间
wi = (ones(NB,1)*duration_time)/8760;       % 每个节点每种负荷水平的权重
pload_i = sum(mpc.bus(:,2:4).*wi,2)/SB;     % 节点 i 负荷有功功率需求
qload_i = pload_i*sqrt(1 - cos_phi^2);      % 节点 i 负荷无功功率需求
Vmax = 1.1^2;                               % 节点电压最大值的平方
Vmin = 0.9^2;                               % 节点电压最小值的平方
Lij_max = 10;                               % 线路 ij 的电流幅值上限
Srate_i = mpc.substation(:,3)/SB/cos_phi;   % 与变电站 i 相连的变压器的额定功率
load('path_cycle_24.mat')                   % 读取24节点配电网中所有环以及根节点之间的所有路径
branch_to_node = zeros(NB,NL);              % 节点的上游支路
branch_from_node = zeros(NB,NL);            % 节点的下游支路
for k = 1:NL
    branch_to_node(mpc.branch(k,3),k) = 1;
    branch_from_node(mpc.branch(k,2),k) = 1;
end
 
%% 优化变量
xij = binvar(NL,1);                         % 0-1 决策变量，表示线路ij是否建设
yi = binvar(NS,1);                          % 0-1 决策变量，表示变电站i是否建设
Lij = sdpvar(NL,1);                         % 线路ij电流的平方
ps_i = sdpvar(NB,1);                        % 变电站节点 i 的从上级电网得到有功功率
qs_i = sdpvar(NB,1);                        % 变电站节点 i 的从上级电网得到无功功率
vi = sdpvar(NB,1);                          % 节点i电压的平方
Pij = sdpvar(NL,1);                         % 线路ij的有功功率
Qij = sdpvar(NL,1);                         % 线路ij的无功功率

🎉3 参考文献

文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。

[1]王颖,许寅,和敬涵,等.基于断线解环思想的配电网辐射状拓扑约束建模方法[J].中国电机工程学报,2021,41(07):2395-2404.DOI:10.13334/j.0258-8013.pcsee.200071.

🌈4 Matlab代码、数据