Inverse Stochastic Moment Analysis of Transient Fow in Randomly Heterogeneous Media

Persistent Link:
http://hdl.handle.net/10150/193932
Title:
Inverse Stochastic Moment Analysis of Transient Fow in Randomly Heterogeneous Media
Author:
Malama, Bwalya
Issue Date:
2005
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:
A geostatistical inverse method of estimating hydraulic parameters of a heterogeneous porous medium at discrete points in space, called pilot points, is presented. In this inverse method the parameter estimation problem is posed as a nonlinear optimization problem with a likelihood based objective function. The likelihood based objective function is expressed in terms of head residuals at head measurement locations in the flow domain, where head residuals are the differences between measured and model-predicted head values. Model predictions of head at each iteration of the optimization problem are obtained by solving a forward problem that is based on nonlocal conditional ensemble mean flow equations. Nonlocal moment equations make possible optimal deterministic predictions of fluid flow in randomly heterogenous porous media as well as assessment of the associated predictive uncertainty. In this work, the nonlocal moment equations are approximated to second order in the standard deviation of log-transformed hydraulic conductivity, and are solved using the finite element method. To enhance computational efficiency, computations are carried out in the complex Laplace-transform space, after which the results are inverted numerically to the real temporal domain for analysis and presentation. Whereas a forward solution can be conditioned on known values of hydraulic parameters, inversion allows further conditioning of the solution on measurements of system state variables, as well as for the estimation of unknown hydraulic parameters. The Levenberg-Marquardt algorithm is used to solve the optimization problem. The inverse method is illustrated through two numerical examples where parameter estimates and the corresponding predictions of system state are conditioned on measurements of head only, and on measurements of head and log-transformed hydraulic conductivity with prior information. An example in which predictions of system state are conditioned only on measurements of log-conductivity is also included for comparison. A fourth example is included in which the estimation of spatially constant specific storage is demonstrated. In all the examples, a superimposed mean uniform and convergent transient flow field through a bounded square domain is used. The examples show that conditioning on measurements of both head and hydraulic parameters with prior information yields more reliable (low uncertainty and good fit) predictions of system state, than when such information is not incorporated into the estimation process.
Type:
text; Electronic Dissertation
Keywords:
Geostatistics; Inverse methods; Stochastic Hydrology
Degree Name:
PhD
Degree Level:
doctoral
Degree Program:
Hydrology; Graduate College
Degree Grantor:
University of Arizona
Committee Chair:
Neuman, Shlomo P.

Full metadata record

DC FieldValue Language
dc.language.isoENen_US
dc.titleInverse Stochastic Moment Analysis of Transient Fow in Randomly Heterogeneous Mediaen_US
dc.creatorMalama, Bwalyaen_US
dc.contributor.authorMalama, Bwalyaen_US
dc.date.issued2005en_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.abstractA geostatistical inverse method of estimating hydraulic parameters of a heterogeneous porous medium at discrete points in space, called pilot points, is presented. In this inverse method the parameter estimation problem is posed as a nonlinear optimization problem with a likelihood based objective function. The likelihood based objective function is expressed in terms of head residuals at head measurement locations in the flow domain, where head residuals are the differences between measured and model-predicted head values. Model predictions of head at each iteration of the optimization problem are obtained by solving a forward problem that is based on nonlocal conditional ensemble mean flow equations. Nonlocal moment equations make possible optimal deterministic predictions of fluid flow in randomly heterogenous porous media as well as assessment of the associated predictive uncertainty. In this work, the nonlocal moment equations are approximated to second order in the standard deviation of log-transformed hydraulic conductivity, and are solved using the finite element method. To enhance computational efficiency, computations are carried out in the complex Laplace-transform space, after which the results are inverted numerically to the real temporal domain for analysis and presentation. Whereas a forward solution can be conditioned on known values of hydraulic parameters, inversion allows further conditioning of the solution on measurements of system state variables, as well as for the estimation of unknown hydraulic parameters. The Levenberg-Marquardt algorithm is used to solve the optimization problem. The inverse method is illustrated through two numerical examples where parameter estimates and the corresponding predictions of system state are conditioned on measurements of head only, and on measurements of head and log-transformed hydraulic conductivity with prior information. An example in which predictions of system state are conditioned only on measurements of log-conductivity is also included for comparison. A fourth example is included in which the estimation of spatially constant specific storage is demonstrated. In all the examples, a superimposed mean uniform and convergent transient flow field through a bounded square domain is used. The examples show that conditioning on measurements of both head and hydraulic parameters with prior information yields more reliable (low uncertainty and good fit) predictions of system state, than when such information is not incorporated into the estimation process.en_US
dc.typetexten_US
dc.typeElectronic Dissertationen_US
dc.subjectGeostatisticsen_US
dc.subjectInverse methodsen_US
dc.subjectStochastic Hydrologyen_US
thesis.degree.namePhDen_US
thesis.degree.leveldoctoralen_US
thesis.degree.disciplineHydrologyen_US
thesis.degree.disciplineGraduate Collegeen_US
thesis.degree.grantorUniversity of Arizonaen_US
dc.contributor.chairNeuman, Shlomo P.en_US
dc.contributor.committeememberYeh, Tian-Chyi J.en_US
dc.contributor.committeememberMyers, Donald E.en_US
dc.contributor.committeememberFerre, Paul A.en_US
dc.identifier.proquest1615en_US
dc.identifier.oclc137356108en_US
All Items in UA Campus Repository are protected by copyright, with all rights reserved, unless otherwise indicated.