Tomographic Imaging of Thermally Induced Fractures in Granite Using Bayesian Inversion

摘要

— The internal structure of rock samples studied in laboratory experiments can be described by a variety of physical parameters. Some of them, like the velocity of acoustic waves, enhanced velocity or quality factor can be reconstructed by means of ultrasonic tomography. This article presents the results of classical velocity tomography imaging, accompanied by the results of attenuation tomography and recently introduced enhanced velocity tomography obtained for a Lac Du Bonnet granite sample subjected to thermal stresses. To invert acoustic data recorded during six heating cycles, a Bayesian inversion scheme accompanied by a genetic algorithm optimization approach and the robust Cauchy norm have been used. To obtain the highest possible spatial resolution of images the inversion was performed in two steps. In the first step a crude parameterization of the sample was used. The result of this stage was next taken as an a priori model for a final inversion with refined parameterization. The choice of parameterization (cell sizes) and damping parameters at both stages was based on an analysis of the resolution operator. Both velocity and enhanced velocity tomography accurately imaged changes in the rock microstructure caused by thermal stresses. However, enhanced velocity tomography gave a much better spatial resolution than velocity tomography. On the other hand, attenuation tomography based on inversion of pulse rise times was able to image only a rough structure of the sample and it has difficulty with reasonable imaging of the crack formed in the sixth heating cycle.