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 127 records · Page 7

EGS Collab Experiment 1: 3D Seismic Velocity Model and Updated Microseismic Catalog from Double-Difference Seismic Tomography

This package contains a 3D Seismic velocity model and an updated microseismic catalog obtained for a double-difference seismic tomography study. The 3D_seismic_velocity_model text file contains x (m), y(m), z(m), P-wave velocity (km/s), P-wave velocity quality indicator (1 for well-constrained; 0 for poorly constrained), S-wave velocity (km/s), and S-wave velocity quality indicator (1 for well-constrained; 0 for poorly constrained). The Updated_MEQ_catalog text file contains event origin time, x(m), y(m), z(m), error in x (m), error in y (m), error in z (m), and RMS misfit (millisecond). The 3D_seismic_P-wave_velocity_model animation file shows slices of the 3D P-wave velocity model. The 3D_seismic_S-wave_velocity_model animation file shows slices of the 3D S-wave velocity model. The Interactive_MEQ_location_comparison API file is an interactive visualization of the updated microseismic event locations and the original seismic catalog. The visualization allows users to view and compare the event locations by dragging, rotating, and zooming in. An updated version of 3D_seismic_velocity_model and associated animations were included, which were calculated with a more strict assumption for quality indicators.

15 GEOTHERMAL ENERGY↗

Large earthquake processes in the northern Vanuatu subduction zone

Abstract The northern Vanuatu (formerly New Hebrides) subduction zone (11°S to 14°S) has experienced large shallow thrust earthquakes with M w > 7 in 1966 ( M S 7.9, 7.3), 1980 ( M w 7.5, 7.7), 1997 ( M w 7.7), 2009 ( M w 7.7, 7.8, 7.4), and 2013 ( M w 8.0). We analyze seismic data from the latter four earthquake sequences to quantify the rupture processes of these large earthquakes. The 7 October 2009 earthquakes occurred in close spatial proximity over about 1 h in the same region as the July 1980 doublet. Both sequences activated widespread seismicity along the northern Vanuatu subduction zone. The focal mechanisms indicate interplate thrusting, but there are differences in waveforms that establish that the events are not exact repeats. With an epicenter near the 1980 and 2009 events, the 1997 earthquake appears to have been a shallow intraslab rupture below the megathrust, with strong southward directivity favoring a steeply dipping plane. Some triggered interplate thrusting events occurred as part of this sequence. The 1966 doublet ruptured north of the 1980 and 2009 events and also produced widespread aftershock activity. The 2013 earthquake rupture propagated southward from the northern corner of the trench with shallow slip that generated a substantial tsunami. The repeated occurrence of large earthquake doublets along the northern Vanuatu subduction zone is remarkable considering the doublets likely involved overlapping, yet different combinations of asperities. The frequent occurrence of large doublet events and rapid aftershock expansion in this region indicate the presence of small, irregularly spaced asperities along the plate interface.

Cleveland, K. Michael↗

Conditioned Simulation of Ground-Motion Time Series at Uninstrumented Sites Using Gaussian Process Regression

Ground-motion time series are essential input data in seismic analysis and performance assessment of the built environment. Because instruments to record free-field ground motions are generally sparse, methods are needed to estimate motions at locations with no available ground-motion recording instrumentation. In this study, given a set of observed motions, ground-motion time series at target sites are constructed using a Gaussian process regression (GPR) approach, which treats the real and imaginary parts of the Fourier spectrum as random Gaussian variables. Model training, verification, and applicability studies are carried out using the physics-based simulated ground motions of the 1906 Mw 7.9 San Francisco earthquake and Mw 7.0 Hayward fault scenario earthquake in northern California. Additionally, the method’s performance is further evaluated using the 2019 Mw 7.1 Ridgecrest earthquake ground motions recorded by the Community Seismic Network stations located in southern California. These evaluations indicate that the trained GPR model is able to adequately estimate the ground-motion time series for frequency ranges that are pertinent for most earthquake engineering applications. The trained GPR model exhibits proper performance in predicting the long-period content of the ground motions as well as directivity pulses.

58 GEOSCIENCES↗

Utah FORGE: Updated Seismic Event Catalogue from the April, 2022 Stimulation of Well 16A(78)-32

These are revised catalogs, related to the April, 2022 well 16A(78)-32 stimulation (phases 1,2, & 3), provided by Geo Energie Suisse (GES) that include additional events at the start of Stage 1 and some tidying up of some locations. These catalogs also include events for additional events that were auto-located to provide a larger dataset for statistical analyses, like b-value calculations. The actual auto-locations have been removed to prevent spurious location plots being created. Times are recorded in UTC (Coordinate Universal Time), and the coordinate reference system is UTM Zone 12N, NAD83.

15 GEOTHERMAL ENERGY↗

Field Validation of MVA Technology for Offshore CCS: Novel Ultra-High-Resolution 3D Marine Seismic Technology (P-Cable) (Final Report)

The objectives of the proposed study were to deploy and validate a specific monitoring technology, high-resolution 3D marine seismic (HR3D), appropriate for large-demonstration and commercial-scale offshore CCS sites. The project accomplished successful acquisition two HR3D seismic surveys. The first HR3D dataset was over the offshore injection site of the Tomakomai, Japan integrated pilot CCS project, which at the time of survey acquisition was actively injecting CO 2 . The first survey also represented a successful international collaboration between the DOE NETL program and Japan’s national CCS program and was the first successful acquisition and use of HR3D over an active CO 2 injection site (Meckel, Feng et al. 2019). The Tomakomai HR3D survey successfully tested a novel 4-streamer HR3D system array in which, for the first time, no cross-cable (aka “P-Cable”) was utilized and only four GeoEel streamers were used instead of the standard 12-streamer configuration. Consequently, this was not, strictly speaking, a deployment of the “P-Cable” system of (Planke and Berndt 2004) but rather a modified version, thereof, and it is the first known demonstration of the modified system configuration. One very positive outcome from the Japanese collaboration earlier in the project was the ability to learn from the Japanese how they used tail buoys with GPS to determine the position of the seismic source and receivers in time and space. Based on that experience, GCCC designed and built six GPS receivers that could be used to position the streamer receivers and the seismic source via tail buoys. A fundamental advance that was made on the original design, was the ability to directly power the tail buoy GPS units and transfer data through the streamers (i.e., vs. the batteries used at Tomakomai). The bulkiness of the GPS batteries caused drag and episodic surging of the buoys, which affected data quality by lifting up the tail end of the streamers so the receivers were not at the same depth. The units were tested onshore for accuracy and functionality, and the design was subsequently and successfully tested in marine acquisition mode during the SLP survey acquisition. The marine acquisition test and survey satisfied Subtasks 2.2.2, Novel Positioning Technology Selection and Subtask 2.2.3, Novel Positioning Technology Deployment. Results of the novel positioning technology selection (Subtask 2.2.2) were considered successful and will be incorporated in future HR3D seismic acquisition projects to reduce costs, improve deployment safety at sea, and integrate both seismic and data recording via a single data transfer through the streamers to the recording system. The project also established a permitting process through NETL NEPA compliance, which included an Environmental Assessment in a marine setting and is required for conducting these types of surveys using Federal funding. The permitting process charted a “boilerplate,” which can allow future surveys related to other funded projects to move forward more expeditiously. Future improvements that could be considered are more robust seals on the GPS module and stronger materials (especially joints) on tail buoy fabrication. These would increase fixed costs, but would be advisable and probably more economic long-term if multiple HR3D surveys are planned. Project Accomplishments include: • Pre-survey Sensitivity Study • Marine geochemistry methods and data analysis • Successful HR3D seismic dataset acquired @ Tomakomai active CO 2 injection marine site • Developed advanced seismic processing techniques • No NRMS anomalies detected in overburden; Demonstration of containment • Repeatability study • Second survey collected @ San Luis Pass, TX • 4D application using positioning techniques developed in the project for monitoring were successful

3D seismic GPS positioning↗

GraphAlign: Graph-Enabled Machine Learning for Seismic Event Filtering

This report summarizes results from a 2 year effort to improve the current automated seismic event processing system by leveraging machine learning models that can operated over the inherent graph data structure of a seismic sensor network. Specifically, the GraphAlign project seeks to utilize prior information on which stations are more likely to detect signals originating from particular geographic regions to inform event filtering. To date, the GraphAlign team has developed a Graphical Neural Network (GNN) model to filter out false events generated by the Global Associator (GA) algorithm. The algorithm operates directly on waveform data that has been associated to an event by building a variable sized graph of station waveforms nodes with edge relations to an event location node. This builds off of previous work where random forest models were used to do the same task using hand crafted features. The GNN model performance was analyzed using an 8 week IMS/IDC dataset, and it was demonstrated that the GNN outperforms the random forest baseline. We provide additional error analysis of which events the GNN model performs well and poorly against concluded by future directions for improvements.

58 GEOSCIENCES↗

Discovering new events beyond the catalogue—application of empirical matched field processing to Salton Sea geothermal field seismicity

Using empirical matched field processing (MFP), we compare 4 yr of continuous seismic data to a set of 195 master templates from within an active geothermal field and identify over 140 per cent more events than were identified using traditional detection and location techniques alone. In managed underground reservoirs, a substantial fraction of seismic events can be excluded from the official catalogue due to an inability to clearly identify seismic-phase onsets. Empirical MFP can improve the effectiveness of current seismic detection and location methodologies by using conventionally located events with higher signal-to-noise ratios as master events to define wavefield templates that could then be used to map normally discarded indistinct seismicity. Since MFP does not require picking, it can be carried out automatically and rapidly once suitable templates are defined. In this application, we extend MFP by constructing local-distance empirical master templates using Southern California Earthquake Data Center archived waveform data of events originating within the Salton Sea Geothermal Field. We compare the empirical templates to continuous seismic data collected between 1 January 2008 and 31 December 2011. The empirical MFP method successfully identifies 6249 additional events, while the original catalogue reported 4352 events. The majority of these new events are lower-magnitude events with magnitudes between M0.2–M0.8. Here, the increased spatial-temporal resolution of the microseismicity map within the geothermal field illustrates how empirical MFP, when combined with conventional methods, can significantly improve seismic network detection capabilities, which can aid in long-term sustainability and monitoring of managed underground reservoirs.

15 GEOTHERMAL ENERGY↗

Utah FORGE: Results of Direct Shear Tests on Saturated Joints in Sierra White Granite

This dataset includes an Excel file with the results of direct shear tests to investigate the mechanical and geophysical response of saturated joints. Tests were conducted on induced tension-fractured Sierra White granite joints in a custom water pressurized chamber. The preparation of granite joints is summarized in the PDF file. The PDF file also includes a description of the direct shear tests, experimental procedure, and results.

15 GEOTHERMAL ENERGY↗

Utah FORGE: Seismic Reflection Data

This is 2D and 3D seismic reflection data from Utah FORGE reprocessed during Phase 2c. The readme file containing an explanation of the data including data formats, software that can be used, processing, and projection and datum used. The Reprocessing document gives the rationale for reprocessing and shows examples of the improvements that were obtained. For all 3D and 2D data the following datasets were created and output in SEG-Y format: - Unmigrated Time - Prestack Time Migration (PSTM), Unenhanced (UnEnh) and Enhanced (Enh) - Prestack Depth Migration (PSDM), Unenhanced (UnEnh) and Enhanced (Enh) - Velocity Model used for Migration

15 GEOTHERMAL ENERGY↗

Utah FORGE: Phase 2C Seismic Observation and Ground Water Well Data and Well Locations

This submission contains the following data associated with Utah FORGE Phase 2C within the Roosevelt Hot Springs geothermal area: - An ArcGIS shapefile with well locations for wells 58-32, 78-32, 68-32, OH-3, OH-4, WOW-2, and WOW-3. Well 58-32 is the deep test well, wells 68-32 and 78-32 are seismic test wells, and the OH and WOW wells were used for ground water monitoring and temperature gradient logging. - Drilling reports, lithologic logs, daily drilling reports, cuttings logs, and drilling photos for seismic monitoring wells 78-32 and 68-32. - Raw and post-processing temperature logs for wells OH-3, OH-4, WOW-2, and WOW-4. The archive also includes a README.txt file, which contains an explanation of the data and how it was acquired, and a text file with coordinates for the logged wells.

15 GEOTHERMAL ENERGY↗

The Modeled Performance of the ORNL Matched Filter Detector for Seismic Signatures of Multi-Axle Tractor Trailers

The MINOS program seeks methods to passively monitor the transfer of special materials, which are reproducible and quantitatively defensible. To date, we have demonstrated that short-period seismometers deployed near vehicle tracks can detect and monitor signatures of multi-axle tractor trailers impacting small road obstacles as means to detect and identify some transfer activities. In this report, we model the seismic velocity signal output from a tractor-trailer impacting a small obstacle (a template waveform) to parameterize a matched filter that identifies signals with similar shapes buried in noise, and that indicate repeating events (target waveforms). We then use this model to perform thousands of synthetic experiments, in which we apply our detector to identify noisy waveforms output from a similarly modeled vehicle data set; we process all data in the same manner as we process real data collected from the Oak Ridge National Laboratory (ORNL) seismic network. Our experiments, while limited, reveal that the matched filter shows a high discrimination capability, and at very low SNR values (at visibility thresholds). This means that template waveforms can match target waveforms that are sourced by models that share the same vehicle speed (target vehicles), and have a low probability of mis-identification. In particular, a matched filter designed to detect a five-axle tractor-trailer will reliably identify that target vehicle down to an SNR of 11.25 dB, and within ±0.1 km/hr, while successfully screening target vehicles with three axles.

42 ENGINEERING↗

WHOLESCALE: Microseismic Event Catalog for San Emidio, Nevada 2022

This submission includes a high-precision seismic event catalog estimated from seismic data collected at San Emidio, Nevada from April to May 2022. The catalog lists the precise time, location, and magnitude of microseismic events recorded during this period. Both the seismic data and microseismic event catalog were produced as part of the Water & Hole Observations Leverage Effective Stress Calculations and Lessen Expenses (WHOLESCALE) project. Attached here are the microseismic event catalog, a link to the raw seismic data, and a detailed description of methods used to create the catalog.

15 GEOTHERMAL ENERGY↗

Lunar seismic profiling experiment natural activity study

The Lunar Seismic Experiment Natural Activity Study has provided a unique opportunity to study the high frequency (4-20 Hz) portion to the seismic spectrum on the moon. The data obtained from the LSPE was studied to evaluate the origin and importance of the process that generates thermal moonquakes and the characteristics of the seismic scattering zone at the lunar surface. The detection of thermal moonquakes by the LSPE array made it possible to locate the sources of many events and determine that they are definitely not generated by astronaut activities but are the result of a natural process on the moon. The propagation of seismic waves in the near-surface layers was studied in a qualitative manner. In the absence of an adequate theoretical model for the propagation of seismic waves in the moon, it is not possible to assign a depth for the scattering layer. The LSPE data does define several parameters which must be satisfied by any model developed in the future.

Duennebier, F. K.↗

SMART Deliverable 6.1.2a: Application of the ORION tool to the IBDP Carbon Storage Site

Forecasting and managing potential induced seismic activity is one of the challenges facing commercialscale geologic carbon sequestration (GCS), as well as other geologic energy extraction and byproduct disposal technologies. Historically, the process to develop robust, science-based forecasts of induced seismicity has required an integrated effort from experts in seismology, geomechanics, and reservoir engineering to manage data, develop and evaluate models of subsurface processes, and to calibrate and interpret the results from a range of models to understand site behavior relative to prescribed standards and in the context of uncertainty in geologic characterization data, forecasting models, and operational scenario uncertainty. The Operational Forecasting of Induced Seismicity (ORION) toolkit is an open-source, observation-based forecasting toolkit that is being co-developed by two U.S. DOE-funded initiatives: the National Risk Assessment Partnership (NRAP) and the Scienceinformed Machine Learning for Accelerating Real Time Decisions in Subsurface Applications (SMART) Initiative. ORION is designed to provide functionality to support decision making about seismic hazard analysis and risk management for GCS stakeholders ranging from the public to site operators to expert seismologists. The tool, which is written as open-source code in the Python programming language, is composed of a desktop graphical user interface (GUI) and an underlying forecasting engine. The forecasting engine uses available reservoir properties, well and fluid injection scenario details, and observed seismic catalog data as inputs to produce a set of temporal and spatio-temporal seismic forecasts.

58 GEOSCIENCES↗

EGS Collab Project: Accomplishments and Plan

To better understand stimulation and inter-well flow in crystalline rock with the goal of developing and improving subsurface heat exchangers for enhanced geothermal energy recovery, the EGS (Enhanced Geothermal Systems) Collab project is developing highly monitored, well- characterized intermediate scale (~10-20 m) field test beds, stimulating, and performing flow tests. These tests provide the ability to evaluate concepts, data to test and validate models, and tools to be used for EGS development, which may be evaluated at DOE’s Frontier Observatory for Research in Geothermal Energy (FORGE) project. Numerical simulation is key to the project, for experiment design and interpretation of results. Testbed characterization included optical and acoustic televiewers, full waveform seismic, electromagnetic, gamma, and temperature. The test block was further characterized using seismic tomography (p- and s-), electrical resistance tomography (ERT), and hydrologic characterization including tracer tests. Coupled process models are constrained by the detailed site characterization data. Installed monitoring systems include: 1) passive seismic/acoustic emissions; 2) continuous active source seismic; 3) ERT, with dynamic imaging using high contrast fluids; and 4) distributed acoustic, temperature, and strain fiber optic sensors. During flow and stimulation tests, fracture aperture strain monitoring is performed using the Step-rate Injection Method for Fracture In-situ Properties (SIMFIP) tool. Pressure and flow in the injection and production boreholes are monitored, and tracer tests and wavefield imaging and inversion are performed to provide additional constraints. Stimulations have been performed at several locations, connecting our injection and production wells. An overview of the EGS Collab project is presented to provide a framework for other EGS Collab presentations, along with results from numerical simulations and experiments performed at the Sanford Underground Research Facility in South Dakota.

Enhanced Geothermal Systems, EGS Collab, stimulati↗

Seismoacoustic Data Fusion for Ground Coupled Airwaves: a FACT Case Study for Calibration and Array Processing

This report covers an inquiry into seismoacoustic array processing using infrasound arrivals combined with resulting Ground Coupled Airwaves (GCA) that are present on collocated seismic sensors. In preparation, data calibration and denoising is completed for a seismoacoustic sensor array that was deployed at the Facility for Acceptance, Calibration, and Testing on Kirtland Airforce Base from August through September of 2021. The events of interest for this study are small, local explosive sources that lead to short duration, impulsive signals on the instruments. The goal is to determine if combining infrasound signals with the corresponding GCAs on collocated seismic sensors can be used to improve the results returned by automated signal detection and characterization (e.g., back azimuth estimates). Preparation for seismic and infrasound data involves removing the instrument response so that sensors have flat power spectra over the frequency range 0.1-10 Hz, where signal from events of interest may be detected. After instrument response removal, deployment conditions specific to this array require a retrospective noise analysis to determine station emplacement characteristics. Once all data is calibrated, a manual search is performed for possible GCA arrivals across the seismoacoustic network. These arrivals are then processed through beamforming and subsequent event identification, resulting in a catalogue of seismoacoustic GCA arrivals with corresponding back azimuth and trace velocity estimations.

47 OTHER INSTRUMENTATION↗

Evaluation of Seismic Artificial Intelligence with Uncertainty

Artificial intelligence has transformed the seismic community with deep learning models (DLMs) that are trained to complete specific tasks within workflows. However, there is still a lack of robust evaluation frameworks for evaluating and comparing DLMs. Here, we address this gap by designing an evaluation framework that jointly incorporates two crucial aspects: performance uncertainty and learning efficiency. To target these aspects, we meticulously construct the training, validation, and test splits using a clustering method tailored to seismic data and enact an expansive training design to segregate performance uncertainty arising from stochastic training processes and random data sampling. The framework’s ability to guard against misleading declarations of model superiority is demonstrated through the evaluation of PhaseNet (Zhu and Beroza, 2018), a popular seismic phase picking DLM, under three training approaches. Our framework helps practitioners choose the best model for their problem and set performance expectations by explicitly analyzing model performance with uncertainty at varying budgets of training data.

58 GEOSCIENCES↗

Deformation analysis of Aceh April 11{sup th} 2012 earthquake using GPS observation data

This research tries to estimate the co-seismic deformation of intraplate earthquake occurred off northern Sumatra coast which is about 100-200 km southwest of Sumatrasubduction zone. The earthquake mechanism was strike-slip with magnitude 8.6 and triggering aftershock with magnitude 8.2 two hours later. We estimated the co-seismic deformation by using the GPS (Global Positioning System) continuous data along western Sumatra coast. The GPS observation derived from Sumatran GPS Array (SuGAr) and Geospatial Information Agency (BIG). For data processing we used GPS Analyze at Massachusetts Institute of Technology (GAMIT) software and Global Kalman Filter (GLOBK) to estimate the co-seismic deformation. From the GPS daily solution, the result shows that the earthquake caused displacement for the GPS stations in Sumatra. GPS stations in northern Sumatra showed the displacement to the northeast with the average displacement was 15 cm. The biggest displacement was found at station BSIM which is located at Simeuleu Island off north west Sumatra coast. GPS station in middle part of Sumatra, the displacement was northwest. The earthquake also caused subsidence for stations in northern Sumatra, but from the time series there was not sign of subsidence was found at middle part of Sumatra. In addition, the effect of the earthquake was worldwide and affected the other GPS Stations around Hindia oceanic.

58 GEOSCIENCES↗