Optimized Model Parameterization using Compact Full Waveform Inversion

2020 ◽  
Author(s):  
Linan Xu ◽  
Edgar Manukyan ◽  
Hansruedi Maurer
Keyword(s):  
Large Scale ◽  
Waveform Inversion ◽  
Gradient Methods ◽  
Significant Loss ◽  
Inverse Model ◽  
Eigenvalue Decomposition ◽  
Model Parameters ◽  
Full Waveform Inversion ◽  
Full Waveform ◽  
Model Parameterization

Summary Seismic Full Waveform Inversion (FWI) has the potential to produce high-resolution subsurface images, but the computational resources required for realistically sized problems can be prohibitively large. In terms of computational costs, Gauss-Newton algorithms are more attractive than the commonly employed conjugate gradient methods, because the former have favorable convergence properties. However, efficient implementations of Gauss-Newton algorithms require an excessive amount of computer memory for larger problems. To address this issue, we introduce Compact Full Waveform Inversion (CFWI). Here, a suitable inverse model parameterization is sought that allows representing all subsurface features, potentially resolvable by a particular source-receiver deployment, but using only a minimum number of model parameters. In principle, an inverse model parameterization, based on the Eigenvalue decomposition, would be optimal, but this is computationally not feasible for realistic problems. Instead, we present two alternative parameter transformations, namely the Haar and the Hartley transformations, with which similarly good results can be obtained. By means of a suite of numerical experiments, we demonstrate that these transformations allow the number of model parameters to be reduced to only a few percent of the original parameterization without any significant loss of spatial resolution. This facilitates efficient solutions of large-scale FWI problems with explicit Gauss-Newton algorithms.

Full waveform inversion with sparse structure constrained regularization

2018 ◽  
Vol 26 (2) ◽  
pp. 243-257 ◽  
Author(s):  
Zichao Yan ◽  
Yanfei Wang
Keyword(s):  
Large Scale ◽  
Regularization Parameter ◽  
Waveform Inversion ◽  
Difference Operators ◽  
Model Parameters ◽  
Full Waveform Inversion ◽  
Structure Information ◽  
Full Waveform ◽  
Ill Posed ◽  
Scale Optimization

AbstractFull waveform inversion is a large-scale nonlinear and ill-posed problem. We consider applying the regularization technique for full waveform inversion with structure constraints. The structure information was extracted with difference operators with respect to model parameters. And then we establish an {l_{p}}-{l_{q}}-norm constrained minimization model for different choices of parameters p and q. To solve this large-scale optimization problem, a fast gradient method with projection onto convex set and a multiscale inversion strategy are addressed. The regularization parameter is estimated adaptively with respect to the frequency range of the data. Numerical experiments on a layered model and a benchmark SEG/EAGE overthrust model are performed to testify the validity of this proposed regularization scheme.

A new parameter set for anisotropic multiparameter full-waveform inversion and application to a North Sea data set

Geophysics ◽  
2016 ◽  
Vol 81 (4) ◽  
pp. U25-U38 ◽  
Author(s):  
Nuno V. da Silva ◽  
Andrew Ratcliffe ◽  
Vetle Vinje ◽  
Graham Conroy
Keyword(s):  
North Sea ◽  
Waveform Inversion ◽  
Real Data ◽  
Model Parameters ◽  
Full Waveform Inversion ◽  
Radiation Patterns ◽  
Data Set ◽  
Full Waveform ◽  
Seismic Image ◽  
Central North Sea

Parameterization lies at the center of anisotropic full-waveform inversion (FWI) with multiparameter updates. This is because FWI aims to update the long and short wavelengths of the perturbations. Thus, it is important that the parameterization accommodates this. Recently, there has been an intensive effort to determine the optimal parameterization, centering the fundamental discussion mainly on the analysis of radiation patterns for each one of these parameterizations, and aiming to determine which is best suited for multiparameter inversion. We have developed a new parameterization in the scope of FWI, based on the concept of kinematically equivalent media, as originally proposed in other areas of seismic data analysis. Our analysis is also based on radiation patterns, as well as the relation between the perturbation of this set of parameters and perturbation in traveltime. The radiation pattern reveals that this parameterization combines some of the characteristics of parameterizations with one velocity and two Thomsen’s parameters and parameterizations using two velocities and one Thomsen’s parameter. The study of perturbation of traveltime with perturbation of model parameters shows that the new parameterization is less ambiguous when relating these quantities in comparison with other more commonly used parameterizations. We have concluded that our new parameterization is well-suited for inverting diving waves, which are of paramount importance to carry out practical FWI successfully. We have demonstrated that the new parameterization produces good inversion results with synthetic and real data examples. In the latter case of the real data example from the Central North Sea, the inverted models show good agreement with the geologic structures, leading to an improvement of the seismic image and flatness of the common image gathers.

FULL-WAVEFORM INVERSION WITH SOURCE AND RECEIVER COUPLING EFFECTS CORRECTION

Geophysics ◽  
2021 ◽  
pp. 1-37
Author(s):  
Xinhai Hu ◽  
Wei Guoqi ◽  
Jianyong Song ◽  
Zhifang Yang ◽  
Minghui Lu ◽  
...  
Keyword(s):  
Waveform Inversion ◽  
Nonlinear Least Squares ◽  
Real Data ◽  
Model Parameters ◽  
Full Waveform Inversion ◽  
Linear Inversion ◽  
Near Surface ◽  
Model Parameter ◽  
Full Waveform ◽  
Coupling Factors

Coupling factors of sources and receivers vary dramatically due to the strong heterogeneity of near surface, which are as important as the model parameters for the inversion success. We propose a full waveform inversion (FWI) scheme that corrects for variable coupling factors while updating the model parameter. A linear inversion is embedded into the scheme to estimate the source and receiver factors and compute the amplitude weights according to the acquisition geometry. After the weights are introduced in the objective function, the inversion falls into the category of separable nonlinear least-squares problems. Hence, we could use the variable projection technique widely used in source estimation problem to invert the model parameter without the knowledge of source and receiver factors. The efficacy of the inversion scheme is demonstrated with two synthetic examples and one real data test.

3D elastic full-waveform inversion of surface waves in the presence of irregular topography using an envelope-based misfit function

Geophysics ◽  
2018 ◽  
Vol 83 (1) ◽  
pp. R1-R11 ◽  
Author(s):  
Dmitry Borisov ◽  
Ryan Modrak ◽  
Fuchun Gao ◽  
Jeroen Tromp
Keyword(s):  
Surface Wave ◽  
Objective Function ◽  
Surface Waves ◽  
Large Scale ◽  
Body Wave ◽  
Waveform Inversion ◽  
Full Waveform Inversion ◽  
S Wave ◽  
Near Surface ◽  
Full Waveform

Full-waveform inversion (FWI) is a powerful method for estimating the earth’s material properties. We demonstrate that surface-wave-driven FWI is well-suited to recovering near-surface structures and effective at providing S-wave speed starting models for use in conventional body-wave FWI. Using a synthetic example based on the SEG Advanced Modeling phase II foothills model, we started with an envelope-based objective function to invert for shallow large-scale heterogeneities. Then we used a waveform-difference objective function to obtain a higher-resolution model. To accurately model surface waves in the presence of complex tomography, we used a spectral-element wave-propagation solver. Envelope misfit functions are found to be effective at minimizing cycle-skipping issues in surface-wave inversions, and surface waves themselves are found to be useful for constraining complex near-surface features.

Efficient 3D elastic full-waveform inversion using wavefield reconstruction methods

Geophysics ◽  
2016 ◽  
Vol 81 (2) ◽  
pp. R45-R55 ◽  
Author(s):  
Espen Birger Raknes ◽  
Wiktor Weibull
Keyword(s):  
Particle Velocity ◽  
Large Scale ◽  
Waveform Inversion ◽  
Transversely Isotropic ◽  
Reconstruction Method ◽  
Full Waveform Inversion ◽  
Reverse Time ◽  
Full Waveform ◽  
Time Migration ◽  
Reconstruction Methods

In reverse time migration (RTM) or full-waveform inversion (FWI), forward and reverse time propagating wavefields are crosscorrelated in time to form either the image condition in RTM or the misfit gradient in FWI. The crosscorrelation condition requires both fields to be available at the same time instants. For large-scale 3D problems, it is not possible, in practice, to store snapshots of the wavefields during forward modeling due to extreme storage requirements. We have developed an approximate wavefield reconstruction method that uses particle velocity field recordings on the boundaries to reconstruct the forward wavefields during the computation of the reverse time wavefields. The method is computationally effective and requires less storage than similar methods. We have compared the reconstruction method to a boundary reconstruction method that uses particle velocity and stress fields at the boundaries and to the optimal checkpointing method. We have tested the methods on a 2D vertical transversely isotropic model and a large-scale 3D elastic FWI problem. Our results revealed that there are small differences in the results for the three methods.

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.

Image-guided sparse-model full waveform inversion

Geophysics ◽  
2012 ◽  
Vol 77 (4) ◽  
pp. R189-R198 ◽  
Author(s):  
Yong Ma ◽  
Dave Hale ◽  
Bin Gong ◽  
Zhaobo (Joe) Meng
Keyword(s):  
Waveform Inversion ◽  
Model Space ◽  
Model Parameters ◽  
Full Waveform Inversion ◽  
Widespread Application ◽  
Low Frequencies ◽  
Sparse Model ◽  
Image Guided ◽  
Full Waveform ◽  
Recorded Data

Multiple problems, including high computational cost, spurious local minima, and solutions with no geologic sense, have prevented widespread application of full waveform inversion (FWI), especially FWI of seismic reflections. These problems are fundamentally related to a large number of model parameters and to the absence of low frequencies in recorded seismograms. Instead of inverting for all the parameters in a dense model, image-guided full waveform inversion inverts for a sparse model space that contains far fewer parameters. We represent a model with a sparse set of values, and from these values, we use image-guided interpolation (IGI) and its adjoint operator to compute finely and uniformly sampled models that can fit recorded data in FWI. Because of this sparse representation, image-guided FWI updates more blocky models, and this blockiness in the model space mitigates the absence of low frequencies in recorded data. Moreover, IGI honors imaged structures, so image-guided FWI built in this way yields models that are geologically sensible.

Simultaneous estimation of velocity and density in acoustic multiparameter full-waveform inversion using an improved scattering-integral approach

Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. R399-R415 ◽  
Author(s):  
Jizhong Yang ◽  
Yuzhu Liu ◽  
Liangguo Dong
Keyword(s):  
Newton Method ◽  
Waveform Inversion ◽  
Integral Approach ◽  
Normal Equation ◽  
Simultaneous Estimation ◽  
Model Parameters ◽  
Full Waveform Inversion ◽  
Trade Off ◽  
Full Waveform ◽  
Diffraction Patterns

Density is known to be difficult to reconstruct in multiparameter full-waveform inversion (FWI). This difficulty results from the similarity of the diffraction patterns of velocity and density at small scattering angles. In addition, the sensitivities of seismic data with respect to velocity and density have different orders of magnitude which make the inversion ill-conditioned. The inverse Hessian has been shown to mitigate the coupling effects and rescale the magnitudes of different parameters, such that reliable updates for all parameters are available. We have investigated the possibility of simultaneous estimations of velocity and density in acoustic media using the truncated Gauss-Newton method. The model updates are calculated using a matrix-free conjugate gradient solution of the Gauss-Newton normal equation. The gradients of the misfit function with respect to the model parameters and the Hessian-vector products are computed using an improved scattering-integral approach. To give some insights into the trade-off effects between velocity and density, and the imaging resolution in FWI, the sensitivity kernels of both parameters are numerically calculated in homogeneous background models, and their spatial distributions and characteristics are analyzed. The synthetic experiments on a canonical inclusion model and the 2004 BP model confirm that, in cases in which the Gauss-Newton approximate Hessian, especially its off-diagonal blocks, is accurately taken into account, the truncated Gauss-Newton method can effectively mitigate the trade-off effects between velocity and density and provide accelerated convergence rate. Hence, well-resolved velocity and density models are expected.

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