Lattice modeling of excavation damage in argillaceous clay formations: Influence of deformation and strength anisotropy

This paper presents modeling of mechanical anisotropy in argillaceous rocks using an irregular lattice modeling approach, namely the rigid-body-spring network. To represent the mechanical anisotropy, new schemes are implemented in the modeling framework. The directionality of elastic deformation is resolved by modifying the element formulation with anisotropic elastic properties. The anisotropy of strength and failure characteristics is facilitated by adopting orientation-dependent failure criteria into the failure model. The verification of the improved modeling procedures is performed against theoretical model predictions for unconfined compression tests with various bedding orientations. Furthermore, excavation damage and fracturing processes in rock formations are simulated for different geomechanical configurations, such as rock anisotropy and tectonic heterogeneity. The simulated excavation damage characteristics are realistic and comparable with the actual field observation at a tunnel located in an argillaceous clay formation. The simulation results provide insights into the excavation damage zone phenomena with an explicit representation of fracturing processes

Development of a Numerical Simulator for Analyzing the Geomechanical Performance of Hydrate-Bearing Sediments

In this paper, we describe the development and application of a numerical simulator that analyzes the geomechanical performance of hydrate-bearing sediments, which may become an important future energy supply. The simulator is developed by coupling a robust numerical simulator of coupled fluid flow, hydrate thermodynamics, and phase behavior in geologic media (TOUGH+HYDRATE) with an established geomechanical code (FLAC3D). We demonstrate the current simulator capabilities and applicability for two examples of geomechanical responses of hydrate bearing sediments during production-induced hydrate dissociation. In these applications, the coupled geomechanical behavior within hydrate-bearing seducements are considered through a Mohr-Coulomb constitutive model, corrected for changes in pore-filling hydrate and ice content, based on laboratory data. The results demonstrate how depressurization-based gas production from oceanic hydrate deposits may lead to severe geomechanical problems unless care is taken in designing the production scheme. We conclude that the coupled simulator can be used to design production strategies for optimizing production, while avoiding damaging geomechanical problems

Analysis of Injection-Induced Micro-Earthquakes in a Geothermal Steam Reservoir, The Geysers Geothermal Field, California

In this study we analyze relative contributions to the cause and mechanism of injection-induced micro-earthquakes (MEQs) at The Geysers geothermal field, California. We estimated the potential for inducing seismicity by coupled thermal-hydrological-mechanical analysis of the geothermal steam production and cold water injection to calculate changes in stress (in time and space) and investigated if those changes could induce a rock mechanical failure and associated MEQs. An important aspect of the analysis is the concept of a rock mass that is critically stressed for shear failure. This means that shear stress in the region is near the rock-mass frictional strength, and therefore very small perturbations of the stress field can trigger an MEQ. Our analysis shows that the most important cause for injection-induced MEQs at The Geysers is cooling and associated thermal-elastic shrinkage of the rock around the injected fluid that changes the stress state in such a way that mechanical failure and seismicity can be induced. Specifically, the cooling shrinkage results in unloading and associated loss of shear strength in critically shear-stressed fractures, which are then reactivated. Thus, our analysis shows that cooling-induced shear slip along fractures is the dominant mechanism of injection-induced MEQs at The Geysers

A New Parameter to Assess Hydromechanical Effect in Single-hole Hydraulic Testing and Grouting

Grouting or filling of the open voids in fractured rock is done by introducing a fluid, a grout, through boreholes under pressure. The grout may be either a Newtonian fluid or a Bingham fluid. The penetration of the grout and the resulting pressure profile may give rise to hydromechanical effects, which depends on factors such as the fracture aperture, pressure at the borehole and the rheological properties of the grout. In this paper, we postulate that a new parameter, {angstrom}, which is the integral of the fluid pressure change in the fracture plane, is an appropriate measure to describe the change in fracture aperture volume due to a change in effective stress. In many cases, analytic expressions are available to calculate pressure profiles for relevant input data and the {angstrom} parameter. The approach is verified against a fully coupled hydromechanical simulator for the case of a Newtonian fluid. Results of the verification exercise show that the new approach is reasonable and that the {angstrom}-parameter is a good measure for the fracture volume change: i.e., the larger the {angstrom}-parameter, the larger the fracture volume change, in an almost linear fashion. To demonstrate the application of the approach, short duration hydraulic tests and constant pressure grouting are studied. Concluded is that using analytic expressions for penetration lengths and pressure profiles to calculate the {angstrom} parameter provides a possibility to describe a complex situation and compare, discuss and weigh the impact of hydromechanical couplings for different alternatives. Further, the analyses identify an effect of high-pressure grouting, where uncontrolled grouting of larger fractures and insufficient (or less-than-expected) sealing of finer fractures is a potential result

Hydromechanical modeling of pulse tests that measure both fluid pressure and fracture-normal displacement of the Coaraze Laboratory site, France

21International audienceIn situ fracture mechanical deformation and fluid flow interactions are investigated through a series of hydraulic pulse injection tests, using specialized borehole equipment that can simultaneously measure fluid pressure and fracture displacements. The tests were conducted in two horizontal boreholes spaced one meter apart vertically and intersecting a near-vertical highly permeable fault located within a shallow fractured carbonate rock. The field data were evaluated by conducting a series of coupled hydromechanical numerical analyses, using both distinct-element and finite-element modeling techniques and both two- and three-dimensional model representations that can incorporate various complexities in fracture network geometry. One unique feature of these pulse injection experiments is that the entire test cycle, both the initial pressure increase and subsequent pressure fall-off, is carefully monitored and used for the evaluation of the in situ hydromechanical behavior. Field test data are evaluated by plotting fracture normal displacement as a function of fluid pressure, measured at the same borehole. The resulting normal displacement-versus-pressure curves show a characteristic loop, in which the paths for loading (pressure increase) and unloading (pressure decrease) are different. By matching this characteristic loop behavior, the fracture normal stiffness and an equivalent stiffness (Young's modulus) of the surrounding rock mass can be back-calculated. Evaluation of the field tests by coupled numerical hydromechanical modeling shows that initial fracture hydraulic aperture and normal stiffness vary by a factor of 2 to 3 for the two monitoring points within the same fracture plane. Moreover, the analyses show that hydraulic aperture and the normal stiffness of the pulse-tested fracture, the stiffness of surrounding rock matrix, and the properties and geometry of the surrounding fracture network significantly affect coupled hydromechanical responses during the pulse injection test. More specifically, the pressure-increase path of the normal displacement-versus-pressure curve is highly dependent on the hydromechanical parameters of the tested fracture and the stiffness of the matrix near the injection point, whereas the pressure-decrease path is highly influenced by mechanical processes within a larger portion of the surrounding fractured rock

A fully coupled three-dimensional THM analysis of the FEBEX in situ test with the ROCMAS code: prediction of THM behavior in a bentonite barrier

Abstract: This paper presents a fully coupled thermal-hydrological-mechanical analysis of FEBEX-a large underground heater test conducted in a bentonite and fractured rock system. System responses predicted by the numerical analysis-including temperature, moisture content, and bentonite-swelling stress-were compared to field measurements at sensors located in the bentonite. An overall good agreement between predicted and measured system responses shows that coupled thermal-hydrological-mechanical processes in a bentonite barrier are well represented by the numerical model. The most challenging aspect of this particular analysis was modeling of the bentonite's mechanical behavior, which at FEBEX turned out to be affected by gaps between prefabricated bentonite blocks. At FEBEX, the swelling pressure did not develop until a few months into the experiment when moisture swelling of bentonite blocks had closed the gaps completely. Moreover, the wetting of the bentonite took place uniformly from the rock and was not impacted by the permeability difference between the Lamprophyres dykes and surrounding rock