1887

Abstract

Summary

The fully implicit method is the most commonly used approach to solve black-oil problems in reservoir simulation. The method requires repeated linearization of large nonlinear systems and produces ill-conditioned linear systems. We present a strategy to reduce computational time that relies on two key ideas: i) a sequential formulation that decouples flow and transport into separate subproblems, and ii) a highly efficient Gauss-Seidel solver for the transport problems. This solver uses intercell fluxes to find all cells that only depend on their upstream neighbors and groups all remaining cells into local clusters of cells that are mutually dependent because of counter-current flow. The single cells and local clusters can then be sorted and solved in sequence, starting from the inflow and moving gradually downstream, since each new cell or local cluster will only depend on upstream neighbors that have already been computed. Altogether, this gives optimal localization and control of the nonlinear solution process.

This method has been successfully applied to real-field problems using the standard first-order finite volume discretization. Here, we extend the idea to first-order dG methods on fully unstructured grids. We also demonstrate proof of concept for the reordering idea by applying it to the full simulation model of the Norne oil field, using a prototype variant of open source OPM Flow simulator.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.201802130
2018-09-03
2024-04-19
Loading full text...

Full text loading...

References

  1. Appleyard, J.R. and Cheshire, I.M.
    [1982] The Cascade Method for Accelerated Convergence in Implicit Simulators. In: European Petroleum Conference, 25–28 October, London, United Kingdom. Society of Petroleum Engineers.
    [Google Scholar]
  2. Bell, J.B., Trangenstein, J.A. and Shubin, G.R.
    [1986] Conservation Laws of Mixed Type Describing Three-Phase Flow in Porous Media. SIAM Journ. Appl. Math., 46(6), 1000–1017.
    [Google Scholar]
  3. Bratvedt, F., Gimse, T. and Tegnander, C.
    [1996] Streamline computations for porous media flow including gravity. Transp. Porous Media, 25(1), 63–78.
    [Google Scholar]
  4. Christie, M.A. and Blunt, M.J.
    [2001] Tenth SPE Comparative Solution Project: A Comparison of Upscaling Techniques. SPE Reservoir Eval. Eng., 4, 308–317. Url: http://www.spe.org/csp/.
    [Google Scholar]
  5. Cockburn, B. and Shu, C.W.
    [1989] TVB Runge-Kutta Local Projection Discontinuous Galerkin Finite Element Method for Conservation Laws II: General. Math. Comp., 52(186), 411–435.
    [Google Scholar]
  6. [1991] The Runge-Kutta local projection P1 - discontinuous-Galerkin finite element method for scalar conservation laws. RAIRO Modél. Math. Anal. Numér., 25(3), 337–361.
    [Google Scholar]
  7. Datta-Gupta, A. and King, M.J.
    [2007] Streamline simulation: theory and practice, SPE Textbook Series, 11. Society of Petroleum Engineers.
    [Google Scholar]
  8. Eikemo, B., Lie, K.A., Dahle, H.K. and Eigestad, G.T.
    [2009] Discontinuous Galerkin methods for transport in advective transport in single-continuum models of fractured media. Adv. Water Resour., 32(4), 493–506. 10.1016/j.advwatres.2008.12.010.
    https://doi.org/10.1016/j.advwatres.2008.12.010 [Google Scholar]
  9. Gries, S., Stüben, K., Brown, G.L., Chen, D. and Collins, D.A.
    [2014] Preconditioning for efficiently applying algebraic multigrid in fully implicit reservoir simulations. SPE J., 19(04), 726–736.
    [Google Scholar]
  10. Jenny, P., Lee, S.H. and Tchelepi, H.A.
    [2006] Adaptive fully implicit multi-scale finite-volume method for multiphase flow and transport in heterogeneous porous media. J. Comput. Phys., 217(2), 627–641.
    [Google Scholar]
  11. Klausen, R.A., Rasmussen, A.F. and Stephansen, A.
    [2012] Velocity interpolation and streamline tracing on irregular geometries. Comput. Geosci., 16, 276–.
    [Google Scholar]
  12. Klemetsdal, ∅., Berge, R., Lie, K.A., Nilsen, H. and Møyner, O.
    [2017] Unstructured gridding and consistent discretizations for reservoirs with faults and complex wells. In: SPE Reservoir Simulation Conference. Society of Petroleum Engineers.
    [Google Scholar]
  13. Klemetsdal, Ø.S., Møyner, O. and Lie, K.A.
    [2018] Robust Nonlinear Newton Solver with Adaptive Interface-Localized Trust Regions. SPE J. Submitted.
    [Google Scholar]
  14. Kwok, F. and Tchelepi, H.
    [2007] Potential-based reduced Newton algorithm for nonlinear multiphase flow in porous media. J. Comput. Phys., 227(1), 706 – 727.
    [Google Scholar]
  15. Lie, K.A.
    [2016] An Introduction to Reservoir Simulation Using MATLAB: User guide for the Matlab Reservoir Simulation Toolbox (MRST). SINTEF Digital, Department of Applied Mathematics.
    [Google Scholar]
  16. Lie, K.A., Krogstad, S., Ligaarden, I.S., Natvig, J.R., Nilsen, H.M. and Skaflestad, B.
    [2012a] Open source MATLAB implementation of consistent discretisations on complex grids. Comput. Geosci., 16, 322–.
    [Google Scholar]
  17. Lie, K.A., Natvig, J.R. and Nilsen, H.M.
    [2012b] Discussion of dynamics and operator splitting techniques for two-phase flow with gravity. Int. J. Numer. Anal. Mod., 9(3), 684–700.
    [Google Scholar]
  18. Lie, K.A., Nilsen, H.M., Rasmussen, A.F. and Raynaud, X.
    [2014] Fast Simulation of Polymer Injection in Heavy-Oil Reservoirs on the Basis of Topological Sorting and Sequential Splitting. SPE J., 19(06), 0991–1004.
    [Google Scholar]
  19. Lie, K.-A.and Møyner, O., Natvig, J.R., Kozlova, A., Bratvedt, K., Watanabe, S. and Li, Z.
    [2017] Successful application of multiscale methods in a real reservoir simulator environment. Comput. Geosci., 21(5), 981–998.
    [Google Scholar]
  20. Møyner, O.
    [2017] Nonlinear solver for three-phase transport problems based on approximate trust regions. Comput. Geosci., 21(5–6), 999–1021.
    [Google Scholar]
  21. Natvig, J.R. and Lie, K.A.
    [2008] Fast computation of multiphase flow in porous media by implicit discontinuous Galerkin schemes with optimal ordering of elements. J. Comput. Phys., 227(24), 10108–10124.
    [Google Scholar]
  22. Natvig, J.R., Lie, K.A., Eikemo, B. and Berre, I.
    [2007] An efficient discontinuous Galerkin method for advective transport in porous media. Adv. Water Resour., 30(12), 2424–2438. 10.1016/j.advwatres.2007.05.015.
    https://doi.org/10.1016/j.advwatres.2007.05.015 [Google Scholar]
  23. OPM
    [2017] The Open Porous Media (OPM) Initative. https://opm-project.org/.
  24. Rasmussen, A.F. and Lie, K.A.
    [2014] Discretization of Flow Diagnostics on Stratigraphic and Unstructured Grids. In: 16th European Conference on the Mathematics of Oil Recovery, Catalania, Sicily, Italy. European Association of Geoscientits and Engineers.
    [Google Scholar]
  25. Shahvali, M. and Tchelepi, H.A.
    [2013] Efficient coupling for non-linear multiphase flow with strong gravity. In: SPE Reservoir Simulation Symposium, 18–20 February, The Woodlands, Texas, USA. Society of Petroleum Engineers.
    [Google Scholar]
  26. Sheth, S.M. and Younis, R.M.
    [2017] Localized Solvers for General Full-Resolution Implicit Reservoir Simulation. In: SPE Reservoir Simulation Conference. Society of Petroleum Engineers.
    [Google Scholar]
  27. Sobolev, S.L. and Vaskevich, V.L.
    [2013] The Theory of Cubature Formulas, 415. Springer Science & Business Media.
    [Google Scholar]
  28. Trangenstein, J.A. and Bell, J.B.
    [1989] Mathematical structure of the black-oil model for petroleum reservoir simulation. SIAM J. Appl. Math., 49(3), 749–783.
    [Google Scholar]
  29. Trottenberg, U., Oosterlee, C.W. and Schuller, A.
    [2000] Multigrid. Academic press.
    [Google Scholar]
  30. Watts, J.
    [1986] A Compositional Formulation of the Pressure and Saturation Equations. SPE Reservoir Engineering, 1(03), 243–252.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.201802130
Loading
/content/papers/10.3997/2214-4609.201802130
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