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.

Patent
   5412206
Priority
Feb 18 1994
Filed
Feb 18 1994
Issued
May 02 1995
Expiry
Feb 18 2014
Assg.orig
Entity
Large
7
46
all paid
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 claim 1, further comprising, prior to step (a), the step of effecting the emission of gamma rays by said gamma emitting element by irradiating said volume with neutrons.
3. A method as recited in claim 1, wherein said ratio is given by
R=(x+h2)2 /(x+h1)2,
wherein x represents the depth of the gamma emitting element.
4. A method as recited in claim 3, further comprising the steps of 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 said plot to infer said estimate of the depth of the gamma emitting element.
5. A method as recited in claim 4, further comprising the step of effecting the emission of gamma rays by said gamma emitting element by irradiating said volume with neutrons.
7. A system as recited in claim 6, further comprising means for effecting the emission of gamma rays by said gamma emitting element by irradiating said volume with neutrons.
8. A system as recited in claim 6, wherein said ratio is given by
R=(x+h2)2 /(x+h1)2,
wherein x represents the depth of the gamma emitting element.
9. A system as recited in claim 8, further comprising means for 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 said plot to infer said estimate of the depth of the gamma emitting element.
10. A system as recited in claim 9, further comprising means for effecting the emission of gamma rays by said gamma emitting element by irradiating said volume with neutrons.
12. A system as recited in claim 11, further comprising a neutron source, attached to said vehicle, for effecting the emission of gamma rays by said gamma emitting element.
13. A system as recited in claim 11, wherein said ratio is given by
R=(x+h2)2 /(x+h1)2,
wherein x represents the depth of the gamma emitting element.
14. A system as recited in claim 13, further comprising means for predetermining a plot of r versus x, and then comparing the ratio r measured for a source at an unknown depth with said plot to infer said estimate of the depth of the gamma emitting element.
15. A system as recited in claim 14, further comprising a neutron source, attached to said vehicle, for effecting the emission of gamma rays by said gamma emitting element.

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 onAssignorAssigneeConveyanceFrameReelDoc
Dec 13 1993RUDDY, FRANK H Westinghouse Electric CorporationASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS 0068960269 pdf
Dec 14 1993SEIDEL, JOHN G Westinghouse Electric CorporationASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS 0068960269 pdf
Dec 21 1993GONZALEZ, JOSEPH L Westinghouse Electric CorporationASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS 0068960269 pdf
Dec 21 1993CONGEDO, THOMAS V Westinghouse Electric CorporationASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS 0068960269 pdf
Feb 18 1994Westinghouse Electric Company(assignment on the face of the patent)
Mar 22 1999CBS CORPORATION FORMERLY KNOWN AS WESTINGHOUSE ELECTRIC CORPORATIONWESTINGHOUSE ELECTRIC CO LLCASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS 0100700819 pdf
Date Maintenance Fee Events
Oct 30 1998M183: Payment of Maintenance Fee, 4th Year, Large Entity.
Sep 24 2002M1552: Payment of Maintenance Fee, 8th Year, Large Entity.
Sep 26 2006M1553: Payment of Maintenance Fee, 12th Year, Large Entity.


Date Maintenance Schedule
May 02 19984 years fee payment window open
Nov 02 19986 months grace period start (w surcharge)
May 02 1999patent expiry (for year 4)
May 02 20012 years to revive unintentionally abandoned end. (for year 4)
May 02 20028 years fee payment window open
Nov 02 20026 months grace period start (w surcharge)
May 02 2003patent expiry (for year 8)
May 02 20052 years to revive unintentionally abandoned end. (for year 8)
May 02 200612 years fee payment window open
Nov 02 20066 months grace period start (w surcharge)
May 02 2007patent expiry (for year 12)
May 02 20092 years to revive unintentionally abandoned end. (for year 12)