楼宇型综合能源服务系统智能优化运行研究

张爱平，赵利兴，刘静

（中国华电科工集团有限公司，北京 100070）

0 引言

随着“碳达峰、碳中和”目标的提出，综合能源服务项目的建设由小范围试点向全国范围推广。综合能源系统是与能源服务升级、多种能源互补发展、可再生能源大规模应用和综合能效提升等相伴而生的。综合能源服务面向能源系统终端，以用户需求为导向，通过集成多种能源系统，提高能源的互补利用，使用户用能满意度得到提升。综合能源服务旨在提供符合能源发展方向、契合用户实际需求、提升综合能源利用率的能源系统解决方案，是推动能源革命的重要举措［1-4］。

欧洲最早提出综合能源系统概念并付诸实施，早在欧盟第五框架（FP5）中就提及［5］，如分布式发电、运输和能源（Distributed Generation Transport and Energy，DGTREN）项目、Energie 项目及Microgrid 项目对综合能源系统进行研究，侧重能源协调、协同优化及靠近用户侧开发。美国注重综合能源智能管理相关理论技术的研发，美国能源部2001 年即提出了综合能源系统、智能管理系统的建设，注重各类能源系统间的协调配合。日本的综合能源管理系统则致力于引导能源结构优化，促进可再生能源规模化开发，这与日本能源严重依赖进口的现状密不可分。

我国综合能源服务尚处于起步阶段，冷、热、电、气各自独立供应，缺乏整体协调管理，存在能效低、运营成本高等问题。综合能源服务系统贴近用户侧，供能负荷需要实时与用户需求相匹配，而用户负荷变化大，导致运行方式复杂［6-10］；同时，综合能源系统机组耦合性强，边际效益条件影响因素多［11-16］。因此，需采用一套智能化控制系统，用高度集成的智能化控制设备取代人工操作，快速响应用户负荷需求变化，保障系统安全、经济运行。

楼宇型综合能源服务系统服务于园区内建筑楼宇智慧用能需求，楼宇用能具有季度性、时段性等特点。优化运行系统要实现供能与楼宇侧用能一体化数据通信，协调优化控制以及源、网、荷友好互动等功能。本文以某楼宇型综合能源服务示范项目为例，对综合能源服务系统优化运行控制策略进行研究，为类似项目提供参考。

1 系统概述

示范项目位于北京市某科技园区，为AB座、CD座、EFG 座、HI 座4 栋独立楼座供冷、热、电、生活热水，建筑规模为：地上约17 万m2，地下约8 万m2。综合能源供能系统包括天然气冷热电三联供系统、太阳能生活热水系统、光储系统及低温烟气余热回收系统。

其中天然气冷热电三联供系统建设2 台3.300 MW 级内燃机+2 台烟气热水型余热溴化锂机组（制冷2×2.908 MW，供热2×1.257 MW（不含缸套水热量））+2台离心式电制冷机（供冷2×1.800 MW）+1台直燃机（供冷3.490 MW，供热2.800 MW）+1 台直燃机（供冷3.490 MW，供热2.800 MW，供生活热水1.100 MW），为园区供冷、热、电及生活热水。太阳能生活热水集热器位于H 座屋顶，共44 组，作为酒店生活热水的补充供应源。光伏装机容量为187.22 kWp，储能装机为100 kW∕240 kW·h，光伏发电直接接入大楼低压配电系统。低温烟气余热回收系统为2 套内燃机组+1 台热泵机组及其辅助设备，每套内燃机组烟道设置1套烟气板式换热器。

园区综合能源服务示范项目自2020 年开始策划实施，计划2021年年底完成建设。通过示范项目的建设，实现综合能源服务系统的优化运行。

2 研究目标

楼宇型综合能源服务系统通过采集园区楼宇侧、管网侧数据建立负荷预测模型，利用大数据分析历史用能数据，找出负荷随季节和时间的变化规律；同时，结合天气预报信息，通过已建立的负荷预测模型，给出冷、热、电、生活热水24 h 负荷预测曲线并根据实时负荷不断进行修正。

在园区能源需求、供能设备运行、外购市电及天然气价格等边界条件的约束下，通过建立多目标函数，实现负荷分配寻优，给出设备组合运行策略。

3 优化运行技术方案设计原则

优化运行技术方案决定了示范项目的运行模式以及各种运行模式下的优化运行策略，设置以下设计原则：（1）天然气供应不受限制，满足机组不同负荷工况下天然气的需求；（2）园区电采用并网上网模式，即园区多余电可以上网，不足电力可以从电网购买市电；（3）在满足园区负荷需求的条件下，尽可能提高可再生能源利用率；（4）提高综合能源利用，余热要尽可能利用；（5）提高园区供能系统运行的经济性；（6）结合项目设备特点，限制部分设备的最低负荷。

4 联供系统介绍

联供设备主要由动力设备、余热利用设备、调峰设备、多能互补系统及相关辅助设备构成。本项目动力设备采用2套燃气内燃发电机组，1台内燃机对应1台烟气热水型溴化锂机组（余热机）。内燃机的缸套水和烟气余热经三通阀连接至余热机，通过调节缸套水、烟气三通阀的开度，可分别调节余热机缸套水及烟气的送入量，即可调节余热机的负荷。2 台余热机冬季可供空调热水，夏季可供空调冷水；同时，余热机集成了生活热水板式换热器，充分吸收烟气的余热，冬季设置低温烟气深度利用系统，利用低温烟气驱动热泵，制取空调热水。H座屋顶太阳能生活热水系统充分利用可再生能源，供给酒店生活热水系统。AB 座屋顶布置的分布式光伏系统为示范光储系统，配置储能系统。调峰设备配置2 台电制冷机组、2 台直燃机组，2 台直燃机组可供空调冷热水，其中一台机组还可供生活热水。系统配置了集成化板式换热器，当室外温度低于5 ℃时，获得空调冷水，对数据机房供冷。

示范项目联供系统结构如图1所示。

图1 联供系统结构Fig.1 Structure of the cogeneration system

5 负荷预测模型设计

针对不同的负荷，结合具体负荷特点进行负荷预测；同时，根据不同负荷的运行时间进行分段和分类预测。

负荷预测模型主要计算步骤如下：（1）获取冷、热、电、生活热水负荷历史记录；（2）分析冷、热、电、生活热水负荷历史记录；（3）确定趋势因子（根据历史记录做出趋势曲线）；（4）确定日内（工作时间段）和季节性方案；（5）预测负荷需求；（6）计算负荷预测误差以及负荷与季节的相关性；（7）将负荷预测转换至计划期段；（8）模型不断修正验算，直至满足预期结果。

负荷预测主要用于预测各类负荷在未来一段时间的变化情况，以便合理安排机组运行方式和机组的负荷调整方式，优化、控制负荷。负荷预测子系统的预测对象有冷、热、电、生活热水负荷，预测内容包括各类负荷的峰谷时刻和持续时间、典型日24 h负荷、短时内的精确负荷预测等，具体负荷预测流程如图2所示。

图2 负荷预测流程Fig.2 Process of load forecasting

5.1 电负荷预测

电负荷预测采用基于线性外推法的短期负荷预测方法，该方法的数学模型为

5.2 冷热负荷预测

冷热负荷要结合楼宇型综合能源的特性来考虑。影响冷热负荷的参数可以分为2 类：第1 类是自由属性参数，包括建筑的朝向、墙体和窗户的热阻、窗墙比等，建筑建造好后这类参数几乎不变；第2 类是运行属性参数，包括建筑的工作时间段、人员安排、人体舒适度要求、内部热源（如照明等用电设备）等。

建立夏冬季单位面积平均冷热负荷数学模型

式中：Qc，Qh为夏季和冬季单位面积平均冷、热负荷，W∕m2；Qc，ba为冷负荷基准值，按6—8 月制冷需求计算，实际基准值可按具体时段进行调整；Tc为工作时段制冷量比例系数之和；θ为建筑正南方向与地理正南方向的夹角；Gc(θ)为制冷季建筑的朝向对制冷的影响；Sc(θ)为制冷时间段辐射系数之和；Sc(t)为制冷时间段温差系数之和；Gc(r)为窗墙比对冷负荷的影响系数；Gc(Rwa)为外墙热阻对冷负荷的影响系数；Gc(Rwi)为窗户热阻对冷负荷的影响系数；Qh，ba为热负荷基准值，按12—2 月热需求计算，实际基准值可按具体时段进行调整；Th为工作时段制热量比例系数之和；Gh(θ)为制冷季建筑朝向对制热的影响；Sh(θ)为制热时段辐射系数之和；Sh(t)为制热时段温差系数之和；Gh(r)为窗墙比对热负荷的影响系数；Gh(Rwa)为外墙热阻对热负荷的影响系数；Gh(Rwi)为窗户热阻对热负荷的影响系数；k1为辐射转换系数；k2，k3为墙体冷、热传导系数。

6 智能优化运行方案设计

6.1 智能优化系统功能

智能优化系统负责信息的采集、通信及优化策略的寻优与显示。系统采集综合能源、空调管网、用户侧的数据，设有负荷预测、机组性能计算、负荷优化指令分配、机组状态监测、机组特性曲线在线学习等功能模块。智能优化系统基本策略是通过采集源、网、荷侧信息，预测未来24 h 负荷需求，完成机组特性曲线的学习、机组性能计算，进一步实现机组负荷的最优分配并下达各负荷指令和设定参数到相关设备。

具体运行优化策略如下：智能优化系统采集综合能源源侧分散控制系统（DCS）和设备信息数据，网络服务器采集空调管网数据、用户用能数据、天气数据等负荷相关数据，智能优化系统通过智能学习算法，完成整个综合能源系统主要设备运行特性的学习；采用智能寻优算法完成机组最优运行模式分析；通过优化控制器向DCS 发送优化控制指令（包括负荷优化分配指令、供热∕供冷方式优化与调控指令），调整机组的运行模式，通过网络服务系统发布综合能源系统及空调管网的最佳运行优化策略给运行值班人员，从而实现综合能源系统的整体优化控制，达到高效低耗、经济匹配、可靠安全的运行目的。

6.2 优化技术路线

在考虑能源市场、能源资源供需平衡和设备可利用性等因素的基础上确定技术路线，建立混合整数 线 性 规 划（Mixed-Integer Linear Programming，MILP）模型，如图3 所示。优化内容包括最优设备组合、最优设备负荷、最优运行策略，以经济成本最小为第1 优化目标，节能和环境指数为第2 优化目标。

图3 综合能源系统智能优化技术路线Fig.3 Roadmap of the intelligent optimization technology for the integrated energy system

6.3 确立优化目标函数

优化可分为单目标优化和多目标优化，综合能源系统中能源的协同优化多为多目标优化，求解多目标优化问题的解法较多，基于分解的多目标进化算法就是一种常用解法。相比于一般的数学优化方法，基于分解的多目标进化算法在解决非线性整数规划问题以及多维目标优化问题上具备很大的优势。一般情况下，多目标优化问题可描述为

式中：X∈Ω，f∈Θ，实现从n维向量到m维向量的转换，如图4所示；Ω，Θ为n及m维向量域。

图4 决策向量转换Fig.4 Transformation of the decision vector

6.4 智能优化运行计算模型

6.4.1 寻优目标函数

（1）全年总运行成本Csys_year，目标是Csys_year最小。

式中：Csys_year为综合能源系统的全年总运行成本；分别为外部能源供应成本、能源转换设备成本、储能总运行成本；Tan为该典型日在每季度所占的天数；n表示第n个季度；l为能源转换设备；j为储能设备。

（2）一天24 h 时序优化目标函数，目标是使供能成本最小，尽可能利用可再生能源循环制取生活热水、光伏发电。

式中：Csys_day为一天24 h供能成本。

6.4.2 优化约束条件

（1）能量平衡约束条件。

1）电力平衡。园区上下网电功率、光伏发电功率、内燃机组发电功率、储能放电功率之和，满足厂用电负荷、园区用电负荷及充电负荷之和。

式中：Pe_net(t，a)为电网上下网电功率，下网为正，上网为负；Ppv(t，a)为园区光伏发电功率；Pmt(t，a)为园区内燃机发电功率；Pdis(t，a)为园区储能充电功率；Pe(t，a)为园区用电负荷；Pchar(t，a)为园区充电负荷，包含储能电池和充电桩；Pauxi(t，a)为厂用电负荷；a表示某天；t表示某时刻。

2）空调热水平衡。

式中：Qh，yr(t，a)为余热机组热负荷；Qh，bh(t，a)为板换热负荷；Qh，zr(t，a)为直燃机组热负荷；Qh(t，a)为园区热负荷需求。

3）空调冷冻水平衡。

式中：Qc，yr(t，a)为余热机组冷负荷；Qc，dz(t，a)为电制冷机组冷负荷；Qc，zr(t，a) 为直燃机组冷负荷；Qc(t，a)为园区冷负荷需求。

4）生活热水平衡。

式中：Qs_sun(t，a)为太阳能生活热水负荷；Qs_sbh(t，a)为生活热水水水板式换热器负荷；Qs_ybh(t，a)为生活热水烟气板式换热器负荷；Qs_h(t，a)为园区生活热水负荷需求。

（2）能源转换设备运行约束条件：每台设备在最低负荷和满负荷之间运行。

（3）能源价格约束：在优化决策过程中，考虑电力峰谷电价和燃气价格进行能源平衡。

6.4.3 优化寻优方法

粒子群算法（Particle Swarm Optimization，PSO）是求解优化问题的常用方法，它具有收敛速度快、计算量小、易于实现等优点，但也有早熟收敛、维数灾难、易陷入局部极值等缺点［17-20］。为避免PSO 的缺点，采用如下解决办法。

（1）调整PSO 的参数来平衡算法的全局探测和局部开采能力。

（2）设计不同类型的拓扑结构，改变粒子学习模式，从而提高种群的多样性。

（3）将PSO 和其他优化算法结合，形成混合PSO算法。

根据园区能源系统的特点，采用混合PSO 算法策略来寻优。算法的基本思想是，采取一种机制，当陷入局部最优或发生早熟时，能够跳出此局部区域并进入解空间的其他区域继续进行搜索，直到找到全局最优解。算法具体流程如图5所示。

图5 混合PSO算法流程Fig.5 Flow chart of the mixed PSO algorithm

引入种群的适应度方差σ2反映粒子的聚集程度：σ2越小，空间中的粒子越收缩于局部；σ2越大，空间中粒子的多样性越好。

式中：n为粒子个数；fi为第i个粒子的适应度；favg为前i个粒子适应度的平均值。

变异机制

式中：gbest(i)为局部最优值；μ为随机变量。

在算法运行过程中，找到最优值后，进行一个是否是局部最优的判断，如果是局部最优就随机改变这个粒子的位置，从而让粒子进入其他区域进行搜索。在搜索过程中，算法可能发现新的个体极值以及全局极值，然后再进行一次是否满足平衡约束的判断，最终找到全局最优解。

实时优化调度模型中，采集实时负荷值并通过优化计算得到各机组出力指令值，每次求解变量相对较少，此时采用混合PSO 算法能够较好地解决非线性少变量的优化问题。混合PSO 优化算法主要包括产生初始解、计算目标函数值、计算约束条件下的函数值、计算适应度值、筛选最优粒子、粒子向最优粒子移动等步骤。能源站空调水温度和发电机出力可调节，使用混合PSO 算法求解，以空调水温度和发电机功率作为决策变量，以经济性为核心目标函数，同时考虑环境温度、机组效率、水泵功耗与供回水温差，得到最优解决方案，即最优控制策略。

7 负荷预测及优化运行

基于园区冷、热、电、生活热水特点建立负荷预测模型，同时结合外界环境温度，绘制全天24 h 负荷预测曲线，如图6—8 所示。由图6—8 可以看出，负荷预测值与实际值的偏差控制在10%以内。

图6 电负荷预测曲线Fig.6 Electricity load forecasting curve

图7 空调水负荷预测曲线Fig.7 Forecasting curve of the load for air conditioning water load

图8 生活热水负荷预测曲线Fig.8 Hot water load forecasting curve

根据冷、热负荷的需求，参考负荷预测值，在满足所考核的能源利用效率的前提下，以实现能源站经济效益最优为目标，对各供能设备及冷热网进行动态负荷优化分配，给出动态优化分配策略，调整设备运行组合方式及负荷参数。初步优化结果显示：通过烟气余热深度利用，可节约直燃机供空调热水的能耗，高峰时段每小时可节约50 m3天然气；通过减少余电上网，可提高能源经济效益；供能初期，通过调整供能参数，实现大温差小流量供能，可节约天然气耗量及厂用电量。2021 年该园区采用优化运行方案后，与2019 年相比，年能源综合效率提升了6.51百分点，自耗电减少了2.24%，见表1。

表1 优化运行前、后数据对比Table 1 Comparison of the data before and after optimized operation

智能优化前期先开环运行，后期逐步闭环运行，直接自动调整、控制设备运行，减少人为干预。后期会进一步结合外界天气情况，优化冷塔运行，节约冷塔耗电量。

8 结论

本文结合楼宇综合能源服务系统的建筑用能特性，考虑季节及时段等特点，同时结合外界天气状况，建立了冷、热、电、生活热水负荷预测模型。源侧根据多能互补系统特点，考虑设备运行特性曲线、能源价格等边界条件，建立优化目标函数，优化寻优算法，实现运行优化方案的寻优。

目前，该优化系统还在不断自学习过程中，后期通过现场实际运行积累的数据，进一步优化调整负荷预测及优化运行算法，逐步细化每个设备的运行参数，为提高项目运行的经济性创造更大空间，以期能为类似项目的优化运行提供借鉴。