Moment Tensor Inversion Research Articles

Complex Multi‐Fault Dynamics in Sikkim Himalaya: New Insights From Local Earthquake Analysis

AbstractAnomalously overturned thrust faults, lineaments and segmentation causing cross‐cutting basement structures characterize the tectonic setting of Sikkim Himalaya. However, its seismotectonics is poorly constrained along with the speculated northward extension of the Dhubri‐Chungthang Fault Zone (DCFZ) causing segmentation. Here, we utilize the precise location of newly acquired local earthquake data and fault plane solutions using full‐waveform moment tensor inversion to better constrain seismically active zones. Transtensional shearing along the Main Himalayan Thrust in central Sikkim is possibly incited by fluid‐rich upper‐crust. Cessation of the mapped 20 km wide mid‐crustal seismogenic zone of DCFZ at Chungthang and its northward discontinuation into the Higher Himalayas is confirmed by the striking variation in focal mechanisms. Earthquakes along imbricated segments in the lower‐crust originate possibly in response to crustal shortening. Extensional shearing along the Moho triggers seismicity to the northwest of Sikkim. Such complex tectonic dynamics instigating persistent seismicity indicates high potential for future great earthquakes in Sikkim Himalaya.

Тензор сейсмічного моменту та механізм вогнища землетрусу у Східній Словаччині 9 жовтня 2023 року

The accuracy of the focal mechanism solution mainly depends on the number of stations used and becomes problematic especially in the case of weak earthquakes and sparse networks. In our study, we retrieve the seismic moment tensor of the M=5 earthquake on 9 October 2023 (18:23:09 UTC, 21.783°E, 49.086°N, depth 11.5 km) in Eastern Slovakia from its records at only four seismic stations. Our seismic moment tensor inversion is based on the point source approach and the use of only direct waves calculated by the matrix method. Displacements on the surface of an elastic, horizontally-layered medium are generated using the frequency-wavenumber integration technique. The advantage of using only direct P- and S- waves in our inversion method is that they are less sensitive to path effects compared to reflected and converted waves, which reduces the impact of an inaccurate velocity structure and improves the accuracy and reliability of the result. Based on the forward modeling, a numerical technique was developed for the inversion of the observed waveforms for the components of the moment tensor M(t) using the generalized inversion solution. Before applying our method to the earthquake of October 9, 2023, it was also tested on the April 23, 2020 earthquake (23:18:26.42 UTC, 21.945°E, 48.781°N, magnitude M=5, depth 9 km), also in Eastern Slovakia, using data from only three stations. The resulting versions of the mechanism compare well with a very reliable version previously determined from the polarities of the first P-waves at a much larger number of stations, which only confirms the reliability of our inversion method and the very possibility of obtaining useful results from data of only limited number of stations.

Gradient-based joint inversion of point-source moment tensor and station-specific time-shifts

SUMMARY The misalignment of the observation and predicted waveforms in regional moment tensor inversion is mainly due to seismic models’ incomplete representation of the Earth's heterogeneities. Current moment tensor inversion techniques, allowing station-specific time-shifts to account for the model error, are computationally expensive. Here, we propose a gradient-based method to jointly invert moment-tensor parameters, centroid depth and unknown station-specific time-shifts utilizing the modern functionalities in deep learning frameworks. A $L_2^2$ misfit function between predicted synthetic and time-shifted observed seismograms is defined in the spectral domain, which is differentiable to all unknowns. The inverse problem is solved by minimizing the misfit function with a gradient descent algorithm. The method's feasibility, robustness and scalability are demonstrated using synthetic experiments and real earthquake data in the Long Valley Caldera, California. This work presents an example of fresh opportunities to apply advanced computational infrastructures developed in deep learning to geophysical problems.

A UK MW catalogue derived from coda envelopes

Summary The United Kingdom (UK) experiences low-to-moderate levels of seismicity; only 12 onshore earthquakes with local magnitude (ML) ≥ 4.0 have been recorded in the past 20 years. It is therefore difficult to estimate moment magnitude (MW) using conventional moment tensor inversion for the majority of UK seismicity, resulting in limited reliable estimates of MW. To address this, we calibrated coda envelopes at 16 broadband seismic stations distributed across the UK, to produce a MW catalogue for 100 events with MW≥2.13 that occurred since 2006. This was achieved using the open-source Coda Calibration Tool, which requires independent source parameter estimates for calibration. For 13 UK events between 2006 and 2022, we used spectral modelling to estimate apparent stress (0.32 to 1.74 MPa), and moment tensor inversion to estimate MW (3.35 to 4.52). These independent source parameters formed a subset of the inputs into the final calibration, which used seismic data from 33 events with coda-derived values of 2.57$\le $Mw$\le $4.49. The resultant coda calibration parameters were applied to 67 further events (MW≥2.13). The coda envelopes exhibit slow seismic coda decay across the UK, with significant energy up to 20 Hz, consistent with other regions of low tectonic activity. This MW catalogue, and the application of the calibration to future UK seismic events, will be useful for both assessing seismic hazard and event characterisation.

Crustal deformation from GNSS measurement and earthquake mechanism along Pieniny Klippen Belt, Southern Poland

The dynamic geological features of the Pieniny Klippen Belt (PKB) in southern Poland are nowadays a focal point of researchers as it is recognised as an active zone of crustal discontinuity. In the present study, we employed long-term analysis of ground-based global navigation satellite system (GNSS) measurements (from 2004 to 2020) to probe the crustal deformation, strain rates, and rotational rates in the PKB unit and the surrounding region, i.e., Magura Nappe (MN) and Podhale Flysch (PF) units. Measured velocities from GNSS observables are modelled by the auto-regressive integrated moving average (ARIMA) method to comprehend the long-term tectonic deformation. Our results showed that the ARIMA-modelled velocity varied from ~ 0.15 to ~ 8.86 mm/yr, indicating about 8.71 mm/yr difference along all units. Such differences suggest that crustal slip along the active thrusts and folds is the major factor causing regional deformation. The strain rates in PKB are also varying from the western to the eastern part. The rotational rates in PKB show a counterclockwise (CCW) pattern similar to the strain rates. These patterns suggesting that the PKB was rotated in the CCW direction with a large angle during the Miocene period. Finally, we analysed the seismicity for a period from 2004 to 2020 by using Bayesian moment tensor inversion and multivariate Bayesian inversion. The Bayesian inversion was applied based on bootstrapping chain analysis to figure out the earthquake mechanism using moment tensor inversion for the mainshock that occurred in Poland on 20 July 2018. The inversion results for the 2018 earthquake resolved a thrusting mechanism with nodal plane-1 having a strike of 346°, dip of 32°, and rake of 92° and a nodal plane-2 with a strike of 163°, dip of 58°, and rake of 89°. Since the seismicity in the Poland region has experienced less significant earthquakes in the last century, it is reasonable to attribute this lower seismic activity to the correspondingly low slip rates discerned through geodetic monitoring efforts.

Full-waveform reconstruction of micro-seismic events via topological derivative approach

Micro-seismic events, naturally occurring within geological formations and quasi-brittle engineered systems, provide a powerful window into the evolving processes of material degradation and failure. Accurate characterization of these events is critical for a comprehensive understanding of the underpinning fracturing mechanisms and potential implications. In this work, we present an algorithm for the spatial reconstruction and characterization of micro-seismic events in a three-dimensional bounded elastic body (with known geometry and nominal material properties) via combined source location and moment tensor inversion. Assuming availability of the full-waveform “acoustic emission” traces whose spectral content can be exposed via Fourier transform, the inverse solution is based on (i) a time-harmonic (forward) elastodynamic model and (ii) the concept of topological derivative as a framework for robust event reconstruction. On exploiting an equivalence between the elastic wavefield generated by the creation of a new micro-surface and that stemming from a suitable set of dipoles and double-couples (whose strengths are synthesized via the seismic moment tensor), we formulate the inverse problem as that for the real (in-phase) and imaginary (out-of-phase) components of the moment tensor at trial grid locations. In this way the optimal solution is obtained via a combinatorial search over a prescribed grid, inherently allowing for successive refinements of event reconstruction over the region(s) of interest. The analysis is illustrated by numerical experiments highlighting the key features of the inversion scheme including the reconstruction of multiple (i.e. contemporaneous) events, localization of the “off-grid” micro-seismic events, and the ability to handle noisy data. The results in particular highlight the utility of multi-frequency event reconstruction toward reducing the demand on the number of sensing locations.

Mechanism of overburden fracture induced earthquakes in coal seam mining

Mining-induced earthquakes are unnatural seismic events that frequently occur in high-position hard and thick rock strata during coal mining. Considering the frequent occurrence of strong mining-induced earthquakes in the Dongtan mining area, this study analysed the fracture migration characteristics of hard and thick rock strata and the focal mechanism of mining-induced earthquakes based on Volasov’s thick-plate and moment tensor inversion theories. The results showed that the main key strata were difficult to break under single-panel mining conditions because of the thick and high-strength rock strata and breakage of the main key strata is caused by multiple-panel mining. Volasov’s thick-plate theoretical calculation indicated an initial fracture span of the main key strata was 314 m, which is consistent with the actual mining distance of the working face. This verified that strong mining-induced earthquakes were induced by the initial fracture of the main key strata. In coal mining, the pure shear failure type of mining-induced earthquakes indicated the highest percentage, and the shear fracture of rock strata was the primary cause of strong mining-induced earthquakes. The dip angle of the focal fracture surface in mining-induced earthquakes was generally within 15°. Through an analysis of the focal mechanism of mining-induced earthquakes, it has a certain guiding role in explaining the mechanism of mining-induced earthquakes.

A new moment tensor catalogue of light to moderate earthquakes in Palu-Koro and Matano faults, Sulawesi, based on the regional network in Indonesia: Fault analysis and insight

Starting from 2009, substantial improvements in the regional seismic network operating in Sulawesi, as part of the Meteorological, Climatological, and Geophysical Agency (BMKG) station deployment, allowed us to determine earthquake mechanism for 140 shallow (5–60 km) light to moderate size (3.8 ≤ Mw ≤ 6.1) earthquakes between 2009 and 2018 along the Palu-Koro and Matano faults, through moment tensor inversion. Ten (10) out of the 140 solutions obtained in this study were found in the International Seismological Centre Bulletin and used as validation. The double-couple mechanism of the best solutions (variance reduction ≥0.4 and condition number ≤ 7.5) is further used to analyze the Palu-Koro and Matano segments. We observed that the Donggala segment of the Palu-Koro fault shows variable earthquake mechanism types. Whereas the two subsequent segments, the Palu and Saluki segments, show predominate oblique strike-slip with normal and reverse mechanisms. The Matano fault is characterized by strike-slip mechanisms, with normal mechanisms also observed. The normal mechanisms depict the pull-apart system in the Matano segment, manifested as Matano Lake. The results of this work improve our understanding of regional tectonics.

Joint Inversion of Regional Waveform, First-Motion Polarity, and Synthetic Aperture Radar Surface Displacement for the Fourth and Sixth North Korean Declared Nuclear Explosions

ABSTRACT This study analyzed the Democratic People’s Republic of Korea’s (DPRK) fourth (DPRK4, 6 January 2016 Mw 4.49) and sixth (DPRK6, 7 September 2017 Mw 5.2) declared nuclear tests, employing a joint seismic and Interferometric Synthetic Aperture Radar (InSAR) inversion to improve understanding of these events and enhance moment tensor (MT) inversion capabilities. The recent efforts have focused on employing seismic waveform and InSAR geodetic deformation data separately to analyze these and the previous nuclear tests (e.g., Chiang et al., 2018; Myers et al., 2018; Wang et al., 2018). Building upon our previous work (Chi-Durán et al., 2021), we performed a joint regional waveform, first-motion (FM) polarity, and surface displacement inversion, which demonstrated improved source-type discrimination, a revised MT solution with reduced scalar moment uncertainty, and an independently constrained location. In this article, we build on the previous results for DPRK6 by including an analysis using a four-layered velocity model with free-surface topography to compute the near-source static deformation Green’s functions. The model consists of a 50 m basalt layer (VP=2.07 km/s, VS=1.2 km/s), a 250 m stratified volcanic deposit layer (VP=1.73 km/s, VS=1.0 km/s), a 700 m weathered granodiorite layer (VP=2.5 km/s, VS=1.3 km/s), and a granodiorite half-space (VP=5.35 km/s, VS=3.09 km/s). The half-space shares the velocity of the regional MDJ2 velocity model (Ford et al., 2010), which has proven effective for waveform inversion in the region. This model considers the range of reported values for various lithologies and weathering effects. Our findings show that using the layered velocity model enhances the recovery of source location and depth for both the explosions by improving fits and reducing uncertainties. The joint inversion also improves source-type discrimination and better constrains the scalar seismic moment necessary for downstream yield estimation.

Seismicity and Focal Mechanism Analysis of The Local Earthquakes in the Zagros Fold-Thrust Belt of the Iraqi Kurdistan Region

During the period of 2020 to May 2022, the Iraqi Kurdistan region experienced several small to moderate size earthquakes. The seismic waveform data of these earthquakes were collected from regional and local broadband seismic stations. Waveform data was analyzed for seventy earthquakes to determine their source parameters, which may contribute to the active tectonics in those areas. Moment tensor inversion is estimated for the detected local earthquakes by analysis of observed seismic waveforms in different selected areas. The waveform inversion determines the seismic moment, depth, and kind of rupturing fault at the source. A grid search is employed as part of the inversion procedure across the strike, dip, and rake angles for each depth between 0.5 and 40 Km in increments of 1 Km. Waveform modeling is used to correlate the observed and created waveforms by applying bandpass filters and a suitable velocity model for Green's function using computer programs in seismology and the ISOLA code" Isolated asperities" for multiple-point source modeling. The results show that the earthquakes with M ≥ 4 have a focal depth ranging between 7 and 20 Km. Furthermore, the regional stresses were identified and the mechanism for the optimum fit was determined. The source mechanisms in the study area are most likely a combination of strike-slip, oblique-slip, and thrust.

Simulation of Shear and Tensile Fractures Using Ductile Phase Field Modelling with the Calibration of P Wave Velocity Measurement and Moment Tensor Inversion

The appropriate understanding and formulation of rock discontinuities via FEM is still challenging for rock engineering, as continuous algorithms cannot handle the discontinuities in rock mass. Also, different failure modes of rock samples, containing tensile and shear failure, need to be computed separately. In this study, a novel double-phase field damage model was introduced with two independent phase field damage variables. The construction of the proposed model follows the thermodynamics framework from the overall Helmholtz free energy, with elastic, plastic and surface damage components. The proposed model is calibrated via traditional damage variables, based on ultrasonic wave velocity measurement and acoustic emission monitoring, and both show great consistency between simulation results and laboratory observations. Then the double-phase field damage model is applied to COMSOL software to simulate microcrack propagation in a pre-fractured rock sample. Both lateral and wing cracks are observed in this study, manifested as shear- and tensile-dominated cracks. We also observed different microcracking mechanisms in the proposed numerical models, such as tensile and shear cracking, the influence of plastic strain and the percolation between tensile and shear microcracks. Overall, this study provides valuable insights into the mechanics of microcracking in rocks, and the proposed model shows promising results in simulating crack propagation.

Apparent Non-Double-Couple Components as Artifacts of Moment Tensor Inversion

Compilations of earthquake moment tensors from global and regional catalogs find pervasive non-double-couple (NDC) componentswith a mean deviation from a double-couple (DC) source of around 20%. Their distributions vary only slightly with magnitude, faulting mechanism, or geologic environments. This consistency suggests thatfor most earthquakes, especially smaller ones whose rupture processes are expected to be simpler, the NDC components are largely artifacts of the moment tensor inversion procedure. This possibility is also supported by the fact that NDC components for individual earthquakes with Mw<6.5 are only weakly correlated betweencatalogs. We explore this possibility by generating synthetic seismograms for the double-couple components of earthquakes around theworld using one Earth model and inverting them with a different Earth model. To match the waveforms with a different Earth model, the inversion changes the mechanisms to include a substantial NDC component while largely preserving the fault geometry (DC component). The resulting NDC components have a size and distribution similar to those reported for the earthquakes in the Global Centroid Moment Tensor (GCMT) catalog. The fact that numerical experiments replicate general features of the pervasive NDC components reported in moment tensor catalogs implies that these components are largely artifacts of the inversions not adequately accounting for the effects of laterally varying Earth structure.

Non‐Double‐Couple Components of Seismic Source: Method and Application to the 2014–2015 Bárðarbunga Volcanic Event Sequence, Iceland

AbstractGenuine non‐double‐couple (non‐DC) components of a seismic source, defined here as the non‐DC components that are not due to summation of pure double‐couple (DC) components, provide important insight into special physical processes in non‐earthquake sources such as explosion, volcano eruption and collapse etc. Yet they remain challenging to be resolved. To address the issue and explore the physical mechanism of those special events, we develop a waveform‐polarity‐based moment tensor (WPMT) inversion method and employ it to study physical process in the 2014–2015 Bárðarbunga volcano event sequence. The WPMT method incorporates P‐wave polarity data and seismic waveforms in the source inversion, designs a source simplicity test to check possible complex rupture in the seismic source, and employs a simulated annealing algorithm to search the best source solution. The simplicity test checks consistency of the source processes in the initiation stage of the event and the major energy release process of the event, thus ensuring that the inferred non‐DC components are genuine to the seismic source. Real event and synthetic tests indicate that the WPMT method can identify and resolve genuine non‐DC components in a seismic source. The WPMT inversions of the Bárðarbunga sequence yield many genuine non‐DC source components and reveal that the eruptions are accompanied by seismic activities in depths of 1–5 km with magma migrations out of chambers, collapses of conduits, failures of normal faults, and a magma recharge at a depth of 9 km accompanied by a failure on a nearby normal fault.

A Source and Ground Motion Study of the Veracruz Earthquakes of 29 October 2009 (Mw5.7) and 4 August 2021 (Mw4.8): Evidence of Strong Azimuthal Variation of Attenuation

We study two moderate earthquakes that occurred offshore the State of Veracruz, in the southwestern Gulf of Mexico, on 29 October 2009 (Mw 5.7) and 4 August 2021 (Mw 4.8). The former was located near the town of Alvarado and latter near the city of Veracruz. The events were well recorded by accelerographs and seismographs at local and regional distances. W-phase regional centroid moment tensor inversion reveals that they had reverse-faulting mechanism, similar to several other earthquakes in the southwestern Gulf of Mexico. Of the seven focal mechanisms now available along the southwestern margin, two are strike slip and the rest are of thrust type, suggesting a heterogeneous stress regime. We take advantage of local and regional recordings produced by these two earthquakes to study the characteristics of the ground motion. Source spectra computed at each station separately (without correcting for the site effect), assuming a reasonable geometrical spreading and Q = 141f 0.63, show remarkably high variability due to difference in path and local site effects. The geometric mean apparent source spectrum (source spectrum including site effects) of both earthquakes may be modeled by an ω2 -Brune source model with a stress drop, Δσ, of 40 MPa. These source spectra, along with the application of stochastic method, yield peak ground acceleration (PGA) and velocity (PGV) as a function of distance in general agreement with the observations. Of greater practical importance is the ground motion at sedimentary sites in the city of Veracruz and at the Laguna Verde Nuclear Power Plant (LVNPP) site, especially from a postulated Mw 6.5 earthquake which is a reasonable scenario event for the region. For the city of Veracruz and LVNPP we estimate site effect with respect to the ω2 -Brune source with Δσ of 2 MPa. There is some support for this Δσ. We apply both stochastic and empirical Green’s functions (EGF) techniques in the estimation of the ground motion. The recording of the 2021 earthquake is taken as the EGF, and Δσ of the EGF and the target event are assumed to be the same and equal to 2 MPa. The predicted PGA and PGV at sedimentary sites in the city of Veracruz and LVNPP above the hypocenter (depth = 20 km) from the postulated Mw 6.5 earthquake are 0.2 g and 10 cm/s and 0.18 g and 3 cm/s, respectively. These results are preliminary as they are based on several assumptions.

Mexico City Earthquake of 11 May 2023 (Mw3.2)

On 11 May 2023 a local earthquake in Mexico City was felt very strongly in Mixcoac, San Angel, and Coyoacán. The event was part of a seismic sequence that had begun about 6 months earlier. Peak Ground Acceleration (PGA) at the closest station (distance ~ 1 km), located in the hill zone, was ~ 0.18 g. Although the response spectrum at short periods at this station exceeded the design spectrum specified in the Mexico City´s Building Code, no structural damage was reported. Moment tensor inversion of bandpass filtered (0.08 – 0.24 Hz) displacement records yields M0 = 6.8x1013 N-m (Mw 3.2), H = 0.7 km, and the likely fault plane characterized by φ = 2700, δ = 760,On 11 May 2023 a local earthquake in Mexico City was felt very strongly in Mixcoac, San Angel, and Coyoacán. The event was part of a seismic sequence that had begun about 6 months earlier. Peak Ground Acceleration (PGA) at the closest station (distance ~ 1 km), located in the hill zone, was ~ 0.18 g. Although the response spectrum at short periods at this station exceeded the design spectrum specified in the Mexico City´s Building Code, no structural damage was reported. Moment tensor inversion of bandpass filtered (0.08 – 0.24 Hz) displacement records yields M0 = 6.8×1013 N-m (Mw 3.2), H = 0.7 km, and the likely fault plane characterized by φ = 2700, δ = 760, λ = -750. These source characteristics are very similar to those estimated for the 17 July 2019 earthquake which occurred during a swarm-like seismic activity about 5 km to the north. Spectral analysis of recordings at 19 sites in the hill zone, 14 in the transition zone, and 41 in the lake-bed zone reveals great variability of the ground motion within each of the zones. Estimated stress drop, Δσ, is 0.5 MPa. A large disparity is found between the observed source spectrum and theoretical source spectrum; their ratio provides an estimation of the amplification of seismic waves as they travel through the layers of decreasing velocity at shallower depth. We denote this ratio as the site effect. Predicted PGA and PGV for an Mw 3.2 earthquake, computed using stochastic technique (Boore 1983, 2003), assuming a Brune ω-2 source, Δσ = 0.5 MPa and including the site effect, are in reasonable agreement with the observations. Expected PGA and PGV at the epicenter of a postulated Mw 5 earthquake are 0.6 g and 60 cm/s at a generic hill-zone site; the expected values are twice as large in the lake-bed zone. These predictions should, however, be taken with caution as they are based on several approximations.

Complex rupture dynamics of the extremely shallow August 2020 M5.1 Sparta, North Carolina earthquake

On August 9, 2020, an Mw 5.1 earthquake ruptured the uppermost crust near the town of Sparta, North Carolina, creating the first co-seismic faulting surface rupture documented in the Eastern United States. Combining deep learning and matched filter earthquake detection, with differential-travel times relocation, we obtain a catalog of 1761 earthquakes, about 5.8 times the number of events listed in the standard USGS/NEIC catalog. The relocated seismicity revealed a complex fault structure with distinct planar alignments, supported by a moment tensor inversion with significant non-double-couple component. The Sparta mainshock with a centroid depth of 1.3 km is interpreted to have nucleated near the intersection of two main fault strands. The mainshock likely ruptured a blind strike-slip fault and a reverse fault associated with the identified surface rupture, both possibly part of a flower structure-like diffuse fault zone. Our observations highlight a complex behavior of extremely shallow earthquakes in stable continental regions.

Using Moment Tensor Inversions for Rapid Seismic Source Detection and Characterization: Application to the North Korean Nuclear Tests

AbstractThe rapid detection and source characterization of any type of seismic events including earthquakes and nuclear explosions is one of the missions of many seismological laboratories. Most often, the techniques used are based on phase picking and amplitude measurements for detecting and locating events, and for estimating magnitudes. From these parameters, especially for anthropogenic sources, event screening is then done empirically, and this may lead to misinterpretations of the source nature. However, it is known that seismic waveform inversion for the determination of the moment tensor has proven to be a reliable source physics-based method for event characterization. Here, we present a technique already used for earthquake monitoring in tectonically active regions, and we test it on very shallow natural and anthropogenic events recorded in the vicinity of the Democratic People’s Republic of Korea (DPRK). From a grid of potential locations and by scanning continuous seismic waveforms, it is possible to implement a rapid detector of seismic events providing the full information of the sources (origin time, location, magnitude, mechanism, and source decomposition). We show its overall performances on all past DPRK nuclear tests and regional earthquakes. From such an approach fast complete event screening is achieved. Source uncertainties can also be estimated. This stand-alone detector and identifier of seismic events may help monitoring seismological agencies to provide a rapid and complete alert for any events within a region of interest before more in-depth discrimination analysis can be run.

The Effects of Layer-Induced Elastic Anisotropy on Microseismic Monitoring Records in Underground Coal Mines

Real-time microseismic observations are extensively used in underground hard rock mines. However, microseismic monitoring in coal mines is not common. In this study, the features of microseismic (MS) event waveforms were discussed under conditions of thin-layered sedimentary media. The possibility of source mechanism evaluation was considered for the MS event recorded in coal mine. Moment tensor (MT) inversion of the waveforms can be challenging due to influence of thin-layered medium structure. Splitting of shear SH- and SV-waves is one of the characteristic features that affect MS records in coal mines. The effect of splitting was demonstrated using real microseismic data collected by an underground seismic array installed in the Polosukhinskaya coal mine, Russia. Velocities of shear SH- and SV-waves were estimated from travel times. The results reveal significant variations in shear wave velocities composed of a radial anisotropy of 75%, thus confirming the importance of a thin-layered model for moment tensor inversion. The thin-layered model was built for geological conditions at the Polosukhinskaya coal mine and synthetic seismograms were generated using QSEIS (Pyrocko OS). A full moment tensor of a single typical event was found employing the layered 1D velocity model.

Improved numerical solution of anisotropic poroelastic wave equation in microseismicity: Graphic process unit acceleration and moment tensor implementation

AbstractThe accuracy and computational efficiency of full waveform forward modelling in poroelastic media are crucial for microseismic monitoring. It enables intuitive, precise and efficient simulation of subsurface responses, thereby improving the reliability of moment tensor inversion and seismic source mechanism interpretation. Additionally, it reflects the role of fluid effects in waveform evolution. In this paper, based on the Biot mechanism, we derived the first‐order velocity–stress equation of poroelastic media and discretized the model using a rotated staggered grid. The rotated staggered grid can well adapt to anisotropic media with high contrast parameters. We provide the finite difference formula based on graphic process unit–acceleration and moment tensor and also provide the graphic process unit workflow for forward modelling of anisotropic poroelastic media. First, two models with different grid sizes were run based on single graphic process unit, 1‐thread Central Processing Unit (CPU) and 16‐thread CPU. The results show that the speedup factors are approximately 14.3 and 3.7, respectively, compared with the running time of 1‐thread CPU and 16‐thread CPU. Then, we compare and evaluate the response of three basic source mechanisms (isotropic expansion, double couple and compensated linear vector dipole) in the model. The comparison of analytical and numerical results verifies the effectiveness of the method. Furthermore, wavefield snapshots of two typical anisotropic media (vertical transversely isotropic and horizontal transversely isotropic) are analysed to correspond to different moment tensor sources. The results showed that the source mechanism does not change the isotropic and anisotropic plane and the wave travel time, but it does change the polarization amplitude of the velocity component. The attenuation of slow qP‐wave increases along with the increase of the value of viscosity. The effect of permeability on wavefield appears with the opposite effect of viscosity. Finally, the seismic waveform differences between multi‐layer elastic media and poroelastic media are compared and analysed. The results showed that the seismic wavefield and waveform of poroelastic media are more complex, the propagation speed of seismic waves is faster, but the attenuation of seismic waves is stronger.

Empirical analysis of unloading creep damage mechanism of salt rock at 110 °C

Underground energy storage in salt caverns involves operating in high-temperature conditions within deep strata, leading to unloading creep in the surrounding salt rock during cavity construction and natural gas extraction. Therefore, understanding the deformation and damage evolution of salt rock in such conditions is vital for ensuring safe construction and prolonged operation of deep salt cavern gas storage. The present study involves graded unloading creep tests and synchronous acoustic emission tests on the salt rock at 110 °C. The moment tensor inversion theory is employed to analyze focal mechanisms and microscopic crack propagation modes based on the test results. The findings reveal that under constant deviatoric stress, decreasing confining pressure exponentially increases the duration of the deceleration creep stage in salt rock. The creep increment and axial steady-state creep rate exhibited a U-shape trend. Instantaneous strain during stress change mainly involves elastic-plastic strain, where the instantaneous strain decreases first and then increases, exhibiting hardening characteristics in the middle process. The evolution of microcracks in salt rock can be categorized into three stages: compaction stage, stable propagation stage, and accelerated propagation stage. The acoustic emission energy, b-value, and fractal dimension of each stage are different. In unloading creep tests, the peak frequency in acoustic emission signal presents zonal distribution in low-frequency and high-frequency regions. It is found that as the confining pressure decreases, the peak frequency of acoustic emission signal and the proportion of high-frequency signal decrease, and the b-value and fractal dimension exhibit an upward-fluctuation-downward trend. The rupture slip angle of the acoustic emission source can effectively characterize the creep slip behavior in salt rock crystals. The slip gradually shifts from vertical to horizontal, while the strike angle and dip angle of rupture orientation are minimally affected by unloading. The obtained results can provide a reference for the safe construction and the design of the operation pressure for ultra-deep salt cavern gas storage.