This invention relates generally to portable x-ray fluorescence (XRF) analyzers.
Portable XRF analyzers are used to detect elements present in a sample. A typical portable XRF analyzer includes an x-ray source for directing x-rays to the sample and a detector responsive to the x-rays emitted from the sample. An analyzer processes the output signals produced by the detector and divides the energy levels of the detected x-ray photons into several energy subranges by counts of the number of x-ray photons detected to produce a graph depicting the x-ray spectrum of the sample. The intensity at different energy levels corresponds to concentrations of different elements.
Portable XRF analyzers are known. See, e.g., the applicants' co-pending applications, U.S. patent application Ser. No. 11/582,038 filed Oct. 17, 2006 entitled “XRF System with Novel Sample Bottle”, and U.S. patent application Ser. No. 11/585,367 filed Oct. 24, 2006 entitled “Fuel Analysis System”, by one or more common inventors hereof and all of the same assignee, incorporated by reference herein. See also U.S. Pat. Nos. 6,501,825; 6,909,770; 6,477,227; and 6,850,592, all of which are incorporated by reference herein. Using a portable XRF analyzer, an operator can detect whether certain elements are present in a sample for use in such applications as, inter alia, alloy, ores and mineral analysis, security and law enforcement, environmental applications, artistic and historic works, biomedical and pharmaceutical applications, process chemistry, and the like. Another key use of portable XRF analyzers is to detect elements listed by the United States Consumer Products Safety Commission (CPSC) such as lead in toys and coatings, and the European Union Directive Restriction on the Use of Certain Hazardous Substances (RoHs). This Directive restricts the use of certain hazardous substances, such as lead (Pb), mercury (Hg), cadmium (Cd), chromium (Cr) and Bromine (Br), and the like, in manufactured electrical and electronic equipment.
Traditionally, it was difficult to analyze the low energy x-rays emitted by the atoms of elements having lower atomic numbers, e.g., elements between sodium (Na) and chlorine (Cl). This is because the lower energy of these x-rays is typically absorbed by the ambient atmosphere (e.g., air) or the material itself. The argon that is naturally present in the air is also fluoresced very efficiently and creates a source of background noise in the same region of the spectrum where some of these low atomic number elements like S and Cl are being measured. Until recently, in order to accurately analyze and detect these lower atomic number elements, the air between the analyzer window and the detector must be removed. This is done either by creating a vacuum or performing a helium (He) purge whereby the helium displaces the air between the analyzer window and the detector. The vacuum or the purge condition prevents the lower energy x-rays from being absorbed into the ambient atmosphere and increases the sensitivity of the XRF analyzer. See U.S. Patent Publication Nos. US 2008/0152079A1 and 2007/0269003 incorporated herein by this reference. In the later reference, there is a pressure sensor and a temperature sensor in the vacuum chamber for determining the density of the air inside the vacuum chamber. A small vacuum pump is used and readings are corrected based on different pressures in the vacuum chamber.
However, recent improvements in detector technologies including the advent of silicon drift detectors (SDDs) in portable XRF devices have allowed for the accurate measurements of low atomic number elements without the need for a vacuum or purge condition. The SDD technology can count at typically 10× higher rates, and have lower intrinsic noise. These detector improvements offset the air absorption problems mentioned previously and thus allow low atomic number elements to be measured more effectively in ambient conditions, i.e. without the need for vacuum or purge conditions. It remains true, however, that using a purge or vacuum condition that removes the intervening air does improve the quality of the analysis even when an SDD detector is employed. Several manufacturers of portable XRF devices now offer commercial portable XRF devices that measure the low atomic number elements with SDD detector technology using only an ambient air environment between the sample window and the detector and/or x-ray source.
A significant problem remains in the use of portable XRF analysis for low atomic number elements in ambient conditions. For the measured x-rays from Na (Z=11) to Ti (Z=22) and, in particular, for important alloy elements such as Mg, Al and Si which are very low in energy, the number of x-rays that travel from the sample to the detector vary based on the ambient air pressure. This is because higher or lower the ambient air pressure means more or less air molecules that can absorb the low energy x-rays from a low atomic number element. The concentration of low atomic number elements reported is proportional to the number of x-rays fluoresced and detected from each low atomic number element in the sample. For example, the reported concentration for Mg is proportional to the number of Mg x-rays detected from the sample. An increase or decrease in air pressure relative to the air pressure at time of factory calibration means that fewer Mg x-rays are detected (for increased air pressure relative to the air pressure at the time of calibration) or more Mg x-rays are detected (for decreased air pressure relative to the air pressure at the time of factory calibration). This means that the reported results for Mg will vary depending upon the local air pressure, thus creating systematic errors in the reported results.
Portable XRF units are typically calibrated at a specific geographical location, usually at the factory, and are therefore subject to only small variations in pressure. They are then used in a variety of climates and geographies where variations in ambient air pressure occur. In some cases, portable XRF units are taken to higher elevations, in the mountains, where the ambient air pressure is substantially lower. In other cases, portable XRF units can be taken hundreds or thousands of feet underground at a mine site. In all these cases, the actual local air pressure can be vastly different from the air pressure at the time of calibration. It is also true that even mild variations in barometric pressure, for example, a low pressure weather front, will affect the measured intensities of the lowest atomic number elements. Portable XRF analyzers are used because they provide fast, accurate quantitative results in the field. It is simply impractical to expect operators to perform manual corrections or recalibrations for changes in barometric pressure due to weather or altitude.
An object of this invention is to provide an improved portable XRF analyzer and such an analyzer which provides an ambient air pressure measurement that is reported to the XRF processor in order to correct the results for low atomic number elements in the field in order to account for the effects of altitude or local barometric pressure. The operator need not correct results or recalibrate the portable XRF unit at the local barometric conditions due to the effects of air pressure on low atomic number elements.
The subject invention, in part, results from the realization that, in one embodiment, an improved portable XRF analyzer automatically corrects results for low atomic number elements by the use of an on-board barometer to yield superior analytical data especially now that it is commercially feasible to employ an SDD to measure these elements without the use of a vacuum or purge condition. A detector behind a window is responsive to the x-rays radiated from the sample. A processor responsive to the detector analyzes the spectrum of emitted x-rays to detect the low atomic number elements and is responsive to a barometric sensor that detects a pressure change at the analyzer and automatically corrects the reported results.
A very compact, low power barometer such as the model MPL115A manufactured by Freescale Semiconductor can be mounted inside a portable XRF device. The barometer measures air pressure independent of type of gas or humidity and reports the value digitally to a processor. The portable XRF device can use the locally measured air pressure at any time and the measured air pressure at the time of calibration to make corrections to the measured intensities of the low atomic number elements for any given measurement. Thus, the reported concentrations for the low atomic number elements are corrected automatically by the portable XRF device for any climate or altitude.
The subject invention features, in one example, a portable analyzer comprising an x-ray source configured to emit x-rays to a sample and a detector subsystem responsive to x-rays irradiated by the sample and outputting the intensities of x-rays detected at different energy levels. A pressure measurement device is disposed to measure the ambient air pressure. A temperature sensor is also typically included. A processing subsystem is responsive to the detector subsystem and the pressure measurement device and is configured to calculate the concentration of at least one low atomic number element in the sample based on the intensity of the x-rays detected by the detector subsystem at an energy level corresponding to the element. The processing subsystem corrects the intensity based on the ambient pressure and temperature.
The detector subsystem typically outputs an intensity I for a known concentration of an element at a calibrated pressure Pc. The processing subsystem then increases the intensity when the gas pressure is greater than Pc and decreases the intensity when the gas pressure is less than Pc. Further included may be stored calibration data (curves or formulas) including at least, for one low atomic number element, an intensity level for a known concentration of the element at a known gas pressure. The processing subsystem is then configured to correct the intensity based on the stored calibration data. In one preferred embodiment, the detector subsystem includes a silicon drift detector and the pressure measurement device is a barometer. A temperature sensor may also be included and used to correct the pressure output by the barometer.
The processor is preferably configured to correct the measured pressure based on the measured temperature by correcting the measured pressure based on the measured temperature and a calibration temperature, to correct the intensity by determining a difference between the corrected pressure and a calibration pressure and using the difference to correct the intensity, and to determine a correction factor which is a function of a constant and the difference correction factor for each element may be empirically determined for a known excitation energy.
A commercial analyzer includes a housing about the x-ray source and the detector subsystem. The housing includes a window through which the x-rays pass to and from the sample. The pressure measurement device is disposed in the housing.
The subject invention also features an XRF analysis method comprising emitting x-rays to a sample, detecting x-rays irradiated by the sample, and measuring the intensities of x-rays detected at different energy levels. The ambient air pressure is measured. The concentration of a typically low atomic number element in the sample is automatically calculated based on the intensity of the x-rays detected at an energy level corresponding to the element. The intensity value is then automatically corrected based on the ambient air pressure before the concentration is calculated.
For an intensity I measured for a known concentration of an element at a calibrated pressure Pc, correcting includes increasing the intensity when the gas pressure is greater than Pc and decreasing the intensity when the gas pressure is less than Pc.
The method may further include storing calibration data including at least, for one element, intensity levels for a known concentration of the element at a different gas pressure. Correcting the measured intensity then includes basing the correcting on the stored calibration data. Temperature information can be used to correct a measured pressure.
This invention also features an XRF analysis method including using an XRF analyzer to induce fluorescence in a calibration sample with a known concentration of an element at a known calibration pressure and temperature. Detecting the fluorescence and storing count rates for the element at the known pressure and temperature in the XRF analyzer. The XRF analyzer is then used to induce fluorescence in a field sample with an unknown concentration of the element at ambient pressure and temperature. The fluorescence is detected and a count rate for the element is determined. The ambient pressure and temperature are measured and the determined count rate is corrected based on the ambient pressure and temperature and the stored count rates for the element at the calibration pressure and temperature.
Other objects, features and advantages will occur to those skilled in the art from the following description of a preferred embodiment and the accompanying drawings, in which:
Aside from the preferred embodiment or embodiments disclosed below, this invention is capable of other embodiments and of being practiced or being carried out in various ways. Thus, it is to be understood that the invention is not limited in its application to the details of construction and the arrangements of components set forth in the following description or illustrated in the drawings. If only one embodiment is described herein, the claims hereof are not to be limited to that embodiment. Moreover, the claims hereof are not to be read restrictively unless there is clear and convincing evidence manifesting a certain exclusion, restriction, or disclaimer.
In one example, portable XRF analyzer 10,
As shown in
Thus, in accordance with the subject invention, analyzer 10,
Processing subsystem 22 thus corrects a measured concentration based on the pressure output by barometer 26. In one example, memory 28 includes stored calibration data (curves or equations) including, at least for one low atomic number element, an intensity level for a known concentration of the element at a known gas pressure. The processing subsystem 22 is then configured to correct the measured intensity and concentration based on the stored calibration data. In general, when the detector system outputs an intensity I for a known concentration of a given element at a calibrated pressure P0, the processor subsystem increases the measured concentration when the measured gas pressure is greater than P0 and decreases the measured concentration when the gas pressure is less than P0.
Memory 28 may include, for an element as shown in
Processing subsystem 22, in one typical embodiment, is configured (e.g., programmed) to energize source 12 upon a command from the user, step 40,
P
m
*T
nom
/T
m where Tnom=300K. (1)
The corrected pressure is then used to correct the intensity value as output from the detector, step 49, using the data in database 28 (either a curve of intensity versus pressure or an equation for the curve). From the pressure corrected intensity value, the concentration of that element can be calculated, this process is then repeated for preferably all detected elements. The results are then typically displayed. It should be noted that to optimize the instrument's real-time measurement performance, this process need not be applied to all detected elements. Accordingly, it may be applied to only the detected elements of concern (e.g., those with a low atomic number) to minimize processing time.
In one example, the difference between the measurement and standard pressure (1013 mbar) (Pnominal) is used. The pressure difference (Pdiff) is:
P
diff
=P
nominal
−P
m(corr) (2)
using this value, the correction factor to apply to the intensity value detected by the device's x-ray detector subsystem is calculated:
correction factor=exp(−CE* Pdiff) (3)
where CE is an empirically determined factor obtained for each element Ei for a known excitation energy Eo. Pdiff=Po−PE. Values of CE greater than 1 results in an increase in intensity and values less than 1 results in decreased intensity.
The corrFactor number can be obtained by processor 22,
As shown at 70 in
Also, if Pm and/or Tm are not available for some reason, processing should not cease and so Pm and Tm can be set to nominal values Po or To (or they can be user defined inputs), steps 74-78. In one example, Po is standard pressure (1013 mbar) and To is standard temperature (300K). If these nominal values are used, as shown at step 79, the correction factor (see equation 3) is 1.
In most cases, Pm and Tm will be available and, in step 80, the effective pressure is calculated as discussed above with reference to equation 1. Step 82 shows how the correction factor is calculated from Po, PE, and CE. Using correction factor returned at step 84, steps 49-51,
One handheld XRF analyzer 10,
In an XRF analyzer, air absorption in the region between the sample, detector and x-ray tube causes a substantial attenuation of the signal for low energy x-rays. If the air pressure in the region changes relative to the air pressure at the time of the factory calibration, the accuracy of the measurement is impacted. The air can attenuate x-rays both from tube to the sample and from the sample to the detector. The amount of attenuation is stronger the lower the energy, so for example magnesium (1.25 keV) will be more heavily affected than aluminum (1.48 keV) or silicon (1.74 keV).
A method for in-field, automatic correction of results based on changing air pressure and temperature is disclosed. One preferred method utilizes a pressure sensor 26,
Absorption A(P, T)=exp(−ρ*K(E)/L) (4)
where K is the energy dependent absorption cross section (in cm2/g) evaluated at energy E, and L is the path length from either the sample 14 to the detector 18, or the tube 12 to the sample 14. The value of K and its energy dependence can be obtained from tabulated material data (one example is the program XCOM available from the National Institute of Science and Technology). The only unknown is the quantity p which is the gas density in units of g/cm3. This value is a function of both temperature and pressure, and for this we substitute in the ideal gas law given by:
ρ=P/(RT) (5)
Here P is the pressure in mBar, T is the temperature in Kelvin (absolute), and R is the gas constant expressed in the right choice of units including the molar mass of air. With this substitution for the equation for the absorption as a function of pressure and temperature it becomes:
A(P, T)=exp(−P*K(E)/(RTL)) (6)
The equation for the total absorption in the region between tube, sample and detector, due to air absorption of both source and fluoresced x-rays is below. This is given by the formula below and accounts for the pressure and temperature dependent absorption from source x-rays that travel from the tube towards the sample, and from the fluoresced x-rays that travel from low atomic number elements fluoresced in the sample to the detector:
In the above equation, the following quantities are defined as follows:
E0=The energy of the incident x-rays in keV. For simplicity assumed to be the average energy of the x-rays from the x-ray source. This assumption may be relaxed at the expense of a more complicated but well known formulation. Other assumptions may also be made.
Ei=The energy of the fluoresced x-ray being detected from element “i.”
LT=The path length from the tube outlet to the sample, in units of cm.
Lp=The path length from the sample to the detector in units of cm.
The final form of the correction is based on the difference between the air absorption during analyzer calibration, at known pressure and temperature values Pc and Tc, and the air absorption that occurs during a given sample measurement at pressure and temperature values Pm and Tm. This is the equation that is used to calculate the correction to the absorption at different pressure and temperature values compared to the time of calibration:
A(final)=exp−{[Pm/Tm*(K(Ei)/LD−K(Eo)/LT)*R]−[Pc/Tc*(K(Ei)/LD−K(E0)/LT)*R]} (8)
The correction can further be reduced if all measurements, including calibrations and unknowns are all corrected to a single nominal pressure and temperature, and thus nominal density, at the time of data processing. The temperature can also be expressed as an effective change in pressure which would result in the same density at the nominal temperature.
The absorption has been measured and compared to predicted values from the above formulas, as shown in
Ic(Mg)*Afinal(Mg), (9)
Where Afinal(Mg) is the absorption factor calculated from the above formula, evaluated at the energy of the magnesium K-alpha x-ray energy, at various values of pressure. The predicted rates for Al, Si and P are given by similar equations except that the absorption factors are determined using the energy of the Al, Si or P K-alpha x-ray energies, respectively. The pressure was provided by the internal pressure sensor included with the unit. As shown in
Assume the detector outputs an intensity of 1000 counts per second for Al (energy level 1.48 keV). The pressure output by the barometer is 100 kPa. The temperature is determined and the pressure of 100 kPa is corrected to 95 kPa based on the temperature. There is at this lower pressure, a 2.16% increased in transmission. The intensity level of 1000 counts is then corrected by 21.6 to 978.4 counts based on a corrected pressure of 950 kPa (1000−1000×21.6%).
In a portable XRF analyzer, the following steps can be executed in order to automatically correct certain elemental results at the time of measurement, due to changes in ambient air pressure, in an example.
The XRF analyzer is factory calibrated for various elements at a known ambient pressure and temperature Pc and Tc respectively, so that the count rates for these various elements Ic (i) for element (i) is known. Practically, this correction need only be applied to the low atomic number elements such as Mg, Al, Si, P and S, but it is recorded for all samples. Thus processor 22,
During field use, the XRF analyzer measures the count rates of various elements “i” Im(i), and also measures the ambient pressure and temperature values Pm and Tm with the onboard pressure 26 and temperature 25 sensors.
Processor 22 calculates the absorption factor Afinal(i) for each element “i” at the new pressure and temperature values Pm and Tm.
Processor 22 then calculates a correction to the calibration intensities of the pure element count rates for each element of interest from the product Im(i)=Ic(i)Afinal. For example, Im(Mg)=Ic(Mg)*Afinal.
The analyzer uses the new calibration intensity Im (i) in order to calculate a corrected concentration of element “i” that includes the effect of a change in pressure and temperature. This is also done for other elements of interest Al, Si and P for example.
Although specific features of the invention are shown in some drawings and not in others, this is for convenience only as each feature may be combined with any or all of the other features in accordance with the invention. The words “including”, “comprising”, “having”, and “with” as used herein are to be interpreted broadly and comprehensively and are not limited to any physical interconnection. Moreover, any embodiments disclosed in the subject application are not to be taken as the only possible embodiments. Other embodiments will occur to those skilled in the art and are within the following claims.
In addition, any amendment presented during the prosecution of the patent application for this patent is not a disclaimer of any claim element presented in the application as filed: those skilled in the art cannot reasonably be expected to draft a claim that would literally encompass all possible equivalents, many equivalents will be unforeseeable at the time of the amendment and are beyond a fair interpretation of what is to be surrendered (if anything), the rationale underlying the amendment may bear no more than a tangential relation to many equivalents, and/or there are many other reasons the applicant can not be expected to describe certain insubstantial substitutes for any claim element amended.