1、1第七章气固相催化反应流化床反应器2 流态化现象：使微粒固体通过与气体或液体接触而转变成类似流体的操作。固体颗粒层与流体接触的不同类型：流体流速增加固定床初始流态化散式流态化聚式流态化腾涌稀相流态化液体气体34流化床的基本概念 当通过床层的流体流量较小时，颗粒受到的升力（浮力与曳力之和）小于颗粒自身重力时，颗粒在床层内静止不动，流体由颗粒之间的空隙通过。此时床层称为固定床。随着流体流量增加，颗粒受到的曳力也随着增大。若颗粒受到的升力恰好等于自身重量时，颗粒受力处于平衡状态，故颗粒将在床层内作上下、左右、前后的激烈运动，这种现象被称为固体的流态化，整个床层称为流化床。5流化床类似液体的性状 轻的
2、物体浮起；表面保持水平；固体颗粒从孔中喷出；床面拉平；床层重量除以截面积等于压强67流化床的优点 1 颗粒流动类似液体，易于处理，控制；2 固体颗粒迅速混合，整个床层等温；3 颗粒可以在两个流化床之间流动、循环，使大量热、质有可能在床层之间传递；4 宜于大规模操作；5 气体和固体之间的热质传递较其它方式高；6 流化床与床内构件的给热系数大。8流化床的缺点 1气体的流动状态难以描述，偏离平推流，气泡使颗粒发生沟流，接触效率下降；2颗粒在床层迅速混合，造成停留时间分布不均匀；3脆性颗粒易粉碎被气流带走；4颗粒对设备磨损严重；5对高温非催化操作，颗粒易于聚集和烧结。9流化床的工业应用 第一次工业应用
3、：1922年 Fritz Winkler获德国专利，1926年第一台高13米，截面积12平方米的煤气发生炉开始运转。目前最重要的工业应用：SOD(Standard Oil Development Company)IV型催化裂化。10散式流化和聚式流化（1）散式流态化 随着流体流量的加大，床层内空隙率增大，颗粒之间间距加大，而颗粒在床层中分布均匀，流体基本上以平推流形式通过床层，人们称这种流化形式为散式流态化。11（2）聚式流态化 在此类流态化形式中，床层明显地分成两部分。其一是乳化相：固体颗粒被分散于流体中，单位体积内颗粒量类似于散式流化床的初始流化状态。其二是气泡相：流体以气泡形式通过床层。
4、12两种流态化的判别 一般认为液固流态化为散式流态化而气固之间的流化状态多为聚式流态化。为散式流态化 为聚式流态化100ReemfPmfmfDLFr100ReemfPmfmfDLFr13弗鲁特数：gduFrP2mfmf 雷诺准数：mfPmfReud umf：初始流化速度 dP：颗粒平均粒径，P：流体密度，颗粒密度 Lmf：初始流化时的浓相段床高 De：流体的扩散系数：流体粘度 14浓相段和稀相段 当流体通过固体床层的空塔速度值高于初始流化速度但低于逸出速度，颗粒在气流作用下悬浮于床层中，所形成的流固混合物称为浓相段。在浓相段上升的气泡在界面上破裂，气泡内颗粒以及受气泡挟带的乳化相中颗粒将被抛向
5、浓相段上方空间。这段空间称为稀相段或称分离段。15颗粒含量床高稀相段浓相段浓相段和稀相段16流态化的不正常现象 沟流：由于流体分布板设计或安装上存在问题，使流体通过分布板进入浓相段形成的不是气泡而是气流，称沟流。沟流造成气体与乳化相之间接触减少，传质与反应效果明显变差。节流（腾涌）1718流化床的工艺计算 1 初始流化速度：颗粒开始流化时的气流速度（气体向上运动时产生的曳力）（床层体积）（固体颗粒分率）（颗粒密度），即：)(1()(1)(gsmfmfgsmfmfttLpLAWAp19 将上式与固定床压降方程(Ergun方程)相结合，可得临界流化速度计算式。Ergun方程：与考虑固定床压降时的方
6、程对照：可以看出所作简化。ps20g3BB2ps03B2B175.11150duduLp s2mg3BBm175.1Re150dddulP20 前一项为粘滞力损失，后一项为动能损失。合并两式并整理：低雷诺数时，粘滞力损失占主导，忽略后一项：2gsg3p2mfgp3mfsmfgp3mf2smf75.11150dudud2gsg3pmfgp3mf2smf1150dud21 解得：高雷诺数时，动能损失占主导，忽略前一项：解得：20Re1150mf3mfgs2p2smfdu 2gsg3p2mfgp3mfs75.1dud1000Re75.13mfggsps2mfdu22 对中等雷诺数，两项都要考虑。计算
7、出临界流化速度后要进行验算，看雷诺数是否在适用范围之内。2 带出速度（终端速度）：当流体对颗粒的曳力与颗粒的重量相等，颗粒会被流体带走：CD-曳力系数2t2pgDgs3p4216udgCd23 对于单颗粒，有半经验公式：21gspD3122gs21Dgs2pD52.520000Re50043.02254500Re2Re10182ReRe24duCdguCgduCp对应对应对应24 以上计算是针对一个颗粒的，在流化床内由于颗粒间有相互影响，故逸出速度由此速度值再加以校正而得。uT=Fu Re10时，Re-F见下图2526 3 反应器内径的计算 VG：气流的体积流量m3s-1 dT：流化床内径m
8、u：气流的空塔流速m.s-1 可见，流化床的内径取决于气流的空塔气速，而流化床的空塔气速应介于初始流化速度（也称临界流化速度）与逸出速度之间。即维持流化状态的最低气速与最高气速之间。uVd4GT27 例8-1 计算萘氧化制苯酐的微球硅胶钒催化剂的起始流化速度和逸出速度 已知催化剂粒度分布如下：催化剂颗粒密度P=1120kg.m-3 气体密度=1.10kg.m-3 气体粘度=0.0302mPas目数120100-12080-10060-8040-60120 100-120 80-100 60-80 40-60 40 dPi 0.121 0.133 0.163 0.208 0.298 0.360
9、Piidx 0.99 0.752 0.797 1.680 0.839 0.139 mm192.0139.0839.0680.1797.0752.099.011PPiidxd7.330408.07.33Re212P3P2mfgd0568.07.331002.381.91.111201.11092.10408.067.1133Re212534mf30 3计算逸出速度（ut）:设Rem21524P2Psm2956.01002.31881.91.111201021.118gdu31 复核Re值 假设Rem2合理。由Re=1.3，Re10可得F=123.11002.31.12956.01021.1Re5
10、4Pud32 浓相段高度的计算 催化剂在床层中堆积高度称静床层高度(L0)。在通入气体到起始流化时，床高LmfL0。若继续加大气量，床层内产生一定量的气泡，浓相段床高(Lf)远大于静床层高度。关于浓相段床高的计算通常用计算床层空隙率(f)来获得。令床层膨胀比Rfmfmff11LLR33 0.2ReP1 1 ReP 200 200ReP500 500ReP n=2.39 则 Lf=RLmfnuu1tf03.0TPRe5.1735.4ddn1.0TPRe1845.4ddn1.0Re45.4n34 稀相段床高的估算 稀相段也称分离段，主要是用来保证床内因气泡破裂而挟带固体颗粒重新回到浓相段所需空间。
11、稀相段床高可由化工原理中非均相分离过程计算而得，也可由下述经验方程估算。1.155.1P032ArRe102.1LL35 例8-2 例8-1中的催化反应过程，若操作气速取12cm.s-1，催化剂装填高度L0=20cm，气体流量为122m3h-1，试估算流化床内径以及浓相段、稀相段床高。解 1计算流化床内径 2计算流化床浓相段床高m6.012.0360012244GTuVdfmfmff11LLRnuu1tf36 当0.2ReP1时03.0TPRe5.1735.4ddn8392.01002.31.112.01092.1Re54P379.48392.06.01092.15.1735.403.04n8
12、139.02965.012.0379.41f687.28139.015.01Rcm74.5320687.2mff RLL37 3计算稀相段床高 4床层总高 L=Lf+L2=53.74+131.36=185.1 cm1.155.1P032ArRe102.1LL7.931002.31.1112081.91.11092.1Ar25342P3Pgdcm36.1317.938392.020102.11.155.132L38流化床的热传递 流化床的热量传递过程大体可分为：固体颗粒之间的热量传递；气体与固体之间的热量传递；床层与床壁（包括换热器）之间的热量传递。由于流化床中颗粒处于高度运动状态，而固体的导热
13、系数较大，因此传热速率很快。床层中温度基本上可以认为是一致的。39 流化床层与器壁的给热系数直到目前为止仍只能通过将实验数据归纳成准数方程而获得。21PPPPPrRe1075.0Nucc40 流化床层与竖放的换热器器壁之间给热系数计算式为 注意：是有单位的，其单位为s.cm-2 床层与横放的换热器器壁之间传热时，给热系数计算式为 43.066.0P8.0PPP43.023.02PrRe110844.1NuccC44.0P3.044.01PrRe66.0Nu41 流化床传热小结 水平管的给热系数比垂直管低515，因此倾向于使用垂直管。颗粒的导热系数和床高对给热系数影响不大；给热系数随颗粒比热的增
14、大而增大，随粒径的增大而降低；流体的导热系数对给热系数hw起最主要的影响，hw与n成正比，n=1/2-2/3。床层直径的影响难于判定；床内管径小时给热系数大；42床层中气泡行为 当气体通过床层时一部分气体与颗粒之间组成乳化相，其余气体以气泡形式通过乳化相。由于气体上升速度与乳化相速度不同，存在明显的速度差异，气泡在上升过程中必然会挟带气泡周围一定量的乳化相物质。气泡在上升时其尾部形成负压，将吸入部分乳化相物质随其上升，这部分称尾涡。43 气泡上升时气泡外侧一定厚度的乳化相将随气泡一起上升，这部分被称为气泡云。尾涡与气泡云统称为气泡晕。44流化床的鼓泡床模型 鼓泡床模型对流化床运动形态作如下简化
15、:（1）认为床层主体部分气泡大小均一且均匀分布于床层之中。（2）床层中乳化相处于起始流化状态，超过起始流化态的气体将以气泡形式通过床层。45（3）床层可分为气泡、气泡晕及乳化相三部分。在气泡、气泡晕和乳化相之间的传质过程是一个串联过程。（4）在 时，进入稀相段的气体只有气泡破裂而逸出的气体，故稀相段气体组成与离开浓相段的气泡中气体组成相同。6mf0uu46反应过程的估算 在流化床的浓相段中，对气体中反应物A而言，存在如下关系：反应消耗反应消耗反应消耗的总消失量乳化相气泡晕气泡相传递传递Aecb47AeAcceAcAbbcccecebbcbcccKccKSkKSkK传递速率相间交换系数AereA
16、crcAbrb31AerAcrAbr31msmolmsmolcfkcfkcfkcfkcfkcfkecb气泡颗粒反应速率气泡体积颗粒体积48 根据此表可得A组分的物料衡算。以单位气体体积为基准（总消失量）=（在气泡中反应的量）+（转移到气泡晕中的量）（转移到气泡晕中的量）=（在气泡晕中反应掉的量）+（转移到乳化相中的量）（转移到乳化相中的量）=（在乳化相中反应掉的量）49 式中 是流化床内总反应速率常数。对该方程进行积分 边值条件为：AbbrAbbAbAbddddddddAcfKlculctltc组分总消失量brKluKcfcddbbrAbAbAfAAfAbfA0AAb00 xxccLlxccl
17、50 已知cA0，cAf（或xAf），利用该式可求得浓相段床高Lf，进而求出催化剂用量。已知cA0，Lf，可求得气体的出口浓度cAf（或转化率xAf）。以反应动力学方程为一级的反应为例：AbbrAbbAbddddAcKlcutc组分总消失量51AereAeAcceAcrcAcAbbcAbrbckccKckccKck在乳化相中反应量转移到乳化相的量在气泡晕中反应量转移到气泡晕的量在气泡中反应量AeAcceAbrbAbbrAbbddccKckcKlcuAeAcceAcrcAcAbbcccKckccKAereAeAcceckccK52 联解此方程，消除cAc，cAe整理后可得ecercbcrbrbr
18、1111KkKkkK53 由 边值条件 代入 若浓相段床高为Lf，则出口气体浓度及转化率为AbbrAbbddcKlcuAfAbf0AAb0ccLlcclbfr0bbrAbAbddfAfA0uLKluKccLccfbbrAfexp1LuKx54 若要求出口转化率为xAf，则需浓相段床高Lf为 b，c，e及Kbc，Kce值由经验公式计算。b的值在0.001-0.01之间。由于该值较小，对计算影响不大。Afbrbf1lnxKuLbwmfmf21bmfmfmfc711.031VVugdu55cbbbmfe11bmfbuuu gud2tb711.0415b2ebmfbc85.55.4dgDduK213b
19、bemfce78.6duDK56 例8-3 计算萘氧化制苯酐的流化床反应器气体出口转化率.已知:1.催化剂：微球硅胶钒催化剂（同例8-1）平均粒径 密度 2气体性质 气体密度 气体粘度 扩散系数m1092.14Pd3Pkg.m1120smPa302.012scm204.0D3kg.m1.157 3流化床特性 静床层高 床层直径 空隙率 操作气速 4反应动力学方程 cm20mf0 LLm6.0Td5.00mf1m.s12.0u14.3AArcR58 解 1计算起始流化速度与逸出速度 2计算操作条件下的空隙率及膨胀比 空隙率 床层膨胀比 浓相段高 稀相段高 13mfm.s1011.8u1tm.s2
20、956.0u8.0f05.3Rcm74.53fLcm36.1312L59 3计算气泡上升速度m01762.081.9711.02956.0711.022tbgud1tbrm.s2956.0 uu1brmfbm.s75.4056.29811.012uuuu60 4计算b，c，e值 取b=0.01 查 图，当时，bwmfmf21bmfmfmfc711.031VVugdubwVVm1092.14Pd47.0bwVV61 代入 式中 代入3221.047.05.000811.001762.081.9711.05.000811.035.0121ccbbbmfe112746.075.40811.012bm
21、fbuuu8435.03221.001.02746.02746.015.01e62 5计算 、值bbcKbceK8488.9762.1981204.085.5762.1811.05.485.55.44152415b2ebmfbcdgDduK435.8762.175.40204.05.078.678.6213213bbemfceduDK63 6计算Kr值4694.28435.01435.84.313221.01849.94.3101.04.31111ecercbcrbrbrKkKkkK64 7计算出口气体中萘的转化率9615.074.5375.404694.2exp1exp1fbbrAfLuKx
