说明书

MATLAB代码：碳交易机制下考虑需求响应的综合能源系统优化运行

关键词：需求响应 碳交易机制 日前优化调度 综合能源系统

参考文档：《碳交易机制下考虑需求响应的综合能源系统优化运行》完美复现

仿真平台：MATLAB yalmip+cplex

主要内容：代码主要做的是碳交易机制下考虑需求响应的综合能源系统优化运行问题，根据负荷响应特性将需求响应分为价格型和替代型 2 类， 分别建立了基于价格弹性矩阵的价格型需求响应模型， 及考虑用能侧电能和热能相互转换的替代型需求响应模型; 其次， 采用基准线法为系统无偿分配碳排放配额，并考虑燃气轮机和燃气锅炉的实际碳排放量，构建一种面向综合能源系统的碳交易机制; 最后，以购能成本、碳交易成本及运维成本之和最小为目标函数，建立综合能源系统低碳优化运行模型，并通过 4 类典型场景对所提模型的有效性进行了验证（四类场景都有对应代码，都实现了）。

%% 碳交易机制下考虑需求响应的综合能源系统优化运行——魏震波
%场景 1: 仅考虑碳交易

clc;clear;close all;% 程序初始化
%% 读取数据
shuju=xlsread('carbon+DR数据.xlsx'); %把一天划分为24小时
load_e=shuju(2,:); %初始电负荷
load_h=shuju(3,:); %初始热负荷
P_PV=shuju(4,:);    %光电预测
P_WT=shuju(5,:);    %风电预测
pe_b=shuju(6,:); %需求响应前电价
pe_a=shuju(7,:); %需求响应电价
ph_b=shuju(8,:); %需求响应前热价
ph_a=shuju(9,:); %需求响应热价

% %% 需求侧定义变量
% Z=zeros(24,24); %需求弹性矩阵
% e_W1=0.5;e_W2=0.3;e_W3=0.15;e_W4=0.05;%约束：固定、可转移、可消减、可替代负荷占比50%，30%，15%，5% 
% h_W1=0.5;h_W2=0.2;h_W3=0.2;h_W4=0.1;%约束：固定、可转移、可消减、可替代负荷占比50%，30%，15%，5%  
% Psl_e=zeros(1,24);%转移电负荷量
% Pcl_e=zeros(1,24);%消减电负荷量
% Prl_e=zeros(1,24);%电负荷被替代量
% Psl_h=zeros(1,24);%转移热负荷量
% Pcl_h=zeros(1,24);%消减热负荷量
% Prl_h=zeros(1,24);%热负荷被替代量
% P2H=1.83; %电转热系数
OP_load_e=zeros(1,24);%优化后的电负荷
OP_load_h=zeros(1,24);%优化后的热负荷
OP_load_e=load_e; %只考虑碳交易，优化后的负荷即为初始电负荷
OP_load_h=load_h; %只考虑碳交易，优化后的负荷即为初始热负荷
%% IES电网交互电价
price_buy_grid=shuju(7,:); %向电网购电价
price_sell_grid=shuju(10,:); %向电网售电价