Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic arrays”

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 OSTI2021

Seismic Array Analysis Using Fiber-Optic Distributed Acoustic Sensing on Small Local and Regional Earthquakes

We examined the feasibility of using fiber-optic Distributed Acoustic Sensing (DAS) as an alternative to a traditional seismic array. Seismic arrays are used routinely globally to measure waveform propagation parameters and signal features including back-azimuth (BAZ) and apparent horizontal velocity through the process of beamforming or frequency-wavenumber (f-k) array analysis with the assumption that the signals travel across the array as a plane-wave. These measurements are useful for identifying signal detections as teleseismic, regional, or local distance seismic phases based on their velocities traveling across the array (e.g., Rost and Thomas, 2001). Signal enhancement is the main benefit of array processesing through the stacking of multiple channels as a phased array. This enhancement increases as the square root of the number of channels relative to a single seismic station (e.g., Rost and Thomas, 2001), which is a factor of 3 from a typical nine-channel array. DAS provides the potential of stacking waveforms from 100’s to 1000’s of channels. However, DAS is a novel technology designed for geophysical exploration and therefore has some limitations which we will explore. We started with the dataset from the PoroTomo project because of its unique experimental layout of colocated geophone array and DAS deployment so the two recording technologies can be compared side-by-side.

58 GEOSCIENCES↗
DOE OSTI2025

STILGAR Seismic Array Data

Two three-component dense seismic arrays were deployed above an active limestone mine in Pleasant Gap, Pennsylvania; one in the fall and one in the spring of 2023. Each campaign lasted approximately four weeks and included 80 to 100 seismic stations with stations spaced between 100 and 400 meters apart. Each campaign also included several small dense seismic arrays. Two types of seismic stations were deployed including FairfieldNodal ZLand 3C All-in-One seismometers and Geospace 3C GS-ONE geophone paired with GSX3 and GSX4 dataloggers. This dataset includes ~ 4 TB of three-component continuously recorded seismic data in the miniseed format data from both campaigns as well as a readme file for the data structure.

58 GEOSCIENCES↗
DOE OSTI2020

Subtask 1.5 – CO2 Injection Monitoring with an Optimized Scalable, Automated, Semipermanent Seismic Array

The scalable, automated, semipermanent seismic array (SASSA) method is a flexible and relatively cost-effective surface geophysical method for regular time-lapse monitoring of the movement of injected carbon dioxide (CO2) in a reservoir for CO2 enhanced oil recovery (EOR) or geologic CO2 storage operations. It has the advantages of a low-environmental-footprint while monitoring regions of a reservoir from the surface without the need for a regular grid distribution of receivers. Automated data collection is possible. As only time-lapse amplitude changes at the reservoir level due to CO2 movement within the reservoir are monitored, the turnaround time to deliver results from the SASSA method can be short, without the need for long, time-consuming data-processing workflows. As data is collected and processed, incremental information can be provided to the field operator. The Energy & Environmental Research Center (EERC) conducted a SASSA field test from September 2018 to November 2020 in a portion of the Bell Creek Field in Montana, which implemented new CO2 EOR field activities during the study period. Lessons learned from a proof-of-concept study were incorporated to improve the data quality of the SASSA method and demonstrate the viability of the technology. The EERC implemented several enhancements to improve data quality, including 1) an iterative survey design, which allowed placing the receivers in strategic locations where the movement of the CO2 in the reservoir could be tracked with minimum interference by the cultural noise in the study area; 2) the use of powerful seismic sources in the form of surface orbital vibrators, and 3) data acquisition during optimal periods. History-matched reservoir simulation was performed to predict gas saturation and pressure response induced by CO2 injection in the study area. The results were compared with the SASSA-measured responses to CO2 injection as a partial validation technique. A match between the two methods was observed for most of the SASSA points predicted to have intersected a CO2 saturation change. The validated results provide confidence that the SASSA method can be used independently as a CO2 saturation monitoring technique. As data are collected and processed, incremental information can be provided to the field operator. The critical components of the SASSA workflow for a successful application of the method are the following: Iterative survey design with information about CO2 injection activities from the oilfield operator. A detailed CO2 injection plan is the key driver to select the strategic monitoring location of the SASSA sensors. After this information is incorporated in the initial distribution of sources and receivers in the study area, high-resolution satellite images are used to identify ground locations not affected by cultural noise sources, such as power lines, pipelines/flow lines, or roadways. In the next iteration of the survey design, a scouting trip to the study area is needed to understand more details of the noise sources identified in the previous step and the intensity of the field activities that can also generate noise during the monitoring. Integrating the information from the scouting trip into the survey design to select the optimum source and receiver locations is the final step. Noise attenuation. The variety of noise types during seismic monitoring of an oil field is enormous. Tailored noise characterization and processing at a node-by-node level can enhance the performance and sensitivity of the SASSA technique. Future advancements that could improve the efficiency and application of the SASSA technology include: Gaining a better understanding of the noise field produced by the seismic source to aid the choice of receiver location. Surface noise from the source can overwhelm the small signal changes due to CO2 that the SASSA method measures. Improved data-processing workflow to automatically analyze and adapt to dynamic noise conditions associated with industrial settings. This subtask was funded through the EERC–DOE Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE- FE0024233.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Imaging the Shallow Structure of the Yucca Flat at the Source Physics Experiment Phase II Site with Horizontal-to-Vertical Spectral Ratio Inversion and a Large- N Seismic Array

The Source Physics Experiment (SPE) is a series of chemical explosions at the Nevada National Security Site (NNSS) with the goal of understanding seismic-wave generation and propagation of underground explosions. To understand explosion source physics, accurate geophysical models of the SPE site are needed. Here, we utilize a large-N seismic array deployed at the SPE phase II site to generate a shallow subsurface model of shear-wave velocity. The deployment consists of 500 geophones and covers an area of, approximately, 2.5 × 2 km. The array is located in the Yucca Flat in the northeast corner of the NNSS, Nye County, Nevada. Using ambient-noise recordings throughout the large-N seismic array, we calculate horizontal-to-vertical spectral ratios (HVSRs) across the array. We obtain 2D seismic images of shear-wave velocities across the SPE phase II site for the shallow structure of the basin. In this work the results clearly image two significant seismic impedance interfaces at ~150–500 and ~350–600 m depth. The shallower interface relates to the contrast between Quaternary alluvium and Tertiary volcanic rocks. The deeper interface relates to the contrast between Tertiary volcanic rocks and the Paleozoic bedrock. The 2D subsurface models support and extend previous understanding of the structure of the SPE phase II site. This study shows that the HVSR method in conjunction with a large-N seismic array is a quick and effective method for investigating shallow structures.

58 GEOSCIENCES↗
DOE OSTI2023

Internal Structure of the Central Garlock Fault Zone From Ridgecrest Aftershocks Recorded by Dense Linear Seismic Arrays

Abstract We provide high‐resolution seismic imaging of the central Garlock fault using data recorded by two dense seismic arrays that cross the Ridgecrest rupture zone (B4) and the Garlock fault (A5). Analyses of fault zone head waves and P ‐wave delay times at array A5 show that the Garlock fault is a sharp bimaterial interface with P waves traveling ∼5% faster in the northern crustal block. The across‐fault velocity contrast agrees with regional tomography models and generates clear P ‐wave reflections in waveforms recorded by array B4. Kirchhoff migration of the reflected waves indicates a near‐vertical fault between 2 and 6 km depth. The P ‐wave delay times imply a ∼300‐m‐wide transition zone near the Garlock fault surface trace beneath array A5, offset to the side with faster velocities. The results provide important constraints for derivations of earthquake properties, simulations of ruptures and ground motion, and future imaging studies associated with the Garlock fault.

58 GEOSCIENCES↗
DOE OSTI2023

Seismoacoustic Analysis of the Large Surface Explosion Coupling Experiment Using a Large- N Seismic Array

Seismoacoustic signals at local distance (<~10 km) are widely used as important constraints on source parameters for near-surface events, yet the seismoacoustic wave generation and energy partitioning are not fully understood. Spatially dense sensors could provide observations in high resolution to capture the full wavefield for better understanding wave propagation and improving source estimation. Recently, spatially dense observations of the local seismoacoustic wavefield produced by a pair of 1-ton surface explosions have been recorded using a large-N seismic array. This large-N array consists of 446 geophones and covers an area of approximately 2×2.5 km 2 . The two surface explosions occur at the same location but at different times with different atmospheric conditions. Both seismic and air–ground coupled acoustic waves from the two surface explosions are well observed. Analyses of signals recorded by the large-N seismic array show different acoustic wave speed and amplitude for the two explosions. A strong spatial variability in acoustic wave speed and amplitude for each explosion is also observed. Finally, the observations suggest the important role of local atmosphere state on wave propagation and source estimation and demonstrate how the use of a large-N capability can improve characterization of the propagation medium and source.

58 GEOSCIENCES↗
DOE OSTI2020

CO 2 Injection Monitoring with a Scalable, Automated Sparse Seismic Array

Conference presentation at American Institute of Chemical Engineers (AIChE) Annual Meeting, Virtual, November 16–20, 2020. The scalable, automated sparse seismic array (SASSA) method has been ideated as a cost-effective alternative or complementary to conventional 3D seismic methods. Based on a sparse array of receivers deployed strategically on Earth’s surface and one or more fixed seismic sources, the SASSA method can be used to study specific portions of the subsurface to detect the movement of CO 2 in a reservoir with lower environmental footprint and shorter turnaround time to deliver results than conventional 3D seismic methods. Energy & Environmental Research Center (EERC) carried out seismic monitoring of an oilfield undergoing CO 2 enhanced oil recovery in Montana, USA, using the SASSA method.

02 PETROLEUM↗
DOE OSTI2021

A Frequency-Domain-Based Algorithm for Detecting Microseismicity Using Dense Surface Seismic Arrays

We propose a new frequency-domain-based algorithm for detecting small-magnitude seismic events using dense surface seismic arrays. Our proposed method takes advantage of the high energy carried by S waves, and approximate known source locations, which are used to rotate the horizontal components to obtain the maximum amplitude. By surrounding the known source area with surface geophones, we achieve a favorable geometry for locating the detected seismic events with the backprojection method. To test our new detection method, we used a dense circular array, consisting of 151 5 Hz three-component geophones, over a 5 km aperture that was in operation at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) in southcentral Utah. We apply the new detection method during a small-scale test injection phase at FORGE, and during an aftershock sequence of an Mw 4.1 earthquake located ~30 km north of the geophone array, within the Black Rock volcanic field. We are able to detect and locate microseismic events (Mw<0) during injections, despite the high level of anthropogenic activity, and several aftershocks that are missing from the regional catalog. By comparing our method with known algorithms that operate both in the time and frequency domain, we show that our proposed method performs better in the case of the FORGE injection monitoring, and equally well for the off-array aftershock sequence. Our new method has the potential to improve microseismic event detections even in extremely noisy environments, and the proposed location scheme serves as a direct discriminant between true and false detections.

Geochemistry & Geophysics↗
DOE OSTI2023

The Variable Continuous Bimaterial Interface in the San Jacinto Fault Zone Revealed by Dense Seismic Array Analysis of Fault Zone Head Waves

Key factors controlling earthquake ruptures include fault geometry, continuity, and seismic velocity structure around the fault. We present a novel tool that better informs deep bimaterial fault geometry embedded in distributed damage and seismicity, associated velocity contrasts across the fault, and their correlations with surface complexities. The method employs fault zone head and direct body waves and is applied to recordings from five spatiotemporally different seismic arrays along the complex San Jacinto fault zone (SJFZ) in southern California. We detect and distinguish these signals based on instantaneous phase coherence and relative energy in a cascading manner from one scale array to another. The analysis reveals a >70-km long continuous bimaterial interface within the SJFZ with several deep northeast dipping fault segments. The northern SJFZ, for instance, locates ~7 km northeast of its surface expression at 18-km depth. P-wave velocity contrasts range from near 0% to >15%, consistent with other bimaterial faults, and differ by a few % depending on fault-array azimuth, implying directional-dependent velocity contrasts. S-wave head waves and velocity contrasts are also imaged for the first time at the southern SJFZ, averaging to 2.9% in agreement with tomography results. The imaged geometry and continuity suggest the SJFZ initiated along remnant tectonic structures and translates to a rupture potential of M > 7.2, i.e., the sizes of its largest paleo-earthquakes. The P and S contrasts, and their ratios, have important implications for earthquake rupture speed, mode, directivity, and frictional heating along the SJFZ and other major faults globally.

58 GEOSCIENCES↗
DOE OSTI2023

Examining Infrasound Propagation at High Spatial Resolution Using a Nodal Seismic Array

Infrasound—acoustic waves in the atmosphere below 20 Hz—is a useful monitoring tool. Topography and atmospheric structure strongly control infrasound propagation, and at common source–receiver distances neither of these effects can be ignored when quantitative source constraints are sought. Detailed spatial measurements of the infrasound wavefield would inform propagation models and improve source estimates. However, the “large-N” deployment strategy now well-known in seismology has not yet been realized for infrasound studies. Here, we use the 900-node seismic array from the 2014 Imaging Magma Under St. Helens (iMUSH) experiment as a proxy for a large-N infrasound network, by leveraging acoustic–seismic coupled arrivals. The active-source component of iMUSH consisted of 23 shallowly buried explosions around Mount Saint Helens volcano; these explosions produced epicentral infrasound recorded on the nodes. We find that the bulk presence of ground-coupled infrasound on the nodes is controlled by wind noise and source–receiver distance, with observed arrivals for eight explosions. Explosions with the most extensive coupling produce complex spatial waveform patterns across the array. These patterns are related to both topographic and atmospheric propagation effects, as well as spatially variable site (coupling) effects. We compare our observations to simple topographic diffraction and high-resolution wind advection models, and full-wave numerical simulations. We find strong spatial correlations between (a) coupled arrival strength and modeled topographic obstruction and (b) coupled arrival time and along-path winds. Our seismoacoustic analyses and results are applicable to other existing and future nodal seismic data sets and can expand the utility of such deployments.

58 GEOSCIENCES↗
DOE OSTI2023

Lateral Variations Across the Southern San Andreas Fault Zone Revealed From Analysis of Traffic Signals at a Dense Seismic Array

Abstract We image the shallow seismic structure across the Southern San Andreas Fault (SSAF) using signals from freight trains and trucks recorded by a dense nodal array, with a linear component perpendicular to SSAF and two 2D subarrays centered on the Banning Fault and Mission Creek Fault (MCF). Particle motion analysis in the frequency band 2–5 Hz shows that the examined traffic sources can be approximated as moving single‐ or multi‐point sources that primarily induce Rayleigh waves. Using several techniques, we resolve strong lateral variations of Rayleigh wave velocities and Q ‐values across the SSAF, including 35% velocity reduction across MCF toward the northeast and strong attenuation around the two fault strands. We further resolve 10% mass density reduction and 45% shear modulus decrease across the MCF. These findings suggest that the MCF is currently the main strand of the SSAF in the area with important implications for seismic hazard assessments.

58 GEOSCIENCES↗
DOE OSTI2020

Internal structure of the San Jacinto fault zone at the Ramona Reservation, north of Anza, California, from dense array seismic data

SUMMARY We image the internal structure of the San Jacinto fault zone (SJFZ) near Anza, California, with seismic data recorded by two dense arrays (RA and RR) from ∼42 000 local and ∼180 teleseismic events occurring between 2012 and 2017. The RA linear array has short aperture (∼470 m long with 12 strong motion sensors) and recorded for the entire analysed time window, whereas the RR is a large three-component nodal array (97 geophones across a ∼2.4 km × 1.4 km area) that operated for about a month in September–October 2016. The SJFZ at the site contains three near-parallel surface traces F1, F2 and F3 from SW to NE that have accommodated several Mw > 6 earthquakes in the past 15 000 yr. Waveform changes in the fault normal direction indicate structural discontinuities that are consistent with the three fault surface traces. Relative slowness from local events and delay time analysis of teleseismic arrivals in the fault normal direction suggest a slower SW side than the NE with a core damage zone between F1 and F2. This core damage zone causes ∼0.05 s delay at stations RR26–31 in the teleseismic P arrivals compared with the SW-most station, and generates both P- and S-type fault zone trapped waves. Inversion of S trapped waves indicates the core damaged structure is ∼100 m wide, ∼4 km deep with a Q value of ∼20 and 40 per cent S-wave velocity reduction compared with bounding rocks. Fault zone head waves observed at stations SW of F3 indicate a local bimaterial interface that separates the locally faster NE block from the broad damage zone in the SW at shallow depth and merges with a deep interface that separates the regionally faster NE block from rocks to the SW with slower velocities at greater depth. The multiscale structural components observed at the site are related to the geological and earthquake rupture history at the site, and provide important information on the preferred NW propagation of earthquake ruptures on the San Jacinto fault.

Geochemistry & Geophysics↗
DOE OSTI2021

A Data-Driven Framework for Automated Detection of Aircraft-Generated Signals in Seismic Array Data Using Machine Learning

Abstract Ground motions associated with aircraft overflights can cover a significant portion of the seismic data collected by shallowly emplaced seismometers, such as new nodal and Distributed Acoustic Sensing systems. This article describes the first published framework for automated detection of aircraft on single channel and multichannel seismic data. The seismic data are converted to spectrograms in a sliding time window and classified as aircraft or nonaircraft in each window using a deep convolutional neural network trained with analyst-labeled data. A majority voting scheme is used to convert the output from the sequence of sliding time windows onto a decision time sequence for each channel and to combine the binary classifications on the decision time sequences across multiple channels. Precision, recall, and F-score are used to quantify the detection performance of the algorithm on nodal data using fourfold time-series cross validation. By applying our framework to data from the Sage Brush Flats nodal array in Southern California, we provide a benchmark performance and demonstrate the advantage of using an array of sensors.

Geochemistry & Geophysics↗
DOE OSTI2020

Dense Seismic Array Study of a Legacy Underground Nuclear Test at the Nevada National Security Site

The complex postdetonation geologic structures that form after an underground nuclear explosion are difficult to constrain because increased heterogeneity around the damage zone affects seismic waves that propagate through the explosion site. Generally, a vertical rubble-filled structure known as a chimney is formed after an underground nuclear explosion that is composed of debris that falls into the subsurface cavity generated by the explosion. Compared with chimneys that collapse fully, leaving a surface crater, partially collapsed chimneys can have remnant subsurface cavities left in place above collapsed rubble. The 1964 nuclear test HADDOCK, conducted at the Nevada test site (now the Nevada National Security Site), formed a partially collapsed chimney with no surface crater. Understanding the subsurface structure of these features has significant national security applications, such as aiding the study of suspected underground nuclear explosions under a treaty verification. In this study, we investigated the subsurface architecture of the HADDOCK legacy nuclear test using hybrid 2D–3D active source seismic reflection and refraction data. The seismic data were acquired using 275 survey shots from the Seismic Hammer (a 13,000 kg weight drop) and 65 survey shots from a smaller accelerated weight drop, both recorded by ~1000 three-component 5 Hz geophones. First-arrival, P-wave tomographic modeling shows a low-velocity anomaly at ~200 m depth, likely an air-filled cavity caused by partial collapse of the rock column into the temporary postdetonation cavity. A high-velocity anomaly between 20 and 60 m depth represents spall-related compaction of the shallow alluvium. Hints of low velocities are also present near the burial depth (~364 m). The reflection seismic data show a prominent subhorizontal reflector at ~300 m depth, a short-curved reflector at ~200 m, and a high-amplitude reflector at ~50 m depth. Comparisons of the reflection sections to synthetic data and borehole stratigraphy suggest that these features correspond to the alluvium–tuff contact, the partial collapse cavity, and the spalled layer, respectively.

58 GEOSCIENCES↗
DOE OSTI2026

Reply to “Comment on ‘The Reduced Detection Rate of Signals That Are Hidden by Earthquakes: Case Studies with Spotlight Detectors That Operate at Seismic Arrays,’ by Joshua D. Carmichael, Brent G. Delbridge, and Richard Alfaro-Diaz” by Paul G. Richards

We reply to a comment that Paul G. Richards (Richards, 2026) directed to an article by Carmichael et al. (2025). Richards expressed concern over two issues. The first issue relates to the size of the explosions discussed in the article. The second relates to the probability that earthquakes and explosions sourced near the same location coincide in time, by chance. We do not address the first issue because the original article does not assign significance to the size of the explosions from the point of view of an experimental party. We do respond to Richards’ comment on the second issue and claim that it can be explained through a clarification about conditional versus joint probabilities. Our present discussion, therefore, describes the difference between the conditional probability that a correlator fails, given that an explosion and earthquake coincidently occur in the same place, and the joint probability that a correlator fails while an earthquake and explosion also coincidentally occur in the same place. This latter, joint probability appears to be rare. We clarify that our original article did not treat the joint probability, only the conditional probability.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗