为探究二次风配风对炉排炉中城市固体垃圾焚烧过程的影响，针对某750t/d垃圾焚烧炉排炉，采用数值模拟的方法对炉膛焚烧过程进行热态模拟，就下二次风投、停运，上二次风布置形式和上二次风风速3个因素进行优化分析.模拟结果显示，通过在炉拱下方增加下二次风能对炉膛前、后炉拱形成包覆作用，阻挡高温烟气冲刷，有利于改善炉拱区域的结渣问题;炉膛上二次风对冲布置或适当增大二次风风速(从45m/s增大至65m/s)均能有效促进烟气混合，提高炉膛烟气的充满度，改善温度分布的均匀性;上二次风对冲布置较错列布置能进一步提高烟气停留时间，降低炉膛出口的CO体积分数，从而提高燃烧效率.
关键词：炉排炉;垃圾(MSW)焚烧;计算流体动力学(CFD);数值模拟;二次风;配风优化
国家统计局2014年统计年鉴[1]显示，2013年我国生活垃圾清运量已达17238.6万吨.如何处理日益围城的生活垃圾成为亟待解决的问题.2012年，国务院发布的《“十二五”全国城镇生活垃圾无害化处理设施建设规划》[2]提出：到2015年，全国城镇垃圾焚烧处理设施处理能力达到无害化处理总能力的35%.在国家政策的大力扶持下，生活垃圾焚烧处理已进入市场化全面发展阶段.炉排炉是当前垃圾焚烧的主要型式之一[3].
目前国内已有企业[4-5]通过引进如德国马丁、比利时西格斯、日本田熊、日立等的先进技术，加以吸收创新，实现了焚烧设备的自主化，但其运行的经济性及稳定性仍有待进一步提高.
基于计算流体动力学(computationalfluiddynamics，CFD)的数值模拟方法作为一种低成本、高效率的研究手段，近年来不断被应用于垃圾焚烧炉的设计及优化研究工作中.Kear等[6]对稻草焚烧炉的炉排顶部温度以及组分分布进行了数值计算.Yin等[7-8]采用数值模拟结合试验的方式对燃烧小麦秸秆的88MW炉排炉进行了研究，得到了不同建模方式对计算结果的影响.马晓茜等[9-12]对炉排的燃烧进行分段处理，讨论富氧条件对炉排燃烧的影响，并且采用CFD手段对炉膛燃烧以及炉膛顶部通过选择性非催化还原(ivenon-catalyticreduction，SNCR)方法脱除NOx进行了研究.胡玉梅等[13-14]采用数值手段对炉排炉二次风的作用以及二次风的布置位置进行了研究，通过调整二次风位置来调节燃烧条件，抑制二恶英生成.
上述研究主要以国外引进炉排为研究对象，且单机处理量处于中小水平，国内鲜见关于750t/d大型炉排炉炉内燃烧模拟以及配风优化的公开研究成果.本文以国内某自主研发的750t/d炉排炉为研究对象，采用数值方法对炉排炉内燃烧进行模拟.下文将从下二次风投、停运，上二次风布置形式以及上二次风风速这3个角度对其二次风配风进行优化分析.
1、研究对象及计算模型
研究对象为国内某公司自主研发的750t/d垃圾焚烧炉排炉，炉排长为11.66m，宽为12.56m，运行速度为7.673m/h，即运行周期为90min.炉膛容积为648.47m3，配风分为三级，一次风由炉排下方由高到低分5级灰斗两列配送，总风量为88700Nm3/h，温度为493K.二次风温度为313K，总流量为38000Nm3/h，其中上、下二次风所占质量分数分别为80%、20%.根据炉膛实际结构尺寸，通过GAMBIT建立三维模型如图1所示，炉膛网格均采用六面体网格，二次风入口采用局部加密处理，网格总量为109.91万，网格质量较好.
图1垃圾焚烧炉膛的几何模型
垃圾燃烧过程按水分蒸发、挥发分析出、挥发分燃烧和焦炭燃烧分为多个阶段进行，垃圾燃烧过程的模拟分为两大部分，即炉排上方垃圾固相的燃烧和固相燃烧析出的气相在炉膛内的燃烧.固相燃烧反应采用FLIC软件进行模拟，床层上固相反应由文献[15-17]的运动模型描述.固相控制方程如式(1)～(4)所示，与气相控制方程类似，同时考虑了垃圾床层的移动.
床层模拟计算得到挥发出的气相温度、速度及各组分质量分数作为入口边界条件导入到炉膛气相燃烧模拟计算中.气相燃烧模拟得到的床层辐射温度将作为床层模拟的边界条件再次迭代计算，直至收敛.炉膛的气相燃烧反应将通过商业软件AN-SYS中的FLUENT组件进行模拟，粘性模型采用标准k-ε湍流模型，壁面函数为标准壁面函数，气相燃烧采用有限速率/涡耗散(EDC)模型，化学反应为以下3步：
辐射模型为DO(discreteordinates)模型，采用Simple算法求解压力速度耦合方程，控制方程的离散形式为二阶迎风.炉膛壁面采用绝热边界条件，设置炉排沿程为速度入口，上、下二次风也均为速度入口，出口设为压力出口.
2、计算工况
由于我国垃圾未经分类处理，成分复杂，其飞灰软化温度较煤粉更低，现运行的垃圾炉普遍存在排炉、炉膛、烟道结焦的现象[18]，尤其在炉膛喉部区域，未燃烬颗粒受重力作用在前、后墙的炉拱上沉积，在燃烧高温的辐射下，熔融结焦.大量焦块脱落严重影响机组运行的稳定性.为改善这一现象，考虑在前、后墙炉拱下方投运下二次风，用以加强炉拱区域气流的扰动、减少积灰，同时减少高温区的传热.为考察下二次风的作用，对停运和投运下二次风时炉膛热态燃烧工况分别进行数值模拟，各工况下二次风总风量保持不变.
炉膛上二次风的作用主要提供后期燃烧所需氧量，加强气流扰动，以确保燃料的燃烬，从而提高锅炉的燃烧效率;同时引导高温烟气流向，减少因高温区贴墙而产生的壁面结焦、结渣现象.垃圾炉上二次风的布置形式一般采用对冲或错列布置，如图2所示，对冲布置即前、后墙风口水平位置相对，风口数相同.错列布置即前、后墙风口水平位置交错，后墙风口数减少一个，同时增大风口面积以保证风量不变.
图2前、后墙二次风布置形式：错列、对冲
为比较错列布置和对冲布置的效果，这里分别计算了上二次风风速vH为55m/s时错列布置和对冲布置下的2个工况，分别标记为CL55和DC55.同时为探究上二次风风速的影响，在保证上二次风角度及对冲布置方式不变的情况下，分别计算上二次风风速为45、55、65m/s的3个工况，记为DC45、DC55、DC65.本文计算的工况如表2所示.
表2计算工况的配风设置Tab.2
3、炉膛热态模拟结果及分析
3.1下二次风投运的影响
图3下二次风变工况(停运、投运)的炉膛温度分布
如图3所示为工况1和工况4的中心截面温度分布及炉拱截面(A-A)的温度分布，可以发现二次风投运后，下二次风由炉拱下方贴壁而上，将喉部炉拱包覆起来，从图3(b)炉拱截面温度分布可以看出下二次风气流在炉拱宽度方向上的间隔吹扫作用，二次风的投运使得火焰趋势在喉部稍稍靠近后墙，这是由于后拱距离上二次风口较远，同时紧邻主燃区，氧量补入较前拱更困难，因此后墙侧下二次风的注入迅速满足其燃烧所需氧量，从而使得喉部燃烧略靠后墙.
为进一步考察下二次风对炉拱的保护作用，量化两工况的差异，表3对前、后墙沿炉拱截面各40mm深度的范围(如图3(b)两侧黑线所示区域)进行了参数统计，表中，T为温度统计值，φ为体积分数统计值.投运下二次风后，前、后炉拱的平均温度均有所降低，前拱降温尤其显著.
除高温外，还原性气氛是导致炉膛结焦结渣的另一大因素[19]，灰分中的Fe2O3被还原成FeO，FeO与SiO2，CaO等形成共晶体后，灰分的熔点大大下降，从而易于熔融结焦.
表3下二次风变工况(停运、投运)炉拱参数统计
从表3可以看出，投运后前拱附近的O2体积分数有所提高，CO体积分数下降一半.后拱附近的O2体积分数虽稍有降低，但CO体积分数减少了一个数量级，因此可以推断出投运后炉膛前、后拱附近区域的还原性气氛均得到大幅改善，综合考虑到下二次风对炉拱附近的吹扫作用以及对高温的隔离作用，下二次风的投运将有利于改善炉膛喉部结焦严重的问题.
3.2、上二次风布置形式的影响
工况2和工况4错列、对冲布置的迹线分布对比如图4所示，对比可以发现，对冲布置方式相比于错列布置方式迹线更为弥散饱满.
图4变上二次风布置形式(错列、对冲)迹线分布
气流在中心碰撞，动量抵消，有利于炉膛气流的均匀扩散，不至于出现气流过分贴壁的情况.因此气流速度分布相对均匀.结合图4中DC55迹线分布可见，该工况在下部迹线呈现竖直向上，下部靠前墙侧的回流区较小，在上部由于烟气向出口转向而偏向前墙，这样的流场一方面能更大限度地利用炉膛容积，在同样锅炉热负荷的时候有效降低容积热负荷，减轻结焦结渣的情况;另一方面，相对于错列布置，烟气在顶部更加偏向前墙也有利于SNCR喷枪布置，为后期脱硝提供便利条件.
通过观察对冲布置和错列布置方式下的炉膛宽度中心截面温度温度分布(见图6(a))，可以发现，两工况下炉排前端为水分蒸发区，因此温度偏低，而中间段为挥发分的燃烧区，温度上升，在中间区域达到最高温度，而且上炉膛部分也正位于该位置上方，高温烟气直接竖直向上，不需要沿着炉拱绕行，这样在一定程度上避免了炉拱壁温太高导致的结渣问题.
对比两工况下炉内温度场可以看出，温度分布并没有发生根本上的改变，仅表现出在上部炉膛的温度场分布稍有不同，说明上二次风的布置方式对于炉膛主燃区的燃烧状况影响甚微.上部炉膛的温度场分布稍有不同，对冲布置下，上炉膛温度降低较不明显，这说明炉膛烟气与上二次风混合更均匀，因此燃烧也就更完全.进一步反映了对冲布置提高燃尽率的作用.
图6变上二次风布置形式(错列、对冲)截面温度分布
观察高度y=10m下的高度截面温度分布图(见图6(b))，可以看出，对冲布置的温度场均匀性更优，在炉壁四周没有局部高温区的存在.这是因为对冲布置方式能使得前、后墙气流在中心汇聚，并且抵消部分动量后再向四周扩散，而不会因为气流残余动量导致高温烟气贴壁而形成局部高温区.
进一步统计出口O2体积分数φ(O2)、CO体积分数φ(O2)、烟气停留时间τ及炉膛喉部平均温度T0如表4所示.根据环保部《生活垃圾焚烧污染控制标准(GB18485-2014)》对炉排炉规定的炉膛设计相关参数[20]，在实际运行中，炉膛温度应≥850℃，炉膛出口φ(O2)≥6%，炉膛出口φ(CO)≤8.0×10-5，烟气在炉膛停留时间应≥2s.
表4变上二次风布置形式(错列、对冲)炉膛参数统计
从表4的统计结果来看，2个计算工况的计算结果均满足国标对垃圾焚烧炉的运行要求，从数值计算的角度验证了该炉型设计方案的可行性.分析组分布以及停留时间变化规律可知，对冲布置下炉膛喉部平均温度更低，出口O2体积分数和CO体积分数更低，烟气停留时间更长，这是由于对冲布置下，前、后墙上二次风动量抵消后，在主燃烧区形成了较大的扰动，整体烟速降低，也较大程度降低了喉部的烟温，烟气更好的混合使得新注入的氧气得到更均匀的分配，提高了燃烧效率，因此出口的CO体积分数较低.
3.3、上二次风风速的影响
计算得到工况3、4、5的炉膛迹线分布如图7所示，对比可以发现，上二次风风速越大，靠前墙侧漩涡尺度越小，涡团越饱满，从而可以看出上二次风速度对炉膛上部混合效果影响较大，风速越高，混合效果越好，炉膛充满度越高.这是因为上二次风风速的提高使得上二次风刚性增强，进入炉膛以后能有效地加剧喉部低速气流的扰动，使得炉膛内的气流分布更加均匀饱满.
图7变上二次风风速(45、55、65m/s)迹线分布
结合图8中y=12m高度截面竖直速度分布曲线可以看出，45m/s速度工况时回流区宽度以及烟气下行速度均最大，65m/s时平面整体呈现无回流区，随着速度的增大，平面竖直速度的峰值由后墙往前墙偏移.但当速度增至65m/s时，烟气上行速度出现明显的整体增大趋势.
图8y=12截面处沿x方向的平均竖直速度分布(45、55、65m/s)
如图9所示为3个工况的截面温度分布，从图9(a)宽度中心截面分布对比可以发现，下炉膛的温度分布没有较大区别，火焰形态也基本相同，而上炉膛随着风速的增大，中心因为二次风注入而形成的低温区域越小，温度分布越均匀.图9(b)高度截面的温度分布也呈现出风速越大越均匀的趋势.
图9变上二次风风速(45、55、65m/s)截面温度分布
因此可以得出和3.2节相似的结论，在保证炉膛总风量即过量空气系数不变的情况下，上二次风风速变化对下炉膛主燃区的燃烧没有较大影响，风速变化的影响仍主要体现在对上炉膛烟气的扰动及导流上.
表5变上二次风风速(45、55、65m/s)炉膛参数统计
如表5所示为出口O2体积分数φ(O2)、CO体积分数φ(CO)、烟气停留时间τ以及炉膛喉部平均温度T0的统计结果.从表中可以看出随着上二次风风速提高，烟气流速会明显增大，因此停留时间相应降低，将对后期SNCR脱销效率造成一定程度的不利影响.从出口CO的体积分数可以看出，上二次风风速提高一定程度上有利于烟气的充分混合，促进了燃料后期的燃尽，而风速增大与风速动量抵消两股因素相互制约，形成对喉部温度影响不确定的现象.
4、结论
(1)下二次的投运能对炉膛前、后墙炉拱区域起到较好的包覆作用，降低了炉壁温度，改善炉拱附近的还原性气氛，同时贴壁的扰流有利于吹扫炉拱的积灰，因此有利于改善炉拱结渣严重的问题.
(2)炉膛上二次风对冲布置方式相对于错列布置方式有利于促进炉膛内气流的混合，提高炉膛的烟气充满度，提高烟气停留时间，同时降低炉膛出口CO体积分数，从而改善燃烧效率.
(3)炉膛上二次风风速从45m/s提高到65m/s以后，适当高速的上二次风能能够加剧气流混合，使得炉膛气流分布更加饱满，温度分布更加均匀.但过高的二次风速不利于降低炉膛烟气停留时间，从而影响后期SNCR脱销效果.
上述研究可为此类大型炉排炉的设计及优化提供理论参考，从而促进该焚烧设备技术的自主化和大型化发展，满足其日益增长的市场需求.由于现有研究对象还处于研发阶段，未能得到相应的试验数据，无法与模型计算结果进行对比.通过与现有焚烧炉运行标准进行对照，验证了设计方案的可行性.限于篇幅，本文未能深入讨论上、下二次风入射角度及位置等因素对炉膛燃烧的影响，进一步的优化研究可以考虑从以上角度展开.