基于Energiron的球团矿气体还原动力学及其机理

张建良，王振阳，邢相栋，刘征建，刘兴乐



基于Energiron的球团矿气体还原动力学及其机理

张建良，王振阳，邢相栋，刘征建，刘兴乐

(北京科技大学 冶金与生态工程学院，北京，100083)

Energiron工艺，采用线性拟合、相关系数对比、扫描电镜能谱分析(SEM-EDS)等研究方法，考察酸性氧化球团矿在特定还原气氛下的气基还原过程。研究结果表明：反应过程可分为3个阶段，第1阶段(0~5 min)控速环节为界面化学反应，表观活化能为28.8 kJ/mol；第2阶段(5~28 min)为界面化学反应和内扩散混合控速；第3阶段(28~34 min)控速环节为内扩散；同时得出不同温度下的有效扩散系数eff和反应速率常数rea，拟合出rea与温度的关系式；最后给出内扩散阻力和界面化学反应阻力随时间的变化趋势。

气基还原；球团矿；Energiron；动力学；线性拟合；微观结构

Energiron是在HYL-ZR基础上建立的气基直接还原工艺[1]，该工艺的特点是还原气的选择具有高度灵活性，可以使用天然气、煤制气或焦炉煤气等，不需要煤气重整炉，还原气可在还原反应炉内依靠金属铁的催化作用进行自重整；增加了CO2分离系统，显著降低了CO2排放量，有数据显示[2]，采用Energiron工艺的DR-EAF流程比传统的BF-BOF流程CO2减排40%~60%。近几年，随着焦煤资源的日渐匮乏和环保要求的逐渐提升，Energiron工艺越来越受到关注。截止到2012年，全球共有5条生产线在建，合计产能为965万t，加之已投产的Energiron工艺，预计总产能将超过1 815万t。我国天然气资源匮乏[3]，曾一度认为不适于发展气基直接还原。但Energiron气基还原工艺较符合我国非焦煤资源丰富、CO2减排压力较大的国情[4]，因此，对Energiron气基还原球团矿的探究很有必要。目前，已有一些研究者对气基还原球团矿的过程进行了动力学分析。王兆才等[5]研究了反应温度、还原气中H2含量对还原速率的影响。Prakash等[6]研究了软质复合球团和自熔性球团在直接还原过程中的动力学限制环节。李永全等[7]研究了氢碳比(H2与CO的体积比)为0.8和1.6的还原气条件下竖炉直接还原球团矿的动力学行为。但根据反应机理，将还原过程分段进行动力学的研究并不多。本文作者通过模拟Energiron工艺的还原条件，采用线性拟合、相关系数对比和SEM-EDS分析的方法，将气体还原球团矿的反应分段进行研究，从而确定各阶段反应机理，并得出一系列动力学参数和关系式。

1 实验

1.1 实验原料

实验所用氧化球团矿来自国内某大型钢铁厂，利用化学分析对本实验球团矿进行检测，其成分见表1。

表1 实验所用氧化球团矿的化学成分(质量分数)

由表1可知：实验所用球团矿全铁约为63%，有害元素S和P质量分数较低，但酸性脉石质量分数较高，SiO2和Al2O3的质量分数接近6.4%，直接还原后酸性脉石将超过8.5%，故经Energiron工艺处理，由该种球团矿生产的直接还原铁(DRI)适合用作转炉炼钢冷却剂或电炉熔分的部分原料。若热送电炉炼钢，需配加一定量的碱性溶剂，会产生较大渣量[8]。由于我国自产铁矿品位普遍偏低，二氧化硅等脉石质量分数较高[9]，故选用表1所示球团矿作为实验原料。

1.2 实验设备与方法

将球团矿置入恒温干燥箱中，在105 ℃下烘干4 h，充分脱除其中的吸附水，烘干后取粒度为10.0~ 12.5 mm的6个球团矿为1组，并称量各组球团矿的初始质量。

模拟Energiron工艺流程[10]，以(H2):(CO)=5:2为还原气配比，分别选取850，900，950和1 000 ℃为还原反应温度。程序还原炉以10 ℃/min的速度进行升温，整个升温过程通入N2，以充分排出恒温区内的O2。将装有球团矿的钼丝篮随炉升温至预定温度，待炉温稳定30 min后，将N2切换至H2和CO混合还原气，还原反应开始。通过电子天平连续测量和计算机记录，系统会每隔12 s记录1次球团矿在实验中的质量损失情况。当质量损失数据趋于稳定时，还原反应基本结束，可得不同温度下气体还原球团矿的质量损失率曲线。取还原温度为1 000 ℃，反应进行至 10 min的球团矿，对其中心切面进行SEM-EDS分析，得到SEM-EDS图。

有研究表明[11]，还原气体流速的增大有利于加快还原反应速率，但是当气流线速度大于0.05 m/s时，还原率随时间的变化不再受到还原气体流速增加的影响，即外部对流传质此时对还原速率不产生影响，实验所用管式炉刚玉管的内径为45 mm，还原气流量选取5 L/min，此时还原气线速度为0.052 m/s，气基还原实验装置见图1。

1—高压气瓶；2—过程控制系统；3—热电偶；4—钼丝篮和球团；5—电阻炉；6—排气烟道；7—电子天平

图1 实验装置示意图

Fig. 1 Experimental device diagram

2 实验结果

干燥后4组球团矿在还原前的初始质量见表2。

表2 不同还原温度对应各组球团矿的初始质量

式中：为还原前球团矿质量，g；为还原后球团矿质量，g；(TFe)为还原前球团矿中全铁质量分数，%；(FeO)为还原前球团矿中氧化亚铁质量分数，%。

通过式(1)，可以计算出各个时刻各组球团矿的还原度，结果见图2。

温度/℃：1—850；2—900；3—950；4—1 000

图2 不同还原温度下(H2):(CO)=5:2时

Fig. 2 Change of reduction degree with time at different temperatures in atmosphere of(H2):(CO)=5:2

由图2可知：当反应时间相同时，温度高的反应还原度较大，如反应进行至20 min时，850，900，950和1 000 ℃的还原度分别为73.29%，80.49%，85.68%和88.67%。但随着还原反应温度的提高，还原度增加的幅度是降低的，说明当温度提高到一定程度时，再增加温度对还原反应的影响不大。另外，高温下的还原度曲线首先达到平缓，说明升温提高了反应速度，缩短了到达反应终点所需的时间，但不同温度下各反应终点处的还原度相差不大，说明在(H2):(CO)=5:2的还原气氛下，温度更多的是影响反应速度，而不是反应接近终点的还原度。

将还原反应温度为1 000 ℃，反应时间为10 min的球团矿中心切面进行SEM-EDS分析[12]，结果见图3。

(a) SEM像；(b)~(f)分别为图3(a)中~点的EDS

图3 还原反应进行10 min时球团矿切面扫描电镜图片及能谱图

Fig. 3 Scanning electron microscopy picture and spectrum diagram of pellet section at 10 min

图3(a)中处为球团矿切面外缘，处为切面中心。由图3(a)可以看出：处较致密，而处较疏松，球团矿从中心到外缘致密度逐渐增大。另外，由图3(b)~(f)可以看出：至处的氧含量逐渐增大，说明球团矿的气体还原反应是从外向内逐步进行的，且还原程度越大的部位越致密，颗粒粒度越小，而还原程度越小的部位则越疏松，颗粒则越粗大。

3 反应动力学分析

定义球团矿的失氧量随时间的变化率d/d为反应的还原速率，计算公式如下：

式中：为还原速率，g/min；∆微小时间内球团矿质量的变化，g；Δ为微小时间，min。

由式(2)可得反应温度为1 000 ℃，还原气氛为(H2):(CO)=5:2的条件下，还原速率随时间的变化曲线及其导函数曲线，见图4。

1—1 000 ℃还原速率；2—1 000 ℃还原速率导数曲线

图4 还原速率与其导数和时间的关系

Fig. 4 Relationship between reduction rate and its derivative with time

由图4可知：在反应温度为1 000 ℃，(H2):(CO)=5:2的还原气氛下，还原速率曲线出现了2次转折，从其导函数曲线可以判断，2次变速发生的时间分别为反应后5 min和28 min。即在反应初始阶段，还原速率较大；在反应开始5 min后，还原速率曲线的斜率发生了变化，曲线下降趋势变缓；在反应进行至28 min时，曲线再次出现转折，还原速率缓慢降低；直到进行至34 min，还原反应基本结束。曲线出现转折说明是气体还原球团矿的控速环节发生了改变。

本实验还原气流量为5 L/min，气流线速度为0.052 m/s，大于外扩散控速的临界流速0.05 m/s，还原率随时间的变化不再受到气流速度增加的影响，即外部对流传质此时对还原速率基本不产生影响，故在本实验条件下，外扩散不应是该气固反应的限制性环节。

当内扩散是限制性环节时，有反应时间与转化率的关系[13]：

L1OCp的伪随机码序列的长度是4 092，周期是8 ms。伪随机码是由一个12级的线性反馈移位寄存器DC1和一个6级的线性移位反馈寄存器DC2模二和产生的，DC1和DC2的时钟速率都是0.511 5 MHz，DC1从第6、8、11个寄存器反馈，DC2从第1、7个寄存器反馈，移位寄存器的移动方向也都是从低位到高位。DC1的初始状态IS1是000011000101，DC2的初始状态IS2根据卫星编号不同从000000到111111。

当界面化学反应是限制性环节时，有反应时间与转化率的关系：

当内扩散与界面化学反应混合控速时，反应时间与转化率的关系为

式中：为反应时间，min；f1为内扩散为限制性环节时的完全反应时间，min；f2为界面化学反应为限制性环节时的完全反应时间，min；f2为完全反应时间，min；为转化率；为界面反应的相对穿透深度。

在0~5，5~28和28~34 min 3个阶段分别作[1−3(1−)2/3+2(1−)]和[1−(1−)1/3]对时间的曲线，并对各个曲线进行线性拟合，比较其线性相关系数，从而确定3个阶段反应的限制性环节[14]。

图5~8中对动力学机理函数线性拟合的相关系数见表3。

1—内扩散；2—化学反应；3—内扩散线性拟合；4—化学反应线性拟合

图5 第1阶段内扩散与界面反应动力学机理函数及其拟合曲线

Fig. 5 Internal diffusion and interface reaction kinetics mechanism function and its fitting curve in the first phase

1—内扩散；2—化学反应；3—内扩散线性拟合；4—化学反应线性拟合

图6 第2阶段内扩散与界面反应动力学机理函数及其拟合曲线

Fig. 6 Internal diffusion and interface reaction kinetics mechanism function and its fitting curve in the second phase

1—内扩散与化学反应混合控速；2—内扩散与化学反应混合控速拟合

图7 第2阶段内扩散与界面反应混合控速函数曲线及其拟合曲线

Fig. 7 Hybrid control function curve of internal diffusion and interface reaction and its fitting curve in the second phase

1—内扩散；2—界面反应；3—内扩散线性拟合；4—界面反应线性拟合

图8 第3阶段内扩散与界面反应动力学机理函数及其拟合曲线

Fig. 8 Internal diffusion and interface reaction kinetics mechanism function and its fitting curve in the third phase

表3 动力学机理函数拟合的相关系数

通过图5所示的界面反应控制与内扩散控制的曲线对比可知：在第1阶段，界面化学反应为限制性环节的动力学机理函数的线性相关度较好，说明在气体还原球团矿的初始阶段，界面化学反应为主要的限制性环节。如图6所示，在第2阶段，内扩散控制曲线与界面化学反应控制曲线均未呈现良好的线性关系，说明在此阶段，还原反应并不存在唯一的限制性环节，即并不是单纯的界面化学反应或是内扩散对反应进行控速。对第2阶段内扩散与界面化学反应混合控速动力学方程进行拟合，结果见图7，发现线性相关度较好，说明在反应进行至第2阶段时，是由内扩散与界面化学反应混合控速。图8则表明：在第3阶段，内扩散为限制性环节的动力学机理函数的线性相关度较好，表明在该阶段的控速环节为内扩散。这是由于在反应初始阶段，气固反应在球团矿的表面进行，此时并没有生成产物铁层，内扩散的阻力为零，界面化学反应是限制性环节。随着还原反应的进行，由H2和CO还原生成的Fe层和还原中间产物层(FeO和Fe3O4)逐渐增厚，导致内扩散的阻力增大，而界面化学反应的阻力相对减少，故此时由界面化学反应和内扩散混合控速。在反应后期，由SEM可知：产物层逐渐增厚且致密度增加，使得还原气H2和CO以及产物气体CO2和H2O扩散更加困难，故此时内扩散成为反应限制性环节。

内扩散为限制性环节时的完全反应时间[15]：

界面化学反应为限制性环节时的完全反应时间为

参照上述1 000 ℃时的计算过程，可以分别得到850，900和950 ℃有效扩散系数eff和化学反应速率常数rea，结果见表4。

表4 不同温度时的Deff和krea

由Arrhenius公式

对式(9)积分可得：

式中：为化学反应速率常数；a为活化能。将不同温度条件下的化学反应速率常数的对数ln对温度的倒数作图，结果见图9。

1—反应速率常数与温度关系；2—反应速率常数与温度关系拟合

图9 反应速率常数与温度关系

Fig. 9 Relationship between reaction rate constant and temperature

由图9拟合可得该反应条件下的表观活化能a=28.8 kJ/mol，由此可得还原气氛为(H2):(CO)=52时化学反应速率常数与温度的关系式，即

在Fe2O3整个还原过程中，FeO还原为Fe这一步最为困难，故在计算反应平衡常数时，若只考虑FeO还原为Fe的反应[16]，则有：

式中：i为内扩散阻力；c为化学反应阻力；0为球团矿初始半径；为反应平衡常数。

由式(12)和(13)，可得各控速环节相对阻力[17]为：

由式(12)~(15)，可得内扩散和化学反应相对阻力随还原度的变化曲线，结果见图10。

1—内扩散相对阻力i相对；2—化学反应相对阻力c相对

图10 还原反应中内扩散与化学反应相对阻力的变化

Fig. 10 Changes of relative resistance of internal diffusion and chemical reaction in reduction reaction

由图10可见：在反应初始阶段(还原度为0~52%)，界面化学反应相对阻力较大，而内扩散相对阻力较小。这是由于在初始阶段，反应在球团矿表层进行，产物层较薄，故初始阶段内扩散阻力较小，界面化学反应为控速环节；随着反应的进行(还原度为52%~96%)，还原产物Fe层逐渐增厚，导致还原气CO和H2以及产物CO2和H2O扩散困难，内扩散阻力逐渐增大，故在此阶段，反应由界面化学反应和内扩散混合控速；在反应后期(还原度为96%~100%)，由于产物Fe层继续增厚且致密度较高，使得还原气CO和H2以及产物CO2和H2O扩散进一步困难，故此时内扩散成为阻碍反应进行的主要原因。

4 结论

1) 基于Energiron法，在(H2):(CO)=5:2的还原气氛下，还原反应可分为3个阶段：第1阶段(0~5 min)控速环节为界面化学反应，活化能为28.8 kJ/mol；第2阶段(5~28 min)为界面化学反应和内扩散混合控速；第3阶段(28~34 min)控速环节为内扩散。

2) 在(H2):(CO)=5:2的还原气氛下，升温可以提高球团矿的还原反应速度，但反应终点处的还原度相差不大，且随着还原反应温度的提高，还原度增加的幅度越来越小。

3) 在(H2):(CO)=5:2的还原气氛下，850，900，950和1 000 ℃的有效扩散系数eff分别为2.81×10−3，4.04×10−3，5.86×10−3和6.77×10−3m2/s；速率常数rea与温度的关系为。

4)随着还原反应的进行，内扩散阻力逐渐增大，而界面化学反应阻力逐渐减少。这是由于产物层逐渐增厚且愈发致密，阻碍了还原气CO和H2以及产物CO2和H2O的扩散。

[1] Volpatti A, Franco B, Martinis A. Energiron direct reduction process: The forefront technology[C]// Association for Iron and Steel Technology. AISTech 2012 Iron and Steel Technology Conference and Exposition. Atlanta: Association for Iron and Steel Technology, 2012: 615−624.

[2] 特诺恩. 联合钢厂采用无碳排放的ENERGIRON直接还原方案减少温室气体的排放[J]. 世界钢铁, 2012(3): 1−8.
TE Nuoen. Decrease of GHG emissions through the carbon free emissions ENERGIRON DR scheme in integrated mills[J]. World Iron & Steel, 2012(3): 1−8.

[3] 齐渊洪, 钱晖, 周渝生, 等. 中国直接还原铁技术发展的现状及方向[J]. 中国冶金, 2013, 23(1): 9−14.
QI Yuanhong, QIAN Hui, ZHOU Yusheng, et al. Present status and trend of direct reduction technology development in China[J]. China Metallurgy, 2013, 23(1): 9−14.

[4] 陈凌, 张涛, 郭敏, 等. 气基竖炉中国钢铁发展的新方向[C]// 2012年全国炼铁生产技术会议暨炼铁学术年会文集. 无锡: 中国金属学会, 2012: 345−350.
CHEN Ling, ZHANG Tao, GUO Min, et al. Gas-based shaft furnace new developing direction of Chinese steel industry[C]// The National Ironmaking Technology Conference and Ironmaking Academic Annual Meeting in 2012. Wuxi: The Chinese Society for Metals, 2012: 345−350.

[5] 王兆才, 储满生, 赵庆杰, 等. 还原气氛和温度对氧化球团气基还原反应进程的影响第七届中国钢铁年会论文集. 北京: 中国金属学会, 2009: 419−427.
WANG Zhaocai, CHU Mansheng, ZHAO Qingjie, et al. Effects of reducting gas composition and temperature on gas reduction reaction process of oxidized pellets[C]// The 7th China Iron & Steel Annual Meeting Proceedings. Beijing: The Chinese Society for Metals, 2009: 419−427.

[6] Prakash S, Goswami M C, Mahapatra A K S, et al. Morphology and reduction kinetics of fluxed iron ore pellets[J]. Ironmaking and Steelmaking, 2000, 27(3): 194−198.

[7] 李永全, 张瑞祥, 杨天钧, 等. 两种还原气条件下球团矿竖炉还原过程的动力学研究[J]. 宝钢技术, 2000(2): 36−40.
LI Yongquan, ZHANG Ruixiang, YANG Tianjun, et al. Kinetic research of pellets reduction in shaft furnace by using two different reducing gases[J]. Bao-Steel Technology, 2000(2): 36−40.

[8] 胡俊鸽, 周文涛, 郭艳玲, 等. 气基竖炉直接还原球团矿技术的发展[J]. 烧结球团, 2012, 37(2): 40−44.
HU Junge, ZHOU Wentao, GUO Yanling, et al. Development of gas-based shaft furnace direct reduction pellets technology[J]. Sinter and Pelletizing, 2012, 37(2): 40−44.

[9] 冯伟, 孙燕, 熊发挥. 浅析我国铁矿石资源可持续发展战略[J]. 四川地质学报, 2012, 30(1): 1−4.
FENG Wei, SUN Yan, XIONG Fahui. The simple analysis on the sustainable development strategy of iron ore resources in China[J]. Sichuan Acta Geologica, 2012, 30(1): 1−4.

[10] 蒙特雷, 布特里奥. 采用ENERGIRON工艺利用焦炉煤气或合成气生产直接还原铁(DRI)[C]// 2007中国钢铁年会. 成都: 中国金属学会, 2007: 171−180.
Alessandro M, Pablo E D. The production of direct reduced iron by using coke oven gas or syngas on ENERGIRON[C]// The National Steel Conference 2007. Chengdu: The Chinese Society for Metals, 2007: 171−180.

[11] 刘建华, 张家芸, 周土平. CO及CO−H2气体还原铁氧化物反应表观活化能的评估[J]. 钢铁研究学报, 2000, 12(1): 5−9.
LIU Jianhua, ZHANG Jiayun, ZHOU Shiping. Assessment of apparent activation energies for reducing iron oxides by CO and CO−H2[J]. Journal of Iron and Steel Research, 2000, 12(1): 5−9.

[12] 陈耀明, 李健. 氧化球团矿中Fe2O3的结晶规律[J]. 中南大学学报(自然科学版), 2007, 38(1): 70−73.
CHEN Yaoming, LIJian. Crystal rule of Fe2O3in oxidized pellet[J]. Journal of Central South University (Science and Technology), 2007, 38(1): 70−73.

[13] Kumar M, Mohapatra P, Patel S K. Studies on the reduction kinetics of hematite iron ore pellets with noncoking coals for sponge iron plants[J]. Mineral Processing and Extractive Metallurgy Review, 2009, 30(4): 372−392.

[14] 黄柱成, 易凌云, 彭虎, 等. 焙烧温度对氧化球团矿性质及其气基直接还原过程的影响[J]. 中南大学学报(自然科学版), 2012, 43(8): 2889−2895.
HUANG Zhucheng, YI Lingyun, PENG Hu, et al. Effect of roast temperature on properties of pellets and its gas-based direct reduction[J]. Journal of Central South University (Science and Technology), 2012, 43(8): 2889−2895.

[15] Long H M, Li J X, Wang P, et al. Reduction kinetics of carbon containing pellets made from metallurgical dust[J]. Ironmaking and Steelmaking, 2012, 39(8): 585−592.

[16] 沈峰满, 施月循. 高炉内气固反应动力学[M]. 北京: 冶金工业出版社, 1996: 167−168.
SHENG Fengman, SHI Yuexun. Gas-solid reaction kinetics in blast furnace[M]. Beijing: Metallurgical Industry Press, 1996: 167−168.

[17] Ono Y, Joko T. Evaluation of mass transfer resistance through the gas phase in reduction kinetics of oxide pellets[J]. Journal of the Iron and Steel Institute of Japan, 1977, 63(12): 1785−1790.

and mechanism of gas-based reduction of pellets based on Energiron

ZHANG Jianliang, WANG Zhenyang, XING Xiangdong, LIU Zhengjian, LIU Xingle

(School of Metallurgical and Ecological Engineering, University of Science and Technology Beijing, Beijing 100083, China)

Based on Energiron process, the research methods of the linear fitting, the correlation coefficient comparing and scanning electron microscope and energy dispersive spectrometer (SEM-EDS) were used to investigate the gas-based reduction process of the acid oxidation pellets under the particular reducing atmosphere. The results show that the reaction process can be divided into three phases: the controlled speed link of first stage (0−5 min) is the interfacial chemical reaction of which the activation energy is 28.8 kJ/mol, the second stage (5−28 min) is the hybrid of interfacial chemical reaction and internal diffusion, and the controlled speed link of the third stage (28−34 mim) is internal diffusion. The effective diffusion coefficienteffand the reaction rate constantreaat different temperatures are calculated at the same time. At last, the variation tendency of resistance of internal diffusion and interfacial with time are given.

gas-based reduction; pellets; Energiron; kinetics; linear fitting; microstructure

TF554

A

1672−7207(2015)01−0041−08

2014−02−13；

2014−04−20

国家重点基础研究发展计划(973计划)项目(2012CB720400) (Project(2012CB720400) supported by Major State Basic Research Development Program of China (973 Program))

王振阳，博士研究生，从事铁矿直接还原、综合利用研究；E-mail: wangzheny@126.com

10.11817/j.issn.1672−7207.2015.01.007

(编辑 杨幼平)