Method and System for Stack Monitoring of Radioactive Nuclides

Abstract

A system and method for monitoring one or more radioactive nuclides present in a stack flow consist of a first detector having a predetermined first sensitivity to gamma radiation and a second detector having a predetermined second sensitivity to gamma radiation and also a predetermined sensitivity to beta radiation. An enclosure proximal to the second detector defines a detection volume and enables the use of calibration factors which are independent of the geometry and material composition of a stack duct. A signal processor with energy window discrimination analyzes the signals from the two detectors. The use of two or more energy windows enables the identification of the nuclide species present in the stack flow and an accurate background-corrected measurement of the released radiation activity concentration for each of the identified nuclide species.

Claims

1. A system for monitoring one or more radioactive nuclides present in a stack flow, comprising: a first detector having a predetermined first sensitivity to gamma radiation; a second detector having a predetermined second sensitivity to gamma radiation and a predetermined sensitivity to beta particles produced by the one or more radioactive nuclides; a flow meter for measuring a flow rate of the stack flow; and an electronic signal processor which receives signals from said first detector, said second detector, and said flow meter; the electronic signal processor comprising energy window discrimination and configured to identify one or more nuclide species present in the stack flow, and to calculate a background-corrected value of radiation activity concentration for each of the identified nuclide species.

2. The system of claim 1 wherein said background-corrected value depends upon a ratio between the first and second sensitivities to gamma radiation.

3. The system of claim 1 further comprising an enclosure which is proximal to said second detector and which defines a detection volume.

4. The system of claim 1 wherein the first detector comprises a scintillator material selected from a group consisting of doped Sodium Iodide, doped Cesium Iodide, and Bismuth Germanate.

5. The system of claim 1 wherein the second detector comprises a scintillator material selected from a group consisting of Anthracene, Stilbene, and Naphthalene.

6. The system of claim 1 wherein the one or more radioactive nuclides comprise a positron-emitting nuclide.

7. The system of claim 1 wherein the one or more radioactive nuclides comprise a nuclide selected from a group consisting of Fluorine-18, Carbon-11, Nitrogen-13, Oxygen-15, and Gallium-68.

8. The system of claim 1 wherein said energy window discrimination comprises two or more energy windows, each defined by a lower limit on kinetic energy.

9. A method for monitoring a stack flow containing one or more radioactive nuclides, comprising the steps of: (a) providing an electronic signal processor with energy window discrimination which receives signals from a first detector having a first sensitivity to gamma radiation; a second detector having a second sensitivity to gamma radiation and a sensitivity to beta particles produced by the one or more radioactive nuclides; and a flow meter for measuring a flow rate of the stack flow; (b) pre-determining sensitivity calibration factors for the first and second detectors and for different nuclides; (c) calculating a background-corrected, beta-only signal for each energy window; (d) identifying one or more nuclide species present in the stack flow by comparing the beta-only signals of different energy windows; and (e) calculating a background-corrected value of the radiation activity concentration for each of the identified nuclide species.

10. The method of claim 9 wherein the calculation of a background-corrected beta-only signal in step (c) depends upon a ratio between the first and second sensitivities to gamma radiation.

11. The method of claim 9 wherein step (a) further comprises providing an enclosure which is proximal to said second detector and which defines a detection volume.

12. The method of claim 9 wherein step (e) further comprises calculating a background-corrected value of a total radiation activity for each of the identified nuclide species.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0022] The invention is herein described, by way of example only, with reference to the accompanying figures.

[0023] FIG. 1: Block diagram of an exemplary stack monitoring system, according to the invention.

[0024] FIGS. 2(a) and 2(b): Cross-sectional drawings of exemplary first and second detectors, respectively, according to the invention.

[0025] FIG. 3: Exemplary energy distribution produced by exposing the second detector to a Strontium 90 (Sr-90) beta-emitting calibration source.

[0026] FIG. 4: An exemplary graph of the positron energy emission distribution of several PET nuclides.

[0027] FIG. 5: A block diagram of an exemplary stack monitoring method, according to the invention.

DETAILED DESCRIPTION

[0028] Referring now to the drawings, FIG. 1 shows a diagram of an exemplary stack monitoring system 100, according to the invention. X-Y-Z reference axes indicate the orientation of the system 100 with respect to the duct walls 120. A fluid (e.g. air) containing radioactive nuclide effluents flows in the direction indicated by arrow 110, which is substantially parallel to the Z-axis. The cross-section of the duct may have a rectangular shape, a circular shape, or some other shape. First detector 130 and second detector 140 are placed back-to-back with their radiation-sensitive surfaces exposed to the fluid flow. Detector 130 detects gamma radiation with high efficiency, and consists for example of a high atomic number (high-Z) scintillator material, such as doped Sodium Iodide. Detector 140 is configured to detect both gamma rays and beta particles, and consists for example of a plastic scintillator material. Detectors 130 and 140 convert scintillator light pulses into electrical output signals 131 and 141, respectively, using optical and electrical components, such as light-guides, photocathodes, and photomultiplier tubes. Flow meter 170 measures the fluid flow velocity, for example in units of meters/second, and generates an output signal 171. Signals 131, 141, and 171 are transmitted to signal processor 180 by means of electrical leads, fiber optic cables, wireless communication, or any combination thereof. System 100 may include integrated power supplies (not shown in the figure) which provide electrical power to detectors 130 and 140, to flow meter 170, and to signal processing unit 180. Alternatively, some or all of the power supplies may be external to system 100.

[0029] The detection volume 160 of detector 140 is defined by an enclosure consisting of an enclosure plate 150 and an enclosure supporting frame 152, which has a negligible effect on air flow into volume 160. The size of volume 160 is proportional to the distance D between the plate 150 and the radiation sensitive surface of detector 140, as shown in FIG. 1. Defining the detection volume 160 as shown, the calibration factors used to calculate activity concentration may be pre-determined in a prototype setup having one type of duct and then normalized for use in other duct types having different size, shape, cross-sectional area, and/or material composition.

[0030] As an alternative to defining the detection volume by means of an enclosure plate and supporting frame, the detection volume can be determined by the duct itself. In this case, the calibration factors are calculated using either simulation (e.g. Monte-Carlo based simulation) or controlled injection of a known activity into the duct. The calibration factors are then functions of the duct's cross-sectional form and area.

[0031] FIG. 2(a) shows a cross-sectional drawing of an exemplary first detector 130, which is configured to detect primarily gamma radiation. One source of the gamma radiation may be from the decay of nuclides trapped on stack duct surfaces or charcoal filters inside the stack. Another, more distant source of gamma radiation may be short-term radiation events originating from outside the stack duct, for example in the chemical hot laboratory where the radioactive isotopes are produced and prepared for shipping. These distant sources of gamma radiation are considered background radiation, since they do not originate from activity that flows through the stack and is released into the environment. Scintillator 132 is preferably made of a high atomic number (high-Z) material, such as doped Sodium Iodide (NaI(Tl)), doped Cesium Iodide (CsI(Tl)), or Bismuth Germanate (BGO). Scintillator 132 is surrounded by a cylindrical detector shell 134, and the entire assembly is protected by a weatherproof detector housing 136. Shell 134 is preferably made of Aluminum having an exemplary thickness of 0.1 to 0.3 mm.; housing 136 is preferably made of stainless steel (e.g. SS-3041) having an exemplary thickness of about 0.03 to 0.08 mm.

[0032] FIG. 2(b) shows a cross-sectional drawing of an exemplary second detector 140, which is configured to detect both gamma radiation and beta particles (in this case, positrons) emitted by a radioactive isotope effluent. Scintillator 142 is preferably made of a low atomic number (low-Z) material with a high beta stopping power, such as organic crystals of Anthracene, Stilbene, or Naphthalene. Insofar as positrons generally have short mean free paths in matter, typically less than 1 mm., maximum detection sensitivity is achieved when scintillator 142 is in direct contact with the air flow. Detector shell 144 and weatherproof detector housing 146 are typically similar in shape and material composition to shell 134 and housing 136 of detector 130. Detector 140 includes a light-tight thin layer 148, which consists, for example, of a thin aluminized Mylar film. Layer 148 enables positrons to pass through to the scintillator while blocking light photons that would otherwise cause scintillator noise.

[0033] FIG. 3 shows an exemplary energy distribution produced by exposing the second detector to a beta-emitting Strontium 90 (Sr-90) calibration source. The detailed distribution was measured using a high-resolution multi-channel analyzer. The scale of the vertical axis is logarithmic.

[0034] When more than one species of nuclide is present in the stack flow, it is necessary to identify the radiation activity concentration of each individual nuclide species. Radioactive isotopes produce beta particles (e.g. positrons) whose kinetic energy spectra are characteristic of the isotope species. For example, FIG. 4 shows normalized positron energy emission distributions for five nuclides commonly used in positron emission tomography (PET). The high energy tails terminate at maximum kinetic energies (Emax) given in the following table.

TABLE-US-00001 TABLE 1 Nuclide Species Nuclide Symbol Emax (keV) Fluorine-18 F-18 634 Carbon-11 C-11 960 Nitrogen-13 N-13 1199 Oxygen-15 O-15 1732 Gallium-68 Ga-68 1920

[0035] Detector 140 is designed to be sensitive to low levels of incident radiation (for example, 1000 Bq/m.sup.3) from the nuclides present in the stack, while being relatively less sensitive to 511 keV gamma rays. The material composition and thickness of scintillator 142 are selected in accordance with previously determined nuclide stopping distances for incident kinetic energies up to Emax.

[0036] Electronic signal processor 180 supports energy windows discrimination, enabling nuclide identification. The energy windows W1-W5 are shown graphically in FIG. 5 by vertical lines. Exemplary lower and upper limits for the energy windows are given in the following table.

TABLE-US-00002 TABLE 2 Energy Windows Energy Lower Limit Upper Limit Nuclide Window (keV) (keV) Identified W1 500 800 F-18 W2 800 1100 C-11 W3 1100 1500 N-13 W4 1500 1800 O-15 W5 1800 open window Ga-68
The ranges of the pre-determined energy windows are chosen so that the Emax values associated with the different nuclides of interest (as given in TABLE 1) all fall in different energy windows.

[0037] FIG. 5 shows a block diagram of an exemplary stack monitoring method, according to the invention. The steps of the method are as follows: [0038] 610: provide a first detector having a first sensitivity to gamma radiation; [0039] 620: provide a second detector having a second sensitivity to gamma radiation and a sensitivity to beta particles produced by one or more radioactive nuclides; [0040] 630: provide a flow meter for measuring a flow rate of a stack flow; [0041] 640: provide an electronic signal processor with energy window discrimination which receives signals from the first detector, the second detector, and the flow meter; [0042] 650: pre-determine sensitivity calibration factors for the first and second detectors (e.g. η's, F, and sensitivity ratio R) and for different nuclides; [0043] 660: calculate a background-corrected, beta-only signal for each energy window; [0044] 670: identify nuclide species present in the stack flow by comparing the beta-only signals of different energy windows; and [0045] 680: calculate a background-corrected value of the radiation activity concentration for each of the identified nuclide species, using the background-corrected beta-only signals of step 660 and the pre-determined sensitivity calibration factors of step 650.
In step 660, the calculation of the background-corrected beta-only signal depends upon a pre-determined ratio, R, between the first and second sensitivities to gamma radiation. Step 620 may optionally include providing an enclosure which is proximal to the second detector, and which defines a detection volume. Furthermore, step 680 may optionally include calculating a background-corrected value of the total released radiation activity for each of the identified nuclide species.

Example

[0046] The following example provides additional specific details of the method of the invention, by way of example only. The exemplary method provides identification and measurement for each radioactive nuclide present in a stack flow.

[0047] Prior to the installation and use of a system 100 inside an isotope manufacturing facility, sensitivity calibration factors are pre-determined for the detectors 130 and 140 (method step 650), as follows. A prototype system is constructed according to the invention and placed inside a calibration duct. Air flowing over a calibrated radiation source, corresponding to one of the nuclides of interest in stack monitoring, such as F-18, is introduced into the calibration duct. Flow meter 170 measures the air velocity in units of m/s. By multiplying the flow rate by the cross-sectional area of the duct, the volume flow rate M is calculated, for example in units of cubic meters per second. When an enclosure plate 150 and enclosure supporting frame 152 are installed, the detection volume 160 is a known constant and the volume flow rate (M) depends only on the air velocity measured by the flow meter. Without the enclosure plate 150 and the frame 152, the detection volume is a detector length dimension multiplied by the duct cross-sectional area.

[0048] Prior to installation of the system 100 inside the duct 3, sensitivity calibration factors are pre-determined for the detectors 130 and 140, including: [0049] i. sensitivity of the first detector to gamma rays, denoted by η.sub.g-g; [0050] ii. sensitivity of the second detector to gamma rays, denoted by η.sub.p-g; and [0051] iii. sensitivity of the second detector to beta particles (e.g. positrons), denoted by η.sub.p-p.
The sensitivities are typically expressed in units of counts per second (cps) per millirad per hour (cps per mR/hr), or as dimensionless percentages (%) which are energy dependent. Note that the first detector is blocked for positrons by the second detector 140, so that the sensitivity of the first detector to positrons is effectively zero.

[0052] An additional calibration factor, F, is calculated as follows. Signal processor 180 receives output signal 131 as well output signals 141 for each of the five energy windows W1 through W5, all in units of cps. A known activity, in units of Bq, is introduced into the duct, and a calibration factor F, is then calculated, in units of Bq/m.sup.3 per cps or nCi/m.sup.3 per cps. The value of F is determined by dividing the known activity by the product of output signal 131 in cps, and the previously determined detection volume in cubic meters. The calibration process is repeated for other calibrated radiation sources, corresponding to other nuclides of interest in stack monitoring.

[0053] The calibration process is repeated using other ducts having different characteristics from those of the calibrating duct, such as size, shape or material composition, in order to determine whether these characteristics influence the values of the sensitivity calibration factors. For example, it may be necessary in some cases to adjust the sensitivity calibration factors by a dimensionless factor which depends upon the ratio between the cross-sectional area of a particular duct and the cross-sectional area of the calibration duct.

[0054] Furthermore, the sensitivity calibration factors of detectors 130 and 140 may require periodic (e.g. annual) maintenance, for example, by inserting a calibrated nuclide source into the duct and checking the accuracy of the measured activity concentrations.

[0055] After calibration of the prototype system, the system is ready to be used as a stack monitor. To calculate a background-corrected beta-only signal (method step 660), the following measurements are acquired in each energy window: [0056] i. background count rate (B.sub.g) of the first detector; [0057] ii. background count rate (B.sub.p) of the second detector; [0058] iii. actual count rate (C.sub.g) of the first detector; and [0059] iv. actual count rate (C.sub.p) of the second detector.
Background count rates are measured under conditions when there is no radiation present; actual count rates are measured with the stack flow on.

[0060] The measurements in each energy window are analyzed to determine which of the following scenarios applies: [0061] i. No background radiation is present and no activity is released. In this case,

C.sub.p≤B.sub.p+K.sub.pσ(B.sub.p)

and

C.sub.g≤B.sub.g+K.sub.gσ(B.sub.g)  where σ(B.sub.g) and σ(B.sub.p) denote standard deviations of the background count rate measured in the first and second detectors, respectively. K.sub.p and K.sub.g are pre-determined constants for adjusting the detection thresholds corresponding to positrons and gamma rays, respectively. Exemplary values are K.sub.p=5 and K.sub.g=3. [0062] ii. Radiation is present but no activity is released. This scenario occurs when the source of the gamma radiation is outside of the stack. For example, this may occur when isotope production is in operation or when a radiation source is being moved in the vicinity of the detectors. In this case,

[00001]$C_p\ge B_p+K_p\sigma \left(B_p\right)$$\mathrm{and}$$\frac{C_p}{C_g}<\frac{{\eta }_{p-g}}{{\eta }_{g-g}}=R$  where R=η.sub.p-g/η.sub.g-g denotes the gamma sensitivity ratio between the two detectors. In the second inequality, forming the ratio C.sub.p/C.sub.g corrects for the effect of background gamma radiation. Note that, if log-units are used the division is replaced by subtraction, that is, log(C.sub.p/C.sub.g)=log(C.sub.p)−log(C.sub.g). The second inequality is satisfied when the ratio of the measured count rates in the two detectors is less than the gamma sensitivity ratio R. [0063] iii. Radiation is present and activity is released. In this case, the nuclide species must be identified and the radiation activity concentration must be measured. The required conditions are as follows:

[00002]$C_p\ge B_p+K_p\sigma \left(B_p\right)$$\mathrm{and}$$\frac{C_p}{C_g}>\frac{{\eta }_{p-g}}{{\eta }_{g-g}}$$\mathrm{and}$$\frac{C_p-C_p.Math.{\eta }_{p-p}}{C_g}=K\frac{{\eta }_{p-g}}{{\eta }_{g-g}}$  where the pre-determined constant K typically has a value close to one. The numerator on the left hand side of the last equation corresponds to the measured count rate in the second detector (C.sub.p) reduced by the count rate (C.sub.P η.sub.p-p) which is contributed by positrons. If the last equation is not satisfied, then there is additional gamma radiation outside the stack, which does not contribute to an activity release, and should therefore be subtracted from the count rate of the first detector.

[0064] When the conditions of scenario iii. are met, there is an activity release from the stack and the nuclide species is identified (method step 670) as follows. A count rate that is only present in energy window W1 indicates the presence of F-18. Similarly, a count rate that is only present in window W2, W3, W4, or W5 indicates the presence of C-11, N-13, O-15, or Ga-68, respectively. The count rate of the identified species is then converted to a radiation activity concentration (in method step 680), denoted by A and having units of Bq/m.sup.3. The value of A is equal to the product of C.sub.g and the calibration factor F. The total released radiation activity for the identified species is denoted by U and has units of Bq. The value of U is equal to the product of the radiation activity concentration, A, the volume air flow rate, M, in units of m.sup.3/sec, and the time duration, in seconds, of the release.

[0065] It will be appreciated that the above descriptions are intended only to serve as examples, and that many other embodiments are possible within the scope of the present invention as defined in the appended claims.