Search

In the 2023 phase of the international collaborative DECOVALEX modeling project, Task E focused on understanding thermal, hydrological, and mechanical (THM) processes related to predicting brine migration in the excavation damaged zone around a heated excavation in salt. Salt is attractive as a disposal medium for radioactive waste because it is self-healing and is essentially impermeable and non-porous in the far field. Investigation of the short-term, near-field behavior is important for radioactive waste disposal because this early period strongly controls the amount of inflowing brine. Brine leads to corrosion of waste forms and waste packages, and possible dissolution of radionuclides with brine transport being a potential transport vector to the accessible environment.The Task was divided into steps. Step 0 included matching unheated brine inflow data from boreholes at the Waste Isolation Pilot Plant (WIPP) and matching temperature observations during a Brine Availability Test in Salt (BATS) heater test. Step 1 included validation of models against a thermo-poroelastic analytical solution, and two-phase flow around an excavation. Finally, Step 2 required all the individual components covered in steps 0 and 1 to come together to match observed brine inflow behavior during the same BATS heater test.There were a range of approaches from the teams, from mechanistic to prescriptive. Given the uncertainties in the problem, some teams used one- or two-dimensional models of the processes, while other teams included more geometrical complexity in three-dimensional models. Task E was a learning experience for the teams involved, and feedback from the modeling teams has led to changes in follow-on BATS experiments at WIPP. The primary Task E lessons learned were the impact of hydrologic initialization methods (wetting up vs. drying down), the difference between confined and unconfined thermal expansion, and the large changes in permeability associated with heating and cooling.

Induced seismicity poses a challenge to the development of Enhanced Geothermal Systems (EGS). Improving monitoring and forecasting techniques is essential to mitigate induced seismicity and thereby fostering a positive perception of EGS projects among local authorities and population. Induced seismicity is the result of complex and coupled thermo-hydro-mechanical-chemical mechanisms. Injection flux and pressure are crucial controlling parameters for both hydraulic stimulation and circulation protocols. We develop a methodology combining a hydro-mechanical model with a seismicity rate model to estimate the magnitude and frequency of mainshocks and aftershocks induced by fluid injection. We apply the methodology to the case of the Basel EGS (2006, Switzerland) to compare the effects of progressive, cyclic and constant injections on the mechanical response of discrete faults. Results from the coupled hydro-mechanical models show that the pore pressure diffusion and consequent enhancement of fault permeability are limited to the vicinity of the injection well during cyclic injection. Additionally, constant injection induces seismicity from the start of the injection but enhances the permeability of most of the faults within a shorter duration, inducing less post-injection seismicity. The methodology can be adapted to any numerical model and allows new projects to be developed by anticipating the safest injection protocol.This article is part of the theme issue 'Induced seismicity in coupled subsurface systems'.

Shearing of fractures and faults is important because it can result in permeability change or even induce seismicity—both are keys for efficient and safe energy recovery and storage in Earth systems. Quantitative analysis of shearing of intersecting fractures and faults is challenging because it can involve dynamic frictional contacts that are complicated by deformation of the rock matrix. To predict the shearing of intersecting fractures/faults, we attempt to answer the question of how intersections impact the shearing of a fracture network and whether we can simplify the description as compared to classical discrete fracture network (DFN) models. To answer these questions, we conducted a series of numerical simulations on scenarios for variable numbers of intersecting fractures. All these examples yield consistent results: the results of using DFNs are consistent with those of using hypothetical major paths. This leads to a new model, which we name simplified discrete fracture network model, to analyze shearing of intersecting fractures/faults using major path(s). We found that the intersections of fractures do not fundamentally change the shearing of two intersecting fractures if the intersecting angles are small. Furthermore, increasing the number of fractures/faults may relax the stress as more fractures/faults become available for shearing and distributing the stress. The simplified DFN model, which can capture efficiently the shearing behavior of each major paths from a large number of intersecting fractures/faults, will be a promising conceptual model that is complementary to existing equivalent continuum and discrete fracture models to analyze shearing of intersecting fractures/faults.

This research paper focuses on the thermo-hydro-mechanical (THM) modeling of brine migration in a heated borehole test conducted as part of the ongoing Brine Availability Test in Salt (BATS) at the Waste Isolation Pilot Plant (WIPP) in New Mexico. It is a component of the international collaboration project DECOVALEX-2023 (DEvelopment of COupled models and their VALidation against EXperiments), which aims to understand the THM processes governing brine flow in heated rock salt repositories through collaborative analysis by multiple research teams. Using the TOUGH–FLAC simulator, THM simulations were performed and compared with data from the BATS phase 1a. This experiment involved two identical horizontal-borehole arrays, one heated and one serving as a control, both equipped with sensor arrays. Analysis of measurements revealed water flow rate surges during heater power transitions, with the highest jump observed during cooling. Acoustic emission activity exhibited distinct patterns in response to heater power changes, suggesting that damage to rock salt is particularly pronounced during the cooling phase. The THM simulations successfully captured these phenomena, highlighting the significance of thermal effects, brine migration, and mechanical behavior in predicting brine availability in heated and damaged rock salt. Our modeling also revealed the critical interplay between heating and cooling-induced damage and its influence on flow properties, particularly affecting brine inflow estimation. Notably, we found that cooling-induced brine inflow spikes result from increased permeability due to tensile dilatancy. These findings have important implications for the development of robust containment strategies and enhance our understanding of the complex processes involved in repository performance.

To ensure the safety of geological disposal of high-level radioactive waste, in-situ experiments have been carried out to examine the behavior of rocks in underground research laboratories (URLs). At the Meuse/Haute-Marne URL in France, the French National Radioactive Waste Management Agency (Andra) has been assessing the Callovo-Oxfordian claystone (COx) as potential host rock of geological disposal by subjecting the COx to in-situ heating mimicking exothermic radioactive waste. Results of the in-situ experiments are used to validate and bolster the numerical simulators for predicting the thermo-hydromechanically (THM) coupled behavior of the COx. The numerical simulators are, however, yet to be tested for predicting the failure and fracture development of the COx during heating, which is of paramount importance to the safety of the geological disposal. In this research, we modelled a recently carried out in-situ experiment at the Meuse/Haute-Marne URL using the TOUGH-FLAC simulator to predict the failure and fracture development of the COx during heating. The objectives are to examine the effects of (i) the weak bedding planes, (ii) the softening rate of matrix/weak plane strengths, and (iii) the stiffness anisotropy of the COx on the development of shear and tensile fractures during heating. Results show that considering failure along the weak planes enabled accurate predictions of fracture development. Also, fracture development intensified at a softening rate beyond a threshold level and the geometry of fractures was significantly affected by the stiffness anisotropy. These results will help boost the reliability of the safety and performance assessment of geological disposal in claystone.

We conducted two-dimensional numerical simulations to investigate the mechanisms underlying the strong spatiotemporal correlation observed between submarine landslides and gas hydrate dissociation due to glacial sea-level drops. Our results suggest that potential plastic deformation or slip could occur at localized and small scales in the shallow-water portion of the gas hydrate stability zone (GHSZ). This shallow-water portion of the GHSZ typically lies within the area enclosed by three points: the BGHSZ–seafloor intersection, the seafloor at ∼600 m below sea level (mbsl), and the base of the GHSZ (BGHSZ) at ∼1,050 mbsl in low-latitude regions. The deep BGHSZ (>1,050 mbsl) could not slip; therefore, the entire BGHSZ was not a complete slip surface. Glacial hydrate dissociation alone is unlikely to cause large-scale submarine landslides. Observed deep-water (much greater than 600 mbsl) turbidites containing geochemical evidence of glacial hydrate dissociation potentially formed from erosion or detachment in the GHSZ pinch-out zone.

To ensure long-term safety and performance, geological nuclear waste repositories require low-permeability barriers such as bentonite buffers and/or shale host rock. Shale is not only known for its low permeability but also for its trend to undergo time-dependent deformation (i.e., creep), which could heal damage, but the effects of shale creep on the long-term performance of nuclear waste repositories have not been clearly understood. In particular, the anisotropic nature of shale (i.e., bedding) could have a significant effect on its creep behavior, and consequently, on the long-term performance of nuclear waste repositories. In this research, numerical simulations were carried out with the objective of showing the effects of anisotropic shale creep on the stress and permeability evolution of a generic geological nuclear waste repository in shale. The TOUGH-FLAC simulator, a thermo-hydromechanically (THM) coupled numerical code, was used for the simulations. To achieve the objective, comparisons were performed between the results of anisotropic shale creep simulations and those of different simulation cases, namely, no creep (i.e., elastic), isotropic creep, and long-term creep shale cases. Results of the comparisons show that the elastic and isotropic creep shale cases respectively led to the overestimation and underestimation of stress and permeability in the repository, whereas the long-term creep shale case, which accumulated greater creep in later periods than in earlier periods, helped to keep large shear and tensile stresses from developing while maintaining compressive spherical stress, resulting in consistently low permeability levels. These results indicate that performance assessments with elastic and isotropic creep formation models will provide the upper and lower bound estimates of stress and permeability, while more reasonable estimates will be given by an anisotropic creep formation model, and that shale with long-term creep characteristics will be beneficial in many aspects of the safety and performance of nuclear waste repositories.

Quick estimates of fluid-induced subsurface deformation are helpful to assess the land uplift/subsidence and to reveal precursors of induced seismicity. Here, we adopt the inclusion theory and Green's function to develop a closed-form solution in a half space for the poroelastic response of a reservoir compartmentalized by an intersecting fault that can be offset and either permeable or impermeable. Simulated results reveal that (1) fault permeability mainly impacts the spatial distribution of displacement while its effect on displacement magnitude is small; (2) ground displacement slightly increases with fault dip while slightly decreases with increasing fault offset; in contrast, reservoir geometry shows a stronger effect than fault geometry: the ground displacement is proportional to the vertical and lateral depth ratios, defined as the ratios of reservoir thickness (h) and width (w) to reservoir depth (D), respectively; (3) the maximum vertical displacement is the double of the horizontal one regardless of fault permeability, fault and reservoir geometries, and mechanical parameters. Comparing the solution in a half space with that in a full space shows that neglecting the free surface underestimates the poroelastic displacement in the overburden. The validity of full-space solutions can be assessed with the product of the lateral and vertical depth ratios, i.e., wh/D2. The full-space solutions become valid when wh/D2 decreases to an intrinsic threshold. This threshold may range from 0.01 to 0.02 for displacement, and its specific value depends on the field background and demands of projects, but can be estimated based on our solution. It is larger for stress than for displacement, and wh/D2 ≤ 0.1 is recommended as a general condition for neglecting the free-surface effects on induced stress. The analytical solution represents a useful tool for estimating ground deformation and for gaining insights of reservoir and fault geometries by analyzing surface deformation patterns.

Pressure solution, a mechanism that involves tight coupling between the geometry and thermal-hydro-mechanical-chemical (THMC) processes, plays an important role in diagenesis. In this study, we make the first attempt to conduct microscale THMC modeling to understand and quantify the impacts of geometry and temperature on pressure solution, taking natural salt rock as an example. This modeling capability is achieved by expanding a novel MC code that we developed previously (Hu et al. J Geophys Res: Solid Earth 126:e2021JB023112, 2021) to include temperature effects. We first conduct a simulation of an example that involves a single brine inclusion within a single halite grain and find that the temperature impact is limited for that case. We then extract geometry from an image of a natural salt rock and conduct simulations with different cases: (A) only temperature and no stress, (B) only stress and no temperature, and (C) with both stress and temperature. These different cases result in quite different phenomena. In case A, dissolution and precipitation occur across the entire system due to isolated pore space reaching a localized mass balance between dissolution, precipitation, and diffusion. In case B, intense geometric features (e.g., major asperities, inclusions) in one area undergo stress concentration, thus dominating pressure solution in that area. In case C, pressure solution is spread out at contacting highly stressed geometric features close to the hotter side. We conclude that geometric features dominate stress distribution, thus dominating pressure solution in a natural salt rock that may be affected by the temperature if a sufficient temperature gradient is applied.

Salt formations are a promising host rock for the disposal of heat-emitting nuclear waste due to their ability to heal fractures, high thermal conductivity, and near-zero permeability. This study investigates the concept of waste emplacement in salt formations by placing waste packages in drifts and backfilling them with crushed salt. The aim is to minimize fracturing and restore the initial properties of the host rock. One of the key issues in assessing the safety and performance of nuclear waste repositories in salt is determining the reconsolidation rate of the crushed salt backfill, which has been traditionally predicted using the C-WIPP model. However, this model was calibrated using data from uniaxial oedometer lab tests, where specimens with high initial porosity values have been compacted to medium porosity values. The observed compaction behavior has then been extrapolated to the low porosity range, leading to a prediction of rapid full reconsolidation of the backfill in less than 20 years. In contrast, the new C-WIPP/TUC model used in this study was calibrated for the medium porosity range, between 17% and 8%, based on more suitable triaxial lab tests. Using coupled Thermo-Hydro-Mechanical (THM) simulations with the TOUGH-FLAC simulator, it was found that the C-WIPP/TUC model predicts a slower reconsolidation rate of the backfill compared to the C-WIPP model, with an average porosity of the backfill remaining roughly 10 times that of natural salt after 10,000 years. The effects of the relatively slow reconsolidation of the crushed salt backfill on the THM behavior of the repository are also examined, highlighting the importance of accurately capturing the behavior of crushed salt for the study of the long-term integrity of a nuclear waste repository in salt.

The Caney shale is an emerging hydrocarbon play located in southwest Oklahoma, USA. Within the Caney shale exist facies which were initially dubbed “reservoir” and “ductile” based on evaluation of well logging data. While past work has shown the distinction of “brittle” and “ductile” is not mechanically justifiable according to formal definitions, the current work shows some important differences between nominally ductile and reservoir zones. First, the “ductile” zones are more clay rich and have textural differences which can be expected to lead to differences in mechanical properties. One important impact of these differences is observed in triaxial creep experiments showing the “ductile” zones are more prone to creep deformation. Numerical simulations predict the “reservoir” zones will experience very little proppant embedment due to creep deformation of hydraulic fractures around proppant particles. On the other hand, “ductile” zones can be expected to undergo creep-driven proppant embedment leading to loss of fracture aperture ranging up to 100% loss, depending upon the spatial density of the proppant distribution. Hence, this research shows the identification of nominally “ductile” zones from well logs, while a misnomer, can be useful in finding clay-rich, creep-prone zones which will be the most prone to proppant embedment and hence vulnerable to greater production decline over time.

Rock salt has a self-sealing capacity, low permeability, and high thermal conductivity, making it a potential host for heat-generating nuclear waste. The feasibility of nuclear waste disposal within salt formations has been investigated mostly for small-sized canisters. Geologic disposal of larger-sized canisters originally designed for spent fuel storage and transportation has lately been examined as a cost-effective alternative. This raises questions about their long-term vertical movement due to their weight and high decay heat. Low-stress creep governs this movement; however, most salt constitutive models do not incorporate it. In this paper, the Norton and the WIPP creep models are compared with the Lux/Wolters/Lerche (LWL) model and a simpler model that combines linear and Norton creep laws (named combined creep model). The LWL and combined creep models consider pressure solution creep, though all incorporate dislocation creep. The models are first applied to creep tests under various stress levels. The LWL and the combined creep models results fit the experimental data well in both high and low stress ranges, whereas the Norton and WIPP models results only fit in higher-stress ranges. The different models are further applied for analyzing long-term canister movement. A sinking rate of - 4.4 × 10 - 7 mm/year was predicted using the Norton and WIPP models versus - 2.1 × 10 - 2 mm/year and - 3.1 × 10 - 2 mm/year using the LWL and the combined creep models, respectively. This comparative study confirms that creep models calibrated exclusively against high-deviatoric stress data might result in an inaccurate estimation of waste packages sinking rate in salt formations.

The current worldwide energy supply is insufficient to meet the rising demand. As a result, the energy prices are expected to keep soaring despite the recent increases in a variety of renewable energy resources. Although not renewable, shale oil and gas — “unconventional” hydrocarbon resources are relatively clean forms of energy resources, which still hold a vast share of the energy market. For many oil and gas companies, meeting profitable production goals from shale reservoirs is sometimes challenging, due to the loss of fracture conductivity and premature declines in the production. In this paper we investigate the stress-dependent changes in the hydraulic conductivity of proppant-filled fractures and mechanical fracture–proppant interactions in Caney Shale, a calcareous, organic-rich mudrock, through laboratory experiments and numerical modeling. American Petroleum Institute (API) fracture conductivity tests were conducted using 2% KCl on five locations within the Caney Shale that consisted of selecting three brittle (reservoir) zones and two ductile zones. Confining pressures ranged from 1,000 psi to 12,000 psi at 210 °F. Conductivity, permeability as well as embedment were measured during the test. Also, an additional, laboratory in-situ visualization test was conducted to examine the detailed proppant-shale matrix interaction under elevated stress (3,920 psi effective stress) and temperature (252 °F), with a synthetic reservoir fluid. Our experimental results have confirmed that improved fracture conductivity is attributed to proppant size, and that the increase in porosity of the proppant pack, closure pressure changes and the reduction in fracture conductivity are a function of many factors such as fracture closure stress.

Geologic storage of carbon dioxide can efficiently contribute to reduce greenhouse gas emissions to the atmosphere. Two major hazards of this technology are leakage towards the ground surface and fluid-induced seismicity. While major faults may be detected and avoided during site characterization, undetected subseismic faults could be encountered once injection has started. This paper investigates leakage and reactivation of undetected faults through coupled thermal–hydraulic–mechanical modeling. The simulations are performed using a recently developed sequential simulator, TOUGH-Pylith, that allows to accurately account for the thermodynamics of brine-CO2 mixtures, and to model faults as surfaces of discontinuity using state-of-the-art fault friction laws. The simulator is benchmarked against well-known analytical solutions and subsequently applied to investigate two cases of CO2 injection close to undetected faults under normal and strike-slip faulting regimes. Although the scenarios are generic and represent unfavorable conditions, they suggest that leakage could occur and that undetected faults could trigger minor seismic events. Therefore, careful site characterization and continuous monitoring during operations should be always performed.

The recent development of the TOUGH3 code allows for a faster and more reliable fluid flow simulator. At the same time, new versions of FLAC3D are released periodically, allowing for new features and faster execution. In this paper, we present the first implementation of the coupling between TOUGH3 and FLAC3Dv6/7, maintaining parallel computing capabilities for the coupled fluid flow and geomechanical codes. We compare the newly developed version with analytical solutions and with the previous approach, and provide some performance analysis on different meshes and varying the number of running processors. Finally, we present two case studies related to fault reactivation during CO2 sequestration and nuclear waste disposal. The use of parallel computing allows for meshes with a larger number of elements, and hence more detailed understanding of thermo-hydro-mechanical processes occurring at depth.

Safety assessment of geological nuclear waste repositories is essential for the permanent disposal of spent nuclear fuels and high-level radioactive waste. The long-term integrity of the host rock as well as the engineered barrier system (e.g., bentonite buffer) surrounding nuclear waste canisters is of particular importance as decay heat from nuclear waste canisters, which will last over thousands of years, significantly disturbs the thermo-hydromechanical (THM) state of the repository. In this study, THM coupled simulations were carried out to investigate the effect of time-dependent deformation (i.e., creep) of shale on the long-term integrity of a generic subsurface nuclear waste repository. The Norton-Bailey creep model, which is also known as the Lemaitre-Menzel-Schreiner model, was employed to simulate the power-law type creep that is observed in shales. The TOUGH-FLAC simulator was employed for the THM coupled modeling of the repository. The objective of this study is to assess the effect of creep in different shales (i.e., high creep shale vs. low creep shale) on long-term stress and permeability changes in the repository. Results show potential advantages of constructing repositories in the high creep shale, as deviatoric stress levels in the formation decreased to zero in 100 years since the emplacement of nuclear waste canisters and the permeability also decreased to the undamaged, intrinsic levels in 10,000 years. Also, mean effective stress levels in the bentonite buffer increased by 100% in the high creep shale case relative to the low creep shale case at 10,000 years due to creep-induced contraction of the nuclear waste disposal tunnel. However, in earlier periods (e.g., 1000 years), the stress levels in the bentonite buffer were twice smaller in the high creep shale case than in the low creep shale case, which shows a tradeoff between the intermediate- (∼1000 years) and long-term (>10,000 years) compaction levels in the bentonite buffer depending on the creep characteristics of the host shale.

The greatest challenges of rigorously modeling coupled hydro-mechanical processes in fractured rocks at different scales are associated with computational geometry. In addition, selections of continuous or discontinuous models, physical laws, and coupling priorities at different scales based on different geometric features determine the applicability of a numerical model for a certain type of problem. In this study, we present our multi-scale modeling capabilities that have been developed based on the numerical manifold method for analyzing coupled hydro-mechanical processes in fractured rocks. Based on their geometric features, the fractures are modeled as continua—finite-thickness porous zones, and discontinua—discontinuous interfaces and microscale asperities and granular systems. Different governing equations, physical laws, coupling priorities, and approaches for addressing fracture intersections and shearing are then applied to describe these. We applied these models to simulate coupled processes in fractured rocks using realistic geometry obtained from rock images at different scales. We first calculated shearing of a single fracture with different models and demonstrated the impacts of asperities on shearing. We then applied the continuous and discontinuous models to simulate a network of rough fractures, demonstrating that contact dynamics contribute significantly to the geometric, multi-physical evolution of systems where rough fractures are not mineral filled. For a discrete fracture network, our coupled processes modeling demonstrates that shearing of the discrete fractures can have a major impact on stress and pore pressure distribution. Lastly, we applied the discontinuous granular model to simulate evolution of a complex granular system with a deformation band, demonstrating that the deformation band can dominate contact dynamics, the structural and the stress evolution of the granular system.

Injection of CO2 for geologic carbon sequestration into deep sedimentary formations involves fluid pressure increases that engage hydromechanical processes that can cause seismicity by activation of existing faults. In this work, we use a coupled multiphase fluid flow and geomechanical simulator to model spatiotemporal fluid pressure and stress changes in order to study the poroelastic effect of CO2 injection on faults in crystalline basement rock below the injection zone. The seismicity rate along features interpreted to be basement faults is modeled using Dieterich's rate-and-state earthquake nucleation model. The methodology is applied to microseismicity detected during CO2 injection into the Mount Simon formation during the Illinois Basin—Decatur Project. The modeling accurately captures an observed reduction in seismicity rate when the injection in the second well was into a slightly shallower zone above the base of the Mount Simon formation. Moreover, the modeling shows that it is important to consider poroelastic stress changes, in addition to fluid pressure changes for accurately modeling of the observed seismicity rate.

Deep geothermal energy (DGE) represents an opportunity for a sustainable and carbon-free energy supply. One of the main concerns of DGE is induced seismicity that may produce damaging earthquakes, challenging its widespread exploitation. It is widely believed that the seismicity risk can be controlled by using doublet systems circulating water to minimize the injection-induced pressure changes. However, cold water reinjection may also give rise to thermal stresses within and beyond the cooled region, whose potential impacts on fault reactivation are less well understood. Here, we investigate by coupled thermo-hydro-mechanical modeling the processes that may lead to fault reactivation in a hot sedimentary aquifer (HSA) in which water is circulated through a doublet. We show that thermal stresses are transmitted much ahead of the cooled region and are likely to destabilize faults located far away from the doublet. Meanwhile, the fault permeability mainly controls the fault reactivation timing, which entails the importance of employing appropriate characterization methods. This investigation is crucial for understanding the mechanisms controlling induced seismicity associated with DGE in a HSA and allows the success of future DGE projects.

Subsidence induced by groundwater depletion is a grave problem in many regions around the world, leading to a permanent loss of groundwater storage within an aquifer and even producing structural damage at the Earth's surface. California's Tulare Basin is no exception, experiencing about a meter of subsidence between 2015 and 2020. However, understanding the relationship between changes in groundwater volumes and ground deformation has proven difficult. We employ surface displacement measurements from Interferometric Synthetic Aperture Radar (InSAR) and gravimetric estimates of terrestrial water storage from the Gravity Recovery and Climate Experiment (GRACE) satellite pair to characterize the hydrological dynamics within the Tulare basin. The removal of the long-term aquifer compaction from the InSAR time series reveals coherent short-term variations that correlate with hydrological features. For example, in the winter of 2018-2019 uplift is observed at the confluence of several rivers and streams that drain into the southeastern edge of the basin. These observations, combined with estimates of mass changes obtained from the orbiting GRACE satellites, form the basis for imaging the monthly spatial variations in water volumes. This approach facilitates the quick and effective synthesis of InSAR and gravimetric datasets and will aid efforts to improve our understanding and management of groundwater resources around the world.

As a direct outcome of economic development coupled with an increase in population, global energy demand will continue to rise in the coming decades. Although renewable energy sources are increasingly investigated for optimal production, the immediate needs require focus on energy sources that are currently available and reliable, with a minimal environmental impact; the efficient exploration and production of unconventional hydrocarbon resources is bridging the energy needs and energy aspirations, during the current energy transition period. The main challenges are related to the accurate quantification of the critical rock properties that influence production, their heterogeneity and the multiscale driven physico-chemical nature of rock–fluid interactions. A key feature of shale reservoirs is their low permeability due to dominating nanoporosity of the clay-rich matrix. As a means of producing these reservoirs in a cost-effective manner, a prerequisite is creation of hydraulic fracture networks capable of the highest level of continued conductivity. Fracturing fluid chemical design, formation brine geochemical composition, and rock mineralogy all contribute to swelling-induced conductivity damage. The Caney Shale is an organic-rich, often calcareous mudrock. Many studies have examined the impact that clay has on different kinds of shale productivity but there is currently no data reported on the Caney Shale in relation to horizontal drilling; all reported data on the Caney Shale is on vertical wells which are shallow, compared to an emerging play that is at double the depth. In this work we develop geochemical–geomechanical integration of rock properties at micro-and nanoscales that can provide insights into the potential proppant embedment and its mitigation. The novel methodology amalgamates the following: computed X-ray tomography, scanning electron microscopy, energy dispersive spectroscopy, micro-indentation, and Raman spectroscopy techniques. Our results show that due to the multiscale heterogeneity in the Caney Shale, these geochemical and structural properties translate into a variation in mechanical properties that will impact interaction between the proppant and the host shale rock.

In the study, we analyze changes in groundwater pressure observed in several boreholes drilled in and around the Mizunami Underground Research Laboratory (MIU) induced by the 2011 off the Pacific coast of Tohoku Earthquake (Mw 9.0). The aim of this project is a development of methodology to evaluate systematic fault activity by numerical analysis. To reach this goal we investigate the behavior of the fault zones present in the area during the passing of seismic waves. We built a simplified hydrogeological model of the MIU site and performed a series of fluid flow simulations with TOUGH2 flow numerical code. We investigate how changes in permeability along three faults present in the study area: the Tsukiyoshi Fault, the Hiyoshi Fault and the Main-Shaft Fault may have influence the groundwater level monitored in boreholes intervals. We also test the influence of the cone of depression at the MIU site and the hydraulic connectivity between the sedimentary cover and the granite aquifers. Our results suggest that two main mechanisms are responsible for the observed changes in groundwater pressure: (1) crustal dilation induced by the Tohoku earthquake causing a groundwater recharge from the sedimentary aquifers to the Toki granite aquifer where the sedimentary cover is thick; and (2) permeability increase along faults critically oriented for shear reactivation and oriented in the direction of the passing seismic wave. In this case, the seismic wave increases the shear stress acting on the fault promoting slip and a change in permeability through a mechanism of slip-induced dilation. Faults not critically stressed and faults critically oriented for shear reactivation but oriented perpendicular to the passing seismic wave are not reactivated.

Microscale numerical modeling is potentially an effective approach for understanding salt creep, but it is quite challenging because natural salt grain boundaries can have complex and evolving geometry, and because contacts between these deformable salt mineral grains are dynamic as a result of compaction, chemical reaction, fluid flow, and heat transfer. In this study, we have overcome these challenges and developed a new microscale mechanical-chemical (MC) model to analyze creep of salt at the microscale, accounting for coupled deformation, dynamic contacts, and chemical reaction in granular systems with realistic geometric representations. The MC model was realized by linking a new microscale mechanical code based on the numerical manifold method (NMM) to a reactive transport code named Crunch. We simulated the processes of reorganization of the salt grains, microfracturing, and pressure solution that contribute to creep at larger scales. Based on this first quantitative microscale model, we found that sharp corners of mineral grains can dominate the contact dynamics, microfracturing, and pressure solution when salt is compacted, thus governing the structural changes and porosity loss of the system. We found that pressure solution, which preferentially dissolves sharp corners and edges, can lead to relatively high porosity loss in the system, thus playing an important role in the creep of salt. Our analysis shows that the dynamic changes of salt granular systems involving grain relocation and pressure solution can occur repeatedly and continuously, thus contributing to longer-term creep of salt at larger scales.

In support of the U.S. DOE SubTER Crosscut initiative, we established a field test facility in a deep mine and designed and carried outin situ hydraulic fracturing experiments relevant to enhanced geothermal systems (EGS) in crystalline rock to characterize the stressfield, understand the effects of rock fabric on fracturing, and gain experience in monitoring using geophysical methods. The project alsoincluded pre- and post-fracturing simulation and analysis, and laboratory measurements and experiments. The kISMET (permeability(k) and Induced Seismicity Management for Energy Technologies) site was established in the West Access Drift of the SanfordUnderground Research Facility (SURF) 4757 ft (1450 m) below ground (on the 4850 ft level (4850L)) in phyllite of the PrecambrianPoorman Formation. We drilled and continuously cored five near-vertical boreholes in a line on 3 m (10 ft) spacing, deviating the twooutermost boreholes slightly to create a five-spot pattern around the test borehole centered in the test volume 40 m below the drift invert(floor) at a total depth of ~1490 m (4890 ft). Laboratory measurements of core from the center test borehole showed P-wave velocityheterogeneity along each core indicating strong, fine-scale (~1 cm or smaller) changes in the mechanical properties of the rock. Fieldmeasurements of the stress field by hydraulic fracturing showed that the minimum horizontal stress at the kISMET site averages 21.7MPa (3146 psi) trending approximately N-S (356 degrees azimuth) and plunging slightly NNW at 12°. The vertical and horizontalmaximum stresses are similar in magnitude at 42-44 MPa (6090-6380 psi) for the depths of testing, which averaged approximately 1530m (5030 ft). Hydraulic fractures were remarkably uniform suggesting core-scale and larger rock fabric did not play a role in controllingfracture orientation. Analytical solutions suggest that the fracture radius of the large fracture (stimulation test) was more than 6 m (20ft), depending on the unknown amount of leak-off.

This paper provides a comprehensive review on the application of proppants to maintain fracture permeability over the lifetime of a well based on published observations from experiments and modeling. The review identifies and describes important processes occurring during proppant embedment, during hydraulic fracturing, laboratory testing of fracture conductivity, proppant embedment and modeling of proppant embedment. Finally, this paper identifies the challenges and knowledge gaps that also provide future avenues of research and opportunities for collaborative technological development which requires an interdisciplinary approach of science, engineering in academia, government, and private sector.

Observations on tens-of-meter scale experiments of fault activation by fluid injection conducted in shales allow exploring how aseismic and seismic events may jeopardize the integrity of a sealing caprock overlying a CO2 sequestration reservoir. We contrast the behavior of shale faults with another set of experiments conducted in carbonates. Significant fluid leakage occurs along the initially low-permeability shale faults when rupture is activated. Most of the leakage pathway closes when fluid injection ceases and fluid pressure drops. Dilatant slip on the fault plane alone does not explain the observed leakage behavior, which is also caused by fault opening favored by the softness of the shale, and by the structure of the fault zone that prevents fluids from diffusing into the adjacent damage zone. Experiments show a large amount of aseismic deformation. Small-magnitude seismicity (Mw < -2.5) is observed outside the pressurized leakage patch. Stress transferred from this aseismic deformation patch can build up to stress-criticality and favor seismicity. Thus, in terms of fault activation in caprocks, aseismic fault slip leading to increased permeability and a loss of seal integrity is of great concern.

Oil and gas shales are a system where understanding the mechanical properties at the microscale is of paramount importance, e.g. to better understand the behavior of proppant-shale contacts and their role in the evolution of propped fractures in unconventional reservoirs. This work shows for the first time an experiment coupling indentation testing with in situ X-ray imaging in a Green River shale sample. A full compliance curve has been measured with the sample in water, allowing to visualize the indentation of the sample in function of axial load, in a purpose-built system for combined in situ indentation and X-ray imaging. A series of 3D datasets were used for a digital volume correlation study to obtain local strain fields. This analysis has been complemented with the analysis of cracks. Finally, geomechanical modeling has been carried out to replicate and generalize the observed behavior in the shale. This study validated this experimental approach, providing a breakthrough in understanding micro-mechanics in shales, and demonstrates how this class of experiments can be important for studies involving the prediction of the evolution of propped fractures in shale reservoirs, with possible applications in a much larger number of application fields (geothermal, materials science, etc.)

The complex geometric features of subsurface fractures at different scales makes mesh generation challenging and/or expensive. In this paper, we make use of neural style transfer (NST), a machine learning technique, to generate mesh from rock fracture images. In this new approach, we use digital rock fractures at multiple scales that represent ‘content’ and define uniformly shaped and sized triangles to represent ‘style’. The 19-layer convolutional neural network (CNN) learns the content from the rock image, including lower-level features (such as edges and corners) and higher-level features (such as rock, fractures, or other mineral fillings), and learns the style from the triangular grids. By optimizing the cost function to achieve approximation to represent both the content and the style, numerical meshes can be generated and optimized. We utilize the NST to generate meshes for rough fractures with asperities formed in rock, a network of fractures embedded in rock, and a sand aggregate with multiple grains. Based on the examples, we show that this new NST technique can make mesh generation and optimization much more efficient by achieving a good balance between the density of the mesh and the presentation of the geometric features. Finally, we discuss future applications of this approach and perspectives of applying machine learning to bridge the gaps between numerical modeling and experiments.

This paper presents a coupled thermal-hydraulic-mechanical (THM) model for simulating failure processes in geological media. The coupling between stress/strain and heat transfer describing how thermal expansion affects stress is based on the theory of thermo-elasticity. Additionally, the coupling between stress/strain and fluid flow describing how pore-pressure affects stress is governed by Biot's consolidation theory. Stress/strain-dependent porosity/permeability/capillary-pressure defines the influence of stress on fluid flow and thermal fields. The THM model is implemented in the solver module that links the multiphase, multicomponent fluid flow simulator TOUGH2 with the geomechanical simulator RFPA (Rock Failure Process Analysis), which is based on the concept that heterogeneity leads to non-linearity and progressive failure behavior. Capability of TOUGH-RFPA simulator was verified by two benchmark problems related to hydraulic-mechanical and thermal-mechanical response in a hollow cylinder. Then the applicability of TOUGH-RFPA was demonstrated for modeling cooling-induced damage around a very deep geothermal well with mesoscopic heterogeneous properties, which showed complex damage evolution and fracturing during well cooling and fluid pressurization.

Experimental studies show that compacted bentonite used as a backfill material for nuclear waste repository experiences strong coupling between chemical and mechanical processes. In this paper, we use a dual-structure expansive soil model, referred to as the Barcelona Expansive Model (BExM), to predict the behavior of bentonite buffer in subsurface emplacement tunnels for high-level radioactive waste. After emplacement, the bentonite is subjected to complex and coupled Thermal-Hydraulic-Mechanical-Chemical (THMC) processes. The BExM constitutive model is implemented in a multi-phase reactive transport and geomechanics simulator, TOUGHREACT-FLAC3D, and the model is verified with one benchmark test on FEBEX bentonite. We utilize a one-way Chemo-Mechanical (C-M) coupling approach, in which chemical changes affect the mechanical behavior of bentonite through the BExM linked with the evolution of mass fraction of smectite, exchangeable cation concentration, and ionic strength via osmotic suction. The parameters of the new coupled model for FEBEX bentonite are calibrated against a series of laboratory experiments with various salinity solutions. Finally, coupled THMC modeling is conducted for a generic argillite repository with bentonite buffer under high temperature, focusing on the long-term chemical change and its effect on the mechanical process.

Fault reactivation due to water injection is assessed within the scope of nuclear waste disposal design. A model using the distinct element method is applied to reproduce the fault reactivation during an experiment carried out at the Mont Terri Underground Research Laboratory in Switzerland. A conceptual model is first presented to understand the hydro-mechanical coupling behavior between water pressure and rock joint movement. The model simulations show that the dominant factors on fault slip are shear stress and frictional resistance. Moreover, modeling shows that fault reversible opening in the normal direction occurs first at a lower pressure, whereas shear displacement as a result of shear slip is produced once a sufficiently high pressure is reached. We demonstrate that a coupled numerical analysis of the rock displacement trend and fluid pressure measured at the injection point allow an in-situ estimation of the principal stresses.

This study presents comprehensive coupled hydro-mechanical numerical modeling of the first and second hydraulic stimulations at the Pohang enhanced geothermal system (EGS) site in order to improve the understanding on the key stimulation mechanisms in the fractured reservoir. Two models for PX-2 and PX-1 wells were developed including fracture zones of different orientations and permeability. Shear dilation and frictional plastic strain-softening were implemented in a non-linear stress dependent fracture model. A history matching of the wellhead pressure curves was carried out in the early days of each stimulation. The PX-2 model with a fracture jacking mechanism achieved the reproduction of highly reversible non-linear aperture changes observed in situ. The PX-1 model successfully simulated a wellhead pressure peak and drop associated with shear slip and dilation by applying frictional plastic strain softening. Simulated critical wellhead pressure for shear slip near the PX-1 well was greater than the prediction by a simple slip potential analysis due to local stress changes by the poroelastic effect. The numerical modeling confirmed that the combination of shear dilation and jacking was adequate to capture the pressure evolution during the increasing step-rate test at PX-1. Besides, spatio-temporal changes in pressure, total stress, and permeability evaluated in the fracture zones greatly enhanced the understanding on the coupled behavior caused by the hydraulic stimulations. The possibility of zones with lower permeability away from the PX-1 well was suggested by an alternative model. The current numerical study demonstrates that the key hydro-mechanical processes of shear slip and dilation and hydraulic jacking observed in the fractured reservoir can be successfully reproduced.

Geologic disposal is a promising solution for a safe permanent isolation of accumulated high-level nuclear waste from nuclear powerplants. The behavior of host rock is highly coupled thermally, hydromechanically, and chemically. Numerical simulations of such coupled phenomena for the extremely long term (>100,000 years) and large length scale (>1 km) of geologic disposal remain to be computationally challenging. In this study, a methodology has been developed to approximate stress and stress-induced permeability change in host rock using only thermo-hydrological (TH) variables. A coupled thermo-hydromechanical (THM) simulation is carried out using TOUGH-FLAC simulator to model THM behaviors of a generic nuclear waste repository, in order to evaluate the performance of the developed methodology, which is implemented in a coupled TH simulation. Results show that stress and permeability change estimated by the developed methodology in the TH-coupled simulation match those calculated in the THM-coupled simulation over the simulated timeframe of over 10,000 years. Details about errors in stress and permeability estimates accrued by the developed methodology are discussed in this paper. The developed methodology will help incorporate stress-induced permeability change into existing TH simulators for the long-term radionuclide transport in geologic disposal.

A coupled multiphase fluid flow and discrete fracturing model is applied to simulate bench-scale gas migration experiments on compacted bentonite. The numerical modeling is based on the linking of the multiphase fluid flow simulator TOUGH2 with a Rigid-Body-Spring Network model, which enables a discrete (lattice) representation of elasticity and individual fractures. The evolution of a complex network of dilatant flow paths is modeled through opening and breakage of lattice interface bonds between porous-elastic matrix elements. To achieve a good match with the experimental results, including an abrupt gas breakthrough along with pressure and stress responses, it was necessary to calibrate model parameters for (1) air-entry pressure, (2) shear and tensile failure of lattice interface bonds, (3) moisture swelling/shrinkage effects on stress, and (4) aperture-dependent permeability of dilatant flow paths. Our best-fit conceptual model considers a pervasive network of discrete flow paths propagating from the gas injection point, whereas some of the experimental data indicate the potential for heterogeneous and unstable flow paths.

Geodetic observations, often in conjunction with other data, provide a cost-effective means for identifying and characterizing geothermal resources. A review of the various methods reveals how the technology for measuring deformation has advanced considerably in the past few decades. Currently, interferometric synthetic aperture radar is the method of choice for monitoring deformation at a geothermal field. A discussion of geodetic monitoring at The Geysers geothermal field, California, illustrates some of the progress made and the challenges that remain.

We present results from an international model comparison study involving a series of controlled fault activation experiments in Opalinus Clay at the Mont Terri Laboratory, Switzerland. The fault activation experiments were conducted in situ by water injection at variable pressure from boreholes targeting different parts of the Main Fault crossing several tunnels and galleries of the Mont Terri Laboratory. The model simulations focused on (1) an experiment activating a discontinuity (fracture or minor fault) within the damage zone of the main fault, and (2) an experiment activating discontinuities close to the core of the Main Fault. The experimental data consist of coupled hydraulic and mechanical responses monitored at an injection borehole and a monitoring borehole located a few meters away. After overcoming several modeling issues along with necessary model developments, a reasonably good agreement was achieved between the modeling results and the field observations. The fault activation experiments displayed an abrupt flow rate increase associated with a sudden fracture opening and rupture propagation after the injection pressure reached above the estimated normal stress on the fracture. This was followed by an abrupt flow rate decrease, indicating hydraulic closing, once the injection pressure decreases to an estimated 1–2 MPa below the stress normal to the opened fracture. The models were able to capture this abrupt hydromechanical behavior, including an observed dominant opening behavior along with the rupture propagation, while the spatial extent of the shear rupture and the quantity of peak injection flow were the most challenging to predict.

We present coupled Thermo–Hydro–Mechanical (THM) modeling of geologic nuclear waste disposal in argillaceous claystone, focusing on thermally-induced pressure changes and the potential for such pressure changes to induce hydro-fracturing. To investigate this possibility, we first conduct a three-dimensional repository scale model, with host rock properties, repository design, and nuclear waste decay heat functions derived from the French concept of geologic disposal in argillaceous claystone. The model simulations show that the highest potential for hydro-fracturing occurs between emplacement micro-tunnels (cells) at the center rather than at the edge of the repository. We further investigate the use of a two-dimensional single cell model as a simpler surrogate for a full three-dimensional model. Our results reveal that such geometric simplification is reasonably accurate for modeling coupled THM processes at the center of the repository, though it overestimates the likelihood for hydro-fracturing and substantially overpredicts ground surface uplift. A parameter study shows the importance of adjacent higher permeability geological layers that play a significant role in dissipating overpressure in the host rock layer. Finally, the study shows the importance of the spacing between emplacement tunnels, which if too short results in a higher fluid pressure and a strongly increased potential for hydro-fracturing. Overall, the study suggests, the limiting factor in the thermal management and design of a repository is not necessarily the maximum temperature in the engineered barrier system near the waste packages, but rather the more modest host rock temperature between emplacement tunnels, due to the potential for thermal damage.

Micromechanical modeling of geomaterials is challenging because of the complex geometry of discontinuities and potentially large number of deformable material bodies that contact each other dynamically. In this study, we have developed a numerical approach for micromechanical analysis of deformable geomaterials with dynamic contacts. In our approach, we detect contacts among multiple blocks with arbitrary shapes, enforce different contact constraints for three different contact states of separated, bonded, and sliding, and iterate within each time step to ensure convergence of contact states. With these features, we are able to simulate the dynamic contact evolution at the microscale for realistic geomaterials having arbitrary shapes of grains and interfaces. We demonstrate the capability with several examples, including a rough fracture with different geometric surface asperity characteristics, settling of clay aggregates, compaction of a loosely packed sand, and failure of an intact marble sample. With our model, we are able to accurately analyze (1) large displacements and/or deformation, (2) the process of high stress accumulated at contact areas, (3) the failure of a mineral cemented rock samples under high stress, and (4) post-failure fragmentation. The analysis highlights the importance of accurately capturing (1) the sequential evolution of geomaterials responding to stress as motion, deformation, and high stress; (2) large geometric features outside the norms (such as large asperities and sharp corners) as such features can dominate the micromechanical behavior; and (3) different mechanical behavior between loosely packed and tightly packed granular systems.

We apply thermoporoelasticity and a sequentially coupling technique for modeling thermally-driven coupled Thermo-Hydro-Mechanical (THM) processes in tight claystone. A THM benchmark case with a corresponding analytic solution for thermoporoelasticity under a constant heat loading verifies the model. Thereafter, two in situ heating experiments are simulated for model validation: a smaller-scale heating experiment (TED experiment) and a larger-scale experiment (ALC experiment) in Callovo-Oxfordian (COx) claystone at the Meuse/Haute-Marne underground research laboratory in France. The model exhibits good performance to match the observed temperature and pore pressure evolution for the smaller-scale TED experiment. For the larger-scale ALC experiment, general trends of thermal-pressurization are captured in the modeling, but pressure is underestimated at some monitoring points during cool-down. This indicates that the THM response in the field may be affected by the variability of rock's properties or irreversible or time-dependent mechanical processes that are not included in the current thermoporoelastic model. The main contributions of this work are as follows: (1) we verify and validate the numerical simulator, TOUGH-FLAC, to be a valuable coupled THM modeling tool; (2) prove that the laboratory determined material parameters can be used as reference values for upscaling experiments. However, to better identify and quantify THM processes with modeling of in situ tests, more emphasize should be dedicated to obtaining high-quality mechanical deformation data.

In this study, the combined effects of geometrical distribution and geomechanical deformation of fracture networks on fluid flow through fractured geological media are investigated numerically. We consider a finite-sized model domain in which the geometry of fracture systems follows a power-law length scaling. The geomechanical response of the fractured rock is simulated using a hybrid finite-discrete element model, which can capture the deformation of intact rocks, the interaction of matrix blocks, the displacement of discrete fractures and the propagation of new cracks. Under far-field stress loading, the locally variable stress distribution in the fractured rock leads to a stress-dependent variable aperture field controlled by compression-induced closure and shear-induced dilatancy of rough fractures. The equivalent permeability of the deformed fractured rock is calculated by solving for the fracture-matrix flow considering the cubic relationship between fracture aperture and flow rate at each local fracture segment. We report that the geometrical connectivity of fracture networks plays a critical role in the hydromechanical processes in fractured rocks. A well-connected fracture system under a high stress ratio condition exhibits intense frictional sliding and large fracture dilation/opening, leading to greater rock mass permeability. However, a disconnected fracture network accommodates much less fracture shearing and opening, and has much lower bulk permeability. We further propose an analytical solution for the relationship between the equivalent permeability of fractured rocks and the connectivity metric (i.e. percolation parameter) of fracture networks, which yields an excellent match to the numerical results. We infer that fluid flow through a well-connected system is governed by traversing channels (forming an “in parallel” architecture) and thus equivalent permeability is sensitive to stress loading (due to stress-dependent fracture permeability), whilst fluid flow through a disconnected system is more ruled by matrix (linking isolated clusters “in series”) and has much less stress dependency.

The greatest challenges of rigorously modeling coupled hydro-mechanical (HM) processes in fractured geological media at different scales are associated with computational geometry. These challenges include dynamic shearing and opening of intersecting fractures at discrete fracture scales as a result of coupled processes, and contact alteration along rough fracture surfaces that triggers structural and physical changes of fractures at micro-asperity scale. In this paper, these challenges are tackled by developing a comprehensive modeling approach for coupled processes in fractured geological media based on numerical manifold method (NMM) at multiple scales. Based on their distinct geometric features, fractures are categorized into three different scales: dominant fracture, discrete fracture, and discontinuum asperity scales. Here the scale is relative, that of the fracture relative to that of the research interest or domain. Different geometric representations of fractures at different scales are used, and different governing equations and constitutive relationships are applied. For dominant fractures, a finite thickness zone model is developed to treat a fracture as a porous nonlinear domain. Nonlinear fracture mechanical behavior is accurately modeled with an implicit approach based on strain energy. For discrete fractures, a zero-dimensional model was developed for analyzing fluid flow and mechanics in fractures that are geometrically treated as boundaries of the rock matrix. With the zero-dimensional model, these fractures can be modeled with arbitrary orientations and intersections. They can be fluid conduits or seals, and can be open, bonded or sliding. For the discontinuum asperity scale, the geometry of rough fracture surfaces is explicitly represented and contacts involving dynamic alteration of contacts among asperities are rigorously calculated. Using this approach, fracture alteration caused by deformation, re-arrangement and sliding of rough surfaces can be captured. Our comprehensive model is able to handle the computational challenges with accurate representation of intersections and shearing of fractures at the discrete fracture scale and rigorously treats contacts along rough fracture surfaces at the discontinuum asperity scale. With future development of three-dimensional (3D) geometric representation of discrete fracture networks in porous rock and contacts among multi-body systems, this model is promising as a basis of 3D fully coupled analysis of fractures at multiple scales, for advancing understanding and optimizing energy recovery and storage in fractured geological media.

Coupled thermal-hydro-mechanical (THM) numerical modeling is conducted to study thermal management associated with geologic disposal of spent nuclear fuel (SNF) in large dual-purpose canisters (DPCs). DPCs, are containers designed for SNF storage and transportation and if determined to be feasible for permanent geological disposal could provide a cost effective disposal solution. However, one of the challenges to direct disposal of DPCs is thermal management to avoid overheating of the Engineered Barrier System (EBS), including bentonite-backfill used as a protective buffer. The model simulations show that the use of a backfill that is thermally engineered for high heat conduction can reduce the EBS temperature to acceptable levels for disposal of large waste canisters in backfilled tunnels. On the other hand, the use of high heat conduction backfill will not reduce the far field rock peak temperature that can occur several thousand years after closure of the repository. This longer term host rock peak temperature generates thermal-poro-elastic stress and geomechanical changes that must be considered in the thermal management and design of a repository.