Stochastic full waveform inversion with deep generative prior for uncertainty quantification

Read original: arXiv:2406.04859 - Published 6/10/2024 by Yuke Xie, Herv'e Chauris, Nicolas Desassis
Total Score

0

Stochastic full waveform inversion with deep generative prior for uncertainty quantification

Sign in to get full access

or

If you already have an account, we'll log you in

Overview

  • This paper proposes a stochastic full waveform inversion (FWI) method that uses a deep generative prior to quantify uncertainty in seismic inversion problems.
  • The method combines a stochastic optimization approach with a pre-trained deep generative model to efficiently explore the space of possible subsurface models.
  • By leveraging the powerful representation capabilities of deep generative models, the approach can capture complex geological structures and provide meaningful uncertainty estimates.

Plain English Explanation

The paper presents a new way to analyze seismic data and create models of the Earth's subsurface. Seismic inversion is the process of using seismic waves to infer the properties of the Earth's interior, such as the density and stiffness of the rocks. This is similar to how doctors use medical imaging to get information about the human body.

However, seismic inversion is a complex problem with a lot of uncertainty. The paper's new method tries to address this by combining two key ideas:

  1. Stochastic optimization: Instead of trying to find a single "best" model, the method explores many possible models in a random (stochastic) way. This allows it to better capture the uncertainty in the solution.

  2. Deep generative models: The method uses a pre-trained neural network that can generate realistic-looking models of the subsurface. This is similar to how AI can generate images of imaginary scenes. By leveraging the power of these deep models, the method can efficiently explore a wide range of possible subsurface structures.

Together, these two ideas allow the method to provide not just a single model, but a whole range of possible models that capture the uncertainty in the seismic inversion problem. [This can be very useful for applications like oil and gas exploration or carbon storage, where understanding the uncertainty is crucial.](https://aimodels.fyi/papers/arxiv/principled-probabilistic-imaging-using-diffusion-models-as)

Technical Explanation

The core of the proposed method is a stochastic optimization approach that uses a deep generative model as a prior distribution for the subsurface properties. The deep generative model is pre-trained on a dataset of realistic geological structures, allowing it to capture the complex spatial relationships and heterogeneities that characterize the Earth's subsurface.

During the inversion process, the method iteratively updates a set of latent variables that parameterize the deep generative model. By sampling from the posterior distribution of these latent variables, the method can generate an ensemble of plausible subsurface models that are consistent with the observed seismic data. This allows for a principled quantification of the uncertainty in the inversion results.

The authors demonstrate the effectiveness of their approach on synthetic and field data examples, showing that it can outperform traditional deterministic FWI methods in terms of both accuracy and uncertainty quantification.

Critical Analysis

The proposed method represents an interesting and promising approach to addressing the challenge of uncertainty quantification in seismic inversion problems. By leveraging the power of deep generative models, the method can capture complex geological structures more effectively than traditional parametric approaches.

However, the method also has some potential limitations. The performance of the deep generative model is highly dependent on the quality and representativeness of the training data, which may be difficult to obtain in practice. Additionally, the computational cost of the stochastic optimization process may be prohibitive for large-scale problems, particularly when dealing with high-resolution seismic data.

Further research is needed to investigate the scalability of the method and to explore ways to improve its computational efficiency. It would also be valuable to see the method applied to a wider range of real-world seismic inversion problems to assess its practical utility and robustness.

Conclusion

The paper presents a novel stochastic FWI method that uses a deep generative prior to quantify uncertainty in seismic inversion problems. By combining stochastic optimization with the powerful representation capabilities of deep generative models, the method can efficiently explore the space of possible subsurface models and provide meaningful uncertainty estimates.

This approach has the potential to significantly improve the reliability and interpretability of seismic inversion results, which could have important implications for applications such as oil and gas exploration, carbon storage, and geohazard assessment. While the method has some limitations, the paper represents an important step forward in the field of uncertainty quantification for seismic imaging.



This summary was produced with help from an AI and may contain inaccuracies - check out the links to read the original source documents!

Follow @aimodelsfyi on 𝕏 →

Related Papers

Stochastic full waveform inversion with deep generative prior for uncertainty quantification
Total Score

0

Stochastic full waveform inversion with deep generative prior for uncertainty quantification

Yuke Xie, Herv'e Chauris, Nicolas Desassis

To obtain high-resolution images of subsurface structures from seismic data, seismic imaging techniques such as Full Waveform Inversion (FWI) serve as crucial tools. However, FWI involves solving a nonlinear and often non-unique inverse problem, presenting challenges such as local minima trapping and inadequate handling of inherent uncertainties. In addressing these challenges, we propose leveraging deep generative models as the prior distribution of geophysical parameters for stochastic Bayesian inversion. This approach integrates the adjoint state gradient for efficient back-propagation from the numerical solution of partial differential equations. Additionally, we introduce explicit and implicit variational Bayesian inference methods. The explicit method computes variational distribution density using a normalizing flow-based neural network, enabling computation of the Bayesian posterior of parameters. Conversely, the implicit method employs an inference network attached to a pretrained generative model to estimate density, incorporating an entropy estimator. Furthermore, we also experimented with the Stein Variational Gradient Descent (SVGD) method as another variational inference technique, using particles. We compare these variational Bayesian inference methods with conventional Markov chain Monte Carlo (McMC) sampling. Each method is able to quantify uncertainties and to generate seismic data-conditioned realizations of subsurface geophysical parameters. This framework provides insights into subsurface structures while accounting for inherent uncertainties.

Read more

6/10/2024

Elastic Full-Waveform Inversion : How the physics of problem improves data-driven techniques?
Total Score

0

Elastic Full-Waveform Inversion : How the physics of problem improves data-driven techniques?

Vahid Negahdari, Seyed Reza Moghadasi, Mohammad Reza Razvan

Full-Waveform Inversion (FWI) is a nonlinear iterative seismic imaging technique that, by reducing the misfit between recorded and predicted seismic waveforms, can produce detailed estimates of subsurface geophysical properties. Nevertheless, the strong nonlinearity of FWI can trap the optimization in local minima. This issue arises due to factors such as improper initial values, the absence of low frequencies in the measurements, noise, and other related considerations. To address this challenge and with the advent of advanced machine-learning techniques, data-driven methods, such as deep learning, have attracted significantly increasing attention in the geophysical community. Furthermore, the elastic wave equation should be included in FWI to represent elastic effects accurately. The intersection of data-driven techniques and elastic scattering theories presents opportunities and challenges. In this paper, by using the knowledge of elastic scattering (Physics of problem) and integrating it with deep learning techniques, we propose methods for the solution of time-harmonic FWI to enhance accuracy compared to pure data-driven approaches. Moreover, by modifying the structure of the Variational Autoencoder, we introduce a probabilistic deep learning method based on the physics of the problem that enables us to explore the uncertainties of the solution. According to the limited availability of datasets in this field and to assess the performance and accuracy of the proposed methods, we create a comprehensive dataset close to reality and conduct a comparative analysis of the presented approaches to it.

Read more

6/11/2024

Physics-guided Full Waveform Inversion using Encoder-Solver Convolutional Neural Networks
Total Score

0

Physics-guided Full Waveform Inversion using Encoder-Solver Convolutional Neural Networks

Matan Goren, Eran Treister

Full Waveform Inversion (FWI) is an inverse problem for estimating the wave velocity distribution in a given domain, based on observed data on the boundaries. The inversion is computationally demanding because we are required to solve multiple forward problems, either in time or frequency domains, to simulate data that are then iteratively fitted to the observed data. We consider FWI in the frequency domain, where the Helmholtz equation is used as a forward model, and its repeated solution is the main computational bottleneck of the inversion process. To ease this cost, we integrate a learning process of an encoder-solver preconditioner that is based on convolutional neural networks (CNNs). The encoder-solver is trained to effectively precondition the discretized Helmholtz operator given velocity medium parameters. Then, by re-training the CNN between the iterations of the optimization process, the encoder-solver is adapted to the iteratively evolving velocity medium as part of the inversion. Without retraining, the performance of the solver deteriorates as the medium changes. Using our light retraining procedures, we obtain the forward simulations effectively throughout the process. We demonstrate our approach to solving FWI problems using 2D geophysical models with high-frequency data.

Read more

5/29/2024

Inversion-DeepONet: A Novel DeepONet-Based Network with Encoder-Decoder for Full Waveform Inversion
Total Score

0

Inversion-DeepONet: A Novel DeepONet-Based Network with Encoder-Decoder for Full Waveform Inversion

Zekai Guo, Lihui Chai, Shengjun Huang, Ye Li

Full waveform inversion (FWI) plays a crucial role in the field of geophysics. There has been lots of research about applying deep learning (DL) methods to FWI. The success of DL-FWI relies significantly on the quantity and diversity of the datasets. Nevertheless, existing FWI datasets, like OpenFWI, where sources have fixed locations or identical frequencies, provide limited information and do not represent the complex real-world scene. For instance, low frequencies help in resolving larger-scale structures. High frequencies allow for a more detailed subsurface features. %A single source frequency is insufficient to describe subsurface structural properties. We consider that simultaneously using sources with different frequencies, instead of performing inversion using low frequencies data and then gradually introducing higher frequencies data, has rationale and potential advantages. Hence, we develop three enhanced datasets based on OpenFWI where each source have varying locations, frequencies or both. Moreover, we propose a novel deep operator network (DeepONet) architecture Inversion-DeepONet for FWI. We utilize convolutional neural network (CNN) to extract the features from seismic data in branch net. Source parameters, such as locations and frequencies, are fed to trunk net. Then another CNN is employed as the decoder of DeepONet to reconstruct the velocity models more effectively. Through experiments, we confirm the superior performance on accuracy and generalization ability of our network, compared with existing data-driven FWI methods.

Read more

8/16/2024