Stochastic inverse problems: Models and metrics
Abstract
In past work, we introduced modelbased inverse methods, and applied them to problems in which the anomaly could be reasonably modeled by simple canonical shapes, such as rectangular solids. In these cases the parameters to be inverted would be length, width and height, as well as the occasional probe liftoff or rotation. We are now developing a formulation that allows more flexibility in modeling complex flaws. The idea consists of expanding the flaw in a sequence of basis functions, and then solving for the expansion coefficients of this sequence, which are modeled as independent random variables, uniformly distributed over their range of values. There are a number of applications of such modeling: 1. Connected cracks and multiple halfmoons, which we have noted in a POD set. Ideally we would like to distinguish connected cracks from one long shallow crack. 2. Cracks of irregular profile and shape which have appeared in cold work holes during bolthole eddycurrent inspection. One side of such cracks is much deeper than other. 3. L or C shaped crack profiles at the surface, examples of which have been seen in bolthole cracks. By formulating problems in a stochastic sense, we are able to leverage the stochasticmore »
 Authors:
 Victor Technologies, LLC, Bloomington, IN 474077706 (United States)
 Computational Tools, Gurnee, IL 60031 (United States)
 Statistical Engineering, Palm Beach Gardens, FL 33418 (United States)
 Air Force Research Laboratory (AFRL/RXCA), Wright Patterson AFB, OH 454337817 (United States)
 Publication Date:
 OSTI Identifier:
 22391239
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: AIP Conference Proceedings; Journal Volume: 1650; Journal Issue: 1; Conference: 41. Annual Review of Progress in Quantitative Nondestructive Evaluation, Boise, ID (United States), 2025 Jul 2014; Other Information: (c) 2015 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; 36 MATERIALS SCIENCE; ALGORITHMS; COMPUTERIZED SIMULATION; CRACKS; EDDY CURRENTS; FLEXIBILITY; INVERSE SCATTERING PROBLEM; MATHEMATICAL MODELS; METRICS; MINIMIZATION; RANDOMNESS; ROTATION; SENSITIVITY; STOCHASTIC PROCESSES; SURFACES
Citation Formats
Sabbagh, Elias H., Sabbagh, Harold A., Murphy, R. Kim, Aldrin, John C., Annis, Charles, and Knopp, Jeremy S.. Stochastic inverse problems: Models and metrics. United States: N. p., 2015.
Web. doi:10.1063/1.4914812.
Sabbagh, Elias H., Sabbagh, Harold A., Murphy, R. Kim, Aldrin, John C., Annis, Charles, & Knopp, Jeremy S.. Stochastic inverse problems: Models and metrics. United States. doi:10.1063/1.4914812.
Sabbagh, Elias H., Sabbagh, Harold A., Murphy, R. Kim, Aldrin, John C., Annis, Charles, and Knopp, Jeremy S.. 2015.
"Stochastic inverse problems: Models and metrics". United States.
doi:10.1063/1.4914812.
@article{osti_22391239,
title = {Stochastic inverse problems: Models and metrics},
author = {Sabbagh, Elias H. and Sabbagh, Harold A. and Murphy, R. Kim and Aldrin, John C. and Annis, Charles and Knopp, Jeremy S.},
abstractNote = {In past work, we introduced modelbased inverse methods, and applied them to problems in which the anomaly could be reasonably modeled by simple canonical shapes, such as rectangular solids. In these cases the parameters to be inverted would be length, width and height, as well as the occasional probe liftoff or rotation. We are now developing a formulation that allows more flexibility in modeling complex flaws. The idea consists of expanding the flaw in a sequence of basis functions, and then solving for the expansion coefficients of this sequence, which are modeled as independent random variables, uniformly distributed over their range of values. There are a number of applications of such modeling: 1. Connected cracks and multiple halfmoons, which we have noted in a POD set. Ideally we would like to distinguish connected cracks from one long shallow crack. 2. Cracks of irregular profile and shape which have appeared in cold work holes during bolthole eddycurrent inspection. One side of such cracks is much deeper than other. 3. L or C shaped crack profiles at the surface, examples of which have been seen in bolthole cracks. By formulating problems in a stochastic sense, we are able to leverage the stochastic global optimization algorithms in NLSE, which is resident in VIC3D®, to answer questions of global minimization and to compute confidence bounds using the sensitivity coefficient that we get from NLSE. We will also address the issue of surrogate functions which are used during the inversion process, and how they contribute to the quality of the estimation of the bounds.},
doi = {10.1063/1.4914812},
journal = {AIP Conference Proceedings},
number = 1,
volume = 1650,
place = {United States},
year = 2015,
month = 3
}

The stochastic collocation (SC) and stochastic Galerkin (SG) methods are two wellestablished and successful approaches for solving general stochastic problems. A recently developed method based on stochastic reduced order models (SROMs) can also be used. Herein we provide a comparison of the three methods for some numerical examples; our evaluation only holds for the examples considered in the paper. The purpose of the comparisons is not to criticize the SC or SG methods, which have proven very useful for a broad range of applications, nor is it to provide overall ratings of these methods as compared to the SROM method.more »Cited by 4

Stochastic spectral methods for efficient Bayesian solution of inverse problems.
We present a reformulation of the Bayesian approach to inverse problems, that seeks to accelerate Bayesian inference by using polynomial chaos (PC) expansions to represent random variables. Evaluation of integrals over the unknown parameter space is recast, more efficiently, as Monte Carlo sampling of the random variables underlying the PC expansion. We evaluate the utility of this technique on a transient diffusion problem arising in contaminant source inversion. The accuracy of posterior estimates is examined with respect to the order of the PC representation, the choice of PC basis, and the decomposition of the support of the prior. The computationalmore » 
Stochastic spectral methods for efficient Bayesian solution of inverse problems
We present a reformulation of the Bayesian approach to inverse problems, that seeks to accelerate Bayesian inference by using polynomial chaos (PC) expansions to represent random variables. Evaluation of integrals over the unknown parameter space is recast, more efficiently, as Monte Carlo sampling of the random variables underlying the PC expansion. We evaluate the utility of this technique on a transient diffusion problem arising in contaminant source inversion. The accuracy of posterior estimates is examined with respect to the order of the PC representation, the choice of PC basis, and the decomposition of the support of the prior. The computationalmore » 
A scalable framework for the solution of stochastic inverse problems using a sparse grid collocation approach
Experimental evidence suggests that the dynamics of many physical phenomena are significantly affected by the underlying uncertainties associated with variations in properties and fluctuations in operating conditions. Recent developments in stochastic analysis have opened the possibility of realistic modeling of such systems in the presence of multiple sources of uncertainties. These advances raise the possibility of solving the corresponding stochastic inverse problem: the problem of designing/estimating the evolution of a system in the presence of multiple sources of uncertainty given limited information. A scalable, parallel methodology for stochastic inverse/design problems is developed in this article. The representation of the underlyingmore » 
Solving geosteering inverse problems by stochastic Hybrid Monte Carlo method
The inverse problems arise in almost all fields of science where the realworld parameters are extracted from a set of measured data. The geosteering inversion plays an essential role in the accurate prediction of oncoming strata as well as a reliable guidance to adjust the borehole position on the fly to reach one or more geological targets. This mathematical treatment is not easy to solve, which requires finding an optimum solution among a large solution space, especially when the problem is nonlinear and nonconvex. Nowadays, a new generation of loggingwhiledrilling (LWD) tools has emerged on the market. The socalled azimuthalmore »