05基于多能互补的热电联供型微网优化运行

说明书

MATLAB代码：基于多能互补的热电联供型微网优化运行

关键词：多能互补 综合需求响应  热电联产 微网 优化调度

参考文档：《基于多能互补的热电联供型微网优化运行》基本完全复现；

仿真平台：MATLAB+yalmip+cplex

主要内容：代码主要做的是多能互补的热电联供型微网优化运行模型，在需求侧对负荷类型进行分类，利用电负荷的弹性和系统供热方式的多样性，构建含电负荷时移、削减响应及热负荷供能方式响应的综合能源需求响应模型，并提出响应补偿机制。在此基础上，以系统运行成本与响应补偿成本之和最小为目标，综合考虑供需双侧设备运行和可调度负荷资源约束，建立基于多能互补的 CHP-MG 优化运行数学模型。同时，为了体现模型的准确性，代码对比了热负荷参与、电负荷参与以及电热负荷均参与或者均不参与四种常见下的模型调度结果，体现了所构模型的经济型。

%方式3：
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];
bb=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,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<=12
        Cgas(i)=1.57;
    elseif  i>=19&&i<=22
        Cgas(i)=1.57;
    elseif i>=13&&i<=18
        Cgas(i)=2.05;
    else
        Cgas(i)=2.05;
    end
end
 %%一天分为24小时，时间步长取1小时/60min
%%%% 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]';