Recover certain low-frequency information for full waveform inversion

Author(s):  
Xiao-Bi Xie
2019 ◽  
Vol 16 (6) ◽  
pp. 1017-1031 ◽  
Author(s):  
Yong Hu ◽  
Liguo Han ◽  
Rushan Wu ◽  
Yongzhong Xu

Abstract Full Waveform Inversion (FWI) is based on the least squares algorithm to minimize the difference between the synthetic and observed data, which is a promising technique for high-resolution velocity inversion. However, the FWI method is characterized by strong model dependence, because the ultra-low-frequency components in the field seismic data are usually not available. In this work, to reduce the model dependence of the FWI method, we introduce a Weighted Local Correlation-phase based FWI method (WLCFWI), which emphasizes the correlation phase between the synthetic and observed data in the time-frequency domain. The local correlation-phase misfit function combines the advantages of phase and normalized correlation function, and has an enormous potential for reducing the model dependence and improving FWI results. Besides, in the correlation-phase misfit function, the amplitude information is treated as a weighting factor, which emphasizes the phase similarity between synthetic and observed data. Numerical examples and the analysis of the misfit function show that the WLCFWI method has a strong ability to reduce model dependence, even if the seismic data are devoid of low-frequency components and contain strong Gaussian noise.


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.


2021 ◽  
Author(s):  
Yao Liu ◽  
Baodi Liu ◽  
Jianping Huang ◽  
Jun Wang ◽  
Honglong Chen ◽  
...  

2021 ◽  
Author(s):  
J. Pruessmann ◽  
G. Eisenberg Klein ◽  
E. Schuenemann ◽  
E. Verschuur ◽  
S. Qu

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 ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. R409-R423
Author(s):  
Polina Zheglova ◽  
Alison Malcolm

Vector-acoustic full-waveform inversion (VAFWI) directly inverts vector-acoustic (VA) data, which consist of pressure and particle displacement components, at the cost of conventional acoustic full-waveform inversion (FWI). VA data contain information about the direction of arrival of the recorded seismic waves. In VAFWI, this directional information is taken into account by introducing an appropriate data weighting. With this weighting, in the geometry of a marine seismic experiment, the VAFWI adjoint calculation approximates inverse wavefield extrapolation, resulting in the natural separation of up- and downgoing recorded waves. If the free-surface effects are modeled during the inversion, the wave separation leads to (1) suppression of surface-related artifacts, (2) constructive interference of receiver ghosts with their primaries leading to preservation of the low-frequency content in the adjoint fields, and (3) compensation for insufficient spatial wavefield sampling on the receiver side. The horizontal displacement component helps interpolate the missing data. Synthetic examples demonstrate that for undersampled data, VAFWI consistently recovers the subsurface properties with higher resolution and fewer artifacts than conventional FWI.


2019 ◽  
Vol 7 (2) ◽  
pp. SB43-SB52 ◽  
Author(s):  
Adriano Gomes ◽  
Joe Peterson ◽  
Serife Bitlis ◽  
Chengliang Fan ◽  
Robert Buehring

Inverting for salt geometry using full-waveform inversion (FWI) is a challenging task, mostly due to the lack of extremely low-frequency signal in the seismic data, the limited penetration depth of diving waves using typical acquisition offsets, and the difficulty in correctly modeling the amplitude (and kinematics) of reflection events associated with the salt boundary. However, recent advances in reflection FWI (RFWI) have allowed it to use deep reflection data, beyond the diving-wave limit, by extracting the tomographic term of the FWI reflection update, the so-called rabbit ears. Though lacking the resolution to fully resolve salt geometry, we can use RFWI updates as a guide for refinements in the salt interpretation, adding a partially data-driven element to salt velocity model building. In addition, we can use RFWI to update sediment velocities in complex regions surrounding salt, where ray-based approaches typically struggle. In reality, separating the effects of sediment velocity errors from salt geometry errors is not straightforward in many locations. Therefore, iterations of RFWI plus salt scenario tests may be necessary. Although it is still not the fully automatic method that has been envisioned for FWI, this combined approach can bring significant improvement to the subsalt image, as we examine on field data examples from the Gulf of Mexico.


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.


Sign in / Sign up

Export Citation Format

Share Document