CO2汽提塔的计算与分析

廖容波，刘利利

(北京众联盛化工工程有限公司 北京 100101)

CO2汽提塔的计算与分析

廖容波，刘利利

(北京众联盛化工工程有限公司 北京 100101)

利用液体穿孔的流量公式，推导出CO2汽提塔液体分布器上的液体分布孔流速与液位高度的关联公式，可计算出汽提塔单管负荷。根据湿壁塔内液膜厚度公式，计算出尿素-甲铵液在不同类型CO2汽提塔中的停留时间。通过热力学计算，归纳总结了不同生产负荷下CO2汽提塔的热负荷和传热系数。计算结果表明，理论计算值与实际生产数据基本吻合。

CO2汽提塔 尿素 液体分布器 热负荷

CO2汽提法尿素装置的高压设备由合成塔、高压甲铵冷凝器、高压洗涤器和汽提塔组成。NH3和CO2在高压甲铵冷凝器中快速合成氨基甲酸铵(以下简称甲铵)，然后进入合成塔中缓慢脱水转化成尿素，未脱水的甲铵随尿液进入汽提塔中等压分解为NH3和CO2后循环回收利用，其中CO2汽提塔是影响尿素装置操作弹性的重要因素之一，直接关系到装置的稳定运行。

1 设备简介及规格

CO2汽提塔为降膜式单管程列管换热器，管内为工艺介质，管外的饱和蒸汽为甲铵分解提供所需的热量。出合成塔的反应液在180～185 ℃下进入汽提塔上部，经液体分布器均布后进入列管，成膜状沿管壁下流，CO2从汽提塔下部逆流汽提反应液，促使未转化的甲铵分解为NH3和CO2，分解气和原料CO2气体从汽提塔塔顶逸出后进入高压甲铵冷凝器。

目前，CO2汽提塔已全部实现国产化，主流设计有3种规格，匹配的尿素装置产能分别为300，520和800 kt/a，其设备参数见表1。

2 液体分布器

根据汽提原理，尿液必须均匀分配至每根换热管中并在内壁形成一层液膜，才会有较高的汽提效率，因此，液体分布器是汽提塔的关键部件。如果液体分布不均，不仅会造成汽提效率下降，而且会影响防腐空气的均匀通过，造成列管局部过热和缺氧腐蚀[1]。

表1 CO2汽提塔设备参数

项 目产能/(t·d-1)外形尺寸/mm换热面积/m2列管数/根液体分布孔直径/mm国产520kt/a1740Φ2400×1259016802875Φ2.5荷型480kt/a162012252600Φ2.3传统300kt/a1000Φ1900×1303010121733Φ2.5湖北三宁公司1450Φ1900×1303010121733Φ2.5国产800kt/a3200Φ2900×1457024544200

注：1)列管规格均为Φ31 mm×3 mm×6 000 mm，单管液体分布孔数均为3个

液体分布器由升气管和液体分布头组成，液体分布器与分布头焊接在一起并插在汽提管上，中间由聚四氟乙烯垫圈进行密封。每个液体分布器上设有3个Φ2.3～2.5 mm互成120°的液体分布孔，均匀排布，其结构示意见图1。

图1 液体分布器结构示意

升气管类似喇叭状，上部直径为Φ25 mm，下部直径为Φ29 mm，焊接于液体分布头上；顶部有

Φ8.0 mm的升气孔，其作用是限流，保持液体分布器内外有一定的压差，进而可维持液体分布器外侧有一定液位，确保尿液在一定的压头下能够均匀分布。液体分布器总高622 mm，其中分布孔距底端35 mm。目前，通用CO2汽提塔液体分布器上的液体分布孔有Φ2.3 mm和Φ2.5 mm 2种规格。

3 计算

3.1 单管负荷

尿素-甲铵液进入受液槽后，沿圆周方向流动，平稳地流入分布器空间。液体在分布器小孔上保持一定的液位，保持液体均匀进入小孔。液体穿孔的流量按下述公式进行计算[2]：

(1)

(2)

(3)

式中：Q——液体流量，m3/s；C0——孔流系数，与小孔的形状和液体的雷诺数有关；

A0——流通面积，m2；

n——孔数；

H——液体高度，m；

d0——孔径，m；

u——液体经过小孔的流速，m/s。

由于汽提塔液体分布器液体分布孔表面粗糙度为3.2，经统计分析，C0取下限0.6。结合式(1)和式(2)，可以推导出液体分布器液体分布孔的流速与受液槽中尿素-甲铵液高度有关：

(4)

根据上述公式，结合2种液体分布器规格，推算出的CO2汽提塔单管理论负荷见表2。

表2 CO2汽提塔单管理论负荷

液体分布孔直径/mm液体密度/(kg·m-3)液位/mm正常最大小孔流速/(m·s1)正常最大单管负荷/(kg·h-1)正常最大Φ2.510503005871.452.0080.9113.0Φ2.310503005501.451.9768.392.8

根据CO2汽提法尿素装置的理论设计和现场反馈，CO2汽提塔单管负荷见表3。

从表2和表3可知，CO2汽提塔单管负荷的理论计算值与实际反馈值基本吻合。由于CO2汽提法尿素装置的生产负荷不能低于设计值的70%，则可计算出汽提塔液体分布器小孔流速最

表3 CO2汽提塔单管设计和实际负荷

项 目生产能力/(t·d-1)列管数量/根进塔液体量/(kg·h-1)单管负荷/(kg·h-1)国产520kt/a1740287522961079.90荷型480kt/a1620260020160077.50传统300kt/a1000173313127775.75湖北三宁公司14501733191342110.40

小值和最低液位分别为1.0 m/s和142 mm。当液体分布器上的液位低于142 mm时，汽提塔中的尿素-甲铵液因液位梯度原因会导致分布不均，造成有些加热管出现干管和腐蚀加剧的现象。

由CO2汽提工艺的操作手册和生产现场反馈可知，汽提塔的操作液位高度应保持在200～400 mm[3]，可计算出单管负荷分别为66.2 kg/h和93.5 kg/h，对应的弹性负荷为82.6%～117.0%，与实际生产数据基本吻合。

3.2 尿液停留时间

CO2汽提塔的壳程采用2.0 MPa蒸汽加热，为尿液中甲铵分解提供热量，但高温会导致尿素水解产生缩二脲。正常情况下，汽提塔后的尿液中缩二脲质量分数增加0.2%～0.3%[3]。

尿液在CO2汽提塔内的停留时间(理论设计值<1 min)可根据湿壁塔内液膜厚度按下式计算(表4)：

(5)

式中：m——液膜厚度，cm； w——单位长度周边的液体质量速度，g/(s·cm)；

μ——液体黏度，mPa·s；

r——液体密度，g/cm3；

g——重力加速度，cm/s2。

表4 CO2汽提塔单管液膜厚度理论计算值

项 目液体密度/(g·cm-3)进口出口液体黏度/(mPa·s)进口出口液体质量速度/(g·s-1·cm-1)进口出口液膜厚度/cm进口出口平均正常情况下停留时间/min国产520kt/a0.991.130.7520.9132.901.700.1900.1550.1730.76荷型480kt/a0.991.130.7520.9132.801.600.1870.1550.1710.77传统300kt/a0.991.130.7520.9132.681.550.1850.1500.1680.76湖北三宁公司0.991.130.7520.9133.912.260.2090.1700.1900.87

注：1)停留时间=汽提塔液膜总体积/[0.5×(进塔液体体积流量+出塔液体体积流量)]

由表4可知：尿素-甲铵液在CO2汽提塔内的停留时间均不超过1 min；湖北三宁公司的尿素装置为传统国产300 kt/a设计，其最大产量为1 450 t/d，但其尿液在汽提塔中停留时间仅有0.87 min，说明该产量下的停留时间基本上不会对汽提塔出口尿液中的缩二脲含量有太大影响。

3.3 热负荷和传热系数

汽提塔实际上是1台多管降膜式加热器，由管外2.0 MPa左右的中压饱和蒸汽供给热量，促使管内甲铵分解反应不断进行。CO2汽提塔单管设计和实际热负荷见表5。

表5 CO2汽提塔单管设计和实际热负荷

项 目生产能力/(t·d-1)换热面积/m2平均温差/℃热负荷/(kJ·h-1)传热系数/(kJ·h-1·m-2·℃-1)国产520kt/a17401680.040.01.33×108473.0荷型480kt/a16201224.638.01.23×108633.0传统300kt/a10001012.738.00.76×108472.6湖北三宁公司14501012.740.91.10×108636.6

由表5可看出：汽提塔的传热系数可满足生产负荷的弹性要求，通过适当提高汽提塔壳侧蒸汽压力，蒸汽温度随之提高，从而可提高CO2汽提塔的处理能力。饱和蒸汽的温度受管侧使用材料限制，而且温度过高的蒸汽会导致尿液温度偏高，从而使缩二脲产生量和尿素的水解量相应增加，因此，应根据尿素装置的产量和汽提塔负荷来调节控制汽提塔的壳侧蒸汽压力。

4 结语

(1)根据CO2汽提塔的理论设计和实际操作条件，利用液体穿孔的流量公式，推导出小孔流速与液位高度的关联公式，由液位高度可准确计算出汽提塔液体分布孔流速和单管负荷，进而可进一步推算出汽提塔的生产负荷。

(2)利用湿壁塔内液膜厚度的公式，计算出尿素-甲铵液在不同类型CO2汽提塔中的停留时间，其计算结果均满足设计要求。

(3)通过热力学计算，归纳总结了不同生产负荷下CO2汽提塔的热负荷和传热系数。

(4)单管负荷、停留时间和热负荷的理论计算值与实际生产数据基本吻合，其计算结果对设计和实践生产有较好的指导意义。

[1] 邸彩霞，崔海生，李进.汽提塔带油原因分析与对策[J].化肥设计，2006(5)：37- 39.

[2] 王汉松.石油化工设计手册：第3卷 化工单元过程[M].北京：化学工业出版社，2002.

[3] 大连工学院.日产千吨合成氨厂——二氧化碳气提法尿素生产工艺[M].北京：石油化学工业出版社，1978.

Calculation and Analysis of CO2Stripper

LIAO Rongbo, LIU Lili

(Beijing Unity Engineering CO., Ltd. Beijing 100101)

Use of flux formula for orifice flow, correlation formula has been deduced of flow velocity in liquid distribution hole at liquid distributor of CO2stripper and height of liquid level, the load of individual pipe of the stripper can be calculated. By formula of liquid film thickness in wetted wall column, the residence time of urea-ammonium carbamate solution in different CO2strippers can be calculated. By thermodynamic calculation, heat load and heat transfer coefficient of CO2stripper at different production load are generalized and summarized. Calculation results show that theoretical calculated value basically agrees with actual production data.

CO2stripper urea liquid distributor heat load

廖容波(1981—)，男，工程师，从事化工工艺设计工作；liao12837@163.com。

TQ441.41

A

1006- 7779(2016)06- 0062- 03

2015- 04- 27)