FR 3075980 A1 20190628 FR 1762989 A 20171222 Procédé d'analyse à l'aide d'un détecteur de particules alpha La présente invention concerne un procédé d'analyse, à l'aide d'un détecteur de particules alpha, d'un échantillon comprenant au moins un radionucléide émetteur d'une pluralité de particules alpha. Les particules alpha ou rayons alpha sont une forme de rayonnement émis par des noyaux instables de grande masse atomique. Elles sont constituées de deux protons et de deux neutrons combinés en une particule identique au noyau d'hélium 4 (hélion). Les particules alpha sont émises par des noyaux radioactifs (radionucléides), comme l'uranium ou le radium, au cours d'un processus nommé « désintégration alpha ». Durant ce phénomène, le radionucléide éjecte une particule alpha et se transforme en un noyau de nombre de masse inférieur. L'analyse du rayonnement alpha est par exemple utilisée dans le domaine de la géologie pour caractériser la présence d'uranium 238 ou de thorium 232 et de leurs descendants. On trouve également des applications dans le domaine médical et plus généralement dans tous les domaines des sciences du vivant ou de l'environnement, par exemple pour imager des tissus animaux ou végétaux, ou détecter la présence de radionucléides tels que le polonium 210. La spectrométrie alpha est une technique d'analyse utilisée pour déterminer la nature du radioélément émetteur alpha avec une grande précision. Ces mesures sont par exemple réalisées à l'aide d'un détecteur à diodes de silicium ou à scintillateur liquide. Cependant, les mesures réalisées ne permettent pas d'obtenir une distribution spatiale de l'activité des radioéléments présents dans l'échantillon. L'autoradiographie est une autre technique d'analyse d'un échantillon qui permet de caractériser, par exemple à l'aide d'un écran phosphore, la distribution spatiale des radionucléides émetteurs alpha sans toutefois permettre de déterminer la nature des radionucléides. Ainsi, la caractérisation de la nature et de la distribution spatiale des radionucléides d'un échantillon fait appel à deux techniques d'analyse distinctes et à des détecteurs de nature différente. Par conséquent, la caractérisation de l'échantillon est longue et fastidieuse. De plus, les résolutions des deux méthodes sont souvent différentes l'une de l'autre car elles dépendent également de la nature du détecteur utilisé, ce qui rend difficile l'intégration et l'interprétation conjointe de ces deux types de mesures. L'invention a pour but de pallier ces inconvénients en proposant un procédé d'analyse à l'aide d'un détecteur de particules alpha qui permette une caractérisation 2 simple et rapide de la nature et de la distribution spatiale des radionucléides présents dans l'échantillon. Pour ce faire, l'invention a pour objet un procédé d'analyse à l'aide d'un détecteur de particules alpha, d'un échantillon comprenant au moins un radionucléide émetteur 5 d'une pluralité de particules alpha, le détecteur comprenant un milieu de détection et une pluralité de cellules de mesure adaptées pour mesurer au moins un signal incident généré par une interaction de la particule alpha avec ledit milieu de détection, le détecteur étant conçu pour fournir une image autoradiographique dudit échantillon, le procédé comprenant les étapes suivantes : 10 - El) dans laquelle on détermine une énergie initiale de chaque particule alpha, ladite étape de détermination (El) ) comprenant une répétition des sous-étapes suivantes pour chaque particule alpha : - D1) dans laquelle on détermine une position d'une première interaction de la particule alpha avec le milieu de détection du détecteur en utilisant au 15 moins une cellule de mesure, - D2) dans laquelle on détermine une énergie déposée par la particule alpha à l'intérieur du détecteur en utilisant le signal incident mesuré, - D3) dans laquelle on détermine, en utilisant l'énergie déposée déterminée à la sous-étape (D2), l'énergie initiale de la particule alpha avant la 20 première interaction de la particule alpha avec le milieu de détection du détecteur, - E2) dans laquelle pour au moins une zone de l'image autoradiographique, on construit un spectre en énergie en utilisant les énergies initiales déterminées des particules alpha correspondantes. 25 Ainsi, le procédé selon l'invention permet de caractériser la nature et la distribution spatiale des radionucléides de l'échantillon en utilisant un unique détecteur, simplement et rapidement. Dans le cadre de l'invention, le détecteur peut être notamment choisi parmi un détecteur gazeux, un détecteur à semi-conducteur ou un détecteur scintillateur. Suivant des modes particuliers de réalisation, le procédé selon l'invention 30 comporte l'une ou plusieurs des caractéristiques suivantes, prises isolément ou selon toutes les combinaisons techniquement possible : - la sous-étape (D1) comprend au moins une opération d'amplification (01, 02) d'une interaction de la particule alpha avec le milieu de détection du détecteur à l'aide d'au moins un champ électrique, 3 - avant l'étape (El), on réalise une étape de calibration (E0) de ladite énergie initiale de la particule alpha en utilisant une source de calibration, - le détecteur de particules alpha est un détecteur gazeux et le signal incident est une charge en mouvement, 5 - la position de la première interaction de la particule alpha avec le milieu de détection du détecteur est déterminée à partir d'une répartition d'une pluralité de signaux incidents mesurés par les cellules de mesure, - la sous-étape (D2) comprend une opération de mesure (04) de charges générées successivement le long de la trajectoire de la particule alpha à l'intérieur du 10 détecteur, et une opération de mesure (06) de la longueur de la trajectoire parcourue par la particule alpha à l'intérieur du détecteur de particules alpha, l'énergie initiale de la particule alpha étant déterminée en utilisant l'énergie déposée déterminée et la longueur de la trajectoire mesurée, 15 - le procédé comprend préalablement à l'opération de mesure (06) de la longueur de la trajectoire de la particule alpha, une opération d'identification (05) des particules alpha ayant traversé le détecteur de particules alpha de part en part, étant entendu que seules les particules alpha ayant traversé le détecteur de part en part sont considérées à l'étape (E2) de construction des spectres en énergie, 20 - l'opération d'identification (05) comprend un calcul d'un paramètre géométrique d'une distribution des charges mesurées dans le détecteur de particules alpha selon une direction du détecteur de particules alpha, - le paramètre géométrique est un coefficient d'aplatissement de la distribution des charges mesurées, 25 - l'opération d'identification (05) comprend l'application d'un seuil sur le coefficient d'aplatissement de la distribution des charges mesurées et la sélection des particules alpha ayant un coefficient d'aplatissement inférieur audit seuil, - l'énergie initiale de la particule alpha est déterminée en utilisant un ratio entre l'énergie déposée par la particule alpha et la longueur de la trajectoire parcourue par 30 ladite particule alpha, - l'énergie initiale de la particule alpha est déterminée en reliant ladite énergie initiale au ratio entre l'énergie déposée par la particule alpha et la longueur de la trajectoire parcourue par ladite particule alpha par une relation de type exponentiel, - la relation de type exponentiel répond à la formule : 4 - Ede?, étant l'énergie déposée par la particule alpha (3B) considérée, - d étant la longueur de la trajectoire parcourue par la particule alpha (3B), - a et b étant des coefficients, - le détecteur est un détecteur à semi-conducteur ou un détecteur scintillateur et le 5 signal incident mesuré est respectivement une charge en mouvement ou un photon. D'autres aspects et avantages de l'invention apparaitront à la lecture de la description qui suit, donnée à titre d'exemple et faite en référence aux dessins annexés parmi lesquels : - la figure 1 est une représentation schématique d'un détecteur gazeux utilisé pour 10 mettre en oeuvre le procédé selon un mode de réalisation ; - la figure 2 est un diagramme représentant les différentes étapes, sous-étapes et opérations d'un procédé selon l'invention ; - la figure 3 est une représentation graphique des charges déposées par une particule alpha dans le détecteur de la figure 1 selon l'axe X de l'anode du détecteur ; 15 - la figure 4 est une représentation graphique de la distribution de la valeur du kurtosis utilisée pour sélectionner les particules alpha ayant traversées le détecteur de part en part ; - la figure 5 est une représentation graphique de la distribution du ratio de l'énergie déposée sur la longueur de la trajectoire parcourue par la particule alpha pour une source 20 de calibration au radium 226 ; - la figure 6 est une représentation graphique du spectre en énergie de la source de calibration au radium 226, après application de la calibration ; - la figure 7 est une représentation graphique d'un spectre en énergie obtenu avec le procédé selon l'invention. 25 Le procédé selon l'invention concerne un procédé d'analyse à l'aide d'un détecteur 1 de particules alpha 3A, 3B, d'un échantillon 5 comprenant au moins un radionucléide émetteur d'une pluralité de particules alpha 3A, 3B (figure 1). Le détecteur 1 est conçu pour fournir une image autoradiographique de l'échantillon 5. 30 Une image autoradiographique est une image de a distribution spatiale des radionucléides émetteurs alpha au sein de l'échantillon 5. Selon un mode de réalisation préféré, le détecteur 1 de particules alpha 3A, 3B est un détecteur gazeux tel que représenté sur la figure 1. Le détecteur 1 est sensiblement similaire au détecteur gazeux décrit dans le 35 document WO 2011/039473 Al. Dans ce qui suit, seules les caractéristiques du détecteur 5 1 nécessaires à la compréhension du procédé d'analyse selon l'invention sont précisées en se référant à la figure 1. Le détecteur 1 comprend une enceinte 7 définie par deux parois principales 9, 11 opposées et parallèles entre elles et des parois latérales 12. 5 Les parois principales 9, 11 ont par exemple une surface comprise entre 1 cm x 1 cm et 30 cm x 30 cm, par exemple 23 cm x 23 cm. L'enceinte 7 contient un milieu de détection adapté pour émettre des électrons primaires sous l'effet du rayonnement alpha émis par l'échantillon 5. Le milieu de détection est typiquement constitué d'un mélange gazeux circulant 10 dans l'enceinte 7 entre une entrée et une sortie (non représentées). Le mélange gazeux est par exemple un mélange entre un gaz rare, par exemple du xénon, et un gaz polyatomique, par exemple du CH4. Le mélange gazeux est avantageusement à une pression comprise entre 0,1 bar et 5 bars, par exemple 1 bar. 15 L'enceinte 7 a typiquement une épaisseur comprise entre 0,5 cm et 20 cm, par exemple 1 cm. Le détecteur 1 comprend une anode 13, une cathode 15, une première électrode 17 et une deuxième électrode 19. Typiquement, l'anode 13, la cathode 15, la première électrode 17 et la deuxième 20 électrode 19 sont parallèles entre elles et parallèles aux deux parois principales 9, 11 de l'enceinte 7. La cathode 15 est par exemple un adhésif conducteur, par exemple un adhésif de cuivre, collé sur une face d'une lame de verre 21, l'échantillon 5 étant disposé sur une face opposée de la lame de verre 21. 25 La cathode 15 et la première électrode 17 définissent entre elles un premier espace d'amplification 23. La première électrode 17 et la deuxième électrode 19 définissent entre elles un espace de dérive 25. La deuxième électrode 19 et l'anode 13 définissent entre elles un deuxième 30 espace d'amplification 27. Les phénomènes d'amplification et de dérive sont décrits dans la suite de la description. La première électrode 17 et la deuxième électrode 19 sont par exemple des microgrilles de type MICROMEGAS telles que décrites dans WO 2011/039473 Al. 6 Avantageusement, des moyens de polarisation (non représentés) sont connectés à la cathode 15, à l'anode 13, à la première électrode 17 et à la deuxième électrode 19. Ils permettent de porter la cathode 15 à un potentiel V1, l'anode 13 à un potentiel V2, la première électrode 17 à un potentiel V3 et la deuxième électrode 19 à un potentiel V4. 5 Typiquement, ces potentiels vérifient la relation V2> V4> V3> V1. Les moyens de polarisation permettent ainsi de créer des champs électriques E1, E2, E3 respectivement dans le premier espace d'amplification 23, l'espace de dérive 25 et le deuxième espace d'amplification 27. Les champs électriques El et E3 dans le premier espace d'amplification 23 et le 10 deuxième espace d'amplification 27 sont typiquement supérieurs à 3 kV/cm. Le champ électrique E2 dans l'espace de dérive 25 est par exemple inférieur à 3 kV/cm. Les potentiels V1, V2, V3 et V4 sont avantageusement inférieurs à quelques kV par exemple 3 kV. 15 Le détecteur 1 comprend une pluralité de cellules de mesure 29 adaptées pour mesurer au moins un signal incident généré par une interaction de la particule alpha 3A, 3B avec le milieu de détection du détecteur 1. Dans le cas du détecteur 1 gazeux décrit ci-dessus, les cellules de mesure 29 sont formées par des anodes élémentaires ou pavés 31 de l'anode 13. 20 L'anode 13 est par exemple du type décrite dans le document WO 2011/039473 Al. Les pavés 31 sont typiquement séparés les uns des autres par des espaces 33 et formant un réseau bidimensionnel de mesure en damier dont les rangées sont alignées selon des axes de coordonnées X et Y perpendiculaires. 25 Chaque pavé 31 forme par exemple un carré de moins d'un millimètre de côté, par exemple de 650 iim. Les pavés 31 sont alternativement affectés à la lecture de l'une ou l'autre des coordonnées X et Y. Deux pavés 31 voisins ne mesurent pas la même position selon la même coordonnée. L'espace 33 entre les pavés 31 est le plus faible possible, tout en permettant d'isoler chaque pavé 31 du pavé 31 voisin. Par exemple, l'espace 33 30 est inférieur ou égal à 100 iim. La structure de chaque pavé 31 est par exemple similaire à celle décrite dans le document WO 2011/039473 Al. Le signal mesuré par les pavés 31 est typiquement une charge électronique incidente. 7 Les pavés 31 permettent de convertir une charge électronique incidente en un courant électrique proportionnel à ladite charge. L'échantillon 5 peut comporter un seul type de radionucléide ou un mélange de radionucléides différents. 5 Les radionucléides sont typiquement d'origine naturelle ou artificielle. Par exemple, les radionucléides sont l'uranium 238 et/ou le thorium 232 et leurs descendants. En variante, les radionucléides sont des actinides artificiels comme par exemple le plutonium. 10 L'échantillon 5 est par exemple un échantillon de roche, de tissu végétal, animal ou humain. L'échantillon 5 est typiquement une couche mince ayant une épaisseur comprise entre 5 lm et 100 11m, par exemple 20 L'échantillon 5 a par exemple des dimensions latérales avantageusement 15 inférieures au champ de vue du détecteur 1, i.e. les dimensions de l'anode 13. Le champ de vue du détecteur 1 correspond à l'angle solide au travers duquel le détecteur 1 est sensible au rayonnement des particules alpha 3A, 3B. L'échantillon 5 est typiquement déposé sur la lame de verre 21, de manière à maintenir l'échantillon sensiblement parallèle à l'anode 13. 20 La lame de verre 21 a typiquement une épaisseur d'un millimètre. En variante, l'échantillon 5 est déposé sur une pièce métallique. Avantageusement, la surface de l'échantillon 5 est sensiblement lisse et sensiblement parallèle à l'électrode 17. Avantageusement, la préparation de l'échantillon 5 comprend une étape de 25 métallisation de la surface de l'échantillon 5, par exemple à l'or/palladium, permettant d'assurer une polarisation homogène de la surface de l'échantillon 5. Cette préparation est typiquement identique à celle, couramment utilisée, d'échantillons pour la microscopie électronique à balayage. Un procédé d'analyse (figure 2) mis en oeuvre au moyen d'un détecteur 1 gazeux 30 tel que décrit plus haut va maintenant être détaillé ci-dessous. En référence à la figure 1, une particule alpha 3A, 3B est émise par l'échantillon 5 au niveau de la face principale 9 de l'enceinte 7 du détecteur 1 gazeux. La particule alpha 3A, 3B traverse l'enceinte du détecteur 1 et interagit avec le milieu gazeux en générant des électrons primaires qui diffusent et se multiplient dans l'enceinte 7 du détecteur 1 35 depuis la cathode 15 jusqu'à l'anode 13. 8 Le procédé comprend une étape de détermination El d'une énergie initiale de chaque particule alpha 3A, 3B émise. L'étape de détermination El comprend une répétition de trois sous-étapes D1, D2 et D3, pour chaque particule alpha 3A, 3B. La première sous-étape D1 comprend une détermination d'une position 5 correspondant à une première interaction de la particule alpha 3A, 3B avec le milieu de détection du détecteur 1 de particules alpha 3A, 3B en utilisant au moins une cellule de mesure 29. La première sous-étape D1 permet la localisation du radionucléide dans l'échantillon 5 qui est à l'origine de l'émission de la particule alpha 3A, 3B. 10 Avantageusement, la première sous-étape D1 comprend au moins une opération d'amplification de l'interaction de la particule alpha 3A, 3B avec le détecteur 1 de particules alpha 3A, 3B, à l'aide d'au moins un champ électrique. Selon le mode de réalisation préféré dans lequel le détecteur 1 de particules alpha 3A, 3B est un détecteur gazeux, la première sous-étape D1 comprend l'amplification de 15 l'interaction de la particule alpha 3A, 3B avec les atomes de gaz présents dans l'enceinte 7 du détecteur 1. Ainsi, dans le mode de réalisation préféré, la première sous-étape D1 comprend une première opération d'amplification 01 dans le premier espace d'amplification 23. Le champ électrique El qui règne à l'intérieur du premier espace d'amplification 23 20 accélère les électrons créés par l'interaction entre la particule alpha 3A, 3B et le milieu gazeux. Chaque électron entre en collision avec un atome de gaz et l'ionise avec la libération d'un nouvel électron. Les deux électrons ionisent à nouveau deux atomes de gaz qui libèrent à leur tour un électron chacun. Ainsi, par effet d'avalanche, de plus en plus d'électrons sont créés entre la cathode 15 et la première électrode 17. 25 Avantageusement, la première opération d'amplification 01 est réalisée au plus proche de l'échantillon 5. Ceci permet de maximiser l'effet d'avalanche électronique et d'obtenir un pic de charge électronique important correspondant aux charges matérialisées par un trapèze 36 sur la figure 1. Comme visible sur la figure 1, la particule alpha 3A, 3B interagit successivement 30 avec plusieurs atomes de gaz dans le premier espace d'amplification 23. Chacune de ces interactions est amplifiée par effet d'avalanche électronique et participe au pic de charge matérialisé par un trapèze 36 sur la figure 1. Les électrons 35A, 35B et 35C générés par le phénomène d'avalanche électronique dans le premier espace d'amplification 23 diffusent ensuite dans l'espace de 9 dérive 25. L'espace de dérive 25 permet de transférer les charges électroniques 39 et 37 jusqu'à l'anode 13. Avantageusement, la première sous-étape D1 comporte une deuxième opération d'amplification 02 dans le deuxième espace d'amplification 27 par effet d'avalanche 5 électronique, de manière similaire à ce qui a été décrit plus haut. La sous-étape D1 comprend alors une opération de détermination 03 de la position moyenne (XO, YO), selon l'axe X et l'axe Y de l'anode 13, de la première interaction entre la particule alpha 3A, 3B et le milieu gazeux, à l'aide des pavés 31 de l'anode 13. 10 La position d'entrée de la particule alpha 3A, 3B dans l'espace d'amplification 23 est typiquement déterminée en utilisant les pavés 31, respectivement affectés à la mesure selon l'axe X et selon l'axe Y, recevant la charge électronique en sortie de l'espace d'amplification 27. La position de la première interaction de la particule alpha 3A, 3B avec le milieu de 15 détection du détecteur 1 est alors avantageusement déterminée à partir d'une répartition des signaux incidents mesurés par les cellules de mesure 29. Dans l'exemple présenté ci-dessus, la position moyenne (XO, YO) est déterminée à partir de la répartition des charges électroniques incidentes mesurées par les pavés 31 de l'anode 13. 20 La position moyenne ainsi déterminée correspond alors à la position au sein de l'échantillon 5 du radionucléide à l'origine de l'émission de la particule alpha 3A, 3B. La mise en oeuvre de la sous-étape D1 suppose que la trajectoire des électrons créés par l'interaction entre la particule alpha 3A, 3B et par le phénomène d'avalanche électronique s'effectue en ligne droite selon la direction des champs électriques dans les 25 espaces d'amplification 23, 27 et de dérive 25, sans déviation d'angle importante. Après avoir traversé le premier espace d'amplification 23, la particule alpha 3A, 3B poursuit son parcours dans l'espace de dérive 25. La particule alpha 3A, 3B interagit alors successivement avec des atomes du milieu gazeux de l'espace de dérive et génère des électrons 39 qui diffusent dans l'espace de dérive 25 jusqu'au deuxième espace 30 d'amplification 27 où les charges électroniques 41 sont amplifiées par effet d'avalanche électronique de manière similaire à ce qui a été décrit précédemment puis détectées par les pavés 31 de l'anode 13. Le procédé comprend une sous-étape de détermination D2 de l'énergie déposée par la particule alpha 3A, 3B à l'intérieur du détecteur 1 de particules alpha 3A, 3B en 35 utilisant le signal incident mesuré. 10 Typiquement, la sous-étape D2 comprend une opération 04 de mesure des charges 41, 36 générées successivement le long d'une trajectoire de la particule alpha 3A, 3B à l'intérieur du détecteur 1. Les charges électroniques 41 et 36 sont typiquement diffusées dans l'espace de 5 dérive 25 et amplifiées dans les deux espaces d'amplification 23 et 27. Ces mesures sont effectuées en utilisant les pavés 31 de l'anode. La figure 3 présente un exemple de distribution de charges électroniques mesurées par l'anode 13 selon l'axe X. La position du pic de charge électronique créé dans le premier espace 10 d'amplification 23 du détecteur 1, correspondant à la première interaction de la particule alpha 3A, 3B avec le milieu de détection du détecteur 1, est localisée en X=0. L'énergie déposée est reliée au nombre d'électrons primaires par la relation : Edep = W X nb primaires (1) W est l'énergie moyenne d'ionisation dans le gaz et nb primaires est le nombre 15 d'électrons primaires. L'énergie moyenne d'ionisation est dépendante du mélange gazeux. La valeur moyenne typiquement retenue pour des mélanges gazeux vaut 30 eV. La charge totale Ctotalei mesurée le long de l'axe i de l'anode 13 s'écrit : Ctotalei = niaia2e + n2a2e (2) n1 est un nombre de charges primaires générées dans le premier espace d'amplification 23 par l'interaction de la particule alpha 3A, 3B avec le milieu de détection. 20 n2 est un nombre de charges primaires générées dans l'espace de dérive 25 par l'interaction de la particule alpha 3A, 3B avec le milieu de détection. al est un gain d'amplification associé à l'espace d'amplification 23. a2 est un gain d'amplification associé à l'espace d'amplification 27. e est la charge d'un électron. 25 i=x ou i=y, selon que les calculs sont effectués avec la distribution de charges mesurée selon l'axe X ou l'axe Y de l'anode 13. En référence à la figure 3, la charge associée à la première partie de la distribution de charges Clis'écrit : = n2a2e (3) La charge associée à la deuxième partie de la distribution de charges Cpici 30 correspondant à la charge mesurée au niveau du pic dû aux charges déposées dans le premier espace d'amplification 23 ayant subi l'amplification du premier et du deuxième espace d'amplification 23 et 27 s'écrit : Cpici = n1a1a2e (4) 11 La C charge - - extrapolationi déposée dans le premier espace d'amplification 23 qui serait mesurée si le gain al dans le premier espace d'amplification 23 était unitaire, s'écrit : Cextrapolationi = n1a2e (5) Cptei peut alors s'écrire : Cpici = Cextrapolationial (6) 5 Si al = 1 (gain unitaire) alors Cptei s'écrit : Cpici = Cextrapolationi (7) Ctotalei peut donc s'écrire sous la forme : Ctotalei = Cpici + C1. (8) Par définition, la charge déposée par la particule alpha 3A, 3B dans le milieu de détection Cdepi peut s'écrire sous la forme : (n1a2e + n2a2e) (C extrapolationi C1i) Cdepi = nie + n2e = = (9) a2 a2 Donc, on obtient : 10 C1. = a2Cdepi Cextrapolationi (10) et Ctotalei = Cpici+Cdepia2 - Cextrapolationi (11) La charge déposée Cdepi s'écrit alors : Cdepi = (Ctotalei Cextrapolationi Cpici) a2 (12) Les charges Ctotalei et Cptei sont typiquement calculées en sommant les valeurs de chaque colonne verticale correspondante de l'histogramme de distribution de charge tel 15 que par exemple représenté sur la figure 3. La charge déposée totale Crde'totale par la particule alpha 3A, 3B est obtenue en sommant la charge déposée Cdepx selon l'axe X et la charge déposée Cdepy selon l'axe Y. La charge déposée totale Ce totale s'écrit également : = e x nbprimaires (13) CdePtotaie On en déduit donc l'équation suivante pour l'énergie déposée totale par la 20 particule alpha 3A, 3B dans le milieu de détection du détecteur 1 : CdePtotaie X W EdePtotaie = (14) Avantageusement, la sous-étape D2 comprend une opération d'identification 05 des particules alpha 3B traversant le détecteur 1 de particules alpha 3A, 3B de part en part. 12 En effet, comme visible sur la figure 1, la nature et la géométrie du détecteur 1 de particules alpha 3A, 3B induit deux types d'évènements : la particule alpha 3B traverse le détecteur 1 de part en part ou la particule alpha 3A est calorimétrisée dans le gaz et s'arrête à l'intérieur de l'enceinte 7 du détecteur 1. 5 Typiquement, l'opération d'identification 05 comprend un calcul d'un paramètre géométrique de la distribution de l'énergie déposée par la particule alpha 3A, 3B dans le détecteur 1 de particules alpha 3A, 3B selon une direction du détecteur 1 de particules alpha 3A, 3B. Par exemple, le paramètre géométrique est un coefficient d'aplatissement de la 10 distribution de l'énergie déposée selon l'axe X ou l'axe Y, encore appelé kurtosis. En effet, une particule alpha 3B qui traverse le détecteur 1 génère une distribution de charges selon l'axe X ou Y de l'anode 13, qui est plus « plate » que celle d'une particule alpha 3A qui est stoppée à l'intérieur de l'enceinte 7 du détecteur 1. Cette caractéristique est liée au fait que le dépôt d'énergie par unité de longueur est plus 15 important en fin de parcours dans le gaz. Le kurtosis correspond à une mesure du coefficient d'aplatissement d'une fonction ou d'une distribution. Plus cette dernière est plate et plus la valeur du kurtosis est faible. Avantageusement, le kurtosis est calculé sur une plage de mesures excluant le pic de charge dû aux charges amplifiées dans le premier espace d'amplification 23. En effet, 20 les charges correspondant au pic, par exemple sur la figure 3 autour de la position X = 0, subissent deux amplifications successives dans les espaces d'amplification 23 et 27. Le kurtosis est ainsi calculé seulement pour les charges qui subissent seulement l'amplification de l'espace d'amplification 27. Par exemple, en se référant à la figure 3, le kurtosis est calculé en tenant compte 25 de la distribution de charges comprises entre X=-1,3 cm et X=-0,15 cm, en excluant le pic de charge centré sur X=0 cm. La figure 4 présente un exemple de distribution de la valeur du kurtosis pour une pluralité de particules alpha émises par l'échantillon, dont l'énergie déposée est mesurée par l'anode 13. 30 Deux distributions 43, 45 de tendance gaussienne sont visibles. La première distribution 43 centrée sur -1,2 correspond aux particules alpha 3B qui ont traversé le détecteur 1 de part en part. La deuxième distribution 45 centrée sur -1,08 correspond aux particules alpha 3A qui ont été calorimétrisées à l'intérieur de l'enceinte 7 du détecteur 1. 13 Comme visible sur la figure 4, les deux distributions 43, 45 se superposent et par conséquent pour certaines particules alpha 3A, 3B, il n'est pas possible à partir de la mesure du kurtosis de déterminer si les particules 3A, 3B ont traversé le détecteur 1 ou si elles ont été calorimétrisées à l'intérieur du détecteur 1. 5 Ainsi, typiquement, l'opération d'identification 05 des particules alpha 3A, 3B comprend l'application d'un seuil sur la valeur du kurtosis et la sélection des particules alpha 3A, 3B dont la valeur du kurtosis est inférieure ou égale au seuil. Dans l'exemple de la figure 4, le seuil est par exemple fixé à -1,12. Ainsi, seules les particules alpha 3A, 3B présentant une valeur de kurtosis inférieure ou égale à -1,12 10 seront sélectionnées pour la suite du procédé. La sous-étape D2 comprend alors typiquement une opération de mesure 06 de la longueur de la trajectoire parcourue par la particule alpha 3B à l'intérieur du détecteur 1 de particules alpha 3A, 3B, seules les particules alpha 3B traversant le détecteur 1 de part en part étant considérées durant cette opération. 15 Ainsi, pour chacune des distributions de charges mesurées selon l'axe X et l'axe Y de l'anode 13, comme par exemple la distribution présentée sur la figure 3, la plus petite coordonnée xmm, ymm et la plus grande coordonnée xmax, ymax, pour lesquelles au moins une charge a été mesurée par l'anode 13 sont déterminées. Les distances parcourues dx et dy, respectivement selon l'axe X et l'axe Y sont 20 alors calculées selon les relations : dx = xmax - xmin et dy = Ymax Ymin (15) La longueur de la trajectoire d est alors calculée en utilisant la relation : d = dx2 dy2 dz2 (16) La longueur dz correspond à l'épaisseur du gaz traversée. Pour les particules 25 alpha 3B ayant traversées le détecteur de part en part, cette épaisseur est connue car elle dépend directement des dimensions du détecteur 1. Par exemple, pour le détecteur 1 décrit plus haut, l'épaisseur dz vaut 1 cm. Le procédé comprend alors une sous-étape de détermination D3, en utilisant l'énergie déposée déterminée, d'une énergie initiale de la particule alpha 3B avant la 30 première interaction de la particule alpha 3A, 3B avec le milieu de détection du détecteur 1 de particules alpha 3A, 3B. Dans le mode de réalisation préféré, l'énergie initiale de la particule alpha 3B est déterminée en utilisant l'énergie déposée mesurée et la longueur de la trajectoire mesurée. 14 Avantageusement, l'énergie initiale de la particule alpha 3B est déterminée en utilisant un ratio entre l'énergie déposée par la particule alpha 3B et la longueur de la trajectoire parcourue par ladite particule alpha 3B. L'énergie initiale de la particule alpha 3B est par exemple déterminée en utilisant 5 une relation de type exponentiel. L'énergie initiale E de la particule alpha 3B est typiquement déterminée en utilisant la relation suivante : E = a x (-Ede)b (17) d d Edep est l'énergie déposée par la particule alpha 3B considérée. 10 d est la longueur de la trajectoire parcourue par la particule alpha. a et b sont deux coefficients. Avantageusement, le procédé comprend, avant l'étape El, une étape de calibration E0 de ladite énergie initiale de la particule alpha 3B en utilisant une source de calibration. 15 Dans le cas où l'énergie initiale de la particule alpha 3B est reliée au ratio entre l'énergie déposée par la particule alpha 3B et la longueur de la trajectoire parcourue par ladite particule alpha 3B selon la relation (17), l'étape de calibration E0 permet de déterminer les coefficients a et b. La calibration est typiquement réalisée pour chaque acquisition, c'est-à-dire avant 20 l'analyse de chaque échantillon 5. La source de calibration contient typiquement au moins deux pics énergétiques Ecall et Ecal2 de rayonnement alpha distincts, par exemple entre 4000 keV et 8000 keV. Par « distinct », on entend que les deux pics énergétiques soient au moins séparés de 1000 keV. 25 La source de calibration est par exemple une source de radium 226. Le spectre énergétique de la source de calibration au radium 226 comprend un premier pic énergétique Ecali à 4800 keV (correspondant au radium 226) et un deuxième pic énergétique Ecan à 7700 keV (correspondant au polonium 214, descendant du radium 226 dans la chaîne de désintégration). 30 Pour déterminer les coefficients a et b, il convient alors de résoudre le système d'équations suivant : Ça x xe_ = Ecan = 4800 keV (18) (a x xe = Ecan = 7700 keV 15 x1 et x2 correspondent respectivement aux valeurs de cEcieP associées aux énergies Ecallet Eca12. La figure 5 représente la distribution de EcdeP pour la source de calibration de radium 226 mesurée avec le détecteur 1. 5 La distribution présente deux pics de `eP distincts, correspondant aux pics énergétiques à 4800 keV et à 7700 keV de la source de calibration de radium 226. La formule de Bethe et la courbe de Bragg de dépôt d'énergie des particules alpha dans la matière, connues de l'homme du métier, traduisent le fait que plus l'énergie de la particule alpha augmente, plus le dépôt d'énergie par unité de longueur est faible. Ainsi, 10 chacun des deux pics deEc÷P peut être associé au pic énergétique correspondant. Les valeurs x1 et x2 sont déterminées après un ajustement de deux fonctions gaussiennes sur la distribution. Dans l'exemple de la figure 5, les valeurs x1 et x2 correspondant respectivement Ede aux valeurs de associées aux énergies Ecanet Eca12 sont x1 = 1295 keV/cm et d 15 x2 = 912 keV /cm. Ces valeurs permettent alors de déterminer les coefficients a et b qui valent dans l'exemple considéré a = 7,51.107 et b = -1,34 par résolution du système d'équations mentionné plus haut. La figure 6 présente le spectre en énergie de la source de calibration de radium 20 après application des coefficients de calibration. Les pics énergétiques sont alors centrés sur &ail = 4800 keV et Ecan = 7700 keV. Le procédé comprend alors pour au moins une zone de l'image autoradiographique, une étape de construction E2 d'un spectre en énergie en utilisant les énergies initiales déterminées des particules alpha 3B correspondantes. 25 Le spectre en énergie représente la distribution des énergies initiales de l'ensemble des particules alpha 3B déterminées au cours de l'étape El du procédé. Pour chaque particule alpha 3B détectée, une position moyenne X0,, YO, appartenant au champ de vue du détecteur 1 et une énergie initiale E, associée sont déterminées, selon le procédé décrit précédemment. 30 Typiquement, le procédé comprend, pour une série de K particules alpha 3B détectées durant une acquisition, le stockage sous forme matricielle des positions moyennes X0,, YO, et des énergies initiales E, associées. 16 La zone de l'image autoradiographique pour laquelle le spectre en énergie est construit au cours de l'étape E2 est typiquement une portion du champ de vue du détecteur 1, couramment dénommée zone d'intérêt. La construction du spectre est alors effectuée en utilisant un sous-ensemble des K 5 particules alpha 3B pour lesquelles les positions moyennes X0,, YO, respectives sont comprises dans la zone d'intérêt. La zone d'intérêt est typiquement une surface à deux dimensions d'aire non nulle appartenant à l'anode 13, dans le cas du détecteur 1 gazeux décrit précédemment. La zone d'intérêt a typiquement une forme quelconque, par exemple circulaire ou 10 ellipsoïdale. Avantageusement, l'étape de construction E2 est réalisée pour toute zone appartenant au champ de vue du détecteur 1. En variante, le détecteur 1 est un détecteur à semi-conducteur. Le milieu de détection est désormais un semi-conducteur monolithique ou 15 composé de plusieurs briques agencées en matrices. L'échantillon 5 et le milieu de détection sont typiquement séparés par un espace de vide avantageusement le plus fin possible pour favoriser la localisation spatiale du point d'émission en réduisant l'effet de parallaxe. Par conséquent, les particules alpha 3A, 3B émises par l'échantillon perdent très peu d'énergie avant l'interaction avec le détecteur 20 et la totalité de l'énergie des particules alpha est déposée dans le milieu semi-conducteur. Les signaux sont mesurés sur les pavés 31. La localisation de leur point d'émission est simultanément possible. Toutes les particules alpha 3A, 3B interagissent avec le détecteur 1. Ainsi, aucune opération de sélection des particules alpha n'est requise. 25 Le signal incident mesuré par les pavés 31 est une charge électronique en mouvement. Par polarisation du milieu semi-conducteur, un courant électrique proportionnel à l'énergie déposée par la particule alpha incidente est induit sur les pavés 31. En variante encore le détecteur est un scintillateur sous forme solide, liquide ou 30 gazeux. L'échantillon 5 est idéalement mis au contact du détecteur scintillateur. La particule alpha émise par l'échantillon 5 dépose son énergie dans le milieu scintillant (gaz, liquide, solide). Des photons de scintillations sont alors émis. Les cellules de mesure 29 comprennent typiquement des photo-détecteurs, par 35 exemple des photo-détecteurs SiPM. 17 Le signal incident mesuré par les cellules de mesure 29 comprend au moins un photon. Les photo-détecteurs sont adaptés pour convertir les photons émis par l'interaction entre la particule alpha 3A, 3B et le milieu scintillant en un courant électrique. 5 Ainsi, grâce au procédé selon l'invention, il est possible à la fois de cartographier la surface de l'échantillon 5 et d'identifier à partir des spectres en énergie construits la nature du radioélément émetteur alpha et sa position au sein de l'échantillon. REVENDICATIONS1.- Procédé d'analyse à l'aide d'un détecteur (1) de particules alpha (3A, 3B), d'un échantillon (5) comprenant au moins un radionucléide émetteur d'une pluralité de particules alpha (3A, 3B), le détecteur (1) comprenant un milieu de détection et une pluralité de cellules de mesure (29) adaptées pour mesurer au moins un signal incident généré par une interaction de la particule alpha (3A, 3B) avec ledit milieu de détection, le détecteur (1) étant conçu pour fournir une image autoradiographique dudit échantillon (5), le procédé comprenant les étapes suivantes : - El) dans laquelle on détermine une énergie initiale de chaque particule alpha (3A, 3B), ladite étape de détermination (El) comprenant une répétition des sous-étapes suivantes pour chaque particule alpha (3A, 3B) : - D1) dans laquelle on détermine une position d'une première interaction de la particule alpha (3A, 3B) avec le milieu de détection du détecteur (1) en utilisant au moins une cellule de mesure (29), - D2) dans laquelle on détermine une énergie déposée par la particule alpha (3A, 3B) à l'intérieur du détecteur (1) en utilisant le signal incident mesuré, - D3) dans laquelle on détermine, en utilisant l'énergie déposée déterminée à la sous-étape (D2), l'énergie initiale de la particule alpha (3A, 3B) avant la première interaction de la particule alpha (3A, 3B) avec le milieu de détection du détecteur (1), - E2) dans laquelle pour au moins une zone de l'image autoradiographique, on construit un spectre en énergie en utilisant les énergies initiales déterminées des particules alpha (3A, 3B) correspondantes. 2.- Procédé selon la revendication 1, dans lequel la sous-étape (D1) comprend au moins une opération d'amplification (01, 02) d'une interaction de la particule alpha (3A, 3B) avec le milieu de détection du détecteur (1) à l'aide d'au moins un champ électrique. 3.- Procédé selon la revendication 1 ou 2, dans lequel, avant l'étape (El), on réalise une étape de calibration (E0) de ladite énergie initiale de la particule alpha (3B) en utilisant une source de calibration. 19 4.- Procédé selon l'une quelconque des revendications 1 à 3, dans lequel le détecteur (1) de particules alpha (3A, 3B) est un détecteur gazeux et le signal incident est une charge en mouvement. 5 5.- Procédé selon la revendication 4, dans lequel la position de la première interaction de la particule alpha (3A, 3B) avec le milieu de détection du détecteur (1) est déterminée à partir d'une répartition d'une pluralité de signaux incidents mesurés par les cellules de mesure (29). 10 6.- Procédé selon la revendication 4 ou 5, dans lequel la sous-étape (D2) comprend une opération de mesure (04) de charges (41, 36) générées successivement le long de la trajectoire de la particule alpha (3A, 3B) à l'intérieur du détecteur (1), et une opération de mesure (06) de la longueur de la trajectoire parcourue par la particule alpha (3A, 3B) à l'intérieur du détecteur (1) de particules alpha (3A, 3B), 15 l'énergie initiale de la particule alpha (3A, 3B) étant déterminée en utilisant l'énergie déposée déterminée et la longueur de la trajectoire mesurée. 7.- Procédé selon la revendication 6, comprenant préalablement à l'opération de mesure (06) de la longueur de la trajectoire de la particule alpha (3A, 3B), une opération 20 d'identification (05) des particules alpha (3B) ayant traversé le détecteur (1) de particules alpha (3A, 3B) de part en part, étant entendu que seules les particules alpha (3B) ayant traversé le détecteur (1) de part en part sont considérées à l'étape (E2) de construction des spectres en énergie. 25 8.- Procédé selon la revendication 7, dans lequel l'opération d'identification (05) comprend un calcul d'un paramètre géométrique d'une distribution des charges mesurées dans le détecteur (1) de particules alpha (3A, 3B) selon une direction du détecteur (1) de particules alpha (3A, 3B). 30 9.- Procédé selon la revendication 8, dans lequel le paramètre géométrique est un coefficient d'aplatissement de la distribution des charges mesurées. 10.- Procédé selon la revendication 9, dans lequel l'opération d'identification (05) comprend l'application d'un seuil sur le coefficient d'aplatissement de la distribution des 20 charges mesurées et la sélection des particules alpha ayant un coefficient d'aplatissement inférieur audit seuil. 11.- Procédé selon l'une quelconque des revendications 6 à 10, dans lequel 5 l'énergie initiale de la particule alpha (3B) est déterminée en utilisant un ratio entre l'énergie déposée par la particule alpha (3B) et la longueur de la trajectoire parcourue par ladite particule alpha (3B). 12.- Procédé selon la revendication 11, dans lequel l'énergie initiale de la particule 10 alpha (3A, 3B) est déterminée en reliant ladite énergie initiale au ratio entre l'énergie déposée par la particule alpha (3B) et la longueur de la trajectoire parcourue par ladite particule alpha (3B) par une relation de type exponentiel. 13.- Procédé selon la revendication 12, dans lequel la relation de type exponentiel 15 répond à la formule : E (cEd) = a x (cEd)b - Edep étant l'énergie déposée par la particule alpha (3B) considérée, - d étant la longueur de la trajectoire parcourue par la particule alpha (3B), - a et b étant des coefficients. 20 14.- Procédé selon l'une quelconques des revendications 1 à 3, dans lequel le détecteur (1) est un détecteur à semi-conducteur ou un détecteur scintillateur et le signal incident mesuré est respectivement une charge en mouvement ou un photon. 25