Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “explosion monitoring”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

Quantification of spatial and seasonal trends in the atmosphere and construction of statistical models for infrasonic propagation

SUMMARY Infrasonic waves are influenced by variations in the density, pressure and temperature as well as the ambient winds. Modelling infrasonic propagation can be challenging due to the dynamic nature of the atmosphere as well as the sparseness of measurements which result in variability and notable uncertainty. A framework is presented to quantify spatial and seasonal trends in atmospheric structure via analysis of the effective sound speed profile and identification of temporal trends in the middle atmospheric waveguide produced by the circumpolar vortex winds. Seasonal definitions identifying typical atmospheric structures during the summer, winter and spring/fall transition periods are identified using atmospheric data from 2010 through 2020. Seasonal trend analysis is conducted for a number of locations across the contiguous United States to quantify spatial variations in atmospheric structure that impact infrasonic propagation. For each season and location, empirical orthogonal function analysis is used to reduce the historical archive of atmospheric data into a smaller representative set that can be analysed using numerical tools more efficiently. Infrasonic ray tracing and finite-frequency modal propagation analyses are applied to construct propagation path geometry and transmission loss statistics which are useful in localization and yield estimation for infrasonic sources, respectively. An example application is detailed in which transmission loss statistics are combined with an explosive source model and noise statistics to quantify the capability of a network to detect nearby sources.

58 GEOSCIENCES↗

From Livingstone to SMV: Formal Verification for Autonomous Spacecrafts

To fulfill the needs of its deep space exploration program, NASA is actively supporting research and development in autonomy software. However, the reliable and cost-effective development and validation of autonomy systems poses a tough challenge. Traditional scenario-based testing methods fall short because of the combinatorial explosion of possible situations to be analyzed, and formal verification techniques typically require a tedious, manual modelling by formal method experts. This paper presents the application of formal verification techniques in the development of autonomous controllers based on Livingstone, a model-based health-monitoring system that can detect and diagnose anomalies and suggest possible recovery actions. We present a translator that converts the models used by Livingstone into specifications that can be verified with the SMV model checker. The translation frees the Livingstone developer from the tedious conversion of his design to SMV, and isolates him from the technical details of the SMV program. We describe different aspects of the translation and briefly discuss its application to several NASA domains.

Pecheur, Charles↗

Sample Delivery and Computer Control Systems for Detecting Leaks in the Main Engines of the Space Shuttle

The main engines of the Space Shuttle use hydrogen and oxygen as the fuel and oxidant. The explosive and fire hazards associated with these two components pose a serious danger to personnel and equipment. Therefore prior to use the main engines undergo extensive leak tests. Instead of using hazardous gases there tests utilize helium as the tracer element. This results in a need to monitor helium in the ppm level continuously for hours. The major challenge in developing such a low level gas monitor is the sample delivery system. This paper discuss a system developed to meet the requirements while also being mobile. Also shown is the calibration technique, stability, and accuracy results for the system.

Griffin, Timothy P.↗

Contrasting Chemical Kinetic Parameters for Near In‐Service‐Temperature Aging and High‐Temperature Thermal Decomposition of Triaminotrinitrobenzene Formulations

Triaminotrinitrobenzene-based high explosives such as LX-17 offer high energy density and exceptional safety, yet their long-term aging behavior at low temperatures remains poorly understood. In this study, several legacy and new production lots of LX-17 were subjected to accelerated aging experiments below 100°C, during which the formation rate of the initial degradation product, monofurazan (F1), was monitored. Kinetic analysis was performed using a sample-age-aware computational approach, yielding activation energy estimates of 82–91 kJ mol −1 for the low-temperature initiation step—markedly lower than the ∼200 kJ mol −1 associated with high-temperature thermal decomposition. Extrapolation from established cookoff models supports the conclusion that the dominant degradation mechanism at low temperatures differs from that at high temperatures. In conclusion, our results provide a unified kinetic framework that bridges the gap between in-service conditions and high-temperature damage.

Chemistry - Chemical explosives↗

Pantex Plant Ogallala Aquifer and Perched Groundwater Contingency Plan

The Pantex Plant Ogallala Aquifer and Perched Groundwater Contingency Plan has been developed in accordance with the requirements identified in the: • Interagency Agreement for the Pantex Superfund Site, Article 8.5 Work to be Performed, • Compliance Plan Provision of Hazardous Waste Permit No. 50284, and • Record of Decision for Groundwater, Soil, and Associated Media, Pantex Plant. A Long‐Term Monitoring System Design has been designed to monitor conditions in the perched groundwater including changes in the perched aquifer as a result of implementing the response actions. Monitoring is required for verifying the effectiveness of perched groundwater response actions (i.e., conditions in the perched aquifer are being affected as intended) and for confirming that the perched aquifer and Ogallala Aquifer characterization as defined in the Resource Conservation and Recovery Act Facility Investigation Report and the Corrective Measure Studies/Feasibility Study remains accurate. If monitoring results obtained through the monitoring network identify an unexpected condition or deviation, contingent actions will be considered and implemented as necessary to ensure continued protection of the Ogallala Aquifer and human health and the environment. Potential deviations to expected technology performance may be encountered for each of the four primary response actions that compose the selected remedy for perched groundwater; Playa 1 Pump and Treat System, Southeast Area Pump and Treat System, Southeast Area In‐Situ Bioremediation System (comprised of the Southeast In‐Situ Bioremediation System Original System, Southeast Area In‐Situ Bioremediation System Extension System, Offsite In‐Situ Bioremediation System, Perchlorate/Chromium ISB, Northeast ISB and County Road 8 ISB), and Zone 11 In‐Situ Bioremediation System. Monitoring will also be conducted to determine if there are deviations to the expected characterization, e.g., contaminants not expected as a result of the RCRA Facility Investigation characterization. Deviations to expected conditions in the Ogallala Aquifer could also be encountered if the response actions in the perched groundwater are not performing as expected, i.e., preventing contaminants from migrating to the Ogallala Aquifer. Currently, Pantex has begun investigation of detections of high explosives above groundwater protection standards in wells on the Texas Tech University property and a plume that is moving to the northeast from that area. Due to those detections, this Plan recognizes the fact that future detections in the Ogallala will be focused on first‐ time detections of analytes. After a remedy is determined, this Plan will require modification to address-deviations and contingent actions. This Plan was developed to identify the contingent actions necessary to mitigate impacts resulting from deviations to site conditions or response action performance. The Plan defines the environmental problem being addressed by the response actions, clarifies the expected conditions and objectives of the response actions, and identifies the potential deviations to the response actions (due to site conditions or technology performance) that could be encountered. The deviations were evaluated to determine the likelihood of occurrence, potential impact, and time to respond to avoid impact. The Plan also identifies the monitoring outlined in the Long‐Term Monitoring System Design Report (Consolidated Nuclear Security, 2024) and Sampling Analysis Plan (PanTeXas Deterrence, 2024) that will be used to detect the deviations. Lastly, the Plan specifies the contingent actions that could be implemented in response to the deviations. Because each response focuses on a discrete portion of the perched aquifer and contaminant plume, each response action has a different set of expected conditions, and therefore differing impacts from deviations to the site and technology expectations. As a result, the contingent actions are identified for each response action and potential deviation including specific constituents, location, and conditions. If deviations are encountered that impact the ability of the response action to meet performance objectives, the contingent actions will be focused on ensuring the response action can meet the performance objective. Contingent actions may be implemented as interim actions (ISMs/removal actions) in accordance with the Record of Decision, Interagency Agreement, and Hazardous Waste Permit‐50284, if warranted by the specific circumstances. For deviations to site characterization expected conditions, the contingent action will focus on determination of the source of the deviation, determination of the appropriate response, and evaluation of additional work to be completed. However, if the deviation to characterization impacts the performance of the response action, the contingent action will again focus on ensuring performance objectives can be met. Early source term removals and cleanup actions have been implemented to protect the Ogallala Aquifer. Because of these actions and based on modeling results, the expected conditions in the Ogallala Aquifer are that constituents of concern will not be detected above the Groundwater Protection Standards (GWPSs) nor will they reach potential points of exposure above the GWPS. The primary deviation of concern for the Ogallala is if constituents are detected in the Ogallala Aquifer near or above GWPSs. If it occurs, this change in expected conditions would require further evaluation of site and contaminant characteristics to determine an appropriate course of action. The evaluation would include additional monitoring, source identification, implementation of interim protective measures (if necessary), and delineation of extent. These evaluations are necessary to determine an appropriate response action for the Ogallala. The primary goal of the Plan is to provide for the continued protection of the Ogallala Aquifer and the health of its consumers. In recognition, this Plan presents a flexible and rational approach for making future decisions associated with confirming the change in perched and Ogallala aquifer conditions and identifying a response (technical activities, changes to response actions, regulatory oversight, and public involvement).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Los Alamos National Laboratory Floodplain Assessment for the High Explosive Transfer Facility Blast Radius Fence Project at Technical Area 08

The National Nuclear Security Administration (NNSA), a semi-autonomous agency within the U.S. Department of Energy (DOE), is proposing to take action at Los Alamos National Laboratory (LANL) within the upper Pajarito Canyon 100-year (yr.) floodplain at Technical Area (TA) 08. This action consists of installation of a 3-strand smooth wire fence to create an operational boundary in a radius of 1250 feet (ft.) from the TA-08 High Explosive (HE) Transfer Facility (Figure 1). The proposed fence is intended to provide a barrier between the public and direct access to the HE Transfer Facility.

54 ENVIRONMENTAL SCIENCES↗

Discrimination of Seismic Events (2006–2020) in North Korea Using P/Lg Amplitude Ratios from Regional Stations and a Bivariate Discriminant Function

Two events of magnitude (m b ) 3.6–3.8 occurred in southern North Korea (NK) on 27 June 2019 and 11 May 2020. Although these events were located ~330–400 km from the known nuclear test site, the fact that they occurred within the territory of NK, a country with a recent history of underground nuclear tests, made them events of interest for the monitoring community. In this work, we used P/Lg ratios from regional stations to categorize seismic events that occurred in NK from 2006 to May 2020, including these two recent events, the six declared NK nuclear tests, and the cavity collapse and triggered earthquakes that followed the 3 September 2017 nuclear explosion. We were able to separate the cavity collapse from the population of nuclear explosions. However, based on P/Lg ratios, the distinction between the earthquakes and the cavity collapse is ambiguous. The performed discriminant analyses suggest that combining Pg/Lg and Pn/Lg ratios results in improved discriminant power compared with any of the ratio types alone. We used the two ratio types jointly in a quadratic discriminant function and successfully classified the six declared nuclear tests and the triggered earthquakes that followed the September 2017 explosion. Our analyses also confirm that the recent southern events of June 2019 and May 2020 are both tectonic earthquakes that occurred naturally.

58 GEOSCIENCES↗

Water Cycle Science Enabled By the GRACE and GRACE-FO Satellite Missions

Satellite observations of the time-variable gravity field revolutionized the monitoring of large-scale water storage changes, beginning with the 2002 launch of the Gravity Recovery and Climate Experiment (GRACE) mission. Most hydrologists were sceptical of the satellite gravimetry approach at first, but validation studies assuaged their concerns and high-profile, GRACE-based groundwater depletion studies caused an explosion of interest. The importance of GRACE observations for hydrologic and cryospheric science became so great that GRACE Follow-On (GRACE-FO) jumped the National Aeronautics and Space Administration’s Earth science mission queue and launched in 2018. A third mass change mission is currently under development. Here, we review key milestones in satellite gravimetry’s progression from the fringes of hydrology to being a staple of large-scale water cycle and water resources studies and the sole source of observations of what is now an ‘essential climate variable’, terrestrial water storage.

climate change↗

A Multi-Analysis Approach for Estimating Regional Health Impacts from the 2017 Northern California Wildfires

In the evening of October 8 and early hours of October 9, 2017, high winds in Northern California downed trees and power lines, igniting some of the most devastating wildfires the state had seen, and within hours unhealthy air quality impacted millions of people. We simulated these air quality conditions using fire detection information from the MODIS, VIIRS, and GOES-16 ABI satellite instruments, and applying a set of three WRF–CMAQ simulations, one data fusion, and three machine learning methods. We investigated using the 5-min available GOES-16 fire detection data to simulate timing of fire activity to allocate emissions hourly for the WRF-CMAQ air quality modeling system. Interestingly, this approach did not necessarily improve results compared to the baseline case, which used a default time profile. However, this approach was key to simulating the initial 12-hr explosive fire activity and smoke impacts. The WRF-CMAQ simulations compared well with observational data for the October 8-15 time period and tended to overestimate concentrations October 16-20. To improve these results, we applied three machine learning algorithms. We also had a unique opportunity to evaluate results with temporary monitors deployed specifically for wildfires, and performance was markedly different. For example, at the permanent monitoring locations, the WRF-CMAQ simulations had a Pearson correlation of 0.65, and the data fusion approach improved this (Pearson correlation = 0.95), while at the temporary monitor locations across the WRF-CMAQ, data fusion, and machine learning datasets, the best Pearson correlation was 0.5. The data fusion and machine learning results were biased low and WRF-CMAQ results were biased high. Finally, we applied the optimized PM2.5 exposure estimate in a short-term exposure-response function. Total estimated mortality attributable to PM2.5 exposure during the smoke episode was 83 (95% confidence interval: 0, 196) with 47% of these deaths attributable to wildland fire smoke.

O'Neill, Susan↗

Panchromatic Observations of SN2011dh Point to a Compact Progenitor Star

We report the discovery and detailed monitoring of X-ray emission associated with the Type IIb SN2011dh using data from the Swift and Chandra satellites, placing it among the best studied X-ray supernovae to date. We further present millimeter and radio data obtained with the SMA, CARMA, and EVLA during the first three weeks after explosion. Combining these observations with early optical photometry, we show that the panchromatic dataset is well-described by non-thermal synchrotron emission (radio/mm) with inverse Compton scattering (X-ray) of a thermal population of optical photons. We derive the properties of the shockwave and the circumstellar environment and find a time-averaged shock velocity of v approximately equals 0.1c and a progenitor mass loss rate of M-dot approximately equals 6 X 10 (exp 5) Solar M/ yr (wind velocity, v(sub w) = 1000 km/s). We show that these properties are consistent with the sub-class of Type IIb supernovae characterized by compact progenitors (Type cIIb) and dissimilar from those with extended progenitors (Type eIIb). Furthermore, we consider the early optical emission in the context of a cooling envelope model to estimate a progenitor radius of R(sub star) approximately equals 10(exp 11) cm, in line with the expectations for a Type cIIb supernova. Together, these diagnostics suggest that the putative yellow supergiant progenitor star identified in archival HST observations is instead a binary companion or unrelated to the supernova. Finally, we searched for the high energy shock breakout pulse using X-ray and gamma-ray observations obtained during the purported explosion date range. Based on the compact radius of the progenitor, we estimate that the shock breakout pulse was detectable with current instruments but likely missed due to their limited temporal/ spatial coverage. Future all-sky missions will regularly detect shock breakout emission from compact SN progenitors enabling prompt follow-up observations of the shockwave with the EVLA and ALMA.

Soderberg, A. M.↗

Seismic Event Characterization Using Full Moment Tensors on the Hypersphere

Moment tensor solutions provide insights into the deformation that has occurred in the source region of a seismic event and are therefore of great value in identifying different types of seismic sources, such as when monitoring for underground nuclear tests. Despite this utility, inversion of waveforms recorded by seismometers for their full seismic moment tensor is not yet routine, and development of robust methods to classify events based on this information is in its infancy. Here, we assemble an inventory of 1405 full moment tensor solutions that include explosive, earthquake, and collapse events, and investigate the use of anisotropic probability distribution functions on the 5D hypersphere to discriminate between these sources. Using a Bayesian classifier, we obtain optimal success rates of 98.4% across all events and demonstrate that modification of the prior probabilities provides a natural way to alter the balance between not missing desirable events (such as explosions) versus misclassifying large numbers of undesired events (such as earthquakes). The approach is specifically designed to progress from traditional, bipolar event screening metrics to more generalized event identification across multiple types of seismic sources. Despite current databases containing insufficient numbers of events to definitively demonstrate at present, we also find intriguing evidence of subgroupings within individual source populations on the hypersphere, for example, between chemical and nuclear explosions, raising the potential possibility of discriminating between these event types in the future.

Geosciences↗

Ionospheric Signature of Surface Mine Blasts from Global Positioning System Measurements

Sources such as atmospheric or buried explosions and shallow earthquakes are known to produce infrasonic pressure waves in the atmosphere. Because of the coupling between neutral particles and electrons at ionospheric altitudes, these acoustic and gravity waves induce variations of the ionospheric electron density. The Global Positioning System (GPS) provides a way of directly measuring the total electron content in the ionosphere and, therefore, of detecting such perturbations in the upper atmosphere. In July and August 1996, three large surface mine blasts (1.5 Kt each) were detonated at the Black Thunder coal mine in eastern Wyoming. As part of a seismic and acoustic monitoring- experiment, we deployed five dual-frequency GPS receivers at distances ranging from 50 to 200 km from the mine and were able to detect the ionospheric perturbation caused by the blasts. The perturbation starts 10 to 15 min after the blast, lasts for about 30 min, and propagates with an apparent horizontal velocity of 1200 meters per second. Its amplitude reaches 3 x 10 (exp 14) el per square meters in the 7-3 min period band, a value close to the ionospheric perturbation caused by the M = 6.7 Northridge earthquake. The small signal-to-noise ratio of the perturbation can be improved by slant-stacking the electron content time-series recorded by the different GPS receivers taking into account the horizontal propagation of the perturbation. The energy of the perturbation is concentrated in the 200 to 300 second period band, a result consistent with previous observations and numerical model predictions. The 300 second band probably corresponds to gravity modes and shorter periods to acoustic modes, respectively. Using a 1-D stratified velocity model of the atmosphere we show that linear acoustic ray tracing fits arrival times at all GPS receivers. We interpret the perturbation as a direct acoustic wave caused by the explosion itself. This study shows that even relatively small subsurface events can produce ionospheric perturbations that are above the detection threshold of the GPS technique. By sensing derivative signals, which can be detected over a relatively broad region, it appears that GPS might be particularly useful for source characterization within the first acoustic quiet zone where infrasound would probably be ineffective. This suggests that dual-frequency GPS monitoring could contribute to Comprehensive Nuclear Test Ban Treaty verification.

Calais, Eric↗

Crush Test Abuse Stand

The purpose of this system is to simulate an internal short on battery cells by causing deformation (a crushing force) in a cell without penetration. This is performed by activating a hydraulic cylinder on one side of a blast wall with a hydraulic pump located on the other. The operator can control the rate of the crush by monitoring a local pressure gauge connected to the hydraulic cylinder or a load cell digital display located at the hydraulic pump control area. The internal short simulated would be considered a worst-case scenario of a manufacturer fs defect. This is a catastrophic failure of a cell and could be a very destructive event. Fully charged cells are to have an internal short simulated at the center of the length of the cell (away from terminals). The crush can be performed with a .- to 1-in. (.0.6- to 2.5-cm) rod placed crossways to the cell axis, causing deformation of the cell without penetration. The OCV (open-circuit voltage) and temperature of the cells, as well as the pressure and crushing force, are recorded during the operation. Occurrence of an internal short accompanied by any visible physical changes such as venting, fires, or explosions is reported. Typical analytical data examined after the test would be plots of voltage, temperature, and pressure or force versus time. The rate of crushing force can be increased or decreased based on how fast the operator pumps the hydraulic pump. The size of cylinder used to compress the battery cell can be easily changed by adding larger or smaller fittings onto the end of the hydraulic cylinder based on the battery/cell size being tested. The cell is crushed remotely and videotaped, allowing the operator to closely monitor the situation from a safe distance.

Collins, Jacob↗

A real-time search for Type Ia Supernovae with late-time interactions with circumstellar material in ZTF data

While it is generally accepted that Type Ia supernovae (SNe Ia) are the terminal explosions of white dwarfs (WDs), the nature of their progenitor systems and the mechanisms that lead up to these explosions remain widely debated. In rare cases, the SN ejecta interact with circumstellar material (CSM) that had previously been ejected from the progenitor system. The longer the delay between the creation of the CSM and the SN explosion, the greater the distance between the SN explosion site and the CSM and the later the onset of the interaction. The unknown distance between the CSM and SN explosion site makes it impossible to predict when the interaction will start. If the time between the SN explosion and the onset of the CSM interaction is of the order of several months to years, the SN has generally faded and it is no longer actively followed up on. This makes it even more difficult to detect the interaction while it is happening. In this work, we report on a real-time monitoring programme running between 13 November 2023 and 9 July 2024. It monitored 6914 SNe Ia for signs of late-time rebrightening using the Zwicky Transient Facility (ZTF). Flagged candidates were rapidly followed up on with photometry and spectroscopy to confirm the late-time excess and its position. We report the discovery of a ∼50 day rebrightening event in SN 2020qxz around 1200 rest-frame days after the peak of its light curve. SN 2020qxz exhibited signs of an early CSM interaction, but had faded from view over two years before its reappearance. Initial follow-up spectroscopy revealed the presence of four emission lines, while later follow-up spectroscopy showed that these had faded shortly after the end of the ZTF-detected rebrightening event. Our best match for these emission lines are H β (blueshifted by ∼5900 km s −1 ) and Ca II λ8542 , N I λ8567 , and K I λλ8763, 8767 (all blueshifted by 5100 km s −1 ; although we note that the line identifications are uncertain). This shows that catching and following up on late-time interactions as they occur can offer new clues on the nature of the progenitor systems that produce these SNe by putting constraints on the possible type of donor star. The only way to do this systematically is to use large sky surveys such as ZTF and the upcoming Vera C. Rubin Observatory’s Legacy Survey of Space and Time (LSST) to monitor a large sample of objects for the rare events that reappear long after the object has faded from view.

circumstellar matter↗

Real-time applications for remotely sensed lightning observations

Real-time monitoring of lightning activity, including the estimation of the location where the flash strikes the earth, provides a basis for the assessment of the impact of both the lightning hazard and severe weather on daily test and operation schedules. Applications to range meteorology are considered, taking into account the Space Transportation System, the handling of explosive materials, and test operations. Attention is given to details regarding forecasting and nowcasting, a low-cost lightning display system for electrical storm tracking and the support of lightning research activities, two versions of employed software, an example of the storm diagnostics and display products which will be available from the next generation weather radar, and a summary of the evolution of 12 New Mexico thunderstorms.

Goodman, S. J.↗

Subsurface spectroscopy of heterogeneous materials using optical wavefront shaping

Plastic-bonded explosives and propellants consist of energetic molecular crystals embedded in a polymer matrix that also includes additives, such as taggants, plasticizers, grit, antioxidants, etc. Furthermore, this heterogeneous composition renders these materials optically opaque, limiting optical characterization techniques (e.g., Raman spectroscopy) to probing chemical reactions at the sur-face. However, many reactions of interest are believed to occur inside the material where current optical techniques can not probe. To address this challenge, we have developed a new optical technique that utilizes feedback assisted wavefront shaping (using spatial light modulators) to focus probe light inside a heterogeneous material, such that we can perform subsurface spectroscopy. Recently, we demonstrated this technique by monitoring subsurface photodegradation and thermal degradation of Eu-doped molecular crystals using both fluorescence and Raman spectroscopy. Based off these successes we are now looking into improvements to this technique to allow for greater time resolution, including replacing our spatial light modulators with fast digital mirror devices, using ns optical phase conjugation, and self-healing optical beams.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

On various mechanisms of shear wave generation from underground chemical explosions in hard rocks

Significant shear waves observed during underground explosions can make them look similar to natural earthquakes. In the near field, they can be caused by the source geometry and emplacement conditions, heterogeneous and anisotropic properties of the rock mass, relaxation of a non-uniform in situ stress around the source after the departure of the shock wave as well as by activation of nearby faults. Here, the stress relaxation mechanism is important at depth when the difference in main principal stresses becomes significant. Another mechanism of shear motion generation is due to explosive gas flow from the source cavity into the cracks initiated by the shock wave and joints connected to the cavity. Also, the shock wave from the explosion can cause weakening of the rock either due to permanent damage or shock-induced liquefaction at saturated joints. Numerical modelling of underground chemical explosions in hard rock at various confinement stresses has been conducted to compare various mechanisms of shear wave generation and their seismic characteristics with that of a simple fault rupture process observed during an earthquake.

58 GEOSCIENCES↗

Differential Muon Tomography to Continuously Monitor Changes in the Composition of Subsurface Fluids

Muon tomography has been used to seek hidden chambers in Egyptian pyramids and image subsurface features in volcanoes. It seemed likely that it could be used to image injected, supercritical carbon dioxide as it is emplaced in porous geological structures being used for carbon sequestration, and also to check on subsequent leakage. It should work equally well in any other application where there are two fluids of different densities, such as water and oil, or carbon dioxide and heavy oil in oil reservoirs. Continuous monitoring of movement of oil and/or flood fluid during enhanced oil recovery activities for managing injection is important for economic reasons. Checking on leakage for geological carbon storage is essential both for safety and for economic purposes. Current technology (for example, repeat 3D seismic surveys) is expensive and episodic. Muons are generated by high- energy cosmic rays resulting from supernova explosions, and interact with gas molecules in the atmosphere. This innovation has produced a theoretical model of muon attenuation in the thickness of rock above and within a typical sandstone reservoir at a depth of between 1.00 and 1.25 km. Because this first simulation was focused on carbon sequestration, the innovators chose depths sufficient for the pressure there to ensure that the carbon dioxide would be supercritical. This innovation demonstrates for the first time the feasibility of using the natural cosmic-ray muon flux to generate continuous tomographic images of carbon dioxide in a storage site. The muon flux is attenuated to an extent dependent on, amongst other things, the density of the materials through which it passes. The density of supercritical carbon dioxide is only three quarters that of the brine in the reservoir that it displaces. The first realistic simulations indicate that changes as small as 0.4% in the storage site bulk density could be detected (equivalent to 7% of the porosity, in this specific case). The initial muon flux is effectively constant at the surface of the Earth. Sensitivity of the method would be decreased with increasing depth. However, sensitivity can be improved by emplacing a greater array of particle detectors at the base of the reservoir.

Coleman, Max↗