APP下载

基于FLIC软件的生物质层燃数值模拟

2021-04-08杨海瑞

关键词:风量生物质组分

肖 凡, 胡 南, 张 缦, 金 燕, 杨海瑞

(1. 太原理工大学 电气与动力工程学院, 太原 030024;2.清华大学能源与动力工程系, 北京 100084; 3. 长春工程学院, 长春 130012)

1 引 言

生物质能是指可以进行光合作用的生命体将太阳能以化学能的形式固定在其体内的能量形式,和太阳能一样都属于可再生能源. 从广义范围来讲, 生物质涵盖所有的植物、微生物及将其作为食物的动物以及其生产出的废物. 如城市及工业有机废弃物料、农作物及其废弃物料、木材和森林废弃物、以及动物粪便等都可以称得上是代表性的生物质[1-2].

中国生物质能能源极其丰富,各类秸秆的年产量可达6亿吨,薪柴等产量大约为2亿吨.每年仍有较大量的人畜生产粪便和工业废弃物料等. 目前,除了很少一部分被用作动物口粮或者家用燃料使用外,其余大多数均被露天焚烧或者就地填埋,还有另外一部分丢在田间被生物降解.

随着环境污染和能源短缺等问题的日益突出,人们愈加重视生物质能利用,许多国家都开始行动起来,积极投身于生物质能的开发及利用. 在多种利用方式中,直燃利用是一种比较重要的规模化利用方式. 生物质直燃代替化石能源,不仅能促进节能减排,还能带动农民增加收入,因此合理有效利用生物质能既具环境效益又兼具社会效益,对促进社会经济增长和改良生态环境极具意义[3].

目前,对生物质进行规模燃烧的最主要的方式包括层燃燃烧和流化床燃烧[4]. 其中,比较常见的是层燃燃烧,一般规模较小,但更适合生物质资源的分布规律. 在东北三省,常采用生物质层燃热风炉,利用热空气对粮食进行干燥加工,可较好地实现生物质的综合利用. 生物质层燃可以分为四个过程[5],依次为干燥、热解、燃烧和燃烬阶段. Yang等[6-8]对垃圾层燃进行了全面系统地理论以及试验研究,开发出一个二维模拟平台FLIC(Fluid Dynamic Incinerator Code),并综合考虑了固定燃烧和移动燃烧,该平台目前已经成为研究层燃的最常用工具. 利用 FLIC软件针对床层进行数值计算,不仅可以得到用于炉膛稀相气体燃烧的边界条件,还能得到炉篦上各成分的燃烧状况及其热解的过程,为炉排炉的设计和优化提供了指导.

对于生物质燃料应用,在实验方面,李晓翔等[9]对农林生物质进行热重实验,结果表明, 生物质热解可大致分为干燥、预热、热分解氧化以及碳化还原4个阶段. 在数值计算方面,赵芳芳[10]用Fluent软件对某电厂改造后的48 t/h生物质炉排炉进行研究;余昭胜等[11]采用FLIC和Fluent软件对秸秆直燃锅炉在特定负荷,同时改变一、二次风比,对锅炉燃烧特性进行研究. 与煤相比,生物质燃料挥发分含量高,固定碳含量低,然而目前很多生物质层燃炉设计和运行依旧沿用燃煤的思维, 不仅影响利用效率, 甚至会引发爆燃事故. 本文采用数值模拟方法研究生物质燃料在炉篦上的燃烧过程,其结果可对该类锅炉的结构以及运行提供优化指导.

2 热风炉结构及燃料特性

本研究对象为东北三省某生物质热风炉,额定热功率为7 MW. 该炉结构简图见图1,炉排长度为4.3 m,宽度为2.5 m,一次风从炉排下部进入炉膛,二次风喷口设在前墙和炉顶,风口内径为27 mm,间隔200 mm.

1.再循环风口 2.二次风口 3.隔墙

该炉燃用的是由玉米秸秆成型燃料,该燃料元素分析、工业分析见表1,其低位热值为12 790 kJ/kg.

表1 生物质元素分析和工业分析

3 基于FLIC平台的层燃模拟

3.1 FLIC软件参数设置

FLIC是谢菲尔德大学开发的床层燃烧反应模拟平台,主要应用于生物质、垃圾等高水、高挥发分等燃料的层燃燃烧模拟. 用户给定燃料的工业分析、元素分析、低位发热量、一次风分配、炉排尺寸等边界条件,然后通过求解炉排上方气相和固相的连续性、动量、能量和组分输运方程,即可得到沿炉排运动方向的速度分布、温度分布以及组分分布,并以此做为稀相空间数值计算的边界条件.

该算法基于四阶龙格库塔法[9],假定挥发分气体成分为CmHn(33.2%)、CO(17.54%),H2(10.22%)和CO2(31.06%). 燃料由辐射传热引燃. 计算时燃料粒径为4 cm×4 cm×10 cm,其形状因子为1.35. 堆积密度为 450 kg/m3,其孔隙率为0.719,料层初始高度为250 mm. 一次风量为209.2 m3/min,温度为273 K,四个风箱中风均匀分配. 床层移动速度为8.53 m/h.

3.2 模拟结果

燃烧计算表明:燃料水分蒸发率为98.25%,挥发分析出率100%,焦炭燃烧率为71.65%,燃料总的质量损失是84.37%,床层计算表明其燃尽程度较高.

图2a 炉排上方气体温度分布

图2、图3分别表示炉篦上方气、固体的温度分布. 由于一次风温较低,燃料随炉排进入炉膛区域以后主要依靠炉膛中烟气辐射进行加热. 由图2b可知,在沿炉排前进0.53 m之前,燃料处于干燥阶段,床层厚度变化不明显,温度较低变化也不大. 从0.53 m位置开始,挥发分开始逸出并燃烧,床层上方气温也很快随之上升,水分蒸发加快,燃料此时由于热解所以不断被耗尽,厚度逐渐减薄. 在0.6~1.0 m之间,床层上方气温基本保持不变,与此同时床层高度下降,这是由于该段有大量挥发分析出,而热解也会吸收热量. 1 m之后,前段逸出的挥发分强烈燃烧,使炉篦上燃料温度达到最高,为1 590 K. 气相最高温度出现在沿炉篦前进方向2~2.2 m处,最高可达1 348 K. 2.2 m之后,挥发分放热基本结束故温度开始变低,并且此时固定碳开始燃烧,但是由于固定碳着火温度较高且此处床层较薄,故对一次风升温程度有限且一次风温度较低,所以固定碳只燃烧大概71%便不再燃烧,所以固定碳燃烧燃尽率并不是很高.

图3 固相温度分布Fig.3 The temperature profile of solid phase

图4所示为燃料主要的组分变化趋势. 水分变化见图4(a),上部燃料受到辐射传热,水分从燃料中脱出,并呈抛物线状沿高度方向呈现.图4(b)为挥发分变化,与图3对比发现,当其表面温度高达600 K时,伴随着挥发分逸处之后燃料体积缩小,床层厚度明显变薄. 固定碳含量变化如图4(c)所示,由于其着火温度较高,在刚燃烧初始阶段其含量保持不变,受到辐射及自身燃烧放热加热后开始燃烧. 固定碳含量在两个阶段会发生突变: 一是挥发分逸出阶段, 此时其含量从最初的13.23%增大至接近54.1%; 二是当固定碳氧化燃烧时,在朝着炉篦前进方向大概1~2.7 m 处,其含量从54.1%降到25.6%以下.

图4 燃料组分含量变化

图5、图6分别为沿炉篦前进方向燃料各组分的变化率和烟气各组分质量分数变化. 在朝炉篦前进方向0~1.9 m 处水分脱出. 从0.5 m 处开始,温度持续上升,到600 K挥发分开始逸出,当析出的可燃气体氧化燃烧时,氧气的质量分数急剧下降至接近0%,生成产物CO2和H2O. 固定碳燃烧所需温度更高,在沿炉排前进方向大概1 m处,床层表面温度达到了1 040 K,开始燃烧并生成产物 CO2. 在1.5 m 处,燃烧率达到最高,温度同样也达到最高1 590 K. 到朝炉篦前进方向2.8 m处,基本完成燃烧,床温下降到966 K. 从图5还可也看出在燃烬后段所需氧气很少,因此消耗量少,故其质量分数接近空气中氧气的质量分数.图5显示挥发分组分变化过程中会产生两个峰值: 第一个峰是由于受热后挥发分逸出,逸出的同时会吸收热量,而且热量尚未完全传递到内部,所以逸出量逐渐减少; 此后,由于里层燃料受热温度升高,之前处于内层挥发分随之逸出,因此第二个逸出峰也随之出现. 本文假定挥发分主要组分是CmHn、CO以及H2,由图6知,CmHn、CO 和 H2也出现相同逸出规律. 从图6中可以看出水的质量分数变化出现三个峰: 第一个峰是由于温度升高燃料中的水分蒸发; 第二个峰是由于挥发分中H2析出与O2反应生成了H2O; 而第三个峰的出现也是由于热传导使得内层燃料中挥发分逸出后与氧气发生反应而后形成.

图5 燃料中各组分的变化率

图6 烟气各组分质量分数变化

3.3 运行优化分析

为了避免因配风原因影响燃烧,所以在实际运行中, 应根据各个阶段不同的风量需求供应不同风量. 例如本文中,在朝炉篦前进方向0~0.5 m段水分脱出,此时并没有氧气消耗,理论上无需供风;朝炉篦前进方向0.3~3 m 为挥发分逸出及固定碳燃烧段,需要大量消耗氧气,此段供风量应达到总风量的80%~90%[11];朝炉篦前进方向3~4.3 m,属燃烬段,也只需少量氧气供应,此段风量供应可以调整为总风量的10%~20%,按需配风,并以此提高燃烬率,同时也可避免炉膛因供风量过大导致烟气热损增加. 针对朝炉篦前进方向0.5~2.1 m处高温烟气,可以在锅炉进行结构设计时适当改变前、后拱角以及其覆盖范围,以此增强辐射,加速燃烧. 鉴于本文一次风温度未经过预热,挥发分析出后料层变薄,同时一次风温度较低使得料层温度降低,进而达不到固定碳燃烧温度,使得固定碳燃尽率只有71.65%,所以为了提高燃尽率,可以适当提高一次风温度.

4 结 论

本文对我国东北三省某用于粮食加工的层燃热风炉进行数值计算,并获得朝炉篦前进方向固相、气相温度分布,以及烟气中的主要成分在朝炉篦前进方向的变化规律. 结果显示,层燃过程可分成水分脱出、挥发分及固定碳逸出燃烧3个阶段. 根据不同阶段特点,按需配风,从而减少烟气热损,保障燃烧效率,并建议根据烟气温度设计拱角以及覆盖长度,增强炉拱辐射强度,加速燃烧进程. 针对一次风温低导致固定碳燃烬率低的情况,建议适当提高一次风温来提高固定碳燃尽率. 本文模拟结果对生物质层燃锅炉的运行与结构优化设计提供了参考.

猜你喜欢

风量生物质组分
生物质水暖炉 农村节能减排好帮手
数据中心间接蒸发冷却空调机组二/一次风量比
近红外定标法分析黏/锦/氨三组分纤维含量
某乘用车冷却系统进风量仿真及优化
生物质炭添加对土壤氧化亚氮排放影响的研究进展
组分分发管理系统在天然气计量的应用
生物质发电:秸秆变燃料环保又增收
定风量调节阀在通风空调系统中的应用
超超临界660 MW机组二次风量异常下降分析与研究
神奇的地下世界