Solute transport along a single fracture in a porous rock: a simple analytical solution and its extension for modeling velocity dispersion

Abstract

A simple and robust solution is developed for the problem of solute transport along a single fracture in a porous rock. The solution is referred to as the solution to the single-flow-path model and takes the form of a convolution of two functions. The first function is the probability density function of residence-time distribution of a conservative solute in the fracture-only system as if the rock matrix is impermeable. The second function is the response of the fracture-matrix system to the input source when Fickian-type dispersion is completely neglected; thus, the effects of Fickian-type dispersion and matrix diffusion have been decoupled. It is also found that the solution can be understood in a way in line with the concept of velocity dispersion in fractured rocks. The solution is therefore extended into more general cases to also account for velocity variation between the channels. This leads to a development of the multi-channel model followed by detailed statistical descriptions of channel properties and sensitivity analysis of the model upon changes in the model key parameters. The simulation results obtained by the multi-channel model in this study fairly well agree with what is often observed in field experiments—i.e. the unchanged Peclet number with distance, which cannot be predicted by the classical advection-dispersion equation. In light of the findings from the aforementioned analysis, it is suggested that forced-gradient experiments can result in considerably different estimates of dispersivity compared to what can be found in natural-gradient systems for typical channel widths.