山西建设执业资格注册管理中心网站,深圳市建设交易中心网站首页,那些网站是做金融行业,专门网页制作工具前言
目前#xff0c;尽管不断地追逐可再生能源全额消纳方式#xff0c;大幅减小弃风弃光电量#xff0c;但是若考虑风电、光伏发电的随机属性#xff0c;全额消纳可能造成电网峰谷差、调峰难度及调峰调频等辅助服务费用的剧增#xff0c;引起电网潜在运行风险。因此尽管不断地追逐可再生能源全额消纳方式大幅减小弃风弃光电量但是若考虑风电、光伏发电的随机属性全额消纳可能造成电网峰谷差、调峰难度及调峰调频等辅助服务费用的剧增引起电网潜在运行风险。因此如果在风光大发与负荷低谷重叠时段选择弃掉一部分的风光尖峰电力不仅能够降低电网运行风险而且有助于提升电网调度的灵活性和整体经济效益。若系统中配置了一定容量的储能则能够在保证总体经济效益最大的前提下提高系统对于可再生能源的消纳能力实现最优弃风弃光率和弃电总量的双降减小净负荷曲线的波动程度和电网的调峰难度使火电机组出力更平稳。因此研究考虑最优弃能率的风光火储联合发电系统的优化调度方法对于实现风光的最经济消纳具有重要意义。
调峰过程及弃风弃光对电网调峰的影响
火电机组调峰过程及调峰成本模型
根据机组的出力状态火电机组的调峰阶段通常可划分为基本调峰阶段和深度调峰阶段而深度调峰阶段按照调峰深度和燃烧介质又可划分为投油深度调峰阶段和不投油深度调峰阶段如下图所示。其中Pmax为机组出力最大值Pmin为基本调峰阶段机组出力最小值Pn为不投油深度调峰阶段机组出力最小值Pm为投油深度调峰阶段机组出力最小值。 合理弃风弃光对电网调峰的影响分析
下图为1天24h内风电和光伏出力曲线。由下图可知风电出力随机性和波动性特征明显在日间负荷高峰时段出力偏低在凌晨和夜间这些时段出力较大具有明显的反调峰特征。若风电全额并网消纳将会增大受端电网净负荷峰谷差增加火电机组的调峰成本。然而光伏出力的特征则与风电完全相反凌晨和夜间时段的光伏出力为零而日间的光伏出力先增大后减小在正午时段出力达到最大值。若光伏全额并网消纳则会降低受端电网正午时段的净负荷及其峰谷差从而降低电网调峰难度增大火电机组调度的灵活性。 下图给出了不同情景下受端电网净负荷曲线及风光并网功率的变化情况。由图可知由于风电和光伏在时间尺度上具有较为明显的互补特性两者之和往往在一定范围内上下浮动。因此风电和光伏两者同时全额并网消纳并不一定能够减小火电机组的调峰成本实现电力系统的最优运行。若在某些时段弃掉一定比例的风光电力尤其是低谷负荷时段则会明显减小净负荷曲线的峰谷差和机组参与深度调峰的时段数。此外由于储能可在低负荷时段充电在峰荷时段放电具有较好的能量时移特性因此在消纳风光电力的同时配置一定容量的储能可以进一步降低受端电网净负荷的峰谷差使净负荷曲线更为平滑从而进一步降低火电的调峰成本。 考虑最优弃能率的双层优化调度模型
程序以寻求各时段的最优弃风弃光率为切入点并利用储能系统良好的能量时移特性构建了考虑最优弃能率的风光火储联合系统双层优化调度模型其结构如下图所示。 上层优化模型以受端电网净负荷方差最小为目标函数在满足弃风弃光率约束及储能电站运行约束的基础上对各时段的弃风弃光率及储能系统充放电功率进行优化由此得出各时段最优风光并网功率及储能充放电功率。同时将考虑最优弃能率的风光并网功率曲线及储能系统充放电功率曲线传递至下层优化模型。
下层优化模型基于上层模型传递的最优风光并网功率曲线及储能系统充放电功率曲线并结合净负荷曲线以受端电网运行成本最小为目标对各时段火电机组出力进行优化。目标函数中同时考虑了多项成本包括火电机组调峰成本、风电、光伏和储能电站运行成本、系统旋转备用成本、环境保护税、风光并网消纳及储能系统的环境收益通过优化求解得出受端电网火电机组的最优出力及最小运行成本。
程序介绍
程序提出了一种考虑最优弃能率的分层优化调度策略主要研究合理弃风弃光对电网调峰难度及整体经济性的影响。为此首先分析了火电机组调峰过程及合理弃风弃光对电网调峰的影响其次综合考虑电网调峰难度及整体经济效益提出了最优弃能率的约束模型并构建了以电网净负荷方差最小为上层目标函数、电网运行成本最小为下层目标函数的双层优化调度模型选取IEEE-30节点系统进行算例分析。程序中算例丰富注释清晰干货满满创新性和可扩展性很高足以撑起一篇高水平论文下面对程序做简要介绍
程序适用平台MatlabYalmipGurobi
程序结果 部分程序
%% 定义决策变量 para2;R_wt sdpvar(1,24);R_pv sdpvar(1,24); P_c sdpvar(1,24); P_d sdpvar(1,24); P_net sdpvar(1,24); %净负荷P_wt1sdpvar(1,24);%并网风力发电P_pv1sdpvar(1,24);%并网光伏发电 %其他变量 Q_s sdpvar(1,24); % 蓄电池剩余电量
% P_Gt1sdpvar(1,24); % 燃气发电机1提供的功率
% P_Gt2sdpvar(1,24); % 燃气发电机2提供的功率
%% 目标函数
P_net(t)P_f(t)-(1-R_wt(t))*P_wt(t)-(1-R_pv(t))*P_pv(t)...
h_c*P_c(t)-h_d*P_d(t);
sumsumP_f(t)-(1-R_wt(t))*P_wt(t)-(1-R_pv(t))*P_pv(t)...
sumsum./24;
%负荷平方差
objectiveobjective((P_f(t)-(1-R_wt(t))*P_wt(t)-(1-R_pv(t))*P_pv(t)...
h_c*P_c(t)-h_d*P_d(t))-sum).^2;objectiveobjective./24;
%% 约束条件
constraint[constraint,0R_wt(t)1];
constraint[constraint,0R_pv(t)1];
for t1:24 %(6-1) 储能系统约束
constraint[constraint, 40Q_s(t)360];
constraint[constraint, Q_s(t)Q_s_initP_c(t)-P_d(t)];
constraint[constraint, Q_s(t)-Q_s(t-1)-P_c(t)P_d(t)0];
% 其他不等式约束
constraint[constraint, 0P_wt(t)400] ;
AAP_wt(t)*R_wt(t);BBP_pv(t)*R_pv(t);
CCP_wt(t);DDP_pv(t);
% 弃风弃光约束%储能系统充放电约束for t1:24 %(6-2) 蓄电池不等式约束constraint[constraint,0P_c(t)P_c_max];
% (6-3) 电池储能初始和最终状态相等约束constraint [constraint, Q_s_init-Q_s(24)0 ];
%% 求解和显示
option sdpsettings(verbose,1,solver,gurobi);%调用gurobi求解器
optimize(constraint,objective,option)%Yalmip求解的命令
部分内容源自网络侵权联系删除
欢迎感兴趣的小伙伴关注并私信获取完整版代码小编会不定期更新高质量的学习资料、文章和程序代码为您的科研加油助力
