1主要内容
该程序参考文献《考虑火电机组储热改造的电力系统低碳经济调度》，利用原文献火电机组储热改造方案建立模型，在传统火电机组的基础上加装热能存储系统，根据储热装置的充热、放热运行机理，构建火电机组储热改造出力模型，以煤耗成本、风光储运行成本、备用成本、环境收益、碳交易成本以及机组启动成本等综合成本为目标，以系统平衡约束、高/低温储热系统出力和容量约束、旋转备用约束、机组爬坡约束等作为模型约束条件，采用matlab+gurobi作为求解平台，程序注释清晰，方便研究学习！
储热改造原理
火电机组储热改造原理如下图：
具体的改造原理不再赘述，但是有一点值得我们重视，就是创新点，这种机器设备改造方面的创新点要比方法上的创新点容易的多，但是这种创新点不管是从创新的质量和实用性上都不错，所以对于搞研究的同学们，下沉到最一线考虑创新点，比仅仅关注方法层面的改进要容易很多。
约束条件
1.高低温储热系统约束
这个约束条件和常见的储能约束条件是一致的，将一个常规模型应用于设备改造上面就成为该文献的最大创新点，并且理论水平和求解过程并不复杂，这就是我们该从中学到的一个重点。
2.火电机组内部能流约束
3.功率平衡约束
由于本模型未考虑弃负荷和弃风光，因此功率平衡和该公式略有出入，但是整体思路是一致的。
4.旋转备用约束
5.其他约束
除此之外，还包括机组最小启停时间约束、爬坡约束、高/低温储热系统储/放热功率、状态约束等。
除约束外，程序还需根据变量运算进行线性化处理，以保证求解的准确性。
2部分程序
%%系统功率平衡约束for t=1:24constraint=[constraint,0.4*(P_1_stgen(t)+P_2_stgen(t)+P_3_stgen(t)+P_4_stgen(t))-P_net(t)==0];%为什么*0.4，燃烧热效率end%%系统旋转备用功率约束for t=1:24constraint=[constraint,u_d*P_D(t)+u_w*P_wt(t)+u_pv*P_pv(t)-(200-0.4*P_1_stgen(t))-(300-0.4*P_2_stgen(t))-(300-0.4*P_3_stgen(t))-(600-0.4*P_4_stgen(t))<=0];constraint=[constraint,u_d*P_D(t)+u_w*P_wt(t)+u_pv*P_pv(t)-(0.4*P_1_stgen(t)-0.4*200)-(0.4*P_2_stgen(t)-0.5*300)-(0.4*P_3_stgen(t)-0.5*300)-(0.4*P_4_stgen(t)-0.3*600)<=0];end%%火电机组爬坡约束for t=1:24if t==1constraint=[constraint,0.4*P_1_stgen(t)<=R_gup(1)];constraint=[constraint,0.4*P_2_stgen(t)<=R_gup(2)];constraint=[constraint,0.4*P_3_stgen(t)<=R_gup(3)];constraint=[constraint,0.4*P_4_stgen(t)<=R_gup(4)];else火电机组1爬坡约束constraint=[constraint,-R_gdomn(1)<=0.4*P_1_stgen(t)-0.4*P_1_stgen(t-1)<=R_gup(1)];
3程序结果
3下载链接
扫码或阅读原文即可跳转程序页面：
想了解更多程序？按照下面教程查看程序目录！
