TY - JOUR
T1 - Theoretical analysis of non-Gaussian heterogeneity effects on subsurface flow and transport
AU - Riva, Monica
AU - Guadagnini, Alberto
AU - Neuman, Shlomo P.
N1 - Funding Information:
Funding from MIUR (Italian Ministry of Education, University and Research, Water JPI, WaterWorks 2014, Project: WE-NEED- Water NEEDs, availability, quality and sustainability) and from the European Union's Horizon 2020 Research and Innovation program (Project ?Furthering the knowledge Base for Reducing the Environmental Footprint of Shale Gas Development? FRACRISK, grant agreement 640979) is acknowledged. This work was supported in part through a contract between the University of Arizona and Vanderbilt University under the Consortium for Risk Evaluation with Stakeholder Participation (CRESP) III, funded by the U.S. Department of Energy. We are grateful to the Associate Editor, Alberto Bellin, for his insightful comments about our study. All data used in the paper will be retained by the authors for at least 5 years after publication and will be available to the readers upon request.
PY - 2017/4/1
Y1 - 2017/4/1
N2 - Much of the stochastic groundwater literature is devoted to the analysis of flow and transport in Gaussian or multi-Gaussian log hydraulic conductivity (or transmissivity) fields, Y(x)=In K(x) (x being a position vector), characterized by one or (less frequently) a multiplicity of spatial correlation scales. Yet Y and many other variables and their (spatial or temporal) increments, ΔY, are known to be generally non-Gaussian. One common manifestation of non-Gaussianity is that whereas frequency distributions of Y often exhibit mild peaks and light tails, those of increments ∆Y, are generally symmetric with peaks that grow sharper, and tails that become heavier, as separation scale or lag between pairs of Y values decreases. A statistical model that captures these disparate, scale-dependent distributions of Y and ∆Y in a unified and consistent manner has been recently proposed by us. This new “generalized sub-Gaussian (GSG)” model has the form Y(x)=U(x)G(x) where G(x) is (generally, but not necessarily) a multiscale Gaussian random field and U(x) is a nonnegative subordinator independent of G. The purpose of this paper is to explore analytically, in an elementary manner, lead-order effects that non-Gaussian heterogeneity described by the GSG model have on the stochastic description of flow and transport. Recognizing that perturbation expansion of hydraulic conductivity K=ey diverges when Y is sub-Gaussian, we render the expansion convergent by truncating Y's domain of definition. We then demonstrate theoretically and illustrate by way of numerical examples that, as the domain of truncation expands, (a) the variance of truncated Y (denoted by Yt) approaches that of Y and (b) the pdf (and thereby moments) of Yt increments approach those of Y increments and, as a consequence, the variogram of Yt approaches that of Y. This in turn guarantees that perturbing Kt=eyt to second order in σYt (the standard deviation of Yt) yields results which approach those we obtain upon perturbing K=ey to second order in σY even as the corresponding series diverges. Our analysis is rendered mathematically tractable by considering mean-uniform steady state flow in an unbounded, two-dimensional domain of mildly heterogeneous Y with a single-scale function G having an isotropic exponential covariance. Results consist of expressions for (a) lead-order autocovariance and cross-covariance functions of hydraulic head, velocity, and advective particle displacement and (b) analogues of preasymptotic as well as asymptotic Fickian dispersion coefficients. We compare these theoretically and graphically with corresponding expressions developed in the literature for Gaussian Y. We find the former to differ from the latter by a factor K=‹U2›/ ‹U›2 ( ‹› denoting ensemble expectation) and the GSG covariance of longitudinal velocity to contain an additional nugget term depending on this same factor. In the limit as Y becomes Gaussian, k reduces to one and the nugget term drops out.
AB - Much of the stochastic groundwater literature is devoted to the analysis of flow and transport in Gaussian or multi-Gaussian log hydraulic conductivity (or transmissivity) fields, Y(x)=In K(x) (x being a position vector), characterized by one or (less frequently) a multiplicity of spatial correlation scales. Yet Y and many other variables and their (spatial or temporal) increments, ΔY, are known to be generally non-Gaussian. One common manifestation of non-Gaussianity is that whereas frequency distributions of Y often exhibit mild peaks and light tails, those of increments ∆Y, are generally symmetric with peaks that grow sharper, and tails that become heavier, as separation scale or lag between pairs of Y values decreases. A statistical model that captures these disparate, scale-dependent distributions of Y and ∆Y in a unified and consistent manner has been recently proposed by us. This new “generalized sub-Gaussian (GSG)” model has the form Y(x)=U(x)G(x) where G(x) is (generally, but not necessarily) a multiscale Gaussian random field and U(x) is a nonnegative subordinator independent of G. The purpose of this paper is to explore analytically, in an elementary manner, lead-order effects that non-Gaussian heterogeneity described by the GSG model have on the stochastic description of flow and transport. Recognizing that perturbation expansion of hydraulic conductivity K=ey diverges when Y is sub-Gaussian, we render the expansion convergent by truncating Y's domain of definition. We then demonstrate theoretically and illustrate by way of numerical examples that, as the domain of truncation expands, (a) the variance of truncated Y (denoted by Yt) approaches that of Y and (b) the pdf (and thereby moments) of Yt increments approach those of Y increments and, as a consequence, the variogram of Yt approaches that of Y. This in turn guarantees that perturbing Kt=eyt to second order in σYt (the standard deviation of Yt) yields results which approach those we obtain upon perturbing K=ey to second order in σY even as the corresponding series diverges. Our analysis is rendered mathematically tractable by considering mean-uniform steady state flow in an unbounded, two-dimensional domain of mildly heterogeneous Y with a single-scale function G having an isotropic exponential covariance. Results consist of expressions for (a) lead-order autocovariance and cross-covariance functions of hydraulic head, velocity, and advective particle displacement and (b) analogues of preasymptotic as well as asymptotic Fickian dispersion coefficients. We compare these theoretically and graphically with corresponding expressions developed in the literature for Gaussian Y. We find the former to differ from the latter by a factor K=‹U2›/ ‹U›2 ( ‹› denoting ensemble expectation) and the GSG covariance of longitudinal velocity to contain an additional nugget term depending on this same factor. In the limit as Y becomes Gaussian, k reduces to one and the nugget term drops out.
KW - analytical solutions
KW - heterogeneity
KW - non-Gaussian
KW - stochastic
KW - subsurface flow and transport
UR - http://www.scopus.com/inward/record.url?scp=85017520257&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85017520257&partnerID=8YFLogxK
U2 - 10.1002/2016WR019353
DO - 10.1002/2016WR019353
M3 - Article
AN - SCOPUS:85017520257
VL - 53
SP - 2998
EP - 3012
JO - Water Resources Research
JF - Water Resources Research
SN - 0043-1397
IS - 4
ER -