CN104360984A - 一种基于两点边值测量下冷却塔填料热力性能计算方法 - Google Patents
一种基于两点边值测量下冷却塔填料热力性能计算方法 Download PDFInfo
- Publication number
- CN104360984A CN104360984A CN201410724779.1A CN201410724779A CN104360984A CN 104360984 A CN104360984 A CN 104360984A CN 201410724779 A CN201410724779 A CN 201410724779A CN 104360984 A CN104360984 A CN 104360984A
- Authority
- CN
- China
- Prior art keywords
- filler
- cooling tower
- thermal performance
- boundary value
- temperature
- 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
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明公开了一种基于两点边值测量下冷却塔填料热力性能计算方法,该方法包括:将填料中的热质传递用四个常微分方程来描述;采集冷却塔现场运行参数和环境参数,确定填料底部(z=0)的空气温度Ta和湿度比X,填料顶部(z=H)的水温度Tw和质量流率mw边值;假设和迭代调整填料底部(z=0)未知的边界条件,从而使得填料顶部(z=H)满足已知边界条件;基于填料热力性能线性模型建立冷却塔填料热力性能计算模型;代入上述计算得到填料顶部(z=H)的条件,计算冷却塔热力性能。本方法基于实时测量值进行模型建立,有利于实现冷却塔热力性能实时在线评估,计算模型在基于两点边值测量下得到简化,相比已有计算方法,在保证计算精度前提下,较大减少计算时间。
Description
技术领域
本发明涉及冷却塔热力性能计算方法,尤其涉及一种基于两点边值测量下冷却塔填料热力性能计算方法。
背景技术
随着工业飞速发展,工业系统对冷却系统要求日益增加。冷却塔作为主要末端冷却装置,被广泛用于排除电厂、暖通空调系统、钢铁厂产生的大量废热,是循环冷却系统重要组成部分。湿式冷却塔因效率高而成为冷却塔主流,其通过空气与高温循环冷却水传热传质过程,将高温高湿空气排出塔外,带走热量,实现冷却水循环利用,因此冷却塔热力性能的优劣不仅直接影响冷却系统节水节能效果,对工业系统正常稳定运行也具有重大意义。
随着节能环保技术被日益重视,冷却系统的节能节水已成为国内外冷却技术研究的重要课题,由于冷却塔热力性能与系统节水节能直接相关,国务院《节能减排“十二五”规划》中指出:要“开展大型公共建筑采暖、空调、通风、照明等节能改造”。国家针对冷却塔具有极大节能潜力,应提高冷却塔能源利用效率问题,提出下一步冷却塔产品必须做节水节能认证,中国质量认证中心出台了《冷却塔节水认证规则》,并即将出台《开式冷却塔节能认证技术规范》予以保证,冷却塔加装热力性能在线监测仪有助于这项工作开展。此外,由于利益驱使,部分冷却塔生产商以次充好配置冷却系统,影响冷却效果,2012年国内某公司就因冷却塔能效设计引起诉讼,这种情况普遍存在将导致用户无法验收,出现能耗过大或故障,对冷却塔加装监测系统显得更加重要。
目前,经典的冷却塔热力性能计算模型有Merkel模型、e-NTU模型和Poppe模型。Merkel模型基于假设建立,模型简单,计算误差大,且需要采用数值积分或迭代法求解,e-NTU模型引入传热单元数,避免数值积分或迭代计算,但和Merkel模型具有相同等级的计算误差,Poppe模型计算结果准确,但模型复杂,求解过程涉及多重迭代,计算时间较长,经典的冷却塔热力性能计算模型无法同时满足高精度和计算量小的要求。随着计算机技术的发展,人工智能算法被引入冷却塔热力性能评估中,避免冷却塔机理建模过程,但需要大量训练样本,且评估时间较长。
发明内容
为解决上述冷却塔填料热力性能计算模型和方法的问题与缺陷,该方法通过将填料中的热质传递问题用一组求解未知函数导数的四个常微分方程来描述,然后将模型简化为两点边值问题,简化冷却塔热力性能计算模型,减少计算时间。本发明提供了一种基于两点边值测量下冷却塔热力性能计算方法。
本发明是通过以下技术方案实现的:
一种基于两点边值测量下冷却塔填料热力性能计算方法,包括四个步骤:
A通过一组求解未知函数导数的四个常微分方程描述填料中的热质传递;
B采集冷却塔现场运行参数和环境参数,确定填料底部即z=0的空气温度Ta和湿度比X,填料顶部即z=H的水温度Tw和质量流率mw边值;
C假设和迭代调整填料底部z=0未知的边界条件,从而使得填料顶部z=H满足已知边界条件;
D基于填料热力性能线性模型建立冷却塔填料热力性能计算模型;
E代入上述计算得到填料顶部z=H的条件,计算冷却塔热力性能。
本发明提供的技术方案的有益效果是:
基于实时测量值进行模型建立,有利于实现冷却塔热力性能实时在线评估,计算模型在基于两点边值测量下得到简化,相比已有评估方法,在保证计算精度前提下,较大减少计算时间。本方法实时性好、适应性强特点,可广泛应用于机械式逆流冷却塔现场热力性能实时评估中。
附图说明
图1是基于两点边值测量下冷却塔填料热力性能计算方法流程图;
图2是逆流方塔的热力性能计算系统安装示意图;
图3是基于两点边值测量下冷却塔填料热力性能计算模型图;
图4是冷却塔填料热力性能线性模型图;
图5是基于两点边值测量下冷却塔填料热力性能计算模型焓湿图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图,以CEF-370不锈钢逆流方塔的热力性能计算为例,对本发明实施方式作进一步地详细描述:
图1是根据本发明实施例的冷却塔填料热力性能计算方法流程图,所述方法包括以下步骤:
步骤10将填料中的热质传递用一组形式如下的求解未知函数导数的四个常微分方程来描述:
已知的四个条件可归纳为:
Ta(z=0)=Tai
X(z=0)=Xi
Tw(z=H)=Twi
mw(z=H)=mwi
步骤20采集冷却塔现场运行参数和环境参数,确定填料底部(z=0)的空气温度Ta和湿度比X,填料顶部(z=H)的水温度Tw和质量流率mw边值;
启动图2所示CEF-370不锈钢逆流方塔的热力性能计算系统,对冷却塔现场进水量mw、进风量ma、进水温度Twi、出水温度Two、出风温度Tao等运行参数和大气压力P0、湿球温度Twb、进风干球温度Tai等环境参数进行采集。
步骤30假设和迭代调整填料底部(z=0)未知的边界条件,从而使得填料顶部(z=H)满足已知边界条件;
假设和迭代调整填料底部(z=0)未知的边界条件 Twb为入口空气的湿球温度,为在出风口温度近似值Tao下计算出的饱和湿度比,从而使得填料顶部(z=H)满足已知边界条件。已知边界条件是:将填料底部(z=0)的未知两个条件(或者说缺少的条件)进行假设,然后根据这个假设去递推填料顶部(z=H)的四个条件,而且递推得到的这些条件与已知的两个条件满足。
步骤40基于填料热力性能线性模型(如图4)建立冷却塔填料热力性能计算模型,可以看出MeTPBVP=MeTPBVP(z=H);
步骤50根据步骤30计算的结果和实时测量值,确定填料顶部(z=H)气温、水温、大气压力和湿球温度Tao、Twi、P0、Twb;
将冷却塔填料热力性能抽象为式MeTPBVP=f(P0,Twb,Ta,Tw),Tai≤Ta≤Tao,Two≤Tw≤Twi,,再利用公式MeTPBVP=f(P0,Twb,Tao,Twi)计算冷却塔热力性能值。
上述z是代表填料高度的坐标,填料厚度为H,z=0代表填料最底层,z=H代表填料上表面。
本实施例CEF-370不锈钢逆流方塔的热力性能计算系统安装参阅图2。为保证采集运行参数和环境参数的准确性,热力性能计算系统传感器网络应严格遵守冷却塔国标GB/T 7190.1-2008规定安装于冷却塔周围。本实施例监测装置21通过人机界面22监测传感器,传感器网络包括进水流量计、进水温度传感器27、出水温度传感器31、干球温度传感器23、湿球温度传感器24、干球温度传感器23、储水盘温度传感器25、温湿传感器26、大气压传感器28、进水量传感器29、风速传感器32、出风温度传感器33等传感器。
基于两点边值测量下冷却塔填料热力性能计算模型如图3所示,冷却塔内填料区的水和空气间的热质交换可抽象为塔内空气薄膜和水膜间的传热传质过程,布水器洒水后高温冷却水掉落过程中与低温低湿空气接触,并向空气传递热量和蒸发水分,使自身温度Tw下降,同时空气焓值ha增加,湿度X上升,体现为气温Ta上升,其模型焓湿图如图5所示。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (6)
1.一种基于两点边值测量下冷却塔填料热力性能计算方法,其特征在于,所述方法包括以下步骤:
A通过一组求解未知函数导数的四个常微分方程描述填料中的热质传递;
B采集冷却塔现场运行参数和环境参数,确定填料底部即z=0的空气温度Ta和湿度比X,填料顶部即z=H的水温度Tw和质量流率mw边值,H为填料的厚度;
C假设和迭代调整填料底部z=0未知的边界条件,从而使得填料顶部z=H满足已知边界条件;
D基于填料热力性能线性模型建立冷却塔填料热力性能计算模型;
E代入上述计算得到填料顶部z=H的条件,计算冷却塔热力性能。
2.根据权利要求1所述的基于两点边值测量下冷却塔填料热力性能计算方法,其特征在于,所述四个常微分方程为:
3.根据权利要求1所述的基于两点边值测量下冷却塔填料热力性能计算方法,其特征在于,所述运行参数包括:进水量mw、进风量ma、进水温度Twi、出水温度Two、出风温度Tao;所述环境参数包括:大气压力P0、湿球温度Twb、进风干球温度Tai。
4.根据权利要求1所述的基于两点边值测量下冷却塔填料热力性能计算方法,其特征在于,所述步骤B具体包括:利用监测装置的传感网络采集模块实时采集冷却塔现场运行参数和环境参数,确定进水温度Twi,、进水质量流率mwi、进风干球温度Tai、进风湿度比Xi边值。
5.根据权利要求1所述的基于两点边值测量下冷却塔填料热力性能计算方法,其特征在于,所述步骤C具体包括:假设和迭代调整填料底部z=0未知的边界条件 Twb为入口空气的湿球温度,为在出风口温度近似值Tao下计算出的饱和湿度比,从而使得填料顶部z=H满足已知边界条件。
6.根据权利要求1所述的基于两点边值测量下冷却塔填料热力性能计算方法,其特征在于,所述步骤E具体包括:根据步骤C计算的结果和实时测量值,确定填料顶部即z=H气温、水温、大气压力和湿球温度Tao、Twi、P0、Twb;
将冷却塔填料热力性能抽象为式MeTPBVP=f(P0,Twb,Ta,Tw),Tai≤Ta≤Tao,Two≤Tw≤Twi,,再利用公式MeTPBVP=f(P0,Twb,Tao,Twi)计算冷却塔热力性能值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410724779.1A CN104360984B (zh) | 2014-12-03 | 2014-12-03 | 一种基于两点边值测量下冷却塔填料热力性能计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410724779.1A CN104360984B (zh) | 2014-12-03 | 2014-12-03 | 一种基于两点边值测量下冷却塔填料热力性能计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104360984A true CN104360984A (zh) | 2015-02-18 |
CN104360984B CN104360984B (zh) | 2017-06-09 |
Family
ID=52528247
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410724779.1A Active CN104360984B (zh) | 2014-12-03 | 2014-12-03 | 一种基于两点边值测量下冷却塔填料热力性能计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104360984B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107607579A (zh) * | 2017-09-13 | 2018-01-19 | 山东大学 | 一种湿式冷却塔的热力性能在线监测系统及其方法 |
CN107631894A (zh) * | 2017-09-13 | 2018-01-26 | 山东大学 | 高位集水冷却塔的塔内性能参数采集测点布设及测试方法 |
CN109443075A (zh) * | 2018-09-11 | 2019-03-08 | 中国南方电网有限责任公司超高压输电公司广州局 | 一种用于闭式冷却塔的热力性能评估方法及监测系统 |
CN109827747A (zh) * | 2019-03-12 | 2019-05-31 | 上海理工大学 | 用于研究横流填料成型片水流特性的装置及方法 |
CN110260705A (zh) * | 2019-05-17 | 2019-09-20 | 华电电力科学研究院有限公司 | 基于风路系统精准调节的逆流式冷却塔淋水填料测试装置及方法 |
CN112162007A (zh) * | 2020-09-25 | 2021-01-01 | 大连斯频德环境设备有限公司 | 一种冷却塔填料性能实验检测系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103293013A (zh) * | 2013-05-13 | 2013-09-11 | 华南理工大学 | 一种冷却塔热力性能评估方法及系统 |
-
2014
- 2014-12-03 CN CN201410724779.1A patent/CN104360984B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103293013A (zh) * | 2013-05-13 | 2013-09-11 | 华南理工大学 | 一种冷却塔热力性能评估方法及系统 |
Non-Patent Citations (4)
Title |
---|
A.A.BRIN: ""Optimization of the mechanical draft cooling tower performance"", 《HTTPS://WWW.RESEARCHGATE.NET/PUBLICATION/277670336》 * |
V.D. PAPAEFTHIMIOU等: ""Thermodynamic study of the effects of ambient air conditions on the thermal performance characteristics of a closed wet cooling tower"", 《APPLIED THERMAL ENGINEERING》 * |
V.D.PAPAEFTHIMIOU等: ""Thermodynamic study of wet cooling tower performance"", 《INTERNATIONAL JOURNAL OF ENERGY RESEARCH》 * |
刘桂雄等: ""基于TBVMC湿式冷却塔热力性能快速评估方法"", 《中国测试》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107607579A (zh) * | 2017-09-13 | 2018-01-19 | 山东大学 | 一种湿式冷却塔的热力性能在线监测系统及其方法 |
CN107631894A (zh) * | 2017-09-13 | 2018-01-26 | 山东大学 | 高位集水冷却塔的塔内性能参数采集测点布设及测试方法 |
CN107607579B (zh) * | 2017-09-13 | 2020-01-17 | 山东大学 | 一种湿式冷却塔的热力性能在线监测系统及其方法 |
CN109443075A (zh) * | 2018-09-11 | 2019-03-08 | 中国南方电网有限责任公司超高压输电公司广州局 | 一种用于闭式冷却塔的热力性能评估方法及监测系统 |
CN109827747A (zh) * | 2019-03-12 | 2019-05-31 | 上海理工大学 | 用于研究横流填料成型片水流特性的装置及方法 |
CN109827747B (zh) * | 2019-03-12 | 2024-05-31 | 上海理工大学 | 用于研究横流填料成型片水流特性的装置及方法 |
CN110260705A (zh) * | 2019-05-17 | 2019-09-20 | 华电电力科学研究院有限公司 | 基于风路系统精准调节的逆流式冷却塔淋水填料测试装置及方法 |
CN110260705B (zh) * | 2019-05-17 | 2024-04-19 | 华电电力科学研究院有限公司 | 基于风路系统精准调节的逆流式冷却塔淋水填料测试装置的工作方法 |
CN112162007A (zh) * | 2020-09-25 | 2021-01-01 | 大连斯频德环境设备有限公司 | 一种冷却塔填料性能实验检测系统 |
Also Published As
Publication number | Publication date |
---|---|
CN104360984B (zh) | 2017-06-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104360984A (zh) | 一种基于两点边值测量下冷却塔填料热力性能计算方法 | |
CN103293013B (zh) | 一种冷却塔热力性能评估方法及系统 | |
CN104748807B (zh) | 一种基于流量校正的电站主蒸汽流量在线计算方法 | |
Gao et al. | Performance prediction of wet cooling tower using artificial neural network under cross-wind conditions | |
CN104749999B (zh) | 装配湿式冷却塔的汽轮发电机组冷端系统优化运行精确指导系统 | |
CN106022624A (zh) | 一种电-热耦合多能流网络状态估计方法 | |
CN104730107B (zh) | 滨海电厂温排水数学模型表面散热系数实验室观测方法 | |
CN108800303B (zh) | 一种基于模式识别的精细化节能供热方法 | |
CN201697888U (zh) | 一种单一材料围护结构热惰性指标的检测装置 | |
Dang et al. | Crosswind influence on cooling capacity in different zones for high level water collecting wet cooling towers based on field test | |
CN109812902A (zh) | 基于红外热源监测的地铁电气设备房空调节能系统及方法 | |
Young et al. | Performance characteristics of a thermosyphon solar domestic hot water system | |
CN202631283U (zh) | 一种冷却塔热力性能在线智能监测系统 | |
CN102494895B (zh) | 对电站汽轮机组节能及优化的分析方法 | |
Ou et al. | Dynamic modeling and validation of a liquid desiccant cooling and dehumidification system | |
CN104298888B (zh) | 一种基于流量冷量关系模型的风机盘管冷量计量方法 | |
CN116894148B (zh) | 一种冷却塔成雾概率计算方法及基于其的消雾控制系统 | |
CN108595723A (zh) | 一种锅炉暖风器回热量计算方法及装置 | |
CN104361148A (zh) | 一种横流式冷却塔的设计方法 | |
CN111738556A (zh) | 一种基于机头微气候评估风机发电能力的方法 | |
CN207528582U (zh) | 一种铝片结霜性能测试装置 | |
CN107133468B (zh) | 一种间接空冷塔冷却扇段进风量的在线软测量方法 | |
CN108628175B (zh) | 基于空冷系统散热量测量的发电厂优化控制系统及方法 | |
Andria et al. | Design and characterization of solar-assisted heating plant in domestic houses | |
CN108692832B (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 |