zotero/storage/QF6GYQDY/.zotero-ft-cache

1229 lines
238 KiB
Plaintext
Raw Permalink Normal View History

2024-08-27 21:48:20 -05:00
REVIEW ARTICLE
10.1029/2019RG000667
Key Points: • We present a comprehensive
summary of the current status of waveformbased seismic source location methods • We illustrate applications of waveformbased methods to natural and induced seismicity ranging from laboratory to field and regional scales • This study outlines the advances and challenges associated with waveformbased location methods
Correspondence to: J. Tan, tanjingqiang@aliyun.com
Citation: Li, L., Tan, J., Schwarz, B., Staněk, F., Poiata, N., Shi, P., et al. (2020). Recent advances and challenges of waveformbased seismic location methods at multiple scales. Reviews of Geophysics, 58, e2019RG000667. https://doi.org/10.1029/2019RG000667
Received 27 SEP 2019 Accepted 11 JAN 2020 Accepted article online 15 JAN 2020 Corrected 27 APR 2020
This article was corrected on 27 APR 2020. See the end of the full text for details.
©2020. The Authors. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.
Recent Advances and Challenges of WaveformBased Seismic Location Methods at Multiple Scales
Lei Li1,2,3 , Jingqiang Tan1,2,3 , Benjamin Schwarz4 , František Staněk5 , Natalia Poiata6,7 , Peidong Shi8 , Leon Diekmann9, Leo Eisner10 , and Dirk Gajewski11
1Key Laboratory of Metallogenic Prediction of Nonferrous Metals and Geological Environment Monitoring, Ministry of Education, Central South University, Changsha, China, 2Hunan Key Laboratory of Nonferrous Resources and Geological Hazard Exploration, Changsha, China, 3School of Geosciences and InfoPhysics, Central South University, Changsha, China, 4GFZ German Research Centre for Geosciences, Potsdam, Germany, 5Institute of Rock Structure and Mechanics, Czech Academy of Sciences, Prague, Czech Republic, 6Institut de Physique du Globe de Paris, Université de Paris, UMR 7154 CNRS, Paris, France, 7National Institute for Earth Physics, Magurele, Romania, 8Institut des Sciences de la Terre (ISTerre), Université Grenoble Alpes, Saint Martin d'Hères, France, 9Department of Earth Sciences, Utrecht University, Utrecht, Netherlands, 10Seismik s.r.o., Prague, Czech Republic, 11Institute of Geophysics, University of Hamburg, Hamburg, Germany
Abstract Source locations provide fundamental information on earthquakes and lay the foundation for
seismic monitoring at all scales. Seismic source location as a classical inverse problem has experienced significant methodological progress during the past century. Unlike the conventional traveltimebased location methods that mainly utilize kinematic information, a new category of waveformbased methods, including partial waveform stacking, time reverse imaging, wavefront tomography, and full waveform inversion, adapted from migration or stacking techniques in exploration seismology has emerged. Waveformbased methods have shown promising results in characterizing weak seismic events at multiple scales, especially for abundant microearthquakes induced by hydraulic fracturing in unconventional and geothermal reservoirs or foreshock and aftershock activity potentially preceding tectonic earthquakes. This review presents a comprehensive summary of the current status of waveformbased location methods, through elaboration of the methodological principles, categorization, and connections, as well as illustration of the applications to natural and induced/triggered seismicity, ranging from laboratory acoustic emission to field hydraulic fracturinginduced seismicity, regional tectonic, and volcanic earthquakes. Taking into account recent developments in instrumentation and the increasing availability of more powerful computational resources, we highlight recent accomplishments and prevailing challenges of different waveformbased location methods and what they promise to offer in the near future.
Plain language summary Earthquakes are a common physical phenomenon involving ground
shaking and rupturing of the surface of the Earth. In addition to the wellknown catastrophic tectonic earthquakes, similar vibration sources also appear at various scales in engineering fields, such as acoustic emissions resulting from microcracks in building walls and bridges, rock bursts in mines, microseismic events generated by mining and fluid injection/extraction, and microseisms caused by crustal activity. Seismic information provides a powerful tool for geophysical and engineering surveys. The source location describes the spatial and temporal extent of an earthquake and lays the foundation for seismic monitoring. Seismic location methods have made significant progress over the last century. Specifically, a category of new waveformbased location methods has emerged as a counterpart of conventional traveltimebased inversion. These methods directly utilize the notion of a wavefield and, very similar to an optical lens, aim at spatially focusing a source's emitted energy. Waveformbased methods have provided robust and effective source location results at various scales. We summarize the development history and current state of waveformbased location methods and discuss the advantages and challenges through their applications for seismic source location at multiple scales.
1. Introduction
Seismic source location is a classical inverse problem in seismology and geoengineering. A seismic source location specifies the origin time and spatial coordinates of seismic energy release generated by a seismic source (earthquake), and it also refers to the process of locating a seismic event (e.g., Lomax et al., 2009).
LI ET AL.
1 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Reviews of Geophysics
10.1029/2019RG000667
LI ET AL.
Figure 1. Scaling of corner frequency and source radius versus moment magnitude and detectable distance associated to differentscale seismic events. Several assumptions are adopted here in this schematic diagram, for example, peak energy approximation for corner frequency and pure shear source model for source radius, and average attenuation is assumed for magnitude calculation. Detectable distance is roughly indicated as a critical distance beyond which the event of a given magnitude cannot be detected with a generic receiver. The illustration considers values in a general sense and for an illustrative purpose only; different geological settings may result in different relationships between these parameters. Modified after Hardy (2003), Kwiatek et al. (2011), Zoback and Gorelick (2012), and Eaton et al. (2018).
This review concerns only the location of individual events under the point source assumption; the contents of location and imaging of a finite source or fault are not covered. Observed seismicity is either naturally occurring in seismogenic areas (e.g., Shearer, 2009) or is triggered/induced by anthropogenic/industrial operations, including mining, water reservoirs impoundment, fluid injection/extraction operations, oil and gas field exploration, and geothermal exploitation (McGarr et al., 2002; Stump, 2002; Zang et al., 2014; Elsworth et al., 2016; Grigoli et al., 2017; Keranen & Weingarten, 2018; Foulger et al., 2018; L. Li et al., 2019; Zoback & Kohli, 2019). In particular, several cases of felt earthquakes with Mw > 3 related with hydraulic fracturing or fluid injection have been reported worldwide, including the United States, the United Kingdom, Western Canada, and Southwestern China (Atkinson et al., 2016; Bao & Eaton, 2016; Clarke et al., 2014; Holland, 2013; Lei et al., 2017; Meng et al., 2019), as well as a recently reported event of Mw 5.5 near an enhanced geothermal system site in South Korea (Grigoli et al., 2018; Kim et al., 2018). In addition to the aforementioned field or exploration scale scenarios, seismicity also occurs at broader scales, varying from laboratory to local and global scales (Eaton, 2018; Hardy, 2003; Shearer, 2009). In earthquake seismology (cover local to global scales), seismic waves and related attributes (e.g., first arrival time) of tectonic earthquakes or tectonic and volcanic tremors are used to estimate velocitydepth models for the Earth's crust and the deeper interior and provide earthquake location and its focal mechanism estimation. Seismic activity associated with millimeterscale fractures and kilometerscale faults share highly similar physical processes (Goodfellow & Young, 2014; Kwiatek et al., 2011), although the dominant frequency content and energy release (magnitude) of laboratory acoustic emissions (AEs) are different from those of globally felt tectonic earthquakes. An illustration of differences in scales of seismic activities with different ranges of corner frequency, magnitude, source radius, and detectable distances is shown in Figure 1. The seismic source radius, magnitude, and the detectable distance or range decrease with corner frequency. Laboratory studies are generally based on rock compression or fracturing tests and high frequency AE monitoring, and they contribute to revealing the smallscale deformation and failure process of brittle rock (Lockner et al., 1991; Scholz, 1968). Seismic waves are used for delineating the nearsurface or subsurface (down to 10 km) structures, which are associated with exploration and development of metals,
2 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Reviews of Geophysics
10.1029/2019RG000667
LI ET AL.
Figure 2. Historical evolution of seismic source location methods, with the most significant conceptual cornerstones highlighted. Red colors indicate the direct involvement of waveform information in the location process. Only several prominent methods are explicitly mentioned; closer investigation is provided in the main text.
minerals, and hydrocarbons, while regional and global earthquakes are used in earthquake seismology to characterize fault zones and the deeper interior (down to 100 km or more; Gibowicz & Kijko, 1994; Stein & Wysession, 2003; Yilmaz, 2015; Ikelle & Amundsen, 2018). For applications at all scales, earthquake source location plays one of the most fundamental roles in the workflow. It provides fundamental information on the earthquake that creates detected seismic waves and lays the foundation for subsequent seismological analysis and characterization, such as magnitude calculation, seismicity pattern delineation, source zone monitoring, source mechanism, and stress field estimation.
Seismic source location problem has experienced significant methodological progress during the past century (e.g., Thurber & Rabinowitz, 2000; Lomax et al., 2009; see also Figure 2 for an overview). The earliest seismic source location method was based on manual triangulation (e.g., Milne, 1886; Pujol, 2004), making use of the arrival time differences between compressional wave (P wave) and shear wave (S wave) for calculating the distance between a source and a station. Once arrival time differences for three or more stations are available, the hypocenter can be determined by intersecting equidistant lines from the stations. This procedure is also known as the graphical approach. Alternatively, Geiger (1912) introduced an iteratively linearized traveltime inversion that requires phase picking, and this algorithm laid the foundation for subsequent numerical location methods. Since the 1980s, the digital seismic recordings and computer technologies enabled seismic source location being processed based on computational or numerical methods, rather than manual approaches (e.g., Sambridge & Kennett, 1986; Thurber & Rabinowitz, 2000). Latter, nonlinear methods that did not involve linearization emerged, and direct (including regular or stochastic) searches, such as the genetic algorithm (Kennett & Sambridge, 1992), and the Monte Carlo searches (Sambridge & Mosegaard, 2002) over the target space were utilized to determine the source location. All these methods use picked arrival times, and in principle, they involve searching for the minimum misfit between the theoretical and observed traveltimes by applying either an inversion or a grid search. Traveltimebased location methods require phase picking and only make use of the traveltime information. They are also known as picking or raybased methods (Z. Li & van der Baan, 2016; Pesicek et al., 2014). Many modifications have been introduced to improve the performance of the traveltimebased methods, most of which lie in the construction of misfit function and inversion strategies, such as joint hypocenter location method (e.g., Douglas, 1967; Pujol, 1992), relative location method (e.g., Fitch, 1975; Grechka et al., 2015; Spence, 1980), plain grid search (e.g., Eisner et al., 2010), master station or stationpair double difference method (e.g., Font et al., 2004; Zhang et al., 2010; Zhou, 1994), joint location and velocity inversion (e.g., Block et al., 1994; Jansky et al., 2010; Z. Zhang et al., 2017; Diekmann et al., 2019), double difference relocation method (Waldhauser & Ellsworth, 2000; H. Zhang
3 of 47
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
& Thurber, 2003; Hauksson & Shearer, 2005; Waldhauser & Schaff, 2008; Kwiatek et al., 2013), and clusterbased relocation methods (G. Lin et al., 2007; P. Shearer et al., 2005; Matoza et al., 2013; Trugman & Shearer, 2017).
As an alternative of the conventional traveltimebased methods, modern waveformbased seismic source location methods do not require phase identification nor phase picking and can detect and locate events with low signaltonoise ratio (SNR) by exploiting migration and stacking techniques commonly used in exploration seismology (e.g., Beskardes et al., 2018; Cesca & Grigoli, 2015; Gajewski et al., 2007; Pesicek et al., 2014; Zhebel & Eisner, 2015). The basic principle of waveformbased methods is locating the source by focusing or reconstructing the source energy into discrete grid of points with a certain migration or imaging operator, which uses traveltime information. It is worth noting that there are several different terminologies for the waveformbased location methods, such as waveform migration or stacking, backprojection imaging, beamforming, and coherence scanning (Cesca & Grigoli, 2015; Gajewski et al., 2007; Maxwell, 2014). Here, for simplicity and avoiding ambiguity, the term “waveformbased” is used throughout the manuscript for the type of location methods using waveforms and/or waveform techniques, not only exploiting the timing of seismic phase arrivals as traveltimebased methods do. For waveformbased methods, no phase picking is required. Every possible source location is firstly considered through stacking and imaging process; then, a detection and picking criterion (or imaging condition) is utilized to select the final location. The procedure can be summarized as “locate and then pick” (Eaton, 2018). In contrast, the general procedure of traveltime based methods is phase picking first and then locating the events through inversion strategies, namely, “pick and then locate.”
The origin of waveformbased methods can be traced back to around 1990s. In 1982, McMechan proposed a pioneering recipe for locating earthquake events through wavefield extrapolation (McMechan, 1982). Then, seismic emission tomography was developed by utilizing maxima of the semblance (a measure of coherence) over space, time, and channels to detect and locate microseismic events (Kiselevitch et al., 1991). An advance on this method, socalled passive seismic emission tomography, has also been proposed to image microseismic emission energy with surface monitoring arrays (Duncan, 2005). This method efficiently produced promising results of source location. The basic idea behind these approaches is the same as diffraction stacking in exploration seismology, namely, stacking the waveforms from individual receivers along oneway traveltime moveout for all possible source locations and origin times. In addition to the Kirchhoff diffraction stacking operator, crosscorrelation stacking operator is another wellestablished, interferometrictype migration operator for source location (Schuster et al., 2004). Crosscorrelation stacking shares the same principle of exploiting differential traveltimes at pairwise stations from common events, with masterstation method and stationpair method mentioned above. This method has been applied to microseismic monitoring in the oil and gas industry (Eisner et al., 2008; Grandi & Oates, 2009; L. Li et al., 2015; Trojanowski & Eisner, 2017; Shi et al., 2019), volcanic tremor data (Droznin et al., 2015; K. L. Li et al., 2017; K. L. Li et al., 2017), mininginduced seismicity (Dales, 2018; L. Li et al., 2018), and slow (tectonic tremor and lowfrequency earthquakes) and regular earthquake (Poiata et al., 2016; Ruigrok et al., 2017; Shi et al., 2019). In this study, we will call it interferometrictype method. Most of these methods use only direct P waves, as the information of the velocity models is usually limited and does not allow use of other arrivals. Therefore, they are also known as partial waveform stacking (PWS; e.g., Pesicek et al., 2014). Beside aforementioned waveform stacking techniques, time reverse imaging (TRI) and full waveform inversion (FWI) are also wellknown waveformbased methods for source location. TRI locates the seismic source by performing a numerical backpropagation of the recorded full waveforms through a fixed velocity model. This approach was initially proposed by McMechan (1982) and was further developed by, for example, Gajewski and Tessmer (2005), Larmat et al. (2008), Artman et al. (2010), and T. Zhu (2014). The basic principle of FWI is to find a model (e.g., velocities and source parameters) that can generate synthetic waveforms that match the observed waveforms best. Though FWI involves high nonlinearity and requires expensive computational effort, it is increasingly used not only for subsurface velocity estimation in exploration seismology (e.g., Virieux & Operto, 2009), but also for seismic source location (e.g., Kaderli et al., 2015; Shekar & Sethi, 2019). More recently, a largely data driven technique, known as wavefront tomography (WT), combines kinematic information estimated
4 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Reviews of Geophysics
10.1029/2019RG000667
LI ET AL.
Figure 3. Physical and methodological classification of different prominent location methods that are discussed in more detail in the text. Like in Figure 2, red colors indicate a stronger (more explicit) utilization of fullwaveform information in the location procedure. Taking into account implementational aspects, a precise positioning of methods is rarely achievable. Possible deviations from the chosen positions are denoted by arrows. The color of these arrows indicates how means of data abstraction or data availability let traditionally separated techniques approach each other.
from wavefront attributes with timereverse modeling to jointly estimate source location and the velocity structure (Diekmann et al., 2019).
So far, there is no universal name or categorization for these relatively modern location methods involving waveform techniques. There are alternative ways of categorizing seismic location methods, besides the concise but rough categorization as traveltime and waveformbased methods (see Figures 2 and 3 for more details). For instance, conventional traveltime inversion involves data domain operation only, and waveform stacking techniques generally determine a source location from imaging function in the image domain, while FWI is an example of hybrid domain methods. Another possible classification strategy is based on the distinction between raybased modeling (e.g., traveltime inversion and stacking methods involving theoretical traveltime) and full waveform modeling (e.g., TRI and FWI). Without loss of generality, here we name those methods simply as waveformbased methods, though some of them also adopt kinematic information, for example, traveltime and wavefront attributes, and involve both data and image domain operations.
Table 1 lists studies and applications of waveformbased location methods applied to seismic sources observed at different scales. A detailed description on transition of traveltimebased to waveformbased methods, and the increasing interest of waveformbased methods at multiple scales is presented in the next section. So far, PWS is still the dominant approach for data at local to global scales. Although TRI and FWI have already been developed for more than a decade, the study concerning their applicability and potentials in seismic source location is still ongoing.
In this work, we aim to provide insights of waveformbased location methods by a comprehensive review of their recent advances and challenges. The methodological theories for different categories of waveform based methods are followed by representative applications ranging from laboratory AE to field hydraulic fracturinginduced seismicity, regional tectonic, and volcanic earthquakes. These examples illustrate the success of applications to various cases of seismic activity covering wide range of scales. We also discuss current challenges and perspectives associated with waveformbased methods. By elaborating the methodological principles, presenting categorization, and connections, as well as illustrating their applications at
5 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Reviews of Geophysics
10.1029/2019RG000667
Table 1 Summary of Studies Using WaveformBased Methods for Seismic Source Location at Different Scales
Scale
Method
Application
Uncertainty estimation
References
Laboratory & small Field & local
PWS TRI PWS
TRI
WT FWI
Regional & global
PWS
TRI FWI
Laboratory AE events HFinduced AE events Laboratory AE events HFinduced events
Mininginduced events
Microearthquakes Synthetic events
Microtremors HFinduced events Synthetic events and
earthquakes Synthetic events
HFinduced events Regular earthquakes and
tremors
Volcanictectonic earthquakes
Aftershocks Glacial earthquakes Synthetic events Synthetic events Gas extraction induced
earthquakes
/ Jackknife resampling /
Gaussian distribution /
Probability density function
Percentage and standard deviation
Jackknife resampling
Jackknife and bootstrap resampling
/
/ / /
/
Velocity uncertainty Percentage /
Probability density function
Percentage and standard deviation
/ Percentage / / / Error ellipsoid
Gu (2016); C. Zhang et al. (2019) LópezComino et al. (2017) Saenger et al. (2011); Nguyen et al. (2018);
T. Zhu et al. (2019) Drew et al. (2005, 2013) Gajewski et al. (2007); Rentsch et al. (2007);
Chambers et al. (2014); Zeng et al. (2014); Zhebel and Eisner (2015); L. Li et al. (2015); Liang et al. (2016); Trojanowski and Eisner (2017); S. Wu et al. (2018) Anikiev et al. (2014)
Gharti et al. (2010)
Grigoli et al. (2013); L. Li et al. (2018); L. Li, Tan, Xie, et al. (2019)
Pesicek et al. (2014)
Gajewski and Tessmer (2005); Douma et al. (2015); Hansen and Schmandt (2015); Z. Li and van der Baan (2016)
Steiner et al. (2008) T. Zhu (2014); T. Zhu et al. (2019) Diekmann et al. (2019)
Michel and Tsvankin (2014, 2017); Kaderli et al. (2018); Wang and Alkhalifah (2018); Igonin and Innanen (2018); Shekar and Sethi (2019)
Huang et al. (2017a, 2017b) Kao and Shan (2004, 2007) Ishii et al. (2007); Kiser and Ishii, (2013);
Grigoli et al. (2016); Poiata et al. (2016, 2018); K. L. Li, Sgattoni, et al. (2017); K. L. Li, Sadeghisorkhani, et al. (2017); Neale et al. (2018) Hendriyana et al. (2018)
Langet et al. (2014)
Shi, Nowacki, et al. (2019) Beskardes et al. (2018) Larmat et al. (2008) Werner and Saenger (2018) Y. Wu and McMechan (1996) Willacy et al. (2019)
Note. PWS = partial waveform stacking; TRI = time reverse imaging; WT = wavefront tomography; FWI = full waveform inversion; HF = hydraulic fracturing; AE = acoustic emission. The slash (/) means the location uncertainty is not considered or estimated.
LI ET AL.
multiple scales, this review is aimed on building a framework for better understanding the idea of waveform based seismic source location and conducting further research.
2. The Rise of WaveformBased Location Methods
Starting from graphical approaches like the circle method, seismic source location has significantly matured in the 1960s, spurred by the rise of modernday computers. Figure 2 illustrates the historical evolution and
6 of 47
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
major milestones of quantitative location methods. With computational infrastructure available, the first practical applications of Geiger's least squares traveltime inversion (Geiger, 1912) became possible, leading to a rapid development of quantitative traveltimebased methods. As a result, nonlinear extensions, capable of honoring, for example, seismic anisotropy, were subsequently introduced in the 1980s. At the same time, owing to the increasing availability of large data volumes, relational techniques, like the doubledifference, the masterstation, or the masterevent method, were developed, which allowed for improved earthquake location estimates in wellinstrumented or seismically active areas. After the 1980s, the installation of dense station networks and broadband instruments then sparked the ongoing trend of utilizing fullwaveform information for seismic location.
As a more or less natural extension of established traveltime inversion, FWI aiming at utilizing amplitude and phase information in the location framework was introduced. At roughly the same time, following the first experiments in the 1980s and early 1990s, the 2000s saw a rapid transfer of waveformbased imaging techniques from hydrocarbon exploration. Tailored to the location problem, methods like coherencebased WT or PWS (Kirchhoff migration) originally emerged in a controlledsource context and were built on the diffraction concept. Likewise, other emerging arrayprocessing techniques such as TRI or matchedfield processing have their origin in the adjacent fields of (hydro)acoustics. Both, FWI and the aforementioned source imaging techniques have in common that their successful application and further development continuously relies on improvements of data acquisition systems. Although in this paper we focus on emerging waveformbased techniques, it should be noted that traveltimebased methods remain a very trusted and continuously developing approach for seismic source location, with improvements concentrating on the reliable and automated identification and picking of seismic phases.
Because of the sparse distribution of early seismological monitoring stations, the single or small number channel techniques, such as those phasepickingbased methods, were the favorite methods. Traveltime based location methods have been continuously improved and undoubtedly dominated the location routines during the 20th century. Phase picking required by these techniques is time consuming, error prone, and sometimes subjective, especially for weak aftershocks or microseismic data, which typically have relatively low SNR (e.g., Chambers et al., 2010; Duncan & Eisner, 2010; Liao et al., 2012). Though advanced and automated procedures can ensure an efficient picking process (e.g., Akram & Eaton, 2016; Grigoli et al., 2018; Ross et al., 2018; Withers et al., 1998), robust phase identification and picking are still challenging for low SNR events, and missed arrival times and picking errors may directly cause location errors. Other nonnegligible limitations of most traveltime or raybased methods consist of their incapability of utilizing multichannel coherency (Poiata et al., 2016; Shi, Angus, et al., 2019), or considering secondary features (like scattering) or finitefrequency effects (Dahlen et al., 2000; H. Zhang & Thurber, 2003).
The limitations imposed by a low SNR on the picking process can be considered one of the main incentives for the development of alternative location methods. The advent of dense local networks for aftershocks and seismic swarm observation, as well as the interest of industry in passive seismic methods for reservoir and fracture monitoring with even denser networks naturally, induced the interest in exploiting the redundancy in the recorded seismic data. This leads to the rise of multichannel approaches for seismic source location, which are well known in reflection seismology and formed the foundation for the “waveformbased methods” (beamforming has been already used in seismology for some time). Specifically, an increasing demand of noiserobust and efficient location approaches arose since early 2000s, due to the rapid development of microseismic monitoring for hydraulic fracturing in unconventional oil and gas fields. Recent interest in microseismic monitoring has been driven by the global expansion in the development of unconventional oil and gas resources, especially shale and tight sands from North America. Microseismic monitoring is a passive seismic technique enabling analysis of microearthquakes, which are used to map fracture reactivation/creations, and fluid movements. Microseismic monitoring is used in mining, fluid injection and extraction, hydraulic fracturing, and other geoengineering operations (Grechka & Heigel, 2017; Maxwell, 2014; Shapiro, 2015). Passive seismic monitoring has been applied in the mining industry and geothermal reservoirs for more than 60 years (Gibowicz & Kijko, 1994; Obert & Duvall, 1945). Currently, microseismic monitoring has become frequently applied in stimulations of unconventional resources (Maxwell, 2014; Rutledge & Phillips, 2003). Seismic source location is the
7 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
most crucial step of microseismic processing, and it is needed for subsequent source mechanism inversion and geomechanical analysis. A multitude of studies have been dedicated to the improvement of waveformbased microseismic location during the past decade (e.g., Anikiev et al., 2014; Cesca & Grigoli, 2015; Trojanowski & Eisner, 2017; L. Li et al., 2018; Diekmann et al., 2019; see also Table 1 for a more representative list of references). The main motivation of these studies is to improve locations of the weak microseismic events as the GutenbergRichter (Gutenberg & Richter, 1956) empirical law implies exponential increase of the number of events if we can locate weaker events. Since the ability to detect and locate a larger number of events offers possibility for a better characterization of the stimulated rock, the methodologies are focusing on locating as weak events as possible. This eventually leads to the use of waveformbased location methods allowing to detect and locate weaker events by stacking of the waveforms.
With the increasing development and application of hydraulic fracturing in unconventional oil and gas fields, AE monitoring of laboratory hydraulic fracturing tests has drawn attention. AE location methods are directly adapted from earthquake location, and the traveltime inversion acts as the dominant approach (e.g., Ge, 2003a, 2003b). In recent years, a category of signal or waveformbased analysis techniques has been developed in AE location and source mechanism inversion (Grosse & Ohtsu, 2008; Ishida et al., 2017), benefiting from the improvement of computation, communication, and storage technologies. Both PWS and TRI have shown promising results in locating highfrequency AE events (e.g., LópezComino et al., 2017; Saenger et al., 2011). In addition, densification of seismic networks and increased volumes of recorded seismic data at regional and global scale is another important reason for increased interest in waveformbased methods, which are beneficial for tectonic earthquake location problem. Earthquakes are predominantly related to failures, faults, or tectonic boundaries across all scales of the crust, and the most reliable structural information on these features has been gained through the location of connected seismicity. Waveformbased methods like seismic migration are still the primary processes for the construction of structural images in controlled source studies. Thus, waveformbased location methods promise to eventually provide a more direct, datadriven characterization of source and structure.
Nevertheless, the distinction between traveltimebased and waveformbased methods is not as strict as it appears. Stacking techniques require traveltimes (obtained from the velocity model), and TRI works with velocities as well. WT provides velocities (inversion) and source locations (folding it with the waveform provides a source image). FWI stands out since it is the only technique combining inversion and imaging in one step. The clearest distinction seems to be the singlechannel and multichannel style of approach, though phase picking can also be enhanced by multichannel techniques (through modern techniques like machine learning). It should be acknowledged that continued advances in traveltimebased methods are emerging that we do not discuss here. Upon closer inspection, it is no surprise that PWS underpins most of the successful applications of waveformbased location in the literature. It can be viewed as a direct descendant of Kirchhoff migration, which—owing to its flexibility, efficiency, and robustness—can still be considered the most trusted controlledsource imaging method to date (Etgen et al., 2009).
3. Methodologies
As many factors fuel the current rise of waveformbased seismic source location techniques, a systematic and consistent classification of different location methods is desirable yet challenging. Figure 3 presents an attempt to perform a categorization in physical or informationtheoretical terms, that is, on how explicitly they make use of the full information content encoded in a seismic time series at all considered stations. The four different poles of the displayed diagram are interrelated either via means of abstraction or instrumentation. The source imaging techniques all have their origin in the fields, where spatial wavefield sampling is reliably dense (controlledsource seismology, optics, or acoustics), whereas more indirect but flexible inversion algorithms like traveltime and fullwaveform inversion are less restricted in that regard. TRI and WT represent direct counterparts and can only be distinguished by the degree of data abstraction. Just like a source location alone, traveltimes and wavefronts can be viewed as idealized abstractions and systematically neglect additional information about the source. The utilization of full waveform information in turn, despite the valuable additional details contained, remains practically challenging and still often leads to unreliable results. It is a combination of the strengths, namely, the
8 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Reviews of Geophysics
10.1029/2019RG000667
LI ET AL.
Figure 4. Waveformbased methods locate seismic source through backpropagation of seismic energy radiated by the source. Model parameters (e.g., velocity model and sourcereceiver geometry) are needed for calculation of traveltime in, for example, PWS and for waveform modeling in, for example, TRI and FWI. Migration or imaging operators can transform waveforms in the data domain to spatial image of the source in the image domain, which can be subsequently used to characterize the faults and fractures in the model domain.
insightfulness of waveform information and the stability of data abstractions like traveltimes, that make hybrid approaches like PWS arguably the most successful candidate so far. Indicated by the arrows in Figure 3, improved instrumentation and advances in wavefield reconstruction, as well as intermediate forms of abstraction, such as characteristic functions (CFs), or waveform fitting lead to the different mindsets approaching each other.
The common principle of waveformbased location methods is focusing or reconstructing the seismic energy radiated by the source into a point with a certain migration or imaging operator. From a theoretical perspective, these methods are based on the KirchhoffHelmholtz formula and Huygens' principle, which are fundamental representation theorems in elastodynamics (Kirchhoff, 1883; Pao & Varatharajulu, 1976; Born & Wolf, 1980; Aki & Richards, 2002; B. B. Baker & Copson, 2003). These theorems provide a theoretical framework for seismic imaging and state the elastic wavefield at any point within a subsurface region that can be completely reconstructed through a weighted integration of the wavefield and its normal derivative measured around a closed surface. This is also the theoretical foundation of Kirchhoff migration (e.g., Schneider, 1978), a routinely used processing technique in exploration seismology. As opposed to the forward propagation of seismic waves from the source to the receiver, source location through waveformbased methods can be regarded as backpropagation of seismic energy from the receivers to the source zone (Figure 4). Thus, the waveformbased methods are based on reciprocity theorem in a general sense. Currently, most waveformbased methods provide absolute locations, which is different from relative locations that are inverted based on known master events. Relative locations basically measure the source locations relative to master events and are often used in highprecision relocation procedures. Though a few pioneering studies (see discussion in later parts of section 3.1) on
9 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Reviews of Geophysics
10.1029/2019RG000667
LI ET AL.
Figure 5. Sketch of the principle for partial waveform stacking methods. (a) Kirchhofftype method, which is based on oneway traveltimes from the source to receivers. (b) Interferometrictype method, which is based on differential traveltimes at pairwise stations from common events. As indicated in the migration operators, Kirchhofftype methods search both source location and origin time, while interferometrictype methods stack the crosscorrelograms of waveforms from pairwise receivers, without inverting the decoupled origin time. The generalized equations of partial waveform stacking and those of commonly used characteristic functions (CFs) are summarized in equations (1) and (2).
combining waveform techniques and relative location idea have shown promising results, these new and hybrid methods are still in the early stage of development. Here we focus on waveformbased, absolute location methods only.
In this section, we present a systematic introduction to the theory and development of waveformbased location methods, including PWS, TRI, WT, and FWI. These methods have all in common that they either directly or indirectly utilize waveform information in the location process (cf. Figure 3). Specific aspects of implementation, ranging from the data abstraction through CFs to velocity model building or the generation of traveltime tables, will be discussed in the following section. While PWS has already proven to be successful in challenging environments, other methods, such as WT, TRI, or FWI, either just appeared on the radar or still face significant practical challenges that need to be addressed. Accordingly, the following section emphasizes implementational aspects where recent advances have been made and discusses the potential and challenges of each method in more detail.
3.1. Partial Waveform Stacking
PWS methods are currently the most widely used among waveformbased location methods, especially for detection and location of seismic events at local and regional scales. These methods usually make use of primary phases only. This category mainly includes diffraction stacking and crosscorrelation stacking, which corresponds to Kirchhofftype and interferometrictype stacking, respectively.
Kirchhoff migration is a weighted backprojection of wavefields onto an isochron surface (a surface with equal traveltime) that passes through the scattering object (Esmersoy & Miller, 1989). When the scattering object turns into a point or an edge, the diffraction theory works (e.g., KlemMusatov, 1994; B. Schwarz, 2019). The seismic source can be approximated by a point in most scenarios ranging from laboratory to global scales of seismic monitoring, because the sourcereceiver distance usually exceeds more than five prevailing wavelengths (e.g., see discussion on dominant frequency of microseismic sources in Eisner et al., 2013). Therefore, the diffraction stacking can be naturally applied to source imaging, and numerous studies have provided development and application of Kirchhofftype methods for source location, which are based on stacking of the waveforms from individual receivers along a theoretical traveltime moveout curve (Figure 5a; e.g., Kao & Shan, 2004; T. Baker et al., 2005; Ishii et al., 2005; Gajewski et al., 2007; Zhebel & Eisner, 2015). In earthquake seismology, Kao and Shan (2004) firstly proposed the sourcescanning algorithm to identify seismic event location and origin time with a brightness function, which is constructed by stacking the recorded seismograms along a moveout theoretical traveltimes. Gajewski et al. (2007) proposed to use diffraction stacking for source location and demonstrated its automatism and feasibility with numerical examples in passive seismic scenarios. They indicated that the location uncertainty was highly dependent on the
10 of 47
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
frequency content of the event. This Kirchhofftype method is also referred to as Kirchhoff reconstruction (T. Baker et al., 2005), backprojection (Kiser & Ishii, 2013), and waveform migration or stacking (Gharti et al., 2010).
Another wellestablished and interferometrictype migration operator for source location is the cross correlation stacking (Schuster et al., 2004), which originated from seismic interferometric imaging and is based on differential traveltimes at pairwise stations from common events (Figure 5b; e.g., Eisner et al., 2008; L. Li et al., 2015). Instead of stacking waveforms along oneway traveltime curves as in diffraction stacking, crosscorrelograms along differential traveltime curves are stacked taking into account the correlation among the receiver channels. Actually, the idea of utilizing differential traveltimes on pairwise receivers from common events has emerged as the masterstation method in earthquake seismology since the 1990s (e.g., H. Zhou, 1994). Font et al. (2004) extended the masterstation method to maximum intersection method by incorporating differential traveltimes from all unique receiver pairs. Analogous to double difference method (Waldhauser & Ellsworth, 2000), Zhang et al. (2010) proposed the stationpair double difference method and used it to locate nonvolcanic tremors successfully. Recently, several studies implemented this idea into waveformbased methods with crosscorrelation techniques from seismic interferometry (e.g., Dales et al., 2017; Droznin et al., 2015; Grandi & Oates, 2009; Poiata et al., 2016; Ruigrok et al., 2017; Shi, Angus, et al., 2019). By including higherorder cross correlations, K. L. Li, Sadeghisorkhani, et al. (2017) and K. L. Li, Sgattoni, et al. (2017) extended the previous single cross correlation stacking to a doublecorrelation method and a more generalized highorder correlation method, which showed superior performance in lateral imaging resolution for locating tremors with sparse station networks. Shi, Angus, et al. (2019) proposed a novel multichannel coherency migration method for locating earthquakes using continuous seismic data based on stacking of Pearson correlation coefficients for waveforms from different receiver pairs or groups. Synthetic and field microseismic data sets demonstrated the multichannel coherency could produce improved source energy focusing.
Besides the imaging operator, the PWS methods also differ in preprocessing of the input waveforms and the event detection and location detection criteria. There are various CFs adopted to improve the SNR and compensate the side effects of source radiation patterns; for example, the input data are converted to
1. envelope (Gharti et al., 2010; Liao et al., 2012; Zeng et al., 2014); 2. semblance (Furumoto et al., 1990; W. Zhang & Zhang, 2013; Staněk et al., 2015; C. Zhang et al., 2019); 3. shortterm average to longterm average ratio (STA/LTA; Drew et al., 2013; Grigoli et al., 2014; L. Li et al.,
2018); 4. kurtosis (Langet et al., 2014; Poiata et al., 2016); 5. multichannel coherency (Shi, Angus, et al., 2019).
Different CFs can exhibit different performance of imaging resolution and noise suppression. Cesca and Grigoli (2015) investigated the location performance of different CFs with both synthetic and field mininginduced seismicity and concluded that kurtosis traces are sensitive with respect to noise, envelope traces produce more stable results but lower imaging resolution, and STA/LTA traces have a better overall performance of location resolution and noise resistance. Trojanowski and Eisner (2017) tested the performance of several CFs with diffraction stacking and crosscorrelation stacking operators, showed the semblance and crosscorrelation stacking have better enhancement of the quality of the stacking profiles, and demonstrated that accounting for the source mechanism by polarization correction improves most efficiently the capability of noise suppression and retrieving weak microseismic events. Beskardes et al. (2018) compared the performance of envelope, STA/LTA, and kurtosis with raw waveforms for local earthquakes detection and location. The comparative study showed for an aftershock sequence that stacking of kurtosis waveforms produced the most robust results for detecting the weak events, while stacking of raw waveforms had the best spatial resolution and retrieved magnitude. Therefore, a twophase backprojection imaging process is suggested (Beskardes et al., 2018): backprojection of kurtosis waveforms for event detection and location within large data volumes, followed by backprojection of raw seismograms of the previously detected events to improve the location accuracy and obtain more reliable magnitudes. A sketch of the principle for PWS methods is shown in Figure 6.
Generally, PWS can be written in a simplified form as creating an image function for all possible origin times t0:
11 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Reviews of Geophysics
10.1029/2019RG000667
Figure 6. The simplified sketch of migration and stacking process for partial waveform stacking. (a) Input waveform generation; (b) subsurface grid discretization; (c) source location picking after migration and stacking of the input waveforms. The raw data are processed to generate input waveforms, which are then stacked as the imaging value of each potential grid points; lastly, the source location is determined by picking the grid point(s) with the maximum imaging value. The illustration presents the general form of the theory and stacking process, neglecting some details (e.g., the origin time issue) and probable difference and variations for different methods. Hot colors of the imaging profiles indicate the increased energy due to constructive interference.
Sðx; t0Þ ¼ ∑ CFðtÞ·MðtÞ;
(1)
N
8 qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
CFðtÞ
¼
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>< >>>>>>>>>>>>>>>>>>>>>>>>>>>>>:
u2 þ H2ðuÞ; envelope
 2
∑ ∑u
nt N
N∑ ∑
u2
;
nt N
semblance
nl∑u
ns
ns∑
u
;
STA=LTA
nl
1 ∑ðuuÞ4

nt 1
nt
2; ðuuÞ2
nt nt
ÂÃ
∑½uiuiŠ ujuj
nt
ðntσiσj
;
kurtosis 2D multichannel
coherency
(2)
LI ET AL.
12 of 47
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
where x is the source coordinates; t0 is the origin time; t is the time sample in the waveform data u; N is the number of receivers; CF(t) is the characteristic function; M is the migration or imaging operators, for example, δ[t (t0+ti,x)] and δ[t (ti,x tj,x)] are imaging operators for Kirchhofftype and interferometrictype stacking methods; δ is the Dirac delta function; and ti,x is the theoretical traveltime from receiver i to source x. Please note that equation (1) is a generalized form of PWS method. Different combinations of CFs and migration operators lead to different specific formulae. The above stackingbased location problem can be regarded as a global optimization problem, in which source coordinates and origin time corresponding to the maximum imaging value max{S(x, t0)} are selected as the inverted source parameters. There are five representative CFs in equation (2), that is, envelope, semblance, STA/LTA, kurtosis, and multichannel coherency, where u is the input waveform, and it is not necessarily the raw waveform, H indicates the Hilbert transform function, nt is the length of selected time window for semblance and multichannel coherency, ns and nl are lengths of short and long time windows for STA/LTA, u is the mean value of u for nt samples, i,j are the station index, and σ is the standard deviation of the corresponding input waveforms. As the formulae denote, the envelope is computed as the amplitude after Hilbert transform, semblance represents ratio of signal energy to total energy for the signal, STA/LTA is the shortterm average to longterm average ratio, kurtosis is defined as the fourth normalized central moment higherorder statistics, and the two dimensional (2D) multichannel coherency is based on the Pearson correlation coefficients of the pairwise groups of records.
Currently, utilizing combination of direct P and S phases for seismic source location has drawn attention, due to the potential of introducing more constraints (e.g., Gharti et al., 2010; Grigoli et al., 2013; Pesicek et al., 2014). However, major challenges remain in building and calibrating S wave velocity model suitable for waveform stacking. Furthermore, the interferences from multiple phases might overweight their additional constraints when assumed multiple phases are submerged in noisy data, or inaccurate velocity models are used (e.g., L. Li et al., 2018). In theory, besides the primary phases, scattering or coda waves (including reflected and converted waves) can also be incorporated into the stacking process to enhance the coherence of the source energy. However, the traveltimes of these secondary phases depend even more strongly on the velocity model than the primary phases, and using such phases may introduce additional bias.
3.2. Time Reverse Imaging
TRI is based on the time reversal (TR) theory and technique, which was first established in medical field and nondestructive evaluation (e.g., Fink, 1999). TR has robust stability in focusing seismic source energy in both space and time. Currently, it is evolving into a standard approach for seismic source location and has been demonstrated to be a reliable location method in various fields, such as ocean acoustics and seismology (e.g., Artman et al., 2010; Douma et al., 2015; Gajewski & Tessmer, 2005; Kuperman et al., 1998). In TRI, the recorded full waveforms are extrapolated backward in time to focus the source energy, and the theoretical foundation is the reciprocity of the elastic wave equation. In principle, this approach shares the same essence with reverse time migration (RTM) technique in exploration seismology, which is based on the principle of timereversal invariance of wavefields, enabling optimal focusing of the source energy if a good velocity model is used (e.g., Baysal et al., 1983; Kees Wapenaar et al., 2010; H.W. Zhou et al., 2018). The inputs of RTM include both forward and reversed wavefield, and the goal is to image the structures, while TRI only backpropagates the recorded waveforms to focus the source energy. The general workflow of TRI consists of three steps (Figure 7): The recorded seismograms at all receivers are first reversed in time; then numerical modeling (e.g., finite difference) is conducted with timereversed seismograms as source wavelets; finally, similar to RTM, wavefield backpropagation is followed by an imaging condition or a focusing criterion (equations (3)(5)), which allows to determine source origin time and location (Figure 7).
This category of location method has been referred to as time reverse modeling (Gajewski & Tessmer, 2005; Steiner et al., 2008; T. Zhu, 2014), TR/reverse imaging (Artman et al., 2010; Kawakatsu & Montagner, 2008; Larmat et al., 2008; Sun et al., 2015), TR extrapolation (Z. Li & van der Baan, 2016), RTM (Nakata & Beroza, 2016), and reverse time imaging (Zou et al., 2014). Here, we use the term “time reverse imaging” for simplicity. During the past decade, many studies have shown the feasibility of TRI in locating seismic sources at laboratory and regional scales (e.g., Larmat et al., 2008; Steiner et al., 2008; Saenger et al., 2011; see
13 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Reviews of Geophysics
10.1029/2019RG000667
Figure 7. Illustration of time reverse imaging. The recorded seismograms are reversed in time and then backpropagated through the model, focused at the source location using an imaging condition. Except the velocity model, no other prior information (e.g., identification of events and the origin time) is necessary when locating events with time reverse imaging. The righttop subfigure illustrates the underlying Huygens principle behind time reverse imaging and why it performs better for reasonably dense and regular networks. The circles indicate the progression of the timereversed wavefields reinjected at the respective receiver locations (marked by yellow circles). The superposition of these elementary circular waves then represents the timereversed source wavefield.
LI ET AL.
Table 1 for more references). The application of TRI for seismic source location was first proposed by McMechan (1982), in which it was used to locate earthquakes. Gajewski and Tessmer (2005) further demonstrated its feasibility in locating lowSNR events with 2D and threedimensional (3D) synthetic acoustic examples and revealed its great potential for very weak events with large passive seismic monitoring networks. Steiner et al. (2008) applied TRI to synthetic and real lowfrequency microtremors originated from a small underground area in the oil reservoirs. They showed that the approach was feasible for microtremor location and could also provide a possible technique for reservoir location, though the characteristics of source radiation significantly affect location accuracy. T. Zhu (2014) improved the quality of TRI by compensating for attenuation and dispersion effects in attenuating media. Z. Li and van der Baan (2016) proposed a representationtheorembased reverse time extrapolation scheme for microseismic location, which can combine the threecomponent particle velocities (displacements) and the pressure wavefield and eliminate the ghost event in downhole monitoring. Werner and Saenger (2018) investigated the significance of station distributions, complex velocity models, and SNRs with respect to location accuracy with 3D numerical examples. They obtained several prerequisites for applying TRI reliably at field scale: Interstation distance should be smaller than the source depth, and the aperture of the array should be at least twice the source depth; an increased number of stations could alleviate the influence of noise; additional scattering of wavefields resulting from complex velocity model could improve location accuracy.
Theoretically, the imaging conditions proposed for RTM can be adapted to TRI in seismic source location. Among the various imaging conditions for RTM, the crosscorrelation imaging condition is a commonly used implementation and enables reliable amplitude reconstruction and good imaging resolution (e.g., Chattopadhyay & McMechan, 2008). Artman et al. (2010) introduced physically meaningful elastic time reverse imaging conditions for passive seismic data. Sun et al. (2015) extended the traditional TRI to a hybrid multiplicative TR imaging algorithm, which is based on multiplication between wavefields from receiver groups (equation (5)) and applies a causal integration over time to obtain the source evolution process of multiple asynchronous sources. Both synthetic and field data examples revealed the feasibility and superiority of the method in imaging the spatial and temporal distribution of microseismic events (T. Zhu et al., 2019). Nakata and Beroza (2016) proposed the geometricmean imaging condition by multiplying all
14 of 47
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
receiver wavefields at each space and time sample (equation (4)) and then take a summation over the time axis. Compared with conventional arithmeticmean imaging condition (equation (3)), this imaging condition could reduce the computational cost and increase the spatial resolution of the source location image. These imaging conditions for TRI in source location are shown in equations (3)(5).
N
Iðx; tÞ ¼ ∑ Riðx; tÞ;
(3)
i¼1
N
Iðx; tÞ ¼ ∏ Riðx; tÞ;
(4)
i¼1
g nj
Iðx; tÞ ¼ ∏ ∑ Rðj×njþkðx; tÞ;
(5)
j¼1 k¼1
where x is the source coordinates and t is the time sample, N is the number of receivers, Ri is the back
propagated
wavefields
g
at
receiver
i,
g
is
the
total
number
of
groups,
and
nj
is
the
number
of
receivers
in
each
group j , that is, ∑ nj ¼ N. Similar to equation (1), the source coordinates and time corresponding to the
j¼1
maximum imaging value max{I(x,t)} are determined as the inverted source parameters. Equation (3)
denotes summation over all the receiver wavefields, while equation (4) indicates multiplication of these wavefields. On one hand, the summation will result in nonzero values across the wave propagation paths, and multiplication will lead to nonzeros only at the focused sources (Sun et al., 2015; Sun et al., 2016). The imaging resolution corresponding to multiplication is supposed to be higher than that of summation. On the other hand, the summation operation in equation (3) is implicit, since the entire receiver wavefields can be backpropagated at once. However, the multiplication has to be conducted explicitly, which increases the computational cost. In order to leverage the efficiency and resolution of TRI, a hybrid multiplicative TRI was proposed by Sun et al. (2015), as shown in equation (5). The summation and multiplicationbased imaging conditions can also be termed as arithmeticmean and geometricmean imaging condition, respectively (Nakata & Beroza, 2016). If the entire receiver wavefields are set as one single group, the hybrid multiplicative imaging condition reduces to the autocorrelation imaging condition as in Artman et al. (2010). The performance of the imaging conditions for different situations (e.g., complex velocity models, irregular array geometry, complex source mechanisms, and earthquake clusters) still needs further investigation.
This methodology is extremely sensitive to the correctness of the velocity model, as focusing and defocusing of the energy are required not only for direct waves but also for all multiply scattered and refracted waves. This makes this methodology suitable only for a very simple medium or more likely to extremely well known velocity model. This poses significant challenges for microseismicity where the wavelengths of the propagating waves range from a few meters (downhole) to tens (surface) of meters and hence velocity models must be very well constructed and is also often problematic for natural seismicity at regional scale where only very general knowledge of the underground structure is available. This makes this methodology more popular and applicable in theoretical studies.
3.3. Wavefront Tomography
Seismic source location methods can be divided into three broad categories based on the role of the velocity model (Grechka & Heigel, 2017): (i) Geigertype inversion with a priori fixed velocities, (ii) tomographic techniques that iteratively update the velocity structure in the location process, and (iii) relative location methods that can position events relative to each other and relative to the selected velocity distribution. Due to the lack of sufficient constraints and related significant computational demands, the challenge of locating and characterizing an unknown source in a largely unknown medium is commonly divided into isolated subproblems. While this separate treatment of source and medium properties led to significant success and can still be considered common practice, there also exist techniques that aim at a joint estimation of the source location and the velocity distribution, thereby accounting for the coupled nature of the full inverse problem (cf. Grechka & Heigel, 2017). A recent attempt at confronting the practical challenges of the underdetermined joint inversion is based on an analysis of the local wavefield coherence (B. Schwarz et al., 2016; Diekmann et al., 2019). The method builds on the notion of wavefronts, which can be
15 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Reviews of Geophysics
10.1029/2019RG000667
Figure 8. Sketch of the wavefrontbased workflow for the joint inversion of velocity structure and source locations (Diekmann et al., 2019). Coherence analysis is
commonly based on the maximization of the semblance norm and results in estimates of reference onset time (indicated by a red dot), the tilt, and the curvature of seismic event as it is recorded at the surface. In addition to an improved signaltonoise ratio, new (imaginary) station responses can be generated. The estimated wavefront attributes then allow for multiple estimates of the source location and a joint inversion of the longwavelength velocity structure. Location estimates and the waveformbased uncertainty are indicated by black dots and a black circle, respectively.
LI ET AL.
reliably tracked, when a sufficient number of stations in a spatial vicinity are considered. There are multiple analytical functions that relate the geometrical attributes of wavefronts to traveltimes (e.g., Höcht et al., 1999; Jäger et al., 2001; Benjamin Schwarz & Gajewski, 2017). The most natural extension of the planar wavefront assumption that is commonly employed in beamforming is the secondorder traveltime expansion around a reference station,
t≈t0
þ
pðxx0Þ
þ
1 2
ðxx0ÞTHðxx0Þ;
(6)
where p is the horizontal slowness, x x0 is the lateral distance vector to the reference location, and t0
denotes the reference onset time. In addition to the firstorder term, the front curvatures of the locally coherent wavefield. For every regularized
2 × 2 matrix central data
pHoienntcÀoxd0e;st0thÁ,ewwitahvien-
a predefined local aperture, the maximization of a coherence measure like the semblance norm (cf. equa-
tion (2)) leads to a maximally coherent stack and local estimates of the coefficient vector p and the coefficient
matrix H (Jäger et al., 2001; Neidell & Taner, 1971). Equation (6) acts as a multidimensional fitting operator,
along which a coherence measure is estimated. The wavefront attributes contained in p and H represent
local approximations of the wavefield and can be traced back to the source. As is illustrated in Figure 8,
for sufficiently dense arrays, this approximation can be used not only to enhance waveforms through coher-
ent summation but also to estimate new (imaginary) station responses, as they would have been recorded at
an intermediate location. For a fixed central location and when the array aperture becomes small in compar-
ison with its distance to the source, the process reduces to conventional beamforming and results in esti-
mates of the directionality of incoming seismic radiation. In the general case of larger apertures, it
additionally includes the secondorder moveout term, which contains information on the emerging wave-
front's curvature. Crucially, in contrast to Kirchhoff migration, event and phaseconsistent traveltimes need
not be computed using a predefined velocity model but rather are automatically extracted from the data. The
fact that coherence analysis can be fully driven by collective wavefield patterns in turn makes it very suitable
for velocity inversion.
Given the nearsurface velocity, the estimated slope and curvature values can be used to characterize the actual geometry of the wavefront, as it arrives at the surface array. This is the underlying idea of WT (Bauer et al., 2017; Duveneck, 2004), which bears close resemblance to eikonal tomography (F.C. Lin et al., 2009) but is applicable for body waves and unknown source locations. Just like conventional traveltime tomography, it is based on a least squares misfit of picked and modeled data, where, however, in addition to traveltime information, the inversion also takes slopes and curvatures into account. As these represent collective properties of the wavefield, they directly relate and stabilize adjacent measurements. With the nearsurface constant velocity backprojection as a starting guess, trial source locations are iteratively updated with the velocity structure. In coherence analysis attribute fields describe the whole
16 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Reviews of Geophysics
10.1029/2019RG000667
Figure 9. Sketch of full waveform inversion for source parameters, for example, source location, origin time, source mechanism, and velocity model. In contrast to other waveformbased location methods, FWI relies on extensive forward modeling starting with a reasonable guess of the longwavelength velocity structure and source characteristics to arrive at full synthetic waveforms that are repeatedly compared with the observations.
recorded event locally, and WT typically delivers multiple location estimates for a single event (illustrated in Figure 8). This redundancy can in turn be used to additionally constrain the inversion by linking phase consistent contributions (Bauer et al., 2019; Ross et al., 2019). Due to the intrinsic possibility of data enhancement and regularization and because a robust estimation of the velocity structure becomes feasible, the wavefront picture naturally complements all other waveformbased seismic location methods, including FWI. Improved locations can usually be obtained by combining the estimated velocity model and reversetime imaging for the enhanced data set (Diekmann et al., 2019). The method appears to be efficient and accurate as long as the data coverage is dense and the underlying velocity structure captured with a single local aperture is only moderately complex. Similar to PWS, kinematic summation trajectories (here wavefronts) are combined with the utilization of waveform coherence, which makes this approach a hybrid method. Also, in analogy to PWS, the concept of diffraction is naturally honored (B. Schwarz, 2019), and the same workflow has been successfully applied to controlledsource backscattering (Bauer et al., 2017, 2019). Owing to different instrumentation and the scalerelated differences in wavefield complexity, the arraybased detection of directional information has dominated in earthquake studies (Vespagrams, slownessback azimuth analysis; e.g., Rost & Thomas, 2002), whereas in exploration seismology, wavefront tilt is often neglected, and curvatures have traditionally been at the center of attention (controlledsource velocity spectra; Taner & Koehler, 1969). In the framework of WT, where slopes and curvatures are combined, both approaches are unified. With the continuing rise of dense nodal deployments and the emergence of quasicontinuous strain measurements with groundcoupled fiberoptic cables, the utilization of collective properties of the wavefield is possible on very short length scales and can be expected to continue and intensify in the foreseeable future.
LI ET AL.
3.4. Full Waveform Inversion
FWI is a tool for geophysical imaging based on the waveform matching technique through minimizing the residuals between the synthetic waveforms and the observed waveforms (e.g., Virieux & Operto, 2009; see also Figure 9). Aside from many impressive applications in earthquake seismology, FWI has led to first highly resolved regional and global images of the Earth's mantle (e.g., A. Fichtner et al., 2009; Bozdağ et al., 2016). The different methods of FWI differ in the way of measuring difference between observed and synthetic waveforms (misfit) and inversion strategy of minimizing this misfit, as well as the variety of inversion parameters, for example, source location, moment tensor, and velocities. The misfit function is also known as objective function, although there is little objective about particular choice of the misfit functions. Enhanced computer resources and a method that enables efficient computation of the gradient of the misfit function, for example, the adjointstate method (Fichtner et al., 2006; Plessix, 2006), are the two crucial factors that make FWI applicable for field data sets. Several existing misfit functions of FWI for source location (and velocity inversion) are
17 of 47
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
FðmÞ
¼
1 2
N
i¼1
usiynðmÞuoi bs2;
(7)
FðmÞ
¼
1 2
N
i¼1
usiynðmÞuoi bs2
þ
ckmk1;
(8)
Fðm; vÞ
¼
N
usi yn ðm;
vÞ*uorebfsuoi bs*usryefnðm;
2;
(9)
i¼1
where m is the source model vector, including the spatial and temporal source components, that is, s(x,y,z,t) or f(x,y,z)w(t); usiyn and uoi bs are synthetic and observed data at receiver i; and N is the number of receivers. c is a positive scalar coefficient that controls the weight of the sparse model s(x,y,z,t). v is the velocity model, * denotes the convolution operator, and uref is the data of a reference trace. Equations (7) and (8) are commonly used misfit functions (e.g., Kaderli et al., 2018; Michel & Tsvankin, 2014, 2017; Shekar & Sethi, 2019), and equation (9) is recently proposed by Wang and Alkhalifah, (2018) and the unknown source origin time is decoupled by convolving reference traces with the observed and modeled traces. When FWI is applied to source location problem, the source location (and maybe together with velocity model) can be inverted and imaged in one step using above misfit functions.
In earthquake seismology, FWI has already been used to solve for the event location and source mechanism since more than two decades (e.g., Y. Wu & McMechan, 1996; RamosMartinez & McMechan, 2001; Sokos & Zahradnik, 2008). Y. Wu and McMechan (1996) used FWI to determine the spatial coordinates and origin time associated with a synthetic doublecouple source in 2D scenario, assuming the velocity model and the source time function are well known. The principal advantage of FWI application to microseismicity is in the short duration of the microseismic event and therefore no need to invert the source time function, which can be approximated by a delta function. Recent studies are trying to include the velocity model in the inversion. In the context of passive seismic monitoring, O Michel and Tsvankin (2014, 2017) utilized FWI to invert for the source parameters (including source location, origin time, and moment tensor), and velocity model (see Figure 9) in a sequential fashion, in which the adjointstate method is adopted to compute the gradient. Synthetic examples based on 2D layered VTI (transversely isotropic with a vertical symmetry axis) media have illustrated the feasibility of the sequential inversion for source parameters and VTI model. Kaderli et al. (2015) tested an acoustic FWI method to invert for the spatialtemporal source characteristics with a homogeneous velocity model, and the spatial and temporal components of the source are updated alternately. Sun et al. (2016) combined the FWI technique with timereversal imaging to obtain a least squares timereversal imaging method, which could provide accurate source location and source function after a sufficient number of linear iterations. Then these inverted source parameters could be further incorporated into the FWI for velocity model updates as shown in simple 2D acoustic model. Huang et al. (2017a, 2017b) developed an acousticwaveequationbased FWI method to locate microseismic events. They proposed to evaluate Frechet derivatives with respect to location parameters and apply the truncated Gauss Newton method to accelerate the inversion process. Synthetic and field microseismic data from above studies have been used to show the practical feasibility of the proposed method. Wang and Alkhalifah (2018) developed a source function independent FWI method to invert for the microseismic source location, source function, and the velocity model simultaneously. The effect of the unknown source origin time is mitigated by convolving reference traces with the observed and modeled traces (equation (9)). Hence, this methodology is extremely sensitive to selecting reference trace. Shekar and Sethi (2019) formulated an FWI scheme for determining the spatiotemporal source function of microseismic events, and an additional regularization term in the misfit function (8) is used by exploring the sparsity inherent in the microseismic location problem. Synthetic examples based on 2D SEG/EAGE overthrust model have illustrated the possibility of the FWI applications (Shekar & Sethi, 2019). Most recent studies of FWI for source location and velocity are limited to 2D synthetic acoustic scenarios, which are still far from field application, and therefore, it is hard to judge if they will be practical because microseismic signals are dominated by S wave energy (in the vicinity of the source) and this is not captured in the acoustic models. Furthermore, radiation complexity is simplified in most of these studies by isotropic (explosive) sources. While using FWI for a fixed velocity medium is already quite successful in source location, further studies on how to improve the performance (both accuracy and efficiency) of FWI including the velocity model are still needed.
18 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
3.5. Comparison and Relation of Different Methods
In a general sense, waveformbased methods locate the seismic source through backpropagation of source energy: TRI uses the reciprocity principle explicitly through backpropagation of seismic wavefields, while other methods make use of it by means of generating waveforms and/or the traveltime tables. The backprojection process of waveformbased methods can also be regarded as a generalized type of beamforming with an image function. Although the idea of waveform stacking is quite the same, the PWS methods distinguish themselves from the beamforming method in at least two major aspects (e.g., Rost & Thomas, 2002; Liao et al., 2012; Poiata et al., 2016; L. Li et al., 2018):
1. First, the nature of problems solved by beamforming is determination of the direction from which the source energy arrives, instead of determining the spatial and temporal distribution of seismic sources.
2. Second, PWS methods stack only waveforms of selected phases (primary phases in most cases) within specific time intervals, while beamforming method sums up all observed phases through a summation over time. More specifically, the general beamforming result is equal to the sum of autocorrelations and approximately delayed cross correlations of the waveforms, and the crosscorrelation term dominates the beamforming value.
An advantage for waveformbased methods is the ability to suppress noise. However, varying degrees of redundancy have been exploited for these methods. The classical diffraction stacking method using single channelbased CFs (e.g., envelope and STA/LTA) does not consider the correlation among the receiver channels, while multichannelbased waveforms (e.g., semblance) and the crosscorrelation stacking operator have explored mutual channel coherency. For TRI and FWI, the entire waveforms including multiples and scattered waves in addition to the direct waves are utilized; thus, in theory more accurate source locations and higherresolution velocity models can be achieved because the scattered waves add illumination points for the source. However, these methods require much more accurate velocity model, which is often not available. Therefore, various imaging functions and conditions are adopted to overcome poor focusing due to inaccurate model. The principal difference between aforementioned methods also lies in the distinction between raybased modeling and wave modeling techniques. The former corresponds to traveltime calculation for PWS, and latter corresponds the full waveform modeling of TRI and FWI. In general, the PWS is more robust because it requires the velocity model to be correct on average between the source and receivers along the whole ray, while the TRI and FWI require the inversion model to be accurate not only on average but to reproduce all reflectors and refractions, a rare achievement. Furthermore, joint velocity inversion and location may eventually become generally applicable, although it is extremely hard to identify the well constrained parts of the velocity model and avoid overfitting data and projecting model errors into source locations.
As discussed above, the critical issue for all methods is the velocity model, since most of these methods need accurate velocity models to locate seismic sources correctly. As explained above, the PWS methods need the velocity model to compute the theoretical traveltime table; this means that a smooth velocity model that predicts accurately average velocity along the path between the source and receiver is required. This is a great advantage of PWS methods. The TRI and FWI also require velocity model to simulate full seismic waveforms that need to realistically fit the observed data. Given the highfrequency nature of the microseismic events (see Figure 1), this may require velocity model reproducing synthetic seismograms at hundreds to several thousands of Hz in downhole monitoring and tens of Hz in surface monitoring. The FWI of active seismic data has partial success if lowfrequency waveforms can be inverted initially to fit long wavelengths of the velocity model. Unfortunately, none of current FWI or TRI method shows this feature with microseismic data sets. The joint inversion of velocity and location is an alternative solution. The wavefrontbased method provides an alternative for simultaneous location and velocity inversion through WT of long wavelength parts of the model. A better understanding and estimation of the role of the velocity model in waveform based location are still a challenge and perspective for future study, which will be discussed in detail in section 5.1. Another common shortcoming of TRI and FWI methods is the relatively expensive computational efforts, due to the storage and transmission of large volume of waveforms during the waveform backprojection process. Moreover, current TRI usually does not enable the flexible implementation of CFs as in PWS methods described above, and the implementation of the TRI generally requires the data to be recorded with a denser and more regular array, which may not be feasible for all configurations. Sparse monitoring
19 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Reviews of Geophysics
10.1029/2019RG000667
Table 2 Advantages and Disadvantages of Discussed WaveformBased Location Methods
Method
Advantages
Disadvantages
PWS TRI WT FWI
Flexible characteristic functions
Joint source mechanism inversion Requires only a smooth velocity model Full waveform exploitation
Either more accurate locations or smaller acquisition system
Local wavefield coherence exploitation
Data enhancement and regularization Full waveform exploitation Either more accurate locations
or smaller acquisition system
Requires a calibrated model to stack the selected part of the waveform
Uses only part of the waveform requires dense array
Need accurate velocity models Need good spatial coverage Computationally demanding Inability to handle complex source mechanisms Computationally demanding High nonlinearity Lack of source mechanism determination Computationally expensive Very high nonlinearity
arrays or any limited spatial coverage will lead to insufficient spatial sampling of the waveformbased location method. For TRI, FWI, and WT, the high nonlinearity of the inversion process associated with multiple source parameters (e.g., source coordinates and source mechanisms) and the velocity model is a challenge to get stable and reliable location results.
Advantages and disadvantages of waveformbased location methods discussed above are summarized in Table 2. It is also worth noting that different methods could be appropriate for different types of data inputs, depending on several factors, for example, application scale (frequencies), geometry of monitoring array, and monitoring targets and objectives. Table 3 lists several open source codes of waveformbased location methods.
Waveformbased methods seldom rely on purely waveform information; instead, they also exploit traveltime information in an explicit or implicit way. For instance, WT combines kinematic information from wavefront attributes and dynamic information through waveform coherence, to obtain a source location image. In this sense, the method can be also regarded as a hybrid method. Another example of hybrid method is an explicit combination of an automatic waveform stacking and manual picking corrections (Bardainne, 2011). This method involves a twostep procedure: First, a microseismic source location is imaged using a fully automated migration technique, then the unreliable picking or detections are manually verified and
Table 3 Available Open Source Codes of WaveformBased Location Methods
Method
Scenario
References
Address
Continuous kurtosisbased migration for seismic event detection and location
A combination of waveform inversion and reversetime modeling method
Multiband array detection and location of seismic sources
Automated full waveform detection and location algorithms
Locating tremor using stacked products of correlations
Diffraction and crosscorrelation stacking with stochastic algorithms
Multichannel Coherency Migration location
3D isotropic acoustic 2D isotropic (visco)
elastic 3D isotropic elastic 3D isotropic elastic 3D isotropic elastic
3D isotropic elastic 3D isotropic elastic
Langet et al. (2014)
Köhn et al. (2016)
Poiata et al. (2016, 2018)
LopezComino et al. (2017)
K. L. Li, Sadeghisorkhani, et al. (2017); K. L. Li, Sgattoni, et al. (2017)
L. Li et al. (2018); L. Li, Tan, Xie, et al. (2019)
Shi, Angus, et al., 2019; Shi, Nowacki, et al., 2019
https://github.com/amaggi/waveloc https://github.com/danielkoehn/DENISEBlackEdition https://github.com/BackTrackBB/backtrackbb/ https://gitext.gfzpotsdam.de/heimann/lassie https://doi.org/10.13140/RG.2.2.36158.08005
https://github.com/leileely/microseismic_stochastic https://github.com/speedshi/seisloc
Note. Data are accessed in December 2019.
LI ET AL.
20 of 47
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
corrected, and a refined location is obtained by combining the waveformbased stacking function and the conventional traveltimebased function. The method delivers more reliable results but also increase the processing time; thus, it provides a tradeoff between location efficiency and reliability.
In addition to the abovementioned methods, there are several variants of empirical waveformbased seismic detection and/or location methods based on correlation of waveforms, such as the template matching method and matched field processing. The main difference between them and abovediscussed techniques is the use of empirical waveforms, instead of a complex velocity model. Though these methods also exploit the waveform information in seismic detection and/or location process, they are slightly different from waveformbased methods discussed here and cannot be strictly classified into one of the above four subcategories.
Template matching method, also known as matched filtering algorithm, is an efficient approach for detecting weak events from continuous waveforms using cross correlation against a reference, the template event (Eisner et al., 2008; Gibbons & Ringdal, 2006; Peng & Zhao, 2009; Shelly et al., 2007). The method detects weak events with waveforms similar to the template event assuming that those occur at similar location and have similar mechanism as the template event. This fully automatic approach has been extensively applied to the detection of regular and lowfrequency earthquakes (e.g., Beaucé et al., 2019; Frank & Shapiro, 2014; Janská & Eisner, 2012; Ross et al., 2019; Shelly et al., 2007), foreshock and aftershock activity (e.g. Peng & Zhao, 2009; Sugan et al., 2014; M. Zhang & Wen, 2015), and microseismic events induced by hydraulic fracturing recorded by small networks (e.g., Caffagni et al., 2016; Eisner et al., 2008). The method can approach PWS by selecting template waveform only around P or S wave arrival and using corresponding velocity models for relative stacking between template (master) and detected events (e.g., Cieplicki et al., 2012) or using the FWI by selecting much larger template interval (e.g., the whole P to S wave interval). This choice makes the methodology not being able to strictly fit into PWS discussed above.
For large data sets collected by dense or regional networks comprising long time series (e.g., 110 years), highperformance computing exploiting multiple central processing units (CPU; Vuan et al., 2018), and/or graphics processing units (GPU; Beaucé et al., 2018) to achieve feasible computing times. Alternatively, more efficient frequency domain operation (Senobari et al., 2019) may be applied. Template matching is commonly used to detect small events or lowSNR signals (SNR < 1) requiring a prior knowledge of template events in the target area without relatively locating them. If the template matching is used also for relative location, the stacking process involved in template matching is applied to crosscorrelation function, and it is designed to produce high stack function as a detection function over all possible relative locations:
CCðx;
y;
z;
¼
1 N
N
i¼1
NCC½tτiΔtiðx;
y;
Š2;
(10)
where CC is stack of cross correlation and NCC is normalized cross correlation between template (master) waveform and detected event waveform. The (x,y,z) represents spatial coordinates where we search for the detected event. Δti is relative traveltime between the master and detected events. τ stands for delay between origin times. The above equation is normalized to be independent of the number of receivers (N) to have maximum value of 1. The advantage of this approach is low sensitivity to the velocity model in the overburden; the drawback is the need for strong master events.
There is a continuous ongoing improvement of the method, and several studies have extended it to a detectionandlocation approach by incorporating waveform stacking techniques (e.g., Withers et al., 1999; Cieplicki et al., 2012; M. Zhang & Wen, 2015; Arrowsmith et al., 2016; Caffagni et al., 2016; Beaucé et al., 2019). For example, Cieplicki et al. (2012) proposed to stack the crosscorrelation waveforms according to the traveltime difference between the template and potential small events. Song et al. (2014) and Caffagni et al. (2016) adopted the matched filtering algorithm for analyzing microseismic events recorded by downhole monitoring arrays. More recently, Beaucé et al. (2019) presented a detection and location approach that combines a PWS approach (Frank & Shapiro, 2014), for estimating the target seismic activity used as template events and the template matching for small event extraction, followed by a relative relocation step. Such composite approaches provide an efficient tool for generating detailed seismic catalogs.
21 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
Matched field processing is an arrayprocessing technique originated from ocean acoustics to track submarines or marine mammals (Baggeroer et al., 1993). It is an extension of plane wave beamforming that aims at locating coherent noise sources through correlation of phase information from the measured wavefield and a synthetic wavefield. The approach considers each grid node as a potential source position and yields phase matches as corresponding imaging values, which indicate the estimation of the source location. Recently, matched field processing has been tested with success for seismic source location in the field of environmental and exploration seismology, such as hydrothermal vents of geysers (Cros et al., 2011; Legaz et al., 2009), hydrocarbon reservoirs (Chmiel et al., 2019; Corciulo et al., 2012), and CO2 degassing zones (Umlauft & Korn, 2019). A similar approach based on timereverse principle and using strain Green's tensor as a matched filter was applied for locating nonimpulsive seismic events along the midoceanic ridge system west of Mexico (Solano et al., 2017). Compared to above waveformbased methods, which mainly work in the time domain, matched field processing is a frequency domain method that takes into account the complex phase, instead of the amplitude of the waveforms. An important advantage of the method is the capability of detecting longduration tremorlike events in the frequency domain, which might not produce a constructive stack in the time domain.
Although very efficient for weak event location, both template matching and matched field processing have a practical limitation, namely, that the seismic sources (i.e., the master event and the target event) must be located close to each other regarding the distance between the master event and the receiver array. In addition, the source mechanisms must be similar to make the matched filter or field effective. Thus, their applicability is limited by the available set of template (master) events and therefore is limited to post processing (Gibbons & Ringdal, 2006; J. Wang et al., 2015). The big advantage of these methods is low sensitivity to the velocity model, which makes it very attractive compared to all abovementioned methods.
4. Applications
Recent studies and applications at different scales have already revealed the performance of waveform based methods and could provide a guidance on their applicability and limits. Next, we will elaborate on several representative examples of waveformbased location methods, ranging from laboratory scale to field and regional scales. So far, only PWS are applied as processing routines for field applications; other waveformbased methods, e.g., TRI and FWI, are still at their infant stage for field application, perhaps due to their high computational expenses, as well as many other issues listed in Table 3. Therefore, most examples discussed here are based on PWS, which has been successfully applied at multiple scales.
4.1. Laboratory and Small Scale 4.1.1. AE Monitoring Experiments AE monitoring under controlled experimental conditions has been widely used to investigate rock fracturing (e.g., Lockner et al., 1991; Sellers et al., 2003; Johnson et al., 2008; Benson et al., 2008; Goebel et al., 2013; Bolton et al., 2019). Traditional AE technology is limited by data acquisition systems and computer storage capabilities, and it is mainly based on statistical parameters (e.g., event count, location, type of failure, and energy) that characterize the activity of AE (e.g., Ishida et al., 2017; Mair et al., 2007). These conventional approaches are called as parametric AE analysis (Grosse & Ohtsu, 2008; Hardy, 2003). A conspicuous shortcoming of parameter analysis is that AE parameters can only provide limited analysis of the fracturing mechanism. In order to improve the reliability of AE technology and further quantify the processing results, the researchers proposed another category of signalbased techniques by making use of the recorded AE waveforms (Grosse & Ohtsu, 2008; Ishida et al., 2017; Stanchits et al., 2015). Waveformbased AE location is a relatively new application of signalbased techniques. There are at least four main challenges that make reliable source characterization of AE events difficult, including lowSNR and high frequency content of AE data, high occurrence rate of AE events, and unreliable sensor calibration (Gu, 2016). Another specific challenge for waveformbased location and characterization of AE is extremely short interval of different phases and reflected waves from sample boundaries. For this reason, PWS methods are thought to be more preferential for AE studies in the foreseeable future.
Saenger et al. (2011) applied TRI to locate and characterize AE events with a numerical model of concrete block. The forward and back propagation of the wavefields is simulated with the rotated staggered finite
22 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
difference method, and the maximum value for each grid point throughout the entire time of modeling was stored to image the event location. Through the comparison between source TRI with limited knowledge of boundary values and full TRI where a complete set of boundary conditions, they demonstrated the feasibility of TRI in locating source areas and characteristics of AEs with a limited number of sensors and effective elastic properties.
Gu (2016) implemented the sourcescanning algorithm—a PWS method—along with two traveltimebased location methods to laboratory AE data. The laboratory loading experiment was conducted on a cylindrical Berea sandstone sample (diameter 36.43 mm and height 76.7 mm), and eight piezoelectric sensors were attached to the rock sample and connected to the data acquisition system to collect AE data. The AE data segment was divided into small windows, and then sourcescanning algorithm was applied to each window assuming only one event was contained in each window. The detection and location were achieved simultaneously by setting a threshold of the imaging value. This study showed that sourcescanning algorithm could provide the initial guess of the AE locations and delineate the fracture shape roughly. 4.1.2. SmallScale Hydraulic Fracturing Experiments Smallscale in situ hydraulic fracturing experiments bridge the gap between laboratory studies and field scale stimulation. Hydraulic fracturing generates massive microseismic or AE data set that should be quickly detected and located. A field test of hydraulic fracture monitoring in hard rock was conducted at the Äspö Hard Rock Laboratory (Sweden) during May and June 2015, and the hydraulic fracturing experiments were operated to compare fracturing growth and induced AE activity under different fluid injection schemes (LópezComino et al., 2017; Zang et al., 2017). Four boreholes were drilled at the experimental tunnel TASN, the orientation of the boreholes was subparallel to the orientation of minimum horizontal compressive stress, allowing hydraulic fractures to propagate perpendicular to the treatment borehole.
The in situ monitoring network consisted of 11 AE sensors and four accelerometers. In this experiment, the acquisition system with the piezoelectric sensors was improved to record continuous signals with 1MHz sampling rate. Homogeneous velocity models obtained from active ultrasonic transmission tests were used for an envelopebased diffraction stacking method to locate AE events (Figure 10). Continuous waveforms along with traveltimes of P and S wave were used in the waveform stacking process. Stacking of CFs (Figure 10b) yielded the preliminary location (Figure 10c) and a quantified detector function (Figure 10d).
LópezComino et al. (2017) presented a novel study on implementation of sufficiently robust and automatic algorithms to detect and locate massive microseismic or AE activity. More than 4,000 AE events were located using the STA/LTAbased waveform stacking and coherence analysis techniques. The absolute location of the largest 258 events (with a detector level threshold of 1,500, which is calculated with normalized amplitude envelope as the CF) delineated the fracture geometry. A low detector threshold can enhance the catalog completeness, but a higher detector threshold may better exclude false detections. Moreover, the location results were further refined by a master event stacking approach (F. Grigoli et al., 2016), which produced spatially clustered locations that outlined the main fracture plane. The Jackknife resampling of STA/LTA configuration was utilized to estimate the location uncertainty, and the median uncertainty was reduced from 1.43 m to 0.64 m through incorporating a master event of high accuracy and maximal coherence value. In summary, smallscale hydraulic fracturing experiments generate highrate, highfrequency AE events. The PWS techniques allow automatic AE events detection and location using continuous waveforms.
4.2. Field and Local Scale 4.2.1. MiningInduced Seismicity Seismicity associated with underground mining has been observed since early last century (e.g., Cook, 1976; S. J. Gibowicz, 1990) and is tightly related to the safety and productivity of mining. For instance, rockbursts are very often the major cause of fatalities in mines. Seismic monitoring can provide a powerful tool for the detection and evaluation of seismic events occurring around underground mines (e.g., Gibowicz & Kijko, 1994; T. Li et al., 2007; Ma et al., 2019). Two examples of waveformbased location methods applied to mininginduced seismic events are shown in this section. A temporary microseismic monitoring network (HAMNET) was deployed from June 2006 to July 2007 in order to record the coal mininginduced seismicity in the Ruhr region, Germany (Bischoff et al., 2010; Grigoli et al., 2013; WehlingBenatelli et al., 2013). The HAMNET network consisted of 15 threecomponent surface seismic stations, deployed in a region of about 3 km × 2 km (Figure 11a). The monitored longwall panel has a spatial extent of about 1 km × 0.3 km and is
23 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Reviews of Geophysics
10.1029/2019RG000667
Figure 10. Detection and location of a sample AE event using envelopebased diffraction stacking (from LópezComino et al. (2017), with permission from Springer Nature). (a) Waveforms of a sample AE event recorded by the 15 sensors. (b) Characteristic function (i.e., normalized amplitude envelopes) for each trace, which are used for stacking with traveltimes of P wave (red line) and S wave (green line). The markers indicate the (best fit) synthetic arrival time of the respective phases at each receiver. (c) Coherence map for the search region. Dark colors denote high coherence values. The location of the detected event is marked by a white star, and receivers are denoted by black triangles. (d) Detector function in a broader time window around the origin time of the detected event. The cutout time window used for the coherence map in (c) is shown in gray color. White stars indicate two detected events within the same processing time window, exceeding a detector level threshold of 500.
located at a depth of about 1.1 km. The exploitation of the coal seam was performed from August 2006 to April 2007. More than 7,000 events have been detected and located using P wave onsets and a homogenous velocity model.
Accompanied by testing on 200 selected events from the above described data set, L. Li et al. (2018) provided a systematic analysis of correlationbased location methods and unified the formulae by describing cross correlation stacking with beamforming, and L. Li et al. (2019) compared the performance of two PWS methods (i.e., diffraction stacking and crosscorrelation stacking) with different stochastic optimization algorithms. The waveforms for all the events had already been truncated as separate files, and the raw waveforms were bandpass filtered (530 Hz). The STA/LTA traces, which were calculated with the energy of filtered waveform (equation (2)), were used as the input waveforms. Only P waves were used to locate these events at first and then compared the result with that of using multiple phases. Figure 10b shows the seismograms and STA/LTA traces of a sample field event, which has a ML magnitude of 0.8. The lateral imaging resolution for the two methods is of reasonable quality due to the proper lateral coverage of the surface monitoring array (Figures 11c and 11d). The unknown origin time is canceled out in the correlation process of crosscorrelation stacking, and the basic imaging pattern is hyperbolae (hyperboloids) intersections in 2D (3D) scenario. The decoupling of the origin time from the location process alleviates the origin time depth tradeoff, which is an inherent characteristic of diffraction stacking (e.g., Anikiev et al., 2014), but significantly lower the vertical imaging resolution.
Another example is waveformbased microseismic event location in the SchlemaAlberoda mining area (Germany). Some of the microseismic events that occurred between 2005 and 2012 in this uranium mine area were located through Kirchhofftype approach (PWS method) to investigate the longterm role of mininginduced seismicity (Hassani et al., 2018). A homogeneous and a 3D velocity model were compared for seismic migration. The resulting images were comparable and exhibited minor differences to the subsurface reflectors. However, the waveformbased location method was more sensitive to the lateral and vertical
LI ET AL.
24 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Reviews of Geophysics
10.1029/2019RG000667
Figure 11. Waveforms and location results of an example field event associated with mininginduced seismicity monitored by the HAMNET network. (a) The location of the field mininginduced events by traveltime inversion. White triangles are stations (the station HMES is noted as the 15th station HM15), and circles are events. Colors of the events indicate the origin time and the size of the circle scales with ML magnitude. (b) The vertical component (left) and STA/LTA traces (right) for a sample field event. (c) The imaging result for the sample event by diffraction stacking, the reference point (east, north) = (0, 0) corresponds to east = 411117 and north = 5721611 in the UTM coordinate system. (d) The imaging result for the sample event by crosscorrelation stacking.
variations in the velocity model than seismic migration algorithm. Waveformbased methods can provide highprecision results using only P waves when an accurate velocity model is available, but a homogeneous or average velocity model may also lead to unreliable locations. Comparison between located seismic sources and the 3D reflection seismic image of the area confirmed that some of reflectors can be directly correlated with the induced seismicity. 4.2.2. Hydraulic FracturingInduced Local Seismicity Hydraulic fracturing stimulation is a methodology used in the oil and gas industry to produce hydrocarbons from unconventional reservoirs or to create heat exchangers in geothermal exploration. Hydraulic fracturing induces microearthquakes, which are located and characterized to map the hydraulic fractures. The monitoring arrays can be installed either in the nearby monitoring wells or deployed at the (near)surface. In contrast to the limited number of receivers (usually smaller than 20) in the well, the surface arrays may consist of several hundreds of stations in different layouts such as circles, regular grids, patches, or starlike patterns. A high number of receivers and spatially dense recordings represent ideal conditions for waveformbased location methods.
In the context of unconventional gas reservoirs, processing of data acquired at large array configurations was initially tailored to mapping high energy radiation (Duncan, 2005). It was then refined to map microseismic events as proposed in academic studies, for example, detect high semblance over space, time, and stations
LI ET AL.
25 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Reviews of Geophysics
10.1029/2019RG000667
Figure 12. Joint source location and source mechanism inversion example for microseismic events induced by hydraulic fracturing (modified after Anikiev et al., 2014, with permission from Oxford University Press). (a) Threedimensional P wave velocity model, wells, and receiver geometry. (b) Example of data with one strong event. The times of the arriving P wave are highlighted by green lines. (c) Obtained stack function for the time interval of 3 min, result of STA/LTA detection applied on the obtained stack, and the semblance function. The example event is detected at 76.8s recording time. (d) The automatically taken P wave arrival amplitudes for the example event. The red and blue circles represent amplitudes with negative and positive polarity, respectively. The size of the circles is proportional to the amplitudes. (e) Locations of the detected events.
LI ET AL.
(Kiselevitch et al, 1991); evaluate brightness function created by stacking of absolute amplitudes in source scanning algorithm (e.g., Kao & Shan, 2004, 2007; Liao et al., 2012); and stack envelopes of waveforms along bodywave moveouts (Gharti et al., 2010) or using semblanceweighted stacking (W. Zhang & Zhang, 2013). The previously referenced methods account for polarity changes arising from the source mechanism's radiation pattern, while only providing location estimates. Later, many authors started to develop methods to determine not only event location but also source mechanism. Gharti et al. (2011) combined TRI with source mechanism estimation; Rodriguez et al. (2012) presented a method for the simultaneous determination of the source origin time, location, and the radiation pattern. Zhebel and Eisner (2015) developed a workflow that combines migrationtype location capabilities with the estimation of the source's full moment tensor. With further improvements, Anikiev et al. (2014) applied this method to a competitive data set acquired during hydraulic fracturing in shales in Oklahoma. Staněk et al. (2015) improved the detection by adding semblance computation of stacked amplitudes along the moveout as another criterion to traditional STA/LTA. It allowed them to eliminate false detections and detect more weak events (Figure 12), which were further interpreted with a bedding plane slip model (Staněk & Eisner, 2017).
26 of 47
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
4.3. Regional and Global Scale 4.3.1. Natural Events Seismic source location and characterization are a key challenge in the analysis of both natural and anthropogenic seismicity. Detailed seismic catalogs are backbone of many studies in seismology such as characterizing the crustal stress field and estimating the Earth structure. Their practical values are directly related to the seismic monitoring and earthquake forecasting. The development and application of waveformbased detection and location methods for the analysis of seismic activity at regional and global scales were pushed forward by densification of seismic instrument deployments as permanent networks (e.g., Figure 13d) or temporary observations (e.g., Figure 13a) as well as accessibility of computational resources. Although traveltimebased methods continue to be a standard procedure for earthquakes detection and location used by seismological centers, waveformbased methods exploiting the coherency of seismic wavefield characteristics are attracting more attention due to their capability to efficiently process large volumes of seismic data in a fully automatic way. Here we illustrate application of waveformbased crosscorrelation stacking to the detection and location of subduction earthquakes in central Chile (Poiata et al., 2016) and the low frequency earthquakes (Katsumata & Kamaya, 2003), shorttransient events hidden in tectonic tremor signal, in western Shikoku, Japan (Poiata et al., 2018), as well as volcanotectonic earthquakes at Uturuncu volcano (Shi, Nowacki, et al., 2019). Compared to diffraction stacking which directly utilizes waveforms of individual stations (e.g. on N traces), the crosscorrelation stacking exploiting multichannel coherency of different station pairs (e.g. N(N1)/2 station pairs), and can extract more information from continuous seismic recordings.
The selected applications cover what are currently the main areas of interest and concern in seismic source seismology extending to the problems of seismic activity monitoring and source characterization. Automatic detection and location of earthquakes during the foreshock/aftershock and swarm sequences or regular background activity in seismically active environments such as subduction zones or active faults provide important information about the state of stress field in the area. The ability to highlight the evolution of seismic activity in greater detail (down to low magnitudes) during, for example, foreshock sequences can potentially lead to the information about the initiation stage of the larger earthquakes. Lowfrequency earthquakes and tectonic tremor are a class of seismic slow earthquakes (K. Obara, 2002) that contain a unique source of information about the mechanical properties of subduction zones during the seismic cycle. Observed in the transition zones surrounding the regions of large (megathrust) earthquakes, slow earthquakes might represent a unique source of information for improving the seismic hazard and risk assessment (K. Obara & Kato, 2016). Thus, detailed monitoring and characterization of slow earthquakes activity represent an important topic in seismic source seismology. In the seismic environments associated to the volcanoes, volcanic earthquakes such as volcanotectonic events and longperiod events are generally associated to the movement of magma below the Earth surface. Volcanoseismic monitoring is an important tool for prediction of volcanic eruptions, as, in many cases, precursory increase of the seismic activity rate often precedes large eruptions. Currently, there are about 1,500 potentially active volcanoes worldwide, and more than 200 volcanoes are seismically monitored (Thompson, 2014). Volcano hazards pose great danger to daily life of nearby populated regions and can affect extended area far from a volcano. Therefore, rapid (automatic) and accurate detection and location of associated seismic events (e.g., volcanotectonic earthquakes and volcanic tremors) are essential for earlywarning and hazard assessment purposes as well as studying the mechanism and structure of volcano system.
With the increased volumes of data available from networks and arrays installed for natural seismic activity observations (in both tectonic and volcanic environments), waveformbased detection and location methods become an increasingly important tool for efficient and fully automated processing of the waveform data having the ability to provide detailed and uniform catalogs. The examples presented here provide an illustration supporting this.
An example application of this approach to the detection and location of subduction earthquakes in central Chile recorded by 22 selected regional seismic stations is shown in Figures 13a13c on a 3minlong velocity record containing two events. Only vertical components are considered under the assumption of P wave propagation in a homogeneous velocity model. Simultaneous detection and location are performed when the value of the image function, normalized by the number of station pairs, exceeds a threshold value of 0.7 (Figure 13c). In this implementation, the CFs (Figure 13b) correspond to the timefrequency kurtosis
27 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Reviews of Geophysics
10.1029/2019RG000667
LI ET AL.
Figure 13. Detection and location example for regular (ac) and lowfrequency (dg) earthquakes (modified after Poiata et al. (2016, 2018), with permission from Oxford University Press). (a) Map view of central Chile; triangles indicate seismic stations. Black open rectangle highlights the investigated area in (c). (b) Original (gray) and kurtosis CFs (black) traces for a M 5.2 earthquake and a preceding M 3.2 event. Green shaded rectangles indicate time windows corresponding to image functions in (c); red ticks show estimated P wave arrival times. (c) Crosscorrelation imaging function for the two time windows in (b). White triangles indicate seismic station. (d) Map view of western Shikoku; black squares show Hinet permanent seismic stations (Okada et al., 2004). Colored dots show location of lowfrequency earthquakes from automatic analysis of 9day continuous data set with crosscorrelation stacking method of Poiata et al. (2016). (e) Example of a 1hr tectonic tremor record (NS component). (f) Waveforms (gray) and kurtosis CF (black) of the portion of seismogram marked by green vertical lines in (e). Red ticks mark estimated S waves arrivals. (g) One horizontal and two vertical sections through the final estimated crosscorrelation imaging function corresponding to the lowfrequency earthquake in (f). White triangles indicate stations used for the location.
28 of 47
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
estimated by the multiband filter approach (Poiata et al., 2016) providing an efficient enhancement of signals like small aftershocks hidden in the coda of large events. This increases the number of the detected and located small earthquakes. An extended application of the same PWS method (Poiata et al., 2016) to a larger data set in the subduction zone of central Chile is shown in the study of Ruiz et al. (2017), analyzing the precursory (foreshock) activity of the 2017 Valparaiso earthquake using the recordings of the temporary and permanent stations of the National Seismological Center of the University of Chile.
Since the discovery of tectonic tremor by K. Obara (2002), one standard technique for the detection and location of such seismic events is based on envelope crosscorrelation stacking (K. Obara, 2002; Kao & Shan, 2004; Wech & Creager, 2008; Ide, 2010). A typical signal of deep tectonic tremor recorded at the surface stations (Figure 13d) is depicted as waveforms energy arrivals that are coherent across the stations (Figure 13e). Thus, envelope CFs estimated in the predominant frequency band of tremor signals (210 Hz) are used in the PWS location schemes assuming S waves as the dominating phase. The tectonic tremor signals may contain mixed smallmagnitude impulsive events called lowfrequency earthquakes (Figure 13f; Katsumata & Kamaya, 2003) that are mixed within the tremor signal (Figures 13e and 13f). Detection and location of these events are highly challenging, since their waveforms are characterized by very low SNR. Poiata et al. (2018) demonstrated that the automatic PWS detection and location scheme based on timefrequency kurtosis and the waveformbased crosscorrelation stack allows for efficient automatic extraction of such events. Figures 13d13g show an example of the detection and location of lowfrequency earthquakes using this approach applied to the 9day continuous records of tremor activity recorded by Hinet (Okada et al., 2004) seismic stations in western Shikoku (Japan). Theoretical S wave traveltimes were calculated assuming a onedimensional velocity model of the area defined at a 3D grid covering the region of interest (Figures 13d and 13g). The maximum likelihood source location is estimated by stacked stationpair cross correlated kurtosis CFs projected according to the theoretical time delays (Figure 13f). Event detection is declared each time the normalized imaging function exceeds the threshold value of 0.65 (Figure 13f). The spatial location corresponds to the maximum of the determined image function, and the timing is calculated according to the provided location and automatically picked S wave arrivals at the stations (Figure 13f). This analysis provided a catalog of 2,212 lowfrequency earthquakes showing details of migration and clustering along the strike of the subduction slab. In comparison, the lowfrequency catalog, based on manual picking and conventional traveltimebased methods, provided by the Japan Meteorological Agency contains 341 events for the same period of time (see Poiata et al., 2018, for more details).
Although there have been many efforts to increase network coverage and station density (Dougherty et al., 2019), seismic monitoring array with poor coverage and sparse deployment is still very common especially in barren areas. Another example illustrating the performance of crosscorrelation stacking is locating volcanictectonic earthquakes recorded with sparse stations at Uturuncu volcano (Shi, Nowacki, et al., 2019). Figure 14a shows a regional elevation map at Uturuncu volcano, Bolivia. As shown, 15 three component seismometers have been temporarily deployed surrounding the inflating Uturuncu to monitor the volcanic activities (Jay et al., 2012). The largest intersection spacing between stations is about 40 km. Four different CFs, that is, envelope (Kao & Shan, 2007), STA/LTA (Drew et al., 2013), kurtosis (Langet et al., 2014), and Pearson correlation coefficients (Shi, Angus, et al., 2019), are compared using the same example of volcanictectonic event. Figure 14b shows the migration profiles of the four different migration methods; the catalog location of the example event is displayed as the red star in Figure 14a. The source energy of the envelope and STA/LTA migration is not well focused due to strong scattering in the recorded wavefield. The horizontal locations of this seismic event using the kurtosis and correlation methods are well consistent with the catalog location. However, the located event depths are deeper than the event depth in the catalog (1.72 and 1.92 km deeper, respectively). Nevertheless, compared to the horizontal location of the seismic event, the event depth is often not well constrained, especially for surface arrays. As discussed (e.g., in section 4.1.2), the tradeoff between event depth and origin time often makes event depth determination problematic and more difficult (e.g., Eisner et al., 2010). The crosscorrelation stacking method is then applied to locate triggered volcanotectonic events on 4 hr of continuous seismic data. The recorded waveform data at a station are shown in Figure 14c, in which many small magnitude events are very difficult and time consuming to pick manually. As the triggered earthquakes start immediately after the surface wave train of the Mw 8.8 Maule earthquake, many events occurred in a short time period with very close or overlapping waveforms. Figure 14d shows the variation of the stacked maximum coherency values at different
29 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Reviews of Geophysics
10.1029/2019RG000667
Figure 14. Detection and location example for volcanic earthquakes (modified after Shi, Nowacki, et al., 2019, with permission from Oxford University Press). (a) Location of the seismic stations and Uturuncu volcano (UTM zone: 19K). The stations are represented by gray reversed triangles. A local volcanotectonic earthquake in the catalog is highlighted by the red star. The color in the figure represents elevation relative to the sea level. The lower left part exhibits a regional map, in
which the red rectangle shows the research area. The lower right part exhibits the velocity model used in the event location, in which the red and blue lines show the P and S wave velocities. The white rectangle shows the imaging area for the 4hr continuous data. (b) Imaging profiles through the maximum migrated value for the volcanotectonic earthquake displayed in (a). The black stars show the corresponding seismic event in the catalog obtained by manual picking. The four columns show results of envelope, STA/LTA, kurtosis, and correlation stacking. The upper and lower rows show XZ and XY profiles, respectively. (c) The recorded Zcomponent waveforms at station UTCA. The instrument response has been removed and the waveforms are filtered using a bandpass filter of 4.221.6 Hz. (d) The maximum coherency value corresponds to origin time for the continuous data. The time interval is 0.08 s. The red points show the 322 verified seismic events. (e) The earthquake locations on the one horizontal and two vertical profiles. Black dots show the 114 event locations in the existing catalog. The colorcoded dots show the verified 322 event locations by crosscorrelation stacking. Colors of the dots indicate the origin times of the earthquakes.
LI ET AL.
30 of 47
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
trial origin times. There are many local peaks rising from the background noise in Figure 14d, which correspond to potential seismic events. In total, 560 seismic events are identified in Figure 14d by tagging the maximum value of each local peak (the number depends of the applied threshold criterion). Finally, 322 seismic events have been verified, which have clear phase arrivals. The verified seismic events are shown as red dots in Figure 14d. Figure 14e shows the locations of the 322 verified seismic events (colored dots) and the 114 current catalog events (black dots). The distribution of automatically located seismic events is consistent with that of the events in the catalog. There are two main earthquake clusters. One is located in the northern part of the study area and close to the volcano. This earthquake cluster occurred earlier (68 a.m.), and the events are mainly triggered by the surface waves (Love and Rayleigh waves) of the Maule earthquake. The other earthquake cluster occurred from 8 to 10 a.m. and is located in the southern part of the study area. The seismic events are mainly triggered by the surface wave overtones of the Maule earthquake. 4.3.2. Induced Events Many human activities can induce seismic activity, the magnitudes and patterns of which vary from continued microseismicity (e.g., mines, hydrocarbon, and geothermal fields), intense microseismicity after a certain event (e.g., hydraulic fracturing), moderatemagnitude seismic activity (e.g., hydrocarbon fields and mines), to catastrophic activity of large magnitude (e.g., reservoir induced seismicity; Suckale, 2009). Induced seismicity monitoring is important for both industrial operations and public safety. Monitoring of the associated microseismicity is essential for reservoir characterization, for example, fracture geometry delineation, fluid migration, and reservoir geomechanical analysis, while largermagnitude earthquakes are mainly exploited for potential geohazard management and mitigation (L. Li, Tan, Wood, et al., 2019). Several cases of local and regional earthquakes associated with hydraulic fracturing have been reported worldwide, such as in southcentral Oklahoma, the United States (Holland, 2013); the Blackpool Area, the United Kingdom (Clarke et al., 2014); the West Canadian Sedimentary Basin, Canada (Atkinson et al., 2016; Bao & Eaton, 2016); and the Sichuan Basin, China (Lei et al., 2017).
Caffagni et al. (2016) adopted a new Matched Filtering Algorithm for detecting and analyzing microseismic events recorded by downhole monitoring arrays. The method required a set of welllocated events that serve as templates to detect and analyze new events. The method involves stacking multichannel cross correlation of the continuous waveform data and waveforms of master events. Besides detection of weak events as conventional template matching techniques do, this new algorithm can also locate events by maximizing the amplitude of stacked envelope functions surrounding the template event. The moveout correction of the envelop functions is determined using a lookup table generated with 2D cylindrical coordinates, defined by radial distance and depth. This approach has been successfully applied in several field data sets associated with hydraulic fracturing in western Canada, such as the multistage treatment at the Garrington field in Alberta (Caffagni et al., 2016), regional earthquakes associated with hydraulic fracturing in the Duvernay shale play (Bao & Eaton, 2016), and the Tony Creek Dual Microseismic Experiment in the KaybobDuvernay region of Alberta (Eaton et al., 2018; Zhang et al., 2019). This waveformbased Matched Filtering Algorithm was demonstrated to produce a more complete catalog, which enables more robust interpretation of the induced microseismicity and potentially provides novel insights into corresponding dynamic rupture processes. A more generalized form of waveform stacking referring to template waveforms in 3D Cartesian coordinates has also been reported in F. Grigoli et al., (2016). Other variants of waveformbased relative location methods involve crosscorrelation waveform stacking based on differential traveltime between template and target events (e.g., M. Zhang & Wen, 2015; L. Li et al., 2016).
Verdon et al. (2017) introduced another example of a waveformbased method applied to induced (micro) seismicity associated with a multiwell, multistage hydraulic fracturing treatment in Oklahoma. The surface monitoring array consisted of 17 threecomponent broadband seismometers. The target depth for the stimulation was approximately 3.8 km. A total of 155 events were detected and located using a STA/LTAbased PWS method, including the 17 events with relatively higher SNR and larger magnitude (ranged from 0.0 < MW < 1.0) detected and located using automated phase picking. Waveformbased location methods enabled a lowering of the detection threshold and could be used at industrial sites to maximize the information of induced seismic events recorded by small, sparse seismometer arrays deployed at the surface.
31 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
5. Challenges and Perspectives
5.1. The Role of the Velocity Model
The velocity model is necessary for both seismic forward modeling and inversion. Waveformbased location methods need velocity models to compute the theoretical traveltime table or generate synthetic waveforms. For example, TRI and FWI require waveform modeling and extrapolation with a good velocity model, and inaccuracies in the velocity model will be accumulated and thus lead to location errors. The most common type of seismic data for velocity model calibration are controlled shots (usually perforation shot), whose source locations and origin time are known. As opposed to weak elastic anisotropy of the earth at the global scale, the anisotropy of tight reservoirs (e.g., kerogenrich shales) at the explorationscale may be strong. The high magnitude of observed anisotropy is another challenge that impedes reliable field microseismic application, and the smallscale natural and hydraulic fractures in these formations make the anisotropy even more complicated (Gajewski et al., 2009; Grechka, 2007; Grechka et al., 2011). Moreover, industrial activities (e.g., hydraulic stimulation) may produce or extend fractures that change the local seismic velocity (Fehler et al., 1998), which in turn affects seismic location. Another example illuminating the influence of the velocity models involves the role of multiple phases in the imaging process. Waveformbased location is sensitive to the velocity model, especially when using multiply scattered and reflected waves. While these later arrivals increase illumination of the source, they may also increase sensitivity of the source location to these scatters and diffractors in the model. A fundamental problem of downhole waveform location is the fact that the use of highfrequency radiation of the source is not accompanied by sufficient smallscale complexity of the utilized velocity model. Controlledsource fullwaveform investigations make use of lowfrequency information to stabilize the inversion process. However, currently, such lowfrequency contributions are not available in microseismic monitoring. Perhaps future research will focus on overcoming instrumental limitations.
The velocity model can be improved by taking advantage of both active and passive seismic data. Both the traveltime and waveform amplitude of recorded data contain information on velocity structure and can be utilized to constrain the inversion either tomographically or through migration principles (Grechka & Heigel, 2017). The joint estimation of both source characteristics and the seismic velocity distribution is expected to confront model uncertainties in a more effective manner (see sections 3.3 and 3.4). A potential problem of joint location and velocity inversion is the tradeoff between the source location, origin time, and velocity model, which might lead to significant bias in the inversion results if it is not treated properly (Jansky et al., 2010). In addition, in the context of microseismic monitoring, it is difficult to obtain an accurate tomographic model due to the poor coverage of the receiver arrays (e.g., downhole monitoring) or a limited number of welllocated microseismic events (Bardainne & Gaucher, 2010). Reservoir imaging using passive seismic data and seismic migration techniques can help resolve fieldscale velocity variations, and this emerging area still faces many challenges (Grechka et al., 2017; Huang & Zhu, 2019, see Section 5.4). The relative location method (Fitch, 1975; Spence, 1980) has been demonstrated to be an effective alternative method, which can compensate the implicit velocity errors and anomalies and create correct relative location of the located events. Waveformbased relative location efforts are likewise on the rise and promise to provide additional, experiencedriven constraints by relying on a catalog of formerly welllocated events. By introducing master events, F. Grigoli et al. (2016) improved the standard diffraction stacking to masterevent waveform stacking with additional traveltime correction terms estimated at each station. L. Li et al. (2016, 2018) proposed two novel waveformbased relative location methods, relative correlation stacking and hybrid correlation stacking, which are based on the differential traveltime from an event pair (i.e., a master event and a target event) to a receiver and double differential traveltime from an event pair to a receiver pair, respectively. The combination of relative location methods and waveformbased methods has great potential since it inherits advantages from both sides: Information from master events can help to compensate the negative effect of velocity model uncertainty and is also applicable to weak seismic event through waveform stacking, though the location results and selection criteria of master events need to be further investigated.
5.2. Resolution and Uncertainty
Location uncertainty estimates generally describe location as a Gaussian distribution with mean and standard deviation. Factors contributing to location uncertainty include the input data quality (e.g., SNR),
32 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
location algorithms (e.g., the imaging operator), the velocity model, and the coverage of the monitoring networks. The location uncertainty estimation for waveformbased methods is still in its early stages, far from being standardized, though there exist several attempts on this topic (see Table 1). For example, the imaging resolution or distribution can also provide preliminary knowledge about the location uncertainty (Anikiev et al., 2014; Kao & Shan, 2004), and the Bootstrap or Jackknife method can help to extract some statistical information of the location results from waveformbased methods (Grigoli et al., 2013; L. Li et al., 2018). According to several existing studies for PWS, the estimated location uncertainty ranges from tens of meters to about 200 m, which are comparable to the corresponding prevailing wavelengths (Grigoli et al., 2013; Pesicek et al., 2014; L. Li et al., 2018).
It may be argued that the image function, introduced in section 3 as a central ingredient of most waveform based location methods, can to some extent be viewed as a probability distribution of the source location. In analogy with the Gaussian model, the detected peak of the image function represents the most likely source location for the given observations. This previously noted “locate, then pick” mentality clearly distinguishes the source imaging methods from conventional inversion, in that the method itself leads to a landscape of possible solutions and their likeliness, rather than a specific solution only. As can be observed in section 4 for essentially all considered examples, the sharpness of the image function is strongly affected by limited and sparse apertures and directional biases, clearly reflecting respective location uncertainties.
All current methods of uncertainty assessment are limited to testing a fixed velocity model. We argue that the estimation of location uncertainty should also consider the uncertainty of the velocity model, which may be a main contributor (Gesret et al., 2015; Huang et al., 2017b). For traveltimebased methods, the velocity uncertainty can be implicitly contained in the observed or theoretical traveltime, while it is not straightforward to estimate velocity uncertainty with waveforms. One potential solution is adding perturbations to the theoretical traveltime table. Given the short development course of waveformbased methods and the difference of location uncertainty estimation between the two types of methods, more systematic benchmarking studies of seismic source location algorithms are in demand, though several recent comparisons of microseismic waveformbased and traveltimebased location algorithms can be found in recent literature (e.g., Pesicek et al., 2014; Trojanowski & Eisner, 2017; L. Li et al., 2018; Wuestefeld et al., 2018; Grigoli, Scarabello, et al., 2018). Note that full waveformbased methods have more sources of uncertainty from, for example, density, S wave velocity, and attenuation models. Compared to the high nonlinearity and requirement of a more accurate velocity model for other methods (see Table 2), we can expect that PWS methods have greater potential to provide an improved uncertainty when model uncertainty is included. This is also partially demonstrated with field applications in section 4. In addition, bandlimited waveforms are acquired and utilized for field application. Thus, the prevailing wavelength or the size of the Fresnel volume of the wavefield, along with the grid spacing, may roughly determine the range of the source imaging resolution and/or location uncertainty for waveformbased methods (Gajewski & Tessmer, 2005). Recent emergence of denser acquisitions across the scales may help solve the prevailing problem of undersampling and resolution. With this ongoing transition from waveforms to wavefields, the increasing volume of dense measurements promises to arrive at a more constrained and thus less uncertain reconstruction properties of a seismic source and the surrounding medium. Compared to local or regional scales, waveforms of small scale events and from the laboratory experiments are characterized by their highfrequency content, resulting in shorter prevailing wavelengths. In these scenarios, an increased image resolution can principally be achieved. However, the related wavefield complexity and increased sensitivity of the records act as additional sources of uncertainty and reasonably accurate, and detailed velocity information is required. This may be regarded as one of the major current challenges; waveformbased methods will have to overcome in order to be successful in AE investigations.
5.3. Complex Source Mechanisms
In contrast to controlled sources in active seismic exploration, both natural and induced/triggered (micro) earthquakes have more complex source mechanisms. Seismicity at different scales may be governed by different physical processes and involves different source mechanisms. Tectonic earthquakes are likely to rupture as shear failure, while volcanictectonic events maybe characterized by other processes, such as collapses, rock bursts, tensile cracks, and explosions (e.g., Cesca & Grigoli, 2015). The underlying mechanism of slow earthquakes, on the other hands, is still not fully understood. Seismicity induced by industrial
33 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
operations, for example, fluid injection or extraction, may involve more complex source mechanisms (from pure shear, over different combinations of shear and tensile to pure tensile) due to pore pressure and stress disturbances (e.g., Eaton et al., 2014; Eisner et al., 2010). Laboratory AE measurements of rock materials also contain information about hybrid failure mechanisms, besides simple tensile and shear type of failure (e.g., Ramsey & Chester, 2004; Q.Z. Zhu, 2017). Complex source mechanisms in anisotropic medium may result in apparent nonshear components of the source radiation. Consequently, the polarities of the amplitudes can be either positive or negative, depending on the direction of takeoff angles for each pair of the source and receiver positions. For example, a pure doublecouple strikeslip or dipslip type of source has such radiation patterns that the stacking value equals to zero at the true source position, if no polarity correction is utilized and the receivers are distributed symmetrically (e.g., Staněk et al., 2015; Zhebel & Eisner, 2015). As mentioned in section 3.1, various CFs have been applied to remove the effects of waveform polarity changes. However, a transformation of the original waveforms with CFs results in losing source mechanism information and reducing imaging resolution, and stacking may also be less effective in terms of noise elimination; that is, the detectability of weak events may be lower. The issues of complex source mechanisms also inhibit the applicability of those methods under the assumption of constant source characteristics (sources are nearly colocated), such as template matching method and matched field processing.
An alternative approach to compensate the effects of waveform polarity changes is joint inversion of location and source mechanism. The diffraction stacking method presented in Gajewski et al. (2007) has been further improved by Anikiev et al. (2014), Chambers et al. (2014), and Zhebel and Eisner (2015) through accounting for the source mechanism. They proposed to include polarity corrections in diffraction stacking through a joint location and moment tensor inversion. The full moment tensor representing source mechanism is inverted, and subsequently, amplitude polarities are corrected at every time before stacking. The method has been applied to field microseismic data acquired with a dense starlike surface array during hydraulic fracturing in shale gas reservoirs. Location quality was consistent with manual processing, and stability of inverted source mechanisms for the large surface arrays was proved by Staněk et al. (2014). They demonstrated a high quality of the source mechanisms even with a high level of noise in the data. Afterward, Liang et al. (2016) and Yu et al. (2018) proposed a similar joint inversion strategy by extending the sourcescanning algorithm to the joint sourcescanning algorithm, which scans both locations and focal mechanism. They carried out a theoretical and numerical analysis of the feasibility of the new method and demonstrated its superior performance in robustness and efficiency with field data from surface microseismic monitoring in the Sichuan Basin.
5.4. From Source Imaging to Reservoir Imaging
Unlike using manmade active seismic sources in exploration seismology, passive sources, such as natural and induced seismicity, ocean waves, and ambient seismic noise, are utilized in passive seismic imaging to image geological structure and locate seismic sources and underground oil, gas, or other resources (Daneshvar et al., 1995; Schuster et al., 2004). Once a seismic event is accurately located, it can resemble an active source used in exploration seismology if unknown source parameters are inverted. Then the advanced seismic migration techniques, along with the weak diffractions and reflections, can be used to image subsurface structures. For example, Dyer et al. (2008) and Reshetnikov et al. (2015) demonstrated the feasibility of seismic reflection imaging techniques in imaging fractures and faults in the vicinity of microseismic clouds in geothermal reservoirs, although uncertainty of these inversions remains questionable.
With the rise of interest in passive seismic monitoring of oilfields and production sites, the concepts of globalscale earthquake characterization techniques have gained increasing interest in the exploration seismology community. A technology exchange between earthquake seismology and exploration seismology has already emerged. It is quite promising to extend the waveformbased methods to reservoir imaging through linking active and passive seismic techniques, for example, using TRI methods. In exploration seismology, the process of RTM has been routinely used to image subsurface structures. With passive seismic records, TRI can be used to image the passive seismic source (e.g., extended finitefault source characteristics), as well as the surrounding fractures and structures. L. Li et al. (2015) applied TRI to synthetic microseismic events monitored with a surface array and a known velocity model, and it could simultaneously image seismic sources and surrounding reservoir structures. Lin and Zhang (2016) and Lin et al. (2018)
34 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
applied RTM to microseismic records containing strongly scattered seismic waves for a synthetic downhole monitoring system. Zhu and Sun (2016) and T. Zhu et al. (2019) proposed a workflow of datadriven and sourceindependent fracture imaging for passive seismic data. In this method, no source information is required, and zerolag cross correlations of the backpropagated direct and scattered wavefields are used to image the scatters and fractures. Numerical examples and an ultrasonic data set have indicated the approach's improved performance in imaging fractures compared to standard RTM, while a field microseismic data from the Marcellus Shale demonstrated the feasibility of the proposed method in imaging subsurface natural fractures (T. Zhu, 2019; Huang & Zhu, 2019).
Though only several exploratory studies have been conducted, passive seismic reservoir imaging bears great potential for providing new opportunities in reservoir characterization (Grechka & Heigel, 2017; B. Schwarz et al., 2017). However, there are still several technical issues that need further study, for example, the influence of irregular spatial distribution of sources and receivers and reliable wavefield separation from low SNR passive seismic data (Grechka & Heigel, 2017).
5.5. Computational Efficiency
Waveformbased location methods involve waveform storage and transmission, which raises new challenges for realtime application, especially for large monitoring arrays and large target regions (Pesicek et al., 2014; Xue et al., 2015; Li et al., 2017). For example, it took tens of minutes to several hours for PWS when hundreds of receivers and millions of grid points were considered (Li, Xie, et al., 2017), let alone TRI and FWI, which are more computational demanding (see Table 2). Currently, there are two solutions to this issue. One solution is taking advantage of highperformance computing technologies, and adopting parallelization for both waveform modeling and inversion process. Parallelization techniques based on software improvements (e.g., OpenMP and message passing interface) and hardware improvements (e.g., GPU) have been utilized in seismic forward modeling since the early 2000s (e.g., Bohlen, 2002; Michéa & Komatitsch, 2010; Thorbecke & Draganov, 2011). For most waveformbased methods, the discretized grid points representing potential source points are independent from each other. Thus, corresponding simulations and image function evaluations lend themselves naturally to parallel processing. Xue et al. (2015), Xue et al., 2016) described the performance of two new GPU features (HyperQ and GPU Direct, which are technologies introduced to improve the GPU performance and usability) applied in microseismic location. According to tests of a realistic 3D microseismic data set, they demonstrated that GPUaccelerated (require only several seconds with two NVidia Tesla C2070 [Fermi] graphics cards) TRI can be 30 times faster compared to an execution on a CPU (require several minutes with a personal workstation with an 8core Intel Xeon E2650 CPU). Though several minutes even hours are acceptable for academic or seismic hazards monitoring purposes, while (near)realtime process is highly required for field industry application and will also benefit the methodological study and nonindustry applications for natural seismicity monitoring (e.g., monitoring of seismic activity at volcanoes).
Another promising solution is incorporating optimization algorithms with waveformbased methods. Conventional gradientbased optimization algorithms are not very suitable for waveformbased methods, since the imaging function has local extrema and a decent initial model cannot be ensured. Fortunately, stochastic optimization algorithms are found to be powerful alternatives. Stochastic algorithms are generally inspired by natural laws or empirical principles and can be classified into global optimization algorithms due to their ability to search the extremum over all the search space, though they are not guaranteed to converge to the global optimum consistently. These strategies are derivativefree and only need random models as initial models; they include evolutionary algorithms (e.g., genetic algorithm), swarm intelligence (e.g., particle swarm optimization), and simulated annealing (e.g., Yang, 2010). Several algorithms have been proven to be effective for solving seismic source location problems. These optimization methods enable efficient convergence and reserve the event clustering by ensuring a highlevel of location accuracy. For example, Kennett et al. (2000) utilized the neighborhood algorithm to obtain reliable hypocenter and source mechanism estimation for a real earthquake. Gharti et al. (2010) completed pioneering work on incorporating differential evolution into waveformbased methods. Pesicek et al. (2014) and Verdon et al. (2017) successfully utilized a covariance matrix adaptation evolution strategy and neighborhood algorithm in diffraction stackingbased microseismic location, respectively. L. Li, Tan, Xie, et al. (2019) implemented three stochastic algorithms to improve the efficiency of waveformbased methods and proposed a parameter tuning
35 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
workflow comprising two types of repeated tests, which is supposed to enhance the algorithmic performance through evaluating the tradeoff between accuracy and efficiency.
The availability of waveformbased location methods toward realtime approaches can be improved with stateoftheart highperformance computing and advanced inversion strategies. On one hand, high performance computing facilities are not universally available, nor are the cloud computing services. On the other hand, the performance of optimization algorithms is not only problem dependent but also data dependent. Besides the algorithmic aspects, many other issues, such as the SNR, the sourcereceiver geometry, and the velocity model, also affect the final performance of optimization algorithms and need further investigation.
6. Conclusions
We present a comprehensive review of the current status of waveformbased seismic location methods and their applications at different scales, ranging from laboratory AEs to regional earthquakes. Seismic monitoring has been widely used in the wider areas of geophysics, seismology, AE, and engineering. In seismic monitoring, source location plays a central role, in that it provides fundamental information about crustal failure and active faults. Waveformbased location methods enable detection and location of weak seismic events recorded during underground activity monitoring, which are beneficial for characterizing tectonically and volcanologically active areas and unconventional resources development.
We elaborate on the evolution, principles, and categorization of waveformbased methods. Here, we classify them into four different categories as PWS, TRI, WT, and FWI. The framework of waveformbased location methods builds on the same theoretical foundation of seismic migration and imaging, and the location process can be regarded as backpropagation of seismic energy from the receivers to the source region. Specifically, the source is localized by focusing or reconstructing the source energy for discrete grid points with a certain migration or imaging operator. A significant advantage of waveformbased methods constitutes their noise resistance. Different preprocessing workflows have been introduced to address the challenges posed by wavefield interference, high noise levels, and polarity changes due to the radiation characteristics of the source. For example, conventional diffraction stacking generally utilizes singlechannelbased CFs, while crosscorrelation stacking explores mutual channel coherency. Waveform stacking methods generally make use of direct waves, while TRI and FWI take advantage of the entire waveforms including multiples and scattered waves. WT in turn integrates collective kinematic attributes (of wavefronts) and waveform information in a hybrid framework. Due to additional constraints from entire waveforms, more accurate source locations with higher certainty may be expected with reliable velocity models and phase recognitions. However, as the underlying problem is highly nonlinear, reliable starting conditions, that is, a priori information on the source and the velocity model, are required. Currently, it appears unlikely that passive seismic FWI or TRI may exist as a standalone tool but it may coexist with other techniques discussed in this review providing these initial conditions.
The feasibility and potential of waveformbased methods have been revealed through representative applications at multiple scales associated with natural and induced seismicity. So far, most field applications of waveformbased methods have been focused on PWS and at exploration and local scale, for example, for the monitoring of hydraulic fracturing in unconventional oil and gas reservoirs and local earthquakes. Meanwhile, these methods were effectively adapted to detect and locate lowfrequency seismic events recorded at local and regional scales. There also exist emerging applications of other methods in rockburst and AE events with comparably highfrequency content, though larger challenges of location uncertainty and velocity reliability are involved. As a general remark, the waveformbased theories and methodologies introduced in this review can be directly applied and/or easily adapted to seismic source location problems at different scales.
Finally, we discuss perspectives and summarize challenges that need to be addressed in future investigations. Arguably, the velocity model currently represents the Achilles heel of most seismic location methods, since it directly determines the reliability of location results. Strategies jointly inverting the location and velocity structure may help to confront this problem, though many challenges remain when field applications are concerned. Recent work on FWI and wavefront attributes of passive seismic events show promising results in constraining the velocity model while locating the source. The imaging resolution achievable by
36 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
waveformbased methods is largely determined by the prevailing wavelength of the events under consideration. The high computational efforts can be alleviated by making use of highperformance computing facilities and advanced optimization algorithms. Passive seismic reservoir imaging is another promising perspective that could provide new opportunities in reservoir characterization. Because of the increased demand for active CO2 reduction from the atmosphere, sequestration technologies may provide an important contribution to meet climate change targets in the future. Passive seismic monitoring is expected to play an important role in reaching public acceptance of such facilities and waveformbased methods most likely resemble the primary means of realizing such monitoring systems.
Glossary
Acoustic emission (AE): Seismic wave radiation from dynamic formation of a microcrack. Due to the small crack lengths, the frequencies of the radiated waves are generally in the kHz (acoustic or ultrasonic) range.
Backprojection: Use of the wave equation to track waves backward in time. Also called backpropagation.
Beamforming: A procedure for forming an array beam. Time delays are applied to the seismograms of an array to bring the desired signal recorded at different receivers or stations into phase, followed by summing of the outputs with appropriate weights. Also called delayandsum processing.
Characteristic function (CF): Nonlinear transformation of the seismogram, providing information that can be exploited by a picking algorithm for event detection and phase picking or directly by a stacking algorithm for event location.
Coherency: The property of two or more waveform traces having an inphase relationship. Also, a measure of the similarity of two or more functions.
Data domain: The data domain or space denotes the set of all possible values of the data.
Diffraction theory: A theory describing the principle and phenomenon of diffraction, which is defined as bending of waves around the corners of an obstacle or an aperture and propagating into the region of the obstacle's geometric shadow. The diffracting object or aperture effectively becomes a secondary source of the propagating wave.
Fault: A discontinuous surface separating two rock masses across which there has been significant displacement as a result of rock mass movement. Fault is also the place where earthquakes often happen.
Full waveform inversion (FWI): The process of using the whole waveform data instead of their particular parameter characteristics to estimate the model parameters of Earth structure and seismic source, by minimizing the misfit between the observed and synthetic seismograms.
Hydraulic fracturing (HF): The process of pumping fluid into a wellbore or the formation at high injection rate and pressure to fracture the rocks.
Image domain: The image domain or space denotes the set of all possible imaging values.
Inverse problem: The problem of determining the parameters of a physical system from measured data. The process of solving an inverse problem is called inversion, which involves searching for a distribution of parameters which fits the observed set of measurements. For instance, geophysical inversions generally derive a model to describe the subsurface that is consistent with the field data.
KirchhoffHelmholtz theory: A theory combines the Helmholtz equation (a timeindependent form of wave equation) with the Kirchhoff integral theorem (the wave function at a certain point can be represented as the integral/sum of wave contributions from the surroundings). It is expressed as KirchhoffHelmholtz representation/integral, which states that the elastic wavefield at any point within a subsurface region that can be reconstructed through a weighted integration of the wavefield and its normal derivative measured around a closed surface.
Kirchhoff migration: A migration technique integrates along diffraction curves and places the results at the crests of the diffraction curves.
Matched field processing: An arrayprocessing technique that originated from underwater acoustics. It aims to locate seismic noise sources through comparison of phase matches between measured and synthetic wavefields over a certain frequency bandwidth.
37 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
Microseismic event: The earthquake that has a magnitude less than 3. See also Seismic event.
Migration: An inversion operation involving rearrangement of seismic information elements so that reflections and diffractions are placed at their true locations. Also called imaging, the transformation of seismic data recorded as waveforms into a scaled version of the true geometry of subsurface geologic features that produced the recorded seismic waveforms.
Model domain: The model domain or space denotes the set of all possible values of model parameters.
Objective function: A function expressing the quality of all points in the model space. Usually objective functions are defined in terms of misfit functions, which quantify the disagreement between observed and theoretical data.
Origin time: The time of occurrence of initial energy release of a seismic event. Also called excitation time.
Partial waveform stacking (PWS): A subcategory of waveformbased seismic source location methods by stacking partial waveforms (generally the primary phases).
Passive seismic monitoring: Seismic investigations using passive monitoring or listening capabilities only, and no controlled seismic energy is stimulated actively by the investigator. Earthquake seismology, ambient noise studies, and microseismic monitoring are examples of passive seismic techniques, which are used for detection and characterization of fractures and faults, amplitude spectra of ground noise, and so forth.
Receiver: An observation point or equipment where ground motion is detected and a seismogram is recorded. Also called station.
Resolution: The ability to separate two features that are close together. The minimum separation of two bodies before their individual identities are lost. For instance, a high imaging resolution of the seismic source means the source energy is highly focused and well distinctive in the imaging profile.
Reverse time migration (RTM): A seismic migration method through backward propagation of the wavefield in time. The image of the subsurface is obtained by seeking the best match between the extrapolation of timereversed waveform data and the theoretical waveforms based on estimated velocity model and source parameters.
Seismic event: A coherent lineup on traces or an arrival of a new seismic wave, usually accompanied by a phase change and an increase in amplitude on a seismic record. A seismic event may indicate a reflection, refraction, diffraction, or other type of wavefront, as well as the movement of the earth resulting from the abrupt release of accumulated strain and energy, for example, tremors, natural earthquakes, microseismic events, and acoustic emissions, which may cover different magnitudes and frequencies and may happen naturally or triggered/induced artificially. See also seismicity and triggered/ induced seismicity.
Seismic source location: The spatial position and time of occurrence for a seismic event. Also refers to the process of locating a seismic event.
Seismic phase: An event on a seismogram section marking the arrival of a new wave, indicated by a change of period and amplitude.
Seismic tremor: A weak vibration of the ground usually lasting for a long time period. Tremor can be observed in a variety of tectonically active areas around the world and can be classified into different types according to their origin, such as tectonic tremor and volcanic tremor. See also Seismic event.
Seismicity: The likelihood and distribution of earthquakes. Also a general term for the activity (e.g., the frequency and intensity) of earthquakes.
Semblance: A measure of multichannel coherence, and it is basically the energy of the stack normalized by the mean energy of the components of the stack.
Shortterm average to longterm average ratio (STA/LTA): The average energy (including both signal and noise) in a shorttime window divided by the average energy of waveforms in a longtime window. Signaltonoise ratio (SNR): The energy or amplitude of the signal divided by all remaining (noise) energy or amplitude at the same channels and time window.
38 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
Reviews of Geophysics
10.1029/2019RG000667
Acknowledgments We thank the editors and three anonymous reviewers for their careful and constructive comments, which helped to improve the quality of the manuscript. We acknowledge Dirk Becker from University of Hamburg for providing Figure 11a. The work is sponsored by the National Natural Science Foundation of China (Grant 41872151), the National Thousand Young Talents Program (Jingqiang Tan, Grant 413170110), the Innovation driven Project of Central South University (Grant 502501005), Hunan Provincial Natural Science Foundation of China (Grant 2019JJ50762), the China Postdoctoral Science Foundation (Grant 2019M652803), and the Open Research Fund Program of State Key Laboratory of Acoustics, Chinese Academy of Sciences (Grant SKLA201911). Benjamin Schwarz received financial support by Geo. X, the Research Network for Geosciences in Berlin and Potsdam. Leon Diekmann is financially supported by the Utrecht Consortium for Subsurface Imaging (UCSI). All data are properly cited and can be accessed in the following references: Grigoli et al. (2013), Anikiev et al. (2014), Staněk et al. (2015), Poiata et al. (2016), LópezComino et al. (2017), L. Li et al. (2018), Poiata et al. (2018), and Shi, Nowacki, et al. (2019).
Source mechanism: The mechanism describes the deformation (motion) and state of the stress field in the source region that generate the recorded seismic waveforms. Also known as focal mechanism. Two different forms are generally used to express the source mechanism. First is the fault plane solution, which describes the initial slippage of the faulting with strike, dip, and rake. The second is the moment tensor, which represents the movement on a fault during an earthquake with a second order tensor, which consists of nine singlecouple forces in the Cartesian coordinate system.
Stacking: The process of summing several seismic traces with corrected with corrected time delays (as known as moveout).
Template matching: A filter that maximizes the output in response to a template signal. Filtering with a matched filter is equivalent to cross correlating with the signal, and it is the very powerful for identifying the presence of a given waveform or event in the presence of noise.
Time reverse imaging (TRI): An imaging technique for source or structure imaging through backpropagating the recorded full waveforms in time.
Tomography: The science and technology of creating images of the interior of a body. For instance, seismic tomography refers to an imaging method using a set of seismic traveltime, amplitude, or waveform data to deduce the velocity and/or attenuation/scattering structure of the Earth's interior.
Traveltime: The time of seismic waves to propagate from one point to another.
Triggered/Induced seismicity: Minor earthquakes and tremors that are caused by human activity that alters the stresses and strains in the Earth's crust. See also Seismicity.
Uncertainty: Random variation in the values assumed by a variable. For instance, the location uncertainty is embodied by a set of locations resulted by random variation of the input data.
Waveform: The complete analog or sufficiently dense sampled digital representation of a continuous wave group (e.g., of a seismic phase) or of a whole wave train. The seismic waveforms of a whole wave train are also called seismograms.
Wavefront: The surface with the same traveltime from a common source point.
Wavefront tomography (WT): A specific category of tomography techniques utilizing wavefront attributes (e.g., wavefront curvature and slopes).
References
Aki, K., & Richards, P. G. (2002). Quantitative seismology (Vol. 1). Sausalito: University Science Books. Akram, J., & Eaton, D. W. (2016). A review and appraisal of arrivaltime picking methods for downhole microseismic data. Geophysics,
81(2), KS71KS91. Anikiev, D., Valenta, J., Staněk, F., & Eisner, L. (2014). Joint location and source mechanism inversion of microseismic events:
Benchmarking on seismicity induced by hydraulic fracturing. Geophysical Journal International, 198(1), 249258. https://doi.org/ 10.1093/gji/ggu126 Arrowsmith, S., Young, C., Ballard, S., Slinkard, M., & Pankow, K. (2016). Pickless event detection and location: The Waveform Correlation EventDetection System (WCEDS) Revisited. Bulletin of the Seismological Society of America, 106(5), 20372044. https://doi.org/10.1785/ 0120160018 Artman, B., Podladtchikov, I., & Witten, B. (2010). Source location using timereverse imaging. Geophysical Prospecting, 58(5), 861873. Atkinson, G. M., Eaton, D. W., Ghofrani, H., Walker, D., Cheadle, B., Schultz, R., et al. (2016). Hydraulic fracturing and seismicity in the Western Canada Sedimentary Basin. Seismological Research Letters, 87(3), 631647. https://doi.org/10.1785/0220150263 Baggeroer, A. B., Kuperman, W. A., & Mikhalevsky, P. N. (1993). An overview of matched field methods in ocean acoustics. IEEE Journal of Oceanic Engineering, 18(4), 401424. https://doi.org/10.1109/48.262292 Baker, B. B., & Copson, E. T. (2003). The mathematical theory of Huygens principle. Providence: AMS Chelsea Publishing. Baker, T., Granat, R., & Clayton, R. W. (2005). Realtime earthquake location using Kirchhoff reconstruction. Bulletin of the Seismological Society of America, 95(2), 699707. Bao, X., & Eaton, D. W. (2016). Fault activation by hydraulic fracturing in western Canada. Science, 354(6318), 14061409. https://doi.org/ 10.1126/science.aag2583 Bardainne, T. (2011). Semiautomatic migration location of microseismic events: A good tradeoff between efficiency and reliability. Presented at the Third Passive Seismic Workshop Actively Passive 2011, Athens, Greece. https://doi.org/10.3997/22144609.20145331 Bardainne, T., & Gaucher, E. (2010). Constrained tomography of realistic velocity models in microseismic monitoring using calibration shots: Constrained tomography of realistic velocity models. Geophysical Prospecting, 58(5), 739753. https://doi.org/10.1111/j.1365 2478.2010.00912.x Bauer, A., Schwarz, B., & Gajewski, D. (2017). Utilizing diffractions in wavefront tomography. Geophysics, 82(2), R65R73. https://doi.org/ 10.1190/geo20160396.1
LI ET AL.
39 of 47
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
Bauer, A., Schwarz, B., Werner, T., & Gajewski, D. (2019). Unsupervised event identification and tagging for diffraction focusing. Geophysical Journal International, 217(3), 21652176. https://doi.org/10.1093/gji/ggz106
Baysal, E., Kosloff, D. D., & Sherwood, J. W. C. (1983). Reverse time migration. Geophysics, 48(11), 15141524. https://doi.org/10.1190/ 1.1441434
Beaucé, E., Frank, W. B., Paul, A., Campillo, M., & Hilst, R. D. (2019). Systematic detection of clustered seismicity beneath the Southwestern Alps. Journal of Geophysical Research Solid Earth, 124, 2019JB018110. https://doi.org/10.1029/2019JB018110
Beaucé, E., Frank, W. B., & Romanenko, A. (2018). Fast matched filter (FMF): An efficient seismic matchedfilter search for both CPU and GPU architectures. Seismological Research Letters, 89(1), 165172. https://doi.org/10.1785/0220170181
Benson, P. M., Vinciguerra, S., Meredith, P. G., & Young, R. P. (2008). Laboratory simulation of volcano seismicity. Science, 322(5899), 249252. https://doi.org/10.1126/science.1161927
Beskardes, G. D., Hole, J. A., Wang, K., Michaelides, M., Wu, Q., Chapman, M. C., et al. (2018). A comparison of earthquake backprojection imaging methods for dense local arrays. Geophysical Journal International, 212(3), 19862002. https://doi.org/10.1093/gji/ ggx520
Bischoff, S. T., Fischer, L., WehlingBenatelli, S., Fritschen, R., Meier, T., & Friederich, W. (2010). Spatiotemporal characteristics of mining induced seismicity in the eastern Ruhrarea. In Presented at the Cahiers du Centre Europe en de Geodynamics et de Seismologie. The European Center for: Geodynamics and Seismology (ECGS).
Block, L. V., Cheng, C. H., Fehler, M. C., & Phillips, W. S. (1994). Seismic imaging using microearthquakes induced by hydraulic fracturing. Geophysics, 59(1), 102112. https://doi.org/10.1190/geo19920156
Bohlen, T. (2002). Parallel 3D viscoelastic finite difference seismic modelling. Computers & Geosciences, 28(8), 887899. https://doi.org/ 10.1016/S00983004(02)000067
Bolton, D. C., Shokouhi, P., RouetLeduc, B., Hulbert, C., Rivière, J., Marone, C., & Johnson, P. A. (2019). Characterizing acoustic signals and searching for precursors during the laboratory seismic cycle using unsupervised machine learning. Seismological Research Letters, 90(3), 10881098. https://doi.org/10.1785/0220180367
Born, M., & Wolf, E. (1980). Principles of optics: Electromagnetic theory of propagation, interference and diffraction of light (6th ed). In Oxford (p. 836). New York: Pergamon Press.
Bozdağ, E., Peter, D., Lefebvre, M., Komatitsch, D., Tromp, J., Hill, J., et al. (2016). Global adjoint tomography: Firstgeneration model. Geophysical Journal International, 207(3), 17391766. https://doi.org/10.1093/gji/ggw356
Caffagni, E., Eaton, D. W., Jones, J. P., & van der Baan, M. (2016). Detection and analysis of microseismic events using a Matched Filtering Algorithm (MFA). Geophysical Journal International, 206(1), 644658.
Cesca, S., & Grigoli, F. (2015). Chapter twofull waveform seismological advances for microseismic monitoring. Advances in Geophysics, 56, 169228.
Chambers, K., Dando, B. D. E., Jones, G. A., Velasco, R., & Wilson, S. A. (2014). Moment tensor migration imaging. Geophysical Prospecting, 62(4), 879896. https://doi.org/10.1111/13652478.12108
Chambers, K., Kendall, J.M., BrandsbergDahl, S., & Rueda, J. (2010). Testing the ability of surface arrays to monitor microseismic activity: Testing the ability of surface arrays. Geophysical Prospecting, 58(5), 821830. https://doi.org/10.1111/j.1365 2478.2010.00893.x
Chattopadhyay, S., & McMechan, G. A. (2008). Imaging conditions for prestack reversetime migration. Geophysics, 73(3), S81S89. https:// doi.org/10.1190/1.2903822
Chmiel, M., Roux, P., & Bardainne, T. (2019). Highsensitivity microseismic monitoring: Automatic detection and localization of subsurface noise sources using matchedfield processing and dense patch arrays. Geophysics, 84(6), KS211KS223. https://doi.org/10.1190/ geo20180537.1
Cieplicki, R., Eisner, L., & Abbott, D. (2012). Correlation detection and location for microseismic events induced by hydraulic fracturing (p. D018). Presented at the 74th EAGE Conference and Exhibition incorporating EUROPEC 2012, Copenhagen, Denmark. https://doi.org/ 10.3997/22144609.20148208
Clarke, H., Eisner, L., Styles, P., & Turner, P. (2014). Felt seismicity associated with shale gas hydraulic fracturing: The first documented example in Europe. Geophysical Research Letters, 41(23), 83088314.
Cook, N. G. W. (1976). Seismicity associated with mining. Engineering Geology, 10(24), 99122. https://doi.org/10.1016/00137952(76) 900156
Corciulo, M., Roux, P., Campillo, M., Dubucq, D., & Kuperman, W. A. (2012). Multiscale matchedfield processing for noisesource localization in exploration geophysics. Geophysics, 77(5), KS33KS41. https://doi.org/10.1190/geo20110438.1
Cros, E., Roux, P., Vandemeulebrouck, J., & Kedar, S. (2011). Locating hydrothermal acoustic sources at Old Faithful Geyser using Matched Field Processing: Old Faithful viewed through array processing. Geophysical Journal International, 187(1), 385393. https://doi.org/ 10.1111/j.1365246X.2011.05147.x
Dahlen, F. A., Hung, S.H., & Nolet, G. (2000). Fréchet kernels for finitefrequency traveltimesI. Theory. Geophysical Journal International, 141(1), 157174. https://doi.org/10.1046/j.1365246X.2000.00070.x
Dales, P. (2018). Interferometric methods for seismic monitoring in industrial environments. Ottawa, CA: University of Ottawa. Retrieved from. http://ruor.uottawa.ca/handle/10393/38303
Dales, P., Audet, P., Olivier, G., & Mercier, J.P. (2017). Interferometric methods for spatio temporal seismic monitoring in underground mines. Geophysical Journal International, 210(2), 731742. https://doi.org/10.1093/gji/ggx189
Daneshvar, M. R., Clay, C. S., & Savage, M. K. (1995). Passive seismic imaging using microearthquakes. Geophysics, 60(4), 11781186. https://doi.org/10.1190/1.1443846
Diekmann, L., Schwarz, B., Bauer, A., & Gajewski, D. (2019). Source localisation and joint velocity model building using wavefront attributes. Geophysical Journal International, 219(2), 9951007. https://doi.org/10.1093/gji/ggz342
Dougherty, S. L., Cochran, E. S., & Harrington, R. M. (2019). The LArgen Seismic Survey in Oklahoma (LASSO) experiment. Seismological Research Letters, 90(5), 20512057.
Douglas, A. (1967). Joint epicentre determination. Nature, 215(5096), 4748. https://doi.org/10.1038/215047a0 Douma, J., Niederleithinger, E., & Snieder, R. (2015). Locating events using time reversal and deconvolution: Experimental application and
analysis. Journal of Nondestructive Evaluation, 34(1), 19. https://doi.org/10.1007/s109210150276x Drew, J., Leslie, H. D., Armstrong, P. N., & Michard, G. (2005). Automated microseismic event detection and location by continuous spatial
mapping (p. Paper 95513). Presented at the SPE Annual Technical Conference and Exhibition. Drew, J., White, R. S., Tilmann, F., & Tarasewicz, J. (2013). Coalescence microseismic mapping. Geophysical Journal International, 195(3),
17731785. https://doi.org/10.1093/gji/ggt331
40 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
Droznin, D. V., Shapiro, N. M., Droznina, S. Y., Senyukov, S. L., Chebrov, V. N., & Gordeev, E. I. (2015). Detecting and locating volcanic tremors on the Klyuchevskoy group of volcanoes (Kamchatka) based on correlations of continuous seismic records. Geophysical Journal International, 203(2), 10011010.
Duncan, P. M. (2005). Is there a future for passive seismic? First Break, 23(6), 111115. Duncan, P. M., & Eisner, L. (2010). Reservoir characterization using surface microseismic monitoring. Geophysics, 75(5), 75A13975A146. Duveneck, E. (2004). Velocity model estimation with dataderived wavefront attributes. Geophysics, 69(1), 265274. https://doi.org/
10.1190/1.1649394 Dyer, B. C., Schanz, U., Ladner, F., Häring, M. O., & Spillman, T. (2008). Microseismic imaging of a geothermal reservoir stimulation. The
Leading Edge, 27(7), 856869. Eaton, D. W. (2018). Passive seismic monitoring of induced seismicity: Fundamental principles and application to energy technologies, (1st
ed.). Cambridge: Cambridge University Press. https://doi.org/10.1017/9781316535547 Eaton, D. W., Igonin, N., Poulin, A., Weir, R., Zhang, H., Pellegrino, S., & Rodriguez, G. (2018). Induced seismicity characterization during
hydraulicfracture monitoring with a shallowwellbore geophone array and broadband sensors. Seismological Research Letters, 89(5), 16411651. https://doi.org/10.1785/0220180055 Eaton, D. W., van der Baan, M., Birkelo, B., & Tary, J.B. (2014). Scaling relations and spectral characteristics of tensile microseisms: Evidence for opening/closing cracks during hydraulic fracturing. Geophysical Journal International, 196(3), 18441857. Eisner, L., Abbott, D., Barker, W. B., Thornton, M. P., & Lakings, J. (2008). Noise suppression for detection and location of microseismic events using a matched filter (pp. 14311435). Presented at the 2008 SEG Annual Meeting, Society of Exploration Geophysicists. Eisner, L., Gei, D., Hallo, M., Opršal, I., & Ali, M. Y. (2013). The peak frequency of direct waves for microseismic events. Geophysics, 78(6), A45A49. https://doi.org/10.1190/geo20130197.1 Eisner, L., Hulsey, B. J., Duncan, P., Jurick, D., Werner, H., & Keller, W. (2010). Comparison of surface and borehole locations of induced seismicity. Geophysical Prospecting, 58(5), 809820. https://doi.org/10.1111/j.13652478.2010.00867.x Eisner, L., WilliamsStroud, S., Hill, A., Duncan, P., & Thornton, M. (2010). Beyond the dots in the box: Microseismicityconstrained fracture models for reservoir simulation. The Leading Edge, 29(3), 326333. https://doi.org/10.1190/1.3353730 Elsworth, D., Spiers, C. J., & Niemeijer, A. R. (2016). Understanding induced seismicity. Science, 354(6318), 13801381. https://doi.org/ 10.1126/science.aal2584 Esmersoy, C., & Miller, D. (1989). Backprojection versus backpropagation in multidimensional linearized inversion. Geophysics, 54(7), 921926. https://doi.org/10.1190/1.1442722 Etgen, J., Gray, S. H., & Zhang, Y. (2009). An overview of depth imaging in exploration geophysics. Geophysics, 74(6), WCA5WCA17. https://doi.org/10.1190/1.3223188 Fehler, M., House, L., Phillips, W. S., & Potter, R. (1998). A method to allow temporal variation of velocity in traveltime tomography using microearthquakes induced during hydraulic fracturing. Tectonophysics, 289(13), 189201. Fichtner, A., Bunge, H.P., & Igel, H. (2006). The adjoint method in seismology I. Theory. Physics of the Earth and Planetary Interiors, 157, 86104. https://doi.org/10.1016/j.pepi.2006.03.016 Fichtner, A., Kennett, B. L. N., Igel, H., & Bunge, H.P. (2009). Full seismic waveform tomography for uppermantle structure in the Australasian region using adjoint methods. Geophysical Journal International, 179(3), 17031725. https://doi.org/10.1111/j.1365 246X.2009.04368.x Fink, M. (1999). Timereversed acoustics. Scientific American, (November), 9197. Fitch, T. J. (1975). Compressional velocity in source regions of deep earthquakes: an application of the master earthquake technique. Earth and Planetary Science Letters, 26(2), 156166. Font, Y., Kao, H., Lallemand, S., Liu, C.S., & Chiao, L.Y. (2004). Hypocentre determination offshore of eastern Taiwan using the maximum intersection method. Geophysical Journal International, 158(2), 655675. Foulger, G. R., Wilson, M., Gluyas, J., Julian, B. R., & Davies, R. (2018). Global review of humaninduced earthquakes. EarthScience Reviews, 178, 438514. Frank, W. B., & Shapiro, N. M. (2014). Automatic detection of lowfrequency earthquakes (LFEs) based on a beamformed network response. Geophysical Journal International, 197(2), 12151223. https://doi.org/10.1093/gji/ggu058 Furumoto, M., Kunitomo, T., Inoue, H., Yamada, I., Yamaoka, K., Ikami, A., & Fukao, Y. (1990). Twin sources of highfrequency volcanic tremor of IzuOshima Volcano, Japan. Geophysical Research Letters, 17(1), 2527. https://doi.org/10.1029/GL017i001p00025 Gajewski, D., Anikiev, D., Kashtan, B., & Tessmer, E. (2007). Localization of seismic events by diffraction stacking. In SEG Technical Program Expanded Abstracts 2007 (pp. 12871291). Society of Exploration Geophysicists. Gajewski, D., Sommer, K., Vanelle, C., & Patzig, R. (2009). Influence of models on se
41 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
Goodfellow, S. D., & Young, R. P. (2014). A laboratory acoustic emission experiment under in situ conditions. Geophysical Research Letters, 41(10), 34223430.
Grandi, S., & Oates, S. J. (2009). Microseismic event location by crosscorrelation migration of surface array data for permanent reservoir monitoring (Vol. Extended Abstracts, p. X012). Presented at the 71st EAGE Conference & Technical Exhibition, EAGE.
Grechka, V. (2007). Multiple cracks in VTI rocks: Effective properties and fracture characterization. Geophysics, 72(5), D81D91. Grechka, V., De La Pena, A., SchisseléRebel, E., Auger, E., & Roux, P.F. (2015). Relative location of microseismicity. Geophysics, 80(6),
WC1WC9. Grechka, V., & Heigel, W. M. (2017). Microseismic monitoring. Tulsa: Society of Exploration Geophysicists. Grechka, V., Li, Z., Howell, B., Garcia, H., & Wooltorton, T. (2017). Highresolution microseismic imaging. The Leading Edge, 36(10),
822828. https://doi.org/10.1190/tle36100822.1 Grechka, V., Singh, P., & Das, I. (2011). Estimation of effective anisotropy simultaneously with locations of microseismic events.
Geophysics, 76(6), WC143WC155. Grigoli, F., Cesca, S., Amoroso, O., Emolo, A., Zollo, A., & Dahm, T. (2014). Automated seismic event location by waveform coherence
analysis. Geophysical Journal International, 196(3), 17421753. https://doi.org/10.1093/gji/ggt477 Grigoli, F., Cesca, S., Krieger, L., Kriegerowski, M., Gammaldi, S., Horalek, J., et al. (2016). Automated microseismic event location using
MasterEvent Waveform Stacking. Scientific Reports, 6, 25744. Grigoli, F., Cesca, S., Priolo, E., Rinaldi, A. P., Clinton, J. F., Stabile, T. A., et al. (2017). Current challenges in monitoring, discrimination,
and management of induced seismicity related to underground industrial activities: A European perspective. Reviews of Geophysics, 55 (2), 310340. https://doi.org/10.1002/2016RG000542 Grigoli, F., Cesca, S., Rinaldi, A. P., Manconi, A., LópezComino, J. A., Clinton, J. F., et al. (2018). The November 2017 M w 5.5 Pohang earthquake: A possible case of induced seismicity in South Korea. Science, 360(6392), 10031006. https://doi.org/10.1126/science. aat2010 Grigoli, F., Cesca, S., Vassallo, M., & Dahm, T. (2013). Automated seismic event location by traveltime stacking: An application to mining induced seismicity. Seismological Research Letters, 84(4), 666677. Grigoli, F., Scarabello, L., Böse, M., Weber, B., Wiemer, S., & Clinton, J. F. (2018). Pickand waveformbased techniques for realtime detection of induced seismicity. Geophysical Journal International, 213(2), 868884. Grosse, C., & Ohtsu, M. (2008). Acoustic emission testing. Berlin, Heidelberg: Springer Berlin Heidelberg. https://doi.org/10.1007/9783540 699729 Gu, C. (2016). Ground motions and source mechanisms of earthquakes in multiscales: Microseismicity to macroseismicity. Massachusetts Institute of Technology. Gutenberg, B., & Richter, C. F. (1956). Earthquake magnitude, intensity, energy, and acceleration: (Second paper). Bulletin of the Seismological Society of America, 46(2), 105145. Hansen, S. M., & Schmandt, B. (2015). Automated detection and location of microseismicity at Mount St. Helens with a largeN geophone array. Geophysical Research Letters, 42(18), 73907397. https://doi.org/10.1002/2015GL064848 Hardy, H. R. (2003). Acoustic emission/microseismic activity. Lisse: Balkema. Hassani, H., Hloušek, F., Alexandrakis, C., & Buske, S. (2018). Migrationbased microseismic event location in the SchlemaAlberoda mining area. International Journal of Rock Mechanics and Mining Sciences, 110, 161167. https://doi.org/10.1016/j. ijrmms.2018.08.002 Hauksson, E., & Shearer, P. (2005). Southern California hypocenter relocation with waveform crosscorrelation, Part 1: Results using the doubledifference method. Bulletin of the Seismological Society of America, 95(3), 896903. https://doi.org/10.1785/ 0120040167 Hendriyana, A., Bauer, K., Muksin, U., & Weber, M. (2018). AICbased diffraction stacking for local earthquake locations at the Sumatran Fault (Indonesia). Geophysical Journal International, 213(2), 952962. https://doi.org/10.1093/gji/ggy045 Höcht, G., de Bazelaire, E., Majer, P., & Hubral, P. (1999). Seismics and optics: Hyperbolae and curvatures. Journal of Applied Geophysics, 42(34), 261281. https://doi.org/10.1016/S09269851(99)000403 Holland, A. A. (2013). Earthquakes triggered by hydraulic fracturing in southcentral Oklahoma. Bulletin of the Seismological Society of America, 103(3), 17841792. Huang, C., Dong, L., Liu, Y., & Yang, J. (2017a). Acoustic waveequation based fullwaveform microseismic source location using improved scatteringintegral approach. Geophysical Journal International, 209(3), 14761488. https://doi.org/10.1093/gji/ggx087 Huang, C., Dong, L.G., Liu, Y.Z., & Yang, J.Z. (2017b). Waveformbased source location method using a source parameter isolation strategy. Geophysics, 82(5), KS85KS97. https://doi.org/10.1190/geo20170062.1 Huang, C., & Zhu, T. (2019). Passive seismic imaging of subsurface natural fractures: Application to Marcellus shale microseismic data. Geophysical Journal International, Accepted Manuscript.. https://doi.org/10.1093/gji/ggz214 Ide, S. (2010). Striations, duration, migration and tidal response in deep tremor. Nature, 466, 356. Igonin, N., & Innanen, K. A. (2018). Analysis of simultaneous velocity and source parameter updates in microseismic FWI. In SEG Technical Program Expanded Abstracts 2018 (pp. 10331037). Anaheim, California: Society of Exploration Geophysicists. https://doi. org/10.1190/segam20182992496.1 Ikelle, L. T., & Amundsen, L. (2018). Introduction to petroleum seismology. Society of Explora
42 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
Jay, J. A., Pritchard, M. E., West, M. E., Christensen, D., Haney, M., Minaya, E., et al. (2012). Shallow seismicity, triggered seismicity, and ambient noise tomography at the longdormant Uturuncu Volcano, Bolivia. Bulletin of Volcanology, 74(4), 817837. https://doi.org/ 10.1007/s0044501105687
Johnson, P. A., Savage, H., Knuth, M., Gomberg, J., & Marone, C. (2008). Effects of acoustic waves on stickslip in granular media and implications for earthquakes. Nature, 451(7174), 5760. https://doi.org/10.1038/nature06440
Kaderli, J., McChesney, M. D., & Minkoff, S. E. (2015). Microseismic event estimation in noisy data via full waveform inversion. In SEG Technical Program Expanded Abstracts 2015 (pp. 11591164). New Orleans, Louisiana: Society of Exploration Geophysicists. https://doi. org/10.1190/segam20155867154.1
Kaderli, J., McChesney, M. D., & Minkoff, S. E. (2018). A selfadjoint velocitystress fullwaveform inversion approach to microseismic source estimation. Geophysics, 83(5), R413R427. https://doi.org/10.1190/geo20170557.1
Kao, H., & Shan, S. J. (2004). The sourcescanning algorithm: Mapping the distribution of seismic sources in time and space. Geophysical Journal International, 157(2), 589594.
Kao, H., & Shan, S. J. (2007). Rapid identification of earthquake rupture plane using SourceScanning Algorithm. Geophysical Journal International, 168(3), 10111020.
Katsumata, A., & Kamaya, N. (2003). Lowfrequency continuous tremor around the Moho discontinuity away from volcanoes in the southwest Japan: DEEP TREMOR IN THE SOUTHWEST JAPAN. Geophysical Research Letters, 30(1), 201204. https://doi.org/ 10.1029/2002GL015981
Kawakatsu, H., & Montagner, J.P. (2008). Timereversal seismicsource imaging and momenttensor inversion. Geophysical Journal International, 175(2), 686688. https://doi.org/10.1111/j.1365246X.2008.03926.x
Kennett, B. L. N., MarsonPidgeon, K., & Sambridge, M. S. (2000). Seismic Source characterization using a neighbourhood algorithm. Geophysical Research Letters, 27(20), 34013404. https://doi.org/10.1029/2000gl011559
Kennett, B. L. N., & Sambridge, M. S. (1992). Earthquake location—Genetic algorithms for teleseisms. Physics of the Earth and Planetary Interiors, 75(13), 103110. https://doi.org/10.1016/00319201(92)90121B
Keranen, K. M., & Weingarten, M. (2018). Induced seismicity. Annual Review of Earth and Planetary Sciences, 2018, 149174. Kim, K.H., Ree, J.H., Kim, Y., Kim, S., Kang, S. Y., & Seo, W. (2018). Assessing whether the 2017 Mw 5.4 Pohang earthquake in South
Korea was an induced event. Science, 360(6392), 10071009. https://doi.org/10.1126/science.aat6081 Kirchhoff, G. (1883). Zur Theorie der Lichtstrahlen. Annalen der Physik, 254(4), 663695. https://doi.org/10.1002/andp.18832540409 Kiselevitch, V. L., Nikolaev, A. V., Troitskiy, P. A., & Shubik, B. M. (1991). Emission tomography: Main ideas, results, and prospects. In SEG
Technical Program Expanded Abstracts. Society of Exploration Geophysicists., 1991, 16021602. Kiser, E., & Ishii, M. (2013). Hidden aftershocks of the 2011 M w 9.0 Tohoku, Japan earthquake imaged with the backprojection method:
Hidden aftershocks. Journal of Geophysical Research Solid Earth, 118(10), 55645576. https://doi.org/10.1002/2013JB010158 KlemMusatov, K. (1994). Theory of seismic diffractions. In F. Hron & L. Lines (Eds.). (p. 431). Tulsa: Society of Exploration Geophysicists.
https://doi.org/10.1190/1.9781560802617 Köhn, D., De Nil, D., al Hagrey, S. A., & Rabbel, W. (2016). A combination of waveform inversion and reversetime modelling for micro-
seismic event characterization in complex salt structures. Environmental Earth Sciences, 75(18), 115. https://doi.org/10.1007/s12665 01660324 Kuperman, W. A., Hodgkiss, W. S., Song, H. C., Akal, T., Ferla, C., & Jackson, D. R. (1998). Phase conjugation in the ocean: Experimental demonstration of an acoustic timereversal mirror. The Journal of the Acoustical Society of America, 103(1), 2540. https://doi.org/ 10.1121/1.423233 Kwiatek, G., Bohnhoff, M., MartínezGarzón, P., Bulut, F., & Dresen, G. (2013). High resolution reservoir characterization using induced seismicity and state of the art waveform processing techniques. First Break, 31(7), 8188. Kwiatek, G., Plenkers, K., Dresen, G., & JAGUARS Research Group (2011). Source parameters of picoseismicity recorded at Mponeng deep gold mine, South Africa: Implications for scaling relations. Bulletin of the Seismological Society of America, 101(6), 25922608. Langet, N., Maggi, A., Michelini, A., & Brenguier, F. (2014). Continuous Kurtosisbased migration for seismic event detection and location, with application to Piton de la Fournaise Volcano, La Reunion. Bulletin of the Seismological Society of America, 104(1), 229246. https:// doi.org/10.1785/0120130107 Larmat, C., Tromp, J., Liu, Q., & Montagner, J.P. (2008). Time reversal location of glacial earthquakes. Journal of Geophysical Research, 113, B09314. https://doi.org/10.1029/2008JB005607 Legaz, A., Revil, A., Roux, P., Vandemeulebrouck, J., Gouédard, P., Hurst, T., & Bolève, A. (2009). Selfpotential and passive seismic monitoring of hydrothermal activity: A case study at Iodine Pool, Waimangu geothermal valley, New Zealand. Journal of Volcanology and Geothermal Research, 179(12), 1118. https://doi.org/10.1016/j.jvolgeores.2008.09.015 Lei, X., Huang, D., Su, J., Jiang, G., Wang, X., Wang, H., et al. (2017). Fault reactivation and earthquakes with magnitudes of up to Mw4. 7 induced by shalegas hydraulic fracturing in Sichuan Basin, China. Scientific Reports, 7(1), 7971. https://doi.org/10.1038/s41598017 08557y Li, K. L., Sadeghisorkhani, H., Sgattoni, G., Gudmundsson, O., & Roberts, R. (2017). Locating tremor using stacked products of correlations. Geophysical Research Letters, 44(7), 31563164. Li, K. L., Sgattoni, G., Sadeghisorkhani, H., Roberts, R., & Gudmundsson, O. (2017). A doublecorrelation tremorlocation method. Geophysical Journal International, 208(2), 12311236. https://doi.org/10.1093/gji/ggw453stack Li, L., Becker, D., Chen, H., Wang, X., & Gajewski, D. (2018). A systematic analysis of correlationbased seismic location methods. Geophysical Journal International, 212(1), 659678. https://doi.org/10.1093/gji/ggx436 Li, L., Chen, H., & Wang, X. (2016). Relative elastic interferometric imaging for microseismic source location. Journal of Geophysics and Engineering, 13(5), 733744. Li, L., Chen, H., & Wang, X.M. (2015). Weightedelasticwave interferometric imaging of microseismic source location. Applied Geophysics, 12(2), 221234. https://doi.org/10.1007/s117700150479z Li, L., Tan, J., Wood, D. A., Zhao, Z., Becker, D., Lyu, Q., et al. (2019). A review of the current status of induced seismicity monitoring for hydraulic fracturing in unconventional tight oil and gas reservoirs. Fuel, 242, 195210. https://doi.org/10.1016/j.fuel.2019.01.026 Li, L., Tan, J., Xie, Y., Tan, Y., Walda, J., Zhao, Z., & Gajewski, D. (2019). Waveformbased microseismic location using stochastic optimization algorithms: A parameter tuning workflow. Computers & Geosciences, 124, 115127. https://doi.org/10.1016/j.cageo.2019.01.002 Li, L., Xie, Y. J., Chen, H., Wang, X. M., & Gajewski, D. (2017). Improving the efficiency of microseismic imaging with particle swarm optimization. In 79th EAGE Conference and Exhibition 2017 (p. We B4 11). Li, T., Cai, M. F., & Cai, M. (2007). A review of mininginduced seismicity in China. International Journal of Rock Mechanics and Mining Sciences, 44(8), 11491171. https://doi.org/10.1016/j.ijrmms.2007
43 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
Li, Z., & van der Baan, M. (2016). Microseismic event localization by acoustic time reversal extrapolation. Geophysics, 81(3), KS123KS134. https://doi.org/10.1190/geo20150300.1
Liang, C., Yu, Y., Yang, Y., Kang, L., Yin, C., & Wu, F. (2016). Joint inversion of source location and focal mechanism of microseismicity. Geophysics, 81(2), KS41KS49.
Liao, Y.C., Kao, H., Rosenberger, A., Hsu, S.K., & Huang, B.S. (2012). Delineating complex spatiotemporal distribution of earthquake aftershocks: an improved SourceScanning Algorithm: Improved SourceScanning Algorithm. Geophysical Journal International, 189(3), 17531770. https://doi.org/10.1111/j.1365246X.2012.05457.x
Lin, F.C., Ritzwoller, M. H., & Snieder, R. (2009). Eikonal tomography: Surface wave tomography by phase front tracking across a regional broadband seismic array. Geophysical Journal International, 177(3), 10911110. https://doi.org/10.1111/j.1365246X.2009.04105.x
Lin, G., Shearer, P. M., & Hauksson, E. (2007). Applying a threedimensional velocity model, waveform cross correlation, and cluster analysis to locate southern California seismicity from 1981 to 2005. Journal of Geophysical Research, 112(B12), B12309. https://doi.org/ 10.1029/2007JB004986
Lin, Y., & Zhang, H. (2016). Imaging hydraulic fractures by microseismic migration for downhole monitoring system. Physics of the Earth and Planetary Interiors, 261, 8897.
Lin, Y., Zhang, H., & Jia, X. (2018). Targetoriented imaging of hydraulic fractures by applying the staining algorithm for downhole microseismic migration. Journal of Applied Geophysics, 150, 278283.
Lockner, D., Byerlee, J. D., Kuksenko, V., Ponomarev, A., & Sidorin, A. (1991). Quasistatic fault growth and shear fracture energy in granite. Nature, 350(6313), 3942.
Lomax, A., Michelini, A., & Curtis, A. (2009). Earthquake location, direct, globalsearch methods. In R. A. Meyers (Ed.), Encyclopedia of Complexity and Systems Science (pp. 24492473). New York, NY: Springer New York. https://doi.org/10.1007/9780387304403_150
LópezComino, J. A., Cesca, S., Heimann, S., Grigoli, F., Milkereit, C., Dahm, T., & Zang, A. (2017). Characterization of hydraulic fractures growth during the Äspö Hard Rock Laboratory Experiment (Sweden). Rock Mechanics and Rock Engineering, 50(11), 29853001. https:// doi.org/10.1007/s0060301712850
Ma, J., Dong, L., Zhao, G., & Li, X. (2019). Ground motions induced by mining seismic events with different focal mechanisms. International Journal of Rock Mechanics and Mining Sciences, 116, 99110. https://doi.org/10.1016/j.ijrmms.2019.03.009
Mair, K., Marone, C., & Young, R. P. (2007). Rate dependence of acoustic emissions generated during shear of simulated fault gouge. Bulletin of the Seismological Society of America, 97(6), 18411849. https://doi.org/10.1785/0120060242
Matoza, R. S., Shearer, P. M., Lin, G., Wolfe, C. J., & Okubo, P. G. (2013). Systematic relocation of seismicity on Hawaii Island from 1992 to 2009 using waveform cross correlation and cluster analysis. Journal of Geophysical Research Solid Earth, 118(5), 22752288. https://doi. org/10.1002/jgrb.50189
Maxwell, S. C. (2014). Microseismic imaging of hydraulic fracturing. Tulsa: Society of Exploration Geophysicists. McGarr, A., Simpson, D., Seeber, L., & Lee, W. H. K. (2002). Case histories of induced and triggered seismicity. In H. K. L. William, K.
Hiroo, C. J. Paul, & K. Carl (Eds.), International Handbook of Earthquake and Engineering Seismology (Vol. 81A, pp. 647664). San Francisco, Calif: Academic Press. McMechan, G. A. (1982). Determination of source parameters by wavefield extrapolation. Geophysical Journal International, 71(4), 613628. Meng, L., McGarr, A., Zhou, L., & Zang, Y. (2019). An investigation of seismicity induced by hydraulic fracturing in the Sichuan Basin of China based on data from a temporary seismic network. Bulletin of the Seismological Society of America, 109(1), 348357. https://doi.org/ 10.1785/0120180310 Michéa, D., & Komatitsch, D. (2010). Accelerating a threedimensional finitedifference wave propagation code using GPU graphics cards: Accelerating a wave propagation code using GPUs. Geophysical Journal International, 182(1), 389402. https://doi.org/10.1111/j.1365 246X.2010.04616.x Michel, O., & Tsvankin, I. (2014). Gradient calculation for waveform inversion of microseismic data in VTI media. Journal of Seismic Exploration, 23(3), 201217. Michel, O., & Tsvankin, I. (2017). Waveform inversion for microseismic velocity analysis and event location in VTI media. Geophysics, 82(4), WA95WA103. https://doi.org/10.1190/geo20160651.1 Milne, J. (1886). Earthquakes and other earth movements (Vol. 56). K. Paul, Trench. Nakata, N., & Beroza, G. C. (2016). Reverse time migration for microseismic sources using the geometric mean as an imaging condition. Geophysics, 81(2), KS51KS60. https://doi.org/10.1190/geo20150278.1 Neale, J., Harmon, N., & Srokosz, M. (2018). Improving microseismic P wave source location with multiple seismic arrays. Journal of Geophysical Research Solid Earth, 123(1), 476492. https://doi.org/10.1002/2017JB015015 Neidell, N. S., & Taner, M. T. (1971). Semblance and other coherency measures for multichannel data. Geophysics, 36(3), 482497. https:// doi.org/10.1190/1.1440186 Nguyen, L. T., Kocur, G. K., & Saenger, E. H. (2018). Defect mapping in pipes by ultrasonic wavefield crosscorrelation: A synthetic verification. Ultrasonics, 90, 153165. https://doi.org/10.1016/j.ultras.2018.06.014 Obara, K. (2002). Nonvolcanic deep tremor associated with subduction in southwest Japan. Science, 296(5573), 16791681. https://doi.org/ 10.1126/science.1070378 Obara, K., & Kato, A. (2016). Connecting slow earthquakes to huge earthquakes. Science, 353(6296), 253257. https://doi.org/10.1126/ science.aaf1512 Obert, L., & Duvall, W. I. (1945). The microseismic method of predicting rock failure in underground mining, Part 1, General method. Bureau of Mines Report of Investigations, 3797. Okada, Y., Kasahara, K., Hori, S., Obara, K., Sekiguchi, S., Fujiwara, H., & Yamamoto, A. (2004). Recent progress of seismic observation networks in Japan —Hinet, Fnet, KNET and KiKnet—. Earth, Planets and Space, 56(8), xvxxviii. https://doi.org/10.1186/ BF03353076 Pao, Y., & Varatharajulu, V. (1976). Huygens principle, radiation conditions, and integral formulas for the scattering of elastic waves. The Journal of the Acoustical Society of America, 59(6), 13611371. https://doi.org/10.1121/1.381022 Peng, Z., & Zhao, P. (2009). Migration of early aftershocks following the 2004 Parkfield earthquake. Nature Geoscience, 2(12), 877881. https://doi.org/10.1038/ngeo697 Pesicek, J. D., Child, D., Artman, B., & Cieślik, K. (2014). Picking versus stacking in a modern microearthquake location: Comparison of results from a surface passive seismic monitoring array in Oklahoma. Geophysics, 79(6), KS61KS68. Plessix, R.E. (2006). A review of the adjointstate method for computing the gradient of a functional with geophysical applications. Geophysical Journal International, 167(2), 495503. https://doi.org/10.1111/j.1365246X.2006.02978.x
44 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
Poiata, N., Satriano, C., Vilotte, J.P., Bernard, P., & Obara, K. (2016). Multiband array detection and location of seismic sources recorded by dense seismic networks. Geophysical Journal International, 205(3), 15481573.
Poiata, N., Vilotte, J.P., Bernard, P., Satriano, C., & Obara, K. (2018). Imaging different components of a tectonic tremor sequence in southwestern Japan using an automatic statistical detection and location method. Geophysical Journal International, 213(3), 21932213. https://doi.org/10.1093/gji/ggy070
Pujol, J. (2004). Earthquake location tutorial: Graphical approach and approximate epicentral location techniques. Seismological Research Letters, 75(1), 6374. https://doi.org/10.1785/gssrl.75.1.63
Pujol, Jose. (1992). Joint hypocentral location in media with lateral velocity variations and interpretation of the station corrections. Physics of the Earth and Planetary Interiors, 75(13), 724. https://doi.org/10.1016/00319201(92)90114B
RamosMartinez, J., & McMechan, G. A. (2001). Sourceparameter estimation by full waveform inversion in 3D heterogeneous, viscoelastic, anisotropic media. Bulletin of the Seismological Society of America, 91(2), 276291. https://doi.org/10.1785/0120000017
Ramsey, J. M., & Chester, F. M. (2004). Hybrid fracture and the transition from extension fracture to shear fracture. Nature, 428(6978), 6366. https://doi.org/10.1038/nature02333
Rentsch, S., Buske, S., Lüth, S., & Shapiro, S. A. (2007). Fast location of seismicity: A migrationtype approach with application to hydraulic fracturing data. Geophysics, 72(1), S33S40.
Reshetnikov, A., Kummerow, J., Asanuma, H., Häring, M., & Shapiro, S. A. (2015). Microseismic reflection imaging and its application to the Basel geothermal reservoir. Geophysics, 80(6), WC39WC49.
Rodriguez, I. V., Sacchi, M., & Gu, Y. J. (2012). Simultaneous recovery of origin time, hypocentre location and seismic moment tensor using sparse representation theory. Geophysical Journal International, 188(3), 11881202.
Ross, Z. E., Meier, M., Hauksson, E., & Heaton, T. H. (2018). Generalized seismic phase detection with deep learning. Bulletin of the Seismological Society of America, 108(5A), 28942901. https://doi.org/10.1785/0120180080
Ross, Z. E., Trugman, D. T., Hauksson, E., & Shearer, P. M. (2019). Searching for hidden earthquakes in Southern California. Science, 364(6442), 767771. https://doi.org/10.1126/science.aaw6888
Ross, Z. E., Yue, Y., Meier, M., Hauksson, E., & Heaton, T. H. (2019). PhaseLink: A deep learning approach to seismic phase association. Journal of Geophysical Research Solid Earth, 124(1), 856869. https://doi.org/10.1029/2018JB016674
Rost, S., & Thomas, C. (2002). Array seismology: Methods and applications. Reviews of Geophysics, 40(3), 21. Ruigrok, E., Gibbons, S., & Wapenaar, K. (2017). Crosscorrelation beamforming. Journal of Seismology, 21(3), 495508. https://doi.org/
10.1007/s1095001696126 Ruiz, S., AdenAntoniow, F., Baez, J. C., Otarola, C., Potin, B., del Campo, F., et al. (2017). Nucleation phase and dynamic inversion of the
M w 6.9 Valparaíso 2017 earthquake in Central Chile: M w 6.9 Valparaíso 2017 earthquake. Geophysical Research Letters, 44(20), 10,29010,297. https://doi.org/10.1002/2017GL075675 Rutledge, J. T., & Phillips, W. S. (2003). Hydraulic stimulation of natural fractures as revealed by induced microearthquakes, Carthage Cotton Valley gas field, east Texas. Geophysics, 68(2), 441452. Saenger, E. H., Kocur, G. K., Jud, R., & Torrilhon, M. (2011). Application of time reverse modeling on ultrasonic nondestructive testing of concrete. Applied Mathematical Modelling, 35(2), 807816. https://doi.org/10.1016/j.apm.2010.07.035 Sambridge, M., & Kennett, B. L. N. (1986). A novel method of hypocentre location. Geophysical Journal International, 87(2), 679697. https://doi.org/10.1111/j.1365246X.1986.tb06644.x Sambridge, M., & Mosegaard, K. (2002). Monte Carlo methods in geophysical inverse problems. Reviews of Geophysics, 40(3), 10091038. https://doi.org/10.1029/2000RG000089 Schneider, W. A. (1978). Integral formulation for migration in two and three dimensions. Geophysics, 43(1), 4976. https://doi.org/10.1190/ 1.1440828 Scholz, C. H. (1968). Microfracturing and the inelastic deformation of rock in compression. Journal of Geophysical Research, 73(4), 14171432. Schuster, G. T., Yu, J., & Sheng, J. (2004). Interferometric/daylight seismic imaging. Geophysical Journal International, 157(2), 838852. Schwarz, B. (2019). An introduction to seismic diffraction. Advances in Geophysics, 60, 164. https://doi.org/10.1016/bs.agph.2019.05.001 Schwarz, B., Bauer, A., & Gajewski, D. (2016). Passive seismic source localization via commonreflectionsurface attributes. Studia Geophysica et Geodaetica, 60(3), 531546. https://doi.org/10.1007/s112000150493x Schwarz, B., Bauer, A., & Gajewski, D. (2017). A unified inversion scheme for diffractions and passive events (p. WS05 B02). Presented at the 79th EAGE Conference and Exhibition 2017Workshops. Schwarz, B., & Gajewski, D. (2017). A generalized view on normal moveout. Geophysics, 82(5), V335V349. https://doi.org/10.1190/ geo20170159.1 Sellers, E. J., Kataka, M. O., & Linzer, L. M. (2003). Source parameters of acoustic emission events and scaling with mininginduced seismicity. Journal of Geophysical Research Solid Earth, 108(B9), 2418. https://doi.org/10.1029/2001jb000670 Senobari, N. S., Funning, G. J., Keogh, E., Zhu, Y., Yeh, C. M., Zimmerman, Z., & Mueen, A. (2019). SuperEfficient CrossCorrelation (SECC): A fast matched filtering code suitable for desktop computers. Seismological Research Letters, 90(1), 322334. https://doi.org/ 10.1785/0220180122 Shapiro, S. A. (2015). Fluidinduced seismicity. Cambridge: Cambridge University Press. Shearer, P., Hauksson, E., & Lin, G. (2005). Southern California hypocenter relocation with waveform crosscorrelation, Part 2: Results using sourcespecific station terms and cluster analysis. Bulletin of the Seismological Society of America, 95(3), 904915. https://doi.org/ 10.1785/0120040168 Shearer, P. M. (2009). Introduction to seismology (2nd ed.). Cambridge: Cambridge University Press. Shekar, B., & Sethi, H. S. (2019). Fullwaveform inversion for microseismic events using sparsity constraints. Geophysics, 84(2), KS1KS12. https://doi.org/10.1190/geo20170822.1 Shelly, D. R., Beroza, G. C., & Ide, S. (2007). Nonvolcanic tremor and lowfrequency earthquake swarms. Nature, 446(7133), 305307. https://doi.org/10.1038/nature05666 Shi, P., Angus, D., Rost, S., Nowacki, A., & Yuan, S. (2019). Automated seismic waveform location using multichannel coherency migration (MCM)I: theory. Geophysical Journal International, 216(3), 18421866. https://doi.org/10.1093/gji/ggy132 Shi, P., Nowacki, A., Rost, S., & Angus, D. (2019). Automated seismic waveform location using Multichannel Coherency Migration (MCM) —II. Application to induced and volcanotectonic seismicity. Geophysical Journal International, 216(3), 16081632. https://doi.org/ 10.1093/gji/ggy507 Sokos, E. N.,
45 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
Solano, E. A., Hjörleifsdóttir, V., & Liu, Q. (2017). Fullwaveform detection of nonimpulsive seismic events based on timereversal methods. Geophysical Journal International, 211(3), 13961413. https://doi.org/10.1093/gji/ggx372
Song, F., Warpinski, N. R., Toksöz, M. N., & Kuleli, H. S. (2014). Fullwaveform based microseismic event detection and signal enhancement: An application of the subspace approach: Microseismic detection and signal enhancement. Geophysical Prospecting, 62(6), 14061431. https://doi.org/10.1111/13652478.12126
Spence, W. (1980). Relative epicenter determination using Pwave arrivaltime differences. Bulletin of the Seismological Society of America, 70(1), 171183.
Stanchits, S., Burghardt, J., & Surdi, A. (2015). Hydraulic fracturing of heterogeneous rock monitored by acoustic emission. Rock Mechanics and Rock Engineering, 48(6), 25132527. https://doi.org/10.1007/s0060301508481
Staněk, F., Anikiev, D., Valenta, J., & Eisner, L. (2015). Semblance for microseismic event detection. Geophysical Journal International, 201(3), 13621369.
Staněk, F., & Eisner, L. (2017). Seismicity induced by hydraulic fracturing in shales: A bedding plane slip model. Journal Of Geophysical Research: Solid Earth, 122(10), 79127926.
Staněk, F., Eisner, L., & Jan Moser, T. (2014). Stability of source mechanisms inverted from Pwave amplitude microseismic monitoring data acquired at the surface: Stability of source mechanisms inverted from Pwave amplitude data. Geophysical Prospecting, 62(3), 475490. https://doi.org/10.1111/13652478.12107
Stein, S., & Wysession, M. (2003). An introduction to seismology, earthquakes, and earth structure. Malden, MA: Blackwell Pub. Steiner, B., Saenger, E. H., & Schmalholz, S. M. (2008). Time reverse modeling of lowfrequency microtremors: Application to hydrocarbon
reservoir localization. Geophysical Research Letters, 35, L03307. https://doi.org/10.1029/2007GL032097 Stump, B. W. (2002). Characterization of mining explosions at regional distances: Implications with the International Monitoring System.
Reviews of Geophysics, 40(4), 1011. https://doi.org/10.1029/1998RG000048 Suckale, J. (2009). Induced seismicity in hydrocarbon fields. Advances in Geophysics, 51, 55106. https://doi.org/10.1016/S00652687(09)
051073 Sugan, M., Kato, A., Miyake, H., Nakagawa, S., & Vuan, A. (2014). The preparatory phase of the 2009 M w 6.3 L'Aquila earthquake by
improving the detection capability of lowmagnitude foreshocks. Geophysical Research Letters, 41(17), 61376144. https://doi.org/ 10.1002/2014GL061199 Sun, J., Xue, Z., Zhu, T., Fomel, S., & Nakata, N. (2016). Fullwaveform inversion of passive seismic data for sources and velocities. In SEG Technical Program Expanded Abstracts 2016 (pp. 14051410). Dallas, Texas: Society of Exploration Geophysicists. https://doi.org/ 10.1190/segam201613959115.1 Sun, J., Zhu, T., Fomel, S., & Song, W.Z. (2015). Investigating the possibility of locating microseismic sources using distributed sensor networks. In SEG Technical Program Expanded Abstracts 2015 (pp. 24852490). New Orleans, Louisiana: Society of Exploration Geophysicists. https://doi.org/10.1190/segam20155888848.1 Taner, M. T., & Koehler, F. (1969). Velocity spectradigital computer derivation applications of velocity functions. Geophysics, 34(6), 859881. https://doi.org/10.1190/1.1440058 Thompson, G. (2014). Seismic Monitoring of Volcanoes. In M. Beer, I. A. Kougioumtzoglou, E. Patelli, & I. S.K. Au (Eds.), Encyclopedia of Earthquake Engineering (pp. 125). Berlin, Heidelberg, Berlin Heidelberg: Springer. https://doi.org/10.1007/9783642361975_411 Thorbecke, J. W., & Draganov, D. (2011). Finitedifference modeling experiments for seismic interferometry. Geophysics, 76(6), H1H18. https://doi.org/10.1190/geo20100039.1 Thurber, C. H., & Rabinowitz, N. (2000). Advances in seismic event location. Dordrecht: Springer. Trojanowski, J., & Eisner, L. (2017). Comparison of migrationbased location and detection methods for microseismic events. Geophysical Prospecting, 65(1), 4763. Trugman, D. T., & Shearer, P. M. (2017). GrowClust: A hierarchical clustering algorithm for relative earthquake relocation, with application to the Spanish Springs and Sheldon, Nevada, Earthquake Sequences. Seismological Research Letters, 88(2A), 379391. https://doi. org/10.1785/0220160188 Umlauft, J., & Korn, M. (2019). 3D fluid channel location from noise tremors using matched field processing. Geophysical Journal International, 219(3), 15501561. https://doi.org/10.1093/gji/ggz385 Verdon, J. P., Kendall, J., Hicks, S. P., & Hill, P. (2017). Using beamforming to maximise the detection capability of small, sparse seismometer arrays deployed to monitor oil field activities. Geophysical Prospecting, 65(6), 15821596. Virieux, J., & Operto, S. (2009). An overview of fullwaveform inversion in exploration geophysics. Geophysics, 74(6), WCC1WCC26. https://doi.org/10.1190/1.3238367 Vuan, A., Sugan, M., Amati, G., & Kato, A. (2018). Improving the detection of lowmagnitude seismicity preceding the Mw 6.3 L'Aquila earthquake: Development of a scalable code based on the cross correlation of template earthquakes. Bulletin of the Seismological Society of America, 108(1), 471480. https://doi.org/10.1785/0120170106 Waldhauser, F., & Ellsworth, W. L. (2000). A doubledifference earthquake location algorithm: Method and application to the northern Hayward fault, California. Bulletin of the Seismological Society of America, 90(6), 13531368. Waldhauser, F., & Schaff, D. P. (2008). Largescale relocation of two decades of Northern California seismicity using crosscorrelation and doubledifference methods. Journal of Geophysical Research Solid Earth, 113, B08311. https://doi.org/10.1029/2007JB005479 Wang, H., & Alkhalifah, T. (2018). Microseismic imaging using a source function independent full waveform inversion method. Geophysical Journal International, 214(1), 4657. https://doi.org/10.1093/gji/ggy121 Wang, J., Templeton, D. C., & Harris, D. B. (2015). Discovering new events beyond the catalogueapplication of empirical matched field processing to Salton Sea geothermal field seismicity. Geophysical Journal International, 203(1), 2232. https://doi.org/10.1093/gji/ggv260 Wapenaar, K., Draganov, D., & Snieder, R. (2010). Tutorial on seismic interferometry: Part 1—Basic principles and applications. Geophysics, 75(5), 75A19
46 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
LI ET AL.
Reviews of Geophysics
10.1029/2019RG000667
Withers, M., Aster, R., & Young, C. (1999). An automated local and regional seismic event detection and location system using waveform correlation. Bulletin of the Seismological Society of America, 89(3), 657669.
Withers, M., Aster, R., Young, C., Beiriger, J., Harris, M., Moore, S., & Trujillo, J. (1998). A comparison of select trigger algorithms for automated global seismic phase and event detection. Bulletin of the Seismological Society of America, 88(1), 95106.
Wu, S., Wang, Y., Zheng, Y., & Chang, X. (2018). Microseismic source locations with deconvolution migration. Geophysical Journal International, 212(3), 20882115. https://doi.org/10.1093/gji/ggx518
Wu, Y., & McMechan, G. A. (1996). Elastic fullwaveform inversion for earthquake source parameters. Geophysical Journal International, 127(1), 6174. https://doi.org/10.1111/j.1365246X.1996.tb01535.x
Wuestefeld, A., Näsholm, S. P., & Greve, S. M. (2018). How the choice of location algorithm affect results: A synthetic comparison. Presented at the Seventh EAGE Workshop on Passive Seismic 2018.
Xue, Q., Wang, Y., & Chang, X. (2016). Fast 3D elastic microseismic source location using new GPU features. Physics of the Earth and Planetary Interiors, 261, 2435. https://doi.org/10.1016/j.pepi.2016.08.001
Xue, Q., Wang, Y., Zhan, Y., & Chang, X. (2015). An efficient GPU implementation for locating microseismic sources using 3D elastic wave timereversal imaging. Computers & Geosciences, 82, 8997.
Yang, X.S. (2010). Natureinspired metaheuristic algorithms. Luniver press.stack the selected part of the waveformYilmaz, Ö. (2015). Engineering seismology with applications to geotechnical engineering. society of exploration geophysicists. https://doi.org/10.1190/ 1.9781560803300
Zang, A., Oye, V., Jousset, P., Deichmann, N., Gritto, R., McGarr, A., et al. (2014). Analysis of induced seismicity in geothermal reservoirs— An overview. Geothermics, 52, 621. https://doi.org/10.1016/j.geothermics.2014.06.005
Zang, A., Stephansson, O., Stenberg, L., Plenkers, K., Specht, S., Milkereit, C., et al. (2017). Hydraulic fracture monitoring in hard rock at 410 m depth with an advanced fluidinjection protocol and extensive sensor array. Geophysical Journal International, 208(2), 790813. https://doi.org/10.1093/gji/ggw430
Zeng, X., Zhang, H., Zhang, X., Wang, H., Zhang, Y., & Liu, Q. (2014). Surface Microseismic monitoring of hydraulic fracturing of a shale gas reservoir using shortperiod and broadband seismic sensors. Seismological Research Letters, 85(3), 668677.
Zhang, C., Qiao, W., Che, X., Lu, J., & Men, B. (2019). Automated microseismic event location by amplitude stacking and semblance. Geophysics, 84(6), KS191KS210. https://doi.org/10.1190/geo20180409.1
Zhang, H., Eaton, D. W., Rodriguez, G., & Jia, S. Q. (2019). Sourcemechanism analysis and stress inversion for hydraulicfracturing induced event sequences near Fox Creek, Alberta. Bulletin of the Seismological Society of America, 109(2), 636651. https://doi.org/ 10.1785/0120180275
Zhang, H., Nadeau, R. M., & Toksoz, M. N. (2010). Locating nonvolcanic tremors beneath the San Andreas fault using a stationpair double difference location method. Geophysical Research Letters, 37(13), L13304.
Zhang, H., & Thurber, C. H. (2003). Doubledifference tomography: The method and its application to the Hayward Fault, California. Bulletin of the Seismological Society of America, 93(5), 18751889. https://doi.org/10.1785/0120020190
Zhang, M., & Wen, L. (2015). An effective method for small event detection: Match and locate (M&L). Geophysical Journal International, 200(3), 15231537. https://doi.org/10.1093/gji/ggu466
Zhang, W., & Zhang, J. (2013). Microseismic migration by semblanceweighted stacking and interferometry (pp. 20452049). Presented at the SEG Technical Program Expanded Abstracts 2013. https://doi.org/10.1190/segam20130970.1
Zhang, Z., Rector, J. W., & Nava, M. J. (2017). Simultaneous inversion of multiple microseismic data for event locations and velocity model with Bayesian inference. Geophysics, 82(3), KS27KS39. https://doi.org/10.1190/geo20160158.1
Zhebel, O., & Eisner, L. (2015). Simultaneous microseismic event localization and source mechanism determination. Geophysics, 80(1), KS1KS9.
Zhou, H. (1994). Rapid threedimensional hypocentral determination using a master station method. Journal of Geophysical Research, 99(B8), 15439. https://doi.org/10.1029/94JB00934
Zhou, H.W., Hu, H., Zou, Z., Wo, Y., & Youn, O. (2018). Reverse time migration: A prospect of seismic imaging methodology. Earth Science Reviews, 179, 207227. https://doi.org/10.1016/j.earscirev.2018.02.008
Zhu, T. (2014). Timereverse modelling of acoustic wave propagation in attenuating media. Geophysical Journal International, 197(1), 483494. https://doi.org/10.1093/gji/ggt519
Zhu, T. (2019). Passive seismic imaging of subwavelength natural fractures: Theory and 2D synthetic and ultrasonic data tests. Geophysical Journal International, 216(3), 18311841. https://doi.org/10.1093/gji/ggy529
Zhu, T., & Sun, J. (2016). Datadriven diffraction imaging of fractures using passive seismic data. In SEG Technical Program Expanded Abstracts. Society of Exploration Geophysicists., 2016, 26792683.
Zhu, T., Sun, J., Gei, D., Carcione, J. M., Cance, P., & Huang, C. (2019). Hybrid multiplicative timereversal imaging reveals the evolution of microseismic events: Theory and field data test. Geophysics, 160. https://doi.org/10.1190/geo20180662.1
Zoback, M. D., & Gorelick, S. M. (2012). Earthquake triggering and largescale geologic storage of carbon dioxide. Proceedings of the National Academy of Sciences, 109(26), 10,16410,168. https://doi.org/10.1073/pnas.1202473109
Zoback, M. D., & Kohli, A. H. (2019). Unconventional reservoir geomechanics: Shale gas, tight oil, and induced seismicity. Cambridge: Cambridge University Press. https://doi.org/10.1017/9781316091869
Zou, Z., Zhou, H., & Gurrola, H. (2014). Reversetime imaging of a doublet of microearthquakes in the Three Gorges Reservoir region. Geophysical Journal International, 196(3), 18581868. https://doi.org/10.1093/gji/ggt499
Erratum
In the originally published version of this article, several typographical errors were present. As a result, Equation 2 published incorrectly. Additionally, equations 79 published erroneously as 911, and the section
numbers in section 4 were not in the appropriate order. The article has since been corrected, and this version
may be considered the authoritative version of record.
47 of 47
19449208, 2020, 1, Downloaded from https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019RG000667, Wiley Online Library on [29/02/2024]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License