Seismic SH Full Waveform Inversion: A tool for high-resolution near-surface characterization

Author(s):  
D. Köhn ◽  
M. Zolchow ◽  
R. Mecking ◽  
D. Wilken ◽  
T. Wunderlich ◽  
...  
Geophysics ◽  
2020 ◽  
pp. 1-57
Author(s):  
Daniele Colombo ◽  
Ernesto Sandoval ◽  
Diego Rovetta ◽  
Apostolos Kontakis

Land seismic velocity modeling is a difficult task largely related to the description of the near surface complexities. Full waveform inversion is the method of choice for achieving high-resolution velocity mapping but its application to land seismic data faces difficulties related to complex physics, unknown and spatially varying source signatures, and low signal-to-noise ratio in the data. Large parameter variations occur in the near surface at various scales causing severe kinematic and dynamic distortions of the recorded wavefield. Some of the parameters can be incorporated in the inversion model while others, due to sub-resolution dimensions or unmodeled physics, need to be corrected through data preconditioning; a topic not well described for land data full waveform inversion applications. We have developed novel algorithms and workflows for surface-consistent data preconditioning utilizing the transmitted portion of the wavefield, signal-to-noise enhancement by generation of CMP-based virtual super shot gathers, and robust 1.5D Laplace-Fourier full waveform inversion. Our surface-consistent scheme solves residual kinematic corrections and amplitude anomalies via scalar compensation or deconvolution of the near surface response. Signal-to-noise enhancement is obtained through the statistical evaluation of volumetric prestack responses at the CMP position, or virtual super (shot) gathers. These are inverted via a novel 1.5D acoustic Laplace-Fourier full waveform inversion scheme using the Helmholtz wave equation and Hankel domain forward modeling. Inversion is performed with nonlinear conjugate gradients. The method is applied to a complex structure-controlled wadi area exhibiting faults, dissolution, collapse, and subsidence where the high resolution FWI velocity modeling helps clarifying the geological interpretation. The developed algorithms and automated workflows provide an effective solution for massive full waveform inversion of land seismic data that can be embedded in typical near surface velocity analysis procedures.


2019 ◽  
Vol 7 (1) ◽  
pp. T141-T154 ◽  
Author(s):  
Md. Iftekhar Alam

Seismic imaging of the shallow subsurface (approximately 5 m) can be very challenging when reflections are absent and the data are dominated by ground roll. I analyzed the transmission coda to produce fine-scale, interpretable vertical and horizontal component seismic velocity ([Formula: see text] and [Formula: see text]) models using full-waveform inversion (FWI). Application of FWI is tested through imaging two buried targets. The first target is a pair of well-documented utility pipes with known diameters (0.8 m) and burial depths (approximately 1.5 m). The second target is a poorly documented former location of the pipe(s), which is now a backfilled void. Data are acquired along a 23 m 2D profile using a static array with single-component vertical and horizontal geophones. Our results indicate considerable velocity updates in the [Formula: see text] and [Formula: see text] models across the pipes and backfill. The pipes appear as negative velocity updates in the final inverted [Formula: see text] and [Formula: see text] models, whereas the backfilled area represents negative and positive velocity updates in the [Formula: see text] and [Formula: see text] models, respectively. Variations of the polarities in the inverted models ([Formula: see text] and [Formula: see text]) across the backfill can be indicative of the medium, which respond differently to the vertical and horizontal component seismic waves. The attenuation models show a general decreasing trend with increasing depth. Therefore, simultaneous applications of vertical ([Formula: see text]) and horizontal ([Formula: see text]) component seismic velocity modeling can be an effective tool to understand the subsurface medium in near-surface characterization.


2017 ◽  
Vol 5 (4) ◽  
pp. SR23-SR33 ◽  
Author(s):  
Xin Cheng ◽  
Kun Jiao ◽  
Dong Sun ◽  
Zhen Xu ◽  
Denes Vigh ◽  
...  

Over the past decade, acoustic full-waveform inversion (FWI) has become one of the standard methods in the industry to construct high-resolution velocity fields from the seismic data acquired. While most of the successful applications are for marine acquisition data with rich low-frequency diving or postcritical waves at large offsets, the application of acoustic FWI on land data remains a challenging topic. Land acoustic FWI application faces many severe difficulties, such as the presence of strong elastic effects, large near-surface velocity contrast, and heterogeneous, topography variations, etc. In addition, it is well-known that low-frequency transmitted seismic energy is crucial for the success of FWI to overcome sensitivity to starting velocity fields; unfortunately, those are the parts of the data that suffer the most from a low signal-to-noise ratio (S/N) in land acquisition. We have developed an acoustic FWI application on a land data set from North Kuwait, and demonstrated our solutions to mitigate some of the challenges posed by land data. More specifically, we have developed a semblance-based high-resolution Radon (HR-Radon) inversion approach to enhance the S/N of the low-frequency part of the FWI input data and to ultimately improve the convergence of the land FWI workflow. To mitigate the impact of elastic effects, we included only the diving and postcritical early arrivals in the waveform inversion. Our results show that, with the aid of HR-Radon preconditioning and a carefully designed workflow, acoustic FWI has the ability to derive a reliable high-resolution near-surface model that could not be otherwise recovered through traditional tomographic methods.


2021 ◽  
Author(s):  
Rahul Dixit ◽  
Pavel Vasilyev ◽  
Ivica Mihaljevic ◽  
Michelle Tham ◽  
Denes Vigh ◽  
...  

Abstract Full-waveform inversion (FWI) has become a well-established method for obtaining a detailed earth model suitable for improved imaging, near-surface characterization and pore-pressure prediction. FWI for onshore data has always been challenging and has seen limited application (Vigh et al, 2018). It requires a dedicated data processing approach related to the lower signal-to-noise ratio, accounting for variable topography and complex near-surface related effects. During the past few years, ADNOC has been acquiring and processing one of the world's largest combined 3D onshore and offshore seismic surveys in the Emirate of Abu Dhabi. The modern acquisition parameters that were implemented enabled the acquisition of broadband onshore seismic data rich in low frequencies that could benefit the initial stages of the FWI workflow. Sand dunes and sabkha layers at the surface, and high-velocity carbonate and dolomite layers in the subsurface pose a significant challenge for near-surface modeling in the UAE. The purpose of this work is to evaluate FWI application onshore UAE for near-surface characterization. We will compare the FWI results with conventional approaches for the near-surface model building that has been used routinely on land datasets in UAE, such as data-driven image-based statics (DIBS, Zarubov et al, 2019). One of the main challenges is data preconditioning, as onshore seismic data typically exhibits high levels of noise. It is imperative to denoise gathers sufficiently prior to the FWI process. A well sonic velocity function with large smoothing was used to build the starting velocity model for FWI. The process aims to minimize the least-squared difference between predicted and observed seismic responses by means of updating the model on which the prediction is based. As the predicted and seismic responses are functions of model parameters as well as source signature, a good estimate of the source wavelet is important for update and convergence in FWI. During this FWI work, source wavelet inversion was done as a separate step and used in subsequent FWI passes. FWI inversion started with adjustive FWI (Kun et al, 2015) on lower frequencies, moving to higher frequencies where both adjustive and least square objective functions were used. We will further show assessment of the anisotropy, initial conditions, usage of geological constraints, and comparisons to the conventional solutions. A comparison of results shows that FWI has successfully added velocity details to the near-surface model that follow the geological trend and conforms to well information while producing a plausible static solution. We have demonstrated the application of FWI onshore UAE for near-surface modeling. Although turnaround time (TAT) has increased compared to the conventional approach, the learning that was gained during this trial will decrease TAT for the future FWI work.


Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE101-VE117 ◽  
Author(s):  
Hafedh Ben-Hadj-Ali ◽  
Stéphane Operto ◽  
Jean Virieux

We assessed 3D frequency-domain (FD) acoustic full-waveform inversion (FWI) data as a tool to develop high-resolution velocity models from low-frequency global-offset data. The inverse problem was posed as a classic least-squares optimization problem solved with a steepest-descent method. Inversion was applied to a few discrete frequencies, allowing management of a limited subset of the 3D data volume. The forward problem was solved with a finite-difference frequency-domain method based on a massively parallel direct solver, allowing efficient multiple-shot simulations. The inversion code was fully parallelized for distributed-memory platforms, taking advantage of a domain decomposition of the modeled wavefields performed by the direct solver. After validation on simple synthetic tests, FWI was applied to two targets (channel and thrust system) of the 3D SEG/EAGE overthrust model, corresponding to 3D domains of [Formula: see text] and [Formula: see text], respectively. The maximum inverted frequencies are 15 and [Formula: see text] for the two applications. A maximum of 30 dual-core biprocessor nodes with [Formula: see text] of shared memory per node were used for the second target. The main structures were imaged successfully at a resolution scale consistent with the inverted frequencies. Our study confirms the feasibility of 3D frequency-domain FWI of global-offset data on large distributed-memory platforms to develop high-resolution velocity models. These high-velocity models may provide accurate macromodels for wave-equation prestack depth migration.


Geophysics ◽  
2021 ◽  
pp. 1-37
Author(s):  
Xinhai Hu ◽  
Wei Guoqi ◽  
Jianyong Song ◽  
Zhifang Yang ◽  
Minghui Lu ◽  
...  

Coupling factors of sources and receivers vary dramatically due to the strong heterogeneity of near surface, which are as important as the model parameters for the inversion success. We propose a full waveform inversion (FWI) scheme that corrects for variable coupling factors while updating the model parameter. A linear inversion is embedded into the scheme to estimate the source and receiver factors and compute the amplitude weights according to the acquisition geometry. After the weights are introduced in the objective function, the inversion falls into the category of separable nonlinear least-squares problems. Hence, we could use the variable projection technique widely used in source estimation problem to invert the model parameter without the knowledge of source and receiver factors. The efficacy of the inversion scheme is demonstrated with two synthetic examples and one real data test.


Sign in / Sign up

Export Citation Format

Share Document