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

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

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

仿真平台：MATLAB yalmip+cplex

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

     

%% 定义决策变量
% 输入参数
parameter_cplex;
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表示放热
% 初始化
objective=0;
constraint=[];

%% 写个风光区间约束
constraint=[constraint,0.8*Ppv<=Ppv_1<=1*Ppv];
constraint=[constraint,0.8*Pwt<=Pwt_1<=1*Pwt];

%% 燃气轮机模型
constraint=[constraint, Pmth == ((Pmt.*(1 - eta_mt - eta_l))./eta_mt) * eta_h * Coph];
%% 电锅炉模型
constraint=[constraint, Pheb == Peb * eta_ah];