9固定床吸附过程的计算

1、第四节 固定床吸附过程的计算固定床吸附器结构简单，但由于气体吸附过程是气固传质，对任一时间 或任一颗粒来说都是不稳定过程，因此固定床吸附器的吸附操作是非稳态的， 计算过程非常复杂， 一般要涉及到物料衡算方程、 吸附等温线方程和传热速率 方程及热量衡算。 而在气态污染物的吸附净化设计中， 由于所涉及到的物系是 低浓度的气态混合物， 且气量一般比较大， 吸附热相对较小， 因此可近似地按 等温过程处理，可不考虑传热速率方程和热量方程（升温脱附除外） 。这样在 设计过程中可采用简化了的方法进行近似计算， 计算时往往提出如下假设：（1） 气相中吸附质浓度低；（2）吸附操作在等温下进行；（3）传质区通过整

2、个床层 时长度保持不变；（ 4）床层长度比传质区长度大得多。 这些简化限制条件对目 前工业上应用的吸附器来说， 一般是符合的。 设计中较常采用的是希洛夫近似 计算法和透过曲线计算法。 计算过程一般是在吸附剂的选择、 吸附设备的选择 和吸附效率确定之后进行的。设计计算的任务是求出吸附器的床层直径和高 度，吸附剂的用量，吸附器的一次循环工作时间， 床层压降等。下面首先介绍 固定床吸附器的吸附过程。一、固定床吸附器的吸附过程 在固定床吸附器的吸附操作中，一般是混合气体从床层的一端进入，净化 了的气体从床层的另一端排出。 因此，首先吸附饱和的应是靠近进气口一端的 吸附剂床层。随着吸附的进行， 整个床层

3、会逐渐被吸附质饱和， 床层末端流出 污染物，此时吸附应该停止， 完成了一个吸附过程。 为了描述吸附过程，提出 了以下概念。（一）吸附负荷曲线与透过曲线1. 吸附负荷曲线 在实际操作中，对于一个固定床吸附器，气体以等速进入床层，气体中的 吸附质就会按某种规律被吸附剂所吸附。 吸附一定时间后， 吸附质在吸附剂上 就会有一定的浓度， 我们把这一定的浓度称为该时刻的吸附负荷。 如果把这一 瞬间床层内不同截面上的吸附负荷对床层的长度 （高度）作一条曲线， 即得吸 附负荷曲线。也就是说，吸附负荷曲线是吸附床层内吸附质浓度 x 随床层长度 z 变化的曲线。z图3 12 理想吸附负荷曲线图3-13吸附负荷曲线

4、在理想状态下，若床层完全没有阻力，吸附会在瞬间达到平衡，即吸附速 率无穷大，则在床层内所有断面上的吸附负荷均为一个相同的值，吸附负荷曲线将是一个直角形的折线，如图 3-12所示。但实际上是不可能的，在实际操 作中由于床层中存在着阻力，在某一瞬间床层内各个截面上的吸附负荷会有差 异，我们这时所绘制的曲线将是图 3-13所示的吸附负荷曲线。图中把曲线分 成了三个区域：饱和区（所有吸附剂已经达到了饱和）、传质区（有一部分吸 附剂还正在吸附）和未用区（所有吸附剂上均未有吸附质）。如果经过一段时 间的吸附，绘制另一时刻的吸附负荷曲线时，会发现曲线前进到了II线的位置，所以我们又形象地把吸附负荷曲线称为吸

5、附波或吸附前沿。当吸附波的下端到达床层未端时，说明已有吸附质漏出，这时床层被穿透，当床层被穿透的 这个时刻，称为破点。此时流出气体中吸附质的浓度称为破点浓度。在实际工作中，由于吸附剂中吸附质的浓度（即吸附负荷）不易测定，故 目前许多场合，曲线的纵坐标都以床层中混合气体的浓度 c来表示。因此，吸 附负荷曲线又可定义为在稳定吸附状态下，床层中气相中吸附质的浓度随床层 高度（长度）变化的曲线。由于床层的阻力不同，吸附负荷曲线会有不同的形状。床层阻力愈大，某 一时刻床层内各截面上浓度差别越大， 吸附负荷曲线也就变得越平缓，这当然 是我们不希望出现的情况。2. 透过曲线吸附负荷曲线表达了床层中浓度分布的

6、情况， 可直观地了解床层内操作的 状况。但要从床层中各部位采出吸附剂样品进行分析是相当困难的，这样易破坏床层的稳定。因此通常改用在一定的时间间隔内。 分析床层流出物中吸附质 浓度的变化，以流出物中吸附质浓度y为纵坐标。时间t为横坐标，则随时间 的推移可画出一条T -y曲线。如图3-14所示，开始时流出物中吸附质浓度为 yB，它是与吸附剂中的XB浓度相平衡的（XB为破点时床层出口端的吸附负荷）。 流出物中吸附质浓度开始上升，到T E时升到yE,即接近床层进口浓度，这时床 层已完全没有吸附能力，吸附波的末端也离开床层了。于是在T -y图上，从T B到T E呈现一个S型曲线，这条曲线称“透过曲线”。

7、它的形状与吸附负荷 曲线是完全相似的，只是方向相反。由于它与吸附负荷曲线成镜面对称相似, 所以也称吸附负荷曲线为“吸附波”或“传质前沿”。图3-14吸附床吸附过稈中流出物分析曲线由于透过曲线易于测定和标绘出来，因此也用它来反映床层内吸附负荷曲 线的形状，而且也能准确地求出破点。如果透过曲线比较陡，说明吸附过程比 较快，反之则速度较慢。如果透过曲线是一条竖直的直线，则说明吸附过程是 飞快的，是理想的吸附波。（二）保护作用时间保护作用时间是固定床吸附器的有效工作时间。它定义为从吸附操作开始 到床层被穿透所经历的时间称为保护作用时间，如图3-14所示的由T 0到T B所经历的时间，到达T B时，床层

8、内吸附剂还没有完全饱和。图中的 yB0,是 根据排放标准规定出的一个值。图3-14还出现一个点，即T E,时间到达T E时，吸附波整个移出床层， 说明床层内的吸附剂已完全饱和，完全失去了吸附能力，这一点称为耗竭点或 称干点，到达干点时，床层内流出的气体中，吸附质浓度基本回复到进口浓度。在实际操作中，一旦达到了破点，就应停止操作，切换到另一吸附床，穿 透了的吸附床转入脱附再生。（三）传质区高度把一个吸附波所占据的床层高度称为传质区高度，用Za表示。从理论上讲，传质区高度应是流出气体中溶质浓度从 0变到yo这个区间内吸附波在Z 轴上占据的长度，但实际上再生后的吸附剂中还残留一定量的吸附质（一般为初

9、始浓度yo的5%10%），而吸附剂完全达到饱和的时间又太长，所以一般 把由破点时间T B对应的气体浓度yB到干点时间T E对应的气体浓度yE这段时 间内吸附波在Z轴上所占据的长度称为传质区高度。为了使吸附操作比较可靠，就必须使床层有足够的长度，起码要包含一个稳定的传质区。而形成一个稳定的传质区需要一定时间。如果吸附器床层长度 比传质区长度还短，那就不能出现一个稳定的传质区， 操作不稳定，出现破点 的时间会比计算的来得快，为避免此点，吸附器床层长度一定要比传质区长度 长。例如实验室内所用吸附柱高度就规定应至少是传质区长度的两倍，而吸附柱直径最少应是最大吸附剂颗粒直径的 10倍。（四）传质区吸附饱

10、和率（度）和剩余饱和能力分率这两个概念可用下式表示：剩余饱和吸附能力分率Za内吸附剂仍具有吸附能 力 一 Za内吸附剂达到饱和时的 吸附能力吸附饱和率度=Za内吸附剂实际的吸附量Za内吸附剂达饱和时的吸_附量这也是量度固定吸附床操作性能的两个指标，吸附饱和率越大，剩余饱和 吸附能力分率越小，说明吸附床的操作性能越好。二、希洛夫近似计算法（一）希洛夫公式在理想状态下，在理想保护作用时间 T*内通过吸附床的吸附质将全部被 吸附，即通过床层的吸附质的量一定等于床层内所吸附的量，即：Gs b Ac。= ZA bXt（3-19）式中 Gs气体通过床层的速率，kg/（m2 s）; A 吸附床层截面积，m2

11、;xt吸附剂的静活性（平衡吸附量），kg/kg ；Tb理想保护作用时间，min;co气体中污染物初始浓度，kg/m3； p b吸附剂堆积密度，kg/m3；Z床层长度，m。由上式可得： :?bXtG s C0(3-20)15:，BXtGsC对于一定的吸附系统和操作条件，p B、XT、Gs、CO均已确定，因此可令(3-21)则（3-21）式可变成：Tb=KZ（3-22）但对一个实际的操作过程，由于床层存在阻力，因此实际上的保护作用时 间T B要比理想保护作用时间 TB短，我们把被缩短的这段时间称为保护作用 时间损失，用TO来表示。阻力越大，T O越大。三个时间的关系可表示如下:(3-23)T B=

12、 tP T o将(3-22)式代入上式，即得：B 二 KZ - o 二 K(Z Zo)(3-24)(3-24)式即为具有实用价值的希洛夫公式，Zo可以称为床层长度损失。T o和Zo均可由实验求得。(二)利用希洛夫公式的简化计算在吸附净化的设计中，常利用希洛夫公式进行简化计算。简化计算还是以 实验作基础，利用希洛夫公式求出 K与T。，再根据生产要求的操作周期求出 吸附床层长度，并根据气速，求出所需床层半径或截面积。具体步骤简述如下：1. 选择吸附剂，确定操作条件，包括温度、压力和流速。固定吸附床的气 体流速一般掌握在0.20.6m/s之间；2. 规定出合适的破点浓度；3. 在一定气速u下，测不同

13、床层长度Z的保护作用时间T B，作出T B-Z 直线，求出K和T 0 ；4. 定出操作周期T B，化为min;5. 将K、T 0、t b代入希洛夫公式，求出 乙若Z过长可以分层。6. 用下式计算床层直径：(3-25)(3-26)7、求吸附剂用量W:W式中 p b吸附剂堆积密度，kg/m3。为避免装填损失，可多取10%装填量。例3-1用活性炭固定床吸附器吸附净化废气。常温常压下废气流量为 1000m3/h，废气中四氯化碳初始浓度为 2000mg/m3，选定空床气速为20m/min。 活性炭平均粒径为3mm，堆积密度为450kg/m3,操作周期为40小时。在上述 条件下，进行动态吸附实验取得如下数

14、据：床层高度 Z/m0.10.150.20.250.30.35透过时间 t B/min109231310462550650请计算固定床吸附器的直径、高度和吸附剂用量。解：以Z为横坐标，T B为纵坐标将上述实验数据描绘在坐标图上得一直 线(例3-1附图)。依据图，求出直线的斜率即为K，截距即为-T 0,得K=2143 (min/m) t 0=95(min)将K、T 0、T B代入希洛夫公式得：40 60 952143=1.164(m)取Z=1.20mt采用立式园柱床进行吸附，计算出吸附床直径:鮎需罟o=1.03(m)可取D=1m。所需吸附剂量 2W 二AZ：、b=1.2 12 450 = 423

15、.9（kg）4考虑装填损失，所需吸附剂量 W为：423.9 x 1.1=466kg三、透过曲线计算法透过曲线计算方法与希洛夫近似计算法相比要复杂一些，但还是要假定吸附体系是一个很简单的恒温体系，混合气体中只有一种可被吸附的吸附质，该体系得到的仅有一个吸附波或传质区。此时固定床吸附器计算的主要内容为传 质区高度Za,保护作用时间t b和全床饱和度S。（一）传质区高度的确定图3-15为一理想透过曲线。气体的初始浓度为y（kg溶质/kg无溶质气体）， 气体流过床层的质量流速为Gs（kg/m2 h），经过一段时间后流出物总量为 W （kg无溶质气体/m3）。透过曲线是比较陡的，流出物中溶质的浓度从基本

16、上 为零迅速上升到进口浓度。以 yB作为破点的浓度，并认为流出物浓度升到接 近y0某一浓度值yE时，吸附剂基本上已耗竭。在破点处流出物量为Wb，而到吸附剂耗竭时，流出物的量为 We。这样，在透过曲线出现期间所积累的流 出物量Wa=WE-WB。把浓度由yB变化剂yE这部分的床层高度称为一个吸附区 或称传质区高度。当吸附波形成后，随着混合气体的不断通入，传质区沿床层不断移动，令 T a为吸附波移动一个传质区高度所需的时间，贝U：,e-wb）=W（3-27）GsGs又令T E为由通气开始至床层耗竭所需要的时间，即传质区形成和移出床层所 需的时间之和，贝U：业（3-28）E Gs设传质区形的时间为T

17、F，则T E- T F应是自吸附波形成开始到移出床层的时 间。在稳定操作时，当吸附波形成后，其前进的距离和所需要的时间之比（即 吸附波前进的速度）应是一个常数。设吸附床高度为乙传质区高度为Za，则:得出传质区高度为：ZaZa E BZa 二Z3-29)设气体在传质区中，从破点到床层完全耗竭所吸附的吸附质的量为U （kg/m2床层截面积），即为图3-15中阴影的面积weU = （y。-y）dWWb（3-30）若传质区中所有的吸附剂均为吸附质所饱和，则其吸附容量应为 yoWakg/m2。但实际情况是，当达到破点时，传质区内仍旧具有一部分吸附容 量，其值为U，式（3-30）所示。若以E代表到达破点时

18、传质区内仍具有的吸附能力与该区内吸附剂总的吸附能力之比，即前述剩余吸附能力分率，则3-31)UyoWa(yo -y)dWWbyWa很显然，（1-E）即代表了吸附区的饱和程度，E愈大，说明吸附区的饱和程度 愈低，形成传质区所需的时间愈短。当 E=0时，说明吸附波形成后，吸附区 内的吸附剂已全部达到饱和，此情况下，吸附形成的时间应与移动一个吸附波 长度的距离所需时间相等：T F= T a若E=1，即传质区内吸附剂基本上不含吸附质，则传质区形成的时间基本上等于零。据此两种极端情况，应有：=（1 -E）3-32）将式（3-32）代入式（3-29）得：3-33)因为T a=Wa/Gs, T E=WE/G

19、S，代入上式即得:Za =ZWaWe _(1_E)Wa (3-34)由式（3-34）可知，要确定传质区的高度Za，必须通过实验得出透过曲线 的形状，从而确定Wa、We和E的值。在实际吸附计算中，E 一般取0.40.6(二)吸附床饱和度设吸附床横截面积为A，吸附床高度为Z,其中吸附剂的堆积密度为p b, 则吸附剂的总量应为ZA p b。若床层全部被饱和，吸附剂与污染物进口浓度yo成平衡的静活性为xt,则此时吸附剂所吸附的污染物的量为：Q =ZA PrXt(3-35)实际操作中，达到破点时，总会有一部分吸附剂未达饱和，此时吸附床中 实际吸附量应为饱和区的吸附量与传质区的吸附量之和。其中：饱和区吸附

20、污染物的量=(Z-Za) A p bxt传质区吸附污染物的量=Za(1-E)A p bxt 于是整个吸附床的饱和度 S为：S _(Z -Za)ABXT Za(1 -E)ABXtZAPbXtZ -ZaE(3-36)Z式中 S假想的吸附剂流量 即yXt(三)传质区中传质单元数和传质单元高度的计算吸附操作过程中，随着吸附的进行，床层内的传质区沿气流方向移动，移 动的速度会远比气流通过的速度慢。为了分析问题的方便，假定传质区移动方 向与气流方向相反，则可把传质区认为是固定在某一高度，如图3-16所示。假设床层高度ZX，则在床层顶面气固相达到平衡状态。可对整个床层作物料衡算：GS(y。-0) =S(Xt

21、 -0)(3-37)(3-38)为操作线斜率(图3-16)。对于床层任上式可以看作是操作线方程, 一截面，则可有如下关系：GS y 二 Sx(3-39)设床层截面积为1,则对床层中微元高度dz作物料衡算，得:Gsdy =Kyap(y -y )dz(3-40)此物料衡算式表示的是单位时间单位面积的dz高度内，气体中溶质的减少量等于吸附剂固体中吸附的吸附质的量。式中 Gs气体流量，kg无溶质气体/m2 h；Kyap气相体积传质总系数，kg吸附质/ m3 h；y*与x成平衡的气相浓度，kg吸附质/kg无溶质气体。将(3-40)式整理并在传质区内积分，即得传质区高度：Za=KGs r dy =HogN

22、og(3-41)Kyap yB y -y上式中Hog、Nog可以称为传质区内的传质单元高度和传质单元数。与处理吸收计算相类似，传质单元数可用图解积分法求取。当平衡线接 近直线时，也可用下式近似计算：Nog=3Aym(3-42)式中 ym3-43)对数平均推动力-ym =(y1 _ y1 ) 一(y2 - y2 )Iny1 -yy2 *2对于低浓度气体，有时也可以用算术平均推动力。下面通过例题讲解有关计算。例3-2用硅胶固定吸附床净化含苯废气。废气初始浓度为y0=0.025kg苯/kg空气，操作温度为 T=298K，p=202.7kPa混合气体密度p v=2.38kg/m3, 动力粘度卩v=1.

23、8x 10-5kg/(m.s)。气流速度为1m/s，吸附周期为90min，破点 浓度yB=0.0025kg苯/kg空气，排放浓度yE=0.020kg苯/kg空气。硅胶堆积密 度p B=650kg/m3,平均粒径dp=6mm。比表面积ap=600m2/m3。给定条件下的 平衡关系为y*=0.167x15,传质单元高度1.42/ d pG S 0 .51。 Hog试计算床层高度。解：设床层截面积为1m2，以简化计算。废气流量 GS 二UgA：v =1 1 2.38kg/s ;传质单元高度Hog1.42 dpGs丁 FpV1.42600 (0.006 2381.8 1Q-、0.51)=0.071 m

24、根据平衡关系式y* =，当y0=0.025时，xt=0.282,过平衡线上该点B作操 作线，并按附表所列逐项计算列出。表中第1栏为yB和yE之间选取的y值；第2栏是自操作线上各点的y所对应的平衡线上y*的值。依次计算出之后， 用y值作横座标，以1/（ y-y*）作纵座标，绘出曲线，并在 yB与yE之间进行 图解积分，即可得到第5栏的值，即Nog的值。据此可得到对应于yE的Nog=5.925。于是得到传质区高度Za：将式（3-31）变换：E =We(yo y)dWWbyoWayyo)d(W -WbWaZa=HoG Nog=0.071 X 5.925=0.42(m) 例3-2附表*1/edyw _

25、wbyyyy-yy-yk y _yWWaV?1234567yB=0.00250.00090.0016625oo0.10.0050.00220.00283581.13750.1920.2o0.00420.00333041.90000.3210.30.0070.00630.00372702.61250.4410.450.00890.00362783.30000.5560.50.0100.01160.00342944.01250.6760.6o0.01480.00273704.83750.8150.70.0120.01800.00205005.92501.0000.850.015oo0.01755y

26、E=0.0200F面根据剩余吸附能力分率E的概念，计算吸附床层高度Z22W -Wb按照上式，若以y/yo为纵座标，以wa为横座标，绘出曲线，曲线与y/yo=1水平线、（W WB）/Wa=1的垂线之间的面积（见附图b），即为E。图解 积分可得：E=0.55。根据物料衡算关系式：(3-44)1 X Z p bSxt=Gs T Byo 将（3-36）式代入上式，得：ZbZ -乙EZXt =Gs Byo(3-45)将已知数代入上式，得:Z 650 Z 一0.42 0.55 0.282 =2.38 90 60 0.025解出 Z:Z=2.0(m)有了床层高度 乙 计算一下全床饱和度S:Z - ZaEZ2

27、 -0.42 0.552二 0.8845可见，全床饱和度接近90%，说明设计基本合理。四、经验估算法用吸附法净化气态污染物时会碰到多种情况，有时会缺乏前述理论计算时 所需要的数据，此时可由生产中或实验测得吸附剂的吸附容量值，用来估算吸附剂的用量，然后根据操作周期和经验气体流速(一般 0.20.6m/s),即可 计算出吸附床高度。例3-3拟用活性炭吸附器回收废气中所含的三氯乙烯。已知废气排放条件为294K，1.38XlO5Pa,废气中含三氯乙烯的体积分数为 2.0X10-3,流量为 12700m3/h,要求三氯乙烯的回收率为 99.5%。测得所要采用的活性炭对三氯 乙烯的吸附容量为 0.28kg

28、三氯乙烯蒸气/kg活性炭，活性炭的堆积密度为 576.7kg/m3，其吸附周期为4小时。操作气速根据经验，取 0.5m/s,求固定吸 附床咼度。解：三氯乙烯体积流量为：12700X 2.0X10-3=25.4m3/h。将此换算成标准状况下的体积：25.451.38 1051.01 1052732943= 32.226(m3/h)4D2700,3.14 0.5 3600计算得三氯乙烯摩尔质量为131.5，则得三氯乙烯质量流量为:131.5 32226 =189.18(kg/h)22.4在四小时内所要吸附的三氯乙烯的量为:189.18X 4X 0.995=752.94kg则所需活性炭体积为：752

29、.94/0.2834.66m 576.7按操作气速0.5m/s计，所需吸附床直径为:= 3.0(m)昇 222得床层截面积：A D2 =0.785 3 -7.065(m2) 于是得床层高度为：4 Z=m/A=4.66/7.065=0.66(m)五、固定床吸附器床层压降估算流体在固定床中的流动情况较之在空管中的流动要复杂得多，固定床中流体是在颗粒间的空隙中流动的，颗粒间空隙形成的孔道是弯曲且相互交错的， 孔道数和孔道截面积沿流向也在不断改变。 所以流体流过床层的压力降，主要 是由于流体与颗粒表面间的摩擦阻力和流体在孔道中的收缩、扩大和再分布等局部阻力引起，当流动状态为层流时，以摩擦阻力为主，当流动状态为湍流时， 以局部阻力为主；由于影响压力降的因素很多，目前尚无一个较完善的通用计 算公式。在设计固定床吸附器时，多根据实际情况，结合相关条件采用经验公 式进行计算，或采用实测数据，实测数据可从同类型的工业装置中获得， 也可 由一定规模的实验装置中获得。下面介绍几种可用于估算固定床吸附器的压力 降的经验公式。（一）气体通过静止吸附剂颗粒层的压力降当气体通过静止的吸附剂颗粒层时，由于床层内堆积的是大量的粒径和形 状不同的颗粒，颗粒之间的空隙结构毫无规则可言， 因而造成气体流动的通道曲折复杂，难以进行理论