1、1第四节 固定床吸附过程的计算固定床吸附器结构简单，但由于气体吸附过程是气固传质，对任一时间或任一颗粒来说都是不稳定过程，因此固定床吸附器的吸附操作是非稳态的，计算过程非常复杂，一般要涉及到物料衡算方程、吸附等温线方程和传热速率方程及热量衡算。而在气态污染物的吸附净化设计中，由于所涉及到的物系是低浓度的气态混合物，且气量一般比较大，吸附热相对较小，因此可近似地按等温过程处理，可不考虑传热速率方程和热量方程（升温脱附除外） 。这样在设计过程中可采用简化了的方法进行近似计算，计算时往往提出如下假设：（1）气相中吸附质浓度低；（2）吸附操作在等温下进行；（3）传质区通过整个床层时长度保持不变；（4）
2、床层长度比传质区长度大得多。这些简化限制条件对目前工业上应用的吸附器来说，一般是符合的。设计中较常采用的是希洛夫近似计算法和透过曲线计算法。计算过程一般是在吸附剂的选择、吸附设备的选择和吸附效率确定之后进行的。设计计算的任务是求出吸附器的床层直径和高度，吸附剂的用量，吸附器的一次循环工作时间，床层压降等。下面首先介绍固定床吸附器的吸附过程。一、固定床吸附器的吸附过程在固定床吸附器的吸附操作中，一般是混合气体从床层的一端进入，净化了的气体从床层的另一端排出。因此，首先吸附饱和的应是靠近进气口一端的吸附剂床层。随着吸附的进行，整个床层会逐渐被吸附质饱和，床层末端流出污染物，此时吸附应该停止，完成了
3、一个吸附过程。为了描述吸附过程，提出了以下概念。（一）吸附负荷曲线与透过曲线1. 吸附负荷曲线在实际操作中，对于一个固定床吸附器，气体以等速进入床层，气体中的吸附质就会按某种规律被吸附剂所吸附。吸附一定时间后，吸附质在吸附剂上就会有一定的浓度，我们把这一定的浓度称为该时刻的吸附负荷。如果把这一瞬间床层内不同截面上的吸附负荷对床层的长度（高度）作一条曲线，即得吸附负荷曲线。也就是说，吸附负荷曲线是吸附床层内吸附质浓度 x 随床层长度 z 变化的曲线。2在理想状态下，若床层完全没有阻力，吸附会在瞬间达到平衡，即吸附速率无穷大，则在床层内所有断面上的吸附负荷均为一个相同的值，吸附负荷曲线将是一个直角
4、形的折线，如图 3-12 所示。但实际上是不可能的，在实际操作中由于床层中存在着阻力，在某一瞬间床层内各个截面上的吸附负荷会有差异，我们这时所绘制的曲线将是图 3-13 所示的吸附负荷曲线。图中把曲线分成了三个区域：饱和区（所有吸附剂已经达到了饱和） 、传质区（有一部分吸附剂还正在吸附）和未用区（所有吸附剂上均未有吸附质） 。如果经过一段时间的吸附，绘制另一时刻的吸附负荷曲线时，会发现曲线前进到了 II线的位置，所以我们又形象地把吸附负荷曲线称为吸附波或吸附前沿。当吸附波的下端到达床层未端时，说明已有吸附质漏出，这时床层被穿透，当床层被穿透的这个时刻，称为破点。此时流出气体中吸附质的浓度称为破
5、点浓度。在实际工作中，由于吸附剂中吸附质的浓度（即吸附负荷）不易测定，故目前许多场合，曲线的纵坐标都以床层中混合气体的浓度 c 来表示。因此，吸附负荷曲线又可定义为在稳定吸附状态下，床层中气相中吸附质的浓度随床层高度（长度）变化的曲线。由于床层的阻力不同，吸附负荷曲线会有不同的形状。床层阻力愈大，某一时刻床层内各截面上浓度差别越大，吸附负荷曲线也就变得越平缓，这当然是我们不希望出现的情况。2. 透过曲线吸附负荷曲线表达了床层中浓度分布的情况，可直观地了解床层内操作的状况。但要从床层中各部位采出吸附剂样品进行分析是相当困难的，这样易破坏床层的稳定。因此通常改用在一定的时间间隔内。分析床层流出物中
6、吸附质浓度的变化，以流出物中吸附质浓度 y 为纵坐标。时间 为横坐标，则随时间的推移可画出一条 -y 曲线。如图 3-14 所示，开始时流出物中吸附质浓度为 yB，它是与吸附剂中的 xB 浓度相平衡的（x B 为破点时床层出口端的吸附负荷） 。流出物中吸附质浓度开始上升，到 E 时升到 yE,即接近床层进口浓度，这时床层已完全没有吸附能力，吸附波的末端也离开床层了。3于是在 -y 图上，从 B 到 E 呈现一个 S 型曲线，这条曲线称“透过曲线”。它的形状与吸附负荷曲线是完全相似的，只是方向相反。由于它与吸附负荷曲线成镜面对称相似，所以也称吸附负荷曲线为“吸附波”或“传质前沿” 。由于透过曲线
7、易于测定和标绘出来，因此也用它来反映床层内吸附负荷曲线的形状，而且也能准确地求出破点。如果透过曲线比较陡，说明吸附过程比较快，反之则速度较慢。如果透过曲线是一条竖直的直线，则说明吸附过程是飞快的，是理想的吸附波。（二）保护作用时间保护作用时间是固定床吸附器的有效工作时间。它定义为从吸附操作开始到床层被穿透所经历的时间称为保护作用时间，如图 3-14 所示的由 0 到 B 所经历的时间，到达 B 时，床层内吸附剂还没有完全饱和。图中的yB0，是根据排放标准规定出的一个值。图 3-14 还出现一个点，即 E，时间到达 E 时，吸附波整个移出床层，说明床层内的吸附剂已完全饱和，完全失去了吸附能力，这
8、一点称为耗竭点或称干点，到达干点时，床层内流出的气体中，吸附质浓度基本回复到进口浓度。在实际操作中，一旦达到了破点，就应停止操作，切换到另一吸附床，穿透了的吸附床转入脱附再生。（三）传质区高度把一个吸附波所占据的床层高度称为传质区高度，用 Za 表示。从理论上讲，传质区高度应是流出气体中溶质浓度从 0 变到 y0 这个区间内吸附波在 Z轴上占据的长度，但实际上再生后的吸附剂中还残留一定量的吸附质（一般为初始浓度 y0 的 5%10%） ，而吸附剂完全达到饱和的时间又太长，所以一4般把由破点时间 B 对应的气体浓度 yB 到干点时间 E 对应的气体浓度 yE 这段时间内吸附波在 Z 轴上所占据的
9、长度称为传质区高度。为了使吸附操作比较可靠，就必须使床层有足够的长度，起码要包含一个稳定的传质区。而形成一个稳定的传质区需要一定时间。如果吸附器床层长度比传质区长度还短，那就不能出现一个稳定的传质区，操作不稳定，出现破点的时间会比计算的来得快，为避免此点，吸附器床层长度一定要比传质区长度长。例如实验室内所用吸附柱高度就规定应至少是传质区长度的两倍，而吸附柱直径最少应是最大吸附剂颗粒直径的 10 倍。（四）传质区吸附饱和率（度）和剩余饱和能力分率这两个概念可用下式表示：这也是量度固定吸附床操作性能的两个指标，吸附饱和率越大，剩余饱和吸附能力分率越小，说明吸附床的操作性能越好。二、希洛夫近似计算法
10、（一）希洛夫公式在理想状态下，在理想保护作用时间 B 内通过吸附床的吸附质将全部被吸附，即通过床层的吸附质的量一定等于床层内所吸附的量，即：(3-19)式中 GS气体通过床层的速率， kg/(m2s)；A吸附床层截面积，m 2；xT吸附剂的静活性（平衡吸附量） ，kg/kg； B理想保护作用时间，min；c0气体中污染物初始浓度，kg/m 3； B吸附剂堆积密度，kg/m 3；Z床层长度，m。由上式可得：(3-20)对于一定的吸附系统和操作条件， B、x T、G s、co 均已确定，因此可令(3-21)则（3-21 ）式可变成： B=KZ (3-22)吸 附 能 力内 吸 附 剂 达 到 饱
11、和 时 的 力内 吸 附 剂 仍 具 有 吸 附 能剩 余 饱 和 吸 附 能 力 分 率 附 量内 吸 附 剂 达 饱 和 时 的 吸内 吸 附 剂 实 际 的 吸 附 量度吸 附 饱 和 率 ZaTBoBSxZAc/ ZcSTB0/KcS05但对一个实际的操作过程，由于床层存在阻力，因此实际上的保护作用时间 B 要比理想保护作用时间 B 短，我们把被缩短的这段时间称为保护作用时间损失，用 o 来表示。阻力越大， o 越大。三个时间的关系可表示如下： B= B o (3-23)将(3-22)式代入上式，即得：(3-24)（3-24 ）式即为具有实用价值的希洛夫公式，Z o 可以称为床层长度损
12、失。 o 和 Zo 均可由实验求得。（二）利用希洛夫公式的简化计算在吸附净化的设计中，常利用希洛夫公式进行简化计算。简化计算还是以实验作基础，利用希洛夫公式求出 K 与 o，再根据生产要求的操作周期求出吸附床层长度，并根据气速，求出所需床层半径或截面积。具体步骤简述如下：1. 选择吸附剂，确定操作条件，包括温度、压力和流速。固定吸附床的气体流速一般掌握在 0.20.6m/s 之间；2. 规定出合适的破点浓度；3. 在一定气速 u 下，测不同床层长度 Z 的保护作用时间 B，作出 B-Z直线，求出 K 和 0；4. 定出操作周期 B，化为 min；5. 将 K、 0、 B 代入希洛夫公式，求出
13、Z，若 Z 过长可以分层。6. 用下式计算床层直径：(3-25)7、求吸附剂用量 W：(3-26)式中 B吸附剂堆积密度，kg/m 3。为避免装填损失，可多取 10%装填量。例 3-1 用活性炭固定床吸附器吸附净化废气。常温常压下废气流量为1000m3/h，废气中四氯化碳初始浓度为 2000mg/m3，选定空床气速为20m/min。活性炭平均粒径为 3mm，堆积密度为 450kg/m3，操作周期为 40小时。在上述条件下，进行动态吸附实验取得如下数据：床层高度 Z/m 0.1 0.15 0.2 0.25 0.3 0.35透过时间 B/min 109 231 310 462 550 650请计算
14、固定床吸附器的直径、高度和吸附剂用量。解：以 Z 为横坐标， B 为纵坐标将上述实验数据描绘在坐标图上得一直线（例 3-1 附图） 。依据图，求出直线的斜率即为 K，截距即为- 0，得)(00ZKZuVD4BAZ6K=2143（min/m） 0=95(min)将 K、 0、 B 代入希洛夫公式得：取 Z=1.20m。 采 用 立 式 园 柱 床 进 行 吸 附 ， 计 算 出 吸 附 床 直 径：可取 D=1m。所需吸附剂量考虑装填损失，所需吸附剂量 W 为：423.91.1=466kg三、透过曲线计算法透过曲线计算方法与希洛夫近似计算法相比要复杂一些，但还是要假定吸附体系是一个很简单的恒温体
15、系，混合气体中只有一种可被吸附的吸附质，该体系得到的仅有一个吸附波或传质区。此时固定床吸附器计算的主要内容为传质区高度 Za,保护作用时间 B 和全床饱和度 S。（一）传质区高度的确定图 3-15 为一理想透过曲线。气体的初始浓度为 y0（kg 溶质/kg 无溶质气体） ，气体流过床层的质量流速为 Gs(kg/m2 h)，经过一段时间后流出物总量为 W（kg 无溶质气体/m 3） 。透过曲线是比较陡的，流出物中溶质的浓度从基本上为零迅速上升到进口浓度。以 yB 作为破点的浓度，并认为流出物浓度升到接近 y0 某一浓度值 yE 时，吸附剂基本上已耗竭。在破点处流出物量为WB，而到吸附剂耗竭时，流
16、出物的量为 WE。这样，在透过曲线出现期间所积累的流出物量 Wa=WE-WB。把浓度由 yB 变化剂 yE 这部分的床层高度称为一个吸附区或称传质区高度。当吸附波形成后，随着混合气体的不断通入，传质区沿床层不断移动，令 a 为吸附波移动一个传质区高度所需的时间，则：(3-27)又令 E 为由通气开始至床层耗竭所需要的时间，即传质区形成和移出床层所需的时间之和，则：（3-28）设传质区形的时间为 F，则 E- F 应是自吸附波形成开始到移出床层的时)(164.239504m)(0.1624uVD)(9.435014.2kgAZWBSaSBEaG)(SEGWBEaZ7间。在稳定操作时，当吸附波形成
17、后，其前进的距离和所需要的时间之比（即吸附波前进的速度）应是一个常数。设吸附床高度为 Z，传质区高度为Za，则：得出传质区高度为：(3-29)设气体在传质区中，从破点到床层完全耗竭所吸附的吸附质的量为U（kg/m 2 床层截面积） ，即为图 3-15 中阴影的面积(3-30)若传质区中所有的吸附剂均为吸附质所饱和，则其吸附容量应为yoWakg/m2。但实际情况是，当达到破点时，传质区内仍旧具有一部分吸附容量，其值为 U，式（3-30）所示。若以 E 代表到达破点时传质区内仍具有的吸附能力与该区内吸附剂总的吸附能力之比，即前述剩余吸附能力分率，则(3-31)很显然， （1-E ）即代表了吸附区的
18、饱和程度， E 愈大，说明吸附区的饱和程度愈低，形成传质区所需的时间愈短。当 E=0 时，说明吸附波形成后，吸附区内的吸附剂已全部达到饱和，此情况下，吸附形成的时间应与移动一个吸附波长度的距离所需时间相等： F= a若 E=1，即传质区内吸附剂基本上不含吸附质，则传质区形成的时间基本上等于零。据此两种极端情况，应有：(3-32)将式（3-32 ）代入式（3-29）得：(3-33)因为 a=Wa/GS， E=WE/GS，代入上式即得：FEaaZEBWdyU)(0aWaydyUB00)(aE)1(aaZ)1(8(3-34)由式（3-34 ）可知，要确定传质区的高度 Za，必须通过实验得出透过曲线的
19、形状，从而确定 Wa、W E 和 E 的值。在实际吸附计算中， E 一般取0.40.6。（二）吸附床饱和度设吸附床横截面积为 A，吸附床高度为 Z，其中吸附剂的堆积密度为 B，则吸附剂的总量应为 ZA B。若床层全部被饱和，吸附剂与污染物进口浓度 y0 成平衡的静活性为 xT，则此时吸附剂所吸附的污染物的量为：(3-35)实际操作中，达到破点时，总会有一部分吸附剂未达饱和，此时吸附床中实际吸附量应为饱和区的吸附量与传质区的吸附量之和。其中：饱和区吸附污染物的量=（Z-Z a）A BxT 传质区吸附污染物的量=Z a(1-E)A BxT 于是整个吸附床的饱和度 S 为：(3-36)（三）传质区中
20、传质单元数和传质单元高度的计算吸附操作过程中，随着吸附的进行，床层内的传质区沿气流方向移动，移动的速度会远比气流通过的速度慢。为了分析问题的方便，假定传质区移动方向与气流方向相反，则可把传质区认为是固定在某一高度，如图 3-16 所示。假设床层高度 Z，则在床层顶面气固相达到平衡状态。可对整个床层作物料衡算：(3-37)式中 S假想的吸附剂流量。即： EaaZ)1(TxQTtBaZxAE)1()(Sa)0()(0TSxyGTSxy09(3-38)上式可以看作是操作线方程， 为操作线斜率（图 3-16） 。对于床层任一截面，则可有如下关系：(3-39)设床层截面积为 1，则对床层中微元高度 dz
21、 作物料衡算，得：(3-40)此物料衡算式表示的是单位时间单位面积的 dz 高度内，气体中溶质的减少量等于吸附剂固体中吸附的吸附质的量。式中 GS气体流量， kg 无溶质气体/m 2h；Kyap气相体积传质总系数，kg 吸附质/ m 3h；y* 与 x 成平衡的气相浓度，kg 吸附质/kg 无溶质气体。将（3-40 ）式整理并在传质区内积分，即得传质区高度：(3-41)上式中 HOG、N OG 可以称为传质区内的传质单元高度和传质单元数。与处理吸收计算相类似，传质单元数可用图解积分法求取。当平衡线接近直线时，也可用下式近似计算：(3-42)式中 y m对数平均推动力(3-43)对于低浓度气体，
22、有时也可以用算术平均推动力。下面通过例题讲解有关计算。例 3-2 用硅胶固定吸附床净化含苯废气。废气初始浓度为 y0=0.025kg苯/kg 空气，操作温度为 T=298K，p=202.7kPa ，混合气体密度 v=2.38kg/m3，动力粘度 V=1.810-5kg/(m.s)。气流速度为 1m/s，吸附周期为 90min，破点浓度 yB=0.0025kg 苯/kg 空气，排放浓度 yE=0.020kg 苯/kg空气。硅胶堆积密度 B=650kg/m3，平均粒径 dp=6mm。比表面积ap=600m2/m3。给定条件下的平衡关系为 y*=0.167x1.5，传质单元高度 SGxydzaKdp
23、yS)(OGypSaNHdKGZEBmOGyN21211ln)()(yym51.0)(42.1VSpOGaH10。试计算床层高度。解：设床层截面积为 1m2，以简化计算。废气流量 ；传质单元高度根据平衡关系式 y*=0.671x1.5 绘出吸附等温线（如附图 a） 。由平衡关系可知，当 y0=0.025 时，x T=0.282，过平衡线上该点 B 作操作线，并按附表所列逐项计算列出。表中第 1 栏为 yB 和 yE 之间选取的 y 值；第 2 栏是自操作线上各点的 y 所对应的平衡线上 y*的值。依次计算出之后，用 y 值作横座标，以 1/（y-y *）作纵座标，绘出曲线，并在 yB 与 yE
24、 之间进行图解积分，即可得到第 5 栏的值，即 NOG 的值。据此可得到对应于 yE 的 NOG=5.925。于是得到传质区高度 Za：Za=HOGNOG=0.0715.925=0.42(m)例 3-2 附表y y* y-y*1 2 3 4 5 6 7yB=0.00250.00500.00750.01000.01250.01500.0175yE=0.02000.00090.00220.00420.00630.00890.01160.01480.01800.00160.00280.00330.00370.00360.00340.00270.0020625358304270278294370500
25、01.13751.90002.61253.30004.01254.83755.925000.1920.3210.4410.5560.6760.8151.0000.10.20.30.40.50.60.70.8下面根据剩余吸附能力分率 E 的概念，计算吸附床层高度 Z。将式（3-31 ）变换：按照上式，若以 y/y0 为纵座标，以 为横座标，绘出曲线，曲线与 y/y0=1 水平线、 (WW B)/Wa=1 的垂线之间的面积（见附图 b） ，即为 E。图解积分可得：E=0.55。根据物料衡算关系式：1Z BSxT=GS By0 (3-44)将（3-36 ）式代入上式，得：skgAuGVgS/38.1 mGdaHVSpOG 071.)8.1306(42.)(451.51.0 )()1()(00aBaWWdyydB aB0GxEBsTaBy1ydEBABW0y
