CN102788517A - 超大型逆流式自然通风湿式冷却塔的设计方法及系统 - Google Patents

Abstract

本发明提供一种超大型逆流式自然通风湿式冷却塔的设计方法及系统，该方法包括以下步骤：获取进风口相对高度ε以及淋水填料阻力系数ξ_f；根据所述进风口相对高度ε、淋水填料阻力系数ξ_f并采用如下公式计算进风口区域阻力系数ξ_a：

根据所述进风口区域阻力系数进行冷却塔的最终特性计算，获得冷却塔的运行特性数据；根据所述运行特性数据进行冷却塔的设计。本发明的方法及系统能够有效掌握冷却塔的实际运行状况，对机组运行稳定性和安全性提供了可靠保障，有效避免了由于计算选型的冷却塔偏小所带来的安全性问题以及计算选型的冷却塔偏大所带来的初期投资浪费问题，因此具有较好的市场应用前景。

Description

超大型逆流式自然通风湿式冷却塔的设计方法及系统

技术领域

本发明涉及一种超大型逆流式自然通风湿式冷却塔的设计方法以及一种超大型逆流式自然通风湿式冷却塔的设计系统。

背景技术

逆流式自然通风湿式冷却塔(以下简称自然塔)广泛的应用于国民经济的诸多部门，包括电力、石油、化工等，其作用是将挟带废热的冷却水在塔内与空气直接接触进行热交换，使废热传输给空气并散入大气。这种类型冷却塔如图1所示，通风筒常采用双曲线形，用钢筋混凝土浇制，塔筒底部为进风口，空气从进风口进入塔体，穿过填料下的雨区，和热水流动成相反方向流过填料，再从塔筒出口流出。

自然塔淋水面积是指冷却塔内“填料区”顶部的断面面积，按淋水面积的大小，冷却塔可初步划分为以下几种：

小型塔：A＜4000m²；

中型塔：4000m²≤A＜7000m²；

大型塔：7000m²≤A＜12000m²；

超大型塔：A≥12000m²(塔底部直径D＞110m、进风口高度h＞11m)。

自然塔空气动力计算是冷却塔进行工艺设计的核心，也是再循环供水系统优化设计和热力设计的依据，而其中阻力计算又是空气动力计算的重要组成部分。自然通风逆流式冷却塔的阻力系数计算一共包括了三个部分，其中比较重要的一部分是进风口区域除雨滴阻力外的人字柱、气流转弯、支柱、填料、收水器与配水装置阻力系数(简称进风口区域阻力系数)计算。可见，进风口区域阻力系数的计算将直接影响冷却塔的设计，尤其是对于底部直径大于110m、进风口高度大于11m、淋水面积大于12000m²的超大型塔来说，如果进风口区域阻力系数的计算结果不准确，将导致最终设计的冷却塔的设计工况与实际运行工况有偏差，即计算选型的冷却塔有可能出现偏大或者偏小的情况：偏大就意味着建设规模过大，造成极大浪费，还可能会使机组出现过冷现象；偏小就会出现冷却能力不足的现象，可能会导致机组发电功率不足，甚至影响运行稳定性和安全性，出现运行事故。

鉴于进风口区域阻力系数对超大型自然塔的设计所具有的重要地位，因此有必要进行自然塔的进风口区域阻力系数研究。传统的超大型逆流式自然通风湿式冷却塔的设计方法一般采用如下公式来计算进风口区域阻力系数：

${\mathrm{\ξ}}_a=(1-3.47\mathrm{\ϵ}+3.65{\mathrm{\ϵ}}^2)(85+2.51{\mathrm{\ξ}}_f-0.206{\mathrm{\ξ}}_f^2+0.00962{\mathrm{\ξ}}_f^3)---\left(1\right)$

式中，ξ_a为进风口区域阻力系数；ξ_f为淋水填料阻力系数(含配水系统收水器等阻力系数)；ε为进风口相对高度(即进风口面积与塔壳底面积比)。

从上面的公式(1)可知，进风口区域阻力系统的计算公式中的变量为淋水填料阻力系数(含收水器和配水系统)和相对进风口高度(进风口面积与塔壳底面积比)，该计算公式对于目前已有常规自然通风冷却塔阻力计算是基本适用的，并且已经被冷却塔设计规范所采用。但由于当时冷却塔淋水规模及试验条件等原因限制，该公式在适用范围上还是具有一定的局限性：

一、公式(1)中，淋水填料阻力系数变化范围为10～20，然而超大型冷却塔的淋水填料高度将较常规冷却塔有所提高，比如填料高度可能达到20m或以上，超出公式(1)的适用范围；

二、公式(1)中的进风口相对高度的适用范围为0.35～0.45，常规冷却塔的运行中冷却塔的相对进风口一般为0.36；然而超大型冷却塔的面积大约为同样装机容量火电冷却塔的80％～100％，进风口的相对高度仍采用0.36势必造成冷却塔的进风口高度较高，增大了厂用电。而且另一方面核电冷却塔高度较高，填料高度有所增大，进风口的高度可相对低些，比如德国的依萨核电16500m²的逆流式自然通风冷却塔的相对进风口高度为0.3。所以，公式(1)的适用性再次受到局限。

近年来，随着电力行业尤其是核电的迅速发展，冷却塔淋水面积不断增大，内陆核电的冷却水系统配套的冷却塔的淋水面积通常都大于12000m²、进风口高度通常大于11m，对于这种规模的超大型冷却塔，进风口区域阻力系数发生较大的变化，再采用传统的计算方法和计算手段，可能或造成比较大的偏差甚至错误，因此研究新的冷却塔进风口区域阻力数值计算技术是十分必要的。

发明内容

基于此，有必要针对上述传统计算进风口区域阻力系数的方法具有一定局限性，导致计算选型的冷却塔偏大或者偏小的问题，提供一种超大型逆流式自然通风湿式冷却塔的设计方法及系统。

一种超大型逆流式自然通风湿式冷却塔的设计方法，包括以下步骤：

获取进风口相对高度ε以及淋水填料阻力系数ξ_f；

根据所述进风口相对高度ε、淋水填料阻力系数ξ_f并采用如下公式计算进风口区域阻力系数ξ_a： ${\mathrm{\ξ}}_a=(76.65-261.33\mathrm{\ϵ}+264.85{\mathrm{\ϵ}}^2)e^{(-0.0217+0.230\mathrm{\ϵ}-0.2191{\mathrm{\ϵ}}^2){\mathrm{\ξ}}_f}-1.5;$

根据所述进风口区域阻力系数进行冷却塔的最终特性计算，获得冷却塔的运行特性数据；

根据所述运行特性数据进行冷却塔的设计。

一种超大型逆流式自然通风湿式冷却塔的设计系统，包括：

参数获取模块，用于获取进风口相对高度ε以及淋水填料阻力系数ξ_f；

进风口区域阻力系数计算模块，用于根据所述进风口相对高度ε、淋水填料阻力系数ξ_f并采用如下公式计算进风口区域阻力系数ξ_a： ${\mathrm{\ξ}}_a=(76.65-261.33\mathrm{\ϵ}+264.85{\mathrm{\ϵ}}^2)e^{(-0.0217+0.230\mathrm{\ϵ}-0.2191{\mathrm{\ϵ}}^2){\mathrm{\ξ}}_f}-1.5;$

运行特性数据计算模块，用于根据所述进风口区域阻力系数进行冷却塔的最终特性计算，获得冷却塔的运行特性数据；

设计模块，用于根据所述运行特性数据进行冷却塔的设计。

由以上方案可以看出，本发明的超大型逆流式自然通风湿式冷却塔的设计方法及系统，为底部直径大于110m、进风口高度大于11m、淋水面积大于12000m²的超大型逆流式自然通风湿式冷却塔进风口区域阻力数值计算提供了解决方法，避免了继续沿用传统的数值计算技术进行冷却塔的进风口区域阻力系数计算所导致的冷却塔的设计工况与实际运行工况有所偏离的问题，具有更加广泛的适用性，能够为超大型冷却塔设计中相关阻力计算提供参考和依据，并最终满足超大型逆流式自然通风湿式冷却塔的设计需要。另外本发明的方法及系统有效掌握了冷却塔的实际运行状况，对机组运行稳定性和安全性提供了可靠保障，有效避免了由于计算选型的冷却塔偏小所带来的安全性问题以及计算选型的冷却塔偏大所带来的初期投资浪费问题，因此具有较好的市场应用前景。

附图说明

图1为逆流式自然通风湿式冷却塔结构示意图；

图2为一种超大型逆流式自然通风湿式冷却塔的设计方法流程示意图；

图3为一种超大型逆流式自然通风湿式冷却塔的设计系统结构示意图。

具体实施方式

下面结合附图以及具体的实施例，对本发明的技术方案作进一步的描述。

参见图1所示，一种超大型逆流式自然通风湿式冷却塔的设计方法，包括以下步骤：

步骤S101，获取进风口相对高度(记为ε)以及淋水填料阻力系数(记为ξ_f)，然后进入步骤S102。

步骤S102，根据所述进风口相对高度(ε)、淋水填料阻力系数(ξ_f)并采用如下公式计算进风口区域阻力系数(记为ξ_a)：

${\mathrm{\ξ}}_a=(76.65-261.33\mathrm{\ϵ}+264.85{\mathrm{\ϵ}}^2)e^{(-0.0217+0.230\mathrm{\ϵ}-0.2191{\mathrm{\ϵ}}^2){\mathrm{\ξ}}_f}-1.5---\left(2\right)$

下面描述本发明进风口区域阻力系数的计算公式的推导过程，其包括了如下几个步骤：

一、构建模拟试验装置

装置试验在试验室进行，需能进行“不同淋水面积、不同进风口高度、不同填料阻力系数、不同填料断面风速工况”组合条件下的进风口区域阻力特性试验。可以通过建立一个冷态模型试验系统来了解超大型冷却塔相关阻力特性，模型将对塔内各主要结构进行模拟，包括淋水填料、进风口高度、塔筒支柱以及塔内梁柱等。具体构建模拟试验装置的过程本发明不作赘述。

二、多工况的试验测试及数据整理

分步骤进行不同淋水面积阻力特性试验、不同塔体结构组成部分阻力特性试验和进风口高度变化对于塔内阻力性能影响试验，然后对试验数据进行整理。具体细分为以下两个步骤：

(1)、试验工况

本试验重点研究超大型自然通风冷却塔不同淋水面积规模的进风口区域(不含雨区)阻力特性。模型需模拟和测试10000m²、16000m²和20000m²原型塔进风口区域阻力特征，同时还比较和分析了塔筒支柱、塔内梁柱等局部结构的阻力影响，以及不同进风口高度与塔进风口区域阻力系数的变化关系。试验中，塔内淋水填料分别采用三种孔板模拟。据此，试验工况布置见表1，共计10组。进风口高度影响试验时，模型淋水面积不变，通过调节进风口高度来实现进风面积与淋水面积的变化。

表1试验工况表

表2淋水填料模拟孔板设计参数

孔板编号	孔板1	孔板2	孔板3
				设计阻力系数	10	25	40
板厚L(mm)	8	8	8
				孔径(mm)	10	10	10
孔板断面直径(m)	1.04	1.04	1.04
				过流面积(m2)	0.28	0.19	0.16
孔隙率(过流面积/断面面积)	0.27	0.18	0.16
				模型通风量(103m3/h)	2～20	2～20	2～20
断面风速(m/s)	0.7～7.0	0.7～7.0	0.7～7.0

断面雷诺数(105)	0.4～4.0	0.4～4.0	0.4～4.0
				孔口雷诺数(104)	0.11～1.1	0.16～1.6	0.20～2.0

(2)、试验数据整理

a)阻力测试

喉部断面各测点阻力等于该点与塔外大气的全压差值，取各点算数平均值作为塔进口至喉部断面的平均阻力，可以用下式计算：

$\mathrm{\ΔP}=\frac{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}\mathrm{\Δ}{P}_i}{n}---\left(3\right)$

式中：

ΔP---塔进口至喉部断面平均阻力，Pa；

ΔP_i---喉部断面测点与外界全压差，Pa。

b)阻力系数

冷却塔的阻力一般可表示为填料断面处气流速度头与冷却塔总阻力系数的积，即如下式：

${\mathrm{\ΔP}}_R={\mathrm{\ξ\ρ}}_m\frac{V_f^2}{2}---\left(4\right)$

式中：ΔP_R---冷却塔的气流阻力，Pa；

ξ---冷却塔的总阻力系数；

ρ---塔内的平均空气密度，kg/m³；

V_f---填料断面平均风速，m/s。

冷却塔中气流受到的阻力主要是受进风口及气流转向、淋水雨区、支撑系统(塔筒支柱、塔内梁柱)、淋水填料、配水系统、除水器、塔出口几部分综合影响的结果。一般将前6部分作为一个量来考虑，公式(4)中的总阻力系数可写为：

$\mathrm{\ξ}={\mathrm{\ξ}}_{1-6}+\mathrm{\α}{\left(\frac{F_f}{F_o}\right)}^2---\left(5\right)$

式中：

ξ---塔总阻力系数；

ξ_1-6---进风口及气流转向、淋水雨区、支撑系统、填料、配水系统、除水器的阻力系数总和；

F_o---冷却塔出口面积，m²；

F_f---淋水面积，m²；

α---出口流速不均匀系数，一般取为1。

其中，本试验未模拟淋水雨区及喉部以上塔出口阻力。填料、配水系统及除水器(塔芯材料)的总阻力由孔板来综合模拟。因此，模型中从进风口到喉部断面的阻力系数即为塔进风口区域阻力系数，如下：

${\mathrm{\ξ}}_t={\mathrm{\ξ}}_1+{\mathrm{\ξ}}_3+{\mathrm{\ξ}}_6=\frac{\mathrm{\ΔP}}{\frac{1}{2}{\mathrm{\ρV}}_f^2}---\left(6\right)$

式中：

ξ_t---塔进风口区域综合阻力系数；

ξ₁---塔进风口及气流转向阻力系数；

ξ₃---支撑系统阻力系数；

ξ₆---除水器(塔芯材料)阻力系数；

ΔP---从塔进口至喉部断面阻力，Pa。

其他符号同前。

三、数值计算公式的提炼

上述工况的模型试验结果见下表3。

表3不同相对进风高度不同填料阻力系数进风口区域阻力系数试验结果

根据上表不同进风口与塔底壳面积比值下的塔进风口区域阻力系数测试结果，即可拟合推导出进风口区域(不含雨区)阻力系数与淋水填料阻力系数及进风口相对高度(进风口面积与塔底壳面积比值)之间的关系表达式如下：

${\mathrm{\ξ}}_a=(76.65-261.33\mathrm{\ϵ}+264.85{\mathrm{\ϵ}}^2)e^{(-0.0217+0.230\mathrm{\ϵ}-0.2191{\mathrm{\ϵ}}^2){\mathrm{\ξ}}_f}-1.5$

式中：ε为进风口相对高度(进风口面积与塔底壳面积比)；

ξ_f为淋水填料阻力系数(含收水器与配水系统)。

步骤S103，根据所述进风口区域阻力系数进行冷却塔的最终特性计算，获得冷却塔的运行特性数据，然后进入步骤S104。

由于进风口区域阻力计算是全塔阻力计算的一部分，因此在得到进风口区域阻力系数后即可进行全塔阻力计算，然后利用冷却塔总阻力等于总抽力的方程式进行冷却塔空气动力学计算，再结合冷却塔热力计算进行冷却塔的最终特性计算，就可以获得冷却塔的运行特性数据(出塔水温等)。上述过程属于公知技术，本发明中不作赘述。

步骤S104，根据所述运行特性数据(如出塔水温等)进行冷却塔的设计。可根据计算得到的精确的运行特性数据设计选择得到合适的冷却塔型号，避免计算选型的冷却塔偏大或偏小的情况出现。

作为一个较好的实施例，本发明中淋水填料阻力系数取值范围则可以为10～25，进风口面积与塔底壳面积比值取值范围可以0.30～0.50，具有更加广泛的适用性，进一步解决了传统计算方法中淋水填料阻力系数以及进风口相对高度的范围具有局限性的问题。

另外，与上述一种超大型逆流式自然通风湿式冷却塔的设计方法相对应的，本发明还提供一种超大型逆流式自然通风湿式冷却塔的设计系统，如图3所示，包括：

参数获取模块101，用于获取进风口相对高度ε以及淋水填料阻力系数ξf；

进风口区域阻力系数计算模块102，用于根据所述进风口相对高度ε、淋水填料阻力系数ξ_f并采用如下公式计算进风口区域阻力系数ξ_a： ${\mathrm{\ξ}}_a=(76.65-261.33\mathrm{\ϵ}+264.85{\mathrm{\ϵ}}^2)e^{(-0.0217+0.230\mathrm{\ϵ}-0.2191{\mathrm{\ϵ}}^2){\mathrm{\ξ}}_f}-1.5;$

运行特性数据计算模块103，用于根据所述进风口区域阻力系数进行冷却塔的最终特性计算，获得冷却塔的运行特性数据；

设计模块104，用于根据所述运行特性数据进行冷却塔的设计。

作为一个较好的实施例，所述进风口相对高度ε的取值范围可以为0.30～0.50。

作为一个较好的实施例，所述淋水填料阻力系数ξ_f的取值范围可以为10～25。

本发明的一种超大型逆流式自然通风湿式冷却塔的设计系统的其它技术特征与上述一种超大型逆流式自然通风湿式冷却塔的设计方法相同，此处不予赘述。

通过以上方案可以看出，本发明的超大型逆流式自然通风湿式冷却塔的设计方法及系统，为底部直径大于110m、进风口高度大于11m、淋水面积大于12000m²的超大型逆流式自然通风湿式冷却塔进风口区域阻力数值计算提供了解决方法，避免了继续沿用传统的数值计算技术进行冷却塔的进风口区域阻力系数计算所导致的冷却塔的设计工况与实际运行工况有所偏离的问题，具有更加广泛的适用性，能够为超大型冷却塔设计中相关阻力计算提供参考和依据，并最终满足超大型逆流式自然通风湿式冷却塔的设计需要。另外本发明的方法及系统有效掌握了冷却塔的实际运行状况，对机组运行稳定性和安全性提供了可靠保障，有效避免了由于计算选型的冷却塔偏小所带来的安全性问题以及计算选型的冷却塔偏大所带来的初期投资浪费问题，因此具有较好的市场应用前景。

以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。

Claims (6)

1.一种超大型逆流式自然通风湿式冷却塔的设计方法，其特征在于，包括以下步骤：

获取进风口相对高度ε以及淋水填料阻力系数ξ_f；

根据所述进风口相对高度ε、淋水填料阻力系数ξ_f并采用如下公式计算进风口区域阻力系数ξ_a： ${\mathrm{\ξ}}_a=(76.65-261.33\mathrm{\ϵ}+264.85{\mathrm{\ϵ}}^2)e^{(-0.0217+0.230\mathrm{\ϵ}-0.2191{\mathrm{\ϵ}}^2){\mathrm{\ξ}}_f}-1.5;$

根据所述进风口区域阻力系数进行冷却塔的最终特性计算，获得冷却塔的运行特性数据；

根据所述运行特性数据进行冷却塔的设计。

2.根据权利要求1所述的超大型逆流式自然通风湿式冷却塔的设计方法，其特征在于，所述进风口相对高度ε的取值范围为0.30～0.50。

3.根据权利要求1或2所述的超大型逆流式自然通风湿式冷却塔的设计方法，其特征在于，所述淋水填料阻力系数ξ_f的取值范围为10～25。

4.一种超大型逆流式自然通风湿式冷却塔的设计系统，其特征在于，包括：

参数获取模块，用于获取进风口相对高度ε以及淋水填料阻力系数ξ_f；

进风口区域阻力系数计算模块，用于根据所述进风口相对高度ε、淋水填料阻力系数ξ_f并采用如下公式计算进风口区域阻力系数ξ_a： ${\mathrm{\ξ}}_a=(76.65-261.33\mathrm{\ϵ}+264.85{\mathrm{\ϵ}}^2)e^{(-0.0217+0.230\mathrm{\ϵ}-0.2191{\mathrm{\ϵ}}^2){\mathrm{\ξ}}_f}-1.5;$

运行特性数据计算模块，用于根据所述进风口区域阻力系数进行冷却塔的最终特性计算，获得冷却塔的运行特性数据；

设计模块，用于根据所述运行特性数据进行冷却塔的设计。

5.根据权利要求4所述的超大型逆流式自然通风湿式冷却塔的设计系统，其特征在于，所述进风口相对高度ε的取值范围为0.30～0.50。

6.根据权利要求4或5所述的超大型逆流式自然通风湿式冷却塔的设计系统，其特征在于，所述淋水填料阻力系数ξ_f的取值范围为10～25。

Images