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 population of 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.

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.

A B S T R A C TModeling naturally fractured reservoirs requires a detailed understanding of the three-dimensional (3D) fracture-network characteristics, whereas generally only one-dimensional (1D) data, often suffering from sampling artifacts, are available as inputs for modeling. Additional fracture properties can be derived from outcrop analogs with the scanline method, but it does not capture their full two-dimensional (2D) characteristics. We propose an improved workflow based on a 2D field-digitizing tool for mapping and analyzing fracture parameters as well as relations to bedding. From fracture data collected along 11 vertical surface outcrops in a quarry in southeast France, we quantify uncertainties in modeling fracture networks. The fracture-frequency distribution fits a Gaussian distribution that we use to evaluate the intrinsic fracture density variability within the quarry at different observation scales along well-analog scanlines. Excluding well length as a parameter, we find that 30 wells should be needed to fully (i.e., steady variance) capture the natural variability in fracture spacing. This illustrates the challenge in trying to predict fracture spacing in the subsurface from limited well data. Furthermore, for models with varying scanline orientations we find that Terzaghi-based spacing corrections fail when the required correction angle is more than 60°. We apply the 1D well Kevin Bisdom holds a M.Sc. degree in petroleum engineering and geosciences from the Delft University of Technology, Netherlands. He is currently a Ph.D. candidate in the section of applied geology at the Delft University of Technology working on geomechanical and fluid flow modeling of fracture networks in folded subsurface structures using outcrop analogs in central Tunisia. is a senior geologist and geophysicist at Total. He is currently head of the naturally fractured reservoir (NFR) study team and scientific adviser for NFR operational and research and development projects. Bertrand holds a Ph.D. in structural geology from the Pierre and Marie Curie University of Paris and has 25 years of experience in the oil and gas industry with Shell and Total. received his master's degree at the University of Pisa (Italy) and the Ph.D. at the ETH-Zurich with a thesis on the tectonosedimentary evolution of the South-Alpine rifted margin. He was then at the VU Amsterdam studying passive margins and foredeep basins. Since 2010, he has been a professor for applied geology at the Delft University of Technology working mainly on fractured reservoirs.

No abstract

General rightsCopyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. Users may download and print one copy of any publication from the public portal for the purpose of private study or research.  You may not further distribute the material or use it for any profit-making activity or commercial gain  You may freely distribute the URL identifying the publication in the public portal If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim.