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

本文主要是介绍Matlab|孤岛划分|弹性配网故障划分模型，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

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:NConstraints=[Constraints,x_ui_square(1,opt_num)==12.66^2];for k=2:33node_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)];%线路有功约束
%% 支路欧姆定律