CN102175703A - Method for correcting peak shift - Google Patents
Method for correcting peak shift Download PDFInfo
- Publication number
- CN102175703A CN102175703A CN 201110004270 CN201110004270A CN102175703A CN 102175703 A CN102175703 A CN 102175703A CN 201110004270 CN201110004270 CN 201110004270 CN 201110004270 A CN201110004270 A CN 201110004270A CN 102175703 A CN102175703 A CN 102175703A
- Authority
- CN
- China
- Prior art keywords
- peak
- road number
- energy
- spectrum
- carbon
- 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
Images
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
The invention discloses a method for correcting peak shift, wherein after the initial functional relation between the energy and the number of channels in a test spectrum is determined via a carbon peak, a hydrogen peak and a leak peak among characteristic peaks, multiple secondary trinomial fittings are performed by utilizing other peaks among the characteristic peaks, thus the final functional relation between the energy and the number of channel in the test spectrum is obtained; and then the characteristic peaks in the test spectrum are enabled to be similar to the characteristic peaks in a standard spectrum to the greatest extent via matching with the standard spectrum, thus the accuracy of the correction is improved.
Description
Technical field
The present invention relates to the pulse neutron coal analysis technical field, particularly relate to a kind of peak drift correction method.
Background technology
The pulse neutron coal on-line analysis is a kind of application of prompt fission gammas's neutron activation analysis technology.It is with the traditional coal property test mode and the difference of method in the past, do not need to coal to be measured sample, sample preparation and sampling operation, as long as checkout equipment is installed on the belt of carrying coal, in to the course of conveying of coal, just finished the analytical work of coal to be measured simultaneously.Because this method has been saved sampling, sample preparation and sampling process,, be a kind of means of coal property test efficiently so can save chemical examination time and chemical examination cost significantly.This method is utilized neutron bombardment coal stream, when colliding nucleic in the coal stream, neutron will react, discharge gamma rays, when not simultaneously by the nucleic that bombarded, the gamma rays that discharges also can be variant on energy, can judge in the stream that produces coal have which nucleic on earth by the difference of these energy.In addition, for same nucleic, the gamma rays that radiate its content in coal stream and it be counted as proportionate relationship.The counting that therefore only need write down the corresponding energy of nucleic just can be known the content of specific species in the coal stream.
The energy of gamma rays and quantity information can use detector and multichannel analyzer to collect.Detector can detect gamma rays, and output is by the gamma rays power spectrum of energy and counting sign.But when temperature and environment changed, detector deviation can occur in detection process, and the energy in the gamma rays power spectrum that causes collecting drifts about, i.e. peak drift.This has just influenced the precision that the ature of coal nucleic is measured to a great extent.Therefore need revise the peak drift, solve the measuring accuracy problem that the peak drift brings with this.
At present, correction to the peak drift is mainly passed through: will test spectrum according to experience and be divided into a plurality of zones, in each zone, seek characteristic peak respectively, afterwards with standard spectrum in the response characteristic peak compare, difference between the road number that the road number of characteristic peak correspondence in the test spectrum is corresponding with characteristic peak in the standard spectrum is coefficient by way of compensation, and the test spectrum is revised.
In the process that realizes the application, the inventor finds that there are the following problems at least in the prior art: the prerequisite that said method is realized is that the peak drift degree of test spectrum is less, when the peak drift is big, the characteristic peak appointed area that may drift about out, enter other zones adjacent with the appointed area, when seeking characteristic peak according to the method described above, will in this appointed area, carry out the peak-seeking operation, can think other peaks by mistake it is the characteristic peak that will seek this moment, the erroneous judgement at occurrence characteristics peak, revise according to the difference between the road number of this Characteristics of Fault peak correspondence road number corresponding afterwards, reduced the accuracy of peak drift correction with characteristic peak in the standard spectrum.In addition, assemble a fairly large number of peak in the high energy district of test spectrum, each peak-to-peak close together, the very easily erroneous judgement at occurrence characteristics peak this moment reduces the accuracy of peak drift correction.
Summary of the invention
In view of this, the invention provides a kind of peak drift correction method, can improve the accuracy of peak drift correction.
Its concrete scheme is as follows:
A kind of peak drift correction method comprises:
Test spectrum to multichannel analyzer output is carried out smooth treatment, and the test that obtains after smooth is composed;
Described test spectrum after smooth is carried out the peak-seeking operation, determine the actual road number corresponding respectively with carbon peak, hydrogen peak and plumbous peak;
With the peak energy at described carbon peak, hydrogen peak and plumbous peak and actual road number as current known quantity group;
Utilize described current known quantity group to carry out the second degree trinomial expression match, determine the road number of test spectrum and the current functional relation between the energy;
Determine in the characteristic peak other any peak, utilize the peak energy at described peak and the theoretical road number that described current functional relation is determined described peak;
With the theoretical road number at described peak is that the peak-seeking operation is carried out at the center in first preset range, determines the actual road number at described peak, and the peak energy and the actual road number at described peak is added into current known quantity group;
Utilize described current known quantity group to carry out the second degree trinomial expression match, determine the functional relation between number and the energy, and be defined as current functional relation;
Judge whether the number of times that described second degree trinomial expression match is carried out reaches preset times, if then utilize current functional relation to determine the pairing energy of each road number, and described test spectrum and standard spectrum are mated; If not, then turn to the step of determining other any peak in the characteristic peak.
Preferably, in said method, describedly determine that respectively the process of the actual road number corresponding with carbon peak, hydrogen peak and plumbous peak is specially:
In described standard spectrum, determine the standard track number of described carbon peak, hydrogen peak and plumbous peak correspondence successively;
Standard track number with described carbon peak, hydrogen peak and plumbous peak correspondence is the center successively in described test spectrum, searches the peak maximum in second preset range;
With pairing road, the peak maximum that finds out successively number respectively as the actual road number of described carbon peak, hydrogen peak and plumbous peak correspondence.
Preferably, in said method, described first preset range be multichannel analyzer total road number ± 0.5%, described second preset range be described multichannel analyzer total road number ± 2%.
Preferably, in said method, described preset times is 5 times.
Preferably, in said method, described characteristic peak comprises: carbon peak, hydrogen peak, plumbous peak, carbon list are escaped peak, oxygen peak, first iron peak and second iron peak.
Preferably, in said method, least square method is adopted in described second degree trinomial expression match.
As described in above-mentioned scheme, in drift correction method in peak disclosed by the invention, after the initial function relation between energy and the road number in the definite test spectrum in the carbon peak in the characteristic peak, hydrogen peak and plumbous peak, by repeatedly second degree trinomial expression match, obtain the energy of test spectrum and the final funtcional relationship between the road number, by mating, make that characteristic peak and the characteristic peak in the standard spectrum in the test spectrum is approaching as far as possible afterwards, improved the accuracy of revising with standard spectrum.
Description of drawings
In order to be illustrated more clearly in the embodiment of the invention or technical scheme of the prior art, to do to introduce simply to the accompanying drawing of required use in embodiment or the description of the Prior Art below, apparently, accompanying drawing in describing below only is some embodiments of the present invention, for those of ordinary skills, under the prerequisite of not paying creative work, can also obtain other accompanying drawing according to these accompanying drawings.
Fig. 1 is the process flow diagram of a kind of peak drift correction method disclosed by the invention;
Fig. 2 is the process flow diagram of method of the actual road number at definite carbon disclosed by the invention peak, hydrogen peak and plumbous peak.
Embodiment
Drift correction method in peak disclosed by the invention, the utilization least square method is carried out the second degree trinomial expression match, and carry out repeatedly iteration, the gamma-ray spectrometry that multichannel analyzer collects in the paired pulses neutron ature of coal online detection instrument is handled, to solve the problem of easily occurrence characteristics peak erroneous judgement in the prior art, correction accuracy.
Below in conjunction with the accompanying drawing in the embodiment of the invention, the technical scheme in the embodiment of the invention is clearly and completely described, obviously, described embodiment only is the present invention's part embodiment, rather than whole embodiment.Based on the embodiment among the present invention, those of ordinary skills belong to the scope of protection of the invention not making the every other embodiment that is obtained under the creative work prerequisite.
The invention discloses a kind of peak drift correction method, can improve the accuracy of peak drift correction.
Referring to Fig. 1, Fig. 1 is the process flow diagram of a kind of peak drift correction method disclosed by the invention.Comprise:
Step S1: the test spectrum to multichannel analyzer output is carried out smooth treatment, and the test that obtains after smooth is composed.
Neutron bombardment coal stream, nucleic in the coal stream and neutron react and produce gamma rays, and detector detects gamma rays, produces a pulse electrical signal, and this electric signal is imported in the multichannel analyzer through after amplifying.Multichannel analyzer is sampled to this pulse signal, gets access to the pulse height of this pulse signal, can convert digital signal to after the processing, and the ray counting of different-energy is placed in the specific road number, and satisfies quadratic function relation between road number and the energy substantially:
E=ax
2+ bx+c+ ε (formula 1)
Wherein, E is the gamma rays energy that is input to multichannel analyzer, and x represents number, and a, b and c are the coefficients (a>0 be function opening direction up) of quadratic polynomial, and ε is an error.
The test spectrum that multichannel analyzer output characterizes with road number and counting is carried out smooth operation to the test spectrum of multichannel analyzer output, can use 5 SMOOTHING FORMULAE methods in force, and the test spectrum is carried out 2 to 5 smooth operations.
Step S2: described test spectrum after smooth is carried out the peak-seeking operation, determine the actual road number corresponding respectively with carbon peak, hydrogen peak and plumbous peak.
Carry out in the smooth operation test spectrum afterwards, the characteristic peak of some nucleic can be clearly, the present patent application people is big and to disturb less characteristic peak be that principle is chosen 3 peaks in the test spectrum to obtain amplitude, as the energy of determining the test spectrum and the iteration original bulk of the functional relation between the road number (being formula 1).The present patent application people finds through test of many times: carry out peak-seeking with mentioned above principle, 3 peaks that select are respectively carbon peak, hydrogen peak and plumbous peak, these three peaks are the most obvious in the test spectrum, determine carbon peak, the pairing actual road in hydrogen peak and plumbous peak number in the test spectrum respectively.
The specific implementation of determining carbon peak, the pairing actual road in hydrogen peak and plumbous peak number as shown in Figure 2.Comprise:
Step S21: the standard track number of in standard spectrum, determining carbon peak, hydrogen peak and plumbous peak correspondence successively;
Step S22: the standard track number with carbon peak, hydrogen peak and plumbous peak correspondence is the center successively in the test spectrum, searches the peak maximum in second preset range;
Step S23: the road number of the peak maximum correspondence that will find out successively is respectively as the actual road number of carbon peak, hydrogen peak and plumbous peak correspondence.
Guaranteeing that detector operation environment and temperature meet under the test condition of industry standard, with the gamma rays that obtains can spectrogram as standard spectrum, the standard of the institute of peak drift correction just reference.In standard spectrum, determine pairing road, carbon peak number, as the standard track number at carbon peak.In the test spectrum, seek the carbon peak afterwards, be specially, standard track number with fixed carbon peak is that the peak maximum is searched in the center in second preset range, promptly searches the point of amplitude maximum in one section test spectrum being determined by the standard track number and second preset range at carbon peak, and this point is the carbon peak.Pairing road, carbon peak number is the actual road number at carbon peak.Obtain the actual road number at hydrogen peak and plumbous peak afterwards according to the method described above successively.
In force, second preset range is total road number of multichannel analyzer ± 2%.If multichannel analyzer is 4096 roads, then second preset range is ± 80 roads.
Step S3: with the peak energy at described carbon peak, hydrogen peak and plumbous peak and actual road number as current known quantity group.
Step S4: utilize described current known quantity group to carry out the second degree trinomial expression match, determine the road number of test spectrum and the current functional relation between the energy.
The peak energy at carbon peak, hydrogen peak and plumbous peak is a standard volume, as current known quantity group, utilizes current known quantity group to carry out the second degree trinomial expression match peak energy at carbon peak, hydrogen peak and plumbous peak and actual road number, determines coefficient a, b and c in the formula 1 with this.
Utilize principle of least square method to carry out the second degree trinomial expression match in force, detailed process is:
According to principle of least square method, the quadratic sum minimum of each road statistics error, promptly
Wherein, x
iBe the actual road number of nucleic characteristic peak, E
iBe the energy of nucleic characteristic peak, ε
iError for nucleic characteristic peak correspondence.
Ask local derviation to get to it:
Utilize least square method to ask undetermined coefficient a, b, c makes formula 1=0, formula 2=0, formula 3=0 can get system of equations:
Can convert thereof into matrix for ease of calculating:
Bring the peak energy and the actual road number at the carbon peak in the current known quantity group, hydrogen peak and plumbous peak into above-mentioned matrix, can solve a by calculating, b, the value of c obtains testing the road number of spectrum and the current functional relation between the energy, and promptly formula 1.
Step S5: determine in the characteristic peak other any peak, utilize the peak energy at described peak and the theoretical road number that described current functional relation is determined described peak.
Step S6: with the theoretical road number at described peak is that the peak-seeking operation is carried out at the center in first preset range, determines the actual road number at described peak, and the peak energy and the actual road number at described peak is added into current known quantity group.
Step S7: utilize described current known quantity group to carry out the second degree trinomial expression match, determine the functional relation between number and the energy, and be defined as current functional relation.
Step S8: judge whether the number of times that described second degree trinomial expression match is carried out reaches preset times, if then turn to step S9, otherwise turn to step S5.
Step S9: utilize current functional relation to determine the pairing energy in each road, and described test spectrum and standard spectrum are mated.
Step S5 to step S7 be an iterative process.Comprise multiple nucleic in the coal stream, another peak except carbon peak, hydrogen peak and plumbous peak in the selected characteristic peak in step S5, the peak energy at this peak is brought in the formula 1, a, b, c and E in this up-to-date style 1 are known quantity, ε is zero in addition, can determine the value of x, promptly determine the road number of this peak correspondence, because formula 1 is to determine through curve fitting, it is not actual test spectral curve, so the x that determine this moment is the theoretical road number of this peak correspondence, need in the test spectrum, carry out the peak-seeking operation, determine the actual road number of this peak correspondence.Detailed process is:
Theoretical road number with this peak is the center, carries out the peak-seeking operation in first preset range, promptly searches the point of amplitude maximum in one section test spectrum being determined by the theoretical road number and first preset range at this peak, and the road number that this point is corresponding is the actual road number at this peak.At this moment, obtained one group of peak energy and the actual road number corresponding again, should organize data and add current known quantity group with this peak.Utilize four groups of data in the current known quantity group to carry out second degree trinomial expression match once more afterwards, determine the road number of new test spectrum and the current functional relation between the energy, promptly determine coefficient a, b and c in the formula 1 based on principle of least square method.Through after iteration, the curve of the test spectrum that formula 1 is represented more approaches actual test spectral curve.
In force, first preset range be multichannel analyzer total road number ± 0.5%, if multichannel analyzer is 4096 roads, then second preset range is ± 20 roads.
Can be clear and definite be, the number of times that the second degree trinomial expression match is carried out is many more, the curve of the test of acquisition spectrum is the test spectral curve of approaching more reality just, preset times is set to 5 in the present invention.
After an iteration was finished, in the selected characteristic peak original another peak, and the peak energy at this peak is brought in the formula 1 that obtains through previous iteration, and the process of iteration obtains new functional relation.Repeat the second degree trinomial expression match, the number of times that carries out up to the second degree trinomial expression match reaches preset times, and the functional relation that last iteration is obtained is as the corresponding relation between road number and the energy in the final test spectrum.
The nucleic of coal stream mainly comprises carbon, hydrogen, lead, oxygen and iron, and in carrying out iterative process, the order of choosing of nucleic is preferably: the carbon list is escaped peak, oxygen peak, first iron peak and second iron peak, and is as shown in the table:
The principle that above-mentioned nucleic is chosen is, the amplitude of choosing in original characteristic peak one by one is big, disturb less peak.
After the second degree trinomial expression match of finishing preset times, utilize current functional relation (be last iterative process determine functional relation) to determine the pairing energy in each road, and will test spectrum and standard spectrum mates.The purpose of coupling is in order to guarantee to test total road number of spectrum and total road number of standard spectrum is consistent, and when total road number of test spectrum is counted less than total road of standard spectrum, will test to compose and extend; When total road of counting the overgauge spectrum when total road of test spectrum is counted, will test to compose and compress.The detailed process of coupling is:
At first calculate the test spectrum and carry out the pseudo-road sign indicating number X in each road after the peak drift correction
i,
Wherein, a, b, c are the coefficient of standard spectrum, E
iBe the energy in i road, E
I+1Be the energy in i+1 road, X
iIt is the pseudo-road sign indicating number in i road.
Calculate the counting Y of test spectrum afterwards through each the road correspondence after the peak drift correction
i,
Wherein, X
I+1Be the pseudo-road sign indicating number in i+1 road, y
iBe the counting in the i road before revising, y
I+1Counting for the i+1 road before revising.
In drift correction method in peak disclosed by the invention, after the initial function relation between energy and the road number in the definite test spectrum in the carbon peak in the characteristic peak, hydrogen peak and plumbous peak, by repeatedly second degree trinomial expression match, obtain the energy of test spectrum and the final funtcional relationship between the road number, afterwards by mating with standard spectrum, make that characteristic peak and the characteristic peak in the standard spectrum in the test spectrum are approaching as far as possible, improved the accuracy of revising.
Each embodiment adopts the mode of going forward one by one to describe in this instructions, and what each embodiment stressed all is and the difference of other embodiment that identical similar part is mutually referring to getting final product between each embodiment.For the disclosed device of embodiment, because it is corresponding with the embodiment disclosed method, so description is fairly simple, relevant part partly illustrates referring to method and gets final product.
The professional can also further recognize, the unit and the algorithm steps of each example of describing in conjunction with embodiment disclosed herein, can realize with electronic hardware, computer software or the combination of the two, for the interchangeability of hardware and software clearly is described, the composition and the step of each example described prevailingly according to function in the above description.These functions still are that software mode is carried out with hardware actually, depend on the application-specific and the design constraint of technical scheme.The professional and technical personnel can use distinct methods to realize described function to each specific should being used for, but this realization should not thought and exceeds scope of the present invention.
The method of describing in conjunction with embodiment disclosed herein or the step of algorithm can directly use the software module of hardware, processor execution, and perhaps the combination of the two is implemented.Software module can place the storage medium of any other form known in random access memory (RAM), internal memory, ROM (read-only memory) (ROM), electrically programmable ROM, electrically erasable ROM, register, hard disk, moveable magnetic disc, CD-ROM or the technical field.
To the above-mentioned explanation of the disclosed embodiments, make this area professional and technical personnel can realize or use the present invention.Multiple modification to these embodiment will be conspicuous concerning those skilled in the art, and defined herein General Principle can realize under the situation that does not break away from the spirit or scope of the present invention in other embodiments.Therefore, the present invention will can not be restricted to these embodiment shown in this article, but will meet and principle disclosed herein and features of novelty the wideest corresponding to scope.
Claims (6)
1. a peak drift correction method is characterized in that, comprising:
Test spectrum to multichannel analyzer output is carried out smooth treatment, and the test that obtains after smooth is composed;
Described test spectrum after smooth is carried out the peak-seeking operation, determine the actual road number corresponding respectively with carbon peak, hydrogen peak and plumbous peak;
With the peak energy at described carbon peak, hydrogen peak and plumbous peak and actual road number as current known quantity group;
Utilize described current known quantity group to carry out the second degree trinomial expression match, determine the road number of test spectrum and the current functional relation between the energy;
Determine in the characteristic peak other any peak, utilize the peak energy at described peak and the theoretical road number that described current functional relation is determined described peak;
With the theoretical road number at described peak is that the peak-seeking operation is carried out at the center in first preset range, determines the actual road number at described peak, and the peak energy and the actual road number at described peak is added into current known quantity group;
Utilize described current known quantity group to carry out the second degree trinomial expression match, determine the functional relation between number and the energy, and be defined as current functional relation;
Judge whether the number of times that described second degree trinomial expression match is carried out reaches preset times, if then utilize current functional relation to determine the pairing energy of each road number, and described test spectrum and standard spectrum are mated; If not, then turn to the step of determining other any peak in the characteristic peak.
2. method according to claim 1 is characterized in that, describedly determines that respectively the process of the actual road number corresponding with carbon peak, hydrogen peak and plumbous peak is specially:
In described standard spectrum, determine the standard track number of described carbon peak, hydrogen peak and plumbous peak correspondence successively;
Standard track number with described carbon peak, hydrogen peak and plumbous peak correspondence is the center successively in described test spectrum, searches the peak maximum in second preset range;
With pairing road, the peak maximum that finds out successively number respectively as the actual road number of described carbon peak, hydrogen peak and plumbous peak correspondence.
3. method according to claim 2 is characterized in that, described first preset range be multichannel analyzer total road number ± 0.5%, described second preset range be described multichannel analyzer total road number ± 2%.
4. method according to claim 1 is characterized in that, described preset times is 5 times.
5. method according to claim 3 is characterized in that, described characteristic peak comprises: carbon peak, hydrogen peak, plumbous peak, carbon list are escaped peak, oxygen peak, first iron peak and second iron peak.
6. method according to claim 1 is characterized in that, least square method is adopted in described second degree trinomial expression match.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100042706A CN102175703B (en) | 2011-01-10 | 2011-01-10 | Method for correcting peak shift |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100042706A CN102175703B (en) | 2011-01-10 | 2011-01-10 | Method for correcting peak shift |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102175703A true CN102175703A (en) | 2011-09-07 |
CN102175703B CN102175703B (en) | 2012-11-21 |
Family
ID=44518908
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011100042706A Expired - Fee Related CN102175703B (en) | 2011-01-10 | 2011-01-10 | Method for correcting peak shift |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102175703B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102819034A (en) * | 2012-09-13 | 2012-12-12 | 成都理工大学 | Energy disperse spectroscopy |
CN105628645A (en) * | 2015-12-28 | 2016-06-01 | 贵州中烟工业有限责任公司 | Infrared analysis model transfer method |
CN105717140A (en) * | 2016-01-29 | 2016-06-29 | 成都理工大学 | Chang'e 2 gamma spectrum and nonlinear spectrum drift calibration method based on energy correspondence |
CN109212583A (en) * | 2017-07-06 | 2019-01-15 | 核工业北京地质研究院 | A kind of multiple tracks gamma energy spectrometer temperature-compensation method based on fitting of a polynomial |
CN110208251A (en) * | 2019-06-20 | 2019-09-06 | 安徽创谱仪器科技有限公司 | Plasma emission spectroscopy interference correction method |
CN111443371A (en) * | 2020-04-20 | 2020-07-24 | 山东省科学院海洋仪器仪表研究所 | Method for judging peak drift of radionuclide in seawater |
CN112099073A (en) * | 2020-09-16 | 2020-12-18 | 北京华力兴科技发展有限责任公司 | Nuclide identification spectrometer |
CN115616010A (en) * | 2022-12-19 | 2023-01-17 | 合肥金星智控科技股份有限公司 | Material component detection method and detection device based on cross-belt neutron activation analysis |
CN117826224A (en) * | 2024-03-06 | 2024-04-05 | 山东省科学院海洋仪器仪表研究所 | Temperature determination method for radionuclide peak drift in seawater |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2219854A (en) * | 1985-09-18 | 1989-12-20 | Cogent Ltd | Coal analysis system |
CN1588020A (en) * | 2004-10-15 | 2005-03-02 | 南京大陆中电科技股份有限公司 | Coal element analytic method and on-line detecting equipment based on spectral library least square method |
CN101509872A (en) * | 2009-03-20 | 2009-08-19 | 清华大学 | Coal quality on-line detecting analytical method based on regression analysis |
-
2011
- 2011-01-10 CN CN2011100042706A patent/CN102175703B/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2219854A (en) * | 1985-09-18 | 1989-12-20 | Cogent Ltd | Coal analysis system |
CN1588020A (en) * | 2004-10-15 | 2005-03-02 | 南京大陆中电科技股份有限公司 | Coal element analytic method and on-line detecting equipment based on spectral library least square method |
CN101509872A (en) * | 2009-03-20 | 2009-08-19 | 清华大学 | Coal quality on-line detecting analytical method based on regression analysis |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102819034A (en) * | 2012-09-13 | 2012-12-12 | 成都理工大学 | Energy disperse spectroscopy |
CN105628645A (en) * | 2015-12-28 | 2016-06-01 | 贵州中烟工业有限责任公司 | Infrared analysis model transfer method |
CN105717140A (en) * | 2016-01-29 | 2016-06-29 | 成都理工大学 | Chang'e 2 gamma spectrum and nonlinear spectrum drift calibration method based on energy correspondence |
CN109212583B (en) * | 2017-07-06 | 2020-10-23 | 核工业北京地质研究院 | Multi-channel gamma energy spectrometer temperature compensation method based on polynomial fitting |
CN109212583A (en) * | 2017-07-06 | 2019-01-15 | 核工业北京地质研究院 | A kind of multiple tracks gamma energy spectrometer temperature-compensation method based on fitting of a polynomial |
CN110208251A (en) * | 2019-06-20 | 2019-09-06 | 安徽创谱仪器科技有限公司 | Plasma emission spectroscopy interference correction method |
CN111443371A (en) * | 2020-04-20 | 2020-07-24 | 山东省科学院海洋仪器仪表研究所 | Method for judging peak drift of radionuclide in seawater |
CN111443371B (en) * | 2020-04-20 | 2022-06-17 | 山东省科学院海洋仪器仪表研究所 | Method for judging peak drift of radionuclide in seawater |
CN112099073A (en) * | 2020-09-16 | 2020-12-18 | 北京华力兴科技发展有限责任公司 | Nuclide identification spectrometer |
CN112099073B (en) * | 2020-09-16 | 2023-08-11 | 北京华力兴科技发展有限责任公司 | Nuclide identification spectrometer |
CN115616010A (en) * | 2022-12-19 | 2023-01-17 | 合肥金星智控科技股份有限公司 | Material component detection method and detection device based on cross-belt neutron activation analysis |
CN117826224A (en) * | 2024-03-06 | 2024-04-05 | 山东省科学院海洋仪器仪表研究所 | Temperature determination method for radionuclide peak drift in seawater |
CN117826224B (en) * | 2024-03-06 | 2024-05-10 | 山东省科学院海洋仪器仪表研究所 | Temperature determination method for radionuclide peak drift in seawater |
Also Published As
Publication number | Publication date |
---|---|
CN102175703B (en) | 2012-11-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102175703B (en) | Method for correcting peak shift | |
Halzen et al. | Charge asymmetry of weak boson production at the LHC and the<? format?> charm content of the proton | |
Bernuzzi et al. | Horizon-absorption effects in coalescing black-hole binaries:<? format?> An effective-one-body study of the nonspinning case | |
CN111045072B (en) | Be applicable to CeBr 3 Gamma energy spectrum iteration spectrum solving method of detector | |
Kavanagh et al. | Probing WIMP particle physics and astrophysics with direct detection and neutrino telescope data | |
Castorina et al. | Excursion set peaks: the role of shear | |
Uher et al. | Peak fitting and identification software library for high resolution gamma-ray spectra | |
Florido et al. | New calibrations for estimating the N/O ratio in H ii regions | |
Enqvist et al. | Measurement and simulation of neutron/gamma-ray cross-correlation functions from spontaneous fission | |
Boelaert | Dijet Angular Distributions in Proton-Proton Collisions: At√ s= 7 TeV and√ s= 14 TeV | |
CN104392084A (en) | Standard substance characteristic-database establishing method and device on basis of terahertz frequency band | |
Horen et al. | Lifetime of the 3 1− state and octupole collectivity in Zr 96 | |
Schlattauer et al. | Calibration of gamma-ray detectors using Gaussian photopeak fitting in the multichannel spectra with a LabVIEW-based digital system | |
Calmet et al. | Characteristic limits in radioactivity measurement: From Currie’s definition to the international standard ISO-11929 publication | |
Ivanov et al. | Single spin asymmetry in open charm photoproduction and decay as a test of pQCD | |
Papp | A critical analysis of the experimental L‐shell Coster–Kronig and fluorescence yields data | |
Gunnink | Plutonium isotopic analysis of nondescript samples by gamma-ray spectrometry | |
Niechciol et al. | Fγ: A new observable for photon-hadron discrimination in hybrid air shower events | |
Blaauw et al. | Comparison of alternative methods for multiplet deconvolution in the analysis of gamma-ray spectra | |
Ruhter et al. | Application of cadmium-zinc-telluride detectors in 235U enrichment measurements | |
Geller et al. | Experiment to measure the increase in electron mass with velocity | |
Hobbs | DIS physics at the EIC and LHeC and connections to the future LHC and $\nu $ A programs | |
Ogren et al. | Demonstration of INCC6 for advanced list-mode data acquisition and analysis using ALMM | |
Liu et al. | A two parameter iterative approach of multiplication correction for assaying plutonium metal parts with neutron multiplicity counting | |
Biselli et al. | Studies of GPDs at Jefferson Lab: results and future experiments |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121121 Termination date: 20200110 |
|
CF01 | Termination of patent right due to non-payment of annual fee |