1 主要内容
程序主要复现《A New Model for Resilient Distribution Systems by Microgrids Formation》，建立灾害情况下配网优化孤岛划分方案，通过虚拟潮流的方式优化得到弹性配网网格划分策略，充分保证配网每个网格的联络性和辐射性，为配电网恢复重构提供重要支撑，程序以33节点系统为例进行实现，验证了方法的有效性。
1.1 DistFlow 模型
1.2 虚拟潮流
1.3 目标函数
以恢复负荷最多为目标函数。
2 部分代码
%% 常量定义 delta_T=1; N=1; nl=37;nc=5; xb=[33 12 22 0.1 0.1;34 8 21 0.1 0.1;35 15 9 0.1 0.1;36 25 29 0.1 0.1;37 18 33 0.1 0.1]; Branch=[Branch;xb]; r_ij=Branch(:,4); x_ij=Branch(:,5); %考虑8:00时间点孤岛划分 p_Solar1=p_Solar(:,8); p_Load1=p_Load(:,8); p_Wind=zeros(33,N); pg_st=[7,13,27];%分布式电源位置 q_Solar=p_Solar1.*0.05;%无功按照有功的0.05倍 %% 定义决策变量 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'); Zij=binvar(nl,N,'full');%网架结构 %% 定义约束条件 Constraints=[]; %% 定义目标函数 f=sum(sum(p_Load1))+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 node_out=find(Branch(:,2)==k); node_in=find(Branch(:,3)==k); [na,nb]=size(node_out); [nc,nd]=size(node_in); Constraints=[Constraints,sum(x_pij(node_in,opt_num)-... r_ij(node_in).*x_Iij_square(node_in,opt_num))+p_Solar1(k,opt_num)+... p_Wind(k,opt_num)-lamda(k,opt_num)*p_Load1(k,opt_num)==sum(x_pij(node_out,opt_num))]; Constraints=[Constraints,sum(x_qij(node_in,opt_num)-... x_ij(node_in).*x_Iij_square(node_in,opt_num))-... q_load(k)==sum(x_qij(node_out,opt_num))]; end 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)];%线路有功约束 %% 支路欧姆定律
3 程序结果
大家“来找茬”，看看能找到以7/13/27为DG的网格吗？
