Search

Full waveform inversion (FWI) aims to reconstruct high-resolution subsurface models from the full wavefield, which includes diving waves, post-critical reflections and short-spread reflections. Most successful applications of FWI are driven by the information carried by diving waves and post-critical reflections to build the long-to-intermediate wavelengths of the velocity structure. Alternative approaches, referred to as reflection waveform inversion (RWI), have been recently revisited to retrieve these long-to-intermediate wavelengths from shortspread reflections by using some prior knowledge of the reflectivity and a scale separation between the velocity macromodel and the reflectivity. This study presents a unified formalism of FWI, named as Joint FWI, whose aim is to efficiently combine the diving and reflected waves for velocity model building. The two key ingredients of Joint FWI are, on the data side, the explicit separation between the short-spread reflections and the wide-angle arrivals and, on the model side, the scale separation between the velocity macromodel and the shortscale impedance model. The velocity model and the impedance model are updated in an alternate way by Joint FWI and waveform inversion of the reflection data (least-squares migration), respectively. Starting from a crude velocity model, Joint FWI is applied to the streamer seismic data computed in the synthetic Valhall model. While the conventional FWI is stuck into a local minimum due to cycle skipping, Joint FWI succeeds in building a reliable velocity macromodel. Compared with RWI, the use of diving waves in Joint FWI improves the reconstruction of shallow velocities, which translates into an improved imaging at deeper depths. The smooth velocity model built by Joint FWI can be subsequently used as a reliable initial model for conventional FWI to increase the high-wavenumber content of the velocity model. [ABSTRACT FROM AUTHOR]

ABSTRACT Building an accurate initial velocity model for full waveform inversion (FWI) is a key issue to guarantee convergence of full waveform inversion towards the global minimum of a misfit function. In this study, we assess joint refraction and reflection stereotomography as a tool to build a reliable starting model for frequency-domain full waveform inversion from long-offset ( i.e., wide-aperture) data. Stereotomography is a slope tomographic method that is based on the inversion of traveltimes and slopes of locally-coherent events in a data cube. One advantage of stereotomography compared to conventional traveltime reflection tomography is the semi-automatic picking procedure of locally-coherent events, which is easier than the picking of continuous events, and can lead to a higher density of picks. While conventional applications of stereotomography only consider short-offset reflected waves, we assess the benefits provided by the joint inversion of reflected and refracted arrivals. Introduction of the refracted waves allows the construction of a starting model that kinematically fits the first arrivals, a necessary requirement for full waveform inversion. In a similar way to frequency-domain full waveform inversion, we design a multiscale approach of stereotomography, which proceeds hierarchically from the wide-aperture to the short-aperture components of the data, to reduce the non-linearity of the stereotomographic inversion of long-offset data. This workflow which combines stereotomography and full waveform inversion, is applied to synthetic and real data case studies for the Valhall oil-field target. The synthetic results show that the joint refraction and reflection stereotomography for a 24-km maximum offset data set provides a more reliable initial model for full waveform inversion than reflection stereotomography performed for a 4-km maximum offset data set, in particular in low-velocity gas layers and in the deep part of a structure below a reservoir. Application of joint stereotomography, full waveform inversion and reverse-time migration to real data reveals that the FWI models and the reverse-time migration images computed from the stereotomography model shares several features with FWI velocity models and migrated images computed from an anisotropic reflection-traveltime tomography model, although stereotomography was performed in the isotropic approximation. Implementation of anisotropy in joint refraction and reflection stereotomography of long-offset data is a key issue to further improve the accuracy of the method. [ABSTRACT FROM AUTHOR]

The development of dense networks of broad-band seismographs makes teleseismic data amenable to full-waveform inversion (FWI) methods for high-resolution lithospheric imaging. Compared to scattered-field migration, FWI seeks to involve the full seismic wavefield in the inversion. We present a parametric analysis of 2-D frequency-domain FWI in the framework of lithospheric imaging from teleseismic data to identify the main factors that impact on the quality of the reconstructed compressional (P)-wave and shear (S)-wave speed models. Compared to controlled-source seismology, the main adaptation of FWI to teleseismic configuration consists of the implementation with a scattered-filed formulation of plane-wave sources that impinge on the base of the lithospheric target located below the receiver network at an arbitrary incidence angle. Seismic modelling is performed with a hp-adaptive discontinuous Galerkin method on unstructured triangular mesh. A quasi-Newton inversion algorithm provides an approximate accounting for the Hessian operator, which contributes to reduce the footprint of the coarse acquisition geometry in the imaging. A versatile algorithm to compute the gradient of the misfit function with the adjoint-state method allows for abstraction between the forward-problem operators and the meshes that are during seismic modelling and inversion, respectively. An approximate correction for obliquity is derived for future application to real teleseismic data under the two-dimension approximation. Comparisons between the characteristic scales involved in exploration geophysics and in teleseismic seismology suggest that the resolution gain provided by full waveform technologies should be of the same order of magnitude for both applications. We first show the importance of the surface-reflected wavefield to dramatically improve the resolving power of FWI by combining tomography-like and migration-like imaging through the incorporation of the forward-scattered and the backscattered wavefields in the inversion. The resolution of FWI is assessed through checkerboard tests and confirms a resolution of the order of the wavelength for both the P and S speeds, when the full wavefield is incorporated in the inversion. Secondly, we show that computationally efficient strategies, which consist of decimating the number of frequency components involved in the inversion, do not apply to teleseismic acquisitions, because the scattering-angle bandwidth sampled by plane-wave sources can be narrow and coarsely sampled, compared to that provided by dense profiles of point sources in exploration seismology. The waveform inversion is less sensitive to the band of incidence angles spanned by the plane-wave sources and to the sampling of this band. However, the deficit of vertically propagating plane waves hampers the vertical resolution of planar layers. Aliasing artefacts created by coarse arrays of receivers are illustrated. We show how taking into account the Hessian in the inversion and the suitable management of frequencies in the inversion help to mitigate these artefacts. Acceptable reconstructions are shown for both the P- and S-wave speeds for a receiver spacing of up to 20 km in the 0.1–0.4 Hz frequency range. Building a reliable initial model for FWI is a highly non-linear problem in exploration seismology. We show how the low-frequency content of the teleseismic sources allow us to build accurate P- and S-wave speed models starting from simple vertical-gradient velocity models. All of these results are derived using a 2-D realistic synthetic experiment that was performed with noise-free data. Most of conclusions of this parametric study should apply in three dimensions. The impact of noise and the footprint of other experimental parameters ... [ABSTRACT FROM AUTHOR]

Double-difference tomographic methods use directly accurate time delays computed between similar signals. Such methods are designed to image very heterogeneous media, such as volcanoes or fault zones. In seismological applications, similar signals are recorded at a given station from earthquakes sharing similar and close-by sources. In seismic exploration experiments, similar signals are often recorded at neighbouring receivers. After a brief presentation of the tomographic algorithm used, a seismological application is summarized. The potential and limits of double-difference tomographic methods are explored using various numerical experiments. They show that two effects are competing in double-difference tomography: (i) the degradation of the stability of the inversion due to the geometrical proximity of the rays used in the differentiation and (ii) the decrease in modelling error, which allows improving the stability of the inversion and using smaller quantities of a priori information when data are sufficiently accurate. The best resolution is obtained for an optimal value of the inter-source or inter-receiver distance. For optimal values of these distances and a priori information, tomography using traveltime differences provides significantly better resolved results than using traveltimes. [ABSTRACT FROM AUTHOR]

Wavelet decomposition of the slowness model allows a multiscale description of the seismic first-arrival time tomography. We propose the introduction of the so-called second generation wavelets that could be used for any mesh structure and do not require a number of samples, such as the power of two in each direction for fast wavelet transform. A linearized procedure for inverting delayed travel-times considering either slowness coefficients or wavelet coefficients has been set up with an efficient ray tracing at each iteration of the inversion procedure. Wavelet decomposition over constant patches (Haar wavelet) or over linear patches (Battle-Lemarie wavelet) of coefficients at different scales are inverted as unknowns of the tomographic linearized system. Reconstruction of these coefficients depends dynamically on the local resolution when considering dense ray coverage. On simple synthetic examples, it has been found necessary to perform a local resolution analysis for specifying wavelet coefficients to be inverted. This resolution analysis could be performed for an initial smooth reconstructed medium and by designing a bit mask operator it allows fine scales to be inverted in specific areas of the model where the resolution is high while not being inverted in other areas where the resolution is poor: the wavelet decomposition will ease the multiscale reconstruction. A few synthetic examples, such as crosshole tomography or surface-surface tomography illustrate the multiscale feature of wavelet tomography. The second generation wavelet approach seems to be a flexible and rather promising tool for controlling the resolution variation of seismic first-arrival tomography. [ABSTRACT FROM AUTHOR]