考虑阶梯式碳交易与供需灵活双响应的综合能源系统优化调度

1 主要内容

该程序方法复现《考虑阶梯式碳交易与供需灵活双响应的综合能源系统优化调度》，提出了供需灵活双响应机制，供应侧引入有机朗肯循环实现热电联产机组热电输出的灵活响应，需求侧在考虑电热气负荷均具备时间维度上需求响应的同时，提出了3种负荷之间具备可替代性;最后构建了以碳排放成本、购能成本、弃风成本、需求响应成本最小为目标的优化调度模型，运用cplex求解器进行求解。

• 目标函数

购能成本+碳交易成本+弃风成本+ DR 补偿成本：

• ​模型：

2 部分代码

%% 决策变量
P_GT_e=sdpvar(1,24); %GT的输出电功率
P_GT_h=sdpvar(1,24); %GT的输出热功率
P_GT_g=sdpvar(1,24); %GT消耗天然气功率
P_h_WHB=sdpvar(1,24); %GT供给WHB的热功率
P_h_ORC=sdpvar(1,24); %GT供给ORC的热功率
P_WHB_h=sdpvar(1,24); %WHB输出热功率
P_ORC_e=sdpvar(1,24); %ORC输出电功率
P_CHP_e=sdpvar(1,24); %chp的输出电功率
P_CHP_h=sdpvar(1,24); %chp的输出热功率
P_DG=sdpvar(1,24); %风电消纳功率
P_g_GB=sdpvar(1,24); %输入GB设备的天然气功率
P_GB_h=sdpvar(1,24); %GB设备输出的热功率
P_e_P2G=sdpvar(1,24); %P2G耗电功率
P_P2G_g=sdpvar(1,24); %P2G输出的天然气功率
%电热气储能
P_ES1_cha=sdpvar(1,24);P_ES2_cha=sdpvar(1,24);P_ES3_cha=sdpvar(1,24); %充放功率
P_ES1_dis=sdpvar(1,24);P_ES2_dis=sdpvar(1,24);P_ES3_dis=sdpvar(1,24);
S_1=sdpvar(1,24);S_2=sdpvar(1,24);S_3=sdpvar(1,24); %各储能的实时容量状态
B_ES1_cha=binvar(1,24);B_ES2_cha=binvar(1,24);B_ES3_cha=binvar(1,24); %充标志
B_ES1_dis=binvar(1,24);B_ES2_dis=binvar(1,24);B_ES3_dis=binvar(1,24); %放标志
P_e_buy=sdpvar(1,24); %购电功率
P_g_buy=sdpvar(1,24); %购气功率
%需求响应负荷
P_e_load=sdpvar(1,24); %经过需求响应后实际的电负荷
P_e_tran=sdpvar(1,24);  %参与DR的可转移电负荷
P_e_tranin=sdpvar(1,24); %可转移电负荷的转入功率
P_e_tranout=sdpvar(1,24); %可转移电负荷的转出功率
Zeta_etranin=binvar(1,24);%可转移电负荷的转入功率二进制变量
Zeta_etranout=binvar(1,24);%可转移电负荷的转出功率二进制变量
P_e_re=sdpvar(1,24);  %参与DR的可替代电负荷
P_e_rein=sdpvar(1,24); %可替代电负荷的转入功率
P_e_reout=sdpvar(1,24); %可替代电负荷的转出功率
Zeta_erein=binvar(1,24);%可替代电负荷的转入功率二进制变量
Zeta_ereout=binvar(1,24);%可替代电负荷的转出功率二进制变量
P_h_load=sdpvar(1,24); %经过需求响应后实际的热负荷
P_h_tran=sdpvar(1,24);  %参与DR的可转移热负荷
P_h_tranin=sdpvar(1,24); %可转移热负荷的转入功率
P_h_tranout=sdpvar(1,24); %可转移热负荷的转出功率
Zeta_htranin=binvar(1,24);%可转移热负荷的转入功率二进制变量
Zeta_htranout=binvar(1,24);%可转移热负荷的转出功率二进制变量
P_h_re=sdpvar(1,24);  %参与DR的可替代热负荷
P_h_rein=sdpvar(1,24); %可替代热负荷的转入功率
P_h_reout=sdpvar(1,24); %可替代热负荷的转出功率
Zeta_hrein=binvar(1,24);%可替代热负荷的转入功率二进制变量
Zeta_hreout=binvar(1,24);%可替代热负荷的转出功率二进制变量
P_g_load=sdpvar(1,24); %经过需求响应后实际的热负荷
P_g_tran=sdpvar(1,24);  %参与DR的可转移热负荷
P_g_tranin=sdpvar(1,24); %可转移热负荷的转入功率
P_g_tranout=sdpvar(1,24); %可转移热负荷的转出功率
Zeta_gtranin=binvar(1,24);%可转移热负荷的转入功率二进制变量
Zeta_gtranout=binvar(1,24);%可转移热负荷的转出功率二进制变量
P_g_re=sdpvar(1,24);  %参与DR的可替代热负荷
P_g_rein=sdpvar(1,24); %可替代热负荷的转入功率
P_g_reout=sdpvar(1,24); %可替代热负荷的转出功率
Zeta_grein=binvar(1,24);%可替代热负荷的转入功率二进制变量
Zeta_greout=binvar(1,24);%可替代热负荷的转出功率二进制变量

3 程序结果