Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismic data processing”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

2D Seismic Reflection Data across Central Illinois

In a continuing collaboration with the Midwest Geologic Sequestration Consortium (MGSC) on the Evaluation of the Carbon Sequestration Potential of the Cambro-Ordovician Strata of the Illinois and Michigan Basins project, Schlumberger Carbon Services and WesternGeco acquired two-dimensional (2D) seismic data in the Illinois Basin. This work included the design, acquisition and processing of approximately 125 miles of (2D) seismic reflection surveys running west to east in the central Illinois Basin. Schlumberger Carbon Services and WesternGeco oversaw the management of the field operations (including a pre-shoot planning, mobilization, acquisition and de-mobilization of the field personnel and equipment), procurement of the necessary permits to conduct the survey, post-shoot closure, processing of the raw data, and provided expert consultation as needed in the interpretation of the delivered product. Three 2D seismic lines were acquired across central Illinois during November and December 2010 and January 2011. Traversing the Illinois Basin, this 2D seismic survey was designed to image the stratigraphy of the Cambro-Ordovician sections and also to discern the basement topography. Prior to this survey, there were no regionally extensive 2D seismic data spanning this section of the Illinois Basin. Between the NW side of Morgan County and northwestern border of Douglas County, these seismic lines ran through very rural portions of the state. Starting in Morgan County, Line 101 was the longest at 93 miles in length and ended NE of Decatur, Illinois. Line 501 ran W-E from the Illinois Basin – Decatur Project (IBDP) site to northwestern Douglas County and was 25 miles in length. Line 601 was the shortest and ran N-S past the IBDP site and connected lines 101 and 501. All three lines are correlated to well logs at the IBDP site. Originally processed in 2011, the 2D seismic profiles exhibited a degradation of signal quality below ~400 millisecond (ms) which made interpretation of the Mt. Simon and Knox sections difficult. The data quality also gradually decreased moving westward across the state. To meet evolving project objectives, in 2012 the seismic data was re-processed using different techniques to enhance the signal quality thereby rendering a more coherent seismic profile for interpreters. It is believed that the seismic degradation could be caused by shallow natural gas deposits and Quaternary sediments (which include abandoned river and stream channels, former ponds, and swamps with peat deposits) that may have complicated or changed the seismic wavelet. Where previously limited by seismic coverage, the seismic profiles have provided valuable subsurface information across central Illinois. Some of the interpretations based on this survey included, but are not limited to: - Stratigraphy generally gently dips to the east from Morgan to Douglas County. - The Knox Supergroup roughly maintains its thickness. There is little evidence for faulting in the Knox. However, at least one resolvable fault penetrates the entire Knox section. - The Eau Claire Formation, the primary seal for the Mt. Simon Sandstone, appears to be continuous across the entire seismic profile. - The Mt. Simon Sandstone thins towards the western edge of the basin. As a result, the highly porous lowermost Mt. Simon section is absent in the western part of the state. - Overall basement dip is from west to east. - Basement topography shows evidence of basement highs with on-lapping patterns by Mt. Simon sediments. - There is evidence of faults within the lower Mt. Simon Sandstone and basement rock that are contemporaneous with Mt. Simon Sandstone deposition. These faults are not active and do not penetrate the Eau Claire Shale. It is believed that these faults are associated with a possible failed rifting event 750 to 560 million years ago during the breakup of the supercontinent Rodinia.

Smith, Valerie↗
DOE OSTI2015

IDC Reengineering Phase 2 Project Scope

The Comprehensive Nuclear-Test-Ban Treaty Organization (CTBTO) operates the International Data Centre (IDC) to support treaty verification activities by Member States. The IDC collects, stores, processes, and distributes data from seismic, hydroacoustic, infrasound, and radionuclide stations. A significant component of the IDC is the seismic, hydroacoustic, and infrasound (SHI) processing system. The IDC has recognized the need to reengineer the SHI processing system and has defined a three-phase reengineering effort. IDC Reengineering Phase 1 (RP1) started in 2011 with the goal to significantly modernize SHI automatic and interactive data processing software, focusing on enhancements to individual components of the system. Reengineering Phase 2 (RP2), started in 2014, addresses the task of specifying a unified architecture for all SHI software, across processing stages, to pave the way for further software development and sustainment in the future. A key objective of RP2 is an architecture sufficient to provide a basis for a cost estimate for the development or enhancement of the software components and subsystems. This report provides a brief description of the scope of the IDC Reengineering Phase 2 project. It describes the goals and objectives of reengineering, the system definition, and the technical scope of the system.This report expands upon the “Scope for IDC Reengineering Phase 2”, July 2014, ECS-DIS-WGB-43-PTS-MATERIAL-12-ADD-1, produced by the IDC.

58 GEOSCIENCES↗
DOE OSTI2019

Frequency-Bessel Transform Method for Effective Imaging of Higher-Mode Rayleigh Dispersion Curves From Ambient Seismic Noise Data

It has been widely recognized that the cross-correlation function of ambient seismic noise data recorded at two stations approximates to the part of Greens function between two stations. Therefore, the cross-correlation function should include higher modes, aside from the fundamental mode. However, the problem of measuring or extracting overtones from ambient seismic noise data remains. In this paper, we propose the frequency-Bessel transform method (F-J method) for extracting the dispersion curves of higher modes from ambient seismic noise data. We then assess the validity, accuracy, and applicability of the F-J method by conducting extensive numerical simulations and processing the observed ambient seismic noise data of the USArray. As demonstrated in this study, the F-J method is a convenient, practical, and accurate method for extracting the dispersion curves of multimodes from ambient seismic noise data and therefore has significant potentiality in the field of ambient seismic noise tomography.

58 GEOSCIENCES↗
DOE OSTI2021

Machine Learning Predicts the Timing and Shear Stress Evolution of Lab Earthquakes Using Active Seismic Monitoring of Fault Zone Processes

Abstract Machine learning (ML) techniques have become increasingly important in seismology and earthquake science. Lab‐based studies have used acoustic emission data to predict time‐to‐failure and stress state, and in a few cases, the same approach has been used for field data. However, the underlying physical mechanisms that allow lab earthquake prediction and seismic forecasting remain poorly resolved. Here, we address this knowledge gap by coupling active‐source seismic data, which probe asperity‐scale processes, with ML methods. We show that elastic waves passing through the lab fault zone contain information that can predict the full spectrum of labquakes from slow slip instabilities to highly aperiodic events. The ML methods utilize systematic changes in P‐wave amplitude and velocity to accurately predict the timing and shear stress during labquakes. The ML predictions improve in accuracy closer to fault failure, demonstrating that the predictive power of the ultrasonic signals improves as the fault approaches failure. Our results demonstrate that the relationship between the ultrasonic parameters and fault slip rate, and in turn, the systematically evolving real area of contact and asperity stiffness allow the gradient boosting algorithm to “learn” about the state of the fault and its proximity to failure. Broadly, our results demonstrate the utility of physics‐informed ML in forecasting the imminence of fault slip at the laboratory scale, which may have important implications for earthquake mechanics in nature.

58 GEOSCIENCES↗
DOE OSTI2013

Rotation-Enabled 7-Degree of Freedom Seismometer for Geothermal Resource Development. Phase 1 Final Report

Under this Department of Energy (DOE) grant, A-Tech Corporation d.b.a. Applied Technology Associates (ATA), seeks to develop a seven-degree-of-freedom (7-DOF) seismic measurement tool for high-temperature geothermal applications. The Rotational-Enabled 7-DOF Seismometer includes a conventional tri-axial accelerometer, a conventional pressure sensor or hydrophone, and a tri-axial rotational sensor. The rotational sensing capability is novel, based upon ATA's innovative research in rotational sensing technologies. The geothermal industry requires tools for high-precision seismic monitoring of crack formation associated with Enhanced Geothermal System (EGS) stimulation activity. Currently, microseismic monitoring is conducted by deploying many seismic tools at different depth levels along a 'string' within drilled observation wells. Costs per string can be hundreds of thousands of dollars. Processing data from the spatial arrays of linear seismometers allows back-projection of seismic wave states. In contrast, a Rotational-Enabled 7-DOF Seismometer would simultaneously measure p-wave velocity, s-wave velocity, and incident seismic wave direction all from a single point measurement. In addition, the Rotational-Enabled 7-DOF Seismometer will, by its nature, separate p- and s-waves into different data streams, simplifying signal processing and facilitating analysis of seismic source signatures and geological characterization. By adding measurements of three additional degrees-of-freedom at each level and leveraging the information from this new seismic observable, it is likely that an equally accurate picture of subsurface seismic activity could be garnered with fewer levels per hole. The key cost savings would come from better siting of the well due to increased information content and a decrease in the number of confirmation wells drilled, also due to the increase in information per well. Improved seismic tools may also increase knowledge, understanding, and confidence, thus removing some current blocks to feasibility and significantly increasing access to potential geothermal sites. During the Phase 1 effort summarized in this final report, the ATA Team modeled and built two TRL 3 proof-of-concept test units for two competing rotational sensor technologies. The two competing technologies were based on ATA's angular rate and angular displacement measurement technologies; Angular rate: ATA's Magnetohydrodynamic Angular Rate Sensor (Seismic MHD); and Angular displacement: ATA's Low Frequency Improved Torsional Seismometer (LFITS). In order to down-select between these two technologies and formulate a go / no go decision, the ATA Team analyzed and traded scientific performance requirements and market constraints against sensor characteristics and components, acquiring field data where possible to validate the approach and publishing results from these studies of rotational technology capability. Based on the results of Phase 1, the ATA Team finds that the Seismic MHD (SMHD) technology is the best choice for enabling rotational seismometry and significant technical potential exists for micro-seismic monitoring using a downhole 7-DOF device based on the SMHD. Recent technical papers and field data confirm the potential of rotational sensing for seismic mapping, increasing confidence that cost-reduction benefits are achievable for EGS. However, the market for geothermal rotational sensing is small and undeveloped. As a result, this report recommends modifying the Phase 2 plan to focus on prototype development aimed at partnering with early adopters within the geothermal industry and the scientific research community. The highest public benefit will come from development and deployment of a science-grade SMHD rotational seismometer engineered for geothermal downhole conditions and an integrated test tool for downhole measurements at active geothermal test sites.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Utah FORGE - Development of a Reservoir Seismic Velocity Model and Seismic Resolution Study

This is data from and a final report on the development of a 3D velocity model for the larger FORGE area and on the seismic resolution in the stimulated fracture volume at the bottom of well 16A-32. The velocity model was developed using RMS velocities of the seismic reflection survey and seismic velocity logs from borehole measurements as an input model. To improve the accuracy of the model in the shallow subsurface, travel times phase arrivals of the direct propagating P-waves were determined from the seismic reflection data, using PhaseNet, a deep-neural-network-based seismic arrival time picking method. The travel times were subsequently inverted using the input velocity model. The seismic resolution study used borehole and surface seismic sensors as well as the seismicity observed during the April 2022 stimulation experiment to estimate the seismic resolution in the activated fracture reservoir. The data contain a 3D P- and S-wave velocity model for the larger FORGE area.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

DEEPEN Leapfrog Geodata Model Cleaned and Reformatted Exploration Datasets from Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. As part of the DEEPEN 3D play fairway analysis (PFA) conducted at Newberry Volcano for multiple play types (conventional hydrothermal, superhot EGS, and supercritical), existing geoscientific exploration datasets needed to be acquired, cleaned, reformatted, and assembled in Leapfrog Geothermal. This GDR submission includes all of the cleaned and reformatted (X (m), Y (m), elevation (m), processed data values) datasets used to build the Leapfrog Geodata model. Existing datasets were acquired from the GDR, from AltaRock, and from other sources. This yielded the following datasets: - Digital elevation model produced from LiDAR data by Ramsey and Bard, 2016 - MT surveys from 2006, 2011, 2014, and 2017 (including single inversions) - Gravity surveys from 2006, 2007, and 2011 (including single) - Earthquake catalogs from PNSN, LLNL, and the Newberry EGS Demonstration project - Seismic velocity model from Templeton et al., 2014 - The Frone, 2015 temperature model and a new one produced through extrapolating downhole temperature measurements and the SMU temperature at depth maps. Two versions of the new model are provided: 250 m spacing and 500 m spacing - EarthVision geologic model with alteration from Moser et al., 2016 - Well data from EGS well 55-29, deep geothermal wells, coreholes (GEO N-2 through 5) and several thermal gradient holes - "Newberry Well Data:" Location, simple lithology, directional survey data, and temperature data for the 34 wells and coreholes used in the Newberry PFA Although there are additional 2D datasets available in the area, such as aeromagnetic surveys, these were not included in the analysis. While it may be possible to project these datasets into three dimensions by assuming the surface measurements do not vary with depth, this method is associated with high uncertainty. Preexisting inversions of these data were unavailable, and inverting additional geophysical datasets is outside the scope of this project.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Ambient field seismology in critical zone hydrological sciences

Passive ambient noise monitoring is an emerging tool in environmental seismology, leveraging the ambient seismic field to assess temporal variations in shallow subsurface properties. This review focuses on the potential and challenges of using scattered coda waves from noise correlation functions to monitor critical zone dynamics. The sensitivity of seismic velocities to various environmental factors, including precipitation, snowmelt, atmospheric pressure, and groundwater fluctuations, underscores the method’s versatility. While coda waves excel in detecting subtle changes due to their scattered nature, ballistic waves provide higher spatial resolution, albeit with challenges in source stability. Advances in seismic sensing, including distributed acoustic sensing and low-cost geophone networks, have enabled high-resolution monitoring of hydrological processes, subsurface deformation, and seismic hazards. Integrating seismic data with hydrological models provides insights into water storage, pore pressure changes, and soil moisture dynamics. However, limitations in spatial resolution, calibration with ground truth data, and coupled effects between environmental factors remain key challenges. This review emphasizes the importance of interdisciplinary approaches in refining methodologies, enhancing sensor deployments, and addressing data gaps. Passive seismic monitoring offers opportunities to understand critical zone processes and their broader impacts on seismic hazards and environmental sustainability.

58 GEOSCIENCES↗
DOE OSTI2020

Seismology in the Cloud: A New Streaming Workflow

Data-intensive research in seismology is experiencing a recent boom, driven in part by large volumes of available data and advances in the growing field of data science. However, there are significant barriers to processing large data volumes, such as long retrieval times from data repositories, complex data management, and limited computational resources. New tools and platforms have reduced the barriers to entry for scientific cluster computing, including the maturation of the commercial cloud as an accessible instrument for research. Here, we build a customized research cluster in the cloud to test a new workflow for large-scale seismic analysis, in which data are processed as a stream (retrieved on-the-fly and acted upon without storing), with data from the Incorporated Research Institutions for Seismology Data Management Center. We use this workflow to deploy a spectral peak detection algorithm over 5.6 TB of compressed continuous seismic data from 2074 stations of the USArray Transportable Array EarthScope network. Using a 50-node cluster in the cloud, we completed the noise survey in 80 hr, with an average data throughput of 1.7 GB per minute. By varying cluster sizes, we find the scaling of our analysis to be sublinear, due to a combination of algorithmic limitations and data center response times. The cloud-based streaming workflow represents an order-of-magnitude increase in acquisition and processing speed compared to a traditional download-store-process workflow, and offers the additional benefits of employing a flexible, accessible, and widely used computing architecture. It is limited, however, due to its reliance on Internet transfer speeds and data center service capacity, and may not work well for repeated analyses or those for which even higher data throughputs are needed. These research applications will require a new class of cloud-native approaches in which both data and analysis are in the cloud.

58 GEOSCIENCES↗
DOE OSTI2022

FutureGen VSP

Processed SEGY files from the FutureGen2 project. Data includes 5 categories of OVSP data with 15 lines in each category (all in TWT). VSP data was processed from Sterling Seismic.

FutureGen2↗
DOE OSTI2019

Integrated model construction for CO2-EOR monitoring via charged-wellbore casing controlled-source electromagnetics

We present a method for the integrated development of 3D electrical conductivity models for CO2-EOR monitoring with charged wellbore casing controlled-source electromagnetics (CWC-CSEM). The model is constructed through the multiphysics integration of reservoir simulation, seismic, borehole resistivity, and transient EM (TEM) data. The process is performed in two stages. First, a large background conductivity model is constructed from the combination of seismic, borehole, and TEM data. Second, a more detailed and dynamic component of the model is created within the production interval by converting the reservoir simulation parameters to conductivity through Archie’s equation. In this presentation, we demonstrate the complete workflow for building these 3D conductivity site models that can be updated throughout production using the Bell Creek oil field as an example. We then show application of the developed site model to simulating the predicted CWC-CSEM data as one step in the larger problem of reservoir imaging and monitoring of injected CO2 during enhanced oil recovery.

CO2, Sequestration, Monitoring, time-lapse, reserv↗
DOE OSTI2015

Crustal and uppermost mantle structure in the Middle East: assessing constraints provided by jointly modelling Ps and Sp receiver functions and Rayleigh wave group velocity dispersion curves

Seismic velocity models are found, along with uncertainty estimates, for 11 sites in the Middle East by jointly modelling Ps and Sp receiver functions and surface (Rayleigh) wave group velocity dispersion. The approach performs a search for models that satisfy goodness-of-fit criteria guided by a variant of simulated annealing and uses statistical tools to assess these products of searches. These tools, a parameter correlation matrix and marginal posterior probability density (PPD) function, allow us to evaluate quantitatively the constraints that each data type imposes on model parameters and to identify portions of each model that are well-constrained relative to other portions. This joint modelling technique, which we call ‘multi-objective optimization for seismology’, does not require a good starting solution, although such a model can be incorporated easily, if available, and can reduce the computation time significantly. Applying the process described above to broadband seismic data reveals that crustal thickness varies from 15 km beneath Djibouti (station ATD) to 45 km beneath Saudi Arabia (station RAYN). A pronounced low velocity zone for both Vp and Vs is present at a depth of ~12 km beneath station KIV located in northern part of greater Caucasus, which may be due to the presence of a relatively young volcano. Similarly, we also noticed a 6-km-thick low velocity zone for Vp beginning at 20 km depth beneath seismic station AGIN, on the Anatolian plateau, while positive velocity gradients prevail elsewhere in eastern Turkey. Beneath station CSS, located in Cyprus, an anomalously slow layer is found in the uppermost mantle, which may indicate the presence of altered lithospheric material. Crustal P- and S-wave velocities beneath station D2, located in the northeastern portion of central Zagros, range between 5.2–6.2 and 3.2–3.8 km s -1 , respectively. Finally, in Oman, we find a Moho depth of 34.0 ± 1.0 km and 25.0 ± 1.0 to 30.0 ± 1.0 km beneath stations S02 and S04, respectively.

58 GEOSCIENCES↗
DOE OSTI2024

Utah FORGE: GES Well 16A(78)-32 and Well 16B(78)-32 Stimulation Seismic Event Catalogs

This dataset contains seismic event catalogs from the hydraulic stimulation of wells 16A(78)-32 and 16B(78)-32 at the Utah FORGE site in April 2024. The data was collected by Geo Energy Suisse (GES) using a variety of seismic monitoring technologies, including 3-component (3C) geophones and distributed acoustic sensing (DAS) systems. These technologies were deployed across several locations, including wells 16A, 16B, and Delano-1, with sensor arrays at multiple depths to capture microseismic activity during the stimulations. The catalogs provide both real-time and manually checked seismic event locations, with detailed parameters such as trigger conditions, velocity models, and data acquisition settings. The dataset includes information on the stimulation stages, event rates, and hydraulic injection conditions for each well, with a report detailing the data acquisition configuration and seismic event location methodologies. Users will need to reference the included report for a complete understanding of the sensor network, data processing techniques, and accuracy considerations.

15 GEOTHERMAL ENERGY↗
NASA NTRS2000

Advanced Shipboard Communications Demonstrations with ACTS

For ships at sea. satellites provide the only option for high data rate (HDR), long haul communications. Furthermore the demand for HDR satellite communications (SATCOM) for military and commercial ships. and other offshore platforms is increasing. Presently the bulk of this maritime HDR SATCOM connectivity is provided via C-band and X-band. However, the shipboard antenna sizes required to achieve a data rate of, say T 1 (1.544 Mbps) with present C-/X-band SATCOM systems range from seven to ten feet in diameter. This limits the classes of ships to which HDR services can be provided to those which are large enough to accommodate the massive antennas. With its high powered K/Ka-band spot beams, the National Aeronautics and Space Administration's (NASA) Advanced Communications Technology Satellite (ACTS) was able to provide T I and higher rate services to ships at sea using much smaller shipboard antennas. This paper discusses three shipboard HDR SATCOM demonstrations that were conducted with ACTS between 1996 and 1998. The first demonstration involved a 2 Mbps link provided to the seismic survey ship MN Geco Diamond equipped with a 16-inch wide, 4.5-inch tall, mechanically steered slotted waveguide array antenna developed by the Jet Propulsion Laboratory. In this February 1996 demonstration ACTS allowed supercomputers ashore to process Geco Diamond's voluminous oceanographic seismic data in near real time. This capability allowed the ship to adjust its search parameters on a daily basis based on feedback from the processed data, thereby greatly increasing survey efficiency. The second demonstration was conducted on the US Navy cruiser USS Princeton (CG 59) with the same antenna used on Geco Diamond. Princeton conducted a six-month (January-July 1997) Western Hemisphere solo deployment during which time T1 connectivity via ACTS provided the ship with a range of valuable tools for operational, administrative and quality-of-life tasks. In one instance, video teleconferencing (VTC) via ACTS allowed the ship to provide life-saving emergency medical aid, assisted by specialists ashore. to a fellow mariner - the Master of a Greek cargo ship. The third demonstration set what is believed to be the all-time SATCOM data rate record to a ship at sea, 45 Mbps in October 1998. This Lake Michigan (Chicago area) demonstration employed one of ACTS' fixed beams and involved the smallest of the three vessels, the 45-foot Bayliner M/V Entropy equipped with a modified commercial-off-the-shelf one-meter antenna. A variety of multi-media services were provided to Entropy through a stressing range of sea states. These three demonstrations provided a preview of the capabilities that could be provided to future mariners on a more routine basis when K/Ka-band SATCOM systems are widely deployed.

Axford, Roy A.↗
DOE OSTI2025

A Comparison of Machine Learning Methods of Association Tested on Dense Nodal Arrays

The association of phase picks to form events is one of the fundamental components of seismology. Large and dense sensor networks, such as >1000 geophone arrays (and distributed acoustic sensing), offer unique challenges in association due to the vast numbers of observations and high likelihood of errant picks. In addition, the large number of stations can greatly increase the time it takes to perform the association. For this reason, machine learning (ML) methods might provide a more optimal method of association for such networks. In this work, we examine how well ML methods (e.g., Gaussian mixture model association, PhaseLink, and Graph Earthquake Neural Interpretation Engine) can incorporate dense seismic arrays into regional networks and how well they handle the increasing numbers of stations. Here, we test their capabilities on two dense seismic deployments, one within Rock Valley Nevada (52 nodes and a 9-station sparse local network), and the LArge-n Seismic Survey in Oklahoma dense nodal array (>1800 vertical-component geophones). Processing data from these two different styles of dense seismic deployments allows testing of how the ML algorithms can merge array data with a broader regional network, how they deal with poorly picked phases, and how they handle anthropogenic noise. We compare the ML-associated bulletins to those obtained using the Rapid Earthquake Association and Location algorithm, a more traditional method of association. We find that there are very small differences in results between the methods for small networks (<100 stations) with low pick rates. For large networks (>1000), there are enough errant picks that some of the ML methods start to create false events out of noise. We also find that the ML methods vary in computation time significantly but are all faster than the traditional method tested here.

58 GEOSCIENCES↗
DOE OSTI2018

Fallon FORGE: Seismic Reflection Profiles

Newly reprocessed Naval Air Station Fallon (1994) seismic lines: pre-stack depth migrations, with interpretations to support the Fallon FORGE (Phase 2B) 3D Geologic model. Data along seven profiles (>100 km of total profile length) through and adjacent to the Fallon site were re-processed. The most up-to-date, industry-tested seismic processing techniques were utilized to improve the signal strength and coherency in the sedimentary, volcanic, and Mesozoic crystalline basement sections, in conjunction with fault diffractions in order to improve the identification and definition of faults within the study area.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Small Seismic Events in Oklahoma Detected and Located by Machine Learning–Based Models

A complete earthquake catalog is essential to understand earthquake nucleation and fault stress. Following the Gutenberg–Richter law, smaller, unseen seismic events dominate the earthquake catalog and are invaluable for revealing the fault state. The published earthquake catalogs, however, typically miss a significant number of small earthquakes. Part of the reason is due to a limitation of conventional algorithms, which can hardly extract small signals from background noise in a reliable and efficient way. To address this challenge, we utilized a machine learning method and developed new models to detect and locate seismic events. These models are efficient in processing a large amount of seismic data and extracting small seismic events. We applied our method to seismic data in Oklahoma, United States, and detected ~14 times more earthquakes compared with the standard Oklahoma Geological Survey catalog. The rich information contained in the new catalog helps better understand the induced earthquakes in Oklahoma.

58 GEOSCIENCES↗