烟气再循环对生物质炉排炉燃烧影响的数值模拟

烟气再循环对生物质炉排炉燃烧影响的数值模拟 摘要针对生物质锅炉实际运行过程中常出现水冷壁腐蚀严重、屏式过热器积灰多和NOx排放量高等问题，以一台某电厂额定蒸发量为130t/h的生物质往复式水冷炉排炉为研究对象，提出二次风掺混再循环烟气燃烧的方法，采用计算流体力学（CFD）数值模拟技术对炉内燃烧过程进行热态模拟，旨在为锅炉的实际运行操作提供理论指导。计算结果表明，采用烟气再循环可以增强炉膛上部气流扰动，改善炉内温度分布的均匀性，提高燃尽率，同时降低屏区火焰温度，减轻大屏积灰结渣风险；后墙下二次风掺混再循环烟气后，主燃区形成还原性气氛，温度下降，有效抑制热力型NOx的生成。后墙下二次风掺混30再循环烟气的工况炉内气流均匀饱满，高温烟气分布从炉膛深度中心向前、后墙两侧稳定下降，NOx排放质量浓度相对于无再循环烟气时减少了32.1。 我国拥有丰富的生物质能资源，据统计可知，我国秸秆产量约为6亿t/a。目前，除少量生物质被用于农村家庭燃料或饲料外，绝大多数生物质被露天焚烧、填埋，或直接丢弃在田间地头进行生物降解[1-4]。 生物质发电兼具环境效益与社会效益，利用生物质发电可以在替代化石能源、促进环境保护、带动农民增收等方面带来诸多积极影响[5]。根据国家能源局生物质能发展“十三五”规划[6]可知，到2020年，生物质能基本实现商业化和规模化利用。生物质发电总装机容量达到1.5107kW，年发电量为91010kWh，其中农林生物质直燃发电7106kW，城镇生活垃圾焚烧发电7.5106kW，沼气发电5105kW。预计2020年，生物质能产业新增投资1960亿元。其中，生物质发电新增投资约400亿元。 适用于生物质直燃燃料的常见炉型有以下几种形式循环流化床炉、水冷振动炉排炉、往复式炉排炉等[7]。目前，国内应用最广泛的是流化床和水冷振动炉排生物质直燃炉型。 锅炉排放的氮氧化物（NOx）是造成大气严重污染的产物之一，它可以形成酸雨、酸雾甚至光化学烟雾，给人类的生活健康带来严重的威胁[8]。生物质固体成型燃料燃烧生成的NOx来源于原料中的氮元素以及空气中N2的高温氧化[9]，根据其来源，可以将NOx分为3种类型热力型NOx、快速型NOx和燃料型NO[10]。烟气再循环技术是为降低NOx的生成而使用的低空气比燃烧技术，属于脱硝技术的一种[11]。国内外对于烟气再循环技术在燃气锅炉上的应用研究较多。宋少鹏等[12]研究烟气再循环率对天然气低氮燃烧特性的影响，结果表明，在燃气锅炉上，10的烟气再循环率可以减少65的氮氧化物排放。胡满银等[13]对某电站四角切圆燃烧锅炉有、无烟气再循环时的燃烧过程进行数值模拟，结果表明，有烟气再循环时，炉内平均温度分布和最大温度分布整体明显降低；炉膛出口CO2体积分数有一定幅度的上升，O2体积分数有较大幅度的下降；NOx排放质量浓度有较大幅度的降低。目前，有关烟气再循环在生物质直燃炉上应用的研究较少。 针对生物质炉排锅炉实际运行过程中存在的前后拱及侧墙水冷壁腐蚀严重，出现爆管事故、屏式过热器积灰多和NOx生成量高的问题，本文以一台电厂额定蒸发量为130t/h的光大往复式水冷炉排炉为研究对象，采用数值方法对炉内生物质燃烧进行模拟，提出二次风掺混再循环烟气燃烧的方法。考察不同再循环烟气配风布置下锅炉内温度分布、流场分布、燃尽率和排放的污染物质量浓度发生的改变，为生物质炉排炉的设计及优化提供理论参考。 1研究对象及计算模型 1.1物理模型 研究对象是一台额定蒸发量为130t/h的生物质往复式水冷炉排炉，炉排总长为9.3m，炉排宽度为12.0m，炉排运行速度为18.6m/h，整个炉型自出渣口到锅筒顶部高度为20.69m。配风分为4级，一次风由炉排下方由高到低分6级灰斗配送，6个灰斗长度分别为1。47、1.47、1.47、1.47、1.47、1.95m，配风体积比为18∶22∶30∶20∶5∶5。点火风分3列布置，上、中、下3列喷口下倾角分别为70、80、90；前墙下二次风分2列布置，下倾角分别为15和20；后墙下二次风分2列布置，下倾角分别为40和50；前、后墙上二次风喷口均分2列布置，配风方向为水平方向。 出于控制烟气中NOx考虑，模拟计算中总过量空气系数取1.2，其中一次风的总过量空气系数为0.7。再循环烟气量取总烟气质量的30，由上、下二次风分配至主燃区，以降低主燃区的温度，同时增强炉膛上部扰动，提高燃尽率，减少CO的残留。 设计燃料按照黄杆灰杆7∶3的质量比掺混而成，掺混后的生物质燃料元素分析（折算收到基）和工业分析（收到基）如表1所示。由元素分析得到单位质量燃料所需的理论空气量为2.487m/kg。当过量空气系数取1.2时，每小时满负荷所需空气量为123259m/h，燃烧生成理论烟气量为200843m/h。 根据锅炉实际尺寸，使用Solidworks建立锅炉模型，如图1所示，使用Ansysmeshing划分网格，炉膛网格均采用六面体网格，其中，对二次风入口区域进行局部加密处理。构建的计算网格模型最小尺寸为1.43010-8m，最大尺寸为6.10910-3m，二次风入口采用局部加密处理，网格总数为297.89万。 1.2数学模型 生物质燃料在炉排炉的燃烧过程可以分解为固体可燃颗粒在炉排上床层内的燃烧和固相析出可燃物（挥发分、炭粒）在炉膛内稀相空间的燃烧。 固相燃烧反应采用谢菲尔德大学垃圾焚烧中心自主开发的模拟平台FLIC。床层上固相反应由Yang等[14-16]的运动模型描述。考虑层料床层的移动[17]，固相控制方程如下。 炉膛的气相燃烧反应计算应用通过商业软件ANSYSFLUENT进行模拟，床层固相燃烧计算得到的炉排上方气体的温度、速度及组分信息，作为炉膛气相燃烧模拟的入口边界条件以UDF的形式导入，具体炉排及炉膛计算的耦合方法参见文献[17]。气相燃烧反应的黏性模型采用标准k-ε湍流模型，壁面函数为标准壁面函数，组分输运为通用有限速率模型，气相燃烧源相采用有限速率/涡耗散（EDC）模型计算，化学反应如下 2CH43O2-→2CO4H2O,（5） 2COO2-→2CO2,（6） 2H2O2-→2H2O。（7） 辐射模型为DO模型，采用Simple算法求解压力速度耦合方程，控制方程的离散形式为二阶迎风。炉膛壁面采用绝热边界条件，设置底部炉排及上下层二次风均为速度入口，炉膛顶部为压力出口。开始计算时，求解温度场、流场、动量方程、颗粒场、燃烧和辐射换热，当残差小于设定标准且烟道出口计算参数稳定，即认为燃烧计算收敛。 NOx排放计算采用燃烧的后处理计算方法，在燃烧计算收敛后，根据燃烧数值解开启NOx反应的模拟计算，此时流动模型、湍流模型、能量方程、辐射模型和组分计算模型不再进行迭代求解，NOx反应与燃烧过程计算不耦合。NOx生成过程较复杂，在计算过程中重点关注燃料型NOx和热力型NOx的转化过程（见图2），NOx的具体转化模式见文献[18,19]。 采用的床层和