Radiation Source Localization Using Surrogate Models Constructed from 3-D Monte Carlo Transport Physics Simulations

Abstract

Recent research has focused on the development of surrogate models for radiation source localization in a simulated urban domain. We employ the Monte Carlo N-Particle (MCNP) code to provide high-fidelity simulations of radiation transport within an urban domain. The model is constructed to employ a source location ($x,y,z$) as input and return the estimated count rate for a set of specified detector locations. Because MCNP simulations are computationally expensive, we develop efficient and accurate surrogate models of the detector responses. We construct surrogate models using Gaussian processes and neural networks that we train and verify using the MCNP simulations. The trained surrogate models provide an efficient framework for Bayesian inference and experimental design. We employ Delayed Rejection Adaptive Metropolis (DRAM), a Markov Chain Monte Carlo algorithm, to infer the location and intensity of an unknown source. The DRAM results yield a posterior probability distribution for the source’s location conditioned on the observed detector count rates. The posterior distribution exhibits regions of high and low probability within the simulated environment identifying potential source locations. In this manner, we can quantify the source location to within at least one of these regions of high probability in the considered cases. Employing these methods, we are able to reduce the space of potential source locations by at least 60%.

Keywords:

• Radiation detection
• inverse problem
• Bayesian inference
• MCNP
• surrogate modeling

Acknowledgments

This research was supported by the U.S. Department of Energy National Nuclear Security Administration under award number DE-NA0002576 through the Consortium for Nonproliferation Enabling Capabilities.