Gibbs Energy Analysis of Phase Equilibria

1982 ◽  
Vol 22 (05) ◽  
pp. 731-742 ◽  
Author(s):  
L.E. Baker ◽  
A.C. Pierce ◽  
K.D. Luks

Abstract Equations of state are used to predict or to match equilibrium fluid phase behavior for systems as diverse as distillation columns and miscible gas floods of oil reservoirs. The success of such simulations depends on correct predictions of the number and the compositions of phases present at a given temperature, pressure, and overall fluid composition. For example, recent research has shown that three or more phases may exist in equilibrium in CO floods. This paper shows why an equation of state can predict the incorrect number of phases or incorrect phase compositions. The incorrect phase descriptions still satisfy the usual restrictions on equality of chemical potentials of components in each phase and conservation of moles in the system. A new method and its mathematical proof are presented for determining when a phase equilibrium solution is incorrect. Examples of instances where incorrect predictions may be made are described. These include a binary system in which a two-phase solution may be predicted for a single-phase fluid and a multicomponent CO /reservoir oil system in which three or more phases may coexist. Introduction Advances in reservoir oil recovery methods have necessitated advances in methods for prediction of phase equilibria associated with those methods. It was long considered sufficient to approximate the reservoir behavior of oil and gas systems with models in which compositions of the phases in equilibrium were unimportant. In such a model, the amounts and properties of the phases are dependent on pressure and temperature only. Later, experience in production from condensate and volatile oil reservoirs showed that models incorporating compositional effects were required to simulate the phase equilibria adequately. This led to the use of convergence pressure correlations and subsequently to the development of more sophisticated equation of state methods for modeling and predicting phase equilibria. For adequate description of the compositional effects that occur in enhanced oil recovery processes such as CO and rich gas flooding, an equation-of-state approach is a virtual necessity. equilibrium The use of equations of state for phase prediction is not limited to the petroleum industry. Such equations also find wide use in basic chemical and physical research, and in the refining and chemical processing industries. Solution techniques for phase equilibrium problems are varied and depend to some extent on the application and equation of state used however, there are three restrictions that all phase equilibrium solutions must satisfy. First, material balance must be preserved. Second, for phases in equilibrium there must be no driving force to cause a net movement of any component from one phase to any other phase. In thermodynamic parlance, the chemical potentials for each component must be the same in all phases. Third, the system of predicted phases at the equilibrium state must have the lowest possible Gibbs energy at the system temperature and pressure. The requirement that the Gibbs energy of a system. at a given temperature and pressure, must be a minimum is a statement of the second law of thermodynamics, equivalent to the more common version requiring the entropy of an isolated system to be a maximum. The equivalence is demonstrated formally in Ref. 1, for example. If the Gibbs energy of a predicted equilibrium state is greater than that of another state that also satisfies Requirements 1 and 2, the state with the greater Gibbs energy is not thermodynamically stable. Requirements 1 and 2, material balance and equality of chemical potentials, are used commonly as the sole criteria for solution of phase equilibrium problems. SPEJ P. 731^

1932 ◽  
Vol 6 (6) ◽  
pp. 596-604 ◽  
Author(s):  
D. LeB. Cooper ◽  
O. Maass

An equation of state for gases at low densities is developed, using a new function for the change in viscosity with temperature, also developed herein.The gas law equation takes the form[Formula: see text]or V(1 + KT)(PV − RT) = λT − a where a and b are constants corresponding to those of the Van der Waals' equation, and K is a constant derived from the proposed viscosity function which is, for carbon dioxide,[Formula: see text]where K is a constant and η is the viscosity at an absolute temperature T.In the case of carbon dioxide the equation was found to follow density results with an accuracy of from 0.01% to within experimental limits, and the viscosity function was found to agree with Sutherland's (10) results between −78.5 and 20 °C.Comparisons with several other equations of state are made. These show that the new equation is probably more accurate than any other.An expanded form of the new equation, namely:[Formula: see text]permits calculations of the slopes of isothermals for any temperature. Comparisons are made with experimental data.The expanded form of the equation may be solved for K, giving the expression:[Formula: see text]where [Formula: see text] and [Formula: see text] and ξ = Rb0, and since the equation enables the calculation of the molecular radius r, the viscosity may be calculated for any temperature and pressure over which the equation holds.


1978 ◽  
Vol 18 (03) ◽  
pp. 173-182 ◽  
Author(s):  
D.D. Fussell ◽  
J.L. Yanosik

Abstract Phase equilibria equations that incorporate the Redlich-Kwong equation of state are nonlinear and, therefore, must be solved by an iterative method. The method of successive substitutions commonly is used. This method, however, almost always diverges near the critical region for bubble point, dew point, and two-phase calculations. Iterative methods that converge for these calculations are presented. These iterative methods are called presented. These iterative methods are called "minimum variable Newton-Raphson" (MVNR) methods because they try to minimize the number of variable for which simultaneous iteration is required and use the Newton-Raphson method for the correction step. Procedures are given for obtaining starting values for the first iteration and several example problems are discussed. Introduction Reservoir performance predictions for gas condensate and volatile oil reservoirs require a knowledge of the vapor-liquid phase equilbria of the reservoir fluids. A similar knowledge also is required when studying multiple-contact, miscible oil recovery methods that involve injection of hydrocarbons and/or carbon dioxide. Such knowledge is obtained experimentally or calculated from physical properties of the components of the physical properties of the components of the reservoir fluid system. Calculation is desirable because experimental determination is both laborious and expensive. A common basis for calculation of vapor-liquid phase equilibria is the single-stage separation unit. phase equilibria is the single-stage separation unit. This unit represents a PVT cell in which a fluid mixture of known over-all composition is equilibrated at the temperature and pressure of interest. Liquid and vapor compositions and moles of liquid and vapor per mole of fluid mixture are determined. Reliable estimates of other fluid properties (such as phase densities and viscosities) are obtained readily with these properties. The Redlich-Kwong equation of state is used widely in the petroleum industry for phase equilibria calculations. The phase equilibria equations that incorporate this equation of state are nonlinear. As a result, they must be solved by an iterative method. The method of successive substitutions commonly is used. This method, however, almost always diverges for bubble point, dew point, and two-phase calculations near the critical region. This region is extremely important when studying multiple-contact, miscible oil recovery methods involving CO2 or rich-gas injection because the path of the over-all fluid mixture passes through path of the over-all fluid mixture passes through this region. The method of successive substitutions also will diverge for some fluid mixtures near their saturation (bubble point or dew point) pressure at conditions removed from the critical region. This paper presents a reliable iterative sequence that can be used to predict phase equilibria of multiple-contact, miscible oil recovery methods. The method includes sequences for calculation of the saturation pressure and phase equilibria in the two phase region. These MVNR methods rely on minimization of the number of unknowns for which simultaneous iteration is required and use the Newton-Raphson method for the correction step. Minimization is subject to the constraint that all additional unknowns can be calculated by using simple linear equations or, at most, an iteration method applied to one equation in one unknown. MVNR is compared with the method of successive substitutions for a two-phase fluid mixture at various pressures for a fixed temperature. MVNR also is compared with the method of successive substitutions for saturation-envelope calculations near the critical region. DESCRIPTION OF PHYSICAL SYSTEM The single-stage separation unit is the basis for the phase equilibria calculations discussed in this study. This unit represents a PVT cell in which a fluid mixture of known over-all composition is equilibrated at the temperature and pressure of interest. SPEJ P. 173


1981 ◽  
Vol 21 (05) ◽  
pp. 535-550 ◽  
Author(s):  
S.T. Lee ◽  
R.H. Jacoby ◽  
W.H. Chen ◽  
W.E. Culham

Abstract Experimental phase equilibrium data are presented for three reservoir oils at conditions approximating those encountered in in-situ thermal recovery processes. The fluid systems involved consist of three major groups of components: flue gas, water, and crude oil. Data were measured at temperatures from 204.4 to 371.1°C (400 to 700°F) and pressures from 6996.0 to 20785.6 kPa (1,000 to 3,000 psia). Experimental phase equilibrium data were used to develop a correlation of binary interaction coefficients of crude-oil fractions required for the Peng-Robinson equation of state. Phase equilibrium data predicted using the Peng-Robinson equation of state, using our interaction coefficients, are compared with experimental data. Generally, the Peng-Robinson equation of state predictions were in close agreement with the experimental data. Effect of feed gas/oil ratio and water/oil ratio on the equilibrium coefficients was examined through the Peng-Robinson equation of state. A study on the feasibility of representing the crude oil by only two fractions was made also. This study includes a procedure for lumping the crude-oil fractions and examples showing the importance of mixing rules in determining the pseudo critical properties of lumped fractions. Introduction The steady growth of commercial thermal recovery processes1 has created a need for basic data on phase equilibria that involve water and hydrocarbons ranging from methane to high boiling-point fractions. The in-situ thermal recovery processes often are operated at pressures above 6800 kPa (1,000 psia) and temperatures above 200°C (400°F). Experimental data and theoretical correlations on phase equilibria approximating these systems are virtually nonexistent. Early work by White and Brown2 dealt with high boiling-point hydrocarbon phase equilibria. However, the highest pressure studied was 6894.8 kPa (1,000 psia) and the lightest component was pentane. Poettmann and Mayland,3 on the basis of an empirical correlation,4 constructed charts of equilibrium coefficients, or K values, as functions of pressure and temperature for various boiling-point fractions. But the maximum pressure studied was 6894.8 kPa (1,000 psia). Later, Hoffmann et al.5 studied phase behavior of a gas-condensate system with the highest pressure reaching 20 684.3 kPa (3,000 psia) but the highest temperature investigated was only 94.2°C (201°F). In 1963, Grayson and Streed6 reported experimental vapor/liquid equilibrium data for high-temperature and high-pressure hydrocarbon systems. They also extended the Chao-Seader correlation to cover the higher temperature ranges. However, the. major light component in Grayson and Streed's system was hydrogen. Recently, because of the increasing activity in carbon dioxide flooding processes, the phase equilibria of systems involving carbon dioxide and crude oil has received attention. Simon et al.7 studied phase behavior and other properties of carbon-dioxide/reservoir-oil systems. Shelton and Yarborough8 examined phase behavior in porous media during carbon dioxide or rich-gas flooding. No extensive data on equilibrium coefficients were reported in those papers, and the temperature ranges (out of physical reality) were below 93.5°C (200°F). None of these papers surveyed included water as a component.


2020 ◽  
Vol 10 (8) ◽  
pp. 3689-3709
Author(s):  
Shadman H. Khan ◽  
Anupama Kumari ◽  
G. Dixit ◽  
Chandrajit B. Majumder ◽  
Amit Arora

Abstract The present work focuses on developing a framework for accurate prediction of thermodynamic conditions for single-component hydrates, namely CH4, CO2, N2, H2S, and C2H6 (coded in MATLAB). For this purpose, an exhaustive approach is adopted by incorporating eight different equations of states, namely Peng–Robinson, van der Waals, Soave–Redlich–Kwong, Virial, Redlich–Kwong, Tsai-Teja, Patel, and Esmaeilzadeh–Roshanfekr, with the well-known van der Waals–Platteeuw model. Overall, for I–H–V phase region, the Virial and van der Waals equation of state gives the most accurate predictions with minimum AAD%. For Lw–H–V phase region, Peng–Robinson equation of state is found to yield the most accurate predictions with overall AAD of 3.36%. Also, genetic programming algorithm is adopted to develop a generalized correlation. Overall, the correlation yields quick estimation with an average deviation of less than 1%. The accurate estimation yields a minimal AAD of 0.32% for CH4, 1.93% for C2H6, 0.77% for CO2, 0.64% for H2S, and 0.72% for N2. The same correlation can be employed for fitting phase equilibrium data for other hydrates too. The tuning parameter, n, is to be used for fine adjustment to the phase equilibrium data. The findings of this study can help for a better understanding of phase equilibrium and cage occupancy behavior of different gas hydrates. The accuracy in phase equilibria is intimately related to industrial applications such as crude oil transportation, solid separation, and gas storage. To date, no single correlation is available in the literature that can accurately predict phase equilibria for multiple hydrate species. The novelty of the present work lies in both the accuracy and generalizability of the proposed correlation in predicting the phase equilibrium data. The genetic programming generalized correlation is convenient for performing quick equilibrium prediction for industrial applications.


2021 ◽  
Author(s):  
Eleanor Green ◽  
Roger Powell

<p>Phase equilibrium modelling offers a welcome window onto rock-forming processes. It underpins the principles of geothermobarometry, which today is commonly carried out via pseudosection calculations in software such as THERMOCALC and Perple_X. Increasingly, phase equilibrium modelling is combined with complementary approaches such as diffusion or geodynamical calculations, in order to simulate Earth processes.</p><p>However, as anyone with experience of pseudosection calculations will know, it is not always easy to make sense of a rock through phase equilibrium modelling. Problems may relate to: (1) in what way the assumption of thermodynamic equilibrium may, or may not, be applied; (2) uncertainties in compositional analysis; and (3) uncertainties in the composition-dependent equations of state (<em>x</em>-eos). The <em>x</em>-eos are the building blocks of the modelling – one <em>x</em>-eos is needed to represent each of the mineral and fluid phases in the calculation. </p><p>Of the problems listed above, (3) is the most opaque for the user. In this talk I will discuss the uncertainties associated with the <em>x</em>-eos, and the implications of those uncertainties for thermobarometry and the simulation of Earth processes. I will describe two tools, currently in development, for investigating <em>x</em>-eos-derived uncertainty in thermobarometry.</p>


1981 ◽  
Vol 21 (01) ◽  
pp. 89-97 ◽  
Author(s):  
M.P. Leach ◽  
W.F. Yellig

Abstract This paper presents matches, using a fully compositional model, of the performances of seven laboratory CO2 displacements of a 10-component synthetic oil. The criteria for achieving a match of laboratory performance include (1) comparisons of predicted and experimentally determined oil recovery and (2) effluent compositional profiles for each component as functions of hydrocarbon pore volumes (HCPV) of CO2 injected.An equation of state was tuned to predict single-contact (PVT) phase equilibria for CO2/synthetic-oil mixtures. The model incorporates this equation of state to predict the multiple-contact phase equilibria during a CO2 displacement test. Input to the model were independently determined gas/oil relative permeability characteristics and - for each laboratory displacement - injection rate, effluent pressure, pore volume, and temperature.The experimental displacements were conducted in linear Berea core systems using a synthetic (C1-C14) oil at 120 and 150 degrees F. Three displacements at 120 degrees F have been published by Metcalfe and Yarborough. Previously, it was thought that these displacements were conducted at selected pressures so that oil displacements encompassed immiscible, multiple-contact miscible (MCM), and contact miscible mechanisms. However, the model results show that only contact miscible and MCM displacement mechanisms were involved. To confirm the mechanistic understanding at 120 degrees F, three additional laboratory displacements were conducted at 150 degrees F. These encompassed pressures such that the displacement was controlled by an immiscible, an MCM, and a contact miscible mechanism, respectively. The model results at 150 degrees F match the experimental data and confirm the mechanistic understanding. The experimental and numerical results are in agreement with the minimum miscibility pressure theory of Yellig and Metcalfe.The results of this study confirm the importance of experimentally determined effluent compositional profiles and fully compositional models for CO2 mechanism studies. Introduction A mechanistic understanding of oil displacement by CO2 is basic to establishing the CO2 requirements and predicting performance for field projects. The petroleum engineering community relies heavily on fully compositional models and sophisticated laboratory experiments to acquire this mechanistic understanding.Currently, it is not deemed feasible to use fully compositional models to simulate performance of field-wide miscible floods. However, these models should be capable of predicting performances of laboratory-scale displacements. The results of such predictions will identify important variables that control oil recovery and that must be incorporated in mechanistically simpler field performance simulators. Further, confidence in field models will be enhanced greatly by the demonstrated ability to predict laboratory floods.Studies to improve the effectiveness and efficiency of multicomponent compositional simulators have been reported. A cell-to-cell flash model has been used to study mechanisms in rich-gas drives. However, no investigations have been reported previously that demonstrate that a fully compositional model can predict results of rich gas or CO2 laboratory displacements, including prediction of the phases and compositions developed in situ. SPEJ P. 89^


Author(s):  
Boris S. Bokstein ◽  
Mikhail I. Mendelev ◽  
David J. Srolovitz

This chapter addresses the general features of phase equilibria and applies them to single component systems. Before extending our study of phase equilibria to the interesting case of multiphase, multicomponent systems, we examine the special case of single phase, two-component systems (Chapter 3). Phase equilibria in multiphase, multicomponent systems is deferred until Chapter 4. A single substance may exist in different states. For example, H2O can exist as water vapor, liquid water, or any one of several solid phases (ices). Different states can co-exist indefinitely under certain sets of conditions. Under such conditions, the co-existence of these states suggests that they are in equilibrium with respect to one another, that is, phase equilibrium has been established. It is convenient to graphically represent phase equilibria in the form of phase diagrams. An example of such a diagram for a one-component system (with no solid state allotropes) is shown in Fig. 2.1. The AO, OB, and OC lines represent conditions for which two phases are in equilibrium. Since each set of two-phase equilibrium is represented by a one-dimensional surface (i.e. a line), we see that we can vary one parameter (either T or p) without entering a one-phase region of the diagram. For example, if we set the temperature to T1 we can find a saturated vapour pressure p1 such that the liquid and gas co-exist. Three phases simultaneously co-exist at point O, which is called the triple point. Since the three-phase co-existence surface is zero dimensional (i.e. a point), three-phase equilibrium only exists at a specific temperature and pressure, that is, no conditions can be varied. On the other hand, every single-phase region of the diagram is a two-dimensional area and, hence, we can simultaneously, vary two parameters (i.e. both the temperature and pressure) and still remain in the same single-phase region of the diagram. Equations describing the lines of phase equilibria will be derived in Section 2.2, below. Unlike the lines describing the solid–liquid or solid–vapor co-existence, the liquid–vapor co-existence line terminates in a single-phase region of the diagram.


Sign in / Sign up

Export Citation Format

Share Document