CN114781102A - 一种城市综合能源***的能量流建模方法及应用 - Google Patents
一种城市综合能源***的能量流建模方法及应用 Download PDFInfo
- Publication number
- CN114781102A CN114781102A CN202210249213.2A CN202210249213A CN114781102A CN 114781102 A CN114781102 A CN 114781102A CN 202210249213 A CN202210249213 A CN 202210249213A CN 114781102 A CN114781102 A CN 114781102A
- Authority
- CN
- China
- Prior art keywords
- heat
- model
- flow
- energy
- inl
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- 230000008569 process Effects 0.000 claims abstract description 16
- 238000009825 accumulation Methods 0.000 claims abstract description 14
- 238000004088 simulation Methods 0.000 claims abstract description 9
- 238000005485 electric heating Methods 0.000 claims abstract description 8
- 238000012360 testing method Methods 0.000 claims abstract description 4
- 239000012530 fluid Substances 0.000 claims description 50
- 238000012546 transfer Methods 0.000 claims description 26
- 230000008859 change Effects 0.000 claims description 7
- 238000010248 power generation Methods 0.000 claims description 7
- 230000014509 gene expression Effects 0.000 claims description 5
- 238000005338 heat storage Methods 0.000 claims description 5
- 238000010438 heat treatment Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000013519 translation Methods 0.000 claims description 4
- 238000004134 energy conservation Methods 0.000 claims description 3
- 230000017525 heat dissipation Effects 0.000 claims description 2
- 238000011002 quantification Methods 0.000 claims description 2
- 230000001360 synchronised effect Effects 0.000 claims description 2
- 230000005540 biological transmission Effects 0.000 abstract description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000010276 construction Methods 0.000 abstract description 2
- 230000010354 integration Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 230000004888 barrier function Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 230000005611 electricity Effects 0.000 description 3
- 230000000295 complement effect Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- 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
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Power Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种城市综合能源***的能量流建模方法及应用,包括:确定城市综合能源***组成设备和网络拓扑结构;建立城市综合能源***中电力***能量流模型,包括设备模型和网络模型;基于火积耗散理论和热电比拟理论,建立城市综合能源***中热力***能量流模型,包括设备模型和网络模型;基于连接电力、热力***子模型得到***整体模型,实现仿真测试和特性分析。本发明针对热力***中的基础设备:供热管道、换热站和建筑热用户建立了一般形式的能量流模型,准确描述热量在各个部分传输和损失的动态过程,有利于较大规模地模型集成和求解,实现了对电热联合***的统一描述和特性分析,为后续调度和控制问题构建提供模型基础。
Description
技术领域
本发明涉及综合能源***技术领域,尤其涉及一种城市综合能源***能量流建模方法及应用。
背景技术
城市综合能源***集成了风电、光伏、火电机组、热电联产机组、电池、储热罐、电热泵、电锅炉等多种形式的供能、能量转换和储能装置,各种能源形式间存在着复杂耦合关系和特性差异。例如,电能比较容易传输但较难存储,而热能比较易于存储但较难传输,因此电力***和热力***的物理特性具有良好的互补性。综合能源***可以充分利用各类形式能源之间的互补特性实现电、气、热等不同能源之间的互联,推动打破传统单一能源发展的技术壁垒、市场壁垒和体制壁垒,实现多能源的互补互济和协调优化,从而有效提升***整体能源利用效率。
目前,综合能源***的建模分析研究一般是针对电力和热力***分别开展的。电力***一般是基于能量流方法研究电能自身的输运特性,通过欧姆定律和基尔霍夫电压电流定律分别描述元件传输和***拓扑特性构造电力潮流约束。而热力***基于工质流方法同时研究伴随工质流动的热能迁移与工质之间的热能传递,在两者的耦合处以互为边界条件的形式体现两者之间的作用和影响。在对不同能源***模型进行集成时,整体求解的收敛性难以保证,且没有基于电力、热力***的物理机制,从而无法对综合能源***开展精确分析与控制设计。
因此,对城市综合能源***提出一种统一的兼具准确度与简洁性的建模和分析方法,对于***的整体特性分析和控制***设计至关重要。
传统的电热综合能源***中电力、热力子***模型是基于不同方法进行构建的。电力***中只有能量传递,不考虑物质传递,因此可以通过能量流的方法基于电学基本定律,考虑元件特性约束和网络拓扑约束,得到电压和电流之间的网络方程。而热力***中不仅有能量传递,还有物质传递,通过需要考虑机械、流动等多物理过程,其数学模型包含分布参数的偏微分方程和集总参数的非线性方程。在对电力、热力子***模型进行简单集成时,不仅模型的收敛性无法保证,同时模型整体复杂度较高,带来了巨大的计算负担。因此本发明通过提出一种城市综合能源***的能量流建模方法,提供了电热综合能源***统一的建模和分析方法,能够综合考虑模型复杂度和准确性,实现对***整体的精准建模和特性分析。
发明内容
为了解决现有技术中的不足,本发明公开一种城市综合能源***能量流建模方法及应用,其技术方案如下:
一种城市综合能源***的能量流建模方法及应用,其特征为:包括如下步骤:
步骤1:确定城市综合能源***组成设备和网络拓扑结构;
步骤2:建立城市综合能源***中电力***能量流模型,包括设备模型和电力***网络功率平衡模型;
步骤3:基于火积耗散理论和热电比拟思想,建立城市综合能源***中热力***能量流模型,包括设备模型和网络模型;
步骤4:基于连接电力、热力***子模型得到***整体模型,实现仿真测试和特性分析。
本发明还公开一种电热综合能源***,其特征为:采用上述城市综合能源***的能量流建模方法实现所述***整体建模和特性分析。
有益效果
与现有技术相比,本发明可以获得包括以下技术效果:
1、综合能源***中热力***规模庞大,结构复杂,设备繁多。本发明基于能量流方法,针对热力***中的基础设备:供热管道、换热站和建筑热用户,分别建立了一般形式的能量流模型,可以准确描述热量在各个部分传输和损失的动态过程,同时有利于较大规模地模型集成和求解。
2、本发明针对电热综合能源***建立了统一的能量流模型,实现了对电热联合***的统一描述和特性分析,可为后续调度和控制问题构建提供模型基础。
附图说明
图1为本发明城市综合能源***示意图;
图2为本发明等效建筑热用户示意图;
图3为本发明环境温度变化曲线示意图;
图4为本发明建筑物内部温度变化曲线示意图。
具体实施方式
实施例1
一种城市综合能源***的能量流建模方法及应用,其特征为:包括如下步骤:
步骤1:确定城市综合能源***组成设备和网络拓扑结构;
所述城市综合能源***的设备主体包括:火电机组、热电联产机组、风电机组、换热站、供热管道,建筑热用户,参见图1所示。综合能源***中风电机组、火电机组、热电联产机组共同为用户提供电负荷,热电联产机组与供热***连通为用户提供热负荷;所述供热***包括一次换热站、一次热网、二次换热站以及二次热网,首先热电联产机组供热抽气通过一次换热站将热量传递给一次热网,一次热网通过二次换热站将热量传递给二次热网,二次热网再通过建筑内散热器将热量传递给建筑物内部,从而满足用户的热负荷需求。
步骤2:建立城市综合能源***中电力***能量流模型,包括设备模型和电力***网络功率平衡模型:
该步骤涉及电力***基础设备包括火电机组、热电联产机组和风电机组,
其中:火电机组模型为t时刻第i台火电机组输出功率模型如下:
热电联产机组模型为t时刻第i台热电联产机组模型如下:
风力发电机组模型为t时刻第i台风电机组模型如下:
步骤2中涉及的电力***网络功率平衡模型为:
步骤3:基于火积耗散理论和热电比拟理论,建立城市综合能源***中热力***能量流模型,包括设备模型和网络模型。
该步骤中涉及的电力***基础设备包括供热管道、换热站和建筑热用户;所述供热管道具体包括一次热网和二次热网,由于一次热网的长度足够长,因此在模型建立中需要考虑对应的热量迁移和漏热过程。
建立一次热网能量流模型的具体方法为:由于流体在管道中的输运需求一定时间,且随着流体在管道内部输运过程中不断与环境进行热交换,因此沿管长方向流体温度也是动态变化的。因此,根据流体输运特性可以建立管道进出口流体温度之间的关系模型。忽略管壁储热和流体之间导热,对一段长度为ΔL的流体微元建立其热平衡方程为:
K(T-Tamb)dΔLdt=-ρcpSΔLdT
其中:K、T、Tamb、d、ρ、cp、S分别表示换热系数、流体温度、环境温度、管道周长、流体密度、流体比热容、管道截面积;ΔL表示流体微元长度,dt表示时间变化值,dT表示温度变化值;该热平衡方程对应的边界条件为:
T(tinl,tinl)=Tinl(tinl)
其中:Tinl表示管道入口流体温度,tinl表示该流体微元进入管道的时刻,从而可以得到在tinl时刻进入管道的流体随时间变化的温度模型为:
其中:τflow=t-tinl,b=τflow/τloss,τloss=ρcpS/K,θ表示时间积分变量,τflow和τloss分别表示管道流动和漏热过程的特征时间,b是对应比值,反映了热量迁移特性和热量传递特性的主导地位;同时,流体微元通过管道的总时长τflow,t可以用下式来确定:
其中,L表示管道总长,v表示流动速度。τflow,t表示管道流动的延迟时间;流体微元从管道入口到出口的总内能变化表达式为:
ΔU(tinl,τflσw,t)=ρcpSΔl[T(tinl,tinl)-T(tinl,tinl+τflow,t)]
其中,ΔU(tinl,τflow,t)是流体微元的总内能变化;从而可以得到单位时间内在tinl时刻流入管道的流体微元流经管道的总漏热量:
其中:
εh,delay,loss(tinl,,τflow,t)是由于漏热引起的平均传热温度变化,Rh,loss是反映整个管道漏热特性的换热热阻,bp=τflow,t/τloss是管道整体热量迁移特性与热量传递特性之间相对重要程度的度量;
同时,流体出口温度可以表示为
Toutl(tinl+τflow,t)=T(tinl,tinl+τflow,t)=Tinl(tinl)+εloss(tinl,τflow,t)
其中:
在此基础上,引入一个描述时间平移的附加热动势εh,delay,flow将上述异步能量流模型的出口温度时间平移至当前时刻温度,计算式为:
εh,delay,flow(tinl,τflow,t)=T(tinl,tinl+τflow,t)-T(tinl-τflow,t,tinl)
平移时间后,瞬时漏热量Φloss,instan可以得到为
考虑到管道内流体的储热能力,得到管道流体的瞬时热平衡方程为:
其中:
Up=GPSTp=ρcpSLTp
Ch,Ps和Tp分别是管内流体的总热容和平均温度。从而可以得到一次管网的能量流模型,其中εh,trans是Φloss,trans由于异步模型向同步模型时间平移所需要的附加热动势和热流源,表达式分别为
εh,trans=Tinl-Tp
Φloss,trans=φloss-φloss,instan
所建立的一次管网能量流模型完整反映了管道输运过程中热量的传递、迁移和存储特性。
建立换热站能量流模型和热用户能量流模型的具体方法为::涉及建筑热用户能量流模型是一种等效建筑模型,包括两部分,一是建筑模型本身,二是一次热网和等效建筑之间的换热模型。等效建筑热用户示意图如图2所示,具体换热过程为:热量首先通过二次换热站由一次热网传递给二次热网,再经过建筑热用户的散热器由二次热网传递给房屋。由于二次热网一般距离热用户较短,所以不考虑二次热网供热管道的蓄热能力,且二次热网与环境的散热过程,则有:
Hhl1=Hhl2
其中:Hhl1、Hhl2分别表示二次换热站换热量和建筑散热器换热量;
对于二次换热站,假设其采用逆流换热器,基于火积耗散理论可以得到逆流换热器的等效热阻为:
其中:ap,hl=khlAhl/Gp,hl,ahl=khlAhl/Ghl,Gp,hl=mP,hlcp,Ghl=mhlcp。表示二次换热站逆流换热器火积耗散热阻,mp,hl和mhl分别表示一次热网、二次热网热流体质量流量;cp表示流体工质比热容,khl、Ahl分别表示二次换热站逆流换热器的换热系数和换热面积;Thl,in、Thl,out分别表示一次热网流体进入和离开二次换热站的温度;Th,i、Th,o分别表示二次热网流体离开、进入二次换热站的温度。
基于能量守恒定量,有如下关系:
从而可以得到
其中:Ren,hl是二次换热站火积热阻的另一种表示形式,该换热器的换热过程可以用温差除以基于火积耗散的热阻表示;从热电比拟原理来说,Hhl1可看作电路中的电流,而温差可以看作是电势差。
步骤3涉及的热力***二次热网实际运行中,一般通过调节其质量流量来保持建筑物室内温度在人体舒适度范围内相对恒定;因此在建模过程中假设二次热网中配置有比例-积分控制器,该比例-积分控制器输入量为建筑物室内温度与其设定值的差值,经比例-积分环节计算后产生控制增量叠加到二次热网流量设定值一起作为二次热网流量的实际控制量,从而维持建筑物室内温度的相对恒定。
建筑热用户模型主要建立由于房屋温度与环境温度不同导致的制热、制冷需求模型。同样假设建筑热用户散热器为逆流换热器,则基于火积耗散理论可以得到该换热器的等效热阻为:
其中:ahl2=kbAb/Ghl,ah=kbAb/Gh,Ghl=mhlcp,Gh=mhch。mh表示散热器房屋侧流体质量流量,ch表示散热器房屋侧流体比热容。kb、Ab别表示逆流换热器的换热系数和换热面积。在房屋侧假设室内温度保持均匀一致,因此有ch=∞,Gh=∞,ah=0。则换热器等效热阻可简化为:
其中:Ren,b是建筑物散热器火积热阻的另一种表示形式;
对建筑物本身根据能量守恒定量可以建立模型如下:
其中,Hb表示建筑物与环境热交换能量,Cb表示建筑物比热容,Tb表示建筑物内部温度,Rb表示建筑物散热过程的等效热阻。
步骤3中涉及的热力***网络平衡模型为:
Tb,min≤Tb,i,t≤Tb,max
其中:Tb,i,t表示t时刻第i个等效建筑内部温度;Tb,min、Tb,max分别表示符合人体舒适度的室内温度下限和上限。通过将建筑热用户室内温度控制在给定区间内,表明热力***中的热负荷达到供需平衡。
步骤4:基于连接电力、热力***子模型得到***整体模型,实现仿真测试和特性分析。
基于步骤2-3分别得到的电力、热力***能量流模型,按照供能原理连接,以电力、热力***网络平衡方程为约束条件,联立求解上述能量流模型,给定输入参数,如火电机组功率,风力发电机组功率,热电联产机组功率,供热机组抽气流量和温度,供热管网基本参数,换热站换热面积和传热系数,建筑物基本参数等,实现给定参数下的***整体模型的求解仿真。图3表示给定的外界环境温度变化曲线,图4表示对应的三处建筑物内部温度变化曲线。该图表明当外界环境温度随时间变化时,通过二次热网回路中控制器调节二次热网供热质量流量,可以始终保持建筑物室内温度在人体舒适度范围内。
本发明研究电热综合能源***的统一建模,以能量流为核心,建立和电力***模型体系相一致的热力***模型,包括供热管道、换热站和热用户等热力***基础设备,将热力***模型融入现有电力***模型体系中,形成反映电、热不同能流整体平衡特性和输运规律的统一潮流模型。本发明提供的综合能源***统一能量流方法,兼具模型准确度与简洁性,可为后续综合能源***的规划设计和优化控制提供良好的模型基础。
Claims (9)
1.一种城市综合能源***的能量流建模方法,其特征为:包括如下步骤:
步骤1:确定城市综合能源***组成设备和网络拓扑结构;
步骤2:建立城市综合能源***中电力***能量流模型,包括设备模型和网络模型;
步骤3:基于火积耗散理论和热电比拟理论,建立城市综合能源***中热力***能量流模型,包括设备模型和网络模型;
步骤4:基于连接电力、热力***子模型得到***整体模型,实现仿真测试和特性分析。
2.根据权利要求1所述的城市综合能源***的能量流建模方法,其特征为:所述步骤1进一步包括如下内容:所述城市综合能源***的设备主体包括:火电机组、热电联产机组、风电机组、换热站、供热管道,建筑热用户;综合能源***中风电机组、火电机组、热电联产机组共同为用户提供电负荷,热电联产机组与供热***连通为用户提供热负荷;所述供热***包括一次换热站、一次热网、二次换热站以及二次热网,首先热电联产机组供热抽气通过一次换热站将热量传递给一次热网,一次热网通过二次换热站将热量传递给二次热网,二次热网再通过建筑内散热器将热量传递给建筑物内部,从而满足用户的热负荷需求。
3.根据权利要求1所述的城市综合能源***的能量流建模方法,其特征为:所述步骤2进一步包括如下内容:
该步骤涉及电力***基础设备包括火电机组、热电联产机组和风电机组,
其中:火电机组模型为t时刻第i台火电机组输出功率模型如下:
热电联产机组模型为t时刻第i台热电联产机组模型如下:
风力发电机组模型为t时刻第i台风电机组模型如下:
电力***网络功率平衡模型为:
4.根据权利要求1所述的城市综合能源***的能量流建模方法,其特征为:所述步骤3进一步包括如下内容:该步骤中涉及的电力***基础设备包括供热管道、换热站和建筑热用户;所述供热管道具体包括一次热网和二次热网;建立一次热网能量流模型、换热站能量流模型和热用户能量流模型。
5.根据权利要求4所述的所述的城市综合能源***的能量流建模方法,其特征为:所述一次热网能量流模型的具体方法为:
忽略管壁储热和流体之间导热,对一段长度为ΔL的流体微元建立其热平衡方程为:
K(T-Tamb)dΔLdt=-ρcpSΔLdT
其中:K、T、Tamb、d、ρ、cp、S分别表示换热系数、流体温度、环境温度、管道周长、流体密度、流体比热容、管道截面积;ΔL表示流体微元长度,dt表示时间变化值,dT表示温度变化值;该热平衡方程对应的边界条件为:
T(tinl,tinl)=Tinl(tinl)
其中:Tinl表示管道入口流体温度,tinl表示该流体微元进入管道的时刻,从而可以得到在tinl时刻进入管道的流体随时间变化的温度模型为:
其中:τflow=t-tinl,b=τflow/τloss,τloss=ρcpS/K,θ表示时间积分变量,τflow和τloss分别表示管道流动和漏热过程的特征时间,b是对应比值,反映了热量迁移特性和热量传递特性的主导地位;同时,流体微元通过管道的总时长τflow,t可以用下式来确定:
其中,L表示管道总长,v表示流动速度。τflow,t表示管道流动的延迟时间;流体微元从管道入口到出口的总内能变化表达式为:
ΔU(tinl,τflow,t)=ρcpSΔl[T(tinl,tinl)-T(tinl,tinl+τflow,t)]
其中,ΔU(tinl,τflow,t)是流体微元的总内能变化;从而可以得到单位时间内在tinl时刻流入管道的流体微元流经管道的总漏热量:
其中:
εh,delay,loss(tinl,,τflow,t)是由于漏热引起的平均传热温度变化,Rh,loss是反映整个管道漏热特性的换热热阻,bp=τflow,t/τloss是管道整体热量迁移特性与热量传递特性之间相对重要程度的度量;
同时,流体出口温度可以表示为
Toutl(tinl+τflow,t)=T(tinl,tinl+τflow,t)=Tinl(tinl)+εloss(tinl,τflow,t)
其中:
在此基础上,引入一个描述时间平移的附加热动势εh,delay,flow将上述异步能量流模型的出口温度时间平移至当前时刻温度,计算式为:
εh,delay,flow(tinl,τflow,t)=T(tinl,tinl+τflow,t)-T(tinl-τflow,t,tinl)
平移时间后,瞬时漏热量Φloss,instan可以得到为
考虑到管道内流体的储热能力,得到管道流体的瞬时热平衡方程为:
其中:
Up=CPSTp=ρcpSLTp
Ch,PS和Tp分别是管内流体的总热容和平均温度。从而可以得到一次管网的能量流模型,其中εh,trans是Φloss,trans由于异步模型向同步模型时间平移所需要的附加热动势和热流源,表达式分别为
εh,trans=Tinl-Tp
Φloss,trans=Φloss-Φloss,instan。
6.根据权利要求4所述的所述的城市综合能源***的能量流建模方法,其特征为:
建立换热站能量流模型和热用户能量流模型的具体方法为:涉及建筑热用户模型是一种等效建筑模型,包括两部分,一是建筑模型本身,二是一次热网和等效建筑之间的换热模型;具体换热过程为:热量首先通过二次换热站由一次热网传递给二次热网,再经过建筑热用户的散热器由二次热网传递给房屋。由于二次热网一般距离热用户较短,所以不考虑二次热网供热管道的蓄热能力,且二次热网与环境的散热过程,则有:
Hhl1=Hhl2
其中:Hhl1、Hhl2分别表示二次换热站换热量和建筑散热器换热量;
对于二次换热站,假设其采用逆流换热器,基于火积耗散理论可以得到逆流换热器的等效热阻为:
其中:ap,hl=khlAhl/Gp,hl,ahl=khlAhl/Ghl,Gp,hl=mp,hlcp,Ghl=mhlcp。表示二次换热站逆流换热器火积耗散热阻,mp,hl和mhl分别表示一次热网、二次热网热流体质量流量;cp表示流体工质比热容,khl、Ahl分别表示二次换热站逆流换热器的换热系数和换热面积;Thl,in、Thl,out分别表示一次热网流体进入和离开二次换热站的温度;Th,i、Th,o分别表示二次热网流体离开、进入二次换热站的温度;
基于能量守恒定量,有如下关系:
从而可以得到
其中:Ren,hl是二次换热站火积热阻的另一种表示形式,该换热器的换热过程可以用温差除以基于火积耗散的热阻表示;从热电比拟原理来说,Hhl1可看作电路中的电流,而温差可以看作是电势差。
7.根据权利要求4所述的所述的城市综合能源***的能量流建模方法,其特征为:
建筑热用户模型主要建立由于房屋温度与环境温度不同导致的制热、制冷需求模型;同样假设建筑热用户散热器为逆流换热器,则基于火积耗散理论可以得到该换热器的等效热阻为:
其中:ahl2=kbAb/Ghl,ah=kbAb/Gh,Ghl=mhlcp,Gh=mhch。mh表示散热器房屋侧流体质量流量,ch表示散热器房屋侧流体比热容。kb、Ab别表示逆流换热器的换热系数和换热面积。在房屋侧假设室内温度保持均匀一致,因此有ch=∞,Gh=∞,ah=0;则换热器等效热阻可简化为:
其中:Ren,b是建筑物散热器火积热阻的另一种表示形式;
对建筑物本身根据能量守恒定量可以建立模型如下:
其中,Hb表示建筑物与环境热交换能量,Cb表示建筑物比热容,Tb表示建筑物内部温度,Rb表示建筑物散热过程的等效热阻。
8.根据权利要求1所述的城市综合能源***的能量流建模方法,其特征为:所述步骤4进一步包括如下内容:基于步骤2中得到的电力***能量流模型和步骤3中得到的热力***能量流模型,按照供能原理连接,以电力、热力***网络平衡方程为约束条件,联立求解上述能量流模型,给定输入参数,实现对***整体模型的求解仿真,并开展对应的特性分析。
9.一种电热综合能源***,其特征为:采用上述权利要求1-7任一所述城市综合能源***的能量流建模方法实现所述***整体建模和特性分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210249213.2A CN114781102A (zh) | 2022-03-14 | 2022-03-14 | 一种城市综合能源***的能量流建模方法及应用 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210249213.2A CN114781102A (zh) | 2022-03-14 | 2022-03-14 | 一种城市综合能源***的能量流建模方法及应用 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114781102A true CN114781102A (zh) | 2022-07-22 |
Family
ID=82423608
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210249213.2A Pending CN114781102A (zh) | 2022-03-14 | 2022-03-14 | 一种城市综合能源***的能量流建模方法及应用 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114781102A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115791872A (zh) * | 2023-01-30 | 2023-03-14 | 天津地热开发有限公司 | 一种地埋管换热器传热分析方法及*** |
CN117450555A (zh) * | 2023-12-22 | 2024-01-26 | 瑞纳智能设备股份有限公司 | 一种供热*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019200662A1 (zh) * | 2018-04-20 | 2019-10-24 | 东北大学 | 电-热-气综合能源***的稳定评估与静态控制方法 |
CN111829059A (zh) * | 2020-07-08 | 2020-10-27 | 华北电力大学 | 一种供热***动态建模方法、模型及调控*** |
CN111899125A (zh) * | 2020-07-31 | 2020-11-06 | 华中科技大学 | 一种综合能源***的优化建模运行方法、装置和介质 |
-
2022
- 2022-03-14 CN CN202210249213.2A patent/CN114781102A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019200662A1 (zh) * | 2018-04-20 | 2019-10-24 | 东北大学 | 电-热-气综合能源***的稳定评估与静态控制方法 |
CN111829059A (zh) * | 2020-07-08 | 2020-10-27 | 华北电力大学 | 一种供热***动态建模方法、模型及调控*** |
CN111899125A (zh) * | 2020-07-31 | 2020-11-06 | 华中科技大学 | 一种综合能源***的优化建模运行方法、装置和介质 |
Non-Patent Citations (3)
Title |
---|
张勇等: "基于线性规划的电-热互联***参数优化方法", 研究与探讨, vol. 22, no. 02, 30 March 2020 (2020-03-30), pages 8 - 12 * |
徐荆州等: "基于能源路由器的源网荷储优化管理", 能效与负荷管理, vol. 20, no. 06, 30 November 2018 (2018-11-30), pages 16 - 17 * |
杨帅;陈磊;徐飞;孙亮;陈群;贺克伦;姜拓;戴远航;: "基于能量流的电热综合能源***弃风消纳优化调度模型", 电网技术, no. 02, 12 February 2018 (2018-02-12) * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115791872A (zh) * | 2023-01-30 | 2023-03-14 | 天津地热开发有限公司 | 一种地埋管换热器传热分析方法及*** |
CN115791872B (zh) * | 2023-01-30 | 2023-06-13 | 天津地热开发有限公司 | 一种地埋管换热器传热分析方法及*** |
CN117450555A (zh) * | 2023-12-22 | 2024-01-26 | 瑞纳智能设备股份有限公司 | 一种供热*** |
CN117450555B (zh) * | 2023-12-22 | 2024-03-22 | 瑞纳智能设备股份有限公司 | 一种供热*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2020093296A1 (zh) | 电-热互联综合能源***的区间潮流计算方法 | |
CN108921727B (zh) | 考虑热负荷动态特性的区域综合能源***可靠性评估方法 | |
CN105910169B (zh) | 基于机理模型预测控制的城市供热***热网调节方法及*** | |
CN111815042B (zh) | 计及精细化热网模型的电热综合能源***优化调度方法 | |
CN108599137B (zh) | 考虑区域热网暂态传热特性的多能流***优化运行方法 | |
CN114781102A (zh) | 一种城市综合能源***的能量流建模方法及应用 | |
Gao et al. | Numerical and experimental investigation on thermal stratification characteristics affected by the baffle plate in thermal storage tank | |
WO2019071763A1 (zh) | 一种考虑管道动态特性的电-热耦合***状态估计方法 | |
CN109492325B (zh) | 一种基于扩展能源集线器的多能耦合***潮流分析方法 | |
CN113190942B (zh) | 供热/冷***虚拟储能能力计算方法、装置与电子设备 | |
CN110135013B (zh) | 一种构建供热管网动态热力模型的整合方法 | |
CN111829059B (zh) | 一种供热***动态建模方法、模型及调控*** | |
Li et al. | Gradient descent iterative method for energy flow of integrated energy system considering multiple modes of compressors | |
Lan et al. | Modeling of the enthalpy transfer using electric circuit equivalents: Theory and application to transients of multi-carrier energy systems | |
CN113190999A (zh) | 供热管网流量调节提升风电消纳的电热协调方法及装置 | |
Huang et al. | A multi-rate dynamic energy flow analysis method for integrated electricity-gas-heat system with different time-scale | |
Ge et al. | Dynamic modeling and intelligent heating strategies of district heating system based on the standardized thermal resistance | |
CN113255105B (zh) | 双向耦合网络结构的电气热综合能源***潮流计算方法 | |
CN114048699B (zh) | 综合能源***电热联合电磁暂态仿真模型建立和计算方法 | |
CN115234973B (zh) | 一种热网动态供热策略确定方法 | |
CN116611706A (zh) | 基于多能源主体的动态碳排放因子测算方法 | |
CN111219767A (zh) | 一种考虑热负荷需求响应的电—气—热综合能源***调控方法 | |
CN109241676A (zh) | 综合能源***中热网及建筑物时间分辨率的选择方法 | |
CN115374527A (zh) | 一种区域电-热-冷综合能源动态模拟***构建方法 | |
CN114372349A (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 |