本程序参考EI论文《基于多能互补的热电联供型微网优化运行》，文章通过储能设备解耦热电联系，建立基于多能互补的综合能源系统/虚拟电厂/微电网优化运行模型。模型包含系统供给侧的多能互补协调与需求侧的综合能源响应两个方面，使供给侧通过能源转换设备扩充供能能力，使需求侧通过改变用户用能选择来提高响应能力，为微网的热、电生产提供优化空间。最后，采用CPLEX软件进行模型求解，程序中算例丰富、注释清晰、干货满满，复现时预留了扩展设备结构，创新性很高！下面对文章和程序做简要介绍！

 

文章亮点：

1. 利用储能实现热电解耦； 2. 可转移、可削减负荷及热负荷需求响应；3. 多能互补、梯级利用；

文章主要工作：

热电联供型微网(CHP-MG)对实现能源可持续发展和构建绿色低碳社会具有重要的应用价值，而内部复杂的能源结构与设备耦合关系，也对其运行优化带来了挑战。利用供需双侧电、热能的互动互补关系，在供给侧采用储能装置实现联供设备的热电解耦，通过各能源转换设备提升系统多能源的供应能力。在需求侧对负荷类型进行分类，利用电负荷的弹性和系统供热方式的多样性，构建含电负荷时移、削减响应及热负荷供能方式响应的综合能源需求响应模型，并提出响应补偿机制。在此基础上，以系统运行成本与响应补偿成本之和最小为目标，综合考虑供需双侧设备运行和可调度负荷资源约束，建立基于多能互补的 CHP-MG 优化运行数学模型。基于算例的仿真结果和对比分析表明：考虑多能互补的供需双侧协同优化能有效提高系统供能的灵活性以及运行经济性。

文中结果：

 程序结果：

部分程序：

 %方式5：储能解耦电热关系，只有热需求响应；
clc;
clear all;
%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;34,38,33,29,26,13,19,13,38,40,59,78,51,38,54,43,24,21,22,22,25,39,48,30;30,30,27,35,27,13,18,14,39,50,58,84,58,39,47,44,25,20,22,18,28,41,36,46;39,32,32,34,23,14,20,14,43,50,61,73,60,41,48,43,26,21,22,20,29,40,43,30;32,39,28,30,21,12,22,16,37,62,53,77,60,39,50,49,25,20,21,18,30,40,38,48;34,43,26,30,22,18,21,12,41,56,65,78,47,37,50,45,23,21,21,20,27,40,39,40;40,43,32,30,24,16,21,15,41,65,62,69,55,44,53,42,27,21,20,21,28,42,41,50;39,40,26,28,25,14,21,18,39,70,54,68,51,38,56,43,25,17,21,17,28,36,40,35;41,43,30,30,22,13,20,13,46,66,55,79,55,37,48,37,28,18,22,19,27,41,41,48;41,38,30,27,23,15,17,14,39,61,62,78,54,39,46,45,22,20,20,20,28,35,36,52]);
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];
% 购气
for i=1:24                   if i>=7&&i<=12Cgas(i)=1.57;elseif  i>=19&&i<=22Cgas(i)=1.57;elseif i>=13&&i<=18Cgas(i)=2.05;elseCgas(i)=2.05;end
end
%%%% 1台MT机组,1台
LHV=9.75; %电热转换系数
aF=[0.5869 0.3952 ]';%参数
bF=[8.6204 -0.185]';%参数
aH=[1.377]';%参数
bH=[20.38]';%参数
%% PCC
Pgrid_max=[320 320]';
Pgrid_min=[0 0]';
%% Pnas 文中Pess
Pnas_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  初始500
SOC_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');%%%%%时平负荷量
yy=intvar(1,24,'full');

欢迎感兴趣的小伙伴关注，小编会不定期更新高质量的学习资料、文章和程序代码，为您的科研加油助力！