牛顿—拉夫森算法在相平衡计算中的应用研究.pdfVIP

维普资讯 投计 ，翟与畿 币 68 1999年 8月 PETRoLEUM EXPLoRATIoN AND I)EVELoPMENT 26No 4 干 教 牛顿．拉夫森算法在相平衡计算中的应用研究 何 伟 孙 雷 — — — — — 一 ]l 审国石油天然气集团公司石j由勘探开发科学研究院 一『芒 用于混合物计算的压缩因子三次方程 HU 舌 Z 一(1一B )Z +(A 一2B 一3B )Z 一 (A B 一B 一B )=0 随着油气 田勘探开发水平的发展和提高 ，目前找 到的凝析油气藏所 占比重越来越大。凝析油气 田的开 其中 发机制和规律与黑油油藏有着显著区别 ，在开采 中其 A = [口 (T)]／(RT) B = (6 )／RT 烃类体系会发生相态转变和相间传递等复杂物理化学 用于混合物各组分逸度计算方程 变化 ，为了解决这些问题 ，相态研究就成 了必不可少的 研究内容。相平衡计算常用的算法有逐步迭代法 、 fn = (z_1]_In(z—B) Broyden算法、梯度型搜索算法和最小变量牛顿一拉夫 森算法等，其中的牛顿一拉夫森算法作为一种最基本的 A f2 6＼ z+2414B 一 五I 一 ／ _0 算法还未得到广泛应用。本文用牛顿一拉夫森算法对 PR方程和LHSS方程作了相平衡计算，推导了它们的 其 中 Jacobi导数矩阵，通过实例计算绘制了相图，并将牛顿一 = ∑ (。 岛) (1_kij) 拉夫森算法与逐步迭代法和Broyden算法作了对 比和 评价 当计算平衡气液相混合物各组分逸度时，，则分 别表示气液相组成M、z。PR状态方程是 目前国内外 状态方程应用 较为常用的状态方程，它在许多方面都有很多优点，但 它也并非十全十美 ，例如它比较适合于甲烷含量低 、中 热力学仅仅给出了体系达到两相平衡的条件和逸 间烃和重组分含量较多的凝析油气体系，但我国已发 度计算公式 ，没有涉及到系统的PvT关系。在经验基 现和投人开发的凝析油气藏种类繁多，甲烷含量高，因 础上发展起来的状态方程能以可信的精度描述气相的 此为了适应我 国油气藏特点，又提 出了一种新的精度 PvT关系，因而得到广泛应用。虽然 目前发表的状态 更高的LHsS状态方程，下面就此方程作一下说明。 方程有很多种类型，并且各种功能更强、应用范围更广 2 用于多组分混台体系的 LHSS方程结构 的模型还在不断研究，但到 目前为止还没有一种能适