快速加热条件下碳酸钙分解动力学

快速加热条件下碳酸钙分解动力学
李佳容;朱建国;朱书骏;刘敬樟;李百航
【摘要】在自行搭建的多气氛竖直管式高温煅烧实验平台上进行碳酸钙高温煅烧实验,样品在推入炉膛后即升温至较高温度.由对样品失质量曲线的分析可知:在快速加热条件下,温度对碳酸钙分解反应有较大影响,在950～1 000℃内影响较大,在1 000～1 100℃内影响较小;通过对比动力学模型得知:在快速加热条件下,反应气氛为CO2-N2混合气氛(CO2、N2体积分数分别为25％、75％)、纯N2时碳酸钙分解反应最概然模型为三维随机成核和随后生长模型,且2种气氛的表观活化能E 分别为207.741、203.786 kJ/mol,说明CO2对分解炉中碳酸钙分解存在抑制作用,会延长碳酸钙的完全分解时间.
【期刊名称】《中国粉体技术》
【年(卷),期】2018(024)006
【总页数】7页(P1-7)
【关键词】碳酸钙;高温煅烧;最概然模型;活化能
【作者】李佳容;朱建国;朱书骏;刘敬樟;李百航
【作者单位】中国科学院大学,北京100049;中国科学院工程热物理研究所,北京100190;中国科学院大学,北京100049;中国科学院工程热物理研究所,北京100190;中国科学院大学,北京100049;中国科学院工程热物理研究所,北京100190;中国科学院工程热物理研究所,北京100190;中国科学院工程热物理研究所,北京100190
【正文语种】中文
【中图分类】TQ127.1+3
我国水泥产量连续30年居世界第一，水泥生产的主要原料是石灰石和黏土[1]。

石灰石中的主要成分碳酸钙在分解炉中经过高温煅烧生成氧化钙和二氧化碳，然后进入回转窑，进行高温煅烧。

碳酸钙在分解炉内的分解过程对水泥生产至关重要。

目前，针对碳酸钙分解动力学参数的实验研究方法主要有热重实验、管式炉实验、高温气固悬浮炉实验等。

热重分析方法是国内外学者最普遍采用的方法：在不同升温速率下得到样品失质量曲线，计算得出动力学模式函数、活化能和指前因子[2-5]。

针对O2-N2、CO2-N2氛围，反应模型研究集中在收缩圆柱体模型、随机成核模型和随后生长模型，表观活化能为180～210 kJ/mol[6-10]。

但是热重实验的实验条件为较低升温速率（5～20 K/min）和较小样品质量（20～50 mg），与实际分解炉内石灰石进入炉膛瞬间升温的复杂状态有较大差别。

常规水平管式炉实验的特点在于样品质量较大（3～5 g），实验过程中可以保持恒温，样品放入后瞬间到达炉膛温度，与分解炉实际情况较一致，但是只能通过样品的不同停留时间来研究分解特性，无法对一份样品进行持续研究[11-12]。

高温气固反应实验中，生料与煤粉可以在装置中实现类似分解炉中的动态过程[13-15]，但是这种高温气固反应实验台无法得到样品质量随时间的变化曲线，只是通过气氛中CO2的变化间接得到石灰石分解率。

而悬浮态高温一维炉通过假定柱塞流间接计算样品转化率，有一定的局限性[16]。

鉴于升温速率严重影响碳酸钙分解反应，本文中使用了自行搭建的多气氛竖直管式高温煅烧实验系统[17]，利用立式管式炉实现恒温，在将样品推入炉膛时迅速升温至恒定温度，与实际分解炉情况类似，并且使用在线电子天平对样品质量进行实时
记录，得到比较完整的样品失质量曲线，进而得到碳酸钙高温煅烧随着温度的变化特性及在特定气氛的表观活化能。

1 实验
1.1 材料
实验材料选用天津市致远化学试剂有限公司生产的碳酸钙分析纯（CaCO3质量分数≥99.0%），符合 GB/T 15897—1995 标准，粒径分布为 0～45 μm，其中
d50=7.3 μm。

1.2 装置
实验装置示意图如图1所示，主要由立式管式炉（炉膛、保温层、热电偶、坩埚、电阻丝）、计算机、电子天平和配气系统（气瓶组、减压阀、流量计、气体混合器）组成。

炉膛内径为Φ50 mm、高度为1 000 mm，炉膛周围布置4支均匀分布的电阻丝加热，温度由一支布置在炉膛中心的热电偶测定，实验过程中保持恒温。

实验气体由气瓶、转子流量计和一个气体混合器供应，可以保持实验过程中气氛恒定。

实验样品平铺于坩埚内，坩埚由刚玉管推入炉膛，置于电子天平上，电子天平与计算机连接，在线记录样品质量。

图1 实验装置示意图Fig.1 Schematic diagram of test device1—炉膛；2—保
温层；3—热电偶；4—坩埚；5—计算机；6—电子天平；7—气瓶；8—减压阀；9—流量计；10—气体混合器；11—电阻丝；12—热电偶。

1.3 条件
实验条件如表1所示，其中纯N2的体积分数为100%，CO2－N2混合气氛中
CO2的体积分数为 25%、N2的体积分数为75%；每次实验样品质量为
（0.5±0.001）g，保证了较大样品质量和快速升温实验条件。

表1 实验条件Tab.1 Test condition工况实验气氛实验温度/K 1 1 223 2 1 273 3 1 323 4 1 373 5纯N2 1 223 6 1 273 7 1 323 8 1 373 CO2-N2混合气氛
2 结果与讨论
2.1 温度对碳酸钙分解反应的影响
碳酸钙转化率曲线如图2所示。

图2 碳酸钙转化率曲线Fig.2 Calcium carbonate conversion curve
图2 a为实验气氛为纯N2条件下碳酸钙转化率随时间的变化曲线。

可以看出，样品分解速度随炉膛温度升高而增大，当炉膛温度处于1 000℃以上时，样品均在1 200 s内分解完全。

碳酸钙分解属于吸热反应，温度升高促进了分解反应的进行，化学反应式为
图2 b为CO2－N2混合气氛下不同温度的碳酸钙转化率随时间的变化曲线。

实际工程中，分解炉内CO2体积分数为20%～30%，与该实验条件接近。

由图可见，1 000℃以上时样品均在2 000 s内分解完成。

图3 不同煅烧温度下样品比孔容积和孔径变化Fig.3 Changes of pore volume and pore size under different calcination temperatures
实验结果表明：随着温度升高，反应速率增大，在950～1 000℃区间内，温度对反应速率有显著影响。

相同温度下，纯N2气氛完全分解时间较短；1 000℃时，纯氮气气氛碳酸钙在1 000 s内完全分解；而在CO2（体积分数为25%）存在时，完全分解需要1 500 s。

CO2的存在不利于碳酸钙分解反应的正向进行，会延长碳酸钙的完全分解时间。

对CO2－N2混合气氛实验前后的样品进行BET检测，不同煅烧温度下样品孔容
积和孔径变化结果如图3所示。

950℃之前，随着煅烧温度的升高，孔径和孔容积逐渐增大；当煅烧温度进一步增加时，孔径和孔容积逐渐缩小；950℃恒温煅烧时结构最为疏松；随着温度升高至1 000℃以上，达到碳酸钙熔融温度，出现过烧现象，不利于碳酸钙的分解，产物活性也会降低。

与陈海、于艳敏等[18-19]的结论
相近似。

2.2 碳酸钙分解反应的表观活化能
使用等温热分析方法分别对纯N2工况与CO2－N2混合气氛工况实验数据进行分析，其中碳酸钙分解动力学函数为
式中：G（α）为分解反应动力学函数，k为反应速率，t为反应时间，α为分解率，m为t时刻样品质量，mi为样品初始质量，me为样品最终质量。

根据对各种经典反应模型的分析，相边界反应模型Rn和随机成核及随后生长模型An较有可能适用于本实验条件下的碳酸钙分解机理。

5种分解机理函数如表2所示。

5种分解机理函数的函数值分布如表3所示。

使用约化时间图法来确定最概然机理函数。

令
则可得
将上述5个模型作为标准模型曲线，将不同温度下碳酸钙的α～t/t0.5曲线作为实验曲线，绘制在同一张图上。

如图4所示，图中点线为5个标准模型数据，曲线
为4个实验温度数据。

通过对比可以看出，实验中的2种氛围均较为符合随机成
核和随后生长模型。

表2 5种分解机理函数Tab.2 Five decomposition mechanism functions符号
名称函数f（x）函数G（x）R2 二维相界面反应，收缩圆柱体模型 2（1-α）1/2
1-（1-α）1/2 R3 三维相界面反应，收缩球状模型 3（1-α）2/3 1-（1-α）1/3
A2 二维随机成核和随后生长模型 2（1-α）[-ln（1-α）]1/2 [-ln（1-α）]1/2 A3 三维随机成核和随后生长模型 3（1-α）[-ln（1-α）]2/3 [-ln（1-α）]1/3 A4 四
维随机成核和随后生长模型 4（1-α）[-ln（1-α）]3/4 [-ln（1-α）]1/4
表3 5种机理函数的函数值分布Tab.3 Distribution of function values for five mechanism functions转化率αR2（α）R3（α）A2（α）A3（α）A4（α）0.10 0.051 317 0.034 511 0.324 593 0.472 309 0.569 731 0.20 0.105 573 0.071 682 0.472 381 0.606 543 0.687 300 0.30 0.163 340 0.112 096 0.597 223 0.709 182 0.772 802 0.40 0.225 403 0.156 567 0.714 721 0.799 388 0.845 412 0.50 0.292 893 0.206 299 0.832 555 0.884 997 0.912 444 0.60 0.367 544 0.263 194 0.957 231 0.971 280 0.978 382 0.70 0.452 277 0.330 567 1.097 257 1.063 830 1.047 500 0.80 0.552 786 0.415 196 1.268 636 1.171 902 1.126 338 0.90 0.683 772 0.535 841 1.517 427 1.320 500 1.231 839 0.99 0.900 000 0.784 557 2.145 966 1.663 726 1.464 912
图4 不同温度下碳酸钙的α～t/t0.5曲线Fig.4 α～t/t0.5curve of calcium carbonate at different temperatures
使用A2、A3、A4模型，得到2种氛围各个温度下的t～G（α）曲线，如图5所示。

图5 2种气氛不同温度下的t～G（α）曲线Fig.5 t～G（α）curves at different temperatures under two atmospheres conditions
通过线性拟合，其斜率为反应速率k，R2表征拟合度。

综合4种温度的拟合情况，R2表征线性拟合度，拟合度越大证明模型与实验结果符合程度越高。

综合来看，A3模型符合程度最高，因此认为本实验中的最概然模型为三维随机成核和随后生长模型。

表4为A3模型下2组氛围的温度T和反应速率k的对应值。

表4 三维随机成核和随后生长模型参数Tab.4 Parameters of the random nucleation and growth model（n=1/3）氛围t/℃ T/K A3模型反应速率k （×10-9）纯N2 950 1 223 199 949 1 000 1 273 961 417 1 050 1 323 1 440
000 1 100 1 373 1 930 000 CO2－N2混合气氛950 1 223 242 255 1 000 1 273 738 089 1 050 1 323 1 590 000 1 100 1 373 2 470 000
由Arrhenius定律可知，反应速率k为
式中：A为阿累尼乌斯常数，R为气体常数，T为反应温度。

两边分别取自然对数，可得
绘制2种气氛下A3模型的1/T～lnk曲线如图6所示，其斜率为-E/R，截距为lnA。

通过线性拟合，得到2种气氛的活化能E与阿累尼乌斯常数（以对数lnA表示），如表5所示。

纯N2气氛下的活化能E比CO2－N2混合气氛下活化能E小，可
以看出CO2抑制碳酸钙的分解。

图6 2种气氛下A3模型的1/T～lnk曲线Fig.6 1/T～lnk curve of A3model under two atmospheres conditions
表5 2种气氛下活化能与阿累尼乌斯常数Tab.5 Activation energy values and arrhenius constants under two atmospheres conditions气氛活化能E/（kJ·mol-1）lnA纯N2 203.786 11.854 CO2－N2混合气氛 207.741 12.261
3 结论
本文中研究快速加热条件下碳酸钙的分解特性后，主要结论如下：
1）煅烧温度越高，碳酸钙分解速率越快，且在950～1 000℃较低温度区间内温
度对分解速率的影响较大，在1 000～1 100℃温度区间内影响较小，CO2的存在会延长碳酸钙的完全分解时间。

2）通过对分解炉CO2－N2混合气氛下不同煅烧温度产物的BET比孔容积与孔径
的分析，得出950℃最有利于碳酸钙的分解，炉膛温度过高引起碳酸钙快烧现象，不利于分解且降低产物活性。

3）在模拟分解炉CO2－N2混合气氛下与对比工况纯N2中，碳酸钙的分解符合
三维随机成核和随后生长模型。

由于升温速率的差异，因此与热重分析方法所得出的相界面模型不同[8]。

反应气氛为CO2－N2混合气氛时碳酸钙分解反应的表观
活化能为 207.741 kJ/mol，纯 N2气氛下表观活化能为203.786 kJ/mol。

可以
看出，CO2对碳酸钙分解存在一定抑制作用。

参考文献（References）：
【相关文献】
[1]王小龙.水泥生产过程中汞的排放特征及减排潜力研究[D].杭州：浙江大学，2017.
[2]田浩臣.高CO2气氛中模拟分解炉条件下的NO生成特性研究[D].广州：华南理工大学，2017.
[3]JIN D，YU X，YUE L，et al.Decomposition kinetics study of AlOOH coated calcium carbonate[J].Materials Chemistry & Physics，2009，115（1）：418-422.
[4]张保生，刘建忠，周俊虎，等.粒度对石灰石分解动力学影响的热重实验研究[J].中国电机工程学报，2010，30（2）：50-55.
[5]ANDERSON H L，KEMMLER A，HOEHNE G W H，et al.Round robin test on the kinetic evaluation of a complex solid state reaction from 13 european laboratories：part1.kinetic TG-analysis[J].Cheminform，1999，332（1）：33-53.
[6]胡贞武，李相国，马保国，等.水泥生料颗粒热分解机理研究[J].水泥工程，2007（6）：4-6，18.
[7]AR I，DOGU G.Calcination kinetics of high purity limestones[J].Chemical Engineering Journal，2001，83（2）：131-137.
[8]王世杰，陆继东，胡芝娟，等.水泥生料分解动力学的研究[J].硅酸盐学报，2003，31（8）：811-814.
[9]CHEN H，LIU N.Application of non-arrheniusequationsin interpreting calcium carbonate decomposition kinetics：revisited[J].Journal of the American Ceramic Society，2010，93（2）：548-553.
[10]CHEN H，LIU N，ZHAO W.Critical study on the identification of reaction mechanism
by the shape of TG/DTG curves[J].Solid State Sciences，2010，12（4）：455-460.
[11]郑瑛，陈小华，周英彪，等.CaCO3分解动力学的热重研究[J].华中科技大学学报（自然科学版），2002，30（8）：71-72.
[12]XIE Y，SUN S，LIU J，et al.The effect of additives on migration and transformation of gaseous pollutants in the vacuum pyrolysis process of waste printed circuit
boards[J].Waste Management & Research the Journal of the International Solid Wastes & Public Cleansing Association Iswa，2016，35（2）：190-199.
[13]ESCARDINO A，GARCIA-TEN J，FELIU C，et al.Kinetic study of the thermal decomposition process of calcite particles in air and CO2，atmosphere[J].Journal of Industrial & Engineering Chemistry，2013，19（3）：886-897.
[14]吕刚.水泥分解炉内NO生成和还原机理的实验及模拟研究[D].武汉：华中科技大学，2011.
[15]王世杰.水泥预分解窑系统内生料分解、煤粉燃烧与NOx控制研究[D].武汉：华中科技大学，2006.
[16]王俊杰，颜碧兰，汪澜，等.高温悬浮态下碳酸钙分解反应动力学的研究[J].武汉理工大学学报，2013，35（5）：41-44.
[17]刘敬樟，李百航，朱建国，等.焦炭颗粒群氧燃烧速率特性[J].中国科学：技术科学，2017，47（12）：1265-1272.
[18]陈海，张世红，杨海平，等.石灰石煅烧产物活性分析[J].无机盐工业，2013，45（3）：12-14.
[19]于艳敏，毕万利，蔺德忠，等.石灰石晶形结构对煅烧石灰活性的影响[J].辽宁科技大学学报，2015，38（2）：88-92.。