Effective Q-compensated reverse time migration using new decoupled fractional Laplacian viscoacoustic wave equation

Seismic waves are attenuated and distorted during propagation because of the conversion of acoustic energy to heat energy. We focus on intrinsic attenuation, which is caused by [Formula: see text], which is the portion of energy lost during each cycle or wavelength. Amplitude attenuation can decrease the energy of the wavefields, and dispersion effects distort the phase of seismic waves. Attenuation and dispersion effects can reduce the resolution of image, and they can especially distort the real position of interfaces. On the basis of the viscoacoustic wave equation consisting of a single standard linear solid, we have derived a new viscoacoustic wave equation with decoupled amplitude attenuation and phase dispersion. Subsequently, we adopt a theoretical framework of viscoacoustic reverse time migration that can compensate the amplitude loss and the phase dispersion. Compared with the other variable fractional Laplacian viscoacoustic wave equations with decoupled amplitude attenuation and phase dispersion terms, the order of the Laplacian operator in our equation is a constant. The amplitude attenuation term is solved by pseudospectral method, and only one fast Fourier transform is required in each time step. The phase dispersion term can be computed using a finite-difference method. Numerical examples prove that our equation can accurately simulate the attenuation effects very well. Simulation of the new viscoacoustic equation indicates high efficiency because only one constant fractional Laplacian operator exists in this new viscoacoustic wave equation, which can reduce the number of inverse Fourier transforms to improve the computation efficiency of forward modeling and [Formula: see text]-compensated reverse time migration ([Formula: see text]-RTM). We tested the [Formula: see text]-RTM by using Marmousi and BP gas models and compared the [Formula: see text]-RTM images with those without compensation and attenuation (the reference image). [Formula: see text]-RTM results match well with the reference images. We also compared the field data migration images with and without compensation. Results demonstrate the accuracy and efficiency of the presented new viscoacoustic wave equation.

[1]  Yanghua Wang,et al.  A stable and efficient approach of inverse Q filtering , 2002 .

[2]  Kees Wapenaar,et al.  Wavefield extrapolation and prestack depth migration in anelastic inhomogeneous media , 2002 .

[3]  S. Brandsberg-Dahl,et al.  The 2004 BP Velocity Benchmark , 2005 .

[4]  R. Clapp,et al.  Q-model building using one-way wave-equation migration Q analysis — Part 2: 3D field-data test , 2018 .

[5]  T. Zhu,et al.  Image-Based Q Tomography using Reverse Time Q Migration , 2015 .

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

[7]  Hui Zhou,et al.  Two efficient modeling schemes for fractional Laplacian viscoacoustic wave equation , 2016 .

[8]  José M. Carcione,et al.  A generalization of the Fourier pseudospectral method , 2010 .

[9]  Gerard T. Schuster,et al.  Attenuation compensation for least-squares reverse time migration using the viscoacoustic-wave equation , 2014 .

[10]  Carlos F.M. Coimbra,et al.  Mechanics with variable‐order differential operators , 2003 .

[11]  Biondo Biondi,et al.  Q-compensated reverse-time migration , 2014 .

[12]  G. McMechan,et al.  Viscoelastic true-amplitude prestack reverse-time depth migration , 2008 .

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

[14]  José M. Carcione,et al.  Wave Fields in Real Media: Wave Propagation in Anisotropic, Anelastic and Porous Media , 2011 .

[15]  Tieyuan Zhu,et al.  Modeling acoustic wave propagation in heterogeneous attenuating media using decoupled fractional Laplacians , 2014 .

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

[17]  Ketil Hokstad,et al.  Prestack depth migration with compensation for absorption and dispersion , 1995 .

[18]  Robert G. Clapp,et al.  Q-model building using one-way wave-equation migration Q analysis — Part 1: Theory and synthetic test , 2018 .

[19]  Viscoacoustic modeling and imaging using low-rank approximation , 2014 .

[20]  Lowrank approximation for time domain viscoacoustic wave equation with spatially varying order fractional Laplacians , 2014 .

[21]  G. McMechan,et al.  True-amplitude prestack depth migration , 2007 .

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

[23]  José M. Carcione,et al.  Time-domain Modeling of Constant-Q Seismic Waves Using Fractional Derivatives , 2002 .

[24]  Wei Huang,et al.  Absorption decomposition and compensation via a two-step scheme , 2015 .

[25]  Robert G. Clapp,et al.  Wave-equation Migration Q Analysis (WEMQA) , 2013 .

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

[27]  Guo Yong.  High Resolution Seismic Data Processing , 2002 .

[28]  Hui Zhou,et al.  Efficient reverse time migration based on fractional Laplacian viscoacoustic wave equation , 2016 .

[29]  Vicente Grau,et al.  Fractional diffusion models of cardiac electrical propagation: role of structural heterogeneity in dispersion of repolarization , 2014, Journal of The Royal Society Interface.

[30]  Causse Emmanuel,et al.  Viscoacoustic reverse-time migration , 2000 .

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