采气工程第二章烃类流体相态课件.ppt

1、主要内容一、烃类流体相态特性 二、油气体系气液平衡计算 三、实际气体状态方程和热力学性质 四、重馏分特征化 五、相态法气藏类型判别简介l相态的几个基本概念 l多组分烃类流体相态变化和原理 l气液平衡的物料和热力学平衡的基本概念 l气液平衡计算的基本步骤 l实际气体状态方程的基本结构 l重馏分特征化基本原理本章掌握重点一、烃类流体相态特性（一）基本概念 l 体系 l 相 l 组分 l 组成 l 自由度 l 相平衡 l 化学位 l 逸度 1、体 系井口保护器井口装置阀 座 锥形阀 弹 簧 孔板压座 外 壳 芯 子 孔板 水气2、相l体系内任何均一的部分 体系总是由一定数目的相构成 相有气、液、固三

2、种，可以不连续存在（如：干冰可由固相直接到气相）相与相之间存在界面 任何一相都能与体系的其它均匀部分分开3、组分纯水盐水5、自由度l在不改变平衡体系中原有平衡相数的条件下可独立改变的量 天然气中一般为压力、温度、浓度、组成等 当某一参数改变引起相数的改变，则该参数不是独立变量盐水4、组成6、化学位l控制物质传递的强度位能控制物质传递的强度位能 它实际上是体系某种物质的偏摩尔内能（内能确定了分子运动速度）当有多相存在时，它决定了某种组分在一相到另一相的传递强度 概念抽象、难于具体化和计算。7、相平衡8、逸度校正了的压力，它是人们在实现由热力学变量（化学位）向物理测量变量转换的一种概念工校正了的压

3、力，它是人们在实现由热力学变量（化学位）向物理测量变量转换的一种概念工具。它直接将数学上的抽象量与实际的可以测量的普通强度量相联系起来。具。它直接将数学上的抽象量与实际的可以测量的普通强度量相联系起来。F等价化学位 逸度系数 可由状态方程计算（对理想气体混合物其值位1.0）iioRT n fifio ()fy Piii i盐水相平衡时，某组分在各相中的逸度相等为相平衡的另一种热力学判定条件。平衡常数KiiLiV 盐水(二）单组分和多组分烃类l单一烃类组分相态特征 l多组分烃类体系相态特征 l多组分体系的反转凝析现象 l常见气藏的P-T相图1 1、单一烃类组分相态特性、单一烃类组分相态特性l理论

4、和实验研究同时表明，对于同一纯物质存在以下关系 它表明图解上述函数关系需要三维空间。在实际应用中三维空间难于应用实际经常使用的是：压力-温度投影图（简称PT图）压力-体积投影图（简称PV图）f P V T(,)0BCD液气、液两相区气临界点TTcT=TcP单组分体积VBCD液气、液两相区气临界点TTcT=TcP多组分体系体积V1 1、单一烃类组分相态特性、单一烃类组分相态特性l理论和实验研究同时表明，对于同一纯物质存在以下关系 它表明图解上述函数关系需要三维空间。在实际应用中三维空间难于应用实际经常使用的是：压力-温度投影图（简称PT图）压力-体积投影图（简称PV图）f P V T(,)0在天

5、然气工业中普遍使用的是 P T 图。（单组分 P-T 相图）（临界点：两相共存的最高温度、压力点）P P D D C C(临临界界点点)固固相相 液液相相 气气相相 H H 三三相相点点 T T单组分 P-T 相图l单一组分烃类体系的相态研究已经揭示了相平衡的许多重要特性。l无论体系多么复杂，多相共存、蒸气压、临界点、连续单相区等重要概念将会应用到任一体系。临界点：气液相内涵性质相同的点，这时气液两相共存。临界参数：Pc，Tc，c，Zc，Vc 正常相变：P衡定，T升高，液相全变为气相。T衡定，P下降，液相也会全变为气相。特点：无两相区，是一个急变过程。2、多组分烃类体系相态特征u单一组分与多组

6、分体系的基本区别在于：体系中每增加一个组分就要增加一个自由度。我们必须规定附加的变量来确定体系的状态。对于一个含有N个组分的体系，其描述方程为：由此可见，多组分体系不可能用三维空间来描述一个完整的体系。f P T V xxxn(,)120 PTTmPmC12345两相共存的最高压力泡点线露点线等液量线气藏类型-相图法：干气藏：TReseroir Tm 湿气藏：TReservoir Tm，但地面分离条件落在两相区 凝析气藏：TcTReservoirTm，Tc临界凝析温度 Pc临界凝析压力3、多组分体系的反转凝析现象 在纯组分体系中，当等温压缩或等压降温时，都会导致有液体的凝析产生。在多组分体系中

7、则常常出现和这相反的逆变现象，即在等温降压或等压升温的情况下，反而会引起液体的凝析。天然气体系中凝析气体系特有的现象2、多组分烃类体系P-T相图及特征两相共存的最高温度3、常见气藏的P-T相图 P 地层条件 分离器 T 干气气藏P-T相图示意图C湿气藏P-T相图示意图 P 地层条件 分离器 T CsepPHB4 P-T图05101520253035100200300400500600T(K)P(MPa)典型特点：干气在任何环节不产生液烃 湿气在分离器有液烃产生 凝析气在地下、井筒和分离器均有液烃产生二、气液平衡及相图计算(一)相平衡计算的基本概念及应用(二)油气体系气液平衡计算物料平衡方程(三

8、)油气体系气液平衡时热力学平衡方程(四)相平衡计算的常用数学模型(五)相图计算1、相平衡 体系中各组分在相与相之间的迁移传递速度相等时的状态 对任一组分，气相分子进入液相的速度和液相分子进入气相的速度相等。特点：P、T一定，平衡时各相组成不变,是一个动态平衡。液/xi汽/yi2、相平衡常数(平衡分配比)定义：平衡时，组分i在气相和液相中的摩尔分数之比。公式：计算方法：收敛压力法 状态方程法iiiiz,T,PfxyK3、相平衡计算模型的原理 相平衡 相平衡时，某组分在各相中的化学位相等 组分i的气相逸度应与液相逸度相等 fvi=fLo,i=1,2,n 逸度单位：与压力单位相同 用途：判断气液两相

9、是否达到平衡的可测工具或称参数。F等价化学位 逸度系数 可由状态方程计算iioRT n fifio ()fy Piii iQuestion：How to calculate composition of every component in system?（一）气液平衡计算基本概念及应用(二)油气体系气液平衡计算物料平衡方程（1）假设条件 体系处于气液相平衡体系处于气液相平衡 等组成膨胀等组成膨胀，孔隙定容，孔隙定容体系温度保持恒定体系温度保持恒定 不考虑油气接触介质的影响不考虑油气接触介质的影响 液/xi汽/yi（2）数学模型（物料平衡方程）设对于一混合烃类体系，取设对于一混合烃类体系，取1

10、mol的物质作为分析单元，那么当其处于任一气液两相的物质作为分析单元，那么当其处于任一气液两相平衡状态对应有以下特性（数学描述）：平衡状态对应有以下特性（数学描述）：平衡气液相的摩尔分量平衡气液相的摩尔分量V和和L分别在分别在0和和1之间变化，且满足归一化条件：之间变化，且满足归一化条件：V+L=1.0 平衡气、液相的组成应分别满足归一化条件：yxii 11或表示为：()yxii 0平衡气、液相中的各组分的量之和应分别等于各自总的量的大小：y Vx Lziii 平衡气、液相中任一组分的分配比例用平衡常数描述：Kyxiii/烃类混合体系中i组分的摩尔组成(二)油气体系气液平衡计算物料平衡方程目标

11、：求在给定压力和温度条件下得气、液两相的摩尔分量及组成 解法：物料平衡方程是高度非线性方程，一般采用迭代法求解（如：牛顿迭代法）。如何利用物料平衡方程求解相平衡问题？VkkzVfiii)1(1)1()(令2)1(1)1()()(VkkzdVVdfVfiii)()(1VfVfVVKK|)(|VfV的初值：进料为液相，可使初值小于0.5；进料为气相，可使初值大于0.5 收敛判据：物料平衡方程的牛顿迭代法要求得平衡气液相的摩尔分量和组成，关键的参数是气液两相的平衡常数Ki，只要能准确知道气液平衡状态下各组成的平衡常数，就能进行气液平衡计算。那么气液两相平衡时其热力学条件和平衡常数有什么关系呢？1、气

12、液平衡常数与气液相逸度的关系 2、逸度系数与逸度、压力和组成的关系 3、气液相平衡热力学方程组 我们已经知道气液两相平衡时，热力学条件是各组分在气相和液相中的逸度相等（ffigil）。已知气液相的逸度表达式分别为：fyPfxPigi igili il 从而：kyxiiiilig 它表示气液平衡时，各组分在气相和液相中的平衡分配比为其逸度系数之比，这即是热力学平衡理论求解相平衡的出发点。逸度系数与逸度、压力和组成的关系为：lnln()lnln()igigiililify Pfx P 只要应用能同时描述气液两相状态方程计算气液相的逸度，就能对相平衡进行精确的描述；相平衡计算的热力学平衡条件方程组为

13、：FP T xyffFP T xyffFP T xyfflglgnnnnlng1111122222000(,)(,)(,)液/xi汽/yiln()lnggjiggViigV TzfpRTRTdVRTZy pzV（三）油气体系气液平衡时热力学平衡方程l油气体系相平衡的统一数学模型：油气体系相平衡的统一数学模型：物料平衡方程组物料平衡方程组 热力学平衡条件热力学平衡条件其中的露点、泡点和闪蒸模型均其中的露点、泡点和闪蒸模型均可从该模型导出可从该模型导出 通过该模型的变化。可计算体系的露通过该模型的变化。可计算体系的露点线、泡点线、等液量线点线、泡点线、等液量线 kTiirii exp.()(/)P

14、r5 37 111 kyxiiiilig yzkkVii ii 11()xzkViii 11()y VxLziii ()()()yxz kkViiiii 1110 FP T xyffFP T xyffFP T xyfflglgnnnnlng1111122222000(,)(,)(,)液/xi汽/yi解释（四）、相平衡计算的常用数学模型1、露点计算 当体系由气相到液相的相态变化时，如果处于露点（有雾状液存在时），其相平衡的物质量的关系表现为平衡气相的摩尔分量V1；平衡液相摩尔分量L0，由平衡气液相的物料平衡方程可得：yzkkVzxzkVzkiiiiiiiiii 1111()()而露点状态液相的组

15、成归一化条件为：xzkiii 1数数学学方方程程：FPTxyffFPTxyffFPTxyffFPTxyzklglgnnnnlngniiii111112222210001(,)(,)(,)(,)对Pd赋初值 对Ki赋初值 计算组成xi 计算逸度fLi，fVi 计算平衡常数Ki 平衡条件检查，即：|fLi-fVi|检查露点是否收敛汽/yi解释露点计算步骤：2、泡点计算 L 1，V 0：yz kkVz kxzkVziiiiiiiiii 1111()()于是泡点状态气相的组成归一化条件为：yz kiii 1泡点定义：第一个气泡刚好析出时的状态。计算泡点就是计算泡点压力或泡点温度。结合气液平衡的热力学条

16、件，得到泡点计算数学方程：FPTxyffFPTxyffFPTxyffFPTxyzklglgnnnn ln gniiii111112222210001(,)(,)(,)(,)给定T（或P）求P（或T）、气相组成和露点计算类似。液/xi汽/yi计算步骤和露点计算相似3、气液平衡计算（天然气地面工程中又称平衡冷凝）u气液平衡计算方程组建立 u气液平衡计算满足的条件 u气液平衡计算方法 F,zi L,xi V,yi P,.T 对每个Ki赋初值。如用Willson公式赋初值 用牛顿迭代法求解平衡冷凝方程，确定V 计算组成xi，yi 用状态方程计算逸度fLi，fVi 计算平衡常数Ki 对每个组分检查平衡条

17、件是否收敛，条件：|fLi-fVi|输入P、T、zi 初值Ki物质平衡计算V L、V、xi、yi计算Ki计算ZL、ZV计算fLi、fVi|fLi-fVi|?NY气液平衡计算步骤：定义：平衡冷凝是气体冷却到一定P、T下达到的平衡。（1）平衡冷凝方程 总物料平衡方程 F=V+L （Kmol/h）液化率：e=L/F 汽化率：1-e=V/F 组分i的物料平衡方程iiiVyLxFziiiye1exz平衡冷凝方程1Ke1ezxiii1Ke1ezKyiiii01e1k1kzT,P,eFiii 其基本计算步骤和气液平衡计算是一样的，因为液化率相当于液相摩尔分数4、地面工程中的平衡冷凝计算u该方法为便于用矩阵形

18、式采用牛顿迭代法计算，把物质平衡方程组中该方法为便于用矩阵形式采用牛顿迭代法计算，把物质平衡方程组中的第的第n+1个方程用压力的形式表示出来：个方程用压力的形式表示出来：u露点计算：露点计算：u u u泡点计算：泡点计算：u u u等温闪蒸：等温闪蒸：11(,)0nigniiiilfFp T xyp11(,)0nilniiiigfFp T xyp11(/1)(,)01(/1)niiligniiiiligzFp T x yV111112222211(,)0(,)0.(,)0(/1)(,)01(/1)lglgnnnn ln gnii li gniiii li gFPTxyffFPTxyffFPTx

19、yffzFpTxyV（四）常用的第二类相平衡计算数学模型三、实际气体状态方程和热力学性质重点和难点(1)*EOS的压力形式：(2)*EOS的Z因子形式(3)*用状态方程进行逸度计算ln()lnggjiggViigV TzfpRTRTdVRTZy pzV其最大的成功之处在于其最大的成功之处在于:(1)第一次导出了满足临界点条件，且对V是简单的三次方型状态方程状态方程有明确的物理意义(2)通过与（Andrews）实测CO2体系临界等温线对比，首次用状态方程阐明气液两相相态转变的连续性(3)提出了两参数对比态原理(4)建立和发展能同时精确描述平衡气液两相相态行为的状态方程是可能的 存在的问题：(1)

20、仅对理想气体做了简单的修正，引入a和b时忽略了实际分子几何形态和分子不对称性及温度对压力的影响；(2)其理论临界压缩因子为0.375，远大于实测值0.292至0.264，仅适用于简单的球形对称的非吸性分子体系。目前在状态方程改进，主要发展方向：(1)基于统计热力学正则分配函数理论发展(2)由统计热力学刚球扰动理论发展(3)按摩尔密度展开级数并结合统计热力学发展(4)基于溶液活度理论发展(5)根据分子热力学建立偏心硬球模型对Van der Waals方程做半经验半理论改进l前三类有较严密的理论基础，但由于结构复杂，实际应用上存在一些困难 l活度理论在描述气液平衡方面不能令人满意 l第五类方程目前

21、用得最为广泛（主要介绍PR、SRK EOS）（1）压力形式 1961 从分子物理学角度、用非球形不对称分子间相互作用位能（引力和斥力）与简单球形对称非极性分子间的位能偏差来改进：1、SRKEOS()R TaTPVbVVbPRTvbav()2c2c2aPTRaccbPRTb0.42747a0.08664b25.0rT1m12176.0574.148.0m 对含弱极性组分的非烃烃混合体系计算有改善。对含H2S的油气体系误差较大。uZ 的三次方型状态方程 将的代入SRK的压力方程:0ABBBAzzz2232RTPaARTbPB 式中：A=f（P,T,yi），B=f（P,T,yi）c2c2aPTRac

22、cbPRTb0.42747a0.08664b25.0rT1m12176.0574.148.0m(2)SRK方程的Z因子型式u（3）混合体系的SRK方程bVVabVRTPmmm ijij5.0jijijimk1aaxxaiiimbxbkij组分i和j的二元交互作用系数zB1bbBABzlnb1zblnmiimii ijj5.0jijijjik1aaxx ijj5.0jijijik1aaxSRK-EOS逸度系数3220mmmmmmmmzzzABBA B1976年Peng和Robinson等人对VDW、RK、SRK进行了改进，它两个方面进行了思考?(1)VDW方程的斥力项，从简单性和实用性来讲，任然

23、是目前较好的；(2)对引力项和分子密度作出了深入的分析，给出了引力项更好的结构2、PR-EOSPaTgvba ttra c t()(,)并并指指出出适适当当地地选选择择g v b(,)的的函函数数形形式式，可可以以更更好好地地反反映映包包括括偏偏心心硬硬球球分分子子体体系系在在内内的的分分子子密密度度对对引引力力项项的的影影响响，并并使使之之更更适适用用于于临临界界区区的的计计算算。P PR R-E EO OS S 给给出出的的g v b(,)的的具具体体结结构构为为：g v bv vbb vb(,)()()（1）压力形式PRTvbaTv vbb vb ()()()PRTv baTv v bb

24、 v bmmmm ()()()aRTPic ic i 04 5 7 2 409 8 71 04222.(.)bRTPic ic i 00 7 7 8 009 8 71 04.(.)iir imT ().11052miii 03 7 4 6 415 4 2 202 6 9 9 2.混合规则：aTxxaakmijijiji jjnin()()()111bxbmiiin 1(2)PR-EOS的Z因子形式及逸度方程P P R R-E E O O S S的的逸逸度度方方程程为为：nfy PbbZnZBABabbnZBZBiiimmmmmmimimmmmm()()()()(.)122224 1 404 1

25、 4其其中中：（注注意意在在平平衡衡计计算算时时yi分分别别表表示示气气液液相相的的组组成成）iiijiji jjnya ak ()().0511l已知：、i l求：和体积特性(Vm、m、g)l对根的讨论：在气相区存在一个实根 液相区存在一个实根 两相区存在三个实根。对气相组成yi取Vmmax计算zmax，对液相组成xi取Vmmin计算zmin，中间根无意义。lz3型方程的解法：牛顿迭代法 牛顿迭代(如SRK）初值怎样取？对气相yizg：z01;对液相xizl：z000AB)BBA(zzz)z(F223)z(F)z(Fzz3、三次方型状态方程求解四、重馏分特征化问题的提出色谱组成准确分析：N2

26、，CO2，H2S，C1，C2，C3,iC4,nC4,iC5,nC5,C6色谱组成不能准确分析：C7，C8，,C44+处理方法 拟组分化：把不易准确测定的第n个组分以后的所有组分合并成一个假的“单一组分”，称为拟组分拟组分 引出新的问题如何正确确定拟组分的Tc,Pc,等相平衡计算所需的热力学性质=引出重馏分特征化重馏分特征化处理技术（1）经验关联式方法经验关联式方法：根据实验测定的热力学性质数据建立经验关联式，用所建立的经验关联式预测重馏分的热力学性质（2）等效碳原子数关联法等效碳原子数关联法：把烃类同系物的碳原子数在数值上看成是连续分布的，将 Cn+的热力学参数与之等效处理，即把 Cn+等效为

27、相同碳原子数的烃类同系物进行其热力学参数的预测 四、重馏分特征化（3）连续热力学分布函数方法连续热力学分布函数方法：通过高精度的色谱分析测定有代表性的油气烃类体系的组成，用概率分布模型描述 Cn+中各个单碳数组分的组成分布规律，建立其热力学性质预测模型（4）以相态实验数据为目标的最优化拟合方法最优化拟合方法：通过对油气体系相态实验数据的拟合计算，获得与体系相匹配的Cn+热力学参数 1、经验关联式法以Edmister的经验关联式为例 式中：Pa-大气压；Pc-临界压力，MPa。Tb-正常沸点；Tc-临界温度，K11)/lg(73)8.0)(360()8.0)(35.225.353()107227

28、.1(03.129.445.05.04206.09371.04bcnCcnaCcncncncnCcncnCcncncnbcncncncnTTPPMPMTMTM2、基于连续热力学理论的等效碳原子数法等效碳原子数概念等效碳原子数概念 将Cn+重馏分的正常沸点，分子量，比重和临界性质等热力学参数与烃类体系同系物的参数作比较，可得到一个整数或非整数碳原子数的“烃类同系物”与之等效。我们把这一碳原子数称为Cn+的等效碳原子数。等效碳原子数关联方法的思路等效碳原子数关联方法的思路 根据等效碳原子数的概念，先将正构烷烃同系物的热力学参数关联为其碳原子数的经验公式，再设法求出Cn+重馏分的等效碳原子数，然后借

29、用正构烷烃同系物的经验公式即可求得其热力学参数 碳原子数Cn5 10 15 20 25 30 35 40 45560 480 400 320 240 160 8044.291.03cncncnMCn=(M+4.028)/13.997式中：Tc-临界温度，K0 5 10 15 205.0 4.03.02.01.00.0Pc碳原子数Cn320000774.00372.0680.0949.5nnnCCCPc式中：Pc-临界压力，MPa等效碳原子数关联方法等效碳原子数关联方法4332107498.32355.09834.5625.85428.162nnnnCCCCTc0.80.60.40.200 5

30、10 15 20碳原子数Cn20.039790.062140.00087nnCC 2、基于连续热力学理论的等效碳原子数法3、基于连续热力学理论的分布函数方法基本思路 准确测定有代表性的黑油，挥发油和凝析气中 Cn+重馏分中各个单碳数组分的组成含量=组分的延伸分析（C7C45+）把分析结果绘制成分子量（或碳原子数）与组成关系的概率分布曲线，建立相应的概率分布数学模型用组成与分子量之间的概率分布模型预测一般油气体系Cn+重馏分中单碳数组分的延伸及其热力学性质 Cn+中各单碳数组分M与zi关系的分布曲线 1.00.10.010.001100 200 300 400 500M2112108922007

31、cMMiMi+1zi/zn+3、基于连续热力学理论的分布函数方法Cn+中各单碳数组分M与zi关系的概率密度分布曲线 0.010.0100 200 300 400 500M2121108922007cM0.008 0.006 0.0040.002p(M)zi=zc7+AREA P(M)-PROBABILITY DENSITY FUNCTION3、基于连续热力学理论的分布函数方法重馏分中各延伸组分应满足物质平衡关系重馏分中各延伸组分应满足物质平衡关系组成分布平衡：各延伸组分组成含量总和等于重馏分的组成组成分布平衡 分子量分布平衡 密度分布平衡NCCnnnMMCCxMdxM)(1 ()nnNnCCC

32、CMMMxxp M式中重馏分的摩尔组成重馏分中最小碳数组分的分子量重馏分中最大碳数组分的分子量重馏分的摩尔组成概率密度分布函数3、基于连续热力学理论的分布函数方法分子量分布平衡：各延伸组分分子量总和等于重馏分的分子量NCCnnnMMCCMMdxx)(1重馏分的摩尔组成的分子量重馏分中最大碳数组分的分子量重馏分中最小碳数组分重馏分的摩尔组成式中nNnnCCCCMMMx 3、基于连续热力学理论的分布函数方法密度分布平衡：各延伸组分密度总和等于重馏分的密度 NCCnNCCnnMMMMCMdxMdx1)()(的分子量 重馏分中最大碳数组分 的分子量 重馏分中最小碳数组分 重馏分的摩尔组成 式中 N n

33、 n C C C M M 分布模型：Mzi关系满足三参数伽玛概率密度分布模型 式中：p(M)-组成的概率密度分布函数 M-分子量 -Cn+重馏分中最小碳数组分的分子量-确定分布曲线形态的特征参数3、基于连续热力学理论的分布函数方法8888880.010.0100 200 300 400 500M0.0080.0060.0040.002p(M)zi=zc7+AREA之间的关系=式中：Cn-Cn+重馏分中最小起始碳原子数 Mcn+-Cn+重馏分的平均分子量，由实验测定给出 形态特征参数分布函数法Cn+重馏分单碳数组分延伸的表征 累积概率分布函数xdxxpxXP)()(累积概率分布函数的解析式xyj

34、yexXPjjy )1(/)(0其中3、基于连续热力学理论的分布函数方法各延伸组分组成分布1)(iinMMCidMMpzz各延伸组分组成分布的解析式01)1(1jjiyjiyCijyeyezziin式中：i=n,n+1,N(=45)咖吗函数的计算)()1()11(1)()1()(1limnkknnnkn !咖吗函数的递推关系00)()()1(jjjj!结束递推计算的精度设置 1610jj3、基于连续热力学理论的分布函数方法重馏分特征化处理过程对Cn+中各单碳数组分进行延伸处理（Spliting)对各单碳数组分的热力学性质进行关联计算对延伸后的单碳数组分进行拟组分化劈分处理(Lumping)（1

35、 1）Cn+Cn+重馏分的延伸重馏分的延伸 目的：根据概率密度函数分布规律，模拟色谱 分析过程将 Cn+中各单碳数组分从起始碳数开 始延伸到某一高碳数。一般是延伸到C45+。意义：用于预测各个单碳数组分的mole组成和 热力学性质。3、基于连续热力学理论的分布函数方法（1）Cn+重馏分的延伸 方法:根据概率密度函数分布规律，以实测的Cn+重馏分平均分子量Mcn+，平均比重 rcn+和组成 zcn+为目标函数进行拟合计算，以获得各个单碳数组分的组成含量和分子量步 骤A:设两个相邻单碳数组分的分子量分别为Mi和Mi+1,则分子量介于Mi和Mi+1之间的组分的mole含量为：1)(MiMicnidM

36、MPzz式中：i=n,n+1,N,表示碳原子数。其选取视精度要求而定，一般最小碳数n可取为7或11；最高碳数N可取为45。3、基于连续热力学理论的分布函数方法B 延伸后的各单碳数的组成需满足以下目标函数式中：Mi-第i个单碳数组分的分子量 ri-第i个单碳数组分的比重 C 根据Watson特征常数的定义和Riazi-Daubert理论，可得以下关系式：0cncncnNniiiiMzMz（1）Cn+重馏分的延伸 0164.11962.241825.115178.0384573.015178.0106607.1)5579.4()(5556.05579.4ibiicniiicnbicncncnTMK

37、MKTMK 为了求解前述目标函数，需要确定的参数有n,N,Mcn+,cn+,Zcn+,以及Mi,Tbi和i三组数据中的任何一组。其中 Mcn+,cn+,Zcn+,由实验测定，若以Mi为调节变量，则可取正构烷烃的分子量为初值。实际计算中，可采用正割法对Mi,Tbi或ri进行迭代，而使610cncncnNniiiiMzMz的解就可以求出各延伸的单碳数组分的zi/zcn+,Mi,ri和Tbi等参数。（2）各延伸的单碳数组分热力学性质的计算选择相适应的经验关联式选择相适应的经验关联式 与各单碳数组分之间的二元交互作用系数)为 式中：C(0668.0 14.0 1 /)/lg(7 3)472(10 71

38、589.1)472(10 53028.5 0629.19 1 1 1 2448.4 86618.3 11 3596.0 3125.2 6 3596.0 5885.0 cici c c i ci c bi ci a ci i bi i bi ci bi i bi ci i bi ci k k k T T P P K T T P K T T P T T-=-=w （3）各延伸组分的拟组分化处理目目 的的 把延伸的单碳数组分按一定规则划分合并成几个拟组分，以便在保证计算精度的情况下减少相平衡计算的工作量。方方 法法 选择一组拟组分化的劈分规则拟组分个数的确定：)lg(3.31 nNINTm式中：n，

39、N分别为单碳数组分的最小和最大碳数各拟组分划分的分界分子量各拟组分划分的分界分子量),.,2,1()ln(1expmjMMmMMjnNnj用用KayKay规则进行各拟组分热力学参数的计算规则进行各拟组分热力学参数的计算 ijijcijicjcijicjijijijijzzPzzPTzzTwwMzzM)/()/()/(/)/(/1)/(iiwzjj w;z 式中：(4)连续热力学分布函数方法综合评述（1）分布模型的合理性:咖玛函数概率分布模型假定分子量与摩尔分数之间存在连续关系，因此不是反映色谱组成分析结果的真实物理模型，但是在数值上合理（2）拟组分数目的合理性:对于一给定的油气体系，需要划分多

40、少个拟组分才能精确预测其相平衡特性？一般情况下，拟组分数目越多，计算结果应越精确。但并不尽然，因为影响计算精度的除了拟组分的个数外，还有拟组分划分方案的合理性。重点是拟组分划分方案的合理性，这需选择合适的混合规则。通常划分的拟组分数为3个。（3）混合规则的选择:Curtis H.Whitson 在“Characterizing HydrocarbonPlus Fractions”(Soc.Pet.Eng.J.(Aug.1983)683-694)一文中同时给出了Kay规则法和平均沸点法两种混合规则。经对比计算，表明两种混合规则的计算结果相差不大，但Kay规则法较为简单，故常用。（4）基础数据的选

41、择:在进行重馏分单碳数组分的延伸时，可以选择单碳数组分的分子量Mi,正常沸点 Tbi，比重 ri中的任一组数据为基础，各自的延伸结果不同，但对油气体系相平衡计算结果影响不大。4、基于相态实验数据的最优化拟合方法问题的提出 应用状态方程方法进行油气藏流体相态模拟计算，必须能明确给出模拟所需的各烃类组分的计算参数如分子量，临界温度，临界压力，偏心因子，以及状态方程的可调参数 ，二元交互作用系数kij等。根据前面的分析，一般对于低碳原子数组分和非烃类气体，都有准确数据可查并可存储在计算程序中供计算选择；但对于高碳数（如Cn+重馏分）组分，虽经过重馏分特征化处理后可得到初始数据，但直接使用初始数据计算

42、可能会产生较大 的误差。因此在进行完整的相平衡计算之前，需通过对相态实验数据的拟合来调整重馏分的热力学参数，以便使计算结果更好地符合实验数据而达最佳匹配。ba,拟合计算一般采用两种方法。一是自动拟合法，即以实测相态数据为目标函数，采用非线性回归方法如单纯形最优化拟合等方法进行热力学参数的匹配；二是采用经验法，即凭人工对各种热力学参数对计算过程的敏感性的认识，确定参数的调整范围，通过多次试凑计算，以获得与实验数据相匹配的结果。两种方法各有其优劣，实用时常先采用人工法进行快速试凑以逼近收敛域，然后通过非线性回归计算得到精确值。非线性回归目标函数的建立0)(),.,(2121eieicimiinYY

43、YWXXXF式中：F(X1,X2,Xn)-目标函数；X1,X2,Xn-所选择的拟合变量(Tcn+，Pcn+,等）；Wi-第i个拟合变量的权重；Yci,Yei-分别为计算和实测的相态参数(如露点压力，泡点压力，反凝析油饱和度，采收率等)。cn4、基于相态实验数据的最优化拟合方法拟合变量的选择：体系中重馏分的临界温度，临界压力，偏心因子和重馏分与C1H4的二元交互作用系数。状态方程中可调的 等参数dcba,拟合目标变量（实测相态数据）的选择:实测的油气体系露点压力，泡点压力；等组成膨胀，定容衰竭，地面分离器闪蒸等数据。五、气藏类型相态判别法简介0.01.02.03.04.05.06.07.08.0

44、9.01 0.0-1 0 0-8 0-6 0-4 0-2 002 04 06 0温 度，压 力，MPaT f=9 0.8 P f=2 0.0 3 M P aP m=9.2 5 M P aT c=-4 1.2 P c=6.3 7 M P aC0.0 2%0.0 5%0.1 0%0.2 0%T m=4 6.0 0.3 0%皮家PK3井相图（Psept1.62 Mpa，T sept32）0.01.02.03.04.05.06.07.08.09.01 0.01 1.01 2.01 3.01 4.0-1 5 0-1 2 0-9 0-6 0-3 003 06 09 01 2 0温 度，压 力，MPaTf=

45、1 0 4.0 Pf=1 7.1 0 M P aP m=1 3.3 2 M P aT c=-8 2.9 P c=3.9 7 M P aC0.0 1%0.0 3%0.0 5%0.1 0%T m=1 0 0.3 03691 21 51 82 1-1 2 0-1 0 0-8 0-6 0-4 0-2 002 04 06 08 01 0 01 2 0温 度，压 力，MPaTf=8 4.2 Pf=2 0.2 9 M P aP m=1 7.1 2 M P aT c=-8 2.0 P c=1 1.4 M P aC0.0 5%0.2%0.4%0.6%Tm=8 2.9 皮家SN139井相图（Psept2.00 M

46、Pa，T sept30.0）井号 参数名 PK3SN92SN1391 182.60263.5196.15油气藏类型无油环凝析气藏各井1参数值1=（C2/C3）+（C1+C2+C3+C4）/C5+气油比（m3/m3）油罐油密度（g/cm3）油气藏类型0.966重质油藏351250.8250.966普通黑油油藏1253500.8020.825黑油油藏与挥发性油藏过渡带3506250.7600.802挥发性油藏62514250.7600.802挥发性油藏与凝析气藏过渡带142512467106860.739低含凝析油气藏湿气气藏按气油比和油罐油密度判别油气藏类型气油比（m3/m3）油罐油密度（g/c

47、m3）油气藏类型0.966重质油藏351250.8250.966普通黑油油藏1253500.8020.825黑油油藏与挥发性油藏过渡带3506250.7600.802挥发性油藏62514250.7600.802挥发性油藏与凝析气藏过渡带142512467106860.739低含凝析油气藏湿气气藏按气油比和油罐油密度判别油气藏类型 利用气层流体组成的四个参数即C2+，C2/C3，100C2/（C3+C4），100（C2+/C1）的经验数据来判别气藏类型气藏类型气藏类型 参数名参数名 气藏气藏无油环无油环 凝析气藏凝析气藏 带油环带油环 凝析气藏凝析气藏 油藏油藏C2+%0.15.05151030

48、2027C2/C3 41602.26.0130.51.3100.C2/(C3+C4)30010000以上以上1704005020020100100(C2+/C1)0.15.05.015104030600四参数判别法参数范围井号 参数名 PK3SN92SN139PN6C2（）5.013.345.111.91C2/C3 4.764.654.274.09100C2/(C3+C4)333.01305.48284.76324.91100（C2/C1）5.373.505.4613.72气藏类型气藏气藏无油环凝析气藏无油环凝析气藏皮家各井四参数值 根据102个油气藏统计（符合率为86%）C5+含量1.75%，无油环凝析气藏 C5+含量1.75%，带油环凝析气藏 C1/C5+52，带油环凝析气藏 C1/C5+52，无油环凝析气藏井号 参数名 PK3SN92SN139PK6C5（）0.550.381.060.31C1/C5 169.73251.0888.07279.94气藏类型无油环凝析气藏