FracAbut: A python']python toolbox for computing fracture stratigraphy using interface impedance

In Naturally Fractured Reservoirs (NFR) diffuse fractures arrangement results from mechanical stratigraphy and tectonic history during failure. Thus, modelling Discrete Fracture Network (DFN) requires to understand and to account for fracture relationships at bed-interface (abutment or crosscutting) in 3D through time (loading path). However, sampling fractures data meaningfully in subsurface has always been a challenge for geologist due to data scarcity. To better understand and forecast fracture networks in stratified rocks, we study outcrops with a focus on geometric relationships between stratigraphic interfaces and fractures. This paper presents an original python toolbox called FracAbut. It is composed of 1 main and 2 auxiliary codes that quantify the geometric relation between fractures and stratigraphic interfaces from 1D (wells, scan-line) and 2D (digital image, photographs data). We calculate the Interface Impedance (II) that accounts for fracture abutment (crossing or not), persistence (single- or multi-bed) and propagation polarity (upward or downward). For each stratigraphic interface FracAbut provides information on fractures (type, number) and interface sensitivity (coupling strength). First, we apply FracAbut on synthetic case studies, then, on naturally fractured and stratified carbonates in Berat, Albania. Using both 1D scan-line and 2D outcrop photograph, we show that i) a mechanical interface can have different coupling above and below based on propagation polarity, ii) FracAbut results can give useful insight on fracture transmissivity, iii) FracAbut is fast and efficient to quantify fracture patterns and classify mechanical interface impact; iv) they are no relation between bed thickness and fracture propagation.