Brevet FR2641615A1 — publié le 19900713

SYSTEME POUR DETECTER LA PRESENCE D'UN SIGNAL PUR DANS UN SIGNAL BRUITE DISCRET MESURE A TAUX MOYEN DE BRUIT CONSTANT DESCRIPTION La présente invention a pour objet un système pour détecter La présence d'un signal pur dans un signal bruite discret mesure, cette détection étant réalisée avec une probabiLité de fausse détection inférieure à un taux de fausse détection prédéteriiné. L'invention concerne La détection des signaux discrets, c'est-à-dire les signaux produits par des evéncments discontinus. L'invention s'applique notamment à La détection des rayonnements corpuscuLaires tels que Les rayons neutroniques ou Les rayons gamma. La détection d'un signal pur dans un signaL bruite mesuré suppose que l'on peut s'affranchir de la présence du bruit dans ledit signal mesuré. L'invention concerne un procédé de détection lorsque Le taux moyen de bruit varie dans un intervalle ayant une borne supérieure MB. Cette condition est pratiquement toujours réalisée pour les rayonnements de neutrons car il y a très peu de neutrons Libres dans l'environnement. Elle est également souvent réaLisée en un lieu donné pendant une certaine durée de temps pour les rayonnements gamma. La détection d'un rayonnement corpusculaire particulier, par exemple nucléaire, se fait en présence d'un bruit corpusculaire ayant plusieurs causes. Dans le cas d'un rayonnement nucléaire, une partie du bruit est due au rayonnement cosmique et l'autre partie est due à la radioactivité naturelle des terrains et des roches. Il existe aussi un bruit dû au détecteur utilise. Le rayonnement cosmique peut varier au cours du temps et le bruit corpusculaire dû aux roches et au terrain est essentiellement variable, selon le Lieu où se fait la mesure. Un rayonnement particulier ne peut donc être détecte que si L'on connatt le bruit ou bien si l'on en a une estimation suffisante. Le stockage des informations est effectué pendant une durée de temps limitée dt. Ce stockage des informations se fait de façon discrète : Lorsqu'une particule est détectée, Le contenu d'une mémoire est augmenté d'une unité. En général, un ou plusieurs parametres définissent l'événement nucléaire détecte. Dans le cas où cet événement ne dépend que d'un paramètre unique, celui-ci peut être par exempLe l'énergie : c'est le cas typique de la détection d'un rayonnement gamma à l'aide d'un germanium intrinsèque ou d'un scintillateur NaI associé à un photomultipLicateur. Ce paramètre unique peut aussi etre le temps, par exemple lorsque l'on mesure le teops de vol d'une particule. Usuellement, pour savoir si un signal a été détecté à l'aide d'un stockagé effectué durant un intervalle de temps dt, la pratique courante consiste à apprécier de quel Le valeur Le comptage obtehu dépasse l'estimation du bruit dont on dispose. Ce bruit peut etre estimé de différentes manières. IL peut etre connu à L'avance par une mesure faite au préalable dans des circonstances où L'on est sûr qu'il n'y a pas de signal particulier. Il peut être également apprécié en cours d'expérience en L'absence de toute source de signal en extrapolant ou en interpoLant les niveaux du signal détecté dans différentes intervaLLes, par exemple énergétiques. Pour ceLa, il faut que le taux de comptage dans les intervalles où l'on est sOr qu'il n'y a pas de bruit soit suffisant, pour pouvoir faire une estimation de la densité de probabilité du bruit dans les intervalles où L'on cherche à détecter un rayonnement particulier. De manière classique, on estime qu'il y a effectivement un signal lorsque le taux de comptage N est tel que N > B+K. 'CM, où M est le taux moyen de bruit estime et K une constante égale à quelques unités. Ce critère connu de détection n' est pas satisfaisant car il est empirique, Là valeur de K étant choisie de façon non rationnelle en se fondant sur des expériences antérieures. Un but de l'invention est de fournir un système de détection fonte sur un critère rationnel, déduit des lois de probabilité d'apparition du bruit et du rayonnement attendu. La détection selon l'invention pressente ainsi l'avantage de ne dépendre que d'un critère objectif. Dans l'art antérieur, un taux de comptage n'est pris en compte que s'il est supérieur au taux moyen de bruit augmenté d'un certain nombre de fois l'écart typeM. On connait par ailleurs une méthode décrite dans La demande de brevet faite par le déposant et publiée sous le numéro Z 580 819 qui repose sur un principe différent. En effet, dans cette demande, on considère, au contraire, que tout dépassement du taux moyen de bruit constitue une information sur La présence ou L'absence d'un signal pur. Cette information est faible lorsque la différence entre le taux de comptage et le taux moyen de bruit est faible. Plus la différence entre le taux de comptage et le taux moyen de bruit est importante, plus L'information apportée par la mesure est importante et plus la probabilité de fausse détection est faible. Le critere pratique de détection d'un signaL pur est donc l'obtention d'une information suffisante, c'est-à-dire supérieure à un seuil Im, c'est-à-dire une probabilité de fausse détection inférieure à un taux de fausse détection prédéterminé. La détection selon cette méthode consiste donc à comparer la probabilité de fausse détection PF, définie par le taux de comptage N mesuré et une estimation du taux moyen de bruit constant M, à un taux de fausse détection prédéterminé TF. Un signal pur est détecté si ladite probabilité est inférieure audit taux. Dans la mesure où L'on connait la probabilité d'occurrence du nombre d'événements détectés pendant chaque intervalle de temps, La probabilité de fausse détection est exprimée en fonction du taux de comptage N et du taux moyen de bruit estimé M. Cette densité de probabilité est par exempLe celLe de la loi de Poisson pour les rayonnements corpusculaires. L'utilisation directe de la relation entre la probabilité de fausse détection PF, Le taux de comptage N et le taux moyen de bruit M pour décider Si un rayonnement particulier est détecte, conduit à -des calculs longs, car la probabilité de fausse détection PF est sous forme d'une série dépendant du taux de comptage N et du taux moyen de bruit M. La durée du calcul dépend de La valeur du taux de comptage N et est usuellement incompatibLe avec une détection en temps réel, si la durée de chaque comptage est de L'ordre de la seconde par exemple. Pour remédier à cet inconvénient, et permettre une détection en temps réel, il a d'abord été proposé d'utiliser une table contenant pour une pluraliste de valeurs de taux de comptage La valeur du taux moyen de bruit associée, telle que la probabilité de fausse détection PF pour ce taux moyen de bruit et ce taux de comptage soit comparé au taux de fausse détection prédéterminée TF. Cette méthode a été appliquée à la conception d'appareilLages portatifs et a consiste à faire une estimation de l';nformation à partir d'une valeur de -3 seuil correspondant à un taux de fausse alarme de 10 Il s'est avéré que, même pour des vitesses de déplacement de L'ordre de 1 m/s, le taux des alarmes entait excessif, par suite des variations importantes du bruit de fond naturel ; il a alors été décide d'utiliser un algorithme plus adapté qui a été mis en oeuvre dans un appareillage embarquable. Cet algorithme est basé sur le principe d'une "Référence Fluctuante" et on peut en trouver une description dans la demande de brevet effectuée par le déposant et publiée sous le n0 2 579 767. Cet algorithme implique le découpage du spectre en énergie, en quelques fenbtres (en pratique quatre fenêtres ont été choisies) ; dans chacune de ces fenêtres, ta probabiLité de faux signal, selon l'algorithme dit de la "Référence Fluctuante", est évaluée en ligne. Or, avec cet algorithme il existe une limite quant à la possibilité de couvrir tout Le spectre, par suite de La résolution limitée du détecteur, en particulier si ce détecteur est du type NaI. Outre l'ensemble de tous ces inconvénients auxquels la demanderesse s'est heurté, La demanderesse a cherche à faciliter la réalisation des détecteurs portatifs. Pour cela, la demanderesse a cherché à simplifier l'électronique utilisée dans Le détecteur portatif connu, par exemple cet Le qui est décrite dans la demande de brevet déjà évoquée et publiée sous Le no 2 580 819 en supprimant l'utilisation de la table de valeurs, ce qui l'a conduite à supprimer la mémoire dans laquelle cette table a été enregistrée. La suppression de cette table n'a pu avoir lieu que parce que la demanderesse s'est posé le problème de La détection de façon différente. En effet, cette simplification est le résultant d'une mise en forme semi-phénomenologique de la quantité d'informations I obtenues en fonction du comptage. La quantité d'informations I est exprimée selon l'invention par une fonction F liee à La borne superieure du bruit MB et au taux moyen de signal observé pendant une mesure de durée déterminée. Comme cela apparait dans la fonction F, au lieu de choisir dans L'évaluation de l'information le taux moyen de bruit estimé M, La demanderesse a choisi La borne supérieure MB du bruit M. Cette conception nouvelle obtenue par une approche différente de celles déjà connues conduit à la réalisation d'un détecteur plus petit, moins cher et donc à un emploi banalisé. La présente invention a plus précisément pour objet un système pour détecter la présence d'un signal pur dans un signal bruité discret mesuré à taux moyen de bruit constant, avec une probabilité de fausse détection PF inférieure à un taux de fausse détection TF prédéterminé, comprenant en série un détecteur pour convertir le rayonnement reçu en impulsions électriques, un moyen de mise en forme, une échelle de comptage binaire, un moyen de traitement et un moyen de signalisation, le moyen de traitement comprenant un processeur et un circuit d'entrée-sortie, ledit système étant caractérisé en ce qu'il permet de : A) - au préalable, effectuer une phase d'initialisation consistant à :: - mesurer un comptage N0 pendant n périodes de temps #t, successives égales, - estimer une valeur du taux moyen de bruit M à partir de ce comptage N0, - déterminer une borne supérieure MB du bruit telle que ne= /3M, ss étant un nombre entier supérieur à 1, B) - effectuer une phase de détection de signal consistant à - mesurer un comptage N durant au moins une période de temps #t1, - estimer un comptage M1 attribué audit signal tel que M1=N i -MB, - calculer la quantité d'informations telle que :: Ml Tri MB112 avec MB régal à la valeur estimée, - émettre un signal S lorsque le taux de signal Ml est positif, la fréquence de ce signal étant une fonction croissante de la quantité d'informations I. Conformément à l'invention, La fonction F permettant de calculer la quantité d'informations dépend de deux fonctions a et K de La borne supérieure MB du bruit estimé de sorte que : &alpha;(MB) M1 I = K(MB) -----1-2---- (2) 112 MB Dans une variante de réalisation, La phase de détection consiste à : - mesurer un comptage pendant plusieurs périodes At successives en effectuant une division de chaque période #ti en k périodes successives Dtii=1,k, et i i en mesurant un comptage n durant chaque période Dtj, la mesure du comptage N se faisant de sorte que :: Ni = Ni-1 + ni - ni-k dans lequel Ni-1 est le comptage mesuré durant la i-I période prEcedente sst; n entant le comptage -effectué durant la période Dt et n k étant le comptage effectué durant la période Dt i-k Avantageusement, le signal S émis à la fin de chaque comptage obtenu à la fin de chaque durée At est un signal sonore de la forme S=H(I).fo, fo étant la fréquence du signal sonore émis dès que le taux de signal M1 est positif et H(I) une fonction croissante de la quantité d'information I. Avantageusement, dans le cas où l'on procède à un échantillonnage de la durée de comptage, le signal S semis est un signal sonore dont la fréquence varie à l'issue de chaque comptage Dtj, ce signal étant de la forme S=H(I).fo, fo étant la fréquence du signal émis des que le taux de signal Ml est positif et H(i) une fonction croissante de la quantité d'information I. D'autres caractéristiques et avantages ressortiront mieux de La description qui suit et qui est donnée à titre illustratif et non limitatif en référence aux figures annexées, dans Lesquelles - La figure 1 représente le schéma de réalisation d'un système selon l'invention, - la figure 2 représente les différentes retapes de la détection opérée par Le système selon l'invention, - la figure 3 représente une variante de réalisation des étapes de détection dans le cas d'un échantillonnage de la durée de comptage At, - la figure 4 représente les variations des fonctions o et K qui sont des fonctions de la borne supérieure du bruit MB, - la figure 5 représente la quantité d'informations I en fonction du taux moyen de signal Ml, à taux moyen de bruit MB constant. L'invention s'applique à tous les rayonnements discrets, tels que les rayonnements corpusculaires, et concerne plus précisément l'analyse du signal délivre par un détecteur pour décider, avec une probabilité de fausse détection inférieure à un taux de fausse détection prédéterminé, si le signal analyse contient un signal pur dO à la présence d'une source particulière ou si ce signal se compose uniquement de bruit. A titre d'exenple, le procédé de l'invention s'applique notamment à la détection des rayons gamma et des rayons neutroniques. L'appareil représente sur La figure 1 comprend un détecteur 2, un moyen de mise en forme 4, une échelle de comptage binaire 6 pouvant comprendre notamment un convertisseur analogique-numérique, un moyen de traitement 8 et une interface de dialogue homme-machine 10 classique en soi constituée par des boutons de commande (par exemple des interrupteurs), un convertisseur tension-fréquence et un moyen de signalisation (une alarme par exemple). Le détecteur 2 dépend du rayonnement à détecter. Pour la détection de rayons gamma, ce détecteur peut Entre un germanium intrinsèque ou un scintillateur NaI associé à un photomultiplicateur ; pour la détection de neutrons, le détecteur est par exemple un compteur à He3. Le détecteur 2 délivre une impulsion de tension pour chaque événement détecté. Cette impulsion est reçue par le moyen de mise en forme 4. Ce dernier peut comprendre notamment en série un amplificateur, un circuit de mise en forme et un discriminateur. L'échelle de comptage binaire 6 est disposée en aval du moyen de mise en forme 4. Elle compte le nombre d1impulsions délivrées par le moyen 4 pendant des durées de comptage t Le moyen de traitement 8 comprend un processeur 14 et un circuit d'entrée-sortie 18. Ces éléments sont reliés par une voie 20 comprenant des lignes d'adresses, des lignes de données et des lignes de commande. Le moyen de signalisation et de commande 10 est en liaison avec le processeur 14 à travers le circuit d'entree-sortie 18 et une voie 21. Il comprend notamment des moyens d'alerte visuels et sonores qui sont déclenchés lorsqu'un signal pur est détecté et des moyens d'alimentation et des interrupteurs (non représentes) pour envoyer des signaux de commande correspondant à des commandes extérieures C et P. On va maintenant décrire le fonctionnement du système de détection de l'invention en considérant, à titre d'exemple,- le cas de La détection d'un rayonnement corpusculaire, tel qu'un rayonnement gamma ou un rayonnement neutronique. Le nombre de particules détectées pendant une durée de mesure sst; est une variable aléatoire dont la densité de probabilité est celle de la loi de Poisson.La probabilité P de détecter durant Le temps pst;, I événements dus aux rayonnements particuliers et J événements dus aux bruits est de I J -(Mî+M) Ml M P e . ~~~ . ~~ I! J! où Ml et M sont les taux moyens associés aux variables aléatoires représentant le nombre de particules détectées dues au rayonnement particulier et le nombre de particules détectées dues au bruit. La probabilité pour qu'il existe un rayonnement particulier, que le stockage N contienne ou non un événement qui lui soit dO, est alors égale à : la probabilité complémentaire est la probabilité de fausse détection PF qui définit le probabilité pour qu'un rayonnement particulier soit détecte alors qu'aucun rayonnement particulier n'existe. Le systeme de détection de L'invention permet de traduire cette probabilité sous forme d'une information I qui est fonction du taux moyen de signal Ml et de La borne supérieure du bruit estimé MB. Les différentes opérations sont décrites en référence à la figure 2 qui représente les différentes opérations dans La détection effectuée par Le système. A) - Au préalable, l'on effectue une phase d'initialisation consistant à : - mesurer un comptage N durant n périodes de temps at; O i successives de durée égale, - puis estimer une valeur du taux moyen de bruit M, à partir de ce comptage No, - puis déterminer une borne supérieure MB du bruit M telle que MB= ssM, ss étant un nombre entier supérieur à I que l'on se fixe en fonction de ta sensibilité desirée, et ensuite, B) - effectuer une phase de détection du signal consistant à :: - mesurer un comptage N durant une ou plusieurs périodes de temps At , - puis estimer un comptage M1 attribué audit signal tel que : M1=Nj-MB, - puis calculer la quantite d'information I tel Le que : Ml I = F(-----) (1) 112 MB avec MB égal à La valeur estimée, et enfin, - émettre un signal S lorsque Le taux de signal Ml est positif, la fréquence de ce signal étant une fonction croissante de la quantité d'information. L'estimation du taux moyen M du bruit sur tout le spectre du signal reçu se fait de façon classique tel que décrit par exemple dans la demande de brevet français publiee sous le n 2 580 819. La fonction F est préalablement enregistrée en mémoire du processeur, qui exécute le calcul avec les données Ml et MB qui lui sont entrées préalablement. Le bloc I de la figure 2 représente une étape consistant à rentrer en mémoire une fois pour toutes la fonction F, la durée des périodes at;. Le bloc A représente la phase d'initiatisation qui a lieu avant toute détection et qui consiste à effectuer une estimation du taux moyen de bruit M et de retenir sa borne supérieure. Ce bloc A comporte les blocs Il à V correspondant à chacune des retapes de la phase d'initialisation A. Ce bloc B représente la phase de détection du signal consistant à déclarer qu'il y a du signal des que le taux M1 est positif et à émettre un signal S dont la fréquence est une fonction croissante de La quantité d'information, cette quantité I ayant une expression semi-phénoménologique donnée sur La figure 4. Le bloc B comporte les blocs VI à XII correspondant aux différentes retapes ayant lieu lors de la phase B. la suite de cette phase B, l'utilisateur du détecteur peut recommencer autant de fois qu'il le désire les comptages durant des durées #t. Selon une variante de réalisation et afin d'augmenter la portée du détecteur, la phase B comporte quelques modifications dues au fait que chaque période de temps st; est divisée en k périodes de temps Dt i i de comptage, i=1 à k. Cette variante correspond aux déroulements des étapes représentées sur la figure 3. Le nombre k est donc préalablement stocké en mémoire du processeur. D'une façon générale, l'on procède à plusieurs comptages successifs Ni-i, Ni, N1+1, Ni+2' jusqu'à ce que l'utilisateur décide d'arrêter. L'on effectue à l'intérieur de chaque période de temps de durée Dt un comptage n et au bout d'une i i durée #ti, la mesure du comptage N est obtenue en procédant au calcul suivant Ni = Ni-1 + ni - ni-k- Le comptage N correspond au comptage i-I obtenu à la période précédente, n correspond au comptage durant la periode Dt de la période #t1+1 et, ni-k correspond au comptage durant la période Dt1-k. Les étapes correspondant à cette variante sont représentées sur la figure 3. Selon cette réalisation, les blocs représentés sur cette figure s'insèrent entre les blocs V et VIII de la figure 2 à la place des blocs VI-VII. Le bloc XV représente une étape de comptages successifs pendant une durée k.Dt . En effet, selon cette variante, l'on effectue un premier comptage pendant une durée k.Dt Le paramètre k correspondant au nombre de sous-périodes Dt est entré en mémoire. Le bloc XVI représente l'étape de mesure de comptage pendant chaque intervalle de temps Dt. Les comptages N1-1, n1, n1, ni+k, i=1, k sont enregistrés en mémoire afin de calculer pour chaque periode de durée #t1, le comptage N tel que : Ni = Ni-1 + - ni-k. Pour les deux variantes de réalisation un signal S pouvant etre un signal sonore est émis. Ce signal est une fonction croissante H de La quantité d'information I. Pour ceLa, tant que Le taux de signal M1 est inférieur à zéro (bloc X, fig. 2), un signal numérique de commande D correspondant par exemple à un niveau bas (état O) est envoye par Le processeur -aux moyens de signalisation (bloc XI, fig. 2) qui ne vont alors émettre aucun signal. Un nouveau comptage ou une nouvelle série de comptages peut être recommencée. Dans le cas contraire, des que le taux de signal Ml est positif, le signal numérique de commande D a alors un niveau haut (état 1) et est envoyé par le processeur aux moyens de signalisation (bloc XII, fig. 2). Le signal S est par exemple un signal sonore dont la fréquence est une fonction croissante de la quantité I. Le signaL S est par exemple obtenu au moyen d'un convertisseur tension-fréquence qui reçoit la quantité d'information I : S=fo.H(I). Lorsque la quantité obtenue dépasse le seuil zéro, le signal émis S a une fréquence fonction de la fréquence fo du signal remise pour le seuil. Sur la figure 4, on a représenté les variations des fonctions a et K en fonction de la borne supérieure du bruit MB. Ces variations permettent pour une borne MB donnée de connaStre directement par simple lecture sur ces courbes ta valeur de la fonction a et de la fonction K à utiliser pour le calcul de la quantité I qui est mis sous la forme Ml a I = K a 2) (2) 112 MB M1 est connu, MB est connu, a et K sont obtenus par ces courbes ou bien calculés à partir de la courbe de variation de I en fonction de Ml pour la borne MB désirée. En effet, on peut également, au lieu de lire directement les valeurs de a et de K sur Les courbes de la figure 4, calculer ces valeurs à partir de La courbe I que l'on a alors enregistrée en fonction du taux de signal Ml pour la borne MB désirée. Les coordonnées de deux points quelconques P1 et P2 pris sur la courbe enregistrée sont utilisées pour établir un système de deux équations à deux inconnues à résoudre. En effet, les coordonnées Mî(Pî), I(P1)] [M1(P2), I(P2)) de ces points sont reportées dans l'expression I donnee par La relation (2). Les relations ainsi obtenues pour chaque point PI, P2 forment un système d'équation à deux inconnues a et K que l'on résout classiquement pour obtenir a et K. Par ailleurs, on prevoit dans le dispositif un bouton de commande C placé sur L'interface de dialogue homme-machine 10 pour que ces derniers envoient un signal (un niveau de tension donné par exemple) vers une entre de commande du microprocesseur 14 afin de déclencher la phase B de détection. Un autre bouton de commande P est prévu sur l'interface de dialogue homme-machine 10 pour que ces derniers envoient un signal (un niveau de tension donné par exemple) vers une entrée de commande du microprocesseur afin de sélectionner le mode de fonctionnement correspondant à l'une ou à l'autre des variantes décrites. REVENDICATIONS 1. Système pour détecter la présence d'un signal pur dans un signal bruite discret mesuré à taux moyen de bruit constant, avec une probabilité de fausse détection PF inférieure à un taux de fausse détection TF prédéterminé, comprenant en série un détecteur (2) pour convertir le rayonnement reçu en impulsions électriques, un moyen de mise en forme (4), une échelle de comptage binaire (6), un moyen de traitement (8) et un moyen de signalisation (10), le moyen de traitement comprenant un processeur (14) et un circuit d'entrée- sortie (18), ledit système etant caractérisé en ce qu'il permet de A) - au préalable, effectuer une phase d'initialisation consistant à - mesurer un comptage N pendant n périodes de temps #t1, successives égales, - estimer une valeur du taux moyen de bruit M à partir de ce comptage No, - déterminer une borne supérieure MB du bruit telle que MB= ssM, ss étant un nombre entier supérieur à 1, B) - pour détecter un signal : - mesurer un comptage N durant au moins une période de temps At1, - estimer un comptage Ml attribué audit signal tel que M1 =N; -NM, - calculer La quantité d'informations I telle que M1 I = F (-------) (1) MB 1/2 avec MB égal à ta valeur estimée, - émettre un signal S lorsque te taux de signal M1 est positif, la fréquence de ce signal étant une fonction croissante de la quantité d'informations I. 2. Système pour détecter la présence d'un signal pur dans un signal bruite selon la revendication 1, caractérisé en ce que la fonction F dépend de deux fonctions a et K de la borne supérieure MB du bruit estimé de sorte que M1 a(MB) I = K(MB) ------------ (2) 1/2 MB 3.Système pour détecter la présence d'un signal pur dans un signal bruité selon la revendication 1 ou 2, caractérisé en ce que La phase de détection B consiste à : - mesurer un comptage pendant plusieurs périodes t successives en effectuant une division de chaque période #t en k périodes successives Dt i=1,k, et en mesurant un comptage ni durant chaque période Dt, la mesure du comptage N se faisant de sorte que : Ni = Ni-1 + ni - ni-k dans lequel Ni-1 est Le comptage mesuré durant La i-1 periode précédente #t1-1' ni étant le comptage effectué durant la période Dti et ni-k étant le comptage effectué durant la période Dt k. 4. Système pour détecter la présence d'un signal pur dans un signal bruite selon l'une quelconque des revendications 1 à 3, caractérisé en ce que Le signal S émis à la fin de chaque comptage obtenu à la fin de chaque durée #t est un signal sonore de La forme S=H(I).fo, fo étant la fréquence du signal sonore émis dès que le taux de signal M1 est positif, H(I) entant une fonction croissante de la quantité I. 5. Système selon la revendication 3, caractérisé en ce que le signal S émis est un signal sonore dont La fréquence varie à l'issue de chaque comptage de durée Dt1, ce signal étant de la forme S=H(I).fo, fo etant la frequence du signal émis dès que le taux de signal M1 est positif et H(I) une fonction croissante de la quantité d'information I.