Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic data”

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 469 records · Page 26

Source Characterization of the Declared North Korean Nuclear Tests From Regional Distance Coda Wave Spectral Ratios

Abstract Seismic observations of underground nuclear explosions provide crucial data on source yield and depth that cannot easily be estimated from other geophysical methods. However, it is difficult to obtain reliable yield estimates for test sites for which we do not have direct seismic calibration experiments. To obtain source information from uncalibrated sites and paths, local and regional seismic records of six, proximal, declared underground nuclear explosions in North Korea are used to compute spectral ratios of narrow‐band waveform envelopes of body‐wave coda that remove path and site effects to reveal precise, relative source moment. The yields of these explosions are obtained from the observed source ratios by simultaneously fitting the classical source model of Mueller and Murphy (1971), https://doi.org/10.1785/bssa0610061675 to all event pairs. The source model provides an impressive fit to the observations considering that the P phase coda derived source spectral ratios did not require a prior knowledge of the source or regionally calibrated corrections to be applied to the data. However, the observed corner frequencies from S wave coda spectral ratios are lower than the source model predictions, but are well fit by models calculated using the corner frequency consistent with the Fisk conjecture. The results presented here provide novel constraints on the spectral distributions of the energy radiated by the sources of the DPRK test series, and allows for an independent evaluation of existing estimated source model parameters.

58 GEOSCIENCES↗

Tracking Seismic Velocity Perturbations at Ridgecrest Using Ballistic Correlation Functions

Abstract We present results based on data of a dense nodal array composed of 147 stations, deployed in 2022 near the epicenter of the 2019 Mw 7.1 Ridgecrest earthquake to investigate characteristics of the seismic wavefields. Through array analyses, we identified two primary components. First, we observed far-field P waves dominating the 0.5–1.2 Hz frequency range, which are likely primarily generated by wind-driven oceanic swell activity. Second, we detected near-field body waves resulting from anthropogenic activities in the frequency range 2–8 Hz. We examined noise correlation functions derived from data of the dense deployment and regional stations to explore fault-zone seismic velocity changes using ballistic arrivals, with a focus on velocity perturbation shortly before and after the Ridgecrest earthquake sequence. Our findings exhibit distinct behavior compared to results obtained through standard coda-wave interferometry. Particularly, we observed a decrease in P-wave travel time on certain station pairs prior to the 2019 earthquake sequence. Supported by detailed investigation of the local seismic wavefields, we interpret the decreasing P-wave travel time as likely caused by a velocity increase away from the fault, possibly related to fluid migration. However, additional information is necessary to verify this hypothesis.

Geochemistry & Geophysics↗

Forecasting induced seismicity in Oklahoma using machine learning methods

Oklahoma earthquakes in the past decade have been mostly associated with wastewater injection. Here we use a machine learning technique—the Random Forest to forecast induced seismicity rate in Oklahoma based on injection-related parameters. We split the data into training (2011.01–2015.05) and test (2015.06–2020.12) periods. The model forecasts seismicity rate during the test period based on input features, including operational parameters (injection rate and pressure), geological information (depth to basement), and modeled pore pressure and poroelastic stress. The results show overall good match with observed seismicity rate (adjusted R 2 of 0.75). The model shows that pore pressure rate and poroelastic stressing rates are the two most important features in forecasting. The absolute values of pore pressure and poroelastic stress, and the injection rate itself, are less important than the stressing rates. These findings further emphasize that temporal changes of stressing rates would lead to significant changes in seismicity rates.

58 GEOSCIENCES↗

Towards Automated/Semiautomated Extraction of Faults from Lidar Data

The Pajarito fault system is a complex zone of deformation and a seismically active region nestled within the Rio Grande rift in north-central New Mexico. Numerous laterally discontinuous faults and associated folds and fractures interact in a manner that has important implications for seismic hazards and risk mitigation. Previous efforts have established a foundation for the location of lineaments and structures in the Pajarito fault system; however, ensuring the completeness of the current lineament mapping is required for identifying areas for field validation, evaluating the potential for future seismic activity, and better understanding fault interaction. Assistance with this fault-mapping task via automated or semiautomated techniques as applied to lidar data over a large area of interest is highly desirable. A proof-of-concept processing flow which transforms lidar point-cloud data into a raster of surficial fault candidates is described and illustrated herein. We report these initial results hold great promise toward achieving our ultimate goal.

58 GEOSCIENCES↗

Time-Lapse Integration at FWU: Fluids, Rock Physics, Numerical Model Integration, and Field Data Comparison

We present the current status of time-lapse seismic integration at the Farnsworth (FWU) CO2 WAG (water-alternating-gas) EOR (Enhanced Oil Recovery) project at Ochiltree County, northwest Texas. As a potential carbon sequestration mechanism, CO2 WAG projects will be subject to some degree of monitoring and verification, either as a regulatory requirement or to qualify for economic incentives. In order to evaluate the viability of time-lapse seismic as a monitoring method the Southwest Partnership (SWP) has conducted time-lapse seismic monitoring at FWU using the 3D Vertical Seismic Profiling (VSP) method. The efficacy of seismic time-lapse depends on a number of key factors, which vary widely from one application to another. Most important among these are the thermophysical properties of the original fluid in place and the displacing fluid, followed by the petrophysical properties of the rock matrix, which together determine the effective elastic properties of the rock fluid system. We present systematic analysis of fluid thermodynamics and resulting thermophysical properties, petrophysics and rock frame elastic properties, and elastic property modeling through fluid substitution using data collected at FWU. These analyses will be framed in realistic scenarios presented by the FWU CO2 WAG development. The resulting fluid/rock physics models will be applied to output from the calibrated FWU compositional reservoir simulation model to forward model the time-lapse seismic response. Modeled results are compared with field time-lapse seismic measurements and strategies for numerical model feedback/update are discussed. While mechanical effects are neglected in the work presented here, complementary parallel studies are underway in which laboratory measurements are introduced to introduce stress dependence of matrix elastic moduli.

4D↗

Failure prediction and estimation of failure parameters

Machine-learning methods and apparatus are disclosed to determine frictional state or other parameters in an earthquake zone or other failing medium, using acoustic emission, seismic waves, or other detectable indicators of microscopic processes. Predictions of future failures are demonstrated in different regimes. A classifier is trained using time series of acoustic emission data along with historic data of frictional state or failure events. In disclosed examples, random forests and gradient boost trees are used, and grid-search or EGO procedures are used for hyperparameter tuning. Once trained, the classifier can be applied to testing or live data in order to assess a frictional state, assess seismic hazard, or make predictions regarding a future failure event. The technology has been developed in a double direct shear apparatus, but can be widely applied to seismic faults, other terrestrial failures, or failures in man-made structures. Variations are disclosed.

Johnson, Paul Allan↗

Detecting an underground tunnel by applying joint traveltime and waveform inversion

Underground tunnel and void detection is a challenging geophysical problem, and many methods have been proposed. Seismic techniques are promising because of the large seismic velocity contrasts between the air-filled void and the surrounding sediment and concrete. We apply a joint seismic traveltime and waveform inversion method to image a buried tunnel with concrete walls and a void space inside. The joint inversion images the top of the concrete tunnel as a high-velocity anomaly and the void space as a low-velocity anomaly. The location of the velocity anomalies predicted by the method agrees with the known location of the tunnel. As a comparison, the stand-alone full waveform inversion is also applied to the data. The first-arrival traveltime tomography shows weak nonlinearity but fails to image the hidden low-velocity layer. Full waveform inversion is able to image complex near-surface structures, but fitting waveforms may not honor the traveltime fit, especially when the data contain noise. Synthetic and real data tests show that the joint inversion method retains the advantages of both traveltime inversion and full waveform inversion and overcomes their respective drawbacks at the same time. The field example shows the joint inversion provides a better reconstruction of the high-velocity feature representing the top of the concrete wall in terms of its magnitude and location.

58 GEOSCIENCES↗

Journey to Time-Variable Moment Tensors through Inversion of Acoustic and Seismoacoustic Data

We explore the capability of acoustic and seismoacoustic datasets to directly resolve a complex, time-variable source consisting of a buried mechanism, represented as a moment tensor, and a spall mechanism, represented as a vertical force at the surface. Traditionally, each component of a resolved moment tensor assumes one underlying source time function, which likely fails to capture the full evolution of a dynamic source, such as an explosion followed by slip on near-source joints or development of spallation. Specifically, we expand previous work to resolve a time-variable moment tensor using single-modality and joint-modality inversion frameworks through analysis of infrasound and seismoacoustic data recorded as part of the Source Physics Experiment Phase II: Dry Alluvium Geology (DAG). We investigate the impact of including signals from seismic-to-air coupling that are local to each infrasound sensor in comparison to mainly atmosphere-propagating acoustic signals, which occur from coupling of the wavefield from the subsurface to the atmosphere directly above the source. Additionally, we assess the ability of our inversion algorithm to fit observed infrasound data using a variety of time-variable source mechanisms. First, we consider the buried moment tensor source alone, which assumes that the determined Green’s functions incorporate effects from spallation or that the impact from spallation is minimal. Second, we examine the estimated buried moment tensor and vertical surface spallation as terms that must both be resolved in the inversion. Third, we assess the ability for an estimated vertical surface spallation source to fit the acoustic data on its own. Finally, we compare results from the joint inversion of both seismic geophone and infrasound acoustic data for the buried-only source compared to buried and spallation sources. Our results are a preliminary investigation into the applications of the inversion technique to recorded datasets and show the technique has limited capabilities using acoustic data alone. Instead, this method shows promise for seismic and seismoacoustic datasets to resolve the time-variable mechanisms of a buried source.

47 OTHER INSTRUMENTATION↗

Utah FORGE Friction-Permeability-Seismicity Laboratory Experiments with Non-Linear Acoustics

Laboratory experimental data on saw-cut interface of Westerly Granite and Utah Forge granitoid rocks. Experiments include velocity-stepping and fluid pressure stepping experiments. Mechanical data from 3 ISCO pumps connected to a Temco pressure vessel measure axial, confining and fault: fluid pressure (kPa), fluid flow rate (mL/min) and volume remaining in pump (mL). Non-linear acoustic data acquired via Verasonics systems connected to PZTs inside the pressure vessel give the timeshift, amplitude and RMS amplitude of the passing P-wave.

15 GEOTHERMAL ENERGY↗

Utah FORGE: Wells 16A(78)-32 and 16B(78)-32 Stimulation Program Report - May 2024

This dataset consists of a comprehensive report documenting the stimulation program conducted in May 2024 on Utah FORGE wells 16A(78)-32 and 16B(78)-32. It includes detailed accounts of operational and scientific objectives, stimulation methodologies, and testing outcomes. The report describes the hydraulic fracturing activities, equipment setups, proppant and fluid usage, and the execution of a nine-hour circulation test. It also includes observations from seismic and fiber optic monitoring systems, providing data on microseismic activity and fracture propagation. Detailed parameters for each stimulation stage are provided, alongside operational challenges and solutions. The dataset includes analyses of injection and production rates, well temperatures, and pressure data, supported by graphical illustrations and logs.

15 GEOTHERMAL ENERGY↗

Report on the recovery of legacy 1997-1998 Depth-of-Burial / Tunnel-and-Hole-Closure explosion data from the Kazakhstan former Soviet nuclear test site, Semipalatinsk.

In 1996, negotiations between the National Nuclear Center of the Republic of Kazakhstan (NNC), the U.S. Defense Special Weapons Administration (DSWA) and the U.S. Department of Energy led to an agreement for a joint series of experiments at the former Soviet nuclear test site Semipalatinsk (STS). These experiments had a dual purpose: to close holes and tunnels originally created for the purpose of testing nuclear weapons, and to characterize explosive sources at a former nuclear test site. This report covers the recovery of legacy data from the experiments conduct in 1997 and 1998. Two teams collected data. Los Alamos National Laboratory (LANL), in collaboration with the NNC, collected the near-source (local) data out to approximately 20 km. Lamont-Doherty Earth Observatory (LDEO), with funding from Lawrence Livermore Laboratory (LLNL), also in collaboration with NNC, collected regional data, 200 – 1000 km distant, leveraging the KZ regional seismic network that LDEO had recently installed in collaboration with NNC. The two teams also worked together to some extent, regarding instrumentation and planning. This report primarily addresses the near-source data collected by the LANL/NNC team.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Toward improved urban earthquake monitoring through deep-learning-based noise suppression

Earthquake monitoring in urban settings is essential but challenging, due to the strong anthropogenic noise inherent to urban seismic recordings. Here, we develop a deep-learning-based denoising algorithm, UrbanDenoiser, to filter out urban seismological noise. UrbanDenoiser strongly suppresses noise relative to the signals, because it was trained using waveform datasets containing rich noise sources from the urban Long Beach dense array and high signal-to-noise ratio (SNR) earthquake signals from the rural San Jacinto dense array. Application to the dense array data and an earthquake sequence in an urban area shows that UrbanDenoiser can increase signal quality and recover signals at an SNR level down to ~0 dB. Earthquake location using our denoised Long Beach data does not support the presence of mantle seismicity beneath Los Angeles but suggests a fault model featuring shallow creep, intermediate locking, and localized stress concentration at the base of the seismogenic zone.

58 GEOSCIENCES↗

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

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

15 GEOTHERMAL ENERGY↗

Preliminary Transfer Learning Results on Israel Data

In this preliminary report, we use publicly available data recorded in Israel to test and expand upon existing machine learning models for seismic-phase detection and arrival-time measurement. We downloaded 3-years of waveform data from Geofon, and cross referenced the waveforms to Israel bulletin picks (Schardong et al., 2021). The initial results using existing models directly generated ubiquitous false detections and that obscured detections of signals that are clearly visible in the waveforms. However, after applying transfer learning (tuning parameters in the existing ML models using one year of the Israel-network data), the results are encouraging, i.e. ML picks agree within a few tenths of a second with bulletin picks and the number of false detections is greatly reduced. The bulletin picks are a good starting point, but they cannot be considered ground-truth. To test potential improvement in picking using ML we would like to relocate the events using the ML picks to see if the events cluster more tightly at known mine locations. However, in order to constrain event locations, we need ML picks for the whole Israeli-Jordanian network, which requires waveforms that are not publicly available.

58 GEOSCIENCES↗

Seismic Features Predict Ground Motions During Repeating Caldera Collapse Sequence

Abstract Applying machine learning to continuous acoustic emissions, signals previously deemed noise, from laboratory faults and slowly slipping subduction‐zone faults, demonstrates hidden signatures are emitted that describe physical details, including fault displacement and friction. However, no evidence currently exists to demonstrate that similar hidden signals occur during seismogenic stick‐slip on earthquake faults—the damaging earthquakes of most societal interest. We show that continuous seismic emissions emitted during the 2018 multi‐month caldera collapse sequence at the Kı̄lauea volcano in Hawai'i contain hidden signatures characterizing the earthquake cycle. Multi‐spectral data features extracted from 30 s intervals of the continuous seismic emission are used to train a gradient boosted tree regression model to predict the GNSS‐derived contemporaneous surface displacement and time‐to‐failure of the upcoming collapse event. This striking result suggests that at least some faults emit such signals and provide a potential path to characterizing the instantaneous and future behavior of earthquake faults.

58 GEOSCIENCES↗

SUBTASK 1.6 – BASIN ELECTRIC CARBON STORAGE RESEARCH PROJECT: NOVEL MONITORING TECHNIQUES

The Energy & Environmental Research Center (EERC) conducted baseline activities associated with an applied research project at Basin Electric Power Cooperative’s (Basin’s) carbon capture and storage (CCS) site in Beulah, North Dakota, to establish novel carbon storage-monitoring techniques as commercial methods under Cooperative Agreement No. DE-FE0024233, Subtask 1.6. The following report summarizes the baseline activities performed and briefly describes the subsequent (operational monitoring) activities that have been proposed to the U.S. Department of Energy (DOE) as part of the overall project to develop and demonstrate novel monitoring techniques at North America’s largest permitted CCS operation. Dakota Gasification Company (DGC), a wholly owned subsidiary of Basin, owns and operates the Great Plains Synfuels Plant (GPSP) approximately 5 miles northwest of the town of Beulah, North Dakota (Figure 1). In 2023, DGC received approval from the North Dakota Industrial Commission (NDIC) to develop a storage facility on-site for injecting a stream of carbon dioxide (CO2) captured from GPSP. DGC will transport the captured CO2 stream with approximately 6.8 miles of transmission lines that extend north of GPSP and inject >1 million tonnes (MMt) of CO2 annually (>1 MMt/yr) over a 12-year period with up to six underground injection control (UIC) Class VI-compliant injection wells completed in the Broom Creek Formation, a predominantly sandstone reservoir and saline aquifer underlying GPSP. The Broom Creek Formation lies approximately 5900 feet (ft) below ground surface (bgs) at GPSP. The commercial scale (i.e., >1 MMt/yr) of DGC’s permitted carbon storage project is ideal for developing and testing the novel monitoring techniques included within Subtask 1.6. The goals of this project are to demonstrate 1) the cost-effectiveness of novel monitoring technologies included as part of this research, 2) technology capability for tracking the CO2 plume and/or associated pressure response in the subsurface and monitoring out-of-zone migration, and 3) compliance with UIC Class VI program requirements. The research activities proposed for the overall project include 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; 5) advanced wellbore-monitoring methods; 6) deployment of an AIM monitoring network; 7) EM monitoring of CO2 with real-time data processing; 8) continued seasonal drone-based surveillance studies; 9) seismic monitoring with passive and active surveys; and 10) wellbore monitoring with nuclear magnetic resonance (NMR) for near-surface characterization. Completion of Activities 1.0–5.0 (baseline activities) are described in this report. Upon authorization of funding by DOE, the EERC will initiate Activities 6.0– 10.0 (operational monitoring activities). Current state-of-the-art (SOA) carbon storage-monitoring techniques require countless labor hours dedicated to the acquisition of data. Once data are gathered, these SOA techniques often rely on commercial facilities to process raw data from the field. However, it is anticipated that next-generation monitoring techniques, such as those being demonstrated, will lower acquisition footprints, be less operationally intensive, and improve data acquisition efficiencies. These new techniques are more conducive to the application of machine learning, artificial intelligence, and automation, thus providing a pathway for integration into active control systems, informing site operability, and improving the integration of data for future CCS projects across the United States. Additionally, reclaimed and active mining lands are present within the project site, creating a unique opportunity to demonstrate the effectiveness of remote sensing and surface-based geophysics monitoring techniques at similar project sites that may include disturbed, unconsolidated, or actively excavated near-surface environments. The efforts included in the overall project will produce necessary designs, learnings, and data acquired during the baseline and operational monitoring periods that are necessary for time-lapse demonstration and validation of the described monitoring techniques. In addition, it is anticipated that the monitoring technologies included in this study will be compliant with UIC Class VI requirements to enable the potential for implementation at other CCS sites across the United States.

42 ENGINEERING↗

Seismotectonic evidence for subduction beneath the Eastern Greater Caucasus

SUMMARY We generated high-resolution 3-D seismic velocity models as well as a relocated earthquake catalogue across the eastern Greater Caucasus and Kura basins. This work was done using data from the recently upgraded Republic Seismological Survey Center's (RSSC) seismic network. We generated our tomographic images of crustal velocity structure in Azerbaijan using double-difference inversions (i.e. tomoDD and hypoDD). Earthquake catalogues from the RSSC between 2011 and 2016 were used; these catalogues include absolute arrival times of 103 288 P- and 120 952 S-wave traveltime picks for 7574 events recorded at 35 stations in Azerbaijan. Beginning with a layered, 1-D velocity model that was estimated using VELEST, we inverted simultaneously for relative location, Vp and Vs on a 3-D grid with dimensions 670 × 445 × 45 km, with a uniform grid spacing of 55 × 55 × 5 km for all of Azerbaijan. We observe that the relocated hypocentres cluster into two depth ranges, at the surface and at depth, that appear to correspond to major fault zones and the top of a subducting plate. Additionally, we note intermediate depth seismicity (∼50–60 km) beneath the Kura Basin, and a northward deepening of earthquake depths. Seismic velocities vary significantly throughout the study region; we observe very slow velocities throughout the Kura Basin between 5 and 15 km, and elevated velocities at 20–35 km. The wholesale velocity structure and seismic structure of Kura Basin strongly mirrors that of the Caspian Sea, which suggests that the geodynamics of the Caspian continue westwards into Azerbaijan. The key results of this study suggest that the northward subduction observed in the Caspian Sea continues beneath the Eastern Greater Caucasus, as well as provides evidence for active faulting along the southern margin of the mountain range.

Gunnels, Michael↗

Local-Regional Surface Wave Dispersion Analysis using Ambient Noise Data of the Southern Nevada National Security Site Area

This report details seismic ambient noise analysis to improve seismic velocity estimates of the southeastern Nevada National Security Site (NNSS). We compare two different methods for estimating surface wave dispersion curves from ambient noise cross-correlations: frequency time analysis (FTAN) and Aki’s cross-spectral method (XSpec). We find that XSpec performs better for our local-regional dataset and frequency content. Using phase velocity estimates from XSpec, we build a preliminary phase velocity dispersion dataset, which contains data for 1,054 station pairs and 8,905 discrete phase velocity measurements for periods between 0.5 and 13 s. This phase velocity dataset will be incorporated into an updated local-regional P and S wave velocity model of the southeastern NNSS in the future and is expected to improve upon shallow velocity estimates.

47 OTHER INSTRUMENTATION↗