第五章气固反应优秀PPT.ppt

第五章气固反应第五章气固反应第一页，本课件共有108页内内 容容第一节第一节 固定床气固相催化反应器的基本类型和数学模型固定床气固相催化反应器的基本类型和数学模型第二节第二节 固定床流体力学固定床流体力学第三节第三节 固定床热量与质量传递过程固定床热量与质量传递过程第四节第四节 绝热式固定床反应器绝热式固定床反应器第五节第五节 连续换热内冷自热式催化反应器连续换热内冷自热式催化反应器第六节第六节 连续换热外冷及外热管式催化反应器连续换热外冷及外热管式催化反应器第七节第七节 薄床层催化反应器薄床层催化反应器第二页，本课件共有108页第一节第一节固定床气固相催化反应器的基本类型和数学模型固定床气固相催化反应器的基本类型和数学模型1.固定床反应器：固定床反应器：凡是流体通过不动的固体物料所形成的床层而进行反应的装置。2.固定床反应器优缺点：固定床反应器优缺点：催化剂不易磨损；床层内流体平推流，较少量催化剂可获较大生产能力；传热较差。第三页，本课件共有108页3.固定床反应器的分类固定床反应器的分类固定床反应器不同的传热要求和传热方式绝热式二段三段四段冷激式反应特征单段绝热式段间反应气冷却或加热方式中间间接换热式多段绝热式原料气冷激式非原料气冷激式换热式对外换热式自热式加压热水（240）导热油（250300）熔盐（300）反应气的流动方向轴向流动固定床反应器径向流动固定床反应器第四页，本课件共有108页4.数字模型数字模型描述固定床反应器的数字模型分为：拟均相和非均相两类拟均相模型忽略床层中粒子与流体间温度与浓度的差别。非均相模型考虑床层中粒子与流体间温度与浓度的差别。第一节第一节固定床气固相催化反应器的基本类型和数学模型固定床气固相催化反应器的基本类型和数学模型第五页，本课件共有108页1.固定床的物理特性固定床的物理特性(1)颗粒密度p床层密度或堆积密度B(2)当量直径及形状系数 具有相同体积的球粒子直径dv 来表示dv？(1)具有相同外表面积球粒子直径Dp来表示Dp？1.具有相同比外表面积球粒子直径ds来表示ds？形状系数=Ss/Sp 第二节第二节固定床流体力学固定床流体力学第六页，本课件共有108页(3)混合颗粒平均直径(4)固定床当量直径及空隙率 de=4RH=4/Se=2/3*(/1-)*ds(5)空隙率及径向流速 了解即可第二节第二节固定床流体力学固定床流体力学dp=1/n（i1=xdii)第七页，本课件共有108页2.单相流体在固定床中的流动及压降单相流体在固定床中的流动及压降(1)流动(了解即可)(2)压降 第二节第二节固定床流体力学固定床流体力学(3)压降影响因素(了解即可)修正雷诺数：修正雷诺数：um 平均流速(空床气速)；L 床层高度；、流体的密度和粘度；ds 比表面当量直径。第八页，本课件共有108页影响固定床压力降的因素影响固定床压力降的因素流体流体流体的密度流体的密度流体的粘度流体的粘度流体的质量流率流体的质量流率床层床层床层的高度床层的高度床层的空隙率床层的空隙率流通截面积流通截面积颗粒颗粒颗粒的形状颗粒的形状颗粒的粒度颗粒的粒度颗粒的表面粗糙度颗粒的表面粗糙度颗粒的物理特性颗粒的物理特性第九页，本课件共有108页床层压降计算例：在内径为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.15第十页，本课件共有108页第十一页，本课件共有108页第十二页，本课件共有108页1.固定床中的传热固定床中的传热传热包括：粒内传热，颗粒与流体间的传热(给热系数hp p)，床层与器壁的传热(给热系数hw w).以床层的平均温度 与管壁温差为推动力-总给热系数 h hc有效导热系数e轴向轴向es径向有效导热系数 er=e0e0+(er)t第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程传热方式：导热、对流传热、热辐射。传热途径：粒内传热、颗粒与流体间传热、床层与器壁间传热。第十三页，本课件共有108页固定床的径向传热方式床层空隙内部流体的传热颗粒之间通过接触的传热颗粒表面附近流体中的传热空隙内部流体的辐射传热流体混合所引起的径向对流传热颗粒表面之间的热辐射传热第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第十四页，本课件共有108页一、颗粒与流体间传热系数一、颗粒与流体间传热系数1、传热因子适用范围颗粒与流体间的传（给）热系数hpdp可用dv代注意：在参考其他教材时，给热系数多以表示。第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第十五页，本课件共有108页计算传热系数的经验公式有很多，可从有关文献或工具书中查取。应用时要注意公式规定的特性尺寸、特性温度和适用范围。第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第十六页，本课件共有108页式中：单位重量催化剂的外表面积 床层比表面积Se的校正系数主体、外表面温度2、流体与颗粒间传热温差的计算、流体与颗粒间传热温差的计算热量平衡热量平衡第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第十七页，本课件共有108页上式可整理成上式可整理成其中，传热量其中，传热量普兰特数普兰特数根据6-12的关联图，查图可求得不同条件下的t。第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第十八页，本课件共有108页图图6-12固定床中流体与颗粒外表面温差固定床中流体与颗粒外表面温差t0rAHrtG时时t第十九页，本课件共有108页e 是针对拟均相模型提出的综合性传热参数，一般是指径向热导率er二、固定床的有效热导率二、固定床的有效热导率有效导热系数e e轴向有效导热系数es径向有效导热系数 er=e0+(er)t t第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第二十页，本课件共有108页式中，式中，流体的热导率流体的热导率流体静止时床层的热导率流体静止时床层的热导率径向与轴向传质速率之比径向与轴向传质速率之比颗粒间距与粒径比的影响颗粒间距与粒径比的影响值可由图值可由图6-14查取。查取。粒径粒径/管径管径图图求有效热导率求有效热导率e第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第二十一页，本课件共有108页包含床层包含床层空隙和颗粒空隙和颗粒对传热的贡献，由下式计算对传热的贡献，由下式计算式中，分别表示颗粒与流体的热导率；床层的平均温度；空隙的辐射给热系数；颗粒的辐射给热系数；颗粒接触点处流体薄膜对导热的影响。第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第二十二页，本课件共有108页床层空隙的辐射给热系数：hrv式中，粒子表面的热辐射率；床层的平均温度；空隙的辐射给热系数；第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第二十三页，本课件共有108页颗粒的辐射给热系数：颗粒接触点处流体薄膜对导热的影响。由图由图6-15查取。查取。式中，粒子表面的热辐射率；床层的平均温度；颗粒辐射给热系数；第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第二十四页，本课件共有108页由图由图6-15查取。查取。图图6-15求所用的值求所用的值第二十五页，本课件共有108页当颗粒直径甚小，床层温度不是很高，以及含有液体时，空隙和颗粒的辐射传热可忽略，床层空隙和颗粒对传热的贡献可简化为：第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第二十六页，本课件共有108页小结：固定床的优缺点固定床的分类基本概念 颗粒直径的表示方法 颗粒的形状系数 混合颗粒的平均直径 床层空隙率 床层比表面积 床层当量直径固定床的径向流速分布床层压降厄根(Ergun)方程修正雷诺数影响固定床压力降的因素固定床中的传热颗粒与流体间传热系数第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第二十七页，本课件共有108页三、床层与器壁间的给热系数三、床层与器壁间的给热系数h0 一维模型中，床层与器壁间传热速率为 h0可由经验公式计算式中，y 无量纲数（适用范围：y 0.2）：床层平均温度;：器壁温度;第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第二十八页，本课件共有108页b 无量纲数均为无量纲的b的函数hW 壁面处对流给热系数；第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第二十九页，本课件共有108页适用范围：适用范围：由图由图6-16查得。查得。第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第三十页，本课件共有108页图图6-16求求h0例如：当例如：当b12.98时，时，分别对应分别对应5.02和和0.25b第三十一页，本课件共有108页由于催化剂表面存在滞流边界层，气流主体浓度与催化剂颗粒表面浓度存在差异。在滞流层内有浓度差，必然存在扩散。一、颗粒与流体间的传质系数一、颗粒与流体间的传质系数 单位体积（或质量）催化剂上关键组分A的传质速率第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程距离0CACAgCASCACRp第三十二页，本课件共有108页整个传质方程的核心，总包了各种条件对传质的影响。由实验关联式计算。关联式之一：传质因子：传质因子：JD实验关联式：实验关联式：气相密度：气相密度;：气体质量流率：气体质量流率;：气相粘度：气相粘度;：气相分子扩散系数：气相分子扩散系数;JD是雷诺数的函数是雷诺数的函数第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程Sc为施密特准数为施密特准数第三十三页，本课件共有108页um：平均流速(空床气速)；g、g：气体的密度和粘度；ds：比表面当量直径。当修正雷诺数在：Sc为施密特准数为施密特准数第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程比表面当量直径：单颗催化剂的外表面积折合成直径为dS的球形颗粒应有的外表面积。第三十四页，本课件共有108页传质对反应的影响达姆科勒数 第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程为外部效率因子，或称外部有效因子，它是D的函数，反映界面反应物浓度降低对反应速率的影响程度第三十五页，本课件共有108页Da为反应速率与扩散速率的比值，反映了体系中外扩散的影响程度。数值越大，或反应速率越快，外扩散的影响就越大。对一级反应第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第三十六页，本课件共有108页第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第三十七页，本课件共有108页特殊情况：反应速率常数k比传质系数kg大得多，则颗粒外表面处A的浓度为零，属外扩散控制。反应速率常数k比传质系数kg小得多，则颗粒外表面处A的浓度与气相主体浓度相等，属内扩散或动力学控制，外扩散可不予考虑。第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第三十八页，本课件共有108页传热因子传热因子传质因子传质因子：气相密度 ：气体质量流率 ：气相粘度 ：气相分子扩散系数 JD是雷诺数的函数施密特准数第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第三十九页，本课件共有108页由传质和传热的类比原理有由传质和传热的类比原理有上式在缺少数据时用来推算数据十分有用上式在缺少数据时用来推算数据十分有用.无论是传质或是传热系数无论是传质或是传热系数,增加质量速率增加质量速率G都可以加快流体都可以加快流体和催化剂外表面间的传递速率和催化剂外表面间的传递速率,但相应的床层压降也增加但相应的床层压降也增加.外扩散速率外扩散速率,内扩散速率内扩散速率,表面反应速率的相对大小是决定反表面反应速率的相对大小是决定反应控制步骤的关键应控制步骤的关键.第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第四十页，本课件共有108页外扩散过程对表面温度的影响外扩散过程对表面温度的影响由传质计算可得：由传质计算可得：第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第四十一页，本课件共有108页两式相等可得：两式相等可得：第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第四十二页，本课件共有108页上式将流固相的温度差与浓度差联系起来了。上式将流固相的温度差与浓度差联系起来了。进一步简化，前面有：进一步简化，前面有：极为相似。极为相似。JHJD相除相除第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第四十三页，本课件共有108页第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第四十四页，本课件共有108页据实验得：据实验得：第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第四十五页，本课件共有108页例：苯加氢反应器在例：苯加氢反应器在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)。估算该处催化剂的外表面浓度。估算该处催化剂的外表面浓度。已知：气体粘度已知：气体粘度=1.410-4g.cm-1s-1，综合扩散系数，综合扩散系数D=0.267cm2s-1。第三节第三节固定床热量与质量传递过程固定床热量与质量传递过程第四十六页，本课件共有108页解：解：（1）计算催化剂的粒径）计算催化剂的粒径dS。第四十七页，本课件共有108页（2）计算床层中气体的修正雷诺数。）计算床层中气体的修正雷诺数。第四十八页，本课件共有108页（3）计算）计算JD和和kg值。值。第四十九页，本课件共有108页（4）计算）计算cAG和和cAS第五十页，本课件共有108页习题：试计算例题中催化剂的外表面处温度。习题：试计算例题中催化剂的外表面处温度。已已知知反反应应热热为为(-H)=2.135105Jmol-1，气气体体的的定定压压比比热热容容cP=49Jmol-1K-1。第五十一页，本课件共有108页小结：颗粒与流体间的传（给）热系数hp传热因子JH流体与颗粒间传热温差的计算固定床的有效热导率床层与器壁间的给热系数 h0传质因子JD第五十二页，本课件共有108页6-3 6-3 催化反应器的数学模型催化反应器的数学模型1 1，非均相拟均相，非均相拟均相2 2，一维模型二维模型，一维模型二维模型3 3，理想流动非理想流动，理想流动非理想流动拟均相适用情况：拟均相适用情况：1 1，化学动力学控制，化学动力学控制2 2，活性较正系数（无宏观动力学资料），活性较正系数（无宏观动力学资料）一维二维：轴向浓度差、温度差；轴径向浓度差、温度差一维二维：轴向浓度差、温度差；轴径向浓度差、温度差理想流动：不考虑返混（理想流动：不考虑返混（PFR）；非理想流动：考虑返混（扩散）非理想流动：考虑返混（扩散）第五十三页，本课件共有108页6-3 6-3 催化反应器的数学模型催化反应器的数学模型维数维数相数相数返混返混根据相态(拟均相？非均相？），维数（一维？二维？），返混（有返混？无返混？）的不同情况，可以建立八种（）不同的数学模型第五十四页，本课件共有108页表表5 51 1催化反应器数学模型分类催化反应器数学模型分类一维模型一维模型A类类：拟均相模型拟均相模型B类类：非均相模型非均相模型AI：基础模型基础模型BI：基础模型相间分布基础模型相间分布AII：AI+轴向返混轴向返混BII：BI+轴向返混轴向返混二维模型二维模型AIII：AI 径向分布径向分布BIII：BI径向分布径向分布AIV：AIII 轴向返混轴向返混BIV：BIII轴向返混轴向返混第五十五页，本课件共有108页6.3.1等温反应器的计算等温反应器的计算 无径向床壁传热，可做平推流处理平推流处理，只在轴向上有温度和浓度的变化，而在与流向垂直的截面上（径向）是等温等浓度。6.3拟均相一维模型拟均相一维模型第五十六页，本课件共有108页6.3拟均相一维模型拟均相一维模型6.3.2单层绝热床的计算单层绝热床的计算 第五十七页，本课件共有108页6.3拟均相一维模型拟均相一维模型第五十八页，本课件共有108页6.3拟均相一维模型拟均相一维模型第五十九页，本课件共有108页化学反应工程6.3.3多层绝热床的计算多层绝热床的计算（了解，自学）（了解，自学）多层绝热床每一层的计算方法与以前介绍的相同，只是在层间加以中间冷却（或加热）改变了温度和浓度，根据层间进行调节，进行物料和热量衡算。6.3拟均相一维模型拟均相一维模型第六十页，本课件共有108页化学反应工程6.3.4多层床的最优化问题多层床的最优化问题（了解，自学）（了解，自学）对于可逆放热反应，要使r尽可能保持最大，必须随转化率的增高，按温度曲线相应降低温度要使床层温度尽可能接近最优分布，以便使催化剂用量尽可能的少，就必须有尽可能多的层数，但层数越多，效果越微，很少超过四层。6.3拟均相一维模型拟均相一维模型第六十一页，本课件共有108页化学反应工程图图自己换热式反应器热平衡图自己换热式反应器热平衡图6.3.5 自己换热式反应器的设计方法（了解，自学）（了解，自学）6.3拟均相一维模型拟均相一维模型第六十二页，本课件共有108页化学反应工程对于单一反应，以反应产物j为着眼组分，物料平衡式：床层内热量衡算：预热管内热量衡算：)()(21111TTUASHrdtdTncSLjBop-D-=)(2122TTUAdldTncSLop-=-气体总的摩尔流量单位床层截面为基准的-ondlrdZnjBjo)(=6.3拟均相一维模型拟均相一维模型第六十三页，本课件共有108页化学反应工程6.4拟均相二维模型拟均相二维模型6.4.1模型方程模型方程（了解，自学）（了解，自学）绝大多数固定床反应器呈圆柱形结构，空间变量分为径向和轴向，描述这两方向上的浓度和温度变化用偏微分方程。第六十四页，本课件共有108页化学反应工程6.4.2模型方程的解法模型方程的解法（了解，自学）（了解，自学）（1）显示差分法（2）隐式差分法6.4拟均相二维模型拟均相二维模型第六十五页，本课件共有108页6.2绝热气固反应器绝热气固反应器（了解，自学）（了解，自学）定义定义绝热反应器是一种和周围没有热交换的反应器。绝热反应器是一种和周围没有热交换的反应器。因此，如果放热反应是在绝热固定床反应器中进行，则气流的温度将沿着入口到出口的方向增加。另一方面，如果吸热反应是在同样的反应器中进行，则气流的温度将会沿着反应器的长度降低。第六十六页，本课件共有108页一维塞流模型一维塞流模型 绝绝热热气气体体反反应应器器最最简简单单的的模模型型是是一一维维塞塞流流模模型型,为为了了适适应应此此模型，反应器必须符合下列条件模型，反应器必须符合下列条件:反应器和周围环境不应有热交换反应器和周围环境不应有热交换(绝热条件绝热条件)。通通过过固固定定床床催催化化剂剂床床的的流流动动应应该该是是在在床床层层内内任任何何一一点点反反应应物物流的线速度都是相同的。流的线速度都是相同的。不应有轴向扩散。不应有轴向扩散。不应有径向扩散。不应有径向扩散。第六十七页，本课件共有108页一维模型绝热固定床气体反应器的设计一维模型绝热固定床气体反应器的设计 Fyodx=rvdV=rvAcdzFyodx(-H)=miCpidT式中,(-H)整个反应系统的反应热；mi每一种反应物(包括反应物和产物)的摩尔流量，kmo1s；Cpi每一种反应物的克分子热容量，kmol(kmo1.K)。第六十八页，本课件共有108页设计用于放热反应的固定床气体反应器设计用于放热反应的固定床气体反应器 对绝热反应器的限制是由限定的出口温度对绝热反应器的限制是由限定的出口温度t最大最大决定的。决定的。如果达到如果达到t最大最大标志着某些不理想的过程如副反应、选择性不标志着某些不理想的过程如副反应、选择性不良和催化剂严重结垢开始出现。如末达到所要求的转化率良和催化剂严重结垢开始出现。如末达到所要求的转化率x，温度就已上升到限制温度，此时，可采用几个方法中的，温度就已上升到限制温度，此时，可采用几个方法中的一个在温度不超过一个在温度不超过t最大最大的情况下来提高的情况下来提高x使其达到所要求的使其达到所要求的值。值。方法方法1(1)把催化剂体积分成两个或两个以上的床层，以便使把催化剂体积分成两个或两个以上的床层，以便使每一床层都达到每一床层都达到t最大最大；(2)用间接换热器冷却两相邻床层之间的气体反应物使其回用间接换热器冷却两相邻床层之间的气体反应物使其回到第一床层入口温度到第一床层入口温度to第六十九页，本课件共有108页方法方法2如果所用间接换热器反应温度太高的话，则可采用第二如果所用间接换热器反应温度太高的话，则可采用第二个方法，把一种冷的反应物喷入两相邻床层之间的反应物个方法，把一种冷的反应物喷入两相邻床层之间的反应物流中，把反应物冷却到流中，把反应物冷却到To由于补充了一种被消耗的反应物，由于补充了一种被消耗的反应物，这种直接冷却或急冷方法又有使反应朝着所要求的方向进行这种直接冷却或急冷方法又有使反应朝着所要求的方向进行的优点。的优点。第七十页，本课件共有108页方法方法3第三个方法是一种老的德国分段急冷的技术，虽然目前仍第三个方法是一种老的德国分段急冷的技术，虽然目前仍在使用，但其缺点是波动很大，且难于设计。在使用，但其缺点是波动很大，且难于设计。第七十一页，本课件共有108页方法方法4第四个方法是在反应器进料中加入一种惰性气或一种过量第四个方法是在反应器进料中加入一种惰性气或一种过量的反应物，其作用好象一个受热器一样，吸收大部分放出的的反应物，其作用好象一个受热器一样，吸收大部分放出的热量，这样，对给定的转化率而言，就减小了温升。和其他热量，这样，对给定的转化率而言，就减小了温升。和其他方法相比，此法有一严重缺点，即为了容纳惰性气或过量的方法相比，此法有一严重缺点，即为了容纳惰性气或过量的反应物，需要增大反应器的尺寸。反应物，需要增大反应器的尺寸。第四个方法派生出来的一个方法是加入一种惰性液体，第四个方法派生出来的一个方法是加入一种惰性液体，它在反应器的通常压力下，其沸点接近反应温度，反应它在反应器的通常压力下，其沸点接近反应温度，反应放出的热使惰性液体蒸发，这样，其所吸收的大量热量放出的热使惰性液体蒸发，这样，其所吸收的大量热量就满足了潜热的需要。在某些情况下，和采用纯气相方就满足了潜热的需要。在某些情况下，和采用纯气相方没相比，这种方法反应速率比较低。此方法的缺点是需没相比，这种方法反应速率比较低。此方法的缺点是需要的反应器大而贵。要的反应器大而贵。第七十二页，本课件共有108页如如果果没没有有一一个个方方法法能能妥妥善善地地解解决决放放出出的的热热量量问问题题，则则必必须须考虑气体反应在下列反应器中的一种内进行：考虑气体反应在下列反应器中的一种内进行：流化床反应器；流化床反应器；带有一种惰性液体的气带有一种惰性液体的气-液相反应器液相反应器非恒温非绝热固定床气体反应器非恒温非绝热固定床气体反应器第七十三页，本课件共有108页绝热固定床气体反应器设计实例：绝热固定床气体反应器设计实例：甲苯加氢脱烷基反应器的设计甲苯加氢脱烷基反应器的设计 反应速率方程为：反应速率方程为：rv=-dCA/d=kCACB0.5=Aexp(-E/RGT)CACB0.5式中式中,CA甲苯浓度，甲苯浓度，kmo1m3；CB氢浓度，氢浓度，kmo1m3；k甲苯加氢脱烷基的反应速率常数甲苯加氢脱烷基的反应速率常数;A频率因子；频率因子；E活化能活化能148114kJkmo1；H-49974kJkmo1转化的甲苯；转化的甲苯；停留时间，停留时间，s。第七十四页，本课件共有108页发表的活化能数值为发表的活化能数值为35.4kcalg.mol，但得不到发表，但得不到发表的频率固子的频率固子A的数值，该数值取决于具体的催化剂的数值，该数值取决于具体的催化剂,现假定现假定A=573xl06，反应器入口的反应温度假定为，反应器入口的反应温度假定为To600十十2732873.2K，反应器压力可从，反应器压力可从35到到70绝对大气压，绝对大气压，假定约假定约50绝对大气压或绝对大气压或50bar。根据芳烃加氢的一般经。根据芳烃加氢的一般经验，氢与芳烃之比可假定为验，氢与芳烃之比可假定为5.0.氢：氢：Cp=20.786甲烷：甲烷：Cp=0.04414T十十27.87苯：苯：Cp=0.1067T十十103.18甲苯：甲苯：Cp0.03535T十十124.85进入加氢脱烷基装置的混合氢气物流的纯度为进入加氢脱烷基装置的混合氢气物流的纯度为90，其，其余的为甲烷。余的为甲烷。第七十五页，本课件共有108页每一反应物或产物的摩尔分数是：每一反应物或产物的摩尔分数是：甲苯：甲苯：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急冷量急冷量)(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)第七十六页，本课件共有108页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+急急冷量冷量)1.5Fyo(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第七十七页，本课件共有108页现在计算每一床层现在计算每一床层(除最后床层外除最后床层外)末端处的急冷量末端处的急冷量:急冷量急冷量=miCpi(T-To)/(Cpi急冷急冷(T-298.2)Cpi急冷急冷=0.9x20.786+0.1x66.413=25.349kJ/KmolK利用利用Ergun方程计算压降方程计算压降:在在Ergun方程中，方程中，G为质量流率其单位为为质量流率其单位为kg(sm2)以空塔为基础的空塔流速以空塔为基础的空塔流速um的方程为的方程为:第七十八页，本课件共有108页若催化剂的形状是直径和高度相等的圆柱体，令若催化剂的形状是直径和高度相等的圆柱体，令rp颗粒半径，则，颗粒半径，则，若给定若给定2rp3mm，则，则dp0.003m。反应混合物的粘度反应混合物的粘度可用以下的关系式从各个反应物的粘度可用以下的关系式从各个反应物的粘度1，2，3，4计算出来此方程是：计算出来此方程是：第七十九页，本课件共有108页ij项可以是12、13、14、21、23、24、31、32、3441、42、43，Mi和Mj是各组分的分子量。第八十页，本课件共有108页假定内径为假定内径为2mdT/dz=(dx/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所所示：示：第八十一页，本课件共有108页第八十二页，本课件共有108页6.4 6.4 多段间接换热式催化反应器优化设计多段间接换热式催化反应器优化设计（要求：了解即可）（要求：了解即可）确定了多段间接换热式催化反应器的段数、进口气体组成、确定了多段间接换热式催化反应器的段数、进口气体组成、最终转化率、生产能力等。多段间接换热式绝热催化反应器最终转化率、生产能力等。多段间接换热式绝热催化反应器的设计是一个确定各段始末温度、转化率的最佳分布问题，的设计是一个确定各段始末温度、转化率的最佳分布问题，最终结果是使催化剂的的总用量最少。最终结果是使催化剂的的总用量最少。1催化床出口温度不受催化剂耐热温度限制时的优化设计催化床出口温度不受催化剂耐热温度限制时的优化设计第八十三页，本课件共有108页目标函数目标函数V VRTRT具有极小值的必要条件如下：具有极小值的必要条件如下：式中，式中，xA，ii段的出口转化率，即段的出口转化率，即i+1段的入口转化率。段的入口转化率。Tii段的入口温度段的入口温度整理后得：整理后得：第八十四页，本课件共有108页又因为：又因为：上式（6-42）可变为：6-24上一段出口处的反应速率与下一段入口处反应速率之绝对值相等上一段出口处的反应速率与下一段入口处反应速率之绝对值相等第八十五页，本课件共有108页各段进口温度各段进口温度T Ti i下，下，又因为：又因为：所以：即，第八十六页，本课件共有108页用图示的形式表示：用图示的形式表示：又可变为，图6-14 方程式（6-45）用图示的形式表示 当xA轴上方的阴影由积与下方的阴影面积相等时，即为符合式(6-45)的xAi值。第八十七页，本课件共有108页设计思路：设计思路：(1)根据反应气体混合物的初始组成，催化剂动力学参数，作过程的平衡曲线与最佳温度曲线。a平衡曲线：根据r=0，可得 XA，eq =f（Teq）b最佳温度曲线：根据r/T=0，可得 XA=f（Topt）(2)假定第一段入口温度为To1，T=To+（XA-XAO）其中，=NTOyAO（-HR）/NTCP求解出第一段出口转化率与出口温度XA1，T1。第八十八页，本课件共有108页（3 3）根据式）根据式 即，第一段出口速率等于第二段入口速率，可求出第二段即，第一段出口速率等于第二段入口速率，可求出第二段入口温度入口温度TO2。(4)重复（2）步求出第二段出口温度与出口转化率 T2，XA2，及重复（3）步可求出第三段的入口温度To3。(5)检查算出的最后一段的的出口转化率与生产任务是否一致，若不一致，重新设定TO1重复（2）-（5）步，直至直至算出最后一段的的出口转化率与生产任务一致(6)检查计算结果与平衡线及最佳温度线的关系。第八十九页，本课件共有108页2.2.催化剂有起始活性温度限制催化剂有起始活性温度限制设计思路设计思路:(1)据据反反应应气气体体混混合合物物的的初初始始组组成成，催催化化剂剂动动力力学学参参数数，作作过过程程的的平平衡衡曲曲线线与最佳温度曲线。与最佳温度曲线。a平衡曲线：平衡曲线：根据根据r=0，可得，可得XA，eq=f（Teq）b最佳温度曲线：根据最佳温度曲线：根据 r/T=0，可得可得XA=f（Topt）（2）设定第一段入口温度为）设定第一段入口温度为To1，To1值为催化剂活性温度，是一固定值。值为催化剂活性温度，是一固定值。（3）设第一段出口转化率为）设第一段出口转化率为XA1，（4）根据绝热温升线可计算第一段的出口温度）根据绝热温升线可计算第一段的出口温度T1。T=To+（XA-XAO）（6-47）其中，其中，=NTOyAO（-HR）/NT2CP第九十页，本课件共有108页(5)(5)由由即，第一段出口速率等于第二段入口速率，可求出第二段入即，第一段出口速率等于第二段入口速率，可求出第二段入口温度口温度TO2。（6）根据：式（）根据：式（6-46）和绝热温升线（和绝热温升线（6-47）第九十一页，本课件共有108页（7）重重复复（5）、（6）步步算算出出最最后后一一段段的的XAn，检检查查算算出出的的最最后后一一段段的的的的出出口口转转化化率率与与生生产产任任务务是是否否一一致致，若若不不一一致致，重重新新设设定定XA1重重复复（3）-（6）步步，直直至至直直至至算算出出最最后后一一段段的的的的出出口口转转化化率与生产任务一致。率与生产任务一致。（8）检查计算结果与平衡线及最佳温度线的关系。）检查计算结果与平衡线及最佳温度线的关系。T=To+（XA2-XA1）（6-47）其中，=NTOyAO（-HR）/NTCP求解出第二段出口转化率与出口温度XA2，T2。3催化剂有最高耐热温度催化剂有最高耐热温度Tmax限制限制设计思路：设计思路：（1）反应气体混合物的初始组成，催化剂动力学参数，作过程的平衡曲线与最佳温度曲线。第九十二页，本课件共有108页（2）假定第一段入口温度为）假定第一段入口温度为To1，根据：式（根据：式（6-46）和绝热温升线（和绝热温升线（6-47）a平衡曲线：根据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。亦可由 第九十三页，本课件共有108页与与 T=To+T=To+（X XA A-X-XAOAO）（6-476-47）计算出计算出X XA1A1。（3）根据式（）根据式（6-43）即，第一段出口速率等于第二段入口速率，可求出第二段入口温度TO2。（4）重复（2）步求出第二段出口温度与出口转化率 T2，XA2，及重复（3）步可求出第三段的入口温度To3。计算过程中，若Ti 出口温度超过Tmax应以Tmax为出口温度Ti。并用式（6-47）或式（6-47）结合式（6-48）计算出口转化率XAi。（5）检查算出的最后一段的的出口转化率与生产任务是否一致，若不一致，重新设定TO1重复（2）-（5）步，直至算出最后一段的的出口转化率与生产任务一致。（6）检查计算结果与平衡线及最佳温度线的关系。第九十四页，本课件共有108页4 4例题例题 某日产干吨两段间接换热式径向氨合成塔某日产干吨两段间接换热式径向氨合成塔合成氨的速率：合成氨的速率：R=dyNH3/do=kT（1+yNH3）2Kf2fH21.5fN2/fNH3fNH3/fH21.5解：（解：（1）不受催化剂温度限制的影响时，）不受催化剂温度限制的影响时，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 第九十五页，本课件共有108页b.b.由由 即，第一段出口速率等于第二段入口速率，可求出第二段入即，第一段出口速率等于第二段入口速率，可求出第二段入口温度口温度TO2=371.7。c.再由再由T=To+（XA2-XA1）（6-47）其中，=NTOyAO（-HR）/NTCP求解出第二段出口转化率与出口温度XA2，T2。T2=461.4,yNH3,2=0.2146,VR2=28.2155 m3