Brevet FR2485752A1 — publié le 19811231

PROCEDE ET DISPOSITIF DE MESURE-DE RAYONS GAMMA DANS UN SONDAGE L'invention se rapporte aux techniques de diagraphie nucléaire dans les sondages et notamment à un procédé et dispositif stabilisé pour la détection de rayons gamma provenant de formations traversées par un sondage. On connait plusieurs techniques de diagraphie nucléaire dans lesquelles on détecte les rayons gamma provenant des formations. Ces rayons gamma peuvent être dûs à la radioactivité naturelle des formations ou résulter de l'irradiation de ces formations par une source de neutrons ou de rayons gamma. La détection de ces rayons gamma, qu'ils proviennent d'une radioactivité naturelle ou provoquée, peut s'effectuer séparément pour plusieurs plages d'énergie ou fens- tres distinctes afin d'obtenir des informations sur le spectre d'énergie de ces rayons gamma. Le spectre obtenu fournit ensuite des indications sur la lithologie des formations. A titre d'exemple, un dispositif d'analyse spectrale des rayons gamma naturels est décrit dans le brevet des Etats-Unis d'Amérique No. 5 976 878 (P. Chevalier et B. Seeman). Dans ce dispositif, le système de détection comprend un cristal à scintillation, un photomultiplicateur, un amplificateur linéaire et un analyseur de hauteur d'impulsions. A la sortie de l'amplificateur apparaissent des impulsions électriques dont l'ampli- tude est proportionnelle à l'énergie des rayons gamma reçus par le cristal. Dans l'analyseur de hauteur d'impulsions le spectre d'énergie des rayons gamma est divisé en cinq fenêtres au moyen de comparateurs ayant des seuils prédéterminés A1 à A6. Les taux de comptage des différentes fenêtres permettent de calculer les teneurs en thorium, uranium et potassium (T, U, K) des formations. Dans cette technique, comme dans toutes celles qui utilisent une étude spectrale, c'est-à- dire une discrimination par hauteur d'énergie du rayonnement reçu, il est d'une extrême importance pour la précision de la mesure d'équiper le système de détection de moyens de stabilisation de gain. En effet les systèmes de détection de rayons gamma et notamment les photomultiplicateurs présentent des variations importantes de gain dues surtout aux changements de température ou aux taux de comptage. Un procédé de stabilisation fréquemment utilisé et décrit dans le brevet ci-dessus consiste à adjoindre à l'appareil une source auxiliaire de rayons gamma dont l'énergie est inférieure à la plage du spectre détecté. Dans l'exemple-décrit on a choisi une source d'américium dont le pic d'émis- sion se situe à 60 KeV. On mesure deux taux de comptage Ml et N1 dans deux demi-fenêtres de largeur déterminée, disposées de part et d'autre de l'énergie 60,KeV et l'on utilise un signal d'erreur fonction de la différence (M1 - N1) pour commander par une boucle de contre-réac- tion la haute tension du photomultiplicateur. Cette technique de stabilisation permet de réduire efficace- ment les erreurs de mesure, jusqu'à 5% dans des appareils utilisés actuellement. Toutefois, il est souhaitable d'améliorer encore cette précision. Lorsqu'on utilise pour la stabilisation une source auxiliaire qui émet des rayons gamma monoénergétiques en dessous de la plage d'énergie intéressant la mesure, on détecte le décalage du spectre dans une plage d'énergie o ce décalage est faible. En effet, pour une certaine variation de gain, le décalage du spectre est d'autant plus grand que l'énergie ou se produit ce décalage est élevée. La détection du décalage est donc peu précise et des erreurs apparais- sent, dues aux variations statistiques des taux de comptage. De plus, lorsqu'on utilise des cristaux scintillateurs de grande dimension, les rayons gamma à faible énergie de la source auxiliaire n'atteignent qu'une faible partie du cristal et la stabilisation qui en découle est sérieusement affectée par une hétérogénéité du cristal ou des effets de vieillissement. Une solution intéressante est d'avoir recours, comme dans le brevet américain No. 3 829 686 (W.E. Schultz et autres) à un pic naturel du spectre détecté dont on connalt par avance le niveau d'énergie exact. Cette solution se heurte à une difficulté majeure dans le cas de l'étude des formations géologiques. Le taux de comptage d'un pic déterminé est lié à la teneur, dans la formation géologique examinée, d'un élément particulier. Dans le cas de la spectroscopie des rayons gamma naturels par exemple, le pic centré sur l'énergie 1,46 MeV est presque entièrement dû à la présence de potassium dans la formation. Autrement dit, si la teneur d'un élément devient très faible, le pie correspondant disparaît et toute stabilisation basée sur ce dernier devient trop entachée d'erreurs pour être valable. Dans le brevet des Etats-Unis d'Amérique No. 3 101 409 (L.E. Fite) on a également proposé une stabilisation utilisant deux pics provenant de deux sources auxiliaires de rayons gamma. L'un des pics est utilisé pour commandèr la haute tension du photomultiplicateur tandis que l'autre pic est utilisé pour commander le seuil inférieur du discriminateur de hauteur d'impulsions. Les deux boucles de stabili- sation sont indépendantes et ne permettent pas de résoudre le cas o l'un des pics dispara!trait. Dans le brevet ce problème n'est même pas posé car les pics proviennent de deux sources auxiliaires de rayons gamma monoénergétiques spécialement ajoutées à l'appareil pour la stabilisation. L'objet de l'invention se rapporte à une technique de détec- tion de rayons gamma incluant une stabilisation basée sur la détection de plusieurs pics à différents niveaux d'énergie dont un au moins est situé dans la plage intéressant la mesure. L'objet de l'invention concerne en outre une technique de diagraphie de rayons gamma ayant une détection particulièrement bien stabilisée grâce à l'utilisation optimale des rayons gamma provenant des formations traversées. Selon un aspect de l'invention, un dispositif pour la mesure de rayons gamma provenant des formations traversées par un sondage comprend des moyens de détection pour convertir les rayons gamma reçus en impulsions électriques dont l'amplitude est liée à l'énergie des rayons gamma; et des moyens de discrimination pour trier ces impul- sions électriques suivant leurs amplitudes afin de fournir des infor- mations sur le spectre d'énergie des rayons gamma reçus. Une partie des moyens de discrimination sépare les impulsions électriques dont les amplitudes sont comprises dans deux premières fenêtres contiguës situées de part et d'autre d'une première valeur prédéterminée et dans deux deuxièmes fenêtres contig ués situées de part et d'autre d'une deuxième valeur prédéterminée, ces valeurs prédéterminées cor- respondant à deux pics de référence du spectre des rayons gamma reçus. Au moins un des pics de référence est situé dans le spectre des rayons gamma provenant des formations et donc susceptible de disparaÂtre. Le dispositif comprend en outre des moyens pour combiner les impulsions desdites fenêtres afin d'élaborer un signal d'erreur; et des moyens commandés par ce signal d'erreur pour modifier la réponse aux rayons gamma du dispositif de mesure, afin de stabiliser ladite réponse. De préférence l'autre pic de référence est aussi situé dans le spectre de rayons gamma provenant des formations. Pour élaborer le signal d'erreur, le dispositif comprend des moyens pour fournir les taux de comptage des impulsions comprises dans chacune desdites fenê- tres et des moyens pour combiner les taux de comptage afin de fournir ledit signal d'erreur. Les moyens pour combiner les taux de comptage sont adaptés à calculer les sommes et les différences entre les taux de comptage des fenêtres contigués et.à combiner ces sommes et ces différences pour minimiser les effets des variations statistiques des taux de comptage. Cette combinaison est par exemple le rapport entre une première fonction linéaire des différences et une deuxième fonction linéaire des sommes. Dans ces fonctions linéaires, les sommes et les différences des taux de comptage sont pondérées par des facteurs prédéterminés d'après l'énergie centrale du pic correspondant, la largeur des fenêtres correspondant à ce pic et la résolution des moyens de détection utilisés. De préférence, le dispositif comprend en outre une autre boucle de stabilisation basée sur le pic d'émission d'une source auxi- hlaire de rayons gamma, ce pic étant situé en dehors du spectre détecté des formations. Selon un autre aspect de l'invention, un procédé pour la mesure de rayons gamma provenant des formations traversées par un sondage comprend les étapes suivantes: convertir les rayons gamma reçus en impulsions électriques dont l'amplitude est liée à l'énergie des rayons gamma; et séparer ces impulsions électriques suivant leurs amplitudes afin de fournir des informations sur le spectre d'énergie des rayons gamma reçus. En particulier on sépare les impul- sions électriques dont les amplitudes sont comprises dans deux premiè- res fenêtres contiguIs situées de part et d'autre d'une première valeur prédéterminée et dans deux deuxièmes fenêtres contiguIs situées de part et d'autre d'une deuxième valeur prédéterminée, ces valeurs cor- respondant à deux pics de référence du spectre des rayons gamma reçus. Au moins l'un des pics de référence est situé dans le spectre de rayons gamma provenant des formations. On combine les impulsions des fenêtres afin d'élaborer un signal d'erreur et l'on utilise ce signal d'erreur pour modifier la réponse du dispositif de mesure afin de stabiliser cette réponse. Les caractéristiques et avantages de l'invention ressorti- ront d'ailleurs mieux-de la description qui va suivre donnée à titre d'exemple non limitatif en référence aux dessins annexés dans lesquels: la figure 1 est un schéma d'un dispositif de mesure de rayons gamma selon l'invention représenté dans un sondage; - la figure 2 est un schéma d'un premier mode de réalisation d'une partie du dispositif de la figure 1; la figure 3 est une représentation d'un spectre de rayons gamma reçus provenant des formations; et - la figure 4 est un organigramme de la suite des opérations effectuées par un calculateur dans un autre mode de réalisation qui est une variante de celui de la figure 2. En référence à la figure 1, une sonde de diagraphie 11 est suspendue à l'extrémité d'un câble 12 dans un sondage 13 qui traverse des formations 14. La sonde 11 qui comporte à sa partie supérieure une enveloppe étanche 15 contenant des circuits électroniques fait partie d'un dispositif de diagraphie utilisé pour la spectroscopie des rayons gamma naturels provenant des formations 14. Ce dispositif détermine les teneurs T, U, K en thorium, uranium et potassium des formations par une analyse spectrale de la radioactivité naturelle des formations selon les principes décrits dans le brevet des Etats- Unis d'Amérique No. 3 976 878 déjà mentionné. Tbutefois ce dispositif est équipé d'un nouveau système de stabilisation. Comme on l'a vu précédemment, la stabilisation d'un tel dispositif est très importante car elle conditionne la précision des mesures. Selon l'invention, on détecte des décalages du spectre à plusieurs niveaux d'énergie correspondant à des pics du spectre reçu, on calcule un signal d'erreur fonction de ces décalages et l'on stabilise la réponse du dispositif de mesure par une bouele de contre-réaction qui tend à annuler ce signal d'erreur. Au moins l'un des pics est choisi dans le spectre de rayons gamma provenant des formations et peut donc diminuer ou même disparattre pendant la mesure. Le signal d'erreur est une fonction des décalages qui utilise au mieux, pour la stabilisation, le spectre de rayons gamma provenant des formations. On sait que le décalage d'un pie i du spectre peut %tre détecté au moyenr d'un paramètre appelé facteur de forme F donné par l'équation: Mi - N. F =.M+N (1). i Mi + Ni dans laquelle M et N sont les taux de comptage des rayons gamma dans i i deux fenêtres de largeur Wi de part et d'autre d'une valeur prédéter- minée Ei o devrait se trouver l'énergie moyenne du pic. En supposant que la forme du pic étudié est gaussienne, on trouve que le décalage détecté AEi du pic peut être exprimé en fonction du facteur de forme par l'équation: AEi = Hi Fi (2) dans laquelle H. est, pour un pic, une constante fonction de l'énergie du pic El, de la largeur Wi des fenêtres de détection et de la résolution ri du détecteur à ce niveau d'énergie. Le coefficient Hi peut être prédéterminé pour un pic Ei, pour des fenêtres de largeur Wiet pour un détecteur ayant une résolution ri au niveau Ei, par l'équation: [(E -E)2 Ei+ Wi E)2 2 -1exp(- 2,77E-Eî dE (3) i-ri Ei ri Hî - i o 2Eî[i +Wexp 2, 77Wi] ri i Selon l'invention on élabore un signal d'erreur X qui est une fonction linéaire des facteurs de forme de plusieurs pies. X i i A AE (4) Les coefficients de pondération Ai sont choisis de façon à minimiser l'effet des erreurs statistiques. En effet, ces erreurs statistiques pourraient introduire des incohérences dans la stabilisa- tion notamment lorsque l'amplitude d'un des pies détectés devient très faible, auquel cas le facteur de forme correspondant Fi, voisin de o/o n'a plus de sens. Pour satisfaire à cette condition, les coefficients sont déterminés de façon à minimiser la variance a2 de X. Les équations (2) et (4) donnent: X i Ai Hi i (5) Si la variance de F est Fi, la variance de X est d'après i Fi le théorème de l'addition des variances: a x = i Ai H ai (6) x 2i 2 F Le calcul montre que a2 est minimal pour A 1L 1 Ai 7 2 2 7 iOF, I Fi Comme valeur de a'i on peut prendre: ai2 1 /(M + Ni) (8) Donc: MAi + NI Mi + N Ai Hi= M+ i i (9) Hi i En combinant les équations (1), (5) et (9) on trouve: z (Mi - Ni) i x (10) Ii (Mi+ Ni) i On peut aussi calculer la variance a2 en remplaçant dans l'équation (6) les valeurs de A et a des équations (7) et (8) pour obtenir 2 1/ Mi+N (11) D'après cette dernière équation, on voit que si l'un des pics est absent, la quantité (Mi + Ni) correspondant à ce pic est nulle. En d'autres termes, la contribution de ce pic au signal d'erreur est nulle et la stabilisation est assurée en totalité grâce aux décalages détectés sur les autres pics. Par contre, dès que le pic précédemment absent prend une amplitude notable, la variance a diminue et ce pic participe efficacement à la stabilisation. Dans le mode de réalisation de la figure 1 pour l'étude spectrale de la radioactivité des formations, la stabilisation est réalisée au moyen de deux boucles de contre-réaction dont l'une utilise deux pics du spectre de rayons gamma provenant des formations pour produire un signal d'erreur. La figure 3 représente un exemple de spectre de rayons gamma naturels provenant des formations. Dans ce spectre on remarque plusieurs pics et notamment un pic 16 à 1,46 MeV principalement dû à la présence de potassium dans les formations et un pic 17 à 2,61 MeV principalement dû à la présence de thorium. Dans le mode de réalisation décrit, ces deux pics sont utilisés pour une boucle de stabilisation. En référence de nouveau à la figure 1, la sonde 11 contient de façon classique un cristal à scintillation 20, par exemple au iodure de sodium, suivi par un photomultiplicateur 21 pour convertir les rayons gamma reçus en impulsions électriques dont l'amplitude est proportionnelle à l'énergie des rayons gamma. Ces impulsions sont amplifiées par un amplificateur linéaire 22 puis appliquées à un discriminateur de hauteurs d'impulsions 23. Le discriminateur 23 comprend onze comparateurs 24 dont les seuils de référence sont choi- sis comme représentés en S1 à Sil sur la figure 3. Dans l'exemple choisi les seuils S1 à Sil sont pris respectivement égaux à 200, 500, 1100, 1365, 1460, 1590, 2000, 2520, 2615, 2745 et 3000 KeV. Chaque comparateur ne fournit une impulsion de sortie que pour une impulsion d'entrée dont l'amplitude est supérieure au seuil de ce comparateur. Les onze sorties des comparateurs sont reliées à des circuits logiques de sélection 25 ou circuits d'anticolncidence qui trient les impulsions comprises dans neuf fenêtres différentes. Cinq fenêtres W1 à W5 (voir figure 3) servent aux mesures proprement dites comme indi- qué dans le brevet des Etats-Unis d'Amérique No. 3 976 878 déjà cité. Ces fenêtres de mesure sont celles comprises entre les seuils (S1-S2), (s2-s3) , (S3-S6), (S6-S7) et (S7-Sll). Quatre fenêtres servent à la stabilisation. Ce sont deux premières fenêtres (S4-S5) et (S5-S6) de part et d'autre de la valeur 1460 KeV et deux deuxièmes fenêtres (S8-S9) et (S9-S10) de part et d'autre de la valeur 2615 KeV. Pour chaque fenêtre un circuit logique d'anticolncidence sélectionne les impulsions dont l'amplitude est comprise entre deux seuils. Les impulsions ainsi triées par amplitudes en neuf fenêtres sont appliquées à un circuit de transmission 26 pour être envoyées en surface par l'intermédiaire du cable. Ces circuits 26 peuvent être du type décrit dans le brevet français No. 2 379 694 (inventeurs A. Belaigues et autres) pour transmettre des mots de données ou de com- mande entre le fond et la surface. Les impulsions des neuf fenêtres peuvent être transmises par exemple sous forme de neuf mots binaires représentant chacun le nombre d'impulsions d'une fenêtre pendant une durée de 16 millisecondes (1/60 sec.). Les circuits 22 à 26 sont placés dans l'enveloppe étanche 15. Les informations arrivant en surface sont reçues par d'autres circuits de transmission 27 qui appliquent les neuf taux de comptage à un circuit d'addition et de mémoire 28 qui effectue et met en mémoire la somme des impulsions dans chaque fenêtre pendant un temps t réglable de 1 à 4 secondes. Les taux de comptage dans les cinq fenêtres de mesures sont appliqués à un circuit de traitement 30 qui combine ces taux de comp- tage pour élaborer des signaux représentatifs des teneurs T, U, K en thorium, uranium et potassium des formations, comme expliqué dans le brevet des Etats-Unis d'Amérique No. 3 976 878 déjà mentionné. Le circuit de traitement 30 fournit aussi la radioactivité naturelle to- tale GRtot par une addition des taux de comptage des cinq fenêtres. Ces résultats sont enregistrés au moyen d'un enregistreur magnétique et aussi sous forme graphique sur un film. L'enregistrement est effectué en fonction de la profondeur détectée par une roue 32 en contact sur le câble 12. Les taux de comptage dans les quatre fenêtres de stabilisa- tion sont appliqués à un calculateur de stabilisation 33 qui élabore un signal d'erreur X selon l'équation (10). Ce signal d'erreur X est appliqué aux circuits de transmission de surface 27 qui le transmet par les conducteurs du câble aux circuits de transmission de fond 26. Après démodulation, ce signal X est appliqué à un convertisseur numé- rique-analogique 34 puis à un filtre 35 qui peut être par exemple un filtre RC avec une constante de temps de l'ordre de 100 secondes de façon à rendre négligeables les variations statistiques. La sortie du filtre 35 est reliée à un circuit de soustraction 36 qui reçoit par ailleurs une tension provenant d'une alimentation de référence 37 pour faire varier cette-tension en réponse au signal d'erreur. La tension ainsi commandée détermine par l'intermédiaire d'un diviseur de tension les valeurs de seuils des comparateurs 24. Le circuit 36 est adapté à diminuer proportionnellement les valeurs de seuils lors- que le signal d'erreur X indique un décalage du spectre vers les basses énergies. Le dispositif comprend en outre une deuxième boucle de stabi- lisation, à action rapide, qui utilise une source auxiliaire 40 de rayons gamma. Cette source 40 disposée dans la sonde 11 est en améri- cium 241 avec un pic d'émission centré sur la valeur 60 KeV, c'est-à- dire en dessous du spectre détecté des rayons gamma provenant des formations. Les rayons gamma provenant de la source auxiliaire 40 donnent à la sortie de l'amplificateur 22 des impulsions électriques appliquées à un amplificateur 41 de gain dix. La sortie de l'amplifi- cateur 41 est appliquée à un analyseur de hauteurs d'impulsions 42 formé de trois comparateurs 43 et de circuits logiques de sélection 44. Les comparateurs 43 dont les valeurs de seuils proviennent d'un diviseur de tension alimenté par la tension de référence 37 séparent sur trois sorties les impulsions dont les amplitudes sont respective- ment supérieures à trois seuils S12' S13 et S14 correspondant à des énergies de rayons gamma de 40, 60 et 80 KeV, c'est-à-dire en dessous du spectre détecté provenant des formations (voir figure 3). Les circuits logiques de sélection 44 ont pour fonction de trier sur deux sorties d'une part les impulsions dont l'amplitude est comprise entre les deux seuils S12 et S13 et d'autre part celles dont l'amplitude est comprise entre les deux seuils S13 et S14. On sépare ainsi les impulsions d'un pic d'américium en deux fenêtres situées de part et d'autre de l'énergie centrale 60 KeV. Les deux sorties de la logique 44 sont respectivement reliées aux entrées positive et négative d'un compteur-décompteur 45 dont le contenu est ainsi repré- sentatif de la différence MAm - NAm des taux de comptage des impul- sions comprises dans les deux fenêtres de l'américium. Le contenu du compteur-décompteur 45 est appliqué à un convertisseur numérique- analogique 46 puis à un circuit d'intégration 47 qui fournit un signal d'erreur pour commander la haute tension 48 du photomultiplicateur. Le signal d'erreur est appliqué de façon à augmenter le gain du photomultiplicateur lorsque le pic de l'amérieium se décale vers les basses énergies, c'est-à-dire lorsque le taux de comptage MAm de la fenêtre inférieure devient plus élevé que celui N de la fenêtre Amn supérieure. Le circuit 47 intègre le signal d'erreur avec une cons- tante de temps de quelques secondes. Les circuits 34 à 48 sont dispo- sés dans l'enveloppe étanche 15. La deuxième boucle de stabilisation basée sur le pic d'amé- ricium et agissant sur le gain du photomultiplicateur a l'effet correctif le plus important et le plus rapide (quelques secondes). Cette boucle pourrait aussi être conçue de façon à agir sur le gain de l'amplificateur 22 ou sur la tension de référence 37. La première boucle de stabilisation, basée sur les pics de potassium et de thorium, est d'action plus lente (100 secondes par exemple) et affine l'effet de la deuxième. On obtient ainsi une erreur sur la mesure inférieure au pourcent. La figure 2 représente plus en détail un premier mode de réalisation, sous forme de circuits, du calculateur de stabilisation 33. Ce calculateur reçoit les taux de comptage pendant le temps t $N, MT et NT des impulsions dans les deux fenêtres du potassium et les deux fenêtres du thorium. Les taux de comptage N et NT sont appliqués à deux circuits de complément 50 et 51 qui calculent respectivement les nombres N et -NT. La sortie du circuit 50 est appliquée à un circuit d'addition 52 qui reçoit aussi le taux de comptage M pour fournir des signaux représentatifs de la différence (MK - NK). Un autre circuit d'addition 53 reçoit les signaux -M et NK pour calculer la somme (M. + NK). Des corrections de bruit de fond sont ensuite effectuées sur la différence (M - N et la somme (MK + N.) au moyen de deux circuits 54 et 55.. Ces corrections améliorent la précision obtenue par la stabilisation. En effet, le pic du potassium a un bruit de fond pro- venant à la fois du thorium et de l'uranium. De plus, le taux de comptage de la fenêtre inférieure du potassium est trop élevé du fait du bruit de fond par effet Compton provenant du potassium lui- même. Dans le circuit 54 on calcule la différence Bd des bruits de fond entre les deux fenêtres du potassium par l'équation Bd (BM - BnK) = alT + a2U+K (12) Les coefficients a1, a2 et a5 sont déterminés dans un puits d'étalonnage en plaçant successivement le dispositif en face de plusieurs formations artificielles dont on connaît les teneurs TUK. La différence de bruit de. fond est ensuite retranchée de la différence des taux de comptage. De mâme dans le circuit 55 on calcule la somme Bs des bruits de fond dans les deux fenêtres par l'équation B = b1 T + b U + b K (13) s 2 3 dans laquelle les coefficients b1, b2 et b3 sont déterminés comme précédemment dans un puits d'étalonnage, et l'on retranche cette somme B de la somme des taux de comptage (MK + Ns). La différence corrigée (MK - NK)C est appliquée à un circuit de multiplication 56 pour être multipliée par le coefficient prédé- terminé 1/HK. Comme on l'a vu précédemment, ce coefficient de pondération dépend des largeurs de fenêtre, de l'énergie centrale du pic et de la résolution du détecteur utilisé (équation 3). Dans l'exemple décrit, ce coefficient î/HK est de 24,58 pour un détecteur dont la résolution pour le césium était de 10$. La somme corrigée (4K + NK)C est appliquée à un circuit de multiplication 57 pour être multipliée par un coefficient égal à î/H4. De façon symétrique, le taux de comptage MT et le complément du taux de comptage NT sont appliqués à un circuit d'addition 62 qui calcule la différence MT - NT. Les deux taux de comptage MT et NT sont aJoutés dans un circuit d'addition 63. Le pic de thorium a un bruit de fond dû principalement à l'uranium et l'on peut envisager des corrections de bruit de fond sur la somme et la différence des taux de comptage, de la même façon que pour le pic de potassium. Comme ce bruit de fond est beaucoup plus faible que celui du potas- sium on a choisi, dans le mode de réalisation décrit, d'appliquer directement les pondérations par les facteurs 1/HT et 1/H au moyen de circuits de multiplication 64 et 65. Dans l'exemple, avec les largeurs de fenêtre choisies, le facteur 1/HT a été pris égal à 31,51 pour le même détecteur. Les circuits 56 et 64 sont reliés à un circuit d'addition 66 qui calcule la somme pondérée X1 des différences corri- i gées des taux de comptage: =1 - (Mi - Ni)/Hi = ( - NK)C/HK +(MT - NT) C/HT (14) Les circuits 57 et 65 sont reliés à un circuit d'addition 67 qui calcule la somme pondérée 2 des sommes corrigées des taux de comptage: X = X (Mi + Ni)/H = (M +N)/) +(MT + NT)II-4 (15) Les signaux de sortie des circuits 66 et 67 sont appliqués à un cir- cuit de division 70 qui caleule le rapport Yl/ 12' c'est-à-dire le signal d'erreur X selon l'équation (10). Le circuit de division 70 comprend principalement un multiplicateur 71, un comparateur 72 et un compteurdécompteur 73. Le nombre E1 de l'équation (14) et la sortie du multiplicateur 71 sont appliqués à l'entrée du comparateur 72 dont les sorties sont appliquées au compteur-décompteur 73. De cette façon la sortie A du multiplicateur est égale à A = X. Y2 et le comparateur 72 a pour effet d'égaler A et E1. Par conséquent la sortie du compteurdécompteur 73 est bien égale au rapport: X = E1 / E2 Les circuits de multiplication utilisés peuvent être du type MPY-12HJ commercialisé par la société T.R.W. Californie, Etats-Unis d'Amérique. Dans un autre mode de réalisation, le calculateur 33 peut être un ordinateur universel programmé de façon convenable pour fournir le signal X. Un organigramme général des opérations à effec- tuer est représenté sur la figure 4. îIbut d'abord on lit les taux de comptage M., NK, MT, N et les teneurs T, U, K calculées (bloc 80). On corrige ensuite les taux de comptage pour le bruit de fond en soustrayant par exemple à chaque taux de comptage un nombre de comptes calculé par une fonction linéaire des teneurs T, U, K (bloc 81). On calcule ensuite le signal d'erreur X (bloc 82) par l'équation (10) et l'on extrait cette valeur calculée du signal d'erreur pour l'appliquer aux circuits de transmission 27. Ce calcul est effectué en temps réel pour chaque niveau de profondeur auquel se trouve la sonde 11 comme représenté par les blocs 84 et -85. Ce mode de réalisa- tion est équivalent à celui de la figure 2. De la même façon les taux de comptage MA et N du pic Am Am d'américium pourraient tre transmis en surface et le signal d'erreur pour la stabilisation correspondante pourrait être calculé en surface par un ordinateur programmé puis retransmis au fond pour commander le gain du dispositif de mesure. On remarquera que pour chaque pic, potassium ou thorium, les fenêtres de détection contigués sont de largeur inégale. Dans l'exemple ces largeurs sont de 95 KeV pour la fenêtre inférieure et KeV pour la fenêtre supérieure. On choisit la fenêtre supérieure plus large afin de tenir compte de la diminution de résolution du détecteur lorsque l'énergie augmente. De nombreuses variantes peuvent 'tre apportées aux modes de réalisation décrits sans sortir du cadre de l'invention. En parti- culier le signal d'erreur peut ttre calculé à partir d'un nombre de pics de référence supérieur à deux, en utilisant la formule (10). Au _0 lieu de deux boucles de stabilisation, on peut combiner les taux de comptage des rayons gamma provenant de la source auxiliaire monoéner- gétique et les taux de comptage d'un ou plusieurs pics du spectre détecté provenant des formations. Le filtrage réalisé par le circuit pourrait tre effectué par une étape de programme ajoutée à l'organigramme de la figure 4. REVENDICATIONS 1. Dispositif pour la mesure de rayons gamma provenant des formations traversées par un sondage comprenant: - des moyens de détection pour convertir les rayons gamma reçus en impulsions électriques dont l'amplitude est liée à l'éner- gie des rayons gamma; - des moyens de discrimination pour trier ces impulsions électriques suivant leurs amplitudes afin de fournir des informa- tions sur le spectre d'énergie des rayons gamma reçus, lesdits moyens de discrimination comprenant des moyens pour séparer les impulsions électriques dont les amplitudes sont comprises dans deux premières fenêtres contigues situées de part et d:'autre d'une première valeur prédéterminée et dans deux deuxièmes fenêtres con- tigdés situées de part et d'autre d'une deuxième valeur prédéter- minée, ces valeurs prédéterminées correspondant à deux pics de référence du spectre des rayons gamma reçus; caractérisé en ce qu'au moins un des pics de référence étant situé -dans le spectre détecté de rayons gamma provenant des formations et par conséquent susceptible de disparaître, le dispositif comprend en outre - des moyens pour combiner les impulsions desdites pre- mières et deuxièmes fenêtres afin d'élaborer un signal d'erreur; et - des moyens commandés par ledit signal d'erreur pour modifier la réponse du dispositif de mesure afin de stabiliser ladite réponse. 2. Dispositif de mesure selon la revendication 1, caracté- risé en ce que l'autre pic de référence est situé dans le spectre de rayons gamma provenant des formations. 3. Dispositif de mesure selon l'une des revendications 1 et 2, caractérisé en ce que les moyens pour élaborer le signal d'erreur comprennent des moyens pour fournir les taux de comptage des impul- sions comprises dans chacune desdites fenêtres et des moyens pour combiner les taux de comptage afin de fournir ledit signal d'er- reur. 4. Dispositif de mesure selon la revendication 3, caracté- risé en ce que les moyens pour combiner les taux de comptage sont adaptés à calculer les différences entre les taux de comptage des fenêtres contigues; et à combiner lesdites différences de façon à minimiser les effets des variations statistiques desdits taux de comptage. 5. Dispositif de mesure selon la revendication 4, caracté- risé en ce que les moyens pour combiner les taux de comptage sont adaptés à: calculer en outre les sommés des taux de comptage dans les fenêtres contigues; calculer une première fonction linéaire desdites différences, calculer une deuxième fonction linéaire desdites sommes; et calculer le rapport de la première à la deu- xième fonction afin d'obtenir le signal d'erreur. 6. Dispositif de mesure selon la revendication 5, caracté- risé en ce que lesdits moyens pour combiner les taux de comptage sont adaptés à calculer ladite. première fonction en effectuant une somme pondérée desdites différences des taux de comptage et ladite deuxième fonction en effectuant une somme pondérée desdites sommés des taux de comptage. 7. Dispositif de mesure selon la revendication 6, caracté- risé en ce que lesdits moyens pour combiner les taux de comptage sont adaptés à pondérer chacune desdites différences des taux de comptage par un facteur prédéterminé d'après l'énergie centrale du pic de référence correspondant, la largeur des fenêtres corres- pondant à ce pie de référence et la résolution des moyens de détection utilisés. 8. Dispositif de mesure selon la revendication 7, caracté- risé en ce que lesdits moyens pour combiner les taux de comptage sont adaptés à pondérer chacune desdites sommes des taux de comp- tage par un deuxième facteur prédéterminé égal au carré du facteur de pondération appliqué à la différence correspondante des taux de comptage. 9. Dispositif de mesure selon l'une des revendications 3 à 8, caractérisé en ce que les moyens pour élaborer le signal d'erreur comprennent des moyens pour corriger les taux de comptage pour tenir compte du bruit de fond. 10. Dispositif de mesure selon l'une des revendications 1 à 8, caractérisé en ce qu'il comprend de plus une source auxiliaire pour émettre des rayons gamma monoénergétiques en dehors de la plage étudiée de rayons gamma provenant des formations; et des moyens sensibles aux rayons gamma émis par ladite source auxiliaire pour modifier la réponse dudit dispositif d'analyse afin d'effectuer une deuxième stabilisation de ladite réponse. 11. Dispositif de mesure selon la revendication 10, carac- térisé en ce que les moyens de discrimination comprennent des moyens pour séparer les impulsions électriques dont les amplitudes sont comprises dans deux troisièmes fenêtres contiguës situées de part et d'autre d'une troisième valeur déterminée correspondant au pic de rayons gamma émis par la source auxiliaire; et en ce que lesdits moyens pour effectuer la deuxième stabilisation comprennent des moyens pour combiner les impulsions desdites troisièmes fenêtres afin d'élaborer un deuxième signal d'erreur et des moyens commandés par le deuxième signal d'erreur pour modi- fier la réponse du dispositif d'analyse. 12. Dispositif de mesure selon la revendication 11, carac- térisé en ce que les moyens commandés par le deuxième signal d'erreur sont adaptés à modifier le gain desdits moyens de détec- tion. 13. Dispositif de mesure selon la revendication 12, carac- térisé en ce que les moyens commandés par le premier signal d'erreur sont adaptés à modifier certains seuils'de comparaison desdits moyens de discrimination de façon que la modification desdits seuils n'ait pas d'influence sur ledit deuxième signal d'erreur. 14. Procédé pour la mesure de rayons gamma provenant des formations traversées par un sondage comprenant les étapes sui- vantes: - convertir les rayons gamma reçus en impulsions électri- ques dont l'amplitude est liée à l'énergie des rayons gamma; - trier ces impulsions électriques suivant leurs amplitu- des afin de fournir des informations sur le spectre d'énergie des rayons gamma reçus, et notamment séparer les impulsions électriques dont les amplitudes sont comprises dans deux premières fenêtres contiguXés situées de part et d'autre d'une première valeur prédéter- minée et dans deux deuxièmes fenêtres contigués situées de part et d'autre d'une deuxième valeur prédéterminée, ces valeurs prédé- terminées correspondant à deux pics de référence du spectre des rayons gamma reçus; caractérisé en ce qu'au moins un des pics de référence étant situé dans le spectre détecté de rayons gamma provenant des formations et par conséquent susceptible de disparaître, le procédé comprend en outre les étapes suivantes: - combiner les impulsions desdites premières et deuxiè- mes fenêtres afin d'élaborer un signal d'erreur; et - modifier au moins l'une des étapes de mesure au moyen dudit signal d'erreur afin de stabiliser ladite mesure. 15. Procédé de mesure selon la revendication 14, caractérisé en ce que l'autre pic de référence est situé dans le spectre de rayons gamma provenant des formations. 16. Procédé de mesure selon l'une des revendications 14 et , caractérisé en ce que l'étape d'élaboration du signal d'erreur comprend les étapes suivantes: - produire les taux de comptage des impulsions compri- ses dans chacune desdites fenêtres; - calculer les sommes et les différences des taux de comptage des fenêtres contiguIs; et - calculer le rapport entre une première fonction liné- aire desdites différences et une deuxième fonction linéaire des- dites sommes afin d'obtenir le signal d'erreur. a48S752 17. Procédé de mesure selon la revendication 16, caractérisé en ce que l'on calcule ladite première fonction linéaire en effec- tuant une somme pondérée desdites différences des taux de comptage et ladite deuxième fonction linéaire en effectuant une somme pon- dérée desdites sommes des taux de comptage. 18. Procédé de mesure selon la revendication 17, caractérisé en ce que les facteurs de pondération utilisés pour calculer la première fonction linéaire sont des facteurs prédéterminés dépen- dant de l'énergie-centrale du pic de référence correspondant, de la largeur des fenrtres correspondant à ce pic de référence et de la résolution obtenue pour convertir les rayons gamma reçus en impulsions. 19. Procédé de mesure selon la revendication 18, caractérisé en ce que les facteurs de pondération utilisés pour calculer la deuxième fonction linéaire sont pris égaux aux carrés des facteurs de pondération utilisés pour calculer la première fonction linéai- re. 20. Procédé de mesure selon l'une des revendications 16 à 19, caractérisé en ce qu'il comprend de plus une étape consistant à corriger les taux de comptage pour tenir compte du bruit de fond. 21. Procédé de mesure selon l'une des revendications 15 à , caractérisé en ce qu'il comprend de plus les étapes suivantes - émettre des rayons gamma monoénergétiques en dehors de la plage étudiée de rayons gamma provenant des formations; et - modifier au moins l'une des étapes de mesure d'après les rayons gamma monoénergétiques détectés afin d'effectuer une deuxième stabilisation de ladite mesure. 22. Procédé de mesure selon la revendication 21, caractérisé en ce qu'il comprend les étapes suivantes: - séparer les impulsions électriques dont les amplitudes sont comprises dans deux troisièmes fenêtres contigues situées de part et d'autre d'une troisième valeur déterminée correspondant au pic de rayons gamma émis par la source auxiliaire; et - combiner les impulsions desdites troisièmes fenêtres afin d'élaborer un deuxième signal d'erreur pour effectuer ladite deuxième stabilisation. 23. Procédé de mesure selon la revendication 22, caracté- risé en ce qu'il comprend l'étape consistant à modifier, par le deuxième signal d'erreur, ladite étape de conversion des rayons gamma en impulsions électriques. 24. Procédé de mesure selon la revendication 23, caractérisé en ce qu'il comprend l'étape consistant à modifier, par le premier signal d'erreur, certains seuils de comparaison de ladite étape -de tri des impulsions électriques.