【调峰】储能辅助电力系统调峰的容量需求研究（Matlab代码实现）

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

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

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

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

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

---

💥1 概述

文献来源：

合理选取储能容量，使储能系统实现经济效益最优和可再生能源发电利用率最大化，是满足调峰需水的惦衣且日忠能灾量需求评估的方法主要分为三类，如图1-1的优化配置问题进行研究。目前针对储能容量需求评估的方法主要分为三类，如图1-1所示。

理论分析法包括时域分析法和频域分析法两个方面，时域分析是基于概率分析或者一阶滤波的方法对储能的容量需求进行分析，频域分析则采用傅里叶变换的方式进行分析。文献利用离散傅里叶变换对可再生能源输出功率进行频谱分析，考虑储能系统充放电效率、荷电状态及可再生能源发电系统目标功率输出波动率的约束，确定能够平抑新能源发电出力的最小储能系统容量;文献[33]提出了一种光电与建筑一体化(BIPV)的储能容量优化的频谱分析方法，根据平滑电网功率波动所需功率的频谱分析结果，确定储的控制容量需求;文献[34]利用离散傅里叶变换(DFT)对风电场输出功率偏差进行频谱分析，得到风电场对各类时间响应的控制需求量，采取不同的控制机制以补偿功率偏差。理论分析法的求解速度较快，但其不一定能够满足技术性的指标，也未考虑储能系统运行的经济性。

📚2 运行结果

 

 

部分代码：

%% 连续变量
% 决策变量：火电出力, 风电出力, 调峰不足,u储能放电，u储能充电 ，储能容量,
P_fire = sdpvar(1,24);       % 火电出力出力
P_wind = sdpvar(1,24); %风电出力
P_loss = sdpvar(1,24);  %调峰不足

%储能
P_ch = sdpvar(1,24);   % 蓄电池充电功率
P_dis = sdpvar(1,24);   % 蓄电池放电功率
Uchr = binvar(1,24);   % 电池充电状态，1表示充电
U_dis = binvar(1,24);   % 电池放电状态，1表示放电
W_Bt= sdpvar(1,24);     % 蓄电池的荷电状态 
H_SS= sdpvar(1);      % 蓄电池容量 

%% 约束
constraint=[];
for t=1:24 % (1) 电能平衡约束    
    constraint = [constraint, P_fire(t)+P_wind(t)+P_loss(t)+P_dis(t)-P_ch(t)-Pload(t)==0 ]; 
end
constraint=[constraint, 0 <= H_SS <= E_max];
for t=1:24 %  蓄电池不等式约束
    constraint=[constraint, 0<=P_ch(t)<=Uchr(t)*E_max];
    constraint=[constraint, 0<=P_dis(t)<=U_dis(t)*E_max];
    constraint=[constraint, Uchr(t)+U_dis(t)<=1];
    constraint=[constraint, H_SS*soc_min<=W_Bt(t)<=H_SS*soc_max];
end
for t=1:24 % 蓄电池等式约束
     if t==1
        constraint=[constraint, W_Bt(t)==H_SS*soc0+P_ch(t)*eta-P_dis(t)/eta];
    else
        constraint=[constraint, W_Bt(t)==W_Bt(t-1)+P_ch(t)*eta-P_dis(t)/eta];
    end
end
 constraint = [constraint, H_SS*soc0==W_Bt(24) ]; 

%% 火电上下限约束
for t=1:24 
    constraint=[constraint, pf_min<=P_fire(t)<=pf_max] ; 
end
%% 火电爬坡约束
for t=2:24 
    constraint=[constraint, -delta<=P_fire(t)-P_fire(t-1)<=delta] ; 
end
%% 风电上下限约束
for t=1:24 
    constraint=[constraint, 0 <= P_wind(t) <= Pwind(t)] ; 
end

%% 调峰不足约束
for t=1:24
 constraint = [constraint, 0<= P_loss(t) <=Pload]; 
end

🎉3 参考文献

部分理论来源于网络，如有侵权请联系删除。

[1]汪赛. 储能辅助电力系统调峰的容量需求研究[D].广西大学,2018.

🌈4 Matlab代码实现