考虑灵活性供需不确定性的储能优化配置(Matlab代码实现）

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

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

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

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

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码、数据、文档

---

💥1 概述

文献来源：

随着风电渗透率的增加，部分常规电源被风电替代，系统的灵活性调节能力不足，大规模风电的

随机波动性给系统带来较大的调峰压力[1]。储能能够实现对风电随机波动性的互补，提高系统的调峰灵活性。因此，在含大规模风电系统中配置储能，研究规划与运行中的灵活性问题，具有重要意义。

由于储能的规划与运行密切相关，文献[7-8] 提出储能配置的双层优化模型，外层优化储能的配

置方案，内层优化储能的充放电功率。高比例可再生能源系统中，灵活性的量化评估逐渐成为系统运行特性的核心和关键[9]，然而以上储能的优化配置模型主要从经济性角度研究储能的规划、运行问题，未考虑系统的灵活性，忽略了关键因素对储能配置方案的影响，不利于系统的安全、可靠运行。

目前，国内外学者针对系统灵活性及其指标的研究取得一定的进展，文献[10]通过统计法计算系

统的调峰不足概率，计算量较大。文献[11-12]采用蒙特卡罗模拟方法计算调峰不足概率/期望，该方法虽然具有较好的精度，但其时间复杂度较高[13]。文献 [14]考虑灵活性资源的优化调度和灵活性需求不确定性，建立储能的双层优化配置模型。但并未考虑机组的随机强迫停运，也忽略了机组的爬坡率对系统灵活性不足风险及储能配置方案的影响。为全面量化大规模风电并网对系统调峰

灵活 性的影响，本文考虑机组的随机停运、最大/最小出力限制以及不同出力状态下的向上/向下爬坡率等调峰能力的不确定性，提出调峰灵活性评估指标，综合考虑系统的经济性和灵活性，基于有效容量分布的时序随机生产模拟，建立储能辅助调峰的双层优化模型。上层储能配置模型以综合成本最小为目标，考虑风电出力的不确定性，将考虑多个随机场景的不确定性规划模型，转化为简单求解的多个确定性模型。下层运行优化模型在上层储能配置方案的基础上，考虑调峰需求和调峰能力的不确定性，以总调峰能力不足期望最小为目标，并将灵活性不足损失成本返回上层目标函数。

📚2 运行结果

部分代码：

%%系统最小技术出力为460，最大为1660，
P_D=[700;750;850;950;1000;1100;1150;1200;1040;1120;1160;1200;1120;1040;1200;1050;1000;1100;1200;1120;1040;1100;900;800];
delta_l_up=1660-max(P_D)-max(P_D)*0.2;%%调峰裕度
delta_l_do=460-min(P_D)+min(P_D)*0.2;
%%计算风电接入调峰能力的要求
peak_demand=zeros(100,24);
CNT=[];
for n=2:24
    peak_demand(:,n)=Wind_tree(:,n)-Wind_tree(:,n-1)+P_D(n)-P_D(n-1);
    
end

for i=1:100
    for j=1:24
    if peak_demand(i,j)>0&&peak_demand(i,j)>delta_l_up
        
CNT(i,j)=peak_demand(i,j)-delta_l_up;
    elseif peak_demand(i,j)<0&&peak_demand(i,j)<delta_l_do
        CNT(i,j)=peak_demand(i,j)-delta_l_do;
    else
       CNT(i,j)=0; 
    end
    end
end
%%统计调峰能力灵活指标
nn1=CNT;
nn1(nn1>0)=0;%%下调峰不足
nn2=CNT;
nn2(nn2<0)=0;%%上调峰不足
nn1=(nn1);
figure
subplot(2,2,1)
% edges = [0 1:0.25:2 24];
% histogram(abs(nn1),edges)
hist(nn1,100);
title('下调峰不足频率直方图');xlabel('缺额MW');ylabel('f');  %%new!
subplot(2,2,2)
[ni,ak]=hist(nn1,100);
fi=ni/length(nn1);
plot((ak),fi);
title('下调峰不足经验分布函数');xlabel('缺额MW');ylabel('F（x）');  %%new!
legend('经验分布函数');
subplot(2,2,3)
hist(nn2,100);
title('上调峰不足频率直方图');xlabel('缺额MW');ylabel('f');  %%new!
subplot(2,2,4)
[ni,ak]=hist(nn2,100);
fi=ni/length(nn2);
plot((ak),fi);
title('上调峰不足经验分布函数');xlabel('缺额MW');ylabel('F（x）');  %%new!
legend('经验分布函数');

🎉3 参考文献

文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。

[1]孙伟卿,宋赫,秦艳辉等.考虑灵活性供需不确定性的储能优化配置[J].电网技术,2020,44(12):4486-4497.DOI:10.13335/j.1000-3673.pst.2020.0667.

🌈4 Matlab代码、数据、文档