👨🎓个人主页：研学社的博客
💥💥💞💞欢迎来到本博客❤️❤️💥💥
🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。
⛳️座右铭：行百里者，半于九十。
📋📋📋本文目录如下：🎁🎁🎁
目录
💥1 概述
1.1 储能辅助电网调峰原理分析
目前电网的调峰形势为在负荷尖峰时段有足够的旋转备用空间，但在负荷低谷时期，机组的向下
调节灵活性严重不足，从而导致大量弃风产生。储能辅助火电机组调峰基本原理如图 1 所示。
从图 1 中可以看出，储能辅助火电机组深度调峰可以有效改善电网调峰压力，减少弃风产生。而
储能系统产生的调峰效果主要取决于其配置方案，配置越高其调峰效果越好，但成本也随之大幅上升。储能系统的配置应兼顾经济性指标及技术性指标。
1.2 系统调峰优化指标模型
1.3 储能系统容量配置模型
详细数学模型及讲解见第4部分。
📚2 运行结果
2.1 模式设置
基础调峰：50%出力，Pmin=0.5*Pmax
深度调峰：40%出力，Pmin=0.4*Pmax, 600mw可以深度调峰，200-300不可以（50%出力）。
投油，30%出力。
场景设置
2.2 新能源消纳曲线：（模式1-模式3）
2.3 场景3下的运行结果
2.4 储能配置结果
2.5 两个灵敏度分析
部分代码：
%储能电池（对应公式30）
summ2=0;
for l=1:1:24
summ2=summ2+max(x2(1,l),0)*ng+min(x2(1,l),0)/np;
end
c =[c,summ2==0];%容量平衡约束
summ4=E0;
for l=1:1:24
summ4=summ4+max(x2(1,l),0)*ng+np*min(x2(1,l),0);
c =[c,Ehmax*0.1<=summ4<=Ehmax*0.9];%容量上下限约束
end
for l=1:1:24
summ9=max(x2(1,l),0)*ng+np*min(x2(1,l),0);
c =[c,-Pgmax<=summ9<=Pgmax];%功率上下限约束
end
%风机
c =[c,0<=x3<=fj];%可再生能源消纳约束
%日负荷平衡约束
for i=1:24
c =[c,sum(x1(:,i))-sum(x2(:,i))+sum(x3(:,i))==load(i)];%
end
%% 目标函数
%储能设备（公式1）
C1=(CE*1000*Ehmax+CP*1000*Pgmax)*(rs*(1+rs)^Nz)/(((1+rs)^Nz)-1)/(365*Nz);
%储能设备（公式2）
C2=0;
for i=1:24
C2=C2+abs(x2(1,i))*0.0075*1000;%储能充放电成本
end
%火电机组（公式3）
C3=0;
for j=1:6
for i=1:24
C3=C3+ag(j)*x1(j,i)*x1(j,i)+bg(j)*x1(j,i)+cg(j);%机组发电成本
end
end
%火电机组（公式4-5）
C4=0;
for j=1:6
for i=1:24
C4=C4+0.05*1000*abs(x1(j,i));%机组发电成本
end
end
C5=0;
% for j=1:2
% for i=1:24
% C5=C5+b2(j,1)*4.8*6130;%投油成本
% end
% end
%弃风成本（公式6）
C6=0;
for i=1:24
C6=C6+0.61*1000*max(fj(i)-x3(1,i),0);%
end
%储能套利（公式7）
C7=0;
for i=1:24
C7=C7+Pgs(i)*(max(x2(1,i),0)*ng+np*min(x2(1,i),0));
end
%储能补偿（公式8）
C8=0;
for i=1:24
C8=C8+0.4*max(x2(1,i),0)*ng;%调单位补偿
end
%火电机组补偿（公式9）
C9=0;
for j=1:6
for i=1:24
C9=C9+0.01*(Pmin1(j)-x1(j,i));%max(0,Pmin1(j)-x1(j,i)),求解时间变长
end
end
%新能源消纳（公式15）
C10=0;
for i=1:24
C10=C10+0.05*(x3(1,i)-xn(i));%指标归一化
end
%调峰成本（公式28-29）
GD=5000;%基础调峰成本
GC=5000*1.2;%不投油深度调峰成本
GO=4.8*6130;%投油深度调峰成本%储能电池（对应公式30）
summ2=0;
for l=1:1:24
summ2=summ2+max(x2(1,l),0)*ng+min(x2(1,l),0)/np;
end
c =[c,summ2==0];%容量平衡约束
summ4=E0;
for l=1:1:24
summ4=summ4+max(x2(1,l),0)*ng+np*min(x2(1,l),0);
c =[c,Ehmax*0.1<=summ4<=Ehmax*0.9];%容量上下限约束
end
for l=1:1:24
summ9=max(x2(1,l),0)*ng+np*min(x2(1,l),0);
c =[c,-Pgmax<=summ9<=Pgmax];%功率上下限约束
end
%风机
c =[c,0<=x3<=fj];%可再生能源消纳约束
%日负荷平衡约束
for i=1:24
c =[c,sum(x1(:,i))-sum(x2(:,i))+sum(x3(:,i))==load(i)];%
end
%% 目标函数
%储能设备（公式1）
C1=(CE*1000*Ehmax+CP*1000*Pgmax)*(rs*(1+rs)^Nz)/(((1+rs)^Nz)-1)/(365*Nz);
%储能设备（公式2）
C2=0;
for i=1:24
C2=C2+abs(x2(1,i))*0.0075*1000;%储能充放电成本
end
%火电机组（公式3）
C3=0;
for j=1:6
for i=1:24
C3=C3+ag(j)*x1(j,i)*x1(j,i)+bg(j)*x1(j,i)+cg(j);%机组发电成本
end
end
%火电机组（公式4-5）
C4=0;
for j=1:6
for i=1:24
C4=C4+0.05*1000*abs(x1(j,i));%机组发电成本
end
end
C5=0;
% for j=1:2
% for i=1:24
% C5=C5+b2(j,1)*4.8*6130;%投油成本
% end
% end
%弃风成本（公式6）
C6=0;
for i=1:24
C6=C6+0.61*1000*max(fj(i)-x3(1,i),0);%
end
%储能套利（公式7）
C7=0;
for i=1:24
C7=C7+Pgs(i)*(max(x2(1,i),0)*ng+np*min(x2(1,i),0));
end
%储能补偿（公式8）
C8=0;
for i=1:24
C8=C8+0.4*max(x2(1,i),0)*ng;%调单位补偿
end
%火电机组补偿（公式9）
C9=0;
for j=1:6
for i=1:24
C9=C9+0.01*(Pmin1(j)-x1(j,i));%max(0,Pmin1(j)-x1(j,i)),求解时间变长
end
end
%新能源消纳（公式15）
C10=0;
for i=1:24
C10=C10+0.05*(x3(1,i)-xn(i));%指标归一化
end
%调峰成本（公式28-29）
GD=5000;%基础调峰成本
GC=5000*1.2;%不投油深度调峰成本
GO=4.8*6130;%投油深度调峰成本
🎉3 文献来源
部分理论来源于网络，如有侵权请联系删除。
[1]李军徽,张嘉辉,李翠萍等.参与调峰的储能系统配置方案及经济性分析[J].电工技术学报,2021,36(19):4148-4160.DOI:10.19595/j.cnki.1000-6753.tces.200678.