Fractures are ubiquitous in geological formations and can often have an impact on subsurface applications such as geothermal energy, groundwater management or CO2 storage. Quantifying the relationship between the uncertainties inherent to fracture networks and the corresponding flow behaviour for these applications remains an open challenge. Simulation studies that are based on outcrop analogues of fracture networks have yielded many new insights about heat and mass transfer in fractured geological formations but are restricted to a limited number of fracture network realizations, simplified assumptions about fracture network properties or deterministic models, making it difficult to analyse a wide range of uncertainties. This study introduces a flexible workflow that generates ensembles of geologically plausible fracture networks that can be based on statistical data from outcrop analogues. The fracture networks are generated using a computationally efficient approach that combines mechanical and statistical methods. The ensembles are then seamlessly linked to multi-purpose flow and transport simulations where the fractures are represented explicitly in a porous and permeable rock matrix. This approach can enable new uncertainty quantification methods, supported by machine-learning-based emulators, to analyse how fracture network properties, such as fracture intensity, fracture aperture or fracture orientation, influence heat and mass transfer in fractured geological formations. The workflow is illustrated using two classic example applications pertinent to fracture network modelling – one based on outcrop data to assess thermal behaviour in geothermal systems, and one synthetic study to analyse the transition from matrix-dominated to fracture-dominated flow – and released as open-source code.

Characterising fractures in geothermal reservoirs is crucial for understanding heat and fluid flow, as fractures control reservoir permeability. Due to data scarcity, estimating fracture network properties remains uncertain. Dynamic data, such as well tests, provides indirect insights into subsurface properties and workflows have been developed to illustrate how uncertainty in fracture data affects flow behaviour. However, they use simplified, randomly generated fracture geometries limiting their applicability to real-world scenarios. This study presents a machine learning workflow for characterizing fractured reservoirs using transient data, focusing on geothermal reservoirs. A comprehensive dataset of 5000 geologically consistent Discrete Fracture Networks (DFNs) was generated using GeoDFN and directly linked to MRST for simulations. The workflow then applies a k-medoids clustering approach, using dynamic time warping (DTW) as a distance metric, to cluster pressure responses with similar transient behaviour. We identified 18 distinct pressure behaviour. Linking clusters to fracture properties reveals that fracture intensity, aperture, and length have the most significant impact on pressure behaviour, while fracture set type was found to be the least important factor. Future work will extend this workflow to temperature transient data and apply advanced machine learning techniques for both forward and inverse modelling of fractured geothermal reservoirs.