Source-independent full wavefield converted-phase elastic migration velocity analysis

Converted phase (CP) elastic seismic signals are comparable in amplitude to the primary signals recorded at large offsets and have the potential to be used in seismic imaging and velocity analysis. We present an approach for CP elastic wave equation velocity analysis that does not use source information and is applicable to surface-seismic, microseismic, teleseismic and vertical seismic profile (VSP) studies. Our approach is based on the cross-correlation between reflected or transmitted PP and CP PS (and/or SS and CP SP) waves propagated backward in time, and is formulated as an optimization problem with a differential semblance criterion objective function for the simultaneous update of both P- and S-wave velocity models. The merit of this approach is that it is fully data-driven, uses full waveform information, and requires only one elastic backward propagation to form an image rather than the two (one forward and one backward) propagations needed for standard reverse-time migration. Moreover, as the method does not require forward propagation, it does not suffer from migration operator source aliasing when a small number of shots are used. We present a derivation of the method and test it with a synthetic model and field micro-seismic data.

[1]  V. Freudenthaler,et al.  Long-range transport of Saharan dust to northern Europe : The 11-16 October 2001 outbreak observed with EARLINET , 2003 .

[2]  William W. Symes,et al.  Wave equation migration velocity analysis by differential semblance optimization , 2005 .

[3]  Bob A. Hardage,et al.  Multicomponent Seismic Technology , 2011 .

[4]  Converted phase elastic migration velocity analysis , 2013 .

[5]  Guy Purnell.  Imaging beneath a high‐velocity layer using converted waves , 1992 .

[6]  P. Williamson,et al.  On the gradient artifacts in migration velocity analysis based on differential semblance optimization , 2010 .

[7]  Kim B. Olsen,et al.  On the implementation of perfectly matched layers in a three‐dimensional fourth‐order velocity‐stress finite difference scheme , 2003 .

[8]  G. McMechan,et al.  2D and 3D elastic wavefield vector decomposition in the wavenumber domain for VTI media , 2010 .

[9]  Paul Sava,et al.  Adjoint wave-equation velocity analysis , 2006 .

[10]  J. Carcione.  Staggered mesh for the anisotropic and viscoelastic wave equation , 1999 .

[11]  Jia Yan,et al.  Improving the efficiency of elastic wave-mode separation for heterogeneous tilted transverse isotropic media , 2011 .

[12]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[13]  J. Ba,et al.  Polarity reversal correction for elastic reverse time migration , 2012 .

[14]  M. V. Hoop,et al.  Beyond receiver functions: Passive source reverse time migration and inverse scattering of converted waves , 2012 .

[15]  Ari Tryggvason,et al.  Three-dimensional imaging of the P- and S-wave velocity structure and earthquake locations beneath Southwest Iceland , 2002 .

[16]  Cengiz Esmersoy,et al.  Inversion of P and SV waves from multicomponent offset vertical seismic profiles , 1990 .

[17]  Alison E. Malcolm,et al.  Seismic Imaging of Hydraulically-Stimulated Fractures: A Numerical Study of the Effect of the Source Mechanism , 2012 .

[18]  Xiang Xiao,et al.  Local vertical seismic profiling (VSP) elastic reverse-time migration and migration resolution: Salt-flank imaging with transmitted P-to-S waves , 2010 .

[19]  Y. Li.  Image-Guided WEMVA for Azimuthal Anisotropy , 2013 .

[20]  B. Arntsen,et al.  Anisotropic migration velocity analysis using reverse-time migration , 2014 .

[21]  Michael G. Bostock,et al.  Multiparameter two-dimensional inversion of scattered teleseismic body waves 3. Application to the Cascadia 1993 data set , 2001 .

[22]  Paul Sava,et al.  Offset and angle-domain common image-point gathers for shot-profile migration , 2002 .

[23]  Paul Sava,et al.  Wave-Equation Migration Velocity Analysis , 2004 .

[24]  R. James Brown,et al.  Converted‐wave seismic exploration: Applications , 2003 .

[25]  Robert R. Stewart,et al.  Rapid map and inversion of P-SV waves , 1991 .

[26]  Moshe Reshef,et al.  Elastic wave calculations by the Fourier method , 1984 .

[27]  Paul Sava,et al.  Wave-equation migration velocity analysis. I. Theory , 2004 .

[28]  Tariq Alkhalifah,et al.  The space–time domain: theory and modelling for anisotropic media , 2001 .

[29]  William Rodi,et al.  Nonlinear conjugate gradients algorithm for 2-D magnetotelluric inversion , 2001 .

[30]  L. Vinnik,et al.  Heterogeneities in the mantle transition zone from observations of P-to-SV converted waves , 1983 .

[31]  L. P. Vinnik,et al.  Detection of waves converted from P to SV in the mantle , 1977 .

[32]  Subsurface Domain Image Warping by Horizontal Contraction and its Application to Wave-Equation Migration Velocity Analysis , 2013 .

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

[34]  Peng Shen.  An RTM Based Automatic Migration Velocity Analysis in Image Domain , 2012 .

[35]  William W. Symes,et al.  Velocity inversion by differential semblance optimization , 1991 .

[36]  M. L’Heureux,et al.  ENSO prediction one year in advance using western North Pacific sea surface temperatures , 2012 .

[37]  Guojian Shan,et al.  RTM Based Wave Equation Migration Velocity Analysis , 2013 .

[39]  J. Stoer,et al.  Introduction to Numerical Analysis , 2002 .

[40]  Michael G. Bostock,et al.  Multiparameter two-dimensional inversion of scattered teleseismic body waves 1. Theory for oblique incidence , 2001 .

[41]  William W. Symes,et al.  Automatic velocity analysis via shot profile migration , 2008 .

[42]  Y. Li,et al.  VTI Migration Velocity Analysis Using RTM , 2012 .