Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “rock deformation”

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

Experimental Investigation of Hydraulic Fracturing and Stress Sensitivity of Fracture Permeability Under Changing Polyaxial Stress Conditions

Understanding and predicting fracture propagation and subsequent fluid flow characteristics is critical to geoenergy technologies that engineer and/or utilize favorable geological conditions to store or extract fluids from the subsurface. Fracture permeability decreases nonlinearly with increasing normal stress, but the relationship between shear displacement and fracture permeability is less well understood. We utilize the new Geo–Reservoir Experimental Analogue Technology (GREAT cell), which can apply polyaxial stress states and realistic reservoir temperatures and pressures to cylindrical samples and has the unique capability to alter both the magnitude and orientation of the radial stress field by increments of 11.25° during an experiment. We load synthetic analogue materials and real rock samples to stress conditions representative of 500–1,000 m depth, investigate the hydraulic stimulation process, and then conduct flow experiments while changing the fluid pressure and the orientation of the intermediate and minimum principal stresses. High–resolution circumferential strain measurements combined with fluid pressure data indicate fracture propagation can be both stable (no fluid pressure drop) and unstable (fluid pressure drop). The induced fractures exhibit both opening and shear displacements during their creation and/or during fluid flow with changing radial stress states. Flow tests during radial stress field rotation reveal that fracture normal effective stress has first–order control on fracture permeability but increasing fracture offset can lead to elevated permeabilities at maximum shear stress. Here, the results have implications for our conceptual understanding of fracture propagation as well as fluid flow and deformation around fractures.

58 GEOSCIENCES↗
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

Finite Element Analysis (FEA) for Water-Foam Fracturing of Granite Rock

In addition to the foam data that were obtained from literature and that were collected from the current study, simulation data was also generated from finite element analysis (FEA) conducted in this study using COMSOL Multiphysics software. The FEA models were built to simulate the experiments conducted at Oak Ridge National Laboratory (ORNL) on cement and granite samples. In these FEA models, temperature was kept at ambient while the pressure profile resembled the loading conditions during the ORNL experiments, where pressure was either monotonically increased or applied cyclically. The cement material was used as a model material and was used to study Von Mises stress and tensile stress distribution for different bore hole length geometry using a parametric sweep with water as fracturing fluid using solid-fluid interaction module. For the granite material, FEA models were developed for stress analysis of cylindrical samples with water or foam fluids. The solid mechanics module in COMSOL was implemented to solve for Von Mises stress and tensile stress. The fluid-structure interaction module was implemented to solve for water-foam interaction on granite cylinder with addition of fluid-loading on structure, i.e., large deformation in solid mechanics with no impact on fluid deformation. Foam was considered as a pseudo single-phase compressible fluid for which material properties were calculated from water and gas (nitrogen) phases. The density of foam is calculated as a function of the densities of water and nitrogen, while viscosity is a function of temperature. Four types of FEA analyses were modelled: 1. Monotonic injection with water 2. Monotonic injection with foam 3. Cyclic injection with water 4. Cyclic injection with foam All the COMSOL files are converted to a zip file which is save in .mph.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Proppant embedment in coal and shale: Impacts of stress hardening and sorption

During methane production in CBM reservoirs, the influence of proppant embedment and permeability damage cannot be neglected – especially where the wall-rock is soft. Effective stresses are elevated during methane recovery, increasing both normal loading stress and confinement and simultaneously overprinting sorption-induced volumetric strains. Experiments and analytic modeling are conducted to define key mechanisms controlling these competitive effects. We independently measure overall sample compaction (external LVDT) and local strain (strain gauge) in the matrix to deconvolve proppant embedment in a propped fracture for different conditions of confining stress. The results show symptomatic behaviors of elastic (shale) and elastoplastic (coal) responses of embedment. Different from shale, the evolution of embedment is convex upwards with increased stress where indented depth increases more rapidly as loading stress increases under constant confinement. In addition, a stress-hardening effect is found to play a pivotal role in determining the characteristics of indentation, which are examined in terms of evolution profiles, deformation regimes, embedment slopes, curvatures, yield points and irreversible indentations. Based on the experimental observations a semianalytical model predicts indentation and the evolution of propped permeability under recreated in-situ stress conditions. A simplified case study is conducted to further illustrate the evolution of aperture and permeability of a propped fracture in CBM reservoirs. The modeling results suggest that proppant embedment is significantly overestimated if the variable stress-hardening (VSH) effect is neglected, especially when effective stress is large. Moreover, a decrease in indentation depth possibly occurs during late stage methane production, resulting in a reversal/recovery in fracture closure. This is because desorption-induced shrinkage becomes the predominant effect, causing an increase in aperture and a reduction in the indented volume of proppant. The resulting recovery in permeability implies that the propped coal fracture has the potential to optimally facilitate methane production as a pathway, even at high closure stresses generated by methane drainage.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

Experimental data for damage mechanics simulation challenge

While there are many computational approaches for simulating damage in rock and other materials, few have been ground truth tested with either known experimental data or with blind data sets. Here, in this work, we present a bench-mark laboratory data set for a damage mechanics challenge to compare computational approaches on damage evolution in brittle-ductile materials. The samples were fabricated through additive manufacturing to produce repeatable specimens designed to fail in controlled ways. The failure was induced in the samples using a 3-point bending test to produce different Modes such as Mode I and mixed Modes including I-II, I-III and I-II-III Modes to generate a calibration data set and a blind challenge data set. Data collected included spatial and temporal measurements from traditional digital load–displacement sensors, 2D digital image correlation measurement to map surface deformations, 3D X-ray microscopy to ground-truth the crack-failure geometry, and laser profilometry to capture surface roughness. The data sets are available, on a data repository, to the community to advance computational models to improve our ability to predict damage in brittle-ductile materials.

3-point bending↗
DOE OSTI2020

Fracture, Dissolution, and Cementation Events in Ordovician Carbonate Reservoirs, Tarim Basin, NW China

Ordovician carbonate rocks of the Yijianfang Formation in the Tabei Uplift, Tarim Basin, contain deeply buried (>6000 m), highly productive oil and gas reservoirs associated with large cavities (>10 m). Previous workers inferred that large cavities are paleocaves (paleokarst) formed near the surface and subsequently buried. Alternately, caves may have formed by dissolution at depth along faults. Using 227 samples from 16 cores, we document textures and cement compositions bearing on cavity histories with petrographic, high-resolution scanning electron microscopy (SEM), isotopic, and fluid inclusion microthermometric observations. Results show that dissolution occurred at depth and was caused by (1) acidic fluids derived from Middle-Late Silurian and/or Devonian-Permian hydrocarbon generation and maturation, (2) high-temperature fluids, of which some were associated with Late Permian igneous activity, and (3) Mg-rich fluids that accompanied Jurassic-Cretaceous deformation and the formation of partially open fractures and stylobreccias (fault breccias). The relative paragenetic sequence of the structure-related diagenesis suggests seven stages of fracturing, dissolution, and cementation. Mottle fabrics in the Yijianfang Formation contain argillaceous carbonate-rich silt and are bioturbation features formed within the marine environment. Those mottled fabrics differ from clearly karstic features in the overlying Lianglitage Formation, which formed by near-surface dissolution and subsequent infilling of cavities by allochthonous sediment. Mottle fabrics are crosscut by compacted fractures filled with phreatic-vadose marine cements and followed by subsequent generations of cement-filled fractures and vugs indicating that some fractures and vugs became cement filled prior to later dissolution events. Calcite cements in fractures and vugs show progressively depleted values of δ18O documenting cement precipitation within the shallow (~220 m), intermediate (~625 m), and deep (~2000 m) diagenetic environments. Deep (mesogenetic) dissolution associated with fractures is therefore the principal source of the high porosity-permeability in the reservoir, consistent with other pieces of evidence for cavities localized near faults.

58 GEOSCIENCES↗
DOE OSTI2019

Initial Simulations of Empty Room Collapse and Reconsolidation at the Waste Isolation Pilot Plant

The Waste Isolation Pilot Plant (WIPP) is a geologic repository for defense-related nuclear waste. If left undisturbed, the virtually impermeable rock salt surrounding the repository will isolate the nuclear waste from the biosphere. If humans accidentally intrude into the repository in the future, then the likelihood of a radionuclide release to the biosphere will depend significantly on the porosity and permeability of the repository itself. Room ceilings and walls at the WIPP tend to collapse over time, causing rubble piles to form on floors of empty rooms. The surrounding rock formation will gradually compact these rubble piles until they eventually become solid salt, but the length of time for a rubble pile to reach a certain porosity and permeability is unknown. This report details the first efforts to build models to predict the porosity and permeability evolution of an empty room as it closes. Conventional geomechanical numerical methods would struggle to model empty room collapse and rubble pile consolidation, so three different meshless methods, the Immersed Isogeometric Analysis Meshfree, Reproducing Kernel Particle Method (RKPM), and the Conformal Reproducing Kernel method, were assessed. First, the meshless methods and the finite element method each simulated gradual room closure, without ceiling or wall collapse. All three methods produced equivalent room closure predictions with comparable computational speed. Second, the Immersed Isogeometric Analysis Meshfree method and RKPM simulated two-dimensional empty room collapse and rubble pile consolidation. Both methods successfully simulated large viscoplastic deformations, fracture, and rubble pile rearrangement to produce qualitatively realistic results. In addition to geomechanical simulations, the flow channels in damaged salt and crushed salt were measured using micro-computed tomography, and input into a computational fluid dynamics simulation to predict the salt's permeability. Although room for improvement exists, the current simulation approaches appear promising.

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

Linking Domain Structure Evolution at a Grain Boundary to Piezoelectric Response via Nano‐Diffraction

Electric-field-induced domain structure switching in a 1 𝜇⁢m -thick Pb 0.99 (Zr 0.45 Ti 0.55 ) 0.98 Nb 0.02 O 3 (Nb-doped PZT) bicrystal film was characterized in situ via nano-focused synchrotron diffraction. The epitaxial film was deposited on a (100) SrTiO 3 bicrystal substrate. The changes in domain structure were mapped within a 5 𝜇⁢m × 5 𝜇⁢m area at and around an in-plane tilt-type (23.6°) grain boundary with 50 nm spatial resolution. Rocking curves collected at each point of the mapped area provided the ability to reconstruct spatially-varying three-dimensional (3D) reciprocal space maps around the grain boundary. The 3D reciprocal space maps reveal how different tilted 𝑎 -type domain variants interact with the grain boundary as a function of increasing electric field. Initially, a subset of 𝑎 -type domain variants with their polarization vectors largely orthogonal to the grain boundary were found in greater abundance within 570 nm of the grain boundary, possibly due to X-ray beam-induced local increases in the electrical conductivity, but after the coercive field was exceeded, reconfiguration of the ferroelastic domains was observed. The spatially-varying reciprocal space maps also facilitated evaluation of the strain field across the mapped area, along with 𝑑 33,𝑓 . Significant spatial heterogeneity of strain and 𝑑 33,𝑓 are observed, especially at the coercive field, which was attributed to maintaining deformation compatibility and correlated ferroelastic switching.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Assessing continuum plasticity postulates with grain stress and local strain measurements in triaxially compressed sand

Critical state and continuum plasticity theories have been used in research and engineering practice in soil and rock mechanics for decades. These theories rely on postulated relationships between material stresses and strains. Some classical postulates include coaxiality between stress and strain rates, stress–dilatancy relationships, and kinematic assumptions in shear bands. Although numerical and experimental data have quantified the strains and grain kinematics in such experiments, little data quantifying grain stresses are available. Here, we report the first-known grain stress and local strain measurements in triaxial compression tests on synthetic quartz sands using synchrotron X-ray tomography and 3D X-ray diffraction. We use these data to examine the micromechanics of shear banding, with a focus on coaxiality, stress-dilatancy, and kinematics within bands. Our results indicate the following: 1) elevated deviatoric stress, strain, and stress ratios in shear bands throughout experiments; 2) coaxial principal compressive stresses and strains throughout samples; 3) significant contraction along shear bands; 4) vanishing volumetric strain but nonvanishing stress fluctuations throughout samples at all stages of deformation. Our results provide some of the first-known in situ stress and strain measurements able to aid in critically evaluating postulates employed in continuum plasticity and strain localization theories for sands.

42 ENGINEERING↗
DOE OSTI2024

Insights into the thermo-hydraulic properties of compacted MX80 bentonite during hydration under elevated temperature

In high-level radioactive waste geological repositories, compacted bentonite undergoes coupled thermo-hydraulic processes due to heat released from a central waste canister and groundwater imbibition from the surrounding host rock. An understanding of these processes is essential for long-term simulations of radionuclide migration and canister corrosion, which requires an understanding of temperature effects on the coupled thermo-hydraulic properties governing these processes. In this study, a tank-scale radial infiltration test was used to investigate water imbibition processes in compacted bentonite under a central heater temperature of 200 °C that simulates high thermal gradients in a repository. Interpretation of this test focuses on evaluation of the liquid water wetting front during hydration and the interpretation of the temperature-dependent transient soil water retention curve (SWRC), thermal conductivity function (TCF), and hydraulic conductivity function (HCF). The SWRC during imbibition follows a temperature-dependent wetting path. Temperature effects on the hydraulic conductivity of bentonite in saturated conditions had the greatest effect on the shape of the HCF, with minimal temperature effects at higher suctions. The transient thermal conductivity data matched well with a new TCF linked with the SWRC shape. Although the bentonite layer was restrained, local deformations during hydration may have affected the shapes of the TCF and HCF.

Engineering↗
DOE OSTI2022

Frictional properties of Opalinus Clay: influence of humidity, normal stress and grain size on frictional stability

SUMMARY The Opalinus Clay (OPA) is a clay-rich formation considered as a potential host rock for radioactive waste repositories and as a caprock for carbon storage in Switzerland. Its very low permeability (10−19 to 10−21 m2) makes it a potential sealing horizon, however the presence of faults that may be activated during the lifetime of a repository project can compromise the long-term hydrological confinement, and lead to mechanical instability. Here, we have performed laboratory experiments to test the effect of relative humidity (RH), grain size (g.s.) and normal stress on rate-and-state frictional properties and stability of fault laboratory analogues corresponding to powders of OPA shaly facies. The sifted host rock powders at different grain size fractions (<63 μm and 63 < g.s. < 125 μm), at room (∼25 per cent) and 100 per cent humidity, were slid in double-direct shear configuration, under different normal stresses (5–70 MPa). We observe that peak friction, μpeak and steady-state friction, μss, depend on water vapour content and applied normal stress. Increasing relative humidity from ∼25 per cent RH (room humidity) to 100 per cent RH causes a decrease of frictional coefficient from 0.41 to 0.35. The analysis of velocity-steps in the light of rate-and-state friction framework shows that the stability parameter (a–b) is always positive (velocity-strengthening), and it increases with increasing sliding velocity and humidity. The dependence of (a–b) on slip rate is lost as normal stress increases, for each humidity condition. By monitoring the variations of the layer thickness during the velocity steps, we observe that dilation (Δh) is directly proportional to the sliding velocity, decreases with normal stress and is unaffected by humidity. Microstructural analysis shows that most of the deformation is accommodated within B-shear zones, and the increase of normal stress (σn) promotes the transition from strain localization and grain size reduction to distributed deformation on a well-developed phyllosilicate network. These results suggest that: (1) the progressive loss of velocity dependence of frictional stability parameter (a–b) at σn > 35 MPa is dictated by a transition from localized to distributed deformation and (2) water vapour content does not affect the deformation mechanisms and dilation, whereas it decreases steady-state friction (μss), and enhances fault stability.

58 GEOSCIENCES↗
DOE OSTI2020

Experimental Evaluation of Geopolymer, Nano-Modified, and Neat Class H Cement by Using Diametrically Compressive Tests

Abstract Cement is the prime barrier material to maintain zonal isolation during wellbore life-cycle. Different types of the cement have different performances on sealability. The primary objective of this study is evaluation of three types of cement systems, neat Class H, nano-synthetic graphite, and rock-based geopolymer cements, on interfacial bond integrity. Diametrical compressive tests are used on casing–cement–formation systems. Cement performances are tested in a conventional configuration (Berea Sandstone). A universal testing system is utilized to analyze the relationship applied load–axial strain and the digital image correlation (DIC) technique is used to generate the strain field maps. Based on the relationship and DIC results, bond strengths in different cement systems are evaluated. Additionally, microscope images provide insight explanations of the bond from micro-view. Based on the experimental results, failures of all cases initiate and propagate from the casing–cement interface. Then, radial cracks develop through the formation. This is because the mechanical interlocking bonds are formed only at the cement–sandstone interface to enhance the bond strength. The comparison of geopolymer and nano-modified cements with the conventional cement (neat Class H) shows that the unconventional cements have a lower maximum load and higher axial strain before failure of the samples. It indicates that both the two kinds of cements are more ductile than conventional cement and can support more deformation before the initiation of cracks.

Energy & Fuels↗
DOE OSTI2019

Viscoplastic Interpretation of Localized Compaction Creep in Porous Rock

Recent laboratory evidence shows that compaction creep in porous rocks may develop through stages of acceleration, especially if the material is susceptible to strain localization. This paper provides a mechanical interpretation of compaction creep based on viscoplasticity and nonlinear dynamics. For this purpose, a constitutive operator describing the evolution of compaction creep is defined to evaluate the spontaneous accumulation of pore collapse within an active compaction band. This strategy enables the determination of eigenvalues associated with the stability of the response which is able to differentiate decelerating from accelerating strain. This mathematical formalism was linked to a constitutive law able to simulate compaction localization. Material point simulations were then used to identify the region of the stress space where unstable compaction creep is expected, showing that accelerating strains correspond to pulses of inelastic strain rate. Such pulses were also found in full-field numerical analyses of delayed compaction, revealing that they correspond to stages of inception and propagation of new bands across the volume of the simulated sample. These results illustrate the intimate relation between the spatial patterns of compaction and their temporal dynamics, showing that while homogeneous compaction develops with decaying rates of accumulation, localized compaction occurs through stages of accelerating deformation caused by the loss of strength taking place during the formation of a band. In addition, they provide a predictive modeling framework to simulate and explain the spatio-temporal dynamics of compaction in porous sedimentary formations.

58 GEOSCIENCES↗
DOE OSTI2020

Acid Erosion of Carbonate Fractures and Accessibility of Arsenic-Bearing Minerals: In Operando Synchrotron-Based Microfluidic Experiment

Underground flows of acidic fluids through fractured rock can create new porosity and increase accessibility to hazardous trace elements such as arsenic. In this study, we developed a custom microfluidic cell for an in operando synchrotron experiment using X-ray attenuation. The experiment mimics reactive fracture flow by passing an acidic fluid over a surface of mineralogically heterogeneous rock from the Eagle Ford shale. Over 48 h, calcite was preferentially dissolved, forming an altered layer 200–500 μm thick with a porosity of 63–68% and surface area >10$\times$ higher than that in the unreacted shale as shown by xCT analyses. Calcite dissolution rate quantified from the attenuation data was 3 $\times$ 10 –4 mol/m 2 s and decreased to 3 $\times$ 10 –5 mol/m 2 s after 24 h because of increasing diffusion limitations. Erosion of the fracture surface increased access to iron-rich minerals, thereby increasing access to toxic metals such as arsenic. Quantification using XRF and XANES microspectroscopy indicated up to 0.5 wt % of As(-I) in arsenopyrite and 1.2 wt % of As(V) associated with ferrihydrite. This study provides valuable contributions for understanding and predicting fracture alteration and changes to the mobilization potential of hazardous metals and metalloids.

58 GEOSCIENCES↗
DOE OSTI2021

A Critical Review of the Physicochemical Impacts of Water Chemistry on Shale in Hydraulic Fracturing Systems

Hydraulic fracturing of unconventional hydrocarbon resources involves the sequential injection of a high-pressure, particle-laden fluid with varying pH’s to make commercial production viable in low permeability rocks. This process both requires and produces extraordinary volumes of water. The water used for hydraulic fracturing is typically fresh, whereas “flowback” water is typically saline with a variety of additives which complicate safe disposal. As production operations continue to expand, there is an increasing interest in treating and reusing this high-salinity produced water for further fracturing. Here in this paper we review the relevant transport and geochemical properties of shales, and critically analyze the impact of water chemistry (including produced water) on these properties. We discuss five major geochemical mechanisms that are prominently involved in the temporal and spatial evolution of fractures during the stimulation and production phase: shale softening, mineral dissolution, mineral precipitation, fines migration, and wettability alteration. A higher salinity fluid creates both benefits and complications in controlling these mechanisms. For example, higher salinity fluid inhibits clay dispersion, but simultaneously requires more additives to achieve appropriate viscosity for proppant emplacement. In total this review highlights the nuances of enhanced hydrogeochemical shale stimulation in relation to the choice of fracturing fluid chemistry.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

DECOVALEX-2023: Task A Final Report

Task A, also known as HGFrac, examines the fracturing processes that may occur in the Callovo-Oxfordian claystone (COx) in the context of the high-level (HLW) and intermediate-level long-lived (ILW-LL) radioactive waste repository in France. Understanding these processes and improving numerical models to reproduce them will aid in the design, optimization, and safety of the repository. Heat and gas fracturing are studied in two independent subtasks following a stepwise approach: laboratory tests/benchmark exercises, in-situ experiment, and finally, an application case. The in-situ heater experiment aimed to thermally induce a hydraulic fracture; temperature and pore pressure were monitored to detect any evidence of fracturing. The in-situ gas injection experiment aimed to study the effect of the stress orientation and gas injection kinetics on the gas fracturing process. The occurrence of fracturing was monitored by gas pressure measured in the injection interval. In both experiments, the excavation-induced fracture network around the heater/injection boreholes played an important role in the reduction of the compressive stress state, leading to both a tensile and shear failure response of the COx. During the first two years of the project, the research teams working on each task developed and/or proposed numerical approaches for reproducing the occurrence of fracturing in the in-situ experiments. The failure criteria were defined by reproducing the measurements from laboratory extension tests for the heat fracturing subtask. In the gas fracturing subtask, their approaches were used for simulating several benchmark exercises, and an inter-comparison between models was carried out. In both tasks, the developed approaches were compared with a simplified approach considering poro-elasticity for the mechanical behaviour of the COx. Most of the developed approaches are based on a continuous medium that takes into account variations in hydraulic properties due to mechanical degradation, such as plastic deformation or damage. Other approaches implicitly modelled weak planes or embedded discontinuities to reproduce fracture propagation. The potential for fracture initiation was also studied through of a discrete approach. In the second half of the project, the research teams mainly focused on interpretative modelling of two in-situ experiments and a blind prediction exercise to test their respective approaches. The models developed by the research teams were also applied at the repository scale to evaluate fracture initiation in a case study under vi unfavourable conditions, particularly in terms of spacing between High-Level Waste cells. The results showed that the poro-elasticity approach could be an efficient tool for understanding the main processes occurring in the COx. One example is the explicit representation of the excavation-induced fracture network around the boreholes, which yielded acceptable results compared to the measurement data. However, advanced approaches were needed to evaluate the potential increase of the excavation-induced fracture network extend and better understand fracture initiation. The stress analyses carried out by the teams revealed that hydraulic boundary conditions had a strong impact on fracture initiation in the heater experiment. Furthermore, in most cases, the results required higher pore pressure increments to reach fracturing than those measured in the experiment. This implies that the measurements may have been biased by the packer’s capacity to fully isolate the piezometric chambers, leading to lower pressures. On the contrary, there was no agreement on the fracturing mode, as some reported either shear or tensile fracturing, while others reported a combination of the two modes. In the case of gas fracturing, the research teams were limited to the comparison of a single point, which complicated their task. Nonetheless, the numerical results were able to reproduce the measurements and capture processes such as longitudinal gas flow through the excavation-induced fracture network, as suggested by some evidence in the observation piezometric chambers. The numerical models also agreed with the measurements in the sense of higher probability of developing along the injection borehole than radially towards the sound rock. The approaches developed by the research teams showed that they are capable of analysing and reproducing fracture initiation in the COx. However, areas of future work should focus on the fracture propagation and fracture aperture, which were out of the scope of this task. To this end, additional data must be gathered for the parameter characterisation and validation of the numerical models. Nonetheless, various approaches showed promising results as they were able to reproduce fracture development under certain conditions.

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

Concurrent measurement of strain and chemical reaction rates in a calcite grain pack undergoing pressure solution: Evidence for surface-reaction controlled dissolution

Pressure solution is inferred to be a significant contributor to sediment compaction and lithification, especially in carbonate sediments. For a sediment deforming primarily by pressure solution, the compaction rate should be directly related to the rate of calcite dissolution, transport along grain contacts, and calcite reprecipitation. Previous experimental work has shown that there is evidence that deformation in wet calcite grain packs is consistent with control by pressure solution, but considerable ambiguity remains regarding the rate limiting mechanism. We present the results of laboratory compaction experiments designed to directly measure calcite dissolution and precipitation rates (recrystallization rates) concurrently with strain rate to test whether measured rates are consistent with predicted rates both in absolute magnitude and time evolution. Recrystallization rates are measured using trace element chemistry (Sr/Ca, Mg/Ca) and isotopes (87Sr/86Sr) of fluids flowing slowly through a compacting grain pack as it is being triaxially compressed. Imaging techniques are used to characterize the grain contacts and strain effects in the post-experiment grain pack. Our data show that calcite recrystallization rates calculated from all three geochemical parameters are in approximate agreement and that the rates closely track strain rate. The geochemically inferred rates are close to predicted rates in absolute magnitude. Uncertainty in grain contact dimensions makes distinguishing between surface reaction control and diffusion control difficult. Measured reaction rates decrease faster than predicted from standard pressure solution creep flow laws. This inconsistency may indicate that calcite dissolution rates at grain contacts are more complex, and more time-dependent, than suggested by geometric models designed to predict grain contact stresses.

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

Advective gas flow in bentonite: Development and comparison of enhanced multi-phase numerical approaches

Understanding the impact of repository gas, generated from degradation of waste and its interaction with the host rock, is essential when assessing the performance and safety function of long-term disposal systems for radioactive waste. Numerical models based on conventional multi-phase flow theory have historically been applied to predict the outcome and impact of gas flow on different repository components. However, they remain unable to describe the full complexity of the physical processes observed in water-saturated experiments (e.g., creation of dilatant pathways) and thus, the development of novel representations for their description is required when assessing fully saturated clay-based systems. This was the primary focus of Task A within the international cooperative project DECOVALEX-2019 (D-2019) and refinement of these approaches is the primary focus of this study (Task B in the current phase of DECOVALEX-2023). This paper summarises development of enhanced numerical representations of key processes and compares the performance of each model against high-quality laboratory test data. Experimental data reveals that gas percolation in water-saturated compacted bentonite is characterised by four key features: (i) a quiescence phase, followed by (ii) the gas breakthrough, which leads to a (iii) peak value, which is then followed by (iv) a negative decay. Three models based on the multiphase flow theory have been developed. These models can provide good initial values and reasonable responses for gas breakthrough (although some of them still predict a too-smooth response). Peak gas pressure values are in general reasonably well captured, although maximum radial stress differences are observed at 48 mm from the base of the sample. Here, numerical peak values of 12.8 MPa are predicted, whereas experimental values are about 11 MPa. These models are also capable of providing a reasonable representation of the negative pressure decay following peak pressure. However, other key specific features (such as the timing of gas breakthrough) still require a better representation. The model simulations and their comparison with experimental data show that these models need to be further improved with respect to model parameter calibration, the numerical representation of spatial heterogeneities in material properties and flow localisation, and the upscaling of the related physical processes and parameters. To further understand gas flow localisation, a new conceptual model has been developed, which shows that discrete channels can possibly be induced through the instability of gas-bentonite interface during gas injection, thus providing a new perspective for modeling gas percolation in low-permeability deformable media.

58 GEOSCIENCES↗