1、第三章 吸收（传质理论之一）1 概述 一、 传质：物质从一相转移到另一相的过程。如：吸收、蒸馏、萃取、结晶、吸附、蒸发、干燥、冷冻、升华、凝华等。二、吸收：用某种液体(吸收剂S)处理气体混合物，使其中某个目标组分(吸收质A)得到分离的过程，不能被吸收的组分统称为惰性组分B。吸收 解吸 互为逆过程。三、吸收的应用：制取产品，净化气体，回收有用组分。四、吸收分类：,利用混合气体中各组分溶解度的差异分离气体,物理吸收:,特点：弱结合,游离,可逆易解吸 速率：g l扩散速率 极限：达g l相平衡,化学吸收：AS发生化学反应而分离气体,五、吸收流程： 在连续操作的立式吸收塔中，气，液逆流接触进行传质。1
2、 要使气相中A被吸收完全-净化气体：则采用单程(单塔或多塔)吸收。2 希望液相A的浓度高-制取产品：则采用液相循环吸收。,特点：强结合，不可逆 速率：反应或扩散速率 极限：化学平衡,(本章只讨论物理吸收),3 要分离和回收A：则采用吸收解吸 联合流程。,131-16,六、收吸剂的选取原则： 高选择性、低蒸气压、低粘度、低价格。,喷淋塔： 气相连续 液膜塔(填料塔)：气相连续鼓泡塔： 液相连续,4 吸收塔,吸收 解吸 联合流程,洗油,苯,七、要解决的问题： 1. 吸收极限相平衡理论 2. 吸收速率双膜理论。 3. 填料塔吸收面积的计算。4. 强化吸收的途径。八、吸收体系溶解度,易溶：NH3, S
3、O3, NO2, HClH2O中等： SO2 , Cl2 , H2S H2O 难溶：O2, H2, CO, CO2 H2O,2 吸收相平衡 当物理吸收达到两相平衡A(g) A(l)时，溶质A在气液两相中的浓度不变，且服从亨利定律关系。,NA(gl)=NA(lg) (NAmolA/s) pAxA？,亨利定律在温度一定总压不太高(5 atm)时, 大多数中等或难溶体系吸收平衡时的液相浓度与其气相分压成正比。,441-12,气液两相浓度有多种表示法，亨利定律有多种形式。1 以pA xA表示的亨利定律 pA*=E xA 或 xA*=pA/ E 即 xA*pA pA*A的气相平衡分压(与浓度为xA的液相
4、平衡)， xA液相A的摩尔分率 (实际浓度)。 xA*A的液相平衡浓度(与分压为pA的气相平衡)， pA气相A的分压 (实际浓度) 。 E亨利系数Pa或atm, E, xA, 即吸收体系的亨利系数越大, 气体越难溶。附表九可查E值。,2 以pA CA表示的亨利定律 p*A= CA/ H 或 C*A= H pA 即 C*A pA CA液相A的物质的量浓度(实际浓度)kmol/m3. C*A液相A的平衡浓度(与分压为pA的气相平衡),541-16,H溶解度系数kmol/m3*Pa or kmol/m3*atm, E, CA,即吸收体系的溶解度系数越大，气体越易溶。 对一定的吸收体系(质A, 剂S确
5、定)，E,H为常数，都反映A在S中的溶解度大小。EH关系？先找CAxA关系：,即E1/ H ，E，H，气体在该吸收剂中越难溶。,xACA关系：,(),3 以yA xA表示的亨利定律A的气相摩尔分率 yA=pA/P = nA/ n总 y*A=p*A/P=(ExA)/P=mxA 或 x*A= yA/m 即 x*AyA m=E/P相平衡系数(无因次)，E, m,气体越难溶。,例51解：空气一水体系，已知p空气=P总和E，求液相平衡浓度(以x*A 和cm3/100g水溶液表示)x*A=pA/E=101325/6.73*109=1.51*10-5 (molA/mol溶液) 即水中空气含量0.00151
6、mol空气/100 mol H2O即100mol(1800g)水中含空气0.00151mol 0.0015*22.4*103cm3=33.6 cm3(标态)故水中空气含量33.6 cm3/1800g=1.86(标态) cm3/100g水,241-12,例5-2解：SO2-水体系,20时溶液的浓度为2.5%(质量百分比, 一般题给气相或液相浓度需要换算),求气相SO2的平衡分压,查附表九知：H= 0.36107 Pa将质量百分比m%xA,p*A=E*xA=0.36*107*0.00716=25.78 kPa (节1),例5-3：解含NH320%(yA1)的NH3-空气混合气体100m3, 用水恒
7、温吸收至混合体中含氨5%(yA2),求NH3被吸收的体积。 解：由于吸收前后混合气体总体积变化，以吸收前混合气体总体积100m3为基准计算：VNH3吸收=VNH3前-VNH3残吸收前：VNH3前=V总前*20 %(yA1)=100*20% = 20 m3 吸收后: VNH3残=V总后(无法求得)*5% (yA2)= VB*(VNH3残/VB) 而: VB =V空气=100-20 =80 m3 (VB吸收前后不变且易求),141-15,若能求出每份空气含NH3的分数VNH3残/ V空气，则吸收后残存的VNH3残可求。 而 VNH3残/ V空气= (VNH3残/V后总)/ (V空气/V后总)=yA
8、2/ yB2 =yA2/(1-yA2)=5%/(1-5%)=5/95 可求 (份NH3/份空气) VNH3残=5/95*80=4.2 m3 VNH3吸收=20-4.2=15.8m3,上例中利用惰性气体量在吸收前后不变的特点，求出吸收后惰性气体中吸收质A含量分率，再乘以惰气的量(vB-m3B, or nB-molB)，即可求吸收后气体中残存A的量，进而求A被吸收的量。,惰性气体中吸收质A含量分率： vA/vB=nA/nB=(nA/n总)/(nB/n总)=yA/yB= yA/(1-yA),吸收前1：nA1/nB=(nA1/n总1)/(nB/n总1)= yA1/yB1 =YA1 = yA1/(1-y
9、A1)吸收后2：nA2/nB=(nA2/n总2)/(nB/n总2)= yA2 /yB2 =YA2 = yA2/(1-yA2)由此提示我们：要求A被吸收的量，要用到A在惰性气体B中A的含量分率yA/yB，令为YA。 定义：比摩尔分率浓度,YA:每mol惰气中所含A的mol数。起始时-YA1,终了时-YA2,4 以YAXA比摩尔分率表示的亨利定律,XA：每mol吸收剂中所含A的mol数。吸收从XA1XA2,即:,或:,上例用比摩尔分率计算： 吸收前： YA1= yA1/yB1=0.20/(1-0.20)=0.25 吸收后：YA2=yA2/yB2=yA2/(1-yA2)=0.05/0.95=0.05
10、3 被吸收NH3的体积： VNH3VB(YA1-YA2)= 80*(0.25-0.053) =15.8 m3 331-9,YA*=mXA,XA*=YA/m,注意：1 吸收达平衡时：YA*=mXA, XA*= YA/m，但两方程的意 义不同，YA*与YA不同。2 yA=nA/nT YA=yA/yB=nA/nB xA=nA/nT XA=xA/xS=nA/nS 3 溶解度mAg(A)/100g(s)xA，CA，XA关系： xA= nA/ nT= (mA/MA)/( mA/MA+100/MS) XA=nA / nB= (mA/MA)/(100/MS) CA= nA/VT= 液* mA/(100+mA)
11、/MA kmol/m3 4 长期放置的吸收体系是平衡体系，在平衡体系中引入气 体或溶液, 平衡将被破坏, 可能再吸收(pAp*A或xAx*A)。,例5-4 已知总压P=101.325 kPa及 PAmA (mAg(A)/1000g(s)关系，求YAXA关系： 解：先转变成 yAxA关系： yA=nA/n总= pA/ P xA=nA/n总=(mA/MA)/( mA/MA+1000/MS) 再转化成YAXA关系： YA= yA/(1- yA) XA= xA/(1- xA)亦可直接求：YA=pA/pB=pA/P-pA XA=nA/nS=(mA/MA)/(1000/MS),441-13,5-3 吸收速
12、率 吸收速率决定吸收达到平衡的时间，决定吸收操作的生产强度，是吸收设备选型和设备设计的重要依据。,二、扩散1 分子扩散：物质以分子热运动方式穿过静止或滞流流体的传递过程特点：传递速率慢。2 对流扩散：物质以相对运动方式穿过湍流流体的传递过程特点：传递速率快。,一、吸收速率定义：NA= dnA/d对于稳定吸收过程：NA=nA/ mol(A)/s 吸收过程是物质的相际转移过程，通过扩散方式进行。,三、费克分子扩散定律 如图:,四、吸收机理双膜理论,单位时间内穿过面积为A，厚度为的扩散层的吸收质A的物质的量,D扩散系数 m2/s,实际吸收过程一般在喷淋塔、填料塔、鼓泡塔中进行。无论采用那种吸收塔，吸
13、收过程都存在气液相界面，类比间壁传热提出了双膜理论以导出吸收速率。 如图：,分子扩散模型,将对流给热视为壁附近滞流层的热传导过程付立叶定率,将吸收视为A穿过相界面附近滞流双膜的分子扩散过程费克定律 331-18,传热过程,吸收过程,理论实质,理论模型,实际情况,1-高温流体传热膜(分)系数2-低温流体传热膜(分)系数,kg-气膜吸收分系数kl -液膜吸收分系数,(节2),导出速率,给热方程,气膜、液膜分速率方程,1 气液两相主体因湍流使A以对流扩散方式传递，传递速率快, 浓度一致，忽略湍流主体中对流扩散的阻力。2 相界面两侧存在滞流的气膜和液膜，双膜中A以分子扩散方式传递，传递速率慢，浓度下降
14、很快，整个传质过程的阻力全部集中在双膜中。,3 到达相界面的A能无阻力地瞬 时溶于液相，且界面上气液两 相达平衡, 即：pi = ci / H 。,双膜理论要点：,141-15,五、费克分子扩散定律的应用当A分子通过分子热运动的碰撞而传递时，其传递速率：,g ,l 气膜、液膜厚度, 由于在膜层中有部分对流传质, 真正的分子扩散距离g,l。对气膜：用pB/P对g加以较正：,=kgA(p-pi),对液膜：用CS /(CA+CS)=CS /CT 对l 加以较正：,=klA(Ci-C),所以，可用界面附近气膜中的扩散速率：NA=kgA(p-pi)或液膜中的扩散速率：NA=klA(Ci-C) 计算吸收速
15、率。,(材三),kg、ky分别以分压、气相摩尔分率表示的气膜吸收分系数,2 液膜吸收分速率方程.,kl、kx分别以物质的量浓度、摩尔分率表示的液膜吸收分系数,六、吸收速率方程1 气膜吸收分速率方程,或表示为:,或表示为：,3 气相吸收总速率方程pi、Ci 、yi、xi 难求,没法消去,p-pi 气相主体与界面上A的分压差。ci-c 界面上与液相主体中A的浓度差。,p*是与液相主体浓度C平衡的分压,吸收过程中，气相浓度，用某处的实际分压与该气相能下降到的最低值之差表示推动力，而气相所能降到的最低分压是与该气相相接触的液相达平衡的平衡分压pA*。如某处的气液相主体浓度为pACA，用pA-pA*表示
16、推动力，而pA*CA平衡，这样增大，阻力也增大（双膜的阻力）。,541-17,NA=kgA(pA-pi)=klA(ci-cA)=klAH(pi-pA),H,H,H,H,KG,Ky,KY -以分压Pa,摩尔分率xA,比摩尔分率YA表示推动力的气相吸收总系数。 4-14,1/KG = 1/kg + 1/klH,4 液相吸收总速率方程,KL,Kx,KX以c, x, X表示推动力的液相吸总系数 (1-13,3-12),注意：以不同浓度表示推动力，其数值和单位各不相同，但NA总是相同的，单位molA/s，通过K、k (吸收总系数和分系数)来调节NA的数值和单位。,七、强化传质吸收的途径和措施 吸收速率
17、NA=KA1 提高K值,即：总阻力=气膜阻力+液膜阻力,1/KG=1/kg+1/klH,1/KL=H/kg+1/kl,KG/KL=H,331-19,H,H,H,对难溶体系：H小，m大， 则： 1/KL 1/kl，1/Kx1/kx 即：总阻力液膜阻力液膜控制 应选用鼓泡塔吸收，以搅拌液体，使液膜阻力。 又k1/，故加大气液两相的相对速度，Re，k，亦可提高K值。,对易溶体系： H大，m小， 则： 1/KG1/kg，1/Ky1/ky 即：总阻力气膜阻力气膜控制 应选用填料塔或喷淋塔吸收，以搅拌气体，使气膜阻力。,3 增大气液接触面积A、增大气液两相的分散度。、开发和采用比表面积（m2/m3）大的填
18、料。,3-13,2 增大吸收推动力 采用逆流吸收。 若=y-y* or: =Y-Y* y和Y恒定，可减小y*、Y* y*= mx, Y*= m X， 增大吸收剂用量，则x和X，y*和Y*都，。 又 m=E/P=S/PHMS P; 使用MS大的吸收剂，可使m； T，H，m，。,541-18,要求 完成一定吸收任务所需传质面积要用NA=KA求A，必须先求：NA、K、(Y-Y*)X1L。一、填料塔的物料衡算式操作线方程 如图：GmolB/s惰性气体的摩尔流量 （工艺条件所确定）,LmolS/s吸收剂的摩尔流量(人为确定)在稳定吸收过程中：对全塔进行物料衡算得： NA=G(Y1-Y2)=L(X1-X2
19、) .(1) 可求NA 式中：G, Y1, Y2, X2，均为定值，只X1L变。 从塔顶至任一截面M-M作物料衡算得：,54填料吸收塔的计算,141-17,上式是一条直线方程, 在YX相图上为一直线, y=kx+b两端点为(X1,Y1), (X2,Y2) 。斜率L/G液气比。 由于(X1,Y1), (X2,Y2)分别是塔底和塔顶的气液相实际组成，所以，直线上任意一点的坐标都表示塔中某一截面上实际的气液组成操作线方程。,二、L/GX1关系 将吸收体系的相平衡线与操作线，绘在同一XY相图上,平衡线：,操作线：,(过(X1, Y1)，(X2, Y2)两点),气相推动力Y-Y*塔顶 Y2-Y2*;塔底
20、Y1-Y1*,X,Y,X1,液相推动力X*-X 塔顶X2*-X2 塔底X1*-X1,实际上仅由 L 决定：L,L/G, 图, , X1; L, L/G, ，X1 图中 。,即：吸收推动力大小取决于操作线偏离平衡线的程度，由 操作线斜率 L/G 确定。,是操作线到平衡线的纵向距离， 从Y2Y1。,是平衡线到操作线的横向距离，从X2X1。,对操作线方程变形有：,式中G, Y1, Y2, X2, 均由工艺条件确定, 只有X1L变。,L，L/G ，(Y-Y*)和 (X*-X) ，图中线Y1、X1都。,L, L/G , (Y-Y*)和(X*-X),图中线Y1 ，X1都，但X1 。,三、吸收剂用量L的确定
21、,即：X1 和由 L而定，如何确定L？,已达平衡X1X1*, 塔下部无用。即：当X1X1*时：L=LM最小吸收剂用量， LM/G最小液气比 (图中线) 。,L再，L/G再 ，操作线与平衡线相交(切)，则：(Y-Y*)和(X*-X)下降到 0, X1到X1*(=Y1/m) , 达最大值; 而 L下降到(Lmin=)LM。 L再 , 则: 不到塔底, 气液,即：,而X1*=Y1/m，Y1、Y2、X2、G、 m已知，所以LM 或 LM/G 可求。,331-20,实际吸收剂用量,四、液相出口浓度X1的计算,知 L可求X1。,求X1和L步骤小结：,据：,或：,4-15,1 求,2 求X1*： X1*=Y
22、1/m 而m=E/P=CT/PH=S/MSPH 或： 据Y1到平衡线上查得X1* (只有气液平衡数据)。,若qv混合气体体积流量Tm3/h，则：,541-19,3 求最小液气比 LM/G：,4 求实际液气比L/G： L/G=(1.22.0)LM/G 求吸收剂用量：L= L/G*G=(1.22.0)LM/G*G,5 求X1： 据,可求X1。,若不需求LmolS/s，可据,例：5-5 、5-6 、P161,求X1。,作业P186 18 (Y*=26.7X),331-21 141-18 541-19,五、吸收推动力的计算,由前面的讨论知： NA=KYA(Y-Y*)=KXA(X*-X)塔底处气液两相推
23、动力：Y1=Y1-Y1* 和 X1=X1*-X1塔顶处气液两相推动力：Y2=Y2-Y2* 和 X2=X2*-X2,类似于求tm= (Tt )m，先找 d(T-t) dA关系，找d(Y-Y*)dA，然后积分。, (Y1-Y1*)(Y2-Y2*) (X1*-X1)(X2*-X2)故须取平均值，对服从亨利定律的吸收体系：以(Y-Y*)为例求： (Y-Y*)m= Ym,X1,Y2,X2,Y1,如图：在吸收塔中取一微元高度dH，dH内填料的表面积dA, 在dH内作物料衡算：衡算方向：塔顶塔底( 浓度),底,底,顶,化5,上式积分得：,对应,代入上式得：,注意： 1表示塔底的脚标，2表示塔顶的脚标。 计算
24、Ym时，Y1, Y2, Y2*, m, X2 易求, 关键在于 求Y1*，要 求Y1* X1L/GLM/GX1*Y1。 计算Xm时, X1*=Y1/m, X2*=Y2/m, X2=x2/1-x2 易求，关键在于求X1，求法同上。,或直接用公式求：,34,54-19,2-13,1 以气相浓度表示的计算式：,又A=H填料*A0* A0填料塔横截面积，填料的有效比表面积，可近似用代替，填料的比表面积。,气相传质单元高度,气相传质单元高度，是气相吸收质的浓度变化量刚好等于过程平均推动力 Ym时所需传质区域对应的填料层高度,是使浓度变化Y1-Y=Ym所需H填料,六、填料层高度的计算式：,2 液相表示计算
25、式,液相传质单元高度,气相传质单元数是经全塔吸收后气相浓度的变化量相当于过程平均推动力的倍数。,141-19,2解：y1=10/17/1000/22.4=22.4/1700=0.01318 =0.99 L=1.5LM Y*=1.2X X2 =0 G=5000*(1-0.01318)/22.4=220.3 kmol/hY1= y1/ (1- y1)=0.01318/(1- 0.01318)=0.013356 Y2=Y1 (1-)=0.0001336 X1*=Y1/m=0.013356/1.2=0.010984 (L=1.5LM/G*G=?) X1 =X1*/1.5=0.010984/1.5=0.
26、007323 Y1*=m* X1=1.2*0.007323=0.008787 Y2*=m*X2=0Y1=0.013356-0.008787=0.004569 Y2 =0.0001336-0=0.0001336 Y1/Y2=34.2,2. 某吸收塔在常压下用清水吸收焦炉气中的氨。焦炉气在标准状态下含氨10g/m3，焦炉气的处理量为5000m3/h(标态)，吸收率不低于99%。水用量为最小用量的1.5倍。混合气体303K时进入吸收塔，在此条件下的平衡关系为：Y*= 1.2 X，试计算气相传质单元数NOG。(塔径为1.4 m,气相吸收总系数KY=220 kmol / (m3*h) )，并求填料层高度
27、H填料。,3、不服从亨利定律的体系填料层高度计算 前面用过 dNA= G dY = L dX = L d(Y*/m),对不服从亨利定律的体系,则无法寻得d(Y-Y*)dNA 关系, 但仍有：,Ym=Y1-Y2/ln(Y1/Y2)=(0.004569-0.0001336)/ln34.2=0.004435/3.5322=0.001256 nG=Y1-Y2/Ym= (0.01336-0.0001336)/0.001256 =0.01322/0.001256=10.53 (个) hG=G/KYA0=220.3/220*0.785*1.42=0.6508 m H填料= hG*nG=0.6508*10.5
28、3=6.85 m,步骤：题给平衡数据换算为XY数据，在XY坐标图上作出平衡线。 求出Y1, Y2, X2, 据X2查得Y2*, 据Y1查得X1*LM/G L/GX1，据(X1, Y1), (X2, Y2)作出操作线。 在Y1Y2之间取一组Y值，对每一个Y值，在操作线上查得对应的X值，据X到平衡线上查得对应的Y*，计算Y-Y*, 并求得1/Y-Y*数据, 列表：,以Y对1/(Y-Y*)作图得一曲线，曲线下方Y1Y2之间的面积即为nG。,同理，对液相有：,液相传质单元数,同理：在X1X2之间取一组X值, 据X到操作线查得Y,再据Y值到平衡线查得X*值，计算X*-X值和1/X*-X值，列表，作图可求
29、：,材6,习题课： P191 24 26例：在逆流吸收塔中，用洗油吸收焦炉气中的苯，已知焦炉气进气流量为1000 m3/h（常压25吸收），其中含苯3%，要求 =95%，洗油中原来含苯0.005 %吸收条件下Y*=0.113X, L=1.2LM。求每小时耗油洗多少 nG,化20,解:Y1=y1/1-y1=0.03/1-0.03=0.03093 Y2=Y1(1-)=?x2=0.00005 X2x2 Y2*=0.113X20, X1*=Y1/0.113=0.0309/0.113=0.2734 LM/G=Y1-Y2/X1*-X2=0.1074, G=1000*(1-y1)*273/298/22.4=
30、39.67 kmol/h, L=1.2LM= 1.2(LM/G)*G=1.2*0.1074*39.67=5.113 kmol/h, X1=(X1*-X2)/1.2+X2=0.2277, Y1*=X1*0.113=0.02573, Ym=(Y1-Y2)/ln(Y1/Y2)=0.00310,习21:含乙炔y1=0.20的乙炔-空气混合气体在25和101.3kPa下与乙炔水溶液相接触,溶液浓度为0.3g (乙炔)/1000g(水).相平衡关系为:p=0.135106x kPa求传质方向和传质推动力Ym,解：实际浓度与平衡浓度之差决定传质方向，如： pAp*A or：xA x*A则继续吸收，反之则解吸。 pC2H2=101.30.2=20.26 kPa 而p*C2H2=0.135106 xC2H2xC2H2=0.3/26/(0.3/26+1000/18)=0.0002077 p*C2H2=28.04 kPa pA pA * 解吸 或:= 0.0002077 而=pC2H2/0.135106=0.00015 xC2H2 x*C2H2 解吸 Y=pC2H2/p空气=20.26/(101.3-20.26)=0.25 Y*=p*C2H2/p空气=28.04/101.3-28.04=0.3828 Y=Y-Y*= 0.25-0.3828= - 0.133 331-22,
