CN109063292B - 一种考虑散热系数区间的热网区间潮流计算方法 - Google Patents

一种考虑散热系数区间的热网区间潮流计算方法 Download PDF

Info

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
Application number
CN201810804007.7A
Other languages
English (en)
Other versions
CN109063292A (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.)
Tsinghua University
Original Assignee
Tsinghua 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 Tsinghua University filed Critical Tsinghua University
Priority to CN201810804007.7A priority Critical patent/CN109063292B/zh
Priority to PCT/CN2018/112695 priority patent/WO2020015241A1/zh
Publication of CN109063292A publication Critical patent/CN109063292A/zh
Priority to US17/146,767 priority patent/US11415952B2/en
Application granted granted Critical
Publication of CN109063292B publication Critical patent/CN109063292B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B17/00Systems involving the use of models or simulators of said systems
    • G05B17/02Systems involving the use of models or simulators of said systems electric
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F24HEATING; RANGES; VENTILATING
    • F24DDOMESTIC- OR SPACE-HEATING SYSTEMS, e.g. CENTRAL HEATING SYSTEMS; DOMESTIC HOT-WATER SUPPLY SYSTEMS; ELEMENTS OR COMPONENTS THEREFOR
    • F24D19/00Details
    • F24D19/10Arrangement or mounting of control or safety devices
    • F24D19/1006Arrangement or mounting of control or safety devices for water heating systems
    • F24D19/1066Arrangement or mounting of control or safety devices for water heating systems for the combination of central heating and domestic hot water
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal 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)热网耦合约束:
Figure BDA0001737751280000021
其中,i为热网中的节点编号,N为节点总数,定义节点温度为:所有流量方向为流出节点的管道的温度,j为热网中的管道编号,M为管道总数,Tend,j表示管道j的末端温度,Tstart,j表示管道的首端温度,Ta为环境温度,λj为管道j的散热系数,Lj为管道j的长度,Cp为管道流体的比热容系数,
Figure BDA0001737751280000022
表示节点相连的管道中,所有流量方向为流出节点的管道温度乘上对应流量的和,
Figure BDA0001737751280000023
表示节点相连的管道中,所有流量方向为流入节点的管道温度乘上对应流量的和,λi,l、λi,u分别表示散热系数的上下解;
(4-2)热望管道散热系数上下界约束:
Figure BDA0001737751280000024
λ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)热网耦合约束:
Figure BDA0001737751280000031
其中,i为热网中的节点编号,N为节点总数,定义节点温度为:所有流量方向为流出节点的管道的温度,j为热网中的管道编号,M为管道总数,Tend,j表示管道j的末端温度,Tstart,j表示管道的首端温度,Ta为环境温度,λj为管道j的散热系数,Lj为管道j的长度,Cp为管道流体的比热容系数,
Figure BDA0001737751280000041
表示节点相连的管道中,所有流量方向为流出节点的管道温度乘上对应流量的和,由于流量方向为流出节点,故各管道对应温度相等。
Figure BDA0001737751280000042
表示节点相连的管道中,所有流量方向为流入节点的管道温度乘上对应流量的和,λi,l、λi,u分别表示散热系数的上下解;
(4-2)热望管道散热系数上下界约束:
Figure BDA0001737751280000043
λ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)热网耦合约束:
Figure FDA0001737751270000011
Figure FDA0001737751270000012
其中,i为热网中的节点编号,N为节点总数,定义节点温度为:所有流量方向为流出节点的管道的温度,j为热网中的管道编号,M为管道总数,Tend,j表示管道j的末端温度,Tstart,j表示管道的首端温度,Ta为环境温度,λj为管道j的散热系数,Lj为管道j的长度,Cp为管道流体的比热容系数,
Figure FDA0001737751270000013
表示节点相连的管道中,所有流量方向为流出节点的管道温度乘上对应流量的和,
Figure FDA0001737751270000014
表示节点相连的管道中,所有流量方向为流入节点的管道温度乘上对应流量的和,λ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)利用信赖域反射算法、有效集算法、内点法或序列二次规划算法,计算得到所有供水网络节点温度和回水网络节点温度的上界和下界,实现考虑散热系数区间的热网区间潮流计算。
CN201810804007.7A 2018-07-20 2018-07-20 一种考虑散热系数区间的热网区间潮流计算方法 Active CN109063292B (zh)

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)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115036923B (zh) * 2022-08-12 2022-10-25 南京邮电大学 一种考虑多能流时序的电-气综合能源系统状态估计方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
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 清华大学 一种基于双侧等效模型的热网稳态运行状态估计方法

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