Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SEISMIC”

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 253 records · Page 14

Using Ambient Noise Seismic Interferometry and Local and Teleseismic Earthquakes to Determine Crustal Thickness and Moho Structure of the Northwestern Gulf of Mexico Margin

Abstract The northwestern part of the Gulf of Mexico has undergone two episodes of continental rifting and collision and produced structural artifacts that are now buried under many kilometers of sediments, complicating investigations of the region. The deep sedimentary package precludes outcrops and points to a need for the application of seismic techniques, but low rates of seismicity in the region and sparse seismic monitoring limit the utility of traditional seismic methods. We therefore use diverse data to perform two‐dimensional seismic tomography across the dry land portion of the margin. Data are gleaned from teleseismic and regional earthquakes, postcritical SsPmp arrivals, and direct P wave energy identified with seismic interferometry that were recorded by a broadband, three‐component array and partially overlapping short‐period, vertical‐component array. The Pn and postcritical SsPmp phase help constrain the Moho discontinuity, which a previous receiver function study suggested was absent beneath the seaward portion of this transect. A high‐velocity body is observed in the crust at the same location as the Houston Magnetic Anomaly, possibly marking rocks from the Alleghenian continental assembly. The crust thins from NW to SE, indicating that extension occurred mostly to the south of the Ouachita orogeny. Our model indicates that the margin's sediment package reaches a maximum thickness of ~15 km at the coast and becomes unresolvably thin near the Llano Uplift.

58 GEOSCIENCES↗

Identifying Different Classes of Seismic Noise Signals Using Unsupervised Learning

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

Johnson, Christopher W.↗

Modeling Seismic Recordings of High-Frequency Guided Infrasound on Mars

NASA's Interior Exploration using Seismic Investigations, Geodesy and Heat Transport (InSight) mission records several high-frequency (>0.5 Hz) dispersive seismic signals on Mars. These signals are due to the acoustic-to-seismic coupling of infrasound generated by the entry and impact of meteorites. This dispersion property is due to infrasound propagating in a structured atmosphere, and we refer to this dispersive infrasound as guided infrasound. We propose to model the propagation of guided infrasound and the seismic coupling to the ground analytically; we use a 1D layered atmosphere on a three-layer solid subsurface medium. The synthetic ground movements fit the observed dispersive seismic signals well and the fitting indicates that the regolith beneath InSight is about 40-m in thickness. We also examine and validate the previously-published subsurface models derived from InSight ambient seismic vibration data.

58 GEOSCIENCES↗

Probing the evolution of fault properties during the seismic cycle with deep learning

We use seismic waves that pass through the hypocentral region of the 2016 M6.5 Norcia earthquake together with Deep Learning (DL) to distinguish between foreshocks, aftershocks and time-to-failure (TTF). Binary and N-class models defined by TTF correctly identify seismograms in test with > 90% accuracy. We use raw seismic records as input to a 7 layer CNN model to perform the classification. Here we show that DL models successfully distinguish seismic waves pre/post mainshock in accord with lab and theoretical expectations of progressive changes in crack density prior to abrupt change at failure and gradual postseismic recovery. Performance is lower for band-pass filtered seismograms (below 10 Hz) suggesting that DL models learn from the evolution of subtle changes in elastic wave attenuation. Tests to verify that our results indeed provide a proxy for fault properties included DL models trained with the wrong mainshock time and those using seismic waves far from the Norcia mainshock; both show degraded performance. Our results demonstrate that DL models have the potential to track the evolution of fault zone properties during the seismic cycle. If this result is generalizable it could improve earthquake early warning and seismic hazard analysis.

58 GEOSCIENCES↗

Alternate Null Hypothesis Correlation: A New Approach to Automatic Seismic Event Detection

This paper describes a new method of seismic signal detection that improves upon the conventional waveform correlation method. Recent studies suggested that a significant limiting factor in the application of waveform correlation to regional and global scale monitoring is the false alarm rate. The false alarms do not originate from detections on noise but rather from seismic arrivals with unrelated source locations. This article presents results from an approach to waveform correlation that exploits techniques from signal processing and machine learning to improve the accuracy of detecting seismic arrivals. We modify the detection model for waveform correlation such that transient signals from noncollocated seismicity are considered when designing the detectors. The new approach uses waveform templates from known catalog events to train a supervised machine learning algorithm that derives a new set of detectors to represent the unique characteristics of the template waveforms; these new detectors maximize the likelihood of detecting only the desired events, thereby minimizing false alarms. We train a waveform correlation template library for a single three-component seismic monitoring station. We then review results from applying the new detectors, known as alternate null hypothesis correlation (ANCorr) templates, to a test set of seismic waveforms. Finally, we compare ANCorr results with those from application of the conventional waveform correlation matched filter technique.

58 GEOSCIENCES↗

Explosion Discrimination Using Seismic Gradiometry and Spectral Filtering of Data

Here, we present a new method to discriminate between earthquakes and buried explosions using observed seismic data. The method is different from previous seismic discrimination algorithms in two main ways. First, we use seismic spatial gradients, as well as the wave attributes estimated from them (referred to as gradiometric attributes), rather than the conventional three-component seismograms recorded on a distributed array. The primary advantage of this is that a gradiometer is only a fraction of a wavelength in aperture compared with a conventional seismic array or network. Second, we use the gradiometric attributes as input data into a machine learning algorithm. The resulting discrimination algorithm uses the norms of truncated principal components obtained from the gradiometric data to distinguish the two classes of seismic events. Using high-fidelity synthetic data, we show that the data and gradiometric attributes recorded by a single seismic gradiometer performs as well as a conventional distributed array at the event type discrimination task.

58 GEOSCIENCES↗

Relocation of the Seismicity of the Caucasus Region

Our objective is to improve the view of the seismicity in the Caucasus region using instrumental data between 1951 and 2019. To create a comprehensive catalog, we combine the bulletins of local agencies and the International Seismological Centre, and use an advanced single-event location algorithm, iLoc, to obtain better locations. We show that relocations with iLoc, using travel-time predictions from the 3D upper mantle velocity model, Regional Seismic Travel Time, improve the locations. Then, using the iLoc results as initial locations and the ground-truth events identified in the iLoc results as fix points, we apply Bayesloc, a multiple-event location algorithm, to simultaneously relocate the entire seismicity of the Caucasus region. Here we demonstrate that the simultaneous relocation of the seismicity with Bayesloc clarifies the location and geometry of major active structures accommodating ongoing convergence between the Arabian and Eurasian continents between the Black and Caspian Seas. Among our major findings is the confirmation of widespread seismicity in the mantle beneath the northern flank of the Greater Caucasus and central Caspian, resulting from north-dipping subduction of the Kura and South Caspian basins and the identification of a discrete band of crustal seismicity beneath the southern flank of the Greater Caucasus.

58 GEOSCIENCES↗

Multiwell Fiber Optic Sensing Reveals Effects of CO2 Flow on Triggered Seismicity

Induced seismicity is one of the main risks for gigaton-scale geological storage of carbon dioxide (⁠CO 2 ⁠). Thus, passive seismic monitoring is often recommended as a necessary component of the monitoring systems for CO 2 storage projects, with a particular forcus on risk mitigation. We present the first field study, CO2CRC Otway Project Stage 3 (Victoria, Australia), where distributed acoustic sensing (DAS) enabled high-precision tracking of the induced seismicity triggered by a small CO 2 injection and also informed the reservoir models. In 610 days of passive seismic monitoring of the Stage 3 injection, we detected 17 microseismic events (maximum moment magnitude Mw 0.1) using five deep boreholes equipped with enhanced-sensitivity optical fiber. The DAS array has sensitivity sufficient for detection and location of induced events with M w ~ ₋2 in a monitoring borehole located up to 1500 m away. Thanks to the dense spatial sampling by the DAS, we were able to estimate the focal mechanisms for events with M w > ₋1.5; although the monitoring boreholes provided very limited angular coverage. The main cluster of the events has the same location and source mechanism as the one triggered by the previous CO 2 injection at the Otway Project site, Stage 2C. Surprizingly, the Stage 2C and Stage 3 events closely followed the actual movement of the CO 2 saturation plume front (not the pressure front), as observed using controlled-source reflection seismic images. The nature of the plume-fault interaction remains unclear, but some alteration of the fault gouge by CO 2 might be responsible for the faults’ reactivation by the pressure perturbation. Importantly, the seismogenic fault could not be identified in the seismic images and was only revealed by DAS observations, which also demonstrated the signature of fluid–rock interaction, that may control the CO 2 flow.

58 GEOSCIENCES↗

Seafloor Seismic Noise Patterns Across the Pacific Basin

Seismic hazard monitoring and global tomography efforts are improved by recording signals at a variety of distances and azimuths to maximize subsurface sampling. Although seismic networks provide good to excellent coverage on land, seafloor stations are still sparse. Inclusion of ocean-based data would greatly improve the global coverage of seismic networks, but the use of seafloor seismic data to complement land-based detection and characterization of events is complicated by the generally much higher ambient noise level in the ocean compared to that observed on land. This noise is driven primarily by sea surface waves and tides, but how seismic noise levels vary with location in the oceans is not well described. Here, in this work, we analyze the relationship between ocean surface wave height and seismic noise in the 0.4–4 Hz frequency band at ocean-bottom seismometer deployments across the Pacific basin. We find that a noise-to-responsiveness ratio (NRR)—the median noise level at a station divided by its sea surface wave height responsiveness—correlates negatively with detection success for large teleseismic earthquakes. Stations that are close to land, with relatively shallow ocean and low wind speed, often have lower NRR than open-ocean stations, but the connection between geographic location and earthquake detection success is imperfect.

58 GEOSCIENCES↗

Final report of a scalable, automated, semipermanent seismic array (SASSA) method for detecting CO 2 extent during geologic CO 2 injection

The scalable, automated, semipermanent seismic array (SASSA) project led and managed by the Energy & Environmental Research Center (EERC) was designed as a 3-year proof-of-concept study to evaluate and demonstrate an innovative application of the seismic method. The concept was to use a sparse surface array of 96 nodal seismic sensors paired with a single, remotely operated active seismic source at a fixed location to monitor for CO 2 saturation changes in a subsurface reservoir by processing the data for time-lapse changes at individual, strategically chosen reservoir reflection points. The combination of autonomous equipment and modern processing algorithms was used to apply the seismic method in a manner different from the normal paradigm of collecting a spatially dense data set to produce an image. It was used instead to monitor individual, strategically chosen reservoir reflection points for detectable signal character changes that could be attributed to the passing of a CO 2 saturation front or, possibly, changes in reservoir pressure. Data collection occurred over the course of 1 year at an oil field undergoing CO 2 injection for enhanced oil recovery (EOR) and focused on four overlapping “five-spot” EOR injector–producer patterns. Selection, procurement, configuration, installation, and testing of project equipment and collection of five baseline data sets were completed in advance of CO 2 injection within the study area. Weekly remote data collection produced 41 incremental time-lapse records for each of the 96 nodes. Validation was provided by two methods: 1) a conventional 2-D seismic line acquired through the center of the study area before injection started and again after the project ended and processed in a time-lapse manner and 2) by CO 2 saturation maps created from reservoir simulations based on injection and production history matching. Interpreted results were encouraging but mixed, with indications of changes likely due to the presence of CO 2 on some node reflection points where and when effects would be expected and noneffects where no CO 2 was expected, while results at some locations where simulation outputs suggested CO 2 should be present were ambiguous. Acquisition noise impacted interpretation of data at several locations. Many lessons learned were generated by the study to inform and improve results on a follow-up study. The ultimate aim of the project was to evaluate whether deployment of a SASSA technology can provide a useful and cost-effective monitoring solution for future CO 2 injection projects. The answer appears to be affirmative, with the expectation that lessons learned applied to future iterations, together with technology advances, will likely result in significant improvements.

58 GEOSCIENCES↗

Seismic Phase Identification with Speech Recognition Algorithms

Seismic signals are composed of the seismic waves (phases) that reach a sensor, similar to the way speech signals are composed of phonemes that reach a listener's ear. Large/small seismic events near/far from a sensor are similar to loud/quiet speakers with high/low-pitched voices. We leverage ideas from speech recognition for the classification of seismic phases at a seismic sensor. Seismic Phase ID is challenging due to the varying paths and distances an event takes to reach a sensor, but there is consistent structure of the makeup (e.g. ordering) of the different phases arriving at the sensor.

58 GEOSCIENCES↗

Imaging Fracture Networks Using Joint Seismic and Electrical Change Detection Techniques (Final Report)

During the initial phase of this SubTER project, we conducted a series of high resolution seismic imaging campaigns designed to characterize induced fractures. Fractures were emplaced using a novel explosive source, designed at Sandia National Laboratories, that limits damage to the borehole. This work provided evidence that fracture locations could be imaged at inch scales using high-frequency seismic tomography but left many fracture properties (i.e. permeability) unresolved. We present here the results of the second phase of the project, where we developed and demonstrated emerging seismic and electrical geophysical imaging technologies that characterize 1) the 3D extent and distribution of fractures stimulated from the explosive source, 2) 3D fluid transport within the stimulated fracture network through use of a contrasting tracer, and 3) fracture attributes through advanced data analysis. Focus was placed upon advancing these technologies toward near real-time acquisition and processing in order to help provide the feedback mechanism necessary to understand and control fracture stimulation and fluid flow. Results from this study include a comprehensive set of 4D cross-hole seismic and electrical data that take advantage of change detection methodologies allowing for perturbations associated with the fracture emplacement and particulate tracer to be isolated. During the testing the team also demonstrated near real-time 4D electrical resistivity tomography imaging and 4D seismic tomography using the CASSM approach with a temporal resolution approaching 1 minute. All of the data collected were used to develop methods of estimating fracture attributes from seismic data, develop methods of assimilating disparate and transient data sets to improve fracture network imaging resolution, and advance capabilities for near real-time inversion of cross-hole tomographic data. These results are illustrated here. Advancements in these areas are relevant to all situations where fracture emplacement is used for reservoir stimulation (e.g. Enhanced Geothermal Systems (EGS) and tight shale gases).

42 ENGINEERING↗

Seismic Spatial Gradients and Machine Learning-Based Classifiers for Explosion Monitoring (LDRD 218327)

This final report summarizes the work completed under the Laboratory Directed Research and Development (LDRD) project “Seismic Spatial Gradients as a Machine Learning-Based Classifier for Explosion Monitoring.” The overarching goal of the project was to explore the efficacy of using machine learning-based classification algorithms where the input data are the spatial gradient of the seismic wavefield collected at a single point on the Earth’s surface. The methods that I describe here are in direct contrast to conventional methods of seismic discrimination which typically rely on a spatially extended network of instruments and physics-based wavefield attributes such as, for example, the ratio between $\textit{P}$ and $\textit{S}$ waves. Rather, we use the spatial gradient of the seismic wavefield observed at a single point on the Earth’s surface and data processing approaches inspired by the machine learning community. We tested two algorithms, a neural network and a modified version of principal component analysis termed Spectrally Filtered Principal Component Analysis (SFPCA). To test these algorithms, we first conducted a series of numerical tests using synthetic data and then conducted a small-scale controlled field experiment. The tests using synthetic data showed that both algorithms had high success rates on gradiometric data, even when simulated noise was added to the signal. Furthermore, we found that using seismic spatial gradients increased the performance of our discrimination algorithms when compared to using just the traditional translational motion seismic data. The tests with field data also showed a high degree of discriminative success.

58 GEOSCIENCES↗

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↗

How far can it go?: Determining the reach of induced seismicity

Induced seismicity is a unique societal and intellectual challenge in the solid Earth sciences. Societally, induced seismicity provides an opportunity to harness the advances of earthquake science to directly reduce hazard. Intellectually, induced seismicity provides a partially controlled, field-scale experiment that allows us to investigate otherwise inaccessible processes. The dramatic increase in induced earthquakes in the 21st century requires rising to both aspects of the challenge. A key issue is establishing general trends of induced seismicity. Fluid injection activity has led to several recent moderate magnitude earthquakes and many detailed studies are available (Bao and Eaton, 2016; Ellsworth, 2013; Kim et al., 2018; Weingarten et al., 2015; Yeck et al., 2017). However, a generalized understanding has been difficult to reach in part because of the variety of potential observables and variables. Some efforts at generalizing behavior have focused on magnitudes or seismicity rate (Bachmann et al., 2011; Galis et al., 2017; Gischig and Wiemer, 2013; McGarr, 1976, 2014; Shapiro et al., 2010). Others have examined migratory behavior (Bourouis and Bernard, 2007; Shapiro and Dinske, 2009; Shapiro et al., 2002). We expand on these results and focus on the spatial distribution of induced events, which has rarely been examined quantitatively. We find that the spatial decay of earthquakes from injection wells is a particularly powerful tool for accessing physics and comparing the factors that exert the most significant control on the resulting seismicity (Goebel and Brodsky, 2018).

15 GEOTHERMAL ENERGY↗

Assessing the Value of Seismic Amplitude Versus Offset (AVO) Attributes for CO2 Storage Project Using a Bayesian Network Model for Decision Support

Attributes versus offset (AVO) are a set of measurements to analyze how the characteristics of reflected seismic waves change as a function of the offset. It can be useful for monitoring CO2 storage sites because the presence of leaked CO2 into the overlying aquifer can change the properties of the rocks and pore fluids composition that can alter the way seismic waves reflect and their amplitudes. The time-lapse changes in AVO attributes derived from repeat seismic surveys can help identify anomalies or shifts in the subsurface that could potentially be used as an indicator for CO2 leak detection. Our study leverages multiple seismic attributes derived from synthetic seismic data and Bayesian network model to quantify the probability of leak detection in the overlying aquifer above the storage reservoir. It helps to quantify the value of individual seismic attributes at multiple monitoring periods based upon their sensitivities.

Kumar, Abhash↗

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↗

SHALLOW MARINE HIGH-RESOLUTION 3D SEISMICS FOR ABOVE- ZONE CO2MONITORING: 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 E.↗