The positive impact that natural fractures can have on geothermal heat production from low-permeability reservoirs has become increasingly recognised and proven by subsurface case studies. In this study, we assess the potential impact of natural fractures on heat extraction from the tight Lower Buntsandstein Subgroup targeted by the recently drilled NLW-GT-01 well (West Netherlands Basin (WNB)). We integrate: (1) reservoir property characterisation using petrophysical analysis and geostatistical inversion, (2) image-log and core interpretation, (3) large-scale seismic fault extraction and characterisation, (4) Discrete Fracture Network (DFN) modelling and permeability upscaling, and (5) fluid-flow and temperature modelling. First, the results of the petrophysical analysis and geostatistical inversion indicate that the Volpriehausen has almost no intrinsic porosity or permeability in the rock volume surrounding the NLW-GT-01 well. The Detfurth and Hardegsen sandstones show better reservoir properties. Second, the image-log interpretation shows predominately NW-SE-orientated fractures, which are hydraulically conductive and show log-normal and negative-power-law behaviour for their length and aperture, respectively. Third, the faults extracted from the seismic data have four different orientations: NW-SE, N-S, NE-SW and E-W, with faults in proximity to the NLW-GT-01 having a similar strike to the observed fractures. Fourth, inspection of the reservoir-scale 2D DFNs, upscaled permeability models and fluid-flow/temperature simulations indicates that these potentially open natural fractures significantly enhance the effective permeability and heat production of the normally tight reservoir volume. However, our modelling results also show that when the natural fractures are closed, production values are negligible. Furthermore, because active well tests were not performed prior to the abandonment of the Triassic formations targeted by the NLW-GT-01, no conclusive data exist on whether the observed natural fractures are connected and hydraulically conductive under subsurface conditions. Therefore, based on the presented findings and remaining uncertainties, we propose that measures which can test the potential of fracture-enhanced permeability under subsurface conditions should become standard procedure in projects targeting deep and potentially fractured geothermal reservoirs.

Responsibly using the subsurface for geo-energy extraction or storage requires an accurate understanding of the static and dynamic behaviour of the targeted reservoirs. While mostly dependent on intrinsic rock properties (e.g. porosity and permeability), this behaviour is also believed to be significantly impacted by the presence of natural fractures. For example, depending on their opening, connectivity, intensity or cement infill, natural fractures can both channel or inhibit fluid flow, thereby causing significant local flow heterogeneities which can both increase or decrease geo-energy production values. Therefore, assessing the possible characteristics of natural fractures has become an integral component in complex reservoir studies.

In layered materials, the deformation style, orientation, confinement, and 3D connectivity of natural fractures is generally impacted by changes in sedimentary facies and alternations in mechanical properties. In this study we address this effect and perform a numerical sensitivity analysis. Mechanical properties, confining pressures, and interfacial frictions are varied for a three-layered model, to investigate and quantify the relation between contrasting material properties, the principal horizontal stresses and fracture behaviour (i.e. deformation style and orientation). Firstly, the results show that tensile stresses develop in the stiffer layers due to the contrasting elastic parameters. The magnitude of these stresses is dependent on the ratio between the elastic parameters of stiffer and softer layers (i.e. E_stiff/E_soft and ν_stiff/ν_soft). There are no horizontal tensile stresses, when applying a compressive horizontal confining pressure (approx. 1/5 of the applied vertical stress). Implementing an interfacial friction lower than 0.2 will result in decoupling of the layers, resulting in slip on the layer boundaries and no tensile stresses within the stiffer layers. Further, the acquired numerical results are in good agreement with previously conducted laboratory work. Finally, we discuss whether the presented results can be used for better relating contrasting mechanical properties to potential fracture deformation styles and orientations in layered outcrops or subsurface reservoirs.

Production from the Ekofisk Chalk Field in the North Sea is believed to be significantly influenced by the presence of a connected fault and fracture network. In the current study, we create a 3D seismic discontinuity cube which is representative of this network within the southern part of the Ekofisk Field. This is done using a multiscale workflow which integrates seismic fault and fracture detection with borehole image log interpretation from three horizontal well sections. The results show that faults and fractures are prevalent in the Ekofisk Formations. Within the study area, faults are mainly organised in three orientations: 1) WNW-ESE, 2) NNE-SSW and 3) NNW-SSE. Smaller E-W striking faults are also observed. The interpreted fractures show a similar pattern and are organized in four orientation groups: NW-SE, WNW-ESE, ENE-WSW and NE-SW. The analysis of seismic discontinuity data (i.e. faults and fractures detectable on seismic) indicates that most small-scale discontinuities occur in proximity to large faults, and that the Lower Ekofisk Formation is characterized by more widespread – and a higher intensity of small-scale seismic discontinuities. It is also demonstrated that along each studied well section, the extracted seismic discontinuities show a qualitative correlation with the image log interpretation. This correlation suggests that the 3D seismic discontinuity cube can serve as a proxy for the fault and fracture network in the southern part of the Ekofisk Chalk Field. Following from our key findings, we conclude that the presented workflow and results could provide a starting point for future studies assessing the impact of natural fractures in the Ekofisk – and other complex reservoirs.

In carbonate rocks, channelized fluid flow through fracture conduits can result in the development of large and connected karst networks. These cavity systems have been found in multiple hydrocarbon and geothermal reservoirs, and are often associated with high-permeability zones, but also pose significant challenges in drilling and reservoir management. Here, we expand on the observed interplay between fractures, fluid flow and large cave systems, using outcrop analysis, drone imagery and fluid-flow modelling. The studied carbonate rocks are heavily fractured and are part of the Salitre Formation (750–650 Ma), located in central Bahia (NE Brazil). Firstly, the fracture and cave network data show a similar geometry, and both systems depict three main orientations, namely; NNE–SSW, NW–SE and ESE–WNW. Moreover, the two datasets are dominated by the longer NNE–SSW features. These observed similarities suggest that the fractures and caves are related. The presented numerical results further acknowledge this observed correlation. These results show that open fractures act as the main fluid-flow conduits, with the aperture model defining the fracture-controlled flow contribution. Furthermore, the performed modelling highlights that geometrical features such as length, orientation and connectivity play an important role in the preferred flow orientations.

Open natural fractures allow fluids to flow, which is necessary for the production of low-permeable geothermal and petroleum reservoirs. These reservoirs often consist of lithological layers with significant variation in rock strength, which makes it difficult to predict fracture containment within the rock succession. In this paper, fractures are classified as contained, when they do not cross the layer interface of the adjacent layer, so that formation and growth are inhibited in the adjacent layer. This study concerns the impact of the differences in rock strength (i.e., mechanical contrast) between two adjacent brittle layers on the fracture containment in finely-layered rocks. Laboratory deformation tests and 2D finite element modelling were performed on three-layered samples with varying mechanical contrasts to examine the behaviour of fractures at different stress conditions. Fractures initiate as shear fractures in weak layers, and propagate with a steeper angle (tensile fracture) through the adjacent stronger layer. The mechanical contrast within a layered rock does not always act as a containment barrier, meaning that fractures do not cross the layer interface of the adjacent layer. This is due to differences in differential stress between the weakest layer and strongest layer. The mechanical contrast combined with the magnitude of the confining pressure has a significant influence on fracture containment. An increase in mechanical contrast and confining pressure prevents fractures from propagating into stronger layers. The results contribute to predict natural fracture containment in brittle sequences at shallow depth in the subsurface.

Natural fracture networks are commonly observed in tight carbonate and chalk reservoirs and are believed to have significant impact on the effective permeability and potential fluid flow behavior. For instance, production from chalk fields in the North Sea is believed to be aided by the presence of natural fracture systems. Apart from enhancing production, fractures can also result in channelized fluid flow and early water break-through. In this study, we propose a multiscale and data-driven workflow of automated fault extraction, image log interpretation and both inverse and forward modelling, to characterize and quantify potential inter-well fracture network geometries and densities. The workflow is exemplified on the Ekofisk chalk field situated in the Norwegian North Sea. The seismic and well data show that the Ekofisk fault and fracture system forms a connected system, which can be subdivided into four main structural orientations. Inverse modeling suggests that the orientations of the fracture/fault system can be explained by three separate normal faulting events. By implementing the structural data into a forward simulation, we characterize the potential inter-well fracture/fault network highlighting that fractures occur in clustered zones which follow the four main observed orientations.

Orthogonal fracture networks form an arrangement of open well-connected fractures which have perpendicular abutment angles and sometimes show topological relations by which fracture sets abut against each other, thus forming a nested network. Previous modelling studies have shown that orthogonal fractures may be caused by a local stress perturbation rather than a rotation in remote stresses. In this study, we expand on the implications of these local stress perturbations using a static finite element approach. The derived stress field is examined to assess the development of implemented microfractures. The results show that the continuous infill of fractures leads to a gradual decrease in the local tensile stresses and strain energies, and, therefore, results in the development of a saturated network, at which further fracture placement is inhibit. The geometry of this fully developed network is dependent on the remote effective stresses and partly on the material properties. Saturated networks range from: (1) a set of closely spaced parallel fractures; (2) a ladder-like geometry; and (3) an interconnected nested arrangement. Finally, we show that our modelling results at which we apply effective tension, are equivalent to having a uniformly distributed internal pore fluid pressure, when assuming static steady state conditions and no dynamic fluid behaviour.

The potential of shale reservoirs for gas extraction is largely determined by the permeability of the rock. Typical pore diameters in shales range from the μm down to the nm scale. The permeability of shale reservoirs is a function of the interconnectivity between the pore space and the natural fracture network present. We have measured the permeability of the Whitby Mudstone, the exposed counterpart of the Posidonia Shales buried in the Dutch subsurface and a possible target for unconventional gas, using different methods and established a correlation with the microstructures and pore networks present down to the nanometer scale. Whitby Mudstone is a clay rich rock with a low porosity. The permeability of the Whitby Mudstone is in the range of 10⁻¹⁸m²–10⁻²¹m². 2D microstructures of the Whitby Mudstone show no connected pore networks, but isolated pore bodies mainly situated in the clay matrix, whereas 3D data shows that connected pore networks are present in less compacted parts of the rock. A closely spaced interconnected fracture network is often required to speed up transport of fluids from the matrix into a producing well. For fluids within the matrix the nearest natural fracture is on average at a distance of approximately 10cm in the Whitby Mudstone. The combination of the permeability data and the porosity data with natural fracture spacing of the fractures present in outcrops along the Yorkshire coast (UK) resulted in new insights into possible fluid pathways from reservoir to well.

Orthogonal fracture networks (ladder-like networks) are arrangements that are commonly observed in outcrop studies. They form a particularly dense and well connected network which can play an important role in the effective permeability of tight hydrocarbon or geothermal reservoirs. One issue is the extent to which both the long systematic and smaller cross fractures can be simultaneously critically stressed under a given stress condition. Fractures in an orthogonal network form by opening mode-I displacements in which the main component is separation of the two fracture walls. This opening is driven by effective tensile stresses as the smallest principle stress acting perpendicular to the fracture wall, which accords with linear elastic fracture mechanics. What has been well recognized in previous field and modelling studies is how both the systematic fractures and perpendicular cross fractures require the minimum principle stress to act perpendicular to the fracture wall. Thus, these networks either require a rotation of the regional stress field or local perturbations in stress field. Using a mechanical finite element modelling software, a geological case of layer perpendicular systematic mode I opening fractures is generated. New in our study is that we not only address tensile stresses at the boundary, but also address models using pore fluid pressure. The local stress in between systematic fractures is then assessed in order to derive the probability and orientation of micro crack propagation using the theory of sub critical crack growth and Griffith’s theory. Under effective tensile conditions, the results indicate that in between critically spaced systematic fractures, local effective tensile stresses flip. Therefore the orientation of the least principle stress will rotate 90 , hence an orthogonal fracture is more likely to form. Our new findings for models with pore fluid pressures instead of boundary tension show that the magnitude of effective tension in between systematic fractures is reduced but does not remove the occurring stress flip. However, putting effective tension on the boundaries will give overestimates in the reduction of the local effective tensile stress perpendicular to the larger systematic fractures and therefore the magnitude of the stress flip. In conclusion, both model approaches indicate that orthogonal fractures can form while experiencing one regional stress regime. This also means that under these specific loading and locally perturbed stress conditions both sets of orthogonal fractures stay open and can provide a pathway for fluid circulation.