Automated seismic waveform location using multichannel coherency migration (MCM)–I: theory

With the proliferation of dense seismic networks sampling the full seismic wavefield, recorded seismic data volumes are getting bigger and automated analysis tools to locate seismic events are essential. Here, we propose a novel Multichannel Coherency Migration (MCM) method to locate earthquakes in continuous seismic data and reveal the location and origin time of seismic events directly from recorded waveforms. By continuously calculating the coherency between waveforms from different receiver pairs, MCM greatly expands the available information which can be used for event location. MCM does not require phase picking or phase identification, which allows fully automated waveform analysis. By migrating the coherency between waveforms, MCM leads to improved source energy focusing. We have tested and compared MCM to other migration-based methods in noise-free and noisy synthetic data. The tests and analysis show that MCM is noise resistant and can achieve more accurate results compared with other migration-based methods. MCM is able to suppress strong interference from other seismic sources occurring at a similar time and location. It can be used with arbitrary 3D velocity models and is able to obtain reasonable location results with smooth but inaccurate velocity models. MCM exhibits excellent location performance and can be easily parallelized giving it large potential to be developed as a real-time location method for very large datasets.

[1]  G. McMechan,et al.  Elastic full-waveform inversion for earthquake source parameters , 1996 .

[2]  Matthias Ohrnberger,et al.  Automatic hypocenter determination of volcano induced seismic transients based on wavefield coherence — an application to the 1998 eruption of Mt. Merapi, Indonesia , 2001 .

[3]  John C. VanDecar,et al.  Determination of teleseismic relative phase arrival times using multi-channel cross-correlation and least squares , 1990 .

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

[5]  George A. McMechan,et al.  Determination of source parameters by wavefield extrapolation , 1982 .

[6]  Sanyi Yuan,et al.  Directional complex-valued coherence attributes for discontinuous edge detection , 2016 .

[7]  Enrico Priolo,et al.  Automated microseismic event location using Master-Event Waveform Stacking , 2016, Scientific Reports.

[8]  Mickael Tanter,et al.  Time-reversed acoustics , 2000 .

[9]  Claire Birnie,et al.  Analysis and models of pre-injection surface seismic array noise recorded at the Aquistore carbon storage site , 2016 .

[10]  G. McMechan,et al.  Source-Parameter Estimation by Full Waveform Inversion in 3D Heterogeneous, Viscoelastic, Anisotropic Media , 2001 .

[11]  Alberto Michelini,et al.  Continuous Kurtosis‐Based Migration for Seismic Event Detection and Location, with Application to Piton de la Fournaise Volcano, La Réunion , 2014 .

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

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

[14]  David Halliday,et al.  Seismic interferometry, surface waves and source distribution , 2008 .

[15]  J. Hunziker,et al.  Seismic interferometry by crosscorrelation and by multidimensional deconvolution: a systematic comparison , 2010 .

[16]  Brian L. N. Kennett,et al.  Automatic Phase-Detection and Identification by Full Use of a Single Three-Component Broadband Seismogram , 2000 .

[17]  V. Oye,et al.  Automated microearthquake location using envelope stacking and robust global optimization , 2010 .

[18]  H. Kao,et al.  FAST TRACK PAPER: Rapid identification of earthquake rupture plane using Source-Scanning Algorithm , 2007 .

[19]  J. Verdon,et al.  Influence of a velocity model and source frequency on microseismic waveforms: some implications for microseismic locations , 2013 .

[20]  H. Kao,et al.  Delineating complex spatiotemporal distribution of earthquake aftershocks: an improved Source-Scanning Algorithm , 2012 .

[21]  S. Minkoff,et al.  Microseismic event estimation in noisy data via full waveform inversion , 2015 .

[22]  P. Podvin,et al.  Finite difference computation of traveltimes in very contrasted velocity models: a massively parallel approach and its associated tools , 1991 .

[23]  R. Guyer,et al.  Tremor source location using time reversal: Selecting the appropriate imaging field , 2009 .

[24]  R. V. Allen,et al.  Automatic phase pickers: Their present use and future prospects , 1982 .

[25]  Simone Cesca,et al.  Full Waveform Seismological Advances for Microseismic Monitoring , 2015 .

[26]  Mathias Fink,et al.  Time‐reversal imaging of seismic sources and application to the great Sumatra earthquake , 2006 .

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

[28]  Frederik Tilmann,et al.  Coalescence microseismic mapping , 2013 .

[29]  B. Steiner Time reverse modeling of low-frequency tremor signals , 2009 .

[30]  Stavros M. Panas,et al.  PAI-S/K: A robust automatic seismic P phase arrival identification scheme , 2002, IEEE Trans. Geosci. Remote. Sens..

[31]  Francesco Grigoli,et al.  Automated Seismic Event Location by Travel‐Time Stacking: An Application to Mining Induced Seismicity , 2013 .

[32]  Aldo Zollo,et al.  Automated seismic event location by waveform coherence analysis , 2014 .

[33]  Jiwei Liu,et al.  Seismic Waveform Classification and First-Break Picking Using Convolution Neural Networks , 2018, IEEE Geoscience and Remote Sensing Letters.

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

[35]  Karl A. Fox,et al.  Methods of Correlation and Regression Analysis, Linear and Curvilinear , 1959 .

[36]  Leo Eisner,et al.  Simultaneous microseismic event localization and source mechanism determination , 2015 .

[37]  S. Yuan,et al.  Microseismic Full Waveform Modeling in Anisotropic Media with Moment Tensor Implementation , 2018, Surveys in Geophysics.

[38]  D. Angus,et al.  Automated seismic waveform location using Multichannel Coherency Migration (MCM)—II. Application to induced and volcano-tectonic seismicity , 2018, Geophysical Journal International.

[39]  Sanyi Yuan,et al.  Simultaneous multitrace impedance inversion with transform-domain sparsity promotionStructures — Exploring impedance inversion , 2015 .