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.