Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model interpretability”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A search for dark matter produced in association with a dark Higgs boson decaying into a Higgs boson pair in 3b or 4b final states using pp collisions at $\sqrt{s}=13$ TeV with the ATLAS detector

A search is performed for dark matter particles produced in association with a resonant pair of Higgs bosons using 140 fb −1 of proton-proton collisions at a centre-of-mass energy of 13 TeV recorded by the ATLAS detector at the Large Hadron Collider. This signature is expected in some extensions of the Standard Model predicting the production of dark matter particles, and is interpreted in terms of a dark Higgs model containing a Z′ mediator in which the dark Higgs boson s decays into a pair of Higgs bosons. The dark Higgs boson is reconstructed through final states with at least three b-tagged jets, produced by the pair of Higgs boson decays, in events with significant missing transverse momentum consistent with the presence of dark matter. The observed data are found to be in good agreement with Standard Model predictions, constraining scenarios with dark Higgs boson masses within the range of 250 to 400 GeV and Z′ mediators up to 2.3 TeV.

Dark Energy and Dark Matter↗
DOE OSTI2022

Testing and Expertise for Marine Energy (TEAMER) Program Support - Numerical Modeling Assistance for iProTech's "PIP" WEC Device: Cooperative Research and Development Final Report, CRADA Number CRD-20-17303

Develop a time-domain, 3 degree of freedom numerical model of IProTech's PIP device, including power take-off (PTO), that can be used to develop control algorithms and establish a baseline technology performance level (TPL). This is important to IProTech because having a good numerical model allows IProTech to optimize design and control algorithms on an ongoing basis. A baseline TPL enables the PIP technology to be benchmarked against other wave energy converter (WEC) concepts, providing (i) justification (or not) for investment in wave tank model testing, (ii) identifying optimum design parameters and control algorithms for the test model, (iii) specifying test conditions and (iv) interpreting wave tank test results. An accredited WEC-Sim model of the PIP device is an essential tool for further development of the technology. 1) The development of numerical models of the PIP concept will enable iProTech to evaluate the potential performance of this concept, and to investigate how sensitive the device is to different variables within the system (e.g. hydraulic components, geometry dimensions, etc.) 2) The numerical models were developed with WEC-Sim and PTO-Sim -- the de-facto industry standard WEC numerical modeling tools. A detailed model of the system's hydraulic circuit was developed with PTO-Sim. However, WEC-Sim is based on linear hydrodynamics and has several limitations; validation against high fidelity models/physical scale models would help to build confidence in the results.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Search for dark matter produced in association with a single top quark and an energetic W boson in $\sqrt{s}=$ 13 TeV $pp$ collisions with the ATLAS detector

This paper presents a search for dark matter, $χ$, using events with a single top quark and an energetic W boson. The analysis is based on proton–proton collision data collected with the ATLAS experiment at $\sqrt{s}=$ 13 TeV during LHC Run 2 (2015–2018), corresponding to an integrated luminosity of 139 fb -1 . The search considers final states with zero or one charged lepton (electron or muon), at least one $b$-jet and large missing transverse momentum. In addition, a result from a previous search considering two-charged-lepton final states is included in the interpretation of the results. The data are found to be in good agreement with the Standard Model predictions and the results are interpreted in terms of 95% confidence-level exclusion limits in the context of a class of dark matter models involving an extended two Higgs-doublet sector together with a pseudoscalar mediator particle. The search is particularly sensitive to on-shell pro duction of the charged Higgs boson state, $H$ ± , arising from the two-Higgs-doublet mixing, and its semi-invisible decays via the mediator particle, $a$: $H$ ± → $W$ ± $a$(→ $χ$ $χ$). Signal models with $H$ ± masses up to 1.5 TeV and a masses up to 350 GeV are excluded assuming a tan $β$ value of 1. For masses of a of 150 (250) GeV, tan $β$ values up to 2 are excluded for $H$ ± masses between 200 (400) GeV and 1.5 TeV. Signals with tan $β$ values between 20 and 30 are excluded for $H$ ± masses between 500 and 800 GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Modeling the effect of habitat availability and quality on endangered winter-run Chinook salmon ( Oncorhynchus tshawytscha ) production in the Sacramento Valley

Reconciliation between water uses and sustaining aquatic species populations requires an effort to identify and quantify essential habitat characteristics for ecosystem health and valuation methods to predict ecosystem response to restoration actions. This process is particularly challenging for anadromous fish species such as California's Sacramento River winter-run Chinook salmon, due to their limited geographic range and diverse life history habitat requirements. Tools, such as life-cycle models, are needed to manage population dynamics and quantify the composite effects of processes across space and time. Nevertheless, complex institutions can hinder result interpretation and communication, and limit model use in decision-making. This paper focuses on the federally endangered and endemic Sacramento River winter-run Chinook (Oncorhynchus tshawytscha) by developing a Winter-Run Habitat-based Population Model (WRHAP). Additionally, WRHAP is a conceptual, freshwater rearing stage model that includes alternative rearing habitats reported in the literature (e.g., floodplains, off-channel and tributaries), defines rules of habitat use based on instream conditions, and incorporates a juvenile growth module that combines bioenergetics modeling with empirical growth rates. Model outputs reasonably follow observed out-migration patterns and provide a realistic smolt size distribution arriving at the San Francisco Bay. This effort demonstrates the importance of currently “non-critical” habitats (as defined by the Endangered Species Act) for juvenile development, with floodplain habitat contributing to a quarter of out-migrating biomass (despite < 18% availability and <10% of total rearing days), and off-channel growth being one of the most sensitive parameters (explaining ~13% of average juvenile weight variance). The model shows the utility of a simple population model to explore relationships between habitat quality/quantity and juvenile development, and to assist water/environmental management and decision-making processes focused on species recovery.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Analyzing the impact of design factors on solar module thermomechanical durability using interpretable machine learning techniques

Solar modules in utility-scale systems are expected to maintain decades of lifetime to rival conventional energy sources. However, cyclic thermomechanical loading often degrades their long-term performance, highlighting the importance of effective design to mitigate thermal expansion mismatches between module materials. Given the complex composition of solar modules, isolating the impact of individual components on overall durability remains a challenging task. In this work, we analyze a comprehensive data set that comprises bill-of-materials (BOM) and thermal cycling power loss from 251 distinct module designs to identify the predominant design factors and their impacts on the thermomechanical durability of modules. The methodology of our analysis combines machine learning modeling (random forest) and Shapley additive explanation (SHAP) to correlate design factors with power loss and interpret the model’s decision-making. The interpretation reveals that silicon type (monocrystalline or polycrystalline), encapsulant thickness, busbar numbers, and wafer thickness predominantly influence the degradation. With lower power loss of around 0.6% on average in the SHAP analysis, monocrystalline cells present better durability than polycrystalline cells. This finding is further substantiated by statistical testing on our raw data set. The SHAP analysis also demonstrates that while thicker encapsulants lead to reduced power loss, further increasing their thickness over around 0.6 to 0.7 mm does not yield additional benefits, particularly for the front side one. In addition, other important BOM features such as the number of busbars are analyzed. This study provides a blueprint for utilizing explainable machine learning techniques in a complex material system and can potentially guide future research on optimizing the design of solar modules.

14 SOLAR ENERGY↗
DOE OSTI2025

Interpretable ensemble learning unveils main aerosol optical properties in predicting cloud condensation nuclei number concentration

Variations in cloud condensation nuclei number concentration (N CCN ) significantly influence cloud microphysics, yet direct N CCN measurements remain challenging. Here, we present an N CCN ensemble learning (NEL) model utilizing ensemble learning and interpretability analysis on aerosol optical parameters. Validated at two land sites, two ocean sites and one polar site within the Atmospheric Radiation Measurement program, the mean absolute percentage error range of the NEL model across different environments is from 12% to 36%, demonstrating high accuracy. Key findings reveal that aerosol optical parameters can serve as predictors for N CCN . Aerosol scattering and backscattering coefficients, absorption coefficient, backscatter fraction (BSF), and Ångström exponent (AE) are positively correlated with N CCN , while single scattering albedo shows negative correlations. N CCN prediction at land sites is highly sensitive to BSF, largely driven by the backscattering coefficient, as fine particles dominate in these sites. At ocean sites, N CCN prediction is more sensitive to AE, primarily influenced by the scattering coefficient, due to the higher proportion of larger particles. At the polar site, N CCN prediction shows sensitivity to both BSF and AE, mainly driven by the scattering coefficient, as polar sites are cleaner and contain larger particles. These differences reflect the variation in particle size and number concentration across different environments.

Atmospheric science↗
DOE OSTI2026

Structure of high-lying excited states in 47 Ca with ℓ ≥ 3 populated in fast-beam 𝛾-ray-tagged one-neutron pickup reactions

Located just one neutron removed from doubly magic, stable 48 Ca, the adjacent isotope 47 Ca has been studied extensively with neutron-removing transfer reactions and interpreted within a shell-model picture. Neutron-adding transfer reactions onto the stable 46 Ca, however, have not been performed in nearly 60 years and never with the aid of high-resolution 𝛾-ray spectroscopy. Here, we report on the 12 C⁡( 46 , Ca 47 Ca+𝛾)⁢𝑋 one-neutron pickup reaction from 12 C onto 46 Ca at more than 50 MeV/nucleon, a regime of orbital angular momentum mismatch that selectively populates 47 Ca final states based on ℓ ≥ 3 neutron configurations. Newly observed excited states and 𝛾-ray transitions are discussed in comparison to 𝑓⁡𝑝 shell-model calculations in a variety of configuration spaces. Furthermore, it is shown that the nature of the observed spectrum of final states can be well interpreted, albeit with a significant fragmentation of the 𝜈⁢0⁢𝑓 5/2 strength observed in the measurement that is not borne out in the structure calculations.

39 ≤ A ≤ 58↗
DOE OSTI2026

elm-diagnostics

elm-diagnostics is a Python package for computing diagnostic analyses and visualizations for the E3SM Land Model (ELM) component and is meant to support new feature development in ELM. The tool reads model history files and performs quantitative analyses including budget-closure checking, variable transformations, temporal aggregations, and statistical summaries to support model evaluation, validation, and scientific interpretation. The framework is designed for extensibility, with modular architecture enabling straightforward addition of new diagnostic methods, derived variables, analysis types, visualization approaches, and model-specific adaptations

Hoffman, Matt [Los Alamos National Laboratory]↗
DOE OSTI2022

A method for automated modeling of plasma shape control response in tokamaks

Here, a simple method commonly used in modeling of tokamak conductor and plasma response for plasma control has been extended to allow more general modeling of tokamak PF conductor circuits. This extension combines a previously existing method for integrating plasma response and tokamak circuits with a more capable circuit modeling technique derived from the theory that underlies present-day automatic circuit simulation codes. This combined method can incorporate any combination of inductors (e.g., coils), resistors, capacitors, voltage sources, and current sources, while retaining the same model structure and ease of interpretation as the simple method. For many tokamaks, use of this extended method would not be the preferred approach, but for certain devices or for specialized models of a device, it is useful and occasionally essential. For example, the method has enabled fast and automated modeling of the complicated and daily-changing circuit topology used in the DIII-D tokamak. This automated modeling is essential to the support of controller design for the frequently changing requirements for DIII-D plasmas. It is likely to also find application in the workflows for new device designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Discrete-Event Model of WIPP Operations

The Waste Isolation Pilot Plant (WIPP) is the critical component of the Department of Energy's (DOE) Transuranic Radioactive Waste (TRU) disposition infrastructure. Quantifying the operations ofWIPP with a discrete-event model demonstrates the capability to assess efficiency, identity bottlenecks, and improve future operations. Such a model has been developed with the simulation software ExtendSim. This report outlines the structure of that model, summarizes the model's successful reproduction of the annual amount of emplaced waste, and demonstrates that WIPP is successfully receiving and emplacing waste at a rate consistent with the rate at which the waste arrives. The model serves as a first step toward future production enhancements at WIPP. Those enhancements will rely on close collaboration with the Carlsbad Field Office (CBFO), accurate interpretation and incorporation of the model results, and effective planning with other DOE Environmental Management (DOE-EM) entities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Measurement of the inclusive and differential $t\overline{t}γ$ cross sections in the dilepton channel and effective field theory interpretation in proton-proton collisions at $\sqrt{s}$ = 13 TeV

The production cross section of a top quark pair in association with a photon is measured in proton-proton collisions in the decay channel with two oppositely charged leptons (e ± μ ∓ , e + e - , or μ + μ - ). The measurement is performed using 138 fb -1 of proton-proton collision data recorded by the CMS experiment at √s = 13 TeV during the 2016–2018 data-taking period of the CERN LHC. A fiducial phase space is defined such that photons radiated by initial-state particles, top quarks, or any of their decay products are included. An inclusive cross section of 175.2 ± 2.5(stat) ± 6.3(syst) fb is measured in a signal region with at least one jet coming from the hadronization of a bottom quark and exactly one photon with transverse momentum above 20 GeV. Differential cross sections are measured as functions of several kinematic observables of the photon, leptons, and jets, and compared to standard model predictions. The measurements are also interpreted in the standard model effective field theory framework, and limits are found on the relevant Wilson coefficients from these results alone and in combination with a previous CMS measurement of the tt¯γ production process using the lepton+jets final state.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Interpreting strain tensor data to characterize and monitor reservoirs for CO2 storage and other applications

Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and this has opened the door to new opportunities for characterization and monitoring. We have deployed strainmeters and then conducted injection well tests in an underlying reservoir at 530m depth. The resulting data indicated that the horizontal strain at shallow strainmeters (30 to 40m depth) was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 n/d over a few days. The signal at two strainmeters at shallow depth were consistent, although the magnitude of the horizontal strains were different reflecting the different radial directions from the well. The signal at a deep strainmeter deployed at reservoir depth was much different, with tensile vertical strains and compressive horizontal strains. These data can be interpreted by inverting poroelastic forward models developed using numerical and analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. This model is fast and can be inverted to estimate reservoir stiffness and geometry. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. The proxy model is periodically updated and refined using the finite element model to ensure accuracy. This approach significantly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. We have shown that the strain tensor in the caprock is sensitive to pressure in the reservoir, boundaries in the reservoir, and pressure in the caprock caused by leaks. These results indicate that coupling strain tensor data with inversion has the potential to help evaluate reservoirs during initial characterization, and to monitor them during the CO2 injection and storage process.

Murdoch, Larry↗
DOE OSTI2020

Influence of physical and chemical hydrology on bioremediation of a U-contaminated aquifer informed by reactive transport modeling incorporating 238 U/ 235 U ratios

Microbially-catalyzed reductive immobilization of aqueous uranium (U) as a solid phase has been proposed as a U remediation technique. Both laboratory and field experiments have demonstrated that this reduction reaction alters the 238 U/ 235 U ratio, producing a 238 U-enriched U(IV) solid. In contrast, other major U reactive transport processes fractionate these isotopes much less. This suggests the potential to quantify the extent of bioreduction occurring in groundwater containing U using the 238 U/ 235 U ratio, which would substantially improve upon current practices largely relying on U concentration measurements alone. Current reactive transport models for uranium dynamics include only concentration measurements, which are strongly influenced by highly coupled reactive transformation and aqueous transport processes. The complex physical and chemical behavior of U potentially compromises such quantitative analysis of U storage and release. Here we report the first numerical reactive transport model which explicitly incorporates variations in the 238 U/ 235 U ratio of U and demonstrates improved interpretation of the principal chemical reactions and groundwater transport processes affecting the subsurface mobility and distribution of this widespread contaminant. Recent U bioreduction studies performed in a contaminated aquifer in Rifle, Colorado, USA applied Rayleigh distillation models to interpret U stable isotope fractionation observed as a result of acetate amendment. These simplified models were unable to resolve the spatiotemporal pattern of U isotope fractionation recorded in the aqueous solutes. Here, we employ the multi-component, isotope-enabled CrunchTope reactive transport software to interpret these measured U isotope ratios, and demonstrate accurate reproduction of observed trends in both geochemistry and 238 U/ 235 U ratios for two consecutive years of field experiments. Overall, our results indicate that accurately modeling both the U concentration and isotope ratio distributions greatly constrains the parameter space of the model. We find that the transport properties of U in the Rifle aquifer are governed by the presence of low-permeability regions, which the isotopes are uniquely sensitive to. When U reduction is spatially constrained by these low-permeability regions, the shift in the 238 U/ 235 U ratio becomes more muted. Accurate modeling of observed U isotope ratios thus provides a powerful means to better understand bioremediation, and the current study serves to advance the application of this novel method.

238U/235U↗
DOE OSTI2023

Autonomous Inversion of In Situ Deformation Measurement Data for Injection-Induced Stress Change

Geologic carbon storage (GCS) is likely to play a key part of the global effort to dramatically reduce CO2 emissions and perhaps even reduce atmospheric CO2 concentrations through carbon negative operations. A critical part of effort to commercialize and widely deploy this technology is developing the capability to rapidly assimilate real-time monitoring data into a form that will enable site operators to make decisions to manage the safe and efficient operations. Two of the risks associate with GCS are the risk of inducing fractures in the sealing formations that can create leakage pathways and the risk of inducing earthquakes of sufficient magnitude to cause public concern, property damage, or safety risks. To properly manage these risks the site operator needs to know the initial state of stress, the change in stress induced by injection, and the relationship between operational parameters such as injection rate and pressure and the change in stress. Current methods of estimating the change in stress require choosing the type of constitutive model and the model parameters based on core, log, and geophysical data during the characterization phase, with little feedback from operational observations to validate or refine these choices. These characterization methods interrogate the geologic formations using length scales, loading rates or magnitudes that are quite different from those encountered by the actual storage system. It is shown that errors in the assumed constitutive response, even when informed by laboratory tests on core samples, are likely to be common, large, and underestimate the magnitude of stress change caused by injection. Recent advances in borehole-based strain instruments and borehole and surface-based tilt and displacement instruments have now enabled monitoring of the deformation of the storage system throughout its operational lifespan. This data can enable validation and refinement of the knowledge of the geomechanical properties and state of the system, but brings with it a challenge to transform the raw data into actionable knowledge. We demonstrate a method that uses automatic differentiation and a finite-element based geomechanical model perform a gradient-based deterministic inversion of geomechanical monitoring data. This approach allows autonomous integration of the instrument data without the need for time consuming manual interpretation and selection of updated model parameters. Furthermore, only isotropic linear elasticity is considered in this paper, the approach presented is very flexible as to what type of geomechanical constitutive response can be used. The approach is easily adaptable to nonlinear physics-based constitutive models to account for common rock behaviors such as creep and plasticity. The approach also enables training of machine learning-based constitutive models by allowing back propagation of errors through the finite element calculations. This enables strongly enforcing known physics, such as conservation of momentum and continuity, while allowing data-driven models to learn the truly unknown physics such as the constitutive or petrophysical responses.

Burghardt, Jeffrey A.↗
DOE OSTI2021

Inter-well connectivity detection in CO 2 WAG projects using statistical recurrent unit models

Routine well-wise injection and production measurements contain significant information on subsurface structure and properties. Data-driven technology that interprets surface data into subsurface structure or properties can assist operators in making informed decisions by providing a better understanding of field assets. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO 2 EOR projects utilizing the water-alternating-gas (WAG) process. SRU is a special type of recurrent neural network (RNN) that allows for better characterization of temporal trends, by learning various statistics of the input at different time scales. In our application, the complete states (injection rate, pressure and cumulative injection) at injectors and pressure states at producers are fed to SRU as the input and the phase rates at producers are treated as the output. Once the SRU is trained and validated, it is then used to assess the connectivity of each injector to any producer using permutation variable importance method, wherein inputs corresponding to an injector are shuffled and the increase in prediction error at a given producer is recorded as the importance (connectivity metric) of the injector to the producer. This method is tested in both synthetic and field-scale cases. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. This significantly improves confidence in our data-driven procedure. The novelty of this work is that it is purely data-driven method and can directly interpret routine surface measurements to intuitive subsurface knowledge. Furthermore, the streamline-based validation procedure provides physics-based backing to the results obtained from data analytics. This study results in a reliable and efficient data analytics framework that is well-suited for large field applications.

42 ENGINEERING↗
DOE OSTI2023

Predictive understanding of the surface tension and velocity of sound in ionic liquids using machine learning

Knowledge of the physical properties of ionic liquids (ILs), such as the surface tension and speed of sound, is important for both industrial and research applications. Unfortunately, technical challenges and costs limit exhaustive experimental screening efforts of ILs for these critical properties. Previous work has demonstrated that the use of quantum-mechanics-based thermochemical property prediction tools, such as the conductor-like screening model for real solvents, when combined with machine learning (ML) approaches, may provide an alternative pathway to guide the rapid screening and design of ILs for desired physiochemical properties. However, the question of which machine-learning approaches are most appropriate remains. In the present study, we examine how different ML architectures, ranging from tree-based approaches to feed-forward artificial neural networks, perform in generating nonlinear multivariate quantitative structure–property relationship models for the prediction of the temperature- and pressure-dependent surface tension of and speed of sound in ILs over a wide range of surface tensions (16.9–76.2 mN/m) and speeds of sound (1009.7–1992 m/s). The ML models are further interrogated using the powerful interpretation method, shapley additive explanations. We find that several different ML models provide high accuracy, according to traditional statistical metrics. The decision tree-based approaches appear to be the most accurate and precise, with extreme gradient-boosting trees and gradient-boosting trees being the best performers. However, our results also indicate that the promise of using machine-learning to gain deep insights into the underlying physics driving structure–property relationships in ILs may still be somewhat premature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Search for heavy Majorana or Dirac neutrinos and right-handed $W$ gauge bosons in final states with charged leptons and jets in $pp$ collisions at $\sqrt{s}=13$ TeV with the ATLAS detector

A search for heavy right-handed Majorana or Dirac neutrinos N R and heavy right-handed gauge bosons W R is performed in events with energetic electrons or muons, with the same or opposite electric charge, and energetic jets. The search is carried out separately for topologies of clearly separated final-state products (“resolved” channel) and topologies with boosted final states with hadronic and/or leptonic products partially overlapping and reconstructed as a large-radius jet (“boosted” channel). The events are selected from pp collision data at the LHC with an integrated luminosity of 139 fb -1 collected by the ATLAS detector at $\sqrt{s}=13$ TeV. No significant deviations from the Standard Model predictions are observed. The results are interpreted within the theoretical framework of a left-right symmetric model, and lower limits are set on masses in the heavy right handed W R boson and N R plane. The excluded region extends to about m(W R ) = 6.4 TeV for both Majorana and Dirac N R neutrinos at m(N R ) < 1 TeV. N R with masses of less than 3.5 (3.6) TeV are excluded in the electron (muon) channel at m(W R ) = 4.8 TeV for the Majorana neutrinos, and lim its of m(N R ) up to 3.6 TeV for m(W R ) = 5.2 (5.0) TeV in the electron (muon) channel are set for the Dirac neutrinos. These constitute the most stringent exclusion limits to date for the model considered.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Potential Formation Damage: An Integrated Reservoir Characterization Study of the Naturally Fractured Carbonate Middle Duperow Formation at the Kevin Dome, Montana

In this study, we integrate geologic and engineering data of a naturally fractured carbonate reservoir at the Kevin Dome, Montana. Well test data are correlated with core description, geochemical and lithology study to determine the flow behavior and communication within the injection test interval and to the surrounding area. Based on a dual-continuum geologic model, numerical brine injection simulations are carried out to validate the interpretation results from our well test analytical models and forecast the probability of CO 2 injection success using current reservoir properties. As a result, our well test analytical models as well as lithology/core description suggest that fluid flow may be mainly restricted to the injection interval and the assumption of radial (horizontal) flow may be appropriate. The well test models also indicate that there is potentially formation damage with a positive skin factor although prior to brine injection well tests, well stimulation through acid treatment was performed. Our numerical simulation results appear to confirm this formation damage by showing additional pressure buildup in the injection data during later test periods. To explain this, acid may have dissolved dolomite then dolomite or calcite may have been formed again further into the matrix/fracture system. Another possible explanation is mechanical clogging of the fractures due to acid dissolving dolomite and dislodging fine grains. Our work also predicts that if no additional well stimulation is performed, the project will have a lower probability of successfully injecting 1 million tons of CO 2 into the Middle Duperow formation over 4 years.

42 ENGINEERING↗