Brevet FR3141554A1 — publié le 20240503

Procédé d’estimation d’un déséquilibre axial de puissance dans un réacteur nucléaire (11), comprenant les étapes suivantes : - obtention d’une consigne de puissance du réacteur (12), - pour chaque variable d’une pluralité de variables du réacteur, détermination d’une séquence de la variable (19), la séquence représentant des variations futures de la variable estimées, - détermination d’une séquence du déséquilibre axial de puissance (13), en prenant en compte les séquences de la pluralité de variables du réacteur (19), la détermination de la séquence du déséquilibre axial de puissance (13) utilisant un module d’apprentissage automatique (14) préalablement entraîné sur des données historiques de réacteurs (17). Figure pour l’abrégé : Fig. 2 Procédé d’estimation d’une valeur future du déséquilibre axial de puissance dans un réacteur nucléaire DOMAINE DE L'INVENTION L’invention concerne l’estimation d’une valeur future du déséquilibre axial de puissance dans un réacteur nucléaire. ETAT DE LA TECHNIQUE Les centrales nucléaires, en tant que moyens de production «pilotables», représentent l’un des contributeurs principaux à l’équilibre du réseau électrique en France. L’électricité ne pouvant pas (ou peu) être stockée, l’égalité entre la puissance électrique produite par les différents moyens de production et la puissance électrique consommée par les usagers doit être vérifiée à chaque instant. Or, les consommateurs et une partie des moyens de production, en particulier les moyens de production d’électricité dits renouvelables telles que l’éolien et le photovoltaïque, sont non pilotables et intermittents : un moyen de production pilotable, tel que le nucléaire, est donc nécessaire pour compenser les fluctuations de la consommation et de la production non pilotable pour assurer l’équilibre offre/demande Pour cela, les réacteurs nucléaires des centrales nucléaires doivent être techniquement en mesure de faire varier rapidement la puissance électrique produite en fonction du besoin, tout en restant dans les domaines de fonctionnement et de sûreté qui leurs sont assignés. Les domaines de fonctionnement sont déterminés notamment en fonction du déséquilibre axial de puissance et des oscillations Xénon présentées maintenant. Déséquilibre axial de puissance dans un réacteur Au sein d’un réacteur l’énergie thermique destinée à être transformée en énergie électrique est produite par fission. La fission nucléaire est un phénomène par lequel un noyau atomique lourd tel que l'Uranium ou le Plutonium se scinde en deux (fission binaire) ou trois (fission ternaire) nucléides plus légers soit de manière spontanée soit après avoir absorbé un neutron. On parle de neutronique pour désigner l'étude du déplacement des neutrons dans la matière et des réactions qu'ils y induisent, en particulier la génération de puissance par fission. La neutronique est à la base de la conception des cœurs de réacteurs nucléaires à fission contrôlée, tels que les réacteurs à eau pressurisée. En rapport avec la figure 1 illustrant un réacteur 1, on peut définir la direction axiale z comme la direction verticale donnée par un fil à plomb. Cette direction z est orientée verticalement vers le haut dans la figure 1. On s’intéresse à caractériser un écart entre une puissance dans la moitié haute du cœur du réacteur et une puissance dans la moitié basse du cœur du réacteur. La courbe 5 représente la puissance en fonction de la hauteur z dans le réacteur. Cet écart varie sous l'effet des barres de contrôle ou grappes, de la présence de poisons neutroniques tels que des produits de fission comme le xénon ou le samarium, de la température du combustible ou encore de son taux d'épuisement qui croît avec le temps au cours de la campagne de fonctionnement du réacteur. L’écart de puissance entre le haut et le bas du réacteur peut être caractérisé par deux grandeurs physiques : un décalage axial (que l’on peut désigner par sa dénomination anglaise « Axial Offset » abrégée en AO) et un déséquilibre axial de puissance que l’on peut désigner par les notations ou « DPAX »). Ces deux grandeurs sont définies par les formules suivantes : avec la puissance nominale du réacteur. Le terme désigne la puissance totale instantanée dans le réacteur. Le décalage axial est une valeur relative de puissance par rapport à la puissance totale. Le déséquilibre axial « DPAX » est une valeur relative de puissance par rapport à la puissance nominale , puissance nominale qui n’est pas nécessairement égale à la puissance totale instantanée P t . Oscillation Xénon Le Xénon 135 ( 135 Xe) est un élément fortement absorbeur de neutrons pouvant être produit par fission mais aussi et surtout par décroissance de l'Iode 135 ( 135 I), un autre produit de la réaction de fission. La présence du Xénon 135 au sein d'un réacteur entraîne une baisse importante de réactivité, faisant de lui un poison neutronique qui nuit à l'entretien de la réaction en chaîne. Après un temps de fonctionnement du réacteur suffisamment long (environ 60 heures) à puissance constante du réacteur, la concentration en Xénon 135 atteint un premier équilibre dans le réacteur. Cet équilibre est atteint lorsque la production de Xénon 135 devient proportionnelle à la puissance du réacteur et que sa destruction devient proportionnelle à la concentration en Iode 135. A l'échelle macroscopique, si l'on diminue la puissance du réacteur une fois cet équilibre atteint, la concentration en Iode 135 décroît exponentiellement jusqu'à atteindre une nouvelle valeur d'équilibre tandis que la concentration en Xénon 135 passe d'abord par un maximum (on parle de pic xénon) avant de décroître jusqu'à atteindre une nouvelle valeur d'équilibre. Le temps nécessaire pour atteindre ce pic (entre 7 et 8 heures) ainsi que la valeur d'antiréactivité qui lui correspond dépendent de la puissance du réacteur avant et après la baisse. La réactivité ρ est une grandeur sans unité exprimée en pcm (partie pour cent mille) qui se définit comme une variation relative du coefficient effectif de multiplication de la population neutronique k eff lui-même défini comme le rapport de la population de neutrons de la présente génération et de la population de neutrons de la génération précédente : Le signe de la réactivité ρ indique si la population de neutrons est stable (ρ= 0), en augmentation (ρ > 0), ou en diminution (ρ On parlera d’antiréactivité lorsque ρ Inversement, si l'on augmente la puissance, la concentration en Iode 135 croît jusqu'à atteindre une nouvelle valeur d'équilibre tandis que la concentration en Xénon 135 passe d'abord par un minimum (on parle de dépression Xénon) avant d'augmenter jusqu'à atteindre une nouvelle valeur d'équilibre. Le temps nécessaire pour atteindre ce minimum (entre 2 et 3 heures) ainsi que la valeur d'antiréactivité qui lui correspond vont dépendre du niveau de puissance avant et après la hausse. Localement, une variation de puissance, par insertion des grappes, se traduit par un effet axial sur la distribution du flux de neutrons et entraîne une modification de la distribution de l'Iode 135 et du Xénon 135. La illustre l’insertion puis le retrait d’une barre de poison neutronique ou grappe 3 s’étendant selon la direction verticale z. L’effet de cet barre 3 insérée puis retirée peut être décomposé en cinq phases correspondants aux lettres A à E en . Cette manœuvre impacte notamment la concentration en Xénon 9 en fonction de la hauteur z dans le réacteur, et la concentration en Iode 7 en fonction de la hauteur z dans le réacteur L’insertion puis le retrait de la barre 3 peut conduire à des oscillations de la concentration en Xénon 9 entre la moitié haute et la moitié basse du réacteur communément appelées oscillations xénon. Cette oscillation se décompose en cinq phases correspondants aux lettres A à E représentées en figures 1. La phase A correspond à un régime permanent dans lequel la distribution axiale de puissance 3 et les concentrations 7 et 9 en Iode 135 et Xénon 135 sont en phase. La phase B correspond à l'insertion d'un groupe de grappes 3 dans le milieu réactionnel du réacteur 1. Cela induit une modification de la distribution axiale de puissance : la puissance baisse dans la moitié haute mais reste à peu près constante dans la moitié basse. L'équilibre entre la production, par décroissance radioactive de l'Iode 135 et la disparition de Xénon 135 par capture d'un neutron ou par décroissance radioactive est rompu. Dans les premiers instants de ce type de transitoires, les distributions en Iode 135 et Xénon 135 restent inchangées tandis que la puissance dans le bas du cœur augmente et celle en haut diminue. Durant la phase C on observe dans le haut du cœur, du fait de la baisse de flux de neutrons provoquée par l’insertion du groupe de grappes, la concentration 7 en Iode 135 qui diminue et le Xénon 135 qui disparaît plus lentement et commence à s'accumuler : la concentration 9 en Xénon 135 augmente provoquant alors une baisse de réactivité. Dans le bas du cœur, du fait de la hausse de flux, la concentration 7 en Iode 135 augmente et la concentration 9 en Xénon 135 diminue car le Xénon 135disparaît plus rapidement. La réactivité augmente dans la moitié basse. Cette redistribution de la réactivité entre le haut et le bas du cœur tend à amplifier le phénomène de déséquilibre avec un maximum du déséquilibre axial atteint au bout de 7 à 8 heures et ce, même si l’on extrait le groupe de barres inséré La phase D correspond à la période située après que le maximum du déséquilibre axial ait été atteint. Dans le haut du cœur, la diminution de concentration 7 en Iode 135 entraîne un déficit de production du Xénon 135 dont la concentration 9 commence à décroître, provoquant ainsi une augmentation de la réactivité. Dans le bas du cœur, l'accumulation de l'Iode 135 compense la disparition du Xénon 135 dont la concentration 9 commence à augmenter, entraînant ainsi une baisse de la réactivité. Ce processus se poursuit au cours du temps jusqu'à obtenir une répartition spatiale de la puissance proche de la situation initiale au bout de 17 à 18 heures. À ce moment-là, la concentration 7 en Iode 135 et la concentration 9 Xénon 135 sont, elles, fortement déséquilibrées. La courbe du Xénon en phase D se trouve dans une position intermédiaire entre la courbe du Xénon en phase C et la courbe du Xénon en phase E. En phase E, le phénomène s'inverse et le pic de puissance bascule vers le haut du cœur. La période de cette oscillation est de l'ordre de 35 heures. Les déformations de flux de neutron induites par le mouvement des grappes de contrôle et/ou par les oscillations Xénon peuvent être à l’origine de la formation de points chauds à l’intérieur du réacteur conduisant à une élévation excessive de température du combustible pouvant lui faire perdre son intégrité. Pour prévenir ce risque, le déséquilibre axial de puissance doit être obligatoirement contrôlé pour le pilotage des réacteurs. Les opérateurs s'appuient sur un diagramme de pilotage représentant le déséquilibre axial en fonction de la puissance totale instantanée. Ce diagramme illustre en particulier une zone de sécurité ou domaine de fonctionnement du réacteur qui est une zone autour d’une droite de référence. Le réacteur doit fonctionner dans la zone proche de la droite de référence. Pour piloter les réacteurs de manière sûre et efficace, les opérateurs ont donc besoin de disposer d’outils capables de prévoir le déséquilibre axial de puissance DPAX en fonction de la stratégie de pilotage envisagée, leur permettant ainsi d’évaluer cette stratégie en amont et d’en envisager éventuellement une autre qui permettra de limiter plus le déséquilibre axial de puissance. Les solutions actuelles d’estimation de l’état futur du réacteur ne permettent pas une aide au pilotage satisfaisante. Concrètement il est demandé aux opérateurs de modifier des consignes de puissance et de surveiller l’évolution des principaux paramètres clés du cœur du réacteur nucléaire afin de respecter les contraintes de fonctionnement. Comme évoqué, la variation de puissance et les autres variations qui lui sont associées (variations de la température moyenne du cœur, mouvements des barres de contrôle, évolution de la concentration en acide borique dans le réfrigérant du réacteur, etc.) génèrent des perturbations dont l’effet est difficile à anticiper : l’état du réacteur, et donc sa réponse aux sollicitations, évolue sans cesse ; les phénomènes en jeu sont complexes et caractérisés par des temps caractéristiques très hétérogènes, de quelques secondes à plusieurs heures, et les évolutions de certains paramètres clés, tels que la concentration et la répartition spatiale du xénon dans le réacteur, ne sont pas mesurables et donc non connus par l’opérateur. Ainsi, pour que les réacteurs nucléaires des centrales nucléaires puissent être techniquement en mesure de faire varier rapidement la puissance électrique produite en fonction du besoin, et cela plus fréquemment que les pratiques antérieures, tout en restant dans les domaines de fonctionnement et de sûreté qui leurs sont assignés, il existe le besoin de disposer d’un outil d’estimation fiable et rapide pour prévoir les effets futures et retardés de chaque action de l’opérateur, de manière aisée, cela dans un contexte de planning tendu. Un but de l’invention est de proposer un procédé d’estimation d’une valeur future d’un déséquilibre axial de puissance dans un réacteur nucléaire. Le but est atteint dans le cadre de la présente invention grâce à un procédé d’estimation d’une valeur future d’un déséquilibre axial de puissance dans un réacteur nucléaire, le procédé comprenant les étapes suivantes : - obtention d’une séquence de valeurs successives d’une consigne de puissance du réacteur, - pour chaque variable d’une pluralité de variables du réacteur, détermination d’une séquence de valeurs successives de la variable, la séquence représentant des variations futures de la variable, les variations étant estimées en prenant en compte la consigne de puissance, la pluralité de variables du réacteur comprenant une concentration en xénon dans une moitié haute du réacteur et une concentration en xénon dans une moitié basse du réacteur, - détermination d’une séquence de valeurs successives du déséquilibre axial de puissance, en prenant en compte les séquences de la pluralité de variables du réacteur, la détermination de la séquence du déséquilibre axial de puissance utilisant un module d’apprentissage automatique préalablement entraîné sur des données historiques de réacteurs, et la détermination de la séquence de la concentration en xénon dans la moitié haute du réacteur et la détermination de la séquence de la concentration en xénon dans la moitié basse du réacteur utilisant un modèle d’évolution d’une concentration en iode et d’une concentration en xénon en fonction d’un flux de neutrons Un tel procédé est avantageusement et optionnellement complété par les différentes caractéristiques suivantes prises seules ou en combinaison : - pour au moins une variable de la pluralité de variables, la détermination de la séquence temporelle de la variable comprend une mesure de la variable ; - la pluralité de variables du réacteur comprend au moins l’une parmi une position d’un dispositif configuré pour absorber des neutrons dans le réacteur, une température moyenne mesurée de la cuve du réacteur, un débit d’un fluide caloporteur circulant dans le réacteur, une concentration en une espèce chimique dans le fluide caloporteur, l’espèce chimique étant configurée pour absorber des neutrons dans le réacteur, et un taux de combustion d’un combustible nucléaire contenu dans le réacteur ; - la détermination de la séquence du déséquilibre axial de puissance prend également en compte au moins l’une parmi une température moyenne de référence de la cuve, un type de gestion du combustible correspondant à un fractionnement d’un cœur du réacteur lors d’un renouvellement du combustible, un enrichissement de noyaux fissiles du réacteur et un type de noyaux lourds du réacteur ; - la détermination de la séquence de la pluralité de variables du réacteur et la détermination de la séquence du déséquilibre axial de puissance prend en compte un scénario de commande du réacteur correspondant à une séquence de valeurs successives d’au moins une variable contrôlable par un opérateur ; - une étape de détermination d’un score du scénario de commande, le score étant de préférence une valeur d’une grandeur choisie parmi un volume d’effluent produits, un écart moyen au déséquilibre axial de référence et une distance moyenne à des limites d’un domaine de fonctionnement du réacteur ; et - le scénario de commande est un premier scénario de commande, les étapes du procédé étant réalisées une seconde fois en remplaçant le premier scénario par un deuxième scénario de commande du réacteur correspondant à une autre séquence de valeurs successives d’au moins une variable contrôlable, le procédé comprenant de préférence une étape de comparaison des scores du premier scénario et du deuxième scénario. L’invention porte également sur un programme d’ordinateur comprenant des instructions adaptées à la mise en œuvre d’au moins une des étapes du procédé tel que décrit ci-avant lorsque ledit programme est exécuté sur un ordinateur. L’invention porte enfin sur un dispositif d’estimation d’une valeur future d’un déséquilibre axial de puissance dans un réacteur nucléaire, le dispositif comprenant un module de calcul configuré pour mettre en œuvre le procédé tel que décrit ci-avant, le module de calcul étant configuré pour - obtenir une séquence de valeurs successives d’une consigne de puissance du réacteur, - déterminer, pour chaque variable d’une pluralité de variables du réacteur, une séquence de valeurs successives de la variable, la séquence représentant des variations futures de la variable, les variations étant estimées en prenant en compte la consigne de puissance, la pluralité de variables du réacteur comprenant une concentration en xénon dans une moitié haute du réacteur, et une concentration en xénon dans une moitié basse du réacteur, - déterminer une séquence de valeurs successives du déséquilibre axial de puissance en prenant en compte les séquence de la pluralité de variables du réacteur, le module de calcul comprenant un module d’apprentissage automatique préalablement entraîné sur des données historiques de réacteurs, et le module de calcul utilisant un modèle d’évolution d’une concentration en iode et d’une concentration en xénon en fonction d’un flux de neutrons, de sorte à déterminer la séquence de la concentration en xénon dans la moitié haute du réacteur et la séquence de la concentration en xénon dans la moitié basse du réacteur. DESCRIPTION DES FIGURES D'autres caractéristiques et avantages de l'invention ressortiront encore de la description qui suit, laquelle est purement illustrative et non limitative, et doit être lue en regard des dessins annexés sur lesquels : la , déjà évoquée, est une illustration schématique en rapport avec un phénomène d’oscillation Xénon ; la est une illustration schématique en rapport avec un mode de réalisation d’un procédé d’estimation d’un déséquilibre axial de puissance dans un réacteur; les figures 3 et 6 illustrent schématiquement un déséquilibre axial de puissance en fonction du temps, les figures 4 et 7 illustrent schématiquement une consigne de puissance en fonction du temps, et les figures 5 et 8 illustrent schématiquement une consigne de puissance en fonction d’un déséquilibre axial de puissance. DESCRIPTION DETAILLEE DE L'INVENTION Procédé d’estimation d’une valeur future d’un déséquilibre axial de puissance En rapport avec la , un procédé d’estimation d’une valeur future d’un déséquilibre axial de puissance dans un réacteur nucléaire 11 est proposé. Le déséquilibre axial de puissance dans un réacteur nucléaire a été précédemment introduit sous la notation « DPAX », il s’agit d’une valeur relative de puissance par rapport à la puissance nominale du réacteur. Lors d’une première étape, une séquence 12 de valeurs successives d’une consigne de puissance P du réacteur est obtenue. Une séquence de valeurs successives prises par une grandeur est également désignée ci-après par l’expression « série temporelle ». Cette étape correspond notamment à la situation où un opérateur qui pilote un réacteur nucléaire doit réaliser un transitoire de puissance 12 de ce réacteur. L’opérateur entre un programme de charge c’est-à dire une évolution de la puissance électrique P du réacteur en fonction du temps t. Cette évolution peut notamment être donnée sous la forme d’une série temporelle de valeurs d’une puissance à respecter, c’est-à-dire de valeurs d’une consigne de puissance, en fonction du temps. Une consigne de puissance est donc ici entendue comme une grandeur physique homogène à une puissance et qui est la puissance à respecter. Lors d’une deuxième étape 18 du procédé, le programme de charge est traité pour déterminer des variations futures 19 d’une pluralité de variables du réacteur. Ces différentes variables sont appelées variables explicatives au sens où elles sont utilisées pour estimer les valeurs d’une autre variable qui est le déséquilibre axial de puissance. Il est donc déterminé des valeurs hypothétiques que prendra chacune des variables explicatives. Autrement dit, on détermine pour chaque variable d’une pluralité de variables du réacteur, une série temporelle de la variable représentant des variations futures de la variable estimées en prenant en compte la consigne de puissance. La série temporelle peut également comprendre des variations passées de la variable. Parmi les variables explicatives, il y a la concentration en Xénon dans la moitié haute du réacteur et la concentration en Xénon dans la moitié basse du réacteur. Ces deux variables sont calculées à l’aide d’un modèle d’évolution d’une concentration en Iode et d’une concentration en Xénon en fonction d’un flux de neutrons. Ce modèle se fonde en particulier sur les équations d’évolution de la concentration en Iode et en Xénon suivantes déclinées pour chaque moitié, haute et basse, du réacteur. Dans ces équations, désigne le flux de neutrons, C I désigne la concentration en Iode, C Xe désigne la concentration en Xénon, désigne le rendement de fission de l’Iode et le rendement de fission du Xénon respectivement. Un rendement de fission d’un produit correspond à la probabilité ou taux de production de ce produit lorsque la fission se produit. désigne la section efficace macroscopique de fission, c’est-à-dire la probabilité macroscopique de fission des noyaux lourds reliée à la probabilité par unité de temps qu’un neutron rencontre un noyau lourd et qu’une fission se produise. désigne le taux de décroissance radioactive de l’Iode et le taux de décroissance radioactive du Xénon. Enfin, désigne la section efficace microscopique d’absorption de neutrons du Xénon reliée à la probabilité par unité de temps qu’un neutron rencontre un noyau de Xénon 135 et qu’il soit capturé (ce qui ne génère pas de fission). est un terme représentant la production en iode 135. est un terme représentant la disparition en iode 135. est un terme représentant la production en xénon 135. est un terme représentant la disparition en xénon 135. Lors d’une troisième étape du procédé, une série temporelle 13 du déséquilibre axial de puissance DPAX est déterminée en prenant en compte les séries temporelles de la pluralité de variables du réacteur. La série temporelle représente en fonction du temps t des variations futures du déséquilibre axial de puissance estimées en prenant en compte les séries temporelles de la pluralité de variables du réacteur et le programme de charge. Cette détermination de la série temporelle du déséquilibre axial de puissance est effectuée en utilisant un module d’apprentissage automatique 14 préalablement entraîné sur des données historiques 15 de réacteurs. Il est à noter que le module d’apprentissage peut être réentraîné à la fin d’un cycle de fonctionnement du réacteur pour prendre en compte dans l’apprentissage des données d’exploitation plus récentes. Cet entraînement du modèle sur toute la base prend moins de 4 heures. Les données historiques 15 des réacteurs sont structurées au cours d’une étape d’extraction 16 de sorte à obtenir des données historiques 17 structurées d’exploitation. Notamment, au cours de cette étape d’extraction, les données peuvent être enrichies de valeurs historiques de décalage axial AO et/ou de déséquilibre axial DPAX. On peut chercher en particulier à déterminer le déséquilibre axial DPAX à partir du décalage axial AO lorsqu’il s’agit de la seule variable disponible, les mesures du décalage axial AO étant généralement plus bruitées que celles du déséquilibre axial DPAX, notamment lorsque la puissance instantanée du réacteur est faible. Le module d’apprentissage automatique 14 est configuré pour être entrainé sur des sets de données historiques 15 issues d’un ou plusieurs réacteurs. Ce module d’apprentissage automatique peut notamment comprendre deux sous-modules. Le premier sous-module est configuré pour apprendre une relation entre un ensemble de variables explicatives et le déséquilibre axial de puissance. Après avoir appris cette relation à partir des données historiques 17 d’exploitation des réacteurs, le premier sous-module traite des valeurs 19 des variables explicatives récentes et hypothétiques pour fournir une première prévision du déséquilibre axial de puissance sur les prochaines heures. Plusieurs algorithmes peuvent être utilisés pour mettre ce premier sous-module en œuvre. Dans une première option, on utilise des forêts d’arbres de décision, telles que les forêts aléatoires (également désignées sous l’expression anglaise « Random Forest »), ou les algorithmes de stimulation et les algorithmes de stimulation par gradient (également désignés sous l’expression anglaise « Boosting algorithms » ou « Gradient Boosting algorithms »). Dans une deuxième option, on utilise des modèles additifs généralisés (également désignés sous l’expression anglaise « generalized additive model » ou « GAM »). Dans une troisième option, on utilise des réseaux de neurones. Le deuxième sous-module met en œuvre une correction de l’estimation produite par le premier sous-module. En particulier, le deuxième sous-module peut mettre en œuvre une simple correction d’un biais de prédiction à un instant initial, c’est-à-dire une différence entre le déséquilibre prévu et le déséquilibre mesuré à l’instant initial de la prévision. Le deuxième sous-module peut également mettre en œuvre un modèle plus complexe ayant été entraîné à prédire l’erreur de prédiction du premier sous-module en s’appuyant sur les dernières valeurs prédites et mesurées du déséquilibre axial de puissance. Le deuxième sous-module peut notamment alors utiliser un modèle autorégressif ou processus autorégressif qui est un modèle de régression pour séries temporelles dans lequel la série est expliquée par ses valeurs passées plutôt que par d'autres variables. En variante le deuxième sous-module peut utiliser un réseau de neurones de type Long Short Term Memory (LSTM) qui est une architecture particulière de réseau de neurones. L’utilisation d’un module d’apprentissage automatique, et un modèle d’évolution d’une concentration en Iode dans le réacteur et d’une concentration en Xénon dans le réacteur en fonction d’un flux de neutrons dans le réacteur permet d’estimer de manière fiable et rapide le déséquilibre axial de puissance. Les principaux intérêts de ce procédé sont les suivants. Tout d’abord, le procédé permet une plus grande facilité d’utilisation : le modèle ne nécessite aucun calage au-delà de la phase d’apprentissage qui est réalisée en amont du déploiement. Autrement dit, cette étape n’est pas à la charge de l’opérateur qui pilote le réacteur et plus généralement l’utilisateur du modèle. Ensuite, la précision atteinte dans les estimations du déséquilibre axial est satisfaisante sur une plage temporelle de l’ordre des transitoires de puissance observés, de quelques heures à une douzaine d’heures. L’application du procédé sur des données du parc nucléaire d’EDF – Electricité de France - entre 2008 et 2021 permet de comparer les estimations aux mesures historiques. L’erreur absolue constatées entre les prévisions du déséquilibre axial DPAX faites par l’outil selon la présente invention et les valeurs réelles du déséquilibre axial DPAX sont inférieures à 2% pour 95% des données. Enfin, le procédé permet de réduire le temps de calcul notamment en comparaison au temps de calcul nécessaire aux codes de calculs 3D : une prévision sur plusieurs heures est générée en une fraction de secondes par le procédé contre plusieurs minutes pour un code 3D. La rapidité de calcul de l’outil permet de traiter un très grand nombre de données, par exemple une trentaine de minutes suffisent pour évaluer les prévisions de l’outil sur 11700 transitoires de puissance. Avantageusement, en exploitation, pour l’évaluation de la stratégie à suivre lors d’un transitoire, la durée de traitement est de l’ordre du dixième de seconde. Avantageusement, l’opérateur peut au cours d’un transitoire de puissance qui a fait l’objet d’une première estimation de l’évolution du déséquilibre axial de puissance, relancer un deuxième procédé d’estimation avant la fin de l’horizon choisi pour la première estimation. L’opérateur peut ainsi disposer d’une seconde estimation, ce qui permet de réduire les écarts de prévision. Les figures 3, 4 et 5 d’une part et les figures 6, 7 et 8 d’autre part illustrent une application du procédé. Les figures 3 et 6 sont relatives au déséquilibre axial DPAX en fonction du temps t : les courbes 31 et 61 représentent le déséquilibre axial DPAX mesuré, les courbes 32 et 62 représentent le déséquilibre axial DPAX estimé selon le procédé. L’intervalle temporel sur lequel sont tracées les courbes 31 et 61 est identique à l’intervalle temporel sur lequel sont tracées les courbes 32 et 62. Cet intervalle est gradué en heures de journée, la valeur 00 :00 correspondant à minuit qui sépare la fin d’une première journée à gauche et le début de la journée suivante à droite. Les figures 4 et 7 représentent la puissance thermique fournie par le cœur en fonction du temps t : les courbes 41 et 71 représentent la consigne de puissance. L’intervalle temporel sur lequel est tracée la courbe 41 est identique à l’intervalle temporel sur lequel est tracée la courbe 71. Cet intervalle est gradué en heures de journée, la valeur 00 :00 correspondant à minuit qui sépare la fin d’une première journée à gauche et le début de la journée suivante à droite. Cet intervalle temporel est également identique à l’intervalle temporel sur lequel sont tracées les courbes 31, 32, 61 et 62. Les courbes 41 et 71 montrent un transitoire de puissance avec successivement dans le temps une diminution de la puissance depuis une valeur haute à une valeur basse, suivi d’un palier de puissance à la valeur basse et enfin une remontée de puissance depuis la valeur basse à la valeur haute. Les figures 5 et 8 représentent les évolutions du déséquilibre axial DPAX, mesurées et estimées lorsque la puissance thermique du réacteur suit un transitoire selon un profil programmé. C’est une relation paramétrique entre deux variables fonction chacune du temps. Plus précisément, les tracés de la sont les chemins parcourus par un point (P th (t), DPAX(t)) à partir d’une origine des temps 55. Ces chemins se déduisent des figures 3 et 4 de même origine des temps (34,42). Les courbes des figures 3 et 4 sont représentatives des fonctions P th = P th (t) et DPAX= DPAX(t) ; variables temporelles. De même les tracés de la sont les chemins parcourus par un point (P th (t), DPAX(t)) à partir d’une origine des temps 85, chemins se déduisant des figures 6 et 7 de même origine des temps (63,73) et dont les courbes respectives sont représentatives des fonctions P th = P th (t) et DPAX= DPAX(t) ; variables temporelles. Le tracé 51 en et le tracé 81 en représentent l’évolution du déséquilibre axial DPAX mesuré pendant l’évolution de la puissance thermique selon un transitoire de profil programmé. Les flèches indiquent le parcours des traces 51 et 81 dans le sens chronologique. De même, la courbe 53 en et le tracé 83 en représentent le déséquilibre axial DPAX estimé selon le procédé de la présente invention pendant l’évolution de la puissance thermique selon un transitoire de profil programmé. Les figures (3,4,5) et (6,7,8) correspondent à des représentations dans un diagramme de pilotage à disposition des opérateurs en salle de commande. Le procédé permet d’estimer les valeurs futures du déséquilibre axial DPAX pour un nombre d’heures que l’on peut fixer en amont du procédé, en fonction d’un horizon temporel désiré. Les figures 3 à 8 correspondent à un horizon temporel de 8 heures. C’est pourquoi, sur la , la courbe 32 correspondant au déséquilibre axial DPAX estimé s’arrête à l’instant 33 en sortie du transitoire. L’estimation correspondant à la courbe 32 s’étale sur 8 heures d’environ 21h30 jusqu’à 5h30 le lendemain matin. Il est à noter que le nombre d’heures d’estimation peut être diminué ou augmenté par l’opérateur. L’avancement dans le temps est illustré sur la par le curseur 34, sur la par le curseur 42, sur la par le curseur 63 et sur la par le curseur 73. En figures 3 et 4, les curseurs 34 et 42 sont placés temporellement en amont du transitoire tandis que sur les figures 6 et 7 les curseurs 63 et 73 sont placés temporellement sur le palier de puissance à la valeur basse et en amont du la remontée de puissance depuis la valeur basse à la valeur haute. Lors de la mise en œuvre du procédé, l’opérateur entre le transitoire de puissance en mémoire de sorte à réaliser la première étape du procédé. La séquence de valeurs successives de la consigne de puissance du réacteur obtenue peut être affichée sous la forme d’une graphique tel qu’illustré en figures 4 et 7. Lorsque la troisième étape du procédé est réalisée, la série temporelle du déséquilibre axial de puissance DPAX déterminée peut être également affichée tel qu’illustré par les courbes 32 et 62. Il est à noter que si plusieurs estimations sont successivement lancées, notamment lors d’un transitoire, les différentes courbes du déséquilibre axial de puissance DPAX estimées peuvent être affichées simultanément. L’exemple des figures 3 à 8 correspond à un transitoire réel ayant été effectué dans le passé sur un des réacteurs du parc nucléaire français. Les courbes 32 et 62 d’estimation du déséquilibre axial de puissance DPAX peuvent être comparées aux courbes 31 et 61 du déséquilibre axial de puissance DPAX réellement mesuré au cours de ce transitoire. Dans cet exemple, les valeurs hypothétiques sont en fait les véritables valeurs mesurées de la pluralité de variables au cours du transitoire, c’est-à-dire que l’on a utilisé les données historiques réellement mesurées sur site lors du transitoire en question. En option, le procédé peut comprendre une étape de mesure pour au moins une variable de la pluralité de variables. La détermination de la séquence temporelle de cette variable ou ces variables comprend alors une mesure de la variable ou des variables. La série temporelle comprend alors des variations passées de la variable. Parmi la pluralité de variables du réacteur on peut notamment utiliser les variables suivantes. Une première variable est la position d’un dispositif configuré pour absorber des neutrons dans le réacteur. Plus précisément il s’agit de la position axiale de ce dispositif qui peut être maintenu au-dessus de la zone où circulent les neutrons ou bien descendu dans cette zone. Plus le dispositif est positionné bas dans la structure et plus il capture de neutrons. Un tel dispositif peut comprendre différentes unités ou groupes de barres. Un premier groupe de barres forme un groupe de compensation de puissance (que l’on peut désigner sous l’abréviation GCP). Le GCP a pour but de compenser le défaut de puissance, c’est-à dire l’antiréactivité due à la variation de puissance. Le GCP est composé de quatre groupes de grappes qui s’insèrent dans le coeur respectivement l’un après l’autre selon la suite de séquences suivantes : dans une première séquence, le premier groupe s’insère jusqu’à atteindre une certaine profondeur, ensuite dans une deuxième séquence le deuxième groupe s’insère jusqu’à une certaine profondeur avant que, dans un troisième séquence le troisième groupe ne s’insère jusqu’à une certaine profondeur, et pour finir dans une quatrième séquence le quatrième groupe s’insère jusqu’à une certaine profondeur. Chaque séquence ne peut se produire avant que la séquence précédente ne soit terminée. Ainsi les grappes présentent à chaque instant un recouvrement relatif entre deux grappes distinctes prises deux à deux. Le groupe GCP peut comprendre un groupe de grappes grises qui sont relativement peu absorbantes par rapport à des grappes « noires » définies ci-après. Les grappes grises permettent de limiter l’impact sur le déséquilibre de puissance. Un tel dispositif peut comprendre un deuxième groupe de barres qui est un groupe régulation de température (que l’on peut désigner sous l’abréviation groupe R) de barres absorbantes, noires. Le groupe R est un autre groupe permettant de faire de la régulation en combinaison avec les variations de la concentration en Bore. En rapport avec la première variable, il est possible de déterminer les variations futures de celle-ci de la manière suivante. Pour le groupe de compensation de puissance GCP, la position de ce groupe peut être liée à la puissance thermique du réacteur par une relation. Cette relation est établie lors d’un essai périodique dédié au cours duquel on fait varier la puissance du réacteur afin de déterminer la position des GCP permettant de compenser efficacement les variations de puissance. On peut choisir une situation de décalibrage en imposant que la position effectivement commandée des barres du groupe de compensation de puissance GCP ne soit pas exactement la position donnée par la relation. On peut notamment commander la position des barres du groupe de compensation de puissance GCP par un décalibrage constant. Pour le groupe de régulation de température, ou groupe R, la position des barres de ce groupe R peut être choisie constante. Dans ce cas, la criticité dans le réacteur est obtenue en réglant la concentration en Bore dans le fluide caloporteur. La criticité est en lien avec la population de neutrons dans le réacteur. Lorsque la criticité est atteinte, c’est-à-dire lorsque le régime est critique, la population de neutrons est constante au cours du temps. Atteindre la criticité permet d’entretenir la réaction exothermique dans le réacteur. Il est à noter que la première variable étant contrôlable, les variations futures déterminées peuvent être librement modifiée par l’utilisateur. Une deuxième variable est la température moyenne mesurée de la cuve du réacteur. Pour déterminer les variations futures de cette deuxième variable, la température moyenne mesurée de la cuve est supposée égale à la température de référence. Une troisième variable est la température moyenne de référence de la cuve. La température moyenne de référence correspond à la température qui devrait régner au sein de la cuve compte-tenu de la puissance fournie par le réacteur tandis que la température moyenne mesurée correspond à la température réellement mesurée de la cuve. Pour déterminer les variations futures de cette troisième variable, la température moyenne de référence est déterminée à partir d’une relation fonction de la puissance thermique future du réacteur. Une quatrième variable est un débit du fluide caloporteur circulant dans le réacteur. Le fluide caloporteur est le fluide primaire qui permet la thermalisation des neutrons et la transmission de chaleur depuis le cœur de réacteur vers le circuit secondaire. Pour déterminer les variations futures de cette quatrième variable, le débit primaire peut être établi à partir de la puissance thermique future du réacteur. Plus précisément, une valeur « nominale » du début au début du cycle est définie au cours d’un essai périodique dédié. Ensuite cette valeur « nominale » varie en fonction de la puissance thermique selon une loi linéaire empirique. Une cinquième variable est une concentration en une espèce chimique dans le fluide caloporteur, l’espèce chimique étant configurée pour absorber des neutrons dans le réacteur. Par exemple cette espèce chimique est le Bore. Les valeurs futures de la concentration en l’espèce chimique comme le Bore peuvent être déterminées à l’aide d’un deuxième modèle statistique utilisant en entrée la concentration en Xénon et la position du dispositif configuré pour absorber des neutrons dans le réacteur. Ce deuxième modèle statistique est semblable au premier modèle d’estimation de valeurs futures du déséquilibre axial de puissance, à la différence que, dans le premier modèle, la concentration en bore est une variable explicative et que, dans le deuxième modèle, la concentration en bore est la variable expliquée ou variable de réponse. Lorsque les variations futures de la cinquième variable, et notamment la concentration en bore sont imposées, alors il est nécessaire d’utiliser une autre variable d’ajustement, comme notamment le groupe R. Si la concentration en Bore n’est pas fournie par l’opérateur, le modèle de prévision n’utilise pas le Bore comme variable explicative. Une sixième variable est un taux de combustion d’un combustible nucléaire contenu dans le réacteur. Pour déterminer les variations futures de cette sixième variable, le taux de combustion futur du combustible est considéré comme étant égal au produit de la puissance thermique par le temps de fonctionnement, produit que l’on divise ensuite par la masse initiale de combustible. On peut également utiliser l’intégrale de la puissance thermique sur le temps de fonctionnement, intégrale que l’on divise ensuite par la masse initiale de combustible. En outre, la détermination de la série temporelle du déséquilibre axial de puissance peut également être effectuée en prenant en compte un type de gestion du combustible correspondant à un fractionnement d’un cœur du réacteur lors d’un renouvellement du combustible, un enrichissement de noyaux fissiles du réacteur et un type de noyaux lourds du réacteur Le type de gestion du combustible et le type de noyaux lourds sont des variables catégorielles et restent constantes tout au long d’un cycle. Par exemple, le type de noyaux lourds du réacteur correspond aux catégories uranium enrichi, uranium retraité, plutonium, etc… La détermination de la série temporelle de la pluralité de variables du réacteur et la détermination de la série temporelle du déséquilibre axial de puissance peut aussi prendre en compte un scénario de commande du réacteur correspondant à une séquence d’au moins une variable contrôlable par un opérateur. Une variable contrôlable peut être la position d’un dispositif configuré pour absorber des neutrons dans le réacteur ou encore la concentration en une espèce chimique dans le fluide caloporteur. L’opérateur peut définir un scénario de pilotage, c’est-à-dire les différentes actions qu’il envisage de réaliser. Celles-ci sont relatives aux variables contrôlables comme par exemple les positions des barres de contrôle ou la concentration en Bore. L’opérateur choisit pour au moins une variable contrôlable une séquence de valeurs successives que prend la variable au cours des prochaines heures. Par exemple, l’opérateur peut se fonder sur une première estimation de l’évolution du déséquilibre axial de puissance sans qu’un scénario ne soit précisé et sur la base de cette simulation, il détermine un scénario et relance une deuxième estimation de l’évolution du déséquilibre axial de puissance qui cette fois prend en compte le scénario déterminé. Le procédé peut en outre comprendre une étape de détermination d’un score du scénario de commande. Ce score peut être par exemple une valeur d’une grandeur choisie parmi un volume d’effluent produits, un écart moyen au déséquilibre axial de référence et une distance moyenne ou marge moyenne à des limites d’un domaine de fonctionnement du réacteur. D’autres formules de calcul du score peuvent être utilisées. Le déséquilibre axial de référence correspond à la droite de référence dans le diagramme de pilotage qui représente le déséquilibre axial référence en fonction de la puissance totale instantanée. Le domaine de fonctionnement du réacteur et donc les limites de ce domaine sont également définies en rapport avec le diagramme de pilotage. Ce domaine de fonctionnement est une zone de sécurité qui se situe autour de la droite de référence. Il peut être par exemple défini par des droites de déséquilibres axiaux limites dans le diagramme. La distance moyenne aux limites peut être évaluée comme la distance à l’une de ces droites de déséquilibres axiaux limites. Le score du scénario permet de fournir à l’opérateur une évaluation du scénario estimé par rapport à des critères de production d’effluent qui nécessiteront un traitement particulier ou des critères de stabilité du réacteur. Il est possible de mettre en œuvre une estimation de deux scénarios que l’opérateur compare ensuite sur la base de leurs scores respectifs. Dans ce cas : - le scénario de commande mentionné précédemment est alors un premier scénario de commande, - un deuxième scénario est obtenu et transmis en entrée du procédé, ce deuxième scénario correspond à une séquence d’au moins une variable contrôlable par un opérateur, séquence différente de celle qui correspond au premier scénario, et - les étapes du procédé sont réalisées une seconde fois en remplaçant le premier scénario par le deuxième scénario. Avantageusement, le procédé comprend une étape de comparaison des scores du premier scénario et du deuxième scénario. Il est à noter que les différentes courbes correspondant aux différents scénarios estimés peuvent être affichées sur le même graphe pour pouvoir être comparées visuellement par l’opérateur. Lorsque le procédé comprend une étape de comparaison des scores du premier scénario et du deuxième scénario, le procédé peut en outre fournir un classement des scénarios par ordre croissant pour le critère retenu ou si plusieurs critères ont été calculés chacun de ces critères. En référence à la , c’est le dispositif 20 d’aide à la décision qui peut fournir ce classement des scénarios. Par exemple, en rapport avec le critère du volume d’effluents produits, un premier scénario au cours duquel on estime que 20 m3 d’effluents sont produits est mieux classé qu’un deuxième scénario au cours duquel on estime que 60 m3 d’effluents sont produits. Par exemple, en rapport avec le critère d’écart au déséquilibre axial de référence, un premier scénario au cours duquel on estime que l’écart au déséquilibre axial de référence est en moyenne égal à 2% est mieux classé qu’un deuxième scénario au cours duquel l’écart au déséquilibre axial de référence est en moyenne égal à 4%. Par exemple, en rapport avec le critère de marge aux limites du domaine de fonctionnement et le cas échéant à des droites de déséquilibres axiaux limites, un premier scénario au cours duquel on estime que la marge aux limites du domaine de fonctionnement est en moyenne égal à 5% est mieux classé qu’un deuxième scénario au cours duquel la marge aux limites du domaine de fonctionnement est en moyenne égal à 3%. Le dispositif 20 d’aide à la décision peut également déterminer un score global du scénario qui prend en compte une pluralité de scores. Le dispositif 20 peut alors fournir une suggestion de scénario qui optimise les différents critères pris en compte dans le score global. Une fois que le scénario est choisi par l’opérateur et que les différentes consignes sont réellement imposées au réacteur, il est possible de réaliser des estimations successivement au cours du temps. Ainsi, le procédé peut être répété toutes les minutes de sorte à fournir chaque minute une nouvelle estimation de l’évolution temporelle de l’évolution du déséquilibre axial DPAX. La fréquence de répétition du procédé peut être déterminée à l’avance par l’opérateur. Le procédé peut également être utilisé pour lister différents scénarios ou différentes stratégies de pilotage « acceptables ». Cela consiste à trouver des scénarios, au sens de suites d’actions à réaliser par l’opérateur, acceptables, c’est à dire qu’ils permettent de respecter les limites d’insertion et d’extraction des groupes de contrôle et permettant au cœur de rester critique (keff = 1). L’idée est de tester différentes solutions de pilotage possibles que l’on détermine par exemple en perturbant un ou plusieurs scénarios identifiés préalablement comme acceptables. Le procédé peut également classer ces différents scénarios. Le procédé peut également être utilisé pour évaluer l’impact sur le déséquilibre axial DPAX d’un éventuel besoin d’augmentation rapide en puissance du réacteur. Pour ce faire, on peut déterminer des transitoires supplémentaires de puissance hypothétiques simulant des augmentations rapides de puissance du réacteur. Il s’agit de créer des séquences de valeurs successives de la puissance du réacteur. Ces séquences prennent comme premières valeurs des valeurs de la séquence de consigne de puissance imposée au réacteur. Ensuite, les étapes du procédé d’estimation du déséquilibre axial sont réalisées sur la base de ces transitoires supplémentaires. On peut déterminer un paramètre de stabilité fondé sur l’écart au déséquilibre axial de référence la distance à la droite de référence ou bien sur la marge aux limites du domaine de fonctionnement. Une alerte peut être émise vers l’opérateur lorsque le paramètre de stabilité franchit un seuil, notamment en cas de risque de dépassement de la marge dont l’opérateur dispose par rapport aux limites de fonctionnement à respecter, ou d’une marge réduite. L’opérateur peut ensuite adapter sa stratégie en conséquence et notamment s’il souhaite pouvoir satisfaire un éventuel besoin d’augmentation rapide en puissance du réacteur. Cette estimation peut être effectuée à fréquence constante, par exemple toutes les demi-heures. Enfin, le procédé présenté jusqu’à présent peut permettre d’évaluer un impact potentiel d’un téléréglage sur le déséquilibre axial DPAX. On entend par téléréglage un réglage secondaire de la puissance électrique produite par l’installation et conséquemment de la puissance thermique produite par le réacteur qui a pour double objectif d’une part de restaurer la réserve primaire à sa valeur nominale et d’autre part de rétablir des valeurs de puissance prévues aux points d'interconnexion avec le reste du réseau, par exemple le réseau électrique européen. Une téléréglage est envoyé par un gestionnaire du réseau de transport d’électricité, selon un signal qui peut prendre des valeurs comprises entre -1 et +1. Ce téléréglage est envoyé à un ensemble de centrales qui constituent une réserve secondaire afin qu'elles modulent la puissance qu'elles étaient initialement chargées de produire : la valeur -1 du signal de téléréglage correspond à une réduction de puissance de la totalité d’une puissance de réserve. Une puissance de réserve d’une centrale correspond à environ 5% de la puissance nominale ; la valeur de +1 du signal de téléréglage correspond à une augmentation de la totalité de la puissance de réserve. Pour évaluer l’impact potentiel d’un téléréglage sur un réacteur, le procédé peut estimer une variation de puissance correspondant à un téléréglage selon des étapes proches de celles qui ont été décrites précédemment en rapport avec l’évaluation d’un impact sur le déséquilibre axial DPAX d’un éventuel besoin d’augmentation rapide en puissance du réacteur. Le procédé peut comprendre une étape d’alerte de l’opérateur en cas de risque de dépassement de la marge dont l’opérateur dispose, relativement aux limites de fonctionnement à respecter. L’opérateur peut ensuite adapter sa stratégie en conséquence. Un objet de l’invention est un programme d’ordinateur comprenant des instructions adaptées à la mise en œuvre d’au moins une des étapes du procédé tel qu’on a pu le présenter jusqu’à présent lorsque ledit programme est exécuté sur un ordinateur. Un objet de l’invention est un dispositif d’estimation d’une valeur future d’un déséquilibre axial de puissance dans le réacteur nucléaire, le dispositif comprenant un module de calcul configuré pour mettre en œuvre le procédé tel qu’on a pu le présenter jusqu’à présent, le module de calcul étant configuré pour - obtenir une séquence de valeurs successives d’une consigne de puissance du réacteur, - déterminer, pour chaque variable d’une pluralité de variables du réacteur, une séquence de valeurs successives de la variable, la séquence représentant des variations futures de la variable, les variations étant estimées en prenant en compte la consigne de puissance, la pluralité de variables du réacteur comprenant une concentration en xénon dans une moitié haute du réacteur, et une concentration en xénon dans une moitié basse du réacteur, - déterminer une séquence de valeurs successives du déséquilibre axial de puissance en prenant en compte les séquence de la pluralité de variables du réacteur, le module de calcul comprenant un module d’apprentissage automatique préalablement entraîné sur des données historiques de réacteurs, et le module de calcul utilisant un modèle d’évolution d’une concentration en iode et d’une concentration en xénon en fonction d’un flux de neutrons, de sorte à déterminer la séquence de la concentration en xénon dans la moitié haute du réacteur et la séquence de la concentration en xénon dans la moitié basse du réacteur. Procédé d’estimation d’une valeur future d’un déséquilibre axial de puissance dans un réacteur nucléaire (11), le procédé comprenant les étapes suivantes : - obtention d’une séquence de valeurs successives d’une consigne de puissance du réacteur (12), - pour chaque variable d’une pluralité de variables du réacteur, détermination d’une séquence de valeurs successives de la variable (19), la séquence représentant des variations futures de la variable, les variations étant estimées en prenant en compte la consigne de puissance, la pluralité de variables du réacteur comprenant une concentration en xénon dans une moitié haute du réacteur et une concentration en xénon dans une moitié basse du réacteur, - détermination d’une séquence de valeurs successives du déséquilibre axial de puissance (13), en prenant en compte les séquences de la pluralité de variables du réacteur (19), la détermination de la séquence du déséquilibre axial de puissance (13) utilisant un module d’apprentissage automatique (14) préalablement entraîné sur des données historiques de réacteurs (17), et la détermination de la séquence de la concentration en xénon dans la moitié haute du réacteur et la détermination de la séquence de la concentration en xénon dans la moitié basse du réacteur utilisant un modèle d’évolution d’une concentration en iode et d’une concentration en xénon en fonction d’un flux de neutrons. Procédé selon la revendication 1 dans lequel pour au moins une variable de la pluralité de variables, la détermination de la séquence temporelle de la variable comprend une mesure de la variable. Procédé selon l’une quelconque des revendications 1 à 2 dans lequel la pluralité de variables du réacteur comprend au moins l’une parmi - une position d’un dispositif configuré pour absorber des neutrons dans le réacteur, - une température moyenne mesurée de la cuve du réacteur, - un débit d’un fluide caloporteur circulant dans le réacteur, - une concentration en une espèce chimique dans le fluide caloporteur, l’espèce chimique étant configurée pour absorber des neutrons dans le réacteur, et - un taux de combustion d’un combustible nucléaire contenu dans le réacteur. Procédé selon la revendication 3 dans lequel la détermination de la séquence du déséquilibre axial de puissance prend également en compte au moins l’une parmi une température moyenne de référence de la cuve, un type de gestion du combustible correspondant à un fractionnement d’un cœur du réacteur lors d’un renouvellement du combustible, un enrichissement de noyaux fissiles du réacteur et un type de noyaux lourds du réacteur. Procédé selon l’une quelconque des revendications 1 à 4 dans lequel la détermination de la séquence de la pluralité de variables du réacteur et la détermination de la séquence du déséquilibre axial de puissance prend en compte un scénario de commande du réacteur correspondant à une séquence de valeurs successives d’au moins une variable contrôlable par un opérateur. Procédé selon la revendication 5 comprenant en outre une étape de détermination d’un score du scénario de commande, le score étant de préférence une valeur d’une grandeur choisie parmi un volume d’effluent produits, un écart moyen au déséquilibre axial de référence et une distance moyenne à des limites d’un domaine de fonctionnement du réacteur. Procédé selon l’une quelconque des revendications 5 à 6, dans lequel le scénario de commande est un premier scénario de commande, les étapes du procédé étant réalisées une seconde fois en remplaçant le premier scénario par un deuxième scénario de commande du réacteur correspondant à une autre séquence de valeurs successives d’au moins une variable contrôlable, le procédé comprenant de préférence une étape de comparaison des scores du premier scénario et du deuxième scénario. Programme d’ordinateur comprenant des instructions adaptées à la mise en œuvre d’au moins une des étapes du procédé selon l’une quelconque des revendications 1 à 7 lorsque ledit programme est exécuté sur un ordinateur. Dispositif d’estimation d’une valeur future d’un déséquilibre axial de puissance dans un réacteur nucléaire, le dispositif comprenant un module de calcul configuré pour mettre en œuvre le procédé selon l’une des revendications 1 à 7, le module de calcul étant configuré pour - obtenir une séquence de valeurs successives d’une consigne de puissance du réacteur, - déterminer, pour chaque variable d’une pluralité de variables du réacteur, une séquence de valeurs successives de la variable, la séquence représentant des variations futures de la variable, les variations étant estimées en prenant en compte la consigne de puissance, la pluralité de variables du réacteur comprenant une concentration en xénon dans une moitié haute du réacteur, et une concentration en xénon dans une moitié basse du réacteur, - déterminer une séquence de valeurs successives du déséquilibre axial de puissance en prenant en compte les séquence de la pluralité de variables du réacteur, le module de calcul comprenant un module d’apprentissage automatique préalablement entraîné sur des données historiques de réacteurs, et le module de calcul utilisant un modèle d’évolution d’une concentration en iode et d’une concentration en xénon en fonction d’un flux de neutrons, de sorte à déterminer la séquence de la concentration en xénon dans la moitié haute du réacteur et la séquence de la concentration en xénon dans la moitié basse du réacteur.