Monte Carlo Simulations of Shear‐thinning Flow in Geological Fractures

Alessandro Lenci (1, 2) , Yves Méheust (2) , M. Putti (3) , V. Di Federico (1)

1 UNIBO - Alma Mater Studiorum Università di Bologna = University of Bologna
2 GR - Géosciences Rennes
3 Unipd - Università degli Studi di Padova = University of Padua

The hydraulic behaviour of fractured rocks under shear-thinning flow is a challenging topic of interest in several fields, related either to environmental remediation or to natural resources recovery. The compound effect of fluid rheology and medium heterogeneity strongly affects flow and transport in fractured geological formations. Here, a stochastic analysis is conducted via Monte Carlo simulations to investigate shear-thinning fluids behaviour in fractures subjected to natural and forced flow, considering different fracture dimensions, for a spatial correlation of the fracture that is an intrinsic parameter of the geological formation, independent of the fracture size. Under the lubrication approximation, a generalized Reynolds equation for shear-thinning fluids is solved using an ad hoc, finite volume-based, numerical scheme. The influence of the rheology and aperture field heterogeneity on ensemble statistics of the velocity components and magnitude, as well as apparent fracture-scale transmissivity, is quantified over 103 fracture realizations. The probability density functions (PDFs) obtained by averaging over the set of realizations and the relative confidence intervals are analysed to comprehend the apparent transmissivity transition from Newtonian to shear-thinning regime. The autocorrelation functions of velocity components are computed to understand the impact of rheology on spatial correlations of the flow. Velocity components exhibit narrow PDFs with nearly exponential decay. More elevated pressure gradients emphasize the shear-thinning behaviour, inducing a more pronounced flow localization, under otherwise identical conditions. This translates at the scale of the fracture into a larger apparent transmissivity as compared to the same configuration with Newtonian rheology, by orders of magnitude.

Mots clés: Monte Carlo Simulations; Geological Fractures; Shear-thinning fluids; Lubrication Theory; Stokes Flow; Stochastic Analysis

Domaines: Mécanique des fluides [physics.class-ph]

Alessandro Lenci, Yves Méheust, M. Putti, V. Di Federico. Monte Carlo Simulations of Shear‐thinning Flow in Geological Fractures. Water Resources Research, 2022, 58 (9), pp.e2022WR032024. ⟨10.1029/2022wr032024⟩. ⟨insu-03760349v2⟩