化学反应工程固定床反应器教材

1、第六章 固定床反应器凡是流体通过固定的固凡是流体通过固定的固体物料所形成的床层而体物料所形成的床层而进行反应的装置都称作进行反应的装置都称作固定床反应器固定床反应器。如：气如：气-固相催化反应器、固相催化反应器、气气-固相非催化反应器。固相非催化反应器。乙苯乙苯水蒸气水蒸气催化剂催化剂产品产品测测温温口口6-1乙苯脱氢的绝热床反应器乙苯脱氢的绝热床反应器6.1.1 固定床反应器的优缺点 固定床层内的气相流动接近平推流，有利于实现较高的转化率与选择性； 可用较少量的催化剂和较小的反应器容积获得较大的生产能力； 结构简单、催化剂机械磨损小，适合于贵金属催化剂； 反应器的操作方便、操作弹性较大。u相

2、对于流化床反应器，固定床反应器缺点 催化剂颗粒较大，有效系数较低； 催化剂床层的传热系数较小，容易产生局部过热； 催化剂颗粒的更换费事，不适于容易失活的催化剂。 6.1.2 固定床反应器类型u固定床反应器形式多种多样，按床层与外界的传热方式分类，可有以下几类： 绝热式固定床反应器， 多段绝热式固定床反应器， 列管式固定床反应器， 自热式反应器。1.绝热式固定床反应器 u反应器外壳包裹绝热保温层，使催化剂床层与外界没有热量交换。中空圆筒的底部 放置搁板，上面堆放固体催化剂。气体从上而下通过催化剂床层。u结构简单，床层横截面温度均匀。单位体积内催化剂量大，即生产能力大。但只适用于热效应不大的反应。

3、67原料气产物催化剂催化剂绝热式固定床反应器 p绝热式固定床反应器可分为轴向反应器和径向反应器。(1)轴向绝热式固定床反应器u如图(a)所示。这种反应器结构最简单，实际上是一个容器，催化剂均匀堆置于床内，预热到一定温度的反应物料自上而下流过床层进行反应，床层同外界无热交换。8(2)径向绝热式固定床反应器u如图 (b) 所示。径向反应器的结构较轴向反应器复杂，催化剂装载于两个同心圆构成的环隙中，流体沿径向流过床层，可采用离心流动或向心流动。u径向反应器的优点是流体流过的距离较短，流道截面积较大，床层阻力降较小。轴向反应器与径向反应器轴向反应器与径向反应器(a)(b)2.多段绝热式固定床反应器u

4、热效应大，常把催化剂床层分成几段(层)，段间采用间接冷却或原料气（或惰性组分）冷激，以控制反应温度在一定的范围内 。u 图 (c) 是用于 SO2 转化的多段绝热反应器，段间引入冷空气进行冷激。对于这类可逆放热反应过程，通过段间换热形成先高后低的温度变化，提高转化率和反应速率。多段绝热床多段绝热床根据段间反应气体的冷却或加热方式，多根据段间反应气体的冷却或加热方式，多段绝热床又分为中间间接换热式和冷激式。段绝热床又分为中间间接换热式和冷激式。特点特点：催化剂床层的温度波动小。催化剂床层的温度波动小。缺点缺点：结构较复杂，催化剂装卸较困难。结构较复杂，催化剂装卸较困难。间接换热间接换热 原料冷激

5、原料冷激 非原料冷激非原料冷激多段固定床绝热反应器多段固定床绝热反应器原料原料原料原料冷冷激激剂剂原料原料产品产品产品产品产品产品3.列管式固定床反应器u 热效应较大，不宜采用绝热式反应器，可采用换热式固定床反应器。此设备如同列管式换热器，又称为列管式固定床反应器。u 如图(d)所示，反应器由多根反应管并联构成，管径一般为25 30，管数可达万根以上。管内装催化剂，传热介质流经管间进行加热或冷却。列管式固定床反应器 u列管式反应器优点： 传热较好，管内温度较易控制； 返混小、选择性较高； 只要增加管数，便可有把握地进行放大； 对于极强的放热反应，还可用同样粒度的惰性物料来稀释催化剂 适用: 原

6、料成本高，副产物价值低以及分离不是十分容易的情况。4. 自热式反应器u采用反应放出的热量来预热新鲜的进料，达到热量自给和平衡，其设备紧凑，可用于高压反应体系。u但其结构较复杂，操作弹性较小，启动反应时常用电加热。226.1.3 传热介质传热介质的选用根据反应的温度范围决定，其温度与催化床的温差宜小，但又必须移走大量的热，常用的传热介质有：1.1.沸腾水：沸腾水：温度范围100300。使用时需注意水质处理，脱除水中溶解的氧。2.2.联苯醚、烷基萘为主的石油馏分：联苯醚、烷基萘为主的石油馏分：粘度低，无腐蚀，无相变，温度范围200 350233.3.熔盐：熔盐：温度范围300400，由无机熔盐KN

7、O3、NaNO3、NaNO2按一定比例组成，在一定温度时呈熔融液体，挥发性很小。但高温下渗透性强，有较强的氧化性。4.4.烟道气：烟道气：适用于600700的高温反应。24256.2 6.2 固定床的传递特性固定床的传递特性 气体在催化剂颗粒之间的孔隙中流动，较在管内流动更容易达到湍流。 气体自上而下流过床层。6.2.1颗粒尺寸 颗粒尺寸是颗粒体系的重要参数，常用粒径来表示。球形粒子的粒径是其直径，其他形状的粒子粒径则需定义。 颗粒的定型尺寸颗粒的定型尺寸最能代表颗粒性质的尺寸为颗粒的当量直径。对于非球形颗粒，可将其折合成具有相同的体积（或外表面积、比表面积）的球形颗粒，以当量直径表示。如体积

8、、外表面积、比表面积当量直径。2627 体积当量直径：体积当量直径：(非球形颗粒折合成同体积的球形颗粒应当具有的直径） 外表面积当量直径：外表面积当量直径： (非球形颗粒折合成相同外表面积的球形颗粒应当具有的直径）31S3S66VddVV球形体积：S2SSddSa球形外表面积：28 比表面积当量直径：比表面积当量直径： (非球形颗粒折合成相同比表面积的球形颗粒应当具有的直径） 混合粒子的平均直径：混合粒子的平均直径：（各不同粒径的粒子直径的加权平均，xi 是直径为d i 的粒子的质量分率）。SSVS32SSV6666SVSddddVSSiidxd1m296.2.2 床层空隙率 B 单位体积床层

9、内的空隙体积(没有被催化剂占据的体积，不含催化剂颗粒内的体积)。 若不考虑壁效应，装填有均匀颗粒的床层，其空隙率与颗粒大小无关。颗粒密度床层堆积密度，床层体积颗粒体积床层体积空隙体积PBPBBPB111VV 床层空隙率是一个重要的参数，影响因素是颗粒形状及大小、粒度分布、颗粒与床层直径比和颗粒的装填方式。 壁效应壁效应 床层空隙率沿床层 径向分布不同，离 壁面约一个粒子直 径处的床层空隙率 最大。 3031 床层内空隙率径向分布不均匀，引起各处的流速不同，因而床层内各处的传热和停留时间也不一样。 为减少壁效应的影响，设计时要求床层直径至少要大于颗粒直径的 8 8倍倍 以上。6-9填充床的空隙率

10、填充床的空隙率床层空隙率床层空隙率B光滑均一光滑均一光滑非均一光滑非均一粘土粘土球形球形圆柱形圆柱形不规则不规则光滑均一光滑均一刚玉均一尺寸刚玉均一尺寸1 / 4 英寸陶质拉西环英寸陶质拉西环熔融刚玉熔融磁铁铝砂6.2.3 床层压降 气体流动通过催化剂床层的空隙所形成的通道，与孔道周壁摩擦而将产生压降。 压降计算通常利用厄根（Ergun）方程：33Bggmsmm1Re:Reud修正的雷诺数，式中： s2mg3BBm175. 1Re150dddulP34 厄根（Ergun）方程中其它参数： 可用来计算床层压力分布。 如压降不大，床层各处物性变化不大，可视为常数，压降将呈线性分布。床层空隙率气体密

11、度颗粒当量直径床层高度空塔气速平均流速Bgsmdlu 当Re m1000时，厄根(Ergun)方程则变为： 催化剂床层压降还有许多计算式，具体参考有关的资料。35LduP3BBs2mgm1Re150LduP3BBs2mg175. 1 影响床层压力降的最大因素： 床层的空隙率 流体的流速 两者稍有增大，会使压力降产生较大变化。 降低床层压降的方法： 增大床层空隙率，如采用较大粒径的颗粒； 降低流体的流速，但要考虑这会使相间的传质和传热变差，需综合考虑。3637例6.1内径为50mm的管内装有4m高的催化剂层，催化剂为球体，催化剂的粒径分布如表所示。空隙率 B = 0.44。在反应条件下气体的密度

12、 g = 2.46 kgm-3，粘度 g = 2.310-5 kgm-1s-1，气体的质量流速 G = 6.2kg m-2s-1。求床层的压降。粒径粒径 dS /mm3.404.606.90质量分率质量分率 0.600.250.1538解： 求颗粒的平均直径 计算修正雷诺数m1096. 3mm96. 390. 615. 060. 425. 040. 360. 0131Siidxd190644. 01103 . 22 . 61096. 31Re53BgSmGd39 计算床层压降Pa10898. 1444. 044. 0146. 21096. 32 . 675. 11903150175. 1Re1

13、50175. 1Re15053323BBgS2m3BBSg2mmLdGLdup改性的 ZSM-5 催化剂其中催化剂的比表面积为300600m2/g，孔容体积为 0.30.8cm3/g甲醇制烯烃项目：甲醇制烯烃项目：反应器结构的计算 1 催化剂的填充量 反应管长度的计算 管束尺寸和反应管的排列 反应床层压降opmmbmbbpombudRRbafdufHp 112例题 乙烯在银催化剂上氧化制环氧乙烷，年产环氧乙烷1106 kg,采用二段空气氧化法。主要反应为： C2H4+1/2O2 C2H4O H1= -103.4kJ/mol(25) C2H4+3O2 CO2+2H2O H1= 1323kJ/mo

14、l(25)根据下列给出中试的数据，估算第一反应器尺寸。 （1）进入第一反应器的原料气组成为：组成C2H4O2CO2N2C2H4Cl2摩尔分数/%3.56.07.782.8微量（2）进入第一反应器进料温度为210 ，反应温度为250 ，反应压力为980.0665KPa,转化率为20%，选择性为66%，空速为5000/h.（3）进入第一反应器采用列管式固定床反应器，列管为直径272.5mm，管长5m,催化剂填充高度5.7m.（5）催化剂为球形，直径dp为5mm,床层孔隙率eB为0.48。（4）管间采用导热油强制外循环换热。导热油进口温度为230 ，出口温度为235 ，导热油对管外壁传热系数a2可取

15、2721kJ/(m2h )（6）年工作7200，反应后分离，精制过程回收率为90%，第一反应器所产环氧乙烷占总产量的90%。4.2.4 固定床中的传热 固定床中的传热组成 颗粒内传热 颗粒与流体间的传热，传热系数 hP 床层与器壁的传热，传热系数 hw、h049u拟均相模型l将包括颗粒与气相流体的床层看作为均一的固体物质，其传热特性用一个有效热导率 e 来表征。 有效热导率 e 根据轴向和径向的传热，可分为轴向和径向的有效热导率 （ ez 和 er）。50511.流体与颗粒外表面间的传热 与传质相仿，由于层流边界层的存在，造成了气流主体与催化剂颗粒外表面存在温度差。因而，必然存在热量传递。 单

16、位时间内从颗粒外表面传递到气相主体的热量为：GSSPddTTShtQ -单位时间传递的热量，J/h； h P -流体对颗粒的传热系数，J/m2hK； T g -气相主体的温度，K； TS -为催化剂颗粒外表面处温度，K； SS -催化剂外表面积，m2； -颗粒表面利用系数，球体 = 1， 圆柱体 = 0.91， 其他形状 = 0.9052tQdd53u传热系数的计算(通过传热因子法计算)： 式中：32PPHPrJGchsPasmkg,g12H气相粘度气体质量流率无量纲数。传热因子GJ1111PKmWKkgJ气体导热系数气体恒压热容gcuPr -普兰特准数 对于气体 Pr = 0.6 1.0 对

17、于液体 Pr = 2 400 适用范围:54ggPPrc35. 0H/3023. 0/876. 2JGdGdGdPPPB1000010GdPl JH 是雷诺数的函数556000Re300Re28. 1J300Re06. 0Re26. 2J1000Re50Re613. 0J50Re01. 0Re904. 0Jm41. 0mHm51. 0mHm51. 0mHm51. 0mH)1 (6ReBSemGdSG或:2.固定床的有效热导率 e 床层的有效热导率是颗粒与流体间的对流传热、颗粒及流体本身的导热以及床层内的辐射传热等几类作用的综合表现，是流体及颗粒特性及流动状态的函数。 通常固定床的热量主要是在中

18、心与管壁间做径向的传递，所以忽略轴向导热 ez。56u确定有效热导率方法 先测定床层中的温度分布，然后根据传热方程式来反求出 e。 -代表流体静止时的床层的热导率。 上式中第二项代表流体流动混合对径向传热的增值。57PrRe)(0Pgege0e 的物理意义代表横向传质与流动方向传质速度之比。 的物理意义代表颗粒间距与粒径比的影响。 g - 流体的有效热导率。58n 的计算：第一项代表床层空隙部分对传热的贡献；第二项代表颗粒部分的贡献。59oesprsBprvBgedhdh3211110 hrv 和 hrs 分别为空隙和颗粒的辐射给热系数。单位：6031001)1 (211162. 1 mBBB

19、rvTh)/(2KmW31002162. 1mrsTh s 和 分别为颗粒和流体的热导率； 为粒子表面的热辐射率； 代表颗粒接触点处流体薄膜的导热影响。61216. 026. 0)(212B3.床层与器壁间的传热系数hw及h0 一维模型，床层径向温度一致，传热速率为： A-传热面积 tm、tw-床层平均温度、器壁温度。62)(0wmttAhqybadddhetpp)(210 二维模型，需要考虑床层径向温度分布，传热速率为：63)(wRwttAhqpwpwpwpwdhdhdhRePr1110644.2.5 固定床中的传质1.流体与催化剂颗粒外表面之间的传质 前面讨论宏观动力学是采用催化剂表面的温

20、度、浓度，然而催化剂表面的温度浓度却难于测量。 本节讨论催化剂表面的温度浓度与气流主体的温度浓度之间的关系，通过可测量的量建立动力学关系。65 流体与催化剂表面之间的传递属外扩散过程。 催化剂表面处滞流边界层中的阻力，使气流主体与催化剂颗粒表面间存在浓度的差异。 组分A从气相主体达到催化剂颗粒外表面上，必须通过扩散运动，流体通过边界层向催化剂颗粒表面扩散的传质速率方程为：ASAGSgAddccSktn -单位时间传递A 物质的摩尔数，mol/s kg -气相传质系数，cm/s CAG -气相主体中A物质的浓度，mol/cm3 CAS -催化剂外表面处A物质浓度，mol/cm3 SS -催化剂外

21、表面积，cm2 -颗粒表面利用系数，球体 = 1， 圆柱体 = 0.91，其他形状 = 0.9066tnddA67u气相传质系数 kg 整个传质方程的核心，反映了传质过程阻力的大小，总包了各种条件对传质的影响。 传质系数对传质速率影响极其重大，传质系数越小，阻力就越大。 传质系数与颗粒的几何尺寸及形状、流体力学条件、流体的物理性质有关。 传质系数 kg 由实验关联式计算。常用 JD 因子计算法关联： 式中：6832ggDScGkJ施密特准数；气体质量流率；气相密度传质因子；,Sc:Scsmkg:mkg:gg123gDDGJ69.cm,1ReRe19. 1,6000Re300Re10. 2,30

22、0Re3 . 0BsBgsm41. 0mDm51. 0mDmD催化剂床层空隙率颗粒比表面当量直径是雷诺数的函数：dGdJJJ。气相分子扩散系数；气相粘度12gsmsPaD70 比表面当量直径：单颗催化剂的外表面积折合成直径为 dS 的球形颗粒应有的外表面积。n注意 Re 的不同定义。2.传质过程对反应的影响 流体与催化剂颗粒外表面之间存在的层流边界层，因此，造成流体主体的气相浓度与颗粒外表面处的气相浓度不同。 宏观动力学方程是以表面气相浓度为计算基准的，外扩散过程直接影响反应的结果。 连续稳定过程，组分A在单位时间内扩散到颗粒外表面处的量等于催化剂中反应掉的量71即： 上式将 CAG 与 CA

23、S 关联起来， CAG 是气相主体浓度，可以直接测定，因而解决了宏观动力学中 CAS 的计算问题。 若本征动力学方程为：72 ASckfVVRccSktnSSAASAGSgAddAckfrA则：其中： Da为坦克莱(Damkohler)数，其物理意义是反应速率与外扩散速率的比值，反映了体系中外扩散的影响程度。数值越大，反应速率越快，外扩散的影响就越大。 Da0，外扩散可忽略。73SgSDaSkkV)(Da)(ASASSgSASAGcfcfSkkVcc74 一级反应： 二级反应： 其他类型(N级)动力学方程：AGASASASAGASASDa11Da)(cccccccfDa21Da41DaAGAS

24、2ASASAG2ASAScccccccf)(AGAScc 有外扩散影响的宏观动力学方程为：75)()()(AGAGSSgAGASAccVSkckfckfR76特殊情况： 外扩散控制-反应速率常数 k 比传质系数 kg 大得多，则颗粒外表面处 A 的浓度为零。 内扩散或动力学控制-反应速率常数 k 比传质系数 kg 小得多，则颗粒外表面处 A 的浓度与气相主体浓度相等，外扩散可不予考虑。AGSgA6cdkRAGAckfR77例 实验室中进行苯加氢反应，在1013.3 kPa下操作，气体质量速度G=3000 kgm-2h-1，催化剂为 89mm圆柱体，颗粒密度P=0.9gcm-3，床层堆积密度B=

25、0.6gcm-3，在反应器某处气体温度为220，气体组成为10苯，80氢，5环己烷和5甲烷(体积分率)，测得该处宏观反应速率 (-RA)=0.015molh-1g-1(cat) 。试估算该处催化剂的外表面浓度。注：气体粘度 = 1.410-4 gcm-1s-1， 扩散系数 D = 0.267cm2s-1。7879解： 计算催化剂的粒径 dScm8308. 0267. 34524. 066cm4524. 09 . 08 . 044cm267. 38 . 09 . 08 . 0242SSS322S222SSVdLdVLddS80 计算床层中气体的雷诺数7 .814333. 011000036009

26、1. 0104 . 1100030008308. 01Re333. 09 . 06 . 0114BgSmPBBGd81 计算 JD 和 kg 值0762. 07 .81419. 1Re19. 141. 041. 0mDJ33mgg.cm1056. 322027383144 .143 .1013RTpM4 .141605. 08405. 028 . 0781 . 0miiMxM82132343332gggDgcm.s397. 6267. 01056. 3104 . 11056. 31030000762. 0DGJk83 计算 CAG 和 CAS32AAGmkmol10472. 222027331

27、4. 81 . 03 .1013RTpc3538SgSPAASAGmkmol1092. 8cmmol1092. 891. 0267. 3397. 64524. 09 . 03600015. 0SkVRcc8432525AGASmkmol10463. 21092. 810472. 21092. 8 cc853.流体与颗粒外表面间的浓度差和温度差 单位时间内传递的热量必然等于单位时间内反应放出的热量。 单位时间内反应放出的热量为： 传递的热量=反应热量)()(ddSAHVRtQ)()()(GSSPSTTShHVRA上式整理:前面的传质过程中，已知：由(A)和(B)式可得：86)()()(SGSSP

28、AHVTTShR(A) )()(ASAGSSgAccVSkR(B)()()(ASAGSSgSGSSgccVSkHVTTS整理上式，得： 上式关联了颗粒浓度与温度之间的变化关系。 可通过气相主体浓度和温度及颗粒表面浓度求颗粒表面温度。87(C)(ASAGggGSccHkTT88u进一步简化(C)式，前面有： JH 和 JD 相除，得：0.41mHm0.51mHm0.41mDm0.51mDm1.28ReJ6000Re3002.26ReJ300Re0.061.19ReJ6000Re3002.10ReJ300Re0.31.0761.191.28JJ6000Re3001.0762.102.26JJ300

29、Re0.06DHmDHm 由J(H) 和 J(D)的定义： 实验得到： 代入上式： 代入(C)式得：89076. 1ScPr/)/(JJPggg32ggPgDHckGkGc1Sc)(Pr/2/3)/(93. 0/Pggckg)(93. 0ASAGPgGSccHcTT90例试计算例5-6中催化剂的外表面处温度。已知反应热为(-H)=2.135105Jmol-1，气体的定压比热容cP=49Jmol-1K-1。解：C1 .2201 . 0220C1 . 0491056. 31092. 810135. 293. 093. 0oSo385PgASAGgSTcccHTT6.3 拟均相一维模型拟均相一维模型

30、概述，目前描述固定床反应器的数学模型可分为拟均相、概述，目前描述固定床反应器的数学模型可分为拟均相、非均相两大类。非均相两大类。 一、拟均相模型一、拟均相模型 忽略床层中催化剂颗粒与流体之间温度与浓度差别，将气忽略床层中催化剂颗粒与流体之间温度与浓度差别，将气相反应物与催化剂看成均匀连续的均相物系。相反应物与催化剂看成均匀连续的均相物系。（1）一维拟均相模型）一维拟均相模型只考虑沿气体流动方向的温度和浓度变化。根据流动形式只考虑沿气体流动方向的温度和浓度变化。根据流动形式还可分为平推流一维模型和轴向分散一维模型。还可分为平推流一维模型和轴向分散一维模型。（2）二维拟均相模型）二维拟均相模型 同

31、时考虑轴向和径向的温度和浓度分布。同时考虑轴向和径向的温度和浓度分布。 二、非均相模型二、非均相模型非均相非均相反应属于扩散和化学动力学共同控制时，则反应属于扩散和化学动力学共同控制时，则催化剂颗粒表面、内部、外部浓度不均一，传递阻力催化剂颗粒表面、内部、外部浓度不均一，传递阻力或传递与动力学阻力不可忽略，应计及催化剂的存在或传递与动力学阻力不可忽略，应计及催化剂的存在和计算宏观反应速率，称为和计算宏观反应速率，称为“非均相非均相”模型。模型。简言之即考虑颗粒与流体之间的温度差和浓度差。简言之即考虑颗粒与流体之间的温度差和浓度差。一般来说，模型考虑得越全面，对过程模拟越精确，一般来说，模型考虑

32、得越全面，对过程模拟越精确，但计算工作量也越大，甚至无法求解。因此，在工程但计算工作量也越大，甚至无法求解。因此，在工程计算允许的误差范围内应尽可能选用简单模型。计算允许的误差范围内应尽可能选用简单模型。6.3.1 等温反应器的计算等温反应器的计算床层温度均匀一致，反应速率常数为常数，反应速度仅与床层温度均匀一致，反应速率常数为常数，反应速度仅与浓度有关。按一维拟均相处理，设计方法与平推流相似。浓度有关。按一维拟均相处理，设计方法与平推流相似。对右图固定床反应器取一微元段进对右图固定床反应器取一微元段进行物料衡算（以行物料衡算（以cat的质量定义）的质量定义）0A0AxFdlTTAAdFF AF质量衡算质量衡算在管式反应器中垂直于流动方向取一个微元，以这个在管式反应器中垂直于流动方向取一个微元，以这个微元对微元对A组份组份做物料衡算：做物料衡算：输入输入 输出输出= 反应反应 积累积累FA FA+dFA (-r