Optimization of PGNAA device and algorithm for testing basicity index of sinter mixture

Long Zhao? Xu Xu? Jing-Bin Lu ? Ya-Lin Gong ? Qing-Feng Song ?Qing-Min Shang ? Wei Zhang ? De-You Yin ? Yan-Feng Li

Abstract Sinter is the main raw material in the blast furnace iron-making process, and basicity (CaO/SiO2) is an important quality index of sinter. Prompt gamma neutron activation analysis is a multi-elemental online detection technology that has been successfully applied in cement,coal, etc. Compared with cement as a raw material, sinter exhibits poor moderation ability and a large neutron absorption cross section. Therefore, cement detection devices are not suitable for sinter mixture detection.In this study, a prompt gamma neutron activation analysis equipment used for testing cement was re-optimized to render it suitable for measuring a sinter mixture. Using Monte Carlo simulation, the comprehensive detection efficiency of the detection device improved by 71.52%.Because of the gamma-ray self-shielding effect of the sinter mixture, the detection errors of CaO and SiO2 are significant. By applying the gamma-ray self-shielding correction algorithm, the detection accuracies improved,and their linear correlation coefficients R2 exceeded 0.99.Furthermore, by applying an improved analyzer to a sintering plant, the first-grade product rate of the factory increased by 4.64%.

Keywords Prompt gamma neutron activation analysis ·Sinter · Basicity · Monte Carlo simulation

1 Introduction

Prompt gamma neutron activation analysis is a rapid and noncontact multi-elemental analysis technique that has been widely used for elemental detection and analysis in various fields, such as cement, coal, and mineral resource industries [1–11]. The prompt gamma neutron activation analysis (PGNAA) technique is based on the detection of prompt gamma rays emitted through thermal neutron capture (nth, γ) or neutron inelastic scattering (n, n′γ). It can distinguish the elemental categories in the material from the characteristic gamma-ray energy spectrum as well as estimate the elemental content from the intensities of characteristic energy peaks in the spectrum. The short response time of PGNAA renders it appropriate for online process control tasks, such as online quality control in cement factories. Because neutrons have a strong penetrating ability, PGNAA technology can realize the fullsection detection of the conveyor belt, and the measurement results can reflect the overall situation. The neutron capture reaction emits gamma rays of high energy, and its detection results are not affected by environmental dust,steam, or other factors. This technology is suitable for industrial production lines with extremely poor environmental conditions.

Sinter is the main raw material in the blast furnace(BF)iron-making process, and basicity (CaO/SiO2) is an important quality index of sinter [12]. The fluctuation of basicity significantly affects the BF smooth operation.Every 10% reduction in basicity fluctuation results in an increase in production by 1.5% and coke reduction by 1%in BF. Currently, the control process of sintering production is lagging severely. The control lag is approximately 5 h,which results in a significant fluctuation in the basicity index of the sinter.Hence,real-time and accurate detection of the basicity is crucial to improve the sintering process and sinter quality. In this study, PGNAA analyzer model DF-5701, which was custom designed for the cement industry, was re-optimized to realize real-time online basicity detection. Because of the difference in chemical composition and physical parameters between the sinter mixture and cement raw material, the thermal neutron flux in the sinter mixture was only 62.58%of that in the cement raw material; therefore, the detection device must be optimized and upgraded. Using the Monte Carlo method for the optimization calculation, the comprehensive detection efficiency of the detection device increased by 71.52%. Gamma-ray self-absorption was observed in the sinter mixture, which contributed negatively to the detection results.In this study,the energy spectrum was restored by measuring the degree of gamma-ray self-absorption to improve the detection accuracy of CaO and SiO2. The linear correlation coefficient R2of the calibration samples of CaO and SiO2improved from 0.962 and 0.970 to 0.995 and 0.994,respectively.Finally,the upgraded analyzer was applied to the sintering plant. After applying the analyzer,the basicity stability rate of the plant increased by 3.98%,and the first-class product rate of the sinter increased by 4.96%.

2 Instrumentation and methods

In neutron-induced gamma techniques,the characteristic gamma-ray spectrum exhibits the following relationship with the neutron flux [13]:

Any unknown sample spectrum is the sum of the products of the elemental amount and the library spectrum of each element for every pulse-height channel [14].

The Monte Carlo simulation [15–17] model of PGNAA analyzer model DF-5701 is shown in Fig.1a. Two 20 μg252Cf neutron sources were placed in the source chamber.Two 5′′× 5′′N(xiāo)aI(Tl)scintillator detectors were used as the prompt gamma-ray detectors. The height of the analyzer belt gallery was 31 cm. The analyzer was suitable for a 1-m-wide belt. The simulated cement raw material load(belt material load per meter) was 100 kg/m, and the material thickness was 20.75 cm.

The main content of the sinter was iron oxide, which constituted more than 75% of the total material in weight.Iron oxide is a good gamma-ray shielding material; therefore,the gamma-ray energy spectrum of the sinter mixture exhibits self-absorption [18–21]. To quantify the self-absorption of gamma rays,a gamma-ray attenuation detection device was added after the PGNAA device. The system diagram is shown in Fig.1b. A137Cs radioactive source was installed under the belt, and a Φ 2′′× 2′′N(xiāo)aI(Tl)scintillator detector was installed above the material to detect the gamma-ray intensity.

The formula for gamma-ray self-absorption compensation is as follows [22]:

Fig.1 (Color online) Experimental model. a DF-5701 simulation model b system structure diagram

The linear absorption coefficient μ of the material is defined as follows:

The formula to calculate Nelementis as follows:

3 Experiments

3.1 Monte Carlo calculations

The main oxides and components of the cement raw material and sintered mixture are shown in Table 1. The thermal neutron capture cross sections of Ca,Si,Fe,and Al are shown in Table 2. According to Tables 1 and 2, the comprehensive thermal neutron capture cross section of the sinter mixture is 5.5 times that of the raw cement, indicating that the sinter mixture with the same load will consume more thermal neutrons than raw cement. Therefore, the detection device of the cement analyzer cannot directly detect the sinter mixture, and the device must be re-optimized.

The probability that the detector records the gamma ray emitted from position A is defined as the spatial gamma detection probability (SGDP) of position A. Furthermore,the spatial comprehensive detection probability (SCDP) at position A is defined as the product of the spatial thermal neutron distribution probability and spatial gamma detection probability at position A.Using the FMESH tally card[23] of MCNP, the spatial thermo neutron distribution probability (STNDP) and spatial fast neutron distribution probability (SFNDP) of the simulation area can be simulated easily. GTFADE is a Monte Carlo simulation software custom designed for simulating the SGDP. Using the GTFADE, the SGDP of the simulation area can be simulated easily.

The Monte Carlo simulation models are shown in Fig.2.In the model, the X-direction is vertical to the paper surface, representing the belt running direction and the length of the test sample; the Y-direction represents the height of the analyzer and the thickness of the test sample;and the Zdirection represents the width of the analyzer and the sample. A cuboid with a length, width, and height of 90,80,and 15 cm,respectively,was set as the simulation area.The simulation area was evenly divided into 90 × 80 × 15 = 108,000 small lattice cells, and the SGDP, STNDP, and SFNDP were simulated. Two252Cf spontaneous fission neutron sources were used in the model under the test material.The energy distribution of neutrons can be approximated as a Watt distribution [24]. The average neutron energy of the source was approximately3.5 MeV, and the distance between the two sources was 18 cm.

Table 1 Oxide composition and content of cement raw material and sinter mixture

Table 2 Thermal neutron capture cross-section data of four elements

Figure 2a shows the initial model (IM), which changes the test material from the cement raw material to the sinter mixture based on the DF-5701 model.The thickness of the sinter mixture was 15 cm.Because the density of the sinter mixture (1.9 tons per cubic meter) was greater than that of the cement raw material (1.3 tons per cubic meter), the height of the material gallery can be reduced. Figure 2b shows optimization model 1 (OM1), which reduced the height of the material gallery from 31 to 27.4 cm based on the IM. Figure 2c shows optimization model 2 (OM2),which reduced the height of the material gallery from 31 to 23.4 cm based on the IM. The neutron reflection plate in the model can enhance the STNDP of the detection area.Figure 2d shows optimization model 3 (OM3), which changed the material of the insulation board into highdensity polyethylene (HDPE) based on OM2. Figure 2e shows optimization model 4 (OM4), which reduced the thickness of the neutron reflector from 13 to 10 cm based on OM4.According to Kamal Hadad et al.[19],HDPE is a good neutron-moderating material. Figure 2f shows optimization model 5 (OM5), which adds a 2-cm-thick HDPE[25,26]above the gamma shield based on OM4 to enhance the neutron moderation ability of the model.

3.2 Experimental activities

In the experiment,the optimized device was used as the detection device, two 20 μg252Cf spontaneous fission neutron sources were used as the neutron excitation sources, and two sets of Φ 5′′× 5′′N(xiāo)aI(Tl) scintillator detectors were used as gamma-ray detectors.

Six groups of calibration samples were prepared in the experiment. The samples contained Fe2O3, SiO2, CaO,MgO, Cl, and other components. The components of the calibration samples are listed in Table 3. To verify the adaptability of the model, nine groups of verification samples were prepared in the experiment, and the sample composition is shown in Table 4.

The cement raw material and sinter mixture were used to perform different load experiments. The experimental results of the cement raw material are shown in Fig.3a.The analytical coefficient of iron AFeincreased linearly with the belt load, demonstrating that the gamma-ray selfabsorption of the cement as a raw material was weak. The experimental results for the sinter mixture are shown in Fig.3b. The analytical coefficient of iron element AFeexhibited a quadratic curve relationship with an increase in the belt load,and the curve exhibited a turning point,which indicates the strong gamma-ray self-absorption of the sinter mixture. Therefore, a compensation correction operation was required. As shown in Fig.3b, the sinter mixture load exhibited an approximately linear relationship in the range of 60–130 kg/m. The belt load directly affected the detection precision of each component in the experiment.The belt load must be in the linear region and should be as large as possible to increase the effective counting rate of the energy spectrum and improve the precision of the analysis results.In this experiment,110 kg/m was selected as the experimental belt load.

Figure 4a shows the PGNAA energy spectrum of the raw cement meal. As shown from the figure, the characteristic peaks of Ca and Si were significant and no superposition interference of other elemental characteristic peaks occurred. In such an ideal energy spectrum, the contents of CaO and SiO2can be calculated directly from the characteristic peak areas of Ca and Si.Figure 4b shows the PGNAA energy spectrum of the sinter mixture. As shown,the characteristic peaks of Fe were ubiquitous in the energy spectrum,and the characteristic peaks of Ca and Si were suppressed by the characteristic peaks of iron and its Compton plateau. Therefore, the component contents of CaO and SiO2with peak areas could not be calculated in such an energy spectrum.The spectral library least-squares method was used in the experiment to accurately deconvolute the energy spectrum. The spectrum library, which contains the characteristic energy spectra of calcium, silicon, iron,aluminum,magnesium,chlorine, sulfur, sodium,and background, was established prior to the experiment.Through the least-squares operation, we can obtain the contribution of each element in the total spectrum, which corresponds to the content of related elements.

Fig.2 (Color online) Monte Carlo simulation model. a IM b OM1 c OM2 d OM3 e OM4 f OM5

Six groups of calibration samples and nine groups of verification samples were placed in the PGNAA detection area to obtain the energy spectrum, and the energy spectrum collection time was 60 min. Subsequently, the samples were placed into the gamma attenuation detection area successively, and the corresponding the count rate with material on the belt N,the count rate with no material on the belt N0were recorded.

Table 3 Chemical composition of calibration samples (loss-free wt%)

Table 4 Chemical composition of verification samples (lossfree wt%)

Fig.3 Change curves of belt load. a Change curve of cement raw material load b change curve of sinter mixture load

Fig.4 PGNAA spectrum a PGNAA spectrum of cement raw material b PGNAA spectrum of sinter mixture

Fig.5 (Color online) SCDP of different test materials and STNDP with different locations. a SCDP of DF-5701 model when test material is cement raw material. b SCDP of DF-5701 model when test material is sinter mixture. c STNDP with different locations (IM vs. OM1)

Table 5 Different locations of STNDP data

4 Results and discussion

4.1 Monte Carlo calculation result

The simulation requirements included the neutron source following Watt distribution, a cross-section database of ENDF/B-VI [27], F4 cell flux tally, and FMESH card. A total of 1 × 109initial particles were used for better simulation statistics.

Different test materials of the same simulation model will exhibit different SCDPs because the test material is a neutron-moderating material.Figure 5a,b shows the SCDP simulation results for different test materials based on the DF-5701 model. The test material shown in Fig.5a is cement raw material,and the test material shown in Fig.5b is the sinter mixture.In the figure,the abscissa and ordinate represent the material length and width, respectively. As shown,the SCDP of cement is significantly higher than that of the sinter mixture. This is because the main element in the sinter mixture is iron, and its neutron-moderation ability is significantly less than those of light elements such as calcium and silicon;furthermore,the STNDP of cement is higher than that of the sinter mixture for the same neutron source. Table 2 shows that the thermal neutron absorption cross section of iron is much higher than that of calcium. Compared with the raw cement material, the sinter mixture consumed more thermal neutrons;therefore,the STNDP of the sinter mixture is much lower than that of the cement raw material.

As shown in Fig.5c and Table 5, the STNDP of OM1 increased by 36.54% at the top of the material but almost no increase was observed at the bottom near the radioactive source. As the height of the gallery decreased, the neutron reflection ability of the device improved;consequently,the thermal neutron flux on the upper part of the material increased.It can be inferred from Fig.5c that the effective approach to improve the SCDP of the system is to reduce the height of the material gallery; however, this is restricted.The passing capacity caused by the change in material shape must be ensured when the belt is operating.

Fig.6 (Color online) Simulation result curves. a Different detection efficiency curves of model DF-5701. b SCDP curves of different models.c Fast and thermo neutron distribution curves. d STNDP and SCDP comparison curves. e Simulation energy spectrum of IM and OM4 model

Figure 6a shows the distribution curves of the STNDP,SGDP, and SCDP with the material thickness in the DF-5701 model. The abscissa in the figure is the y-axis coordinate of the model,which represents the thickness change of the material. As shown, the STNDP decreased with increasing material thickness,the SGDP increased with the material thickness, and the SCDP increased with the material thickness. A thicker material enabled the absorption of more neutrons and hence a reduced STNDP. At locations closer to the detector, the spatial detection angle is larger; therefore, the SGDP increases with the thickness of the material.

Figure 6b and Table 6 shows the total comprehensive detection efficiency of each simulation model.As shown in the figure and table, OM4 achieved the highest total comprehensive detection efficiency. Figure 6c shows a comparison of the probability curves of fast neutron and thermal neutron distributions under two models: the DF-5701 model for cement material detection and OM4. As shown in the figure, the thermal neutron distribution probability of OM4 is lower than that of the cement model,whereas the fast neutron distribution probability is higher than that of the cement model, indicating that OM4 may exhibit problems in neutron moderation. To solve this problem,OM5 was designed.The comparison results of the OM5 and OM4 simulations are shown in Fig.6d. As shown, the comprehensive detection efficiency of OM5 islower than that of OM4,and that OM4 is the best detection model for this simulation calculation.

Table 6 Counts with different model data

Figure 6d shows that the probability of thermal neutron distribution at the bottom of the material in the OM5 model improved, but the probability of thermal neutron distribution at the top of the material decreased, and the comprehensive detection efficiency of the OM5 model decreased as well. This indicates that the sinter mixture possessed sufficient moderation ability for the fast neutrons emitted by the252Cf neutron source, and no additional moderation materials were required. However, the neutron moderation ability of cement, coal, bauxite, and other bulk materials was higher than that of the sinter mixture.Therefore,when using PGNAA technology to detect these materials, no additional moderation materials are required.

The MCNP code was used to simulate the characteristic gamma spectrum of the IM and OM4 models. A comparison diagram of the energy spectrum is shown in Fig.6e.It is evident that the characteristic peak counts of Fe,Ca,and Si in the characteristic energy spectrum of the OM4 model improved significantly.

4.2 Measurement experiment result

Table 7 shows the experimental data of the SiO2calibration samples. It includes the results of the chemical element test by the sampling method in the laboratory, the online results obtained using the library least-squares methods in PGNAA, and the results corrected by self-absorption. Table 8 shows the experimental data of the CaO before gamma self-absorption correction. b2 Calibration curve of calcium oxide after gamma self-absorption correction. c1 SiO2comparison curve of verification samples c2 CaO comparison curve of verification samples. d Basicity contrast curve calibration samples. Table 10 shows the experimental data of the SiO2verification samples. Table 9 shows the experimental data of the CaO verification samples. Figure 7a1, a2 shows the calibration curves of SiO2. Figure 7b1, b2 shows the calibration curves of CaO.Figure 7c1 shows the SiO2comparison curve of the verification samples. Figure 7c2 shows the CaO comparison curve of the verification samples. After the gamma-ray self-absorption correction, the calibration RMS error of SiO2increased from 0.333 to 0.148% (Table 7), and the linear correlation coefficient R2improved from 0.970 to 0.994 (Fig.7a1,a2).Meanwhile,the calibration RMS error of CaO increased from 0.398 to 0.143% (Table 8),and the linear correlation coefficient R2improved from 0.962 to 0.995 (Fig.7b1,b2).The CaO and SiO2RMS errors of the nine groups of verification samples were 0.152% and 0.205% (Tables 9 and 10), respectively.

Fig.7 (Color online) Curves of measurement experiment. a1 Calibration curve of silicon dioxide before gamma self-absorption correction. a2 Calibration curve of silicon dioxide after gamma self-absorption correction. b1 Calibration curve of calcium oxide

Table 10 Experimental data of SiO2 verification samples

Table 11 Error data of PGNAA analyzer in sintering plant

Table 12 Production data before and after using the analyzer

4.3 Application effect of analyzer

The improved analyzer was applied in a sintering plant and installed in the second mixing process point of the sintering batching workshop. The analyzer installation point was approximately 15 min behind the batching belt;however,it was approximately 4 h and 45 min ahead of the original control method.

In the actual sintering production process, the composition of each element changes slightly.By substituting the content of each element detected by the analyzer as a variable into the calculation formula, the dynamic detection μmEicurve can be realized, and the accuracy of the correction model improved. A comparison of the basicity data curves obtained using the analyzer and experimentally is shown in Fig.7d. It is clear that the basicity trend obtained using the PGNAA analyzer was consistent with that from performing a chemical test in the laboratory.The error data are shown in Table 11. The standard error of basicity was only 0.04,indicating the excellent accuracy of the PGNAA analyzer in detecting the basicity index of the sinter mixture. In August 2017, the sintering plant began using the analyzer to guide batching production. After the analyzer was used, the basicity stability rate (BSR) of the sintering plant improved by 3.98%, and the first-class product rate (FCPR)of the sinter improved by 4.64%.The factory production data are shown in Table 12.

5 Conclusion

In the experiment, the DF-5701 neutron activation analyzer was re-optimized, and the height of the material corridor was reduced from 31 to 23.4 cm. Thermal insulation materials were replaced with HDPE, and the thickness of the neutron reflector was reduced from 13 to 10 cm.The distance between the measured material and the detector was reduced from 43 to 32.4 cm, which effectively increased the comprehensive detection efficiency of the detection area.The comprehensive detection efficiency of the optimized model increased by 71.52%.

Author contributionsAll authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by Long Zhao,Qing-Feng Song and Xu Xu.The first draft of the manuscript was written by Long Zhao and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.