Application of image-guided full waveform inversion to a 2D ocean-bottom cable data set

2012 ◽  
Author(s):  
Yong Ma ◽  
Jianxin (Jerry) Yuan ◽  
Yunqing Shen ◽  
Bin Gong
2013 ◽  
Author(s):  
Yi Wang ◽  
Kirk Wallace ◽  
Houzhu Zhang ◽  
Alexandre Bertrand ◽  
YunQing Shen

Geophysics ◽  
2017 ◽  
Vol 82 (6) ◽  
pp. R299-R308 ◽  
Author(s):  
Antoine Guitton ◽  
Tariq Alkhalifah

Choosing the right parameterization to describe a transversely isotropic medium with a vertical symmetry axis (VTI) allows us to match the scattering potential of these parameters to the available data in a way that avoids a potential tradeoff and focuses on the parameters to which the data are sensitive. For 2D elastic full-waveform inversion in VTI media of pressure components and for data with a reasonable range of offsets (as with those found in conventional streamer data acquisition systems), assuming that we have a kinematically accurate normal moveout velocity ([Formula: see text]) and anellipticity parameter [Formula: see text] (or horizontal velocity [Formula: see text]) obtained from tomographic methods, a parameterization in terms of horizontal velocity [Formula: see text], [Formula: see text], and [Formula: see text] is preferred to the more conventional parameterization in terms of [Formula: see text], [Formula: see text], and [Formula: see text]. In the [Formula: see text], [Formula: see text], and [Formula: see text] parameterization and for reasonable scattering angles (<[Formula: see text]), [Formula: see text] acts as a “garbage collector” and absorbs most of the amplitude discrepancies between the modeled and observed data, more so when density [Formula: see text] and S-wave velocity [Formula: see text] are not inverted for (a standard practice with streamer data). On the contrary, in the [Formula: see text], [Formula: see text], and [Formula: see text] parameterization, [Formula: see text] is mostly sensitive to large scattering angles, leaving [Formula: see text] exposed to strong leakages from [Formula: see text] mainly. These assertions will be demonstrated on the synthetic Marmousi II as well as a North Sea ocean bottom cable data set, in which inverting for the horizontal velocity rather than the vertical velocity yields more accurate models and migrated images.


Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. WCC149-WCC157 ◽  
Author(s):  
René-Édouard Plessix

With the acquisition of wide-aperture seismic data sets, full-waveform inversion is an attractive method for deriving velocity models. Three-dimensional implementations require an efficient solver for the wave equation. Computing 3D time-harmonic responses with a frequency-domain solver is complicated because a large linear system with negative and positive eigenvalues must be solved. Time-domain schemes are an alternative. Nevertheless, existing frequency-domain iterative solvers with an efficient preconditioner are a viable option when full-waveform inversion is formulated in the frequency domain. An iterative solver with a multigrid preconditioner is competitive because of a high-order spatial discretization. Numerical examples illustrated the efficiency of the iterative solvers. Three dimensional full-waveform inversion was then studied in the context of deep-water ocean-bottom seismometer acquisition. Three dimensional synthetic data inversion results showed the behavior of full-waveform inversion with respect to the initial model and the minimum frequency available in the data set. Results on a 3D real ocean-bottom seismometer data set demonstrated the relevance of full-waveform inversion, especially to image the shallow part of the model.


Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. B201-B210 ◽  
Author(s):  
Abdullah AlTheyab ◽  
G. T. Schuster

We have developed an efficient approach for picking first-break wavefronts on coarsely sampled time slices of 3D shot gathers. Our objective was to compute a smooth initial velocity model for multiscale full-waveform inversion (FWI). Using interactive software, first-break wavefronts were geometrically modeled on time slices with a minimal number of picks. We picked sparse time slices, performed traveltime tomography, and then compared the predicted traveltimes with the data in-between the picked slices. The picking interval was refined with iterations until the errors in traveltime predictions fell within the limits necessary to avoid cycle skipping in early arrivals FWI. This approach was applied to a 3D ocean-bottom-station data set. Our results indicate that wavefront picking has 28% fewer data slices to pick compared with picking traveltimes in shot gathers. In addition, by using sparse time samples for picking, data storage is reduced by 88%, and therefore allows for a faster visualization and quality control of the picks. Our final traveltime tomogram is sufficient as a starting model for early arrival FWI.


Geophysics ◽  
2014 ◽  
Vol 79 (2) ◽  
pp. R63-R77 ◽  
Author(s):  
Denes Vigh ◽  
Kun Jiao ◽  
Dave Watts ◽  
Dong Sun

Recent computational improvements allowed us to simulate elastic wavefields in a 3D manner and undertake the challenge of executing elastic full-waveform inversion (EFWI). The 3D SEG/EAGE overthrust synthetic data were used to run the initial tests, which included using all three components for the simulation. The inversion targeted two regions: the channel system and the overthrusted zone, which proved the effectiveness of EFWI to delineate geology in terms of [Formula: see text] and [Formula: see text] velocity fields. For the field data experiment to demonstrate the technologies, we elected to use a Gulf of Mexico ocean bottom cable data set, which allowed us to take advantage of relatively large offsets along with the 4C acquisition. The input data were minimally processed mostly through noise removal, and the initial model was a Gaussian smoothed version of grid tomography output, which is done by a prestack migrated gather flattening process. During EWFI, a multiscale approach was followed to ensure convergence, and the early stages of the [Formula: see text]/[Formula: see text] ratio were constrained by the mud rock-line ratio. When the last sets of inversions were executed, this constraint was eliminated to ensure the simultaneous update of the [Formula: see text] and [Formula: see text] velocity fields. The density was kept constant to keep the inversion at a simple level, which allowed us to draw essential conclusions. The velocity fields were validated through an imaging algorithm of the elastic reverse time migration, and the imaging shows clear structural improvements when inputting the inverted velocities in conjunction with the measurements. If full-waveform inversion can provide multiple earth parameters, the user can use the process to detect gas zones along with sand and shale content of the subsurface, which will further assist the drilling decisions. We achieved this by simulating the earth more accurately with the elastic wave propagation in the algorithms.


Geophysics ◽  
2016 ◽  
Vol 81 (4) ◽  
pp. U25-U38 ◽  
Author(s):  
Nuno V. da Silva ◽  
Andrew Ratcliffe ◽  
Vetle Vinje ◽  
Graham Conroy

Parameterization lies at the center of anisotropic full-waveform inversion (FWI) with multiparameter updates. This is because FWI aims to update the long and short wavelengths of the perturbations. Thus, it is important that the parameterization accommodates this. Recently, there has been an intensive effort to determine the optimal parameterization, centering the fundamental discussion mainly on the analysis of radiation patterns for each one of these parameterizations, and aiming to determine which is best suited for multiparameter inversion. We have developed a new parameterization in the scope of FWI, based on the concept of kinematically equivalent media, as originally proposed in other areas of seismic data analysis. Our analysis is also based on radiation patterns, as well as the relation between the perturbation of this set of parameters and perturbation in traveltime. The radiation pattern reveals that this parameterization combines some of the characteristics of parameterizations with one velocity and two Thomsen’s parameters and parameterizations using two velocities and one Thomsen’s parameter. The study of perturbation of traveltime with perturbation of model parameters shows that the new parameterization is less ambiguous when relating these quantities in comparison with other more commonly used parameterizations. We have concluded that our new parameterization is well-suited for inverting diving waves, which are of paramount importance to carry out practical FWI successfully. We have demonstrated that the new parameterization produces good inversion results with synthetic and real data examples. In the latter case of the real data example from the Central North Sea, the inverted models show good agreement with the geologic structures, leading to an improvement of the seismic image and flatness of the common image gathers.


2010 ◽  
Author(s):  
Yong Ma ◽  
Dave Hale ◽  
Zhaobo (Joe) Meng ◽  
Bin Gong

Sign in / Sign up

Export Citation Format

Share Document