1887

Abstract

Summary

Phase equilibrium calculation is at the heart of compositional reservoir simulation. The conventional example is the isothermal isobaric flash (T,P) which must be solved in every grid-block at each iteration during simulation. This is conventionally solved through the repeated use of stability analysis and phase-split calculation. To accurately represent the fluid properties it is useful to use an equation of state.

Most compositional reservoir simulations use the cubic equations of state (SRK or PR). More advanced equations of state which, for example, take account of association are attractive alternatives (e.g. CPA or PC-SAFT) for some simulations. Each of these models are functions for the Helmholtz energy with the natural variables (T,V,n), using the conventional flash framework it is necessary to solve the equation for volume (at a given pressure) at each iteration. Though simple for the cubic equations of state this is a more significant issue when using advanced equations of state where the association equations must also be solved iteratively at each volume iteration. This increase in computational cost is one reason that the more advanced equations of state are not yet in common use.

An alternative framework to solve the isothermal isobaric flash problem is possible. Instead of solving the equation of state for volume at each iteration the pressure of the phases is matched only at the final equilibrium point. This allows for the natural variables, (T,V,n), of the equation of state to be used to co-solve the equation of state with the equilibrium equations. Using this framework means that the equation of state does not need to be solved for volume at each iteration. This means that the more complex equations of state are only marginally more computationally expensive than the simple cubics.

In this work we will present a method to solve the (T,P) flash problem using the natural variables of the equation of state, (T,V,n). The resulting framework will be used in a multiphase, compositional, 3D reservoir simulator and demonstrated using a number of examples. The computational cost of the proposed method will be compared with the conventional method for solving the (T,P) flash problem when solving the same simulation problem.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.201802114
2018-09-03
2020-05-30
Loading full text...

Full text loading...

References

  1. Baker, L.E., Pierce, A.C. and Luks, K.D.
    [1982] Gibbs Energy Analysis of Phase Equilibria. Society of Petroleum Engineers Journal, 22(5), 1–12.
    [Google Scholar]
  2. Gill, P.E. and Murray, W.
    [1974] Newton-type methods for unconstrained and linearly constrained optimization. Mathematical Programming, 7(1), 311–350.
    [Google Scholar]
  3. Gross, J. and Sadowski, G.
    [2001] Perturbed-Chain SAFT: An Equation of State Based on a Perturbation Theory for Chain Molecules. Industrial & Engineering Chemistry Research, 40(4), 1244–1260.
    [Google Scholar]
  4. Heidemann, R.A. and Michelsen, M.L.
    [1995] Instability of Successive Substitution. Industrial & Engineering Chemistry Research, 34(3), 958–966.
    [Google Scholar]
  5. Huron, M.J. and Vidal, J.
    [1979] New mixing rules in simple equations of state for representing vapourliquid equilibria of strongly non-ideal mixtures.
    [Google Scholar]
  6. Kontogeorgis, G.M., Voutsas, E.C., Yakoumis, I.V. and Tassios, D.P.
    [1996] An Equation of State for Associating Fluids. Industrlal & Engineering Chemistry Research, 35(11), 4310–4318.
    [Google Scholar]
  7. McDonald, C.M. and Floudas, C.A.
    [1995] Global optimization for the phase stability problem. AIChE Journal, 41(7), 1798–1814.
    [Google Scholar]
  8. Michelsen, M.L.
    [1982a] The isothermal flash problem. Part I. Stability. Fluid Phase Equilibria, 9(1), 1–19.
    [Google Scholar]
  9. [1982b] The isothermal flash problem. Part II. Phase-split calculation. Fluid Phase Equilibria, 9(1), 21–40.
    [Google Scholar]
  10. [1994] Calculation of multiphase equilibrium. Computers and Chemical Engineering, 18(7), 545–550.
    [Google Scholar]
  11. Michelsen, M.L., Mollerup, J. and Breil, M.P.
    [2008] Thermodynamic Models: Fundamental & Computational Aspects.Tie-Line publications, 2nd edn.
    [Google Scholar]
  12. Nagarajan, N.R., Cullick, A.S. and Griewank, A.
    [1991] New strategy for phase equilibrium and critical point calculations by thermodynamic energy analysis. Part I. Stability analysis and flash. Fluid Phase Equilibria, 62(3), 191–210.
    [Google Scholar]
  13. Nichita, D.V.
    [2018a] A volume-based approach to phase equilibrium calculations at pressure and temperature specifications.Fluid Phase Equilibria, 461, 70–83.
    [Google Scholar]
  14. [2018b] Volume-based phase stability testing at pressure and temperature specifications.Fluid Phase Equilibria, 458, 123–141.
    [Google Scholar]
  15. Paterson, D., Michelsen, M.L., Stenby, E.H. and Yan, W.
    [2018a] RAND-Based Formulations for Isothermal Multiphase Flash. SPE Journal, 23(02), 535–549.
    [Google Scholar]
  16. Paterson, D., Michelsen, M.L., Yan, W. and Stenby, E.H.
    [2018b] Extension of modified RAND to multiphase flash specifications based on state functions other than (T,P).Fluid Phase Equilibria, 458, 288–299.
    [Google Scholar]
  17. Péneloux, A., Rauzy, E. and Fréze, R.
    [1982] A consistent correction for Redlich-Kwong-Soave volumes. Fluid Phase Equilibria, 8(1), 7–23.
    [Google Scholar]
  18. Peng, D.Y. and Robinson, D.B.
    [1976] A New Two-Constant Equation of State. Industrial & Engineering Chemistry Fundamentals, 15(1), 59–64.
    [Google Scholar]
  19. Pereira, F.E., Galindo, A., Jackson, G. and Adjiman, C.S.
    [2014] On the impact of using volume as an independent variable for the solution of PâĂŞT fluid-phase equilibrium with equations of state.Computers & Chemical Engineering, 71, 67–76.
    [Google Scholar]
  20. Petitfrere, M. and Nichita, D.V.
    [2014] Robust and efficient Trust-Region based stability analysis and multiphase flash calculations.Fluid Phase Equilibria, 362, 51–68.
    [Google Scholar]
  21. Rasmussen, C., Krejbjerg, K., Michelsen, M. and Bjurstrøm, K.
    [2006] Increasing the Computational Speed of Flash Calculations With Applications for Compositional, Transient Simulations.SPE Reservoir Evaluation & Engineering, 9(October 2003), 5–8.
    [Google Scholar]
  22. Soave, G.
    [1972] Equilibrium constants from a modified Redlich-Kwong equation of state. Chemical Engineering Science, 27(6), 1197–1203.
    [Google Scholar]
  23. Von Solms, N., Michelsen, M.L., Kontogeorgis, G.M. and Solms, N.V.
    [2003] Computational and physical performance of a modified PC-SAFT equation of state for highly asymmetric and associating mixtures.Ind. Eng. Chem. Res. Industrial & engineering chemistry research, 42, 1098–1105.
    [Google Scholar]
  24. Yan, W., Michelsen, M.L. and Stenby, E.H.
    [2012] Calculation of Minimum Miscibility Pressure Using Fast Slimtube Simulation. In: SPE Improved Oil Recovery Symposium. Society of Petroleum Engineers.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.201802114
Loading
/content/papers/10.3997/2214-4609.201802114
Loading

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error