CN109212583A - A kind of multiple tracks gamma energy spectrometer temperature-compensation method based on fitting of a polynomial - Google Patents
A kind of multiple tracks gamma energy spectrometer temperature-compensation method based on fitting of a polynomial Download PDFInfo
- Publication number
- CN109212583A CN109212583A CN201710545798.1A CN201710545798A CN109212583A CN 109212583 A CN109212583 A CN 109212583A CN 201710545798 A CN201710545798 A CN 201710545798A CN 109212583 A CN109212583 A CN 109212583A
- Authority
- CN
- China
- Prior art keywords
- temperature
- multiple tracks
- fitting
- energy spectrometer
- polynomial
- 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.)
- Granted
Links
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
Abstract
The invention belongs to radioactivity survey fields, and in particular to a kind of multiple tracks gamma energy spectrometer temperature-compensation method based on fitting of a polynomial.Specifically includes the following steps: setting measurement parameter Step 1: the data based on multiple tracks gamma energy spectrometer acquire；Step 2: establishing temperature correction model, multinomial correction coefficient is determined；Step 3: being based on fitting of a polynomial, chooses mineral samplers and be corrected using the multinomial coefficient in measurement parameter and step 2 in step 1.The present invention establishes temperature correction curve using Digital Signal Processing combination polynomial fitting, realize the real time correction of measurement spectrogram, the related fieldss services such as the qualitative and quantitative analysis ability for improving spectrogram is preferably exploration prospecting, oil gas generaI investigation, regional geologic reconnaissance.
Description
Technical field
The invention belongs to radioactivity survey fields, and in particular to a kind of multiple tracks gamma energy spectrometer temperature based on fitting of a polynomial
Compensation method.
Background technique
Nuclear geophysical prospecting gamma spectrum measuring technique is usually used in reconnoitring radioactive occurrence, and be successfully used for Solid Mineral at present looks for mine
Exploration, oil-gas field reconnaissance, regional geologic reconnaissance, hydrogeology and engineering geology etc..With to low energy γ in natural radiation field
The further investigation of spectroscopy and its application, and adapt to the need of environmental radiation investigation, building materials and building decoration materials radiological measuring
It wants, it is desirable that gamma spectrum system should be able to adapt to measure difference caused by environment temperature.At present for the side of gamma spectrum temperature-compensating
Method mainly has temperature-compensation circuit realization.The system is mainly by main amplifier, peak value of pulse retainer, multichannel pulse amplitude analysis
Device, slide rule, multiple tracks data buffer and the composition such as embedded micro energy lose computer and its peripheral hardware.Using " drawing outside every interior " heat transfer
Technology blocks the variation of external temperature field which by thermal insulation layer, can only slowly, equably pass to the inside of probe, and in popping one's head in
The variation in portion temperature field can be sent to rapidly other positions again, guarantee that the consistency of probe interior temperature distribution realizes temperature with this
Degree compensation more can.
Summary of the invention
The multiple tracks gamma energy spectrometer temperature-compensating based on fitting of a polynomial that the technical problem to be solved in the present invention is to provide a kind of
Method establishes temperature correction curve using Digital Signal Processing combination polynomial fitting, realizes the reality of measurement spectrogram
Shi Jiaozheng improves the qualitative and quantitative analysis ability of spectrogram, is preferably exploration prospecting, oil gas generaI investigation, the phases such as regional geologic reconnaissance
Close industry service.
In order to solve the above technical problems, a kind of multiple tracks gamma energy spectrometer temperature-compensating side based on fitting of a polynomial of the present invention
Method, specifically includes the following steps:
Step 1: the data based on multiple tracks gamma energy spectrometer acquire, measurement parameter is set；
Step 2: establishing temperature correction model, multinomial correction coefficient is determined；
Step 3: being based on fitting of a polynomial, mineral samplers are chosen, using more in measurement parameter and step 2 in step 1
Binomial coefficient is corrected.
In the step 2, relevant voltage threshold value and sampling time are set, in k different temperatures T1,T2,T3,...,TkItem
Measurement standard source under part137The gamma spectrum figure of Cs, under corresponding temperature137The peak position of Cs is P1,P2,P3,...,Pk；
Polynomial of degree n correction coefficient a0、a1、a2、…、anIt is determined by following formula:
Wherein ()-1Representing matrix asks quasi-, ()TRepresenting matrix transposition.
In the step 3, mineral samplers are specially chosen, using multinomial in measurement parameter and step 2 in step 1
Formula coefficient chooses temperature TmFor normal temperature (1≤m≤k), TjAt a temperature of137Cs peak position is Pj, will spectral data school under this condition
Just arrive TmAt a temperature of:
Therefore
Wherein PjTo correct preceding peak position,For peak position after correction.
In the step 1, the measurement parameter of multiple tracks gamma energy spectrometer is sample frequency: 8M HZ.
In the step 1, the measurement parameter of multiple tracks gamma energy spectrometer is time of measuring: 10 hours.
In the step 1, the measurement parameter of multiple tracks gamma energy spectrometer is voltage pulse threshold value: 2050.
In the step 1, the measurement parameter of multiple tracks gamma energy spectrometer is energy calibration coefficient: slope 2.91, intercept :-
5971。
The beneficial technical effect of the present invention lies in: this method is based on Digital Signal Processing, utilizes the temperature of formulation, peak
Bit model establishes calibration curve, is compensated to temperature, it can be achieved that measuring the real time correction of spectrogram.Temperature based on fitting of a polynomial
It spends backoff algorithm and handles spectral data, realize quick, the precise calibration of gamma spectrum figure, surveyed for radioactivity survey, field geology
The related fieldss such as spy, environmental evaluation provide powerful technique guarantee.
Detailed description of the invention
Fig. 1 is one standard source of embodiment137Multiple tracks gamma spectrum figure under the conditions of Cs T=25 °；
Fig. 2 is multiple tracks gamma spectrum figure under the conditions of T=-10 ° of embodiment two；
Fig. 3 is multiple tracks gamma spectrum figure under the conditions of T=1 ° of embodiment two；
Fig. 4 is multiple tracks gamma spectrum figure under the conditions of T=25 ° of embodiment two；
Fig. 5 is multiple tracks gamma spectrum figure under the conditions of T=50 ° of embodiment two.
Specific embodiment
Invention is further described in detail with reference to the accompanying drawings and examples.
Embodiment one
A kind of multiple tracks gamma energy spectrometer temperature-compensation method based on fitting of a polynomial of the present invention, specifically includes the following steps:
Step 1: the data based on multiple tracks gamma energy spectrometer acquire, measurement parameter is set；
NaI crystal counter output signal pulses use Sing plus waveform Spectrum stripping techniques in analysis in pulse amplitude analysis part
The data processing method for replacing traditional peak value sampling to keep completes the complicated hardware capability of tradition with software technology, to improve
All technical, reliability, flexibility and the operability of instrument are carried out one by one by the impulse amplitude generated to gamma-rays
The accurate measurement to radioactive sample content is realized in parsing；
This instrument is aided with pulse former, pulse amplifying circuit, high pressure conversion electricity by combined type NaI crystal counter
The part such as road, A/D sample conversion circuit, data collection and analysis processing software forms.As shown in Fig. 1；
Using the above-mentioned multiple tracks gamma energy spectrometer system acquisition data built；In order to obtain stable gamma spectrum, relevant parameter
It adjusts as follows:
Sample frequency: 8M HZ
Time of measuring: 10 hours
Voltage pulse threshold value: 2050
Energy calibration coefficient: slope 2.91, intercept: -5971；
Step 2: establishing temperature correction model, multinomial correction coefficient is determined, using the measurement parameter in step 1, if
Relevant voltage threshold value and sampling time are set, in k different temperatures T1,T2,T3,...,TkUnder the conditions of measurement standard source137The γ of Cs
Energy spectrum diagram, under corresponding temperature137The peak position of Cs is P1,P2,P3,...,Pk；Spectrogram such as Fig. 1 institute under the conditions of one of temperature
Show；
Temperature T is fitted using polynomial of degree n and peak position P relationship is as follows:
Wherein a0、a1、a2、…、anFor polynomial of degree n correction coefficient.The matrix form of above-mentioned formula is as follows:
Under above-mentioned k temperature experiment condition, following matrix equation is obtained using polynomial of degree n fittingTherefore polynomial of degree n correction coefficient a0、a1、a2、…、anIt is determined by following formula:
Wherein ()-1Representing matrix asks quasi-, ()TRepresenting matrix transposition.
Step 3: being based on fitting of a polynomial, mineral samplers are chosen, using more in measurement parameter and step 2 in step 1
Binomial coefficient chooses temperature TmFor normal temperature (1≤m≤k), TjAt a temperature of137Cs peak position is Pj, will spectral data under this condition
It is corrected to TmAt a temperature of:
Therefore
Wherein PjTo correct preceding peak position,For peak position after correction.
Embodiment two
Step 1, it obtains under normal temperature137The gamma spectrum figure of Cs
Relevant voltage threshold value and sampling time are set, in 4 different temperatures: T1=-10 °, T2=1 °, T3=25 °, T4=
50 lower measurement standard sources137The gamma spectrum figure of Cs is as shown in Fig. 2-Fig. 5.
Step 2, calibration model is established, correction coefficient is obtained
Using spectral data under the conditions of the aforementioned four temperature of 3 order polynomials fitting, following matrix equation is established:
Wherein P1=535.44, P2=582.48, P3=673.62, P4=797.10；T1=-10 °, T2=1 °, T3=
25°、T4=50 °.Therefore correction coefficient are as follows:
Step 3, based on the temperature-compensating of fitting of a polynomial
Identical sampling condition is set, is acquired under different temperatures137Cs peak position is as shown in table 1.With room temperature T in experiment3=25 °
For standard (peak position: 673.62Kev), spectral data in table 1 is corrected at this temperature.Utilize the correction coefficient in step 2
And the peak position such as table 2 after following formula correction:
Wherein Tm=25, a0、a1、a2、a3Value is as shown in step 2.
It can be seen that the multiple tracks gamma spectrum figure for acquiring under condition of different temperatures from table 2, by quasi- based on multinomial
After the temperature compensation algorithm correction of conjunction,137Cs peak position is all corrected near 673Kev substantially, average deviation 0.78Kev, substantially
Meet practical application request.
Under 1 condition of different temperatures of table137Cs peak position: temperature (degree Celsius) peak position (Kev)
| Temperature | -5 | 5 | 10 | 15 | 30 | 35 | 40 | 45 |
| Peak position | 556.67 | 598.01 | 618.37 | 634.63 | 694.56 | 716.14 | 741.72 | 766.27 |
After 2 temperature-compensating of table correction137Cs peak position
Claims (10)
1. a kind of multiple tracks gamma energy spectrometer temperature-compensation method based on fitting of a polynomial, it is characterised in that: specifically include following step
It is rapid:
Step 1: the data based on multiple tracks gamma energy spectrometer acquire, measurement parameter is set；
Step 2: establishing temperature correction model, multinomial correction coefficient is determined；
Step 3: being based on fitting of a polynomial, mineral samplers are chosen, using the multinomial in measurement parameter and step 2 in step 1
Coefficient is corrected.
2. a kind of multiple tracks gamma energy spectrometer temperature-compensation method based on fitting of a polynomial according to claim 1, feature
It is: in the step 2, relevant voltage threshold value and sampling time is set, in k different temperatures T1,T2,T3,...,TkCondition
Lower measurement standard source137The gamma spectrum figure of Cs, under corresponding temperature137The peak position of Cs is P1,P2,P3,...,Pk；
Polynomial of degree n correction coefficient a0、a1、a2、…、anIt is determined by following formula:
Wherein ()-1Representing matrix asks quasi-, ()TRepresenting matrix transposition.
3. a kind of multiple tracks gamma energy spectrometer temperature-compensation method based on fitting of a polynomial according to claim 2, feature
It is: in the step 3, mineral samplers is specially chosen, using the system of polynomials in measurement parameter and step 2 in step 1
Number chooses temperature TmFor normal temperature (1≤m≤k), TjAt a temperature of137Cs peak position is Pj, spectral data will be corrected under this condition
TmAt a temperature of:
Therefore
Wherein PjTo correct preceding peak position,For peak position after correction.
4. a kind of multiple tracks gamma energy spectrometer temperature-compensation method based on fitting of a polynomial according to claim 3, feature
Be: in the step 1, the measurement parameter of multiple tracks gamma energy spectrometer is sample frequency: 8M HZ.
5. a kind of multiple tracks gamma energy spectrometer temperature-compensation method based on fitting of a polynomial according to claim 4, feature
Be: in the step 1, the measurement parameter of multiple tracks gamma energy spectrometer is time of measuring: 10 hours.
6. a kind of multiple tracks gamma energy spectrometer temperature-compensation method based on fitting of a polynomial according to claim 5, feature
Be: in the step 1, the measurement parameter of multiple tracks gamma energy spectrometer is voltage pulse threshold value: 2050.
7. a kind of multiple tracks gamma energy spectrometer temperature-compensation method based on fitting of a polynomial according to claim 6, feature
Be: in the step 1, the measurement parameter of multiple tracks gamma energy spectrometer is energy calibration coefficient: slope 2.91, intercept: -5971.
8. a kind of multiple tracks gamma energy spectrometer temperature-compensation method based on fitting of a polynomial according to claim 1, feature
Be: in the step 1, the measurement parameter of multiple tracks gamma energy spectrometer is sample frequency: 8M HZ, time of measuring: 10 hours, electricity
Pressure pulse threshold: 2050, energy calibration coefficient: slope 2.91, intercept: -5971.
9. a kind of multiple tracks gamma energy spectrometer temperature-compensation method based on fitting of a polynomial according to claim 8, feature
It is: in the step 2, relevant voltage threshold value and sampling time is set, in k different temperatures T1,T2,T3,...,TkCondition
Lower measurement standard source137The gamma spectrum figure of Cs, under corresponding temperature137The peak position of Cs is P1,P2,P3,...,Pk；
Polynomial of degree n correction coefficient a0、a1、a2、…、anIt is determined by following formula:
Wherein ()-1Representing matrix asks quasi-, ()TRepresenting matrix transposition.
10. a kind of multiple tracks gamma energy spectrometer temperature-compensation method based on fitting of a polynomial according to claim 9, feature
It is: in the step 3, mineral samplers is specially chosen, using the system of polynomials in measurement parameter and step 2 in step 1
Number chooses temperature TmFor normal temperature (1≤m≤k), TjAt a temperature of137Cs peak position is Pj, spectral data will be corrected under this condition
TmAt a temperature of:
Therefore
Wherein PjTo correct preceding peak position,For peak position after correction.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710545798.1A CN109212583B (en) | 2017-07-06 | 2017-07-06 | Multi-channel gamma energy spectrometer temperature compensation method based on polynomial fitting |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710545798.1A CN109212583B (en) | 2017-07-06 | 2017-07-06 | Multi-channel gamma energy spectrometer temperature compensation method based on polynomial fitting |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109212583A true CN109212583A (en) | 2019-01-15 |
| CN109212583B CN109212583B (en) | 2020-10-23 |
Family
ID=64993424
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710545798.1A Active CN109212583B (en) | 2017-07-06 | 2017-07-06 | Multi-channel gamma energy spectrometer temperature compensation method based on polynomial fitting |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109212583B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113729751A (en) * | 2021-10-12 | 2021-12-03 | 上海交通大学 | Detection method and equipment for tracer activity and PET scanner |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1045186A (en) * | 1989-02-22 | 1990-09-05 | 成都地质学院 | Digital spectrum-stabilizing device with natural low-energy peak source for referencial use |
| CN101713829A (en) * | 2009-12-22 | 2010-05-26 | 清华大学 | Rapid spectrum stabilization method and device for ray energy spectrometer |
| CN102175703A (en) * | 2011-01-10 | 2011-09-07 | 长沙开元仪器股份有限公司 | Method for correcting peak shift |
| CN102279408A (en) * | 2011-03-28 | 2011-12-14 | 核工业西南物理研究院 | Processing method of multichannel energy spectrometry and apparatus thereof |
| US20140084149A1 (en) * | 2012-09-27 | 2014-03-27 | Christian Stoller | Use of spectral information to extend temperature range of gamma-ray detector |
| CN104777509A (en) * | 2015-04-29 | 2015-07-15 | 东华理工大学 | Aviation gamma energy disperse spectroscopy spectrum stabilization method based on energy spectrum reconstitution technology |
| US20170139063A1 (en) * | 2015-11-18 | 2017-05-18 | Weatherford Technology Holdings, Llc | Gain Stabilization of Radiation Detectors Via Spectrum Analysis |
-
2017
- 2017-07-06 CN CN201710545798.1A patent/CN109212583B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1045186A (en) * | 1989-02-22 | 1990-09-05 | 成都地质学院 | Digital spectrum-stabilizing device with natural low-energy peak source for referencial use |
| CN101713829A (en) * | 2009-12-22 | 2010-05-26 | 清华大学 | Rapid spectrum stabilization method and device for ray energy spectrometer |
| CN102175703A (en) * | 2011-01-10 | 2011-09-07 | 长沙开元仪器股份有限公司 | Method for correcting peak shift |
| CN102279408A (en) * | 2011-03-28 | 2011-12-14 | 核工业西南物理研究院 | Processing method of multichannel energy spectrometry and apparatus thereof |
| US20140084149A1 (en) * | 2012-09-27 | 2014-03-27 | Christian Stoller | Use of spectral information to extend temperature range of gamma-ray detector |
| CN104777509A (en) * | 2015-04-29 | 2015-07-15 | 东华理工大学 | Aviation gamma energy disperse spectroscopy spectrum stabilization method based on energy spectrum reconstitution technology |
| US20170139063A1 (en) * | 2015-11-18 | 2017-05-18 | Weatherford Technology Holdings, Llc | Gain Stabilization of Radiation Detectors Via Spectrum Analysis |
| US9702990B2 (en) * | 2015-11-18 | 2017-07-11 | Weatherford Technology Holdings, Llc | Gain stabilization of radiation detectors via spectrum analysis |
Non-Patent Citations (2)
| Title |
|---|
| 杨焕章: "航空γ能谱探头温度效应及校正技术研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
| 田酉牧: "光纤陀螺温度漂移与补偿方法的研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113729751A (en) * | 2021-10-12 | 2021-12-03 | 上海交通大学 | Detection method and equipment for tracer activity and PET scanner |
| CN113729751B (en) * | 2021-10-12 | 2024-01-30 | 上海交通大学 | Method and equipment for detecting activity of tracer and PET scanner |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109212583B (en) | 2020-10-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| McLean et al. | An algorithm for U‐Pb isotope dilution data reduction and uncertainty propagation | |
| Cummer et al. | The lightning‐TGF relationship on microsecond timescales | |
| Casanovas et al. | Temperature peak-shift correction methods for NaI (Tl) and LaBr3 (Ce) gamma-ray spectrum stabilisation | |
| Fletcher et al. | Matrix effects and calibration limitations in ion probe U–Pb and Th–Pb dating of monazite | |
| Cano-Ott et al. | Pulse pileup correction of large NaI (Tl) total absorption spectra using the true pulse shape | |
| Seljak | Bias, redshift space distortions and primordial nongaussianity of nonlinear transformations: application to Ly-α forest | |
| Cummer et al. | Three years of lightning impulse charge moment change measurements in the United States | |
| Smith et al. | Note on the experimental determination of the relative fast-neutron sensitivity of a hydrogenous scintillator | |
| CN104777509B (en) | Airborne Gamma-ray spectrometry instrument spectrum stabilization method based on power spectrum reconfiguration technique | |
| Abbaszadeh et al. | Characterization of a sub-assembly of 3D position sensitive cadmium zinc telluride detectors and electronics from a sub-millimeter resolution PET system | |
| Peterson et al. | Photoelectrons and solar ionizing radiation at Mars: Predictions versus MAVEN observations | |
| Brodrick et al. | X‐ray magnitude of the 4 November 2003 solar flare inferred from the ionospheric attenuation of the galactic radio background | |
| Pu et al. | VHF radio spectrum of a positive leader and implications for electric fields | |
| Bhardwaj et al. | On the evening time exosphere of Mars: Result from MENCA aboard Mars Orbiter Mission | |
| Hellesen et al. | Impact of digitization for timing and pulse shape analysis of scintillator detector signals | |
| Reeder et al. | Gain calibration of a β/γ coincidence spectrometer for automated radioxenon analysis | |
| Dickson | Recent advances in aerial gamma-ray surveying | |
| CN109212583A (en) | A kind of multiple tracks gamma energy spectrometer temperature-compensation method based on fitting of a polynomial | |
| Romer et al. | A low-power analog front-end amplifier for SiPM based radiation detectors | |
| Kulisek et al. | Real-time airborne gamma-ray background estimation using NASVD with MLE and radiation transport for calibration | |
| VanBommel et al. | Seasonal atmospheric argon variability measured in the equatorial region of Mars by the Mars exploration rover alpha particle X‐ray spectrometers: Evidence for an annual argon‐enriched front | |
| Arabshahi et al. | The energy spectrum of X‐rays from rocket‐triggered lightning | |
| Gardner et al. | Use of an iterative convolution approach for qualitative and quantitative peak analysis in low resolution gamma-ray spectra | |
| Sabbatucci et al. | First principles pulse pile-up balance equation and fast deterministic solution | |
| Herndon | Measurement analysis in uncertainty space |
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 | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Li Boping Inventor after: Guo Dongfa Inventor after: Xia Chenguang Inventor after: Huang Qiuhong Inventor after: Dong Chen Inventor after: Hu Yong Inventor before: Hu Yong Inventor before: Guo Dongfa Inventor before: Li Boping Inventor before: Huang Qiuhong Inventor before: Dong Chen |