【顶刊|修正】多区域综合能源系统热网建模及系统运行优化【复现+延伸】

主要内容

部分代码

结果一览

下载链接

主要内容

该程序复现《多区域综合能源系统热网建模及系统运行优化》模型并进一步延伸，基于传热学的基本原理建立了区域热网能量传输通用模型，对热网热损方程线性化实现热网能量流建模，在冷热电联供系统运行优化模型基础上建立了含有热网的多区域 IES 优化混合整数线性规划模型，以从电网购电费用、向电网售电所得费用、燃气费、弃光成本和热网运行费用之和为目标，包括CCHP、热网和耦合关系约束条件，采用蒙特卡洛和kmeans聚类生成可再生能源出力随机优化场景，并将原确定性模型调整为日前确定、实时鲁棒优化模型。程序采用matlab+gurobi编程，本次对程序进行了修正优化，IES方向的同学可以下载学习！

本次主要对程序做了三方面的修正：

1.解决了求解速度慢的问题，通过优化程序代码，目前求解时间为5min左右。

2.修正了热功率失衡问题。

3.通过addpath命令应用解决了程序对文件夹依赖问题，程序下载后可直接运行，无需修改文件夹路径及名称。

部分代码

%上述模型为 0-1 混合整数线性规划模型，其中决策变量有：CCHP 中各设备出力、储能装置输入输出
% 电网购电、售电量、转换设备输入以及热网中各管段输入输出热功率。目前该模型已有成熟的求解算法
% 导入参数
%居民区 CCHP 1;
%商业区 CCHP 2
%工业区 CCHP 3
%办公区 CCHP 4
%--------------------------------------------------------------------------第一阶段
warning off
addpath('..\data')
addpath('..\子函数')
clc
clear
close all
Params = LoadParams();

NumOfScence = 3;                                              % 场景数
NumOfDistract = 4;                                             % 区域数

Params = WindSceneProduce(Params,NumOfScence,1000);           % 光伏不确定性场景


p = ones(1,NumOfScence) /NumOfScence;

Preal =  Params.CCHP1.P';                 % 真实场景概率
R = 4;               % 循环水泵的数量

StateParams = StateParamsInit(NumOfDistract);
% 主问题热网变量
StateTemData = cell(1,24);
for i = 1:24
StateTemData{i} = TemFlowRate(Params);
end


% 实时情况
DifferentScence = NowParamsInit(NumOfScence,NumOfDistract);
% 子问题热网变量
DifferentTemData = cell(1,NumOfScence);

for j = 1:NumOfScence
    a = cell(1,24);
    for i = 1:24
        a{i} = TemFlowRate(Params);
    end
    DifferentTemData{j} = a;
end
clear a i j