The source mechanism analysis of significant local earthquake Mw 5.6 January 16, 2017 near Mt. Sibayak and Mt. Sinabung using moment tensor inversion of BMKG waveform data

2018 ◽  
Author(s):  
R. A. Prasetyo ◽  
N. Heryandoko ◽  
D. S. J. Sianipar ◽  
I. Suardi ◽  
S. Rohadi
Geophysics ◽  
2011 ◽  
Vol 76 (6) ◽  
pp. WC65-WC75 ◽  
Author(s):  
Jing Du ◽  
Norm R. Warpinski

Although microseismic monitoring of hydraulic fractures has primarily been concerned with the dimensions, complexity, and growth of fractures or fracture systems, there is an ever-increasing desire to extract more information about the hydraulic-fracturing and/or natural fractures from microseismic data. Source mechanism analysis, which is concerned with deducing details of the failure process from the microseismic waveform data, is, therefore, attracting more attention. However, most of the studies focus more on the moment-tensor inversion than on extracting fault-plane solutions (FPSs) from inverted moment tensors. The FPSs can be extracted from the inverted moment-tensor, but there remains a question regarding how errors associated with the inversion of the moment-tensor affect the accuracy of the FPSs. We examine the uncertainties of FPS, given the uncertainties of the amplitude data, by looking into the uncertainty propagation from amplitude data into the moment-tensor and then into the resultant FPS. The uncertainty propagation method will be demonstrated using two synthetic examples.


2021 ◽  
Author(s):  
◽  
Elizabeth de Joux Robertson

<p>The aim of this project is to enable accurate earthquake magnitudes (moment magnitude, MW) to be calculated routinely and in near real-time for New Zealand earthquakes. This would be done by inversion of waveform data to obtain seismic moment tensors. Seismic moment tensors also provide information on fault-type. I use a well-established seismic moment tensor inversion method, the Time-Domain [seismic] Moment Tensor Inversion algorithm (TDMT_INVC) and apply it to GeoNet broadband waveform data to generate moment tensor solutions for New Zealand earthquakes. Some modifications to this software were made. A velocity model can now be automatically used to calculate Green's functions without having a pseudolayer boundary at the source depth. Green's functions can be calculated for multiple depths in a single step, and data are detrended and a suitable data window is selected. The seismic moment tensor solution that has either the maximum variance reduction or the maximum double-couple component is automatically selected for each depth. Seismic moment tensors were calculated for 24 New Zealand earthquakes from 2000 to 2005. The Global CMT project has calculated CMT solutions for 22 of these, and the Global CMT project solutions are compared to the solutions obtained in this project to test the accuracy of the solutions obtained using the TDMT_INVC code. The moment magnitude values are close to the Global CMT values for all earthquakes. The focal mechanisms could only be determined for a few of the earthquakes studied. The value of the moment magnitude appears to be less sensitive to the velocity model and earthquake location (epicentre and depth) than the focal mechanism. Distinguishing legitimate seismic signal from background seismic noise is likely to be the biggest problem in routine inversions.</p>


2021 ◽  
Author(s):  
Andreas Steinberg ◽  
Hannes Vasyura-Bathke ◽  
Peter Gaebler ◽  
Lars Ceranna

&lt;p&gt;&lt;span&gt;Estimating fast earthquakes&amp;#8217; source mechanism is essential for near real time hazard assessments, which are based on shakemaps and further downstream analysis such as physics based aftershock probability calculations. The model and data uncertainties associated to the estimated source mechanism are also crucial. We propose a Baysian Machine Learning algorithm trained on normalized synthetic waveforms for estimating the full moment tensor of earthquakes almost instantaneously with associated source parameter uncertainties. &lt;/span&gt;&lt;/p&gt;&lt;p&gt;&lt;span&gt;A prior assumption is an appropriate location of the earthquakes along with its associated uncertainties. Here, this is obtained by already established Machine learning based algorithms, where the training data set is computed by forward calculations of synthetic waveforms based on Green&amp;#8217;s functions calculated for a specified 1-D velocity model using the Pyrocko software package. The &lt;/span&gt;&lt;span&gt;learned labels, &lt;/span&gt;&lt;span&gt;which are the information learned by the Machine Learning algorithm associated to the data,&lt;/span&gt;&lt;span&gt; are the moment tensor components, described with only five unique parameters. For predefined locations in an area of interest we train a full independent Bayesian Convolutional Neural Network (BNN). &lt;/span&gt;&lt;/p&gt;&lt;p&gt;&lt;span&gt;With variational inference the weights of the network are not scalar but represent a distribution &lt;/span&gt;&lt;span&gt;of weights for the activation of neurons. &lt;/span&gt;&lt;span&gt; Each evaluation of input data into our BNN yields therefore to a set of predictions with associated probabilities. This allows us to evaluate an ensemble of possible source mechanisms for each evaluation of input waveform data. &lt;/span&gt;&lt;/p&gt;&lt;p&gt; &lt;span&gt;As a test set, we trained our models for an area south of the Coso geothermal field in California for a fixed set of broadband stations at maximum 150 km distance. We validate our approach with a subset of earthquakes from the Ridgecrest 2019-2020 sequence. For this data set we compare the results of the estimates of our Machine Learning based approach with independently determined focal mechanism and moment tensors. Overall, we benchmark our approach with&lt;/span&gt;&lt;span&gt; data &lt;/span&gt;&lt;span&gt;unseen during the training process &lt;/span&gt;&lt;span&gt;of the Machine Learning models&lt;/span&gt;&lt;span&gt; and show its capabilities for generating similar source mechanism estimations as independent studies within only a few seconds processing time per earthquake. &lt;/span&gt;&lt;span&gt;We finally apply the method to seismic data of a research network monitoring the area around two south-german geothermal power plants.&lt;/span&gt;&lt;span&gt; Our approach demonstrates the potential of Machine Learning for being implemented in operational frameworks for fast earthquake source mechanism estimation with associated uncertainties.&lt;/span&gt;&lt;/p&gt;&lt;p&gt;&lt;br&gt;&lt;br&gt;&lt;/p&gt;


2021 ◽  
Vol 73 (03) ◽  
pp. 31-33
Author(s):  
Pat Davis Szymczak

Real-time analysis of microseismic events using data gathered during hydraulic fracturing can give engineers critical feedback on whether a particular fracturing job has achieved its goal of increasing porosity and permeability and boosting stimulated reservoir volume (SRV). Currently, no perfect way exists to understand clearly if a fracturing operation has had the intended effect. Engineers collect data, but the methods used to gather it, manually sort it, and analyze it provide an inconclusive picture of what really is happening underground. Daniel Stephen Wamriew, a PhD candidate at the Skolkovo Institute of Science and Technology (Skoltech) in Moscow, said he believes this can change with advances in artificial intelligence and machine learning that can enhance accuracy in determining the location of a microseismic event while obtaining stable source mechanism solutions, all in real time. Wamriew presented his research at the 2020 SPE Russian Petroleum Technology Conference in Moscow in October in paper SPE 201925, “Deep Neural Network for Real-Time Location and Moment Tensor Inversion of Borehole Microseismic Events Induced by Hydraulic Fracturing.” The paper’s coauthors included Marwan Charara, Aramco Research Center, and Evgenii Maltsev, Skolkovo Institute of Science and Technology. Skoltech is a private institute established in 2011 as part of a multiyear partnership with the Massachusetts Institute of Technology. “People in the field mainly want to know if they created more fractures and if the fractures are connected,” Wamriew explained in a recent interview with JPT. “So, we need to know where exactly the fractures are, and we need to know the orientation (the source mechanism).” It Starts With Data “Usually, when you do hydraulic fracturing, a lot of data comes in,” Wamriew said. “It is not easy to analyze this data manually because you have to choose what part of the data you deal with, and, in doing that, you might leave out some necessary data that the human eye has missed.” To solve this problem, Wamriew proposes feeding microseismic data gathered during a fracturing job into a convolutional neural network (CNN) that he is constructing (Fig. 1). Humans discard nothing. Wave signals from actual events along with noise of all kinds goes into a machine, and the CNN delivers valuable information to reservoir engineers who want to understand the likely SRV. Companies today can identify the location of microseismic events, even without the help of artificial intelligence—though the techniques are always open to refinement—but analyzing the orientation (and hence their understanding of whether and how the fractures are connected) is a difficult and often expensive task that is usually left undone. “Current source mechanism solutions are largely inconsistent,” Wamriew said. “One scientist collects data and performs the moment tensor inversion, and another does the same and gets different results, even if they both use the same algorithm. When we handle data manually, we choose the process, and, in doing so, we introduce errors at every step because we are truncating, rounding up, and rounding down. We end up with something far from reality.”


2021 ◽  
Author(s):  
◽  
Elizabeth de Joux Robertson

<p>The aim of this project is to enable accurate earthquake magnitudes (moment magnitude, MW) to be calculated routinely and in near real-time for New Zealand earthquakes. This would be done by inversion of waveform data to obtain seismic moment tensors. Seismic moment tensors also provide information on fault-type. I use a well-established seismic moment tensor inversion method, the Time-Domain [seismic] Moment Tensor Inversion algorithm (TDMT_INVC) and apply it to GeoNet broadband waveform data to generate moment tensor solutions for New Zealand earthquakes. Some modifications to this software were made. A velocity model can now be automatically used to calculate Green's functions without having a pseudolayer boundary at the source depth. Green's functions can be calculated for multiple depths in a single step, and data are detrended and a suitable data window is selected. The seismic moment tensor solution that has either the maximum variance reduction or the maximum double-couple component is automatically selected for each depth. Seismic moment tensors were calculated for 24 New Zealand earthquakes from 2000 to 2005. The Global CMT project has calculated CMT solutions for 22 of these, and the Global CMT project solutions are compared to the solutions obtained in this project to test the accuracy of the solutions obtained using the TDMT_INVC code. The moment magnitude values are close to the Global CMT values for all earthquakes. The focal mechanisms could only be determined for a few of the earthquakes studied. The value of the moment magnitude appears to be less sensitive to the velocity model and earthquake location (epicentre and depth) than the focal mechanism. Distinguishing legitimate seismic signal from background seismic noise is likely to be the biggest problem in routine inversions.</p>


2021 ◽  
Vol 8 (1) ◽  
Author(s):  
Hijrah Saputra ◽  
Wahyudi Wahyudi ◽  
Iman Suardi ◽  
Ade Anggraini ◽  
Wiwit Suryanto

AbstractThis study comprehensively investigates the source mechanisms associated with the mainshock and aftershocks of the Mw = 6.3 Yogyakarta earthquake which occurred on May 27, 2006. The process involved using moment tensor inversion to determine the fault plane parameters and joint inversion which were further applied to understand the spatial and temporal slip distributions during the earthquake. Moreover, coseismal slip distribution was overlaid with the relocated aftershock distribution to determine the stress field variations around the tectonic area. Meanwhile, the moment tensor inversion made use of near-field data and its Green’s function was calculated using the extended reflectivity method while the joint inversion used near-field and teleseismic body wave data which were computed using the Kikuchi and Kanamori methods. These data were filtered through a trial-and-error method using a bandpass filter with frequency pairs and velocity models from several previous studies. Furthermore, the Akaike Bayesian Information Criterion (ABIC) method was applied to obtain more stable inversion results and different fault types were discovered. Strike–slip and dip-normal were recorded for the mainshock and similar types were recorded for the 8th aftershock while the 9th and 16th June were strike slips. However, the fault slip distribution from the joint inversion showed two asperities. The maximum slip was 0.78 m with the first asperity observed at 10 km south/north of the mainshock hypocenter. The source parameters discovered include total seismic moment M0 = 0.4311E + 19 (Nm) or Mw = 6.4 with a depth of 12 km and a duration of 28 s. The slip distribution overlaid with the aftershock distribution showed the tendency of the aftershock to occur around the asperities zone while a normal oblique focus mechanism was found using the joint inversion.


2020 ◽  
Author(s):  
Gene Aaron Ichinose ◽  
Sean Ricardo Ford ◽  
Robert J. Mellors

Sign in / Sign up

Export Citation Format

Share Document