Accurate multiscale finite element methods for two-phase flow simulations

In this paper we propose a modified multiscale finite element method for two-phase flow simulations in heterogeneous porous media. The main idea of the method is to use the global fine-scale solution at initial time to determine the boundary conditions of the basis functions. This method provides a significant improvement in two-phase flow simulations in porous media where the long-range effects are important. This is typical for some recent benchmark tests, such as the SPE comparative solution project [M. Christie, M. Blunt, Tenth spe comparative solution project: a comparison of upscaling techniques, SPE Reser. Eval. Eng. 4 (2001) 308-317], where porous media have a channelized structure. The use of global information allows us to capture the long-range effects more accurately compared to the multiscale finite element methods that use only local information to construct the basis functions. We present some analysis of the proposed method to illustrate that the method can indeed capture the long-range effect in channelized media.

[1]  L. Durlofsky.  Numerical calculation of equivalent grid block permeability tensors for heterogeneous porous media , 1991 .

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

[3]  L. Durlofsky,et al.  A coupled local-global upscaling approach for simulating flow in highly heterogeneous formations , 2003 .

[4]  Thomas Y. Hou,et al.  A Multiscale Finite Element Method for Elliptic Problems in Composite Materials and Porous Media , 1997 .

[5]  Patrick Jenny,et al.  Adaptive Multiscale Finite-Volume Method for Multiphase Flow and Transport in Porous Media , 2005, Multiscale Model. Simul..

[6]  T. Arbogast.  Implementation of a Locally Conservative Numerical Subgrid Upscaling Scheme for Two-Phase Darcy Flow , 2002 .

[7]  E Weinan,et al.  The Heterogeneous Multi-Scale Method , 2002 .

[8]  F. Brezzi.  Interacting with the subgrid world , 2005 .

[9]  Jørg E. Aarnes,et al.  On the Use of a Mixed Multiscale Finite Element Method for GreaterFlexibility and Increased Speed or Improved Accuracy in Reservoir Simulation , 2004, Multiscale Model. Simul..

[10]  T. Hou,et al.  Multiscale Finite Element Methods for Nonlinear Problems and Their Applications , 2004 .

[11]  Louis J. Durlofsky,et al.  Coarse scale models of two phase flow in heterogeneous reservoirs: volume averaged equations and their relationship to existing upscaling techniques , 1998 .

[12]  I. Babuska,et al.  Generalized Finite Element Methods: Their Performance and Their Relation to Mixed Methods , 1983 .

[13]  Giancarlo Sangalli,et al.  Capturing Small Scales in Elliptic Problems Using a Residual-Free Bubbles Finite Element Method , 2003, Multiscale Model. Simul..

[14]  T. Hughes,et al.  The variational multiscale method—a paradigm for computational mechanics , 1998 .

[15]  Zhiming Chen,et al.  A mixed multiscale finite element method for elliptic problems with oscillating coefficients , 2003, Math. Comput..

[16]  C. W. Gear,et al.  Equation-Free, Coarse-Grained Multiscale Computation: Enabling Mocroscopic Simulators to Perform System-Level Analysis , 2003 .

[17]  C. Schwab,et al.  Homogenization via p -FEM for problems with microstructure , 2000 .

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

[19]  Clayton V. Deutsch,et al.  GSLIB: Geostatistical Software Library and User's Guide , 1993 .

[20]  I. Babuska,et al.  Special finite element methods for a class of second order elliptic problems with rough coefficients , 1994 .

[21]  E Weinan,et al.  The Heterognous Multiscale Methods , 2003 .

[22]  Michael Andrew Christie,et al.  Upscaling for reservoir simulation , 1996 .

[23]  Thomas Y. Hou,et al.  Convergence of a Nonconforming Multiscale Finite Element Method , 2000, SIAM J. Numer. Anal..

[24]  J. W. Barker,et al.  A critical review of the use of pseudo-relative permeabilities for upscaling , 1997 .