A homotopy method for well-log constraint waveform inversion

Geophysics ◽  
2007 ◽  
Vol 72 (1) ◽  
pp. R1-R7 ◽  
Author(s):  
Bo Han ◽  
Hongsun Fu ◽  
Hong Liu
Keyword(s):  
Iterative Methods ◽  
Global Minimum ◽  
Waveform Inversion ◽  
Homotopy Method ◽  
Well Log ◽  
Local Minima ◽  
Acoustic Wave Equation ◽  
Low Frequencies ◽  
True Model ◽  
Starting Model

The primary difficulty posed by the method of waveform inversion for the acoustic-wave equation is the presence of local minima of the objective function. Waveform inversion fails to converge to the global minimum unless the wavefield contains very low frequencies, or the starting model is very close to the true model. Constraints can improve the convergence of the method; however, if the starting model is far from the correct model, normal iterative methods will still get trapped in local minima. We designed a homotopy method to improve the robustness of waveform inversion; it makes natural use of constraints, such as well logs. In addition, to further condition the inverse problem, we incorporate standard Tikhonov regularization. We demonstrate through a synthetic example that our method is more likely to find a global minimum than normal iterative methods.

Multiscale seismic waveform inversion

Geophysics ◽  
1995 ◽  
Vol 60 (5) ◽  
pp. 1457-1473 ◽  
Author(s):  
Carey Bunks ◽  
Fatimetou M. Saleck ◽  
S. Zaleski ◽  
G. Chavent
Keyword(s):  
Iterative Methods ◽  
Global Minimum ◽  
Multigrid Method ◽  
Waveform Inversion ◽  
Seismic Inversion ◽  
Local Minima ◽  
Data Set ◽  
Seismic Waveform ◽  
Inversion Methods ◽  
Iterative Inversion

Iterative inversion methods have been unsuccessful at inverting seismic data obtained from complicated earth models (e.g. the Marmousi model), the primary difficulty being the presence of numerous local minima in the objective function. The presence of local minima at all scales in the seismic inversion problem prevent iterative methods of inversion from attaining a reasonable degree of convergence to the neighborhood of the global minimum. The multigrid method is a technique that improves the performance of iterative inversion by decomposing the problem by scale. At long scales there are fewer local minima and those that remain are further apart from each other. Thus, at long scales iterative methods can get closer to the neighborhood of the global minimum. We apply the multigrid method to a subsampled, low‐frequency version of the Marmousi data set. Although issues of source estimation, source bandwidth, and noise are not treated, results show that iterative inversion methods perform much better when employed with a decomposition by scale. Furthermore, the method greatly reduces the computational burden of the inversion that will be of importance for 3-D extensions to the method.

scholarly journals Adaptive waveform inversion: Theory

Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. R429-R445 ◽  
Author(s):  
Michael Warner ◽  
Lluís Guasch
Keyword(s):  
Waveform Inversion ◽  
Computational Cost ◽  
Synthetic Data ◽  
Velocity Model ◽  
Seismic Inversion ◽  
New Method ◽  
Low Frequencies ◽  
Full Waveform ◽  
True Model ◽  
Starting Model

Conventional full-waveform seismic inversion attempts to find a model of the subsurface that is able to predict observed seismic waveforms exactly; it proceeds by minimizing the difference between the observed and predicted data directly, iterating in a series of linearized steps from an assumed starting model. If this starting model is too far removed from the true model, then this approach leads to a spurious model in which the predicted data are cycle skipped with respect to the observed data. Adaptive waveform inversion (AWI) provides a new form of full-waveform inversion (FWI) that appears to be immune to the problems otherwise generated by cycle skipping. In this method, least-squares convolutional filters are designed that transform the predicted data into the observed data. The inversion problem is formulated such that the subsurface model is iteratively updated to force these Wiener filters toward zero-lag delta functions. As that is achieved, the predicted data evolve toward the observed data and the assumed model evolves toward the true model. This new method is able to invert synthetic data successfully, beginning from starting models and under conditions for which conventional FWI fails entirely. AWI has a similar computational cost to conventional FWI per iteration, and it appears to converge at a similar rate. The principal advantages of this new method are that it allows waveform inversion to begin from less-accurate starting models, does not require the presence of low frequencies in the field data, and appears to provide a better balance between the influence of refracted and reflected arrivals upon the final-velocity model. The AWI is also able to invert successfully when the assumed source wavelet is severely in error.

Projection methods and applications for seismic nonlinear inverse problems with multiple constraints

Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. R251-R269 ◽  
Author(s):  
Bas Peters ◽  
Brendan R. Smithyman ◽  
Felix J. Herrmann
Keyword(s):  
Inverse Problems ◽  
Prior Information ◽  
Waveform Inversion ◽  
Projection Methods ◽  
Multiple Constraints ◽  
Local Minima ◽  
Low Frequencies ◽  
Trade Off ◽  
Nonlinear Inverse Problems ◽  
Easy Integration

Nonlinear inverse problems are often hampered by local minima because of missing low frequencies and far offsets in the data, lack of access to good starting models, noise, and modeling errors. A well-known approach to counter these deficiencies is to include prior information on the unknown model, which regularizes the inverse problem. Although conventional regularization methods have resulted in enormous progress in ill-posed (geophysical) inverse problems, challenges remain when the prior information consists of multiple pieces. To handle this situation, we have developed an optimization framework that allows us to add multiple pieces of prior information in the form of constraints. The proposed framework is more suitable for full-waveform inversion (FWI) because it offers assurances that multiple constraints are imposed uniquely at each iteration, irrespective of the order in which they are invoked. To project onto the intersection of multiple sets uniquely, we use Dykstra’s algorithm that does not rely on trade-off parameters. In that sense, our approach differs substantially from approaches, such as Tikhonov/penalty regularization and gradient filtering. None of these offer assurances, which makes them less suitable to FWI, where unrealistic intermediate results effectively derail the inversion. By working with intersections of sets, we avoid trade-off parameters and keep objective calculations separate from projections that are often much faster to compute than objectives/gradients in 3D. These features allow for easy integration into existing code bases. Working with constraints also allows for heuristics, where we built up the complexity of the model by a gradual relaxation of the constraints. This strategy helps to avoid convergence to local minima that represent unrealistic models. Using multiple constraints, we obtain better FWI results compared with a quadratic penalty method, whereas all definitions of the constraints are in terms of physical units and follow from the prior knowledge directly.

scholarly journals Full-waveform inversion with extrapolated low-frequency data

Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. R339-R348 ◽  
Author(s):  
Yunyue Elita Li ◽  
Laurent Demanet
Keyword(s):  
Waveform Inversion ◽  
Low Frequency ◽  
Velocity Model ◽  
Propagation Model ◽  
Full Waveform Inversion ◽  
Frequency Data ◽  
Local Minima ◽  
Bandwidth Extension ◽  
Low Frequencies ◽  
Full Waveform

The availability of low-frequency data is an important factor in the success of full-waveform inversion (FWI) in the acoustic regime. The low frequencies help determine the kinematically relevant, low-wavenumber components of the velocity model, which are in turn needed to avoid convergence of FWI to spurious local minima. However, acquiring data less than 2 or 3 Hz from the field is a challenging and expensive task. We have explored the possibility of synthesizing the low frequencies computationally from high-frequency data and used the resulting prediction of the missing data to seed the frequency sweep of FWI. As a signal-processing problem, bandwidth extension is a very nonlinear and delicate operation. In all but the simplest of scenarios, it can only be expected to lead to plausible recovery of the low frequencies, rather than their accurate reconstruction. Even so, it still requires a high-level interpretation of band-limited seismic records into individual events, each of which can be extrapolated to a lower (or higher) frequency band from the nondispersive nature of the wave-propagation model. We have used the phase-tracking method for the event separation task. The fidelity of the resulting extrapolation method is typically higher in phase than in amplitude. To demonstrate the reliability of bandwidth extension in the context of FWI, we first used the low frequencies in the extrapolated band as data substitute, to create the low-wavenumber background velocity model, and then we switched to recorded data in the available band for the rest of the iterations. The resulting method, extrapolated FWI, demonstrated surprising robustness to the inaccuracies in the extrapolated low-frequency data. With two synthetic examples calibrated so that regular FWI needs to be initialized at 1 Hz to avoid local minima, we have determined that FWI based on an extrapolated [1, 5] Hz band, itself generated from data available in the [5, 15] Hz band, can produce reasonable estimations of the low-wavenumber velocity models.

scholarly journals Full waveform inversion in time and frequency domain of velocity modeling in seismic imaging: FWISIMAT a Matlab code

2018 ◽  
Vol 22 (4) ◽  
pp. 291-300
Author(s):  
Sagar Singh ◽  
Ali Ismet Kanli ◽  
Sagarika Mukhopadhyay
Keyword(s):  
Frequency Domain ◽  
Seismic Imaging ◽  
Waveform Inversion ◽  
Velocity Model ◽  
Full Waveform Inversion ◽  
Initial Model ◽  
Final Model ◽  
Full Waveform ◽  
True Model ◽  
Starting Model

This paper investigates the capability of acoustic Full Waveform Inversion (FWI) in building Marmousi velocity model, in time and frequency domain. FWI is an iterative minimization of misfit between observed and calculated data which is generally solved in three segments: forward modeling, which numerically solves the wave equation with an initial model, gradient computation of the objective function, and updating the model parameters, with a valid optimization method. FWI codes developed in MATLAB herein FWISIMAT (Full Waveform Inversion in Seismic Imaging using MATLAB) are successfully implemented using the Marmousi velocity model as the true model. An initial model is obtained by smoothing the true model to initiate FWI procedure. Smoothing ensures an adequate starting model for FWI, as the FWI procedure is known to be sensitive on the starting model. The final model is compared with the true model to review the number of recovered velocities. FWI codes developed in MATLAB herein FWISIMAT (Full Waveform Inversion in Seismic Imaging using MATLAB) are successfully implemented usingMarmousi velocity model astrue model. An initial model is derived from smoothing the true model to initiate FWI procedure. Smoothing ensures an adequate starting model for FWI, as the FWI procedure is known to be sensitive onstarting model. The final model is compared with the true model to review theamount of recovered velocities. 

scholarly journals Full-waveform inversion using a nonlinearly smoothed wavefield

Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. R117-R127 ◽  
Author(s):  
Yuanyuan Li ◽  
Yunseok Choi ◽  
Tariq Alkhalifah ◽  
Zhenchun Li ◽  
Kai Zhang
Keyword(s):  
Objective Function ◽  
Global Minimum ◽  
Waveform Inversion ◽  
Gradient Methods ◽  
Low Frequency ◽  
Full Waveform Inversion ◽  
Half Cycle ◽  
Low Frequencies ◽  
Full Waveform ◽  
Global Correlation

Conventional full-waveform inversion (FWI) based on the least-squares misfit function faces problems in converging to the global minimum when using gradient methods because of the cycle-skipping phenomena. An initial model producing data that are at most a half-cycle away from the observed data is needed for convergence to the global minimum. Low frequencies are helpful in updating low-wavenumber components of the velocity model to avoid cycle skipping. However, low enough frequencies are usually unavailable in field cases. The multiplication of wavefields of slightly different frequencies adds artificial low-frequency components in the data, which can be used for FWI to generate a convergent result and avoid cycle skipping. We generalize this process by multiplying the wavefield with itself and then applying a smoothing operator to the multiplied wavefield or its square to derive the nonlinearly smoothed wavefield, which is rich in low frequencies. The global correlation-norm-based objective function can mitigate the dependence on the amplitude information of the nonlinearly smoothed wavefield. Therefore, we have evaluated the use of this objective function when using the nonlinearly smoothed wavefield. The proposed objective function has much larger convexity than the conventional objective functions. We calculate the gradient of the objective function using the adjoint-state technique, which is similar to that of the conventional FWI except for the adjoint source. We progressively reduce the smoothing width applied to the nonlinear wavefield to naturally adopt the multiscale strategy. Using examples on the Marmousi 2 model, we determine that the proposed FWI helps to generate convergent results without the need for low-frequency information.

An overview of full-waveform inversion in exploration geophysics

Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. WCC1-WCC26 ◽  
Author(s):  
J. Virieux ◽  
S. Operto
Keyword(s):  
High Performance ◽  
Waveform Inversion ◽  
Real Data ◽  
Quantitative Information ◽  
Full Waveform Inversion ◽  
Low Frequencies ◽  
Fitting Procedure ◽  
Full Waveform ◽  
Resolution Imaging ◽  
Starting Model

Full-waveform inversion (FWI) is a challenging data-fitting procedure based on full-wavefield modeling to extract quantitative information from seismograms. High-resolution imaging at half the propagated wavelength is expected. Recent advances in high-performance computing and multifold/multicomponent wide-aperture and wide-azimuth acquisitions make 3D acoustic FWI feasible today. Key ingredients of FWI are an efficient forward-modeling engine and a local differential approach, in which the gradient and the Hessian operators are efficiently estimated. Local optimization does not, however, prevent convergence of the misfit function toward local minima because of the limited accuracy of the starting model, the lack of low frequencies, the presence of noise, and the approximate modeling of thewave-physics complexity. Different hierarchical multiscale strategies are designed to mitigate the nonlinearity and ill-posedness of FWI by incorporating progressively shorter wavelengths in the parameter space. Synthetic and real-data case studies address reconstructing various parameters, from [Formula: see text] and [Formula: see text] velocities to density, anisotropy, and attenuation. This review attempts to illuminate the state of the art of FWI. Crucial jumps, however, remain necessary to make it as popular as migration techniques. The challenges can be categorized as (1) building accurate starting models with automatic procedures and/or recording low frequencies, (2) defining new minimization criteria to mitigate the sensitivity of FWI to amplitude errors and increasing the robustness of FWI when multiple parameter classes are estimated, and (3) improving computational efficiency by data-compression techniques to make 3D elastic FWI feasible.

scholarly journals Full-intensity waveform inversion

Geophysics ◽  
2018 ◽  
Vol 83 (6) ◽  
pp. R649-R658 ◽  
Author(s):  
Yike Liu ◽  
Bin He ◽  
Huiyi Lu ◽  
Zhendong Zhang ◽  
Xiao-Bi Xie ◽  
...  
Keyword(s):  
Objective Function ◽  
Waveform Inversion ◽  
Low Frequency ◽  
Velocity Model ◽  
Low Frequencies ◽  
Full Waveform ◽  
Seismic Waveform ◽  
Using Data ◽  
Full Intensity ◽  
Starting Model

Many full-waveform inversion schemes are based on the iterative perturbation theory to fit the observed waveforms. When the observed waveforms lack low frequencies, those schemes may encounter convergence problems due to cycle skipping when the initial velocity model is far from the true model. To mitigate this difficulty, we have developed a new objective function that fits the seismic-waveform intensity, so the dependence of the starting model can be reduced. The waveform intensity is proportional to the square of its amplitude. Forming the intensity using the waveform is a nonlinear operation, which separates the original waveform spectrum into an ultra-low-frequency part and a higher frequency part, even for data that originally do not have low-frequency contents. Therefore, conducting multiscale inversions starting from ultra-low-frequency intensity data can largely avoid the cycle-skipping problem. We formulate the intensity objective function, the minimization process, and the gradient. Using numerical examples, we determine that the proposed method was very promising and could invert for the model using data lacking low-frequency information.

Wave-equation reflection traveltime inversion with dynamic warping and full-waveform inversion

Geophysics ◽  
2013 ◽  
Vol 78 (6) ◽  
pp. R223-R233 ◽  
Author(s):  
Yong Ma ◽  
Dave Hale
Keyword(s):  
Wave Equation ◽  
Waveform Inversion ◽  
Velocity Model ◽  
Full Waveform Inversion ◽  
Local Minima ◽  
Low Frequencies ◽  
Traveltime Inversion ◽  
Velocity Models ◽  
Full Waveform ◽  
Dynamic Warping

In reflection seismology, full-waveform inversion (FWI) can generate high-wavenumber subsurface velocity models but often suffers from an objective function with local minima caused mainly by the absence of low frequencies in seismograms. These local minima cause cycle skipping when the low-wavenumber component in the initial velocity model for FWI is far from the true model. To avoid cycle skipping, we discovered a new wave-equation reflection traveltime inversion (WERTI) to update the low-wavenumber component of the velocity model, while using FWI to only update high-wavenumber details of the model. We implemented the low- and high-wavenumber inversions in an alternating way. In WERTI, we used dynamic image warping (DIW) to estimate the time shifts between recorded data and synthetic data. When compared with correlation-based techniques often used in traveltime estimation, DIW can avoid cycle skipping and estimate the time shifts accurately, even when shifts vary rapidly. Hence, by minimizing traveltime shifts estimated by dynamic warping, WERTI reduces errors in reflection traveltime inversion. Then, conventional FWI uses the low-wavenumber component estimated by WERTI as a new initial model and thereby refines the model with high-wavenumber details. The alternating combination of WERTI and FWI mitigates the velocity-depth ambiguity and can recover subsurface velocities using only high-frequency reflection data.

Building good starting models for full-waveform inversion using adaptive matching filtering misfit

Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. U61-U72 ◽  
Author(s):  
Hejun Zhu ◽  
Sergey Fomel
Keyword(s):  
Least Squares ◽  
Waveform Inversion ◽  
Basin Of Attraction ◽  
Low Frequency ◽  
Quadratic Function ◽  
Full Waveform Inversion ◽  
Local Minima ◽  
Full Waveform ◽  
Gradient Based ◽  
Starting Model

We have proposed a misfit function based on adaptive matching filtering (AMF) to tackle challenges associated with cycle skipping and local minima in full-waveform inversion (FWI). This AMF is designed to measure time-varying phase differences between observations and predictions. Compared with classical least-squares waveform differences, our misfit function behaves as a smooth, quadratic function with a broad basin of attraction. These characters are important because local gradient-based optimization approaches used in most FWI schemes cannot guarantee convergence toward true models if misfit functions include local minima or if the starting model is far away from the global minimum. The 1D and 2D synthetic experiments illustrate the advantages of the proposed misfit function compared with the classical least-squares waveform misfit. Furthermore, we have derived adjoint sources associated with the proposed misfit function and applied them in several 2D time-domain acoustic FWI experiments. Numerical results found that the proposed misfit function can provide good starting models for FWI, particularly when low-frequency signals are absent in recorded data.

Sign in / Sign up

Export Citation Format

Share Document