Non-Fickian dispersion in porous media explained by heterogeneous microscale matrix diffusion

Mobile-immobile mass transfer is widely used to model non-Fickian dispersion in porous media. Nevertheless, the memory function, implemented in the sink/source term of the transport equation to characterize diffusion in the matrix (i.e., the immobile domain), is rarely measured directly. Therefore, the question can be posed as to whether the memory function is just a practical way of increasing the degrees of freedom for fitting tracer test breakthrough curves or whether it actually models the physics of tracer transport. In this paper we first present a technique to measure the memory function of aquifer samples and then compare the results with the memory function fitted from a set of field-scale tracer tests performed in the same aquifer. The memory function is computed by solving the matrix diffusion equation using a random walk approach. The properties that control diffusion (i.e., mobile-immobile interface and immobile domain cluster shapes, porosity, and tortuosity) are investigated by X-ray microtomography. Once the geometry of the matrix clusters is measured, the shape of the memory function is controlled by the value of the porosity at the percolation threshold and of the tortuosity of the diffusion path. These parameters can be evaluated from microtomographic images. The computed memory function compares well with the memory function deduced from the field-scale tracer tests. We conclude that for the reservoir rock studied here, the atypical non-Fickian dispersion measured from the tracer test is well explained by microscale diffusion processes in the immobile domain. A diffusion-controlled mobile-immobile mass transfer model therefore appears to be valid for this specific case.