Illumination compensation of shadow zones in extended least squares migrated images by solving the linear inverse problem in tomographic full waveform inversion

2019 ◽  
Author(s):  
Rahul Sarkar ◽  
Biondo Biondi
Keyword(s):  
Inverse Problem ◽  
Least Squares ◽  
Waveform Inversion ◽  
Full Waveform Inversion ◽  
Illumination Compensation ◽  
Linear Inverse Problem ◽  
Full Waveform

scholarly journals Adding Prior Information in FWI through Relative Entropy

Entropy ◽  
2021 ◽  
Vol 23 (5) ◽  
pp. 599
Author(s):  
Danilo Cruz ◽  
João de Araújo ◽  
Carlos da Costa ◽  
Carlos da Silva
Keyword(s):  
Inverse Problem ◽  
Relative Entropy ◽  
Global Minimum ◽  
Prior Information ◽  
Waveform Inversion ◽  
Petroleum Industry ◽  
Synthetic Data ◽  
Full Waveform Inversion ◽  
Full Waveform

Full waveform inversion is an advantageous technique for obtaining high-resolution subsurface information. In the petroleum industry, mainly in reservoir characterisation, it is common to use information from wells as previous information to decrease the ambiguity of the obtained results. For this, we propose adding a relative entropy term to the formalism of the full waveform inversion. In this context, entropy will be just a nomenclature for regularisation and will have the role of helping the converge to the global minimum. The application of entropy in inverse problems usually involves formulating the problem, so that it is possible to use statistical concepts. To avoid this step, we propose a deterministic application to the full waveform inversion. We will discuss some aspects of relative entropy and show three different ways of using them to add prior information through entropy in the inverse problem. We use a dynamic weighting scheme to add prior information through entropy. The idea is that the prior information can help to find the path of the global minimum at the beginning of the inversion process. In all cases, the prior information can be incorporated very quickly into the full waveform inversion and lead the inversion to the desired solution. When we include the logarithmic weighting that constitutes entropy to the inverse problem, we will suppress the low-intensity ripples and sharpen the point events. Thus, the addition of entropy relative to full waveform inversion can provide a result with better resolution. In regions where salt is present in the BP 2004 model, we obtained a significant improvement by adding prior information through the relative entropy for synthetic data. We will show that the prior information added through entropy in full-waveform inversion formalism will prove to be a way to avoid local minimums.

Full-waveform inversion of salt models using shape optimization and simulated annealing

Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. R793-R804 ◽  
Author(s):  
Debanjan Datta ◽  
Mrinal K. Sen ◽  
Faqi Liu ◽  
Scott Morton
Keyword(s):  
Simulated Annealing ◽  
Least Squares ◽  
Level Set ◽  
Waveform Inversion ◽  
Optimization Techniques ◽  
Model Parameters ◽  
Full Waveform Inversion ◽  
Velocity Models ◽  
Full Waveform ◽  
Starting Model

A good starting model is imperative in full-waveform inversion (FWI) because it solves a least-squares inversion problem using a local gradient-based optimization method. A suboptimal starting model can result in cycle skipping leading to poor convergence and incorrect estimation of subsurface properties. This problem is especially crucial for salt models because the strong velocity contrasts create substantial time shifts in the modeled seismogram. Incorrect estimation of salt bodies leads to velocity inaccuracies in the sediments because the least-squares gradient aims to reduce traveltime differences without considering the sharp velocity jump between sediments and salt. We have developed a technique to estimate velocity models containing salt bodies using a combination of global and local optimization techniques. To stabilize the global optimization algorithm and keep it computationally tractable, we reduce the number of model parameters by using sparse parameterization formulations. The sparse formulation represents sediments using a set of interfaces and velocities across them, whereas a set of ellipses represents the salt body. We use very fast simulated annealing (VFSA) to minimize the misfit between the observed and synthetic data and estimate an optimal model in the sparsely parameterized space. The VFSA inverted model is then used as a starting model in FWI in which the sediments and salt body are updated in the least-squares sense. We partition model updates into sediment and salt updates in which the sediments are updated like conventional FWI, whereas the shape of the salt is updated by taking the zero crossing of an evolving level set surface. Our algorithm is tested on two 2D synthetic salt models, namely, the Sigsbee 2A model and a modified SEG Advanced Modeling Program (SEAM) Phase I model while fixing the top of the salt. We determine the efficiency of the VFSA inversion and imaging improvements from the level set FWI approach and evaluate a few sources of uncertainty in the estimation of salt shapes.

scholarly journals ML-descent: An optimization algorithm for full-waveform inversion using machine learning

Geophysics ◽  
2020 ◽  
Vol 85 (6) ◽  
pp. R477-R492 ◽  
Author(s):  
Bingbing Sun ◽  
Tariq Alkhalifah
Keyword(s):  
Neural Network ◽  
Machine Learning ◽  
Inverse Problem ◽  
Optimization Algorithm ◽  
Waveform Inversion ◽  
Slight Modification ◽  
Descent Method ◽  
Quadratic Functions ◽  
Full Waveform Inversion ◽  
Full Waveform

Full-waveform inversion (FWI) is a nonlinear optimization problem, and a typical optimization algorithm such as the nonlinear conjugate gradient or limited-memory Broyden-Fletcher-Goldfarb-Shanno (LBFGS) would iteratively update the model mainly along the gradient-descent direction of the misfit function or a slight modification of it. Based on the concept of meta-learning, rather than using a hand-designed optimization algorithm, we have trained the machine (represented by a neural network) to learn an optimization algorithm, entitled the “ML-descent,” and apply it in FWI. Using a recurrent neural network (RNN), we use the gradient of the misfit function as the input, and the hidden states in the RNN incorporate the history information of the gradient similar to an LBFGS algorithm. However, unlike the fixed form of the LBFGS algorithm, the machine-learning (ML) version evolves in response to the gradient. The loss function for training is formulated as a weighted summation of the L2 norm of the data residuals in the original inverse problem. As with any well-defined nonlinear inverse problem, the optimization can be locally approximated by a linear convex problem; thus, to accelerate the training, we train the neural network by minimizing randomly generated quadratic functions instead of performing time-consuming FWIs. To further improve the accuracy and robustness, we use a variational autoencoder that projects and represents the model in latent space. We use the Marmousi and the overthrust examples to demonstrate that the ML-descent method shows faster convergence and outperforms conventional optimization algorithms. The energy in the deeper part of the models can be recovered by the ML-descent even when the pseudoinverse of the Hessian is not incorporated in the FWI update.

scholarly journals Constrained Full Waveform Inversion for Borehole Multicomponent Seismic Data

Geosciences ◽  
2019 ◽  
Vol 9 (1) ◽  
pp. 45
Author(s):  
Marwan Charara ◽  
Christophe Barnes
Keyword(s):  
Inverse Problem ◽  
Spatial Correlation ◽  
Seismic Data ◽  
Waveform Inversion ◽  
Synthetic Data ◽  
Model Space ◽  
Elastic Model ◽  
Full Waveform Inversion ◽  
Spatial Correlations ◽  
Full Waveform

Full-waveform inversion for borehole seismic data is an ill-posed problem and constraining the problem is crucial. Constraints can be imposed on the data and model space through covariance matrices. Usually, they are set to a diagonal matrix. For the data space, signal polarization information can be used to evaluate the data uncertainties. The inversion forces the synthetic data to fit the polarization of observed data. A synthetic inversion for a 2D-2C data estimating a 1D elastic model shows a clear improvement, especially at the level of the receivers. For the model space, horizontal and vertical spatial correlations using a Laplace distribution can be used to fill the model space covariance matrix. This approach reduces the degree of freedom of the inverse problem, which can be quantitatively evaluated. Strong horizontal spatial correlation distances favor a tabular geological model whenever it does not contradict the data. The relaxation of the spatial correlation distances from large to small during the iterative inversion process allows the recovery of geological objects of the same size, which regularizes the inverse problem. Synthetic constrained and unconstrained inversions for 2D-2C crosswell data show the clear improvement of the inversion results when constraints are used.

scholarly journals Land seismic multiparameter full waveform inversion in elastic VTI media by simultaneously interpreting body waves and surface waves with an optimal transport based objective function

2019 ◽  
Vol 219 (3) ◽  
pp. 1970-1988 ◽  
Author(s):  
Weiguang He ◽  
Romain Brossier ◽  
Ludovic Métivier ◽  
René-Édouard Plessix
Keyword(s):  
Objective Function ◽  
Surface Waves ◽  
Least Squares ◽  
Optimal Transport ◽  
Waveform Inversion ◽  
Body Waves ◽  
Model Parameters ◽  
Full Waveform Inversion ◽  
Initial Model ◽  
Full Waveform

SUMMARY Land seismic multiparameter full waveform inversion in anisotropic media is challenging because of high medium contrasts and surface waves. With a data-residual least-squares objective function, the surface wave energy usually masks the body waves and the gradient of the objective function exhibits high values in the very shallow depths preventing from recovering the deeper part of the earth model parameters. The optimal transport objective function, coupled with a Gaussian time-windowing strategy, allows to overcome this issue by more focusing on phase shifts and by balancing the contributions of the different events in the adjoint-source and the gradients. We first illustrate the advantages of the optimal transport function with respect to the least-squares one, with two realistic examples. We then discuss a vertical transverse isotropic (VTI) example starting from a quasi 1-D isotropic initial model. Despite some cycle-skipping issues in the initial model, the inversion based on the windowed optimal transport approach converges. Both the near-surface complexities and the variations at depth are recovered.

Sign in / Sign up

Export Citation Format

Share Document