A computational framework for the solution of infinite-dimensional Bayesian statistical inverse problems with application to global seismic inversion
Quantifying uncertainties in large-scale forward and inverse PDE simulations has emerged as a central challenge facing the field of computational science and engineering. The promise of modeling and simulation for prediction, design, and control cannot be fully realized unless uncertainties in model...
Main Author: | |
---|---|
Other Authors: | |
Format: | Others |
Language: | en |
Published: |
2015
|
Subjects: | |
Online Access: | http://hdl.handle.net/2152/31374 |
id |
ndltd-UTEXAS-oai-repositories.lib.utexas.edu-2152-31374 |
---|---|
record_format |
oai_dc |
spelling |
ndltd-UTEXAS-oai-repositories.lib.utexas.edu-2152-313742015-09-20T17:33:51ZA computational framework for the solution of infinite-dimensional Bayesian statistical inverse problems with application to global seismic inversionMartin, James Robert, Ph. D.Bayesian inferenceInfinite-dimensional inverse problemsUncertainty quantificationLow-rank approximationOptimalityScalable algorithmsHigh performance computingMarkov chain Monte CarloStochastic NewtonSeismic wave propagationQuantifying uncertainties in large-scale forward and inverse PDE simulations has emerged as a central challenge facing the field of computational science and engineering. The promise of modeling and simulation for prediction, design, and control cannot be fully realized unless uncertainties in models are rigorously quantified, since this uncertainty can potentially overwhelm the computed result. While statistical inverse problems can be solved today for smaller models with a handful of uncertain parameters, this task is computationally intractable using contemporary algorithms for complex systems characterized by large-scale simulations and high-dimensional parameter spaces. In this dissertation, I address issues regarding the theoretical formulation, numerical approximation, and algorithms for solution of infinite-dimensional Bayesian statistical inverse problems, and apply the entire framework to a problem in global seismic wave propagation. Classical (deterministic) approaches to solving inverse problems attempt to recover the “best-fit” parameters that match given observation data, as measured in a particular metric. In the statistical inverse problem, we go one step further to return not only a point estimate of the best medium properties, but also a complete statistical description of the uncertain parameters. The result is a posterior probability distribution that describes our state of knowledge after learning from the available data, and provides a complete description of parameter uncertainty. In this dissertation, a computational framework for such problems is described that wraps around the existing forward solvers, as long as they are appropriately equipped, for a given physical problem. Then a collection of tools, insights and numerical methods may be applied to solve the problem, and interrogate the resulting posterior distribution, which describes our final state of knowledge. We demonstrate the framework with numerical examples, including inference of a heterogeneous compressional wavespeed field for a problem in global seismic wave propagation with 10⁶ parameters.Ghattas, Omar N.2015-09-18T16:27:29Z2015-082015-08-10August 20152015-09-18T16:27:29ZThesistextapplication/pdfhttp://hdl.handle.net/2152/31374en |
collection |
NDLTD |
language |
en |
format |
Others
|
sources |
NDLTD |
topic |
Bayesian inference Infinite-dimensional inverse problems Uncertainty quantification Low-rank approximation Optimality Scalable algorithms High performance computing Markov chain Monte Carlo Stochastic Newton Seismic wave propagation |
spellingShingle |
Bayesian inference Infinite-dimensional inverse problems Uncertainty quantification Low-rank approximation Optimality Scalable algorithms High performance computing Markov chain Monte Carlo Stochastic Newton Seismic wave propagation Martin, James Robert, Ph. D. A computational framework for the solution of infinite-dimensional Bayesian statistical inverse problems with application to global seismic inversion |
description |
Quantifying uncertainties in large-scale forward and inverse PDE simulations has emerged as a central challenge facing the field of computational science and engineering. The promise of modeling and simulation for prediction, design, and control cannot be fully realized unless uncertainties in models are rigorously quantified, since this uncertainty can potentially overwhelm the computed result. While statistical inverse problems can be solved today for smaller models with a handful of uncertain parameters, this task is computationally intractable using contemporary algorithms for complex systems characterized by large-scale simulations and high-dimensional parameter spaces. In this dissertation, I address issues regarding the theoretical formulation, numerical approximation, and algorithms for solution of infinite-dimensional Bayesian statistical inverse problems, and apply the entire framework to a problem in global seismic wave propagation. Classical (deterministic) approaches to solving inverse problems attempt to recover the “best-fit” parameters that match given observation data, as measured in a particular metric. In the statistical inverse problem, we go one step further to return not only a point estimate of the best medium properties, but also a complete statistical description of the uncertain parameters. The result is a posterior probability distribution that describes our state of knowledge after learning from the available data, and provides a complete description of parameter uncertainty. In this dissertation, a computational framework for such problems is described that wraps around the existing forward solvers, as long as they are appropriately equipped, for a given physical problem. Then a collection of tools, insights and numerical methods may be applied to solve the problem, and interrogate the resulting posterior distribution, which describes our final state of knowledge. We demonstrate the framework with numerical examples, including inference of a heterogeneous compressional wavespeed field for a problem in global seismic wave propagation with 10⁶ parameters. |
author2 |
Ghattas, Omar N. |
author_facet |
Ghattas, Omar N. Martin, James Robert, Ph. D. |
author |
Martin, James Robert, Ph. D. |
author_sort |
Martin, James Robert, Ph. D. |
title |
A computational framework for the solution of infinite-dimensional Bayesian statistical inverse problems with application to global seismic inversion |
title_short |
A computational framework for the solution of infinite-dimensional Bayesian statistical inverse problems with application to global seismic inversion |
title_full |
A computational framework for the solution of infinite-dimensional Bayesian statistical inverse problems with application to global seismic inversion |
title_fullStr |
A computational framework for the solution of infinite-dimensional Bayesian statistical inverse problems with application to global seismic inversion |
title_full_unstemmed |
A computational framework for the solution of infinite-dimensional Bayesian statistical inverse problems with application to global seismic inversion |
title_sort |
computational framework for the solution of infinite-dimensional bayesian statistical inverse problems with application to global seismic inversion |
publishDate |
2015 |
url |
http://hdl.handle.net/2152/31374 |
work_keys_str_mv |
AT martinjamesrobertphd acomputationalframeworkforthesolutionofinfinitedimensionalbayesianstatisticalinverseproblemswithapplicationtoglobalseismicinversion AT martinjamesrobertphd computationalframeworkforthesolutionofinfinitedimensionalbayesianstatisticalinverseproblemswithapplicationtoglobalseismicinversion |
_version_ |
1716824574550605824 |