孤岛划分|弹性配网故障划分模型

目录

1 主要内容

1.1 DistFlow 模型

1.2 虚拟潮流

1.3 目标函数

2 部分代码

3 程序结果

4 程序链接

---

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))];     endConstraints = [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的网格吗？