煤粉炉水冷刮板式干排渣技术研究.pdf

Advances in Energy and Power Engineering 电力与能源进展电力与能源进展,2023,11(3),101-111 Published Online June 2023 in Hans.https:/www.hanspub.org/journal/aepe https:/doi.org/10.12677/aepe.2023.113013 文章引用文章引用:李永华,付志.煤粉炉水冷刮板式干排渣技术研究J.电力与能源进展,2023,11(3):101-111.DOI:10.12677/aepe.2023.113013 煤粉炉水冷刮板式干排渣技术研究煤粉炉水冷刮板式干排渣技术研究 李永华，付李永华，付 志志 华北电力大学能源动力与机械工程学院，河北 保定 收稿日期：2023年5月8日；录用日期：2023年5月18日；发布日期：2023年6月30日 摘摘 要要 由于我国水资源短缺国情，干排渣技术得到提倡，但风冷式干排渣导致的炉底漏风问题无法避免。为避由于我国水资源短缺国情，干排渣技术得到提倡，但风冷式干排渣导致的炉底漏风问题无法避免。为避免漏风问题，提出一种新型间接水冷干排渣方式，采用数值模拟的方法研究了新排渣方式下，底渣的运免漏风问题，提出一种新型间接水冷干排渣方式，采用数值模拟的方法研究了新排渣方式下，底渣的运动状态和冷却效果，对比了新排渣方式和传统排渣方式的优缺点，结果表明：新排渣方式在运输灰渣的动状态和冷却效果，对比了新排渣方式和传统排渣方式的优缺点，结果表明：新排渣方式在运输灰渣的同时加强了灰渣颗粒之间的混合，排渣温度更均匀，由于换热面积的降低导致降温效果略微低于钢带式同时加强了灰渣颗粒之间的混合，排渣温度更均匀，由于换热面积的降低导致降温效果略微低于钢带式运输方式，但仍满足运输阶段的排渣冷却，灰渣温度从运输方式，但仍满足运输阶段的排渣冷却，灰渣温度从523.15 K降低到降低到449.9 K。新排渣结构下，冷却。新排渣结构下，冷却水流速高于水流速高于0.5 m/s后排渣温度不再显著降低。后排渣温度不再显著降低。关键词关键词 锅炉排渣锅炉排渣，结构设计结构设计，离散元仿真离散元仿真，计算流体力学计算流体力学 Research on Dry Slagging Technique of Water-Cooled Scraper in Pulverized Coal Furnace Yonghua Li,Zhi Fu College of Energy and Power Machinery Engineering,North China University of Electric Power,Baoding Hebei Received:May 8th,2023;accepted:May 18th,2023;published:Jun.30th,2023 Abstract Due to the shortage of water resources in our country,dry slag technology is advocated,but air leakage caused by air-cooled dry slagging system is unavoidable.In order to avoid air leakage,a new indirect water-cooled dry slagging method is proposed.Numerical simulation is used to study the motion state and cooling effect of bottom slagging under the new slagging method,and the ad-vantages and disadvantages of the new slagging method and the traditional slagging method are 李永华，付志 DOI:10.12677/aepe.2023.113013 102 电力与能源进展 compared.The results show that the new slagging method strengthens the mixing between ash particles while transporting ash,and the slagging temperature is more uniform.The cooling effect is slightly lower than that of the steel belt transportation method due to the reduction of heat ex-change area,but it still satisfies the slagging cooling in the transportation stage,and the ash tem-perature decreases from 523.15 K to 449.9 K.Under the new slag discharge structure,the slag temperature no longer decreases significantly after the cooling water flow velocity is higher than 0.5 m/s.Keywords Boiler Slag Discharge,Structural Design,Discrete Element Simulation,Computational Fluid Dynamics Copyright 2023 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 目前煤粉锅炉常用的排渣系统有湿排渣和干排渣两种。由于冷却介质不同，再由我国水资源短缺的国情出发，干排渣技术得到提倡。但风冷式干排渣中冷渣空气被认为是炉底漏风，属于无序漏风。在控制炉膛出口过剩空气系数不变的情况下，冷却风量的增加导致通过空气预热的一次风量减少，因此这种泄漏影响燃烧条件和排烟温度1。刘铭媛2指出在 6%漏风量情况下，炉底漏风全部由燃尽风弥补时，漏风对燃烧过程影响很小，但由二次风弥补时，会严重影响燃烧过程，火焰中心上升，煤粉燃尽推迟，减温水量增加。许多学者的研究3 4表明，将冷却风量控制在锅炉所需风量的 1%以下可以忽略其对炉膛燃烧的影响。对于锅炉效率而言，何彦鹏5等研究发现：当冷却风量小于 0.45%才能够保证 MAC 系统对锅炉影响较小，当冷却风量比例增加时，锅炉效率也会降低。风冷式排渣系统在大排渣量下灰渣处理能力不足，一部分由于冷却风量不能太大而影响原先设计的烟风系统；一部分由于灰渣堆积，被覆盖在中心的灰渣难以将热量传递给空气。目前也有学者提出排渣废热的回收应用。黄嘉驷6等人分析了冷渣水回水节点对机组整体循环效率的影响。李慧君7等人提出采用有机朗肯循环(ORC)余热发电技术回收排渣废热。更有学者提出与煤气化行业联合8。对于排渣废热的应用研究较多，但如何在煤粉炉中实现这一过程却少有提及。本文提出一种间接水冷的排渣方式。采用刮板运输机运输灰渣颗粒，这样可以避免风冷式钢带运输机与灰渣摩擦力不够而引起的输渣能力不足问题，还可以使灰渣颗粒在运输过程中产生相对位移，更好的被冷却。在回收排渣废热的同时，不影响炉膛燃烧条件和排烟温度。2.模型构建模型构建 思路源于炉膛水冷壁，运输灰渣的设备由钢带式运输机改变为刮板式运输机，运输机的底板由类似于水冷壁结构的钢结构代替，如图 1 所示。拟采用 24 根平行管路来组成 1.32 m 宽的运输轨道，刮板间距 0.324 m，运输速度 0.046 m/s。在排渣量 5.46 t/h 的工况下，单根管子、两个刮板内的灰渣量为 0.4451 kg。灰渣物性均采用真实物性，来源于电厂技术协议和文献9 10数据插值得到，展示于表 1。Open AccessOpen Access李永华，付志 DOI:10.12677/aepe.2023.113013 103 电力与能源进展 (a)整体结构图 (b)截面结构图 Figure 1.Illustrative diagram of model structure 图图 1.模型结构说明图 Table 1.Physical properties of ash and steel structures 表表 1.灰渣和钢结构的物性表 参数(单位)灰渣 钢结构 密度(kg/m3)2300 7800 杨氏模量(Pa)1e8 2e11 泊松比 0.3 0.3 恢复系数 0.5 0.5 静摩擦系数 0.9 0.637 滑动摩擦系数 0.04 0.04 形状 球形 直径(mm)1 粒径分布 单一粒径 初始温度(K)523.15 313.15 比热容(J/(kgK)850.4 574 导热系数(W/(mK)1 40.4 3.数学模型数学模型及网格无关验证及网格无关验证 3.1.数学模型数学模型 3.1.1.颗粒运动及换热计算模型颗粒运动及换热计算模型 EDEM 采用离散元的数值模拟方法，灰渣颗粒之间的作用采用 Hertz-Mindlin(no slip)模型来计算。在这个模型中，法向力成分是基于赫兹接触理论。切向力模型是基于 Mindlin-Deresiewicz 的工作。法向力和切向力都有阻尼成分，其中阻尼系数与恢复系数有关，切向摩擦力遵循库仑摩擦定律模型，滚动摩擦力以独立于接触的方向性恒定扭矩模型实现。具体计算公式如下：法向力是重合度的函数：3*243nnFEr=(1)式中n是法线重合度；E*是有效杨氏模量；r*是平均粒径，定义分别如下：()()22*111ijijwwEEE=+(2)李永华，付志 DOI:10.12677/aepe.2023.113013 104 电力与能源进展 *111ijRRr=+(3)式中 E，w，R 分别为杨氏模量、泊松比和半径。另外法线阻尼力通过下式给出：*526drelnnnFS m v=(4)式中 m*是等效质量；relnv 是相对速度的法线分量，m*和和 Sn的定义如下：1*111immm=+(5)22lnlnee=+(6)*2nnSER=(7)式中 e 是恢复系数。切向力 Ft取决于切向重叠度t和切线刚度 St：tttFS=(8)*8tnSGR=(9)式中 G*是等效剪切模量，切向阻尼由下式给出：*526dreltttFS m v=(10)式中reltv 是相对切向速度，切向力受摩擦力的限制，最后，颗粒总的受力用合力计算：ntFFF=+(11)颗粒滚动的计算，滚动摩擦的计算是通过在壁面加扭矩来体现的，计算公式如下：irniiF R=(12)式中r是滚动摩擦系数，R 是接触点与质心的距离，i是接触点上物体的单位角速度矢量。颗粒的温度会随热通量的计算而改变从而实现温度的传递，定义如下：1 21 2p pcp pQhT=(13)接触面积的作用表现在传热系数中，并由下列公式给出：1/3*12*12434ppNcppk kF rhkkE=+(14)式中 FN是法向力，N。颗粒的温度更新通过以下公式计算：ddppheatTm CQt=(15)3.1.2.等效介质假设法计算模型等效介质假设法计算模型 整个运输过程的换热通过 Fluent 软件实现，主要方程如下 连续性方程：()0t+=v (16)李永华，付志 DOI:10.12677/aepe.2023.113013 105 电力与能源进展 式中代表密度，kg/m3；t 代表时间，s；v是速度矢量，m/s。动量方程：()()()pgt+=+vvv (17)()23I=+vvv (18)式中 p 代表压力，Pa；g 代表重力加速度 m/s2；代表应力张量；是粘度，kg/(ms)。能量方程：()()ieffhijjTHuHpkStxxx+=+(19)Prptefftckk=+(20)式中 H 代表总能，keff表示湍流有效导热系数，W/(mk)；T 表示温度，K；Sh代表源项，t表示湍流粘度，kg/(ms)；Prt表示湍流普朗特数。本文使用 Realizable k-模型进行计算，Y+的范围是 0300，选择增强型壁面处理近壁函数。相关计算公式如下：()()tjkjjkjkkkuGtxxx+=+(21)()()212tjjjjuC SCtxxxkv+=+(22)式中 k 表示湍流动能，m2/s2；表示湍流动能耗散率，m3/s2；Gk表示由于平均速度梯度而产生的湍流动能；S 表示平均应变速率的模数；C2 是常数，其值为 1.9；k和e是 k 和 e 的经验值，其值分别为 1.0和 1.2。对于固体域，只计算能量方程，因为连续性方程和动量方程自动满足。能量方程如下：()()()thhk TSt+=+v (26)式中 h 表示显焓，J/kg；St表示体积热源项，W/m3；左面第二项是对流项，是由于固体移动产生的。在稳态条件下运行时，灰渣温度可以看作不随时间变化，所以模拟采用伪瞬态 Couple 算法。压力使用 PRESTO 格式，动量、能量和湍流耗散率使用 QUICK 格式。二阶迎风格式用于湍流动能。为确保计算结果收敛，监测底灰的温度变化，直到稳定。3.2.网格无关性验证网格无关性验证 网格无关性验证是做模拟不可缺少的一步，根据软件不同验证的需求程度也不同。单纯 EDEM 软件的模拟计算不需要进行网格验证，该软件的网格用来检索灰渣颗粒及其接触状态，一般取颗粒半径的 23倍，本文取 3 倍颗粒半径，时间步长取 1.64e6。所以本节针对使用 Fluent 模拟换热过程的间接水冷结构进行网格无关性验证，如图 2。当网格数量大于 42 万时，冷却水的出口温度和排渣温度不再大幅度变化，所以本文选取 59 万网格作为模拟使用的网格。李永华，付志 DOI:10.12677/aepe.2023.113013 106 电力与能源进展 Figure 2.Pipe heat exchange grid irrelevance verification diagram 图图 2.管道换热网格无关性验证图 4.结果与讨论结果与讨论 4.1.运动过程分析运动过程分析 颗粒运动模拟参考滚筒冷渣器内颗粒运动研究11 12。钢带式运输机和刮板式运输机颗粒运动情况如图 3。Figure 3.Ash particle movement diagram(steel belt conveyor(left)scraper conveyor(right)图图 3.灰渣颗粒运动图(钢带式运输机(左)刮板式运输机(右)前 2.2 s 灰渣掉落堆积，运输机不工作，两种运输方式的颗粒状态相同，2.2 s 颗粒生成完成，运输机开始运行，直到 15.55 s 模拟结束。通过钢带式运输机运输灰渣时，灰渣颗粒之间没有相对位移，一直保持一定的厚度移动；通过刮板运输机运输灰渣的过程中，2.2 s 到 6.5 s 之间，随着刮板的移动，右侧颗粒失去了刮板的支撑，坍落成自然堆积状态，左侧颗粒在刮板的推动作用下，逐渐堆积，颗粒从平铺状态变成左高右低的状态，到 6.5 s 以后，灰渣堆积成固定形状，以固定的堆积高度和堆积角度随着刮板一起向前移动，堆积高度约为 74.10 mm，堆积角度约为 26.92。6.5 s 以后刮板运输机中灰渣运动状态看似规律，但颗粒运动十分复杂。李永华，付志 DOI:10.12677/aepe.2023.113013 107 电力与能源进展 紧靠刮板和底板的区域，这里的颗粒被挤压逐渐向上移动，而灰渣之间的摩擦小于灰渣与两侧挡板的摩擦，加上重力作用，导致这些颗粒向中间靠拢；位于最高位置的颗粒有的保持高度随刮板直线运动，有些被挤压到斜面区域然后滚落到底板。图 4 和图 5 分别显示了刮板式运输机两个挡板中心截面和挡板附近的颗粒运动状态，速度高于 0.08 m/s 的颗粒用红色表示，小于 0.02 m/s 的颗粒用蓝色表示。中心截面的颗粒明显分成三个区域，底板附近的颗粒移动速度小于刮板移动速度，可以看成类似静止区，斜面部位红色区域是滚落区，还有中间和刮板运行速度相差不大的移动区。挡板附近颗粒移动十分混乱，在底板附近有更大范围的类似静止区，在颗粒内部不仅出现了许多速度小于 0.02 m/s 的区域，还出现了个别速度大于 0.08 m/s 的颗粒，这是颗粒受重力，挤压，摩擦共同作用的结果。从 6.5 s 到 15.55 s，颗粒保持这种状态向前移动。Figure 4.Movement of particles in the central section 图图 4.中心截面颗粒的运动状态 Figure 5.Movement of particles in the cross section near the baffle 图图 5.挡板附近截面颗粒的运动状态 李永华，付志 DOI:10.12677/aepe.2023.113013 108 电力与能源进展 4.2.换热分析换热分析 钢带式运输机中灰渣的冷却过程及局部放大图如图 6，最上层灰渣得不到冷却，温度难以传递到底板，所以排渣过程中灰渣的最高温度保持初始温度不变，最下层灰渣一直受凝结水的间接冷却，最低温度随时间一直降低，从 523.15 K 降低到 431.74 K。整体温度分布呈现上高下低的趋势，类似于非稳态导热温度分布。刮板式运输机中灰渣的冷却过程及不同视角局部放大图如图 7。相比于钢带式运输机，刮板式运输机使灰渣之间相互混合，颗粒的运动轨迹不是一条水平线，灰渣颗粒之间的温差更小、温度更均匀，最大温差只有 52.56 K。灰渣颗粒最高温度同样一直不变，一部分是由于模拟时间短；一部分由于颗粒混合效果并不均匀，刮板不能保证所有灰渣颗粒都能够与底板接触。最低温度呈现先降低后升高的趋势，到 11.3 s 时灰渣颗粒达到最低温度，之后颗粒最低温度上升，说明 11.3 s 后，一开始被堆积在下面的灰渣颗粒被挤压到更高的位置，这些灰渣颗粒开始得不到底板的冷却，开始受周围更高温度的灰渣颗粒加热。Figure 6.Temperature distribution of ash particles in steel belt conveyor 图图 6.钢带式运输机灰渣颗粒温度分布图 Figure 7.Temperature distribution of ash particles in scraper conveyor 图图 7.刮板式运输机灰渣颗粒温度分布图 李永华，付志 DOI:10.12677/aepe.2023.113013 109 电力与能源进展 表 2 为刮板式运输机和钢带式运输机排渣冷却效果对比，同样冷却时长，钢带运输机降温 2.35 K，刮板运输机降温 1.51 K，降温效果降低 35.74%。这是由于刮板运输形式使灰渣与底板的接触面积降低，根据模拟结果，刮板式运输机灰渣投影到底板的面积为 90.20 cm2，仅为钢带式运输机灰渣投影到底板面积的 50.62%。数值上来看，刮板运输机在同样换热面积下，具备更好的冷却条件，如果采用合适的排渣方式弥补这一缺陷，刮板运输机将会有更好的冷却效果。刮板推动灰渣运动可以使中心高温颗粒接触冷却底板，增加换热温差，有类似破坏温度边界层的功能。Table 2.Comparison table of cooling effect 表表 2.降温效果对比表 运输机 时刻(s)颗粒平均温度(K)平均温降(K)刮板运输机 2.2 522.76,1.51 15.55 521.25 钢带运输机 2.2 522.76 2.35 15.55 520.41 整体排渣过程的冷却通过 Fluent 软件模拟，进口位置的灰渣温度为某热电厂现场测量温度，通过等效介质假设模拟灰渣在整个运输阶段的降温，灰渣的堆积参数及其他模拟条件如表 3 所示。Table 3.Physical properties parameters 表表 3.物性参数 材料 参数(单位)数值 底灰 密度(kg/m3)900 初始温度(K)523.15 质量流量(kg/s)2.5277 导热系数(W/(m2K)13 0.10420.0003t=+比热容(J/(kgK)13 572.87870.502ct=+斯蒂芬玻尔兹曼常数(W/(m2K4)5.67 108 发射率14 0.94 运输速度(m/s)0.046 水 密度(kg/m3)998 初始温度(K)313.15 导热系数(W/(m2K)0.6 比热容(J/(kgK)4200 流速(m/s)0.1,0.3,0.5,0.7,1 粘度(kg/(ms)0.001 在仅考虑锅炉给水的间接冷却条件下，保守地估计新结构对灰渣的降温效果。同样的冷却时间下，灰渣可以最低被冷却到 408.45 K，灰渣温度降低了 114.7 K，对比校验后应用刮板式运输机可以使灰渣降温 73.25 K，所以排渣温度约为 449.9 K，满足冷却灰渣的要求。图 8 显示了冷却水流速 0.3 m/s 时钢结构的截面随轴向的温度变化情况。Z 变大方向为灰渣的移动方李永华，付志 DOI:10.12677/aepe.2023.113013 110 电力与能源进展 向。钢结构内部的温度分布呈对称趋势，与水冷壁温度分布相似，距离冷却水越近的地方，温度越低。最低温度一直出现在钢结构的最底部，距离冷却水最远的地方温度最高。随着灰渣的移动，钢结构截面的温差越来越小，最后接近冷却水进口温度。Z=0.2 m 处的钢结构的最高温度和最低温度差距最大，因为此处灰渣与底板换热最强烈。Figure 8.Indirect water-cooled steel structure temperature distribution with position 图图 8.间接水冷钢结构温度分布随位置变化图 图 9 说明了冷却水流量与排渣温度的关系，当冷却水流速高于 0.5 m/s 时，增加冷却水流量不会对排渣温度造成明显影响。值得注意的是冷却水流速等于 0.1 m/s 时，冷却水出口温度达到了 360.95 K。这说明排渣可以作为加热厂内生活用水和热网回水的热量来源。这也说明排渣废热有产生蒸汽的能力，可以与有机朗肯循环结合以产生更高品位的能量。Figure 9.Indirect water-cooled cooling effect with cooling water flow rate 图图 9.间接水冷冷却效果随冷却水流速变化图 5.结论结论 1)刮板运输机可以有效的使中心灰渣与冷源接触，排渣温度更均匀，不会产生钢带式运输中出现的温度分层现象，但由于其工作特点，导致灰渣与底板的接触面积减少了 49.38%，所以整体冷却效果不及同冷却条件下的钢带式运输机，在运输段进口温度 523.15 K 的情况下，可以降温 73.25 K，排渣温度约为李永华，付志 DOI:10.12677/aepe.2023.113013 111 电力与能源进展 449.9 K，基本满足运输过程中排渣的冷却要求。2)冷却水管的温度分布各截面趋势相同，底部温度低，肋片处温度高，延灰渣移动方向，冷却水管截面温差越来越小，越来越接近冷却水温度。基金项目基金项目 国家自然科学基金项目(52072081)。参考文献参考文献 1 童家麟,齐晓娟,韩平,等.炉底漏风对某 600 MW 机组锅炉性能的影响J.热力发电,2021,50(4):126-31.2 刘铭媛,陈国庆,戴维葆,等.炉底漏风对炉内燃烧及NOx生成排放特性的影响J.机械工程学报,2019,55(24):210-215+226.3 王轶锋,高飞燕.发电厂干式排渣系统对锅炉效率的影响试验及分析J.内蒙古电力技术,2010,28(2):52-53+56.4 曹旭东.电站锅炉干式排渣系统的选型设计及其对锅炉效率的影响D:硕士学位论文.北京:华北电力(北京)大学,2002.5 何彦鹏,刘宏举,郝峰.干排渣系统降氮节能改造技术研究及应用J.上海节能,2021(6):619-623.6 黄嘉驷,张爱军,魏铜生,等.冷渣器冷却水回水接入点分析及改进方案J.热力发电,2009,38(3):75-77+80.7 李慧君,谷凯娜,马超.机组余热利用的热经济性分析J.热力发电,2014,43(1):6-10.8 Barati,M.,Esfahani,S.and Utigard,T.A.(2011)Energy Recovery from High Temperature Slags.Energy,36,5440-5449.https:/doi.org/10.1016/j.energy.2011.07.007 9 司小东.滚筒冷渣器传热模型研究D:硕士学位论文.北京:清华大学,2011.10 张超.超临界直流锅炉垂直管屏水冷壁的热应力分析J.广东电力,2022,35(7):133-140.11 彭涛,何庆中,王佳,等.组合式滚筒冷渣机中颗粒径向扩散运动研究J.热力发电,2019,48(12):69-74.12 何庆中,赖镜安,王佳,等.基于 EDEM-FLUENT 耦合的滚筒冷渣机颗粒轴向扩散运动J.科学技术与工程,2022,22(3):1004-1010.13 石磊.干式排渣系统传热特性研究及可视化软件开发D:硕士学位论文.北京:华北电力大学,2017.14 Qiu,L.,Sang,D.W.,Li,Y.L.,et al.(2020)Numerical Simulation of Gas-Solid Heat Transfer Characteristics of Porous Structure Composed of High-Temperature Particles in Moving Bed.Applied Thermal Engineering,181,Article ID:115925.https:/doi.org/10.1016/j.applthermaleng.2020.115925