王艺儒,詹淑慧(北京建筑大学 环境与能源工程学院,北京100044)
粗煤气是煤在焦炉中高温干馏产生的混合气,蕴含着大量的显热[1-2]。粗煤气热量回收方法主要有上升管汽化换热、导热油夹套换热、热管式换热、余热锅炉换热等[3-7]。余热锅炉换热技术投资比较少,能有效克服结焦,稳定性和回收效率比较高。本文在显热换热前提下,采用Aspen Plus软件对余热锅炉换热面积进行模拟计算,确定换热管束结构,采用关联式法对余热锅炉换热面积进行校核计算。
① 工艺流程
余热锅炉工艺流程见图1。由焦炉排出的高温粗煤气进入余热锅炉,加热余热锅炉进水,加热后的水可供居民、工业用户使用。
图1 余热锅炉工艺流程
② 粗煤气组成
粗煤气主要成分为干煤气、其他杂质(包括水蒸气)、煤焦油。干煤气(组成见表1)是粗煤气中的轻组分,约占粗煤气质量的49%。其他杂质约占粗煤气质量的40%,且大部分为水蒸气。煤焦油约占粗煤气质量的11%,一般分为两部分[8]:一部分是确定的轻组分;另一部分为无法确定的组分,一般称为沥青。
表1 干煤气的组成
③ 物料物性参数
余热锅炉可视为管壳式换热器,粗煤气走管程,水走壳程。设定粗煤气流速为20 m/s,进口温度为400 ℃,出口温度为200 ℃,定性温度为进出口温度的算术平均值(300 ℃)。设定水的流速为1.0 m/s,进口温度为30 ℃,出口温度为80 ℃,定性温度为进出口水温的算术平均值(55 ℃),物性参数见表2[9]330。采用Aspen Plus物性估算出定性温度下的粗煤气的物性参数,见表2。估算时,将粗煤气视为干煤气,组成见表1。物性方法分为状态方程法、活度系数法,状态方程法主要用于对气体物性分析过程,常见的状态方程有Ideal、SRK、PR等,SRK方程可以比较精确地推算出混合气体物料的真实物性参数。因此,物性方法选取SRK方程。
表2 定性温度下的粗煤气、水的物性参数
① 初步设计
选用规格为D57×3.5的碳钢换热管,选取Aspen Plus软件适用于管壳式换热器的单元模块(HeatX)中的简捷法模型(shortcut)。该模型不需要换热器结构或几何尺寸,仅通过输入少量已知量即可对换热器进行初步设计。输入粗煤气的进出口温度、水的进口温度,设定水的质量流量为10 t/h。运行模拟程序,得出余热锅炉热负荷为644.71 kW,水出口温度为80 ℃,换热面积(以下称为模拟换热面积)为34.71 m2。
管壳式换热器总管数N的计算式为[10]:
(1)
(2)
式中N——管壳式换热器总管数
qm——粗煤气的质量流量,kg/s
u——粗煤气的流速,m/s
ρ——定性温度下粗煤气的密度,kg/m3
do——换热管的外直径,m
Φ——管壳式换热器热负荷,W
cp——定性温度下粗煤气的比定压热容,J/(kg·K)
Δt——粗煤气进出口温差,℃
换热管单管长度L的计算式为:
(3)
(4)
式中L——换热管单管长度,m
A——管壳式换热器换热面积,m2
K——管壳式换热器总传热系数,W/(m2·K)
Δtm——对数平均温差,℃
将已知参数代入式(1)、(2),可计算得到管壳式换热器总管数为50.71,进而由式(3)、(4)可计算得到换热管单管长度为3.94 m。根据GB/T 28712.2—2012《固定管板式换热器型式与基本参数》第4.3.1条规定,换热管单管长度选取4.5 m。
② 换热管束结构设计
换热管束排布方式选取采用最常见的正三角形排列方式(见图2),为满足排布方式需要,管壳式换热器总管数确定为52。
图2 换热管束排布方式
管内表面传热系数hin的关联式为[9]182:
(5)
(6)
(7)
式中hin——管内表面传热系数,W/(m2·K)
λ——定性温度下粗煤气的热导率,W/(m·K)
din——换热管的内直径,m
Re——粗煤气的雷诺数
Pr——粗煤气的普朗特数
n——指数,流体被加热时取0.1,流体被冷却时取0.3,这里取0.3
μ——定性温度下粗煤气的动力黏度,Pa·s
将已知参数代入式(5)~(7),可计算得到管内表面传热系数为114.18 W/(m2·K)。
壳程为强迫对流换热,管外表面传热系数ho的关联式为[9]183:
(8)
(9)
(10)
式中ho——管外表面传热系数,W/(m2·K)
λw——定性温度下水的热导率,W/(m·K)
de——壳程当量直径,m
uw——水的流速,m/s
ρw——定性温度下水的密度,kg/m3
μw——定性温度下水的动力黏度,Pa·s
Prw——水的普朗特数
Nc——位于管束中心线上的管数,为8
cp,w——定性温度下水的比定压热容,J/(kg·K)
将已知参数代入式(8)~(10),可计算得到管外传热系数为1 673.29 W/(m2·K)。
忽略换热管壁的热阻,管壳式换热器总传热系数K的计算式为[9]183:
(11)
式中Ro——管外单位面积污垢热阻,m2·K/W
Rin——管内单位面积污垢热阻,m2·K/W
管外单位面积污垢热阻、管内单位面积污垢热阻按GB/T 151—2014《热交换器》附录E分别选取为3.5×10-4、1.76×10-4m2·K/W。将已知参数代入式(11),可计算得到管壳式换热器总传热系数为89.87 W/(m2·K)。进而由式(4)计算得到管壳式换热器换热面积(以下称为计算换热面积)为31.41 m2,与模拟换热面积(为34.71 m2)十分接近,说明模拟结果可信。
由Aspen Plus软件模拟计算得到的余热锅炉换热面积为34.71 m2,由关联式法计算得到的换热面积为31.41 m2,二者十分接近,模拟结果可信。