A method for determining the depth of a gamma emitting element beneath the surface of a volume of soil is disclosed. The disclosed method includes the steps of detecting gamma rays at a first height h1 above the surface; detecting gamma rays at a second height h2 above the surface; determining a ratio r representative of the ratio of the number of gamma rays detected at h2 to the number of gamma rays detected at h1, or vice versa; and inferring, on the basis of the ratio r, an estimate of the depth of the gamma emitting element beneath the surface.
|
1. A method for determining the depth of a gamma emitting element beneath the surface of a volume of material, comprising the steps of:
(a) detecting gamma rays at a first height h1 above the surface; (b) detecting gamma rays at a second height h2 above the surface; (c) determining a ratio r representative of the ratio of the number of gamma rays detected at h2 to the number of gamma rays detected at h1, or vice versa; and (d) inferring, on the basis of the ratio r, an estimate of the depth of the gamma emitting element beneath the surface.
6. A system for determining the depth of a gamma emitting element beneath the surface of a volume of material, comprising:
(a) detecting means for detecting gamma rays at a first height h1 and a second height h2 above the surface; (b) means coupled to said detecting means for determining a ratio r representative of the ratio of the number of gamma rays detected at h2 to the number of gamma rays detected at h1, or vice versa; and (c) means for inferring, on the basis of the ratio r, an estimate of the depth of the gamma emitting element beneath the surface.
11. A system for determining the depth of a gamma emitting element beneath the surface of a volume of material, comprising:
(a) a gamma detector; (b) a vehicle that is movable over said surface, said vehicle carrying said gamma detector and including means for moving said gamma detector vertically in relation to said surface; and (c) a computer and associated support electronics for determining a ratio r of a number of gamma rays detected at a height h2 to a number of gamma rays detected at another h1, or vice versa, and for deducing, on the basis of the ratio r, an estimate of the depth of the gamma emitting element beneath the surface.
2. A method as recited in
3. A method as recited in
R=(x+h2)2 /(x+h1)2, wherein x represents the depth of the gamma emitting element. 4. A method as recited in
5. A method as recited in
7. A system as recited in
8. A system as recited in
R=(x+h2)2 /(x+h1)2, wherein x represents the depth of the gamma emitting element. 9. A system as recited in
10. A system as recited in
12. A system as recited in
13. A system as recited in
R=(x+h2)2 /(x+h1)2, wherein x represents the depth of the gamma emitting element. 14. A system as recited in
15. A system as recited in
|
The present invention relates generally to methods and apparatus for use in soil remediation, and more particularly to methods and apparatus for determining the depth of a gamma emitting element beneath the surface of a volume of material.
With the growing awareness of contamination of large tracts of land with chemically or radioactively hazardous elements, there is a corresponding international effort to initiate remediation activities to restore affected regions to an environmental status considered acceptable. To this end, soil washing and other methods are being developed. For these methods to be technically efficient and cost effective, it is necessary to accurately identify where the contamination in a field is located.
U.S. patent application Ser. No. 958,215, filed Oct. 8, 1992, titled "Prompt Gamma Neutron Activation Analysis System," which is hereby incorporated by reference, discloses the use of Prompt Gamma Neutron Activation Analysis (PGNAA) in soil remediation. Experimental and analytical data presented therein demonstrate that the disclosed methods and system are capable of measurements of trace element concentrations within a material sample by achieving extremely high signal-to-noise ratios.
PGNAA employs neutron-induced reactions. Such reactions can be divided into two broad categories--threshold reactions and exoergic reactions. For threshold reactions, energy in the form of neutron kinetic energy is required to supply a certain minimum energy to the reaction system before the reaction can proceed. Neutrons with energies below this minimum threshold energy are incapable of inducing the nuclear reaction. For exoergic reactions, the threshold is zero; that is, neutrons with all energies are capable of inducing the reaction.
FIG. 1 illustrates the process of neutron activation at a nuclear level. A neutron of energy E from a source 16 collides with the nucleus of an atom in the sample and initiates a reaction. For a neutron of thermal energy (0.0252 eV), the reaction might be absorption of the neutron into the nucleus, creating the next higher mass isotope of that element. If the neutron is more energetic (e.g., with several mega-electronvolts of kinetic energy), other nuclear reactions are possible. These other reactions include reactions wherein nuclear transmutation to another element occurs. In each of these cases, the residual nucleus is left in a highly excited internal state, and decays to its ground state almost instantaneously, emitting a gamma ray of several mega-electronvolts of energy. The energy of this gamma ray is uniquely characteristic of the quantum structure of the residual nucleus, and thus is a signature of the original target nucleus. The number of atoms of each of the elements of interest in a sample can be estimated by detecting (with a detector 12) and collecting the spectrum of gamma rays emitted by the sample and integrating the appropriate peaks.
Three methods for inferring the depth distribution of a contaminant element in a field of soil or other matrix material are disclosed by the above-cited U.S. patent application Ser. No. 958,215. These methods are summarized below:
1. The first method uses the inherent sensitivity bias of fast neutron-induced nuclear reactions for interrogating shallow depths (several inches) versus the deeper interrogation which is characteristic of thermal neutron-induced nuclear reactions.
2. The second method is based on the fact that, if a neutron-induced reaction on a particular nuclide of a contaminant element gives rise to more than one gamma ray, each with a known branching ratio, and if the energies of two or more such gamma rays are significantly different, the different amounts by which the gamma rays attenuate in propagating through the soil back to a detector will alter the ratio of their detected signals by an amount which is a unique signature of the average depth of the contaminant aggregate.
3. The third method makes use of the successive placement of none, one, or more reflectors above the neutron tube. The reflectors are large masses of highly moderating (low atomic mass constituent) materials. Graphite and polyethylene are examples. Such materials have the effect of utilizing neutrons which otherwise would have travelled upward from the neutron tube, significantly lowering their energy by collisions with the nuclei of the atomic constituents and reflecting them back toward the ground. The result is that the thermal neutron flux can be tailored to be highest in the shallow (several inches) depths of the soil with the reflectors, versus peaking at 6 to 10 inches (15.2 to 25.4 cm) below the soil surface in the absence of a reflector. Thus, even using only thermal neutron-induced nuclear reactions, it becomes possible to infer a depth distribution of a contaminant using thermal neutron PGNAA.
The present invention provides a fourth novel method for inferring the average depth of a contaminant using nuclear counting. The invention is applicable for detecting PGNAA induced gamma rays and naturally emitted gamma rays. The invention takes advantage of the fact that the number of gamma ray signals detected from a finite-size source in a fixed counting time interval decreases in a predictable fashion as the distance between the source and detector is varied.
According to the present invention, a method for determining the depth of a gamma emitting element beneath the surface of a volume of material comprises the steps of detecting gamma rays at a first height h1 above the surface; detecting gamma rays at a second height h2 above the surface; determining a ratio R representative of the ratio of the number of gamma rays detected at h2 to the number of gamma rays detected at h1, or vice versa; and inferring, on the basis of the ratio R, an estimate of the depth of the gamma emitting element beneath the surface.
In one presently preferred embodiment of the invention, the method further comprises the step of effecting the emission of gamma rays by the gamma emitting element by irradiating the volume of material with neutrons.
In addition, in the presently preferred embodiment the ratio is given by R=(x+h2)2 /(x+h1)2, wherein "x" represents the depth of the gamma emitting element. The depth estimation may be carried out by predetermining a plot of R versus x by implanting gamma emitting sources at known depths and measuring corresponding values of R, and then comparing the ratio R measured for a source at an unknown depth with the plot to infer the estimate of the depth of the gamma emitting element.
The present invention also provides a system including means for carrying out the method summarized above. One exemplary embodiment of such a system comprises a gamma detector; a vehicle that is movable over the surface, the vehicle carrying the gamma detector and including means for moving the gamma detector vertically in relation to the surface; and a computer and associated support electronics for determining a ratio R of a number of gamma rays detected at a height h2 to a number of gamma rays detected at another h1, or vice versa, and for deducing, on the basis of the ratio R, an estimate of the depth of the gamma emitting element beneath the surface. The system may advantageously include a neutron source, attached to the vehicle, for effecting the emission of gamma rays by the gamma emitting element.
Other features of the invention are described below.
FIG. 1 is a schematic diagram of a basic PGNAA system as known from the prior art.
FIG. 2 is a schematic diagram of one embodiment of a system for determining the depth of a gamma emitting element in accordance with the present invention.
FIG. 3 is a plot of the ratio of detected signals at 6 and 12 inch (0.15 and 0.30 meter) detector heights as a function of source depth. This plot indicates that the signal ratio is a unique function of source depth and is independent of photon energy.
FIG. 4 is a plot of expected view ratio R as a function of the depth x of a source beneath the surface. The view ratio is plotted for detector heights h1 =6 inches (0.15 meter) and h2 =12 (0.30 meter) and for several different source diameters. FIG. 4 depicts only the solid angle factor effect on the ratio R.
FIG. 5 illustrates the origin of a mild energy dependence in the ratio of yields measured when the detector is placed at different heights above the ground.
FIG. 2 schematically depicts one exemplary embodiment of a system for determining the depth of a gamma emitting element beneath the surface of a soil in accordance with the present invention. The system comprises a vehicle 10; gamma detector 12; support electronics 14, including a computer 14a; and an optional neutron source 16 for PGNAA applications The gamma emitting element, denoted "S", is shown at a depth "x" beneath the surface 18 of a soil (or other matrix material) volume 20. The gamma detector 12 is shown at a height "h1 " above the surface. The vertical double arrow above the gamma detector 12 indicates that the vehicle 10 is capable of moving the gamma detector to different heights in carrying out the inventive method described herein. The individual components depicted, including the vehicle 10, gamma detector 12, support electronics 14, and neutron source 16, are of common construction and, therefore, are not described in detail in this specification. If necessary, the reader may refer to the above-cited U.S. patent application for further information about gamma ray detection and PGNAA.
In presently preferred embodiments of the invention, measurements are made in several steps, with the gamma detector fixed at different heights above the surface. The present inventors have discovered that the average depth of the aggregate of gamma-emitting material can be deduced from the ratio of detected gamma intensity at different detector elevations above ground. It can also be demonstrated that, for localized sources, the intensity ratio versus depth relationship is independent of the energy of the detected photons or the matrix material composition.
The present invention is based on the phenomena experienced by a stream of gamma rays of energy E emitted isotropically by a gamma source S at a depth x beneath the surface of soil or some other matrix material (e.g., concrete, water, or other material). The source S could contain a naturally radioactive isotope or could be emitting prompt gamma rays as a result of PGNAA interrogation (e.g., by the neutron source 16 of the FIG. 2 embodiment).
First, assume that S is a point source, which is the simplest case. For an emission rate of N gammas per second of energy E from source S, the intensity of these gamma rays arriving at a narrow detector, located at an elevation h above the surface, is given by:
φ=N(e-μx /4π(x+h)2),
where μ is the linear attenuation coefficient (with units of inverse length) of the matrix material for photons of energy E.
The exponential in the numerator of the above expression corresponds to the attenuation of the gamma ray intensity as a result of travelling through a distance x of the matrix material, and the denominator (4π (x+h)2) corresponds to the inverse square falloff with distance of the intensity emitted from a point source. (This inverse square falloff is discussed below with reference to FIG. 4.) From this expression, one can see that, if a measurement of the rate of gamma detection is made with the gamma detector 12 placed at each of two detector heights, h1 and h2, the ratio R=φ(h1)/100 (h2) of these rates will be given by:
R=(x+h2)2 /(x+h1)2
This result for a point source depends only on the detector heights h1 and h2 and the source depth x below ground. For the special case of the source S at the surface, x=0 and the ratio is given by the square of the ratio of the detector heights. Conversely, if the ratio is measured for known heights h1 and h2, the depth x can be inferred. In measurements made using two different detector heights, the gamma rays from a localized source reaching the detector 12 have penetrated the same depth x of the matrix material. Therefore, the exponential term cancels when the ratio is formed. In this case, the inventive method is independent of the matrix material and the energy of the gamma ray.
As an example, FIG. 3 presents the expected ratio R for h1 =6 inches (15.2 cm) and h2 =12 inches (30.5 cm), as a function of the depth x of a point source located beneath the surface. The expected ratio is shown for gamma ray energies of 4 MeV and 1 MeV. The ratio is a unique function of x and is independent of the attenuation coefficient μ. This is indicated by the two sets of predictions for different gamma ray energies. The behavior for a source at the surface is supported by the following data, taken in a series of five minute counting intervals in the inventors' laboratory:
TABLE 1 |
______________________________________ |
Measured count rates and their ratios, from a series |
of five minute counting intervals, using a small 60 Co |
source, placed at detector heights of six inches and twelve |
inches. Ratios include minimum standard deviations. These data |
support the inverse square variation of intensity from a |
narrow source with distance to the detector. |
Measured Count |
Measured Count |
Measured Ratio |
60 Co |
Rate (cpm) for |
Rate (cpm) for |
R = CR (h1)/ |
Gamma h1 = 6 h2 = 12 CR(h2) (with |
Energy inches inches Std Dev) |
______________________________________ |
1173 keV |
1868 459 4.07 (± 0.09) |
1332 keV |
1654 424 3.90 (± 0.09) |
______________________________________ |
As the lateral dimensions of a below ground gamma source increase, a lower order dependence on depth follows, as indicated in FIG. 4. FIG. 4 depicts only the solid angle factor (geometric efficiency) effect on the ratio R, for a three-inch diameter detector held at heights of six inches (15.2 cm) and twelve inches (30.5 cm) above the surface. The inverse square dependence characteristic of a point source is also plotted for comparison. In addition, for an extended source, gamma rays from the outer regions will travel different depths in the absorbing medium, for the cases of different detector heights. This will introduce a weak dependence on photon energy and matrix material, and can be accounted for analytically. This is explained below with reference to FIG. 5. In field applications, the measurement of the lateral falloff of count rate with position can indicate which approximate source width should be used for depth inference.
As mentioned above, a weak energy dependence emerges in the ratio of yields measured when the detector 12 (FIG. 1) is placed at different heights above the ground (or surface of some other matrix material). This is a consequence of the fact that, for a source region which is below the soil surface and of finite width (i.e., of finite dimensions parallel to the surface), the gamma rays from the source perimeter region reaching the detector have travelled through differing depths of soil, in the cases of the detector placed at two different elevations above the ground. This is illustrated in FIG. 5, which depicts how a first source region "S1 " close to a vertical line connecting the two detector positions h1, h2 produces gamma rays that travel substantially the same distance through the soil 20 toward the two detector positions, whereas a second source region "S2 " far from the vertical line produces gamma rays that travel distances through the soil that differ by an amount "δ". Moreover, it is apparent that this effect is more severe for gamma sources nearer the surface 18. This effect may be readily modelled in commercially available transport codes that treat the gamma transport problem. The DOT and MCNP codes are two commercially available, internationally recognized examples of such code.
As mentioned, the present invention is not limited to methods and apparatus employing PGNAA, since it can be employed in connection with naturally radioactive sources. In addition, it is apparent that the invention may be employed in determining the depth of a radioactive source implanted in a volume of material other than soil, wherein the detector is placed above, below, or on a side of the volume, as the reference frame of the disclosed system may be rotated in accordance with known techniques. Accordingly, except where they are expressly so limited, the scope of protection of the following claims is not limited to the details specified above.
Congedo, Thomas V., Ruddy, Frank H., Gonzalez, Joseph L., Seidel, John G.
Patent | Priority | Assignee | Title |
11841335, | Apr 12 2018 | Riken | Nondestructive inspection method and apparatus comprising a neutron source and a gamma-ray detection device for determining a depth of a target component in an inspection target |
5608214, | Oct 30 1995 | Core Laboratories LP | Gamma ray spectral tool for well logging |
6518579, | Jun 11 1999 | Rensselaer Polytechnic Institute | Non-destructive in-situ method and apparatus for determining radionuclide depth in media |
7235780, | Apr 26 2002 | Rijksuniversiteit Groningen | Method and system for detecting a property of a pavement by measuring gamma-radiation |
7279683, | Apr 24 2002 | MITSUBISHI HEAVY INDUSTRIES, LTD; DIRECTOR-GENERAL OF THE INSTITUTE OF SPACE AND ASTRONAUTICAL SCIENCE, THE | Distance measurement apparatus of gamma ray source using multilayered ray detector |
7330128, | May 23 2005 | System and method for detecting radiological waste by trash collection vehicles | |
8680477, | Oct 15 2009 | IHI Corporation | Non-destructive inspection method and device |
Patent | Priority | Assignee | Title |
2395322, | |||
2557158, | |||
2562914, | |||
2781453, | |||
2800857, | |||
3008046, | |||
3124684, | |||
3341706, | |||
3354310, | |||
3433310, | |||
3444721, | |||
3463922, | |||
3483371, | |||
3508046, | |||
3715758, | |||
3723727, | |||
3781556, | |||
3801816, | |||
3832545, | |||
3889112, | |||
3918056, | |||
3938146, | Apr 01 1974 | Del Norte Technology, Inc. | Secure encoder for trilateralization locator utilizing very narrow acceptance periods |
4056969, | Apr 28 1975 | Barringer Research Limited | Detection of concealed metalliferous deposits, hydrocarbons and explosives |
4146796, | Mar 18 1977 | Apparatus for radiation source depth determination in a material | |
4232220, | May 04 1979 | Schlumberger Technology Corporation | Background subtraction system for pulsed neutron logging of earth boreholes |
4232317, | Nov 01 1978 | Quantized hyperbolic and inverse hyperbolic object location system | |
4248310, | Jun 28 1979 | Skip row cultivators | |
4278885, | May 04 1978 | Outokumpu Oy | Apparatus for measuring the concentrations of elements in a material by the capture gamma method |
4302285, | Feb 19 1979 | MALOE GOSUDARSTVENNOE PREDPRIYATIE METIST -USSR | Neutron activation analysis installation |
4314155, | Jul 19 1979 | Australian Atomic Energy Commission | Method and apparatus for elemental analysis employing combination of neutron inelastic scattering and γ ray scattering |
4317033, | Jun 15 1979 | Kenting Earth Sciences Limited; Hudson's Bay Oil and Gas Company Limited | Gamma ray prospecting system |
4350887, | Nov 05 1980 | The United States of America as represented by the United States | Method of assaying uranium with prompt fission and thermal neutron borehole logging adjusted by borehole physical characteristics |
4421981, | Apr 29 1980 | FORTHSTAR LIMITED, A BRITISH COMPANY | Method and apparatus for conducting remote surveys of the earth's crust |
4464330, | May 13 1982 | United States of America as represented by the United States Department of Energy | Apparatus for irradiating a continuously flowing stream of fluid |
4483817, | Jan 31 1983 | The United States of America as represented by the Administrator of the | Method and apparatus for mapping the distribution of chemical elements in an extended medium |
4581531, | Nov 16 1983 | Mobil Oil Corporation | Method for identifying hydrothermal alteration areas of the earth |
4645926, | Apr 10 1985 | WESTERN ATLAS INTERNATIONAL, INC , | Method for induced gamma ray logging |
4749858, | Aug 27 1986 | Troxler Electronic Laboratories, Inc. | Nuclear measuring gauge with automatic detection of source depth |
4754136, | Jan 17 1986 | CLEAN SOIL, INC ; CLEAN SOIL INC | Method of detecting underground tank leak |
4825073, | Dec 14 1987 | Halliburton Logging Services Inc. | Method for determining depth of penetration of radioactive tracers in formation fractures |
5025150, | Oct 14 1988 | WASHINGTON GROUP INTERNATIONAL, INC , AN OHIO CORPORATION | Site survey method and apparatus |
5038042, | Jan 16 1990 | Micron Technology, Inc | High resolution scintillation counters |
5133901, | Mar 01 1991 | WESTINGHOUSE ELECTRIC CORPORATION A CORP OF PENNSYLVANIA | System and method for on-line monitoring and control of heavy metal contamination in soil washing process |
5286973, | Nov 13 1991 | ODETICS, INC | Radiation mapping system |
EP81075, | |||
GB2148490, |
Executed on | Assignor | Assignee | Conveyance | Frame | Reel | Doc |
Dec 13 1993 | RUDDY, FRANK H | Westinghouse Electric Corporation | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 006896 | /0269 | |
Dec 14 1993 | SEIDEL, JOHN G | Westinghouse Electric Corporation | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 006896 | /0269 | |
Dec 21 1993 | GONZALEZ, JOSEPH L | Westinghouse Electric Corporation | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 006896 | /0269 | |
Dec 21 1993 | CONGEDO, THOMAS V | Westinghouse Electric Corporation | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 006896 | /0269 | |
Feb 18 1994 | Westinghouse Electric Company | (assignment on the face of the patent) | / | |||
Mar 22 1999 | CBS CORPORATION FORMERLY KNOWN AS WESTINGHOUSE ELECTRIC CORPORATION | WESTINGHOUSE ELECTRIC CO LLC | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 010070 | /0819 |
Date | Maintenance Fee Events |
Oct 30 1998 | M183: Payment of Maintenance Fee, 4th Year, Large Entity. |
Sep 24 2002 | M1552: Payment of Maintenance Fee, 8th Year, Large Entity. |
Sep 26 2006 | M1553: Payment of Maintenance Fee, 12th Year, Large Entity. |
Date | Maintenance Schedule |
May 02 1998 | 4 years fee payment window open |
Nov 02 1998 | 6 months grace period start (w surcharge) |
May 02 1999 | patent expiry (for year 4) |
May 02 2001 | 2 years to revive unintentionally abandoned end. (for year 4) |
May 02 2002 | 8 years fee payment window open |
Nov 02 2002 | 6 months grace period start (w surcharge) |
May 02 2003 | patent expiry (for year 8) |
May 02 2005 | 2 years to revive unintentionally abandoned end. (for year 8) |
May 02 2006 | 12 years fee payment window open |
Nov 02 2006 | 6 months grace period start (w surcharge) |
May 02 2007 | patent expiry (for year 12) |
May 02 2009 | 2 years to revive unintentionally abandoned end. (for year 12) |