Brevet FR2947344A1 — publié le 20101231

PROCÉDÉ DE DÉTECTION DE RAYONNEMENT GAMMA ET/OU NEUTRONIQUE ET DISPOSITIF ASSOCIE L'invention concerne un procédé de détection de rayonnement gamma et/ou neutronique et le dispositif associé. ARRIERE PLAN DE L'INVENTION Dans le domaine nucléaire, le besoin existe de discriminer les parties de rayonnement provenant des neutrons des parties de rayonnement provenant des rayons gamma afin d'identifier la menace NR (NR pour � Nucléaire et Radiologique �), par exemple pour évaluer la gravité d'une menace ou, plus concrètement, pour évaluer une masse de matière fissile. A ce jour, la discrimination à partir de signaux issus de scintillateurs organiques est réalisée par analyse de forme. Deux procédés connus sont utilisés qui sont basés sur un principe consistant à exploiter une différence de population des états énergétiques de molécules selon le pouvoir ionisant des particules incidentes, lequel est lié à leur masse et à leur charge. Un premier procédé consiste à dupliquer un signal issu du scintillateur sur deux lignes. Le signal de la première ligne est intégré puis différencié. Le passage du signal par zéro est caractéristique d'une particule ayant interagi avec le matériau fluorescent du scintillateur. Ensuite, la distance temporelle entre le début dudit signal et le passage par zéro est déterminée. Sur la seconde ligne, le début du signal est déterminé par un discriminateur à fraction constante. Le temps ainsi déterminé est directement lié à la particule qui a interagi avec le scintillateur. Plus ce temps est important, plus la particule est ionisante. Le second procédé consiste à intégrer un signal issu d'un scintillateur sur deux fenêtres temporelles différentes, la première contenant l'ensemble du signal, et la seconde uniquement la partie dite retardée, i.e. celle contenue dans la seconde partie du signal, la partie décroissante. Par analyse biparamétrique, on peut alors distinguer deux parties dans le graphique bidimensionnel affichant la charge totale versus la charge retardée, chacune correspondant à une zone d'intérêt contenant les impulsions soit d'origine neutronique, que soit d'origine gamma. Ces deux procédés comparent deux signaux entre eux. Un inconvénient lié à une comparaison entre deux signaux est que, compte tenu de la variabilité des formes, il est très difficile de discriminer les particules entre elles dans les cas où la différence des signaux est très faible, ce qui est le cas des scintillateurs organiques solides. OBJET DE L'INVENTION L'invention a pour objet un procédé de discrimination de particules adapté pour être utilisé avec des scintillateurs organiques solides n'ayant pas les inconvénients de l'art antérieur. BREVE DESCRIPTION DE L'INVENTION Dans ce but, l'invention propose un procédé de détection de rayonnement gamma et/ou neutronique dans lequel le rayonnement gamma et/ou neutronique est incident sur un scintillateur qui délivre un signal lumineux qui est transmis à un photomultiplicateur, lequel photomultiplicateur effectue une conversion lumière/électrons et délivre un signal électronique qui comprend une composante gamma et/ou une composante neutronique, le procédé de détection comprenant les étapes suivantes mises en oeuvre par un calculateur : ù retarder d'un temps TAU et atténuer d'un coefficient ALPHA le signal électronique pour obtenir un signal retardé et atténué, ù soustraire le signal retardé et atténué du signal électronique pour obtenir un signal de différence qui comprend une composante gamma de différence et/ou une composante neutronique de différence, et ù calculer une première grandeur sigmal et une seconde grandeur sigma2 telles que : ô sigmal = f S3(t) dt Tl Tz sigma2 = f S3(t) dt où les grandeurs T1 et T2 sont des instants qui encadrent un passage par zéro d'une composante gamma de différence et ô est un instant de passage par zéro de la composante gamma de différence, et ù calculer les grandeurs : y=sigmal/A1, et x=sigma2/A2, où Al est l'amplitude du signal de différence pendant la durée 5-T1 et A2 l'amplitude du signal de différence pendant la durée T2-5. BREVE DESCRIPTION DES FIGURES D'autres caractéristiques et avantages de l'invention apparaîtront à la lumière d'un mode de réalisation préférentiel fait en référence aux figures jointes, parmi lesquelles : û La Figure 1 représente un schéma de principe de système de détection de l'invention ; û La figure 2 représente un schéma détaillé d'un circuit qui appartient au système représenté en figure 1 ; û Les figures 3a-3d et 4a-4d représentent différents signaux gamma et neutroniques aptes à décrire différentes étapes du procédé de discrimination de l'invention ; et û La figure 5 représente, à titre d'exemple, l'affichage d'un résultat obtenu par le procédé de discrimination de l'invention. Sur toutes les figures les mêmes références désignent les mêmes éléments. EXPOSE DETAILLE DU MODE DE REALISATION PREFERENTIEL DE L'INVENTION La Figure 1 représente un schéma de principe d'un système de détection de l'invention. Le système comprend un détecteur 1, par exemple un scintillateur organique solide, adapté pour recevoir des rayonnements gamma et/ou des rayonnements neutroniques. La sortie du détecteur 1 est reliée à l'entrée d'un photomultiplicateur 2 polarisé par une haute tension HT. Le photomultiplicateur 2 a une sortie reliée à l'entrée d'un convertisseur analogique/numérique 3 (par exemple un convertisseur 1GHz, 8 bits), dont la sortie est reliée à l'entrée d'un calculateur 4. Le calculateur 4 est décrit en détail en référence à la figure 2. La sortie du calculateur 4 est reliée à un dispositif d'affichage 5 qui affiche le résultat des calculs effectués par le calculateur. Selon le mode de réalisation représenté en figure 1, le système de détection comprend un convertisseur analogique/numérique 3. Il faut cependant noter que l'invention concerne également le cas où le système de détection ne comprend pas de convertisseur analogique/numérique. Le traitement des signaux est alors un traitement analogique. En fonctionnement, des rayonnements gamma et/ou neutroniques interagissent avec le scintillateur organique solide 1, lequel délivre un signal lumineux qui est transmis au photomultiplicateur 2. Le photomultiplicateur 2 effectue une conversion lumière/électrons et délivre un signal électronique qui est transmis au convertisseur analogique/numérique 3 qui convertit les signaux analogiques en signaux numériques. Les signaux numériques sont ensuite envoyés au calculateur 4. Le résultat du calculateur 4 est ensuite affiché sur le dispositif d'affichage 5. La Figure 2 représente une vue détaillée du calculateur 4. Le circuit d'entrée du calculateur 4 comprend préférentiellement un préamplificateur 6. Ce préamplificateur peut toutefois être supprimé dès lors que le système de détection comprend un photomultiplicateur 2 de performances suffisantes pour induire un signal exploitable. La sortie du préamplificateur 6 est reliée, d'une part, à l'entrée d'un dispositif de retard 7 et, d'autre part, à une première entrée d'un opérateur de différence 9. La sortie du dispositif de retard est reliée à l'entrée d'un dispositif d'atténuation 8 dont la sortie est reliée à une seconde entrée de l'opérateur de différence 9. La sortie de l'opérateur de différence 9 est reliée à l'entrée d'un opérateur d'extraction 10 dont une première sortie est reliée à une première entrée d'un dispositif d'acquisition 11 et une seconde sortie est reliée à une seconde entrée du dispositif d'acquisition 11. La sortie du dispositif d'acquisition 11 est reliée à l'entrée du dispositif d'affichage 5. En fonctionnement, un signal SO est fourni en entrée du préamplificateur 6. Le préamplificateur 6 amplifie le signal SO et délivre un signal S1. Le signal S1 comprend une composante gamma et/ou une composante neutronique, la composante neutronique se décomposant elle-même en deux composantes qui se superposent, à savoir une composante rapide due à une désexcitation prompte dans le matériau du scintillateur et une composante retardée due à une désexcitation retardée dans le matériau du scintillateur. Les figures 3a et 3b représentent, respectivement, une composante gamma S1(y) et une composante neutronique S1(n) du signal S1 dans le cas de signaux analogiques. La composante neutronique S1(n) comprend la composante neutronique rapide Sa et la composante neutronique retardée Sb. Le signal S1 qui est délivré par le préamplificateur 6 est alors transmis directement sur la première entrée de l'opérateur différence 9 et, via le dispositif de retard 7 et le dispositif d'atténuation 8, sur la deuxième entrée de l'opérateur 9. Le dispositif de retard 7 retarde le signal d'un temps TAU et le dispositif d'atténuation 8 atténue le signal d'un coefficient ALPHA. Un signal S2 est délivré par le dispositif de retard 8. Les figures 3c et 3d représentent, respectivement, la composante gamma retardée et atténuée S2(y) et la composante neutronique retardée et atténuée S2(n) du signal S2. Les valeurs TAU et ALPHA sont choisies à l'aide d'un processus d'itération. TAU a de préférence une valeur comprise entre 0 et 1Ons et ALPHA a de préférence une valeur comprise entre 0 et 1. Le processus d'itération est effectué au préalable en faisant passer un même signal de référence plusieurs fois dans le dispositif de retard 7 et le dispositif d'atténuation 8. Lors du processus d'itération, plusieurs sources radioactives présentant individuellement ou collectivement des flux gamma et/ou neutron peuvent être utilisées pour calibrer le système. Les valeurs de TAU et de ALPHA qui sont choisies sont celles pour lesquelles les deux lobes visualisés en sortie du dispositif d'acquisition qui représentent, l'un un rayonnement gamma, l'autre un rayonnement neutronique, sont les plus séparés. L'opérateur de différence 9 effectue la différence entre le signal S1 et le signal S2. Le signal S3 délivré en sortie de l'opérateur 9 est alors . S3 = S1 - S2 Le signal S3 comprend une composante gamma S3(y) et/ou une composante neutronique S3(n). Les figures 4a et 4b représentent, respectivement, la composante gamma S3(y) et la composante neutronique S3(n) en fonction du temps. On constate que le signal S3(y) passe par zéro, alors que la courbe S3(n) ne passe jamais par zéro. L'opérateur d'extraction 10 qui reçoit le signal S3 est programmé pour être utilisé à partir d'un instant T1 jusqu'à un instant T2, un paramètre temporel ô situé entre les instants T1 et T2 indiquant l'instant de passage par zéro du signal S3. Le paramètre ô est réglé au préalable, lors d'une étape d'étalonnage, pour chaque détecteur, à l'aide d'un signal d'étalonnage. Une fois réglé, ce paramètre ô ne dépend que de la mise en oeuvre du système de détection et, éventuellement, du vieillissement du matériau scintillateur et/ou du vieillissement du photomultiplicateur (diminution du rendement de fluorescence et/ou du gain du photomultiplicateur). L'opérateur d'extraction 10 calcule des grandeurs sigmal et sigma2 telles que : sigmal = J S3(t) dt Tl Tz sigma2 = f S3(t) dt30 Les figures 4c et 4d illustrent de façon symbolique les quantités sigmal et sigma2 pour chacun des signaux S3 (y) et S3 (n) . Le dispositif d'acquisition 11 reçoit en entrée les grandeurs calculées sigmal et sigma2 et calcule les grandeurs x et y telles que : Y = sigmal/A1, et x = sigma2/A2 où Al est l'amplitude du signal de différence pendant la durée 5-T1 et A2 l'amplitude du signal de différence pendant la durée T2-5. Le résultat du traitement effectué par le dispositif d'acquisition 11 est affiché par le dispositif d'affichage 5. La figure 5 illustre, à titre d'exemple, l'affichage d'un résultat. Avantageusement, selon l'invention, le fait de corréler le signal reçu avec lui-même permet de prendre en compte toute la variabilité interne des impulsions qui constituent le signal. Il est alors possible de calculer une grandeur (sigma2) apte à discriminer un signal gamma d'un signal neutronique. La grandeur sigma2 est en effet négative pour un signal gamma et positive pour un signal neutronique. En affichant sur un même diagramme le couple de valeurs (x,y) pour chaque incidence d'un rayonnement reçu, les valeurs représentant un rayonnement gamma peuvent alors être séparées des valeurs représentant un rayonnement neutronique, les valeurs représentant un rayonnement gamma étant sur la partie gauche du diagramme et les valeurs représentant un rayonnement neutronique étant sur la partie droite. De ce fait, une zone R1 située principalement du côté des x négatifs rassemble l'ensemble des incidences des rayonnements gamma et une zone R2 située principalement du côté des x positifs rassemble l'ensemble des incidences des rayonnements neutroniques (voir figure 5). Les zones R1 et R2 peuvent toutefois se recouvrir partiellement du fait du phénomène connu de diaphonie. Compte tenu de la nature des rayonnements discriminés, le procédé de l'invention ne peut conduire à l'établissement de résultats fiables qu'après intégration d'un grand nombre d'impulsions. En effet, bien qu'une analyse soit effectuée impulsion par impulsion, un résultat pertinent ne peut être considéré que d'un point de vue global, lorsqu'une analyse en nombre a été conduite. Typiquement il faut une analyse comprise entre 100 et 1000 impulsions unitaires (gamma et/ou neutrons) pour avoir un résultat pertinent permettant l'analyse de la nature de la source de rayonnement observée.20 REVENDICATIONS1. Procédé de détection de rayonnement gamma et/ou neutronique dans lequel le rayonnement gamma et/ou neutronique est incident sur un scintillateur (1) qui délivre un signal lumineux qui est transmis à un photomultiplicateur (2), lequel photomultiplicateur effectue une conversion lumière/électrons et délivre un signal électronique (Si) qui comprend une composante gamma (S1(y)) et/ ou une composante neutronique (S1(n)), caractérisé en ce qu'il comprend les étapes suivantes mises en oeuvre par un calculateur : ù retarder (7) d'un temps TAU et atténuer (8) d'un 15 coefficient ALPHA le signal électronique (Si) pour obtenir un signal retardé et atténué (S2), soustraire (9) le signal retardé et atténué (S2) du signal électronique (Si) pour obtenir un signal de différence (S3) qui comprend une composante gamma de 20 différence (S3(y)) et/ou une composante neutronique de différence (S3(n)), et ù calculer (10) une première grandeur sigmal et une deuxième grandeur sigma2 telles que : sigmal = J S3(t) dt Tl Tz 25 sigma2 = f S3(t) dt où les grandeurs Ti et T2 sont des instants qui encadrent un passage par zéro d'une composante gamma de différence (S3(y)) et ô un instant de passage par zéro de la composante gamma de différence (S3(y)), 30 etcalculer des grandeurs y et x telles que : y=sigmal/A1, et x=sigma2/A2, où Al est l'amplitude du signal de différence pendant la durée 5-T1 et A2 l'amplitude du signal de différence pendant la durée T2-5. 2. Procédé selon la revendication 1, dans lequel les valeurs de TAU, de ALPHA et de ô sont déterminées lors d'une étape d'étalonnage qui précède la détection. 3. Procédé selon la revendication 1 ou 2, dans lequel TAU a une valeur comprise entre 0 et 1Ons et ALPHA une valeur comprise entre 0 et 1. 4. Procédé selon l'une quelconque des revendications précédentes, dans lequel le scintillateur est un scintillateur organique solide. 5. Procédé selon l'une quelconque des revendications précédentes, dans lequel les grandeurs Y et x sont transmises à un dispositif d'affichage (5). 25 6. Système de détection de rayonnement gamma et/ou neutronique dans lequel le rayonnement gamma et/ou neutronique est fourni à un scintillateur (1) qui délivre un signal lumineux qui est transmis à un photomultiplicateur (2), lequel photomultiplicateur 30 effectue une conversion lumière/électrons et délivre un signal électronique (Si) qui comprend une composante20gamma (S1(y) et/ ou une composante neutronique (S1(n)), caractérisé en ce qu'il comprend un calculateur qui comprend : ù un dispositif de retard (7) et un dispositif d'atténuation (8) pour, respectivement, retarder d'un temps TAU et atténuer d'un coefficient ALPHA le signal électronique (Si) afin d'obtenir un signal retardé et atténué (S2), ù un opérateur de différence (9) pour soustraire le signal retardé et atténué (S2) du signal électronique (Si) afin d'obtenir un signal de différence (S3) qui comprend une composante gamma de différence (S3(y)) et/ou une composante neutronique de différence (S3 (n)) , et ù une unité de calcul (10) pour calculer une première grandeur sigmal et une deuxième grandeur sigma2 telles que : sigmal = J S3(t) dt Tl Tz sigma2 = f S3(t) dt ù où les grandeurs Ti et T2 sont des instants qui encadrent un passage par zéro d'une composante gamma de différence (S3(y)) et ô un instant de passage par zéro de la composante gamma de différence (S3(y)), et ù une unité de calcul (11) pour calculer des grandeurs y et x telles que : y=sigmal/A1, et x=sigma2/A2, oùAl est l'amplitude du signal de différence pendant la durée 5-T1 et A2 l'amplitude du signal de différence pendant la durée T2-5. 7. Système selon la revendication 6, dans lequel les valeurs de TAU, de ALPHA et de ô sont déterminées lors d'une étape d'étalonnage qui précède la détection. 8. Système selon la revendication 6 ou 7, dans lequel TAU a une valeur comprise entre 0 et 1Ons et ALPHA une valeur comprise entre 0 et 1. 9. Système selon l'une quelconque des revendications 6 à 8, dans lequel le scintillateur est un scintillateur organique solide. 10. Système selon l'une quelconque des revendications 6 à 9, dans lequel un dispositif d'affichage (5) affiche les grandeurs x et y.