Research Article  Open Access
The Feasibility Study for Multigeometries Identification of Uranium Components Using PCALSSVM Based on Correlation Measurements
Abstract
The geometry of uranium components is one of the key characteristics and strictly confidential. The geometry identification of metal uranium components was studied using ^{252}Cf sourcedriven correlation measurement method. For the 3 uranium samples with the same mass and enrichment, there are subtle differences in neutron signals. Even worse, the correlation functions were disturbed by scatter neutrons and include “accidental” coincidence, which is not conductive to the geometry identification. In this paper, we proposed an identification method combining principal component analysis and leastsquare support vector machine (PCALSSVM). The results based on PCALSSVM showed that the training precision was 100% and the test precision was 95.83% of the identification model. The total precision of the identification model was 98.41%, which indicated that the identification model was an effective way to identify the geometry properties with the correlation functions.
1. Introduction
The geometry of uranium components is one of the key characteristics and strictly confidential [1, 2]. As for uranium components, they are very difficult to be detected in passive measurements because of the low neutron yield and the weak energy of gamma ray. Uranium components can be shielded easily. Hence, the active method is considered and widely used.
Correlation measurements used for in the field of nuclear material identification have been developed in recent decades [3–11]. In active measurements, ^{252}Cf sourcesdriven correlation measurement method is a popular one [12, 13]. Based on this method, ORNL (Oak Ridge National Laboratory) had designed a nuclear material identification system (NMIS) aiming to verify the fissile material [14–16]. This method has been used successfully in the past to identify various uranium components, but nearly all of them focused on mass and enrichment [7, 17–19]. Indeed, geometric difference exists practically, but few studies have been done.
In this paper, the research work of the geometry identification to the metal uranium components has been conducted. The uranium components have the same mass, enrichment, and different geometry. Just because the uranium components only have geometric difference, the geometry identification is very difficult by the correlation functions simply. By analyzing the features and differences of the correlation functions, we proposed an identification method combining principal component analysis and leastsquare support vector machine (PCALSSVM) [20, 21] to build an identification model, and the uranium components with different geometries can be effectively distinguished.
2. Methodology
2.1. ^{252}Cf SourcesDriven Correlation Measurement Method
In ^{252}Cf sourcesdriven correlation measurement method, ^{252}Cf is employed as an interrogating neutron source.^{ 252}Cf source, uranium component, neutron detectors, and relevant data acquisition and display system are included in the measurement setup. Figure 1 shows the flowchart of the verification system. What if the uranium component is irradiated by ^{252}Cf source, ^{235}U goes through fission readily within 10^{−14} seconds in average and emits gamma rays, neutrons, and photons [19]. The mutual effects between spontaneous fission neutrons and gamma rays from ^{252}Cf source and the uranium component are chiefly in three possible processes: penetration, scattering, and induced fission. Then the neutrons and gamma rays are detected and recorded by the detector. This method takes advantage of innumerable neutron counts and makes inferences about the properties of uranium components observed by the correlation function between the ^{252}Cf source and detector.
2.2. The Correlation Function
As shown in Figure 1, ^{252}Cf source is placed in an ionization chamber called channel 1 (CH1). The chamber’s operation mode is typically pulse mode so that individual source fission event is observable. Neutrons from the ^{252}Cf source interacting with the uranium component are detected by a scintillator detector, which is placed straightly behind the uranium component and called channel 2 (CH2). The uranium component is placed between the ^{252}Cf source and the scintillator detector.
Thus, we can obtain signals , from CH1 and CH2, respectively, and derive the correlation function between the two signals. The correlation function , which statistically reflects the correlation neutrons characteristics between two random signals, is defined aswhere τ is the delayed time and is the observing time. The correlation functions can depict some properties of uranium components, such as enrichment and mass.
Figure 2 shows the correlation function between the ^{252}Cf source and a detector measured by ORNL [22]. Figure 2 has several meaningful features marking as letters from (a) to (g). (a) denotes the emitting time of ^{252}Cf spontaneous fission. (b) shows the directly transmitted gamma rays recorded by the detector. (c) means the scattered gamma rays and directly transmitted neutrons. (d) corresponds to the scattered neutrons, the neutrons, and gamma rays from induced fission that are detected lately (e). Some of these divisions are not distinct in time but overlap. For a long time, slow decay (f) will reduce to the background level (g).
2.3. PCALSSVM
Principal component analysis (PCA) is a powerful multivariate statistical technique which usually applied to data treatment of high dimensionality. It is based on the reduction of the variable numbers for a lower value of the poriginal ones, in order to represent the characteristics of this data set, however, without loss of information. Therefore, it can cover the combination among the variables as well as grouping of samples. Least squaressupport vector machine (LSSVM) is a method based on the Statistical Learning Theory (SLT) that employs a leastsquare linear system as a cost function, resembling a regularized network. PCALSSVM deals with combination of the feature extraction and classification. Geometry identification of the uranium components is separated into two phases: the feature extraction from the correlation functions and feature reduction by PCA and (2) the classification by LSSVM classifier.
3. Simulations, Results, and Discussion
3.1. Simulation Setup and Results
Considering high radioactivity of highenrichment uranium components, the simulation method is commonly used in the first place instead of making some experiments. In fact, the simulation research will help us further improve the performance of realdata experiments. Correlation measurements were carried out by application of MC simulation. Figure 3 shows the structure of the simulation setup. The distance between the ^{252}Cf source and the detector is fixed for all measurements.
The particle transport code, Geant 4, was used for all simulations. Figure 4 shows the correlation functions acquired by the MC simulation. In order to make the results more accurate, we remove the transmitted gamma rays by our proposed “rejecting window” method [23, 24]. When ^{252}Cf spontaneous fission, directly transmitted neutrons, the scattered neutrons, and the neutrons from induced fission all are detected by the detector, the peaks appeared in the curve.
For the neutrons detected by our detectors, we choose part of them whose energy covers 0~20 MeV as the objects to be processed. The average energy of ^{252}Cf source spontaneous fission neutron is 2.13 MeV. The average energy for the induced neutrons is about 2 MeV. Furthermore, neutron energy decreases when neutrons collide with uranium nuclei, which is beneficial to improving the probability of ^{235}U fission reaction. Thus, neutrons with energy less than 2 MeV account for the vast majority. Neutrons with energy ranging 2~20 MeV are very few. Neutrons with energy larger than 20 MeV are rare. That is why we choose 0~20 MeV.
Compared with Figures 2 and 4, there are several differences: due to eliminating the transmitted gamma rays, the feature (b) does not exist in Figure 4; (2) the experiment setups are different, such as the mass and enrichment of the uranium component. But beyond that, both results are consistent qualitatively in the total correlation functions. Thus it can be seen that the MC code is skilled in correctly simulating the correlation experiments between the ^{252}Cf source and detector.
3.2. Geometry Identification Based on Correlation Measurements
Three kinds of uranium components with the same mass (16.9736 kg) and same enrichment (93.15%) but different geometry are used. Table 1 lists the parameters of 3 samples.

Figure 5 shows the correlation functions for the 3 uranium components in different energy intervals. Here we focus on neutron correlation events and remove the gamma rays firstly; there is no gamma peak in Figure 5. The directly transmitted neutrons are the first particles detected by CH2 and then the scattered neutrons and the fission neutrons induced from the uranium components. The correlation function distribution presents a wide range because neutrons move at different velocities. The distribution curves given in Figure 5 show that, as time goes on, the correlation function decreases with the neutron chain reaction progress. Considering the features of the correlation function and compared with the distribution curves after the neutron peak, Figure 5 shows inconsistently in different energy intervals and cannot realize effectively the geometry identification.
(a)
(b)
(c)
(d)
To analyze the parameters of 3 uranium components, which have the same enrichment, mass and volume, so they have the same neutron net production. The 3 samples have different geometry, so they have different superficial area. Different superficial areas lead to different neutron leakages. The bigger superficial area has more possibility of stimulating the heavy nuclei and to induce easily more fission. Considering the above, the geometry identification of the 3 samples is very difficult, because there are subtle differences in neutron signals. Even worse, the correlation functions include “real” coincidence and “accidental” coincidence as shown in Figure 6. It is worth noting in Figure 6 that the “accidental” coincidence is chiefly from scatter neutrons which is not conductive to the geometry identification. Unfortunately, the detector just records the neutron counts but cannot pick out the fission neutrons or scatter neutrons. So, the geometry identification is very difficult by the correlation functions simply.
3.3. Geometry Identification Model Based on PCALSSVM
In order to identify the geometry of uranium components efficiently, we proposed a method combined principal component analysis and leastsquare support vector machine (PCALSSVM) to build an identification model based on the correlation functions of the 3 different geometry properties. In the above research, the correlation functions were divided into different groups with an energy interval of 0.1 MeV. The sample number of each geometry property was 21. So 21 × 3 = 63 groups of correlation functions were obtained. The 63 groups of correlation functions were shown in Figure 7.
In Figure 7, with the increase of the energy, the correlation functions of the sphere in a high energy level were mixed with the correlation functions of the cylinder in a low level. The similar situation also existed in the correlation functions of cylinder and cuboid, which is consistent with the above results. The geometry property was unable to be identified directly with the correlation functions.
A geometry identification model based on PCALSSVM was constructed. Firstly, PCA was applied to extract the feature components from the 63 groups of correlation functions. The feature components were the linear combination of the correlation functions and they were able to characterize the original data well. According to the extraction through PCA, the redundant information was deleted and the dimension was reduced dramatically compared with the original data. The number of the feature components was related to the including information percentage. Normally, an information percentage of no less than 85% contained in feature components was acceptable. In the experiment, the information percentage was selected as 95% and 3 feature components were obtained. The coefficients of the first feature component were shown in Figure 8.
The coefficients of the first feature component were the weights in characterizing the geometry properties with the correlation functions. When building the identification model of the geometry properties, the feature components were applied to replace the correlation functions as the input of the identification model. The labels of the 3 geometry properties were sets 1, 2, and 3 as output of the model corresponding to sphere, cylinder, and cuboid, respectively. Then LSSVM was applied to train the identification model of the geometry properties. In the model training process, the correlation functions were divided into training set and test set with a proportion of 13 : 8. Radial basis function was selected as the kernel function. The regularization parameter and kernel width of the LSSVM were determined by particle swarm optimization (PSO) method. At last, the identification model was trained and shown in Figure 9. In Figure 9, the distributions of each sample in the principal component feature space were shown. Coordinate axis showed the value of first feature component of each sample while the showed the second. The identification model divided the space into three parts. Each part represented the region where the feature components of the corresponding geometry property distributed.
The training precision was 100% and the test precision was 95.83% of the identification model. The total precision of the identification model was 98.41%, which indicated that the identification model was an effective way to identify the geometry properties with the correlation functions. However, the number of the sample data shown in this paper was not rich enough. In the future work, more data is expected to be achieved to test the reliability and precision of the identification model. Moreover, the identification model can be optimized for the possibility of identifying the geometry properties with different uranium masses and enrichments.
4. Conclusions
This paper describes a set of simulation experiments for the 3 uranium samples which have the the same mass and enrichment but different geometry. The experiments were performed using ^{252}Cf sourcesdriven correlation measurement method, based on the correlation functions. However, the correlation functions are disturbed by “accidental” coincidence chiefly from scatter neutrons and showed disagreement in different energy intervals. Therefore, we cannot realize the geometry identification by the correlation functions directly. We proposed further a method combined principal component analysis and leastsquare support vector machine (PCALSSVM). The results based on PCALSSVM showed that the training precision was 100% and the test precision was 95.83% of the identification model. The total precision of the identification model was 98.41%, which indicated that the identification model was an effective way to identify the geometry properties with the correlation functions. Obviously, this correlation function based identification method should be feasible not only for these three geometries but also for other geometries; it is one of our followup studies.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was partially supported by the National Natural Science Foundation of China (no. 11605017) and the Fundamental Research Funds for the Central Universities (no. 10611CDJXZ238826).
References
 X. W. Du, An elementary, of nuclear arms control technology, National Defense Industry Press, Beijing, China, 1996 (Chinese).
 C. A. Liu and J. Wu, An introduction of verification technology of nuclear arms control, vol. 3, National Defense Industry Press, Beijing, China, 2007.
 T. Uckan, M. Wyatt, J. Mihalczo, T. Valentine, J. Mullens, and T. Hannon, “^{252}Cfsourcecorrelated transmission measurements for uranyl fluoride deposit in a 24inOD process pipe,” Nuclear Instruments and Methods in Physics Research, vol. 422, no. 13, pp. 26–34, 1999. View at: Publisher Site  Google Scholar
 J. L. MuñozCobo, R. B. Perez, T. E. Valentine, Y. Rugama, and J. T. Mihalczo, “Stochastic transport theory of neutron and photon coupled fields: neutron and photon counting statistics in nuclear assemblies,” Annals of Nuclear Energy, vol. 27, no. 12, pp. 1087–1114, 2000. View at: Publisher Site  Google Scholar
 T. Valentine and J. Mihalczo, “MCNPDSP: A neutron and gamma ray Monte carlo calculation of sourcedriven noisemeasured parameters,” Annals of Nuclear Energy, vol. 23, no. 16, pp. 1271–1287, 1996. View at: Publisher Site  Google Scholar
 G. P. Marc, M. Monterial, P. Marleau et al., “Gamma/neutron timecorrelation for special nuclear material detection – Active stimulation of highly enriched uranium,” Annals of Nuclear Energy, vol. 72, pp. 358–366, 2014. View at: Publisher Site  Google Scholar
 S. A. Pozzi and J. Segovia, “^{252}Cf sourcecorrelated transmission measurements and genetic programming for nuclear safeguards,” Nuclear Instruments and Methods in Physics Research A, vol. 491, no. 12, pp. 326–341, 2002. View at: Publisher Site  Google Scholar
 M. Marseguerra, E. Padovani, and S. A. Pozzi, “Use of the MCNPPolimi code for timecorrelation safeguards measurements,” Progress in Nuclear Energy, vol. 43, no. 14, pp. 305–311, 2003. View at: Publisher Site  Google Scholar
 P. Feng, S. Y. Liu, and B. Wei, “Simulation and experimental study of a random neutron analyzing system with Cf252 neutron source,” Nuclear Science and Techniques, vol. 22, pp. 39–46, 2011. View at: Google Scholar
 J. T. Mihalczo and V. K. Paré, “Theory of correlation measurement in time and frequency domains with ^{252}Cf,” Annals of Nuclear Energy, vol. 2, no. 25, pp. 97–105, 1975. View at: Publisher Site  Google Scholar
 J. A. Williams, J. T. Mihalczo, C. W. Ricker, F. L. Glesius, and T. A. Kniss, “A highsensitivity, positionsensitive fission chamber for subcriticality measurements of spent fuel,” Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, vol. 299, no. 13, pp. 187–190, 1990. View at: Publisher Site  Google Scholar
 A. M. Stolle, “Application of stochastic diffusion theory to the interpretation of the ^{252C}f sourcedriven noise analysis method for subcriticality determination,” Annals of Nuclear Energy, vol. 20, no. 3, pp. 167–183, 1993. View at: Publisher Site  Google Scholar
 J. S. Li, Y. ZHang, and Y. Jin, “Development of ^{252}Cf fast fission chamber,” Nuclear Electronics & Detection Technlogy, vol. 21, no. 4, pp. 264–267, 2001. View at: Google Scholar
 S. A. Pozzi, R. B. Oberer, L. G. Chiang, J. K. Mattingly, and J. T. Mihalczo, “Higherorder statistics from NMIS to measure neutron and gamma ray cross talk in plastic scintillators,” Nuclear Instruments and Methods in Physics Research Section A, vol. 481, no. 13, pp. 739–748, 2002. View at: Publisher Site  Google Scholar
 B. R. Grogan and J. T. Mihalczo, “Simulating NMIS imaging and fissionmapping measurements,” INMM 52nd Annual Meeting, 2011. View at: Google Scholar
 J. Mihalczo, J. Mattingly, J. Neal, and J. Mullens, “NMIS plus gamma spectroscopy for attributes of HEU, PU and HE detection,” Nuclear Instruments and Methods in Physics Research Section B, vol. 213, pp. 378–384, 2004. View at: Publisher Site  Google Scholar
 S. Zhang and J. Wu, “Active method of delayed neutrons to determine uranium enrichment,” Nuclear Techniques, vol. 32, no. 12, pp. 922–926, 2009. View at: Google Scholar
 P. Huang, J. S. Li, and et al., “Influence of uranium components geometry on time correlation coincidence mass measurement,” Atomic Energy Science Technology, vol. 46, pp. 53–56, 2012. View at: Publisher Site  Google Scholar
 W. Xie, J. Li, J. Gong, J. Zhu, and P. Huang, “Experimental study on the measurement of uranium casting enrichment by timedependent coincidence method,” Chinese Physics C, vol. 37, no. 10, 2013. View at: Publisher Site  Google Scholar
 Y. Zhang, Y. Qin, Z.Y. Xing, L.M. Jia, and X.Q. Cheng, “Roller bearing safety region estimation and state identification based on LMDPCALSSVM,” Measurement, vol. 46, no. 3, pp. 1315–1324, 2013. View at: Publisher Site  Google Scholar
 D. Çalişir and E. Dogantekin, “A new intelligent hepatitis diagnosis system: PCALSSVM,” Expert Systems with Applications, vol. 38, no. 8, pp. 10705–10708, 2011. View at: Publisher Site  Google Scholar
 J. T. Mihalczo, J. A. Mullens, J. K. Mattingly, and T. E. Valentine, “Physical description of nuclear materials identification system (NMIS) signatures,” Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, vol. 450, no. 2, pp. 531–555, 2000. View at: Publisher Site  Google Scholar
 Y. Ren, B. Wei, J. Li, C. Ye, and P. Feng, “Implementation of ^{252}Cfsourcedriven power spectrum density measurement system,” High Power Laser and Particle Beams, vol. 24, no. 1, pp. 215–219, 2012 (Chinese). View at: Publisher Site  Google Scholar
 B. Wei, J. Jin, Y. Ren, and Y. Tang, “Method of γpeak removal in coincidence measurement system for verification of nuclear arms control,” High Power Laser and Particle Beams, vol. 23, no. 10, pp. 2781–2785, 2011 (Chinese). View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Mi Zhou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.