CN116452009B - 一种考虑时间关联特性的综合能源枢纽灵活域评估方法 - Google Patents
一种考虑时间关联特性的综合能源枢纽灵活域评估方法 Download PDFInfo
- Publication number
- CN116452009B CN116452009B CN202310175488.0A CN202310175488A CN116452009B CN 116452009 B CN116452009 B CN 116452009B CN 202310175488 A CN202310175488 A CN 202310175488A CN 116452009 B CN116452009 B CN 116452009B
- Authority
- CN
- China
- Prior art keywords
- time
- energy
- flexible domain
- hub
- comprehensive energy
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 239000011159 matrix material Substances 0.000 claims abstract description 53
- 230000008878 coupling Effects 0.000 claims abstract description 27
- 238000010168 coupling process Methods 0.000 claims abstract description 27
- 238000005859 coupling reaction Methods 0.000 claims abstract description 27
- 238000011156 evaluation Methods 0.000 claims abstract description 6
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 108
- 239000003345 natural gas Substances 0.000 claims description 54
- 238000004146 energy storage Methods 0.000 claims description 41
- 239000007789 gas Substances 0.000 claims description 28
- 230000036962 time dependent Effects 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000005266 casting Methods 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 abstract description 6
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 abstract description 5
- 230000009286 beneficial effect Effects 0.000 abstract description 5
- 229910052799 carbon Inorganic materials 0.000 abstract description 5
- 238000004134 energy conservation Methods 0.000 abstract description 3
- 230000009467 reduction Effects 0.000 abstract description 3
- 230000005611 electricity Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000007246 mechanism Effects 0.000 description 4
- 230000002776 aggregation Effects 0.000 description 2
- 238000004220 aggregation Methods 0.000 description 2
- 239000000969 carrier Substances 0.000 description 2
- 239000000446 fuel Substances 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 238000005096 rolling process Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000011426 transformation method Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0637—Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/80—Management or planning
- Y02P90/82—Energy audits or management systems therefor
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Mathematical Analysis (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Game Theory and Decision Science (AREA)
- Primary Health Care (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Quality & Reliability (AREA)
- Development Economics (AREA)
- Mathematical Optimization (AREA)
- Operations Research (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种考虑时间关联特性的综合能源枢纽灵活域评估方法。方法包括:建立综合能源枢纽的多能流模型,将效率参数矩阵输入,输出连接关系参数矩阵;建立灵活域模型,将连接关系参数矩阵输入,输出时间无关灵活域和时间耦合灵活域;采用射线投影方法获取灵活域信息,实现综合能源枢纽灵活域的评估。本发明方法基于可获取考虑时间关联特征的综合能源枢纽灵活域,可提升对能源枢纽的灵活资源的利用水平,提升系统的安全性,充分发挥多能替代,有利于节能减碳。
Description
技术领域
本发明涉及了一种综合能源枢纽灵活域评估方法,具体涉及一种考虑时间关联特性的综合能源枢纽灵活域评估方法。
背景技术
不同能源载体之间的相互作用可有效提高能源效率,减少碳排放。作为多能源用户的典型形式,由燃气锅炉、储能设备等组成的综合能源枢纽已被广泛引入,以满足终端负荷的多能源需求(如工业制造和商业建筑)。综合能源枢纽可以调整不同时期的能源消耗,实现不同能源载体的能源替代。在不牺牲终端负荷用能需求的情况下调整多能源输入的能力被定义为综合能源枢纽的运行灵活性。通过多能源协调,综合能源枢纽可以为能源系统提供灵活调节资源,有利于能源系统的安全可靠运行。
综合能源枢纽的灵活域评估是能源系统利用其灵活资源的基础。在传统灵活域评估的过程中,主要研究了独立时间间隔的灵活域,决策变量只与当前系统状态有关,模拟了综合能源枢纽的静态运行特性。然而,除静态约束以外,综合能源枢纽的灵活域还会受到储能设备的充电/放电特征的影响,因此综合能源枢纽的灵活域具有时间关联特征。考虑到时间耦合和多能源相互作用,采用传统方法难以评估不同时间间隔的综合能源枢纽灵活域。因此,亟需提出新的综合能源枢纽灵活域评估方法,以有效考虑综合能源枢纽的时间关联特性,使得能源系统充分挖掘综合能源枢纽的灵活资源并安全经济低碳运行。
发明内容
为了解决背景技术中存在的问题,本发明所提供一种考虑时间关联特性的综合能源枢纽灵活域评估方法。本发明方法通过对考虑时间关联特性的综合能源枢纽灵活域进行表征,可指导综合能源枢纽调度中心向综合能源系统调度中心提供灵活资源,可提升综合能源系统调度中心对综合能源枢纽的灵活资源的利用水平,提升系统的安全性,充分发挥多能替代,有利于节能减碳。
本发明采用的技术方案是:
本发明综合能源枢纽灵活域评估方法包括如下步骤:
步骤1)建立考虑能源设备运行约束的综合能源枢纽的多能流模型,将效率参数矩阵输入多能流模型中,多能流模型输出连接关系参数矩阵。
步骤2)建立综合能源枢纽的灵活域模型,将连接关系参数矩阵输入灵活域模型中,灵活域模型输出综合能源枢纽的时间无关灵活域和时间耦合灵活域。
步骤3)根据时间无关灵活域和时间耦合灵活域,采用射线投影方法获取灵活域信息,即获取综合能源枢纽的电力输入和天然气输入的可调范围,实现综合能源枢纽灵活域的评估。
所述的步骤1)中,综合能源枢纽包括燃气锅炉、电储能设备、电热泵和热电联产机组,综合能源枢纽和综合能源系统以及终端能源负荷连接,综合能源系统提供电力和天然气能源,并通过综合能源枢纽为终端能源负荷提供热力、电力和天然气能源。
燃气锅炉是一种通过消耗天然气来产热的装置;电储能设备是指能将电力储存的装置;电热泵是消耗电能实现热能从低温到高温传递的装置;热电联产机组使用天然气作为燃料产生电力和热能的装置。基于综合能源系统提供的电力和天然气,满足终端能源负荷侧用户的电力、天然气、热等能源需求。
所述的步骤1)中,建立的考虑能源设备运行约束的综合能源枢纽的多能流模型,具体如下:
[GLEH,PLEH]T=Ax
Bx=[HTL,GTL,PTL]T
其中,PLEH和GLEH分别表示从综合能源系统到综合能源枢纽的电力和天然气输入;A表示综合能源枢纽的效率参数矩阵;x表示综合能源枢纽的决策变量;B表示综合能源枢纽的连接关系参数矩阵;HTL、PTL和GTL分别表示从综合能源枢纽到终端能源负荷的热力、电力和天然气输出。
所述的综合能源枢纽的决策变量x具体如下:
其中,HGB、HEHP和HCHP分别表示综合能源枢纽中燃气锅炉、电热泵和热电联产机组的热功率输出集合;PCHP和PES分别表示综合能源枢纽中热电联产机组和电储能设备的电功率输出集合;PN和GN分别表示终端能源负荷的电力负荷和天然气负荷;EES表示电储能设备的荷电状态集合。
所述的综合能源枢纽的效率参数矩阵A具体如下:
其中,ηGB和ηEHP分别表示燃气锅炉和电热泵的转换效率;和/>分别表示热电联产机组的热能和电能输出效率。
第一行中的元素表示决策变量与综合能源枢纽的电力输入之间的关系。第二行中的元素表示决策变量和天然气输入之间的关系。
所述的综合能源枢纽的连接关系参数矩阵B具体如下:
其中,连接关系参数矩阵B中的元素1表示综合能源枢纽中的不同装置之间存在连接关系,元素0表示综合能源枢纽中的不同装置之间无连接关系。
连接关系参数矩阵B中的不同行描述了决策变量和综合能源枢纽的多个能量输出之间的关系。
所述的能源设备运行约束具体如下:
EES,t=EES,t-1-PES,t-1Δt
0≤HGB≤HGB,max
0≤HEHP≤HEHP,max
其中,EES,t和EES,t-1分别表示电储能设备在t时刻和t-1时刻的荷电状态集合;PES,t-1表示电储能设备在t-1时刻的充放电率集合;Δt表示时间间隔;EES,min和EES,max分别表示电储能设备的最小和最大荷电状态集合;PEcS,max表示电储能设备的最大充放电率集合;HGB,max、HEHP,max和HCHP,max分别表示燃气锅炉、热电联产机组和热电联产机组的最大热功率输出集合;PCHP,min和PCHP,max分别表示热电联产机组的最小和最大电功率输出集合。
能源设备运行约束为考虑到时间特征的不同能源设备的运行约束,综合能源枢纽中的一些设备的运行状态不仅受到当前时间状态的约束,还会受到其他时间运行状态的约束,如电储能。综合能源枢纽中的没有时间耦合约束的能源设备包括如燃气锅炉、电热泵和热电联产机组。
多能流模型的紧凑形式可以表示为:
Vαxt≤Vβxt-1+U
其中,Vα和Vβ分别表示第一和第二常数矩阵;xt和xt-1分别表示t时刻和t-1时刻的决策变量;U表示参数向量;I8×8和08×8分别表示的8乘8的单位矩阵和0矩阵;Λα和Λβ分别表示第一和第二常数矩阵;xmax和xmin分别表示决策变量的最大值和最小值;EES,max和EES,min分别表示电储能设备的最大和最小荷电状态。
将多能流划分为两部分,即时间无关部分和时间耦合部分,具体如下:
其中,GLEH和PLEH分别表示综合能源枢纽的天然气输入和电力输入;A1是通过将时间耦合变量的列置零,A2是通过将时间无关变量的列置零,A1和A2均从A中通过矩阵变换方法计算出来。
所述的步骤2)中,建立的综合能源枢纽的灵活域模型具体如下:
其中,Ψ1表示综合能源枢纽的时间无关灵活域;表示闵可夫斯基求和;Ψ2表示综合能源枢纽的时间耦合灵活域;Vα,1表示时间无关灵活域Ψ1的常数矩阵;U1表示时间无关灵活域Ψ1的参数矩阵;GLEH,1和PLEH,1分别表示综合能源枢纽中时间无关灵活域Ψ1的天然气输入和电力输入;A1表示时间无关向量;GLEH,2和PLEH,2分别表示时间耦合灵活域Ψ2的天然气输入和电力输入;A2表示时间有关向量。
Ψ1是以天然气输入和电力输入为轴的二维区域,与综合能源枢纽的多能源替代密切相关。Ψ1中的能量流约束只保留了与时间无关的约束,由Vα,1和U1表示。Ψ2可以表征为与电力输入相对应的一维区间,与储能的充电/放电特征有关。关于时间耦合灵活域Ψ2,其时间耦合约束是指储能约束。因此,我们可以将Ψ2重新表述为一个储能系统模型。PES,t是由能源系统调度中心决定的调度指令的一部分,综合能源枢纽只需要提交容量限制(EES,min,EES,max)和最大充放电比例PEcS,max,这些都是每个时间段的常数。PES,t和EES,t分别代表t时刻的电力输入和储能设备水平的向量。
与时间无关的灵活域Ψ1是一个高维多面体,对能量输入向量的投影是一个二维多边形Θ,其可代表综合能源枢纽与时间无关的灵活性,具体如下:
综合能源枢纽的灵活性可以表示为两个多面体的聚合,即能源系统调度中心可以获得有关灵活域的信息Ψ,用于调度和风险管理等。灵活域模型具体可表示如下:
上述公式表示一个高维多面体,其中电能存在时间耦合约束,这可能使综合能源枢纽在不同时间间隔的运行灵活性相互影响。基于t-1的状态,综合能源枢纽调度中心来评估综合能源枢纽在时间t的灵活性,而综合能源枢纽的运行状态是由能源系统调度中心在能源系统层面调度的。综合能源枢纽在一天中不同时间的灵活性评估和运行调度是一个滚动过程。
综合能源枢纽调度中心是指负责综合能源枢纽运行的机构,可对综合能源枢纽内部设备的运行状态进行调整。综合能源枢纽运行的灵活性主要是指从电-气耦合系统输入到综合能源枢纽的能量的调节能力,它受制于综合能源枢纽中的设备(如热电联产机组)的运行限制和终端负荷的能源需求。由于能源输入包含两种形式,即电力和燃气,因此灵活域被描述为一个二维空间。综合能源枢纽的电力和天然气输入可以在这个区域内进行调节,这可以为电-气耦合系统提供灵活资源。灵活域是指综合能源枢纽的电力输入和天然气输入的可调范围,可表示为以电力输入为横坐标、天然气输入为纵坐标的二维区域。电-气耦合系统调度中心是指负责综合能源枢纽实际调度运行的机构,可以调整综合能源枢纽中各个装置的运行状态。
所述的时间无关灵活域Ψ1的常数矩阵Vα,1具体如下:
Vα,1=[I8×8 -I8×8 Λα,1 -Λα,1]T
Λα,1=[0 0 0 0 0 0 0 1]
其中,I8×8表示8乘8的单位矩阵;Λα,1表示时间无关灵活域Ψ1下的常数矩阵。
所述的时间无关灵活域Ψ1的参数矩阵U1具体如下:
其中,xmax和xmin分别表示决策变量的最大值和最小值;EES,max和EES,min分别表示电储能设备的最大和最小荷电状态。
所述的步骤3)中,根据时间无关灵活域和时间耦合灵活域,采用射线投影方法获取灵活域信息Μex,具体为使用射线投射方法来搜索时间无关灵活域Ψ1内的参考点,以时间无关灵活域Ψ1内任何一个点为起点,投射射线延该点右侧绘制并延伸到无限远,计算投射射线与时间无关灵活域Ψ1边缘相交的次数,若相交次数为奇数,则确定起点位于灵活域内并作为参考点;根据确定的参考点从数学上获取时间无关灵活域Ψ1的最终边界。
灵活域信息Μex具体如下:
其中,表示时间无关灵活域Ψ1内参考点和投影点之间的线段ab的外向法向量;和/>分别表示采用射线投影方法获取的投影点的横纵坐标。
时间无关灵活域Ψ1的最终边界的某一条边上的投影点可以通过做垂直线来找。与时间无关灵活域Ψ1,可以被描述为一个凸形的多边形,它限制了区域内的天然气和电力输入。关于与时间有关的灵活域Ψ2,它由一组线性方程表示,对应于储能设备的充电/放电约束。
所述的时间无关灵活域Ψ1内参考点和投影点之间的线段ab的外向法向量具体如下:
其中,和/>分别表示时间无关灵活域Ψ1内的参考点R的横纵坐标。
本发明的有益效果是:
本发明的方法为在多能耦合下实现综合能源枢纽灵活域评估的新方法,克服以往评估方法只考虑了独立时间间隔的不足;本发明方法基于可获取考虑时间关联特征的综合能源枢纽灵活域,可提升综合能源调度中心对能源枢纽的灵活资源的利用水平,提升系统的安全性,充分发挥多能替代,有利于节能减碳。
附图说明
图1为本发明方法的逻辑框图;
图2为本发明中电-气耦合系统和综合能源枢纽拓扑示意图;
图3的(a)为t=1时刻灵活域评估结果示意图;
图3的(b)为t=6时刻灵活域评估结果示意图;
图3的(c)为t=24时刻灵活域评估结果示意图。
具体实施方式
下面结合附图及具体实施例对本发明作进一步详细说明。
如图1所示,本发明综合能源枢纽灵活域评估方法包括如下步骤:
本发明综合能源枢纽灵活域评估方法包括如下步骤:
步骤1)建立考虑能源设备运行约束的综合能源枢纽的多能流模型,将效率参数矩阵输入多能流模型中,多能流模型输出连接关系参数矩阵。
步骤1)中,综合能源枢纽包括燃气锅炉、电储能设备、电热泵和热电联产机组,综合能源枢纽和综合能源系统以及终端能源负荷连接,综合能源系统提供电力和天然气能源,并通过综合能源枢纽为终端能源负荷提供热力、电力和天然气能源。
燃气锅炉是一种通过消耗天然气来产热的装置;电储能设备是指能将电力储存的装置;电热泵是消耗电能实现热能从低温到高温传递的装置;热电联产机组使用天然气作为燃料产生电力和热能的装置。基于综合能源系统提供的电力和天然气,满足终端能源负荷侧用户的电力、天然气、热等能源需求。
步骤1)中,建立的考虑能源设备运行约束的综合能源枢纽的多能流模型,具体如下:
[GLEH,PLEH]T=Ax
Bx=[HTL,GTL,PTL]T
其中,PLEH和GLEH分别表示从综合能源系统到综合能源枢纽的电力和天然气输入;A表示综合能源枢纽的效率参数矩阵;x表示综合能源枢纽的决策变量;B表示综合能源枢纽的连接关系参数矩阵;HTL、PTL和GTL分别表示从综合能源枢纽到终端能源负荷的热力、电力和天然气输出。
综合能源枢纽的决策变量x具体如下:
其中,HGB、HEHP和HCHP分别表示综合能源枢纽中燃气锅炉、电热泵和热电联产机组的热功率输出集合;PCHP和PES分别表示综合能源枢纽中热电联产机组和电储能设备的电功率输出集合;PN和GN分别表示终端能源负荷的电力负荷和天然气负荷;EES表示电储能设备的荷电状态集合。
综合能源枢纽的效率参数矩阵A具体如下:
其中,ηGB和ηEHP分别表示燃气锅炉和电热泵的转换效率;和/>分别表示热电联产机组的热能和电能输出效率。
第一行中的元素表示决策变量与综合能源枢纽的电力输入之间的关系。第二行中的元素表示决策变量和天然气输入之间的关系。
综合能源枢纽的连接关系参数矩阵B具体如下:
其中,连接关系参数矩阵B中的元素1表示综合能源枢纽中的不同装置之间存在连接关系,元素0表示综合能源枢纽中的不同装置之间无连接关系。
连接关系参数矩阵B中的不同行描述了决策变量和综合能源枢纽的多个能量输出之间的关系。
能源设备运行约束具体如下:
EES,t=EES,t-1-PES,t-1Δt
0≤HGB≤HGB,max
0≤HEHP≤HEHP,max
其中,EES,t和EES,t-1分别表示电储能设备在t时刻和t-1时刻的荷电状态集合;PES,t-1表示电储能设备在t-1时刻的充放电率集合;Δt表示时间间隔;EES,min和EES,max分别表示电储能设备的最小和最大荷电状态集合;表示电储能设备的最大充放电率集合;HGB,max、HEHP,max和HCHP,max分别表示燃气锅炉、热电联产机组和热电联产机组的最大热功率输出集合;PCHP,min和PCHP,max分别表示热电联产机组的最小和最大电功率输出集合。
能源设备运行约束为考虑到时间特征的不同能源设备的运行约束,综合能源枢纽中的一些设备的运行状态不仅受到当前时间状态的约束,还会受到其他时间运行状态的约束,如电储能。综合能源枢纽中的没有时间耦合约束的能源设备包括如燃气锅炉、电热泵和热电联产机组。
多能流模型的紧凑形式可以表示为:
Vαxt≤Vβxt-1+U
其中,Vα和Vβ分别表示第一和第二常数矩阵;xt和xt-1分别表示t时刻和t-1时刻的决策变量;U表示参数向量;I8×8和08×8分别表示的8乘8的单位矩阵和0矩阵;Λα和Λβ分别表示第一和第二常数矩阵;xmax和xmin分别表示决策变量的最大值和最小值;EES,max和EES,min分别表示电储能设备的最大和最小荷电状态。
将多能流划分为两部分,即时间无关部分和时间耦合部分,具体如下:
其中,GLEH和PLEH分别表示综合能源枢纽的天然气输入和电力输入;
A1是通过将时间耦合变量的列置零,A2是通过将时间无关变量的列置零,A1和A2均从A中通过矩阵变换方法计算出来。
步骤2)建立综合能源枢纽的灵活域模型,将连接关系参数矩阵输入灵活域模型中,灵活域模型输出综合能源枢纽的时间无关灵活域和时间耦合灵活域。
步骤2)中,建立的综合能源枢纽的灵活域模型具体如下:
其中,Ψ1表示综合能源枢纽的时间无关灵活域;表示闵可夫斯基求和;Ψ2表示综合能源枢纽的时间耦合灵活域;Vα,1表示时间无关灵活域Ψ1的常数矩阵;U1表示时间无关灵活域Ψ1的参数矩阵;GLEH,1和PLEH,1分别表示综合能源枢纽中时间无关灵活域Ψ1的天然气输入和电力输入;A1表示时间无关向量;GLEH,2和PLEH,2分别表示时间耦合灵活域Ψ2的天然气输入和电力输入;A2表示时间有关向量。
Ψ1是以天然气输入和电力输入为轴的二维区域,与综合能源枢纽的多能源替代密切相关。Ψ1中的能量流约束只保留了与时间无关的约束,由Vα,1和U1表示。Ψ2可以表征为与电力输入相对应的一维区间,与储能的充电/放电特征有关。关于时间耦合灵活域Ψ2,其时间耦合约束是指储能约束。因此,我们可以将Ψ2重新表述为一个储能系统模型。PES,t是由能源系统调度中心决定的调度指令的一部分,综合能源枢纽只需要提交容量限制(EES,min,EES,max)和最大充放电比例这些都是每个时间段的常数。PES,t和EES,t分别代表t时刻的电力输入和储能设备水平的向量。
与时间无关的灵活域Ψ1是一个高维多面体,对能量输入向量的投影是一个二维多边形Θ,其可代表综合能源枢纽与时间无关的灵活性,具体如下:
综合能源枢纽的灵活性可以表示为两个多面体的聚合,即能源系统调度中心可以获得有关灵活域的信息Ψ,用于调度和风险管理等。灵活域模型具体可表示如下:
上述公式表示一个高维多面体,其中电能存在时间耦合约束,这可能使综合能源枢纽在不同时间间隔的运行灵活性相互影响。基于t-1的状态,综合能源枢纽调度中心来评估综合能源枢纽在时间t的灵活性,而综合能源枢纽的运行状态是由能源系统调度中心在能源系统层面调度的。综合能源枢纽在一天中不同时间的灵活性评估和运行调度是一个滚动过程。
综合能源枢纽调度中心是指负责综合能源枢纽运行的机构,可对综合能源枢纽内部设备的运行状态进行调整。综合能源枢纽运行的灵活性主要是指从电-气耦合系统输入到综合能源枢纽的能量的调节能力,它受制于综合能源枢纽中的设备(如热电联产机组)的运行限制和终端负荷的能源需求。由于能源输入包含两种形式,即电力和燃气,因此灵活域被描述为一个二维空间。综合能源枢纽的电力和天然气输入可以在这个区域内进行调节,这可以为电-气耦合系统提供灵活资源。灵活域是指综合能源枢纽的电力输入和天然气输入的可调范围,可表示为以电力输入为横坐标、天然气输入为纵坐标的二维区域。电-气耦合系统调度中心是指负责综合能源枢纽实际调度运行的机构,可以调整综合能源枢纽中各个装置的运行状态。
时间无关灵活域Ψ1的常数矩阵Vα,1具体如下:
其中,I8×8表示8乘8的单位矩阵;Λα,1表示时间无关灵活域Ψ1下的常数矩阵。
时间无关灵活域Ψ1的参数矩阵U1具体如下:
其中,xmax和xmin分别表示决策变量的最大值和最小值;EES,max和EES,min分别表示电储能设备的最大和最小荷电状态。
步骤3)根据时间无关灵活域和时间耦合灵活域,采用射线投影方法获取灵活域信息,即获取综合能源枢纽的电力输入和天然气输入的可调范围,实现综合能源枢纽灵活域的评估。
步骤3)中,根据时间无关灵活域和时间耦合灵活域,采用射线投影方法获取灵活域信息Μex,具体为使用射线投射方法来搜索时间无关灵活域Ψ1内的参考点,以时间无关灵活域Ψ1内任何一个点为起点,投射射线延该点右侧绘制并延伸到无限远,计算投射射线与时间无关灵活域Ψ1边缘相交的次数,若相交次数为奇数,则确定起点位于灵活域内并作为参考点;根据确定的参考点从数学上获取时间无关灵活域Ψ1的最终边界。
灵活域信息Μex具体如下:
其中,表示时间无关灵活域Ψ1内参考点和投影点之间的线段ab的外向法向量;和/>分别表示采用射线投影方法获取的投影点的横纵坐标。
时间无关灵活域Ψ1的最终边界的某一条边上的投影点可以通过做垂直线来找。与时间无关灵活域Ψ1,可以被描述为一个凸形的多边形,它限制了区域内的天然气和电力输入。关于与时间有关的灵活域Ψ2,它由一组线性方程表示,对应于储能设备的充电/放电约束。
时间无关灵活域Ψ1内参考点和投影点之间的线段ab的外向法向量具体如下:
其中,和/>分别表示时间无关灵活域Ψ1内的参考点R的横纵坐标。
本发明的具体实施例如下:
以改进的IEEE33节点和比利时20节点组成的电-气耦合系统为例,结合技术方案和附图详细叙述本发明的具体实施。
电-气耦合系统包括电力系统和天然气系统。电力系统包括41个电力线路和9台发电机组组成,而天然气系统由17条管道、6个气源和2台压缩机组成。在电力系统中,节点5、8和13处的三台发电机被燃气机组取代,燃气机组分别由天然气系统中的天然气节点3、7和20供应。电力系统1、22、25节点新增3个风电机组。有4个综合能源枢纽接入天然气系统的天然气节点3、10、16和19,以及电力系统的电节点9、16、20和29。每个综合能源枢纽都包含电热泵,热电联产装置,燃气锅炉和电储能设备。拓扑图如图2所示。燃气锅炉和电热泵的转换效率分别设置为0.9和0.8。电储能的容量限制设置为6MW。天然气的热值为8.6kW/m3。
分析结果如图3所示。可以发现,综合能源枢纽灵活域以天然气输入为纵坐标,电力输入为横坐标,表示为一个二维的多边形。在此区域内,任一一点均可满足终端的负荷需求。
Claims (6)
1.一种考虑时间关联特性的综合能源枢纽灵活域评估方法,其特征在于:方法包括如下步骤:
步骤1)建立考虑能源设备运行约束的综合能源枢纽的多能流模型,将效率参数矩阵输入多能流模型中,多能流模型输出连接关系参数矩阵;
步骤2)建立综合能源枢纽的灵活域模型,将连接关系参数矩阵输入灵活域模型中,灵活域模型输出综合能源枢纽的时间无关灵活域和时间耦合灵活域;
步骤3)根据时间无关灵活域和时间耦合灵活域,采用射线投影方法获取灵活域信息,即获取综合能源枢纽的电力输入和天然气输入的可调范围,实现综合能源枢纽灵活域的评估;
所述的步骤1)中,综合能源枢纽包括燃气锅炉、电储能设备、电热泵和热电联产机组,综合能源枢纽和综合能源系统以及终端能源负荷连接,综合能源系统提供电力和天然气能源,并通过综合能源枢纽为终端能源负荷提供热力、电力和天然气能源;
所述的步骤1)中,建立的考虑能源设备运行约束的综合能源枢纽的多能流模型,具体如下:
其中,PLEH和GLEH分别表示从综合能源系统到综合能源枢纽的电力和天然气输入;A表示综合能源枢纽的效率参数矩阵;x表示综合能源枢纽的决策变量;B表示综合能源枢纽的连接关系参数矩阵;HTL、PTL和GTL分别表示从综合能源枢纽到终端能源负荷的热力、电力和天然气输出;
所述的能源设备运行约束具体如下:
EES,t=EES,t-1-PES,t-1Δt
0≤HGB≤HGB,max
0≤HEHP≤HEHP,max
其中,EES,t和EES,t-1分别表示电储能设备在t时刻和t-1时刻的荷电状态集合;PES,t-1表示电储能设备在t-1时刻的充放电率集合;Δt表示时间间隔;EES,min和EES,max分别表示电储能设备的最小和最大荷电状态集合;表示电储能设备的最大充放电率集合;HGB,max、HEHP,max和HCHP,max分别表示燃气锅炉、热电联产机组和热电联产机组的最大热功率输出集合;PCHP,min和PCHP,max分别表示热电联产机组的最小和最大电功率输出集合;
所述的步骤2)中,建立的综合能源枢纽的灵活域模型具体如下:
其中,Ψ1表示综合能源枢纽的时间无关灵活域;表示闵可夫斯基求和;Ψ2表示综合能源枢纽的时间耦合灵活域;Vα,1表示时间无关灵活域Ψ1的常数矩阵;U1表示时间无关灵活域Ψ1的参数矩阵;GLEH,1和PLEH,1分别表示综合能源枢纽中时间无关灵活域Ψ1的天然气输入和电力输入;A1表示时间无关向量;GLEH,2和PLEH,2分别表示时间耦合灵活域Ψ2的天然气输入和电力输入;A2表示时间有关向量。
2.根据权利要求1所述的一种考虑时间关联特性的综合能源枢纽灵活域评估方法,其特征在于:所述的综合能源枢纽的决策变量x具体如下:
其中,HGB、HEHP和HCHP分别表示综合能源枢纽中燃气锅炉、电热泵和热电联产机组的热功率输出集合;PCHP和PES分别表示综合能源枢纽中热电联产机组和电储能设备的电功率输出集合;PN和GN分别表示终端能源负荷的电力负荷和天然气负荷;EES表示电储能设备的荷电状态集合。
3.根据权利要求1所述的一种考虑时间关联特性的综合能源枢纽灵活域评估方法,其特征在于:所述的综合能源枢纽的效率参数矩阵A具体如下:
其中,ηGB和ηEHP分别表示燃气锅炉和电热泵的转换效率;和/>分别表示热电联产机组的热能和电能输出效率;
所述的综合能源枢纽的连接关系参数矩阵B具体如下:
其中,连接关系参数矩阵B中的元素1表示综合能源枢纽中的不同装置之间存在连接关系,元素0表示综合能源枢纽中的不同装置之间无连接关系。
4.根据权利要求1所述的一种考虑时间关联特性的综合能源枢纽灵活域评估方法,其特征在于:所述的时间无关灵活域Ψ1的常数矩阵Vα,1具体如下:
Λα,1=[0 0 0 0 0 0 0 1]
其中,I8×8表示8乘8的单位矩阵;Λα,1表示时间无关灵活域Ψ1下的常数矩阵;
所述的时间无关灵活域Ψ1的参数矩阵U1具体如下:
其中,xmax和xmin分别表示决策变量的最大值和最小值;EES,max和EES,min分别表示电储能设备的最大和最小荷电状态。
5.根据权利要求1所述的一种考虑时间关联特性的综合能源枢纽灵活域评估方法,其特征在于:所述的步骤3)中,根据时间无关灵活域和时间耦合灵活域,采用射线投影方法获取灵活域信息Μex,具体为使用射线投射方法来搜索时间无关灵活域Ψ1内的参考点,以时间无关灵活域Ψ1内任何一个点为起点,投射射线延该点右侧绘制并延伸,计算投射射线与时间无关灵活域Ψ1边缘相交的次数,若相交次数为奇数,则确定起点位于灵活域内并作为参考点;
灵活域信息Μex具体如下:
其中,表示时间无关灵活域Ψ1内参考点和投影点之间的线段ab的外向法向量;和/>分别表示采用射线投影方法获取的投影点的横纵坐标。
6.根据权利要求5所述的一种考虑时间关联特性的综合能源枢纽灵活域评估方法,其特征在于:所述的时间无关灵活域Ψ1内参考点和投影点之间的线段ab的外向法向量具体如下:
其中,和/>分别表示时间无关灵活域Ψ1内的参考点R的横纵坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310175488.0A CN116452009B (zh) | 2023-02-15 | 2023-02-15 | 一种考虑时间关联特性的综合能源枢纽灵活域评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310175488.0A CN116452009B (zh) | 2023-02-15 | 2023-02-15 | 一种考虑时间关联特性的综合能源枢纽灵活域评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116452009A CN116452009A (zh) | 2023-07-18 |
CN116452009B true CN116452009B (zh) | 2024-03-12 |
Family
ID=87130981
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310175488.0A Active CN116452009B (zh) | 2023-02-15 | 2023-02-15 | 一种考虑时间关联特性的综合能源枢纽灵活域评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116452009B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110322051A (zh) * | 2019-06-06 | 2019-10-11 | 国网浙江省电力有限公司经济技术研究院 | 考虑n-1安全约束的综合能源系统优化配置方法 |
WO2019200662A1 (zh) * | 2018-04-20 | 2019-10-24 | 东北大学 | 电-热-气综合能源系统的稳定评估与静态控制方法 |
WO2020237700A1 (zh) * | 2019-05-28 | 2020-12-03 | 国电南瑞科技股份有限公司 | 一种多能源系统运行调度方法 |
CN113537559A (zh) * | 2021-06-03 | 2021-10-22 | 浙江华云电力工程设计咨询有限公司 | 一种乡村综合能源规划资源优化系统及方法 |
-
2023
- 2023-02-15 CN CN202310175488.0A patent/CN116452009B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019200662A1 (zh) * | 2018-04-20 | 2019-10-24 | 东北大学 | 电-热-气综合能源系统的稳定评估与静态控制方法 |
WO2020237700A1 (zh) * | 2019-05-28 | 2020-12-03 | 国电南瑞科技股份有限公司 | 一种多能源系统运行调度方法 |
CN110322051A (zh) * | 2019-06-06 | 2019-10-11 | 国网浙江省电力有限公司经济技术研究院 | 考虑n-1安全约束的综合能源系统优化配置方法 |
CN113537559A (zh) * | 2021-06-03 | 2021-10-22 | 浙江华云电力工程设计咨询有限公司 | 一种乡村综合能源规划资源优化系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN116452009A (zh) | 2023-07-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110245878B (zh) | 智能楼宇群体的分布式综合能源需求响应协同优化方法 | |
Li et al. | A hybrid optimization-based scheduling strategy for combined cooling, heating, and power system with thermal energy storage | |
CN105447599B (zh) | 基于储热热电联产机组与电锅炉的弃风消纳协调调度模型 | |
CN110163443B (zh) | 考虑电-气综合需求响应的天然气调压站微能网优化调度方法 | |
Chen et al. | Economic and environmental operation of power systems including combined cooling, heating, power and energy storage resources using developed multi-objective grey wolf algorithm | |
CN111244939B (zh) | 一种计及需求侧响应的多能互补系统两级优化设计方法 | |
Xu et al. | Data mining based plant-level load dispatching strategy for the coal-fired power plant coal-saving: A case study | |
CN111950122A (zh) | 园区综合能源系统运行优化方法 | |
CN115186902A (zh) | 温室综合能源系统的调控方法、装置、终端及存储介质 | |
CN116822683A (zh) | 一种基于碳交易计及不确定性的综合能源系统优化运行方法 | |
CN116432824A (zh) | 基于多目标粒子群的综合能源系统优化方法及系统 | |
Bai et al. | Robust online operation of power systems with advanced adiabatic compressed air energy storage and renewable generation | |
CN115051388A (zh) | 一种基于分布鲁棒的“源-网-荷-储”两阶段调度优化方法 | |
CN111724026B (zh) | 一种多能源网络与配水网耦合运行的优化方法 | |
Varganova et al. | Energy-efficient distribution of heat between the boiler units of industrial power palnts with the use of computers | |
CN116452009B (zh) | 一种考虑时间关联特性的综合能源枢纽灵活域评估方法 | |
CN112214883A (zh) | 一种微能源网综合能源系统运行优化方法和系统 | |
CN111767621A (zh) | 一种基于知识迁移q学习算法的多能源系统优化调度方法 | |
CN111523947A (zh) | 一种虚拟电厂发电成本生成方法 | |
CN116957229A (zh) | 基于Hausdorff距离的微电网两阶段分布鲁棒优化调度方法 | |
CN114387129B (zh) | 考虑天然气管网的综合能源调度方法、装置、设备及介质 | |
CN115758763A (zh) | 一种计及源荷不确定性的多能流系统优化配置方法及系统 | |
CN213783243U (zh) | 一种面向产业园区的综合能源系统运行优化装置 | |
CN114936762A (zh) | 计及柔性电负荷的综合能源系统扩展规划方法 | |
CN115034453A (zh) | 考虑碳交易市场的综合能源系统灵活运行方法及计算平台 |
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 |