CN110082811B - Energy track counting reassignment process and method for gamma spectrum analysis - Google Patents
Energy track counting reassignment process and method for gamma spectrum analysis Download PDFInfo
- Publication number
- CN110082811B CN110082811B CN201810072826.7A CN201810072826A CN110082811B CN 110082811 B CN110082811 B CN 110082811B CN 201810072826 A CN201810072826 A CN 201810072826A CN 110082811 B CN110082811 B CN 110082811B
- Authority
- CN
- China
- Prior art keywords
- energy
- line segment
- channel
- count
- interval
- Prior art date
- Legal status (The legal status 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 status listed.)
- Active
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/36—Measuring spectral distribution of X-rays or of nuclear radiation spectrometry
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Measurement Of Radiation (AREA)
Abstract
The invention relates to an energy track counting reassignment process and method for gamma spectrum analysis, comprising the following steps: determining two mutually parallel line segments with equal length, and respectively naming the two line segments as AB and MN, wherein MN represents an integer energy channel, AB represents a non-integer energy channel and is original data; the method comprises the following steps of (1) partitioning an AB line segment into a plurality of energy intervals, selecting a midpoint in the middle of each energy interval, wherein end points and the midpoint of each interval are nodes, and the nodes are marked by letters; partitioning the MN line segment into a plurality of target intervals, and adopting letter identification at the end points of the target intervals; and selecting the end point of the corresponding target interval as a vertical line towards the line segment AB, adopting letter identification for the foot on the line segment AB, obtaining the length ratio in the energy channel where the foot is located through the position of the foot, and calculating the integral counting of the energy channel in the target interval. And the integral operation is carried out on the non-integral energy channel address, so that the gamma spectrum analysis is conveniently carried out by a deconvolution iterative algorithm based on a response matrix, and the gamma spectrum analysis efficiency is improved.
Description
Technical Field
The invention relates to the fields of gamma spectrometers, gamma energy spectrum analysis, in-situ measurement and in-water radioactivity monitoring, in particular to an energy channel counting redistribution process and method for gamma spectrum analysis.
Background
In the process of gamma spectrometer energy spectrum analysis, the spectrometer is usually calibrated for energy. After calibration, since the energy calibration coefficients are usually fit by polynomial, the coefficients are usually floating point numbers which cannot be predicted in advance, and thus the addresses are changed from integer numbers to floating point numbers through the energy calibration. It is not beneficial to the gamma spectrum analysis based on the deconvolution iterative algorithm of the response matrix, nor to the comparison of the monte carlo simulation result with the experimental result, because the abscissa of the typical monte carlo simulation result is usually the integer energy address boundary. Therefore, after the pulse height spectrum measured by experiments passes through the energy scale, the integer operation is carried out on the non-integer energy channel address, the abscissa of the energy spectrum data after the energy scale is non-integer energy, the non-integer energy channel is converted into the integer energy, and meanwhile, the full spectrum counting rate and the spectrum shape information are guaranteed not to be lost.
Disclosure of Invention
The invention aims to solve the defects in the prior art and provides an energy track counting reassignment process and method for gamma spectrum analysis.
In order to achieve the purpose, the invention adopts the technical scheme that:
a trace count reassignment process and method for gamma spectrum analysis comprising the steps of:
determining two parallel line segments with equal length, and respectively naming the two line segments as AB and MN, wherein MN represents an integer energy channel, AB represents a non-integer energy channel and is original data;
the method comprises the following steps of (1) partitioning an AB line segment into a plurality of energy intervals, selecting a midpoint in the middle of each energy interval, wherein end points and the midpoint of each interval are nodes, and the nodes are marked by letters;
partitioning the MN line segment into a plurality of target intervals, and adopting letter identification at the end points of the target intervals;
and selecting the end point of the corresponding target interval as a vertical line towards the line segment AB, adopting letter identification for the foot on the line segment AB, obtaining the length ratio in the energy channel where the foot is located through the position of the foot, and calculating the integral counting of the energy channel in the target interval.
Further, four nodes D, F, H and J are selected from the line segment AB, the line segment AB is divided into five energy intervals of AD, DF, FH, HJ and JB, and C, E, G, I and K are selected as the midpoints of the five energy intervals respectively;
two nodes O and P are selected on the line segment MN, and the MN is divided into three target intervals of MO, OP and PN.
Furthermore, the nodes O and P respectively make vertical lines towards the line segment AB, and the vertical feet on the line segment AB are respectively W and X.
Furthermore, the energy channel of the target interval MO is an AW energy channel on the AB line segment, and the energy channel of the MP is an AX energy channel on the AB line segment.
Further, the total count within the target interval OP = total count within MP-total count within MO.
Furthermore, the positions of nodes O and P on the MN and four nodes D, F, H and J on the AB are selected randomly.
Furthermore, the MO counting method of the target interval is to directly count the counts in the complete original energy interval AD and add the length ratio of the DW in the energy channel DF to the count in the channel, i.e. the MO counting method of the target interval is to multiply the counts in the energy channel DF by the length ratio of the DW in the energy channel DF, i.e. the MO counting method of the target interval
Can count in the channel.
Further, the target interval PN total count = MN inner total count-MP inner total count.
Further, the method for calculating the total count in the MP is the same as the method for calculating the total count in the MO.
The invention has the beneficial effects that: according to the invention, through the research on the non-integer energy track caused by energy scales in the energy spectrum analysis, the basic property of an integer is fully utilized as a starting point, and through the redistribution of the track address and the counting, a set of effective energy track integer process and method are realized, the integer operation is carried out on the non-integer energy track address, the gamma spectrum analysis is conveniently carried out based on the deconvolution iterative algorithm of a response matrix, the comparison between a Monte Carlo simulation result and an experimental result is facilitated, and the gamma spectrum analysis efficiency is improved.
Drawings
FIG. 1 is a schematic diagram of a geometric scale-normalization method for reassignment of energy track counts according to the present invention;
FIG. 2 is a schematic diagram of an intra-lane count reassignment implementation of the present invention;
FIG. 3 is a comparison graph of a test of energy spectrum of an underwater gamma spectrometer with energy trace count redistribution.
Detailed Description
The invention is further described with reference to the following figures and detailed description.
As shown in fig. 1 and fig. 2, a trace count reassignment process and method for gamma spectrum analysis comprises the following steps:
determining two parallel line segments with equal length, and respectively naming the two line segments as AB and MN, wherein MN represents an integer energy channel, AB represents a non-integer energy channel and is original data;
the method comprises the following steps of (1) partitioning an AB line segment into a plurality of energy intervals, selecting a midpoint in the middle of each energy interval, wherein end points and the midpoint of each interval are nodes, and the nodes are marked by letters;
partitioning the MN line segment into a plurality of target intervals, and adopting letter identification at the end points of the target intervals;
and selecting the end point of the corresponding target interval as a vertical line towards the line segment AB, adopting letter identification for the foot on the line segment AB, obtaining the length ratio in the energy channel where the foot is located through the position of the foot, and calculating the integral counting of the energy channel in the target interval.
The specific operation is as follows:
selecting four nodes D, F, H and J on a line segment AB, dividing the line segment AB into five energy intervals of AD, DF, FH, HJ and JB, and selecting C, E, G, I and K as the midpoints of the five energy intervals respectively;
two nodes O and P are selected on a line segment MN, the MN is divided into three target intervals MO, OP and PN, the nodes O and P respectively make vertical lines towards the line segment AB, and the vertical legs on the line segment AB are respectively W and X.
The positions of nodes O and P on the MN and four nodes D, F, H and J on the AB are selected randomly, the energy channel of the target interval MO is an AW energy channel on the AB line segment, and the energy channel of the MP is an AX energy channel on the AB line segment.
For the interval MO, AD belongs to a complete energy interval, while DW is less than a complete energy interval. For interval MP, there are three complete energy intervals AD, DF, FH, and HX belongs to the incomplete energy region. Fig. 1 shows the case where the original data interval is smaller than an integer. Cases where the data interval is greater than an integer may also be considered similarly.
The primary purpose of the integer is to redistribute the ordinate counts of the spectral data while maintaining the full spectral count and spectral shape information. For the leftmost new energy channel MO, the count should be made by directly counting the counts in the complete original energy interval AD and adding the length ratio of DW in its energy channel DF multiplied by its in-channel count, i.e. the DW length in its energy channel DF is counted
Can count in the channel.
Thus, the process of one-time energy channel integer counting reallocation is completed. For non-left-end intervals such as OP, since the left end O rarely coincides with the original energy track boundary, the subtraction method is used, i.e. total count in OP = total count in MP-total count in MO. The calculation method of the total count in the MP is similar to that of the total count in the MO, and is not described in detail.
Writing a program according to the block diagram shown in fig. 2 can realize the process of counting redistribution in the energy channels, and mainly uses the above geometric proportion principle;
further, the energy spectrum of the simulated underwater gamma spectrometer is tested by energy channel counting redistribution, and a comparison graph before and after the test is shown in fig. 3:
where the solid line is the simulated original spectrum and the abscissa is an integer, but becomes non-integer after the energy calibration, the energy traces are re-converted to integer counts using the method and process described above. The full spectrum count rate difference is only about 0.6%. There is no change in energy resolution. The full spectrum counting rate and the spectrum shape information are ensured.
The energies and counts taken over the energy range 550keV to 572keV are given in the following table:
processing front and back partial data
It can be seen that the new energy interval counts are accumulated indeed, ensuring that the total count rate is unchanged.
According to the invention, through the research on the non-integer energy track caused by energy scales in the energy spectrum analysis, the basic property of an integer is fully utilized as a starting point, and through the redistribution of the track address and the counting, a set of effective energy track integer process and method are realized, the integer operation is carried out on the non-integer energy track address, the gamma spectrum analysis is conveniently carried out based on the deconvolution iterative algorithm of a response matrix, the comparison between a Monte Carlo simulation result and an experimental result is facilitated, and the gamma spectrum analysis efficiency is improved.
It will be apparent to those skilled in the art that various changes and modifications may be made in the present invention without departing from the spirit and scope of the invention. Thus, if such modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is intended to include such modifications and variations.
Claims (9)
1. A trace count reassignment process and method for gamma spectrum analysis comprising the steps of:
a. determining two mutually parallel line segments with equal length, and respectively naming the two line segments as AB and MN, wherein MN represents an integer energy channel, AB represents a non-integer energy channel and is original data;
b. dividing the AB line segment into five energy intervals, selecting a midpoint in the middle of the energy intervals, wherein the end points and the midpoint of each interval are nodes, and the nodes are marked by letters;
c. the MN line segment is divided into three target intervals, letter identification is adopted at the end point of the target interval, namely four nodes are selected from the line segment AB, and two nodes are selected from the line segment MN;
d. selecting the end point of the corresponding target interval as a vertical line towards the line segment AB, marking the foot on the line segment AB by letters, obtaining the length ratio in the energy channel where the foot is located through the position of the foot, calculating the integral counting of the energy channel of the target interval, for a new energy channel with the left end point as the leftmost end M, directly counting the whole original energy interval before the foot corresponding to the right end point, adding the length ratio of the line segment between the right end point of the original energy interval and the foot in the whole energy channel where the line segment is located to multiply the in-channel counting, and for the interval with the left end point not being the leftmost end, adopting the method of subtracting the energy channels with the adjacent left end points as the leftmost end M to calculate the counting of the target interval.
2. The energy track counting reassignment process and method for gamma spectrometry according to claim 1, wherein four nodes D, F, H, J are selected on line segment AB, dividing AB into five energy intervals of AD, DF, FH, HJ and JB, and selecting C, E, G, I, K as the midpoints of the five energy intervals respectively;
two nodes O and P are selected on a line segment MN to divide the MN into three target intervals of MO, OP and PN.
3. The energy track count reassignment process and method for gamma spectroscopy according to claim 2 wherein the nodes O, P are respectively perpendicular to line segment AB, where the foot of the perpendicular is W and X respectively.
4. The energy track counting reassignment process and method for gamma spectrum analysis according to claim 3, wherein the energy track of the target interval MO is an AW energy track on AB line segment, and the energy track of MP is an AX energy track on AB line segment.
5. The channel count reassignment process for gamma spectrometry according to claim 4, wherein the total count in target interval OP = MP total count-MO total count.
6. The energy track count reassignment process and method for gamma spectroscopy according to claim 5, wherein the positions of the nodes O, P on MN and the four nodes D, F, H, J on AB are arbitrarily selected.
7. The method of claim 5, wherein the MO count of the target region is directly counted in the complete original energy interval of AD, and the length ratio of DW in its energy channel DF is multiplied by the count in its channel, i.e. DW length ratio is added to the count in its energy channel DF
Can count in the channel.
8. The channel count reassignment procedure and method for gamma spectrometry according to claim 4, wherein target interval PN inner total count = MN inner total count-MP inner total count.
9. The energy track count reassignment process for gamma spectrometry according to claim 7, wherein the total count in MP is calculated in the same way as the total count in MO.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810072826.7A CN110082811B (en) | 2018-01-25 | 2018-01-25 | Energy track counting reassignment process and method for gamma spectrum analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810072826.7A CN110082811B (en) | 2018-01-25 | 2018-01-25 | Energy track counting reassignment process and method for gamma spectrum analysis |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110082811A CN110082811A (en) | 2019-08-02 |
| CN110082811B true CN110082811B (en) | 2023-04-14 |
Family
ID=67412025
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810072826.7A Active CN110082811B (en) | 2018-01-25 | 2018-01-25 | Energy track counting reassignment process and method for gamma spectrum analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110082811B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112906193B (en) * | 2021-01-21 | 2023-05-12 | 兰州大学 | Gaussian broadening and spectrum resolving method of gamma response function of detector |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| NO954005D0 (en) * | 1994-10-14 | 1995-10-09 | Western Atlas Int Inc | Method for analyzing radiation and apparatus for carrying out the method |
| CN102608649A (en) * | 2012-03-02 | 2012-07-25 | 成都理工大学 | Statistics distributed gamma or X ray energy spectrum unscrambling method |
| CN103135125A (en) * | 2011-12-01 | 2013-06-05 | 中国辐射防护研究院 | Energy spectrum analytical method of X-ray filtering reference radiation |
| JP2016186489A (en) * | 2016-06-01 | 2016-10-27 | コミシリア ア レネルジ アトミック エ オ エナジーズ オルタネティヴズ | Spectrometric analysis method and related device |
| CN106202353A (en) * | 2016-07-06 | 2016-12-07 | 郑州大学 | A kind of visable representation method of time series data |
| TWI600916B (en) * | 2016-04-12 | 2017-10-01 | 行政院原子能委員會核能研究所 | Simple radiation spectrum analyzer |
| CN107229787A (en) * | 2017-05-24 | 2017-10-03 | 南京航空航天大学 | A kind of gamma-ray spectrum analysis method based on approximation coefficient and deep learning |
| CN107526113A (en) * | 2017-08-10 | 2017-12-29 | 中国石油天然气股份有限公司 | Capture the determination method and apparatus of gamma spectra |
-
2018
- 2018-01-25 CN CN201810072826.7A patent/CN110082811B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| NO954005D0 (en) * | 1994-10-14 | 1995-10-09 | Western Atlas Int Inc | Method for analyzing radiation and apparatus for carrying out the method |
| CN103135125A (en) * | 2011-12-01 | 2013-06-05 | 中国辐射防护研究院 | Energy spectrum analytical method of X-ray filtering reference radiation |
| CN102608649A (en) * | 2012-03-02 | 2012-07-25 | 成都理工大学 | Statistics distributed gamma or X ray energy spectrum unscrambling method |
| TWI600916B (en) * | 2016-04-12 | 2017-10-01 | 行政院原子能委員會核能研究所 | Simple radiation spectrum analyzer |
| JP2016186489A (en) * | 2016-06-01 | 2016-10-27 | コミシリア ア レネルジ アトミック エ オ エナジーズ オルタネティヴズ | Spectrometric analysis method and related device |
| CN106202353A (en) * | 2016-07-06 | 2016-12-07 | 郑州大学 | A kind of visable representation method of time series data |
| CN107229787A (en) * | 2017-05-24 | 2017-10-03 | 南京航空航天大学 | A kind of gamma-ray spectrum analysis method based on approximation coefficient and deep learning |
| CN107526113A (en) * | 2017-08-10 | 2017-12-29 | 中国石油天然气股份有限公司 | Capture the determination method and apparatus of gamma spectra |
Non-Patent Citations (4)
| Title |
|---|
| Information retrieval methods for high resolution γ-ray spectra;WU Hexi et al.;《Nuclear Science and Techniques》;20121231;全文 * |
| Reduing noise in gamma-ray spectrometry using spectral component analysis;Minty，B et al.;《EXPLORATION GEOPHYSICS》;20020101;全文 * |
| 低空气比释动能率和窄谱系列过滤X射线参考辐射能谱测量和分析;牛强 等;《辐射防护》;20121031;全文 * |
| 地层自然伽马射线能量峰的识别;王卫华;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20080315;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110082811A (en) | 2019-08-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kaastra et al. | Optimal binning of X-ray spectra and response matrix design | |
| Brenton et al. | Accurate mass measurement: terminology and treatment of data | |
| Aprile et al. | Response of the XENON100 dark matter detector to nuclear recoils | |
| Vidmar et al. | Calculation of true coincidence summing corrections for extended sources with EFFTRAN | |
| CN103913764B (en) | A kind of NaI based on Gaussian response matrix (TI) scintillation detector gamma spectrum high-resolution inversion analysis system and method | |
| CN101819168B (en) | Spectral analysis method | |
| CN103853929A (en) | Low-resolution gamma energy spectrum inversion analysis process and method based on Monte Carlo response matrix | |
| Zhaunerchyk et al. | Theory and simulations of covariance mapping in multiple dimensions for data analysis in high-event-rate experiments | |
| CN110082811B (en) | Energy track counting reassignment process and method for gamma spectrum analysis | |
| CN109471999B (en) | Gamma radiation field data correction calculation method and system for non-uniform source item distribution | |
| CN110515115B (en) | Evaluation method for solving peak area of gamma spectrum | |
| WO2018066587A1 (en) | System and method for accurately quantifying composition of a target sample | |
| CN112906193A (en) | Gauss widening and spectrum-resolving method for gamma response function of detector | |
| Guibert et al. | Non-linear approach of TL response to dose: polynomial approximation | |
| Arinc et al. | Application of the least-squares method to the analysis of XRF spectral intensities from atmospheric particulates collected on filters | |
| CN110121663B (en) | Radioactivity measuring method and radioactivity measuring system | |
| CN110082812B (en) | Improved energy track counting redistribution method for gamma spectrum analysis | |
| JP6882098B2 (en) | Control devices, measurement systems, control methods, and programs | |
| Tsirigotis et al. | Reconstruction efficiency and discovery potential of a Mediterranean neutrino telescope: A simulation study using the Hellenic Open University Reconstruction and Simulation (HOURS) package | |
| Roque | X-ray imaging using 100 µm thick Gas Electron Multipliers operating in Kr-CO2 mixtures | |
| Haefner et al. | Demonstration of Event Position Reconstruction based on Diffusion in the NEXT-White Detector | |
| CN110121664B (en) | Radioactivity measuring method and radioactivity measuring system using data expansion | |
| Lin et al. | Sinc Nyström Method for Singularly Perturbed Love's Integral Equation | |
| Klusoň et al. | Study of errors using scintillation spectra unfolding for the calculation of photon field dosimetric characteristics | |
| Feng et al. | An Integrated Single-Step Gamma Spectrum Analysis Method Based on Bregmanized Variational Deconvolution |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |