微❤关注"电气仔推送"获得资料（专享优惠）
参考文档：
A New Model for Resilient Distribution Systems by Microgrids Formation；
主动配电网故障恢复的重构与孤岛划分统一模型；
同时考虑孤岛与重构的配电网故障恢复运行策略。
运行环境：
matlab+yalmip+cplex
主要内容：
程序主要以33节点为研究对象，编制配电网故障重构模型，采用虚拟潮流约束保证配电网的连通性和辐射性，以网损和负荷损失作为目标函数，包括潮流约束、电压电流约束、辐射性约束等。
主动配电网故障恢复的重构与孤岛划分统一模型 随着分布式电源和储能装置大量接入配电网，配电网 在发生故障之后可以进行重构和孤岛划分，提高故障恢复的水平。建立包含多类型分布式电源、柔性负荷和储能的多 时间段故障动态恢复模型，并计及故障恢复时间和检修次序。利用二阶锥技术将 所建模型转换成混合整数二阶锥规划模型。结果表明，与传统方法相比，所提故障恢复方法计算速度更快，故障恢复率更高。
部分代码：
%% 检修策略约束
for i=1:nl
if ismember(i,F)
Constraints = Constraints, Zij(i,1)==0;%故障线路初始状态
for t=2:N
Constraints = Constraints, Zij(i,t)\>=Zij(i,t-1);%故障线路下一个时刻状态要大于等于上一个时刻
end
end
end
\~,na=size(F);%故障线路条数
for t=1+Tf:N
Constraints = Constraints, sum(Zij(F,t))-sum(Zij(F,t-Tf))\<=h;%每个时段最多只能检修h条故障线路
end
for i=1:6
for t=2:N
Rep(i,t-1)=Zij(F(i),t)-Zij(F(i),t-1);%线路维修状态和线路状态关系
end
%% 定义目标函数
f=sum(sum(repmat(r_ij,1,N).*x_Iij_square))+100*sum(sum(p_Load1))+100*sum(sum(-lamda.*p_Load1));
Constraints=Constraints,0\<=lamda\<=1;
%% 有功与无功平衡
for opt_num=1:N
Constraints=Constraints,x_ui_square(1,opt_num)==12.66\^2;
for k=2:33
Constraints = Constraints, -0.01.\*Zij(33:37,opt_num)\<= x_pij(33:37,opt_num) \<= 0.01.\*Zij(33:37,opt_num);%线路有功约束
Constraints = Constraints, -4.\*Zij(1:32,opt_num) \<= x_pij(1:32,opt_num) \<= 4.\*Zij(1:32,opt_num);%线路有功约束
%% 支路欧姆定律
%欧姆定律约束
m = (1.05*1.05 - 0.95*0.95)*12.66^2;
M = (ones(nl,N) - Zij)*m;
%% 运行电压电流约束
Constraints=Constraints,x_ui_square(:,opt_num)\>=(0.95\*12.66)\^2;
Constraints=Constraints,x*ui*square(:,opt_num)\<=(1.05\*12.66)\^2;
运行结果：
