Many rocks contain planar heterogeneities, in the form of open fractures, veins and/or stylolites, but scarce data exist on how strength and fracture pattern formation is affected by the presence of a singular planar heterogeneity in an otherwise uniform matrix. The mechanics of stylolite-bearing and/or fractured limestone is of interest to several engineering applications, from quarries to subsurface gas or geothermal reservoirs. We have performed Brazilian Disc tests on pre-fractured Indiana limestone samples and Treuchtlinger Marmor discs which contain cohesive stylolites, investigating Brazilian test Strength and the resulting fracture pattern. All experiments were filmed, and where possible analyzed with particle image velocimetry. When viewed in 2D, the planar discontinuity was set at different rotation angles compared to the principal loading direction, where perpendicular to the loading direction is defined as 0⁰. The results show that all samples are weaker than their intact counterparts. For the pre-fractured Indiana limestone, there is 10–75% angle-dependent weakening. However, in the samples with a stylolite, strength is weakened by 35–75%, independent of direction. Several new cracks appeared when fracturing a stylolite-sample, where the orientation is heavily influenced by the stylolite orientation. The fracture pattern and associated stress drops are more complex for high angles. In these samples always more than one fracture formed, whereas in pre-fractured samples usually only one new fracture formed. This suggests a potential for more permeability increase when hydrofracturing a stylolite-rich interval. Comparison with Finite Element Models indicates that this difference in fracture pattern is caused by the strength contrast between the anastomosing stylolite zone and the matrix material, leading to stress concentrations effects. This causes (micro-) fracture nucleation to occur locally, promotes fracture coalescence and fracture growth at lower overall sample-load conditions compared to intact samples.

We present an overview of induced seismicity due to subsurface engineering in the Netherlands. Our overview includes events induced by gas extraction, underground gas storage, geothermal heat extraction, salt solution mining and post-mining water ingress. Compared to natural seismicity, induced events are usually small (magnitudes ≤ 4.0). However, due to the soft topsoils in combination with shallow hypocentres, in the Netherlands events exceeding magnitude 1.5–2.0 may be felt by the public. These events can potentially damage houses and infrastructure, and undermine public acceptance. Felt events were induced by gas production in the north of the Netherlands and by post-mining water ingress in the south-east. Notorious examples are the earthquakes induced by gas production from the large Groningen gas field with magnitudes up to 3.6. Here, extensive non-structural damage incurred and public support was revoked. As a consequence, production will be terminated in 2022 leaving approximately 800 billion cubic metres of gas unexploited. The magnitudes of the events observed at underground gas storage, geothermal heat production and salt solution mining projects have so far been very limited (magnitudes ≤ 1.7). However, in the future larger events cannot be excluded. Project- or industry-specific risk governance protocols, extensive gathering of subsurface data and adequate seismic monitoring are therefore essential to allow sustainable use of the Dutch subsurface now and over the decades to come.

Increasing reservoir connectivity to the wellbore and bypassing the damaged area is crucial in improving the productivity of the wells and enhancing the swept area. This has become feasible by a new technology called radial jet drilling (RJD), in which relatively long, small-diameter laterals can be drilled radially from the main wellbore. In this study, the authors attempt to gain a better understanding of the efficiency of a high-velocity jet drilling on chalk destruction, and also identify parameters controlling the jet drilling. For this purpose, two distinct outcrop chalks from Austin, Texas (US) and Northern Province, Welton (UK) are used in this study, which are analogs to the reservoir chalk in the North Sea. In conjunction with the jet drilling experiments, basic rock mechanics testing is carried out in order to correlate the rock strength and stiffness properties to the jet drilling performance. Jet drilling of boreholes is evaluated not only by varying the fluid and nozzle type and the fluid pressure at the nozzle, but also varying the jet drilling setup under unconfined and also confined stress fields resembling reservoir condition. Results of our study show a clear correlation of the rock strength (and stiffness) on the threshold pressure and specific energy required to break the rock. Tight chalk requires more than 30% higher pump pressure than used in soft chalk for breaking the chalk, having more than twice the strength properties. Soft chalk presents larger borehole size and better rate of penetration, both with water and acid-aided fluid, owing to its higher matrix permeability value, as well as lower mechanical properties that favor diffusion of the jet drilling fluid into the rock and faster erosion/breakage compared with tight chalk. Static nozzles create a larger surface area compared with rotating nozzles. The penetration rate of the nozzle is improved significantly under stress confinement. In addition, jet drilling in the direction of minimum principal stress (σ3) appears to be faster owing to localization of shear failure around the drilled hole induced by the differential stresses compared with the jet drilling in the direction of maximum principal stress (σ1) under isotropic stress or ambient conditions.

In every tight formation reservoir, natural fractures play an important role for mass and energy transport and stress distribution. Enhanced Geothermal Systems (EGS) make no exception, and stimulation aims at increasing the reservoir permeability to enhance fluid circulation and heat transport. EGS development relies upon the complex task of predicting accurate hydraulic fracture propagation pathway by taking into account reservoir heterogeneities and natural or preexisting fractures. In this contribution, we employ the variational phase-field method, which handles hydraulic fracture initiation, propagation, and interaction with natural fractures and is tested under varying conditions of rock mechanical properties and natural fractures distributions. We run bidimensional finite element simulations employing the open-source software OpenGeoSys and apply the model to simulate realistic stimulation scenarios, each one built from field data and considering complex natural fracture geometries in the order of a thousand of fractures. Key mechanical properties are derived from laboratory measurements on samples obtained in the field. Simulations results confirm the fundamental role played by natural fractures in stimulation's predictions, which is essential for developing successful EGS projects.

The mechanical dynamics of volcanic systems can be better understood with detailed knowledge on strength of a volcanic edifice and subsurface. Previous work highlighting this on Mt. Etna has suggested that its carbonate basement could be a significant zone of widespread planar weakness. Here, we report new deformation experiments to better quantify such effects. We measure and compare key deformation parameters using Etna basalt, which is representative of upper edifice lava flows, and Comiso limestone, which is representative of the carbonate basement, under upper crustal conditions. These data are then used to derive empirical constitutive equations describing changes in rocks strength with pressure, temperature, and strain rate. At a constant strain rate of 10^-5 s^-1 and an applied confining pressure of 50 MPa, the brittle-to-ductile transitions were observed at 975 °C (Etna basalt) and 350 °C (Comiso limestone). For the basaltic edifice of Mt. Etna, the strength is described with a Mohr-Coulomb failure criterion with μ ~ 0.704, C = 20 MPa. For the carbonate basement, strength is best described by a power law-type flow in two regimes: a low-T regime with stress exponent n ~ 5.4 and an activation energy Q ~ 170.6 kJ/mol and a high-T regime with n ~ 2.4 and Q ~ 293.4 kJ/mol. We show that extrapolation of these data to Etna's basement predicts a brittle-to-ductile transition that corresponds well with the generally observed trends of the seismogenic zone underneath Mt. Etna. This in turn may be useful for future numerical simulations of volcano-tectonic deformation of Mt. Etna, and other volcanoes with limestone basements.

Radial water jet drilling (RJD) is a method of enhancing heat recovery by accessing and connecting to high permeable zones within geothermal reservoirs. The wall rock geometry behind an advancing water jet borehole under in-situ conditions is largely unknown. Water jet drilling tests were performed on 300 mm cubical blocks of weak porous sandstone under true-triaxial boundary stress conditions at the Delft Technical University (DTU) rock mechanics laboratory. Some of these tests showed distinct breakout features depending on the applied stress field. Geometries of resulting boreholes are recovered using X-Ray CT scans, and are analysed using segmentation software (Avizo). The code Solidity, using a combined finite-discrete element method with a cohesive zone fracture model, simulates stress take-up and wall shearing giving breakouts comparable to the experiments. The results lead to the suggestion that criteria based on Kirsch solutions would be suitable to provide general guidance on in-situ stress and rock strength conditions free of breakouts. FEMDEM models appear well-suited to examine geometries and dimensions that can be sustained by given strengths under deeper in-situ conditions. Wall-rock failure and a process of jet-hole enlargement together with the potential benefits of greater heat recovery arising from larger holes is also briefly discussed.

In geodynamic numerical models of volcanic systems, the volcanic basement hosting the magmatic reservoir is often assumed to exhibit constant elastic parameters with a sharp transition from the host rocks to the magmatic reservoir. We assess this assumption by deriving an empirical relation between elastic parameters and temperature for Icelandic basalts by conducting a set of triaxial compression experiments between 200 °C and 1000 °C. Results show a significant decrease of Young's modulus from ∼38 GPa to less than 4.7 GPa at around 1000 °C. Based on these laboratory data, we develop a 2D axisymmetric finite-element model including temperature-dependent elastic properties of the volcanic basement. As a case study, we use the Snæfellsjökull volcanic system, Western Iceland to evaluate pressure differences in the volcanic edifice and basement due to glacial unloading of the volcano. First, we calculate the temperature field throughout the model and assign elastic properties accordingly. Then we assess unloading-driven pressure differences in the magma chamber at various depths in models with and without temperature-dependent elastic parameters. With constant elastic parameters and a sharp transition between basement and magma chamber we obtain results comparable to other studies. However, pressure changes due to surface unloading become smaller when using more realistic temperature-dependent elastic properties. We ascribe this subdued effect to a transition zone around the magma chamber, which is still solid rock but with relatively low Young's modulus due to high temperatures. We discuss our findings in the light of volcanic processes in proximity to the magma chamber, such as roof collapse, dyke injection, or deep hydrothermal circulation. Our results aim at quantifying the effects of glacial unloading on magma chamber dynamics and volcanic activity.