1主要内容
程序主要以33节点为研究对象，编制配电网故障重构模型，采用虚拟潮流约束保证配电网的连通性和辐射性，以网损和负荷损失作为目标函数，包括潮流约束、电压电流约束、sop约束、辐射性约束等，该方法可实现单条线路故障或者多条线路故障下的重构。
该虚拟潮流方法在《主动配电网故障恢复与孤岛划分模型/多时段》里面有详细模型，如下：
辐射状及联通性约束
引自《A New Model for Resilient Distribution Systems by Microgrids Formation》
英文文献《A New Model for Resilient Distribution Systems by Microgrids Formation》给出了两个比较值得关注的信息：
网络辐射性的充要条件:【每个子图联通】【线路支路数量=节点数量-子图数量】
给出了虚拟潮流的计算方式。
这样就打通了二阶锥模型中配网辐射性和连通性约束实现问题，之前有文献采用图论方式来实现，但是该方法并非充要条件，导致在变换故障线路过程中还需要进一步增加其他约束，限制了方法应用的推广性，而虚拟潮流方式真正解决这个问题。
单条线路故障
虚拟潮流方式严格保证了网络的辐射性和连通性，实现最佳效果。
多条线路故障
设置了6/13/33三条线路故障，得到如下结果。
虚拟潮流方式严格保证了网络的辐射性和连通性，实现最佳效果。
2部分代码
%% 定义决策变量lamda=sdpvar(33,N,'full');x_Iij_square=sdpvar(37,N,'full'); x_ui_square=sdpvar(33,N,'full');x_pij=sdpvar(37,N,'full'); x_qij=sdpvar(37,N,'full');x_p_sop1=sdpvar(2,N,'full'); x_q_sop1=sdpvar(2,N,'full');opt_num=1;beta=binvar(37,37,'full');Zij=binvar(nl,1);%网架结构% Z0=[ones(nl-nc,1);zeros(nc,1)];%初始拓扑% assign(Zij,Z0);% x_sop_L=sdpvar(2,N,'full');%% 定义约束条件Constraints=[];%for opt_num=1:NConstraints = [Constraints, sum(Zij) == 32];Constraints = [Constraints, 0<= x_pij(33:37) <= 0.1.*Zij(33:37)];%% 定义目标函数f=r_ij'*x_Iij_square(:,opt_num)+100*sum(p_load)+100*sum(-lamda.*p_load);Constraints=[Constraints,x_ui_square(1,opt_num)==12.66^2];Constraints=[Constraints,0<=lamda<=1];
4下载链接
扫码或阅读原文即可跳转程序页面：
想了解更多程序？按照下面教程查看程序目录！
