Lattice Boltzmann simulation of fluid flow in fracture networks with rough, self-affine surfaces.

Using the lattice Boltzmann method, we study fluid flow in a two-dimensional (2D) model of fracture network of rock. Each fracture in a square network is represented by a 2D channel with rough, self-affine internal surfaces. Various parameters of the model, such as the connectivity and the apertures of the fractures, the roughness profile of their surface, as well as the Reynolds number for flow of the fluid, are systematically varied in order to assess their effect on the effective permeability of the fracture network. The distribution of the fractures' apertures is approximated well by a log-normal distribution, which is consistent with experimental data. Due to the roughness of the fractures' surfaces, and the finite size of the networks that can be used in the simulations, the fracture network is anisotropic. The anisotropy increases as the connectivity of the network decreases and approaches the percolation threshold. The effective permeability K of the network follows the power law K approximately (beta), where is the average aperture of the fractures in the network and the exponent beta may depend on the roughness exponent. A crossover from linear to nonlinear flow regime is obtained at a Reynolds number Re approximately O(1), but the precise numerical value of the crossover Re depends on the roughness of the fractures' surfaces.

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

[2]  Tamaz Chelidze,et al.  Evidence of fractal fracture , 1990 .

[3]  T. Sanderson,et al.  Characteristics of energetic particle events associated with interplanetary shocks , 1985 .

[4]  S. Redner,et al.  Introduction To Percolation Theory , 2018 .

[5]  L. Luo,et al.  Theory of the lattice Boltzmann method: From the Boltzmann equation to the lattice Boltzmann equation , 1997 .

[6]  Reconstruction and geostatistical analysis of multiscale fracture apertures in a large block of welded tuff , 1992 .

[7]  A. Coates,et al.  Determination of comet Halley Gas emission characteristics from mass loading of the solar wind , 1990 .

[8]  N. Odling,et al.  Natural fracture profiles, fractal dimension and joint roughness coefficients , 1994 .

[9]  Mukhopadhyay,et al.  Scaling properties of a percolation model with long-range correlations. , 1996, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[10]  J. Thovert,et al.  Permeability of Self-Affine Fractures , 2001 .

[11]  James M. Keller,et al.  Lattice-Gas Cellular Automata: Inviscid two-dimensional lattice-gas hydrodynamics , 1997 .

[12]  Muhammad Sahimi,et al.  COARSENING OF HETEROGENEOUS MEDIA : APPLICATION OF WAVELETS , 1997 .

[13]  Pierre M. Adler,et al.  PERMEABILITY OF THREE-DIMENSIONAL FRACTURE NETWORKS , 1998 .

[14]  Pierre M. Adler,et al.  PERCOLATION AND CONDUCTIVITY OF SELF-AFFINE FRACTURES , 1999 .

[15]  R. A. Johns,et al.  Nondestructive measurements of fracture aperture in crystalline rock cores using X ray computed tomography , 1993 .

[16]  P. C. Robinson,et al.  Numerical calculations of critical densities for lines and planes , 1984 .

[17]  Z. Jaeger,et al.  Fluid Flow Through a Crack Network in Rocks , 1983 .

[18]  Stephen R. Brown.  A note on the description of surface roughness using fractal dimension , 1987 .

[19]  Muhammad Sahimi,et al.  Multiresolution wavelet coarsening and analysis of transport in heterogeneous media , 2002 .

[20]  Permeability of self-affine rough fractures , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[21]  Eva Hakami,et al.  Aperture measurements and flow experiments on a single natural fracture , 1996 .

[22]  Yannis C. Yortsos,et al.  Percolation in a fractional Brownian motion lattice , 1996 .

[23]  Brian Berkowitz,et al.  Analysis of fracture network connectivity using percolation theory , 1995 .

[24]  P. A. Witherspoon,et al.  Hydromechanical behavior of a deformable rock fracture subject to normal stress , 1981 .

[25]  Alex Hansen,et al.  Study of fracture permeability using Lattice Gas Automata , 1995 .

[26]  Takaji Inamuro,et al.  A NON-SLIP BOUNDARY CONDITION FOR LATTICE BOLTZMANN SIMULATIONS , 1995, comp-gas/9508002.

[27]  Equivalent channel network model for permeability and electrical conductivity of fracture networks , 2002 .

[28]  Pierre M. Adler,et al.  Solute transport at fracture intersections , 2002 .

[29]  Olivier Bour,et al.  Connectivity of random fault networks following a power law fault length distribution , 1997 .

[30]  Brian Berkowitz,et al.  Application of a percolation model to flow in fractured hard rocks , 1991 .

[31]  M. Sahimi.  Flow phenomena in rocks : from continuum models to fractals, percolation, cellular automata, and simulated annealing , 1993 .

[32]  D. T. Snow,et al.  The frequency and apertures of fractures in rock , 1970 .

[33]  B. L. Cox,et al.  Fractal Surfaces: Measurement and Applications in the Earth Sciences , 1993 .

[34]  M. Isichenko.  Percolation, statistical topography, and transport in random media , 1992 .

[35]  Stéphane Roux,et al.  Field measurements of the roughness of fault surfaces , 1993 .

[36]  Olivier Bour,et al.  On the connectivity of three‐dimensional fault networks , 1998 .

[37]  Muhammad Sahimi,et al.  Non-linear and non-local transport processes in heterogeneous media: from long-range correlated percolation to fracture and materials breakdown , 1998 .

[38]  Qisu Zou,et al.  A improved incompressible lattice Boltzmann model for time-independent flows , 1995 .

[39]  Pierre M. Adler,et al.  Fractures and Fracture Networks , 1999 .

[40]  Muhammad Sahimi,et al.  Fluid flow across mass fractals and self-affine surfaces , 1996 .