Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “rock”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Quantifying Dynamic Water Storage in Unsaturated Bedrock with Borehole Nuclear Magnetic Resonance

Quantifying the volume of water that is stored in the subsurface is critical to studies of water availability to ecosystems, slope stability, and water-rock interactions. In a variety of settings, water is stored in fractured and weathered bedrock as rock moisture. However, few techniques are available to measure rock moisture in unsaturated rock, making direct estimates of water storage dynamics difficult to obtain. Here, we use borehole nuclear magnetic resonance (NMR) at two sites in seasonally dry California to quantify dynamic rock moisture storage. Additionally, we show strong agreement between NMR estimates of dynamic storage and estimates derived from neutron logging and mass balance techniques. The depths of dynamic storage are up to 9 m and likely reflect the depth extent of root water uptake. To our knowledge, these data are the first to quantify the volume and depths of dynamic water storage in the bedrock vadose zone via NMR.

58 GEOSCIENCES↗
DOE OSTI2026

Stress and Strain Heterogeneity and Persistence in Uniaxially‐ and Triaxially‐Loaded Sandstone

Two critical questions in brittle rock mechanics are how rocks developed localized strains and to what extent internal stress heterogeneity controls this localization and subsequent macroscopic failure. Definitive answers have not yet been found, but would provide insight into rock fracture mechanics as relevant to hydrocarbon extraction and sequestration. Here, we use synchrotron X‐ray tomography (XRT) and 3D X‐ray diffraction (3DXRD) during uniaxial and triaxial tests on Nugget and Bentheimer sandstones to examine strain and stress localization prior to mechanical failure. 3DXRD was used to measure intra‐granular lattice strains which were used to compute elastic stress tensors of each grain. Digital volume correlation (DVC) was applied to XRT images to determine the strain field in the sample. Both samples featured marked spatial heterogeneity, localization, and temporal persistence of elevated stresses and strains during their mechanical deformation toward failure. Both samples featured a majority of grains with at least one principal stress component that was tensile, a signature of the influence of heterogeneity on stress transmission. Measurements further revealed that compressive stress orientations and statistics evolved in a similar manner to those of inter‐particle forces in loose granular materials, with triaxially‐compressed rock exhibiting enhanced grain stress heterogeneity compared to uniaxially‐compressed rock. Our results complement recent work by others who employed XRT and scanning 3DXRD to study triaxially‐compressed sandstone, but extend those results to uniaxial compression, sandstones of varied porosity, and grain stress measurements throughout the 3D full extent of the samples rather than in a single layer examined with scanning 3DXRD.

58 GEOSCIENCES↗
DOE OSTI2026

A wide-range thermodynamically consistent constitutive model for dynamic loading of partially saturated porous media

SUMMARY This paper presents a novel thermodynamically consistent constitutive model for partially saturated porous rocks across a wide range of conditions. The material states generated behind the shock wave from an explosive source can vary significantly, ranging from crushed and melted rock near the source to a poroelastic medium in the far field. In the model, rock strength is determined by the effective pressure, which is calculated using two independent equations of state: one for the solid rock and another for the pore fluid. The model accounts for shock-induced liquefaction resulting from fluid pressure buildup in the pore spaces near the explosive source. Simultaneously, it describes the increase in wave propagation speed due to elastic pore contraction in both dry and partially saturated rocks. This model is applied to investigate how fluid saturation affects the amplitude and shape of the generated waves, as well as the residual stress surrounding the cavity formed by spherical explosions.

Vorobiev, Oleg [Lawrence Livermore National Labora↗
DOE OSTI2021

Determining fast-S and slow-S propagation directions with SV-P data produced by buried explosives and recorded with vertical geophones

We have evaluated the concept of practicing S-wave reflection seismology with legacy 3D seismic data generated by a P-wave source and recorded with only vertical geophones. This type of S-wave imaging is based on the principle that seismic P-wave sources not only produce a downgoing illuminating P wavefield, but they also simultaneously produce a downgoing illuminating SV wavefield that, in almost all cases, is suitable for S-wave reflection imaging. The S-mode used in this study is the SV-P, or converted-P, mode. This mode involves a downgoing illuminating SV wavefield and an upgoing reflected P-mode that is recorded by vertical geophones. In flat-layered stratigraphy, the lengths of the SV and P raypaths in SV-P imaging are identical to the lengths of the SV and P raypaths in P-SV imaging with P-sources and 3C geophones. P-SV imaging of deep rocks has been practiced for more than two decades; SV-P imaging is a new concept. SV-P data should provide the same options for investigating deep rocks as do P-SV data. We have determined one of the equivalences between SV-P data extracted from vertical-geophone data and P-SV data extracted from horizontal geophones: that both modes react to azimuth-dependent variations in the S velocity in anisotropic rocks. Azimuthal variations in the SV-P traveltime can be used to define the polarization direction of the fast-S-wave mode, which is also the azimuth of the maximum horizontal stress (SHmax). Our investigation demonstrates a noninvasive method for monitoring changes in the SHmax azimuth across a CO 2 storage reservoir, or any targeted porous rock, as fluids are cycled into, and then out of, that rock’s pore space.

Geochemistry & Geophysics↗
DOE OSTI2020

Crystalline Disposal R&D at LBNL: FY20 Progress Report

Within the Spent Fuel & Waste Science and Technology (SFWST) Program, research work continues further with the goal to better understand long-term performance of disposal systems in three main geologic rock types: clay/shale, salt, and crystalline rock. This report documents LBNL’s research activities related to investigations of crystalline host rock according to the scope of two work packages: SF-20LB01030207 “Crystalline International Collaborations – LBNL,” and SF-20LB01030203 “Crystalline Disposal R&D – LBNL.” These research activities correspond are related to key Features, Events, and Processes (FEPs), ranked in importance from medium to high, as listed in Table 7 of the Used Fuel Disposition Campaign Disposal Research and Development Roadmap (FCR&D-USED-2011-000065 REV0) (Nutt, 2011). Specifically, these research activities address FEP 2.2.01, Excavation Disturbed Zone (EZD). The results of these research activities provide important insights into understanding and predicting flow and transport processes that could occur in low-permeability crystalline rocks, in which fractures might serve as main conduits for fluid flow and radionuclide transport. The evolution of the EDZ during the excavation of the tunnel as well as the evolution of microcrack growth within the EDZ after the emplacement of backfill are critical for predicting the long- term behavior of the EDZ. A number of factors including stress, temperature, water activity, capillary pressure, chemistry, and mineralogy can affect the rock evolution, which require advanced experimental tools to study it.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

KOMPASS: Compaction of crushed salt for the safe containment

In Germany, rock salt formations are a possible host rock taken into account for the safe disposal of heat-emitting radioactive waste. With respect to crushed salt will be used in the repository for backfilling of open cavitied (using dry material). With time, the crushed salt will be compacted by the convergence of the host rock and reaches porosities comparable with the rock salts. The compaction behaviour of crushed salt has been investigated within the last 40 years, however, its behaviour at low porosities and the resulting low permeabilities becomes relevant with the introduction of the approach of the containment providing rock zone. In the current state, the database and process understanding have some important gaps in knowledge referring the material behaviour, existing laboratory and numerical models, especially for the porosity range. The objective of this project was the development of methods and strategies for the reduction of deficits in the prediction of crushed salt compaction leading to an improvement of the prognosis quality. It includes the development of experimental methods for determining crushed salt properties in the range of low porosities, the enhancement of process understanding and the investigation and development of existing numerical models.

Crushed salt compaction↗
DOE OSTI2022

Tuscaloosa Marine Shale Laboratory

The Tuscaloosa Marine Shale (TMS) in Louisiana and Mississippi is an Upper Cretaceous source rock formation sandwiched between the sands of the upper and lower Tuscaloosa sections. The TMS is believed to be the source rock for underlying prolific Tuscaloosa sand formation. The TMS has an unproven estimate of 7,000,000,000 bbls of recoverable oil while its current total average production is about 3,000 bbls of oil per day in 2017. In 2013 and 2014, more than 80 wells were drilled horizontally into the TMS that were fractured using multi-stage fracturing technology. The results from this have been mixed, but recent production for several wells show an appealing initial oil production rate of more than 1000 bbl/day. The preliminary core analysis by industry partners and a few literature studies shows that the TMS is one of the most clay-rich and sensitive shales to water. Due to these and other technical problems, there is high risk for the economic development of TMS compared to other shale plays. The experiences of major industrial players in the TMS show the necessity of open and collaborative efforts to better understand the critical gaps in the development of this challenging and potentially highly economic shale play to enable more cost-efficient and environmentally-sound recovery from this unconventional liquid-rich shale play. The overall objective of this project is to form a consortium of science and industry partners to address the following six major objectives using scientific and technical approaches: 1. To improve wellbore integrity by better understanding the sources of the wellbore instability issues, proposing innovative mud and cement design for the TMS. 2. To improve formation evaluation using laboratory techniques for the evaluation of mineralogical composition, organic content, and produced-water chemistry as well as well log and geophysical analysis. 3. To determine the role of geologic discontinuities on fracture growth and shale creep behavior using digital image correlation technique. 4. To investigate the application of stable CO 2 foam and super-hydrophobic proppants for improved reservoir stimulation. 5. To better understand the nature of water/hydrocarbon/CO 2 flow in clay and organic-rich formation and the role of water/fluid interaction on recovery. 6. To prepare better socio-economic environment for TMS development by community engagement. Subsequently, the TMS virtual laboratory conducted testing and analysis of various properties of rock and formation fluids from the TMS, including but not limited to the following: Analyzing reports and logs to better understand the source of wellbore instability in TMS wells; Experiments to design a customized cement based on TMS requirements; Experiments to obtain the mineralogical and geochemical composition of TMS samples; Seismic analysis of TMS geophysical data to better predict total organic carbon (TOC) content and brittleness in TMS; Well log analysis to better estimate the TOC and geo-mechanical properties of TMS; Experiments on formation water to understand the chemistry of produced water; Experiments to determine the role of lamination and natural fractures on fracture propagation or rock deformation using digital image correlation technique in in-direct tensile tests, semi-circular bend test and creep tests Experiments to determine the stability and rheological properties of nanoparticle-stabilized CO 2 foam in TMS rock samples; Experiments to determine fluid dynamics in un-propped TMS fractures and the role of nano-coating of proppants on fluid dynamics in fractures with proppants; Micro-fluidics experiments to enhance the understanding of fluid dynamics in tight liquidrich pores with high clay content; Socio-economic studies to better engage communities in TMS development.

58 GEOSCIENCES↗
DOE OSTI2024

DECOVALEX-2023: Task C Final Report

The Full-scale Emplacement (FE) heater experiment at the Mont Terri Underground Rock Laboratory (URL) was designed and conducted by Nagra to replicate an emplacement tunnel of Nagra’s reference repository design at 1:1 scale. Alongside testing the technical feasibility of constructing disposal tunnels, emplacing waste containers in the tunnels and then backfilling them, the main goals of the FE experiment are (1) to obtain a better understanding of the coupled effects of induced thermo-hydro-mechanical (THM) processes that may occur and (2) to validate existing coupled THM models (Müller et al., 2017). A key aspect of ensuring safety for repositories located in low-permeability rock involves minimizing any damage to the rock itself, thereby preserving its integrity and promoting a stable environment Amongst a number of processes that could damage the rock is the increase in pore pressure due to thermal loading caused by heat emitted from the waste. To reduce the potential damage of the rock, it is important to analyse the evolution of heat over time due to the heat load of the containers and assess possible consequences by coupled THM models. The aim of Task C of DECOVALEX-2023 was to build 3D numerical models of the FE experiment, focussing in particular on the heating induced pore pressure change in the Opalinus Clay. Data from a large number of sensors were available from the FE experiment for model comparison. These sensors measured temperature and relative humidity in the bentonite around the heaters, and temperature, pressure and displacement/strain in the surrounding Opalinus clay. Data were available from the start of excavation (April 2012) up to August 2020 for most sensors (more than 5 years from the start of heating in December 2014). To fulfil the overall aim of the task, the work was broken down into a number of steps, starting with simpler models to build confidence in each team’s approach and then moving to more complex models that better represent the FE experiment. Step 0 consisted of 2D benchmark models, gradually increasing the number of processes that are represented from thermal (T) only models in Step 0a, to coupled thermal hydraulic (TH) models in Step 0b with a representation of changing porosity, to coupled thermo-hydro-mechanical (THM) models in Step 0c, where porosity changes are calculated by the mechanical model. A detailed specification of processes, parameters, initial and boundary conditions was provided for this step, with the ambition that all teams would work towards close agreement in their model results, thus building confidence in the model implementations. vi It was not straightforward to achieve agreement between the teams, so additional steps (Step 0b2, 0b3, 0c2, 0c3) were added along with derivation of some analytical solutions against which the models could be compared. The reasons for the differences between teams were investigated and found to be caused primarily by different conceptual model assumptions (including temperature dependence of the thermal expansion of water), different model formulations (including porosity evolution) and differences in modelled domain sizes, boundary conditions and grid discretisation. This demonstrates that comparisons between multiple modelling teams and/or comparison with analytical results and experimental data are highly beneficial in providing an indication of uncertainty in model predictions. At the conclusion of Step 0, almost all teams had achieved a close agreement in model results and those that had not achieved an agreement knew the reason for this. Step 1 moved from 2D models to 3D models of the FE experiment without adding technical features like shotcrete or EDZ, and only considering the heating phase. Initially the 3D model was tightly specified to continue to build confidence in the model implementations (Step 1a). The results of Step 1a were compared to the data from the FE-experiment without the teams seeing the data. The teams were then provided with a sub-set of the data from the FE-experiment and invited to consider how best to use the large dataset for model comparison (Step 1b). Teams were then asked to use the data provided to calibrate their models, only changing material property values rather than adding features or processes to their models (Step 1c). In Step 1, teams were asked to only model the heating phase of the experiment, so pressure in the Opalinus Clay was reported as change in pressure since the initial conditions were specified rather than modelled. The change from 2D to 3D models was accompanied by an increase in the dispersion of results between the teams. Some of this was resolved during the task, but some remained and is potentially due to model discretisation. Calibration of parameters was useful in improving the fit of the models to the data but the remaining differences indicated that the models were missing features or processes. In Step 2, the teams were asked to update their models with additional features and processes as well as calibrating parameters to try and improve the fit of the models to the data. Teams were encouraged to represent ventilation of the open FE tunnel prior to backfilling with heaters and bentonite and in Step 2, the absolute pressure in the Opalinus Clay was compared between the teams. Teams took different approaches, but there was consideration of adding shotcrete and an EDZ into the model, representing stress change during excavation and different approaches to modelling ventilation of the FE tunnel. Overall, the documented results showed a very good agreement for temperature. The results for porewater pressure evolution showed a significant improvement for most teams compared to Step 1c with a good agreement to the measurements for several teams whereas some teams overpredicted the pressure increase and others overpredicted the drainage effect especially for the sensors close to the heater. Step 3 was an opportunity for teams to use the models developed in Step 1 and Step 2 to make predictions about the temperature and pressure changes that will be expected at the FE experiment over the next few years in light of the planned changes in thermal output of the heaters.

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

Analyzing the Validity of Brazilian Testing Using Digital Image Correlation and Numerical Simulation Techniques

Characterizing the mechanical behavior of rocks plays a crucial role to optimize the fracturing process in unconventional reservoirs. However, due to the intrinsic anisotropy and heterogeneity in unconventional resources, fracture process prediction remains the most significant challenge for sustainable and economic hydrocarbon production. During the deformation tracking under compression, deploying conventional methods (strain gauge, extensometer, etc.) is insufficient to measure the deformation since the physical attachment of the device is restricted to the size of the sample, monitoring limited point-wise deformation, producing difficulties in data retrieval, and a tendency to lose track in failure points, etc. Where conventional methods are limited, the application of digital image correlation (DIC) provides detailed and additional information of strain evolution and fracture patterns under loading. DIC is an image-based optical method that records an object with a camera and monitors the random contrast speckle pattern painted on the facing surface of the specimen. To overcome the existing limitations, this paper presents numerical modeling of Brazilian disc tests under quasi-static conditions to understand the full-field deformation behaviors and finally, it is validated by DIC. As the direct tensile test has limitations in sample preparation and test execution, the Brazilian testing principle is commonly used to evaluate indirectly the tensile strength of rocks. The two-dimensional numerical model was built to predict the stress distribution and full-field deformation on Brazilian disc under compression based on the assumptions of a homogenous, isotropic and linear elastic material. The uniaxial compression test was conducted using the DIC technique to determine the elastic properties of Spider Berea sandstone, which were used as inputs for the simulation model. The model was verified by the analytical solution and compared with the digital image correlation. The numerical simulation results showed that the solutions matched reasonably with the analytical solutions where the maximum deviation of stress distribution was obtained as 14.59%. The strain evolution (normal and shear strains) and displacements along the central horizontal and vertical planes were investigated in three distinguishable percentages of peak loads (20%, 40%, and 90%) to understand the deformation behaviors in rock. The simulation results demonstrated that the strain evolution contours consistently matched with DIC generated contours with a reasonable agreement. The changes in displacement along the central horizontal and vertical planes showed that numerical simulation and DIC generated experimental results were repeatable and matched closely. In terms of validation, Brazilian testing to measure the indirect tensile strength of rocks is still an issue of debate. The numerical model of fracture propagation supported by digital image correlation from this study can be used to explain the fracturing process in the homogeneous material and can be extended to non-homogeneous cases by incorporating heterogeneity, which is essential for rock mechanics field applications.

58 GEOSCIENCES↗
DOE OSTI2021

Geology and Hydrothermal Alteration of the Raft River Geothermal System, Idaho

The Raft River geothermal system is located in southern Idaho, near the Utah-Idaho state boarder in the Raft River Valley. The field, which is owned and operated by U.S. Geothermal, has been selected as an EGS demonstration site by the U. S. Department of Energy. This paper summarizes ongoing geologic and petrologic investigations being conducted in support of this project. The reservoir is developed in fractured Proterozoic schist and quartzite, and Archean quartz monzonite cut by younger diabase intrusions. The basement complex was deformed during the mid Tertiary and covered by approximately 5000 ft of late Tertiary sedimentary and volcanic deposits. Listric normal faults of Cenozoic age disrupt the Tertiary deposits but do not offset the basement rocks. RRG-9, the target well, was drilled southwest of the main well field to a measured depth (MD) of 6089 ft. The well is deviated to the west and cased to a depth of 2316 ft MD. It penetrated the Proterozoic reservoir rocks at a depth of 5286 ft MD and encountered a maximum temperature of 139 degrees C. Bottomhole temperatures in other deep wells range from 133 to 149 degrees C. X-ray diffraction and thin section analyses are being conducted on samples from 5 deep wells, RRG- 1, 2, 3, 7 and 9, to determine the characteristics of the rock types and hydrothermal alteration within the geothermal system. Thin section analyses of samples from RRG-9 document the presence of strong alteration and brecciation at the contact between the Tertiary and basement complex. The Tertiary rocks consist of ash-flow tuffs, lava flows, tuffaceous siltstone, greywacke, and sandstone. No core is available from RRG-9 but core was obtained from RRG-3C. The sample is a brecciated and altered siltstone from the base of the Tertiary sequence and is similar to rocks at the base of the Tertiary deposits in RRG-9. Here, the results of thermal and quasi-static mechanical property measurements that were conducted on the core sample are presented.

Jones, Clay↗
DOE OSTI2024

Evaluating Heat Extraction Performance of Closed-Loop Geothermal Systems with Thermally Conductive Enhancements in Conduction-Only Reservoirs: Preprint

We investigated the impact of disc-shaped and linear thermal enhancements on performance of closed-loop geothermal systems using steady-state and transient COMSOL numerical simulations. Thermal enhancements refer to thermally conductive material introduced in the rocks surrounding a wellbore to compensate for the relatively low rock thermal conductivity and increase heat extraction. Materials proposed include composite cements with thermally conductive additives such as metals or graphite to obtain overall cement thermal conductivities of 10 to 100x the rock thermal conductivity. Our simulation approach was modeling a subsegment (e.g., 100 m long) of a wellbore with idealized (e.g., perfect disc-shaped) thermal enhancements, and we report results as relative increase in heat extraction with respect to systems without thermal enhancements. We did not evaluate technical feasibility or cost of installing thermal enhancements. Simulations indicate that for 5-mm thick disc-shaped thermal enhancements with 5-m radius, repeated every 1 m along the wellbore and with thermal conductivity of 100x the rock thermal conductivity, the thermal output increases roughly 20% with respect to a closed-loop system without thermal enhancements. For fishbone structure thermal enhancements, we estimate a roughly 7% increase in thermal output for 5-m long radially outward pipes with half the radius of the main wellbore, repeated every 2.5 m, with thermal conductivity of 100x the rock thermal conductivity.

advanced geothermal systems↗
DOE OSTI2021

Acoustic velocity and permeability of acidized and propped fractures in shale

We report that from geochemical reactions to proppant emplacement, hydraulic fracturing induces various chemomechanical fracture alterations in shale reservoirs. Hydraulic fracturing through the injection of a vast amount and variety of fluids and proppants has substantial impacts on fluid flow and hydrocarbon production. There is a strong need to improve our understanding on how fracture alterations affect flow pathways within the stimulated rock volume and develop monitoring tools. We have conducted time-lapse rock-physics experiments on clay-rich (carbonate-poor) Marcellus shales to characterize the acoustic velocity and permeability responses to fracture acidizing and propping. Acoustic P- and S-wave velocities and fracture permeability were measured before and after laboratory-induced fracture alterations along with microstructural imaging through X-ray computed tomography and scanning electron microscopy. Our experiments indicate that the S-wave velocity is an important geophysical observable, particularly the S-wave polarized perpendicular to fractures because it is sensitive to fracture stiffness. The acidizing and propping of a fracture decrease its elastic stiffness. This effect is stronger for acidizing, so it is possible that proppant monitoring will be masked by chemical alteration except when propping is highly efficient (i.e., most fractures are propped). However, fracture permeability is undermined by the softening of fracture surfaces due to acidizing, while being greatly enhanced by propping. These contrasting effects on fluid flow in combination with similar seismic attributes indicate the importance of experiments to improve existing rock-physics models, which must include changes to the rock frame. Such improvements are necessary for a correct interpretation of seismic velocity monitoring of flow pathways in stimulated shales.

58 GEOSCIENCES↗
DOE OSTI2020

The German Site Selection Process - Implementation and Challenges - 20528

Federal Company for Radioactive Waste Disposal (Bundesgesellschaft fuer Endlagerung mbH; hereinafter referred to as BGE) is a state-owned company established in September 2016 following the 'Act on the Organisational Restructuring in the Field of Radioactive Waste Storage'. 'Repository Site Selection Act' (hereinafter referred to as StandAG) defines the site selection procedure, aiming to ensure the best possible safety for storing high-level radioactive waste in Germany for the duration of one million years. The selection procedure is a participatory, transparent, learning and self-questioning process based on scientific expertise. Possible repository sites are sought all over Germany, without excluding or predetermining any potential site at the beginning. The procedure consists of three phases. After each phase, the legislative decides on the regions for further consideration in the next phase of the site selection procedure. The current first phase (Phase I) is divided into two steps. Objective of Step 1 is to determine subareas of favorable geological conditions based on available data. Step 2 of Phase I will lead to identification of specific regions that are considered for further surface-based geoscientific and geophysical exploration. Potential host rocks are rock salt, clay rock and crystalline rock. As part of the site selection procedure, an overview of the status of the planning and implementation of R and D activities is given. Topics of particular interest are, for example, the adaptation of the maximum temperature of 100 deg. C at the container surface, the impact of long-term climate variations as well as minimally invasive exploration methods. (authors)

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

Dating fractures using luminescence

Rock fracturing (cracking) is a universal process that drives and limits chemical degradation, sediment production and erosion, and deterioration of infrastructure. Despite extensive research gains in rock mechanics on one hand and geochronology on the other, there remains a glaring gap in our ability to understand the long term evolution of natural, in situ fractures. Here we develop a novel fracture exposure dating technique, grounded in modern advances in luminescence geochronology. We apply our new dating method to a granitic boulder from a glacial outwash terrace in California, US. We conclude that the longest, clast-splitting E-W fracture appeared shortly after the boulder's deposit, whereas the secondary N-S fracture appeared 5 ka after the deposition, approximately correlating with the Last Glacial Maximum and Younger Dryas periods of the region, respectively. However, dating of the third fracture (<< 50 µm width) which does not fully split the rock, is ambiguous due to negligible daylight penetration and poor determination of fracture width. The fracture dating method presented herein brings with it the potential to decipher relationships that are crucial for the interpretation and modeling of, for example, long-term landscape and atmospheric evolution relating rock weathering to climate change and erosion.

58 GEOSCIENCES↗
DOE OSTI2020

Localization of deformation in fluid-saturated sandstone

Herein we discuss a laboratory study of inelastic deformation and localization was conducted to evaluate microseismicity and fracture in fluid-saturated rock. Plane strain compression experiments were performed on oil- and water-saturated Berea sandstone under dry, drained, undrained, and quasi-unjacketed conditions. Parameters associated with the elastic and inelastic material response were accurately measured and used in a constitutive model to predict localization of deformation under the drained and undrained conditions. Acoustic emission (AE) activity was recorded and the onset of inelastic response in fluid-saturated rock, similar to dry rock, coincided with an increase in AE rate. Fracture initiation and propagation were successfully monitored by AE locations. Clustering of AE events and the start of inhomogeneous deformation inferred from displacement measurements are shown to be well correlated with the model prediction of the onset of localized deformation under various rock-fluid interaction conditions.

42 ENGINEERING↗
DOE OSTI2022

Integrated Experimental and Modeling Study of Geochemical Reactions of Simple Fracturing Fluids with Caney Shale

Interactions between rock minerals and hydraulic fluids directly impact the geochemical and geomechanical properties of shale formations. However, the mechanisms of geochemical reactions in shale unconventional reservoirs remain poorly understood. Here, to investigate the geochemical reactions between shale and hydraulic fracturing fluids, a series of batch reactor experiments were undertaken. Three rock samples with different mineralogical compositions and three fluid samples of different compositions (deionized water, deionized water + 2% potassium chloride (KCl), and deionized water + 0.5% choline chloride (C 5 H 14 ClNO) were used. Experiments were undertaken at reservoir temperature and atmospheric pressure. Elemental compositions of effluents after 1, 3, 7, 14, 28 days were analyzed using Inductively Couple Plasma Mass Spectrometry (ICP-MS). Medical Computed Tomography (CT) scan and X-ray Fluorescence (XRF) were conducted on the entire core run to help upscale results obtained from rock-fluid interaction experiments. Geochemical modeling using a reactive simulator, TOUGHREACT, was undertaken to corroborate experimental results. Results show that lower pH triggered high dissolution rates in the rock samples, especially the carbonate components. As pH increased, the rate of dissolution declined significantly, though for most cases dissolution still continued. Observed dissolved silica concentrations were much higher than the quartz solubility, suggesting that much of the silica originates from more soluble silica polymorphs and possibly desorption from clay mineral exchange sites. Concentration of most elemental species in solution increased but aluminium and magnesium concentrations declined rapidly following initial entry into solution. Geochemical modeling corroborated the conclusions regarding mineral dissolution and precipitation observed from experiments, notably; the dissolution of calcite and oxidation of pyrite in reacted shale samples, the likely presence of silica polymorphs such as opal, chalcedony or amorphous silica in these samples, and the reduction of Al and Mg concentrations in solution by precipitation of secondary aluminosilicate phases. The de-flocculation of clay minerals during reaction implies fines migration after hydraulic fracturing. This is detrimental to reservoir productivity as clay fines are displaces and lodged within the micro and nano-fractures created during fracturing. The immediate consumption of aluminium and magnesium also has implications on blockage of hydrocarbon pathways due to precipitation of new minerals in these locations.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2020

Engineered Barrier System R&D and International Collaborations – LANL (FY20): Spent Fuel and Waste Disposition

The United States Department of Energy’s Spent Fuel and Waste Disposition (SFWD) program is investigating the design and safety function of generic nuclear geologic repositories in a variety of geologic settings (salt, argillite, and crystalline rock). Different configurations and loadings of spent nuclear fuel and waste within disposal canisters are also being investigated, some of which have the potential to generate repository temperatures higher than previously considered (i.e., temperatures >100ºC) by foreign and domestic concepts. This report expands on engineered barrier material stability in a high temperature crystalline rock repository through high temperature hydrothermal experiments. Experiments were designed to develop engineered barrier system (EBS) concepts in a hightemperature crystalline environment in 1) bentonite-Grimsel Granodiorite interactions, 2) bentonite-cement reactions, and 3) interaction between waste canister materials and bentonite. Experiment results are applied to understanding long-term repository performance in terms of radionuclide isolation. One hydrothermal experiment was completed in the rocking autoclaves at LANL in FY-20: IEBS-6 (Grimsel Granodiorite + Wyoming bentonite + cured ordinary Portland cement + Grimsel Granodiorite synthetic groundwater, 250ºC/150 bar, 8 weeks). Several other experiments were planned but were delayed due to the COVID-19 pause in laboratory work. In addition, some characterization of the reaction products of experiments conducted in FY-20 was prevented by COVID-19. Mineral phase chemistry was not measured via electron microprobe analyses for IEBS-6. Quantitative X-ray diffraction results from IEBS-6 and HBT-1 were not completed. The missing analyses and discussion of the results will be included in next year’s report. New characterization that was completed includes scanning electron microscopy of reaction products from IEBS-6, quantitative X-ray diffraction results are presented from IEBS-1 through IBES-5, preliminary scanning electron microscope images and chemical analyses for IEBS-6, X-ray diffraction of the clay fraction from all experiments, measurement/imaging of mineral growth on the surface of steel coupons. Major observations pertaining to bentonite stability in a Grimsel Granodiorite environment include the stability of Na-montmorillonite at 250°C, the formation of trace CSH phases, and the formation of bentonite colloids on experiment cooling. The addition of a cured chip of Portland cement to the bentonite-Grimsel system results in slightly higher pH values and the formation of diverse secondary mineral phases that were not observed in the previous experiments (e.g., analcime, garronite, CSH phases). The new characterization efforts related to the interaction of stainless-steel coupons and bentonite clay focused on thickness and mineralogy of phases that formed at the steel surface. In Wyoming bentonite + Grimsel Granodiorite systems, newly formed minerals at the bentonitesteel coupon interface included alteration of the outermost steel edge to Fe,Cr-oxide phases, followed by Fe-rich phyllosilicates (Fe-saponite, chlorite) and interbedded Fe,Cr,Ni-sulfide phases (pentlandite). Hydrothermal experiments were completed to assess uranium-carbonate complexation at conditions relevant to high-temperature disposal. Autoclave solubility experiments were conducted at 150 to 250ºC with a range of carbonate and uranium concentrations. The experiment results were characterized via situ UV-Visible spectroscopy and synchrotron-based in situ XAS techniques. Results show a significant decrease in the stability of uranyl-carbonate complexes at temperatures above 100ºC along with the precipitation of uranium oxides. Further, at T > 200ºC, results show that uranyl-hydroxyl complexes control solubility of uranium instead of the previously predicted uranyl-carbonate species. These results are significant for understanding the mobility of uranium in the EBS, which will likely contain carbonate-rich fluids. International research efforts focused on three main areas: 1) participation in international conferences, 2) building collaborations with foreign repository programs, and 3) the initiation of an experimental program to complement the full-scale HotBENT test at the Grimsel test site. This experiment included Wyoming bentonite + low carbon steel + Grimsel Granodiorite synthetic groundwater and was run at the planned maximum temperature of the HotBENT test (200°C). Complete characterization of reaction products was hindered by disruptions to laboratory work but will be reported in the next FY. The experimental results obtained in FY-20 continue to document the wide-ranging effects of bulk composition and pressure-temperature conditions in the mineralogical and geochemical evolution of a high-temperature repository environment. Concepts developed will be used to inform models of long-term material stability in a generic crystalline rock-hosted repository.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Water-Weakening and Time-Dependent Deformation of Organic-Rich Chalks

Abstract The Ghareb Formation is a shallowly buried porous chalk in southern Israel that is being considered as a host rock for a geologic nuclear waste repository. Setup and operation of a repository will induce significant mechanical, hydrological and chemical perturbations in the Ghareb. Developing a secure repository requires careful characterization of the rock behavior to different loads. To characterize hydromechanical behavior of the Ghareb, several short- and long-term deformation experiments were conducted. Hydrostatic loading tests were conducted both dry and water-saturated, using different setups to measure elastic properties, time-dependent behavior, and permeability. A set of triaxial tests were conducted to measure the elastic properties and rock strength under differential loading at dry and water-saturated conditions. The hydrostatic tests showed the Ghareb began to deform inelastically around 12–15 MPa, a relatively low effective pressure. Long-term permeability measurements demonstrated that permeability declined with increasing effective pressure and was permanently reduced by ~ 1 order of magnitude after unloading pressure. Triaxial tests showed that water saturation significantly degrades the rock properties of the Ghareb, indicating water-weakening is a significant risk during repository operation. Time-dependent deformation is observed during hold periods of both the hydrostatic and triaxial tests, with deformation being primarily visco-plastic. The rate of deformation and permeability loss is strongly controlled by the effective pressure as well. Additionally, during holds of both hydrostatic and triaxial tests, it is observed that when water-saturated, radial strain surpassed axial strain when above effective pressures of 13–20 MPa. Thus, deformation anisotropy may occur in situ during operations even if the stress conditions are hydrostatic when above this pressure range.

58 GEOSCIENCES↗