High order discontinuous Galerkin method for simulating miscible flooding in porous media

We present a high-order method for miscible displacement simulation in porous media. The method is based on discontinuous Galerkin discretization with weighted average stabilization technique and flux reconstruction post processing. The mathematical model is decoupled and solved sequentially. We apply domain decomposition and algebraic multigrid preconditioner for the linear system resulting from the high-order discretization. The accuracy and robustness of the method are demonstrated in the convergence study with analytical solutions and heterogeneous porous media, respectively. We also investigate the effect of grid orientation and anisotropic permeability using high-order discontinuous Galerkin method in contrast with cell-centered finite volume method. The study of the parallel implementation shows the scalability and efficiency of the method on parallel architecture. We also verify the simulation result on highly heterogeneous permeability field from the SPE10 model.

[1]  Cornelis Vuik,et al.  Fast linear solver for diffusion problems with applications to pressure computation in layered domains , 2014, Computational Geosciences.

[2]  Ludmil T. Zikatanov,et al.  Two‐level preconditioning of discontinuous Galerkin approximations of second‐order elliptic equations , 2006, Numer. Linear Algebra Appl..

[3]  Y. Epshteyn,et al.  CONVERGENCE OF HIGH ORDER METHODS FOR MISCIBLE DISPLACEMENT , 2008 .

[4]  Peter Bastian,et al.  A fully-coupled discontinuous Galerkin method for two-phase flow in porous media with discontinuous capillary pressure , 2013, Computational Geosciences.

[5]  K. Nikitin,et al.  A Monotone Non-linear Finite Volume Method for Advection-diffusion Equations and Multiphase Flows , 2012 .

[6]  Jizhou Li,et al.  Convergence of a high order method in time and space for the miscible displacement equations , 2015 .

[7]  Robert Scheichl,et al.  Algebraic multigrid for discontinuous Galerkin discretizations of heterogeneous elliptic problems , 2012, Numer. Linear Algebra Appl..

[8]  Khalid Aziz,et al.  Use of Voronoi Grid in Reservoir Simulation , 1994 .

[9]  Yu. V. Vassilevski,et al.  A monotone nonlinear finite volume method for diffusion equations on conformal polyhedral meshes , 2009 .

[10]  Minfu Feng,et al.  A stabilized mixed discontinuous Galerkin method for the incompressible miscible displacement problem , 2011 .

[11]  M. R. Todd,et al.  Methods for Increased Accuracy in Numerical Reservoir Simulators , 1972 .

[12]  A. Ern,et al.  A discontinuous Galerkin method with weighted averages for advection–diffusion equations with locally small and anisotropic diffusivity , 2008 .

[13]  Hussein Hoteit,et al.  New two‐dimensional slope limiters for discontinuous Galerkin methods on arbitrary meshes , 2004 .

[14]  Peter Bastian,et al.  Generic implementation of finite element methods in the Distributed and Unified Numerics Environment (DUNE) , 2010, Kybernetika.

[15]  M. Wheeler,et al.  Discontinuous Galerkin methods for coupled flow and reactive transport problems , 2005 .

[16]  Sören Bartels,et al.  Discontinuous Galerkin Finite Element Convergence for Incompressible Miscible Displacement Problems of Low Regularity , 2009, SIAM J. Numer. Anal..

[17]  Bernardo Cockburn,et al.  The Runge-Kutta local projection discontinous Galerkin finite element method for conservation laws , 1990 .

[18]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous galerkin finite element method for conservation laws. II: General framework , 1989 .

[19]  Luciane A. Schuh,et al.  Discontinuous Galerkin approximation of two-phase flows in heterogeneous porous media with discontinuous capillary pressures , 2010 .

[20]  Béatrice Rivière,et al.  A high order method for solving the black-oil problem in porous media , 2015 .

[21]  E. Koval,et al.  A Method for Predicting the Performance of Unstable Miscible Displacement in Heterogeneous Media , 1963 .

[22]  Yuri V. Vassilevski,et al.  A monotone nonlinear finite volume method for diffusion equations and multiphase flows , 2014, Computational Geosciences.

[23]  S. Nicaise,et al.  An accurate H(div) flux reconstruction for discontinuous Galerkin approximations of elliptic problems , 2007 .

[24]  Andreas Dedner,et al.  A generic grid interface for parallel and adaptive scientific computing. Part II: implementation and tests in DUNE , 2008, Computing.

[25]  Huiqin Lu,et al.  Eigenvalue Problem and Unbounded Connected Branch of Positive Solutions to a Class of Singular Elastic Beam Equations , 2011 .

[26]  Guglielmo Scovazzi,et al.  A high-order, fully coupled, upwind, compact discontinuous Galerkin method for modeling of viscous fingering in compressible porous media , 2013 .

[27]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: one-dimensional systems , 1989 .

[28]  Richard E. Ewing,et al.  Efficient Time-Stepping Methods for Miscible Displacement Problems in Porous Media , 1982 .

[29]  D. Moissis,et al.  Simulation of viscous fingering during miscible displacement in nonuniform porous media , 1988 .

[30]  Markus Blatt A Parallel Algebraic Multigrid Method for Elliptic Problems with Highly Discontinuous Coefficients , 2010 .

[31]  B. Rivière,et al.  Superconvergence and H(div) projection for discontinuous Galerkin methods , 2003 .

[32]  Xiaojie Lin,et al.  Existence of Solutions to a Nonlocal Boundary Value Problem with Nonlinear Growth , 2011 .

[33]  H. Huang,et al.  A fully-coupled upwind discontinuous Galerkin method for incompressible porous media flows: High-order computations of viscous fingering instabilities in complex geometry , 2013, J. Comput. Phys..

[34]  Peter Bastian,et al.  Higher Order Discontinuous Galerkin Methods for Flow and Transport in Porous Media , 2003 .

[35]  Todd Arbogast,et al.  A discontinuous Galerkin method for two-phase flow in a porous medium enforcing H(div) velocityand continuous capillary pressure , 2013, Computational Geosciences.

[36]  Andreas Dedner,et al.  A generic grid interface for parallel and adaptive scientific computing. Part I: abstract framework , 2008, Computing.

[37]  Max Jensen,et al.  Stable Crank–Nicolson Discretisation for Incompressible Miscible Displacement Problems of Low Regularity , 2009, 0910.1516.

[38]  Igor Mozolevski,et al.  Accurate velocity reconstruction for Discontinuous Galerkin approximations of two-phase porous media flows , 2009 .

[39]  B. Rivière,et al.  Discontinuous Galerkin methods for flow and transport problems in porous media , 2001 .

[40]  Jizhou Li,et al.  Numerical solutions of the incompressible miscible displacement equations in heterogeneous media , 2015 .