Brevet FR2818427A1 — publié le 20020621

L'invention concerne un procédé de détermination du temps de doublement d'un flux neutronique dans le coeur d'un réacteur nucléaire. Le temps de doublement du flux neutronique dans le coeur d'un réacteur nucléaire est égal à la période de temps mise par le flux pour être multiplié par deux. Ce temps de doublement est donc un indicateur d'une vitesse de croissance du flux à un instant donné qui sert à la surveillance de la puissance du réacteur, en particulier lors de transitoires de fonctionnement normal, incidentel ou accidentel du réacteur nucléaire tels que des démarrages, montées en puissance, baisses de puissance, arrêt normal ou arrêt d'urgence. Dans le cas des réacteurs nucléaires comportant une cuve renfermant le coeur du réacteur nucléaire, et en particulier dans le cas des réacteurs nucléaires refroidis par de l'eau sous pression, le flux neutronique peut être mesuré par des capteurs disposés à l'extérieur de la cuve, à proximité d'une paroi de la cuve dans une zone renfermant le coeur du réacteur nucléaire. De tels capteurs sont généralement réalisés sous la forme de chambres"excore", c'est-à-dire situées à l'extérieur de la cuve, qui comportent un dépôt d'une substance susceptible de libérer des particules sous l'effet d'un bombardement électronique, telle que le bore et qui sont remplies d'un gaz inerte tel que de l'argon. Sous l'effet du bombardement des neutrons provenant du coeur du réacteur nucléaire qui s'échappent de la cuve, le bore libère des particules alpha et l'argon remplissant la chambre se trouve ionisé. Le passage des particules alpha+ et e- entre deux électrodes permet de compter le nombre de particules libérées par unité de temps qui est représentatif du flux neutronique à l'intérieur du coeur du réacteur nucléaire. Les mesures des capteurs sont fournies sous la forme d'un nombre de coups par seconde correspondant au nombre de passages de particules par seconde entre les électrodes. Le nombre de coups par seconde ou taux de comptage Tc est fourni par une carte électronique du système de contrôle du réacteur nucléaire appelée carte ictomètre qui compte le nombre d'impulsions qu'elle reçoit et calcule le taux de comptage en fonction d'un temps de comptage qui est déterminé par le programme de la carte ictomètre. Le taux de comptage est utilisé pour le calcul du temps de doublement Td du flux neutronique qui est exprimé en seconde. Les procédés de calcul utilisés jusqu'ici pour déterminer le temps de doublement à partir des valeurs mesurées des taux de comptage se traduisent par l'obtention d'un signal qui présente de fortes fluctuations au cours du temps. De ce fait, l'opérateur chargé de la surveillance de l'évolution du flux neutronique dans le coeur du réacteur nucléaire depuis la salle de commande du réacteur dispose d'une information sous forme graphique qu'il est très difficile d'interpréter, ce qui complique les tâches de surveillance du réacteur nucléaire. Il est donc apparu nécessaire d'améliorer le signal représentatif du temps de doublement du flux neutronique, pour que celui-ci soit plus facilement exploitable par un opérateur. Pour obtenir un signal plus facilement exploitable, on a proposé d'utiliser une méthode de calcul du temps de doublement par moyenne glissante mais le résultat n'est pas très satisfaisant et ce procédé de calcul par moyenne glissante présente l'inconvénient de nécessiter une très grande capacité de mémoire dans un calculateur. Le but de l'invention est donc de proposer un procédé de détermination du temps de doublement d'un flux neutronique dans le coeur d'un réacteur nucléaire comportant une cuve renfermant le coeur du réacteur, dans lequel on dispose, à l'extérieur et à proximité de la cuve du réacteur nucléaire, au moins une chambre de mesure comprenant un revêtement en une substance capable d'émettre des particules sous l'effet d'un flux neutronique, de manière à soumettre la substance au flux neutronique dans le coeur du réacteur, à travers une paroi de la cuve du réacteur nucléaire, on compte le nombre de particules émises par unité de temps, à des intervalles de temps At pendant le fonctionnement du réacteur nucléaire et on détermine le temps de doublement du taux de comptage par le calcul, ce procédé étant perfectionné de manière qu'on obtienne un signal représentatif du temps de doublement parfaitement stable et facilement utilisable, avec des moyens de calcul de capacité modérée. Dans ce but : - on détermine, pour chaque valeur du taux de comptage Tc, l'intervalle de temps At séparant la mesure du taux de comptage Tc d'une mesure précédente du taux de comptage, - on optimise les paramètres d'un algorithme de calcul utilisant un filtre de Kalman, à partir de l'intervalle de temps At déterminé précédemment, de paramètres du filtre calculés dans une étape précédente de détermination du temps de doublement et du temps de doublement Td déterminé précédemment, - on calcule, au moyen de l'algorithme utilisant le filtre de Kalman avec ses paramètres optimisés, une nouvelle valeur du temps de doublement, et - on met en mémoire les paramètres optimisés et la valeur du temps de doublement. Afin de bien faire comprendre l'invention, on va maintenant décrire, en se référant aux figures jointes en annexe, un exemple de mise en oeuvre du procédé suivant l'invention dans le cas d'un réacteur nucléaire à eau sous pression dont le flux neutronique est évalué par des chambres excore. La figure 1 est une vue schématique en plan de la cuve d'un réacteur nucléaire renfermant le coeur et de chambres excore disposées à l'extérieur et au voisinage de la paroi de la cuve utilisées pour déterminer le temps de doublement du flux neutronique dans le coeur du réacteur nucléaire, au cours d'une phase de fonctionnement du réacteur. La figure 2 est une courbe donnant le taux de comptage Tc en fonction du temps, pendant la phase de fonctionnement du réacteur nucléaire. La figure 3 est un diagramme donnant le temps de doublement Td du flux en fonction du temps déterminé par le procédé suivant l'invention et, de manière comparative, par un procédé suivant l'art antérieur à moyenne glissante. La figure 4 est un schéma montrant le principe du procédé de calcul du temps de doublement selon l'invention. La figure 5 est un diagramme fonctionnel d'un filtre de Kalman utilisé dans le cadre de la mise en oeuvre du procédé suivant l'invention. Sur la figure 1, on a représenté de manière schématique la paroi d'une cuve 1 d'un réacteur nucléaire à eau sous pression qui présente une forme générale cylindrique et qui renferme le coeur 2 du réacteur nucléaire constitué par des assemblages de combustible de forme prismatique droite placés de manière juxtaposée à l'intérieur de la cuve du réacteur nucléaire, avec leurs axes dans la direction axiale de la cuve. Le flux dans le réacteur nucléaire est évalué à l'aide de chambres excore 4 qui sont réparties autour de la cuve du réacteur nucléaire, de manière à fournir des mesures ou évaluations de flux neutronique pour différentes zones du coeur. Par exemple, on peut utiliser quatre chambres excore 4 disposées au voisinage de la paroi 3 de la cuve 1 et placées à 900 l'une de l'autre autour de l'axe de la cuve. Dans ce cas, on obtient, sur chacune des chambres excore 4, des mesures représentatives du flux neutronique dans un quart ou quadrant du coeur du réacteur nucléaire. Les chambres excore 4 qui peuvent être multi-étagées comportent chacune une enveloppe s'étendant suivant une partie au moins de la hauteur du coeur dans la direction axiale de la cuve qui est revêtue par une couche d'une substance telle que le bore susceptible d'émettre des particules alpha sous l'effet du bombardement neutronique provenant du coeur 2 du réacteur nucléaire et traversant la paroi 3 de la cuve. Les chambres excore 4 sont remplies d'un gaz neutre tel que l'argon qui est ionisé sous l'effet du faisceau de neutrons et des particules émises par la couche de bore. Les chambres excore 4 comportent des électrodes entre lesquelles passent les particules alpha et les électrons e- produits par le faisceau neutronique. Les chambre excore 4 sont reliées à une unité de traitement 5 comportant une carte ictomètre permettant d'effectuer le comptage par unité de temps des particules traversant les électrodes des chambres excore, ce comptage se traduisant par un nombre de coups par seconde qui est représentatif du flux neutronique reçu par la chambre excore. Les évaluations du flux à partir du taux de comptage Tc en coups par seconde peuvent être utilisées pour chacune des chambres excore 4 ou en- core sous la forme d'une moyenne des résultats relatifs à chacune des chambres excore. L'enregistrement du taux de comptage relatif à une chambre excore ou à l'ensemble des chambres excore, pendant une phase de fonctionnement du réacteur nucléaire, peut être représenté sous la forme d'un diagramme tel que le diagramme de la figure 2. On voit que, pendant une phase de fonctionnement du réacteur nucléaire d'une durée de l'ordre de 80 secondes (prise à titre d'exemple), le taux de comptage est resté pratiquement stable et égal à 2000 coups par seconde pendant une durée de 60 secondes puis a augmenté de la valeur de 2000 coups par seconde à une valeur proche de 10000 coups par seconde pendant une durée de 20 secondes. Le diagramme de la figure 2 montre donc l'évolution globale du taux de comptage et ne permet pas de déterminer directement la vitesse de variation instantanée du flux neutronique. Sur la figure 3, on a représenté, sous la forme de la courbe 6 en traits fins, le temps de doublement du flux neutronique dans le coeur du réacteur nucléaire calculé par un procédé suivant l'art antérieur à moyenne glissante, avec un temps de cycles de calcul de 300 ms. Du fait des fluctuations au cours du temps du taux de comptage, ces fluctuations étant pratiquement aléatoires sur une durée égale au temps de cycle de calcul, la courbe 6 représentative du temps de doublement calculé par moyenne glissante est fortement perturbée et ne permet pas de déterminer facilement la vitesse de variation du flux neutronique à un instant donné. En particulier, les valeurs du temps de doublement qui traduisent le bruit de la mesure du taux de comptage (voir figure 2) passent de valeurs positives à des valeurs négatives suivant les variations instantanées du taux de comptage. A titre de comparaison, la courbe 7 obtenue par le procédé suivant l'invention qui sera décrit ci-après permet d'obtenir une représentation par- faitement utilisable de l'évolution du temps de doublement, pendant le fonctionnement du réacteur nucléaire. Ce temps de doublement est pratiquement constant et égal à 300 secondes pendant la première minute de fonctionnement du réacteur nucléaire puis décroît rapidement jusqu'à une valeur très faible, proche de zéro, ce qui traduit une augmentation très rapide du flux neutronique. Cette évolution rapide du flux neutronique est parfaitement apparente à un opérateur en salle de commande disposant d'une représentation graphique du temps de doublement, sur un écran. Selon l'invention, on utilise une méthode de calcul optimisée du temps de doublement au moyen d'un algorithme utilisant un filtre de Kalman, de manière à limiter au maximum le bruit de la mesure. La courbe du temps de doublement ainsi obtenue (courbe 7) de la figure 3 peut être ainsi facilement interprétée par l'opérateur dans la salle de commande. Description d'un filtre de Kalman Soit un système linéaire de vecteur d'état X (pour le calcul du temps de doublement, les composantes du vecteur d'état doivent comprendre le temps de doublement ou une expression permettant de le calculer, ainsi que le taux de comptage). On cherche à connaître les valeurs que prend X aux instants ti, t2, t3, On note ces valeurs X (1), X (2), X (3),.... On suppose que l'évolution de X obéit à la relation : X (k+1) = On suppose que U (k) est un bruit blanc centré (c'est-à-dire que U (k) est de moyenne nulle et est indépendant de U (j) si k est différent de j). La matrice de covariance de U (k) est notée Q (k). On suppose ensuite qu'aux instants ti, t2, tus,.... on effectue les mesures y (1), y (2), y (3), etc. et que y est une fonction linéaire de X : y (k) = M (k) + w (k) où : M (k) est la matrice de mesure (ou d'observation), et w (k) le bruit de mesure, supposé être un bruit blanc centré de matrice de covariante R (k). Le filtre de Kalman est un algorithme utilisé pour donner une estimation de X (k) à partir des mesures y. Le coeur de l'algorithme du filtre de Kalman consiste à calculer la matrice de covariance de l'erreur dans l'estimation de X. Cette matrice est appelée P. Elle est calculée à chaque nouvelle observation faite sur la base des mesures. Application du filtre de Kalman au calcul du temps de doublement Le taux de comptage Tc et le temps de doublement Td sont supposés être liés par la relation : Tc (t) = Tco 2trad Tc (t) = Tco exp [ (ln 2) t/Td] où : exp est la fonction exponentielle de base e, et In 2 le logarithme népérien de 2. Une évolution linéaire est obtenue en prenant le logarithme népérien de l'expression : In Tc (t) = In Tco + (In 2) t/Td. Xi 1 On prend le vecteur d'état X avec X2 xi (k) = In TC (tk) x2 (k) = (In 2)/Td (tk). En considérant le temps de doublement comme constant entre deux acquisitions successives du taux de comptage (c'est-à-dire que Td est à variation lente par rapport à la fréquence d'acquisition des observations), on a : xi (k+1) = xi (k) + x2 (k) At (k) + u1 (k) x2 (k+1) =X2 (k) + u2 (k) où : At est la durée entre deux mesures rafraîchies et prises en compte, et u1 et U2 des bruits blancs centrés. Dans le cas où Tc n'est pas rafraîchi, aucune nouvelle valeur n'est calculée mais on effectue une mémorisation de la dernière valeur. Comme c'est le taux de comptage qui est mesuré, l'observation peut s'écrire : y (k) =xi (k) + w (k) où x est le bruit de mesure, bruit blanc centré non corrélé avec Ui et U2. L'application du filtre de Kalman se particularise donc en prenant les matrices : X (k)) le vecteur d'état, (x (k) j n At (k) C (k) = matrice de changement d'état, (0 1 1 U (k)) le bruit traduisant le caractère aléatoire du phénomène U (k) j physique, y (k) la mesure ou l'observation (matrice 1x1), M (k) = (1 0) la matrice de mesure (constante) w (k) le bruit de mesure. L'algorithme de calcul utilisant le filtre de Kalman permet de fournir, en sortie, la valeur X2 = (In 2) /Td. On en déduit, par le calcul, 1/Td = X2/tn 2 et la valeur instantanée du temps de doublement du flux neutronique Td. Cependant, comme il est visible en particulier sur la figure 4, il est nécessaire de réaliser différents calculs préalables et une optimisation des paramètres, avant d'effectuer le calcul du temps de doublement du flux neutronique par l'intermédiaire du filtre de Kalman 8 à un instant donné, à partir de la valeur x2 (k). La carte ictomètre qui fournit le taux de comptage Tc est programmée de manière à fournir ce taux de comptage selon une fréquence qui varie avec le taux de comptage Tc. La période de temps séparant deux déterminations du taux de comptage, appelée période de rafraîchissement des taux de comptage T, cto, varie selon la règle suivante : Pour Tc de 0 à 100 coups par seconde : T, cto = 20, 48 secondes. Pour Tc compris entre 100 et 204800 coups par seconde : T, cto varie de 20,48 secondes à 10 millisecondes, selon la formule : T. cato = 2048/Tc. Pour Tc compris entre 204800 et 6,5. 106 coups par seconde : T, cto = 10 millisecondes. La durée du cycle de traitement de l'unité centrale affectée au calcul du temps de doublement (cette unité est schématisée à la figure 4) est de 300 millisecondes. Si T, cato Dans le cas où Lcto est supérieur à 300 ms, on ne peut effectuer un nouveau calcul avant qu'une nouvelle valeur ne soit donnée par la carte ictomètre. La détermination de l'intervalle de temps At (k) correspondant à l'intervalle de temps de rafraîchissement de la carte ictomètre est donc définie de la manière suivante par l'unité de calcul 9, à partir de la valeur T du nombre d'unités de temps des mesures (625 microsecondes). L'unité 9 calcule, à partir de T (k), T, cato = T x 0,000625, ce qui donne le temps de rafraîchissement en seconde. Ce temps T, cto est comparé à 300 millisecondes. Dans le cas où T, cto > = 300 ms : At = Tcto. Dans le cas où T, cato La valeur de At retenue est entrée dans l'unité de calcul 10 permettant l'optimisation dynamique des paramètres de calcul. L'unité 10 reçoit les paramètres obtenus au cycle précédent et gardés en mémoire et en particulier la valeur du temps de doublement précédemment calculée Td (k-1). Dans l'unité de calcul 11, on réalise, à partir de la valeur instantanée Tc (k) du taux de comptage, le calcul de la valeur y (k) = In (Tc). L'algorithme de calcul 8 utilisant un filtre de Kalman permet de calculer la valeur x2 (k) à partir de la valeur y (k). A partir de la valeur x2 (k), on calcule, dans l'unité de calcul 12, la valeur de Td (k) = ln 2/x2 (k). La valeur instantanée du temps de doublement Td (k) ainsi que les paramètres optimisés de l'algorithme de calcul sont introduits dans la mémoire 13, pour être fournis au cycle suivant à l'unité 10 d'optimisation dynamique des paramètres. L'algorithme de calcul de la valeur courante de l'inverse du temps de doublement 1/Td~n est optimisé par l'ajustement de paramètres qui sont les termes de deux matrices Q et R. Les termes de la matrice Q qui représente la covariance du bruit blanc centré de l'évolution de X peuvent être fixés à 1, par exemple : q11 = q12 = q21 = q22 = 1. La matrice R représente la covariance du bruit de mesure et peut être déterminée à partir de l'intervalle de temps de rafraîchissement At calculé précédemment et de l'inverse de la valeur instantanée précédente du temps de doublement courant Td~n par la formule : R = a/ (At 11/Td~npréc 13). Les données d'entrée en provenance de la mémoire 13 utilisée pour effectuer l'ajustement des paramètres sont constituées, en plus de l'inverse de la valeur courante précédente du temps de doublement, par les valeurs xi et X2 précédentes et les termes pu, P12 et P22 de la matrice P de cova- riance de l'erreur dans l'estimation de X, déterminés lors d'un cycle de calcul précédent. Les valeurs d'initialisation utilisées pour xi et X2 sont les suivantes : xi (0) = ! nTco=0 (Tc= 1) X2 (0) = In 2/Td (O) = 0,002 (Td > 300 s, flux stable). Comme indiqué plus haut, R représente la variance du bruit de mesure ; en conséquence, plus R est grand, plus le filtre est lent et plus les perturbations sont atténuées. Inversement, plus R est petit et plus le filtre est rapide mais moins le bruit est réduit. R joue donc un rôle sur le rapport rapidité/précision du filtre. Pour augmenter la rapidité du filtre à faible taux de comptage (le nombre de points délivrés par la carte ictomètre étant plus faible), on choisit R proportionnel à Ri avec Ri = 1/At. Pour Q et R fixes et pour un bruit d'amplitude proportionnel au taux de comptage, on constate que : - pour un temps de doublement élevé en valeur absolue, il est nécessaire d'augmenter R pour obtenir une précision correcte, - pour un temps de doublement faible en valeur absolue, il est nécessaire de diminuer R pour obtenir une réponse rapide. On ajuste donc R (avec Q fixe) en fonction du temps de doublement Td précédemment calculé. R est proportionnel à R2 = 1/11/Td~préc 13. En définitive, en adoptant le coefficient de proportionnalité a, on retrouve la formule de calcul de R covariance du bruit de mesure donnée plus haut. On rappellera dans le tableau 1 ci-dessous les différentes étapes du calcul : TABLEAU 1 T, cto = T * 0, 000625 Conversion de la donnée d'entrée T en seconde #t = Ticto si Ticto > = 300 ms Calcul de At, le temps entre deux mesures #t = 0,3 si Ticto 00 ms rafraîchies et prises en compte par l'UC R = &alpha;/(#t # 1 / Td~npréc.]3) Calcul de R, covariance du bruit de mesure y = ln(Tc) Calcul de y, l'observation au cycle courant #1 = x1préc. + #tx2préc. Calcul de la première estimation du vec- #2 = x2préc. teur d'état X (équation 1 du filtre de Kalman) #11 = p11préc. + 2 #tp12péc. + #t2p22préc + q11 Calcul de la covariance de l'erreur de la #12 = p12préc. + #t p22préc + q12 première estimation, P p22 = p22préc = q22 k1 = #11 / (#11 + R) Calcul du gain du filtre, K k2 = #12 / (#11 + R) x1 = #1 + k1 (y - #1) Calcul de la deuxième estimation de X x2 = #2 + k1 (y - #1) (nouvelle valeur du vecteur d'état X) p11 = k1 R Calcul de la covariance P de l'erreur de la P12 = k2 R nouvelle estimation de X P22 = P22-k2 P12 1/Td~n = x2 / ln 2 si x2 / ln 2 > 1/300 Calcul de la nouvelle valeur de 1/dn avec 1/Td~n = 1/300 sinon écrêtage Xlpréc = Xi X2préc = X2 Pupréc. = Pli Mémorisation des variables Xi, X2, pu, P12, P12préc = P12 P22 et 1/Td~n pour le cycle suivant P22préc = P22 1/Td n réc= 1/Td n Sur la figure 5, on a représenté un diagramme fonctionnel du filtre de Kalman 8 utilisé pour effectuer le calcul du vecteur d'état X (k) au cycle cou- rant à partir de la valeur mesurée y (k). La signification des différentes notations utilisées sur la figure 5 est donnée ci-dessous : y (k) : mesure ou observation au cycle courant (k), X (kalk) : estimation du vecteur d'état au cycle courant, sur la base des mesures y (1), y (2),..., y (k), X (k-1 lk-1) : estimation du vecteur d'état au cycle précédent, sur la base des mesures y (1), y (2),..., y (k), X (kl-1) estimation du vecteur d'état au cycle courant, sur la ba- se du vecteur d'état précédent, Mém. : mémorisation du vecteur d'état d'un cycle à l'autre, (D (k-1) : matrice de changement d'état, M (k) : matrice de mesure ou d'observation, K (k) : gain du filtre de Kalman. En définitive, le procédé de calcul utilisant un filtre de Kalman permet d'obtenir des valeurs et une représentation graphique du temps de doublement qui peuvent être facilement exploitées par un opérateur chargé de la surveillance du flux neutronique dans le coeur du réacteur nucléaire. La représentation graphique des variations du temps de doublement du flux neutronique dans le coeur d'un réacteur nucléaire est représentée par la courbe 7 sur la figure 3, dans le cas de la phase de fonctionnement d'un réacteur nucléaire qui a été décrite plus haut. REVENDICATIONS 1.-Procédé de détermination du temps de doublement d'un flux neu- tronique dans le coeur d'un réacteur nucléaire comportant une cuve (1) renfermant le coeur (2) du réacteur nucléaire, dans lequel on dispose, à l'extérieur et à proximité de la cuve du réacteur nucléaire, au moins une chambre de mesure (4) comprenant un revêtement en une substance capable d'émettre des particules sous l'effet d'un flux neutronique, de manière à soumettre la substance au flux neutronique dans le coeur du réacteur nucléaire à travers une paroi (3) de la cuve (1) du réacteur nucléaire, on compte le nombre de particules émises par unité de temps, pour obtenir, à des intervalles de temps (At) pendant le fonctionnement du réacteur nucléaire, une valeur mesurée appelée taux de comptage (Tc) et on détermine le temps de doublement du taux de comptage par le calcul, caractérisé par le fait : - qu'on détermine, pour chaque valeur du taux de comptage (Tc) l'in- tervalle de temps (At) séparant la mesure du taux de comptage d'une mesure précédente du taux de comptage (Tc), - qu'on optimise les paramètres d'un algorithme de calcul utilisant un filtre de Kalman, à partir de l'intervalle de temps (At) déterminé précédemment, de paramètres du filtre calculés dans une étape précédente de détermination du temps de doublement et du temps de doublement (Td) déterminé précédemment, - qu'on calcule, au moyen de l'algorithme utilisant le filtre de Kalman avec ses paramètres optimisés, une nouvelle valeur du temps de doublement (Td) et qu'on met en mémoire les paramètres optimisés et la nouvelle valeur du temps de doublement (Td). 2.-Procédé suivant la revendication 1, caractérisé par le fait qu'on prend comme vecteur d'état du filtre de Kalman le vecteur (X) dont les com- posantes sont (xi (k) = In Tc (tk)) et (x2 (k) = (In 2)/Td (tk)). 3.-Procédé suivant l'une quelconque des revendications 1 et 2, caractérisé par le fait qu'on optimise les paramètres à partir d'une matrice (Q) de covariance d'un bruit aléatoire relatif à l'évolution du vecteur d'état (X) et de la covariance (R) du bruit de la mesure. 4.-Procédé suivant la revendication 3, caractérisé par le fait qu'on choisit tous les termes de la matrice (Q) égaux à 1. 5.-Procédé suivant l'une quelconque des revendications 2 à 4, caractérisé par le fait qu'on choisit (R) proportionnel à (Ri x R2) avec (Ri = 1 3 At) et (R2 = 1/!1/Tdpréc . 6.-Procédé suivant l'une quelconque des revendications 1 à 5, ca- ractérisé par le fait qu'on choisit l'intervalle de temps (At) séparant les mesures du taux de comptage (Tc) en fonction de la valeur du taux de comptage (Tc).