简介
本程序参考EI论文《基于多能互补的热电联供型微网优化运行》，文章通过储能设备解耦热电联系，建立基于多能互补的综合能源系统/虚拟电厂/微电网优化运行模型。模型包含系统供给侧的多能互补协调与需求侧的综合能源响应两个方面，使供给侧通过能源转换设备扩充供能能力，使需求侧通过改变用户用能选择来提高响应能力，为微网的热、电生产提供优化空间。最后，采用CPLEX软件进行模型求解，程序中算例丰富、注释清晰、干货满满，复现时预留了扩展设备结构，创新性很高！下面对文章和程序做简要介绍！
创新点及主要工作
创新点
1. 利用储能实现热电解耦；2. 可转移、可削减负荷及热负荷需求响应；
3. 多能互补、梯级利用；
主要工作
热电联供型微网(CHP-MG)对实现能源可持续发展和构建绿色低碳社会具有重要的应用价值，而内部复杂的能源结构与设备耦合关系，也对其运行优化带来了挑战。利用供需双侧电、热能的互动互补关系，在供给侧采用储能装置实现联供设备的热电解耦，通过各能源转换设备提升系统多能源的供应能力。在需求侧对负荷类型进行分类，利用电负荷的弹性和系统供热方式的多样性，构建含电负荷时移、削减响应及热负荷供能方式响应的综合能源需求响应模型，并提出响应补偿机制。在此基础上，以系统运行成本与响应补偿成本之和最小为目标，综合考虑供需双侧设备运行和可调度负荷资源约束，建立基于多能互补的 CHP-MG 优化运行数学模型。基于算例的仿真结果和对比分析表明：考虑多能互补的供需双侧协同优化能有效提高系统供能的灵活性以及运行经济性。
文中结果
程序结果
上述两图片程序中预留了可扩展设备接口！！
部分程序
%方式5：储能解耦电热关系，只有热需求响应；%10个高概率风电功率w=2*fix([38,37,27,29,23,14,21,13,43,76,59,70,49,41,51,41,28,21,18,18,33,41,45,42;Pwind=w(5,:);% 电热负荷load=[178 160 200 220 258 260 286 380 380 385 360 320 280 254 266 285 312 335 326 288 200 185 146 120];%电用户% Phot=[360,353,342,342,331,328,287,265,240,245,244,238,233,229,232,237,249,257,266,270,301,321,341,366];Phot=[350,343,332,332,331,318,317,285,260,245,244,238,233,249,262,277,299,267,266,270,301,321,341,346];bb=[0,0,0,0,10,0,0,0,0,50,0,0,0,0,0,0,0,50,0,0,0,0,0,0]; %这个不用管% 日前计划 可控负荷% 可中断负荷L5=[19,24,24,19,24,29,64,49,64,59,59,59,49,59,59,39,39,64,64,19,44,19,44,19]; %可中断% 可转移负荷L6=[18,8,13,13,18,13,13,18,23,33,23,18,18,18,18,18,18,18,13,13,13,8,8,8];%可平移% 供热选择量--电热量L7=[17,22,22,17,22,26,58,85,88,84,84,84,65,64,74,55,35,58,58,17,40,17,40,17];% 购气%%%% 1台MT机组,1台LHV=9.75; %电热转换系数aF=[0.5869 0.3952 ]';%参数bF=[8.6204 -0.185]';%参数aH=[1.377]';%参数bH=[20.38]';%参数Pgrid_max=[320 320]';Pgrid_min=[0 0]';%% Pnas 文中PessPnas_max=[60 60]'; Pnas_min=[0 0]';%%%燃气锅炉参数H_max=600;H_n=0.98%锅炉容量/转换率Peh_max=60;n_Peh=0.98;%转换设备/转换率H_storage_max=800; h_n=0.98;h_charge=0.9;h_discharge=1;%热储能容量/自损/充热/放热；%%%储能参数% 容量 1000 初始500SOC_max=[800]';SOC_min=[200]';SOC_ini=0.5;%% 峰 平 谷 电价buy=[1.243 0.8934 0.47];%谷/峰 购电价%% 赋值MT=intvar(1,24,'full');Pgrid=intvar(1,24,'full'); % 购 、 售Pnas=intvar(2,24,'full'); % 充、 放H=intvar(1,24,'full');%锅炉Q=intvar(1,24,'full');%可中断电负荷P=intvar(1,24,'full');%转热负荷Hti=intvar(1,24,'full');%充热Hto=intvar(1,24,'full');%放热xx=intvar(1,24,'full');%%%%%时平负荷量
部分内容源自网络，侵权联系删除！
欢迎感兴趣的小伙伴点击文末阅读原文获取完整版代码，小编会不定期更新高质量的学习资料、文章和程序代码，为您的科研加油助力！
更多创新性预测及优化程序请点击公众号首页：《预测优化》|《综合能源》链接！
