Higher-resolution convection schemes for flow in porous media on highly distorted unstructured grids

Higher-resolution schemes are presented for convective flow approximation on highly distorted unstructured grids. The schemes are coupled with continuous full-tensor Darcy-flux approximations. A sequence of non-uniform and distorted grid formulations are developed and compared for a range of unstructured meshes with variable grid spacing. The higher-order schemes are constructed using non-uniform grid slope limiters such that they are stable with a local maximum principle, ensuring that solutions are free of spurious oscillations. Benefits of the resulting schemes are demonstrated for classical test problems in reservoir simulation including cases with full-tensor permeability fields. The test cases involve a range of unstructured grids with variations in grid spacing, orientation and permeability that lead to flow fields that are poorly resolved by standard simulation methods. The higher-order formulations are shown to effectively reduce numerical diffusion, leading to improved resolution of concentration and saturation fronts.

[1]  P. Sweby.  High Resolution Schemes Using Flux Limiters for Hyperbolic Conservation Laws , 1984 .

[2]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[3]  Mayur Pal,et al.  Convergence study of a family of flux‐continuous, finite‐volume schemes for the general tensor pressure equation , 2006 .

[4]  Michael G. Edwards,et al.  Higher‐resolution hyperbolic‐coupled‐elliptic flux‐continuous CVD schemes on structured and unstructured grids in 3‐D , 2006 .

[5]  Michael G. Edwards,et al.  Dynamically Adaptive Godunov Schemes With Renormalization in Reservoir Simulation , 1993 .

[6]  Michael G. Edwards,et al.  M-Matrix Flux Splitting for General Full Tensor Discretization Operators on Structured and Unstructured Grids , 2000 .

[7]  Michael G. Edwards.  Cross flow tensors and finite volume approximation with by deferred correction , 1998 .

[8]  Patrick Jenny,et al.  A Finite-Volume Method with Hexahedral Multiblock Grids for Modeling Flow in Porous Media , 2002 .

[9]  L. Durlofsky.  A triangle based mixed finite element–finite volume technique for modeling two phase flow through porous media , 1993 .

[10]  Ivar Aavatsmark,et al.  Discretization on Unstructured Grids for Inhomogeneous, Anisotropic Media. Part I: Derivation of the Methods , 1998, SIAM J. Sci. Comput..

[11]  Michael G. Edwards,et al.  Physically Based Higher Order Godunov Schemes for Compositional Simulation , 2001 .

[12]  K. Morgan,et al.  A review and comparative study of upwind biased schemes for compressible flow computation. Part III: Multidimensional extension on unstructured grids , 2002 .

[13]  Michael G. Edwards,et al.  Unstructured, Control-Volume Distributed, Full-Tensor Finite-Volume Schemes with Flow Based Grids , 2002 .

[14]  M. Edwards.  Higher order hyperbolic - coupled - elliptic flux-continuous CVD finite volume schemes in two and three dimensions , 2006 .

[15]  A. Jameson.  Artificial diffusion, upwind biasing, limiters and their effect on accuracy and multigrid convergence in transonic and hypersonic flows , 1993 .

[16]  Phillip Colella,et al.  Higher order Godunov methods for general systems of hyperbolic conservation laws , 1989 .

[17]  Martin J. Blunt,et al.  Tenth SPE Comparative Solution Project: a comparison of upscaling techniques , 2001 .

[18]  M. Blunt,et al.  Implicit flux limiting schemes for petroleum reservoir simulation , 1992 .

[19]  Khalid Aziz,et al.  A Control Volume Scheme for Flexible Grids in Reservoir Simulation , 1997 .

[20]  Michael G. Edwards,et al.  Higher‐resolution hyperbolic‐coupled‐elliptic flux‐continuous CVD schemes on structured and unstructured grids in 2‐D , 2006 .

[21]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[22]  Michael G. Edwards,et al.  A higher-order Godunov scheme coupled with dynamic local grid refinement for flow in a porous medium , 1996 .

[23]  Michael Andrew Christie,et al.  Tenth SPE Comparative Solution Project: a comparison of upscaling techniques , 2001 .

[24]  Michael G. Edwards,et al.  Non-Upwind Versus Upwind Schemes for Hyperbolic Conservation Laws in Porous Media , 2005 .