A new algorithm for rachford-rice for multiphase compositional simulation

R. Okuno, Russell Taylor Johns, K. Sepehrnoori

Research output: Chapter in Book/Report/Conference proceedingConference contribution

2 Citations (Scopus)

Abstract

Flash calculations for use in compositional simulation are more difficult and time-consuming as the number of equilibrium phases increase beyond two. Because of its complexity many simulators do not even attempt to incorporate three or more hydrocarbon phases even though such cases are important in many low-temperature gas floods, or for high temperatures where hydrocarbons can partition into water. Multiphase flash algorithms typically use successive substitution (SS) followed by Newton's method. For Np-phase flash calculations, (Np-1) Rachford-Rice (RR) equations are solved in every iteration step in SS, and depending on the choice of independent variables, in Newton's method. Solution of RR equations determines both compositions and amounts of phases for a fixed overall composition and set of K-values. A robust algorithm for RR is critical to obtain convergence in multiphase compositional simulation, and has not been satisfactorily developed unlike the traditional two-phase flash. In this paper, we develop a robust and efficient algorithm for RR equations for multiphase compositional simulation. We focus primarily on systems with three equilibrium hydrocarbon phases, although other applications are possible. We derive a function, whose gradient vector consists of RR equations. This correct solution to the RR equations is formulated as a minimization of the non-monotonic convex function using the independent variables of (Np-1) phase mole fractions. The key to obtaining a robust algorithm is that we specify non-negative constraints for the resulting equilibrium phase compositions, which are described by a very small region with no poles. The minimization uses Newton's direction with a line search technique to exhibit super-linear convergence. We show a case in which a previously developed method cannot converge while our algorithm rapidly converges in a few iterations. We implement the algorithm both in a stand-alone flash code and in UTCOMP, a multiphase compositional simulator, and show that the algorithm is guaranteed to converge when a multiphase region exists as indicated by stability analysis.

Original languageEnglish (US)
Title of host publicationSociety of Petroleum Engineers - SPE Eastern Regional/AAPG Eastern Section Joint Meeting 2008
Pages353-367
Number of pages15
StatePublished - 2008
EventSPE Eastern Regional/AAPG Eastern Section Joint Meeting 2008 - Pittsburgh, PA, United States
Duration: Oct 11 2008Oct 15 2008

Other

OtherSPE Eastern Regional/AAPG Eastern Section Joint Meeting 2008
CountryUnited States
CityPittsburgh, PA
Period10/11/0810/15/08

Fingerprint

rice
simulation
Hydrocarbons
phase equilibrium
Newton-Raphson method
hydrocarbon
simulator
substitution
Substitution reactions
Simulators
Chemical analysis
Phase composition
stability analysis
Phase equilibria
Poles
Gases
Temperature
Water
gas
method

All Science Journal Classification (ASJC) codes

  • Geochemistry and Petrology
  • Geotechnical Engineering and Engineering Geology

Cite this

Okuno, R., Johns, R. T., & Sepehrnoori, K. (2008). A new algorithm for rachford-rice for multiphase compositional simulation. In Society of Petroleum Engineers - SPE Eastern Regional/AAPG Eastern Section Joint Meeting 2008 (pp. 353-367)
Okuno, R. ; Johns, Russell Taylor ; Sepehrnoori, K. / A new algorithm for rachford-rice for multiphase compositional simulation. Society of Petroleum Engineers - SPE Eastern Regional/AAPG Eastern Section Joint Meeting 2008. 2008. pp. 353-367
@inproceedings{1f93a24cd208446e93f249af2f76f01c,
title = "A new algorithm for rachford-rice for multiphase compositional simulation",
abstract = "Flash calculations for use in compositional simulation are more difficult and time-consuming as the number of equilibrium phases increase beyond two. Because of its complexity many simulators do not even attempt to incorporate three or more hydrocarbon phases even though such cases are important in many low-temperature gas floods, or for high temperatures where hydrocarbons can partition into water. Multiphase flash algorithms typically use successive substitution (SS) followed by Newton's method. For Np-phase flash calculations, (Np-1) Rachford-Rice (RR) equations are solved in every iteration step in SS, and depending on the choice of independent variables, in Newton's method. Solution of RR equations determines both compositions and amounts of phases for a fixed overall composition and set of K-values. A robust algorithm for RR is critical to obtain convergence in multiphase compositional simulation, and has not been satisfactorily developed unlike the traditional two-phase flash. In this paper, we develop a robust and efficient algorithm for RR equations for multiphase compositional simulation. We focus primarily on systems with three equilibrium hydrocarbon phases, although other applications are possible. We derive a function, whose gradient vector consists of RR equations. This correct solution to the RR equations is formulated as a minimization of the non-monotonic convex function using the independent variables of (Np-1) phase mole fractions. The key to obtaining a robust algorithm is that we specify non-negative constraints for the resulting equilibrium phase compositions, which are described by a very small region with no poles. The minimization uses Newton's direction with a line search technique to exhibit super-linear convergence. We show a case in which a previously developed method cannot converge while our algorithm rapidly converges in a few iterations. We implement the algorithm both in a stand-alone flash code and in UTCOMP, a multiphase compositional simulator, and show that the algorithm is guaranteed to converge when a multiphase region exists as indicated by stability analysis.",
author = "R. Okuno and Johns, {Russell Taylor} and K. Sepehrnoori",
year = "2008",
language = "English (US)",
isbn = "9781605606736",
pages = "353--367",
booktitle = "Society of Petroleum Engineers - SPE Eastern Regional/AAPG Eastern Section Joint Meeting 2008",

}

Okuno, R, Johns, RT & Sepehrnoori, K 2008, A new algorithm for rachford-rice for multiphase compositional simulation. in Society of Petroleum Engineers - SPE Eastern Regional/AAPG Eastern Section Joint Meeting 2008. pp. 353-367, SPE Eastern Regional/AAPG Eastern Section Joint Meeting 2008, Pittsburgh, PA, United States, 10/11/08.

A new algorithm for rachford-rice for multiphase compositional simulation. / Okuno, R.; Johns, Russell Taylor; Sepehrnoori, K.

Society of Petroleum Engineers - SPE Eastern Regional/AAPG Eastern Section Joint Meeting 2008. 2008. p. 353-367.

Research output: Chapter in Book/Report/Conference proceedingConference contribution

TY - GEN

T1 - A new algorithm for rachford-rice for multiphase compositional simulation

AU - Okuno, R.

AU - Johns, Russell Taylor

AU - Sepehrnoori, K.

PY - 2008

Y1 - 2008

N2 - Flash calculations for use in compositional simulation are more difficult and time-consuming as the number of equilibrium phases increase beyond two. Because of its complexity many simulators do not even attempt to incorporate three or more hydrocarbon phases even though such cases are important in many low-temperature gas floods, or for high temperatures where hydrocarbons can partition into water. Multiphase flash algorithms typically use successive substitution (SS) followed by Newton's method. For Np-phase flash calculations, (Np-1) Rachford-Rice (RR) equations are solved in every iteration step in SS, and depending on the choice of independent variables, in Newton's method. Solution of RR equations determines both compositions and amounts of phases for a fixed overall composition and set of K-values. A robust algorithm for RR is critical to obtain convergence in multiphase compositional simulation, and has not been satisfactorily developed unlike the traditional two-phase flash. In this paper, we develop a robust and efficient algorithm for RR equations for multiphase compositional simulation. We focus primarily on systems with three equilibrium hydrocarbon phases, although other applications are possible. We derive a function, whose gradient vector consists of RR equations. This correct solution to the RR equations is formulated as a minimization of the non-monotonic convex function using the independent variables of (Np-1) phase mole fractions. The key to obtaining a robust algorithm is that we specify non-negative constraints for the resulting equilibrium phase compositions, which are described by a very small region with no poles. The minimization uses Newton's direction with a line search technique to exhibit super-linear convergence. We show a case in which a previously developed method cannot converge while our algorithm rapidly converges in a few iterations. We implement the algorithm both in a stand-alone flash code and in UTCOMP, a multiphase compositional simulator, and show that the algorithm is guaranteed to converge when a multiphase region exists as indicated by stability analysis.

AB - Flash calculations for use in compositional simulation are more difficult and time-consuming as the number of equilibrium phases increase beyond two. Because of its complexity many simulators do not even attempt to incorporate three or more hydrocarbon phases even though such cases are important in many low-temperature gas floods, or for high temperatures where hydrocarbons can partition into water. Multiphase flash algorithms typically use successive substitution (SS) followed by Newton's method. For Np-phase flash calculations, (Np-1) Rachford-Rice (RR) equations are solved in every iteration step in SS, and depending on the choice of independent variables, in Newton's method. Solution of RR equations determines both compositions and amounts of phases for a fixed overall composition and set of K-values. A robust algorithm for RR is critical to obtain convergence in multiphase compositional simulation, and has not been satisfactorily developed unlike the traditional two-phase flash. In this paper, we develop a robust and efficient algorithm for RR equations for multiphase compositional simulation. We focus primarily on systems with three equilibrium hydrocarbon phases, although other applications are possible. We derive a function, whose gradient vector consists of RR equations. This correct solution to the RR equations is formulated as a minimization of the non-monotonic convex function using the independent variables of (Np-1) phase mole fractions. The key to obtaining a robust algorithm is that we specify non-negative constraints for the resulting equilibrium phase compositions, which are described by a very small region with no poles. The minimization uses Newton's direction with a line search technique to exhibit super-linear convergence. We show a case in which a previously developed method cannot converge while our algorithm rapidly converges in a few iterations. We implement the algorithm both in a stand-alone flash code and in UTCOMP, a multiphase compositional simulator, and show that the algorithm is guaranteed to converge when a multiphase region exists as indicated by stability analysis.

UR - http://www.scopus.com/inward/record.url?scp=70349468140&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=70349468140&partnerID=8YFLogxK

M3 - Conference contribution

SN - 9781605606736

SP - 353

EP - 367

BT - Society of Petroleum Engineers - SPE Eastern Regional/AAPG Eastern Section Joint Meeting 2008

ER -

Okuno R, Johns RT, Sepehrnoori K. A new algorithm for rachford-rice for multiphase compositional simulation. In Society of Petroleum Engineers - SPE Eastern Regional/AAPG Eastern Section Joint Meeting 2008. 2008. p. 353-367