未燃煤粉在高炉内可参加铁水渗碳、碳气化反应和还原反应, 但随着高炉喷煤量的提高, 若未燃煤粉堆积量超过了高炉内部反应消耗的量, 未燃煤粉将在高炉内部堆积, 影响高炉透气性, 使压差升高, 严重时影响高炉顺行; 因此未燃煤粉在高炉内的过量堆积是限制高炉喷煤量的重要因素之一[1-3].
为了改善高炉炉况、降低未燃煤粉对高炉的负面影响, 有必要对未燃煤粉在高炉内的分布进行研究.国内外学者对高炉内未燃煤粉的运动分布作了相应的研究.文献[4-5]分别通过物理模拟的方法研究了喷吹一定量煤粉后, 煤粉在高炉模型内的分布及其对压差的影响; Dong等[6]通过填充床模型考察了软熔带对高炉内气体和粉体流动的影响; Nogami等[7]研究了死料柱结构和粉体直径对高炉内气体和粉体流动的影响.然而这些研究并未考虑高炉内各区域空隙度的变化.本文将高炉分为若干空隙率不同的区域, 考察高炉模型内未燃煤粉的堆积分布状况.气固两相流通常采用“Euler-Lagrange”法和“Euler-Euler”法.本文计算未燃煤粉的堆积, 由于粉体体积分数较高, 故采用“Euler-Euler”法建立数学模型, 对未燃煤粉在高炉内的堆积分布进行研究, 考察了不同软熔带形状下高炉内未燃煤粉的堆积分布.
1 数学模型高炉内的过程非常复杂, 对高炉内的气体和煤粉运动进行数值模拟时, 视高炉为散料体填充床, 气体和未燃煤粉为互相贯穿的连续介质[8]; 并作如下假设:①气体、煤粉以均匀速度从风口吹入; ②气体为不可压缩的稳态流动; ③忽略传热传质及化学反应; ④填充颗粒为直径相等、球形度相同的颗粒; ⑤忽略高炉块状带的空隙度变化.
高炉内的未燃煤粉根据其运动状态可分为两部分[9], 一部分为被气体带出炉外的运动煤粉, 另一部分为由于阻力滞留在高炉内部的静止煤粉.Hidaka等[10]对填充床中静止粉体的分布进行了试验研究, 得到了静止粉体体积分数的经验公式, 当气体速度接近于0时, 该公式不能描述高炉内堆积的未燃煤粉实际情况; 因此对该公式作了修正, 修正后的公式如下:
(1) |
式中:φfs为静止煤粉的体积分数; φfmax为未燃煤粉的体积分数最大值, 本文取值为0.24[8]; vg为气体的表观速度, m/s.
根据未燃煤粉在高炉内的流动行为, 高炉内气体、未燃煤粉和填充颗粒的体积分数满足式(2):
(2) |
式中:φg为气体的体积分数; φf为未燃煤粉的体积分数,φf=φfd+φfs; φfd为运动煤粉的体积分数; φp为填充颗粒的体积分数.
计算所采用的连续性方程和动量守恒方程如下:
(3) |
(4) |
式中:ρg, ρf分别为气体、未燃煤粉的密度, kg·m-3; ug, ufd为气体、运动煤粉的速度, m·s-1.
(5) |
(6) |
式中:τg, τf为气体、未燃煤粉的黏性应力; p为气体和未燃煤粉共享的压力, Pa; g为重力加速度, m·s-2; Fgf, Ffg为气体与未燃煤粉间的相互作用力,Fgp为气体与填充颗粒间的相互作用力,Ffp为未燃煤粉与填充颗粒之间的相互作用力, kg·m-2·s-2.
湍流方程为κ-ε方程, 详细方程见文献[11].
1.1 计算模型计算模型为某1 250 m3高炉, 其炉膛截面结构尺寸如图 1所示, 计算区域为二维规则几何体; 为减少计算量, 选取1/2高炉模型进行网格划分, 网格类型为四边形结构性网格.
各边界条件的设置如下:1)炉壁和炉底设置为固体壁面, 采用无滑移边界条件, 各相速度为零; 2)入口采用速度入口边界条件, 在入口处各相的速度由各相流量确定, 速度方向为垂直于边界面; 3)出口采用压力出口边界条件, 出口处各相梯度均为零; 4)对称轴设置为对称边界条件, 在对称边界上, 垂直边界的速度为零, 其他物理量的值在该边界内外是相等的.近壁面采用壁面函数法.
1.3 软熔带参数模型忽略软熔带矿石层中的煤气流量[12], 将软熔带矿石层设定为气体和未燃煤粉不能通过的区域.软熔带参数的定义如图 1所示.倒V, V和W型三种高炉软熔带的参数见表 1.
1) 高炉各区域及空隙度参数.将高炉模型分为块状带、软熔带“焦窗”、滴落带、死料柱和风口回旋区5个区域(见图 2).根据文献[13], 从料面到风口间的高炉平均空隙度在0.231~0.426, 高炉块状带的空隙度在0.35~0.46, 因此将高炉块状带、软熔带“焦窗”、滴落带、死料柱和风口回旋区的空隙度ε1, ε2, ε3, ε4,ε5分别设为0.4, 0.35, 0.3, 0.25, 1.0.
将高炉模型从炉喉沿径向分为三等分(如图 2所示), 即高炉模型分为三个区域:中心区域、中间区域和边缘区域, 以便于比较未燃煤粉在高炉模型各区域的堆积分布.
2) 物性参数.高炉模型中的填充颗粒球形度均设为0.75, 颗粒直径为0.025 m; 气体为1.01×105 Pa下20 ℃的空气; 烟煤和无烟煤的密度分别为1 250~1 500 kg·m-3, 1 360~1 800 kg·m-3[1]; 本研究选取未燃煤粉密度为1 400 kg·m-3, 气体和煤粉的物性参数及模拟高炉的基本操作参数见表 2.
采用1:25比例的高炉二维模型, 通过比较中心线压差来验证数学模型的正确性, 高炉模型几何尺寸及测压点如图 3a所示.
由图 3b可以看出, 计算值与实验值的趋势基本吻合.因此数学模型能较好地模拟高炉填充床模型中气固两相流的运动.
2 计算结果与讨论 2.1 不同软熔带形状下未燃煤粉的分布软熔带的形状是高炉上下部调剂的综合体现, 主要受布料制度和送风制度的影响, 是决定高炉煤气流稳定运行的关键.在相同喷吹条件下, 三种软熔带形状下高炉内未燃煤粉的堆积分布如图 4所示.由图可知:1)高炉块状带的未燃煤粉堆积量较少, 未燃煤粉主要堆积在高炉软熔带下方; 死料柱、风口回旋区下方、软熔区域是未燃煤粉易堆积的区域; 2)软熔带形状对软熔区域及死料柱上的未燃煤粉的分布有重要影响:倒V型软熔带高炉中, 未燃煤粉易在软熔区域的中心堆积; V型软熔带高炉中, 未燃煤粉在软熔区域的边缘容易堆积, 在死料柱上方未燃煤粉的堆积量较多; W型软熔带高炉中, 未燃煤粉在高炉软熔区域的中心和边缘的堆积比较分散, 显著堆积的区域较少, 未燃煤粉在软熔区域的分布相对均匀.
图 5为高炉内运动煤粉的分布图, 由图可见:运动煤粉主要在高炉中间区域运动; 软熔带改变了未燃煤粉的运动方向, 对运动煤粉起到了二次分布的作用.在倒V型软熔带和W型软熔带高炉内, 运动煤粉经过软熔带后流动方向趋向边缘; 在软熔带高炉内, 运动粉体经过软熔带后流动方向趋于高炉中心.运动煤粉通过软熔带“焦窗”向上流动, 在V型和倒V型高炉块状带区, 运动煤粉汇聚成细条状; 在W型高炉块状带内, 运动煤粉汇聚现象有所降低, 说明W型高炉块状带内运动煤粉分布较为均匀, 有利于该区域的气体组成和温度均匀分布.
为比较不同软熔带高炉内未燃煤粉堆积量的多少, 定义了高炉内未燃煤粉的平均体积分数φf, 平均体积分数越大, 则高炉模型内未燃煤粉累积越多, 反之亦然.平均体积分数φf定义如下:
(7) |
式中,V为模型计算域体积, m3.
三种软熔带高炉内未燃煤粉的平均体积分数如图 6所示.由图可见, 相同喷吹条件下, V型软熔带高炉内堆积的未燃煤粉最多, 其次是W型软熔带高炉, 倒V型软熔带高炉堆积的未燃煤粉最少.从容纳未燃煤粉量的多少考虑, V型和W型软熔带高炉是较好的选择.但在实际生产中, V型软熔带高炉在高炉中心区域堆积的未燃煤粉较多, 不利于发展中间气流, 生产指标不好, 对高炉顺行不利.因此从容纳未燃煤粉量和高炉实际操作的角度考虑, W型软熔带在一定全压差条件下可容纳较多煤粉, 且分布较均匀, 有利于高炉顺行, 这是W型软熔带的高炉操作优势之一.
将高炉沿垂直高度划分为若干层, 每层分为中心、中间和边缘三个区域, 通过平均体积分数比较高炉模型各部位的未燃煤粉堆积量.图 7为三种软熔带高炉分层区域的未燃煤粉平均体积分数.由图可见, 三种软熔带高炉的块状带和炉缸区域未燃煤粉平均体积分数基本一致,即未燃煤粉堆积量基本一致.在软熔带及其根部和上部区域, 未燃煤粉平均体积分数相差较大, 说明软熔带形状对该区域未燃煤粉的分布影响较大.在三种软熔带高炉的中心和中间区域, 倒V型高炉堆积的未燃煤粉量较多, 其次是W型和V型高炉; 而在边缘区域, 与之相反.因此比较三种软熔带高炉, W型软熔带高炉的未燃煤粉在软熔带区域径向分布相对均匀, 这也是W型软熔带的高炉操作优势之一.
图 8为三种软熔带高炉在高炉中心线上的压差分布(以高炉出口压力101 325 Pa为基准压力, 以风口水平线为起点).由图可见:1)在高炉块状带区域, 三种软熔带高炉的压差分布无显著差别; 2)高炉软熔带的形状对高炉中心压差的影响主要在软熔带区域, V型软熔带压力分布不均, 这是高炉操作不选择V型软熔带的原因之一.
1) 软熔带形状对高炉内未燃煤粉分布的影响主要在软熔带及其根部和上部区域.未燃煤粉主要堆积在高炉软熔带下方, 死料柱、风口回旋区下方、软熔区域是未燃煤粉最容易堆积的区域.
2) 三种软熔带中, W型软熔带是最佳的高炉操作模式, 具体表现在:①未燃煤粉在高炉内分布相对均匀; ②在一定压差条件下可容纳较多煤粉.
[1] |
Shen Y S, Yu A B, Peter A, et al.
Modelling in-furnace phenomena of pulverized coal injection in ironmaking blast furnace:effect of coke bed porosities[J]. Minerals Engineering, 2012, 33: 54–65.
DOI:10.1016/j.mineng.2011.10.014 |
[2] |
Nomura S, Callcott T G.
Maximum rates of pulverized coal injection in ironmaking blast furnaces[J]. ISIJ International, 2011, 51(7): 1033–1043.
DOI:10.2355/isijinternational.51.1033 |
[3] |
The Technical Society, the Iron and Steel Institute of Japan.
Production and technology of iron and steel in Japan during 2015[J]. ISIJ International, 2016, 56(6): 905–916.
DOI:10.2355/isijinternational.56.905 |
[4] |
刘新, 陈星秋.
未燃煤粉在高炉内的分布特性的实验研究[J]. 东北大学学报(自然科学版), 2000, 21(2): 177–179.
( Liu Xin, Chen Xing-qiu. Experimental investigation on behavior of UPC (unburnt pulverized coal)in blast furnace[J]. Journal of Northeastern University(Natural Science), 2000, 21(2): 177–179. DOI:10.3321/j.issn:1005-3026.2000.02.018 ) |
[5] |
李洪会, 孙亮, 沈峰满, 等.
未燃煤粉在填充床内的分布及对料柱压差的影响[J]. 东北大学学报(自然科学版), 2001, 22(3): 303–305.
( Li Hong-hui, Sun Liang, Shen Feng-man, et al. Distribution of unburned pulverized coal and effects on pressure drop of gas in packed bed[J]. Journal of Northeastern University(Natural Science), 2001, 22(3): 303–305. DOI:10.3321/j.issn:1005-3026.2001.03.019 ) |
[6] |
Dong X F, Zhang S J, Pinson D, et al.
Gas-powder flow and powder accumulation in a packed bed:Ⅰ.experimental study[J]. Powder Technology, 2004, 149(1): 1–9.
DOI:10.1016/j.powtec.2004.09.040 |
[7] |
Nogami H, Austin P R, Yagi J, et al.
Numerical investigation on effects of deadman structure and powder properties on gas and powder flows in lower part of blast furnace[J]. ISIJ International, 2004, 44(3): 500–509.
DOI:10.2355/isijinternational.44.500 |
[8] |
Dong X F, Zhang S J, Pinson D, et al.
Gas-powder flow and powder accumulation in a packed bed:Ⅱ.numerical study[J]. Powder Technology, 2004, 149(1): 10–22.
DOI:10.1016/j.powtec.2004.09.039 |
[9] |
Pintowantoro S, Nogami H, Yagi J.
Numerical analysis of static holdup of fine particles in blast furnace[J]. ISIJ International, 2004, 44(2): 304–309.
DOI:10.2355/isijinternational.44.304 |
[10] |
Hidaka N, Matsumoto T, Kusakabe K, et al.
Flow of fine particles entrained by upward gas flow in packed beds of coarse particles[J]. Journal of Chemical Engineering of Japan, 1999, 32(2): 197–203.
DOI:10.1252/jcej.32.197 |
[11] |
丁智敏, 姜鑫, 张枥, 等.未燃煤粉在高炉中堆积分布的数值模拟[C]//2016年全国高炉炼铁学术年会论文集.邯郸: 中国金属学会, 2016: 576-583.
( Ding Zhi-min, Jiang Xin, Zhang Li, et al.Numerical simulation on the distribution of unburned pulverized coal in blast furnace[C] //The Conference of Blast Furnace Ironmaking of China in 2016.Handan: The Chinese Society for Metals, 2016: 576-583. ) |
[12] |
于海彬, 陈义胜.
软熔带对煤气流动影响的数值模拟[J]. 包头钢铁学院学报, 2006, 25(1): 1–4.
( Yu Hai-bin, Chen Yi-sheng. Numerical simulation of influence on gas flow caused by cohesive zone[J]. Journal of Baotou University of Iron and Steel Technology, 2006, 25(1): 1–4. DOI:10.3969/j.issn.2095-2295.2006.01.001 ) |
[13] |
周传典.
高炉炼铁生产技术手册[M]. 北京: 冶金工业出版社, 2008: 129.
( Zhou Chuan-dian. Technical manuals of blast furnace ironmaking production[M]. Beijing: Metallurgical Industry Press, 2008: 129. ) |
[14] |
解维伟.
煤化学与煤质分析[M]. 北京: 冶金工业出版社, 2012: 116.
( Xie Wei-wei. Coal chemistry and coal quality analysis[M]. Beijing: Metallurgical Industry Press, 2012: 116. ) |