Diffusion of Tracer in Altered Tonalite: Experiments and Simulations with Heterogeneous Distribution of Porosity

Numerical time-domain-diffusion simulations were used for studying the diffusion behavior of tracer molecules in rock matrix with homogeneous and heterogeneous porosity. As the heterogeneous sample in these simulations, a 3D tomographic image of altered tonalite was used, in which the mineral components and the pores resolved by X-ray microtomography were represented by their respective intragranular porosities determined previously by the ¹⁴C-PMMA method. The apparent diffusion coefficient of a tracer in altered tonalite was determined experimentally, and was then used in the simulations. In the altered tonalite analyzed, inclusion of heterogeneity in the porosity increased the diffusion coefficient by 16 %. Altered and pristine feldspar was the main mineral component in the sample (72 %), and it also provided the dominant contribution to tracer diffusion, explaining alone 52 % of the diffusion coefficient. The large pores resolved by microtomography (6 %) and altered and pristine mica (22 %) gave an equal contribution to the diffusion coefficient. The simulation method applied was also validated by comparing the results to both an analytical and a numerical solution to the diffusion equation in a homogenous medium. In addition, the method was compared to discrete-time random-walk simulations in the case of randomly placed overlapping spheres.