Attenuation compensation for least-squares reverse time migration using the viscoacoustic-wave equation

Strong subsurface attenuation leads to distortion of amplitudes and phases of seismic waves propagating inside the earth. Conventional acoustic reverse time migration (RTM) and least-squares reverse time migration (LSRTM) do not account for this distortion, which can lead to defocusing of migration images in highly attenuative geologic environments. To correct for this distortion, we used a linearized inversion method, denoted as Qp-LSRTM. During the leastsquares iterations, we used a linearized viscoacoustic modeling operator for forward modeling. The adjoint equations were derived using the adjoint-state method for back propagating the residual wavefields. The merit of this approach compared with conventional RTM and LSRTM was that Qp-LSRTM compensated for the amplitude loss due to attenuation and could produce images with better balanced amplitudes and more resolution below highly attenuative layers. Numerical tests on synthetic and field data illustrated the advantages of Qp-LSRTM over RTM and LSRTM when the recorded data had strong attenuation effects. Similar to standard LSRTM, the sensitivity tests for background velocity and Qp errors revealed that the liability of this method is the requirement for smooth and accurate migration velocity and attenuation models.

[1]  Einar Kjartansson,et al.  Constant Q-wave propagation and attenuation , 1979 .

[2]  G. Schuster,et al.  Least-squares migration of incomplete reflection data , 1999 .

[3]  S. Bickel,et al.  Plane-wave Q deconvolution , 1985 .

[4]  S. T. Chen,et al.  Subsurface imaging using reversed vertical seismic profiling and crosshole tomographic methods , 1990 .

[5]  Youli Quan,et al.  Seismic attenuation tomography using the frequency shift method , 1997 .

[6]  G. Schuster,et al.  Acoustic wave-equation traveltime and waveform inversion of crosshole seismic data , 1995 .

[7]  Wei Dai,et al.  Least-squares Migration of Simultaneous Sources Data With a Deblurring Filter , 2009 .

[8]  A. Tarantola.  LINEARIZED INVERSION OF SEISMIC REFLECTION DATA , 1984 .

[9]  Efficient iterative viscoacoustic linearized inversion , 1995 .

[10]  Shuki Ronen,et al.  Least-squares Reverse Time Migration/inversion For Ocean Bottom Data: a Case Study , 2011 .

[11]  Bertrand Duquet,et al.  Kirchhoff modeling, inversion for reflectivity, and subsurface illumination , 2000 .

[12]  J. Leveille,et al.  Attenuation Compensation in Viscoacoustic Reserve-Time Migration , 2013 .

[13]  J. Leveille,et al.  Waveform Inversion With Attenuation , 2012 .

[14]  Joakim O. Blanch,et al.  Modeling of a constant Q; methodology and algorithm for an efficient and optimally inexpensive viscoelastic technique , 1995 .

[15]  Yanghua Wang.  Inverse-Q filtered migration , 2008 .

[16]  R. Christensen.  Theory of viscoelasticity : an introduction , 1971 .

[17]  Wei Dai,et al.  Multi‐source least‐squares reverse time migration , 2012 .

[18]  G. McMechan.  MIGRATION BY EXTRAPOLATION OF TIME‐DEPENDENT BOUNDARY VALUES* , 1983 .

[19]  E. Baysal,et al.  Reverse time migration , 1983 .

[20]  N. Hargreaves,et al.  Inverse Q filtering by Fourier transform , 1991 .

[21]  Matthew A. Brzostowski,et al.  3-D tomographic imaging of near‐surface seismic velocity and attenuation , 1992 .

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

[23]  Wei Dai,et al.  Multisource Least-squares Reverse Time Migration , 2012 .

[24]  3D prestack depth migration with compensation for frequency dependent absorption and dispersion , 2009 .

[25]  Yaxun Tang.  Target-oriented wave-equation least-squares migration/inversion with phase-encoded Hessian , 2009 .

[26]  Gerard T. Schuster,et al.  Multisource least‐squares migration of marine streamer and land data with frequency‐division encoding , 2012 .

[27]  N. Chemingui,et al.  Wave Equation Migration With Attenuation And Anisotropy Compensation , 2011 .

[28]  G. West,et al.  Inverse Q Migration , 1994 .

[29]  Joakim O. Blanch,et al.  Viscoelastic finite-difference modeling , 1994 .

[30]  Matthew A. Brzostowski,et al.  3-D Tomographic Imaging of Near-surface Seismic Velocity And Attenuation , 1991 .

[31]  Xin Wang,et al.  Least-squares migration of multisource data with a deblurring filter , 2011 .

[32]  R. Plessix,et al.  Frequency-domain finite-difference amplitude-preserving migration , 2004 .

[33]  Gerard T. Schuster,et al.  Multi-source Wave Equation Least-squares Migration with a Deblurring Filter , 2010 .

[34]  N. Whitmore.  Iterative Depth Migration By Backward Time Propagation , 1983 .

[35]  F. Santosa,et al.  Inverse Problems of Acoustic and Elastic Waves , 1984 .

[36]  Dave Nichols,et al.  Wavepath-consistent Effective Q Estimation for Q-compensated Reverse-time Migration , 2012 .

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

[38]  Robert H. Stolt,et al.  Seismic migration : theory and practice , 1986 .

[39]  M. Carcione,et al.  Wave propagation simulation in a linear viscoacoustic medium , 1997 .

[40]  Joakim O. Blanch,et al.  Viscoelastic finite-difference modeling , 1994 .

[41]  Joe Zhou,et al.  Compensating for Visco-Acoustic Effects in TTI Reverse Time Migration , 2010 .

[42]  Bin Wang,et al.  Compensating Visco-Acoustic Effects in Anisotropic Resverse-Time Migration , 2012 .

[43]  Gerard T. Schuster,et al.  3D Multi-source Least-squares Reverse Time Migration , 2010 .

[44]  Yu Zhang,et al.  Compensating for visco-acoustic effects with reverse time migration , 2010 .