湿法烟气脱硫塔内传递与化学反应过程CFD模拟

doi:10.11949/j.issn.0438-1157.20181368

化工学报 ›› 2019, Vol. 70 ›› Issue (6): 2117-2128.DOI: 10.11949/j.issn.0438-1157.20181368

湿法烟气脱硫塔内传递与化学反应过程CFD模拟

曲江源¹(),齐娜娜¹(),关彦军¹,滕阳¹,徐文青²,朱廷钰²,张锴¹

1. 华北电力大学热电生产过程污染物监测与控制北京市重点实验室，北京 102206
2. 中国科学院过程工程研究所湿法冶金清洁生产技术国家工程实验室，北京 100190

收稿日期:2018-11-18 修回日期:2019-03-10 出版日期:2019-06-05 发布日期:2019-06-05
通讯作者: 齐娜娜
作者简介:<named-content content-type="corresp-name">曲江源</named-content>（1993－），男，博士研究生，<email>qujy_ncepu@163.com</email>
基金资助:
国家重点研发计划项目(2017YFC0210601);青年科学基金项目(51706070)

CFD simulation of transfer and chemical reaction process in wet flue gas desulfurization tower

Jiangyuan QU¹(),Nana QI¹(),Yanjun GUAN¹,Yang TENG¹,Wenqing XU²,Tingyu ZHU²,Kai ZHANG¹

1. Beijing Key Laboratory of Emission Surveillance and Control for Thermal Power Generation, North China Electric Power University, Beijing 102206, China
2. National Engineering Laboratory for Hydrometallurgical Cleaner Production Technology, Institute of Process Engineering, Chinese Academy of Sciences, Beijing 100190, China

Received:2018-11-18 Revised:2019-03-10 Online:2019-06-05 Published:2019-06-05
Contact: Nana QI

摘要/Abstract

摘要：

以某330 MW燃煤发电机组湿法烟气脱硫塔为研究对象，采用欧拉-拉格朗日方法建立了塔内气-液两相流动、热质交换模型，以溶解平衡、质量守恒及电荷守恒描述浆液内13种溶质组分瞬时化学反应特性，通过用户自定义函数实现流动模型与传质模型、化学反应模型的耦合。基于上述模型预测了脱硫塔内气-液两相流动、液滴蒸发与SO₂化学吸收过程，获得了SO₂浓度与浆液pH的径向分布特性，详细分析了气-液流动对化学吸收特性的影响，结果表明局部液气比分布特性是影响SO₂径向分布的关键因素之一，可通过调控近壁区及主管道区的两相流动状态提高脱硫塔的吸收性能；随着气相侧SO₂浓度提高或液滴粒径减小，浆液pH下降速率增大且各化学组分浓度达到稳定状态用时缩短。

关键词: 湿法烟气脱硫, 气液两相流, 热质交换, 吸收, 数值模拟

Abstract:

The Eulerian-Lagrangian method is employed to describe the gas-liquid flow and the heat and mass transfer process in a full-scale wet flue gas desulfurization(WFGD) tower, which is in operation at a 330 MW coal-fired power unit. Instantaneous chemical reactions involving 13 dissolved species in aqueous phase are assumed to be at equilibrium obeying mass balance for the specific chemical compositions and electro neutrality. The gas-liquid flow model is coupled with mass transfer and chemical reactions procedures through user defined functions (UDFs). The gas-liquid hydrodynamics, droplets evaporation and SO₂ chemical absorption is predicted in the spray scrubber. Radial distribution characteristics of SO₂ concentration and pH of the droplets are obtained and the effect of gas-liquid flow behavior on chemical absorption process is analyzed. The results show that the local liquid-gas ratio distribution is one of the key factors affecting the radial distribution of SO₂. Furthermore, the absorption performance of the desulfurization tower can be improved by controlling the two-phase flow pattern in the area around the wall and that beneath the main pipe. The drop rates of pH increase with the raise of SO₂ concentration in flue gas or the decrease of droplet diameter. Besides, the higher concentration of SO₂ in gas phase or the smaller characteristic size of droplet can accelerate the process that the concentration of each chemical component reaches a stable state.

Key words: wet flue gas desulfurization, gas-liquid flow, heat and mass transfer, absorption, numerical simulation

中图分类号:

X 701.3

曲江源, 齐娜娜, 关彦军, 滕阳, 徐文青, 朱廷钰, 张锴. 湿法烟气脱硫塔内传递与化学反应过程CFD模拟[J]. 化工学报, 2019, 70(6): 2117-2128.

Jiangyuan QU, Nana QI, Yanjun GUAN, Yang TENG, Wenqing XU, Tingyu ZHU, Kai ZHANG. CFD simulation of transfer and chemical reaction process in wet flue gas desulfurization tower[J]. CIESC Journal, 2019, 70(6): 2117-2128.

图/表 12

图1 脱硫塔及计算域

Fig.1 Schematic diagram of desulfurization tower and computational domain

表1 液相内化学反应平衡常数

Table 1 Reaction equilibrium constants in droplets

Parameter	Formula	Eq.
$K H 2 O$ /(mol/L)²	[H⁺][OH^-]	(19)
$K S O 2, a q$ /(mol/L)	[H⁺][HSO₃ ^-]/[SO_2,aq]	(20)
$K H S O 3 -$ /(mol/L)	[H⁺][SO₃ ^2-]/[HSO₃ ^-]	(21)
$K H S O 4 -$ /(mol/L)	[H⁺][SO₄ ^2-]/[HSO₄ ^-]	(22)
$K C O 2, a q$ /(mol/L)	[H⁺][HCO₃ ^-]/[CO_2,aq]	(23)
$K H C O 3 -$ /(mol/L)	[H⁺][CO₃ ^2-]/[HCO₃ ^-]	(24)

表1 液相内化学反应平衡常数

Table 1 Reaction equilibrium constants in droplets

Parameter	Formula	Eq.
$K H 2 O$ /(mol/L)²	[H⁺][OH^-]	(19)
$K S O 2, a q$ /(mol/L)	[H⁺][HSO₃ ^-]/[SO_2,aq]	(20)
$K H S O 3 -$ /(mol/L)	[H⁺][SO₃ ^2-]/[HSO₃ ^-]	(21)
$K H S O 4 -$ /(mol/L)	[H⁺][SO₄ ^2-]/[HSO₄ ^-]	(22)
$K C O 2, a q$ /(mol/L)	[H⁺][HCO₃ ^-]/[CO_2,aq]	(23)
$K H C O 3 -$ /(mol/L)	[H⁺][CO₃ ^2-]/[HCO₃ ^-]	(24)

图2 网格划分与不同高度XY截面SO2平均浓度相对偏差分布

Fig.2 Schematic diagram of geometry discretization and relative errors profiles of average SO2 concentration at different XY planes

图3 烟气速度分布云图

Fig.3 Contour of gas velocity distribution

图4 无量纲轴向速度径向分布

Fig.4 Dimensionless axial velocity radial profiles

图5 浆液浓度分布云图

Fig.5 Contour of discrete phase concentration

图6 不同截面液滴粒径分布特性

Fig.6 Diameter distribution character of droplets in different cross sections

表2 温度、水蒸气体积分数实测值与计算值

Table 2 Measured values and simulation results of gas temperature and vapor volume fraction

项目	T _g,in /℃	T _g,out /℃	$V H 2 O, i n$ /%	$V H 2 O, o u t / %$
实测 Ⅰ	125.0	48.0	7.5	13.2
实测 Ⅱ	135.0	50.0	6.2	13.7
实测 Ⅲ	131.0	51.0	7.1	13.4
平均值	130.3	49.7	7.0	13.4
模拟值	130.0	50.0	7.0	13.8
相对误差	—	0.60%	—	2.99%

表2 温度、水蒸气体积分数实测值与计算值

Table 2 Measured values and simulation results of gas temperature and vapor volume fraction

项目	T _g,in /℃	T _g,out /℃	$V H 2 O, i n$ /%	$V H 2 O, o u t / %$
实测 Ⅰ	125.0	48.0	7.5	13.2
实测 Ⅱ	135.0	50.0	6.2	13.7
实测 Ⅲ	131.0	51.0	7.1	13.4
平均值	130.3	49.7	7.0	13.4
模拟值	130.0	50.0	7.0	13.8
相对误差	—	0.60%	—	2.99%

图7 入口烟气温度对塔内烟气温度与相对湿度的影响

Fig.7 Influence of initial gas temperature on temperature and relative humidity of flue gas in tower

图8 烟气中SO2浓度径向分布

Fig.8 Radial profiles of SO2 concentration in gas phase

图9 喷淋浆液空间分布与pH分布

Fig.9 Slurry distribution in spray scrubber colored by pH

图10 喷嘴位置与液滴粒径对pH的影响

Fig.10 Influences of injection location and droplet diameter on pH

参考文献 34

1	史文峥, 杨萌萌, 张绪辉, 等 . 燃煤电厂超低排放技术路线与协同脱除[J]. 中国电机工程学报, 2016, 36(16): 4308-4318.
	Shi W Z , Yang M M , Zhang X H , et al . Ultra-low emission technical route of coal-fired power plants and the cooperative removal[J]. Proceedings of the CSEE, 2016, 36(16): 4308-4318.
2	潘丹萍, 郭彦鹏, 黄荣廷, 等 . 石灰石-石膏法烟气脱硫过程中细颗粒物形成特性[J]. 化工学报, 2015, 66(11): 4618-4625.
	Pan D P , Guo Y P , Huang R T , et al . Formation of fine particles in flue gas desulphurization process using limestone-gypsum[J]. CIESC Journal, 2015, 66(11): 4618-4625.
3	张东辉, 庄烨, 朱润儒, 等 . 燃煤烟气污染物超低排放技术及经济分析[J]. 电力建设, 2015, 36(5): 125-130.
	Zhang D H , Zhuang Y , Zhu R R , et al . Ultra-low air pollutant control technologies for coal-fired flue gas and its economic analysis[J]. Electric Power Construction, 2015, 36(5): 125-130.
4	张军 . 超低排放的湿法高效脱硫协同除尘的机理及模型研究[D]. 杭州: 浙江大学, 2018.
	Zhang J . Mechanism and model study on high efficiency desulfurization and dedust based on ultra-low emission technology [D]. Hangzhou: Zhejiang University, 2018.
5	张军, 郑成航, 张涌新, 等 . 某1000MW燃煤机组超低排放电厂烟气污染物排放测试及其特性分析[J]. 中国电机工程学报, 2016, 36(5): 1310-1314.
	Zhang J , Zheng C H , Zhang Y X , et al . Experimental investigation of ultra-low pollutants emission characteristics from a 1000MW coal-fired power plant[J]. Proceedings of the CSEE, 2016, 36(5): 1310-1314.
6	杨业, 徐超群, 朱燕群, 等 . 臭氧氧化结合硫代硫酸钠溶液喷淋同时脱硫脱硝[J]. 化工学报, 2016, 67(5): 2041-2047.
	Yang Y , Xu C Q , Zhu Y Q , et al . Simultaneous removal of SO₂ and NO _x by combination of ozone oxidation and Na₂S₂O₃ solution spray [J]. CIESC Journal, 2016, 67(5): 2041-2047.
7	Córdoba P . Status of flue gas desulphurization (FGD) systems from coal-fired power plants: overview of the physic-chemical control processes of wet limestone FGDs [J]. Fuel, 2015, 144: 274-286.
8	Angelo J , Lightfoot E , Howard D . Generalization of the penetration theory for surface stretch: application to forming and oscillating drops[J]. AIChE Journal, 1966, 12(4): 751-760.
9	Altwicker E , Lindhjem C . Absorption of gases into drops[J]. AIChE Journal, 1988, 34(34): 329-332.
10	Hsu C , Shih S . Semiempirical equation for liquid-phase mass-transfer coefficient for drops[J]. AIChE Journal, 1993, 39(6): 1090-1092.
11	Brogren C , Karlsson H . Modeling the absorption of SO₂, in a spray scrubber using the penetration theory[J]. Chemical Engineering Science, 1997, 52(18): 3085-3099.
12	Zhong Y , Gao X , Huo W , et al . A model for performance optimization of wet flue gas desulfurization systems of power plants[J]. Fuel Processing Technology, 2008, 89(11): 1025-1032.
13	祝杰, 叶世超, 吴振元, 等 . 石灰石-石膏脱硫喷淋塔吸收区高度计算模型[J]. 化工学报, 2014, 65(8): 2896-2901.
	Zhu J , Ye S C , Wu Z Y , et al . Mathematical modeling for height of absorption region in spray tower for limestone-gypsum wet flue gas desulfurization[J]. CIESC Journal, 2014, 65 (8): 2896-2901.
14	Nygaard H , Kiil S , Johnsson J , et al . Full-scale measurements of SO₂ gas phase concentrations and slurry compositions in a wet flue gas desulphurization spray absorber[J]. Fuel, 2004, 83(9): 1151-1164.
15	Gao X , Huo W , Luo Z Y , et al . CFD simulation with enhancement factor of sulfur dioxide absorption in the spray scrubber[J]. Journal of Zhejiang University-Science A(Applied Physics & Engineering), 2008, 9(11): 1601-1613.
16	Gómez A , Fueyo N , Tomás A . Detailed modeling of a flue-gas desulfurization plant[J]. Computers & Chemical Engineering, 2007, 31(11): 1419-1431.
17	Crowe C , Sharma M , Stock D . The particle-source-in cell (PSI-CELL) model for gas-droplet flow[J]. Journal of Fluids Engineering, 1977, 99(2): 325-332.
18	林永明, 高翔, 施平平, 等 . 大型湿法烟气脱硫喷淋塔内阻力特性数值模拟[J]. 中国电机工程学报, 2008, 28(5): 28-33.
	Lin Y M , Gao X , Shi P P , et al . Numerical simulation on resistance characteristic of large scale wet flue gas desulfurization spraying scrubber[J]. Proceedings of the CSEE, 2008, 28(5): 28-33.
19	钟毅, 高翔, 王惠挺, 等 . 基于CFD技术的湿法烟气脱硫系统性能优化[J]. 中国电机工程学报, 2008, 28(32): 18-23.
	Zhong Y , Gao X , Wang H T , et al . Performance optimization of wet flue gas desulphurization system based on CFD technology[J]. Proceedings of the CSEE, 2008, 28(32): 18-23.
20	Chen Z , Wang H M , Zhuo J K , et al . Experimental and numerical study on effects of deflectors on flow field distribution and desulfurization efficiency in spray towers[J]. Fuel Processing Technology, 2017, 162: 1-12.
21	Marocco L , Inzoli F . Multiphase Euler-Lagrange CFD simulation applied to wet flue gas desulphurization technology[J]. International Journal of Multiphase Flow, 2009, 35(2): 185-194.
22	Marocco L . Modeling of the fluid dynamics and SO₂ absorption in a gas-liquid reactor[J]. Chemical Engineering Journal, 2010, 162(1): 217-226.
23	Montanes C , Gomez S A , Fueyo N , et al . Computational evaluation of wall rings in wet flue-gas desulfurization plants[J]. International Journal of Energy for A Clean Environment, 2009, 10(1/2/3/4): 15-36.
24	Mors S , Alexander A . An investigation of particle trajectories in two-phase flow systems[J]. Journal of Fluid Mechanics, 1972, 55(2): 193-208.
25	Wang Y , James P . Assessment of an eddy-interaction model and its refinements using predictions of droplet deposition in a wave-plate demister[J]. Chemical Engineering Research & Design, 1999, 77(8): 692-698.
26	林瑜, 陈德珍, 尹丽洁 . 喷淋层组合方式对大型脱硫塔内流动和热湿交换过程影响的数值模拟[J]. 中南大学学报(自然科学版), 2017, 48(10): 2572-2582.
	Lin Y , Chen D Z , Yin L J . Numerical simulation of impact of spraying layers scheme on gas-liquid two-phases flow, heat and mass transfer in large scale desulphurization absorption tower[J]. Journal of Central South University (Science and Technology), 2017, 48(10): 2572-2582.
27	Zhu J , Ye S C , Bai J , et al . A concise algorithm for calculating absorption height in spray tower for wet limestone-gypsum flue gas desulfurization[J]. Fuel Processing Technology, 2015, 129: 15-23.
28	Füller E N , Schettler P D , Giddings J C . New method for prediction of binary gas-phase diffusion coefficient[J]. Industrial & Engineering Chemistry, 1966, 58(5): 19-27.
29	Newton G H , Kramlich J , Payne R . Modeling the SO₂-slurry droplet reaction.[J]. AIChE Journal, 1990, 36(12): 1865-1872.
30	Maurer G . On the solubility of volatile weak electrolytes in aqueous solutions—thermodynamics of aqueous systems with industrial applications[J]. ACS Symposium, 1980, 133: 139-172.
31	Neveux T , Moullec Y L . Wet industrial flue gas desulfurization unit: model development and validation on industrial data[J]. Industrial & Engineering Chemistry Research, 2011, 50(12): 7579-7592.
32	Wang P H , Zhuang L W , Dai G C . Synergistic effect of droplet self-adjustment and rod bank internal on fluid distribution in a WFGD spray column[J]. Chemical Engineering Science, 2017, 162: 227-244.
33	Tseng C C , Li C J . Eulerian-Eulerian numerical simulation for a flue gas desulfurization tower with perforated sieve trays[J]. International Journal of Heat and Mass Transfer, 2018, 116: 329-345.
34	Chen Z , You C F , Liu H Z , et al . The synergetic particles collection in three different wet flue gas desulfurization towers: a pilot-scale experimental investigation[J]. Fuel Processing Technology, 2018, 179: 344-350.

[1]	宋嘉豪, 王文. 斯特林发动机与高温热管耦合运行特性研究[J]. 化工学报, 2023, 74(S1): 287-294.
[2]	张思雨, 殷勇高, 贾鹏琦, 叶威. 双U型地埋管群跨季节蓄热特性研究[J]. 化工学报, 2023, 74(S1): 295-301.
[3]	黄琮琪, 吴一梅, 陈建业, 邵双全. 碱性电解水制氢装置热管理系统仿真研究[J]. 化工学报, 2023, 74(S1): 320-328.
[4]	金正浩, 封立杰, 李舒宏. 氨水溶液交叉型再吸收式热泵的能量及分析[J]. 化工学报, 2023, 74(S1): 53-63.
[5]	叶展羽, 山訸, 徐震原. 用于太阳能蒸发的折纸式蒸发器性能仿真[J]. 化工学报, 2023, 74(S1): 132-140.
[6]	张义飞, 刘舫辰, 张双星, 杜文静. 超临界二氧化碳用印刷电路板式换热器性能分析[J]. 化工学报, 2023, 74(S1): 183-190.
[7]	王志国, 薛孟, 董芋双, 张田震, 秦晓凯, 韩强. 基于裂隙粗糙性表征方法的地热岩体热流耦合数值模拟与分析[J]. 化工学报, 2023, 74(S1): 223-234.
[8]	江河, 袁俊飞, 王林, 邢谷雨. 均流腔结构对微细通道内相变流动特性影响的实验研究[J]. 化工学报, 2023, 74(S1): 235-244.
[9]	米泽豪, 花儿. 基于DFT和COSMO-RS理论研究多元胺型离子液体吸收SO₂气体[J]. 化工学报, 2023, 74(9): 3681-3696.
[10]	袁佳琦, 刘政, 黄锐, 张乐福, 贺登辉. 泡状入流条件下旋流泵能量转换特性研究[J]. 化工学报, 2023, 74(9): 3807-3820.
[11]	何松, 刘乔迈, 谢广烁, 王斯民, 肖娟. 高浓度水煤浆管道气膜减阻两相流模拟及代理辅助优化[J]. 化工学报, 2023, 74(9): 3766-3774.
[12]	邢雷, 苗春雨, 蒋明虎, 赵立新, 李新亚. 井下微型气液旋流分离器优化设计与性能分析[J]. 化工学报, 2023, 74(8): 3394-3406.
[13]	韩晨, 司徒友珉, 朱斌, 许建良, 郭晓镭, 刘海峰. 协同处理废液的多喷嘴粉煤气化炉内反应流动研究[J]. 化工学报, 2023, 74(8): 3266-3278.
[14]	高燕, 伍鹏, 尚超, 胡泽君, 陈晓东. 基于双流体喷嘴的磁性琼脂糖微球的制备及其蛋白吸附性能探究[J]. 化工学报, 2023, 74(8): 3457-3471.
[15]	程小松, 殷勇高, 车春文. 不同工质在溶液除湿真空再生系统中的性能对比[J]. 化工学报, 2023, 74(8): 3494-3501.

湿法烟气脱硫塔内传递与化学反应过程CFD模拟

CFD simulation of transfer and chemical reaction process in wet flue gas desulfurization tower

RichHTML

PDF (PC)

可视化

被引次数

摘要/Abstract

引用本文

使用本文

图/表 12

参考文献 34

相关文章 15

编辑推荐

Metrics

本文评价