【代码分享】考虑储能设备容量及季节性的电热综合能源系统优化调度

程序名称：考虑储能设备容量及季节性的电热综合能源系统优化调度

实现平台：matlab-yalmip-cplex/gurobi

代码简介：代码主要做的是电热综合能源系统的优化调度问题，代码注重对设备容量以及季节性差异的算例测试对比，算例中，对比了储能容量、储热容量变化对调度结果的影响，同时对比了夏季以及冬季两个典型日下电热综合能源系统的优化调度结果。结果非常清晰，出图量超多，分析十分到位，效果非常清楚，加以完善即是一篇中核期刊水平。附带参考文献。

参考文献：《含电热联合系统的微电网运行优化》《碳交易机制下考虑需求响应的综合能源系统优化运行》

代码获取方式：【代码分享】考虑储能设备容量及季节性的电热综合能源系统优化调度

专题推荐：论文推荐，代码分享（点击即可跳转）

运行结果展示

%% 市场交易参数
% 购电电价
Crb = [0.17, 0.17, 0.17, 0.17, 0.17, 0.17, 0.17, 0.49, 0.49, 0.49, 0.83, 0.83, 0.83, 0.83, 0.83, 0.49, 0.49, 0.49, 0.83, 0.83, 0.83, 0.49, 0.49, 0.17];
%购电电价
Crs = [0.13, 0.13, 0.13, 0.13, 0.13, 0.13, 0.13, 0.38, 0.38, 0.38, 0.65, 0.65, 0.65, 0.65, 0.65, 0.38, 0.38, 0.38, 0.65, 0.65, 0.65, 0.38, 0.38, 0.13];
Pgrid_min = -60;%购电最大值
Pgrid_max = 60;%售电最大值

%% 燃气轮机
Pmt_min = 15;%燃气轮机出力最小值
Pmt_max = 65;%燃气轮机出力最大值
Rmt_down = -5;%燃气轮机向下爬坡率
Rmt_up = 10;%燃气轮机向上爬坡率
Cm_mt = 0.025;%燃气运行维护成本

eta_l = 0.15;%散热损失率
eta_h = 0.9;%烟气回收率
Coph = 1.2;%溴冷机制热系数

%燃气轮机发热效率
eta_mt = [0.281511485,0.279660704,0.274107223,0.281714324,0.285366323,0.282392617,0.278652065,0.273527594,0.265437288,0.263784876,0.253341649,0.240684731,0.225559863,0.219695875,0.23198922,0.233347706,0.236346234,0.259440614,0.273964246,0.283534984,0.278375758,0.271587717,0.262965241,0.252302695];

%% 电锅炉
Peb_min = 0;%电锅炉出力最小值
Peb_max = 50;%电锅炉出力最大值
Reb_down = -3;  %电锅炉向下爬坡率
Reb_up = 5;%电锅炉向上爬坡率
eta_ah = 0.85; %热电转换效率
Cm_eb = 0.016; %电锅炉运行维护成本
Cst_eb = 2.74;

%% 燃料电池
Pfc_min = 5;%燃料电池出力最小值
Pfc_max = 40;%燃料电池出力最大值
Rfc_down = -2;%燃料电池向下爬坡率
Rfc_up = 2;%燃料电池向上爬坡率
Cm_fc = 0.026;%燃料电池运行维护成本

% 燃料电池效率
eta_fc = [0.633102762,0.630810014,0.626210435,0.622515243,0.619496279,0.615849659,0.614475289,0.609895576,0.606003752,0.602201083,0.599696666,0.595793234,0.595413236,0.593021958,0.595233072,0.599833088,0.604433107,0.60902944,0.613626799,0.618226815,0.62282674,0.627426773,0.632026849,0.636626883];
Cst_fc = 1.2;

%% 热泵
HPmin=0;%热泵出力最小值
HPmax=20;%热泵出力最大值
HP=0.025;%维护成本
COP_HP=4.4;%制热系数

%% 电储能
tau_b = 0.001;% 自损耗系数
eta_bch = 0.9;%充电效率
eta_bdis = 0.9;%放电效率
Eb_init = 30;%初始蓄电量
Eb_min = 0.2*150;%蓄电量最小值
Eb_max = 0.8*150;%蓄电量最大值
Cm_Eb = 0.0018;%电池维护成本
Pbch_min=0;  %电储能充电功率最小值
Pbch_max=37.5;%电储能充电功率最大值
Pbdis_min=0;  % 电储能放电功率最小值
Pbdis_max=37.5;%电储能放电功率最大值

% 热储能
tau_h = 0.01;% 自损耗系数
eta_hch = 0.9;%蓄热效率
eta_hdis = 0.9;%放热效率
Eh_init = 0.2*100;%初始蓄热量
Eh_min = 0.1*100;%蓄热量最小值
Eh_max = 0.8*100/1.5;;%蓄热量最大值
Cm_Eh = 0.0016;%维护成本
Phch_min=0;    %热储能充热功率最小值
Phch_max=25;%热储能充热功率最大值
Phdis_min=0;   % 热储能放热功率最小值
Phdis_max=25;%热储能放热功率最大值

% 单价系数
Cgas = 2.5;%燃气单价
L_gas = 9.7;%天然气低热值
Uinit = 0;
Che = 0.1; % 制热收益系数

%环境成本参数
PP=2.53;%电量的折算系数
e_GT=0.3;%供电效率
h_GT=0.4;%供热效率
epei=0.57;%碳排放配额
epai=0.6101;%实际排放系数


%% 定义决策变量
% 输入参数
Ppv_1 = sdpvar(1,24);% 光伏实际出力
Pwt_1 = sdpvar(1,24);% 风电实际出力
% 电设备
Pmt = sdpvar(1,24);            % 燃气轮机电功率
Peb = sdpvar(1,24);            % 电锅炉电功率 
Pfc = sdpvar(1,24);              % 燃料电池电功率
Pgrid = sdpvar(1,24);             % 电网功率                   >0 买电，<0 卖电
Eb= sdpvar(1,24);               % 电储能的容量
Pbch = sdpvar(1,24);          % 电储能充电功率
Pbdis = sdpvar(1,24);         % 电储能放电功率
% 热设备
Pmth = sdpvar(1,24);          % 燃气轮机热功率
Pheb = sdpvar(1,24);          % 电锅炉热功率
Eh= sdpvar(1,24);               % 热储能的容量
Phch = sdpvar(1,24);          % 热储能储热功率
Phdis = sdpvar(1,24);         % 热储能放热功率
% 辅助变量
Ubch = binvar(1,24);   % 电池充电状态，1表示充电
Ubdis = binvar(1,24);   % 电池放电状态，1表示放电
Uhch = binvar(1,24);   % 热储能储热状态，1表示储热
Uhdis = binvar(1,24);   % 热储能放热状态，1表示放热

%修改-结构修改，加入热泵
P_HP=sdpvar(1,24,'full');%热泵输入功率