Coupling geological and numerical models to simulate groundwater flow and contaminant transport in fractured media

A new modeling approach is presented to improve numerical simulations of groundwater flow and contaminant transport in fractured geological media. The approach couples geological and numerical models through an intermediate mesh generation phase. As a first step, a platform for 3D geological modeling is used to represent fractures as 2D surfaces with arbitrary shape and orientation in 3D space. The advantage of the geological modeling platform is that 2D triangulated fracture surfaces are modeled and visualized before building a 3D mesh. The triangulated fractures are then transferred to the mesh generation software that discretizes the 3D simulation domain with tetrahedral elements. The 2D triangular fracture elements do not cut through the 3D tetrahedral elements, but they rather form interfaces with them. The tetrahedral mesh is then used for 3D groundwater flow and contaminant transport simulations in discretely fractured porous media. The resulting mesh for the 2D fractures and 3D rock matrix is checked to ensure that there are no negative transmissibilities in the discretized flow and transport equation, to avoid unrealistic results. To test the validity of the approach, flow and transport simulations for a tetrahedral mesh are compared to simulations using a block-based mesh and with results of an analytical solution. The fluid conductance matrix for the tetrahedral mesh is also analyzed and compared with known matrix values.

[1]  Siegfried Selberherr,et al.  On the interplay between meshing and discretization inthree-dimensional diffusion simulation , 2000, IEEE Trans. Comput. Aided Des. Integr. Circuits Syst..

[2]  René Therrien,et al.  A method to discretize non‐planar fractures for 3D subsurface flow and transport simulations , 2008 .

[3]  Takeo Taniguchi,et al.  Numerical experiments for 3-dimensional flow analysis in a fractured rock with porous matrix , 1996 .

[4]  Frank W. Letniowski,et al.  Three-Dimensional Delaunay Triangulations for Finite Element Approximations to a Second-Order Diffusion Operator , 1992, SIAM J. Sci. Comput..

[5]  S. Geiger,et al.  Combining finite element and finite volume methods for efficient multiphase flow simulations in highly heterogeneous and structurally complex geologic media , 2004 .

[6]  Olaf Kolditz,et al.  Development and application of a CAD interface for fractured rock , 2005 .

[7]  P. F. Andersen,et al.  Saltwater intrusion in aquifers: Development and testing of a three‐dimensional finite element model , 1987 .

[8]  C. W. Gable,et al.  Automated grid generation from models of complex geologic structure and stratigraphy , 1996 .

[9]  P. Huyakorn,et al.  Techniques for Making Finite Elements Competitve in Modeling Flow in Variably Saturated Porous Media , 1984 .

[10]  Thomas Graf,et al.  Modeling coupled thermohaline flow and reactive solute transport in discretely-fractured porous media , 2006 .

[11]  Olaf Kolditz,et al.  Numerical analysis of coupled hydrosystems based on an object-oriented compartment approach , 2008 .

[12]  Mario Putti,et al.  Finite Element Approximation of the Diffusion Operator on Tetrahedra , 1998, SIAM J. Sci. Comput..

[13]  E. A. Sudicky,et al.  Contaminant transport in fractured porous media: Analytical solution for a single fracture , 1981 .

[14]  E. Sudicky,et al.  Three-dimensional analysis of variably-saturated flow and solute transport in discretely-fractured porous media , 1996 .

[15]  Peter A. Forsyth,et al.  A control volume finite element method for three‐dimensional NAPL groundwater contamination , 1991 .

[16]  Paul E. Allaire.  Basics of the Finite Element Method: Solid Mechanics, Heat Transfer, and Fluid Mechanics , 1985 .

[17]  H. Trease,et al.  Geological applications of automatic grid generation tools for finite elements applied to porous flow modeling , 1996 .

[18]  Sebastian Bauer,et al.  A process-oriented approach to computing multi-field problems in porous media , 2004 .

[19]  G. Zyvoloski,et al.  Grid Resolution Study of Ground Water Flow and Transport , 2005, Ground water.

[20]  Mario Putti,et al.  Accuracy of Galerkin finite elements for groundwater flow simulations in two and three‐dimensional triangulations , 2001 .

[21]  Emil O. Frind,et al.  Contaminant transport in fractured porous media: Analytical solutions for a system of parallel fractures , 1982 .

[22]  P. Huyakorn,et al.  A three‐dimensional finite‐element model for simulating water flow in variably saturated porous media , 1986 .