420 t/h四角切圆锅炉低氮燃烧改造数值模拟研究
0 引 言
随着国家对于环保要求的提高,燃煤电站锅炉的超低排放已普遍实施。其中超低排放明确规定,氮氧化物排放浓度不得高于50 mg/Nm3[1]。为达到严格的超低排放标准,目前国内绝大部分电站锅炉均实施了NOx排放控制技术改造。目前主流氮氧化物控制技术主要采用炉内低NOx燃烧技术和烟气脱硝技术(SCR)相结合的方案。采用低NOx燃烧技术可经济有效地降低炉膛出口NOx排放量,从而减轻SCR运行负担。数值模拟方法是研究煤粉锅炉燃烧特性的经济高效手段[2-3]。Choi等[4]以一台500 MWe四角切圆煤粉炉为对象,利用Fluent软件研究了炉内的流场、颗粒运动、组分和温度分布,并重点考察了NOx的生成,发现空气分级技术能有效抑制NOx的生成。Zhou等[5]运用CFD技术研究了一台1 000 MWe双切圆煤粉炉的燃烧和NOx排放特性,指出燃尽风多组布置比单组布置降氮效果更好。张晓辉等[6]对一台燃用烟煤的220 t/h四角切圆锅炉进行数值模拟,发现水平浓淡燃烧结合燃尽风方式可显著降低NOx排放水平。本文针对一台燃用烟煤的420 t/h四角切圆锅炉,采用水平浓淡燃烧结合空气分级的方式进行改造,并运用数值模拟的方法,重点研究了一次风浓淡比、SOFA风率和SOFA风射流角度等参数对锅炉燃烧和NOx排放的影响规律,从而确定最佳的低氮燃烧改造方案。
1 研究对象
本文针对一台420 t/h四角切圆锅炉进行了低氮燃烧改造,锅炉的炉膛尺寸及燃烧器布置如图1所示。锅炉的原始燃烧器分3层布置在前墙和后墙角区,为直流双通道风包粉结构,一、二次风由双通道燃烧器喷入炉膛。假设切圆直径为800 mm,炉膛截面尺寸为8 650 mm×11 940 mm,燃烧室高度为36 450 mm。改造前锅炉未采用空气分级燃烧,存在燃烧稳定性差和NOx排放浓度高等问题。锅炉燃用煤种的煤质特性见表1。
表1 煤质特性分析
Table 1 Ultimate and proximate analyses of coal
工业分析/%元素分析/%MarAarVarFCarCarHarOarNarSt,arQnet,ar/(MJ·kg-1)13.5023.7323.5939.1849.542.789.030.640.7818.79
改造采用百叶窗式水平浓淡燃烧器替换原直流双通道风包煤燃烧器,如图2所示。该燃烧器通过浓缩器将一次风粉分离成浓淡2股气流,分别从钝体两侧喷入炉膛,其中浓粉气流为富燃料燃烧,造成该区域缺氧,能有效抑制NOx的生成[6];淡粉气流的过量空气系数较高,但由于温度较低,因而NOx生成可得到部分抑制[7-9]。因此,低氮燃烧器在煤粉着火初期便限制了NOx的生成[10-11]。改造后一次风切圆直径由800 mm改为大切圆702.5 mm和小切圆662.5 mm。
图1 锅炉结构及喷口布置
Fig.1 Boiler structure and nozzle arrangement
图2 百叶窗式浓淡燃烧器结构示意
Fig.2 Structure of shutter type rich-lean burner
本改造方案将原主燃烧区部分二次风布置在主燃区上方5.8 m处的3层SOFA喷口送入炉膛(图1),炉膛被分成贫氧的主燃区、还原区和富氧的燃尽区3个燃烧区域。空气分级燃烧是一种可有效降低NOx排放的技术[12-13]。主燃烧区煤粉处在还原性气氛下,减少NOx生成,在燃尽区未燃尽的碳和CO等燃烧完全。
2 数值计算模型与方法
结合研究对象的自身特点,数值模拟过程中所选择的计算模型为:湍流模型采用带旋流修正的k-ε模型,煤粉颗粒运动采用Langrange随机轨道模型来模拟,气相燃烧采用混合分数/概率密度函数模型,采取两部竞争反应模型作为挥发分的析出模型,采用动力/扩散控制反应速率模型模拟焦炭燃烧,采用P-1模型计算辐射换热,热力型NOx生成采用Zeldovich机理[14],燃料型NOx采用总体反应速率模型[15-16]。
3 计算结果与分析
3.1 计算模型的验证
首先对数值计算结果的可靠性进行验证。在改前100% ECR负荷工况下,对锅炉进行数值模拟,将模拟结果和热态试验测量结果进行对比(表2)。发现折焰角处的烟温、NOx排放值和飞灰含碳量虽存在一定的偏差,但尚在可接受范围内,测量结果与模拟结果基本吻合。表明模拟采用模型的可靠性。
表2 测量结果和模拟结果对比
Table 2 Comparison between the test results and simulated results
项目分隔屏进口烟温/K出口NOx含量(6% O2)/(mg·Nm-3)飞灰含碳量/%模拟结果1 3497132.02测量结果1 3457001.81
3.2 一次风浓淡比对低氮燃烧的影响
为了保证单一变量,计算在空气分级燃烧条件下进行,其SOFA风率设为30%。
采用水平浓淡燃烧器后,3种不同浓淡比下炉膛出口温度、氧含量、NOx含量以及飞灰含碳量见表3。可知随着浓淡比的增加,炉膛出口温度逐渐增加,而NOx含量则逐渐降低。浓淡比为4∶1时,飞灰含碳量最低。为了进一步解释上述结果,对3种不同浓淡比下各参数沿炉高的变化规律进行研究,结果如图3所示。
表3 不同浓淡比下模拟结果(折焰角)
Table 3 Simulation results at different rich-lean ratios(furnace arch)
浓淡比温度/KO2浓度/%NOx含量(6% O2)/(mg·Nm-3)飞灰含碳量/%6∶11 3393.772553.454∶11 3333.762713.202∶11 3273.762863.42
由图3可知,3个工况下各参数整体分布规律相似,差别主要集中在主燃区。当一次风浓淡比从2∶1逐渐增大到6∶1时,燃烧器区域截面平均温度呈现先上升后下降的趋势。随着浓淡比增大,一方面浓侧气流的挥发分浓度增加,着火温度降低,着火热有降低趋势;另一方面浓侧煤粉量的增加也导致着火热有增大趋势,上述相反作用使得在某一浓淡比时,煤粉气流的着火热需求最小,着火距离最短,煤粉燃尽率最高,从而炉膛温度最高[17]。O2浓度则随浓淡比的增大呈现先降低后升高的相反趋势,这是由于着火及燃烧性能越好,氧消耗量越大,其中氧浓度场与温度场的变化规律相一致。主燃区CO浓度随浓淡比增大而升高,浓淡比增大,浓粉气流过量空气系数降低,煤粉着火后生成的CO转化为CO2的量减少;虽然淡侧CO生成量有所减少,但其减少量远低于浓侧CO的增加量,故CO浓度呈现升高趋势。随着浓淡比增大,CO浓度升高,增强了主燃区域的还原性,抑制挥发分含氮中间产物氧化成NO,从一方面解释了NOx浓度随浓淡比的增大而显著降低的趋势。另一方面,浓淡比增大使浓煤粉气流挥发分析出速率加快,强化挥发分含氮中间产物HCN和NH3将已生成的NO还原为N2;同时,淡侧气流煤粉浓度低,含氮基团析出量变小,与氧反应生成NO的量减少。在主燃区与燃尽区之间,温度迅速升高并达到最高值;O2浓度迅速降低并达到最低值,而CO浓度则小幅上升。在该区域内高温低氧气氛下,主燃区生成的NOx被大量还原为N2,因而使NOx浓度大幅下降。在燃尽区,随着煤粉与二次风、SOFA风的混合,未燃尽碳、CO与剩余氧反应,3个工况的温度值、O2浓度值、CO浓度值趋于一致,NOx浓度值维持随浓淡比增加而降低的趋势。
图3 不同浓淡比下各参数沿炉膛高度分布
Fig.3 Various parameters along the height of furnace at different rich-lean ratios
综上分析,浓淡比为4∶1时,煤粉的着火及燃烧强度最优,燃尽率最高,NOx排放量也在可接受范围内。因此,选择浓淡比为4∶1的水平浓淡燃烧器作为改造燃烧器,后文在此基础上进行进一步研究。
3.3 SOFA风率对低氮燃烧的影响
3种不同SOFA风率下炉膛出口温度、氧含量、NOx排放量以及飞灰含碳量的变化见表4。可知,炉膛出口烟温、飞灰含碳量随SOFA风率的增加而增加,20% SOFA风率时,NOx浓度较高,SOFA风率由30%增加到40%时,NOx浓度基本保持不变。为了进一步解释上述结果,对3种不同SOFA风率下各参数沿炉高的变化规律进行研究。
表4 不同SOFA风率下模拟结果(折焰角)
Table 4 Simulation results at different SOFA ratios(furnace arch)
SOFA风率/%温度/KO2浓度/%NOx含量(6% O2)/(mg·Nm-3)飞灰含碳量/%401 3983.852713.95301 3333.762713.20201 2913.803042.17
不同SOFA风率下各参数沿炉膛高度分布如图4所示。可知在主燃烧区域随着一、二次风交替喷入,温度呈现锯齿状分布,在主燃区与燃尽区之间达到最高值;随着煤粉的燃尽,炉内温度逐渐降低。主燃区温度与SOFA风率呈非线性相关关系,这是因为主燃区温度受煤粉燃烧放热和二次风吸热双重影响,存在某一SOFA风率下主燃区温度最高。而随着SOFA风率增加,煤粉着火燃烧过程明显延迟,主燃烧区强度降低,未燃烬的煤粉进入燃尽区进行燃烧放热,促使燃烧火焰中心上移,炉膛出口烟温上升。
由图4可知,随着SOFA风率增加,主燃区的O2浓度降低,而CO浓度升高。这是因为,主燃区过量空气系数降低,煤粉缺氧燃烧显著加剧,CO浓度大幅升高;随着燃尽风的喷入,CO浓度急剧降低,且降低幅度随SOFA风率的增加而增大;在炉膛出口附近,SOFA风率增加,煤粉颗粒停留时间缩短,化学不完全燃烧损失增加,CO浓度升高。
图4 不同SOFA风率下各参数沿炉膛高度分布
Fig.4 Various parameters along the height of furnace at different SOFA ratios
主燃区NOx浓度值随SOFA风率的增加而降低,低O2高CO浓度的强还原性气氛抑制了HCN及NH3被氧化成NO,相反促进了其与已生成的NO发生反应生成N2。同时,主燃区产生的CO和焦炭也还原了部分NOx。在高SOFA风风率下,高温区缩小,生成的热力型NOx减少[18]。在燃尽区,NOx浓度有所升高,这是由于煤粉继续燃烧,会重新生成部分NOx,且烟气中的还原性物质基本燃尽,对NOx的还原作用减弱[19]。SOFA风率越高,在燃尽区燃烧的煤粉越多,NOx生成量也越多。在炉膛出口处,20% SOFA风率工况的NOx浓度较高,而30%和40% SOFA风率的NOx浓度相近。
综合考虑炉膛的经济性和降低NOx排放效果,SOFA风率定为30%较为合理。
3.4 SOFA射流角度对低氮燃烧的影响
为便于叙述,定义SOFA喷嘴上扬角度为正角度,下倾角度为负角度,水平方向为0°,分别对-30°、-15°、0°、+15°和+30°五种工况进行对比计算。
不同SOFA射流角度下炉膛出口温度、氧含量、NOx以及飞灰含碳量的变化见表5。可知随着SOFA风射流角度上扬,还原区加长,有利于降低NOx浓度,但燃尽区的火焰中心会上升,煤粉燃尽时间变短,炉膛出口温度和飞灰含碳量上升。对5种SOFA射流角度下各参数沿炉高的变化规律进行深入分析,如图5所示。
表5 不同SOFA风射流角度下模拟结果(折焰角)
Table 5 Simulation results at different SOFA injection angles(furnace arch)
SOFA喷嘴角度/°温度/KO2浓度/%NOx含量(6% O2)/(mg·Nm-3)飞灰含碳量/%-301 3243.973192.35-151 3363.992832.9001 3333.762713.20+151 3483.662503.43+301 3613.572423.97
由图5可知,主燃烧区,除射流角度下倾30°外,其他工况下温度较接近,这表明在一定范围内,改变射流角度对燃烧器区域燃烧影响很小,主要影响煤粉在燃烧器至炉膛出口段的燃烧过程。当射流角度下倾30°时,燃尽区下方温度均较低,因为下倾角度过大,部分燃尽风起到了提前冷却的效果。射流角度在-15°~30°时,主燃烧区至燃尽区之间的温度随射流角度增加而降低,因为SOFA下倾具有“压火”和及时补充氧气的作用,可增加煤粉在此区域的停留时间和燃烧放热。在燃尽区,炉膛温度先降低再增加,燃尽风首先起到冷却作用,然后与未燃尽煤粉混合燃烧放热从而使温度再升高;炉膛出口烟温随射流角度上扬而升高,因为射流角度增加延迟了未燃尽煤粉燃烧,使得高温区上移。
图5 不同SOFA风射流角度下各参数沿炉膛高度分布
Fig.5 Various parameters along the height of furnace at different SOFA injection angles
燃尽风喷入炉膛后,O2浓度先急剧增加而后随着煤粉的燃尽而下降,CO浓度则大幅下降。随射流角度增加,O2浓度降低而CO浓度升高,这是因为,射流角度增大延迟了煤粉燃尽过程,增加了化学不完全燃烧损失;这种低氧高CO的强还原性气氛大大抑制了NOx生成。增大射流角度相当于加长了还原区长度,生成的NOx能更充分地被还原。
实际运行时可根据此规律调整SOFA风射流角度,以达到比较好的降氮效果、燃烧效率和汽温等参数。
3.5 低氮燃烧改造效果分析
根据前文研究,确定所研究锅炉的低氮燃烧改造方案为:将原燃烧器改为浓淡比为4∶1左右的水平浓淡燃烧器;二次风分4层均等布置,与一次风同向切圆,二次风风率为36.5%;SOFA风率为30%,分3层布置,SOFA射流角度上扬15°。热态试验结果见表6。
表6 热态试验结果
Table 6 Results of thermal measurement
项目试验1试验2试验3试验4试验5锅炉蒸发量/(t·h-1)214298373391410主蒸汽压力/MPa13.0513.2813.1913.8213.61主蒸汽温度/℃539542544544541NOx含量/(mg·Nm-3)337284222225256飞灰含碳量/%2.923.073.813.163.92大渣含碳量/%1.350.970.740.485.23
由表6可知,改造后锅炉在低氮运行工况下,能较好地降低NOx排放,保证煤粉燃尽率。负荷360~410 t/h时整体运行状况良好,主蒸汽温度在540~550 ℃,NOx排放量由改造前的750~800 mg/ m3降至220~250 mg/ m3,降幅达65%~70%,飞灰含碳量在3%~4%。中低负荷下,为保证炉内空气动力场和稳定燃烧,通常运行氧体积分数较高,主燃区处于富氧气氛,促进了燃烧初期NOx的生成,导致炉膛出口NOx浓度高于高负荷运行时;通过运行调整,低负荷下其他运行参数变化不大。
4 结 论
1)随浓淡比增加,水平浓淡燃烧器分级程度加大,NOx浓度有所降低,煤粉燃烧强度及飞灰含碳量均在浓淡比为4∶1时最优。
2)随着SOFA风率的增加,主燃区氧含量降低,煤粉燃烧强度减弱,火焰中心随之上移,在炉膛出口处温度、CO浓度、飞灰含碳量增加。在SOFA风率为30%时综合性能最优。
3)随SOFA射流角度上扬,还原区域变长,燃尽区变短,炉膛出口NOx浓度降低,温度、O2浓度、CO浓度及飞灰含碳量增加。
4)热态试验结果表明,改造后各工况下锅炉燃烧状况良好,运行参数稳定,NOx排放达标,燃烧经济性好,表明该改造方案可达到良好的低氮燃烧效果。
[1] 韩文科. 煤电超低排放:机遇与挑战[J]. 环境保护,2016,44(8):39-41.
HAN Wenke .Coal-fired power plants ultra-low emissions-opportunities and challenges[J]. Environmental Protection,2016,44(8):39-41.
[2] KORYTNYI Efim,SAVELIEV Roman,PERELMAN Miron,et al. Computational fluid dynamic simulations of coal-fired utility boilers:An engineering tool[J]. Fuel,2009,88(1):9-18.
[3] 钟礼今,张殿平,田登峰,等. 700 MW四角切圆锅炉低NOx燃烧的数值模拟[J]. 热能动力工程,2015,30(1):118-123.
ZHONG Lijin,ZHANG Dianping,TIAN Dengfeng,et al. Numerical simulation of the low NOx combustion of a 700 MW tangentially-fired boiler[J]. Journal of Engineering for Thermal Energy and Power,2015,30(1):118-123.
[4] CHOI Choeng Ryul,KIM Chang Nyung. Numerical investigation on the flow,combustion and NOx emission characteristics in a 500 MWe tangentially fired pulverized-coal boiler[J]. Fuel,2009,88(9):1720-1731.
[5] ZHOU Hao,MO Guiyuan,SI Dongbo,et al. Numericalsimulation of the NOx emissions in a 1 000 MW tangentially fired pulverized-coal boiler:Influence of the multi-group arrangement of the separated over fire air[J]. Energy & Fuels,2011,25(5):2004-2012.
[6] 张晓辉,孙锐,孙绍增,等. 燃尽风与水平浓淡燃烧联用对NOx生成的影响[J]. 中国电机工程学报,2007,27(29):56-61.
ZHANG Xiaohui,SUN Rui,SUN Shaozeng,et al. Effects ofcombined application of OFA and horizontal bias burner on NOx emission characteristics[J]. Proceedings of the CSEE,2007,27(29):56-61.
[7] 毕明树,张丛,周一卉. 煤粉浓淡空气分级燃烧过程的数值模拟[J]. 动力工程学报,2010,30(8):612-616.
BI Mingshu,ZHANG Cong,ZHOU Yihui . Numerical simulation on processes of pulverized coal bias air-staged combustion[J]. Journal of Chinese Society of Power Engineering,2010,30(8):612-616.
[8] 秦裕琨,李争起,孙锐,等. 风包粉煤粉燃烧原理及实验研究[J]. 中国电机工程学报,2000,20(5):60-63.
QIN Yukun,LI Zhengqi,SUN Rui,et al. Study on experiment and the principle air-surrounding-fuel(ASF)pulverized coal combustion[J]. Proceedings of the CSEE,2000,20(5):60-63.
[9] 孙锐,李争起,孙绍增,等. 四角切圆锅炉炉内煤粉燃烧过程数值模拟[J]. 机械工程学报,2006,42(8):107-113.
SUN Rui,LI Zhengqi,SUN Shaozeng,et al. Numerical simulation on pulverized coal combustion process in a tangential firing furnace[J]. Chinese Journal of Mechanical Engineering,2006,42(8):107-113.
[10] 张抒明. 水平浓淡和竖直浓淡燃烧器的综合性能比较[J]. 华北电力技术,2002(10):45-47.
ZHANG Shuming . Comparison of synthetic behavior of horizontal type conc-weak conditioning coal-burner with vertical type[J]. North China Electric Power,2002(10):45-47.
[11] 周昊,杜黎龙,王正华,等. 水平浓淡燃烧器出口气固浓淡射流的混合特性研究[J]. 中国电机工程学报,2003,23(7):178-182.
ZHOU Hao,DU Lilong,WANG Zhenghua,et al. The diffusion characteristics of the horizontal fuel rich-lean gas-solid burner flow[J]. Proceedings of the CSEE,2003,23(7):178-182.
[12] ZHANG Xiaohui,ZHOU Jue,SUN Shaozeng,et al. Numerical investigation of low NOx combustion strategies in tangentially-fired coal boilers[J]. Fuel,2015,142:215-221.
[13] BELOEVI Srdjan,TOMANOVI Ivan,CRNOMARKOVI Nenad,et al. Numerical study of pulverized coal-fired utility boiler over a wide range of operating conditions for in-furnace SO2/NOx reduction[J]. Applied Thermal Engineering,2016,94:657-669.
[14] ZELDOVICH Y B. The oxidation of nitrogen in combustion and explosions[J].Acta Physicochimica URSS,1946,21(4):577-628.
[15] FENIMORE C P. Formation of nitric oxide from fuel nitrogen in ethylene flames[J]. Combustion & Flame,1972,19(2):289-296.
[16] SHIN Mi-Soo,PARK Mi-sun,JANG Dong-soon. Numericalstudy on the thermal NOx reduction by addition of moisture in LNG flame[J]. Journal of Korean Society of Environmental Engineers,2014,36(12):837-842.
[17] 安恩科,徐通模. 水平浓淡煤粉燃烧器的研究[J]. 热能动力工程,1999(1):39-41.
AN Enke,XU Tongmo . A study of horizontal dense-dilute pulverized coal burners[J]. Journal of Engineering for Thermal Energy & Power,1999(1):39-41.
[18] 赵振宁,吕伟为,张清峰,等. 燃尽风对600MW对冲燃烧锅炉燃烧特性影响的数值模拟[J]. 电站系统工程,2013,29(6):21-23.
ZHAO Zhenning,LYU Weiwei,ZHANG Qingfeng,et al. Numerical simulation on effects of over fire air on the combustion characteristics in 600 MW swirl-opposed boiler[J]. Power System Engineering,2013,29(6):21-23.
[19] 刘敦禹,秦明,刘辉,等. 深度空气分级条件下炉内氮氧化物生成的数值模拟[J]. 燃烧科学与技术,2011,17(5):469-475.
LIU Dunyu,QIN Ming,LIU Hui,et al. Numerical simulation of NO-production under deep air staged combustion in furnace[J]. Journal of Combustion Science and Technology,2011,17(5):469-475.
Numerical simulation of low NOx burner retrofit of a 420 t/h tangentially lean coal-fired boiler
移动阅读
CHEN Tao,YU Jie,ZHAO Peng,et al.Numerical si mulation of low NOx burner retrofit of a 420 t/h tangentially lean coal-fired boiler[J].Clean Coal Technology,2019,25(5):73-79.