Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismic monitoring”

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 163 records · Page 9

Geophysical monitoring using active seismic techniques at the Citronelle Alabama CO 2 storage demonstration site

Between August 2012 and September 2014, about 114,000 metric tonnes of CO 2 was captured from the coal-fired Plant Barry Power Station at Bucks Alabama and injected into the Paluxy Formation above the oil pool in the southeast unit of the Citronelle Oilfield. Various monitoring methods were deployed at land surface and in project wells to measure system performance, comply with permit requirements and test new and innovative monitoring tools. The monitoring program relied heavily on active seismic methods for subsurface imaging of geologic structure and time-lapse seismic techniques to track the CO 2 migration in the injection interval. Both conventional geophone/hydrophone and fiber-optic based Distributed Acoustic Sensing (DAS) arrays were deployed and tested, allowing a side by side comparison of the equipment and techniques. Geophysical imaging of the subsurface was successful using DAS in the offset vertical seismic profile (OVSP) survey configuration. A high resolution OVSP image of the subsurface was obtained in 2014 with DAS, which exceeded project expectations in comparison to a lower resolution image obtained in 2012 using a conventional 80-level geophone array. A time-lapse image of the redistribution of CO 2 after injection ended in September 2014 was obtained with two DAS OVSP surveys from June 2014 and December 2015, thus successfully demonstrating its proof-of-concept. Unfortunately, a pre-injection baseline survey with DAS, which was in its initial stage of technology development in 2012, did not have sufficient quality for use, making it difficult to interpret the acquired DAS time-lapse difference. Additional research in this area has since demonstrated the utility of time-lapse DAS OVSP. DAS data were also acquired during a cross-well seismic survey conducted in 2014. Unfortunately, the DAS technique was not success in the cross-well survey configuration because the system noise level was too high in the crosswell frequency output range (100–1200 Hz) of the piezoelectric source (increasing by a factor of ten compared to VSP frequency band). Additionally, the cross-well geometry causes sub-horizontal (broadside) incidence on the vertical DAS fiber cable, which is known to be problematic. Current research is focused on improving the DAS cable response to broadside acoustic energy. Time-lapse seismic surveys using commercially available conventional arrays were also acquired. In contrast to the DAS acquired data, the cross-well seismic results obtained with the conventional array was highly successful and clearly showed the CO 2 remained in zone at the end of injection. Time-lapse differencing of the OSVP surveys acquired with the conventional arrays proved to be inconclusive. Finally, changes in wellbore conditions between surveys and unavoidable changes in equipment (the array used for the baseline survey was retired) affected data quality, making it difficult to interpret the OVSP results.

58 GEOSCIENCES↗

Adjoint Waveform Tomography for Next Generation Seismic Analyses and Monitoring

The development of methods and capabilities to compute complete waveform simulations in three-dimensional (3D) Earth models along with adjoint methods for computing the fully 3D sensitivity kernels in the 2000's set the stage for new advances in seismic imaging. I believe that the full benefits of adjoint waveform tomography (AWT) are not yet fully realized and this will be an important direction for the future of seismic tomography.

58 GEOSCIENCES↗

Aquifer Monitoring Using Ambient Seismic Noise Recorded With Distributed Acoustic Sensing (DAS) Deployed on Dark Fiber

Abstract Groundwater is a critical resource for human activities worldwide, and a vital component of many natural ecosystems. However, the state and dynamics of water‐bearing aquifers remain uncertain, mostly due to the paucity of subsurface data at high spatial and temporal resolution. Here, we show that analysis of infrastructure‐generated ambient seismic noise acquired on distributed acoustic sensing (DAS) arrays has potential as a tool to track variations in seismic velocities (d v / v ) caused by groundwater level fluctuations. We analyze 5 months of ambient noise acquired along an unused, 23 km‐long telecommunication fiber‐optic cable in the Sacramento Valley, CA, a so‐called “dark fiber." Three array subsections, ∼6 km apart, are processed and the stretching technique is applied to retrieve daily d v / v beneath each location. Near the Sacramento river, d v / v variations in the order of 2%–3% correlate with precipitation events and fluctuations in river stage of ∼1.5 m. In contrast, regions away (2.5 km) from the river do not experience large d v / v variations. These observations reveal short‐scale spatial variability in aquifer dynamics captured by this approach. Dispersion analysis and surface wave inversion of noise gathers reveal that seismic velocity perturbations occur at depths of 10–30 m. Rock physics modeling confirms that observed d v / v are linked to pore pressure changes at these depths, caused by groundwater table fluctuations. Our results suggest that DAS combined with ambient noise interferometry provides a means of tracking aquifer dynamics at high spatial and temporal resolutions at local to regional scales, relevant for effective groundwater resource management.

58 GEOSCIENCES↗

Predicting Seafloor Seismic Noise from Sea Surface Wave Heights and Vice Versa

Seismic hazard monitoring and global seismic tomography efforts are improved by recording seismic signals close to sources and at a variety of azimuths surrounding the source. While seismic networks have been deployed at many locations on land, seafloor stations are still sparse. Inclusion of ocean-based data would greatly improve the global coverage of seismic networks. However, the use of seafloor seismic data to complement land-based detection and characterization of events is complicated by the fact that ambient noise on the seafloor is in general much higher than noise levels on land. This noise is driven primarily by sea surface waves and tides. In this project, I define an empirical relationship between wave height and seafloor seismic displacement, and tie this relationship to the microseism noise theory of Longuet-Higgins and Ursell (1948). I find that while there appears to be a linear relationship between wave height and seismic displacement, wave height values can only be used with low confidence to predict seismic displacement and vice versa, and furthermore, factors other than wave height are significant in determining the level of seismic noise at a particular site.

58 GEOSCIENCES↗

Illinois State Geological Survey (ISGS), Illinois Basin - Decatur Project (IBDP) Seismic Data, July 7, 2021. Midwest Geological Sequestration Consortium (MGSC) Phase III Data Sets. DOE Cooperative Agreement No. DE-FC26-05NT42588.

Seismic data from the IBDP, included primarily under folders: Active_Seismic_Data and Passive_Seismic_Events_Monitoring. The data included under IBDP_Located_Microseismic_Event_Data are a subset of microseismic (Passive Seismic Events Monitoring) data acquired throughout the project from pre-injection to shut-in. These data represent the “located” microseismic events and the two datasets, located in folders Downhole_Geophone_Data and Surface_Seismometer_Data, are correlated.

3D Seismic,Carbon Sequestration,Decatur,Illinois B↗

Assessing Dynamic Time Warping Techniques for Discriminating Seismic Sources at Local and Regional Distances

Effective monitoring of seismic explosions and hazard assessment relies heavily on the accurate discrimination of underground seismic sources. This study investigates the application of novel nonlinear alignment techniques, specifically Dynamic Time Warping (DTW), for event-type discrimination at regional and local distances. Building on prior research that used DTW and Elastic Shape Analysis (ESA) in discrimination at regional distances, we evaluate the performance of recently developed variants of DTW, including a method that employs Pearson cross-correlation as a measure of warping distance and a time distortion coefficient that quantifies the type and degree of time distortion between signals. By analyzing observational datasets that include different source types, we assess the performance of these approaches for realistic monitoring scenarios. Specifically, we consider a dataset recorded at regional distances in the Korean Peninsula and a local-distance subset from the Unconstrained Utah Event Bulletin catalog to evaluate DTW-based discrimination across multiple distance scales. Additionally, we introduce the maximum cross-correlations of warped waveforms as a similarity metric for event classification. Through hierarchical cluster analysis and dendrogram interpretation, we present our findings, highlighting the strengths and limitations of these techniques in seismic event classification.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Applying Compression-Based Metrics to Seismic Data in Support of Global Nuclear Explosion Monitoring

The analysis of seismic data for evidence of possible nuclear explosion testing is a critical global security mission that relies heavily on human expertise to identify and mark seismic signals embedded in background noise. To assist analysts in making these determinations, we adapted two compression distance metrics for use with seismic data. First, we demonstrated that the Normalized Compression Distance (NCD) metric can be adapted for use with waveform data and can identify the arrival times of seismic signals. Then we tested an approximation for the NCD called Sliding Information Distance (SLID), which can be computed much faster than NCD. We assessed the accuracy of the SLID output by comparing it to both the Akaike Information Criterion (AIC) and the judgments of expert seismic analysts. Our results indicate that SLID effectively identifies arrival times and provides analysts with useful information that can aid their analysis process.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

A Wrapper to Use a Machine-Learning-Based Algorithm for Earthquake Monitoring

Seismology is one of the main sciences used to monitor volcanic activity worldwide. Fast, efficient, and accurate seismicity detectors are crucial to assess the activity level of a volcano in near–real time and to issue timely warnings. Traditional real–time seismic processing software uses phase onset pickers followed by a phase association algorithm to declare an event and estimate its location. The pickers typically do not identify whether the detected phase is a P or S arrival, which can have a negative impact on hypocentral location quality and complicates phase association. We implemented the deep–neural–network–based method PhaseNet to identify in real time P and S seismic waves on data from one– and three–component seismometers. We tuned the Earthworm binder_ew associator module to use the phase identification from PhaseNet to detect and locate the events, which we archive in a SeisComP3 database. We assessed the performance of the algorithm by comparing the results with existing catalogs built to monitor seismic and volcanic activity in Mayotte and the Lesser Antilles region. Our algorithm, which we refer to as PhaseWorm, showed promising results in both contexts and clearly outperformed the previous automatic method implemented in Mayotte. As a result, this innovative real–time processing system is now operational for seismicity monitoring in Mayotte and Martinique.

58 GEOSCIENCES↗

Findings on subtask 1.6 – basin electric carbon storage research project: novel monitoring techniques

The Energy & Environmental Research Center (EERC) led a study to validate novel and emerging technologies as commercial monitoring techniques for application in carbon dioxide (CO 2 ) injection operations. This applied research was conducted at Basin Electric Power Cooperative’s (Basin Electric’s) active CO 2 -injection operations in Beulah, North Dakota, in two phases. The EERC previously completed a set of baseline (preinjection) activities in Phase 1, which included 1) design of an automated, integrated, modular (AIM) monitoring station; 2) time-lapse electromagnetic (EM) field surveys; 3) drone-based surveillance studies; 4) time-lapse monitoring with seismic methods; and 5) advanced wellbore-monitoring methods.

54 ENVIRONMENTAL SCIENCES↗

Evaluating the Aftershock Duration of Induced Earthquakes

As the carbon sequestration community prepares to scale up the number and size of commercial operations, the need for tools and methods to assess and mitigate risks associated with these operations becomes increasingly important. One outstanding question is whether aftershocks of induced events decay quickly after injection operations cease or if aftershock activity persists for hundreds of years before returning to background levels more akin to tectonic events (Stein and Liu, 2009). Appropriate estimates of the aftershock duration impact several operational management decisions including mitigation strategies and post-injection monitoring for seismic activity. It is hypothesized that induced earthquake rates may diminish more quickly after injection is stopped, owing to higher stressing rates from injected fluids. Alternatively, it is plausible that only the first event in the sequence is induced by increased fluid overpressures, whereas subsequent events (e.g., aftershocks) respond to the stored tectonic stresses and static and dynamic stress changes due to the mainshock (Keranen et al., 2013). Further, if the aftershock duration can be linked to stressing rates due to injection, then it follows that operational strategies to reduce seismic hazard by reducing injection rates or volumes may be successful. However, if aftershocks of induced events are relieving stored tectonic stresses, then altering injection volumes may not alleviate ongoing seismic activity. Furthermore, knowledge of an aftershock duration could aid in the determination of an appropriate post-injection monitoring period for ongoing seismicity, which is a factor in overall operational costs. In this study, we model induced seismicity sequences in Oklahoma with a coupled Coulomb rate–state earthquake rate model (Dieterich, 1994; Kroll et al., 2017) to estimate aftershocks durations. Results for the current study indicate that elevated rates of aftershock activity following induced mainshocks return to background seismicity rates in less than five years, contrary to the tens to hundreds of years observed for tectonic aftershocks.

58 GEOSCIENCES↗

Seismic Signal Detection on International Monitoring System 3-Component Stations using PhaseNet

In this report we discuss training a deep learning seismic signal detection model on 3-component stations from the International Monitoring System (IMS) using the PhaseNet architecture. Using 14 years of associated signals from the International Data Centre’s (IDC) Late Event Bulletin (LEB), we auto-curated training data consisting of signal windows containing associated arrivals, and noise windows that contain no LEB-associated signals. We trained several models using different waveform window durations (30 seconds and 100 seconds), with and without bandpass filtering. We evaluated the effectiveness of our models using associated signals from the Unconstrained Global Event Bulletin (UGEB) and found that several of our models outperformed the signal detections from the IDC’s Selected Event List 3 (SEL3) arrival table. The SEL3 bulletin evaluated on the UGEB dataset with 100-second waveform windows registered a precision and recall of .15 and .48, respectively, versus .19 and .59 for our filtered-data model. For the 30-second waveform window dataset, the SEL3 bulletin achieved a precision and recall of .31 and .47, respectively, versus .32 and .60 for our filtered-data model. Finally, our models detected signals from all source-to-receiver distances, suggesting it is feasible to use a single PhaseNet model for the IMS network.

58 GEOSCIENCES↗

Report on the LLNL Global Seismic Waveform Tomography Modeling Project

Earth models have important applications for seismic event monitoring including the location and characterization of potential underground nuclear tests. Seismic waves generated by earthquakes and man-made events are altered by variations in Earth’s subsurface properties; and these distortions effect the arrival times of energy packets (phases) and obscure the characteristics of the original energy source (event mechanism). Three-dimensional models of Earth’s seismic properties can help explain and predict the distortions to the seismic wavefield and reveal the properties of a source that generated the waves, including the event location using the model-predicted timing of the phases and source mechanism using the modelpredicted waveform characteristics.

58 GEOSCIENCES↗

Assessing the Effectiveness of Generalized Likelihood Ratio Test Detector Schemes in Seismic Event Detection and the Avoidance of Nontarget Signals

Cross-correlation techniques have played a long-standing and pivotal role in seismic event monitoring. However, the performance of correlation-based detectors is challenged by nuisance seismicity, or nontarget signals. Such detections are a problem when the mission is to automatically map events to the correct source region. Using aftershocks of the 2014 $M_w$ 6.0 South Napa, California, earthquake, we demonstrate the effectiveness of utilizing a dynamic correlation processor framework in a generalized likelihood ratio test (GLRT) detector configuration to minimize nontarget detections. A GLRT maximizes a detection statistic with respect to one or more unknown parameters. In this case, the detection statistic is a template signal match against the waveform in a window sliding over a data stream, and the unknown parameter is an index variable indicating group membership of the template event or events. Detected events are assigned to the event group that yields the largest detection statistic. In this work, our results show that a GLRT detector will outperform a suite of independently operating correlation and subspace detectors in terms of having a lower nontarget detection rate at a given missed detection rate. We also show that a GLRT detector composed of a few high-rank subspace detectors has a slightly higher nontarget detection rate, but a significantly lower missed detection rate, than a GLRT detector composed of many low-rank subspace detectors. The high-rank GLRT configuration produced impressive results even with marginal data (single channel, single station, and very low time bandwidth product), which bodes well for the utility of building efficient aftershock classification systems and global monitoring systems at larger scales. However, future work is required to assess performance at the regional scale and to assess the performance of the system at detecting target events not used in the detector template creation.

58 GEOSCIENCES↗

Exploring Uncertainty in Moment Estimation for Small Earthquakes in Southern Nevada Using the Coda Envelope Method

Compiling source parameter estimates for small earthquakes is important both for our understanding of earthquake physics and for accurately assessing earthquake hazard. Reliable source parameter estimates are difficult to achieve for small earthquakes, in part due to our inability to accurately model the relevant physical processes at high frequencies. The coda envelope methodology developed by Mayeda and Walter (1996) and Mayeda et al. (2003) can mitigate this concern and estimate the moment of small earthquakes by determining the parameters that control the shape of the S-wave coda envelope while eliminating path effects by minimizing the scatter between seismic stations. Here, we use an open-source implementation of this technique called the Coda Calibration Tool (CCT; Barno, 2017) to calculate CCT-based moment magnitude estimates of small earthquakes (M L 0–3) in the Rock Valley, Nevada, region within the Nevada National Security Site. The Rock Valley data set is of particular interest because it allows us to explore the changes in uncertainties of the coda calibration method with earthquake size and depth. We found that a consistent linear relationship exists between the local magnitude M L and our coda-derived M w estimates for earthquakes as small as M L 0–3, but that current CCT workflows do not accurately characterize very shallow events. We also demonstrate that the epistemic uncertainty in the apparent stress value assumed by the CCT algorithm can influence magnitude estimates of small earthquakes. In conclusion, these results provide valuable insight into the seismicity of this region, and inform future analysis and modeling efforts for nuclear monitoring and seismic hazard.

58 GEOSCIENCES↗