石油与天然气化工  2015, Vol. 44 Issue (1): 12-16
磷石膏、NH3联合矿化天然气净化厂尾气中CO2的流程模拟
林婷 , 魏文韫 , 刘凌岭 , 余徽     
四川大学化学工程学院
摘要:针对CCUS新工艺——磷石膏、NH3联合矿化天然气净化厂排放尾气中CO2,应用流程模拟方法,在物性估算、模型建立的基础上,对磷石膏矿化CO2反应的工艺条件进行了优化,提出反应温度为60 ℃、反应时间为100 min的优化方案。对比优化后工艺条件与现有工艺条件的模拟结果可知,优化后尾气中CO2、NH3、磷石膏排放量比现有工艺排放量均减少3%以上。以降低工程能耗为优化目标,通过对全流程进行换热网络的设计,可使热公用工程能耗降低15%,冷公用工程能耗降低5.6%。
关键词磷石膏    NH3    天然气净化厂    CO2    流程模拟    换热网络    
Mineralization process simulation of CO2 in tail gas of natural gas purification plant by phosphor gypsum combined with ammonia
Lin Ting , Wei Wenyun , Liu Lingling , Yu Hui     
School of Chemical Engineering, Sichuan University, Chengdu 610065, China
Abstract: Using CCUS new technology of phosphor gypsum, ammonia combined mineralizing CO2 in tail gas of natural gas purification plant, and based on estimation of physical parameters and model establishment, the process conditions for making use of the phosphorus gypsum mineralization of CO2 were optimized with process simulation method. It proposed the optimized reaction temperature was 60 ℃, and the reaction time was 100 minutes. Compared the simulation results of optimized process conditions with the existing technology conditions, it was concluded that the exhaust gas of CO2, NH3, phosphor gypsum emissions were reduced by more than 3% compared with the existing process emissions. Aimed at the optimization goal of reducing the energy consumption, the heat exchange network design of whole process can reduce the energy consumption of hot utility by 15% and that of cold utility by 5.6%, respectively.

天然气净化厂排放的尾气中含有大量CO2,目前工业上应用最多的脱除技术是有机胺吸收法[1],但因脱除成本高、增加发电成本等弊端影响了其工业化应用[2]。近年来,谢和平等针对CO2捕获、利用及储存(Carbon Capture, Utilization and Storage,以下简称CCUS)工艺提出了实现CO2减排的CO2捕获及利用(Carbon Capture and Utilization,以下简称CCU)发展新理念[3],提出利用天然矿物或工业废料的CO2矿化技术,将CO2作为一种资源,在低能耗、低成本的条件下矿化CO2,转化联产高附加值化工产品,真正实现CO2的高效利用[4]。工业废料因具有很好的矿化潜力而得到广泛的研究。例如磷石膏,作为磷铵工业中的固体废弃物,目前中国每年产出约5 000×104 t,利用其中的Ca2+进行CO2矿化,每年可消耗CO2约1.25×108 t[5]。因此,应用Aspen Plus软件对该工艺流程进行综合模拟、计算和分析,针对矿化反应的主要工艺单元,基于灵敏度分析和设计规定进行模拟优化,并进行全流程换热网络的设计,通过Aspen Energy Analyzer中不同方案的筛选和优化,确定最佳换热方案,以期为工程应用提供工艺参数选择的依据。

1 工艺说明
1.1 工艺原理

磷石膏(以二水合物为主)矿化CO2制备硫酸铵和碳酸钙的反应为复分解反应,其反应方程式见式(Ⅰ)和式(Ⅱ)。

$ 2{\text{N}}{{\text{H}}_3} + {\text{C}}{{\text{O}}_2} + {{\text{H}}_2}{\text{O}} \to {({\text{N}}{{\text{H}}_4})_2}{\text{C}}{{\text{O}}_3} $ (Ⅰ)
$ {\text{CaS}}{{\text{O}}_4}\cdot2{{\text{H}}_2}{\text{O}} + {({\text{N}}{{\text{H}}_4})_2}{\text{C}}{{\text{O}}_3} \to {({\text{N}}{{\text{H}}_4})_2}{\text{S}}{{\text{O}}_4} + {\text{CaC}}{{\text{O}}_3} + 2{{\text{H}}_2}{\text{O}} $ (Ⅱ)

反应推动力来源于CaSO4·2H2O和CaCO3的溶解度差别,二者的溶度积常数相差近3 000倍,CaSO4·2H2O溶解的Ca2+与NH4HCO3产生的CO32-结合,形成更难溶的CaCO3沉淀,使溶液中Ca2+浓度降低,促使CaSO4·2H2O溶解。同时,NH3·H2O的存在,使体系呈碱性,保持足够量的CO32-,使CaSO4·2H2O较完全地转化为CaCO3结晶。整个反应过程的速度与CaSO4·2H2O溶解速度、Ca2+和CO32-扩散速度及形成CaCO3晶体的速度等有关[6]

1.2 基本条件
1.2.1 原料条件

(1) 以中国石化普光天然气净化厂尾气为例,其尾气组成见表 1

表 1    尾气组成 Table 1    Composition of tail gas

(2) 以瓮福集团产生的磷石膏固体废物为例,其组成(干基)见表 2

表 2    磷石膏组成(干基) Table 2    Composition of phosphor gypsum (dry basis)

1.2.2 公用物料规格

该工艺中公用物料规格见表 3

表 3    公用物料规格 Table 3    Specifications of common materials

1.3 拟实现的性能指标

本工艺拟实现尾气CO2矿化率≥75%、磷石膏转化率≥90%、NH3利用率≥99%的指标。

2 流程模拟
2.1 数学模型

低浓度尾气CO2直接转化磷石膏制硫酸铵的总反应见式(Ⅲ)。

$ {\text{C}}{{\text{O}}_2} + 2{\text{N}}{{\text{H}}_3} + {\text{CaS}}{{\text{O}}_4}\cdot2{{\text{H}}_2}{\text{O }} \to {({\text{N}}{{\text{H}}_4})_2}{\text{S}}{{\text{O}}_4} + {\text{CaC}}{{\text{O}}_3} + {{\text{H}}_2}{\text{O}} $ (Ⅲ)

该反应是一个复杂的气-液-固三相反应过程,其间可能发生的化学反应主要有6种[7]。一般认为主要遵循以下连串、可逆反应步骤[8]

$ {\text{C}}{{\text{O}}_2} + 2{\text{N}}{{\text{H}}_3} \rightleftharpoons {\text{N}}{{\text{H}}_2}{\text{CO}}{{\text{O}}^ - } + {\text{NH}}_4^{\text{ + }} + {{\text{Q}}_1} $ (Ⅳ)
$ {\text{N}}{{\text{H}}_2}{\text{CO}}{{\text{O}}^ - } + {{\text{H}}_{\text{2}}}{\text{O}} \rightleftharpoons {\text{NH}}_4^{\text{ + }} + {\text{CO}}_3^{2 - } - {{\text{Q}}_2} $ (Ⅴ)
$ {\text{C}}{{\text{a}}^{2 + }} + {\text{CO}}_3^{2 - } \rightleftharpoons {\text{CaC}}{{\text{O}}_3} + {{\text{Q}}_3} $ (Ⅵ)

当反应达到平衡时,溶液中的[Ca2+]为定值。

$ [{\rm{C}}{{\rm{a}}^{2 + }}] = \frac{{{K_{{\rm{SP, }}}}_{{\rm{CaS}}{{\rm{O}}_{\rm{4}}} \cdot {\rm{2}}{{\rm{H}}_{\rm{2}}}{\rm{O}}}}}{{{\rm{[SO}}_4^{2-}]}} = \frac{{{K_{{\rm{SP, CaC}}{{\rm{O}}_{\rm{3}}}}}}}{{[{\rm{CO}}_3^{2{\rm{-}}}]}} $ (1)

在25℃下,

$ {K_{{\rm{SP, CaS}}{{\rm{O}}_{\rm{4}}} \cdot {\rm{2}}{{\rm{H}}_{\rm{2}}}{\rm{O}}}} = [{\rm{C}}{{\rm{a}}^{2 + }}]{\rm{[SO}}_4^{2-}] = 9.1 \times {10^{ - 6}} $ (2)
$ {K_{{\rm{SP, CaC}}{{\rm{O}}_{\rm{3}}}}} = [{\rm{C}}{{\rm{a}}^{2 + }}][{\rm{CO}}_3^{2{\rm{-}}}] = 2.9 \times {10^{ - 6}} $ (3)

由于CaSO4·2H2O的溶度积比CaCO3的溶度积要大得多,反应主要向右侧进行,故磷石膏转化率一般较大[9]。CaSO4·2H2O在碳铵溶液中的转化反应可分为液相中CO32-扩散至固液界面处,与溶解于界面内的Ca2+结合生成CaCO3微晶析出,SO42-再扩散至溶液中等步骤,反应为阴离子转移及晶体重组过程。反应生成CaCO3微晶聚集于未反应的CaSO4·2H2O固相周围,形成包裹硫酸钙的产物层。转化反应过程可以用缩芯模型来描述,属于化学反应控制。其宏观动力学反应可用以下反应速率方程表达[10]

$ r = {k_1}[{\rm{CO}}_3^{2{\rm{-}}}] $ (4)
$ {k_1} = 5.99 \times {10^{10}}{\rm{exp}}(\frac{{ - 68\;920}}{{RT}}) $ (5)

式中,r为反应速率,mol/(m3·s);k1为速率常数,mol/(m3·s);T为温度,K;R为通用气体常数,8.314 J/mol·K。

由于工艺流程涉及无机电解质体系,故选用ENRTL-RK物性方法对物性进行关联。

2.2 工艺流程

磷石膏制硫酸铵的现有生产工艺[11]过程主要由磷石膏洗涤工段、尾气净化工段、矿化反应体系、固液分离工段、硫酸铵结晶、干燥工段等组成。针对现有生产工艺流程进行必要的简化后,建立硫酸铵生产的工艺流程模型,如图 1所示。

图 1     硫酸铵生产的工艺流程图 Figure 1     Process flow of ammonium sulfate production

2.3 现有工艺操作条件

该工艺中反应流程以NH3为动力源、热源和反应促进剂,在全混流反应器中进行反应,替代传统的液固搅拌和热源,缩短流程、减少设备、降低能耗。反应流程的三相反应器应选用全混流反应器CSTR(Continuous Stirred Tank Reactor)。目前,朱家骅[12]等对磷石膏矿化CO2制备硫酸铵和碳酸钙工艺流程进行了详细研究,并得到气-液-固三相全混流反应器中的操作条件,见表 4

表 4    现有工艺操作条件 Table 4    Operating conditions of existing process

该反应的有效相态为气-液两相,停留时间3 h。

3 模拟条件的优化
3.1 灵敏度分析

温度是影响反应速率的主要因素,温度越高,反应速率越快。但在该反应体系中,温度升高将导致碳铵溶液的分解[13]。在全混流反应器中,对碳铵与磷石膏反应的反应物浓度与生成物浓度在温度为0~150 ℃范围内做灵敏度分析,结果如图 2所示。

图 2     反应在0~150 ℃范围内的灵敏度分析图 Figure 2     Reaction sensitivity analysis at the range of 0-150 ℃

图 2可知,生成的硫酸铵和碳酸钙摩尔流量随温度的升高而增加,但当温度达到60 ℃时,由于碳铵的分解使得反应生成物的摩尔流量逐渐减少。因此,磷石膏与碳铵反应温度若控制在60 ℃,硫酸铵和碳酸钙的产率将达最大。

3.2 停留时间的确定

通过模拟计算得到碳铵反应转化率随停留时间变化的关系,如图 3所示。由图 3可知,随着停留时间的增加,磷石膏的转化率不断增大。当停留时间达到100 min时,转化率基本不变,即反应时间在100 min时就可以达到较高的转化率。通过上述模拟计算分析优化,得到在全混流反应器中的最佳反应条件:温度为60 ℃,停留时间为100 min。

图 3     停留时间对碳铵反应转化率的影响 Figure 3     Influence of retention time on reaction conversion rate of ammonium carbonate

3.3 模拟结果及分析

将上述最佳反应条件带入Aspen Plus流程模拟中进行再次模拟分析,将得到的数据结果与现有工艺操作条件下的数据进行对比,见表 5

表 5    优化前后工艺模拟对比数据表 Table 5    Comparison of process simulation data before and after optimization

表 5所示,反应条件优化之后,转化率仍可达到性能指标。与现有的操作工艺条件相比,不仅降低了反应温度,且尾气、NH3、磷石膏排放量分别减少3.5%、3.7%和4.1%。

4 换热网络的设计

将流程模拟得到所需的物流数据导入Aspen Energy Analyzer,形成反应工段模拟物流信息,见表 6

表 6    反应工段流程模拟物流信息 Table 6    Stream information of reaction process simulation

在换热过程中,只有在热物流的温度超过冷物流时,才能将热量从热物流自发地传递至冷物流。为使设计可行,换热器内部进行换热的每个点热、冷流之间均要存在一定的温差,根据物流信息,取最小传热温差为10 ℃,计算得到冷、热物流夹点图,如图 4所示。

图 4     物流夹点图 Figure 4     Pinch point of the stream

图 4可知,当热物流温度为54.7 ℃,冷物流温度为44.7 ℃时,换热温差最小,换热网络具有最大的热回收量,且需要的公用工程量也最小。通过Aspen Energy Analyzer可确定出所需的最小公用工程量,如表 7所示。

表 7    夹点温度和公用工程 Table 7    Pinch point temperature and utilities

在确定夹点温度和公用工程后,根据夹点原理,通过物流之间的换热匹配,实现最大程度的热量集成。经过Aspen Energy Analyzer中不同换热方案的筛选和优化,可确定以下最优换热方案。

含有CO2的的尾气温度为300 ℃,先用于加热磷石膏料浆,使尾气出口温度降至132.4 ℃,再用于加热NH3,使NH3达到指定换热效果,尾气也降至所需的85 ℃,磷石膏料浆升温至45.26 ℃,再以低压蒸汽加热至所需温度。富氨溶液和碳铵溶液均由公用工程提供热量。

经过换热物流匹配后,得到换热网络的公用工程信息表,见表 8

表 8    流程公用工程信息 Table 8    Utility information of the process

5 结论

(1) 通过对磷石膏、NH3联合矿化工业尾气中CO2工艺全流程的模拟与反应条件的优化,得到最佳反应条件:温度60 ℃;停留时间100 min。最后通过全流程模拟得到尾气转化率88.4%,磷石膏转化率95.9%,NH3转化率99.0%,均符合性能要求。与现有工艺操作条件相比,尾气排放量减少3.5%,NH3排放量减少3.7%,磷石膏排放量减少4.1%。

(2) 通过利用Aspen Energy Analyzer软件进行换热网络的设计,实现了对冷热流股之间的匹配,优化了换热方案,热公用工程耗量可节约15%,冷公用工程耗量可节约5.6%。

参考文献
[1]
张永军, 苑慧敏, 万书宝, 等. 天然气中二氧化碳脱除技术[J]. 化工中间体, 2008(9): 2.
[2]
CIFERNO J P, DIPIETRO P, TARKA T. An economic scoping study for CO2 capture using aqueous ammonia [R]. Washington D. C: US Department of Energy DOC, 2005.
[3]
谢和平, 谢凌志, 王昱飞, 等. 全球二氧化碳减排不应是CCS, 应是CCU[J]. 四川大学学报:工程科学版, 2012, 44(4): 1-5.
[4]
TAI C Y, CHEN W R, SHIH S M. Factors affecting wollastonite carbonation under CO2 supercritical conditions[J]. AIChE Journal, 2006, 52(1): 292-299. DOI:10.1002/(ISSN)1547-5905
[5]
谢和平, 刘虹, 吴刚. 中国未来二氧化碳减排技术应向CCU方向发展[J]. 中国能源, 2012, 34(10): 15-18. DOI:10.3969/j.issn.1003-2355.2012.10.003
[6]
张茂林, 王永秋, 刘清理, 等. 磷石膏转化法生产硫酸钾的工艺研究[J]. 淮北煤师院学报:自然科学版, 2002, 23(2): 32-36.
[7]
VUČAK M, PERI J, MIKI A, et al. A study of carbon dioxide absorption into aqueous monoethanolamine solution containing calcium nitrate in the gas-liquid reactive precipitation of calcium carbonate[J]. Chemical Engineering Journal, 2002, 87(2): 171-179. DOI:10.1016/S1385-8947(01)00239-X
[8]
李希, 戈育芳. 氨盐水碳酸化反应的机理与反应器模型Ⅰ碳酸化反应的动力学模型[J]. 化学反应工程与工艺, 1999, 12(4): 350-356.
[9]
万邦江, 胡武洪. 磷石膏预处理及综合利用[J]. 中国西部科技, 2011, 10(25): 5-7. DOI:10.3969/j.issn.1671-6396.2011.25.003
[10]
鲁厚芳. 杂质对石膏-碳酸铵转化过程的影响[D]. 成都: 四川大学, 2002.
[11]
朱家骅, 谢和平, 夏素兰, 等. 烟气二氧化碳与磷石膏转化一步法节能节水清洁工艺: 中国, 201210223218. 4[P]. 2012-07-02.
[12]
朱家骅, 谢和平, 何浩明, 等. (四川大学)三相全混流二氧化碳氨化矿化反应系统与装置: 中国, 201210237185. 9[P]. 2012-07-10.
[13]
刘芳, 王淑娟, 张曦, 等. 碳酸氢铵溶液解吸二氧化碳研究[J]. 环境科学学报, 2009, 29(9): 1886-1890.