scholarly journals Time-domain full-waveform inversion of exponentially damped wavefield using the deconvolution-based objective function

Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. R77-R88 ◽  
Author(s):  
Yunseok Choi ◽  
Tariq Alkhalifah
Keyword(s):  
Objective Function ◽  
Frequency Band ◽  
Waveform Inversion ◽  
Low Frequency ◽  
Full Waveform Inversion ◽  
Low Frequencies ◽  
Full Waveform ◽  
Adjoint State ◽  
Recorded Data ◽  
Adjoint State Method

Full-waveform inversion (FWI) suffers from the cycle-skipping problem when the available frequency-band of data is not low enough. We have applied an exponential damping to the data to generate artificial low frequencies, which helps FWI to avoid cycle skipping. In this case, the least-squares misfit function does not properly deal with the exponentially damped wavefield in FWI because the amplitude of traces decays almost exponentially with increasing offset in a damped wavefield. Thus, we use a deconvolution-based objective function for FWI of the exponentially damped wavefield. The deconvolution filter includes inherently a normalization between the modeled and observed data; thus, it can address the unbalanced amplitude of a damped wavefield. We specifically normalize the modeled data with the observed data in the frequency-domain to estimate the deconvolution filter and selectively choose a frequency-band for normalization that mainly includes the artificial low frequencies. We calculate the gradient of the objective function using the adjoint-state method. The synthetic and benchmark data examples indicate that our FWI algorithm generates a convergent long-wavelength structure without low-frequency information in the recorded data.

scholarly journals Full-waveform inversion, Part 3: Optimization

2018 ◽  
Vol 37 (2) ◽  
pp. 142-145 ◽  
Author(s):  
Philipp Witte ◽  
Mathias Louboutin ◽  
Keegan Lensink ◽  
Michael Lange ◽  
Navjot Kukreja ◽  
...  
Keyword(s):  
Objective Function ◽  
Wave Equations ◽  
Waveform Inversion ◽  
Full Waveform Inversion ◽  
Acoustic Wave Equation ◽  
Automated Generation ◽  
Domain Specific ◽  
Full Waveform ◽  
Adjoint State ◽  
Adjoint State Method

This tutorial is the third part of a full-waveform inversion (FWI) tutorial series with a step-by-step walkthrough of setting up forward and adjoint wave equations and building a basic FWI inversion framework. For discretizing and solving wave equations, we use Devito ( http://www.opesci.org/devito-public ), a Python-based domain-specific language for automated generation of finite-difference code ( Lange et al., 2016 ). The first two parts of this tutorial ( Louboutin et al., 2017 , 2018 ) demonstrated how to solve the acoustic wave equation for modeling seismic shot records and how to compute the gradient of the FWI objective function using the adjoint-state method. With these two key ingredients, we will now build an inversion framework that can be used to minimize the FWI least-squares objective function.

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.

Decoupled Fréchet kernels based on a fractional viscoacoustic wave equation

Geophysics ◽  
2021 ◽  
pp. 1-42
Author(s):  
Guangchi Xing ◽  
Tieyuan Zhu
Keyword(s):  
Wave Equation ◽  
Fractional Laplacian ◽  
Waveform Inversion ◽  
Full Waveform Inversion ◽  
Physical Processes ◽  
Full Waveform ◽  
Adjoint State ◽  
The Finite Difference Method ◽  
Adjoint State Method ◽  
Laplacian Operators

We formulate the Fréchet kernel computation using the adjoint-state method based on a fractional viscoacoustic wave equation. We first numerically prove that both the 1/2- and the 3/2-order fractional Laplacian operators are self-adjoint. Using this property, we show that the adjoint wave propagator preserves the dispersion and compensates the amplitude, while the time-reversed adjoint wave propagator behaves identically as the forward propagator with the same dispersion and dissipation characters. Without introducing rheological mechanisms, this formulation adopts an explicit Q parameterization, which avoids the implicit Q in the conventional viscoacoustic/viscoelastic full waveform inversion ( Q-FWI). In addition, because of the decoupling of operators in the wave equation, the viscoacoustic Fréchet kernel is separated into three distinct contributions with clear physical meanings: lossless propagation, dispersion, and dissipation. We find that the lossless propagation kernel dominates the velocity kernel, while the dissipation kernel dominates the attenuation kernel over the dispersion kernel. After validating the Fréchet kernels using the finite-difference method, we conduct a numerical example to demonstrate the capability of the kernels to characterize both velocity and attenuation anomalies. The kernels of different misfit measurements are presented to investigate their different sensitivities. Our results suggest that rather than the traveltime, the amplitude and the waveform kernels are more suitable to capture attenuation anomalies. These kernels lay the foundation for the multiparameter inversion with the fractional formulation, and the decoupled nature of them promotes our understanding of the significance of different physical processes in the Q-FWI.

Full-waveform inversion with borehole constraints for elastic VTI media

Geophysics ◽  
2020 ◽  
Vol 85 (6) ◽  
pp. R553-R563
Author(s):  
Sagar Singh ◽  
Ilya Tsvankin ◽  
Ehsan Zabihi Naeini
Keyword(s):  
Objective Function ◽  
Reservoir Characterization ◽  
Waveform Inversion ◽  
Low Frequency ◽  
Rock Physics ◽  
Anisotropic Media ◽  
Full Waveform Inversion ◽  
Full Waveform ◽  
Trade Offs ◽  
Vti Media

The nonlinearity of full-waveform inversion (FWI) and parameter trade-offs can prevent convergence toward the actual model, especially for elastic anisotropic media. The problems with parameter updating become particularly severe if ultra-low-frequency seismic data are unavailable, and the initial model is not sufficiently accurate. We introduce a robust way to constrain the inversion workflow using borehole information obtained from well logs. These constraints are included in the form of rock-physics relationships for different geologic facies (e.g., shale, sand, salt, and limestone). We develop a multiscale FWI algorithm for transversely isotropic media with a vertical symmetry axis (VTI media) that incorporates facies information through a regularization term in the objective function. That term is updated during the inversion by using the models obtained at the previous inversion stage. To account for lateral heterogeneity between sparse borehole locations, we use an image-guided smoothing algorithm. Numerical testing for structurally complex anisotropic media demonstrates that the facies-based constraints may ensure the convergence of the objective function towards the global minimum in the absence of ultra-low-frequency data and for simple (even 1D) initial models. We test the algorithm on clean data and on surface records contaminated by Gaussian noise. The algorithm also produces a high-resolution facies model, which should be instrumental in reservoir characterization.

Tackling cycle skipping in full-waveform inversion with intermediate data

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. R411-R427 ◽  
Author(s):  
Gang Yao ◽  
Nuno V. da Silva ◽  
Michael Warner ◽  
Di Wu ◽  
Chenhao Yang
Keyword(s):  
Objective Function ◽  
Waveform Inversion ◽  
Synthetic Data ◽  
Velocity Model ◽  
Full Waveform Inversion ◽  
Data Set ◽  
Intermediate Data ◽  
Full Waveform ◽  
L2 Norm ◽  
Recorded Data

Full-waveform inversion (FWI) is a promising technique for recovering the earth models for exploration geophysics and global seismology. FWI is generally formulated as the minimization of an objective function, defined as the L2-norm of the data residuals. The nonconvex nature of this objective function is one of the main obstacles for the successful application of FWI. A key manifestation of this nonconvexity is cycle skipping, which happens if the predicted data are more than half a cycle away from the recorded data. We have developed the concept of intermediate data for tackling cycle skipping. This intermediate data set is created to sit between predicted and recorded data, and it is less than half a cycle away from the predicted data. Inverting the intermediate data rather than the cycle-skipped recorded data can then circumvent cycle skipping. We applied this concept to invert cycle-skipped first arrivals. First, we picked up the first breaks of the predicted data and the recorded data. Second, we linearly scaled down the time difference between the two first breaks of each shot into a series of time shifts, the maximum of which was less than half a cycle, for each trace in this shot. Third, we moved the predicted data with the corresponding time shifts to create the intermediate data. Finally, we inverted the intermediate data rather than the recorded data. Because the intermediate data are not cycle-skipped and contain the traveltime information of the recorded data, FWI with intermediate data updates the background velocity model in the correct direction. Thus, it produces a background velocity model accurate enough for carrying out conventional FWI to rebuild the intermediate- and short-wavelength components of the velocity model. Our numerical examples using synthetic data validate the intermediate-data concept for tackling cycle skipping and demonstrate its effectiveness for the application to first arrivals.

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.

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 Extrapolated full-waveform inversion with deep learning

Geophysics ◽  
2020 ◽  
Vol 85 (3) ◽  
pp. R275-R288 ◽  
Author(s):  
Hongyu Sun ◽  
Laurent Demanet
Keyword(s):  
Neural Network ◽  
Numerical Experiments ◽  
Waveform Inversion ◽  
Low Frequency ◽  
Full Waveform Inversion ◽  
Bandwidth Extension ◽  
Low Frequencies ◽  
Full Waveform ◽  
The Neural Network ◽  
Band Limited

The lack of low-frequency information and a good initial model can seriously affect the success of full-waveform inversion (FWI), due to the inherent cycle skipping problem. Computational low-frequency extrapolation is in principle the most direct way to address this issue. By considering bandwidth extension as a regression problem in machine learning, we have adopted an architecture of convolutional neural network (CNN) to automatically extrapolate the missing low frequencies. The band-limited recordings are the inputs of the CNN, and, in our numerical experiments, a neural network trained from enough samples can predict a reasonable approximation to the seismograms in the unobserved low-frequency band, in phase and in amplitude. The numerical experiments considered are set up on simulated P-wave data. In extrapolated FWI (EFWI), the low-wavenumber components of the model are determined from the extrapolated low frequencies, before proceeding with a frequency sweep of the band-limited data. The introduced deep-learning method of low-frequency extrapolation shows adequate generalizability for the initialization step of EFWI. Numerical examples show that the neural network trained on several submodels of the Marmousi model is able to predict the low frequencies for the BP 2004 benchmark model. Additionally, the neural network can robustly process seismic data with uncertainties due to the existence of random noise, a poorly known source wavelet, and a different finite-difference scheme in the forward modeling operator. Finally, this approach is not subject to strong assumptions on signals or velocity models of other methods for bandwidth extension and seems to offer a tantalizing solution to the problem of properly initializing FWI.

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.

Sign in / Sign up

Export Citation Format

Share Document