CN109447323B - 一种计及节点热价的综合能源系统两阶段容量配置方法 - Google Patents
一种计及节点热价的综合能源系统两阶段容量配置方法 Download PDFInfo
- Publication number
- CN109447323B CN109447323B CN201811154415.9A CN201811154415A CN109447323B CN 109447323 B CN109447323 B CN 109447323B CN 201811154415 A CN201811154415 A CN 201811154415A CN 109447323 B CN109447323 B CN 109447323B
- Authority
- CN
- China
- Prior art keywords
- heat
- node
- ith
- power
- stage
- 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 47
- 230000005611 electricity Effects 0.000 claims abstract description 14
- 238000012937 correction Methods 0.000 claims abstract description 11
- 239000007789 gas Substances 0.000 claims description 95
- 238000010438 heat treatment Methods 0.000 claims description 34
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 19
- 238000005485 electric heating Methods 0.000 claims description 16
- 238000003860 storage Methods 0.000 claims description 15
- 238000005338 heat storage Methods 0.000 claims description 14
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 14
- 239000002918 waste heat Substances 0.000 claims description 11
- 238000001816 cooling Methods 0.000 claims description 10
- 238000004519 manufacturing process Methods 0.000 claims description 10
- 238000011084 recovery Methods 0.000 claims description 10
- 230000005540 biological transmission Effects 0.000 claims description 9
- 239000003345 natural gas Substances 0.000 claims description 7
- 238000009434 installation Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000010248 power generation Methods 0.000 claims description 5
- 238000012546 transfer Methods 0.000 claims description 5
- 239000000446 fuel Substances 0.000 claims description 4
- 238000013461 design Methods 0.000 claims description 3
- 238000005286 illumination Methods 0.000 claims description 3
- 238000002347 injection Methods 0.000 claims description 3
- 239000007924 injection Substances 0.000 claims description 3
- 230000005855 radiation Effects 0.000 claims description 3
- 238000011144 upstream manufacturing Methods 0.000 claims description 3
- 238000007599 discharging Methods 0.000 claims 1
- 239000004576 sand Substances 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 4
- 238000005457 optimization Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 6
- 239000000243 solution Substances 0.000 description 3
- 238000004146 energy storage Methods 0.000 description 2
- 230000020169 heat generation Effects 0.000 description 2
- 239000008400 supply water Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- 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
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种计及节点热价的综合能源系统两阶段容量配置方法,包括建立区域热网运行约束模型;第一阶段,基于系统参数及负荷数据,建立包含多个热电联供子系统最优容量配置的目标函数及相应约束条件;建立区域热网温度修正方程,并基于潮流追踪法,建立节点热价方程;第二阶段,基于电价及节点热价,建立考虑需求侧响应的年购能费用目标函数及相应约束条件;建立综合能源系统两阶段规划方法的迭代求解过程,迭代求解最优配置容量。本发明方法首先提出一种节点热价的计算方法,并用于需求侧响应中,从而在规划阶段可以充分考虑系统优化运行及需求侧响应的影响。
Description
技术领域
本发明属于综合能源系统应用领域,具体来说是一种计及节点热价的综合能源系统两阶段容量配置方法。
背景技术
综合能源系统以电力系统为核心,结合供冷、供热、供气等子系统,以实现能源的梯级利用为基本原则,通过联供和储能设备进行各类能量的转化、分配与存储的多能流耦合一体化系统。而热电联产设备可利用发电余热进行制热,充分体现能源梯级利用的原则,其作为综合能源系统中的核心设备,可得到人们广泛的研究。在热电联产设备的容量配置方法,大部分文献主要针对单个热电联产设备进行容量配置研究,忽略了区域热网的协调调度作用。
发明内容
发明目的:为解决现有技术的不足,提供一种计及节点热价的综合能源系统的两阶段规划方法。
技术方案:为实现上述发明目的,本发明采用以下技术方案:
一种计及节点热价的综合能源系统两阶段容量配置方法,该方法包括以下步骤:
(1)建立区域热网运行约束模型;
(2)第一阶段,基于系统参数及负荷数据,建立包含多个热电联供子系统最优容量配置的目标函数及相应约束条件;
(3)建立区域热网温度修正方程组,并基于潮流追踪法,建立节点热价方程组;
(4)第二阶段,基于电价及节点热价,建立考虑需求侧响应的年购能费用目标函数及相应约束条件;
(5)建立综合能源系统两阶段规划方法的迭代求解过程,迭代求解综合能源系统的最优配置容量。
进一步的,步骤(1)包括:
(11)对于任意连接节点j,流入该节点的热媒流量之和等于流出该节点的热媒流量之和,区域热网包括热源节点、换热站节点和连接节点三种类型,具体运行约束方程如下:
(12)采用量调节的方式,事先设定好区域热网中各个管段的参考方向,则对于第i条管段在t时段的热媒流量,满足以下约束:
其中,和分别表示流向分别与管段参考方向相同和相反的热媒流量,mi,min和mi,max分别是管段允许流经热媒流量最小值和最大值,和均是0-1变量,用来表示管段中热媒流量的状态,等于1表示管段中存在该流向的热媒流量;
(13)热媒在传输过程中存在热量的损失,因此第i条管段在t时段的热功率在管段参考方向的始端与末端表示如下:
(14)第j个热源节点、换热站节点与区域热网交换热功率如下:
进一步的,步骤(2)包括:
(21)基于区域热网运行约束模型,建立综合能源系统第一阶段的最优容量配置及优化运行目标函数:
第i个热电联供子系统的年投资及运行费用如下:
第i个热电联供子系统中的第χ种设备的年投资费用表示如下:
其中,χ分别表示燃气轮机、燃气锅炉、蓄电池和蓄热槽,表示第χ种设备的安装容量,表示第χ种设备的单位投资费用,y表示第χ种设备的类型,为0-1变量,等于1表示第i个热电联供子系统安装了第y种类型的χ设备,j表示利率,n表示规划周期;
第i个热电联产子系统年运行成本包括燃气轮机和燃气锅炉所消耗的燃料成本:
其中,cgas是天然气的热值,是第y种类型的燃气轮机在t时段的输出电功率,是第y种类型的燃气锅炉在t时段的输出热功率,和是0-1变量,分别表示第i个热电联产子系统是否安装第y种类型的燃气轮机和燃气锅炉,和分别表示第y种类型的燃气轮机/燃气锅炉的效率,Igt和Igb分别表示备选燃气轮机和燃气锅炉的集合,Y是年规划周期,N是日调度周期,Δt是时间间隔;
(22),基于区域热网模型,建立综合能源系统第一阶段的最优容量配置及优化运行约束条件。
更进一步的,步骤(22)包括:
(a)热电功率平衡约束:
其中,和分别表示第i个热电联供子系统在t时段所产生的电功率和热功率,和分别表示第i个热电联供子系统在t时段换热器、燃气轮机、燃气锅炉、余热回收装置所产生的热功率以及蓄热槽存储和释放的热功率;和分别表示第i个热电联供子系统在t时段燃气轮机的电功率以及蓄电池的充放电功率;表示第i个热电联供子系统中燃气轮机的热电比;ηhe和ηhr分别表示换热器和余热回收装置的效率;
(b)设备约束方程:
进一步的,步骤(3)包括:
(31)建立区域热网的温度修正方程组:
(a)区域热网中第i条管段的始端温度和末端温度关系如下:
其中,τamb表示管段周围环境温度,λi和li分别表示第i条管段的导热系数和管段长度;
(b)对于任意连接节点j,从该节点流出的任意管段i具有相同的始端温度;
(c)对于任意连接节点j,该节点的热功率平衡方程如下:
(d)对于热源节点j,该节点的供水温度应等于设计好的供水温度;
其中,τs表示设计供水温度;
(32)建立区域热网的节点热价方程,对于任意节点k,该点的节点热价包括节点产热成本和节点传输成本:
(a)对于任意节点k,该点的节点产热成本表示如下:
其中,Qk,t表示注入节点k的热功率,是节点i所产生的热功率,从式(17)可知,对于热源节点,该点所产生的热功率等于热电联产子系统所产生的热功率对于其他节点,该点所产生的热功率为零;Au是顺流追踪矩阵,其第i行第j列元素如下:
对于第i个热电联供子系统,其单位产热成本表示如下:
其中,和分别表示第i个热电联供子系统在t时段内燃气轮机和燃气锅炉所产生的热功率;cgas表示天然气的单位热值价格;和分别表示第i个热电联供子系统在t时段内燃气轮机和燃气锅炉生成单位热功率的产热成本;和分别表示第i个热电联供子系统中燃气轮机和燃气锅炉的制热效率;
(b)对于任意节点k,该点的节点传输成本表示如下:
其中,Γp表示所有管段的集合,是节点i的热负荷功率,从式(21)可知,对于换热站节点,该点的热负荷功率等于换热站所消耗的热功率对于其他节点,该点的热负荷功率为零;Ad是逆流追踪矩阵,其第i行第j列元素如下:
其中,β是权重因子,用来调整单位热传输损耗成本的大小;
(c)因此,对于节点k,该点的节点热价表示如下:
进一步的,步骤(4)包括:
(41)基于需求侧响应,建立第二阶段年最优运行目标函数;
(a)第二阶段年运行费用为:
f2=min(Cel+Chl+Cpv) (25);
其中,Cel和Chl分别表示负荷集线器的年购电和购热费用,Cpv是光伏电池的年投资费用;
(b)负荷集线器的年购电和购热费用具体如下:
(c)光伏电池的年投资费用具体如下:
其中,cpv表示光伏电池的单位投资费用,ηpv表示光伏电池的效率,Apv表示光伏电池的安装面积,Θpv表示光伏电池的集合;
(42)建立设备约束方程:
(a)光伏电池在t时段产生的电功率满足以下约束方程:
其中,表示光伏电池在t时段所产生电功率用于供电;表示负荷集线器利用光伏所发电功率用于电制热,供给第i个换热站节点的热负荷;ppv表示光伏电池的单位额定功率;Apv表示光伏电池的总安装面积;ηpv表示光伏电池的效率;Itt表示光照辐射强度;表示负荷集线器通过电制热供给第i个换热站节点的最大热负荷;
(b)在需求侧响应中考虑负荷集线器采用光伏所发电功率进行电制热,因此对于第i个换热站节点的热负荷更新如下:
进一步的,步骤(5)包括:
(51))初始化热网温度,预测电热冷负荷数据;
(52)基于区域热网运行约束模型,结合负荷数据及各种设备参数,求解第一阶段的目标函数,获取热网中热媒流量及相应热功率和热损耗;
(53)基于区域热网温度修正方程组,节点热价方程,求解节点热/冷价;
(54)将节点热价作为输入变量,并结合需求侧响应方法,求解第二阶段的目标函数及并更新电热冷负荷数据;
(55)判断相邻两次迭代的温差是否小于误差要求,若达到要求则终止迭代,得到综合能源系统的最优配置容量;否则将更新的热网温度以及电热冷负荷数据返回步骤(52)继续迭代。
有益效果:与现有技术相比,本发明具有以下优点:
本发明的一种计及节点热价的综合能源系统的两阶段规划方法,首先建立区域热网运行约束模型;然后基于系统参数及负荷数据,建立第一阶段综合能源系统中各设备的最优容量配置的目标函数及相应约束条件;然后建立区域热网温度修正方程及节点热价方程;基于电价及节点热价,建立考虑需求侧响应的年购能费用目标函数及相应约束条件;最后,提出一种两阶段规划方法的迭代求解策略。与现有成果相比,本方法充分考虑的需求侧响应对设备容量配置及系统调度运行的影响,避免了设备配置过量造成的设备闲置问题,进一步提高的系统的经济性。此外,该成果所提出的计算节点热价的方法,可为以后热网的精细化供热提供理论依据。
附图说明
图1是综合能源系统结构图;
图2是典型热电联产系统结构图;
图3是综合能源系统两阶段规划示意图;
图4是迭代算法流程图;
图5是综合能源系统结构图;
图6是代表性节点的节点热价示意图。
具体实施方式
下面结合附图和具体实施例对本发明的技术方案进行详细说明。
如图1所示,综合能源系统包括多个热电联产系统,各个热电联产系统通过区域热网互联。
如图2所示,典型热电联产系统包括:热电联产设备(燃气轮机)、燃气锅炉、余热回收装置以及蓄电池和蓄热槽等储能设备,他们之间通过电网和气网连接;通过消耗天然气来产生电能和热能。
如图3所示,本发明基于区域热网,提出一种计及节点热价的综合能源系统的两阶段规划方法。第一阶段的目标函数为综合能源系统中年设备投资费用和运行费用最低;第二阶段的目标函数为考虑需求侧响应的年购能费用最低。此外,提出一种节点热价的计算方法,并应用于第二阶段的需求侧响应中。
具体包括以下步骤:
步骤10)建立区域热网运行约束模型;
步骤101)对于任意连接节点j,流入该节点的热媒流量之和等于流出该节点的热媒流量之和,区域热网包括热源节点、换热站节点和连接节点三种类型,具体运行约束方程如下:
步骤102)本实施例采用量调节的方式,因此需事先设定好区域热网中各个管段的参考方向,则对于第i条管段在t时段的热媒流量,满足以下约束:
其中,和分别表示流向分别与管段参考方向相同和相反的热媒流量,mi,min和mi,max分别是管段允许流经热媒流量最小值和最大值,和均是0-1变量,用来表示管段中热媒流量的状态,等于1表示管段中存在该流向的热媒流量。
步骤103)热媒在传输过程中存在热量的损失,因此第i条管段在t时段的热功率在管段参考方向的始端与末端表示如下:
步骤104)第j个热源节点/换热站节点与区域热网交换热功率如下:
步骤20)建立基于综合能源系统的最优容量配置及优化运行模型;
该步骤为第一阶段,基于系统参数及负荷数据,建立包含多个热电联供子系统最优容量配置的目标函数及相应约束条件;
步骤201),基于区域热网运行约束模型,建立综合能源系统第一阶段的最优容量配置及优化运行目标函数:
本发明考虑热电联供系统中燃气轮机、燃气锅炉、蓄电池和蓄热槽的配置费用,因此第i个热电联供子系统的年投资及运行费用如下:
第i个热电联供子系统中的第χ种设备的年投资费用表示如下:
其中,χ分别表示燃气轮机、燃气锅炉、蓄电池和蓄热槽,表示第χ种设备的安装容量,表示第χ种设备的单位投资费用,y表示第χ种设备的类型,为0-1变量,等于1表示第i个热电联供子系统安装了第y种类型的χ设备,j表示利率,n表示规划周期。
第i个热电联产子系统年运行成本包括燃气轮机和燃气锅炉所消耗的燃料成本:
其中,cgas是天然气的热值,是第y种类型的燃气轮机在t时段的输出电功率,是第y种类型的燃气锅炉在t时段的输出热功率,和是0-1变量,分别表示第i个热电联产子系统是否安装第y种类型的燃气轮机和燃气锅炉,和分别表示第y种类型的燃气轮机/燃气锅炉的效率,Igt和Igb分别表示备选燃气轮机和燃气锅炉的集合,Y是年规划周期,N是日调度周期,Δt是时间间隔。
步骤202),基于区域热网模型,建立综合能源系统第一阶段的最优容量配置及优化运行约束条件:
(1)热电功率平衡约束:
其中,和分别表示第i个热电联供子系统在t时段所产生的电功率和热功率,和分别表示第i个热电联供子系统在t时段换热器、燃气轮机、燃气锅炉、余热回收装置所产生的热功率以及蓄热槽存储和释放的热功率;和分别表示第i个热电联供子系统在t时段燃气轮机的电功率以及蓄电池的充放电功率;表示第i个热电联供子系统中燃气轮机的热电比;ηhe和ηhr分别表示换热器和余热回收装置的效率。
(2)设备约束方程:
步骤30)建立区域热网温度修正方程,并基于潮流追踪法,建立节点热价方程;
步骤301)建立区域热网的温度修正方程:
(1)区域热网中第i条管段的始段温度和末端温度关系如下:
其中,τamb表示管段周围环境温度,λi和li分别表示第i条管段的导热系数和管段长度。
(2)对于任意连接节点j,从该节点流出的任意管段i具有相同的始段温度。
(3)对于任意连接节点j,该节点的热功率平衡方程如下:
(4)对于热源节点j,该节点的供水温度应等于设计好的供水温度。
其中,τs表示设计供水温度,Θchp表示热源节点的集合。
步骤302)建立区域热网的节点热价方程,对于任意节点k,该点的节点热价包括节点产热成本和节点传输成本:
(1)对于任意节点k,该点的节点产热成本表示如下:
其中,Qk,t表示注入节点k的热功率,是节点i所产生的热功率,从式(17)可知,对于热源节点,该点所产生的热功率等于热电联产子系统所产生的热功率对于其他节点,该点所产生的热功率为零;Au是顺流追踪矩阵,其第i行第j列元素如下:
对于第i个热电联供子系统,其单位产热成本表示如下:
其中,和分别表示第i个热电联供子系统在t时段内燃气轮机和燃气锅炉所产生的热功率;cgas表示天然气的单位热值价格;和分别表示第i个热电联供子系统在t时段内燃气轮机和燃气锅炉生成单位热功率的产热成本;和分别表示第i个热电联供子系统中燃气轮机和燃气锅炉的制热效率。
(2)对于任意节点k,该点的节点传输成本表示如下:
其中,Γp表示所有管段的集合,是节点i的热负荷功率,从式(21)可知,对于换热站节点,该点的热负荷功率等于换热站所消耗的热功率对于其他节点,该点的热负荷功率为零;Ad是逆流追踪矩阵,其第i行第j列元素如下:
其中,β是权重因子,用来调整单位热传输损耗成本的大小。
(3)因此,对于节点k,该点的节点热价表示如下:
步骤40)建立第二阶段负荷集线器的年最优运行目标函数;
第二阶段,基于电价及节点热价,建立考虑需求侧响应的年购能费用目标函数及相应约束条件;
步骤401)基于需求侧响应,建立第二阶段年最优运行目标函数。
(1)第二阶段年运行费用
f2=min(Cel+Chl+Cpv) (25);
其中,Cel和Chl分别表示负荷集线器的年购电和购热费用,Cpv是光伏电池的年投资费用。
(2)负荷集线器的年购电和购热费用具体如下:
(3)光伏电池的年投资费用具体如下:
其中,cpv表示光伏电池的单位投资费用,ηpv表示光伏电池的效率,Apv表示光伏电池的安装面积,Θpv表示光伏电池的集合。
步骤402)建立设备约束方程:
(1)光伏电池在t时段产生的电功率满足以下约束方程:
其中,表示光伏电池在t时段所产生电功率用于供电;表示负荷集线器利用光伏所发电功率用于电制热,供给第i个换热站节点的热负荷;ppv表示光伏电池的单位额定功率;Apv表示光伏电池的总安装面积;ηpv表示光伏电池的效率;Itt表示光照辐射强度;表示负荷集线器通过电制热供给第i个换热站节点的最大热负荷。
(2)在需求侧响应中考虑负荷集线器采用光伏所发电功率进行电制热,因此对于第i个换热站节点的热负荷更新如下:
步骤50)建立综合能源系统两阶段规划方法的迭代求解过程,如图4所示,包括以下步骤:
步骤501)初始化k=0,初始化负荷数据和热网温度T(0),预测电热冷负荷数据;
步骤502)基于区域热网运行约束模型,结合负荷数据及各种设备参数,求解第一阶段的目标函数,获取热网中热媒流量m(k);
步骤503)基于区域热网温度修正方程组求解热网温度T"及相应热功率,基于节点热、冷价方程组,求解节点热价chl(k);
步骤504)将电价cel(k)和节点热价chl(k)作为输入变量,并结合需求侧响应方法,求解第二阶段的目标函数,并更新电热冷负荷数据;
步骤505)判断相邻两次迭代的温差δt是否小于误差要求,若达到要求则终止迭代,得到综合能源系统的最优配置容量;否则将更新的热网温度T(k)=T",以及更新的电热冷负荷数据,k=k+1,返回步骤502)继续迭代。
基于如图5所示某综合能源系统结构图,该系统包含2个CCHP子系统,25个换热站节点,以及50条热网管段,两个子系统分别在节点1和节点51。表1给出两种方案的规划结果,其中,方案一是不考虑节点热价以及需求侧响应的影响,方案二是本发明中的考虑节点热价及需求侧响应的两阶段规划方法。从规划结果可以看出,相比方案一,方案二可有效减少燃气轮机和燃气锅炉的配置容量。此外,图6给出CCHP1、CCHP2和代表性节点的节点热价随时间的变化特征。从图6中可以看出,本发明所提节点热价在时间和空间上具有明显的差异,为综合能源系统的精细化调度提供很好的支撑作用。
表1综合能源系统容量配置结果
本发明提出了一种综合能源系统的两阶段规划方法,该方法考虑了需求侧响应及用户侧屋顶光伏渗透的影响。该综合能源系统包含多个热电联供子系统(CCHP systems),各个子系统通过区域热网进行连接。第一阶段优化目标是将多个热电联供子系统年投资费用和运营费用降到最低,并获得最优的设备类型及容量。在此基础上,提出了一种节点热价格的计算方法,并利用到第二阶段的优化中。第二阶段的目标函数为负荷集线器(LA)的年购能总费用最低,将电价、节点热、价作为输入变量,将新的电、热负荷作为输出变量返回到第一阶段。在第二阶段的需求侧响应中,负荷集线器通过管理用户侧光伏发电进行电制热。本方法的有益效果在于:首先提出一种节点热价的计算方法,并用于需求侧响应中,从而在规划阶段可以充分考虑系统优化运行及需求侧响应的影响。
Claims (6)
1.一种计及节点热价的综合能源系统两阶段容量配置方法,其特征在于,该方法包括以下步骤:
(1)建立区域热网运行约束模型;具体包括以下步骤:
(11)对于任意节点i,流入该节点的热媒流量之和等于流出该节点的热媒流量之和,区域热网包括热源节点、换热站节点和连接节点三种类型,节点i为热源节点、换热站节点或连接节点,具体运行约束方程如下:
(12)采用量调节的方式,事先设定好区域热网中各个管段的参考方向,则对于第lij条管段在t时段的热媒流量,满足以下约束:
(13)热媒在传输过程中存在热量的损失,因此第lij条管段在t时段的热功率在管段参考方向的始端与末端表示如下:
(14)第is个热源节点、换热站节点ih与区域热网交换热功率如下:
其中,分别表示第is个热源节点和换热站节点ih在t时段与区域热网交换的热功率,Θchp和Θsb分别表示热源节点和换热站节点的集合,和分别表示与第is个热源节点/第ih个换热站节点相连的管段集合,第is个热电联供子系统等效为第is个热源节点,表示连接节点is和节点j的管段,表示连接节点ih和节点j的管段;
(2)第一阶段,基于系统参数及负荷数据,建立包含多个热电联供子系统最优容量配置的目标函数及相应约束条件;
(3)建立区域热网温度修正方程组,并基于潮流追踪法,建立节点热价方程组;
(4)第二阶段,基于电价及节点热价,建立考虑需求侧响应的年购能费用目标函数及相应约束条件;
(5)建立综合能源系统两阶段规划方法的迭代求解过程,迭代求解综合能源系统的最优配置容量。
2.根据权利要求1所述的一种计及节点热价的综合能源系统两阶段容量配置方法,其特征在于,步骤(2)包括:
(21)基于区域热网运行约束模型,建立综合能源系统第一阶段的最优容量配置及优化运行目标函数:
第is个热电联供子系统的年投资及运行费用如下:
第is个热电联供子系统中的第χ种设备的年投资费用表示如下:
其中,表示第is个热电联供子系统中的第χ种设备的年投资费用,χ分别表示燃气轮机、燃气锅炉、蓄电池和蓄热槽,表示第χ种设备的安装容量,表示第χ种设备的单位投资费用,y表示第χ种设备的类型,为0-1变量,等于1表示第is个热电联供子系统安装了第y种类型的χ设备,π表示利率,T表示规划周期;
第is个热电联供子系统年运行成本包括燃气轮机和燃气锅炉所消耗的燃料成本:
其中,表示第is个热电联供子系统年运行成本包括燃气轮机和燃气锅炉所消耗的燃料成本,cgas是天然气的热值,是第y种类型的燃气轮机在t时段的输出电功率,是第y种类型的燃气锅炉在t时段的输出热功率,和是0-1变量,分别表示第is个热电联供子系统是否安装第y种类型的燃气轮机和燃气锅炉,和分别表示第y种类型的燃气轮机/燃气锅炉的效率,Igt和Igb分别表示备选燃气轮机和燃气锅炉的集合,Y是年规划周期,N是日调度周期,Δt是时间间隔;
(22),基于区域热网模型,建立综合能源系统第一阶段的最优容量配置及优化运行约束条件。
3.根据权利要求2所述的一种计及节点热价的综合能源系统两阶段容量配置方法,其特征在于,步骤(22)包括:
(a)热电功率平衡约束:
其中,表示第is个热电联供子系统在t时段所产生的电功率, 和分别表示第is个热电联供子系统在t时段换热器、燃气轮机、燃气锅炉、余热回收装置所产生的热功率以及蓄热槽存储和释放的热功率;和分别表示第is个热电联供子系统在t时段燃气轮机的电功率以及蓄电池的充放电功率;表示第is个热电联供子系统中燃气轮机的热电比;ηhe和ηhr分别表示换热器和余热回收装置的效率;
(b)设备约束方程:
4.根据权利要求1所述的一种计及节点热价的综合能源系统两阶段容量配置方法,其特征在于,步骤(3)包括:
(31)建立区域热网的温度修正方程组:
(a)区域热网中第lij条管段的始端温度和末端温度关系如下:
(b)对于任意节点i,从该节点流出的任意管段lij具有相同的始端温度;
(c)对于任意节点i,该节点的热功率平衡方程如下:
(d)对于任意节点i,该节点的供水温度应等于设计好的供水温度;
其中,τs表示设计供水温度;
(32)建立区域热网的节点热价方程,对于任意节点i,该点的节点热价包括节点产热成本和节点传输成本:
(a)对于任意节点i,该点的节点产热成本表示如下:
其中,Qi,t表示注入节点i的热功率,是节点i所产生的热功率,从式(17)可知,对于热源节点,该点所产生的热功率等于热电联产子系统所产生的热功率对于其他节点,该点所产生的热功率为零;Au是顺流追踪矩阵,其第i行第is列元素如下:
对于第is个热电联供子系统,其单位产热成本表示如下:
其中,和分别表示第is个热电联供子系统在t时段内燃气轮机和燃气锅炉所产生的热功率;cgas表示天然气的单位热值价格;和分别表示第is个热电联供子系统在t时段内燃气轮机和燃气锅炉生成单位热功率的产热成本;和分别表示第is个热电联供子系统中燃气轮机和燃气锅炉的制热效率;表示第is个热电联供子系统中燃气轮机的热电比;
(b)对于任意节点i,该点的节点传输成本表示如下:
其中,Γp表示所有管段的集合,是节点i的热负荷功率,从式(21)可知,对于换热站节点,该点的热负荷功率等于换热站所消耗的热功率对于其他节点,该点的热负荷功率为零;Ad是逆流追踪矩阵,其第i行第ih列元素如下:
其中,β是权重因子,用来调整单位热传输损耗成本的大小;
(c)因此,对于节点i,该点的节点热价表示如下:
5.根据权利要求1所述的一种计及节点热价的综合能源系统两阶段容量配置方法,其特征在于,步骤(4)包括:
(41)基于需求侧响应,建立第二阶段年最优运行目标函数;
(a)第二阶段年运行费用为:
f2=min(Cel+Chl+Cpv) (25);
其中,f2表示第二阶段年运行费用;Cel和Chl分别表示负荷集线器的年购电和购热费用,Cpv是光伏电池的年投资费用;
(b)负荷集线器的年购电和购热费用具体如下:
其中,表示第is个热电联供子系统的发电成本,Δt为时间间隔,Y是年规划周期,N是日调度周期,表示第is个热电联供子系统在t时段所产生的电功率;表示换热站节点ih的用热成本;表示换热站节点ih的热功率;
(c)光伏电池的年投资费用具体如下:
其中,cpv表示光伏电池的单位投资费用,ηpv表示光伏电池的效率,Apv表示光伏电池的安装面积,Θpv表示光伏电池的集合,ip表示第ip个光伏电池节点,π表示利率,T表示规划周期;
(42)建立设备约束方程:
(a)光伏电池在t时段产生的电功率满足以下约束方程:
其中,表示光伏电池在t时段所产生电功率用于供电;表示负荷集线器利用光伏所发电功率用于电制热,供给第ih个换热站节点的热负荷;ppv表示光伏电池的单位额定功率;Apv表示光伏电池的总安装面积;ηpv表示光伏电池的效率;Itt表示光照辐射强度;表示负荷集线器通过电制热供给第ih个换热站节点的最大热负荷;
(b)在需求侧响应中考虑负荷集线器采用光伏所发电功率进行电制热,因此对于第ih个换热站节点的热负荷更新如下:
6.根据权利要求1所述的一种计及节点热价的综合能源系统两阶段容量配置方法,其特征在于,步骤(5)包括:
(51))初始化热网温度,预测电热冷负荷数据;
(52)基于区域热网运行约束模型,结合负荷数据及各种设备参数,求解第一阶段的目标函数,获取热网中热媒流量及相应热功率和热损耗;
(53)基于区域热网温度修正方程组,节点热价方程,求解节点热/冷价;
(54)将节点热价作为输入变量,并结合需求侧响应方法,求解第二阶段的目标函数及并更新电热冷负荷数据;
(55)判断相邻两次迭代的温差是否小于误差要求,若达到要求则终止迭代,得到综合能源系统的最优配置容量;否则将更新的热网温度以及电热冷负荷数据返回步骤(52)继续迭代。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811154415.9A CN109447323B (zh) | 2018-09-30 | 2018-09-30 | 一种计及节点热价的综合能源系统两阶段容量配置方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811154415.9A CN109447323B (zh) | 2018-09-30 | 2018-09-30 | 一种计及节点热价的综合能源系统两阶段容量配置方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109447323A CN109447323A (zh) | 2019-03-08 |
CN109447323B true CN109447323B (zh) | 2021-08-20 |
Family
ID=65545689
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811154415.9A Active CN109447323B (zh) | 2018-09-30 | 2018-09-30 | 一种计及节点热价的综合能源系统两阶段容量配置方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109447323B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110059853B (zh) * | 2019-03-11 | 2023-09-05 | 上海电气分布式能源科技有限公司 | 一种冷热电三联供系统中设备的配置方法 |
CN110288152B (zh) * | 2019-06-25 | 2022-08-26 | 河海大学 | 考虑电/热柔性负荷的区域综合能源系统储能配置方法 |
CN110336307A (zh) * | 2019-08-12 | 2019-10-15 | 国网天津市电力公司 | 基于热网蓄热特性电热综合能源系统风电消纳方法及装置 |
CN111815068B (zh) * | 2020-07-22 | 2023-01-06 | 福州大学 | 二阶段约束边界紧缩算法求解城市综合能源网络优化方法 |
CN112581112B (zh) * | 2020-12-04 | 2022-02-22 | 陕西启迪瑞行清洁能源研究院有限公司 | 一种分时间段的梯级式热费计量方法及装置 |
CN112669529B (zh) * | 2020-12-04 | 2021-12-17 | 陕西启迪瑞行清洁能源研究院有限公司 | 一种多梯级品位热能计费装置 |
CN112634525B (zh) * | 2020-12-04 | 2021-11-30 | 陕西启迪瑞行清洁能源研究院有限公司 | 一种基于回水温度及分时段的多梯级品位热能计费方法 |
CN113011655B (zh) * | 2021-03-22 | 2022-10-25 | 昆明理工大学 | 一种基于两阶段随机规划的社区多能源系统的混合储能系统规划方法 |
CN113704931A (zh) * | 2021-06-24 | 2021-11-26 | 中国电力科学研究院有限公司 | 考虑管网动态及设备变效率特性的综合能源系统规划方法 |
CN113487207B (zh) * | 2021-07-16 | 2022-06-28 | 重庆大学 | 基于多智能体系统最优用户集群的多能源社区中的多目标能源管理系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102331110A (zh) * | 2011-08-31 | 2012-01-25 | 北京中科华誉能源技术发展有限责任公司 | 基于吸收式换热的区域热电冷联合能源系统及其方法 |
CN106786753A (zh) * | 2016-12-29 | 2017-05-31 | 上海博翎能源科技有限公司 | 多用户的区域能源互联网的系统及其调节方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101089170B1 (ko) * | 2010-03-30 | 2011-12-02 | 서울과학기술대학교 산학협력단 | 지역난방 네트워크의 고급 프로세스 제어 방법 |
GB201218342D0 (en) * | 2012-10-12 | 2012-11-28 | Univ Leuven Kath | Method and system for distributing and/or controlling an energy flow taking into account constraints relating to the electricity network |
CN106055773B (zh) * | 2016-05-27 | 2019-05-21 | 东南大学 | 结合热网模型的多区域综合能源系统配置模型的建立方法 |
CN106527142B (zh) * | 2016-12-06 | 2019-08-27 | 国网江苏省电力公司徐州供电公司 | 一种有源配电网环境下的冷热电联供系统协调调度方法 |
CN106505634B (zh) * | 2016-12-14 | 2019-03-19 | 东南大学 | 基于两阶段协调优化与控制的冷热电联供型微网运行方法 |
-
2018
- 2018-09-30 CN CN201811154415.9A patent/CN109447323B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102331110A (zh) * | 2011-08-31 | 2012-01-25 | 北京中科华誉能源技术发展有限责任公司 | 基于吸收式换热的区域热电冷联合能源系统及其方法 |
CN106786753A (zh) * | 2016-12-29 | 2017-05-31 | 上海博翎能源科技有限公司 | 多用户的区域能源互联网的系统及其调节方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109447323A (zh) | 2019-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109447323B (zh) | 一种计及节点热价的综合能源系统两阶段容量配置方法 | |
Ma et al. | The optimal structure planning and energy management strategies of smart multi energy systems | |
CN109359839B (zh) | 一种综合能源系统中区域热网的节点热价计算方法 | |
Wang et al. | Synergistic planning of an integrated energy system containing hydrogen storage with the coupled use of electric-thermal energy | |
Li et al. | Power system planning based on CSP-CHP system to integrate variable renewable energy | |
CN106447529A (zh) | 考虑热水管网的分布式能源系统建模和运行优化方法 | |
CN106779471B (zh) | 一种多能互联交直流混合微电网系统及优化配置方法 | |
CN109389248A (zh) | 基于区域热网的综合能源系统热电功率协调调度方法 | |
CN114139837B (zh) | 一种考虑双层碳排放优化分配模型的区域多系统双层分散优化调度方法 | |
Gou et al. | Holistic analysis and optimization of distributed energy system considering different transport characteristics of multi-energy and component efficiency variation | |
CN117081143A (zh) | 促进分布式光伏就地消纳的园区综合能源系统协调优化运行方法 | |
CN112085263A (zh) | 一种用户侧分布式能源系统混合储能优化配置方法和系统 | |
Jia et al. | Multi-objective optimization study of regional integrated energy systems coupled with renewable energy, energy storage, and inter-station energy sharing | |
CN117371669B (zh) | 一种考虑碳交易风险成本的园区综合能源系统运行方法 | |
CN112200695A (zh) | 一种城市级综合能源系统优化调度方法 | |
Pasban-Gajan et al. | Optimal scheduling of renewable-based energy hubs considering time-of-use pricing scheme | |
CN115906456A (zh) | 考虑需求侧响应不确定性的含氢储能ies调度优化模型 | |
CN115860406A (zh) | 基于内部电价激励的园区综合能源系统能量调度方法 | |
CN212784787U (zh) | 一种污水厂可再生能源综合利用系统 | |
Chen et al. | Optimally coordinated operation of a multi-energy microgrid with coupled electrical and heat networks | |
CN116306194A (zh) | 基于多目标优化理论的综合能源系统能效优化配置方法 | |
CN114386256A (zh) | 一种考虑电热设备灵活性约束及热网特性的区域电热系统优化调度方法 | |
CN113762643A (zh) | 区域综合能源系统的储能容量优化配置方法 | |
Xi et al. | Multi-objective optimization strategy for regional multi-energy systems integrated with medium-high temperature solar thermal technology | |
Hao et al. | A Multi-Objective Optimization Planning of Park-level Integrated Energy System Based on Two-tier Planning Model |
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 | ||
TA01 | Transfer of patent application right | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20201203 Address after: 213300 room 428, building a, 218 Hongkou Road, Kunlun Street, Liyang City, Changzhou City, Jiangsu Province (in Zhongguancun Science and Technology Industrial Park, Jiangsu Province) Applicant after: Liyang Research Institute of Southeast University Applicant after: SOUTHEAST University Address before: 211189 Jiangsu Road, Jiangning District, Southeast University, No. 2, No. Applicant before: SOUTHEAST University |
|
GR01 | Patent grant | ||
GR01 | Patent grant |