A new algorithm for rachford-rice for multiphase compositional simulation

R. Okuno, Russell Taylor Johns, K. Sepehrnoori

Research output: Contribution to journalArticle

49 Citations (Scopus)

Abstract

Flash calculations for use in compositional simulation are more difficult and time-consuming as the number of equilibrium phases increases 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 gasfloods 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 an algorithm for RR equations for multiphase compositional simulation that is guaranteed to converge to the correct solution independent of the number of phases for both positive and negative flash calculations. 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 nonmonotonic convex function using the independent variables of (Np-1) phase mole fractions. The key to obtaining a robust algorithm is that we specify nonnegative 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 superlinear 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 standalone flash code and in UTCOMP (Chang et al. 1990), 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)
Pages (from-to)313-325
Number of pages13
JournalSPE Journal
Volume15
Issue number2
DOIs
StatePublished - Jan 1 2010

Fingerprint

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

All Science Journal Classification (ASJC) codes

  • Energy Engineering and Power Technology
  • Geotechnical Engineering and Engineering Geology

Cite this

Okuno, R. ; Johns, Russell Taylor ; Sepehrnoori, K. / A new algorithm for rachford-rice for multiphase compositional simulation. In: SPE Journal. 2010 ; Vol. 15, No. 2. pp. 313-325.
@article{effbc27a5f124a1688f5af38f47462ec,
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 increases 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 gasfloods 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 an algorithm for RR equations for multiphase compositional simulation that is guaranteed to converge to the correct solution independent of the number of phases for both positive and negative flash calculations. 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 nonmonotonic convex function using the independent variables of (Np-1) phase mole fractions. The key to obtaining a robust algorithm is that we specify nonnegative 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 superlinear 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 standalone flash code and in UTCOMP (Chang et al. 1990), 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 = "2010",
month = "1",
day = "1",
doi = "10.2118/117752-PA",
language = "English (US)",
volume = "15",
pages = "313--325",
journal = "SPE Journal",
issn = "1086-055X",
publisher = "Society of Petroleum Engineers (SPE)",
number = "2",

}

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

In: SPE Journal, Vol. 15, No. 2, 01.01.2010, p. 313-325.

Research output: Contribution to journalArticle

TY - JOUR

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

AU - Okuno, R.

AU - Johns, Russell Taylor

AU - Sepehrnoori, K.

PY - 2010/1/1

Y1 - 2010/1/1

N2 - Flash calculations for use in compositional simulation are more difficult and time-consuming as the number of equilibrium phases increases 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 gasfloods 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 an algorithm for RR equations for multiphase compositional simulation that is guaranteed to converge to the correct solution independent of the number of phases for both positive and negative flash calculations. 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 nonmonotonic convex function using the independent variables of (Np-1) phase mole fractions. The key to obtaining a robust algorithm is that we specify nonnegative 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 superlinear 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 standalone flash code and in UTCOMP (Chang et al. 1990), 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 increases 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 gasfloods 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 an algorithm for RR equations for multiphase compositional simulation that is guaranteed to converge to the correct solution independent of the number of phases for both positive and negative flash calculations. 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 nonmonotonic convex function using the independent variables of (Np-1) phase mole fractions. The key to obtaining a robust algorithm is that we specify nonnegative 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 superlinear 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 standalone flash code and in UTCOMP (Chang et al. 1990), 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=77954841836&partnerID=8YFLogxK

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

U2 - 10.2118/117752-PA

DO - 10.2118/117752-PA

M3 - Article

AN - SCOPUS:77954841836

VL - 15

SP - 313

EP - 325

JO - SPE Journal

JF - SPE Journal

SN - 1086-055X

IS - 2

ER -