Abstract
This study presents a comprehensive methodology for the detection and characterization of fractures in geological samples using X-ray computed tomography (CT). By combining convolution-based image processing techniques with advanced neural network-based segmentation, the proposed approach achieves high precision in identifying complex fracture networks. The method was applied to a marly limestone sample from the Maiolica Formation, part of the Umbria-Marche stratigraphic succession (Northern Apennines, Italy), a geological context where fractures often vary in size and contrast and are frequently filled with minerals such as calcite or clays, making their detection challenging. A critical part of the work involved addressing multiple sources of uncertainty that can impact fracture identification and measurement. These included the inherent spatial resolution limit of the CT system (voxel size of 70.69 μm), low contrast between fractures and the surrounding matrix, artifacts introduced by the tomographic reconstruction process (specifically the Radon transform), and noise from both the imaging system and environmental factors. To mitigate these challenges, we employed a series of preprocessing steps such as Gaussian and median filtering to enhance image quality and reduce noise, scanning from multiple angles to improve data redundancy, and intensity normalization to compensate for shading artifacts. The neural network segmentation demonstrated superior capability in distinguishing fractures filled with various materials from the host rock, overcoming the limitations observed in traditional convolution-based methods. Overall, this integrated workflow significantly improves the reliability and accuracy of fracture quantification in CT data, providing a robust and reproducible framework for the analysis of discontinuities in heterogeneous and complex geological materials.