您现在的位置是:首页 >其他 >含sop的配电网重构(含风光|可多时段拓展)网站首页其他

含sop的配电网重构(含风光|可多时段拓展)

电力程序小学童 2024-06-17 10:32:11
简介含sop的配电网重构(含风光|可多时段拓展)

目录

1 主要内容

2 部分程序

3 下载链接 


主要内容

之前分享了很多配电网重构的程序,每个程序针对场景限定性比较大,程序初学者修改起来难度较大,本次分享一个基础程序,针对含sop的配电网重构模型,含风电和光伏,优化了33节点网络电压合理性,任意修改故障线路不仅得到满足辐射和联络性网络重构结果,也能得到完美的节点电压图。

采用IEEE33节点系统,光伏接入7、13、27节点,风机接于5节点,但是由于考虑单时间段(1时刻),光伏出力为0,风机功率为0.8MW,单线路故障得到如下结果图:

10线路故障:

19线路故障:

多线路故障结果:

5、31线路故障:

部分程序

Branch=[Branch;xb];
r_ij=Branch(:,4);   x_ij=Branch(:,5);
% p_Solar=zeros(33,1);
p_Wind=zeros(33,1);
p_Wind(5,1)=0.8;%5节点加风电
p_Load=p_load;
% p_Solar(5)=0.05;p_Solar(11)=0.05;p_Solar(15)=0.05;p_Solar(22)=0.05;p_Solar(31)=0.05;%分布式功率
q_Solar=p_Solar.*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');
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);%网架结构
lx=[5,31];%故障线路编号
% Z0=[ones(nl-nc,1);zeros(nc,1)];%初始拓扑
% assign(Zij,Z0);
% x_sop_L=sdpvar(2,N,'full');
%% 定义约束条件
Constraints=[];
%for opt_num=1:N    
Constraints = [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];%负荷消纳比例
    Constraints=[Constraints,Zij(lx)==0];%设置故障线路
    Constraints=[Constraints,Zij(35)==0];
    Constraints=[Constraints,0<=x_Iij_square<=100*Zij];

3 下载链接 

点击直达! 

风语者!平时喜欢研究各种技术,目前在从事后端开发工作,热爱生活、热爱工作。