Brevet FR3006076A1 — publié le 20141128

Procédé de caractérisation radiologique d'un générateur de vapeur DOMAINE TECHNIQUE DE L'INVENTION La présente invention se rapporte au domaine des centrales électro- nucléaires, et plus particulièrement à un procédé permettant de caractériser radiologiquement un générateur de vapeur de centrale nucléaire. ARRIÈRE-PLAN TECHNOLOGIQUE DE L'INVENTION La caractérisation radiologique permet d'évaluer l'activité d'une source radioactive. Le besoin de caractérisation radiologique d'un générateur de vapeur apparaît lors de la fin de vie de celui-ci afin de gérer au mieux sa mise en stockage. Diverses solutions sont envisageables pour se débarrasser d'un générateur de vapeur usé. Dans une première solution, le générateur de vapeur peut être mis en décharge radioactive par stockage, sans décontamination, après avoir été éventuellement découpé. Dans ce cas, la mise en stockage nécessite une caractérisation radiologique. Dans une deuxième solution, le générateur de vapeur peut être décontaminé pour être mis en stockage définitif dans une catégorie de déchets radioactifs inférieure telle que les déchets très faiblement actifs (TFA) en France. Dans ce cas, la caractérisation radiologique avant décontamination permet de définir la performance de décontamination à réaliser. Après décontamination, la caractérisation radiologique est requise pour la mise en stockage du générateur de vapeur. Dans une troisième solution, le générateur de vapeur peut être décontaminé pour le recyclage des métaux. Dans ce cas, la caractérisation radiologique avant décontamination permet de définir la performance de la décontamination à réaliser. Après décontamination, la caractérisation radiologique est requise pour vérifier que les critères d'acceptation en ce qui concerne la radioactivité sont bien remplis en filière de recyclage. La caractérisation radiologique peut être utilisée également pour caractériser un générateur de vapeur en phase de maintenance. Le générateur de vapeur considéré est notamment un générateur de vapeur utilisé dans les réacteurs nucléaires de type réacteur à eau pressurisée 5 (REP) ou, en anglais, Pressurized Water Reactor (PWR) ou d'autres générateur de vapeur dont la caractérisation radiologique peut être nécessaire. Un générateur de vapeur comprend habituellement des tubes en U véhiculant un fluide primaire ou eau primaire qui permettent d'échanger de la chaleur avec un fluide secondaire ou eau secondaire. Cet échange de chaleur a 10 pour conséquence la vaporisation du fluide secondaire d'où le nom de générateur de vapeur. L'ensemble des tubes en U forme un faisceau tubulaire comprenant plusieurs milliers de tubes, par exemple 3000 tubes. L'essentiel de l'activité radiologique d'un générateur de vapeur est distribuée sur la surface interne de ces tubes en U. 15 Chaque tube en U part d'une boîte à eau dite boîte à eau chaude ou autrement appelée boîte chaude correspondant à l'entrée de l'eau primaire pour arriver dans une autre boîte à eau dite boîte à eau froide ou autrement appelée boîte froide correspondant à la sortie de l'eau primaire. L'ensemble des deux boîtes est approximativement hémisphérique. 20 Une première extrémité de chaque tube en U débouche dans la boîte d'eau froide et l'autre extrémité débouche dans la boîte d'eau chaude. Les deux boîtes sont séparées par une paroi ou plaque de partition. La partie rectiligne de chaque tube en U issue de la boîte chaude est appelée « demi-tube chaud » et la partie rectiligne de chaque tube en U issue 25 de la boîte froide est appelée « demi-tube froid ». Dans le faisceau tubulaire, l'ensemble des demi-tubes froids est séparé de l'ensemble des demi-tubes chauds par un espace vide appelé « ru ». On pose que les demi-tubes chauds sont similaires entre eux et que les demi-tubes froids sont similaires entre eux. Les extrémités des tubes peuvent être disposées selon une maille 30 carrée. Selon la type de générateur de vapeur d'autres types de mailles peuvent être envisagés telles que des mailles triangulaires. Actuellement, plusieurs techniques permettent la caractérisation radiologique. Une technique consiste à prélever un demi-tube, ou la surface interne d'un demi-tube afin d'en mesurer l'activité. La technique est fiable dans la mesure où le prélèvement est bien représentatif. En revanche, cette technique amène quelques difficultés qui sont, par exemple, les suivantes. Il est nécessaire de prélever un grand nombre de demi-tubes avec éventuellement plusieurs prélèvements dans chaque demi-tube afin de garantir la représentativité globale. Même s'il est possible de prélever de la surface interne des demi-tubes pour chaque boîte d'eau, il est nécessaire d'en prélever toute la surface et rien que la surface en question. En effet, l'activité est distribuée dans les oxydes déposés ou créés à la surface sur quelques micromètres de métal de base, par exemple, entre 5 um et 20 um. Il faut être sûr de prélever toute cette épaisseur sans percer le tube. Comme le résultat à obtenir est exprimé en activité surfacique, il faut avoir une connaissance précise de la surface prélevée pour que le résultat soit fiable. En outre, le prélèvement d'un demi-tube est difficile, voire impossible à cause de l'inaccessibilité du demi-tube à découper au milieu du faisceau tubulaire. Au total, cette technique de prélèvement est destructive, ce qui interdit, par exemple, une décontamination plus avant s'il s'avère que les objectifs ne sont pas atteints. Cette technique est aussi onéreuse et délicate à mettre en oeuvre. De plus, elle permet de ne mesurer que l'activité à l'intérieur des tubes alors qu'il peut exister une activité du côté du tube d'eau secondaire, distribuée sous forme de points chauds. Une autre technique se base sur la mesure de rayonnement gamma, ou autrement appelé comptage gamma, soit global, soit spectrométrique. Dans le cas des mesures spectrométriques, les mesures sont réalisées par la détection de certaines énergies d'émissions caractéristiques. Par exemple, il est possible de détecter les raies 1173 keV et 1332 keV du 60Co ou la raie 662 keV du 137Cs. Le choix d'une mesure du rayonnement gamma global ou spectroscopique dépend du temps de refroidissement du générateur de vapeur, qui correspond de la durée depuis l'arrêt du réacteur ou depuis le découplage du générateur de vapeur du circuit primaire, et des performances attendues de la caractérisation radiologique. Les principes généraux de cette technique sont les suivants : - mesure du rayonnement gamma depuis l'extérieur du générateur de vapeur et/ou dans les tubes et/ou ailleurs dans le générateur de vapeur. Pour l'intérieur des tubes, on utilise une sonde suffisamment petite pour être introduite dans les tubes qui ont un diamètre, par exemple, de 15 mm. Il faut noter qu'avec ces petites dimensions, les performances d'une mesure spectrométrique sont faibles et conduisent à des temps d'acquisition importants. La mesure de rayonnement gamma globale est mieux adaptée. - exploitation des mesures gamma à l'aide de modèles numériques et/ou physiques permettant de prendre en compte les propagations, atténuations et diffusions des rayonnements entre leur lieu d'émission et lieu de détection. - quelques analyses par prélèvement de manière à déterminer le spectre de radioéléments avec la proportion d'activité de chaque élément. À partir de mesures purement gamma, il est possible de déduire l'activité de radioéléments peu ou faiblement émetteurs de rayonnement gamma tels que les émetteurs bêta purs ou le plutonium. Il est aussi possible de distribuer une activité gamma globale mesurée en activité par émetteurs de rayonnement gamma. Plusieurs avantages découlent de la mesure de rayonnement gamma : - La mesure est non destructive. - La mesure est non locale. En effet, le rayonnement gamma se propage et traverse la matière, ce qui fait que, par exemple, un sondage dans un tube émetteur de rayonnement gamma est également représentatif des tubes voisins. Le nombre de sondage nécessaires est donc moindre que celui des prélèvements nécessaires dans une technique de prélèvements. - La mesure est rapide et peu onéreuse. La difficulté de la technique provient de la combinaison de deux caractéristiques. La première caractéristique est le fait que le système de générateur de vapeur est extrêmement complexe par son milieu inhomogène comprenant des milliers de tubes. La deuxième caractéristique est le fait que la distribution d'activité dans les tubes, parfois hors des tubes, est inhomogène. Il peut en effet exister plusieurs types de tubes et de demi-tubes. Il y a des différences entre les demi-tubes chauds et les demi-tubes froids. Il y a le cas aussi de tubes particuliers tels que les tubes pastillés. Les tubes pastillés sont des tubes qui ont été neutralisés par des pastilles car ils présentaient des défauts tels que des fuites. DESCRIPTION GÉNÉRALE DE L'INVENTION La présente invention a pour premier objet de permettre la caractérisation radiologique d'un générateur de vapeur à tube de circulation de fluide, quelque soit le type de tube, par mesure de rayonnement gamma. La présente invention a également pour objet de permettre la caractérisation radiologique d'un générateur de vapeur par mesure de rayonnement gamma en prenant en compte des distributions inhomogènes de l'activité causées par exemple par un certain nombre de tubes particuliers et un certain nombre de zones différenciées avec pour chaque zone une homogénéité de l'activité des tubes non particuliers et avec des différences d'activité entre chaque zone. À cet effet, l'invention concerne un procédé de caractérisation radiologique d'un générateur de vapeur comprenant une pluralité de tubes en U, chaque tube en U étant formé d'au moins deux demi-tubes dont les coordonnées d'au moins un demi-tube sont enregistrées dans un fichier de mesure stocké dans un moyen de stockage numérique, le générateur de vapeur comprenant des tubes dits non-pastillés et éventuellement des tubes pastillés, le procédé étant caractérisé en ce qu'il comprend au moins les étapes suivantes : - détermination de familles de demi-tubes parmi les demi-tubes ; - au moins une mesure par une sonde du taux de comptage Ni du rayonnement gamma dans au moins un demi-tube par famille, les mesures étant enregistrées dans le fichier de mesures en association avec les coordonnées du ou des demi-tubes sondés ; - détermination pour chaque famille par le module de calcul du ou des coefficients eij d'influence d'une matrice d'influence pour obtenir une matrice ayant une dimension inférieure ou égale au nombre de familles, un coefficient d'influence correspondant à l'influence du taux de comptage d'un terme source de rayonnement gamma sur le taux de comptage du demi-tube sondé, les coefficients d'influence étant enregistrés dans un fichier de calcul stocké dans le moyen de stockage numérique ; - détermination par le module de calcul de l'inverse de la matrice d'influence dont les coefficients sont enregistrés dans le fichier de calcul stocké dans le moyen de stockage numérique ; - détermination par le module de calcul d'une valeur d'activité Ai du générateur de vapeur à partir des mesures de taux de comptage, la valeur d'activité étant enregistrée dans un fichier de résultats stocké dans le module de stockage numérique. Selon une autre particularité, l'étape de détermination de l'inverse de la matrice d'influence par le module de calcul comprend une étape dans laquelle le module de calcul calcule l'inverse la matrice d'influence en considérant que l'activité des demi-tubes dans une famille est la même que celle mesurée pour au moins un demi-tube de cette famille. Selon une autre particularité, l'étape de détermination de l'inverse de la matrice d'influence par le module de calcul comprend une étape dans laquelle le module de calcul calcule l'inverse de la matrice d'influence en considérant que l'activité des demi-tubes varie en fonction des coordonnées des demi-tubes. Selon une autre particularité, le ou les coefficients d'influences sont déterminés à partir d'une fonction d'interpolation eij = f(diJ) exécutée par le module de calcul en utilisant les valeurs des mesures de taux de comptage enregistrées dans le fichier de mesures et au moins des informations représentatives des distances entre le tube sondé et les termes sources enregistrées dans un fichier représentatif du générateur de vapeur mémorisé dans le moyen de stockage numérique. Selon une autre particularité, le ou les coefficients d'influence sont déterminés en outre à partir d'une fonction d'interpolation eij = f(dii,autres_paramètres) exécutée par le module de calcul en utilisant les valeurs des mesures de taux de comptage enregistrés dans le fichier de mesures et au moins un autre paramètre déterminé à partir de la moyenne d'un échantillonnage d'au moins une caractéristique structurelle du générateur de vapeur. Selon une autre particularité, la caractéristique structurelle du générateur de vapeur est une information représentative de l'épaisseur moyenne EM de matière traversée selon un plan perpendiculaire au tube sondé entre le terme source et le centre du demi-tube sondé enregistré dans un fichier représentatif du générateur de vapeur mémorisé dans le moyen de stockage numérique. Selon une autre particularité, le terme source est la partie interne d'un demi-tube et/ou une source autre que la partie interne d'un demi-tube. Selon une autre particularité, l'étape de détermination de familles comprend : - la constitution d'au moins une famille de demi-tubes pastillés comprenant chacune un seul demi-tube pastillé, et - la constitution d'au moins une famille de tubes non pastillés comprenant tous les demi-tubes non pastillés similaires en activité radioactive. Selon une autre particularité, l'étape de détermination de familles 5 comprend : - la constitution d'au moins une famille de demi-tubes pastillés comprenant chacune un seul demi-tube pastillé, et - la constitution d'au moins une famille de demi-tubes non pastillés de telle manière que chaque famille de demi-tubes non pastillés comprend 10 des demi-tubes non pastillés adjacents et ayant au moins un critère de classification semblable. Selon une autre particularité, un moyen d'affichage du module de calcul affiche chaque demi-tube mesuré du générateur de vapeur en association de ses coordonnées, la constitution des familles étant effectuée par l'utilisateur qui 15 utilise des moyens interactifs du module de calcul pour sélectionner sur le moyen d'affichage, pour chaque famille, les demi-tubes à affecter à cette famille à l'aide d'un moyen de sélection. Selon une autre particularité, lors de l'étape de détermination par le module de calcul de l'inverse de la matrice d'influence comprend une étape de 20 factorisation de type Cholesky de la matrice d'influence par le module de calcul. L'invention, avec ses caractéristiques et avantages, ressortira plus clairement à la lecture de la description faite en référence aux dessins annexés dans lesquels : La figure 1 est un graphe représentant l'influence relative d'un demi- s tube en fonction de la distance entre le demi-tube influent et le demi-tube sondé. La figure 2 est un graphe représentant la contribution cumulée relative en fonction de la distance entre le demi-tube influent et le demi-tube sondé. La figure 3 représente une coupe transversale du générateur de vapeur 10 avec la représentation de parcours d'un rayonnement gamma selon la direction de ce rayonnement. La figure 4 représente la coupe transversale de trois demi-tubes dont un demi-tube sondé, un demi-tube atténuant et un demi-tube influent. 15 vapeur. La figure 5 représente une coupe longitudinale d'un générateur de vapeur. La figure 6 représente une coupe transversale du générateur de DESCRIPTION DES MODES DE RÉALISATION PRÉFÉRÉS DE L'INVENTION 20 La suite de la description fera référence aux figures citées ci-dessus. La caractérisation radiologique d'un générateur (0) de vapeur est un problème difficile à résoudre. Le système permettant de mettre en oeuvre le procédé comprend au moins un module de calcul tel qu'un processeur d'un dispositif informatique. Le 25 système comprend également un moyen de stockage numérique connecté au module de calcul permettant de mémoriser et d'enregistrer des fichiers. Le système peut comprendre également un moyen d'affichage connecté au module de calcul tel qu'un écran et/ou un moyen de saisie connecté au module de calcul tel qu'un clavier ou une souris. Une sonde (12a) est connectée à un moyen de traitement permettant d'interpréter le signal obtenu par les mesures pour l'afficher sur un moyen d'affichage. La valeur affichée sur le moyen d'affichage est alors saisie par le moyen de saisie du dispositif informatique afin de la stocker dans le moyen de stockage numérique. Comme nous l'avons exprimé ci-dessus, plusieurs méthodes existent déjà. La construction d'un modèle réaliste, physique ou numérique, d'un 10 générateur (0) de vapeur apparaît comme impossible. Pour un générateur (0) de vapeur de 3000 tubes, le problème à résoudre peut se décrire sous forme matricielle : /No\ /e0,0 e0,1 --- eon\ //41\ [1] N1 1 e1,0 e1,1 el,n A2 _ - . . - - el - - . . 1 \1\ni \en,0 en,1 --- en,ni \Ani où : 15 i prend les valeurs de 0 à n, n est le nombre total de demi-tubes (11, 12) moins un demi-tube, Par exemple, i = 0 à n = 5999 pour un générateur (0) de vapeur à 3000 tubes. Ai est l'activité dans le demi-tube n°i exprimée, par exemple, en 20 Bq/cm2, Ni est le taux de comptage de rayonnement gamma mesuré en sondant le demi-tube n°i, eii est le coefficient d'influence de l'activité dans le demi-tube n°j sur le résultat du sondage gamma dans le demi-tube n°i. Les tubes du générateur (0) de vapeur peuvent avoir, par exemple et de façon non limitative, les dimensions suivantes : un diamètre interne de 15 mm environ et un diamètre externe de 20 mm environ. Dans la présente demande, la longueur des tubes est supposée comme infinie, mais en réalité un demi-tube a une longueur comprise, par exemple et de façon non limitative, entre 10 m et 15 m. Le générateur (0) de vapeur a un ru (3) de largeur, par exemple et de façon non limitative, de 95 mm. La matière des tubes peut être, par exemple, de l'acier inoxydable 316 10 possédant une densité de 8000 kg/m3. Chaque tube en U part d'une boîte à eau dite boîte à eau chaude ou autrement appelée boîte (13c) chaude correspondant à l'entrée de l'eau primaire pour arriver dans une autre boîte à eau dite boîte à eau froide ou autrement appelée boîte (13f) froide correspondant à la sortie de l'eau primaire. 15 La partie rectiligne de chaque tube en U issue de la boîte (13c) chaude est appelée demi-tube (11c) chaud et la partie rectiligne de chaque tube en U issue de la boîte (13f) froide est appelée demi-tube (11f) froid. Le faisceau (1) formé par l'ensemble des demi-tubes (11f) froids est le faisceau (1) froid. Le faisceau (2) formé par l'ensemble des demi-tubes (11c) chauds est le faisceau (1) 20 chaud. Le ru (3) sépare la zone où débouchent les demi-tubes (11f) froids de la zone où débouchent les demi-tubes chauds (11c). Il peut exister des termes sources d'activité supplémentaires ailleurs que sur la surface (12b) interne des tubes. Les termes sources supplémentaires se traduisent par autant de dimensions supplémentaires à la matrice, avec, par 25 exemple, le demi-tube n°n représentant, en fait, l'activité distribuée du côté du circuit (14) d'eau secondaire pour un lieu de sondage gamma associé. Dans la suite de la description par souci de simplification, nous ne parlerons que de la partie interne (12b) des demi-tubes. Mais, il faut comprendre que la description est valable à la fois pour la surface interne des demi-tubes et pour les termes sources supplémentaires ailleurs que sur la surface interne des tubes. Pour un générateur (0) de vapeur à 3000 tubes, la matrice d'efficacité des demi-tubes correspondant à la matrice des coefficients d'influence est de 5 dimension 6000 x 6000 et possède donc 36 000 000 coefficients. Ce nombre n(n+1) 5999(5999+1) peut être réduire à 2 2 n + 1 = 5999 = 17 997 001 grâce aux symétries eii = ej,i et ei,i = e0,0. Ces symétries ne s'appliquent pas aux termes sources supplémentaires qui se trouvent hors des surfaces internes des tubes. 10 Les valeurs Ai sont les inconnues et les valeurs N1 sont des données obtenues à partir des résultats de mesure par sondage. Le problème consiste donc à inverser une matrice très grande dont il faut d'abord calculer les coefficients. Le calcul d'un seul coefficient ou de dix coefficients ou de cent 15 coefficients est tout à fait réalisable avec des codes de calculs disponibles tels que le code de transport Monte-Carlo à N particules (Monte-Carlo N-Particle transport code en anglais ou MCNP) ou tout autre code permettant à un module de calcul de réaliser ce calcul. Il est possible d'envisager d'automatiser le calcul de coefficients plus 20 nombreux mais pour 17 997 001 coefficients, le temps de calcul machine est très important. En effet, il faut en moyenne 1 h de calcul pour le calcul d'un seul coefficient avec une méthode de type Monte-Carlo qui est la seule méthode à même de calculer les coefficients relatifs au comptage gamma global. Ainsi, pour 17 997 001 coefficients, il faudrait plus de 2000 ans de calcul. Même avec 25 les méthodes dites « en ligne directe » qui ne calculent que des coefficients relatifs à des comptages spectroscopiques et qui sont capables de réaliser le calcul d'un coefficient en 1 min, il faudrait 34 ans de temps de calcul machine. En considérant que l'influence de l'activité d'un demi-tube sur le sondage dans un autre demi-tube diminue avec la distance, on peut penser qu'à partir d'une certaine distance cette influence devient négligeable. Mais le nombre de demi-tubes présents à une certaine distance d'un demi-tube donné augmente avec la distance, ce qui contrebalance l'effet de diminution. Cette diminution n'est pas totale mais elle est suffisante pour que la distance à partir de laquelle il est justifié de négliger l'influence soit comparable au rayon du générateur (0) de vapeur. La figure 1 représente l'ordre de grandeur de l'influence d'un demi-tube par rapport à la distance partant de ce demi-tube d'origine. La figure 2 représente l'accumulation cumulée relative par rapport à la distance partant de ce demi-tube d'origine. Ces figures montrent plusieurs choses. D'une part, l'influence d'un demi-tube pour le sondage sur un autre demi-tube décroît rapidement avec la distance entre les demi-tubes. Par exemple, il ne reste que 1/100 de l'influence du demi-tube sondé sur lui-même dès une distance de 20 cm. D'autre part, si l'on se permet 5% de sous-estimation, il faut calculer les coefficients jusqu'à au moins 40 cm dans le cas de demi-tubes ayant tous la même activité. Nous allons maintenant étudier les différents principes et méthodes matricielles de résolution des systèmes linéaires. Deux groupes de méthodes sont généralement utilisés pour résoudre des systèmes linéaires de grande taille : des méthodes directes et des méthodes itératives. L'utilisation de l'un ou l'autre de ces groupes dépend généralement de la structure du système à résoudre qui peut être complet ou incomplet, mais aussi du problème à étudier. On s'accorde sur le fait qu'il n'y a pas de solution évidente sur le choix d'une approche ou d'une autre si l'on considère en même temps les deux propositions que sont la structure et le problème à résoudre. Ici, de façon générale et très synthétique, nous présentons les contraintes à prendre en compte pour privilégier une orientation d'utilisation de méthode de résolution. De façon particulière, nous présentons les méthodes numériques les plus utilisées à l'heure actuelle. Ainsi, avant d'aborder la résolution matricielle nécessaire à la solution de notre problème, nous aborderons de façon générale les méthodes d'analyse numérique existantes et les contraintes associées. De façon générale, il ressort que les méthodes directes sont utilisées pour résoudre des problèmes à une ou deux dimensions alors que les méthodes itératives sont privilégiées pour résoudre des problèmes à trois dimensions. S'ajoute également dans le choix de l'utilisation d'une méthode ou d'une autre, un coût de résolution qui pourrait être exprimé en tant que paramètre ou fonction. Ce coût peut s'apparenter au temps nécessaire à l'obtention des solutions. Et comme tout paramètre ou fonction de coût, il s'agit de la minimiser. La minimisation passe classiquement par deux approches complémentaires : d'un côté, la réduction de la complexité du système à résoudre par l'analyse numérique qui vise à poser le problème pour réduire le nombre d'opérations élémentaires à effectuer, de l'autre côté, l'augmentation des performances de traitement. On parle ici, d'un côté, d'algorithme ou de systématisation de l'enchaînement des calculs allié de l'autre côté à une architecture adaptée, c'est-à-dire, les fonctions de traitement en séries ou répartis associées aux moyens informatiques. De façon particulière, concernant l'analyse numérique, comme il a été noté précédemment, si l'ordre d'une matrice dépasse quelques unités, résoudre un système linéaire dont la représentation matricielle est N=E-A (voir équation [I]) devient délicat, car le nombre d'opérations à effectuer est trop important. Cela est le cas si on pratique le calcul E-1 puis A = E-I- - N avec une méthode comme le déterminant de Cramer, poussé par une application « littérale » de la définition de l'inverse d'une matrice : E - E-I- = E-1.E = I tcom(E) obtenue par la relation E-1 = 1E1 . Avec I la matrice identité, tCOM(E) la transposée de la comatrice E etlElle déterminant de E. En fait, les méthodes directes de résolution de système linéaires visent à déterminer une matrice M inversible telle queM-N=M-E-A. Il s'agit en fait de se ramener à la résolution de systèmes triangulaires qui présentent des propriétés intéressantes. Ces propriétés sont exploitées dans les méthodes de Gauss, de factorisation dites LU (LowerUpper). En fonction de la structure de la matrice, la factorisation LU s'avère encore plus intéressante par un nombre d'opérations moindre par des matrices diagonales, des matrices bandes, des matrices tridiagonales, etc. Ainsi, d'autre factorisations peuvent être envisagées : factorisation QR, Householder, factorisation de Cholesky, etc. Les méthodes directes visent à donner des solutions exactes et de manière générale, il s'agit de diminuer la complexité de la résolution du système linéaire. À titre d'exemple, on présente ci-dessous une évaluation en fonction des différentes méthodes : - la méthode du déterminant de Cramer aboutit à une complexité en 0(n!), - la méthode de Gauss-Jordan aboutit à une complexité en 0(2n3), - la méthode de Householder ou QR aboutit à une complexité en 0((4n3)/3), - la méthode de décomposition de Cholesky aboutit à une complexité en 0(n3/3). Les méthodes directes reposent sur un nombre déterminé et fini d'opérations. Ce nombre est d'autant plus élevé que la taille du système est grande. La propagation d'éventuelles erreurs de calculs va de pair avec ce nombre. De toute façon, il est difficile de tenir compte et surtout de quantifier des erreurs avec ces méthodes directes. Aussi, la solution n'est que rarement exacte. Quant aux méthodes itératives, elles visent à donner des solutions approchées avec la possibilité d'estimer les erreurs. Le principe de ces méthodes itératives est toujours le même. Dans un système d'équations de type b = A. x, on décompose A = M - N de telle façon que M soit inversible et que l'une des deux conditions soit remplie. Si A est symétrique définie positive, M + tN est aussi définie positive et le rayon spectrale qui correspond au maximum des valeurs de M - 1. N est petit, typiquement inférieur à 1. Dès lors, les méthodes de choix de M et N sont alors connues : des méthodes stationnaires (les méthodes de Jacobi, Gauss-Seidel, sur-relaxation successive (AOR), etc.), des méthodes non stationnaires (les méthodes de gradient conjugué, Quasi Minimale Résiduelle de Freund et Nachtigal, de Krylov, etc.) Toutes ces méthodes sont stables, mais les conditions de convergences sont souvent difficiles à mettre en oeuvre. On estime qu'une méthode itérative ne sera compétitive que si elle converge avec un nombre d'itérations indépendants de la taille du système à résoudre ou bien croissant de manière sous-linéaire. Dans le cadre de l'état de l'art, la mesure gamma n'est possible que moyennant des hypothèses extrêmement réductrices, comme, par exemple, une activité homogène sur tous les tubes et un faisceau tubulaire de générateur (0) de vapeur modélisé par un milieu homogène. Aussi, il est à privilégier l'usage d'une méthode directe qu'il est nécessaire d'adapter. Le procédé de caractérisation radiologique d'un générateur (0) de vapeur selon l'invention comprend une étape de détermination de famille de demi-tubes parmi l'ensemble des demi-tubes. Le procédé comprend également une étape de sondage. Le sondage est réalisé par une sonde (12a). Cette sonde est, par exemple et de façon non limitative, une sonde CZT qui est une diode de composition chimique Cc10,9Zno,1Te. Le volume sensible de cette diode est par exemple de façon non limitative, un parallélépipède de dimension d'environs 5 mm x 5 mm x 2,5 mm et d'une densité de 5780 kg/m3. La sonde est introduite dans un tube (11s) de telle manière que l'une de ses grandes dimensions, par exemple 5 mm, soit dans l'axe du tube. La sonde peut être utilisée en comptage globale, c'est-à-dire qu'elle détecte les événements à partir d'un seuil, par exemple de 10 keV, et compte une impulsion pour tout événement ayant déposé une énergie supérieure à ce seuil dans le volume sensible. De ce fait, la sonde est sensible non seulement au rayonnement direct n'ayant pas interagit entre son émission et sa détection, mais aussi au rayonnement diffusé. Elle peut être aussi utilisée en comptage spectrométrique. La sonde est utilisée, de préférence, en comptage globale. En effet, si elle était utilisée en comptage spectrométrique, son comptage ne serait pas influencé par le rayonnement diffusé. Cela simplifierait les calculs mis en oeuvre par le module de calcul mais le petit volume de la sonde, imposé par le diamètre des tubes, rend le comptage spectrométrique trop peu efficace pour une mesure avec des temps d'acquisition raisonnables. L'activité à déterminer est déposée sur la surface interne de chaque tube. Le dépôt est supposé uniforme dans chaque demi-tube. L'activité est constituée, par exemple, de 60Co. Les résultats des mesures sont exprimés en impulsions détectées par désintégration par seconde dans un demi-tube, par exemple, de 5 m. Une désintégration de 60Co produisant deux photons à 1173 keV et 1332 keV, le résultat est donc la somme des efficacités en impulsion par photon émis pour les deux niveaux d'énergie. L'étape de sondage ne revient pas à sonder tous les demi-tubes du générateur mais seulement un échantillonnage de tube. Ainsi, par exemple pour un générateur comprenant 6000 demi-tubes, la matrice de dimension 6000 x 6000 est réduite à une dimension plus petite. Cette dimension est inférieure ou égale au nombre de sondages sur la base de familles de demi-tubes dont les activités sont corrélées ou équivalentes. Tout d'abord, il est nécessaire de définir qu'un tube pastillé est un tube qui n'est plus utilisé dans le générateur (0) de vapeur, c'est-à-dire qu'il a été bouché afin que l'eau primaire ne circule à l'intérieur. Généralement, un tube est pastillé par une pastille suite à un endommagement du tube, tel qu'une fuite. Un tube non pastillé est un tube qui est utilisé. Pour illustrer la notion de famille, considérons un cas de base, avec un générateur (0) de vapeur contenant quelques tubes pastillés, des demi-tubes (11) non pastillés chauds ayant une homogénéité d'activité et des tubes (12) non pastillés froids ayant une homogénéité d'activité avec des différences d'homogénéité d'activité entre les demi-tubes chauds et les demi-tubes froids. Les demi-tubes froids sont donc considérés comme similaires en activité radioactive entre eux et les demi-tubes chauds sont considérés comme similaires en activité radioactive entre eux. Ici, « similaire » signifie soit égal, soit compris dans un intervalle qui est fonction de l'incertitude recherchée pour la détermination d'activité du générateur de vapeur, soit corrélé. Par activité radioactive corrélée entre eux veut dire que l'on peut établir une loi permettant de déduire l'activité d'un demi-tube à partir de la position de ce demi-tube et de l'activité d'un autre demi-tube de la famille, la position de cet autre demi-tube étant connue. Nous avons alors, une famille pour les tubes chauds, une famille pour les tubes froids et une famille pour chaque tube pastillé. Par exemple, s'il y a 10 tubes pastillés, il y a donc 20 demi-tubes pastillés. Il y a alors 22 familles, l'activité des demi-tubes étant la même entre les membres d'une même famille : - chacun des demi-tubes pastillés constitue à lui seul une famille. Donc avec 10 tubes pastillés, il y a 20 familles de demi-tubes pastillés, - l'ensemble des demi-tubes chauds non pastillés constitue une famille, - l'ensemble des demi-tubes froids non pastillés constitue une famille. Notons que l'hypothèse d'homogénéité d'activité à l'intérieur d'une famille est la plus simple. Mais, il est possible de travailler sur des hypothèses plus complexes où, par exemple, l'activité varie continûment avec la position géographique du demi-tube. Il est possible également que deux tubes ayant des demi-tubes froids faisant partie d'une même famille aient des demi-tubes chauds de faisant pas partie d'une même famille, c'est-à-dire qu'il n'y a pas de symétrie de familles entre les demi-tubes froids et les demi-tubes chauds. Dans l'hypothèse où l'activité varie continûment avec la position géographique du demi-tube, au moins une famille de demi-tubes non pastillés est constituée de telle manière que chaque famille de demi-tubes non pastillés comprend des demi-tubes non pastillés adjacents et ayant un critère de classification semblable. Un critère de classification semblable peut être, de façon non limitative, un taux de comptage similaire. Ce critère de classification peut être une limite de distance à ne pas dépasser par rapport à une source de rayonnement gamma. D'autres critères de classification peuvent être utilisés seuls ou en combinaison avec d'autres critères. Les sondages sont réalisés par au moins une sonde. La sonde mesure le taux de comptage à l'endroit où elle est disposée et transforme cette mesure en un signal électrique. Ce signal électrique est interprété par le module de traitement en une valeur représentative de ce signal. La valeur représentative de ce signal est alors ajoutée dans un fichier de mesures enregistré dans le moyen de stockage numérique. Le fichier de mesure peut associer également chaque mesure avec les coordonnées de la sonde au moment de la mesure. Ces coordonnées peuvent être le numéro du demi-tube où est réalisée la mesure auquel sont associées les coordonnées cartésiennes du tube. Le numéro du tube peut être saisi par des moyens de saisie par un utilisateur à chaque mesure. Le fichier de mesure peut associer chaque mesure effectuée en fonction de la position du demi-tube sondé définies sur une représentation graphique des demi-tubes d'un générateur (0) de vapeur enregistrée dans le moyen de stockage numérique du module de calcul ou définie par le moyen de saisie par l'utilisateur. Les demi-tubes peuvent être disposés selon une maille carrée ou d'autres types de mailles telles que des mailles triangulaires. Chaque demi-tube est repéré par des coordonnées cartésiennes, elles-mêmes exprimées en nombre de mailles. Au moins un sondage doit être réalisé dans chaque famille. En reprenant l'exemple précédent avec 22 familles, la matrice initiale de 6000 x 6000 se réduit alors à une nouvelle matrice de 22 x 22. Cette matrice est plus facile à inverser mais le problème de calcul des coefficients, qui sont de 17 997 002 dans cet exemple, n'est pas résolu parce que les coefficients de la nouvelle matrice, qui sont de 484 dans cet exemple, sont chacun une somme d'une partie des coefficients de la matrice initiale. Il faut donc toujours calculer tous les coefficients de la matrice initiale avant de pouvoir calculer les coefficients de la nouvelle matrice. Calcul des coefficients Les coefficients sont calculés par le module de calcul du dispositif informatique. Le principe du calcul des coefficients réside dans l'établissement d'une fonction d'interpolation permettant d'exprimer le coefficient d'influence de l'activité du demi-tube j sur le sondage dans le demi-tube i en fonction, dans un premier temps de la distance entre ces deux demi-tubes : eii = f (dl, ), où f( ) est la fonction analytique à établir et du est la distance entre les deux demi-tubes i et j. Pour cela, le module de calcul détermine cette fonction d'interpolation en calculant un échantillonnage de coefficients à partir des mesures de taux de comptage enregistrées dans le fichier de mesures. Le module de calcul corrèle ensuite les valeurs des mesures de taux de comptage à celles des informations représentatives des distances. La distance est calculée par le module de calcul à partir des coordonnées associées dans le fichier de mesures à la mesure en question ou des coordonnées stockées dans le fichier représentatif du générateur (0) de vapeur décrit plus loin. L'avantage de cette fonction est d'envisager l'automatisation du calcul complet des toutes les valeurs de coefficients de la matrice initiale, puis celui des coefficients de la matrice réduite avec un temps de calcul machine très court, de l'ordre de la seconde. Le calcul par le module de calcul des coefficients échantillonnés est plus long mais acceptable lorsque l'échantillonnage est raisonnable. Il n'est effectué qu'une seule fois par type de générateur (0) de vapeur. L'inconvénient est le caractère approximatif de cette fonction. En effet, deux caractéristiques rendent délicate l'interpolation. La première caractéristique provient de l'existence de ru (3) entre le demi-faisceau de demi-tubes chauds et le demi-faisceau de demi-tubes froids. Pour une même distance entre demi-tubes, l'épaisseur de matière traversée est différente selon que l'on traverse ou non le ru (3). Or, l'atténuation du rayonnement gamme est en relation exponentielle avec l'épaisseur de matière traversée, ce qui a des conséquences directes sur les coefficients d'influence. La deuxième caractéristique provient de la périodicité spatiale de la disposition des demi-tubes dans le plan perpendiculaire aux demi-tubes, ce qui a pour effet de générer des directions favorisées et défavorisées. La figure 3 illustre ce phénomène. Plusieurs essais ont montré que interpolation par régression linéaire est, préférentiellement mais de manière non limitative, l'interpolation la mieux 20 adaptée. Dans cette figure 3, la point A correspond au demi-tube sondé. Les points B et C correspondent aux deux demi-tubes influents sur la mesure du demi-tube A. La direction allant de A à B est une direction défavorisée car elle 25 traverse six demi-tubes lors du trajet du demi-tube influent B vers le demi-tube sondé A. La direction allant de A vers C est une direction favorisée car elle traverse quatre demi-tubes lors du trajet du demi-tube influent C vers le dei-tube sondé A alors que la distance entre les demi-tubes A et B est voisine de la distance entre les demi-tubes A et C. Le module de calcul peut, par un échantillonnage adapté des coefficients à calculer, faire la moyenne des effets de ces phénomènes. Ainsi, le ou les coefficients d'influence sont déterminés en outre à partir d'une fonction d'interpolation exécutée par le module de calcul en utilisant les valeurs des mesures de taux de comptage enregistrés dans le fichier de mesures et au moins un autre paramètre déterminé à partir de la moyenne d'un échantillonnage d'au moins une caractéristique structurelle du générateur (0) de vapeur.de la moyenne d'un échantillonnage d'au moins une caractéristique structurelle du générateur (0) de vapeur. La fonction f peut être alors dépendante d'autre paramètres : eii = f (du, autres_paramètres). Les caractéristiques structurelles peuvent être le fait que le rayonnement traverse le ru (3) ou non, l'angle du vecteur joignant le demi-tube sondé et le demi-tube influent par rapport à un repère lié au générateur (0) de vapeur, le nombre de demi-tubes traversés lors du trajet, etc. Préférentiellement, la caractéristique structurelle est une information représentative de l'épaisseur moyenne (EM) de matière traversée selon un plan perpendiculaire au tube (11s) sondé entre la surface (12b) interne du demi-tube (11i) influent et le centre (12a) du demi-tube (11s) sondé dans une géométrie à deux dimensions. Le centre (12a) du demi-tube sondé est le lieu où se situe la sonde comme représenté par la figure 4. Il est possible également que l'information représentative de l'épaisseur moyenne de matière traversée puisse être déterminée dans une géométrie à trois dimensions. Dans cette figure 4, le rayonnement (13) mesuré par la sonde (12a) et provenant du demi-tube (11i) influent est atténué par un demi-tube (11a) intermédiaire. Une seconde fonction g est alors déterminée par le module de calcul pour calculer les coefficients d'atténuation : eii = f(dii) . g(Em). Le moyen de stockage numérique comprend un fichier représentatif du générateur (0) de vapeur sondé. Ce fichier comprend des valeurs correspondant aux coordonnées de chaque demi-tube et/ou le maillage des demi-tubes dans le plan de coupe perpendiculaire aux tubes, la forme des mailles et leurs dimensions. Le fichier représentatif du générateur (0) de vapeur peut comprendre aussi des valeurs correspondant aux dimensions de chaque demi-tube telles que leur section interne et externe, leur épaisseur et leur taille ainsi qu'à la dimension du ru (3). Le fichier représentatif du générateur (0) de vapeur peut comprendre aussi des valeurs correspondant à la nature du matériau des tubes, à la densité du matériau des demi-tubes, à la position respective des demi-tubes sondés et des demi-tubes influents. À partir du fichier représentatif du générateur (0) de vapeur, le module de calcul calcule une intégrale à une dimension sur le pourtour interne du demi- tube influent dans le plan de coupe perpendiculaire aux tubes et passant par la sonde. Le module de calcul calcule alors la somme des épaisseurs traversées dans l'épaisseur du demi-tube influent, dans l'épaisseur de tous les demi-tubes traversés et dans l'épaisseur du demi-tube sondé comme représenté sur la figure 4. Le module de calcul détermine ainsi l'épaisseur moyenne à partie du fichier représentatif du générateur (0) de vapeur par les valeurs de forme des mailles, le diamètre interne des tubes, le diamètre externe des tubes, la nature du matériau des tubes, la densité du matériau des tubes, la position respective des demi-tubes sondés et des demi-tubes influents. Les valeurs d'épaisseurs sont alors enregistrées dans le fichier de mesures en étant associées au demi- tube sondé en question. Ainsi, le ou les coefficients d'influence sont déterminés à partir d'une fonction d'interpolation établie par le moyen de calcul à partir des mesures de taux de comptage enregistrés dans le fichier de mesures et d'au moins l'épaisseur moyenne de matière traversée selon un plan perpendiculaire au tube sondé entre le terme source qui peut être la surface interne du demi-tube influent et le centre du demi-tube sondé et/ou une source autre que la surface interne d'un demi-tube. Ainsi, comme pour la fonction f, à partir d'un échantillonnage des valeurs de mesures et les valeurs d'épaisseurs calculées enregistrées dans le fichier de mesures, il est possible d'interpoler la fonction g. Plusieurs essais ont montré que préférentiellement mais de manière non limitative, une interpolation polynomiale l'interpolation la mieux adaptée. Ce calcul de coefficients colle au plus près aux phénomènes physiques qui sont des effets de distance pure et d'atténuation lors de la traversée de matière, tout en restant simple et rapide à calculer. Ce calcul prend en compte également la traversée ou non du ru (3). Les coefficients de la matrice initiale calculés par le module de calcul sont stockés dans un fichier mémorisé dans le moyen de stockage numérique. Résolution du système Après le calcul des coefficients par le module de calcul, la matrice des coefficients est inversée. La méthodologie utilisée permet une résolution directe des différentes 20 équations correspondant aux taux de comptage, constituant le système à résoudre par le module de calcul. Par convention, un tube dépastillé est nommé « P » et un tube non pastillé est nommé « NP » ; la boîte froide est nommée « BF » et la boîte 25 chaude est nommée « BC ». Par convention, on attribue l'indice 0 au demi-tube sondé. Les demi-tubes dépastillés sont des demi-tubes pastilles qui ont été dépastillés afin de pouvoir introduire la sonde à l'intérieur et faire des mesures. Les mesures par la sonde permettent d'observer un gradient des valeurs d'activités sur les tubes NP. En complément de la notion de famille pour les tubes NP, lorsqu'un tube P considéré n'est pas entouré uniquement de tubes NP dans sa zone d'influence, l'activité propre de chaque tube P situé dans cette zone doit être prise en compte dans l'équation qui peut alors s'écrire de la manière suivante : .ek + --- + A1.3, .ek k J -2,1-2,1 m = eoo + A r2 i2 + Al,, 1_,1 + A 1, \\ fri;. ek1,NFI1,NF k1,NF 11,NF Nrmim - Armim.e00 + + --- + A1.3 + AiP ek + ANP NF e kM,NFIM,NF kM,NF .M,NF avec M : le nombre total de demi-tubes P en BF et en BC. N im F i- m. : . le comptage mesuré par la sonde en coups par seconde du mième tube P considéré, situé aux coordonnées(im, jm) en BF ou en BC. ekrn,n/rn,n : l'efficacité calculée par le module de calcul du mième tube P en BF ou en BC situé en (im, Jin) par rapport au nième tube P considéré, situé 5 en (in, jn) en BF ou en BC. e00 : l'efficacité propre du tube sondé calculée par le module de calcul. NF : le nombre de familles en BF et en BC, c'est-à-dire NF est égal en principe à deux fois le nombre de familles dans chaque boîte. ANiP : l'activité moyenne des tubes NP de la première famille 10 déterminée par le module de calcul à partir de la ou des mesures faites pour chaque famille par la sonde. Ainsi si une seule mesure a été prise par la sonde pour cette première famille, AlP est égale à cette mesure. Si plusieurs mesures ont été prises pour cette première famille, AlP est égale à la moyenne de ces mesures qui est calculée par le module de calcul. 15 ANNFP : l'activité moyenne des tubes NP de la NFième famille déterminée par le module de calcul à partir de la ou des mesures faites pour chaque famille par la sonde. Ainsi si une seule mesure a été prise par la sonde pour cette NFième famille, N ANFP est égale à cette mesure. Si plusieurs mesures ont été prises pour cette NFième famille, ANNPF est égale à la moyenne de ces mesures 20 qui est calculée par le module de calcul. km,NF lm,NF : l'ensemble des indices des tubes NP non mesurés de la NFième famille, par rapport au mième tube P considéré, situé en (im, Jm) en BF ou en BC. L'ensemble des 2 x M équations pour les tubes P peut être écrit sous la forme d'une matrice carrée suivante : [ e00 - - - e00 ekellei Ariii min/ A /mn/ y y Al P Elc11 E111 ek11111 All\\TTIci,NFILAIF eki,NFILNF + - - - + Al Y, NP -kei E tel ekeile A yie i NF --M,NFIM,NF ekM,NFIM,NF z [2] Nous obtenons une équation [2] de la forme x = y + z représentant l'activité des tubes P en fonction du comptage des demi-tubes P et des demi-tubes NP. Avec : - le terme x qui correspond au vecteur des comptages des M tubes P mesuré par la sonde, - le terme y qui se compose de la matrice d'influence des tubes P sur les tubes P dont les coefficients ont été calculés par le module de calcul et du vecteur d'activité des tubes P qui sont les inconnues de l'équation à résoudre - le terme z qui se compose de NF vecteurs des influences des tubes NP sur le comptage des tubes P. Dès lors, une factorisation de Cholesky est envisageable sur ce système pour déterminer le vecteur de l'activité des demi-tubes P. La AF 'n'Ifni. factorisation est réalisée par le module de calcul à partir des valeurs du fichier de mesure stocké dans les moyens de mémoire. Une hypothèse est posée sur les demi-tubes NP. Cette hypothèse consiste à estimer que toutes les activités des demi-tubes NP sont identiques dans une même famille, donc que Al = A0 quelque soit i. Dans cette hypothèse, la relation peut s'écrire pour une seule famille : No = Ao. eoi j=0 avec : - N0 : le taux de comptage mesure dans le demi-tube sondé, - jErl=0 : la somme des coefficients de l'ensemble des demi-tubes, y compris du demi-tube sondé, - A0 : l'activité du demi-tube sondé. Ainsi l'activité pour les demi-tubes NP pour une famille peut être calculée par le module de calcul par la relation : A0 = j=0 e0,1 Cette relation est utilisée pour chaque famille. Le choix de No peut être déterminé par des techniques statistiques, par exemple, de type bootstrap. Le module de calcul détermine ainsi une valeur d'activité Ai du générateur (0) de vapeur. La valeur d'activité est enregistrée dans un fichier de résultats stocké dans le module de stockage numérique. No Malgré la dimension initiale du système à résoudre, la méthode de résolution mise en oeuvre par l'invention dans les moyens de calcul et la diminution de la complexité permettant d'aboutir à des temps de calcul de l'ordre de la seconde. La description fait référence à des tubes ou demi-tubes car c'est généralement la configuration des éléments de tube d'un générateur (0) de vapeur, mais il doit être clair dans l'esprit du lecteur que la méthode de l'invention s'applique quelque soit le type de tube de générateur, qu'ils soient tubulaires ou non. Il doit être évident pour les personnes versées dans l'art que la présente invention permet des modes de réalisation sous de nombreuses autres formes spécifiques sans l'éloigner du domaine d'application de l'invention comme revendiqué. Par conséquent, les présents modes de réalisation doivent être considérés à titre d'illustration, mais peuvent être modifiés dans le domaine défini par la portée des revendications jointes, et l'invention ne doit pas être limitée aux détails donnés ci-dessus. REVENDICATIONS1. Procédé de caractérisation radiologique d'un générateur (0) de vapeur comprenant une pluralité de tubes en U, chaque tube en U étant formé d'au moins deux demi-tubes (11f, 11c) dont les coordonnées d'au moins un demi-tube sont enregistrées dans un fichier de mesure stocké dans un moyen de stockage numérique, le générateur (0) de vapeur comprenant des tubes dits non-pastillés et éventuellement des tubes pastillés, le procédé étant caractérisé en ce qu'il comprend au moins les étapes suivantes : - détermination de familles de demi-tubes parmi les demi-tubes ; - au moins une mesure par une sonde (12a) du taux de comptage Ni du rayonnement gamma dans au moins un demi-tube par famille, les mesures étant enregistrées dans le fichier de mesures en association avec les coordonnées du ou des demi-tubes sondés (11s) ; - détermination pour chaque famille par le module de calcul du ou des coefficients eii d'influence d'une matrice d'influence pour obtenir une matrice ayant une dimension inférieure ou égale au nombre de familles, un coefficient d'influence correspondant à l'influence du taux de comptage d'un terme source de rayonnement gamma sur le taux de comptage du demi-tube sondé, les coefficients d'influence étant enregistrés dans un fichier de calcul stocké dans le moyen de stockage numérique ; - détermination par le module de calcul de l'inverse de la matrice d'influence dont les coefficients sont enregistrés dans le fichier de calcul stocké dans le moyen de stockage numérique ; - détermination par le module de calcul d'une valeur d'activité Ai du générateur (0) de vapeur à partir des mesures de taux de comptage, la valeur d'activité étant enregistrée dans un fichier de résultats stocké dans le module de stockage numérique. 2. Procédé selon la revendication 1, caractérisé en ce que l'étape de détermination de l'inverse de la matrice d'influence par le module de calcul comprend une étape dans laquelle le module de calcul calcule l'inverse la matrice d'influence en considérant que l'activité des demi-tubes dans une famille est la même que celle mesurée pour au moins un demi-tube de cette famille. 3. Procédé selon la revendication 1, caractérisé en ce que l'étape de détermination de l'inverse de la matrice d'influence par le module de calcul comprend une étape dans laquelle le module de calcul calcule l'inverse de la matrice d'influence en considérant que l'activité des demi-tubes varie en fonction des coordonnées des demi-tubes. 4. Procédé selon l'une des revendications 1 à 3, caractérisé en ce que le ou les coefficients d'influences sont déterminés à partir d'une fonction d'interpolation eii = f (dii) exécutée par le module de calcul en utilisant les valeurs des mesures de taux de comptage enregistrées dans le fichier de mesures et au moins des informations représentatives des distances entre le tube sondé et les termes sources enregistrées dans un fichier représentatif du générateur (0) de vapeur mémorisé dans le moyen de stockage numérique. 5. Procédé selon l'une des revendications 1 à 3, caractérisé en ce que le ou les coefficients d'influence sont déterminés en outre à partir d'une fonction d'interpolation eij = f (dij, autres_paramètres) exécutée par le module de calcul en utilisant les valeurs des mesures de taux de comptage enregistrés dans le fichier de mesures et au moins un autre paramètre déterminé à partir de la moyenne d'un échantillonnage d'au moins une caractéristique structurelle du générateur (0) de vapeur. 6. Procédé selon la revendication 5, caractérisé en ce que la caractéristique structurelle du générateur (0) de vapeur est une information représentative de l'épaisseur moyenne EM de matière traversée selon un plan perpendiculaire au tube sondé entre le terme source (12b) et le centre (12a) du demi-tubesondé (11s). enregistré dans un fichier représentatif du générateur (0) de vapeur mémorisé dans le moyen de stockage numérique. 7. Procédé selon l'une des revendications 1 à 6, caractérisé en ce que le terme source (12b) est la partie interne d'un demi-tube (11i) et/ou une source autre que la partie interne d'un demi-tube. 8. Procédé selon l'une des revendications 1 à 7, caractérisé en ce que l'étape de détermination de familles comprend : - la constitution d'au moins une famille de demi-tubes pastillés comprenant chacune un seul demi-tube pastillé, et - la constitution d'au moins une famille de tubes non pastillés comprenant tous les demi-tubes non pastillés similaires en activité radioactive. 9. Procédé selon l'une des revendications 1 à 8, caractérisé en ce que l'étape de détermination de familles comprend : la constitution d'au moins une famille de demi-tubes pastillés comprenant chacune un seul demi-tube pastillé, et - la constitution d'au moins une famille de demi-tubes non pastillés de telle manière que chaque famille de demi-tubes non pastillés comprend des demi-tubes non pastillés adjacents et ayant au moins un critère de classification semblable. 10. Procédé selon l'une des revendications 1 à 9, caractérisé en ce qu'un moyen d'affichage du module de calcul affiche chaque demi-tube mesuré du générateur (0) de vapeur en association de ses coordonnées, la constitution des familles étant effectuée par l'utilisateur qui utilise des moyens interactifs du module de calcul pour sélectionner sur le moyen d'affichage, pour chaque famille, les demi-tubes à affecter à cette famille à l'aide d'un moyen de sélection. 11. Procédé selon au moins l'une des revendications 1 à 10, caractérisé en ce que lors de l'étape de détermination par le module de calcul de l'inverse dela matrice d'influence comprend une étape de factorisation de type Cholesky de la matrice d'influence par le module de calcul.