1. 项目背景与核心价值
冷热电多能互补综合能源系统(Integrated Energy System, IES)是当前能源领域的研究热点,它通过电、热、冷、气等多种能源形式的协同优化,实现能源的高效利用。传统IES优化调度往往只关注经济性指标,而忽略了用户体验这一关键因素。这套MATLAB代码的创新之处在于,首次将预测平均投票数(PMV)这一国际通用的热舒适度指标引入优化模型,使系统调度不仅考虑经济成本和碳排放,还能主动保障用户的环境舒适度。
PMV指标由丹麦学者P.O. Fanger提出,它综合考虑了环境温度、湿度、空气流速、平均辐射温度、人体活动强度和衣着情况六个因素,通过-3(冷)到+3(热)的7级标度量化人体热感受。在建筑能耗中,供暖和制冷占总能耗的40%-60%,而用户对温度的主观感受直接影响设备使用行为。通过PMV约束,系统可以动态调整供能策略,在舒适与节能之间找到最佳平衡点。
2. 系统架构与关键技术
2.1 整体架构设计
系统采用"源-网-荷-储"协同框架,主要包含以下组件:
- 能源生产层:光伏阵列、风力发电机、燃气轮机(CHP)
- 能源转换层:余热锅炉、吸收式制冷机、电制冷机、P2G(电转气)设备
- 能源存储层:隐含考虑了建筑本身的热惯性储能特性
- 负荷需求层:电负荷、热负荷、冷负荷、气负荷
各组件通过能量流网络连接,形成如图1所示的耦合关系。特别值得注意的是,系统通过热惯性模型将建筑结构作为虚拟储能设备,利用墙体、家具等的蓄热特性平抑负荷波动。
2.2 关键技术实现
2.2.1 用户舒适度建模
PMV值与室内温度的关系可表示为:
code复制PMV = (0.303e^(-0.036M) + 0.028)×[(M-W) - 3.05×10^-3(5733-6.99(M-W)-Pa) - 0.42((M-W)-58.15) - 1.7×10^-5M(5867-Pa) - 0.0014M(34-Ta) - 3.96×10^-8fcl((Tcl+273)^4-(Tr+273)^4) - fclhc(Tcl-Ta)]
其中M为代谢率,W为机械功,Pa为水蒸气分压,fcl为服装面积系数,Tcl为服装表面温度,hc为对流换热系数,Tr为平均辐射温度。
在代码实现中,通过固定其他参数(设定典型办公场景:代谢率1.2met,服装热阻1clo),将PMV简化为室内温度的函数,建立约束条件:
matlab复制% PMV约束示例
if PMV_set == 1 % 中等舒适要求
constraints = [constraints, 20 <= T_in <= 24];
elseif PMV_set == 2 % 较高舒适要求
constraints = [constraints, 22 <= T_in <= 23];
end
2.2.2 热惯性动态模型
建筑热动态采用RC网络模型,将墙体简化为热阻R和热容C的等效电路:
code复制C*dT/dt = (T_out - T_in)/R + Q_HVAC
离散化后得到代码中的递推公式:
matlab复制T_in(t) = T_in(t-1) + dt/(R*C)*(T_out(t) - T_in(t-1)) + Q_HVAC(t)/C
其中dt为时间步长(1小时),Q_HVAC为供热/冷功率。
2.2.3 多能流耦合约束
系统需要同时满足四种能量平衡:
- 电平衡:
code复制P_PV + P_WT + P_MT + P_net = P_load + P_EC/COP_EC + P_P2G/η_P2G - 热平衡:
code复制P_EH + P_GH = Q_load + P_AC/COP_AC - 冷平衡:
code复制P_AC + P_EC = Q_cool - 气平衡:
code复制G_buy + G_P2G = G_MT + G_GH + G_load
3. 优化模型构建
3.1 决策变量定义
使用YALMIP工具箱定义优化变量:
matlab复制% 连续变量
P_MT = sdpvar(24,1); % 燃气轮机出力
P_EH = sdpvar(24,1); % 余热锅炉出力
P_GH = sdpvar(24,1); % 燃气锅炉出力
P_AC = sdpvar(24,1); % 吸收式制冷机
P_EC = sdpvar(24,1); % 电制冷机
P_P2G = sdpvar(24,1); % P2G设备
P_buy = sdpvar(24,1); % 电网购电
P_sell = sdpvar(24,1);% 电网售电
% 二进制变量
Temp_net = binvar(24,1); % 购售电状态标志
3.2 目标函数
提供两种优化目标:
-
经济成本最小化:
matlab复制Cost_gas = 0.25 * sum(G_buy); % 天然气成本 Cost_elec = 0.7 * sum(P_buy) - 0.2 * sum(P_sell); % 电网交互成本 Objective = Cost_gas + Cost_elec; -
碳排放最小化:
matlab复制Emission_gas = 0.4483 * sum(G_buy); % 天然气碳排放 Emission_elec = 0.805 * sum(P_buy); % 电网购电碳排放 Objective = Emission_gas + Emission_elec;
3.3 约束条件
典型设备约束示例:
matlab复制% 燃气轮机约束
constraints = [constraints, 9 <= P_MT <= 18]; % 出力上下限
constraints = [constraints, P_EH == 0.68*0.6*P_MT]; % 余热回收关系
% 购售电逻辑约束
for t = 1:24
constraints = [constraints, P_buy(t) <= 1000*Temp_net(t)];
constraints = [constraints, P_sell(t) <= 1000*(1-Temp_net(t))];
end
4. 代码实现与使用指南
4.1 运行环境配置
- 安装MATLAB R2017b或更高版本
- 添加YALMIP工具箱(建议R20200930版)
- 安装IBM ILOG CPLEX求解器(学术版免费)
4.2 数据准备
示例数据文件shuju_data.xlsx应包含以下工作表:
- 电负荷:24小时数值,单位kW
- 气负荷:24小时数值,单位m³/h
- 光伏预测:预期出力,单位kW
- 风电预测:预期出力,单位kW
- 室外温度:每小时实测值,单位℃
4.3 典型运行流程
matlab复制% 主脚本示例
clc; clear; close all;
% 读取数据
data = readtable('shuju_data.xlsx');
% 设置PMV等级 (1或2)
PMV_level = 1;
% 选择优化目标 (1:经济, 2:碳排放)
opt_mode = 1;
if opt_mode == 1
[cost, dispatch] = economic_opt(data, PMV_level);
else
[emission, dispatch] = emission_opt(data, PMV_level);
end
% 结果可视化
plot_results(dispatch, data);
5. 应用案例与结果分析
5.1 不同PMV设置对比
设置两种PMV等级运行经济优化模型:
| 指标 | PMV=1(宽松) | PMV=2(严格) | 变化率 |
|---|---|---|---|
| 总成本(元) | 5824 | 6127 | +5.2% |
| 碳排放(kg) | 2865 | 3012 | +5.1% |
| 燃气轮机运行小时数 | 18 | 21 | +16.7% |
分析表明:舒适度要求每提高1个PMV等级,系统运行成本增加约5%,主要源于燃气轮机更频繁启停和更高出力的需求。
5.2 双目标优化对比
固定PMV=1时两种优化目标的结果:
| 设备 | 经济最优出力 | 低碳最优出力 | 差异 |
|---|---|---|---|
| 燃气轮机 | 142 kW | 98 kW | -31% |
| 光伏利用率 | 78% | 92% | +18% |
| 电网购电量 | 256 kWh | 312 kWh | +22% |
结果显示:低碳模式会抑制燃气轮机使用,优先消纳可再生能源,但可能增加电网依赖。
6. 工程实践建议
-
参数校准建议:
- 建筑热阻R和热容C应通过实测数据辨识
- PMV参数需根据实际人群特征调整(如老年人偏好更高温度)
- 设备效率参数应使用厂家提供的实测曲线
-
扩展应用方向:
- 增加电储能、储热罐等物理储能设备
- 考虑需求响应机制,引入电价弹性负荷
- 结合机器学习改进风光出力预测精度
-
常见问题排查:
-
问题1:求解器报"infeasible"
- 检查负荷数据是否超出设备容量总和
- 验证PMV约束的温度范围是否合理
-
问题2:结果出现不合理波动
- 确认时间步长是否与数据采样率匹配
- 检查热惯性模型参数单位是否一致
-
这套代码为综合能源系统设计提供了兼顾经济性、环保性和舒适性的决策工具,特别适合区域能源站、绿色建筑等场景的优化运行。通过调整PMV参数,运营商可以在用户体验和运营成本之间找到最佳平衡点。
