Engineering topics
Guglielmi, Yves
Publications and source records attributed to Guglielmi, Yves.
Hydraulic Modeling of Induced and Propagated Fractures: Analysis of Flow and Pressure Data From Hydromechanical Experiments in the COSC-1 Deep Borehole in Crystalline Rock Near Åre, Sweden
To characterize the coupled hydromechanical behavior of rock fractures, the step-rate injection method for fracture in-situ properties (SIMFIP) was conducted with a specialized downhole probe developed by Guglielmi et.al. In June 2019, a field campaign was carried out near Åre, Sweden, where the SIMFIP probe was applied in the Collisional Orogeny in the Scandinavian Caledonides-1 scientific borehole to understand the dynamics of injection-induced fracture initiation, fracture opening, and shearing due to water injection-withdrawal in a borehole interval isolated by two packers. Three intervals were investigated at ~500 m depth: (a) an unfractured section (intact rock), (b) a section with non-conductive fractures, and (c) a section with hydraulically conductive fractures. Pressure, injection flow rate, and borehole wall displacement were simultaneously measured during the tests. In the present study, the geometry of the induced fracture and deformation of existing fractures at different time stages of the tests are determined based on a hydrologic model by using the measured pressure and flow data during each time stage of the experiment. A numerical model for the fluid flow within the fracture and the packed-off borehole interval is implemented within COMSOL Multiphysics. By matching model simulations with observed data for all three sections, estimates of the induced and propagated fractures' radius and aperture at successive time stages have been obtained in each case. We could also determine the non-linear relationship between fracture aperture and pressure for values above fracture opening pressures. The model results provide insights for the understanding of pressure-induced fracture initiation and propagation in crystalline rock.
Recommended Practices for Managing Induced Seismicity Risk Associated with Geologic Carbon Storage
The geologic storage of carbon dioxide (CO 2 ) is one method to help reduce or eliminate atmospheric CO 2 emissions. The sequestered CO 2 is originally captured from the atmosphere or from a stationary industrial source and subsequently injected into a deep subsurface porous rock formation. To facilitate the successful deployment of field scale carbon storage projects, the U.S. Department of Energy (DOE) is developing tools and protocols for defensible, science-based frameworks to quantify and mitigate risks associated with the long-term storage of CO 2 . This protocol specifically addresses the risk of induced seismicity due to injection in a geologic carbon storage (GCS) site. This integrated and risk-based protocol is a product of the U.S. DOE Fossil Energy’s National Risk Assessment Partnership (NRAP), a multi-year collaborative research effort of Los Alamos National Laboratory (LANL), Lawrence Berkeley National Laboratory (LBNL), Lawrence Livermore National Laboratory (LLNL), National Energy Technology Laboratory (NETL), and Pacific Northwest National Laboratory (PNNL). These recommended practices describe a set of 7 steps to evaluate, manage, communicate, and mitigate the risk of induced seismicity at GCS sites. The base methodology of the recommended practices follows a framework similar to the Protocol for Addressing Induced Seismicity Associated with Enhanced Geothermal Systems (Majer et al., 2012), developed for the Geothermal Technology Office of the U.S. DOE. These recommended practices present a framework to systematically assess the induced seismicity risk and quantify the associated uncertainties. These recommendations are based on current research and are sufficiently general to allow for modification and application to a variety of different types of sites. The substance of the recommended practices contained herein includes both technical and non-technical issues, and covers all operational stages of the GCS project lifecycle. They start at the preliminary risk assessment phase, continue through site assessment and characterization, include best practice communication and seismic monitoring plan methodologies, discuss the evaluation and mitigation of seismic hazard and risk, and closes with an exploration of operational management plans, which conclude when the induced seismicity risk abates back to background level. The focus of these recommendations is on actively managing the risks associated with induced seismicity by developing an actionable risk management plan that starts at the project proposal stage and continues through site closure through an iterative assessment and improvement process. The audience of this document is expected to include all interested stakeholders (e.g., operators, project developers, regulators, and the general public) and is expressly written to be accessible to this broad range of partners. This document is intended to disseminate knowledge gained through recent advances in the science of induced seismicity hazard and risk assessments, to provide updates based on recent experience gained by similar corollary injection-induced seismicity cases, and most importantly to establish a uniform framework to carry out a successful induced seismicity risk management plan for carbon storage projects in the future. These recommendations do not directly address any domestic or international regulations or standards. A complementary NRAP report makes recommendations for the assessment and management of environmental subsurface risks associated with unwanted fluid migration at GCS sites (Thomas et al., 2021) and should be referred to in order to address those additional GCS site risks.
Field-scale fault reactivation experiments by fluid injection highlight aseismic leakage in caprock analogs: Implications for CO 2 sequestration
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 CO 2 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 (M w < -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.
Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies (M2SF-21SN010301072)
The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress in modeling and experimental approaches towards the characterization of chemical and physical phenomena that could impact the long-term safety assessment of heatgenerating nuclear waste disposition in deep-seated clay/shale/argillaceous rock. International collaboration activities such as heater tests, continuous field data monitoring, and postmortem analysis of samples recovered from these have elucidated key information regarding changes in the engineered barrier system (EBS) material exposed to years of thermal loads. Chemical and structural analyses of sampled bentonite material from such tests as well as experiments conducted on these are key to the characterization of thermal effects affecting bentonite clay barrier performance and the extent of sacrificial zones in the EBS during the thermal period. Thermal, hydrologic, and chemical data collected from heater tests and laboratory experiments has been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches (e.g., continuum and discrete) to tackle issues related to flow and transport at various scales of the host-rock, its interactions with barrier materials, and EBS design concept.
Close Observation of Hydraulic Fracturing at EGS Collab Experiment 1: Fracture Trajectory, Microseismic Interpretations, and the Role of Natural Fractures
Abstract Despite the wide application of hydraulic fracturing in enhanced geothermal system (EGS) development and unconventional hydrocarbon production, our understanding of the many factors affecting the propagation of hydraulic fractures has relied on circumstantial evidence. This is partly due to the scarcity of direct observations in the subsurface. The EGS Collab project attempts to address these issues in the context of EGS research by performing intermediate‐scale (∼10 m) hydraulic stimulation experiments in a thoroughly characterized and heavily instrumented underground testbed. This paper analyzes the data collected from the first suite of hydraulic fracturing tests in this testbed, consisting of seven stimulation episodes. High‐quality microseismic data delineate five planar features very clearly. Combining fracture‐wellbore intersection observations from distributed temperature sensing (DTS) and visual observations from an open‐hole well, as well as prior in‐situ stress measurements, we conclude with high certainty that the four larger planes were hydraulic fractures. The growth of the hydraulic fractures was temporarily halted by a prominent, open natural fracture in the testbed but they eventually crossed and slightly reoriented under continued stimulation. Mineral‐filled (i.e., healed) natural fractures, though prevalent in the testbed, did not have perceivable effects on hydraulic fracture propagation. The high‐quality, mutually corroborating data sets allowed conclusions to be drawn with high confidence and attests to the advantage of intermediate‐scale experiments in subsurface research.
In-situ observation of pre-, co- and post-seismic shear slip at 1.5 km depth
Understanding the initiation and arrest of earthquakes is one of the long-standing challenges of seismology. Here we report on direct observations of borehole displacement by a meter-sized shear rupture induced by pressurization of metamorphic rock at 1.5 km depth. We observed the acceleration of sliding, followed by fast co-seismic slip and transient afterslip. Total displacements were about 7, 5.5 and 9.5 micrometers, respectively for the observed pre-slip, co-seismic slip and afterslip. The observed pre-slip lasted about 0.4 seconds. Co-seismic slip was recorded by the 1 kHz displacement recording and a 12-component array of 3-C accelerometers sampled at 100 kHz. The observed afterslip is consistent with analytical models of arrest in a velocity-strengthening region and subsequent stress relaxation. < p >The observed slip vector agrees with the activation of a bedding plane within the phyllite, which is corroborated by relocated seismic events that were observed during the later stages of the injection experiment.
Modelling the water injection induced fault slip and its application to in-situ stress estimation
Not Available
An international model comparison study of controlled fault activation experiments in argillaceous claystone at the Mont Terri Laboratory
Not Available
Hydraulic fracturing experiments at 1500 m depth in a deep mine: Highlights from the kISMET project
In support of the U.S. DOE SubTER Crosscut initiative, we established a field test facility in a deep mine and designed and carried out in situ hydraulic fracturing experiments relevant to enhanced geothermal systems (EGS) in crystalline rock to characterize the stress field, understand the effects of rock fabric on fracturing, and gain experience in monitoring using geophysical methods. The project also included 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 Sanford Underground Research Facility (SURF) 4757 ft (1450 m) below ground (on the 4850 ft level (4850L)) in phyllite of the Precambrian Poorman Formation. We drilled and continuously cored five near-vertical boreholes in a line on 3 m (10 ft) spacing, deviating the two outermost 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 velocity heterogeneity along each core indicating strong, fine-scale (~1 cm or smaller) changes in the mechanical properties of the rock. Field measurements of the stress field by hydraulic fracturing showed that the minimum horizontal stress at the kISMET site averages 21.7 MPa (3146 psi) trending approximately N-S (356 degrees azimuth) and plunging slightly NNW at 12°. The vertical and horizontal maximum stresses are similar in magnitude at 42-44 MPa (6090-6380 psi) for the depths of testing, which averaged approximately 1530 m (5030 ft). Hydraulic fractures were remarkably uniform suggesting core-scale and larger rock fabric did not play a role in controlling fracture orientation. Analytical solutions suggest that the fracture radius of the large fracture (stimulation test) was more than 6 m (20 ft), depending on the unknown amount of leak-off.
Decovalex-2019 (Executive Summary)
The DECOVALEX Project is an on-going international research collaboration, established in 1992, to advance the understanding and modeling of coupled Thermal (T), Hydrological (H), Mechanical (M) and Chemical (C) processes in geological in geological systems. DECOVALEX was initially motivated by the recognition that prediction of these coupled effects is an essential part of the performance and safety assessment of geologic disposal systems for radioactive waste and spent nuclear fuel. Later it was realized that these processes also play a critical role in other subsurface engineering activities, such as subsurface CO 2 storage, enhanced geothermal systems, and unconventional oil and gas production through hydraulic fracturing. Research teams from many countries (e.g., Canada, China, Czech Republic, Finland, France, Germany, Japan, Republic of Korea, Spain, Sweden, Switzerland, Taiwan, United Kingdom, and the United States) various institutions have participated in the DECOVALEX Project over the years, providing a wide range of perspectives and solutions to these complex problems. These institutions represent radioactive waste management organizations, national research institutes, regulatory agencies, universities, as well as industry and consulting groups. The overall aim of DECOVALEX-2019 was to increase the understanding of various thermo-hydro-mechanical-chemical processes of importance for radionuclide release and transport from a repository to the biosphere and how they can be described and modelled using mathematical models. The scientific and technical objectives are: to increase the basic understanding of T-H-M-C coupled processes in fractured rocks (crystalline, sedimentary, argillaceous) and buffer materials; to investigate the predictive capabilities of different codes to field experiments and to perform verification of codes; to exchange experimental data, and improve the understanding of the constitutive behavior of crystalline and argillaceous rock masses and buffer materials; and to perform THMC calculations in a performance/safety assessment context.
DECOVALEX-2019 (Task B Final Report)
The DECOVALEX Project is an on-going international research collaboration, established in 1992, to advance the understanding and modeling of coupled Thermal (T), Hydrological (H), Mechanical (M) and Chemical (C) processes in geological in geological systems. DECOVALEX was initially motivated by the recognition that prediction of these coupled effects is an essential part of the performance and safety assessment of geologic disposal systems for radioactive waste and spent nuclear fuel. This document is the final report of Task B which was proposed and coordinated by the Swiss Federal Nuclear Safety Inspectorate (ENSI) and Lawrence Berkeley National Laboratory (LBNL), presenting the technical definitions of the problems studied, approaches applied, achievements made and outstanding issues for future research. Task B focused mainly on modelling of fault activation experiments (FS experiment) recently performed at the Mont Terri underground research laboratory, Switzerland. The FS experiment explores the coupling between fault reactivation in a clay host rock and the potential-enhanced fluid displacement through a previously low-permeability formation. Based on this experiment, the DECOVALEX-2019 Task B addresses how the change in permeability induced by the fault activation and the resulting fluid flow within the fault can be simulated including the dependence on its mechanical behavior. This will support the understanding of the processes during fault activation itself and it can help to determine consequences for contaminant transport due to potentially created permeable flow paths in otherwise low permeability argillaceous rocks.
Slip Tendency Analysis of Fracture Networks to Determine Suitability of Candidate Testbeds for the EGS Collab Hydroshear Experiment
Experiment 2 of the EGS Collab project is aimed at testing stimulation by hydro-shearing of existing natural fractures, versus Experiment 1, which is focused on hydraulic fracturing a rock mass to enhance permeability. The main criterion for the testbed selection in Experiment 2 is the presence of an interconnected network of fractures, at least 10 meters in extent, and with orientations such that shear slip can be induced at injection pressures less than the minimum horizontal stress (Shmin). The feasibility analysis for this experiment requires a well-constrained stress state along with well-characterized fracture networks. The fracture systems at two candidate locations in the Sanford Underground Research Facility (SURF), the 4850 Level (number refers to depth below ground surface in feet) and the 4100 Level, have been characterized to different extents. The 4850 Level has a well-characterized Discrete Fracture Network (DFN) from borehole (drilled for an experimental test bed) and drift observations. In contrast, the 4100 level has a fracture network characterized only by observations from the drift wall, which provides little constraint on fracture extents. This paper will present assessments for the interpreted slip potential for the natural fractures at the two locations. Data uncertainties are addressed by performing a probabilistic analysis that takes into account the uncertainty in the stress state, uncertainty in the fracture properties, and preliminary borehole locations. There is strong evidence of shear stimulation of some natural fractures in Experiment 1 on the 4850 level. We use examples of natural fractures strongly linked to shear stimulation in Experiment 1 from multiple monitoring indicators to test the consistency of the stress model and guide the uncertainty interpretation. This analysis is intended to guide the site selection process for Experiment 2 by highlighting the fracture orientations that are likely to be shear-stimulated in a majority of modeled realizations.
EGS Collab Experiment 1: SIMFIP Notch-164 GRL Paper
Characterizing the stimulation mode of a fracture is critical to assess the hydraulic efficiency and the seismic risk related to deep fluid manipulations. We have monitored the three-dimensional displacements of a fluid-driven fracture during water injections in a borehole at ~1.5 km depth in the crystalline rock of the Sanford Underground Research Facility (USA). The fracture initiates at 61% of the minimum horizontal stress by micro-shearing of the borehole on a foliation plane. As the fluid pressure increases further, borehole axial and radial displacements increase with injection time highlighting the opening and sliding of a new hydrofracture growing ~10 m away from the borehole, in accordance with the ambient normal stress regime and in alignment with the microseismicity. Our study reveals how fluid-driven fracture stimulation can be facilitated by a mixed-mode process controlled by the complex hydromechanical evolution of the growing fracture. The data presented in this submission refer to the SIMFIP measurements and analyses of the stimulation tests conducted on the 164 ft (50 m) notch of the Sanford Underground Research Facility (SURF), during the EGS-Collab test 1. In addition to the datafiles, there is the draft of a manuscript submitted to Geophysical Research Letters (GRL).
Crystalline Disposal R&D at LBNL: FY20 Progress Report
Within the Spent Fuel & Waste Science and Technology (SFWST) Program, research work continues further with the goal to better understand long-term performance of disposal systems in three main geologic rock types: clay/shale, salt, and crystalline rock. This report documents LBNL’s research activities related to investigations of crystalline host rock according to the scope of two work packages: SF-20LB01030207 “Crystalline International Collaborations – LBNL,” and SF-20LB01030203 “Crystalline Disposal R&D – LBNL.” These research activities correspond are related to key Features, Events, and Processes (FEPs), ranked in importance from medium to high, as listed in Table 7 of the Used Fuel Disposition Campaign Disposal Research and Development Roadmap (FCR&D-USED-2011-000065 REV0) (Nutt, 2011). Specifically, these research activities address FEP 2.2.01, Excavation Disturbed Zone (EZD). The results of these research activities provide important insights into understanding and predicting flow and transport processes that could occur in low-permeability crystalline rocks, in which fractures might serve as main conduits for fluid flow and radionuclide transport. The evolution of the EDZ during the excavation of the tunnel as well as the evolution of microcrack growth within the EDZ after the emplacement of backfill are critical for predicting the long- term behavior of the EDZ. A number of factors including stress, temperature, water activity, capillary pressure, chemistry, and mineralogy can affect the rock evolution, which require advanced experimental tools to study it.
Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies
The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress in modeling and experimental approaches towards the characterization of chemical and physical phenomena that could impact the long-term safety assessment of heat-generating nuclear waste disposition in deep clay/shale/argillaceous rock. International collaboration activities such as heater tests and postmortem analysis of samples recovered from these have elucidated key information regarding changes in the engineered barrier system (EBS) material exposed to years of thermal loads. Chemical and structural analyses of sampled bentonite material from such tests has as well as experiments conducted on these are key to the characterization of thermal effects affecting bentonite clay barrier performance and the extent of sacrificial zones in the EBS during the thermal period. Thermal, hydrologic, and chemical data collected from heater tests and laboratory experiments has been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches (e.g., continuum vs. discrete) to tackle issues related to flow and transport at various scales of the host-rock and EBS design concept. Consideration of direct disposal of large capacity dual-purpose canisters (DPCs) as part of the back-end SNF waste disposition strategy has generated interest in improving our understanding of the effects of elevated temperatures on the EBS design. This is particularly important for backfilled repository concepts where temperature plays a key role in the EBS behavior and long-term performance. This report describes multiple R&D efforts on disposal in argillaceous geologic media through development and application of coupled THMC process models, experimental studies on clay/metal/cement barrier and host-rock (argillite) material interactions, molecular dynamic (MD) simulations of water transport during (swelling) clay dehydration, first-principles studies of metaschoepite (UO 2 corrosion product) stability, and advances in thermodynamic plus surface complexation database development. Drift-scale URL experiments provides key data for testing hydrological-chemical (HC) model involving strong couplings of fluid mixing and barrier material chemical interactions. The THM modeling focuses on heater test experiments in argillite rock and gas migration in bentonite as part of international collaboration activities at underground research laboratories (URLs). In addition, field testing at an URL involves in situ analysis of fault slip behavior and fault permeability. Pore-scale modeling of gas bubble migration is also being investigated within the gas migration modeling effort. Interaction experiments on bentonite samples from heater test under ambient and elevated temperatures permit the evaluation of ion exchange, phase stability, and mineral transformation changes that could impact clay swelling. Advances in the development, testing, and implementation of a spent nuclear fuel (SNF) degradation model coupled with canister corrosion focus on the effects of hydrogen gas generation and its integration with Geologic Disposal Safety Assessment (GDSA). GDSA integration activities includes evaluation of groundwater chemistries in shale formations.
Modeling of fluid injection-induced fault reactivation using coupled fluid flow and mechanical interface model
The present study is aimed at developing a numerical model to reproduce coupled hydro-mechanical processes associated with fault reactivation by fluid injection in low permeability rock, as part of the DECOVALEX-2019 project Task B. We proposed a modeling approach for simulating the processes using the TOUGH-FLAC simulator, and modeled a fault reactivation experiment conducted at Mont Terri Rock Laboratory in Switzerland. The first step of the study involved benchmark calculations considering a simplified fault plane and geometry. Fluid flow along a fault was modeled using elements of aperture-sized thickness on the basis of Darcy's law and the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by zero-thickness interface elements in FLAC3D upon which a slip and/or separation is allowed. A methodology to connect a TOUGH2 volume element to a FLAC3D interface element was developed for handling the hydro-mechanical interactions on the fault during fluid injection. Two different fault models for describing the evolutions of hydraulic aperture by elastic fracture opening and failure-induced aperture increase were considered in the benchmark calculations. In the coupling process, the changes in geometrical features and hydrological properties induced by mechanical deformation were continuously updated. The transient responses of the fault and host rock to stepwise pressurization were examined during the simulation. The hydro-mechanical behavior, including the injection flow rate, pressure distribution around the borehole, stress conditions, and displacements in normal and shear directions were monitored in the surrounding rock and along the fault. The results of benchmark calculations suggest that the developed model reasonably represents the hydro-mechanical behavior of a fault and the surrounding rock. This modeling approach was applied to the fault reactivation experiment of the Mont Terri Rock Laboratory. In this interpretive modeling, a parametric study was conducted to examine the effects of input parameters regarding in situ stress and fault properties on the hydro-mechanical responses of the fault to water injection. Then, an optimal parameter set to reproduce the field experiment results was chosen by trial-and-error. The injection flow rate and pressure response during fault reactivation closely matched those obtained at the site, which indicates the capability of the model to appropriately capture the progressive pathway evolution during fault reactivation tests at the site. Finally, the anchor displacements were overestimated by the model, but a fair agreement was obtained in terms of the order of magnitude and the variation tendency.
Estimating perturbed stress from 3-D borehole displacements induced by fluid injection in fractured or faulted shales
SUMMARY Hydrofracturing stress measurements in fractured and anisotropic shales are notoriously difficult, because opening of existing geological features tends to prevent the creation of a pure hydraulic fracture perpendicular to the least compressive principal stress. Here we show how adding 3-D borehole-displacement measurements while conducting the hydraulic injection test helps to better constrain the principal stress orientations and magnitudes. We developed a 3-D fully coupled hydromechanical numerical model to analyse the displacement, fluid pressure and injection flow-rate data measured during an injection pressure-step-rate test conducted to activate a faulted borehole interval in the Mont Terri Opalinus Clay (Switzerland). We find that injected fluids can only penetrate the fault when it is at or above the Coulomb failure pressure. Borehole displacement orientations are sensitive to a ∼15° variation in the stress–tensor orientation and a 1 MPa stress magnitude variation. Although some dispersion occurs while rupture is propagating along the fault plane ∼4 m away from the borehole, the maximum density of displacement orientations consistently informs about the stress orientation. Thus, an extended injection step-rate approach coupled with an accurate in situ measurement of the borehole wall displacements can be used to better constrain the local stress field perturbations in fractured shales and in heterogeneous rock in general.