Open Access

Research on dose correction method of vehicle-borne environmental radiation measurement equipment


Cite

Introduction

The Fukushima accident in Japan in March 2011 caused severe radioactive nuclear leakage and spread to surrounding areas, which had a large-scale impact on the environment in a short period of time [1, 2]. An important measurement and the key to ensuring quality in the environmental radiation monitoring of nuclear accidents is the gamma radiation dose rate of the ambient surface [3]. The most commonly used technical monitoring is vehicle-borne environmental radiation measurement equipment due to its high mobility, rapid response, and the ability to obtain the distribution data of the environmental radiation field.

Since the beginning of the 1980s, developed countries such as the United States, Canada, and Japan have successfully developed vehicle-borne large-volume NaI gamma spectrometer detection systems for environmental radiation monitoring, nuclear accident emergency, and radioactive source detection [4, 5]. In recent years, China's environmental protection departments have actively conducted capacity building in environmental radiation monitoring. Domestic enterprises have also developed and produced a number of vehicle-borne radiation detection systems to assess the ambient dose rate and the dose rates of detected gamma-emitting nuclides. In 2007, the Ministry of Environmental Protection issued the “National standards for the construction of radiation environmental monitoring and supervision institutions”, which included the vehicle-borne radioactive detection system as the configuration object [6].

After the Fukushima nuclear accident in Japan, the stationary station nearby was destroyed by tsunami and other natural disasters. Thus, the vehicle-borne environmental radiation measurement equipment has a crucial role in environmental dose rate monitoring. The measurement principle of this equipment is similar to that of the stationary station. The difference is that the former has a strong mobile performance than the latter.

Vehicle-borne radiation detection system has been used to measure environmental dose rate [6, 7]. However, the use of this system for accurate monitoring of the environment and radioactive material radiation dose remains a difficult problem. The use of different measurements for the environmental dose value of a point can lead to inconsistent or even widely different results. At present, China has not released the relevant standards and technical specifications for the dose calibration of a vehicle-borne spectrometer. This study presents a new method to accurately calculate the surface environmental dose using vehicle-borne radiation measuring equipment. A near-ground reference radiation field is built on the basis of the typical radionuclides of the Fukushima nuclear accident. The Monte Carlo code FLUKA is used to simulate the ambient dose under different reference radiation fields. All simulated results are converted into air dose 1 m above the ground surface, as required by the International Commission on Radiological Protection (ICRP 1996) [8], to establish the dose correction factor (CF) for the vehicle-borne radiation environment measurement equipment.

Simulation methods and model construction

The FLUKA code [9] is used to simulate the environmental dose of vehicle-borne environmental radiation measurement equipment in the reference radiation fields. The incident gamma photons interact with the scintillator material in large-volume detectors, and part of the energy of gamma photons is deposited in the detector. The energy of gamma rays deposited in NaI (Tl) crystals is recorded using the USRTRACK card and converted into dose rate. The activity of each nuclide is 1 Ci at the initial time of all simulation data, and the simulation platform is Dell Server Precision T7920, CPU: Intel(R) Core (TM) i7-9750H @ 2.6 GHz, RAM 64.0 GB.

Model construction of vehicle-borne radiation environment measurement equipment

Vehicle-borne environmental radiation measurement technology measures the type and content of radionuclides on the ground and the air-absorbed dose rates along the road using the gamma spectrometer mounted on a car [4]. The Monte Carlo code FLUKA is used to simulate the ambient dose rate of vehicle-borne environmental radiation measurement equipment. A schematic of the FLUKA model used for simulations is illustrated in Figs. 1 and 2. The 4 L NaI (Tl) crystal is covered by a 0.5-mm-thick MgO reflector, which is surrounded by aluminium with a thickness of 2 mm, and by a layer of 2-mm-thick SiO2 on the photomultiplier coupling side. Given that NaI (Tl) crystals are susceptible to deliquescence, an aluminium box is used to seal them. The dry MgO powder between the crystal and the aluminium box acts as a reflective layer to reflect the light emitted in the scintillator in all directions to the crystal, and SiO2 can effectively transmit the light to the photomultiplier tube. The main component of the on-board radiation measurement equipment is a 4 L NaI (Tl) detector mounted on the roof (in horizontal placement) of a Cheetah Land Cruiser LBA6482LQ4. The point source is assumed to emit 1.332 MeV gamma rays, and the distance between the isotropic gamma-ray point source and the 10 cm × 10 cm surface of the NaI (Tl) crystal is 100 cm. The intrinsic detection efficiency of a NaI detector is defined as the ratio of the number of pulses deposited in the NaI (Tl) crystal to the incident radiation source on the surface of the detector. The simulation shows that the intrinsic detection efficiency of NaI (Tl) for 1.332 MeV gamma ray is 65.30%.

Fig. 1

NaI (Tl) detector model.

Fig. 2

Simulation model of vehicle-borne radiation environmental measurement equipment.

The car body structure is simplified to improve the efficiency of simulation calculations. Except for the chassis and tires, the car body adopts a box structure, while other scattered and lighter parts are ignored to simulate the real car to the maximum extent. Table 1 shows the detailed information, such as the structure and materials, of each component in the studied equipment.

Detailed information of the structure and materials of each component in the generalized model of vehicle-borne environmental radiation measurement equipment

Structure Material Component Density (g/cm3) Thickness (cm)
Vehicle chassis Alloy steel 3.45% C, 96.55% Fe 7.85 10.00
Vehicle tires Rubber 11.84% H, 88.16% C 0.92
Carriage Aluminium skin 100% Al 2.70 0.12
Polyurethane foam 54.4% C, 12.1% N, 33.5% O 0.04
Glass 46.0% O, 33.6% Si, 10.7% Ca, 9.7% Na 2.40 1.00
Aluminium alloy 97.43% Al, 1.55% C, 1.02% Cu 2.80 3.76
Detector crystal Sodium iodide 15.34% Na, 84.66% I 3.67
Magnesium oxide 39.7% O, 60.3% Mg 2.58 0.05
Aluminium 100% Al 2.70 0.20
Silicon dioxide 53.2% Si, 46.8% O 2.32 0.20
Construction of a near-ground reference radiation field model

A nuclear accident produces a complex radiation field composed of hundreds of radionuclides with varying distributions over time [10, 11]. Therefore, using the traditional ground radiation detection method is insufficient. In this study, the radiation field distribution in the radioactive source area of a nuclear accident is simulated using the Monte Carlo program by sampling the radiation rays and energy directions. According to the national standard GB/T17982-2018 [12], nuclear accidents are divided into time periods or phases: the early stage is the initial moment of decay (nuclear accident); the middle stage starts 1 day after the accident; and the later stage starts 100 days after the accident. The cylindrical source model is used to approximate the infinite radiation field after the nuclear accident, as shown in Fig. 3. In the early phase, radioactive substances are released into the atmosphere and dispersed in the air (Fig. 3a). These radionuclides are assumed to be uniformly distributed in the form of the plume. The cylinder has a radius of 1000 cm and a height of 2000 cm. The entire space is filled with air, and the radiation field volume is V = 6280 m3. In the mid-late nuclear phase of nuclear accidents, the radionuclides rapidly deposited in the soil after a nuclear accident, and the early radioactive materials are gradually deposited on the soil to form a radiation field of ground radioactive fallout (Fig. 3b). Similarly, the cylinder has a radius of 1000 cm and a height of 2000 cm. The bottom of the column is filled with a 1-cm-thick layer of soil [13], and the air above this surface is used as an air model. The soil density [14] is 1.7 g/cm3, and the soil mass is M = 4710 kg. Table 2 shows the air and soil content.

Fig. 3

Reference radiation field model: (a) early radiation field and (b) mid-late radiation field.

Mass fraction percentage of main components of air and soil

Air

C N O Ar
0.02 75.50 23.20 1.28

Soil

H O Al Si Fe
6.67 48.33 18.00 22.50 4.50
Source term analysis of radiation field

The consequence assessment and environmental radiation monitoring of the Fukushima accident is a complex and arduous task, and the source item is an important basis for accident classification [15, 16]. Numerous researchers of the Fukushima accident assess the total amount of radioactive materials released into the atmosphere to evaluate the environmental impact [16, 17, 18, 19, 20]. However, in daily training, nuclear emergency equipment is used in a radiation field environment formed by a single nuclide. For this reason, this study specifically selects five representative radionuclides released after nuclear accidents (Table 3) as a single-nuclide reference to study the distribution of radiation fields in the early phases of nuclear accidents [16].

Representative nuclides in single-nuclide reference radiation field

Nuclides 60Co 132Te 137Cs 131I 133Xe
Half-life 5.27 a 76.80 h 30.07 a 8.02 d 5.25 d
Gamma energy (keV) 1173.24, 1332.50 228.16 661.66 364.49 81.00

The Fukushima accident released various types of characteristic radionuclides. However, simultaneous analyses of the natural radionuclides in the environment are necessary considering the background effect. Table 4 shows the selected nuclides required for the background reference radiation field in the “ 2018 Annual Report of the National Radiation Environment” [21]. The aerosol activity is the air model for the early radiation field and the soil activity is the soil model for the mid-late stage. The complex radiation field distribution composed of hundreds of nuclides is affected by various factors, such as meteorological and geological conditions. The source term distribution of radiation field needs to be obtained from a specific time and location. Thus, this study uses the nuclide data during the Fukushima accident as the source term of the mixed reference radiation field.

Environmental background composition of the reference radiation field

Distribution Aerosol Soil
Nuclides 7Be 40K 137Cs 210Pb 137Cs 226Ra
Energy (keV) 477.61 1460.83 661.66 46.54 661.66 186.21
Background (Bq) 5.66 × 101 2.36 × 100 1.37 × 10−2 2.84 × 101 3.91 × 105 2.98 × 104

Tables 5 and 6 include the main nuclides, radioactive concentration, and activity of the mixed reference radiation field in the early and middle phases of the Fukushima accident, respectively. The radioactivity calculation of the nuclides in the reference radiation fields is shown in Eqs. (1) and (2). A1=C1×V {A_1} = {C_1} \times V A2=C2×M {A_2} = {C_2} \times M where A is the radionuclide activity, C is the radionuclide concentration, 1 and 2 are the different reference radiation fields, V is the volume of the early air reference radiation field, and M is the mass of the mid-late soil reference radiation field.

Nuclide composition of the reference radiation field in the air of early stages of a nuclear accident

Nuclides 106Ru 129Te 131I 132Te 132I 133I 134Cs 136Cs 137Cs 144Ce
Energy (keV) 511.85 459.60 364.49 228.16 667.72 529.87 604.72 818.50 661.66 696.51
Concentration (Bq/m3) 1.24 × 102 1.30 × 103 3.23 × 103 3.22 × 102 1.23 × 103 3.80 × 101 3.55 × 101 5.00 × 100 3.75 × 101 5.00 × 103
Activity (Bq) 7.79 × 105 8.16 × 106 2.53 × 107 2.02 × 106 9.26 × 106 2.39 × 105 3.99 × 105 3.17 × 104 4.26 × 105 3.14 × 107

Nuclide composition of the reference radiation field in the soil of middle and late stages of a nuclear accident

Nuclides 89Sr 131I 134Cs 136Cs 137Cs 238Pu 241Am
Energy (keV) 583.00 364.49 604.72 818.50 661.66 42.49 59.54
Concentration (Bq/kg) 1.31 × 102 3.44 × 104 8.63 × 104 1.11 × 103 8.53 × 104 2.40 × 10−3 1.21 × 10−3
Activity (Bq) 6.17 × 105 1.62 × 108 4.06 × 108 5.30 × 106 4.02 × 108 1.13 × 101 3.17 × 104
Vehicle-borne radiation system measurement

When different measuring platforms are used to assess the ground deposition of nuclides, all measurement results must be corrected to maintain data consistency from the different material structures and placement height of the measuring instrument [22, 23]. In this study, the detector is installed on the top of the car, parts of the ambient gamma rays are shielded by the body material, and the height attenuation caused by the detection equipment affects the measurement results. Therefore, CF needs to match the dose measured by the vehicle-borne radiation ambient measurement equipment with the air dose. The Monte Carlo simulation is used to study the air absorption dose rate detected by NaI (Tl) detection crystal at two positions inside and outside the vehicle. Thus, the dose CF is established, as shown in Eq. (3). CF=(D˙carD˙air) {\rm{CF}} = \left({{{{{\dot D}_{{\rm{car}}}}} \over {{{\dot D}_{{\rm{air}}}}}}} \right) where CF is the dose CF of the measured dose in the vehicle-borne radiation measurement equipment and and are the air absorption dose rates measured by the detector mounted on the car roof in the car and 1 m above the ground in the air, respectively. Under the same conditions, the simulated detector is placed on the car roof at different dose rates to verify the dose CF. In the end, the air dose rate at 1 m above the ground can be obtained by applying the dose attenuation CF to the vehicle-borne radiation dose measurement, as shown in Eq. (4). D˙1m=D˙carCF {\dot D_{1m}} = {{{{\dot D}_{{\rm{car}}}}} \over {{\rm{CF}}}} where is the ambient dose rate at 1 m above the ground and is the ambient dose rate measured by the vehicle-borne radiation measurement equipment.

Results and discussion
Ambient dose rate of the single-nuclide reference radiation field

The Monte Carlo program is used to simulate the dose rate of 60Co, 131I, 132Te, 133Xe, and 137Cs radionuclides with activity of 1 Ci in the early reference radiation field, as shown in Fig. 4. The ambient dose rate increases exponentially with the increase in radionuclide energy in the single-nuclide reference radiation field. In determination radionuclide activity, the dose rate change is only related to nuclide energy.

Fig. 4

Trend of ambient dose rate with gamma energy in the single reference radiation field.

Table 7 shows the ambient dose rates at the calibration locations. The rates range from 2.8 × 10−5 μGy/s to 1.7 × 10−2 μGy/s as measured at the height of 1 m above the ground, and the rates range from 2.5 × 10−5 μGy/s to 1.5 × 10−2 μGy/s as measured by the detector at the car roof. Under the same conditions, the trends of environmental dose rates detected by the NaI (Tl) detector inside and outside the vehicle are roughly the same. The dose rate measured outside is larger than that measured by the vehicle-borne detector because of the shielding effect of the vehicle and its contents.

Dose rates measured at the height of 1 m above the ground and at the car roof of single-nuclide radiation field (μGy/s)

Nuclides Energy (keV) Car (μGy/s) Air (μGy/s)
133Xe 81.00 2.5149 × 10−5 2.8241 × 10−5
132Te 228.16 6.8513 × 10−4 7.8970 × 10−4
131I 364.49 2.4188 × 10−3 2.7888 × 10−3
137Cs 661.66 7.0965 × 10−3 8.1073 × 10−3
60Co 1173.24, 1332.50 1.5559 × 10−2 1.7436 × 10−2
Ambient dose rate of mixed-nuclide reference radiation fields

Tables 5 and 6 show various radionuclides and activity in the mixed-nuclide reference radiation fields. For the early phase of nuclear accidents, this study considers 2 h as the step length. For the mid-late stage, 100 days as a time node is used to calculate 12 different periods (1 day, 3, 5, 7, 9, 11, 13, 20, 40, 60, 80, and 100 days) to obtain the change curve of atomic number and activity of the radioactive decay progeny, as shown in Fig. 5. The radionuclide activity in two different reference fields decreases exponentially over time. The results show that these decreases are not related to the radiation field model but only related to the decay constant.

Fig. 5

Results of radionuclide activity over time: (a) early radiation field and (b) mid-late radiation field.

Figures 6a and 6b show the trend of changes in the environmental dose rate measured by the NaI (Tl) detector inside and outside of the vehicle, respectively. In this study, the early reference radiation field is the air body source, and the air absorption dose rate is evaluated to only consider the shielding effect of the car body. Thus, the dose rates inside and outside the vehicle slightly differ. By comparison, the mid-late reference radiation field is the soil body source, and the NaI (Tl) detector is placed on the roof, which essentially increases the source range. Thus, the NaI (Tl) detector receives a large difference in the ambient dose rate inside and outside the vehicle. The dose rate changes less and tends to flatten due to the long half-life and slow decay of certain radionuclides, such as 137Cs.

Fig. 6

Ambient dose rates of the mixed-nuclide radiation field varying with time: (a) early radiation field and (b) mid-late radiation field.

Correction factor of the single-nuclide reference radiation field

By convention, the ICRP 1996 [9] stipulates measurement of ambient dose rate at 1 m height to ensure consistent results when using different vehicle radiation measurement systems. In this study, only the shielding effect of the car body and the influence of the detector distance from the ground are considered to evaluate the air absorption dose rate. The CF of the vehicle-borne radiation environment measurement equipment under different radiation fields can be obtained using Eq. (3).

The ambient dose rate of the single-nuclide reference radiation field is simulated using the Monte Carlo program. Table 8 shows the results of the dose CF obtained by Eq. (3). The mean dose CF is 0.8786, and the error is within 2%. The results show that energy change negligibly affects the dose CF.

Dose CF of single-nuclide reference radiation field

Nuclides 133Xe 132Te 131I 137Cs 60Co
CF 0.8905 0.8676 0.8673 0.8753 0.8923

CF, correction factor.

Correction factor of the mixed-nuclide reference radiation field

The dose CF of the mixed-nuclide reference radiation field using the vehicle-borne environmental radiation measurement equipment is listed in Tables 9 and 10. The early reference field CF is regarded as a fixed value of 0.8830, and the mid-late mixed-nuclide reference field dose CF fluctuates at approximately 0.6711. The dose rate measured by the vehicle-borne equipment is applied to 1 m above the ground in air using Eq. (4).

Dose CF of early mixed-nuclide reference radiation field

Time (h) 0 2 4 6 8 10 12 14 16 18 20 22 24
CF 0.883 0.883 0.883 0.884 0.884 0.884 0.884 0.884 0.884 0.884 0.884 0.883 0.883

CF, correction factor.

Dose CF of mid-late mixed-nuclide reference radiation field

Time (d) 1 3 5 7 9 11 13 20 40 60 80 100
CF 0.669 0.669 0.670 0.670 0.670 0.671 0.671 0.672 0.672 0.672 0.672 0.673

CF, correction factor.

Effects of detector position

The dose rate measurements are repeated at several locations with different ambient dose levels and detector installation locations to verify the validity of the dose CF. Table 11 shows the average results of the detector placed at different positions. The data are divided into dose rate and CF according to Eq. (3).

Detector located at different locations of the dose CF

Central place Car (μGy/s) Air (μGy/s) CF Error
97.5, 160, 160 7.0442 × 10−3 7.8729 × 10−3 0.8947 0.8813 + 0.015
97.5, 195, 160 7.0905 × 10−3 7.9737 × 10−3 0.8892 0.8813 + 0.009
97.5, 230, 160 7.0965 × 10−3 8.1073 × 10−3 0.8753 0.8813 − 0.007
47.5, 240, 160 7.0088 × 10−3 8.0049 × 10−3 0.8756 0.8813 − 0.006
147.5, 240, 160 7.0541 × 10−3 8.0180 × 10−3 0.8798 0.8813 − 0.002
97.5, 265, 160 7.1126 × 10−3 8.1010 × 10−3 0.8780 0.8813 − 0.004
97.5, 300, 160 6.8867 × 10−3 7.8778 × 10−3 0.8742 0.8813 − 0.008
97.5, 335, 160 7.1724 × 10−3 8.2174 × 10−3 0.8728 0.8813 − 0.010
47.5, 370, 160 7.0047 × 10−3 7.9291 × 10−3 0.8834 0.8813 + 0.002
97.5, 370, 160 6.9351 × 10−3 7.8866 × 10−3 0.8794 0.8813 − 0.002
147.5, 370, 160 7.2050 × 10−3 8.1576 × 10−3 0.8832 0.8813 + 0.002
97.5, 405, 160 6.9543 × 10−3 7.8149 × 10−3 0.8899 0.8813 + 0.010

CF, correction factor.

Conclusion

This study establishes a near-ground reference radiation field based on typical radionuclides of the Fukushima accident as the source term data. The Monte Carlo program simulation is used to study the environmental dose of vehicle-borne environmental radiation measurement system. The dose CF is established by simulating the ambient dose rate of the Cheetah Land Cruiser LBA6482LQ4 with a 4 L NaI (Tl) detector in the early and mid-late reference radiation fields of nuclear accidents. The main conclusions are as follows:

In the early reference radiation field, the radionuclides are uniformly distributed in the air, the ambient dose rate is only related to the vehicle material structure, and the average dose CF is 0.8813.

In the mid-late reference radiation field, the ambient dose rate is highly correlated with the body material structure and the height of the detector from the surface under the assumption that radio-nuclides all sank to the surface to form soil body source. The average value of the environmental dose CF is 0.6711.

The results of the ambient dose CF are independent of the detector position, but this case is only in the established reference radiation field. This paper presents the study of the dose correction method for the monitoring process. The dynamic measurement of this process will be presented in a subsequent paper.

eISSN:
1508-5791
Language:
English
Publication timeframe:
4 times per year
Journal Subjects:
Chemistry, Nuclear Chemistry, Physics, Astronomy and Astrophysics, other