FR 3072794 A1 20190426 FR 1759892 A 20171020 PROCEDE DE DETERMINATION DES PARAMETRES D'UN MODELE SIMPLIFIE D'UN SYSTEME DE STOCKAGE DE L'ENERGIE, PROCEDE DE PILOTAGE UTILISANT UN TEL MODELE ET DISPOSITIF ASSOCIE DOMAINE TECHNIQUE DE L'INVENTION Le domaine technique de l'invention est celui de la gestion de l'énergie. La présente invention concerne un procédé de détermination des paramètres d'un modèle physique et en particulier un procédé de détermination des paramètres d'un modèle simplifié d'un système de stockage de l'énergie. ARRIERE-PLAN TECHNOLOGIQUE DE L'INVENTION Dans le cadre des énergies renouvelables, il est primordial d'optimiser le stockage de l'énergie et l'utilisation de cette énergie stockée pour atténuer le caractère intermittent de certaines sources d'énergie telles que les éoliennes ou les panneaux photovoltaïques. Le stockage est, de manière générale, assuré par un système de stockage de l'énergie comprenant un dispositif de stockage en charge du stockage de l'énergie proprement dit et d'un dispositif de conversion qui va assurer la charge ou la décharge du dispositif de stockage ainsi que la conversion de l'énergie stockée sous une forme adaptée à une utilisation domestique ou industrielle, par exemple en convertissant un courant continu en courant alternatif. Il est à noter qu'un système de stockage inclut souvent des auxiliaires qui peuvent être associés au dispositif de stockage ou au dispositif de conversion. Dans la suite du texte, lorsque l'on parle de système de stockage, ces systèmes auxiliaires sont considérés comme inclus dans ledit système de stockage. Afin de piloter un tel système de stockage de manière efficace, il est primordial de prévoir le comportement de ce dernier de sorte à anticiper la charge, la décharge ainsi que la puissance que le dispositif peut fournir ou absorber. Ce pilotage s'effectue par la mise à jour de consignes de fonctionnement, ladite mise à jour étant effectuée à intervalles réguliers, en général de l'ordre de la minute voire de l'heure. Afin de déterminer les consignes adaptées au pilotage de l'installation, il est 2 donc nécessaire de simuler le comportement du système de stockage avant l'établissement de chaque nouvelle consigne. Pour cela, un modèle représentant le dispositif de stockage et un modèle représentant le dispositif de conversion sont assemblés de sorte à constituer un modèle complexe. Afin de piloter le système de 5 stockage de l'énergie, il peut sembler raisonnable d'effectuer une simulation avec un pas de temps du même ordre de grandeur que le pas temps de pilotage (i.e. de quelques minutes voire quelques heures), la simulation ainsi effectuée permettant de générer une consigne. Cependant, il est en général nécessaire d'adopter un pas de temps bien inférieur au pas de temps de pilotage (i.e. de l'ordre de la seconde 10 voire de la milliseconde) afin d'obtenir une précision suffisante pour générer une consigne fiable. Ainsi, pour chaque nouvelle consigne, il est nécessaire d'effectuer une simulation gourmande en ressource de calcul et en mémoire. Il existe donc un besoin d'un procédé de pilotage dans lequel les consignes sont générées en ayant recours à une simulation dont le pas de temps est du même 15 ordre de grandeur que le pas de temps de pilotage lui-même, et qui nécessite ainsi une puissance de calcul moindre (ou un temps de calcul moindre pour une puissance de calcul donnée) tout en garantissant une précision suffisante. Il existe donc également un besoin d'un modèle simplifié permettant de mettre en oeuvre un tel procédé de pilotage. 20 RESUME DE L'INVENTION L'invention offre une solution aux problèmes évoqués précédemment, en permettant d'obtenir un modèle simplifié dont le temps de simulation est limité. En effet, le modèle simplifié obtenu à l'aide d'un procédé selon l'invention permet d'effectuer une simulation d'un système de stockage de l'énergie avec un pas de 25 temps du même ordre de grandeur que le pas temps de pilotage dudit système, tout en garantissant une précision des résultats simulés similaire à celle obtenue avec une simulation ayant un pas de temps beaucoup plus petit que le pas de temps du pilotage lui-même. 3 Un premier aspect de l'invention concerne un procédé de détermination des paramètres d'un modèle simplifié d'un système de stockage de l'énergie, ledit système comprenant un dispositif de stockage de l'énergie et un dispositif de conversion, ledit système pouvant être modélisé au moyen d'un modèle complexe 5 incluant un modèle du dispositif de stockage de l'énergie et un modèle du dispositif de conversion ; ledit modèle complexe recevant en entrée une puissance de consigne P . ac_sp et un état charge SOC, et fournissant en sortie l'état de charge SOC du dispositif de stockage ainsi que la puissance Pac en sortie du dispositif de stockage ; ledit procédé étant caractérisé en ce qu'il comprend : 10 - une première étape de mise en oeuvre d'une pluralité de simulations du système de stockage de l'énergie à l'aide du modèle complexe, chaque simulation étant par exemple effectuée avec un pas de temps At; - une deuxième étape de calcul à partir des résultats obtenus lors de la première étape : 15 - d'une table de la variation temporelle de l'état de charge du système en fonction de la puissance de consigne P . ac_sp et de l'état de charge SOC ; - d'une table de puissance maximale fournie en fonction de l'état de charge SOC ; 20 - d'une table de puissance minimale fournie en fonction de l'état de charge SOC. Le modèle simplifié obtenu permet d'attribuer en fonction d'une puissance de consigne P . ac_sp et d'un état de charge SOC p fournis en entrée, et à partir des tables déterminées lors de la deuxième étape, une puissance Pac et un état de charge du 25 système SOC. Dans la suite, par convention, la puissance Pac lors de la charge est considérée comme négative et la puissance Pac lors de la décharge est considérée comme positive. Grâce à l'invention, il n'est plus nécessaire de choisir entre un pas de temps de simulation du même ordre de grandeur que le pas de temps de pilotage (i.e. la durée séparant deux mises à jour des consignes de pilotage) conduisant à une précision insatisfaisante, et un pas de temps de simulation bien plus 4 faible permettant d'obtenir une simulation certes précise, mais gourmande en ressources de calcul et en mémoire. En effet, le modèle simplifié obtenu à l'aide d'un procédé selon un premier aspect de l'invention permet d'adopter un pas de temps de simulation du même ordre de grandeur que le pas de temps de pilotage 5 tout en conservant une précision suffisante pour ledit pilotage. Outre les caractéristiques qui viennent d'être évoquées dans le paragraphe précédent, le procédé selon un premier aspect de l'invention peut présenter une ou plusieurs caractéristiques complémentaires parmi les suivantes, considérées individuellement ou selon toutes les combinaisons techniquement possibles. 10 Avantageusement, chaque valeur de la table de la variation de l'état de charge du système est obtenue à l'aide d'une simulation d'une durée tESS effectuée pour un état de charge SOC(j) et une puissance de consigne Pac_sp\., M donnés et appartenant à un premier sous-ensemble de la pluralité de simulations, ladite valeur étant : - égale à la valeur de la variation de l'état de charge ASOC obtenue lors de 15 ladite simulation sur le temps tESS si la puissance moyenne au cours de ladite simulation est égale à la puissance de consigne P . ac_sp. - dans le cas contraire, égale à: interp f[P' sp(i - 1), ],[dsocd(it-i,j) ,PSOtECs(si,j)1, D (; ac sp Y') avec : dSOC(i.-1,j) 20 - la variation de charge obtenue ou calculée lors de la dt simulation correspondant à l'état de charge SOC(j) et à une puissance de consigne Pac_sp(i-1) ; - ASOC(i, j), la variation de charge obtenue lors de la simulation correspondant à l'état de charge SOC(j) et à une puissance de 25 consigne P . ac_spv, ; - , la valeur moyenne de la puissance Pac lors de la simulation correspondant à l'état de charge SOC(j) et à une puissance de consigne P ac_spv, ; - interpaxo, xl], [yo, x} est la fonction qui détermine la valeur de y 30 correspondant à la valeur de x par interpolation à partir des valeurs de 5 x0,x1,y0 et Ainsi, il est possible d'établir la variation temporelle de l'état de charge correspondant à une simulation donnée même lorsque cette dernière n'est pas constante lors de ladite simulation. On entend par égale que la puissance Pac est 5 égale à la puissance de consigne Pw_sp à plus ou moins 5%, voire plus ou moins 2%, de préférence plus ou moins 1%. De manière alternative, chaque valeur de la table de la variation de l'état de charge du système est obtenue à l'aide d'une simulation d'une durée tESS effectuée pour un état de charge SOC(j) et une puissance de consigne Pac_sp,., M donnés et appartenant 10 à un premier sous-ensemble de la pluralité de simulations, ladite valeur étant égale à la valeur de la variation de l'état de charge ASOC obtenue lors de ladite simulation sur le temps tEss. Ainsi, il est possible d'établir la variation temporelle de l'état de charge correspondant à une simulation donnée même lorsque cette dernière n'est pas 15 constante lors de ladite simulation en ayant recours à une approximation peu exigeante en ressource de calcul. En effet, les inventeurs ont mis en évidence que, de manière surprenante, cette approximation permettait d'obtenir de très bons résultats. De manière encore alternative chaque valeur de la table de la variation de l'état de 20 charge du système est obtenue à l'aide d'une simulation d'une durée tESS effectuée pour un état de charge SOC(j) et une puissance de consigne Pac_sp,., M donnés et appartenant à un premier sous-ensemble de la pluralité de simulations, ladite valeur étant : - égale à la valeur de la variation de l'état de charge ASOC obtenue lors de 25 ladite simulation sur le temps tESS si la puissance moyenne au cours de ladite simulation est égale à la puissance de consigne P . ac_sp. Sd SOC(m,j) , d OC(m,j) dt dt - dans le cas contraire, égale à ou est la variation temporelle de charge obtenue lors de la simulation correspondant à l'état de charge SOCi(j) et la puissance de consigne Pac_sp(M), Pac_sp(M) étant la puissance 6 de consigne la plus proche de la puissance de consigne Pac_sp(i) pour laquelle la puissance moyenne au cours de ladite simulation est égale à la puissance de consigne P ( . ac_sp\r11). Ainsi, il est possible d'établir la variation temporelle de l'état de charge 5 correspondant à une simulation donnée même lorsque cette dernière n'est pas constante lors de ladite simulation en ayant recours à une approximation peu exigeante en ressource de calcul. Avantageusement, chaque simulation du premier sous-ensemble de simulations est effectuée pour un état de charge SOC(j) et une puissance de consigne P (i) . ac_sp, 10 données et en ce que chaque état de charge SOC(j) est séparé du précédent SOC(j-1) et/ou du suivant SOC(j+1) par un pas d'état de charge adaptatif et/ou chaque puissance de consigne ac_sp P (i) est séparée de la précédente Pac_sp(i-1) . , et/ou de la suivante Pac_spy (i+ 1) d'un pas de puissance de consigne adaptatif. Ainsi, le nombre de simulations effectuées afin de constituer la table de la variation 15 temporelle de l'état de charge est réduit. Avantageusement, l'étape de simulation est répétée pour une pluralité de durées tESS. Ainsi, il est possible d'intégrer une variable temporelle au modèle simplifié. Avantageusement, la pluralité de simulations comporte une simulation effectuée 20 avec un état de charge initial SOC ini égal à l'état de charge maximum SOCmax, une durée égale à la durée nécessaire au déchargement complet du système de stockage tDch et une puissance de consigne P . ac_sp infinie positive, et le calcul de la table de puissance maximale fournie en fonction de l'état de charge SOC comporte : 25 - une étape de détermination de périodes de durée tESS au sein de ladite simulation, chaque période étant identifiée au moyen d'un entier positif k ; - pour chacune de ces périodes, une étape de détermination de l'état de charge SOCk au début de la période k et de la puissance 7 moyenne k durant la période k. Ainsi, la pluralité de couples (SOCk, k) vient constituer la table de puissance maximale fournie en fonction de l'état de charge SOC. On entend par une puissance de consigne P ac_sp infinie positive une puissance de consigne P . ac_sp largement 5 supérieure, en valeur absolue à la puissance que le système peut fournir. Pour rappel, par convention, une valeur positive de la puissance Pac est associée à la décharge du système de stockage. Avantageusement, la pluralité de simulations comporte une simulation effectuée avec un état de charge initial SOC ini égal à l'état de charge minimum SOCmin, une 10 durée égale à la durée nécessaire au chargement complet du système de stockage tch et une puissance de consigne P . ac_sp infinie négative, et le calcul de la table de puissance minimale fournie en fonction de l'état de charge SOC comporte : - une étape de détermination de périodes de durée tESS au sein de ladite simulation, chaque période étant identifiée au moyen d'un entier 15 positif k' ; - pour chacune de ces périodes, une étape de détermination de l'état de charge SOCk' au début de la période k' et de la puissance moyenne k' durant la période k' ; Ainsi, la pluralité de couples (SOCk', k) vient constituer la table de puissance 20 minimale fournie en fonction de l'état de charge SOC. On entend par une puissance de consigne P ac_sp infinie négative une puissance de consigne P . ac_sp largement supérieure, en valeur absolue à la puissance que le système peut recevoir lors de la charge. Pour rappel, par convention, une valeur négative de la puissance Pac est associée à la charge du système de stockage. 25 Un deuxième aspect de l'invention concerne un procédé de pilotage d'un système de stockage de l'énergie calculant des consignes de fonctionnement dudit système à partir d'un modèle dudit système, ledit modèle étant obtenu à l'aide d'un procédé selon un premier aspect de l'invention. 8 Un troisième aspect de l'invention concerne un dispositif de pilotage d'un système de stockage de l'énergie comprenant des moyens pour envoyer des consignes de fonctionnement au système de stockage de l'énergie, des moyens pour recevoir des données concernant le fonctionnement du système de stockage de l'énergie et 5 des moyens pour mettre en oeuvre un procédé de pilotage selon un deuxième aspect de l'invention. Un quatrième aspect de l'invention concerne un produit programme d'ordinateur comprenant des instructions qui, lorsque le programme est exécuté par un ordinateur, conduisent celui-ci à mettre en oeuvre les étapes du procédé selon un 10 premier aspect de l'invention. Un cinquième aspect de l'invention concerne un produit programme d'ordinateur comprenant des instructions qui conduisent le dispositif de pilotage selon un troisième aspect de l'invention à exécuter les étapes du procédé selon un deuxième aspect de l'invention. 15 Un sixième aspect de l'invention concerne un support lisible par ordinateur, sur lequel est enregistré le programme d'ordinateur selon un quatrième ou un cinquième aspect de l'invention. L'invention et ses différentes applications seront mieux comprises à la lecture de la description qui suit et à l'examen des figures qui l'accompagnent. 20 BREVE DESCRIPTION DES FIGURES Les figures sont présentées à titre indicatif et nullement limitatif de l'invention. - La figure 1 montre un ordinogramme d'un mode de réalisation d'un procédé selon un premier aspect de l'invention. - La figure 2 montre une représentation schématique d'un système de 25 stockage de l'énergie. - La figure 3 montre un schéma de principe d'un procédé selon un premier aspect de l'invention. 9 - Les figures 4A et 4B montrent une simulation mise en jeu dans un procédé selon un premier aspect de l'invention. - Les figures 5A et 5B montrent une simulation mise en jeu dans un procédé selon un premier aspect de l'invention. 5 - Les figures 6A et 6B montrent une illustration en 3D d'une table de la variation temporelle de l'état de charge du système selon un premier aspect de l'invention. - Les figures 7A et 7B montrent une simulation mise en jeu dans un procédé selon un premier aspect de l'invention. 10 - Les figures 8A et 8B montrent une simulation mise en jeu dans un procédé selon un premier aspect de l'invention. - La figure 9 montre un graphique illustrant les performances en temps et en précision d'un modèle obtenu à l'aide d'un procédé selon un premier aspect de l'invention. 15 DESCRIPTION DETAILLEE D'AU MOINS UN MODE DE REALISATION DE L'INVENTION Sauf précision contraire, un même élément apparaissant sur des figures différentes présente une référence unique. Un premier mode de réalisation d'un procédé selon un premier aspect de l'invention 20 illustré à la figure 1 concerne un procédé 100 de détermination des paramètres d'un modèle simplifié MS d'un système de stockage de l'énergie ESS. Le système de stockage de l'énergie ESS illustré à la figure 2 comprend un dispositif de stockage de l'énergie DSE, par exemple une batterie, et un dispositif de conversion DC et peut être modélisé au moyen d'un modèle complexe MC illustré à la figure 3 incluant 25 un modèle du dispositif de stockage de l'énergie MCS et un modèle du dispositif de conversion MCC. Le modèle complexe MC reçoit en entrée une puissance de consigne P . ac_sp et un état de charge SOC, et fournit en sortie l'état de charge SOC du dispositif de stockage ainsi que la puissance Pac en sortie du système de stockage de l'énergie. Le procédé 100 selon l'invention comprend une première 30 étape 101 de mise en oeuvre d'une pluralité de simulations du système de stockage 10 de l'énergie SSE à l'aide du modèle complexe MC, chaque simulation étant par exemple effectuée avec un pas de temps At. Lors de ces simulations, l'état de charge SOC p fourni en entrée du modèle complexe est l'état de charge initial du système pour la première itération puis, pour les itérations suivantes, l'état de 5 charge calculé lors de l'itération précédente. Le procédé selon un premier aspect de l'invention comprend également une deuxième étape 102 de calcul à partir des résultats obtenus lors de la première étape 101 : - d'une table SOCV_TC de la variation temporelle ci dS0tC - de l'état de charge du 10 système en fonction de la puissance de consigne P . ac_sp et de l'état de charge SOC; - d'une table PAC _ MAX_ TC de puissance maximale en fonction de l'état de charge SOC ; - d'une table PAC _ MIN _TC de puissance minimale en fonction de l'état de 15 charge SOC. Pour rappel, par convention, la puissance Pac lors de la charge est considérée comme négative car elle est absorbée et la puissance Pac lors de la décharge est considérée comme positive car elle est fournie à l'extérieur du système considéré. Ainsi, pour un état de charge donnée, la puissance minimale constitue la puissance 20 minimale P . ac_min (négative) que peut absorbée le système pour ledit état de charge et la puissance maximale fournie constitue la puissance maximale Pac_max (positive) que peut fournir le système pour ledit état de charge. Le modèle simplifié obtenu permet d'attribuer en fonction d'une puissance de consigne P . ac_sp et d'un état de charge SOC p fournis en entrée, et à partir des 25 tables SOCV_TC, PAC_MAX_TC, PAC_MIN_TC déterminées lors de la deuxième étape 102, une puissance Pac et un état de charge du système SOC. En effet, à partir d'une puissance de consigne P . ac_sp ainsi que d'un état de charge SOC p fourni en entrée du modèle, il est possible de déterminer la puissance minimale à l'aide de la table PAC _ MIN _TC de puissance minimale et la puissance maximale à l'aide 30 de la table PAC _ MAX_ TC de puissance maximale. Si la puissance de consigne 11 Pc_sp est située dans cet intervalle, alors la puissance délivrée par le système est égale à ladite consigne P . ac_sp, sinon elle est égale à la valeur limite la plus proche de ladite consigne. Il est également possible, à partir de la puissance de consigne Pac_sp une fois saturée et de l'état de charge SOC p fournis en entrée, de dSOC 5 déterminer la variation temporelle de l'état de charge à l'aide de la dt table SOCV_TC de la variation temporelle de l'état de charge ci dS0tC - du système de stockage de l'énergie ESS. L'état de charge SOC p fourni en entrée pourra par exemple correspondre à un état de charge d'une itération précédente lorsque le modèle simplifié est utilisé pour effectuer une simulation. Ainsi, le modèle obtenu 10 permet de modéliser un système de stockage de l'énergie ESS de manière rapide sans pour autant remettre en cause la précision des simulations permettant de modéliser ledit système. Dans un mode de réalisation, chaque valeur de la table SOCV_TC de la variation dSOC temporelle de l'état de charge -dt du système est obtenue à l'aide d'une simulation 15 effectuée pour un état de charge SOC(j) choisi comme état de charge initial SOCini et une puissance de consigne ac_sp,., P M donnés et appartenant à un premier sous- ensemble de la pluralité de simulations, ladite simulation étant effectuée sur une durée tESS. Cette simulation va permettre d'associer à chaque puissance de consigne Pac_sp(i) et à chaque état de charge SOC(j) une variation temporelle de M 20 l'état de charge -ddSOCt J )-Dans un mode de réalisation, le pas de temps At utilisé pour les simulations à l'aide du modèle complexe est choisi de sorte que 102AktEss, de préférence 103AktEss. Les figures 4A-4B et 5A-5B illustrent deux simulations du premier sous ensemble de la pluralité de simulations effectuées pour deux états de charge SOC différents 25 et pour un même temps de simulation tEss. Les figures 4A et 5A reproduisent la puissance Pac en sortie du système de gestion de l'énergie en fonction du temps et les figures 4B et 5B reproduisent l'état de charge SOC du système de gestion de l'énergie en fonction du temps. 12 La simulation illustrée aux figures 4A et 4B décrit une simulation dans laquelle la puissance Pac en sortie du système de stockage de l'énergie sur toute la durée tESS de la simulation est égale à la puissance de consigne P . ac_sp. On entend par égale que la puissance Pac en sortie du le système ESS est égale à la puissance de 5 consigne Pw_sp à plus ou moins 5%, voire plus ou moins 2%, de préférence plus ou moins 1%. Lors de cette simulation, l'état de charge SOC décroît de manière constante de sorte qu'il est aisé, à partir de cette simulation, d'extraire une valeur dSOC de la variation temporelle de l'état de charge -dt. Les figures 5A et 5B décrivent quant à elles une simulation dans laquelle la puissance Pac en sortie du système de 10 stockage de l'énergie sur toutes la durée tESS de la simulation n'est pas constante, mais varie au cours de la simulation. Cette variation peut par exemple s'expliquer par le fait que le dispositif de stockage DS du système de stockage de l'énergie ESS est presque totalement chargé (ou déchargé). Lors de cette simulation, la variation dSOC temporelle de l'état de charge -dt du système de stockage ESS n'est plus 15 constante, l'état de charge SOC ayant tendance à former un plateau sur la fin de la simulation pour les raisons qui viennent d'être évoquées. Dans ce cas de figure, il peut être plus difficile d'attribuer une valeur représentant correctement la variation temporelle de l'état de charge dSOC du système de stockage de l'énergie associée à la puissance de consigne Pac_sp pour un état de charge donné SOC. 20 Afin de prendre en compte ces deux cas de figure, dans un mode de réalisation dit « avec interpolation », le procédé selon un premier aspect de l'invention prévoit que chaque valeur de la table SOCV_TC de correspondance est : - égale à la valeur de la variation de l'état de charge ASOC obtenue lors de ladite simulation sur le temps tESS si la puissance moyenne au cours 25 de ladite simulation est égale à la puissance de consigne Pac_sp (situation illustrée aux figures 2A et 2B) ; - dans le cas contraire (situation illustrée aux figures 3A et 3B), égale à: [dSOC(i - 1, , j) ASOC(i,j)1 interp {[P' sp(i - 1), G Pac(i,j) >1, .,_ , Pac sp(i)} dt LEss avec : 13 dSOC(i-1,j - dt ) la variation temporelle de charge obtenue lors de la simulation ou calculée correspondant à l'état de charge SOC(j) et à un puissance de consigne Pac_sp(i-1) ; - ASOC(i,j), la variation de charge obtenue lors de la simulation 5 correspondant à l'état de charge SOC(j) et à une puissance de consigne P ac_sp(i) ; - , la valeur moyenne de la puissance Pac lors de la simulation correspondant à l'état de charge SOC(j) et à une puissance de consigne P ac_sp(i) ; 10 - interpax0,x1], [yo, x} est la fonction qui détermine la valeur de y correspondant à la valeur de x par interpolation à partir des valeurs de Xo, Xi, y 0 et Autrement dit, on réalise une interpolation à l'aide des simulations précédemment effectuées. Afin de réaliser cette interpolation, il est par exemple envisageable de 15 réaliser, à partir des simulations du premier sous-ensemble, un premier tableau dans lequel les colonnes représentent des états de charges SOC différents, les lignes représentent des puissances de consigne P ac_sp différentes et dont les cases contiennent la puissance moyenne au cours de ladite simulation (cf. tableau 1) et un deuxième tableau dans lequel les colonnes représentent des états 20 de charges SOC différents, les lignes représentent des puissances de consigne P ac_sp différentes et dont les cases contiennent la variation totale de l'état de charge ASOC au cours de la simulation (cf. tableau 2). Plus particulièrement, pour chaque état de charge SOC, une simulation est effectuée pour une puissance de consigne Pac_sp nulle puis les simulations suivantes sont effectuée pour des 25 puissances de consignes P ac_sp allant croissantes. Ainsi, lorsque l'on atteint une puissance P ac_sp qui ne peut pas être maintenue lors de la simulation et qui nécessite donc une interpolation, la variation temporelle de l'état de charge dSOC obtenue à la simulation précédente peut être utilisée pour ladite interpolation. Il est procédé de la même façon pour les valeurs négatives de la 30 puissance de consigne P ac_sp. Autrement dit, si l'on souhaite calculer la table dt 14 dSOC variation temporelle de l'état de charge -dt pour des puissances de consigne Pac_sp allant de -Ps à Ps, on effectue les simulations en commençant par une puissance de consigne Pac_sp nulle que l'on fait croitre entre chaque simulation jusqu'à une valeur de puissance de consigne égale à Ps. Puis, on effectue ensuite les 5 simulations en commençant par une puissance de consigne Pac_sp nulle que l'on fait décroitre entre chaque simulation jusqu'à une valeur de puissance de consigne égale à -Ps. SOC(0) ... SOC(j-1) SOC(j) SOC(j+1) ... SOC(W) Pac_sp(k) (k,O) (k,j-1) (k,j) (k,j+1) (k,k') ... Pac_sp(i+1) (i+1 ,0) (i+1 ,j-1) (i+1 ,j) (i+1 ,j+1 ) (i+1 ,k') Pac_sp(i) (i3O) (i,j-1) (i,j) (i,j+1) (i,k') Pac_sp(i1) (i-1 ,0) (i-1 ,j-1) (i-1 ,j) (i-1 ,j+1 ) (i-1 ,k') ... Pac_sp(0) (0,0) (0,j-1) (0,j) (0,j+1) (0,k') Tableau 1 15 SOC(0) .. SOC(j-1) SOC(j) SOC(j+1) .. SOC(k') Pac_sp(k) ASOC(k,O) ASOC(k,j-1) ASOC(k,j) ASOC(k,j+1) ASOC(k,k') - - - Pac_sp(i+1) ASOC(i+1 ,0) ASOC(i+1 ,j-1) ASOC(i+1 ,j) ASOC(i+1 ,j+1 ) ASOC(i+1 ,k') Pac_sp(i) ASOC(i3O) ASOC(i,j-1) ASOC(i,j) ASOC(i,j+1) ASOC(i,k') Pac_sp(i-1) ASOC(i-1 ,0) ASOC(i-1 ,j-1) ASOC(i-1 ,j) ASOC(i-1,j+1) ASOC(i-1 ,k') - - - Pac_sp(0) ASOC(0,0) ASOC(0,j-1) ASOC(0,j) ASOC(0,j+1) ASOC(0,k') Tableau 2 Les figures 6A et 6B illustrent en 3D la table de la variation de l'état de charge SOCV_TC du système (selon l'axe z) en fonction de la puissance de consigne Pac_sp et de l'état de charge SOC. On voit clairement apparaitre sur ces deux figures deux zones séparées par deux lignes noires qui correspondent aux limites imposées par 5 la puissance minimale et la puissance maximale, ces valeurs limites variant en fonction de l'état de charge SOC comme cela sera décrit plus loin. Cette représentation permet de faire apparaitre un aspect important concernant la table de la variation de l'état de charge SOCV_TC du système : elle permet d'attribuer une valeur à la variation temporelle de l'état de charge cl÷tc pour des puissances de 10 consigne Pac_sp supérieures (en valeur absolue) aux valeurs limites que sont la puissance maximale et la puissance minimale pour un état de charge donné. Ceci n'est rendu possible que par le recours à l'interpolation comme cela vient d'être décrit ou bien encore en utilisant deux autres variantes dites « sans interpolation » ou « avec la valeur précédente » que nous allons décrire maintenant. Il est à noter 15 que l'objectif de cette interpolation n'est pas de renseigner la table pour des valeurs qui dans tous les cas ne seront pas utilisées car non atteignables, mais de rendre plus précis les points à la limite des zones précédemment décrites. 16 Dans un mode de réalisation alternatif dit « sans interpoltaEstsion », chaque valeur de la table SOCV TC de correspondance est égale àASOC(i,j), avec ASOC(i,j) la variation de charge obtenue lors de la simulation correspondant à l'état de charge SOC(j) et à une puissance de consigne Pac_sp(i). En effet, les inventeurs ont mis en 5 évidence que, de manière surprenante, cette approximation permettait d'obtenir de très bons résultats pour un coût moindre en termes de puissance de calcul (cet aspect de l'invention sera illustré dans la suite). Dans un mode de réalisation alternatif dit « avec la valeur précédente », chaque valeur de la table de la variation de l'état de charge SOCV_TC du système est égale 10 à: - la valeur de la variation de l'état de charge ASOC obtenue lors de ladite simulation sur le temps tESS si la puissance moyenne au cours de ladite simulation est égale à la puissance de consigne Pac_sp. Sd SOC(m,j) d OC(m,j) dt dt - dans le cas contraire, ou est la variation de charge obtenue 15 lors de la simulation correspondant à l'état de charge SOC(j) et la puissance de consigne Pac_sp(M), Pac_sp(M) étant la puissance de consigne la plus proche de la puissance de consigne Pac_sp(i) pour laquelle la puissance moyenne au cours de ladite simulation est égale à la puissance de consigne Pac_sp(M). 20 Ainsi, lorsque la puissance de consigne Pac_sp ne peut pas être maintenue lors de la simulation, les résultats d'une simulation précédente sont utilisés afin d'attribuer dsdot c une valeur à la variation temporelle de l'état de charge - du système de stockage de l'énergie ESS. Cet aspect permet de limiter les calculs nécessaires à l'obtention dSOC de la table SOCV_TC de la variation temporelle de l'état de charge -dt du système 25 avec cependant certaines limitations abordées dans la suite et illustrées à la figure 9. Dans un mode de réalisation illustré aux figures 7A et 7B, afin d'obtenir la table PAC _ MAX_ TC de puissance maximale en fonction de l'état de charge SOC, la pluralité de simulations comporte une simulation effectuée avec un état de charge 17 initial SOCini égal à l'état de charge maximum SOCmax, une durée égale à la durée nécessaire au déchargement complet du système de stockage tDch (c'est-à-dire jusqu'à que l'état de charge atteigne la valeur de l'état de charge minimum SOCmin) et une puissance de consigne P ac_sp infinie positive. On entend par une puissance 5 de consigne P ac_sp infinie positive une puissance de consigne P ac_sp largement supérieure à la puissance que le système peut fournir. Ainsi, la puissance moyenne k sur un intervalle donné peut être considérée comme la puissance maximale pour l'état de charge SOCk correspondant audit intervalle. De plus, le calcul de la table PAC MAX TC de puissance maximale en fonction de l'état de charge SOC 10 comporte : - une étape de détermination de périodes de durée tESS au sein de ladite simulation, chaque période étant identifiée au moyen d'un entier positif k; - pour chacune de ces périodes, une étape de détermination de l'état de charge SOCk au début de la période k et de la puissance 15 moyenne k durant la période k. La pluralité de couples (SOCk, k) vient alors constituer la table PAC_MAX_TC de puissance maximale en fonction de l'état de charge SOC de sorte que pour chaque état de charge SOCk donné, il est possible d'attribuer une puissance maximale (égale à k). 20 Dans un mode de réalisation illustré aux figures 8A et 8B, afin d'obtenir la table PAC _ MIN _TC de puissance minimale en fonction de l'état de charge SOC, la pluralité de simulations comporte une simulation effectuée avec un état de charge initial SOC ini égal à l'état de charge minimale SOCmin, une durée égale à la durée nécessaire au chargement complet du système de stockage tch (c'est-à-dire jusqu'à 25 que l'état de charge atteigne la valeur de l'état de charge maximum SOCmax) et une puissance de consigne P ac_sp infinie négative. On entend par une puissance de consigne P ac_sp infinie négative une puissance de consigne P ac_sp largement supérieure à la puissance que le système peut accepter lors de la charge. Ainsi, la puissance moyenne k' sur un intervalle donné peut être considérée comme la 30 puissance minimale pour l'état de charge SOCk' correspondant audit intervalle. De plus, le calcul de la table PAC_MIN_TC de puissance minimale en fonction de l'état 18 de charge SOC comporte : - une étape de détermination de périodes de durée tESS au sein de ladite simulation, chaque période étant identifiée au moyen d'un entier positif k' ; - pour chacune de ces périodes, une étape de détermination de l'état de 5 charge SOCK au début de la période k' et de la puissance moyenne k' durant la période k'. La pluralité de couples (SOCK, k') vient alors constituer la table PAC_MIN_TC de puissance minimale en fonction de l'état de charge SOC de sorte que pour chaque état de charge SOCk' donné, il est possible d'attribuer une puissance 10 minimale (égale à k'). Dans un mode de réalisation, chaque simulation d'un premier sous-ensemble de simulations est effectué pour un état de charge SOC(j) et une puissance de consigne ac_spv, P M donnés et chaque état de charge SOC(j) est séparé du précédent . SOC(j-1) et/ou du suivant SOC(j+1) par un pas d'état de charge adaptatif et/ou 15 chaque puissance de consigne ac_spv,M P est séparée de la précédente Pac_sp(i-1) . et/ou de la puissance suivante ac_spv+ P (i 1 ) d'un pas de puissance de consigne . adaptatif. Dans un mode de réalisation, l'étape de simulation est répétée pour une pluralité de durées tESS. Ainsi, l'on constitue une table de correspondance par durée de 20 simulation tESS, la table de correspondance utilisée par le système de pilotage étant choisie en fonction de la précision nécessaire. La figure 9 permet de mettre en avant les avantages du modèle obtenu à l'aide d'un procédé selon un premier aspect de la présente invention. Le premier histogramme HMC correspond au temps de simulation lorsque la simulation est effectuée à l'aide 25 d'un système complexe et avec un pas de temps At égal à 1 seconde. Le graphique comprend également cinq groupements H1, H2, H3, H4, H5 de quatre histogrammes, chaque groupement correspondant à une simulation identique à celle effectuée avec le modèle complexe, mais effectuée à l'aide d'un modèle simplifié obtenu par un procédé selon un premier aspect de l'invention et pour un 30 temps de simulation tESS donné. Autrement dit, chaque groupement correspond à 19 une simulation effectuée avec un modèle simplifié et avec un pas de temps de simulation égal au temps de simulation tEss utilisé lors du procédé de détermination des paramètres dudit modèle simplifié. Dans chaque groupement, le premier histogramme concerne le temps de simulation 5 et les deuxième, troisième et quatrième histogrammes concernent : - pour le deuxième histogramme, l'erreur RMSE entre la simulation utilisant le modèle complexe et une simulation utilisant un modèle obtenu à l'aide d'un procédé selon un premier aspect de l'invention dans le mode de réalisation dit « avec interpolation » ; 10 - pour le troisième histogramme, l'erreur RMSE entre la simulation utilisant le modèle complexe et une simulation utilisant un modèle obtenu à l'aide d'un procédé selon un premier aspect de l'invention dans le mode de réalisation dit « sans interpolation » ; - pour le quatrième histogramme, l'erreur RMSE entre la simulation utilisant le 15 modèle complexe et une simulation utilisant un modèle obtenu à l'aide d'un procédé selon un premier aspect de l'invention dans le mode de réalisation dit « avec la valeur précédente ». Le groupement H1 concerne un cas où le temps tESS est égal à 926 secondes. Dans ce cas de figure, la précision obtenue avec les méthodes « avec interpolation » et 20 « sans interpolation » est sensiblement identique tandis que la précision obtenue avec la méthode « avec la valeur précédente » est inférieure. De plus, le temps de simulation à l'aide d'un modèle obtenu par un procédé selon l'invention est largement inférieur au temps de simulation à l'aide d'un modèle complexe selon l'état de l'art antérieur. 25 Le groupement H2 concerne un cas où le temps tESS est égal à 463 secondes. On constate que la précision des méthodes « avec interpolation » et « sans interpolation » est améliorée de manière significative alors que la précision avec la méthode « avec la valeur précédente » s'est nettement dégradée. Les groupements H3 concernant un cas où le temps tESS est égal à 93 secondes, 30 H4 concernant un cas où le temps tESS est égal à 47 secondes et H5 concernant un 20 cas où le temps tESS est égal à 10 secondes ne font que confirmer les tendances observées dans les deux groupements précédents à savoir que l'erreur diminue lorsque le temps tESS diminue et que les méthodes « avec interpolation » et « sans interpolation » permettent d'obtenir des résultats similaires tandis que la méthode « 5 avec la valeur précédente » est toujours moins performante. Afin d'illustrer le déroulé du procédé selon un premier aspect de l'invention, nous allons dans la suite présenter un modèle complexe comprenant un modèle relatif à un dispositif de conversion ainsi qu'un modèle relatif à un dispositif de stockage. Il est important de noter que le procédé selon un premier aspect de l'invention ne 10 dépend pas du type de modèle complexe utilisé et le modèle complexe qui va suivre est donné à titre purement illustratif. Il permettra de faire apparaitre l'avantage en termes de simplification de la modélisation (et de donc du pilotage) d'un système de stockage de l'énergie que procure le modèle obtenu à l'aide d'un procédé selon un premier aspect de l'invention. 15 Dans la suite, les valeurs d'entrée sont initialisées lors de la première itération puis mises à jour lors des itérations suivantes. Les valeurs de sortie sont quant à elles calculées en fonction des valeurs d'entrée ainsi que des paramètres du modèle. Les paramètres du modèle sont, en règle générale, fournis par le constructeur du dispositif de conversion ou du dispositif de stockage, mais peuvent également être 20 déterminés à partir de tests et de mesures expérimentales. Par ailleurs, on notera Xt la valeur de la grandeur X à l'itération t. De plus, la fonction liée à une table de correspondance Y sera notée fy. Il est également supposé que lorsqu'une valeur n'est pas directement disponible dans une table de correspondance, cette dernière est obtenue par interpolation, par exemple interpolation linéaire, à partir des valeurs 25 disponibles dans ladite table de correspondance (il s'agit d'une méthode standard d'utilisation d'une table de correspondance). L'utilisation du modèle relatif à un dispositif de conversion va maintenant être détaillée. Les entrées, paramètres et sorties relatives audit modèle seront introduits au fur et à mesure. Le modèle permet tout d'abord de calculer l'angle de déphasage 30 41) entre la partie active P . ac_sp et la partie réactive n -ac_sp de la puissance de consigne. 21 La partie active P . ac_sp et la partie réactive n -ac_sp de la puissance de consignes sont des entrées du modèle. L'angle de déphase s'obtient simplement à l'aide de la relation suivante : Qac = tan-'( ) Pac sp 5 Toutefois, afin de simplifier l'exposé qui va suivre, la puissance réactive de consigne Qac_sp, est considérée nulle sur l'ensemble de la simulation et donc la valeur de l'angle de déphasage 0 ne dépend pas de l'itération. Une fois l'angle de déphase connu, le modèle permet de calculer la puissance maximale apparente Smax_LUT à l'aide d'une table de correspondance Smax_LUT 10 fournissant la puissance maximale apparente Smax en fonction de l'angle de déphasage 0 de sorte que : Smax - fsmaxLuT (cl)) Le modèle permet ensuite de calculer les puissances active P . ac_sat et réactive Qac_stat maximales que le système de conversion peut absorber. Ce calcul s'effectue 15 à l'aide d'une fonction de saturation fonction fsat qui modélise la saturation appliquée aux consignes de puissances active et réactive. Autrement dit : [Pac sat,Qac satl = fsat(Smax,d),Pac sp,Qac_sp) Comme mentionné précédemment, la puissance réactive de consigne Qac_sp est considérée nulle sur l'ensemble de la simulation et donc, il en est de même de 20 Qac_sat. La relation précédente peut donc être réécrite de la manière suivante : Pac sat = min(Smax, Pac sp) Cette condition permet de s'assurer que même si la puissance de consigne P . ac_sp est supérieure à la puissance maximale Smax pouvant être fournie par le dispositif de conversion, la puissance réellement fournie est limitée par la puissance 25 maximale que le dispositif de conversion peut fournir. 22 Le modèle permet ensuite de calculer la puissance continue Pdc correspondant à la puissance de consigne P ac_sp à l'aide de la table de correspondance Pac_LUT donnant la puissance alternative Pac en fonction de la puissance continue Pdc, la puissance réactive Qac, de la tension Udc reçue en entrée et de la tension Uac ; 5 autrement dit : Pac = fPac LUT (P dc, Q ac' Udc, Uac) Il convient ensuite d'utiliser la relation précédente afin de déterminer Pdc. Pour cela, les valeurs de Pac et Qac sont choisies de sorte que Pac = Pac_sp et Qac = Qac_sp, la tension Uac étant connue et égale à la tension nominale du système de 10 conversion Unom fournie en entrée, et la méthode de Newton-Raphson est utilisée pour résoudre la minimisation suivante : Min(lPac_sp fPac LUT (Pdc, Qac sp, U dc,U ac)I OÙ Pdc est la valeur initiale prise égale à P ac_sp à partir de laquelle est effectué la minimisation, ladite minimisation permettant d'obtenir la puissance continue Pdc. 15 Le modèle permet ensuite d'envoyer le courant de consigne de la batterie Idc_sp au modèle du dispositif de stockage après que ce dernier est été calculé à l'aide de la relation suivante : Pdc dc sp = U dc Le modèle permet par ailleurs de calculer la puissance active Pac à l'aide de la table 20 de correspondance Pac_LUT introduite auparavant de sorte que : Pac = fPac LUT (P dc, Qac' Udc, Uac) Comme mentionné précédemment, la puissance réactive de consigne Qac_sp est considérée nulle sur l'ensemble de la simulation et donc, il en est de même de Qac. Il est également important de noter que la tension de sortie Uac est imposée par le 25 système de sorte à ne pas avoir une tension fluctuante et sa valeur est donnée par 23 la tension alternative nominale du système de conversion Unom, cette dernière étant un paramètre du modèle. Dans ce même exemple de réalisation, on considère un modèle relatif à un dispositif de stockage prenant en entrées le courant de consigne Idc_sp (fourni par le dispositif 5 de conversion) et la température. Ce modèle est couplé au modèle détaillé précédemment pour constituer le modèle complexe du système de stockage de l'énergie ESS. Dans la suite, on considère que la température du dispositif de stockage est fournie par un modèle externe au modèle complexe et fournie en entrée de ce dernier. De plus, l'hypothèse est faite que la température n'influe pas 10 sur le comportement de la batterie ou que la température est maintenue par un dispositif externe. Dans ce cas la température en entrée reste constante et égale à la température initiale. Il ne s'agit cependant que d'un exemple de réalisation et il est également envisageable de coupler le modèle complexe avec un modèle thermique lorsque les hypothèses précédentes ne sont plus respectées. 15 Le modèle permet tout d'abord de calculer la puissance perdue dans le système auxiliaire à l'aide d'une table de correspondant dc P_aux_ LUT fournissant la puissance . perdue Pac_aux dans le système auxiliaire du dispositif de stockage en fonction du courant de consigne Idc_sp reçu par le dispositif de stockage de sorte que : PL aux = fPdc auxLUT (1f1c sp) 20 En considérant que les éléments du système auxiliaire (les ventilateurs par exemple) fonctionnent à leur tension optimale connue de l'utilisateur et égale à Udc_aux, le courant consommé par ledit système auxiliaire, le modèle permet d'obtenir le courant auxiliaire Idc_aux à l'aide de la relation suivante : pt It - dc aux dc aux - `-' dc aux 25 Une fois le courant consommé par ledit système auxiliaire Idc_aux connu, le modèle permet d'évaluer l'influence de ce dernier sur le courant de consigne Idc_sp afin d'obtenir un courant de consigne modifié lac' à l'aide de la relation suivante : 24 t c = sp 1(ti _c. aux Autrement dit, lorsque le dispositif de stockage est chargé avec un courant Idc_sp, le courant servant réellement à ladite charge est égal à lac'. De la même manière, lorsque le dispositif de stockage est déchargé par un courant Idc_sp, le courant que 5 le système de stockage peut réellement fournir est égal à lac'. Le modèle permet également de prendre en compte le fait que le courant lac' est limité par un courant maximum Ichmax lors de la charge du système de stockage et un courant maximum IDchmax lors de la décharge dudit système. Pour cela, le courant maximum Ichmax lors de la charge est déterminé à l'aide d'une table de 10 correspondance Ichmax_LUT fournissant le courant maximum Ichmax lors de la charge en fonction de l'état de charge SOC, de l'état de santé SOH et de la température T : 1=ax fichmax_LuT(S0Ct,5'0Ht,r) De même, le courant maximum IDchmax lors de la décharge est obtenu à l'aide de table de correspondance IDchmax_LUT fournissant le courant maximum IDchmax lors 15 de la décharge en fonction de l'état de charge SOC, de l'état de santé SOH et de la température T: fiDchmax_LUT(S°Ct,S0Ht,Tt) Autrement dit, le courant lac peut être déterminé à l'aide de la relation suivante : ct Si lihMax REVENDICATIONS1. Procédé (100) de détermination des paramètres d'un modèle simplifié d'un système de stockage de l'énergie (ESS), ledit système (ESS) comprenant un dispositif de stockage de l'énergie (DSE) et un dispositif de conversion (DC), ledit système (ESS) pouvant être modélisé au moyen d'un modèle complexe (MC) incluant un modèle du dispositif de stockage de l'énergie (MCS) et un modèle du dispositif de conversion (MCC) ; ledit modèle complexe (MC) recevant en entrée une puissance de consigne P . ac_sp et un état de charge SOC, et fournissant en sortie l'état de charge SOC du dispositif de stockage ainsi que la puissance Pac en sortie du dispositif de stockage ; ledit procédé étant caractérisé en ce qu'il comprend : - une première étape (101) de mise en oeuvre d'une pluralité de simulations du système de stockage de l'énergie (ESS) à l'aide du modèle complexe (MC) ; - une deuxième étape (102) de calcul à partir des résultats obtenus lors de la première étape (101) : - d'une table de la variation temporelle de l'état de charge (SOCV_TC) du système en fonction de la puissance de consigne P . ac_sp et de l'état de charge SOC ; - d'une table de puissance maximale (PAC_MAX_TC) en fonction de l'état de charge SOC ; - d'une table de puissance minimale (PAC_MIN_TC) en fonction de l'état de charge SOC ; le modèle simplifié obtenu permettant d'attribuer en fonction d'une puissance de consigne P . ac_sp et d'un état de charge SOC p fournis en entrée, et à partir des tables (SOCV_TC, PAC_MAX_TC, PAC_MIN_TC) déterminées lors de la deuxième étape (102), une puissance Pac et un état de charge du système SOC. 2. Procédé (100) selon la revendication précédente caractérisé en ce que chaque valeur de la table de la variation de l'état de charge (SOCV_TC) du système est obtenue à l'aide d'une simulation d'une durée tESS effectuée 33 5 pour un état de charge SOC(j) et une puissance de consigne P M donnés ac_spv, et appartenant à un premier sous-ensemble de la pluralité de simulations, ladite valeur étant : - égale à la valeur de la variation de l'état de charge ASOC obtenue lors de ladite simulation sur le temps tESS si la puissance moyenne au cours de ladite simulation est égale à la puissance de consigne P . ac_sp. - dans le cas contraire, égale à: alac '(i 1), 1, rdsoc(i-i,j) ,âsoc(i,i) (n} dt tESS interp rac splt 10 avec : dSOC - la variation de charge obtenue lors de la simulation dt correspondant à l'état de charge SOC(j) et à un puissance de consigne Pac_sp(i-1) ; - ASOC(i, j), la variation de charge obtenue lors de la simulation correspondant à l'état de charge SOC(j) et à une puissance de consigne P . ac_spv, ; - Pac (i,j) >, la valeur moyenne de la puissance Pac lors de la simulation correspondant à l'état de charge SOC(j) et à une puissance de consigne P ac_sp ./. 20 3. Procédé (100) selon la revendication 1 caractérisé en ce que chaque valeur de la table de la variation de l'état de charge (SOCV_TC) du système est obtenue à l'aide d'une simulation d'une durée tESS effectuée pour un état de charge SOC(j) et une puissance de consigne ac_sp P M donnés et appartenant . " à un premier sous-ensemble de la pluralité de simulations, ladite valeur étant 25 égale à la valeur de la variation de l'état de charge ASOC obtenue lors de ladite simulation sur le temps tEss. 4. Procédé (100) selon la revendication 1 caractérisé en ce que chaque valeur de la table de la variation temporelle de l'état de charge (SOCV_TC) du système est obtenue à l'aide d'une simulation d'une durée tESS effectuée 30 pour un état de charge SOC(j) et une puissance de consigne P M donnés . ac_sp" 15 34 et appartenant à un premier sous-ensemble de la pluralité de simulations, ladite valeur étant : - égale à la valeur de la variation de l'état de charge ASOC obtenue lors de ladite simulation sur le temps tESS si la puissance moyenne 5 au cours de ladite simulation est égale à la puissance de consigne P ac_sp. d SSOC(m,j) d OC(m,j) dt dt - dans le cas contraire, égale à ou est la variation de charge obtenue lors de la simulation correspondant à l'état de charge SOC(j) et la puissance de consigne Pac_sp(M), Pac_sp(M) étant la 10 puissance de consigne la plus proche de la puissance de consigne Pac_sp(i) pour laquelle la puissance moyenne au cours de ladite simulation est égale à la puissance de consigne Pac_sp(M). 5. Procédé (100) selon l'une des trois revendications précédentes caractérisé en ce que chaque état de charge SOC(j) est séparé du précédent SOC(j-1) 15 et/ou du suivant SOC(j+1) par un pas d'état de charge adaptatif et/ou chaque puissance de consigne Pac_sp(i)- est séparée de la précédente Pac_sp, 1, et/ou de la suivante Pac_spv (i+ 1 ) d'un pas de puissance de consigne adaptatif. 6. Procédé (100) selon l'une des revendications 3 à 5 caractérisé en ce que l'étape de simulation est répétée pour une pluralité de durées tESS. 20 7. Procédé (100) selon l'une des revendications précédentes caractérisé en ce que la pluralité de simulations comporte une simulation effectuée avec un état de charge initial SOC ini égal à l'état de charge maximum SOCmax, une durée égale à la durée nécessaire au déchargement complet du système de stockage tDch et une puissance de consigne Pac_sp infinie positive et en ce 25 que le calcul de la table PAC_MAX_TC de puissance maximale en fonction de l'état de charge SOC comporte : - une étape de détermination de périodes de durée tESS au sein de ladite simulation, chaque période étant identifiée au moyen d'un entier positif k; 30 - pour chacune de ces périodes, une étape de détermination de l'état de charge SOCk au début de la période k et de la puissance moyenne k durant la période k; la pluralité de couple (SOCk, k) venant constituer la table PAC MAX TC. 5 8. Procédé (100) selon l'une des revendications précédentes caractérisé en ce que la pluralité de simulations comporte une simulation effectuée avec un état de charge initial SOC ini égal à l'état de charge minimum SOCmin, une durée égale à la durée nécessaire au chargement complet du système de stockage tch et une puissance de consigne P ac_sp infinie négative et en ce 10 que le calcul de la table PAC_MIN_TC de puissance minimale en fonction de l'état de charge SOC comporte : - une étape de détermination de périodes de durée tESS au sein de ladite simulation, chaque période étant identifiée au moyen d'un entier positif k' ; 15 - pour chacune de ces périodes, une étape de détermination de l'état de charge SOCk' au début de la période k' et de la puissance moyenne k' durant la période k' ; la pluralité de couple (SOCk', k) venant constituer la table PAC MIN TC. 20 9. Procédé de pilotage d'un système de stockage de l'énergie (ESS) calculant des consignes de fonctionnement dudit système (ESS) à partir d'un modèle dudit système (ESS) caractérisé en ce que ledit modèle est obtenu à l'aide d'un procédé (100) selon l'une des revendications précédentes. 10. Dispositif de pilotage d'un système de stockage de l'énergie (ESS) 25 comprenant des moyens pour envoyer des consignes de fonctionnement au système de stockage de l'énergie (ESS), des moyens pour recevoir des données concernant le fonctionnement du système de stockage de l'énergie (ESS) caractérisé en ce qu'il comporte des moyens pour mettre en oeuvre un procédé de régulation selon la revendication précédente. 30