《气液相平衡方面的基础知识汇总.docx》由会员分享,可在线阅读,更多相关《气液相平衡方面的基础知识汇总.docx(25页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、 气 液 相 平 衡 相是指系统的某一局部具有相的物理和化学性质,具有一样的组成,并且与另外的相以确定的边界隔开。出现在世有产品里的大多数的相相是液态烃和气相。水也是普遍存在的另一个液态相。在给定的系统里,当描述变更的变量随着时间和位置的变更而保持恒定时,液态烃, 气相和水这些相将会平衡共存。而确定平衡状态的主要变量是系统的温度,压力和组分。 对于设计外表别离装置和改良组分模型,不同相能够共存的条件是特殊值得考虑并且具有实践价值的重要问题。这些计算的类型是以平衡常数的原理为根底的。一、 平衡常数 一给定组分的平衡常数是指组分的气相的摩尔分数与液相摩尔分数之比。在数学上,它们的关系是 (51)
2、在 EMBED Equation.3 (100磅/)的压力下,和的定律为对于志向溶液方法供应了一个预料平衡常数的简化的方法。定律是指多组分系统中单一组分所产生的局部压力等于它的液相摩尔分数与该祖父的气相压力的乘积即= (52) 组分i的局部压力,磅/) 组分i的气相压力,磅/) 组分i的液相摩尔分数 定律是指某一组分的局部压力等于它的气相摩尔分数与系统的总压力的乘积,即 EMBED Equation.3 (53)系统总压力,磅/ 在平衡状态下,依据上面的定律可知,被某一组分气相作用产生的局部压力必需与该液相作用产生的局部压力平衡。因此,将描述两大定律的方程结合得 EMBED Equation.
3、3 ,整理以上关系式并代入平衡常数定义式得 EMBED Equation.3 EMBED Equation.3 (54) 这个方程说明对志向溶液,不管烃类混合物的组分如何,平衡常数仅仅有表达系统压力和温度的作用。第一章说明,组分的气相压力有表达温度的作用。 现阶段,介绍并给以下的术语下定义是很有必要的。 给定的烃类混合物中组分i的摩尔分数, n 烃类混合物的总摩尔数, 液相的摩尔总数, 气相的摩尔总数。 由定义可知, (55) 这个方程说明系统总的摩尔数与气相的摩尔数和液相的摩尔数之和是相等的。 组分的物质平衡的结果: EMBED Equation.3 EMBED Equation.3 EMB
4、ED Equation.3 (56) 在这里 EMBED Equation.3 组分i在系统里的总摩尔数, EMBED Equation.3 组分i在液相里的总摩尔数, EMBED Equation.3 组分i在气相里的总摩尔数。 也依据摩尔百分数的定义,我们可以写成 =1 (57) =1 (58) =1 (59) 以烃类混合物的一摩尔为根底进展全部的相平衡计算是便利的,即1。该假设简化方程55和56为 + =1 (510) EMBED Equation.3 (511) 联合方程54和511以便从方程511中消掉得到 EMBED Equation.3 EMBED Equation.3 ,解出的
5、值 = /(+) (512) 通过联合511和54方程消掉也能从511中求出的值,得到 EMBED Equation.3 (+) 513 再联合512和58及513和59方程得 EMBED Equation.3 EMBED Equation.3 (+)=1 (514) EMBED Equation.3 EMBED Equation.3 (+)=1 (515) EMBED Equation.3 - =0,即 /(+) EMBED Equation.3 (+)=0,那么 EMBED Equation.3 /(+)=0,以(1-)代替得= EMBED Equation.3 (516)二、 闪蒸计算闪
6、蒸计算是指全部油层和加工工程计算的总和。在确定的压力和温度下,闪蒸计算须要了解在一个油箱或者在海上钻井装置中共存的液态烃和气态烃的数量,无论什么时候这样做都是很有必要的。它也可以用来确定共存烃类相态的组成。闪蒸计算前须要确定的是:1气相的摩尔分数;2液相的摩尔分数;3液相的组成;4气相的组成. 计算的步骤是确定全部组分的烃类混合物的, ,和,即可确定,以及确定一系列可计算的平衡常数,如,然后把它们记录下来。 步骤:1的计算 用迭代法可解出方程516中的。下面是这种迭代法的运用:的值。好的假设值是靠以下的关系式估算的 () , A = (-1), B = (-1)/ b.依据方程(516)通过运
7、用假设值而得到的值估算出的结果。 确定值比预定的限度小例如10,那么假设的就是所须要的解。的确定值比预定的限度大,那么新的就要有以下的表达式估算:() EMBED Equation.3 EMBED Equation.3 ()新的的值 关于的的第一个导数,可由 = (-1) /(-1)+1 得到。 e.以上步骤始终迭答下去直到一个收敛被获得。 2计算 从给出的(510)方程可以计算出液相的摩尔数=1- 3计算 通过运用方程512可计算液相的组成4计算 由方程513可以确定气相的组成例51 由如下成分组成的碳水化合物在50磅/和100F条件下进展别离。假设一个志向的解决方法,求别离后的各组成? 解
8、:1从图表中见第一章可知蒸汽压力并依据方程54计算出平衡常数,组 成100F 的压力 C 190 i 0.10 72.2 1.444 C 0.10 51.6 i 0.20 20.44 0.4088 0.20 0.20 4.956 0.09912 2通过解方程516得到迭代次数 0 12 3 4 0.0 3求解 组 成方程512方程513 Ci C Ci C C4求解和 三, 平衡常数 平衡常数是组分在液相与气相中的比例。依据方程54来计算,还须要蒸汽压力和系统压力。字运用方程54时有一些根本的假设:1气相是志向气体并满足达西定理2液相有一个志向的溶解度并定理以上的假设是不契合实际的,并且在高压
9、条件 预料的平衡常数是不精确的。 为了得到真实解,平衡常数不仅仅是压力和温度的函数,而且还是炭氢合物组成的函数。这个关系可用以下的的数学式表达式表示:= 有许多方法被建议用来进展平衡常数的计算,一些是简洁的数学表达式。一些是困难的表达式包含几种独立的不同的组成。这些方法以后再介绍。 威尔逊关系 威尔逊在1968年提出了一个简洁的热力学表达式来求K值这个表达式如下: EMBED Equation.3 5.37(1+)(1) 517 成分i的临界压力 P 系统压力 成分i的临界温度 T 系统温度 成分i的热容系数 以上这个关系在低压状况下所计算出的平衡常数是你叫精确的。 斯坦丁关系 斯坦丁在197
10、9年提出了一套计算平衡常数的公式,在1000以下和温度在200F以下的状况比拟适合。 成分的特征因素在给定的压力 状况下呈现直线关系,其数学表达式如下: =(11/) 518 其中,(/14.7)/ (11/) 519 成分i的特征系数 成分i的沸点 七直线关系的根本方程为: (P) 或 =10 520这里系数a, c为直线的截距和斜率。 通过对等压线和18个方程的估计系数a, c与压力的关系由下式给出 EMBED Equation.3 (521) (522) 斯坦丁指出:N, , , 到的烃类的平衡常数,通过这些成分的沸点和相关参数能得到很好的证明。作者建议如下的的修正值:成 分 EMBED
11、 Equation.3 RN4701096521941136331 30094C1145303C1779416i C237471 C2153491i C2368542 C24805572738610 当进展闪蒸计算时,对附加局部除以上局部的烃类的平衡常数的问题总是被提出。凯特在1937年提出C以上的烃类的K值有15%的能用C的K值代替. 斯坦丁提出了一个可供选择的方法,通过试验与方程520的集合来求解平衡常数C以上局部,并计算出了特征系数。 斯坦丁建议运用如下的计算机程序来确定参数b和: 第一步:由下面的关系式来确定正常石蜡烃类的碳原子数 7.30+0.0075(460)+0.0016P (5
12、23) 第二步:由下面的表达式来计算参数b和沸点 (524) EMBED Equation.3 (525) 以上的计算值可代入方程518来求。 例 5-2 烃饿里混合物的组成如下,在1000和150F的状况下进展闪蒸别离: 成 分 NCCi Ci C CC 假设C 解:a. 威尔逊关系 第一步:通过运用第一章例2-6的方法来计算临界压力, 临界温度和, EMBED Equation.3 263.67, 第二步:运用方程517求解成 分 EMBED Equation.3 1071N493CCCi i EMBED Equation.3 CC 第一步:计算系数a和c由521和512得到 第三步:运用
13、方程524和525计算参数b和沸点, 3833.369, EMBED Equation.3 R 第四步:运用方程520的结果 成分652194N470109C30094C1145303C1799416i 20374712153491i 23685422480557C27386100C 集中收敛方法 在早期的高压相平衡探讨中我们发觉,当定组分的烃类混合物在恒温状况下,随着压力的增加,全部的组分的平衡值都会向一点收敛。我们把这个压力叫做改混合物的收敛压力P。在求组分的平衡常数时,收敛压力是特殊有必要的。图51:平衡常数与压力的相对关系图如上。 收敛压力理论在图中表示得特殊明确,图表示了在常温状况下
14、,在双对数坐标纸上,烃类混合物的平衡常数与压力的关系。对全部组合在一个详细的压力即收敛压力状况下,从图上可以看出平衡常数逐渐的趋向于一个常量=1。对不同的烃类混合物可能会出现不同的收敛压力。 自然气程序运用协会作出了K值与温度, 压力等组成和收敛压力的函数关系。这个关系图可从的工程数据丛书中获得,书中还包括了许多收敛压力状况下的K值,详细为800,1000,1500,2000,3000,5000和10000。 几位专家调查发觉,烃类混合物的收敛压力在4000或者以上时,平衡常数的值与烃类混合类的系统压力在1000以下的时候特殊相像。这个结论说明白,在1000以下时烃类混合物的组成对平衡常数的影
15、响很小。 在运用自然气程度协会举荐的平衡常数的解法时存在的问题就是在选择适合的图之前就必需知道收敛压力。下面介绍三种求解收敛压力的方法。 A:哈顿方法:哈顿在1956年提出了一个求解烃类混合物收敛压力的方法。这个方法是基于“假设和“等价两个概念,运用如图52的收敛压力图。方法可依据为如下几步: 第一步:为收敛压力估计一个值; 第二步:从一个适合的平衡常数图上读取一个K值,这个K值是依据每个成分目前的系统压力和温度而定的; 第三步:利用估计的K值和系统组成进展闪蒸计算; 第四步:在0.1液相当中找出最高的烃类成分; 第五步:将液体的摩尔组成转化为质量组成; 第六补:除开最高成分后,重复第四步,并
16、求出主要成分的标准质量组成。 第七步:计算平均重量,临界温度和压力;由下式计算 = = EMBED Equation.3 其中是成分i的标准质量组成。 第八步:由假设成分所得的临界参数返回到图5-2; 第九步:读取新的收敛压力的值; 第十步:假设新的收敛压力的值是不够志向的重复步骤二到步骤九。见附图5-2 这里应当留意:假设计算出的新的收敛压力的值是介于图中所得的值之间时,对图的说明特殊必要的。假设K值的变更不大而P时,此时值就可以作为真实值被用。 B:斯坦丁方法:斯坦丁在1977年提出收敛压力与分子量有一个粗略的直线关系,后来,威特斯顿用如下的关系式来表达了这一思想:=60M*-4200 其
17、中,M*表示重烃C以上的局部的分子量。 C:雷撒方法:雷撒在1952年提出了一个对轻烃的收敛压力的一个简洁的递推关系。他运用温度和C以上局部的分子量和重度作为相关参数。这一关系可用图5-3来表示附图5-3。关系式可用如下的数学表达式来表示: 2381.8542+46.341487(M*) EMBED Equation.3 (M*) /(460) (527) 其中:M* C的分子量 C的重度 EMBED Equation.3 = 为如下的常数 以上的数学表达式适合于求解烃类混合物的温度在50F到300F时的收敛压力。 托普和威特森关系 托普和威特森在1981年对威尔逊方程517做了改良,使它能够
18、 用于高压计算。威尔逊的方程被修改,增加了收敛压力在方程中,关系如下: = (/) ()*5.37A(1+) (1) 1- (528) 其中:P 系统压力 收敛压力 T 系统温度 成分i的热容系数例5-3重复了5-2的工作,并且计算了平衡常数运用了托普方法: 解:第一步:从方程527得收敛压力 第二步:由方程528计算平衡常数成 分NCCCi CC 纹恩关系 纹恩在1954年提出了求解收敛压力在5000烃类混合物的平衡常数的方法。如右图5-4该图由温度, 压力的网格, 收敛压力刻度, 平衡常数刻度组成。烃类混合物的平衡常数可以干脆从图上读得,K值的刻度是通过志向成分的直线与系统温度和压力的网格
19、的交点。 对于收敛压力不在5000的烃类混合物,纹恩供应了一个校正图版,它表示了系统压力, 收敛压力和网格压力之间的关系见图5-5。网格压力是用于连接图5-4获得的K值不在5000时,校正过程可简洁为以下几步: 第一步:在图5-4中描出系统温度, 压力的点; 第二步:在图5-4中把刚刚作的点与平衡常数刻度线上的1.0点用直线连接起来; 第三步:在这条直线上找出临界点和图5-5中相关的压力直线的交点就是不同于系统温度, 压力的网格点; 第四步:读取平衡常数的刻度。它是通过从第二步的网格点直线和志向组成的交点。 状态方程 总体来说,上面所述的方程为计算平衡常数供应了更多的阅历和技术理论根底,我们将
20、在第六章对它进展更加详细的论述。四, C以上局部的平衡常数 C以上局部的平衡常数与其它成分的值经常是不同的,这是因为C的局部自身就是一种混合物。 有几种技术可以用来估算C以上局部的K值以下介绍几种: 坎波贝儿的方法 坎波贝儿在1973年发觉对任何烃类系统,K的对数值与任何成分的 EMBED Equation.3 成直线关系。因此坎波贝儿提出作一条好的回来直线。通过这条直线C以上局部的K值就可以得到。坎波贝儿还指出,K的对数值与1/也呈直线关系。通过这条直线可以得到C以上局部在它的平均泡点压力时的平衡常数。坎波贝儿的思想在图5-6和5-7中得以说明。 纹恩的方法 纹恩在1954年提出了对重烃局部
21、且露点在210F以上时平衡常数的计算式: EMBED Equation.3 ( EMBED Equation.3 ) (529) 其中: C以上局部的K值; 在系统压力, 温度和收敛压力下的K值 b 粘度 纹恩将重质组分中的挥发组分的相关性质,如露点,绘成图表见图5-8。其数学表达式如下: EMBED Equation.3 ( -460)(-460) (-460) (-460) (530) 其中:露点 EMBED Equation.3 EMBED Equation.3 EMBED Equation.3 EMBED Equation.3 卡茨的方法 卡茨(1957)建议对于丁烷及其以上的重质组分
22、,其平衡比例许再乘以才显得比拟合理,相关数学表达式见如下方程: 0.15K K 531 其中: EMBED Equation.3 的平衡常数 K 丁烷的平衡常数 霍夫曼等的方法 霍夫曼, 克努波, 郝考特等指出平衡常数的对数与系统的确定压力有关,而系统的确定压力与特征系数的反比关系。特征系数定义为如下的关系式: = 532 其中,成分i的露点 成分i的临界温度 T 系统温度 例5-4 烃类系统的温度和压力分别为201和1000,有试验得局部的K值外其它的K值如下所示:组 分i i 0142 局部的分子量和详细重度为198.7和0.841。假设混合物的收敛压力为5000,试用以下方法计算的平衡常
23、数: a: 运用斯坦丁关系;b: 运用托普关系;c: 纹恩方程;d: 卡茨方程;e: 霍夫曼方法.解:首先计算临界参数,运用第二章的方程21,求,和 EMBED Equation.3 , EMBED Equation.3 , 运用上面所得的参数再代入方程221计算得:a: 运用斯坦丁关系:第一步:运用方程521和522得参数a和c,第二步:由方程520计算K EMBED Equation.3 b: 运用托普关系: 第一步:运用方程528计算A值: 1- 第二步:由方程528计算K,得c: 纹恩方程: 第二步:求K和 K,得1.52 , 第三步:由方程529计算K,即 EMBED Equatio
24、n.3 d: 卡茨方程:代入方程531得, EMBED Equation.3 e: 霍夫曼方法 第一步:计算()和特征函数组 分47691515650i 360295i 165142第二步:计算()/找出直线上对应的点图5-9如右,读出即第三步:求解K, 相平衡模型组成:1描述平衡液相组成, 物质量摩尔数及平衡常数与温度, 压力关系的物料平衡条件方程组; 2描述平衡液相组成, 物质量, 平衡常数与逸度关系的热力学平衡条件方程组及用语相平衡计算的状态方程机构体系。一, 油气体系体系气液相平衡计算烃类气液平衡体系,气相摩尔重量V,液相摩尔重量L,1,平衡常数是指组分的气相的摩尔分数与液相摩尔分数之
25、比, EMBED Equation.3 =1, EMBED Equation.3 =1。在数学上,它们的关系是 。平衡气相质量守恒方程为: EMBED Equation.3 气相组成方程为: = EMBED Equation.3 气象平衡方程组为: EMBED Equation.3 EMBED Equation.3 1液相组成方程为:=液相物料平衡方程组:1气液两相总物料平衡方程为:-= EMBED Equation.3 二, 油气体系相态计算热力学平衡依据热力学相平衡理论确定出组分气, 液相的逸度系数,进而确定平衡常数。依据流体热力学平衡理论,当油气烃类体系到达气液两相平衡时,体系中各组分在
26、气液相中的逸度和应满足=。体系气液相逸度的表达式分别为 = EMBED Equation.3 P = EMBED Equation.3 P联合以上两式得, EMBED Equation.3 EMBED Equation.3 EMBED Equation.3 其中 EMBED Equation.3 , EMBED Equation.3 平衡气液相中组分的逸度系数。三, 常用的两种类型相平衡计算数学模型1、 第一种类型的相平衡计算数学模型在油气藏烃类体系相平衡计算中,第一类数学模型是基于适合接受逐步迭代算法或牛顿干脆算法而建立的。而依据油气体系相态分析的须要,可分别建立适合于油气体系露点, 泡点和
27、等温闪蒸给定温度,局部汽化或局部液化的相态变更过程计算的相平衡条件方程组。1) 适合于露点计算的相态模型当体系相态变更处于露点状态时,起平衡的物质量平衡关系表现为平衡气相的摩尔重量V1;平衡液相的摩尔重量L0,有 EMBED Equation.3 EMBED Equation.3 EMBED Equation.3 EMBED Equation.3 归一化条件那么为 EMBED Equation.3 EMBED Equation.3 1 联合以上方程可得到露点计算目标函数方程组为 EMBED Equation.3 0 EMBED Equation.3 0 EMBED Equation.3 0 1
28、 2) 适合泡点计算的相态模型 同理可得,V0,L1 EMBED Equation.3 EMBED Equation.3 EMBED Equation.3 EMBED Equation.3 EMBED Equation.3 1于是,得到泡点计算平衡条件目标方程组为 EMBED Equation.3 0 EMBED Equation.3 0 EMBED Equation.3 0 EMBED Equation.3 EMBED Equation.3 13适合于等温闪蒸计算的相态模型 油气体系的相态变更处于局部汽化和液化的状态时,平衡气液相的摩尔重量分别在0和1之间变更,由以上方程组合可得 EMBED
29、 Equation.3 0 EMBED Equation.3 0 EMBED Equation.3 0 EMBED Equation.3 EMBED Equation.3 =0因此,可分别求解所对应相态中的1个变量。 用压力和热力学逸度表示物料热力平衡条件方程,更直观的反映出热力学参数对相平衡的影响: 对露点计算 EMBED Equation.3 EMBED Equation.3 EMBED Equation.3 0对泡点计算 EMBED Equation.3 EMBED Equation.3 EMBED Equation.3 0 牛顿迭代法求解非线性方程组过程中的物料平衡方程为 EMBED
30、Equation.3 EMBED Equation.3 =0 平衡过程计算 平衡计算是自然气加工工程中经常用到的过程计算,这种计算包括以下四类根本问题:1) 确定压力P或温度T下,液相组成x,确定泡点温度或泡点压力和平衡的气相组成y,称为泡点计算。2) 确定压力P或温度T下,气相组成y,确定泡点温度或泡点压力和平衡的气相组成x,称为露点计算。3) 料液的量F和组成,计算在确定压力P和温度T闪蒸得到的气液两相的量V,L和组成,称为等温闪蒸计算。4) 料液的量F和组成,其压力为P,温度为T,在绝热状况下瞬间降压至,系统温度降至,这个过程为绝热闪蒸过程。其中,泡点压力和露点压力的计算比拟简洁,只须要
31、求得规定温度下各组分得饱和蒸汽压,即可干脆计算出P和或。泡点计算和露点计算相像,只涉及相平衡关系,而闪蒸计算除涉及到相平衡关系外,海涉及到物料衡算和热量衡算。一、 泡点计算给定的液相组成x和总压P,计算确定液体刚开场沸腾的温度即泡点和产生的气相平衡组成y。计算用到的方程有:1. 相平衡关系 = EMBED Equation.3 ,i2. 气液相组成加和归一化方程 EMBED Equation.3 EMBED Equation.3 1, EMBED Equation.3 =13. 相平衡常数关联式 =, i共有方程数22个方程,变量数为32,故自由度数=32-22, 现给定P和1个,那么上述方程
32、组有唯一解。因为上述方程和T式非线性的需以迭代法求解。二、 露点计算露点温度式确定组成y的气相在确定压力P下冷凝,当开场凝出第一滴液珠时的温度。计算用到的方程有:1. 相平衡关系 = EMBED Equation.3 ,i 2. 气液相组成加和归一化方程 EMBED Equation.3 EMBED Equation.3 1, EMBED Equation.3 =1 3. 相平衡常数关联式 =, i 三, 闪蒸计算 在泡点计算和露点计算之间的温度, 压力下,存在着两相混合物,每一相的量和组成取决于对体系有影响的条件,这些条件最普遍的组合就是固定温度T或固定H和P,通过物料平衡和相平衡来进展计算
33、。1) 等温闪蒸的计算气相混合物的局部冷凝和液体混合物的局部汽化均属于等温闪蒸,此时进料的组成和量是的,要求计算在指定压力和温度下产生的气液两相的组成和量。三、 对每一组分的物料平衡算和平衡条件为 EMBED Equation.3 ,= EMBED Equation.3 式中: F 进入的原料量,; L 液相的量,; V 气相的量,; 进料的摩尔组成; 液相的摩尔组成; 气相的摩尔组成; 引入汽化分率,闪蒸条件即目标函数就变为 1 EMBED Equation.3 1 EMBED Equation.3 0 其相应的算法格式为 EMBED Equation.3 通过逐次靠近法求得e之后,以1为初值,用这种方法常可求得一个收敛的解,= 和运用了一个常引起更快收敛而对e的初值无特殊要求的目标函数: EMBED Equation.3 EMBED Equation.3 0 闪蒸方程:= EMBED Equation.3 其算法格式为 EMBED Equation.3 在应用以上两种算法格式进展汽化, 冷凝计算前宜先推断混合物在指定的温度压力下是否处于两相区,其判据如下: =1 进料处于泡点,e = 0 EMBED Equation.3 1 T e 0 1 T1 T e 1 进料为过热蒸汽
限制150内