专题综述
文章编号:1009 9425(2003) 02 0001 05
精馏计算中相对挥发度 的选择
徐 文 灏
(北京西城区月坛北街铁三区5栋2门3号 100045)
摘要:推演了在精馏计算中相对挥发度 的选择对最少理论塔板数N m 及最小回流比R m 的决定影响。从基础理论上阐明了 的选用原则。并且对几个经典理论公式的相互关系作出了解释。
关键词:精馏; 相对挥发度; 制氩; 馏分; 回流比; 理论塔板数
中图分类号:TQ028 13 文献标识码:A
Choosing relative volatility in calculation of distillation
Xu Wen hao
(Designing and Research Institute o f Beijing Capital Iron and Steel Com pany (pre retirement ) , No . 3, Unit 2, Bldg . 5, Tie San Area , North Moon Temple Street , West U rban District , Beijing 100045)
Abstract :The decisive influence of relative volatility chosen in calculation of distillation on the minimum number of theoretical plates N m and the minimum reflux ratio R m is given by derivation and deduction. The principle of choosing is expounded on the basis of basic theory. In addition the c orrelation of several classical theoretical equations is ex plained with c ogenc y and relevance.
Keywords :Distillation; Relative volatility; Argon production; Fraction; Reflux ratio; Theoretical plate
笔者在论文 全精馏制氩中回流比R 的抉择与讨论
[1]
择, 亦都不加说明。这与未能深入研究有关系。笔者在空分精馏计算的多篇论文中阐述了许多有意思的结果及问题, 将借助以往的计算, 在基础理论上
作出阐述, 使之充实精馏基础理论。
中得出结论:用入料馏分层的 值和塔顶馏
出产品的 值的平均值( lh ) av 计算的理论塔板数是精馏任务最低要求的值; 而仅用入料馏分层的( lh ) F 值求得的理论塔板数却是需要的最高极限值。此概念可推至其他精馏, 这为相对挥发度 的选择提出了一个明确的原则
[1]
1 最少理论塔板数(芬斯克公式1932年)
公式推演与探讨
以l 代表轻关键组分, h 代表重关键组分。
y l x l (P l /P ) x l P l x l
=== lh y h x h x h (P h /P ) x h P h y l x l 即(y ) 1=( lh ) 1(x ) w
h
h
[2]
。
早年, 笔者曾浏览过一些经典公式的原著论文, 但因没有带着问题, 未曾注意到公式对 值的要求, 当然原著对各自公式之间的关系也很少涉及。以后, 在工作中亦都延用了塔上下的( lh ) a v 平均值, 因为各种书籍和手册亦都是不加阐述, 更无讨论。从未见有专论深入探讨此问题。虽然亦偶有不同的选
收稿日期:2003 01 03
(1)
由物料平衡知:
V 1 y l (1)=L 2 x l (2) -W x l (w )
作者简介:徐文灏(1933 ) , 男, 教授级高工。1955年毕业于浙江大学化工系, 北京首钢设计院工作, 1958年以来, 先后主持过国产3350、日本6000、法国6500、德国30000m 3/h 空分装置的工程设计、引进建设投产工作。著有40多篇论文, 获多次奖励。
本文是在作者1999年, 2001年一系列低温全精馏制氩工艺计算演绎论文上对精馏基础理论的一个研究成果。
深冷技术 2003年第2期 Cryogenic Technology 2 2003
1
专题综述
在回流比为 , 即全回流时, V 1=L 2, W =0, 则y l (1) =x l (2) , 同理y h (1) =x h (2) , 故(y l /y h ) 1=(x l /x h ) 2, 代入(1) 式, 则:
(x l /x h ) 2=( lh ) 1(x l /x h ) w
徐文灏:精馏计算中相对挥发度 的选择
(2) N m 值除与塔顶塔底的馏分组成要求有关外, 直接受制于( lh ) a v 。就是说不管是哪一组多元系统精馏, 只要入料及馏出的原料及产品组成相同, N m 就只决定于( lh ) a v 。当然, ( lh ) av 亦是由该多元系统的入料及馏出组分所决定。
(3) 由此延伸, 不管什么计算精馏塔板数的方法, 对同一任务其计算值只有准确度的区别, 但应该
是相近的。因为计算方法的理论都是正确的。这就是说, 凡是计算中用了( lh ) a v 值, 所求得的N 值都应该是精馏所要求的最少理论塔板数, 而且是相近的N m 值。
(4) 正如前述, 上述N m 公式是全回流下的 静态 精馏, 所得的N m 值与用各种方法计算的N 值是不一样的。各种方法是利用( lh ) a v 求得的, 实际生产情况下, 只有馏出产品的真正的 动态 精馏, 其N 值要比N m 大许多, 甚至达到一倍半。但因为用的是( lh ) av , 所以求得的塔板数比起用大于(
lh ) a v 值计算的塔板数, 还应该是最少, 如图1所示, 是实际生产计算中的最少理论塔板要求值。
(2)
式中:P 、P L 、P h 代表平衡液面的总压和分压; V 代表上升蒸汽量; L 代表回流液体量; W 代表釜底放出液体量(以W 、1、2 分别代表釜底、第一层、第二层塔板 ) 。
同样, 由平衡关系可得:(y l /y h ) 2=( lh ) 2 (x l /x h ) 2
由物料平衡可得:(y l /y h ) 2=(x l /x h ) 3故: (x l /x h ) 3=( lh ) 2(x l /x h ) 2=
( lh ) 2( lh ) 1(x l /x h ) w
依此类推, 直至塔顶, 可得:(x l /x h ) D =( lh ) n ( lh ) n -1( lh ) n -2
( lh ) 2( lh ) 1(x l /x h ) w
如果各理论板的 lh 均相等, 则( lh ) n ( lh ) 1=( lh )
n
(3)
显然, 因为是全回流, 所得的n 值一定是最少, 故n 即以N m 表示之。但即使是全回流下, 各层的 lh 亦是不等的, 但可假定为:
( lh ) n ( lh ) 1=( lh ) a v
式中:( lh ) a v =
Nm
Nm
( lh ) n ( lh ) 1
由于各层中的( lh ) i 值都是未知的, 只有入料层(F ) 、塔顶(D ) 、釜底(W ) 的 lh 是可知的(其液相或汽相组成可以由物料平衡求出) 。
因此,
2
可以
3
( lh ) D ( lh ) F ( lh ) w 或
图1 纯氩塔精馏N-R 曲线图
( lh ) D ( lh ) F 代表( lh ) a v , 这样(3) 式即变为:N m =
=
lg[(x l /x h ) D /(x l /x h ) w ]
lg( lh ) a v lg[(x l /x h ) D (x h /x l ) w ]lg( lh ) av
(4)
图1曲线说明, 就是再加大R, 为了满足塔顶D 的产质产量, N 值都不会有大的降低。这样求得的N 值就是借助( lh ) a v 所求得的生产所需最低塔板值。这一概念即成为讨论R m 中 选值的讨论基础。
由于是逐层计算而得到的全回流下的理论塔板数, 所以N m 应是全回流下精馏需要的最少理论塔板数。任何有D 排出的情况都将带来回流量的减少, 回流比由 开始降低, 引起塔板数的增多。
结合以往计算的实例, 将最少理论塔板数的概念扩展阐述如下:
(1) 毫无疑问, 计算N m 必须采用( lh ) av 值。
2 最小回流比(柯尔本公式1941年)
式推演及探讨
精馏段一般操作方程式
[3]
[3, 4]
公
为:
(5)
y =xR /(R +1) +x D /(R +1) 式中 x D 塔顶产品组成
R 回流比
2 2003年第期
徐文灏:精馏计算中相对挥发度 的选择
专题综述
R m 值。而精馏段的R m 就只能用入料层的 lh (F ) 值和入料层及塔顶产品组成求得, 从而所求得的N 值才是真正的最大值。果然, 在文献[1]中计算的185层以上各层的y 值基本没有变化, 185层已是可读的最大值, 实际就相当于 了。如果采用小于这算得的R m 值, 将会得到N = 的结果, 见图1。
(3) 化学工程手册
(6)
[6]
当回流比为最小时, 塔板数即为 , 此时精馏段某一层的两关键组分的操作线方程式为:
R m x l (x l ) D
y l =(R +1) +(R +1) ,
m m R m x h (x h ) D
y h =(R +1) +(R +1)
m m 故:
y l R m x l +(x l ) D x l ==( lh ) y h R m x h +(x h ) D x h
将公式(6) 中 lh 明确指
(x l ) D 1lh (x h ) D
即R m =[-]
( lh -1) x l x h
定为进料层的 lh (F ) , 这是正确的。但未作讨论和阐明, 更未提及其他文献采用( lh ) a v 的可否与比较。必须指出, 若是真正的最小R m 值, 它求得的塔板数N 即使不是最高的, 亦是偏高的。因为入料层以上的各
层的 只有将实际采用的R =lh 值都会比 lh (F ) 小。AR m 中的A 取较大值, 才能求出合适的考虑到建设费及运行费的N 值。而用( lh ) a v 所求出的R m 值及取较小A 值的R 计算的N 值, 可以视作完成精馏任务的底限值。如图1中所示氩精馏中的150层即是一种明确的答案, 显然计算的N 值是最低实际需要
[5]
值。在氩塔计算中所引美、德专利中都说 最少需要150层 , 最好需要180层 或 180~200层 。其实, 180层以上到200层精馏效果几乎是没有进展的; 150到180层, 对提高纯氩质量亦无多少改进, 因
-6
为都已经达到氧含量低于10, 增加塔板数只会徒然加大空气压力。
(4) 关于恩特乌德法
[1、2]
结合以往实例, 将最小回流比的概念扩展阐解如下:
(1) 在许多书籍及论文中, 对 lh 都未加以讨论, 即延用了( lh ) a v 的求法:是最近出版的新书
[4]
特别( lh ) D ( lh ) F 值。
明确规定用此方法, 条件是只
要入料和塔顶的 lh 值相差不多。但是, 由本文1 (3) 中阐述, 只要用( lh ) a v 值, 不管用什么计算方法所求得的理论塔板数都应该是生产所需的最低限度塔板数。以氩塔计算实例
[5]
说明, 用上述R m 公式, 采用
( lh ) a v 值, 逐层塔板计算法算得的塔板数大抵都在150层左右, 如图1所示。这使人恍然理解, 由于( lh ) a v 值是偏小的, 造成R m 偏大, 再加上采用的R =AR m (A >1) , 促成N 值计算就接近不变的最小值。
(2) 从公式(6) 来看, 显然 lh 是指某一层的 lh 。如果各层的 lh 非常接近, 自然用( lh ) av 计算, 偏差尚不大, 但往往入料和塔顶的 lh 值差别不可忽视, 在精细的塔板计算中将会有较大误差。
而且由公式(6) 看, 对应 lh 的是这一层的组分汽液相值, 采用( lh ) a v 与这一层组分不对应, 亦必然会导致误差, 在理论上亦说不通。按理, 最小回流比是指采用小于此回流比值, 塔板数就将无穷大, 这不符合有限塔板的精馏情况。有读者已经发现在1999年第2期 深冷技术
[5]
计算R m 公式(1945)
(7a) (7b)
i (x D , i ) m
-Q =R m +1
i
i x F , i )
=1-q
i -Q
中最难挥发的组分为基准)
式中 i 组分i 的相对挥发度(以进料混合物
q 进料的液相分率
x F , i 进料混合物中组分i 的摩尔分率(x D , i ) m 最小回流比下馏出液中组分i
的摩尔分率
Q 方程式的根(对于有C 个组分的系统
将有C 个根, 只取 l
一个根)
计算比较麻烦, 首先要按芬斯克公式由关键组分的分配比及全回流下估算出馏出液的组成, 然后用试差法由式7a 求Q 代入式7b 求R m , 例题见文献[1]中。
3中的表5中只算到氩馏分含
氩8%, 而不象表6中算到7%Ar?这是因为算到7%时, R m 就大于所采用的R =55, 即55小于R m , 应该是理论板数无穷大, 但计算的N 值仍然是可读值, 略高于150, 当时不理解亦无法解释。为了不影响课题主攻方向, 只好把它悬置而略去。现在从理论上详细分析就明白了:用( lh ) a v 所求得的R m 根本不能是这一层的最小值, 更不是入料层以上精馏塔的最小深冷技术 2003年第2期
Cryogenic Technology 2 2003
专题综述
有一点亦是明确的, 所用 i 都是进料组分的 值, 因为只有进料层的x i 、 这与上述认i 是可知的。定R m 应按入料层 反过来说, 只要lh 计算是一致的。( lh ) a v 与( lh ) F 比较接近, 用( lh ) a v 代入计算所得的R m 亦会偏大, 而进一步求得的N 值亦将是低限值。
徐文灏:精馏计算中相对挥发度 的选择
3 简捷法精馏计算 线
[2、4]
吉利兰曲
(1940) 分析与探讨
[4]
吉利兰根据61座两组分和多组分精馏塔的逐层计算结果, 整理得出关联曲线, 利用此曲线可以
用N m 、R m 、R 求出理论塔板数N , 此即简捷法。如图2
所示。
图2 吉利兰曲线
图2中曲线左端延线表示在最小回流比下的操作情况, 此时(R -R m ) /(R +1) 接近零, 而(N -N m ) /(N +1) 接近1(即N = ) 。而曲线右端表示在全回流下操作状况, 此时(R -R m ) /(R +1) 接近1(即R = ) , (N -N m ) /(N +1) 接近零, 即N N m
[7]
所归结的:只要用( lh ) a v 来作计算, 不管什么计算方法, 求得的理论塔板数总是低限值, 而低值总是接近的, 正如图1所示。结果一致的现象由此得到解释。反过来说, 计算结果一致正好证明了用( lh ) a v 去计算塔板数所得的都是完成精馏的低限值, 见表1。
表1 法国6500m 3/h 空分下塔精馏计算结果汇总表
简捷法
下塔段位 - 段 - 段
逐层法812 48
(用各段的R m 计算) 7 5813 08
简捷法(以全塔的R m 计算) 7 88
12 62
实际塔板数(由此可估计板效率)
1414
[3]
。
结合以往计算实例, 可以得出如下一些概念:(1) 既然曲线中用了N m , 则其所用的一定是
( lh ) a v 。在同一计算中, 既然N m 用了( lh ) av , 则可以判定R m 中亦一致地用了( 因此所求得的N 值lh ) av 。必然是完成精馏任务的低限值, 这是以上讨论顺理成章的理解。
(2) 或者正是因为上述吉利兰曲线计算N m 、R m
时都采用( lh ) a v , 就造成了习惯性的计算通用, 再加上前面提到的通常的如果塔顶塔底的 lh 值相差不大时可以用( , 就使概念混淆不lh ) a v 的明确解释清, 掩盖了应该用( 又lh ) F 来计算R m 的正确概念。因为用( lh ) a v 求得的R m 及N 亦是能完成任务的低限值, 不会出偏差。所以就这样长期使用下去, 亦不为过。
(3) 在以往的空分计算中, 用逐层法和简捷法
[3]
计算的结果非常接近。当时虽然很满意, 但亦非常惊奇, 无法解释其一致性。现在看来, 正如以上讨论
[4]
表2 法国6500m 3/h 空分上塔精馏段计算结果汇总表
上塔段位 段(由上向下)
段 段
段 ~ 段全塔
逐层法简捷法87471826
7 967 694 246 3718 326 26
实际塔板数
15
128446479
塔板效率/%
53645014 5
(因有膨胀过热气影响)
28 633 2
(4) 十分明显, 在吉利兰曲线中亦包含了R 与R m 之间有一个比值系数A 。由于用( lh ) av 所求得的R m 已足够大, 如文献[1、5]所计算的那样, 因此此系数值不会很大, 大亦无益。据文献[8]中所写: 根据
2003年第期
2
徐文灏:精馏计算中相对挥发度 的选择
专题综述
间取几个R m 值计算出各自的N 值, 绘成如图1的曲线, 就可以很科学地决定选用的R 及N 值来。因为高限N 值与低限N 值有的相差还是较大的, 不能轻信塔顶与顶底 lh 值相差不大时可采用( lh ) a v 的说法。
(4) 对于精密精馏, 只能用逐层法计算; 简捷法只适用于分离度要求不高, 满足于求得需要的塔板最低限度值的任务。但简捷法由于使用了R m 、N m 等经典公式, 却能带来精馏概念的廓清。
(5) 所有这些讨论都是泛指所有的精馏工况。因此这些结论都适用于各精馏物系。计算实例很好地佐证了这些结论。
自1978年以来, 笔者对空分及氩精馏作了工艺理论计算的系统探讨, 发表了系列论文, 形成完整的计算方法与实例, 其中亦廓清了许多概念。其中N 2 Ar O 2汽液平衡关联式及本文还对基础理论作了发展, 它们都已验证于实验数据与计算实用中。近年来同业还有同志作了另外的开拓, 如在操作线物料平衡中引入了热量平衡, 即进行不等分子回流的计算。精馏是空分装置的核心与主体, 美、德、法、英、日以及俄罗斯的大公司都有自己的计算程序和数据库, 各公司几十年来在精馏上所取得的进展亦是有目共睹, 突飞猛进, 塔板数一加再加, 从而促进和带动了塔板结构、设置的变化乃至稀有气体提取的进展。深切希望同业抓住这个中心继续消化吸收, 独立钻研, 赶上先进并创造出我们自己的理论与实践成果, 带动行业的发展。
参考文献:
[1]徐文灏. 全精馏制氩中回流比R 的抉择与讨论. 深冷技术,
2001(4):1~12.
[2]陈洪舫编. 基本有机化工分离工程. 化工版, 1981:120~121,
125.
[3]徐文灏. 全低压空分装置相对挥发度法精馏计算(上) 逐层
法.
(下) 简捷法. 深冷技术, 1979(2) :1~24; 1979(3) :12~28.
[4]邓修, 吴俊生编著. 化工分离工程. 科学版, 2002:55~56. [5]徐文灏. 低温全精馏制氩工艺计算演绎之一、二. 深冷技术,
1999(1)1~12; 1999(2) 1~9; 1999(3) :16~18. [6]化学工程手册:第三卷. 化工版. P -74、75.
[7]贾绍仪. 蔡诚敬主编. 化工传质与分离工程. 化工版. 2001年
1月第1版:179.
[8]刘家琪主编. 分离工程. 化工版. 2002年5月第一版:107.
Pair J. R 和Bolles W. l: Chem. Eng. 1969(75) :56 的研究结果, R /R m 的最优化值约为1 05, 但是在比值稍大的一定范围内都接近最佳条件。 对于每一物系都可按此考虑, 吉利兰曲线显示在(R -R m ) /(R +1) 中包含了这一因素。从氩精馏来看, 1 05亦是合适的。这与文献[5]中美国专利4994098所采用的1 04~1 06是一致的。
(5) 吉利兰线图绘制的依据是八种物系在广泛的精馏条件下由逐板计算法得出的结果。组分数目2~11个, 进料热状态包括冷料至过热蒸汽等五种情况, R m 为0 53~7 0, 组分间相对挥发度为1 26~4 05, 理论板数为2 4~43 1。可见吉利兰线图适用于石油、苯、酚系等蒸馏要求分离度不高的物系。因为它所求得的亦正是逐层法所求得的最低限度需要塔板数, 此值在R 一定大时, 是趋于一致的。在空分上能成功计算即是因此。用简捷法还解决了过去粗氩塔计算中塔顶95%Ar 含量的粗氩中氮、氧含量随进料组成而变的假定值, 这亦是简捷法的贡献。但对于精密精馏的纯氩计算则只能用逐层法, 唯逐层法才能分别用( lh ) F 和( lh ) av 判别任务需要的最高和最低塔板极限
[1]
[5][5]
[7]
。
4 结语
通过这样以实际计算结果佐证的对理论的探讨, 是否可以归结为以下结论:
(1) 由于全回流(R ) 下求N m 必须采用( lh ) a v 。因此可以推论:用( lh ) a v 在逐层法或简捷法中计算的R m 值已足够大, 所计算的N 值必然是精馏所需的最低限度塔板数。但此值不是全回流下的N m 值, 而是实际生产所需最低值。
(2) 用入料层的( lh ) F 求得的R m 才是该精馏任务的最小回流比, 小于此R m , N , 精馏不能进行。采用略大于R m 的R 计算出来的塔板数N 才是实际需要的最多塔板数, 实际上此值是偏高的, 因为进料层以上塔板的( lh ) i 都比( lh ) F 小, 这导致计算中最上的多层塔板纯度虽稍有提高, 但徙增塔板阻力。
(3) 对于纯产品的精密精馏, 应该用( lh ) F 计算出R m 、R (R 可采用1 05R m 左右) , 再用逐层法计算出高限N 值来; 再用( R , 计算出低限lh ) a v 求R m 、N 值来。借助计算机程序计算, 可以在R m 与R m 之深冷技术 2003年第2期
Cryogenic Technology 2 2003
5
专题综述
文章编号:1009 9425(2003) 02 0001 05
精馏计算中相对挥发度 的选择
徐 文 灏
(北京西城区月坛北街铁三区5栋2门3号 100045)
摘要:推演了在精馏计算中相对挥发度 的选择对最少理论塔板数N m 及最小回流比R m 的决定影响。从基础理论上阐明了 的选用原则。并且对几个经典理论公式的相互关系作出了解释。
关键词:精馏; 相对挥发度; 制氩; 馏分; 回流比; 理论塔板数
中图分类号:TQ028 13 文献标识码:A
Choosing relative volatility in calculation of distillation
Xu Wen hao
(Designing and Research Institute o f Beijing Capital Iron and Steel Com pany (pre retirement ) , No . 3, Unit 2, Bldg . 5, Tie San Area , North Moon Temple Street , West U rban District , Beijing 100045)
Abstract :The decisive influence of relative volatility chosen in calculation of distillation on the minimum number of theoretical plates N m and the minimum reflux ratio R m is given by derivation and deduction. The principle of choosing is expounded on the basis of basic theory. In addition the c orrelation of several classical theoretical equations is ex plained with c ogenc y and relevance.
Keywords :Distillation; Relative volatility; Argon production; Fraction; Reflux ratio; Theoretical plate
笔者在论文 全精馏制氩中回流比R 的抉择与讨论
[1]
择, 亦都不加说明。这与未能深入研究有关系。笔者在空分精馏计算的多篇论文中阐述了许多有意思的结果及问题, 将借助以往的计算, 在基础理论上
作出阐述, 使之充实精馏基础理论。
中得出结论:用入料馏分层的 值和塔顶馏
出产品的 值的平均值( lh ) av 计算的理论塔板数是精馏任务最低要求的值; 而仅用入料馏分层的( lh ) F 值求得的理论塔板数却是需要的最高极限值。此概念可推至其他精馏, 这为相对挥发度 的选择提出了一个明确的原则
[1]
1 最少理论塔板数(芬斯克公式1932年)
公式推演与探讨
以l 代表轻关键组分, h 代表重关键组分。
y l x l (P l /P ) x l P l x l
=== lh y h x h x h (P h /P ) x h P h y l x l 即(y ) 1=( lh ) 1(x ) w
h
h
[2]
。
早年, 笔者曾浏览过一些经典公式的原著论文, 但因没有带着问题, 未曾注意到公式对 值的要求, 当然原著对各自公式之间的关系也很少涉及。以后, 在工作中亦都延用了塔上下的( lh ) a v 平均值, 因为各种书籍和手册亦都是不加阐述, 更无讨论。从未见有专论深入探讨此问题。虽然亦偶有不同的选
收稿日期:2003 01 03
(1)
由物料平衡知:
V 1 y l (1)=L 2 x l (2) -W x l (w )
作者简介:徐文灏(1933 ) , 男, 教授级高工。1955年毕业于浙江大学化工系, 北京首钢设计院工作, 1958年以来, 先后主持过国产3350、日本6000、法国6500、德国30000m 3/h 空分装置的工程设计、引进建设投产工作。著有40多篇论文, 获多次奖励。
本文是在作者1999年, 2001年一系列低温全精馏制氩工艺计算演绎论文上对精馏基础理论的一个研究成果。
深冷技术 2003年第2期 Cryogenic Technology 2 2003
1
专题综述
在回流比为 , 即全回流时, V 1=L 2, W =0, 则y l (1) =x l (2) , 同理y h (1) =x h (2) , 故(y l /y h ) 1=(x l /x h ) 2, 代入(1) 式, 则:
(x l /x h ) 2=( lh ) 1(x l /x h ) w
徐文灏:精馏计算中相对挥发度 的选择
(2) N m 值除与塔顶塔底的馏分组成要求有关外, 直接受制于( lh ) a v 。就是说不管是哪一组多元系统精馏, 只要入料及馏出的原料及产品组成相同, N m 就只决定于( lh ) a v 。当然, ( lh ) av 亦是由该多元系统的入料及馏出组分所决定。
(3) 由此延伸, 不管什么计算精馏塔板数的方法, 对同一任务其计算值只有准确度的区别, 但应该
是相近的。因为计算方法的理论都是正确的。这就是说, 凡是计算中用了( lh ) a v 值, 所求得的N 值都应该是精馏所要求的最少理论塔板数, 而且是相近的N m 值。
(4) 正如前述, 上述N m 公式是全回流下的 静态 精馏, 所得的N m 值与用各种方法计算的N 值是不一样的。各种方法是利用( lh ) a v 求得的, 实际生产情况下, 只有馏出产品的真正的 动态 精馏, 其N 值要比N m 大许多, 甚至达到一倍半。但因为用的是( lh ) av , 所以求得的塔板数比起用大于(
lh ) a v 值计算的塔板数, 还应该是最少, 如图1所示, 是实际生产计算中的最少理论塔板要求值。
(2)
式中:P 、P L 、P h 代表平衡液面的总压和分压; V 代表上升蒸汽量; L 代表回流液体量; W 代表釜底放出液体量(以W 、1、2 分别代表釜底、第一层、第二层塔板 ) 。
同样, 由平衡关系可得:(y l /y h ) 2=( lh ) 2 (x l /x h ) 2
由物料平衡可得:(y l /y h ) 2=(x l /x h ) 3故: (x l /x h ) 3=( lh ) 2(x l /x h ) 2=
( lh ) 2( lh ) 1(x l /x h ) w
依此类推, 直至塔顶, 可得:(x l /x h ) D =( lh ) n ( lh ) n -1( lh ) n -2
( lh ) 2( lh ) 1(x l /x h ) w
如果各理论板的 lh 均相等, 则( lh ) n ( lh ) 1=( lh )
n
(3)
显然, 因为是全回流, 所得的n 值一定是最少, 故n 即以N m 表示之。但即使是全回流下, 各层的 lh 亦是不等的, 但可假定为:
( lh ) n ( lh ) 1=( lh ) a v
式中:( lh ) a v =
Nm
Nm
( lh ) n ( lh ) 1
由于各层中的( lh ) i 值都是未知的, 只有入料层(F ) 、塔顶(D ) 、釜底(W ) 的 lh 是可知的(其液相或汽相组成可以由物料平衡求出) 。
因此,
2
可以
3
( lh ) D ( lh ) F ( lh ) w 或
图1 纯氩塔精馏N-R 曲线图
( lh ) D ( lh ) F 代表( lh ) a v , 这样(3) 式即变为:N m =
=
lg[(x l /x h ) D /(x l /x h ) w ]
lg( lh ) a v lg[(x l /x h ) D (x h /x l ) w ]lg( lh ) av
(4)
图1曲线说明, 就是再加大R, 为了满足塔顶D 的产质产量, N 值都不会有大的降低。这样求得的N 值就是借助( lh ) a v 所求得的生产所需最低塔板值。这一概念即成为讨论R m 中 选值的讨论基础。
由于是逐层计算而得到的全回流下的理论塔板数, 所以N m 应是全回流下精馏需要的最少理论塔板数。任何有D 排出的情况都将带来回流量的减少, 回流比由 开始降低, 引起塔板数的增多。
结合以往计算的实例, 将最少理论塔板数的概念扩展阐述如下:
(1) 毫无疑问, 计算N m 必须采用( lh ) av 值。
2 最小回流比(柯尔本公式1941年)
式推演及探讨
精馏段一般操作方程式
[3]
[3, 4]
公
为:
(5)
y =xR /(R +1) +x D /(R +1) 式中 x D 塔顶产品组成
R 回流比
2 2003年第期
徐文灏:精馏计算中相对挥发度 的选择
专题综述
R m 值。而精馏段的R m 就只能用入料层的 lh (F ) 值和入料层及塔顶产品组成求得, 从而所求得的N 值才是真正的最大值。果然, 在文献[1]中计算的185层以上各层的y 值基本没有变化, 185层已是可读的最大值, 实际就相当于 了。如果采用小于这算得的R m 值, 将会得到N = 的结果, 见图1。
(3) 化学工程手册
(6)
[6]
当回流比为最小时, 塔板数即为 , 此时精馏段某一层的两关键组分的操作线方程式为:
R m x l (x l ) D
y l =(R +1) +(R +1) ,
m m R m x h (x h ) D
y h =(R +1) +(R +1)
m m 故:
y l R m x l +(x l ) D x l ==( lh ) y h R m x h +(x h ) D x h
将公式(6) 中 lh 明确指
(x l ) D 1lh (x h ) D
即R m =[-]
( lh -1) x l x h
定为进料层的 lh (F ) , 这是正确的。但未作讨论和阐明, 更未提及其他文献采用( lh ) a v 的可否与比较。必须指出, 若是真正的最小R m 值, 它求得的塔板数N 即使不是最高的, 亦是偏高的。因为入料层以上的各
层的 只有将实际采用的R =lh 值都会比 lh (F ) 小。AR m 中的A 取较大值, 才能求出合适的考虑到建设费及运行费的N 值。而用( lh ) a v 所求出的R m 值及取较小A 值的R 计算的N 值, 可以视作完成精馏任务的底限值。如图1中所示氩精馏中的150层即是一种明确的答案, 显然计算的N 值是最低实际需要
[5]
值。在氩塔计算中所引美、德专利中都说 最少需要150层 , 最好需要180层 或 180~200层 。其实, 180层以上到200层精馏效果几乎是没有进展的; 150到180层, 对提高纯氩质量亦无多少改进, 因
-6
为都已经达到氧含量低于10, 增加塔板数只会徒然加大空气压力。
(4) 关于恩特乌德法
[1、2]
结合以往实例, 将最小回流比的概念扩展阐解如下:
(1) 在许多书籍及论文中, 对 lh 都未加以讨论, 即延用了( lh ) a v 的求法:是最近出版的新书
[4]
特别( lh ) D ( lh ) F 值。
明确规定用此方法, 条件是只
要入料和塔顶的 lh 值相差不多。但是, 由本文1 (3) 中阐述, 只要用( lh ) a v 值, 不管用什么计算方法所求得的理论塔板数都应该是生产所需的最低限度塔板数。以氩塔计算实例
[5]
说明, 用上述R m 公式, 采用
( lh ) a v 值, 逐层塔板计算法算得的塔板数大抵都在150层左右, 如图1所示。这使人恍然理解, 由于( lh ) a v 值是偏小的, 造成R m 偏大, 再加上采用的R =AR m (A >1) , 促成N 值计算就接近不变的最小值。
(2) 从公式(6) 来看, 显然 lh 是指某一层的 lh 。如果各层的 lh 非常接近, 自然用( lh ) av 计算, 偏差尚不大, 但往往入料和塔顶的 lh 值差别不可忽视, 在精细的塔板计算中将会有较大误差。
而且由公式(6) 看, 对应 lh 的是这一层的组分汽液相值, 采用( lh ) a v 与这一层组分不对应, 亦必然会导致误差, 在理论上亦说不通。按理, 最小回流比是指采用小于此回流比值, 塔板数就将无穷大, 这不符合有限塔板的精馏情况。有读者已经发现在1999年第2期 深冷技术
[5]
计算R m 公式(1945)
(7a) (7b)
i (x D , i ) m
-Q =R m +1
i
i x F , i )
=1-q
i -Q
中最难挥发的组分为基准)
式中 i 组分i 的相对挥发度(以进料混合物
q 进料的液相分率
x F , i 进料混合物中组分i 的摩尔分率(x D , i ) m 最小回流比下馏出液中组分i
的摩尔分率
Q 方程式的根(对于有C 个组分的系统
将有C 个根, 只取 l
一个根)
计算比较麻烦, 首先要按芬斯克公式由关键组分的分配比及全回流下估算出馏出液的组成, 然后用试差法由式7a 求Q 代入式7b 求R m , 例题见文献[1]中。
3中的表5中只算到氩馏分含
氩8%, 而不象表6中算到7%Ar?这是因为算到7%时, R m 就大于所采用的R =55, 即55小于R m , 应该是理论板数无穷大, 但计算的N 值仍然是可读值, 略高于150, 当时不理解亦无法解释。为了不影响课题主攻方向, 只好把它悬置而略去。现在从理论上详细分析就明白了:用( lh ) a v 所求得的R m 根本不能是这一层的最小值, 更不是入料层以上精馏塔的最小深冷技术 2003年第2期
Cryogenic Technology 2 2003
专题综述
有一点亦是明确的, 所用 i 都是进料组分的 值, 因为只有进料层的x i 、 这与上述认i 是可知的。定R m 应按入料层 反过来说, 只要lh 计算是一致的。( lh ) a v 与( lh ) F 比较接近, 用( lh ) a v 代入计算所得的R m 亦会偏大, 而进一步求得的N 值亦将是低限值。
徐文灏:精馏计算中相对挥发度 的选择
3 简捷法精馏计算 线
[2、4]
吉利兰曲
(1940) 分析与探讨
[4]
吉利兰根据61座两组分和多组分精馏塔的逐层计算结果, 整理得出关联曲线, 利用此曲线可以
用N m 、R m 、R 求出理论塔板数N , 此即简捷法。如图2
所示。
图2 吉利兰曲线
图2中曲线左端延线表示在最小回流比下的操作情况, 此时(R -R m ) /(R +1) 接近零, 而(N -N m ) /(N +1) 接近1(即N = ) 。而曲线右端表示在全回流下操作状况, 此时(R -R m ) /(R +1) 接近1(即R = ) , (N -N m ) /(N +1) 接近零, 即N N m
[7]
所归结的:只要用( lh ) a v 来作计算, 不管什么计算方法, 求得的理论塔板数总是低限值, 而低值总是接近的, 正如图1所示。结果一致的现象由此得到解释。反过来说, 计算结果一致正好证明了用( lh ) a v 去计算塔板数所得的都是完成精馏的低限值, 见表1。
表1 法国6500m 3/h 空分下塔精馏计算结果汇总表
简捷法
下塔段位 - 段 - 段
逐层法812 48
(用各段的R m 计算) 7 5813 08
简捷法(以全塔的R m 计算) 7 88
12 62
实际塔板数(由此可估计板效率)
1414
[3]
。
结合以往计算实例, 可以得出如下一些概念:(1) 既然曲线中用了N m , 则其所用的一定是
( lh ) a v 。在同一计算中, 既然N m 用了( lh ) av , 则可以判定R m 中亦一致地用了( 因此所求得的N 值lh ) av 。必然是完成精馏任务的低限值, 这是以上讨论顺理成章的理解。
(2) 或者正是因为上述吉利兰曲线计算N m 、R m
时都采用( lh ) a v , 就造成了习惯性的计算通用, 再加上前面提到的通常的如果塔顶塔底的 lh 值相差不大时可以用( , 就使概念混淆不lh ) a v 的明确解释清, 掩盖了应该用( 又lh ) F 来计算R m 的正确概念。因为用( lh ) a v 求得的R m 及N 亦是能完成任务的低限值, 不会出偏差。所以就这样长期使用下去, 亦不为过。
(3) 在以往的空分计算中, 用逐层法和简捷法
[3]
计算的结果非常接近。当时虽然很满意, 但亦非常惊奇, 无法解释其一致性。现在看来, 正如以上讨论
[4]
表2 法国6500m 3/h 空分上塔精馏段计算结果汇总表
上塔段位 段(由上向下)
段 段
段 ~ 段全塔
逐层法简捷法87471826
7 967 694 246 3718 326 26
实际塔板数
15
128446479
塔板效率/%
53645014 5
(因有膨胀过热气影响)
28 633 2
(4) 十分明显, 在吉利兰曲线中亦包含了R 与R m 之间有一个比值系数A 。由于用( lh ) av 所求得的R m 已足够大, 如文献[1、5]所计算的那样, 因此此系数值不会很大, 大亦无益。据文献[8]中所写: 根据
2003年第期
2
徐文灏:精馏计算中相对挥发度 的选择
专题综述
间取几个R m 值计算出各自的N 值, 绘成如图1的曲线, 就可以很科学地决定选用的R 及N 值来。因为高限N 值与低限N 值有的相差还是较大的, 不能轻信塔顶与顶底 lh 值相差不大时可采用( lh ) a v 的说法。
(4) 对于精密精馏, 只能用逐层法计算; 简捷法只适用于分离度要求不高, 满足于求得需要的塔板最低限度值的任务。但简捷法由于使用了R m 、N m 等经典公式, 却能带来精馏概念的廓清。
(5) 所有这些讨论都是泛指所有的精馏工况。因此这些结论都适用于各精馏物系。计算实例很好地佐证了这些结论。
自1978年以来, 笔者对空分及氩精馏作了工艺理论计算的系统探讨, 发表了系列论文, 形成完整的计算方法与实例, 其中亦廓清了许多概念。其中N 2 Ar O 2汽液平衡关联式及本文还对基础理论作了发展, 它们都已验证于实验数据与计算实用中。近年来同业还有同志作了另外的开拓, 如在操作线物料平衡中引入了热量平衡, 即进行不等分子回流的计算。精馏是空分装置的核心与主体, 美、德、法、英、日以及俄罗斯的大公司都有自己的计算程序和数据库, 各公司几十年来在精馏上所取得的进展亦是有目共睹, 突飞猛进, 塔板数一加再加, 从而促进和带动了塔板结构、设置的变化乃至稀有气体提取的进展。深切希望同业抓住这个中心继续消化吸收, 独立钻研, 赶上先进并创造出我们自己的理论与实践成果, 带动行业的发展。
参考文献:
[1]徐文灏. 全精馏制氩中回流比R 的抉择与讨论. 深冷技术,
2001(4):1~12.
[2]陈洪舫编. 基本有机化工分离工程. 化工版, 1981:120~121,
125.
[3]徐文灏. 全低压空分装置相对挥发度法精馏计算(上) 逐层
法.
(下) 简捷法. 深冷技术, 1979(2) :1~24; 1979(3) :12~28.
[4]邓修, 吴俊生编著. 化工分离工程. 科学版, 2002:55~56. [5]徐文灏. 低温全精馏制氩工艺计算演绎之一、二. 深冷技术,
1999(1)1~12; 1999(2) 1~9; 1999(3) :16~18. [6]化学工程手册:第三卷. 化工版. P -74、75.
[7]贾绍仪. 蔡诚敬主编. 化工传质与分离工程. 化工版. 2001年
1月第1版:179.
[8]刘家琪主编. 分离工程. 化工版. 2002年5月第一版:107.
Pair J. R 和Bolles W. l: Chem. Eng. 1969(75) :56 的研究结果, R /R m 的最优化值约为1 05, 但是在比值稍大的一定范围内都接近最佳条件。 对于每一物系都可按此考虑, 吉利兰曲线显示在(R -R m ) /(R +1) 中包含了这一因素。从氩精馏来看, 1 05亦是合适的。这与文献[5]中美国专利4994098所采用的1 04~1 06是一致的。
(5) 吉利兰线图绘制的依据是八种物系在广泛的精馏条件下由逐板计算法得出的结果。组分数目2~11个, 进料热状态包括冷料至过热蒸汽等五种情况, R m 为0 53~7 0, 组分间相对挥发度为1 26~4 05, 理论板数为2 4~43 1。可见吉利兰线图适用于石油、苯、酚系等蒸馏要求分离度不高的物系。因为它所求得的亦正是逐层法所求得的最低限度需要塔板数, 此值在R 一定大时, 是趋于一致的。在空分上能成功计算即是因此。用简捷法还解决了过去粗氩塔计算中塔顶95%Ar 含量的粗氩中氮、氧含量随进料组成而变的假定值, 这亦是简捷法的贡献。但对于精密精馏的纯氩计算则只能用逐层法, 唯逐层法才能分别用( lh ) F 和( lh ) av 判别任务需要的最高和最低塔板极限
[1]
[5][5]
[7]
。
4 结语
通过这样以实际计算结果佐证的对理论的探讨, 是否可以归结为以下结论:
(1) 由于全回流(R ) 下求N m 必须采用( lh ) a v 。因此可以推论:用( lh ) a v 在逐层法或简捷法中计算的R m 值已足够大, 所计算的N 值必然是精馏所需的最低限度塔板数。但此值不是全回流下的N m 值, 而是实际生产所需最低值。
(2) 用入料层的( lh ) F 求得的R m 才是该精馏任务的最小回流比, 小于此R m , N , 精馏不能进行。采用略大于R m 的R 计算出来的塔板数N 才是实际需要的最多塔板数, 实际上此值是偏高的, 因为进料层以上塔板的( lh ) i 都比( lh ) F 小, 这导致计算中最上的多层塔板纯度虽稍有提高, 但徙增塔板阻力。
(3) 对于纯产品的精密精馏, 应该用( lh ) F 计算出R m 、R (R 可采用1 05R m 左右) , 再用逐层法计算出高限N 值来; 再用( R , 计算出低限lh ) a v 求R m 、N 值来。借助计算机程序计算, 可以在R m 与R m 之深冷技术 2003年第2期
Cryogenic Technology 2 2003
5