The Lilstock outcrop in the southern Bristol Channel provides exceptional exposures of several limestone beds displaying stratabound fracture networks, providing the opportunity to create a very large, complete, and ground-truthed fracture model. Here we present the result of automated fracture extraction of high-resolution photogrammetric images (0.9 cm/pixel) of the full outcrop, obtained using an unmanned aerial vehicle, to obtain a spatially extensive, full-resolution map of the complete fracture network with nearly 350,000 ground-truthed fractures. We developed graph-based functions to resolve some common issues that arise in automatic fracture tracing such as incomplete traces, incorrect topology, artificial fragmentation, and linking of fracture segments to generate geologically significant trace interpretations. The fracture networks corresponding to different regions within the outcrop are compared using several network metrics and the results indicate both inter- and intra-network (layer to layer) structural variabilities. The dataset is a valuable benchmark in the study of large-scale natural fracture networks and its extension to stochastic network generation in geomodelling. The dataset also highlights the intrinsic spatial variation in natural fracture networks that can occur even in weakly-deformed rocks over relatively short length scales of tens of metres.

Representing fractures explicitly using a discrete fracture network (DFN) approach is often necessary to model the complex physics that govern thermo-hydro-mechanical-chemical processes (THMC) in porous media. DFNs find applications in modelling geothermal heat recovery, hydrocarbon exploitation, and groundwater flow. It is advantageous to construct DFNs from the photogrammetry of fractured outcrop analogues as the DFNs would capture realistic, fracture network properties. Recent advances in drone photogrammetry have greatly simplified the process of acquiring outcrop images, and there is a remarkable increase in the volume of image data that can be routinely generated. However, manually digitizing fracture traces is time-consuming and inevitably subject to interpreter bias. Additionally, variations in interpretation style can result in different fracture network geometries, which, may then influence modelling results depending on the use case of the fracture study. In this paper, an automated fracture trace detection technique is introduced. The method consists of ridge detection using the complex shearlet transform coupled with post-processing algorithms that threshold, skeletonize, and vectorize fracture traces. The technique is applied to the task of automatic trace extraction at varying scales of rock discontinuities, ranging from 10° to 102m. We present automatic trace extraction results from three different fractured outcrop settings. The results indicate that the automated approach enables the extraction of fracture patterns at a volume beyond what is manually feasible. Comparative analysis of automatically extracted results with manual interpretations demonstrates that the method can eliminate the subjectivity that is typically associated with manual interpretation. The proposed method augments the process of characterizing rock fractures from outcrops.

In this paper, we present a two-dimensional numerical model for modelling of hydraulic fracturing in anisotropic media. The numerical model is based on extended finite element method. The saturated porous medium is modelled using Biot’s theory of poroelasticity. An enhanced local pressure model is used for modelling the pressure within the fracture, taking into account the external fluid injection and the leak-off. Directional dependence of all the rock properties, both elastic and flow related, is taken into account. A combination of the Tsai–Hill failure criterion and Camacho–Ortiz propagation criterion is proposed to determine the fracture propagation. We study the impact on fracture propagation (in both magnitude and direction) due to anisotropies induced by various parameters, namely ultimate tensile strength, Young’s modulus, permeability and overburden pressure. The influence of several combinations of all these anisotropies along with different grain orientations and initial fracture directions on the fracture propagation direction is studied. Different regimes are identified where the fracture propagation direction is controlled by the degree of material anisotropy instead of the stress anisotropy.