High-resolution seismic attenuation imaging from wide-aperture onshore data by visco-acoustic frequency-domain full-waveform inversion

SUMMARY Here we assess the potential of the visco-acoustic frequency domain full-waveform inversion (FWI) to reconstruct P-wave velocity (VP) and P-wave attenuation factor (Q) from surface onshore seismic data. First, we perform a sensitivity analysis of the FWI based upon a grid search analysis of the misfit function and several synthetic FWI examples using velocity and Q models of increasing complexity. Subsequently, we applied both the acoustic and viscoacousticFWItorealsurfacewide-apertureonshoreseismicdatafromthePolishBasin,wherea strongattenuationoftheseismicdataisobserved.Thesensitivityanalysisofthevisco-acoustic FWI suggests that the FWI can jointly reconstruct the velocity and the attenuation factor if the signatureoftheattenuationissufficientlystronginthedata.Asyntheticexamplecorresponding to a homogeneous background model with an inclusion shows a reliable reconstruction of VP and Q in the inclusion, when Q is as small as 90 and 50 in the background model and in the inclusion, respectively. A first application of acoustic FWI to real data shows that a heuristic normalization of the data with offset allows us to compensate for the effect of the attenuation in the data and reconstruct a reliable velocity model. Alternatively, we show that visco-acoustic FWI allows us to reconstruct jointly both a reliable velocity model and a Q model from the true-amplitude data. We propose a pragmatical approach based upon seismic modelling and source wavelet estimation to infer the best starting homogeneous Q model for visco-acoustic FWI. We find the source wavelet estimation quite sensitive to the quality of the velocity and attenuation models used for the estimation. For example, source-to-source wavelets are significantly more consistent when computed in the final FWI model than in the initial one. A good kinematic and amplitude match between the early-arriving phases of the real and time-domain synthetic seismograms computed in the final FWI model provides an additional evidence of the reliability of the final FWI model. We find the recovered velocity and attenuation models consistent with the expected lithology and stratigraphy in the study area. We link high-attenuation zones with the increased clay content and the presence of the mineralized fluids.

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

[2]  Angus I. Best,et al.  The relationships between the velocities, attenuations and petrophysical properties of reservoir sedimentary rocks1 , 1994 .

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

[4]  Stephen J. Wright,et al.  Numerical Optimization , 2018, Fundamental Statistical Inference.

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

[6]  J. Dvorkin,et al.  Modeling attenuation in reservoir and nonreservoir rock , 2006 .

[7]  Stéphane Operto,et al.  Quantitative imaging of the Permo-Mesozoic complex and its basement by frequency domain waveform tomography of wide-aperture seismic data from the Polish Basin , 2008 .

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

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

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

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

[12]  A. Tarantola.  Inversion of seismic reflection data in the acoustic approximation , 1984 .

[13]  Theodoros Klimentos,et al.  Relationships among compressional wave attenuation, porosity, clay content, and permeability in sandstones , 1990 .

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

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

[16]  W. Mulder,et al.  Migration for velocity and attenuation perturbations , 2010 .

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

[18]  S. Operto,et al.  Mixed‐grid and staggered‐grid finite‐difference methods for frequency‐domain acoustic wave modelling , 2004 .

[19]  A. Tarantola.  Inverse problem theory : methods for data fitting and model parameter estimation , 1987 .

[20]  Klaus Bauer,et al.  Waveform tomography images of velocity and inelastic attenuation from the Mallik 2002 crosshole seismic surveys , 2005 .

[21]  Theodoros Klimentos,et al.  Attenuation of P- and S-waves as a method of distinguishing gas and condensate from oil and water , 1995 .

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

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

[24]  Yanghua Wang,et al.  Seismic Inverse Q Filtering , 2008 .

[25]  C. Mccann,et al.  Sonic to ultrasonic Q of sandstones and limestones: Laboratory measurements at in situ pressures , 2009 .

[26]  Graham. J. Hicks,et al.  Arbitrary source and receiver positioning in finite‐difference schemes using Kaiser windowed sinc functions , 2002 .

[27]  Romain Brossier,et al.  Robust elastic frequency‐domain full‐waveform inversion using the L1 norm , 2009 .

[28]  Z. M. Song,et al.  Frequency-domain acoustic-wave modeling and inversion of crosshole data; Part II, Inversion method, synthetic experiments and real-data results , 1995 .

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

[30]  Wenzhi Zhao,et al.  Nonlinear process control of wave-equation inversion and its application in the detection of gas , 2007 .

[31]  S. Seneviratne,et al.  A regional perspective on trends in continental evaporation , 2009 .

[32]  Detecting near-surface objects with seismic waveform tomography , 2009 .

[33]  Fuchun Gao,et al.  Waveform tomography at a groundwater contamination site: VSP-surface data set , 2006 .

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

[35]  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 .

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

[37]  G. McMechan,et al.  Multifrequency viscoacoustic modeling and inversion , 1996 .

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

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

[40]  S. Assefa,et al.  Attenuation of P‐ and S‐waves in limestones , 1999 .

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

[42]  Effective sub-Zechstein salt imaging using low-frequency seismics — Results of the GRUNDY 2003 experiment across the Variscan front in the Polish Basin , 2007 .

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

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

[45]  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.

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

[47]  P. Williamson,et al.  Frequency-domain acoustic-wave modeling and inversion of crosshole data; Part 1, 2.5-D modeling method , 1995 .

[48]  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 .

[49]  Stéphane Rondenay,et al.  Effects of surface scattering in full-waveform inversion , 2009 .

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

[51]  George A. McMechan,et al.  2.5D full‐wavefield viscoacoustic invercion1 , 1995 .

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