目录
主要内容
该模型构建了含冷、热、电负荷的冷热电联供型综合能源系统优化调度模型,设备包括燃气轮机、电制冷机、锅炉、风光出力以及购售电交易,以用户购电购热冷量的成本、CCHP收益以及碳排放成本为目标,采用权值多目标方式将多目标转化为单目标求解,同时将约束采用罚函数形式进行表达,实现系统的经济运行,程序采用matlab编写,注释清晰,方便学习。
部分代码
PGnom1=1200;PGnom2=2800;%电功率基准值为1200和2800 Copac=0.80;Copec=3; %-------------% for t=1:24 ygg(1,t)=(Aa+Bb*(X((t-1)*D+1)/PGnom1)+Cc*(X((t-1)*D+1)/PGnom1)^2+Dd*(X((t-1)*D+1)/PGnom1)^3)/100; %X((t-1)*D+1)在对应的该时刻的该燃气发电机的电功率 ygg(2,t)=(Aa+Bb*(X((t-1)*D+2)/PGnom2)+Cc*(X((t-1)*D+2)/PGnom2)^2+Dd*(X((t-1)*D+2)/PGnom2)^3)/100; %ygg对应的是第i台燃气发电机的发电效率 end %-------------% for t=1:24 %R是燃气机发电模型的λ FG(t)=R*(X((t-1)*D+1)/ygg(1,t)+X((t-1)*D+2)/ygg(2,t)); %燃气发电机在单时段的燃气耗量 %R疑似是锅炉模型的λ FBoi(t)=R*(X((t-1)*D+3)+X((t-1)*D+4))/yBoi; %锅炉在单时段消耗的燃料量 FE(t)=FG(t)+FBoi(t); %t时刻的燃气量 if X((t-1)*D+5)>=0 Rr(t)=rElc(t); PG(t)=X((t-1)*D+5); %向主电网交互的电量,此时是买电 else Rr(t)=RElc(t); PG(t)=0; %向主电网交互的电量,此时是卖电 end CE(t)=Rr(t)*X((t-1)*D+5); %单时段与主电网间能量交互成本 HJ(t)=ue*PG(t)+uf*FE(t); %二氧化碳的排放量 %---单时段功率平衡约束---% h(t)=abs(X((t-1)*D+1)+X((t-1)*D+2)+Pwind(t)+Ppv(t)-X((t-1)*D+6)+X((t-1)*D+5)-PD(t)-X((t-1)*D+8));%电负荷平衡约束 end %---用的燃气量不同,得到的价格也不同---% if sum(FE)<250 rLrg=5.257; elseif 250<=sum(FE)<=4167 rLrg=3.25; else rLrg=2.814; end %-----------------------------------% f1=rLrg*sum(FE)+sum(CE); %燃气成本和向电网购电或是售电后的花费 f2=W*sum(HJ); %碳排放消费费用