以垃圾焚烧厂脱酸塔为模拟对象，建立了脱酸塔气固液之间的传质模型及化学反应模型，对旋转产生的喷雾建立了浆滴模型以分析其蒸发和传质过程，通过对已有实验模型与脱酸塔实际运行结果进行对比，验证了所建化学反应模型和浆滴模型的准确性，模拟了某垃圾焚烧厂脱酸塔的内流特征，得到其压力、温度、速度和酸性气体质量分数等参数的分布及脱除率。
结果表明：脱酸塔中轴线上部是脱酸反应主要发生部位;由于烟气回流区的作用，部分没有蒸发完全的浆滴被吹向脱酸塔壁面，容易造成脱酸塔壁面的结渣和腐蚀，导致脱酸塔壁面所受的热应力增大，不利于脱酸塔长期稳定运行;浆滴进入脱酸塔5s内完全蒸发。
垃圾焚烧技术已成为目前解决日益严重的垃圾污染问题、缓解能源危机的重要措施。但焚烧过程产生大量的酸性气体(如SO2和HCl等)，使得垃圾焚烧烟气净化技术的研究与设备优化备受关注。
其烟气净化技术主要借鉴燃煤烟气净化方法，一般采用半干式喷雾干燥吸收法。但垃圾焚烧厂与燃煤电厂烟气净化的主要区别在于脱除对象不同：燃煤电厂主要脱除SO2，而垃圾焚烧厂主要脱除HCl和SO2，且在脱除过程中HCl与SO2之间存在复杂的相互影响和竞争关系。
脱酸塔作为烟气净化的重要设备，气液固在塔内的流动、传热传质过程对脱酸塔的设计和运行具有重要的影响。臧仁德等研究了垃圾焚烧中喷雾干燥塔内烟气的温度分布，指出增加烟气流速可降低石灰浆滴(以下简称浆滴)对SO2、HCl的脱除率，且SO2脱除率降低幅度更显著。
李午申对双流体喷雾干燥塔内酸性气体净化进行了研究，指出SO2和HCl在半干法烟气净化中的净化效率随吸收剂用量的增大而提高。Kiil等开展了脱酸塔装置的实验研究，分析了不同高度上的脱酸效率，指出石灰石剩余量越大，HCl脱除率越高。
因垃圾焚烧中脱除气体不同，所以在模型建立与气体有关部分的处理与燃煤电厂将存在较大差异，如气体浓度场、离散相与连续相的传质等过程，因此不能简单套用脱硫模型分析垃圾烟气脱酸过程。
燃煤电厂烟气脱硫过程大多采用双模理论模拟吸收过程，将传质过程简化为通过两层薄膜的扩散过程，其在反映客观实际方面仍有缺点和局限性，且传质过程大多用单纯液滴(H2O)或固体颗粒作为离散相来完成吸收过程，很少分析水的蒸发过程。这显然与实际的SO2脱除过程存在明显区别。
因此，合理建立脱酸过程的数学模型对准确评估塔内脱酸过程至关重要。目前，关于垃圾焚烧中的脱酸过程数学模型研究较少，且均是简单计算脱酸的传质和反应过程，未考虑气相的速度场、温度场、压力场、气体浓度场和脱酸塔内的喷淋对脱酸效率的影响，这势必影响模拟结果的准确性。
因此，为准确模拟脱酸塔中的气液固流动和化学反应，须充分考虑气固、气液两相间、固体颗粒间以及固体颗粒与壁面间的相互作用。
此外，大多计算流体动力学(CFD)软件的优势主要在于纯气相流动与传热的计算，而对高浓度气固两相流的计算并不十分有效，这主要表现在：未充分考虑气固两相间、固体颗粒间、颗粒与壁面间的相互作用;计算中气固两相的耦合性弱，与真实两相流存在差异;对求解对象中固体颗粒的浓度有限制。
为解决上述问题和满足工业领域气液固三相流的数值计算要求，笔者采用Barracuda软件对某脱酸塔内气液固流动和脱酸反应进行模拟，其特点是采用计算颗粒流体动力学(CPFD)方式求解气液固耦合的流动问题，基于Arrhenius理论计算化学反应速率。
针对脱酸过程，笔者建立了离散相与连续相的化学反应模型和浆滴模型，通过Hill等实验模型验证所建模型的正确性;在此基础上，模拟实际脱酸塔内的脱酸过程，并与实际工况进行对比，进一步验证所建理论模型的可靠性。
1数值模拟方法
Barracuda软件是一款基于MP-PIC(multi-phase particle-in-cell)方式的数值计算方法，能有效解决流体和大量颗粒在三维空间内的耦合计算问题以及在三维空间内耦合求解颗粒和流体动量方程的问题。
在CPFD计算中，所用数值计算颗粒并非物理意义上的颗粒，而是由一定数量、具有相同性质(成分、尺寸、密度和温度等)的颗粒组成的“颗粒团”或“颗粒簇”，由此可将含有数以亿计颗粒的气固系统简化为仅含有数以百万计颗粒的系统，进而对任何固体颗粒实现全尺寸分布的模拟。
CPFD可实现颗粒浓度从很小到紧密堆积的模拟，而无需事先确定颗粒浓度范围，从而得到颗粒的质量、动量、传热和磨损等完整的拉格朗日意义下的信息，为快速准确捕捉脱酸塔中气固液流动的动态特征提供了一种有效可行的手段。
1.1基本假设
针对烟气脱酸过程进行以下基本假设：
(1)连续相与离散相的反应速率受HCl、SO2和Ca(OH)2的扩散速率控制。
(2)随着蒸发进行，浆滴表面的液膜厚度逐渐减小，气相与离散相化学反应速率逐渐增大，液膜消失时化学反应终止。
(3)忽略雾化盘转动对流场的影响。
(4)忽略传质过程中的反应热和脱硫剂的溶解热。
(5)浆滴和反应完全的石灰颗粒均为球形，且两者在整个反应过程中始终保持球形。
(6)浆滴与SO2的反应过程不影响水分蒸发。
1.2连续相(即气相)控制方程
脱酸塔内将烟气视为连续介质，其控制方程包括连续性方程、动量方程、能量方程、组分方程及k-ε湍流方程，通用形式为
2物理模型及边界条件
以某垃圾焚烧厂额定工况为研究条件，在Bar-racuda软件中建立3D脱酸塔物理模型(见图3)。
图3脱酸塔结构示意图
该脱酸塔高21m，截面直径为8.5m。锅炉满负荷运行时，垃圾焚烧处理量为330t/d，烟气量约为77881m3/h;锅炉出口烟气温度为194℃，脱酸塔出口压力为-300Pa，温度为155.5℃，工业水系统用量为110.6L/h，石灰浆供应量为2.3m3/h。表1给出了某电厂实际运行中烟气含量的工业分析。
表1烟气成分分析
基于上文涉及到的化学反应速率方程，编写脱硫和脱酸程序分别嵌入到CPFD软件中计算Hill等搭建的实验模型和该垃圾焚烧厂的脱酸塔模型。
模拟中两模型的网格数目分别约为2.7万和24万，与Fluent生成网格的区别在于CPFD生成的网格为笛卡尔网格，复杂模型最大网格数目控制在100万以内，可较精确地进行模拟。
计算时间步长设为10-3s，此步长需与化学反应速率相匹配，过小将导致计算效率下降，而过大则导致部分网格出现单步长内超温的错误。模拟在8核处理器上单机运行，计算50s的动态过程大约需要CPU运行60h。
3结果与分析
3.1理论模型验证
Hill等建立了喷雾干燥器两相流动中单个液滴的传质传热模型，并搭建如图4(a)所示内径为0.16m、高1.5m的喷雾干燥实验台，需要指出SO2在喷雾干燥吸收过程中的吸收效率取决于吸收过程与干燥过程的叠加，因吸收和反应只发生在水存在的情况下，吸收效率很大程度上取决于蒸发速率与吸收速率之比，且其模拟结果与实验值吻合较好。
为验证所建理论模型的正确性，首先模拟了Hill等搭建的实验装置，结果如图4(b)所示。
由图4(b)可知，模拟结果与实验值平均偏差为3%，标准偏差为1.6%，而王乃华运用双膜理论所得结果平均偏差为7.43%，标准偏差为5.6%，由此可见根据本文所建模型计算所得结果更接近实验值，保证了模型的准确性。
图4物理模型及SO2脱除率
由图4(b)还可知，SO2脱除率随n(Ca)/n(S)的增加而提高，且当n(Ca)/n(S)>2.0时SO2脱除率升高速度减缓。
这是因为随着吸收剂增加，相当于脱硫反应方程式中一种反应物浓度增加，从而加快了脱硫反应;而当n(Ca)/n(S)增至一定数值后，继续增加此值，SO2脱除率的升高速度将变缓，这主要是因为酸性气体的扩散速率达到一定时不再增大，再增加n(Ca)/n(S)对SO2脱除率的影响较小。
3.2与实际电厂运行参数的对比
为进一步评估上述模型对脱酸塔内脱酸过程模拟的准确性，表2给出了模拟结果与实际运行数据的对比。
表2模拟结果与实际运行数据的对比
由表2可知，脱酸塔出入口压力、温度和速度等模拟结果与实际运行数据吻合良好，如脱酸塔入口烟气温度相差0.05K，入口压力相差2Pa，出口烟气温度相差4.35K，SO2和HCl脱除率平均偏差小于3%，由此说明了本文所建模型在合理评估脱酸塔内浆滴的蒸发、SO2和HCl的脱除等过程或反应方面的可靠性。
为阐述脱酸塔内的内流特征对脱酸过程的影响，图5～图9给出了脱酸塔内的流场、SO2和HCl脱除率及蒸发速度。图5分别为Y=0截面和垂直于Z轴不同截面上的压力分布。
由图5可知，受脱酸塔入口段至塔体的扩容效果影响，烟气在进入塔体后压力有所降低，但分布较均匀;塔体下部的锥形结构相对于塔体中部逐渐变窄，压力又逐渐升高。
图5脱酸塔压力分布
由图6可知，在脱酸塔中轴线附近的温度较高，尤其是中轴线上部，这主要是由于烟气进入塔体后主要沿Z轴负方向流动，且未与浆滴进行换热，由于在脱酸塔上部雾化盘喷洒出的浆滴温度较低，当烟气与浆滴发生化学反应并进行换热后，烟气温度逐渐降低，并趋于平衡。
脱酸塔壁面附近温度梯度较小，但右侧壁面局部温度较低，致使当地热应力增大，容易产生热变形，不利于脱酸塔长期稳定运行。
图6脱酸塔温度分布
模拟中忽略雾化盘转动影响，将其简化为液滴沿雾化盘出口切向喷出，雾化盘转速为10000rmin，对应出口速度约为113m/s。
图7(a)为脱酸塔含湿颗粒的质量分数分布。烟气与浆滴均为切向进入脱酸塔，且二者沿圆周运动方向相反，使得大部分浆滴由烟气携带向下运行，一小部分未蒸发完全的浆滴被吹向脱酸塔壁面。图7(b)中颗粒被吹向脱酸塔壁面位置附近，也就是图7(a)标记的截面处的速度矢量图。
由图7(b)可知，右侧存在顺时针方向的回流区，且雾化盘雾化的浆滴进入塔体后初始速度垂直于Z轴，由于这2个原因的共同影响，浆滴被烟气吹向了脱酸塔壁面，增大了脱酸塔壁面结渣和腐蚀的可能性。
图7烟气流场速度矢量和含湿颗粒质量分数分布
图8 Y=0截面SO2和HCl质量分数分布
图8为SO2和HCl在Y=0截面上的质量分数分布。结合图6和图8得出，受入口结构作用，烟气进入塔体后，首先沿脱酸塔中轴线逆时针旋转向下流动，然后又与上部雾化盘产生的顺时针方向的喷雾发生传质传热过程，导致脱酸塔中轴线上部温度梯度和酸性气体质量分数梯度较大。
由图8可知，塔体上部被吸收的酸性气体占SO2和HCl总脱除量的60%和90%，是脱酸反应主要发生部位。图8(a)右侧SO2质量分数低是因为一部分浆滴受烟气作用向右偏斜。
由图8(b)可知，HCl与浆滴的化学反应速率比SO2与浆滴的化学反应速率快很多，在塔体上部HCl几乎已经全部被反应殆尽。图9给出了计算区间内主要特征参数随时间的变化，图中正负值分别代表流入塔内和塔内流出。
图9主要特征参数随时间的变化
由图9可知，40s后曲线波动较小，因此选取40～50s这一段时间作为参考时间段进行分析，上述结果及平均值是基于该段时间内相应物理量的时均值。
由图9(a)可知，由于浆滴的加入出口烟气质量流量相对于入口烟气质量流量略微有所增大。由图9(b)可知，出口HCl质量流量有略微浮动，与入口HCl质量流量相差很大，这是因为约97%的HCl被石灰浆液吸收，仅3%的HCl从出口流出。
与图9(a)、图9(b)和图9(d)相比，图9(c)中出口SO2质量流量随时间的变化最为明显，其主要原因是浆滴与HCl反应较快，且SO2与浆滴的化学反应速率远小于HCl与浆滴的化学反应速率，因此在相同条件下，可认为浆滴先与HCl气体反应，再与SO2反应。
另外，因脱酸塔内温度较高，在浆滴还未与SO2反应完全时就已蒸干失去活性，因此SO2脱除率相对较低。
由图9(d)可知，出口水蒸气质量流量在前5s内小幅增大，其原因在于浆滴中液态水的蒸发使气体中水蒸气质量流量增大。由于浆滴蒸发时间是衡量化学反应时间的一个标准，因此浆滴发生脱酸反应的时间主要在其进入脱酸塔5s内。
表3脱酸实验、实际运行工况与模拟结果脱除率的对比
表3对比了相关文献所做脱酸实验与某垃圾焚烧厂实际运行及本文模拟所得脱酸结果。由表3可知，受实验装置和反应条件等多种因素影响，所得结果差距显著，但本文所得模拟结果与实际运行结果吻合良好。
其中臧仁德等、解海卫等在入口烟气温度为220℃和200℃，李午申在烟气温度为150℃，n(Ca)/n(S)为0.5～3.5的实验条件下，所得SO2和HCl脱除率分别约为50%和90%，而在该电厂实际运行条件下，n(Ca)/n(S)为12时，SO2和HCl脱除率才达到55.63%和56.68%。其原因可能是在实验状况下，模型较小，烟气流速较低，并且实验为理想状态。
Liu搭建的实验台较简单，200℃时SO2脱除率最低仅为23%;而Kiil等在搭建的脱酸实验平台测试的SO2和HCl脱除率最高为67%和91%，虽然与脱酸塔实际运行时的脱酸效率较为接近，但其入口烟气温度仅为50℃，降低温度将显著缩短浆滴的蒸发时间，加快反应速率。
因此，该实验结果并不能准确反映出脱酸塔内真实的脱酸过程。另外，受多种因素和条件的制约，所搭建的实验台仅为实验室规模的缩小模型，其内部流场、化学反应时间和浆滴蒸发时间都有可能与实际装置存在显著区别，因而导致所得结果存在显著区别。
而本文研究对象为原型实际装置，通过合理考虑脱酸塔实际运行中的速度场、温度场、压力场和质量分数场，所得模拟结果能合理反映其实际脱酸过程，可为脱酸塔的运行及优化提供一种可以借鉴的理论模型。
4结论
(1)针对垃圾焚烧厂脱酸塔的脱酸过程，通过求解连续相和离散相的控制方程，充分考虑脱酸塔内速度场、温度场、压力场和酸性气体质量分数场对脱酸过程的影响，建立了脱酸塔气固液之间的传质模型及化学反应模型，并在充分考虑其蒸发及传质的基础上对旋转产生的喷雾建立了浆滴模型。
通过与实验装置和现场运行数据对比，验证了本文提出的气液固的传质模型、浆滴模型和化学反应模型的正确性。
(2)脱酸塔中轴线上部温度梯度和酸性气体质量分数梯度较大，是脱酸反应主要发生部位;一小部分没有蒸发完全的浆滴由于烟气回流区作用被吹向脱酸塔壁面，容易造成脱酸塔壁面的结渣和腐蚀，不利于脱酸塔长期稳定运行。
烟气脱酸主要发生在脱酸塔的上部，脱酸塔中轴线附近温度较高，温度梯度、SO2和HCl质量分数梯度较大，且总脱除量60%的SO2和90%的HCl在上部被吸收，并且HCl与浆滴的化学反应速率远大于SO2与浆滴的化学反应速率;浆滴进入脱酸塔5s内被全部蒸干，即浆滴发生脱酸反应时间主要在5s内。