Inferring fault structures and overburden depth in 3D from geophysical data using machine learning algorithms – A case study on the Fenelon gold deposit, Quebec, Canada

We apply a machine learning approach to automatically infer two key attributes – the location of fault or shear zone structures and the thickness of the overburden – in an 18 km² study area within and surrounding the Archean Fenelon gold deposit in Quebec, Canada. Our approach involves the inversion of carefully curated borehole lithological and structural observations truncated at 480 m below the surface, combined with magnetic and Light Detection and Ranging survey data. We take a computationally low-cost approach in which no underlying model for geological consistency is imposed. We investigated three contrasting approaches: (1) an inferred fault model, in which the borehole observations represent a direct evaluation of the presence of fault or shear zones; (2) an inferred overburden model, using borehole observations on the overburden-bedrock contact; (3) a model with three classes – overburden, faulted bedrock and unfaulted bedrock, which combines aspects of (1) and (2). In every case, we applied all 32 standard machine learning algorithms. We found that Bagged Trees, fine K-nearest neighbours and weighted K-nearest neighbour were the most successful, producing similar accuracy, sensitivity and specificity metrics. The Bagged Trees algorithm predicted fault locations with approximately 80% accuracy, 70% sensitivity and 73% specificity. Overburden thickness was predicted with 99% accuracy, 77% sensitivity and 93% specificity. Qualitatively, fault location predictions compared well to independently construct geological interpretations. Similar methods might be applicable in other areas with good borehole coverage, providing that criteria used in borehole logging are closely followed in devising classifications for the machine learning training set and might be usefully supplemented with a variety of geophysical survey data types.