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.

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.

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.