计及建筑热负荷弹性的综合能源系统调度方法
Optimal Dispatch of Integrated Energy System Based on Flexibility of Thermal Load
通讯作者: 邵常政,助理研究员;E-mail:cshao@cqu.edu.cn.
责任编辑: 孙伟
收稿日期: 2021-12-28 修回日期: 2022-02-18 接受日期: 2022-04-2
基金资助: |
|
Received: 2021-12-28 Revised: 2022-02-18 Accepted: 2022-04-2
作者简介 About authors
胡博(1983-),教授,博士生导师,从事电力系统可靠性评估及综合能源系统优化运行研究.
具有一定弹性的建筑热负荷被视为电-热综合能源系统运行优化的重要调节资源.考虑建筑热负荷具有规模大、单体容量小的特点,非侵入式的数据驱动方法成为量化建筑热负荷弹性的有效手段.然而,由于数据不足或模型精度不够,该方法将不可避免地产生误差,给电-热综合能源系统的优化调度带来认知不确定性.因此,提出一种考虑建筑热负荷弹性并兼容相关认知不确定性的电-热综合能源系统优化调度方法.分析基于数据驱动的建筑热负荷需求弹性评估方法,将评估过程中产生的误差建模为认知不确定性,并通过改进的D-S证据理论对多源误差进行融合;采用拉丁超立方抽样方法生成表征热负荷弹性认知不确定性的场景,并通过模糊聚类法进行场景削减;将构造的场景集嵌入电-热综合能源系统的协调优化调度中,实现对建筑热负荷弹性及相关认知不确定性的综合考虑.算例仿真结果表明,考虑建筑热负荷需求弹性及认知不确定性对减少弃风、提高电热综合能源系统的运行灵活性至关重要.
关键词:
The flexibility of thermal loads of buildings is a valuable balancing resource for operation of the heat and electricity integrated energy system (HE-IES). Considering the characteristics of large scale and small single load capacity of the themal load, the non-intrusive data-driven method has become an effective means to quantify the flexibility of building thermal load. However, due to the inaccuracy of the model or the lack of data, this method inevitably produces errors and brings epistemic uncertainty to the optimal dispatch of the HE-IES. An optimal dispatch model of the HE-IES that is compatible with the epistemic uncertainty of demand flexibility in the thermal loads of buildings is proposed. First, a data-driven flexible demand assessment method for building thermal load is described. The measurement errors are modeled as epistemic uncertainty and the multiple error sources are combined by using the D-S evidence theory. Then, the representative scenarios are selected to represent the epistemic uncertainty of the demand flexibility based Latin hypercube sampling(LHS) method, and the scenarios are reduced by the fuzzy clustering method. Finally, the representative scenarios are embedded in the coordinated and optimized dispatch of the HE-IES to realize the comprehensive consideration of the thermal load flexibility and related epistemic uncertainty of the building. The results demonstrate that considering the epistemic uncertainties of the thermal load demand is crucial for reducing the wind power curtailments and improving the operational flexibility of HE-IES.
Keywords:
本文引用格式
胡博, 程欣, 邵常政, 黄威, 孙悦, 谢开贵.
HU Bo, CHENG Xin, SHAO Changzheng, HUANG Wei, SUN Yue, XIE Kaigui.
为量化建筑热负荷的弹性,文献[9]中将建筑热负荷建模为“储能电池”,以电池的充放电代表建筑热负荷功率的调节过程.然而,这类建模方法往往忽略了建筑本身的传热特性,使模型本身失去建筑传热特性物理层面的可解释性.此外,文献[10-11]中提出采用实验方法对建筑传热特性及潜在灵活性进行研究.但工作量较大,限制了该方法的应用[12].考虑上述方法的缺点,文献[12-13]中提出利用数据驱动的方法量化建筑热负荷的弹性.文献[14⇓-16]中基于大规模监测数据和专家信息,构建量化热负荷弹性的关键参数与其影响因素的关系模型,进而实现对建筑热负荷弹性的预测.为提高供热系统的智能化程度,文献[17]中采用混合数据驱动模型对未来1 h的热负荷进行预测建模,文献[18]中建立了基于深度长短期记忆网络的综合能源多负荷联合预测模型.然而,数据驱动方法不可避免地存在数据不足和训练模型不精确的问题,随之带来热负荷弹性量化测量误差,这种误差可能会导致对建筑热负荷需求弹性的过度高估,从而对电-热综合能源系统的协调优化调度构成风险.
基于该背景,本文提出一种考虑建筑热负荷弹性并兼容相关认知不确定性的电-热综合能源系统协调优化调度方法.首先,分析基于数据驱动的建筑热负荷需求弹性评估方法,将评估过程中产生的误差建模为认知不确定性,并通过改进的D-S证据理论对多源误差进行融合,获得描述建筑热负荷需求弹性评估的综合性误差值.然后,采用拉丁超立方抽样方法和模糊聚类法生成代表性场景以表征热负荷弹性及认知不确定性.最后,实现计及建筑热负荷弹性以及相关认知不确定性的电-热综合能源系统协调优化调度.
1 建筑热负荷需求弹性的量化与表征
1.1 基于一阶等效热参数模型的热负荷需求建模
式中:
图1
1.2 基于数据驱动的建筑热负荷需求弹性评估
(1) 提取过程.基于智能监控设备,收集历史建筑热负荷数据和气象参数,并基于历史数据从中提取kp,r,t和Cp,r,t作为后续研究样本.
(2) 训练过程.基于数据驱动算法建立kp,r,t和Cp,r,t与气象参数之间的关系,并利用气象预测数据对kp,r,t和Cp,r,t进行预测,为建筑热负荷需求弹性评估提供指导.
然而,在上述提取和训练过程中,针对kp,r,t和Cp,r,t预测的测量误差不可避免.这些测量误差均可归结为认知不确定性引起的误差.认知不确定性是由于模型不准确或缺乏数据而导致对问题缺乏理解,进而产生的认知偏差[23].
在上述基于数据驱动的预测方法中,导致kp,r,t和Cp,r,t预测误差的认知不确定性主要来源于提取和训练两个过程.因此,预测kp,r,t和Cp,r,t时有必要考虑多源认知不确定性带来的后果.如图2所示,提出基于改进D-S证据理论的多源误差融合模型解决上述问题,其中
图2
图2
建筑热参数kp,r,t多源误差融合示意图
Fig.2
Illustration of multi-source measurement errors fusion of kp,r,t
2 多源误差认知不确定性建模
合理处理kp,r,t和Cp,r,t预测过程中的多源误差及误差相互之间的关系极其重要.基于改进的D-S证据理论,全面考虑认知不确定性导致的多源误差对kp,r,t和Cp,r,t预测的影响.表1为kp,r,t和Cp,r,t的多误差来源及各误差数量级分析.可知,kp,r,t和Cp,r,t不同来源误差所属数量级存在差异,而针对同一处理过程,kp,r,t和Cp,r,t对应的误差所属数量级也并不相同.这是因为kp,r,t和Cp,r,t对应的建筑热负荷体量存在差异,且受气象参数影响的程度不同,导致同一模型求解所得的kp,r,t和Cp,r,t预测误差存在差异.
表1 kp,r,t和Cp,r,t多源误差来源及数量级分析
Tab.1
参数 | 误差来源(误差数量级) |
---|---|
kp,r,t | 提取过程(10-3) |
训练过程(10-2) | |
Cp,r,t | 提取过程(10-6) |
训练过程(10-3) |
2.1 多源误差认知不确定性建模
设非空集合Θ为识别框架,定义映射m: 2Θ→ [0, 1]满足
式中:m为基本概率分配函数即mass函数,可以理解为对子集X的支持度或信任度.如果m(X)>0,则称X为焦元.
信任函数定义如下:
根据D-S证据理论的Dempster融合规则,融合后焦元X的mass函数如下:
式中:A=X1∩X2∩…∩Xn;m1, m2, …, mn 为识别框架Θ中不同证人针对n个焦元X提供的mass函数,X1, X2, …, Xn为对应焦元.
由表1可知,在提取和训练的过程中,kp,r,t和Cp,r,t平均相对误差的数量级具有差异性.根据传统证据理论,将kp,r,t和Cp,r,t的多源误差视作事件集,当各误差出现的概率表征mass函数时,由于kp,r,t和Cp,r,t多源误差所属不同数量级,将导致构造的多源误差对应的事件集无重合区,即不同证人只对指定区域的事件提供专家信息.然而,每位证人只有知悉全局,并给出整个事件集的支持度,才可以提供完整准确的专家信息.因此,上述做法不合逻辑,必然会影响证据融合的结果.
针对以上问题,将认知不确定性导致的kp,r,t和Cp,r,t多源测量误差视作是对同一事件集的多方面证据描述.基于改进的D-S证据理论,多源误差相融合,得到一个描述kp,r,t和Cp,r,t预测误差的综合性结果,具体改进如下.
(1) 构建识别框架:
式中:Xd为kp,r,t和Cp,r,t多源相对误差对应的同一等效事件;n为等效事件的数量.
(2) 构建mass函数.kp,r,t和Cp,r,t的测量误差主要来自两个方面,即对于同一事件Xd,存在两个证人提供对应的mass函数值,直接将多源误差视作多位证人对同一等效事件的信任度,具体关系如下:
利用式(7)对kp,r,t多源相对误差进行融合,得到m(Xd)便可表征kp,r,t的综合性相对误差.本文采用
图3所示用以说明算法改进前后各方面性能的提升.可知,基于上述方法改进后,在处理kp,r,t和Cp,r,t多源误差融合时,等效事件集数量的减少使得计算负担同步下降;同时各位证人对整个事件集提供专家信息,保证了专家信息的完整性和有效性.
图3
2.2 基于信息熵融合的改进证据理论
识别框架Θ的mass函数对所有命题都提供支持度或称信任度.Pignistic概率转换方法转换的目的是对系统已获得的各命题信任值进行重新分配,以获得更可靠的决策依据,具体定义如下:
式中:θi为焦元Xd包含的信息;|Xd|为焦元Xd中的元素个数.
当信任函数进行信息转移时,子命题的支持度也随之变化,为此定义条件支持度的概念度量信任函数信息转移过程中对子命题的支持度的变化,即:
式中:θi∈H, H⊆Θ, ∀θj(θj∈H); SH(θj|θi)和LH(θj|θi)分别为mass函数在集合H上提取信息θi对θj的条件支持度和剩余支持度,并且可以进一步得到以下关系:
式(15)说明mass函数在H上对子命题的条件支持度是相互的.式(16)表示当对子命题的无条件支持度一定时,子命题的条件支持度越大,其剩余支持度越小.从式(17)可知,mass函数在H上目标子命题对子命题的剩余支持度即为mass函数在H上提取目标子命题信息后对子命题的无条件支持度.
采用一种新的信任函数信息转移算法,该算法基于Pignistic转换,逐步提取子命题的条件支持度,并将其重点放在无条件支持度大的子命题上,具体算法如下:
(1) 输入事件集Θ={X1,X2,…,Xd}和对应的mass函数值{m(X1), m(X2), …, m(Xd)}.
(2) 由式(18)计算每一个自命题的无条件支持度,即
(3) 令pi=SΘ(θi),得到分布(p1, p2, …, pn)为信任函数的初始状态.
(4) 通过式(13)和(14)计算条件支持度与剩余支持度.
(5) 选择信息融合下条件支持度最大的子命题,重新计算每个原子命题的条件支持度,直到所有信息融合完毕.
(6) 修正后的(p1, p2, …, pn)被认为是识别框架Θ的新mass函数值.
基于上述讨论即可获得修正后描述kp,r,t 和 Cp,r,t综合性相对误差的mass函数值.
3 基于代表性场景的电热综合能源系统协调调度模型
3.1 计及建筑热负荷需求弹性及认知不确定性的代表性场景生成
重点介绍计及热负荷需求弹性及认知不确定性的场景生成方法,代表日的选取方法参考文献[28].
3.1.1 考虑kp,r, t和Cp,r,t认知不确定性的场景生成
建筑热负荷需求弹性主要由kp,r,t和Cp,r,t刻画,因此通过抽样kp,r,t和Cp,r,t的误差场景表示各自相对误差的不确定性.目前,蒙特卡罗(MC)抽样和拉丁超立方(LHS)抽样是使用较为广泛的抽样方法[29].MC抽样技术是完全随机的,在输入分布的范围内,样本可以落在任何位置,但样本更有可能从高发生概率的分布区域中抽取.而LHS抽样根据输入变量的累积分布函数(CDF)生成随机样本数[30],它从随机变量的整个分布中进行抽样,可以确保小概率事件在模拟的输出中被准确地代表,从而提供更稳定和更精确的估计结果.同时LHS抽样迭代次数少,在抽样效率和运行时间方面的表现均更为优异.因此,采用LHS抽样方法模拟kp,r,t和Cp,r,t的不确定性,应用过程如下:
(1) 确定抽样样本数Ns,将[0,1]分成Ns层,并在每层中随机产生一个随机数.
(2) 根据样本数NLHS,设置NLHS个区间.
(3) 判断步骤(1)中抽取的随机数所处NLHS个区间中的位置.
(4) 若随机数处于第q个区间,则该随机数对应的kp,r,t和Cp,r,t的相对误差可在此区间上利用线性法计算,由此产生Ns个描述kp,r,t和Cp,r,t不确定性的场景.
3.1.2 基于模糊聚类算法的场景缩减技术
LHS抽样方法虽然在一定程度上可以提高建筑热负荷需求弹性不确定性评估的准确性,但也为后续的调度过程带来巨大的计算负担.对此,使用模糊聚类算法进行场景缩减,并生成代表性场景来描述建筑热负荷需求弹性及不确定性.
场景缩减技术根据场景的相似程度将场景分为不同的集群.模糊控制是自动控制领域的经典方法,模糊集引入了隶属度的概念,不遵循经典数学中“非0即1”的约定,而是提供了更灵活的聚类结果[31].
模糊聚类算法的目标函数是最小化聚类场景与聚类簇中聚类中心之间的欧氏距离,具体如下:
式中:U 为z×g隶属度矩阵,z为设定的聚类类别数,g为待选场景的数量;W为聚类中心矩阵;G为候选聚类场景的一维矩阵;uvw 为场景xv 对聚类 w 的隶属度;xw为待聚类的第w个场景;dvw 为xw和av之间的欧氏距离;av为每个聚类的聚类中心;模糊系数q'一般由经验决定,q'值过小或过大均会导致聚类算法的效果变差.各参数需满足的约束条件及算法具体流程见文献[32].
基于上述模型即可生成描述建筑热负荷需求弹性及认知不确定性的代表性场景,并应用于电热协调调度模型中.
3.2 电热综合能源系统协调调度模型
基于上述方法,综合考虑建筑热负荷需求弹性及认知不确定性,将建筑热负荷的灵活需求响应模型纳入电热综合能源系统优化调度模型中,以解决促进风电消纳等能源问题.
融合建筑热负荷灵活需求响应模型的电-热协调调度模型简述如下:
式中:
基于上述讨论,提出计及建筑热负荷需求弹性和认知不确定性的电-热综合能源系统协调调度平台如图4所示,其主要分为以下4个子模块,数据交汇中心为各个子模块间数据传递提供桥梁.
图4
图4
建筑热负荷需求弹性及认知不确定性处理方法及应用
Fig.4
Processing method and application of flexible demand and epistemic uncertainty of building thermal loads
(1) 数据预处理模块.基于历史热负荷数据和气象参数,提取传热参数k和C,并基于数据驱动训练算法,建立k、C与气象参数间的关系,从而实现未来气象条件下k和C的预测.
(2) 误差融合模块.将提取和训练过程中产生的k、C多源测量误差建模为认知不确定性,并采用改进的证据理论将多源测量误差融合以获得描述k和C不确定性的综合性误差值.
(3) 场景生成模块.采用LHS抽样和模糊聚类法生成描述建筑热负荷需求弹性及认知不确定性的代表性场景,用于后续的电热协调调度模块.
(4) 电-热综合能源系统协调调度模块.建筑热负荷的灵活需求响应模型成功嵌入电-热综合能源系统优化调度中,实现计及建筑热负荷需求弹性及认知不确定性的电-热综合能源系统协调调度.
4 算例仿真
构建算例仿真以验证本文方法的准确性和有效性.首先,讨论用于评估建筑热负荷需求弹性及认知不确定性方法的有效性.然后,给出考虑kp,t和Cp,t认知不确定性前后建筑热负荷评估结果.此外,使用简单的测试系统与已有调度方法对比,证明融合认知不确定性模型后电热综合能源系统调度灵活性的提升.模型仿真均在 MATLAB R2018b环境下进行,线性优化部分由Gurobi求解器求解.
4.1 kp,t和Cp,t的多源误差融合
基于智能监控设备获取300组历史热负荷及天气参数数据集,从中提取kp,t和Cp,t并求取对应的认知不确定性导致的偏差,同时给出多源误差的融合结果.
图5
图6
图7
图7
多源误差融合后Cp,t的综合性相对误差
Fig.7
Relative error of Cp,t after considering error fusion
图8
图8
多源误差融合后kp,t的综合性相对误差
Fig.8
Relative error of kp,t after considering error fusion
4.2 多源误差融合的精确性分析
设置两个子算例Case 1和Case 2,通过比较是否考虑kp,t和Cp,t多源误差融合时建筑热负荷的估计值,详细说明考虑多源误差融合的必要性,验证认知不确定性处理方法的准确性.其中,Case 1的训练误差作为传热系数kp,t和Cp,t的估计误差评估建筑热负荷的需求弹性;Case 2在评价建筑热负荷需求弹性时,考虑kp,t和Cp,t多源误差融合.
图9
图10
4.3 应用
为了说明建筑热负荷的需求弹性可以为系统提供更大的调节空间,将提出的建筑热负荷灵活响应模型应用于实际的电热综合能源系统中.系统结构及参数主要来自文献[34],并在此基础上,增加一台功率为600 MW的风力涡轮机和一台1000 kW的燃气发电机组,电能通过配电网络中每条馈线处的11/0.433 kV 变压器提供给负载;在区域供热系统中,设置3个节点提供热能.
使用LHS抽样生成更多描述建筑热负荷需求弹性及认知不确定性的场景.为了减小计算量,通过模糊聚类技术进行场景缩减.图11为具有代表性的弃风场景结果,其中P为弃风功率.可知,弃风现象主要发生在夜间,这是因为夜晚热负荷需求量增长,热电联产机组“以热定电”的运行方式使得夜间热电联产机组电出力随之增长,而夜间电负荷需求下降,限制了风电上网空间,造成大量的弃风.
图11
准确评估建筑热负荷的潜在弹性可以提升风电消纳效率.图12为是否考虑建筑热负荷需求弹性及不确定性的调度过程中风电消纳量(P1)的比较.其中情况1为考虑热负荷需求及不确定性;情况2为考虑热负荷需求但不考虑不确定性;情况3为不考虑热负荷需求弹性.可知,当考虑建筑热负荷需求弹性及不确定性,并将建筑热负荷灵活需求响应模型嵌入电-热综合能源系统优化调度中时,在大部分时间段内系统风电消纳量明显增加.相较于情况2和情况3,情况1在1 d内的弃风率分别降低了2.6%和5.1%.由此得出结论,考虑建筑热负荷弹性可以显著提高风电消纳,进一步证明了本文的热负荷需求弹性不确定性评估模型对提高系统和风电消纳的灵活性是必要且有效的.
图12
5 结语
本文研究了建筑热负荷需求弹性及认知不确定性的处理方法及应用,以促进电热综合能源系统协调优化调度中可变风电的消纳.将基于数据驱动的建筑热负荷需求弹性评估方法过程中产生的多源误差建模为认知不确定性,并通过改进的D-S证据理论对多源误差进行融合;采用拉丁超立方抽样方法和模糊聚类法生成表征热负荷需求弹性及认知不确定性的代表性场景;将构造的场景集嵌入电-热综合能源系统协调调度中,实现对建筑热负荷需求弹性及相关认知不确定性的综合考虑.在后续研究中,将针对建筑热负荷弹性不确定性在能源市场中的应用提出解决方案,即对充分考虑能价激励、用户行为下,建筑热负荷弹性不确定性评估及应用展开进一步研究.
参考文献
综合能源系统多时间尺度动态时域仿真关键技术
[J].
Key technologies for multi-time-scale dynamic time-domain simulation of integrated energy system
[J].
含多能微网群的区域电热综合能源系统分层自治优化调度
[J].
Hierarchical autonomous optimal dispatching of district integrated heating and power system with multi-energy microgrids
[J].
Economic dispatch of integrated energy systems with robust thermal comfort management
[J].DOI:10.1109/TSTE.5165391 URL [本文引用: 1]
Dispatch model of combined heat and power plant considering heat transfer process
[J].DOI:10.1109/TSTE.2017.2671744 URL [本文引用: 1]
计及储热备用效益的电热综合能源系统优化调度模型
[J].
Optimal dispatch model of electricity-heat integrated energy system considering reserved benefits of heat storage
[J].
考虑风电消纳的区域综合能源系统源荷协调经济调度
[J].
Source-load coordination economic dispatch method for regional integrated energy system considering wind power accommodation
[J].
计及P2H的电-热互联综合能源系统概率能量流分析
[J].
Analysis of probabilistic energy flow for integrated electricity-heat energy system with P2H
[J].
Dispatch model for CHP with pipeline and building thermal energy storage considering heat transfer process
[J].DOI:10.1109/TSTE.2018.2829536 URL [本文引用: 1]
Flexible dispatch of a building energy system using building thermal storage and battery energy storage
[J].DOI:10.1016/j.apenergy.2019.03.187 URL [本文引用: 2]
Experimental investigation of demand response potential of buildings: Combined passive thermal mass and active storage
[J].DOI:10.1016/j.apenergy.2020.115956 URL [本文引用: 1]
Estimating demand response potential under coupled thermal inertia of building and air-conditioning system
[J].DOI:10.1016/j.enbuild.2018.10.022 URL [本文引用: 1]
Multi-step ahead forecasting of heat load in district heating systems using machine learning algorithms
[J].DOI:10.1016/j.energy.2019.116085 URL [本文引用: 2]
数据驱动的综合能源系统负荷预测综述
[J].
Review of data-driven load forecasting for integrated energy system
[J].
Smart meter data-driven evaluation of operational demand response potential of residential air conditioning loads
[J].DOI:10.1016/j.apenergy.2020.115708 URL [本文引用: 1]
A simplified prediction model for energy use of air conditioner in residential buildings based on monitoring data from the cloud platform
[J].DOI:10.1016/j.scs.2020.102194 URL [本文引用: 1]
Optimal operation strategy for combined heat and power system based on solid electric thermal storage boiler and thermal inertia
[J].DOI:10.1109/ACCESS.2019.2958877 URL [本文引用: 1]
Demand response potential evaluation for residential air conditioning loads
[J].DOI:10.1049/gtd2.v12.19 URL [本文引用: 1]
Transmission-constrained unit commitment considering combined electricity and district heating networks
[J].DOI:10.1109/TSTE.2015.2500571 URL [本文引用: 1]
计及热力不确定性的智能建筑电-热联合鲁棒经济调度
[J].
Robust dispatch for electrical-thermal combined intelligent building considering impacts of uncertainties on thermal side
[J].
Graduation formula: A new method to construct belief reliability distribution under epistemic uncertainty
[J].DOI:10.23919/JSEE.5971804 URL [本文引用: 1]
D-S证据理论综述
[J].
A summary of D-S evidence theory
[J].
多源信息融合中一种新的证据合成算法
[J].
A novel combination rule of evidence theory in multi-source information fusion
[J].
基于大焦元的子焦元的信任函数逼近方法
[J].
The approximation method of the belief function based on the sub focal elements of the large focal elements
[J].
基于融合信息熵性质的信任函数概率逼近
[J].
A probability approximations of belief function based on fusion of the properties of information entropy
[J].
Clustering representative days for power systems generation expansion planning: Capturing the effects of variable renewables and energy storage
[J].DOI:10.1016/j.apenergy.2019.113603 URL [本文引用: 1]
Probabilistic evaluation of available load supply capability for distribution system
[J].DOI:10.1109/TPWRS.2013.2245924 URL [本文引用: 1]
Comparison of simulation methods for power system reliability indexes and their distributions
[J].DOI:10.1109/TPWRS.2008.919425 URL [本文引用: 1]
模糊聚类理论发展及应用的研究进展
[J].
Research progress on the development and application of fuzzy clustering theory
[J].
Risk-oriented renewable energy scenario clustering for power system reliability assessment and tracing
[J].DOI:10.1109/Access.6287639 URL [本文引用: 1]
Market-based integrated generation expansion planning of electric power system and district heating systems
[J].DOI:10.1109/TSTE.5165391 URL [本文引用: 1]
Combined analysis of electricity and heat networks
[J].DOI:10.1016/j.apenergy.2015.01.102 URL [本文引用: 1]
/
〈 |
|
〉 |
