CN109063292B - 一种考虑散热系数区间的热网区间潮流计算方法 - Google Patents
一种考虑散热系数区间的热网区间潮流计算方法 Download PDFInfo
- Publication number
- CN109063292B CN109063292B CN201810804007.7A CN201810804007A CN109063292B CN 109063292 B CN109063292 B CN 109063292B CN 201810804007 A CN201810804007 A CN 201810804007A CN 109063292 B CN109063292 B CN 109063292B
- Authority
- CN
- China
- Prior art keywords
- supply network
- interval
- heat dissipation
- pipeline
- 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.)
- Active
Links
- 230000017525 heat dissipation Effects 0.000 title claims abstract description 39
- 238000004364 calculation method Methods 0.000 title claims abstract description 19
- 230000008878 coupling Effects 0.000 claims abstract description 11
- 238000010168 coupling process Methods 0.000 claims abstract description 11
- 238000005859 coupling reaction Methods 0.000 claims abstract description 11
- 238000010438 heat treatment Methods 0.000 claims abstract description 8
- 238000000034 method Methods 0.000 claims abstract description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 18
- 238000005457 optimization Methods 0.000 claims description 6
- 238000013461 design Methods 0.000 claims description 3
- 239000012530 fluid Substances 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims description 3
- 230000005855 radiation Effects 0.000 abstract description 6
- 230000032683 aging Effects 0.000 abstract description 2
- 238000009413 insulation Methods 0.000 abstract description 2
- 235000006679 Mentha X verticillata Nutrition 0.000 description 2
- 235000002899 Mentha suaveolens Nutrition 0.000 description 2
- 235000001636 Mentha x rotundifolia Nutrition 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000005485 electric heating Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012954 risk control Methods 0.000 description 1
- 230000033772 system development Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B17/00—Systems involving the use of models or simulators of said systems
- G05B17/02—Systems involving the use of models or simulators of said systems electric
-
- 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
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F24—HEATING; RANGES; VENTILATING
- F24D—DOMESTIC- OR SPACE-HEATING SYSTEMS, e.g. CENTRAL HEATING SYSTEMS; DOMESTIC HOT-WATER SUPPLY SYSTEMS; ELEMENTS OR COMPONENTS THEREFOR
- F24D19/00—Details
- F24D19/10—Arrangement or mounting of control or safety devices
- F24D19/1006—Arrangement or mounting of control or safety devices for water heating systems
- F24D19/1066—Arrangement or mounting of control or safety devices for water heating systems for the combination of central heating and domestic hot water
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Thermal Sciences (AREA)
- Mechanical Engineering (AREA)
- Combustion & Propulsion (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Automation & Control Theory (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Control Of Temperature (AREA)
Abstract
本发明涉及一种考虑散热系数区间的热网区间潮流计算方法,属于电‑热耦合系统运行和控制技术领域。本方法考虑了热网管道参数中,厂家的散热系数和真值存在偏差,同时真实的散热系数会随着管道和保温层的老化而发生改变的情况,因此将散热系数描述为一个区间,求解集中供热系统在参数区间下的潮流解。这样的潮流解可以为系统运行控制提供更为全面的信息,方便运行者在控制和调节中考虑到散热系数与真实值偏差带来的安全上的风险。
Description
技术领域
本发明涉及一种考虑散热系数区间的热网区间潮流计算方法,属于电-热耦合系统运行和控制技术领域。
背景技术
为提高能源总体使用效率,扩大可再生能源消纳能力,对多类能源互联集成和互补融合的需求增长迅速,综合能源系统成为能源系统发展的重要趋势之一。其中,电热耦合系统是综合能源系统的一个重要组成部分。传统的供热系统对于控制的精确要求不高,而电热耦合系统则对原有的控制提出了新的要求。在传统供热系统中,往往采用调节温度的方式来进行调节,通过采集能源站、负荷处包括流量、压力、节点温度在内的信息,结合厂家所给的管道长度、内径、关闭粗糙度和散热系数等参数,求解热网潮流。然而,厂家所给的参数往往会与真实值存在一定的偏差,同时,随着运行时间和运行工况的变化,管道的一些参数也会随着发生改变,传统供热系统由于控制要求不高,往往忽略了这一问题。电热耦合系统由于需要进行更多的调控和风险控制,对于参数的不确定性提出了新的要求。其中,散热系数是影响热网温度分布的一个重要参数,因此,需要研究其不确定性对于热网潮流的影响。
发明内容
本发明的目的是提出一种考虑散热系数区间的热网区间潮流计算方法,以克服现有利用厂家所给参数计算热网潮流存在误差的不足,满足目前快速发展的综合能源系统多能联供的质量要求。
本发明提出的考虑散热系数区间的热网区间潮流计算方法,包括以下步骤:
本发明提出的考虑散热系数区间的热网区间潮流计算方法,包括以下步骤:
(1)从热网的数据采集与监视控制系统中,获得热网系统的状态估计结果,包括:管道的流量、能源站的出水温度、负荷的出水温度;
(2)从厂家所给的设计手册中,获取电-热耦合多能流系统中的管道散热系数λi;
(3)设定所有的管道散热系数位于给定的区间[λi,l,λi,u],区间[λi,l,λi,u]由厂家提供的管道散热系数乘上集中供热系统考虑的散热系数波动比σl和σu得到:
λi,l=λi×σl
λi,u=λi×σu
其中波动比σl和σu为设定值,满足:
0<σl<1
σu>1
(4)建立热网的区间潮流计算模型,包括:
(4-1)热网耦合约束:
其中,i为热网中的节点编号,N为节点总数,定义节点温度为:所有流量方向为流出节点的管道的温度,j为热网中的管道编号,M为管道总数,Tend,j表示管道j的末端温度,Tstart,j表示管道的首端温度,Ta为环境温度,λj为管道j的散热系数,Lj为管道j的长度,Cp为管道流体的比热容系数,表示节点相连的管道中,所有流量方向为流出节点的管道温度乘上对应流量的和,表示节点相连的管道中,所有流量方向为流入节点的管道温度乘上对应流量的和,λi,l、λi,u分别表示散热系数的上下解;
(4-2)热望管道散热系数上下界约束:
λi,l,λi,u由的取值由步骤(3)决定。
(4-3)建立热网区间潮流计算模型的优化目标:
maxTsupply(return),i,i=1,2,…,N
minTsupply(return),i,i=1,2,…,N
其中,Tsupply,i、Treturn,i分别表示热网管道i的供水节点温度和回水节点温度。该目标函数为多目标的多次优化问题,每次只优化一个变量的最大值或最小值;
(5)利用信赖域反射算法、有效集算法、内点法或序列二次规划算法,计算得到所有供水网络节点温度和回水网络节点温度的上界和下界,实现考虑散热系数区间的热网区间潮流计算。
本发明提出的考虑是散热系数区间的热网区间潮流计算方法,其优点是:
本发明方法考虑了热网管道参数中,厂家的散热系数和真值存在偏差,同时真实的散热系数会随着管道和保温层的老化而发生改变的情况,因此将散热系数描述为一个区间,求解集中供热系统在参数区间下的潮流解。这样的潮流解可以为系统运行控制提供更为全面的信息,方便运行者在控制和调节中考虑到散热系数与真实值偏差带来的安全上的风险。本发明方法可以克服现有利用厂家所给参数计算热网潮流存在误差的不足,得到更为保守、可靠的热网潮流解,有利于热网的安全运行与控制。
具体实施方式
本发明提出的考虑散热系数区间的热网区间潮流计算方法,包括以下步骤:
(1)从热网的数据采集与监视控制系统中,获得热网系统的状态估计结果,包括:管道的流量、能源站的出水温度、负荷的出水温度;
(2)从厂家所给的设计手册中,获取电-热耦合多能流系统中的管道散热系数λi;
(3)设定所有的管道散热系数位于给定的区间[λi,l,λi,u],区间[λi,l,λi,u]由厂家提供的管道散热系数乘上集中供热系统考虑的散热系数波动比σl和σu得到:
λi,l=λi×σl
λi,u=λi×σu
其中波动比σl和σu为设定值,满足:
0<σl<1
σu>1
(4)建立热网的区间潮流计算模型,包括:
(4-1)热网耦合约束:
其中,i为热网中的节点编号,N为节点总数,定义节点温度为:所有流量方向为流出节点的管道的温度,j为热网中的管道编号,M为管道总数,Tend,j表示管道j的末端温度,Tstart,j表示管道的首端温度,Ta为环境温度,λj为管道j的散热系数,Lj为管道j的长度,Cp为管道流体的比热容系数,表示节点相连的管道中,所有流量方向为流出节点的管道温度乘上对应流量的和,由于流量方向为流出节点,故各管道对应温度相等。表示节点相连的管道中,所有流量方向为流入节点的管道温度乘上对应流量的和,λi,l、λi,u分别表示散热系数的上下解;
(4-2)热望管道散热系数上下界约束:
λi,l,λi,u由的取值由步骤(3)决定。
(4-3)建立热网区间潮流计算模型的优化目标:
maxTsupply(return),i,i=1,2,…,N
minTsupply(return),i,i=1,2,…,N
其中,Tsupply,i、Treturn,i分别表示热网管道i的供水节点温度和回水节点温度。该目标函数为多目标的多次优化问题,每次只优化一个变量的最大值或最小值;
(5)利用信赖域反射算法(Trust region reflective)、有效集算法(Activeset)、内点法(Interior point)或序列二次规划算法(SQP),计算得到所有供水网络节点温度和回水网络节点温度的上界和下界,实现考虑散热系数区间的热网区间潮流计算。
Claims (1)
1.一种考虑散热系数区间的热网区间潮流计算方法,其特征在于该方法包括以下步骤:
本发明提出的考虑散热系数区间的热网区间潮流计算方法,包括以下步骤:
(1)从热网的数据采集与监视控制系统中,获得热网系统的状态估计结果,包括:管道的流量、能源站的出水温度、负荷的出水温度;
(2)从厂家所给的设计手册中,获取电-热耦合多能流系统中的管道散热系数λi;
(3)设定所有的管道散热系数位于给定的区间[λi,l,λi,u],区间[λi,l,λi,u]由厂家提供的管道散热系数乘上集中供热系统考虑的散热系数波动比σl和σu得到:
λi,l=λi×σl
λi,u=λi×σu
其中波动比σl和σu为设定值,满足:
0<σl<1
σu>1
(4)建立热网的区间潮流计算模型,包括:
(4-1)热网耦合约束:
其中,i为热网中的节点编号,N为节点总数,定义节点温度为:所有流量方向为流出节点的管道的温度,j为热网中的管道编号,M为管道总数,Tend,j表示管道j的末端温度,Tstart,j表示管道的首端温度,Ta为环境温度,λj为管道j的散热系数,Lj为管道j的长度,Cp为管道流体的比热容系数,表示节点相连的管道中,所有流量方向为流出节点的管道温度乘上对应流量的和,表示节点相连的管道中,所有流量方向为流入节点的管道温度乘上对应流量的和,λi,l、λi,u分别表示散热系数的上下解;
(4-2)热望管道散热系数上下界约束:
λi,l ≤λi≤λi,u
λi,l,λi,u由的取值由步骤(3)决定。
(4-3)建立热网区间潮流计算模型的优化目标:
max Tsupply(return),i,i=1,2,…,N
min Tsupply(return),i,i=1,2,…,N
其中,Tsupply,i、Treturn,i分别表示热网管道i的供水节点温度和回水节点温度。该目标函数为多目标的多次优化问题,每次只优化一个变量的最大值或最小值;
(5)利用信赖域反射算法、有效集算法、内点法或序列二次规划算法,计算得到所有供水网络节点温度和回水网络节点温度的上界和下界,实现考虑散热系数区间的热网区间潮流计算。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810804007.7A CN109063292B (zh) | 2018-07-20 | 2018-07-20 | 一种考虑散热系数区间的热网区间潮流计算方法 |
PCT/CN2018/112695 WO2020015241A1 (zh) | 2018-07-20 | 2018-10-30 | 一种考虑散热系数区间的热网区间潮流计算方法 |
US17/146,767 US11415952B2 (en) | 2018-07-20 | 2021-01-12 | Method, apparatus, and storage medium for controlling heating system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810804007.7A CN109063292B (zh) | 2018-07-20 | 2018-07-20 | 一种考虑散热系数区间的热网区间潮流计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109063292A CN109063292A (zh) | 2018-12-21 |
CN109063292B true CN109063292B (zh) | 2022-12-23 |
Family
ID=64817673
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810804007.7A Active CN109063292B (zh) | 2018-07-20 | 2018-07-20 | 一种考虑散热系数区间的热网区间潮流计算方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US11415952B2 (zh) |
CN (1) | CN109063292B (zh) |
WO (1) | WO2020015241A1 (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115036923B (zh) * | 2022-08-12 | 2022-10-25 | 南京邮电大学 | 一种考虑多能流时序的电-气综合能源系统状态估计方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2965019A4 (en) * | 2013-03-04 | 2017-01-11 | Greensleeves LLC | Energy management systems and methods of use |
CN106022624B (zh) * | 2016-05-27 | 2019-07-26 | 清华大学 | 一种电-热耦合多能流网络状态估计方法 |
CN106067677B (zh) * | 2016-05-27 | 2018-07-06 | 清华大学 | 一种电-热耦合多能流系统静态安全分析方法 |
CN106056478B (zh) * | 2016-06-12 | 2019-05-21 | 清华大学 | 一种电-热耦合系统中热网的区间潮流计算方法 |
CN106447529A (zh) * | 2016-08-30 | 2017-02-22 | 上海交通大学 | 考虑热水管网的分布式能源系统建模和运行优化方法 |
CN106532782B (zh) * | 2016-12-12 | 2021-09-21 | 中国电力科学研究院有限公司 | 一种提升风电利用率的热电联产机组运行优化方法及装置 |
CN107563674B (zh) * | 2017-10-09 | 2020-07-10 | 清华大学 | 一种考虑管道动态特性的电-热耦合系统状态估计方法 |
CN107817681B (zh) * | 2017-10-16 | 2020-04-24 | 清华大学 | 一种基于双侧等效模型的热网稳态运行状态估计方法 |
-
2018
- 2018-07-20 CN CN201810804007.7A patent/CN109063292B/zh active Active
- 2018-10-30 WO PCT/CN2018/112695 patent/WO2020015241A1/zh active Application Filing
-
2021
- 2021-01-12 US US17/146,767 patent/US11415952B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN109063292A (zh) | 2018-12-21 |
US11415952B2 (en) | 2022-08-16 |
US20210132558A1 (en) | 2021-05-06 |
WO2020015241A1 (zh) | 2020-01-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112128841B (zh) | 一种基于负荷预测和室温反馈修正的全网平衡调节方法 | |
CN106895478B (zh) | 一种供热管网调节方法和系统 | |
CN107817681B (zh) | 一种基于双侧等效模型的热网稳态运行状态估计方法 | |
CN111222213B (zh) | 一种热力网络动态仿真方法及装置 | |
Zhou et al. | Dynamic modeling of thermal conditions for hot-water district-heating networks | |
CN108920866B (zh) | 基于滚动时域估计理论的热网动态调节运行参数估计方法 | |
CN104281057A (zh) | 一种应用于变压器冷却系统的复合pid模糊控制方法 | |
CN111928335B (zh) | 一种基于智能阀二次网水力平衡方法 | |
CN111829059B (zh) | 一种供热系统动态建模方法、模型及调控系统 | |
CN114491817B (zh) | 一种船舶核动力管网仿真中的母管模拟方法及系统 | |
CN103925641A (zh) | 一种供热管网水力平衡自动调节方法 | |
JP6033674B2 (ja) | 熱供給制御装置、熱供給システム及び熱供給制御方法 | |
CN109063292B (zh) | 一种考虑散热系数区间的热网区间潮流计算方法 | |
CN111898816A (zh) | 一种综合能源系统动态状态估计方法 | |
Wang et al. | Investigation on the use of pumps in HVAC systems for providing ancillary services in smart grids | |
CN113111555B (zh) | 一种基于叠加解耦法的质调节热力系统能流快速计算方法 | |
CN101726038A (zh) | 一种集中供热控制系统 | |
Mu et al. | A novel air flowrate control method based on terminal damper opening prediction in multi-zone VAV system | |
CN112100778A (zh) | 一种综合能源系统的调度方法和系统 | |
CN106839083B (zh) | 一种智能用热管理系统及其管理方法 | |
CN209484698U (zh) | 一种基于热力网和水力发电机的综合供热系统 | |
IT202100004994A1 (it) | Metodo e sistema di controllo di un impianto fluidico | |
Dabis et al. | Control of Domestic Hot Water production is instantaneous heating system with a speed controlled pump | |
CN115906587B (zh) | 一种综合能源系统供热管网动态耦合模型建立方法 | |
Vasek et al. | Pulled plug-flow model for 4th generation district heating |
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 |