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: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 程序结果

4 下载链接