Fluid flow in naturally fractured reservoirs is often controlled by subseismic-scale fracture networks. Although the fracture network can be partly sampled in the direct vicinity of wells, the inter-well scale network is poorly constrained in fractured reservoir models. Outcrop analogues can provide data for populating domains of the reservoir model where no direct measurements are available. However, extracting relevant statistics from large outcrops representative of inter-well scale fracture networks remains challenging. Recent advances in outcrop imaging provide high-resolution datasets that can cover areas of several hundred by several hundred meters, i.e. the domain between adjacent wells, but even then, data from the high-resolution models is often upscaled to reservoir flow grids, resulting in loss of accuracy. We present a workflow that uses photorealistic georeferenced outcrop models to construct geomechanical and fluid flow models containing thousands of discrete fractures covering sufficiently large areas, that does not require upscaling to model permeability. This workflow seamlessly integrates geomechanical Finite Element models with flow models that take into account stress-sensitive fracture permeability and matrix flow to determine the full permeability tensor. The applicability of this workflow is illustrated using an outcropping carbonate pavement in the Potiguar basin in Brazil, from which 1082 fractures are digitised. The permeability tensor for a range of matrix permeabilities shows that conventional upscaling to effective grid properties leads to potential underestimation of the true permeability and the orientation of principal permeabilities. The presented workflow yields the full permeability tensor model of discrete fracture networks with stress-induced apertures, instead of relying on effective properties as most conventional flow models do.

Regional natural fracture networks often show variations on a scale below that captured by seismic reflection data. This variability is not considered in most reservoir models, but likely impacts uncertainties in permeability. We quantify this uncertainty using a database of 13,000 fractures in nine outcrops digitised in the carbonate Jandaíra Formation (Potiguar basin, Brazil). Distance between outcrops is on average 11 km, with a minimum of 300 m, which is comparable to the distance between wells in naturally fractured reservoirs. In between outcrops, significant variations exist in orientation, intensity, length and topology. Using discrete fracture-matrix flow models, we model the permeability of each deterministic pattern and find that small changes in geometry and topology result in permeability variations that are not captured by connectivity-based analyses such as percolation probabilities, particularly when the matrix is permeable. The permeability variations associated with subseismic-scale fracture variability are not captured in conventional stochastic models, but can be captured using deterministic outcrop models with flow through discrete fractures. The deterministic models provide a permeability range associated with subseismic fracture variability, that can be assigned to grid cells of fractured reservoir flow models, as an alternative to assuming constant permeability in the absence of subseismic-scale deformation.

The shallow-marine carbonate rocks of the Jandaíra Formation have been subject to significant permeability variations through time due to various events of fracturing and calcite cementation. As a consequence, the Jandaíra Formation accommodated fluid flow only during specific moments in time. We reconstructed these episodes of fluid flow based on isotope characterizations and microscope characteristics of calcite veins and host rock cements. The Jandaíra Formation, which belongs to the post-rift sequence of the Potiguar Basin in northeast Brazil, was deposited from the Turonian onward until a marine regression exposed it in the Campanian. Due to the subaerial exposure, meteoric waters flushed out marine connate waters, leading to an event of early diagenesis and full cementation of the Jandaíra Formation. Fluid flow through the resulting impermeable carbonate formation appears to be closely related to fracturing. Fracturing in the Late Cretaceous induced a drastic increase in permeability, giving rise to extensive fluid circulation. Host rock dissolution associated to the circulating fluids led to calcite vein cementation within the fracture network, causing it to regain an impermeable and sealing character. In the research area, fluid flow occurred during early burial of the Jandaíra Formation at estimated depths of 400–900 m. This study documents the first application of fluid inclusion isotope analysis on vein precipitates, which allowed full isotopic characterization of the paleo-fluids responsible for calcite vein cementation. The fluid inclusion isotope data indicate that upwelling of groundwater from the underlying Açu sandstones provided the fluids to the fracture network. In Miocene times, renewed tectonic compression of a lower intensity created a secondary fracture network in the Jandaíra Formation. The density of this fracture network, however, was too low to induce a new episode of fluid circulation. As a result, this tectonic event is associated with the development of barren extensional fractures.

Pervasive fracture networks are common in many reservoir-scale carbonate bodies even in the absence of large deformation and exert a major impact on their mechanical and flow behaviour. The Upper Cretaceous Jandaíra Formation is a few hundred meters thick succession of shallow water carbonates deposited during the early post-rift stage of the Potiguar rift (NE Brazil). The Jandaíra Formation in the present onshore domain experienced <1.5 km thermal subsidence and, following Tertiary exhumation, forms outcrops over an area of >1000 km². The carbonates have a gentle, <5⁰, dip to the NE and are affected by few regional, low displacement faults or folds. Despite their simple tectonic history, carbonates display ubiquitous open fractures, sub-vertical veins, and sub-vertical as well as sub-horizontal stylolites. Combining structural analysis, drone imaging, isotope studies and mathematical modelling, we reconstruct the fracturing history of the Jandaíra Formation during and following subsidence and analyse the impact fractures had on coeval fluid flow. We find that Jandaíra carbonates, fully cemented after early diagenesis, experienced negligible deformation during the first few hundreds of meters of subsidence but were pervasively fractured when they reached depths >400-500 m. Deformation was accommodated by a dense network of sub-vertical mode I and hybrid fractures associated with sub-vertical stylolites developed in a stress field characterised by a sub-horizontal σ₁ and sub-vertical σ₂. The development of a network of hybrid fractures, rarely reported in the literature, activated the circulation of waters charged in the mountainous region, flowing along the porous Açu sandstone underlying the Jandaíra carbonates and rising to the surface through the fractured carbonates. With persisting subsidence, carbonates reached depths of 800-900 m entering a depth interval characterised by a sub-vertical σ₁. At this stage, sub-horizontal stylolites developed liberating calcite which sealed the sub-vertical open fractures transforming them in veins and preventing further flow. During Tertiary exhumation, several of the pre-existing veins and stylolites opened and became longer, and new fractures were created typically with the same directions of the older features. The simplicity of our model suggests that most rocks in passive margin settings might have followed a similar evolution and thus display similar structures. Basin Research

Predicting equivalent permeability in fractured reservoirs requires an understanding of the fracture network geometry and apertures. There are different methods for defining aperture, based on outcrop observations (power law scaling), fundamental mechanics (sublinear length-aperture scaling), and experiments (Barton-Bandis conductive shearing). Each method predicts heterogeneous apertures, even along single fractures (i.e., intrafracture variations), but most fractured reservoir models imply constant apertures for single fractures. We compare the relative differences in aperture and permeability predicted by three aperture methods, where permeability is modeled in explicit fracture networks with coupled fracture-matrix flow. Aperture varies along single fractures, and geomechanical relations are used to identify which fractures are critically stressed. The aperture models are applied to real-world large-scale fracture networks. (Sub)linear length scaling predicts the largest average aperture and equivalent permeability. Barton-Bandis aperture is smaller, predicting on average a sixfold increase compared to matrix permeability. Application of critical stress criteria results in a decrease in the fraction of open fractures. For the applied stress conditions, Coulomb predicts that 50% of the network is critically stressed, compared to 80% for Barton-Bandis peak shear. The impact of the fracture network on equivalent permeability depends on the matrix hydraulic properties, as in a low-permeable matrix, intrafracture connectivity, i.e., the opening along a single fracture, controls equivalent permeability, whereas for a more permeable matrix, absolute apertures have a larger impact. Quantification of fracture flow regimes using only the ratio of fracture versus matrix permeability is insufficient, as these regimes also depend on aperture variations within fractures.

Aperture has a controlling impact on porosity and permeability and is a source of uncertainty in modeling of naturally fractured reservoirs. This uncertainty results from difficulties in accurately quantifying aperture in the subsurface and from a limited fundamental understanding of the mechanical and diagenetic processes that control aperture. In the absence of cement bridges and high pore pressure, fractures in the subsurface are generally considered to be closed. However, experimental work, outcrop analyses and subsurface data show that some fractures remain open, and that aperture varies even along a single fracture. However, most fracture flow models consider constant apertures for fractures. We create a stress-dependent heterogeneous aperture by combining Finite Element modeling of discrete fracture networks with an empirical aperture model. Using a modeling approach that considers fractures explicitly, we quantify equivalent permeability, i.e. combined matrix and stress-dependent fracture flow. Fracture networks extracted from a large outcropping pavement form the basis of these models. The results show that the angle between fracture strike and σ₁ has a controlling impact on aperture and permeability, where hydraulic opening is maximum for an angle of 15°. At this angle, the fracture experiences a minor amount of shear displacement that allows the fracture to remain open even when fluid pressure is lower than the local normal stress. Averaging the heterogeneous aperture to scale up permeability probably results in an underestimation of flow, indicating the need to incorporate full aperture distributions rather than simplified aperture models in reservoir-scale flow models.