3D rock fabric analysis using micro-tomography:​ an introduction to the open-source TomoFab MATLAB code

Petri, Benoît; Almqvist, Bjarne S.G.; Pistone, Mattia

Affiliations

¹Institute of Earth Sciences, University of Lausanne, Géopolis, 1015, Lausanne, Switzerland

²Université de Strasbourg, CNRS, IPGS UMR 7516, F-67000, Strasbourg, France

³Department of Earth Sciences, Uppsala University, Villavägen 16, 752 36, Uppsala, Sweden

³Department of Geology, Franklin College of Arts and Sciences, University of Georgia, 210 Field Street, Athens, GA, 30602-2501, USA

Abstract and keywords

The study of rock fabric properties (orientation, planar, linear, anisotropy) is key to unravelling the geological processes that generated them. With advancements in data acquisition and treatment, X-ray micro-computed tomography (μXCT) represents a powerful method to analyse the shape preferred orientation (SPO) of rock-forming elements, including minerals, aggregates, and pores, in the three-dimensional space. After reconstruction and segmentation of μXCT images, we developed a novel protocol to construct and analyse the fabric tensor, a second-rank symmetric tensor constructed using the orientation and the length of the three characteristic axes of each grain (simplified to a best fit ellipsoid). The analysis of the fabric tensor permits calculation of mean principal directions and associated confidence ellipses, and quantifies the degree of anisotropy (P′) and the shape (T) of the fabric ellipsoid by eigenvalue and eigenvector analysis.

We implement this method in the TomoFab open-source MATLAB package. The code integrates a graphical user interface (GUI) that allows the visualisation of the full set of ellipsoid orientation, shape, and size. Density plots and contouring can be utilised to identify fabrics graphically, and a full set of fabric parameters can be calculated based on the analysis of the fabric tensor and/or the analysis of each principal direction orientation tensor.

We demonstrate the versatility of TomoFab with synthetic datasets and a field- and laboratory-based investigation of a sample presenting a magmatic foliation and lineation, collected in the Mafic Complex within the lower crustal section of the Ivrea-Verbano Zone (North Italy). In the light of these developments, we stress that μXCT represents a pertinent tool for rock fabric analysis to characterise the SPO of rock components. This approach can be performed parallel or complementary to other rock fabric quantification methods (e.g., AMS, EBSD) and applied to various rock types. TomoFab is freely available for download at https://github.com/benpetri/tomofab.

Keywords: Microfabric analysis; Tensor analysis; Micro-tomography; X-ray computed tomography; MATLAB GUI

Links

DOI: 10.1016/j.cageo.2020.104444

WoS: 000526797800010

History

Received: 2019-03-25

Revised: 2020-02-11

Accepted: 2020-02-11

Online: 2020-02-15

Cited Works (5)

Year	Entry
1984	Harrigan TP, Mann RW. Characterization of microstructural anisotropy in orthotropic materials using a second rank tensor. J Mater Sci. March 1984;19(3):761-767.
1974	Whitehouse WJ. The quantitative morphology of anisotropic trabecular bone. J Microsc (Oxford). July 1974;101:153-168.
1997	Odgaard A. Three-dimensional methods for quantification of cancellous bone architecture. Bone. 1997;20(4):315-328.
1998	Smit TH, Schneider E, Odgaard A. Star length distribution: a volume-based concept for the characterization of structural anisotropy. J Microsc (Oxford). September 1998;191(3):249-257.
1997	Odgaard A, Kabel J, van Rietbergen B, Dalstra M, Huiskes R. Fabric and elastic principal directions of cancellous bone are closely related. J Biomech. 1997;30(5):487-495.