1887
  1. Home
  2. Conferences
  3. Conference Proceedings
  4. ECMOR XIV - 14th European Conference on the Mathematics of Oil Recovery
  5. Conference paper

Abstract

Summary

The Fully-Implicit method (FIM) is often the method of choice for the temporal discretization of the partial differential equations governing multiphase flow in porous media, especially when nonlinearity is severe. The FIM offers unconditional stability, but requires the solution of large coupled systems of nonlinear algebraic equations. Newton-based methods – often with damping heuristics - are employed to solve the nonlinear systems. However, Newton-based solvers can suffer from convergence problems, especially for large time steps in the presence of highly nonlinear flow physics. To overcome such convergence problems, the timestep is reduced, and the Newton steps are restarted from the solution of the previous (converged) time step. Recently, potential ordering and the reduced-Newton method were used to solve immiscible three-phase flow in the presence of buoyancy and capillary effects (e.g., ). Here, we extend the potential-ordering method to interphase mass transfer. Specifically, we deal with the black-oil model with variable bubble-point pressure. The convergence properties and the computational efficiency of the potential-ordering nonlinear solver are superior to existing damped Newton methods.

The nonlinear iteration process starts with the latest pressure field. Here, we use the Algebraic MultiGrid (AMG) from Fraunhofer (Stueben, International Multigrid Conference, 1983). Based on the latest pressure (potential) distribution, a directed graph is formed, in which nodes represent grid cells (control volumes) and edges represent phase fluxes between cells. As proposed by Natvig & Lie (ECMOR, 2008), and Shahvali & Tchelepi ( ), Tarjan’s strongly connected components algorithm is used to order the nonlinear discrete system into a block triangular form. For the transport step, the potential ordering is used to update the saturations/compositions, one cell (control volume) at a time – from highest to lowest phase potential. The algorithm deals effectively with mass transfer between the liquid and gas phases, including phase disappearance (e.g., gas going back in solution) and reappearance (e.g., gas leaving solution), as a function of pressure and composition.

The new nonlinear ordering-based approach was implemented in Stanford’s general-purpose research simulator (AD-GPRS), and was applied to challenging Black-Oil problems using highly heterogeneous permeability fields (e.g., layers of the SPE10 test case). Detailed robustness and performance comparisons of the potential based solver with state-of-the-art nonlinear/linear solvers (e.g. damped Newton with CPR-based AMG) are presented for variable bubble-point black-oil problems using highly detailed 3D heterogeneous models. The results show that for large time steps (e.g., corresponding to fluid throughput – CFL – numbers on the order of hundreds to thousands), our nonlinear ordering-based solver reduces the number of nonlinear iterations significantly, which also leads to gains in the overall computational cost.

© EAGE Publications BV
Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.20141767
2014-09-08
2024-04-20

  • From This Site

    /content/papers/10.3997/2214-4609.20141767
    dcterms_title,dcterms_subject,pub_keyword
    -contentType:Journal -contentType:Contributor -contentType:Concept -contentType:Institution
    10
    5
Loading full text...

Full text loading...

References

  1. Appleyard, J.R., Cheshire, I.M. et al.
    [1982] The Cascade method for accelerated convergence in implicit simulators. European Petroleum Conference, Society of Petroleum Engineers, doi:10.2118/12804‑MS.
     https://doi.org/10.2118/12804-MS [Google Scholar]
  2. Aziz, K. and Settari, A.
    [1979] Petroleum reservoir simulation, vol. 476. Applied Science PublishersLondon.
     [Google Scholar]
  3. Christie, M., Blunt, M. et al.
    [2001] Tenth SPE comparative solution project: A comparison of upscaling techniques. SPE Reservoir Evaluation & Engineering, 4(04), 308–317, doi:10.2118/72469‑PA.
     https://doi.org/10.2118/72469-PA [Google Scholar]
  4. Coats, K.H. et al.
    [2000] A note on IMPES and some IMPES-based simulation models. SPE Journal, 5(03), 245–251, doi:10.2118/65092‑PA.
     https://doi.org/10.2118/65092-PA [Google Scholar]
  5. Deuflhard, P.
    [2011] Newton methods for nonlinear problems: affine invariance and adaptive algorithms, vol. 35. Springer.
     [Google Scholar]
  6. Iranshahr, A., Voskov, D.V. and Tchelepi, H.A.
    [2013] Tie-simplex based compositional space parameterization: Continuity and generalization to multiphase systems. AIChE Journal, 59(5), 1684–1701.
     [Google Scholar]
  7. Jenny, P., Tchelepi, H.A. and Lee, S.H.
    [2009] Unconditionally convergent nonlinear solver for hyperbolic conservation laws with S-shaped flux functions. Journal of Computational Physics, 228(20), 7497–7512, doi: 10.1016/j.jcp.2009.06.032.
     https://doi.org/10.1016/j.jcp.2009.06.032 [Google Scholar]
  8. Killough, J.
    [1995] Ninth SPE comparative solution project: a reexamination of Black-Oil simulation. 1995 SPE Symposium on Reservoir Simulation, The Woodlands, Texas, USA, Society of Petroleum Engineers, doi:10.2118/29110‑MS.
     https://doi.org/10.2118/29110-MS [Google Scholar]
  9. Kwok, F. and Tchelepi, H.
    [2007] Potential-based reduced-Newton algorithm for nonlinear multiphase flow in porous media. Journal of Computational Physics, 227(1), 706–727, doi:10.1016/j.jcp.2007.08.012.
     https://doi.org/10.1016/j.jcp.2007.08.012 [Google Scholar]
  10. Li, B., Benson, S. and Tchelepi, H.
    [2014] Nonlinearity analysis for efficient modelling of long-term CO2 storage.
     [Google Scholar]
  11. Lie, K.A., Nilsen, H.M., Rasmussen, A.F. and Raynaud, X.
    [2013] Fast simulation of polymer injection in heavy-oil reservoirs based on topological sorting and sequential splitting. 2013 SPE Reservoir Simulation Symposium, The Woodlands, Texas, USA, Society of Petroleum Engineers, doi:10.2118/163599‑MS.
     https://doi.org/10.2118/163599-MS [Google Scholar]
  12. 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. Journal ofComputational Physics, 227(24), 10108– 10124, doi:10.1016/j.jcp.2008.08.024.
     https://doi.org/10.1016/j.jcp.2008.08.024 [Google Scholar]
  13. Ortega, J.M. and Rheinboldt, W.C.
    [2000] Iterative solution of nonlinear equations in several variables, vol. 30. Siam.
     [Google Scholar]
  14. Peaceman, D.W.
    [2000] Fundamentals ofnumerical reservoir simulation. Elsevier.
     [Google Scholar]
  15. Russell, T. et al.
    [1989] Stability analysis and switching criteria for adaptive implicit methods based on the CFL condition. 1989 SPE Symposium on Reservoir Simulation, The Woodlands, Texas, USA, Society of Petroleum Engineers, doi:10.2118/18416‑MS.
     https://doi.org/10.2118/18416-MS [Google Scholar]
  16. Saad, Y. and Schultz, M.H.
    [1986] GMRES: A generalized minimal residual algorithm for solving nonsymmetric linear systems. SIAM Journal on scientific and statistical computing, 7(3), 856–869.
     [Google Scholar]
  17. Shahvali, M. and Tchelepi, H.
    [2013] Efficient coupling for nonlinear multiphase flow with strong gravity. 2013 SPE Reservoir Simulation Symposium, The Woodlands, Texas, USA, Society of Petroleum Engineers, doi:10.2118/163659‑MS.
     https://doi.org/10.2118/163659-MS [Google Scholar]
  18. Stüben, K.
    [2001] A review of algebraic multigrid. Journal of Computational and Applied Mathematics, 128(1), 281–309.
     [Google Scholar]
  19. Tarjan, R.
    [1972] Depth-first search and linear graph algorithms. SIAM journal on computing, 1(2), 146–160.
     [Google Scholar]
  20. Voskov, D.V. and Tchelepi, H.A.
    [2012] Comparison of nonlinear formulations for two-phase multi-component eos based simulation. Journal of Petroleum Science and Engineering, 82, 111–.
     [Google Scholar]
  21. Wallis, J., Kendall, R., Little, T. et al.
    [1985] Constrained residual acceleration of conjugate residual methods. 1985 SPE Reservoir Simulation Symposium, The Woodlands, Texas, USA, Society of Petroleum Engineers, doi:10.2118/13536‑MS.
     https://doi.org/10.2118/13536-MS [Google Scholar]
  22. Wang, X. and Tchelepi, H.A.
    [2013] Trust-region based solver for nonlinear transport in heterogeneous porous media. Journal of Computational Physics, 253, 114–137, doi:10.1016/j.jcp.2013.06.041.
     https://doi.org/10.1016/j.jcp.2013.06.041 [Google Scholar]
  23. Younis, R., Tchelepi, H. and Aziz, K.
    [2010] Adaptively localized continuation-Newton method–nonlinear solvers that converge all the time. SPE Journal, 15(2), 526–544, doi:10.2118/119147‑PA.
     https://doi.org/10.2118/119147-PA [Google Scholar]
  24. Younis, R.M.
    [2011] Modern advances in software and solution algorithms for reservoir simulation. Ph.D. thesis, Stanford University.
     [Google Scholar]
  25. Zhou, Y.
    [2012] Parallel general-purpose reservoir simulation with coupled reservoir models and multisegment wells. Ph.D. thesis, Stanford University.
     [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.20141767
Loading
Ordering-based Nonlinear Solver for Fully-implicit Simulation
Conference Proceedings 2014, 1 (2014); https://doi.org/10.3997/2214-4609.20141767
/content/papers/10.3997/2214-4609.20141767
/content/papers/10.3997/2214-4609.20141767
Loading

Data & Media loading...

Most Cited This Month Most Cited RSS feed

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