高磷鲕状赤铁矿气化脱磷动力学研究.pdf
高磷鲕状赤铁矿气化脱磷动力学研究 ① 张 伟, 王 辉, 邢宏伟, 田铁磊, 李 杰 (河北联合大学 冶金与能源学院 河北省现代冶金技术重点实验室,河北 唐山 063009) 摘 要 为探索烧结过程中高磷鲕状赤铁矿的脱磷机理,采用综合热分析仪,在升温速率分别为 10、15、20 ℃ / min 的条件下,通过 与 Fe2O3的对比试验,对高磷铁矿进行了气化脱磷动力学研究。 结果表明气化脱磷反应在第2 失重阶段发生,且温度为850 ℃时, Ca5(PO4)3F 和脱磷剂反应开始,1 050 ℃左右,脱磷反应最剧烈。 采用 Ozawa 法计算了高磷铁矿反应的第 1、2 阶段和 Fe2O3反应的 第 2 阶段活化能,分别为 104.71,250.55 和 168.80 kJ/ mol,脱磷反应过程中克服能垒需要更高能量;气化脱磷反应机理函数符合二 维扩散 Valensi 方程。 关键词 鲕状赤铁矿; 气化脱磷; 活化能; 反应机理 中图分类号 TF046文献标识码 Adoi10.3969/ j.issn.0253-6099.2015.03.022 文章编号 0253-6099(2015)03-0079-04 Kinetic Study for Gasification Dephosphorization of High⁃phosphorus Oolitic Hematite ZHANG Wei, WANG Hui, XING Hong⁃wei, TIAN Tie⁃lei, LI Jie (College of Metallurgy and Energy, Hebei United University, Hebei Key Laboratory of Modern Metallury Technology, Tangshan 063009, Hebei, China) Abstract In order to explore the dephosphorization mechanism in the sintering process of high⁃phosphorus oolitic hematite, the kinetics of gasification dephosphorization of high⁃phosphorus iron ore were studied with heating rate at 10, 15 and 20 ℃ / min respectively, based on the contrast experiment with Fe2O3by using simultaneous thermal analyzer. The results showed that the gasification dephosphorization occurred at the second stage of weight loss, and Ca5(PO4)3F and dephosphorization agent began to react when the temperature increased to 850 ℃, with the reaction reaching a maximum at 1 050 ℃. The Ozawa method was used to calculate the activation energy at two stages of reaction of high⁃ phosphorus iron ore and the second stage of Fe2O3reaction, leading to the values of 104.71, 250.55 and 168.80 kJ/ mol. It is found that more energy will be taken to overcome the energy barrier in the gasification dephosphorization process and the reaction mechanism is in consistent with the two⁃dimensional diffusion Valensi equation. Key words oolitic hematite; gasification dephosphorization; activation energy; reaction mechanism 我国高磷鲕状赤铁矿资源丰富,已探明储量达 40 亿吨,其中鄂西⁃湘西北成矿区是我国“宁乡式”鲕状赤 铁矿的主要富集区,仅鄂西地区已探明储量就达 16.87 亿吨[1]。 高磷赤铁矿被认为是国内外难选矿石之一, 磷元素以氟磷灰石的形式存在于高磷鲕状赤铁矿 中[2-3],其含磷量高,脉石嵌布粒度极细,分选难度大, 目前无法大规模开发利用[4-5]。 国内外越来越多的科 研工作者投身到高磷矿的开发利用上,但综合考虑提 铁降磷、节能减排和成本等因素,效果均不太理想,高 磷鲕状赤铁矿脱磷始终是冶金行业公认的难题[6-10]。 铁矿粉的烧结过程一般不能脱磷,但考虑烧结原 料配碳后会存在局部的还原性气氛,含磷矿物有可能 被还原,如果在卤族元素添加剂作用下,形成含磷气体 产物,则气化脱磷在理论上是可行的[11-12]。 本文以烧结气化脱磷理论分析和热力学计算为基 础,采用综合热分析仪,根据 TG 和 DTA 曲线,对气化 脱磷反应动力学进行了探索性研究,计算得到了反应 活化能和脱磷反应机理函数。 ①收稿日期 2014-12-09 基金项目 国家自然科学基金资助项目(51274081) 作者简介 张 伟(1977-),女,河北沧州人,实验师,硕士,主要从事冶金节能与资源优化方面研究。 第 35 卷第 3 期 2015 年 06 月 矿矿 冶冶 工工 程程 MINING AND METALLURGICAL ENGINEERING Vol.35 №3 June 2015 1 试 验 1.1 试验原料 试验用铁矿粉为湖北宜昌高磷鲕状赤铁矿,煤粉 为阳泉无烟煤,Fe2O3、SiO2、SiC、Na2SO4、CaCl2为分析 纯。 高磷鲕状赤铁矿和煤粉的成分及工业分析见 表 1~2。 表 1 高磷鲕状赤铁矿化学成分(质量分数) / % TFeFe2O3FeOSiO2CaOMgOAl2O3PS 53.9645.748.316.984.050.255.101.370.074 表 2 阳泉煤的工业分析结果/ % 灰分(Aad)挥发分(Vad)固定碳(FCad)水分(Mad) 16.619.5773.822.44 1.2 试验装置及条件 采用北京恒久 HCT-4 综合热分析仪(TG⁃DTA), 样品质量约 80.00 mg(高磷铁矿 88.9%、SiO21.41%、 SiC 1.57%、Na2SO41.34%、CaCl21.36%、配碳量 4%), 混合均匀后放入 Φ6 mm 8 mm 刚玉坩埚中,通入氮 气和氧气的混合气体(N2∶ O2= 4 ∶1),气体流量为 55 mL/ min,采用不同的升温速率(10、15、20 ℃ / min) 由室温升至 1 350 ℃,同步记录 TG⁃DTA 曲线,进行动 力学分析。 为了更准确分析气化脱磷反应阶段动力 学,选取纯 Fe2O3代替高磷矿,其他条件不变,进行对 比试验(注原料为高磷铁矿的实验以下简称方案 1, 对比试验称为方案 2)。 1.3 试验原理 原料中添加 CaCl2后,经过烧结,赋存在氟磷灰石 中的磷转化为 PCl3气态形式,在抽风负压作用下外 排,从而起到脱磷效果。 用 Factsage 软件对气化脱磷 反应进行热力学计算 CaCl2 +SO 2(g)+O2(g) → CaSO4 +Cl 2(g) (1) T<1 144 ℃,ΔGθ (1)<0 2Ca5(PO4)3F+15C+9Cl2+9SiO2→ 9CaSiO3+6PCl3(g)+15CO(g)(2) T开=810 ℃ 2 试验结果与分析 2.1 差热和热重分析 对试验样品(方案 1)分别以不同的升温速率进行 差热⁃热重分析,所得 TG 曲线、TG⁃DTA 曲线如图 1~2 所示。 0 -2 -4 -6 -8 -10 -12 -14 300 500 700 900 1100 1300 TGmg , ;520 min ;515 min ;510 min 图 1 不同升温速率下反应的 TG 曲线 300 500 700 900 1100 1300 TGmg , 0 -2 -4 -6 -8 -10 -12 -14 0 -4 -8 -12 -16 -20 -24 DTAμν TG DTA 图 2 升温速率为 15 ℃/ min 时反应的 TG⁃DTA 曲线 从图 1 可见,随着升温速率增大,热重曲线很有规 律地向高温区偏移。 该反应有明显两段失重,第 1 个 阶梯段失重较缓慢,温度区间跨度大,其主要原因是高 磷鲕状赤铁矿中 Fe2O3被 C 还原为 Fe3O4,产生 CO 气 体并排出,且铁氧化物的还原伴随整个过程。 当温度 升到 850 ℃以上时,第 1 个阶梯反应没有完全结束,第 2 个阶梯段失重紧接发生,失重速度明显快于第 1 个 阶梯段。 由热力学知识可知,Ca5(PO4)3F 的脱磷反应 在此阶段发生,产生大量 PCl3气体,与 CO、CO2一起 排除,所以失重迅速,最大失重速率发生在 1 050 ℃ (由于升温速率较大,开始反应温度 850 ℃ 较理论计 算值 810 ℃有所偏移,属于正常现象)。 当温度达到 1 200 ℃后,失重曲线逐渐平坦,失重反应结束。 由图 2 可知,当升温速率为 15 ℃ / min 时,在第 1 个失重阶段,有一个尖锐的吸热峰,是 Fe2O3还原吸热 所致,在 1 050 ℃左右,第 2 个失重阶段出现一个更大 更尖锐的吸热峰,原因是在 1 050 ℃,氟磷灰石的脱磷 反应剧烈,吸收大量热量,产生 PCl3气体。 2.2 动力学分析 根据升温方式不同,动力学研究可分为等温法和 非等温法。 本试验采用非等温法,非等温法能在开始 到结束的整个温度范围内计算动力学参数,且试验样 品用量少。 因此,根据所得实验数据分别采用 Ozawa 法和 Kissinger 法处理,求解反应活化能。 08矿 冶 工 程第 35 卷 2.2.1 Ozawa 法 非等温、非均相条件下反应的动力 学方程为 dα dT = A β exp - Ea RT f(α)(3) 式中 α 为反应的转化率;β 为升温速率,K/ min;f(α)为反 应机理函数的微分形式;Ea为反应活化能,kJ/ mol;A 为 频率因子;R 为气体常数,8.314 J/ mol;T 为绝对温度,K。 Ozawa 法是对式(3)进行近似变换而得 lgβ = lg AEa RG(α) - 2.315 - 0.4567 Ea RT (4) 式中 G(α) 为f(α) 的积分形式,定义为 G(α)=∫ α 0 dα f(α); Ea可用如下方法求得在不同 β 下,选择相同转化率 α,则 lg AEa RG(α) 是定值。 选取不同升温速率下相同转化率的数据作图,用 最小二乘法进行线性拟合,根据直线斜率计算反应活 化能 Ea,方案 1 和方案 2 的活化能见表 3。 表 3 不同转化率下的反应活化能 转化率 α 方案1(第1个阶梯)方案1(第2个阶梯)方案2(第2个阶梯) Ea Ea Ea Ea Ea Ea 0.1092.01169.62225.16 0.20113.41168.55217.82 0.30138.00178.88215.61 0.40148.75181.73175.63 0.50142.36223.46157.65 0.60140.04104.71256.45250.55156.39168.80 0.70123.34336.58142.26 0.8083.02340.49135.75 0.8561.24350.98143.09 0.9057.17298.35139.84 0.9552.53242.01147.60 注Ea为反应活化能平均值;反应活化能的单位为 kJ/ mol。 从表 3 可以看出,高磷铁矿气化脱磷反应的第 1 个失重阶段 Ea为 104.71 kJ/ mol,第 2 个阶段 Ea为 250.55 kJ/ mol,第 2 阶段的反应活化能明显高于第 1 阶段;又根据纯 Fe2O3第 2 阶段 Ea为 168.80 kJ/ mol, 可以判断在失重第 2 阶段发生了气化脱磷反应。 由气 化脱磷反应的热力学计算可知,气化脱磷反应为吸热 反应,克服能垒使反应发生需要更高的能量,因此方案 1 第 2 失重阶段反应活化能是比较高的。 2.2.2 Kissinger 法 对式(3)进行Kissinger 法处理可得 ln β T2 max = ln RA Ea - Ea R 1 Tmax (5) 在方案 1 的 DTA 曲线上,每条曲线都有一些峰, 峰顶所在温度为 Tmax,不同升温速率下,Tmax会有所偏 移,将所得数据作图,见表 4 和图 3,用最小二乘法线 性拟合成一条直线,根据斜率即可计算出反应活化 能 Ea。 表 4 不同升温速率下的 ln(β/ Tmax)与 1/ Tmax值(方案 1 第 2 个阶梯) 升温速率 / (℃min -1 ) Tmax / K T -1 max / K -1 ln(β/ Tmax) 101 3190.000758150-12.06667321 151 3440.000744048-11.69876084 201 3570.000736920-11.43033105 -11.4 -11.7 -12.0 0.0007360.0007440.0007520.000760 ln β/Tmax T-1 K-1 max 图3 ln(β/ Tmax)与1/ Tmax线性拟合结果(方案1 第2 个阶梯) 根据 ln(β/ Tmax)与 1/ Tmax线性拟合所得直线斜率 计算可得,Ea为 244.67 kJ/ mol,与 Ozawa 法计算出的 活化能 250.55 kJ/ mol 基本保持一致。 2.2.3 反应机理函数 反应的动力学模型主要有化学 反应、扩散反应、相界反应和成核与生长等,常采用的微 分和积分形式的动力学机理函数见表 5[13]。 本文选取 失重最剧烈的第 2 阶段(脱磷反应阶段),用试验数据选 取表 5 中不同的动力学机理函数进行线性拟合。 表 5 常用化学反应机理函数 模型G(α)f(α)反应机理 化学 反应 F1-ln(1-α)1-α一级反应 F1.52[(1-α)-1/ 2 -1] (1-α)3/ 21.5 级反应 F2(1-α) -1-1 (1-α)2二级反应 扩散 控制 D1α21/2α 一维扩散, Parabolic 法则 D2(1-α)ln(1-α)+α-ln(1-α)-1 二维扩散, Valensi 方程 D3[1-(1-α)1/ 3]21.5(1-α)2/ 3[1-(1-α)1/ 3]-1 三维扩散, Jander 方程 相界 反应 R1α1一维 R21-(1-α)1/ 22(1-α)1/ 2二维,收缩圆柱体 R31-(1-α)1/ 33(1-α)2/ 3三维,收缩球体 成核 生长 A2[-ln(1-α)]1/ 22(1-α)[-ln(1-α)]1/ 2 随机核化, Aevrami⁃Erofeev 方程Ⅰ A3[-ln(1-α)]1/ 33(1-α)[-ln(1-α)]2/ 3 随机核化, Aevrami⁃Erofeev 方程Ⅱ 18第 3 期张 伟等 高磷鲕状赤铁矿气化脱磷动力学研究 Bagchi 等[14]提出用微分法和积分法相结合的方 法对非等温动力学数据进行处理以确定反应机理,常 见的非等温处理方法有 Freeman⁃Correll 法(式 6) 和 Coats⁃Redfern 法(式 7) ln dα/ dt f(α) = lnA - Ea RT (6) ln G(α) T2 = ln AR βEa - Ea RT (7) 把各种反应机理函数的微分式和积分式分别代入 式(6) ~(7),代入试验数据计算并进行线性拟合,结 果见表 6。 表 6 不同动力学机理函数的动力学参数 序号机理函数 Freeman⁃Correll 法Coats⁃Redfern 法 EarEar 1F1232.760.97182192.900.99177 2F1.5304.560.99012231.350.97675 3F2376.300.9833274.850.95766 4D1242.100.83819284.190.99707 5D2261.520.90723272.010.99962 6D385.200.42086361.810.99797 7R1127.780.75634131.320.99642 8R2161.000.85082159.570.99931 9R3184.940.91327170.100.99784 10A2125.530.8828785.660.99002 11A389.790.7704749.940.98761 据文献[12],当两种方法求解的 Ea和计算所得 Ea相近,且相关系数较大(r≈1),则可由此判断反应 机理函数。 根据表 6,反应机理函数 D2,即二维扩散 Valensi 方程所得活化能 261.52 和 272.01 kJ/ mol 与前 面计算所得反应活化能 250.55 和 244.67 kJ/ mol 接 近,相关系数较大,分别为 0.90723 和 0.99962。 由此 判断该反应的动力学机理函数为二维扩散 Valensi 方 程G(α)= (1-α)ln(1-α)+α,f(α)= -ln(1-α)-1。 3 结 论 1) 高磷铁矿在整个反应过程中有两段明显失重, 通过与 Fe2O3的对比试验可知,气化脱磷反应发生在 第 2 失重阶段,为强吸热反应,脱磷产物为 PCl3。 2) Ca5(PO4)3F 在 850 ℃左右发生脱磷反应,随 升温速率增加,失重曲线向高温区偏移,反应为强吸热 反应, 因而克服能垒需要 更 高 能 量。 Ozawa 法 和 Kissinger 法计算得到的反应活化能分别为 250.55 和 244.67 kJ/ mol。 3) 根据 Freeman⁃Correll 法和 Coats⁃Redfern 法可 知,气化脱磷反应机理函数为G(α)= (1-α)ln(1-α)+α, f(α)= -ln(1-α)-1。 参考文献 [1] 赵一鸣,毕承思. 宁乡式沉积铁矿床的时空分布和演化[J]. 矿床 地质,2004,19(4)350-362. [2] 雷 婷. 高磷鲕状赤铁矿中磷矿物的还原焙烧行为及铁磷分离技 术[D]. 长沙中南大学资源加工与生物工程学院,2012. [3] 赵延成. 建议不再沿用“胶磷矿”这一矿物惯用名称[J]. 化工地 质,1994(3)215. [4] 孙炳泉. 近年我国复杂难选铁矿石选矿技术进展[J]. 金属矿山, 2006(3)11-13. [5] 印万忠,丁亚卓. 我国年选铁矿石资源利用现状[J]. 有色金属, 2006(8)163-168. [6] Xu C Y,Sun T C. Mechanism of phosphorus removal in beneficiation of high phosphorous oolitic hematite by direct reduction roasting with dephosphorization agent[J]. Transactions of Nonferrous Metals Society of China, 2012,22(11)2806-2812. [7] 江礼科,梁斌,邱礼有,等. 氟磷灰石热炭还原动力学研究[J]. 成 都科技大学学报,1995(5)1-8. [8] Tang H Q,Guo Z C,Zhao Z L. Phosphorus Removal of High Phosphorus Iron Ore by Gas⁃Based Reduction and Melt Separation[J]. Journal of Iron and Steel Research International, 2010,17(9)1-6. [9] Xia W T,Ren Z D,Gao Y F. Removal of Phosphorus from High Phos⁃ phorus Iron Ores by Sslective HCI Leaching Method[J]. Journal of Iron and Steel Research International, 2011,18(5)1-4. [10] Lonkov K,Gaydardzhiev S,Correa A,et al. Amenability for Process⁃ ing of Oolitic Iron Ore Concentrate for Phosphorus Removal [ J]. Minerals Engineering, 2013,46119-127. [11] 田铁磊,张玉柱,李 杰. 钢渣还原脱磷技术研究现状及展望 [J]. 甘肃冶金,2012(5)4-6. [12] 张 伟,刘卫星,田铁磊,等. 脱磷剂对高磷钢渣气化脱磷率影响 的研究[J]. 烧结球团,2014(1)47-50. [13] 胡荣祖,史启祯. 热分析动力学分析[M]. 北京科学出版社, 2001. [14] Bagchi T P,Sen P K. Kinetics of densification of powder compacts during the initial of sintering with constant rates of heating.A thermal analysis approach ( Ⅱ) Haematite compacts [ J]. Thermochimica Acta, 1981(15)175-178. 28矿 冶 工 程第 35 卷