高岭土热分解动力学

第 卷第 期 张爱华等：高岭土的热分解动力学 · 1 ·

审稿专家[1]的主要意见:

该文介绍了采用热分析方法研究高岭石热分解过程的动力学特征.研究思路清晰,方法合理,公式应用正确,依据明确,具有一定的可读性.但理论推导的结果均存在一定偏差,如果作者能通过对多个地区高岭石进行比较研究,将更能说明问题. 审稿专家[2]的主要意见:

该文介绍了采用热分析方法研究高岭石热分解过程的动力学特征.研究思路清晰,方法合理,公式应用正确,依据明确,具有一定的可读性.但理论推导的结果均存在一定偏差,如果作者能通过对多个地区高岭石进行比较研究,将更能说明问题.

高岭土热分解动力学

张爱华，何明中，秦芳芳，严慧

（中国地质大学材料科学与化学工程学院，武汉 430074）

摘 要：采用综合热分析仪在动态空气气氛条件下研究了高岭土的热分解过程，利用热重分析数据对高岭土的热分解过程进行了动力学分析。用迭代的等转化率方法获取了准确的活化能，将得到的活化能应用到Malek 方法中推测其反应机理，并进一步求得了指前因子A 。结果发现：在400~700 ℃高岭土脱去羟基，生成了结晶度较差的偏高岭土，该过程遵循的化学反应模型，其微分和积分表达式分别为n

f )1()(αα-=，n

G n

---=-1)1(1)(1αα其

中n =2.1，表观活化能为182.428 kJ/mol ，指前因子A 的范围为：(4.566~4.635) ×1011 s -1。

关键词：高岭土；热分解；动力学；迭代法；Malek 法

中图分类号：TQ170.1 文献标识码：A 文章编号：0454-5648(2009)12

K INETICS FOR T HERMAL D ECOMPOSITION OF K AOLINITE

ZHANG Aihua, HE Mingz fang, YAN Hui

(Faculty of Material Science and Chemical Engineering, China University of Geosciences,

Wuhan 430074, China)

Abstract: The thermal decomposition processes of kaolinite have been studied in dynamic air using simultaneous thermal analysis instrument. Thermogravimetricdata was used to carry out the kinetic analysis. Firstly, the exact activation energy was calculated by the interative method; Secondly, using the activation energy in Malek method to confer the reaction mechanism; Lastly, using the above data to get the pre-exponential factor A . The hydroxylation of kaolinite in the temperature range 400–700 ℃ followed the model of chemical reaction. The function were ()()n f αα-=1，()n

G n

---=

-111)(1ααwhere n =2.1, the apparent activation energy was

182.428 kJ/mol, the range of the A was from 4.566×1011 to 4.635×1011 s -1.

Key words: kaolinite; thermal decomposition; kinetic; iterative method; Malek method

高岭土是一种天然矿物，在我国有着丰富的储量,现已成为造纸、陶瓷、橡胶、化工、涂料、医药和国防等几十个行业所必需的矿物原料。油漆涂料和造纸是国内优质煅烧高岭土的最主要的消费领域。[1] 因此，对高岭土的热分解过程进行动力学分析是十分重要和迫切的。通过动力学分析可以得到高岭土热分解的活化能，指前因子和反应机理函数，从而为生产实践提供指导。但是目前国内对高岭土的热分解过程进行动力学分析的很少，同时在

反应级数方面与国外报道的也有所不同。[2-5]

实验利用迭代的等转化率方法和Malek 方法研究了广东茂名产高岭土的热分解过程，计算出了热分解过程的动力学参数。

1 动力学数据的处理方法

1.1 迭代的等转化率法求取活化能E a

根据等温动力学理论，固体分解反应的动力学 方程可表示为：

第 卷 第 期2009年月

硅 酸 盐 学 报

JOURNAL OF THE CHINESE CERAMIC SOCIETY

V ol. ，No. , 2009

· 2 · 硅 酸 盐 学 报

)(exp d d a ααf RT E A t ⋅⎪⎭

⎫

⎝⎛-=

(1) 其中α为转化率，)(αf 为动力学模型的微分形式，A 为指前因子，E a 为活化能，T 为反应温度，R 为普适气体常数，t 为时间。

若在非等温条件下，升温速率β恒定，t T d /d =β，则

)(exp d d a ααf RT E βA T ⋅⎪⎭

⎫

⎝⎛-=

(2) 分离变量并进行积分变形可以得到Ozawa 方程和KAS 方程。

Ozawa 方程为：

RT E

R G AE βa a 0516.1)(0048.0ln ln -⎥⎦

⎤⎢⎣⎡=α (3)

收稿日期： 修改稿收到日期： 第一作者：张爱华(1983~)，女，硕士研究生 通讯作者：何明中(1954~)，男，教授。

Received date: Approved date :

First author:ZHANG Aihua(1983~),female,graduate student for master degree.

E-mail ：zhangaihua525@.

Correspondent author: HE Mingzhong(1954~),male,professor. E-mail: hmzhong@

KAS 方程为：

RT E R G AE T β

2a a )(ln ln -

⎥⎦

⎤⎢⎣⎡=α (4)

其中，)(αG 为动力学模型的积分形式，在上述的

计算中引入了温度积分近似公式，对计算结果产生了一定的误差。

迭代法[6]即让所求近似的E a 值逐步逼近真实的E a 值，其方程为：

RT E R G AE x H β

a a 0516.1)(0048.0ln )(ln -⎥⎦

⎤⎢⎣⎡=α (5) RT E R G AE T x h β

a a 2)(ln )(ln

-

⎥⎦

⎤⎢⎣⎡=α (6)

H (x )和h (x )随着x [)/(a RT E x =]的变化而变化，它

们分别为方程(5)和(6)进行迭代的定义函数，其定义式为：

)0516.1exp(0048.0/)()exp()(2

x x x h x x H --= (7)

12024012020968818)(234

234+++++++=x x x x x

x x x x h (8) 迭代过程可分为3步：(1) 设H (x )=1或h (x )=1并估算初始的活化能E a ，即为常规的等转化率法所求得的E a1；(2) 由E a1值以及不同升温速率下所对应的T ，代入x 以及H (x )或h (x )的定义式中，求取x 及相应的H (x )或h (x )的值，再将H (x )或h (x )的值代入方程(5)或(6)中，以)](/ln[x H β对1/T 或

2

)(ln T

x h β

对1/T 进行线性回归，通过斜率可求得新的E a2；(3) 以E a2代替E a1，重复(2)，直至E a i -E a(i -1)<0.01 kJ/mol 时，即可认为所求得的E a 值为真实值。[7]

在不同的α处采取相同的步骤得到α与E a 的关系，如果在主要的α范围内（0.2≤α≤0.8），[8]E a 的数值基本不变，则认为该反应可以用单一的动力

学模型来描述。

1.2 Malek 法判定最可能的机理函数

由动力学方程(2)可知： T RT E A f d exp )(d a ⋅⎪⎭⎫

⎝⎛-=βαα (9)

⎰

⎰⎰⋅⎪⎭

⎫ ⎝⎛-≈⋅⎪⎭

⎫

⎝⎛-==T T T T

RT E A T

RT E A G f 00d exp d exp )()(d a a 0ββαααα (10)

因此

)()(a

x P R

AE G βα= (11)

)(x P 称为温度积分近似式，其表示形式如下：

x x x P x

x

d e )(2

⎰

∞--=

(12) 其中)/(a RT E x =，由于)(x P 是不收敛积分，得不到精确解析解，采用唐万军提出的近似公式来代替：[9]

87391198

.100198882.11

)exp()(+⋅

-=x x x x P 当α=0.5时

)()5.0(5.0a

x P R

AE G β= (13)

式(11)除以式(13)可得，

)

5.0()

()()(5.0G G x P x P α= (14)

等号左边与实验数据有关，将实验数据代入左端。而等号右边与定义的理论机理函数有关，将α（84.016.0≤≤α，间隔0.02）和45种机理[10]对应的机理函数代入右端。然后计算实验值和理论值之间的标准偏差，根据标准偏差的大小可以找到合适的机理函数。