CN110535128B - 基于用能舒适度的多区域综合能源系统协同调度方法 - Google Patents

基于用能舒适度的多区域综合能源系统协同调度方法 Download PDF

Info

Publication number
CN110535128B
CN110535128B CN201910839302.0A CN201910839302A CN110535128B CN 110535128 B CN110535128 B CN 110535128B CN 201910839302 A CN201910839302 A CN 201910839302A CN 110535128 B CN110535128 B CN 110535128B
Authority
CN
China
Prior art keywords
heat
temperature
user
load
power
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910839302.0A
Other languages
English (en)
Other versions
CN110535128A (zh
Inventor
刘鑫蕊
姜潇
孙秋野
张化光
黄博南
赵芯莹
孙丝萝
孟媛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northeastern University China
Original Assignee
Northeastern University China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Northeastern University China filed Critical Northeastern University China
Priority to CN201910839302.0A priority Critical patent/CN110535128B/zh
Publication of CN110535128A publication Critical patent/CN110535128A/zh
Application granted granted Critical
Publication of CN110535128B publication Critical patent/CN110535128B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/067Enterprise or organisation modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06Electricity, gas or water supply
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/12Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
    • H02J3/14Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by switching loads on to, or off from, network, e.g. progressively balanced loading
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02BCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
    • Y02B70/00Technologies for an efficient end-user side electric power management and consumption
    • Y02B70/30Systems integrating technologies related to power network operation and communication or information technologies for improving the carbon footprint of the management of residential or tertiary loads, i.e. smart grids as climate change mitigation technology in the buildings sector, including also the last stages of power distribution and the control, monitoring or operating management systems at local level
    • Y02B70/3225Demand response systems, e.g. load shedding, peak shaving
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S20/00Management or operation of end-user stationary applications or the last stages of power distribution; Controlling, monitoring or operating thereof
    • Y04S20/20End-user application control systems
    • Y04S20/222Demand response systems, e.g. load shedding, peak shaving

Abstract

本发明提供一种基于用能舒适度的多区域综合能源系统协同调度方法,包括:将电负荷分为可中断负荷、可调节负荷、可转移负荷三种类型;根据区域用能特性差异,将区域划分为低热电比的商业区和高热电比的居民区与工业区;根据采暖负荷建立基于气象预报的室内温度模型,结合热舒适度确定热需求斑马线;在传统燃料成本最优基础上附加环境成本与用能满意度成本建立优化目标;引入区域协同电力、热力供应机制实现区域间的整体热供需平衡。本发明基于用能舒适度提高了机组运行的灵活性,充分利用区域间负荷特性与用能舒适度进行联合调度达到能量管理的目的。

Description

基于用能舒适度的多区域综合能源系统协同调度方法
技术领域
本发明涉及综合能源系统运行优化领域,具体而言,涉及一种基于用能舒适度的多区域综合能源系统协同调度方法。
背景技术
随着能源枯竭以及环境污染的情况加剧,人们明显意识到传统的能源利用形式不符合现今社会的能源清洁高效发展的科学理念。综合能源系统可利用清洁能源实现能量的梯级利用,并凭借能效高、环境效益好等优势得到了领域内的广泛关注,是未来能源系统发展的重要方向。
目前综合能源系统优化通常以单区域为研究对象,但单区域的负荷特性比较单一、局限性强,从而制约了综合能源系统的优化结果。此外随着生活水平的提高,用户的用能舒适度也是能源调度应该关注的重点,而现有能源调度解决方案通常只针对电或者只针对热进行舒适度研究,缺乏基于用能舒适度的多区域综合能源系统协同调度方法。
发明内容
根据上述提出单区域研究负荷特性单一、局限性强,且未考虑用户舒适度的技术问题,而提供一种基于用能舒适度的多区域综合能源系统协同调度方法。本发明主要结合用户用能舒适度可实现某区域富裕热量或电量消纳,实现不同区域设备的协同调度,实现最优工作配合,提高综合能源系统经济性。
本发明采用的技术手段如下:
一种基于用能舒适度的多区域综合能源系统协同调度方法,包括:
S1、将电负荷分为可中断负荷、可调节负荷、可转移负荷三种类型,分别以断电时间、功率、转移时刻为衡量指标判断各类型电负荷用户的用电舒适度,将热负荷分为采暖用热、生产工艺用热两种类型,对于采暖用热以温度为衡量指标判断用户的用热舒适度;
S2、根据区域用能特性差异,将区域划分为低热电比的商业区和高热电比的居民区与工业区,并考虑热网动态特性建立传输延时与温度损耗模型;
S3、根据采暖负荷建立基于气象预报的室内温度模型,结合热舒适度确定热需求斑马线;
S4、在传统燃料成本最优基础上附加环境成本与用能满意度成本建立优化目标,所述用能满意度成本与用能舒适度相关;
S5、引入区域协同电力、热力供应机制,使电量富裕的区域将富余电量供给正处于用电高峰的区域,并采用热环网的传输形式进行热能交互,从而实现区域间的整体热供需平衡。
较现有技术相比,本发明具有以下优点:
本发明通过热环网将子区域连接形成多区域综合能源系统,将余热注入热网,也可以从热网获取热量。而多余的电量也可通过区域内电网进行协同消纳,基于用能舒适度也可提高机组运行的灵活性,由此可充分利用区域间负荷特性与用能舒适度进行联合调度达到能量管理的目的。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图做以简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例中多区域综合能源系统结构图。
图2是本发明实施例采用的CHP系统结构图。
图3是本发明实施例中热网传输延时的管道流量示意图。
图4是本发明实施例的区域能量协同机制流程图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
如图1-4所示,本发明提供了一种基于用能舒适度的多区域综合能源系统协同调度方法,包括:
S1、将电负荷分为可中断负荷、可调节负荷、可转移负荷三种类型,分别以断电时间、功率、转移时刻为衡量指标判断各类型电负荷用户的用电舒适度,将热负荷分为采暖用热、生产工艺用热两种类型,对于采暖用热以温度为衡量指标判断用户的用热舒适度。
具体地,可中断负荷模型为:
Figure GDA0003927632240000031
式中:μes,z为可中断负荷每次断电时产生的用户不舒适度;T为断电时长;Tes,best、Tes,max分别为最佳断电时长与最长断电时间。
可调节负荷模型为:
Figure GDA0003927632240000032
式中:μea,j为可调节负荷每次调节产生的用户不舒适度;Pa为调节后运行功率;α为调节裕度参数;Pea,N为可调负荷额定功率。
可转移负荷模型为:
Figure GDA0003927632240000041
式中:μet,k为可转移负荷每次转移产生的用户不舒适度;t为工作时刻;Tet,best、Tet,on、Tet,off分别为最佳启动时间、负荷可转移启、止时刻。
热舒适度模型为:
以温度为衡量指标判断用户对房间温度的热舒适度,采用模糊隶属度函数给出热舒适度模型。
Figure GDA0003927632240000042
式中:μh,i为用户i的热舒适度;Tn,i为用户i的室内温度;a由科学依据定为最适宜人们日常生活的25℃,上下偏差2℃为用户热舒适温度范围。
S2、根据区域用能特性差异,将区域划分为低热电比的商业区和高热电比的居民区与工业区,并考虑热网动态特性建立传输延时与温度损耗模型。
具体地,传输延时模型:
将热水管道出水温度定义为在Δt为管时间内对应的各部分流入温度的线性加权。
Figure GDA0003927632240000043
式中:mi,t为t时刻管道i中的水流量;mi,tΔt为从t时刻到t+Δt时刻内注入管道i的热水质量;αi,t为t-1时刻流出管道i的热水与其注入管道i时刻的时间差值;βi,t为t时刻流出管道i的热水与其注入管道i时刻的时间差值;t-αi,t为t-1时刻流出管道i的那部分热水注入管道i的时刻;t-βi,t为t时刻流出管道i的那部分热水注入管道i的时刻;
Figure GDA0003927632240000044
为t时刻出水管温度;
Figure GDA0003927632240000045
为t-αi,t时刻入水管温度;
Figure GDA0003927632240000046
为t-βi,t时刻入水管温度;
Figure GDA0003927632240000047
为k时刻入水管温度;Ai为管道i的横截面积;Li为管道i的长度。
其中:
Figure GDA0003927632240000048
Figure GDA0003927632240000051
温度损耗模型:
考虑到热水在管道传输的过程中会与管道外壁进行热交换,所以采用苏霍夫温降公式进行管道出口温度修正,修正温度如下:
Figure GDA0003927632240000052
其中Ta,t为t时刻的环境温度;Ji,t为温度损失系数。
计算方法为:
Figure GDA0003927632240000053
其中λi为管道的导热系数。
S3、根据采暖负荷建立基于气象预报的室内温度模型,结合热舒适度确定热需求斑马线。
具体地,用户室内温度模型:
Figure GDA0003927632240000054
式中为Qra,t=KSFSNS(Tav-Tn,t),
Figure GDA0003927632240000055
且Qra,t需满足Qra,t=cwmj(Tin-Tout);进而得到散热器功率与供水温度和气象温度的关系:
Figure GDA0003927632240000056
Figure GDA0003927632240000057
式中:Tn,t、Ta,t、Ts、Tav分别为用户在t时段的室内、气象温度和供水温度、散热器内热水的平均温度;Qra,t为t时段散热器功率;δ为空气热传导系数;ΔT为调度时间段时长;Tc为调度周期,单位为h;Ura,t表示散热器在t时段的开关状态,Ura,t=1表示t时段散热器开启;Ura,t=0表示t时段散热器关停;Ks为散热器的传热系数;Fs为散热器的散热面积;Ns为散热器组装片数、连接形式和安装形式的修正系数;h为传热系数;l为管道长度;cw为水的比热容;mj为与用户j相连的上级热水质量流。
热需求斑马线:
用户热舒适度μh,i与室内温度相关,在满足规定舒适度下限μ-h基础上,以0.2为μh的变化梯度,制定用户舒适度等级,每个等级对应的热需求形成斑马线,由此称之为热需求斑马线,用户可选择其任意斑马线作为热舒适度阈值。
Lh,i=cwmj(Ti μ-Tb)+Lh,loss
Figure GDA0003927632240000061
式中:Lh,i、Lh,loss、ΔLh,i分别为用户i的热需求、热损失、热需求变化区间,称之为基于热舒适度的热需求斑马带;Ti μ、Tb
Figure GDA0003927632240000062
分别为用户设定舒适度所对应的温度、热水管道水源基准温度、最舒适温度。
S4、在传统燃料成本基础上附加环境成本与用能满意度成本建立优化目标,所述用能满意度成本与用能舒适度相关。
具体地,以每日总费用、环境成本、用户不适度共同组成优化的目标函数。其数学表达式如下:
Figure GDA0003927632240000063
Figure GDA0003927632240000064
Figure GDA0003927632240000065
Figure GDA0003927632240000066
Figure GDA0003927632240000067
式中:f为供能成本函数;
Figure GDA0003927632240000068
为购电花费;
Figure GDA0003927632240000069
为购气花费;
Figure GDA00039276322400000610
为环境费用;
Figure GDA00039276322400000611
为用户满意度费用;pe为上级电网电价;
Figure GDA00039276322400000612
为从上级电网的购电功率;pg为天然气单价;Vmt为燃气轮机耗气量;Vgb为燃气锅炉耗气量;ζgrid、ζmt、ζgb分别为电力、燃气轮机、燃气锅炉的折算后的CO2排放因子。
Figure GDA00039276322400000613
分别为燃气轮机和燃气锅炉的输出功率;
Figure GDA00039276322400000614
分别为用户电、热不满意度折算成本。
约束条件包括:
Figure GDA00039276322400000615
Figure GDA00039276322400000616
Figure GDA00039276322400000617
Figure GDA0003927632240000071
Figure GDA0003927632240000072
Figure GDA0003927632240000073
Figure GDA0003927632240000074
Figure GDA0003927632240000075
Figure GDA0003927632240000076
Figure GDA0003927632240000077
X=(μes,zea,jet,kh,i)
式中:
Figure GDA0003927632240000078
为燃气轮机每小时的天然气耗量;Hg为天然气热值;
Figure GDA0003927632240000079
为燃气轮机的发电效率;ηloss为能量损耗率;
Figure GDA00039276322400000710
为热回收效率;
Figure GDA00039276322400000711
为燃气锅炉产热功率;
Figure GDA00039276322400000712
为燃气锅炉单位小时消耗的天然气体积;
Figure GDA00039276322400000713
为燃气锅炉的产热效率;
Figure GDA00039276322400000714
分别为燃气轮机发电的上下限;
Figure GDA00039276322400000715
Figure GDA00039276322400000716
分别为燃气轮机产热的上下限;
Figure GDA00039276322400000717
分别为购电的上下限:
Figure GDA00039276322400000719
Figure GDA00039276322400000720
分别为燃气锅炉产热的上下限;X为舒适度向量。
S5、引入区域协同电力、热力供应机制,使电量富裕的区域将富余电量供给正处于用电高峰的区域,并采用热环网的传输形式进行热能交互,从而实现区域间的整体热供需平衡。
本实施例中,综合能源系统以热电联供系统为例,多区域综合能源系统结构如图1所示,每个子区域建有CHP系统,CHP系统结构图如图2所示,其中电热主要来源于燃气轮机,燃气锅炉作为补燃装置投入,并通过热环网将子区域连接形成多区域综合能源系统。CHP可将余热注入热网,也可以从热网获取热量,CHP多余的电量也可通过区域内电网进行协同消纳,基于用能舒适度也可提高机组运行的灵活性,由此可充分利用区域间负荷特性与用能舒适度进行联合调度达到能量管理的目的。图3所示为热网传输延时的管道流量示意图。
具体地,由于区域可进行协调配合所以可对产生能量剩余的区域采取区域协同机制促进能量消纳,增加经济性。图4为能量协同分配流程图,其协同分配收益度如下:
(1)电协同收益度
因为区域不可能富余大量电完全供给另一区域,所以有
Figure GDA00039276322400000721
但由于需满足区域热舒适度,所以若电余量过多会出现不能完全消纳的情况,区域消纳电的上限为
Figure GDA0003927632240000081
且规划时不允许所有区域在某一时刻在以热定电的同时全部有电剩余。
Figure GDA0003927632240000082
式中:P′e为区域剩余电量之和;
Figure GDA0003927632240000083
为区域购电量。
Figure GDA0003927632240000084
为改变机组出力而能增加的耗电上限。
首先定义函数
Figure GDA0003927632240000085
然后进行按如下计方案进行热分配计算,计算顺序有先后,但最终按计算总结果进行统一电余量协同分配:
1)
Figure GDA0003927632240000086
1.当
Figure GDA0003927632240000087
时,富裕电量可完全消纳
2.当
Figure GDA0003927632240000088
时,
Figure GDA0003927632240000089
式中:
Figure GDA00039276322400000810
为区域消纳电量。
2)
Figure GDA00039276322400000811
1.当
Figure GDA00039276322400000812
时,电协同收益度
Figure GDA00039276322400000813
如下:
Figure GDA00039276322400000814
式中:
Figure GDA00039276322400000815
为区域Ωi在t时刻的电价。
先满足电协同收益度高区域的购电量。
2.当
Figure GDA00039276322400000816
时,补充区域购电后仍剩余的电
Figure GDA00039276322400000817
然后按基于调动机组的电收益较高区域优先获得电量的顺序,进行区域最大电量消耗,若有剩余继续供给下一区域。
基于调动机组的收益判据如下:
Figure GDA0003927632240000091
Figure GDA0003927632240000092
Figure GDA0003927632240000093
Figure GDA0003927632240000094
Figure GDA0003927632240000095
Figure GDA0003927632240000096
Figure GDA0003927632240000097
Figure GDA0003927632240000098
式中:Ωi为划分的区域;
Figure GDA0003927632240000099
为电协同收益度;
Figure GDA00039276322400000910
为协同前t时刻区域运行成本;
Figure GDA00039276322400000911
为协同后t时刻区域运行成本;
Figure GDA00039276322400000912
Figure GDA00039276322400000913
分别为电热成本参数;
Figure GDA00039276322400000914
为协同后区域发电功率;
Figure GDA00039276322400000915
为协同后区域发热功率;
Figure GDA00039276322400000916
为协同后区域燃气轮机功率;
Figure GDA00039276322400000917
为热电比系数;
Figure GDA00039276322400000918
为燃气轮机出力上限。
3)
Figure GDA00039276322400000919
首先按如下电协同收益度进行,购电的补余:
Figure GDA00039276322400000920
然后若有电剩余,进行按基于调动机组的电收益较高区域优先获得电量的顺序,进行区域最大电量消耗,若有剩余继续供给下一区域。
基于调动机组的收益判据上文已提及。
4)
Figure GDA00039276322400000921
按基于调动机组的电收益较高区域优先获得电量的顺序,进行区域最大电量消耗,若有剩余继续供给下一区域。
基于调动机组的收益判据上文已提及。
(2)热协同收益度
因为各区域均采用一热定电的模式,且供热设备都可满足本区域热需求,在某区域用电高峰需要向电网购电时,协同优化为了减少购电会增加机组出力,此时会产生协同热剩余,前提不存在三个区域同时热剩余。
然后进行按如下计方案进行热分配计算,计算顺序有先后,但最终按计算总结果进行统一热余量协同分配:
1)只有一个区域有热剩余
首先在自区域以提升区域热舒适度的方式将热储存起来;当该区域热舒适度达到上限时仍有热余量,对其余区域进行如下判断:
1.若只有一个区域存在燃气轮机补燃,则优先对该区域进行燃气轮机补燃,若有热余量然后寻找热损失最小路径通过增加区域热舒适度进行热消纳,当该区域舒适度达到上限仍有热余量,增加下一区域热舒适度进行热消纳。
2.若有两个区域存在燃气轮机补燃,先寻找热损失最小路径进行热分配,若有热余量继续对另一区域进行热分配减小燃气轮机出力,若仍有热余量,寻找热损失最小路径通过增加区域热舒适度进行热消纳,当该区域舒适度达到上限仍有热余量,增加下一区域热舒适度进行热消纳。
3.若无区域存在燃气轮机补燃,寻找热损失最小路径通过增加区域热舒适度进行热消纳,当该区域舒适度达到上限仍有热余量,增加下一区域热舒适度进行热消纳。
2)有两个区域热剩余
首先在自区域以提升区域热舒适度的方式将热储存起来;当该区域热舒适度达到上限时仍有热余量,对剩余区域进行如下操作:
寻找热损失最小路径通过增加区域热舒适度进行热消纳,若存在热余量,则停止,进行弃热。
(3)热损最小路径确定方法
1)首先对三个区域进行编号分别为:Ω1、Ω2、Ω3
2)定义函数
Figure GDA0003927632240000111
3)定义集合:A和B。
Figure GDA0003927632240000112
时,i∈A;
Figure GDA0003927632240000113
时,j=i∈B;
4)定义管道热传输优先度:
Figure GDA0003927632240000114
式中:lij为有余热区域管道到无余热区域管道的长度。
5)进行管道热传输优先级排序,优先度大的管道首先进行热传输计算。
采用CPLEX进行模型求解,给出最优能量协同调度策略。
上述实施例首先建立了包含热传输延时、传输损耗、用能舒适度等约束的综合能源系统模型,其中将舒适度范围内房屋蓄热温度作为热负荷的指标,并提出在舒适度范围内热负荷的可调节斑马带以提升机组的灵活性,在日内调度中可有效降低系统运行成本。其次,建立多区域能量协同配合机制对区域剩余能量进行协调配合,提高能源利用率,当子区域处于用电高峰时也可采用燃气轮机满发状态减小向上级电网购电量,若产生余热可进行多区域的能量协同配合,进而提升系统运行效率。
下面通过具体的应用实例对本发明的效果进行进一步说明。
以某综合能源区域为例,如图1所示,该区域划分为商业区、居民区和工业区三个子区域,其区域热管道长度分别为0.8km、1.2km、0.6km,区域电价分别为0.9元/kWh、0.5元/kWh、0.7元/kWh,天然气价格为2.3元/m3。为便于比较分析:设置案例1为基本案例,该基本案例既不考虑用能舒适度特性,也不考虑区域间的能量协同配合;案例2仅考虑用能舒适度特性;案例3为本发明实施例所提模型。各案例运行结果如表1所示:
表1运行结果
Figure GDA0003927632240000121
由表1可知,案例2相对于案例1节省了0.545万元,节省率为2.84%。其中,案例2由于考虑了用能舒适度,即电负荷可在舒适度内调节、转移或中断用电器以减少用电量,从而减少峰时购电;而基于热舒适度的热负荷在热需求斑马带内是可调的,由此在用电高峰期燃气轮机多发电,从而减少购电费用。案例3相对于案例2节省了0.849万元,节省率为4.43%;相对于案例1节省了1.394万元,节省率为7.27%。因为区域能量协同调度的加入可将富裕能量分配给其它区域提高能量利用率,又因为各区域的协调配合可以使得各区域机组运行灵活性更大,进而提升经济效益。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (4)

1.一种基于用能舒适度的多区域综合能源系统协同调度方法,其特征在于,包括:
S1、将电负荷分为可中断负荷、可调节负荷、可转移负荷三种类型,分别以断电时间、功率、转移时刻为衡量指标判断各类型电负荷用户的用电舒适度,将热负荷分为采暖用热、生产工艺用热两种类型,对于采暖用热以温度为衡量指标判断用户的用热舒适度;用户的用热舒适度模型为:
Figure FDA0003927632230000011
其中,μh,i为用户i的热舒适度;Tn,i为用户i的室内温度;a为最适宜人们日常生活温度,其上下偏差2℃为用户热舒适温度范围;
S2、根据区域用能特性差异,将区域划分为低热电比的商业区和高热电比的居民区与工业区,并考虑热网动态特性建立传输延时与温度损耗模型;
S3、根据采暖负荷建立基于气象预报的室内温度模型,结合热舒适度确定热需求斑马线,包括:
S31、建立用户室内温度模型,进而得到散热器功率与供水温度和气象温度的关系,所述用户室内温度模型为:
Figure FDA0003927632230000012
Qra,t=KSFSNS(Tav-Tn),
Figure FDA0003927632230000013
Qra,t=cwmj(Tin-Tout)
式中:Tn,t为t时段的室内温度,Ta,t为t时段的气象温度,Tav为t时段的散热器内热水的平均温度,Qra,t为t时段散热器功率,δ为空气热传导系数,ΔT为调度时间段时长;Tc为调度周期,Ura,t表示散热器在t时段的开关状态,Ura,t=1表示t时段散热器开启,Ura,t=0表示t时段散热器关停,Ks为散热器的传热系数,Fs为散热器的散热面积,Ns为散热器组装片数、连接形式和安装形式的修正系数,cw为水的比热容,mj为与用户j相连的上级热水质量流,
S32、在满足规定舒适度下限μ h基础上,以固定步长为μh的变化梯度,制定用户舒适度等级,每个等级对应的热需求形成斑马线:
Lh,i=cwmj(Ti μ-Tb)+Lh,loss
Figure FDA0003927632230000024
式中:Lh,i为用户i的热需求,Lh,loss为用户i的热损失,ΔLh,i为热需求变化区间,基于热舒适度的热需求斑马带;Ti μ为用户设定舒适度所对应的温度,Tb为水管道水源基准温度、
Figure FDA0003927632230000021
为最舒适温度;
S4、在传统燃料成本最优基础上附加环境成本与用能满意度成本建立优化目标,所述用能满意度成本与用能舒适度相关;
S5、引入区域协同电力、热力供应机制,使电量富裕的区域将富余电量供给正处于用电高峰的区域,并采用热环网的传输形式进行热能交互,从而实现区域间的整体热供需平衡。
2.根据权利要求1所述的多区域综合能源系统协同调度方法,其特征在于,可中断负荷用户用电舒适度模型为:
Figure FDA0003927632230000022
其中,μes,z为可中断负荷每次断电时产生的用户不舒适度,T为断电时长,Tes,best为最佳断电时长,Tes,max为最长断电时间;
可调节负荷用户用电舒适度模型为:
Figure FDA0003927632230000023
其中,μea,j为可调节负荷每次调节产生的用户不舒适度,Pa为调节后运行功率,α为调节裕度参数,Pea,N为可调负荷额定功率;
可转移负荷用户用电舒适度模型为:
Figure FDA0003927632230000031
其中:μet,k为可转移负荷每次转移产生的用户不舒适度,t为工作时刻,Tet,best为最佳启动时间,Tet,on为负荷可转移启动时刻,Tet,off为负荷可转移终止时刻。
3.根据权利要求1或2所述的多区域综合能源系统协同调度方法,其特征在于,传输延时与温度损耗模型为:
传输延时:
Figure FDA0003927632230000032
Figure FDA0003927632230000033
Figure FDA0003927632230000034
其中,mi,t为t时刻管道i中的水流量;mi,tΔt为从t时刻到t+Δt时刻内注入管道i的热水质量;αi,t为t-1时刻流出管道i的热水与其注入管道i时刻的时间差值;βi,t为t时刻流出管道i的热水与其注入管道i时刻的时间差值;
Figure FDA0003927632230000035
为t时刻出水管温度;
Figure FDA0003927632230000036
为t-αi,t时刻入水管温度;
Figure FDA0003927632230000037
为t-βi,t时刻入水管温度;
Figure FDA0003927632230000038
为k时刻入水管温度;Ai为管道i的横截面积;Li为管道i的长度;
温度损耗:
Figure FDA0003927632230000039
其中,Ta,t为t时刻的环境温度;Ji,t为温度损失系数,且
Figure FDA00039276322300000310
其中,λi为管道的导热系数。
4.根据权利要求2所述的多区域综合能源系统协同调度方法,其特征在于,步骤S4中优化目标为:
Figure FDA0003927632230000041
Figure FDA0003927632230000042
Figure FDA0003927632230000043
Figure FDA0003927632230000044
Figure FDA0003927632230000045
其中:f为供能成本函数,
Figure FDA0003927632230000046
为购电花费,
Figure FDA0003927632230000047
为购气花费,
Figure FDA0003927632230000048
为环境费用,
Figure FDA0003927632230000049
为用户满意度费用,pe,t为上级电网电价,
Figure FDA00039276322300000410
为从上级电网的购电功率,pg为天然气单价,Vmt为燃气轮机耗气量,Vgb为燃气锅炉耗气量;ζgrid为电力折算后的CO2排放因子,ζmt为燃气轮机折算后的CO2排放因子,ζgb为燃气锅炉的折算后的CO2排放因子,
Figure FDA00039276322300000411
为燃气轮机的输出功率,
Figure FDA00039276322300000412
为燃气锅炉的输出功率,
Figure FDA00039276322300000413
为用户电不满意度折算成本,
Figure FDA00039276322300000414
为用户热不满意度折算成本。
CN201910839302.0A 2019-09-05 2019-09-05 基于用能舒适度的多区域综合能源系统协同调度方法 Active CN110535128B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910839302.0A CN110535128B (zh) 2019-09-05 2019-09-05 基于用能舒适度的多区域综合能源系统协同调度方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910839302.0A CN110535128B (zh) 2019-09-05 2019-09-05 基于用能舒适度的多区域综合能源系统协同调度方法

Publications (2)

Publication Number Publication Date
CN110535128A CN110535128A (zh) 2019-12-03
CN110535128B true CN110535128B (zh) 2022-12-13

Family

ID=68667188

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910839302.0A Active CN110535128B (zh) 2019-09-05 2019-09-05 基于用能舒适度的多区域综合能源系统协同调度方法

Country Status (1)

Country Link
CN (1) CN110535128B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112418477B (zh) * 2020-06-23 2023-07-07 深圳职业技术学院 一种基于区域能源中心的居民用能双目标优化方法
CN111967728B (zh) * 2020-07-27 2022-05-10 国网河北省电力有限公司营销服务中心 计及用能舒适度时变性的商场楼宇调峰容量评估方法
CN111952980B (zh) * 2020-08-11 2022-12-09 合肥工业大学 考虑需求侧响应不确定性的综合能源系统优化方法及系统
CN112594761B (zh) * 2020-11-11 2021-12-24 华电电力科学研究院有限公司 一种集中式区域供能方法
CN112836849A (zh) * 2020-12-21 2021-05-25 北京华能新锐控制技术有限公司 一种考虑风电不确定性的虚拟电厂调度方法
CN113489063B (zh) * 2021-07-01 2022-07-15 上海发电设备成套设计研究院有限责任公司 一种区域综合能源的联合调度系统及方法
CN113610316B (zh) * 2021-08-19 2023-10-03 东北电力大学 不确定环境下考虑综合需求响应的园区综合能源系统优化调度方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107665384A (zh) * 2017-10-27 2018-02-06 天津大学 一种含多区域能源站的电力‑热力综合能源系统调度方法
CN108629470A (zh) * 2017-03-17 2018-10-09 华北电力大学 基于非合作博弈的多能互补系统能量管理与优化运行
CN109193628A (zh) * 2018-09-18 2019-01-11 华北电力大学 一种基于一致性的综合能源系统能量管理方法
CN110046773A (zh) * 2019-05-30 2019-07-23 东北大学 综合能源系统多主体合作优化运行与成本效益分配方法
CN110188492A (zh) * 2019-06-04 2019-08-30 南通大学 一种考虑热网特性的冷热电联供微网优化调度方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9651929B2 (en) * 2014-09-29 2017-05-16 International Business Machines Corporation HVAC system control integrated with demand response, on-site energy storage system and on-site energy generation system
CN106998079B (zh) * 2017-04-28 2020-05-05 东南大学 一种热电联合优化调度模型的建模方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108629470A (zh) * 2017-03-17 2018-10-09 华北电力大学 基于非合作博弈的多能互补系统能量管理与优化运行
CN107665384A (zh) * 2017-10-27 2018-02-06 天津大学 一种含多区域能源站的电力‑热力综合能源系统调度方法
CN109193628A (zh) * 2018-09-18 2019-01-11 华北电力大学 一种基于一致性的综合能源系统能量管理方法
CN110046773A (zh) * 2019-05-30 2019-07-23 东北大学 综合能源系统多主体合作优化运行与成本效益分配方法
CN110188492A (zh) * 2019-06-04 2019-08-30 南通大学 一种考虑热网特性的冷热电联供微网优化调度方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Distributed_Optimal_Energy_Management_for_Energy_Internet;Huaguang Zhang;《IEEE TRANSACTIONS ON INDUSTRIAL INFORMATICS》;20171231;第3081-3097页 *
Thermal Profiling of Residential Energy Use;Adrian Albert;《IEEE TRANSACTIONS ON POWER SYSTEMS》;20150331;第602-611页 *
参与双重市场的区域综合能源系统日前优化调度;黄伟;《电力系统自动化》;20190625;第68-75页 *
面向综合能源系统的静态电压稳定性分析;黄雨佳;《中国电机工程学报》;20190831;第44-53页 *

Also Published As

Publication number Publication date
CN110535128A (zh) 2019-12-03

Similar Documents

Publication Publication Date Title
CN110535128B (zh) 基于用能舒适度的多区域综合能源系统协同调度方法
CN108258679B (zh) 考虑供热网储热特性的电-热综合能源系统优化调度方法
WO2018130231A1 (zh) 一种基于热网和房屋热惯性的综合能源系统优化方法
CN109858759A (zh) 一种工业园区综合能源平衡调度方法
CN111738502A (zh) 促进富余风电消纳的多能互补系统需求响应运行优化方法
CN110209135B (zh) 基于微型热电联产多时间尺度的家庭能源优化调度方法
CN109447323A (zh) 一种计及节点热价的综合能源系统两阶段容量配置方法
CN108416472A (zh) 一种区域多源供冷供热系统优化调度方法
CN110968827B (zh) 一种多区域综合能源系统优化配置方法
CN112165122A (zh) 一种综合能源系统的运行方法及系统
CN111553544B (zh) 基于一致性算法的工业园区分布式综合需求响应方法
CN107392366A (zh) 计及用户行为分析的综合能源系统供需双侧综合优化方法
CN107749645A (zh) 一种控制高压大容量储热供暖装置的方法
CN115857348A (zh) 考虑两联供热泵舒适供能的分布式能量系统容量优化方法
CN114154328A (zh) 提升灵活性的电热综合能源系统两阶段随机优化调度方法
CN112488363B (zh) 基于广义储能的多能源电力系统优化调度方法
Pang et al. Collaborative power tracking method of diversified thermal loads for optimal demand response: A MILP-Based decomposition algorithm
CN116341881B (zh) 计及热网灵活性的电-热系统鲁棒超前调度方法及系统
Gan-yun et al. Optimal scheduling of regional integrated energy system considering integrated demand response
CN113128799A (zh) 能源管控方法、装置、电子设备及计算机存储介质
CN113091124B (zh) 一种基于区域供能能力的电热互补供热系统分区调控方法
CN113761727B (zh) 一种含分布式电热泵的热电联合系统优化调度模型构建方法
CN111241459B (zh) 一种考虑自趋优的工业园区分布式综合需求响应方法
CN113435708A (zh) 一种基于综合需求响应的楼宇综合能源管理方法与系统
Zhou et al. Day-ahead and intraday optimization models for demand response of the aggregators of thermostatically controlled loads

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant