1、基于 CFD技术的镁法脱硫塔入口烟道流场优化分析要:以某 300MW 机组镁法脱硫塔为研究对象, 针对入口烟道内存在的气流冲壁现象,利用计算流体动力学软件 Ansys13.0,对其进行了流场数值模拟。塔内烟气采用 Reynolds 时均 Navier-Stokes 方程描述,离散相采用拉格朗日模型描述。通过对布置不同导流板的方案分析表明,在烟道内布置若干数量的导流板, 可以有效地缓解入口烟道内壁的气流冲刷现象;同时导流板的布置使得塔内烟气流场分布更加均匀, 气液两相混合更的均匀,并提高了吸收效率。综合考虑入口导流板数量以及塔内流场分布, 布置三块导流板最为合理;另外针对入口处布置不同数量导流板
2、的模拟结果对镁法喷淋塔的现场运行以及结构优化提供理论依据。关键词:镁法烟气脱硫;导流板;优化;Ansys13.0中国是一个以煤炭为主要能源的国家,并且大部分煤炭品位低, 原煤中灰分和硫含量较高1,煤炭燃烧生成的 SO2 已成为中国大气污染的主要污染物。火电厂排放的 SO2 约占全国 SO2 排放量的 1/32。我国在 1995 年 SO2 排放量严重超出环境的自净能力, 其总量达到 2370 万吨,超越欧洲和美国等地区。目前我国发电企业仍以火电为主,控制 SO2 排放成为我国发电企业面临的重大问题。自 2002 年,中国在电力行业内开展了大规模的 SO2 治理工程。为控制以 SO2 排放为主造
3、成酸雨污染的恶化趋势, 脱硫工艺的发展也越来越成熟 ,这其中应用比较广泛的是镁法、钙法、氨法等工艺3-4。镁法 脱硫技术的优点有:占地面积小、综合效益好、运行费用低、安全可靠、脱硫效率高、没有二次污染、副产品能产生一定的经济效益等,同时又避免了大型湿法脱硫的诸多缺点,因此镁法脱硫将逐步得到更为广泛的应用5。多数湿法喷淋塔在设计和运行时均没有充分考虑烟道的优化,喷淋塔内气液两相流场以及进口烟道内的流场分析对于整个系统的可靠性及安全运行具有重要的意义。本文针对某 300MW 机组镁法脱硫塔入口烟道存在的气流冲刷内壁面的现象进行分析,采取喷淋塔入口烟道处布置不同导流板的方案作为应对策略, 运用 An
4、sys13.0 软件,对其导流板的数量以及塔内流场进行了三维数值模拟,其模拟结果与原塔进行对比分析, 选出最优方案, 期望对系统的实际运行及优化提供理论依据。1 喷淋塔内部烟气流场数学模型的建立1.1 烟气相流场的控制方程在脱硫塔内,烟气向上流动,液体向下流动,属于三维非定常气液两相流动,气液两相在塔内的运动过程中伴随着质量和动量方面的交换,这些交换是一直存在的。欧拉坐标系下的烟气流场方程为式中: 为通用变量(可以是湍动能、耗散率和速度等的求解变量); 为广义扩散系数;U 为速度矢量;Sp 为液滴对烟气的作用源项。在此方程的基础上气相的模拟采用标准 k- 双方程模型实现方程的封闭6-7。同时
5、,近壁面选用壁面函数法来修正8。因为 Launder 等9 人提出的 k- 双方程模型有较高的计算精度且计算量较小的优点, 所以雷诺应力项的封闭采用 k- 双方程。1.2 液相控制方程由于 Lagrangian 颗粒轨道模型可以完整地考虑颗粒与流体间的相互作用, 采用随机轨道方法进行修正,并且在计算过程中跟踪计算颗粒沿轨道的动量增加与损失, 与此同时所得计算结果方便用于后续的连续相计算中,因此液相模拟采用 Lagrangian 颗粒随机轨道模型。液滴相控制方程式中:dp 为粒径;u 为速度,m/s; 为动力粘度,Pas; 为密度,kg/m3;d 为粒径,m;CD 为同 ReD 有关的函数。2
6、脱硫塔内烟气流场的模拟计算2.1 建立脱硫塔模型及假设本文模拟的喷淋塔模型如图 1 所示( 简化后模型),其模型为逆流喷淋塔。烟气从脱硫塔入口进入塔内,塔内喷嘴向下喷淋液滴 ,烟气与液滴互相接触发生化学反应进而达到净化烟气的目的,被净化后的烟气从塔内上部的除雾器排出喷淋塔。图 1 喷淋塔结构示意图由于采取以下假设能使建模以及计算过程得到简化10-13,同时结果的精确性受到略微的影响。因为本文主要的研究是对脱硫塔内的流场作定性分析而非定量分析,略微的误差是可以接受的。针对本文的模拟,对烟气与液滴两相作如下假设 :(1)喷淋塔模拟区域为除雾器以下至浆液面部分 ;(2)烟气可看作是不可压缩流体 ;(
7、3)暂不考虑塔内喷嘴等部件对烟气流场的影响 ;(4)脱硫塔内进行的过程复杂多变 ,其内部的流场是非稳态的,本次模拟假定为稳态。2.2 气液两相流场参数及边界条件烟气与液滴的参数及边界条件,如表 1 所示。2.3 网格的无关性讨论本文分别选取了 67 万、82 万、109 万以及 134 万数量的网格分别对四类几何模型进行了模拟,结果发现随着网格数增加 ,相邻网格数间相对误差逐渐减小 ,直至网格数达到 109万左右,再增加网格数目,其模拟结果基本不变,经考虑最终选取 109 万网格方案, 由于 4 种几何模型相差不大,故均采用 109 万网格, 这种做法不但运算简洁,结果的精确性也得到了保证。3
8、 计算结果与讨论通过对在喷淋塔入口布置不同数量的导流板进行模拟,其模拟结果与原塔进行对比分析,选出一个最优的布置方案。顶层喷淋层上方 1m 处的截面的速度云图, 如图 2 所示。图 2 顶层喷淋层上方 1m 处的截面的速度云图从图 2(a)中可以看出,原塔截面的速度不规整 ,中部左侧速度比右侧大 ,说明经过喷淋层的烟气与喷淋液滴接触不充分, 直接影响吸收效率; 而图 2(b)、图 2(c)、图 2(d)中布置导流板后, 速度分布逐步得到了改善 ,图 2(d)中改善最为明显,速度分布趋于均匀。比较分析得出,脱硫塔入口处布置三块导流板的方案中烟气经过喷淋层到达出口时速度比较均匀,说明气液两相在这之
9、前接触的比较充分 ,这种布置方案对于增强浆液与烟气的传质交换,优化气液分布均匀性 ,提高系统 SO2 的吸收率都有很大的意义。不同入口条件下的喷淋塔系统的气流速度云图,如图 3 所示。图 3 不同入口条件下的气流速度云图从图 3(b)中看出布置一块导流板时,入口内壁的气流冲刷现象没有得到改善, 还是存着冲刷现象;图 3(c)中布置两块导流板时较之前的方案内壁面冲刷现象有所缓解;而图 3(d)中布置三块导流板时,冲刷现象得到了明显改善。比较分析得出, 布置三块导流板在很大程度上对流经该区域的气流进行了疏导,减少了气流对该区域的冲刷。图 4 给出了底层喷淋层下方 1m 处截面的速度等值图,图中数值
10、代表速度 ,此处的截面比较具有代表性,由于该截面是烟气与液相的主要反应区域,SO2 吸收的好坏跟该截面有很大关系。因此,烟气在此横截面上均匀性以及气液混合均匀程度的好坏均可由该截面上的速度分布的好坏来反映。另外,从图 4(a)、图 4(b)、图 4(c)中可以看出,截面中间区域主流速度为 1m/s,与塔的设计烟气速度 3m/s-5m/s 相差很大 ,反而越靠近壁面速度越大 ,与设计速度比较接近;而图4(d)中布置三块导流板的截面速度分布比较均匀,整个截面的速度基本在 3m/s-5m/s 左右,这是最有利于脱硫反应的速度。图 4 底层喷淋层下方 1m 处截面的速度等值图通过上述分析得出,在喷淋塔
11、入口处布置三块导流板时 ,在同一截面里 ,烟气速度分布较为均匀,因此有利于气液两相很好的混合 ,增加接触时间,从而提高 SO2 吸收率。对于脱硫系统来说,研究脱硫塔阻力特性是比较有意义的。不同入口条件下吸收塔内压力沿塔高的变化曲线,如图 5 所示。由图 5 可以看出,三种方案的曲线的下降趋势与原塔大致相同,没有出现突然下降和或者突然上升的趋势 ,也说明布置导流板的几种方案都是可行的, 影响压降的主要因素是喷淋层的层数 14,本次为三层喷淋层。考虑到之前讨论的合理性,认为布置三块导流板是最优方案。图 5 不同入口条件下吸收塔内压力沿塔高的变化曲线塔内喷淋的效果图,如图 6 所示。一般情况下, 喷
12、嘴喷出的浆液有很多种形态,此次模拟中,喷淋塔喷嘴为实体锥形,基本上喷淋形成的锥体可以认为是以喷嘴位置为顶点的锥体。这样的喷淋能全面覆盖烟气, 使浆液与烟气能更好的接触,利于 SO2 吸收,从而提高其脱除率。图 6 中空锥形喷嘴喷淋效果图4 结论对比结果可知,在喷淋塔入口烟道内合理地布置导流板 ,能够有效的改善入口烟道以及塔内流场分布状况,有利于两侧的烟气平衡 ,使得气液两相在塔内更好地接触 ,提高 SO2 吸收率;同时,从压力曲线也可以看出, 几种方案的下降趋势和原塔相近,说明几种方案都是可行的,但相比其他方案,入口处布置三块导流板的优化结果更为明显。与此同时,几何模型中也布置了四块甚至更多导流板 ,并对其进行数值模拟 ,结果发现,布置四块或更多导流板时,其模拟结果与三块导流板的结果相差无几, 非常接近;另外,综合考虑成本的问题,认为方案 d(布置三块导流板 )为最优方案。期望本次模拟的结果能对系统的实际运行及优化提供理论依据。东北电力大学学报作者：鞠铠阳 1,李威 2,王欣 2,洪文鹏；1. 天津国投津能发电有限公司,2.东北电力大学能源与动力工程学院
