DEVICE FOR ESTIMATING A DIRECTION POINTING TOWARD AN IRRADIATING SOURCE
20250102610 · 2025-03-27
Assignee
Inventors
- Yoann MOLINE (Gif-Sur-Yvette, FR)
- Gwenolé CORRE (Gif-Sur-Yvette, FR)
- Jordan GAMEIRO (Gif-Sur-Yvette, FR)
- Clément LYNDE (Gif-Sur-Yvette, FR)
- Romuald WOO (Gif-Sur-Yvette, FR)
Cpc classification
International classification
G01T1/29
PHYSICS
Abstract
A device for estimating a position of a radioactive source, comprising: a holder; and at least two radiation detectors configured to generate a detection signal, the device being such that, each detector being assigned a position, the detectors are placed around an iso-centroid of each position, i.e., a centroid of each position equally weighted. The device also includes processing circuitry configured to compute a centroid of each position weighted by the detection signal measured by each detector; and depending on the centroid, estimate a direction pointing toward the source, corresponding to a direction between the device and the source.
Claims
1. A device for estimating a direction pointing toward a radioactive source, the radioactive source emitting radiation composed of gamma particles or X-rays or neutrons, the device comprising: a holder; and at least two radiation detectors, each detector being configured to generate a detection signal, depending on a number of particles detected by the detector, wherein each detector is assigned a position, the detectors being placed around a center corresponding to an iso-centroid of each position, a centroid of each position being equally weighted, the device further comprising processing circuitry configured to compute a centroid of each position weighted by the detection signal measured by the detector occupying said position; depending on the centroid, estimate a direction pointing toward the radioactive source, corresponding to a direction between the device and the radioactive source; store the detection signals of the detectors at a first time; store the detection signals of the detectors at a second time, after the first time; take into account a position of each detector at the first time and at the second time; and compute the centroid of the respective positions of the detectors at each time, weighted by the detection signals measured by each detector at each time, respectively.
2. The device as claimed in claim 1, further comprising at least three non-aligned detectors.
3. The device as claimed in claim 1, further comprising an orienting unit configured to estimate a change in the position of the detectors between the first time and the second time.
4. The device as claimed in claim 1, wherein the device further comprises a display, connected to the processing circuitry, and configured to display the direction pointing toward the radioactive source.
5. The device as claimed in claim 1, wherein at least one detector or each detector comprises a detection scintillator.
6. The device as claimed in claim 5, wherein at least one detector or each detector comprises an inorganic detection scintillator.
7. The device as claimed in claim 1, and wherein the processing circuitry is further configured to: take into account a distance of movement of the device between two successive measurement times, in the direction pointing toward the radioactive source; measure the detection signals generated by a detector of the device at the two successive measurement times; and estimate a distance between the device and the source depending on the distance of movement and on the detection signals measured at each time.
8. The device as claimed in claim 7, further comprising a geolocating unit configured to measure the distance of movement of the device between the two successive measurement times.
9. A method for determining a direction pointing toward a radioactive source using the device as claimed in claim 1, the method comprising: a) measuring detection signals with each detector; b) computing a centroid of each position weighted by the signal measured by the detector occupying said position; and c) estimating the direction pointing toward the source depending on the computed centroid.
10. The method as claimed in claim 9, wherein: the step a) comprises: measuring a detection signal resulting from at least two detectors at a first time; and measuring the detection signal resulting from said detectors at a second time; the step b) comprises taking into account a position of each detector at the first time and at the second time; and the step c) comprises computing a centroid of the respective positions of the detectors at each time, weighted by the detection signals measured by each detector at each time, respectively.
11. The method as claimed in claim 10, wherein between the first time and the second time the device rotates about the iso-centroid, the angle of rotation being predetermined and/or measured.
12. The method as claimed in claim 9, further comprising displaying the direction pointing toward the radioactive source.
13. The method as claimed in claim 9, further comprising: d) moving the device between two positions along the direction pointing toward the source, by a distance of movement (d); e) at each position, measuring a detection signal by means of a detector; and f) taking into account the distance of movement and the detection signal measured at each position to estimate a distance between the radioactive source and the device.
14. The method as claimed in claim 13, further comprising: g) comparing the distance between the source and the device, estimated in step f), with a predetermined threshold distance; and h) depending on the comparison, validating or invalidating the estimated distance between the radioactive source and the device.
Description
FIGURES
[0057]
[0058]
[0059]
[0060]
[0061]
[0062]
[0063]
[0064]
[0065]
[0066]
DESCRIPTION OF PARTICULAR EMBODIMENTS
[0067]
[0068] Use of the device presupposes that the radioactive source is placed at a distance such that, in view of its activity, the radiation incident on the device is sufficient to be detected.
[0069] In the examples shown, the device is intended to be borne by a user, so as to indicate the direction pointing toward the source. It is therefore a device intended to be used nomadically. Alternatively, as described below, the device may be mounted on a vehicle, in particular in the context of deployment in hostile environments.
[0070] The device 1 comprises a holder 2 on which a plurality of detectors 10.sub.1, 10.sub.2 are placed. Each detector is configured to detect ionizing radiation, in particular neutrons or X-ray or gamma photons.
[0071] Each detector is formed from a detection material connected to a detecting circuit. Under the effect of an interaction in the detection material, a particle forming the ionizing radiation (neutron or X-ray or gamma photon) generates an electrical pulse, the latter being detected by the detecting circuit. When the detection material is a scintillator, the electrical pulse is generated by a photodetector coupled to the detection material. When the detection material is a semiconductor, the electrical impulse results from collection of charge carriers created in the material under the effect of the interaction. The detecting circuit allows a count rate to be established, which corresponds to a number of pulses detected per unit of time. For example, a number of pulses per second, usually referred to as the counts per second, is established. Alternatively, the detecting circuit allows a count to be established, which corresponds to a number of pulses detected during a given length of measurement time.
[0072] Generally, the detecting circuit allows a detection signal to be formed, the latter corresponding either to a count rate or to a count. In the examples described below, the detection signal is a count rate.
[0073] Preferably, the detection material of each detector is solid. It may in particular be a scintillator. When the particles forming the radiation are photons, inorganic scintillator materials are preferred, due to their high atomic number, which means they attenuate well. Among inorganic scintillators, mention may be made, non-limitingly, of NaI(TI), CsI(TI), GBO (bismuth germanate), LaBr.sub.3, YSO(Ce) (yttrium orthosilicate: ceriumY.sub.2SiO.sub.5:Ce), and Srl.sub.2(Eu). When the particles forming the radiation are neutrons, organic scintillator materials are preferred, because of their high light-atom content.
[0074] In the example illustrated in
[0075] The device also comprises: [0076] a processing unit 20, connected to each detector. The processing unit may comprise a computer, for example a field-programmable gate array (FPGA) or an application specific integrated circuit (ASIC) or a microprocessor, programmed to perform computations using the count rates generated by each respective detector. The processing unit is programmed to establish a direction, pointing toward the source, based on the count rates. [0077] an orienting unit 21, configured to measure an orientation of the device. For example, it may be a question of an inertial measurement unit. The orienting unit is optional. [0078] a display unit 22, for displaying the direction pointing toward the source. It may be a question of a screen or of a set of diodes. The display unit is optional. [0079] a geolocating unit 23, configured to establish a position of the device. For example, it may be a question of a GPS system when the device is used outdoors. It may be a question of a system using beacons emitting a short-range (Wi-Fi for example) electromagnetic signal when the device is used indoors. The geolocating unit is optional.
[0080]
[0081] Based on the position of each detector on the holder 2, an iso-centroid O may be established. The iso-centroid O is a geometrical point located equidistant from each detector. To compute the iso-centroid, each detector 10.sub.i may be assigned a discrete position P.sub.i, which corresponds to the center of the detector. i is a natural number designating each detector.
[0082] The processing unit 20 is configured to compute a centroid B of the position P.sub.i of each detector 10.sub.i, weighted by the count rate measured by each detector, respectively. If P.sub.i designates a point on the holder associated with each detector 10.sub.i, which may for example be the center of each detector, the centroid B is such that:
N.sub.i corresponds to the count rate measured by the detector 10.sub.i.
[0083] The vector {right arrow over (OB)} points toward the source S. Thus, computation of the vector {right arrow over (OB)} makes it possible to estimate a direction, pointing toward the source, from the device. The direction pointing toward the source may be displayed on the display unit 22. The device then operates in a manner equivalent to a compass. The display unit may comprise an arrow pointing toward the source.
[0084] Using only two detectors, a single measurement is not enough, because the direction obtained is necessarily parallel to the direction in which the detectors are aligned.
[0085] However, by rotating the device about an axis perpendicular to the holder and passing through the iso-centroid O, the vector {right arrow over (OB)} may be computed a number of times and the orientation of the device for which the norm {right arrow over (OB)} is maximum may be determined. In such an orientation, the detectors are aligned in the direction pointing toward the source.
[0086] When the detectors 10.sub.1, 10.sub.2 and the source S are aligned, the detector 10.sub.1 forms a front detector and the detector 10.sub.2 forms a rear detector. The front detector is interposed between the source and the rear detector. In this arrangement, the radiation emitted by the source and reaching the rear detector is attenuated by the front detector. Since the detection materials of the detectors are absorbent, in particular when they are inorganic detectors, the attenuation of the radiation by the front detector is significant. Thus, in this configuration, the discrepancy between the respective count rates measured by each detector is maximum. This maximizes the norm {right arrow over (OB)} .
[0087] The display unit 22 may display an arrow the size of which is correlated to the norm {right arrow over (OB)}. The user may thus visually determine the orientation of the device maximizing the norm {right arrow over (OB)}. Alternatively, a numerical value, corresponding to the norm, is displayed on the display unit.
[0088]
[0089] The processing unit 20 is configured to store the respective count rates at each measurement time and to compute a position of the centroid B.
where N.sub.i(t.sub.j) is the count rate of detector 10.sub.i, occupying position P.sub.i(t.sub.j), at measurement time t.sub.j, with j varying between 1 and 2 in this example.
[0090]
[0091] According to one possibility, the device comprises an orienting unit 21 such as described above. In this case, the orienting unit 21 may be configured to detect, during the rotation of the device, when an angle of rotation of 90 with respect to the first time t.sub.1 has been reached, and to carry out the measurements directly.
[0092] According to another possibility, the orienting unit 21 determines an angle of rotation between the first time and the second time. This allows a more precise computation of the position of the centroid B and avoids the need for the user to rotate through an angle close to 90. The user may actuate a switch, so as to signal the first and second times to the processing unit. The orienting unit determines the angle between the first time and the second time and delivers it to the processing unit. The actual value of the angle is taken into account to compute the centroid. The rotation between the two times may thus be different from 90. It may be 45 or 110. However, it is preferable for the angle of rotation to be close to 90.
[0093]
[0094] Whatever the configuration, when the radiation to be detected is gamma radiation, use of detectors based on inorganic materials maximizes absorption of some of the radiation emitted by the source by the detectors closest the source. This increases the contrast between the count rates measured by each detector. The centroid B is thus further away from the iso-centroid O, this increasing the accuracy of the direction estimated by the device. In order for the absorption resulting from one detector to be able to significantly influence the count rate of another detector, the distance between detectors is preferably less than 30 or 40 cm. In order for the effect of absorption to be more significant, it is preferable for the distance between the detectors to be less than 20 cm.
[0095]
[0096] The inventor has carried out an experimental trial in which four detectors were arranged in the same horizontal plane, in a configuration similar to the one shown in
[0101] The difference between the various count rates made it possible to determine, without difficulty, the direction pointing toward the source.
[0102] A configuration such as shown in
TABLE-US-00001 TABLE 1 5 m 4 m 3 m 2 m 1 m 10.sub.3 (left) 187.61 293.15 521.15 1172.59 4690.35 10.sub.1 (front) 205.37 320.90 570.48 1283.58 5134.33 10.sub.4 (right) 190.29 297.33 528.59 1189.33 4757.32 10.sub.2 (rear) 98.72 154.25 274.22 616.99 2467.95
[0103] It may be seen that the contrast observed between the front (10.sub.1) and rear (10.sub.2) detectors is high whatever the distance.
[0104] In another model, a .sup.137Cs source such as described above placed 5 meters from the device was modeled. The material forming each detector was varied. Regardless of the material, each detector was cylindrical in size, and had a diameter and height equal to 25.4 mm. Table 2 shows the count rates estimated for each detector as a function of the detection materials, the latter being NaI(TI), BGO, LaBr.sub.3(Ce), and Srl.sub.2(Eu).
TABLE-US-00002 TABLE 2 Nal(Tl) BGO LaBr.sub.3(Ce) Srl.sub.2(Eu) 10.sub.3 (left) 1958.49 3326.15 2454.65 2247.79 10.sub.1 (front) 2064.44 3501.76 2577.38 2375.98 10.sub.4 (right) 1987.16 3362.26 2486.10 2281.82 10.sub.2 (rear) 963.66 568.30 909.89 949.59
[0105] It may be seen that the contrast between the front (10.sub.1) and rear (10.sub.2) detectors is particularly high for BGO or, to a lesser extent, LaBr.sub.3(Ce), because of their high atomic number, which is conducive to a strong attenuation of the photons.
[0106]
[0107] Operation of the device according to the embodiment described with reference to
[0111] The embodiments described above make it possible to determine a direction pointing toward a radioactive source. The directional information may then be used to estimate a distance between the device and the source.
[0112] Assuming a point source, the intensity of the radiation decreases proportionally to the inverse of the square of the distance. If r designates a distance between a radioactive source and the device, the intensity of the radiation, between the device and the source, varies as 1/r.sup.2. The intensity may correspond to a flux (number of particles per unit time) or to a fluence rate (number of particles per unit time and area). The count rate of each detector varies linearly as a function of the intensity of the radiation emitted by the source. Knowing the direction in which the source is located, by implementing one of the embodiments described above, it is possible to estimate a distance between the device and the source, via movement, by a known distance, with respect to the source, and via measurement of the variation in the count rate of at least one detector. The movement with respect to the source is made in the direction pointing toward the source. Knowing the distance traveled, combined with measurement of the count rates before and after the movement, makes it possible to evaluate the distance between the device and the source.
[0113] Thus, when the device is moved, in a rectilinear direction, pointing toward the source, between a position r.sub.1 and a position r.sub.2, see
where N.sub.1 and N.sub.2 are the count rates measured, by a detector (and preferably the same detector) when the device occupies the positions r.sub.1 and r.sub.2 with respect to the source, respectively. If d is a distance between r.sub.1 and r.sub.2, then r.sub.2=r.sub.1+d.
[0114] d is a signed distance, which is negative if r.sub.2<r.sub.1 and positive if r.sub.2>r.sub.1.
[0115] Development of (3) leads to:
N.sub.2 and N.sub.1 are measured. d may be predetermined or measured. When d is predetermined, the user bearing the device moves, along the direction pointing toward the source, by a predetermined distance, for example 50 or 60 cm, this corresponding to one increment. Alternatively, the device is equipped with a geolocating unit 23, allowing the distance d traveled between the positions r.sub.2 and r.sub.1 to be measured.
[0116] The processing unit 20 is programmed to estimate r.sub.1 by solving (4).
[0117] Count rates follow a Poisson distribution. Expression (4) may be employed provided that the difference between N.sub.1 and N.sub.2 is greater than a threshold, so as to be representative from a statistical point of view. Assuming that N.sub.1 and N.sub.2 are high enough to follow a Gaussian distribution, and a risk of false positive of 2.5%, the difference between N.sub.1 and N.sub.2 is considered statistically significant when
with 1=1.96.
[0118] Beyond the statistical threshold (1){square root over (2N.sub.2)}, the difference between N.sub.1 and N.sub.2 has a 97.5% probability of being statistically representative and of not being due to statistical fluctuations.
[0119] A feasibility trial has been carried out using a BGO inorganic scintillator detector of 25.4 mm25.4 mm (1 inch by 1 inch). The scintillator was moved, in increments of 20 cm, in a straight line, with respect to a point .sup.137Cs source. The .sup.137Cs source generated an equivalent dose rate of 10 Sv.h.sup.1, and was placed 5 meters from the detector. The distance between the detector and the source varied between 10 m and 2 m.
[0120] The relative difference
was also computed, the scale of the latter being plotted on the right-hand y-axisunit %. Curve b) corresponds to the variation in the relative difference as a function of distance.
[0121] At each measurement point, the statistical threshold as defined in (5) was computed. Curve c) shows the variation in the statistical threshold as a function of distance.
[0122] The difference between the count rates N.sub.2N.sub.1 may be considered to be statistically representative when N.sub.2N.sub.1>(1){square root over (2N.sub.2)}. In this example, N.sub.2N.sub.1 crosses the statistical threshold at a threshold distance d.sub.s of 6.5 m from the source. The differential measurement N.sub.2N.sub.1 may thus be considered to allow an estimation of the distance between the detector and source over a distance range of between 0 and 6.5 m, i.e. up to the threshold distance. This corresponds to a relative difference of about 6%.
[0123]
(grayscale levels) as a function of the variation in distance r.sub.1r.sub.2 (x-axisunit meterr.sub.2<r.sub.1), for various distances r.sub.1 with respect to the source (y-axisunit meter). The more r.sub.1 increases, the greater the distance variation r.sub.1r.sub.2 must be to observe a relative difference greater than statistical fluctuations. This figure shows that the distance r.sub.1r.sub.2 is preferably greater than 20 cm, to obtain a relative difference greater than 6%, the latter corresponding to the statistical threshold: see
[0124]
[0125] To determine the distance between the source and the device, the detector closest to the source will preferably be used, since said detector is not subject to attenuation of the radiation emitted by the source by the other detectors.
[0126]
[0127] Step 100: measuring the count rates of the detectors forming the device.
[0128] When the device comprises more than 2 detectors, step 120 may be implemented directly.
[0129] Step 110: this step is implemented in particular when the device comprises only 2 detectors. The device is rotated around iso-centroid O and perpendicular to the holder. The count rate of each detector is measured.
[0130] Step 120: Depending on the measurements performed in step 100 and the optional step 110, computing the centroid of the positions of each detector weighted by the count rates measured by each respective detector.
[0131] Step 130: Depending on the position of the centroid resulting from step 120, determining a direction pointing toward the source.
[0132] Step 140: moving the device in the direction resulting from step 130 by a predetermined distance or any distance.
[0133] Step 150: when the distance of movement of step 140 is arbitrarily chosen, estimating the distance by means of the geolocating unit or by means of the inertial measurement unit, beginning from a starting position.
[0134] Step 160: measuring the count rate of a detector of the device.
[0135] Step 170: based on the measurement resulting from step 160 and on the measurement made by the detector before the movement, for example in step 100 (or 110), estimating a distance between the source and the device.
[0136] Step 180: comparing the estimated distance with the threshold distance d.sub.s, so as to validate or invalidate the measurement. The threshold distance was determined beforehand.
[0137] The device described above may be used while being carried by a user. It may also be mounted on an autonomous vehicle. In this case, the direction determined by the device is transmitted to a control unit of the vehicle, so that the vehicle may be steered depending on this direction. Use of a vehicle is recommendable when operating in hostile environments, for example in the presence of contamination or high levels of radiation.