Full waveform inversion method using envelope objective function without low frequency data

2014 ◽  
Vol 109 ◽  
pp. 36-46 ◽  
Author(s):  
Benxin Chi ◽  
Liangguo Dong ◽  
Yuzhu Liu
Geophysics ◽  
2021 ◽  
pp. 1-54
Author(s):  
Milad Bader ◽  
Robert G. Clapp ◽  
Biondo Biondi

Low-frequency data below 5 Hz are essential to the convergence of full-waveform inversion towards a useful solution. They help build the velocity model low wavenumbers and reduce the risk of cycle-skipping. In marine environments, low-frequency data are characterized by a low signal-to-noise ratio and can lead to erroneous models when inverted, especially if the noise contains coherent components. Often field data are high-pass filtered before any processing step, sacrificing weak but essential signal for full-waveform inversion. We propose to denoise the low-frequency data using prediction-error filters that we estimate from a high-frequency component with a high signal-to-noise ratio. The constructed filter captures the multi-dimensional spectrum of the high-frequency signal. We expand the filter's axes in the time-space domain to compress its spectrum towards the low frequencies and wavenumbers. The expanded filter becomes a predictor of the target low-frequency signal, and we incorporate it in a minimization scheme to attenuate noise. To account for data non-stationarity while retaining the simplicity of stationary filters, we divide the data into non-overlapping patches and linearly interpolate stationary filters at each data sample. We apply our method to synthetic stationary and non-stationary data, and we show it improves the full-waveform inversion results initialized at 2.5 Hz using the Marmousi model. We also demonstrate that the denoising attenuates non-stationary shear energy recorded by the vertical component of ocean-bottom nodes.


Geophysics ◽  
2020 ◽  
Vol 85 (6) ◽  
pp. R553-R563
Author(s):  
Sagar Singh ◽  
Ilya Tsvankin ◽  
Ehsan Zabihi Naeini

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.


Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. R77-R88 ◽  
Author(s):  
Yunseok Choi ◽  
Tariq Alkhalifah

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.


Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. R339-R348 ◽  
Author(s):  
Yunyue Elita Li ◽  
Laurent Demanet

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.


Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. R117-R127 ◽  
Author(s):  
Yuanyuan Li ◽  
Yunseok Choi ◽  
Tariq Alkhalifah ◽  
Zhenchun Li ◽  
Kai Zhang

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.


2014 ◽  
Vol 1 (2) ◽  
pp. 1757-1802
Author(s):  
C. Huang ◽  
L. Dong ◽  
Y. Liu ◽  
B. Chi

Abstract. Low frequency is a key issue to reduce the nonlinearity of elastic full waveform inversion. Hence, the lack of low frequency in recorded seismic data is one of the most challenging problems in elastic full waveform inversion. Theoretical derivations and numerical analysis are presented in this paper to show that envelope operator can retrieve strong low frequency modulation signal demodulated in multicomponent data, no matter what the frequency bands of the data is. With the benefit of such low frequency information, we use elastic envelope of multicomponent data to construct the objective function and present an elastic envelope inversion method to recover the long-wavelength components of the subsurface model, especially for the S-wave velocity model. Numerical tests using synthetic data for the Marmousi-II model prove the effectiveness of the proposed elastic envelope inversion method, especially when low frequency is missing in multicomponent data and when initial model is far from the true model. The elastic envelope can reduce the nonlinearity of inversion and can provide an excellent starting model.


Geophysics ◽  
2020 ◽  
Vol 85 (6) ◽  
pp. A37-A43
Author(s):  
Jinwei Fang ◽  
Hui Zhou ◽  
Yunyue Elita Li ◽  
Qingchen Zhang ◽  
Lingqian Wang ◽  
...  

The lack of low-frequency signals in seismic data makes the full-waveform inversion (FWI) procedure easily fall into local minima leading to unreliable results. To reconstruct the missing low-frequency signals more accurately and effectively, we have developed a data-driven low-frequency recovery method based on deep learning from high-frequency signals. In our method, we develop the idea of using a basic data patch of seismic data to build a local data-driven mapping in low-frequency recovery. Energy balancing and data patches are used to prepare high- and low-frequency data for training a convolutional neural network (CNN) to establish the relationship between the high- and low-frequency data pairs. The trained CNN then can be used to predict low-frequency data from high-frequency data. Our CNN was trained on the Marmousi model and tested on the overthrust model, as well as field data. The synthetic experimental results reveal that the predicted low-frequency data match the true low-frequency data very well in the time and frequency domains, and the field results show the successfully extended low-frequency spectra. Furthermore, two FWI tests using the predicted data demonstrate that our approach can reliably recover the low-frequency data.


Sign in / Sign up

Export Citation Format

Share Document