Search NASA⌕ Search

SEARCH · Search NASA

Results for “microearthquakes”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Machine Learning‐Assisted Microearthquake Location Workflow for Monitoring the Newberry Enhanced Geothermal System

Abstract Enhanced geothermal systems (EGS) offer a sustainable energy source but face challenges in accurately locating microearthquakes induced during reservoir stimulation. Locating these microearthquakes provides reliable feedback on the stimulation progress. Current deep learning methods for locating earthquakes require extensive data sets for training, which is problematic as detected microearthquakes are often limited. To address the scarcity of training data, we propose a practical workflow using probabilistic multilayer perceptron (PMLP) which predicts microearthquake locations from cross‐correlation time lags in waveforms. Utilizing a 3D velocity model of Newberry site derived from ambient noise interferometry, we generate numerous synthetic microearthquakes and 3D acoustic waveforms for PMLP training. Accurate synthetic tests prompt us to apply the trained network to the 2012 and 2014 stimulation field waveforms. To enhance the accuracy of source localization, we carefully handpick the P‐arrival times. Predictions on the 2012 stimulation data set show major microseismic activity at depths of 0.5–1.2 km, correlating with a known casing leakage scenario. In the 2014 data set, the majority of predictions concentrate at 2.0–2.9 km depths, consistent with results obtained from conventional physics‐based inversion, and align with the presence of natural fractures from 2.0 to 2.7 km. We validate our findings by comparing the synthetic and field picks, demonstrating a satisfactory match for the first arrivals. By combining the benefits of quick inference speeds and accurate location predictions, we demonstrate the feasibility of using realistic synthetic data set to locate microseismicity for EGS monitoring.

15 GEOTHERMAL ENERGY↗

Deep learning forecasts the spatiotemporal evolution of fluid-induced microearthquakes

Microearthquakes generated by subsurface fluid injection record the evolving stress state and permeability of reservoirs. Forecasting their spatiotemporal evolution is therefore critical for applications such as enhanced geothermal systems, carbon dioxide sequestration and other geoengineering applications. Here we propose a transformer neural network model that ingests hydraulic stimulation history and prior microearthquake observations to forecast four key quantities: cumulative microearthquake count, cumulative logarithmic seismic moment, and the 50th- and 95th-percentile extents of the microearthquake cloud. Applied to the EGS Collab Experiment 1 dataset, the model achieves R2 > 0.98 for the 1-s forecast horizon and R2 > 0.88 for the 15-s forecast horizon across all targets, and supplies uncertainty estimates through a learned standard deviation term. These accurate, uncertainty-quantified forecasts enable real-time inference of fracture propagation and permeability evolution, demonstrating the strong potential of deep-learning approaches to improve seismic-risk assessment and guide mitigation strategies in future fluid-injection operations.

Chung, Jaehong↗

Data Arrays for Microearthquake (MEQ) Monitoring using Deep Learning for the Newberry EGS Sites

The 'Machine Learning Approaches to Predicting Induced Seismicity and Imaging Geothermal Reservoir Properties' project looks to apply machine learning (ML) methods to Microearthquake (MEQ) data for imaging geothermal reservoir properties and forecasting seismic events, in order to advance geothermal exploration and safe geothermal energy production. As part of the project, this submission provides data arrays for 149 microearthquakes between the year 2012 and 2013 at the Newberry EGS Site for use with the Deep Learning Algorithm that has been developed. The data provided includes raw waveform data, location data, normalized waveform data, and processed waveform data. Penn State Geothermal Team has shared the following files from the project: - 149 microearthquakes (MEQs) between 2012 and 2013 at Newberry EGS sites, 'Normalized Waveform Inputs.npz' are normalized waveforms. - labels of 149 MEQs: Processed Waveform Inputs.npz - location labels of 149 MEQs: Location Data.npz Note: .npz is the python file format by NumPy that provides storage of array data.

15 GEOTHERMAL ENERGY↗

Seismic monitoring of EGS fracture stimulations at Utah FORGE (Part 1): Time-lapse variations of b-values and Shear-wave Splitting Rates of induced microearthquakes

Characterizing geologic formations near geothermal wells is crucial for efficient extraction of subsurface geothermal energy resources. To monitor effective fracture stimulation of enhanced geothermal systems (EGS), we analyze a catalog of induced microearthquakes to characterize seismicity induced by Phase 2C fracture stimulations of the Utah FORGE project and derive the b-values for 424 induced microearthquakes in three stages of fracture stimulation. Here, we find the time-lapse variations of these b-values. We apply eigenvector rotation and waveform correlation to 3C geophone data to separate fast S 1 and slow S 2 waves, determine their directions, measure their lag time, and calculate the S-wave splitting rate (SSR) for six induced microearthquakes in Stage 1 and 21 in Stage 3. Both SSR values and b-values increase with time. The increase in these values may be related to heavily fractured zones created by fracture stimulations during Stage 2 injections. Large and highly variable SSR values in Stage 3 may reflect heterogeneities in the fracture zones. Time-lapse variations of SSR values and b-values may be one of the most useful approaches to monitoring and characterizing the effectiveness and efficiency of various types of fracture stimulations.

58 GEOSCIENCES↗

Crustal permeability generated through microearthquakes is constrained by seismic moment

We link changes in crustal permeability to informative features of microearthquakes (MEQs) using two field hydraulic stimulation experiments where both MEQs and permeability evolution are recorded simultaneously. The Bidirectional Long Short-Term Memory (Bi-LSTM) model effectively predicts permeability evolution and ultimate permeability increase. Our findings confirm the form of key features linking the MEQs to permeability, offering mechanistically consistent interpretations of this association. Transfer learning correctly predicts permeability evolution of one experiment from a model trained on an alternate dataset and locale, which further reinforces the innate interdependency of permeability-to-seismicity. Models representing permeability evolution on reactivated fractures in both shear and tension suggest scaling relationships in which changes in permeability (Δk) are linearly related to the seismic moment (M) of individual MEQs as Δk ∝ M. This scaling relation rationalizes our observation of the permeability-to-seismicity linkage, contributes to its predictive robustness and accentuates its potential in characterizing crustal permeability evolution using MEQs.

15 GEOTHERMAL ENERGY↗

Quantifying rupture characteristics of microearthquakes in the Parkfield Area using a high-resolution borehole network

It is well known that large earthquakes often exhibit significant rupture complexity such as well separated subevents. With improved recording and data processing techniques, small earthquakes have been found to exhibit rupture complexity as well. Studying these small earthquakes offers the opportunity to better understand the possible causes of rupture complexities. Specifically, if they are random or are related to fault properties. We examine microearthquakes (M < 3) in the Parkfield, California, area that are recorded by a high-resolution borehole network. We quantify earthquake complexity by the deviation of source time functions and source spectra from simple circular (omega-square) source models. We establish thresholds to declare complexity, and find that it can be detected in earthquakes larger than magnitude 2, with the best resolution above M2.5. Comparison between the two approaches reveals good agreement (>90 per cent), implying both methods are characterizing the same source complexity. For the two methods, 60–80 per cent (M 2.6–3) of the resolved events are complex depending on the method. The complex events we observe tend to cluster in areas of previously identified structural complexity; a larger fraction of the earthquakes exhibit complexity in the days following the M w 6 2004 Parkfield earthquake. Ignoring the complexity of these small events can introduce artefacts or add uncertainty to stress drop measurements. Focusing only on simple events however could lead to systematic bias, scaling artefacts and the lack of measurements of stress in structurally complex regions.

58 GEOSCIENCES↗

Induced microearthquakes predict permeability creation in the brittle crust

Predicting the evolution of permeability accurately during stimulation at the reservoir scale and at the resolution of individual fractures is essential to characterize the fluid transport and the reactive/heat-transfer characteristics of reservoirs where stress exerts significant control. Here, we develop a hybrid machine learning (ML) model to visualize in situ permeability evolution for an intermediate-scale (~10 m) hydraulic stimulation experiment. This model includes an ML model that was trained using the well history of flow rate and wellhead pressure and MEQ data from the first three stimulation episodes to predict average permeability from the statistical features of the MEQs alone for later episodes. Moreover, a physics-inspired model is integrated to estimate in situ fracture permeability spatially. This method relates fracture permeability to fracture dilation and scales dilation to the equivalent MEQ magnitude, according to laboratory observations. The seismic data are then applied to define incremental changes in permeability in both space and time. Our results confirm the excellent agreement between the ground truth and model-predicted permeability evolution. The resulting permeability map defines and quantifies flow paths in the reservoir with the averaged permeability comparing favorably with the ground truth of permeability.

induced seismicity↗

Tectonic implications of the microearthquake seismicity and fault plane solutions in southern Peru

Because the contortion in the seismic zone in southern Peru is aligned approximately parallel to the direction of relative plate motion, rather than perpendicular to the coast of Peru, the position of the contortion need not migrate with respect to the overriding South American plate as the Nazca plate subducts beneath it, and the flow in the surrounding asthenosphere could be in a steady state. In addition, the position of the contortion defines the northern boundary of the volcanic arc in southern Peru. The inference that a wedge of asthenospheric material must overlie the downgoing slab for subduction-related volcanism to occur is thereby strengthened.

Grange, F.↗

Tectonic deformation of the Andes and the configuration of the subducted slab in central Peru: Results from a micro-seismic experiment

The vast majority of the microearthquakes recorded occurred to the east: on the Huaytapallana fault in the Eastern Cordillera or in the western margin of the sub-Andes. The sub-Andes appear to be the physiographic province subjected to the most intense seismic deformation. Focal depths for the crustal events here are as deep as 50 km, and the fault plane solutions, show thrust faulting on steep planes oriented roughly north-south. The Huaytapallana fault in the Cordillera Oriental also shows relatively high seismicity along a northeast-southwest trend that agrees with the fault scarp and the east dipping nodal plane of two large earthquakes that occurred on this fault in 1969. The recorded microearthquakes of intermediate depth show a flat seismic zone about 25 km thick at a depth of about 100 km. This agrees with the suggestion that beneath Peru the slab first dips at an angle of 30 deg to a depth of 100 km and then flattens following a quasi-horizontal trajectory. Fault plane solutions of intermediate depth microearthquakes have horizontal T axes oriented east-west.

Suarez, G.↗

Induced seismicity during the 2012 Newberry EGS stimulation: Assessment of two advanced earthquake detection techniques at an EGS site

To improve the microearthquake detection capability during the 2012 Newberry enhanced geothermal system (EGS) demonstration, we applied two advanced earthquake detection algorithms to the continuous seismic data. These techniques have the power to detect microearthquakes within the seismic noise that would otherwise be missed using traditional techniques. We compare the efficacy of a matched field processing algorithm with a multi-channel waveform correlation algorithm. Both methods identified between 60%–80% more events, with each slightly outperforming the other in different seismic noise contexts. Incorporating either template-matching technique into a seismic processing routine could therefore effectively increase the sensitivity of the seismic network without changes to the existing hardware. Additionally, lowering the magnitude of completeness would significantly add to the number of events that could subsequently inform future decision making at the reservoir level.

15 GEOTHERMAL ENERGY↗

Machine Learning Approaches to Predicting Induced Seismicity and Imaging Geothermal Reservoir Properties

This project developed machine learning (ML) methods, lab data sets, and field data to advance geothermal exploration and geothermal energy production. The work had three focus areas. One involved the development of ML methods to use microearthquakes (MEQs) for imaging geothermal reservoir properties and improving subsurface characterization – most importantly the evolution of permeability within the evolving reservoir. This part of the work included development of ML approaches for automated MEQ location, focal mechanism determination and identification of earthquake precursors. The second area focused on using MEQ signals generated by geothermal exploration and production to predict the relationship between fluid injection and seismicity. Here, we extended to reservoir scale our success in using ML to predict laboratory earthquakes and fault zone stress state. The third focus area was on lab experiments. Here, we developed new ML models for lab earthquake prediction and identification of precursors to failure to improve earthquake forecasting and early warning in geothermal settings. Major outcomes of our work include ML models that learn from MEQ signals during geothermal exploration and production to predict induced seismicity. MEQs occur naturally in connection with drilling and energy production. We developed ML methods to use the seismic waves from these events to characterize the elastic, hydraulic and poromechanical properties of reservoirs. Our work illuminated fracture geometry and the evolution of fracture permeability by incorporating seismic coda wave analysis and ML methods to relate fluid injection and seismicity. We significantly expanded laboratory earthquake prediction to include methods that use both passive measurements of microearthquakes within the lab fault zones and also active source acoustic measurements of fault zone elastic properties. These methods can now predict fault zone stress state, time to failure and the magnitude of lab earthquakes. Our work showed that repetitive stick- slip failure events during frictional sliding (the lab equivalent of earthquakes) are preceded by a cascade of micro-failure events that radiate energy in a manner that foretells unstable failure – manifest as laboratory MEQs. We documented a mapping between fracture properties and statistical attributes of elastic radiation. We extended existing works to geothermal reservoir scale and developed ML methods to determine reservoir permeability, fracture properties, and their evolution during geothermal energy production. An attractive feature of ML algorithms is their ability to handle big datasets and reveal patterns and correlations that may remain invisible to conventional analyses. Our work connected data from field, laboratory and intermediate scales to study permeability, stress, strength, fracture stiffness and geometry. At the field scale we used data from the Newberry Volcano field site, UtahFORGE, EGS Collab, and also the Bedretto underground research lab in Switzerland. These data sets are bridging the gap between the lab scale, theory, and reservoir scale. Our work produced plain language summaries to improve public understanding of DOE research. We also developed openly distributed ML and seismicity datasets for use by all researchers and we published connections between induced seismicity in geothermal areas and reservoir properties including permeability, fracture properties, and stress state. Our models are designed for the large data sets of induced seismicity typically associated with geothermal sites. We produced labeled event catalogs and used them on geothermal data to assess how ML can facilitate geothermal production and exploration. All datasets are available on the GDR Productivity: The project produced 32 publications in peer reviewed journals (two are in review). It supported the work of 6 PhD students, 40 conference presentations, 6 keynote talks at national meetings, and mentoring and professional development for 4 postdoctoral fellows.

15 GEOTHERMAL ENERGY↗

Strain buildup and release, earthquake prediction and selection of VLB sites for the margins of the North Pacific

Microearthquake studies of the southern section of the San Andreas fault system were conducted in order to elucidate the role of known but little studied complications in the fault system which could affect the SAFE measurement. Data reduction from microearthquake studies in Baja California, Mexico, was completed. In addition, styles of deformation on the Alpine fault in New Zealand and on the San Andreas fault were compared. Sections of these faults with comparable physical characteristics seem to deform in a similar manner, indicating that deformation style is controlled by certain fundamental relations and will continue in the future. The geology of eastern New Guinea is also discussed. Deformation in the area is proposed to be related to collision of a series of island arcs with central New Guinea, which is part of the Australian continent. Among other unusual properties, eastern New Guinea contains the highest anomaly of the gravimetrically determined geoid based on GEM 6 and surface gravity data.

Scholz, C. H.↗

Complementary Subsurface Characterization Methods to Develop a Geologic Model for the EGS Collab Experiment, Sanford Underground Research Facility

The EGS (Enhanced Geothermal Systems) Collab project was performed within the Sanford Underground Research Facility (SURF) with a goal of understanding processes and evaluation of models related to hydraulic stimulation of rock at depth. The present work deals with the development of Testbed 2 where experiments were conducted at a depth of 1.25 km and were located within a well-characterized testbed in a metamorphic, amphibolite host rock. A total of eleven boreholes varying in length between 10.6 m and 81.2 m were continuously cored to develop the testbed. In addition to the continuous coring of the amphibolite host rock, geophysical instrumentation supporting electrical resistivity tomography (ERT), microearthquake (MEQ) detection, and optical fiber providing distributed temperature sensing (DTS) and distributed acoustic sensing (DAS) were installed in the monitoring boreholes, all of which produced a comprehensive complementary suite of characterization and monitoring technologies. Pre-stimulation characterization of the groundwater conditions identified only two hydraulically significant fractures, neither of which transects the central portion of the testbed. Flow and pressure monitoring indicated low preexisting pore pressure conditions likely affected by the mine openings.

15 GEOTHERMAL ENERGY↗

Complementary Subsurface Characterization Methods to Develop a Geologic Model for the EGS Collab Experiment, Sanford Underground Research Facility

The EGS (Enhanced Geothermal Systems) Collab project was performed within the Sanford Underground Research Facility (SURF) with a goal of understanding processes and evaluation of models related to hydraulic simulation of rock at depth. The present work deals with the development of Testbed 2 where experiments were conducted at a depth of 1.25 km and were located within a well-characterized testbed in a metamorphic, amphibolite host rock. A total of eleven boreholes varying in length between 10.6 m and 81.2 m were continuously cored to develop the testbed. In addition to the continuous coring of the amphibolite host rock, geophysical instrumentation supporing electrical resisivity tomography (ERT), microearthquake (MEQ) detection, and optical fiber providing distributed temperature sensing (DTS) and distributed acoustic sensing (DAS) were installed in the monitoring boreholes, all of which produced a comprehensive complementary suite of characterization and monitoring technologies.

15 GEOTHERMAL ENERGY↗

Three-dimensional paganica fault morphology obtained from hypocenter clustering (L'Aquila 2009 seismic sequence, Central Italy)

In seismic modelling, fault planes are normally assumed to be flat due to the lack of data which can constrain fault morphology. However, incorporating 3D fault morphology is important for modelling several phenomena, for example calculating mainshock induced stress changes. Here we utilize a data-analytical method to unveil the 3D rupture morphology of faults using unsupervised clustering techniques applied to earthquake hypocenters in seismic sequences. We apply this method to the 2009 L'Aquila seismic sequence which involved a M W 6.1 mainshock on April 6th. We use a dataset of about 50,000 relocated events, mostly microearthquakes, reaching magnitude of completeness equal to 0.7. Clustering distinguishes the earthquakes as occurring in three main clusters along with other minor fault segments. We then represent the morphology of the main Paganica fault system (responsible for the largest mainshock) using splines. This method shows promise as a step toward robustly and quickly obtaining 3D rupture morphologies where earthquake sequences have been monitored. The 3D model is presented interactively online, and the processing is presented in an interactive Jupyter Notebook (https://bit.ly/2MnCFdj).

58 GEOSCIENCES↗

Dynamics, Radiation, and Overall Energy Budget of Earthquake Rupture With Coseismic Off-Fault Damage

Earthquake ruptures dynamically activate coseismic off–fault damage around fault cores. Systematic field observation efforts have shown the distribution of off–fault damage around main faults, while numerical modeling using elastic–plastic off–fault material models have demonstrated the evolution of coseismic off–fault damage during earthquake ruptures. Laboratory–scale microearthquake experiments have pointed out the enhanced high–frequency radiation due to the coseismic off–fault damage. However, the detailed off–fault fracturing mechanisms, subsequent radiation, and its contribution to the overall energy budget remain to be fully understood because of limitations of current observational techniques and model formulations. Here, we constructed a new physics–based dynamic earthquake rupture modeling framework, based on the combined finite–discrete element method, to investigate the fundamental mechanisms of coseismic off–fault damage, and its effect on the rupture dynamics, the radiation and the overall energy budget. We conducted a 2–D systematic case study with depth and showed the mechanisms of dynamic activation of the coseismic off–fault damage. We found the decrease in rupture velocity and the enhanced high–frequency radiation in near field due to the coseismic off–fault damage. We then evaluated the overall energy budget, which shows a significant contribution of the coseismic off–fault damage to the overall energy budget even at depth, where the damage zone width becomes narrower. Here, the present numerical framework for the dynamic earthquake rupture modeling thus provides new insights into earthquake rupture dynamics with the coseismic off–fault damage.

58 GEOSCIENCES↗

Identifying Different Classes of Seismic Noise Signals Using Unsupervised Learning

Abstract Proper classification of nontectonic seismic signals is critical for detecting microearthquakes and developing an improved understanding of ongoing weak ground motions. We use unsupervised machine learning to label five classes of nonstationary seismic noise common in continuous waveforms. Temporal and spectral features describing the data are clustered to identify separable types of emergent and impulsive waveforms. The trained clustering model is used to classify every 1 s of continuous seismic records from a dense seismic array with 10–30 m station spacing. We show that dominate noise signals can be highly localized and vary on length scales of hundreds of meters. The methodology demonstrates the complexity of weak ground motions and improves the standard of analyzing seismic waveforms with a low signal‐to‐noise ratio. Application of this technique will improve the ability to detect genuine microseismic events in noisy environments where seismic sensors record earthquake‐like signals originating from nontectonic sources.

Johnson, Christopher W.↗

Preliminary Analysis of Source Physics Experiment Explosion-Triggered Microseismicity Using the Back-Projection Method

A series of four chemical explosions were detonated in a deep borehole within the Yucca Flat Dry Alluvium Geology (DAG) at the Nevada National Security Site between 2018 and 2019. The two larger chemical explosions of 50 tons (DAG-2) and 10 tons (DAG-4) TNT equivalent yield triggered energetic aftershock sequences numbering 1392 and 347 microearthquakes, respectively, within the first 10 days. No significant aftershock activity was observed for the two smaller 1-ton explosions (DAG-1 and DAG-3). Here, we used a back-projection method based on travel-time migration and stacking of signal-to-noise ratio traces to detect, associate and locate aftershocks from a subset of 22-geophones within a larger 2 urn:x-wiley:21699313:media:jgrb54910:jgrb54910-math-0001 2 km seismic array surrounding the borehole. The aftershocks located within 300 m of the borehole and the depths were above the working points of 300 and 50 m depths of DAG-2 and DAG-4, respectively, ruling out triggering slip on geologic faults or disturbances beneath neighboring collapse craters. DAG-2 and DAG-4 aftershocks decayed at similar rates, with power-law exponents of p = 1.48 and p = 1.49, respectively. These decay rates are comparable to aftershocks sequences triggered by earthquakes and historical nuclear explosions at Yucca Flat. A smooth power-law aftershock decay within the first 10 days suggests a triggering mechanism from explosion generated stress relaxation due to the diffusion of high gas pressures in the cavity and radial fractures. A more random and episodic aftershock rate would be expected due to cavity collapse or falling rubble in chimney formation.

58 GEOSCIENCES↗