WO2024057514A1 - 放射能分析装置及び放射能分析方法 - Google Patents
放射能分析装置及び放射能分析方法 Download PDFInfo
- Publication number
- WO2024057514A1 WO2024057514A1 PCT/JP2022/034696 JP2022034696W WO2024057514A1 WO 2024057514 A1 WO2024057514 A1 WO 2024057514A1 JP 2022034696 W JP2022034696 W JP 2022034696W WO 2024057514 A1 WO2024057514 A1 WO 2024057514A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- radiation
- signal
- inverse problem
- radioactivity
- problem calculation
- Prior art date
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 65
- 230000005855 radiation Effects 0.000 claims abstract description 239
- 238000004364 calculation method Methods 0.000 claims abstract description 68
- 238000001514 detection method Methods 0.000 claims abstract description 68
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 58
- 238000005316 response function Methods 0.000 claims abstract description 56
- 238000005259 measurement Methods 0.000 claims abstract description 54
- 239000000126 substance Substances 0.000 claims description 17
- 230000010365 information processing Effects 0.000 description 21
- 238000010586 diagram Methods 0.000 description 20
- 239000011159 matrix material Substances 0.000 description 16
- 238000000034 method Methods 0.000 description 16
- 230000006870 function Effects 0.000 description 14
- 238000012545 processing Methods 0.000 description 11
- 239000004065 semiconductor Substances 0.000 description 10
- 239000013078 crystal Substances 0.000 description 8
- 238000010801 machine learning Methods 0.000 description 8
- 230000005250 beta ray Effects 0.000 description 6
- 230000005260 alpha ray Effects 0.000 description 5
- 230000005251 gamma ray Effects 0.000 description 5
- 230000006399 behavior Effects 0.000 description 4
- 230000003321 amplification Effects 0.000 description 3
- 230000005255 beta decay Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000003199 nucleic acid amplification method Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000000691 measurement method Methods 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 239000000941 radioactive substance Substances 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- OYEHPCDNVJXUIW-FTXFMUIASA-N 239Pu Chemical compound [239Pu] OYEHPCDNVJXUIW-FTXFMUIASA-N 0.000 description 1
- GUTLYIVDDKVIGB-OUBTZVSYSA-N Cobalt-60 Chemical compound [60Co] GUTLYIVDDKVIGB-OUBTZVSYSA-N 0.000 description 1
- LXQXZNRPTYVCNG-YPZZEJLDSA-N americium-241 Chemical compound [241Am] LXQXZNRPTYVCNG-YPZZEJLDSA-N 0.000 description 1
- 230000009118 appropriate response Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- NIWWFAAXEMMFMS-OIOBTWANSA-N curium-244 Chemical compound [244Cm] NIWWFAAXEMMFMS-OIOBTWANSA-N 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 229910052732 germanium Inorganic materials 0.000 description 1
- GNPVGFCGXDBREM-UHFFFAOYSA-N germanium atom Chemical compound [Ge] GNPVGFCGXDBREM-UHFFFAOYSA-N 0.000 description 1
- 239000012212 insulator Substances 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000004020 luminiscence type Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 229910044991 metal oxide Inorganic materials 0.000 description 1
- 150000004706 metal oxides Chemical class 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- OYEHPCDNVJXUIW-VENIDDJXSA-N plutonium-238 Chemical compound [238Pu] OYEHPCDNVJXUIW-VENIDDJXSA-N 0.000 description 1
- 230000005258 radioactive decay Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 229910052716 thallium Inorganic materials 0.000 description 1
- BKVIYDNLLOSFOA-UHFFFAOYSA-N thallium Chemical compound [Tl] BKVIYDNLLOSFOA-UHFFFAOYSA-N 0.000 description 1
- 229910052727 yttrium Inorganic materials 0.000 description 1
- VWQVUPCCIRVNHF-UHFFFAOYSA-N yttrium atom Chemical compound [Y] VWQVUPCCIRVNHF-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/167—Measuring radioactive content of objects, e.g. contamination
Definitions
- the present disclosure relates to a radioactivity analysis device and a radioactivity analysis method.
- a radiation detector measures the energy distribution of radiation emitted from the measurement target, and based on the measurement results, the radionuclide that is the radioactive substance contained in the measurement target is determined. identification and calculation of the radioactivity intensity [Bq (becquerel)] of the radionuclide.
- Non-Patent Document 1 describes a mathematical method for calculating the energy spectrum of incident radiation using a response function of a radiation detector based on the output wave height distribution (i.e., pulse height distribution) of a pulse signal output from a radiation detector. It explains the unfolding method. Furthermore, Patent Document 1 describes a method of improving the accuracy of estimation results by increasing the number of response functions used in inverse problem calculations.
- the radioactivity analysis method described above has a problem in that the accuracy of radioactivity analysis is low when the measurement target contains radionuclides that emit multiple types of radiation.
- the present disclosure has been made to solve the above problems, and aims to provide a radioactivity analysis device and a radioactivity analysis method that have high accuracy in radioactivity analysis.
- the radioactivity analyzer of the present disclosure includes a radiation detection unit that detects radiation emitted from a measurement target and outputs a first signal based on the radiation, and indicates the type of the radiation based on the first signal.
- a radiation line type determination unit that determines the radiation line type and outputs a second signal including information indicating the radiation line type;
- an inverse problem calculation algorithm selection unit that selects an inverse problem calculation algorithm based on the inverse problem calculation algorithm; and an inverse problem calculation using the selected inverse problem calculation algorithm and a response function of the radiation detection unit stored in advance in a storage device.
- the method further includes an inverse problem calculation unit that identifies a radionuclide included in the measurement target and calculates the radioactivity intensity of the radionuclide from the first signal and the second signal.
- the radioactivity analysis method of the present disclosure is a method executed by a computer, in which the first signal is output from a radiation detection unit that detects radiation emitted from a measurement target and outputs a first signal based on the radiation. a step of determining a radiation line type indicating the type of radiation based on the first signal and outputting a second signal including information indicating the radiation line type, and storing the information in advance in a storage device.
- radioactivity analysis device and radioactivity analysis method of the present disclosure highly accurate radioactivity analysis can be performed.
- FIG. 1 is a block diagram schematically showing the configuration of a radioactivity analyzer according to Embodiment 1.
- FIG. FIG. 3 is a diagram showing the wave height (that is, energy) distribution of pulse signals output from the detector of the radiation detection unit for each radiation line type indicating the type of radiation.
- 1 is a diagram showing a hardware configuration of a radioactivity analyzer according to Embodiment 1.
- FIG. 3 is a flowchart showing a radioactivity analysis method according to Embodiment 1.
- FIG. FIG. 2 is a block diagram schematically showing the configuration of a radioactivity analyzer according to Embodiment 2.
- FIG. FIG. 3 is a block diagram schematically showing the configuration of a radioactivity analyzer according to Embodiment 3.
- FIG. 3 is a diagram showing the rise and decay characteristics of the amplitude of a pulse signal output from a detector of a radiation detection unit for each type of radiation line.
- FIG. 3 is a block diagram schematically showing the configuration of a radioactivity analyzer according to Embodiment 4.
- FIG. (A) to (C) are schematic diagrams showing trajectory shapes in a detector of a radiation detection unit for each radiation line type.
- FIG. 3 is a block diagram schematically showing the configuration of a radioactivity analyzer according to Embodiment 5.
- FIG. FIG. 3 is a diagram showing the relationship between the brightness value and the spread of a locus shape in a position detection type detector of a radiation detection unit.
- FIG. 1 is a block diagram schematically showing the configuration of a radioactivity analyzer 1 according to the first embodiment.
- the radioactivity analyzer 1 detects radiation emitted from a measurement target (also referred to as “measurement material” or “measurement object”) 70, and outputs a first signal D1 based on the detected radiation. section 10, and an information processing section that identifies radioactive substances (also referred to as “radionuclides” or “nuclides”) contained in the measurement target 70 and calculates the radioactivity intensity [Bq (Becquerel)] of the radionuclides. 1a.
- the information processing unit 1a is, for example, a computer serving as an information processing device, and is capable of implementing the radioactivity analysis method according to the first embodiment.
- the information processing section 1a of the radioactivity analyzer 1 includes a radiation line type determination section 20, an inverse problem calculation algorithm selection section 30, and an inverse problem calculation section 50.
- the inverse problem calculation algorithm selection unit 30 selects the radiation radiation indicated by the second signal D2 from a plurality of inverse problem calculation algorithms (for example, stored as the inverse problem calculation algorithm database 31) stored in advance in the storage device.
- the inverse problem calculation algorithm is selected based on the radiation line type indicating the type.
- the inverse problem calculation unit 50 performs inverse problem calculation using the selected inverse problem calculation algorithm and the response function of the radiation detection unit 10 that is stored in advance in the storage device (for example, stored as the response function database 41).
- the radionuclide contained in the measurement target 70 is identified and the radioactivity intensity [Bq] of the radionuclide is calculated from the first signal D1 and the second signal D2.
- the inverse problem calculation section 50 identifies each of the plurality of radionuclides and determines the radioactivity strength [Bq (becquerel)] of each of the plurality of radionuclides. ] can be calculated.
- the inverse problem calculation unit 50 calculates the radiation intensity (for example, radiation fluence (flux) [1/(cm 2 ⁇ s)]) may be calculated.
- the calculated radioactivity may be expressed as radioactivity concentration [Bq/g] (or [Bq/cm 3 ], etc.), which is the intensity of radioactivity per unit weight (or unit volume) of the measurement target 70. good.
- the storage device in which the inverse problem calculation algorithm database (also referred to as "inverse problem calculation algorithm DB") 31 is stored and the storage device in which the response function database (also referred to as "response function DB") 41 is stored are as follows. They may be the same device or different devices.
- the storage device in which the inverse problem calculation algorithm DB 31 is stored and the storage device in which the response function DB 41 is stored may be part of the radioactivity analysis device 1, but may be other devices that can communicate with the radioactivity analysis device 1. It may be part of a device (eg, a computer on a network).
- the radioactivity analyzer 1 is connected to a display unit 60 as an information output device that presents the calculation results of the inverse problem calculation unit 50 to the user.
- the display section 60 may be a part of the radioactivity analyzer 1.
- the display unit 60 is, for example, a display device such as a liquid crystal display that displays images.
- another information output device for example, a printer, etc. may be provided.
- the radiation detection unit 10 includes a detector 11 that detects radiation emitted from a measurement target 70, an amplifier 12 that amplifies the output pulse output from the detector 11, and a waveform that shapes the output pulse output from the amplifier 12. It has a shaper 13 and a pulse height analysis section 14 that generates a pulse height distribution of output pulses output from the waveform shaper 13.
- the detector 11 has a function of outputting signal strength (for example, pulse amplitude) according to the energy given to the detector by radiation.
- Examples of the detector 11 that are widely used include a scintillation detector that combines a scintillator crystal and a photomultiplier, a scintillation detector that combines a scintillator crystal and a semiconductor photodetector, a semiconductor detector, or a gas detector. etc. are known.
- the semiconductor detector is, for example, a germanium semiconductor detector.
- a gas detector detects the number of electrons (i.e., a value proportional to the amount of energy loss of the charged particles) produced by the collision of charged particles with electrons in the gas molecules as radiation passing through the gas in the container. This is a detector that measures .
- Gas detectors include Geiger-Mueller (GM) counters and detectors filled with gas and using an ionization chamber.
- the amplifier 12 is a circuit that has a function of amplifying the amplitude of the output pulse, which is the signal output from the detector 11.
- the amplifier 12 is used to amplify the output pulse to a signal strength (for example, a sufficiently large pulse amplitude) sufficient for detection by the pulse height analysis section 14 and the radiation type determination section 20 at the subsequent stage. If the output of the detector 11 is sufficiently large, there is no need to provide the amplifier 12.
- the waveform shaper 13 is a circuit that has the function of removing noise from the signal output from the amplifier 12 or the detector 11 and adjusting the amplification factor.
- the waveform shaper 13 is formed of an analog circuit, it is generally configured by combining an integrating circuit and a differentiating circuit.
- the waveform shaper 13 removes noise and adjusts the amplification factor, and forms the output pulse into a waveform that allows the pulse height analysis section 14 to easily detect the pulse height.
- the waveform shaper 13 When the waveform shaper 13 is formed by a digital circuit, it digitally converts the signal output from the previous stage by analog-to-digital conversion, removes noise and adjusts the amplification factor by a digital filter, and converts the signal into a digital waveform-shaped signal. Perform pulse processing.
- the pulse height analyzer 14 creates a frequency distribution of the peak values of the output pulses output from the waveform shaper 13 or a frequency distribution of the time widths of the output pulses.
- the pulse height analyzer 14 acquires the wave height of the output pulse by combining a circuit that holds the peak value of the output pulse and analog-to-digital conversion.
- the wave height of the output pulse is obtained by combining a circuit that converts the rise and fall times of the output pulse into voltage and analog-to-digital conversion.
- the pulse height analyzer 14 obtains the wave height of the output pulse from the amplitude or time width of the output pulse processed by the digital filter.
- the pulse height analysis unit 14 creates a frequency distribution from the wave heights of the acquired output pulses.
- FIG. 2 is a diagram showing the wave height (that is, energy) distribution of the pulse signal output from the radiation detection unit 10 for each radiation line type.
- the measured wave height distribution may be converted into an energy distribution.
- the radiation type is one or more of ⁇ rays, ⁇ rays, and ⁇ rays.
- the radiation type determining unit 20 has a function of detecting the type of radiation (that is, the radiation type) emitted from the measurement object 70 from the pulse height distribution or energy distribution obtained by the pulse height analysis unit 14.
- alpha rays ( ⁇ rays) detected by the radiation detection unit 10 are generally emitted from alpha ray emitting nuclides and have energy of 4 MeV or more.
- plutonium-238 ( 238 Pu) emits ⁇ -rays with an energy of 5.50 MeV and ⁇ -rays with an energy of 5.46 MeV.
- Plutonium-239 ( 239 Pu) emits ⁇ -rays with an energy of 5.16 MeV.
- Americium-241 ( 241 Am) emits ⁇ -rays with an energy of 5.49 MeV.
- Curium-244 ( 244 Cm) emits alpha rays with an energy of 5.95 MeV.
- ⁇ rays having high energy are detected by the detector 11 so that a pulse height distribution is formed in the high energy region, as shown by the solid line in FIG.
- the energy of the alpha rays emitted from the measurement object 70 has an asymmetrical peak structure due to energy loss in the air, but exhibits a characteristic in which a peak structure appears on the wave height distribution.
- the radiation type determination unit 20 can determine whether the measurement target 70 contains an ⁇ -ray emitting nuclide.
- gamma rays ( ⁇ rays) emitted from radionuclides generally have energy specific to the nuclide, and the wave height distribution of ⁇ rays has a peak structure close to a normal distribution. It has a continuous component that reflects the process of losing energy due to scattering inside or around the detector 11.
- the energy of gamma rays of general radionuclides is 3 MeV or less, which is lower than the peak of the pulse height distribution of alpha rays.
- thallium (Tl) emits gamma rays with an energy of 2.6 MeV.
- Yttrium (Y) emits gamma rays with an energy of 1.8 MeV.
- Cobalt-60 ( 60 Co) emits two gamma rays with energies of 1.3 MeV and 1.1 MeV.
- the energy of ⁇ -rays emitted from ⁇ -ray-emitting nuclides is lower than the energy of ⁇ -rays emitted from ⁇ -ray-emitting nuclides, and the shape of the peak structure is also different (it has a steep peak structure and a continuous component). Therefore, the radiation type determining unit 20 can determine whether the measurement target 70 contains a gamma-ray emitting nuclide from the pulse height distribution.
- Beta rays ( ⁇ rays) emitted from radionuclides are emitted during the process of beta decay, which is a type of radioactive decay of atomic nuclei.
- beta decay electrons (ie, beta rays) and antielectron neutrinos are released.
- beta decay the electron and antielectron neutrino share energy, so the pulse height distribution of the ⁇ ray, which is the energy of the emitted electron (i.e. ⁇ ray), is as shown by the dashed line in Figure 2. , a continuous distribution.
- the pulse height distribution of ⁇ -rays depends on the structure of the detector 11, a peak structure appears in a low energy region, but unlike ⁇ -rays and ⁇ -rays, it does not have a peak structure in a specific energy region. Similar to the pulse height distribution of gamma rays, the pulse height distribution of ⁇ rays is a continuous distribution in the low energy region. Furthermore, in the low energy region, the pulse height distribution of ⁇ rays is larger than the pulse height distribution of ⁇ rays. Since the shape of the wave height distribution of ⁇ -rays is different from that of other radiation types, the radiation type determining unit 20 can determine whether the measurement target 70 contains a ⁇ -ray emitting nuclide from the wave height distribution.
- the characteristics of each radiation line type in the wave height distribution can be detected, for example, by analytically detecting the above characteristics, by verifying the degree of agreement with a previously assumed wave height distribution, or by detecting the characteristics of each radiation line type based on a previously assumed wave height distribution. This can be achieved by either a method of determining the radiation line type by machine learning or a method of learning to determine the radiation line type by unsupervised machine learning.
- the wave height distribution assumed in advance can be created by calculating the wave height distribution under various conditions. The calculation of the wave height distribution includes, for example, calculation based on radiation behavior analysis.
- the function of performing detection using the pulse height distribution of the pulse height analysis unit 14 may be performed by the radiation type determination unit 20. In that case, it is possible not to include the pulse height analysis section 14.
- the radiation type output from the detector 11 by determining the radiation type output from the detector 11 and selecting an inverse problem calculation algorithm to be used according to the detected radiation type, it is possible to detect a plurality of radiation types in the measurement object 70. Even if it is included, it is possible to determine the radioactivity intensity of each radiation type.
- the characteristics of each radiation line type in the wave height distribution can be detected, for example, by analytically detecting the above characteristics, by verifying the degree of agreement with a previously assumed wave height distribution, or by detecting the characteristics of each radiation line type based on a previously assumed wave height distribution. This can be achieved by learning to determine radiation line types using machine learning or unsupervised machine learning.
- the wave height distribution assumed in advance can be created by calculating the wave height distribution under various conditions. The calculation of the wave height distribution includes, for example, calculation based on radiation behavior analysis.
- the inverse problem calculation unit 50 identifies the radionuclide included in the measurement target 70 from the wave height distribution and calculates the radioactivity intensity [Bq] of the radionuclide.
- the radioactivity strength [Bq] may be expressed as radioactivity concentration.
- the nuclide of the radionuclide contained in the measurement object 70, the radioactivity intensity [Bq] of the radionuclide, and the radioactivity concentration are known, and the measurement is performed in a state where the wave height distribution is obtained as a forward problem. be done.
- R(E) is a response function
- S(E) is the intensity of radioactivity (or radioactivity concentration) for each energy of the measurement target 70.
- algorithms for performing inverse problem calculations include inverse matrix calculation, pseudo-inverse matrix calculation, deconvolution, deconvolution, and successive approximation.
- the response function can be prepared by using, for example, a method such as radiation behavior analysis that can perform calculations that reflect the physical process of radiation.
- M(E) and S(E) can be created by using radiation behavior analysis in the same manner as above.
- S(E) and M(E ) can also be used to solve inverse problem operations.
- ⁇ Inverse problem calculation algorithm selection unit 30> Next, the inverse problem calculation algorithm selection unit 30 will be explained. A plurality of inverse problem calculation algorithms that can perform inverse problem calculations as described above are prepared in advance. In some cases, it may not be possible to obtain a solution due to the nature of the matrix that is the response function R(E). The response function R(E) differs depending on the radiation type.
- the ⁇ -ray response function R(E) is shown below as equation (2).
- the following response function R(E) is represented by a matrix whose diagonal elements have values greater than zero.
- the ⁇ -ray response function R(E) is shown below as equation (3).
- the response function R(E) below is represented by a matrix whose diagonal elements are zero.
- an analytical calculation algorithm such as an inverse matrix calculation, a pseudo-inverse matrix calculation, or a successive approximation method can be used.
- the response function for a radiation type that is attenuated in the air, such as ⁇ -rays is a matrix that does not have diagonal components, as shown in equation (3).
- the response function is represented by a matrix without diagonal elements, an inverse matrix operation cannot be used because an inverse matrix does not exist.
- pseudo-inverse matrix calculations may have problems with calculation accuracy, such as the existence of solutions that are less than or equal to 0.
- the response function can be created using the estimated spectrum from the assumed radionuclide.
- the condition is that the response function is not a square matrix, and there is no inverse matrix, so inverse problem calculation using a network using machine learning may be able to estimate the optimal result rather than analytical inverse problem calculation.
- the inverse problem calculation algorithm selection unit 30 optimally selects one or more calculation algorithms based on the detected radiation line type so that the solution converges or the calculation speed becomes faster.
- the inverse problem calculation algorithm DB 31 stores a plurality of calculation algorithms that can be used in the inverse problem calculation unit 50 as described above.
- the response function database 41 has parameters of a response matrix or machine learning network structure used for inverse problem calculation.
- matrices such as equations (2) and (3) are held. Further, it may have information regarding a network for performing inverse problem calculation using machine learning, which is a method other than analytical calculation.
- FIG. 3 is a diagram showing an example of the hardware configuration of the radioactivity analyzer 1 according to the first embodiment.
- the hardware configuration of FIG. 3 is also applied to the radioactivity analyzers according to the second to fifth embodiments.
- the radioactivity analyzer 1 includes a radiation detection section 10 and an information processing section 1a.
- the information processing unit 1a includes a processor 101 such as a CPU (Central Processing Unit), a memory 102 that is a volatile storage device, and a nonvolatile storage device 103 such as a hard disk drive (HDD) or a solid state drive (SSD). have.
- the memory 102 is, for example, a semiconductor memory such as a RAM (Random Access Memory).
- the storage device 103 may store the database shown in FIG.
- Each function of the information processing unit 1a is realized by, for example, a processing circuit.
- the processing circuitry may be dedicated hardware or may be a processor 101 executing a program stored in memory 102.
- the processor 101 may be any one of a processing device, an arithmetic device, a microprocessor, a microcomputer, and a DSP (Digital Signal Processor).
- the processing circuit is dedicated hardware, the processing circuit is, for example, a single circuit, a composite circuit, a programmed processor, a parallel programmed processor, an ASIC (Application Specific Integrated Circuit), an FPGA (Field-Programmable Gate Arra), etc. y ), or a combination of any of these.
- ASIC Application Specific Integrated Circuit
- FPGA Field-Programmable Gate Arra
- the radioactivity analysis program executed by the information processing unit 1a is realized by software, firmware, or a combination of software and firmware.
- the radioactivity analysis program is installed in the information processing unit 1a via a network or from a recording medium.
- Software and firmware are written as programs and stored in memory 102.
- the processor 101 can realize the functions of each section shown in FIG. 1 by reading and executing the display control program stored in the memory 102.
- the information processing unit 1a may be partially realized by dedicated hardware and partially realized by software or firmware. In this way, the processing circuit can implement each of the above-mentioned functions using hardware, software, firmware, or any combination thereof.
- FIG. 4 is a flowchart showing the radioactivity analysis method according to the first embodiment.
- This radioactivity analysis method is executed by a computer serving as the information processing section 1a.
- This radioactivity analysis method includes a step S1 of receiving a first signal D1 from a radiation detection section 10 that detects radiation emitted from a measurement target 70 and outputs a first signal D1 based on the radiation; A step S2 of determining a radiation line type indicating the type of radiation based on the information and outputting a second signal D2 including information indicating the radiation line type, and a step S2 of determining a radiation line type indicating the type of radiation based on Step S3 of selecting an inverse problem calculation algorithm based on the radiation type, and executing the inverse problem calculation using the selected inverse problem calculation algorithm and the response function of the radiation detection unit 10 stored in advance in the storage device.
- the method includes a step S4 of identifying a radionuclide included in the measurement target 70 and calculating the radioactivity intensity of the radionuclide from the first signal D1 and the second signal D2. . Further, in step S4, the intensity of radiation (for example, the fluence of radiation) may be calculated.
- the intensity of radiation for example, the fluence of radiation
- FIG. 5 is a block diagram schematically showing the configuration of the radioactivity analyzer 2 according to the second embodiment.
- components that are the same as or correspond to those shown in FIG. 1 are given the same reference numerals as those shown in FIG.
- the radioactivity analysis device 2 according to the second embodiment differs from the radioactivity analysis device 1 according to the first embodiment in that the information processing section 2a includes a response function selection section 40.
- the radioactivity analyzer 2 includes a radiation detection unit 10 that detects radiation emitted from a measurement target 70 and outputs a first signal D1 based on the detected radiation, and a radioactivity detector 10 that detects radiation emitted from a measurement target 70 and outputs a first signal D1 based on the detected radiation. It has an information processing unit 2a that specifies the nuclide and calculates the radioactivity strength [Bq] of the radionuclide.
- the information processing unit 2a is, for example, a computer, and is capable of implementing the radioactivity analysis method according to the second embodiment.
- a plurality of response functions are stored in advance in a storage device (for example, the storage device 103 in FIG. 3) as the response function DB 41.
- the response function selection unit 40 selects a response function to be used from a plurality of response functions stored in the response function DB 41 based on the second signal D2 indicating the radiation type determined by the radiation type determination unit 20. select.
- the inverse problem calculation unit 50 performs the first inverse problem calculation by using the inverse problem calculation algorithm selected by the inverse problem calculation algorithm selection unit 30 and the response function selected by the response function selection unit 40.
- the radionuclide contained in the measurement target 70 is identified and the radioactivity intensity [Bq] of the radionuclide is calculated from the signal D1 and the second signal D2. Further, the inverse problem calculation unit 50 may calculate the intensity of radiation from the first signal D1 and the second signal D2 by performing an inverse problem calculation.
- a suitable inverse problem calculation algorithm is selected from a plurality of inverse problem calculation algorithms, and a suitable inverse problem calculation algorithm is selected from a plurality of response functions.
- the second embodiment is the same as the first embodiment.
- FIG. 6 is a block diagram schematically showing the configuration of the radioactivity analyzer 3 according to the third embodiment.
- the radioactivity analyzer 3 according to the third embodiment has a radiation detection unit 10b including a waveform discriminator 15 that discriminates a detection signal for each type of radiation detected by the detector 11 and outputs a discrimination signal. This is different from the radioactivity analyzer 1 according to the first embodiment.
- the waveform discriminator 15 is a circuit that allows the pulse shapes of the output of the detector 11 and the output signal output from the amplifier 12 to differ depending on the type of radiation line.
- the radiation line type determining unit 20b determines the radiation line type based on the first signal D1 including the discrimination signal, and outputs a second signal D2 including information indicating the determined radiation line type.
- the radioactivity analyzer 3 includes a radiation detection unit 10b that detects radiation emitted from the measurement target 70 and outputs a first signal D1 including a discrimination signal indicating the wave height distribution of the radiation, It has an information processing unit 3a that specifies the nuclide and calculates the radioactivity strength [Bq] of the radionuclide.
- the information processing unit 3a is, for example, a computer, and is capable of implementing the radioactivity analysis method according to the third embodiment.
- FIG. 7 is a diagram showing the rise and attenuation characteristics of the amplitude of the pulse signal output from the detector 11 of the radiation detection section 10b for each radiation line type.
- a semiconductor is used as the substance that interacts with radiation in the detector 11
- electrons and holes are generated inside the semiconductor by the passage of radiation, but the mobility of electrons and the mobility of holes are different. is different.
- the pulse rises in a gentle shape see Figure 7). (indicated by a solid line).
- the waveform discriminator 15 can detect a specific radiation line type with high accuracy by discriminating the radiation line type using such a difference in waveform shape.
- the inverse problem calculation unit 50 solves the first problem by executing the inverse problem calculation using the inverse problem calculation algorithm selected by the inverse problem calculation algorithm selection unit 30 and the response function. From the signal D1 and the second signal D2, the radionuclide contained in the measurement target 70 is identified and the radioactivity intensity [Bq] of the radionuclide is calculated. Further, the inverse problem calculation unit 50 may calculate the intensity of radiation from the first signal D1 and the second signal D2 by performing an inverse problem calculation.
- the radioactivity analysis device 3 and the radioactivity analysis method according to the third embodiment by selecting and using a suitable inverse problem calculation algorithm from a plurality of inverse problem calculation algorithms, It is possible to improve the accuracy of problem calculations.
- the waveform discriminator 15 determines a specific radiation line type by the waveform discriminator 15 and combining the determination result with detection from the wave height distribution by the radiation line type determination unit 20, detection accuracy can be increased. Further, by correlating the detection result of a specific radiation type with the pulse height distribution, it is possible to extract the pulse height distribution of the specific radiation type, thereby separating the pulse height distribution and solving the inverse problem calculation.
- Embodiment 3 is the same as Embodiment 1 or 2.
- FIG. 8 is a block diagram schematically showing the configuration of the radioactivity analyzer 4 according to the fourth embodiment.
- the radioactivity analysis device 4 according to the fourth embodiment has two points: the detector 11c of the radiation detection section 10c has a function of detecting a position where energy is imparted by radiation, and the radiation detection section 10c has a trajectory shape detection section 16. This is different from the radioactivity analyzer 1 according to the first embodiment.
- the radioactivity analyzer 4 includes a radiation detection unit 10c that detects radiation emitted from the measurement target 70 and outputs a first signal D1 based on the detected radiation and a trajectory signal indicating the trajectory shape of the radiation, and
- the information processing unit 4a specifies the radionuclides included in the radionuclide 70 and calculates the radioactivity strength [Bq] of the radionuclides.
- the information processing unit 4a is, for example, a computer, and is capable of implementing the radioactivity analysis method according to the fourth embodiment.
- the radiation detection unit 10c includes a detector 11c having a substance that interacts with radiation, and a circuit that detects the position of a locus to which radiation imparts energy within the interacting substance and outputs a locus shape signal indicating the shape of the locus. It has a trajectory shape detection section 16.
- the radiation line type determining unit 20c determines the radiation line type based on the first signal D1 including the trajectory shape signal, and outputs a second signal D2 including information indicating the determined radiation line type.
- the detector 11c has a function of detecting a position energized by radiation.
- the detector 11c is a device in which a plurality of scintillator crystals are arranged in an array and a photodetector is placed on each of the plurality of scintillator crystals arranged in an array (for example, see Patent Document 2).
- the detector 11c includes a plurality of regularly arranged scintillator crystals and a camera (for example, a CMOS (Complementary Metal Oxide Semiconductor) camera element, SOI (Silicon-On-Insulator)) that photographs the light generated by the scintillator crystals. An image sensor, etc.) may also be used.
- CMOS Complementary Metal Oxide Semiconductor
- SOI Silicon-On-Insulator
- the detector 11c may be a device having a semiconductor radiation detector (for example, a Ge crystal) and a plurality of strip electrodes regularly arranged at intervals so as to sandwich the semiconductor radiation detector (for example, a Ge crystal). good. Further, the detector 11c may be a device that detects a position energized by radiation using a gas-type radiation detector having a plurality of wires or a plurality of regularly arranged gas-type radiation detectors.
- FIGS. 9(A) to 9(C) are schematic diagrams showing the locus shapes in the detector 11c of the radiation detection unit 10c for each radiation line type.
- FIGS. 9A to 9C show trajectory shapes within the detector 11c depending on the radiation type (when viewed from the side and when viewed from above).
- radiation having a short range and giving large energy such as ⁇ -rays
- the ⁇ ray has a long range and passes through the detector in a linear manner, so it has an elongated shape.
- the ⁇ -rays react at one point, so they have a point-like shape with no spread.
- the radioactivity analysis device 4 and the radioactivity analysis method according to the fourth embodiment by selecting and using a suitable inverse problem calculation algorithm from a plurality of inverse problem calculation algorithms, It is possible to improve the accuracy of problem calculations.
- detection accuracy can be increased by determining a specific radiation line type by the trajectory shape detection unit 16 and combining the determination result with detection from the pulse height distribution by the radiation line type determination unit 20. Further, by correlating the detection result of a specific radiation type with the pulse height distribution and extracting the pulse height distribution of the specific radiation type, it is possible to separate the pulse height distribution and solve the inverse problem calculation.
- Embodiment 4 is the same as Embodiment 1 or 2.
- FIG. 10 is a block diagram schematically showing the configuration of the radioactivity analyzer 5 according to the fifth embodiment.
- components that are the same as or correspond to those shown in FIG. 1 are given the same reference numerals as those shown in FIG.
- the radioactivity analysis device 5 according to the fifth embodiment differs from the radioactivity analysis device 1 according to the first embodiment in that the radiation detection section 10d is constituted by a position detection type detector 11d.
- the radioactivity analyzer 5 detects radiation emitted from the measurement target 70, and outputs a signal D5 indicating the brightness value and trajectory spread according to the radiation line type as a first signal based on the detected radiation. It has a radiation detection section 10d and an information processing section 5a that specifies the radionuclide included in the measurement target 70 and calculates the radioactivity intensity [Bq] of the radionuclide.
- the information processing unit 5a is, for example, a computer, and is capable of implementing the radioactivity analysis method according to the fifth embodiment.
- the radiation detection unit 10d includes a substance that interacts with the incident radiation, and generates a signal indicating the position where the radiation imparts energy to the substance and the amount of energy imparted to the substance within the interacting substance as a first signal D5. It has a position detection type detector 11d that outputs.
- the radiation type determination unit 20d determines the radiation type based on a first signal D5 including a position where energy is applied to the substance and a signal indicating the amount of energy applied to the substance, and determines the determined radiation type.
- a second signal D2 containing information indicating .
- the position detection type detector 11d has a function of detecting a position to which energy is applied by the incident radiation, and outputs a pulse having an amplitude corresponding to the applied energy.
- the position detection type detector 11d has a function of detecting a position to which energy is applied by the incident radiation, and has a function of outputting a luminance value of light generated at the position to which energy is applied.
- the position detection type detector 11d outputs information on the position where energy is applied and the pulse amplitude.
- the position detection type detector 11d outputs an image including the position where energy is applied and brightness value information.
- the radiation type determination unit 20d determines the radiation type based on the position information and the pulse height value output from the position detection type detector 11d, or the image including the position information and the brightness value.
- the position detection type detector 11d can be realized by the same method as the detector 11c in the fourth embodiment.
- FIG. 11 is a diagram showing the relationship between the brightness value in the position detection type detector 11d of the radiation detection unit 10d and the spread of the trajectory shape.
- the relationship between the brightness value and the spread of the locus shape differs depending on the radiation line type.
- the luminance value of the light emission in the ⁇ -ray trajectory is higher than that of other radiation line types, and the spread of the brightness value (width in the horizontal axis direction in FIG. 11) is larger than that of other radiation line types.
- the spread of the ⁇ -ray locus is higher than that of other radiation line types, and the luminance value of the emission in the ⁇ -ray locus is lower than that of other radiation line types.
- the spread of the gamma-ray trajectory is lower than that of other radiation types, and the luminance value of the luminescence in the gamma-ray trajectory is lower than that in the alpha-ray trajectory.
- the radiation type determination unit 20d can determine the radiation type included in the radiation emitted from the measurement object 70 and generate the second signal D2. .
- the radioactivity analysis device 5 and the radioactivity analysis method according to the fifth embodiment by selecting and using a suitable inverse problem calculation algorithm from a plurality of inverse problem calculation algorithms, It is possible to improve the accuracy of problem calculations.
- the pulse height distribution of the specific radiation type can be extracted, and the wave height distribution can be separated and the inverse problem can be solved. Can solve operations.
- Embodiment 5 is the same as Embodiment 1 or 2.
- 1 to 5 radioactivity analyzer 1a to 5a information processing section, 10, 10b, 10c, 10d radiation detection section, 11, 11b, 11c detector, 11d position detection type detector, 12 amplifier, 13 waveform shaper, 14 Wave height analysis unit, 15 Waveform discriminator, 16 Trajectory shape detection unit, 20, 20b, 20c, 20d Radiation line type determination unit, 30 Inverse problem calculation algorithm selection unit, 31 Inverse problem calculation algorithm DB, 40 Response function selection unit, 41 Response function DB, 50 inverse problem calculation section, 60 display section, 70 measurement object.
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Measurement Of Radiation (AREA)
Abstract
放射能分析装置(1)は、測定対象(70)から放出される放射線を検出し、放射線に基づく第1の信号(D1)を出力する放射線検出部(10)と、第1の信号(D1)に基づいて放射線の種類を示す放射線線種を判定し、放射線線種を示す情報を含む第2の信号(D2)を出力する放射線線種判定部(20)と、記憶装置に予め保存されている複数の逆問題演算アルゴリズムから放射線線種に基づいて逆問題演算アルゴリズムを選択する逆問題演算アルゴリズム選択部(30)と、選択された逆問題演算アルゴリズムと記憶装置に予め保存されている放射線検出部(10)の応答関数とを用いる逆問題演算を実行することによって、第1の信号(D1)と第2の信号(D2)とから測定対象(70)に含まれる放射性核種の特定と放射性核種の放射能の強さの計算とを行う逆問題演算部(50)とを有する。
Description
本開示は、放射能分析装置及び放射能分析方法に関する。
一般に、測定対象の放射能を分析する放射能分析では、放射線検出器によって測定対象から放出される放射線のエネルギー分布が計測され、計測の結果に基づいて測定対象に含まれる放射性物質である放射性核種の特定と放射性核種の放射能の強さ[Bq(ベクレル)]の計算とが行われる。
このように、放射能分析では、結果を観測して原因を推定する逆問題演算が用いられる。例えば、非特許文献1は、放射線検出器から出力されたパルス信号の出力波高分布(すなわち、パルス波高分布)に基づいて、放射線検出器の応答関数を用いて入射放射線のエネルギースペクトルを求める数学的手法であるアンフォールディングを説明している。また、特許文献1は、逆問題演算に用いられる応答関数の数を増やすことで推定結果の精度を向上させる方法を説明している。
原子力規制庁監視情報、放射能測定法シリーズNo.17、「連続モニタによる環境γ線測定法」、原子力規制委員会、平成29年12月改定
しかしながら、上記放射能分析方法では、測定対象に複数の種類の放射線を放出する放射性核種が含まれている場合に、放射能分析の精度が低いという課題があった。
本開示は、上記課題を解決するためになされたものであり、放射能分析の精度が高い放射能分析装置及び放射能分析方法を提供することを目的とする。
本開示の放射能分析装置は、測定対象から放出される放射線を検出し、前記放射線に基づく第1の信号を出力する放射線検出部と、前記第1の信号に基づいて前記放射線の種類を示す放射線線種を判定し、前記放射線線種を示す情報を含む第2の信号を出力する放射線線種判定部と、記憶装置に予め保存されている複数の逆問題演算アルゴリズムから前記放射線線種に基づいて逆問題演算アルゴリズムを選択する逆問題演算アルゴリズム選択部と、前記選択された逆問題演算アルゴリズムと記憶装置に予め保存されている前記放射線検出部の応答関数とを用いる逆問題演算を実行することによって、前記第1の信号と前記第2の信号とから前記測定対象に含まれる放射性核種の特定と前記放射性核種の放射能の強さの計算とを行う逆問題演算部とを有することを特徴とする。
本開示の放射能分析方法は、コンピュータによって実行される方法であって、測定対象から放出される放射線を検出し、前記放射線に基づく第1の信号を出力する放射線検出部から前記第1の信号を受け取るステップと、前記第1の信号に基づいて前記放射線の種類を示す放射線線種を判定し、前記放射線線種を示す情報を含む第2の信号を出力するステップと、記憶装置に予め保存されている複数の逆問題演算アルゴリズムから前記放射線線種に基づいて逆問題演算アルゴリズムを選択するステップと、前記選択された逆問題演算アルゴリズムと記憶装置に予め保存されている前記放射線検出部の応答関数とを用いる逆問題演算を実行することによって、前記第1の信号と前記第2の信号とから前記測定対象に含まれる放射性核種の特定と前記放射性核種の放射能の強さの計算とを行うステップとを有することを特徴とする。
本開示の放射能分析装置及び放射能分析方法によれば、精度の高い放射能分析を行うことができる。
以下に、実施の形態に係る放射能分析装置及び放射能分析方法を、図面を参照しながら説明する。以下の実施の形態は、例にすぎず、実施の形態を適宜組み合わせること及び各実施の形態を適宜変更することが可能である。
《1》実施の形態1
〈放射能分析装置1〉
図1は、実施の形態1に係る放射能分析装置1の構成を概略的に示すブロック図である。放射能分析装置1は、測定対象(「測定資料」又は「被測定物」ともいう。)70から放出される放射線を検出し、検出された放射線に基づく第1の信号D1を出力する放射線検出部10と、測定対象70に含まれる放射性物質(「放射性核種」又は「核種」ともいう。)の特定と放射性核種の放射能の強さ[Bq(ベクレル)]の計算とを行う情報処理部1aとを有している。情報処理部1aは、例えば、情報処理装置としてのコンピュータであり、実施の形態1に係る放射能分析方法を実施することができる。
〈放射能分析装置1〉
図1は、実施の形態1に係る放射能分析装置1の構成を概略的に示すブロック図である。放射能分析装置1は、測定対象(「測定資料」又は「被測定物」ともいう。)70から放出される放射線を検出し、検出された放射線に基づく第1の信号D1を出力する放射線検出部10と、測定対象70に含まれる放射性物質(「放射性核種」又は「核種」ともいう。)の特定と放射性核種の放射能の強さ[Bq(ベクレル)]の計算とを行う情報処理部1aとを有している。情報処理部1aは、例えば、情報処理装置としてのコンピュータであり、実施の形態1に係る放射能分析方法を実施することができる。
放射能分析装置1の情報処理部1aは、放射線線種判定部20と逆問題演算アルゴリズム選択部30と逆問題演算部50とを有している。逆問題演算アルゴリズム選択部30は、記憶装置に予め保存されている複数の逆問題演算アルゴリズム(例えば、逆問題演算アルゴリズムデータベース31として保存されている。)から、第2の信号D2が示す放射線の種類を示す放射線線種に基づいて逆問題演算アルゴリズムを選択する。逆問題演算部50は、選択された逆問題演算アルゴリズムと記憶装置に予め保存されている放射線検出部10の応答関数(例えば、応答関数データベース41として保存されている。)とを用いる逆問題演算を実行することによって、第1の信号D1と第2の信号D2とから測定対象70に含まれる放射性核種の特定と放射性核種の放射能の強さ[Bq]の計算とを行う。測定対象70に複数の放射性核種が含まれている場合には、逆問題演算部50は、複数の放射性核種の各々の特定と複数の放射性核種の各々の放射能の強さ[Bq(ベクレル)]の計算とを行うことができる。また、逆問題演算部50は、逆問題演算を実行することによって、第1の信号D1と第2の信号D2とから放射線の強さ(例えば、放射線のフルエンス(フラックス)[1/(cm2・s)])の計算を行ってもよい。
計算される放射能は、測定対象70の単位重量(又は単位体積)あたりの放射能の強さである放射能濃度[Bq/g](又は[Bq/cm3]など)で表されてもよい。また、逆問題演算アルゴリズムデータベース(「逆問題演算アルゴリズムDB」とも記す。)31が保存される記憶装置と応答関数データベース(「応答関数DB」とも記す。)41が保存される記憶装置とは、同じ装置又は互いに異なる装置のいずれであってもよい。逆問題演算アルゴリズムDB31が保存される記憶装置と応答関数DB41が保存される記憶装置とは、放射能分析装置1の一部であってもよいが、放射能分析装置1と通信可能な他の装置(例えば、ネットワーク上のコンピュータ)の一部であってもよい。
また、放射能分析装置1は、逆問題演算部50の演算結果をユーザに提示する情報出力装置としての表示部60に接続されている。ただし、表示部60は、放射能分析装置1の一部であってもよい。表示部60は、例えば、画像を表示する液晶ディスプレイなどの表示装置である。表示部60の代わりに、又は、表示部60に加えて、他の情報出力装置(例えば、プリンタ、など)を備えることもできる。
〈放射線検出部10〉
放射線検出部10は、測定対象70から放出される放射線を検出する検出器11と、検出器11から出力される出力パルスを増幅する増幅器12と、増幅器12から出力される出力パルスを整形する波形整形器13と、波形整形器13から出力される出力パルスの波高分布を生成する波高分析部14とを有している。
放射線検出部10は、測定対象70から放出される放射線を検出する検出器11と、検出器11から出力される出力パルスを増幅する増幅器12と、増幅器12から出力される出力パルスを整形する波形整形器13と、波形整形器13から出力される出力パルスの波高分布を生成する波高分析部14とを有している。
まず、検出器11について説明する。検出器11は、放射線が検出器に与えたエネルギーに応じた信号強度(例えば、パルス振幅)を出力する機能を持つ。検出器11として広く用いられるものとして、シンチレータ結晶と光電子増倍管とを組み合わせたシンチレーション検出器、シンチレータ結晶と半導体の光検出器とを組み合わせたシンチレーション検出器、半導体検出器、又はガス検出器、などが知られている。半導体検出器は、例えば、ゲルマニウム半導体検出器である。ガス検出器は、容器内のガス中を通過する放射線としての荷電粒子とガス分子の中の電子との衝突によって生成された電子の数(すなわち、荷電粒子のエネルギー損失の量に比例する値)を計測する検出器である。ガス検出器としては、ガイガー・ミュラー(GM)計数管、ガスが充填され電離箱を用いた検出器がある。
次に、増幅器12について説明する。増幅器12は、検出器11から出力された信号である出力パルスの振幅を増幅する機能を有する回路である。増幅器12は、出力パルスを後段の波高分析部14及び放射線線種判定部20での検出に十分な信号強度(例えば、十分に大きなパルス振幅)に増幅するために用いられる。検出器11の出力が十分に大きい場合には、増幅器12を備える必要はない。
次に、波形整形器13について説明する。波形整形器13は、増幅器12又は検出器11から出力される信号のノイズ除去、及び増幅率の調整を行う機能を有する回路である。波形整形器13は、アナログ回路で形成される場合には、一般的に積分回路及び微分回路を組み合わせて構成される。波形整形器13は、ノイズ除去と増幅率調整とを行い、出力パルスを、波高分析部14における波高検出が容易になるような形状の波形にする。波形整形器13は、デジタル回路で形成される場合には、アナログデジタル変換によって前段から出力される信号をデジタル変換し、デジタルフィルタによってノイズ除去と増幅率調整とを行い、波形整形したデジタルでのパルス処理を行う。
次に、波高分析部14について説明する。波高分析部14は、波形整形器13から出力される出力パルスのピーク値の頻度分布又は出力パルスの時間幅の頻度分布を作成する。波形整形器13で前述のようなアナログ回路の処理が行われる場合には、波高分析部14は、出力パルスのピーク値を保持する回路と、アナログデジタル変換を組み合わせて出力パルスの波高を取得する、又は、出力パルスの立上りと立下りの時間を電圧に変換する回路とアナログデジタル変換を組み合わせることで出力パルスの波高を取得する。波形整形器13で前述のようなデジタル回路の処理が行われる場合には、波高分析部14は、デジタルフィルタによって処理された出力パルスの振幅又は時間幅から出力パルスの波高を取得する。波高分析部14は、取得した出力パルスの波高から頻度分布を作成する。
図2は、放射線検出部10から出力されるパルス信号の波高(すなわち、エネルギー)分布を放射線線種ごとに示す図である。図2に示されるように、計測される波高分布は、エネルギー分布に変換されても良い。波高分布とエネルギー分布とは、一対一に対応している。例えば、放射線線種は、α線、β線、γ線のうちの1つ以上である。
〈放射線線種判定部20〉
次に、放射線線種判定部20について説明する。放射線線種判定部20は、波高分析部14で得られた波高分布又はエネルギー分布から測定対象70から放出される放射線の種類(すなわち、放射線線種)を検出する機能を有する。
次に、放射線線種判定部20について説明する。放射線線種判定部20は、波高分析部14で得られた波高分布又はエネルギー分布から測定対象70から放出される放射線の種類(すなわち、放射線線種)を検出する機能を有する。
図2に示されるように、一般的に、放射線検出部10で検出されるアルファ線(α線)は、α線放出核種から放出され、4MeV以上のエネルギーを持つ。例えば、プルトニウム238(238Pu)は、5.50MeVのエネルギーを持つα線と5.46MeVのエネルギーを持つα線とを放出する。プルトニウム239(239Pu)は、5.16MeVのエネルギーを持つα線を放出する。アメリシウム241(241Am)は、5.49MeVのエネルギーを持つα線を放出する。キュリウム244(244Cm)は、5.95MeVのエネルギーを持つα線を放出する。高いエネルギー(一般には、4MeV以上のエネルギー)持つα線は、検出器11によって、図2の実線のように、高いエネルギーの領域に波高分布が形成されるよう検出される。測定対象70から放出されるα線のエネルギーは、空気中でのエネルギー損失により非対称なピーク構造になるが、波高分布上にピーク構造が現れる特徴を示す。放射線線種判定部20は、このピーク構造の有無を検知することで、測定対象70にα線放出核種が含まれるかどうかを判定することができる。
図2に点線で示されるように、一般的に、放射性核種から放出されるガンマ線(γ線)は、核種特有のエネルギーを持ち、γ線の波高分布は、正規分布に近い形状のピーク構造と検出器11の内部又は周辺での散乱によってエネルギーを失う過程を反映した連続成分とを有する。一般的な放射性核種のγ線のエネルギーは、3MeV以下であり、α線の波高分布のピークより低い。例えば、タリウム(Tl)は、2.6MeVのエネルギーを持つγ線を放出する。イットリウム(Y)は、1.8MeVのエネルギーを持つγ線を放出する。コバルト60(60Co)は、1.3MeVと1.1MeVのエネルギーを持つ2本のγ線を放出する。γ線放出核種から放出されるγ線のエネルギーは、α線放出核種から放出されるα線のエネルギーよりも低く、ピーク構造の形状も異なる(急峻なピーク構造と、連続成分とを有する。)ため、放射線線種判定部20は、波高分布から測定対象70にγ線放出核種が含まれるかどうかを判定することができる。
放射性核種から放出されるベータ線(β線)は、原子核の放射性崩壊の一種であるベータ崩壊の過程で放出される。ベータ崩壊では、電子(すなわち、β線)と反電子ニュートリノとが放出される。ベータ崩壊では、電子と反電子ニュートリノとがエネルギーを分け合うために、放出される電子(すなわち、β線)のエネルギーであるβ線の波高分布は、は、図2に一点鎖線で示されるように、連続的な分布となる。β線の波高分布は、検出器11の構造によるが、ピーク構造が低いエネルギーの領域に現れるが、α線及びγ線のように特定のエネルギーの領域にピーク構造を持たない。β線の波高分布は、γ線の波高分布と同様に、低いエネルギーの領域で連続的な分布になる。また、低いエネルギーの領域で、β線の波高分布は、γ線の波高分布より大きい。β線の波高分布の形状は、他の放射線線種と異なるので、放射線線種判定部20は、波高分布から測定対象70にβ線放出核種が含まれるかどうかを判定することができる。
波高分布での放射線線種ごとの特徴の検出は、例えば、上記の特徴を解析的に検出する方法、事前に想定した波高分布との一致度検証によって検出する方法、事前に想定した波高分布の機械学習によって放射線線種を判定する方法、教師無し機械学習によって放射線線種を判定するように学習する方法のいずれかによって実現可能である。事前に想定した波高分布は、様々な条件での波高分布を計算することで作成することができる。波高分布の計算には、例えば、放射線挙動解析による計算がある。
波高分析部14の波高分布を用いて検出を行う機能は、放射線線種判定部20が行ってもよい。その場合には、波高分析部14を備えないことが可能である。
従来の技術では、放射能分析を行うために、特許文献1のように応答関数を調整することで精度を上げることが提案されてきた。実施の形態1の手法では、測定対象70から放出される放射線の種類が1種類である場合が想定されており、複数の放射線線種が測定対象70に含まれる場合には、応答関数の形状が変わることは考慮されていなかった。この結果、応答関数が変わった場合には、単一の逆問題演算アルゴリズムでは、複数の放射線線種が測定される条件で正しく逆問題演算を行うことができなかったため、正しく放射線の強度と放射能の強さを測定することができなかった。
実施の形態1では、検出器11から出力される放射線線種を判定し、検出された放射線線種に応じて使用する逆問題演算アルゴリズムを選ぶことで、測定対象70に複数の放射線線種が含まれる場合においても、それぞれの放射線線種の放射能の強さを求めることが可能となる。
波高分布での放射線線種ごとの特徴の検出は、例えば、上記の特徴を解析的に検出する方法、事前に想定した波高分布との一致度検証によって検出する方法、事前に想定した波高分布の機械学習によって放射線線種を判定する方法、教師無し機械学習によって放射線線種を判定するように学習することで実現可能である。事前に想定した波高分布は、様々な条件での波高分布を計算することで作成することができる。波高分布の計算には、例えば、放射線挙動解析による計算がある。
〈逆問題演算部50〉
次に、逆問題演算部50について説明する。逆問題演算部50は、波高分布から測定対象70に含まれる放射性核種の核種の特定及び放射性核種の放射能の強さ[Bq]の計算を行う。なお、放射能の強さ[Bq]は、放射能濃度で表されてもよい。一般的な測定では、測定対象70に含まれる放射性核種の核種、放射性核種の放射能の強さ[Bq]、放射能濃度が既知の状態で測定を行い、波高分布を得ることを順問題とされる。放射能分析では、測定対象70の情報が不明で測定結果である波高分布のみ知ることができるため順問題を逆に解く逆問題を解く必要がある。放射能分析で逆問題演算を行うためには、一般的に、波高分布と、測定対象70の放射能の強さ又は放射能濃度とを関係づける応答関数が用いられる。波高Eの頻度を測定して得られた波高分布M(E)は、以下の式(1)で求めることができる。
次に、逆問題演算部50について説明する。逆問題演算部50は、波高分布から測定対象70に含まれる放射性核種の核種の特定及び放射性核種の放射能の強さ[Bq]の計算を行う。なお、放射能の強さ[Bq]は、放射能濃度で表されてもよい。一般的な測定では、測定対象70に含まれる放射性核種の核種、放射性核種の放射能の強さ[Bq]、放射能濃度が既知の状態で測定を行い、波高分布を得ることを順問題とされる。放射能分析では、測定対象70の情報が不明で測定結果である波高分布のみ知ることができるため順問題を逆に解く逆問題を解く必要がある。放射能分析で逆問題演算を行うためには、一般的に、波高分布と、測定対象70の放射能の強さ又は放射能濃度とを関係づける応答関数が用いられる。波高Eの頻度を測定して得られた波高分布M(E)は、以下の式(1)で求めることができる。
ここで、R(E)は応答関数であり、S(E)は測定対象70のエネルギーごとの、放射能の強さ(又は放射能濃度)である。この場合に、逆問題演算を行うためのアルゴリズムとしては、例えば、逆行列演算、擬似逆行列演算、デコンボリューション、逆畳み込み積分、及び逐次近似法がある。応答関数は、例えば、放射線挙動解析のように放射線の物理プロセスを反映した演算をすることができるものを用いることで用意することができる。また、式(1)のような解析的な演算ではなく、様々な条件のM(E)とS(E)の組み合わせから、M(E)からS(E)を推定できるように学習させた機械学習を用いることでも演算が実現できる。M(E)とS(E)は、上記と同様に、放射線挙動解析を用いることで作成することができる。
上記のように、式(1)のように解析的な関係ではなく、既知のS(E)とM(E)との関係から機械学習で用いられるネットワーク構造によるS(E)とM(E)の関係から逆問題演算を解くこともできる。
〈逆問題演算アルゴリズム選択部30〉
次に、逆問題演算アルゴリズム選択部30について説明する。上記のように逆問題演算を行うことができる複数の逆問題演算アルゴリズムが予め用意されている。演算は、応答関数R(E)である行列の性質によって解を求めることができない場合もある。応答関数R(E)は、放射線線種によって異なる。
次に、逆問題演算アルゴリズム選択部30について説明する。上記のように逆問題演算を行うことができる複数の逆問題演算アルゴリズムが予め用意されている。演算は、応答関数R(E)である行列の性質によって解を求めることができない場合もある。応答関数R(E)は、放射線線種によって異なる。
以下に、式(2)として、γ線の応答関数R(E)を示す。以下の応答関数R(E)は、対角成分が0より大きい値を持つ行列で表される。
以下に、式(3)として、α線の応答関数R(E)を示す。以下の応答関数R(E)は、対角成分が0である行列で表される。
γ線のように、空気中でのエネルギーの連続的な減衰が無い放射線線種についての応答関数は、式(2)のように、対角成分を持つ(すなわち、対角成分が0ではない)行列が応答関数となる。このような応答関数を持つ放射線線種に対して逆問題演算を行う場合には、逆行列演算、擬似逆行列演算、又は逐次近似法などの解析的な演算アルゴリズムを用いることができる。
α線のように、空気中での減衰がある放射線線種についての応答関数は、式(3)のように、対角成分を持たない行列となる。例えば、対角成分を持たない行列で応答関数が表される場合には、逆行列が存在しないため逆行列演算は用いることができない。また、擬似逆行列演算は、0以下の解が存在するなど演算精度に問題があることがある。
β線の逆問題演算を行う場合には、応答関数を想定される放射性核種からの推定スペクトルで作ることができる。この場合には、応答関数が正方行列では無い条件となり、逆行列が存在しないため、解析的な逆問題演算よりも機械学習によるネットワークによる逆問題演算が最適な結果を推定できる場合がある。逆問題演算アルゴリズム選択部30は、検出した放射線線種に基づいて解の収束又は演算速度の速くなるように1つ又は複数の演算アルゴリズムを検出した放射線線種に応じて最適な選択を行う。
逆問題演算アルゴリズムDB31は、上記のように逆問題演算部50で用いることができる複数の演算アルゴリズムを格納している。
〈応答関数データベース41〉
次に、応答関数データベース41について説明する。応答関数DB41は、逆問題演算に用いる応答行列又は機械学習のネットワーク構造のパラメータを持つ。解析的に逆問題演算を行う場合には、式(2)及び式(3)のような行列を保持する。また、解析的な演算以外の方法である機械学習を用いた逆問題演算を行うためのネットワークに関する情報を持っても良い。
次に、応答関数データベース41について説明する。応答関数DB41は、逆問題演算に用いる応答行列又は機械学習のネットワーク構造のパラメータを持つ。解析的に逆問題演算を行う場合には、式(2)及び式(3)のような行列を保持する。また、解析的な演算以外の方法である機械学習を用いた逆問題演算を行うためのネットワークに関する情報を持っても良い。
〈表示部60〉
逆問題演算で導出された測定対象70に含まれる放射性核種、放射能の強さ[Bq]、放射線の強さ、及び放射能濃度は、表示部60にて表示する。
逆問題演算で導出された測定対象70に含まれる放射性核種、放射能の強さ[Bq]、放射線の強さ、及び放射能濃度は、表示部60にて表示する。
図3は、実施の形態1に係る放射能分析装置1のハードウェア構成の例を示す図である。図3のハードウェア構成は、実施の形態2~5に係る放射能分析装置にも適用される。図3に示されるように、放射能分析装置1は、放射線検出部10と情報処理部1aとを有している。情報処理部1aは、CPU(Central Processing Unit)などのプロセッサ101と、揮発性の記憶装置であるメモリ102と、ハードディスクドライブ(HDD)又はソリッドステートドライブ(SSD)などの不揮発性の記憶装置103とを有している。メモリ102は、例えば、RAM(Random Access Memory)などの半導体メモリである。記憶装置103は、図1に示されるデータベースを記憶するものであってもよい。
情報処理部1aの各機能は、例えば、処理回路により実現される。処理回路は、専用のハードウェアであってもよく、又はメモリ102に格納されるプログラムを実行するプロセッサ101であってもよい。プロセッサ101は、処理装置、演算装置、マイクロプロセッサ、マイクロコンピュータ、及びDSP(Digital Signal Processor)のいずれであってもよい。
処理回路が専用のハードウェアである場合、処理回路は、例えば、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ASIC(Application Specific Integrated Circuit)、FPGA(Field-Programmable Gate Array)、又はこれらのうちのいずれかを組み合わせたものである。
処理回路がプロセッサ101である場合、情報処理部1aによって実行される放射能分析プログラムは、ソフトウェア、ファームウェア、又はソフトウェアとファームウェアとの組み合わせにより実現される。放射能分析プログラムは、ネットワークを経由して又は記録媒体から情報処理部1aにインストールされる。ソフトウェア及びファームウェアは、プログラムとして記述され、メモリ102に格納される。プロセッサ101は、メモリ102に記憶された表示制御プログラムを読み出して実行することにより、図1に示される各部の機能を実現することができる。
なお、情報処理部1aは、一部を専用のハードウェアで実現し、一部をソフトウェア又はファームウェアで実現するようにしてもよい。このように、処理回路は、ハードウェア、ソフトウェア、ファームウェア、又はこれらのうちのいずれかの組み合わせによって、上述の各機能を実現することができる。
図4は、実施の形態1に係る放射能分析方法を示すフローチャートである。この放射能分析方法は、情報処理部1aとしてのコンピュータによって実行される。この放射能分析方法は、測定対象70から放出される放射線を検出し、放射線に基づく第1の信号D1を出力する放射線検出部10から第1の信号D1を受け取るステップS1と、第1の信号に基づいて放射線の種類を示す放射線線種を判定し、放射線線種を示す情報を含む第2の信号D2を出力するステップS2と、記憶装置に予め保存されている複数の逆問題演算アルゴリズムから放射線線種に基づいて逆問題演算アルゴリズムを選択するステップS3と、選択された逆問題演算アルゴリズムと記憶装置に予め保存されている前記放射線検出部10の応答関数とを用いる逆問題演算を実行することによって、第1の信号D1と前記第2の信号D2とから測定対象70に含まれる放射性核種の特定と前記放射性核種の放射能の強さの計算とを行うステップS4とを有している。また、ステップS4において、放射線の強さ(例えば、放射線のフルエンス)の計算を行ってもよい。
〈効果〉
以上に説明したように、実施の形態1に係る放射能分析装置1及び放射能分析方法によれば、複数の逆問題演算アルゴリズムから、適した逆問題演算アルゴリズムを選択して使用することで、逆問題演算の精度を上げることができる。
以上に説明したように、実施の形態1に係る放射能分析装置1及び放射能分析方法によれば、複数の逆問題演算アルゴリズムから、適した逆問題演算アルゴリズムを選択して使用することで、逆問題演算の精度を上げることができる。
《2》実施の形態2
図5は、実施の形態2に係る放射能分析装置2の構成を概略的に示すブロック図である。図5において、図1に示される構成と同一又は対応する構成には、図1に示される符号と同じ符号が付されている。実施の形態2に係る放射能分析装置2は、情報処理部2aが応答関数選択部40を有する点において、実施の形態1に係る放射能分析装置1と異なる。
図5は、実施の形態2に係る放射能分析装置2の構成を概略的に示すブロック図である。図5において、図1に示される構成と同一又は対応する構成には、図1に示される符号と同じ符号が付されている。実施の形態2に係る放射能分析装置2は、情報処理部2aが応答関数選択部40を有する点において、実施の形態1に係る放射能分析装置1と異なる。
放射能分析装置2は、測定対象70から放出される放射線を検出し、検出された放射線に基づく第1の信号D1を出力する放射線検出部10と、測定対象70に含まれる放射性物質を示す放射性核種の特定と放射性核種の放射能の強さ[Bq]の計算とを行う情報処理部2aとを有している。情報処理部2aは、例えば、コンピュータであり、実施の形態2に係る放射能分析方法を実施することができる。
実施の形態2においては、記憶装置(例えば、図3の記憶装置103)に応答関数DB41として予め複数の応答関数が応答関数DB41として保存されている。応答関数選択部40は、放射線線種判定部20において判定された放射線線種を示す第2の信号D2に基づいて、応答関数DB41に保存されている複数の応答関数から、使用する応答関数を選択する。逆問題演算部50は、逆問題演算アルゴリズム選択部30によって選択された逆問題演算アルゴリズムと応答関数選択部40によって選択された応答関数とを用いて、逆問題演算を実行することによって、第1の信号D1と第2の信号D2とから測定対象70に含まれる放射性核種の特定と放射性核種の放射能の強さ[Bq]の計算とを行う。また、逆問題演算部50は、逆問題演算を実行することによって、第1の信号D1と第2の信号D2とから放射線の強さの計算を行ってもよい。
以上に説明したように、実施の形態2に係る放射能分析装置2及び放射能分析方法によれば、複数の逆問題演算アルゴリズムから、適した逆問題演算アルゴリズムを選択し、複数の応答関数から適した応答関数を選択し、これらを使用することで、逆問題演算の精度を上げることができる。
上記以外に関し、実施の形態2は、実施の形態1と同じである。
《3》実施の形態3
図6は、実施の形態3に係る放射能分析装置3の構成を概略的に示すブロック図である。図6において、図1に示される構成と同一又は対応する構成には、図1に示される符号と同じ符号が付されている。実施の形態3に係る放射能分析装置3は、放射線検出部10bが、検出器11によって検出された放射線の放射線種別ごとに検出信号を弁別して弁別信号を出力する波形弁別器15とを有する点において、実施の形態1に係る放射能分析装置1と異なる。波形弁別器15は、検出器11の出力及び増幅器12から出力する出力信号のパルス形状が放射線線種によって差が出るようにする回路である。
図6は、実施の形態3に係る放射能分析装置3の構成を概略的に示すブロック図である。図6において、図1に示される構成と同一又は対応する構成には、図1に示される符号と同じ符号が付されている。実施の形態3に係る放射能分析装置3は、放射線検出部10bが、検出器11によって検出された放射線の放射線種別ごとに検出信号を弁別して弁別信号を出力する波形弁別器15とを有する点において、実施の形態1に係る放射能分析装置1と異なる。波形弁別器15は、検出器11の出力及び増幅器12から出力する出力信号のパルス形状が放射線線種によって差が出るようにする回路である。
放射線線種判定部20bは、弁別信号を含む第1の信号D1に基づいて放射線線種を判定し、判別された放射線線種を示す情報を含む第2の信号D2を出力する。放射能分析装置3は、測定対象70から放出される放射線を検出し、放射線の波高分布を示す弁別信号を含む第1の信号D1を出力する放射線検出部10bと、測定対象70に含まれる放射性核種の特定と放射性核種の放射能の強さ[Bq]の計算とを行う情報処理部3aとを有している。情報処理部3aは、例えば、コンピュータであり、実施の形態3に係る放射能分析方法を実施することができる。
図7は、放射線検出部10bの検出器11から出力されるパルス信号の振幅の立ち上がり及び減衰特性を放射線線種ごとに示す図である。例えば、検出器11の、放射線と相互作用する物質として半導体を用いた場合、半導体の内部では、放射線の通過によって電子と正孔とが発生するが、電子の移動度と正孔の移動度とは異なる。α線のように検出器11の表面付近で反応する放射線線種では、正孔の移動距離を長くするように検出器11を設定することでパルスの立上りが緩やかな形状になる(図7に実線で示される)。また、γ線及びβ線のように検出器11の内部(表面よりも深い位置)で反応する放射線線種では、立上りが急峻な形状となる(図7に破線で示される)。波形弁別器15は、このような波形形状の違いを用いて、放射線線種を弁別することで特定の放射線線種を精度良く検出できる。
実施の形態3においては、逆問題演算部50は、逆問題演算アルゴリズム選択部30によって選択された逆問題演算アルゴリズムと、応答関数とを用いて、逆問題演算を実行することによって、第1の信号D1と第2の信号D2とから測定対象70に含まれる放射性核種の特定と放射性核種の放射能の強さ[Bq]の計算とを行う。また、逆問題演算部50は、逆問題演算を実行することによって、第1の信号D1と第2の信号D2とから放射線の強さの計算を行ってもよい。
以上に説明したように、実施の形態3に係る放射能分析装置3及び放射能分析方法によれば、複数の逆問題演算アルゴリズムから、適した逆問題演算アルゴリズムを選択し使用することで、逆問題演算の精度を上げることができる。
また、波形弁別器15によって特定の放射線線種を判定し、判定の結果を放射線線種判定部20での波高分布からの検出を組み合わせることで、検出精度を上げることができる。また、特定の放射線線種の検出結果と波高分布とを対応付けることで特定の放射線線種の波高分布を取り出すことで波高分布を分離して逆問題演算を解くことができる。
上記以外に関し、実施の形態3は、実施の形態1又は2と同じである。
《4》実施の形態4
図8は、実施の形態4に係る放射能分析装置4の構成を概略的に示すブロック図である。図8において、図1に示される構成と同一又は対応する構成には、図1に示される符号と同じ符号が付されている。実施の形態4に係る放射能分析装置4は、放射線検出部10cの検出器11cが放射線によってエネルギー付与された位置を検出する機能を有する点及び放射線検出部10cが軌跡形状検出部16を有する点において、実施の形態1に係る放射能分析装置1と異なる。
図8は、実施の形態4に係る放射能分析装置4の構成を概略的に示すブロック図である。図8において、図1に示される構成と同一又は対応する構成には、図1に示される符号と同じ符号が付されている。実施の形態4に係る放射能分析装置4は、放射線検出部10cの検出器11cが放射線によってエネルギー付与された位置を検出する機能を有する点及び放射線検出部10cが軌跡形状検出部16を有する点において、実施の形態1に係る放射能分析装置1と異なる。
放射能分析装置4は、測定対象70から放出される放射線を検出し、検出された放射線に基づく第1の信号D1と放射線の軌跡形状を示す軌跡信号を出力する放射線検出部10cと、測定対象70に含まれる放射性核種の特定と放射性核種の放射能の強さ[Bq]の計算とを行う情報処理部4aとを有している。情報処理部4aは、例えば、コンピュータであり、実施の形態4に係る放射能分析方法を実施することができる。
放射線検出部10cは、放射線と相互作用する物質を有する検出器11cと、相互作用する物質内において放射線がエネルギーを付与した軌跡の位置を検出し、軌跡の形状を示す軌跡形状信号を出力する回路である軌跡形状検出部16とを有している。放射線線種判定部20cは、軌跡形状信号を含む第1の信号D1に基づいて放射線線種を判定し、判定された放射線線種を示す情報を含む第2の信号D2を出力する。
実施の形態4では、検出器11cが放射線によってエネルギー付与された位置を検出する機能を有する。例えば、検出器11cは、アレイ状に複数のシンチレータ結晶を並べ、アレイ状の並べられた複数のシンチレータ結晶の各々の上に光検出器を配置した装置である(例えば、特許文献2参照)。また、検出器11cは、規則的に配列された複数のシンチレータ結晶と、シンチレータ結晶で発生した光を撮影するカメラ(例えば、CMOS(Complementary Metal Oxide Semiconductor)カメラ素子、SOI(Silicon-On-Insulator)撮像素子、など)とからなる装置であってもよい。また、検出器11cは、半導体放射線検出器(例えば、Ge結晶体)と、これを挟むように、且つ、間隔を開けて規則的に配列された複数のストリップ電極とを持つ装置であってもよい。また、検出器11cは、複数のワイヤーを持つガス式放射線検出器又は規則的に配列された複数のガス式放射線検出器によって放射線によってエネルギー付与された位置を検出する装置であってもよい。
図9(A)から(C)は、放射線検出部10cの検出器11cにおける軌跡形状を放射線線種ごとに示す概略図である。図9(A)から(C)は、放射線線種による検出器11c内における軌跡形状(横から見た場合と上から見た場合)を示す。図9(A)に示されるように、α線のように飛程が短く付与エネルギーが大きい放射線は、広がった点状の形状が軌跡形状となる。図9(B)に示されるように、β線は、飛程が長い特性により、線状に検出器を通過するため細長い形状になる。図9(C)に示されるように、γ線は、一点で反応するため広がりのない点状の形状となる。この形状の違いを放射線線種判定部20cで検出することで、測定対象70から放出される放射線に含まれる放射線線種を判定し、第2の信号D2を生成できる。
以上に説明したように、実施の形態4に係る放射能分析装置4及び放射能分析方法によれば、複数の逆問題演算アルゴリズムから、適した逆問題演算アルゴリズムを選択し使用することで、逆問題演算の精度を上げることができる。
また、軌跡形状検出部16によって特定の放射線線種を判定し、判定の結果を放射線線種判定部20での波高分布からの検出を組み合わせることで、検出精度を上げることができる。また、特定の放射線線種の検出結果と波高分布とを対応付け特定の放射線線種の波高分布を取り出すことで波高分布を分離して逆問題演算を解くことができる。
上記以外に関し、実施の形態4は、実施の形態1又は2と同じである。
《5》実施の形態5
図10は、実施の形態5に係る放射能分析装置5の構成を概略的に示すブロック図である。図10において、図1に示される構成と同一又は対応する構成には、図1に示される符号と同じ符号が付されている。実施の形態5に係る放射能分析装置5は、放射線検出部10dが位置検出型検出器11dで構成されている点において、実施の形態1に係る放射能分析装置1と異なる。
図10は、実施の形態5に係る放射能分析装置5の構成を概略的に示すブロック図である。図10において、図1に示される構成と同一又は対応する構成には、図1に示される符号と同じ符号が付されている。実施の形態5に係る放射能分析装置5は、放射線検出部10dが位置検出型検出器11dで構成されている点において、実施の形態1に係る放射能分析装置1と異なる。
放射能分析装置5は、測定対象70から放出される放射線を検出し、検出された放射線に基づく第1の信号として、放射線線種に応じた輝度値と軌跡の拡がりと示す信号D5を出力する放射線検出部10dと、測定対象70に含まれる放射性核種の特定と放射性核種の放射能の強さ[Bq]の計算とを行う情報処理部5aとを有している。情報処理部5aは、例えば、コンピュータであり、実施の形態5に係る放射能分析方法を実施することができる。
放射線検出部10dは、入射する放射線と相互作用する物質を備え、相互作用する物質内において放射線が物質にエネルギーを付与した位置と物質に付与したエネルギーの量を示す信号を第1の信号D5として出力する位置検出型検出器11dを有している。放射線線種判定部20dは、物質にエネルギーを付与した位置と物質に付与したエネルギーの量を示す信号とを含む第1の信号D5に基づいて放射線線種を判定し、判定された放射線線種を示す情報を含む第2の信号D2を出力する。
位置検出型検出器11dは、入射する放射線によってエネルギー付与された位置を検出する機能を有し、付与されたエネルギーに対応する振幅を持つパルスを出力する。或いは、位置検出型検出器11dは、入射する放射線によってエネルギー付与された位置を検出する機能を有し、エネルギー付与された位置で発生する光の輝度値を出力する機能を有する。位置検出型検出器11dは、エネルギー付与された位置とパルス振幅の情報を出力する。或いは、位置検出型検出器11dは、エネルギー付与された位置と輝度値情報とを含む画像を出力する。放射線線種判定部20dは、位置検出型検出器11dから出力された位置情報とパルス波高値、又は位置情報と輝度値とを含む画像に基づいて、放射線線種を判定する。位置検出型検出器11dは、実施の形態4における検出器11cと同様の方法で実現可能である。
図11は、放射線検出部10dの位置検出型検出器11dにおける輝度値と軌跡形状の広がりとの関係を示す図である。図11に示されるように、放射線線種によって、輝度値と軌跡形状の広がりとの関係が異なる。α線軌跡における発光の輝度値は、他の放射線線種のものよりも高く、輝度値の広がり(図11の横軸方向の幅)は、他の放射線線種のものよりも大きい。β線軌跡の広がりは、他の放射線線種のものよりも高く、β線軌跡における発光の輝度値は、他の放射線線種のものよりも低い。γ線軌跡の広がりは、他の放射線線種のものよりも低く、γ線軌跡における発光の輝度値は、α線軌跡におけるものよりも低い。これらの違いを放射線線種判定部20dで検出することで、放射線線種判定部20dは、測定対象70から放出される放射線に含まれる放射線線種を判定し、第2の信号D2を生成できる。
以上に説明したように、実施の形態5に係る放射能分析装置5及び放射能分析方法によれば、複数の逆問題演算アルゴリズムから、適した逆問題演算アルゴリズムを選択し使用することで、逆問題演算の精度を上げることができる。
また、放射線線種の判定精度を上げることができ、特定の放射線線種の検出結果と波高分布とを対応付けることで特定の放射線線種の波高分布を取り出すことで波高分布を分離して逆問題演算を解くことができる。
上記以外に関し、実施の形態5は、実施の形態1又は2と同じである。
《6》変形例
実施の形態2で説明した、応答関数選択部40を実施の形態3、4、5のいずれかに組み合わせることが可能である。
実施の形態2で説明した、応答関数選択部40を実施の形態3、4、5のいずれかに組み合わせることが可能である。
1~5 放射能分析装置、 1a~5a 情報処理部、 10、10b、10c、10d 放射線検出部、 11、11b、11c 検出器、 11d 位置検出型検出器、 12 増幅器、 13 波形整形器、 14 波高分析部、 15 波形弁別器、 16 軌跡形状検出部、 20、20b、20c、20d 放射線線種判定部、 30 逆問題演算アルゴリズム選択部、 31 逆問題演算アルゴリズムDB、 40 応答関数選択部、 41 応答関数DB、 50 逆問題演算部、 60 表示部、 70 測定対象。
Claims (8)
- 測定対象から放出される放射線を検出し、前記放射線に基づく第1の信号を出力する放射線検出部と、
前記第1の信号に基づいて前記放射線の種類を示す放射線線種を判定し、前記放射線線種を示す情報を含む第2の信号を出力する放射線線種判定部と、
記憶装置に予め保存されている複数の逆問題演算アルゴリズムから前記放射線線種に基づいて逆問題演算アルゴリズムを選択する逆問題演算アルゴリズム選択部と、
前記選択された逆問題演算アルゴリズムと記憶装置に予め保存されている前記放射線検出部の応答関数とを用いる逆問題演算を実行することによって、前記第1の信号と前記第2の信号とから前記測定対象に含まれる放射性核種の特定と前記放射性核種の放射能の強さの計算とを行う逆問題演算部と
を有することを特徴とする放射能分析装置。 - 応答関数選択部をさらに有し、
前記予め保存されている前記応答関数として、複数の応答関数が保存されており、
前記応答関数選択部は、前記複数の応答関数から前記第2の信号が示す放射線線種に基づいて応答関数を選択し、
前記逆問題演算部は、前記選択された逆問題演算アルゴリズムと前記選択された応答関数とを用いて前記逆問題演算を実行することによって、前記第1の信号と前記第2の信号とから前記測定対象に含まれる放射性核種の特定と前記放射性核種の放射能の強さの計算とを行う
ことを特徴とする請求項1に記載の放射能分析装置。 - 前記逆問題演算部は、前記逆問題演算を実行することによって、前記第1の信号と前記第2の信号とから前記放射線の強さの計算を行う
ことを特徴とする請求項1又は2に記載の放射能分析装置。 - 前記放射線検出部は、前記放射線を検出して検出信号を出力する検出器と、前記放射線の放射線種別ごとに前記検出信号を弁別して弁別信号を出力する波形弁別器とを有し、
前記放射線線種判定部は、前記弁別信号を含む前記第1の信号に基づいて前記放射線線種を判定し、前記放射線線種を示す情報を含む前記第2の信号を出力する
ことを特徴とする請求項1から3のいずれか1項に記載の放射能分析装置。 - 前記放射線検出部は、前記放射線と相互作用する物質を有する検出器と、前記物質内において前記放射線がエネルギーを付与した軌跡の位置を検出し、前記軌跡の形状を示す軌跡形状信号を出力する軌跡形状検出部とを有し、
前記放射線線種判定部は、前記軌跡形状信号を含む前記第1の信号に基づいて前記放射線線種を判定し、前記放射線線種を示す情報を含む前記第2の信号を出力する
ことを特徴とする請求項1から3のいずれか1項に記載の放射能分析装置。 - 前記放射線検出部は、前記放射線と相互作用する物質を備え、前記物質内において前記放射線が前記物質にエネルギーを付与した位置と前記物質に付与した前記エネルギーの量を示す信号を前記第1の信号として出力する位置検出型検出器を有し、
前記放射線線種判定部は、前記物質にエネルギーを付与した位置と前記物質に付与した前記エネルギーの量とを示す信号を含む前記第1の信号に基づいて前記放射線線種を判定し、前記放射線線種を示す情報を含む前記第2の信号を出力する
ことを特徴とする請求項1から3のいずれか1項に記載の放射能分析装置。 - 前記放射線線種は、α線、β線、γ線のうちの1つ以上である
ことを特徴とする請求項1から6のいずれか1項に記載の放射能分析装置。 - コンピュータによって実行される放射能分析方法であって、
測定対象から放出される放射線を検出し、前記放射線に基づく第1の信号を出力する放射線検出部から前記第1の信号を受け取るステップと、
前記第1の信号に基づいて前記放射線の種類を示す放射線線種を判定し、前記放射線線種を示す情報を含む第2の信号を出力するステップと、
記憶装置に予め保存されている複数の逆問題演算アルゴリズムから前記放射線線種に基づいて逆問題演算アルゴリズムを選択するステップと、
前記選択された逆問題演算アルゴリズムと記憶装置に予め保存されている前記放射線検出部の応答関数とを用いる逆問題演算を実行することによって、前記第1の信号と前記第2の信号とから前記測定対象に含まれる放射性核種の特定と前記放射性核種の放射能の強さの計算とを行うステップと
を有することを特徴とする放射能分析方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2022/034696 WO2024057514A1 (ja) | 2022-09-16 | 2022-09-16 | 放射能分析装置及び放射能分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2022/034696 WO2024057514A1 (ja) | 2022-09-16 | 2022-09-16 | 放射能分析装置及び放射能分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2024057514A1 true WO2024057514A1 (ja) | 2024-03-21 |
Family
ID=90274692
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2022/034696 WO2024057514A1 (ja) | 2022-09-16 | 2022-09-16 | 放射能分析装置及び放射能分析方法 |
Country Status (1)
Country | Link |
---|---|
WO (1) | WO2024057514A1 (ja) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014157132A (ja) * | 2013-02-18 | 2014-08-28 | Mitsubishi Electric Corp | 放射能分析装置および放射能分析方法 |
JP2015537215A (ja) * | 2012-11-23 | 2015-12-24 | クロメック リミテッドKromek Limited | スペクトルデータの検出・操作方法 |
-
2022
- 2022-09-16 WO PCT/JP2022/034696 patent/WO2024057514A1/ja unknown
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015537215A (ja) * | 2012-11-23 | 2015-12-24 | クロメック リミテッドKromek Limited | スペクトルデータの検出・操作方法 |
JP2014157132A (ja) * | 2013-02-18 | 2014-08-28 | Mitsubishi Electric Corp | 放射能分析装置および放射能分析方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4160275B2 (ja) | エネルギー測定方法及び測定装置 | |
US8084748B2 (en) | Radioactive material detecting and identifying device and method | |
JP4766407B2 (ja) | 放射線線量計および放射線線量計算プログラム | |
US20070290136A1 (en) | Pulse shape discrimination method and apparatus for high-sensitivity radioisotope identification with an integrated neutron-gamma radiation detector | |
EP2883085B1 (en) | Gamma-ray spectrometer | |
WO2013066882A1 (en) | Neutron spectrometer | |
JP4061367B2 (ja) | ZnS(Ag)シンチレーション検出器 | |
JP2014066518A (ja) | 放射能分析装置 | |
Parks | The Compton effect-Compton scattering and gamma ray spectroscopy | |
JPH01260388A (ja) | エネルギー補正オフセットにより像を一様にしたガンマカメラ | |
JP6524484B2 (ja) | 放射線計測方法及び放射線計測装置 | |
JPH05341047A (ja) | 効果的なα及びβ(γ)線同時測定法及びその検出器 | |
Richards et al. | High count rate pulse shape discrimination algorithms for neutron scattering facilities | |
JP7085964B2 (ja) | 放射線検出器 | |
WO2024057514A1 (ja) | 放射能分析装置及び放射能分析方法 | |
JPH11194170A (ja) | 放射性物質検査装置及び放射性廃棄物検査システム | |
JP2002006055A (ja) | 放射線測定方法および装置 | |
KR102663201B1 (ko) | 방사선 측정 장치 및 장치의 동작 방법 | |
KR102280128B1 (ko) | 방사성 핵종 판별 장치 및 방법 | |
JPH09211133A (ja) | 放射線モニタ | |
JP3595375B2 (ja) | ポジトロン検出装置 | |
JP2972936B2 (ja) | シンチレーションパルス波高データの組合わせ方法および装置 | |
RU2390800C2 (ru) | Способ и устройство для измерения спектральной и интегральной плотности потока нейтронов | |
JP7378377B2 (ja) | 放射線分析装置およびダストモニタ装置 | |
JP7352963B2 (ja) | 放射線検出器 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 22958834 Country of ref document: EP Kind code of ref document: A1 |