第五章气固反应.ppt

《第五章气固反应.ppt》由会员分享，可在线阅读，更多相关《第五章气固反应.ppt（108页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第五章第五章 固定床气固定床气-固相催固相催 化反应工程化反应工程覃吴覃吴内内 容容第一节第一节 固定床气固相催化反应器的基本类型和数学模型固定床气固相催化反应器的基本类型和数学模型第二节第二节 固定床流体力学固定床流体力学第三节第三节 固定床热量与质量传递过程固定床热量与质量传递过程第四节第四节 绝热绝热式式固定床反应器固定床反应器第五节第五节 连续换热内冷自热式催化反应器连续换热内冷自热式催化反应器第六节第六节 连续换热外冷及外热管式催化反应器连续换热外冷及外热管式催化反应器第七节第七节 薄床层催化反应器薄床层催化反应器第一节第一节固定床气固相催化反应器的基本类型和数学模型固定床气固相催化

2、反应器的基本类型和数学模型1.固定床反应器：固定床反应器：凡是流体通过不动的固体物料所形成的床层而进行反应的装置。2.固定床反应器优缺点：固定床反应器优缺点：催化剂不易磨损；床层内流体平推流，较少量催化剂可获较大生产能力；传热较差。3.固定床反应器的分类固定床反应器的分类固定床反应器不同的传热要求和传热方式绝热式二段三段四段冷激式反应特征单段绝热式段间反应气冷却或加热方式中间间接换热式多段绝热式原料气冷激式非原料气冷激式换热式对外换热式自热式加压热水（240）导热油（250300）熔盐（300）反应气的流动方向轴向流动固定床反应器径向流动固定床反应器4.数字模型数字模型描述固定床反应器的数字模

3、型分为：拟均相和非均相两类拟均相模型忽略床层中粒子与流体间温度与浓度的差别。非均相模型考虑床层中粒子与流体间温度与浓度的差别。第一节第一节固定床气固相催化反应器的基本类型和数学模型固定床气固相催化反应器的基本类型和数学模型1.固定床的物理特性固定床的物理特性(1)颗粒密度p床层密度或堆积密度B(2)当量直径及形状系数 具有相同体积的球粒子直径dv 来表示dv？具有相同外表面积球粒子直径Dp来表示Dp？具有相同比外表面积球粒子直径ds来表示ds？形状系数=Ss/Sp 第二节第二节固定床流体力学固定床流体力学(3)混合颗粒平均直径(4)固定床当量直径及空隙率 de=4RH=4/Se=2/3*(/1

4、-)*ds(5)空隙率及径向流速 了解即可第二节第二节固定床流体力学固定床流体力学dp=1/n（i 1=xdii)2.单相流体在固定床中的流动及压降单相流体在固定床中的流动及压降(1)流动(了解即可)(2)压降 第二节第二节固定床流体力学固定床流体力学(3)压降影响因素(了解即可)修正雷诺数：修正雷诺数：um 平均流速(空床气速)；L 床层高度；、流体的密度和粘度；ds 比表面当量直径。影响固定床压力降的因素影响固定床压力降的因素流体流体流体的密度流体的密度流体的粘度流体的粘度流体的质量流率流体的质量流率床层床层床层的高度床层的高度床层的空隙率床层的空隙率流通截面积流通截面积颗粒颗粒颗粒的形状

5、颗粒的形状颗粒的粒度颗粒的粒度颗粒的表面粗糙度颗粒的表面粗糙度颗粒的物理特性颗粒的物理特性床层压降计算例：在内径为50mm的管内装有4m高的催化剂层，催化剂的粒度分布如下表。催化剂为球体，空隙率B0.44。在反应条件下的气体密度=2.46kg/m3，粘度=2.310-5kg/(ms)，气体的质量流速G=6.2kg/(m2s)。求床层压降。粒度粒度 d ds s/mm/mm3.403.404.604.606.906.90质量分率质量分率 w w0.600.600.250.250.150.151.固定床中的传热固定床中的传热传热包括：粒内传热，颗粒与流体间的传热(给热系数h hp p)，床层与器壁

6、的传热(给热系数h hw w).).以床层的平均温度 与管壁温差为推动力-总给热系数 h hc c有效导热系数有效导热系数 e e轴向轴向 eses径向有效导热系数径向有效导热系数 erer=e0e0+(+(erer)t t第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程传热方式：导热、对流传热、热辐射。传热途径：粒内传热、颗粒与流体间传热、床层与器壁间传热。固定床的径向传热方式床层空隙内部流体的传热颗粒之间通过接触的传热颗粒表面附近流体中的传热空隙内部流体的辐射传热流体混合所引起的径向对流传热颗粒表面之间的热辐射传热第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程一、

7、颗粒与流体间传热系数一、颗粒与流体间传热系数1、传热因子适用范围颗粒与流体间的传（给）热系数hpdp可用dv代注意：在参考其他教材时，给热系数多以表示。第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程计算传热系数的经验公式有很多，可从有关文献或工具书中查取。应用时要注意公式规定的特性尺寸、特性温度和适用范围。第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程式中：单位重量催化剂的外表面积 床层比表面积Se的校正系数主体、外表面温度2、流体与颗粒间传热温差的计算、流体与颗粒间传热温差的计算热量平衡热量平衡第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程上式可整理

8、成上式可整理成其中，传热量其中，传热量普兰特数普兰特数根据6-12的关联图，查图可求得不同条件下的t。第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程图图6-12固定床中流体与颗粒外表面温差固定床中流体与颗粒外表面温差t0rAHrtG时时te 是针对拟均相模型提出的综合性传热参数，一般是指径向热导率er二、固定床的有效热导率二、固定床的有效热导率有效导热系数有效导热系数 e e轴向有效导热系数轴向有效导热系数 eses径向有效导热系数径向有效导热系数 erer=e e0 0+(+(erer)t t第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程式中，式中，流体的热导率流

9、体的热导率流体静止时床层的热导率流体静止时床层的热导率径向与轴向传质速率之比径向与轴向传质速率之比颗粒间距与粒径比的影响颗粒间距与粒径比的影响值可由图值可由图6-14查取。查取。粒径粒径/管径管径图图求有效热导率求有效热导率e第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程包含床层包含床层空隙和颗粒空隙和颗粒对传热的贡献，由下式计算对传热的贡献，由下式计算式中，分别表示颗粒与流体的热导率；床层的平均温度；空隙的辐射给热系数；颗粒的辐射给热系数；颗粒接触点处流体薄膜对导热的影响。第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程床层空隙的辐射给热系数：hrv式中，粒子表面的

10、热辐射率；床层的平均温度；空隙的辐射给热系数；第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程颗粒的辐射给热系数：颗粒接触点处流体薄膜对导热的影响。由图由图6-15查取。查取。式中，粒子表面的热辐射率；床层的平均温度；颗粒辐射给热系数；第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程由图由图6-15查取。查取。图图6-15求所用的值求所用的值当颗粒直径甚小，床层温度不是很高，以及含有液体时，空隙和颗粒的辐射传热可忽略，床层空隙和颗粒对传热的贡献可简化为：第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程小结：固定床的优缺点固定床的分类基本概念 颗粒直径的表示方

11、法 颗粒的形状系数 混合颗粒的平均直径 床层空隙率 床层比表面积 床层当量直径固定床的径向流速分布床层压降厄根(Ergun)方程修正雷诺数影响固定床压力降的因素固定床中的传热颗粒与流体间传热系数第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程三、床层与器壁间的给热系数三、床层与器壁间的给热系数h0 一维模型中，床层与器壁间传热速率为 h0可由经验公式计算式中，y 无量纲数（适用范围：y 0.2）：床层平均温度;：器壁温度;第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程b 无量纲数均为无量纲的b的函数hW 壁面处对流给热系数；第三节第三节固定床热量与质量传递过程固定床热

12、量与质量传递过程适用范围：适用范围：由图由图6-16查得。查得。第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程图图6-16求求h0例如：当例如：当b12.98时，时，分别对应分别对应5.02和和0.25b由于催化剂表面存在滞流边界层，气流主体浓度与催化剂颗粒表面浓度存在差异。在滞流层内有浓度差，必然存在扩散。一、颗粒与流体间的传质系数一、颗粒与流体间的传质系数 单位体积（或质量）催化剂上关键组分A的传质速率第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程距离0CACAgCASCACRp整个传质方程的核心，总包了各种条件对传质的影响。由实验关联式计算。关联式之一：传质因

13、子：传质因子：JD实验关联式：实验关联式：气相密度：气相密度;：气体质量流率：气体质量流率;：气相粘度：气相粘度;：气相分子扩散系数：气相分子扩散系数;JD是雷诺数的函数是雷诺数的函数第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程Sc为施密特准数为施密特准数um：平均流速(空床气速)；g、g：气体的密度和粘度；ds：比表面当量直径。当修正雷诺数在：Sc为施密特准数为施密特准数第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程比表面当量直径：单颗催化剂的外表面积折合成直径为dS的球形颗粒应有的外表面积。传质对反应的影响达姆科勒数 第三节第三节固定床热量与质量传递过程固定床

14、热量与质量传递过程为外部效率因子，或称外部有效因子，它是D的函数，反映界面反应物浓度降低对反应速率的影响程度Da为反应速率与扩散速率的比值，反映了体系中外扩散的影响程度。数值越大，或反应速率越快，外扩散的影响就越大。对一级反应第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程特殊情况：反应速率常数k比传质系数kg大得多，则颗粒外表面处A的浓度为零，属外扩散控制。反应速率常数k比传质系数kg小得多，则颗粒外表面处A的浓度与气相主体浓度相等，属内扩散或动力学控制，外扩散可不予考虑。第三节第三节固定床热量与质量传递过程固定床热量与

15、质量传递过程传热因子传热因子传质因子传质因子：气相密度 ：气体质量流率 ：气相粘度 ：气相分子扩散系数 JD是雷诺数的函数施密特准数第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程由传质和传热的类比原理有由传质和传热的类比原理有上式在缺少数据时用来推算数据十分有用上式在缺少数据时用来推算数据十分有用.无论是传质或是传热系数无论是传质或是传热系数,增加质量速率增加质量速率G都可以加快流都可以加快流体和催化剂外表面间的传递速率体和催化剂外表面间的传递速率,但相应的床层压降也但相应的床层压降也增加增加.外扩散速率外扩散速率,内扩散速率内扩散速率,表面反应速率的相对大小是决表面反应速率的相

16、对大小是决定反应控制步骤的关键定反应控制步骤的关键.第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程外扩散过程对表面温度的影响外扩散过程对表面温度的影响由传质计算可得：由传质计算可得：第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程两式相等可得：两式相等可得：第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程上式将流固相的温度差与浓度差联系起来了。上式将流固相的温度差与浓度差联系起来了。进一步简化，前面有：进一步简化，前面有：极为相似。极为相似。JHJD相除相除第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第三节第三节固定床热量与质量传递过程固定

17、床热量与质量传递过程据实验得：据实验得：第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程例：苯加氢反应器在例：苯加氢反应器在1013.3kPa下操作，气体质量速度下操作，气体质量速度G=3000kg.m-2h-1，催化剂为，催化剂为89mm圆柱体，颗粒密度圆柱体，颗粒密度P=0.9g.cm-3，床层堆积密度，床层堆积密度B=0.6g.cm-3，在反应器某处，在反应器某处气体温度为气体温度为220，气体组成为，气体组成为10苯，苯，80氢，氢，5环己环己烷和烷和5甲烷甲烷(体积分率体积分率)，测得该处宏观反应速率，测得该处宏观反应速率(-rA)=0.015mol.h-1g-1(cat

18、)。估算该处催化剂的外表面浓度。估算该处催化剂的外表面浓度。已知：气体粘度已知：气体粘度=1.410-4g.cm-1s-1，综合扩散系数，综合扩散系数D=0.267cm2s-1。第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程解：解：（1）计算催化剂的粒径）计算催化剂的粒径dS。（2）计算床层中气体的修正雷诺数。）计算床层中气体的修正雷诺数。（3）计算）计算JD和和kg值。值。（4）计算）计算cAG和和cAS习题：试计算例题中催化剂的外表面处温度。习题：试计算例题中催化剂的外表面处温度。已已知知反反应应热热为为(-H)=2.135105Jmol-1，气气体体的的定定压压比比热热容容

19、cP=49Jmol-1K-1。小结：颗粒与流体间的传（给）热系数hp传热因子JH流体与颗粒间传热温差的计算固定床的有效热导率床层与器壁间的给热系数 h0传质因子JD6-3 6-3 催化反应器的数学模型催化反应器的数学模型1 1，非均相拟均相，非均相拟均相2 2，一维模型二维模型，一维模型二维模型3 3，理想流动非理想流动，理想流动非理想流动拟均相适用情况：拟均相适用情况：1 1，化学动力学控制，化学动力学控制2 2，活性较正系数（无宏观动力学资料），活性较正系数（无宏观动力学资料）一维二维：轴向浓度差、温度差；轴径向浓度差、温度差一维二维：轴向浓度差、温度差；轴径向浓度差、温度差理想流动：不考

20、虑返混（理想流动：不考虑返混（PFR）；非理想流动：考虑返混（扩散）非理想流动：考虑返混（扩散）6-3 6-3 催化反应器的数学模型催化反应器的数学模型维数维数相数相数返混返混根据相态(拟均相？非均相？），维数（一维？二维？），返混（有返混？无返混？）的不同情况，可以建立八种（）不同的数学模型表表5 51 1催化反应器数学模型分类催化反应器数学模型分类一维模型一维模型A类类：拟均相模型拟均相模型B类类：非均相模型非均相模型AI I：基础模型基础模型基础模型基础模型BI I：基础模型相间分布基础模型相间分布基础模型相间分布基础模型相间分布AIIII：AI+I+轴向返混轴向返混轴向返混轴向返混BI

21、I II：BI+I+轴向返混轴向返混轴向返混轴向返混二维模型二维模型AIIIIII：AI I 径向分布径向分布径向分布径向分布BIII III：BI I径向分布径向分布径向分布径向分布AIVIV：AIII III 轴向返混轴向返混轴向返混轴向返混BIV IV：BIIIIII轴向返混轴向返混轴向返混轴向返混6.3.1等温反应器的计算等温反应器的计算 无径向床壁传热，可做平推流处理平推流处理，只在轴向上有温度和浓度的变化，而在与流向垂直的截面上（径向）是等温等浓度。6.3拟均相一维模型拟均相一维模型6.3拟均相一维模型拟均相一维模型6.3.2单层绝热床的计算单层绝热床的计算 6.3拟均相一维模型拟

22、均相一维模型6.3拟均相一维模型拟均相一维模型化学反应工程6.3.3多层绝热床的计算多层绝热床的计算（了解，自学）（了解，自学）多层绝热床每一层的计算方法与以前介绍的相同，只是在层间加以中间冷却（或加热）改变了温度和浓度，根据层间进行调节，进行物料和热量衡算。6.3拟均相一维模型拟均相一维模型化学反应工程6.3.4多层床的最优化问题多层床的最优化问题（了解，自学）（了解，自学）对于可逆放热反应，要使r尽可能保持最大，必须随转化率的增高，按温度曲线相应降低温度要使床层温度尽可能接近最优分布，以便使催化剂用量尽可能的少，就必须有尽可能多的层数，但层数越多，效果越微，很少超过四层。6.3拟均相一维模

23、型拟均相一维模型化学反应工程图图自己换热式反应器热平衡图自己换热式反应器热平衡图6.3.5 自己换热式反应器的设计方法（了解，自学）（了解，自学）6.3拟均相一维模型拟均相一维模型化学反应工程对于单一反应，以反应产物j为着眼组分，物料平衡式：床层内热量衡算：预热管内热量衡算：)()(21111TTUASHrdtdTncSLjBop-D-=)(2122TTUAdldTncSLop-=-气体总的摩尔流量单位床层截面为基准的-ondlrdZnjBjo)(=6.3拟均相一维模型拟均相一维模型化学反应工程6.4拟均相二维模型拟均相二维模型6.4.1模型方程模型方程（了解，自学）（了解，自学）绝大多数固定

24、床反应器呈圆柱形结构，空间变量分为径向和轴向，描述这两方向上的浓度和温度变化用偏微分方程。化学反应工程6.4.2模型方程的解法模型方程的解法（了解，自学）（了解，自学）（1）显示差分法（2）隐式差分法6.4拟均相二维模型拟均相二维模型6.2绝热气固反应器绝热气固反应器（了解，自学）（了解，自学）定义定义绝热反应器是一种和周围没有热交换的反应器。绝热反应器是一种和周围没有热交换的反应器。因此，如果放热反应是在绝热固定床反应器中进行，则气流的温度将沿着入口到出口的方向增加。另一方面，如果吸热反应是在同样的反应器中进行，则气流的温度将会沿着反应器的长度降低。一维塞流模型一维塞流模型 绝绝热热气气体体

25、反反应应器器最最简简单单的的模模型型是是一一维维塞塞流流模模型型,为为了了适适应应此模型，反应器必须符合下列条件此模型，反应器必须符合下列条件:反应器和周围环境不应有热交换反应器和周围环境不应有热交换(绝热条件绝热条件)。通通过过固固定定床床催催化化剂剂床床的的流流动动应应该该是是在在床床层层内内任任何何一一点点反反应物流的线速度都是相同的。应物流的线速度都是相同的。不应有轴向扩散。不应有轴向扩散。不应有径向扩散。不应有径向扩散。一维模型绝热固定床气体反应器的设计一维模型绝热固定床气体反应器的设计 Fyodx=rvdV=rvAcdz Fyodx(-H)=miCpidT式中,(-H)整个反应系统

26、的反应热；mi每一种反应物(包括反应物和产物)的摩尔流量，kmo1s；Cpi每一种反应物的克分子热容量，kmol(kmo1.K)。设计用于放热反应的固定床气体反应器设计用于放热反应的固定床气体反应器 对绝热反应器的限制是由限定的出口温度对绝热反应器的限制是由限定的出口温度t最大最大决定的。决定的。如果达到如果达到t最大最大标志着某些不理想的过程如副反应、选择标志着某些不理想的过程如副反应、选择性不良和催化剂严重结垢开始出现。如末达到所要求的性不良和催化剂严重结垢开始出现。如末达到所要求的转化率转化率x，温度就已上升到限制温度，此时，可采用几，温度就已上升到限制温度，此时，可采用几个方法中的一个

27、在温度不超过个方法中的一个在温度不超过t最大最大的情况下来提高的情况下来提高x使其使其达到所要求的值。达到所要求的值。方法方法1 (1)把催化剂体积分成两个或两个以上的床层，以便使把催化剂体积分成两个或两个以上的床层，以便使每一床层都达到每一床层都达到t最大最大；(2)用间接换热器冷却两相邻床层之间的气体反应物用间接换热器冷却两相邻床层之间的气体反应物使其回到第一床层入口温度使其回到第一床层入口温度to 方法方法2如果所用间接换热器反应温度太高的话，则可采用如果所用间接换热器反应温度太高的话，则可采用第二个方法，把一种冷的反应物喷入两相邻床层之间的第二个方法，把一种冷的反应物喷入两相邻床层之间

28、的反应物流中，把反应物冷却到反应物流中，把反应物冷却到To由于补充了一种被消由于补充了一种被消耗的反应物，这种直接冷却或急冷方法又有使反应朝着耗的反应物，这种直接冷却或急冷方法又有使反应朝着所要求的方向进行的优点。所要求的方向进行的优点。方法方法3 第三个方法是一种老的德国分段急冷的技术，虽然目第三个方法是一种老的德国分段急冷的技术，虽然目前仍在使用，但其缺点是波动很大，且难于设计。前仍在使用，但其缺点是波动很大，且难于设计。方法方法4第四个方法是在反应器进料中加入一种惰性气或一第四个方法是在反应器进料中加入一种惰性气或一种过量的反应物，其作用好象一个受热器一样，吸收大种过量的反应物，其作用好

29、象一个受热器一样，吸收大部分放出的热量，这样，对给定的转化率而言，就减小部分放出的热量，这样，对给定的转化率而言，就减小了温升。和其他方法相比，此法有一严重缺点，即为了了温升。和其他方法相比，此法有一严重缺点，即为了容纳惰性气或过量的反应物，需要增大反应器的尺寸。容纳惰性气或过量的反应物，需要增大反应器的尺寸。第四个方法派生出来的一个方法是加入一种惰性液第四个方法派生出来的一个方法是加入一种惰性液体，它在反应器的通常压力下，其沸点接近反应温度，体，它在反应器的通常压力下，其沸点接近反应温度，反应放出的热使惰性液体蒸发，这样，其所吸收的大量反应放出的热使惰性液体蒸发，这样，其所吸收的大量热量就满

30、足了潜热的需要。在某些情况下，和采用纯气热量就满足了潜热的需要。在某些情况下，和采用纯气相方没相比，这种方法反应速率比较低。此方法的缺点相方没相比，这种方法反应速率比较低。此方法的缺点是需要的反应器大而贵。是需要的反应器大而贵。如如果果没没有有一一个个方方法法能能妥妥善善地地解解决决放放出出的的热热量量问问题题，则则必必须考虑气体反应在下列反应器中的一种内进行：须考虑气体反应在下列反应器中的一种内进行：流化床反应器；流化床反应器；带有一种惰性液体的气带有一种惰性液体的气-液相反应器液相反应器非恒温非绝热固定床气体反应器非恒温非绝热固定床气体反应器 绝热固定床气体反应器设计实例：绝热固定床气体反

31、应器设计实例：甲苯加氢脱烷基反应器的设计甲苯加氢脱烷基反应器的设计 反应速率方程为：反应速率方程为：rv=-dCA/d=kCACB0.5=Aexp(-E/RGT)CACB0.5 式中式中,CA甲苯浓度，甲苯浓度，kmo1m3；CB氢浓度，氢浓度，kmo1m3；k甲苯加氢脱烷基的反应速率常数甲苯加氢脱烷基的反应速率常数;A频率因子；频率因子；E活化能活化能 148114kJkmo1；H-49974kJkmo1转化的甲苯；转化的甲苯；停留时间，停留时间，s。发表的活化能数值为发表的活化能数值为35.4kcalg.mol，但得不到，但得不到发表的频率固子发表的频率固子A的数值，该数值取决于具体的催化

32、剂的数值，该数值取决于具体的催化剂,现假定现假定A=573 x l06，反应器入口的反应温度假定为，反应器入口的反应温度假定为To600十十2732873.2K，反应器压力可从，反应器压力可从35到到70绝绝对大气压，假定约对大气压，假定约50绝对大气压或绝对大气压或50 bar。根据芳烃加。根据芳烃加氢的一般经验，氢与芳烃之比可假定为氢的一般经验，氢与芳烃之比可假定为5.0.氢：氢：Cp=20.786 甲烷：甲烷：Cp=0.04414T十十27.87 苯：苯：Cp=0.1067T十十103.18 甲苯：甲苯：Cp0.03535T十十124.85 进入加氢脱烷基装置的混合氢气物流的纯度为进入加

33、氢脱烷基装置的混合氢气物流的纯度为90，其余的为甲烷。，其余的为甲烷。每一反应物或产物的摩尔分数是：每一反应物或产物的摩尔分数是：甲苯：甲苯：Fyo(1-x)苯：苯：Fyox氢：氢：Fyo(5-x)+0.9急冷量急冷量甲烷：甲烷：Fyo(5/9+x)+0.1急冷量急冷量 摩尔分数：摩尔分数：Fyo(1+5+59)+急冷量，或急冷量，或6.5555Fyo+急冷量急冷量每一组分的摩尔分数为：每一组分的摩尔分数为：甲苯：甲苯：y1Fyo(1-x)(6.5555Fyo+急冷量急冷量)yA苯：苯：y2Fyox(6.5555Fy。十急冷量。十急冷量)氢：氢：y3(Fyo(5一一x)+0.9急冷量急冷量)(

34、6.5555Fyo+急冷量急冷量)=yB甲烷：甲烷：y4=(Fyo(59+x)+0.1急冷量急冷量)(6.5555Fyo+急冷量急冷量)浓度项可写成：浓度项可写成：C=N/V=Np/(NRG*T)=p/(RG*T)因为因为:p=Py 所以：所以：C=yP/(RG*T)CA=PFyo(1-x)(RG*T(6.5555Fyo+急冷量急冷量)CB=P(Fyo(5-x)+0.9急冷量急冷量)(RG*T(6.5555Fyo+急冷量急冷量)若把这些若把这些CA和和CB项代入式项代入式(6-32)，则：，则：rv=-dCA/d=kCACB0.5=Aexp(-E/RGT)(P(RG*T(6.5555Fyo+急

35、急冷量冷量)1.5 Fyo(1-x)(Fyo(5-x)+0.9急冷量急冷量)0.5将上式与（将上式与（6-8）合并：）合并：dx/dz=AcAexp(-E/RGT)(P/RG*T(6.5555Fyo+急急 冷冷 量量)1.5(1-x)(Fyo(5-x)+0.9急冷量急冷量)0.5 现在计算每一床层现在计算每一床层(除最后床层外除最后床层外)末端处的急冷量末端处的急冷量:急冷量急冷量=miCpi(T-To)/(Cpi急冷急冷(T-298.2)Cpi 急冷急冷=0.9x20.786+0.1x66.413=25.349 kJ/Kmol K 利用利用Ergun方程计算压降方程计算压降:在在Ergun方

36、程中，方程中，G为质量流率其单位为为质量流率其单位为kg(s m2)以空塔为基础的空塔流速以空塔为基础的空塔流速um的方程为的方程为:若催化剂的形状是直径和高度相等的圆柱体，令若催化剂的形状是直径和高度相等的圆柱体，令rp颗粒半径，则，颗粒半径，则，若给定若给定2rp3mm，则，则dp0.003m。反应混合物的粘度反应混合物的粘度可用以下的关系式从各个反应物的粘可用以下的关系式从各个反应物的粘度度1，2，3，4计算出来此方程是：计算出来此方程是：ij项可以是12、13、14、21、23、24、31、32、3441、42、43，Mi和Mj是各组分的分子量。假定内径为假定内径为2mdT/dz=(d

37、x/dz)Fyo(-H)/(miCpi)dx/dz=AcAexp(-/RGT)(P/RG*T(6.5555Fyo+急急冷冷量量)1.5(1-x)(Fyo(5-x)+0.9急冷量急冷量)0.5联联立立解解微微分分方方程程可可求求出出完完成成一一定定生生产产任任务务的的反反应应器器高高度度z。电电算算结结果果如如表表6-1所示：所示：6.4 6.4 多段间接换热式催化反应器优化设计多段间接换热式催化反应器优化设计（要求：了解即可）（要求：了解即可）确定了多段间接换热式催化反应器的段数、进口气体确定了多段间接换热式催化反应器的段数、进口气体组成、最终转化率、生产能力等。多段间接换热式绝热催组成、最终

38、转化率、生产能力等。多段间接换热式绝热催化反应器的设计是一个确定各段始末温度、转化率的最佳化反应器的设计是一个确定各段始末温度、转化率的最佳分布问题，最终结果是使催化剂的的总用量最少。分布问题，最终结果是使催化剂的的总用量最少。1催化床出口温度不受催化剂耐热温度限制时的优化设计催化床出口温度不受催化剂耐热温度限制时的优化设计目标函数目标函数V VRTRT具有极小值的必要条件如下：具有极小值的必要条件如下：式中，式中，xA，i i段的出口转化率，即段的出口转化率，即i+1段的入口转化率。段的入口转化率。Ti i 段的入口温度段的入口温度 整理后得：整理后得：又因为：又因为：上式（6-42）可变为

39、：6-24上一段出口处的反应速率与下一段入口处反应速率之绝对值相等上一段出口处的反应速率与下一段入口处反应速率之绝对值相等各段进口温度各段进口温度T Ti i下，下，又因为：又因为：所以：即，用图示的形式表示：用图示的形式表示：又可变为，图6-14 方程式（6-45）用图示的形式表示 当xA轴上方的阴影由积与下方的阴影面积相等时，即为符合式(6-45)的xAi值。设计思路：设计思路：(1)根据反应气体混合物的初始组成，催化剂动力学参数，作过程的平衡曲线与最佳温度曲线。a平衡曲线：根据r=0，可得 XA，eq =f（Teq）b最佳温度曲线：根据r/T=0，可得 XA=f（Topt）(2)假定第一

40、段入口温度为To1，T=To+（XA-XAO）其中，=NTOyAO（-HR）/NTCP求解出第一段出口转化率与出口温度XA1，T1。（3 3）根据式）根据式 即，第一段出口速率等于第二段入口速率，可求出第二段即，第一段出口速率等于第二段入口速率，可求出第二段入口温度入口温度TO2。(4)重复（2）步求出第二段出口温度与出口转化率 T2，XA2，及重复（3）步可求出第三段的入口温度To3。(5)检查算出的最后一段的的出口转化率与生产任务是否一致，若不一致，重新设定TO1重复（2）-（5）步，直至直至算出最后一段的的出口转化率与生产任务一致(6)检查计算结果与平衡线及最佳温度线的关系。2.2.催化

41、剂有起始活性温度限制催化剂有起始活性温度限制设计思路设计思路:(1)据据反反应应气气体体混混合合物物的的初初始始组组成成，催催化化剂剂动动力力学学参参数数，作作过过程程的的平平衡衡曲曲线线与最佳温度曲线。与最佳温度曲线。a平衡曲线：平衡曲线：根据根据r=0，可得，可得 XA，eq =f（Teq）b最佳温度曲线：根据最佳温度曲线：根据 r/T=0，可得可得 XA=f（Topt）（2）设定第一段入口温度为）设定第一段入口温度为To1，To1值为催化剂活性温度，是一固定值。值为催化剂活性温度，是一固定值。（3）设第一段出口转化率为）设第一段出口转化率为XA1，（4）根据绝热温升线可计算第一段的出口温

42、度）根据绝热温升线可计算第一段的出口温度T1。T=To+（XA-XAO）（6-47）其中，其中，=NTOyAO（-HR）/NT2CP(5)(5)由由即，第一段出口速率等于第二段入口速率，可求出第二段即，第一段出口速率等于第二段入口速率，可求出第二段入口温度入口温度TO2。（6）根据：式（）根据：式（6-46）和绝热温升线（和绝热温升线（6-47）（7）重重复复（5）、（6）步步算算出出最最后后一一段段的的XAn，检检查查算算出出的的最最后后一一段段的的的的出出口口转转化化率率与与生生产产任任务务是是否否一一致致，若若不不一一致致，重重新新设设定定XA1重重复复（3）-（6）步步，直直至至直直至

43、至算算出出最最后后一一段段的的的出口转化率与生产任务一致。的出口转化率与生产任务一致。（8）检查计算结果与平衡线及最佳温度线的关系。）检查计算结果与平衡线及最佳温度线的关系。T=To+（XA2-XA1）（6-47）其中，=NTOyAO（-HR）/NTCP求解出第二段出口转化率与出口温度XA2，T2。3催化剂有最高耐热温度催化剂有最高耐热温度Tmax限制限制设计思路：设计思路：（1）反应气体混合物的初始组成，催化剂动力学参数，作过程的平衡曲线与最佳温度曲线。（2）假定第一段入口温度为）假定第一段入口温度为To1，根据：式（根据：式（6-46）和绝热温升线（和绝热温升线（6-47）a平衡曲线：根据

44、r=0，可得 XA，eq =f（Teq）b最佳温度曲线：根据r/T=0，可得 XA=f（Topt）T=To+（XA-XAO）（6-47）其中，=NTOyAO（-HR）/NTCP求解出第一段出口转化率与出口温度XA1，T1。若T1超过Tmax，则以Tmax为T1，并由式（6-47）算出XA1。亦可由 与与 T=To+T=To+（X XA A-X-XAOAO）（6-476-47）计算出计算出X XA1A1。（3）根据式（）根据式（6-43）即，第一段出口速率等于第二段入口速率，可求出第二段入口温度TO2。（4）重复（2）步求出第二段出口温度与出口转化率 T2，XA2，及重复（3）步可求出第三段的入

45、口温度To3。计算过程中，若Ti 出口温度超过Tmax应以Tmax为出口温度Ti。并用式（6-47）或式（6-47）结合式（6-48）计算出口转化率XAi。（5）检查算出的最后一段的的出口转化率与生产任务是否一致，若不一致，重新设定TO1重复（2）-（5）步，直至算出最后一段的的出口转化率与生产任务一致。（6）检查计算结果与平衡线及最佳温度线的关系。4 4例题例题 某日产干吨两段间接换热式径向氨合成塔某日产干吨两段间接换热式径向氨合成塔合成氨的速率：合成氨的速率：R=dyNH3/do=kT（1+yNH3）2Kf2fH21.5fN2/fNH3 fNH3/fH21.5解：（解：（1）不受催化剂温度

46、限制的影响时，）不受催化剂温度限制的影响时，a.TO1=340,初始组成已经确定，由下式，初始组成已经确定，由下式，T=To+（XA-XAO）（6-47）其中，=NTOyAO（-HR）/NTCP求解出第一段出口转化率与出口温度XA1，T1。T1=516.4,yNH3,1=0.1487,VR1=11.7367 m3 b.b.由由 即，第一段出口速率等于第二段入口速率，可求出第二段即，第一段出口速率等于第二段入口速率，可求出第二段入口温度入口温度TO2=371.7。c.再由再由 T=To+（XA2-XA1）（6-47）其中，=NTOyAO（-HR）/NTCP求解出第二段出口转化率与出口温度XA2，

47、T2。T2=461.4,yNH3,2=0.2146,VR2=28.2155 m3（2 2）受催化剂温度限制的影响时，）受催化剂温度限制的影响时，a.a.第一段出口温度第一段出口温度T T1 1计算值超过了催化剂耐热温度计算值超过了催化剂耐热温度Tmax,TTmax,T1 1=Tmax=Tmax=510,510,并由式（并由式（6-476-47）算出）算出X XA1A1。亦可由。亦可由 与与 T=To+（XA-XAO）计算出计算出XA1。T1=510,yNH3,1=0.1432,VR1=10.2508 m3b.由由 即即，第第一一段段出出口口速速率率等等于于第第二二段段入入口口速速率率，可可求求

48、出出第第二二段段入口温度入口温度TO2=368。c.再由再由 T=To+（XA2-XA1）（6-47）其中，=NTOyAO（-HR）/NTCP求解出第二段出口转化率与出口温度XA2，T2。T2=465.3,yNH3,2=0.2144,VR2=30.6867 m3 6.5 连续换热式固定床催化反应器优化设计连续换热式固定床催化反应器优化设计（要求：了解即可）（要求：了解即可）连连续续换换热热式式催催化化床床可可分分为为催催化化剂剂装装在在管管间间的的内内冷冷自自热热式式、催催化化剂剂装装在在列列管管内内的的外外冷冷列列管管式式和和外外部部供供热热管管式式。本本节节主主要要以以三三套套管管并并流流

49、式式合合成成氨氨催催化化反反应应器器为为例例讨讨论论连连续续换换热热式式固固定定床床催催化化反反应器优化设计问题应器优化设计问题.图6-15 三套管并流式冷管催化床温度分布及操作情况 1.内冷自热式氨合成反应器的数学模型内冷自热式氨合成反应器的数学模型1.1.绝热段绝热段:(1)物料衡算物料衡算氨合成为变摩尔数反应氨合成为变摩尔数反应,合成氨反应的动力学方程为:将式将式(6-49)(6-49)、(6-50)(6-50)代入式（代入式（6-516-51）整理后得：）整理后得：其中，其中，再因为，Wo为校标准状态下氨分解基气体体积流率及床层截面积计算的虚拟线速度，m3(m2s)。对于绝热层，应按绝

50、热层床层截面积来计算Wo，以Woh表示，再计入活性校正系数COR及表征催化剂使用时间对宏观反应速率影响的寿命因子TF，最后可得绝热层中氨摩尔分率随床高l变化的表达式，式（6-52）变化为：(2)热量衡算热量衡算 T Tbobo（设定）、（设定）、y yioio已知时，联立式（已知时，联立式（6-546-54）与式（）与式（6-566-56）可计算组成可计算组成y yi i、床层温度、床层温度T Tb b随床层高度的变化。随床层高度的变化。1.2 冷却段冷却段(1)物料衡算物料衡算冷却段的温度分布与流冷却段的温度分布与流体流向如图体流向如图6-17所示所示图6-17 冷却段的温度分布与流体流向（