Uncertainty in FPSs from moment-tensor inversion

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):  
◽  
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>


1982 ◽  
Vol 72 (2) ◽  
pp. 439-456
Author(s):  
Thorne Lay ◽  
Jeffrey W. Given ◽  
Hiroo Kanamori

Abstract The seismic moment and source orientation of the 8 November 1980 Eureka, California, earthquake (Ms = 7.2) are determined using long-period surface and body wave data obtained from the SRO, ASRO, and IDA networks. The favorable azimuthal distribution of the recording stations allows a well-constrained mechanism to be determined by a simultaneous moment tensor inversion of the Love and Rayleigh wave observations. The shallow depth of the event precludes determination of the full moment tensor, but constraining Mzx = Mzy = 0 and using a point source at 16-km depth gives a major double couple for period T = 256 sec with scalar moment M0 = 1.1 · 1027 dyne-cm and a left-lateral vertical strike-slip orientation trending N48.2°E. The choice of fault planes is made on the basis of the aftershock distribution. This solution is insensitive to the depth of the point source for depths less than 33 km. Using the moment tensor solution as a starting model, the Rayleigh and Love wave amplitude data alone are inverted in order to fine-tune the solution. This results in a slightly larger scalar moment of 1.28 · 1027 dyne-cm, but insignificant (&lt;5°) changes in strike and dip. The rake is not well enough resolved to indicate significant variation from the pure strike-slip solution. Additional amplitude inversions of the surface waves at periods ranging from 75 to 512 sec yield a moment estimate of 1.3 ± 0.2 · 1027 dyne-cm, and a similar strike-slip fault orientation. The long-period P and SH waves recorded at SRO and ASRO stations are utilized to determine the seismic moment for 15- to 30-sec periods. A deconvolution algorithm developed by Kikuchi and Kanamori (1982) is used to determine the time function for the first 180 sec of the P and SH signals. The SH data are more stable and indicate a complex bilateral rupture with at least four subevents. The dominant first subevent has a moment of 6.4 · 1026 dyne-cm. Summing the moment of this and the next three subevents, all of which occur in the first 80 sec of rupture, yields a moment of 1.3 · 1027 dyne-cm. Thus, when the multiple source character of the body waves is taken into account, the seismic moment for the Eureka event throughout the period range 15 to 500 sec is 1.3 ± 0.2 · 1027 dyne-cm.


2020 ◽  
Vol 110 (5) ◽  
pp. 2112-2123 ◽  
Author(s):  
Bernard Dost ◽  
Annemijn van Stiphout ◽  
Daniela Kühn ◽  
Marloes Kortekaas ◽  
Elmer Ruigrok ◽  
...  

ABSTRACT Recent developments in the densification of the seismic network covering the Groningen gas field allow a more detailed study of the connection between induced seismicity and reactivated faults around the gas reservoir at 3 km depth. With the reduction of the average station distance from 20 km to 4–5 km, a probabilistic full-waveform moment tensor inversion procedure could be applied, resulting in both improved hypocenter location accuracy and full moment tensor solutions for events of M≥2.0 recorded in the period 2016–2019. Hypocenter locations as output from the moment tensor inversion are compared to locations from the application of other methods and are found similar within 250 m distance. Moment tensor results show that the double-couple (DC) solutions are in accordance with the known structure, namely normal faulting along 50°–70° dipping faults. Comparison with reprocessed 3D seismic sections, extended to a depth of 6–7 km, demonstrate that (a) most events occur along faults with a small throw and (b) reactivated faults in the reservoir often continue downward in the Carboniferous underburden. From non-DC contributions, the isotropic (ISO) component is dominant and shows consistent negative values, which is expected in a compacting medium. There is some indication that events connected to faults with a large throw (&gt;70  m) exhibit the largest ISO component (40%–50%).


EKSPLORIUM ◽  
2021 ◽  
Vol 42 (2) ◽  
pp. 111
Author(s):  
Priyobudi Priyobudi ◽  
Mohamad Ramdhan

ABSTRAK. Keberadaan sesar aktif dengan pola sesar naik di daerah Plampang berhasil diungkap dari sebaran hiposenter terelokasi, hasil inversi momen tensor, dan pemodelan perubahan tegangan Coulomb. Studi ini juga berhasil mengungkap sumber gempa pada sesar aktif tersebut dengan kedalaman relatif dangkal yang bisa menjadi ancaman di Pulau Sumbawa jika magnitudo maksimumnya rilis di masa yang akan datang. Hasil relokasi hiposenter menunjukkan sebaran episenter berarah barat daya–timur laut. Hal ini didukung juga oleh hasil inversi momen tensor yang menunjukkan bidang sesar berarah barat daya–timur laut (N2240E) dengan dip cukup curam (490). Penampang seismisitas vertikal pada arah dip menunjukkan adanya pola sesar naik yang semakin landai seiring bertambahnya kedalaman. Bidang sesar yang landai menunjukkan struktur decollement pada kedalaman 10–15 km dan berangsur menjadi curam sebagai struktur splay fault pada kedalaman 0–10 km. Hal tersebut konsisten dengan hasil inversi momen tensor yang menunjukkan mekanisme pergerakan sesar naik terjadi pada kedalaman 7 km. Pemodelan perubahan tegangan Coulomb menunjukkan adanya penambahan stress di luar area bidang sesar sehingga memicu terjadinya aftershocks. Sebaran gempa susulan menunjukkan adanya bidang sesar hipotetik dengan panjang 19 km dan lebar 12 km. Sesar sebesar ini berpotensi membangkitkan gempa dengan kekuatan Mw 6,4. Gempa Sumbawa 13 Juni 2020 dengan magnitudo M 5,3 disebabkan oleh sebagian kecil aktivitas dari bidang sesar tersebut.ABSTRACT. The existence of an active fault with a reverse fault mechanism in the Plampang area is successfully delineated from the distribution of the relocated hypocenter, the moment tensor inversion, and the Coulomb stress changes. This study also reveals the source of the earthquake in the active fault with a relatively shallow depth which can be a threat on Sumbawa Island if the maximum magnitude is released in the future. Seismicity from hypocenter relocation shows the distribution of the epicenter with a southwest–northeast direction. It is also supported by the moment tensor inversion result which shows the fault plane trending southwest–northeast (N2240E) with a steep dip (490). The vertical section of seismicity in the dip direction shows that the slope of the plane has a lower angle with increasing depth. The lower angle of a fault plane shows a decollement structure at a depth of 10–15 km and gradually becomes steep as a splay fault structure at a depth of 0–10 km. It is consistent with the result of moment tensor inversion which shows the mechanism of a reverse fault that occurred at a depth of 7 km. The Coulomb stress changes show the stress increasing outside the fault plane area, which triggers aftershocks. The distribution of aftershocks shows a hypothetical fault plane of 19 km long and 12 km wide. A fault of this size has the potential to generate an earthquake with a magnitude maximum of Mw 6.4. The Sumbawa earthquake on June 13, 2020, having M 5.3 was caused by a small part of the activity from the fault.


Kappa Journal ◽  
2021 ◽  
Vol 5 (1) ◽  
pp. 31-39
Author(s):  
Sri Rizki Eka Putri ◽  
◽  
Hiden Hiden ◽  
Suhayat Minardi ◽  
◽  
...  

Nusa Tenggara is one of the areas with a high level of seismic activity in Indonesia because this area is located between the Indian Ocean plate, which moves northward and pushes the Eurasian plate. One method that is often used to determine an earthquake's epicentre is the Tensor Moment Inversion method. This study aims to determine the moment tensor magnitude of each earthquake event and determine earthquake characteristics based on the earthquake focus mechanism in Nusa Tenggara from 2018 to 2019. The earthquake with a magnitude of ≥ 5.7 SR and to find out the fault parameters, namely strike, dip, and rake using waveform data. One method that is often used to determine an earthquake's epicentre is the Tensor Moment Inversion method. The results showed that the fault planes formed were reverse faults and oblique faults. It has been calculated the moment tensor for each of the six components, namely Mxx, Myy, Mzz, Mxy, Myz and Mxz. From the results of the focal analysis of the 2018-2019 Nusa Tenggara earthquake mechanism, the values of the fault plane orientation parameters such as strike, dip and rake are obtained. For strikes in Nusa Tenggara on area 1, namely: 73° to 122°, Dip: 20° to 72° and Rake: 53° to 139°. While in field 2 for a strike, it is 232° to 280°, Dip 28° to 75°, Rake: 52° to 102°.


Sign in / Sign up

Export Citation Format

Share Document