目录

1 主要内容

2 部分代码

3 程序结果

4 下载链接

---

1 主要内容

之前分享了含sop的配电网优化模型，链接含sop的配电网优化,很多同学在咨询如何增加储能约束，并进行多时段的优化，本次拓展该部分功能，在原代码的基础上增加储能模型，以IEEE33节点为研究对象，在8节点增加储能，并且在7、10、24、27节点配置光伏，在13和30节点配置风电，储能不仅满足功率约束，同时需要满足soc约束条件。

2 部分代码

clear;clc;
run ieee_33_node_system.m;
run DG_Load.m
t1=clock;
%% 常量定义
delta_T=1;  N=24;%时间
S_sopi1=0.5;    S_sopj1=0.5;
S_sopi2=0.5;    S_sopj2=0.5;    %SOP容量
A_sop=0.02;     %SOP损耗
q_CB=60/1000;   C_loss=0.08;%与有功功率损耗相关的成本系数
WL = 0.833;     Wv = 0.167;
k12_0=1;
​
r_ij=Branch(:,4);   x_ij=Branch(:,5);  %%第四，五列
% p_Solar=zeros(33,24);p_Wind=zeros(33,24);p_Load=zeros(33,24);
% q_Solar=zeros(33,24);q_Wind=zeros(33,24);q_Load=zeros(33,24);
result=zeros(1,24);
Horizon=N;
%储能参数
capmax=400;
EESmax=100;
EESmin=0;
socmax=0.9;
socmin=0.2;
theta=0.01;%自放电率
yita=0.95;
%% 定义决策变量p_Solar
x_Iij_square=sdpvar(32,N,'full');   x_ui_square=sdpvar(33,N,'full');%线路电流 节点电压
x_pij=sdpvar(32,N,'full');      x_qij=sdpvar(32,N,'full');%线路传输有功   无功
x_p_sop1=sdpvar(2,N,'full');    x_q_sop1=sdpvar(2,N,'full');%sop1有功 无功
x_p_sop2=sdpvar(2,N,'full');    x_q_sop2=sdpvar(2,N,'full');%sop2有功 无功
x_p_sop1_loss=sdpvar(2,N,'full');   x_p_sop2_loss=sdpvar(2,N,'full');%sop1损耗 sop2损耗  
K12_t=intvar(1,N,'full');     k12_t=sdpvar(1,N,'full');%OLTC抽头选择 k12_t没用到
N_CB=intvar(1,N,'full');      Q_CB=sdpvar(1,N,'full');%CB数量 CB无功
x_AUX=sdpvar(33,N,'full');%电压偏差
x_b_OLTC=binvar(11,N,'full');  x_K121_t=intvar(1,N,'full');%抽头选择变量（只有一个是1） 
x_K122_t=intvar(1,N,'full'); %OLTC正负变化辅助变量  
x_N1_CB=intvar(1,N,'full');%CB变化辅助变量
K_t=intvar(1,N,'full');     N_t=intvar(1,N,'full');%变量没有用到
x_N2_CB=intvar(1,N,'full');   x_vc_k=sdpvar(11,N,'full');%CB变化辅助变量 oltc电压方值
x_ui_square_OLTC=sdpvar(1,N,'full');
​
x_P_ch = sdpvar(1, Horizon);%充电
x_P_dis = sdpvar(1, Horizon);%放电
x_u_ch = binvar(1, Horizon);%充电状态
x_u_dis = binvar(1, Horizon);%放电状态

3 程序结果

sop有功

sop无功

4 下载链接

点击直达！