Reverse time migration for microseismic sources using the geometric mean as an imaging condition

ABSTRACTTime reversal is a powerful tool used to image directly the location and mechanism of passive seismic sources. This technique assumes seismic velocities in the medium and propagates time-reversed observations of ground motion at each receiver location. Assuming an accurate velocity model and adequate array aperture, the waves will focus at the source location. Because we do not know the location and the origin time a priori, we need to scan the entire 4D image (3D in space and 1D in time) to localize the source, which makes time-reversal imaging computationally demanding. We have developed a new approach of time-reversal imaging that reduces the computational cost and the scanning dimensions from 4D to 3D (no time) and increases the spatial resolution of the source image. We first individually extrapolate wavefields at each receiver, and then we crosscorrelate these wavefields (the product in the frequency domain: geometric mean). This crosscorrelation creates another imaging condition, and focusi...

[1]  G. Beroza,et al.  Non-volcanic tremor and low-frequency earthquake swarms , 2007, Nature.

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

[3]  Charles F. Richter,et al.  Magnitude and energy of earthquakes , 2010 .

[4]  Jose Pujol,et al.  Earthquake Location Tutorial: Graphical Approach and Approximate Epicentral Location Techniques , 2004 .

[5]  G. Beroza,et al.  Complex evolution of transient slip derived from precise tremor locations in western Shikoku, Japan , 2007 .

[6]  Jamieson H Steidl,et al.  Report on the August 2012 Brawley Earthquake Swarm in Imperial Valley, Southern California , 2013 .

[7]  G. Arfken Mathematical Methods for Physicists , 1967 .

[8]  Haruko Sekiguchi,et al.  Rupture process of the 2011 Tohoku‐Oki mega‐thrust earthquake (M9.0) inverted from strong‐motion data , 2011 .

[9]  Ted Urbancic,et al.  Microseismic moment tensors: A path to understanding frac growth , 2010 .

[10]  E. Tessmer,et al.  Reverse modelling for seismic event characterization , 2005 .

[11]  Paul Sava,et al.  Wave‐equation migration velocity analysis. II. Subsalt imaging examples , 2004 .

[12]  Roel Snieder,et al.  Focusing of elastic waves for microseismic imaging , 2014 .

[13]  G. Papanicolaou,et al.  Theory and applications of time reversal and interferometric imaging , 2003 .

[14]  K. Sabra,et al.  Detecting remotely triggered temporal changes around the Parkfield section of the San Andreas fault , 2010 .

[15]  J. Claerbout Toward a unified theory of reflector mapping , 1971 .

[16]  Paul Sava,et al.  Isotropic angle-domain elastic reverse-time migration , 2008 .

[17]  Brad Artman,et al.  Source location using time‐reverse imaging , 2010 .

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

[19]  Jeroen Tromp,et al.  Time reversal location of glacial earthquakes , 2008 .

[20]  George A. McMechan,et al.  Imaging conditions for prestack reverse-time migration , 2008 .

[21]  H. Kao,et al.  The Source‐Scanning Algorithm: mapping the distribution of seismic sources in time and space , 2004 .

[22]  Paul Sava,et al.  Probabilistic microearthquake location for reservoir monitoring , 2010 .

[23]  Moon Gi Kang,et al.  Super-resolution image reconstruction: a technical overview , 2003, IEEE Signal Process. Mag..

[24]  Paul Sava Micro-earthquake monitoring with sparsely sampled data , 2011 .

[25]  Robust 3D scalar imaging condition for elastic RTM , 2017 .

[26]  Paul Sava,et al.  Extended imaging conditions for wave‐equation migration , 2011 .

[27]  T. Moser,et al.  Diffraction imaging by focusing‐defocusing: An outlook on seismic superresolution , 2004 .

[28]  F. Poletto,et al.  Synthesis and composition of virtual-reflector (VR) signals , 2010 .

[29]  Andrew W. Hill,et al.  Beyond the dots in the box: Microseismicity-constrained fracture models for reservoir simulation , 2010 .

[30]  Paul Sava,et al.  Angle-domain common-image gathers by wavefield continuation methods , 2003 .

[31]  Toshifumi Matsuoka,et al.  Acceleration of computation speed for elastic wave simulation using a Graphic Processing Unit , 2011 .

[32]  N. Nakata,et al.  Introduction to microseismic source mechanisms , 2015 .

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

[34]  Clifford H. Thurber,et al.  Double-Difference Tomography: The Method and Its Application to the Hayward Fault, California , 2003 .

[35]  Michael Fehler,et al.  Petroleum reservoir characterization using downhole microseismic monitoring , 2010 .

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

[37]  Mathias Fink,et al.  Time-reversal acoustics in complex environments , 2006 .

[38]  I. Vasconcelos,et al.  Imaging condition for nonlinear scattering-based imaging: Estimate of power loss in scattering , 2012 .

[39]  F. Waldhauser,et al.  A Double-Difference Earthquake Location Algorithm: Method and Application to the Northern Hayward Fault, California , 2000 .