FR 3118216 A1 20220624 FR 2013754 A 20201218 L’invention concerne un procédé de gestion coordonnée, sur une période de temps, d’unités de production et de consommation. Le procédé comprend, lors d’une phase de programmation: - recevoir (S1) des simulations de consommation de chaque unité de consommation couvrant la période de temps et paramétrées par une probabilité de transition nominale caractérisant la probabilité de passer d’un état à un autre, - déterminer (S2), en fonction des simulations, une consommation moyenne cible des unités de consommation et un signal de contrôle, - transmettre (S3), à chaque unité de consommation, le signal de contrôle pour mettre à jour (S4) la probabilité de transition nominale sur la base du signal de contrôle reçu et d’une fonction locale caractérisant les spécificités de fonctionnement de l’unité de consommation, et - transmettre (S3), à chaque unité de production, un programme de production cohérent avec la consommation moyenne cible Figure de l’abrégé : Figure 2. Programmation conjointe de flexibilités de production et de consommation Le domaine de l’invention se rapporte à la programmation et à la gestion coordonnées de flexibilités de production et de consommation, notamment le pilotage d’unités de consommation, alimentées en énergie électrique par des unités de production, aptes à moduler leur consommation à la demande d’un opérateur de flexibilité. La flexibilité de consommation est la capacité d’une unité de consommation, par exemple dans le secteur de l’industrie, à moduler sa consommation en énergie électrique à la demande d’un opérateur, appelé opérateur de flexibilité. L’opérateur peut ainsi piloter à distance le fonctionnement d’équipements industriels ou domestiques au niveau d’une ou plusieurs unités de consommation ou retarder un tel fonctionnement. La demande de l’opérateur de flexibilité à une unité de consommation peut être une demande de réduction de sa consommation électrique en période de forte tension sur le système électrique. L’unité de consommation peut alors répondre par l’arrêt du fonctionnement d’équipements industriels pendant quelques minutes voire plusieurs heures. On parle alors d’effacement de consommation électrique. A l’inverse, l’opérateur de flexibilité peut aussi demander une hausse de la consommation électrique pour absorber une production excédentaire d’énergie électrique. La programmation et la gestion des flexibilités de consommation visent ainsi à résorber ou à tirer profit d’un déséquilibre entre l’offre et la demande en énergie électrique au sein d’un système comprenant des unités de consommation alimentées par des unités de production. L’enjeu des flexibilités de consommation est donc de réaliser des économies d’énergie électrique, d’optimiser la performance énergétique et plus généralement de permettre à un tel système de s’adapter aux fluctuations de production et de consommation électriques. Aujourd’hui, les flexibilités de consommation sont parfois ignorées au profit d’autres problématiques liées à la programmation des unités de production, comme la question du dimensionnement ou la prise en compte des contraintes techniques des moyens de production. Lorsque les flexibilités de consommation sont prises en compte, celles-ci sont modélisées par un stockage virtuel agrégeant l’ensemble des flexibilités sans distinction des contraintes individuelles des différents équipements des unités de consommation. Un contrôle automatique des unités de consommation peut également être mis à place, par exemple d’un jour sur l’autre, avec la détermination d’un programme prévisionnel fourni en amont à chaque unité de consommation, laquelle est alors asservie en temps réel à la consigne indiquée dans ce programme prévisionnel. De manière générale, on distingue habituellement entre une approche locale, donc une approche du point de vue d’une unité de consommation gérant en temps réel les différents équipements industriels pour réduire leur consommation électrique, et une approche globale, dans laquelle l’opérateur de flexibilité agrège les flexibilités offertes par l’ensemble des unités de consommation pour optimiser une consommation électrique globale. Toutefois, de telles méthodes ne permettent pas de s’assurer qu’une consommation cible déterminée par l’opérateur de flexibilité est atteignable par l’ensemble d’unités de consommation. De plus, les modèles prédictifs proposant de mettre à jour en temps réel une programmation prévisionnelle et de recalculer la consigne de consommation à chaque événement sont très coûteux en ressources et en temps de calcul. Résumé La présente invention vient améliorer la situation. A ce titre, la présente invention concerne un procédé, mis en œuvre par des moyens informatiques, de gestion coordonnée, sur une période de temps prédéterminée échantillonnée en une succession d’instants, d’un ensemble d’une ou plusieurs unités de production et d’un ensemble d’une ou plusieurs unités de consommation alimentées en énergie électrique par l’ensemble d’une ou plusieurs unités de production. Le procédé comprend, lors d’une phase de programmation: - recevoir au moins une simulation de consommation de chaque unité de consommation, chaque simulation de consommation d’une unité de consommation couvrant la période de temps et étant paramétrée par une probabilité de transition nominale caractérisant la probabilité de l’unité de consommation de passer, à un instant donné de la période de temps, d’un état de fonctionnement à un autre, - déterminer, en fonction de la ou des simulations de consommation reçues, une consommation moyenne cible de l’ensemble d’une ou plusieurs unités de consommation et un signal de contrôle échantillonnés sur la période de temps, - transmettre, à chaque unité de consommation, le signal de contrôle pour mettre à jour, au niveau de chaque unité de consommation, la probabilité de transition nominale sur la base du signal de contrôle reçu et d’une fonction locale caractérisant les spécificités de fonctionnement de l’unité de consommation considérée, et - transmettre, à chaque unité de production, un programme de production sur la période de temps cohérent avec la consommation moyenne cible. Dans un ou plusieurs modes de réalisation, le programme de production comprend au moins, pour chaque unité de production, un profil de production déterminé conjointement à la consommation moyenne cible et au signal de contrôle. Avantageusement, le signal de contrôle échantillonné est déterminé lors de la phase de programmation selon un algorithme itératif défini comme suit : [Math. 1] avec : où : - λ p (t k ) est la valeur à l’instant t k du signal de contrôle à l’itération p, le signal de contrôle étant initialisé à un signal de contrôle initial λ 0 prédéterminé ; - ρ p est un réel positif ; - f 0 est une fonction pénalisant un écart entre un profil de production de l’ensemble d’une ou plusieurs unités de production et un profil de consommation de l’ensemble d’une ou plusieurs unités de consommation ; - T est le nombre d’instants de la succession d’instants de la période de temps prédéterminée ; - est un profil de défaillance à l’itération p, (t k ) est la valeur à l’instant t k du profil de défaillance ; - est un réel caractérisant une contrainte de déviation en énergie nulle de l’ensemble d’une ou plusieurs unités de consommation par rapport à une consommation nominale ; - n est le nombre d’unités de production ; - A j est un ensemble de contraintes techniques constituant l’ensemble des états techniques dans lesquels la j-ème unité de production peut fonctionner ; - p j (x j ) est le profil de production de la j-ème unité de production lorsque celle-ci fonctionne dans un état technique x j , p j (x j )(t k ) étant la valeur à l’instant t k du profil de production p j (x j ) ; - f j (x j ) est un coût de production de la j-ème unité de production lorsque celle-ci fonctionne dans un état technique x j ; - N est le nombre d’unités de consommation ; - M est le nombre de simulations de consommation de chaque unité de consommation ; - κ est une constante prédéterminée ; - exp(·) est la fonction exponentielle ; et - U i,q (t k ) est la valeur à l’instant t k de la consommation de la i-ème unité de consommation selon la q-ième simulation de consommation, l’algorithme itératif étant interrompu lorsqu’un critère d’arrêt prédéfini est vérifié. Le critère d’arrêt est par exemple défini comme suit : où : - est une norme ; et - ε est un réel positif. Typiquement, le signal de contrôle initial permettant d’initialiser le signal de contrôle préalablement à la mise en œuvre de l’algorithme itératif est le signal de contrôle obtenu à l’issue d’une période de temps précédent la période de temps considérée. Avantageusement, le signal de contrôle déterminé à l’issue de l’algorithme itératif est défini comme suit : où p * est le rang de l’itération, parmi toutes les itérations jusqu’à la dernière itération vérifiant le critère d’arrêt, pour laquelle un saut de dualité est minimal, soit : [Math. 8] avec : Avantageusement, la consommation moyenne cible est définie comme suit : [Math. 10] où : - r(t k ) est la valeur à l’instant t k de la consommation moyenne cible de l’ensemble d’une ou plusieurs unités de consommation ; et - p * est le rang de l’itération, parmi toutes les itérations jusqu’à la dernière itération vérifiant le critère d’arrêt, pour laquelle un saut de dualité est minimal, soit : [Math. 11] avec : Dans un ou plusieurs modes de réalisation, la probabilité de transition nominale est mise à jour comme suit lors de la phase de programmation : [Math. 13] avec : [Math. 15] où : - X 1 , X 2 et Y sont des états possibles d’une unité de consommation, dX 1 , dX 2 et dY sont les différentielles de tels états, - P 0 i (t k ,X 1 ,dX 2 ) est la probabilité de transition nominale caractérisant la probabilité de la i-ème unité de consommation de passer, à l’instant t k , d’un état de fonctionnement X 1 à un autre état de fonctionnement X 2 ; - Q 0 i (t k ,X 1 ,dX 2 ) est la probabilité de transition nominale mise à jour lors de la phase initiale de programmation caractérisant la probabilité de la i-ème unité de consommation de passer, à l’instant t k , d’un état de fonctionnement X 2 à un autre état de fonctionnement X 2 ; - V i (t k ,X) est la valeur à l’instant t k de la fonction locale caractérisant les spécificités de fonctionnement de la i-ème unité de consommation et de caractéristiques techniques desdits équipements, lorsque la i-ème unité de consommation fonctionne dans l’état X ; - exp(·) est la fonction exponentielle ; - N est le nombre d’unités de consommation ; - κ est une constante prédéterminée ; - λ * (t k ) est la valeur à l’instant t k du signal de contrôle déterminé ; et - c i étant une fonction qui, à un état de la i-ème unité de consommation, associe un profil de consommation, c i (Y) étant la valeur de la consommation de la i-ème unité de consommation fonctionnant dans l’état Y . La fonction locale dépend par exemple d’une utilisation programmée d’équipements industriels de l’unité de consommation, de caractéristiques techniques desdits équipements industriels et d’aléas locaux. Dans un ou plusieurs modes de réalisation, le procédé comprend en outre, suite à la phase de programmation, une phase de gestion mise en œuvre sur la période de temps comprenant la mise à jour en temps réel du signal de contrôle en fonction de la consommation moyenne cible ; et la mise à jour en temps réel de la probabilité de transition nominale de chaque unité de consommation sur la base du signal de contrôle mis à jour. Avantageusement, dans un tel cas, la probabilité de transition nominale est par exemple mise à jour comme suit lors de la phase de gestion : [Math. 15] où : - est la probabilité de transition nominale mise à jour lors de la phase de gestion caractérisant la probabilité de la i-ème unité de consommation de passer, à l’instant t k , d’un état de fonctionnement X 1 à un autre état de fonctionnement X 2 ; - est la valeur à l’instant t k du signal de contrôle mis à jour. La phase de gestion est par exemple mise en œuvre comme suit : - prévoir, à chaque instant de la période de temps, une consommation moyenne à l’instant suivant de l’ensemble d’une ou plusieurs unités de consommation en fonction de la probabilité de transition nominale mise à jour de chaque unité de consommation, - mettre à jour, à chaque instant, le signal de contrôle sur la base d’une comparaison entre la consommation moyenne à l’instant suivant prévue et de la consommation moyenne cible, et - transmettre, à chaque instant, le signal de contrôle mis à jour à chaque unité de consommation pour mettre à jour, au niveau de chaque unité de consommation, la probabilité de transition nominale sur la base du signal de contrôle mis à jour reçu. Avantageusement, la mise à jour du signal de contrôle lors de la phase de gestion est mise en œuvre comme suit : - construire un modèle de consommation prédictif permettant d’approximer, en fonction du signal de contrôle à un instant donné, la consommation moyenne prévue à l’instant suivant, puis, à chaque instant de la période de temps prédéterminée : - comparer une consommation moyenne réelle et une consommation moyenne prévue à cet instant de l’ensemble d’une ou plusieurs unités de consommation, - adapter par apprentissage, à chaque instant, le modèle de consommation prédictif sur la base de la comparaison, et - mettre à jour, à chaque instant, le signal de contrôle par application réciproque du modèle de consommation prédictif à la consommation moyenne cible visée à l’instant suivant. Avantageusement, au cours de la phase de gestion mise en œuvre sur la période de temps, chaque unité de consommation fonctionne selon la probabilité de transition nominale mise à jour en temps réel, tandis que chaque unité de production suit le programme de production reçu. L’invention concerne en outre un programme informatique, comprenant des instructions pour la mise en œuvre du procédé décrit précédemment lorsque les instructions sont mises en œuvre par au moins un processeur. Enfin, l’invention concerne une unité de traitement pour la gestion coordonnée, sur une période de temps prédéterminée échantillonnée en une succession d’instants, d’un ensemble d’une ou plusieurs unités de production et d’un ensemble d’une ou plusieurs unités de consommation alimentées en énergie électrique par l’ensemble d’une ou plusieurs unités de production. L’unité de traitement est configurée pour, lors d’une phase de programmation: - recevoir au moins une simulation de consommation de l’unité de consommation, chaque simulation de consommation d’une unité de consommation couvrant la période de temps et étant paramétrée par une probabilité de transition nominale caractérisant la probabilité de l’unité de consommation de passer, à un instant donné de la période de temps, d’un état de fonctionnement à un autre, - déterminer, en fonction de la ou des simulations de consommation reçues, une consommation moyenne cible de l’ensemble d’une ou plusieurs unités de consommation et un signal de contrôle échantillonnés sur la période de temps, - transmettre, à chaque unité de consommation, le signal de contrôle pour mettre à jour, au niveau de chaque unité de consommation, la probabilité de transition nominale sur la base du signal de contrôle reçu et d’une fonction locale caractérisant les spécificités de fonctionnement de l’unité de consommation considérée, et - transmettre, à chaque unité de production, un programme de production sur la période de temps cohérent avec la consommation moyenne cible. D’autres caractéristiques, détails et avantages apparaîtront à la lecture de la description détaillée ci-après, et à l’analyse des dessins annexés, sur lesquels : illustre schématiquement un système comprenant un ensemble d’une ou plusieurs unités de production, un ensemble d’une ou plusieurs unités de consommation et une unité de traitement selon l’invention. illustre une phase de programmation du procédé de gestion selon l’invention ; et illustre une phase de gestion en temps réel du procédé de gestion selon l’invention. La illustre schématiquement un système SYS. Le système SYS comprend un ensemble d’une ou plusieurs unités de production UP 1 , UP 2 , UP n , un ensemble d’une ou plusieurs unités de consommation UC 1 , UC 2 , UC N et une unité de traitement UNT. Dans la suite de la description, on considère que le système SYS comprend avantageusement plusieurs unités de production, ici un nombre n d’unités de production UP 1 , UP 2 , UP n . De même, on considère que le système SYS comprend avantageusement plusieurs unités de consommation, par exemple un nombre N d’unités de consommation UC 1 , UC 2 , UC N . Dans le contexte de l’invention, le système SYS peut être assimilé à un portefeuille mixte composé de moyens de production conventionnels, à savoir les unités de production UP 1 , UP 2 , UP n , et de flexibilités de consommation, à savoir les unités de consommation UC 1 , UC 2 , UC N . La gestion d’un tel portefeuille vise à limiter le déséquilibre entre l’offre et la demande en énergie électrique au sein du système SYS dans le but de réaliser des économies d’énergie électrique et de s’adapter aux fluctuations de production et de consommation électriques. Les unités de production UP 1 , UP 2 , UP n sont agencées pour alimenter en énergie électrique les unités de consommation UC 1 , UC 2 , UC N . En d’autres termes, les unités de production UP 1 , UP 2 , UP n sont électriquement connectées aux unités de consommation UC 1 , UC 2 , UC N au sein du système SYS. Chaque unité de production peut être contrôlée en temps réel ou suivre une programmation prédéterminée. Le fonctionnement de chaque unité de production peut être caractérisé par un profil de production et un coût de production. Le profil de production et le coût de production dépendent de l’état technique, aussi appelée état ou mode de fonctionnement, dans lequel l’unité de production fonctionne. On peut ainsi définir, pour chaque unité de production, des contraintes techniques qui constituent l’ensemble des états techniques possibles de sorte que, à un instant donné, l’unité de production considérée fonctionne selon l’un de ces états techniques. Plus spécifiquement, la j-ème unité de production peut être caractérisée par une fonction p j qui, à chaque état technique x j des contraintes techniques A j de la j-ème unité de production, associe un profil de production p j (x j ). Le profil de production p j (x j ) permet de déterminer, à un instant t k donné de la période de temps prédéterminée à venir, la production en énergie électrique p j (x j )(t k ) de la j-ème unité de production fonctionnant dans l’état technique x j . Par ailleurs, il est possible de définir un profil de production de l’ensemble d’une ou plusieurs d’unités de production UP 1 , UP 2 , UP n en fonction du profil de production de chaque unité de production. Le profil de production de l’ensemble d’une ou plusieurs d’unités de production UP 1 , UP 2 , UP n permet de déterminer, à un instant donné de la période de temps, la production totale en énergie électrique de l’ensemble d’une ou plusieurs d’unités de production UP 1 , UP 2 , UP n . Ce profil de production dépend par conséquent des états techniques respectifs des unités de production UP 1 , UP 2 , UP n . Les unités de consommation UC 1 , UC 2 , UC N sont agencées pour être alimentées en énergie électrique par les unités de production UP 1 , UP 2 , UP n . Chaque unité de consommation comprend typiquement un ou plusieurs équipements industriels ou domestiques dont le fonctionnement requiert une alimentation en énergie électrique. De tels équipements peuvent être par exemple un chauffe-eau électrique à accumulation, un réfrigérateur ou encore un système de climatisation. Comme expliqué précédemment, les unités de consommations UC 1 , UC 2 , UC N sont flexibles. En d’autres termes, les unités de consommations UC 1 , UC 2 , UC N sont aptes à moduler leur consommation en énergie électrique respective. Dans le système SYS décrit ici les unités de consommation UC 1 , UC 2 , UC N modulent leur consommation en énergie électrique à la demande de l’unité de traitement UNT, aussi appelée opérateur de flexibilité. Le fonctionnement de chaque unité de consommation est caractérisé par un profil de consommation et une probabilité de transition nominale. Le profil de consommation dépend de l’état de fonctionnement de l’unité de consommation. La probabilité de transition nominale caractérise la probabilité de l’unité de consommation de passer, à un instant donné, d’un état de fonctionnement à un autre. En temps réel, le fonctionnement de chaque unité de consommation peut être modélisé par un comportement aléatoire basé sur la probabilité de transition nominale. Les états de chaque unité de consommation peuvent être supposés mutuellement indépendants et suivent une distribution de probabilité connue. Plus spécifiquement, la i-ème unité de consommation peut être caractérisée par une fonction c i qui, à chaque état ou mode de fonctionnement X i de la i-ème unité de consommation, associe un profil de consommation c i (X i ). Le profil de consommation c i (X i ) permet de déterminer, à un instant t k donné de la période de temps prédéterminée à venir, la consommation en énergie électrique c i (X i )(t k ) de la i-ème unité de consommation fonctionnant dans l’état X i . Par ailleurs, il est possible de définir un profil de consommation de l’ensemble d’une ou plusieurs d’unités de consommation UC 1 , UC 2 , UC N en fonction du profil de consommation de chaque unité de consommation. Le profil de consommation de l’ensemble d’une ou plusieurs d’unités de consommation UC 1 , UC 2 , UC N permet de déterminer, à un instant donné de la période de temps, la consommation totale en énergie électrique de l’ensemble d’une ou plusieurs d’unités de consommation UC 1 , UC 2 , UC N . Ce profil de consommation dépend par conséquent des états ou modes de fonctionnement respectifs des unités de consommation UC 1 , UC 2 , UC N . Les spécificités de fonctionnement de chaque unité de consommation peuvent également être décrites par une fonction locale. Une telle fonction locale dépend par exemple d’une utilisation programmée d’équipements de l’unité de consommation, de caractéristiques techniques de ces équipements et d’aléas locaux. Par ailleurs, comme expliqué dans la suite de la description, une partie au moins de la fonction locale associée à chaque unité de consommation est commune aux fonctions locales respectives des autres unités de consommation. Sur demande de l’unité de traitement UNT, chaque unité de consommation peut réduire sa consommation électrique en période de forte tension sur le système SYS. Cette réduction de la consommation en énergie électrique peut être obtenue par l’arrêt du fonctionnement d’équipements industriels ou domestiques au niveau de l’unité de consommation pendant quelques minutes voire plusieurs heures. On parle alors d’effacement de consommation électrique. Inversement, chaque unité de consommation peut augmenter sa consommation électrique pour absorber une production excédentaire d’énergie électrique. Dans le contexte de la présente invention, on cherche à déterminer, dans le cadre de la gestion des flexibilités de consommation, une programmation des unités de consommation UC 1 , UC 2 , UC N sur une période de temps prédéterminée. Plus précisément, la période de temps prédéterminée est échantillonnée en une succession d’instants de sorte que le fonctionnement de chaque unité de consommation est déterminé pour chaque instant de la période de temps prédéterminée. Typiquement, la période de temps prédéterminée représente plusieurs heures voire un ou plusieurs jours, tandis que le pas de temps est d’une ou plusieurs minutes. Une telle programmation est déterminée en amont de la période de temps prédéterminée et vise à déterminer une consommation moyenne cible des unités de consommation UC 1 , UC 2 , UC N atteignable. Par ailleurs, la présente invention permet aussi une gestion en temps réel, plus exactement d’un instant de la période de temps prédéterminée à l’instant suivant, des flexibilités de consommation pour s’adapter aux contraintes et aléas locaux et s’assurer d’atteindre la consommation moyenne cible calculée en amont. En outre, la présente invention propose la détermination d’un programme de production, destiné aux unités de production UP 1 , UP 2 , UP n cohérent avec les besoins en énergie électrique des unités de consommation UC 1 , UC 2 , UC N . Il est donc proposé une articulation entre la programmation, antérieure à la période de temps prédéterminée, et la gestion en temps réel des flexibilités de consommation. Cette articulation est rendue possible par l’unité de traitement UNT. Au sein du système SYS, l’unité de traitement UNT est un opérateur de flexibilité et a donc pour rôle la gestion des flexibilités de production et de consommation du système SYS. Par conséquent, l’unité de traitement UNT est configurée pour programmer conjointement le fonctionnement, sur une période de temps prédéterminée échantillonnée, des unités de production UP 1 , UP 2 , UP n et des unités de consommation UC 1 , UC 2 , UC N . L’unité de traitement peut être configurée également pour gérer en temps réel, au cours de la période de temps prédéterminée, le fonctionnement des unités de consommation UC 1 , UC 2 , UC N . Plus spécifiquement, dans le cas présent, l’unité de traitement UNT est configurée pour déterminer, lors d’une phase de programmation, un signal de contrôle commun aux unités de consommations UC 1 , UC 2 , UC N ainsi qu’une consommation moyenne cible des unités de consommations UC 1 , UC 2 , UC N . Cette consommation moyenne cible permet par ailleurs de déterminer un programme de production cohérent. L’unité de traitement UNT est également configurée pour transmettre, à chaque unité de consommation, le signal de contrôle afin de mettre à jour, au niveau de chaque unité de consommation, la probabilité de transition nominale sur la base du signal de contrôle reçu et de la fonction locale caractérisant les spécificités de fonctionnement de l’unité de consommation considérée. En particulier, une partie au moins de la fonction locale associée à une unité de consommation dépend du signal de contrôle reçu et est donc commune à tous les unités de consommation UC 1 , UC 2 , UC N . L’unité de traitement UNT est configurée en outre pour transmettre, à chaque unité de production, un programme de production sur la période de temps cohérent avec la consommation moyenne cible. Un tel programme de production permet ainsi aux unités de production UP 1 , UP 2 , UP n d’être coordonnées avec la consommation moyenne cible des unités de consommation UC 1 , UC 2 , UC N sur la période de temps à venir et d’adapter leur production en conséquence pour assurer un équilibre entre l’offre et la demande en énergie électrique. Comme expliqué précédemment, le fonctionnement d’une unité de consommation est caractérisé par sa probabilité de transition nominale de passer, à un instant donné, d’un état de fonctionnement à un autre. Ainsi, la détermination puis la transmission, aux unités de consommation UC 1 , UC 2 , UC N , du signal de contrôle permet une programmation du fonctionnement de chaque unité de consommation. Typiquement, la période de temps sur laquelle le fonctionnement de chaque unité de consommation est programmé correspond à une journée et chaque programmation est donc réalisée d’un jour sur l’autre. Par ailleurs, l’unité de traitement UNT est configurée en outre pour, lors d’une phase de gestion en temps réel mise en œuvre sur la période de temps prédéterminée, mettre à jour en temps réel le signal de contrôle en fonction de la consommation moyenne cible. Ainsi, au cours de cette phase de gestion en temps réel, l’unité de traitement UNT est en outre configurée pour mettre à jour en temps réel la probabilité de transition nominale de chaque unité de consommation sur la base du signal de contrôle mis à jour. Comme illustrée sur la , l’unité de traitement UNT comprend une mémoire MEM, un processeur PROC et un module de communication COM. La mémoire MEM est agencée pour stocker les instructions d’un programme informatique dont la mise en œuvre par le processeur PROC se traduit par le fonctionnement de l’unité de traitement UNT. Le processeur PROC permet ainsi, lors de la phase de programmation, de déterminer le signal de contrôle et la consommation moyenne cible de l’ensemble d’une ou plusieurs unités de consommation UC 1 , UC 2 , UC N pour mettre à jour la probabilité de transition nominale de chaque unité de consommation. Le processeur PROC permet également, lors de la phase de gestion mise en œuvre sur la période de temps, de mettre à jour en temps réel le signal de contrôle en fonction de la consommation moyenne cible et de mettre à jour en temps réel la probabilité de transition nominale de chaque unité de consommation sur la base du signal de contrôle ainsi mis à jour. Le module de communication COM est agencé pour communiquer avec les unités de production UP 1 , UP 2 , UP n et les unités de consommation UC 1 , UC 2 , UC N . Par exemple, le module de communication COM permet de recevoir, en provenance de chaque unité de production, les contraintes techniques et le coût de production caractérisant l’unité de production considérée. Le module de communication COM permet aussi de transmettre, à chaque unité de consommation, le signal de contrôle déterminé par le processeur PROC. De même, le module de communication COM permet de transmettre, à chaque unité de production, le programme de production cohérent avec la consommation moyenne cible. Lors de la phase de gestion en temps réel au cours de la période de temps prédéterminée, le module de communication COM permet de transmettre à chaque unité de consommation le signal de contrôle mis à jour à chaque instant. Il est connu de la personne du métier, qu'il existe de nombreux types différents de réseaux de communication de données, par exemple des réseaux de radiocommunication, cellulaires ou non cellulaires, et qu’en fonction du mode de réalisation, le module de communication COM pourra intégrer un ou plusieurs sous-modules de communication, par exemple de communication radiofréquence et être configuré pour l’émission et la réception de signaux radiofréquences, selon une ou plusieurs technologies, telles que TDMA, FDMA, OFDMA, CDMA, ou un ou plusieurs standards de radiocommunication, tels que GSM, EDGE, CDMA, UMTS, HSPA, LTE, LTE-A, WiFi (IEEE 802.11) et WiMAX (IEEE 802.16), ou leurs variantes ou évolutions, actuellement connus ou développés ultérieurement. Le fonctionnement de l’unité de traitement UNT, et plus généralement du portefeuille mixte modélisé par le système SYS, va à présent être décrit en référence à la qui illustre la phase de programmation puis en référence à la qui illustre la phase de gestion, mise en œuvre suite à la phase de programmation. Dans le cadre du procédé décrit ci-après, l’unité de traitement UNT est chargée de la gestion des flexibilités de production et de consommation du système SYS pour une période de temps prédéterminée échantillonnée en une succession d’instants. En d’autres termes, l’unité de traitement UNT a pour rôle la programmation conjointe, pour une gestion coordonnée, de la production des unités de production et de la consommation des unités de consommation. Lors d’une étape S1, l’unité de traitement UNT reçoit au moins une simulation de consommation de chaque unité de consommation du système SYS. L’unité de traitement UNT peut recevoir également des données ou informations nécessaires pour la construction du programme, pour la période de temps à venir, des unités de production UP 1 , UP 2 , UP n et des unités de consommation UC 1 , UC 2 , UC N . Chaque simulation de consommation d’une unité de consommation couvre la période de temps prédéterminée. Chaque simulation relative à une unité de consommation est paramétrée par une probabilité de transition nominale caractérisant la probabilité d’unité de consommation de passer, à un instant donné de la période de temps, d’un état de fonctionnement à un autre. Avantageusement, l’unité de traitement UNT reçoit un nombre élevé de simulations de consommation de chaque unité de traitement. Dans la suite de la description du procédé, l’unité de traitement UNT reçoit le même nombre M de simulations de consommation pour chaque unité de consommation. En d’autres termes, dans le système SYS illustré sur la , l’unité de traitement UNT reçoit M simulations de consommation de l’unité de consommation UC 1 , M simulations de consommation de l’unité de consommation UC 2 et M simulations de consommation de l’unité de consommation UC N . Les simulations de consommation sont par exemple transmises directement par les unités de consommation UC 1 , UC 2 , UC N à l’unité de traitement UNT, laquelle les reçoit via le module de communication COM. Lors d’une étape S2, l’unité de traitement UNT détermine, en fonction de la ou des simulations de consommation reçues, une consommation moyenne cible de l’ensemble d’une ou plusieurs unités de consommation UC 1 , UC 2 , UC N et un signal de contrôle. Plus exactement, la consommation moyenne cible et le signal de contrôle sont déterminés par le processeur PROC de l’unité de traitement UNT en exécutant les instructions du programme informatique stocké dans la mémoire MEM. Par ailleurs, avantageusement, lors de l’étape S2, au moins un profil de production est déterminé pour chaque unité de production. Le ou les profils de production sont déterminés conjointement à la consommation moyenne cible et au signal de contrôle. Le profil de production d’une unité de production dépend de l’état technique dans lequel l’unité de production fonctionne et permet de déterminer, à chaque instant, la quantité d’énergie électrique produite par l’unité de production considérée. Un ensemble de contraintes techniques possibles peut être déterminé pour chaque unité de production et être fourni à l’unité de traitement UNT. De même, le profil de consommation d’une unité de consommation dépend de l’état de l’état de fonctionnement de l’unité de consommation et permet de déterminer, à chaque instant, la quantité d’énergie électrique consommée par l’unité de consommation considérée. Comme expliqué précédemment, la période de temps prédéterminée sur laquelle est programmée le fonctionnement des unités de production UP 1 , UP 2 , UP n et des unités de consommation UC 1 , UC 2 , UC N est échantillonnée en une succession d’instants. Par conséquent, la consommation moyenne cible et le signal de contrôle sont également échantillonnés, leurs valeurs respectives étant connues pour chaque échantillon de temps, donc pour chaque instant. Il est à noter que la consommation cible et le signal de contrôle sont prolongeables par continuité sur la période de temps prédéterminée, par exemple sous la forme de fonctions en escalier. La programmation vise globalement à minimiser le critère suivant qui permet de caractériser le déséquilibre entre l’offre et la demande en énergie électrique au sein d’un système SYS et auquel doit remédier l’unité de traitement UNT par la gestion des flexibilités de consommation : où : - E est l’espérance mathématique ; - f 0 est une fonction pénalisant un écart entre le profil de production de l’ensemble d’une ou plusieurs unités de production et le profil de consommation de l’ensemble d’une ou plusieurs unités de consommation ; - n est le nombre d’unités de production ; - p j (x j ) est le profil de production de la j-ème unité de production lorsque celle-ci fonctionne dans un état technique x j ; - N est le nombre d’unités de consommation ; - c i étant une fonction qui, à un état de la i-ème unité de consommation, associe un profil de consommation, c i (X i ) étant la valeur de la consommation de la i-ème unité de consommation fonctionnant dans l’état X i ; - f j (x j ) est un coût de production de la j-ème unité de production lorsque celle-ci fonctionne dans un état technique x j ; - κ est une constante prédéterminée, permettant notamment de pénaliser les déviations de la consommation par rapport au profil nominal ; - D( , ) est la pseudo distance de Kullback-Leibler entre μ i , qui correspond à la distribution de probabilité contrôlée de la i-ème unité de consommation, et , qui correspond à la distribution de probabilité nominale, donc en l’absence de contrôle, de la i-ème unité de consommation. Dans le critère précédent à minimiser, le premier terme correspond à une pénalité de défaillance, le deuxième terme correspond à un coût de production total de l’ensemble d’une ou plusieurs unités de production UP 1 , UP 2 , UP n et le troisième terme correspond à une pénalité d’écart à la consommation nominale de l’ensemble d’une ou plusieurs unités de consommation UC 1 , UC 2 , UC N . Comme illustré sur la , l’étape S2 de détermination de la consommation moyenne cible et du signal de contrôle est mise en œuvre selon un algorithme itératif. En entrée de cet algorithme, le rang p d’itération est initialisé à 1 et le signal de contrôle est également initialisée à un signal de contrôle initial 0 . Typiquement, le signal de contrôle initial 0 permettant d’initialiser le signal de contrôle préalablement à la mise en œuvre de l’algorithme itératif est le signal de contrôle obtenu à l’issue d’une période de temps précédent la période de temps considérée. En d’autres termes, le dernier signal de contrôle obtenu lors de la période de temps précédant la période de temps à venir pour laquelle une programmation est requise peut être utilisé pour initialiser le signal de contrôle. L’algorithme itératif est défini comme suit : avec : où : - λ p (t k ) est la valeur à l’instant t k du signal de contrôle à l’itération p, le signal de contrôle étant initialisé à un signal de contrôle initial λ 0 prédéterminé ; - ρ p est un réel positif ; - f 0 est une fonction pénalisant un écart entre un profil de production de l’ensemble d’une ou plusieurs unités de production et un profil de consommation de l’ensemble d’une ou plusieurs unités de consommation ; - T est le nombre d’instants de la succession d’instants de la période de temps prédéterminée ; - est un profil de défaillance à l’itération p, (t k ) est la valeur à l’instant t k du profil de défaillance ; - est un réel caractérisant une contrainte de déviation en énergie nulle de l’ensemble d’une ou plusieurs unités de consommation par rapport à une consommation nominale ; - n est le nombre d’unités de production ; - A j est un ensemble de contraintes techniques constituant l’ensemble des états techniques dans lesquels la j-ème unité de production peut fonctionner ; - p j (x j ) est le profil de production de la j-ème unité de production lorsque celle-ci fonctionne dans un état technique x j , p j (x j )(t k ) étant la valeur à l’instant t k du profil de production p j (x j ) ; - f j (x j ) est un coût de production de la j-ème unité de production lorsque celle-ci fonctionne dans un état technique x j ; - N est le nombre d’unités de consommation ; - M est le nombre de simulations de consommation de chaque unité de consommation ; - κ est une constante prédéterminée ; - exp(·) est la fonction exponentielle ; et - U i,q (t k ) est la valeur à l’instant t k de la consommation de la i-ème unité de consommation selon la q-ième simulation de consommation Par ailleurs, l’algorithme itératif est interrompu lorsqu’un critère d’arrêt prédéfini est vérifié. Dans l’exemple illustré en , le critère d’arrêt est défini comme suit : [Math. 22] où : - est une norme ; et - ε est un réel positif. Bien entendu, la personne du métier comprend que d’autres critères d’arrêt peuvent être définis. Par exemple, le critère d’arrêt de l’algorithme itératif peut être l’atteinte ou le dépassement d’un nombre prédéterminé d’itérations. La formule du critère d’arrêt proposé plus haut permet de caractériser que la différence entre le profil de défaillance relatif notamment à la pénalité d’écart entre profil de production et profil de consommation, d’une part, et la différence entre le profil de production et la consommation, d’autre part, doit être arbitrairement aussi faible que possible. Lorsque le critère d’arrêt est vérifié, il est mis fin à l’algorithme itératif qui s’arrête donc à une itération p donnée. L’algorithme a alors permis de générer autant de signaux de contrôle qu’il y a eu d’itération, sans compter le signal de contrôle initial λ 0 . Le signal de contrôle cherché est alors à chercher parmi l’ensemble de ces signaux de contrôle. Par exemple, si un nombre α d’itérations ont été nécessaires, le signal de contrôle déterminé à l’issue de l’étape S2 est le signal de contrôle initial λ 0 ou un des signaux de contrôle parmi le nombre de signaux de contrôle générés au cours de l’algorithme itératif. Par exemple, le signal de contrôle sélectionné correspond au signal de contrôle, parmi tous les signaux de contrôle générés jusqu’à la dernière itération vérifiant le critère d’arrêt, pour lequel un saut de dualité est minimal. En d’autres termes, le rang p * du signal de contrôle sélectionné vérifie : [Math. 23] avec : La fonction Gap correspond au saut de dualité du signal de contrôle λ p obtenu à l’itération p de l’algorithme. Le signal de contrôle déterminé * correspond ainsi au signal de contrôle déterminé à l’itération de rang p * , soit : Le rang p * de l’itération à laquelle le signal de contrôle permet d’atteindre, parmi tous les signaux de contrôle générés, un saut de dualité minimal permet également de déterminer la consommation moyenne cible de l’ensemble d’une ou plusieurs unités de consommation UC 1 , UC 2 , UC N . Comme expliqué précédemment, la consommation moyenne cible permet de déterminer non seulement la consommation des unités de consommations UC 1 , UC 2 , UC N mais aussi le programme de production des unités de production UP 1 , UP 2 , UP n . Plus précisément, le signal de contrôle transmis aux unités de consommation UC 1 , UC 2 , UC N vise à atteindre, pour chaque unité de consommation, la consommation moyenne cible sur la période de temps prédéterminée. La consommation moyenne cible est, elle, utilisée pour déterminer un programme de production des unités de production UP 1 , UP 2 , UP n cohérent avec la consommation moyenne cible de manière à ce que celles-ci adaptent leur production, sur la période de temps prédéterminée, aux besoins en énergie électrique des unités de consommation UC 1 , UC 2 , UC N . La consommation moyenne cible est donnée par la formule ci-après : [Math. 26] où : - r(t k ) est la valeur à l’instant t k de la consommation moyenne cible de l’ensemble d’une ou plusieurs unités de consommation ; et - p * est le rang de l’itération, parmi toutes les itérations jusqu’à la dernière itération vérifiant le critère d’arrêt, pour laquelle un saut de dualité est minimal. Lors d’une étape S3, l’unité de traitement UNT transmet le signal de contrôle déterminé à chaque unité de consommation. Le signal de contrôle transmis est ainsi commun aux unités de consommation UC 1 , UC 2 , UC N du système SYS. Le signal de contrôle est transmis à chaque unité de consommation dans le but d’atteindre, sur la période de temps à venir, la consommation moyenne cible. Le signal de contrôle est par exemple transmis par le module de communication COM de l’unité de traitement UNT. Toujours lors de cette étape S3, l’unité de traitement transmet, également via le module de communication COM, un programme de production à chaque unité de production du système SYS. En effet, étant donné que les unités de consommation UC 1 , UC 2 , UC N sont configurées pour suivre le signal de contrôle et atteindre la consommation moyenne cible, un programme de production complémentaire, conforme à cette consommation moyenne cible, doit être ainsi transmis. Ainsi, le système SYS peut atteindre, sur la période de temps à venir, un équilibre entre l’offre et la demande d’énergie électrique. L’écart entre la quantité d’énergie électrique produite et la quantité d’énergie électrique consommée est ainsi minimisé. Avantageusement, le programme de production comprend en outre, pour chaque unité de production, au moins un profil de production. Comme expliqué précédemment, le programme de production transmis à la j-ème unité de production peut comprendre une fonction p j qui, à chaque état technique x j des contraintes techniques A j de la j-ème unité de production, associe un profil de production p j (x j ). Le profil de production p j (x j ) permet de déterminer, à un instant t k donné de la période de temps prédéterminée à venir, la production en énergie électrique p j (x j )(t k ) de la j-ème unité de production fonctionnant dans l’état technique x j . Lors d’une étape S4, la probabilité de transition nominale est mise à jour, au niveau de chaque unité de consommation, sur la base du signal de contrôle reçu et de la fonction locale caractérisant les spécificités de fonctionnement de l’unité de consommation considérée. Comme expliqué précédemment, la fonction locale associée à une unité de consommation dépend par exemple d’une utilisation programmée d’équipements de l’unité de consommation, de caractéristiques techniques de ces équipements et d’aléas locaux. Par ailleurs, la fonction locale d’une unité de consommation dépend en outre du signal de contrôle reçu, de sorte que la fonction locale ne peut être déterminée que sur réception du signal de contrôle en provenance de l’unité de traitement UNT, et notamment du module de communication COM. Typiquement, la probabilité de transition nominale de la i-ème unité de consommation est mise à jour comme suit lors de la phase de programmation : avec : [Math. 28] où : - X 1 , X 2 et Y sont des états possibles d’une unité de consommation, dX 1 , dX 2 et dY sont les différentielles de tels états, - P 0 i (t k ,X 1 ,dX 2 ) est la probabilité de transition nominale caractérisant la probabilité de la i-ème unité de consommation de passer, à l’instant t k , d’un état de fonctionnement X 1 à un autre état de fonctionnement X 2 ; - Q 0 i (t k ,X 1 ,dX 2 ) est la probabilité de transition nominale mise à jour lors de la phase initiale de programmation caractérisant la probabilité de la i-ème unité de consommation de passer, à l’instant t k , d’un état de fonctionnement X 2 à un autre état de fonctionnement X 2 ; - V i (t k ,X) est la valeur à l’instant t k de la fonction locale caractérisant les spécificités de fonctionnement de la i-ème unité de consommation et de caractéristiques techniques desdits équipements, lorsque la i-ème unité de consommation fonctionne dans l’état X ; - exp(·) est la fonction exponentielle ; - N est le nombre d’unités de consommation ; - κ est une constante prédéterminée ; - λ * (t k ) est la valeur à l’instant t k du signal de contrôle déterminé ; et - c i étant une fonction qui, à un état de la i-ème unité de consommation, associe un profil de consommation, c i (Y) étant la valeur de la consommation de la i-ème unité de consommation fonctionnant dans l’état Y. Il est à noter que les états des unités de consommation UC 1 , UC 2 , UC N ne sont pas nécessairement des états discrets et peuvent être continus. Ainsi, certains équipements industriels ou domestiques peuvent fonctionner selon des états discrets avec un premier état correspondant au mode « marche » et un second état correspondant au mode « arrêt ». Toutefois, certains équipements industriels ou domestiques peuvent fonctionner selon des états continus, comme par exemple un système de climatisation ou un chauffe-eau électrique à accumulation qui, lors du mode de chauffe, peut fonctionner à des températures différentes. Par conséquent, il est pertinent dans la formule précédente de mise à jour de la probabilité de transition nominale de considérer la différentielle des états de fonctionnement pour tenir compte du fait que les états de fonctionnement des unités de consommation UC 1 , UC 2 , UC N peuvent être continus et non discrets. L’étape S4, avec la mise à jour de la probabilité de transition nominale de chaque unité de consommation clôt la phase de programmation du procédé de gestion selon l’invention. A ce stade, chaque unité de consommation dispose d’une probabilité de transition nominale pour la période de temps prédéterminée à venir. Cette probabilité de transition nominale varie, d’un instant donnée de la période de temps prédéterminée à un autre, en fonction du signal de contrôle, lui-même échantillonnée sur la période de temps prédéterminée, et de la fonction locale, dépendant également du signal de contrôle. Comme expliqué précédemment, le fonctionnement en temps réel de chaque unité de consommation peut être modélisé par un comportement aléatoire basé sur la probabilité de transition nominale. Ainsi le fonctionnement de chaque équipement industriel ou domestique au niveau de chaque unité de consommation dépend de la probabilité de transition nominale de l’unité de consommation. L’ensemble des unités de consommation UC 1 , UC 2 , UC N vise à atteindre la consommation moyenne cible, également transmise lors de l’étape S3 par le module de communication COM de l’unité de traitement UNT, déterminée à l’issue de l’algorithme itératif. Par ailleurs, dans le cadre du procédé de gestion selon l’invention, la phase de programmation peut avantageusement être complétée par une phase de gestion en temps réel, mise en œuvre au cours de la période de temps prédéterminée et donc postérieurement à la phase de programmation. Une telle phase de gestion permet une articulation cohérente avec la phase de programmation. A noter que, en l’absence de phase de gestion en temps réel, la probabilité de transition nominale de chaque unité de consommation peut être programmée comme suit sur la période de temps prédéterminée : Une telle probabilité de transition nominale vise à calculer une déviation préliminaire en minimisant la quantité suivante : [Math. 29] avec: La phase de gestion mise en œuvre sur la période de temps prédéterminée est décrite ci-après en référence à la . Une telle phase de gestion vise à mettre à mettre à jour en temps réel signal de contrôle pour mettre à jour, également en temps réel, la probabilité de transition nominale au niveau de chaque unité de consommation du système SYS. Lors de la phase de gestion, l’ensemble d’une ou plusieurs unités de production UP 1 , UP 2 , UP n suit le programme de production transmis par l’unité de traitement UNT. La production de chaque unité de production est ainsi programmée sur la période de temps prédéterminée à venir. Lors d’une étape S5, l’unité de traitement UNT, et plus précisément le processeur PROC, met à jour en temps réel le signal de contrôle en fonction de la consommation moyenne cible. Pour ce faire, l’unité de traitement UNT prévoit, à chaque instant de la période de temps prédéterminée, une consommation moyenne à l’instant suivant de l’ensemble d’une ou plusieurs unités de consommation UC 1 , UC 2 , UC N en fonction de la probabilité de transition nominale mise à jour de chaque unité de consommation. L’unité de traitement UNT met ensuite à jour, à chaque instant, le signal de contrôle sur la base d’une comparaison entre la consommation moyenne à l’instant suivant prévue et de la consommation moyenne cible. Dans le mode de réalisation illustré sur la , la mise à jour du signal de contrôle lors de l’étape S5 de la phase de gestion est mise en œuvre comme suit : L’unité de traitement UNT construit un modèle de consommation prédictif permettant d’approximer, en fonction du signal de contrôle à un instant donné, la consommation moyenne prévue à l’instant suivant. En d’autres termes, une fonction F est générée par l’unité de traitement UNT. Cette fonction F prend en entrée la valeur du signal de de contrôle à un instant donné et renvoie la consommation moyenne de l’ensemble d’une ou plusieurs unités de consommation UC 1 , UC 2 , UC N prévue pour l’instant suivant. Une fois le modèle de consommation prédictif construit, celui-ci est modifié par apprentissage à chaque instant de la période de temps prédéterminée. A chaque instant de la période de temps prédéterminée, l’unité de traitement UNT compare une consommation moyenne réelle et une consommation moyenne prévue à cet instant de l’ensemble d’une ou plusieurs unités de consommation UC 1 , UC 2 , UC N . Une telle comparaison permet de mesurer la différence ou la distance entre la consommation moyenne approximée et la consommation moyenne réelle à cet instant. Ensuite, l’unité de traitement UNT, et plus exactement le processeur PROC, adapte par apprentissage, à chaque instant, le modèle de consommation prédictif, donc la fonction F, sur la base de la comparaison effectuée. Enfin, l’unité de traitement met à jour, à chaque instant, le signal de contrôle par application réciproque du modèle de consommation prédictif à la consommation moyenne cible visée à l’instant suivant. L’application réciproque F -1 permet ainsi de s’adapter aux contraintes, aléas et spécificités des unités de consommation UC 1 , UC 2 , UC N du système SYS. En d’autres termes, à l’issue de l’étape S5, le signal de contrôle est mis à jour comme suit : où ζ (t k ) est la valeur à l’instant t k du signal de contrôle mis à jour. Lors d’une étape S6, l’unité de traitement UNT transmet, à chaque instant de la période de temps prédéterminée, le signal de contrôle mis à jour à chaque unité de consommation. Comme expliqué précédemment, le signal de contrôle mis à jour est par exemple transmis par le module de communication COM. Il est à noter que le signal de contrôle mis à jour au cours de la phase de gestion en temps réel est commun à tous les unités de consommation UC 1 , UC 2 , UC N de même que le signal de contrôle obtenu à l’issue de la phase de programmation. Enfin, lors d’une étape S7, la probabilité de transition nominale est mise à jour, au niveau de chaque unité de consommation, sur la base du signal de contrôle mis à jour reçu. Par exemple, la probabilité de transition nominale peut être mise à jour comme suit lors de la phase de gestion : [Math. 32] où : - Q i (t k ,X 1 ,dX 2 ) est la probabilité de transition nominale mise à jour lors de la phase de gestion caractérisant la probabilité de la i-ème unité de consommation de passer, à l’instant t k , d’un état de fonctionnement X 1 à un autre état de fonctionnement X 2 ; - (t k ) est la valeur à l’instant t k du signal de contrôle mis à jour. La personne du métier comprend que, au cours de la période de temps prédéterminée pendant laquelle est mise en œuvre la phase de gestion en temps réel, chaque unité de consommation du système SYS fonctionne selon la probabilité de transition nominale mise à jour en temps réel. Ainsi, à chaque instant de la période de temps, le signal de contrôle est mis à jour par le processeur PROC au niveau de l’unité de traitement UNT. Le signal de contrôle ainsi modifié est transmis à chaque unité de consommation. La probabilité de transition nominale est mise à jour, également à chaque instant, en fonction du signal de contrôle reçu. Le fonctionnement de chaque unité de consommation est paramétré par la probabilité de transition nominale. En d’autres termes, l’unité de consommation passe d’un état ou mode de fonctionnement à un autre conformément à la probabilité de transition nominale mise à jour. Au niveau de la i-ème unité de consommation, la fonction c i permet de déterminer, pour chaque état de fonctionnement, le profil de consommation c i (X i ) associé. Le profil de consommation c i (X i ) permet de déterminer, à chaque instant t k de la période de temps, la consommation c i (X i )(t k ) de la i-ème unité de consommation. De même, les unités de production UP 1 , UP 2 , UP n suivent le programme de production, sur la période de temps, transmis par l’unité de traitement UNT. Comme expliqué précédemment, ce programme de production est cohérent avec la consommation moyenne cible et permet d’améliorer l’équilibre entre l’offre et la demande en énergie électrique. Ainsi, les équipements domestiques ou industriels localisés au niveau de chaque unité de consommation modulent leur consommation en énergie électrique, fournie par l’ensemble d’une ou plusieurs unités de production UP 1 , UP 2 , UP n en modifiant leur état ou mode de fonctionnement en temps réel en suivant la probabilité de transition nominale mise à jour à chaque réception du signal de contrôle, également mis à jour en temps réel. L’utilisation du modèle prédictif permet de se rapprocher de la consommation moyenne cible déterminée lors de la phase de programmation tout en étant compte, pour chaque unité de consommation, des données locales spécifiques. Procédé, mis en œuvre par des moyens informatiques, de gestion coordonnée, sur une période de temps prédéterminée échantillonnée en une succession d’instants, d’un ensemble d’une ou plusieurs unités de production (UP 1 , UP 2 , UP n ) et d’un ensemble d’une ou plusieurs unités de consommation (UC 1 , UC 2 , UC N ) alimentées en énergie électrique par ledit ensemble d’une ou plusieurs unités de production, ledit procédé comprenant, lors d’une phase de programmation: - recevoir (S1) au moins une simulation de consommation de chaque unité de consommation, chaque simulation de consommation d’une unité de consommation couvrant ladite période de temps et étant paramétrée par une probabilité de transition nominale caractérisant la probabilité de ladite unité de consommation de passer, à un instant donné de ladite période de temps, d’un état de fonctionnement à un autre, - déterminer (S2), en fonction de la ou des simulations de consommation reçues, une consommation moyenne cible de l’ensemble d’une ou plusieurs unités de consommation et un signal de contrôle échantillonnés sur ladite période de temps, - transmettre (S3), à chaque unité de consommation, ledit signal de contrôle pour mettre à jour (S4), au niveau de chaque unité de consommation, la probabilité de transition nominale sur la base du signal de contrôle reçu et d’une fonction locale caractérisant les spécificités de fonctionnement de l’unité de consommation considérée, et - transmettre (S3), à chaque unité de production, un programme de production sur ladite période de temps basée sur ladite consommation moyenne cible. Procédé de gestion selon la revendication 1, dans lequel le programme de production comprend au moins, pour chaque unité de production, un profil de production déterminé conjointement à la consommation moyenne cible et au signal de contrôle. Procédé selon la revendication 2, dans lequel le signal de contrôle échantillonné est déterminé lors de la phase de programmation selon un algorithme itératif défini comme suit : [Math. 32] avec : [Math. 33] [Math. 34] [Math. 35] [Math. 36] où : λ p (t k ) est la valeur à l’instant t k du signal de contrôle à l’itération p, le signal de contrôle étant initialisé à un signal de contrôle initial 0 prédéterminé ; - ρ p est un réel positif ; - f 0 est une fonction pénalisant un écart entre un profil de production de l’ensemble d’une ou plusieurs unités de production et un profil de consommation de l’ensemble d’une ou plusieurs unités de consommation ; - T est le nombre d’instants de la succession d’instants de la période de temps prédéterminée ; - est un profil de défaillance à l’itération p, (t k ) est la valeur à l’instant t k du profil de défaillance ; - est un réel caractérisant une contrainte de déviation en énergie nulle de l’ensemble d’une ou plusieurs unités de consommation par rapport à une consommation nominale ; - n est le nombre d’unités de production ; - A j est un ensemble de contraintes techniques constituant l’ensemble des états techniques dans lesquels la j-ème unité de production peut fonctionner ; - p j (x j ) est le profil de production de la j-ème unité de production lorsque celle-ci fonctionne dans un état technique x j , p j (x j )(t k ) étant la valeur à l’instant t k du profil de production p j (x j ) ; - f j (x j ) est un coût de production de la j-ème unité de production lorsque celle-ci fonctionne dans un état technique x j ; - N est le nombre d’unités de consommation ; - M est le nombre de simulations de consommation de chaque unité de consommation ; - κ est une constante prédéterminée ; - exp(·) est la fonction exponentielle ; et - U i,q (t k ) est la valeur à l’instant t k de la consommation de la i-ème unité de consommation selon la q-ième simulation de consommation, l’algorithme itératif étant interrompu lorsqu’un critère d’arrêt prédéfini est vérifié. Procédé selon la revendication 3, dans lequel le critère d’arrêt est défini comme suit : [Math. 37] où : - est une norme ; et - ε est un réel positif Procédé selon la revendication 3 ou 4, dans lequel le signal de contrôle initial permettant d’initialiser le signal de contrôle préalablement à la mise en œuvre de l’algorithme itératif est le signal de contrôle obtenu à l’issue d’une période de temps précédent la période de temps considérée. Procédé selon l’une des revendications 3 à 5, dans lequel le signal de contrôle déterminé à l’issue de l’algorithme itératif est défini comme suit : [Math. 38] où p * est le rang de l’itération, parmi toutes les itérations jusqu’à la dernière itération vérifiant le critère d’arrêt, pour laquelle un saut de dualité est minimal, soit : [Math. 39] avec : [Math. 40] Procédé selon l’une des revendications 3 à 6, dans lequel la consommation moyenne cible est définie comme suit : [Math. 41] où : - r(t k ) est la valeur à l’instant t k de la consommation moyenne cible de l’ensemble d’une ou plusieurs unités de consommation ; et - p * est le rang de l’itération, parmi toutes les itérations jusqu’à la dernière itération vérifiant le critère d’arrêt, pour laquelle un saut de dualité est minimal, soit : [Math. 42] avec : [Math. 43] Procédé selon l’une des revendications précédentes, dans lequel la probabilité de transition nominale est mise à jour comme suit lors de la phase de programmation : [Math. 44] avec : [Math. 45] où : - X 1 , X 2 et Y sont des états possibles d’une unité de consommation, dX 1 , dX 2 et dY sont les différentielles de tels états, - P 0 i (t k ,X 1 ,dX 2 ) est la probabilité de transition nominale caractérisant la probabilité de la i-ème unité de consommation de passer, à l’instant t k , d’un état de fonctionnement X 1 à un autre état de fonctionnement X 2 ; - Q 0 i (t k ,X 1 ,dX 2 ) est la probabilité de transition nominale mise à jour lors de la phase initiale de programmation caractérisant la probabilité de la i-ème unité de consommation de passer, à l’instant t k , d’un état de fonctionnement X 2 à un autre état de fonctionnement X 2 ; - V i (t k ,X) est la valeur à l’instant t k de la fonction locale caractérisant les spécificités de fonctionnement de la i-ème unité de consommation et de caractéristiques techniques desdits équipements, lorsque la i-ème unité de consommation fonctionne dans l’état X ; - exp(·) est la fonction exponentielle ; - N est le nombre d’unités de consommation ; - κ est une constante prédéterminée ; - λ * (t k ) est la valeur à l’instant t k du signal de contrôle déterminé ; et - c i étant une fonction qui, à un état de la i-ème unité de consommation, associe un profil de consommation, c i (Y) étant la valeur de la consommation de la i-ème unité de consommation fonctionnant dans l’état Y . Procédé selon l’une des revendications précédentes, dans lequel la fonction locale caractérisant les spécificités de fonctionnement d’une unité de consommation dépend d’une utilisation programmée d’équipements industriels de ladite unité de consommation, de caractéristiques techniques desdits équipements industriels et d’aléas locaux. Procédé selon l’une des revendications précédentes, comprenant en outre, suite à la phase de programmation, une phase de gestion mise en œuvre sur ladite période de temps comprenant la mise à jour (S5) en temps réel du signal de contrôle en fonction de ladite consommation moyenne cible ; et la mise à jour (S7) en temps réel de la probabilité de transition nominale de chaque unité de consommation sur la base du signal de contrôle mis à jour. Procédé selon la revendication 8, comprenant en outre, suite à la phase de programmation, une phase de gestion mise en œuvre sur ladite période de temps comprenant la mise à jour en temps réel du signal de contrôle en fonction de ladite consommation moyenne cible ; et la mise à jour en temps réel de la probabilité de transition nominale de chaque unité de consommation sur la base du signal de contrôle mis à jour, dans lequel la probabilité de transition nominale est mise à jour comme suit lors de la phase de gestion : [Math. 46] où : - Q i (t k ,X 1 ,dX 2 ) est la probabilité de transition nominale mise à jour lors de la phase de gestion caractérisant la probabilité de la i-ème unité de consommation de passer, à l’instant t k , d’un état de fonctionnement X 1 à un autre état de fonctionnement X 2 ; - (t k ) est la valeur à l’instant t k du signal de contrôle mis à jour. Procédé selon la revendication 10 ou 11, dans lequel la phase de gestion est mise en œuvre comme suit : - prévoir, à chaque instant de ladite période de temps, une consommation moyenne à l’instant suivant de l’ensemble d’une ou plusieurs unités de consommation en fonction de la probabilité de transition nominale mise à jour de chaque unité de consommation, - mettre à jour, à chaque instant, le signal de contrôle sur la base d’une comparaison entre la consommation moyenne à l’instant suivant prévue et de la consommation moyenne cible, et - transmettre (S6), à chaque instant, ledit signal de contrôle mis à jour à chaque unité de consommation pour mettre à jour, au niveau de chaque unité de consommation, la probabilité de transition nominale sur la base du signal de contrôle mis à jour reçu. Procédé selon la revendication 12, dans laquelle la mise à jour du signal de contrôle lors de la phase de gestion est mise en œuvre comme suit : - construire (S51) un modèle de consommation prédictif permettant d’approximer, en fonction du signal de contrôle à un instant donné, la consommation moyenne prévue à l’instant suivant, puis, à chaque instant de la période de temps prédéterminée : - comparer une consommation moyenne réelle et une consommation moyenne prévue audit instant de l’ensemble d’une ou plusieurs unités de consommation, - adapter (S52) par apprentissage, à chaque instant, ledit modèle de consommation prédictif sur la base de ladite comparaison, et - mettre à jour (S53), à chaque instant, le signal de contrôle par application réciproque du modèle de consommation prédictif à la consommation moyenne cible visée à l’instant suivant. Procédé selon l’une des revendications 10 à 13, dans lequel, au cours de la phase de gestion mise en œuvre sur ladite période de temps, chaque unité de consommation fonctionne selon la probabilité de transition nominale mise à jour en temps réel, tandis que chaque unité de production suit le programme de production reçu. Programme informatique, comprenant des instructions pour la mise en œuvre du procédé selon l’une des revendications précédentes lorsque lesdites instructions sont mises en œuvre par au moins un processeur (PROC). Unité de traitement (UNT) pour la gestion coordonnée, sur une période de temps prédéterminée échantillonnée en une succession d’instants, d’un ensemble d’une ou plusieurs unités de production (UP 1 , UP 2 , UP n ) et d’un ensemble d’une ou plusieurs unités de consommation (UC 1 , UC 2 , UC N ) alimentées en énergie électrique par ledit ensemble d’une ou plusieurs unités de production, ladite unité de traitement étant configurée pour, lors d’une phase de programmation: - recevoir au moins une simulation de consommation de chaque unité de consommation, chaque simulation de consommation d’une unité de consommation couvrant ladite période de temps et étant paramétrée par une probabilité de transition nominale caractérisant la probabilité de ladite unité de consommation de passer, à un instant donné de ladite période de temps, d’un état de fonctionnement à un autre, - déterminer, en fonction de la ou des simulations de consommation reçues, une consommation moyenne cible de l’ensemble d’une ou plusieurs unités de consommation et un signal de contrôle échantillonnés sur ladite période de temps, - transmettre, à chaque unité de consommation, ledit signal de contrôle pour mettre à jour, au niveau de chaque unité de consommation, la probabilité de transition nominale sur la base du signal de contrôle reçu et d’une fonction locale caractérisant les spécificités de fonctionnement de l’unité de consommation considérée, et - transmettre, à chaque unité de production, un programme de production sur ladite période de temps basée sur ladite consommation moyenne cible.