水氯镁石的热解机理及动力学

第6卷第5期过程工程学报Vol 6 No52006年10月The Chinese Journal of Process EngineeringOct.2006水氯镁石的热解机理及动力学黄小芳,吴玉龙,杨明德,胡湖生,党杰,张建安清华大学核能与新能源技术设计研究院,北京100084)摘要:采用热重差热质谱联用技术(TG-DIA-MS和热分析技术(rG-DA)研究了水氯镁石热解过程的反应机理和动力学行为.结果显示,水氯镁石热解分4个阶段,前2个阶段共脱除4个结晶水,后2步脱水和水解并存,第1步的分解反应属于球形对称相界面反应为控制步骤的R3机理,后3步属于 Avrami- Erofeev的成核及核成长为控制步骤的A1机理,4步表观活化能E1分别为81.13,1254,2304和1641kJ/mol,频率因子分别为295×1010,2.70×1043,4.07×1026和475×106s-1关键词:水氯镁石;热解机理;动力学中图分类号:O642.13文献标识码:A文章编号:1009-606X(2006)05-0729-051前言材料制备等过程中得到了较好的应用验证山.本工作采用热重差热联用仪(TG-DTA),用 Doyle方程和 Coats-镁是一种极为重要的有色金属,自然界的镁资源主 Redfern方程,代入16个常见固相热分解机理的动力学要以水氯镁石MgCl26HO形式存在,海水、井矿盐卤机制函数,研究了水氯镁石热分解反应机理及其动力学水及盐湖卤水都含有大量的水氯镁石.我国的镁资源十行为,为水氯镁石脱水的研究和工业化提供了基础数据金属镁主要通过电解无水氯化镁得到,水氯镁石必须先2实验脱水形成无水氯化镁后才能用来生产金属镁,因此水氯2.1试样选取镁石脱水成了镁资源利用中十分重要的一环.水氯镁石试样1为北京益利精细化学品有限公司提供的分不能直接受热脱水至无水氯化镁,前4个结晶水容易失析纯六水氯化镁,平均粒径0.31mm.试样2由试样1去,但脱除后2个结晶水时易发生水解.因此关于水氯研细得到,平均粒径0.8mm镁石脱水,包括机理及动力学的研究对镁资源的合理利2.2实验方法用及我国西部盐湖资源的开发至关重要.试样1的TG-DTA-MS测定采用德国 Netzsch公司目前国内外在水氯镁石热解脱水的机理和动力学生产的STA429C综合热分析仪与瑞士 Balzers公司生产方面做了很多研究陈建军等用示差扫描量热法的 ThermoStarTM质谱仪联机耦合进行.实验条件如下:研究了青海察尔汘盐湖水氯镁石脫水过程的动力学;测定气氛氮气,温度范围室温至300,升温速率Kish等四通过DTA,TG,DTG等方法研究发现,水氯镁5/min,试样用量31450mg,TG量程200mg,DTA石脱水反应的表观反应级数和表观活化能随着氯化镁量程1000uV,质量扫描速度1amus,系统真空度水合数的减少而增加; Novikov等研究了伏尔加矿区3.6×10+Pa.预先在各通道中设定不同的质量数,预定水氯镁石的脱水机理和热力学这些工作对热解过程的的检测离子见表总机理大都是根据毎步热分解阶段的失重推测得到,并表1预设定的检测正离子没有对热解过程中逸出的气体作进一步定性考察,因此Table 1 Expected positive ions for mass spectra所得结论的可靠性未知.此外,在研究热解失重动力学时,对失重的动力学机制仅假设为最简单的幂函数形式而未作广泛的检验,因此所得每步分解机理和动力学参数准确性不佳.本研究采用热重差热质谱联用仪(TG-DIA-MS)研究了水氯镁石的热解机理,定性检测试样2中国煤化工青工 EXSTA6000了每步热分解的逸出气体,从而更准确地解释了各步热型热重差热仪CNMH(氩气,温度范围解机理.该方法已在κrO2粉制备、碳纳米管SiO2复合室温至300,升温速率5/min,样品用量23.324mng,收稿日期:2005-11-09,修回日期:2006-01-13作者简介:黄小芳(1980-,男,浙江省瑞安市人,硕士研究生,化学工程专业;杨明德,通讯联系人,Te:010-89796088,E-mail:yangmd@mail.tsinghua.edu.cn730过程工程学报第6卷灵敏度TG1mg,DTA1mV3.2TGDA热分解曲线分析3结果与讨论试样2的TG-DTA曲线如图3所示DTA曲线的形状及各峰顶温度与文献基本相似,共有5个吸热峰,3.1 TG DIA-M测定结果及分析其中第2个吸热峰较弱, Kirsh等根据实验结果认为该试样1的TG-DTA曲线如图1所示,热解逸出气峰为氯化镁的熔化过程.图3的TG曲线在50-119之体的正离子质谱图见图2.由图可见,试样分解分3个间出现1个台阶,质量损失约为168%,与MgC26H2O台阶,第1台阶50-210,质量损失约为34.0%,与加热失去2个结晶水的质量损失理论值(17.7%)相当,相MgCl26H2O脱除4个结晶水的质量损失理论值(35.4%)应的DTA曲线在104出现1个吸热峰;在119-161相当,相应的DIA曲线出现吸热峰,MS图出现HO之间出现1个台阶,质量损失为195%,也与MgCl26H2O和H2O峰,而没有CI与HCI峰,由此表明这一台阶失失去2个结晶水相近,相应的DTA曲线在148出现1去4个结晶水,而不存在水解.第2台阶210-240个吸热峰;在161-195之间出现1个台阶,质量损失质量损失为12.49%,DIA曲线在相应位置岀现吸热峰,为10.3%,相应的DTA曲线在181出现1个吸热峰;MS图在此区间出现H2O,OH与HC峰,其中H2O在195-256之间出现1个台阶,质量损失为10.4%,离子强度是HCI离子强度的30多倍,可认为这一过程相应的DA曲线在235出现1个吸热峰.以脱水为主,存在微量水解.第3台阶240-289,质量损失为10.59%,DTA曲线出现吸热峰,MS图出现10H2O,OH与HCI峰,且H2O的强度与HC在同一数量级,由此说明这一过程脱水反应与水解反应共存三3工100图3试样2的热重分析曲线Fig 3 TG-DTA curves of the sample 2-0.6100150200250300根据试样的TG-DTA-MS实验和TG-DTA实验测Temperature(℃)定结果,可以推测六水氯化镁热解分4个阶段,前2个图1试样1的热重分析曲线阶段每步失去2个结晶水,第3阶段以脱水为主,存在Fig. I TG-DTA curves of the sample 1微量的水解,第4阶段水解和脱水并存.以上结果与Kirsh等以及陈建军等的结论基本一致,唯一的差异6是本研究的第3阶段存在微量水解.因本研究采用热重:4a差热质谱联用仪研究了水氯镁石的热解过程,用质谱仪定性检测了每步热分解的逸出气体,从而更准确地解释了各步热解过程0由图1和3可知,试样1与2的TG曲线的失重率相当,分别为57.3%和57.6%,形状也相近,两者的DTA曲线峰也都有5个吸热峰.但图1的DTA曲线的吸热峰250300温度有滞后现象,原因可能是由于样品1的用量和粒度大于样品2中国煤化工实验条件也有所图2试样1热分解逸出气体的质谱图不同.根据CNMHG响,样品2的Fig 2 Mass spectrograms of gases evaporated fromTG-DTA实验结果更准确thermal decomposition of the sample 1第5期黄小芳等:水氯镁石的热解机理及动力学7314热分解反应机理及动力学分析为从试样2的单条TG曲线对固相分解反应非等温动力学进行研究,本研究采用 Doyle机理方程根据非等温反应动力学理论,线性升温条件下固相5.3305-1.052物质的分解反应动力学方程为In F(a)=In AEf(a)(1)式中[∥,由式(2)可知mFa对mT作图应是条直线.取16个常见固体热分解反应机理121的F(a式中β为线性升温速率(min,a为在温度TK时的反值(见表2)分别代入式(2),用最小二乘法求其线性相关应分解分数,A为频率因子mn),E为活化能(Jmo),系数同时,从直线的斜率求得反应活化能E,从截距f0x).动力学机制函数求得反应频率因子A表2常见固体热分解反应机理Table 2 Most frequently used mechanisms of solid state processes2-Dimensional diffusion (Valensi Eq.[In(1-a)x+(1-)ln(1-∞)3-Dimensional diffusion(sphere, Jander Eq.1.5(1-∞31-(1-ax)31-(1-ax)13]23-Dimensional diffusion(cylinder, G-B Eq.5[(1-∞)-1(1-20/3)-(1-0)2a(1-)Nucleationand nuclei growth( Avrami-Erofeev Eq,n=1)-ln(1-x)h(Avrami-Erofeev Eq, n=1.5)A1.51.5(1-∞)-n(l-n(1-x)Nucleation and nuclei growth( Avrami-Erofeev Eq, n=2)2(1-∝)-n(l-ln(1-α)ExponentialExponential nucleation(Mample Eq.Exponential nucleation(Mample Eq.Power law(n=1.5)2(1-a)32(1-ax)-2为了进行对比,本研究也采用了 Coats- Redfern方核及核成长为控制步骤的Al机理.程求取动力学参数按相关性较好的 Doyle方程式(2)求取动力学参数,[fat]=In ar1-2rTIb(3)4个阶段的计算结果为:E1=8113kJ/mo,A1=2.95×10s-,E2=1254kJ/mol,A2=270×1013s-,E3=2304k/mo,式中2RTE的值若随温度变化很小,则在程序升温过程A3=4.07×103°s-1,E4=1641kJ/mol,A4=4.75×10s-.求中可视为常数,故lnF(α)T对1作图应是一条直线,得的热分解活化能在热分解温度范围内满足 Doyle方程从斜率可求得反应活化能,取反应过程的平均温度,从的假设20≤E/RD)≤60,说明采用 Doyle方程是合理的,截距中求得反应频率因子A.故相应的动力学方程式分别为应用以上两种积分方法是直接由实验数据α和T求da1dT=2.95×101°e813x0m3(1-a)算动力学参数,避开了由dadT求算可能引入的计算误算,是较为合理的day/dT=270×1015e-12540(1-a2)采用 Doyle方程和 Coats- Redfern方程代入和Tda3/dT=407×1026c23040m(1-a3)的实验数据,分别对MgCl26HO热分解反应动力学进行了参数计算,结果见表3,其中s为平均标准偏差da/dT=4.75×101e104108(1-a4)判断固相分解反应机理一般以相关性为主要判据,文献[5,刀]都假设热解反应机理为最简单的幂函数如相关性均较接近,可选取平均标准偏差及相关经验为形式f)=k(1-x)来计算动力学参数,而本工作采用辅助判据.表4为2种计算方法的最优分解反应机理及 Doyle方程和 Coats- Redfern方程,代入16个常见固相活化能与频率因子.可以看岀,ρ种计算方法结果较接热分解机理氯镁石热分解反近,但由于应用 Doyle方程的条件比 Coats- Redfern方应机理及其中国煤化最优的机理和动程更严格,拟合的相关性较Coas- Redfern法普遍有所力学参数,并LCNMHy,N1x们、ts- Redfern方程提高.从表可见,第一阶段的机理是以球形对称相界面的约束条件,确定了最优的机理和动力学参数,因此本反应为控制步骤的R3机理,后3个阶段的机理是以成工作得出的结果更为合理732过程工程学报第6卷表3据式(2)拟合的各种qz6HO热分解反应固相反应机理Syms Table 3 Calculated results using Eq (2)for different solid state reaction mechanisms for decomposition of MgC12-6H,0st pnase0.16103L56.51×100.98310.2820l80.10.99800.1070145.04.08×10190.23802.50×1020.16700.99890.08156.7×109207.67.11×10240.9993101×10100.99370.10100.023881.130.99610.0833l12.0142×10132.28×10121254Al.50.99860.037060.281.19×1080.99910.0293931×10350.02200.01830.148.22×100.99910.014641.8699×100.994665.7390.080.99460.04021.47×0.98310.070133×1030.99460.98310.04700.99460.020116.430.98310.035322.5242.521.56×1050.89250.34802.68×1030.92650.290085.04606×100Third phaseFourth phase0.96190.422032746.33×100.9850236.59.62×100.9832410.31.36×1040.1590293.6300×10300.993497870.1850193.76.24×1020.99449374×10139832205.1107×102307951.87×10140.14604.07×l0.99864.75×10160.990409730.998615.27.79×10120.99869.97×10486152×103PI0.2110163.72.58×10480.98500.1330l18.28l.843.74×10985096190.0704576×1030.98500.04443.76×103961940.92C1.50.117080.40120×10°0.9180.15202.98×11053l10.0表4据式(2)和(3)拟合的水氯镁石最优热分解机理及活化能Table 4 Calculated results using Eqs. (2)and(3) for the best fit mechanisms ofbischofite thermal decomposition and apparent activation energyEquation Phase Mechanism Frequency factor(s Correlation coefficient Standard deviation Apparent activation energy (k J/mol)First2.95×10100.999702381254407×10360.99040.146030.44.75×101.37×100024l0.9990(3)Second2.40×1015Thirdl.04×100.9984应表观活化能分别为2304和164.1kJ/mol,频率因子分5结论别为4.07×1026和475×106s在惰性气体保护下,温度区间由室温到300、线参考文献性升温速率为5/min时,六水氯化镁热解过程分4个1郭起海镁盐的开发与应用应用化工,199(3:1415阶段,前2步热分解过程失去4个结晶水,其机理分别2]沈朴,兰旭光.水氯镁石的利用:制备工业氧化镁卩.环境保护,为球形对称相界面反应为控制步骤的R3机理与成核及(3]杨宝刚,于佩志等.MgO在氯化镁电解过程中的电泳行核成长为控制步骤的AI机理,分解反应表观活化能分为[J.轻金中国煤化工别为8113和1254kJ/mo,频率因子分别为2.95×101014陈建军,陈CNMHG化脱水过程的动力学研究卩J和2.70×103s-.后2步热解过程中脱水和水解共存,其5陈建军,马培华,陈关城,等,示差扫描量热法研究水氯镁石脱机理都是成核及核成长为控制步骤的A1机理,分解反水过程的动力学门海湖盐与化工,200,29(5):4-7第5期黄小芳等:水氯镁石的热解机理及动力学733[6 Moscowtz H, Lando D, Cohen H, et al. Study on the Dehydrate10陆昌伟,张玉峰.SiC纤维/LAS复合材料的TG-DTA-MS研究JBischofite J). Ind Eng. Chem. Prod, 1978, 17(2): 156-160复合材料学报,1994,l1(3):38-42.「刀] Kirsh Y, Yariv s, Shoval s. Kinetic Analysis of Thermal Dehydration[于惠梅,陆昌伟,宁金威,溶胶-凝胶法制备碳纳米管SiO2复合nd Hydrolysis of MgCl2 6H,o by DTA and tG []. Thermal. Anal材料过程的TG-DSC-MS研究[.质谱学报,2004,25(2):65-681987,32:392-408l2]李峰,何静,段雪,等.α-磷酸锆的制备及热分解非等温动力学[8 Novikov V L, Abbakumov V G, vikhman S I Kinetic Parameters of研究J无机化学学报,1999,15(1:55-6Chemically Enriched Magnesium Oxide Sintering [] Ogneupory, [13] Gao X, Dollimore D. The Thermal Decomposition of Oxalates 26: A1993,4:30-34Kinetic Study of the Thermal Decomposition of Manganess (f)[9]陆昌伟,施剑林,奚同庚,等.化学法制备超细ZrO粉热分解过Oxalate Dihydrate [J]. Thermochim. Acta, 1993. 215: 47-63程的TG-DTA-MS研究门.材料研究学报,1994,8(2):144-148Mechanism and Kinetics of Thermal decomposition of BischofiteHUANG Xiao-fang, wU Yu-long, YANG Ming-de, HU Hu-sheng DANG Jie, ZHANG Jian-an(Institute of Nuclear and New Energy Technology, TSinghua University, Beijing 100084, China)Abstract: Reaction mechanism and kinetic behavior of bischofite dehydration were studied by means of the TG-DTA-Ms couplingtechnique and the TG-DTA technique. The results show that there are four steps in the thermal decomposition of bischofite, in the firsto steps, four crystallized waters were lost and in the last two steps combined dehydration and hydration took place. The decompositionaction of the first step belongs to the R3 mechanism with 3-dimensional phase boundary reaction as the control step, last three stepsbelongs to the Al mechanism that Avrami-Erofeev nucleation and nuclei growth as control steps. The apparent active energy of foursteps is: 81.13, 125.4, 230.4 and 164.1 kJ/mol respectively, and the frequency factor is 2.95x10 0, 2.70x10 5, 4.07x1026 and 4.75x106S-'Key words: bischofite; thermal decomposition mechanism; kinetics中国煤化工CNMHG