【EI复现】参与调峰的储能系统配置方案及经济性分析（Matlab代码实现）

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

1.1 储能辅助电网调峰原理分析

1.2 系统调峰优化指标模型

1.3 储能系统容量配置模型

📚2 运行结果

2.1 模式设置

2.2 新能源消纳曲线：（模式1-模式3） 

2.3 场景3下的运行结果 

2.4 储能配置结果 

2.5 两个灵敏度分析 

🎉3 文献来源

🌈4 Matlab代码、数据、文章讲解

---

​

💥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.

🌈4 Matlab代码、数据、文章讲解