Bayesian dual inversion of experimental telescope acceptance and integrated flux for geophysical muon tomography

N. Lesparre (1) , D. Gibert (1) , J. Marteau (2)

1 IPGP - Institut de Physique du Globe de Paris
2 IPNL - Institut de Physique Nucléaire de Lyon

Density tomography of rock volumes with cosmic muons involves telescopes equipped with pixelized matrices of scintillator strips able to simultaneously measure the flux of muons in hundredths of directions. The resulting muon radiography images are a measure of the amount of matter integrated along each line of sight inside the geological target. This information constitutes the primary data at the root of muon density 3-D tomography. Before being used for either interpretation or tomography inversion, the radiographies must be corrected from artefacts due to imperfect detection capacity of the detection matrices. We present a correction method based on a Bayesian inversion to construct a probabilistic model of the distorted telescope acceptance from which undistorted radiographies may be obtained. The method also allows to simultaneously derive a stochastic model for the incident flux of muons. The resulting non-linear inverse problem is solved with the Metropolis-annealing algorithm, which allows to easily implement symmetry constraints to reduce the non-uniqueness. An inversion of real data acquired with one of our field muon telescopes is presented and discussed.

Domaines: Instrumentations et Détecteurs [physics.ins-det]; Planète et Univers [physics]

N. Lesparre, D. Gibert, J. Marteau. Bayesian dual inversion of experimental telescope acceptance and integrated flux for geophysical muon tomography. Geophysical Journal International, 2012, 188, pp.490-497. ⟨10.1111/j.1365-246X.2011.05268.x⟩. ⟨in2p3-00666938⟩