CN106055773A - 结合热网模型的多区域综合能源系统配置模型的建立方法 - Google Patents
结合热网模型的多区域综合能源系统配置模型的建立方法 Download PDFInfo
- Publication number
- CN106055773A CN106055773A CN201610362848.8A CN201610362848A CN106055773A CN 106055773 A CN106055773 A CN 106055773A CN 201610362848 A CN201610362848 A CN 201610362848A CN 106055773 A CN106055773 A CN 106055773A
- Authority
- CN
- China
- Prior art keywords
- represent
- website
- formula
- unit
- heat
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
Abstract
本发明公开了一种结合热网模型的多区域综合能源系统配置模型的建立方法,该方法包括以下步骤:步骤10)分别建立热网模型和CCHP模型;步骤20)结合建立的热网模型和CCHP模型,建立配置模型。该方法能够很好的利用各区域内冷热电负荷存在的明显的峰谷交错现象,实现多种能源的供需互补,降低容量配置,解决单独规划时出现的配置浪费、设备利用率低的问题,有效提升系统整体经济性。
Description
技术领域
本发明属于综合区域能源系统规划领域,具体来说,涉及一种结合热网模型的多区域综合能源系统配置模型的建立方法。
背景技术
能源危机和环境污染的双重压力促使人们对现有的能源消费模式进行反思,并开始对电、气、热等各种形式能源的综合利用进行研究。我国近期提出的“互联网+”智慧能源的行动计划,描绘了未来能源互联网愿景,指出充分考虑电、气、冷、热等不同形式能源在生产、传输、消费各环节的协同耦合是构建能源互联网的重要基础;打破各能源分开规划、独立运行的既有模式,发展在规划、运行、建设各阶段有机协调的综合能源系统是实现能源可持续发展的必经之路。
冷热电联供系统凭借其能效高、环境效益好等优势已成为综合能源系统的主要研究方向之一。目前对CCHP系统的研究的大多属于独立型CCHP系统,其规模小、终端用户品种少,对冷热电能的需求时间段单一,在规划或运行时,只能按本区域负荷峰值单独设计和建设,不可避免的导致设备能效低下,经济效益略有欠缺。结合热网(heating network,HN)的区域型CCHP系统正是解决上述问题的有效方法之一,其终端用户品种多样,各子区域内冷热电负荷存在明显的峰谷交错现象,运行时可协调各CCHP系统的热出力以匹配该区域热需求,实现能源的优化调度。
发明内容
技术问题:本发明所要解决的技术问题是:提出一种结合热网模型的多区域综合能源系统配置模型的建立方法,该方法能够很好的利用各区域内冷热电负荷存在的明显的峰谷交错现象,实现多种能源的供需互补,降低容量配置,解决单独规划时出现的配置浪费、设备利用率低的问题,有效提升系统整体经济性。
技术方案:为解决上述技术问题,本发明提出了一种结合热网模型的多区域综合能源系统配置模型的建立方法,该方法包括以下步骤:
步骤10)分别建立热网模型和CCHP模型;
步骤20)结合建立的热网模型和CCHP模型,建立配置模型。
作为优选例,所述的步骤10)中,建立热网模型的过程为:
步骤101)建立目标函数,具体包括步骤1011)—步骤1012):
步骤1011)建立管道年投资成本函数:
式中:Cpi,inv为管道年投资成本,R为管道成本年折算系数,Kpi,fix为管道固定费用,Nl为热网管道总数;Lk为热网中第k段管道的长度,单位:m;Kpi,var为与管道容量相关的可变费用,单位:元/kW·m;为第k段管道流过热能的最大值,单位:kW;
步骤1012)建立水泵运行费用函数:
式中:Cpu,ope为水泵运行费用;Nc为CCHP联供系统总数,ηehr,i为第i个CCHP系统中,水泵的耗电输热比,表示传输单位热量时耗费的电量,单位:kW;QHex,i,t为在t时刻,第i个CCHP系统与热网交互的热能,单位:kW;Ke,i,t为第i个CCHP系统t时刻购电单价,单位:元/kWh;△t为时间间隔,Nt表示全年时段数,Nt=8760/△t;
步骤102)建立约束条件,具体包括步骤1021)—步骤1022):
步骤1021)建立如式(4)所示的节点流量平衡约束条件:
[qi,t+∑qij,t]i=0 式(4)
式中:qi,t为t时刻第i个节点与第i个CCHP系统之间的流量;节点是CCHP系统与热网相连的点,节点与CCHP系统相对应设置;qij,t为t时刻与第i个节点相连接的各管段流量;
步骤1022)建立如式(5)所示的热能-流量约束条件:
qi,t=QHex,i,t/c(Ts-Tr) 式(5)
其中:c为水的比热容,Ts为给水温度,Tr为回水温度;
将式(5)带入式(4),得到式(6):
[QHex,i,t+∑QHij,t]i=0 式(6)
式中:QHij,t为在t时刻,与第i个节点相连接的各管段热量;
步骤1023)建立热损平衡约束:
设管道k与管道k-1之间的节点为节点i,QHk,t,1为流入管道k的热量,QHk,t,2为流出管道k的热量;由式(6)可知,对于节点i有:
QHex,i,t=QHk-1,t,2-QHk,t,1 式(7)
式中:QHk,t,1为t时刻流入管道k的热量,即为节点i流向下一节点j的热量,QHk-1,t,2为t时刻流出管道k-1的热量,即为上一节点i-1流向节点i的热量;QHex,i,t为在t时刻,第i个CCHP系统与热网交互的热能;
对于管段k有:
式中:QHk,t,1和QHk,t,2表示管道k两端的热能;δ为单位长度管道热损率;lk为第k段管道的长度。
作为优选例,所述的Δt=1h。
作为优选例,所述的步骤10)中,建立CCHP模型的过程为:
步骤111),建立目标函数,具体包括步骤1111)至步骤1112):
步骤1111):建立如式(9)所示的机组年投资成本函数:
式中:Ccchp,inv表示机组年投资成本;ΩGT表示燃气轮机型号集合,ΩGB表示燃气锅炉型号集合,ΩAC表示吸收式制冷机型号集合,ΩEC表示电制冷机型号集合;r1表示燃气轮机型号,r2表示燃气锅炉型号,r3表示吸收式制冷机型号,r4表示电制冷机型号;表示r1型燃气轮机的单位容量价格,表示r2型燃气锅炉的单位容量价格,表示r3型吸收式制冷机的单位容量价格,表示r4型电制冷机的单位容量价格,单位:元/kW;表示第i站点r1型燃气轮机的容量,表示第i站点r2型燃气锅炉的容量,表示第i站点r3型吸收式制冷机的容量,表示i站点r4型电制冷机的容量,单位:kW;表示在第i站点安装r1型燃气轮机,表示在第i站点安装r2型燃气锅炉,表示在第i站点安装r3型吸收式制冷机,表示在第i站点安装r4型电制冷机;
步骤1112):建立如式(10)所示的燃料费用函数:
式中:Cf,ope表示燃料费用,表示第i站点r1型燃气轮机在t时刻的发出的电量,单位:kW;ηgt,i为燃气轮机的发电效率,表示第i站点r2型燃气锅炉在t时刻发出的热能,单位:kW;ηgb,i为燃气锅炉的发热效率;Kf表示燃气单价,单位:元/m3;Hng为天然气的热值;
步骤1113):建立如式(11)所示的购电费用函数:
式中,Ce,ope表示系统全年运行电费,单位:元;Ke,i,t为第i个CCHP系统t时刻购电单价,单位:元/kWh;Pex,i,t表示第i站点t时刻与上级电网交互的电能,单位:kW;
步骤112)建立约束条件,具体包括步骤1121)至步骤1123):
步骤1121):建立如式(12)所示的冷热电功率平衡约束条件:
式中:Pd,i,t为第i个区域t时刻用户电负荷需求,表示第i站点r1型燃气轮机在t时刻的发出的电量,单位:kW,Pex,i,t表示第i站点t时刻与上级电网交互的电能,单位:kW;表示第i站点t时刻r4型电制冷机吸收的电量,QHd,i,t为第i区域t时刻用户热负荷需求,ηhr,i表示热回收器的回收效率,ηhe,i为热交换器效率,表示第i站点t时刻r3型吸收式制冷机机吸收的热量,QCd,i,t为第i区域t时刻用户冷负荷需求,COPac,i为第i站点吸收式制冷机的制冷效率,COPec,i为第i站点电制冷机的制冷效率;
步骤1122)建立设备容量约束条件:
建立如式(13)所示的燃气轮机与燃气锅炉容量与设备型号及其安装状态匹配条件:
其中,Wgt,i表示第i站点燃气轮机的总容量,Wgb,i表示第i站点燃气锅炉的总容量,Wac,i表示第i站点吸收式制冷机的总容量,Wec,i表示第i站点电制冷机的总容量;
步骤1123)建立如式(14)所示的运行上下限约束:
式中:u表示设备运行状态变量,表示第i站点t时刻r1型燃气轮机运行状态,当时,表示第i站点t时刻r1型燃气轮机运行;表示r1型燃气轮机功率下限,表示r1型燃气轮机功率上限;表示第i站点t时刻r2型燃气锅炉运行状态,当时,表示第i站点t时刻r2型燃气锅炉运行;为r2型燃气锅炉功率上限,为r2型燃气锅炉功率下限;表示第i站点t时刻r3型吸收式制冷机运行状态,当时,表示第i站点t时刻r3型吸收式制冷机运行;为r3型吸收式制冷机功率上限,为r3型吸收式制冷机功率下限;表示第i站点t时刻r4型电制冷机运行状态,当时,表示第i站点t时刻r4型电制冷机运行;表示r4型电制冷机功率上限,表示r4型电制冷机功率下限,为第i站点购电上限,为第i站点购电下限;为第i站点与热网交互上限,为第i站点与热网交互下限;
步骤1124)建立如式(15)所示的状态变量约束:
作为优选例,所述的步骤20)中,建立如式(16)所示的配置模型:
其中,ξ表示优化变量,其包含0-1变量与连续变量,f(ξ)为配置模型的目标函数,即式(1)、(3)、(9)—(11);h(ξ)为等式约束,即式(7)-(8)、(12)-(13);g(ξ)为不等式约束,即式(14)-(15);obj.表示目标函数;s.t.表示约束条件。
有益效果:与现有技术相比,本发明具有以下优点:
本发明实施例提出的一种结合热网模型的多区域综合能源系统配置模型的建立方法,首先针对各区域CCHP系统之间的环状热网,建立了考虑节点流量平衡、热能-流量约束及热损平衡约束的热网模型。结合CCHP系统能量平衡约束和热网模型,建立了多区域CCHP系统容量协同优化配置的混合整数线性规划模型。较各CCHP系统单独规划而言,本文提出的结合热网的多区域CCHP系统协同配置模型在规划阶段,提高了燃气轮机容量,降低了燃气锅炉容量,设备年投资成本较高,同时需要铺设热网管道,故在规划时经济性不明显;在运行阶段,由于考虑了热网,故能灵活调节各CCHP系统机组出力以匹配该区域的负荷,大幅降低运行费用,经济性十分显著。
由于多区域终端用户品种多样,结合热网后,能有效利用各子区域内冷热电负荷存在的明显的峰谷交错现象,通过热网转移负荷,在综合能源系统侧表现为有效提高系统运行时间,解决了单独规划时出现的燃气轮机利用率较低,燃气轮机配置容量浪费等问题。同时,在各区域热能交互时,本发明实施例提出的热网模型充分考虑了热网节点热量约束以及热网热损约束,故在无弃热的情况下,热能传输总是选择热损最小的路径,从而进一步降低系统运行费用。
附图说明
图1是本发明实施例中含热网的多区域综合能源系统结构图;
图2是本发明实施例的模型结构图;
图3是本发明实施例中简化的热网结构图;
图4是本发明实施例中热网管道热量分布图;
图5是本发明实施例采用的CCHP系统结构图;
图6是本发明实施案例的区域划分图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施案例对本发明进行深入地详细说明。应当理解,此处所描述的具体实施案例仅仅用以解释本发明,并不用于限定发明。
综合能源系统以冷热电联供系统为例,含热网的多区域综合能源系统结构如图1所示。每个子区域内设有一CCHP系统用以供给该区域的冷、热、电能。CCHP系统从燃气网获取燃气,燃烧产生热能与电能,部分热能与电能经过吸收式制冷与电制冷产生冷能,不足的电能由电网提供。每个CCHP系统可通过热网交互热能,第i个CCHP系统与热网交互的热能为QHex,i,t,但各CCHP系统之间不交互电能,电能只与上级电网交互。
如图2所示,本发明实施例的一种结合热网模型的多区域综合能源系统配置模型的建立方法,包括以下步骤:
步骤10)分别建立热网模型(即HN模型)和CCHP模型;
步骤20)结合建立的热网模型和CCHP模型,建立配置模型。
在上述实施例中,所述的步骤10)中,建立热网模型的过程为:
步骤101)建立目标函数,具体包括步骤1011)—步骤1012):
步骤1011)建立管道年投资成本函数:
式中:Cpi,inv为管道年投资成本。R为管道成本年折算系数。Kpi,fix为管道固定费用。Kpi,fix与管道容量无关,而与管道挖掘铺设费用相关,单位:元/m。Nl为热网管道总数;Lk为热网中第k段管道的长度,单位:m;Kpi,var为与管道容量相关的可变费用,单位:元/kW·m。为第k段管道流过热能的最大值,单位:kW。管道成本年折算系数R依据式(2)确定:
式中:y为投资回收年限,m为年利率。
步骤1012)建立水泵运行费用函数:
式中:Cpu,ope为水泵运行费用;Nc为CCHP联供系统总数,ηehr,i为第i个CCHP系统中,水泵的耗电输热比,表示传输单位热量时耗费的电量,单位:kW;QHex,i,t为在t时刻,第i个CCHP系统与热网交互的热能,单位:kW;Ke,i,t为第i个CCHP系统t时刻购电单价,单位:元/kWh;△t为时间间隔,Nt表示全年时段数,Nt=8760/△t。作为优选,Δt=1h。
步骤102)建立约束条件,具体包括步骤1021)—步骤1022)。
热网为各CCHP系统之间热网,属“源-源”网。热网对传输的安全可靠性要求极高,故优选采用环状网的布置形式。环状网在任一管道损坏时,闸阀可将它与其余管线隔开检修,安全可靠性较高,但同一区域,环状网的管线长度要长于树状网,造价也更高。为便于分析,将区域热网简化如图3。黑色实线条表示管道,黑色虚线表示热水流动正方向。图3主要体现管道中的热水走向。对于任一节点,均有三个方向的流量。以节点2为例,由上一节点流入的流量q12,t,流向CCHP系统的流量q2,t,流向下一节点的流量q23,t。其中,流向CCHP系统的流量q2,t所携带的热量为QHex,2,t。
步骤1021)对任何形式的网络,都须满足节点流量平衡约束,即对网络中任一节点,流向该节点的流量一定等于流出该节点的流量。建立如式(4)所示的节点流量平衡约束条件:
[qi,t+∑qij,t]i=0 式(4)
式中:qi,t为t时刻第i个节点与第i个CCHP系统之间的流量;节点是CCHP系统与热网相连的点,节点与CCHP系统相对应设置。即第1个CCHP系统对应第1个节点,CCHP系统的数量与节点的数量相等。qij,t为t时刻与第i个节点相连接的各管段流量。
步骤1022)建立如式(5)所示的热能-流量约束条件:
qi,t=QHex,i,t/c(Ts-Tr) 式(5)
其中:c为水的比热容。通常,c=4187J/(kg·℃)。Ts为给水温度,Tr为回水温度。
将式(5)带入式(4),得到式(6):
[QHex,i,t+∑QHij,t]i=0 式(6)
式中:QHij,t为在t时刻,与第i个节点相连接的各管段热量。
步骤1023)建立热损平衡约束:
由于热能在管网中传输不可避免的有热损,故热网模型还需包含热损平衡约束。由于存在热损,故对于某管道,其从一端流入的热能与从另一端流出的热能并不相同。如图4所示,设管道k与管道k-1之间的节点为节点i,QHk,t,1为流入管道k的热量,QHk,t,2为流出管道k的热量。
由式(6)可知,对于节点i有:
QHex,i,t=QHk-1,t,2-QHk,t,1 式(7)
式中:QHk,t,1为t时刻流入管道k的热量,即为节点i流向下一节点j的热量,QHk-1,t,2为t时刻流出管道k-1的热量,即为上一节点i-1流向节点i的热量;QHex,i,t为在t时刻,第i个CCHP系统与热网交互的热能;
对于管段k有:
式中:QHk,t,1和QHk,t,2表示管道k两端的热能;δ为单位长度管道热损率;lk为第k段管道的长度。设节点i的下一节点是节点j。若热水由节点i流向节点j,即QHk,t,1>0时,QHk,t,1>QHk,t,2,即ηk<1。若热水由节点j流入节点i,即QHk,t,1<0时,QHk,t,1<QHk,t,2,即ηk>1。
所述的步骤10)中,建立CCHP模型的过程为:
步骤111),建立目标函数,具体包括步骤1111)至步骤1112):
步骤1111):建立如式(9)所示的机组年投资成本函数:
式中:Ccchp,inv表示机组年投资成本;ΩGT表示燃气轮机型号集合,ΩGB表示燃气锅炉型号集合,ΩAC表示吸收式制冷机型号集合,ΩEC表示电制冷机型号集合;r1表示燃气轮机型号,r2表示燃气锅炉型号,r3表示吸收式制冷机型号,r4表示电制冷机型号;表示r1型燃气轮机的单位容量价格,表示r2型燃气锅炉的单位容量价格,表示r3型吸收式制冷机的单位容量价格,表示r4型电制冷机的单位容量价格,单位:元/kW;表示第i站点r1型燃气轮机的容量,表示第i站点r2型燃气锅炉的容量,表示第i站点r3型吸收式制冷机的容量,表示i站点r4型电制冷机的容量,单位:kW;表示在第i站点安装r1型燃气轮机,表示在第i站点安装r2型燃气锅炉,表示在第i站点安装r3型吸收式制冷机,表示在第i站点安装r4型电制冷机。
步骤1112):建立如式(10)所示的燃料费用函数:
式中:Cf,ope表示燃料费用,表示第i站点r1型燃气轮机在t时刻的发出的电量,单位:kW;ηgt,i为燃气轮机的发电效率,表示第i站点r2型燃气锅炉在t时刻发出的热能,单位:kW;ηgb,i为燃气锅炉的发热效率;Kf表示燃气单价,单位:元/m3;Hng为天然气的热值。
步骤1113):建立如式(11)所示的购电费用函数:
式中,Ce,ope表示系统全年运行电费,单位:元;Ke,i,t为第i个CCHP系统t时刻购电单价,单位:元/kWh;Pex,i,t表示第i站点t时刻与上级电网交互的电能,单位:kW。
步骤112)建立约束条件,具体包括步骤1121)至步骤1123):
如图5所示,HR为热回收器,HE为热交换器。第i个CCHP系统在t时刻,燃气轮机(GT)消耗燃气产生的电能为Pgt,i,t,产生的热能经热回收器(HR)回收的热能为QHgt,i,t;电能一部分直接供给i区域t时刻用户电需求Pd,i,t,另一部分Pec,i,t经过电制冷机(EC)产生冷能QCec,i,t,不足的电能由电网购电Pex,i,t;热能一部分经热交换器(HE)供给i区域t时刻用户热需求QHd,i,t,,另一部分QHac,i,t经吸收式制冷机(AC)产生冷能QCac,i,t,不足的热能由燃气锅炉(GB)燃烧燃气补充QHgb,i,t;i区域t时刻用户冷需求QCd,i,t由QCec,i,t与QCac,i,t组成。
步骤1121):建立如式(12)所示的冷热电功率平衡约束条件:
式中:Pd,i,t为第i个区域t时刻用户电负荷需求。本实施例为多区域综合能源系统协同配置模型,应用对象为一片综合型的用户区域,首先可将此综合型区域按照冷热电负荷特性分为若干个区域,每个区域配置一个CCHP系统,即第i个CCHP系统为第i个区域供冷热电能,其与环状热网相连的节点为第i个节点。表示第i站点r1型燃气轮机在t时刻的发出的电量,单位:kW;Pex,i,t表示第i站点t时刻与上级电网交互的电能,单位:kW;表示第i站点t时刻r4型电制冷机吸收的电量;QHd,i,t为第i区域t时刻用户热负荷需求,ηhr,i表示热回收器的回收效率;ηhe,i为热交换器效率;表示第i站点t时刻r3型吸收式制冷机机吸收的热量;QCd,i,t为第i区域t时刻用户冷负荷需求;COPac,i为第i站点吸收式制冷机的制冷效率;COPec,i为第i站点电制冷机的制冷效率。
步骤1122)建立设备容量约束条件:
建立如式(13)所示的燃气轮机与燃气锅炉容量与设备型号及其安装状态匹配条件:
其中,Wgt,i表示第i站点燃气轮机的总容量,Wgb,i表示第i站点燃气锅炉的总容量,Wac,i表示第i站点吸收式制冷机的总容量,Wec,i表示第i站点电制冷机的总容量。
步骤1123)建立如式(14)所示的运行上下限约束:
式中:u表示设备运行状态变量,表示第i站点t时刻r1型燃气轮机运行状态,当时,表示第i站点t时刻r1型燃气轮机运行;表示r1型燃气轮机功率下限,表示r1型燃气轮机功率上限;表示第i站点t时刻r2型燃气锅炉运行状态,当时,表示第i站点t时刻r2型燃气锅炉运行;为r2型燃气锅炉功率上限,为r2型燃气锅炉功率下限;表示第i站点t时刻r3型吸收式制冷机运行状态,当时,表示第i站点t时刻r3型吸收式制冷机运行;为r3型吸收式制冷机功率上限,为r3型吸收式制冷机功率下限;表示第i站点t时刻r4型电制冷机运行状态,当时,表示第i站点t时刻r4型电制冷机运行;表示r4型电制冷机功率上限,表示r4型电制冷机功率下限,为第i站点购电上限,为第i站点购电下限;为第i站点与热网交互上限,为第i站点与热网交互下限。
步骤1124)建立如式(15)所示的状态变量约束:
上述实施例中,所述的步骤20)中,建立如式(16)所示的配置模型:
其中,ξ表示优化变量,其包含0-1变量与连续变量,f(ξ)为配置模型的目标函数,即式(1)、(3)、(9)—(11);h(ξ)为等式约束,即式(7)-(8)、(12)-(13);g(ξ)为不等式约束,即式(14)-(15);obj.表示目标函数;s.t.表示约束条件。
该实施例中,配置模型的目标函数与约束均为线性,故此模型为典型的混合整数线性配置模型。该实施例中,建立CCHP模型的目的是确定在考虑热能交互时,设备的安装状态、最优容量配置及运行情况。需要配置的设备有燃气轮机(Gas Turbine,GT)、燃气锅炉(Gas Boiler,GB)、吸收式制冷机(Absorption Chiller,AC)、电制冷机(ElectricalChiller,EC)。建立HN模型的目的是确定热网中的实时流量,选取管道型号,实现水泵与闸阀的自动控制。
该实施例首先针对各区域CCHP系统之间的环状热网,建立了考虑节点流量平衡、热能-流量约束及热损平衡约束的热网模型。结合CCHP系统能量平衡约束和热网模型,建立了多区域CCHP系统容量协同优化配置的混合整数线性配置模型。本发明实施例在运行阶段,能够通过热网实现多区域综合能源系统的供需匹配,提升燃气轮机利用率,大幅降低运行费用,综合经济性十分显著。同时,利用热网模型充分考虑了热网节点热量约束以及热网热损约束,故在无弃热的情况下,热能传输总是选择热损最小的路径,从而进一步降低系统运行费用
下面例举一实施例。
以天津市某一综合区域为例,用户区域划分如图6所示。其中居民区、商业区、办公区与工业区的面积分别约为12万m2、12万m2、18万m2与30万m2,各区域CCHP系统之间热网管道长度依次为1.2km、0.65km、1.5km、0.6km。集中供热用天然气价格为2.37元/m3,居民电价0.48元/kWh,商业电价0.9元/kWh,工业电价0.71元/kWh。设定规划年限为10年,年利率为5%。本实施例以LINGO11.0求解结果为例。为便于比较分析,将未结合热网模型的独立规划模型作为对比例,将结合热网模型的协同规划模型作为本实施例。对比例和本实施例中,各CCHP系统优化配置结果及费用如表1所示。
表1
通过表1可以看出,应用本实施例模型后,低热电比用户(办公与商业)燃气轮机配置上升,高热电比用户(居民与工业)燃气锅炉配置下降。这是因为当考虑热网时,负荷热电比较低的用户在提高燃气轮机出力以满足自身电需求、降低购电量的同时,多余的热能可通过热网补贴其他负荷热电比较高的用户,故在本发明实施例下,办公与商业的燃气轮机配置增加。同时,热电比较高的用户在对比例中,不足的热能只能通过燃气锅炉补燃,而结合热网后,热能可由办公与商业区补给,从而减轻燃气锅炉的出力,降低燃气锅炉的配置。
本实施例相对于对比例节省了燃气锅炉投资成本,增加了燃气轮机投资成本,本实施例需考虑管道铺设费用,故总的设备年投资成本较对比例增加了768.71万元;本实施例中,商业与办公燃气轮机出力增大,燃气费用略有上升,但同时其购电费用显著降低,且考虑到商业与办公均属于电费较高的商业用电,故本实施例在运行时的经济性十分明显,年运行费用较对比例降低了2665.75万元。综上,本发明实施例运行成本低,经济性十分显著。
以上所述仅是本发明的优选实施方式,应该指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (5)
1.一种结合热网模型的多区域综合能源系统配置模型的建立方法,其特征在于,该方法包括以下步骤:
步骤10)分别建立热网模型和CCHP模型;
步骤20)结合建立的热网模型和CCHP模型,建立配置模型。
2.按照权利要求1所述的结合热网模型的多区域综合能源系统配置模型的建立方法,其特征在于,所述的步骤10)中,建立热网模型的过程为:
步骤101)建立目标函数,具体包括步骤1011)—步骤1012):
步骤1011)建立管道年投资成本函数:
式中:Cpi,inv为管道年投资成本,R为管道成本年折算系数,Kpi,fix为管道固定费用,Nl为热网管道总数;Lk为热网中第k段管道的长度,单位:m;Kpi,var为与管道容量相关的可变费用,单位:元/kW·m;为第k段管道流过热能的最大值,单位:kW;
步骤1012)建立水泵运行费用函数:
式中:Cpu,ope为水泵运行费用;Nc为CCHP联供系统总数,ηehr,i为第i个CCHP系统中,水泵的耗电输热比,表示传输单位热量时耗费的电量,单位:kW;QHex,i,t为在t时刻,第i个CCHP系统与热网交互的热能,单位:kW;Ke,i,t为第i个CCHP系统t时刻购电单价,单位:元/kWh;△t为时间间隔,Nt表示全年时段数,Nt=8760/△t;
步骤102)建立约束条件,具体包括步骤1021)—步骤1022):
步骤1021)建立如式(4)所示的节点流量平衡约束条件:
[qi,t+Σqij,t]i=0 式(4)
式中:qi,t为t时刻第i个节点与第i个CCHP系统之间的流量;节点是CCHP系统与热网相连的点,节点与CCHP系统相对应设置;qij,t为t时刻与第i个节点相连接的各管段流量;
步骤1022)建立如式(5)所示的热能-流量约束条件:
qi,t=QHex,i,t/c(Ts-Tr) 式(5)
其中:c为水的比热容,Ts为给水温度,Tr为回水温度;
将式(5)带入式(4),得到式(6):
[QHex,i,t+∑QHij,t]i=0 式(6)
式中:QHij,t为在t时刻,与第i个节点相连接的各管段热量;
步骤1023)建立热损平衡约束:
设管道k与管道k-1之间的节点为节点i,QHk,t,1为流入管道k的热量,QHk,t,2为流出管道k的热量;由式(6)可知,对于节点i有:
QHex,i,t=QHk-1,t,2-QHk,t,1 式(7)
式中:QHk,t,1为t时刻流入管道k的热量,即为节点i流向下一节点j的热量,QHk-1,t,2为t时刻流出管道k-1的热量,即为上一节点i-1流向节点i的热量;QHex,i,t为在t时刻,第i个CCHP系统与热网交互的热能。
对于管段k有:
式中:QHk,t,1和QHk,t,2表示管道k两端的热能;δ为单位长度管道热损率;lk为第k段管道的长度。
3.按照权利要求2所述的结合热网模型的多区域综合能源系统配置模型的建立方法,其特征在于,所述的Δt=1h。
4.按照权利要求1所述的结合热网模型的多区域综合能源系统配置模型的建立方法,其特征在于,所述的步骤10)中,建立CCHP模型的过程为:
步骤111),建立目标函数,具体包括步骤1111)至步骤1112):
步骤1111):建立如式(9)所示的机组年投资成本函数:
式中:Ccchp,inv表示机组年投资成本;ΩGT表示燃气轮机型号集合,ΩGB表示燃气锅炉型号集合,ΩAC表示吸收式制冷机型号集合,ΩEC表示电制冷机型号集合;r1表示燃气轮机型号,r2表示燃气锅炉型号,r3表示吸收式制冷机型号,r4表示电制冷机型号;表示r1型燃气轮机的单位容量价格,表示r2型燃气锅炉的单位容量价格,表示r3型吸收式制冷机的单位容量价格,表示r4型电制冷机的单位容量价格,单位:元/kW;表示第i站点r1型燃气轮机的容量,表示第i站点r2型燃气锅炉的容量,表示第i站点r3型吸收式制冷机的容量,表示i站点r4型电制冷机的容量,单位:kW;表示在第i站点安装r1型燃气轮机,表示在第i站点安装r2型燃气锅炉,表示在第i站点安装r3型吸收式制冷机,表示在第i站点安装r4型电制冷机;
步骤1112):建立如式(10)所示的燃料费用函数:
式中:Cf,ope表示燃料费用,表示第i站点r1型燃气轮机在t时刻的发出的电量,单位:kW;ηgt,i为燃气轮机的发电效率,表示第i站点r2型燃气锅炉在t时刻发出的热能,单位:kW;ηgb,i为燃气锅炉的发热效率;Kf表示燃气单价,单位:元/m3;Hng为天然气的热值;
步骤1113):建立如式(11)所示的购电费用函数:
式中,Ce,ope表示系统全年运行电费,单位:元;Ke,i,t为第i个CCHP系统t时刻购电单价,单位:元/kWh;Pex,i,t表示第i站点t时刻与上级电网交互的电能,单位:kW;
步骤112)建立约束条件,具体包括步骤1121)至步骤1123):
步骤1121):建立如式(12)所示的冷热电功率平衡约束条件:
式中:Pd,i,t为第i个区域t时刻用户电负荷需求,表示第i站点r1型燃气轮机在t时刻的发出的电量,单位:kW,Pex,i,t表示第i站点t时刻与上级电网交互的电能,单位:kW;表示第i站点t时刻r4型电制冷机吸收的电量,QHd,i,t为第i区域t时刻用户热负荷需求,ηhr,i表示热回收器的回收效率,ηhe,i为热交换器效率,表示第i站点t时刻r3型吸收式制冷机机吸收的热量,QCd,i,t为第i区域t时刻用户冷负荷需求,COPac,i为第i站点吸收式制冷机的制冷效率,COPec,i为第i站点电制冷机的制冷效率;
步骤1122)建立设备容量约束条件:
建立如式(13)所示的燃气轮机与燃气锅炉容量与设备型号及其安装状态匹配条件:
其中,表示第i站点燃气轮机的总容量,表示第i站点燃气锅炉的总容量,表示第i站点吸收式制冷机的总容量,表示第i站点电制冷机的总容量;
步骤1123)建立如式(14)所示的运行上下限约束:
式中:u表示设备运行状态变量,表示第i站点t时刻r1型燃气轮机运行状态,当时,表示第i站点t时刻r1型燃气轮机运行;表示r1型燃气轮机功率下限,表示r1型燃气轮机功率上限;表示第i站点t时刻r2型燃气锅炉运行状态,当时,表示第i站点t时刻r2型燃气锅炉运行;为r2型燃气锅炉功率上限,为r2型燃气锅炉功率下限;表示第i站点t时刻r3型吸收式制冷机运行状态,当时,表示第i站点t时刻r3型吸收式制冷机运行;为r3型吸收式制冷机功率上限,为r3型吸收式制冷机功率下限;表示第i站点t时刻r4型电制冷机运行状态,当时,表示第i站点t时刻r4型电制冷机运行;表示r4型电制冷机功率上限,表示r4型电制冷机功率下限,为第i站点购电上限,为第i站点购电下限;为第i站点与热网交互上限,为第i站点与热网交互下限;
步骤1124)建立如式(15)所示的状态变量约束:
5.按照权利要求1所述的结合热网模型的多区域综合能源系统配置模型的建立方法,其特征在于,所述的步骤20)中,建立如式(16)所示的配置模型:
其中,ξ表示优化变量,其包含0-1变量与连续变量,f(ξ)为配置模型的目标函数,即式(1)、(3)、(9)—(11);h(ξ)为等式约束,即式(7)-(8)、(12)-(13);g(ξ)为不等式约束,即式(14)-(15);obj.表示目标函数;s.t.表示约束条件。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610362848.8A CN106055773B (zh) | 2016-05-27 | 2016-05-27 | 结合热网模型的多区域综合能源系统配置模型的建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610362848.8A CN106055773B (zh) | 2016-05-27 | 2016-05-27 | 结合热网模型的多区域综合能源系统配置模型的建立方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106055773A true CN106055773A (zh) | 2016-10-26 |
CN106055773B CN106055773B (zh) | 2019-05-21 |
Family
ID=57175461
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610362848.8A Active CN106055773B (zh) | 2016-05-27 | 2016-05-27 | 结合热网模型的多区域综合能源系统配置模型的建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106055773B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107642772A (zh) * | 2017-09-11 | 2018-01-30 | 国家电网公司 | 冷热电联产系统同时满足负荷需求运行控制方法 |
CN108197768A (zh) * | 2018-04-02 | 2018-06-22 | 厦门大学 | 一种能源系统与管网布局联合优化方法 |
CN108694483A (zh) * | 2018-07-28 | 2018-10-23 | 天津大学 | 考虑电热气耦合的综合能源系统设备选型及容量规划方法 |
CN109447323A (zh) * | 2018-09-30 | 2019-03-08 | 东南大学 | 一种计及节点热价的综合能源系统两阶段容量配置方法 |
CN109449918A (zh) * | 2018-09-13 | 2019-03-08 | 北京国电通网络技术有限公司 | 区域内综合能源系统的构建方法和规划装置 |
CN110110904A (zh) * | 2019-04-17 | 2019-08-09 | 华北电力大学 | 考虑经济性、独立性和碳排放的综合能源系统优化方法 |
CN110188492A (zh) * | 2019-06-04 | 2019-08-30 | 南通大学 | 一种考虑热网特性的冷热电联供微网优化调度方法 |
CN110619425A (zh) * | 2019-08-06 | 2019-12-27 | 国网山东省电力公司经济技术研究院 | 一种考虑源网荷储差异特性的多功能区综合能源系统协同规划方法 |
CN110968827A (zh) * | 2019-11-13 | 2020-04-07 | 国家电网有限公司 | 一种多区域综合能源系统优化配置方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102829575A (zh) * | 2012-09-10 | 2012-12-19 | 陈戈 | 多源区域型冷热电联供能源网系统及方法 |
CN105576710A (zh) * | 2016-02-18 | 2016-05-11 | 东南大学 | 一种综合能源系统中的分布式电源的配置方法 |
-
2016
- 2016-05-27 CN CN201610362848.8A patent/CN106055773B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102829575A (zh) * | 2012-09-10 | 2012-12-19 | 陈戈 | 多源区域型冷热电联供能源网系统及方法 |
CN105576710A (zh) * | 2016-02-18 | 2016-05-11 | 东南大学 | 一种综合能源系统中的分布式电源的配置方法 |
Non-Patent Citations (3)
Title |
---|
HAICHAO WANG,ET AL.: ""Modeling and optimization of CHP based district heating system with renewable energy production and energy storage"", 《APPLIED ENERGY》 * |
王成山 等: ""冷热电联供微网优化调度通用建模方法"", 《中国电机工程学报》 * |
王进 等: ""与电力系统协同区域型分布式冷热电联供能源系统集成方案"", 《电力系统自动化》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107642772A (zh) * | 2017-09-11 | 2018-01-30 | 国家电网公司 | 冷热电联产系统同时满足负荷需求运行控制方法 |
CN108197768A (zh) * | 2018-04-02 | 2018-06-22 | 厦门大学 | 一种能源系统与管网布局联合优化方法 |
CN108197768B (zh) * | 2018-04-02 | 2021-07-27 | 厦门大学 | 一种能源系统与管网布局联合优化方法 |
CN108694483A (zh) * | 2018-07-28 | 2018-10-23 | 天津大学 | 考虑电热气耦合的综合能源系统设备选型及容量规划方法 |
CN109449918A (zh) * | 2018-09-13 | 2019-03-08 | 北京国电通网络技术有限公司 | 区域内综合能源系统的构建方法和规划装置 |
CN109447323A (zh) * | 2018-09-30 | 2019-03-08 | 东南大学 | 一种计及节点热价的综合能源系统两阶段容量配置方法 |
CN110110904A (zh) * | 2019-04-17 | 2019-08-09 | 华北电力大学 | 考虑经济性、独立性和碳排放的综合能源系统优化方法 |
CN110188492A (zh) * | 2019-06-04 | 2019-08-30 | 南通大学 | 一种考虑热网特性的冷热电联供微网优化调度方法 |
CN110619425A (zh) * | 2019-08-06 | 2019-12-27 | 国网山东省电力公司经济技术研究院 | 一种考虑源网荷储差异特性的多功能区综合能源系统协同规划方法 |
CN110619425B (zh) * | 2019-08-06 | 2022-03-04 | 国网山东省电力公司经济技术研究院 | 一种考虑源网荷储差异特性的多功能区综合能源系统协同规划方法 |
CN110968827A (zh) * | 2019-11-13 | 2020-04-07 | 国家电网有限公司 | 一种多区域综合能源系统优化配置方法 |
CN110968827B (zh) * | 2019-11-13 | 2023-06-27 | 国家电网有限公司 | 一种多区域综合能源系统优化配置方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106055773B (zh) | 2019-05-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106055773B (zh) | 结合热网模型的多区域综合能源系统配置模型的建立方法 | |
Söderman et al. | Structural and operational optimisation of distributed energy systems | |
Marques et al. | Thermodynamic analysis of tri-generation systems taking into account refrigeration, heating and electricity load demands | |
CN108960503B (zh) | 基于内点法的综合能源系统多场景优化分析方法 | |
CN106447113B (zh) | 一种基于运行优化模型的多区域综合能源系统运行方法 | |
CN103728881B (zh) | 一种多楼宇冷热电联供系统的优化运行方法 | |
WO2011134432A1 (zh) | 实现能源优化利用的泛能网及提供能源交易和服务的方法 | |
CN102359739A (zh) | 零能源损耗率热电厂的燃蒸循环热电冷三联供系统与方法 | |
CN111899120A (zh) | 一种综合能源规划与运行联合优化模型建立方法 | |
CN202267113U (zh) | 一种零能源损耗率的热电厂的燃蒸循环热电冷三联供系统 | |
CN110968827B (zh) | 一种多区域综合能源系统优化配置方法 | |
CN107784382A (zh) | 基于能源路由器的用户侧能源互联网规划方法 | |
CN105258384A (zh) | 一种集成热化学过程的热电冷多联产系统 | |
CN110502791B (zh) | 基于能源集线器的综合能源系统稳态建模方法 | |
CN110535128A (zh) | 基于用能舒适度的多区域综合能源系统协同调度方法 | |
CN109472413A (zh) | 考虑热管网传输特性的园区综合能源系统优化调度方法 | |
CN106940033A (zh) | 基于弃风电能的多设备联合高低温独立蓄热供热系统 | |
WO2017188083A1 (ja) | マイクログリッドの運転計画装置および方法、並びにマイクログリッドの運転計画装置で使用される地域エネルギー管理装置およびエネルギー管理装置 | |
Mo et al. | Design method of combined cooling, heating, and power system coupled with cascaded latent heat thermal energy storage based on supply-demand energy-exergy matching | |
CN107269391A (zh) | 基于燃气轮机的医院综合供能系统及其应用 | |
CN109162810A (zh) | 一种用于分布式能源站的多能流测点系统 | |
CN207111237U (zh) | 一种应用于大型公建的分布式能源系统 | |
CN102373998B (zh) | 一种化石燃料与地热相结合的生活热水-电联供系统及方法 | |
Wu et al. | Multi-parameter cooperative optimization and solution method for regional integrated energy system | |
CN205156426U (zh) | 一种集成热化学过程的热电冷多联产系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |