Multi-scale modeling of three-phase three-component processes in heterogeneous porous media

Abstract Flow and transport processes in porous media occur on different spatial and temporal scales and may also be locally different. Additionally, the structure of the porous medium itself generally shows a high dependence on the spatial scale. As an example, the contamination of the unsaturated zone with a light non-aqueous phase liquid is studied, corresponding to a domain with randomly distributed heterogeneities where complex three-phase–three-component processes are relevant only in a small (local) subdomain. This subdomain needs fine resolution as the complex processes are governed by small-scale effects. For a comprehensive fine-scale model taking into account three-phase–three-component processes as well as heterogeneities in the whole (global) model domain, data collection is expensive and computational time is long. Therefore, we developed a general multi-scale concept where on the one hand, the global flow field influences the local three-phase–three-component processes on the fine-scale. On the other hand, a coarse-scale saturation equation is solved where the effects of the fine-scale multi-phase–multi-component processes in the subdomain are captured by source/sink terms and the effects of fine-scale heterogeneities by a macrodispersion term. It turned out that the new multi-scale algorithm represents a flexible and extendable tool for incorporating processes of different complexity occurring at different locations in one model domain, while reducing the amount of required data. Using a simplified numerical example, it could be shown that the multi-scale algorithm functions very well. However, more research has to be done in order to further improve its computational efficiency.

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

[2]  Rainer Helmig,et al.  Numerical simulation of non-isothermal multiphase multicomponent processes in porous media. 2. Applications for the injection of steam and air , 2002 .

[3]  Chi-Wang Shu,et al.  The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. IV. The multidimensional case , 1990 .

[4]  I. Babuska,et al.  A DiscontinuoushpFinite Element Method for Diffusion Problems , 1998 .

[5]  Yalchin Efendiev,et al.  Numerical modeling of subgrid heterogeneity in two phase flow simulations , 2002 .

[6]  Michel Fortin,et al.  Mixed and Hybrid Finite Element Methods , 2011, Springer Series in Computational Mathematics.

[7]  Ivan Yotov,et al.  Advanced Techniques and Algorithms for Reservoir Simulation, IV: Multiblock Solvers and Preconditioners , 2002 .

[8]  Karsten Pruess,et al.  A Fickian diffusion model for the spreading of liquid plumes infiltrating in heterogeneous media , 1996 .

[9]  S. E. Buckley,et al.  Mechanism of Fluid Displacement in Sands , 1942 .

[10]  Holger Class.  Theorie und numerische Modellierung nichtisothermer Mehrphasenprozesse in NAPL-kontaminierten porösen Medien , 2001 .

[11]  L. Young,et al.  A Generalized Compositional Approach for Reservoir Simulation , 1983 .

[12]  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 .

[13]  Rainer Helmig,et al.  Multi‐scale modelling of two‐phase–two‐component processes in heterogeneous porous media , 2006, Numer. Linear Algebra Appl..

[14]  Cass T. Miller,et al.  Multiphase flow and transport modeling in heterogeneous porous media: challenges and approaches , 1998 .

[15]  Ronald W. Falta,et al.  Modeling sub‐grid‐block‐scale dense nonaqueous phase liquid (DNAPL) pool dissolution using a dual‐domain approach , 2003 .

[16]  J. Oden,et al.  A discontinuous hp finite element method for convection—diffusion problems , 1999 .

[17]  A. Bellin,et al.  On the use of block‐effective macrodispersion for numerical simulations of transport in heterogeneous formations , 2003 .

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

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

[20]  D. W. Pollock.  Semianalytical Computation of Path Lines for Finite‐Difference Models , 1988 .

[21]  J. Gómez-Hernández,et al.  Upscaling hydraulic conductivities in heterogeneous media: An overview , 1996 .

[22]  Chi-Wang Shu,et al.  The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V , 1998 .

[23]  Karsten Pruess,et al.  Numerical modeling of steam injection for the removal of nonaqueous phase liquids from the subsurface: 2. Code validation and application , 1992 .

[24]  Jennifer Niessner,et al.  Multi-scale modeling of multi-phase-multi-component processes in heterogeneous porous media , 2006 .

[25]  William G. Gray,et al.  Mechanics and thermodynamics of multiphase flow in porous media including interphase boundaries , 1990 .

[26]  Richard E. Ewing,et al.  Resource recovery, confinement, and remediation of environmental hazards , 2002 .

[27]  Yalchin Efendiev,et al.  Modeling of subgrid effects in coarse‐scale simulations of transport in heterogeneous porous media , 2000 .

[28]  Du Plessis,et al.  Fluid Transport in Porous Media , 1997 .

[29]  A. Bellin,et al.  The concept of block-effective macrodispersivity and a unified approach for grid-scale- and plume-scale-dependent transport , 1999, Journal of Fluid Mechanics.

[30]  P. Renard,et al.  Calculating equivalent permeability: a review , 1997 .

[31]  C. L. Farmer,et al.  Upscaling: a review , 2002 .

[32]  Rainer Helmig,et al.  Heat transport in the unsaturated zone – comparison of experimental results and numerical simulations , 1998 .

[33]  A. V. Balakrishnan,et al.  Vistas in applied mathematics : numerical analysis, atmospheric sciences, immunology , 1986 .

[34]  D. D. Fussell,et al.  AN ITERATIVE TECHNIQUE FOR COMPOSITIONAL RESERVOIR MODELS , 1979 .

[35]  William G. Gray,et al.  Toward an improved description of the physics of two-phase flow , 1993 .

[36]  E Weinan,et al.  Heterogeneous multiscale method: A general methodology for multiscale modeling , 2003 .

[37]  R. Helmig.  Multiphase Flow and Transport Processes in the Subsurface: A Contribution to the Modeling of Hydrosystems , 2011 .

[38]  Steven E. Koch,et al.  A Survey of Nested Grid Techniques and Their Potential for Use Within the MASS Weather , 1987 .

[39]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

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

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

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

[43]  Sophie Didierjean,et al.  Analytical solutions of one-dimensional macrodispersion in stratified porous media by the quadrupole method: convergence to an equivalent homogeneous porous medium , 2004 .

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

[45]  Rainer Helmig,et al.  Numerical simulation of non-isothermal multiphase multicomponent processes in porous media.: 1. An efficient solution technique , 2002 .

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

[47]  Magne S. Espedal,et al.  Macrodispersion for two-phase, immiscible flow in porous media , 1994 .

[48]  Jean-Louis Auriault,et al.  UPSCALING HETEROGENEOUS MEDIA BY ASYMPTOTIC EXPANSIONS , 2002 .

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

[50]  P. Witherspoon,et al.  Numerical modeling of steam injection for the removal of nonaqueous phase liquids from the subsurface. 1. Numerical formulation , 1992 .

[51]  W. G. Gray,et al.  Recent advances in theories of two-phase flow in porous media , 1997 .

[52]  Jack C. Parker,et al.  A parametric model for constitutive properties governing multiphase flow in porous media , 1987 .

[53]  William G. Gray,et al.  General conservation equations for multi-phase systems: 3. Constitutive theory for porous media flow. , 1980 .