CN102621170B - Method for automatically determining measurement time in detection of energy spectrometer - Google Patents
Method for automatically determining measurement time in detection of energy spectrometer Download PDFInfo
- Publication number
- CN102621170B CN102621170B CN 201210061557 CN201210061557A CN102621170B CN 102621170 B CN102621170 B CN 102621170B CN 201210061557 CN201210061557 CN 201210061557 CN 201210061557 A CN201210061557 A CN 201210061557A CN 102621170 B CN102621170 B CN 102621170B
- Authority
- CN
- China
- Prior art keywords
- spectral line
- energy spectrometer
- stability
- measuring time
- data
- 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
Abstract
The invention discloses a method for automatically determining measurement time in the detection of an energy spectrometer. The method is characterized by comprising the following steps of: (1) tracking and monitoring data collected by the energy spectrometer in real time for a sample to be measured and caching the data; (2) carrying out data pretreatment on the spectral line data tracked and cached in the step (1); (3) acquiring a stability parameter according to the spectral line data tracked currently and the spectral line data treated in the step (2); (4) carrying out stability judgment, automatically judging whether the collection of a current spectral line meets a stability condition or not; and automatically stopping the measurement sampling of the sample if the condition is met, otherwise continuing the measurement sampling and tracking the spectral line data if the condition is not met. By using the method, the measurement time can be automatically determined in the process that the sample is measured by the energy spectrometer; the measurement efficiency of the sample is effectively improved; the influence of environmental change on a measurement result is effectively reduced; and meanwhile, the uncertainty of the measurement time judged manually is reduced.
Description
Technical field
The present invention relates to a kind of method of automatically determining Measuring Time in energy spectrometer is surveyed.
Background technology
High energy particle is surveyed the energy spectrometer, as the NaI energy spectrometer, is widely adopted in radioactivity of building materials detection and soil radioactivity detection, and becomes one of thermometrically equipment of GB requirement.The NaI energy spectrometer is general adopts abundant vitriol chamber as shielding, and NaI scintillation crystal coupling photomultiplier consists of gamma-ray detector, and the high-energy photon of sample radiation is surveyed and energy discrimination.Because the NaI energy spectrometer has good environment alpha ray shield condition, the NaI scintillation crystal possesses again higher detection efficiency and energy resolution preferably, thus the NaI detector can measure comparatively accurately the sample of low-activity (radium-226 in sample (
226Ra), thorium-232 (
232Th), kalium-40 (
40When K) the specific radioactivity sum surpassed 37Bq/kg, when spreading factor was 1, uncertainty of measurement was not more than 20%).Because sample activity is very low, in order to obtain enough detection counting, need long Measuring Time: it is generally acknowledged that Measuring Time is longer, the statistical counting that obtains is more, spectral line is carried out the statistical error that the spectrum unscrambling analysis brings less.Yet Measuring Time is long, has a strong impact on the one hand the detection efficiency of sample, also can be subjected to the impact of the factors such as environment temperature, stability of instrument simultaneously, thereby affect measurement result.
In present blind sample detects, measure the experience that often depends on the energy spectrometer user of choosing of duration, this can increase the uncertainty of measurement result.
And the needed Measuring Time of the sample of different activity is also different, the sample of low-activity is with respect to the sample of high activity, often need measurement sampling time of more growing in measurement, yet because sample activity in blind sample detection is unknown, the experimenter tends to adopt the long as far as possible time to sample, and can waste unnecessary Measuring Time like this, affects the detection efficiency of sample, simultaneously, also can increase the probability of the unfavorable factor introducing measurement results such as environmental change.And, if find that when sample measurement is completed Measuring Time is inadequate, generally can require this sample to remeasure, efficient is measured in impact more.
Summary of the invention
Technical matters to be solved by this invention is that a kind of method of automatically determining Measuring Time in energy spectrometer is surveyed is provided for the deficiency of above-mentioned existence.The method can not only effectively improve the measurement efficient of sample, reduces environmental change to the impact of measurement result, and can reduce the artificial uncertainty of judging Measuring Time.。
The present invention addresses the above problem the technical scheme that adopts to be:
A kind of method of automatically determining Measuring Time in energy spectrometer is surveyed is characterized in that comprising the following steps:
A: in real time sample is carried out the tracking and monitoring energy spectrometer and measure spectral line data and the buffer memory that gathers;
B: to the spectral line data of trace cache in steps A, carry out data pre-processing;
C: obtain stability parameter according to the spectral line data of processing in the spectral line data of current tracking and step B;
D: carry out judgement of stability, whether the collection of the current spectral line of automatic decision satisfies stable condition; If satisfy condition, automatically stop the measurement collection of sample, if do not satisfy condition, continue to measure and sample and follow the tracks of spectral line data.
Press such scheme, adopt the energy spectrometer acquisition module to carry out the spectral line collection to sample, adopt the energy spectrometer control module that spectral line data is followed the tracks of detection and buffer memory, data pre-processing, obtained stability parameter and spectral line data is carried out judgement of stability.
Press such scheme, the method that follow the tracks of to detect image data in steps A is for constant duration or obtain the spectral line data of current collection according to the certain time interval cycle, and the mode of opening up special data space or setting up file is carried out buffer memory with this image data.
Press such scheme, described constant duration is 20s-800s.
Press such scheme, the data pre-processing described in step B comprises filtering processing and/or normalized.
Press such scheme, described filtering is treated to location, the road length according to energy spectrometer, and image data is carried out smothing filtering, and filter length is at 0 to 300 point.
Press such scheme, described normalized is time normalization processing and/or sample activity normalized; Described time normalization is treated to the collection spectral line data divided by current Measuring Time, obtains the image data in the unit interval; Described is with gathering the gross-count rate of spectral line data divided by spectral line data, obtaining the collection spectral line data in unit activity to the activity normalized.
Press such scheme, the gross-count rate of described spectral line data is the total counting of spectral line divided by the Measuring Time of track record data for the first time.
Press such scheme, the method for obtaining stability parameter described in step C comprises gathering dynamically carrying out stability calculating and/or the situation of change of sampling spectral line resolving spectra being carried out stability calculate of spectral line shape.
Press such scheme, described dynamic change to spectral line shape carry out stability be calculated as current spectrum wire shaped and front collection spectral line shape of preserving once or several times poor, with the curve of difference dynamic change as stability parameter.
Press such scheme, then described difference dynamic changing curve asks the average of the quadratic sum of sequence for all collection spectral line data and last tracking spectral line data corresponding subtracting each other successively.
Press such scheme, described situation of change to the spectral line resolving spectra is carried out stability and is calculated as current spectral line and the front spectral line of preserving are once or several times done poor contrast, separate and to obtain that in spectral line, each gathers content and the inside and outside Irradiation index of element after analysis of spectrum, with the curve of every resolving spectra dynamic change as stability curve.
By technique scheme, described solution analysis of spectrum adopts least square-direct matrix.
Press such scheme, the judgement of stability described in step D adopts the method judgement of a subthreshold.
Press such scheme, the judgement of stability described in step D adopts the method judgement of sets of threshold values.
Compared with prior art, the present invention has obtained obvious technique effect:
1, automatic monitoring spectral line calculates and differentiates by stability, and definite sample is reasonably measured acquisition time automatically.
2, under the inconsistent prerequisite of Measuring Time under different activity samples, automatically determine to measure the sampling duration and avoided the long-time time waste that brings of measuring of traditional employing, guarantee to have improved the efficient of sample measurement under the prerequisite of measurement result.
3, automatically determine rationally effectively to measure duration, avoided the interference of the adverse effect factors such as environmental change to measurement result, increased the accuracy of measurement result.
Description of drawings
Fig. 1 is the method flow diagram of one embodiment of the invention.
Fig. 2 is sample measurement spectrum in the embodiment of the present invention.
Fig. 3 is spectral line shape stability curve in the embodiment of the present invention.
Fig. 4 is radium element spectral line resolving spectra stability curve in embodiment in the present invention.
Fig. 5 is thorium element spectral line resolving spectra stability curve in embodiment in the present invention.
Fig. 6 is potassium element spectral line resolving spectra stability curve in embodiment in the present invention.
Embodiment
Below in conjunction with accompanying drawing, embodiments of the present invention is further illustrated.
As shown in Figure 1, the present invention comprises the following steps successively:
A: in real time sample is carried out the tracking and monitoring energy spectrometer and measure spectral line data and the buffer memory that gathers;
B: to the spectral line data of trace cache in steps A, carry out data pre-processing;
C: obtain stability parameter according to the spectral line data of processing in the spectral line data of current tracking and step B;
D: carry out judgement of stability, whether the collection of the current spectral line of automatic decision satisfies stable condition; If satisfy condition, automatically stop the measurement collection of sample, if do not satisfy condition, continue to measure and sample and follow the tracks of spectral line data.
Wherein, adopt the energy spectrometer acquisition module to carry out the spectral line collection to sample, the energy spectrometer control module is followed the tracks of detection and buffer memory, data pre-processing, is obtained stability parameter and spectral line data is carried out judgement of stability spectral line data.
Described acquisition module and control module may reside in an equipment or instrument, perhaps are separated from each other, and connect with agreement and communicate by letter by certain communication interface.Acquisition module can be under the control of control module collection signal, perhaps stop collection signal.And can carry out in control module inside corresponding data and process.Described control module can be flush bonding processor, computing machine or other control and computing equipment.What control module realized Measuring Time determines to be divided into four steps automatically, all completes in control module inside, is followed successively by that spectral line is traced and monitored, buffer memory to gathering, the spectral line data pre-treatment, and spectral line stability is calculated and the spectral line determination of stability.
For to rise when the energy spectrometer acquisition module begins to sample to sample, control module starts timer, at interval of 600 seconds, preserves current spectral line data in the method for tracing and monitoring energy spectrometer measurement image data described in steps A.Store method is the spatial cache P (j) that opens up 4096 length; J=1,2,3 ..., 4096.Store the numerical value in 4096 roads in current spectral line.
As shown in Figure 2, in figure, four curves are respectively the spectral line of the 7th secondary tracking, the 14th secondary tracking, the 21st secondary tracking and the 28th secondary tracking.
Spectral line pre-treatment in step B has comprised spectral line filtering and time normalization and three processing modes of sample activity normalization.Spectral line data P (j) formula below consisting of carries out filtering.Filtering is treated to location, the road length according to energy spectrometer, and sampled data is carried out smothing filtering, and filter length is at 0 to 300 point.Obtain filtered data P
SM(j); J=1,2,3 ..., 4096.
The measurement sampling needle of spectral line is to P
SM(j) spectral line is carried out time normalization.Normalized computing formula is as follows:
P
ST(j)＝P
SM(j)/t(i)；j＝1，2，3，...，4096
It is the i time surveillance and tracking that i represents current, and the time is t (i)=600 * i, and unit is second.
Obtain the spectral line P of time normalization
ST(j).
For P
ST(j) spectral line is carried out activity normalization, wherein activity adopts gross-count rate Cts (1) expression of spectral line when following the tracks of for the first time.
Cts (1) calculates with following formula:
The normalized formula of activity is as follows:
P
SA(j)＝P
ST(j)/Cts(1)；j＝1，2，3，...，4096
Obtain the normalized spectral line P of activity
SA(j);
Get stable calculating in described control module and comprise two computing method, be respectively the calculating of spectral line shape stability and the calculating of spectral line resolving spectra stability.
When adopting spectral line shape stability to calculate, need current tracking and normalized spectral line P
SA(j) and last the tracking and spectral line P that normalization obtains
SA_Last(j), adopt following formula computational stability parameter:
Q is the stability parameter of obtaining.
Fig. 3 calculates for each the spectral line that gained spectral line shape stability parameter obtains, until the 10th time, stability parameter satisfies condition, and stops sampled measurements.
When adopting the stability calculating of spectrum unscrambling result of calculation, need current tracking and filtering spectral line P afterwards
SM(j) and last the tracking and the spectral line P of filtering
SM_Last(j): respectively these two spectral lines are separated analysis of spectrum, spectrum unscrambling adopts the algorithm of least square-inverse matrix to obtain radium in spectral line, thorium, three kinds of nucleic of potassium to divide other content, be expressed as C
Ra, C
ThAnd C
KFollow the tracks of in conjunction with last the C that spectral line obtains
Ra_LastC
Th_LastAnd C
L_LastBelow adopting, formula carries out stability calculating:
Q
R＝|C
Ra-C
Ra_Last|
Q
T＝|C
Th-C
Th_Last|
Q
K＝|C
K-C
K_Last|
Q
R, Q
TAnd Q
KBe the stability parameter of calculating.
Described stability criterion adopts the way of a subthreshold, in the spectral line shape determination of stability, we adopt threshold value Threshold＜0.04 as stable condition, in spectral line resolving spectra stability criterion, we adopt sets of threshold values, satisfy simultaneously following three conditions:
Fig. 4, Fig. 5, Fig. 6 are the spectral line that the stability parameter of each resolving spectra obtains radium, thorium, three kinds of nucleic of potassium, until the 15th time, stability condition satisfies, and stops sampled measurements.
With current stability parameter and threshold ratio, if stability parameter satisfies the threshold condition of setting, namely stop measuring sampling.If do not satisfy, continue to get back to the first step, the surveillance and tracking spectral line.
Claims (15)
1. automatically determine the method for Measuring Time for one kind in the energy spectrometer detection, it is characterized in that comprising the following steps:
A: in real time sample is carried out the tracking and monitoring energy spectrometer and measure spectral line data and the buffer memory that gathers;
B: to the spectral line data of trace cache in steps A, carry out data pre-processing;
C: obtain stability parameter according to the spectral line data of processing in the spectral line data of current tracking and step B;
D: carry out judgement of stability, whether the collection of the current spectral line of automatic decision satisfies stable condition; If satisfy condition, automatically stop the measurement collection of sample, if do not satisfy condition, continue to measure and sample and follow the tracks of spectral line data.
2. the method for automatically determining Measuring Time in energy spectrometer is surveyed according to claim 1, it is characterized in that adopting the energy spectrometer acquisition module to carry out the spectral line collection to sample, adopt the energy spectrometer control module that spectral line data is followed the tracks of detection and buffer memory, data pre-processing, obtained stability parameter and spectral line data is carried out judgement of stability.
3. the method for automatically determining Measuring Time in energy spectrometer is surveyed according to claim 1 and 2, it is characterized in that in steps A following the tracks of the method that detects image data for obtaining the spectral line data of current collection with the constant duration cycle, and the mode of opening up special data space or setting up file is carried out buffer memory with this image data.
4. the method for automatically determining Measuring Time in energy spectrometer is surveyed according to claim 3, is characterized in that described constant duration is 20s-800s.
5. the method for automatically determining Measuring Time in energy spectrometer is surveyed according to claim 1 and 2, is characterized in that the data pre-processing described in step B comprises filtering processing and/or normalized.
6. the method for automatically determining Measuring Time in energy spectrometer is surveyed according to claim 5, is characterized in that described filtering is treated to location, the road length according to energy spectrometer, and image data is carried out smothing filtering, and filter length is at 0 to 300 point.
7. the method for automatically determining Measuring Time in energy spectrometer is surveyed according to claim 5, is characterized in that described normalized is time normalization processing and/or sample activity normalized; Described time normalization is treated to the collection spectral line data divided by current Measuring Time, obtains the image data in the unit interval; Described sample activity normalized obtains the collection spectral line data in unit activity for gathering the gross-count rate of spectral line data divided by spectral line data.
8. the method for automatically determining Measuring Time in energy spectrometer is surveyed according to claim 7, the gross-count rate that it is characterized in that described spectral line data are the total counting of spectral line divided by the Measuring Time of track record data for the first time.
9. the method for automatically determining Measuring Time in energy spectrometer is surveyed according to claim 1 and 2, it is characterized in that the method for obtaining stability parameter described in step C comprises carrying out dynamically to what gather spectral line shape that stability is calculated and/or the situation of change of sampling spectral line resolving spectra is carried out stability calculating.
10. the method for automatically determining Measuring Time in energy spectrometer is surveyed according to claim 9, it is characterized in that described to gather spectral line shape dynamically carry out stability be calculated as the current spectrum wire shaped poor with the collection spectral line shape of front preservation once or several times, with the curve of difference dynamic change as stability parameter.
11. the method for automatically determining Measuring Time in energy spectrometer is surveyed according to claim 10, it is characterized in that described difference dynamic changing curve is that all collection spectral line data and last time are followed the tracks of spectral line data corresponding subtracting each other successively, then asks the average of the quadratic sum of sequence.
12. the method for automatically determining Measuring Time in energy spectrometer is surveyed according to claim 9, it is characterized in that described situation of change to sampling spectral line resolving spectra carries out stability and be calculated as current spectral line and the front spectral line of preserving are once or several times done poor the contrast, separate and to obtain that in spectral line, each gathers content and the inside and outside Irradiation index of element after analysis of spectrum, with the curve of every resolving spectra dynamic change as stability curve.
13. the method for automatically determining Measuring Time in energy spectrometer is surveyed according to claim 12 is characterized in that described solution analysis of spectrum adopts least square-direct matrix.
14. the method for automatically determining Measuring Time in energy spectrometer is surveyed according to claim 10 is characterized in that
Judgement of stability described in step D adopts the method judgement of a subthreshold.
15. the method for automatically determining Measuring Time in energy spectrometer is surveyed according to claim 12 is characterized in that the judgement of stability described in step D adopts the method judgement of sets of threshold values.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201210061557 CN102621170B (en) | 2012-03-09 | 2012-03-09 | Method for automatically determining measurement time in detection of energy spectrometer |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201210061557 CN102621170B (en) | 2012-03-09 | 2012-03-09 | Method for automatically determining measurement time in detection of energy spectrometer |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102621170A CN102621170A (en) | 2012-08-01 |
| CN102621170B true CN102621170B (en) | 2013-11-06 |
Family
ID=46561212
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 201210061557 Active CN102621170B (en) | 2012-03-09 | 2012-03-09 | Method for automatically determining measurement time in detection of energy spectrometer |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102621170B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104239709A (en) * | 2014-09-05 | 2014-12-24 | 西安奥华电子仪器有限责任公司 | Method for determining yield of stratum element logging instrument by virtue of spectrum unfolding |
| CN104239710A (en) * | 2014-09-05 | 2014-12-24 | 西安奥华电子仪器有限责任公司 | Method for confirming yield through spectrum unfolding analysis of stratum element logging instrument |
| CN108457644B (en) * | 2017-12-26 | 2020-04-10 | 中国石油天然气股份有限公司 | Gamma-ray energy spectrum unscrambling method and device for element capture energy spectrum logging |
| CN111551978B (en) * | 2020-05-08 | 2022-09-27 | 中国辐射防护研究院 | Method for evaluating uncertainty of radioactive soil sample sampling process |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3123860B2 (en) * | 1993-07-20 | 2001-01-15 | 日本電子株式会社 | Elemental analyzer using wavelength dispersive spectrometer and energy dispersive spectrometer |
| JP2000030657A (en) * | 1998-07-13 | 2000-01-28 | Anelva Corp | Appearance mass analyzing method |
| JP3790643B2 (en) * | 1999-07-07 | 2006-06-28 | 日本電子株式会社 | Surface analyzer with energy dispersive X-ray detector |
| CN1115571C (en) * | 1999-08-18 | 2003-07-23 | 上海伽玛射线技术有限公司 | Specific radioactivity measuring method and equipment |
| JP4822826B2 (en) * | 2005-12-06 | 2011-11-24 | 日本電子株式会社 | Particle analysis method and apparatus |
| DE112008001646T5 (en) * | 2007-06-20 | 2010-06-10 | Thermo Electron Scientific Instruments LLC, Madison | Processing spectrometer pile-up events |
| CN101806910B (en) * | 2010-05-11 | 2012-07-18 | 西北核技术研究所 | Method for absolutely measuring for radioactive gas nuclide activity |
-
2012
- 2012-03-09 CN CN 201210061557 patent/CN102621170B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN102621170A (en) | 2012-08-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107229787B (en) | Gamma energy spectrum analysis method based on approximation coefficient and deep learning | |
| CN103837558B (en) | Multielement composition and content detection device and detection method in a kind of aqueous solution based on PGNAA technology | |
| CN103424766B (en) | A kind of nuclide rapid identification method based on pattern-recognition | |
| CN105093261B (en) | A kind of neutron, the γ discrimination systems of digitlization neutron spectrometer | |
| CN103605148B (en) | A kind of gamma spectrometry method under high count rate | |
| CN111221030A (en) | Neutron-gamma detector based on physical integration and neutron-gamma online screening method | |
| CN102621170B (en) | Method for automatically determining measurement time in detection of energy spectrometer | |
| CN104360376A (en) | Gamma camera having function of identifying radioactive source, namely nuclide, and nuclide identification method | |
| US20150025852A1 (en) | Signal data processing for radiation detection | |
| CN105425275B (en) | Anti-coincidence output discrimination system | |
| Li et al. | Review of recent gamma spectrum unfolding algorithms and their application | |
| CN103852475A (en) | Multichannel potassium measuring instrument based on gamma ray | |
| CN206020674U (en) | A kind of portable Low background α, β measuring instrument | |
| CN108008440A (en) | low background α, β measuring device | |
| RU2004105034A (en) | ENVIRONMENTAL RADIOACTIVITY MONITOR | |
| CN204789406U (en) | A signal processing system that is used for beta to penetrate line method measuring atmospheric particulates concentration | |
| CN208334663U (en) | low background α, β measuring device | |
| KR101282962B1 (en) | neutron puluse-height spectrum analysis method and the neutron detector system using thereof | |
| CN104111482A (en) | Double-detector X-ray fluorescence logging probe tube and method | |
| CN111239797B (en) | Collector based on radiation particle event and rapid nuclide identification method | |
| CN203285434U (en) | Logging instrument and system | |
| EP3676640B1 (en) | Methods and systems for calibration of particle detectors | |
| CN107422362B (en) | A kind of noise-reduction method of energy spectrum and time spectrum | |
| CN105628715A (en) | Potassium measuring instrument of sodium iodide crystals of raw material potassium ores | |
| CN106405609A (en) | Multiparticle event capturing method and device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |