Un survol d'hélicoptère de six secondes identifie des sources radioactives mixtes

Un hélicoptère qui traverse un site contaminé à environ 10 m/s ne passe que quelques secondes à proximité d'une source donnée. Une équipe de l'Institut Paul Scherrer (PSI), en Suisse, et de l'Université du Michigan affirme que cela suffit. Dans un article publié le 26 octobre 2026 dans Nature Communications, les chercheurs exploitent 6 s de données issues d'un seul passage, validées par rapport à des mesures indépendantes in situ et en laboratoire. Ils rapportent une identification décisive de mélanges de sources (>5σ), une localisation à l'échelle du mètre et une restitution des activités au pourcent près. Le résumé de l'université chiffre ces résultats plus concrètement : sources localisées à moins de deux mètres et activités estimées à moins de deux points de pourcentage de leurs valeurs réelles.

L'avancée est avant tout computationnelle. Le détecteur est un réseau commercial de NaI(Tl). Ce qui a changé, c'est la façon dont l'équipe construit, et évalue rapidement, le modèle physique auquel sont comparés les comptages épars.

01Pourquoi un spectre de survol est difficile à lire

Chaque radionucléide émet des rayons gamma à des énergies caractéristiques. En principe, le spectre mesuré permet d'identifier l'isotope et sa quantité. En pratique, chaque matière possède une empreinte constituée de son spectre d'énergie gamma, et cette empreinte se dégrade en vol. Les photons sont diffusés dans l'air et dans la structure de l'aéronef, et perdent de l'énergie avant d'atteindre les cristaux. Les comptages sont en outre rares. Selon le PSI, sur environ 10 000 photons émis par une source au sol, un seul atteint l'hélicoptère, et celui-là est encore diffusé par les surfaces métalliques à l'intérieur de l'appareil.

Les auteurs qualifient d'ill-posed, c'est-à-dire mal posée, l'inversion qui en résulte. La spectrométrie gamma mobile est la référence pour l'identification des sources, mais les approches actuelles manquent de précision, car les données sont rares et l'inversion est mal posée. Le remède classique est l'analyse en spectre complet. Au lieu d'intégrer quelques fenêtres d'énergie, elle ajuste l'ensemble du spectre mesuré à des gabarits simulés. Le problème, c'est le coût. Des gabarits fidèles exigent un transport Monte-Carlo à travers l'air et la cellule de l'appareil, ce qui peut prendre des heures, voire des jours, pour chaque géométrie de vol, chaque atmosphère et chaque position de source. Personne ne peut effectuer ce calcul en pleine situation d'urgence.

02Scinder le problème du transport en deux

L'astuce centrale est la factorisation. L'équipe sépare le trajet d'un photon en deux volets : le transport atmosphérique, et la réponse de l'hélicoptère et du détecteur. D'après le résumé du Michigan, les calculs atmosphériques tiennent compte de la température et de l'humidité de l'air, tandis que le volet aéronef repose sur des fonctions précalculées propres au modèle d'hélicoptère. La section méthodes de l'article décrit une « banque de flux gamma » générée avec le module de désintégration radioactive du code FLUKA. Elle décrit aussi des simulations distinctes de la plateforme, qui enregistrent le dépôt d'énergie dans chaque cristal et le réduisent à des spectres d'amplitude d'impulsions attendus.

Le gros du travail Monte-Carlo est donc réalisé une seule fois, en amont. Le preprint fait état d'environ 10² heures-cœur par source, avec 2 × 10⁹ particules primaires réparties sur 10² exécutions indépendantes. Une géométrie de vol donnée ne nécessite ensuite qu'une consultation et une étape de combinaison. Selon Breitenmoser, c'est cette réorganisation, plutôt qu'une puissance de calcul accrue, qui a ramené à des microsecondes des simulations qui prenaient des heures ou des jours.

Le modèle de la plateforme s'appuie sur des travaux antérieurs du PSI. Le système suisse comprend quatre cristaux prismatiques de NaI(Tl) d'un volume total de 16,8 L, installés dans la soute d'un Super Puma. Dans une étude de validation antérieure, l'équipe avait rapporté une erreur spectrale relative médiane inférieure à 10 % pour la plupart des radionucléides dans le modèle de réponse du détecteur.

03Inférence bayésienne sur des milliers de gabarits

Avec des gabarits aussi peu coûteux, l'inférence peut se permettre d'être exhaustive. Le cadre méthodologique compare des milliers de spectres prédits au spectre mesuré et met à jour en continu sa confiance statistique quant à ce qui se trouve au sol. Les mots-clés du preprint désignent la méthode de Monte-Carlo par chaînes de Markov comme échantillonneur. Les inconnues sont estimées conjointement : les nucléides présents, la position de la source et son activité. Le bruit de fond naturel dû au potassium, au thorium et à l'uranium du sol est lui aussi ajusté.

Le résultat est une distribution a posteriori plutôt qu'un unique meilleur ajustement. Cela compte sur le plan opérationnel. Un intervenant obtient une estimation d'activité assortie d'un intervalle de crédibilité, ainsi qu'une préférence quantifiée pour « Cs-137 plus Ba-133 » par rapport à « Cs-137 seul ». Le résultat à >5σ sur le mélange est un énoncé de ce type.

04Ce qui a réellement été testé

Les données de terrain proviennent de l'exercice ARM22, sur le terrain d'entraînement militaire de Thoune, le 16 juin 2022. Selon le preprint, les deux configurations utilisaient les deux mêmes sources scellées, du Cs-137 et du Ba-133, avec une acquisition à 1 Hz. Les configurations étaient les suivantes :

• Vol stationnaire : environ 5 min à environ 90 m, utilisé pour un étalonnage contrôlé des activités et du bruit de fond.
• Passage unique : environ 30 m d'altitude et environ 10 m/s, sur trois scénarios : Cs-137 seul, Ba-133 seul, et une double source colocalisée avec moins de 20 cm d'écart.

Les valeurs de référence étaient indépendantes. Les activités des sources provenaient des mesures de laboratoire du fabricant, et le bruit de fond naturel était contraint par une spectroscopie in situ à 23 positions sur le terrain d'entraînement. À 10 m/s, six secondes correspondent à une trajectoire d'environ 60 m (mon calcul). L'inférence repose donc sur un court tronçon de vol autour de la source.

À titre de comparaison, les travaux bayésiens antérieurs sur les levés par véhicule routier ne sont pas directement comparables. Une étude réalisée à 50 km/h avec un échantillonnage à 1 s a constaté que les écarts prédits sur l'activité et la position restaient autour de 30 %, en raison de l'incertitude des données. Cette étude employait des détecteurs et une géométrie différents. Elle illustre néanmoins l'écart entre un modèle approximatif de source ponctuelle et un modèle en spectre complet propre à la plateforme.

05Pourquoi c'est important, et ce qui reste à démontrer

En cas d'accident de réacteur ou de détonation nucléaire, les aéronefs peuvent être le seul moyen d'atteindre la zone touchée. L'équipe soutient que les données acquises en vol peuvent éclairer les décisions d'évacuation, de mise à l'abri et de distribution d'iodure de potassium. Pour la sécurité nucléaire, la propriété utile est la résolution des mélanges. Distinguer une source médicale de Cs-137 d'un inventaire mixte ou industriel, et le faire en un seul passage, change ce que fait ensuite une équipe de recherche. Les auteurs citent également la surveillance courante d'installations, la cartographie géochimique et même l'abondance de radionucléides extraterrestres parmi les applications. Les documents de presse y ajoutent la remédiation environnementale et la prospection minière.

Les limites apparaissent dans la conception même de l'étude :

• Peu de nucléides : les essais ont porté sur deux sources ponctuelles scellées. Un rejet réel implique des dizaines d'isotopes, des panaches étendus et des dépôts au sol ; la bibliothèque de gabarits et le problème d'inférence s'en trouvent tous deux élargis.
• Une seule plateforme et un seul site : les fonctions de l'aéronef sont propres à un modèle d'hélicoptère. Une autre cellule nécessite son propre modèle Monte-Carlo validé et son propre précalcul.
• Géométrie favorable : les trajectoires de vol passaient directement au-dessus des sources, dont la présence était connue. Une recherche à l'aveugle sur une zone inconnue constituerait un test plus difficile.
• Pourcentages : les chiffres de « deux points de pourcentage » et de « deux mètres » proviennent des communiqués de presse. L'article lui-même parle de « percent-level » et de « meter-scale », c'est-à-dire d'un niveau du pourcent et d'une échelle du mètre. Les lecteurs devraient consulter les tableaux par scénario de l'article avant de citer un chiffre unique.

Les prochaines étapes ne sont pas détaillées dans les sources que j'ai examinées. Il serait raisonnable d'envisager davantage de nucléides, des sources étendues et distribuées, ainsi que des essais en conditions opérationnelles. Les données sont archivées sur Zenodo, de sorte que des groupes indépendants peuvent tester l'approche sur les mêmes spectres. Le résultat principal est qu'une factorisation fondée sur la physique a rendu un traitement bayésien rigoureux assez rapide pour être appliqué aux données d'un vol de quelques secondes.