LNG液化天然气Chapter03气液相平衡课件

,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,College of Chemistry & Chemical Engineering,SWPU,*,Southwest Petroleum University,3 Vapour Liquid Equilibrium,Phase Behavior of Natural Gas,Equation of State,Gas-liquid Equilibrium,1,College of Chemistry & Chemical Engineering,SWPU,3 Vapour Liquid EquilibriumPh,3.1 Phase Behavior of Natural Gas,一、基本概念,1、相(Phase),内部物理性质和化学性质完全均匀的体系。,在不同相之间，有明显的相分界面；,同一相物质可成片存在，也可以孤立的泡状、滴状等存在；,一相中可含有多种组分。,水(H,2,O),Oil,(C,7,、C,8,、),2,College of Chemistry & Chemical Engineering,SWPU,3.1 Phase Behavior of Natural,3.1 Phase Behavior of Natural Gas,2、相态,含一种或多种物质的单一介质，,p,、,T,、,V,之间的关系可用连续方程表达，则该介质为同一相态。,利用克拉贝伦方程解释相态转变：,式中：,V,为体积改变量；,H,为相变的热焓。,当压力、温度恒定而体积和热焓发生有限变化时，相变发生。,3,College of Chemistry & Chemical Engineering,SWPU,3.1 Phase Behavior of Natural,3.1 Phase Behavior of Natural Gas,3、相平衡(Phase Equilibrium),在,p,、,T,一定时：,A相分子进入B相的量B相分子进入A相的量时的状,4、自由度(Degree of Freedom),在不改变平衡体系中原有平衡相数的条件下可独立改变的量(如,p,、,T,和浓度等)的数目。,5、相态种类及其特点,(1)气相(Gas Phase),4,College of Chemistry & Chemical Engineering,SWPU,3.1 Phase Behavior of Natural,3.1 Phase Behavior of Natural Gas,特点是其占据容器整个空间且易被压缩。,(2)液相(Liquid Phase),形状与容器底部相同，难于被压缩，仅在极高压力下可被压缩。,(3)固相(Solid Phase),有固定形状而不能被压缩。,可能存在几种液相和固相，但气相只有一种。天然气处理与加工中通常只需考虑液相和气相。,5,College of Chemistry & Chemical Engineering,SWPU,3.1 Phase Behavior of Natural,3.1 Phase Behavior of Natural Gas,6、饱和蒸汽压(Vapor Pressure),蒸汽,液体,在一定,T,下，液体与在液面上的蒸气呈平衡状态时，此蒸气所产生的压力。,体现为气相分子对器壁的压力。,7、沸点(Boiling Point),当液体的蒸汽压力等于外加压力时，液体内部产生气泡，并不断地冒出而沸腾的温度。,8、泡点(Bubble Point),开始从液相分离出,第一个气泡,的气液共存态。,6,College of Chemistry & Chemical Engineering,SWPU,3.1 Phase Behavior of Natural,3.1 Phase Behavior of Natural Gas,泡点温度：,在,p,一定的时，开始从液相中出现第一个气泡的温度。,9、露点(Dew Point),开始从气相中凝结出第一滴液滴的气液共存态。,露点温度：,在,p,一定的时，开始从气相中凝结出第一滴液滴的温度。,二、相态研究的工具,1、相律(Phase Rule),7,College of Chemistry & Chemical Engineering,SWPU,3.1 Phase Behavior of Natural,3.1 Phase Behavior of Natural Gas,对,N,相,C,个组分的烃类非均匀封闭体系，自由度F为：,F,C,N,2,2、相图(Phase Diagram),直观的相态研究和表示方法。,表示相平衡态与物系,x,、,T,和,p,等状态变量之间的关系图，亦称相平衡状态图。,3、相平衡方程(Equation of Phase Equilibrium),泡点方程、露点方程和闪蒸方程。,8,College of Chemistry & Chemical Engineering,SWPU,3.1 Phase Behavior of Natural,3.1 Phase Behavior of Natural Gas,三、烃类体系的相图,烃类体系的相态不仅与体系中烃类物质的组成有关，而且还取决于体系所处的,T,、,p,和所占体积。,1、单组分体系,理论及实验的研究同时表明，纯物质的相态存在下述关系式：,f,(,p,V,T,)0,(1)纯组分,p,-,V,-,T,相图,9,College of Chemistry & Chemical Engineering,SWPU,3.1 Phase Behavior of Natural,3.1 Phase Behavior of Natural Gas,由一系列面组成，表示单相或由两相组成的混合相。,固,固,液,液,液-汽,汽,气,临界点,三相线,固- 汽,气,临界点,液-汽,液,固,固- 汽,三相线,汽,10,College of Chemistry & Chemical Engineering,SWPU,3.1 Phase Behavior of Natural,3.1 Phase Behavior of Natural Gas,将三维立体相图进行的投影可得到二维平面相图。,11,College of Chemistry & Chemical Engineering,SWPU,3.1 Phase Behavior of Natural,3.1 Phase Behavior of Natural Gas,(2)纯组分,p,-,T,相图,C点：Critical Point,表示纯物质汽-液两相可以共存的最高温度,T,C,和最高压力,P,C,。,高于,T,C,和,P,C,，由虚线隔开的区域为,密流,区。,密流区的流体称超临界流体，在这个区域流体的属性不同与气体也不同于液体，它具有特殊的属性。,12,College of Chemistry & Chemical Engineering,SWPU,3.1 Phase Behavior of Natural,3.1 Phase Behavior of Natural Gas,当,p,p,C,，,T,T,C,时：,两相性质相同。,当,T,T,C,，,p,p,C,时,处于密流区，不符合气体定义，也不符合液体定义。,当,p,p,C,，,T,T,C,时,等温加压或等压降温均可液化，属于,汽体,；,当,p,p,C,，,T,T,C,时,等温加压可变为流体，等压降温可液化，属于,气体,；,13,College of Chemistry & Chemical Engineering,SWPU,3.1 Phase Behavior of Natural,3.1 Phase Behavior of Natural Gas,(3)纯组分,p,-,V,相图,在单相区，等温线为光滑的曲线或直线。,高于临界温度的等温线光滑无转折点；,低于临界温度的等温线有转折点，由三部分组成。,等温线在临界点处出现水平拐点，该点一阶导数和二阶导数皆为零。,14,College of Chemistry & Chemical Engineering,SWPU,3.1 Phase Behavior of Natural,3.1 Phase Behavior of Natural Gas,2、多组分体系,与多组分体系的差别：体系中每增加一个组分就多增加一个自由度。,f,(,p,V,T,x,1,x,2, ,x,C,)0,p,V,B,D,两相区,C,汽,液,T,T,C,T,=,T,C,T,T,C,C,临界点,(1)多组分体系,p,-,V,相图,等温相变伴随着压力改变。,凝析过程(,D,B,),压力升高。,蒸发过程(,B,D,),压力降低。,15,College of Chemistry & Chemical Engineering,SWPU,3.1 Phase Behavior of Natural,3.1 Phase Behavior of Natural Gas,(2)多组分体系,p,-,T,相图,泡点线,露点线,液气,临界点,纯组分泡点线和露点线重合为单一蒸气压曲线，而多组分则存在露点线和泡点线；,混合物临界点C，不是气液相能够共存的,最高压力，也不是气-液相能够平衡共存的最高温度。,16,College of Chemistry & Chemical Engineering,SWPU,3.1 Phase Behavior of Natural,2,3.1 Phase Behavior of Natural Gas,3、天然气的反凝析现象,A,B,C,T,m,p,m,温度T,压力,p,气液两相区,液,气,3,1,4,(1)等温反凝析,(2)等温反蒸发,(3)等压反凝析,(4)等压反蒸发,将上述四种反常现象统称作“反凝析现象”。,17,College of Chemistry & Chemical Engineering,SWPU,23.1 Phase Behavior of Natural,3.1 Phase Behavior of Natural Gas,天然气等温反凝析过程分析如下：,设体系原始态为A，等温降压AF,AB(露点)降压,相变：气相开始出现液相；,BD 降压：,相变：BB,1,B,2,B,3,D(反常相变),液相：010203040。,CDC,T,BC为反凝析区,D E (露点) 降压：,相变：DD,3,D,2,D,1,E(正常相变),液相：403020100。,EF降压：单一气相,18,College of Chemistry & Chemical Engineering,SWPU,3.1 Phase Behavior of Natural,3.1 Phase Behavior of Natural Gas,4、天然气反凝析现象的利用,在天然气开采和油气藏类型识别方面具有重要意义。,天然气井的开采多属于等温降压过程，当天然气从高压井中喷出时，由于压力降低，出现反凝析现象，会凝析出液体烃。,油气藏开发等温反凝析凝析气藏；,研究意义：指导凝析气藏开发，减少凝析油，在地层中的损失。,19,College of Chemistry & Chemical Engineering,SWPU,3.1 Phase Behavior of Natural,3.2 Equation of State,状态方程：描述物质存在状态的基本关系式。,状态方程的发展先后有三次具有重要意义的突破：,范德华(Van der Waals)方程,Redlich-Kwong(RK)方程,基于普遍化范德华配分函数发展的扰动硬链理论(PHCT),目前，,开发的状态方程有数百个，但既能用于非极性和极性化合物，又有较高精度且形式简单、计算方便的状态方程则不多见。,20,College of Chemistry & Chemical Engineering,SWPU,3.2 Equation of State 状态方程：描述物,3.2 Equation of State,一、范德华(van der Waals)方程,理想气体状态方程是最为简单和古老的状态方程：,pV,RT,van der Waals于1873年首次提出了用于真实气体的状态方程,式中：,a,和,b,称为范德华常数，,a,/,V,2,为压力修正项，,b,为体积修正项，,a,和,b,与气体的临界常数有关。,21,College of Chemistry & Chemical Engineering,SWPU,3.2 Equation of State 一、范德华(va,3.2 Equation of State,应用：,不能精确描述流体的性质，不适于工程应用，但它在状态方程的发展史上却具有重要意义。,二、Redlich-Kwong(RK)方程,Redlich和Kwong于1949年在范德华方程的基础上进一步修正后得到一个两参数状态方程。,1、方程标准形式,式中：,a,反映分子间吸收力的大小，,b,反映分子的大小。,22,College of Chemistry & Chemical Engineering,SWPU,3.2 Equation of State 应用：不能精确描,3.2 Equation of State,a,、,b,可由临界等温线所具有的特性定出,a,=,a,R,2,T,C,2.5,/,p,C,=0.42748,R,2,T,C,2.5,/,p,C,b,=,b,RT,C,/,p,C,=0.08664,RT,C,/,p,C,2、多项式形式,RK方程通常又通过压缩因子表达如下：,Z,3,Z,2,(,A,B,B,2,),A,AB,0,方程参数：,A,=,ap,/(,R,2,T,2,),B,=,bp,/(,RT,),23,College of Chemistry & Chemical Engineering,SWPU,3.2 Equation of State a、b可由临界等,3.2 Equation of State,RK方程应用于混合物时，参数,a,m,和,b,m,可由纯组分的相应参数,a,i,和,b,i,按某种规则计算。,式中：,k,ij,为二元交互作用参数。通常由二元气-液相平衡实验数据回归得到。,24,College of Chemistry & Chemical Engineering,SWPU,3.2 Equation of State RK方程应用于混,3.2 Equation of State,三、Soave-Redlich-Kwong方程(SRK方程),1972年Soave将偏心因子作为第三个参数引入RK状态方程。,将RK状态方程中,a,/,T,0.5,一项改用具有普遍意义的温度和偏心因子的函数,a,(,T,)来代替。,1、方程标准形式,25,College of Chemistry & Chemical Engineering,SWPU,3.2 Equation of State 三、Soave-,3.2 Equation of State,2、多项式形式,Z,3,Z,2,(,A,B,B,2,),Z,AB,0,方程参数：,A,=,a,(,T,),p,/(,RT,),2,B,=,bp,/(,RT,),3、对纯组分方程参数计算式,26,College of Chemistry & Chemical Engineering,SWPU,3.2 Equation of State 2、多项式形式2,3.2 Equation of State,在SRK方程中,a,(,T,)可表示成,4、对混合物方程参数计算式,在处理混合物时，可采用以下混合规则,27,College of Chemistry & Chemical Engineering,SWPU,3.2 Equation of State 在SRK方程中a,3.2 Equation of State,5、SRK方程的热力学计算式,(1)逸度系数计算式,(2)焓差计算式,(3) 熵差计算式,28,College of Chemistry & Chemical Engineering,SWPU,3.2 Equation of State 5、SRK方程的,3.2 Equation of State,四、Peng-Robinson状态方程(PR方程),Peng Robinson指出，SRK虽取得明显改进，但仍有不足之处，如预测液相密度欠准确：普遍小于实验数据。,Peng和Robinson在Soave模型基础上作了一些改进。,1、方程标准形式,29,College of Chemistry & Chemical Engineering,SWPU,3.2 Equation of State 四、Peng-R,3.2 Equation of State,2、多项式形式,Z,3,(1,B,),Z,2,(A2,B,3,B,2,),Z,(,AB,B,2,B,3,) 0,方程参数：,A,=,a,(,T,),p,/(,RT,),2,B,=,bp,/(,RT,),3、对纯组分方程参数计算式,在PR方程中,a,(,T,)可表示成,30,College of Chemistry & Chemical Engineering,SWPU,3.2 Equation of State,3.2 Equation of State,4、对混合物方程参数计算式,在处理混合物时，可采用以下混合规则,5、PR方程的热力学计算式,(1)逸度系数计算式,31,College of Chemistry & Chemical Engineering,SWPU,3.2 Equation of State 4、对混合物方程,3.2 Equation of State,(2)焓差计算式,(3)熵差计算式,32,College of Chemistry & Chemical Engineering,SWPU,3.2 Equation of State (2)焓差计算式,3.2 Equation of State,六、维里(virial)状态方程,Kamerling Onnes于是1901年提出的virial状态方程可以表示成以下两种形式：,(1)将压缩因子,Z,表示为摩尔体积倒数1/,V,的幂级数,式中：,B,、,C,、,D,、分别表示为第二、第三、第四、维里系数，均与压力和密度无关。,对纯组分仅为温度的函数，对混合物尚与组成有关。,33,College of Chemistry & Chemical Engineering,SWPU,3.2 Equation of State 六、维里(vir,3.2 Equation of State,(2)将压缩因子,Z,表示为压力,p,的幂级数,式中：,B,、,C,、,D,、的意义同前。,可以证明两种表达式系数之间具有如下的关系：,维里方程具有可靠的理论基础，不存在任何假设。维里系数与分子间的作用力直接相关。,34,College of Chemistry & Chemical Engineering,SWPU,3.2 Equation of State (2)将压缩因子,3.2 Equation of State,第二维里系数考虑了两个分子间的相互作用所造成的对理想气体行为的偏差。,第三维里系数则考虑了三个分子间的相互作用所产生的影响。,更高的维里系数具有类似的物理意义。,维里系数的确定以及由维里方程推导的逸度和焓熵计算式请参见有关手册。,35,College of Chemistry & Chemical Engineering,SWPU,3.2 Equation of State 第二维里系数考虑,3.3 Gas-liquid Equilibrium,一、相平衡计算内容,泡点计算：,一定,p,(或,T,)下，已知液相组成,x,，确定泡点,T,b,(或,p,b,)和平衡的汽相组成,y,；,露点计算：,一定,p,(或,T,)下，已知汽相组成,y,，确定露点,T,d,(或,p,d,)和平衡的液相组成,x,；,等温闪蒸：,给定料液的量,F,和组成,z,F,，计算在一定,p,和,T,下闪蒸得到的汽液两相的量(,V,、,L,)和组成(,y,，,x,)。,二、相平衡常数(Phase Equilibrium Constant),36,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium一、相平,3.3 Gas-liquid Equilibrium,1、定义,在一定的,T,和,p,条件下，气液两相处于平衡时，组分,i,在气相中的摩尔分数,y,i,和液相中的摩尔分数,x,i,的比值。,2、公式,3、理想平衡常数,对于理想气体和液体，由拉乌尔定律得：,37,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium1、定义,3.3 Gas-liquid Equilibrium,对于理想气体和液体，由道尔顿定律得,由拉乌尔定律和道尔顿定律有,则理想平衡常数为：,适用范围：低压，且远离临界压力。,38,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium对于理想,3.3 Gas-liquid Equilibrium,三、相平衡模型的构成,1、物料平衡条件,描述气液相组成、物质量及平衡常数之间的关系。,汽,相,y,1,，,y,2,，,，,y,n,T,，,p,，,V,液,相,x,1,，,x,2,，,，,x,n,L,气液量,L,V,F,各组分在气液相中量的衡算,y,i,V,x,i,L,Z,i,F,气液相组成归一化,y,i,1.0 ,x,i,1.0,39,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium三、相平,3.3 Gas-liquid Equilibrium,经处理，可得气、液相组成和物料平衡的基本方程。,平衡组分分配比：,K,i,y,i,/,x,i,物料的量衡算：,y,i,V,x,i,L,Z,i,F,气相组成方程：,(气化率),液相组成方程：,40,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium经处理，,3.3 Gas-liquid Equilibrium,2、热力学平衡方程,关键在于能否准确确定气液两相达到平衡后各组分的平衡常数,K,i,，,K,i,一般是,T,、,p,和组成的函数。,由热力学知识，气相逸度可表示为：,液相逸度可表示为：,于是，相平衡常数,K,可表示为：,41,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium2、热力,3.3 Gas-liquid Equilibrium,热力学平衡条件方程组：,42,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium热力学平,3.3 Gas-liquid Equilibrium,四、泡点(Bubble Point)计算,已知液相组成,x,i,和压力,p,(或温度,T,)，确定液体刚开始沸腾时的温度,T,b,(或压力,p,)以及产生气相的平衡组成,y,i,。,1、泡点计算方程,(1)相平衡关系：,y,i,K,i,x,i,(2)归一化方程：,(3)相平衡常数：,43,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium四、泡点,3.3 Gas-liquid Equilibrium,2、简化计算法,设,K,仅与,T,和,p,有关，可采用手算。以泡点温度为例：,设,T,P,给定,由,K,-图查读,K,i,是,否,调整,T,结束,其中,为试差的允许偏差，手算中一般取,0.001。,T,调整：,若(1-,K,i,x,i,)0，则,K,i,偏大，原设,T,b,0,过大(,K,i,随,T,而,)，重新假设时应降低,T,b,0,；否则应增大,T,b,0,。,44,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium2、简化,3.3 Gas-liquid Equilibrium,例题：,由正丁烷(1)、正戊烷(2)和正己烷(3)组成的混合物加入到压力为0.2MPa的精馏塔的加料板。,混合物的组成为,x,1,=0.15，,x,2,=0.4，,x,3,=0.45，要求呈饱和液体加入塔中。试计算料液的泡点温度是多少？,温度，,K,1,K,2,K,3,50,2.50,0.76,0.28,60,3.10,0.98,0.38,58,3.00,0.96,0.36,45,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium例题：由,3.3 Gas-liquid Equilibrium,解：,假设温度为50，则,(1-,K,i,x,i,)1(2.50.150.760.40.450.28),0.1950,说明假设的温度过小，重新假设温度为60，则,(1-,K,i,x,i,)1(3.10.150.980.40.450.38),-0.0280,说明假设的温度过大，重新假设温度为58，则,(1-,K,i,x,i,)1(3.00.150.960.40.450.36),0.0040,说明假设的温度合适，则泡点温度为58。,46,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium解：假设,3.3 Gas-liquid Equilibrium,3、严格计算,K,i,应是组成、,p,和,T,的函数，无法通过简单的迭代过程完成泡点计算。,泡点计算框图如右图：内循环为,K,循环，外循环为,T,循环。,47,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium3、严格,3.3 Gas-liquid Equilibrium,五、露点(D,ew Point)计算,已知液相组成,y,i,和压力,p,(或温度,T,)，确定气体刚开始结露时的温度,T,d,(或压力,p,)以及产生液相的平衡组成,x,i,。,1、露点计算方程,(1)相平衡关系：,x,i,y,i,/,K,i,(2)归一化方程：,(3)相平衡常数：,48,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium五、露点,3.3 Gas-liquid Equilibrium,2、简化计算法,设,K,仅与,T,和,p,有关，可采用手算。以露点温度为例：,设,T,P,给定,由,K,-图查读,K,i,是,否,调整,T,结束,其中,为试差的允许偏差，手算中一般取,0.001。,T,调整：,若(,y,i,/,K,i,-1.0)0，则,K,i,偏大，原设,T,d,0,过大(,K,i,随,T,而,)，重新假设时应降低,T,d,0,；否则应增大,T,d,0,。,49,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium2、简化,3.3 Gas-liquid Equilibrium,3、严格计算,K,i,应是组成、,p,和,T,的函数，无法通过简单的迭代过程完成露点计算。,泡点计算框图如右图：内循环为,K,循环，外循环为,T,循环。,50,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium3、严格,3.3 Gas-liquid Equilibrium,六、等温闪蒸,起平衡分离作用，是一种单级平衡分离过程。如部分汽化和部分冷凝过程。,1、等温闪蒸过程,V,y,i,F,Z,i,加热器或冷却器,Q,L,x,i,T,p,2、变量数,该过程涉及的基本变量有,p,、,T,、,F,、,V,、,L,、,z,i,、,x,i,、,y,i,，变量数共为3c5(c为进料的总组分数)。,51,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium六、等温,3.3 Gas-liquid Equilibrium,3、基本方程,(1)相平衡方程,(2)物料平衡方程,(3)归一化方程,4、独立变量数和计算任务,52,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium3、基本,3.3 Gas-liquid Equilibrium,基本方程数为2,c,3，其独立变量数,N,f,应为,需指定(,c,2)个独立变量方能求解。指定,F,、,p,、,T,(或,e,)、,z,i,作为独立变量，共,c,2个。,任务：,汽、液相的量,V,、,L,和组成,y,i,、,x,i,。,5、目标方程和迭代方程,等温闪蒸单元过程的目标方程可由物科平衡方程与相平衡方程相结合得出。,(1)目标方程,53,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium基本方程,3.3 Gas-liquid Equilibrium,(2)迭代方程,采用Newton-Raphson算法。,54,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium(2)迭,3.3 Gas-liquid Equilibrium,(3)相平衡常数,(4)平衡判据,55,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium(3)相,3.3 Gas-liquid Equilibrium,6、闪蒸过程判据,作等温闪蒸计算前，需判断物系在给定,T,、,p,下是否处于两相区。,1,T,T,b,，进料处于泡点，,e,0,1,T,T,b,，,e,0,1,T,T,b,，进料为过冷液体,1,T,T,d,，进料处于露点，,e,1,1,T,T,d,，,e,1,1,T,T,d,，进料为热蒸汽,仅,K,i,z,i,和,z,i,/,K,i,均1，才处于两相区(0,e,1)。,56,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium6、闪蒸,3.3 Gas-liquid Equilibrium,7、简化计算,简化计算时假设,K,仅与,T,和,p,有关，即,K,=,f,(,T,p,)。,判断混合物是否处于两相区；,由物料衡算方程和相平衡方程建立目标方程；,采用Newton迭代法求解目标方程，得到汽化率,e,；,根据,e,和,K,i,求气、液相组成,x,i,、,y,i,。,8、严格计算,结合状态方程模型进行,闪蒸计算框图,如下。,57,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium7、简化,3.3 Gas-liquid Equilibrium,七、绝热闪蒸,是节流阀在绝热情况下减压而导致温度变化的过程。,1、绝热闪蒸过程,2、变量数,基本变量：,F,、,V,、,L,、,z,i,、,x,i,、,y,i,、,p,I,、,p,F,、,T,I,、,T,F,、,H,I,M,、,H,F,M,、,H,V,、,H,L,。,变量总数：3c11,F,、 z,i,；,p,1,、,T,1,、,H,I,M,节流阀,Q,0,T,F,p,F,H,F,M,V,、,y,i,、,H,V,L,、,x,i,、,H,L,58,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium七、绝热,3.3 Gas-liquid Equilibrium,3、基本方程,(1)相平衡方程,(2)物料平衡方程,(3)归一化方程,(4)焓计算方程：,H,I,M,、,H,F,M,、,H,V,、,H,L,共4个方程。,59,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium3、基本,3.3 Gas-liquid Equilibrium,(5)焓平衡方程,4、独立变量数和计算任务,方程数共2,c,8个，则独立变量数,N,f,应为：,通常指定以下,c,3个变量作为独立交量：,F,、,z,i,、,p,I,、,T,I,、,p,F,(或,T,F,)。,任务：,节流后的温度,T,F,(或,p,F,)，平衡汽、液相的量,V,(或,e,)、,L,，组成,y,i,、,x,i,和焓,H,V,、,H,L,等性质。,60,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium(5)焓,3.3 Gas-liquid Equilibrium,5、计算方程,与等温闪蒸相比，增加了热平衡方程和若干计算有关物流焓的方程。,(1)热平衡方程,(2)温度调整式(正割法),61,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium5、计算,3.3 Gas-liquid Equilibrium,(3)汽化率目标方程,(4)迭代方程,(5)相平衡常数,62,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium(3)汽,3.3 Gas-liquid Equilibrium,6、严格计算框图,结合状态方程的,绝热闪蒸计算框图,如下图所示。,(1)内层,e,-循环：,在,K,i,为常数的情况下用Newton迭代法求解气化分率,e,。,(2)中层,K,-循环：,根据,e,-循环求出的,y,i,、,x,i,分别计算,f,i,V,、,f,i,L,，判断是否满足|,f,i,L,/,f,i,V,1|10,-4,的要求。,若不足则由,K,i,计算式调整,K,i,。,(3)外层,T,F,-循环：,当,K,-循环收敛后，计算,H,V,、,H,L,，由,e,求,H,F,M,，判断是否满足|(,H,F,M,H,I,M,)/,H,I,M,|10,-4,的要求。,63,College of Chemistry & Chemical Engineering,SWPU,3.3 Gas-liquid Equilibrium6、严格,相平衡计算上机作业,针对天然气组成特性，编写天然气组分气液平衡计算的计算程序。,基本要求：,建立天然气组分物性数据库，能够取物性数据；,输入、输出数据能够保存；,模块化结构，如计算压缩因子、计算逸度系数、计算气化率、计算焓值等。,64,College of Chemistry & Chemical Engineering,SWPU,相平衡计算上机作业针对天然气组成特性，编写天然气组分气液平衡,