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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Locating underground features with seismic data processing

Methods are presented for determining the location of underground features (e.g., CO 2 ). One method includes capturing, by sensors distributed throughout a region, seismic traces associated with seismic signals generated by a seismic source. For multiple sensors, active noise is identified or passive noise is measured within each seismic trace and values for attributes associated with the active or passive noise are determined. Further, an unsupervised machine-learning model, based on the values of the attributes, is utilized to determine noise characteristics for multiple sensors. The sensors are grouped in clusters based on the noise characteristics for each sensor. For multiple clusters, a noise filter is created based on the noise characteristics of the sensors in the cluster, and the noise filter of the cluster is applied, for multiple sensors, to the seismic traces of the sensor. Additionally, the filtered seismic traces are analyzed to determine a location of CO 2 underground.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Rock Valley Accelerated Weight Drop Seismic Data Processing and Picking of P-wave and S-wave Arrival Times

Rock Valley, in the southern end of the Nevada National Security Site, hosts a fault system that was responsible for a shallow (< 3 km below surface ) magnitude 3.7 earthquake in May 1993. In order to better understand this system, seismic properties of the shallow subsurface need to be better constrained. In April and May of 2021, accelerated weight drop (AWD) active-source seismic data were recorded in order to measure P- and S-wave travel-times for the area. This report describes the processing and phase picking of the recorded seismic waveforms. In total, we picked 7,982 P-wave arrivals at offsets up to ~2500 m, and 4,369 S-wave arrivals at offsets up to ~2200 m. These travel-time picks can be inverted for shallow P-wave and S-wave velocity structure in future studies.

58 GEOSCIENCES↗
DOE OSTI2015

Multiple attenuation to reflection seismic data using Radon filter and Wave Equation Multiple Rejection (WEMR) method

Separation between signal and noise, incoherent or coherent, is important in seismic data processing. Although we have processed the seismic data, the coherent noise is still mixing with the primary signal. Multiple reflections are a kind of coherent noise. In this research, we processed seismic data to attenuate multiple reflections in the both synthetic and real seismic data of Mentawai. There are several methods to attenuate multiple reflection, one of them is Radon filter method that discriminates between primary reflection and multiple reflection in the τ-p domain based on move out difference between primary reflection and multiple reflection. However, in case where the move out difference is too small, the Radon filter method is not enough to attenuate the multiple reflections. The Radon filter also produces the artifacts on the gathers data. Except the Radon filter method, we also use the Wave Equation Multiple Elimination (WEMR) method to attenuate the long period multiple reflection. The WEMR method can attenuate the long period multiple reflection based on wave equation inversion. Refer to the inversion of wave equation and the magnitude of the seismic wave amplitude that observed on the free surface, we get the water bottom reflectivity which is used to eliminate the multiple reflections. The WEMR method does not depend on the move out difference to attenuate the long period multiple reflection. Therefore, the WEMR method can be applied to the seismic data which has small move out difference as the Mentawai seismic data. The small move out difference on the Mentawai seismic data is caused by the restrictiveness of far offset, which is only 705 meter. We compared the real free multiple stacking data after processing with Radon filter and WEMR process. The conclusion is the WEMR method can more attenuate the long period multiple reflection than the Radon filter method on the real (Mentawai) seismic data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Modeling of time-lapse seismic monitoring using CO 2 leakage simulations for a model CO 2 storage site with realistic geology: Application in assessment of early leak-detection capabilities

Time-lapse surface seismic surveys have been widely used at carbon sequestration sites for site characterization, monitoring subsurface CO 2 plume migration, and detecting potential CO 2 leakage from a storage reservoir. Monitoring in the first permeable unit directly above the primary seal is important for early detection of CO 2 leakage. Forward modeling of time-lapse seismic data can be used to assess the utility of the seismic method for CO 2 leakage detection. We develop a workflow for forward modeling of time-lapse seismic data, including constructing seismic velocity models using flow simulation outputs, modeling of pre-stack and post-stack synthetic seismic data following seismic data processing sequence and analysis of processed synthetic time-lapse seismic data. We apply the forward modeling and analysis workflow to assessing the detectability and the earliest detection time of seismic monitoring using the hypothetical CO 2 leakage scenarios for a model geologic storage site with realistic geology. We then derive the detection thresholds using the simulated normalized root-mean-square (NRMS) differences for the no-leakage case at a range of signal-to-noise ratios, representing the background noise levels in seismic data. We then compare NRMS differences triggered by the CO 2 leakage to the detection thresholds at each time step to quantify the detectability and the earliest detection time of seismic monitoring. We analyze the effects of the acquisition parameters and elastic parameters on the produced synthetic seismic data and earliest detection time. Our modeling results indicate that high signal-to-noise ratio is needed to detect the CO 2 leakage at the model site. Minimizing the background noise in seismic data is crucial for improving the detectability of the seismic method. Increasing the shot density or increasing the dominant frequency of the source wavelet is likely to increase the possibility of the leakage detection and reduce the time needed for the detection. The elastic parameters used in the rock physics modeling have significant effects on the resultant seismic velocity models and synthetic seismic data, highlighting their critical roles in understanding and interpreting time-lapse seismic reflection data associated with CO 2 monitoring, verification and accounting activities.

58 GEOSCIENCES↗
DOE OSTI2021

WHOLESCALE: Seismic Survey Data from San Emidio Nevada 2021

This dataset includes raw and processed seismic data from the 2021 seismic survey at the San Emidio geothermal field in Nevada. In April and May 2021, 37 tri-axial short period seismographs were deployed in a 1.8km diameter cluster centered on 40.367278 N, 119.409019 W. The first data record started at 2021-04-06T07:09:10Z UTC and the last record ended at 2021-05-11T02:58:52Z UTC. The 37 stations included 29 SmartSolo IGU-16HR 3C all-in-one 5 Hz seismographs and 8 DataCube seismographs with 4.5 Hz HGS HG-6(B coil) tri-axial geophones. The raw format (level 0) data includes 353 GB of 500 sps SmartSolo data in native DLD format, 113 GB of 400 sps DataCube data in native DataCube format, and 3.4 GB of GPS data collected during the RTK GPS survey. The SAC data (level 1) files include 564 GB of hourly SAC files. The experiment was designed to monitor seismic activity before, during, and after the planned three-day plant maintenance shutdown April 19-21, 2021. The pumping stopped at 2021/04/19 12:51:45 UTC and resumed about 2021/04/21 21:00:00 UTC. The dataset is hosted in an AWS data lake, with an associated GDR data set providing the metadata. Links to both of these resources are included below. Additionally, this collection features data and metadata from a 2016 seismic survey at the same site for comparative analysis.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

New Architecture to Support Integration and Processing of Seismic Data from Heterogeneous Sources

The Geophysical Monitoring Program (GMP) at Lawrence Livermore National Lab (LLNL) maintains a database and supporting infrastructure for geophysical data used in support of the Nuclear Detonation Detection mission. This database includes data from multiple sources, many of which do not distribute data to the public or for which there is no automated means of access. For example, Figure 1 shows (left) the distribution of waveform data in our database by source. The Incorporated Research Institutions for Seismology Data Management Center (IRISDMC) is our major source of waveform data and those data may be retrieved at will using the Federated Digital Seismograph Networks FDSN web Application Programming Interface (API). However, the next 6 most important sources of waveform data have no or only limited automated access to waveforms. As Figure 1 (right) shows, it is very common for waveform records associate with an event in our database to come from two or more sources, and in some cases data come from 10 sources. This diversity of data sources drives our need for efficient and correct integration of metadata, parametric data, and waveform data.

58 GEOSCIENCES↗
DOE OSTI2023

A Novel Approach to Map Permeability Using Passive Seismic Emission Tomography

Newly acquired magnetotelluric data and passive seismic data collected with tightly spaced geophone arrays are combined with historic drilling, active seismic, and potential fields data to generate 3-D permeability maps. A cooperative inversion methodology has been developed using active seismic, magnetotelluric, and gravity data in order to produce more robust velocity models for passive seismic data processing without requiring expensive 3-D active seismic surveys. The cooperative inversion estimates velocities from other geophysical data where no prior seismic velocity information is available at two geothermal sites in Nevada: San Emidio and Crescent Valley.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Cardinal: Seismic and Geoacoustic Array Processing

Data collected via seismic and infrasound array deployments are leveraged in the geosciences to detect and characterize a myriad of natural and anthropogenic sources. These deployments consist of numerous sensors placed in a predetermined configuration to amplify signal strength and improve the efficacy of array processing techniques used to measure signal directionality and waveform coherence. High‐fidelity feature extraction is often predicated on interstation distance as well as the frequency content and wavelength of an incident signal. Numerous array processing softwares analyze data in sequential frequency bands to obtain a more detailed characterization of a signal. However, current algorithms are limited in their ability to determine optimal array configuration for each band. We introduce an open‐source Python code, called Cardinal, to process seismic and infrasound array data in discretized time–frequency space with the option of applying an adaptive array design to determine optimal subarray configuration for each frequency band. To reduce computational time, the array processing step can be run in parallel using multithreading. Furthermore, the software has the capability to aggregate array processing results from different time–frequency pixels to produce separate sets of detections, or families, with added utility via the application of an adaptive semblance threshold, which aids in isolating signals‐of‐interest from coherent background noise. Upon appropriate configuration, Cardinal exhibits the potential to combine distinct seismic and infrasound phases into separate families.

Adaptive Array↗
DOE OSTI2016

Large-scale seismic waveform quality metric calculation using Hadoop

Here in this work we investigated the suitability of Hadoop MapReduce and Apache Spark for large-scale computation of seismic waveform quality metrics by comparing their performance with that of a traditional distributed implementation. The Incorporated Research Institutions for Seismology (IRIS) Data Management Center (DMC) provided 43 terabytes of broadband waveform data of which 5.1 TB of data were processed with the traditional architecture, and the full 43 TB were processed using MapReduce and Spark. Maximum performance of ~0.56 terabytes per hour was achieved using all 5 nodes of the traditional implementation. We noted that I/O dominated processing, and that I/O performance was deteriorating with the addition of the 5th node. Data collected from this experiment provided the baseline against which the Hadoop results were compared. Next, we processed the full 43 TB dataset using both MapReduce and Apache Spark on our 18-node Hadoop cluster. We conducted these experiments multiple times with various subsets of the data so that we could build models to predict performance as a function of dataset size. We found that both MapReduce and Spark significantly outperformed the traditional reference implementation. At a dataset size of 5.1 terabytes, both Spark and MapReduce were about 15 times faster than the reference implementation. Furthermore, our performance models predict that for a dataset of 350 terabytes, Spark running on a 100-node cluster would be about 265 times faster than the reference implementation. We do not expect that the reference implementation deployed on a 100-node cluster would perform significantly better than on the 5-node cluster because the I/O performance cannot be made to scale. Finally, we note that although Big Data technologies clearly provide a way to process seismic waveform datasets in a high-performance and scalable manner, the technology is still rapidly changing, requires a high degree of investment in personnel, and will likely require significant changes in other parts of our infrastructure. Nevertheless, we anticipate that as the technology matures and third-party tool vendors make it easier to manage and operate clusters, Hadoop (or a successor) will play a large role in our seismic data processing.

58 GEOSCIENCES↗
DOE OSTI2013

Advanced Seismic Data Analysis Program (The Hot Pot Project), DOE Award: DE-EE0002839 (Phase 1 Report)

A five-line (23 mile) reflection- seismic survey was conducted at the Hot Pot geothermal prospect area in north-central Nevada under the USDOE (United States Department of Energy) Geothermal Technologies Program. The project objective was to utilize innovative seismic data processing, integrated with existing geological, geophysical and geochemical information, to identify high-potential drilling targets and to reduce drilling risk. Data acquisition and interpretation took place between October 2010 and April 2011. The first round of data processing resulted in large areas of relatively poor data, and obvious reflectors known from existing subsurface information either did not appear on the seismic profiles or appeared at the wrong depth. To resolve these issues, the velocity model was adjusted to include geologic input, and the lines were reprocessed. The resulting products were significantly improved, and additional detail was recovered within the high-velocity and in part acoustically isotropic basement. Features visible on the improved seismic images include interpreted low angle thrust faults within the Paleozoic Valmy Formation, which potentially are reactivated in the current stress field. Intermediate-depth wells are currently targeted to test these features. The seismic images also suggest the existence of Paleogene sedimentary and volcanic rocks which potentially may function as a near-surface reservoir, charged by deeper structures in Paleozoic rocks.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Imagery Applications for Advanced Event Analytics (WBS 24.3.1.3.3-IDC FY202 Final Project Report)

Accurate event locations and replicability of location analyses are essential for assessing the nature of an event, its context, ambient site conditions, and proximity to relevant facilities and infrastructure. Additionally, accurate event locations provide valuable information that reduce uncertainties, improve confidence in event analyses, and inform in-field verification activities. However, event location/relocation and replicability are difficult due to a number of factors, including spatially-sparse network coverage in some areas of the globe and variability in seismic data processing. This team proposed that the incorporation of high-fidelity imagery as a data backbone to the analytical assessment of a suspected underground explosion and/or an advanced seismic event bulletin produced by the International Data Centre (IDC) of the Preparatory Commission for the Comprehensive Nuclear-Test-Ban Treaty Organization (CTBTO PrepCom) could reduce uncertainties and improve confidence in analyses. Specifically, temporally-separated images can reduce uncertainty by identifying areas where change has occurred (e.g., building construction or demolition, road or facilities improvements). The primary goal of this project was to develop an automated geospatial processing script for imagery change detection to better reflect needs of the technical community (including the IDC) and to make the use of such a tool accessible in a variety of settings across platforms. Technical experts at Los Alamos National Laboratory successfully built GAIA: the Geospatial Automated Imagery Analysis tool, to fill this need. GAIA combines five tool components to produce orthorectified time separated imagery and imagery change detection maps. Our toolkit (1) reduces error by providing a standardized workflow for image analyses and (2) significantly reduces processing time from between 7 and 24+ hours to approximately 5 minutes. Technical experts at Sandia National Laboratories supported GAIA via beta-testing and by introducing a web-based system approach for increased applicability. To test the function, performance, broad application, and ease-of-use of GAIA, we applied it to four separate test cases. The results of this preliminary investigation show promise in reducing uncertainty in seismic event locations: if satellite imagery can show regions where operations that produce seismic activity likely occurred, then pursuing imagery to locate epicenters of seismic nuclear events could reduce the time needed to find the true epicenter location.

58 GEOSCIENCES↗
DOE OSTI2020

Imagery Applications for Advanced Event Analytics (WBS 24.3.1.3.3-IDC FY2020 Final Project Report)

Accurate event locations and replicability of location analyses are essential for assessing the nature of an event, its context, ambient site conditions, and proximity to relevant facilities and infrastructure. Additionally, accurate event locations provide valuable information that reduce uncertainties, improve confidence in event analyses, and inform in-field verification activities. However, event location/relocation and replicability are difficult due to a number of factors, including spatially-sparse network coverage in some areas of the globe and variability in seismic data processing. This team proposed that the incorporation of high-fidelity imagery as a data backbone to the analytical assessment of a suspected underground explosion and/or an advanced seismic event bulletin produced by the International Data Centre (IDC) of the Preparatory Commission for the Comprehensive Nuclear-Test-Ban Treaty Organization (CTBTO PrepCom) could reduce uncertainties and improve confidence in analyses. Specifically, temporally-separated images can reduce uncertainty by identifying areas where change has occurred (e.g., building construction or demolition, road or facilities improvements). The primary goal of this project was to develop an automated geospatial processing script for imagery change detection to better reflect needs of the technical community (including the IDC) and to make the use of such a tool accessible in a variety of settings across platforms. Technical experts at Los Alamos National Laboratory successfully built GAIA: the Geospatial Automated Imagery Analysis tool, to fill this need. GAIA combines five tool components to produce orthorectified time-separated imagery and imagery change detection maps. Our toolkit (1) reduces error by providing a standardized workflow for image analyses and (2) significantly reduces processing time from between 7 and 24+ hours to approximately 5 minutes. Technical experts at Sandia National Laboratories supported GAIA via beta-testing and by introducing a web-based system approach for increased applicability. To test the function, performance, broad application, and ease-of-use of GAIA, we applied it to four separate test cases. The results of this preliminary investigation show promise in reducing uncertainty in seismic event locations: if satellite imagery can show regions where operations that produce seismic activity likely occurred, then pursuing imagery to locate epicenters of seismic nuclear events could reduce the time needed to find the true epicenter location.

58 GEOSCIENCES↗
DOE OSTI2019

Global- and Local-Scale High-Resolution Event Catalogs for Algorithm Testing

During the development of new seismic data processing methods, the verification of potential events and associated signals can present a nontrivial obstacle to the assessment of algorithm performance, especially as detection thresholds are lowered, resulting in the inclusion of significantly more anthropogenic signals. Here, we present two 14 day seismic event catalogs, a local–scale catalog developed using data from the University of Utah Seismograph Stations network, and a global–scale catalog developed using data from the International Monitoring System. Each catalog was built manually to comprehensively identify events from all sources that were locatable using phase arrival timing and directional information from seismic network stations, resulting in significant increases compared to existing catalogs. Here, the new catalogs additionally contain challenging event sequences (prolific aftershocks and small events at the detection and location threshold) and novel event types and sources (e.g., infrasound only events and long–wall mining events) that make them useful for algorithm testing and development, as well as valuable for the unique tectonic and anthropogenic event sequences they contain.

58 GEOSCIENCES↗
DOE OSTI2021

Seismically Detecting Nuclear Reactor Operations Using a Power Spectral Density (PSD) Misfit Detector

To explore the ability to indirectly detect and attribute various operations conducted at a nuclear reactor using waveform data, we investigated the seismic signals recorded near the High Flux Isotope Reactor (HFIR) located at Oak Ridge National Laboratory in Oak Ridge, Tennessee. Specifically, we processed seismic data collected from a single seismoacoustic station, WACO, near the HFIR facility, and employed a power spectral density misfit detector to identify signals of interest and associate the detections with operational events. Initial results suggest that this method provides a promising means of regularly detecting at least 19 unique operations. Furthermore, with additional station deployment and more comprehensive data logs, we anticipate that future analysis will offer an additional means to seismically monitor nuclear reactors (such as HFIR) health and performance more accurately.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Shallow Marine High-Resolution 3D Seismics for Above-Zone CO2 Monitoring: A Case Study in Offshore Japan

Shallow marine high-resolution 3D (HR3D) seismic utilizes a novel seismic acquisition technology for imaging fine geological features not revealed by conventional 3D surveys. HR3D technology has been previously used for shallow subsurface studies, however is rarely implemented for monitoring injected fluids, and has not previously been used to monitor an active CO2 storage project. In this study, we for the first time demonstrate the acquisition, processing and initial interpretation of a HR3D seismic survey in Tomakomai, offshore Japan, where a large-scale CO2 storage project is being undertaken. As a seismic monitoring tool, the HR3D survey images the geologic overburden (< 1000 ms), aiming to detect potential CO2 leakage pathways above the injection reservoir. The HR3D acquisition system comprises a 210 in3 GI air gun and 4 streamers that are 25 m in length with 10 m inline separation. Each streamer has 8 channels with a 3.125 m group interval, yielding a very small final bin size of 3 x 3 m. The short offset and low-fold coverage of the HR3D data creates new challenges in seismic data processing, which demands a different data processing workflow tailored to the HR3D acquisition. We present a data processing workflow involving multiple software processing packages, which tackled high-resolution related static issues and significantly reduced various noise elements. Dip-steering and structural filtering were employed on the stacked volume to further enhance the lateral coherence and to remove random noise and suppress acquisition footprints. The final volume shows flat and coherent above-zone stratigraphy, suggesting a good seal quality above the CO2 reservoir. The successful implementation of HR3D acquisition and the customized processing workflow verify the feasibility of HR3D seismics for the above- zone monitoring of offshore CO2 injection and storage projects.

Feng, Ye Erickson↗