高炉煤粉燃烧率通用模型

第36卷第2期 北京科技大学学报 Vol.36 No.2 2014年2月 Journal of University of Science and Technology Beijing Feb.2014 高炉煤粉燃烧率通用模型 陈川区，程树森 北京科技大学治金与生态工程学院，北京100083 ☒通信作者，E-mail:chenchuan(606@l63.com 摘要建立了煤粉燃烧率通用模型，模型可以根据煤粉的工业分析值计算燃烧动力学参数并预测煤粉燃烧率.通过对比前 人的实验数据，验证了模型的准确性，同时研究了影响高炉煤粉燃烧率的若干因素。研究结果表明：在高炉喷煤过程中，煤粉 颗粒在2ms左右就可以达到热风速度，由于煤粉颗粒在直吹管内停留时间短并且温度较低，因此在直吹管内煤粉不会发生燃 烧.煤粉进入风口回旋区后，挥发分瞬间全部析出，并且颗粒粒径越小，挥发分开始析出时间越早，降低煤粉粒径和增加氧气 体积分数均有利于提高煤粉燃烧率.氧气体积分数每增加1%，燃烧率提高2%.随着喷煤量的增加，煤粉燃烧率逐渐降低 当提高煤粉喷吹量时，为了保证较高的燃烧率，实际操作过程中应提高富氧率并适当降低煤粉粒径 关键词高炉：煤粉燃烧：数学模型：因素分析 分类号TF531 General model of pulverized coal combustion ratio in blast furnaces CHEN Chuan,CHENG Shu-sen School of Metallurgical and Ecological Engineering,University of Science and Technology Beijing,Beijing 100083,China Corresponding author,E-mail:chenchuan606@163.com ABSTRACT A general model of pulverized coal combustion ratio was put forward,which can calculate the kinetic parameters and combustion ratio based on coal industrial analysis.The accuracy of this model was verified by comparing with previous experimental data.Some factors influencing the coal combustion ratio were discussed.It is found that during the process of pulverized coal injection, the pulverized coal particle can reach the velocity of hot blast within 2 ms.In the blowpipe,due to the pulverized coal particle with short residence time and lower temperature,combustion reaction does not occur.When the pulverized coal particle enters into the race- way,the volatile is released instantly.The smaller the coal diameter is,the earlier the volatile will be released.The coal combustion ratio can be improved by decreasing the coal diameter or increasing the oxygen enrichment.When the oxygen enrichment increases 1% the combustion ratio increases about 2%.With increasing coal injection rate,the coal combustion ratio decreases.Therefore,when the coal injection rate increases,we should increase the oxygen enrichment and decrease the particle diameter appropriately to maintain the high combustion ratio. KEY WORDS blast furnaces:coal combustion:mathematical models:factor analysis 随着高炉大型化的发展，焦炭已经不能满足高 炉况不顺，影响高炉的正常生产.因此，提高高炉内 炉“高产、低耗”的生产要求，通过喷吹煤粉置换部 煤粉的燃烧率，对促进喷吹，提高煤比，保证高炉的 分焦炭己成为炼铁节能降耗的重要技术措施，而煤 稳定顺行具有不可忽视的作用B-) 粉在风口前的燃烧率是限制煤粉喷吹量的重要因 回旋区内煤粉的燃烧是在有限空间内的高温、 素口.风口前煤粉燃烧率过低将会导致大量的未燃 高压、高速流动问题和颗粒燃烧反应的多相流问题 煤粉，从而降低煤粉燃烧率，影响煤焦置换比可.同 的结合。煤粉燃烧速率的快慢是衡量煤粉燃烧性能 时，未燃煤粉过多还会恶化炉缸透气透液性，进而使 的一个重要参数.目前有经验的因、半经验的因、数 收稿日期：2012-11-17 基金项目：国家自然科学基金资助项目(61271303) DOI:10.13374/j.issn1001-053x.2014.02.014:http://journals.ustb.edu.cn

第 36 卷 第 2 期 2014 年 2 月 北京科技大学学报 Journal of University of Science and Technology Beijing Vol． 36 No． 2 Feb． 2014 高炉煤粉燃烧率通用模型 陈 川，程树森 北京科技大学冶金与生态工程学院，北京 100083  通信作者，E-mail: chenchuan606@ 163． com 摘 要 建立了煤粉燃烧率通用模型，模型可以根据煤粉的工业分析值计算燃烧动力学参数并预测煤粉燃烧率． 通过对比前 人的实验数据，验证了模型的准确性，同时研究了影响高炉煤粉燃烧率的若干因素． 研究结果表明: 在高炉喷煤过程中，煤粉 颗粒在 2 ms 左右就可以达到热风速度，由于煤粉颗粒在直吹管内停留时间短并且温度较低，因此在直吹管内煤粉不会发生燃 烧． 煤粉进入风口回旋区后，挥发分瞬间全部析出，并且颗粒粒径越小，挥发分开始析出时间越早． 降低煤粉粒径和增加氧气 体积分数均有利于提高煤粉燃烧率． 氧气体积分数每增加 1% ，燃烧率提高 2% ． 随着喷煤量的增加，煤粉燃烧率逐渐降低． 当提高煤粉喷吹量时，为了保证较高的燃烧率，实际操作过程中应提高富氧率并适当降低煤粉粒径． 关键词 高炉; 煤粉燃烧; 数学模型; 因素分析 分类号 TF 531 General model of pulverized coal combustion ratio in blast furnaces CHEN Chuan ，CHENG Shu-sen School of Metallurgical and Ecological Engineering，University of Science and Technology Beijing，Beijing 100083，China  Corresponding author，E-mail: chenchuan606@ 163． com ABSTＲACT A general model of pulverized coal combustion ratio was put forward，which can calculate the kinetic parameters and combustion ratio based on coal industrial analysis． The accuracy of this model was verified by comparing with previous experimental data． Some factors influencing the coal combustion ratio were discussed． It is found that during the process of pulverized coal injection， the pulverized coal particle can reach the velocity of hot blast within 2 ms． In the blowpipe，due to the pulverized coal particle with short residence time and lower temperature，combustion reaction does not occur． When the pulverized coal particle enters into the race￾way，the volatile is released instantly． The smaller the coal diameter is，the earlier the volatile will be released． The coal combustion ratio can be improved by decreasing the coal diameter or increasing the oxygen enrichment． When the oxygen enrichment increases 1% the combustion ratio increases about 2% ． With increasing coal injection rate，the coal combustion ratio decreases． Therefore，when the coal injection rate increases，we should increase the oxygen enrichment and decrease the particle diameter appropriately to maintain the high combustion ratio． KEY WOＲDS blast furnaces; coal combustion; mathematical models; factor analysis 收稿日期: 2012--11--17 基金项目: 国家自然科学基金资助项目( 61271303) DOI: 10． 13374 /j． issn1001--053x． 2014． 02． 014; http: / /journals． ustb． edu． cn 随着高炉大型化的发展，焦炭已经不能满足高 炉“高产、低耗”的生产要求，通过喷吹煤粉置换部 分焦炭已成为炼铁节能降耗的重要技术措施，而煤 粉在风口前的燃烧率是限制煤粉喷吹量的重要因 素［1］． 风口前煤粉燃烧率过低将会导致大量的未燃 煤粉，从而降低煤粉燃烧率，影响煤焦置换比［2］． 同 时，未燃煤粉过多还会恶化炉缸透气透液性，进而使 炉况不顺，影响高炉的正常生产． 因此，提高高炉内 煤粉的燃烧率，对促进喷吹，提高煤比，保证高炉的 稳定顺行具有不可忽视的作用［3--4］． 回旋区内煤粉的燃烧是在有限空间内的高温、 高压、高速流动问题和颗粒燃烧反应的多相流问题 的结合． 煤粉燃烧速率的快慢是衡量煤粉燃烧性能 的一个重要参数． 目前有经验的［5］、半经验的［6］、数

第2期 陈川等：高炉煤粉燃烧率通用模型 ·227· 值求解)等方法来计算煤粉的燃烧速率.无论用 什么方法计算，都将涉及煤粉表面氧化反应动力学 me dt 8 TCap delum-ugl (u-u), 参数问题.杨天钧等o、岑可法等m、Shen 等2)、Guo等d计算煤粉燃烧速率时，都假设煤 -gcpl-1a。-w+ me dt 粉动力学参数为常数或者仅是温度的函数，并没有 1 考虑煤粉成分对动力学参数的影响.傅维镳的通 6 wgd (p-pa). 过分析大量的实验数据，找到了煤粉燃烧的通用规 律，只需要知道煤粉的工业分析值，就可以统一计算 式中：F和F,分别为煤粉颗粒运动过程在水平方 煤粉燃烧反应动力学参数，从而解决了长期困扰人 向和竖直方向受到的空气阻力，N;F,为重力和浮力 们在研究煤焦燃烧反应动力学参数时所出现的无规 的合力，N;m。为煤粉颗粒的质量，kg;d,为颗粒直 律状况 径，mp,和p分别为碳粒和气体的密度，kgm-3;“ 本文根据煤粉燃烧通用规律计算煤粉燃烧动力 为空气速度，m·s1;u,为煤粉颗粒速度，m·sl;Ca 学参数，结合高炉喷煤工艺，建立了高炉煤粉燃烧率 为颗粒在空气中运动的阻力系数，C:=。 24 + 通用模型，采用四阶龙格库塔法对模型进行离 散，并编写了煤粉燃烧率计算软件.计算过程中考 6 +0.4,Re为雷诺数，Re=4,-,,L为 1+Re L 虑了煤粉成分对燃烧动力学参数的影响，使模型更 加符合实际情况.通过对比模型预测值与前人的实 气体动力黏度，m2s. 验数据，验证了模型的准确性，同时研究了煤粉粒 1.2粒径方程 假设煤粉燃烧速率为G,则煤粉质量方程为 径、煤粉成分、鼓风温度、氧气体积分数以及喷吹量 对煤粉燃烧率的影响，结合理论分析对如何进一步 提高煤粉燃烧率进行了探讨. dims=Gmdp dt 1煤粉燃烧模型 不考虑挥发分析出对煤粉粒径的影响，则煤粉 在燃烧过程中粒径变化方程为 回旋区内煤粉燃烧数学模型基本假设条件： (1)回旋区内气体温度、压力保持稳定：(2)煤粉比 -26 (2) dt 焦炭优先燃烧：(3)高温煤粉的燃烧，主要以 2C+02=2C0反应为主；(4)煤粉颗粒内部温度均 1.3热解方程 匀：(5)不考虑挥发分析出对煤粉粒径的影响. 煤粉的热解方程采用 1.1运动方程 根据牛顿定律，煤粉颗粒二维运动方程如下： 出=化.-0kep(-后) (3) du匹=fds' me di 本文中的热解方程主要是对已有的热解方程进 行改进.以往的热解方程中动力学参数的确定大部 dug =Fay+F: (1) me dt 分是依据经验，本文给出煤粉热解过程中动力学参 数的表达式。通过对以往实验数据的进行拟合得 山，=√m+u 到热解动力学参数k,和E的表达式： expB.132exp(9.113×10-5T)+0.02399exp(0.003663T）], T≤1600： T-2173)2 T>1600. 10×76560exp(-0.01381T)）+0.7519exp(0.001966T）],T<1850: E= 110×{30.84exp- /T-20491 630.8 T≥1850. 热解模型中V,表示挥发分最大产量，它是煤粉 热解通用模型中的重要参数，与煤种、加热条件、煤 v.a(a 粒尺寸等因素有关. 式中：'为煤的可燃基挥发分：T。为煤粉的初温

第 2 期 陈 川等: 高炉煤粉燃烧率通用模型 值求解［7--9］等方法来计算煤粉的燃烧速率． 无论用 什么方法计算，都将涉及煤粉表面氧化反应动力学 参 数 问 题． 杨 天 钧 等［10］、岑 可 法 等［11］、Shen 等［12--13］、Guo 等［14］计算煤粉燃烧速率时，都假设煤 粉动力学参数为常数或者仅是温度的函数，并没有 考虑煤粉成分对动力学参数的影响． 傅维镳［15］通 过分析大量的实验数据，找到了煤粉燃烧的通用规 律，只需要知道煤粉的工业分析值，就可以统一计算 煤粉燃烧反应动力学参数，从而解决了长期困扰人 们在研究煤焦燃烧反应动力学参数时所出现的无规 律状况． 本文根据煤粉燃烧通用规律计算煤粉燃烧动力 学参数，结合高炉喷煤工艺，建立了高炉煤粉燃烧率 通用模型，采用四阶龙格库塔法［16］对模型进行离 散，并编写了煤粉燃烧率计算软件． 计算过程中考 虑了煤粉成分对燃烧动力学参数的影响，使模型更 加符合实际情况． 通过对比模型预测值与前人的实 验数据，验证了模型的准确性，同时研究了煤粉粒 径、煤粉成分、鼓风温度、氧气体积分数以及喷吹量 对煤粉燃烧率的影响，结合理论分析对如何进一步 提高煤粉燃烧率进行了探讨． 1 煤粉燃烧模型 回旋区内煤粉燃烧数学模型基本假设条件: ( 1) 回旋区内气体温度、压力保持稳定; ( 2) 煤粉比 焦炭 优 先 燃 烧; ( 3 ) 高温煤粉的燃烧，主 要 以 2C + O2 = 2CO反应为主; ( 4) 煤粉颗粒内部温度均 匀; ( 5) 不考虑挥发分析出对煤粉粒径的影响． 1. 1 运动方程 根据牛顿定律，煤粉颗粒二维运动方程如下: mp dupx dt = Fdx， mp dupy dt = Fdy + Fg， up = u2 px + u2 槡 py        ． ( 1) mp dupx dt = 1 8 πCd ρgd2 p | upx － ugx | ( ugx － upx ) ， mp dupy dt = 1 8 πCd ρgd2 p | upy － ugy | ( ugy － upy ) + 1 6 πgd3 p ( ρp － ρg ) ． 式中: Fdx和 Fdy分别为煤粉颗粒运动过程在水平方 向和竖直方向受到的空气阻力，N; Fg为重力和浮力 的合力，N; mp 为煤粉颗粒的质量，kg; dp 为颗粒直 径，m; ρp和 ρg分别为碳粒和气体的密度，kg·m － 3 ; ug 为空气速度，m·s － 1 ; up为煤粉颗粒速度，m·s － 1 ; Cd 为颗粒在空气中运动的阻力系数，Cd = 24 Ｒe + 6 1 + 槡Ｒe + 0. 4，Ｒe 为雷诺数，Ｒe = | up － ug | dp L ，L 为 气体动力黏度，m2 ·s － 1 ． 1. 2 粒径方程 假设煤粉燃烧速率为 Gc，则煤粉质量方程为 dmp dt = Gc，mp = 1 6 πd3 p ρp ． 不考虑挥发分析出对煤粉粒径的影响，则煤粉 在燃烧过程中粒径变化方程为 ddp dt = － 2Gc πρp d2 p ． ( 2) 1. 3 热解方程 煤粉的热解方程采用 dV dt = ( V∞ － V) k0 ( exp － E ＲT ) p ． ( 3) 本文中的热解方程主要是对已有的热解方程进 行改进． 以往的热解方程中动力学参数的确定大部 分是依据经验，本文给出煤粉热解过程中动力学参 数的表达式． 通过对以往实验数据［15］进行拟合得 到热解动力学参数 k0和 E 的表达式: k0 = exp［8. 132exp( 9. 113 × 10 － 5T) + 0. 02399exp( 0. 003663T) ］， T≤1600; { [ ( exp 22. 9exp － T － 2173 ) 409. 5 ] 2 + 15. 37exp － [ ( T － 1674 ) 331. 3 ] } 2 { ， T ＞ 1600． E = 104 ×［76560exp( － 0. 01381T) + 0. 7519exp( 0. 001966T) ］， T ＜ 1850; 104 × 30. 84exp － { [ ( T － 2049 ) 630. 8 ] } 2 { ， T≥1850． 热解模型中 V∞ 表示挥发分最大产量，它是煤粉 热解通用模型中的重要参数，与煤种、加热条件、煤 粒尺寸等因素有关． V∞ = 1. 2( Vdaf ) 0. 8 [ exp － 2 × 106 Ｒ( T∞ － T0 ) 2 ] ． 式中: Vdaf为煤的可燃基挥发分; T0 为煤粉的初温， · 722 ·

·228+ 北京科技大学学报 第36卷 K;T为环境温度，K 炉壁与碳粒之间的辐射换热Q为 1.4炭粒温度 QR=5 S (T -T)=Ts do(T-T) 煤粉的热量交换主要包括三部分：(1)与热风 之间进行对流换热；(2)碳粒燃烧反应放热以及挥 式中，S。为煤粒的表面积，T.为壁温，e。为黑度系 发分析出吸热：(3)与炉壁之间辐射换热.对流传热 数，σ为玻尔兹曼常数 方程如下： 煤粉总的传热方程为， Qc=TdA (2+0.6Re.5Pro.33)(T.-T) :-Qc+Qu-Qw mpcpe dt (4) 式中，T.和T。分别为气体表面温度和颗粒的温度， 入。为气体导热系数 式中：cpe为煤粉颗粒比热容，J小(kgK)-1 热量被碳粒和气体吸收，反应热换Q为 1.5燃烧速率 =治74尝 影响煤粉燃烧状态的因素很多，例如煤粉粒径、 QM=dt 温度、环境氧气体积分数、氧的扩散系数等都将影响 式中：m。为碳粒质量，kg;Hc为碳粒的反应热，包 煤粉的燃烧状态。为了求解煤粉燃烧表面动力学参 括煤粉的热解和气化反应吸热和碳的氧化反应放 数，根据傅维镳的提出的煤燃烧通用规律，可得煤 热，Hnc=9797kJkg1. 粉燃烧反应动力学参数的通用表达式为： rko.h=4.109(Fz+27)8.%×10-20-(0.8363+0.7082b+0.2150F2+0.0267Fb3+0.00107Fb)ep(-Fb)], Fh>-2: ka-4109(+2)×i0a[a1637+05en(0g)】 Fh≤-2. (5) 式中：k。为反应频率因子，ms:Fz为煤质指数， 面气体扩散系数，m2·s-1;E为活化能，180kJ· Fz=(Va+Ma)2Ca×100:Fb为燃烧状态参数)， mol-;Nu为努赛尔数：R为气体常数.采用迭代法 当Fh-2; G.=2214×109'-1.057×10-6.81x107-2.2×10-+4.06×10-F2+1.953× 10-4Fb-9.551×10-4F3-3.179×10-3Fh2+2.254×10-2Fb+0.07724,D>10,Fb>-2; 2m是m-忌 Fh≤-2. D.Nu (7) 式中：D为空间反应邓克尔数，D=4.95×105× Y=Y。-(M,×Vom/)×(1-G), (8) rp/Nu:颗粒燃烧速率G.=2πd,p.D.G。,kgsl; G=(v+)/(1-A). (9) 1.6氧气体积分数以及燃烧率的计算 式中：Y。为初始氧气体积分数；G为煤粉燃烧率； 鼓风中的氧气主要用于煤粉和焦炭的燃烧，煤 Vo为单位质量煤粉耗氧量，0.933m3·kg1:M,为 粉在燃烧过程中周围的氧气体积分数逐渐变小，假 吨铁喷吹煤粉量，kg;V为吨铁鼓风量，计算过程中 设煤粉比焦炭优先燃烧，在煤粉未燃烬前氧气体积 为1250m3;V为释放的挥发分的质量分数；w.为己 分数的变化规律为 燃烧的碳质量分数；A,为煤粉灰分的质量分数

北 京 科 技 大 学 学 报 第 36 卷 K; T∞ 为环境温度，K． 1. 4 炭粒温度 煤粉的热量交换主要包括三部分: ( 1) 与热风 之间进行对流换热; ( 2) 碳粒燃烧反应放热以及挥 发分析出吸热; ( 3) 与炉壁之间辐射换热． 对流传热 方程如下: QC = πdpλg ( 2 + 0. 6Ｒe0. 5Pr0. 33 ) ( Tg － Tp ) ． 式中，Tg 和 Tp 分别为气体表面温度和颗粒的温度， λg 为气体导热系数． 热量被碳粒和气体吸收，反应热换 QM 为 QM = dmp dt Hreac = 1 2 ρp d2 p ddp dt Hreac ． 式中: mp 为碳粒质量，kg; Hreac 为碳粒的反应热，包 括煤粉的热解和气化反应吸热和碳的氧化反应放 热，Hreac = 9797 kJ·kg － 1 ． 炉壁与碳粒之间的辐射换热 QＲ 为 QＲ = εp Spσ( T4 p － T4 w ) = πεp d2 pσ( T4 p － T4 w ) ． 式中，Sp 为煤粒的表面积，Tw 为壁温，εp 为黑度系 数，σ 为玻尔兹曼常数． 煤粉总的传热方程为， mp cp，c dTp dt = QC + QM － QＲ． ( 4) 式中: cp，c为煤粉颗粒比热容，J·( kg·K) － 1 ． 1. 5 燃烧速率 影响煤粉燃烧状态的因素很多，例如煤粉粒径、 温度、环境氧气体积分数、氧的扩散系数等都将影响 煤粉的燃烧状态． 为了求解煤粉燃烧表面动力学参 数，根据傅维镳［15］提出的煤燃烧通用规律，可得煤 粉燃烧反应动力学参数的通用表达式为: k0，ch = 4. 109( Fz + 27) 18. 98 × 10 － 22［1 － ( 0. 8363 + 0. 7082Fb + 0. 2150Fb2 + 0. 0267Fb3 + 0. 00107Fb4 ) exp( － Fb) ］， Fb ＞ － 2; k0，ch = 4. 109( Fz + 27) 18. 98 × 10 [ － 22 ( 0. 1637 + 0. 86exp 0. 51 ) ] Fb + 1. 59 { ， Fb≤ － 2． ( 5) 式中: k0，ch为反应频率因子，m·s － 1 ; Fz 为煤质指数， Fz = ( Vad + Mad ) 2 Cad × 100; Fb 为燃烧状态参数［17］， 当 Fb ＜ － 2. 3 时燃烧反应属于动力控制，当 － 2. 3 ＜ Fb ＜ 8 时属于内孔扩散--动力控制，当 Fb≥8 时则属 于扩散控制． 定义如下: Fb = ln [ k0，ch dp DsNu Yo，∞ ( exp － E ＲT ) ] p ． ( 6) 式中: Yo，∞ 为质量分数，Yo，∞ = 2. 75Yo，∞ ; Ds为颗粒表 面气体 扩 散 系 数，m2 ·s － 1 ; E 为 活 化 能，180 kJ· mol － 1 ; Nu 为努赛尔数; Ｒ 为气体常数． 采用迭代法 计算 k0，ch和 Fb，直至| ( k0，ch ) n － ( k0，ch ) n － 1 |在允许的 误差范围内，初值( k0，ch ) 0 = 4. 109 × ( Fz + 27) 18. 98 × 10 － 22 ． 由于燃烧速率与 Fb 和 Dag有关，通过拟合实验 数据得到量纲一的通用燃烧速率公式: Gc = 1. 414 × 10 － 9Fb9 － 1. 483 × 10 － 8Fb8 － 6. 382 × 10 － 7Fb7 － 1. 167 × 10 － 6Fb6 + 4. 427 × 10 － 5Fb5 + 1. 26 × 10 － 4Fb4 － 1. 238 × 10 － 3Fb3 － 2. 441 × 10 － 3Fb2 + 0. 2748Fb + 0. 0829， Dag≤104 ，Fb ＞ － 2; 2. 214 × 10 － 9Fb9 － 1. 057 × 10 － 9Fb8 － 6. 811 × 10 － 7Fb7 － 2. 92 × 10 － 6Fb6 + 4. 06 × 10 － 5Fb5 + 1. 953 × 10 － 4Fb4 － 9. 551 × 10 － 4Fb3 － 3. 179 × 10 － 3Fb2 + 2. 254 × 10 － 2Fb + 0. 07724， Dag ＞ 104 ，Fb ＞ － 2; ( 12 /44) k0，ch dp DsNu Yo，∞ ( exp － E ＲT ) p ， Fb≤ － 2          ． ( 7) 式中: Dag为空间反应邓克尔数，Dag = 4. 95 × 1013 × r 2 p ρ 1. 5 s /Nu; 颗粒燃烧速率 Gc = 2πdp ρsDsGc，kg·s － 1 ; 1. 6 氧气体积分数以及燃烧率的计算 鼓风中的氧气主要用于煤粉和焦炭的燃烧，煤 粉在燃烧过程中周围的氧气体积分数逐渐变小，假 设煤粉比焦炭优先燃烧，在煤粉未燃烬前氧气体积 分数的变化规律为 Y = Y0 － ( Mp × V( O2) /Vb ) × ( 1 － G) ， ( 8) G = ( Vad + wc ) /( 1 － Aad ) ． ( 9) 式中: Y0 为初始氧气体积分数; G 为煤粉燃烧率; V( O2) 为单位质量煤粉耗氧量，0. 933 m3 ·kg － 1 ; Mp为 吨铁喷吹煤粉量，kg; Vb为吨铁鼓风量，计算过程中 为 1250 m3 ; Vad为释放的挥发分的质量分数; wc为已 燃烧的碳质量分数; Aad为煤粉灰分的质量分数． · 822 ·

第2期 陈川等：高炉煤粉燃烧率通用模型 ·229· 2 模型计算及验证 y(3)+K3)(4)+K4, 根据式(1)~(9)可得煤粉燃烧通用模型， y6)+2K,(5): 如下： =4,小 元=7（①+）y2)+K2), =,dT) y(3)+2K3)(4)+2K4, dt dTp y65)+2,5)）小： dt =f(V,d,Tp,u）, =f(,d,T,4,）, 元=7)+1y2)+K2), dt =6(W,dyTxur“n） y(3)+人3)④+2K④， dt 令子=于(d。,T。u)=（听5jff), (5)+2K(5): 应用四阶龙格库塔法对模型进行离散，离散格式 y=寸+(K1+2K2+2K3+K)h/6,t=ih. 如下： (ii)）output(x,y）stop. ()取时间步长h=10-5,t=0时，初始化各个 以上算法通过C语言编程实现，可以计算颗粒 参数， 在不同时刻的V,d,T,“r,“pm值，当i=n或者燃烧 y=(V(0),d.(0),T.(0),(0),u(0). 率G>0.99时计算结束. (i)Fori=l,2,…,n; 根据煤粉燃烧通用模型计算不同煤粉颗粒的燃 =j(y(1),y(2),y(3),y(4),y(5): 烧速率，将计算值与Smith和Tyler-9、Zhang 等0、Fu和Zhang0给出的实验数据进行比较，如 =)+K10y2)+2K2, 图1所示. ■预测值(d=89μm) ■预测值(d=1004um)） ■预测值(d=150μm) ●实验值Smith等，褐煤 ●实验值F如等.烟煤) ●实验值(Zhang等，无烟煤) ■ 00 10 ● 5 1000 1500 20001400 1600 1800 20001400 16001800 2000 温度K 温度K 温度K 图1煤粉燃烧速率预报值与实验值比较 Fig.I Comparison between the predicted results and previous experiment results 图1分别给出了褐煤、烟煤和无烟煤在不同 3结果分析 温度下的燃烧速率预报值和实验值.煤粉颗粒粒 径分别为89、100和150um,煤粉温度从1111~ 由于热风速度很高，因此煤粉在回旋区内的停 1873K,计算值和实验值吻合较好.因此，只需要 留时间很短.梁建国等回对高炉下部气固湍流进 知道煤粉的工业分析值，就能统一的确定煤粉的 行了研究，结果表明在风口回旋区内煤粉颗粒跟随 动力学参数，从而比较准确的预测煤粉的燃烧 热风作循环运动，煤粉颗粒在风口回旋区内的停留 速率. 时间在10~25ms,粒径较小的颗粒在回旋区做循环

第 2 期 陈 川等: 高炉煤粉燃烧率通用模型 2 模型计算及验证 根据式 ( 1 ) ～ ( 9 ) 可 得 煤 粉 燃 烧 通 用 模 型， 如下: dV dt = f1 ( V，dp，Tp ) ， ddp dt = f2 ( V，dp，Tp ) ， dTp dt = f3 ( V，dp，Tp，up ) ， dupx dt = f4 ( V，dp，Tp，up，upy ) ， dupy dt = f5 ( V，dp，Tp，up，upy )              ． 令 f  = f  ( V，dp，Tp，up，upy ) = ( f1，f2，f3，f4，f5 ) ， 应用四阶龙格库塔法对模型进行离散，离散格式 如下: ( i) 取时间步长 h = 10 － 5，t = 0 时，初始化各个 参数， y  = ( V( 0) ，dp ( 0) ，Tp ( 0) ，upx ( 0) ，upy ( 0) ) ． ( ii) For i = 1，2，…，n; K  1 = f  ( y( 1) ，y( 2) ，y( 3) ，y( 4) ，y( 5) ) ; K  2 = f (  y( 1) + h 2 K1 ( 1) ，y( 2) + h 2 K1 ( 2) ， y( 3) + h 2 K1 ( 3) ，y( 4) + h 2 K1 ( 4) ， y( 5) + h 2 K1 ( 5 ) ) ; K  3 = f (  y( 1) + h 2 K2 ( 1) ，y( 2) + h 2 K2 ( 2) ， y( 3) + h 2 K2 ( 3) ，y( 4) + h 2 K2 ( 4) ， y( 5) + h 2 K2 ( 5 ) ) ; K  4 = f (  y( 1) + h 2 K4 ( 1) ，y( 2) + h 2 K4 ( 2) ， y( 3) + h 2 K4 ( 3) ，y( 4) + h 2 K4 ( 4) ， y( 5) + h 2 K4 ( 5 ) ) ; y  = y  + ( K  1 + 2K  2 + 2K  3 + K  4 ) h /6，t = i·h． ( iii) output( x，y) stop． 以上算法通过 C 语言编程实现，可以计算颗粒 在不同时刻的 V，dp，Tp，upx，upy值，当 i = n 或者燃烧 率 G ＞ 0. 99 时计算结束． 根据煤粉燃烧通用模型计算不同煤粉颗粒的燃 烧 速 率，将 计 算 值 与 Smith 和 Tyler［18--19］、Zhang 等［20］、Fu 和 Zhang［21］给出的实验数据进行比较，如 图 1 所示． 图 1 煤粉燃烧速率预报值与实验值比较 Fig． 1 Comparison between the predicted results and previous experiment results 图 1 分别给出了褐煤、烟煤和无烟煤在不同 温度下的燃烧速率预报值和实验值． 煤粉颗粒粒 径分别为 89、100 和 150 μm，煤粉温度从 1111 ～ 1873 K，计算值和实验值吻合较好． 因此，只需要 知道煤粉的工业分析值，就能统一的确定煤粉的 动力学 参 数，从而比较准确的预测煤粉的燃烧 速率． 3 结果分析 由于热风速度很高，因此煤粉在回旋区内的停 留时间很短． 梁建国等［22］对高炉下部气固湍流进 行了研究，结果表明在风口回旋区内煤粉颗粒跟随 热风作循环运动，煤粉颗粒在风口回旋区内的停留 时间在 10 ～ 25 ms，粒径较小的颗粒在回旋区做循环 · 922 ·

·230 北京科技大学学报 第36卷 运动的时间可达到1、。为了研究煤粉粒径、煤粉成 250 分、风温、氧气体积分数以及喷吹量对煤粉燃烧率的 200 影响，计算过程中假设煤粉粒径分布均匀，不考虑煤 粉粒径对停留时间的影响，煤粉在回旋区内的停留 150 -d.=64um 时间均为25ms.计算过程中用到的参数如表1 d=74μm 7-d=94um 所示. 50 表1煤粉燃烧率计算参数 4 6 1 Table 1 Parameters in calculating the combustion ratio 时间ms 参数 数值 图2煤粉粒径对速度的影响 挥发分析出量，()。 0 Fig.2 Influence of coal diameter on the particle velocity 粒径，(d.)o/μm 多 799和694K,进入回旋区后煤粉颗粒达到理论燃烧 纵向速度，(4w)。/(ms) 0 温度后保持不变.从图中可以看出煤粉颗粒的升温 鼓风温度，TK 1473 速率非常大，当煤粉粒径为94um时，在8ms内，颗 初始氧气体积分数，Y 0.23 粒温度由300K上升至2373K,平均升温速率高达 挥发分的质量分数，V 0.1445 2.6×10K·s1,并且煤粉颗粒越小，升温速率 固定碳的质量分数，C 0.6220 越高 黑度系数，6p 0.8 气体常数，R/(J.mol-1) 8.3145 2500 反应热，H/(kJkg1) 9797 2000 吨铁煤粉喷吹量，Mkg 190 1500 喷枪位置，PLlm 0.5 量一d=64um 1000 d=74 um 煤粉初温，(T,)。K 300 d=94μm 500 横向速度，(umo/(ms) 5 4 6 10 鼓风风速，山，/(ms) 230 时间/ms 理论燃烧温度，TK 2373 图3不同煤粉颗粒温度变化 煤粉密度，p。/(kgm3) 1600 Fig.3 Change in temperature with time for different particles 水分的质量分数，Ma 0.0096 图4为煤粉颗粒粒径对挥发分析出量（质量分 灰分的质量分数，Ad 0.2239 数)的影响.从图中可以看出煤粉挥发分几乎瞬间 普朗特数，PA 0.7 析出，三种粒度的煤粉颗粒最大挥发分析出量相同， 非氧化反应活化能，E/(k·mol-l) 180 但是挥发分析出的时间不同，粒径为64、74和94μm 碳的热值，q。/kg) 3.34×107 的煤粉颗粒挥发分开始析出时间分别为2.7、3.2和 煤粉停留时间，r/ms 25 4.2ms,煤粉颗粒越小，挥发分开始析出时间越早. 注：喷枪位置PL表示喷枪与风口前端之间的距离 计算过程中喷枪出口位置距离风口前端为0.5m, 3.1煤粉粒径对燃烧率的影响 当煤粉颗粒粒径为64、74和94m时，煤粉颗粒在 煤粉颗粒的速度和温度直接影响颗粒的燃烧 直吹管的停留时间分别为2.31、2.35和2.43ms,如 率，因此首先研究了煤粉粒径对煤粉速度和温度的 图5所示.颗粒粒径越小，在直吹管的停留时间越 影响，如图2和图3所示. 短.由于煤粉颗粒在直吹管内的停留时间小于挥发 由图2可以看出，煤粉颗粒在很短的时间内就 分开始析出时间，因此煤粉颗粒在直吹管内几乎不 达到了热风速度.由于小颗粒的煤粉跟随性好，所 析出挥发分. 以煤粉粒径越小，达到热风速度所需的时间越短 图6为燃烧状态参数及反应频率因子随时间变 在2ms之后，不同粒径的煤粉颗粒速度均为热风速 化规律.从图中可以看出：在燃烧前期，Fb<-2.3, 度.图3给出了不同粒径的煤粉颗粒温度随时间的 煤粉燃烧属于动力控制：燃烧反应进行10s之后， 变化曲线.在直吹管内，环境温度为1473K,粒径为 -2.3<Fb<8,燃烧反应属于内孔扩散一动力控制， 64、74和94μm的煤粉颗粒最高温度分别为876、 Fb和k。.变化幅度非常小，近似为常数；当煤粉粒

北 京 科 技 大 学 学 报 第 36 卷 运动的时间可达到 1 s． 为了研究煤粉粒径、煤粉成 分、风温、氧气体积分数以及喷吹量对煤粉燃烧率的 影响，计算过程中假设煤粉粒径分布均匀，不考虑煤 粉粒径对停留时间的影响，煤粉在回旋区内的停留 时间均为 25 ms． 计算过程中用到的参数如表 1 所示． 表 1 煤粉燃烧率计算参数 Table 1 Parameters in calculating the combustion ratio 参数 数值 挥发分析出量，( V) 0 0 粒径，( dp ) 0 /μm 74 纵向速度，( upy ) 0 /( m·s － 1 ) 0 鼓风温度，Tf /K 1473 初始氧气体积分数，Y 0. 23 挥发分的质量分数，Vad 0. 1445 固定碳的质量分数，Cad 0. 6220 黑度系数，εp 0. 8 气体常数，Ｒ /( J·mol － 1 ) 8. 3145 反应热，Hreac /( kJ·kg － 1 ) 9797 吨铁煤粉喷吹量，Mp /kg 190 喷枪位置，PL /m 0. 5 煤粉初温，( Tp ) 0 /K 300 横向速度，( upx ) 0 /( m·s － 1 ) 5 鼓风风速，ug /( m·s － 1 ) 230 理论燃烧温度，T /K 2373 煤粉密度，ρp /( kg·m － 3 ) 1600 水分的质量分数，Mad 0. 0096 灰分的质量分数，Aad 0. 2239 普朗特数，Pr 0. 7 非氧化反应活化能，E /( kJ·mol － 1 ) 180 碳的热值，qc /( J·kg － 1 ) 3. 34 × 107 煤粉停留时间，τ / ms 25 注: 喷枪位置 PL 表示喷枪与风口前端之间的距离． 3. 1 煤粉粒径对燃烧率的影响 煤粉颗粒的速度和温度直接影响颗粒的燃烧 率，因此首先研究了煤粉粒径对煤粉速度和温度的 影响，如图 2 和图 3 所示． 由图 2 可以看出，煤粉颗粒在很短的时间内就 达到了热风速度． 由于小颗粒的煤粉跟随性好，所 以煤粉粒径越小，达到热风速度所需的时间越短． 在 2 ms 之后，不同粒径的煤粉颗粒速度均为热风速 度． 图 3 给出了不同粒径的煤粉颗粒温度随时间的 变化曲线． 在直吹管内，环境温度为 1473 K，粒径为 64、74 和 94 μm 的煤粉颗粒最高温度分别为 876、 图 2 煤粉粒径对速度的影响 Fig． 2 Influence of coal diameter on the particle velocity 799 和 694 K，进入回旋区后煤粉颗粒达到理论燃烧 温度后保持不变． 从图中可以看出煤粉颗粒的升温 速率非常大，当煤粉粒径为 94 μm 时，在 8 ms 内，颗 粒温度由 300 K 上升至 2373 K，平均升温速率高达 2. 6 × 105 K·s － 1，并且煤粉颗粒越小，升温 速 率 越高． 图 3 不同煤粉颗粒温度变化 Fig． 3 Change in temperature with time for different particles 图 4 为煤粉颗粒粒径对挥发分析出量( 质量分 数) 的影响． 从图中可以看出煤粉挥发分几乎瞬间 析出，三种粒度的煤粉颗粒最大挥发分析出量相同， 但是挥发分析出的时间不同，粒径为 64、74 和94 μm 的煤粉颗粒挥发分开始析出时间分别为 2. 7、3. 2 和 4. 2 ms，煤粉颗粒越小，挥发分开始析出时间越早． 计算过程中喷枪出口位置距离风口前端为 0. 5 m， 当煤粉颗粒粒径为 64、74 和 94 μm 时，煤粉颗粒在 直吹管的停留时间分别为 2. 31、2. 35 和 2. 43 ms，如 图 5 所示． 颗粒粒径越小，在直吹管的停留时间越 短． 由于煤粉颗粒在直吹管内的停留时间小于挥发 分开始析出时间，因此煤粉颗粒在直吹管内几乎不 析出挥发分． 图 6 为燃烧状态参数及反应频率因子随时间变 化规律． 从图中可以看出: 在燃烧前期，Fb ＜ － 2. 3， 煤粉燃烧属于动力控制; 燃烧反应进行 10 ms 之后， － 2. 3 ＜ Fb ＜ 8，燃烧反应属于内孔扩散--动力控制， Fb 和 k0，ch变化幅度非常小，近似为常数; 当煤粉粒 · 032 ·

第2期 陈川等：高炉煤粉燃烧率通用模型 ·231· 粉越细安全问题越大，所以必须选择合适的煤粉 15 粒度. 10 10 25 5 量一d=64m 0 20 d =74 um ---74m =94μm -10 -94m 0 20 h 4 6 10 时间/m“ 5 30 图4不同煤粉颗粒挥发分析出量 1015 20 25 0 Fig.4 Released volatile amount for different particles 时间/ms 图6燃烧状态参数及反应频率因子随时间变化规律 3.0 Fig.6 Variation of combustion state parameter and kinetic parameter with time 100 量-d=64山m 80 -d=74m d=94 pm 60 40 64μm 74 pm 94μm 煤粉粒径 20 图5煤粉颗粒在直吹管内的停留时间 Fig.5 Residence time of pulverized coal particles in the blowpipe 0 5 10152025 时间/ms 径分别为74um和94um时，燃烧状态参数Fb分别 图7粒径对煤粉颗粒燃烧率的影响 为2.4和2.7，反应频率因子分别为1.1×10°和 Fig.7 Influence of particle diameter on the combustion ratio 1.2×10°m·s-1,因此煤粉粒径对燃烧反应动力学 参数影响较小.图7给出了不同粒度煤粉颗粒的燃 3.2煤粉成分对燃烧率的影响 烧率变化曲线.从图中可以看出，在2.5ms以内，煤 图8为挥发分含量对燃烧状态参数的影响.从 粉燃烧率近似为0，由于煤粉在直吹管内停留时间 图中可以看出：在燃烧前期，Fb<-2.3,煤粉燃烧属 小于2.5ms,因此煤粉颗粒在直吹管内不会燃烧. 于动力控制；在燃烧后期，-2.3<Fb<8,燃烧反应 Bortz和Flament在风温为1100~1200℃，氧体积 属于内孔扩散一动力控制.挥发分含量越高，反应状 分数为1%~2%，风速为200m·s-的条件下研究 态参数和反应频率因子越大，说明煤粉成分对燃烧 煤粉燃烧，发现无烟煤在3s内煤粉几乎没有燃 动力学参数影响很大，在计算煤粉燃烧率时必须考 烧.模型计算结果与前人观察到的现象相吻合.当 虑煤种对反应速率系数的影响.图9分别计算了挥 煤粉颗粒从直吹管进入风口回旋区后，环境温度瞬 发分质量分数为4.5%、14.5%和24.5%的煤粉颗 间达到2100℃左右，挥发分迅速析出，因此燃烧率 粒的燃烧率.由于煤粉在直吹管内燃烧速率很低， 突然增大.图7中曲线的斜率表示煤粉燃烧速率与 所以煤粉在直吹管内的燃烧率近似为零；当煤粉颗 煤粉颗粒燃烧前固定碳以及挥发分质量之和的比 粒进入风口回旋区后，煤粉燃烧率迅速提高。在相 值，由于10ms之后煤粉燃烧状态参数Fb变化很小 同时刻，煤粉挥发分含量越高，煤粉燃烧率越高.当 (见图6)，根据式(7)可得，煤粉在燃烧过程中燃烧 燃烧反应进行5ms之后，由于燃烧状态参数Fb近 速率变化不大.因此，煤粉燃烧率与时间近似呈线 似为常数，所以煤粉燃烧率随时间呈线性关系。并 性关系.计算过程中三种粒径的煤粉成分相同，由 且通过对比图中虚线可以看出，挥发分含量高的煤 于大颗粒煤粉挥发分以及固定碳的质量大，所以煤 粉燃烧率曲线斜率较大，说明煤粉挥发分含量越高， 粉颗粒越大，煤粉燃烧率曲线的斜率越小.在25s 煤粉燃烧速率越大.当燃烧反应进行25s时，上述 时上述三种粒径煤粉的燃烧率分别93.5%、72.5% 三种煤粉颗粒的燃烧率分别为53.6%、72.5%和 和51.2%，粒径平均每减小10m,燃烧率提高 91.2%.煤粉挥发分质量分数每提高10%，煤粉燃 10%~20%.但是考虑到磨煤的动力消耗，以及煤 烧率提高18%~19%.因此，当煤粉在回旋内不能

第 2 期 陈 川等: 高炉煤粉燃烧率通用模型 图 4 不同煤粉颗粒挥发分析出量 Fig． 4 Ｒeleased volatile amount for different particles 图 5 煤粉颗粒在直吹管内的停留时间 Fig． 5 Ｒesidence time of pulverized coal particles in the blowpipe 径分别为 74 μm 和 94 μm 时，燃烧状态参数 Fb 分别 为 2. 4 和 2. 7，反应频率因子分别为 1. 1 × 106 和 1. 2 × 106 m·s － 1，因此煤粉粒径对燃烧反应动力学 参数影响较小． 图 7 给出了不同粒度煤粉颗粒的燃 烧率变化曲线． 从图中可以看出，在 2. 5 ms 以内，煤 粉燃烧率近似为 0，由于煤粉在直吹管内停留时间 小于 2. 5 ms，因此煤粉颗粒在直吹管内不会燃烧． Bortz 和 Flament［23］在风温为 1100 ～ 1200 ℃，氧体积 分数为 1% ～ 2% ，风速为 200 m·s － 1 的条件下研究 煤粉燃烧，发现无烟煤在 3 ms 内煤粉几乎没有燃 烧． 模型计算结果与前人观察到的现象相吻合． 当 煤粉颗粒从直吹管进入风口回旋区后，环境温度瞬 间达到 2100 ℃左右，挥发分迅速析出，因此燃烧率 突然增大． 图 7 中曲线的斜率表示煤粉燃烧速率与 煤粉颗粒燃烧前固定碳以及挥发分质量之和的比 值，由于 10 ms 之后煤粉燃烧状态参数 Fb 变化很小 ( 见图 6) ，根据式( 7) 可得，煤粉在燃烧过程中燃烧 速率变化不大． 因此，煤粉燃烧率与时间近似呈线 性关系． 计算过程中三种粒径的煤粉成分相同，由 于大颗粒煤粉挥发分以及固定碳的质量大，所以煤 粉颗粒越大，煤粉燃烧率曲线的斜率越小． 在 25 ms 时上述三种粒径煤粉的燃烧率分别 93. 5% 、72. 5% 和 51. 2% ，粒径平均每减小 10 μm，燃 烧 率 提 高 10% ～ 20% ． 但是考虑到磨煤的动力消耗，以及煤 粉越细安全问题越大，所以必须选择合适的煤粉 粒度． 图 6 燃烧状态参数及反应频率因子随时间变化规律 Fig． 6 Variation of combustion state parameter and kinetic parameter with time 图 7 粒径对煤粉颗粒燃烧率的影响 Fig． 7 Influence of particle diameter on the combustion ratio 3. 2 煤粉成分对燃烧率的影响 图 8 为挥发分含量对燃烧状态参数的影响． 从 图中可以看出: 在燃烧前期，Fb ＜ － 2. 3，煤粉燃烧属 于动力控制; 在燃烧后期，－ 2. 3 ＜ Fb ＜ 8，燃烧反应 属于内孔扩散--动力控制． 挥发分含量越高，反应状 态参数和反应频率因子越大，说明煤粉成分对燃烧 动力学参数影响很大，在计算煤粉燃烧率时必须考 虑煤种对反应速率系数的影响． 图 9 分别计算了挥 发分质量分数为 4. 5% 、14. 5% 和 24. 5% 的煤粉颗 粒的燃烧率． 由于煤粉在直吹管内燃烧速率很低， 所以煤粉在直吹管内的燃烧率近似为零; 当煤粉颗 粒进入风口回旋区后，煤粉燃烧率迅速提高． 在相 同时刻，煤粉挥发分含量越高，煤粉燃烧率越高． 当 燃烧反应进行 5 ms 之后，由于燃烧状态参数 Fb 近 似为常数，所以煤粉燃烧率随时间呈线性关系． 并 且通过对比图中虚线可以看出，挥发分含量高的煤 粉燃烧率曲线斜率较大，说明煤粉挥发分含量越高， 煤粉燃烧速率越大． 当燃烧反应进行 25 ms 时，上述 三种煤粉颗粒的燃烧率分别为 53. 6% 、72. 5% 和 91. 2% ． 煤粉挥发分质量分数每提高 10% ，煤粉燃 烧率提高 18% ～ 19% ． 因此，当煤粉在回旋内不能 · 132 ·

·232· 北京科技大学学报 第36卷 充分燃烧时，可以通过提高煤粉挥发分含量或采用 78 混合喷吹提高煤粉燃烧率 16 10 0 Fb 74 0 72 -20 ---V=45% =24.5% 20 0 00 1000 1100 1200 300 风温℃ 0 10152025 图10风温对燃烧率的影响 时间ms Fig.10 Influence of blast temperature on the combustion ratio 图8挥发分含量对燃烧状态参数及反应频率因子的影响 Fig.8 Influence of volatile content on the combustion state parameter and kinetic parameter 85 100 量一V=45% 80 -V-14.5% 80 24.59% 60 70 21 23 25 27 20 20H 氧气体积分数/% 图11氧气体积分数对燃烧率的影响 0 5 10152025 Fig.11 Influence of oxygen volume fraction on the combustion ratio 时间/ms 图9挥发分含量对煤粉燃烧率的影响 的影响比较显著.安阳1号高炉生产实践表明：吨 Fig.9 Influence of volatile content on the combustion ratio 铁煤粉喷煤量增加20kg,相应的氧体积分数应提高 3.3操作参数对燃烧率的影响 1%,如此可维持较高的煤粉燃烧率的.因此，在大 图10为风温对煤粉燃烧率的影响.从图中可 喷煤量的条件下，为了保证较高的燃烧率，实际操作 以看出，风温越高，煤粉燃烧率越高，当风温从 过程中应采用富氧喷吹 900℃提高至1300℃时，燃烧率提高了5%.但是， 95 煤粉燃烧率与风温之间并不是线性关系.当风温低 90 于1100℃时，燃烧率受风温的影响比较显著：当风 85 温超过1100℃时，风温对煤粉燃烧率的影响逐渐减 弱.图11为氧气体积分数对煤粉燃烧率的影响.燃 烧率随氧气体积分数增加近似呈线性增长.氧气体 积分数平均增加1%，燃烧率提高2%左右.许莹和 100 130 150170 190 胡宾生通过实验发现，当氧气体积分数从0提高 吨铁喷吹量从g 到3%时，煤粉的燃烧率提高6%~7%左右.计算 图12喷吹量对煤粉燃烧率的影响 结果与前人实验结果吻合较好.由于在风口回旋区 Fig.12 Influence of injection amount on the combustion ratio 内，固定碳的燃烧速率受氧在碳颗粒表面传质速度 的限制，因此鼓风中氧气体积分数的增加有助于提 4 结论 高固定碳的燃烧速率，改善煤粉在高炉风口区的燃 (1)在高炉喷煤过程中，煤粉颗粒在2ms左右 烧率 就可以达到热风速度，由于煤粉颗粒在直吹管内停 图12为喷吹量对煤粉燃烧率的影响.从图中 留时间短并且温度较低，因此在直吹管内煤粉不会 可以看出，随着喷吹量的增加，煤粉燃烧率逐渐降 发生燃烧煤粉进入风口回旋区后，挥发分瞬间全 低.当吨铁煤粉喷吹量由100kg提高到190kg时， 部析出，并且颗粒粒径越小，挥发分开始析出时间 煤粉燃烧率下降了19%，因此喷煤量对煤粉燃烧率 越早

北 京 科 技 大 学 学 报 第 36 卷 充分燃烧时，可以通过提高煤粉挥发分含量或采用 混合喷吹提高煤粉燃烧率． 图 8 挥发分含量对燃烧状态参数及反应频率因子的影响 Fig． 8 Influence of volatile content on the combustion state parameter and kinetic parameter 图 9 挥发分含量对煤粉燃烧率的影响 Fig． 9 Influence of volatile content on the combustion ratio 3. 3 操作参数对燃烧率的影响 图 10 为风温对煤粉燃烧率的影响． 从图中可 以看出，风 温 越 高，煤粉燃烧率越高，当 风 温 从 900 ℃提高至 1300 ℃ 时，燃烧率提高了 5% ． 但是， 煤粉燃烧率与风温之间并不是线性关系． 当风温低 于 1100 ℃时，燃烧率受风温的影响比较显著; 当风 温超过 1100 ℃时，风温对煤粉燃烧率的影响逐渐减 弱． 图11 为氧气体积分数对煤粉燃烧率的影响． 燃 烧率随氧气体积分数增加近似呈线性增长． 氧气体 积分数平均增加 1% ，燃烧率提高 2% 左右． 许莹和 胡宾生［24］通过实验发现，当氧气体积分数从 0 提高 到 3% 时，煤粉的燃烧率提高 6% ～ 7% 左右． 计算 结果与前人实验结果吻合较好． 由于在风口回旋区 内，固定碳的燃烧速率受氧在碳颗粒表面传质速度 的限制，因此鼓风中氧气体积分数的增加有助于提 高固定碳的燃烧速率，改善煤粉在高炉风口区的燃 烧率． 图 12 为喷吹量对煤粉燃烧率的影响． 从图中 可以看出，随着喷吹量的增加，煤粉燃烧率逐渐降 低． 当吨铁煤粉喷吹量由 100 kg 提高到 190 kg 时， 煤粉燃烧率下降了 19% ，因此喷煤量对煤粉燃烧率 图 10 风温对燃烧率的影响 Fig． 10 Influence of blast temperature on the combustion ratio 图 11 氧气体积分数对燃烧率的影响 Fig． 11 Influence of oxygen volume fraction on the combustion ratio 的影响比较显著． 安阳 1 号高炉生产实践表明: 吨 铁煤粉喷煤量增加 20 kg，相应的氧体积分数应提高 1% ，如此可维持较高的煤粉燃烧率［25］． 因此，在大 喷煤量的条件下，为了保证较高的燃烧率，实际操作 过程中应采用富氧喷吹． 图 12 喷吹量对煤粉燃烧率的影响 Fig． 12 Influence of injection amount on the combustion ratio 4 结论 ( 1) 在高炉喷煤过程中，煤粉颗粒在 2 ms 左右 就可以达到热风速度，由于煤粉颗粒在直吹管内停 留时间短并且温度较低，因此在直吹管内煤粉不会 发生燃烧． 煤粉进入风口回旋区后，挥发分瞬间全 部析出，并且颗粒粒径越小，挥发分开始析出时间 越早． · 232 ·

第2期 陈川等：高炉煤粉燃烧率通用模型 ·233· (2)煤粉成分对燃烧动力学参数影响很大，在 (2):163 计算煤粉燃烧率时必须考虑煤种对反应动力学参数 (张丽丽，周慎杰，陈举华.高炉三维气固湍流和煤粉燃烧过 程数值模拟.南京理工大学学报，2007,31(2)：163) 的影响.煤粉挥发分含量每提高10%，煤粉燃烧率 Zhang S F,Wen L Y,Bai C G,et al.Numerical model of tuyere 提高18%~19%.因此，当煤粉在回旋内不能充分 pulverized coal and raceway coke combustion in BF.Chongqing 燃烧时，可以通过提高煤粉挥发分含量或采用混合 Unig Nat Sci Ed,2005,28(8):42 喷吹提高煤粉燃烧率。 (张生富，温良英，白晨光，等.高炉风口煤粉及回旋区焦炭 (3)提高风温有利于提高煤粉燃烧率，但是随 燃烧过程数学模型.重庆大学学报：自然科学版，2005,28 着风温的提高，风温对煤粉燃烧率的影响逐渐减弱. (8):42) (4)降低煤粉粒径、增加氧气体积分数可以有 [10]YangT J,Cang D Q.Ding Y L Oxygen Enrichment Pulrerized Coal in Blast Furnace.Beijing:Metallurgical Industry Press, 效地提高煤粉燃烧率.煤粉颗粒粒径每减小10m, 1996 燃烧率提高10%~20%.氧气体积分数每增加 (杨天钧，苍大强，丁玉龙.高炉宫氧煤粉喷吹.北京：治金 1%,燃烧率提高2%.随着喷煤量的增加，煤粉燃烧 工业出版社，1996) 率逐渐降低.因此，当提高煤粉喷吹量时，为了保证 [11]Cen K F,Yao Q,Luo Z Y,et al.Combustion Theory and Pollu- 较高的燃烧率，实际操作过程中应提高富氧率并适 tion Control.Beijing:Mechanical Industry Press,2008 (岑可法，姚强，骆仲泱，等.燃烧理论与污染控制.北京： 当降低煤粉粒径 机械工业出版社，2008) [12]Shen Y S,Guo B Y,Yu A B,et al.Model study of the effects of 参考文献 coal properties and blast conditions on pulverized coal combus- tion.1SJt,2009,49(6):819 [Kong D W,Zhang JL,Lin X H,et al.Combustion characteris- [13]Shen Y S,Guo B Y,Yu A B,et al.Three-dimensional model- tics and kinetic analysis of blended coal in BF.J fron Steel Res, ling of coal combustion in blast furnace.IS/J Int,2008,48(6): 2011,23(11):4 777 (孔德文，张建良，林祥海，等.高炉喷吹混合煤的燃烧特性 [14]Guo B Y,Zulli P.Rogers H,et al.Three-dimensional simula- 及动力学分析.钢铁研究学报，2011,23(11)：4) tion of flow and combustion for pulverised coal injection.IS/ Raygan S,Abdizadeh H,Eskandari R A.Evaluation of four coals lnt,2005,45(9):1272 for blast fumace pulverized coal injection.J Iron Steel Res Int, 15]Fu W B.Coal Combustion Theory and Macroscopic General Rule. 2010,17(3):8 Beijing:Tsinghua University Press,2003 B]Shen Y S,Guo B Y,Yu A B,et al.A three-dimensional numeri- (傅维镳.煤燃烧理论及其宏观通用规律。北京：清华大学 cal study of the combustion of coal blends in blast furnace.Fuel, 出版社，2003) 2009,88(2):255 [16]Han X L.Numerical Analysis.Beijing:Higher Education Press, 4]Jiang H B.Zhang JL,Fu JX,et al.Properties and structural op- 2011 timization of pulverized coal for blast fumace injection.J Iron Steel (韩旭里.数值分析.北京：高等教育出版社，2011) Resm,2011,18(3):6 [17]Fu W B,Zhang B L.A simple and general method for calculat- Tang G L,Na S R.The calculation of pulverized coal in raceway ing the buming rate of coal char/carbon particles in air.Combust cavity and the behaviour of the unbumnt pulverized coal.Baotou Sci Technol,1993,89:405 Univ Iron Steel Technal,1996,15(3):38 [18]Smith I W,Tyler R J.The reactivity of a porous brown coal char (唐光临，那树人.煤粉燃烧率的计算及未燃煤粉的行为.包 to oxygen between 630 and 1812 K.Combust Sci Technol,1974 头钢铁学院学报，1996,15(3)：38) 9(314):87 [6]Xu W R,Wu K,Zhang LL,et al.Research on utilization factor [19]Smith I W,Tyler R J.Internal burning of pulverized semi-an- of injected into BF pulverized coal at large PCI rate.fron Steel, thracite:the relation between particle structure and reactivity. 2006,41(4):10 Fuel,1972,51(4):312 (徐万仁，吴铿，张龙来，等.高炉大喷煤时煤粉利用率的研 [20] Zhang X W,Lu PS,Cao P F,et al.The kinetics of combustion 究.钢铁，2006,41(4)：10) of pulverized anthracite char in the temperature range 1200 to Zhang H B,Qin G T,Ji R S,et al.Numerical simulation of pul- 2000 K.J Eng Thermophys,1991,12(1):78 verized-coal combustion.J Beijing Univ Aeronaut Astronaut, (张喜文，路不思，曹鹏飞，等.粉状无烟煤炭粒在1200- 2009,35(5):536 2000K温度范围内燃烧的化学动力学特性.工程热物理学 (张宏博，秦国彤，纪任山，等.煤粉燃烧过程的数值模拟 报，1991,12(1)：78) 北京航空航天大学学报，2009,35(5)：536) 1]Fu W B.Zhang B L.A generalized relationship between the ki- [8]Zhang LL,Zhou S J,Chen J H.Numerical simulations of three- netic parameters of char combustion and coalproperties.Combust dimensional gas-particle turbulent flow and pulverized coal com- Sci Technol,1997,3(1)1 bustion in blast furnace.J Nanjing Unin Sci Technol,2007,31 (傅维标，张百立.煤焦燃烧反应动力学参数与煤种的通用

第 2 期 陈 川等: 高炉煤粉燃烧率通用模型 ( 2) 煤粉成分对燃烧动力学参数影响很大，在 计算煤粉燃烧率时必须考虑煤种对反应动力学参数 的影响． 煤粉挥发分含量每提高 10% ，煤粉燃烧率 提高 18% ～ 19% ． 因此，当煤粉在回旋内不能充分 燃烧时，可以通过提高煤粉挥发分含量或采用混合 喷吹提高煤粉燃烧率． ( 3) 提高风温有利于提高煤粉燃烧率，但是随 着风温的提高，风温对煤粉燃烧率的影响逐渐减弱． ( 4) 降低煤粉粒径、增加氧气体积分数可以有 效地提高煤粉燃烧率． 煤粉颗粒粒径每减小 10 μm， 燃烧 率 提 高 10% ～ 20% ． 氧 气 体 积 分 数 每 增 加 1% ，燃烧率提高 2% ． 随着喷煤量的增加，煤粉燃烧 率逐渐降低． 因此，当提高煤粉喷吹量时，为了保证 较高的燃烧率，实际操作过程中应提高富氧率并适 当降低煤粉粒径． 参 考 文 献 ［1］ Kong D W，Zhang J L，Lin X H，et al． Combustion characteris￾tics and kinetic analysis of blended coal in BF． J Iron Steel Ｒes， 2011，23( 11) : 4 ( 孔德文，张建良，林祥海，等． 高炉喷吹混合煤的燃烧特性 及动力学分析． 钢铁研究学报，2011，23( 11) : 4) ［2］ Ｒaygan S，Abdizadeh H，Eskandari Ｒ A． Evaluation of four coals for blast furnace pulverized coal injection． J Iron Steel Ｒes Int， 2010，17( 3) : 8 ［3］ Shen Y S，Guo B Y，Yu A B，et al． A three-dimensional numeri￾cal study of the combustion of coal blends in blast furnace． Fuel， 2009，88( 2) : 255 ［4］ Jiang H B，Zhang J L，Fu J X，et al． Properties and structural op￾timization of pulverized coal for blast furnace injection． J Iron Steel Ｒes Int，2011，18( 3) : 6 ［5］ Tang G L，Na S Ｒ． The calculation of pulverized coal in raceway cavity and the behaviour of the unburnt pulverized coal． J Baotou Univ Iron Steel Technol，1996，15( 3) : 38 ( 唐光临，那树人． 煤粉燃烧率的计算及未燃煤粉的行为． 包 头钢铁学院学报，1996，15( 3) : 38) ［6］ Xu W Ｒ，Wu K，Zhang L L，et al． Ｒesearch on utilization factor of injected into BF pulverized coal at large PCI rate． Iron Steel， 2006，41( 4) : 10 ( 徐万仁，吴铿，张龙来，等． 高炉大喷煤时煤粉利用率的研 究． 钢铁，2006，41( 4) : 10) ［7］ Zhang H B，Qin G T，Ji Ｒ S，et al． Numerical simulation of pul￾verized-coal combustion． J Beijing Univ Aeronaut Astronaut， 2009，35( 5) : 536 ( 张宏博，秦国彤，纪任山，等． 煤粉燃烧过程的数值模拟． 北京航空航天大学学报，2009，35( 5) : 536) ［8］ Zhang L L，Zhou S J，Chen J H． Numerical simulations of three￾dimensional gas-particle turbulent flow and pulverized coal com￾bustion in blast furnace． J Nanjing Univ Sci Technol，2007，31 ( 2) : 163 ( 张丽丽，周慎杰，陈举华． 高炉三维气固湍流和煤粉燃烧过 程数值模拟． 南京理工大学学报，2007，31( 2) : 163) ［9］ Zhang S F，Wen L Y，Bai C G，et al． Numerical model of tuyere pulverized coal and raceway coke combustion in BF． J Chongqing Univ Nat Sci Ed，2005，28( 8) : 42 ( 张生富，温良英，白晨光，等． 高炉风口煤粉及回旋区焦炭 燃烧过程数学模型． 重庆大学学报: 自然科学版，2005，28 ( 8) : 42) ［10］ Yang T J，Cang D Q，Ding Y L． Oxygen Enrichment Pulverized Coal in Blast Furnace． Beijing: Metallurgical Industry Press， 1996 ( 杨天钧，苍大强，丁玉龙． 高炉富氧煤粉喷吹． 北京: 冶金 工业出版社，1996) ［11］ Cen K F，Yao Q，Luo Z Y，et al． Combustion Theory and Pollu￾tion Control． Beijing: Mechanical Industry Press，2008 ( 岑可法，姚强，骆仲泱，等． 燃烧理论与污染控制． 北京: 机械工业出版社，2008) ［12］ Shen Y S，Guo B Y，Yu A B，et al． Model study of the effects of coal properties and blast conditions on pulverized coal combus￾tion． ISIJ Int，2009，49( 6) : 819 ［13］ Shen Y S，Guo B Y，Yu A B，et al． Three-dimensional model￾ling of coal combustion in blast furnace． ISIJ Int，2008，48( 6) : 777 ［14］ Guo B Y，Zulli P，Ｒogers H，et al． Three-dimensional simula￾tion of flow and combustion for pulverised coal injection． ISIJ Int，2005，45( 9) : 1272 ［15］ Fu W B． Coal Combustion Theory and Macroscopic General Ｒule． Beijing: Tsinghua University Press，2003 ( 傅维镳． 煤燃烧理论及其宏观通用规律． 北京: 清华大学 出版社，2003) ［16］ Han X L． Numerical Analysis． Beijing: Higher Education Press， 2011 ( 韩旭里． 数值分析． 北京: 高等教育出版社，2011) ［17］ Fu W B，Zhang B L． A simple and general method for calculat￾ing the burning rate of coal char /carbon particles in air． Combust Sci Technol，1993，89: 405 ［18］ Smith I W，Tyler Ｒ J． The reactivity of a porous brown coal char to oxygen between 630 and 1812 K． Combust Sci Technol，1974， 9( 3 /4) : 87 ［19］ Smith I W，Tyler Ｒ J． Internal burning of pulverized semi-an￾thracite: the relation between particle structure and reactivity． Fuel，1972，51( 4) : 312 ［20］ Zhang X W，Lu P S，Cao P F，et al． The kinetics of combustion of pulverized anthracite char in the temperature range 1200 to 2000 K． J Eng Thermophys，1991，12( 1) : 78 ( 张喜文，路丕思，曹鹏飞，等． 粉状无烟煤炭粒在 1200 － 2000 K 温度范围内燃烧的化学动力学特性． 工程热物理学 报，1991，12( 1) : 78) ［21］ Fu W B，Zhang B L． A generalized relationship between the ki￾netic parameters of char combustion and coalproperties． Combust Sci Technol，1997，3( 1) : 1 ( 傅维标，张百立． 煤焦燃烧反应动力学参数与煤种的通用 · 332 ·

·234· 北京科技大学学报 第36卷 关系.燃烧科学与技术，1997,3(1)：1) 24]Xu Y.Hu BS.Proper combustion of oxygen enrichment percent- 2]Liang JG,Yan HJ,Zhang J Y.Numerical simulation and opti- age with large amount of coal powder injection in blast furnace.J mization of gas particle turbulent flow and pulverized coal com- Iron Steel Res,2006,18(1)6 bustion in the lower zone of blast furace.Res fron Steel,2012, (许莹，胡宾生.高炉大喷煤量与富氧率的合理搭配.钢铁 40(4):14 研究学报，2006,18(1)：6) (梁建国，闫红杰，张家元·高炉下部气固湍流和煤粉燃烧 25]Zhao Z Q,Guo Y Q.Oxygen enrichment and coal injection in- 的数值模拟与优化研究.钢铁研究，2012,40(4)：14) dustrial experiment on No.I BF of ironmaking plant of Anyang 3]Bortz S,Flament G.Measurements of pulverized coal combustion Iron Steel Group Company.Res Iron Steel,2005(2):1 under conditions simulating a blast furnace environment.Iron- (赵正清，郭永强.安钢炼铁厂1号高炉宫氧喷煤工业试验 making Steelmaking,1982,10(5)222 钢铁研究，2005(2)：1)

北 京 科 技 大 学 学 报 第 36 卷 关系． 燃烧科学与技术，1997，3( 1) : 1) ［22］ Liang J G，Yan H J，Zhang J Y． Numerical simulation and opti￾mization of gas particle turbulent flow and pulverized coal com￾bustion in the lower zone of blast furnace． Ｒes Iron Steel，2012， 40( 4) : 14 ( 梁建国，闫红杰，张家元． 高炉下部气固湍流和煤粉燃烧 的数值模拟与优化研究． 钢铁研究，2012，40( 4) : 14) ［23］ Bortz S，Flament G． Measurements of pulverized coal combustion under conditions simulating a blast furnace environment． Iron￾making Steelmaking，1982，10( 5) : 222 ［24］ Xu Y，Hu B S． Proper combustion of oxygen enrichment percent￾age with large amount of coal powder injection in blast furnace． J Iron Steel Ｒes，2006，18( 1) : 6 ( 许莹，胡宾生． 高炉大喷煤量与富氧率的合理搭配． 钢铁 研究学报，2006，18( 1) : 6) ［25］ Zhao Z Q，Guo Y Q． Oxygen enrichment and coal injection in￾dustrial experiment on No． 1 BF of ironmaking plant of Anyang Iron ＆ Steel Group Company． Ｒes Iron Steel，2005( 2) : 1 ( 赵正清，郭永强． 安钢炼铁厂 1 号高炉富氧喷煤工业试验． 钢铁研究，2005( 2) : 1) · 432 ·