一种氨法脱硫喷淋塔系统及SO₂吸收传质系数的确定方法

一种氨法脱硫喷淋塔系统及SO₂吸收传质系数的确定方法
【专利摘要】本发明公开了一种氨法脱硫喷淋塔系统及SO2吸收传质系数的确定方法，属于氨法烟气脱硫技术领域。本发明基于喷淋塔氨法烟气脱硫工艺的特点，选取喷淋塔内浆液pH值、烟气流速ug和液气比L/Q为参数，先计算出不同pH值、烟气流速ug及液气比L/Q条件下的氨法脱硫SO2吸收传质速率ky，然后采用数值拟合的方法重新构建出具有普遍适用性的氨法脱硫SO2吸收传质系数ky′的表达式。采用本发明的方法得到的氨法脱硫SO2传质系数能够全面、准确的反映工艺操作条件对SO2传质吸收过程的影响，且具有普遍适用性，能够为优化氨法脱硫吸收塔的设计和运行提供理论指导，具有较强的实际意义和工程应用价值，且其计算结果与实验结果的误差小于14.77％。
【专利说明】
一种氨法脱硫喷淋塔系统及S〇2吸收传质系数的确定方法
技术领域
[0001 ]本发明属于氨法烟气脱硫技术领域，更具体地说，涉及一种氨法脱硫喷淋塔系统 及S〇2吸收传质系数的确定方法。
【背景技术】
[0002] 二氧化硫(S02)是我国目前危害仅次于可吸入颗粒物的主要大气污染物，但它同 时也是一种重要资源，目前我国脱硫方法主要采用石灰石-石膏法脱硫工艺，一是因为脱硫 剂石灰石来源丰富且价格便宜，二是因为石灰石-石膏法技术在国外相当成熟且公开，获取 容易。但是该方法所用系统较复杂、投资多、占地面积大，且容易产生二次污染、运行费用 高，中小热电厂难以承受，从而限制了该法的应用。
[0003] 氨法脱硫是以化学反应活性较高的氨(NH3)为脱硫剂吸收烟气中的S02,系统动力 消耗低，在高效脱除so 2的同时，还能实现硫资源的回收，得到硫酸铵副产品（该过程涉及的 主要反应如下（1)~（5)式），且无废水和废渣的排放，符合我国发展循环经济、创建节约型 社会的政策需求，成为缓解上述石灰石-石膏法脱硫工艺存在问题的有效途径，该技术目前 在我国大中型燃煤烟气脱硫项目上也逐步得到了广泛应用，发展前景良好。
[0004] 2NH3+H20+S〇2^ (NH4) 2SO3 (1)
[0005] (NH4)2S〇3+S〇2+H2〇-2NH4HS〇3 (2)
[0006] NH4HS〇3+NH3-(NH4)2S〇3 (3)
[0007] (NH4) 2 S03+1 /2〇2^ (NH4) 2 S04 (4)
[0008] NH4HS〇3+l/2〇2^NH4HS〇4 (5)
[0009] 吸收塔是氨法脱硫工艺的核心设备，工程上普遍采用的是喷淋吸收塔。目前，我国 氨法烟气脱硫喷淋塔的设计和运行主要依据流场模拟对气、液流速分布进行优化，以及根 据经验控制pH值和停留时间等参数来实现，系统在运行过程中很难达到技术性和经济性的 最优组合，且容易发生明显的氨的挥发性损耗和硫酸铵气溶胶逃逸等问题。另外，2012年1 月1日开始实施的《火电厂大气污染物排放标准》(GB13223-2011)进一步收紧了 S02和颗粒 物的排放浓度限值，这就对现有氨法脱硫技术提出了更高的要求。研究并确立氨法烟气脱 硫S0 2吸收传质系数的表达式，是对氨法脱硫喷淋塔工程的设计和运行进行有效优化的关 键，对实现该"绿色"脱硫技术的尚效、稳定运彳丁，提尚氣法脱硫效率以及抑制氣的挥发性损 耗和气溶胶粒子的逃逸具有重要的意义。
[0010] 目前，有关氨法脱硫s〇2吸收传质系数方面的研究尚未见报道，国内外学者通过对 与之相关的S02吸收反应动力学和吸收过程建模进行了一些研究，比较典型的有:Gao X等 (Applied Energy，2010,87(8):2647~2651)依据氨法脱硫实质上是以（NH4)2S03来吸收S〇2 这一特点，对(NH 4) 2S03吸收S02的反应动力学过程进行了研究，结果表明当(NH 4) 2S03浓度在 0.05mol/L以下时，吸收过程由气膜和液膜共同控制，当（NH4) 2S03浓度在0.05mol/L以上时 由气膜控制，且吸收反应速率与(NH4) 2S〇3浓度无关，与初始S02浓度呈0.6级响应。陈梅倩等 对氨法脱硫化学反应动力学进行了实验研究，并建立了氨法脱硫反应的Arrhenius方程(环 境科学学报，2005,25(7):886~890)。刘国荣（化工学报，2010,61(9):2463~2467)和彭健 (高校化学工程学报，2011，25(6): 1073~1077)等基于气相传质系数和液相传质系数相等 这一原理，建立数学模型分别对喷淋塔和湿壁塔氨法脱硫过程进行了数值模拟。上述氨法 脱硫吸收反应动力学及模型的研究为氨法脱硫技术的应用和发展提供了很多有价值的参 考。然而，氨法脱硫过程中，吸收浆液的pH值沿吸收段高度方向上逐渐下降，液相中+4价硫 浓度及30 2、出503、11303与8032_之间的平衡均不断变化，50 2吸收是一个动态的过程，现有氨 法脱硫S02吸收模型的气相传质推动力的计算未从气-液动态平衡的角度出发，且对S0 2吸收 传质系数的简化程度较大，考察不全面，一般仅考虑了浆液的物理化学特性，从而导致不能 够真实、准确地反映氨法脱硫工艺中的传质过程，进而导致对氨法脱硫喷淋塔及脱硫工艺 的优化设计不理想，不能够有效地用于指导实际生产过程。
[0011] Codolo M C等(International Journal of Heat and Mass Transfer，2013,66: 80~89)对喷淋塔NaOH溶液吸收S02的体积传质系数进行了研究，考察了烟气流速、液体流 速、喷嘴直径等因素的影响。Kaji R等（Journal of Chemical Engineering of Japan， 1985，18:169~172)采用准数关联的方式研究了NaOH溶液吸收S02的传质系数。然而，上述 文献均是对喷淋塔NaOH溶液吸收S0 2的体积传质系数的研究，其体积传质系数是在特定的 实验条件下获得的，实验环境及条件与氨法脱硫工艺悬殊较大，因此获得的体积传质系数 在其它吸收设备中也就不具有通用性，不能够反映氨法脱硫工艺中的传质过程。此外，由于 考虑不全面，上述研究中得出的体积传质系数也不能真实地反映其传质过程。
[0012] 因此，研究并建立具有普遍适用性的氨法脱硫s〇2吸收传质系数的表达式，且能够 全面、准确的反映出不同工艺条件下s〇 2的吸收速率，对指导并优化氨法脱硫喷淋塔的设计 和运行，提高脱硫效率具有重要的理论意义和工程应用价值。

【发明内容】

[0013] 1.发明要解决的技术问题
[0014] 本发明的目的在于克服现有氨法脱硫喷淋塔实验系统易发生明显的氨的挥发性 损耗和硫酸铵气溶胶逃逸等问题，不能真实地模拟氨法脱硫工艺过程，以及现有氨法脱硫 S0 2吸收传质系数不能够真实、准确地反映氨法脱硫工艺中S02的吸收过程，从而导致不能用 于有效地指导氨法脱硫喷淋塔工程的设计和运行的不足，提供了一种氨法脱硫喷淋塔系统 及S0 2吸收传质系数的确定方法。本发明的氨法脱硫喷淋塔实验系统能够真实地模拟工业 中的氨法脱硫工艺过程，有效避免氨的挥发性损耗及硫酸铵气溶胶的逃逸现象;采用本发 明的s〇 2吸收传质系数的确定方法对s〇2吸收传质系数进行优化，得到的s〇2吸收传质系数能 够真实、准确地反映氨法脱硫工艺中s〇 2的吸收过程，且具有普遍适用性，从而能够为氨法 脱硫喷淋塔的精确设计和运行控制提供理论依据，有利于推动氨法烟气脱硫技术的发展。
[0015] 2 ?技术方案
[0016] 为达到上述目的，本发明提供的技术方案为：
[0017] 本发明的一种氨法脱硫喷淋塔系统，包括喷淋塔本体，该喷淋塔本体内部自上而 下依次包括除雾区、喷淋区和浆液槽，除雾区设有除雾器，除雾器的上方设有烟气出口，所 述喷淋区设有喷淋管，喷淋管的进水口通过循环管道及循环栗与浆液槽底部相连通，还包 括气压式加氨系统，该气压式加氨系统包括氨水罐和气栗，氨水罐的进气口通过管道与气 栗相连，且氨水罐的出液口通过管道与循环栗的输出口相连。
[0018] 更进一步地，所述喷淋塔本体侧壁底部的烟气进口通过进气管与混合罐相连，该 混合罐分别通过管道与S02钢瓶及氧化风机相连，来自S0 2钢瓶的S02气源与来自氧化风机的 空气经混合罐混合均匀后通过喷淋塔本体侧壁底部的烟气进口管道进入喷淋塔本体内部。
[0019] 本发明的一种氨法烟气脱硫S02吸收传质系数的计算方法，该方法基于喷淋塔氨 法烟气脱硫工艺的特点，选取喷淋塔内浆液pH值、烟气流速udP液气比L/Q为参数，先计算 出不同pH值、烟气流速％及液气比L/Q条件下的氨法脱硫S0 2吸收传质速率ky，然后采用数值 拟合的方法重新构建出具有普遍适用性的氨法脱硫S0 2吸收传质系数k/的表达式。
[0020] 更进一步地，所述喷淋塔内浆液pH值的范围为4.5~6.0，烟气流速％的范围为2.0 ~4.0m/s，液气比L/Q的范围为2.0~4.0L/m 3。
[0021] 更进一步地，该计算方法的具体步骤如下：
[0022]步骤一、确定如下参数:氨法脱硫系统喷淋塔的D、G、L、pH、P、Q、T、yi、y2、Z、P、y、v、 〇;
[0023]其中，D为喷淋塔直径，通过直接测量得出，单位m;L为浆液喷淋量，通过流量计测 量得出，单位为L/h;pH为喷淋塔内浆液pH值，由pH计测量得到;P为烟气总压力，近似与大气 压相等，为l〇1325Pa; Q为烟气流量，通过流量计测量得出，单位m3/h; T为烟气温度，通过温 度传感器测量得出，单位K;G为烟气摩尔流率，根据公式G=1(T3 ? PQ/(RTjtD2)/4计算得出， 单位为kmol/(m2 ? h) ;yi、y2分别为喷淋塔进口、出口烟气中S〇2的摩尔分率，由烟气分析仪 测得，单位均为ppm; Z为喷淋塔吸收区高度，直接测量得到，单位m; P为喷淋塔内浆液密度， 通过密度计测量得到，单位kg/m3;浆液的粘度y和表面张力〇分别利用粘度计和表面张力仪 测得，其单位分别为Pa ? s和N/m;v为浆液的运动粘度，根据公式V = y/P计算得出，单位m2/ s；
[0024]步骤二、氨法烟气脱硫S02吸收传质速率的计算：启动氨法脱硫喷淋塔系统，改变 浆液pH值、烟气流速ug和浆液喷淋量L，测定不同pH值、烟气流速％和液气比L/Q条件下喷淋 塔进、出口烟气中S0 2的浓度yi、y2，通过对液滴运动和塔壁液膜流动状态的计算，确定气-液 传质接触面积a，进而计算出不同实验条件下的氨法脱硫S0 2吸收传质速率ky;
[0025]步骤三、氨法脱硫S02吸收传质系数表达式的构建：
[0026]选取pH值、烟气流速也和液气比L/Q为参数，构建S02吸收传质系数的表达式如下：
[0027] :QT' ⑴，
[0028] 上式中ko为系数，X1、X2和X3均为指数，对上式两边取对数可得：
[0029] In k/y=ln ko+xiln pH+X2ln ug+X3ln(L/Q) (2)，
[0030]将步骤二中通过测量计算所得的氨法脱硫S〇2吸收传质速率ky与pH、ug及L/Q按公 式(2)进行多元线性拟合，确定系数k〇及指数X1、X2和X3的值，进而重新构建出氨法脱硫50 2吸 收传质系数k/的表达式。
[0031]更进一步地，步骤二中不同条件下的氨法脱硫s〇2吸收传质速率ky根据如下公式进 行计算：
[0033]式(3)中，a为喷淋塔内气-液传质接触面积，单位为m2/m3，包含两部分:一是喷淋塔 内液滴的比表面积Ad/ (( JtD2/4 ) Z)，m2/m3;二是喷淋塔壁液膜的比表面积Am/ (( jtD2/4) Z)，m2/ m3,因此公式(3)可表述为如下形式：
[0035]更进一步地，所述喷淋塔微元内液滴的表面积Ad根据如下公式进行计算：
[0037]上式(5)中，Ld为以液滴形式存在的浆液喷淋量，单位为L/h; Ud为液滴的下落速度， 单位m/s; d为液滴的Sauter平均直径，单位为m，液滴直径d根据如下公式(6)进行计算： [0038] d = 133.0.(D/2).We-0.74 (6)，
[0039]上式中，韦博数We根据公式(7)进行计算：
[0041]上式中，4为喷嘴衆液的喷出速度，单位为m/ s，根据动量守丨旦方程/hi 〃丨;喝?! 2进 行计算，其中，nu为喷嘴喷出浆液的质量流量，单位为kg/s，通过测量得到为喷嘴运行压 力，单位为MPa，通过测量得到;d k为喷嘴孔径，单位m; 〇为浆液表面张力，N/m;
[0042] 公式(5)中液滴下落速度Ud表述为：
[0044] 该公式中，g为重力加速度，单位为m/s2;Pd、Pg分别为液滴密度和气体密度，单位均 为kg/m 3，其中液滴密度Pd即为喷淋塔内浆液密度P; ug为气体速度，单位m/s，通过下式计算 得到：
[0045] ug = Q/(360〇JiD2)/4 (9)，
[0046]公式(8)中Cd为曳力系数，根据下式进行计算：
[0048]上式中，雷诺数Red根据公式（11)进行计算：
[0050]液滴的运动方程为：
[0051 ] uAdt = dz ( 12)?
[0052]联立公式(5)~（12)，沿喷淋塔吸收区高度方向上积分即计算出液滴表面积Ad; [0053]所述喷淋塔壁液膜的面积^按照如下公式进行计算：
[0054] Am=^(D-2hm) ? dz (13)，
[0055] 式(13)中，hm为喷淋塔壁液膜厚度，单位为m，根据如下公式进行计算：
[0057]上式中，Rem为液膜流动雷诺数，根据如下公式进行计算：
[0058] Rem=4LmP/jidli (15)，
[0059] 公式（15)中，Lm为塔壁液体流量，L/h;联立公式（13)~（15)，沿喷淋塔吸收区高度 方向上积分即计算出喷淋塔壁液膜面积A m。
[0060] 更进一步地，塔壁液体流量Lm和以液滴形式存在的浆液量Ld分别近似为总浆液流 量L 的 8%、92%。
[0061 ]如表1所示为本发明中的参数说明。
[0062]表1本发明中的参数说明

[0065] 3 ?有益效果
[0066] 采用本发明提供的技术方案，与现有技术相比，具有如下显著效果：
[0067] (1)本发明的一种氨法脱硫喷淋塔系统，通过气压式加氨系统对向系统内补充吸 收剂順 3,从而可以有效克服采用高位槽和栗来加氨带来的氨的挥发问题。通过该喷淋塔系 统对工业中的氨法脱硫工艺进行模拟，能够更加真实地反映氨法脱硫工艺过程，从而为工 业中氨法脱硫喷淋塔的优化设计和运行提供可靠依据。
[0068] (2)本发明的氨法烟气脱硫S02吸收传质系数的确定方法，该方法基于喷淋塔氨法 烟气脱硫工艺的特点，选取喷淋塔内浆液pH值、烟气流速udP液气比L/Q为参数，先计算出 不同pH值、烟气流速％及液气比L/Q条件下的氨法脱硫S0 2吸收传质速率ky，然后采用数值拟 合的方法重新构建出具有普遍适用性的氨法脱硫S0 2吸收传质系数k/的表达式。采用本发 明的方法得到的氨法烟气脱硫S02吸收传质系数k/能够更加真实、准确地反映氨法脱硫工 艺中的S0 2传质吸收过程，从而能够直接、有效地用于指导氨法脱硫喷淋塔及脱硫工艺的优 化设计与运行，有利于提高氨法脱硫工艺的脱硫效率，更具有实际意义和工程应用价值。
[0069] (3)本发明的氨法烟气脱硫S02吸收传质系数的确定方法，通过对喷淋塔内液滴运 动和塔壁液膜流动的理论分析，在计算氨法脱硫S0 2吸收传质速率的过程中全面考虑了喷 淋液滴面积和塔壁液膜面积的影响，从而保证建立的传质系数表达式不只是针对特定的喷 淋吸收塔有效，具有普遍适用性。
[0070] (4)本发明的氨法烟气脱硫S02吸收传质系数的确定方法，通过该方法得到的氨法 烟气脱硫S0 2吸收传质系数的表达式更加准确，对该表达式进行试验验证，结果表明采用本 发明获得的表达式进行计算得到的计算值与实验结果之间的误差较小，小于14.77%。
【附图说明】
[0071 ]图1为本发明的氨法脱硫喷淋塔系统的结构示意图；
[0072]图2为本发明的喷淋塔氨法脱硫的物理模型图；
[0073]图3为本发明的喷淋塔氨法脱硫S02吸收物料的平衡示意图。
[0074] 不意图中的标号说明：
[0075] 1、S02钢瓶;2、减压阀；3、稳流器;4、风机;5、流量计;6、取样口； 7、混合罐;8、气栗； 9、除雾器;10、喷淋塔;11、浆液槽;12、循环栗;13、控制柜;14、氨水罐;15、烟气分析仪。
【具体实施方式】
[0076] 为进一步了解本发明的内容，现结合附图和实施例对本发明作详细描述。
[0077] 实施例1
[0078] 1、氨法脱硫S02吸收传质系数的确定
[0079]本实施例的氨法脱硫喷淋塔系统如图1所示，喷淋塔直径D为0.3m，吸收区高度Z为 2.0m。该氨法脱硫喷淋塔系统包括喷淋塔本体10，该喷淋塔本体10内部自上而下依次包括 除雾区、喷淋区和浆液槽11，除雾区设有除雾器9，除雾器9的上方设有烟气出口，该烟气出 口与尾气排放管道相连，且尾气排放管道上设有控制柜13。所述喷淋区设有喷淋管，该喷淋 管的进水口通过循环管道及循环栗12与浆液槽11底部相连通。浆液槽11底部还通过管道与 气栗8相连。来自S0 2钢瓶1的S02气源与来自氧化风机4的空气经混合罐7混合均匀后通过喷 淋塔本体10侧壁底部的烟气进口管道进入喷淋塔本体10内部，由喷淋塔下部向上流过喷淋 区并与喷淋而下的液滴逆流接触，以(NH4) 2S〇3-NH4HS〇3混合溶液为基础对s〇2进行吸收。 [0080]浆液中，对s〇 2起吸收作用的是(NH4)2S〇3，NH4HS〇3不具有吸收s〇2的能力。随着吸收 反应的进行，吸收液中NH 4HS〇3的浓度逐渐增大，浆液吸收s〇2的能力开始下降，为了维持浆 液的吸收能力，需要向系统内补充吸收剂順 3,使部分NH4HS03转化为(NH4)2S〇3。吸收剂NH3的 补给采用气压式加氨系统，该气压式加氨系统包括氨水罐14和气栗8,氨水罐14的进气口通 过管道与气栗8相连，且氨水罐14的出液口通过管道与循环栗12的输出口相连，利用气栗8 提供的空气的压力将氨水(含量25~28wt. %)从氨水罐14中压出，经氨水罐14的出液管道 上的流量计5调节流量后，在浆液循环栗12入口处与浆液进行混合，有效地克服了采用高位 槽和栗来加氨带来的氨的挥发问题。氧化空气量根据单位时间吸收的S0 2的量确定，由气栗 8来提供。本实施例中，所述S02钢瓶1与混合罐7之间的管道上设有减压阀2和稳流器3,所述 氧化风机4的出气管道、气栗8与浆液槽11之间的管道以及氨水罐14的出液管道均设有流量 计5。所述喷淋塔本体10侧壁底部的烟气进口管道、与喷淋管相连的循环管道以及喷淋塔本 体10顶部的尾气排放管道上均设有取样口 6,通过烟气分析仪15对喷淋塔本体10的进口和 出口烟气进行分析。通过该喷淋塔系统对工业中的氨法脱硫工艺进行模拟，能够更加真实 地反映氨法脱硫工艺过程，从而为工业中氨法脱硫喷淋塔的优化设计和运行提供可靠依 据。
[0081 ]本实施例的氨法烟气脱硫S02吸收传质系数的确定方法，该方法基于喷淋塔氨法 烟气脱硫工艺的特点，选取喷淋塔内浆液pH值、烟气流速udP液气比L/Q为参数，先计算出 不同pH值、烟气流速也及液气比L/Q条件下的氨法脱硫S0 2吸收传质速率ky，然后采用数值拟 合的方法重新构建出具有普遍适用性的氨法脱硫s〇2吸收传质系数k/的表达式，该方法适 用于喷淋塔内浆液pH值的范围为4.5~6.0，烟气流速ug的范围为2.0~4. Om/S，液气比L/Q 的范围为2.0~4. OL/m3的情况，其具体步骤为：
[0082]步骤一、确定如下参数:氨法脱硫系统喷淋塔的D、G、L、pH、P、Q、T、yi、y2、Z、P、y、v、 〇;
[0083]其中，D为喷淋塔直径，通过直接测量得出，单位m;L为浆液喷淋量，通过流量计测 量得出，单位为L/h;pH为喷淋塔内浆液pH值，由pH计测量得到;P为烟气总压力，近似与大气 压相等，为l〇1325Pa; Q为烟气流量，通过流量计测量得出，单位m3/h; T为烟气温度，通过温 度传感器测量得出，单位K;G为烟气摩尔流率，根据公式G=1(T3 ? PQ/(RTjtD2)/4计算得出， 单位为kmol/(m2 ? h)，该式中R为气体状态常数，单位为J/(mol ? K) ;yi、y2分别为喷淋塔进 口、出口烟气中S02的摩尔分率，由烟气分析仪测得，单位均为ppm;Z为喷淋塔吸收区高度， 直接测量得到，单位m; P为喷淋塔内浆液密度，通过密度计测量得到，单位kg/m3;浆液的粘 度y和表面张力〇分别利用粘度计和表面张力仪测得，其单位分别为Pa ? s和N/m;v为浆液的 运动粘度，根据公式v=lVf?计算得出，单位m2/s。
[0084] 步骤二、氨法烟气脱硫S02吸收传质速率的计算：
[0085] 启动氨法脱硫喷淋塔系统，改变浆液pH值、烟气流速ug和浆液喷淋量L，测定不同 pH值、烟气流速％和液气比L/Q条件下喷淋塔进、出口烟气中S0 2的浓度yi、y2，通过对液滴运 动和塔壁液膜流动状态的计算，确定气-液传质接触面积a，进而计算出不同实验条件下的 氨法脱硫S0 2吸收传质速率ky，如表2所示为不同pH值、烟气流速udP液气比L/Q条件下计算 得到的S0 2吸收传质速率ky。具体的，本实施例以表2中序号1实验为例，详细描述氨法脱硫 S0 2吸收传质速率ky的计算过程。控制实验条件浆液pH为5.5，烟气流速％为1.76m/s，浆液喷 淋量L为1350L/h，液气比3L/m 3，T为323.151(，实验测得进口302浓度71为36(^111，出口30 2浓 度y2为 184ppm。
[0086]氨法脱硫S02吸收传质速率ky根据如下公式进行计算：
[0088]式(3)中，a为喷淋塔内气-液传质接触面积，单位为m2/m3,因此要想计算得到吸收 传质速率ky，必须要先计算出喷淋塔内气-液传质接触面积a，喷淋塔内气-液传质接触面积 a包含两部分:一是喷淋塔内液滴的比表面积Ad/(UD 2/4)Z)，m2/m3;二是喷淋塔壁液膜的比 表面积Am/((3iD 2/4)Z)，m2/m3,因此公式(3)可表述为如下形式：
[0090]沿喷淋塔高度方向上将吸收塔划分为多个微元，由于微元高度很小，假定微元内 液滴下落的速度不变，另外，忽略液滴的凝并和破碎，液滴直径采用Sauter平均直径。因此， 喷淋塔微元内液滴的表面积Ad可表述为：
[0092]上式(5)中，Ld为以液滴形式存在的浆液喷淋量，单位为L/h，近似为总浆液流量L 的92% dd为液滴的下落速度，单位m/s; d为液滴的Sauter平均直径，单位为m，液滴直径d根 据如下公式(6)进行计算：
[0093] d = 133.0.(D/2).We-0.74 (6)，
[0094] h式中，韦博数We枏据公式（7)讲行计筧：
[0096]上式中，为喷嘴衆液的喷出速度，单位为m/s，根据动量守丨旦方程wi.4=户进 行计算，其中，nu为喷嘴喷出浆液的质量流量，单位为kg/s，通过测量得到;Pi为喷嘴运行压 力，单位为MPa，通过测量得到;d k为喷嘴孔径，单位m，直接测量得到；〇为浆液表面张力，N/ m。本实施例中，浆液密度P为1273 ? 9kg ? nf3，浆液表面张力〇为0 ? 073N/m，喷嘴运行压力PiS 0.4MPa，喷嘴孔径dk为4 X 10-3m，计算得到4为15m/s，联立式(6)和(7)，计算出液滴的直径d 为1.3 X l(T3m，考虑到液滴的凝并，取液滴直径d为2.0 X l(T3m，从而减小计算误差。
[0097] 公式(5)中液滴下落速度Ud表述为：
[0099] 该公式中，g为重力加速度，单位为m/s2;Pd、Pg分别为液滴密度和气体密度，单位均 为kg/m 3，其中液滴密度Pd即为喷淋塔内浆液密度P; ug为气体速度，单位m/s，通过下式计算 得到：
[0100] ug = Q/(360〇JiD2)/4 (9)，
[0101]公式(8)中Cd为曳力系数，根据下式进行计算：
[0103]上式中，雷诺数Red根据公式(11)进行计算：
[0105] 上式中，气体粘度yg通过查《化学化工物性手册》获得;液滴的运动方程为：
[0106] ￡\由=tfe (12),
[0107] 通过查询《化学化工物性手册》得到烟气密度Pg取1.23kg/m3,烟气流速u gSl .76m/ 8，烟气粘度1%为1.81\10-午&%，微元高度如取0.01111，联立公式(5)~（12)，沿喷淋塔吸收 区高度方向上积分即计算出液滴表面积Ad为0.199m 2。
[0108] 喷淋塔氨法脱硫过程中，部分浆液落在喷淋塔壁上并形成液膜，塔壁液膜在与烟 气接触的过程中参与S02吸收反应，喷淋塔壁液膜的面积^按照如下公式进行计算：
[0109] Am=^(D-2hm) ? dz (13)，
[0110] 式（13)中，hm为喷淋塔壁液膜厚度，单位为m，喷淋塔壁液膜按层流流动考虑，其厚 度hm可按Feind公式计算：
[0112] 上式中，Rem为液膜流动雷诺数，根据如下公式进行计算：
[0113] Rem=4LmP/3idy (15)，
[0114] 公式（15)中，Lm为塔壁液体流量，L/h，近似为总浆液流量L的8%，浆液运动粘度V 为0.79X10-6m2/s，浆液粘度y为1.005 X 10-3Pa ? s，微元高度为0.01m，联立公式（13)~ (15)，沿喷淋塔吸收区高度方向上积分即计算出喷淋塔壁液膜厚度为0.0031m，喷淋塔壁液 膜面积Am为1.84m 2。
[0115] 将实验测得的喷淋塔进、出口 S02浓度ydPy2,以及计算得到的液滴表面积Ad和塔 壁液膜面积Am带入式(4)，可以计算得到S0 2吸收传质速率ky为3.11kmol/(m2 ? h)。同理，改 变pH、烟气流速％和浆液喷淋量L，计算出不同实验条件下的氨法脱硫S02吸收传质速率k y。
[0116] 本实施例中，氨法脱硫S02吸收传质速率ky的计算式即公式(3)的推导过程如下：
[0117] 氨法脱硫过程中，稳态条件下S02的传质速率可表述为：
[0118] nSB2 =kap-p：〇,) (16)，
[0119] 其中，Kg为气相总传质速率，p和!^分别为S〇2在气相主体和相界面上的分压，根据 双膜理论，气相总传质速率Kg与气、液相传质分速率kt；、ki的关系可表述为：
[0121] 上式中，H为S〇2溶解度系数，单位为mol/m3。
[0122] 由于S02易溶于水，S02溶解度系数较大，氨法脱硫S02液相传质阻力可以忽略不计， 传质总阻力等于气相传质分阻力，式(17)可表达为：
[0124] 因此，式（16)可转化为：
[0125] Nm： =K(p-pi〇,) = /f, (,v-.v,) (19)，
[0126] 在吸收塔高度方向上dz微元内对S02进行物料衡算(如附图2、附图3所示），得到如 下方程：
[0128] 式中，a为传质接触面积，m2/m3;y为气相主体中S02的摩尔分率，yi为气液界面处S0 2 的平衡摩尔分率;G为烟气摩尔流率，kmol/(m2 ? h)。
[0129] 式(20)可转化为：
[0131] 由溶液中H2S03、HSO;和SOf的分布系数（公开号：CN 105032154A)可知，在氨法脱 硫的pH范围（pH为4.0~6.0)内，溶液中的硫主要以1 ISa形式存在，H2S03的量极少，几乎没 有游离态的S02,因而气相中的S0 2的体积分数远远大于气液界面的S02的体积分数，即y>> yi。另外，燃煤烟气中S〇2的浓度较低，式(21)可以简化为：
[0132] Gdy = kyaydz (22),
[0133] 式(22)可转化为：
[0135]沿吸收塔高度上从上向下对式(23)两边进行积分(如附图3所示）：
[0137]式中，yi、y2分别为喷淋吸收塔进、出口烟气中S〇2的摩尔分率(体积比浓度）；Z为喷 淋塔吸收区高度，m;由积分式(24)计算可得：
[0139] 因此，由公式(25)即可推导出气相分传质速率ky的表达式：
[0140] 表2不同实验条件下氨法脱硫S02吸收传质速率

[0143] 步骤三、氨法脱硫S02吸收传质系数表达式的构建：
[0144] 选取pH值、烟气流速％和液气比L/Q为参数，构建S02吸收传质系数的表达式如下：
[0145] -p/gp (丄），
[0146] 上式中k〇为系数，XI、X2和x3均为指数，对上式两边取对数可得：
[0147] In k/y=ln k〇+xi In pH+X2 In ug+X3 ln(L/Q) (2)，
[0148] 将步骤二中通过测量计算所得的氨法脱硫S02吸收传质速率ky与pH、ug&L/Q(具体 数据见表2所示)按公式(2)进行多元线性拟合，确定系数k Q及指数X1、X2和X3的值，进而重新 构建出氨法脱硫s〇2吸收传质速率ky的表达式，本实施例重新构建的氨法脱硫s〇2吸收传质 速率k/的表达式为：6
[0149] 鉴于现有技术中氨法脱硫喷淋塔及脱硫工艺的优化设计效果在很大程度上依赖 于氨法烟气脱硫S02吸收传质系数表达式的准确性，
[0150] 由于现有技术中氨法烟气脱硫S02吸收传质系数的表达式不能够真实、准确地反 映S02吸收传质过程，从而不能有效用于指导氨法脱硫喷淋塔工程的优化设计和运行，导致 氨法脱硫工艺的脱硫效率不能得到有效提高。基于此，发明人一直致力于氨法烟气脱硫工 艺的研究，发明人通过大量实验研究及实践经验发现，氨法烟气脱硫工艺中S0 2吸收传质过 程主要受喷淋塔内浆液pH值、烟气流速％和液气比L/Q的影响较大，因此其选用这三个参数 为计算参数，重新构建出氨法烟气脱硫S0 2吸收传质速率与pH值、ug&L/Q之间的计算表达 式。采用本发明的方法得到的氨法烟气脱硫S0 2吸收传质系数能够更加真实、准确地反映氨 法脱硫工艺中的S02传质吸收过程，从而能够直接、有效地用于指导氨法脱硫喷淋塔及脱硫 工艺的优化设计与运行，更具有实际意义和工程应用价值。
[0151] 此外，本发明的氨法烟气脱硫S02吸收传质系数的计算方法，通过对喷淋塔内液滴 运动和塔壁液膜流动的理论分析，在计算氨法脱硫S0 2吸收传质系数的过程中全面考虑了 喷淋液滴面积和塔壁液膜面积的影响，建立的传质系数表达式不只是针对特定的喷淋吸收 塔有效，具有普遍适用性。
[0152] 2、氨法脱硫S02吸收传质系数表达式的实验验证
[0153] 采用如图1所示的喷淋塔氨法脱硫实验系统，控制浆液pH为5.5，烟气流速ug为1.76m/ s，液气比为3L/m 3，实验测得喷淋塔进口 S02浓度为360ppm，喷淋塔出口 S02浓度为184ppm，联立式 (3)~(15)计算得到的S02吸收传质系数为3.11kmol/(m 2 ? h)。将实验条件pH为5.5，Uf^1.76m/S， 液气比为3L/m3代入建立的氨法脱硫S02传质速率表达式：x續严，计 算得到S02吸收传质系数为3.08kmol/(m2 ?!〇,二者之间的误差为1.1%。因此，采用本发明 的方法得到的S02吸收传质系数较为准确，能够较好地反映的S0 2吸收传质过程。
[0154] 实施例2
[0155] 采用如实施例1中的喷淋塔氨法脱硫实验系统，控制浆液pH为5.5，烟气流速ug为 3.34m/s，液气比为3L/m 3，实验测得喷淋塔进口 S02浓度为360ppm，喷淋塔出口 S02浓度为 233ppm，联立式(3)~（15)计算得到的S〇2吸收传质系数为3.545kmol/(m 2 ? h)，计算方法同 实施例1。将实验条件pH为5.5，％为3.34m/s，液气比为3L/m3代入建立的氨法脱硫S0 2传质速 率表达式，计算得到S02吸收传质系数为3.54kmol/(m2 ? h)，二者之间的误差仅为0.13%。
[0156] 实施例3
[0157] 采用如实施例1中的喷淋塔氨法脱硫实验系统，控制浆液pH为5.68，烟气流速％为 2.16m/s，液气比为3L/m3,实验测得喷淋塔进口 S02浓度为360ppm，喷淋塔出口 S02浓度为 183.7ppm，联立式(3)~（15)计算得到的S〇2吸收传质系数为3.75kmol/(m 2 ? h)，计算方法同 实施例1。将实验条件pH为5.68，％为2.16m/s，液气比为3L/m3代入建立的氨法脱硫S0 2传质 速率表达式，计算得到S02吸收传质系数为3.60kmol/(m2 ? h)，二者之间的误差为3.94%。
[0158] 实施例4
[0159]采用如实施例1中的喷淋塔氨法脱硫实验系统，控制浆液pH为5.95，烟气流速％为 2.16m/s，液气比为3L/m3,实验测得喷淋塔进口 S02浓度为360ppm，喷淋塔出口 S02浓度为 162ppm，联立式(3)~（15)计算得到的S〇2吸收传质系数为4.45kmol/(m2 ? h)，计算方法同实 施例1。将实验条件pH为5.95，％为2.16m/s，液气比为3L/m3代入建立的氨法脱硫S0 2传质速 率表达式，计算得到S〇2吸收传质系数为4.23kmol/(m2 ? h)，二者之间的误差为4.86%。 [0160] 实施例5
[0161 ]采用如实施例1中的喷淋塔氨法脱硫实验系统，控制浆液pH为5.5，烟气流速ug为 2.16m/s，液气比为2.5L/m3，实验测得喷淋塔进口 S02浓度为360ppm，喷淋塔出口 S02浓度为 235ppm，联立式(3)~（15)计算得到的S〇2吸收传质系数为2.41kmol/(m 2 ? h)，计算方法同实 施例1。将实验条件pH为5.5，ugS2.16m/s，液气比为2.5L/m 3代入建立的氨法脱硫S02传质速 率表达式，计算得到S02吸收传质系数为2.83kmol/(m 2 ? h)，二者之间的误差为14.6%。
[0162] 实施例6
[0163] 采用如实施例1中的喷淋塔氨法脱硫实验系统，控制浆液pH为5.5，烟气流速ug为 2.16m/s，液气比为3L/m 3,实验测得喷淋塔进口 S02浓度为360ppm，喷淋塔出口 S02浓度为 220ppm，联立式(3)~（15)计算得到的S〇2吸收传质系数为2.74kmol/(m 2 ? h)，计算方法同实 施例1。将实验条件pH为5.5，ugS2.16m/s，液气比为3L/m3代入建立的氨法脱硫S0 2传质速率 表达式，计算得到302吸收传质系数为3.221〇11〇1/( 1112?11)，二者之间的误差为14.77%。
【主权项】
1. 一种氨法脱硫喷淋塔系统，包括喷淋塔本体（10)，该喷淋塔本体（10)内部自上而下 依次包括除雾区、喷淋区和浆液槽(11)，除雾区设有除雾器(9)，除雾器(9)的上方设有烟气 出口，所述喷淋区设有喷淋管，喷淋管的进水口通过循环管道及循环栗（12)与浆液槽（11) 底部相连通，其特征在于:还包括气压式加氨系统，该气压式加氨系统包括氨水罐(14)和气 栗(8)，氨水罐（14)的进气口通过管道与气栗(8)相连，且氨水罐（14)的出液口通过管道与 循环栗（12)的输出口相连。2. 根据权利要求1所述的一种氨法脱硫喷淋塔系统，其特征在于:所述喷淋塔本体（10) 侧壁底部的烟气进口通过进气管与混合罐(7)相连，该混合罐(7)分别通过管道与SO 2钢瓶 (1)及氧化风机(4)相连，来自SO2钢瓶（1)的SO2气源与来自氧化风机(4)的空气经混合罐(7) 混合均匀后通过喷淋塔本体(10)侧壁底部的烟气进口管道进入喷淋塔本体(10)内部。3. -种氨法烟气脱硫SO2吸收传质系数的确定方法，其特征在于:该方法基于喷淋塔氨 法烟气脱硫工艺的特点，选取喷淋塔内浆液PH值、烟气流速％和液气比L/Q为参数，先计算 出不同PH值、烟气流速％及液气比L/Q条件下的氨法脱硫SO 2吸收传质速率ky，然后采用数值 拟合的方法重新构建出具有普遍适用性的氨法脱硫SO 2吸收传质系数k/的表达式。4. 根据权利要求3所述的一种氨法烟气脱硫SO2吸收传质系数的确定方法，其特征在于： 所述喷淋塔内浆液pH值的范围为4.5~6.0，烟气流速％的范围为2.0~4. Om/s，液气比L/Q 的范围为2.0~4.0L/m3。5. 根据权利要求3或4所述的一种氨法烟气脱硫SO2吸收传质系数的确定方法，其特征在 于:其具体步骤如下： 步骤一、确定如下参数:氨法脱硫系统喷淋塔的D、G、L、pH、P、Q、T、yi、y2、Z、P、y、v、 〇; 其中，D为喷淋塔直径，通过直接测量得出，单位m; L为浆液喷淋量，通过流量计测量得 出，单位为L/h;pH为喷淋塔内浆液pH值，由pH计测量得到;P为烟气总压力，近似与大气压相 等，为101325Pa; Q为烟气流量，通过流量计测量得出，单位m3/h; T为烟气温度，通过温度传 感器测量得出，单位K;G为烟气摩尔流率，根据公式G=HT3 · PQ/(RTjtD2)/4计算得出，单位 为kmol/(m2 别为喷淋塔进口、出口烟气中SO 2的摩尔分率，由烟气分析仪测得， 单位均为ppm; Z为喷淋塔吸收区高度，直接测量得到，单位m; P为喷淋塔内浆液密度，通过密 度计测量得到，单位kg/m3;浆液的粘度μ和表面张力σ分别利用粘度计和表面张力仪测得， 其单位分别为Pa · s和N/m;v为楽:液的运动粘度，根据公式ν=μ/ρ计算得出，单位m2/s; 步骤二、氨法烟气脱硫SO2吸收传质速率的计算： 启动氨法脱硫喷淋塔系统，改变浆液pH值、烟气流速Ug和浆液喷淋量L，测定不同pH值、 烟气流速％和液气比L/Q条件下喷淋塔进、出口烟气中SO2的浓度yi、y2，通过对液滴运动和 塔壁液膜流动状态的计算，确定气-液传质接触面积a，进而计算出不同实验条件下的氨法 脱硫SO 2吸收传质速率ky; 步骤三、氨法脱硫SO2吸收传质系数表达式的构建： 选取pH值、烟气流速％和液气比L/Q为参数，构建SO2吸收传质系数的表达式如下： -a/gF (1)， 上式中k〇为系数，Xl、X2和X3均为指数，对上式两边取对数可得： Ink7 y= lnk〇+xilnpH+X2lnug+X3ln(L/Q) (2)， 将步骤二中通过测量计算所得的氨法脱硫SO2吸收传质速率ky与pH、ug&L/Q按公式(2) 进行多元线性拟合，确定系数k〇及指数XI、X2和X3的值，进而重新构建出氨法脱硫S〇2吸收传 质系数k/的表达式。6. 根据权利要求5所述的一种氨法烟气脱硫SO2吸收传质系数的确定方法，其特征在于： 步骤二中不同条件下的氨法脱硫SO 2吸收传质速率ky根据如下公式进行计算：式(3)中，a为喷淋塔内气-液传质接触面积，单位为m2/m3，包含两部分:一是喷淋塔内液 滴的比表面积Ad/(( JtD2/4 )Z)，m2/m3;二是喷淋塔壁液膜的比表面积Am/(( jtD2/4) Z)，m2/m3，因 此公式(3)可表述为如下形式：7. 根据权利要求6所述的一种氨法烟气脱硫SO2吸收传质系数的确定方法，其特征在于： 所述喷淋塔微元内液滴的表面积Ad根据如下公式进行计算：上式(5)中，Ld为以液滴形式存在的浆液喷淋量，单位为L/h; Ud为液滴的下落速度，单位 m/s; d为液滴的Sauter平均直径，单位为m，液滴直径d根据如下公式(6)进行计算： d=133.0 · (D/2) · We-Q.74 (6)， 上式中，韦博数We根据公式(7)进行计算：上式中，为喷嘴衆液的喷出速度，单位为m/s，根据动量守丨旦方程=P1TicZ42进行计 算，其中，πη为喷嘴喷出浆液的质量流量，单位为kg/s，通过测量得到;^为喷嘴运行压力，单 位为MPa，通过测量得到;d k为喷嘴孔径，单位m; 〇为浆液表面张力，N/m; 公式(5)中液滴下落速度Ud表述为：C8), 该公式中，g为重力加速度，单位为m/s2;Pd、Pg分别为液滴密度和气体密度，单位均为 kg/m3，其中液滴密度Pd即为喷淋塔内浆液密度P; Ug为气体速度，单位m/s，通过下式计算得 到： ug = Q/(360〇3iD2)/4 (9)， 公式(8)中Cd为曳力系数，根据下式进行计算：上式中，雷诺数Red根据公式（11)进行计算：(12)， 输逾的;云云fl卞*?呈为： 联立公式(5)~（12)，沿喷淋塔吸收区高度方向上积分即计算出液滴表面积Ad; 所述喷淋塔壁液膜的面积照如下公式进行计算：式(13)中，hm为喷淋塔壁液膜厚度，单位为m，根据如下公式进行计算：上式中，Rem为液膜流动雷诺数，根据如下公式进行计算： Rem - 4LmP/兀(15)， 公式（15)中，Lm为塔壁液体流量，L/h;联立公式（13)~（15 )，沿喷淋塔吸收区高度方向 上积分即计算出喷淋塔壁液膜面积Am。8.根据权利要求7所述的一种氨法烟气脱硫SO2吸收传质系数的确定方法，其特征在于： 塔壁液体流量Lm和以液滴形式存在的浆液量Ld分别近似为总浆液流量L的8%、92%。
【文档编号】B01D53/50GK106000041SQ201610384391
【公开日】2016年10月12日
【申请日】2016年5月30日
【发明人】贾勇, 殷李国, 吴胜华, 陈 光, 李智芳, 丁希楼, 陈宜华, 戴波
【申请人】安徽工业大学