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.

@en

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.

@en

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.

@en

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.

@en

Many widely used numerical models of rock fracture based on mesoscale laboratory test characterisation of effective ‘intact’ strength parameters neglect microstructure effects. They therefore cannot explain grain boundary and pore effects on crack propagation and consequently are inadequate for models of rock destruction that exploit point and indentation stresses. Understanding deep drilling processes involving drill-bit buttons and/or water-jetting where rock loading is concentrated in domains with fewer mineral grains will therefore require models with microstructure. To investigate microscale failure mechanisms of granular rocks in diverse scenarios, we target a porous sandstone and introduce a novel workflow consisting of a computerized tomography (CT) based microstructure construction approach and a complementary mechanical numerical approach. The construction approach extracts the realistic rock microstructure and transforms the large voxel number CT-scan data into significantly fewer triangular elements. The finite-discrete element method (FDEM) with grain-based model (GBM) is adopted to solve the mechanics. The microscale failure mechanism of sandstone during the Brazilian test was thoroughly analysed using the numerical results together with the post failure CT-scan test data. The build-up of compressive and tensile stress chains, micro-crack nucleation, local relaxation, chain switching and final crack-path development exploiting pores was illustrated, revealing the micro-to-macro failure mechanism in time and space. Fracture paths in the specimens during Brazilian tensile test were dominated by the pores and the inter-grain boundaries. The tensile strength of the inter-grain joints was estimated to be at least 3.67 times the mesoscale specimen's intact tensile strength, while the pores account for 72.76% of the fracture path. The influence of the cementation distribution and microscale discontinuities was investigated with numerical cases.@en

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.

@en

We study the potential formation damage effects due to radial jet drilling on Bentheim sandstone, a homogeneous sandstone that is well established in terms of mechanical properties and widely used as reservoir analogue. We quantify how properties degrade with distance from a jetted hole on μm scale (porosity) and cm scale (compressive, tensile and acoustical properties). Moreover, we perform true-triaxial compression tests on samples with and without a jetted hole. It is concluded that, for this material, jetting has no direct influence on the surrounding area. No significant changes compared to intact material is found, nor can a significant change be detected with respect to distance to a jetted hole. Differences fall within the intra-block variability, and differences between blocks can be well explained by block-to-block variation. Differences in results of the true-triaxial compression tests can be attributed to a different sample geometry, which is corroborated by a numerical simulation. We conclude that the stress field around the jetted hole can therefore be well approximated by the Kirsch equations, modified for compression, and thereby accurately describe the stability of the lateral borehole.@en

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.

@en

Two 0.3 m × 0.3 m × 0.3 m cubic blocks of shale and coal were used for hydraulic fracturing experiments under true tri-axial stress conditions. The shale block used was highly homogeneous and without visible fractures, while the coal block contained a host of natural fractures. The mechanical and hydraulic properties of both rocks were characterized through multi-stage triaxial tests, Brazilian disk tests, and porosity and permeability measurements. A true tri-axial rock testing machine equipped with loading, pump and acoustic systems was used in the experiment. The acoustic system uses 48 transducers with active sources to repetitively generate and receive ultrasonic P/S wave pulses to reveal fracture initiation and growth. Before the experiment, initial seismic response of both blocks was recorded under hydrostatic stress conditions to characterize anisotropy and heterogeneity of the blocks as reference. Silicon oil was injected centrally into both blocks to create a hydrofracture under deviatoric stress conditions and the load, displacement, pump pressure and volume, and seismic response during the injection process were recorded. Results from two blocks are being compared in terms of hydrofracture geometry and seismic features.@en

Cretaceous Sandstone bodies in the subsurface of western Netherlands are already used for heating some of the greenhouses in that area. The reservoirs used are typically at depths between 1500 and 3000m, with temperatures generally <100 °C. For higher temperature applications deeper reservoirs are required. However, deeper reservoirs are subjected to higher effective pressures due to more overburden, which can lead to more compacted rocks, and thereby reduced permeability. We assess the effects of effective pressure on Triassic Buntsandstein, a formation targeted to act as a deep geothermal reservoir in the western Netherlands. Rock samples are acquired from laterally equivalent quarries and prepared for permeability measurements within a tri-axial apparatus. To determine anisotropy, cores are drilled both perpendicular and parallel to bedding. Experiments are conducted by maintaining hydrostatic confining pressure, stepwise increasing up to 700 bar (if still permeable enough for accurate measurements) and a pore pressure of 25 bar. At each step the permeability is assessed by imposing a number of constant flow rates and continuous measurement of the pore pressure difference between up and downstream reservoirs. Throughout the experiment the sample strain is measured in radial and axial directions, such that elastic constants can be determined and micromechanical mechanisms may be observed. In addition to measurements on in-tact rock samples, we also assess the effect of induced fracturing on permeability by similar measurements. First, rock samples are fractured within the tri-axial cell with normal jacketing to evaluate the stress conditions of failure. Secondly, the experiment is repeated using relatively strong jackets which remain sealing after sample failure, allowing for permeability measurements. Preliminary results show that an increase of confining pressure leads to a decrease of permeability by three orders of magnitude, from 1e-13 to 1e-16 m2. Anisotropy results in permeability parallel to bedding to be roughly one order of magnitude higher than perpendicular to it. Based on the collected data, the validity of the available exponential permeability-porosity-stress relationship is assessed and the model parameters with the best fitting characteristic is chosen for the selected formation. The established relationship is then used as an input for field scale numerical simulation of cold fluid circulation in Buntsandstein formation to predict the reservoir behavior over longer term of fluid circulation. The Finite Element Method is used to evaluate the reservoir behaviour during injection/production of the cold/hot fluid in a fully coupled poro-thermo-elastic environment. Weighted residual method is used for deriving the weak formulation of the mass-, momentum- and energy balance equations. Consequently the standard Galerkin approach is used for spatial discretization of the weak formulas. Temporal discretization is also carried out in a fully implicit manner to avoid the time-stepping limitation. The preliminary results of this study show a promising capacity of heat extraction from the Buntsandstein formation as a geothermal reservoir within western Netherlands. @en

Thermal fracturing is considered to be a potential mechanism to create additional fractures in geothermal fields. The injected cold water into the hot host rock suddenly cools down the host rock, causing a considerable shrinkage of the material and thus potentially increased local stresses that may potentially lead to the formation of cooling related fractures. This is likely to happen in the near wellbore environment or along existing faults or fractures, ie. areas where the hot rocks juxtaposed to cold fluids. In this research, we experiment with thermal fracturing by exposing heated granitic and basaltic samples with cold water to see the extend of the thermal microfracturing inside the samples at different temperatures. Before and after the heat treatment, the micro CT-scanner is used to get high-resolution 3D images of fracture planes and fracture network connectivity. Moreover, the porosity is measured before and after treatment by using the pycnometer to see the effect of the different temperatures. In addition, the changes in geomechanical behaviour are tested by using an unconfined compressive strength (UCS) apparatus on heat treated and non-heat treated samples. We compare the changes in Young Modulus, Poisson’s Ratio and ultimate strength of the various samples and record the influence of the thermal fractures on the stressdriven fracturing behaviour in the UCS test.@en

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.

@en