多相催化动力学1-5章 北京化工大学

1、1,多相催化反应动力学理论与实践,学时：40 讲解： 工业催化与反应器研究室 综合楼517 电话64436787,2,化学反应动力学建立与应用的引入, 过程工程科学的核心理论 三传一反（热量、质量、动量传递，反应）, 过程工程问题解决的基础条件 方向、限度（物理化学、热力学） 速率（热质传递、化学反应等）, 化学反应动力学通式（过程源、汇项的解决）,3, 化学反应动力学如何获得？,对于化学反应： a+b c, 基元反应 幂律型方程,a、b和c竞争吸附双曲型方程, 非基元反应 幂律型方程,a、b和c竞争吸附双曲型方程, 动力学参数k0、e、na、nb、kj如何确定？,4,涉及/解决的问题, 多相

2、催化的含义 动力学的概念，化学反应动力学内涵 化学反应动力学建立过程、实施方案 反应过程热力学分析（平衡、网络结构） 催化剂性能考察（本体性能、附加性能） 反应过程动力学特性分析 动力学实验设计、实施 动力学模型选取、建立及其模型参数估计 化学反应动力学模型适定性、可信性检验,5, 课程目的 掌握建立工程实用的多相催化反应动力学方程（模型）的原理和方法（建模）（动力学模型建立、数据获取和模型识别等）。, 催化剂制备与应用 粉体制备成型（催化剂工程设计）装填反应器 反应器分析与设计, 多相催化的含义 均相催化、多相催化（催化剂、反应物系与催化剂或反应物系间非同一物相。气固、气液、液固、气液固）。

3、,6, 课程意义 为催化剂性能改进、工艺过程开发和反应器优化设计与操作提供重要的基础应用工具。, 课程特点 系统性、理论与实践的结合性和工程实用性。实用性（技巧性）很强（以工业应用为目的），对数学、基本催化作用原理和计算机编程均有较高的要求。, 课程学习方法 超大型作业和倒叙的方式如何建立实用的动力学方程？,7, 课程安排（40学时） chap 1 绪论（4小时） chap 2 多相催化反应动力学基础理论（8小时） chap 3 多相催化反应动力学实验研究（10小时） chap 4 数据处理和动力学模型识别（6小时） chap 5 非定态反应动力学（3小时） chap 6 多相催化宏观反应动力

4、学及其描述（3小时） chap 7 多相催化反应动力学研究的实例分析（4小时） 结课测验 （2小时）, 课程考核方式 考试 70%，课程论文 30%,8, 参考书 1 berty j m. experiments in catalytic reaction engineering. elesiver, 1999 2 smith j m. chemical engineering kinetics （3th）, mcgraw-hill book company, 1981 3 张继炎等译，贝伦斯 m，霍夫曼 h，林肯 a著，化学反应工程学，中国石化出版社. 1994 4 李绍芬. 化学与催化反应

5、工程，化学工业出版社. 1994 5 朱炳辰. 化学反应工程学，化学工业出版社. 1998 6 陈甘堂. 化学反应工程学，化学工业出版社. 1981 7 李成岳译. 马特路斯著. 催化反应器中的非定态过程，科学出版社. 1994,9,1 绪论,1.1 化学反应动力学的概念 动力学含义：研究在一定推动力作用下，某个或数个具有特定物理含义的变量与其影响因素之间的相互关系。 传递动力学（石头击水水波纹传递过程描述） 气体动力学（飞机风洞实验 ） 流体流动动力学 摩擦动力学 化学反应动力学等,10, 化学反应动力学 定义：定量描述化学反应速率及其影响因素之间的函数关系。实际为有可能定量地描述化学反应随

6、时间的变化对过程条件的依赖关系。, 影响因素：温度、压力、组分浓度、催化剂种类与浓度、传递过程（流动、扩散等） 。, 化学反应速率的表达：,11,反应程度的概念：,多相催化中常见的化学反应速率表达式：,误区说明：化学反应速度（矢量）和化学反应速率 （标量）？, 化学反应动力学的分类：本征动力学（微观动力学、化学动力学）、宏观动力学 约定：只要未特别说明，均指本征动力学,12,1.2 多相催化动力学的研究内容,1.2.1 化学反应速率定量关系的建立 考察多相催化反应的速率变化规律，建立多相催化反应速率与反应条件（c、p、t等）之间的经验（半经验或理论关系）关系，给出定量的函数关系式。,1.2.2

7、 化学反应图式或反应机理的印证 动力学对机理的提出具有否定性判据效果，但无肯定性效果 表观动力学方程 工程动力学,13,1.2.3 化学反应动力学特征参数与热力学性质 和物质结构特征的关系 利用动力学的研究结果和物质结构理论，分析化合物结构与其反应能力之间的关系；建立多相催化反应体系的动力学特征参数与其热力学性质和物质结构特征之间的经验或半经验关系。,1.2.4 基元反应速率常数的理论计算 从理论上（即从量子化学及统计物理）计算基元反应速率常数，并对复杂反应给以定量描述 。,1.2.11.2.3已进行，1.2.4内容仅处于初级阶段,14,1.3 常见的动力学型式,1.3.1 幂律型,1.3.2

8、 双曲型,15,1.3.3 分段型,16,1.3.4 机理型 朗格谬尔（langmuir）吸附机理模型 焦姆金（temkn -ea、eda ） 弗朗得里希（freundlich-ea、edlna）,由于缺乏对催化过程的足够认识，多数动力学方程都是半经验型的（如级数ni=0.784），是对实验数据的合理拟合。因此，外推能力很差，且同种催化剂、同种反应过程，因不同人、不同方式和不同时域获得的动力学方程会有相当大的差异。使用时一定要注意动力学方程的测试条件。现今多相催化动力学是一个实验性特强的工作！,17,1.4 化学反应动力学的应用, 反应条件选取：选择性（平行、连串） 催化剂改进方向 催化剂设计

9、 催化剂工程设计 反应过程分析与优化,18,1.5 非均相催化反应动力学涉及学科, 催化剂制备与应用 原料预处理原粉制备焙烧成型 预活化销售与使用 以金属催化剂为例 催化剂合成（制备）： 固体状态化学、表面胶 体化学、无机化学 、有机金属化学 催化剂表征：结构、性质、光谱、衍射 非均相催化动力学：反应动力学（宏观水平） 反应机理（微观水平） 工程应用,19,1.6 建立化学反应动力学的步骤,1.6.1 对反应过程的认识 热力学分析、热力学图式（方向、限度、可逆性），可能的反应产物（往往工艺开发分析中已经解决）,1.6.2 基本物料衡算关系 化学反应计量学化学计量独立反应（独立反应数和计量独立反

10、应的选取） 关键组分（加合性原则） 相态变化（恒容、变容、相态）,20,1.6.3 实验装置和分析方法 预备实验、计量标定等,1.6.4 动力学图式的确定 动力学独立数和动力学独立反应,1.6.5 反应条件的确定 考察的目的 内外扩散检验,21,1.6.6 实验设计 单因素法、正交实验法、序贯法、均匀设计等,1.6.7 实验数据处理（将采集量转化为常规表示） a. 将实验获得数据（色谱、容量分析等）转化为常规（习惯）数据 b. 图、表给出数据细节（完整） c. 根据需要给出图表（简化细节） d. 动力学处理需要数据 一次数据（实验获得的最原始数据），二次数据（a），三次数据（b），四次数据（c

11、），五次数据（d）,22,1.6.8 动力学模型及参数估计 模型的选取 数值计算方法的确定（目标函数、数值方法、计算机程序等） 估值参数的初步筛选,1.6.9 动力学模型检验 与实验数据方面的检验（故意留一些不参与估值的数据，待模型估值后用以检验 物理化学准则检验 秩和检验、参差分布及f（或2）统计检验 工业应用检验,23,1.7 小结,1.7.1 化学反应动力学分类,化学反应动力学,宏观动力学 （过程动力学传递过程化学过程）,吸附动力学,表面反应动力学,脱附动力学,内扩散化学动力学,外扩散化学动力学,化学动力学 内扩散外扩散,本征动力学 （微观动力学/化学动力学）,活化动力学,失活动力学,宏

12、观失活动力学,常见本征动力学涵盖内容,稳（定）态动力学,瞬态动力学,24,1.7.2 化学反应动力学模型建立方法,化学反应动力学模型建立,物种分析（gc-ms等表征）,网络结构判识 （反应图式）,热力学分析,动力学分析,实验装置 测试方法与措施,催化剂性能考察、动力学测试区域、装置性能、实验技术、措施等,实验设计、实验安排、数据分析,动力学模型建立,动力学模型筛选/推导,参数优化估计,动力学模型检验,物理/化学准则检验,统计检验,应用检验,25,26,化学反应动力学模型建立方法,化学反应动力学模型建立,物种分析（gc-ms等表征）,网络结构判识 （反应图式）,热力学分析,动力学分析,实验装置

13、测试方法与措施,催化剂性能考察、动力学测试区域、装置性能、实验技术、措施等,实验设计、实验安排、数据分析,动力学模型建立,动力学模型筛选/推导,参数优化估计,动力学模型检验,物理/化学准则检验,统计检验,应用检验,27,多相催化反应动力学建立步骤,（1）对反应过程的认识 （2）基本物料衡算关系 （3）实验装置和分析方法 （4）动力学图式的确定 （5）反应条件的确定 （6）实验设计 （7）实验数据处理 （8）动力学模型及参数估计 （9）动力学模型检验,28,2 多相催化反应动力学基础理论,2.1 催化剂基本组成及其催化特征 2.1.1 催化剂分类：主要分为四类,29, 良好的催化活性（尤指低温活

14、性） 良好的选择性 较长的使用寿命 适宜的物理织构（sg、vg、孔径分 布、活性组分分布、形状、尺寸） 较强的抗毒能力 较高的机械强度（不易磨损、破碎）,2.1.2 工业催化剂的使用要求,30, 催化剂定义 参与反应，但反应前后的组成和特性不变，对正逆反应均等倍加速或减速,2.1.3 工业催化剂的性质和特点, 催化剂性能, 直接性能（活性、选择性和稳定性） 附加性能（强度、热稳定性和寿命等）,31,32,2.1.4 催化剂组成及各组分的作用, 活性组分 金属类 半导体类 酸碱类 绝缘体类iiiv、iva和va金属或非金属氧化 物，卤化物等, 助催化剂（助剂） 1 结构型助剂 作用：增加主活性组

15、分活性构造的稳定性和活性组分的比表面积,33,按作用机理分类： (a) 分隔（离）活性组分细小微晶高熔点、难还 原的氧化物比较适宜 (b) 与活性组分作用生成高熔点的化合物或固（熔） 体。如fe3o4中加入少量的al2o3（结构助剂）,2 电子型助剂 k2o、na2o不改变催化剂本性，但改变活性组分电子结构，从而使反应分子的化学吸附能力和反应活化能发生变化,3 扩散型助剂 改善催化剂孔结构（造孔剂）如分子筛中加入有机酸等（织构结构）,34,4 毒化型助剂 毒化催化剂某些有害的活性中心，如为防止结焦，可加少量碱性物质毒化引起裂解生焦的酸性中心, 载体 作用：改变催化剂机械强度、导热性和热稳定性，

16、增大活性表面和提高适宜孔结构，提供活性中心减少催化剂活性组分用量 常见类型（起稳定剂和分散剂的作用）：al2o3、mgo、硅胶、硅藻土等, 抑制剂 （乙烯环氧化二氯乙烷的作用）,35,2.1.5 催化剂制备与各制备过程的作用（目的）,1 原粉制备 催化剂使用 5个使用区域 原粉制备方法 （1）浸渍法载体放入活性组分溶液中（盐、碱、 酸等）浸渍 （2）共沉淀法 （3）（机械）混（共）合法,原料预处理原粉制备焙烧成型预活化销售,36,（4）沥滤法骨架催化剂制备。骨架ni、骨架co、 骨架fe、骨架cu。如ni催化剂，将 含ni和al各50%的铝镍合金破碎、 过筛后用20%naoh溶液溶解al， 留

17、下具有高比表面积的ni （5）其它方法：热熔融法、电解法、粘结剂法等,2 焙烧过程描述 焙烧目的 (a) 稳定催化剂组成和活性使催化剂具有一定化学组成和稳定的活性 (b) 使催化剂具有（保持）一定的晶型、粒度大小、孔结构和比表面积,37,(c) 增强或减弱机械强度使催化剂具有一定机械强度, 焙烧过程变化 (a) 化学变化 eg. agno3 ag +1/2o2 + no2 (b) 催化剂结构变化 (c) 颗粒变化（晶型转变等） eg. 水软铝石转化为al2o3 (d) 孔结构变化,38,3 成型 压缩成型（打片） 挤出成型（挤条） 转动成型（糖衣机） 喷雾成型（喷雾干燥原理） 其它成型方法（油

18、中成型法、喷动成型、 蜂窝及纤维催化剂成型）,39,2.2 吸附作用基本原理,2.2.1 工业多相催化反应历程（以固体催化 剂为例）动力学模型选定时应考虑的过程 目的是获得动力学方程，但应首先清晰研究过程的目的（本征、宏观），对催化过程的认识是必须的，而吸附是发生多相催化反应的重要基础。催化剂多孔介质，外表面积多孔介质丰富的内部微孔。 例如ag催化剂上乙烯氧化制环氧乙烷c2h4 + o2 c2h4o，深层氧化主要是生成co和co2，因而采用非多孔性介质即，尽量控制反应在催化剂表面上进行，以提高环氧乙烷的选择性。,40,多孔介质上的多相催化反应历程,41,多孔介质上的多相催化反应历程 外扩散反应

19、组分由主体物流 内扩散反应组分由催化剂外表面 吸附反应组分在催化剂内表面活性中心上吸 附（以化学吸附为主，物理吸附为附） 表面化学反应 脱附产物 内扩散产物由催化剂内表面 外扩散产物由催化剂外表面 其中，步骤（3）、（4）和（5）为化学动力学过程（本征动力学）,42,2.2.2 固体催化剂颗粒的主要结构参数 催化剂性能主要由活性、选择性和寿命组成，这些均与催化剂结构参数直接或间接联系（如分子筛选择性、寿命和内表面利用率等）。 主要结构参数：比表面积、孔体积、孔体积（孔径）分布、固体密度、颗粒密度、孔隙率、孔内扩散描述的非特征参数曲节因子（迷宫因子）,1 比表面积 sg（m2/g）单位质量催化剂

20、所具有的表面积，多层吸附理论bet法测定。一般为51000 m2/g,43,2 孔体积（孔容积） vg（m3/g）单位质量催化剂所具有的微孔体积 测定： (a) 水煮法 (b) 氦汞法 (b1) 氦法：测定的是固体物质占据的体积vt (b2) 汞法：测定的是固体物质 + 微孔占据的 体积vs,3 固体密度（真密度）不含微孔时的密度 固体质量 mp,44,4 颗粒密度（假密度）注意与床层堆密度的区别 固体质量 mp,5 颗粒孔隙率微孔所占总体积的分率,注意：床层空隙率 b,45,6 孔径分布测定 压汞法测大孔 1210-8 m （即 r 1210-8 m 100200 以上） 氮吸附法测微孔 1

21、10-9 m r 210-8 m （即 10200 ） 1 = 10-8 cm ,46,2.2.3 多孔介质吸附脱附理论,1 吸附类别 物理吸附,(a) 依靠范德华力（分子间引力），类似于蒸汽的凝结或气体的液化，无明显选择性 (b) 多层吸附 (c) 吸附热大致与被吸附组分的冷凝潜热相当 2.121 kj/mol（0.55 kcal/mol） (d) 吸附与脱附能快速达到平衡，可逆，随温度升高，吸附量降低,47, 化学吸附,(a)由固体表面与吸附分子间的化学键力所致，形成吸附化学键，产生表面活化络合物 (b)多为单分子层吸附，吸附满后不再吸附（即吸附量有上限） (c)吸附热与化学反应热相当（同

22、一数量级）21418 kj/mol（5100 kcal/mol） (d)随温度升高，吸附量增加,48, 物理与化学吸附对比 （1）物理吸附尽管能使物质吸附在催化剂表面，但不改变活化能，因而，不是构成催化过程的主要原因； （2）化学吸附因化学键的作用，使被吸附分子发生变化，更显活泼,性，是多相催化过程的主要原因； （3）低温下以物理吸附为主，高温下以化学吸附为主。,49,2 影响吸附（脱附）速率的因素,吸附过程：a + a活性中心 脱附过程：a a + 表面覆盖度（组分吸附速率）： 空位率： 单一组分吸附时： 多组分吸附时： ,50,1 影响吸附速率因素 (a) 与吸附分子对表面的碰撞数成正比

23、气体（流体）对固体表面碰撞频率越快，吸附速率越快。据分子运动论，组分a碰撞频率/时间与其分压成正比,paa气体分压；maa气体分子量；kbboltzmann常数；t温度 k,(b) 与吸附分子碰撞在空吸附位（自由表面）上的几率成正比。其为覆盖率的函数,51,(c) 与吸附活化能成指数关系。占有高于吸附吸附活化能ea的那些分子才能被吸附，它们在全部气体分子中所占分率为exp(-ea/rt)（有效分率）, 吸附速率,按化学反应 a + a,凝聚系数，代表具有能量超过ea的分子撞击在空吸附位上而被吸附的几率,52,2 影响脱附速率因素 (a) 与已经吸附了气体分子的吸附位占总吸附位的分率成正比。此为

24、覆盖率的函数,(b) 与脱附活化能ed呈指数关系。而ed与a有关 气体（流体）对固体表面碰撞频率越快，吸附速率越快。据分子运动论，组分a碰撞频率/时间与其分压成正比。 脱附速率,53,按化学反应 a a + ,净吸附速率,吸附热 ed=ea+c。 多数化学反应为放热过程，所以，大多数情况下c为负值,54,3 吸附等温方程（吸附量的确定）,吸附平衡时 r=0,吸附平衡常数,ea和ed确定困难，有多种模型。ka（吸附平衡常数）表征固体表面对气体分子吸附能力的强弱，是一个吸附特征值，与温度的关系为（一般qc 0 ）,55,基于对qc的不同处理，得出不同吸附模型（实际为ea和ed与a的关系）。 在温度

25、恒定时，达到吸附平衡条件下，吸附量（常以体积计）与压力的关系曲线叫等温吸附线。尚有：等压吸附线和吸附等容（量）线。 吸附量 ： 等温吸附模型 (a) 理想吸附(langmuir)均匀表面吸附理论：qc、ea和ed恒定 (b) 真实吸附模型（不均匀表面吸附模型）,56,b1 temkin（焦姆金）线性关系 ea=ea0+a ed=ed0-a qc=qc0-a,b2 freundlich（弗朗得里希）对数关系 ea=ea0+lna ed=ed0-lna qc=qc0-flna,57,4 均匀表面吸附理论（langmuir模型）,理想吸附模型/均匀表面吸附模型 基本假定： 1 吸附剂表面均一。即，整

26、个催化剂表面具有均匀的吸附能力。 2 吸附和脱附活化能与表面吸附的程度无关（ea和ed恒定） 3 单分子层吸附（每个活性中心仅能吸附一个气体分子） 4 被吸附分子间没有相互作用 note：吸附和脱附在一定条件下可达平衡,58,（a）单分子吸附,吸附速率： ra=ka pav 脱附速率： rd=kd a 净吸附速率：r= ra rd= ka pav kd a 而： av 1 v 1- a r= ka pa(1- a) kd a 吸附平衡时： r= 0,吸附平衡分压,59,（b）单分子解离吸附,吸附速率： ra=ka pav2 脱附速率： rd=kd a2 净吸附速率：r= ra rd= ka p

27、av 2 kd a2 而： av 1 v 1- a r= ka pa(1- a)2 kd a2 吸附平衡时： r= 0,60,（c）双分子吸附,a的净吸附速率：ra= raa rda= kaa pav kdaa 吸附平衡时： a kapa*v b的净吸附速率：rb= rab rdb= kab pbv kdbb 吸附平衡时： b kbpb*v 而 ab + v 1,61,（d）多分子吸附,i的净吸附速率：ri= rai rdi= kai piv kdii 吸附平衡时： i kipi*v 而 iv 1,62,吸附极限/端情况（渐近线）, 稀疏覆盖 单组分：if kapa*1 a 1,63,64,3

28、 热力学分析在动力学研究中的应用, 热力学分析目的： 方向、限度、热效应和热力学图式, 封闭体系 化学反应热,生成热 燃烧热 总热效应,65,例题3-1：一实验装置采用羰基合成工艺生产乙醇，化学反应方程式为：2co(g)+4h2(g)c2h5oh(g)+h2o(g)。已知下列化合物的生成热及燃烧热：,生成热 （下标f表示）,燃烧热 （下标c表示）,试求1atm、25oc时，乙醇合成反应的化学反应热。,高热值、低热值概念？,66,解：（1）通过气态甲醇的燃烧热求气态甲醇的生成热 气态乙醇完全燃烧的化学反应方程式为： c2h5oh(g)3o22co2(g)3h2o(l) 则，该反应1atm、25o

29、c时下的化学反应热,67,（2）乙醇合成反应的化学反应热的计算 化学反应方程式为： 2co(g)+4h2(g)c2h5oh(g)+h2o(g) 则，该反应1atm、25oc时下的化学反应热,68, 化学平衡常数与温度的关系,vant hoff 方程,计算gibbs自由标准生成焓,经验公式 计算法（内插、外推） 增量法（基团贡献法）,对称数。被认为是刚性的分子围绕一个（或多个）任意对称轴旋转时所能占据的同等位置的数目，如，对甲烷=12；对乙烷=2,69, 化学平衡组成计算 依据：化学平衡 元素平衡,将多反应体系分成若干个单一反应（即在不同反应器内进行不同的反应）,70,例3-2 试从化学热力学角

30、度出发，分析下列两个反应在800k时哪个更易进行？并计算相应的化学平衡常数。 c6h12c2h4+c4h8 （1） c6h12c6h6+3h2 （2） 已知800k时， goc6h12=319.2kj/mol ； goc2h4=103.9kj/mol goc6h6=221.5kj/mol ； goc4h8=209.3kj/mol,解1,71,解2,72,例3-3：甲烷与水蒸汽混合物中水碳比为5，设转化过程没有炭黑析出，试计算1atm和600oc条件下的平衡组成。已知该条件下的主反应及平衡常数如下： 反应1：ch4+h2o=co+3h2 kp1=0.9618 反应2：co+h2o=co2+h2

31、kp2=2.7,解：以1mol入料ch4为基准，设x为按反应1转化了的ch4摩尔数，y为按反应2变换了的co摩尔数，m为水碳比，则平衡物料关系可写为：,反应1：ch4+h2o=co+3h2 始态： 1 m 0 0 终态： 1-x m-x x 3x,反应2：co+h2o=co2+h2 始态： x m-x 0 3x 终态： x-y m-x-y y 3x+y,结果：ch4 1-x；h2o m-x-y ；co x-y；h2 3x+y； co2 y 系统总摩尔数 1+m+2x,73,反应1：ch4+h2o=co+3h2 始态： 1 m 0 0 终态： 1-x m-x x 3x,反应2：co+h2o=co

32、2+h2 始态： x m-x 0 3x 终态： x-y m-x-y y 3x+y,结果：ch4 1-x；h2o m-x-y ；co x-y；h2 3x+y； co2 y 系统总摩尔数 1+m+2x；总压 p,74,经整理，即可得到如下的平衡组成,迭代序列 构建技巧,75,例3-4 乙烯在银催化剂上用过量空气部分氧化制取环氧乙烷，化学反应方程式为：主反应：2c2h4+o2=2c2h4o；副反应：c2h4+3o2=2co2+2h2o。反应器进口混合气体的流量为45 000nm3h-1，原料组成mol分率为 c2h4 3.5%，n2 82.0%，o2 14.5%。乙烯转化率为32%，反应选择性为69

33、，则反应器出口气体的流量为（ ）kgh-1。 （1）57 107.5； （2）57415.2； （3）56 981.95； （4）11 289.9,解：该题干扰因素明显，表面上看似乎化学反应对过程计算非常必要且难以计算，但实际上，由于题目要求计算的是反应器出口混合气体的质量，而题设条件除反应器进出口流股外无其它信息，因此，根据质量守恒定律可知，反应器进出口的混合气体质量流量相等，而化学反应的存在仅仅是元素的重排。,76,解：（1）反应器进口平均分子量计算,（2）反应器进、出口质量流量计算,故：反应器出口气体的流量为57415.2 kgh-1,77,化学反应计量学是研究化学反应进行时，反应混合物

34、组成变化规律的科学，其基础是化学反应计量式,4 化学反应计量学和动力学反应图式的判识,4.1 化学反应计量学的概念和依据,甲醇合成,习惯表示方法co2h2 =ch3oh（1）,物种摩尔数变化关系,化学计量式co2h2 ch3oh 0（2）,反应物为，产物为。计量上，（1）和（2）无本质差别,78,不失一般性，对n个ai物种参与的k个反应的复杂反应系统（j反应），化学反应方程式为：,同样移项得化学计量关系（文献常见形式）,ji计量系数，反应物为“”，生成物为“+”（得到为正，失去为负）,79,需要解决的问题 如何对这一复杂反应体系进行描述？ 描述的途径是什么？（组成，动力学，路径等）, 静态问题

35、： 确定混合物组成时，至少知道几种物质的量 关键组分数及其选取 确定组分变化量时，至少需要几个化学反应 化学计量独立反应数及其选择, 动态问题： 组分变化量与各个反应进行程度之间如何关联 nij 化学计量独立反应是否为动力学独立反应 动力学独立反应数及其选择,80,4.2关键组分与化学计量独立反应（原子矩阵法）,封闭体系中，l个元素组成的n个物质的化学反应过程遵从原子守恒定律（元素平衡）,i=1n，h=1l（元素），bhh元素的原子数,但化学反应是按摩尔数变化进行的，ni=ni-ni0,hi的矩阵b有（ll）个矩阵元，设其矩阵的秩为r。则该体系的自由未知数（自由度）为,r=n-r（r关键组分）

36、 一般，rv=r（rv化学计量独立反应数）,化学计量学基础,81,例1关键组分和化学计量独立反应确定,甲烷水蒸气转化制合成气。 物质ch4、h2o、co、co2、h2、c； 元素c、h、o；物种数i=16，元素数h=13。确定并选取关键组分和化学计量独立反应。, 化学反应式确定（热力学分析）,82,关键组分数和化学独立反应数确定 元素原子守恒,经初等变换可得秩为：r=3,以及,83,关键组分数和化学独立反应数为,根据反应原料组成和目的产物，选取如下的化学计量独立反应： ch4+h2o=co+3h2(1) co+h2o=co2+h2 (6) ch4=c+2h2 (7),无析炭反应发生时，r=3，

37、r=5-3=2,选取： ch4+h2o=co+3h2(1) co+h2o=co2+h2 (6),84,例2化学计量独立反应和物料衡算,甲醇合成。原料 co/co2/h2，产物： co/co2/h2 ，ch3oh，h2o和n2+ch4。确定并选取化学计量独立反应，并进行基本物料衡算。, 主反应 co+2h2=ch3oh co2+3h2= ch3oh +h2o co2+h2 =co+h2o,秩：r=3 ；关键组分数和化学独立反应数为：,85, 化学计量独立反应 co+2h2=ch3oh co2+3h2= ch3oh +h2o, 选取co和co2为关键组分进行物料衡算,m表示甲醇,86,m表示甲醇,

38、co+2h2=ch3oh co2+3h2= ch3oh +h2o,87,88,作业：试推导选取ch3oh和co2为关键组分时的结果。答案如下：,89,4.3组分变化量与反应程度的关系,反应程度：,广义量,单一反应：化学计量式 反应程度 反应速率,90,复合反应： 化学计量式 反应程度 反应速率,j=1m（m动力学独立反应），即rv=n-r1，组分ai摩尔数变化为多个反应j造成的，多数情况下，mrv。,91,4.4计量学与反应动力学之间的关系,化学动力学的目的：根据反应系数随时间的变化导出反应进行程度或反应机理的、有物理和化学意义的规律性。相反，大多数从计量学角度考虑的反应没有动力学意义（随机选

39、定） 动力学独立反应数是否一定等于化学计量独立反应数（或说关键组分数）？即，是否一定存在,回答：rv不一定总等于m,92,动力学图式分析 条件：a1、a2和a3为同一分子式组成的同分异构体，原料中有大量的a1和少量a2和a3，判别可能的动力学反应路径（前三个：实际反应数等于关键组分数）,途径1 途径2 途径3 途径4,93,一个特别的例子蔗糖的转化反应 （实际反应数小于关键组分数）,rv计算：相应的元素物种矩阵的秩为2，计量学上应该有两个（4-2=2）线性无关的反应方程式。第二个关键反应也许是异构化反应：,c12h22o11（蔗糖）+h2o=c6h12o6（葡萄糖）+ c6h12o6（果糖）,

40、c6h12o6（葡萄糖）c6h12o6（果糖）,但从化学上看，这种异构反应是不可能的,因此，计量学只为关键反应数提供了上限,94,总结：无论是为了选择化学上很有意义的关键反应，还是为了确定它们的数目，都需要进行动力学实验。计量学只是对系统的数学处理划出了外部范围只为关键反应数提供了上限。,95,4.5动力学路径分析实例,例：无定形c（a1）石墨（a2）或金刚石（a3），提出问题（1）三种物质均由c构成，计量独立化学反应数？（2）以纯原料a1进行变空速实验发现，低空速下产物中只有a2，而无a3生成；高空速下，a1、a2和a3并存，动力学图式是什么？,解：（1）计量独立化学反应数？ 方法一：对于c

41、、h、o体系，关键组分数 物种数组成这些物种的元素个数3-1=2,方法二：,96,4.5动力学路径分析实例,（2）动力学路途判别 可能的路径,途径1 途径2 途径3 途径4,低空速下只有产物a2存在反应,高空速下a2和a3并存，存在反应,结合低空速实验结果，a3也可能生成a2，因此，可能的反应路径为,或,97,例 苯与异丙醇气相烷基化合成异丙苯,实验确定苯与异丙醇气相烷基化合成异丙苯的动力学反应网络结构。原料由苯（b）和异丙醇（ipa）组成,硕士论文题目 ysbh-02催化剂上苯与异丙醇气相烷基化 合成异丙苯工艺条件及本征动力学研究,98, 实验流程,产物分布测定和热力学分析 gc-ms测定表

42、明，产物主要由丙烯（pr）、水、异丙苯（ipb）、少量二异丙苯（dipb）、微量正丙苯（npb）和微量三异丙苯（tipb）组成,99,例 苯与异丙醇气相烷基化合成异丙苯,主反应：ipa=h2o+pr（1） pr+b=ipb（2） ipa+b=ipb+h2o（3） 副反应：2ipb=dipb+b（4） ipb=npb（5） ipb+pr=dipb（6） ipb+ipa=dipb+h2o（7） prdipb=tipb （8） 除反应（2）和反应（5）为热力学可逆反应外，其它6个反应均可看作正向不可逆反应,100,异丙醇与苯烷基化反应网络结构的提出,ipa异丙醇，pr丙烯，b苯，ipb异丙苯，dip

43、b二异丙苯，dipb三异丙苯,停留时间对产物分布的影响,特别注意: 动力学图式反应机理,101, 过程组成：外扩散过程、内扩散过程和化学动 力学过程三部分组成（共7个步骤）,5 物理传递过程对化学反应速率的影响, 过程描述的关键基础 1 控制步骤存在与否 2 操作状态（定态或非定态） 3 过程描述基础理论（传递和反应）,一、多相催化过程中的物理传递现象及其描述,102, 外扩散控制 内扩散控制 化学动力学控制（吸附、脱附、表面化学反应） 共同控制,二、多相催化过程控制步骤判断,三、多相催化过程宏观速率的描述 效率因子的概念, 过程速率定量方法（宏观动力学） 1 宏观活性实验测定及其模型化 2 依据本征动力学和反应-扩散模型计算,103,5.1 多相催化过程中的物理传递现象及其描述,104,外扩散控制 cag cas cac ca* 内扩散控制 cag cas cac ca* 动力学控制 cag cas cac ca* 共同控制？,105,5.1.1 外扩散过程,类型：总体流动、粘滞流流动、返混和扩散等 定态速率：(ra)g=kgse(cag-cas) gglobal；kg外扩散传质系数，m/s； se单位体积催化床中颗粒的外表面积，m2/m3； (ra)g