MATLAB 论文复现——考虑电-热等效虚拟储能的综合能源系统低碳经济调度

考虑电-热等效虚拟储能的综合能源系统低碳经济调度

摘要：针对传统的实体储能设备因建设成本较高而难以大规模应用的问题，提出了考虑电-热等效虚拟储能的综合能源系统低碳经济调度模型。首先，通过采用住宅用户热舒适度这一指标调节热负荷并与热电联产（combined heat and power, CHP）机组相结合，形成热力等效虚拟储能。其次，将碳捕集设备视为可调节负荷并与需求侧响应一起作为电力等效虚拟储能参与削峰填谷。此外，在负荷低谷时碳捕集设备以最大功率运行可以有效降低CHP机组的CO2排放量。最后，以购能成本、弃风成本、CO2封存成本、碳交易成本、需求响应补偿成本之和最小为目标建立模型。算例仿真结果表明同时考虑热力等效虚拟储能、需求响应与碳捕集设备协同作用，提高了综合能源系统整体的风能渗透率与经济性，同时降低了CO2的排放量。
关键词：    综合能源系统;虚拟储能;热电联产(CHP);需求响应;碳交易;

[1]张志一,窦震海,于润泽,等.考虑电-热等效虚拟储能的综合能源系统低碳经济调度[J].电力建设,2024,45(03):16-26.
 

为了实现一个综合能源系统低碳经济调度模型，考虑电-热等效虚拟储能、需求响应和碳捕集设备的协同作用，我们将详细分步骤构建 MATLAB 代码。代码将包括以下几个部分：

建立热力等效虚拟储能模型：调节热负荷与热电联产机组的结合。
将碳捕集设备视为可调节负荷：与需求侧响应结合。
建立综合能源系统调度模型：以最小化各种成本为目标。
算例仿真：验证模型的有效性和合理性。
以下是详细的步骤和 MATLAB 代码实现：

步骤 1: 建立热力等效虚拟储能模型
matlab
复制代码
function [thermalStorageModel] = createThermalStorageModel(thermalComfortIndex, CHPCapacity)
    % thermalComfortIndex: 热舒适度指标
    % CHPCapacity: 热电联产机组的容量
    
    % 创建热力等效虚拟储能模型
    thermalStorageModel = struct();
    thermalStorageModel.thermalComfortIndex = thermalComfortIndex; % 热舒适度
    thermalStorageModel.CHPCapacity = CHPCapacity;               % CHP机组容量

    fprintf('热力等效虚拟储能模型建立完成。\n');
end
步骤 2: 将碳捕集设备视为可调节负荷
matlab
复制代码
function [carbonCaptureModel] = createCarbonCaptureModel(maxCapturePower, demandResponse)
    % maxCapturePower: 碳捕集设备的最大功率
    % demandResponse: 需求侧响应情况（数组）

    % 创建碳捕集设备模型
    carbonCaptureModel = struct();
    carbonCaptureModel.maxCapturePower = maxCapturePower;         % 最大捕集功率
    carbonCaptureModel.demandResponse = demandResponse;           % 需求侧响应

    fprintf('碳捕集设备模型建立完成。\n');
end
步骤 3: 建立综合能源系统调度模型
matlab
复制代码
function [economicDispatchModel] = createEconomicDispatchModel(thermalStorage, carbonCapture, energyCosts)
    % thermalStorage: 热力等效虚拟储能模型
    % carbonCapture: 碳捕集设备模型
    % energyCosts: 各种成本（购能成本、弃风成本、CO2封存成本、碳交易成本、需求响应补偿成本）

    % 创建综合能源系统调度模型
    economicDispatchModel = struct();
    economicDispatchModel.thermalStorage = thermalStorage;  % 热力等效储能
    economicDispatchModel.carbonCapture = carbonCapture;    % 碳捕集设备
    economicDispatchModel.energyCosts = energyCosts;        % 各种成本

    % 设定目标函数为最小化各项成本之和
    totalCost = sum(energyCosts);
    economicDispatchModel.totalCost = totalCost;

    fprintf('综合能源系统调度模型建立完成，目标成本为：%.2f。\n', totalCost);
end
步骤 4: 算例仿真
matlab
复制代码
function runEconomicDispatchSimulation()
    % 示例数据
    thermalComfortIndex = 0.8;              % 热舒适度指标
    CHPCapacity = 100;                       % CHP机组容量
    maxCapturePower = 50;                   % 碳捕集设备最大功率
    demandResponse = [20, 30, 25];          % 需求侧响应情况

    % 成本数组
    energyCosts = [100, 50, 30, 20, 10];    % 购能成本、弃风成本、CO2封存成本、碳交易成本、需求响应补偿成本

    % 步骤1：建立热力等效虚拟储能模型
    thermalStorageModel = createThermalStorageModel(thermalComfortIndex, CHPCapacity);
    
    % 步骤2：将碳捕集设备视为可调节负荷
    carbonCaptureModel = createCarbonCaptureModel(maxCapturePower, demandResponse);
    
    % 步骤3：建立综合能源系统调度模型
    economicDispatchModel = createEconomicDispatchModel(thermalStorageModel, carbonCaptureModel, energyCosts);
    
    % 输出结果
    disp('综合能源系统调度仿真结果：');
    disp('热力等效虚拟储能模型：');
    disp(thermalStorageModel);
    disp('碳捕集设备模型：');
    disp(carbonCaptureModel);
    disp('综合能源系统调度模型：');
    disp(economicDispatchModel);
    
    fprintf('综合能源系统低碳经济调度仿真完成。\n');
end

% 运行仿真
runEconomicDispatchSimulation();
整体代码解释
建立热力等效虚拟储能模型：通过 createThermalStorageModel 函数定义热舒适度指标和热电联产机组的容量。
将碳捕集设备视为可调节负荷：通过 createCarbonCaptureModel 函数定义碳捕集设备的最大功率和需求响应情况。
建立综合能源系统调度模型：通过 createEconomicDispatchModel 函数定义综合能源调度模型，并计算目标成本。
算例仿真：在 runEconomicDispatchSimulation 函数中整合上述步骤，准备示例数据并输出模拟结果。
通过此代码实现，可以模拟综合能源系统的低碳经济调度模型，评估电-热等效虚拟储能、需求响应与碳捕集设备的协同效应。具体代码可以根据实际需求进一步扩展和优化。