Abstract:Fractured crystalline rocks have been chosen or are under consideration by several countries as host rock formations for deep geological repositories for spent nuclear fuel. In such geological formations, flow and solute transport are mostly controlled by a network of connected natural fractures, each of them being characterised by internal heterogeneity, also denoted as roughness. Fractures are, in turn, subject to variable load caused by various factors, such as the presence of thick ice sheets formed during glaciation periods. Understanding how coupled hydro-mechanical (HM) processes affect flow and transport at the scale of a single natural fracture is crucial for a robust parameterisation of large-scale discrete fracture network models, which are not only used for nuclear waste disposal applications but are also of interest to problems related to geothermics, oil and gas production or groundwater remediation. In this work, we analyse and model an HM experiment carried out in a single natural fracture and use the results of both, the experimental and the modelling work, to get insights into fundamental questions such as the applicability of local cubic law or the effect of normal load on channeling. The initial fracture aperture was obtained from laser scanning of the two fracture surfaces and an equivalent initial aperture was then defined by moving the two fracture surfaces together and comparing the results obtained using a Navier–Stokes based computational fluid dynamics (CFD) model with the experimental flowrate obtained for unloaded conditions. The mechanical effect of the different loading stages was simulated using a high-resolution contact model. The different computed fracture apertures were then used to run groundwater flow simulations using a modified Reynolds equation. The results show that, without correction, local cubic law largely overestimates flowrates. Instead, we show that by explicitly acknowledging the difference between the mechanical aperture and the hydraulic aperture and setting the latter equal to 1/5 of the former, cubic law provides a very reasonable approximation of the experimental flowrates over the entire loading cycle. A positive correlation between fluid flow channeling and normal load is also found.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: How to understand the mechanism of water flowing through natural rough fractures? Specifically, the researchers are concerned with how the changes in fracture apertures affect water flow and solute transport under different normal loading conditions. This problem is crucial for the application of deep geological disposal of radioactive waste, and is also relevant to fields such as geothermal energy utilization, oil and gas extraction, and groundwater remediation. ### Background of the Paper The paper points out that some countries have chosen or are considering using fractured crystalline rocks as host rocks for deep geological disposal of radioactive waste. In such geological formations, water flow and solute transport are mainly controlled by a network of interconnected natural fractures, and each fracture has internal heterogeneity, that is, roughness. These fractures are also affected by variable loads caused by various factors, such as thick ice layers formed during the glacial period. ### Research Objectives In order to more accurately simulate and predict the water flow behavior in these fractures, the researchers conducted a series of experiments and numerical analyses, aiming to explore the following key issues: 1. **Applicability of the Local Cubic Law (LCL)**: The local cubic law is a simplified model that assumes parallel and smooth fracture walls. However, this assumption does not hold in natural rough fractures, so its applicability needs to be verified. 2. **The influence of normal load on water flow channelization**: Study the changes in the distribution and flow paths of water in fractures under different normal loading conditions. ### Overview of Methods The researchers adopted an integrated approach, combining experimental data, fracture surface scanning, a computational fluid dynamics (CFD) model based on the Navier - Stokes equations, an elastic - plastic contact model, and a model based on the modified Reynolds equation. The specific steps include: 1. **Experimental Setup**: Conduct a hydromechanical experiment on a rock sample containing a natural fracture, and record the flow rate and pressure under different loading conditions. 2. **Fracture Surface Scanning**: Use laser scanning technology to obtain detailed topological data of the fracture surface. 3. **CFD Simulation**: Reconstruct the fracture volume based on the fracture surface data and perform numerical simulation using the Navier - Stokes equations. 4. **Mechanical Model**: Use a high - resolution contact model (HRCM) to simulate the changes in fracture apertures under different loading conditions. 5. **Modified Reynolds Equation Model**: Combine the results of the experiment and the mechanical model, and use the modified Reynolds equation to simulate water flow. ### Key Findings 1. **Limitations of the Local Cubic Law**: The unmodified local cubic law significantly overestimates the flow rate. By explicitly considering the difference between the mechanical aperture and the hydraulic aperture and setting the hydraulic aperture to 1/5 of the mechanical aperture, the local cubic law can provide a reasonable approximation of the flow rate. 2. **The influence of normal load on water flow channelization**: It was found that there is a positive correlation between the normal load and water flow channelization. ### Formula Summary - **Navier - Stokes Equations**: \[ \rho \left( \frac{\partial \mathbf{u}}{\partial t}+(\mathbf{u} \cdot \nabla)\mathbf{u} \right)=-\nabla p+\mu \nabla^{2}\mathbf{u} \] \[ \nabla \cdot \mathbf{u} = 0 \] where \(\mathbf{u}\) is the velocity of water, \(\rho\) is the density of water, \(\mu\) is the viscosity of water, and \(p\) is the pressure of water. - **Modified Local Cubic Law**: \[ T=(A \cdot \delta_{m})^{B}\frac{\rho g}{12\mu} \] where \(T\) is the permeability, \(\delta_{m}\) is the mechanical aperture, and \(A\) and \(B\) are constants. - **Modified Reynolds Equation**: \[ \nabla \cdot(\delta_{m}^{B}\nabla h)

Experimental and numerical analysis of flow through a natural rough fracture subject to normal loading

Laboratory Simulation of Flow Through Rough-Walled Microfractures under High Hydraulic Gradient

Numerical Simulation of Unsaturated Flow in Single Fracture Considering the Water Exchange Between Fracture and Matrix

Experimental and numerical study of the geometrical and hydraulic characteristics of a single rock fracture during shear

Quantitative Assessments on Fluid Flow Through Fractures Embedded in Permeable Host Rocks: Experiments and Simulations

Mesoscopic Deformation of a Hydraulic and Mechanical Aperture of a Single Fracture under Normal Stress

Reduced‐Order Models Unravel the Joint Impact of Aperture Heterogeneity and Shear‐Thinning Rheology on Fracture‐Scale Flow Metrics

Modeling the impact of natural roughness of tension joints on heat transport

Channeling flow and anomalous transport due to the complex void structure of rock fractures

Modeling of Flow Characteristics in 3D Rough Rock Fracture with Geometry Changes under Confining Stresses

Experimental investigation on macroscopic behavior and microfluidic field of nonlinear flow in rough-walled artificial fracture models

A Mechanical-Hydraulic-solute Transport Model for Rough-Walled Rock Fractures Subjected to Shear under Constant Normal Stiffness Conditions

Analytical solution of coupled stress-flow-transport processes in a single rock fracture

An Experimental Setup to Investigate Non-Newtonian Fluid Flow in Variable Aperture Channels

Radial Flow Behaviors of a Rough Beishan Granite Fracture under Normal and Thermal Loadings

A model for saturated–unsaturated flow with fractures acting as capillary barriers

The Role of Fracture Surface Roughness in Macroscopic Fluid Flow and Heat Transfer in Fractured Rocks

Effect of shear-induced contact area and aperture variations on nonlinear flow behaviors in fractal rock fractures

Modified Governing Equation and Numerical Simulation of Seepage Flow in a Single Fracture with Three-Dimensional Roughness

A Hydromechanical Model for a Single Rock Fracture Subjected to Shearing

Linear and nonlinear fluid flow responses of connected fractures subject to shearing under constant normal load and constant normal stiffness boundary conditions