High frequency full waveform inversion as an interpretation solution

2019 ◽  
Vol 59 (2) ◽  
pp. 904
Author(s):  
Laurence Letki ◽  
Matt Lamont ◽  
Troy Thompson

The amount of available data to help us characterise the subsurface is ever increasing. Large seismic surveys, long offsets, multi- and full-azimuth datasets, including 3D and 4D, marine, ocean-bottom nodes and extremely high fold land surveys, are now common. In parallel, computing power is also increasing and, in combination with better data, this enables us to develop better tools and to use better physics to build models of the subsurface. Wave-equation based techniques, such as full waveform inversion (FWI), have therefore become a lot more practical. FWI uses the entire wavefield, including refractions and reflections, primaries and multiples, to generate a refined, high resolution Earth model. This technique is now commonly used at lower frequencies (up to 12 Hz) to derive more accurate models for improved seismic imaging and reduced depth conversion uncertainty. By including higher frequencies in FWI, we can attempt to resolve for finer and finer details. FWI models using the entire bandwidth of the seismic data constitute an interpretation product in itself, with applications in both structural interpretation and reservoir characterisation. Incorporating more physics within the FWI implementation, combined with modern supercomputer facilities, promises to increase the focus on very high frequency FWI in the coming years. In this paper, through a series of field examples, we illustrate the applications and rewards of high frequency FWI: from improved imaging, improved quantitative interpretation and depth conversion to a direct interpretation of the FWI models.

2021 ◽  
Vol 40 (5) ◽  
pp. 335-341
Author(s):  
Denes Vigh ◽  
Xin Cheng ◽  
Kun Jiao ◽  
Wei Kang ◽  
Nolan Brand

Full-waveform inversion (FWI) is a high-resolution model-building technique that uses the entire recorded seismic data content to build the earth model. Conventional FWI usually utilizes diving and refracted waves to update the low-wavenumber components of the velocity model. However, updates are often depth limited due to the limited offset range of the acquisition design. To extend conventional FWI beyond the limits imposed by using only transmitted energy, we must utilize the full acquired wavefield. Analyzing FWI kernels for a given geology and acquisition geometry can provide information on how to optimize the acquisition so that FWI is able to update the velocity model for targets as deep as basement level. Recent long-offset ocean-bottom node acquisition helped FWI succeed, but we would also like to be able to utilize the shorter-offset data from wide-azimuth data acquisitions to improve imaging of these data sets by developing the velocity field with FWI. FWI models are heading toward higher and higher wavenumbers, which allows us to extract pseudoreflectivity directly from the developed velocity model built with the acoustic full wavefield. This is an extremely early start to obtaining a depth image that one would usually produce in much later processing stages.


2020 ◽  
Author(s):  
Suyang Chen ◽  
Netifetu Usman-Lecky ◽  
Tamunoemi Okougbo ◽  
Emmanuel Saragoussi ◽  
Jean-Patrick Mascomère ◽  
...  

2019 ◽  
Vol 38 (3) ◽  
pp. 204-213 ◽  
Author(s):  
Ping Wang ◽  
Zhigang Zhang ◽  
Jiawei Mei ◽  
Feng Lin ◽  
Rongxin Huang

Full-waveform inversion (FWI), proposed by Lailly and Tarantola in the 1980s, is considered to be the most promising data-driven tool for automatically building velocity models. Many successful examples have been reported using FWI to update shallow sediments, gas pockets, and mud volcanoes. However, successful applications of FWI to update salt structures had almost only been seen on synthetic data until recent progress at the Atlantis Field in the Gulf of Mexico. We revisited some aspects of FWI algorithms to minimize cycle-skipping and amplitude discrepancy issues and derived an FWI algorithm that is able to build complex salt velocity models. We applied this algorithm to a variety of data sets, including wide-azimuth and full-azimuth (FAZ) streamer data as well as ocean-bottom-node data, with different geologic settings in order to demonstrate the effectiveness of the method for salt velocity updates and to examine some fundamentals of the salt problem. We observed that, in multiple cases, salt velocity models from this FWI algorithm produced subsalt images of superior quality. We demonstrate with one FAZ streamer data example in Keathley Canyon that we do not necessarily need very high frequency in FWI for subsalt imaging purposes. Based on this observation, we envision that sparse node for velocity acquisition may provide appropriate data to handle large and complex salt bodies with FWI. We believe the combination of advanced FWI algorithms and appropriate data acquisition will bring a step change to subsalt imaging.


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 ◽  
2021 ◽  
pp. 1-52
Author(s):  
Yuzhu Liu ◽  
Xinquan Huang ◽  
Jizhong Yang ◽  
Xueyi Liu ◽  
Bin Li ◽  
...  

Thin sand-mud-coal interbedded layers and multiples caused by shallow water pose great challenges to conventional 3D multi-channel seismic techniques used to detect the deeply buried reservoirs in the Qiuyue field. In 2017, a dense ocean-bottom seismometer (OBS) acquisition program acquired a four-component dataset in East China Sea. To delineate the deep reservoir structures in the Qiuyue field, we applied a full-waveform inversion (FWI) workflow to this dense four-component OBS dataset. After preprocessing, including receiver geometry correction, moveout correction, component rotation, and energy transformation from 3D to 2D, a preconditioned first-arrival traveltime tomography based on an improved scattering integral algorithm is applied to construct an initial P-wave velocity model. To eliminate the influence of the wavelet estimation process, a convolutional-wavefield-based objective function for the preprocessed hydrophone component is used during acoustic FWI. By inverting the waveforms associated with early arrivals, a relatively high-resolution underground P-wave velocity model is obtained, with updates at 2.0 km and 4.7 km depth. Initial S-wave velocity and density models are then constructed based on their prior relationships to the P-wave velocity, accompanied by a reciprocal source-independent elastic full-waveform inversion to refine both velocity models. Compared to a traditional workflow, guided by stacking velocity analysis or migration velocity analysis, and using only the pressure component or other single-component, the workflow presented in this study represents a good approach for inverting the four-component OBS dataset to characterize sub-seafloor velocity structures.


Sign in / Sign up

Export Citation Format

Share Document