【微电网调度】基于粒子群实现对含并网蓄电池微电网的优化调度附matlab代码
智能优化算法 神经网络预测 雷达通信 无线传感器 电力系统
信号处理 图像处理 路径规划 元胞自动机 无人机
⛄ 内容介绍
随着能源革命的进行,利用风,光资源的开发与微电网应用的研究受到了各国科学家的关注.由于分布式发电技术本身在电力电子技术方面不成熟且发电的不稳定性并不遵循正常的人工调度,因此微电网发电技术深受人们诟病.文章以风力发电机,光伏机组和燃气轮机为基础,以蓄电池为储能装置,建立了一系列的微电网系统,并使用粒子群算法对发电成本进行迭代从而实现最优解.
⛄ 部分代码
cost=rid(PGT)
%适应值函数
%--------------------------------------------------------------------------
T=24;
N=5;
% 输入原始数据,包括各时段负荷大小，发电机有功输出上下限，发电机耗量成本系数，各时段风电场预测的平均输出功率
pmax(1)=115;pmin(1)=0;
pmax(2)=40;pmin(2)=10;
pmax(3)=40;pmin(3)=10;
pmax(4)=30;pmin(4)=0;
% 污染物排放成本
% aa表示NOX，bb表示so2,cc表示co2,从左到右依次为Grid,DG，MT，FC
aa=[3.6,21.8,0.03,0.44];
bb=[4.54,0.454,0.006,0.008];
cc=[23,1.432,1.078,1.596];
% 柴油发电机的耗量参数
a=0.4333;b=0.2333;c=0.0071;
% 微型燃气轮机,P_mt表示燃气轮机发出的功率，Xl_mt表示燃气轮机的效率
Xl_mt=0.5;
=0.39;
L=9.7; % 天然气的低热值，单位是kW.h/m3,千瓦时每立方米
% Xl_mt=0.0753*(P_mt/65)^3-0.3095*(P_mt/65)^2+0.1068;
% 燃料电池
Xl_fc=0.6;
=0.39;
=0.07; % 电网平均电价
=0.07; % 微网多余电能销售给电网的电价
% =[0.01 0.01 0.01 0.01 0.01 0.04 0.04 0.04 0.04 0.04 0.04 0.1 0.1 0.1 0.04 0.04 0.04 0.1 0.1 0.1 0.1 0.04 0.01 0.01];
% =[0.03 0.03 0.03 0.03 0.03 0.06 0.06 0.06 0.06 0.06 0.06 0.09 0.09 0.09 0.06 0.06 0.06 0.09 0.09 0.09 0.09 0.06 0.03 0.03];
%%%%%%%%%%=[ 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24]
%--------------------------------------------------------------------------
C=zeros(1,T);
Cost=0;
for jj=1:T
% s=a+b*PGT(1,jj)+c*PGT(1,jj)*PGT(1,jj)+*PGT(2,jj)/Xl_fc+*PGT(3,jj)/(0.0753*(PGT(3,jj)/65)^3-0.3095*(PGT(3,jj)/65)^2+0.4174*PGT(3,jj)/65+0.1068);
% % 目标1
if PGT(1,jj)>=0
s(jj)=a+b*PGT(4,jj)+c*PGT(4,jj)*PGT(4,jj)+(/L)*PGT(2,jj)/Xl_fc+(/L)*PGT(3,jj)/Xl_mt+*PGT(1,jj);
else
s(jj)=a+b*PGT(4,jj)+c*PGT(4,jj)*PGT(4,jj)+(/L)*PGT(2,jj)/Xl_fc+(/L)*PGT(3,jj)/Xl_mt+*PGT(1,jj);
end
Cost=Cost+s(jj);
% 目标2
% if PGT(1,jj)>=0
% s(jj)=0.0042*(aa(1)* PGT(1,jj)+aa(2)*PGT(4,jj)+aa(3)*PGT(3,jj)+aa(4)*PGT(2,jj))+0.00099*(bb(1)*PGT(1,jj)+bb(2)*PGT(4,jj)+bb(3)*PGT(3,jj)+bb(4)*PGT(2,jj))+0.*(cc(1)*PGT(1,jj)+cc(2)*PGT(4,jj)+cc(3)*PGT(3,jj)+cc(4)*PGT(2,jj));
% else
% s(jj)=0.0042*(aa(2)*PGT(4,jj)+aa(3)*PGT(3,jj)+aa(4)*PGT(2,jj))+0.00099*(bb(2)*PGT(4,jj)+bb(3)*PGT(3,jj)+bb(4)*PGT(2,jj))+0.*(cc(2)*PGT(4,jj)+cc(3)*PGT(3,jj)+cc(4)*PGT(2,jj));
% end
% Cost=Cost+s(jj);
% % 目标3
% if PGT(1,jj)>=0
% s1(jj)=a+b*PGT(4,jj)+c*PGT(4,jj)*PGT(4,jj)+(/L)*PGT(2,jj)/Xl_fc+(/L)*PGT(3,jj)/Xl_mt+*PGT(1,jj);
% s2(jj)=0.0042*(aa(1)* PGT(1,jj)+aa(2)*PGT(4,jj)+aa(3)*PGT(3,jj)+aa(4)*PGT(2,jj))+0.00099*(bb(1)*PGT(1,jj)+bb(2)*PGT(4,jj)+bb(3)*PGT(3,jj)+bb(4)*PGT(2,jj))+0.*(cc(1)*PGT(1,jj)+cc(2)*PGT(4,jj)+cc(3)*PGT(3,jj)+cc(4)*PGT(2,jj));
% else
% s1(jj)=a+b*PGT(4,jj)+c*PGT(4,jj)*PGT(4,jj)+(/L)*PGT(2,jj)/Xl_fc+(/L)*PGT(3,jj)/Xl_mt+*PGT(1,jj);
% s2(jj)=0.0042*(aa(2)*PGT(4,jj)+aa(3)*PGT(3,jj)+aa(4)*PGT(2,jj))+0.00099*(bb(2)*PGT(4,jj)+bb(3)*PGT(3,jj)+bb(4)*PGT(2,jj))+0.*(cc(2)*PGT(4,jj)+cc(3)*PGT(3,jj)+cc(4)*PGT(2,jj));
% end
% Cost=Cost+s1(jj)+s2(jj);
end
cost=Cost;
⛄ 运行结果
⛄ 参考文献
[1] 程宇旭.基于改进粒子群算法的微电网能量优化调度研究及实现[D].中南大学,2014.
