CN109447323A - 一种计及节点热价的综合能源系统两阶段容量配置方法 - Google Patents

一种计及节点热价的综合能源系统两阶段容量配置方法 Download PDF

Info

Publication number
CN109447323A
CN109447323A CN201811154415.9A CN201811154415A CN109447323A CN 109447323 A CN109447323 A CN 109447323A CN 201811154415 A CN201811154415 A CN 201811154415A CN 109447323 A CN109447323 A CN 109447323A
Authority
CN
China
Prior art keywords
node
heat
indicate
pipeline section
caloric value
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.)
Granted
Application number
CN201811154415.9A
Other languages
English (en)
Other versions
CN109447323B (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.)
Southeast University
Liyang Research Institute of Southeast University
Original Assignee
Southeast University
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 Southeast University filed Critical Southeast University
Priority to CN201811154415.9A priority Critical patent/CN109447323B/zh
Publication of CN109447323A publication Critical patent/CN109447323A/zh
Application granted granted Critical
Publication of CN109447323B publication Critical patent/CN109447323B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06312Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
    • 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy 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,流入该节点的热媒流量之和等于流出该节点的热媒流量之和,区域热网包括热源节点、换热站节点和连接节点三种类型,具体运行约束方程如下:
其中,mi,t表示第i条管段在t时段的热媒流量,表示所有与连接节点j相连的管段的集合,Θin表示所有连接节点组成的集合;
(12)采用量调节的方式,事先设定好区域热网中各个管段的参考方向,则对于第i条管段在t时段的热媒流量,满足以下约束:
其中,分别表示流向分别与管段参考方向相同和相反的热媒流量,mi,min和mi,max分别是管段允许流经热媒流量最小值和最大值,均是0-1变量,用来表示管段中热媒流量的状态,等于1表示管段中存在该流向的热媒流量;
(13)热媒在传输过程中存在热量的损失,因此第i条管段在t时段的热功率在管段参考方向的始端与末端表示如下:
其中,c是水的比热容,分别是供水管段在管段参考方向始端与末端的热功率,分别是供水管段在管段参考方向始端与末端的热媒温度,τr是回水管段温度;
(14)第j个热源节点、换热站节点与区域热网交换热功率如下:
其中,分别表示第j个热源节点和换热站节点在t时段与区域热网交换的热功率,Θchp和Θsb分别表示热源节点和换热站节点的集合,分别表示第j个热源节点/换热站节点相连的管段集合。
进一步的,步骤(2)包括:
(21)基于区域热网运行约束模型,建立综合能源系统第一阶段的最优容量配置及优化运行目标函数:
其中,表示第i个热电联供子系统的年投资及运行费用,Θchp表示热电联供子系统的集合;
第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)设备约束方程:
其中,分别表示第i个热电联供子系统在t时段余热回收装置和换热器所产生的最大热功率,分别表示第i个热电联供子系统燃气轮机、燃气锅炉、蓄电池和蓄热槽的安装容量。
进一步的,步骤(3)包括:
(31)建立区域热网的温度修正方程组:
(a)区域热网中第i条管段的始端温度和末端温度关系如下:
其中,τamb表示管段周围环境温度,λi和li分别表示第i条管段的导热系数和管段长度;
(b)对于任意连接节点j,从该节点流出的任意管段i具有相同的始端温度;
其中,表示t时段内节点j处的热媒温度;
(c)对于任意连接节点j,该节点的热功率平衡方程如下:
(d)对于热源节点j,该节点的供水温度应等于设计好的供水温度;
其中,τs表示设计供水温度;
(32)建立区域热网的节点热价方程,对于任意节点k,该点的节点热价包括节点产热成本和节点传输成本:
(a)对于任意节点k,该点的节点产热成本表示如下:
其中,Qk,t表示注入节点k的热功率,是节点i所产生的热功率,从式(17)可知,对于热源节点,该点所产生的热功率等于热电联产子系统所产生的热功率对于其他节点,该点所产生的热功率为零;Au是顺流追踪矩阵,其第i行第j列元素如下:
其中,表示与节点i直接连接的上游节点组成的集合;
对于第i个热电联供子系统,其单位产热成本表示如下:
其中,分别表示第i个热电联供子系统在t时段内燃气轮机和燃气锅炉所产生的热功率;cgas表示天然气的单位热值价格;分别表示第i个热电联供子系统在t时段内燃气轮机和燃气锅炉生成单位热功率的产热成本;分别表示第i个热电联供子系统中燃气轮机和燃气锅炉的制热效率;
(b)对于任意节点k,该点的节点传输成本表示如下:
其中,Γp表示所有管段的集合,是节点i的热负荷功率,从式(21)可知,对于换热站节点,该点的热负荷功率等于换热站所消耗的热功率对于其他节点,该点的热负荷功率为零;Ad是逆流追踪矩阵,其第i行第j列元素如下:
其中,是表示与节点i直接连接的下游节点组成的集合。是管段j的单位热传输损耗成本,表示如下:
其中,β是权重因子,用来调整单位热传输损耗成本的大小;
(c)因此,对于节点k,该点的节点热价表示如下:
进一步的,步骤(4)包括:
(41)基于需求侧响应,建立第二阶段年最优运行目标函数;
(a)第二阶段年运行费用为:
f2=min(Cel+Chl+Cpv) (25);
其中,Cel和Chl分别表示负荷集线器的年购电和购热费用,Cpv是光伏电池的年投资费用;
(b)负荷集线器的年购电和购热费用具体如下:
其中,表示第i个热电联供子系统的发电成本,Δt为时间间隔;
(c)光伏电池的年投资费用具体如下:
其中,cpv表示光伏电池的单位投资费用,ηpv表示光伏电池的效率,Apv表示光伏电池的安装面积,Θpv表示光伏电池的集合;
(42)建立设备约束方程:
(a)光伏电池在t时段产生的电功率满足以下约束方程:
其中,表示光伏电池在t时段所产生电功率用于供电;表示负荷集线器利用光伏所发电功率用于电制热,供给第i个换热站节点的热负荷;ppv表示光伏电池的单位额定功率;Apv表示光伏电池的总安装面积;ηpv表示光伏电池的效率;Itt表示光照辐射强度;表示负荷集线器通过电制热供给第i个换热站节点的最大热负荷;
(b)在需求侧响应中考虑负荷集线器采用光伏所发电功率进行电制热,因此对于第i个换热站节点的热负荷更新如下:
其中,表示第i个换热站节点在t时段的初始预测热负荷,COPh,air表示中央空调制热系数。
进一步的,步骤(5)包括:
(51))初始化热网温度,预测电热冷负荷数据;
(52)基于区域热网运行约束模型,结合负荷数据及各种设备参数,求解第一阶段的目标函数,获取热网中热媒流量及相应热功率和热损耗;
(53)基于区域热网温度修正方程组,节点热价方程,求解节点热/冷价;
(54)将节点热价作为输入变量,并结合需求侧响应方法,求解第二阶段的目标函数及并更新电热冷负荷数据;
(55)判断相邻两次迭代的温差是否小于误差要求,若达到要求则终止迭代,得到综合能源系统的最优配置容量;否则将更新的热网温度以及电热冷负荷数据返回步骤(52)继续迭代。
有益效果:与现有技术相比,本发明具有以下优点:
本发明的一种计及节点热价的综合能源系统的两阶段规划方法,首先建立区域热网运行约束模型;然后基于系统参数及负荷数据,建立第一阶段综合能源系统中各设备的最优容量配置的目标函数及相应约束条件;然后建立区域热网温度修正方程及节点热价方程;基于电价及节点热价,建立考虑需求侧响应的年购能费用目标函数及相应约束条件;最后,提出一种两阶段规划方法的迭代求解策略。与现有成果相比,本方法充分考虑的需求侧响应对设备容量配置及系统调度运行的影响,避免了设备配置过量造成的设备闲置问题,进一步提高的系统的经济性。此外,该成果所提出的计算节点热价的方法,可为以后热网的精细化供热提供理论依据。
附图说明
图1是综合能源系统结构图;
图2是典型热电联产系统结构图;
图3是综合能源系统两阶段规划示意图;
图4是迭代算法流程图;
图5是综合能源系统结构图;
图6是代表性节点的节点热价示意图。
具体实施方式
下面结合附图和具体实施例对本发明的技术方案进行详细说明。
如图1所示,综合能源系统包括多个热电联产系统,各个热电联产系统通过区域热网互联。
如图2所示,典型热电联产系统包括:热电联产设备(燃气轮机)、燃气锅炉、余热回收装置以及蓄电池和蓄热槽等储能设备,他们之间通过电网和气网连接;通过消耗天然气来产生电能和热能。
如图3所示,本发明基于区域热网,提出一种计及节点热价的综合能源系统的两阶段规划方法。第一阶段的目标函数为综合能源系统中年设备投资费用和运行费用最低;第二阶段的目标函数为考虑需求侧响应的年购能费用最低。此外,提出一种节点热价的计算方法,并应用于第二阶段的需求侧响应中。
具体包括以下步骤:
步骤10)建立区域热网运行约束模型;
步骤101)对于任意连接节点j,流入该节点的热媒流量之和等于流出该节点的热媒流量之和,区域热网包括热源节点、换热站节点和连接节点三种类型,具体运行约束方程如下:
其中,mi,t表示第i条管段在t时段的热媒流量,表示所有与连接节点j相连的管段的集合,Θin表示所有连接节点组成的集合。
步骤102)本实施例采用量调节的方式,因此需事先设定好区域热网中各个管段的参考方向,则对于第i条管段在t时段的热媒流量,满足以下约束:
其中,分别表示流向分别与管段参考方向相同和相反的热媒流量,mi,min和mi,max分别是管段允许流经热媒流量最小值和最大值,均是0-1变量,用来表示管段中热媒流量的状态,等于1表示管段中存在该流向的热媒流量。
步骤103)热媒在传输过程中存在热量的损失,因此第i条管段在t时段的热功率在管段参考方向的始端与末端表示如下:
其中,c是水的比热容,分别是供水管段在管段参考方向始端与末端的热功率,分别是供水管段在管段参考方向始端与末端的热媒温度,τr是回水管段温度;在本专利中视为常数。
步骤104)第j个热源节点/换热站节点与区域热网交换热功率如下:
其中,分别表示第j个热源节点和换热站节点在t时段与区域热网交换的热功率,Θchp和Θsb分别表示热源节点和换热站节点的集合,分别表示第j个热源节点/换热站节点相连的管段集合。
步骤20)建立基于综合能源系统的最优容量配置及优化运行模型;
该步骤为第一阶段,基于系统参数及负荷数据,建立包含多个热电联供子系统最优容量配置的目标函数及相应约束条件;
步骤201),基于区域热网运行约束模型,建立综合能源系统第一阶段的最优容量配置及优化运行目标函数:
其中,表示第i个热电联供子系统的年投资及运行费用,Θchp表示热电联供子系统的集合。
本发明考虑热电联供系统中燃气轮机、燃气锅炉、蓄电池和蓄热槽的配置费用,因此第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)设备约束方程:
其中,分别表示第i个热电联供子系统在t时段余热回收装置和换热器所产生的最大热功率,分别表示第i个热电联供子系统燃气轮机、燃气锅炉、蓄电池和蓄热槽的安装容量。
步骤30)建立区域热网温度修正方程,并基于潮流追踪法,建立节点热价方程;
步骤301)建立区域热网的温度修正方程:
(1)区域热网中第i条管段的始段温度和末端温度关系如下:
其中,τamb表示管段周围环境温度,λi和li分别表示第i条管段的导热系数和管段长度。
(2)对于任意连接节点j,从该节点流出的任意管段i具有相同的始段温度。
其中,表示t时段内节点j处的热媒温度。
(3)对于任意连接节点j,该节点的热功率平衡方程如下:
(4)对于热源节点j,该节点的供水温度应等于设计好的供水温度。
其中,τs表示设计供水温度,Θchp表示热源节点的集合。
步骤302)建立区域热网的节点热价方程,对于任意节点k,该点的节点热价包括节点产热成本和节点传输成本:
(1)对于任意节点k,该点的节点产热成本表示如下:
其中,Qk,t表示注入节点k的热功率,是节点i所产生的热功率,从式(17)可知,对于热源节点,该点所产生的热功率等于热电联产子系统所产生的热功率对于其他节点,该点所产生的热功率为零;Au是顺流追踪矩阵,其第i行第j列元素如下:
其中,表示与节点i直接连接的上游节点组成的集合。
对于第i个热电联供子系统,其单位产热成本表示如下:
其中,分别表示第i个热电联供子系统在t时段内燃气轮机和燃气锅炉所产生的热功率;cgas表示天然气的单位热值价格;分别表示第i个热电联供子系统在t时段内燃气轮机和燃气锅炉生成单位热功率的产热成本;分别表示第i个热电联供子系统中燃气轮机和燃气锅炉的制热效率。
(2)对于任意节点k,该点的节点传输成本表示如下:
其中,Γp表示所有管段的集合,是节点i的热负荷功率,从式(21)可知,对于换热站节点,该点的热负荷功率等于换热站所消耗的热功率对于其他节点,该点的热负荷功率为零;Ad是逆流追踪矩阵,其第i行第j列元素如下:
其中,是表示与节点i直接连接的下游节点组成的集合。是管段j的单位热传输损耗成本,表示如下:
其中,β是权重因子,用来调整单位热传输损耗成本的大小。
(3)因此,对于节点k,该点的节点热价表示如下:
步骤40)建立第二阶段负荷集线器的年最优运行目标函数;
第二阶段,基于电价及节点热价,建立考虑需求侧响应的年购能费用目标函数及相应约束条件;
步骤401)基于需求侧响应,建立第二阶段年最优运行目标函数。
(1)第二阶段年运行费用
f2=min(Cel+Chl+Cpv) (25);
其中,Cel和Chl分别表示负荷集线器的年购电和购热费用,Cpv是光伏电池的年投资费用。
(2)负荷集线器的年购电和购热费用具体如下:
其中,表示第i个热电联供子系统的发电成本,Δt为时间间隔。
(3)光伏电池的年投资费用具体如下:
其中,cpv表示光伏电池的单位投资费用,ηpv表示光伏电池的效率,Apv表示光伏电池的安装面积,Θpv表示光伏电池的集合。
步骤402)建立设备约束方程:
(1)光伏电池在t时段产生的电功率满足以下约束方程:
其中,表示光伏电池在t时段所产生电功率用于供电;表示负荷集线器利用光伏所发电功率用于电制热,供给第i个换热站节点的热负荷;ppv表示光伏电池的单位额定功率;Apv表示光伏电池的总安装面积;ηpv表示光伏电池的效率;Itt表示光照辐射强度;表示负荷集线器通过电制热供给第i个换热站节点的最大热负荷。
(2)在需求侧响应中考虑负荷集线器采用光伏所发电功率进行电制热,因此对于第i个换热站节点的热负荷更新如下:
其中,表示第i个换热站节点在t时段的初始预测热负荷,COPh,air表示中央空调制热系数。
步骤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 (7)

1.一种计及节点热价的综合能源系统两阶段容量配置方法,其特征在于,该方法包括以下步骤:
(1)建立区域热网运行约束模型;
(2)第一阶段,基于系统参数及负荷数据,建立包含多个热电联供子系统最优容量配置的目标函数及相应约束条件;
(3)建立区域热网温度修正方程组,并基于潮流追踪法,建立节点热价方程组;
(4)第二阶段,基于电价及节点热价,建立考虑需求侧响应的年购能费用目标函数及相应约束条件;
(5)建立综合能源系统两阶段规划方法的迭代求解过程,迭代求解综合能源系统的最优配置容量。
2.根据权利要求1所述的一种计及节点热价的综合能源系统两阶段容量配置方法,其特征在于,步骤(1)包括:
(11)对于任意连接节点j,流入该节点的热媒流量之和等于流出该节点的热媒流量之和,区域热网包括热源节点、换热站节点和连接节点三种类型,具体运行约束方程如下:
其中,mi,t表示第i条管段在t时段的热媒流量,表示所有与连接节点j相连的管段的集合,Θin表示所有连接节点组成的集合;
(12)采用量调节的方式,事先设定好区域热网中各个管段的参考方向,则对于第i条管段在t时段的热媒流量,满足以下约束:
其中,分别表示流向分别与管段参考方向相同和相反的热媒流量,mi,min和mi ,max分别是管段允许流经热媒流量最小值和最大值,均是0-1变量,用来表示管段中热媒流量的状态,等于1表示管段中存在该流向的热媒流量;
(13)热媒在传输过程中存在热量的损失,因此第i条管段在t时段的热功率在管段参考方向的始端与末端表示如下:
其中,c是水的比热容,分别是供水管段在管段参考方向始端与末端的热功率,分别是供水管段在管段参考方向始端与末端的热媒温度,τr是回水管段温度;
(14)第j个热源节点、换热站节点与区域热网交换热功率如下:
其中,分别表示第j个热源节点和换热站节点在t时段与区域热网交换的热功率,Θchp和Θsb分别表示热源节点和换热站节点的集合,分别表示第j个热源节点/换热站节点相连的管段集合。
3.根据权利要求1所述的一种计及节点热价的综合能源系统两阶段容量配置方法,其特征在于,步骤(2)包括:
(21)基于区域热网运行约束模型,建立综合能源系统第一阶段的最优容量配置及优化运行目标函数:
其中,表示第i个热电联供子系统的年投资及运行费用,Θchp表示热电联供子系统的集合;
第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),基于区域热网模型,建立综合能源系统第一阶段的最优容量配置及优化运行约束条件。
4.根据权利要求3所述的一种计及节点热价的综合能源系统两阶段容量配置方法,其特征在于,步骤(22)包括:
(a)热电功率平衡约束:
其中,分别表示第i个热电联供子系统在t时段所产生的电功率和热功率,分别表示第i个热电联供子系统在t时段换热器、燃气轮机、燃气锅炉、余热回收装置所产生的热功率以及蓄热槽存储和释放的热功率;分别表示第i个热电联供子系统在t时段燃气轮机的电功率以及蓄电池的充放电功率;表示第i个热电联供子系统中燃气轮机的热电比;ηhe和ηhr分别表示换热器和余热回收装置的效率;
(b)设备约束方程:
其中,分别表示第i个热电联供子系统在t时段余热回收装置和换热器所产生的最大热功率,分别表示第i个热电联供子系统燃气轮机、燃气锅炉、蓄电池和蓄热槽的安装容量。
5.根据权利要求1所述的一种计及节点热价的综合能源系统两阶段容量配置方法,其特征在于,步骤(3)包括:
(31)建立区域热网的温度修正方程组:
(a)区域热网中第i条管段的始端温度和末端温度关系如下:
其中,τamb表示管段周围环境温度,λi和li分别表示第i条管段的导热系数和管段长度;
(b)对于任意连接节点j,从该节点流出的任意管段i具有相同的始端温度;
其中,表示t时段内节点j处的热媒温度;
(c)对于任意连接节点j,该节点的热功率平衡方程如下:
(d)对于热源节点j,该节点的供水温度应等于设计好的供水温度;
其中,τs表示设计供水温度;
(32)建立区域热网的节点热价方程,对于任意节点k,该点的节点热价包括节点产热成本和节点传输成本:
(a)对于任意节点k,该点的节点产热成本表示如下:
其中,Qk,t表示注入节点k的热功率,是节点i所产生的热功率,从式(17)可知,对于热源节点,该点所产生的热功率等于热电联产子系统所产生的热功率对于其他节点,该点所产生的热功率为零;Au是顺流追踪矩阵,其第i行第j列元素如下:
其中,表示与节点i直接连接的上游节点组成的集合;
对于第i个热电联供子系统,其单位产热成本表示如下:
其中,分别表示第i个热电联供子系统在t时段内燃气轮机和燃气锅炉所产生的热功率;cgas表示天然气的单位热值价格;分别表示第i个热电联供子系统在t时段内燃气轮机和燃气锅炉生成单位热功率的产热成本;分别表示第i个热电联供子系统中燃气轮机和燃气锅炉的制热效率;
(b)对于任意节点k,该点的节点传输成本表示如下:
其中,Γp表示所有管段的集合,是节点i的热负荷功率,从式(21)可知,对于换热站节点,该点的热负荷功率等于换热站所消耗的热功率对于其他节点,该点的热负荷功率为零;Ad是逆流追踪矩阵,其第i行第j列元素如下:
其中,是表示与节点i直接连接的下游节点组成的集合。是管段j的单位热传输损耗成本,表示如下:
其中,β是权重因子,用来调整单位热传输损耗成本的大小;
(c)因此,对于节点k,该点的节点热价表示如下:
6.根据权利要求1所述的一种计及节点热价的综合能源系统两阶段容量配置方法,其特征在于,步骤(4)包括:
(41)基于需求侧响应,建立第二阶段年最优运行目标函数;
(a)第二阶段年运行费用为:
f2=min(Cel+Chl+Cpv) (25);
其中,Cel和Chl分别表示负荷集线器的年购电和购热费用,Cpv是光伏电池的年投资费用;
(b)负荷集线器的年购电和购热费用具体如下:
其中,表示第i个热电联供子系统的发电成本,Δt为时间间隔;
(c)光伏电池的年投资费用具体如下:
其中,cpv表示光伏电池的单位投资费用,ηpv表示光伏电池的效率,Apv表示光伏电池的安装面积,Θpv表示光伏电池的集合;
(42)建立设备约束方程:
(a)光伏电池在t时段产生的电功率满足以下约束方程:
其中,表示光伏电池在t时段所产生电功率用于供电;表示负荷集线器利用光伏所发电功率用于电制热,供给第i个换热站节点的热负荷;ppv表示光伏电池的单位额定功率;Apv表示光伏电池的总安装面积;ηpv表示光伏电池的效率;Itt表示光照辐射强度;表示负荷集线器通过电制热供给第i个换热站节点的最大热负荷;
(b)在需求侧响应中考虑负荷集线器采用光伏所发电功率进行电制热,因此对于第i个换热站节点的热负荷更新如下:
其中,表示第i个换热站节点在t时段的初始预测热负荷,COPh,air表示中央空调制热系数。
7.根据权利要求1所述的一种计及节点热价的综合能源系统两阶段容量配置方法,其特征在于,步骤(5)包括:
(51))初始化热网温度,预测电热冷负荷数据;
(52)基于区域热网运行约束模型,结合负荷数据及各种设备参数,求解第一阶段的目标函数,获取热网中热媒流量及相应热功率和热损耗;
(53)基于区域热网温度修正方程组,节点热价方程,求解节点热/冷价;
(54)将节点热价作为输入变量,并结合需求侧响应方法,求解第二阶段的目标函数及并更新电热冷负荷数据;
(55)判断相邻两次迭代的温差是否小于误差要求,若达到要求则终止迭代,得到综合能源系统的最优配置容量;否则将更新的热网温度以及电热冷负荷数据返回步骤(52)继续迭代。
CN201811154415.9A 2018-09-30 2018-09-30 一种计及节点热价的综合能源系统两阶段容量配置方法 Active CN109447323B (zh)

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 true CN109447323A (zh) 2019-03-08
CN109447323B 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)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110059853A (zh) * 2019-03-11 2019-07-26 上海电气分布式能源科技有限公司 一种冷热电三联供系统中设备的配置方法
CN110288152A (zh) * 2019-06-25 2019-09-27 河海大学 考虑电/热柔性负荷的区域综合能源系统储能配置方法
CN110336307A (zh) * 2019-08-12 2019-10-15 国网天津市电力公司 基于热网蓄热特性电热综合能源系统风电消纳方法及装置
CN111815068A (zh) * 2020-07-22 2020-10-23 谢仕炜 二阶段约束边界紧缩算法求解城市综合能源网络优化方法
CN112581112A (zh) * 2020-12-04 2021-03-30 陕西启迪瑞行清洁能源研究院有限公司 一种分时间段的梯级式热费计量方法及装置
CN112634525A (zh) * 2020-12-04 2021-04-09 陕西启迪瑞行清洁能源研究院有限公司 一种基于回水温度及分时段的多梯级品位热能计费方法
CN112669529A (zh) * 2020-12-04 2021-04-16 陕西启迪瑞行清洁能源研究院有限公司 一种多梯级品位热能计费装置
CN113011655A (zh) * 2021-03-22 2021-06-22 昆明理工大学 一种基于两阶段随机规划的社区多能源系统的混合储能系统规划方法
CN113487207A (zh) * 2021-07-16 2021-10-08 重庆大学 基于多智能体系统最优用户集群的多能源社区中的多目标能源管理系统
CN113704931A (zh) * 2021-06-24 2021-11-26 中国电力科学研究院有限公司 考虑管网动态及设备变效率特性的综合能源系统规划方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20110109039A (ko) * 2010-03-30 2011-10-06 서울과학기술대학교 산학협력단 지역난방 네트워크의 고급 프로세스 제어 방법
CN102331110A (zh) * 2011-08-31 2012-01-25 北京中科华誉能源技术发展有限责任公司 基于吸收式换热的区域热电冷联合能源系统及其方法
CN106055773A (zh) * 2016-05-27 2016-10-26 东南大学 结合热网模型的多区域综合能源系统配置模型的建立方法
CN106505634A (zh) * 2016-12-14 2017-03-15 东南大学 基于两阶段协调优化与控制的冷热电联供型微网运行方法
CN106527142A (zh) * 2016-12-06 2017-03-22 国网江苏省电力公司徐州供电公司 一种有源配电网环境下的冷热电联供系统协调调度方法
CN106786753A (zh) * 2016-12-29 2017-05-31 上海博翎能源科技有限公司 多用户的区域能源互联网的系统及其调节方法
US10063052B2 (en) * 2012-10-12 2018-08-28 Katholieke Universiteit Leuven Method and system for distributing and/or controlling an energy flow taking into account constraints relating to the electricity network

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20110109039A (ko) * 2010-03-30 2011-10-06 서울과학기술대학교 산학협력단 지역난방 네트워크의 고급 프로세스 제어 방법
CN102331110A (zh) * 2011-08-31 2012-01-25 北京中科华誉能源技术发展有限责任公司 基于吸收式换热的区域热电冷联合能源系统及其方法
US10063052B2 (en) * 2012-10-12 2018-08-28 Katholieke Universiteit Leuven Method and system for distributing and/or controlling an energy flow taking into account constraints relating to the electricity network
CN106055773A (zh) * 2016-05-27 2016-10-26 东南大学 结合热网模型的多区域综合能源系统配置模型的建立方法
CN106527142A (zh) * 2016-12-06 2017-03-22 国网江苏省电力公司徐州供电公司 一种有源配电网环境下的冷热电联供系统协调调度方法
CN106505634A (zh) * 2016-12-14 2017-03-15 东南大学 基于两阶段协调优化与控制的冷热电联供型微网运行方法
CN106786753A (zh) * 2016-12-29 2017-05-31 上海博翎能源科技有限公司 多用户的区域能源互联网的系统及其调节方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
CHENYU WU 等: "Combined Economic Dispatch Considering the Time-Delay of District Heating Network and Multi-Regional Indoor Temperature Control", 《IEEE TRANSACTIONS ON SUSTAINABLE ENERGY》 *
顾伟 等: "多区域综合能源系统热网建模及系统运行优化", 《中国电机工程学报》 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110059853A (zh) * 2019-03-11 2019-07-26 上海电气分布式能源科技有限公司 一种冷热电三联供系统中设备的配置方法
CN110059853B (zh) * 2019-03-11 2023-09-05 上海电气分布式能源科技有限公司 一种冷热电三联供系统中设备的配置方法
CN110288152B (zh) * 2019-06-25 2022-08-26 河海大学 考虑电/热柔性负荷的区域综合能源系统储能配置方法
CN110288152A (zh) * 2019-06-25 2019-09-27 河海大学 考虑电/热柔性负荷的区域综合能源系统储能配置方法
CN110336307A (zh) * 2019-08-12 2019-10-15 国网天津市电力公司 基于热网蓄热特性电热综合能源系统风电消纳方法及装置
CN111815068A (zh) * 2020-07-22 2020-10-23 谢仕炜 二阶段约束边界紧缩算法求解城市综合能源网络优化方法
CN111815068B (zh) * 2020-07-22 2023-01-06 福州大学 二阶段约束边界紧缩算法求解城市综合能源网络优化方法
CN112581112A (zh) * 2020-12-04 2021-03-30 陕西启迪瑞行清洁能源研究院有限公司 一种分时间段的梯级式热费计量方法及装置
CN112581112B (zh) * 2020-12-04 2022-02-22 陕西启迪瑞行清洁能源研究院有限公司 一种分时间段的梯级式热费计量方法及装置
CN112669529A (zh) * 2020-12-04 2021-04-16 陕西启迪瑞行清洁能源研究院有限公司 一种多梯级品位热能计费装置
CN112634525A (zh) * 2020-12-04 2021-04-09 陕西启迪瑞行清洁能源研究院有限公司 一种基于回水温度及分时段的多梯级品位热能计费方法
CN113011655A (zh) * 2021-03-22 2021-06-22 昆明理工大学 一种基于两阶段随机规划的社区多能源系统的混合储能系统规划方法
CN113011655B (zh) * 2021-03-22 2022-10-25 昆明理工大学 一种基于两阶段随机规划的社区多能源系统的混合储能系统规划方法
CN113704931A (zh) * 2021-06-24 2021-11-26 中国电力科学研究院有限公司 考虑管网动态及设备变效率特性的综合能源系统规划方法
CN113487207A (zh) * 2021-07-16 2021-10-08 重庆大学 基于多智能体系统最优用户集群的多能源社区中的多目标能源管理系统
CN113487207B (zh) * 2021-07-16 2022-06-28 重庆大学 基于多智能体系统最优用户集群的多能源社区中的多目标能源管理系统

Also Published As

Publication number Publication date
CN109447323B (zh) 2021-08-20

Similar Documents

Publication Publication Date Title
CN109447323A (zh) 一种计及节点热价的综合能源系统两阶段容量配置方法
CN111445090B (zh) 一种离网型综合能源系统双层规划方法
CN111738502B (zh) 促进富余风电消纳的多能互补系统需求响应运行优化方法
CN106786753B (zh) 多用户的区域能源互联网的系统及其调节方法
CN105869075A (zh) 一种冷热电联供型微型能源网经济优化调度方法
CN110535128B (zh) 基于用能舒适度的多区域综合能源系统协同调度方法
CN110110913A (zh) 大型园区综合能源系统能源站优化配置方法
CN106447529A (zh) 考虑热水管网的分布式能源系统建模和运行优化方法
CN110458353A (zh) 电热联合微网能量梯级优化方法及系统
CN106779471A (zh) 一种多能互联交直流混合微电网系统及优化配置方法
CN109389248B (zh) 基于区域热网的综合能源系统热电功率协调调度方法
CN114077934B (zh) 一种综合能源微网互联系统及其调度方法
CN115170343A (zh) 一种区域综合能源系统分布式资源和储能协同规划方法
CN111553544B (zh) 基于一致性算法的工业园区分布式综合需求响应方法
CN111062598B (zh) 综合能源系统的分布式优化调度方法及系统
CN112085263A (zh) 一种用户侧分布式能源系统混合储能优化配置方法和系统
CN116308881A (zh) 利用供热管网储热的综合能源系统多时间尺度调度方法
CN117081143A (zh) 促进分布式光伏就地消纳的园区综合能源系统协调优化运行方法
CN112883630B (zh) 用于风电消纳的多微网系统日前优化经济调度方法
CN117035202B (zh) 考虑需求响应的电热综合能源系统双层协同扩展规划方法
CN111724026B (zh) 一种多能源网络与配水网耦合运行的优化方法
CN113128868A (zh) 一种区域综合能源系统调度优化方法及装置
CN116502921A (zh) 一种园区综合能源系统优化管理系统及其协调调度方法
CN114757469B (zh) 考虑多能流设备的区域综合能源系统日前优化调度方法
CN115906456A (zh) 考虑需求侧响应不确定性的含氢储能ies调度优化模型

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

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

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant