气固多相催化反应动力学基础课件.ppt

1、气固多相催化反应动力学马双忱1研究方式：两个分支宏观和微观1、宏观化学动力学：从宏观变量(浓度、温度、压力等)出发，研究反应速率。研究具有普遍意义的规律、特性参数速率常数、级数、活化能、指前因子等。宏观变量cB,T,p宏观动力学宏观动力学特征k,n,Ea,A微观动力学微观变量d,t,r,v,e2、微观动力学：从物质微观特性(如分子尺寸、几何构型、以及分子各运动能级)出发，研究基元反应的速率。化学动力学的任务和目的2反应动力学的重要参量：速率常数：可用来比较催化剂的活性；活化能：判断活性中心的异同；指前因子：用于求取活性中心的数目。气固多相催化反应动力学特点：反应速率与反应物的表面浓度/覆盖度有

2、关；总包反应动力学带有吸附和脱附动力学的特征，甚至内扩散的影响。3催化剂反应物气气相相主主体体4催化反应分七步：反应物从气流主体扩散到固体催化剂的外表面（外扩散过程）反应物进一步向催化剂的微孔内扩散（内扩散过程）反应物在催化剂的表面上被吸附（吸附过程）吸附的反应物转化为反应的生成物（表面反应过程）反应物从催化剂的内表面上脱附下来（脱附过程）脱附下来的生成物分子从微孔内向外扩散到催化剂外表面处（内扩散过程）生成物分子从催化剂外表面扩散到主体气流中被带出（外扩散过程）宏观动力学：宏观动力学：包含传递过程和化学反应的气-固催化反应动力学。本征动力学：本征动力学：不考虑内外扩散传递过程影响的动力学特性

3、(只涉及到上只涉及到上述过程的第述过程的第、V V三步三步)。5多相催化的扩散过程多相催化的扩散过程6 吸附相中的化学反应吸附相中的化学反应+7分解反应物分解反应物HIN2ONH3CH4非催化活化能非催化活化能Ea/kJmol-1183245335335催化剂催化剂AuPtAuPtWOs 锇Pt催化活化能催化活化能Ea/kJmol-110558121136163197230251kS/k(常温常温A不变时不变时)1013102110211019103010241014多相催化反应与气相非催化反应的速率比较实际过程中，Ea降低时，A下降是不可忽略的。指前因子可以代表碰撞频率，催化剂活性位减少可能

4、是导致指前因子减小的原因。8气气-固催化非均相反应动力学固催化非均相反应动力学 气气-固催化固催化非均相非均相反应动力学反应动力学均匀表面吸附均匀表面吸附动力学方程动力学方程 存在控制步骤存在控制步骤 无控制步骤无控制步骤 非均匀表面吸附动力学方程非均匀表面吸附动力学方程 9稳定浓度法（定态法）：稳定浓度法（定态法）：当体系达到稳定状当体系达到稳定状态时，表面中间态物种浓度不随时间变化态时，表面中间态物种浓度不随时间变化平衡浓度法：存在平衡浓度法：存在控制步骤过程的总速率取控制步骤过程的总速率取决于某个控制步骤的速率，而其它步骤处于决于某个控制步骤的速率，而其它步骤处于平衡平衡 1.1 反应速

5、率反应速率1、反应速率定义：反应速率表示反应的快慢，通常定义为参与反应的某种反应物或产物i的量随反应时间的变化率，是反应空间。对于均相催化反应，是反应体系的体积；对于固体催化剂的气固多相催化反应，可以是催化剂的体积V、表面积S或质量W。对于基元反应：有dtdnrii1BArr11BA10如果一个化学反应，反应物分子在碰撞中如果一个化学反应，反应物分子在碰撞中相互作用直接转化为生成物的分子。相互作用直接转化为生成物的分子。2、转换数或转换频率（turnover number or turnover frequency，TOF）：单位时间内每个活性中心引发的总包反应的次数。（总包反应是所包含的所有

6、基元反应的集合，亦称非基元反应或简称总反应亦称非基元反应或简称总反应）1.2 速率方程与动力学参数速率方程与动力学参数微分速率方程：幂式速率方程：k：速率常数；，：反应级数 双曲线式速率方程：k：速率常数；i：常数)P,P,P,P,P(2121 fr BAPkPr 22111PPPkrii11基元过程一般服从Arrhenius定律表观活化能表观活化能：形式上遵循Arrhenius定律的总包反应所对应的活化能。但其具体物理意义视情况而定。动力学参数动力学参数：速率常数、反应级数、指前因子和活化能。1.3 速率控制步骤速率控制步骤催化反应由许多基元反应构成，如果其总速率由其中一步的速率决定，这一步

7、就称为速率控制步骤速率控制步骤。从速控步骤的假定可以进一步推论，在定态时，除速控步骤之外的其它各步都近似地处于平衡状态。E/RT)Aexp(k指前因子指前因子活化能活化能1213内外扩散阻滞效应的判断实例内外扩散阻滞效应的判断实例 In order to be sure that the reaction is working under a chemical regime,preliminary experiments were performed to evaluate the influence of external and internal diffusion.In order to

8、 evaluate the influence of external diffusion,the reactions were run at a constant contact time,at different volumes of catalyst and therefore at different space velocities.It was found that the external diffusion did not play an important role at contact times up to 8 s The influence of internal di

9、ffusion was studied by running the reaction under the same reaction conditions and by changing the particle size of the catalyst.It was found that the internal diffusion did not play an important role at particle size smaller than 0.250 mmApplied Catalysis A:General 192(2000)71791.4 表面质量作用定律表面质量作用定律

10、表面质量作用定律：挪威化学家古德贝格和瓦格（挪威化学家古德贝格和瓦格（Guldberg和和Weage）在前人试验的基础上提出，）在前人试验的基础上提出，理想吸附层中的表面基元反应，其速率与反应物在表面上的浓度成正比，而表面浓度的幂是化学计量方程的计量系数。如对反应：其速率吸附物种表面经常用覆盖度来代替，A是A在催化剂表面的覆盖率，它表示表面反应物A-S的浓度大小，因此也称表面浓度。所以 kaaBAaaBACkCr BAkr 15质量作用定律只适用于基元反应质量作用定律只适用于基元反应 获取速率方程的方法：机理模型法机理模型法和经验模型法经验模型法。机理速率方程机理速率方程：假定一个机理，借助于

11、吸附、脱附、以及表面反应的规律推导出的速率方程。分两种情况：理想吸附模型的速率方程；实际吸附模型的速率方程。2.1 理想吸附模型的速率方程理想吸附模型的速率方程假定吸附、脱附均符合Langmuir模型；表面反应则应用表面质量作用定律。1、表面反应为速控步骤时的速率方程：吸附、脱附步骤处于准平衡。16（1）单分子反应例1、设一反应机理模型为A*A*A*B*k *Ark*代表催化剂，总反应速率由速控步骤表面反应速率表面反应速率决定，由表决定，由表面质量作用面质量作用定律：由Langmuir吸附等温式，*1AAAAA AAAAPPkr1低压或A的吸附很弱时，总包反应表现为一级反应；高压或A的吸附很强

12、时，总包反应表现为零级反应。17实验结果表明PH3在钨表面的分解可能按以上机理进行：例2、设一反应机理模型为*C*C*B*B*C*B*A*A*Ak2*)1(CCBBAAAAAPPPPkkr总反应速率18（2）双分子反应Langmuir-Hinshelwood机理：表面反应发生在两个吸附物种表面反应发生在两个吸附物种间，且此步骤为速控步骤间，且此步骤为速控步骤。（。（Reaction between adsorbates）例1、设反应A+B C按以下机理进行*C*C*C*B*A*B*B*A*Ak2)1(CCBBAABABABAPPPPPkkr总反应速率pAQ QA AfastfastACBpBQ

13、 QB B19例2、设反应A2 B按以下机理进行*2B*2A*2A*2A2k22)1(2121AAAAAPPkkr总反应速率Rideal机理：吸附物种和气相分子间的反应为速控步骤。（Reaction between adsorbate and incoming molecule）例1、设反应A+B C按以下机理进行*C*C*CB*A*A*AkpAQ QA AfastfastRDSkhetAB20CCAABAABAPPPPkPkr1总反应速率金属表面CO氧化反应机理：LH机理See:T.Engel,G.Ertl,Molecular beam studies on the catalytic ox

14、idation of carbon monoxide on a Pd(111)surface,J.Chem.Phys.69(1978):1267-1281.212、吸附或脱附为速控步骤时的速率方程吸附或脱附为速控步骤时的速率方程例、设反应A B的反应机理包括以下三步第一步吸附是速控步骤，其它各步都处于近似平衡，则总反应速率等于第一步的净反应速率假设A相对应的平衡压力为PA*，则A*A*A*B*B*B*-kk AAkPkr0BBAABBAAAAAPPPPP*11*1*022BBAAAAAPPPkPkr*1*总反应平衡常数K=PB/PA*，A是吸附系数，实际上是吸附是吸附系数，实际上是吸附平衡常数

15、，在一定温度下为一常数。平衡常数，在一定温度下为一常数。A=k+/k-最后有 其中BBAPkKPPkr1)(BAKk吸附吸附-化学平衡化学平衡：由吸附平衡和表面化学反应平衡组成的总平衡。233、无速控步骤时的速率方程无速控步骤时的速率方程-稳定态处理法处理模型：在催化反应的连续系列中，如各步骤速率相近和远离平衡的情况下，没有速控步骤。稳态近似：假定各步骤速率相近，从而中间物种浓度在较长时间内恒定，即：其中，i为表面中间物种浓度。例、反应A B的反应机理包括以下两步0dtdi1122A*A*A*B*-kkkk 24根据稳态近似假定，0*dtdA*2*10201AABAkkPkPk因为1*0A21

16、21210212121*kkPkPkkkkkPkPkPkPkBABABAA总反应速率r21212121kkPkPkPkkPkkrBABA这是一种很重要的处理方法，在化学中有广泛运用。25理想吸附模型的机理速率方程通式：n)(（吸附项）推动力项动力学项速率 264、表观活化能和补偿效应以按Rideal机理进行的双分子反应A+B C为例：*C*C*CB*A*A*Ak速控步骤CCAABAAPPPPkr1总反应速率在低覆盖度时，1CCAAPP则BAappBAAPPkPPkr表观速率常数Aappkk27dTddTkddTkdAapplnlnlnk是表面反应速率常数)/exp(RTEAk表面反应活化能A是

17、A的吸附平衡常数)/exp(0RTHA吸附热HEEapp表观活化能表观活化能补偿效应补偿效应：在一系列催化剂上进行某一反应，或在不同条件下处理的同种催化剂上进行一系列反应，将得到的k用Arrhenius方程处理，有时出现活化能E和指前因子A同时增加或同时减小的情况，这样是其中一个的作用抵消了另一个的作用，E和A这种同时同方向的变化称为补偿效应补偿效应。28甲苯在掺碱金属甲苯在掺碱金属K的的V2O5/TiO2催化剂上氧化的补偿效应催化剂上氧化的补偿效应K为电子助剂，降低反应活化能，但K吸附在V单层上，降低了活性中心-表面V物种的浓度，从而降低了指前因子。29 cRTEAdeslnCompensa

18、tion effectA compensation effect can result from a number of sources such as a heterogeneous surface that contains adsorption sites with a range of binding energies（结合能）;lateral interactions（侧向干扰）,in particular if they are strong enough to give rise to coverage dependent phase changes in the adsorbe

19、d layer;Adsorbate-induced changes in the substrate structures（载体结构）.302.2 实际吸附模型的反应速率方程实际吸附模型的反应速率方程假定以反应有如下机理：*DCB*A*A*Ak速控步骤根据表面质量作用定律ABkPr假定吸附热随覆盖度的变化是对数关系利用Freundlich方程ln0 qqnAAKP1nABPkKPr1当吸附或脱附为速控步骤时，可以用Elovich或管孝男的吸附、脱附速率方程代表反应的速率方程。31合成氨反应机理：22223N2*2N*H2*2H*H*N*(NH)*(NH)*H*(NH)*(NH)*NH2*铁催化

20、剂上，认为N2的吸附是合成氨的速控步骤，N2的脱附是氨分解的速控步骤假定吸附能量与覆盖度按线性关系变化，根据Elovich吸附速率方程，反应的速率方程为假设与N2平衡的压力为PN2*，按Temkin等温方程，有22N*(N)*adkk )/exp(22RTPkrNNaa)ln(1*202NNPaf32利用合成氨反应总平衡3223*2HNHNKPPP)ln(1322302HNHNKPPaf合成氨的速率方程)ln(1exp322302HNHNaaKPPafRTPkr令RTf23322230322NHHNNHHNaaPPPkPaKPPkr其中，,0aKkka33同理，氨分解的速率方程13223322

21、30HNHHNHddPPkKPPakr其中，10Kakkd所以合成氨的净反应速率方程当=0.5，由该方程得到的结果与实验吻合。1322323322HNHNHHNdaPPkPPPkrrr34三、经验模型法建立速率方程三、经验模型法建立速率方程直接用某种函数去表达动力学数据，建立速率方程。最常选用的函数是冥函数，对不可逆反应：对可逆反应：参数确定方法：尝试法、孤立法、线性回归法。imiikPrimiimiiiPkPkr35四、动力学方法与反应机理四、动力学方法与反应机理 通常先测定动力学数据，然后用这些数据检验代表不同机理 的速率方程，在检验的基础上，提出反应可能遵循的机理。对对于反应机理的确定，

22、动力学的证明是必要的，但不是充分的于反应机理的确定，动力学的证明是必要的，但不是充分的。4.1 动力学数据的测定动力学数据的测定测定动力学数据应当在内、外扩散不成为速控步骤的情况下测定动力学数据应当在内、外扩散不成为速控步骤的情况下进行。进行。动力学数据测定的主要内容是测定速率。361、连续进料搅拌槽式反应器 反应速率Fm：物料的质量流率；V：反应体积；：停留时间N0，nf:单位质量进料和出料中目的组分的摩尔数。/00fmfnnFVnnr372、柱塞流管式反应器 对于反应体积为V的均匀截面反应管，当反应物料质量流率为Fm（kg/s)时，体积元dV内的物料恒算式为其中反应速率r是单位时间、单位反

23、应体积内反应物转化的摩尔数(mol/m3.s)，x是单位质量进料中目的组分转化的摩尔数(mol/kg)。其积分形式为有时用催化剂质量W代替V。积分反应器：要求高转化率：x-V/Fm图微分反应器：要求转化率低于1%。dxFrdVm)/(mFVddxr 出口入口xxmrdxFV)/(mFVxr384.2 建立速率方程和拟定机理实例建立速率方程和拟定机理实例1、SO2在Pt/Al2O3上的氧化反应根据已有实验结果，假设该反应为Rideal机理，反应步骤为I氧的解离反应II表面反应III产物脱附2O2*2O*吸附系数O223O*SOSO*-kk 表面反应平衡常数KII33SO*SO*脱附系数1/SO3

24、3932SOIIOKK322SOSOO2132SOOSOIIkPkrr总反应总反应平衡常数表面反应为速控步骤应用Langmuir等温方程332233333222221212121212111SOSOOOSOSOSOSOSOOOOOOPPPPPP有33223322212121211)1(SOSOOOSOSOOOPPPKPPkr40利用3221,SOIIOIIKKkkK332233222121212111SOSOOOSOSOOOPPPKPPkr3332121SOSOSOOBPAPKPPrPO2近似为常数，rPKPPBPASOSOOSO33231213SOBPAR或者令可以看出，R和PSO3间有线性

25、关系。41左图说明从动力学数据而言，上述假定机理模型与实验数据符合。从截距和斜率中可以求出，A=-120.6（kPa）1.5hgmol-1B 150.3（kPa）1.5hgmol-1因此SO2氧化速率方程为6.1203.15026.71333221SOSOSOOPPPPr422、2-丁烯氧化脱氢制丁二烯假定速率方程形式为22loglogloglogOBbOaBPbPakrPkPr利用先行回归法求出 a=0.7 b=0.1 c=1.310-3 mmol/sg(kPa)0.8所以可能的经验速率方程形式为1.027.03103.1OBPPr434.3 机理模型研究的几种动力学方法机理模型研究的几种动

26、力学方法 （略讲略讲）基本准则：（1）速率常数和吸附平衡常数为正；（2）速率常数和吸附平衡常数的温度系数合理；（3）速率常数符合Arrhenius定律，吸附平衡常数符合vant Hoff定律。对于速率常数，活化能和指前因子为正；对于吸附平衡，吸附焓和吸附熵通常为负；（4）同系物进行同一反应，其相应平衡常数在相近的温度下有接近的数值。441、压力检定法：从实验上观测速率与反应物压力的不同关系，可以帮助判断反应的可能机理模型。例、对 A+B C，判断其遵循LH机理还是Rideal机理？1）让PA=PB=（1/2）Pt2）让转化率趋于0，速率即为初始速率r0 LH机理：表面反应为速控机理：表面反应为

27、速控步骤步骤LH机理：机理：A的吸附为速控步骤的吸附为速控步骤220)1(ttbPaPr20)1(ttPbPar452、程序升温技术（(TPAT)）：检测样品加热过程中从催化剂表面脱附分子的信号，以脱附分子信号对样品表面温度作图，即得脱附谱。可以推测吸附物种的种数吸附物种的种数和表面覆盖度表面覆盖度，吸附物吸附物种的脱附活化能种的脱附活化能、指前因子指前因子和脱附级数脱附级数，从脱附活化能可以求出吸附活化能吸附活化能。RTEvNdTdNtTTRTEvNdtdNdmaatdmaaexpexp0采用线性升温：脱附速率方程：以温度表示的脱附速率方程：46假设v和Ed与覆盖度无关，当脱附速率最大（脱附

28、峰）时，0)(dTdTdNd pdpapdpdpdadmaddmaamaRTEvNRTERTEvRTEdTdNRTENmvRTERTENvRTENvdTdexp)(2expexpexp0exp2212m=1:m=2:47结论:(1)一级脱附：Tp与表面覆盖度无关；(2)二级脱附：表面覆盖度增加，Tp向低温方向位移；(3)用 ln(Tp2)对 1/Tp作图，可以求出脱附动力学参数。一级反应：二级反应：假设（Na）p=（1/2）N0RERTETdpdplnln1ln202lnln1lnNRERTETdpdp当表面反应为速控步，脱附往往反映了表面反应的动力学信息。当表面反应为速控步，脱附往往反映了表

29、面反应的动力学信息。48TPD：O2/Rh493、过渡应答法：对一个达到稳态的反应体系进行某种扰动（最常见的是浓度），收集该体系向新的稳定态过渡时所显示的各种应答信号。通过对应答信号的分析可以获得有关反应动力学和机理的知识。反应：A+B C*C*C*CB*A*A*A332211-kkkkkk反应机理（1）A A应答：反应物进料中只有A，没有B。当体系达到稳定态后，对A的浓度进行扰动：一是体系立即达到新的稳定态，这表明催化剂对A的吸附已经饱和；二是体系不能立即达到新的稳定态，存在时间过渡态，这表明在过渡时间内A不是发生吸附就是发生脱附。50当反应物入口处浓度由0突跃至CA0：在t到t1内，催化剂

30、上A的吸附增量为图中阴影面积：FV为物料的体积流率。经过时间t1，可继续吸附的量：dtFCCVVtAAA100)(dtFCCVVtAAA1)(0（2）A A应答：反应物进料中有A和B。当体系达到稳定态后，突然升高入口处的A浓度，则体系呈现时间过渡态。出口处A浓度沿曲线II逐渐变到CA0。由于A在步骤II中被消耗，使曲线II位于曲线I的下方。达到稳定态后，A的净吸附速率等于表面反应的净速率。51（3）A C应答：反应物A突跃对产物C应答曲线的影响。如果A发生突跃，产物C在瞬间达到稳定态，则表明在反应途径上任何一种中间物种的吸附或脱附在动力学上都不是重要的，表面反应起速控步骤的作用；在多数情况下，

31、反应物-产物的应答都不是即时的，表现出时间过渡态。如果A、B的浓度同时突跃为0，表面反应没有逆过程：曲线I；表面反应有逆过程：曲线II。52突跃前：0.2atm A+0.8atm He突跃：引入0.2 atm B曲线Q：k1 k3 k2;S组曲线：k1 k3 k2;R组曲线：k1很小。53同位素标记与瞬态应答ssNNTOFNRNROHCHHCO1/3242反应速率与表面物种的平均寿命以及表面物种的数目。可以用相同压力的同位素来迅速反应物，这样可以避免反应压力的突变对参数的影响。5455565758催化反应动力学催化反应动力学：表面吸附-表面反应-表面脱附。由此建立的动力学方程即为催化反应。控制总rr 0dtxd控控制制步步骤骤法法定定态态法法建立方法建立方法Return59 Over!60