Transition from non-Fickian to Fickian longitudinal transport through 3-D rough fractures: Scale-(in)sensitivity and roughness dependence.

Understanding transport in rough fractures from non-Fickian to Fickian regimes and the prediction of non-Fickian transport is critical for the development of new transport theories and many practical applications. Through computational experiments that fall within the macrodispersion regime, we first simulated and analyzed solute transport through synthetic rough fractures with stationary geometrical properties (i.e., fracture roughness σb/〈b〉 and correlation length λ, where b refers to aperture with its standard deviation σb and arithmetic mean 〈b〉) across increasing fracture longitudinal transport domain length L, with L/λ ranging from 2.5 to 50. The results were used to determine how solute transport behavior evolves with increasing scale in the longitudinal direction. Moreover, a set of correlated fractures with aperture fields following normal and log-normal distributions was created to further identify and quantify the dependence of non-Fickian transport on roughness. We found that although persistent intermittent velocity structures were present, the breakthrough curves (BTCs) and residence time distributions showed diminishing early arrival and tailing, features of non-Fickian transport, with increasing longitudinal L/λ, ultimately converging to a Fickian transport regime given σb/〈b〉 remained constant. Inverse analysis of the experimental BTCs with the advection-dispersion equation (ADE) model showed that the dispersion coefficient (D) was non-trivially scale-dependent. Simulation results for rough fractures with varying σb/〈b〉 and L/λ indicated that the ratio of fluid velocity to transport velocity fitted to the ADE model depends on σb/〈b〉 and L/λ. The continuous time random walk (CTRW) performed much better across all transport scales, and resulted in scale-independent fitted parameters, i.e., β in the memory function. The fitted β is proportional to σb/〈b〉but is insensitive to L/λ. Therefore, bulk longitudinal solute transport across the pre-asymptotic and asymptotic regimes can be estimated based on the CTRW model parameterized by measurable fracture physical properties.

[1]  Stephen R. Brown,et al.  Fluid flow through rock joints: The effect of surface roughness , 1987 .

[2]  Yiping Guo,et al.  Influence of aperture field heterogeneity and anisotropy on dispersion regimes and dispersivity in single fractures , 2009 .

[3]  Daniel M. Tartakovsky,et al.  Perspective on theories of non-Fickian transport in heterogeneous media , 2009 .

[4]  M. Cardenas,et al.  Non‐Fickian transport through two‐dimensional rough fractures: Assessment and prediction , 2014 .

[5]  Richard A. Ketcham,et al.  Modification of the Local Cubic Law of fracture flow for weak inertia, tortuosity, and roughness , 2015 .

[6]  Christopher C. Pain,et al.  Non-linear regimes of fluid flow in rock fractures , 2004 .

[7]  B. Berkowitz,et al.  Interpretation and nonuniqueness of CTRW transition distributions: Insights from an alternative solute transport formulation , 2014 .

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

[9]  Pierre M. Adler,et al.  Permeability of a Single Fracture; Validity of the Reynolds Equation , 1995 .

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

[11]  Xiaoxian Zhang,et al.  Persistence of anomalous dispersion in uniform porous media demonstrated by pore‐scale simulations , 2007 .

[12]  Tanguy Le Borgne,et al.  Modeling preasymptotic transport in flows with significant inertial and trapping effects – The importance of velocity correlations and a spatial Markov model , 2014 .

[13]  B. Berkowitz Characterizing flow and transport in fractured geological media: A review , 2002 .

[14]  Paul W. J. Glover,et al.  Fluid flow through rough fractures in rocks. II: A new matching model for rough rock fractures , 2006 .

[15]  Li Li,et al.  Solute transport in low‐heterogeneity sandboxes: The role of correlation length and permeability variance , 2014 .

[16]  Brian Berkowitz,et al.  Anomalous transport in correlated velocity fields. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[17]  Andreas Englert,et al.  Mixing, spreading and reaction in heterogeneous media: a brief review. , 2011, Journal of contaminant hydrology.

[18]  Seung Hyun Lee,et al.  Tail shortening with developing eddies in a rough‐walled rock fracture , 2015 .

[19]  R. Glass,et al.  Solute transport in variable‐aperture fractures: An investigation of the relative importance of Taylor dispersion and macrodispersion , 2000 .

[20]  Ruben Juanes,et al.  Pore‐scale intermittent velocity structure underpinning anomalous transport through 3‐D porous media , 2014 .

[21]  Tanguy Le Borgne,et al.  Lagrangian statistical model for transport in highly heterogeneous velocity fields. , 2008, Physical review letters.

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

[23]  N. Spycher,et al.  Impact of fluid-rock chemical interactions on tracer transport in fractured rocks. , 2013, Journal of contaminant hydrology.

[24]  Karen Mair,et al.  Roughness of fault surfaces over nine decades of length scales , 2012 .

[25]  Brian Berkowitz,et al.  Time behavior of solute transport in heterogeneous media: transition from anomalous to normal transport , 2003 .

[26]  Clint Dawson,et al.  Upscaling transport of a reacting solute through a peridocially converging-diverging channel at pre-asymptotic times. , 2015, Journal of contaminant hydrology.

[27]  Arturo A. Keller,et al.  Effect of fracture aperture variations on the dispersion of contaminants , 1999 .

[28]  Y. Tsang Usage of “Equivalent apertures” for rock fractures as derived from hydraulic and tracer tests , 1992 .

[29]  Stephen E. Silliman,et al.  An interpretation of the difference between aperture estimates derived from hydraulic and tracer tests in a single fracture , 1989 .

[30]  C. Chrysikopoulos,et al.  Contaminant transport in a fracture with spatially variable aperture in the presence of monodisperse and polydisperse colloids , 2003 .

[31]  Daniel M. Tartakovsky,et al.  Semi‐analytical solutions for solute transport and exchange in fractured porous media , 2012 .

[32]  D. Bolster,et al.  Substrate size and heterogeneity control anomalous transport in small streams , 2014 .

[33]  Ruben Juanes,et al.  Spatial Markov model of anomalous transport through random lattice networks. , 2011, Physical review letters.

[34]  Brian Berkowitz,et al.  Numerical simulation of non‐Fickian transport in geological formations with multiple‐scale heterogeneities , 2002 .

[35]  Olivier Bour,et al.  Persistence of incomplete mixing: a key to anomalous transport. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[36]  M. Fourar,et al.  Non-Fickian dispersion in a single fracture. , 2008, Journal of contaminant hydrology.

[37]  R. Glass,et al.  Experimental and simulated solute transport in a partially‐saturated, variable‐aperture fracture , 2002 .

[38]  Haiyan Zhou,et al.  Steady-state saturated groundwater flow modeling with full tensor conductivities using finite differences , 2010, Comput. Geosci..

[39]  C. Chrysikopoulos,et al.  Transport of polydisperse colloids in a saturated fracture with spatially variable aperture , 2000 .

[40]  Solute transport in periodical heterogeneous porous media: Importance of observation scale and experimental sampling , 2015 .

[41]  M. Cardenas,et al.  An efficient quasi-3D particle tracking-based approach for transport through fractures with application to dynamic dispersion calculation. , 2015, Journal of contaminant hydrology.

[42]  Brian Berkowitz,et al.  Non-Fickian Transport in Transparent Replicas of Rough-Walled Rock Fractures , 2013, Transport in Porous Media.

[43]  P. Bennett,et al.  Theory for dynamic longitudinal dispersion in fractures and rivers with Poiseuille flow , 2012 .

[44]  Brian Berkowitz,et al.  Theory of anomalous chemical transport in random fracture networks , 1998 .

[45]  M. Cardenas Direct simulation of pore level Fickian dispersion scale for transport through dense cubic packed spheres with vortices , 2009 .

[46]  David J. Brush,et al.  Fluid flow in synthetic rough‐walled fractures: Navier‐Stokes, Stokes, and local cubic law simulations , 2003 .

[47]  Richard A. Ketcham,et al.  Effects of inertia and directionality on flow and transport in a rough asymmetric fracture , 2009 .

[48]  B. Berkowitz,et al.  Anomalous Transport in “Classical” Soil and Sand Columns , 2004, Soil Science Society of America Journal.

[49]  D. Bolster,et al.  Natural Organic Matter Transport Modeling with a Continuous Time Random Walk Approach. , 2014, Environmental engineering science.

[50]  Yi‐Feng Chen,et al.  The Friction Factor in the Forchheimer Equation for Rock Fractures , 2016, Rock Mechanics and Rock Engineering.

[51]  D. Bolster,et al.  The impact of inertial effects on solute dispersion in a channel with periodically varying aperture , 2012 .

[52]  Yiping Guo,et al.  On the appropriate “equivalent aperture” for the description of solute transport in single fractures: Laboratory‐scale experiments , 2008 .

[53]  G. Taylor Dispersion of soluble matter in solvent flowing slowly through a tube , 1953, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[54]  Pierre M. Adler,et al.  Taylor dispersion in porous media. Determination of the dispersion tensor , 1993 .

[55]  Alexandre M Tartakovsky,et al.  Flow intermittency, dispersion, and correlated continuous time random walks in porous media. , 2013, Physical review letters.

[56]  Richard A. Ketcham,et al.  Navier‐Stokes flow and transport simulations using real fractures shows heavy tailing due to eddies , 2007 .

[57]  John F. Brady,et al.  A non-local description of advection-diffusion with application to dispersion in porous media , 1987, Journal of Fluid Mechanics.