Horizontal Water Injection Wells: Injectivity and Containment of Injection-Induced Fractures

Horizontal water injection wells have the capacity to inject larger volumes of water and have a smaller surface footprint than vertical wells. We present a new quantitative analysis on horizontal well injectivity, injection scheme (matrix vs. fracturing), and fracture containment. To precisely predict injector performance and delineate safe operating conditions, we simulate particle plugging, thermo-poro-elastic stress changes, thermal convection and conduction and fracture growth/containment in reservoirs with multiple layers. Simulation results show that matrix injection in horizontal wells continues over a longer time than vertical injectors as the particle deposition occurs slowly on the larger surface area of horizontal wellbores. At the same time, heat loss occurs uniformly over a longer wellbore length to cause less thermal stress reduction and delay fracture initiation. As a result, the horizontal well length and the injection rates are critical factors that control fracture initiation and long-term injectivity of horizontal injectors. To predict fracture containment accurately, thermal conduction in the caprock and associated thermal stresses are found to be critical factors. We show that ignoring these factors underestimates fracture height growth. Based on our simulation analysis, we suggest strategies to maintain high injectivity and delay fracture initiation by controlling the injection rate, temperature, and water quality. We also provide several methods to design horizontal water injectors to improve fracture containment considering wellbore orientation relative to the local stress orientations. Well placement in the local maximum horizontal stress direction induces longitudinal fractures with better containment and less fracture turning than transverse fractures. When the well is drilled perpendicular to the maximum horizontal stress direction, the initiation of transverse fractures is delayed compared with the longitudinal case. Flow control devices are recommended to segment the flow rate and the wellbore. This helps to ensure uniform water placement and helps to keep the fractures contained.

Temporal Changes in Seismic Anisotropy as a New Eruption Forecasting Tool

<p>The orientation of crustal anisotropy changed by ~80 degrees in association with the 1995/96 eruption of Mt. Ruapehu volcano, New Zealand. This change occurred with a confidence level of more than 99.9%, and affects an area with a radius of at least 5 km around the summit. It provides the basis for a new monitoring technique and possibly for future mid-term eruption forecasting at volcanoes. Three deployments of seismometers were conducted on Mt. Ruapehu in 1994, 1998 and 2002. The fast anisotropic direction was measured by a semi-automatic algorithm, using the method of shear wave splitting. Prior to the eruption, a strong trend for the fast anisotropic direction was found to be around NW-SE, which is approximately perpendicular to the regional main stress direction. This deployment was followed by a moderate phreatomagmatic eruption in 1995/96, which ejected material with an overall volume of around 0.02-0.05 km3. Splitting results from a deployment after the eruption (1998) suggested that the fast anisotropic direction for deep earthquakes (>55 km) has changed by around 80 degrees, becoming parallel to the regional stress field. Shallow earthquakes (<35 km) also show this behaviour, but with more scatter of the fast directions. Another deployment (2002) covered the exact station locations of both the 1994 and the 1998 deployments and indicates further changes. Fast directions of deep events remain rotated by 80 degrees compared to the pre-eruption direction, whereas a realignment of the shallow events towards the pre-eruption direction is observed. The interpretation is that prior to the eruption, a pressurised magma dike system overprinted the regional stress field, generating a local stress field and therefore altering the fast anisotropic direction via preferred crack alignment. Numerical modelling suggests that the stress drop during the eruption was sufficient to change the local stress direction back to the regional trend, which was then observed in the 1998 experiment. A refilling and pressurising magma dike system is responsible for the newly observed realignment of the fast directions for the shallow events, but is not yet strong enough to rotate the deeper events with their longer delay times and lower frequencies. These effects provide a new method for volcano monitoring at Mt. Ruapehu and possibly at other volcanoes on Earth. They might, after further work, serve as a tool for eruption forecasting at Mt. Ruapehu or elsewhere. It is therefore proposed that changes in anisotropy around other volcanoes be investigated.</p>

Temporal Changes in Seismic Anisotropy as a New Eruption Forecasting Tool

<p>The orientation of crustal anisotropy changed by ~80 degrees in association with the 1995/96 eruption of Mt. Ruapehu volcano, New Zealand. This change occurred with a confidence level of more than 99.9%, and affects an area with a radius of at least 5 km around the summit. It provides the basis for a new monitoring technique and possibly for future mid-term eruption forecasting at volcanoes. Three deployments of seismometers were conducted on Mt. Ruapehu in 1994, 1998 and 2002. The fast anisotropic direction was measured by a semi-automatic algorithm, using the method of shear wave splitting. Prior to the eruption, a strong trend for the fast anisotropic direction was found to be around NW-SE, which is approximately perpendicular to the regional main stress direction. This deployment was followed by a moderate phreatomagmatic eruption in 1995/96, which ejected material with an overall volume of around 0.02-0.05 km3. Splitting results from a deployment after the eruption (1998) suggested that the fast anisotropic direction for deep earthquakes (>55 km) has changed by around 80 degrees, becoming parallel to the regional stress field. Shallow earthquakes (<35 km) also show this behaviour, but with more scatter of the fast directions. Another deployment (2002) covered the exact station locations of both the 1994 and the 1998 deployments and indicates further changes. Fast directions of deep events remain rotated by 80 degrees compared to the pre-eruption direction, whereas a realignment of the shallow events towards the pre-eruption direction is observed. The interpretation is that prior to the eruption, a pressurised magma dike system overprinted the regional stress field, generating a local stress field and therefore altering the fast anisotropic direction via preferred crack alignment. Numerical modelling suggests that the stress drop during the eruption was sufficient to change the local stress direction back to the regional trend, which was then observed in the 1998 experiment. A refilling and pressurising magma dike system is responsible for the newly observed realignment of the fast directions for the shallow events, but is not yet strong enough to rotate the deeper events with their longer delay times and lower frequencies. These effects provide a new method for volcano monitoring at Mt. Ruapehu and possibly at other volcanoes on Earth. They might, after further work, serve as a tool for eruption forecasting at Mt. Ruapehu or elsewhere. It is therefore proposed that changes in anisotropy around other volcanoes be investigated.</p>

Study on the Influence of In Situ Stress Distribution on the Stability of Roadway Surrounding Rock

In order to understand the instability characteristics of surrounding rock in the process of deep roadway excavation, a three-dimensional numerical model was established by FLAC3D to systematically analyze the influence of roadway surrounding rock stability under different in situ stress distribution forms, and the environmental coefficient of mining-induced stress η was defined, the larger the environmental coefficient of mining-induced stress is, the larger the surrounding rock stress environment is, and the range where the η coefficient is greater than 0.2 is called with the destruction-danger zone. When the initial vertical stress is maximum principal stress and minimum principal stress, by comparing the roadway along the middle ground stress direction and minimum or maximum in-situ stress direction, the variation characteristics of displacement, failure zone and failure hazard zone of roadway surrounding rock are obtained, which provides theoretical basis for the treatment of disaster accidents such as roadway surrounding rock instability and rock burst caused by deep high in-situ stress.

Experiment and Failure Analysis of Rock-Like Material Affected by Different Excavation Depths

In order to increase the understanding of the strength and failure mechanism of rock mass during tunnel excavation, a series of uniaxial compression tests were conducted on mortar specimen with cracks and holes by using a rock mechanics servo-controlled testing system. And by monitoring the experimental process, the initiation, propagation, and coalescence process of cracks were observed and characterized. According to the experimental results, the influences of the excavation depth on the mechanical parameters and fracture characteristics of mortar specimens with single hole and the ones with single-hole crack were analyzed in detail. In the specimens with single hole, the peak strength decreases with the increase of hole depth, but the peak strain and elastic modulus have no obvious linear correlation with the hole depth. And the position and angle of initial crack change differently with the increase of the hole depth. The position of initial crack moves from the side of the hole to the top of the hole. When the hole depth exceeds 50%, the crack initiation angle is no longer inclined to the axial stress direction, but parallel to the axial stress direction. In the specimens with single prefabricated crack, the wing-shaped secondary cracks are generated at the tip of the precrack, and the antiwing-shaped secondary cracks are generated at the tip when approaching the peak stress. However, in the specimens with single-hole crack, no antiwing-shaped crack appears. And when the hole depth reaches 80%, two wing-shaped cracks appear at the precrack tip. One of the new wing-shaped cracks appears in the direction of the extension line of the precrack.

Crustal Fault Zones (CFZ) as Geothermal Power Systems: A Preliminary 3D THM Model Constrained by a Multidisciplinary Approach

The Pontgibaud crustal fault zone (CFZ) in the French Massif Central provides an opportunity to evaluate the high-temperature geothermal potential of these naturally permeable zones. Previous 2D modeling of heat and mass transfer in a fault zone highlighted that a subvertical CFZ concentrates the highest temperature anomalies at shallow depths. By comparing the results of these large-scale 2D numerical models with field data, the depth of the 150°C isotherm was estimated to be at a depth of 2.5 km. However, these results did not consider 3D effects and interactions between fluids, deformation, and temperature. Here, field measurements are used to control the 3D geometry of the geological structures. New 2D (thin-section) and 3D (X-ray microtomography) observations point to a well-defined spatial propagation of fractures and voids, exhibiting the same fracture architecture at different scales (2.5 μm to 2 mm). Moreover, new measurements on porosity and permeability confirm that the highly fractured and altered samples are characterized by large permeability values, one of them reaching 10-12 m2. Based on a thermoporoelastic hypothesis, a preliminary 3D THM numerical model is presented. A first parametric study highlights the role of permeability, stress direction, and intensity on fluid flow. In particular, three different convective patterns have been identified (finger-like, blob-like, and double-like convective patterns). The results suggest that vertical deformation zones oriented at 30 and 70° with respect to the maximum horizontal stress direction would correspond to the potential target for high-temperature anomalies. Finally, a large-scale 3D numerical model of the Pontgibaud CFZ, based on THM coupling and the comparison with field data (temperature, heat flux, and electrical resistivity), allows us to explore the spatial geometry of the 150°C isotherm. Although simplified hypotheses have been used, 3D field data have been reproduced.