Nonlinear acceleration of sequential fully implicit (SFI) method for coupled flow and transport in porous media

Abstract The sequential fully implicit (SFI) method was introduced along with the development of the multiscale finite volume (MSFV) framework, and has received considerable attention in recent years. Each time step for SFI consists of an outer loop to solve the coupled system, in which there is one inner Newton loop to implicitly solve the pressure equation and another loop to implicitly solve the transport equations. Limited research has been conducted that deals with the outer coupling level to investigate the convergence performance. In this paper we extend the basic SFI method with several nonlinear acceleration techniques for improving the outer-loop convergence. Specifically, we consider numerical relaxation, quasi-Newton (QN) and Anderson acceleration (AA) methods. The acceleration techniques are adapted and studied for the first time within the context of SFI for coupled flow and transport in porous media. We reveal that the iterative form of SFI is equivalent to a nonlinear block Gauss–Seidel (BGS) process. The effectiveness of the acceleration techniques is demonstrated using several challenging examples. The results show that the basic SFI method is quite inefficient, suffering from slow convergence or even convergence failure. In order to better understand the behaviors of SFI, we carry out detailed analysis on the coupling mechanisms between the sub-problems. Compared with the basic SFI method, superior convergence performance is achieved by the acceleration techniques, which can resolve the convergence difficulties associated with various types of coupling effects. We show across a wide range of flow conditions that the acceleration techniques can stabilize the iterative process, and largely reduce the outer iteration count.

[1]  Hamdi A. Tchelepi,et al.  Compositional Multiscale Finite-Volume Formulation , 2014 .

[2]  H. L. Stone Estimation of Three-Phase Relative Permeability And Residual Oil Data , 1973 .

[3]  W. Wall,et al.  Fixed-point fluid–structure interaction solvers with dynamic relaxation , 2008 .

[4]  J. W. Watts,et al.  A Compositional Formulation of the Pressure and Saturation Equations , 1986 .

[5]  Yalchin Efendiev,et al.  Hybrid upwind discretization of nonlinear two-phase flow with gravity , 2015 .

[6]  Yousef Saad,et al.  Two classes of multisecant methods for nonlinear acceleration , 2009, Numer. Linear Algebra Appl..

[7]  Homer F. Walker,et al.  Anderson acceleration and application to the three-temperature energy equations , 2017, J. Comput. Phys..

[8]  Patrick Jenny,et al.  Multiscale finite-volume method for compressible multiphase flow in porous media , 2006, J. Comput. Phys..

[9]  Shingo Watanabe,et al.  A Stable Multi-phase Nonlinear Transport Solver with Hybrid Upwind Discretization in Multiscale Reservoir Simulator , 2016 .

[10]  Jostein R. Natvig,et al.  Constrained pressure residual multiscale (CPR-MS) method for fully implicit simulation of multiphase flow in porous media , 2015, J. Comput. Phys..

[11]  Bo Lu,et al.  Iteratively coupled reservoir simulation for multiphase flow in porous media , 2008 .

[12]  Jiamin Jiang,et al.  Efficient C1-continuous phase-potential upwind (C1-PPU) schemes for coupled multiphase flow and transport with gravity , 2017 .

[13]  H. Tchelepi,et al.  Multi-scale finite-volume method for elliptic problems in subsurface flow simulation , 2003 .

[14]  Yalchin Efendiev,et al.  Accurate multiscale finite element methods for two-phase flow simulations , 2006, J. Comput. Phys..

[15]  Hamdi A. Tchelepi,et al.  Analysis of Hybrid Upwinding for Fully-Implicit Simulation of Three-Phase Flow with Gravity , 2016, SIAM J. Numer. Anal..

[16]  Yalchin Efendiev,et al.  C1-Continuous relative permeability and hybrid upwind discretization of three phase flow in porous media , 2016 .

[17]  Patrick Jenny,et al.  Unconditionally convergent nonlinear solver for hyperbolic conservation laws with S-shaped flux functions , 2009, J. Comput. Phys..

[18]  K.-A. Lie,et al.  Successful application of multiscale methods in a real reservoir simulator environment , 2016, Computational Geosciences.

[19]  Jan Vierendeels,et al.  Partitioned simulation of the interaction between an elastic structure and free surface flow , 2010 .

[20]  Hamdi A. Tchelepi,et al.  Nonlinear analysis of multiphase transport in porous media in the presence of viscous, buoyancy, and capillary forces , 2015, J. Comput. Phys..

[21]  Bradley T. Mallison,et al.  Implicit Hybrid Upwind scheme for coupled multiphase flow and transport with buoyancy , 2016 .

[22]  Patrick Jenny,et al.  Iterative multiscale finite-volume method , 2008, J. Comput. Phys..

[23]  A. Kozlova,et al.  A Real-Field Multiscale Black-Oil Reservoir Simulator , 2015, ANSS 2015.

[24]  Hamdi A. Tchelepi,et al.  Adaptively Localized Continuation-Newton Method—Nonlinear Solvers That Converge All the Time , 2010 .

[25]  K. Aziz,et al.  Petroleum Reservoir Simulation , 1979 .

[26]  L. Baker Three-Phase Relative Permeability Correlations , 1988 .

[27]  Olav Møyner,et al.  A Multiscale Restriction-Smoothed Basis Method for Compressible Black-Oil Models , 2016 .

[28]  Rainer Helmig,et al.  EfficientModeling of Flow and Transport in Porous Media Using Multiphysics andMultiscale Approaches , 2010 .

[29]  Yalchin Efendiev,et al.  Generalized multiscale finite element methods (GMsFEM) , 2013, J. Comput. Phys..

[30]  Olav Møyner,et al.  A multiscale two-point flux-approximation method , 2014, J. Comput. Phys..

[31]  Patrick Jenny,et al.  Adaptive fully implicit multi-scale finite-volume method for multi-phase flow and transport in heterogeneous porous media , 2006, J. Comput. Phys..

[32]  Hamdi A. Tchelepi,et al.  Trust-region based solver for nonlinear transport in heterogeneous porous media , 2013, J. Comput. Phys..

[33]  Alexander Düster,et al.  Accelerated staggered coupling schemes for problems of thermoelasticity at finite strains , 2012, Comput. Math. Appl..

[34]  Patrick Jenny,et al.  Modified sequential fully implicit scheme for compositional flow simulation , 2017, J. Comput. Phys..

[35]  Jostein R. Natvig,et al.  A Multiscale Mixed Finite Element Solver for Three Phase Black Oil Flow , 2009 .

[36]  Keith H. Coats,et al.  A Note on Impes and Some Impes-Based Simulation Models , 1999 .

[37]  Patrick Jenny,et al.  Multiscale finite-volume method for parabolic problems arising from compressible multiphase flow in porous media , 2009, J. Comput. Phys..

[38]  Hamdi A. Tchelepi,et al.  Algebraic multiscale solver for flow in heterogeneous porous media , 2014, J. Comput. Phys..

[39]  Homer F. Walker,et al.  Anderson Acceleration for Fixed-Point Iterations , 2011, SIAM J. Numer. Anal..

[40]  H. Tchelepi,et al.  Multiscale finite-volume formulation for multiphase flow in porous media: black oil formulation of compressible, three-phase flow with gravity , 2008 .

[41]  P. Jenny,et al.  Multiscale finite-volume method for density-driven flow in porous media , 2008 .

[42]  Roland Masson,et al.  Accelerating the convergence of coupled geomechanical‐reservoir simulations , 2007 .

[43]  Zhangxin Chen,et al.  The sequential method for the black-oil reservoir simulation on unstructured grids , 2003 .