The Dynamic Fracture Model Generator algorithm was developed by Michael Welch and Mikael Lüthje by combining Linear Elastic Fracture Mechanics (LEFM) with Subcritical Fracture Propagation Theory (SFP) to calculate the rate of fracture nucleation, propagation and interaction as a function of mechanical properties, deformation history and the geometry of the fractured layer.


It is based on a simple conceptual model of fracture growth:


  • We start with the simple case of an intact planar, homogeneous layer. This contains an initial population of small, circular microfractures (shown in red)
  • As strain is applied to the layer, a local stress concentration builds up around the fracture edges, and they start to grow. Eventually, they become big enough to span the entire layer.
  • At this point, they start to propagate laterally and can be modelled as elongated rectangular layer-bound fractures (shown in blue).
  • As more layer-bound fractures nucleate, they start to interact in two ways:


A. By stress shadow interaction between parallel fractures

B. By intersection with perpendicular or oblique fractures.

We can model the growth of the fracture network in two ways:


  • By modelling the nucleation and growth of individual fractures, giving an explicit Discrete Fracture Network (DFN) model. This is ideal for detailed geomechanical and flow simulations, as dilation and fluid flow can be simulated within each fracture but is computationally expensive.
  • By modelling the evolution of functions representing the fracture density and size distribution for the entire network. This allows easier simulation of very large networks containing potentially millions of fractures.

Advantages of DFN Generator

Current DFN tools


Dynamic Fracture Model Generator


Accuracy:


Since they are calibrated against limited and incomplete (1D) data, conventional tools often generate geologically unrealistic fracture models, and at best there will be a high level of uncertainty in the results.

Accuracy:


Since it simulates the physical processes of fracture nucleation and propagation, the model output will automatically honour the geology, geomechanical properties and structural evolution of the reservoir. The detail of the output will reflect the detail of the input data, but even using basic input data the outputs will be geologically consistent.


Ease of use:


Since there is very little constraint on the input parameters, a “trial and error” approach with a lot of tweaking is required to get reasonable results. It can therefore be very time consuming to generate models with the current tools.

Ease of use:


Since this is a deterministic model it is easy to set up and run. It can typically be run with data already available in most models (e.g. curvature of top horizon, default mechanical properties for the known lithology) – a first pass model can be generated in c.½ day. However it is easy to build more detailed models, e.g. taking into account lateral variability in mechanical property due to facies variations, using existing Petrel functionality.


Uncertainty and risk analysis:


Since the input parameters are very poorly constrained and often do not represent geological attributes at all, it is difficult to systematically generate multiple realistic realisations or to rank them in likelihood. It is of course easy to generate many realisations by changing input parameters but it is difficult to analyse these systematically to quantify the range of possible outcomes.

Uncertainty and risk analysis:


Since it simulates fracture growth, it automatically generates multiple realisations at different stages of fracture development, which will accurately reflect the distribution of fractures we would expect for a low to high intensity of deformation. Where there is uncertainty in e.g. mechanical properties, we can use the Petrel workflow functionality to generate multiple realisations very quickly.


Usefulness:


Since the output from the models is limited, and often merely reflects the input data put into the models, they are only able to provide limited input into quantifying parameters of interest (e.g. connected fracture volume, permeability anisotropy etc.)

Usefulness:


The model generates (or can easily be modified to generate) many parameters and properties (e.g. fracture size distributions, fracture connectivity) that are required for complex flow modelling calculations. These properties come out of the simulation, and cannot be calculated directly from the model inputs.


For more information, see the following references:


- Welch, M.J., Lüthje, M., Oldfield, S.J. 2000. Modelling the Evolution of Natural Fracture Networks. Springer Nature Switzerland AG., https://doi.org/10.1007/978-3-030-52414-2


- Welch, M.J., Lüthje, M., Glad, A.C. 2019. Influence of fracture nucleation and propagation rates on fracture geometry: Insights from geomechanical modelling. Petroleum Geoscience, 25, 470-489, https://doi.org/10.1144/petgeo2018-161