以锅炉的选择性催化还原（Selective Catalytic Reduction，SCR）脱硝系统为研究对象，运用数值模拟的方法分析原喷氨格栅结构下烟气与氨气的混合效果，对其结构和布置形式进行优化调整，发现缩小喷氨圆管上喷氨孔的直径并采用两侧布置大孔径中间布置小孔径的形式，可增强氨气射流的穿透力，NH3摩尔浓度的变异系数Cv最高下降20%，烟气与氨气混合均匀性得到大幅提升。
关键词：喷氨格栅；数值模拟；变异系数；混合均匀性
0 引言
选择性催化还原（Selective Catalytic Reduction，SCR）脱硝是指在催化剂作用下，喷入还原剂氨或尿素，将烟气中的NOx还原为N2和H2O。烟气氨氮分布均匀性被视为SCR脱硝性能评价的一个重要指标，作为SCR脱硝系统结构的一部分，喷氨格栅可促使氨气和烟气在进入SCR反应器前充分混合，喷氨装置设计不合理将直接造成氨氮混合不均匀，进而影响到进入催化剂层的反应。只有烟气与氨具有良好的混合均匀性，才能保证催化剂层达到最佳的催化反应和氮脱除效率。国内外常用喷氨格栅进行多点喷氨，使氨均匀地分布在整个反应器截面上。越来越严的排放标准对SCR反应器内的速度场、浓度场、喷氨格栅喷射三者之间的耦合提出了更高要求，系统均流与混合是脱硝系统运行优化的关键之一。以锅炉的SCR脱硝系统为研究对象，采用数值模拟计算方法，在分析原喷氨格栅结构烟气与氨气的混合效果的基础上，对其结构和布置形式进行优化调整，为脱硝喷氨格栅结构参数的优化设计提供参考。
1 模拟对象与方法
1.1模拟对象的几何结构及边界条件
脱硝还原剂采用氨气，原始 SCR 喷氨格栅主要由氨气风道和烟道组成， 计算区域的几何模型如图1（a）所示 ，氨 与 空 气混合 稀 释 后 经氨气入口 进 入 环形氨气风道，并从喷氨圆管的管壁圆孔喷出；烟气从高温烟气入口自上而下流动， 并在方形段烟道内与氨气混合，最终从烟气出口流出。氨气风道为矩形，布置在烟道周边，两侧与喷氨圆管连通，烟道内共布置 5 根喷氨圆管，烟道内每根喷氨圆管中心线上，均设置有对称布置的喷氨孔， 喷口开孔方向与烟气流向、喷氨圆管中心线垂直。SCR 喷氨格栅模型网格划分如图 1（b）所示，运用 ANSYS MESH 软件对三维几何模型进行网格划分， 采用六面体与四面体混合网格，对喷氨圆管网格进行局部加密，最终的网格数量控制在 100 万左右。
图 1 SCR 喷氨格栅几何模型及网格划分
SCR 喷氨格栅入口参数见表 1， 对部分参数进行了简化处理，如用高温空气代替高温烟气，用纯氨气代替氨气与空气的混合气体， 其他参数保持与实际情况一致。
表 1 SCR 喷氨格栅入口参数
1.2 物理模型
对 SCR 喷氨格栅区域进行流场优化模拟是基于 N-S 流动控制方程的求解。采用标准 k-ε 模型模拟气体湍流流动。采用 Species 物质输运模型模拟 NH3在烟气中的混合与扩散， 但不涉及化学反应。开启能量方程，考虑空气与氨气的换热。本模拟假设烟气为单相气体， 不考虑高温烟气中粉尘对流场的影响，将烟气视为不可压缩流体，且为定常流动；假设高温烟气入口和氨气入口的速度分布均匀。烟道入口采用速度进口边界条件， 烟道出口为 Outflow 边界条件；喷氨入口为速度入口，喷射角度与烟气流动方向垂直；喷氨圆管及其他边界设为绝热壁面条件，采用标准壁面方程，无滑移边界条件。
2 模拟结果与分析
2.1原始SCR喷氨格栅的混合分析
原始 SCR 喷氨格栅共设置有 5 根喷氨圆管，每根圆管管壁上开有圆形喷氨孔，其布置如图 2 所示：喷氨孔水平方向上双侧对称布置，间距均为 20 mm，孔直径为 7 mm，每根喷氨圆管布置 20 个喷氨孔，从喷氨孔出来氨气的流向垂直于烟气流向。
图 2 原始圆管上喷氨孔布置示意
通过建立现有 SCR 喷氨格栅区域的全尺度三维模型， 并利用 Fluent 18.0 进行数值模拟计算，获得了现有 SCR 喷氨格栅烟道内的温度和 NH3质量分数分布。图 3 为原始喷氨格栅的温度分布，喷氨入口截面的温度分布如图 3（a）所示，氨气风道的温度较低，方形段烟道的温度较高，这是由于氨气初始温度为 150 ℃，而高温烟气初始温度为 370 ℃。5根喷氨圆管均出现两侧到中间，温度明显逐渐升高的现象，最高温升达 180 ℃。由于壁面均已设置为绝热，所以排除导热造成管内氨气温度升高，这可能是由于通过喷氨孔部分高温空气混入了喷氨圆管中。烟气出口温度分布如图 3（b）所示，总体上看出口的温度分布并不十分均匀，两侧存在局部低温区。
图 3 原始喷氨格栅温度分布
图 4 为原始喷氨格栅的 NH3质量分数分布，喷氨入口截面的 NH3质量分数分布如图 4（a）所示，氨气风道的 NH3质量分数分布为 1，方形段烟道的为 0。5 根喷氨圆管均出现两侧到中间，NH3质量分数分布逐渐降低的现象。而模拟过程中只有 NH3和空气两种组分，这说明随着 NH3在喷氨圆管中流动，方形烟道中部分空气通过喷氨孔进入到圆管中。烟气出口NH3质量分数分布如图 4（b）所示，总体上看出口NH3的分布并不十分均匀，存在中间浓度低，两侧浓度高的现象。
无论从温度还是 NH3质量分数的分布来看，采用原始的喷氨格栅结构都存在高温烟气与氨气混合均匀性较差的问题， 即烟道出口两侧氨气浓度高，中间浓度低的情况。这可能是由于氨气沿着圆管由两侧向中间流动时，其流量在逐渐减小；且喷氨孔是水平布置，高温空气垂直流动；并最终导致喷氨圆管的中间位置高温空气更容易通过喷氨孔进入圆管， 而氨气则更难从喷氨圆管的喷氨孔流入方形烟道。因此，优化设计时还因考虑在工艺允许的情况下， 进一步缩小圆管中间段喷氨孔的直径。
图 4 原始喷氨格栅 NH3质量分数分布
2.2 优化后 SCR 喷氨格栅的混合分析
对原始 SCR 喷氨格栅进行了优化设计，其结构如图 5 所示。喷氨圆管上喷氨孔还是以 20 mm 等间距布置， 有 D5.5 mm、D5.0 mm、D4.5 mm 及 D4.0 mm4 种规格，具体地，两侧布置大孔径中间布置小孔径，喷氨孔的数量和原始喷氨圆管一样，在水平方向上双侧布置，每根喷氨圆管布置 20 个喷氨孔，从喷氨孔出来氨气的流向垂直于烟气流向。通过数值计算获得了优化后 SCR 喷氨格栅烟道内的温度和 NH3质量分数分布。
图 6 为优化后喷氨格栅的温度分布， 其温度标尺和图 3 原始喷氨格栅的温度标尺保持一致。喷氨入口截面的温度分布如图 6（a）所示，氨气风道的温度较低，方形段烟道的温度较高， 这同样是由于氨气和空气的初始温度不一致。5 根喷氨圆管在方形烟道内温度稍微升高了一点， 最高温升不超过 30 ℃，并未出现 原始结构 两 侧 到 中 间 温 度 明 显 升 高 的 现象。烟气出口温度分布如图 6（b）所示，虽然出口还存在小范围的局部低温区， 但总体上看其温度分布还是比较均匀， 相比较于原始喷氨格栅出口的温度分布，局部低温区范围大大较小，温度均匀性明显提升。
图 7 为优化后喷氨格栅的 NH3质量分数分布，其质量分数标尺和图 4 原始喷氨格栅的质量分数标尺保持一致。喷氨入口截面的 NH3质量分数分布如图 7（a）所示，氨气风道的 NH3质量分数分布为 1，方形段烟道为 0。5 根喷氨圆管在方形烟道中 NH3质量分数均出现了小幅降低， 这说明有少量空气通过喷氨孔进入圆管中。但相较于原始喷氨格栅，混入喷氨圆管的空气大幅减少。烟气出口 NH3质量分数分布如图 7（b）所示，总体上看出口 NH3的分布比较均匀，仅存在小范围的低浓度区。
图 5 优化后圆管上喷氨孔的布置
图 6 优化后喷氨格栅温度分布
图 7 优化后喷氨格栅 NH3质量分数分布
2.3 优化前后 NH3分布均匀性对比分析
为进一步了解喷氨格栅优化前后 NH3的分布均匀性，将对 NH3的摩尔浓度进行定量分析。采用变异系数 Cv这一参数作为衡量浓度均匀性的标准，并将其定义为
图 8 优化前后混合烟道各流通截面 NH3摩尔浓度 Cv值的对比
如图 1（a）所示，沿着混合烟道气流方向由上而下分别截取 x=0.1 m、0.2 m、0.3 m、0.4 m、0.5 m 和0.6 m 6 个截面， 并对其 NH3摩尔浓度的 Cv值进行比较分析。图 8 为优化前后混合烟道各流通截面NH3摩尔浓度 Cv值的对比， 可以看出无论优化前还是优化后，NH3摩尔浓度的变异系数都是随着 x 值增大而减小，说明随着烟气与 NH3的不断掺混，NH3的分布越来均匀；且经过 0.6 m 的混合距离，两种结构下的 NH3变异系数 Cv值均减小一半，均匀性均提高了一倍。然而无论哪个截面，优化后的 Cv值均明显小于优化前，下降幅度在 10%～20%之间，说明仅通过调整喷氨孔径来优化喷氨格栅结构，NH3分布的均匀性就能大幅提高。
3 结语
基于原有的SCR喷氨格栅结构进行模拟分析，发现其布置并不合理，喷氨入口截面和烟气出口均存在中间NH3质量分数较低，而两侧较高的现象，烟气与氨气混合均匀性较差。通过缩小喷氨圆管上喷氨孔的直径并采用两侧布置大孔径中间布置小孔径的形式，增强了氨气射流的穿透力，使烟气与氨气混合均匀性得到大幅提升，并最终确立了较优化的喷氨格栅结构。
