Initial Report on the Development of a Monte CarloMarkov Chain Joint Inversion Approach for Geothermal Exploration
Abstract
Geothermal exploration and subsequent characterization of potential resources typically employ a variety of geophysical, geologic and geochemical techniques. However, since the data collected by each technique provide information directly on only one or a very limited set of the many physical parameters that characterize a geothermal system, no single method can be used to describe the system in its entirety. Presently, the usual approach to analyzing disparate data streams for geothermal applications is to invert (or forward model) each data set separately and then combine or compare the resulting models, for the most part in a more or less ad hoc manner. However, while each inversion may yield a model that fits the individual data set, the models are usually inconsistent with each other to some degree. This reflects uncertainties arising from the inevitable fact that geophysical and other exploration data in general are to some extent noisy, incomplete, and of limited sensitivity and resolution, and so yield nonunique results. The purpose of the project described here is to integrate the different model constraints provided by disparate geophysical, geological and geochemical data in a rigorous and consistent manner by formal joint inversion. The objective is to improve the fidelity ofmore »
 Authors:
 Publication Date:
 Research Org.:
 Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
 Sponsoring Org.:
 USDOE
 OSTI Identifier:
 920472
 Report Number(s):
 UCRLTR230755
TRN: US200818%%91
 DOE Contract Number:
 W7405ENG48
 Resource Type:
 Technical Report
 Country of Publication:
 United States
 Language:
 English
 Subject:
 58 GEOSCIENCES; 15 GEOTHERMAL ENERGY; DRILLING; EXPLORATION; GEOTHERMAL EXPLORATION; GEOTHERMAL SYSTEMS; RESOLUTION; SENSITIVITY; TARGETS
Citation Formats
Foxall, W, Ramirez, A, Carlson, S, Dyer, K, and Sun, Y. Initial Report on the Development of a Monte CarloMarkov Chain Joint Inversion Approach for Geothermal Exploration. United States: N. p., 2007.
Web. doi:10.2172/920472.
Foxall, W, Ramirez, A, Carlson, S, Dyer, K, & Sun, Y. Initial Report on the Development of a Monte CarloMarkov Chain Joint Inversion Approach for Geothermal Exploration. United States. doi:10.2172/920472.
Foxall, W, Ramirez, A, Carlson, S, Dyer, K, and Sun, Y. Wed .
"Initial Report on the Development of a Monte CarloMarkov Chain Joint Inversion Approach for Geothermal Exploration". United States.
doi:10.2172/920472. https://www.osti.gov/servlets/purl/920472.
@article{osti_920472,
title = {Initial Report on the Development of a Monte CarloMarkov Chain Joint Inversion Approach for Geothermal Exploration},
author = {Foxall, W and Ramirez, A and Carlson, S and Dyer, K and Sun, Y},
abstractNote = {Geothermal exploration and subsequent characterization of potential resources typically employ a variety of geophysical, geologic and geochemical techniques. However, since the data collected by each technique provide information directly on only one or a very limited set of the many physical parameters that characterize a geothermal system, no single method can be used to describe the system in its entirety. Presently, the usual approach to analyzing disparate data streams for geothermal applications is to invert (or forward model) each data set separately and then combine or compare the resulting models, for the most part in a more or less ad hoc manner. However, while each inversion may yield a model that fits the individual data set, the models are usually inconsistent with each other to some degree. This reflects uncertainties arising from the inevitable fact that geophysical and other exploration data in general are to some extent noisy, incomplete, and of limited sensitivity and resolution, and so yield nonunique results. The purpose of the project described here is to integrate the different model constraints provided by disparate geophysical, geological and geochemical data in a rigorous and consistent manner by formal joint inversion. The objective is to improve the fidelity of exploration results and reservoir characterization, thus addressing the goal of the DOE Geothermal Program to improve success in exploration for economically viable resources by better defining drilling targets, reducing risk, and improving exploration/drilling success rates.},
doi = {10.2172/920472},
journal = {},
number = ,
volume = ,
place = {United States},
year = {Wed Apr 25 00:00:00 EDT 2007},
month = {Wed Apr 25 00:00:00 EDT 2007}
}

We explore the full parameter space of Minimal Supergravity (mSUGRA), allowing all four continuous parameters (the scalar mass m{sub 0}, the gaugino mass m{sub 1/2}, the trilinear coupling A{sub 0}, and the ratio of Higgs vacuum expectation values tan {beta}) to vary freely. We apply current accelerator constraints on sparticle and Higgs masses, and on the b {yields} s{gamma} branching ratio, and discuss the impact of the constraints on g{sub {mu}}2. To study dark matter, we apply the WMAP constraint on the cold dark matter density. We develop Markov Chain Monte Carlo (MCMC) techniques to explore the parameter regions consistentmore »

MultiResolution MarkovChainMonteCarlo Approach for System Identification with an Application to FiniteElement Models
Estimating unknown system configurations/parameters by combining system knowledge gained from a computer simulation model on one hand and from observed data on the other hand is challenging. An example of such inverse problem is detecting and localizing potential flaws or changes in a structure by using a finiteelement model and measured vibration/displacement data. We propose a probabilistic approach based on Bayesian methodology. This approach does not only yield a single bestguess solution, but a posterior probability distribution over the parameter space. In addition, the Bayesian approach provides a natural framework to accommodate prior knowledge. A Markov chain Monte Carlo (MCMC)more » 
Stochastic Inversion of Electrical Resistivity Changes Using a Markov Chain, Monte Carlo Approach
We describe a stochastic inversion method for mapping subsurface regions where the electrical resistivity is changing. The technique combines prior information, electrical resistance data and forward models to produce subsurface resistivity models that are most consistent with all available data. Bayesian inference and a Metropolis simulation algorithm form the basis for this approach. Attractive features include its ability to: (1) provide quantitative measures of the uncertainty of a generated estimate and, (2) allow alternative model estimates to be identified, compared and ranked. Methods that monitor convergence and summarize important trends of the posterior distribution are introduced. Results from a physicalmore » 
Joint reconstructions of CO2 plumes using a Markov Chain Monte Carlo approach
We describe a stochastic inversion method for mapping subsurface regions where CO{sub 2} saturation is changing. The technique combines prior information with measurements of injected CO{sub 2} volume, reservoir deformation and electrical resistivity. Bayesian inference and a Metropolis simulation algorithm form the basis for this approach. The method can (a) jointly reconstruct disparate data types such as surface or subsurface tilt, electrical resistivity, and injected CO{sub 2} volume measurements, (b) provide quantitative measures of the result uncertainty, (c) identify competing models when the available data are insufficient to definitively identify a single optimal model and (d) rank the alternative modelsmore » 
Stochastic Engine Final Report: Applying Markov Chain Monte Carlo Methods with Importance Sampling to LargeScale DataDriven Simulation
Accurate prediction of complex phenomena can be greatly enhanced through the use of data and observations to update simulations. The ability to create these datadriven simulations is limited by error and uncertainty in both the data and the simulation. The stochastic engine project addressed this problem through the development and application of a family of Markov Chain Monte Carlo methods utilizing importance sampling driven by forward simulators to minimize time spent search very large state spaces. The stochastic engine rapidly chooses among a very large number of hypothesized states and selects those that are consistent (within error) with all themore »