WO2024057514A1 - 放射能分析装置及び放射能分析方法 - Google Patents

放射能分析装置及び放射能分析方法 Download PDF

Info

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
Application number
PCT/JP2022/034696
Other languages
English (en)
French (fr)
Inventor
理 笹野
真照 林
Original Assignee
三菱電機株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 三菱電機株式会社 filed Critical 三菱電機株式会社
Priority to PCT/JP2022/034696 priority Critical patent/WO2024057514A1/ja
Publication of WO2024057514A1 publication Critical patent/WO2024057514A1/ja

Links

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

放射能分析装置（１）は、測定対象（７０）から放出される放射線を検出し、放射線に基づく第１の信号（Ｄ１）を出力する放射線検出部（１０）と、第１の信号（Ｄ１）に基づいて放射線の種類を示す放射線線種を判定し、放射線線種を示す情報を含む第２の信号（Ｄ２）を出力する放射線線種判定部（２０）と、記憶装置に予め保存されている複数の逆問題演算アルゴリズムから放射線線種に基づいて逆問題演算アルゴリズムを選択する逆問題演算アルゴリズム選択部（３０）と、選択された逆問題演算アルゴリズムと記憶装置に予め保存されている放射線検出部（１０）の応答関数とを用いる逆問題演算を実行することによって、第１の信号（Ｄ１）と第２の信号（Ｄ２）とから測定対象（７０）に含まれる放射性核種の特定と放射性核種の放射能の強さの計算とを行う逆問題演算部（５０）とを有する。

Description

放射能分析装置及び放射能分析方法
　本開示は、放射能分析装置及び放射能分析方法に関する。
　一般に、測定対象の放射能を分析する放射能分析では、放射線検出器によって測定対象から放出される放射線のエネルギー分布が計測され、計測の結果に基づいて測定対象に含まれる放射性物質である放射性核種の特定と放射性核種の放射能の強さ［Ｂｑ（ベクレル）］の計算とが行われる。
　このように、放射能分析では、結果を観測して原因を推定する逆問題演算が用いられる。例えば、非特許文献１は、放射線検出器から出力されたパルス信号の出力波高分布（すなわち、パルス波高分布）に基づいて、放射線検出器の応答関数を用いて入射放射線のエネルギースペクトルを求める数学的手法であるアンフォールディングを説明している。また、特許文献１は、逆問題演算に用いられる応答関数の数を増やすことで推定結果の精度を向上させる方法を説明している。
原子力規制庁監視情報、放射能測定法シリーズＮｏ．１７、「連続モニタによる環境γ線測定法」、原子力規制委員会、平成２９年１２月改定
特許第５８３２４０４号公報
　しかしながら、上記放射能分析方法では、測定対象に複数の種類の放射線を放出する放射性核種が含まれている場合に、放射能分析の精度が低いという課題があった。
　本開示は、上記課題を解決するためになされたものであり、放射能分析の精度が高い放射能分析装置及び放射能分析方法を提供することを目的とする。
　本開示の放射能分析装置は、測定対象から放出される放射線を検出し、前記放射線に基づく第１の信号を出力する放射線検出部と、前記第１の信号に基づいて前記放射線の種類を示す放射線線種を判定し、前記放射線線種を示す情報を含む第２の信号を出力する放射線線種判定部と、記憶装置に予め保存されている複数の逆問題演算アルゴリズムから前記放射線線種に基づいて逆問題演算アルゴリズムを選択する逆問題演算アルゴリズム選択部と、前記選択された逆問題演算アルゴリズムと記憶装置に予め保存されている前記放射線検出部の応答関数とを用いる逆問題演算を実行することによって、前記第１の信号と前記第２の信号とから前記測定対象に含まれる放射性核種の特定と前記放射性核種の放射能の強さの計算とを行う逆問題演算部とを有することを特徴とする。
　本開示の放射能分析方法は、コンピュータによって実行される方法であって、測定対象から放出される放射線を検出し、前記放射線に基づく第１の信号を出力する放射線検出部から前記第１の信号を受け取るステップと、前記第１の信号に基づいて前記放射線の種類を示す放射線線種を判定し、前記放射線線種を示す情報を含む第２の信号を出力するステップと、記憶装置に予め保存されている複数の逆問題演算アルゴリズムから前記放射線線種に基づいて逆問題演算アルゴリズムを選択するステップと、前記選択された逆問題演算アルゴリズムと記憶装置に予め保存されている前記放射線検出部の応答関数とを用いる逆問題演算を実行することによって、前記第１の信号と前記第２の信号とから前記測定対象に含まれる放射性核種の特定と前記放射性核種の放射能の強さの計算とを行うステップとを有することを特徴とする。
　本開示の放射能分析装置及び放射能分析方法によれば、精度の高い放射能分析を行うことができる。
実施の形態１に係る放射能分析装置の構成を概略的に示すブロック図である。 放射線検出部の検出器から出力されるパルス信号の波高（すなわち、エネルギー）分布を放射線の種類を示す放射線線種ごとに示す図である。 実施の形態１に係る放射能分析装置のハードウェア構成を示す図である。 実施の形態１に係る放射能分析方法を示すフローチャートである。 実施の形態２に係る放射能分析装置の構成を概略的に示すブロック図である。 実施の形態３に係る放射能分析装置の構成を概略的に示すブロック図である。 放射線検出部の検出器から出力されるパルス信号の振幅の立ち上がり及び減衰特性を放射線線種ごとに示す図である。 実施の形態４に係る放射能分析装置の構成を概略的に示すブロック図である。 （Ａ）から（Ｃ）は、放射線検出部の検出器における軌跡形状を放射線線種ごとに示す概略図である。 実施の形態５に係る放射能分析装置の構成を概略的に示すブロック図である。 放射線検出部の位置検出型検出器における輝度値と軌跡形状の広がりとの関係を示す図である。
　以下に、実施の形態に係る放射能分析装置及び放射能分析方法を、図面を参照しながら説明する。以下の実施の形態は、例にすぎず、実施の形態を適宜組み合わせること及び各実施の形態を適宜変更することが可能である。
《１》実施の形態１
〈放射能分析装置１〉
　図１は、実施の形態１に係る放射能分析装置１の構成を概略的に示すブロック図である。放射能分析装置１は、測定対象（「測定資料」又は「被測定物」ともいう。）７０から放出される放射線を検出し、検出された放射線に基づく第１の信号Ｄ１を出力する放射線検出部１０と、測定対象７０に含まれる放射性物質（「放射性核種」又は「核種」ともいう。）の特定と放射性核種の放射能の強さ［Ｂｑ（ベクレル）］の計算とを行う情報処理部１ａとを有している。情報処理部１ａは、例えば、情報処理装置としてのコンピュータであり、実施の形態１に係る放射能分析方法を実施することができる。
　放射能分析装置１の情報処理部１ａは、放射線線種判定部２０と逆問題演算アルゴリズム選択部３０と逆問題演算部５０とを有している。逆問題演算アルゴリズム選択部３０は、記憶装置に予め保存されている複数の逆問題演算アルゴリズム（例えば、逆問題演算アルゴリズムデータベース３１として保存されている。）から、第２の信号Ｄ２が示す放射線の種類を示す放射線線種に基づいて逆問題演算アルゴリズムを選択する。逆問題演算部５０は、選択された逆問題演算アルゴリズムと記憶装置に予め保存されている放射線検出部１０の応答関数（例えば、応答関数データベース４１として保存されている。）とを用いる逆問題演算を実行することによって、第１の信号Ｄ１と第２の信号Ｄ２とから測定対象７０に含まれる放射性核種の特定と放射性核種の放射能の強さ［Ｂｑ］の計算とを行う。測定対象７０に複数の放射性核種が含まれている場合には、逆問題演算部５０は、複数の放射性核種の各々の特定と複数の放射性核種の各々の放射能の強さ［Ｂｑ（ベクレル）］の計算とを行うことができる。また、逆問題演算部５０は、逆問題演算を実行することによって、第１の信号Ｄ１と第２の信号Ｄ２とから放射線の強さ（例えば、放射線のフルエンス（フラックス）［１／（ｃｍ２・ｓ）］）の計算を行ってもよい。
　計算される放射能は、測定対象７０の単位重量（又は単位体積）あたりの放射能の強さである放射能濃度［Ｂｑ／ｇ］（又は［Ｂｑ／ｃｍ３］など）で表されてもよい。また、逆問題演算アルゴリズムデータベース（「逆問題演算アルゴリズムＤＢ」とも記す。）３１が保存される記憶装置と応答関数データベース（「応答関数ＤＢ」とも記す。）４１が保存される記憶装置とは、同じ装置又は互いに異なる装置のいずれであってもよい。逆問題演算アルゴリズムＤＢ３１が保存される記憶装置と応答関数ＤＢ４１が保存される記憶装置とは、放射能分析装置１の一部であってもよいが、放射能分析装置１と通信可能な他の装置（例えば、ネットワーク上のコンピュータ）の一部であってもよい。
　また、放射能分析装置１は、逆問題演算部５０の演算結果をユーザに提示する情報出力装置としての表示部６０に接続されている。ただし、表示部６０は、放射能分析装置１の一部であってもよい。表示部６０は、例えば、画像を表示する液晶ディスプレイなどの表示装置である。表示部６０の代わりに、又は、表示部６０に加えて、他の情報出力装置（例えば、プリンタ、など）を備えることもできる。
〈放射線検出部１０〉
　放射線検出部１０は、測定対象７０から放出される放射線を検出する検出器１１と、検出器１１から出力される出力パルスを増幅する増幅器１２と、増幅器１２から出力される出力パルスを整形する波形整形器１３と、波形整形器１３から出力される出力パルスの波高分布を生成する波高分析部１４とを有している。
　まず、検出器１１について説明する。検出器１１は、放射線が検出器に与えたエネルギーに応じた信号強度（例えば、パルス振幅）を出力する機能を持つ。検出器１１として広く用いられるものとして、シンチレータ結晶と光電子増倍管とを組み合わせたシンチレーション検出器、シンチレータ結晶と半導体の光検出器とを組み合わせたシンチレーション検出器、半導体検出器、又はガス検出器、などが知られている。半導体検出器は、例えば、ゲルマニウム半導体検出器である。ガス検出器は、容器内のガス中を通過する放射線としての荷電粒子とガス分子の中の電子との衝突によって生成された電子の数（すなわち、荷電粒子のエネルギー損失の量に比例する値）を計測する検出器である。ガス検出器としては、ガイガー・ミュラー（ＧＭ）計数管、ガスが充填され電離箱を用いた検出器がある。
　次に、増幅器１２について説明する。増幅器１２は、検出器１１から出力された信号である出力パルスの振幅を増幅する機能を有する回路である。増幅器１２は、出力パルスを後段の波高分析部１４及び放射線線種判定部２０での検出に十分な信号強度（例えば、十分に大きなパルス振幅）に増幅するために用いられる。検出器１１の出力が十分に大きい場合には、増幅器１２を備える必要はない。
　次に、波形整形器１３について説明する。波形整形器１３は、増幅器１２又は検出器１１から出力される信号のノイズ除去、及び増幅率の調整を行う機能を有する回路である。波形整形器１３は、アナログ回路で形成される場合には、一般的に積分回路及び微分回路を組み合わせて構成される。波形整形器１３は、ノイズ除去と増幅率調整とを行い、出力パルスを、波高分析部１４における波高検出が容易になるような形状の波形にする。波形整形器１３は、デジタル回路で形成される場合には、アナログデジタル変換によって前段から出力される信号をデジタル変換し、デジタルフィルタによってノイズ除去と増幅率調整とを行い、波形整形したデジタルでのパルス処理を行う。
　次に、波高分析部１４について説明する。波高分析部１４は、波形整形器１３から出力される出力パルスのピーク値の頻度分布又は出力パルスの時間幅の頻度分布を作成する。波形整形器１３で前述のようなアナログ回路の処理が行われる場合には、波高分析部１４は、出力パルスのピーク値を保持する回路と、アナログデジタル変換を組み合わせて出力パルスの波高を取得する、又は、出力パルスの立上りと立下りの時間を電圧に変換する回路とアナログデジタル変換を組み合わせることで出力パルスの波高を取得する。波形整形器１３で前述のようなデジタル回路の処理が行われる場合には、波高分析部１４は、デジタルフィルタによって処理された出力パルスの振幅又は時間幅から出力パルスの波高を取得する。波高分析部１４は、取得した出力パルスの波高から頻度分布を作成する。
　図２は、放射線検出部１０から出力されるパルス信号の波高（すなわち、エネルギー）分布を放射線線種ごとに示す図である。図２に示されるように、計測される波高分布は、エネルギー分布に変換されても良い。波高分布とエネルギー分布とは、一対一に対応している。例えば、放射線線種は、α線、β線、γ線のうちの１つ以上である。
〈放射線線種判定部２０〉
　次に、放射線線種判定部２０について説明する。放射線線種判定部２０は、波高分析部１４で得られた波高分布又はエネルギー分布から測定対象７０から放出される放射線の種類（すなわち、放射線線種）を検出する機能を有する。
　図２に示されるように、一般的に、放射線検出部１０で検出されるアルファ線（α線）は、α線放出核種から放出され、４ＭｅＶ以上のエネルギーを持つ。例えば、プルトニウム２３８（２３８Ｐｕ）は、５．５０ＭｅＶのエネルギーを持つα線と５．４６ＭｅＶのエネルギーを持つα線とを放出する。プルトニウム２３９（２３９Ｐｕ）は、５．１６ＭｅＶのエネルギーを持つα線を放出する。アメリシウム２４１（２４１Ａｍ）は、５．４９ＭｅＶのエネルギーを持つα線を放出する。キュリウム２４４（２４４Ｃｍ）は、５．９５ＭｅＶのエネルギーを持つα線を放出する。高いエネルギー（一般には、４ＭｅＶ以上のエネルギー）持つα線は、検出器１１によって、図２の実線のように、高いエネルギーの領域に波高分布が形成されるよう検出される。測定対象７０から放出されるα線のエネルギーは、空気中でのエネルギー損失により非対称なピーク構造になるが、波高分布上にピーク構造が現れる特徴を示す。放射線線種判定部２０は、このピーク構造の有無を検知することで、測定対象７０にα線放出核種が含まれるかどうかを判定することができる。
　図２に点線で示されるように、一般的に、放射性核種から放出されるガンマ線（γ線）は、核種特有のエネルギーを持ち、γ線の波高分布は、正規分布に近い形状のピーク構造と検出器１１の内部又は周辺での散乱によってエネルギーを失う過程を反映した連続成分とを有する。一般的な放射性核種のγ線のエネルギーは、３ＭｅＶ以下であり、α線の波高分布のピークより低い。例えば、タリウム（Ｔｌ）は、２．６ＭｅＶのエネルギーを持つγ線を放出する。イットリウム（Ｙ）は、１．８ＭｅＶのエネルギーを持つγ線を放出する。コバルト６０（６０Ｃｏ）は、１．３ＭｅＶと１．１ＭｅＶのエネルギーを持つ２本のγ線を放出する。γ線放出核種から放出されるγ線のエネルギーは、α線放出核種から放出されるα線のエネルギーよりも低く、ピーク構造の形状も異なる（急峻なピーク構造と、連続成分とを有する。）ため、放射線線種判定部２０は、波高分布から測定対象７０にγ線放出核種が含まれるかどうかを判定することができる。
　放射性核種から放出されるベータ線（β線）は、原子核の放射性崩壊の一種であるベータ崩壊の過程で放出される。ベータ崩壊では、電子（すなわち、β線）と反電子ニュートリノとが放出される。ベータ崩壊では、電子と反電子ニュートリノとがエネルギーを分け合うために、放出される電子（すなわち、β線）のエネルギーであるβ線の波高分布は、は、図２に一点鎖線で示されるように、連続的な分布となる。β線の波高分布は、検出器１１の構造によるが、ピーク構造が低いエネルギーの領域に現れるが、α線及びγ線のように特定のエネルギーの領域にピーク構造を持たない。β線の波高分布は、γ線の波高分布と同様に、低いエネルギーの領域で連続的な分布になる。また、低いエネルギーの領域で、β線の波高分布は、γ線の波高分布より大きい。β線の波高分布の形状は、他の放射線線種と異なるので、放射線線種判定部２０は、波高分布から測定対象７０にβ線放出核種が含まれるかどうかを判定することができる。
　波高分布での放射線線種ごとの特徴の検出は、例えば、上記の特徴を解析的に検出する方法、事前に想定した波高分布との一致度検証によって検出する方法、事前に想定した波高分布の機械学習によって放射線線種を判定する方法、教師無し機械学習によって放射線線種を判定するように学習する方法のいずれかによって実現可能である。事前に想定した波高分布は、様々な条件での波高分布を計算することで作成することができる。波高分布の計算には、例えば、放射線挙動解析による計算がある。
　波高分析部１４の波高分布を用いて検出を行う機能は、放射線線種判定部２０が行ってもよい。その場合には、波高分析部１４を備えないことが可能である。
　従来の技術では、放射能分析を行うために、特許文献１のように応答関数を調整することで精度を上げることが提案されてきた。実施の形態１の手法では、測定対象７０から放出される放射線の種類が１種類である場合が想定されており、複数の放射線線種が測定対象７０に含まれる場合には、応答関数の形状が変わることは考慮されていなかった。この結果、応答関数が変わった場合には、単一の逆問題演算アルゴリズムでは、複数の放射線線種が測定される条件で正しく逆問題演算を行うことができなかったため、正しく放射線の強度と放射能の強さを測定することができなかった。
　実施の形態１では、検出器１１から出力される放射線線種を判定し、検出された放射線線種に応じて使用する逆問題演算アルゴリズムを選ぶことで、測定対象７０に複数の放射線線種が含まれる場合においても、それぞれの放射線線種の放射能の強さを求めることが可能となる。
　波高分布での放射線線種ごとの特徴の検出は、例えば、上記の特徴を解析的に検出する方法、事前に想定した波高分布との一致度検証によって検出する方法、事前に想定した波高分布の機械学習によって放射線線種を判定する方法、教師無し機械学習によって放射線線種を判定するように学習することで実現可能である。事前に想定した波高分布は、様々な条件での波高分布を計算することで作成することができる。波高分布の計算には、例えば、放射線挙動解析による計算がある。
〈逆問題演算部５０〉
　次に、逆問題演算部５０について説明する。逆問題演算部５０は、波高分布から測定対象７０に含まれる放射性核種の核種の特定及び放射性核種の放射能の強さ［Ｂｑ］の計算を行う。なお、放射能の強さ［Ｂｑ］は、放射能濃度で表されてもよい。一般的な測定では、測定対象７０に含まれる放射性核種の核種、放射性核種の放射能の強さ［Ｂｑ］、放射能濃度が既知の状態で測定を行い、波高分布を得ることを順問題とされる。放射能分析では、測定対象７０の情報が不明で測定結果である波高分布のみ知ることができるため順問題を逆に解く逆問題を解く必要がある。放射能分析で逆問題演算を行うためには、一般的に、波高分布と、測定対象７０の放射能の強さ又は放射能濃度とを関係づける応答関数が用いられる。波高Ｅの頻度を測定して得られた波高分布Ｍ（Ｅ）は、以下の式（１）で求めることができる。
Figure JPOXMLDOC01-appb-M000001
　ここで、Ｒ（Ｅ）は応答関数であり、Ｓ（Ｅ）は測定対象７０のエネルギーごとの、放射能の強さ（又は放射能濃度）である。この場合に、逆問題演算を行うためのアルゴリズムとしては、例えば、逆行列演算、擬似逆行列演算、デコンボリューション、逆畳み込み積分、及び逐次近似法がある。応答関数は、例えば、放射線挙動解析のように放射線の物理プロセスを反映した演算をすることができるものを用いることで用意することができる。また、式（１）のような解析的な演算ではなく、様々な条件のＭ（Ｅ）とＳ（Ｅ）の組み合わせから、Ｍ（Ｅ）からＳ（Ｅ）を推定できるように学習させた機械学習を用いることでも演算が実現できる。Ｍ（Ｅ）とＳ（Ｅ）は、上記と同様に、放射線挙動解析を用いることで作成することができる。
　上記のように、式（１）のように解析的な関係ではなく、既知のＳ（Ｅ）とＭ（Ｅ）との関係から機械学習で用いられるネットワーク構造によるＳ（Ｅ）とＭ（Ｅ）の関係から逆問題演算を解くこともできる。
〈逆問題演算アルゴリズム選択部３０〉
　次に、逆問題演算アルゴリズム選択部３０について説明する。上記のように逆問題演算を行うことができる複数の逆問題演算アルゴリズムが予め用意されている。演算は、応答関数Ｒ（Ｅ）である行列の性質によって解を求めることができない場合もある。応答関数Ｒ（Ｅ）は、放射線線種によって異なる。
　以下に、式（２）として、γ線の応答関数Ｒ（Ｅ）を示す。以下の応答関数Ｒ（Ｅ）は、対角成分が０より大きい値を持つ行列で表される。
Figure JPOXMLDOC01-appb-M000002
　以下に、式（３）として、α線の応答関数Ｒ（Ｅ）を示す。以下の応答関数Ｒ（Ｅ）は、対角成分が０である行列で表される。
Figure JPOXMLDOC01-appb-M000003
　γ線のように、空気中でのエネルギーの連続的な減衰が無い放射線線種についての応答関数は、式（２）のように、対角成分を持つ（すなわち、対角成分が０ではない）行列が応答関数となる。このような応答関数を持つ放射線線種に対して逆問題演算を行う場合には、逆行列演算、擬似逆行列演算、又は逐次近似法などの解析的な演算アルゴリズムを用いることができる。
　α線のように、空気中での減衰がある放射線線種についての応答関数は、式（３）のように、対角成分を持たない行列となる。例えば、対角成分を持たない行列で応答関数が表される場合には、逆行列が存在しないため逆行列演算は用いることができない。また、擬似逆行列演算は、０以下の解が存在するなど演算精度に問題があることがある。
　β線の逆問題演算を行う場合には、応答関数を想定される放射性核種からの推定スペクトルで作ることができる。この場合には、応答関数が正方行列では無い条件となり、逆行列が存在しないため、解析的な逆問題演算よりも機械学習によるネットワークによる逆問題演算が最適な結果を推定できる場合がある。逆問題演算アルゴリズム選択部３０は、検出した放射線線種に基づいて解の収束又は演算速度の速くなるように１つ又は複数の演算アルゴリズムを検出した放射線線種に応じて最適な選択を行う。
　逆問題演算アルゴリズムＤＢ３１は、上記のように逆問題演算部５０で用いることができる複数の演算アルゴリズムを格納している。
〈応答関数データベース４１〉
　次に、応答関数データベース４１について説明する。応答関数ＤＢ４１は、逆問題演算に用いる応答行列又は機械学習のネットワーク構造のパラメータを持つ。解析的に逆問題演算を行う場合には、式（２）及び式（３）のような行列を保持する。また、解析的な演算以外の方法である機械学習を用いた逆問題演算を行うためのネットワークに関する情報を持っても良い。
〈表示部６０〉
　逆問題演算で導出された測定対象７０に含まれる放射性核種、放射能の強さ［Ｂｑ］、放射線の強さ、及び放射能濃度は、表示部６０にて表示する。
　図３は、実施の形態１に係る放射能分析装置１のハードウェア構成の例を示す図である。図３のハードウェア構成は、実施の形態２～５に係る放射能分析装置にも適用される。図３に示されるように、放射能分析装置１は、放射線検出部１０と情報処理部１ａとを有している。情報処理部１ａは、ＣＰＵ（Ｃｅｎｔｒａｌ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ）などのプロセッサ１０１と、揮発性の記憶装置であるメモリ１０２と、ハードディスクドライブ（ＨＤＤ）又はソリッドステートドライブ（ＳＳＤ）などの不揮発性の記憶装置１０３とを有している。メモリ１０２は、例えば、ＲＡＭ（Ｒａｎｄｏｍ　Ａｃｃｅｓｓ　Ｍｅｍｏｒｙ）などの半導体メモリである。記憶装置１０３は、図１に示されるデータベースを記憶するものであってもよい。
　情報処理部１ａの各機能は、例えば、処理回路により実現される。処理回路は、専用のハードウェアであってもよく、又はメモリ１０２に格納されるプログラムを実行するプロセッサ１０１であってもよい。プロセッサ１０１は、処理装置、演算装置、マイクロプロセッサ、マイクロコンピュータ、及びＤＳＰ（Ｄｉｇｉｔａｌ　Ｓｉｇｎａｌ　Ｐｒｏｃｅｓｓｏｒ）のいずれであってもよい。
　処理回路が専用のハードウェアである場合、処理回路は、例えば、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ　Ｓｐｅｃｉｆｉｃ　Ｉｎｔｅｇｒａｔｅｄ　Ｃｉｒｃｕｉｔ）、ＦＰＧＡ（Ｆｉｅｌｄ－Ｐｒｏｇｒａｍｍａｂｌｅ　Ｇａｔｅ　Ａｒｒａｙ）、又はこれらのうちのいずれかを組み合わせたものである。
　処理回路がプロセッサ１０１である場合、情報処理部１ａによって実行される放射能分析プログラムは、ソフトウェア、ファームウェア、又はソフトウェアとファームウェアとの組み合わせにより実現される。放射能分析プログラムは、ネットワークを経由して又は記録媒体から情報処理部１ａにインストールされる。ソフトウェア及びファームウェアは、プログラムとして記述され、メモリ１０２に格納される。プロセッサ１０１は、メモリ１０２に記憶された表示制御プログラムを読み出して実行することにより、図１に示される各部の機能を実現することができる。
　なお、情報処理部１ａは、一部を専用のハードウェアで実現し、一部をソフトウェア又はファームウェアで実現するようにしてもよい。このように、処理回路は、ハードウェア、ソフトウェア、ファームウェア、又はこれらのうちのいずれかの組み合わせによって、上述の各機能を実現することができる。
　図４は、実施の形態１に係る放射能分析方法を示すフローチャートである。この放射能分析方法は、情報処理部１ａとしてのコンピュータによって実行される。この放射能分析方法は、測定対象７０から放出される放射線を検出し、放射線に基づく第１の信号Ｄ１を出力する放射線検出部１０から第１の信号Ｄ１を受け取るステップＳ１と、第１の信号に基づいて放射線の種類を示す放射線線種を判定し、放射線線種を示す情報を含む第２の信号Ｄ２を出力するステップＳ２と、記憶装置に予め保存されている複数の逆問題演算アルゴリズムから放射線線種に基づいて逆問題演算アルゴリズムを選択するステップＳ３と、選択された逆問題演算アルゴリズムと記憶装置に予め保存されている前記放射線検出部１０の応答関数とを用いる逆問題演算を実行することによって、第１の信号Ｄ１と前記第２の信号Ｄ２とから測定対象７０に含まれる放射性核種の特定と前記放射性核種の放射能の強さの計算とを行うステップＳ４とを有している。また、ステップＳ４において、放射線の強さ（例えば、放射線のフルエンス）の計算を行ってもよい。
〈効果〉
　以上に説明したように、実施の形態１に係る放射能分析装置１及び放射能分析方法によれば、複数の逆問題演算アルゴリズムから、適した逆問題演算アルゴリズムを選択して使用することで、逆問題演算の精度を上げることができる。
《２》実施の形態２
　図５は、実施の形態２に係る放射能分析装置２の構成を概略的に示すブロック図である。図５において、図１に示される構成と同一又は対応する構成には、図１に示される符号と同じ符号が付されている。実施の形態２に係る放射能分析装置２は、情報処理部２ａが応答関数選択部４０を有する点において、実施の形態１に係る放射能分析装置１と異なる。
　放射能分析装置２は、測定対象７０から放出される放射線を検出し、検出された放射線に基づく第１の信号Ｄ１を出力する放射線検出部１０と、測定対象７０に含まれる放射性物質を示す放射性核種の特定と放射性核種の放射能の強さ［Ｂｑ］の計算とを行う情報処理部２ａとを有している。情報処理部２ａは、例えば、コンピュータであり、実施の形態２に係る放射能分析方法を実施することができる。
　実施の形態２においては、記憶装置（例えば、図３の記憶装置１０３）に応答関数ＤＢ４１として予め複数の応答関数が応答関数ＤＢ４１として保存されている。応答関数選択部４０は、放射線線種判定部２０において判定された放射線線種を示す第２の信号Ｄ２に基づいて、応答関数ＤＢ４１に保存されている複数の応答関数から、使用する応答関数を選択する。逆問題演算部５０は、逆問題演算アルゴリズム選択部３０によって選択された逆問題演算アルゴリズムと応答関数選択部４０によって選択された応答関数とを用いて、逆問題演算を実行することによって、第１の信号Ｄ１と第２の信号Ｄ２とから測定対象７０に含まれる放射性核種の特定と放射性核種の放射能の強さ［Ｂｑ］の計算とを行う。また、逆問題演算部５０は、逆問題演算を実行することによって、第１の信号Ｄ１と第２の信号Ｄ２とから放射線の強さの計算を行ってもよい。
　以上に説明したように、実施の形態２に係る放射能分析装置２及び放射能分析方法によれば、複数の逆問題演算アルゴリズムから、適した逆問題演算アルゴリズムを選択し、複数の応答関数から適した応答関数を選択し、これらを使用することで、逆問題演算の精度を上げることができる。
　上記以外に関し、実施の形態２は、実施の形態１と同じである。
《３》実施の形態３
　図６は、実施の形態３に係る放射能分析装置３の構成を概略的に示すブロック図である。図６において、図１に示される構成と同一又は対応する構成には、図１に示される符号と同じ符号が付されている。実施の形態３に係る放射能分析装置３は、放射線検出部１０ｂが、検出器１１によって検出された放射線の放射線種別ごとに検出信号を弁別して弁別信号を出力する波形弁別器１５とを有する点において、実施の形態１に係る放射能分析装置１と異なる。波形弁別器１５は、検出器１１の出力及び増幅器１２から出力する出力信号のパルス形状が放射線線種によって差が出るようにする回路である。
　放射線線種判定部２０ｂは、弁別信号を含む第１の信号Ｄ１に基づいて放射線線種を判定し、判別された放射線線種を示す情報を含む第２の信号Ｄ２を出力する。放射能分析装置３は、測定対象７０から放出される放射線を検出し、放射線の波高分布を示す弁別信号を含む第１の信号Ｄ１を出力する放射線検出部１０ｂと、測定対象７０に含まれる放射性核種の特定と放射性核種の放射能の強さ［Ｂｑ］の計算とを行う情報処理部３ａとを有している。情報処理部３ａは、例えば、コンピュータであり、実施の形態３に係る放射能分析方法を実施することができる。
　図７は、放射線検出部１０ｂの検出器１１から出力されるパルス信号の振幅の立ち上がり及び減衰特性を放射線線種ごとに示す図である。例えば、検出器１１の、放射線と相互作用する物質として半導体を用いた場合、半導体の内部では、放射線の通過によって電子と正孔とが発生するが、電子の移動度と正孔の移動度とは異なる。α線のように検出器１１の表面付近で反応する放射線線種では、正孔の移動距離を長くするように検出器１１を設定することでパルスの立上りが緩やかな形状になる（図７に実線で示される）。また、γ線及びβ線のように検出器１１の内部（表面よりも深い位置）で反応する放射線線種では、立上りが急峻な形状となる（図７に破線で示される）。波形弁別器１５は、このような波形形状の違いを用いて、放射線線種を弁別することで特定の放射線線種を精度良く検出できる。
　実施の形態３においては、逆問題演算部５０は、逆問題演算アルゴリズム選択部３０によって選択された逆問題演算アルゴリズムと、応答関数とを用いて、逆問題演算を実行することによって、第１の信号Ｄ１と第２の信号Ｄ２とから測定対象７０に含まれる放射性核種の特定と放射性核種の放射能の強さ［Ｂｑ］の計算とを行う。また、逆問題演算部５０は、逆問題演算を実行することによって、第１の信号Ｄ１と第２の信号Ｄ２とから放射線の強さの計算を行ってもよい。
　以上に説明したように、実施の形態３に係る放射能分析装置３及び放射能分析方法によれば、複数の逆問題演算アルゴリズムから、適した逆問題演算アルゴリズムを選択し使用することで、逆問題演算の精度を上げることができる。
　また、波形弁別器１５によって特定の放射線線種を判定し、判定の結果を放射線線種判定部２０での波高分布からの検出を組み合わせることで、検出精度を上げることができる。また、特定の放射線線種の検出結果と波高分布とを対応付けることで特定の放射線線種の波高分布を取り出すことで波高分布を分離して逆問題演算を解くことができる。
　上記以外に関し、実施の形態３は、実施の形態１又は２と同じである。
《４》実施の形態４
　図８は、実施の形態４に係る放射能分析装置４の構成を概略的に示すブロック図である。図８において、図１に示される構成と同一又は対応する構成には、図１に示される符号と同じ符号が付されている。実施の形態４に係る放射能分析装置４は、放射線検出部１０ｃの検出器１１ｃが放射線によってエネルギー付与された位置を検出する機能を有する点及び放射線検出部１０ｃが軌跡形状検出部１６を有する点において、実施の形態１に係る放射能分析装置１と異なる。
　放射能分析装置４は、測定対象７０から放出される放射線を検出し、検出された放射線に基づく第１の信号Ｄ１と放射線の軌跡形状を示す軌跡信号を出力する放射線検出部１０ｃと、測定対象７０に含まれる放射性核種の特定と放射性核種の放射能の強さ［Ｂｑ］の計算とを行う情報処理部４ａとを有している。情報処理部４ａは、例えば、コンピュータであり、実施の形態４に係る放射能分析方法を実施することができる。
　放射線検出部１０ｃは、放射線と相互作用する物質を有する検出器１１ｃと、相互作用する物質内において放射線がエネルギーを付与した軌跡の位置を検出し、軌跡の形状を示す軌跡形状信号を出力する回路である軌跡形状検出部１６とを有している。放射線線種判定部２０ｃは、軌跡形状信号を含む第１の信号Ｄ１に基づいて放射線線種を判定し、判定された放射線線種を示す情報を含む第２の信号Ｄ２を出力する。
　実施の形態４では、検出器１１ｃが放射線によってエネルギー付与された位置を検出する機能を有する。例えば、検出器１１ｃは、アレイ状に複数のシンチレータ結晶を並べ、アレイ状の並べられた複数のシンチレータ結晶の各々の上に光検出器を配置した装置である（例えば、特許文献２参照）。また、検出器１１ｃは、規則的に配列された複数のシンチレータ結晶と、シンチレータ結晶で発生した光を撮影するカメラ（例えば、ＣＭＯＳ（Ｃｏｍｐｌｅｍｅｎｔａｒｙ　Ｍｅｔａｌ　Ｏｘｉｄｅ　Ｓｅｍｉｃｏｎｄｕｃｔｏｒ）カメラ素子、ＳＯＩ（Ｓｉｌｉｃｏｎ－Ｏｎ－Ｉｎｓｕｌａｔｏｒ）撮像素子、など）とからなる装置であってもよい。また、検出器１１ｃは、半導体放射線検出器（例えば、Ｇｅ結晶体）と、これを挟むように、且つ、間隔を開けて規則的に配列された複数のストリップ電極とを持つ装置であってもよい。また、検出器１１ｃは、複数のワイヤーを持つガス式放射線検出器又は規則的に配列された複数のガス式放射線検出器によって放射線によってエネルギー付与された位置を検出する装置であってもよい。
特表２０１６－５１４８３５号公報
　図９（Ａ）から（Ｃ）は、放射線検出部１０ｃの検出器１１ｃにおける軌跡形状を放射線線種ごとに示す概略図である。図９（Ａ）から（Ｃ）は、放射線線種による検出器１１ｃ内における軌跡形状（横から見た場合と上から見た場合）を示す。図９（Ａ）に示されるように、α線のように飛程が短く付与エネルギーが大きい放射線は、広がった点状の形状が軌跡形状となる。図９（Ｂ）に示されるように、β線は、飛程が長い特性により、線状に検出器を通過するため細長い形状になる。図９（Ｃ）に示されるように、γ線は、一点で反応するため広がりのない点状の形状となる。この形状の違いを放射線線種判定部２０ｃで検出することで、測定対象７０から放出される放射線に含まれる放射線線種を判定し、第２の信号Ｄ２を生成できる。
　以上に説明したように、実施の形態４に係る放射能分析装置４及び放射能分析方法によれば、複数の逆問題演算アルゴリズムから、適した逆問題演算アルゴリズムを選択し使用することで、逆問題演算の精度を上げることができる。
　また、軌跡形状検出部１６によって特定の放射線線種を判定し、判定の結果を放射線線種判定部２０での波高分布からの検出を組み合わせることで、検出精度を上げることができる。また、特定の放射線線種の検出結果と波高分布とを対応付け特定の放射線線種の波高分布を取り出すことで波高分布を分離して逆問題演算を解くことができる。
　上記以外に関し、実施の形態４は、実施の形態１又は２と同じである。
《５》実施の形態５
　図１０は、実施の形態５に係る放射能分析装置５の構成を概略的に示すブロック図である。図１０において、図１に示される構成と同一又は対応する構成には、図１に示される符号と同じ符号が付されている。実施の形態５に係る放射能分析装置５は、放射線検出部１０ｄが位置検出型検出器１１ｄで構成されている点において、実施の形態１に係る放射能分析装置１と異なる。
　放射能分析装置５は、測定対象７０から放出される放射線を検出し、検出された放射線に基づく第１の信号として、放射線線種に応じた輝度値と軌跡の拡がりと示す信号Ｄ５を出力する放射線検出部１０ｄと、測定対象７０に含まれる放射性核種の特定と放射性核種の放射能の強さ［Ｂｑ］の計算とを行う情報処理部５ａとを有している。情報処理部５ａは、例えば、コンピュータであり、実施の形態５に係る放射能分析方法を実施することができる。
　放射線検出部１０ｄは、入射する放射線と相互作用する物質を備え、相互作用する物質内において放射線が物質にエネルギーを付与した位置と物質に付与したエネルギーの量を示す信号を第１の信号Ｄ５として出力する位置検出型検出器１１ｄを有している。放射線線種判定部２０ｄは、物質にエネルギーを付与した位置と物質に付与したエネルギーの量を示す信号とを含む第１の信号Ｄ５に基づいて放射線線種を判定し、判定された放射線線種を示す情報を含む第２の信号Ｄ２を出力する。
　位置検出型検出器１１ｄは、入射する放射線によってエネルギー付与された位置を検出する機能を有し、付与されたエネルギーに対応する振幅を持つパルスを出力する。或いは、位置検出型検出器１１ｄは、入射する放射線によってエネルギー付与された位置を検出する機能を有し、エネルギー付与された位置で発生する光の輝度値を出力する機能を有する。位置検出型検出器１１ｄは、エネルギー付与された位置とパルス振幅の情報を出力する。或いは、位置検出型検出器１１ｄは、エネルギー付与された位置と輝度値情報とを含む画像を出力する。放射線線種判定部２０ｄは、位置検出型検出器１１ｄから出力された位置情報とパルス波高値、又は位置情報と輝度値とを含む画像に基づいて、放射線線種を判定する。位置検出型検出器１１ｄは、実施の形態４における検出器１１ｃと同様の方法で実現可能である。
　図１１は、放射線検出部１０ｄの位置検出型検出器１１ｄにおける輝度値と軌跡形状の広がりとの関係を示す図である。図１１に示されるように、放射線線種によって、輝度値と軌跡形状の広がりとの関係が異なる。α線軌跡における発光の輝度値は、他の放射線線種のものよりも高く、輝度値の広がり（図１１の横軸方向の幅）は、他の放射線線種のものよりも大きい。β線軌跡の広がりは、他の放射線線種のものよりも高く、β線軌跡における発光の輝度値は、他の放射線線種のものよりも低い。γ線軌跡の広がりは、他の放射線線種のものよりも低く、γ線軌跡における発光の輝度値は、α線軌跡におけるものよりも低い。これらの違いを放射線線種判定部２０ｄで検出することで、放射線線種判定部２０ｄは、測定対象７０から放出される放射線に含まれる放射線線種を判定し、第２の信号Ｄ２を生成できる。
　以上に説明したように、実施の形態５に係る放射能分析装置５及び放射能分析方法によれば、複数の逆問題演算アルゴリズムから、適した逆問題演算アルゴリズムを選択し使用することで、逆問題演算の精度を上げることができる。
　また、放射線線種の判定精度を上げることができ、特定の放射線線種の検出結果と波高分布とを対応付けることで特定の放射線線種の波高分布を取り出すことで波高分布を分離して逆問題演算を解くことができる。
　上記以外に関し、実施の形態５は、実施の形態１又は２と同じである。
《６》変形例
　実施の形態２で説明した、応答関数選択部４０を実施の形態３、４、５のいずれかに組み合わせることが可能である。
　１～５　放射能分析装置、　１ａ～５ａ　情報処理部、　１０、１０ｂ、１０ｃ、１０ｄ　放射線検出部、　１１、１１ｂ、１１ｃ　検出器、　１１ｄ　位置検出型検出器、　１２　増幅器、　１３　波形整形器、　１４　波高分析部、　１５　波形弁別器、　１６　軌跡形状検出部、　２０、２０ｂ、２０ｃ、２０ｄ　放射線線種判定部、　３０　逆問題演算アルゴリズム選択部、　３１　逆問題演算アルゴリズムＤＢ、　４０　応答関数選択部、　４１　応答関数ＤＢ、　５０　逆問題演算部、　６０　表示部、　７０　測定対象。

Claims (8)

  1. 　測定対象から放出される放射線を検出し、前記放射線に基づく第１の信号を出力する放射線検出部と、
    　前記第１の信号に基づいて前記放射線の種類を示す放射線線種を判定し、前記放射線線種を示す情報を含む第２の信号を出力する放射線線種判定部と、
    　記憶装置に予め保存されている複数の逆問題演算アルゴリズムから前記放射線線種に基づいて逆問題演算アルゴリズムを選択する逆問題演算アルゴリズム選択部と、
    　前記選択された逆問題演算アルゴリズムと記憶装置に予め保存されている前記放射線検出部の応答関数とを用いる逆問題演算を実行することによって、前記第１の信号と前記第２の信号とから前記測定対象に含まれる放射性核種の特定と前記放射性核種の放射能の強さの計算とを行う逆問題演算部と
    　を有することを特徴とする放射能分析装置。
  2. 　応答関数選択部をさらに有し、
    　前記予め保存されている前記応答関数として、複数の応答関数が保存されており、
    　前記応答関数選択部は、前記複数の応答関数から前記第２の信号が示す放射線線種に基づいて応答関数を選択し、
    　前記逆問題演算部は、前記選択された逆問題演算アルゴリズムと前記選択された応答関数とを用いて前記逆問題演算を実行することによって、前記第１の信号と前記第２の信号とから前記測定対象に含まれる放射性核種の特定と前記放射性核種の放射能の強さの計算とを行う
    　ことを特徴とする請求項１に記載の放射能分析装置。
  3. 　前記逆問題演算部は、前記逆問題演算を実行することによって、前記第１の信号と前記第２の信号とから前記放射線の強さの計算を行う
    　ことを特徴とする請求項１又は２に記載の放射能分析装置。
  4. 　前記放射線検出部は、前記放射線を検出して検出信号を出力する検出器と、前記放射線の放射線種別ごとに前記検出信号を弁別して弁別信号を出力する波形弁別器とを有し、
    　前記放射線線種判定部は、前記弁別信号を含む前記第１の信号に基づいて前記放射線線種を判定し、前記放射線線種を示す情報を含む前記第２の信号を出力する
    　ことを特徴とする請求項１から３のいずれか１項に記載の放射能分析装置。
  5. 　前記放射線検出部は、前記放射線と相互作用する物質を有する検出器と、前記物質内において前記放射線がエネルギーを付与した軌跡の位置を検出し、前記軌跡の形状を示す軌跡形状信号を出力する軌跡形状検出部とを有し、
    　前記放射線線種判定部は、前記軌跡形状信号を含む前記第１の信号に基づいて前記放射線線種を判定し、前記放射線線種を示す情報を含む前記第２の信号を出力する
    　ことを特徴とする請求項１から３のいずれか１項に記載の放射能分析装置。
  6. 　前記放射線検出部は、前記放射線と相互作用する物質を備え、前記物質内において前記放射線が前記物質にエネルギーを付与した位置と前記物質に付与した前記エネルギーの量を示す信号を前記第１の信号として出力する位置検出型検出器を有し、
    　前記放射線線種判定部は、前記物質にエネルギーを付与した位置と前記物質に付与した前記エネルギーの量とを示す信号を含む前記第１の信号に基づいて前記放射線線種を判定し、前記放射線線種を示す情報を含む前記第２の信号を出力する
    　ことを特徴とする請求項１から３のいずれか１項に記載の放射能分析装置。
  7. 　前記放射線線種は、α線、β線、γ線のうちの１つ以上である
    　ことを特徴とする請求項１から６のいずれか１項に記載の放射能分析装置。
  8. 　コンピュータによって実行される放射能分析方法であって、
    　測定対象から放出される放射線を検出し、前記放射線に基づく第１の信号を出力する放射線検出部から前記第１の信号を受け取るステップと、
    　前記第１の信号に基づいて前記放射線の種類を示す放射線線種を判定し、前記放射線線種を示す情報を含む第２の信号を出力するステップと、
    　記憶装置に予め保存されている複数の逆問題演算アルゴリズムから前記放射線線種に基づいて逆問題演算アルゴリズムを選択するステップと、
    　前記選択された逆問題演算アルゴリズムと記憶装置に予め保存されている前記放射線検出部の応答関数とを用いる逆問題演算を実行することによって、前記第１の信号と前記第２の信号とから前記測定対象に含まれる放射性核種の特定と前記放射性核種の放射能の強さの計算とを行うステップと
    　を有することを特徴とする放射能分析方法。
PCT/JP2022/034696 2022-09-16 2022-09-16 放射能分析装置及び放射能分析方法 WO2024057514A1 (ja)

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)

* Cited by examiner, † Cited by third party
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 クロメック リミテッドＫｒｏｍｅｋ Ｌｉｍｉｔｅｄ スペクトルデータの検出・操作方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015537215A (ja) * 2012-11-23 2015-12-24 クロメック リミテッドＫｒｏｍｅｋ Ｌｉｍｉｔｅｄ スペクトルデータの検出・操作方法
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
EP2773982A1 (en) Neutron spectrometer
JP4061367B2 (ja) ＺｎＳ（Ａｇ）シンチレーション検出器
Parks The Compton effect-Compton scattering and gamma ray spectroscopy
JP3597973B2 (ja) ダストモニタ
JPH01260388A (ja) エネルギー補正オフセットにより像を一様にしたガンマカメラ
JP6524484B2 (ja) 放射線計測方法及び放射線計測装置
JPH05341047A (ja) 効果的なα及びβ（γ）線同時測定法及びその検出器
JP7085964B2 (ja) 放射線検出器
WO2024057514A1 (ja) 放射能分析装置及び放射能分析方法
JPH11194170A (ja) 放射性物質検査装置及び放射性廃棄物検査システム
JP2002006055A (ja) 放射線測定方法および装置
KR102663201B1 (ko) 방사선 측정 장치 및 장치의 동작 방법
KR102280128B1 (ko) 방사성 핵종 판별 장치 및 방법
JPH09211133A (ja) 放射線モニタ
JP2012103179A (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