Brevet FR3090128A1 — publié le 20200619

Description Titre de l'invention : PROCEDE DE CARACTERISATION D'UN SPECTROMETRE, PRODUIT PROGRAMME D'ORDINATEUR ET CALCULATEUR ASSOCIES Domaine technique 100011 La présente invention concerne un procédé d'étalonnage d'un spectromètre comprenant un détecteur, le spectromètre étant associé à une plage d'énergie de détection comprenant n, canaux de détection. 100021 L'invention concerne également un produit programme d'ordinateur et un cal- culateur. 100031 L'invention s'applique au domaine de la physique nucléaire, notamment à la caracté- risation de spectromètres dédiés à l'instruction nucléaire, plus particulièrement à l'étalonnage en rendement desdits spectromètres. Technique antérieure [0004] De façon classique, un spectromètre X et/ou y est un dispositif de mesure nucléaire comprenant, d'une part, une chaîne de détection comportant un détecteur et un organe d'acquisition, et, d'autre part, un analyseur configuré pour analyser un signal de détection généré par la chaîne de détection, l'analyseur étant associé à une plage d'énergie de détection prédéterminée, subdivisée en un nombre entier prédéterminé n, de canaux de détection. La correspondance entre chaque canal de détection de l'analyseur et l'intervalle correspondant de la plage d'énergie de détection est établie au cours d'un étalonnage en énergie. [0005] Le spectromètre est configuré pour fournir, suite à l'exposition du détecteur à un rayonnement ionisant, un spectre du rayonnement, c'est-à-dire un histogramme de la population de photons détectés en fonction de l'énergie déposée dans le détecteur par lesdits photons, conformément à la subdivision en n' de canaux de détection de la plage d'énergie de détection prédéterminée. [0006] Un tel spectromètre, pour pouvoir remplir une fonction de quantification de l'activité radiologique, requiert, outre son étalonnage en énergie, un étalonnage dit « en rendement de mesure en absorption totale », encore appelé « étalonnage en rendement ». [0007] Un tel étalonnage en rendement consiste à associer, à chacun des n, canaux de détection, et pour une configuration expérimentale donnée (définie, notamment, par une position de la source par rapport au détecteur, une géométrie de la source, et une influence de l'environnement), un rendement correspondant de l'ensemble du spectromètre. En d'autres termes, une fois l'étalonnage réalisé, si la configuration expé- 2 rimentale demeure inchangée, la détermination du nombre d'impulsions détectées dans chaque canal de détection donne accès à l'activité d'une source quelconque. [0008] Toutefois, il est malaisé d'établir un étalonnage en rendement d'un spectromètre donné qui décrive exhaustivement, pour un radio-isotope donné, l'ensemble des configurations expérimentales. C'est la raison pour laquelle est généralement utilisé un modèle numérique du spectromètre. [0009] Pour établir un tel modèle numérique, il est connu d'avoir recours à des logiciels qui, moyennant la fourniture en entrée de valeurs mesurées, lors de la mise en oeuvre de configurations expérimentales prédéterminées, de paramètres prédéterminés, analysent la réponse du spectromètre dans chacune desdites configurations expérimentales prédéterminées afin de déterminer l'impact de chaque paramètre et générer le modèle numérique recherché à partir de solutions préenregistrées. [0010] Néanmoins, une telle approche ne donne pas entière satisfaction. [0011] En effet, une telle approche n'offre pas à un utilisateur possédant une expertise et des outils de calcul adéquats la possibilité de modéliser lui-même l'environnement et l'exposition. Or, il est parfois difficile de faire correspondre une configuration expérimentale à une configuration expérimentale prédéterminée attendue. En outre, une telle approche suppose que les valeurs fournies des paramètres prédéterminés soient justes et précises. Or les résultats de la procédure d'étalonnage peuvent se révéler extrêmement sensibles à toute variation desdits paramètres, et les paramètres internes du spectromètre (par exemple, des dimensions du détecteur) ne peuvent, dans la plupart des cas, être déterminés qu'avec de larges incertitudes relatives. [0012] Un but de l'invention est donc de proposer un procédé d'étalonnage qui soit plus flexible pour l'utilisateur, et qui puisse être mis en oeuvre même si les valeurs de pa- ramètres d'un modèle ne peuvent pas être obtenues avec une grande précision. Exposé de l'invention [0013] A cet effet, l'invention a pour objet un procédé d'étalonnage du type précité, comprenant les étapes suivantes : - pour chacune parmi n, configurations expérimentales dans lesquelles le détecteur est exposé à au moins une source de rayonnements ionisants, et pour chacun des n, canaux de détection, détermination, en fonction d'un signal de détection obtenu au moyen du spectromètre, d'un rendement de mesure expérimental correspondant ; - mise en oeuvre, successivement au cours du temps, d'une étape de calcul comprenant : - une phase d'élaboration de n', modèles du spectromètre, l'élaboration de chacun des n', modèles comportant le tirage, pour chaque paramètre prédéterminé parmi n' paramètres prédéterminés, d'une valeur respective suivant une loi de probabilité courante 3 correspondant au paramètre prédéterminé ; - pour chacun des n', modèles élaborés, une phase de détermination, pour chacune des n, configurations expérimentales, d'un rendement de mesure simulé de chacun des n, canaux de détection ; - pour chacun des n', modèles élaborés, une phase de calcul, à partir du rendement de mesure simulé et du rendement de mesure expérimental associé à chacune des n, configurations expérimentales et à chacun des n, canaux de détection, d'un critère de mérite correspondant ; - pour chacun des np paramètres prédéterminés, une phase d'adaptation des grandeurs caractéristiques de la loi dc probabilité courante correspondante, à partir de l'évolution de la valeur du critère de mérite en fonction des n', valeurs prises par le paramètre prédéterminé ; - lorsqu'une condition prédéterminée est remplie, étape de choix, pour chacun des np paramètres prédéterminés, d'une valeur étalonnée respective à partir de la dernière loi de probabilité courante correspondante ; n' n' n, et n', étant des entiers prédéterminés. [0014] En effet, avec un tel procédé d'étalonnage, l'utilisateur est en mesure d'utiliser toute configuration expérimentale que les circonstances ou que son expérience le poussent à mettre en place, de sorte qu'il n'est plus contraint par des configurations expérimentales prédéterminées potentiellement difficiles à mettre en oeuvre. [0015] En outre, un tel procédé d'étalonnage autorise un plus grande tolérance dans les valeurs renseignées des paramètres du modèle, dans la mesure où la valeur de chaque paramètre est mise à jour (par l'intermédiaire de la mise à jour de la loi de probabilité correspondante) durant l'exécution du procédé objet de l'invention, ladite mise à jour étant guidée par une optimisation du critère de mérite obtenu à partir des résultats fournis par le modèle et des résultats simulés (en l'occurrence, le rendement de mesure simulé et le rendement de mesure expérimental). [0016] Le procédé objet de l'invention est donc plus flexible pour l'utilisateur, et est sus- ceptible d'être mis en oeuvre même si les valeurs de paramètres d'un modèle ne peuvent pas être obtenues avec une grande précision, tout en conduisant à un résultat [0017] Suivant d'autres aspects avantageux de l'invention, le procédé comporte une ou plusieurs des caractéristiques suivantes, prise(s) isolément ou suivant toutes les combinaisons techniquement possibles : [0018] - la phase d'adaptation comporte, pour chacun des np paramètres prédéterminés : - pour chacun des n,' modèles élaborés, une association de la valeur correspondante du paramètre prédéterminé et du critère de mérite correspondant ; - une modélisation, à partir des na, valeurs prises par le paramètre prédéterminé et des n', critères de mérite correspondants, de l'évolution du critère de mérite en fonction de la valeur du paramètre prédéterminé ; - à partir du résultat de la modélisation, une modification des grandeurs caractéristiques de la loi de probabilité courante correspondante ; [0019] - pour chacun des np paramètres prédéterminés, si la valeur d'un facteur de forme dépendant d'un résultat de la modélisation correspondante est supérieure à une première valeur prédéterminée, la modification des grandeurs caractéristiques de la loi de probabilité courante correspondante est une réduction de l'écart-type de la loi de probabilité et/ou une modification de l'espérance de la loi de probabilité ; [0020] - pour chacun des np paramètres prédéterminés, si la valeur d'un facteur de forme dépendant d'un résultat de la modélisation correspondante est inférieure à une deuxième valeur prédéterminée, la modification de la loi de probabilité courante correspondante est une augmentation de l'écart-type de la loi de probabilité ; [0021] - la condition prédéterminée est remplie si l'étape de calcul est mise en oeuvre un nombre prédéterminé de fois, et/ou si, pour chacun des n1, paramètres prédéterminés, la variation relative de toute ou partie des grandeurs caractéristiques de la loi de probabilité correspondante entre deux mises en oeuvres successives de l'étape de calcul est inférieure ou égale à un seuil prédéterminé ; [0022] - pour chacun des n1, paramètres prédéterminés, la valeur étalonnée respective choisie est l'espérance de la dernière loi de probabilité courante correspondante. [0023] L'invention a également pour objet un produit programme d'ordinateur comportant des instructions de code de programme qui, lorsqu'elles sont exécutées par un ordinateur, mettent en oeuvre le procédé tel que défini ci-dessus. [0024] En outre, l'invention a pour objet un calculateur pour l'étalonnage d'un spectromètre comprenant un détecteur, le spectromètre étant associé à une plage d'énergie de détection comprenant n, canaux de détection, le calculateur étant configuré pour : - pour chacune parmi n, configurations expérimentales dans lesquelles le détecteur est exposé à au moins une source de rayonnements ionisants, et pour chacun des n, canaux de détection, déterminer, en fonction d'un signal de détection obtenu au moyen du spectromètre, un rendement de mesure expérimental correspondant ; - mettre en oeuvre, successivement au cours du temps, une étape de calcul comprenant : - une phase d'élaboration de n,' modèles du spectromètre, l'élaboration de chacun des n,' modèles comportant le tirage, pour chaque paramètre prédéterminé parmi np paramètres prédéterminés, d'une valeur respective suivant une loi de probabilité courante correspondant au paramètre prédéterminé ; - pour chacun des mn modèles élaborés, une phase de détermination, pour chacune des n, configurations expérimentales, d'un rendement de mesure simulé de chacun des n' canaux de détection ; - pour chacun des nr, modèles élaborés, une phase de calcul, à partir du rendement de mesure simulé et du rendement de mesure expérimental associé à chacune des n, configurations expérimentales et à chacun des n, canaux de détection, d'un critère de mérite correspondant ; - pour chacun des ni, paramètres prédéterminés, une phase d'adaptation des grandeurs caractéristiques de la loi de probabilité courante correspondante, à partir de l'évolution de la valeur du critère de mérite en fonction des n,' valeurs prises par le paramètre prédéterminé ; - lorsqu'une condition prédéterminée est remplie, choisir, pour chacun des n, paramètres prédéterminés, une valeur étalonnée respective à partir de la dernière loi de probabilité courante correspondante ; n,' n,, n, et n', étant des entiers prédéterminés. Brève description des dessins [0025] L'invention sera mieux comprise à l'aide de la description qui va suivre, donnée uniquement à titre d'exemple non limitatif et faite en se référant aux dessins annexés sur lesquels : [0026] [fig.1] est une représentation schématique d'un spectromètre associé à un système d'étalonnage selon l'invention ; [0027] If ig.21 est un ordinogramme illustrant schématiquement le procédé d'étalonnage selon l' invention ; et [0028] Ifig.31 est un graphique représentant l'évolution du biais moyen de modèles d'un spec- tromètre en fonction der indice d'itération d'une étape de calcul du procédé d'étalonnage selon l'invention. [0029] EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS [0030] Un spectromètre 2 est représenté sur la figure I. Le spectromètre 2 est associé à un système d'étalonnage 4 selon l'invention. [0031] Le spectromètre 2 est configure pour fournir, suite à son exposition à un rayonnement, notamment ionisant, issu d'une source 5, un spectre dudit rayonnement. En outre, le système d'étalonnage 4 est configure pour réaliser un étalonnage en rendement du spectromètre 2. Le système d'étalonnage 4 est assimilable à un calculateur. [0032] Par « calculateur », il est par exemple entendu, au sens de la présente invention, un ordinateur, ou encore un circuit de type ASIC (de l'anglais « application-specific in- 6 tegrated circuit », signifiant circuit intégré propre à une application) ou un circuit logique programmable (en anglais, « programmable logic device »), tel qu'un circuit FPGA (de l'anglais «fie/d-programmable gaie array », signifiant réseau de portes programmables in situ). [0033] Le spectromètre 2 est, par exemple, un spectromètre X et/ou y. [0034] De façon classique, le spectromètre 2 comporte une chaîne de détection 6 et un analyseur 8 relié en sortie de la chaîne de détection 6. [0035] La chaîne de détection 6 est configurée pour générer un signal de détection repré- sentatif de la détection de particules rayonnées par la source 5. En outre, l'analyseur 8 est configure pour analyser le signal de détection généré par la chaîne de détection 6. [0036] La chaîne de détection 6 comporte un détecteur 10 et un organe d'acquisition 12 relié en sortie du détecteur 10. [0037] Le détecteur 10 est configure pour détecter les particules du rayonnement issu de la source 5, et pour délivrer un signal électrique représentatif de ladite détection. Dans le cas où la source 5 émet des photons y, le détecteur 10 est, par exemple, une diode HPGe (de l'anglais « high-purity germanium », signifiant « germanium de haute pureté »). [0038] L'organe d'acquisition 12 est configure pour générer le signal de détection à partir du signal électrique délivré par le détecteur 10. Par exemple, l'organe d'acquisition 12 est configure pour amplifier et/ou filtrer et/ou numériser le signal électrique délivré par le détecteur 10 afin de générer le signal de détection. [0039] Le détecteur 10 et l'organe d'acquisition 12 sont connus et ne seront pas décrits davantage. [0040] L'analyseur 8 est configure pour analyser le signal de détection délivré par la chaîne de détection 6, pour établir un spectre du rayonnement émis par la source 5, ou encore pour mesurer l'activité radiologique de la source 5. [0041] L'analyseur 8 est associé à une plage d'énergie de détection prédéterminée, subdivisée en un nombre entier prédéterminé n, de canaux de détection. [0042] Chaque canal de détection est repéré par un rang noté x. Par la suite, l'expression canal x » désignera un canal de rang x donné. [0043] En outre, chaque canal de détection x de l'analyseur 8 est associé à un intervalle d'énergie correspondant de la plage d'énergie de détection, noté E,. Par la suite, l'expression « intervalle d'énergie E' » désignera l'intervalle d'énergie associé au canal x. [0044] Au sens de la présente invention, le spectre est défini comme un histogramme de la population de particules détectées, en fonction de l'énergie déposée dans le détecteur 10 par lesdites particules. Plus précisément, à chaque canal x est associé le nombre de particules qui ont déposé, dans le détecteur 10, une énergie comprise dans l'intervalle 7 d'énergie L. [0045] Un tel analyseur 8 est classiquement connu et ne sera pas d'écrit davantage. [0046] Da façon classique, le fonctionnement du spectromètre 2 est susceptible d'être décrit par un modèle physique, un tel modèle physique du spectromètre comportant np paramètres prédéterminés. Le modèle physique du spectromètre est, notamment, destine à être utilise pour déterminer l'activité radiologique d'une source quelconque en fonction d'une réponse du spectromètre 2 lors de l'exposition du détecteur 10 du spectromètre 10 à ladite source. [0047] Un tel modèle physique est, comme cela sera d'écrit ultérieurement, stocké dans le système d'étalonnage 4. [0048] Le système d'étalonnage 4 est configure pour déterminer, durant un étalonnage du spectromètre 2, une valeur, dite « valeur étalonnée », de chacun parmi les np paramètres prédéterminés du modèle physique prédéterminé du spectromètre 2. [0049] Chaque paramètre est repéré par un rang note p. Par la suite, l'expression « paramètre p » désignera un paramètre prédéterminé de rang p donné. [0050] Le système d'étalonnage 4 comprend une mémoire 14, un processeur 16 relie à la mémoire 14 et, par exemple, une interface homme/machine 18. [0051] La mémoire 14 comporte un emplacement mémoire de configuration 20, un em- placement mémoire de traitement 22 et un logiciel d'étalonnage 24. [0052] L'emplacement mémoire de configuration 20 est configure pour stocker des données relatives au spectromètre 2. [0053] En particulier, l'emplacement mémoire de configuration 20 est configure pour stocker le modèle physique du spectromètre 2. Comme énoncé précédemment, le modèle physique du spectromètre 2 comprend ni, paramètres. De tels paramètres sont, par exemple, les dimensions du détecteur 10, les masses volumiques des matériaux constituant le spectromètre 2 (en particulier, constituant le détecteur 10), ou encore les concentrations des différents isotopes présents dans lesdits matériaux. [0054] L'emplacement mémoire de configuration 20 est également configure pour stocker l'intervalle d'énergie E. associé à chacun des n, canaux de détection. [0055] L'emplacement mémoire de configuration 20 est, en outre, configure pour stocker des informations relatives à des configurations expérimentales mises en oeuvre pour l'étalonnage du spectromètre 2, et pour chacune desquelles le détecteur 10 est exposé à au moins une source de rayonnements ionisants. Les configurations expérimentales sont deux à deux distinctes. [0056] De telles informations comportent, par exemple, un nombre n, de configurations ex- périmentales mises en oeuvre pour l'étalonnage du spectromètre 2. [0057] Chaque configuration expérimentale est repérée par un rang noté z. Par la suite, l'expression « configuration expérimentale z » désignera une configuration expé- 8 ri mentale de rang z donné. [0058] En outre, de telles informations relatives aux configurations expérimentales mises en oeuvre pour l'étalonnage du spectromètre 2 et stockées dans l'emplacement mémoire de configuration 20 comportent, par exemple, la nature de la ou des source(s) utilisée(s) lors de chaque configuration expérimentale z, ainsi que la position de la ou des source(s) par rapport au spectromètre 2. [0059] En outre, l'emplacement mémoire de configuration 20 est configuré pour stocker, en relation avec chacune des n, configurations expérimentales, le signal de détection correspondant obtenu au moyen du spectromètre 2. [0060] L'emplacement mémoire de configuration 20 est, en outre, configuré pour stocker une première valeur prédéterminée strictement positive, appelée « premier seuil de significativité ». [0061] L'emplacement mémoire de configuration 20 est également configuré pour stocker une deuxième valeur prédéterminée strictement négative, appelée « deuxième seuil de significativité », par exemple égale à l'opposé du premier seuil de significativité. [0062] En outre, l'emplacement mémoire de configuration 20 est configuré pour stocker, pour chacun des n' paramètres du modèle physique du spectromètre 2, une loi de probabilité correspondante. Une telle loi de probabilité est, par exemple, une loi normale. [0063] L'emplacement mémoire de configuration 20 est également configuré pour stocker un premier coefficient réel strictement supérieur à 1, et un deuxième coefficient réel strictement supérieur à un I. [0064] L'emplacement mémoire de traitement 22 est configuré pour stocker, pour chacun des n,, paramètres du modèle physique du spectromètre 2, les valeurs de grandeurs caractéristiques de la loi de probabilité correspondante. A titre d'exemple, dans le cas où ladite loi de probabilité est une loi normale, de telles grandeurs caractéristiques sont son espérance et son écart-type. [0065] Par la suite, l'expression « loi de probabilité courante » désignera la donnée d'une loi de probabilité et de la valeur courante de ses grandeurs caractéristiques. [0066] Comme cela sera décrit ultérieurement, la valeur des grandeurs caractéristiques de chaque loi de probabilité est mise à jour au cours du temps, durant la mise en oeuvre du procédé d'étalonnage objet de l'invention. [0067] Le logiciel d'étalonnage 24 est configuré pour calculer la valeur étalonnée de chacun des n,, paramètres du modèle physique du spectromètre 2. Pour ce faire, le logiciel d'étalonnage 24 est configuré mettre en oeuvre une étape 25 de détermination de rendement de mesure expérimental, une étape de calcul 26, une étape de vérification 27 et une étape de choix 28. [0068] Le logiciel d'étalonnage 24 est configuré de façon à déterminer, au cours de l'étape de détermination 25, pour chaque configuration expérimentale z, un rendement de 9 mesure expérimental correspondant à chaque canal x, et ce à partir du signal de détection fourni par le spectromètre 2 dans ladite configuration expérimentale z. [0069] Plus précisément, le logiciel d'étalonnage 24 est configure pour déterminer le rendement de mesure expérimental, noté fr,P, par la mise en oeuvre de la relation (1) suivante : [0070] [Math.1] eNp S, f = (1) z A [0071] où fe.m) est le rendement de mesure expérimental associé à la configuration expérimentale z, pour le canal x; [0072] S', est le nombre d'impulsions (obtenu après déduction du bruit de fond) cor- respondant au dépôt, dans le détecteur 10, d'une énergie appartenant à l'intervalle d'énergie Et, pour la configuration expérimentale z ; [0073] A, est l'activité, connue, de la ou des source(s) dans la configuration expérimentale z; [0074] L est le nombre de photons d'énergie E, émis par désintégration du ou des radio- isotopc(s), accessible dans les bases dc données nucléaires ; [0075] Ta est un temps actif de la mesure, défini comme la différence entre un temps total de la mesure et un temps mort de la mesure ; et [0076] T, est un temps réel de la mesure, défini comme le temps total de la mesure. [0077] L'écart-type correspondant au rendement de mesure expérimental Pm'. noté G. est donné par la relation (2) suivante : [0078] [Math.2] G( rxzex,P ) AIr 'D (2) a [0079] Le logiciel d'étalonnage 24 est également configure pour mettre en oeuvre, après l'étape de détermination 25, l'étape de calcul 26. Plus précisément, le logiciel d'étalonnage 24 est configure pour mettre en oeuvre l'étape de calcul 26 de façon itérative, c'est-à-dire à plusieurs reprises, successivement dans le temps, jusqu'à ce qu'une condition prédéterminée soit remplie. [0080] L'étape de calcul 26 comporte une phase 30 d'élaboration de modèles, une phase 32 de détermination de rendement de mesure simulé, une phase 34 de calcul de critère de mérite et une phase d'adaptation 36. [0081] La description des phases 30, 32, 34 et 36 va être faite en référence à une itération donnée quelconque, notée k, de l'étape de calcul 26. [0082] Le logiciel d'étalonnage 24 est configuré pour générer,au cours de la phase d'élaboration 30, n', modèles du spectromètre 2. [0083] Chaque modèle du spectromètre est repéré par un rang noté ni. Par la suite, l'expression « modèle ni » désignera un modèle de rang m donné du spectromètre 2. [0084] Plus précisément, le logiciel 24 est configuré pour, au cours de la phase d'élaboration 30, générer un modèle ni quelconque à partir du modèle physique du spectromètre 2 stocké dans l'emplacement mémoire de configuration 20, en choisissant, pour chacun des np paramètres prédéterminés du modèle m, une valeur respective dépendant de la loi de probabilité courante correspondante. [0085] En d'autres ternies, pour une configuration expérimentale z donnée, et pour un modèle ni donné, la valeur d'un paramètre p est obtenue par un tirage mettant en oeuvre la loi de probabilité courante associée. [0086] En outre, le logiciel d'étalonnage 24 est configuré pour déterminer, après la phase d'élaboration 30, et au cours de la phase 32 de détermination de rendement de mesure simulé, pour chaque configuration expérimentale z, et pour chaque modèle ni élaboré (c'est-à-dire généré), un rendement de mesure simulé de chacun des n, canaux de détection. Le rendement de mesure simulé est noté fuui [0087] Plus précisément, pour chaque configuration expérimentale z et pour chaque modèle ni, le rendement de mesure simulé Pm est obtenu au moyen dudit modèle m, pour ladite configuration z. [0088] Le logiciel d'étalonnage 24 est également configuré pour calculer, après la phase 32 de détermination de rendement de mesure simulé, et au cours de la phase de calcul 34, un critère de mérite associé à chaque modèle m généré durant la phase d'élaboration 30. [0089] Plus précisément, pour un modèle m quelconque de l'itération k, le logiciel d'étalonnage 24 est configuré pour calculer le critère de mérite correspondant, noté B,, .k. à partir du rendement de mesure simulé th' associé à chacune des n, configurations expérimentales et à chacun des n, canaux de détection, et du rendement de mesure expérimental t'ex') associés à chacune des n, configurations expérimentales, à chacun des n , canaux de détection, et à l'itération courante k. [0090] Par exemple, le critère de mérite B associé à un modèle m quelconque de l'itération k est défini comme un biais, par exemple obtenu au moyen de la formule (3) suivante : [0091] [Math.3] I 2 ( rsim \ ,'sizn,i' 2 2 fexP il a lxink. xr G , z, feXp Y. kfCXI) \ 4 V X, / esim 1) _feK TS fexp x,7 rix (3) pt ) "n -11,,nz 7=- 100921 où 11 f- situ X, zan.,k est le rendement de mesure simulé associé au canal x, dans la configuration expérimentale z et pour le modèle m, à l'itération k courante ; et [0093] 1 etS1111 I X, Z,111,k ) est l'écart-type du rendement de mesure simulé m fsi x, zan,k [0094] En outre, le logiciel d'étalonnage 24 est configure pour, après la phase de calcul 34. et au cours de la phase d'adaptation 36, mettre à jour la valeur courante des grandeurs caractéristiques de la loi de probabilité associée à chacun des np paramètres prédéterminés. En d'autres termes, le logiciel d'étalonnage 24 est configure pour. au cours de la phase d'adaptation 36 de l'itération k, écrire dans l'emplacement mémoire de traitement 22, pour la loi de probabilité associée à chacun des np paramètres prédéterminés, les valeurs des grandeurs caractéristiques correspondantes qui seront à mettre en oeuvre au cours de l'itération suivante k+1. [0095] Plus précisément, le logiciel d'étalonnage 24 est configure pour, au cours de la phase d'adaptation 36, mettre à jour la valeur des grandeurs caractéristiques de la loi de probabilité associée à chaque paramètre p à partir de l'analyse de l'évolution de la valeur du critère de mérite en fonction des n' valeurs tirées pour le paramètre p au cours de la phase d'élaboration 30. [0096] Pour ce faire, le logiciel d'étalonnage 24 est, de préférence, configure pour, au cours de la phase d'adaptation 36, associer la valeur du paramètre p et la valeur du critère de mérite B11,k obtenue pour chaque modèle m. En d'autres termes, pour chaque paramètre p, le logiciel d'étalonnage 24 est configure de façon à associer la valeur du paramètre p tirée pour chaque modèle m et la valeur du critère de mérite correspondant audit modèle m. Il en résulte, pour chaque paramètre p, np, couples associant la valeur du paramètre p tirée pour chaque modèle m généré et la valeur du critère de mérite correspondant audit modèle na. [0097] Dans ce cas, le logiciel d'étalonnage 24 est également configure pour calculer, pour chaque paramètre p, et à partir des n, valeurs prises par ledit paramètre p et des n', critères de mérite correspondants, une approximation de l'évolution du critère de mérite Bk en fonction de la valeur du paramètre p. [0098] Une telle approximation est, par exemple, fondée sur une fonction d'approximation quadratique, notée y,' de la forme y1, = ap.p2 + b,'p + c, où a, b1, et c1, sont des réels associés, pour l'itération k, au paramètre p et qui sont déterminés par le logiciel d'étalonnage 24 durant la modélisation décrite ci-dessus, par exemple au moyen de la 12 méthode des moindres carrés. [0099] En outre, dans ce cas, le logiciel d'étalonnage 24 est configuré pour modifier la valeur courant des grandeurs caractéristiques de la loi de probabilité correspondant à chaque paramètre p à partir du résultat de la modélisation. [0100] Pour ce faire, le logiciel d'étalonnage 24 est, par exemple, configuré pour calculer, pour chaque paramètre p, un facteur de forme, noté V', représentatif de la convexité de la fonction d'approximation y,. [0101] Par exemple, pour un paramètre p donné quelconque, le logiciel d'étalonnage 24 est configuré pour calculer le facteur de forme V, correspondant par la mise en oeuvre de la relation (4) suivante : [0102] [Math.4] a var (y ) Pl