《武汉工程大学学报》  2017年02期 134-140   出版日期：2017-05-04   ISSN:1674-2869   CN:42-1779/TQ

---

碱式碳酸锆的热分解非等温动力学及机理研究


YUAN H， LI H F， REN Y G. Non-isothermal kinetics and mechanism of thermal decomposition of basic zirconium（Ⅳ） carbonate［J］. Journal of Wuhan Institute of Technology，2017，39（2）：134-140.氧化锆因具高硬度、高耐磨、低摩擦系数、高离子电导率、高弹性模量、高化学惰性、低热导率、高熔点等物理化学性能［1］，而广泛应用于先进陶瓷、固体电解质、热障涂层、电-光材料、气体传感器、催化剂等领域［2］. 常压下，氧化锆存在3种晶型［3］：单斜相（m-ZrO2，P21/c，＜1 170 ℃）、四方相（t-ZrO2，P42/nmc，1 170 ℃~2 370 ℃）、立方相（c-ZrO2，Fm3m，＞2 370 ℃）. 氧化锆粉体的粒径、形状、粒度分布等对氧化锆陶瓷制品的性能有重要影响［4］. 但纳米粉体因具很高的比表面积而致使晶粒易团聚. 研究表明：团聚体的相对强度及密度可很大程度地影响粉体成型密实程度及烧结体性能［5］. 因此，制备纳米氧化锆粉体除要求具有高活性、高比表面积、高纯度、热处理温度相对低等，还须防止超细晶粒团聚. 即使发生团聚，团聚体晶粒间也须具较弱的束缚力以使粉末成型过程容易被破坏. 制备纳米氧化锆粉体的方法通常有溶胶-凝胶法［6］、水热或熔剂热法［7］、微乳液沉淀法［8］、表面活性剂诱导的非水解/非溶胶-凝胶热解合成［9］、微波辐射法［10］、超声辅助/水热合成法［11］、化学气相沉积［12］、水解法［13］、热分解法［14］等. 其中，热分解是制备纳米氧化锆比较方便、实用、有效的方法. 该法具有工艺简单，成本低，热处理温度低等优点. 文献［2］报道了热分解Zr2O3（OH）2?H2O、Zr（OH）4、Zr3（C6H5O7）4?4H2O、ZrO2-x（OH）2x、C6H5O7Zr、ZrOCl2?8H2O、ZrOx（OH）4-2x等含锆化合物制备纳米氧化锆粉体，其动力学研究仅限于无定形ZrO2晶化动力学及亚稳的t-ZrO2→m-ZrO2相变动力学. 碱式碳酸锆（Basic Zirconium（Ⅳ） Carbonate， ZrCO3（OH）2?ZrO2，CAS：57219-64-4）是一种重要的含锆化合物，广泛应用于建材、皮革、涂料、造纸、化妆品、油漆、陶瓷、脂肪酸盐等诸多领域. 它既是制备碳酸锆铵、碳酸锆钾等衍生化合物的重要原料，也是制备纳米氧化锆粉体比较经济的前驱体原料. 所以，其热分解机理研究具有重要意义，但至今未见碱式碳酸锆热分解动力学及机理的研究. 本研究采用热分析技术研究碱式碳酸锆热分解非等温动力学，解析其热分解过程，研究其非等温热分解动力学及热分解反应机理. 1　实验部分1.1　试剂与材料碱式碳酸锆（ZrCO3（OH）2?ZrO2，阿拉丁试剂（上海）有限公司，ZrO2质量分数≥40%），综合热分析仪（德国耐驰仪器制造有限公司，NETZSCH STA 409 PC/PG型）. 温度范围30 ℃ ~1 200 ℃，动态氮气与氧气吹扫的流量分别为20 mL/min、10 mL/min. 采用Al2O3刚玉坩埚，以1 550 ℃煅烧5 h的a-Al2O3为参比试样. 2.2　实验方法取10 mg左右碱式碳酸锆粉末样品，采用综合热分析仪对样品进行TG-DTG-DSC（thermal gravity-differential thermal gravity-differential scanning calorimeter）分析，非等温热分解升温速率β分别为10 ℃/min、20 ℃/min、30 ℃/min，终点温度为1 200 ℃，根据固相反应动力学理论及TG-DTG数据，计算热分解动力学相关参数. 2　结果与讨论2.1　TG-DTG分析图1为碱式碳酸锆不同升温速率的TG-DTG曲线. ZrCO3（OH）2?ZrO2的非等温热分解过程分为4个阶段：室温~229 ℃ （第一阶段）、229 ℃~503 ℃（第二阶段）、503 ℃~642 ℃ （第三阶段）、642 ℃~ 1 200 ℃（第四阶段）. 以10 ℃/min、20 ℃/min、30 ℃/min升温速度的失重率分别为19.27%、20.47%、20.89%，与碱式碳酸锆［ZrCO3（OH）2?ZrO2］理论失重率20.09%相吻合. 2.2　非等温动力学研究2.2.1　迭代法求反应活化能　根据非等温反应动力学理论，线性升温条件下的固态物质热分解动力学方程［15］为：[dαdT=Aβ?exp(-EaRT)?f(α)] （1）式中，[α]为转化率（%），[β]为加热速度（[β]=dT/dt，℃?min-1），[Ea]为表观活化能（kJ?mol-1），[A]为指前因子，[f(α)]为转化函数，R为摩尔气体常数（R=8.314 J?mol-1?K-1），T为对应温度（K）. 根据公式（1），经系列变形可得Kissinger-Akahira-Sunose （KAS）方程［16］与Flynn-Wall-Ozawa （FWO）方程［16］，分别为：KAS方程 [lnβT2=lnARg(α)Ea-EaRT] （2）FWO方程[lnβ=ln0.004 8AEag(α)R-1.051 6EaRT] （3）为求得更接近实际值的表观活化能[Ea]，通常采用迭代法求解. KAS方程与FWO方程对应迭代方程［17］分别表示为：KAS迭代方程 [lnβh(x)T2=lnARg(α)Ea-EaRT] （4）FWO迭代方程 [lnβH(x)=ln0.004 8AEag(α)R-1.051 6EaRT] （5）式中[x=Ea/RT]，[g(a)]为最可机函数积分式. 其中[h(x)]、[H(x)]表达式［17］为：[h(x)=x4+18x3+88x2+96xx4+20x3+120x2+240x+120] （6）[H(x)=exp(-x)h(x)/x20.004 8exp(-1.051 6x)] （7）迭代的步骤为：（ⅰ） 假定h（x）=1或H（x）=1，判断活化能初始值[Ea1]；（ⅱ） 利用[Ea1]，根据[ln[β/H(x)]~1/T]或者[ln[β/(h(x)T2)]~1/T]关系计算活化能初始值[Ea2]；（ⅲ） 重复步骤（ⅱ），利用[Ea2]代替[Ea1]，直到满足[Eai-Ea(i-1)≤0.01]kJ?mol-1为止，[Eai]就是固相热分解表观活化能的准确数值；由于表观活化能的估值无需选择特定动力学模型［即特定[g(α)]函数］，故迭代法求得的表观活化能[Ea]较传统等转化率法更可信. 根据公式（2）与公式（3），ZrCO3（OH）2?ZrO2非等温热分解各阶段不同转化率α对应[lnβ~1 000/T]（FWO法）、[ln[β/T2]~1 000/T]（KAS法）关系见图2与图3所示. 利用最小二乘法线性回归分析得ZrCO3（OH）2?ZrO2非等温热分解各阶段、不同转化率[α]条件下活化能[E1a]（FWO法）、[E2a]（KAS法）. 再分别以[E1a]、[E2a]为初始值，利用公式（4）与公式（5），采用迭代法求得相应转化率α条件下活化能[E3a]（FWO迭代法）、[E4a]（KAS迭代法），所求[E1a]、[E2a]、[E3a]、[E4a]结果见表1所示. 利用不同计算方法得到的ZrCO3（OH）2?ZrO2非等温热分解各阶段不同转化率[α]条件下活化能的相对误差REai （%）见表1所示. 由表1结果可以看出，FWO法计算所得活化能值[E1a]稍高，而KAS法计算所得活化能值[E2a]与迭代结果比较接近. 同时可知，对于同一转化率α、以不同活化能初始值[E1a]或[E2a]，利用[ln[β/H(x)]~1/T]或者[ln[β/(h(x)T2)]~1/T]迭代求得活化能[E3a]、[E4a]是相等的. 说明KAS法及迭代法计算得活化能值更具可靠性. 故，0.2≤[α]≤0.8条件下，ZrCO3（OH）2?ZrO2非等温热分解各阶段活化能为94.963 kJ?mol-1 （第一阶段） 、309.781 kJ?mol-1 （第二阶段）、362.591 kJ?mol-1 （第三阶段）. [E1a]由FWO方程求得. [E2a]由KAS 方程求得. [E3a]由FWO方程所得迭代结果. [E4a]由KAS 方程所得迭代结果. [REia] = 100*（[Eia]-[E4a]）/[E4a]. BKAS由公式（4）所得截距. BFWO 由 公式（5）所得截距. lnAKAS 由 公式（4）所得指前因子. lnAFWO 由公式（5）所得指前因子. 3.2.2　最可机函数判定　通过方程（8）［18］判定ZrCO3（OH）2?ZrO2热分解动力学的最可机函数[g(α)]：[lng(α)=[lnAEaR+lne-xx2+lnh(x)]-lnβ] （8）将各升温速率[β]下，对应同一温度的转化率α带入方程（8），30种机理函数[g(α)]代数表达式参见文献［18］. 作[lng(α)~lnβ]关系图，利用最小二乘法线性回归，如反应机制符合某一特定[g(α)]函数，则直线斜率应趋近等于-1.000 0，同时其线性相关系数R2较高. ZrCO3（OH）2?ZrO2非等温热分解加热速度[β]=20 ℃?min-1时各阶段以[lng(α)~lnβ]关系线性回归求得的截距B、斜率k及相关系数R2见表2所示. 由表2可知，ZrCO3（OH）2?ZrO2非等温热分解各阶段所对应最可机理函数分别为：第Ⅰ阶段，室温~229 ℃，最可机理函数为No.27，g（α）=［-ln（1-α）］3，Avrami-Erofeev方程，随机成核和随后生长机制；第Ⅱ阶段，229 ℃~503 ℃，最可机理函数为No.28，g（[α]）=［-ln（1-[α]）］4，Avrami-Erofeev方程，随机成核和随后生长机制；第Ⅲ阶段，503 ℃~642 ℃，最可机理函数为No.21， g（[α]）=（1-[α]）-1-1，化学反应机制. 2.2.3　求解指前因子[A]　将2.2.2求得最可机函数[g(α)]分别带入公式（4）、公式（5），根据[ln[β/H(x)]~1/T]或[ln[β/(h(x)T2)]~1/T]直线截距B计算指前因子A，见表1所示. 由表1可知，ZrCO3（OH）2?ZrO2非等温热分解各阶段所对应指前因子分别为：第一阶段，lnA=27.407 6~27.864 1；第二阶段，lnA=57.730 4~58.873 6 ；第三阶段，lnA=53.268 8~55.034 6. 3　结　语采用热分析详细研究了碱式碳酸锆在空气中的热分解非等温动力学. ZrCO3（OH）2?ZrO2非等温热分解分为4阶段. 第Ⅰ阶段表观活化能Ea=94.963 kJ?mol-1， lnA=27.407 6~27.864 1，最可机理函数g（[α]）=［-ln（1-[α]）］3，随机成核和随后生长机制；第Ⅱ阶段表观活化能Ea=309.781 kJ?mol-1， lnA=57.730 4~58.873 6，最可机理函数g（[α]）=［-ln（1-[α]）］4，随机成核和随后生长机制；第Ⅲ阶段表观活化能Ea=362.591 kJ?mol-1， lnA=53.268 8~55.034 6，最可机理函数g（α）=（1-α）-1-1，化学反应机制.