Brevet FR2652651A1 — publié le 19910405

La présente invention concerne un système de détection de substances et en particulier d'expLosifs, par irradiation neutronique de ceux-ci. Elle s'applique notamment au contrôle de bagages, en particuLier dans les aéroports. On connaît déjà des systèmes de détection d'explosifs, utiLisant une irradiation de ceux-ci par des neutrons. On peut consulter notamment à ce sujet les documents FR-A-2201765 et EP-A-0227497. On rappelle à ce propos les notions suivantes : IL existe trois processus de production de photons gamma par interaction de neutrons avec des noyaux : - La capture d'un neutron par un noyau, engendrant un isotope qui se désexcite par émission d'un rayonnement gamma, appelé rayonnement gamma de capture (capture gamma radiation dans Les publications en Langue anglaise), cette interaction ayant essentiellement lieu avec des neutrons thermiques ("thermaL neutrons1 dans Les publications en langue angLaise), - La diffusion inélastique d'un neutron sur un noyau qui se désexcite par émission d'un rayonnement gamma appelé rayonnement gamma prompt ("prompt gamma radiation" dans Les publications en Langue angLaise), cette interaction existant seulement avec des neutrons rapides ("fast neutrons" dans Les publications en tangue angLaise) ayant assez d'énergie (c'est-à-dire une énergie au moins égaLe à celle du rayonnement gamma prompt), et - l'activation d'un noyau par un neutron thermique ou rapide, qui crée un noyau radioactif ayant une certaine durée de vie et se désintégrant en émettant un rayonnement gamma appelé rayonnement gamma d'activation ("activation gamma radiation" dans les publications en Langue anglaise). Par ailleurs, les explosifs usuels comportent des atomes de carbone, d'hydrogène, d'oxygène et d'azote et tes éLéments Les plus intéressants pour détecter La présence d'un explosif sont l'azote et l'oxygène, notamment à cause de La séparabilité des raies (ou pics) gammas, créés par diffusion inélastique de neutrons rapides, l'énergie des photons gamma se répartissant entre 2,3129 MeV et 6,1304 MeV. L'utilisation du rayonnement gamma de capture de l'azote à 10,83 MeV est également intéressante du fait du bon rapport signal/bruit dans L'intervalle d'énergie proche de 10 MeV. Les systèmes de détection connus1 mentionnés pLus haut, présentent l'inconvénient de donner Lieu à une probabilité trop élevée de défaut de détection (de -2 l'ordre de 10 ou plus) lors de contrôles à cadences élevées et pour de faibles quantités d'explosif. La présente invention a pour but de remédier à cet inconvénient en proposant un système de détection de substances et en particulier d'explosifs, qui permet d'abaisser cette probabilité de défaut de détection et même d'atteindre des valeurs très faibLes pour celle-ci -4 -10 -10 (10 , 10 ou même moins de 10 ) en utilisant des moyens de détection de haute résolution pour détecter le rayonnement gamma De façon précise, La présente invention a pour objet un système de détection d'une substance susceptible d'être contenue dans un objet, système caractérisé en ce qu' iL comprend - des moyens d'irradiation de l'objet par des neutrons (thermiques et/ou rapides), - des moyens de détection du rayonnement gamma susceptible d'être alors émis par L'objet, et - des moyens éLectroniques de traitement des signaux fournis par Les moyens de détection, ces moyens électroniques de traitement étant prévus pour compter les photons gamma correspondant à chaque raie d'un ensemble de raies caractéristiques d'au moins un élément chimique de la substance, déterminer, pour chaque raie i, une probabilité de fausse détection de l'élément chimique associe à cette raie, c'est-à-dire la probabilité PFi pour que le signal détecté, correspondant à cette raie i, soit dû à un bruit de fond, déterminer le produit de ces probabilités de fausse détection1 comparer ce produit à un seuil fixé par Les utilisateurs du système, et avertir ces utilisateurs si Le produit est inferieur au seuil fixé par eux, l'objet étant alors présumé contenir la substance. Des éléments importants de L'invention, permettant de diminuer ta probabilité de défaut de détection par rapport aux systèmes connus, sont donc : - Le choix d'un ou d'une pLuralité d'éLéments chimiques que contient La substance que l'on veut détecter, par exemple Azote seul ou l'azote et l'oxygène ou encore l'azote, l'oxygène et le carbone dans le cas d'un explosif, - le choix d'une pluralité de raies qui sont caractéristiques du ou des éléments choisis et dont on compte tes photons gamma correspondants, par exemple, dans le cas d'un explosif, Les raies à 5,106 MeV et à 10,83 MeV de L'azote ou bien L'une et/ou L'autre des raies précédentes ainsi qu'une ou plusieurs raies caractéristiques de L'oxygène et, éventuellement, en plus, des raies respectivement caractéristiques d'éléments tels que le carbone, L'hydrogène, Le chlore (explosif chloré), et - L'évaluation de la probabilité de fausse détection relative à chacune de ces raies (ou d'une quantité d'information liée à cette probabilité). Chaque probabilité de fausse détection PFi peut être déterminée par la formule : où Ni représente le comptage, pendant un temps Dt, qui correspond à la raie i et qui est dû à des réactions nucléaires induites par des neutrons et aussi à un bruit de fond M'2i relatif à la raie i, pendant Le temps Dt. Chaque bruit de fond M'2i peut être déterminé en faisant la somme d'un bruit de fond physique M2i, reLatif à la raie i, et de La borne supérieure Ci de comptages relatifs à cette raie i, pendant le temps Dt, sur des objets qui sont susceptibles de contenir La substance mais dont on est sûr qu'ils n'en contiennent pas. SeLon un premier mode de réalisation particulier du système objet de l'invention, Les moyens électroniques de traitement sont prévus pour : déterminer, pour chaque raie i, une quantité dtinformation li relative à cette raie i et définie par ta formule PFi li = - log - PFi 1 - PFi (après avoir déterminé Ni puis PFi), déterminer la somme de ces quantités d'information, comparer cette somme à un seuil fixé par les utilisateurs du système, et avertir ces utilisateurs si cette somme est supérieure au seuil auquel elLe est comparée, l'objet étant alors présumé contenir la substance. Selon un deuxième mode de réalisation particulier, permettant un temps de calcul pLus court et donc un temps de contrôle plus court de L'objet, Les moyens électroniques de traitement sont prévus pour déterminer, pour chaque raie i, une quantité d'information Ii reLative à cette raie i et définie par La formule (après avoir déterminé Ni mais sans avoir à déterminer PFi), formule où K et a sont des fonctions mémorisées d'un bruit de fond M'2i relatif à La raie i et Nîi est la partie entière de la différence Ni-M'2i, Ni représentant le comptage, pendant un temps Dt, qui correspond à La raie i et qui est dû à des réactions nucléaires induites par des neutrons et aussi au bruit de fond M'2i relatif à la raie i, pendant Le temps Dt, déterminer La somme de ces quantités d'information, comparer cette somme à un seuiL fixé par les utiLisateurs du système, et avertir ces utilisateurs si cette somme est superieure au seuil auquel elLe est comparée, L'objet étant alors présumé contenir la substance. Le système objet de L'invention s'appLique aussi bien à une configuration "fermée" qu'à une configuration "ouverte" : Les moyens d'irradiation peuvent comprendre une source de neutrons rapides et une enceinte de thereaLisation de ces neutrons rapides, dans laquelle se trouve la source et qui est prévue pour recevoir L'objet, ou peuvent comprendre une source de neutrons rapides qui est prévue pour irradier L'objet (directement, sans enceinte de thermalisation). Afin de pouvoir contrôler L'ensemble de L'objet, Les moyens de détection peuvent comprendre une pluralité de détecteurs de rayonneaent gamma, tes moyens électroniques de traitement comprenant une pluralité de chaînes de détection, respectivement associées aux détecteurs, la configuration de ces détecteurs pouvant être iatricielle. Selon un mode de réalisation préféré du système objet de L'invention, permettant de diminuer fortement La probabilité de défaut de détection par rapport aux systèmes connus, Les moyens de détection du rayonnement gamma sont des moyens de détection de haute résolution. A cet effet, on peut par exemple utiLiser un ou une pLuraLité de détecteurs en germanium de haute pureté. De préférence, les moyens d'irradiation comprenant une source de neutrons, les moyens de détection sont protégés du rayonnesent direct de la source et collimatés vers L'objet. Dans une réalisation préférée de L'invention, les moyens d'irradiation comprennent une source pulsée de neutrons, prévue pour fournir des bouffées de neutrons, et les moyens électroniques de traitement coopérent avec Les moyens de détection pour faire des mesures dans des intervalles de temps où l'on est sûr de détecter seulement L'une des catégories de photons gamma engendrés lors de l'irradiation de L'objet par des neutrons. Avec cette source de neutrons pulsée, on réalise par exemple une détection synchrone des photons gamma prompts et une détection des photons gamma de capture entre Les bouffées de neutrons. Ceci permet d'obtenir, pour chaque type de détection, un rapport signaL/bruit nettement supérieur à celui que L'on obtient avec une source continue de neutrons. Dans La présente invention, les moyens d'irradiation comprennent, de préférence, une source de neutrons de 14 MeV engendrés par des réactions de fusion En effet, une source de neutrons d'une telle énergie permet de contrôler la plupart des éléments chimiques. Enfin, dans un mode de réalisation particuLier de L'invention, approprié à La détection d'explosif, Les moyens électroniques de traitement sont prévus pour compter Les photons gamma d'un ensemble de raies qui sont respectivement caractéristiques de L'azote, de L'oxygène et éventuellement du carbone, dans le cas d'explosifs nitrés, ou d'autres éléments dans Le cas d'autres explosifs. La présente invention sera mieux comprise à La Lecture de la description suivante, d'exemples de réalisation donnés à titre purement indicatif et nullement Limitatif, en référence aux dessins annexés sur lesquels - la figure 1 est un graphique illustrant les variations de La quantité d'information en fonction du taux moyen de signal, à taux de bruit physique constant et pour divers taux de bruit physique, - La figure 2 est un graphique permettant de comparer l'une à l'autre deux expressions de la quantité d'information, - la figure 3 est un graphique montrant la décroissance monotone, en fonction du taux de bruit physique, de chacune de ces deux expressions de la quantité d'information, - la figure 4 est une vue en perspective schéiatique d'un mode de réalisation particulier du système objet de l'invention, utilisant une pluralité de détecteurs fixes de rayonnement gamma et une enceinte de thermalisation, - la figure 5 illustre schématiquement des moyens électroniques de traitement qui sont utilisables pour traiter les impulsions fournies par ces détecteurs, - La figure 6 illustre schémat;ouement d'autres moyens électroniques de traitement qui sont également utilisables pour traiter les impulsions fournies par les détecteurs représentés -sur la figure 4, - la figure 7 est une vue schématique et partielle d'un autre mode de réalisation particulier du système objet de L'invention, utilisant un détecteur mobile de rayonnement gamma, et - ta figure 8 est une vue schématique d'un autre mode de réalisation particulier du système objet de l'invention, de configuration "ouverte". Dans ce qui suit, on suppose que l'on veut contrôler des bagages qui sont susceptibles de contenir des expLosifs, par exempLe des explosifs nitrés. On ne peut disposer par bagage que d'une durée de temps Limitée afin de pouvoir assurer un débit suffisant, mesuré en nombre de bagages par heure. Cette durée de temps est par exemple de l'ordre de quelques secondes. Un bagage est donc amené à un système conforme à L'invention où it est irradié par des neutrons et l'on effectue un comptage, pendant une durée Dt, de photons gamma parmi lesquels figurent, bien entendu, des photons gamma émis par L'objet à La suite de L'irradiation de ce dernier. A cette durée Dt correspond un nombre moyen d'événements M2 qui est dû à diverses causes physiques (mouvement propre du ou des détecteurs de rayonnement gamma, rayonnement du milieu environnant, rayonnement cosmique ...) et que l'on appelle "bruit de fond physique - Pour cette même durée Dt, il existe un comptage moyen Ml, que l'on peut attribuer à la présence d'un produit suspect, de la matière azotée par exemple. Des essais effectués sur un grand nombre de bagages montrent qu'au bagage typique du voyageurimoyen' correspond un comptage maximum Ci sous un pic i, comptage maximum qui est attribuable par exemple à de la matière azotée (Le pic i étant, par exemple, Le pic à 10,83 MeV de L'azote) sans que ce bagage contienne une quelconque quantité d'explosif. On attribue alors à la présence d'explosif tout comptage qui est supérieur à la somme M2i+Ci, M2i étant la valeur du comptage due au bruit de fond physique sous le pic i pendant Le temps Dt. On considère alors, pour La détection d'une quantité anormale d'azote attribuable à un explosif, qu'il existe un bruit de fond M'2i tel que M'2i = M2i + Ci (1) Conformément à La présente invention, on choisît à l'avance d'examiner un certain nombre de raies i, i allant de 1 à n, dans Le spectre de rayonnement gamma et l'on caractérise la présence de chaque raie i par La probabilité PFi pour que le signal détecté soit dû au bruit de fond ou à une quantité de matière suspecte inférieure à la borne supérieure de la quantité normale, en d'autres termes, pour qu'il y ait fausse détection. On fait correspondre au comptage Ni "brut", sous Le pic i, et au comptage moyen de bruit M'2i, une probabilité de fausse détection PFi teLle que De plus, on considère la valeur Nîi égale à La partie entière de la différence entre le comptage brut Ni sous le pic i et M'2i, c'est-à-dire N1i = E(Ni - M'2i) (3) en notant E(x) la partie entière du nombre x. On fait correspondre à La valeur Nti et à la valeur MX2i du bruit de fond, une probabilité de défaut de détection PDi qui caractérise le défaut inévitable du système de ne pas détecter ce qu'il devrait détecter, cette probabilité de défaut de détection étant teLle que On peut montrer que, à bruit de fond M'2i croissant, PDi et PFi tendent vers une mêse valeur, PDi restant supérieur à PFi, La probabilité PDi constituant un élément d'appréciation fondamental pour connaStre l'aptitude du système à ne pas manquer un bagage qui comporterait une charge d'expLosif. Cependant, dans la présente invention, on utilise la probabilité PFi étant donné son importance "opérationnelle" (et Le fait que PDi et PFi tendent vers ladite même valeur). Dans la présente invention, on peut également mesurer L'information obtenue pendant le temps Dt au moyen d'une fonction appelée "quantité d'information" et notée li pour La raie i, cette quantité Ii dépendant de la probabilité de fausse détection PFi obtenue au cours d'un stockage de photons gamma. Chaque quantité d'information satisfait aux conditions suivantes : - additivité des quantités d'information obtenues tors de deux stockages séparés, effectués dans les mêmes conditions, et - nullité de la quantité d'information avant tout stockage. Les probabilités de fausse détection PF' et PF" sur un phénomène particulier, au cours de deux comptages distincts, sont indépendantes, si bien que la probabiLité de fausse détection PF résultant des deux comptages considérés globalement, a pour expression : PF = PF' x PF" (5) On choisit alors l'expression suivante pour la quantité d'informations Ii : PFi li = - log ------- (6) 1-PFi Avant ta mesure, les probabilités de vraie et de fausse détection sont a priori égales à 0,5 et l'expression (6) de La quantité d'information est alors nulle. Elle est également nulle dans le cas où le comptage Ni est égal au comptage moyen de bruit de fond M2i car PFi est alors égal à 0,5. La quantité d'information Ii est positive si le comptage Ni est supérieur à M2i et négative si Ni est inférieur à M2i. IL est à noter qu'on pourrait changer le signe de La quantité d'information en changeant de critère (obtention d'une quantité d'information Ii' sur La non-existence d'un phénomène particulier, avec Ii'= li). Au sujet de cette notion de quantité d'information, on se reportera aux documents suivants : - Article de M.L. RAMBAUT, intitulé "Détection d'un rayonnement corpusculaire et information" et publié dans Nuclear Instruments and Methods in Physics Research A248 (1986) p.471-482 North Holland, Amsterdam. - EP-A-0200596 - EP-A-0202980 - demande de brevet français n08712970 du 18 septembre 1987. Sur la figure 1, on a tracé la courbe des variations de la quantité d'information I en fonction du taux moyen de signal M1 exprimé en nombre d'événements pendant Le temps Dt, et à bruit, noté M2, constant et exprimé en nombre d'événements pendant Le temps Dt, pour plusieurs valeurs de M2 (10, 50, 100, 4 199, 316, 501, 10 et 10 ), la base des logarithmes étant égale à 10 et la quantité d'information étant exprimée en "digits". Le réseau de courbes de la figure 1 montre bien quantitativement l'avantage qu'il y a à placer le ou les détecteurs de rayonnement gamma Là où le bruit M2 est minimal. Le système objet de L'invention peut comprendre une source de neutrons de fission, par exemple une source de Californium 252. Cependant, ce système comprend avantageusenent une source pulsée de neutrons de fusion de 14 MeV ; en effet, si la mesure n'est faite que dans L'intervalle de temps où L'on est sûr de détecter L'une des trois catégories de photons gamma mentionnées pLus haut, on peut réduire proportionnellement la durée de stockage et donc La valeur du bruit de fond. Au lieu d'utiliser La formule (6) pour La quantité d'information Ii, on peut utiliser une formule semi-phénomenologique que L'on peut déduire des formules (21 et (6) et qui est la suivante où K et a sont des fonctions de M'2i. Cette formule (7), qui est fondée, par l'intermédiaire de la formule (2), sur le calcul des probabilités, facilite Le dialogue entre Le système objet de L'invention et les utilisateurs de ce dernier. D'autre part, eLLe est préférable à une formule du genre : En effet, en considérant la figure 2, sur laquelle une quantité d'information I(7) selon la formule (7) et une quantité d'information I(8) selon la formule (8), exprimées en "digits", sont représentées en fonction d'un taux de signal M1, pour un bruit M2 fixé, égaL à 199 (événements dûs au bruit physique pendant te temps Dt), on voit que l'expression de la formule (7 > fournit plus d'information que L'expression de La formule (8). Dans La présente invention, Le choix des moyens de détection du rayonnement gamma est important. On pourrait choisir un ou des détecteurs de faible résolution, utilisant par exemple des cristaux d'iodure de sodium, ou de très faible résolution, utiLisant par exemple des cristaux de germanate de bismuth, de coût relativement faible. Cependant, de tels détecteurs ne fournissent pas une garantie suffisante aux utilisateurs du système, étant donné que Leur emploi conduit à une probabilité de défaut de détection d'explosif PD trop élevée. C1est pourquoi on utilise de préférence, dans la présente invention, un ou des détecteurs de haute résolution, comportant par exemple un cristal de geriianiua intrinsèque. De tels détecteurs permettent de séparer tous les pics de photons gamma utiles (dont certains ne seraient même pas visibLes sur un spectre obtenu avec un détecteur de faible résolution) et de disposer sous un pic gamma, d'un rapport signal/bruit maximum et donc d'obtenir une quantité d'information maximum dans un temps donné Dt. Comme on le verra par la suite, L'utilisation de tels detecteurs de haute résolution, en association avec l'utilisation des probabiLités de fausse détection (ou des quantités d'information du genre de celles de La formule (6) ou de la formule (7)), permet d'abaisser considérablement La probabilité de défaut de détection de l'explosif, jusqu'à des niveaux bien inférieurs à 0 Certes, on pourrait utiliser des détecteurs de faible résolution tout en augmentant le flux de neutrons par unité de temps. Cependant, ceci conduirait à une augmentation de l'importance des problèmes de radio-protection inhérents aux systèmes de détection d'explosif par irradiation neutronique et conduirait surtout à des empilements dans les circuits électroniques analogiques de traitement des signaux fournis par Les détecteurs de rayonnement gamma. Pour diminuer PF et PD et pour augmenter la quantité d'information, on augmente donc la résolution des moyens de détection, ce qui entraine une réduction du bruit dans La même proportion. Si lton suppose que La densité de probabilité sous un pic i, en fonction de L'énergie des photons gamma, est gaussienne, l'écart-type s du comptage des événements est proportionnel au bruit de fond physique M2i sous te pic.Donc, en supprimant L'indice i pour pLus de simplicité et en considérant deux détecteurs repérés seuLement par des indices A et B, on peut écrire sA/sB = M2A/M2B (9) En conséquence, les efficacités des détecteurs étant supposées identiques et Les comptages dûs au signaL étant aussi, de ce fait, identiques dans les deux détecteurs, si les quantités d'information IA et IB, correspondant respectivement aux détecteurs A et B, sont exprimées par la formule (8), on obtient que chacune de ces quantités d'information IA et IB, obtenues pendant le même temps Dt, est inversement proportionnelle à la racine carrée de l'écart-type correpondant. En revanche, si l'on utilise la formule (7) pour calculer tes quantités d'information lA et IB, si le taux de signal MIA est égal au taux de signal M1B et si Le bruit M2A est inférieur au bruit M2B, on obtient que le rapport IA/IB est inférieur au rapport sB /sA . Ceci résulte du fait que la quantité dtinformation exprimée par la formule (7) décroît de façon monotone en fonction du bruit tout en restant supérieure à la quantité d'information exprimée par la formule (8), comme on le voit sur La figure 3 où l'on a choisi une valeur particulière du comptage dans la voie signal, à savoir M1=100 événements durant le temps Dt. A titre d'exemple, on considère Le cas de la raie à 10,83 MeV de l'azote, on prend un écart-type de 500 keV environ pour le détecteur d'indice B, ce qui caractérise typiquement La faible résolution d'un détecteur en germanate de bismuth ; on considère que La largeur de la bande d'énergie autour de 10,83 MeV, dans laquelle sont détectés les photons gamma, est égale à 1,3 MeV ; pour Le détecteur d'indice A, on prend un écart-type de 10 keV et une bande d'énergie de largeur 20 keV autour de 10,83 MeV, ce qui correspond typiquement à un détecteur en germanium intrinsèque ; on obtient alors un rapport IA/IB supérieur à 8 ; plus précisément, si L'on applique La formule (2), à des -2 valeurs de 10 environ pour PF et PD (M2=120 et M1=130 événements pendant le temps Dt), valeurs qui sont accessibles à l'aide d'un détecteur de faible résolution, il correspond, pour un détecteur de haute -26 résolution, une valeur proche de 10 et par conséquent une sécurité de détection quasiment absolue. On notera également que des moyens de détection de haute resolution permettent d'augmenter le débit des bagages par unité de temps, en réduisant la durée d'emmagasinage, en vue d'obtenir une valeur -2 entre 10-2 -26 intermédiaire entre 10 et 10 26 pour PF et pour PD. La figure 4 est une vue en perspective schématique et partielle d'un système conforme à l'invention, destiné à contrôler des bagages pour s'assurer qu il ne contiennent pas d'explosif. Le système représenté sur la figure 4 comprend une enceinte 2 de thermalisation de neutrons et, dans cette enceinte, au milieu de cette dernière par exemple, une source 4 de neutrons rapides, par exemple une source de neutrons de 14 MeV. L'enceinte 2 est faite d'une matière apte à thermaliser Les neutrons rapides, par exemple une matière hydrogénée telle que le polystyrène. Pour gagner du temps, le système est prévu pour permettre le contrôle de deux bagages 6 en même temps et comprend à cet effet deux tapis roulants 8 qui traversent l'enceinte 2 et qui permettent de faire circuler Les bagages 6 de part et d'autre de La source 4, l'enceinte 2 comportant, pour chaque tapis roulant 8, une entrée 10 permettant L'entrée des bagages dans l'enceinte et une sortie 12 permettant aux bagages de quitter l'enceinte. Les tapis roulants 8 sont faits d'un matériau mince qui présente une faible absorption vis-à-vis du rayonnement gamma et contient le moins possible d'éléments parmi ceux que l'on veut détecter. Le système représenté sur la figure 4 comprend également des moyens de détection du rayonnement gamma émis par Les bagages lorsqu'ils sont irradiés par des neutrons. Ces moyens de détection comprennent, pour chaque tapis roulant 8, une matrice de détecteurs 14, de préférence à haute résolution. Chaque matrice de détecteurs 14 est placée sous te tapis roulant 8 correspondant, au niveau de la source 4, de telle façon que chaque bagage 6 à controler, qui circule sur ce tapis roulant et y repose de préférence par sa plus grande surface, puisse être mis en présence de La source 4 et de la matrice de détecteurs correspondant à ce tapis roulant. Dans une variante non représentée, les matrices de détecteurs sont placées, non pas sous leurs tapis roulants respectifs, mais de part et d'autre de La source, entre les deux tapis. On voit sur La figure 4 que chaque détecteur 14 est protégé du rayonnement direct de la source de neutrons et cotlimaté vers le bagage 6 à contrôler. A cet effet, on utilise un bloc ou blindage 16 qui est fait d'un matériau apte à arrêter les neutrons et dans lequel sont prévus des Logements qui accueillent chacun un détecteur 14. La figure 5 illustre schématiquement des moyens éLectroniques de traitement des impulsions fournies par Les détecteurs 14 (qui fonctionnent simuLtanément). Ces moyens électroniques de traitement comprennent une pluralité de chaînes de détection. Chaque détecteur 14 est associé à une chaîne de détection qui comprend successivement, à la suite du détecteur 14 associé, un préamplificateur 18, un amplificateur proportionnel 20, un discriminateur à portes 22 et un codeur d'amplitude 24. La Linéarité intégrale et la Linéarité différentielLe du codeur d'ampLitude doivent être suffisantes pour exploiter pleinement les possibilités que permettent L'emploi d'un détecteur de haute résolution. Le nombre et La disposition des détecteurs dépend de La forme et des dimensions des bagages à contrôler. Le nombre de détecteurs et donc de chaînes de détection dépend aussi de l'angle solide observable par un détecteur, de la valeur maximum de PF choisie pour caractériser La présence d'un explosif et du nombre de neutrons émis par la source par unité de temps. Comme on le voit sur La figure 5, les codeurs d'amplitude sont reliés à un calculateur 26 qui est lui-même reLié à des moyens de signalisation 28 (aptes à afficher les résultats d'un contrôle et/ou à émettre un signaL sonore s 'il y a présomption d'explosif dans un bagage contrôlé), le calculateur recevant les valeurs des codages d'amplitude relatives à chaque channe. De préférence, la source 4 est une source puLsée, apte à fournir des bouffées de neutrons ("neutron bursts" dans les publications en langue angLaise) et, comme on Le voit sur les figures 4 et 5, cette source puLsée est reliée à chacun des discriminateurs à porte 22 par une ligne de synchronisation 30, en vue d'effectuer la détection des photons gamma prompts en corncidence avec les bouffées de neutrons rapides et la détection des photons gamma de capture entre ces bouffées de neutrons. Le système représenté sur les figures 4 et 5 fonctionne de la façon suivante : un bagage, placé par un tapis roulant 6, pénètre dans L'enceinte de thermalisation 2 et est mis en présence de La source 4 et de La matrice de détecteurs 14 correspondant au tapis roulant. La source émet des neutrons rapides qui sont thermalisés par tes parois de l'enceinte de sorte que le bagage est aussi irradié par des neutrons thermiques et émet essentiel Lestent des photons gamma prompts, de capture et d'activation. Les détecteurs 14 de la matrice reçoivent ces photons gamma et fournissent des signaux aux moyens éLectroniques de traitement qui Les traitent de la façon suivante On suppose qu'on a choisi d'examiner des raies de L'azote et des raies de l'oxygène par exemple. Le caLculateur a en mémoire La quantité M'2i correspondant à la raie i, pour chacune des n raies choisies. Chaque chaîne de détection permet au calculateur 26 d'acquérir une valeur de comptage Ni pour chaque raie i. Le calculateur peut ainsi caLcuLer, grâce à La formule (2), la probabilité PFi pour chaque raie i et pour chaque chaîne de détection. Le calculateur 26 calcule alors, pour chaque chaîne de détection, la quantité PF égale au produit des quantités PFi conformément à la formule suivante : PF étant La probabilité pour que l'ensemble des comptages effectués sous tous les pics estimés pertinents pour la détection d'explosif soit dû uniquement au bruit ou à une quantité d'explosif inférieure à la quantité maximum présente dans des bagages ne contenant aucun explosif. Ensuite, pour chaque chaine de détection, le calculateur 26 compare PF à un seuil S1 fixé par Les utilisateurs. Si, pour l'une au moins des chaînes de detection, PF est inférieur à S1, Le bagage est présumé contenir un explosif et le calculateur commande alors Les moyens de signalisation 28. Ce bagage suspect est alors dirigé jusqu'à une zone de stockage des bagages suspects pour y être examiné. Si, au contraire, PF est au moins égal à S1 pour chacune des chaines de détection, le bagage est considéré comme ne contenant aucun explosif et quitte L'enceinte de thermalisation 2 par La sortie correspondant au tapis roulant sur lequel le bagage se trouve. Dans une variante de réalisation du système représenté sur Les figures 4 et 5, le calculateur calcule, pour chaque chaîne de détection, La probabiLité PFi pour chaque raie i puis, par la formule (6), la quantité d'information li correspondante. Ensuite, Le calculateur calcule la somme de ces quantites d'information (pour chaque channe de détection). Chaque somme est ensuite comparée à un seuil S2 fixé par Les utiLisateurs. Si la somme est supérieure à S2 pour l'une quelconque des chaines de détection, le bagage est présumé contenir un explosif, les utilisateurs sont avertis de ceci et Le bagage est dirigé vers La zone des bagages suspects. Si, au contraire, la somme est inférieure ou égale à Sv pour chacune des channes de détection, Le bagage est considéré comme ne contenant aucun explosif et quitte L'enceinte de thermalisation. Dans une autre variante de réalisation, au lieu dtutiliser la formule (6) pour calculer la quantité 'information Ii, on utilise la formule (7), Le calculateur ayant en mémoire les fonctions K et a, ce qui conduit à un temps de calcul plus court. La figure 6 ilLustre schématiquement d'autres moyens électroniques de traitement des impulsions fournies par les détecteurs 14. Ces autres moyens éLectroniques de traitement diffèrent de ceux qui sont représentes sur la figure 5 par le fait que la calculateur 26 est remplacé par une pluralité de processeurs 32 qui sont respectivement associés aux chaînes de détection (comportant chacune les éLéments 18, 20, 22 et 24), chaque processeur 32 étant placé à ta suite du codeur 24 correspondant et relié à ce dernler. Dans le cas des moyens électroniques de traitement représentés sur La figure 6, chaque processeur 32 calcule le produit PF (ou la somme des quantités d'informations) pour La chaine de détection à laquelle il est associe. Tous tes processeurs 32 sont reliés à des moyens de signaLisation 34 qui permettent d'avertir les utilisateurs du système lorsqu'il s'avère, grâce à l'une au moins des chaînes de détection, qu'il y a présomption de présence d'explosif. Sur la figure 7, on a représenté schématiqueient et partiellement un autre système de détection conforme à l'invention, qui diffère de celui qui est représenté sur La figure 4 par Le fait qu'il utiLise un seul détecteur de haute résolution 36 au lieu d'un matrice de tels détecteurs, ce détecteur 36 étant placé comme précédemment et logé dans des moyens 38 de cotlimation et de protection vis-à-vis du rayonnement direct de la source 4. Des moyens non représentés permettent le dépLacement du détecteur 36, associé aux moyens 38, suivant un trajet 40 lui permettant de contrôler la totalité du bagage 6, ce dernier étant immobilisé. Les signaux fournis par le détecteur 36 sont envoyés à une chaîne de détection composée, comme précédemment, d'un préamplificateur 18, d'un ampLificateur proportionnel 20, d'un discriminateur à porte 22 et d'un codeur d'amplitude 24, ce dernier étant relié à un calculateur 42, Lui-même relié à des moyens de signalisation 44. Les moyens de déplacement du détecteur 36 permettent d'immobiliser celui-ci sous le bagage 6 en plusieurs points et, si aucun explosif n'a été détecté en chacun de ces points, te bagage est considéré comme ne contenant pas d'explosif alors que si pour l'un des points, il y a présomption de présence d'explosif, le bagage est dirigé vers La zone des bagages suspects pour y être examiné. Sur La figure 8, on a représenté schématiquement et partiellement un autre système de détection conforme à l'invention. Le système représenté sur la figure 8 diffère de cetui qui est représenté sur La figure 4 par le fait qu'il a une configuration "ouverte" et donc qu'il ne comporte pas d'enceinte de thermalisation des neutrons. Le système représenté sur la figure 8 comprend encore une source puLsée 46 de neutrons rapides, mise en présence du bagage à contrôler. Le système représenté sur la figure 8 comprend également un détecteur 50 de haute résolution pour détecter le rayonnement gamma et ce détecteur est, comice précédemment, relié à une chaîne de détection du genre de celle qui est représentée sur la figure 7, cette chaîne de détection étant reliée à un calcuLateur, Lui-m ême relié à des moyens de signalisation, Comme précédemment, la source 46 est reliée par une ligne de synchronisation au discriminateur à porte de La chaîne pour les raisons expliquées pLus haut. De plus, Le détecteur 50 est protégé du rayonnement direct de la source 46 et collimaté, par des moyens appropriés 52, vers le bagage 6 à contrôler Ainsi, chaque bagage à contrôleur se trouve en regard du détecteur 50 collimaté, la source 46 étant placée de façon à pouvoir irradier Le bagage 6 par Les neutrons rapides qu'elle engendre, de sorte que te bagage 6 émet essentiellement un rayonnement gamma qui est dû aux neutrons rapides et qui arrive au détecteur 50. tes impulsions émises par ce dernier sont traitées par les moyens électroniques de traitement associés à ce détecteur 50. Bien entendu, au Lieu d'un seul détecteur, on pourrait utiliser une matrice de détecteurs, chaque détecteur de la matrice étant associé à une chaîne de détection comme on L'a expliqué en se référant aux figures 5 et 6. La présente invention n'est pas limitée à la détection d'expLosifs nitrés. Elle s' applique également, par exemple, à La détection d'explosifs non nitrés et même à la détection d'un grand nombre de composés chimiques, en expLoitant les signatures gamma reLatives à des éLéments chimiques de ces composés. L'invention s'applique en particulier à la détection de produits stupéfiants, en exploitant notamment les pics du carbone et de L'azote. En prenant l'exemple d'un contrôLe de bagages pour s'assurer qu'ils ne contiennent pas d'explosif, on explique ci-aprés L'obtention préalable des paramètres et fonctions (M2i, Ci 14kan, K et a) qu'iL convient de mémoriser pour faire les divers calculs. Le bruit de fond physique M2i sous le pic i est estimé par extrapolation de La densité de probabilité de part et d'autre de ce pic. La contribution maximum (Ci) au comptage des matériaux non explosifs du bagage typique est mesurabLe si l'on dispose d'un nombre suffisant de bagages de voyageurs, remplis de produits habituels (vêtements, objets de toilette ...) : on fait passer de tels bagages (sans explosif) dans le système de la figure 4 (ou on les place successivement devant le système de la figure 8), d'où divers comptages dont on fait la moyenne pour obtenir Ci. Il est à noter que le calcul de PFi peut être fait directement, à partir de la formule (2), ou bien en tabulant d'une façon ou d'une autre les valeurs de PFi en fonction de Ni et de M2i, par exemple comme cela est expliqué dans le document EP-A-0200596. Les valeurs de K(M'2i) et de a(M'2i) peuvent être tabulées dans les intervalles utiles de M'2i. Pour faire une telle tabulation préalable, on détermine des courbes du genre de celles de la figure 1 : en fixant les paramètres M1 et M2, on calcule PF (formule (2)) puis I (formule (6)) et ce, pour diverses valeurs du coup le (M1, M2). A partir des courbes ainsi obtenues et pour diverses valeurs du paramètre M2, on détermine K et a à une valeur fixée de M2 correspond une courbe du réseau de courbes obtenu ; on choisit deux points PI et P2 de cette courbe d'où deux couples (I(P1), M1(P1)) et (I(P2 ), M1(P2)) respectivement associés à ces deux points ; à partir de la formule (7) on obtient alors un système de deux équations à deux inconnues dont la solution donne K(M2) et a(M2) pour la valeur fixée de M2. REVENDICATIONS 1. Système de détection d'une substance susceptible d'être contenue dans un objet (6), système caractérisé en ce qu'il comprend - des moyens (4, 46) d'irradiation de l'objet par des neutrons, - des moyens (14, 36, 50) de détection du rayonnement gamma susceptible d'être alors émis par l'objet, et - des moyens électroniques (18, 2û, 22, 24, 26, 28, 32, 34, 42, 44) de traitement des signaux fournis par Les moyens de détection, ces moyens électroniques de traitement étant prévus pour compter les photons gamma correspondant à chaque raie d'un ensemble de raies caractéristiques d'au moins un élément chimique de la substance, déterminer, pour chaque raie i, une probabilité de fausse détection de L'élément chimique associé à cette raie, c'est-à-dire la probabilité PFi pour que le signal détecté, correspondant à cette raie i, soit dû à un bruit de fond, déterminer le produit de ces probabilités de fausse détection, comparer ce produit à un seuil fixé par les utilisateurs du système, et . avertir ces utilisateurs si le produit est inférieur au seuil fixé par eux, L'objet étant alors présumé contenir la substance. 2. Système selon la revendication 1, caractérisé en ce que chaque probabilité PFi est déterminée par la formule : où Ni représente Le comptage, pendant un temps Dt, qui correspond à La raie i et qui est dû à des réactions nucléaires induites par des neutrons et aussi à un bruit de fond M'2i reLatif à la raie i, pendant le temps Dt. 3. Système selon la revendication 2, caractérisé en ce que chaque bruit de fond M'2i est déterminé en faisant la somme d'un bruit de fond physique M2i, reLatif à La raie i, et de la borne supérieure Ci de comptages relatifs à cette raie i, pendant le temps Dt, sur des objets qui sont susceptibles de contenir ta substance mais dont on est sûr qu'ils n'en contiennent pas. 4. Système selon l'une quelconque des revendications 1 à 3, caractérisé en ce que les moyens électroniques de traitement (18, 20, 22, 24, 26, 28, 32, 34, 42, 44) sont prévus pour : déterminer, pour chaque raie i, une quantité d'information Ii relative à cette raie i et définie par la formule PFi Ii = - Log ------- 1 - PF1 déterminer La somme de ces quantités d'information, comparer cette somme à un seuil fixé par les utilisateurs du système, et avertir ces utilisateurs si cette somme est supérieure au seuil auquel elle est comparee, l'objet étant alors présumé contenir La substance. 5. Système selon l'une quelconque des revendications 1 à 3, caractérisé en ce que les moyens électroniques de traitement (18, 20, 22, 24, 26, 28, 32, 34, 42, 44) sont prévus pour déterminer, pour chaque raie i, une quantité dtinformation Ii relative à cette raie i et définie par la foreuse :: où K et a sont des fonctions mémorisées d'un bruit de fond M'2i relatif à la raie i et Nîi est la partie entière de La différence Ni-M'2i, Ni représentant le comptage, pendant un temps Dt, qui correspond à la raie i et qui est dû à des réactions nucléaires induites par des neutrons et aussi au bruit de fond M'2i relatif à la raie i, pendant le temps Dt, déterminer la somme de ces quantités d'information, . comparer cette somme à un seuil fixé par les utilisateurs du système, et avertir ces utilisateurs si cette somme est supérieure au seuil auquel elle est comparée, l'objet étant alors présumé contenir la substance. 6. Système selon l'une quelconque des revendications 1 à 5, caractérisé en ce que les moyens d'irradiation comprennent : - une source C4) de neutrons rapides, et - une enceinte (2) de thermalisation de ces neutrons rapides, dans Laquelle se trouve La source et qui est prévue pour recevoir l'objet (6). 7. Système selon t'une quelconque des revendications 1 à 5, caractérisé en ce que les moyens d'irradiation comprennent une source (46) de neutrons rapides qui est prévue pour irradier l'objet (6). 8. Système selon L'une quelconque des revendications 1 à 7, caractérisé en ce que les moyens de détection comprennent une pluralité de détecteurs (14) de rayonnement gamma et en ce que les moyens électroniques de traitement comprennent une pluralité de chaînes de détection (18, 20, 22, 24), respectivement associées aux détecteurs. 9. Système selon t'une quelconque des revendications 1 à 8, caractérisé en ce que les moyens de détection (14, 36, 50) du rayonnement gamma sont des doyens de détection de haute résoLution. 10. Système selon l'une quelconque des revendications î à 9, caractérisé en ce que, les moyens d'irradiation comprenant une source (4, 46) de neutrons, les moyens de détection (14, 36, 50) sont protégés du rayonnement direct de la source et collimatés vers l'objet. 11. Système selon l'une quelconque des revendications 1 à 10, caractérisé en ce que les moyens d'irradiation comprennent une source pulsée (4, 46) de neutrons, prévue pour fournir des bouffées de neutrons et en ce que Les moyens électroniques de traitement coopérent avec les moyens de détection pour faire des mesures dans des intervalles de temps où l'on est sûr de detecter seuLement l'une des catégories de photons gamma engendres lors de l'irradiation de l'objet par des neutrons. 12. Système selon l'une quelconque des revendications 1 à 11, caractérisé en ce que les moyens d'irradiation comprennent une source (4, 46) de neutrons de 14 MeV engendrés par des réactions de fusion.