Significance of higher order moments to the completecharacterization of the travel time pdf in heterogeneous porous media using the maximum entropyprinciple
2010 (English)In: Water resources research, ISSN 0043-1397, E-ISSN 1944-7973, Vol. 46, no W05502Article in journal (Refereed) Published
The travel time formulation of advective transport in heterogeneous porous media is of interest both conceptually, e. g., for incorporating retention processes, and in applications where typically the travel time peak, early, and late arrivals of contaminants are of major concern in a regulatory or remediation context. Furthermore, the travel time moments are of interest for quantifying uncertainty in advective transport of tracers released from point sources in heterogeneous aquifers. In view of this interest, the travel time distribution has been studied in the literature; however, the link to the hydraulic conductivity statistics has been typically restricted to the first two moments. Here we investigate the influence of higher travel time moments on the travel time probability density function (pdf) in heterogeneous porous media combining Monte Carlo simulations with the maximum entropy principle. The Monte Carlo experimental pdf is obtained by the adaptive Fup Monte Carlo method (AFMCM) for advective transport characterized by a multi-Gaussian structure with exponential covariance considering two injection modes (in-flux and resident) and lnK variance up to 8. A maximum entropy (MaxEnt) algorithm based on Fup basis functions is used for the complete characterization of the travel time pdf. All travel time moments become linear with distance. Initial nonlinearity is found mainly for the resident injection mode, which exhibits a strong nonlinearity within first 5I(Y) for high heterogeneity. For the resident injection mode, the form of variance and all higher moments changes from the familiar concave form predicted by the first-order theory to a convex form; for the in-flux mode, linearity is preserved even for high heterogeneity. The number of moments sufficient for a complete characterization of the travel time pdf mainly depends on the heterogeneity level. Mean and variance completely describe travel time pdf for low and mild heterogeneity, skewness is dominant for lnK variance around 4, while kurtosis and fifth moment are required for lnK variance higher than 4. Including skewness seems sufficient for describing the peak and late arrivals. Linearity of travel time moments enables the prediction of asymptotic behavior of the travel time pdf which in the limit converges to a symmetric distribution and Fickian transport. However, higher-order travel time moments may be important for most practical purposes and in particular for advective transport in highly heterogeneous porous media for a long distance from the source.
Place, publisher, year, edition, pages
2010. Vol. 46, no W05502
FULL-BAYESIAN APPROACH; SOLUTE FLUX APPROACH; MEAN UNIFORM-FLOW; GROUNDWATER-FLOW; INVERSE PROBLEM; CONSERVATIVE SOLUTES; MASS-TRANSFER; TRANSPORT; INFORMATION; STATISTICS
IdentifiersURN: urn:nbn:se:kth:diva-10651DOI: 10.1029/2009WR008220ISI: 000277266300007ScopusID: 2-s2.0-77952264451OAI: oai:DiVA.org:kth-10651DiVA: diva2:222646
QC 201007142009-06-092009-06-092011-01-03Bibliographically approved