1. 引言
煤炭在我国能源需求和消耗中占据重要位置,而燃煤所产生的氮氧化物是污染大气的主要有害物质之一 [1],所以燃煤电厂NOx的排放成为了我国烟气污染物治理的重点。目前常用的脱硝技术有低氮燃烧、选择性非催化还原(SNCR)及选择性催化还原(SCR)技术,其中SCR技术因具有受锅炉运行条件影响较小,反应条件易于控制,脱硝效率高(可达90%以上)等优点得到了较广泛的应用 [2] [3]。SCR脱硝系统内部烟气流场不均匀时,易造成还原剂与烟气混合程度低及氨逃逸过量等情况,对脱硝效率产生消极影响。因此,为保证脱硝系统的高效运行,对其内部流场进行优化就显得十分重要。
目前CFD软件在脱硝系统优化中很受欢迎,该软件通过建模、迭代计算及分析结果,可以达到优化设计研究的目的,具有可操作性强,低成本,可缩短产品周期及计算结果直观的优点,是脱硝系统流场优化的重要工具。
本文以煤粉炉为例,从SCR系统流场优化角度,列举了目前SCR系统数值模拟优化的进展。
2. SCR反应概述
选择性催化还原(SCR)技术是在烟气中加入氨和尿素等还原剂,在320℃~420℃的温度窗口,还原剂与烟气中的NOx在催化剂表面发生反应,生产氮气和水。该反应是氧化反应和还原反应相互竞争和作用的结果,反应示意图如图1所示。
主要反应如下:
SCR反应受温度、氨氮浓度等因素影响,浓度低或者混合不充分都会影响脱硝效果及运行的经济性,因此对SCR系统进行流场模拟就显得十分重要。流场模拟主要分为确定求解方程、前处理(建模及绘制网格)、Fluent参数设置,求解及后处理几个部分。
3. SCR数值模拟概述
3.1. 控制方程
根据SCR运行特点,在建模时可进行如下假设和简化:
1) 将烟气视为不可压缩流体;2) 假设进口烟气速度分布均匀;3) 催化剂层及整流装置等多孔结构采用多孔介质模型进行计算。基于以上假设及简化,SCR系统烟气流场的连续性方程、动量方程、能量方程及组分守恒方程的通用形式可表示如下 [4]
(1)
式中:S为广义源项;
为广义扩散系数;
为通用的因变量;
为流体密度;
为时间;
为速度。
Fluent多孔介质模型就是在定义为多孔介质的区域结合了一个根据经验假设为主的流动阻力。本质上是在动量方程上叠加了一个动力源项,源项由两部分组成:粘性损失项和内部损失项。且在主流方向和非主流方向相差不超过1000倍。
(2)
式中:
为流体粘度;
为i方向(
,
或
)的速度分量;
为i方向(
,
或
)动量源项;
为内部阻力因子,
为介质渗透性。
湍流模型采用
模型,其连续性方程为 [5]
(3)
式中:
为流体密度;
为时间;
为位移矢量;
为i组分的速度矢量物质输运模型。
绝大多数SCR流场模拟时,采用无反应的混合组分,通过第i种物质对流扩散方程预估其质量分数
[6]
(4)
式中:
为物质i的扩散通量,由浓度梯度产生;
为化学反应净产生速率;
为离散相和定义的源相导致的额外产生速率。此处
和
均取0。无反应的混合组分,只需考虑i种物质的湍流扩散通量
[7]。
式中:
是混合物中第i种物质的扩散系数。
是湍流施密特数,
。
是流体粘性系数,
是扩散系数。
3.2. 建模及网格划分
流场模拟前处理软件可选择Gambit,Workbench,Icem等。建模时根据设备结构特点对内部装置进行合理简化,可大大提高计算效率,例如喷氨喷嘴可简化为喷口,忽略对流场影响较小的装置(构架、梁等)。
在弯头或者结构较复杂处建议采用非结构性网格。为同时兼顾计算的准确性及高效性,网格数量可通过网格无关性认证来确定。
3.3. Fluent设置
在Fluent设置中,进出口烟道一般设置为速度入口及压力出口,烟气与氨气的混合选择多组分输运模型,催化剂、整流格栅等多孔结构可采用多孔介质设置,孔隙度根据结构参数确定。模拟时,不考虑氨和NOx的化学反应,只考虑二者的混合。
4. SCR内部装置模拟优化
通过数值模拟优化SCR系统的内部装置,例如导流板,喷氨格栅及整流格栅等,可改善SCR脱硝系统内烟气流场,从而提高脱硝效率。
4.1. 导流板的优化
陈东林等人 [8] 对湖南一台300 MW燃煤机组进行导流板优化模拟,发现在喷氨后布置导流板对首层催化剂速度流场产生较大影响,而喷氨前的导流板布置对首层催化剂层氨浓度场分布影响较大,对导流板位置优化提供了一定参考。
王汉青等人 [9] 对某300 MW机组脱硝装置进行模拟,通过分析速度相对偏差等五个参数,发现距离反应器最近处的导流板尺寸对流场影响最大,而且在反应器顶部安装整流格栅,可使流场更均匀。为使某330 MW机组首层催化剂入口速度偏差小于15%,改善系统压降,杨松等 [10] 发现反应器与水平烟道连接处安装直板、两个转弯烟道安装弧形导流板时脱硝装置速度场及浓度场分布较均匀。
大部分学者进行流场模拟时都选择忽略了脱硝系统中的灰尘,为使模拟结果更有现实指导意义,朱天宇等 [11] 对某600 MW燃煤机组脱硝装置进行模拟时加入了飞灰颗粒的离散项,发现该装置流场受飞灰影响不大,但考虑到催化剂堵塞及使用寿命,在灰量较大的SCR脱硝系统中建议安装吹灰装置。还得出了弧形板后连接一段竖直板导流效果更好的结论,该结论与文献 [10] 都提出了弧形板与竖直板的配合,可为导流板的结构优化提供新的设计依据。
于玉真等 [12] 在前人的基础上,以某660 MW燃煤超临界锅炉为研究对象,利用正交实验法,证明了安装导流板并不会增加脱硝系统的压力损失,氨烟气混合器会增加一定的压降,但其整流效果显著。毛剑宏等 [13] 发现添加导流板后脱硝系统的压降减少很多,并将模拟结果与冷态试验数据进行对比分析,进一步验证了模拟结果,具有一定的参考价值。综上,可以发现导流板不仅能使流场更均匀,还能降低系统压降,因此合理添加导流板就显得尤为重要。
4.2. 喷氨格栅及整流格栅的优化
吴学智等 [14] 提出不均匀喷氨法,将喷氨格栅划分成若干个区域,根据各区域的氨浓度分布来调节喷氨格栅的开度,使氨氮混合更充分,增强了反应效率。在实际工况中,大型燃煤锅炉省煤器出口烟气流场呈不均匀分布,但当前大部分学者对SCR系统进行模拟时采用了均匀入口边界,故周英贵 [15] 决定对非均匀入口条件的SCR氨喷射做模拟研究,针对脱硝入口烟气的不均匀性调整喷氨量,更符合工程实际。王乐乐等 [16] 使用数值模拟的方法对某电厂的喷氨格栅结构优化方案进行了验证,得出了该双向调节(横向烟道与纵向烟道)方案具有可行性的结论,该评估方法对火电厂脱硝系统优化具有一定的借鉴作用。
张巍等 [17] 对某600 MW机组整流格栅进行了优化模拟,发现该系统整流格栅的压力损失为100 Pa时,首层催化剂入口速度偏差及浓度偏差最小。陶立新等 [18] 将整流格栅结构优化为阶梯式,发现整流效果增强。徐劲等 [19] 发现,整流格栅及首层催化剂高度都会对首层催化剂入口烟气速度及氨浓度分布造成影响。格栅高度变化时,需相应改变催化剂的安装高度。高建强等 [20] 通过模拟研究确定了某1000 MW燃煤机组SCR脱硝反应器整流格栅的最佳间距以及整流格栅与催化剂安装高度的比值对流场的影响。整流格栅与催化剂安装高度的耦合研究,进一步为脱硝反应器流场优化提供了新的思路和方向。
5. 结论
1) 导流板的安装位置,数量及结构对导流效果存在影响,且距离反应器最近处的导流板尺寸对导流效果影响最大。安装位置一般位于烟道转角处,安装数量可参考相近功率机组的脱硝系统,结构上采用弧形导流板配合竖直直板效果更好。对导流板进行模拟优化发现该装置不仅能够改善流场分布,提高脱硝效率,而且合适的导流优化措施会减小系统压降。
2) 不均匀喷氨技术的应用进一步增强了还原剂与NOx的反应,提高了脱硝效率,降低了氨逃逸率。
3) 对导流板、喷氨格栅及整流格栅的进行优化时,要考虑到各结构之间的耦合优化。在安装导流板的基础上,同时在反应器入口配以整流格栅,整流格栅间距、高度和结构的优化及整流格栅与催化剂安装高度的耦合都会使流场更加均匀。