Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “UNDERGROUND EXPLOSIONS”

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 OSTI2024

Fracture Network Influence on Rock Damage and Gas Transport following an Underground Explosion

Simulations of rock damage and gas transport following underground explosions that omit preexisting fracture networks in the subsurface cannot fully characterize the influence of geo-structural variability on gas transport. Previous studies do not consider the impact that fracture network structure and variability have on gas seepage. In this study, we develop a sequentially coupled, axi-symmetric model to look at the damage pattern and resulting gas breakthrough curves following an underground explosion given different fracture network realizations. We simulate 0.327 and 0.164 kT chemical explosives with burial depths of 100 m for 90 stochastically generated fracture networks. Gases quickly reach the surface in 30% of the higher yield simulations and 5% of the lower yield simulations. The fast breakthrough can be attributed to the formation of connected pathways between fractures to the surface. The formation of a connected damage pathway to the surface is not clearly correlated with the fracture intensity (P32) in our simulations. Breakthrough curves with slower transport are highly variable depending on the fracture network sample. The variability in the breakthrough behavior indicates that ignoring the influence of fracture networks on rock damage, which strongly influences the hydraulic properties following an underground explosion, will likely lead to a large underestimation of the uncertainty in the gas transport to the surface. This work highlights the need for incorporation of fracture networks into models for accurately predicting gas seepage following underground explosions.

58 GEOSCIENCES↗
DOE OSTI2017

Two Sources of Nonisotropic Radiation From Underground Explosions in Granite

Significant tangential ground motion observed during underground explosions makes it difficult to distinguish them from natural earthquakes. Such motion can be generated by the source geometry and emplacement conditions, by the heterogeneous nature of the rock mass (mechanical properties may vary in space due to the presence of cracks, joints, faults, and various geologic layers) and also by the nonuniform in situ stress state. The last mechanism is increasingly important with depth when the difference in main principal stresses becomes significant. This paper is focused on the role of material strength of the rock mass in generation of nonradial motion during explosions in prestressed media. Numerical modeling of underground chemical explosions in granite at various depths has been conducted to compare two possible mechanisms of shear wave generation. The first, caused by rock mass anisotropy, is important at shallow depth. The second is related to elastic-plastic relaxation around the cavity created by the explosion. As a result, tangential motions for these two mechanisms have different signatures.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The influence of physical and algorithmic factors on simulated far-field waveforms and source–time functions of underground explosions using unsupervised machine learning

SUMMARY Characterizing explosion sources and differentiating between earthquake and underground explosions using distributed seismic networks becomes non-trivial when explosions are detonated in cavities or heterogeneous ground material. Moreover, there is little understanding of how changes in subsurface physical properties affect the far-field waveforms we record and use to infer information about the source. Simulations of underground explosions and the resultant ground motions can be a powerful tool to systematically explore how different subsurface properties affect far-field waveform features, but there are added variables that arise from how we choose to model the explosions that can confound interpretation. To assess how both subsurface properties and algorithmic choices affect the seismic wavefield and the estimated source functions, we ran a series of 2-D axisymmetric non-linear numerical explosion experiments and wave propagation simulations that explore a wide array of parameters. We then inverted the synthetic far-field waveform data using a linear inversion scheme to estimate source–time functions (STFs) for each simulation case. We applied principal component analysis (PCA), an unsupervised machine learning method, to both the far-field waveforms and STFs to identify the most important factors that control variance in the waveform data and differences between cases. For the far-field waveforms, the largest variance occurs in the shallower radial receiver channels in the 0–50 Hz frequency band. For the STFs, both peak amplitude and rise times across different frequencies contribute to the variance. We find that the ground equation of state (i.e. lithology and rheology) and the explosion emplacement conditions (i.e. tamped versus cavity) have the greatest effect on the variance of the far-field waveforms and STFs, with the ground yield strength and fracture pressure being secondary factors. Differences in the PCA results between the far-field waveforms and STFs could possibly be due to near-field non-linearities of the source that are not accounted for in the estimation of STFs and could be associated with yield strength, fracture pressure, cavity radius and cavity shape parameters. Other algorithmic parameters are found to be less important and cause less variance in both the far-field waveforms and STFs, meaning algorithmic choices in how we model explosions are less important, which is encouraging for the further use of explosion simulations to study how physical Earth properties affect seismic waveform features and estimated STFs.

58 GEOSCIENCES↗
DOE OSTI2020

Seismoacoustic analysis of underground explosions using the Rayleigh integral

SUMMARY A coupled seismoacoustic model is developed for the analysis of acoustic signals produced by underground explosive events with an aim to develop a means of improving estimated depth and yield for explosion monitoring. A ground spall model is used to predict surface motion characteristics produced by an underground explosion and the Rayleigh integral is applied to relate the surface motion to the acoustic signal some distance from surface ground zero. The low-frequency component of the ground motion associated with the prolonged free fall of lofted material during spall is found to dominate the acoustic signal propagating away from surface ground zero at shallow angles. The model is applied to observed ground motion and acoustic signals recorded during the Source Physics Experiment (SPE) with promising results. In addition to accurately predicting characteristics of the observed acoustic signals from several of the SPE events, the model provides a means of explaining the lack of signals observed during several events in the SPE due to the directionality of the higher frequency acoustic signals associated with the uplift and closure components as well as the lack of a spatially localized, longer duration dwell.

58 GEOSCIENCES↗
DOE OSTI2020

Low-frequency Electromagnetic Detection Limits of Underground Nuclear Explosions

Underground nuclear explosions produce a broadband electromagnetic signal that can propagate to the surface at some distance. We perform a logistic regression on historical electromagnetic detections to calculate a probabilistic detection curve for prediction of the propagation distance for any given underground nuclear explosion. The curve predicts positive detections at scaled ranges of 2414.6, 2025.7, 1908.3, and 1672.4 m/ktonne 1/3 at probabilities of 50, 90, 95 and 99%, respectively. For a theoretical 100 ktonne underground nuclear explosion these scaled ranges translate to absolute ranges of 11.2, 9.4, 8.9, and 7.8 km at probabilities of 50, 90, 95 and 99%, respectively. We also motivate a need for improved electromagnetic sensor design and deployment by extending the detection analysis using the quoted signal-to-noise ratio of two and noise amplitude of 10 pT along with attenuation relationships that go as 1/r to 1/r 3 . We find that reducing the noise floor by a factor often could extend the detection limit to approximately 5 to 25 km/ktonne 1/3 , which translates to 24 to 112 km for a theoretical 100 ktonne underground nuclear explosion.

58 GEOSCIENCES↗
DOE OSTI2017

Detecting Surface Changes from an Underground Explosion in Granite Using Unmanned Aerial System Photogrammetry

Efficient detection and high-fidelity quantification of surface changes resulting from underground activities are important national and global security efforts. In this investigation, a team performed field-based topographic characterization by gathering high-quality photographs at very low altitudes from an unmanned aerial system (UAS)-borne camera platform. The data collection occurred shortly before and after a controlled underground chemical explosion as part of the United States Department of Energy’s Source Physics Experiments (SPE-5) series. The high-resolution overlapping photographs were used to create 3D photogrammetric models of the site, which then served to map changes in the landscape down to 1-cm-scale. Separate models were created for two areas, herein referred to as the test table grid region and the nearfield grid region. The test table grid includes the region within ~40 m from surface ground zero, with photographs collected at a flight altitude of 8.5 m above ground level (AGL). The near-field grid area covered a broader area, 90–130 m from surface ground zero, and collected at a flight altitude of 22 m AGL. The photographs, processed using Agisoft Photoscan® in conjunction with 125 surveyed ground control point targets, yielded a 6-mm pixel-size digital elevation model (DEM) for the test table grid region. This provided the ≤3 cm resolution in the topographic data to map in fine detail a suite of features related to the underground explosion: uplift, subsidence, surface fractures, and morphological change detection. The near-field grid region data collection resulted in a 2-cm pixel-size DEM, enabling mapping of a broader range of features related to the explosion, including: uplift and subsidence, rock fall, and slope sloughing. This study represents one of the first works to constrain, both temporally and spatially, explosion-related surface damage using a UAS photogrammetric platform; these data will help to advance the science of underground explosion detection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Application of Active Seismic and Electrical Methods to Detect and Characterize Subsurface Effects of an Underground Explosion

During the period of underground nuclear testing in the U.S., now ended for over 20 years, the focus of work was to conduct the test and analyze the data as quickly as possible; there was little concern about application of technologies that could be used to detect and characterize the location and depth of the explosion. Thus, there is little experience with applying geophysical methods in places where an underground nuclear explosion (UNE) has taken place. On the other hand, application of geophysical methods to the detection of shallow artifacts and water table depth is commonly applied in the environmental remediation and site characterization realm of commercial geophysical surveys. For nuclear nonproliferation treaty applications which may incorporate both shallow and deep geophysical probing methods (such as on-site inspection (OSI) under the Comprehensive Nuclear Test Ban Treaty [CTBT]), there thus remains a need to better understand detection capabilities of the deeper targets: anomalies related to an explosion cavity or rubble zone that may be located several hundred meters below the ground surface. The purpose of this Coordinated Inspection Technologies (CIT) project was to gain experience in applying active seismic methods and electrical methods at a site in Nevada where an underground nuclear test (UNE) had been conducted in the past to get a better understanding of the capabilities of these methods for detection and characterization of deep subsurface effects of the explosion.

58 GEOSCIENCES↗
DOE OSTI2014

Modelling infrasound signal generation from two underground explosions at the Source Physics Experiment using the Rayleigh integral

For this study, we use the Rayleigh integral (RI) as an approximation to the Helmholtz–Kirchoff integral to model infrasound generation and propagation from underground chemical explosions at distances of 250 m out to 5 km as part of the Source Physics Experiment (SPE). Using a sparse network of surface accelerometers installed above ground zero, we are able to accurately create synthetic acoustic waveforms and compare them to the observed data. Although the underground explosive sources were designed to be symmetric, the resulting seismic wave at the surface shows an asymmetric propagation pattern that is stronger to the northeast of the borehole. This asymmetric bias may be attributed to the subsurface geology and faulting of the area and is observed in the acoustic waveforms. We compare observed and modelled results from two of the underground SPE tests with a sensitivity study to evaluate the asymmetry observed in the data. This work shows that it is possible to model infrasound signals from underground explosive sources using the RI and that asymmetries observed in the data can be modelled with this technique.

58 GEOSCIENCES↗
DOE OSTI2015

On the generation of horizontal shear waves by underground explosions in jointed rocks

This paper describes computational studies of non-spherical ground motions generated by spherical explosions in a heavily jointed granite formation. Various factors affecting the shear wave generation are considered, including joint spacing, orientation, persistence and properties. Simulations are performed both in 2D for a single joint set to elucidate the basic response mechanisms, and in 3D for multiple joint sets to realistically represent in situ conditions in a realistic geologic setting. The joints are modeled explicitly using both contact elements and weakness planes in the material. Simulations are performed both deterministically and stochastically to quantify the effects of geologic uncertainties on near field ground motions. The mechanical properties of the rock and the joints as well as the joint spacing and orientation are taken from experimental test data and geophysical logs corresponding to the Climax Stock granitic outcrop, which is the geologic setting of the Source Physics Experiment (SPE). Agreement between simulation results and near field wave motion data from SPE enables newfound understanding of the origin and extent of non-spherical motions associated with underground explosions in fractured geologic media.

58 GEOSCIENCES↗
DOE OSTI2022

Combining translational and rotational seismic motions to invert local-scale seismic data for time-variable moment tensors: do rotational motions help for high-frequency seismic data produced by underground explosions?

SUMMARY We present an analysis of combining translational and rotational seismic data in an inversion for the time-variable source time functions corresponding to the components of the seismic moment tensor. We conduct a series of numerical experiments where the data are simulated by a combination of an underground explosion and a co-located double couple shear source and recorded on surface-mounted seismometers within 1–2 km of the source. The experiments are designed to mimic explosion seismology experiments, and thus the data are in the 1–10 Hz frequency range and contain very few surface waves. We use a Monte Carlo method to propagate Earth model uncertainty into the estimates of seismic source parameters. In our experiments, we find that the uncertainty of the estimated seismic source parameters increases when we add rotational seismic motions to the inversion when using a constant number of data channels. In this case, the increased degree of uncertainty in the final results is most likely due to the near-surface Earth model uncertainty that we introduce in our simulations. However, for a fixed number of seismic stations, adding rotational seismic motions to the inversion acts to decrease the uncertainty of the estimated seismic source parameters, most likely due to the increase in the number of data channels used in the inversion.

Poppeliers, Christian (ORCID:0000000159526849)↗
DOE OSTI2020

Free-field Ground Motion Induced by Underground Explosions at Aqueduct Mesa with Predictions for Physical Experiment One (PE1)

We fit standard power-law models for ground acceleration and particle velocity using scaled ground motion observations from underground nuclear explosions at Aqueduct Mesa on the Nevada National Security Site. The models are then validated using small-scale high-explosive (HE) test data from nearby Rainier Mesa. The comparison between model predictions and observations is good and can be made more favorable if the assumed yield for the HE tests are doubled. Further validation is made with a numerical experiment that additionally shows a transition in ground motion attenuation to elastic propagation (1 = r). An extended model is provided to incorporate the long-range transition in attenuation, which results in increased ground motions at farther range. We then make predictions for the field experiment, Physical Experiment One (PE1), to take place in P-tunnel on Aqueduct Mesa, and suggest that the models could be used in planning various operations for PE1.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2020

The Complicated Link Between Material Properties and Microfracture Density for an Underground Explosion in Granite

Geophysical techniques are often implemented as quick and inexpensive ways to locate and characterize fractures in the subsurface, which is important for a number of geoscience fields. Seismic velocities are the most widely used proxies for identification of fractures, but the correlation is not always well-defined. In this study we present material property data: unconfined compressive strength (UCS), bulk density (ρ), Young's modulus (E), Poisson's ratio (ν), P wave velocity (V p ), and S wave velocity (V s ), in conjunction with microfracture densities measured on samples of granite collected before and after underground chemical explosions. Results indicate the relationship between fractures and material properties is complex, even in this single-lithology environment. We interpret that this complexity arises from varying fracture mechanisms (e.g. dilation-inducing fractures vs compression-inducing fractures) in different parts of the core, due to differences in stress conditions. Additional complexity may result from chemical interactions between the fresh fractures and the fluids in the area. Water content appears to have a significant, if not dominant, role in the unconfined compressive strength (UCS) of the samples. We suggest caution when using elastic property measurements as a proxy for fracturing in areas of explosion-induced damage, or in other areas where a variety of mechanisms induce fracturing.

58 GEOSCIENCES↗
DOE OSTI2014

Source Characterization of Underground Explosions from Combined Regional Moment Tensor and First-Motion Analysis

Here in this study, we investigate the 14 September 1988 U.S.–Soviet Joint Verification Experiment nuclear test at the Semipalatinsk test site in eastern Kazakhstan and two nuclear explosions conducted less than 10 years later at the Chinese Lop Nor test site. These events were very sparsely recorded by stations located within 1600 km, and in each case only three or four stations were available in the regional distance range. We have utilized a regional distance seismic waveform method fitting long-period, complete, three-component waveforms jointly with first-motion observations from regional stations and teleseismic arrays. The combination of long-period waveforms and first-motion observations provides a unique discrimination of these sparsely recorded events in the context of the Hudson et al. (1989) source-type diagram. We demonstrate through a series of jackknife tests and sensitivity analyses that the source type of the explosions is well constrained. One event, a 1996 Lop Nor shaft explosion, displays large Love waves and possibly reversed Rayleigh waves at one station, indicative of a large F-factor. We show the combination of long-period waveforms and P-wave first motions are able to discriminate this event as explosion-like and distinct from earthquakes and collapses. We further demonstrate the behavior of network sensitivity solutions for models of tectonic release and spall-based tensile damage over a range of F-factors and K-factors.

58 GEOSCIENCES↗
DOE OSTI2023

Blue Canyon Dome: Development of a Small-Scale Testbed for Monitoring Underground Explosions

This report documents the development of the Blue Canyon Dome (BCD) testbed, including test site selection, development, instrumentation, and logistical considerations. The BCD testbed was designed for small-scale explosive tests (~5 kg TNT equivalence maximum) for the purpose of comparing diagnostic signals from different types of explosives, the assumption being that different chemical explosives would generate different signatures on geophysical and other monitoring tools. The BCD testbed is located at the Energetic Materials Research and Testing Center near Socorro, New Mexico. Instrumentation includes an electrical resistivity tomography array, geophones, distributed acoustic sensing, gas samplers, distributed temperature sensing, pressure transducers, and high-speed cameras. This SAND report is a reference for BCD testbed development that can be cited in future publications.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗