On the footprint of anisotropy on isotropic full waveform inversion: the Valhall case study

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 inconsis- tency 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 illus- trated 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 hierar- chical 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 δ andwere previously developed by anisotropic reflection traveltime tomography. The normal moveout (VNMO = V0 √ 1 + 2δ) and horizontal (Vh = V0 √ 1 + 2� ) 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 tomog- raphy (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 reconstruc- tions 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.

[1]  Dong-Joo Min,et al.  Frequency-domain elastic full waveform inversion for VTI media , 2010 .

[2]  P. Jaiswal,et al.  Seismic imaging of the Naga thrust using multiscale waveform inversion , 2009 .

[3]  L. Sirgue,et al.  The Importance of Low Frequency and Large Offset in Waveform Inversion , 2006 .

[4]  Sensitivity Study of Seismic Reflection/Refraction Data , 1995 .

[5]  Romain Brossier,et al.  Two-dimensional Seismic Imaging of the Valhall Model From Synthetic OBC Data By Frequency Domain Elastic Full-waveform Inversion , 2009 .

[6]  Romain Brossier,et al.  Two-dimensional frequency-domain visco-elastic full waveform inversion: Parallel algorithms, optimization and performance , 2011, Comput. Geosci..

[7]  G. H. F. Gardner,et al.  FORMATION VELOCITY AND DENSITY—THE DIAGNOSTIC BASICS FOR STRATIGRAPHIC TRAPS , 1974 .

[8]  R. Shipp,et al.  Seismic waveform inversion in the frequency domain, Part 2: Fault delineation in sediments using crosshole data , 1999 .

[9]  Albert Tarantola,et al.  Inverse problem theory - and methods for model parameter estimation , 2004 .

[10]  E. Polak,et al.  Note sur la convergence de méthodes de directions conjuguées , 1969 .

[11]  P. Mora Elastic wave‐field inversion of reflection and transmission data , 1988 .

[12]  Marwan Charara,et al.  The domain of applicability of acoustic full-waveform inversion for marine seismic data , 2009 .

[13]  William W. Symes,et al.  Migration velocity analysis and waveform inversion , 2008 .

[14]  R. Plessix Three-dimensional frequency-domain full-waveform inversion with an iterative solver , 2009 .

[15]  R. Plessix A review of the adjoint-state method for computing the gradient of a functional with geophysical applications , 2006 .

[16]  A. Tarantola,et al.  Two‐dimensional nonlinear inversion of seismic waveforms: Numerical results , 1986 .

[17]  Jean Virieux,et al.  Building starting model for full waveform inversion from wide-aperture data by stereotomography , 2010 .

[18]  Jean Virieux,et al.  Crustal seismic imaging from multifold ocean bottom seismometer data by frequency domain full waveform tomography: Application to the eastern Nankai trough , 2006 .

[19]  Jean-Yves L'Excellent,et al.  FWT2D: A massively parallel program for frequency-domain full-waveform tomography of wide-aperture seismic data - Part 1: Algorithm , 2009, Comput. Geosci..

[20]  Jean Virieux,et al.  An overview of full-waveform inversion in exploration geophysics , 2009 .

[21]  Alan G. Green,et al.  Effects of near‐surface waveguides on shallow high‐resolution seismic refraction and reflection data , 1996 .

[22]  Olav I. Barkved,et al.  Thematic Set: Full waveform inversion: the next leap forward in imaging at Valhall , 2010 .

[23]  Patrick R. Amestoy,et al.  MUltifrontal Massively Parallel Solver (MUMPS Version 4.3) Users' guide , 2003 .

[24]  Jean-Yves L'Excellent,et al.  FWT2D: A massively parallel program for frequency-domain full-waveform tomography of wide-aperture seismic data - Part 2: Numerical examples and scalability analysis , 2009, Computational Geosciences.

[25]  Hicks,et al.  Gauss–Newton and full Newton methods in frequency–space seismic waveform inversion , 1998 .

[26]  Antoine Guitton,et al.  A Preconditioning Scheme For Full Waveform Inversion , 2010 .

[27]  R. Shipp,et al.  Two-dimensional full wavefield inversion of wide-aperture marine seismic streamer data , 2002 .

[28]  R. Pratt,et al.  Reconciliation of crosshole seismic velocities with well information in a layered sedimentary environment , 1996 .

[29]  Ying Rao,et al.  Reflection seismic waveform tomography , 2009 .

[30]  R. Pratt,et al.  INVERSE THEORY APPLIED TO MULTI‐SOURCE CROSS‐HOLE TOMOGRAPHY.: PART 1: ACOUSTIC WAVE‐EQUATION METHOD1 , 1990 .

[31]  Florian Bleibinhaus,et al.  Structure of the California Coast Ranges and San Andreas Fault at SAFOD from seismic waveform inversion and reflection imaging , 2007 .

[32]  R. Gerhard Pratt,et al.  Efficient waveform inversion and imaging: A strategy for selecting temporal frequencies , 2004 .

[33]  R. Plessix,et al.  Thematic Set: Full waveform inversion of a deep water ocean bottom seismometer dataset , 2010 .

[34]  C. Bunks,et al.  Multiscale seismic waveform inversion , 1995 .

[35]  René-Édouard Plessix,et al.  Application of Acoustic Full Waveform Inversion to a Low-frequency Large-offset Land Data Set , 2010 .

[36]  Tariq Alkhalifah,et al.  Velocity analysis for transversely isotropic media , 1995 .

[37]  P. Mora Nonlinear two-dimensional elastic inversion of multioffset seismic data , 1987 .

[38]  R. Pratt Seismic waveform inversion in the frequency domain; Part 1, Theory and verification in a physical scale model , 1999 .

[39]  Jean-Luc Mari,et al.  Signal Processing for Geologists and Geophysicists , 1999 .

[40]  René-Édouard Plessix,et al.  A parametrization study for surface seismic full waveform inversion in an acoustic vertical transversely isotropic medium , 2011 .

[41]  Olav I. Barkved,et al.  Initial Experience Operating a Permanent 4C Seabed Array For Reservoir Monitoring At Valhall. , 2004 .

[42]  J. Munns,et al.  The Valhall Field: a geological overview , 1985 .

[43]  R. G. Pratt,et al.  Two‐dimensional velocity models from wide‐angle seismic data by wavefield inversion , 1996 .

[44]  Ilya Tsvankin,et al.  Normal moveout from dipping reflectors in anisotropic media , 1995 .

[45]  Jean Virieux,et al.  Multiscale imaging of complex structures from multifold wide-aperture seismic data by frequency-domain full-waveform tomography: application to a thrust belt , 2004 .

[46]  Denes Vigh,et al.  3D Full waveform inversion on a Gulf of Mexico WAZ data set , 2010 .

[47]  S. Operto,et al.  Which data residual norm for robust elastic frequency-domain full waveform inversion? , 2010 .

[48]  Guy Chavent,et al.  Nonlinear Least Squares for Inverse Problems , 2010 .

[49]  Jean Virieux,et al.  Two-dimensional Acoustic Anisotropic (VTI) Full Waveform Inversion: the Valhall Case Study , 2011 .

[50]  J. Scales,et al.  Bayesian seismic waveform inversion: Parameter estimation and uncertainty analysis , 1998 .

[51]  L. Thomsen Weak elastic anisotropy , 1986 .

[52]  P. Podvin,et al.  Finite difference computation of traveltimes in very contrasted velocity models: a massively parallel approach and its associated tools , 1991 .

[53]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[54]  신창수 Apparatus, methods and recording medium for imaging a subsurface using the waveform inversion in the laplace-fourier domain , 2008 .

[55]  Jean Virieux,et al.  Finite-difference frequency-domain modeling of viscoacoustic wave propagation in 2D tilted transversely isotropic (TTI) media , 2009 .

[56]  Romain Brossier,et al.  Parsimonious finite-volume frequency-domain method for 2-D P–SV-wave modelling , 2008 .

[57]  J. Sheng,et al.  Early arrival waveform tomography on near-surface refraction data , 2006 .

[58]  S. Singh,et al.  Sensitivity Study of Seismic Reflection/Refraction Data , 1995 .

[59]  R. G. Pratt,et al.  Full waveform tomography for lithospheric imaging: results from a blind test in a realistic crustal model , 2007 .

[60]  Ilya Tsvankin,et al.  Seismic Signatures and Analysis of Reflection Data in Anisotropic Media , 2001 .

[61]  R. G. Pratt,et al.  P092 SEISMIC WAVEFORM TOMOGRAPHY: THE EFFECT OF LAYERING AND ANISOTROPY , 2001 .

[62]  Jan H. Kommedal,et al.  3D Waveform Inversion on Valhall Wide-azimuth OBC , 2009 .

[63]  R. Pratt Inverse theory applied to multisource cross-hole tomography, Part2 : Elastic wave-equation method , 1990 .

[64]  S. Operto,et al.  Seismic imaging of complex onshore structures by 2D elastic frequency-domain full-waveform inversion , 2009 .

[65]  Jean Virieux,et al.  Acoustic VTI Full Waveform Inversion: Sensitivity Analysis And Realistic Synthetic Examples , 2011 .