气液相平衡相是指系统的某一部分具有相的物理和化学性质,具有相同的组成,并且与另外的相以一定的边界隔开。
出现在世有产品里的大多数的相相是液态烃和气相。
水也是普遍存在的另一个液态相。
在给定的系统里,当描述变化的变量随着时间和位置的改变而保持恒定时,液态烃、气相和水这些相将会平衡共存。
而决定平衡状态的主要变量是系统的温度,压力和组分。
对于设计表面分离装置和改进组分模型,不同相能够共存的条件是非常值得考虑并且具有实践价值的重要问题。
这些计算的类型是以平衡常数的原理为基础的。
一、平衡常数一给定组分的平衡常数是指组分的气相的摩尔分数yi与液相摩尔分数xi之比。
在数学上,它们的关系是(5----1) 在(100磅/)的压力下,Raoult和Dalton的定律为对于理想溶液方法提供了一个预测平衡常数的简化的方法。
Raoult定律是指多组分系统中单一组分所产生的局部压力等于它的液相摩尔分数与该祖父的气相压力的乘积即=(5----2) ------组分i的局部压力,磅/)------组分i的气相压力,磅/)------组分i的液相摩尔分数Dalton定律是指某一组分的局部压力等于它的气相摩尔分数与系统的总压力的乘积,即=P (5----3) P------系统总压力,磅/在平衡状态下,根据上面的定律可知,被某一组分气相作用产生的局部压力必须与该液相作用产生的局部压力平衡。
因此,将描述两大定律的方程结合得=P ,整理以上关系式并代入平衡常数定义式得/=/P=(5----4) 这个方程表明对理想溶液,不管烃类混合物的组分如何,平衡常数仅仅有体现系统压力和温度的作用。
(第一章表明,组分的气相压力有体现温度的作用。
)现阶段,介绍并给以下的术语下定义是很有必要的。
------给定的烃类混合物中组分i的摩尔分数,n ------烃类混合物的总摩尔数,------液相的摩尔总数,------气相的摩尔总数。
由定义可知,(5----5) 这个方程表明系统总的摩尔数与气相的摩尔数和液相的摩尔数之和是相等的。
已知组分的物质平衡的结果:=+ (5----6) 在这里------组分i在系统里的总摩尔数,------组分i在液相里的总摩尔数,------组分i在气相里的总摩尔数。
也根据摩尔百分数的定义,我们可以写成=1 (5----7)=1 (5----8)=1 (5----9) 以烃类混合物的一摩尔为基础进行全部的相平衡计算是方便的,即n=1。
该假设简化方程(5----5)和(5----6)为+=1 (5----10)=+ (5----11) 联合方程(5----4)和(5----11)以便从方程(5----11)中消掉得到=+ ,解出的值= /(+) (5----12) 通过联合(5----11)和(5----4)方程消掉也能从(5----11)中求出的值,得到=/(+) (5----13)再联合(5----12)和(5----8)及(5----13)和(5----9)方程得=/(+)=1 (5----14)=/(+)=1 (5----15) 已知-=0,即/(+)-/(+)=0,那么/(+)=0,以(1-)代替得=(5----16) 二、闪蒸计算闪蒸计算是指所有油层和加工工程计算的总和。
在一定的压力和温度下,闪蒸计算需要了解在一个油箱或者在海上钻井装置中共存的液态烃和气态烃的数量,无论什么时候这样做都是很有必要的。
它也可以用来确定共存烃类相态的组成。
闪蒸计算前需要确定的是:1>气相的摩尔分数;2>液相的摩尔分数;3>液相的组成;4>气相的组成.计算的步骤是确定已知全部组分的烃类混合物的, ,和,即可确定,以及确定一系列可计算的平衡常数,如,然后把它们记录下来。
步骤:1>的计算用Newton---Raphson迭代法可解出方程(5----16)中的。
下面是这种迭代法的运用:a.假定的值。
好的假设值是靠以下的关系式估算的=A/(A-B) ,A =[(-1)],B =[(-1)/ ]b.按照方程(5----16)通过使用假设值而得到的值估算出的结果。
c.如果绝对值比预定的限度小(例如10),那么假设的就是所需要的解。
d.如果的绝对值比预定的限度大,那么新的就要有以下的表达式估算:()=-/()----新的的值-----关于的的第一个导数,可由={ (-1)/[(-1)+1]} 得到。
e.以上步骤一直迭答下去直到一个收敛被获得。
2>计算从给出的(5----10)方程可以计算出液相的摩尔数=1-3>计算通过运用方程(5----12)可计算液相的组成4>计算由方程(5----13)可以决定气相的组成例5—1由如下成分组成的碳水化合物在50磅/和100F条件下进行分离。
假设一个理想的解决方法,求分离后的各组成?解:1>从图表中(见第一章)可知蒸汽压力并根据方程(5----4)计算出平衡常数,组成100F 的压力C0.20 190 3.80i ----C0.10 72.2 1.444n---- C0.10 51.6 1.032i ----C0.20 20.44 0.4088n----C0.20 15.57 0.31140.20 4.956 0.099122>通过解方程(5----16)得到迭代次数0 0.08196579 0.073E-021 0.1079687 8.894E-042 0.1086363 7.60E-073 0.1086368 1.49E-084 0.1086368 0.03>求解=1-0.1086368=0.89136314>求解和组成方程(5----12)方程(5----13)C0.1534 0.5827i ---- C0.0954 0.1378n---- C0.0997 0.1028i ---- C0.2137 0.0874n---- C0.2162 0.06730.2216 0.0220三、平衡常数平衡常数是组分在液相与气相中的比例。
根据方程(5----4)来计算,还需要蒸汽压力和系统压力。
字运用方程(5----4)时有一些基本的假设:1>气相是理想气体并满足达西定理2>液相有一个理想的溶解度并Raoult定理以上的假设是不契合实际的,并且在高压条件预测的平衡常数是不准确的。
为了得到真实解,平衡常数不仅仅是压力和温度的函数,而且还是炭氢合物组成的函数。
这个关系可用以下的的数学式表达式表示:=有许多方法被建议用来进行平衡常数的计算,一些是简单的数学表达式。
一些是复杂的表达式包含几种独立的不同的组成。
这些方法以后再介绍。
威尔逊关系威尔逊在1968年提出了一个简单的热力学表达式来求K值这个表达式如下:=EXP[5.37(1+)(1-/T)] (5----17)------成分i的临界压力P ------系统压力------成分i的临界温度T ------系统温度------成分i的热容系数以上这个关系在低压情况下所计算出的平衡常数是你叫精确的。
斯坦丁关系斯坦丁在1979年提出了一套计算平衡常数的公式,在1000pisa以下和温度在200F以下的情况比较适合。
成分的特征因素在给定的压力情况下呈现直线关系,其数学表达式如下:=(1/-1/) (5----18)其中,=log(/14.7)/ (1/-1/) (5----19)------成分i的特征系数------成分i的沸点七直线关系的基本方程为: log(P)=a+C或=10(5----20)这里系数a、c为直线的截距和斜率。
通过对等压线和18个方程的估计系数a、c与压力的关系由下式给出a=1.2+0.00045P+P(5----21)c=0.89-0.00017P-(5----22)斯坦丁指出:N、CO、H S 、到的烃类的平衡常数,通过这些成分的沸点和相关参数能得到很好的证明。
作者建议如下的的修正值:成分RN470 109CO652 194H S 1136 331300 94C1145 303C1779 416i ---- C237 471n---- C2153 491i ---- C2368 542n---- C2480 5572738 610当进行闪蒸计算时,对附加部分(除以上部分的烃类)的平衡常数的问题总是被提出。
凯特在1937年提出C以上的烃类的K值有15%的能用C的K 值代替.斯坦丁提出了一个可供选择的方法,通过实验与方程(5----20)的集合来求解平衡常数(C以上部分),并计算出了特征系数。
斯坦丁建议使用如下的计算机程序来决定参数b和:第一步:由下面的关系式来决定正常石蜡烃类的碳原子数n=7.30+0.0075(T-460)+0.0016P (5----23) 第二步:由下面的表达式来计算参数b和沸点b=1013+324n-4.256(5----24)=301+59.85n-0.971(5----25) 以上的计算值可代入方程(5----18)来求。
例5-2烃饿里混合物的组成如下,在1000pisa和150F的情况下进行闪蒸分离:成分CO0.009N0.0030.535C0.115C0.088i ----C0.023n---- C 0.023i ---- C0.015n---- C0.0150.015C0.159假设C以上烃类的分子量和比重为198.71和0.8527计算平衡常数分别使用:a.威尔逊关系; b.斯坦丁关系解:a. 威尔逊关系第一步:通过使用第一章例2-6的方法来计算临界压力、临界温度和,=1294.1R,=263.67psia, =0.5346第二步:运用方程(5----17)求解成分,psia ,RCO1071 0.225 547.9 2.0923N493 0.040 227.6 16.343C667.8 0.0104 343.37 7.155C707.8 0.0986 550.09 1.263C616.3 0.1524 666.01 0.349i ----C529.1 0.1848 734.98 0.144n----C550.7 0.2010 765.65 0.106i ----C490.4 0.2223 829.1 0.046n----C488.6 0.2539 845.7 0.036C436.9 0.3007 913.7 0.013C263.67 0.5346 1294.1 0.0255b.斯坦丁关系第一步:计算系数a和c[由(5----21)和(5----12)得到]a=1.65, c=0.685第二步:用(5----23)计算碳原子数,n=10.025第三步:运用方程(5----24)和(5----25)计算参数b和沸点Tb, b=3833.369, =803.41R第四步:运用方程(5----20)的结果成分CO652 194 2.292 1.660N470 109 3.541 11.901C300 94 2.700 3.159C1145 303 1.902 0.897C1799 416 1.375 0.391i ----C2037 471 0.985 0.211n----C2153 491 0.855 0.172i ----C2368 542 0.487 0.096n----C2480 557 0.387 0.082C2738 610 0 0.045C3833.369 803.41 -1.513 0.004集中收敛方法在早期的高压相平衡研究中我们发现,当定组分的烃类混合物在恒温情况下,随着压力的增加,所有的组分的平衡值都会向一点收敛。