dfnWorks: A discrete fracture network framework for modeling subsurface flow and transport

Abstract dfn W orks is a parallelized computational suite to generate three-dimensional discrete fracture networks (DFN) and simulate flow and transport. Developed at Los Alamos National Laboratory over the past five years, it has been used to study flow and transport in fractured media at scales ranging from millimeters to kilometers. The networks are created and meshed using dfn G en , which combines fram (the feature rejection algorithm for meshing) methodology to stochastically generate three-dimensional DFNs with the L a G ri T meshing toolbox to create a high-quality computational mesh representation. The representation produces a conforming Delaunay triangulation suitable for high performance computing finite volume solvers in an intrinsically parallel fashion. Flow through the network is simulated in dfn F low , which utilizes the massively parallel subsurface flow and reactive transport finite volume code pflotran . A Lagrangian approach to simulating transport through the DFN is adopted within dfn T rans to determine pathlines and solute transport through the DFN. Example applications of this suite in the areas of nuclear waste repository science, hydraulic fracturing and CO2 sequestration are also included.

[1]  Satish Karra,et al.  Three-phase numerical model for subsurface hydrology in permafrost-affected regions (PFLOTRAN-ICE v1.0) , 2014 .

[2]  Satish Karra,et al.  Effect of advective flow in fractures and matrix diffusion on natural gas production , 2015 .

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

[4]  Cesare Marchetti,et al.  The future of natural gas: A Darwinian analysis , 1987 .

[5]  J. Erhel,et al.  A mixed hybrid Mortar method for solving flow in discrete fracture networks , 2010 .

[6]  Jean-Raynald de Dreuzy,et al.  Influence of spatial correlation of fracture centers on the permeability of two‐dimensional fracture networks following a power law length distribution , 2004 .

[7]  C. C. Law,et al.  ParaView: An End-User Tool for Large-Data Visualization , 2005, The Visualization Handbook.

[8]  P C Lichtner,et al.  High resolution numerical investigation on the effect of convective instability on long term CO2 storage in saline aquifers , 2007 .

[9]  Charles Hansen,et al.  The Visualization Handbook , 2011 .

[10]  G. Marsily,et al.  Modeling fracture flow with a stochastic discrete fracture network: calibration and validation: 1. The flow model , 1990 .

[11]  S Pacala,et al.  Stabilization Wedges: Solving the Climate Problem for the Next 50 Years with Current Technologies , 2004, Science.

[12]  David M. Mount,et al.  A point-placement strategy for conforming Delaunay tetrahedralization , 2000, SODA '00.

[13]  Vladimir Cvetkovic,et al.  Numerical and analytical modeling of advective travel times in realistic three‐dimensional fracture networks , 2011 .

[14]  Jean-Raynald de Dreuzy,et al.  A Generalized Mixed Hybrid Mortar Method for Solving Flow in Stochastic Discrete Fracture Networks , 2012, SIAM J. Sci. Comput..

[15]  H. Coxeter,et al.  Introduction to Geometry. , 1961 .

[16]  Laurence C. Hull,et al.  Streamline routing through fracture junctions , 1986 .

[17]  Velimir V. Vesselinov,et al.  MODEL ANALYSIS AND DECISION SUPPORT (MADS) FOR COMPLEX PHYSICS MODELS , 2012 .

[18]  Kanti V. Mardia,et al.  A Connectivity Index for Discrete Fracture Networks , 2007 .

[19]  George A Zyvoloski,et al.  An Investigation of Numerical Grid Effects in Parameter Estimation , 2006, Ground water.

[20]  Satish Karra,et al.  Three-phase numerical model for subsurface hydrology in permafrost-affected regions , 2014 .

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

[22]  Scott L. Painter,et al.  Time domain particle tracking methods for simulating transport with retention and first‐order transformation , 2008 .

[23]  Martin Vohralík,et al.  Numerical simulation of fracture flow with a mixed-hybrid FEM stochastic discrete fracture network model , 2005 .

[24]  Wenqing Wang,et al.  Geometric modelling and object-oriented software concepts applied to a heterogeneous fractured network from the Grimsel rock laboratory , 2007 .

[25]  Chuan Lu,et al.  Simulating subsurface flow and transport on ultrascale computers using PFLOTRAN , 2007 .

[26]  Jean-Raynald de Dreuzy,et al.  Flow Simulation in Three-Dimensional Discrete Fracture Networks , 2009, SIAM J. Sci. Comput..

[27]  William Gropp,et al.  PETSc Users Manual Revision 3.4 , 2016 .

[28]  Lee Hartley,et al.  Approaches and algorithms for groundwater flow modeling in support of site investigations and safety assessment of the Forsmark site, Sweden , 2013 .

[29]  Stefano Berrone,et al.  The virtual element method for discrete fracture network simulations , 2014 .

[30]  Scott L. Painter,et al.  Calculation of resident groundwater concentration by post-processing particle-tracking results , 2012, Computational Geosciences.

[31]  K. Pruess,et al.  TOUGH2 User's Guide Version 2 , 1999 .

[32]  Georg Kosakowski,et al.  Transport behavior in three‐dimensional fracture intersections , 2003 .

[33]  Carl W. Gable,et al.  Pathline tracing on fully unstructured control-volume grids , 2012, Computational Geosciences.

[34]  Derek Elsworth,et al.  A hybrid boundary element-finite element analysis procedure for fluid flow simulation in fractured rock masses , 1986 .

[35]  Chuan Lu,et al.  PFLOTRAN: Reactive Flow & Transport Code for Use on Laptops to Leadership-Class Supercomputers , 2012 .

[36]  Satish Karra,et al.  PFLOTRAN User Manual A Massively Parallel Reactive Flow and Transport Model for Describing Surface and Subsurface Processes , 2015 .

[37]  Satish Karra,et al.  Shale gas and non-aqueous fracturing fluids: Opportunities and challenges for supercritical CO2 , 2015 .

[38]  Géraldine Pichot,et al.  Influence of fracture scale heterogeneity on the flow properties of three-dimensional discrete fracture networks (DFN) , 2012 .

[39]  G E Hammond,et al.  Evaluating the performance of parallel subsurface simulators: An illustrative example with PFLOTRAN , 2014, Water resources research.

[40]  S. Bachu.  Sequestration of CO2 in geological media in response to climate change: road map for site selection using the transform of the geological space into the CO2 phase space , 2002 .

[41]  Jim Ruppert,et al.  A Delaunay Refinement Algorithm for Quality 2-Dimensional Mesh Generation , 1995, J. Algorithms.

[42]  Jeffrey D. Hyman,et al.  Conforming Delaunay Triangulation of Stochastically Generated Three Dimensional Discrete Fracture Networks: A Feature Rejection Algorithm for Meshing Strategy , 2014, SIAM J. Sci. Comput..

[43]  Peter Jackson,et al.  Multi-scale groundwater flow modeling during temperate climate conditions for the safety assessment of the proposed high-level nuclear waste repository site at Forsmark, Sweden , 2014, Hydrogeology Journal.

[44]  Satish Karra,et al.  Progress Towards Coupled Simulation of Surface/Subsurface Hydrologic Processes and Terrestrial Ecosystem Dynamics Using the Community Models PFLOTRAN and CLM , 2012 .

[45]  C. Fidelibus,et al.  Derivation of equivalent pipe network analogues for three‐dimensional discrete fracture networks by the boundary element method , 1999 .

[46]  Adler,et al.  Fractured Porous Media , 2012 .

[47]  Satish Karra,et al.  Influence of injection mode on transport properties in kilometer‐scale three‐dimensional discrete fracture networks , 2015 .

[48]  Satish Karra,et al.  Particle tracking approach for transport in three-dimensional discrete fracture networks , 2015, Computational Geosciences.

[49]  Stefano Berrone,et al.  A PDE-Constrained Optimization Formulation for Discrete Fracture Network Flows , 2013, SIAM J. Sci. Comput..

[50]  Glenn E. Hammond,et al.  Field‐scale model for the natural attenuation of uranium at the Hanford 300 Area using high‐performance computing , 2010 .

[51]  Nelson L. Max,et al.  A contract based system for large data visualization , 2005, VIS 05. IEEE Visualization, 2005..

[52]  Roussos Dimitrakopoulos,et al.  Discretizing two‐dimensional complex fractured fields for incompressible two‐phase flow , 2011 .

[53]  K. Lipnikov,et al.  On the Reconstruction of Darcy Velocity in Finite-Volume Methods , 2012, Transport in Porous Media.

[54]  Kassem Mustapha,et al.  A New Approach to Simulating Flow in Discrete Fracture Networks with an Optimized Mesh , 2007, SIAM J. Sci. Comput..

[55]  Jean-Raynald de Dreuzy,et al.  Transport and intersection mixing in random fracture networks with power law length distributions , 2001 .