Monte Carlo full waveform inversion of tomographic crosshole data using complex geostatistical a priori information

2010 ◽  
Author(s):  
Knud S. Cordua ◽  
Thomas M. Hansen ◽  
Klaus Mosegaard
Keyword(s):  
Monte Carlo ◽  
A Priori ◽  
Waveform Inversion ◽  
Full Waveform Inversion ◽  
A Priori Information ◽  
Full Waveform ◽  
Priori Information

scholarly journals Monte Carlo full-waveform inversion of crosshole GPR data using multiple-point geostatistical a priori information

Geophysics ◽  
2012 ◽  
Vol 77 (2) ◽  
pp. H19-H31 ◽  
Author(s):  
Knud Skou Cordua ◽  
Thomas Mejer Hansen ◽  
Klaus Mosegaard
Keyword(s):  
Inverse Problem ◽  
Monte Carlo ◽  
A Priori ◽  
Waveform Inversion ◽  
Multiple Point ◽  
Full Waveform Inversion ◽  
A Posteriori ◽  
A Priori Information ◽  
Full Waveform ◽  
Priori Information

We present a general Monte Carlo full-waveform inversion strategy that integrates a priori information described by geostatistical algorithms with Bayesian inverse problem theory. The extended Metropolis algorithm can be used to sample the a posteriori probability density of highly nonlinear inverse problems, such as full-waveform inversion. Sequential Gibbs sampling is a method that allows efficient sampling of a priori probability densities described by geostatistical algorithms based on either two-point (e.g., Gaussian) or multiple-point statistics. We outline the theoretical framework for a full-waveform inversion strategy that integrates the extended Metropolis algorithm with sequential Gibbs sampling such that arbitrary complex geostatistically defined a priori information can be included. At the same time we show how temporally and/or spatiallycorrelated data uncertainties can be taken into account during the inversion. The suggested inversion strategy is tested on synthetic tomographic crosshole ground-penetrating radar full-waveform data using multiple-point-based a priori information. This is, to our knowledge, the first example of obtaining a posteriori realizations of a full-waveform inverse problem. Benefits of the proposed methodology compared with deterministic inversion approaches include: (1) The a posteriori model variability reflects the states of information provided by the data uncertainties and a priori information, which provides a means of obtaining resolution analysis. (2) Based on a posteriori realizations, complicated statistical questions can be answered, such as the probability of connectivity across a layer. (3) Complex a priori information can be included through geostatistical algorithms. These benefits, however, require more computing resources than traditional methods do. Moreover, an adequate knowledge of data uncertainties and a priori information is required to obtain meaningful uncertainty estimates. The latter may be a key challenge when considering field experiments, which will not be addressed here.

Improvements to elastic full-waveform inversion using cross-gradient constraints

Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. R105-R115 ◽  
Author(s):  
Edgar Manukyan ◽  
Hansruedi Maurer ◽  
André Nuber
Keyword(s):  
Field Data ◽  
A Priori ◽  
Waveform Inversion ◽  
Structural Similarity ◽  
Model Parameters ◽  
Full Waveform Inversion ◽  
A Priori Information ◽  
Individual Parameter ◽  
Full Waveform ◽  
Priori Information

Seismic full-waveform inversion (FWI) is potentially a powerful method for obtaining high-resolution subsurface images, but the results are often distorted by nonlinear effects and parameter trade-offs. Such distortions can be particularly severe in the case of multiparameter FWI, such as elastic FWI, in which inversion is performed simultaneously for P- and S-wave velocities and density. The problem can be alleviated by adding constraints in the form of plausible a priori information. A usually well-justified constraint includes the structural similarity of different model parameters; i.e., an anomalous body likely exhibits variations in all elastic properties, although their magnitudes may be different. To consider such types of a priori information, we have developed a structurally constrained elastic FWI, which is based on minimization of the cross products of gradients of different model parameters. Our synthetic 2D experiments show that structurally constrained FWI can significantly improve model reconstruction. It is also demonstrated that our approach still leads to improved results, even when the structural similarity between the individual parameter types is not exactly met. Inversions of field data show that in comparison to conventional FWI, structurally constrained FWI is able to match the field data equally well while requiring less structural complexity of the subsurface.

scholarly journals A synthetic study to assess the applicability of full-waveform inversion to infer snow stratigraphy from upward-looking ground-penetrating radar data

Geophysics ◽  
2016 ◽  
Vol 81 (1) ◽  
pp. WA213-WA223 ◽  
Author(s):  
Lino Schmid ◽  
Jürg Schweizer ◽  
John Bradford ◽  
Hansruedi Maurer
Keyword(s):  
Water Content ◽  
Liquid Water ◽  
A Priori ◽  
Waveform Inversion ◽  
Liquid Water Content ◽  
A Priori Information ◽  
Full Waveform ◽  
Snow Stratigraphy ◽  
Priori Information ◽  
Ground Penetrating

Snow stratigraphy and liquid water content are key contributing factors to avalanche formation. Upward-looking ground-penetrating radar (upGPR) systems allow nondestructive monitoring of the snowpack, but deriving density and liquid water content profiles is not yet possible based on the direct analysis of the reflection response. We have investigated the feasibility of deducing these quantities using full-waveform inversion (FWI) techniques applied to upGPR data. For that purpose, we have developed a frequency-domain FWI algorithm in which we additionally took advantage of time-domain features such as the arrival times of reflected waves. Our results indicated that FWI applied to upGPR data is generally feasible. More specifically, we could show that in the case of a dry snowpack, it is possible to derive snow densities and layer thicknesses if sufficient a priori information is available. In case of a wet snowpack, in which it also needs to be inverted for the liquid water content, the algorithm might fail, even if sufficient a priori information is available, particularly in the presence of realistic noise. Finally, we have investigated the capability of FWI to resolve thin layers that play a key role in snow stability evaluation. Our simulations indicate that layers with thicknesses well below the GPR wavelengths can be identified, but in the presence of significant liquid water, the thin-layer properties may be prone to inaccuracies. These results are encouraging and motivate applications to field data, but significant issues remain to be resolved, such as the determination of the generally unknown upGPR source function and identifying the optimal number of layers in the inversion models. Furthermore, a relatively high level of prior knowledge is required to let the algorithm converge. However, we feel these are not insurmountable and the new technology has significant potential to improve field data analysis.

A Backward Monte Carlo method for fast-ion-loss simulations

2021 ◽  
Author(s):  
Filippo Zonta ◽  
Lucia Sanchis ◽  
Eero Hirvijoki
Keyword(s):  
Monte Carlo ◽  
Monte Carlo Method ◽  
A Priori ◽  
Target Surface ◽  
A Priori Information ◽  
Ion Distributions ◽  
Plasma Facing Components ◽  
Priori Information ◽  
Fast Ion

Abstract This paper presents a novel scheme to improve the statistics of simulated fast-ion loss signals and power loads to plasma-facing components in fusion devices. With the so-called Backward Monte Carlo method, the probabilities of marker particles reaching a chosen target surface can be approximately traced from the target back into the plasma. Utilizing the probabilities as {\it a priori} information for the well-established Forward Monte Carlo method, statistics in fast-ion simulations are significantly improved. For testing purposes, the scheme has been implemented to the ASCOT suite of codes and applied to a realistic ASDEX Upgrade configuration of beam-ion distributions.

scholarly journals Two-grid full-waveform Rayleigh-wave inversion via a genetic algorithm — Part 1: Method and synthetic examples

Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. R805-R814 ◽  
Author(s):  
Zhen Xing ◽  
Alfredo Mazzotti
Keyword(s):  
Genetic Algorithm ◽  
Rayleigh Wave ◽  
A Priori ◽  
Computing Time ◽  
A Priori Information ◽  
Fine Grid ◽  
Near Surface ◽  
Reference Models ◽  
Full Waveform ◽  
Priori Information

When reliable a priori information is not available, it is difficult to correctly predict near-surface S-wave velocity models from Rayleigh waves through existing techniques, especially in the case of complex geology. To tackle this issue, we have developed a new method: two-grid genetic-algorithm Rayleigh-wave full-waveform inversion (FWI). Adopting a two-grid parameterization of the model, the genetic algorithm inverts for unknown velocities and densities at the nodes of a coarse grid, whereas the forward modeling is performed on a fine grid to avoid numerical dispersion. A bilinear interpolation brings the coarse-grid results into the fine-grid models. The coarse inversion grid allows for a significant reduction in the computing time required by the genetic algorithm to converge. With a coarser grid, there are fewer unknowns and less required computing time, at the expense of the model resolution. To further increase efficiency, our inversion code can perform the optimization using an offset-marching strategy and/or a frequency-marching strategy that can make use of different kinds of objective functions and allows for parallel computing. We illustrate the effect of our inversion method using three synthetic examples with rather complex near-surface models. Although no a priori information was used in all three tests, the long-wavelength structures of the reference models were fairly predicted, and satisfactory matches between “observed” and predicted data were achieved. The fair predictions of the reference models suggest that the final models estimated by our genetic-algorithm FWI, which we call macromodels, would be suitable inputs to gradient-based Rayleigh-wave FWI for further refinement. We also explored other issues related to the practical use of the method in different work and explored applications of the method to field data.

Anisotropic full-waveform inversion with tilt-angle recovery

Geophysics ◽  
2017 ◽  
Vol 82 (3) ◽  
pp. R135-R151 ◽  
Author(s):  
Herurisa Rusmanugroho ◽  
Ryan Modrak ◽  
Jeroen Tromp
Keyword(s):  
Tilt Angle ◽  
Transverse Isotropy ◽  
A Priori ◽  
Waveform Inversion ◽  
Transversely Isotropic ◽  
Seismic Wave Propagation ◽  
Full Waveform Inversion ◽  
Full Waveform ◽  
Anisotropic Strength ◽  
Global Seismology

By allowing spatial variations in the direction of the anisotropic fast axis, tilted transverse isotropy (TTI) helps to image complex or steeply dipping structures. Without a priori geologic constraints, however, recovery of all the anisotropic parameters can be nontrivial and nonunique. We adopt two methods for TTI inversion with tilt-angle recovery: one based on the familiar Voigt parameters, and another based on the so-called Chen and Tromp parameters known from regional and global seismology. These parameterizations arise naturally in seismic wave propagation and facilitate straightforward recovery of the tilt angle and anisotropic strength. In numerical experiments with vertical transversely isotropic starting models and TTI target models, we find that the Voigt as well as the Chen and Tromp parameters allow quick and robust recovery of steeply dipping anticlinal structures.

Sign in / Sign up

Export Citation Format

Share Document