分离工程-II--06-1-多级分离计算一逐板法

1、第六章 多组分多级分离的严格计算 逐板计算法,平衡级分离过程种类繁多，分离设备结构各异，但主要受以下关系式所制约： 即： 物料平衡关系式(M)， 热量平衡或焓平衡关系式(H)， 相平衡关系式(E) 摩尔分率加和式(S)， 也就是说，一个平衡级分离过程可用上述四种关系式(MSEH)描述。,平衡级分离的严格计算，就是根据已知条件求解上述关系式组成的方程组，获得未知诸变量问题。上述方程组是非线性的，特别是当用很多平衡级分离非理想性很强的多组分时，所建立的方程组求解起来非常复杂，所以，如何求解方程组是多年来许多学者持续研究的课题。目前已经开发出几种适用于各种情况的多组分多级分离的计算方法，大体可以解决

2、当前化工分离中的计算问题。,多级分离计算的类型，按求解问题的特点可分为设计型和操作型两种。设计型计算是从一定的分离要求出发，求所需要的平衡级数；操作型计算是从给定的平衡级出发求，在一定条件下的分离性能。 我们所讲的多组分多级分离的严格计算主要以研究操作型计算为主。 首先，我们来学习多组分多级分离的逐级计算法。这种方法既有设计型算法，也可用于操作型算法。,6.1 多组分精馏的逐级计算法,逐级计算法用于双组分多级分离的计算，化工原理中已作过详细介绍。采用逐级计算用于多组分多级分离过程时，必须逐个组分从顶端或底端开始，交替利用相平衡关系及物料平衡逐级计算，直到某平衡级的组成满足预定要求为止，从确定平

3、衡级数。由于非关键组分的浓度分布难以准确确定，采用逐级计算十分繁琐，因而对于多组分多级分离过程，主要研究操作型计算方法。,5,xF,6.1.1 设计型算法,逐级计算法作为设计型算法，在用于双组分多级分离计算是很成功的，梯级做图法就属于此法。但对于多组分逐级计算，在计算终点处，我们需要对非关键组分浓度的合理分布进行反复试算。作为设计变量，在多组分分离中只规定了两个轻重关键组分的分离要求，由此，仅能定出塔底塔顶产品中两关键组分的浓度，其他非关键组分的浓度都是未知的，不可任意规定。然而，逐级计算要求先假定上下任一端的全部组分浓度，以便逐级向另一端计算。,然而，对非关键组分浓度的任何不准确的假设，都会

4、导致计算结果产生偏差而不能契合。因此，逐级计算法用于多组分分离时，只有当分离物系中除轻重关键组分外，只含有比重关键组分更重的重组分，而不含或几乎不含较轻关键组分更轻的轻组分时，则塔底产品中除规定的轻重关键组分之外，必定几乎集中了全部的重组分，这样，我们设计塔底组成才会相当准确，逐级计算至塔顶，不致有过大的偏差，从而大大地减少试算量，甚至不必试差也可得到比较满意的结果。,在计算既有轻组分，又有重组分的分离物系时，常采用从两端向中间逐级计算，当计算至轻、重关键组分的计算值于进料相符时，还需考察非关键组分计算值是否于进料实际组成契合一致，以判明原设定的塔顶塔底组成是否恰当，契合不好时，要重新设定，重

5、复计算，直至完全契合。 它的计算步骤时，先确定各组分在塔顶、塔底的分配，然后分别从塔底及塔顶逐级计算到进料级，如果气相或液相的关键组分之比与与原料之比相等，则可估计所需平衡级数及进料级位置。,逐级计算法用于恒沸精馏及萃取精馏是有效的，双组分溶液进行恒沸精馏或萃取精馏时须加入夹带剂或溶剂，属于三组分物系，不会同时存在轻重非关键组分，可以用逐级计算法进行计算。,6.1.2 操作型算法,采用操作型算法计算多组分多级分离过程，首先规定平衡级数及进料级位置，计算可能的分离程度。计算时，先假定平衡级各级温度（初值），从塔的顶底两端分别向进料级计算组成，校正温度改变两端产品组成，直到满足进料级契合要求。它的

6、特点是先规定级数及进料位置，用于已有的精馏塔进行查定考核是方便的。当 用此法做精馏塔设计计算时，需用简洁法估算平衡级数及进料位置，才能做逐级计算，由计算结果是否满足设计要求，判断原来固定的级数是否合适，必要时，需重新设定级数，重新计算，直至满足要求。,操作型逐级计算公式中诸变量常以组分流率表示，各变量定义如下： ln,i=Lnxn,i 第n级i组分液相流率 vn,i=Vnyn,i 第n级i组分汽相流率 fi=Fzi 进料中i组分流率 di=Dxdi 塔顶产品中i组分流率 wi=Wxwi 塔底产品中i组分流率 An,i=Ln/(Vnkn,i) 吸收因子（n级，i组组分） Sn,i=Vnkn,i/

7、Ln 蒸出因子（n级，i组组分）,将组分流率定义式代入相平衡关系式yni=knixni中消去yni及xni可导出以组分流率为的变量的相平衡关系式： ln,i=An,ivn,i vn,i=Sn,iln,i 逐级计算公式中以比流率为变量，即以某组分在各级上的比流率与该组分在塔顶或塔底出料流率之比为变量，因而又称为比流率法。,1). 直接迭代法,取塔级排序自下而上，以塔底再沸器为第1级，顺次向上，进料级序号为M, 塔顶第1个平衡级为第N级，塔顶冷凝器为N+1级，自塔顶冷凝器（第N+1级）至精馏段任一级作i组分的物料平衡， vn,i=ln+1,i+di 等号两侧乘以Ani/di An,ivn,i/di

8、=An,i(ln+1,i/di+1) 因ln,i=An,ivn,i, 故 ln,i/di=An,i(ln+1,i/di+1),(6-1),(6-2),上式即为精馏段逐级计算公式，由任一级的比流率值（ln+1,i/di）可计算出其下一级的比流率值（ln,i/di），如果已知各级的吸收因子且能确定塔顶冷凝器的比流率值（lN+1,i/di），可从塔顶开始逐级向下利用逐级计算公式计算出精馏段各级的比流率值。塔顶冷凝器的比流率可丛回流与产品的关系进行计算。,塔顶冷凝器为全凝器时， lN+1,i/di= LN+1xN+1,i/(Dxd,i) =LN+1/D=R 塔顶冷凝器为分凝器时， lN+1,i/di=

9、 LN+1,ixN+1,i/(Dyd,i) = R/kN+1,i 计算出lN+1,i/di后，就可利用逐级计算公式计算lNi/di，进而计算精馏段各级的比流率值，最后计算出lM+1,i/di。,(6-3),(6-4),自塔底再沸器（第一级）至提馏段任 一级（如第n级）作i组分的物料平衡 ln+1,i=vn,i+wi 等号两侧除以wi ln+1,i/wi=vn,i/wi+1 因vn,i=Sn,iln,i,故 ln+1,i/wi=Sn,i(ln,i/wi)+1 上式为提馏段逐级计算公式,(6-3),(6-4),由任一级的比流率值(ln,i/wi)可计算其上一级的比流率值(ln+1,i/wi)，如果

10、已知各级的蒸出因子且能确定塔底再沸器的比流率值（l1,i/wi），可从塔底开始逐级向上利用逐级计算公式计算出提馏段各级的比流率值。,塔底再沸器的出料流率等于离开第一级的液相流率，故 l1,i/wi=1 将结果代入计算公式可计算出l2,i/wi l2,i/wi=S1,i(l1,i/wi)+1=S1,i+1 反复运用逐级计算公式，可计算出提馏段各级的比流率值，最后计算出 lM+1,i/wi。,用逐级计算法分别从塔顶和塔底计算至进料级，所计算出的比流率值应满足全塔的物料平衡，如不满足说明原假设条件有误，对假设条件进行调整重新进行计算，直到满足要求，即所谓进料契合。判别是否满足契合要求的判据定为塔顶产

11、品流率的计算值与规定值之差是否小于规定的允许误差。不满足要求时，将计算结果直接用作下一轮迭代的假设值，故称为直接迭代法。,首先，通过物料平衡关系关联塔底产品和塔顶产品各组分的流率比。 精馏段组分物料平衡 VM,i=lM+1,i+di 等式两边除以di, VM,i/di=lM+1,i/di+1 提馏段组分物料平衡 VM,i=lM+1,i-wi 等式两边除以wi, VM,i/wi=lM+1,i/wi-1,(6-5),(6-6),全塔组分组分物料平衡 fi=wi+di 等式两边除以di, fi/di=wi/di+1 作进料级物料平衡，关联精馏段和提馏段组分物料关系 VMi=VM,i+fVi 等式两边

12、除以di, VMi/di=VMi/di+fVi/di 上式可转化为：,(6-9),(6-8),(6-7),将式（6-7），（6-8）和（6-9）代入上式并整理得： 当进料为汽液两相时，可由精馏段及提馏段逐级计算所得lM+1,i/di及lM+1,i/wi,用上式计算出塔底塔顶产品组分流率比wi/di，,(6-10),若进料为饱和液体或过冷液体时，则fV i=0, 式（6-10）化为： 若进料为饱和蒸汽或过热蒸汽时，则fV,i=fi，式（6-10）化为,(6-11),(6-12),用wi/di值通过全塔组分物料平衡，计算塔顶产品率: fi=di+wi=di (1+wi/di) di=fi/ (1+

13、wi/di) wi=(wi/di)di Dcal为塔顶产品流率计算值,(6-13),如果Dcal与原规定的产品流率值D接近相等，满足精度要求时，说明原设的温度与流率值合适，否则说明原设不当，需重新设定进行计算，直至满足要求，获得正确结果。,用直接迭代法求得各级各组分液相组成，校正各级温度。 新一轮迭代计算所用的各级温度用泡点法求算，首先直接用逐级计算所得比流率值求出各级液相组成，据此求算各级泡点即新一轮迭代温度。 精馏段各级各组分液相组成为：,(6-14),提馏段各级各组分液相组成为： 有了各级液相组成，可计算在操作压力下的泡点温度。如果需要对流率进行校正，要靠热量平衡计算。,(6-15),操

14、作型逐级计算步骤总结如下： 1）输入进料条件F, zi, q, 操作压力P, 塔顶产品流 率D, 回流比R，平衡级数N及进料位置M,求 出各级汽液相流率V, L, L； 2）假定各级温度Tn，求出各级各组分平衡常数 kj,i； 3）求精馏段吸收因子Ani及提馏段蒸出因子Sni； 4）精馏段逐级计算，计算各级各组分的比流率 值lni/di，计算至进料级的上一级lM+1,i/di；,5）提馏段逐级计算，计算各级各组分的比流率 值lni/wi，计算至进料级的上一级lM+1,i/wi； 6）进料级契合，求出塔底塔顶组分的比流率值 wi/di, 塔顶产品流率di，塔底产品流率wi； 7）判断塔顶产品流率

15、计算值Dcal=与规定值D之 差是否 满足精度要求ABS(D-Dcal)?,如果 满足则计算结束，输出结果； 8）否则，计算各级各组分液相摩尔分率xni，求 各级泡点温度Tn及各级各组分的平衡常数 kni； 9）返回第(3)步。,习题： 某精馏塔有9个平衡级（包括塔釜），进料在第5级（自下数），处理含A，B，C，D四组分的饱和溶液，进料量100摩尔/小时，塔顶设部分冷凝器，回流比为4，B、C组分的回收率为98%，按恒摩尔流假定确定流率，操作条件下各组分的平衡常数按下式计算。 原料组成及各组分相对挥发度,Ki=ai+bit (t: ) 常数ai, bi列表如下： 试用操作型算法逐级算出各级温度及

16、组成。,2)、 收敛法,直接迭代算法直接利用组成的计算值计算泡点温度，使计算在远离收敛解处盘桓，导致收敛速度慢，甚至产生发散或震荡而得不到结果，如果能对组成的计算值作合理修正，使其接近收敛要求，会大大改善收敛性质。收敛法是一种加速收敛的方法，其计算步骤于直接迭代法基本相同，只是增加一个搜寻修正系数 的步骤，用来修正组成的计算值。,收敛法和直接迭代法一样，先分别由塔顶和塔底计算比流率值，lni/di及lni/wi，通过进料级契合计算出塔底塔顶产品组分的比流率值wi/di，此时不是用di=fi/ (1+wi/di)、wi=(wi/di)di直接求出塔顶、塔底产品组分流率di及wi，而是通过下面的方

17、法求出校正后的di及wi的值。,具体做法是： 对全塔做组分的物料衡算： fi=di+wi 由此可得： 显然，对于正确的di，应满足：,(6-16),(6-17),(6-18),(6-19),以上式子中的 值，可有前面学过的(6-10)式求出： 为了得到正确的di，用 乘上一个系数之后，再计算di:,(6-10),(6-20),所以，wi为： 下标co表示校正后的值，下标ca表示计算值。 将式(6-20)带入式(6-19)中，得到：,(6-21),(6-22),上式可写成如下函数形式 : 上式中，除变量外均为已知量，是一个关于的一元方程，可通过牛顿迭代法求解，g()函数的图形如图所示。,(6-23),图6-1 g()函数,由图可见，取值范围在0之间 在0区域里， g()函数是单调下降的，取=0为初值进行牛顿迭代，必定收敛于g()=0的根= r处， g()之导数式为：,(6-24),(6-25),(6-26),由第r次迭代值求第r+1次迭代值时，按下列公式进行计算 通常收敛判据用相对偏差形式： 在未满足 之前，即未达=1之前，(di)co及(wi)co仍有近似值，作为下一轮迭代计算之用。,(6-27),收敛法用下列公式修正其液相组成： 精馏段液相组成(摩尔