Brevet FR2950703A1 — publié le 20110401

PROCEDE DE DETERMINATION DE RAPPORT ISOTOPIQUE DE CHAMBRE A FISSION DESCRIPTION Domaine technique et art antérieur L'invention concerne un procédé de détermination de rapport isotopique de chambre à fission. La mesure par chambre à fission permet de remonter à la puissance d'un réacteur nucléaire ou d'obtenir des taux de fission caractéristiques d'un flux neutronique à certains emplacements du coeur sur une gamme énergétique donnée (cas d'une chambre à fission miniature utilisée dans un réacteur expérimental pour mesurer une atténuation de flux en réflecteur par exemple). Une chambre à fission est constituée d'une enceinte étanche remplie d'un gaz inerte à fort pouvoir d'ionisation. De la matière fissile est déposée à l'intérieur de la chambre à fission sur un élément formant anode, l'enceinte étanche formant cathode. Lorsque les neutrons dont le flux est à mesurer entrent en contact avec la matière fissile, la fission se produit et il y a émission de deux fragments de fission qui partent à 180° l'un de l'autre. L'un des deux fragments est stoppé dans l'anode et l'autre ionise le gaz, avant d'être arrêté par la cathode. L'ionisation du gaz conduit à la création d'électrons dans la chambre à fission. Les électrons créés par ionisation sont alors collectés par la cathode et transmis à une chaîne de traitement électronique via un câble coaxial. Le nombre d'électrons créés est directement proportionnel au flux de neutrons qui est à l'origine de la fission. La connaissance des impuretés présentes dans une chambre à fission est cruciale pour une interprétation correcte des mesures. Dans le cas où la matière fissile est, par exemple, de l'Uranium 238U, il est possible que plusieurs dizaines de pourcents du signal détecté provienne d'impuretés d'Uranium 235U. Il est ainsi nécessaire de connaître le rapport isotopique 235U/238U avec la plus grande précision afin de pouvoir extraire le signal utile du signal détecté. L'incertitude associée à la connaissance de ce rapport isotopique doit donc être la plus faible possible, idéalement de l'ordre de quelques pourcents. A ce jour, plusieurs méthodes sont connues pour déterminer ce rapport isotopique. Une première méthode consiste à dissoudre le matériau concerné ou un échantillon de ce matériau et à effectuer une analyse chimique du matériau dissous. Par cette méthode, il est possible d'atteindre la connaissance de l'isotopie du matériau avec une précision inférieure au pourcent. Un inconvénient majeur de cette méthode est son caractère destructeur du matériau. En plus de son caractère destructeur, un autre inconvénient de cette méthode réside dans la génération d'effluents qu'il faut - in fine - traiter. Une autre méthode est l'interrogation neutronique. L'interrogation neutronique présente l'avantage de ne pas être destructive. Un inconvénient de l'interrogation neutronique est toutefois de ne pas pouvoir être réalisée sur une faible quantité de matériau (quelques centaines de microgrammes dans le cas d'une chambre à fission miniature). D'autres méthodes sont également connues, comme, par exemple, l'étalonnage par rapport à une chambre étalon ou l'utilisation de spectres différents et de matériaux ayant des caractéristiques neutroniques particulières. Ces autres méthodes présentent l'inconvénient d'être difficiles à mettre en oeuvre et de fournir des résultats relativement imprécis. L'invention ne présente pas les inconvénients mentionnés ci-dessus. L'invention propose en effet de mesurer le rapport isotopique en mettant en oeuvre un procédé qui ne détruit pas la chambre à fission et dont l'incertitude est de l'ordre de quelques pourcents (typiquement inférieure à 10%). Exposé de l'invention L'invention concerne un procédé de détermination de rapport isotopique de chambre à fission contenant de la matière fissile ayant un isotope principal X et au moins un isotope d'impureté Y, les isotopes X et Y ayant une décroissance radioactive selon les deux équations suivantes : X -> X' caractérisée par XX, FX, et Y -> Y' caractérisée par Xy, Fy, où X' et Y' sont des isotopes « fils » respectifs des isotopes X et Y, la décroissance de l'isotope X (respectivement Y) étant caractérisée par l'émission d'une particule gamma par l'isotope fils X' (respectivement Y') à une énergie E1 (respectivement E2) suivant une probabilité d'émission I7(El) (respectivement I7(E2)) , les grandeurs À et À étant, respectivement, la constante de décroissance radioactive de l'isotope principal X et la constante de décroissance radioactive de l'isotope d'impureté Y, et FX et Fy étant, respectivement, un facteur de branchement de décroissance d'isotope utilisé pour une mesure d'activité de l'isotope principal et un facteur de branchement de décroissance d'isotope utilisé pour une mesure d'activité de l'isotope d'impureté. Le procédé est caractérisé en ce qu'il comprend les étapes suivantes : mesure, à l'aide d'un banc de spectrométrie placé dans une configuration de mesure donnée, d'une surface nette S(E1) d'un premier pic gamma de la matière fissile à une première énergie E1 et d'une surface nette S(E2) d'un deuxième pic gamma de la matière fissile à une deuxième énergie E2, détermination, à l'aide de sources ponctuelles de référence, dans la configuration de mesure donnée, d'un rendement d'absorption totale de référence RoP(E1) à la première énergie E1 et d'un rendement d'absorption totale de référence RoP (E2) à la deuxième énergie E2, calcul, à l'aide d'un calculateur, dans la configuration de mesure donnée, d'un transfert de rendement intégral T(E1) de la matière fissile à la première énergie E1 et d'un transfert de rendement intégral T(E2) de la matière fissile à la deuxième énergie E2, et calcul, à l'aide d'un calculateur, du rapport isotopique de chambre à fission R à l'aide de l'équation : R-xXS(E2)X1 (El)XRo(E1)XT(El)x X . S(E1) Ir (E2 ) Ro (E2 ) T(E2) Fy Par banc de spectrométrie placé dans une configuration de mesure donnée, il faut entendre que les différents éléments qui constituent le banc de spectrométrie sont placés, les uns par rapport aux autres, selon une géométrie qui n'est pas modifiée d'une mesure à l'autre. Les mesures de surface nette sont effectuées soit avant que la matière fissile ne soit introduite dans la chambre à fission, soit après introduction de la matière fissile dans la chambre à fission. Dans le premier cas, les mesures sont effectuées soit sur un échantillon du bloc de matière fissile, soit sur le bloc de matière fissile dans son entier. Dans le deuxième cas, les mesures sont effectuées sur la chambre à fission dans son entier. Le procédé de l'invention permet avantageusement une détermination facile et performante du rapport isotopique d'une chambre à fission sur la base de mesures de rayonnements gamma et de données de bibliothèque. De façon avantageuse également, la mesure des surfaces nettes est effectuée pendant une durée suffisamment longue qui permet d'obtenir une faible incertitude de comptage (la durée de mesures des surfaces nettes est en effet comprise entre une heure et quelques semaines, par exemple entre une heure et dix semaines). Brève description des figures D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture d'un mode de réalisation préférentiel fait en référence aux figures jointes parmi lesquelles : - La figure 1 représente un schéma de principe de dispositif utilisé pour la mise en oeuvre du procédé de l'invention ; - La figure 2 représente un schéma de principe du procédé de l'invention ; - La figure 3 représente une vue de détail d'une première étape du procédé de l'invention ; - La figure 4 représente une vue de détail d'une deuxième étape du procédé de l'invention; et - La figure 5 représente une vue de détail d'une troisième étape du procédé de l'invention. Exposé détaillé d'un mode de réalisation préférentiel de l'invention Selon le mode de réalisation préférentiel de l'invention décrit ci-dessous, le matériau fissile est l'Uranium 238U et les impuretés sont constituées d'Uranium 23'U. L'invention concerne cependant tout autre matériau fissile ainsi que toute impureté associée à cet autre matériau fissile. Le matériau fissile peut être, par exemple 240Pu, 242Pu, 243Am2232Th et l'impureté associée peut être, par exemple, 239Pu, 241Am2233Th. Dans le cas où le matériau fissile est l'Uranium 238U et où les impuretés sont constituées d'Uranium 235U, le rapport isotopique R est donné par l'équation (1) ci-dessous: R= U238 X 4U235 XU235 Au238 où Au238 et AU235 sont, respectivement, l'activité de l'Uranium 238U et l'activité de l'Uranium 235U et Au238 et AU235 sont, respectivement, la constante de décroissance radioactive de l'Uranium 238U et la constante de décroissance radioactive de l'Uranium 235U Le schéma de décroissance de l'Uranium 2 et de l'uranium 235U s'écrit : 224U 0,9985 2 a 3'Th L:i4m Pa T1/2=1, 17 minutes T1/2=4, 468*10" ans T1/2=24, 09 jours FPa = 0,0015 4 2~~U 2'Th T1/2=7.038*10 ans où le symbole « a » représente la décroissance par émission de particules a, le symbole « R » la décroissance par émission de particules R, T1/2 la demi-vie de l'isotope concerné, Fu le facteur de branchement de décroissance de l'isotope 234mPa vers 234U, c'est-à-dire la probabilité pour le Protactinium 234mPa de se désintégrer en Uranium 234U et FPa le facteur de branchement de décroissance de l'isotope 234mPa vers 234Pa, c'est-à-dire la probabilité pour le Protactinium 234mPa de se désintégrer en Protactinium 234Pa (Fu+FPa=1) . (1) A l'équilibre séculaire (également appelé équilibre nucléaire), le schéma de décroissance de l'Uranium 238U conduit donc à exprimer l'activité AU238 de l'Uranium 238U par rapport à l'activité APa234m du Protactinium 234mPa comme suit : 1 X APa234m F ( 2 ) U D' où : R = U238 X AU235 X Fu ?LU235 APa234m Les deux raies gamma traditionnellement mesurées pour quantifier les activités de l'Uranium 23'U et du Protactinium 234mPa sont les raies à 185,7keV (intensité de 57,0 %) et à 1001keV (intensité de 0,839%). Ce sont ces deux raies qui sont mesurées dans le cadre du mode de réalisation préférentiel de l'invention. L'activité A d'un isotope quelconque calculée à partir de la mesure d'une de ses raies d'émission gamma d'énergie Eo est donnée par la formule suivante : Ir(Eo)RoP (Eo)T(Eo)Ataet (4) où . S(Eo) est la surface nette du pic gamma à l'énergie Eo, dans la configuration de mesure donnée; Ccoinc (Eo) est la correction de coïncidences vraies (détection simultanée de rayonnements émis lors d'une même désintégration radioactive), dans la configuration de mesure donnée ; AU238 (3) A = S(Eo)Ccoine(E0) I7(Eo) est la probabilité d'émission ou intensité de la raie gamma d'énergie Eo pour l'isotope fils de l'isotope considéré; RoP(Eo) est le rendement d'absorption totale de référence à l'énergie Eo, qui est obtenu par étalonnage du détecteur gamma à l'aide d'une source étalon ponctuelle, dans la configuration de mesure donnée ; T(Eo) est le transfert de rendement intégral qui permet de corriger les effets d'angle solide et de matrice dans la configuration de mesure donnée (hétérogénéité des sources et autoabsorption gamma dus à la mesure d'un objet dans une configuration différente de celle utilisée lors de l'étalonnage) ; Atact est la durée de mesure active, c'est-à-dire la durée de mesure corrigée du temps mort électronique. De l'expression de A ci-dessus, il découle que le rapport isotopique R pour l'Uranium 238U s'écrit : R= U238 x S(185.7keV) x Ir (1001keV) x Ro (1001keV) x T (1001keV) x F (5) 41235 S(1001keV) Ir (185,7keV) Ro (185,7keV) T(185,7keV) U où . - Au235 est une donnée connue dans les bases de données internationales (Au235=ln (2) / (2,22102 1016) s-1) ; - Au238 est une donnée connue dans les bases de données internationales (Au238=ln (2) / (1,40996 1017) s-1) ; S(185,7keV) est une grandeur mesurée; S(1001keV) est une grandeur mesurée; - IY(1001keV) est une donnée connue dans les bases de données internationales (IY(1001kev)=0,00839) ; - IY(185,7keV) est une donnée connue dans les bases de données internationales (IY(185,7kev)=0,570) ; Ro''(1001keV) est une grandeur mesurée par étalonnage; Ro''(185,7keV) est une grandeur mesurée par étalonnage ; T(1001keV) est une grandeur calculée par modélisation ; T(185,7keV) est une grandeur calculée par modélisation ; Fu est une donnée connue dans les bases de données internationales. L'équation (5) ci-dessus concerne le cas où le matériau fissile est l'Uranium 238U (isotope principal) et l'impureté l'Uranium 235U (isotope d'impureté associé), les raies gamma utilisées pour quantifier les activités de l'isotope principal et de l'isotope d'impureté étant respectivement les raies à 1001keV et 185,7keV. Dans le cas général d'un matériau fissile X (isotope principal) et d'une impureté associée Y (isotope d'impureté), les raies gamma utilisées pour quantifier les activités de l'isotope principal et de l'isotope d'impureté étant, respectivement, une raie d'énergie E1 et une raie d'énergie E2, le rapport isotopique R s'écrit : R-~xXS(E2)XIr(Ei)XRo(Ei)XT(El)XFx (6) X . S (E1) I r (E2 ) Ro (E2 ) T (E2 ) Fy où . À est la constante de décroissance radioactive de l'isotope principal de la matière fissile destinée à équiper la chambre à fission, À est la constante de décroissance radioactive de l'isotope d'impureté associé à l'isotope principal, ù S(E1) est la surface nette d'un premier pic gamma de la matière fissile à l'énergie E1 et S(E2) la surface nette d'un deuxième pic gamma de la matière fissile à l'énergie E2, ù RoP(E1) est le rendement d'absorption totale de référence à l'énergie E1 et RoP (E2) le rendement d'absorption totale de référence à l'énergie E2, ù T(E1) est le rendement intégral de la matière fissile à l'énergie E1 et T(E2) le rendement intégral de l'échantillon de la matière fissile à l'énergie E2, ù IY (E1) et IY (E2) sont, respectivement, la probabilité d'émission de raie gamma du « fils » de l'isotope principal à l'énergie E1 et la probabilité d'émission de raie gamma du « fils » de l'isotope d'impureté à l'énergie E2, FX est le facteur de branchement de décroissance d'un isotope utilisé pour mesurer l'activité de l'isotope principal, et ù Fy est le facteur de branchement de décroissance d'un isotope utilisé pour mesurer l'activité de l'isotope d'impureté. Comme cela a été mentionné précédemment avec l'exemple de l'Uranium 238U et de l'Uranium 235U, les facteurs de branchement de décroissance sont des probabilités connues liées au mode de désintégration des isotopes. Dans le cas où l'isotope utilisé pour mesurer l'activité se désintègre de manière univoque, le facteur de branchement de décroissance est égal à 1 et, dans le cas où la désintégration n'est pas univoque, le facteur de branchement est égal à la probabilité associée à la désintégration choisie. La figure 1 représente un schéma de principe de dispositif utilisé pour la mise en oeuvre du procédé de l'invention. Le dispositif comprend un banc Sp de spectrométrie gamma et un calculateur K. Le banc Sp de spectrométrie gamma délivre des signaux de mesure Sm qui sont transmis au calculateur K. Le banc de spectrométrie Sp est constitué d'une enceinte E dans laquelle est placé un détecteur D. Le détecteur D, préférentiellement une diode en germanium de grande pureté, est un détecteur de rayons gamma qui présente une bonne efficacité sur un large domaine d'énergie, par exemple des énergies allant de 50keV à 2MeV. De façon préférentielle, le détecteur D présente également une aptitude à limiter l'effet Compton. De la matière fissile est placée dans l'enceinte E. La référence B sur la figure 1 représente cette matière fissile : il s'agit soit d'un échantillon du bloc de matière fissile destiné à équiper la chambre à fission, soit du bloc de matière fissile lui-même, soit de la chambre à fission dans son entier équipée du bloc de matière fissile. L'interaction des rayons gamma issus de la matière fissile B avec le détecteur D conduit à une ionisation des atomes du milieu détecteur. Les charges ainsi crées sont collectées par une alimentation continue haute tension, par exemple une alimentation de quelques milliers de volts (non représentée sur la figure 1). L'enceinte E est préférentiellement une enceinte circulaire, constituée de matériaux permettant d'avoir un niveau très faible de bruit de fond, par exemple, un empilement de plomb basse activité, de plomb ultra basse activité, d'étain et de cuivre. Le refroidissement du détecteur est assuré par un système cryo-cycle qui permet un refroidissement hybride électrique/azote liquide. Ce système permet des durées d'acquisition de plusieurs semaines sans perturbation liée à l'apport en azote liquide. Les signaux de mesure Sm construits à partir du rayonnement gamma sont transmis au calculateur K pour traitement. Le calculateur K met en oeuvre des procédés de traitement aptes à délivrer les grandeurs souhaitées. Le procédé de l'invention va maintenant être décrit en référence aux figures 2-5. Le procédé de l'invention comprend (cf. figure 2) une étape 1 de détermination des surfaces nettes S(185,7keV) et S(1001keV) de la matière fissile, une étape 2 de détermination des rendements d'absorption totale de référence aux énergies de 185,7keV et de 1001keV, Ro''(185,7keV) et Ro''(1001keV), une étape 3 de détermination des transferts de rendement intégral aux énergies de 185,7keV et de 1001keV, T(185,7keV) et T(1001keV), et une étape 4 de calcul du rapport isotopique R à partir des résultats des étapes 1, 2 et 3 et de données D connues issues de bibliothèques de données. Les données D sont les données À 2m, À 2m, IY(185,7keV), IY(1001keV) et Fu définies précédemment. L'étape 1 de détermination des surfaces nettes S(185,7keV) et S(1001keV) comprend une étape la de mesure de spectre gamma de la matière fissile B. La mesure du spectre gamma est effectuée, à l'aide du banc de mesure Sp, sur une bande d'énergie qui comprend les raies utiles à 185,7keV et 1001keV. L'étape la (Figure 3) est suivie d'une étape lb d'extraction, à l'aide du calculateur K, des surfaces nettes des raies à 185,7keV et 1001keV. Une mesure de bruit de fond gamma est également effectuée, à l'aide du banc de mesure Sp, sur la même bande d'énergie lors de l'étape 1N, pour mesurer les éventuels pics gamma du bruit de fond aux énergies E1 et E2 se superposant aux pics déterminés dans la mesure de l'étape la. A partir des résultats obtenus à l'issue des étapes lb et 1N, le calculateur K calcule, à l'étape 1c, les surfaces nettes S(185,7keV) et S(1001keV). Le logiciel de calcul utilisé pour le calcul des surfaces nettes est, par exemple, le logiciel GENIE2000 dont le manuel d'utilisation est disponible sur le réseau Internet à l'adresse électronique suivante : http://www.ipnas.ulg.ac.be/garnir/pdf/genie2000.pdf, Tout autre type de logiciel connu de l'homme du métier et susceptible de conduire au résultat souhaité peut également être utilisé. La figure 4 illustre les différentes étapes élémentaires qui composent l'étape 2 de détermination des rendements d'absorption totale de référence aux énergies de 185,7keV et de 1001keV, Ro''(185,7keV) et RoP (1001keV) . L'étape 2 de détermination des rendements d'absorption totale de référence Ro''(185,7keV) et Ro''(1001keV) est une étape connue en soi. L'étape 2 commence par une étape 2a de choix de sources ponctuelles dont les activités sont bien connues. Chaque source ponctuelle j est caractéristique d'une émission de particules gamma à une énergie Ei choisie dans une gamme énergétique connue, par exemple la gamme 50keV-2MeV. Chaque source ponctuelle j est placée dans l'enceinte pour la mesure de son activité (étape 2b). L'activité AmP(Ej) de chaque source ponctuelle j à l'énergie Ej est alors mesurée à l'aide du banc de mesure Sp et du calculateur K. Le rendement d'absorption totale de référence R'oP(Ej) de chaque source ponctuelle à l'énergie Ej est ensuite calculé, à l'étape 2c, à partir de l'activité mesurée AmP (Ei) et de l'activité de la source ponctuelle qui est connue a priori AoP (Ei) . Il vient . R' of' (Fi) = AmP(Fi)/ AoP(Fi), On obtient alors une courbe des rendements d'absorption totale de référence R'oP(Ej) pour les différentes énergies E. La courbe des rendements R'oP(Ej) est ensuite ajustée, à l'aide d'une fonction analytique Fa, par le calculateur K (étape 2d). La fonction analytique Fa consiste par exemple à exprimer le logarithme népérien de R'oP en fonction d'un polynôme de logarithme népérien de l'énergie. De la courbe ajustée délivrée à l'étape 2d sont déduites, à l'étape 2e, les valeurs des rendements d'absorption totale de référence aux énergies souhaitées de 185,7keV et 1001keV, RoP (185, 7keV) et RoP (1001keV) . La figure 5 illustre les différentes étapes élémentaires qui composent l'étape 3 de détermination des transferts de rendement intégraux. L'étape 3 de détermination des transferts de rendement intégraux comprend, tout d'abord, une étape 3a de modélisation du banc de mesure équipé du détecteur D, une étape 3b de modélisation de l'élément mesuré lors de l'étape 1 de détermination des surfaces nettes (chambre à fission CH ou échantillon de bloc de matière fissile) et une étape 3c de modélisation des conditions d'utilisation de la source ponctuelle utilisée lors de la mesure des rendements d'absorption dans la même configuration de mesure du banc gamma que dans l'étape 3b. A partir des données de modélisation délivrées à l'issue des étapes 3a et 3b, une étape de calcul 3d calcule, par exemple par la méthode de Monte Carlo pour la résolution de l'équation de transport des photons, l'activité relative de l'élément modélisé déterminée par le détecteur D modélisé dans le calcul aux énergies respectives de 185,7keV et 1001keV, à savoir les grandeurs AccH (185, 7keV) et AccH (1001keV) . De même, à partir des données de modélisation délivrées à l'issue des étapes 3a et 3c, l'étape de calcul 3d calcule, par exemple par la méthode de Monte Carlo, l'activité relative de la source ponctuelle déterminée par le détecteur D modélisé dans le calcul aux énergies de 185,7keV et 1001keV, à savoir les grandeurs AcP (185, 7keV) et AcP (1001keV) . A l'issue de l'étape 3d, une étape de calcul 3e calcule les grandeurs T(185,7keV) et T(1001keV) à partir des activités calculées, à savoir : 1(185, 7keV) =ACCH (185, 7keV) / Aj (185, 7keV) T (1001keV) =AccH (1001keV) / AcP (1001keV) A partir de l'ensemble des grandeurs déterminées comme décrit ci-dessus et des données de bibliothèque mentionnées précédemment, il est alors possible de déterminer le rapport isotopique R tel que : R _ U238 X S(185.7keV) Ir (1001keV) X Ro (1001keV) X T(1001keV) Xrr235 S(1001keV) Ir (185,7keV) Ro (185,7keV) T(185,7keV ou, dans le cas général : R = ~x x S(Ez) Ir(E1) x RoP(Ei) X T(E1) X FX . X S(E1) Ir (E) Ro (E2) T(E2) FY REVENDICATIONS1. Procédé de détermination de rapport isotopique de chambre à fission contenant de la matière fissile ayant un isotope principal X et au moins un isotope d'impureté Y, les isotopes X et Y ayant une décroissance radioactive selon les deux équations suivantes . X -> X' caractérisée par XX, FX, et Y -> Y' caractérisée par Xy, Fy, où X' et Y' sont des isotopes « fils » respectifs des isotopes X et Y, la décroissance de l'isotope X (respectivement Y) étant caractérisée par l'émission d'une particule gamma par l'isotope fils X' (respectivement Y') à une énergie El (respectivement E2) suivant une probabilité d'émission I7(El) (respectivement I7(E2)) , les grandeurs À et À étant, respectivement, la constante de décroissance radioactive de l'isotope principal X et la constante de décroissance radioactive de l'isotope d'impureté Y, et Fy et Fy étant, respectivement, un facteur de branchement de décroissance d'isotope utilisé pour une mesure d'activité de l'isotope principal et un facteur de branchement de décroissance d'isotope utilisé pour une mesure d'activité de l'isotope d'impureté, caractérisé en ce qu'il comprend les étapes suivantes : mesure, à l'aide d'un banc de spectrométrie placé dans une configuration de mesure donnée, d'une surface nette S(E1) d'un premier pic gamma de la matière fissile à une première énergie E1 et d'unesurface nette S(E2) d'un deuxième pic gamma de la matière fissile à une deuxième énergie E2, détermination, à l'aide de sources ponctuelles de référence, dans la configuration de mesure donnée, d'un rendement d'absorption totale de référence RoP(E1) à la première énergie E1 et d'un rendement d'absorption totale de référence RoP (E2) à la deuxième énergie E2, calcul, à l'aide d'un calculateur, dans la configuration de mesure donnée, d'un transfert de rendement intégral T(E1) de la matière fissile à la première énergie E1 et d'un transfert de rendement intégral T(E2) de la matière fissile à la deuxième énergie E2, et calcul, à l'aide d'un calculateur, du rapport isotopique de chambre à fission R à l'aide de l'équation : R= X S(E2) X Ir (E1) X Ro (E1) X T (Ei) x Fx X . S(E1) Ir (E2 ) Ro (E2 ) T(E2) Fy 2. Procédé selon la revendication 1, dans lequel les mesures de surface nette de la matière fissile sont effectuées, avant introduction de la matière fissile dans la chambre à fission, sur un échantillon du bloc de la matière fissile destiné à équiper la chambre à fission. 3. Procédé selon la revendication 1, dans lequel les mesures de surface nette de la matière fissile sont effectuées, avant introduction de la matière fissile dans la chambre à fission, sur le blocde matière fissile destiné à équiper la chambre à fission. 4. Procédé selon la revendication 1, dans lequel les mesures de surface nette de la matière fissile sont effectuées, après introduction de la matière fissile dans la chambre à fission, sur la chambre à fission dans son entier. 5. Procédé selon l'une quelconque des revendications précédentes, dans lequel l'isotope principal est l'Uranium 238U et l'isotope d'impureté est l'Uranium 23"U 6. Procédé selon l'une quelconque des revendications précédentes, dans lequel la durée des étapes de mesure de surface nette est comprise entre une heure et dix semaines.