Structural controls on anomalous transport in fractured porous rock

Anomalous transport is ubiquitous in a wide range of disordered systems, notably in fractured porous formations. We quantitatively identify the structural controls on anomalous tracer transport in a model of a real fractured geological formation that was mapped in an outcrop. The transport, determined by a continuum scale mathematical model, is characterized by breakthrough curves (BTCs) that document anomalous (or “non-Fickian”) transport, which is accounted for by a power-law distribution of local transition times ψ(t) within the framework of a continuous time random walk (CTRW). We show that the determination of ψ(t) is related to fractures aligned approximately with the macroscopic direction of flow. We establish the dominant role of fracture alignment, and assess the statistics of these fractures by determining a concentration-visitation weighted residence time histogram. We then convert the histogram to a probability density function (pdf) that coincides with the CTRW ψ(t) and hence anomalous transport. We show that the permeability of the geological formation hosting the fracture network has a limited effect on the anomalous nature of the transport; rather, it is the fractures transverse to the flow direction that play the major role in forming the long BTC tail associated with anomalous transport. This is a remarkable result, given the complexity of the flow field statistics as captured by concentration transitions. This article is protected by copyright. All rights reserved.

[1]  R. Helmig,et al.  Black-Oil Simulations for Three-Component -- Three-Phase Flow in Fractured Porous Media , 2007 .

[2]  M. Dentz,et al.  Modeling non‐Fickian transport in geological formations as a continuous time random walk , 2006 .

[3]  Martha Lien,et al.  A 3D Computational Study of Effective Medium Methods Applied to Fractured Media , 2013, Transport in Porous Media.

[4]  S. Geiger,et al.  Numerical simulation of water injection into layered fractured carbonate reservoir analogs , 2006 .

[5]  K. S. Schmid,et al.  Higher order FE-FV method on unstructured grids for transport and two-phase flow with variable viscosity in heterogeneous porous media , 2013, J. Comput. Phys..

[6]  Sebastian Geiger,et al.  A Novel Multi-rate Dual-porosity Model for Improved Simulation of Fractured and Multi-porosity Reservoirs , 2011 .

[7]  Brian Berkowitz,et al.  Computing “Anomalous” Contaminant Transport in Porous Media: The CTRW MATLAB Toolbox , 2005, Ground water.

[8]  Sean Andrew McKenna,et al.  Tracer tests in a fractured dolomite: 2. Analysis of mass transfer in single‐well injection‐withdrawal tests , 1999 .

[9]  J. Long,et al.  From field data to fracture network modeling: An example incorporating spatial structure , 1987 .

[10]  Alberto Guadagnini,et al.  Origins of anomalous transport in heterogeneous media: Structural and dynamic controls , 2014 .

[11]  F. Schwartz,et al.  An Analysis of the Influence of Fracture Geometry on Mass Transport in Fractured Media , 1984 .

[12]  Brian Berkowitz,et al.  ANOMALOUS TRANSPORT IN RANDOM FRACTURE NETWORKS , 1997 .

[13]  S. Geiger,et al.  Non‐Fourier thermal transport in fractured geological media , 2010 .

[14]  D. Coumou,et al.  Hydrothermal, multiphase convection of H2O‐NaCl fluids from ambient to magmatic temperatures: a new numerical scheme and benchmarks for code comparison , 2014 .

[15]  M. Dentz,et al.  Impact of velocity correlation and distribution on transport in fractured media: Field evidence and theoretical model , 2015 .

[16]  Donald M. Reeves,et al.  Transport of conservative solutes in simulated fracture networks: 2. Ensemble solute transport and the correspondence to operator‐stable limit distributions , 2008 .

[17]  Sean A. McKenna,et al.  Tracer tests in a fractured dolomite: 3. Double‐porosity, multiple‐rate mass transfer processes in convergent flow tracer tests , 2001 .

[18]  K. Stüben A review of algebraic multigrid , 2001 .

[19]  Ruben Juanes,et al.  Predictability of anomalous transport on lattice networks with quenched disorder. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[20]  Ruben Juanes,et al.  Anomalous transport on regular fracture networks: Impact of conductivity heterogeneity and mixing at fracture intersections. , 2015, Physical review. E, Statistical, nonlinear, and soft matter physics.

[21]  Kenneth Stuart Sorbie,et al.  Semianalytical solutions for cocurrent and countercurrent imbibition and dispersion of solutes in immiscible two‐phase flow , 2011 .

[22]  G. Bertotti,et al.  Fracture-network analysis of the Latemar Platform (northern Italy): integrating outcrop studies to constrain the hydraulic properties of fractures in reservoir models , 2014 .

[23]  K. Bisdom,et al.  The impact of different aperture distribution models and critical stress criteria on equivalent permeability in fractured rocks , 2016 .

[24]  P. King,et al.  Prediction of vein connectivity using the percolation approach: model test with field data , 2006 .

[25]  D. Benson,et al.  Application of a fractional advection‐dispersion equation , 2000 .

[26]  S. Gorelick,et al.  Multiple‐Rate Mass Transfer for Modeling Diffusion and Surface Reactions in Media with Pore‐Scale Heterogeneity , 1995 .

[27]  Georg Kosakowski,et al.  Anomalous transport of colloids and solutes in a shear zone. , 2004, Journal of contaminant hydrology.

[28]  Jean-Raynald de Dreuzy,et al.  Particle-tracking simulations of anomalous transport in hierarchically fractured rocks , 2013, Comput. Geosci..

[29]  S. Geiger,et al.  Multiscale fracture network characterization and impact on flow: A case study on the Latemar carbonate platform , 2015 .

[30]  Kenneth Stuart Sorbie,et al.  Analytical solutions for co- and counter-current imbibition of sorbing, dispersive solutes in immiscible two-phase flow , 2012, Computational Geosciences.

[31]  F. Afșar,et al.  How facies and diagenesis affect fracturing of limestone beds and reservoir permeability in limestone–marl alternations , 2014 .

[32]  X. Sanchez‐Vila,et al.  On matrix diffusion: formulations, solution methods and qualitative effects , 1998 .

[33]  Sebastian Geiger,et al.  A Novel Multi-Rate Dual-Porosity Model for Improved Simulation of Fractured and Multiporosity Reservoirs , 2013 .

[34]  X. Sanchez‐Vila,et al.  On the formation of multiple local peaks in breakthrough curves , 2015 .

[35]  M. Belayneh Palaeostress orientation inferred from surface morphology of joints on the southern margin of the Bristol Channel Basin, UK , 2004, Geological Society, London, Special Publications.

[36]  Martin J. Blunt,et al.  Streamline‐based dual‐porosity simulation of reactive transport and flow in fractured reservoirs , 2004 .

[37]  Vladimir Cvetkovic,et al.  Upscaling particle transport in discrete fracture networks: 2. Reactive tracers , 2007 .

[38]  M. Belayneh,et al.  Finite Element - Node-Centered Finite-Volume Two-Phase-Flow Experiments With Fractured Rock Represented by Unstructured Hybrid-Element Meshes , 2007 .

[39]  Jens T. Birkholzer,et al.  Continuous time random walk analysis of solute transport in fractured porous media , 2008 .

[40]  T. Le Borgne,et al.  Inferring transport characteristics in a fractured rock aquifer by combining single‐hole ground‐penetrating radar reflection monitoring and tracer test data , 2012, 1701.01877.

[41]  S. Geiger,et al.  Upscaling solute transport in naturally fractured porous media with the continuous time random walk method , 2009 .

[42]  B. Berkowitz,et al.  Measurement and analysis of non-Fickian dispersion in heterogeneous porous media. , 2003, Journal of contaminant hydrology.

[43]  K. St A review of algebraic multigrid , 2001 .

[44]  Stefan M. Luthi,et al.  Multi-scale fracture network analysis from an outcrop analogue: A case study from the Cambro-Ordovician clastic succession in Petra, Jordan , 2012 .

[45]  Igor M. Sokolov,et al.  ANOMALOUS TRANSPORT IN EXTERNAL FIELDS : CONTINUOUS TIME RANDOM WALKS AND FRACTIONAL DIFFUSION EQUATIONS EXTENDED , 1998 .

[46]  M. Belayneh,et al.  Fluid flow partitioning between fractures and a permeable rock matrix , 2004 .