3D Digital Outcrop Model-based quantification of fracture intensity: the Discontinuity Intensity Calculator and Estimator (DICE) open-source application

In outcrop-based fracture studies, the quantification of fracture intensity is often limited by the limitations of the manual sampling technique, characterized by punctual measurements (e.g. sampling spot, scanline, scanwindow) and moderate biases (e.g. fracture length truncation, technical and personal errors). The proximal remote sensing technologies, as terrestrial or Uncrafted Aerial Vehicle (UAV)-based LiDAR and photogrammetry, can help to overcome these limitations due to the possibility to obtain high-resolution and accurate quantitative data from the digital twin of the outcrop, the so-called Digital Outcrop Model (DOM). The DOMs can be very useful in outcrop-based fracture studies because their analysis allows to obtain several quantitative information with manual and/or automatic methods and with continuity in each position of the outcrop, increasing the accuracy of the fracture intensity estimations. However, due to the novelty of DOM technology and the lack of well-defined DOM-based fracture sampling procedures, these huge fracture datasets are often difficult to study and interpret, and therefore, the benefits of the DOM cannot be fully exploited. 

For this reason we present a complete workflow based on the DICE (Discontinuity Intensity Calculator and Estimator) open-source MATLAB© application that allows to quantitatively characterize the fractures of rocky outcrops from the 3D Digital Outcrop Models (DOMs). The proposed workflow consists in the following steps: (1) fracture mapping onto the 3D DOMs; (2) calculation of the fractures dimension, position and orientation; (iii) determination by DICE algorithm of the discontinuity parameters (persistence/dimension, distribution, spacing and intensity) using different 3D sampling techniques (3D scanline, 3D circular scan window and spherical scan volume). The differences of these sampling techniques and the fracture intensity parameters that can be obtained (p10, p21, p32) are discussed, showing the advantages and limitations of each DICE method.