Robust Nonlinear Newton Solver With Adaptive Interface-Localized Trust Regions
- Øystein S. Klemetsdal (Norwegian University of Science and Technology) | Olav Møyner (Norwegian University of Science and Technology) | Knut-Andreas Lie (SINTEF Digital)
- Document ID
- Society of Petroleum Engineers
- SPE Journal
- Publication Date
- August 2019
- Document Type
- Journal Paper
- 1,576 - 1,594
- 2019.Society of Petroleum Engineers
- Newton's method, flux function, multiphase flow, trust region, transport solver
- 8 in the last 30 days
- 142 since 2007
- Show more detail
- View rights & permissions
The interplay of multiphase-flow effects and pressure/volume/temperature behavior encountered in reservoir simulations often provides strongly coupled nonlinear systems that are challenging to solve numerically. In a sequentially implicit method, many of the essential nonlinearities are associated with the transport equation, and convergence failure for the Newton solver is often caused by steps that pass inflection points and discontinuities in the fractional-flow functions. The industry-standard approach is to heuristically chop timesteps and/or dampen updates suggested by the Newton solver if these exceed a predefined limit. Alternatively, one can use trust regions (TRs) to determine safe updates that stay within regions that have the same curvature for numerical flux. This approach has previously been shown to give unconditional convergence for polymer- and waterflooding problems, also when property curves have kinks or near-discontinuous behavior. Although unconditionally convergent, this method tends to be overly restrictive. Herein, we show how the detection of oscillations in the Newton updates can be used to adaptively switch on and off TRs, resulting in a less-restrictive method better suited for realistic reservoir simulations. We demonstrate the performance of the method for a series of challenging test cases ranging from conceptual 2D setups to realistic (and publicly available) geomodels such as the Norne Field and the recent Olympus model from the Integrated Systems Approach for Petroleum Production (ISAPP) optimization challenge.
|File Size||1 MB||Number of Pages||19|
Christie, M. A. and Blunt, M. J. 2001. Tenth SPE Comparative Solution Project: A Comparison of Upscaling Techniques. SPE Res Eval & Eng 4 (4): 308–317. SPE-72469-PA. https://doi.org/10.2118/72469-PA.
Demidov, D. 2018. AMGCL—A C++ Library for Solving Large Sparse Linear Systems With Algebraic Multigrid Method, March 2018, https://github.com/ddemidov/amgcl (accessed 24 May 2019).
Demidov, D. and Rossi, R. 2017. Subdomain Deflation and Algebraic Multigrid: Combining Multiscale With Multilevel. arXiv preprint arXiv:1710.03940.
Fritsch, F. N. and Carlson, R. E. 1980. Monotone Piecewise Cubic Interpolation. SIAM J. Numer. Anal. 17 (2): 238–246. https://doi.org/10.1137/0717021.
Gries, S., Stu¨ben, K., Brown, G. L. et al. 2014. Preconditioning for Efficiently Applying Algebraic Multigrid in Fully Implicit Reservoir Simulations. SPE J. 19 (4): 726–736. SPE-163608-PA. https://doi.org/10.2118/163608-PA.
Hilden, S. T., Møyner, O., Lie, K.-A. et al. 2016. Multiscale Simulation of Polymer Flooding With Shear Effects. Transp Porous Media 113 (1): 111–135. https://doi.org/10.1007/s11242-016-0682-2.
Integrated Systems Approach for Petroleum Production (ISAPP). 2017. The Olympus Benchmark Case, http://www.isapp2.com/optimization-challenge/reservoir-model-description.html (accessed 24 May 2019).
Jenny, P., Lee, S. H., and Tchelepi, H. A. 2006. Adaptive Fully Implicit Multi-Scale Finite-Volume Method for Multi-Phase Flow and Transport in Heterogeneous Porous Media. J Comput Phys 217 (2): 627–641. https://doi.org/10.1016/j.jcp.2006.01.028.
Jenny, P., Tchelepi, H. A., and Lee, S. H. 2009. Unconditionally Convergent Nonlinear Solver for Hyperbolic Conservation Laws With S-Shaped Flux Functions. J Comput Phys 228 (20): 7497–7512. https://doi.org/10.1016/j.jcp.2009.06.032.
Jiang, J. and Tchelepi, H. A. 2018. Nonlinear Acceleration of Sequential Fully Implicit (SFI) Method for Coupled Flow and Transport in Porous Media. arXiv preprint arXiv:1810.02326.
Killough, J. E. and Wheeler, M. F. 1987. Parallel Iterative Linear Equation Solvers: An Investigation of Domain Decomposition Algorithms for Reservoir Simulation. Presented at the SPE Symposium on Reservoir Simulation, San Antonio, Texas, 1–4 February. SPE-16021-MS. https://doi.org/10.2118/16021-MS.
Klemetsdal, Ø. S., Rasmussen, A. F., Møyner, O. et al. 2018. Nonlinear Gauss-Seidel Solvers With Higher Order For Black-Oil Models. Presented at ECMOR XVI–16th European Conference on the Mathematics of Oil Recovery, Barcelona, Spain, 3–6 September. https://doi.org/10.3997/2214-4609.201802130.
Krogstad, S., Lie, K.-A., Møyner, O. et al. 2015. MRST-AD—An Open-Source Framework for Rapid Prototyping and Evaluation of Reservoir Simulation Problems. Presented at the SPE Reservoir Simulation Symposium, Houston, 23–25 February. SPE-173317-MS. https://doi.org/10.2118/173317-MS.
Krogstad, S., Lie, K.-A., Nilsen, H. M. et al. 2017. Efficient Flow Diagnostics Proxies for Polymer Flooding. Computat Geosci 21 (5): 1203–1218. https://doi.org/10.1007/s10596-017-9681-9.
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. https://doi.org/10.1016/j.jcp.2007.08.012.
Li, B. and Tchelepi, H. A. 2014. Unconditionally Convergent Nonlinear Solver for Multiphase Flow in Porous Media Under Viscous Force, Buoyancy, and Capillarity. Energy Procedia 59: 404–411. https://doi.org/10.1016/j.egypro.2014.10.395.
Lie, K.-A. 2019. An Introduction to Reservoir Simulation Using MATLAB/GNU Octave: User Guide for the MATLAB Reservoir Simulation Toolbox (MRST). Cambridge, UK: Cambridge University Press.
Lie, K.-A., Nilsen, H. M., Rasmussen, A. F. et al. 2014. Fast Simulation of Polymer Injection in Heavy-Oil Reservoirs on the Basis of Topological Sorting and Sequential Splitting. SPE J. 19 (6): 991–1004. SPE-163599-PA. https://doi.org/10.2118/163599-PA.
Lie, K.-A., Møyner, O., Natvig, J. R. et al. 2017. Successful Application of Multiscale Methods in a Real Reservoir Simulator Environment. Computat Geosci 21 (5): 981–998. https://doi.org/10.1007/s10596-017-9627-2.
Lu, B., Alshaalan, T., and Wheeler, M. F. 2007. Iteratively Coupled Reservoir Simulation for Multiphase Flow. Presented at the SPE Annual Technical Conference and Exhibition, Anaheim, California, 11–14 November. SPE-110114-MS . https://doi.org/10.2118/110114-MS.
MATLAB is a registered trademark of The MathWorks, Inc., 1 Apple Hill Drive, Natick, Massachusetts 01760-2098.
Moncorgé, A., Tchelepi, H., and Jenny, P. 2017. Modified Sequential Fully Implicit Scheme for Compositional Flow Simulation. J Comput Phys 337 (15 May): 98–115. https://doi.org/10.1016/j.jcp.2017.02.032.
Moncorgé, A., Tchelepi, H., and Jenny, P. 2018. Sequential Fully Implicit Formulation for Compositional Simulation Using Natural Variables. J Comput Phys 371 (15 October): 690–711. https://doi.org/10.1016/j.jcp.2018.05.048.
Møyner, O. 2017. Nonlinear Solver for Three-Phase Transport Problems Based on Approximate Trust Regions. Computat Geosci 21 (5–6): 999–1021. https://doi.org/10.1007/s10596-017-9660-1.
Møyner, O. and Lie, K.-A. 2016. A Multiscale Restriction-Smoothed Basis Method for Compressible Black-Oil Models. SPE J. 21 (6): 2079–2096. SPE-173265-PA. https://doi.org/10.2118/173265-PA.
Møyner, O. and Moncorgé, A. 2018. Nonlinear Domain Decomposition Scheme For Sequential Fully Implicit Formulation of Compositional Multiphase Flow. Presented at ECMOR XVI–16th European Conference on the Mathematics of Oil Recovery, Barcelona, Spain, 3–6 September. https://doi.org/10.3997/2214-4609.201802128.
Møyner, O. and Tchelepi, H. 2018. A Mass-Conservative Sequential Implicit Multiscale Method for Isothermal Equation-of-State Compositional Problems. SPE J. 23 (6): 2376–2393. SPE-182679-PA. https://doi.org/10.2118/182679-PA.
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. https://doi.org/10.1016/j.jcp.2008.08.024.
Odeh, A. S. 1981. Comparison of Solutions to a Three-Dimensional Black-Oil Reservoir Simulation Problem. J Pet Technol 33 (1): 13–25. SPE-9723-PA. https://doi.org/10.2118/9723-PA.
Open Porous Media Initiative. 2015. The Norne Dataset, https://github.com/OPM/opm-data (accessed 24 May 2019).
Raynaud, X., Lie, K.-A., Nilsen, H. M. et al. 2016. The Single-Cell Transport Problem for Two-Phase Flow With Polymer. Computat Geosci 20 (3): 495–507. https://doi.org/10.1007/s10596-015-9502-y.
Schlumberger. 2013. ECLIPSE: Technical Description 2013.2. Houston: Schlumberger.
Sheth, S. M. and Younis, R. M. 2017. Localized Solvers for General Full-Resolution Implicit Reservoir Simulation. Presented at the SPE Reservoir Simulation Conference, Montgomery, Texas, 20–22 February. SPE-182691-MS. https://doi.org/10.2118/182691-MS.
Tarjan, R. 1972. Depth-First Search and Linear Graph Algorithms. SIAM J. Comput. 1 (2): 146–160. https://doi.org/10.1137/0201010.
Todd, M. R. and Longstaff, W. J. 1972. The Development, Testing, and Application of a Numerical Simulator for Predicting Miscible Flood Performance. J Pet Technol 24 (7): 874–882. SPE-3484-PA. https://doi.org/10.2118/3484-PA.
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. https://doi.org/10.1137/0149044.
Trottenberg, U., Oosterlee, C. W., and Schuller, A. 2000. Multigrid, first edition. Cambridge, Massachusetts: Academic Press.
Voskov, D. V. and Tchelepi, H. 2011. Compositional Nonlinear Solver Based on Trust Regions of the Flux Function Along Key Tie-Lines. Presented at the SPE Reservoir Simulation Symposium, The Woodlands, Texas, 21–23 February. SPE-141743-MS. https://doi.org/10.2118/141743-MS.
Wallis, J. R. 1983. Incomplete Gaussian Elimination as a Preconditioning for Generalized Conjugate Gradient Acceleration. Presented at the SPE Reservoir Simulation Symposium, San Francisco, 15–18 November. SPE-12265-MS. https://doi.org/10.2118/12265-MS.
Wang, X. and Tchelepi, H. A. 2013. Trust-Region Based Solver for Nonlinear Transport in Heterogeneous Porous Media. J Comput Phys 253 (15 November): 114–137. https://doi.org/10.1016/j.jcp.2013.06.041.
Watts, J. 1986. A Compositional Formulation of the Pressure and Saturation Equations. SPE Res Eng 1 (3): 243–252. SPE-12244-PA. https://doi.org/10.2118/12244-PA.
Younis, R., Tchelepi, H. A., and Aziz, K. 2010. Adaptively Localized Continuation-Newton Method—Nonlinear Solvers That Converge All the Time. SPE J. 15 (2): 526–544. SPE-119147-PA. https://doi.org/10.2118/119147-PA.