SEARCH · Engineering Papers

Results for “rock finding”

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 181 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Nanoconfinement matters in humidified CO 2 interaction with metal silicates

With enigmatic observations of enhanced reactivity of wet CO 2 -rich fluids with metal silicates, the mechanistic understanding of molecular processes governing carbonation proves critical in designing secure geological carbon sequestration and economical carbonated concrete technologies. Here, we use the first principle and classical molecular simulations to probe the impact of nanoconfinement on physicochemical processes at the rock–water–CO 2 interface. We choose nanoporous calcium–silicate–hydrate (C–S–H) and forsterite (Mg 2 SiO 4 ) as model metal silicate surfaces that are of significance in cement chemistry and geochemistry communities, respectively. We show that while a nanometer-thick interfacial water film persists at unsaturated conditions consistent with in situ infrared spectroscopy, the phase behavior of the water–CO 2 mixture changes from its bulk counterpart depending on the surface chemistry and nanoconfinement. We also observe enhanced solubility at the interface of water and CO 2 phases, which could amplify the CO 2 speciation rate. Additionally, through free energy calculations, we show that CO 2 could be found in a metastable state near the C–S–H surface, which can potentially react with surface water and hydroxyl groups to form carbonic acid and bicarbonate. These findings support the explicit consideration of nanoconfinement effects in reactive and non-reactive pore-scale processes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Integrated Modeling Report for CO 2 Storage with Enhanced Oil Recovery in Northern Michigan

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership initiative. Since its founding, MRCSP has made significant strides toward making CCUS a viable option for states in the region. The public/private consortium, funded through the DOE Regional Carbon Sequestration Initiative, brings together nearly 40 industry partners and 10 states. Battelle, as the project lead, oversees research, development and operations and coordinates activities among the partners. The incremental, phased approach has built a valuable knowledge base for the industry and paved the way for commercial-scale adoption of CCUS technologies. From 2008 to 2020, MRCSP Phase III focused on the development of large-scale injection projects. This report is part of a series of reports prepared under the Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase). These reports summarize and detail the findings of the work conducted under the Phase III project. MRCSP has been investigating various reservoir characterization and modeling technologies related to Carbon Capture, Utilization, and Storage (CCUS) in conjunction with carbon dioxide enhanced oil recovery (CO 2 -EOR) operations in multiple Silurian-age (Niagaran), oil-bearing, carbonate pinnacle reefs in northern Michigan, USA. This report provides a comprehensive discussion of reservoir modeling studies that were conducted for tracking oil production, forecasting CO 2 plume migration, and estimating associated storage in a number of reefs that were at different stages of their CO 2 -EOR life cycle. The modeling process for simulating oil production, CO 2 injection, and associated storage in these reefs entails two phases. The first phase, geologic framework modeling, integrates all pertinent geological and geophysical data (from logs, cores and seismic surveys) about reservoir structure, geometry, rock types, and property distributions (porosity, permeability, water saturation) into a 3-D distributed grid-based static earth model (SEM). The second phase, dynamic reservoir modeling, uses the SEM as a platform to simulate the movement of oil, gas, water, and CO 2 within the reservoir during primary hydrocarbon production, as well as during subsequent phases such as CO 2 -injection assisted EOR, plume migration, and associated storage. In addition, an assessment of coupled process effects is also carried out, where the impacts of geochemical and geo-mechanical processed induced by CO 2 injection are studied.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2023

An adjoint-based optimization method for jointly inverting heterogeneous material properties and fault slip from earthquake surface deformation data

SUMMARY Analysis of tectonic and earthquake-cycle associated deformation of the crust can provide valuable insights into the underlying deformation processes including fault slip. How those processes are expressed at the surface depends on the lateral and depth variations of rock properties. The effect of such variations is often tested by forward models based on a priori geological or geophysical information. Here, we first develop a novel technique based on an open-source finite-element computational framework to invert geodetic constraints directly for heterogeneous media properties. We focus on the elastic, coseismic problem and seek to constrain variations in shear modulus and Poisson’s ratio, proxies for the effects of lithology and/or temperature and porous flow, respectively. The corresponding nonlinear inversion is implemented using adjoint-based optimization that efficiently reduces the cost function that includes the misfit between the calculated and observed displacements and a penalty term. We then extend our theoretical and numerical framework to simultaneously infer both heterogeneous Earth’s structure and fault slip from surface deformation. Based on a range of 2-D synthetic cases, we find that both model parameters can be satisfactorily estimated for the megathrust setting-inspired test problems considered. Within limits, this is the case even in the presence of noise and if the fault geometry is not perfectly known. Our method lays the foundation for a future reassessment of the information contained in increasingly data-rich settings, for example, geodetic GNSS constraints for large earthquakes such as the 2011 Tohoku-oki M9 event, or distributed deformation along plate boundaries as constrained from InSAR.

Geochemistry & Geophysics↗

DOE OSTI2024

Report series: finding of effect and mitigation documentation for the mercury solar photovoltaic array and battery energy storage system, area 23, nevada national security site, nye county, nevada

The U.S. Department of Energy (DOE), National Nuclear Security Administration Nevada Field Office (NNSA/NFO) proposes to install solar photovoltaic (PV) power generation arrays and an associated battery energy storage system (BESS) for the town of Mercury at the Nevada National Security Site (NNSS) in Nye County, Nevada (Figure 1). The purposes of the development of this facility are to support long-term efforts to modernize Mercury and to provide energy-resilient infrastructure and address climate adaptation needs. Because it is within the boundary of the Mercury Historic District (MHD), it is subject to the terms of the Programmatic Agreement Between the National Nuclear Security Administration Nevada Field Office and the Nevada State Historic Preservation Officer Regarding Modernization and Operational Maintenance of the Nevada National Security Site, at Mercury in Nye County, Nevada (hereafter referred to as the Mercury PA). An identification and evaluation report prepared for this undertaking determined that a contributing element to the MHD, the Mercury airstrip (26NY15777), is within the APE (Haynes 2024). The Mercury airstrip was developed following the closure of Camp Desert Rock in 1957 and used until late 1963 or 1964 when the Camp Desert Rock Airport was renovated, and the Mercury Bypass road constructed. The town of Mercury and the immediate surrounding area have been determined eligible for listing in the National Register of Historic Places (NRHP) as the MHD (SHPO Resource No. D230) under Criteria A and C for its importance in supporting nuclear testing and scientific research from 1951 through 1992 (Reed 2019). Originally recorded in 2016, 26NY15777 was determined individually not eligible for listing in the NRHP because it lacked sufficient integrity to convey its significance (Palmer 2016). It was subsequently determined in 2018 to be a contributing element of the MHD in an architectural survey of the district (Palmer 2018) and identified in Appendix C of the Mercury PA as a Category I contributing element. However, as per Stipulation IV.B.2, because the airstrip had already been formally evaluated and determined not to be individually eligible for the NRHP in consultation with the SHPO, this categorization was an error. NNSA/NFO reported this to the SHPO in Haynes 2024 and the SHPO concurred on June 11, 2024 (Reed). Accordingly, the airstrip is a Category II Property for the purposes of complying with the Mercury PA. The Mercury airstrip is a historic property for the purposes of compliance with Section 106 of the National Historic Preservation Act (NHPA) and subject to the stipulations of the Mercury PA.

25 ENERGY STORAGE↗

DOE OSTI2023

Radionuclide Transport Simulations Supporting Proposed Borehole Waste Disposal in Israel

A scientific collaboration between the U.S. and Israel is underway to assess the suitability of a potential site for subsurface radioactive waste disposal in the Negev Desert, Israel. The Negev Desert has several favorable attributes for geologic disposal, including an arid climate, a deep vadose zone, interlayered low-permeability lithologies, and carbonate rocks with high uranium-sorption potential. These features may provide a robust natural barrier to radionuclide migration. Geologic and laboratory characterization data from the Negev Desert are incorporated into multiphase flow and transport models, solved using PFLOTRAN, to aid in site characterization and risk analysis that will support decision-making for waste disposal in an intermediate-depth borehole design. The lithology with the greatest uranium sorption potential at the site is phosphorite. We use modeling to evaluate the ability of this layer to impact uranium transport around a proposed disposal borehole. The current objective of the simulations is focused on characterizing hypothetical leakage from waste canisters and subsequent uranium migration under three infiltration scenarios. Here, we describe a hydrogeologic model based on data from a local exploratory borehole and present results for uranium flow and transport simulations under varying infiltration scenarios. We find that under the current climate conditions, it is likely that uranium will remain in the near-field of the borehole for thousands of years. However, under a hypothesized extreme climate scenario representing an increase in infiltration by a factor of 300x above present-day values, uranium may break through the phosphorite layer and exit the base of the model domain (~200 m above the water table) within 1000 years. Simulation results have direct implications for the planning of nuclear waste disposal in the Negev Desert, and specifically in intermediate-depth boreholes.

58 GEOSCIENCES↗

DOE OSTI2022

Utah FORGE Well 16A(78)-32 Stimulation DFN Fracture Plane Evaluation and Data

This dataset includes files used to fit planar fractures through the preliminary earthquake catalogs of the three stages of the April 2022 well 16A(78)-32 stimulation which is linked bellow. These planar features have been used to update the FORGE reference Discrete Fracture Network (DFN) model. The files are provided to encourage other modelers to use additional workflows to find additional/alternative features. To this end, the dataset includes the cleaned earthquake catalog data translated to the FORGE reference model global reference frame, the well trajectory of 16A(78)-32 in those same coordinates, the fit 15 planar features in csv format, and a pdf file with slides illustrating the process used to fit the features. A recorded presentation of this material is available from the October 2022 FORGE Modeling and Simulation Forum which is also linked below.

15 GEOTHERMAL ENERGY↗

DOE OSTI2023

Organic matter sulfurization and organic carbon burial in the Mesoproterozoic

Throughout the Proterozoic Era, sedimentary organic carbon burial helped set the pace of global oxygenation and acted as a major modulator of atmospheric CO 2 and climate. Although Proterozoic rocks generally contain low concentrations of organic matter (OM), there are key exceptions to this rule, including the relatively OM-rich Arctic Bay shales from Baffin Island, Canada (Bylot Supergroup, Borden Basin, ~1.05 Ga). The mechanisms driving elevated OM concentrations in these and other Proterozoic shales remain poorly understood. In the Mesozoic and Cenozoic, organic matter sulfurization can be a major driver of enhanced OM burial across a range of redox conditions comparable to those inferred for many Proterozoic environments. Therefore, in this study, we evaluate the role of sulfurization in driving OM preservation in the Mesoproterozoic Borden Basin and discuss its relevance to Proterozoic systems in general. We present the first evidence for syngenetic-to-early-diagenetic OM sulfurization in a Proterozoic basin, which begins to fill a several-billion-year gap in our record of organic S across Earth history. We find that OM sulfurization was particularly extensive in shales from a relatively shallow-water section (Alpha River) but less extensive in shales deposited in deeper water (Shale Valley), which is consistent with models that infer sulfidic ‘wedges’ or O 2 -minimum-zone-type structures on shelf margins at least intermittently at this time. At the shallower site, organic S and pyrite are similarly 34 S-depleted and thus likely formed at roughly the same time near the sediment–water interface under conditions previously interpreted to have been ferruginous to intermittently sulfidic. In contrast, at the deeper-water site, large S-isotope differences between pyrite and organic S along with low apparent OM sulfurization intensities indicate that pyrite formation was favored over OM sulfurization during early sedimentation under variable but primarily ferruginous conditions. Although Mesoproterozoic biomass can be substantially sulfurized, indicators of sulfurization intensity are not correlated with OM concentrations, and therefore sulfurization does not appear to have been the primary driver of enhanced OM concentrations in Arctic Bay Formation shales. The link between sulfurization and total OM preservation may have been modulated during the deposition of Arctic Bay Formation shales by differences in iron availability, nutrient cycling, and particle dynamics in the Mesoproterozoic.

58 GEOSCIENCES↗

DOE OSTI2026

A parametric study of slow dynamic nonlinear elasticity with comparisons to models

Several phenomenological models that aspire to quantitative description of anomalous nonlinear mesoscopic elasticity are reviewed and compared with laboratory measurements. This class of nonlinearity, best known perhaps for slow dynamics and aging, is seen widely in imperfectly consolidated granular solids but is not well understood. Typical slow dynamic tests show that a modest conditioning oscillatory "pump" strain depresses material stiffness, which then recovers like the logarithm of time after conditioning ceases. Several phenomenological models based on physical arguments have been proposed that predict the material stiffness response to arbitrary pump strain histories during conditioning and recovery. Approximate closed form and numerical solutions to the models are presented that predict the quantitative influence of three key pump parameters: the pump's strain amplitude, the pump's strain rate, and the pump’s duration. Laboratory measurements on Berea sandstone, concrete and a confined single aluminum bead find that slow dynamic responses are linear in pump strain and independent of pump frequency. Measurements also show that, after pump-off, stiffness recovers over times far longer than the pump duration. These observations and others are compared to model predictions. One of the considered models, based on a picture of fast brittle damage and slow healing, successfully matches all these behaviors.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Influence of curing temperature on belite cement hydration: A comparative study with Portland cement

Belite cements (BCs) could be a more sustainable binder than Portland cements (PCs) but adequate knowledge of the hydration features has still to be built. In particular, the mild curing hydration temperature effects have been extensively studied for PCs but not for BCs. This research was triggered by a previous work reporting improved mechanical strengths of BCs at higher curing temperatures. Here, we report the hydration characteristics of a BC at 20, 40 and 60 °C and compared to those of a typical PC. We have corroborated previous findings and used a multi-technique approach including Rietveld phase analysis, thermal analysis, calorimetry, silicon MAS-NMR, mercury intrusion porosimetry and chiefly synchrotron X-ray microtomography, to thoroughly understand the different behavior. In a nutshell, the improved mechanical performances at mild curing temperatures for BCs are mainly due to a much larger belite degree of reaction, with lower porosity coarsening of belite cements also playing a role.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Coarse-grained simulation of colloidal self-assembly, cation exchange, and rheology in Na/Ca smectite clay gels

Knowledge Gap: The aggregation of clay minerals—layered silicate nanoparticles—strongly impacts fluid flow, solute migration, and solid mechanics in soils, sediments, and sedimentary rocks. Experimental and computational characterization of clay aggregation is inhibited by the delicate water-mediated nature of clay colloidal interactions and by the range of spatial scales involved, from 1 nm thick platelets to flocs with dimensions up to micrometers or more. Simulations: Using a new coarse-grained molecular dynamics (CGMD) approach, we predicted the microstructure, dynamics, and rheology of hydrated smectite (more precisely, montmorillonite) clay gels containing up to 2,000 clay platelets on length scales up to 0.1 μm. Further, simulations investigated the impact of simulation time, platelet diameters (6 to 25nm), and the ratio of Na to Ca exchangeable cations on the assembly of tactoids (i.e., stacks of parallel clay platelets) and larger aggregates (i.e., assemblages of tactoids). We analyzed structural features including tactoid size and size distribution, basal spacing, counterion distribution in the electrical double layer, clay association modes, and the rheological properties of smectite gels. Findings: Our results demonstrate new potential to characterize and understand clay aggregation in dilute suspensions and gels on a scale of thousands of particles with explicit representation of counterion clouds and with accuracy approaching that of all-atom molecular dynamics (MD) simulations. For example, our simulations predict the strong impact of Na/Ca ratio on clay tactoid formation and the shear-thinning rheology of clay gels.

42 ENGINEERING↗

DOE OSTI2024

Optimizing Deep Geothermal Drilling for Energy Sustainability in the Appalachian Basin

This study investigates the geological and geomechanical characteristics of the MIP 1S geothermal well in the Appalachian Basin to optimize drilling and address the wellbore stability issues encountered. Data from well logs, sidewall core analysis, and injection tests were used to derive elastic and rock strength properties, as well as stress and pore pressure profiles. A robust 1D-geomechanical model was developed and validated, correlating strongly with wellbore instability observations. This revealed significant wellbore breakout, widening the diameter from 12 ¼ inches to over 16 inches. Advanced technologies like Cerebro Force™ In-Bit Sensing were used to monitor drilling performance with high accuracy. This technology tracks critical metrics such as bit acceleration, vibration in the x, y, and z directions, Gyro RPM, stick-slip indicators, and bending on the bit. Cerebro Force™ readings identified hole drag caused by poor hole conditions, including friction between the drill string and wellbore walls and the presence of cuttings or debris. This led to higher torque and weight on bit (WOB) readings at the surface compared to downhole measurements, affecting drilling efficiency and wellbore stability. Optimal drilling parameters for future deep geothermal wells were determined based on these findings.

Environmental Sciences & Ecology↗

DOE OSTI2025

Enhanced Weathering May Benefit From Co‐Application With Organic Amendments

Abstract Enhanced weathering has emerged as a promising natural climate solution that has the potential to remove billions of tons of carbon from the atmosphere if widely adopted in agricultural settings. Despite this potential, few field trials have been published that verify the carbon dioxide removal (CDR) potential of enhanced weathering in croplands and, until now, none had been published in grazing lands. Anthony et al. (2025, https://doi.org/10.1029/2024AV001480 ) conducted the first trial of enhanced weathering in a California rangeland and showed weathering of ground silicate rocks despite drought conditions in an already dry climate. Co‐application of inorganic (silicate rocks) with organic (biochar and compost) amendments revealed not just additive, but synergistic effects whereby organic amendments increased rates of weathering. This is important because field CDR rates were <10% of the theoretical maximum (i.e., the rate if basalt was completely weathered); thus, methods to improve weathering rates will be necessary for this practice to scale in a meaningful way. Multi‐carbon pool measurements revealed not only how co‐application of soil amendments heightened net carbon benefits, but also how soil amendments complemented each other to produce net benefits for soil carbon, biomass growth, and greenhouse gas emission reductions. Anthony et al. (2025, https://doi.org/10.1029/2024AV001480 ) produce new insights toward our understanding of enhanced weathering as well as introduce paths for future research concerning combined amendment applications, synergistic mechanisms for carbon storage, and deployment in various agricultural contexts. While questions remain about the fate of weathering products in arid regions, Anthony et al. (2025, https://doi.org/10.1029/2024AV001480 ) present novel findings on the potential for significant weathering to occur even under suboptimal conditions.

Almaraz, Maya [Yale Center for Natural Carbon Capt↗

DOE OSTI2022

Seeking Repeating Anthropogenic Seismic Sources: Implications for Seismic Velocity Monitoring at Fault Zones

Abstract Seismic velocities in rocks are highly sensitive to changes in permanent deformation and fluid content. The temporal variation of seismic velocity during the preparation phase of earthquakes has been well documented in laboratories but rarely observed in nature. It has been recently found that some anthropogenic, high‐frequency (>1 Hz) seismic sources are powerful enough to generate body waves that travel down to a few kilometers and can be used to monitor fault zones at seismogenic depth. Anthropogenic seismic sources typically have fixed spatial distribution and provide new perspectives for velocity monitoring. In this work, we propose a systematic workflow to seek such powerful seismic sources in a rapid and straightforward manner. We tackle the problem from a statistical point of view, considering that persistent, powerful seismic sources yield highly coherent correlation functions (CFs) between pairs of seismic sensors. The algorithm is tested in California and Japan. Multiple sites close to fault zones show high‐frequency CFs stable for an extended period of time. These findings have great potential for monitoring fault zones, including the San Jacinto Fault and the Ridgecrest area in Southern California, Napa in Northern California, and faults in central Japan. However, extra steps, such as beamforming or polarization analysis, are required to determine the dominant seismic sources and study the source characteristics, which are crucial to interpreting the velocity monitoring results. Train tremors identified by the present approach have been successfully used for seismic velocity monitoring of the San Jacinto Fault in previous studies.

58 GEOSCIENCES↗

DOE OSTI2020

Connecting Hindered Transport in Porous Media across Length Scales: From Single-Pore to Macroscopic

Hindered mass transport is widely observed in various porous media, including filtration membranes and porous rocks. However, there is no universal model capable of predicting transport in porous media due to the heterogeneity of porous structures and the complexity of the underlying microscopic mechanisms, including electrostatic and hydrodynamic effects. Here, we used a highly ordered porous medium (i.e. a silica inverse opal film) as a model system to directly explore the effects of geometric parameters (i.e. pore size, pore throat size and tracer particle size) and microscopic interaction parameters (e.g. controlled by ionic strength) on nanoparticle transport in porous environments using three-dimensional single-particle tracking. We found a linear scaling relation between the macroscopic diffusion coefficient and microscopic diffusion behaviors and a combination of parameters associated with microscopic features and phenomena, including both geometric effects and particle-level interactions. The proportionality coefficient relating micro and macro behaviors was complex and related to the connectivity of the matrix and pore-size variation, which could lead to tortuous diffusion pathways, hindering macroscopic transport. Furthermore, these findings provide important insights into nanoparticles and macromolecules transport in various porous media.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

A dynamic microbial sulfur cycle in a serpentinizing continental ophiolite

Summary Serpentinization is the hydration and oxidation of ultramafic rock, which occurs as oceanic lithosphere is emplaced onto continental margins (ophiolites), and along the seafloor as faulting exposes this mantle‐derived material to circulating hydrothermal fluids. This process leads to distinctive fluid chemistries as molecular hydrogen (H 2 ) and hydroxyl ions (OH − ) are produced and reduced carbon compounds are mobilized. Serpentinizing ophiolites also serve as a vector to transport sulfur compounds from the seafloor onto the continents. We investigated hyperalkaline, sulfur‐rich, brackish groundwater in a serpentinizing continental ophiolite to elucidate the role of sulfur compounds in fuelling in situ microbial activities. Here we illustrate that key sulfur‐cycling taxa, including Dethiobacter , Desulfitispora and ‘ Desulforudis ’, persist throughout this extreme environment. Biologically catalysed redox reactions involving sulfate, sulfide and intermediate sulfur compounds are thermodynamically favourable in the groundwater, which indicates they may be vital to sustaining life in these characteristically oxidant‐ and energy‐limited systems. Furthermore, metagenomic and metatranscriptomic analyses reveal a complex network involving sulfate reduction, sulfide oxidation and thiosulfate reactions. Our findings highlight the importance of the complete inorganic sulfur cycle in serpentinizing fluids and suggest sulfur biogeochemistry provides a key link between terrestrial serpentinizing ecosystems and their submarine heritage.

Sabuda, Mary C.↗

DOE OSTI2025

Discrete fracture network model benchmarks developed and applied in a DECOVALEX-2023 repository performance assessment study

This study presents newly developed benchmarks for modeling flow and transport within discrete fracture networks (DFNs) and useful methods for analyzing the results. The new benchmarks are designed to test modeling approaches for use in probabilistic performance assessment models of deep geologic repositories in fractured rock. The benchmarks simulate flow and transport through a 1 km 3 block of fractured rock. The first simulates migration of a short pulse of tracer through a simple network of four intersecting fractures. The second adds 1089 stochastically generated fractures. The third changes the pulse to a continuous point source. Evaluation of model performance relies on moment analysis and comparison of the results of different models. The expected nondimensional first moment of the conservative tracer for each benchmark is 1. The benchmarks were simulated by teams from Canada, Czechia, Germany, Korea, Sweden, Taiwan, and the United States as part of a DECOVALEX-2023 study (decovalex.org). The teams used various approaches, including explicit DFN modeling, DFN upscaling to an equivalent continuous porous medium (ECPM), and a combination of both methods. Transport mechanisms are modeled using either the advection-dispersion equation or particle tracking. Results demonstrate strong agreement among the models in breakthrough behavior up to the 75th percentile. Significant deviations in first moments and well-clustered outputs led to the identification of inaccuracies in several models. Such findings exemplify the benefit of exercising these benchmarks and using the presented methods to test DFN flow and transport models.

Benchmark↗

DOE OSTI2025

Development of an ERT‐Based Framework for Bentonite Buffers Monitoring From Laboratory Tests: 1. Characterizing Thermal–Hydrological–Mechanical Processes

Abstract Bentonite clay is widely used in engineered barrier systems for the permanent disposal of high‐level radioactive waste due to its low permeability, high swelling capacity, and thermal stability. However, the complex thermal‐hydrological‐mechanical (THM) processes induced by heating from decaying radioactive waste and hydration from surrounding rock can lead to heterogeneous changes that are difficult to measure and predict. This study develops an Electrical Resistivity Tomography (ERT)‐based framework for monitoring THM processes, progressing from sample‐scale to bench‐scale tests, to inform field‐scale applications. Sample‐scale tests analyzed small bentonite samples under controlled variations in water content, temperature, and porosity to establish fundamental resistivity relationships. Bench‐scale tests involved larger bentonite columns subjected to heating (up to 200°C) and hydration under controlled pressure, simulating repository conditions. ERT measurements, complemented by X‐ray CT imaging, temperature monitoring, and tracing sensors, revealed coupled THM processes, such as hydration‐induced compression, swelling, and thermal gradients, leading to complex resistivity patterns. The results demonstrate the potential of ERT for capturing THM‐induced resistivity changes, though challenges remain in upscaling and quantitative analysis. This study evaluates laboratory test capabilities and proposes future improvements for understanding THM‐induced resistivity responses. A conceptual framework for ERT implementation in field‐scale monitoring is presented, synthesizing findings from both scales and exploring how ERT data can inform long‐term modeling and reduce prediction uncertainties. Overall, this ERT‐based framework offers a robust method for monitoring bentonite buffers, aiding in early issue detection and supporting the safe long‐term disposal of radioactive waste in geological repositories, while highlighting the need for future development. Plain Language Summary Bentonite clay is crucial in engineered barrier systems (EBS) for containing high‐level radioactive waste due to its ability to absorb water, swell, seal and remain stable under high temperatures. When bentonite absorbs water and heats up from radioactive decay, it experiences complex changes in its physical and mechanical properties. Understanding these changes is important for ensuring the long‐term safety and effectiveness of EBS. This study used Electrical Resistivity Tomography (ERT), a non‐invasive method that measures electrical conductivity to monitor these changes during laboratory experiments. The ERT data revealed significant variations in resistivity corresponding to changes in water content, temperature, and density, providing detailed spatial and temporal insights into the behavior of bentonite. These findings enhance our ability to predict the long‐term performance of bentonite barriers, ensuring the safe containment of radioactive waste. By improving our understanding of bentonite's behavior, this research supports the development of more reliable and effective barrier systems for radioactive waste disposal, protecting the environment and public health. Key Points ERT monitoring was employed to capture resistivity changes in bentonite during controlled heating and hydration experiments, providing insights into THM processes ERT data reveal significant resistivity changes correlated with water content, temperature, and mechanical effects, enhancing the understanding of THM dynamics in bentonite This study explores the potential of the framework for application in field‐scale EBS monitoring, emphasizing the need for integrating additional geophysical methods for comprehensive subsurface imaging

Chen, Hang↗

DOE OSTI2024

Fe Kα XANES, Fe Kβ HERFD XANES and EPMA flank method determinations of the oxidation state of Fe in garnet

The ferric to total iron ratios (Fe 3+ /ΣFe) of garnets can be paired with thermodynamic mineral activity models to quantify the oxygen fugacity of garnet-bearing rocks. However, techniques with a high analytical and spatial resolution are necessary to distinguish differences in garnet Fe 3+ /ΣFe ratios at the percent level and to accurately measure garnets that are zoned or contain inclusions. We acquired conventional Fe Kα and high-resolution energy fluorescence detection (HERFD) Fe Kβ X-ray absorption near edge structure (XANES) spectra and electron microprobe flank method analyses on a suite of 27 peridotitic and eclogitic garnets with Fe 3+ /ΣFe ratios previously determined by Mössbauer spectroscopy to evaluate the precision of each technique. We examined variations in the energy and intensity of three XANES spectral features as a function of Fe 3+ /ΣFe ratios: 1) the intensity ratio of two-post edge features (I-ratio; Fe Kα only); 2) the energy of the Fe edge at 90% normalized intensity (E 0.9 ; Fe Kα only) and 3) the pre-edge centroid energy (Fe Kα and HERFD Fe Kβ). In accordance with previous work, we find the energies of garnet pre-edge centroids are relatively insensitive to Fe 3+ /ΣFe ratios. The I-ratios of peridotitic and eclogitic garnets are offset from each other at low Fe 3+ /ΣFe ratios (≤0.13); I-ratio garnet XANES calibrations are composition-specific. The E 0.9 feature is independent of garnet major element composition in spectra that have been corrected for the effects of self-absorption. We produce two Fe Kα garnet XANES calibrations based on variations in the E0.9 feature; one calibration with all garnet reference materials included (Fe 3+ /ΣFe up to 1.0; “all garnet calibration”) and another calibration specific to garnets with low Fe 3+ /ΣFe ratios (“low ferric calibration”). Fe 3+ /ΣFe ratios calculated from the mean of up to 25 flank method measurements on eight garnet reference materials fall within 4% absolute of a one-to-one correlation with Fe 3+ /ΣFe ratios measured by Mössbauer. The standard error of the mean Fe 3+ /ΣFe ratio calculated from flank method approaches the Mössbauer-determined Fe 3+ /ΣFe ratio within estimated error (3%) after three analyses. Flank method precision is enhanced at higher beam current; however, the precision of the flank method does not approach the precision of XANES under any microprobe analytical condition tested here. Garnet reference materials detailed here are available by request to the Smithsonian Institution.

58 GEOSCIENCES↗