SEARCH · Search NASA

Results for “Single event effects”

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.

SourceResearch connectionFrom yearThrough year

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

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2024

A review of criticality dosimetry at the Y-12 National Security Complex and practical importance of dose accuracy in emergency response

A nuclear criticality results in the emission of both neutron and gamma radiation and can produce doses to personnel near the event that exceed 0.1 Gy (10 rad). The primary purpose of nuclear accident dosimetry is to rapidly identify affected personnel in need of prompt medical treatment and to reassure personnel who have been only minimally exposed. While accurate dosimetry is desired, it must be recognized that dose determinations made from whole-body dosimeters or simple triage methods are very rough estimates and contain significant uncertainties. Even when accounting for factors like varying neutron energy spectra, mean photon energies, body orientation within the radiation field, and transient effects on dosimeter response, etc., the end value is a dosimetric quantity defined for very specific radiological conditions and determined within a simple phantom usually at a single depth. Of more importance is the biological response to the radiation, which will vary by person and can be affected by the individual’s radiation sensitivity, age, gender, mass, and underlying health conditions. The overall biological, person-specific response to a given dose cannot be precisely determined except by patient symptom observation and individual biological dosimetry (e.g. chromosome analysis, lymphocyte ratios, etc.). This work describes and discusses the criticality accident dosimetry program at the Y-12 National Security Complex, a United States Department of Energy National Nuclear Security Administration facility. In conclusion, the primary goals of the Y-12 accident dosimetry program are, among others, the rapid identification of significantly exposed persons, prompt routing of exposed workers for medical evaluation and treatment, and the ultimate processing of dosimeters to assign doses to personnel.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2025

Search for Higgs boson decays into a Z boson and a light hadronically decaying resonance in pp collisions at s =13 TeV with the ATLAS detector

A search for decays of the Higgs boson into a Z boson and a light resonance, with a mass of 0.5–3.5 GeV, is performed using the full 140 fb−1 dataset of 13 TeV proton–proton collisions recorded by the ATLAS detector during LHC Run 2. Leptonic decays of the Z boson and hadronic decays of the light resonance are considered. The resonance can be interpreted as a J / ψ or η c meson, an axion-like particle, or a light pseudoscalar predicted in two-Higgs-doublet models. Due to its low mass, this resonance is produced with a high Lorentz boost in the laboratory frame and therefore reconstructed as a single small-radius jet of hadrons. A neural network is used to correct the Monte Carlo simulation of the total expected background using data from sideband regions. Two additional neural networks are used to distinguish signal from background, enhancing the purity of the signal region. A binned profile-likelihood fit is performed on the final-state invariant mass distribution. No significant excess of events relative to the expected background is observed, and upper limits at 95% confidence level are set on the Higgs boson's branching fraction to a Z boson and a light resonance. The exclusion limit is ∼10% for the lower masses, and increases for higher masses. Upper limits on the effective coupling C Z H eff / Λ of an axion-like particle to a Higgs boson and Z boson are also set at 95% confidence level, and range from 0.9 to 2 TeV−1.

Aad, G↗

DOE OSTI2021

Pulse-shape discrimination against low-energy Ar-39 beta decays in liquid argon with 4.5 tonne-years of DEAP-3600 data

The DEAP-3600 detector searches for the scintillation signal from dark matter particles scattering on a 3.3 tonne liquid argon target. The largest background comes from 39 Ar beta decays and is suppressed using pulse-shape discrimination (PSD). We use two types of PSD estimator: the prompt-fraction, which considers the fraction of the scintillation signal in a narrow and a wide time window around the event peak, and the log-likelihood-ratio, which compares the observed photon arrival times to a signal and a background model. We furthermore use two algorithms to determine the number of photons detected at a given time: (1) simply dividing the charge of each PMT pulse by the mean single-photoelectron charge, and (2) a likelihood analysis that considers the probability to detect a certain number of photons at a given time, based on a model for the scintillation pulse shape and for afterpulsing in the light detectors. The prompt-fraction performs approximately as well as the log-likelihood-ratio PSD algorithm if the photon detection times are not biased by detector effects. We explain this result using a model for the information carried by scintillation photons as a function of the time when they are detected.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

The Promise of Mutation Resistant Drugs for SARS-CoV-2 That Interdict in the Folding of the Spike Protein Receptor Binding Domain

In recent work, we proposed that effective therapeutic drugs aimed at treating the SARS-CoV-2 infection could be developed based on interdicting in the early steps of the folding pathway of key viral proteins, including the receptor binding domain (RBD) of the spike protein. In order to provide for a drug target on the protein, the earliest contact-formation event along the dominant folding pathway of the RBD spike protein was predicted employing the Sequential Collapse Model (SCM). The segments involved in the predicted earliest contact were suggested to provide optimal folding interdiction target regions (FITRs) for potential therapeutic drugs, with a focus on folding interdicting peptides (FIPs). In this paper, we extend our analysis to include 13 known single mutations of the RBD spike protein as well as the triple mutation B1.351 and the recent double mutation B1.617.2. The results show that the location of the FITR does not change in any of the 15 studied mutations, providing for a mutation-resistant drug design strategy for the RBD-spike protein.

Bergasa-Caceres, Fernando↗

DOE OSTI2024

Seismic Signal Detection on International Monitoring System 3-Component Stations using PhaseNet

In this report we discuss training a deep learning seismic signal detection model on 3-component stations from the International Monitoring System (IMS) using the PhaseNet architecture. Using 14 years of associated signals from the International Data Centre’s (IDC) Late Event Bulletin (LEB), we auto-curated training data consisting of signal windows containing associated arrivals, and noise windows that contain no LEB-associated signals. We trained several models using different waveform window durations (30 seconds and 100 seconds), with and without bandpass filtering. We evaluated the effectiveness of our models using associated signals from the Unconstrained Global Event Bulletin (UGEB) and found that several of our models outperformed the signal detections from the IDC’s Selected Event List 3 (SEL3) arrival table. The SEL3 bulletin evaluated on the UGEB dataset with 100-second waveform windows registered a precision and recall of .15 and .48, respectively, versus .19 and .59 for our filtered-data model. For the 30-second waveform window dataset, the SEL3 bulletin achieved a precision and recall of .31 and .47, respectively, versus .32 and .60 for our filtered-data model. Finally, our models detected signals from all source-to-receiver distances, suggesting it is feasible to use a single PhaseNet model for the IMS network.

58 GEOSCIENCES↗

DOE OSTI2024

FY24 Report on Water NSTF Testing: Lower Tank Inlet Piping Configuration

The following report serves as a summary of the accomplishments and testing results by the Natural convection Shutdown heat removal Test Facility (NSTF) experimental program over the past 12-month period. A major activity included reconfiguration of the chimney piping geometry which altered the discharge position into the tank from the original 50% tank height to a new lower position at a 10.9% tank height. This modification increases the volume of available coolant thus extending long-term operating capacity, however also results in a decreased liquid driving head which has the potential to reduce the natural circulation efficiency and decrease overall performance. Examination of the tradeoffs for these two configurations is an area of interest for RCCS designers and drove planned test activities this year. Ten matrix tests were performed in FY24, totaling 195 hours of heated operations and 12.9 MWh of electrical heating throughout the year, with eight classified as Accepted per NQA-1 and two classified as Trending. Testing prior to the facility reconfiguration examined the facility response to throttling at the tank inlet, demonstrating increased sensitivity of flow instabilities to throttling within the two-phase region when compared to throttling at the lower sensitivity single-phase inlet region. Testing after the chimney reconfiguration began with a baseline test at conditions of 80% inventory fill and prototypic thermal power input of 2.1 MW t . In addition to establishing a reference for nominal system behavior and performance, repeat testing with multiple subsequent runs demonstrated strong repeatability between tests performed at the same conditions in this new configuration. These tests also examined if influences would occur to system behavior with installation of new higher-resolution instrumentation within the upper chimney. In this critical region where boiling and flashing phenomena dominate, the test results provided confidence that the new instrumentation does not uncharacteristically influence the observed behaviors. An inventory parametric series was then initiated to examine system behavior with the new chimney configuration at six varying initial inventory levels, ranging from high 80% to low 20% fill. Two-phase oscillations with similar peak and mean flow rates were observed when comparing to the previous mid tank configuration. Generally, similar system response trends with inventory were also observed, such as two phase oscillations suppressing as inventories were lowered. However, in one absolute inventory comparison at the highest fill of 80%, flow oscillations saw a gradual growth over the 4 hours of two-phase operation in the new lower tank inlet configuration, opposite of the gradual dampening observed in the mid tank inlet configuration. This can be attributed to, in part, a greater hydrostatic head pressure above the two-phase discharge region where the boiling front is developed. Furthermore, the change resulted in greater sensitivity to liquid degassing phenomena during single-phase heating, causing loop instabilities to form which trigged moderate flow degradation during the period approaching saturation and boiling conditions. This behavior had been observed previously under some conditions but was a common occurrence in recent testing with the newer lower tank inlet configuration. Initial observations from these first data sets suggest an overall larger window of stability for the mid tank inlet configuration compared to the lower tank inlet. Lastly, the lowest inventory fill test was repeated over an extended testing window to examine depletion behavior. Natural circulation flow and effective heat removal performance were observed during most of the testing period; only after the tank became fully drained (0% fill) did flow stagnate and violent geysering events occur. This early observation confirms one relative advantage over the mid tank inlet configuration, which stagnated under comparable conditions at ~20% inventory remaining in the tank.

42 ENGINEERING↗

DOE OSTI2023

Simulations of Water Flow in Relation to a Steady Criticality in an Unsaturated Alluvial Repository

Oak Ridge National Laboratory (ORNL) is examining the conditions under which a criticality event could be ruled out for large packages placed in an unsaturated alluvial repository. This work specifically examines the water level in a stylized waste package that is assumed (1) to fail on the top, and (2) to empty as water is driven away—not by boiling—but by fluid flow phenomena driven by thermal gradients. This analysis uses the PFLOTRAN code package. Specifically, this year’s work includes examination of the heat rate at which water is driven away for various percolation rates using the same assumptions listed in the 2022 annual report. The STANDARDS (formerly UNF-ST&DARDS) package has predicted dual-purpose canister (DPC) neutron multiplication factors (k-effectives) ranging above 1.1 (unitless) with failed components. The DPC designs were not created with the intent to use for disposal purposes. Emplacing such packages may be worth consideration if it can be shown that criticality is either extremely unlikely or of very low consequence to the repository’s safety case. This work aims to help establish hydrologic boundary conditions for a highly hypothetical, highly stylized subsurface criticality event. The approach to critical is assumed to be slow, producing steady heating rates under a few kilowatts. For a DPC containing capillary media, the simulated power required for dry-out was proportional to the percolation rate but was unaffected by the size of a large breech on top of the DPC. Results for single cases with no capillary media in the DPC suggest that higher infiltration rates may be required in those cases to cause flooding. Bounding mechanisms and assumptions must be identified and characterized further by performing additional perturbations to various input parameters to demonstrate robustness of the results. This might be accomplished by examining responses to power pulses to determine the system transfer function required to enable a reduced order feedback model.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

Alpha-decay induced shortening of fission tracks simulated by in situ ion irradiation

The diffusion of defects (e.g., vacancies and interstitials) and elements used for dating (e.g., He and Pb) in a mineral structure is a thermal process: This is the primary assumption used to determine the age and thermal history of minerals. For instance, thermal history reconstruction, through the number and length distribution of tracks produced by spontaneous fission of 238 U, is obtained by assuming a thermal event to be the only energy source for shortening of fission tracks. Here, we report a new, non-thermal energy source that induces additional shortening of fission tracks by the irradiation of alpha-recoils from the alpha-decay of 238 U and 232 Th. We simulate alpha-decay induced track-shortening by combining ion accelerator irradiations with transmission electron microscopy. This allows for the first observation of track-shrinkage during in situ ion irradiation. We show that rather than alpha-particles, alpha-recoils induce a significant shortening of fission tracks by nuclear-collisions. The shortening of track-length can be quantified as a function of alpha-decay event dose. However, apatite is less sensitive than zircon to this non-thermal process. Furthermore, the findings exemplify the interactions among different types of self-irradiation from alpha-particles, alpha-recoils and fission-fragment nuclei in single mineral grains and have important implications for the use of zircon and apatite for radiometric dating and thermochronology.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Search for new physics in the final state with a single photon and large missing transverse momentum in proton-proton collisions at √𝑠 =13 TeV

A search for new physics in events featuring a single photon and missing transverse momentum is presented, using proton-proton √𝑠 =13 TeV collision data corresponding to an integrated luminosity of 101 fb −1 collected by the CMS experiment at the CERN LHC between 2017 and 2018. This analysis, combined with a previous study of 36 fb −1 of 2016 data (totaling 137 fb −1 ), reveals no significant deviations from standard model expectations. The results are then used to establish 95% confidence level limits on parameters in theoretical models involving dark matter and large extra dimensions. Compared to the 2016-only analysis, this search achieves up to a 14% improvement in exclusion reach for mediator masses in simplified dark matter models, along with 11% and 10% enhancements in the limits on the effective field theory suppression scale and the fundamental Planck scale, respectively. These results are the most stringent constraints on these parameters to date.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Measurement of $\nu_\mu$ CC Interactions With Two-Proton Final State in MINERvA

This dissertation presents a measurement of charged–current (CC) muon–neutrino interactions with exactly two protons and no pions in the final state (CC~$2p\,0\pi$), using data collected by the MINERvA detector in the NuMI medium–energy beam at Fermilab. Such two–proton topologies are a sensitive probe of nuclear dynamics in the few–GeV regime, including multi–nucleon correlations (npnh, notably $2p2h$) and intranuclear final–state interactions (FSI) such as pion absorption and nucleon rescattering. A precise experimental characterization of these processes is essential both for neutrino–interaction theory and for reducing systematic uncertainties in oscillation experiments that rely on accurate modeling of neutrino–nucleus interactions. Events are selected by requiring a $\nu_\mu$ CC interaction with a reconstructed $\mu^-$ and two proton tracks originating from a common vertex in MINERvA’s finely segmented scintillator tracker, with no reconstructed mesons. Muon charge and momentum are constrained by matching to the MINOS Near Detector, while proton identification exploits energy–loss profiles and stopping–proton features. Backgrounds from pion–producing channels that enter the signal region through FSI or reconstruction effects are constrained with data–driven sidebands (Michel–electron and isolated–cluster “blob” samples) and tuned via a simultaneous fit across signal and sideband regions. To correct detector resolution and acceptance effects, the analysis employs iterative Bayesian unfolding with extensive validation: statistical pseudo–experiments, and robustness checks against generator systematic “universes” and additional strong shape warps. Single–differential cross sections are reported for three observables tailored to the two–proton final state: the opening–angle cosine $\cos\!\left(\theta_{pp}\right)$, the leading–proton momentum, and the subleading–proton momentum. Systematic uncertainties include contributions from neutrino flux, interaction modeling (e.g., npnh and resonance parameters, pion FSI), and detector response (calibration, reconstruction efficiencies). The resulting distributions provide targeted constraints on the interplay of multi–nucleon dynamics and FSI that shape CC~$2p\,0\pi$ final states on hydrocarbon. Comparisons to modern GENIE–based simulations highlight kinematic regions where model components require refinement. These measurements thus inform generator tuning and improve the reliability of neutrino–energy reconstruction strategies for current and future long–baseline oscillation programs.

Syrotenko, Vladyslav S. [Tufts U.]↗

DOE OSTI2024

Measurement of the 1-jettiness event shape observable in deep-inelastic electron-proton scattering at HERA

The H1 Collaboration reports the first measurement of the 1-jettiness event shape observable $τ^b_1$ in neutral-current deep-inelastic electron-proton scattering (DIS). The observable $τ^b_1$ is equivalent to a thrust observable defined in the Breit frame. The data sample was collected at the HERA ep collider in the years 2003–2007 with center-of-mass energy of $\sqrt{s}$ = 319 GeV, corresponding to an integrated luminosity of 351.1 pb -1 . Triple differential cross sections are provided as a function of $τ^b_1$, event virtuality $Q^2$, and inelasticity y, in the kinematic region $Q^2$ > 150 GeV 2 . Single differential cross section are provided as a function of $τ^b_1$ in a limited kinematic range. Double differential cross sections are measured, in contrast, integrated over $τ^b_1$ and represent the inclusive neutral-current DIS cross section measured as a function of $Q^2$ and y. The data are compared to a variety of predictions and include classical and modern Monte Carlo event generators, predictions in fixed-order perturbative QCD where calculations up to $\mathcal{O}$($α^3_s$) are available for $τ^b_1$ or inclusive DIS, and resummed predictions at next-to-leading logarithmic accuracy matched to fixed order predictions at $\mathcal{O}$($α^2_s$). These comparisons reveal sensitivity of the 1-jettiness observable to QCD parton shower and resummation effects, as well as the modeling of hadronization and fragmentation. Within their range of validity, the fixed-order predictions provide a good description of the data. Monte Carlo event generators are predictive over the full measured range and hence their underlying models and parameters can be constrained by comparing to the presented data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

A Deep Exon Cryptic Splice Site Promotes Aberrant Intron Retention in a Von Willebrand Disease Patient

A translationally silent single nucleotide mutation in exon 44 (E44) of the von Willebrand factor (VWF) gene is associated with inefficient removal of intron 44 in a von Willebrand disease (VWD) patient. This intron retention (IR) event was previously attributed to reordered E44 secondary structure that sequesters the normal splice donor site. We propose an alternative mechanism: the mutation introduces a cryptic splice donor site that interferes with the function of the annotated site to favor IR. We evaluated both models using minigene splicing reporters engineered to vary in secondary structure and/or cryptic splice site content. Analysis of splicing efficiency in transfected K562 cells suggested that the mutation-generated cryptic splice site in E44 was sufficient to induce substantial IR. Mutations predicted to vary secondary structure at the annotated site also had modest effects on IR and shifted the balance of residual splicing between the cryptic site and annotated site, supporting competition among the sites. Further studies demonstrated that introduction of cryptic splice donor motifs at other positions in E44 did not promote IR, indicating that interference with the annotated site is context dependent. We conclude that mutant deep exon splice sites can interfere with proper splicing by inducing IR.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Unified 0.25-degree gridded infrastructure-critical extreme weather for the United States from 1979 to 2100

Extreme weather events can severely disrupt critical infrastructure, triggering cascading effects on power, transportation, and essential services. However, standard weather and climate datasets often lack specialized variables necessary for hazard assessments. We present a unified dataset of infrastructure-critical weather and climate variables across the United States at 0.25° resolution, covering daily or sub-daily intervals from 1979 to 2100. The dataset includes temperature, dew point, wind gusts, precipitation partitioned by rain, snow, and freezing rain or ice pellets, lightning, and wildfire metrics. Historical conditions (1979-2023) are synthesized from observations and reanalysis products, while future projections are derived from 14 CMIP6 global climate models (historical, SSP245, and SSP585 experiments). Physically based and data-driven methods are used to estimate variables not directly provided by existing models. By integrating these variables into a single unified dataset, we enable consistent, high-resolution assessments of weather-related infrastructure risks across past and future periods, supporting wide-ranging applications in energy, transportation, water resources, emergency management, and beyond.

Climate and Earth system modelling↗

DOE OSTI2021

Effect of Glacial/Interglacial Recharge Conditions on Flow of Meteoric Water Through Deep Orogenic Faults: Insights Into the Geothermal System at Grimsel Pass, Switzerland

Many meteoric-recharged, fault-hosted geothermal systems in amagmatic orogenic belts have been active through the Pleistocene glacial/interglacial climate fluctuations. The effects of such climate-induced recharge variations on fluid flow patterns and residence times of the thermal waters are complex and may influence how the geothermal and mineralization potential of the systems are evaluated. We report systematic thermal–hydraulic simulations designed to reveal the effects of recharge variations, using a model patterned on the orogenic geothermal system at Grimsel Pass in the Swiss Alps. Fault-bound circulation of meteoric water through the Grimsel fault is driven to depths of ~10 km by the high alpine topography. Simulations suggest that the current single pass flow is typical of interglacial periods, during which 1) meteoric recharge into the fault is high (above tens of centimeters per year), 2) conditions are at or somewhat below the critical Rayleigh number, and 3) the hydraulic connectivity along the fault plane is extensive (an extent of at least 10 km into increasingly higher terrain is required to explain the 10 km penetration depth). The subcritical condition constrains the bulk fault permeability to < 1e-14 m2. In contrast, the limited recharge during the numerous Pleistocene glaciation events likely induced a layered flow system with single pass flow confined to shallow depths while non-Rayleigh convection occurred in the deep fault. The same layering can be observed at low aspect ratios (length/depth) of the fault plane, when the available recharge area limits the flux through the fault.

58 GEOSCIENCES↗

DOE OSTI2022

NF-κB perturbation reveals unique immunomodulatory functions in Prx1 + fibroblasts that promote development of atopic dermatitis

Skin is composed of diverse cell populations that cooperatively maintain homeostasis. Up-regulation of the nuclear factor κB (NF-κB) pathway may lead to the development of chronic inflammatory disorders of the skin, but its role during the early events remains unclear. Here, through analysis of single-cell RNA sequencing data via iterative random forest leave one out prediction, an explainable artificial intelligence method, we identified an immunoregulatory role for a unique paired related homeobox-1 (Prx1) + fibroblast subpopulation. Disruption of Ikkb–NF-κB under homeostatic conditions in these fibroblasts paradoxically induced skin inflammation due to the overexpression of C-C motif chemokine ligand 11 (CCL11; or eotaxin-1) characterized by eosinophil infiltration and a subsequent T H 2 immune response. Because the inflammatory phenotype resembled that seen in human atopic dermatitis (AD), we examined human AD skin samples and found that human AD fibroblasts also overexpressed CCL11 and that perturbation of Ikkb–NF-κB in primary human dermal fibroblasts up-regulated CCL11. Monoclonal antibody treatment against CCL11 was effective in reducing the eosinophilia and T H 2 inflammation in a mouse model. Together, the murine model and human AD specimens point to dysregulated Prx1 + fibroblasts as a previously unrecognized etiologic factor that may contribute to the pathogenesis of AD and suggest that targeting CCL11 may be a way to treat AD-like skin lesions.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

Observation of double J / ψ meson production in p Pb collisions at s NN = 8.16 TeV

The first observation of the concurrent production of two J / ψ mesons in proton-nucleus collisions is presented. The analysis is based on a proton-lead ( p Pb ) data sample recorded at a nucleon-nucleon center-of-mass energy of 8.16 TeV by the CMS experiment at the CERN LHC and corresponding to an integrated luminosity of 174.6 nb − 1 . The two J / ψ mesons are reconstructed in their μ + μ − decay channels with transverse momenta p T > 6.5 GeV and rapidity | y | < 2.4 . Events where one of the J / ψ mesons is reconstructed in the dielectron channel are also considered in the search. The p Pb → J / ψ J / ψ + X process is observed with a significance of 5.3 standard deviations. The measured inclusive fiducial cross section, using the four-muon channel alone, is σ ( p Pb → J / ψ J / ψ + X ) = 22.0 ± 8.9 ( stat ) ± 1.5 ( syst ) nb . A fit of the data to the expected rapidity separation for pairs of J / ψ mesons produced in single (SPS) and double (DPS) parton scatterings yields σ SPS p Pb → J / ψ J / ψ + X = 16.5 ± 10.8 ( stat ) ± 0.1 ( syst ) nb and σ DPS p Pb → J / ψ J / ψ + X = 5.4 ± 6.2 ( stat ) ± 0.4 ( syst ) nb , respectively. This latter result can be transformed into a lower bound on the effective DPS cross section, closely related to the squared average interparton transverse separation in the collision, of σ eff > 1.0 mb at 95% confidence level. © 2024 CERN, for the CMS Collaboration 2024 CERN

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Identification of Faults Susceptible to Induced Seismicity (Final Report)

Central to the work documented in this report is the capability of geocellular models to represent the geologic conceptual model updated with fault identification from machine learning and joint inversion modeling of microseismic data measured and recorded as a consequence of CO 2 injection at a field demonstration site: the Illinois Basin - Decatur Project (IBDP). This work required seven unique geocellular models with 100s of simulated variations to gain a very high degree of confidence in the identification of geologic features present that contributed to induced microseismicity at IBDP. All forward modeling: pressure modeling, stress modeling, and seismic modeling used the same geologic conceptual model and representations of that model at different scales. The pressure modeling and poroelastic modeling created “snapshots” of pore pressure and stress field changes at different times during CO 2 injection, in which microseismic events were clustered (in time). These pressure and stress snapshots, within the framework and architecture of the geologic conceptual model via the geocellular model, informed the single fault and fault network models to ascertain the likelihood of fault movement (seismic or aseismic). The outcomes of the pressure, stress, and fault/fault network (seismic) modeling confirmed that the faults in the geologic conceptual model in Task 2 were likely the source of microseismic events measured at IBDP and acted as conduits for pressure to be transmitted from the injection interval into the Precambrian crystalline basement rock. This closely coordinated and integrated unique modeling approach was conducted to prove the viability of our proposed workflow 1) to better resolve crystalline basement faults, 2) detect subseismic faults that could be activated by injection, 3) increase the certainty in fault detection and their susceptibility to release seismic energy, and 4) understand transmission of pressure vertically from the well to the underlying fractured crystalline basement. The proposed methodology was effective in guiding an iterative process of calibrating forward modeling results based on similar geocellular models while honoring the geologic conceptual model (i.e., characterization data and knowledge of regional geology); this led to higher level of certainty in the identification of fault/faults zones to control seismicity and transmission of pressure to the regions of recorded and located injection induced seismicity.

58 GEOSCIENCES↗

DOE OSTI2023

Testing and Design of Discriminants for Local Seismic Events Recorded during the Redmond Salt Mine Monitoring Experiment

The Redmond Salt Mine (RSM) Monitoring Experiment in Utah was designed to record seismoacoustic data at distances less than 50 km for algorithm testing and development. During the experiment from October 2017 to July 2019, six broadband seismic stations were operating at a time, with three of them having fixed locations for the duration, whereas the three other stations were moved to different locations every one-and-half to two-and-half months. RSM operations consist of nighttime underground blasting several times per week. The RSM is located in proximity to a belt of active seismicity, allowing direct comparison of natural and anthropogenic sources. Using the recorded data set, we built 1373 events with local magnitude (M L ) of -2.4 and lower to 3.3. For 75 blasts (RMEs) from the Redmond Salt Mine and 206 tectonic earthquakes (EQs), both M L and the coda duration magnitude (MC) are well constrained. We used these events to test and design discriminants that separate the RMEs from the EQs and are effective at local distances. The discriminants consist of M L -M C , low-frequency Sg to high-frequency Sg, Pg/Sg phase-amplitude ratios, and Rg/Sg spectral amplitude ratios, as well as different combinations of two or more of these classifiers. The areas under the receiver operating characteristic curves (AUCs) of 0.92–1.0 for M L -M C , low-frequency Sg to high-frequency Sg, and Rg/Sg indicate that these discriminants are very effective. Conversely, the AUC of only 0.57 for Pg/Sg suggests that this discriminant is only slightly better than a random classifier. Among the effective classifiers, Rg/Sg, shows the lowest likelihood of misclassification (4.3%) for the populations. In conclusion, results of joint discriminant analyses suggest that even the arguably ineffective single classifier, like Pg/Sg in this case, can provide some value when used in combination with others.

58 GEOSCIENCES↗