On the role of density and attenuation in three-dimensional multiparameter viscoacoustic VTI frequency-domain FWI: an OBC case study from the North Sea

3-D frequency-domain full waveform inversion (FWI) is applied on North Sea wide-azimuth ocean-bottom cable data at low frequencies (≤10 Hz) to jointly update vertical wave speed, density and quality factor Q in the viscoacoustic VTI approximation. We assess whether density and Q should be viewed as proxy to absorb artefacts resulting from approximate wave physics or are valuable for interpretation in the presence of soft sediments and gas cloud. FWI is performed in the frequency domain to account for attenuation easily. Multiparameter frequency-domain FWI is efficiently performed with a few discrete frequencies following a multiscale frequency continuation. However, grouping a few frequencies during each multiscale step is necessary to mitigate acquisition footprint and match dispersive shallow guided waves. Q and density absorb a significant part of the acquisition footprint hence cleaning the velocity model from this pollution. Low Q perturbations correlate with low-velocity zones associated with soft sediments and gas cloud. However, the amplitudes of the Q perturbations show significant variations when the inversion tuning is modified. This dispersion in the Q reconstructions is however not passed on the velocity parameter suggesting that cross-talks between first-order kinematic and second-order dynamic parameters are limited. The density model shows a good match with a well log at shallow depths. Moreover, the impedance built a posteriori from the FWI velocity and density models shows a well-focused image with however local differences with the velocity model near the sea bed where density might have absorbed elastic effects. The FWI models are finally assessed against time-domain synthetic seismogram modelling performed with the same frequency-domain modelling engine used for FWI.

[1]  A. Calvert,et al.  Seismic waveform tomography across the Seattle Fault Zone in Puget Sound: resolution analysis and effectiveness of visco-acoustic inversion of viscoelastic data , 2013 .

[2]  Joint inversion of velocity and density in preserved-amplitude full-waveform inversion , 2016 .

[3]  A. Green,et al.  Guided waves in near‐surface seismic surveys , 1998 .

[4]  Ludovic Métivier,et al.  Efficient 3-D frequency-domain mono-parameter full-waveform inversion of ocean-bottom cable data: application to Valhall in the visco-acoustic vertical transverse isotropic approximation , 2015 .

[5]  M. Warner,et al.  Anisotropic 3D full-waveform inversion , 2013 .

[6]  A. Fichtner,et al.  Block-diagonal Approximate Hessian for Preconditioning in Full Waveform Inversion , 2013 .

[7]  K. Marfurt Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations , 1984 .

[8]  Ludovic Métivier,et al.  A guided tour of multiparameter full-waveform inversion with multicomponent data: From theory to practice , 2013 .

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

[10]  Estimating seismic velocities below the sea-bed using surface waves , 2004 .

[11]  Yuzhu Liu,et al.  Simultaneous estimation of velocity and density in acoustic multiparameter full-waveform inversion using an improved scattering-integral approach , 2016 .

[12]  Wim A. Mulder,et al.  An ambiguity in attenuation scattering imaging , 2009 .

[13]  R. Pratt,et al.  Reflection waveform inversion using local descent methods: Estimating attenuation and velocity over a gas-sand deposit , 2001 .

[14]  S. Operto,et al.  Fast full waveform inversion with source encoding and second-order optimization methods , 2015 .

[15]  V. Akcelik,et al.  Full Waveform Inversion for Seismic Velocity and Anelastic Losses in Heterogeneous Structures , 2007 .

[16]  Ludovic Métivier,et al.  Computationally efficient three-dimensional acoustic finite-difference frequency-domain seismic modeling in vertical transversely isotropic media with sparse direct solver , 2014 .

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

[18]  Ludovic Métivier,et al.  The SEISCOPE optimization toolbox: A large-scale nonlinear optimization library based on reverse communication , 2016 .

[19]  Scott A. Morton,et al.  Time-lapse full-waveform inversion with ocean-bottom-cable data: Application on Valhall field , 2016 .

[20]  Walter I. Futterman,et al.  Dispersive body waves , 1962 .

[21]  H. Kolsky,et al.  LXXI. The propagation of stress pulses in viscoelastic solids , 1956 .

[22]  M. Toksöz,et al.  Diffraction tomography and multisource holography applied to seismic imaging , 1987 .

[23]  Bastien Dupuy,et al.  A downscaling strategy from FWI to microscale reservoir properties from high-resolution images , 2016 .

[24]  Gregory Beylkin,et al.  A new slant on seismic imaging: Migration and integral geometry , 1987 .

[25]  D. Köhn,et al.  The Role of Density in Acoustic Full Waveform Inversion of Marine Reflection Seismics , 2012 .

[26]  D. Min,et al.  Full waveform inversion strategy for density in the frequency domain , 2012 .

[27]  Jean Virieux,et al.  Asymptotic viscoacoustic diffraction tomography of ultrasonic laboratory data: a tool for rock properties analysis , 2000 .

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

[29]  S. Operto,et al.  Which parameterization is suitable for acoustic vertical transverse isotropic full waveform inversion? Part 1: Sensitivity and trade-off analysis , 2013 .

[30]  Patrick R. Amestoy,et al.  Fast 3D frequency-domain full-waveform inversion with a parallel block low-rank multifrontal direct solver: Application to OBC data from the North Sea , 2016 .

[31]  D. Boiero,et al.  Surface- and guided-wave inversion for near-surface modeling in land and shallow marine seismic data , 2013 .

[32]  R. G. Pratt,et al.  Inversion strategies for visco-acoustic waveform inversion , 2013 .

[33]  Denes Vigh,et al.  Elastic full-waveform inversion application using multicomponent measurements of seismic data collection , 2014 .

[34]  Ludovic Métivier,et al.  Acoustic Multi-Parameter FWI for the Reconstruction of P-Wave Velocity, Density and Attenuation: Preconditioned Truncated Newton Approach , 2015 .

[35]  Stéphane Operto,et al.  3D ray+Born migration/inversion—Part 1: Theory , 2003 .

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

[37]  Paul G. Richards,et al.  Quantitative Seismology: Theory and Methods , 1980 .

[38]  S. Operto,et al.  3D finite-difference frequency-domain modeling of visco-acoustic wave propagation using a massively parallel direct solver: A feasibility study , 2007 .

[39]  Romain Brossier,et al.  Multiparameter full waveform inversion of multicomponent ocean-bottom-cable data from the Valhall field. Part 1: imaging compressional wave speed, density and attenuation , 2013 .

[40]  D. Köhn,et al.  Acoustic full waveform tomography in the presence of attenuation: a sensitivity analysis , 2013 .

[41]  Andrew Ratcliffe,et al.  Anisotropic 3 D full-waveform inversion , 2013 .

[42]  Wim A. Mulder,et al.  Seismic attenuation imaging with causality , 2011 .

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

[44]  C. Shin,et al.  Improved amplitude preservation for prestack depth migration by inverse scattering theory , 2001 .

[45]  Stéphane Operto,et al.  3 D ray + Born migration / inversion — Part 1 : Theory , 2003 .

[46]  L. Neil Frazer,et al.  Practical aspects of reflectivity modeling , 1987 .

[47]  Stéphane Operto,et al.  High-resolution seismic attenuation imaging from wide-aperture onshore data by visco-acoustic frequency-domain full-waveform inversion , 2011 .

[48]  Pengliang Yang,et al.  A review on the systematic formulation of 3-D multiparameter full waveform inversion in viscoelastic medium , 2016 .