天然气工程-烃类流体相态.pptVIP

十五863子课题验收汇报 3、气液平衡计算（天然气地面工程中又称平衡冷凝） 气液平衡计算方程组建立 气液平衡计算满足的条件 气液平衡计算方法 第二十五页，共七十二页。 十五863子课题验收汇报 ①对每个Ki赋初值。如用Willson公式赋初值 ②用牛顿迭代法求解平衡冷凝方程，确定V ③计算组成xi，yi ④用状态方程计算逸度fLi，fVi ⑤计算平衡常数Ki ⑥对每个组分检查平衡条件是否收敛，条件：| fLi- fVi |＜ε 输入P、T、zi 初值Ki 物质平衡计算V L、V、xi、yi 计算Ki 计算ZL、ZV 计算fLi、fVi | fLi- fVi |＜ε? N Y 气液平衡计算步骤： 第二十六页，共七十二页。 十五863子课题验收汇报 定义：平衡冷凝是气体冷却到一定P、T下达到的平衡。 （1）平衡冷凝方程 ①总物料平衡方程 F = V+L （Kmol/h） 液化率：e = L/F 汽化率：1-e = V/F ②组分i的物料平衡方程 ③平衡冷凝方程 其基本计算步骤和气液平衡计算是一样的，因为液化率相当于液相摩尔分数 4、地面工程中的平衡冷凝计算 第二十七页，共七十二页。 十五863子课题验收汇报 该方法为便于用矩阵形式采用牛顿迭代法计算，把物质平衡方程组中的第n+1个方程用压力的形式表示出来： 露点计算： 泡点计算： 等温闪蒸： （四）常用的第二类相平衡计算数学模型 第二十八页，共七十二页。 十五863子课题验收汇报 三、实际气体状态方程和热力学性质 重点和难点 (1)*EOS的压力形式： (2)*EOS的Z因子形式 (3)*用状态方程进行逸度计算 第二十九页，共七十二页。 十五863子课题验收汇报 其最大的成功之处在于: (1)第一次导出了满足临界点条件，且对V是简单的三次方型状态方程状态方程有明确的物理意义 (2)通过与（Andrews）实测CO2体系临界等温线对比，首次用状态方程阐明气液两相相态转变的连续性 (3)提出了两参数对比态原理 (4)建立和发展能同时精确描述平衡气液两相相态行为的状态方程是可能的 存在的问题： (1)仅对理想气体做了简单的修正，引入a和b时忽略了实际分子几何形态和分子不对称性及温度对压力的影响； (2)其理论临界压缩因子为0.375，远大于实测值0.292至0.264，仅适用于简单的球形对称的非吸性分子体系。 第三十页，共七十二页。 十五863子课题验收汇报 目前在状态方程改进，主要发展方向： (1)基于统计热力学正则分配函数理论发展 (2)由统计热力学刚球扰动理论发展 (3)按摩尔密度展开级数并结合统计热力学发展 (4)基于溶液活度理论发展 (5)根据分子热力学建立偏心硬球模型对Van der Waals方程做半经验半理论改进 前三类有较严密的理论基础，但由于结构复杂，实际应用上存在一些困难 活度理论在描述气液平衡方面不能令人满意 第五类方程目前用得最为广泛（主要介绍PR、SRK –EOS） 第三十一页，共七十二页。 十五863子课题验收汇报 （1）压力形式 1961 ＳＲＫ从分子物理学角度、用非球形不对称分子间相互作用位能（引力和斥力）与简单球形对称非极性分子间的位能偏差来改进： 1、SRK—EOS 对含弱极性组分的非烃－烃混合体系计算有改善。 对含H2S的油气体系误差较大。 第三十二页，共七十二页。 十五863子课题验收汇报 Z 的三次方型状态方程 将ＰＶ＝ｚＲＴ的Ｖ代入SRK的压力方程: 式中： A = f（P,T,yi），B = f（P,T,yi） (2) SRK方程的Z因子型式 第三十三页，共七十二页。 十五863子课题验收汇报 （3）混合体系的SRK方程 kij——组分i和j的二元交互作用系数 SRK-EOS逸度系数 第三十四页，共七十二页。 十五863子课题验收汇报 1976年Peng和Robinson等人对VDW、RK、SRK进行了改进，它两个方面进行了思考? (1)VDW方程的斥力项，从简单性和实用性来讲，任然是目前较好的； (2)对引力项和分子密度作出了深入的分析，给出了引力项更好的结构 2、PR-EOS 第三十五页，共七十二页。 十五863子课题验收汇报 （1）压力形式 第三十六页，共七十二页。 十五863子课题验收汇报 (2)PR-EOS的Z因子形式及逸度方程 第三十七页，共七十二页。 十五863子课题验收汇报 已知：Ｐ、Ｔ、ｙi 求：ｚ和体积特性(Vm、ν、ρm、ρg) 对根的讨论： 在气相区存在一个实根 液相区存在一个实根 两相区存在三个实根。对气相组成yi取Vmmax计算zmax，对液相组成xi取Vmmin计算zmin，中间根无意义。 z3型方程的解法：牛顿迭代法 牛顿