精华公布栏

发信人: realplayer (演绎真实人生), 信区: Energy
标  题: 燃用宽筛分煤循环流化床锅炉燃烧模拟计算
发信站: 哈工大紫丁香 (2002年03月22日21:41:46 星期五), 站内信件

原文参考：http://www.china5e.com/dissertation/cfb/0058.htm

燃用宽筛分煤循环流化床锅炉燃烧模拟计算

刘文铁　李炳熙　赵广播　谭　袖　陆慧林

-----------------------------------------------------------------------
---------

［摘要］　建立了循环流化床锅炉炉膛颗粒燃烧和脱硫反应模型。该模型考虑了炉
膛下部为高颗粒浓度的密相区和上部为低颗粒浓度的稀相区的特征。模拟计算给出
了烟气温度、热流密度和各气体成分(O2、C2O、CO、H2O和S2O)的轴向分布。模拟
计算结果的趋势是合理的。
关键词：循环流化床锅炉；数值模拟计算；燃烧

1　前言
　　循环流化床锅炉具有燃用煤种广，燃烧效率高且又能满足环境排放标准，显示
出明显的经济效益、环境效益和社会效益。在过去的十多年里，该技术取得了长足
发展，但由于该过程涉及到流动、传热和燃烧等多因素非线性耦合的相互影响，依
然存在着如何设计和放大循环流化床锅炉的挑战。数学模拟方法能够较全面地反映
多种因素的相互影响，提供一个综合性的定量评估，甚至揭示出难以自觉而有价值
的信息。这对提高设计和放大的可靠性有重大意义。因而循环流化床锅炉的燃烧数
学模型被逐渐建立和发展。Weib et al.(1987)和Maggio et al.(1995)采用单元法
建立了包括炉膛、分离器、回料系统和外置式热交换器等的数学模型。Hyppanen
et al.(1990)建立的三维循环流化床锅炉炉膛模型能够预测烟气温度、气体成分分
布等参数。然而几乎所有的上述数学模型都末能很好地考虑炉膛下部高颗粒浓度区
的颗粒流动和燃烧。Johnsson & Leckner(1995)研究结果表明在炉膛下部的平均颗
粒密度大约为800～1200　kg/m3，他们认为该区域具有鼓泡流化床的特征，而被称
为密相区。Montat et al.(1995)在一台电功率250MW循环流化床锅炉试验表明密相
区的颗粒浓度约700～1000 kg/m3，说明炉膛下部颗粒流动与湍动鼓泡流化床相似
。密相区高度约1.0米。实际运行和设计经验也表明燃用宽筛分煤循环流化床锅炉
炉膛下部密相区必须运行在鼓泡流化床操作范围内。从而沿炉膛高度形成不同的流
动、传热和燃烧过程，建立的数学模型应该考虑该特征对其过程的影响。

2　数学模型
　　假设循环流化床锅炉炉膛内的颗粒流动、传热和燃烧过程为稳态过程，沿炉膛
高度分为下部密相区和上部稀相区。
2.1　密相区
　　假设：(1)密相区由气泡相和乳化相组成，乳化相空隙率取为临界流化空隙率
，超过临界流化速度的气体通过气泡相。(2)沿径向各相参数均匀分布，沿轴向为
逆流返混模型。(3)煤的燃烧分为三个过程：加热蒸发、挥发份析出和焦碳燃烧。
当颗粒温度达到100℃时蒸发过程终止。挥发份析出过程同时受化学反应速率和质
扩散速率控制，当颗粒温度达到密相区温度时该过程结束。碳粒燃烧服从缩核模型
。(4)颗粒为宽筛分，粒子为球形。
　　煤燃烧主要生成C2O，CO和H2O等气体组分。沿轴向微元高度气泡相和乳化相内
各气体组分的质量守恒方程如下。Lu等(1998)给出了上述方程的详细推导过程和各
项的计算方法。

（1a）

　

（1b）

（2a）

　　

（2b）

　

（3a）

　

（3b）

　　对单一颗粒，能量守恒方程是：

　　

（4）

其中：

（4a）

（4b）

　　对惰性颗粒(石灰石颗粒和循环颗粒)，能量守恒方程简化为：

　　

（5）

　　颗粒粒径d(i)的质量守恒方程可表示如下：
　　［来自稀相区的颗粒流量］－［流出密相区的颗粒流量］＋［来自颗粒d（i+
1）燃烧的流量］－［来自颗粒d(i－1)燃烧的流量］＋［颗粒d(i)燃烧减少的颗粒
流量］＝0，或者表示为

　　Wfall,i－Ei＋Wi+1－Wi－Wburn,i=0

(6)

2.2　稀相区
　　假设：(1)忽略径向的不均匀性，燃烧反应按一维轴向考虑。(2)不考虑煤粒的
干燥加热过程。(3)忽略气相和颗粒的径向交换。
　　各气体组分的质量守恒方程如下：

　　

（7a）

　

（7b）

　　

（7c）

　　能量守恒方程可表示如下：

（8）

　　轴向颗粒浓度分布按Kunii & Levenspiel(1991)关联式计算：

　　　　

（9）

　　其中系数a取为3.0/u，εsd为可调节参数。
2.3　石灰石脱硫反应
　　石灰石脱硫反应过程为：

　　

　　反应速率按下式计算(Saraiva et al., 1993):

（10）

　　

（11）

　　SR＝5.6×104－38.4Tp　　1253K

（12a）

　　SR＝－36.7×104＋35.9Tp　　1253K

（12b）

　　其中λ是石灰石活性系数。
2.4　边界条件
　　在密相区入口，气体组分满足：

　　CO2，b＝CO2，e＝CO2，inlet;　CH2O＝CCO2＝CCO＝0

（13）

　　在密相区出口，各气体组分的平均值按下式计算：

　

（14）

3　数值模拟计算和分析
　　上述描述密相区和稀相区颗粒流动、传热和燃烧微分方程可采用龙格—库塔法
求解，可确定在不同操作条件下的密相区和稀相区各气体组分、温度、热流密度、
颗粒温度、颗粒粒度分布和颗粒含碳量等参数。
　　模拟计算是某台58　MW高温热水循环流化床祸炉。为简化计算，取平均炉膛横
截面为2.5m×3.1m,高为8.4m。分离器分离效率为88.5%。表1给出燃用燃料特性。
图1表示所用燃料颗粒的粒度分布。平均煤粒粒径为2.38mm。图中也给出密相区颗
粒的颗粒粒度分布，平均颗粒粒径为1.86mm。计算结果表明随着颗粒直径增大，颗
粒达到操作温度所需要的时间增大。图2表示燃料、密相区排出的溢流颗粒和炉膛
出口扬析颗粒碳流量分布。密相区颗粒平均含碳量为2.45%。图3表示烟气温度沿炉
膛高度分布。密相区平均烟气温度为946℃。计算结果表明沿炉膛高度烟气温度逐
渐下降。图4表示热流密度沿炉膛高度的变化。由图可见沿炉膛高度热流密度逐渐
下降，入口空气流速对热流密度影响很小，其变化趋势与烟气温度变化是相似的。

表1　燃用燃料特性

C(%dry)
70.2
N(&dry)
1.4
A(%)
16.5

H(%dry)
4.9
S(%dry)
0.71
C(%)
31.9

O(%dry)
6.1
W(%)
3.9
Qdw(kJ/kg)
26.83

图1　燃料和密相区颗粒粒度分布

图2　燃料、扬析和溢流颗粒碳流量分布

图3　烟气温度沿炉膛高度分布

图4　热流密度沿炉膛高度的分布

图5　各气体组份沿炉膛高度的变化

　　图5表示各气体组分O2、C2O、CO和H2O沿炉膛高度的变化。计算中一次风与二
次风比例为75/25(%)。二次风位置在1.2 m处。由图可见各气体组分沿炉膛高度下
降，较高的O2浓度和较低的CO2浓度表明炉膛稀相区燃烧强度低于密相区。图6表示
Ca/S模尔比对S2O浓度分布的影响。计算中石灰石平均粒径为150微米。正如大量实
验和模拟计算所表示的S2O排放量随着Ca/S模尔比增大而较低。模拟预测的趋势是
合理的。

图6　Ca/S模尔比对SO2分布的影响

　　主要符号　　　　　　
　　A＝面积；p＝压力；
　　C＝浓度，比热；R＝总反应率；
　　d＝颗粒粒径；T＝温度；
　　H＝发热值；U＝速度；
　　h＝传热系数；V＝体积；
　　m＝颗粒质量；Z＝高度
　　希腊　　　　　
　　ε＝空隙率；
　　ρ＝密度；
　　下脚标
    f＝煤颗粒；
    g＝气体；
    mf＝最小流化；
　　α＝放热系数
　　s＝惰性颗粒；
　　w＝壁面　　

致谢：本项目研究得到清华大学煤清洁燃烧技术国家重点实验室基金的资助

本文联系人：陆慧林(1943～．)　男　教授　150001　哈尔滨工业大学热能工程教
研室

作者单位：(哈尔滨工业大学)　

参考文献

［1］Hyppanen T, lee Y Y & Rainio A. A three dimensional model for circulating
fluidized bed combustion. Circulating Fluidized Bed Technology Ⅲ， 1990
:563～568.
［2］Johnsson F & Leckner B. Vertical distribution of soilds in a CFB furnace
, 13th Int. Conf. On FBC, Orlando, USA∶1993.
［3］Kunni D & Levenspiel O. Entrainment of solids from fluidized beds.
Powder Technol, 1990,61:193～206.
［4］Leckner B, Golriz M R, Zhang W, Andersson B A & Johnsson F. Boundary
layers-first measurements in the 12 MW CFB research plant at Chalmers
University. Proc. 11th　Int. Conf. on Fluidized Bed Combustion, ASME, 1991
:771～776.
［5］Lu Huilin, Zhao Guobuo, Bie Ruoshan, Yang Lidan. Numerical computation
of a circulating fluidized bed combustor, Intern.J. of Energy Research
, accepted, 1998.
［6］Montat D & Dimaggio T. 1D two phase description of the thermal hydraulic
behaviour of the furnace of E.Huchet 125 MW CFB bolier, 5th Int. Conf.
CFB, Bejing, China∶1995.
［7］Saraiva P C, Azevedo J　L　T　 & Garvalho M G. Mathematical simulation
of a circulating fluidized bed combustor.Combst. Sci. and Tech, 1993,93
:223～243.
［8］Weib V, Fett F N, Helmrich H & Janssen K. Mathematical modeling of
circulating fluidized bed reactors by reference to a solids decomposition
reaction and coal combustion.　Chem. Eng. Process, 1987,22:79～90.
［9］Werther J & Wein J. Expansion behavior of gas fluidized beds in the
turbulent regime. AIChE Symp., 1994,90:31～49.

　

--
生活需要realplayer!

※ 来源:·哈工大紫丁香 bbs.hit.edu.cn·[FROM: 202.118.243.253]
※ 修改:·realplayer 於 03月22日21:42:15 修改本文·[FROM: 202.118.243.253]