Origins of anomalous transport in heterogeneous media: Structural and dynamic controls

Abstract

Anomalous (or "non-Fickian") transport is ubiquitous in the context of tracer migration in geological formations. We quantitatively identify the origin of anomalous transport in a representative model of a heterogeneous porous medium under uniform (in the mean) flow conditions; we focus on anomalous transport which arises in the complex flow patterns of lognormally distributed hydraulic conductivity (K) fields, with several decades of K values. Transport in the domains is determined by a particle tracking technique and characterized by breakthrough curves (BTCs). The BTC averaged over multiple realizations demonstrates anomalous transport in all cases, which is accounted for entirely by a power law distribution ∼t-1-β of local transition times. The latter is contained in the probability density function ψ(t) of transition times, embedded in the framework of a continuous time random walk (CTRW). A unique feature of our analysis is the derivation of ψ(t) as a function of parameters quantifying the heterogeneity of the domain. In this context, we first establish the dominance of preferential pathways across each domain, and characterize the statistics of these pathways by forming a particle-visitation weighted histogram, Hw(K), of the hydraulic conductivity. By converting the ln(K) dependence of Hw(K) into time, we demonstrate the equivalence of Hw(K) and ψ(t), and delineate the region of Hw(K) that forms the power law of ψ(t). This thus defines the origin of anomalous transport. Analysis of the preferential pathways clearly demonstrates the limitations of critical path analysis and percolation theory as a basis for determining the origin of anomalous transport. Furthermore, we derive an expression defining the power law exponent β in terms of the Hw(K) parameters. The equivalence between Hw(K) and ψ(t) is a remarkable result, particularly given the nature of the K heterogeneity, the complexity of the flow field within each realization, and the statistics of the particle transitions. Key Points Quantitative connection between CTRW parameters and conductivities is determined Dynamic controls are critical factors to determine key transport features Transport is not explained only by structural knowledge of the disordered medium © 2014. American Geophysical Union. All Rights Reserved.