02含电热联合系统的微电网运行优化

说明书

MATLAB代码：含电热联合系统的微电网运行优化

关键词：微网 电热联合系统 优化调度

参考文档：《含电热联合系统的微电网运行优化》完全复现

仿真平台：MATLAB yalmip+cplex

主要内容：提出基于电热联合调度的区域并网型微电网运行优化模型。综合网内储能特性、分时电价、电热负荷与分布式电源的时序特征，以包含风机、光伏电池、热电联产系统、电锅炉、 燃料电池和储能系统的并网型微电网为例， 采用 Cplex优化软件求得调度周期内各微电源最佳出力及总运行成本，并与两种常见电热调度方式进行比较。仿真算例表明：联合调度模型能实现电热统一协调调度并降低微电网运行成本。该模型可为电热之间能源互联及规划运营提供参考。

     

%% 参数
% 风光热电
Ppv = [0, 0, 0, 0, 0, 0, 2.52, 10.08, 16.38, 25.2, 28.98, 30.24, 28.98, 26.46, 21.42, 16.38, 11.34, 2.52, 0 , 0, 0, 0, 0, 0];
Pwt = [32.76, 37.8, 41.57, 39.06, 35.28, 15.12, 17.64, 21.42, 25.2, 15.12, 7.56, 11.34, 6.3, 5.04, 17.64, 21.42, 34.02, 34.02, 22.68, 18.9, 20.16, 20.16, 30.24, 39.06];
Phe = [128.5, 132.28, 139.84, 129.76, 134.8, 127.24, 123.47, 119.69, 119.69, 112.13, 100.79, 89.45, 79.37, 66.77, 76.85, 85.67, 104.57, 114.65, 127.24, 146.14, 141.1, 134.8, 131.02, 129.76];
Pl = [62.99, 57.95, 55.43, 51.65, 50.39, 62.99, 81.89, 97.01, 117.17, 127.24, 123.47, 108.35, 99.53, 93.23, 97.01, 109.61, 118.4, 122.21, 129.76, 136.06, 125.98, 105.83, 89.45, 66.77];
Cm_pv = 0.0235;
Cm_wt = 0.0196;

% 电网买卖电价
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];
Cst_mt = 1.94;

% 电锅炉
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;

% 电储能
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;
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;          % 制热收益系数