Premium
Matrix Diffusion in Fractured Media: New Insights Into Power Law Scaling of Breakthrough Curves
Author(s) -
Hyman Jeffrey D.,
Rajaram Harihar,
Srinivasan Shriram,
Makedonska Nataliia,
Karra Satish,
Viswanathan Hari,
Srinivasan Gowri
Publication year - 2019
Publication title -
geophysical research letters
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 2.007
H-Index - 273
eISSN - 1944-8007
pISSN - 0094-8276
DOI - 10.1029/2019gl085454
Subject(s) - pareto distribution , power law , advection , scaling , diffusion , statistical physics , matrix (chemical analysis) , continuous time random walk , physics , pareto principle , random walk , distribution (mathematics) , mechanics , mathematics , mathematical analysis , geometry , statistics , materials science , thermodynamics , composite material
We develop a theoretical model for power law tailing behavior of transport in fractured rock based on the relative dominance of the decay rate of the advective travel time distribution, modeled using a Pareto distribution (with tail decaying as ∼ time −(1+ α ) ), versus matrix diffusion, modeled using a Lévy distribution. The theory predicts that when the advective travel time distribution decays sufficiently slowly ( α <1), the late‐time decay rate of the breakthrough curve is −(1+ α /2) rather than the classical −3/2. However, if α >1, the −3/2 decay rate is recovered. For weak matrix diffusion or short advective first breakthrough times, we identify an early‐time regime where the breakthrough curve follows the Pareto distribution, before transitioning to the late‐time decay rate. The theoretical predictions are validated against particle tracking simulations in the three‐dimensional discrete fracture network simulator dfnWorks, where matrix diffusion is incorporated using a time domain random walk.