Models and simulations of variable-density flow in fractured porous media

We develop a numerical technique for variable-density flow in fractured porous media, in which fractures are (d − 1)-dimensional manifolds, with d being the dimension of the ambient space. The PDEs of variable-density flow are firstly presented in the same form for both the fractures and the enclosing medium. Then, the equations defined in the fractures are averaged along the fracture width and formulated in (d − 1)-dimensions. The resulting PDEs are solved together with those defined in the enclosing medium, which maintain their d-dimensional form. The discretisation of the coupled system of d-and (d − 1)-dimensional PDEs follows a finite-volume method requiring a special construction of the discretisation grid, obtained by the algorithm explained in this paper. The accuracy of our technique is tested by comparing the produced results with those obtained in simulations in which the fractures maintain dimension d. In all simulations the fractured medium is three-dimensional.

[1]  Roussos Dimitrakopoulos,et al.  An efficient method for discretizing 3D fractured media for subsurface flow and transport simulations , 2011 .

[2]  Joke Blom,et al.  Computing brine transport in porous media with an adaptive‐grid method , 1992 .

[3]  T. Graf,et al.  Grid convergence of variable-density flow simulations in discretely-fractured porous media , 2011 .

[4]  Klaus Gärtner,et al.  Meshing Piecewise Linear Complexes by Constrained Delaunay Tetrahedralizations , 2005, IMR.

[5]  W. Hackbusch.  Iterative Solution of Large Sparse Systems of Equations , 1993 .

[6]  J. Bear.  Hydraulics of Groundwater , 1979 .

[7]  W. Kinzelbach,et al.  Numerical simulation of three-dimensional saltwater-freshwater fingering instabilities observed in a porous medium , 2006 .

[8]  G. Wittum,et al.  The saltpool benchmark problem – numerical simulation of saltwater upconing in a porous medium , 2002 .

[9]  Alexander Fuchs.  Almost Regular Triangulations of Trimmend NURBS-Solids , 2001, Engineering with Computers.

[10]  S. P. Neuman,et al.  Trends, prospects and challenges in quantifying flow and transport through fractured rocks , 2005 .

[11]  H. Borouchaki,et al.  Geometric surface mesh optimization , 1998 .

[12]  A Two-Dimensional Areal Model for Density Dependent Flow Regime , 2001 .

[13]  Richard E. Ewing,et al.  Numerical Simulation of Multiphase Flow in Fractured Porous Media , 2000 .

[14]  L. Paul Chew,et al.  Constrained Delaunay triangulations , 1987, SCG '87.

[15]  Richard Barrett,et al.  Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods , 1994, Other Titles in Applied Mathematics.

[16]  J. Sharp,et al.  Heterogeneity effects on possible salinity‐driven free convection in low‐permeability strata , 2009 .

[17]  Stefan Lang,et al.  Large‐scale density‐driven flow simulations using parallel unstructured Grid adaptation and local multigrid methods , 2005, Concurr. Pract. Exp..

[18]  Jens Lang,et al.  THREE-DIMENSIONAL ADAPTIVE COMPUTATION OF BRINE TRANSPORT IN POROUS MEDIA , 2002 .

[19]  Jonathan Richard Shewchuk,et al.  Constrained Delaunay Tetrahedralizations and Provably Good Boundary Recovery , 2002, IMR.

[20]  Zhiqiang Cai,et al.  On the finite volume element method , 1990 .

[21]  Peter Frolkovič,et al.  Maximum principle and local mass balance for numerical solutions of transport equation coupled with variable density flow. , 1998 .

[22]  C. Oldenburg,et al.  Layered Thermohaline Convection in Hypersaline Geothermal Systems , 1998 .

[23]  H. Murphy.  Convective instabilities in vertical fractures and faults , 1979 .

[24]  E. Sudicky,et al.  Density-dependent solute transport in discretely-fractured geologic media: is prediction possible? , 1998 .

[25]  William G. Gray,et al.  Derivation of conditions describing transport across zones of reduced dynamics within multiphase systems , 1989 .

[26]  Gabriel Wittum,et al.  Simulation of density-driven flow in fractured porous media , 2010 .

[27]  J. Bear.  Dynamics of Fluids in Porous Media , 1975 .

[28]  Jesús Carrera,et al.  A methodology to interpret cross-hole tests in a granite block , 2006 .

[29]  G. Wittum,et al.  Simulation of salinity- and thermohaline-driven flow in fractured porous media , 2012 .

[30]  R. Therrien,et al.  Variable-density groundwater flow and solute transport in porous media containing nonuniform discrete fractures , 2005 .

[31]  Olaf Kolditz,et al.  Variable-density flow and transport in porous media: approaches and challenges , 2002 .

[32]  Philippe Angot,et al.  ASYMPTOTIC AND NUMERICAL MODELLING OF FLOWS IN FRACTURED POROUS MEDIA , 2009 .

[33]  T. Clement,et al.  Improving the worthiness of the Henry problem as a benchmark for density‐dependent groundwater flow models , 2004 .

[34]  R. Therrien,et al.  Stable-unstable flow of geothermal fluids in fractured rock , 2009 .

[35]  S. M. H. Karimian,et al.  Pressure-based control-volume finite element method for flow at all speeds , 1995 .

[36]  Chin-Fu Tsang,et al.  Flow and Contaminant Transport in Fractured Rock , 1993 .