×

Nonlinear solver for three-phase transport problems based on approximate trust regions. (English) Zbl 1396.76061

Summary: Implicit transport solvers used in reservoir simulation can take longer time steps than explicit solvers, but for long time steps, the commonly used Newton-Raphson’s method will often fail to converge. The convergence issues may manifest themselves as oscillating residuals even though the implicit discretization itself is stable. This behavior occurs because the fractional flow-type flux functions often change between convex and concave during long time steps, resulting in multiple contraction regions for the Newton-Raphson solver. The common strategy to overcome this is to set limits on the saturation changes during the nonlinear iteration, but such a limit has to be determined on a case by case basis, excess iterations may be required, and practical convergence is not guaranteed for a given problem. Previous work on this problem by multiple authors has resulted in solvers based on trust regions, where unconditional convergence can be obtained for incompressible two-phase flow provided a priori analytical knowledge of the flux function exists. The goal of our work is to extend this methodology to a solver where inflection points demarking the different contraction regions do not need to be explicitly known. Instead, these values are estimated during the solution process, giving improved convergence by a local computation for each interface in the simulation model. By systematically reducing updates over regions known to produce convergence issues, it is possible to greatly reduce the computational expense, making the same formulation suitable for an arbitrary number of components. We present a series of numerical results, including arbitrary time-step lengths for two and three-phase gravity segregation, as well as three-dimensional gas and water injection problems with wells and a mixture of both viscous and gravity-dominated flow regimes. The test cases are a systematic validation on a wide variety of both analytical and tabulated relative permeability curves.

MSC:

76M12 Finite volume methods applied to problems in fluid mechanics
65M08 Finite volume methods for initial value and initial-boundary value problems involving PDEs
76S05 Flows in porous media; filtration; seepage
76Txx Multiphase and multicomponent flows
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Baker, L.: Three-phase relative permeability correlations. In: SPE Enhanced Oil Recovery Symposium. Society of Petroleum Engineers (1988)
[2] Bratvedt, F; Gimse, T; Tegnander, C, Streamline computations for porous media flow including gravity, Transp. Porous Media, 25, 63-78, (1996) · Zbl 0925.76753 · doi:10.1007/BF00141262
[3] Brenier, Y; Jaffré, J, Upstream differencing for multiphase flow in reservoir simulation, SIAM J. Numer. Anal., 28, 685-696, (1991) · Zbl 0735.76071 · doi:10.1137/0728036
[4] Brooks, RH; Corey, AT, Hydraulic properties of porous media and their relation to drainage design, Trans. ASAE, 7, 26-0028, (1964) · doi:10.13031/2013.40684
[5] Christie, M.A., Blunt, M.J.: Tenth SPE comparative solution project: a comparison of upscaling techniques. In: SPE Reservoir Simulation Symposium. Society of Petroleum Engineers (2001)
[6] Fritsch, FN; Carlson, RE, Monotone piecewise cubic interpolation, SIAM J. Numer. Anal., 17, 238-246, (1980) · Zbl 0423.65011 · doi:10.1137/0717021
[7] Gmelig Meyling, RHJ; Guerillot, D (ed.); Guillon, O (ed.), A characteristic finite element method for solving non-linear convection-diffusion equations on locally refined grids, 255-262, (1990), Arles
[8] Gmelig Meyling, RHJ, Numerical methods for solving the nonlinear hyperbolic equations of porous media flow, 503-517, (1991), Lund · Zbl 0767.76050
[9] Gries, S; Stüben, K; Brown, GL; Chen, D; Collins, DA; etal., Preconditioning for efficiently applying algebraic multigrid in fully implicit reservoir simulations, SPE J., 19, 726-736, (2014) · doi:10.2118/163608-PA
[10] Hamon, FP; Tchelepi, HA, Analysis of hybrid upwinding for fully-implicit simulation of three-phase flow with gravity, SIAM J. Numer. Anal., 54, 1682-1712, (2016) · Zbl 1382.65269 · doi:10.1137/15M1020988
[11] IO Center NTNU: The Norne benchmark case url: http://www.ipt.ntnu.no/ norne/wiki/doku.php (2012)
[12] Jenny, P; Lee, SH; Tchelepi, HA, Multi-scale finite-volume method for elliptic problems in subsurface flow simulation, J. Comput. Phys., 187, 47-67, (2003) · Zbl 1047.76538 · doi:10.1016/S0021-9991(03)00075-5
[13] Jenny, P; Tchelepi, HA; Lee, SH, Unconditionally convergent nonlinear solver for hyperbolic conservation laws with S-shaped flux functions, J. Comput. Phys., 228, 7497-7512, (2009) · Zbl 1391.76553 · doi:10.1016/j.jcp.2009.06.032
[14] Killough, J., Wheeler, M., et al.: Parallel iterative linear equation solvers: an investigation of domain decomposition algorithms for reservoir simulation. In: SPE Symposium on Reservoir Simulation. Society of Petroleum Engineers (1987)
[15] Killough, J., et al.: Ninth SPE comparative solution project: a reexamination of black-oil simulation. In: SPE Reservoir Simulation Symposium. Society of Petroleum Engineers (1995)
[16] Kippe, V; Aarnes, JE; Lie, KA, A comparison of multiscale methods for elliptic problems in porous media flow, Comput. Geosci., 12, 377-398, (2008) · Zbl 1259.76047 · doi:10.1007/s10596-007-9074-6
[17] Kwok, F; Tchelepi, H, Potential-based reduced Newton algorithm for nonlinear multiphase flow in porous media, J. Comput. Phys., 227, 706-727, (2007) · Zbl 1388.76258 · doi:10.1016/j.jcp.2007.08.012
[18] Lee, S; Efendiev, Y, C1-continuous relative permeability and hybrid upwind discretization of three phase flow in porous media, Adv. Water Resour., 96, 209-224, (2016) · doi:10.1016/j.advwatres.2016.07.017
[19] Lee, S; Efendiev, Y; Tchelepi, H, Hybrid upwind discretization of nonlinear two-phase flow with gravity, Adv. Water Resour., 82, 27-38, (2015) · doi:10.1016/j.advwatres.2015.04.007
[20] Li, B; Tchelepi, HA, Unconditionally convergent nonlinear solver for multiphase flow in porous media under viscous force, buoyancy, and capillarity, Energy Procedia, 59, 404-411, (2014) · doi:10.1016/j.egypro.2014.10.395
[21] Li, B; Tchelepi, HA, Nonlinear analysis of multiphase transport in porous media in the presence of viscous, buoyancy, and capillary forces, J. Comput. Phys., 297, 104-131, (2015) · Zbl 1349.76820 · doi:10.1016/j.jcp.2015.04.057
[22] Lie, KA; Natvig, JR; Nilsen, HM, Discussion of dynamics and operator splitting techniques for two-phase flow with gravity, Int. J Numer. Anal. Mod. (Special issue in memory of Magne Espedal), 9, 684-700, (2012) · Zbl 1277.76119
[23] Lie, KA; Nilsen, HM; Rasmussen, AF; Raynaud, X; etal., Fast simulation of polymer injection in heavy-oil reservoirs on the basis of topological sorting and sequential splitting, SPE J., 19, 991-1, (2014) · doi:10.2118/163599-PA
[24] Møyner, O; Lie, KA, A multiscale restriction-smoothed basis method for high contrast porous media represented on unstructured grids, J. Comput. Phys., 304, 46-71, (2016) · Zbl 1349.76824 · doi:10.1016/j.jcp.2015.10.010
[25] Natvig, JR; Lie, KA, Fast computation of multiphase flow in porous media by implicit discontinuous Galerkin schemes with optimal ordering of elements, J. Comput. Phys., 227, 10108-10124, (2008) · Zbl 1218.76029 · doi:10.1016/j.jcp.2008.08.024
[26] Odeh, AS, Comparison of solutions to a three-dimensional black-oil reservoir simulation problem, J. Petrol. Tech., 33, 13-25, (1981) · doi:10.2118/9723-PA
[27] Schlumberger: ECLIPSE 2013.2 Technical Description (2013)
[28] Sheth, S; Younis, R, Localized computation of Newton updates in fully-implicit two-phase flow simulation, Procedia Comput. Sci., 80, 1392-1403, (2016) · doi:10.1016/j.procs.2016.05.445
[29] Stone, H, Probability model for estimating three-phase relative permeability, J. Petrol. Tech., 22, 214-218, (1970) · doi:10.2118/2116-PA
[30] Stone, H.: Estimation of three-phase relative permeability and residual oil data. J. Can. Petrol. Technol. 12(4) (1973)
[31] The Open Porous Media Initiative: The Norne dataset url: https://github.com/OPM/opm-data (2015)
[32] Thiele, M.R., Batycky, R.P., Blunt, M.J.: A streamline-based 3D field-scale compositional reservoir simulator. In: SPE Annual Technical Conference and Exhibition. San Antonio (5-8 October 1997), sPE 38889
[33] Trangenstein, JA; Bell, JB, Mathematical structure of the black-oil model for petroleum reservoir simulation, SIAM J. Appl. Math., 49, 749-783, (1989) · Zbl 0669.76125 · doi:10.1137/0149044
[34] Trottenberg, U., Oosterlee, C.W., Schuller, A.: Multigrid. Academic, New York (2000)
[35] Voskov, D.V., Tchelepi, H.A., et al.: Compositional nonlinear solver based on trust regions of the flux function along key tie-lines. In: SPE Reservoir Simulation Symposium. Society of Petroleum Engineers (2011)
[36] Voskov, D., Tchelepi, A.: Tie-simplex-based nonlinear solver for mass-variables compositional formulation. In: ECMOR XIII-13th European Conference on the Mathematics of Oil Recovery (2012)
[37] Wallis, J., et al.: Incomplete Gaussian elimination as a preconditioning for generalized conjugate gradient acceleration. In: SPE Reservoir Simulation Symposium. Society of Petroleum Engineers (1983)
[38] Wang, X; Tchelepi, HA, Trust-region based solver for nonlinear transport in heterogeneous porous media, J. Comput. Phys., 253, 114-137, (2013) · Zbl 1349.76406 · doi:10.1016/j.jcp.2013.06.041
[39] Watts, JW, A compositional formulation of the pressure and saturation equations, SPE J., 1, 243-252, (1986)
[40] Younis, R; Tchelepi, HA; Aziz, K; etal., Adaptively localized continuation-Newton method—nonlinear solvers that converge all the time, SPE J., 15, 526-544, (2010) · doi:10.2118/119147-PA
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.