Stochastic approach to steady state flow in nonuniform geologic media

Persistent Link:
http://hdl.handle.net/10150/191175
Title:
Stochastic approach to steady state flow in nonuniform geologic media
Author:
Orr, Shlomo.
Issue Date:
1993
Publisher:
The University of Arizona.
Rights:
Copyright © is held by the author. Digital access to this material is made possible by the University Libraries, University of Arizona. Further transmission, reproduction or presentation (such as public display or performance) of protected items is prohibited except with permission of the author.
Abstract:
This dissertation considers the effect of measuring randomly varying local hydraulic conductivities K(x) on one's ability to predict steady state flow within a bounded domain, driven by random source and boundary functions. That is, the work concerns the prediction of local hydraulic head h(x) and Darcy flux q(x) by means of their unbiased ensemble moments (h(x))(κ) and (q(x))(κ) conditioned on measurements of K(x). These predictors satisfy a deterministic flow equation in which (q(x))(κ) = -(κ)(x)∇(h(x))(κ) + r(κ)(x) where κ(x) is a relatively smooth unbiased estimate of K(x) and r(κ)(x) is a "residual flux." A compact integral expression is derived for r(κ)(x) which is rigorously valid for a broad class of K(x) fields, including fractals. It demonstrates that (q(x))(κ) is nonlocal and non-Darcian so that an effective hydraulic conductivity does not generally exist. It is shown analytically that under uniform mean flow the effective conductivity may be a scalar, a symmetric or a nonsymmetric tensor, or a set of directional scalars which do not form a tensor. For cases where r(κ)(x) can neither be expressed nor approximated by a local expression, a weak (integral) approximation (closure) is proposed, which appears to work well in media with pronounced heterogeneity and improves as the quantity and quality of K(x) measurements increase. The nonlocal deterministic flow equation can be solved numerically by standard methods; the theory here shows clearly how the scale of grid discretization should relate to the scale, quantity and quality of available data. After providing explicit approximations for the prediction error moments of head and flux, some practical methods are discussed to compute κ(x) from noisy measurements of K(x) and to calculate required second moments of the associated estimation errors when K(x) is log normal. Nonuniform mean flow is studied by conducting high resolution Monte Carlo simulations of two dimensional seepage to a point sink in statistically homogeneous and isotropic log normal K(x) fields. These reveal the existence of radial effective hydraulic conductivity which increases from the harmonic mean of K(x) near interior and boundary sources to geometric mean far from such sources for σ^2/Υ (the variance of ln K) at least as large as 4. They suggest the possibility of replacing r(κ)(x) by a local expression at distances of few conditional integral scales from the interior and boundary sources. Special attention is paid to the "art" of random field generation, and comparisons are made between four alternative methods with five different random number generators.
Type:
Dissertation-Reproduction (electronic); text
Keywords:
Hydrology.; Groundwater flow -- Mathematical models.; Monte Carlo method.
Degree Name:
Ph. D.
Degree Level:
doctoral
Degree Program:
Hydrology and Water Resources; Graduate College
Degree Grantor:
University of Arizona
Committee Chair:
Neuman, Shlomo P.

Full metadata record

DC FieldValue Language
dc.language.isoenen_US
dc.titleStochastic approach to steady state flow in nonuniform geologic mediaen_US
dc.creatorOrr, Shlomo.en_US
dc.contributor.authorOrr, Shlomo.en_US
dc.date.issued1993en_US
dc.publisherThe University of Arizona.en_US
dc.rightsCopyright © is held by the author. Digital access to this material is made possible by the University Libraries, University of Arizona. Further transmission, reproduction or presentation (such as public display or performance) of protected items is prohibited except with permission of the author.en_US
dc.description.abstractThis dissertation considers the effect of measuring randomly varying local hydraulic conductivities K(x) on one's ability to predict steady state flow within a bounded domain, driven by random source and boundary functions. That is, the work concerns the prediction of local hydraulic head h(x) and Darcy flux q(x) by means of their unbiased ensemble moments (h(x))(κ) and (q(x))(κ) conditioned on measurements of K(x). These predictors satisfy a deterministic flow equation in which (q(x))(κ) = -(κ)(x)∇(h(x))(κ) + r(κ)(x) where κ(x) is a relatively smooth unbiased estimate of K(x) and r(κ)(x) is a "residual flux." A compact integral expression is derived for r(κ)(x) which is rigorously valid for a broad class of K(x) fields, including fractals. It demonstrates that (q(x))(κ) is nonlocal and non-Darcian so that an effective hydraulic conductivity does not generally exist. It is shown analytically that under uniform mean flow the effective conductivity may be a scalar, a symmetric or a nonsymmetric tensor, or a set of directional scalars which do not form a tensor. For cases where r(κ)(x) can neither be expressed nor approximated by a local expression, a weak (integral) approximation (closure) is proposed, which appears to work well in media with pronounced heterogeneity and improves as the quantity and quality of K(x) measurements increase. The nonlocal deterministic flow equation can be solved numerically by standard methods; the theory here shows clearly how the scale of grid discretization should relate to the scale, quantity and quality of available data. After providing explicit approximations for the prediction error moments of head and flux, some practical methods are discussed to compute κ(x) from noisy measurements of K(x) and to calculate required second moments of the associated estimation errors when K(x) is log normal. Nonuniform mean flow is studied by conducting high resolution Monte Carlo simulations of two dimensional seepage to a point sink in statistically homogeneous and isotropic log normal K(x) fields. These reveal the existence of radial effective hydraulic conductivity which increases from the harmonic mean of K(x) near interior and boundary sources to geometric mean far from such sources for σ^2/Υ (the variance of ln K) at least as large as 4. They suggest the possibility of replacing r(κ)(x) by a local expression at distances of few conditional integral scales from the interior and boundary sources. Special attention is paid to the "art" of random field generation, and comparisons are made between four alternative methods with five different random number generators.en_US
dc.description.notehydrology collectionen_US
dc.typeDissertation-Reproduction (electronic)en_US
dc.typetexten_US
dc.subjectHydrology.en_US
dc.subjectGroundwater flow -- Mathematical models.en_US
dc.subjectMonte Carlo method.en_US
thesis.degree.namePh. D.en_US
thesis.degree.leveldoctoralen_US
thesis.degree.disciplineHydrology and Water Resourcesen_US
thesis.degree.disciplineGraduate Collegeen_US
thesis.degree.grantorUniversity of Arizonaen_US
dc.contributor.chairNeuman, Shlomo P.en_US
dc.contributor.committeememberMaddock, Thomasen_US
dc.contributor.committeememberYeh, Jim T-Cen_US
dc.contributor.committeememberWarrick, Art W.en_US
dc.contributor.committeememberWierrenga, Peter J.en_US
dc.identifier.oclc226969115en_US
All Items in UA Campus Repository are protected by copyright, with all rights reserved, unless otherwise indicated.