Search

Multiparameter full waveform inversion (FWI) is a challenging quantitative seismic imaging method for lithological characterization and reservoir monitoring. The difficulties in multiparameter FWI arise from the variable influence of the different parameter classes on the phase and amplitude of the data, and the trade-off between these. In this framework, choosing a suitable parametrization of the subsurface and designing the suitable FWI workflow are two key methodological issues in non-linear waveform inversion. We assess frequency-domain visco-acoustic FWI to reconstruct the compressive velocity (VP), the density (ρ) or the impedance (IP) and the quality factor (QP), from the hydrophone component, using a synthetic data set that is representative of the Valhall oil field in the North Sea. We first assess which of the (VP, ρ) and (VP, IP) parametrizations provides the most reliable FWI results when dealing with wide-aperture data. Contrary to widely accepted ideas, we show that the (VP, ρ) parametrization allows a better reconstruction of both the VP, ρ and IP parameters, first because it favours the broad-band reconstruction of the dominant VP parameter, and secondly because the trade-off effects between velocity and density at short-to-intermediate scattering angles can be removed by multiplication, to build an impedance model. This allows for the matching of the reflection amplitudes, while the broad-band velocity model accurately describes the kinematic attributes of both the diving waves and reflections. Then, we assess different inversion strategies to recover the quality factor QP, in addition to parameters VP and ρ. A difficulty related to attenuation estimation arises because, on the one hand the values of QP are on average one order of magnitude smaller than those of VP and ρ, and on the other hands model perturbations relative to the starting models can be much higher for QP than for VP and ρ during FWI. In this framework, we show that an empirical tuning of the FWI regularization, which is adapted to each parameter class, is a key issue to correctly account for the attenuation in the inversion. We promote a hierarchical approach where the dominant parameter VP is reconstructed first from the full data set (i.e. without any data preconditioning) to build a velocity model as kinematically accurate as possible before performing the joint update of the three parameter classes during a second step. This hierarchical imaging of compressive wave speed, density and attenuation is applied to a real wide-aperture ocean-bottom-cable data set from the Valhall oil field. Several geological features, such as accumulation of gas below barriers of claystone and soft quaternary sediment are interpreted in the FWI models of density and attenuation. The models of VP, ρ and QP that have been developed by visco-acoustic FWI of the hydrophone data can be used as initial models to perform visco-elastic FWI of the geophone data for the joint update of the compressive and shear wave speeds. [ABSTRACT FROM AUTHOR]

Multiparameter elastic full waveform inversion (FWI) is a promising technology that allows inferences to be made on rock and fluid properties, which thus narrows the gap between seismic imaging and reservoir characterization. Here, we assess the feasibility of 2-D vertical transverse isotropic visco-elastic FWI of wide-aperture multicomponent ocean-bottom-cable data from the Valhall oil field. A key issue is to design a suitable hierarchical data-driven and model-driven FWI workflow, the aim of which is to reduce the nonlinearity of the FWI. This nonlinearity partly arises because the shear (S) wavespeed can have a limited influence on seismic data in marine environments. In a preliminary stage, visco-acoustic FWI of the hydrophone component is performed to build a compressional (P)-wave velocity model, a density model and a quality-factor model, which provide the necessary background models for the subsequent elastic inversion. During the elastic FWI, the P and S wavespeeds are jointly updated in two steps. First, the hydrophone data are inverted to mainly update the long-to-intermediate wavelengths of the S wavespeeds from the amplitude-versus-offset variations of the P–P reflections. This S-wave velocity model is used as an improved starting model for the subsequent inversion of the better-resolving data recorded by the geophones. During these two steps, the P-wave velocity model is marginally updated, which supports the relevance of the visco-acoustic FWI results. Through seismic modelling, we show that the resulting visco-elastic model allows several P-to-S converted phases recorded on the horizontal-geophone component to be matched. Several elastic quantities, such as the Poisson ratio, and the ratio and product between the P and S wavespeeds, are inferred from the P-wave and S-wave velocity models. These attributes provide hints for the interpretation of an accumulation of gas below lithological barriers. [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]

SUMMARY The validity of isotropic approximation to perform acoustic full waveform inversion (FWI) of real wide-aperture anisotropic data can be questioned due to the intrinsic kinematic inconsistency between short- and large-aperture components of the data. This inconsistency is mainly related to the differences between the vertical and horizontal velocities in vertical-transverse isotropic (VTI) media. The footprint of VTI anisotropy on 2-D acoustic isotropic FWI is illustrated on a hydrophone data set of an ocean-bottom cable that was collected over the Valhall field in the North Sea. Multiscale FWI is implemented in the frequency domain by hierarchical inversions of increasing frequencies and decreasing aperture angles. The FWI models are appraised by local comparison with well information, seismic modelling, reverse-time migration (RTM) and source-wavelet estimation. A smooth initial VTI model parameterized by the vertical velocity V0 and the Thomsen parameters δ and ε were previously developed by anisotropic reflection traveltime tomography. The normal moveout () and horizontal () velocity models were inferred from the anisotropic models to perform isotropic FWI. The VNMO models allows for an accurate match of short-spread reflection traveltimes, whereas the Vh model, after updating by first-arrival traveltime tomography (FATT), allows for an accurate match of first-arrival traveltimes. Ray tracing in the velocity models shows that the first 1.5 km of the medium are sampled by both diving waves and reflections, whereas the deeper structure at the reservoir level is mainly controlled by short-spread reflections. Starting from the initial anisotropic model and keeping fixed δ and ε models, anisotropic FWI allows us to build a vertical velocity model that matches reasonably well the well-log velocities. Isotropic FWI is performed using either the NMO model or the FATT model as initial model. In both cases, horizontal velocities are mainly reconstructed in the first 1.5 km of the medium. This suggests that the wide-aperture components of the data have a dominant control on the velocity estimation at these depths. These high velocities in the upper structure lead to low values of velocity in the underlying gas layers (either equal or lower than vertical velocities of the well log), and/or a vertical stretching of the structure at the reservoir level below the gas. This bias in the gas velocities and the mispositioning in depth of the deep reflectors, also shown in the RTM images, are required to match the deep reflections in the isotropic approximation and highlight the footprint of anisotropy in the isotropic FWI of long-offset data. Despite the significant differences between the anisotropic and isotropic FWI models, each of these models produce a nearly-equivalent match of the data, which highlights the ill-posedness of acoustic anisotropic FWI. Hence, we conclude with the importance of considering anisotropy in FWI of wide-aperture data to avoid bias in the velocity reconstructions and mispositioning in depth of reflectors. Designing a suitable parameterization of the VTI acoustic FWI is a central issue to manage the ill-posedness of the FWI. [ABSTRACT FROM AUTHOR]