创刊于1987年, 双月刊
主管:

江西理工大学

主办:

江西理工大学
江西省有色金属学会

ISSN:1674-9669
CN:36-1311/TF
CODEN YJKYA9

应用正交设计研究火焰原子吸收分析中的干扰效应

邱春

邱春. 应用正交设计研究火焰原子吸收分析中的干扰效应[J]. 有色金属科学与工程, 2000, 14(3): 33-35.
引用本文: 邱春. 应用正交设计研究火焰原子吸收分析中的干扰效应[J]. 有色金属科学与工程, 2000, 14(3): 33-35.
QIU Chun. Computational eliminating the interference effects of flame atomic absorption analysis by orthogonal design[J]. Nonferrous Metals Science and Engineering, 2000, 14(3): 33-35.
Citation: QIU Chun. Computational eliminating the interference effects of flame atomic absorption analysis by orthogonal design[J]. Nonferrous Metals Science and Engineering, 2000, 14(3): 33-35.

应用正交设计研究火焰原子吸收分析中的干扰效应

详细信息
    作者简介:

    邱春(1962-), 男, 江西泰和人, 赣州有色冶金研究所分析实验室工程师, 主要从事分析仪器与分析方法研究

  • 中图分类号: O433.5+1

Computational eliminating the interference effects of flame atomic absorption analysis by orthogonal design

  • 摘要: 采用正交多项式回归设计方法, 试验和研究了稀土氧化物中钙的火焰原子吸收光谱法(FAAS)测定时的干扰效应。结果表明:该方法可以有效地研究和消除干扰效应。用回归方程修正的结果和实际浓度一致。
    Abstract: The interference effects for FAAS determining calcium in RE oxides were studied by using orthogonal polynomial regression design.In this paper, it is shown that the method is efficient to study and eliminate interference effects in FAAS and the results corrected by using polynomial regressive equation are consistent with the actual concentrations.
  • 氧化镝用作制取金属镝的原料、玻璃添加剂,还用于金属卤素灯、磁光记忆材料、钇铁或钇铝石榴石、原子能工业中.因此,氧化镝的应用范围十分广泛,应用价值也很高.制取氧化镝通常有2种方法,一种是草酸镝盐热分解制得稀土氧化物,另一种是碳酸镝盐热分解制得稀土氧化物.草酸镝盐是目前制备氧化镝,尤其是制备具有可控粒度的氧化镝常用的前驱体,它具有沉淀物晶型好,易于过滤分解等优点.因此,研究草酸镝盐的热分解动力学,对掌握草酸镝盐热分解过程的规律以及控制产物的组织形貌和粒度有着重大的意义.

    热分析技术被广泛用于热分解机理研究和动力学的研究中[1-5],许多学者用此法对草酸盐的热分解过程进行了研究,如李靖华等[6]研究了草酸锰的热分解,潘云祥等[7]研究了草酸镍脱水过程的非等温动力学,邹华红等[8]研究了不同样品量、升温速度、粒径大小和反应气氛下水合草酸钙热分解差热分析数据的区别,分析这些因素对热分析数据的影响.关于稀土草酸盐水合物的热分解也有过较多研究[9-10].如倪兆艾等[11]采用程序式升温的热分析法研究稀土草酸盐的热分解机理,得出脱水峰值温度随镧系元素原子序数的增大而升高.郑晓玲等[12]用热分析技术研究了重稀土草酸盐在氢气气氛中的热分解规律,并提出了热分解机理.张克立等[13]用热分析法研究了十水草酸镧的热分解机理.顾翼东等[14]研究了Pr2(C2O43﹒10H2O热分解过程,得出其DTA曲线有3个峰.方正东等[15]用热重-差热分析法研究了十水草酸铈在程序升温下的热分解行为,TG曲线表明它有5个分解阶段,前3个阶段为脱水过程,后2个阶段为无水草酸铈热分解过程,最终产物为CeO2.他们都只是研究了热分解机理,而未研究热分解动力学.关于稀土草酸盐的动力学研究,胡艳宏等[16]利用DTA法研究了草酸铈前驱体的热分解动力学,詹光等[17]研究了十水草酸镧在室温至900 ℃下的热分解动力学,张湘辉等[18]研究了La2(CO33﹒3.4H2O在空气中热脱水的非等温动力学.在国外研究中,AYODELE B V等[19]采用热重分析法研究了轻稀土硝酸盐水合物热分解的非等温动力学及机理,分析其分解机理,得出了硝酸盐前驱体热分解的各动力学参数,KOVÁCS T N等[20]研究了2种不同的钼酸铵在惰性气体和氧气气氛下的热分解行为.但关于稀土草酸盐的许多研究都只是研究了分解机理,对动力学参数的研究较少.

    为了研究草酸镝的热分解过程机理和动力学,拟用热重-差热分析法研究草酸镝在程序升温下的热分解行为,通过不同升温速率下草酸镝盐分解的TG-DTA数据对其进行动力学分析,以确定草酸镝盐分解的分解反应级数、分解活化能和分解机理,为草酸镝盐分解制备良好形貌和粒度的氧化镝提供一定的理论依据.

    称取5 g氧化镝Dy2O3加水溶解,将12 mol/L的HCl配制成6 mol/L的HCl溶液,再称取58.006 g固体H2C2O4﹒2H2O并加水溶解配制成20 %的草酸溶液,取该溶液70 mL分别加入氧化镝Dy2O3水溶液中,置于磁力子搅拌器中,在60 ℃下持续搅拌30 min.搅拌后的溶液趁热过滤,将过滤后的滤渣置于表面皿中,放在烘干室内,100 ℃下烘干3 h.

    取适量的待分析样品置于北京博渊同步热分析仪DTU-3B的坩埚,以5 ℃/min、10 ℃/min、15 ℃/min、20 ℃/min、25 ℃/min、30 ℃/min的升温速度加热样品,初始温度设定30 ℃,终止温度设定1 000 ℃.通过同步热分析仪记录热分析数据,用Origin作图软件作出不同的升温速度下的TG曲线和DTA曲线,根据不同数据处理法算出动力学参数.

    Kissinger法[21]用到以下方程:

    $$\ln \left( {\frac{\beta }{{{\mathit{\Gamma }^2}}}} \right) = \ln \left( {\frac{{RA}}{{Ea}}} \right) - \frac{{Ea}}{R} \cdot \frac{1}{\mathit{\Gamma }}$$ (1)

    进行差热分析实验时,DTA会得到峰值温度T,可以计算$\ln \left( {\frac{\beta }{{{\mathit{\Gamma }^2}}}} \right)$和$\frac{1}{\mathit{\Gamma }}$,在每个升温速度下,将$\frac{1}{\mathit{\Gamma }}$作为横坐标,$\ln \left( {\frac{\beta }{{{\mathit{\Gamma }^2}}}} \right)$作为纵坐标,可得到1条直线,斜率为$ - \frac{{Ea}}{R}$,截距为$\ln \left( {\frac{{RA}}{{Ea}}} \right)$,此时,便可以很方便地求得化学反应表观活化能Ea和指前因子A.

    Ozawa法[22]由于不知道机理函数,因此只能求出化学反应的表观活化能Ea. Ozawa法的基本方程:

    $$\lg \beta = \lg \left[ {\frac{{A \cdot Ea}}{{R \cdot G\left( a \right)}}} \right] - 2.315 - 0.4567\frac{{Ea}}{{RT}}$$ (2)

    如式(2),将不同升温速度下的TG曲线拟合在同一坐标系中,做1条转化率α相同的线,与每条TG线都会有1个交点,读出每个交点的温度T.再将各点的lgβ和$\frac{1}{\mathit{\Gamma }}$作图,可得到1条直线,该直线的斜率为-0.456 7$\frac{{Ea}}{R}$,采用上述相同转化率法多计算几次取平均值,就可以最终确定化学反应表观活化能Ea.

    Crane法[23]的动力学方程为:

    $$\frac{{{\rm{d}}\left( {\ln \beta } \right)}}{{{\rm{d}}\frac{1}{{{T_{\max }}}}}} = - \left( {\frac{{Ea}}{{n \cdot R}} + {T_{\max }}} \right)$$ (3)

    找出DTA曲线中在不同升温速度β下,每个峰的对应的峰值温度Tmax,求出lnβ和$\frac{1}{{{T_{\max }}}}$,将lnβ对$\frac{1}{{{T_{\max }}}}$作图,通过最小二乘法拟合,拟合直线的斜率为$ - \frac{{Ea}}{{n \cdot R}}$,把Ozawa法求得的表观活化能Ea值代入求出反应级数n.

    Coats-Redfern法[24]是根据式(4):

    $$\ln \left[ {\frac{{f\left( a \right)}}{{{\mathit{\Gamma }^2}}}} \right] = \ln \left( {\frac{{A \cdot R}}{{\beta \cdot Ea}}} \right) - \frac{{Ea}}{{R \cdot T}}$$ (4)

    把常见的机理函数代入式(4)中,将$\ln \left[ {\frac{{f\left( a \right)}}{{{\mathit{\Gamma }^2}}}} \right]$对$\frac{1}{T}$作图,并通过最小二乘法进行线性拟合,该拟合直线的斜率为$ - \frac{{Ea}}{R}$,截距为ln($\frac{{A \cdot R}}{{\beta \cdot Ea}}$),就可以求得每种反应机理函数下的表观活化能Ea和指前因子A.选取拟合直线相关系数r最好的一组,该组表观活化能Ea和指前因子A即为所求动力学参数,其代入的机理函数则为该反应的机理函数.

    为了研究六水草酸镝热分解的动力学参数,先对其分解机理进行分析,确定分解机理后,采用Kissinger法、Ozawa法和Crane法进一步确定热分解反应的表观活化能Ea、指前因子A和反应级数n,再根据Coats-Redfern法,确定其热分解反应的机理函数gα).

    为了分析草酸镝的热分解机理,需要对草酸镝TG-DTA曲线进行分析,研究DTA曲线上的吸热峰和放热峰,结合TG失重量来确定分解机理.现选取20 K/min的TG-DTA-DTG曲线进行分析如图 1所示.

    图  1  20 K/min六水草酸镝热分解TG-DTA曲线
    Figure  1.  TG-DTA curves of dysprosium oxalate with six water thermal decomposition reaction in 20 K/min

    图 1可见,TG失重线共有2个失重台阶,因此,草酸镝的热分解过程为2步分解过程.第1个失重过程的温度在450~500 K之间,在DTA曲线上是一个吸热峰,故认为它是一个失水的过程,根据TG线算出第1个台阶的失重率约为16.16 %,与失去6个水分子的理论失重率15.77 %相当,仅相差0.39 %.因此,草酸镝在这个过程中应是失去6个结晶水.对应第2个平台,DTA曲线上有一个放热峰,根据TG线得到该过程的失重率为30.11 %,与失去3个CO2分子和3个CO分子的理论失重率31.53 %相近,仅相差1.42 %. DTA最后一个放热峰应是产物的晶型转变过程,因为这时在DTG线上没有出现峰,说明此时没有发生质量变化.且样品草酸镝在505 ℃(778.16 K)和800 ℃(1 073.15 K)焙烧1 h后的XRD物相分析均为Dy2O3.(图 2图 3所示),只是对应不同的晶型结构而已,这充分说明DTA上的第3个放热峰确实为相变过程所致.样品DTA峰值温度和TG失重率如表 1所列.通过查阅相关文献[25],重稀土草酸盐结晶水个数一般为6个,结合实验TG失重率,可以确定本实验草酸镝含6个结晶水,为六水草酸镝.

    图  2  六水草酸镝热分解505 ℃产物XRD物相分析
    Figure  2.  XRD pattern of the product in 505 ℃ of dysprosium oxalate with six water thermal decomposition reaction
    图  3  六水草酸镝热分解最终产物XRD物相分析
    Figure  3.  XRD pattern of the final product of dysprosium oxalate with six water thermal decomposition reaction
    表  1  六水草酸镝热分解机理数据
    Table  1.  Data of thermal decomposition mechanism of dysprosium oxalate with six water
    300
    /(K·min-1)
    DTA峰值
    温度/K
    TG失
    重率/%
    理论
    失重率/%
    相差
    /%
    第1个平台480.743 416.1615.770.39
    第2个平台684.255 830.1131.531.42
    下载: 导出CSV 
    | 显示表格

    根据以上分析,可以确定六水草酸镝(Dy2(C2O43·6H2O)热分解过程如下:

    第1步:(Dy2(C2O43·6H2O→(Dy2(C2O43+6H2O↑

    第2步:(Dy2(C2O43→Dy2O3+3CO2↑+3CO↑

    为了研究草酸镝水合物在分解过程中组织形貌的演变过程,对草酸镝原样、300 ℃和800 ℃下的灼烧样进行了SEM检测,结果如图 4所示.

    图  4  草酸镝的SEM形貌
    Figure  4.  SEM images of dysprosium oxalate

    图 4可见,原样及2个灼烧温度下的试样,在相同的放大倍数下,都呈片层状,说明分解过程样品的形貌具有继承性.并且,随着温度的升高,试样的粒度明显减小,片层状更加明显.

    为分析草酸镝的分解过程动力学,现将不同升温速度下的DTA曲线拟合在一起,得到图 5.

    图  5  不同升温速度下六水草酸镝热分解DTA曲线
    Figure  5.  DTA curves of dysprosium oxalate with six water thermal decomposition reaction at different heating rates

    图 5可见在草酸镝的热分解过程中,DTA曲线中一共有3个峰,其中第1个峰是吸热峰,后两个峰是放热峰,但是根据前文反应机理的推断,最后一个峰为反应过程中发生的相变等过程,因此,只需分析前两个峰.

    六水草酸镝热分解2步反应的DTA数据经过Kissinger法处理后,求得各个反应过程的动力学参数如表 2所列.

    表  2  六水草酸镝热分解过程分步反应动力学参数
    Table  2.  Kinetic parameters of dysprosium oxalate with six water thermal decomposition reaction
    第1步反应第2反应
    斜率-8 679.927 26-15 087.927 1
    截距8.552 2412.085 01
    指前因子A4.49×1072.67×109
    活化能Ea/(kJ·mol-1)72.16125.44
    下载: 导出CSV 
    | 显示表格

    经Kissinger法处理差热分析,六水草酸镝((Dy2(C2O43﹒6H2O)热分解过程中发生的第1步脱水反应表观活化能Ea和指前因子A分别是72.16 kJ/mol和4.49×107,第2步分解反应表观活化能Ea和指前因子A分别是125.44 kJ/mol和2.67×109.

    将不同升温速度下草酸镝热分解TG线拟合在同一个坐标系中,得到图 6,利用等转化率法对其进行分析.这6个升温速度下,草酸镝热分解过程的TG线重合性都较好,TG线上共有2个失重平台.

    图  6  不同升温速度下六水草酸镝热分解TG曲线
    Figure  6.  TG curves of dysprosium oxalate with six water thermal decomposition at different heating rates

    根据Ozawa法,求得2个分解反应的活化能值如表 3所列.

    表  3  六水草酸镝热分解反应不同转化率活化能值
    Table  3.  Activation energy in different conversion rate of dysprosium oxalate with six water thermal decomposition reaction
    α第1步反应
    Ea / (kJ·mol-1)
    第2步反应
    Ea / (kJ·mol-1)
    0.145.08113.64
    0.275.37118.86
    0.378.66118.76
    0.477.94118.89
    0.576.37119.58
    0.673.01113.66
    0.772.06104.90
    0.876.5898.30
    0.978.97110.24
    Ea平均值
    / (kJ·mol-1)
    72.67112.98
    下载: 导出CSV 
    | 显示表格

    实验采用Ozawa法处理差热分析数据,六水草酸镝((Dy2(C2O43﹒6H2O)热分解过程中发生的第1步脱水反应表观活化能Ea是72.67 kJ/mol,第2步分解反应表观活化能Ea是112.98 kJ/mol.

    为求分解反应的反应级数,采用Crane法进行求解.Crane法是用和Kissinger法相同的峰顶温度,把lnβ对每个峰的峰顶温度$\frac{1}{{{T_{\max }}}}$作图.将2个峰的lnβ对$\frac{1}{{{T_{\max }}}}$图分别作出,通过最小二乘法分别拟合,并拟合在图 7中,根据Crane法,需要代入其他方法求得的活化能,才能求出反应级数,因此,将Ozawa法求得的活化能代入,每条拟合直线的斜率和Crane法求出的反应级数如表 4所列.

    图  7  Crane法拟合曲线
    Figure  7.  The graph fitted by crane method
    表  4  Crane法计算反应级数
    Table  4.  Reaction order calculated by the Crane method
    参数第1步反应第2步反应
    斜率-9 627.423 40-16 425.345 05
    Ozawa求得的活化能Ea /(J • mol-1)72 670.776 36112 979.814 6
    Crane法求得的反应级数n0.910.83
    下载: 导出CSV 
    | 显示表格

    因此,采用Crane法处理六水草酸镝热分析数据,得出六水草酸镝第1步脱水反应的反应级数n为0.91,近似为1级反应,第2步分解反应的反应级数n为0.83,也近似为1级反应.

    根据Coats-Redfern法,在确定机理函数时,把12种常见的机理函数代入Coats-Redfern方程中,求出每个反应机理函数下的$\ln \left[ {\frac{{g\left( a \right)}}{{{\mathit{\Gamma }^2}}}} \right]$.

    将$\ln \left[ {\frac{{g\left( a \right)}}{{{\mathit{\Gamma }^2}}}} \right]$对$\frac{1}{\mathit{\Gamma }}$作图,并通过最小二乘法进行线性拟合,得到2步反应的相关系数如表 5所列.

    表  5  六水草酸镝热分解反应数据
    Table  5.  Data of dysprosium oxalate with six water thermal decomposition reaction
    机理函数g(a)表观活化能
    Ea /(kJ·mol-1)
    指前因子A第1步反应
    相关系数r
    表观活化能
    Ea / (kJ·mol-1)
    指前因子A第2步反应
    相关系数r
    A227.141.70×1020.893 3927.571.03 ×1010.581 06
    A315.365.75×1000.854 1514.266.16×10-10.434 54
    A444.811.93×1040.911 947.554.46×1020.654 84
    D190.437.29×1080.815 4294.545.84×1050.548 =1
    D2101.958.83×1090.854 23108.524.40×1060.603 59
    D3116.911.17×10110.896 51127.083.36×1070.669 16
    D4106.867.52×1090.869 72114.603.12×1060.626 84
    F162.481.84×1060.919 5367.531.60×1040.685 80
    F293.078.37×1090.985 59106.422.79×1070.842 29
    F3131.342.59×10140.998 96155.522.74×10110.921 38
    R250.683.32×1040.858 4152.794.29×1020.580 16
    R354.356.28×1040.880 7557.35=.13×1020.61668
    下载: 导出CSV 
    | 显示表格

    结合其他方法求得的活化能数据和指前因子数据,再根据求出的相关系数,找出最匹配的机理函数即为Coats-Redfern法求得的机理函数.Kissinger法求得第1步反应的表观活化能为72.16 kJ/mol,指前因子为4.49×107,与之最匹配的机理函数为F1;第2步反应的表观活化能为125.44 kJ/mol,指前因子为2.67×109,与之最匹配的机理函数为F2.所以,六水草酸镝热分解第1步反应机理函数为F1,说明反应的控制步骤为化学反应,该反应表观活化能为62.48 kJ/mol,指前因子为1.84×106;第2步反应的机理函数为F2,说明反应的控制步骤为化学反应,该反应表观活化能为106.42 kJ/mol,指前因子为2.79×107.故六水草酸镝热分解2步反应的反应机理函数分别为:

    第1步:ga)=-ln(1-a

    第2步:ga)=(1-a-1-1

    其中,a为反应某时刻转化率.

    文中用3种方法计算了草酸镝热分解活化能,结果汇总于表 6.

    表  6  六水草酸镝热分解反应表观活化能值
    Table  6.  Apparent activation energy of dysprosium oxalate with six water thermal decomposition reaction
    样品
    名称
    反应
    步骤
    Kissinger法
    Ea /(kJ·mol-1)
    Ozawa法
    Ea /(kJ·mol-1)
    Coats-Redfern法
    Ea /(kJ·mol-1)
    草酸镝第1步72.1612.6162.48
    第2步125.44112.98106.42
    下载: 导出CSV 
    | 显示表格

    根据表 6可以看出3种方法计算的草酸镝热分解过程活化能值相差不大,说明计算较为合理.

    在空气中六水草酸镝((Dy2(C2O43·6H2O)的热分解发生2步反应,分解过程为:

    第1步:(Dy2(C2O43·6H2O→(Dy2(C2O43+6H2O↑

    第2步:(Dy2(C2O43→Dy2O3+3CO2↑+3CO↑

    第1步反应机理函数为F1,表观活化能为62.48 kJ/mol,指前因子为1.84×106;第2步反应的机理函数为F2,该反应表观活化能为106.42 kJ/mol,指前因子为2.79×107. Crane法计算得到的六水草酸镝热分解2步反应的反应级数均近似为1.

    六水草酸镝热分解2步反应的反应机理函数分别为:

    第1步为F1:g(a)=-ln(1-a

    第2步为F2:g(a)=(1-a-1-1

    其中,a为反应某时刻转化率.

  • 表  1   试验安排与结果

    下载: 导出CSV

    表  2   试验结果的方差分析

    下载: 导出CSV

    表  3   回归系数及其变差平方和

    下载: 导出CSV

    表  4   回归系数的方差分析

    下载: 导出CSV

    表  5   R计算值与实测值以及结果误差比较

    下载: 导出CSV
  • [1] 金曼一, 张会强, 周欣磊.火焰原子吸收光谱分析中干扰效应计算消除方法的验证〔J〕.分析试验室, 1990, 9(3):38-39. http://www.oalib.com/paper/4558623
    [2]

    L Pszonicki, A Lechotycki, M Krupinski.Computational Elimination of Interferent Effects in Flame Atomic-Absorption Spectrometry〔J〕.Talanta, 1988, 35(6):465-472. doi: 10.1016/0039-9140(88)80109-7

    [3]

    Michael Thompson, Stephen J Walton, Shirley J Wood.Statistical Appraisal of Interference Effects by Atomic-absorption Spectrophotometry in Applied Geochemistry〔J〕.Analyst, 1979, 104(4):299-312.

    [4] 邓勃, 秦建候, 王小芹.分析测试中试验设计和优化方法二、正交多项式回归设计〔J〕.分析试验室, 1985, 4(11):47-56. https://www.wenkuxiazai.com/doc/a9ef88213b3567ec112d8a27...
    [5] 邓勃, 王小芹.利用正交多项式回归设计研究石墨炉原子吸收法测定铂的条件初探〔J〕.分析试验室, 1986, 5(4):11-14. http://www.cnki.com.cn/Article/CJFDTotal-HBGT201504029.htm
    [6] 李述信.原子吸收光谱分析中干扰及消除方法〔M〕.北京:北京大学出版社, 1987.1-23.
表(5)
计量
  • 文章访问数:  69
  • HTML全文浏览量:  15
  • PDF下载量:  2
  • 被引次数: 0
出版历程
  • 收稿日期:  2000-05-28
  • 发布日期:  2000-09-29
  • 刊出日期:  2000-09-19

目录

/

返回文章
返回