主要内容
该程序复现《多区域综合能源系统热网建模及系统运行优化》模型并进一步延伸，基于传热学的基本原理建立了区域热网能量传输通用模型，对热网热损方程线性化实现热网能量流建模，在冷热电联供系统运行优化模型基础上建立了含有热网的多区域 IES 优化混合整数线性规划模型，以从电网购电费用、向电网售电所得费用、燃气费、弃光成本和热网运行费用之和为目标，包括CCHP、热网和耦合关系约束条件，采用蒙特卡洛和kmeans聚类生成可再生能源出力随机优化场景，并将原确定性模型调整为日前确定、实时鲁棒优化模型。程序采用matlab+gurobi编程，本次对程序进行了修正优化，IES方向的同学可以下载学习！
本次主要对程序做了三方面的修正：
1.解决了求解速度慢的问题，通过优化程序代码，目前求解时间为5min左右。
2.修正了热功率失衡问题。
3.通过addpath命令应用解决了程序对文件夹依赖问题，程序下载后可直接运行，无需修改文件夹路径及名称。
部分代码
%上述模型为 0-1 混合整数线性规划模型，其中决策变量有：CCHP 中各设备出力、储能装置输入输出% 电网购电、售电量、转换设备输入以及热网中各管段输入输出热功率。目前该模型已有成熟的求解算法% 导入参数%居民区 CCHP 1;%商业区 CCHP 2%工业区 CCHP 3%办公区 CCHP 4%--------------------------------------------------------------------------第一阶段warning offaddpath('..\data')addpath('..\子函数')clcclearclose allParams = 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:24StateTemData{i} = TemFlowRate(Params);end% 实时情况DifferentScence = NowParamsInit(NumOfScence,NumOfDistract);% 子问题热网变量DifferentTemData = cell(1,NumOfScence);for j = 1:NumOfScencea = cell(1,24);for i = 1:24a{i} = TemFlowRate(Params);endDifferentTemData{j} = a;endclear a i j
结果一览
主场景
因出图结果较多，其他场景出图结果不再展示，有兴趣的可以获取程序源码进一步了解。
“阅读原文”获取程序源码
