内循环流化床冷渣器内颗粒混合特性数值模拟.docx

1、内循环流化床冷渣器内颗粒混合特性数值模拟内循环流化床冷渣器内颗粒混合特性数值模拟摘 要本文分析了目前国内普遍使用的各种循环流化床锅炉冷渣器的特点，并且根据国内冷渣器的使用现状，选用简化的冷渣器模型作为研究对象，采用数值模拟计算方法，利用计算流体动力学(CFD)的FLUENT软件对渣器进行研究。针对冷渣器内的气固两相流体动力特性，文中选用欧拉模型并结合颗粒动力学理论，研究冷渣器内颗粒流动与混合特性以及颗粒在冷渣器内的循环流动规律等;最后将计算结果与实验结果进行分析、比较和修正。研究结果表明:此模型能够较好地反映冷渣器内部的气固两相流体动力特性，该研究结果为冷渣器结构的优化设计提供了可靠的依据，所

2、得的数据和运行规律为冷渣器的可靠运行提供参考依据和操作指导。关键词：循环流化床 冷渣器 数值模拟 气固流动Circulating Fluidized Bed Ash Cooler Numerical Simulation of Particle MixingABSTRACTIt analyzes the features of many kinds of slag cooler of commonly used Circulating Fluidized Bed (CFB) boilers, and according to the domestic use of slag cooler st

3、atus, uses a simplified model of slag cooler as the research object. The method of Numerical simulation and FLUENT software of CFD was used to research slag cooler. For the fluid dynamic characteristics of the gas-solid in slag cooler, it takes Eulerian model combined with the kinetic theory of gran

4、ular flow (KTGF) to research the flow and mixing characteristics and circulation flow pattern of the particle in the slag cooler. Finally, the results show that this model can reflect gas-solid flow characteristics well in the slag cooler through comparing calculated results with experimental result

5、s. The results would provide a reliable basis for the optimization of the structure of slag cooler, the dates and operation discipline would provide reference parameters and operational guidance for credible working of slag cooler.Key Words: CFB boiler; slag cooler; numerical simulation; gas-solid f

6、low第一章 引言1.1 选题背景和意义能源工业是支撑国民经济可持续发展的基础性行业，关系到社会生产、人民生活和国家安全。进入二十一世纪以来, 我国的工业化与城市化进程加快, 能源需求正以前所未有的速度增长。“ 七五” 、“ 八五” 、“ 九五” 三个计划期, 我国能源消费量年均增长速度分别是6.5%,6.0%，-0.07% , 能源消费弹性系数始终小于“1 ” 。进入“ 十五” 以后, 能源消费速度开始加快, 并且逐步超过同期国内生产总值的增长速度, 平均能源消费增长速度和平均能源消费弹性分别达到10%和1.05。2003年、2004年能源消费弹性系数高达1.53和1.59。2005年能源消

7、费量为22.33亿吨标煤,是2000年的6倍, 能源消费绝对量增长了8.5亿吨标煤。其中, 煤炭消费增长了64%, 石油消费增长了45%, 电力消费增长了85%1。作为燃煤排放主要污染物的CO2导致全球的气候变暖，造成温室效应，SO2 带来全球大面积的酸雨，导致植被破坏等严重危害，这些均已引起全人类的关注2-5。我国每年排入大气的87%的SO2和67%的NOX均来自煤的直接燃烧6，由此可见，煤的高效率、低污染燃烧技术对社会的可持续发展具有重要的意义。循环流化床燃烧技术是近40年发展起来的一种高效、低成本的清洁燃烧技术，也是目前商业化程度最好的洁净煤燃烧技术，发展迅速。它具有以下几个优点：（1）

8、燃烧效率高，接近或达到同容量煤粉锅炉的效率水平；（2）燃料适应性强，不仅可以燃用烟煤等优质煤，而且可燃用各种劣质燃料，如褐煤、贫煤、洗中煤、泥煤、石焦油甚至工业废弃物和城市垃圾等；（3）符合调节比宽，在25%的额定负荷下仍能稳定燃烧；（4）负荷调节方便快捷，负荷连续变化速率可达（7%-12%）/min；（5）污染排放物少，低温燃烧及分级送风使NOX生成量少，可用石灰石做脱硫添加剂，低成本实现炉内脱硫；（6）灰渣便于综合利用6。特别在燃用高灰份、低挥发份或高硫份等其它燃烧设备难以适应的劣质燃料以及低负荷要求较高的调峰电厂和负荷波动较大的自备电站中，循环流化床锅炉是最佳选择，是当前公认的燃煤技术的

9、重大创新。我国是世界上少数几个以煤为主要能源的国家之一，发展高效、清洁的循环流化床锅炉，具有长远的战略意义。综上所述，循环流化床比较适合我国国情，因此发展循环流化床技术有利于我国能源工业以至社会的可持续发展。随着循环流化床燃烧技术的迅猛发展，对其高温炉渣的冷却热回收技术的研究日益得到国内外研究者的重视，冷渣器作为保证流化床锅炉安全、高效运行的重要部件，它的不正常工作是导致CFB机组停炉和减负荷运行的主要原因之一7。冷渣器的开发应用;对降低高温炉渣的热污染;减轻排渣的劳动强度;回收高温渣的热量以及提高锅炉热效率都起到了重要作用。尤其为燃用低热值、高灰份劣质煤的流化床锅炉配用合适的冷渣器就显的更加

10、重要。当循环流化床燃用石煤、煤研石等低热值、高灰分劣质燃料后，所排出的热渣量大，不仅造成了热损失同时也污染了环境。它不仅恶化了运行现场的工作条件，灰渣中残留的硫、氮仍可以在炉外以SO2和NOx形式排放出来。同时，高温热渣给运输和处理造成了困难;使得机械化操作难以进行。一般循环流化床锅炉的底渣温度在800-950左右8，这不仅造成了大量的热损失，降低了锅炉效率，而且对环境也会造成一定的影响。目前，冷渣器运行过程中普遍存在许多问题，如冷渣器进口堵塞，冷渣器内出现结焦，排渣温度高造成排渣口结焦或堵塞等。这些故障与冷渣器内部的结构和流体动力特性有着密切的关系。因此，对冷渣器内流体动力特性的试验研究以及

11、模拟冷渣器内颗粒流动特性对于改善冷渣器性能，提高循环流化床锅炉热效率等有着重要意义。1.2 国内外研究现状1.2.1 冷渣器研究现状与国外CFB锅炉运行情况相比，国内CFB锅炉底渣排放特点有:底渣排放量大，粒度及排放量变化范围大9。目前，国内使用的冷渣器设备10主要有：风水联合冷渣器、滚筒冷渣器、钢带风冷式冷渣器、水冷螺旋式冷渣器等。在国外，不少循环流化床电厂选用螺旋式冷渣器即水冷绞龙冷渣11，此类冷渣器材料选用较好但造价很高。美国Foster Wheeler公司开发出流化床分选冷渣器，冷却效果较好，但存在埋管磨损问题；B&W公司正在研制开发的冷渣器与浙江大学首先提出的流化移动叠置式灰渣冷却技

12、术的概念较为相似12,13。下面，简单的介绍一下各类冷渣器的优缺点。（1）风水联合式冷渣器布置于炉膛两侧，整体采用绝热钢板结构，外部采用保温材料进行绝热处理，内部铺设耐磨耐火浇注料以防止灰渣磨损钢板。风水联合式冷渣器优点主要包括：煤质适应性强，冷渣器的设计留有较大裕量；低速流化，使冷渣器始终保持在鼓泡床状态运行，可在保证换热前提下大大减弱材料的磨损，延长了冷渣器使用寿命；独立的流化风机和钟罩式风帽，冷渣器配有专用的罗茨风机，可保证冷渣器流化风具有稳定的压力和流量；倾斜的床面及多点排渣，冷渣器各风室底部耐磨耐火材料敷设成倾斜状，分别向排渣口倾斜3-4度，可使冷渣器底部的粗渣自动流向粗渣排渣口14

13、,15。但从长期的运行中发现，该冷渣器还存在许多问题：灰渣复燃结焦；对进渣粒度要求较严格(一般小于20mm 时排渣顺畅)，处理大渣的能力不够，有时会出现堵渣；因细灰未能有效的被分离而导致热风管道堵塞；床内埋管磨损；送风系统设计不足，造成调节困难。（2）滚筒式冷渣器是目前国内使用较多的一类冷渣器。这类冷渣器的冷却主体是由内外两个滚筒组成的双层密封套筒，中间的空隙通冷却水。有的还在内筒的内壁上焊接固定的螺旋叶片，供物料导向前进。冷渣器由电动机驱动做旋转运动，热渣进入内筒由螺旋叶片推动前进，同时与外筒逆流的冷却水接触降温。冷却风是由出料口进入由入料口排出，逆向与热渣直接接触。在风水联合的滚筒冷渣器中

14、，冷却风一般是辅助性的16。滚筒冷渣器也是目前比较适合我国国情的，用于CFB锅炉较多的冷渣设备，但主要用于中、小型CFB锅炉，出渣量小，出渣温度较高。其特点包括：安全可靠，设有低水压和超水温报警装置；冷却能力强，灰渣的热回收率高达90%以上：冷却水温升大，进出水温差达60；调整负荷简单方便；设备密封性好，极大的改善了工作环境；低磨损，使用寿命长。（3）钢带冷渣器主要部件是超强钢带，它是由耐高温钢制成，并与其它部件密闭在壳体中。冷却介质是冷却风，钢带兼具冷渣及输送功能。整套系统的启动与停运由PLC控制，设置超温保护、零转速保护及电气设备保护。钢带冷渣器的工作过程为：超强钢带由头部主动滚筒驱动，将

15、热渣从入渣口输送到出渣口，与此同时与逆向的冷却风充分接触来降低灰渣的温度。该类冷渣器一般分两级，每级可以由一条或多条钢带组成；每条钢带的出口处都设有紧急喷淋装置，防止底渣中可燃物含量超出设计值14。钢带冷渣器主要靠国外引进，价格较高，目前国内很少使用这种冷渣器，据调查使用效果较好，确保了机组的稳定运行。但风冷换热效果差，降低了锅炉热效率。（4）螺旋式冷渣器即水冷绞龙结构主要有单螺旋轴或多螺旋轴，冷却方式为外壳、叶片、和轴水冷，能够在输送过程中实现高温灰渣的间接冷却。采用水冷螺旋式冷渣器除渣系统的主要优点是系统结构及控制简单，通过控制调速电机的转速来调节锅炉的排渣量；缺点主要是冷却能力小，对煤种

16、的适应性差，设备抗热冲击能力差及磨损严重。有的电厂将水冷螺旋式冷渣器和纯风冷冷渣器串联使用，取得了较好的效果，既解决了设备的热冲击问题，又提高了系统的除渣能力；但其系统结构及控制都比较复杂17,18,19。此外，国外曾在80年代初提出并试验了许多越过高压力位差直接输送粒状固体的方法，但至今仍未成功。压力下的排渣系统主要是应用在PFBC高温炉渣的排放方面。英国Leatherhead PFBC的排渣冷渣系统，高温炉渣连续的由布风板底部的排渣管引出，依靠一气控排渣机构来控制排渣速率；热渣用气力输送方式，经过带水夹套的冷却管冷却后送入床料收集斗，再排入下面的变压锁斗卸压；上述方案的显著优点在于排渣锁斗

17、在空间布置上拥有了很大的自由20。德国Babcock PFBC的排渣系统，高温炉渣从流化床床层底部向下移动时，得到一小部分流化风和水冷壁的有效冷却，到出渣口处，渣温已冷至300左右；压力壳外有一个在加压环境下工作的叶轮加料器，用以控制排渣速率；冷渣排入锁斗系统，最后气力输送至灰仓；这种排渣冷渣方式突破了以往的思路，使冷渣过程在炉内完成，从而提高了空间利用率，并省去了昂贵的高温隔离阀21。我国“八五”攻关的PFBC-CC中试装置22，其高温炉渣连续排放，冷却系统采用非机械方式控制排渣速率；高温炉渣在压力下以浅床流化方式进行冷却；冷渣热风回送至炉内沸腾段作二次风；而经过冷却后的炉渣则由变压罐卸压后

18、排入常压环境；东南大学热能工程研究所建立了高温炉渣连续排放及冷却系统的冷态试验台，并进行了冷态实验23,24,25，灰渣冷却器的设计采用了内循环浅床的形式26。 1.2.2 试验研究概况试验研究流化床冷渣器内颗粒运动主要采用两种方法：追踪冷渣器内单个颗粒的运动轨迹与测定冷渣器内示踪粒子的混合特性。前者通常采用放射性的示踪粒子或比较容易观察的异于物料的粒子27，后者因需要测量示踪粒子的浓度大小，一般采用具有相同物理性质的两种粒子进行混合，然后通过某种不同性质来区分它们，如（1）化学性质区分法；（2）放射性粒子法；（3）磁性粒子法；（4）升华法；（5）热粒子法。方法（1）-（4）都有一个共同的缺点

19、：从物料中分离示踪粒子相当困难，而且难以达到快速取样，更为重要的一点是冷渣器内气泡运动过程具有一定的随机性，两个相同工况的运行过程是不完全相同的。采用加热的物料颗粒作为示踪粒子的方法具有以下优点：（1）示踪性质与物料相同；（2）流化过程中可在线测量示踪粒子的浓度；（3）同时测量冷渣器内不同位置的示踪粒子的浓度。热粒子与物料在粒径、形状、密度方面完全相同运行后不必分离，几分钟后便与物料充分混合，炉渣的性质与实验前后完全相同，在此短暂的时间内可进行多组测量。另外还满足了由于气泡和颗粒运动具有一定的随机性，同一混合实验需要重复若干次后，才能获得准确且有代表性的实验数据。Valenzuela认为28，

20、29由于固体物料与流化空气的热容量之比相当高，在较短的时间里空气带走的热量与示踪颗粒带入的热量相比可忽略不计。空气的对流换热系数较高，空气的热容量又较低，可以近似认为床内任一点空气与颗粒具有相同的温度，微元内温度的升高正比于示踪粒子的浓度。沈来宏27采用此法在床体尺寸为130092060mm3的流化床上，研究床内颗粒混合并建立流化床内颗粒混合的对流-扩散模型，实验数据与模型数值结果吻合良好。任有中和Gibbs30采用此法研究床内颗粒横向扩散系数和热扩散系数，热颗粒与床料温差30，在较短时间内（小于60s），忽略流化空气穿过床层所带走的热量，理论解与实测到的床层温度能吻合一致，但随着时间的延长，

21、两者的相差会越来越大。杨海瑞31，32通过分析微元内的质量能量守恒，从理论上证明了在一定的条件下，微元内的温升与示踪剂的浓度近似成正比。1.2.3 数值模拟概况合理的物理和数学模型是模拟的关键。流态化气固流动系统有以下两个特点：(1)颗粒相的体积份额较大，它对流体相的运动规律有较大的影响，之间存在强烈耦合；(2)颗粒之间的相互作用大大加强，并对整个系统的运律有较大的影响。常用的颗粒轨道、无滑移(或单流体模型)和小滑移等稀相模型对流化床系统不再适用。近年来，对气固两相流的模拟主要有两种基于欧拉方法的颗粒相拟流体模型、基于拉格朗日方法的颗粒轨道模型。1.2.3.1 颗粒相拟流体模型拟流体模型属于流

22、体力学方法中的连续介质模型。应用于两相流动时，也称为双流体模型(Two Fluid Model,TFM)，最早由是Jackson33，Murray34和Anderson和Jackson35等人提出的，主要思想是求解局部平均的Navier-Stocks方程。在颗粒相拟流体模型中，在把流体相视为连续介质的同时，把固体相也视为一种类似流体的连续介质(拟流体)，并根据守恒原理分别建立气固相的瞬时、局部动量控制方程组。在这类模型中，空间的每一点都是各相所共有，即两相相互渗透且各相之间存在着相互作用。近二、三十年来基于颗粒相拟流体模型的流态化过程数值模拟取得了很大的进展，其中最为著名的是Gidaspow等

23、人的研究小组。在文献36,37中，他们对鼓泡流化床内气泡的形成、长大及破裂过程进行了成功的预测(见图1)，通过计算也得到了流化床内流体及固体颗粒相的速度与床层中空隙率的分布。在文献38,39中，他们对PYROFLOW型循环流化床锅炉的简化模型进行了数值模拟，通过模拟得到了整个循环流化床回路内气固两相的速度及固体颗粒浓度的瞬时分布。Kuipers等人则采用拟流体模型研究了三维流化床上升段内的流体动力学特性40。另外，Gidaspow等人也将拟流体模型从单纯计算流化床的流体动力学特性推广到计算流化床内壁面与床层之间的传热特性41。Enwald等人在文献42中曾对拟流体模型在流态化领域的应用做了较为

24、全面的总结。然而，自从该模型提出时开始就存在着一些难以解决的问题，主要包括：(1)连续介质假设的争议传统流体动力学应用了两类变量来描述流动现象。一类用以表征流动现象，即速度、密度、温度和压力等物理量；另一类用以表征力和前一类变量之间的关系，如粘度、比热、导热系数等。根据连续介质假设，这些参数在流体所占据的空间形成连续场，且变化得足够缓慢，服从局域平衡假设，可用流体力学微分方程组描述力和流之间的关系。然而，由于存在非线性、非平衡性，实际系统未必满足这一假设。为将颗粒相视为连续介质，双流体模型中的微元体必须远大于颗粒间的平均距离及颗粒的平均自由程。但在这样的尺度上，结构和秩序，或者说非平衡现象可能

25、已经相对明显了。由于拟流体模型是从连续的观点来处理问题，忽略了颗粒的个体离散行为。因此所得计算结果也只能反应多个离散颗粒在局部区域内的平均特性而无法在颗粒的层次上获得其运动、燃烧等动态信息，这也是双流体模型的一大局限。(2)颗粒相本构方程的建立将实际的离散颗粒相处理成连续介质时会遇到很多困难，这些困难主要集中在如何确定颗粒相拟流体的粘性应力和静压力。目前常用的有两种方法，一是采用“固体弹性模量”(Solid elasticity modulus)来描写颗粒相应力。它是通过变换颗粒应力张量而导出的称之为“颗粒压力”的一个经验表达式43。另一种方法是采用“颗粒动理论”(Kinetic theory

26、 granular media)来获得44,45。在该理论中，通过与气体分子运动论相比拟，引入了一个类似气体温度的物理量，即所谓的“颗粒温度”。固体粘性应力和静压力是“颗粒温度”的函数。借助于这一理论，实际固体颗粒间的法向相互作用力被表达为固体颗粒相的静压力和体积粘性，而切向作用力则归结为剪切粘性。但是该理论对颗粒间保持双体碰撞的情况适用，对于实际流化床中常存在的多颗粒碰撞和颗粒间长时间相互接触的情况则不太适用。(3)模型控制方程组仍存在争论Lyczkowski等人发现基于Anderson-Jackson方程的微分方程组(他们称为模型A)在计算中对初值问题适定性不好，进而他们提出了称为模型B的

27、改进方程组46。Boemer等人则根据颗粒动力学理论发现，虽然模型B的结果较接近实验值，但模型A在理论上要比模型B更准确些47。另外，对颗粒相拟流体方程的边界条件的确定也需要一定的经验和技巧48。迄今为止，采用颗粒相拟流体模型时，如何获得收敛的、正确的结果仍然是一个非常有挑战性的问题。1.2.3.2 颗粒轨道模型与颗粒相拟流体模型不同，在颗粒轨道模型中，虽然仍将流体相看作连续介质，固体颗粒相则按离散体系处理。颗粒轨道模型在气固两相流研究中有着悠久的历史，早期的颗粒轨道模型也称为单颗粒动力学模型，它不考虑固体颗粒之间及固体颗粒对流体相的作用，只计算流体对离散固体颗粒的作用。由于上述假设，单颗粒动

28、力学模型只是在稀相流动，即固体颗粒浓度较小的气固两相流，模拟中近似适用，对流态化过程这样的稠密气固两相流则无能为力。后来的研究者对单颗粒动力学模型进行了改进使之能够充分考虑颗粒与流体之间的相互作用，并且也计及了与颗粒扩散无关的、两相间的大速度滑移和温度滑移。目前，颗粒轨道模型是稀相气固两相流模拟中最为广泛应用的模型。但该模型对于固体颗粒之间的相互作用仍未考虑（对确定性轨道模型）或考虑不充分（对随机轨道模型）。另外，对于颗粒浓度较高的情况，计算量较大且不容易获得收敛的结果，因而该种模型也不适于流态化过程的模拟。近年来，属于颗粒轨道模型范畴的离散单元法模拟(Discrete Element Met

29、hod，DEM)在流态化模拟中取得了很大的成功。DEM方法最早在1979年由Cundall等人提出并用于模拟地质学上的岩石运动规律49。Tsuji等人在1992年首次将该方法应用于浓相水平管道气力输送研究中50，1993年又将该模型与计算流体力学模型(Computational Fluid Dynamics，CFD)相结合，成功用于流化床内气固两相流动过程的模拟51。目前，DEM已发展成为粉体系统和流化颗粒系统模拟的一种非常重要的方法，它处理离散颗粒相的思路是：在拉格朗日坐标系中对系统的每一个颗粒的受力和轨迹进行跟踪研究，颗粒之间的相互作用力根据有关模型直接确定。这样就避免了颗粒相拟流体模型中

30、建立颗粒相输运系数本构方程及固体颗粒静压力表达式的困难。对于颗粒与流体之间的相互作用，则与拟流体及一般颗粒轨道模型类似，采用有关的经验关系式确定。由于DEM同时考虑了颗粒之间及颗粒与流体间相互作用，因此该模型非常适合于流化床这样的浓相气固两相流动。DEM方法与颗粒相拟流体模型相比，其不足之处是由于要跟踪系统中的每一个粒子的受力和运动情况、并且时间步长又受到限制，因此其计算量较大。但由于DEM模拟能够在颗粒水平上给出流态化系统动态变化过程的非常丰富的受力和运动信息，在流态化系统微观机理研究领域显示出巨大的潜能，在国际学术界备受关注。目前大多数的研究人员都是通过实验的方法来研究冷渣器内流体动力特性

31、，如排渣速率、流化风速、配风比和床层压降等。这种主要依靠实验结果再加上相似性原理整理数据而得到的经验或半经验公式方法，也有不少缺陷，例如它经验性强、范围窄、结果多依赖于实验条件以及实验条件不能适应广泛变化等，测量也会有误差，且在大多数情况下，这种全比例的实验是极其昂贵的，往往也是不可能的，而缩小比例的实验结果外推到全比例的设备上时又多会失真。由于实验条件的限制，单纯通过实验来研究冷渣器的性能不仅周期长而且费用高。另一方面随着现代电子计算机、数值分析学、计算流体力学、计算传热学及计算燃烧学等学科的产生和发展，自上世纪七十年代以来，逐步形成了两相流动和燃烧数值模拟7。综上所述，建立合理的物理及数学

32、模型，利用欧拉方法、拉格朗日方法、颗粒动力学理论等，并借助于计算机模拟技术来对冷渣器的流动细节进行研究，对内流体流动规律进行分析，进而优化冷渣器的结构，可以大大缩短研发周期，具有重要的工程应用价值。1.4 研究思路和技术方法1.4.1 研究目标、研究内容高压高温灰渣的连续排放和冷却装置（冷渣器）是内循环流化床的关键技术之一。能否实现良好的气固混合和物料平衡是实现灰渣的稳定连续排放和冷却技术的重要环节，采用内循环流化床冷渣器系统可较好的解决这一问题。到目前为止，内循环流化床冷渣器内的颗粒流动和混合规律的掌握仍然较少。本课题利用数值模拟的途径研究冷渣器内颗粒混合特性。通过研究使之得出冷渣器进渣量与进风量等的最佳操作工况。1.4.2 拟解决的关键问题1）深入理解冷渣器内气体和固体的流动规律；2）建立合理的物理模型及气固湍流数学模型；3）重点在于获得可靠的结果和优化参数。1.4.3 课题的研究方法和技术路线分析本文力求利用成熟的计算流体力学商业软件Fluent，并选用合适的计算模型对冷渣中的气固两相流动进行数值模拟，为进一步改进现有