FR 2939895 A1 20100618 FR 0858564 A 20081215 PROCÉDÉ DE DÉTECTION NON INTRUSIVE D'ELEMENT CHIMIQUE Domaine technique et art antérieur L'invention concerne un procédé de détection non intrusive d'au moins un élément chimique présent dans un liquide ou un solide. L'invention s'applique, par exemple, à la détection d'explosif(s), d'arme(s) chimique(s) ou de drogue(s) dans le domaine de la sécurité civile (inspection de bagages, conteneurs de fret, etc.), à la recherche de produits de contrebande (contrôles douaniers), au déminage (sécurité civile et applications militaires) ou encore à la caractérisation des matériaux constitutifs des déchets radioactifs issus de l'industrie nucléaire. Différentes méthodes sont connues pour la détection non intrusive d'éléments chimiques. Certaines de ces méthodes sont basées sur la mise en oeuvre de mesures neutroniques. Parmi les méthodes basées sur la mise en oeuvre de mesures neutroniques, certaines identifient les éléments chimiques par détection de rayonnements gamma de capture (n, y) induits par des neutrons thermiques, ou par détection de rayonnements gamma de diffusion inélastique (n, n'y) ou encore par détection de rayonnements gamma faisant suite à d'autres réactions induites par des neutrons rapides, à savoir (n, p) , (n, d) , (n, t) , (n,î), etc.. Dans la suite de la description, l'expression � rayonnement gamma de type inélastique � est utilisée pour tout rayonnement gamma induit par des neutrons rapides. La détection de rayonnements de capture induits par des neutrons thermiques est plus sensible en utilisant un générateur neutronique pulsé pour l'émission des neutrons car le rapport signal (rayonnement de type capture) sur bruit (rayonnement de type inélastique) est maximum entre les impulsions. Une illustration de ce type de détection est fournie dans le brevet US 6 393 085 B1. Dans ce document on décrit aussi la détection possible de rayonnements gamma de type inélastique, cette fois-ci pendant les impulsions, mais une telle méthode présente une sensibilité de détection nettement moindre qu'en utilisant un tube à particule associée (TPA), qui permet de localiser spatialement l'objet inspecté à partir d'une mesure de temps de vol neutronique dont le principe est décrit en détail dans la demande de brevet US 2007/0241283 Al. En bref, chaque neutron est émis simultanément à une particule a et dans la direction opposée. La détection de la particule î permet de déterminer la direction d'émission du neutron et la mesure du temps séparant la détection de la particule î et du rayonnement y de type inélastique induit par le neutron permet de déterminer le temps de vol neutronique. On dispose donc d'un moyen de localiser en 3D le lieu d'interaction et donc d'éliminer les événements produits en dehors de la zone d'intérêt, d'où une nette amélioration du rapport signal sur bruit par rapport à l'interrogation pulsée. Cependant, comme l'émission du TPA est continue, il s'ensuit que, avec cette méthode, la détection des rayonnements de capture n'est pas optimale en terme de rapport signal sur bruit. Un problème réside en ce que ni l'une ni l'autre des deux méthodes, l'une basée sur le tube pulsé, l'autre sur le TPA, ne permet de détecter simultanément avec la meilleure sensibilité possible les rayonnements de capture et ceux de type inélastique, et par suite d'identifier de manière satisfaisante tous les éléments recherchés. Une fois une première méthode effectuée, il faut ensuite effectuer la deuxième méthode pour élargir le champ des éléments détectables, ou améliorer le niveau de confiance sur la présence de ceux détectables par les deux méthodes. La durée totale d'une opération de détection est en conséquence relativement longue, ce qui représente un inconvénient majeur en terme de sécurité. L'invention ne présente pas cet inconvénient. Exposé de l'invention En effet, l'invention concerne un procédé de détection de présence d'un élément chimique dans un objet par émission de neutrons sur l'objet, caractérisé en ce que l'émission de neutrons sur l'objet est constituée, d'une part, par une émission continue de neutrons provenant d'un générateur neutronique à particule associée et, d'autre part, par une émission d'impulsions neutroniques qui se superposent à l'émission continue de neutrons, les impulsions neutroniques provenant d'un générateur neutronique pulsé qui génère des impulsions neutroniques de durée d'impulsion T2, deux impulsions neutroniques successives étant séparées d'une durée T4, les émissions continue et pulsée de neutrons sur l'objet produisant un rayonnement gamma de capture et un rayonnement gamma de type inélastique. Selon une caractéristique supplémentaire de l'invention, le procédé comprend : - une détection, à l'aide de la technique de la particule associée, pendant les durées T4 qui séparent les impulsions de durées T2, du rayonnement gamma de type inélastique et du rayonnement gamma de capture, û une formation , dans une plage de temps de vol neutronique relative à l'objet, d'un premier spectre gamma à partir du rayonnement gamma de type inélastique et du rayonnement gamma de capture détectés à l'étape de détection à l'aide de la technique de la particule associée, - une formation, dans une plage temporelle disjointe de la plage de temps de vol neutronique relative à l'objet ou de toute autre plage de temps de vol neutronique relative à tout autre objet susceptible d'environner l'objet, d'un spectre de bruit de fond aléatoire relatif aux coïncidences fortuites, à partir du rayonnement gamma de type inélastique et du rayonnement gamma de capture détectés à l'étape de détection à l'aide de la technique de la particule associée, - une formation, à partir du spectre de bruit de fond aléatoire relatif aux coïncidences fortuites, d'un spectre de bruit de fond aléatoire normalisé sur la plage temporelle relative à l'objet, - une soustraction du premier spectre gamma et du spectre de bruit de fond aléatoire normalisé pour former un deuxième spectre gamma, et - une étape de recherche de la présence d'un élément chimique dans l'objet à partir du deuxième spectre gamma et de spectres d'éléments chimiques de référence susceptibles d'être détectés pour délivrer une première donnée apte à traduire la présence ou l'absence de l'élément chimique dans l'objet. Selon une autre caractéristique supplémentaire de l'invention, le procédé comprend, en outre : - une détection, à l'aide de la technique d'interrogation neutronique pulsée, du rayonnement gamma de capture et du rayonnement gamma de type inélastique, d'une part pendant les durées d'impulsions T2 et, d'autre part, pendant les durées T4 qui séparent les impulsions de durée T2. - une formation d'un troisième spectre gamma, à partir du rayonnement gamma de capture et du rayonnement gamma de type inélastique détectés, pendant les durées T4, à l'aide de la technique d'interrogation neutronique pulsée, - une formation, à partir du troisième spectre gamma, d'une première donnée représentative d'au moins une raie de capture de l'élément chimique, - une formation d'un quatrième spectre gamma constitué d'une pluralité de canaux d'énergie, à partir du rayonnement de type inélastique détecté pendant les durées d'impulsion T2, à l'aide de la technique d'interrogation neutronique pulsée, - une formation d'un spectre gamma normalisé à partir du quatrième spectre gamma, la formation du spectre gamma normalisé étant effectuée en multipliant le contenu de chaque canal d'énergie du quatrième spectre gamma par un facteur de normalisation F sensiblement proportionnel au nombre de neutrons émis pendant une durée T4 et sensiblement inversement proportionnel au nombre de neutrons émis pendant une durée T2, - une formation, à partir du spectre gamma normalisé, d'une deuxième donnée représentative d'un bruit de fond de rayonnement gamma de type inélastique, - une étape de soustraction de la première donnée et de la deuxième donnée, une étape de validation de la présence ou de l'absence de l'élément chimique dans l'objet en fonction du résultat de l'étape de soustraction, et si la présence de l'élément chimique est validée dans l'objet, une étape supplémentaire qui délivre une deuxième donnée apte à confirmer ou infirmer la validation de la présence de l'élément chimique sur la base du résultat de l'étape de validation et de données de poids relatives aux différentes raies gamma de rayonnement de capture de l'élément chimique, ces poids étant attribués en fonction d'une probabilité de détection et d'un risque d'interférence associés à chaque raie. Selon encore une autre caractéristique supplémentaire de l'invention, le facteur de normalisation F est donné par l'équation : F = (I4xT4)/(I2xT2)x[(1-TM(T4))/(1-TM(T2))], où TM(T4) et TM(T2) sont, respectivement, un temps mort d'une électronique de traitement qui traite les signaux électriques issus de la détection des rayonnements gamma pendant les durées T4 et un temps mort d'une électronique de traitement qui traite les signaux électriques issus de la détection des rayonnements gamma détectés pendant les durées T2. Selon encore une caractéristique supplémentaire du procédé de l'invention, une étape de décision délivre un signal d'alerte si la première donnée ou la deuxième donnée confirme la présence de l'élément chimique. Selon encore une caractéristique supplémentaire de l'invention, dès lors que la présence de plusieurs éléments chimiques est confirmée, il comprend une étape supplémentaire de calcul de proportions relatives de tout ou partie des éléments chimiques dont la présence est confirmée. Selon encore une caractéristique supplémentaire du procédé de l'invention, la détection, à l'aide de la technique de la particule associée, du rayonnement gamma de capture et du rayonnement gamma de type inélastique pendant les durées T4 qui séparent les impulsions de durée T2 est effectuée par un détecteur préférentiellement adapté pour la détection de rayonnement gamma de type inélastique. Selon encore une caractéristique supplémentaire du procédé de l'invention, la détection, à l'aide de la technique d'interrogation neutronique pulsée, du rayonnement gamma de capture et du rayonnement gamma de type inélastique, d'une part, pendant les durées d'impulsion T2 et, d'autre part, pendant les durées T4 qui séparent les impulsions de durée T2 est effectuée par un détecteur préférentiellement adapté pour la détection de rayonnement gamma de capture. Une utilisation simultanée d'un générateur neutronique à particule associé G1 qui émet des neutrons en continu et d'un générateur neutronique pulsé G2 qui émet des impulsions neutroniques présente - a priori - un inconvénient pour la détection des rayonnements de capture du fait de la présence d'une émission neutronique résiduelle entre les impulsions du générateur G2, laquelle émission résiduelle est due à l'émission continue du générateur G1. Cette émission résiduelle dégrade en effet le rapport signal (rayonnements de capture) sur bruit (rayonnements de type inélastique). De façon avantageuse, le procédé de l'invention propose des moyens pour supprimer cet inconvénient. Dans le cadre de l'invention, l'utilisation combinée des générateurs G1 et G2 permet alors d'obtenir, en une seule acquisition, des mesures exploitables en terme de sensibilité, pour les rayonnements gamma de capture détectés par interrogation neutronique pulsée et pour les rayonnements gamma de type inélastique détectés par la méthode de la particule associée. Le champ des éléments détectables s'en trouve élargi sans augmenter le temps d'acquisition. De plus, pour les éléments détectables par les deux méthodes, le niveau de confiance de la détection est amélioré. Avec le système de détection de l'invention, il est avantageusement possible d'utiliser une large fraction de la période T des impulsions neutroniques (typiquement 95%) pour détecter simultanément les deux types de rayonnements induits par les neutrons thermiques (capture) et par les neutrons rapides (diffusion inélastique et réactions (n,î), (n,p), etc.) Brève description des figures D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture 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 non intrusive qui met en oeuvre le procédé de l'invention ; - la figure 2 représente un spectre d'émission de neutrons conforme à l'invention ; - la figure 3 représente les étapes d'un cycle de traitement du procédé de détection non intrusif de l'invention ; - la figure 4 représente un exemple de spectre temporel pour la technique de la particule associée ; - la figure 5 représente des étapes élémentaires d'une étape particulière du cycle de traitement représenté en figure 3 ; - la figure 6 représente un pic de détection d'une raie gamma dans un spectre gamma ; Sur toutes les figures les mêmes références désignent les mêmes éléments. Description détaillée de modes de réalisation préférentiels de l'invention La figure 1 représente un schéma de principe de système de détection non intrusive qui met en oeuvre le procédé de l'invention. Autour de l'objet 0 susceptible de contenir des éléments chimiques à détecter, le système de détection de l'invention comprend un générateur neutronique à particule associée G1, un générateur neutronique pulsé G2, deux détecteurs Dl et D2, un circuit électronique de traitement E qui traite et met en forme les signaux issus des détecteurs Dl et D2, des bases de données Bo, B1, B2 qui stockent les données délivrées par le circuit électronique de traitement E, des bases de données de référence BRC et BS qui stockent, respectivement, des données relatives à des raies de capture de référence et à des spectres de référence de rayonnement gamma de type inélastique et un calculateur C qui calcule un signal qui traduit la présence ou l'absence d'élément(s) chimique(s) ou de produits dangereux à partir des données contenues dans les bases de données Bo, B1, B2, BRC et B. Les générateurs neutroniques G1 et G2 émettent des neutrons dont une partie se dirige vers l'objet à étudier. Les générateurs neutroniques G1 et G2 fonctionnent simultanément (la figure 2 représente l'évolution temporelle de l'intensité d'émission des neutrons émis simultanément par les générateurs neutroniques G1 et G2). Un rayonnement gamma y1 résulte de l'interaction des neutrons émis par le générateur G1 avec l'objet à étudier et un rayonnement gamma y2 résulte de l'interaction des neutrons émis par le générateur G2 avec l'objet à étudier. Chacun des rayonnements gamma y1 ou y2 est susceptible d'être un rayonnement gamma de capture ou un rayonnement gamma de type inélastique. Chacun des détecteurs Dl ou D2 est apte à détecter un rayonnement de capture et/ou un rayonnement de type inélastique. Toutefois, le détecteur Dl est préférentiellement adapté pour la détection de rayonnement gamma de type inélastique par la méthode de la particule associée (le détecteur Dl est ainsi choisi pour sa forte efficacité de détection et sa bonne résolution temporelle) et le détecteur D2 est préférentiellement adapté pour la détection de rayonnement gamma de capture par la méthode d'interrogation neutronique pulsée (le détecteur D2 est ainsi choisi pour sa bonne résolution en énergie). Le détecteur Dl est, par exemple, un cristal scintillateur de forte efficacité et le détecteur D2, par exemple, un cristal semi-conducteur à haute résolution en énergie. La figure 1 représente le cas où un noyau du composé chimique X1 subit une interaction avec un neutron du générateur G1, provoquant l'émission d'un rayonnement y1 qui est détecté par le détecteur Dl et où un noyau du deuxième composé chimique X2 subit une interaction avec un neutron du générateur G2, provoquant l'émission d'un rayonnement y2 qui est détecté par le détecteur D2. Un neutron issu du générateur G1 induit l'émission d'un rayonnement y1 dans l'objet 0, ici de type inélastique, lequel rayonnement y1 est détecté dans le détecteur Dl, et un neutron émis par le générateur G2 se ralentit du fait de multiples collisions dans l'objet 0 puis induit l'émission d'un rayonnement gamma y2 de capture, lequel rayonnement y2 est détecté dans le détecteur D2. L'invention concerne cependant tous les autres cas possibles de détection comme, par exemple, les cas où les noyaux d'un même élément chimique sont détectés à l'aide des deux rayonnements gamma y1 et y2, ou seulement par l'un des deux rayonnements y1 ou y2. De façon connue en soi, le générateur neutronique à particule associée G1 émet des neutrons par action d'un faisceau continu F1 de deutérium sur une cible de tritium T1. L'émission de neutrons s'effectue par une réaction dans laquelle une particule a est produite en même temps qu'un neutron. La particule a est détectée par un détecteur à localisation L qui appartient au générateur G1. Le détecteur à localisation L est, par exemple, constitué d'un ensemble de N détecteurs élémentaires dk (k=1, 2, N) assemblés sous forme de matrice. La connaissance du détecteur élémentaire dk qui, parmi les N détecteurs, détecte la particule a permet de situer avec précision la direction dans laquelle la particule a a été émise. Dans un autre mode de réalisation de l'invention, le détecteur à localisation L n'est pas constitué d'un ensemble de détecteurs élémentaires mais d'un seul bloc détecteur sur lequel sont positionnés plusieurs capteurs, la localisation de la particule a étant alors effectuée par un calcul de barycentre. Particule a et neutron étant émis dans des directions sensiblement opposées, la connaissance de la direction dans laquelle la particule a est émise permet de connaître la direction dans laquelle le neutron est émis. Dès qu'une particule a est détectée par le détecteur L, un signal S1 qui traduit la détection de la particule a est transmis au circuit électronique E. Le signal S1 est un signal de synchronisation qui autorise la prise en compte, par le circuit électronique E, d'une impulsion délivrée par le détecteur Dl. Le signal S1 comprend également une information qui identifie la direction dans laquelle le neutron a été émis (détecteur élémentaire dk qui a détecté la particule a dans le cas d'un détecteur matriciel ou position de la détection dans le cas d'un bloc détecteur). Selon le mode de réalisation décrit ci-dessus, le signal S1 de synchronisation est délivré par le détecteur L. Dans un autre mode de réalisation de l'invention, le signal de synchronisation S1 est délivré par le détecteur D1, dès lors que celui-ci a reçu une impulsion de rayonnement gamma. Le générateur neutronique pulsé G2 émet des neutrons sous forme d'impulsions neutroniques intenses par action d'un faisceau pulsé F2 de deutérium sur une cible de tritium T2. Le générateur neutronique pulsé G2 délivre également un signal impulsionnel électrique S2 qui reproduit la forme temporelle des impulsions neutroniques. Le signal S2 est transmis au circuit électronique de traitement E. Le signal S2 commande le circuit électronique de traitement E afin de prendre en compte et d'identifier les impulsions délivrées par le détecteur D2, d'une part, durant chaque impulsion neutronique intense et, d'autre part, entre deux impulsions neutroniques intenses successives. Le signal S2 commande également le circuit électronique de traitement E afin de prendre en compte les impulsions délivrées par le détecteur Dl uniquement entre les impulsions neutroniques intenses successives. En effet, l'acquisition de signal par la technique de la particule associée ne doit pas se faire pendant l'émission des impulsions neutroniques intenses, car il en résulterait alors une forte dégradation du rapport signal sur bruit. La base de données Bo stocke les événements EVO de type inélastique détectés par le détecteur D2, pendant les impulsions neutroniques intenses, sous la forme d'une liste des énergies EO déposées dans le détecteur D2, lesquelles énergies sont proportionnelles à l'amplitude des impulsions électriques générées par le détecteur D2 sous l'action des rayonnements gamma. Dans la suite de la description, ces énergies seront nommées "énergies des impulsions détectées" dans un but de simplification. La base de données B1 stocke les événements EV1 de type inélastique et de capture détectés par le détecteur D2, entre les impulsions neutroniques intenses, sous la forme d'une liste des énergies El des impulsions qui sont générées par le détecteur D2 sous l'action des rayonnements gamma. La base de données B2 stocke les événements EV2 de type inélastique et de capture résiduels qui sont détectés, à l'aide de la technique de la particule associée, par le détecteur Dl, entre les impulsions neutroniques intenses. Comme cela sera précisé ultérieurement, chaque événement EV2 rassemble un ensemble d'informations, à savoir : une énergie gamma, un temps de coïncidence entre la détection de la particule î et celle du rayonnement y associé et une information de position de détection de la particule a. Le calculateur C calcule le signal d'alerte S(alerte) sur la base des données stockées dans les bases de données Bo, B1r B2 et de données de référence stockées dans les bases de données BRC et BS comme cela sera décrit ultérieurement, en référence à la figure 3. La figure 2 représente l'intensité I des neutrons qui sont émis simultanément par les générateurs neutroniques G1 et G2. Des neutrons sont émis en continu par le tube G1 avec une intensité I4 et des impulsions neutroniques intenses sont émises avec une période T par le générateur neutronique G2. Pendant une impulsion neutronique intense, l'intensité globale des neutrons émis par les tubes G1 et G2 évolue de I4 à I2 pendant une durée T1 (par exemple égale à lus), conserve sensiblement la valeur I2 pendant une durée T2 (par exemple quelques dizaines de microsecondes) et diminue de I2 à I4 pendant une durée T3 sensiblement égale à T1. Ensuite, une durée T4 sépare la fin d'une impulsion neutronique intense du début de l'impulsion neutronique intense suivante. Du fait du caractère simultané de l'émission de neutrons par les tubes G1 et G2, il faut souligner que, contrairement à l'art antérieur connu dans le domaine de l'interrogation neutronique pulsée, le nombre de neutrons émis ne revient pas à zéro entre les impulsions neutroniques intenses mais conserve une valeur non nulle. Pendant la durée d'émission T2 d'une impulsion neutronique intense et pendant la durée T4 qui sépare deux impulsions neutroniques intenses, le signal de synchronisation S2 place le circuit électronique E en attente de réception de rayonnement gamma provenant du détecteur D2. Si un rayonnement gamma est détecté pendant la phase de durée T2 d'émission neutronique intense, il est dû, de façon la plus probable, aux réactions des neutrons rapides qui induisent un rayonnement gamma de type inélastique et, avec une faible probabilité, aux neutrons thermiques. Si un tel rayonnement est détecté, son énergie est mesurée puis convertie, préférentiellement, sous la forme d'une donnée numérique, par le circuit électronique E. Les données numériques représentatives des énergies ainsi mesurées constituent les événements EVO (i.e. les énergies EO des impulsions) qui sont stockés dans la base de données Bo. Si un rayonnement gamma est détecté pendant la durée T4 qui sépare deux impulsions neutroniques intenses, il peut provenir d'une réaction de capture ou d'une réaction de type inélastique. L'énergie d'un tel rayonnement est mesurée puis convertie, préférentiellement, sous la forme d'une donnée numérique, par le circuit électronique E. Les données numériques représentatives des énergies ainsi détectées constituent les événements EV1 (i.e. les énergies El des impulsions) qui sont stockés dans la base de données B. Par ailleurs, dès lors qu'une particule a est détectée pendant la durée T4, le signal S1 place le circuit électronique E en attente d'un rayonnement gamma détecté par le détecteur D1. Si un rayonnement gamma est détecté en coïncidence avec la particule a, l'énergie du rayonnement gamma est mesurée puis convertie, préférentiellement sous la forme d'une donnée numérique, par le circuit électronique E. Le circuit électronique E calcule également la durée Tay qui sépare la détection de la particule a de la détection du rayonnement gamma. La base de données B2 stocke les événements EV2 qui sont détectés par le détecteur D1 sous la forme de données numériques. Un événement EV2 rassemble l'énergie E2 du rayonnement gamma détecté, la durée Tay et une information ik qui identifie la position à laquelle la particule a a été le détecteur à localisation L. Une ces événements correspond à des gamma de capture. Cette fraction est préférentiellement supprimée du correspondant à l'objet observé, comme cela sera décrit ultérieurement en référence à la figure 3. L'ensemble des mesures effectuées sur une 25 période T (T=T1+T2+T3+T4) constitue un cycle de mesures. Après Nm cycles de mesures, par exemple 600000 cycles de mesure correspondant à une durée totale de mesure d'environ dix minutes, un cycle de traitement est effectué. 30 La figure 3 représente les différentes étapes mises en oeuvre dans un cycle de traitement de détectée sur fraction de rayonnements d'événements spectre gamma20 l'invention. Un cycle de traitement est mis en oeuvre par le calculateur C. Un cycle de traitement comprend tout d'abord un ensemble d'étapes de construction de spectres, à savoir : une étape 1 de construction d'un spectre gamma So à partir des données stockées dans la base de données Bo, une étape 2 de construction d'un spectre gamma S1 à partir des données stockées dans la base de données B1r et une étape 3 de construction d'un spectre gamma Sn dépourvu de bruit aléatoire à partir des données stockées dans la base de données B2, l'étape 3 étant elle-même constituée d'une étape 4 de sélection des événements stockés dans la base de données B2 pour déterminer, parmi ces événements, ceux qui correspondent à l'objet étudié, d'une étape 5 de construction d'un spectre gamma brut Sb de l'objet étudié à partir des événements sélectionnés à l'issue de l'étape 4, d'une succession d'étapes el, e2 pour former un spectre de bruit aléatoire normalisé Sa(n) à partir des événements stockés dans la base de données B2 et d'une étape 6 de soustraction du spectre gamma brut Sb et du spectre de bruit aléatoire normalisé Sa(n) pour former le spectre Sn dépourvu de bruit aléatoire qui ne comprend alors que les événements de type inélastique induits par les neutrons rapides dans l'objet étudié. Les étapes 1, 2 et 3 de construction de spectres sont suivies d'étapes de traitement de spectres, à savoir : - une étape 7 de recherche d'élément(s) chimique(s) à partir des spectres So et SI, et - une étape 8 de recherche(s) d'élément(s) chimique(s) à partir du spectre Sn. Les étapes 7 et 8 de traitement de spectres sont, quant à elles, suivies d'une étape de prise de décision 9. Par étape de construction de spectre à partir de données stockées, il faut entendre une étape, connue en soi, de construction des histogrammes des énergies des événements associés à ces données. L'étape 3 de construction du spectre Sn dépourvu de bruit aléatoire à partir des données stockées dans la base de données B2 va maintenant être décrite. Les événements de type inélastique et de capture EV2 qui sont stockés dans la base de données B2 ne sont constitués que d'événements qui se produisent entre les impulsions neutroniques. Le circuit électronique E comprend, à cette fin, des moyens qui sélectionnent, sur la base du signal S2 qui provient du générateur pulsé G2, les événements détectés par le détecteur Dl dans les plages temporelles souhaitées. Comme cela a été mentionné plus haut, l'étape 3 comprend une étape 4 de sélection des événements stockés dans la base de données B2. L'étape 4 de sélection des événements stockés dans la base de données B2 repose sur la technique connue de la particule associée. Cette étape est une étape connue en soi qui permet de déterminer, parmi tous les événements de la base B2, ceux qui concernent l'objet étudié. A partir des événements ainsi sélectionnés à l'étape 4, on construit le spectre gamma brut de l'objet étudié Sb lors de l'étape 5. L'étape 5 est une étape connue en soi, à l'image des étapes 1 et 2 mentionnées précédemment. A l'étape 5 succède une étape 6 qui soustrait du spectre Sb le spectre de bruit aléatoire normalisé Sa(n). Le spectre de bruit aléatoire normalisé Sa(n) est construit par les étapes successives el et e2, à partir de données spectrales stockées dans la base B2. Un spectre de bruit aléatoire Sa est tout d'abord formé à partir des données spectrales stockées dans la base B2 (étape el) et le spectre de bruit aléatoire normalisé Sa(n) est alors formé à partir du spectre de bruit aléatoire Sa (étape e2). La formation du spectre de bruit aléatoire normalisé Sa(n) est précisée ci-dessous en référence à la figure 4. La figure 4 représente un exemple de spectre temporel obtenu à l'aide de la technique de la particule associée. Ce spectre représente les temps séparant la détection de la particule alpha et du rayonnement gamma, temps à partir desquels on peut déduire les temps de vol neutroniques. Le spectre de la figure 4 comprend, à titre d'exemple, un niveau de bruit b représentatif du niveau des coïncidences fortuites, niveau au-dessus duquel se dessinent un pic p1 de très faible amplitude qui correspond aux photons induits dans le générateur G1, un pic p2 de forte amplitude qui représente les coïncidences vraies de l'objet étudié et un pic p3 d'amplitude moyenne qui représente les événements ne correspondant pas à l'objet étudié. La largeur temporelle At du pic p2 est la plage temporelle relative à l'objet étudié. L'étape el consiste, par exemple, à sélectionner les événements du spectre temporel dans une plage temporelle At(Sa) négative par rapport à la plage temporelle relative à l'objet étudié et disjointe de la plage temporelle relative à l'objet étudié, puis à construire le spectre Sa de bruit de fond aléatoire dû aux coïncidences fortuites à partir des événements ainsi sélectionnés. La durée At(Sa) est préférentiellement supérieure ou égale à la plage temporelle At de l'objet étudié afin d'obtenir une meilleure précision statistique. L'étape e2 consiste ensuite à normaliser le spectre Sa à la largeur At de la plage temporelle de l'objet étudié. Le facteur de normalisation est égal à At/At(Sa). Le spectre normalisé Sa(n) délivré à l'issue de l'étape e2 est alors retranché du spectre brut Sb à l'étape 6. Il en résulte le spectre gamma dépourvu de bruit aléatoire Sn de l'objet étudié. La présence d'élément(s) chimique(s) est alors recherchée, à partir du spectre Sn, à l'aide d'un algorithme de traitement qui est mis en oeuvre, de façon connue en soi, par l'étape 8 (cf. figure 3). La recherche d'élément(s) chimique(s) est effectuée à l'aide de données de référence répertoriées dans une liste préétablie contenue dans une base de données Bs. La base de données Bs contient des données de spectre complet pour chacun des éléments chimiques susceptibles d'être détectés. La résolution en énergie du détecteur Dl ne permet pas un traitement par raie contrairement au détecteur D2. L'étape 8 met donc en oeuvre une déconvolution de l'ensemble du spectre Sn pour extraire la contribution relative de chaque élément présent. La déconvolution est réalisée par un algorithme d'ajustement à l'aide des données contenues dans la base Bs. La fraction d'événements dans Sn associée à un élément X est ensuite normalisée au nombre total de coups dans Sn. On calcule alors le pourcentage FX de l'élément X dans le spectre complet avec la variance 6(FX) qui lui est associée. Un élément X est considéré comme détecté si, par exemple, l'inéquation suivante est vérifiée : 6(FX) / FX La figure 5 illustre les étapes élémentaires mises en oeuvre par l'algorithme Al. L'algorithme Al détermine la présence d'éléments chimiques qui correspondent à des raies gamma répertoriées dans la base de données de référence BRC. L'algorithme Al inclut la recherche d'interférences possibles des rayonnements de type inélastique via l'analyse du spectre Sc. De façon connue en soi, tout élément X recherché est caractérisé par un ensemble de n raies gamma de capture d'énergies Elx, E2x, ..., EnX caractéristiques des isotopes naturels de l'élément X. Dans le cadre de l'invention, il existe un risque d'interférence des raies de capture qui correspondent à des éléments recherchés avec des raies de capture d'éléments autres que les éléments recherchés ou avec des raies de type inélastique dues à l'émission résiduelle du générateur G1 entre les impulsions intenses délivrées par le générateur G2. L'algorithme Al a pour but de sélectionner, dans une liste de raies associées à chaque élément X d'intérêt, celles présentant une intensité suffisante et un risque d'interférence faible. L'algorithme Al comprend une étape ES1 de traitement du spectre Si, une étape N de normalisation du spectre Sc, une étape ES2 de traitement du spectre So normalisé, une étape de soustraction 16 qui soustrait le résultat de l'étape de traitement ES2 du résultat de l'étape de traitement ES1, une étape 17 de validation de raie et une étape 18 de détermination de présence d'un élément. L'étape de traitement ES1 comprend, successivement, une étape 10 de recherche de pics, une étape 11 d'extraction d'aires nettes et une étape 12 de validation de détection de raies. Chacune de ces étapes est une étape connue en soi. L'étape 10 de recherche de pics est suivie, pour chaque pic détecté, d'une étape 11 d'extraction de l'aire nette A du pic. La figure 6 représente un pic de spectre gamma. L'aire nette A du pic représenté en figure 6 est la surface du pic qui dépasse le niveau de bruit de fond Compton Bc. Une raie Exi d'un élément X est considérée comme détectée dans le spectre si l'aire nette As1 (Exi) de la raie Exi vérifie, par exemple, l'équation : [Asi (Exi) ] /Asi (Exi) Le facteur F peut également être obtenu de façon empirique en mesurant, lors d'une phase d'étalonnage réalisée dans des conditions similaires à celles d'une mesure d'inspection, pendant les durées T2 et T4, les intensités de raies gamma de type inélastique. Les raies gamma utilisées pour le calcul du facteur F sont alors précisément choisies pour leur intensité, pour leur énergie et parce qu'elles ne présentent pas d'interférence connue avec des rayonnements de capture. Le spectre gamma normalisé So(n) est alors formé en multipliant le nombre de coups présent dans chaque canal d'énergie du spectre gamma So par le facteur F. L'étape de traitement ES2 succède à l'étape N de normalisation. L'étape de traitement ES2 comprend, de façon connue en soi, une recherche des pics présents dans le spectre normalisé So(n) (étape 13), une extraction des aires nettes Aso (nor) (EXi) des pics du spectre normalisé (étape 14) et une étape de validation de détection de raie (étape 15). De chaque aire nette de capture An (EXi) de raie EXi délivrée suite de l'étape 12, on soustrait, à l'étape 16, l'aire nette normalisée Aso (nor) (Ex�) délivrée suite à l'étape 15. Dans le cas où il n'y a pas de pic détecté à l'étape 13 (ce qui est le cas le plus fréquent), l'aire nette normalisée Aso (nor) (EXi) relative à l'énergie EXi est définie comme étant égale à zéro. A l'issue de l'étape 16, on obtient alors une aire nette Ans1(Exi) dépourvue des contributions de type inélastique, à savoir : AnS1 (Exi) ù An (Exi) ù ASO (nor) (Exi) La raie de capture d'un élément X est confirmée à l'étape 17, si la relation suivante est, par exemple, vérifiée : a [Anse (EXi) ] /Anse (EXi) REVENDICATIONS1. Procédé de détection de présence d'un élément chimique (X) dans un objet par émission de neutrons sur l'objet, caractérisé en ce que l'émission de neutrons sur l'objet est constituée, d'une part, par une émission continue de neutrons provenant d'un générateur neutronique à particule associée (G1) et, d'autre part, par une émission d'impulsions neutroniques qui se superposent à l'émission continue de neutrons, les impulsions neutroniques provenant d'un générateur neutronique pulsé (G2) qui génère des impulsions neutroniques de durée d'impulsion T2, deux impulsions neutroniques successives étant séparées d'une durée T4, les émissions continue et pulsée de neutrons sur l'objet produisant un rayonnement gamma de capture et un rayonnement gamma de type inélastique. 2. Procédé selon la revendication 1, caractérisé en ce qu'il comprend : - une détection, à l'aide de la technique de la particule associée, pendant les durées T4 qui séparent les impulsions de durées T2, du rayonnement gamma de type inélastique et du rayonnement gamma de capture, û une formation (4, 5), dans une plage de temps de vol neutronique relative à l'objet, d'un premier spectre gamma (Sb) à partir du rayonnement gamma de type inélastique et du rayonnement gamma de capture détectés à l'étape de détection à l'aide de la technique de la particule associée, une formation (e1), dans une plage temporelle disjointe de la plage de temps de vol neutronique relative à l'objet ou de toute autre plage de temps de vol neutronique relative à tout autre objet susceptible d'environner l'objet, d'un spectre de bruit de fond aléatoire (Sa) relatif aux coïncidences fortuites, à partir du rayonnement gamma de type inélastique et du rayonnement gamma de capture détectés à l'étape de détection à l'aide de la technique de la particule associée (4), - une formation (e2), à partir du spectre de bruit de fond aléatoire relatif aux coïncidences fortuites (Sa), d'un spectre de bruit de fond aléatoire normalisé (Sa(n)) sur la plage temporelle relative à l'objet, - une soustraction (6) du premier spectre gamma (Sb) et du spectre de bruit de fond aléatoire normalisé (Sa(n)) pour former un deuxième spectre gamma (Sn), et - une étape de recherche (8) de la présence d'un élément chimique dans l'objet à partir du deuxième spectre gamma (Sn) et de spectres d'éléments chimiques de référence susceptibles d'être détectés (Bs) pour délivrer une première donnée (dix) apte à traduire la présence ou l'absence de l'élément chimique (X) dans l'objet. 3. Procédé selon la revendication 2, caractérisé en ce qu'il comprend, en outre: - une détection, à l'aide de la technique d'interrogation neutronique pulsée, du rayonnementgamma de capture et du rayonnement gamma de type inélastique, d'une part pendant les durées d'impulsions T2 et, d'autre part, pendant les durées T4 qui séparent les impulsions de durée T2. - une formation (2) d'un troisième spectre gamma (Si), à partir du rayonnement gamma de capture et du rayonnement gamma de type inélastique détectés, pendant les durées T4, à l'aide de la technique d'interrogation neutronique pulsée, - une formation (ES1), à partir du troisième spectre gamma (Si), d'une première donnée (As1 (E X)) représentative d'au moins une raie de capture (Evv, E2x, ..., E,x) de l'élément chimique (X), - une formation (1) d'un quatrième spectre 15 gamma (So) constitué d'une pluralité de canaux d'énergie, à partir du rayonnement de type inélastique détecté pendant les durées d'impulsion T2, à l'aide de la technique d'interrogation neutronique pulsée, - une formation (N) d'un spectre gamma 20 normalisé (S(nor)) à partir du quatrième spectre gamma (So), la formation du spectre gamma normalisé (So(nor) ) étant effectuée en multipliant le contenu de chaque canal d'énergie du quatrième spectre gamma par un facteur de normalisation (F) sensiblement proportionnel 25 au nombre de neutrons émis pendant une durée T4 et sensiblement inversement proportionnel au nombre de neutrons émis pendant une durée T2, une formation (ES2), à partir du spectre gamma normalisé (So (nor)) , d'une deuxième donnée 30 (Aso (nor) (Ex)) représentative d'un bruit de fond de rayonnement gamma de type inélastique, 10- une étape de soustraction (16) de la première donnée (As1(EiX)) et de la deuxième donnée (Asonor (EiX)) , - une étape de validation (17) de la présence ou de l'absence de l'élément chimique (X) dans l'objet en fonction du résultat de l'étape de soustraction, et si la présence de l'élément chimique (X) est validée dans l'objet, une étape supplémentaire (18) qui délivre une deuxième donnée (djX) apte à confirmer ou infirmer la validation de la présence de l'élément chimique (X) sur la base du résultat de l'étape de validation (17) et de données de poids (W (E1X) , W (E2X) , W(EnX)) relatives aux différentes raies gamma de rayonnement de capture de l'élément chimique (X), ces poids étant attribués en fonction d'une probabilité de détection et d'un risque d'interférence associés à chaque raie. 4. Procédé selon la revendication 3, dans lequel le facteur de normalisation (F) est donné par l'équation : F = (I4xT4)/(I2xT2)x[(1-TM(T4))/(1-TM(T2))], où TM(T4) et TM(T2) sont, respectivement, un temps mort d'une électronique de traitement (E) qui traite les signaux électriques issus de la détection des rayonnements gamma pendant les durées T4 et un temps mort d'une électronique de traitement (E) qui traite les signaux électriques issus de la détection des rayonnements gamma détectés pendant les durées T2. 5. Procédé selon la revendication 3 ou 4, dans lequel une étape de décision (9) délivre un signal d'alerte (S(alerte)) si la première donnée (dix) ou la deuxième donnée (djx) confirme la présence de l'élément chimique. 6. Procédé selon la revendication 5, dans lequel, dès lors que la présence de plusieurs éléments chimiques est confirmée, il comprend une étape supplémentaire de calcul de proportions relatives de tout ou partie des éléments chimiques dont la présence est confirmée. 7. Procédé selon l'une quelconque des revendications 2 à 6, dans lequel la détection, à l'aide de la technique de la particule associée, du rayonnement gamma de capture et du rayonnement gamma de type inélastique pendant les durées T4 qui séparent les impulsions de durée T2 est effectuée par un détecteur (Dl) préférentiellement adapté pour la détection de rayonnement gamma de type inélastique. 8. Procédé selon l'une quelconque des revendications 3 à 7, dans lequel la détection, à l'aide de la technique d'interrogation neutronique pulsée, du rayonnement gamma de capture et du rayonnement gamma de type inélastique, d'une part, pendant les durées d'impulsion T2 et, d'autre part, pendant les durées T4 qui séparent les impulsions de durée T2 est effectuée par un détecteur (D2)préférentiellement adapté pour la détection de rayonnement gamma de capture.