American Geophysical Union, Water Resources Research, 10(51), p. 8273-8293, 2015
DOI: 10.1002/2015wr017645
Full text: Unavailable
Both flow field heterogeneity and mass transfer between mobile and immobile domains have been studied separately for explaining observed anomalous transport. Here we investigate non-Fickian transport using high-resolution 3-D X-ray microtomographic images of Berea sandstone containing microporous cement with pore size below the setup resolution. Transport is computed for a set of representative elementary volumes and results from advection and diffusion in the resolved macroporosity (mobile domain) and diffusion in the microporous phase (immobile domain) where the effective diffusion coefficient is calculated from the measured local porosity using a phenomenological model that includes a porosity threshold (φθ) below which diffusion is null and the exponent n that characterizes tortuosity-porosity power-law relationship. We show that both flow field heterogeneity and microporosity trigger anomalous transport. Breakthrough curve (BTC) tailing is positively correlated to microporosity volume and mobile-immobile interface area. The sensitivity analysis showed that the BTC tailing increases with the value of φθ, due to the increase of the diffusion path tortuosity until the volume of the microporosity becomes negligible. Furthermore, increasing the value of n leads to an increase in the standard deviation of the distribution of effective diffusion coefficients, which in turn results in an increase of the BTC tailing. Finally, we propose a continuous time random walk upscaled model where the transition time is the sum of independently distributed random variables characterized by specific distributions. It allows modeling a 1-D equivalent macroscopic transport honoring both the control of the flow field heterogeneity and the multirate mass transfer between mobile and immobile domains. © 2015. American Geophysical Union. All Rights Reserved. ; The authors wish to thank PaulTafforeau for his help in the X-raymicrotomography data acquisition atthe European Synchrotron Radiationfacility (ESRF, France). All data used inthis work are available upon from thecorresponding author, while detailsand codes for OpenFOAM are availableat http://openfoam.org/. MDacknowledges the support of theEuropean Research Council (ERC)through the project MHetScale(617511). ; Peer reviewed