Multiphase isenthalpic flash: General approach and its adaptation to thermal recovery of heavy oil

Duncan Paterson, Wei Yan*, Michael L. Michelsen, Erling H. Stenby

*Corresponding author for this work

Research output: Contribution to journalJournal articleResearchpeer-review

2 Downloads (Pure)

Abstract

Inthalpic flash is a basic equilibrium calculation in process simulation. The recent interest in isenthalpic multiphase flash is mainly caused by the need for simulating thermal recovery of heavy oil. We present here systematic solutions to multiphase isenthalpic flash with full thermodynamics (such as EoS models) or with correlations for K‐factors, and discuss how to tailor the general methods to systems encountered in thermal recovery of heavy oil. First, for the general situation with full thermodynamics we recommend a solution strategy which uses Newton's method for rapid convergence in the majority of cases and Q‐function maximization to safeguard convergence when Newton's method fails. The solution procedure is a generalization of Michelsen's state function based two‐phase flash to multiple phases. The general solution does not need special considerations for the components in the system and is not limited to the selected thermodynamic models and the number of phases. For thermal recovery processes where gas, oil, and aqueous phases are typically involved, the stability analysis and initialization steps are tailored to improve the efficiency. Second, as it is quite common in thermal reservoir simulators to describe phase equilibrium and heat properties with temperature‐dependent K‐factors and separate correlations for heat capacities, we propose a formulation as an extension of the ideal solution isothermal flash formulation to solve such problems. It uses a Newton–Raphson procedure to converge in the majority of cases and a nested loop procedure with the outer loop for a temperature search as a fallback approach for convergence. If the correlations for K‐factors and for heat capacities are thermodynamically consistent, the outer loop can be treated as a maximization. Finally, we present systematic tests of the proposed algorithms using examples with full thermodynamics or K‐factor based thermodynamics. The algorithms prove robust and efficient even in challenging cases including a narrow‐boiling system, a degenerate system, and a four‐phase system. The additional computational cost relative to the corresponding isothermal flash is modest and would be suitable for the purpose of thermal reservoir simulation.
Original languageEnglish
JournalAIChE Journal
Volume65
Issue number1
Pages (from-to)281-293
Number of pages13
ISSN0001-1541
DOIs
Publication statusPublished - 2019

Keywords

  • Phase equilibrium
  • Oil shale/tar sands
  • Thermodynamics/classical

Cite this

Paterson, Duncan ; Yan, Wei ; Michelsen, Michael L. ; Stenby, Erling H. / Multiphase isenthalpic flash: General approach and its adaptation to thermal recovery of heavy oil. In: AIChE Journal. 2019 ; Vol. 65, No. 1. pp. 281-293.
@article{28766ca0b2414435bc17cfc06d8a3e73,
title = "Multiphase isenthalpic flash: General approach and its adaptation to thermal recovery of heavy oil",
abstract = "Inthalpic flash is a basic equilibrium calculation in process simulation. The recent interest in isenthalpic multiphase flash is mainly caused by the need for simulating thermal recovery of heavy oil. We present here systematic solutions to multiphase isenthalpic flash with full thermodynamics (such as EoS models) or with correlations for K‐factors, and discuss how to tailor the general methods to systems encountered in thermal recovery of heavy oil. First, for the general situation with full thermodynamics we recommend a solution strategy which uses Newton's method for rapid convergence in the majority of cases and Q‐function maximization to safeguard convergence when Newton's method fails. The solution procedure is a generalization of Michelsen's state function based two‐phase flash to multiple phases. The general solution does not need special considerations for the components in the system and is not limited to the selected thermodynamic models and the number of phases. For thermal recovery processes where gas, oil, and aqueous phases are typically involved, the stability analysis and initialization steps are tailored to improve the efficiency. Second, as it is quite common in thermal reservoir simulators to describe phase equilibrium and heat properties with temperature‐dependent K‐factors and separate correlations for heat capacities, we propose a formulation as an extension of the ideal solution isothermal flash formulation to solve such problems. It uses a Newton–Raphson procedure to converge in the majority of cases and a nested loop procedure with the outer loop for a temperature search as a fallback approach for convergence. If the correlations for K‐factors and for heat capacities are thermodynamically consistent, the outer loop can be treated as a maximization. Finally, we present systematic tests of the proposed algorithms using examples with full thermodynamics or K‐factor based thermodynamics. The algorithms prove robust and efficient even in challenging cases including a narrow‐boiling system, a degenerate system, and a four‐phase system. The additional computational cost relative to the corresponding isothermal flash is modest and would be suitable for the purpose of thermal reservoir simulation.",
keywords = "Phase equilibrium, Oil shale/tar sands, Thermodynamics/classical",
author = "Duncan Paterson and Wei Yan and Michelsen, {Michael L.} and Stenby, {Erling H.}",
year = "2019",
doi = "10.1002/aic.16371",
language = "English",
volume = "65",
pages = "281--293",
journal = "A I Ch E Journal",
issn = "0001-1541",
publisher = "JohnWiley & Sons, Inc.",
number = "1",

}

Multiphase isenthalpic flash: General approach and its adaptation to thermal recovery of heavy oil. / Paterson, Duncan; Yan, Wei; Michelsen, Michael L.; Stenby, Erling H.

In: AIChE Journal, Vol. 65, No. 1, 2019, p. 281-293.

Research output: Contribution to journalJournal articleResearchpeer-review

TY - JOUR

T1 - Multiphase isenthalpic flash: General approach and its adaptation to thermal recovery of heavy oil

AU - Paterson, Duncan

AU - Yan, Wei

AU - Michelsen, Michael L.

AU - Stenby, Erling H.

PY - 2019

Y1 - 2019

N2 - Inthalpic flash is a basic equilibrium calculation in process simulation. The recent interest in isenthalpic multiphase flash is mainly caused by the need for simulating thermal recovery of heavy oil. We present here systematic solutions to multiphase isenthalpic flash with full thermodynamics (such as EoS models) or with correlations for K‐factors, and discuss how to tailor the general methods to systems encountered in thermal recovery of heavy oil. First, for the general situation with full thermodynamics we recommend a solution strategy which uses Newton's method for rapid convergence in the majority of cases and Q‐function maximization to safeguard convergence when Newton's method fails. The solution procedure is a generalization of Michelsen's state function based two‐phase flash to multiple phases. The general solution does not need special considerations for the components in the system and is not limited to the selected thermodynamic models and the number of phases. For thermal recovery processes where gas, oil, and aqueous phases are typically involved, the stability analysis and initialization steps are tailored to improve the efficiency. Second, as it is quite common in thermal reservoir simulators to describe phase equilibrium and heat properties with temperature‐dependent K‐factors and separate correlations for heat capacities, we propose a formulation as an extension of the ideal solution isothermal flash formulation to solve such problems. It uses a Newton–Raphson procedure to converge in the majority of cases and a nested loop procedure with the outer loop for a temperature search as a fallback approach for convergence. If the correlations for K‐factors and for heat capacities are thermodynamically consistent, the outer loop can be treated as a maximization. Finally, we present systematic tests of the proposed algorithms using examples with full thermodynamics or K‐factor based thermodynamics. The algorithms prove robust and efficient even in challenging cases including a narrow‐boiling system, a degenerate system, and a four‐phase system. The additional computational cost relative to the corresponding isothermal flash is modest and would be suitable for the purpose of thermal reservoir simulation.

AB - Inthalpic flash is a basic equilibrium calculation in process simulation. The recent interest in isenthalpic multiphase flash is mainly caused by the need for simulating thermal recovery of heavy oil. We present here systematic solutions to multiphase isenthalpic flash with full thermodynamics (such as EoS models) or with correlations for K‐factors, and discuss how to tailor the general methods to systems encountered in thermal recovery of heavy oil. First, for the general situation with full thermodynamics we recommend a solution strategy which uses Newton's method for rapid convergence in the majority of cases and Q‐function maximization to safeguard convergence when Newton's method fails. The solution procedure is a generalization of Michelsen's state function based two‐phase flash to multiple phases. The general solution does not need special considerations for the components in the system and is not limited to the selected thermodynamic models and the number of phases. For thermal recovery processes where gas, oil, and aqueous phases are typically involved, the stability analysis and initialization steps are tailored to improve the efficiency. Second, as it is quite common in thermal reservoir simulators to describe phase equilibrium and heat properties with temperature‐dependent K‐factors and separate correlations for heat capacities, we propose a formulation as an extension of the ideal solution isothermal flash formulation to solve such problems. It uses a Newton–Raphson procedure to converge in the majority of cases and a nested loop procedure with the outer loop for a temperature search as a fallback approach for convergence. If the correlations for K‐factors and for heat capacities are thermodynamically consistent, the outer loop can be treated as a maximization. Finally, we present systematic tests of the proposed algorithms using examples with full thermodynamics or K‐factor based thermodynamics. The algorithms prove robust and efficient even in challenging cases including a narrow‐boiling system, a degenerate system, and a four‐phase system. The additional computational cost relative to the corresponding isothermal flash is modest and would be suitable for the purpose of thermal reservoir simulation.

KW - Phase equilibrium

KW - Oil shale/tar sands

KW - Thermodynamics/classical

U2 - 10.1002/aic.16371

DO - 10.1002/aic.16371

M3 - Journal article

VL - 65

SP - 281

EP - 293

JO - A I Ch E Journal

JF - A I Ch E Journal

SN - 0001-1541

IS - 1

ER -