WO2017018923A1 - Method for calibrating an x-ray fluorescence spectrometer - Google Patents
Method for calibrating an x-ray fluorescence spectrometer Download PDFInfo
- Publication number
- WO2017018923A1 WO2017018923A1 PCT/SE2016/050730 SE2016050730W WO2017018923A1 WO 2017018923 A1 WO2017018923 A1 WO 2017018923A1 SE 2016050730 W SE2016050730 W SE 2016050730W WO 2017018923 A1 WO2017018923 A1 WO 2017018923A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- calculated
- radiation
- background
- sample
- intensity
- Prior art date
Links
- 238000000034 method Methods 0.000 title abstract description 61
- 238000004876 x-ray fluorescence Methods 0.000 title description 4
- 230000005855 radiation Effects 0.000 claims abstract description 152
- 238000005259 measurement Methods 0.000 claims abstract description 75
- 238000013178 mathematical model Methods 0.000 claims abstract description 20
- 239000000463 material Substances 0.000 claims abstract description 14
- 239000000470 constituent Substances 0.000 claims abstract description 10
- 230000003595 spectral effect Effects 0.000 claims description 19
- 238000009826 distribution Methods 0.000 description 15
- 239000000203 mixture Substances 0.000 description 12
- 238000004364 calculation method Methods 0.000 description 10
- 238000012937 correction Methods 0.000 description 6
- 239000010405 anode material Substances 0.000 description 5
- 238000000611 regression analysis Methods 0.000 description 5
- 238000001228 spectrum Methods 0.000 description 5
- 238000013459 approach Methods 0.000 description 4
- 238000000691 measurement method Methods 0.000 description 4
- 239000010948 rhodium Substances 0.000 description 4
- 238000000441 X-ray spectroscopy Methods 0.000 description 3
- 238000007792 addition Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 239000011324 bead Substances 0.000 description 2
- 229910052703 rhodium Inorganic materials 0.000 description 2
- MHOVAHRLVXNVSD-UHFFFAOYSA-N rhodium atom Chemical compound [Rh] MHOVAHRLVXNVSD-UHFFFAOYSA-N 0.000 description 2
- 238000004611 spectroscopical analysis Methods 0.000 description 2
- 239000010421 standard material Substances 0.000 description 2
- 239000011573 trace mineral Substances 0.000 description 2
- 235000013619 trace mineral Nutrition 0.000 description 2
- 230000005461 Bremsstrahlung Effects 0.000 description 1
- 241001137251 Corvidae Species 0.000 description 1
- ZOKXTWBITQBERF-UHFFFAOYSA-N Molybdenum Chemical compound [Mo] ZOKXTWBITQBERF-UHFFFAOYSA-N 0.000 description 1
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 description 1
- 239000005864 Sulphur Substances 0.000 description 1
- 229920002678 cellulose Polymers 0.000 description 1
- 239000001913 cellulose Substances 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000010790 dilution Methods 0.000 description 1
- 239000012895 dilution Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052750 molybdenum Inorganic materials 0.000 description 1
- 239000011733 molybdenum Substances 0.000 description 1
- 235000015108 pies Nutrition 0.000 description 1
- 229920000136 polysorbate Polymers 0.000 description 1
- 230000005477 standard model Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000004846 x-ray emission Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/22—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by measuring secondary emission from the material
- G01N23/223—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by measuring secondary emission from the material by irradiating the sample with X-rays or gamma-rays and by measuring X-ray fluorescence
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/443—Emission spectrometry
Definitions
- the present invention relates to a method for calibrating an X-ray fluorescence spectrometer.
- the invention also relates to a method for measuring the composition of a sample by the use of an X-ray fluorescence spectrometer, which measurement method comprises such a calibration in an initial step.
- X-ray fluorescence spectroscopy is in general used to analyse the composition of matter, more specifically the contents of different elements in a sample.
- the sample is radiated with high-energy radiation, such as X-rays or gamma rays, from a radiation source.
- a radiation source is a radiation tube, which is arranged to accelerate electrons from a cathode towards an anode, using a known potential difference of about 10-100 keV.
- the electron impact will generate X-ray radiation, both a relatively specific part due to photoelectrical interaction with the anode material and a radiation part with relatively wide wave- length distribution as a result of bremsstrahlung.
- the radiation emitted from the radiation source impacts the sample, it in turn emits radiation of different origins.
- a respective characteristic pattern of photoelectrical ⁇ induced radiation will result, in the form of a respective characteristic spectral line pattern.
- a background radiation component which is due to scattering within the sample.
- the radiation emitted by the sample is then measured by a detector of a spectrometer.
- a detector of a spectrometer will detect radiation within a particular wavelength interval, but will also be prone to detecting radiation of neighbouring wavelengths, to a varying degree depending on the particular design of the spectrometer and its detector (line overlapping).
- the goal is to detect the relative amounts of different elements in the sample based on the intensity of the different photoelectric radiation peaks emitted therefrom. For elements present in the sample in larger proportions, this can be done with relatively high precision and accuracy. However, for elements with lower proportions, or even trace elements, the presence of the background radiation component will severely affect the measurement precision and accuracy.
- SE 0600946-8 describes a method in which a sample is investigated using spectrometry.
- J. E. Fernandez who is referred to in the said article, has developed a mathematical model allowing one to, given a set of particular characteristics concerning a combination of a radiation source and a spectrometer, calculate an expected detected radiation across a specified wavelength interval from a sample of known composition.
- This mathematical model is based upon the current knowledge of the physics relevant to the above described type of spectroscopy, which is detailed enough to provide reliable results.
- the Fernandez model is one example of such mathematical models, which are all ultimately based upon the standard model of physics.
- the invention relates to a method for calibrating a piece of measurement equipment for spectrometric measurement of constituent elements in a material sample, which piece of equipment comprises a source of radiation and a spectrometer detector, which method comprises, for different respective wavelength intervals, measuring the radiation intensity emitted from a number of known samples, each with a respective known proportion of a particular respective element, as well as calculating a respective expected detected radiation intensity across the same wavelength interval of each of said samples using a mathematical model, which method is characterised in that the method furthermore comprises the steps a) given a set of a plurality of different such samples, and for a particular background wavelength interval not comprising a spectral line of any element present to a meas- urable extent in the sample, determining a background relationship between measured and calculated radiation intensities; b) for each of the samples used in step a) and said background wavelength interval, calculating a respective relative difference between the calculated radiation intensity and a corresponding measured intensity calculated using said background relationship; and c) for a
- the invention relates to a system for measuring the composition of a material sample using a piece of measurement equipment for spectrometric measurement of constituent elements in such a material sample, which piece of equipment comprises a source of radiation and a spectrometer detector, which system is arranged to receive meas- urement values from the said spectrometer detector of the radiation intensity emitted from a number of known samples, each with a respective known proportion of a particular respective element, and for different respective wavelength intervals, as well as to calculate a respective expected detected radiation intensity across the same wavelength interval of each of said samples using a mathematical model, which system is characterised in that the system is arranged to, given a plurality different such samples, and for a particular background wavelength interval not comprising a spectral line of any element present to a measurable extent in the sample, determine a background relationship between measured and calculated radiation intensities; for each of the samples in said set and said background wavelength interval, calculate a respective relative difference between the calculated radi- ation intensity and a corresponding measured intensity calculated using said
- Figure la is a flowchart illustrating a calibration method according to the present invention
- Figure lb is a continuation of the flowchart of figure la, further showing a subsequent measurement part of the same method
- Figure 2 is a schematic overview of a setup with which a method according to the present invention can be used, wherein a system 100 according to the present invention is a part of said setup;
- ⁇ ⁇ wavelength interval comprising characteristic photoelectric spectrum line for basic element n.
- ⁇ ⁇ wavelength interval for background point x, not comprising any characteristic photoelectric spectrum line for any basic element being present, or at least not expected to be present in measurable quantities, in a sample for which the background point in question is selected, preferably not for any basic element.
- T m measured intensity
- T c calculated intensity
- T m (l) k ⁇ ⁇ T c ( ⁇ ) + ⁇ ⁇ linear relation describing the relation between measured and calculated intensity for the element or background point corresponding to wavelength interval ⁇ .
- IB c ( ⁇ ) calculated background intensity for wavelength interval ⁇ .
- T c ( ⁇ ) IB c ( ⁇ ) + IP c ( ⁇ ).
- IS c ( ⁇ ) component in calculated background intensity for waveguide interval ⁇ which is due to primary radiation scattering in a material sample.
- component in calculated background intensity for waveguide interval ⁇ which is
- FIG. 1 illustrates a method according to the present invention for calibrating a piece of measurement equipment for spectrometric measurement of constituent elements in a material sample, which method will now be discussed in detail.
- the piece of equipment comprises a source of radiation and a spectrometer detector.
- the source of radiation is a radiation tube, preferably using rhodium as the anode material.
- the spectrometer and its detector may be of any suitable type. Since such equipment is as such conventional and well-known, it is not described in detail herein.
- the calibration method comprises, for different respective wavelength intervals ⁇ , measuring the radiation intensity emitted from a number of known samples s pn , each with a respective known proportion p of a particular respective element n, as well as calculating a respective expected detected radiation intensity across the same wavelength interval ⁇ of each of said samples using a particular selected mathematical model.
- the said mathematical model is selected for calculating an expected detected radiation intensity by the spectrometer detector across a particular wavelength interval ⁇ , which radiation is emitted from a particular sample as a result of subjection of the sample to certain high-energy radiation from the radiation source.
- the high-energy ra- diation is preferably X-ray radiation.
- input parameters to the mathematical model comprise characteristics of the radiation source, such as radiation wavelength distribution; sample composition in terms of basic element proportions; spectrometer detector characteristics, such as spectrometer and spectrometer detector geometry; and wavelength detection interval.
- characteristics of the radiation source such as radiation wavelength distribution; sample composition in terms of basic element proportions; spectrometer detector characteristics, such as spectrometer and spectrometer detector geometry; and wavelength detection interval.
- spectrometer detector characteristics such as spectrometer and spectrometer detector geometry
- wavelength detection interval e.g., spectrometer and spectrometer detector geometry
- Such a model can be selected based on current knowledge about electromagnetic radiation, photoelectric effect, etc.
- the model may take into consideration the geometry and anode material of the radiation source; the geometry of the sample and the sample chamber; as well as the geometry and specifications of the spectrometer and its detector.
- the present calibration method is based upon the use of calibration standard samples of known and accurate composition with respect to basic element proportions. Such calibration samples are also well-known as such.
- the calibration method is commenced.
- one sample s pn at a time is subjected to high-energy radiation of known (such as calculated) characteristics, and the radiation intensity T m ( ⁇ ) emitted from the sample Spn in question is measured for a number of respective wavelength intervals ⁇ .
- These wave- length intervals preferably comprise at least one wavelength interval ⁇ ⁇ comprising a respective characteristic photoelectric spectrum line of one respective element n contained, to a known proportion p, in the sample s pn in question.
- the said wavelength intervals also preferably comprise at least one, preferably several different, wavelength intervals ⁇ ⁇ for a respective background point, not comprising such a characteristic photoelectric spectrum line for any of the elements n used in the calibration, preferably not comprising such a spectrum line for any element. It is realized that these measurements also can be performed later during the calibration method, as needed. However, it is preferred that all measurements on the known standards s P n are made at the offset of the calibration method, and that the method thereafter only comprises calculation steps. This provides for an efficient method, in which an initial measuring step can be proceeded by a calculation step using a high-performance computer. Preferably, the measurement on each known standard s pn are made in parallel in the above described piece of equipment.
- a first calibration sequence is initiated in order to, given a set of a plurality of different samples s pn , and for a particular background wavelength interval ⁇ ⁇ not comprising a spectral line of any element which is substantially present to a measurable extent in the sample, determine a background relationship (in terms of between measured and
- Output variables for this calibration step comprise, for each such background point ⁇ , regression constants describing a relationship between meas-
- interval ⁇ ⁇ constitutes an output variable of this calibration step, in the preferred case in which the interval is varied as described below during the calibration.
- a background point B x is selected, together with a corresponding wavelength interval comprising the background point B x
- a sample s pn is selected, and the intensity is calculated using the
- This step is repeated using a plurality of different samples s pn , preferably having a different respective known proportion p of one and the same basic element.
- samples s pn Preferably, more than three different samples s pn , with different proportions p and/or elements n, are used, such as at least 10 different samples s pn , preferably with at least two different elements n.
- a regression analysis is performed, with the aim of finding an optimal relationship, preferably a linear relationship, between measured and calculated intensities for the background point ⁇ ⁇ wavelength interval currently considered:
- regression analysis may be conventional as such, for instance using
- Output variables from the regression analysis comprise regression constants (a background regression constant representing a linear relationship between
- the wavelength interval ⁇ ⁇ in other words the pair of lower and upper wavelength limits, associated with the highest correlation between measured and calculated radiation intensities is selected as the wavelength interval ⁇ ⁇ for the background point B x in question to be used throughout the rest of the method.
- the wavelength interval ⁇ ⁇ in other words the pair of lower and upper wavelength limits, associated with the highest correlation between measured and calculated radiation intensities is selected as the wavelength interval ⁇ ⁇ for the background point B x in question to be used throughout the rest of the method.
- at least ten different intervals are tested.
- a "try all possibilities" approach is used, across a selection of potential lower and upper wavelength interval limits in the vicinity of the background point B x , and the best combination is then selected.
- a second calibration sequence is performed, according to the following.
- a respective relative difference is calculated between the calculated radiation intensity and a corresponding background point B x (preferably the same background point or background points as used in the above described first calibration sequence) and all samples s P n used in the first calibration sequence.
- an optimal, preferably linear, relationship is determined between measured T m and corrected calculated T mc intensities for each calibrated basic element n, preferably along with a corresponding optimal wavelength interval comprising a
- this third calibration sequence for a set of at least two different samples s pn comprising the same certain element n, and for a particular element wavelength interval ⁇ ⁇ comprising a spectral line of the certain element n, an element relationship (in terms of is determined between measured radiation intensities and
- a respective element n a background point B x a wave- length interval ⁇ ⁇ comprising a spectral line for the element n, and a sample s pn are selected.
- the intensity for the wavelength interval ⁇ ⁇ is calculated, using
- the calculated intensity T c is corrected, within the domain of the selected mathematical model, in the sense that the component of the
- tion intensities is not performed with respect to spectrometer detector measurement er- rors resulting from line overlaps.
- the selected mathematical model is hence used to calculate a respective background radiation component IB C of the various radiation intensities as a sum of radiation due to scattering in the sample IS C and spectral line overlap IO c as detected by the spectrometer detector.
- a respective corrected intensity T c is calculated in this way for at least two samples s pn comprising different proportions of one and the same element n, preferably for all samples s pn used in the second calibration sequence comprising known proportions of the element n in question currently considered in the second calibration sequence.
- a regression analysis is performed for the element n in question, with the aim of determining an optimal linear relationship between measured and corrected calculated intensities Note that k Xn is an element regression constant
- This regression analysis is, again, performed in a suitable way, such as using least mean squares, and is preferably performed for a number of different values for the wavelength interval ⁇ ⁇ , in a way which may be similar to the above described determination of an optimal value for the respective background points intervals ⁇ ⁇ .
- the pair of lower and upper wavelength limits associated with the highest correlation between measured and corrected calculated intensities is selected as the wavelength interval ⁇ ⁇ for the element n in question to be used throughout the rest of the method.
- at least ten different intervals are tested.
- a "try all possibilities" approach is used, across a selection of potential lower and upper wavelength interval limits in the vicinity of the element n spectral line, and the best combination (the one yielding the highest correlation) is then selected.
- the best combination the one yielding the highest correlation
- Such calibration analysis is performed for all elements n that are to be measured in subsequent measuring steps of unknown samples s u , and preferably for more than one background point ⁇ for each element n.
- the difference D is calculated as a weighted sum of the differences D Bx for each such background point are the weight fac
- the first and second calibration sequences are preferably performed for several different background wavelength intervals resulting in one respective difference
- the third calibration sequence is preferably performed using several different relative weight allocations for the said different background wavelength intervals ⁇ ⁇ , noted, ⁇ ⁇ , and the above described element relationship is determined for the particular relative weight allocation for which the correlation between the measured and corrected radiation intensities is the highest.
- the radiation source is a radiation tube
- the background wavelength interval ⁇ ⁇ comprises a peak of the Rayleigh or Compton scattered radiation from the sample of a characteristic line of the anode material of the source of radiation.
- two different background points ⁇ ⁇ and ⁇ ⁇ are used, preferably represented by non-overlapping respective wavelength intervals not comprising any photoelectric line of any element, and one of the background points preferably comprises a peak of such Rayleigh radiation, while the other one comprises a peak of such Compton radiation.
- one of the background points preferably comprises a peak of such Rayleigh radiation, while the other one comprises a peak of such Compton radiation.
- the determined calibration constants are stored, as applica
- the measurement on sample s u is performed using a piece of measurement equipment for spectrometric measurement of constituent elements in a material sample, which piece of equipment comprises a source of radiation and a spectrometer detector. More specifically, the measurement method comprises first performing a calibration method according to the above described, which calibration is performed with the same piece of equipment and then performing said measurement using a respective element relationship, in terms of as calibration constant for each one of a number of elements n.
- the measurement method comprises measuring the radiation intensity of the unknown sample s u , for at least one of the at least one background wavelength intervals ⁇ ⁇ used in the above described first calibration sequence, preferably the respective optimal background wavelength interval ⁇ ⁇ determined in that respective sequence run- through, and for at least two elements (in other words, at least two respective wavelength intervals ⁇ n) used in the above described third calibration step.
- This background point ⁇ ⁇ is selected.
- the radiation intensity is measured for the selected background
- the radiation intensity is calculated using the same se
- the measured radiation intensity is translated to calcu
- a relative difference is calculated according to I
- This first measurement sequence is performed for at least one background point pref erably for all background points used in the above described calibration sequences.
- n is selected.
- the radiation intensity for element n is calculated, and the result is corrected using the calculated difference
- r C ' (A n , s u ) constitutes an expected detected radiation intensity as corrected using the difference for the background point
- weight factors it is preferred that the same weight factors are used for the calculated difference D in this step as described above.
- the measured radiation intensity This is preferably performed using the element
- proportion p of the currently considered element n is calculated.
- this derivative is independent of the background radiation in the selected mathematical model.
- the final distribution is the element distribution determined by the method at this finishing point.
- the invention furthermore relates to a system 100, illustrated in figure 2, for measuring the composition of a material sample 210 of the above described kind, arranged in a measuring chamber 200, using the above-described piece of measurement equipment for spectromet- ric measurement of constituent elements in such a material sample.
- the piece of equipment comprises a radiation source 201 and a spectrometer detector 202, as described above.
- the equipment is controlled by a per se conventional control device 203.
- Such a system 100 is preferably connected to the piece of equipment, or at least to the spectrometer detector 202, and arranged to receive measurement values from the spectrometer detector 202 over an electric interface 102, preferably a standard interface, such as a digital data interface.
- the system 100 is arranged both to perform the calibration step sequences described above, as well as the measurement step sequences described above.
- the system comprises a user interface, such as in the form of a screen/keyboard interface 320 and/or an API (Application Programming Interface) 101 in order to output calibration- and measurement results, and to receive parameter input.
- An API 101 is, for instance, arranged to communicate with a server 310.
- the invention also relates to a computer software product arranged to be used for calibrating the above-described piece of measurement equipment.
- the software product is executable on a piece of hardware, such as on the said system 100, which is connected to said piece of measurement equipment and arranged to receive measurement values from the said spectrometer detector of the said radiation intensity.
- a software prod- uct is arranged to implement all the above described functionality in terms of calibration and measurement steps, and preferably also to operate the said user interface 320 and/or API 101.
- the present invention allows for precise and quick calibration of a piece of X-ray spectrom- etry equipment, and also for subsequent measurement in a corresponding precise and quick manner of low-proportion component elements in a sample.
- the method allows such precise calibration and measurement to be performed using initial measurement and subsequent calculations, which admits the use of parallel measurements across several channels and use of high-speed computing circuitry.
- MSXRF a mathematical model-based piece of software called MSXRF
- a setup similar to the one illustrated in figure 2 to perform parts of a method according to the present invention radiation intensities of Molybdenum (Mo) and scatter from a Rhodium (Rh) tube were measured for 26 geological standards.
- the standards contained very low concentrations of Mo, 0.094-15.5 ppm.
- the standards were prepared as beads containing 50% weight of standard material and 50% of cellulose. Intensities of both Rayleigh scatter and Compton scatter by the sample of the Rh k-alpha tube line were measured.
- the Compton scattered peak coincided with a Th line (L2-G1, 0.653 A. This line was used for the MSXRF calls as a line within the window is required by MSXRF).
- Optimum MSXRF wavelength window widths were determined in accordance with the above described, and intensity-intensity calibration lines were produced for three meas- ured intensities. Data of measured and calculated intensities produced during the calibration was transferred to a file which was imported to Excel for further calculations.
- the calculated background concentrations for Mo were now adjusted using the relative deviations from their calibration line of Rh_r and Rh_c for the corresponding standards.
- a slide control was included in the Excel sheet to be able to interactively find the optimum division of the correction between Rh_r and Rh_c. In this case, the best result was achieved when 30% of the Rh_r deviation and 70% of the Rh_c deviation was used.
- Regression calculations were performed to find a relation (calibration line) be- tween the adjusted calculated Mo intensities and the measured Mo intensities. Using this method, the standard deviation of the calculated intensities of Mo for the 26 standards from the calibration line was reduced to 1.2%. The standard deviation of the corresponding Mo concentrations from the certified values was 3.0 ppm.
- Table 1 shows the results for the different standards before and after correction. Namely, Table 1 shows the relative difference of the calculated intensity from the calibration line and this difference converted to ppm using the calculated intensity per cone derivative of the standard. The concentration deviations listed are the deviations for the standard material ("the sample”) so the deviation calculated for the measured bead has been multiplied with 2 (the dilution factor).
- Figures 3a and 3b show the measured/calculated intensity correlation before and after the Mo calculated background intensities had been adjusted using the deviations of Rh_r and Rh_c, respectively.
- the method can use other types of measurement equipment so long as the calibration principles described herein are used.
- a radiation tube instead of a radiation tube, another type of high energy radiation source can be employed.
- a calibration according to the present invention can be performed first, followed by several measurements according to the above on several different samples, preferably samples containing overlapping element compositions, preferably sam- pies prepared for measurement in the same way as the samples used during calibration.
- Method for calibrating a piece of measurement equipment for spectrometric measurement of constituent elements in a material sample which piece of equipment comprises a source of radiation (201) and a spectrometer detector (202), which method comprises, for different respective wavelength intervals ( ⁇ ), measuring the radiation intensity emitted from a number of known samples (s pn ), each with a respective known proportion (p) of a particular respective element (n), as well as calculating a respective expected detected radiation intensity across the same wavelength interval ( ⁇ ) of each of said samples using a mathematical model, c h a r a c t e r i s e d i n that the method furthermore comprises the steps
- the determination of the background relationship comprises calculating a background regression constant representing a linear relationship between measured and calculated
- step c) the determination of the element relationship comprises calculating an element regression constant representing a linear relationship between measured radiation intensities
- step a) is performed several times, wherein the background wavelength interval used for the calculation of the radiation intensities from the samples is modified to be
- step a different for different performances of step a), and in that the background relationship is determined for the particular modified background wavelength interval for which the
- step c) is performed several times, wherein the element wavelength interval ( ⁇ ⁇ ) used for the calculation of the radiation intensities from the samples is modified to be different for different performances of step c), and in that the element relationship is determined for the particular modified element wavelength interval ( ⁇ ⁇ ) for which the correlation between measured and corrected calculated radiation intensities is the highest.
- step c) is performed using several different relative weight allocations (yi,...yivi) for the said different background wavelength intervals ⁇ and in that the element relationship is determined for the relative weight allocation for which the correlation between the measured and corrected
- the source of radiation (201) is a radiation tube
- the background wavelength interval comprises a peak of the Rayleigh or Compton radiation
- the method comprises the sub steps i) assuming a particular element distribution of the sample
- step iv) is modified based upon the said derivative, and in that the iteration in step iv) is
- a piece of measurement equipment for spectrometric measurement of constituent elements in such a material sample which piece of equipment comprises a source
Abstract
Method for calibrating a piece of spectrometric measurement equipment of constituent elements in a material sample, which piece of equipment comprises a source of radiation and a spectrometer detector, which method comprises, for different wavelength intervals, measuring the radiation intensity emitted from a number of known samples, and calculating a corresponding expected detected radiation using a mathematical model. The method comprises the steps a) for a plurality of different samples, and a particular background wavelength, determining a background relationship between measured and calculated radiation intensities; b) for each sample and said background wavelength, calculating a difference between the calculated and calculated measured radiation intensity;and c) for two different samples comprising the same element, and for a particular element wavelength, determining an element relationship between measured and corrected calculated radiation intensities; wherein the element relationship is used as calibration constant. The invention also relates to a system and a computer software product.
Description
Method for calibrating an X-ray fluorescence spectrometer
The present invention relates to a method for calibrating an X-ray fluorescence spectrometer. The invention also relates to a method for measuring the composition of a sample by the use of an X-ray fluorescence spectrometer, which measurement method comprises such a calibration in an initial step.
X-ray fluorescence spectroscopy is in general used to analyse the composition of matter, more specifically the contents of different elements in a sample. The sample is radiated with high-energy radiation, such as X-rays or gamma rays, from a radiation source. One example of such a radiation source is a radiation tube, which is arranged to accelerate electrons from a cathode towards an anode, using a known potential difference of about 10-100 keV. The electron impact will generate X-ray radiation, both a relatively specific part due to photoelectrical interaction with the anode material and a radiation part with relatively wide wave- length distribution as a result of bremsstrahlung.
As the radiation emitted from the radiation source impacts the sample, it in turn emits radiation of different origins. For all elements contained in the sample, a respective characteristic pattern of photoelectrical^ induced radiation will result, in the form of a respective characteristic spectral line pattern. Apart from this photoelectric radiation, there will also be a background radiation component, which is due to scattering within the sample.
The radiation emitted by the sample is then measured by a detector of a spectrometer. Such spectrometer detector will detect radiation within a particular wavelength interval, but will also be prone to detecting radiation of neighbouring wavelengths, to a varying degree depending on the particular design of the spectrometer and its detector (line overlapping).
When using such a setup for analysing a sample, the goal is to detect the relative amounts of different elements in the sample based on the intensity of the different photoelectric radiation peaks emitted therefrom. For elements present in the sample in larger proportions, this can be done with relatively high precision and accuracy. However, for elements
with lower proportions, or even trace elements, the presence of the background radiation component will severely affect the measurement precision and accuracy.
In order to improve precision and accuracy, it has been proposed to carefully measure the background component radiated from the sample. However, such measurements are time consuming. For trace elements, careful background radiation measurements are very important for being able to correctly determine proportions of such elements. Typically, background radiation will be measured on both sides of a particular photoelectric peak for a particular element.
Hence, it would be desirable to be able to precisely calibrate a piece of X-ray spectrometry equipment in a precise and quick way. Furthermore, it would be desirable to be able to then perform accurate and precise measurements of low-proportion component elements in a sample in a less time-consuming way than what has previously been possible.
SE 0600946-8 describes a method in which a sample is investigated using spectrometry.
In the article "Application of J. E. Fernandez algorithms in the evaluation of X-ray intensities measured on fused glass discs for a set of international standards and a proposed calibra- tion procedure", J Malmqvist; X-RAY SPECTROMETRY; X-Ray Spectrom 2001;30:83-92, a method is described in which calibration is performed using a known standard sample.
The present invention solves the above described problems. J. E. Fernandez, who is referred to in the said article, has developed a mathematical model allowing one to, given a set of particular characteristics concerning a combination of a radiation source and a spectrometer, calculate an expected detected radiation across a specified wavelength interval from a sample of known composition. This mathematical model is based upon the current knowledge of the physics relevant to the above described type of spectroscopy, which is detailed enough to provide reliable results. The Fernandez model is
one example of such mathematical models, which are all ultimately based upon the standard model of physics.
The present invention relies on the use of such a mathematical model. However, since such models as such can be devised by those knowledgeable in the relevant physics, and since the example constituted by Fernandez' model is available, these models are not described in closer detail herein than what is necessary in order to understand the present invention. In the following, they are simply referred to as "models". Hence, the invention relates to a method for calibrating a piece of measurement equipment for spectrometric measurement of constituent elements in a material sample, which piece of equipment comprises a source of radiation and a spectrometer detector, which method comprises, for different respective wavelength intervals, measuring the radiation intensity emitted from a number of known samples, each with a respective known proportion of a particular respective element, as well as calculating a respective expected detected radiation intensity across the same wavelength interval of each of said samples using a mathematical model, which method is characterised in that the method furthermore comprises the steps a) given a set of a plurality of different such samples, and for a particular background wavelength interval not comprising a spectral line of any element present to a meas- urable extent in the sample, determining a background relationship between measured and calculated radiation intensities; b) for each of the samples used in step a) and said background wavelength interval, calculating a respective relative difference between the calculated radiation intensity and a corresponding measured intensity calculated using said background relationship; and c) for a set of at least two different samples comprising the same certain element, and for a particular element wavelength interval comprising a spectral line of the certain element, determining an element relationship between measured radiation intensities and calculated radiation intensities for the certain element, which calculated background radiation intensities have first each been corrected using said difference for the particular sample in question; wherein the hence determined element relationship is used as calibration constant for the certain element in question.
Furthermore, the invention relates to a system for measuring the composition of a material sample using a piece of measurement equipment for spectrometric measurement of constituent elements in such a material sample, which piece of equipment comprises a source of radiation and a spectrometer detector, which system is arranged to receive meas- urement values from the said spectrometer detector of the radiation intensity emitted from a number of known samples, each with a respective known proportion of a particular respective element, and for different respective wavelength intervals, as well as to calculate a respective expected detected radiation intensity across the same wavelength interval of each of said samples using a mathematical model, which system is characterised in that the system is arranged to, given a plurality different such samples, and for a particular background wavelength interval not comprising a spectral line of any element present to a measurable extent in the sample, determine a background relationship between measured and calculated radiation intensities; for each of the samples in said set and said background wavelength interval, calculate a respective relative difference between the calculated radi- ation intensity and a corresponding measured intensity calculated using said background relationship; and, for a set of at least two different samples comprising the same certain element, and for a particular element wavelength interval comprising a spectral line of the certain element, determine an element relationship between measured radiation intensities and calculated radiation intensities for the certain element, which calculated radiation intensities have first each been corrected using said difference for the particular sample in question, and in that the system is arranged to use the hence determined element relationship as calibration constant for the certain element in question during subsequent measurement on an unknown sample using the said piece of equipment. The invention also relates to a computer software product arranged to be executed on a piece of hardware, such as a system according to the invention, and arranged to perform a method according to the present invention.
In the following, the invention will be described in detail, with reference to exemplifying embodiments of the invention and to the enclosed drawings, wherein:
Figure la is a flowchart illustrating a calibration method according to the present invention; Figure lb is a continuation of the flowchart of figure la, further showing a subsequent measurement part of the same method;
Figure 2 is a schematic overview of a setup with which a method according to the present invention can be used, wherein a system 100 according to the present invention is a part of said setup; and
Figures 3a and 3b are respective diagrams showing experimental results of a run-through of a method according to the present invention. Below, the following variables and relationships will be used to describe various exemplifying embodiments of the present invention: n = a basic element (such as a metal or sulphur).
Bx = background point x
λη = wavelength interval comprising characteristic photoelectric spectrum line for basic element n.
λΒ = wavelength interval for background point x, not comprising any characteristic photoelectric spectrum line for any basic element being present, or at least not expected to be present in measurable quantities, in a sample for which the background point in question is selected, preferably not for any basic element.
Tm = measured intensity.
Tc = calculated intensity.
Tm(l) = kλ■ Tc(λ) + Ιλ linear relation describing the relation between measured and calculated intensity for the element or background point corresponding to wavelength interval λ.
IBc(λ) = calculated background intensity for wavelength interval λ.
/Pc(/l) = calculated photoelectric intensity for wavelength interval λ.
Tc(λ) = IBc(λ) + IPc(λ).
ISc(λ) = component in calculated background intensity for waveguide interval λ which is due to primary radiation scattering in a material sample.
component in calculated background intensity for waveguide interval λ which is
due to intensity additions from overlapping spectral lines.
tion. One example is the "background relationship" described below.
measured intensity is translated to calculated intensity via the linear relationship Tm (X) =
as determined during calibration for the wavelength interval λ.
spn = known standard sample with known proportion p of element n.
sp = known standard sample with known element distribution p.
su = unknown sample.
sa = assumed composition of unknown sample, in terms of the relative proportions of a set of elements in the sample.
Figure 1 illustrates a method according to the present invention for calibrating a piece of measurement equipment for spectrometric measurement of constituent elements in a material sample, which method will now be discussed in detail. The piece of equipment comprises a source of radiation and a spectrometer detector. Preferably, the source of radiation is a radiation tube, preferably using rhodium as the anode material. The spectrometer and its detector may be of any suitable type. Since such equipment is as such conventional and well-known, it is not described in detail herein.
As will be detailed in the following, the calibration method according to the present inven- tion comprises, for different respective wavelength intervals λ, measuring the radiation intensity emitted from a number of known samples spn, each with a respective known proportion p of a particular respective element n, as well as calculating a respective expected detected radiation intensity across the same wavelength interval λ of each of said samples using a particular selected mathematical model.
Thus, in a first step, the said mathematical model is selected for calculating an expected detected radiation intensity by the spectrometer detector across a particular wavelength interval λ, which radiation is emitted from a particular sample as a result of subjection of the sample to certain high-energy radiation from the radiation source. The high-energy ra- diation is preferably X-ray radiation. Preferably, input parameters to the mathematical model comprise characteristics of the radiation source, such as radiation wavelength distribution; sample composition in terms of basic element proportions; spectrometer detector characteristics, such as spectrometer and spectrometer detector geometry; and wavelength detection interval. Such a model can be selected based on current knowledge about electromagnetic radiation, photoelectric effect, etc. One example of such a model is the one developed by J. E. Fernandez and referred to in the above mentioned article. For instance, the model may take into consideration the geometry and anode material of the radiation source; the geometry of the sample and the sample chamber; as well as the geometry and specifications of the spectrometer and its detector.
The present calibration method is based upon the use of calibration standard samples of known and accurate composition with respect to basic element proportions. Such calibration samples are also well-known as such. In a next step, the calibration method is commenced.
In a next step, one sample spn at a time is subjected to high-energy radiation of known (such as calculated) characteristics, and the radiation intensity Tm( λ) emitted from the sample Spn in question is measured for a number of respective wavelength intervals λ. These wave- length intervals preferably comprise at least one wavelength interval λη comprising a respective characteristic photoelectric spectrum line of one respective element n contained, to a known proportion p, in the sample spn in question. The said wavelength intervals also preferably comprise at least one, preferably several different, wavelength intervals λΒ for a respective background point, not comprising such a characteristic photoelectric spectrum line for any of the elements n used in the calibration, preferably not comprising such a spectrum line for any element.
It is realized that these measurements also can be performed later during the calibration method, as needed. However, it is preferred that all measurements on the known standards sPn are made at the offset of the calibration method, and that the method thereafter only comprises calculation steps. This provides for an efficient method, in which an initial measuring step can be proceeded by a calculation step using a high-performance computer. Preferably, the measurement on each known standard spn are made in parallel in the above described piece of equipment. Then, a first calibration sequence is initiated in order to, given a set of a plurality of different samples spn, and for a particular background wavelength interval λΒ not comprising a spectral line of any element which is substantially present to a measurable extent in the sample, determine a background relationship (in terms of between measured and
calculated radiation intensities.
Hence, in this first calibration sequence an optimal linear relationship is determined between measured and calculated intensities for at least one of the above described background points Bx. Output variables for this calibration step comprise, for each such background point Βχ, regression constants describing a relationship between meas-
ured and calculated intensities across the background wavelength interval
interval λΒ constitutes an output variable of this calibration step, in the preferred case in which the interval is varied as described below during the calibration.
Hence, in a first step of said first calibration sequence, a background point Bx is selected, together with a corresponding wavelength interval comprising the background point Bx
in question.
selected mathematical model for the particular sample spn in question and the selected
background point Bx. This step is repeated using a plurality of different samples spn, preferably having a different respective known proportion p of one and the same basic element. Preferably, more than three different samples spn, with different proportions p and/or elements n, are used, such as at least 10 different samples spn, preferably with at least two different elements n.
In a next step, a regression analysis is performed, with the aim of finding an optimal relationship, preferably a linear relationship, between measured and calculated intensities for the background point λΒχ wavelength interval currently considered:
least square means. Output variables from the regression analysis comprise regression constants (a background regression constant representing a linear relationship between
measured and calculated radiation intensities) and Ιλ . Note that the intensity
was measured in the above described measurement step.
For each considered background point Bx, it is preferred that such regression constants are determined for several different wavelength intervals λΒχ, all comprising the background point Bx. Then, in a next step, the wavelength interval λΒχ, in other words the pair of lower and upper wavelength limits, associated with the highest correlation between measured and calculated radiation intensities is selected as the wavelength interval λΒχ for the background point Bx in question to be used throughout the rest of the method. Preferably, at least ten different intervals are tested. According to one preferred embodiment, a "try all possibilities" approach is used, across a selection of potential lower and upper wavelength interval limits in the vicinity of the background point Bx, and the best combination is then selected. However, it is also possible to use a more elaborate iterative goal-seeking approach to numerically find the optimal interval limits. Since no additional measurement is needed, and since these tests can be performed in full using only calculations, and without having to perform any additional measurements, the fine-tuning of the wavelength interval λΒ can be performed quickly.
These steps are repeated for at least one background point Bx, preferably for at least two background points.
Thereafter, for at least one background point Bx (preferably the same background point or background points as used in the above described first calibration sequence) and all samples sPn used in the first calibration sequence, a second calibration sequence is performed, according to the following. In this second calibration sequence, a respective relative difference is calculated between the calculated radiation intensity and a corresponding
measured intensity calculated using the background relationship determined in the first calibration sequence.
background point wavelength interval and sample combination in question is translated into a corresponding calculated intensity , using the previously established re-
In a third calibration sequence, an optimal, preferably linear, relationship is determined between measured Tm and corrected calculated Tmc intensities for each calibrated basic element n, preferably along with a corresponding optimal wavelength interval comprising a
respective photoelectric spectral line for each such element n.
According to the invention, in this third calibration sequence, for a set of at least two different samples spn comprising the same certain element n, and for a particular element wavelength interval λη comprising a spectral line of the certain element n, an element relationship (in terms of is determined between measured radiation intensities and
calculated radiation intensities for the said certain element n, which calculated background
radiation intensities have first each been corrected using the above described calculated relative difference for the particular sample spn in question.
Hence, in next subsequent steps, a respective element n, a background point Bx a wave- length interval λη comprising a spectral line for the element n, and a sample spn are selected.
the selected mathematical model. The calculated intensity Tc is corrected, within the domain of the selected mathematical model, in the sense that the component of the
calculated background intensity Tc constituted by primary radiation scattered by the sample Spn currently used is corrected using the difference previously calculated with respect to the background point Bx currently used: —
tion intensities is not performed with respect to spectrometer detector measurement er- rors resulting from line overlaps. In this preferred case, the selected mathematical model is hence used to calculate a respective background radiation component IBC of the various radiation intensities as a sum of radiation due to scattering in the sample ISC and spectral line overlap IOc as detected by the spectrometer detector.
A respective corrected intensity Tc is calculated in this way for at least two samples spn comprising different proportions of one and the same element n, preferably for all samples spn used in the second calibration sequence comprising known proportions of the element n in question currently considered in the second calibration sequence.
In a next step, a regression analysis is performed for the element n in question, with the aim of determining an optimal linear relationship between measured and corrected calculated intensities Note that kXn is an element regression constant
representing a linear relationship between measured radiation intensities and corrected calculated radiation intensities. This regression analysis is, again, performed in a suitable
way, such as using least mean squares, and is preferably performed for a number of different values for the wavelength interval λη, in a way which may be similar to the above described determination of an optimal value for the respective background points intervals λΒ . Hence, the pair of lower and upper wavelength limits associated with the highest correlation between measured and corrected calculated intensities is selected as the wavelength interval λη for the element n in question to be used throughout the rest of the method. Preferably, at least ten different intervals are tested. According to one preferred embodiment, a "try all possibilities" approach is used, across a selection of potential lower and upper wavelength interval limits in the vicinity of the element n spectral line, and the best combination (the one yielding the highest correlation) is then selected. However, it would also be possible to use a more elaborate goal-seeking approach, in order to numerically find the optimal interval limits.
Such calibration analysis is performed for all elements n that are to be measured in subsequent measuring steps of unknown samples su, and preferably for more than one background point Βχ for each element n.
In the preferred case that several such background points Bx are used, it is preferred that the difference D is calculated as a weighted sum of the differences DBx for each such background point are the weight fac
tors. Then, different sets of such weight factors are tried, and the set resulting in the highest correlation between measured and corrected calculated intensities is selected for each element n. This set of weight factors is then used for the calculation in the third calibration sequence for that element n.
Hence, the first and second calibration sequences are preferably performed for several different background wavelength intervals resulting in one respective difference
the third calibration sequence is preferably performed using several different relative weight allocations for the said different background wavelength intervals λΒι,
[...], λΒχ, and the above described element relationship is determined for the particular relative weight allocation for which the correlation between the measured and corrected radiation intensities is the highest. In a particularly preferred embodiment, the radiation source is a radiation tube, and, in the first calibration sequence, the background wavelength interval λΒ comprises a peak of the Rayleigh or Compton scattered radiation from the sample of a characteristic line of the anode material of the source of radiation. Preferably, two different background points λΒ and λΒ are used, preferably represented by non-overlapping respective wavelength intervals not comprising any photoelectric line of any element, and one of the background points preferably comprises a peak of such Rayleigh radiation, while the other one comprises a peak of such Compton radiation. This makes it possible to measure over a relatively short time period in the initial measuring step of the samples spn, and still to obtain measurement values with acceptable statistical accuracy.
ble, for use in subsequent measurement steps.
After the above described calibration sequences, the hence determined element relation- ship, in terms of preferably in combination with the determined background
relation, in terms of and/or the determined optimal wavelength intervals
and/or λη, may be used as calibration constants for each considered respective element n in a preferred method according to the invention for measuring the composition of an unknown sample su. The method steps of this preferred aspect of the invention is described in the following (see figure lb).
Like is the case for the above described calibration sequences, the measurement on sample su is performed using a piece of measurement equipment for spectrometric measurement of constituent elements in a material sample, which piece of equipment comprises a source of radiation and a spectrometer detector. More specifically, the measurement method comprises first performing a calibration method according to the above described, which
calibration is performed with the same piece of equipment and then performing said measurement using a respective element relationship, in terms of as calibration constant for each one of a number of elements n.
In general, the measurement method comprises measuring the radiation intensity of the unknown sample su, for at least one of the at least one background wavelength intervals λΒχ used in the above described first calibration sequence, preferably the respective optimal background wavelength interval λΒχ determined in that respective sequence run- through, and for at least two elements (in other words, at least two respective wavelength intervals λn) used in the above described third calibration step.
In the following, a particularly preferred such measurement method is described in detail.
Hence, in a first step of a first measurement sequence, a particular element ni, [...], ΠΝ dis- tribution pi, [...], pw of the sample su is assumed. This first element fraction pn is derived from the relative measured intensity of the element. The relative measured intensity is calculated as the element measured intensity of the current sample divided by the measured intensity corresponding to the calculated 100% intensity of the element (this calculated intensity is converted to a measured intensity via the element calibration constants).
selected from the background point or points already used in the above described first calibration sequence. In a next step, the radiation intensity is measured for the selected background
point and the sample
lected mathematical model as used in the calibration sequences.
In a next method step, the measured radiation intensity is translated to calcu
lated intensity using the above described relation
for the particular background point λΒ in question in the above described first calibration sequence.
Preferably, a relative difference is calculated according to I
This first measurement sequence is performed for at least one background point pref
erably for all background points used in the above described calibration sequences.
After this first measurement sequence, a high precision measure of the background radiation error is hence known, for the particular sample su and the particular measurement equipment used. In a subsequent second measurement sequence, the respective radiation intensity for each analysed element n is measured. It is realised that this measurement
can be performed beforehand, for instance in connection to the measurement of the background points in the first measurement sequence, preferably using a parallel measure
ment.
In a next step of the second measurement sequence, a particular element n is selected.
In a next step, the radiation intensity
for element n is calculated, and the result is corrected using the calculated difference
culated background intensity constituted by primary radiation scattered by the sample which is corrected, and additions from line overlap are not corrected for. Hence, rC' (An, su) constitutes an expected detected radiation intensity as corrected using the difference for the background point In case several background points were used, as
described above with weight factors
it is preferred that the same weight factors are used for the calculated difference D in this step as described above.
as described above:
In a next step, a difference is calculated between the corrected calculated radiation intensity
and the corresponding translated measured radiation intensity
proportion p of the currently considered element n is calculated. Preferably, this derivative is independent of the background radiation in the selected mathematical model.
bly according to
These second measurement sequence steps are performed for all considered elements n.
Then, in a next step, the method iterates back to the step in the second measurement sequence in which an element n is selected, now using the modified element distribution
The iteration is then repeated until the distribution is stable for all analysed ele
ments n. The final distribution is the element distribution determined by the method at this finishing point.
The invention furthermore relates to a system 100, illustrated in figure 2, for measuring the composition of a material sample 210 of the above described kind, arranged in a measuring chamber 200, using the above-described piece of measurement equipment for spectromet- ric measurement of constituent elements in such a material sample. The piece of equipment comprises a radiation source 201 and a spectrometer detector 202, as described above. The equipment is controlled by a per se conventional control device 203. Such a system 100 is preferably connected to the piece of equipment, or at least to the spectrometer detector 202, and arranged to receive measurement values from the spectrometer detector 202 over an electric interface 102, preferably a standard interface, such as a digital data interface. Preferably, the system 100 is arranged both to perform the calibration step sequences described above, as well as the measurement step sequences described above. Preferably, the system comprises a user interface, such as in the form of a screen/keyboard interface 320 and/or an API (Application Programming Interface) 101 in order to output calibration- and measurement results, and to receive parameter input. An API 101 is, for instance, arranged to communicate with a server 310. The invention also relates to a computer software product arranged to be used for calibrating the above-described piece of measurement equipment. The software product is executable on a piece of hardware, such as on the said system 100, which is connected to said piece of measurement equipment and arranged to receive measurement values from the said spectrometer detector of the said radiation intensity. Preferably, such a software prod- uct is arranged to implement all the above described functionality in terms of calibration
and measurement steps, and preferably also to operate the said user interface 320 and/or API 101.
The present invention allows for precise and quick calibration of a piece of X-ray spectrom- etry equipment, and also for subsequent measurement in a corresponding precise and quick manner of low-proportion component elements in a sample. In particular, the method allows such precise calibration and measurement to be performed using initial measurement and subsequent calculations, which admits the use of parallel measurements across several channels and use of high-speed computing circuitry.
In the following, an example is described, in which a mathematical model-based piece of software called MSXRF was used with a setup similar to the one illustrated in figure 2 to perform parts of a method according to the present invention. Hence, radiation intensities of Molybdenum (Mo) and scatter from a Rhodium (Rh) tube were measured for 26 geological standards. The standards contained very low concentrations of Mo, 0.094-15.5 ppm. The standards were prepared as beads containing 50% weight of standard material and 50% of cellulose. Intensities of both Rayleigh scatter and Compton scatter by the sample of the Rh k-alpha tube line were measured. For the instrument used, the Compton scattered peak coincided with a Th line (L2-G1, 0.653 A. This line was used for the MSXRF calls as a line within the window is required by MSXRF).
Optimum MSXRF wavelength window widths were determined in accordance with the above described, and intensity-intensity calibration lines were produced for three meas- ured intensities. Data of measured and calculated intensities produced during the calibration was transferred to a file which was imported to Excel for further calculations.
The results for Mo using conventional algorithms, not using the background point correction according to the present invention, was analysed in Excel. The standard deviation of the calculated Mo intensities from the calibration line was 3.5%. Using the intensity/concentration derivatives calculated during calibration and the calculated intensities for Mo,
the resulting standard deviation of the concentrations from the certified concentrations was 9.0 ppm.
The calculated background concentrations for Mo (overlap intensities excluded) were now adjusted using the relative deviations from their calibration line of Rh_r and Rh_c for the corresponding standards. A slide control was included in the Excel sheet to be able to interactively find the optimum division of the correction between Rh_r and Rh_c. In this case, the best result was achieved when 30% of the Rh_r deviation and 70% of the Rh_c deviation was used. Regression calculations were performed to find a relation (calibration line) be- tween the adjusted calculated Mo intensities and the measured Mo intensities. Using this method, the standard deviation of the calculated intensities of Mo for the 26 standards from the calibration line was reduced to 1.2%. The standard deviation of the corresponding Mo concentrations from the certified values was 3.0 ppm. During the measurements of the Mo intensities, around 60,000 pulses were captured for each sample. This means that the standard deviation of the intensity measurements due to the stochastics of the photon production was around 0.4%. The number of pulses captured for Rh_r, around 250,000 pulses and for Rh_c, around 500,000 pulses, was high enough to give just neligible contribution to the Mo intensity deviation. A longer Mo intensity meas- urement would reduce the 1.2% intensity deviation to some extent.
Table 1 shows the results for the different standards before and after correction. Namely, Table 1 shows the relative difference of the calculated intensity from the calibration line and this difference converted to ppm using the calculated intensity per cone derivative of the standard. The concentration deviations listed are the deviations for the standard material ("the sample") so the deviation calculated for the measured bead has been multiplied with 2 (the dilution factor).
Figures 3a and 3b show the measured/calculated intensity correlation before and after the Mo calculated background intensities had been adjusted using the deviations of Rh_r and Rh_c, respectively.
Above, preferred embodiments have been described. However, it is apparent to the skilled person that many modifications can be made to the disclosed embodiments without departing from the basic idea of the invention. For instance, the order and timing of measurements can be another one than as described above, both for the calibration sequences and the measurement sequences, as long as the measured values are available when to be used in the various method steps.
Furthermore, the method can use other types of measurement equipment so long as the calibration principles described herein are used. For instance, instead of a radiation tube, another type of high energy radiation source can be employed.
Although a quite detailed description of preferred algorithm steps have been described above, it is realized that the detailed implementation may differ in details as long as the principles as claimed are used.
Furthermore, it is realized that a calibration according to the present invention can be performed first, followed by several measurements according to the above on several different samples, preferably samples containing overlapping element compositions, preferably sam- pies prepared for measurement in the same way as the samples used during calibration.
Hence, the invention is not limited to the described embodiments, but can be varied within the scope of the enclosed claims.
C L A I M S
1. Method for calibrating a piece of measurement equipment for spectrometric measurement of constituent elements in a material sample (210), which piece of equipment comprises a source of radiation (201) and a spectrometer detector (202), which method comprises, for different respective wavelength intervals (λ), measuring the radiation intensity emitted from a number of known samples (spn), each with a respective known proportion (p) of a particular respective element (n), as well as calculating a respective expected detected radiation intensity across the same wavelength interval (λ) of each of said samples using a mathematical model, c h a r a c t e r i s e d i n that the method furthermore comprises the steps
a) given a set of a plurality of different such samples (spn), and for a particular background wavelength interval (λΒχ) not comprising a spectral line of any element present to a measurable extent in the sample, determining a background relationship between measured and calculated radiation intensities;
b) for each of the samples (spn) used in step a) and said background wavelength interval (λΒχ), calculating a respective relative difference between
the calculated radiation intensity and a corresponding measured intensity calculated using said background relationship;
c) for a set of at least two different samples (spn) comprising the same certain element (n), and for a particular element wavelength interval (λη) comprising a spectral line of the certain element (n), determining an element relationship between measured radiation intensities and calculated radiation intensities for the certain element (n), which calculated background radiation intensities have first each been corrected using said difference for the particular sample
[spn) in question;
wherein the hence determined element relationship is used as calibration constant for the certain element (n) in question.
2. Method according to claim 1, c h a r a c t e r i s e d i n that, in step a), the determination of the background relationship comprises calculating a background regression constant representing a linear relationship between measured and calculated
radiation intensities, and in that, in step b), said measured intensities are corrected using said background regression constant
3. Method according to claim 1 or 2, c h a r a c t e r i s e d i n that, in step c), the determination of the element relationship comprises calculating an element regression constant representing a linear relationship between measured radiation intensities
and corrected calculated radiation intensities.
4. Method according to any one of the preceding claims, c h a r a c t e r i s e d i n that step a) is performed several times, wherein the background wavelength interval used for the calculation of the radiation intensities from the samples is modified to be
different for different performances of step a), and in that the background relationship is determined for the particular modified background wavelength interval for which the
correlation between measured and calculated radiation intensities is the highest.
5. Method according to any one of the preceding claims, c h a r a c t e r i s e d i n that, in step c), the correction of the calculated radiation intensities is not performed with respect to spectrometer detector measurement errors resulting from line overlaps.
6. Method according to any one of the preceding claims, c h a r a c t e r i s e d i n that step c) is performed several times, wherein the element wavelength interval (λη) used for the calculation of the radiation intensities from the samples is modified to be different for different performances of step c), and in that the element relationship is determined for the particular modified element wavelength interval (λη) for which the correlation between measured and corrected calculated radiation intensities is the highest. 7. Method according to any one of the preceding claims, c h a r a c t e r i s e d i n that steps a) and b) are performed for several different background wavelength inter-
vals resulting in one respective difference for each such back
ground wavelength interval and each used sample, in that step c) is performed using several different relative weight allocations (yi,...yivi) for the said different background wavelength intervals
Μ and in that the element relationship is determined for the relative weight allocation for which the correlation between the measured and corrected
radiation intensities is the highest.
8. Method according to any one of the preceding claims, c h a r a c t e r i s e d i n that the model is used to calculate a background radiation component of the
various radiation intensities as a sum of radiation due to scattering in the sample (210)
and spectral line overlap as detected by the spectrometer detector
(202).
9. Method according to any one of the preceding claims, c h a r a c t e r i s e d i n that the source of radiation (201) is a radiation tube, and in that, in step a), the background wavelength interval comprises a peak of the Rayleigh or Compton radiation
from the sample (210) of a characteristic line of the anode material of the source of radiation (201). 10. Method for measuring the composition of an unknown sample (210, su) using a piece of measurement equipment for spectrometric measurement of constituent elements in a material sample (210), which piece of equipment comprises a source of radiation (201) and a spectrometer detector (202), c h a r a c t e r i s e d i n that the method comprises first performing a calibration method according to any one of the preceding claims of said piece of equipment and then performing said measurement using a respective element relationship as calibration constant for each one of a number of elements (n).
11. Method according to claim 10, c h a r a c t e r i s e d i n that the method comprises measuring, for at least one of the at least one background wavelength interval used in step a) and for at least two element wavelength intervals (λη) used in step c),
the radiation intensity of the unknown sample (210, su), wherein the method comprises the sub steps
i) assuming a particular element distribution of the sample
(210);
ii) for at least one particular of said background wavelength intervals used in
step a), calculating the radiation intensity using the model and correcting the calculated radiation intensity using the background relationship determined in step a), and calculating an difference between the calculated radia
tion intensity and the calculated intensity corresponding to the measured intensity;
responding element wavelength interval (λη), calculate an expected detected radiation intensity as corrected using the said difference determine a
calculated radiation intensity corresponding to the measured radiation intensity, and calculate a difference between the said corrected calculated radiation intensity and the said corresponding calculated radiation intensity; and iv) modify the assumed element distribution and iterate back to
step i).
12. Method according to claim 10 or 11, c h a r a c t e r i s e d i n that, in step ii), the correction of the calculated radiation intensity is not performed with respect to spectrometer detector measurement errors resulting from line overlaps.
13. Method according to any one of claims 10-12, c h a r a c t e r i s e d i n that, in step iii), for each element (n) in said element distribution the
derivative of the calculated radiation intensity with respect to the contents of the element (n) in question is calculated, in that, in step iv), the assumed element distribution
performed until the assumed element distribution stabilizes.
a piece of measurement equipment for spectrometric measurement of constituent elements in such a material sample which piece of equipment comprises a source
of radiation (201) and a spectrometer detector (202), which system (100) is arranged to
Claims
receive measurement values from the said spectrometer detector (202) of the radiation intensity emitted from a number of known samples (spn), each with a respective known proportion (p) of a particular respective element (n), and for different respective wavelength intervals (A), as well as to calculate a respective expected detected radiation intensity across the same wavelength interval (λ) of each of said samples using a mathematical model, c h a r a c t e r i s e d i n that the system (100) is arranged to, given a set of a plurality of different such samples (spn), and for a particular background wavelength interval (λΒχ) not comprising a spectral line of any element present to a measurable extent in the sample, determine a background relationship between measured a nd calculated radiation intensities; for each of the samples (spn) in said set and said background wavelength interval (λΒχ), calculate a respective relative difference between the calculated radia
tion intensity and a corresponding measured intensity calculated using said background relationship; and, for a set of at least two different samples (spn) comprising the same certain element (n), and for a particular element wavelength interval (λη) comprising a spectral line of the certain element (n), determine an element relationship between measured radiation intensities and calculated radiation intensities for the certain element (n), which calculated background radiation intensities have first each been corrected using said difference for the particular sample (spn) in question, and in that the system (100) is ar¬
ranged to use the hence determined element relationship as calibration constant for the certain element (n) in question during subsequent measurement on an unknown sample (su) using the said piece of equipment.
15. Computer software product arranged to be used for calibrating a piece of measurement equipment for spectrometric measurement of constituent elements in a material sample, which piece of equipment comprises a source of radiation (201) and a spectrometer detector (202), which software product is executable on a piece of hardware (100) which is connected to said piece of measurement equipment and arranged to receive measurement values from the said spectrometer detector (202) of the radiation intensity emitted from a number of known samples (spn), each with a respective known proportion (p) of a particular respective element (n), and for different respective wavelength intervals (λ), as well as to calculate a respective expected detected radiation intensity across the same wavelength
interval (λ) of each of said samples using a mathematical model, c h a r a c t e r i s e d i n that the software product is arranged to, given a set of a plurality of different such samples (spn), and for a particular background wavelength interval (λΒχ) not comprising a spectral line of any element present to a measurable extent in the sample, determine a background relationship between measured and calculated radiation intensities; for each of the samples (spn) in said set and said background wavelength interval (λΒχ), calculate a respective relative difference between the calculated radiation intensity and a
corresponding measured intensity calculated using said background relationship; and, for a set of at least two different samples (spn) comprising the same certain element (n), and for a particular element wavelength interval (λη) comprising a spectral line of the certain element (n), determine an element relationship between measured radiation intensities and calculated radiation intensities for the certain element (n), which calculated background radiation intensities have first each been corrected using said difference for the
particular sample (spn) in question, and in that the software product is arranged to make available the hence determined element relationship for use as calibration constant for the certain element (n) in question during subsequent measurement on an unknown sample (su) using the said piece of equipment.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
SE1551044A SE539014C2 (en) | 2015-07-28 | 2015-07-28 | Method for calibrating and X-ray fluorescence spectrometer |
SE1551044-9 | 2015-07-28 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2017018923A1 true WO2017018923A1 (en) | 2017-02-02 |
Family
ID=57886848
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/SE2016/050730 WO2017018923A1 (en) | 2015-07-28 | 2016-07-26 | Method for calibrating an x-ray fluorescence spectrometer |
Country Status (2)
Country | Link |
---|---|
SE (1) | SE539014C2 (en) |
WO (1) | WO2017018923A1 (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5982847A (en) * | 1996-10-28 | 1999-11-09 | Utah State University | Compact X-ray fluorescence spectrometer for real-time wear metal analysis of lubrucating oils |
WO2007126371A1 (en) * | 2006-04-28 | 2007-11-08 | Xrf Analytical Ab | Method at spectrometry for investigation of samples, where the sample contains at least two elements |
WO2010082897A1 (en) * | 2009-01-19 | 2010-07-22 | Xrf Analytical Ab | Method for spectrometry for investigation of samples containing at least two elements. |
US20150051877A1 (en) * | 2013-08-19 | 2015-02-19 | Kla-Tencor Corporation | Metrology Tool With Combined XRF And SAXS Capabilities |
-
2015
- 2015-07-28 SE SE1551044A patent/SE539014C2/en unknown
-
2016
- 2016-07-26 WO PCT/SE2016/050730 patent/WO2017018923A1/en active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5982847A (en) * | 1996-10-28 | 1999-11-09 | Utah State University | Compact X-ray fluorescence spectrometer for real-time wear metal analysis of lubrucating oils |
WO2007126371A1 (en) * | 2006-04-28 | 2007-11-08 | Xrf Analytical Ab | Method at spectrometry for investigation of samples, where the sample contains at least two elements |
WO2010082897A1 (en) * | 2009-01-19 | 2010-07-22 | Xrf Analytical Ab | Method for spectrometry for investigation of samples containing at least two elements. |
US20150051877A1 (en) * | 2013-08-19 | 2015-02-19 | Kla-Tencor Corporation | Metrology Tool With Combined XRF And SAXS Capabilities |
Non-Patent Citations (1)
Title |
---|
SITKO, R.: "Influence if X-ray tube spectral distribution on uncertainty of calculated fluorescent radiation intensity", SPECTROCHIMICA ACTA PART B, vol. 62, 2007, pages 777 - 786, XP002644867 * |
Also Published As
Publication number | Publication date |
---|---|
SE539014C2 (en) | 2017-03-21 |
SE1551044A1 (en) | 2017-01-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11210366B2 (en) | Analysis of X-ray spectra using fitting | |
US10921267B2 (en) | X-ray fluorescence analysis method, X-ray fluorescence analysis program, and X-ray fluorescence spectrometer | |
Donovan et al. | A new EPMA method for fast trace element analysis in simple matrices | |
Pella et al. | An analytical algorithm for calculation of spectral distributions of x‐ray tubes for quantitative x‐ray fluorescence analysis | |
Boumans | Detection limits and spectral interferences in atomic emission spectrometry | |
Alitti et al. | An improved determination of the ratio of W and Z masses at the CERN pp collider | |
JP2008309807A (en) | Method of determining background corrected count of radiation quantum in x-ray energy spectrum | |
US11402343B2 (en) | X-ray fluorescence spectrometer | |
US9188686B2 (en) | Radiation detector assembly and sample analyzer | |
CN108603849B (en) | Device for determining the effective energy spectrum of an X-ray tube | |
JPH10123071A (en) | Method and equipment for x ray analysis | |
JP3889187B2 (en) | X-ray fluorescence analysis method and apparatus | |
JP2005513478A5 (en) | ||
US10746675B2 (en) | Image processing device, analysis device, and image processing method for generating an X-ray spectrum | |
WO2017018923A1 (en) | Method for calibrating an x-ray fluorescence spectrometer | |
JP2006071311A (en) | X-ray fluorescence spectrometer and program used for it | |
JP5874108B2 (en) | X-ray fluorescence analyzer | |
Ménesguen et al. | Experimental determination of L fluorescence yields of gadolinium | |
Wegrzynek et al. | Evaluation of the uncertainty of element determination using the energy‐dispersive x‐ray fluorescence technique and the emission–transmission method | |
US11435303B2 (en) | X-ray analysis device and peak search method | |
Rackwitz et al. | Calculation of X‐ray tube spectra by means of photon generation yields and a modified Kramers background for side‐window X‐ray tubes | |
JP4523958B2 (en) | X-ray fluorescence analyzer and program used therefor | |
Šmit et al. | Parametrization of a tabletop micro‐XRF system | |
Liu et al. | Quantitative energy-dispersive X-ray fluorescence analysis for unknown samples using full-spectrum least-squares regression | |
Kaur et al. | Measurements of mass attenuation coefficients for 51Sb over energy region 4 keV–14 keV using synchrotron radiation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 16830924 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 16830924 Country of ref document: EP Kind code of ref document: A1 |