Computation of multiphase equilibrium phenomena with an equation of state

Abstract This paper presents a general approach to multiphase equilibrium computations using an equation of state. Various types of calculations are described: flash calculations, saturation-pressure and saturation-temperature calculations, and calculations of the feed composition corresponding to a given saturation pressure and temperature. A stagewise procedure for flash calculations using the QNSS method allows efficient computation of phase equilibria for systems with any number of phases. The initial guess of the composition of an extra phase which is required to start the iteration is obtained from an analysis of the Gibbs-energy surface. Newton's method is used to construct multiphase boundaries on various types of diagrams (pressure-temperature, pressure-composition and temperature-composition). Multiphase critical points on these boundaries are also estimated by interpolation. The algorithms developed are tested for a typical reservoir-oil-CO 2 mixture which exhibits three-phase liquid 1-liquid 2-vapor (L 1 L 2 V) separation. The predictions obtained using the Peng-Robinson of state are satisfactory.

[1]  Robert A. Heidemann.  Non-uniqueness in phase and reaction equilibrium computations , 1978 .

[2]  M. Michelsen.  The isothermal flash problem. Part I. Stability , 1982 .

[3]  Khalid Aziz,et al.  An accelerated successive substitution algorithm , 1983 .

[4]  D. Peng,et al.  A rigorous method for predicting the critical properties of multicomponent systems from an equation of state , 1977 .

[5]  D. Peng,et al.  A New Two-Constant Equation of State , 1976 .

[6]  Michael L. Michelsen,et al.  Calculation of phase envelopes and critical points for multicomponent mixtures , 1980 .

[7]  R. Heidemann,et al.  Computation of high pressure phase equilibria , 1983 .

[8]  R. Heidemann,et al.  The calculation of critical points , 1980 .

[9]  K. Aziz,et al.  Computation of Multiphase Equilibrium For Compositional Simulation , 1982 .

[10]  M. Michelsen,et al.  Calculation of critical points from cubic two-constant equations of state , 1981 .

[11]  L. Yarborough,et al.  Multiple Phase Behavior in Porous Media During CO2 or Rich-Gas Flooding , 1977 .

[12]  Yau-Kun Li,et al.  The Development of a General Phase Envelope Construction Algorithm for Reservoir Fluid Studies , 1982 .

[13]  L. X. Nghiem,et al.  A Robust Iterative Method for Flash Calculations Using the Soave-Redlich-Kwong or the Peng-Robinson Equation of State , 1983 .

[14]  R. Risnes,et al.  Equilibrium Calculations for Coexisting Liquid Phases , 1984 .

[15]  F. Orr,et al.  Interpretation of Pressure-Composition Phase Diagrams for CO2/Crude-Oil Systems , 1984 .

[16]  K. D. Luks,et al.  Partial miscibility behavior in cryogenic natural gas systems , 1983 .

[17]  L. X. Nghiem,et al.  A New Approach to Quasi-Newton Methods With Application to Compositional Modeling , 1983 .

[18]  Franklin M. Orr,et al.  The Effect of Phase Behavior on CO2-Flood Displacement Efficiency , 1981 .

[19]  H. H. Rachford,et al.  Procedure for Use of Electronic Digital Computers in Calculating Flash Vaporization Hydrocarbon Equilibrium , 1952 .

[20]  F. Orr,et al.  Phase Behavior of CO2 and Crude Oil in Low-Temperature Reservoirs , 1981 .

[21]  M. Michelsen.  The isothermal flash problem. Part II. Phase-split calculation , 1982 .

[22]  L. E. Baker,et al.  Gibbs energy analysis of phase equilibria , 1982 .