Muography applied to nuclear waste storage sites ANNALS OF GEOPHYSICS, 60, 1, 2017, S0106; doi:10.4401/ag-7381 S0106 Muography applied to nuclear waste storage sites Raffaello D’Alessandro1,2,*, Fabio Ambrosino3,4, Lorenzo Bonechi2, Luigi Cimmino1,2, David G. Ireland5, Ralf Kaiser5, David F. Mahon5, Nicola Mori1,2, Pasquale Noli4, Giulio Saracino3,4, Craig Shearer6, Lorenzo Viliani1,2, Guangliang Yang5 1 Università di Firenze, Dipartimento di Fisica e Astronomia, Sesto Fiorentino, Florence, Italy 2 Istituto Nazionale di Fisica Nucleare (INFN), Sezione di Firenze, Sesto Fiorentino, Florence, Italy 3 Università di Napoli Federico II, Dipartimento di Fisica, Naples, Italy 4 Istituto Nazionale di Fisica Nucleare (INFN), Sezione di Napoli, Naples, Italy 5 SUPA, School of Physics and Astronomy, University of Glasgow, Glasgow, Scotland, UK 6 National Nuclear Laboratory, Central Laboratory, Sellafield, Seascale, Cumbria, England, UK ABSTRACT Legacy storage sites for nuclear waste can pose a serious environmental problem. In fact, since certain sites date from the middle of the last cen- tury when safety protocols had not been properly established and strict bookkeeping was not enforced, a situation has evolved where the content of storage silos is basically known only with a large uncertainty both on quantity and quality. At the same time maintenance work on old storage structures is becoming ever more urgent and yet this work requires exactly that information which is now lacking on the type of waste that was stored inside. Because of the difficulty in accessing the storage silos and the near impossibility of making visual inspections inside, techniques have to be developed which can determine the presence or absence of heavy elements (i.e. uranium) within the structures. Muography is a very prom- ising technique which could allow the survey of previously inaccessible structures. We have begun an evaluation performing feasibility studies using simulations based on real case scenarios. This paper will outline the storage site scenarios and then present some of the results obtained from the Monte Carlo simulations. 1. Introduction The monitoring of legacy storage sites for nuclear waste is growing the interest of the nuclear industry. The limited knowledge of the content of some storage facilities coming from an epoch when there were no strict rules for bookkeeping impairs the ability to per- form the needed maintenance operations. Investigation campaigns on these sites are strongly constrained by factors like the difficult access to high-radiation envi- ronments and the impossibility to make visual inspec- tions. Given the critical importance of this topic in the frame of environment conservation, new inspection methods which cope with these difficulties must be de- veloped. Muon radiography is a technique that is attracting much attention in many fields related to the investiga- tion of the internal density structure of a given target body. Its application to the detection and monitoring of high-Z materials (such as uranium) inside various kind of containers ranging from a nuclear reactor [Morris et al. 2014, Miyadera et al. 2013] to cargos [Blackwell and Kudryavtsev 2013, Riggi et al. 2013] has been studied recently. Monte Carlo simulations also show that muon radiography can be successful in de- tecting cm-sized uranium samples inside small storage canisters filled with concrete [Clarkson et al. 2015]. This result, however, does not provide any firm evidence that the technique can scale up to storage facilities with typical dimensions of several meters. When enlarging the size of the inspected struc- ture, two kinds of problems arise. The first one is con- nected to the geometric acceptance of the detection device and thus to the data acquisition time needed to collect a reasonable amount of muons. Typically, the presence of high-Z materials is investigated by meas- uring the scattering position inside the target and the scattering angle of an incoming muon: a large number of muons will be scattered at large angles when tra- versing dense regions inside the inspected structure. This technique requires at least two tracking devices, one upstream of the structure, to measure the direc- Article history Received July 5, 2016; accepted October 5, 2016. Subject classification: Muons, Radiography, Nuclear waste. tion of the incoming muon, and one downstream, to measure the direction after traversing the structure. The distance between the trackers when they are placed at opposite sides of a large storing facility will reduce the geometric acceptance of the detection de- vice, making the inspection of large volumes unfeasible (although the monitoring of a localized sample with known position would still be possible by means of a proper arrangement of the detectors). The second problem with large facilities is the increased back- ground due to multiple scattering inside the structure itself. A shield of concrete which is several meters thick can add a non-negligible contribution to the scattering angle given by the passage through the high-Z region, and increase the probability of a large-angle scattering in the concrete itself. These effects must be assessed by means of dedicated studies. In this paper we present the preliminary results of a Monte Carlo study of the performance of a muon ra- diography technique based on muon absorption when applied to large storage volumes. Muon absorption is widely used in investigating large geological structures [Tanaka et al. 2001, Ambrosi et al. 2011, Cârloganu et al. 2013] thanks to its large acceptance given by the usage of a single, downstream detector. In the next sec- tions we describe the measurement scenario (storage silo, detector, measurement principle) and the data analysis techniques and then we briefly discuss the results. 2. Measurement scenario In order to obtain results which are relevant to the real world, a real storage scenario has been considered. The storage facility is a cylindrical silo filled with con- crete, with radius 3.5 m and height 4 m: this kind of fa- cility is present e.g. in the storage site managed by Sellafield Ltd. (UK). The hypothesized uranium con- tent of the real silos is an unknown number of cm-sized uranium debris in unknown locations. A Monte Carlo simulation based on Geant4 [Agostinelli et al. 2003] for the above described scenario has been set up. The simulated silo contains various uranium cubes with sizes of 2, 5 and 10 cm at different locations, and a cubic air bubble with a size of 10 cm in order to assess the effects of eventual low-density re- gions inside the concrete. The detector is made of two parallel planes spaced by 50 cm; each plane is (2×2)m2 wide and consists of two orthogonal layers of scintil- lating bars made of poly-vinyl-toluene. The assumed spatial resolution of the detector is 0.3 cm. These figures have been chosen in order to resemble a realistic, scaled- up muon detector for geological surveys [Ambrosi et al. 2011]. The simulation setup is shown in Figure 1. A total number of 1.2×109 events inside the de- tector acceptance have been simulated, which corre- sponds to a data acquisition time of about 160 days. The energy and angular spectra of the muon flux used for random particle generation have been parameterized from ground measurements performed with a mag- netic spectrometer [Bonechi et al. 2005]. The simulated energy range is 0.7 - 130 GeV, while the polar angle ranges from 0 (i.e. vertical muon) to 80 degrees. 3. Data analysis The absorption technique seeks for missing muons due to enhanced absorption when traversing high-den- sity regions inside the target. This requires an estimate of the expected number of detected events coming from any direction when no density fluctuation is present, in order to compute the amount of missing events from a given direction which will be proportional to the column density along that direction. For large targets like geo- logical structures this method can only provide infor- mation about the column density of the target along a line of sight, since the detector is point-like when com- pared to the dimension of the target, thus resulting in a two-dimensional density map. However, in the storage silo scenario the dimension of the detector is not negli- gible, and this fact can be exploited to obtain a full three-dimensional estimate of the target density. As shown in Figure 2, missing muon tracks due to a local- ized high-density region inside the target are detected along different directions, allowing for a stereoscopic re- construction of the position of the absorption point. De- tected tracks can be back-extrapolated to the silo region, the intersections with YZ planes at a given X can be com- puted and from those a YZ histogram of intersection points can be built. This histogram (the “back projection map”) will show a deficit of events with respect to the D’ALESSANDRO ET AL. 2 Figure 1. Schematic representation of the simulated scenario. The direction of arrival of detected muons is from right to left. 3 homogeneous target case in the portion of the his- togram (the “signal region”) where the YZ plane inter- sects the high-density region inside the silo. Together with the X position of the back-projection plane, the YZ position of the deficit will give the full three-dimensional location of the high-density region. If the plane does not intersect the high density region then the signal region will be wider and the signal-to-noise will be less, analo- gously to what happens to optical images outside the focal plane of the optical system. Many back-projection layers at different X coordinates can thus be considered, and the X position of the high-density region will be that of the layer giving the smallest signal region with the highest signal-to-noise ratio. This method is described in details in Bonechi et al. [2015]. The above described method requires the knowl- edge of the expected event count for a silo containing no uranium debris. This information cannot be meas- ured directly since the content of the silos is unknown and thus there is no known uranium-free silo. However, since uranium debris are much smaller than the silo and their positions are scattered randomly in the silo vol- ume, it seems reasonable to assume that the position of a debris inside a silo will not match any of the posi- tions of the debris inside another silo, so that the sig- nal regions of the back projection maps at a given X for the two silos will never overlap. With this assumption, subtracting the back projection map of silo A from that of silo B will show signal regions for B as regions with a deficit of events (direct signal), and those for A as an excess of events (inverse signal). In order to avoid the necessity of acquiring data for two silos, the mirror symmetry about the Z axis of the measurement sce- nario can be exploited. Assuming that no uranium de- bris will lie exactly on the Z axis, one can mirror the back projection map about the Z axis and subtract it from the direct map. In this way, every direct signal re- gion at coordinate (Y, Z) will have a mirrored inverse signal region at (−Y, Z). To investigate the presence of a signal from ura- nium debris and study its variation with acquisition time, simulated data have been analyzed by placing back projection layers so that they intersect the debris inside the silo. In a realistic scenario this would not be possible, since the positions of the debris would be un- known. A full stack of planes with a suitable regular spacing along X would be required to perform a tomo- graphic scan of the silo, and suitable algorithms to find the plane containing the signal region with minimum extension and maximum signal-to-noise would have to be implemented in order to determine the position of the debris. This is far beyond the scope of this work, which only aims at assessing the eventual presence of a signal and its time dependence. So for this goal we re- lied on the knowledge of the positions given by the Monte Carlo truth. The simulated data have been processed with the mirror map technique in order to produce event count difference maps. Each bin of the difference maps is computed as the difference of two independent Poisson random variables (i.e. event counts), and thus follows a Skellam distribution whose RMS is given by the square root of the sum of the ex- pectation values of the two Poisson variables, which have been estimated with the event count themselves. Signal-to-noise maps have thus been built by dividing each bin of the difference maps by its Skellam RMS. An example of these maps is shown in Figure 3. The back projection layer is set at X = −100 cm, where MUOGRAPHY APPLIED TO NUCLEAR WASTE STORAGE SITES Figure 2. Measurement principle: missing tracks (dashed lines) due to absorption in a high-density region (black box) “hit” the detector with different angles, allowing for a stereoscopic view of the ab- sorption point. Figure 3. Signal-to-noise map for an acquisition time of about 160 days computed at X = −100 cm. The direct (red) and inverse (blue) signal regions are clearly visible at Z ~ 350 cm. a 5 cm uranium cube is present at coordinates (Y, Z) = (−15, 350) cm (see Figures 1 and 2 for the details of the reference frame). The signal is clearly visible with a sig- nal-to-noise value for the single bin around 6, thanks to the uranium sample being placed near to the detector and in the upper part of the silo. This condition in- crease the statistics since the muon flux increase with decreasing polar angle, and makes the effect of multi- ple scattering in the concrete less effective due to the reduced thickness of traversed material. Indeed, the same 5 cm sample placed near the center of the silo does not produce a detectable signal. On the contrary, a 10 cm sample in the same position produce a signal- to-noise of about 7. The detection capability is thus strongly influenced by the position of the uranium de- bris, other than their size. To put a lower limit on this last parameter, 2 cm samples have been placed in vari- ous positions, but none of them resulted to give a de- tectable signal with the available statistics. Also, the 10 cm air bubble placed near the center of the silo resulted to be not visible as well. In order to define the lower signal-to-noise thresh- old for detection, the probability of having a statistical fluctuation in a single bin giving a signal-to-noise greater than a given value has been studied. When tak- ing into account the look-elsewhere effect given by the ~2000 bins corresponding to the silo region in the sig- nal-to-noise map, it turns out that the probability of having one or more bins with a signal-to-noise ratio greater than 4 due to statistical fluctuations is about 9%, decreasing to 0.08% when considering a threshold of 5. A signal-to-noise value of 5 have thus been as- sumed as a safe detection threshold. To assess the ac- quisition time needed to reliably detect a uranium sample, the analysis has been repeated by considering subsets of the whole simulation data set. A plot show- ing the dependence of the signal-to-noise ratio on the acquisition time for the 5 cm and the 10 cm visible sam- ples is shown in Figure 4. An acquisition time of 1-2 months result to be sufficient to detect both the 5 cm and the 10 cm samples. 4. Discussion In this work the feasibility of searching legacy nu- clear waste storage silos for uranium presence leverag- ing the absorption of atmospheric muons has been investigated. This is a novel application of a well estab- lished and developed technique in the field of the mon- itoring of geological structures, which present many peculiar aspects when applied to the nuclear storage scenario. This preliminary analysis based on Monte Carlo simulations shows that in a realistic case there is a good discovery potential for well-positioned uranium debris with a typical size greater than 5 cm. The re- quired data acquisition time is of the order of a couple of months. The study of systematic effects given by e.g. non-uniform concrete density, the presence of other heavy debris like steel rods, surrounding buildings, dif- ferent parameterizations of the muon flux etc. are be- yond the scope of this preliminary, and will be addressed in a future work. References Agostinelli, S., al. (2003). Geant4 - a simulation toolkit, Nucl. Instrum. Meth. A, 506, 250-303. Ambrosi, G., et al. (2011). The MU-RAY project: Vol- cano radiography with cosmic-ray muons, Nucl. In- strum. Meth. A, 628, 120-123. Blackwell, T.B., and V.A. Kudryavtsev (2013). Identifi- cation of Nuclear Materials in Cargo Containers Using Cosmic Rays, Nuclear Science Symposium and Medical Imaging Conference (NSS/MIC), 1-8. Bonechi, L., M. Bongi, D. Fedele, M. Grandi, S. Riccia- rini and E. Vannuccini (2005). Development of the ADAMO detector: test with cosmic rays at different zenith angles, Proceedings of the 29th International Cosmic Ray Conference, 9, 283-286. Bonechi, L., R. D’Alessandro, N. Mori and L. Viliani (2015). A projective reconstruction method of un- derground or hidden structures using atmospheric muon absorption data, J. Instrum., 10, P02003. Cârloganu, C., et al. (2013). Towards a muon radiogra- phy of the Puy de Dôme, Geosci. Instrum. Meth., 2 (1), 55-60. Clarkson, A., D.J. Hamilton, M. Hoek, D.G. Ireland, J.R. Johnstone, R. Kaiser, T. Keri, S. Lumsden, D.F. Mahon, D’ALESSANDRO ET AL. 4 Figure 4. Signal-to-noise ratios for the 10 cm (open squares) and the 5 cm (full squares) samples as functions of the data acquisition time. Error bars represent the statistical errors. 5 B. McKinnon, M. Murray, S. Nutbeam-Tuffs, C. Shearer, G. Yang and C. Zimmerman (2015). Char- acterising encapsulated nuclear waste using cosmic- ray muon tomography, J. Instrum., 10, P03020. Miyadera, H., K.N. Borozdin, S.J. Greene, E.C. Milner, C.L. Morris, Z. Luki , K. Masuda and J.O. Perry (2013). Imaging Fukushima Daiichi reactors with muons, AIP Advances, 3, 052133. Morris, C.L., et al. (2014). Analysis of muon radiogra- phy of the Toshiba nuclear critical assembly reactor, Appl. Phys. Lett., 104, 024110. Riggi, S., V. Antonuccio-Delogu, M. Bandieramonte, U. Becciani, A. Costa, P. La Rocca, P. Massimino, C. Petta, C. Pistagna, F. Riggi, E. Sciacca and F. Vitello (2013). Muon tomography imaging algorithms for nuclear threat detection inside large volume con- tainers with the Muon Portal detector, Nucl. In- strum. Meth. A, 728, 59-68. Tanaka, H., K. Nagamine, N. Kawamura, S.N. Naka- mura, K. Ishida and K. Shimomura (2001). Devel- opment of the Cosmic-Ray Muon Detection System for Probing Internal-Structure of a Volcano, Hyperf. Int., 138 (1), 521-526. *Corresponding author: Raffaello D’Alessandro, Istituto Nazionale di Fisica Nucleare (INFN), Sezione di Firenze, Sesto Fiorentino, Florence, Italy; email: candi@fi.infn.it. © 2017 by the Istituto Nazionale di Geofisica e Vulcanologia. All rights reserved. MUOGRAPHY APPLIED TO NUCLEAR WASTE STORAGE SITES << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles false /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Warning /CompatibilityLevel 1.3 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.1000 /ColorConversionStrategy /LeaveColorUnchanged /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams true /MaxSubsetPct 100 /Optimize false /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments false /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile (None) /AlwaysEmbed [ true /AndaleMono /Apple-Chancery /Arial-Black /Arial-BoldItalicMT /Arial-BoldMT /Arial-ItalicMT /ArialMT /CapitalsRegular /Charcoal /Chicago /ComicSansMS /ComicSansMS-Bold /Courier /Courier-Bold /CourierNewPS-BoldItalicMT /CourierNewPS-BoldMT /CourierNewPS-ItalicMT /CourierNewPSMT /GadgetRegular /Geneva /Georgia /Georgia-Bold /Georgia-BoldItalic /Georgia-Italic /Helvetica /Helvetica-Bold /HelveticaInserat-Roman /HoeflerText-Black /HoeflerText-BlackItalic /HoeflerText-Italic /HoeflerText-Ornaments /HoeflerText-Regular /Impact /Monaco /NewYork /Palatino-Bold /Palatino-BoldItalic /Palatino-Italic /Palatino-Roman /SandRegular /Skia-Regular /Symbol /TechnoRegular /TextileRegular /Times-Bold /Times-BoldItalic /Times-Italic /Times-Roman /TimesNewRomanPS-BoldItalicMT /TimesNewRomanPS-BoldMT /TimesNewRomanPS-ItalicMT /TimesNewRomanPSMT /Trebuchet-BoldItalic /TrebuchetMS /TrebuchetMS-Bold /TrebuchetMS-Italic /Verdana /Verdana-Bold /Verdana-BoldItalic /Verdana-Italic /Webdings ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 150 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.10000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 150 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.10000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.08250 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile (None) /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName (http://www.color.org) /PDFXTrapped /Unknown /CreateJDFFile false /SyntheticBoldness 1.000000 /Description << /ENU (Use these settings to create PDF documents with higher image resolution for high quality pre-press printing. The PDF documents can be opened with Acrobat and Reader 5.0 and later. These settings require font embedding.) /JPN /FRA /DEU /PTB /DAN /NLD /ESP /SUO /NOR /SVE /KOR /CHS /CHT /ITA >> >> setdistillerparams << /HWResolution [2400 2400] /PageSize [595.000 842.000] >> setpagedevice