FR 3059845 A1 20180608 FR 1661779 A 20161201 1 PROCÉDÉ D'EQUILIBRAGE D'ETATS DE CHARGE D'UNE PLURALITE DE DISPOSITIFS DE STOCKAGE D'ENERGIE ELECTRIQUE DOMAINE TECHNIQUE DE L'INVENTION La présente invention concerne un procédé d'équilibrage d'état de charge d'une pluralité de dispositifs de stockage d'énergie électrique appartenant à un système de stockage d'énergie électrique. i.o ARRIÈRE-PLAN TECHNOLOGIQUE DE L'INVENTION Les dispositifs de production d'énergie électrique basés sur les énergies renouvelables produisent une quantité d'énergie qui dépend directement de leur environnement. Dans le cas par exemple d'une centrale éolienne, plus il y a de vent, 15 plus la centrale produit d'énergie électrique, alors qu'en l'absence de vent, la centrale ne produit pas d'énergie. Un tel dispositif ne peut pas adapter sa production à la demande en énergie d'un réseau électrique auquel il est connecté. Par conséquent, ce dispositif de 20 production nécessite un système de stockage permettant d'emmagasiner l'énergie produite qui n'est pas immédiatement consommée et de la restituer ultérieurement afin de répondre à des sollicitations du réseau électrique. De manière classique, un système de stockage d'énergie comporte plusieurs 25 batteries connectées en parallèle et/ou en série. Initialement, les batteries présentent un même état de charge. Au cours du temps, les états de charge des batteries dérivent les uns par rapport aux autres, par exemple à cause de sollicitations qui sont différentes d'une batterie à l'autre suite à l'application d'une gestion d'énergie basique, ou parce que les batteries ne vieillissent pas de façon homogène. Le 30 système de stockage est alors limité soit par la batterie qui a l'état de charge le plus élevé lorsque le système est en charge, soit par la batterie qui a l'état de charge le plus faible lorsque le système est en décharge. Les capacités du système de stockage sont tirées vers le bas. 3059845 2 La disparité des états de charge accentue davantage l'inhomogénéité du vieillissement des batteries, ce qui se répercute sur les états de charges, et ainsi de suite. Par ailleurs, la limitation du système de stockage conduit à une perte d'un réservoir énergétique non négligeable. Pour tenter de pallier ces problèmes, le 5 système de stockage est généralement surdimensionné, ce qui augmente inutilement le coût du système. Actuellement, il existe des algorithmes de gestion de tels systèmes qui assurent la mise en sécurité des batteries par exemple en déconnectant une batterie endommagée ou dont la température est trop élevée. Cependant ces algorithmes de gestion ne permettent pas d'assurer un équilibrage des états de charge entre les différentes batteries que comporte le système de stockage. Il apparaît donc un besoin d'améliorer la gestion d'énergie du système de 15 stockage pour permettre d'homogénéiser dans le temps les états de charge des batteries. RÉSUMÉ DE L'INVENTION 20 La présente invention vise à équilibrer des états de charge entre plusieurs dispositifs de stockage d'énergie appartenant à un même système de stockage d'énergie. Selon un premier aspect de l'invention, ce but est atteint en prévoyant un 25 procédé d'équilibrage d'états de charge d'une pluralité de dispositifs de stockage d'énergie électrique formant un système de stockage d'énergie, les dispositifs de stockage étant connectés à un bus de distribution par l'intermédiaire d'éléments de connexion, le système de stockage présentant un courant de fonctionnement traversant le bus de distribution. Le courant de fonctionnement appartient à une 30 plage de courant s'étendant entre un courant de charge maximal, de signe négatif, et un courant de décharge maximal, de signe positif. Cette plage de courant est divisée en au moins deux intervalles de courant : un premier intervalle s'étendant, en valeurs absolues, entre le courant de charge maximal et un courant nul (phase de charge), et un deuxième intervalle s'étendant entre un courant nul et le courant de décharge 3059845 3 maximal (phase de décharge). Chaque élément de connexion présente au moins deux états de conduction électrique. Le procédé d'équilibrage comporte les étapes suivantes : - (a) obtenir l'état de charge des dispositifs de stockage par des moyens de 5 mesures de l'état de charge ; - (b) calculer, à partir des états de charge obtenus au cours de l'étape (a), une cartographie des états de conduction des éléments de connexion associée à chaque intervalle de courant ; - (c) déterminer l'intervalle de courant auquel appartient le courant de 10 fonctionnement du système de stockage ; - (d) appliquer aux éléments de connexion la cartographie associée à l'intervalle de courant déterminé au cours de l'étape (c) ; les étapes (a) et (b) étant répétées périodiquement pour mettre à jour les cartographies des états de conduction des éléments de connexion. On entend par « cartographie des états de conduction » une représentation indiquant dans quel état de conduction placer chaque élément de connexion. Par conséquent, lorsque la cartographie est appliquée aux éléments de connexion, chacun d'entre eux se place dans l'état de conduction qui lui est attribué dans la cartographie, en fonction d'un courant de sollicitation auquel est soumis le système de stockage d'énergie. Le procédé de fabrication selon l'invention peut également comporter une ou plusieurs caractéristiques parmi les suivantes, considérées individuellement ou selon les combinaisons techniquement possibles : - l'étape (b) de calcul des cartographies comporte les opérations suivantes : o construire un problème d'optimisation défini par des variables de décision, des contraintes, et une fonction objectif visant à minimiser l'écart entre les états de charge des dispositifs de stockage ; o résoudre le problème d'optimisation ; les dispositifs de stockage sont connectés en parallèle au bus de distribution par l'intermédiaire de contacteurs, le problème d'optimisation étant non-linéaire ; 3059845 4 le problème d'optimisation est linéarisé afin de garantir l'obtention d'une solution optimale avec un temps de calcul réduit ; les dispositifs de stockage sont connectés en parallèle au bus de distribution par l'intermédiaire de résistances variables, le problème d'optimisation étant 5 non-linéaire ; dans ce cas, le problème d'optimisation ne peut pas être linéarisé et la cartographie sera alors sous optimale ; le procédé comporte une étape préalable de caractérisation d'au moins un dispositif de stockage pour déterminer, en fonction de l'état de charge du dispositif de stockage, les variations de chacun des paramètres 10 caractéristiques du dispositif de stockage suivants : la tension de circuit ouvert, la résistance interne, le courant de charge maximal et le courant de décharge maximum ; le premier intervalle de courant et/ou le deuxième intervalle de courant sont subdivisés en un nombre prédéterminé d'intervalles de courant ; 15 le nombre prédéterminé d'intervalles de courant est paramétré par un utilisateur. Un deuxième aspect de l'invention concerne un produit programme d'ordinateur comprenant des instructions pour la mise en oeuvre d'un procédé selon 20 le premier aspect de l'invention, lorsque ces instructions sont exécutées par un processeur. BRÈVE DESCRIPTION DES FIGURES L'invention et ses différentes applications seront mieux comprises à la lecture 25 de la description qui suit et à l'examen des figures qui l'accompagnent, parmi lesquelles : - la figure 1 représente schématiquement un premier exemple de contexte de mise en oeuvre du procédé selon l'invention ; - la figure 2 représente un diagramme fonctionnel d'un mode de mise oeuvre du 30 procédé selon l'invention ; - les figures 3A, 3B et 3C représentent pour les modes de charge et de décharge, en fonction de l'état de charge de la batterie, les variations respectivement 3059845 5 de la tension de circuit ouvert, de la résistance interne, et du courant de fonctionnement maximal ; - les figures 4A et 4B représentent respectivement un exemple de profil de sollicitation en courant appliqué au système de stockage de la figure 1 et à un 5 système de stockage de l'art antérieur, et les états de charges initiaux de leurs batteries; - la figure 5 représente les variations au cours du temps des états de charge minimaux et maximaux du système de stockage de la figure 1 et du système de stockage de l'art antérieur lorsque le profil de la figure 4A est appliqué ; 10 - la figure 6A représente les variations au cours du temps des courants minimaux et maximaux fournis ou absorbés par les batteries du système de stockage de la figure 1 et du système de stockage de l'art antérieur lorsque le profil de la figure 4A est appliqué ; - la figure 6B représente les variations au cours du temps de la différence entre 15 les courants minimaux et maximaux de la figure 6A, pour le système de stockage de la figure 1 et pour le système de stockage de l'art antérieur ; - les figures 7A et 7B représentent les variations au cours du temps du nombre de batteries connectées respectivement dans le système de stockage de la figure 1 et dans le système de stockage de l'art antérieur lorsque le profil de la figure 4A est 20 appliqué ; - la figure 8 représente schématiquement un deuxième exemple de contexte de mise en oeuvre du procédé selon l'invention. Les figures ne sont présentées qu'à titre indicatif et nullement limitatif de l'invention. 25 Pour plus de clarté, les éléments identiques ou similaires sont repérés par des signes de référence identiques sur toutes les figures. DESCRIPTION DÉTAILLÉE DE MODES DE MISE EN OEUVRE DE L'INVENTION 30 La figure 1 représente schématiquement un premier exemple de contexte de mise en oeuvre du procédé selon l'invention. La figure 1 montre un système de stockage d'énergie électrique 100 comportant N dispositifs de stockage d'énergie 120 montés en parallèle. Chaque dispositif de stockage 120 est connecté électriquement à un bus de distribution 140 par l'intermédiaire d'un élément de 3059845 6 connexion 160. Le bus de distribution 140 est connecté électriquement à un réseau électrique (non représenté). Le réseau électrique comporte un ou plusieurs dispositifs de production 5 d'énergie électrique permettant de fournir de l'énergie au système de stockage 100, tels que des centrales photovoltaïques ou éoliennes, ou des générateurs diesel. Le réseau électrique comporte également un ou plusieurs consommateurs d'énergie électrique. 10 Le système de stockage 100 présente deux modes de fonctionnement : un mode de charge et un mode de décharge. Dans le mode de charge, le système de stockage 100 emmagasine l'énergie électrique fournie par le réseau électrique. Inversement, dans le mode de décharge, le système de stockage 100 restitue au réseau électrique l'énergie emmagasinée pour alimenter les consommateurs en 15 énergie électrique. Le système de stockage 100 présente une tension de fonctionnement prise sur le bus de distribution 140, et un courant de fonctionnement traversant le bus de distribution 140. Le courant de fonctionnement est un courant de charge lorsque le 20 système de stockage 100 est dans le mode de charge. Inversement, le courant de fonctionnement est un courant de décharge lorsque le système de stockage 100 est dans le mode de décharge. Par convention, le courant de charge est de signe négatif et le courant de décharge est de signe positif. 25 La connexion entre le bus de distribution 140 et le réseau électrique est effectuée par l'intermédiaire d'un convertisseur tension alternative - tension continue 180 qui a pour fonction d'adapter la forme et l'amplitude de la tension de fonctionnement. 30 Le procédé selon l'invention peut s'appliquer à tous les types de dispositifs de stockage d'énergie tels que des batteries électrochimiques ou des piles à combustible. Dans le système de stockage 100 de la figure 1, les dispositifs de stockage 120 sont des batteries électrochimiques identiques, et les éléments de connexions 160 sont des contacteurs 160. 3059845 7 Le contacteur 160 peut être placé dans deux états de conduction électrique : un état dit « ON » lorsque le contacteur 160 est fermé et un état dit « OFF » lorsque le contacteur 160 est ouvert. Lorsque le contacteur est à l'état ON, le contacteur 160 5 conduit l'électricité, et lorsque le contacteur est à l'état OFF, le contacteur 160 ne conduit pas l'électricité. La présente invention a pour objectif de piloter les contacteurs 160 de manière à ce que le système de stockage 100 réponde aux sollicitations du réseau électrique 10 tout en faisant en sorte que les batteries 120 conservent toutes un état de charge sensiblement identique au cours du temps. Les contacteurs 160 sont pilotés en utilisant une cartographie qui indique dans quel état de conduction (ON ou OFF) chaque contacteur 160 doit être placé en fonction du courant de sollicitation. 15 Le courant de fonctionnement du système de stockage 100 appartient à une plage de courant qui s'étend entre le courant de charge maximal 'climax (de signe négatif) et le courant de décharge maximal Idismax (de signe positif) du système de stockage 100. Le courant de fonctionnement est borné par ces courants 'climax et Idismax qui sont des paramètres intrinsèques du système de stockage 100. 20 La plage de courant est divisée en au moins deux intervalles de courant, un premier intervalle s'étendant entre le courant de charge maximal Ichmax et un courant nul, et un deuxième intervalle s'étendant entre le courant nul et le courant de décharge maximal Idismax. Ainsi, les modes de charge et de décharges sont dissociés 25 pour permettre le calcul de la cartographie des états de conductions des contacteurs 160 pour chaque intervalle de courant, comme décrit plus en détail par la suite. Une première cartographie est calculée pour le premier intervalle correspondant au mode de charge, et une deuxième cartographie est calculée pour le deuxième intervalle correspondant au mode de décharge. 30 Avantageusement, le premier intervalle de courant et le deuxième intervalle de courant peuvent être également divisés. Les intervalles de courant ainsi formés ont une taille plus réduite ce qui permet d'obtenir une cartographie plus adaptée pour chaque intervalle. Plus le nombre d'intervalles de courant est élevé, plus l'équilibrage 3059845 8 des états de charge des batteries 120 est efficace et converge rapidement. Le nombre et les bornes des intervalles de courant sont par exemple paramétrés par un utilisateur. 5 La figure 2 représente un diagramme fonctionnel d'un mode de mise en oeuvre du procédé selon l'invention. Au cours d'une étape (a), l'état de charge de chaque batterie 120 est obtenu par des moyens de mesures associés à chaque batterie 120. Ensuite, à partir de ces états de charges, la cartographie des états de conduction des contacteurs 160 associée à chaque intervalle de courant est 10 calculée. Pour le système de stockage 100 de la figure 1, le calcul de la cartographie est un problème d'optimisation de type combinatoire, non-linéaire et non-convexe, également appelé MINLP pour « Mixed Integer Non-Linear Programming en anglais. 15 De manière générale, un problème d'optimisation se définit par un ensemble de variables de décision prenant chacune leur valeur parmi un ensemble appelé domaine, et par un ensemble de relations logiques, appelées contraintes, qui restreignent les combinaisons de valeurs permises. Un tel problème se définit 20 également par une fonction objectif qui caractérise la qualité d'une solution. Une solution à un tel problème est une affectation d'une valeur de son domaine à chaque variable de décision qui satisfait toutes les contraintes. L'ensemble des solutions optimales, dans le cas d'un problème de minimisation, correspond aux solutions ayant la plus petite valeur de la fonction objectif parmi l'ensemble des solutions 25 existantes. Dans le contexte de l'invention, la fonction objectif consiste à équilibrer les états de charge SOC des N batteries 120 du système de stockage 100, ce qui revient à minimiser la différence entre les états de charge des batteries 120 prises 30 deux à deux. Cette différence peut se traduire mathématiquement par la matrice suivante : 3059845 9 /SOC, - SOC1 SOC1 - SOC2 SOC1 - SOCN_, SOC1 - SOCN \ SOC2 - SOC1 SOC2 - SOC2 SOC2 - SOCN_, SOC2 - SOCN SOC3 - SOC1 SOC3 - SOC2 SOC3 - SOCN_, SOC3 - SOCN \SOCN - soc, SOCN - SOC2 SOCN - SOCN_, SOCN - SOCNI Cette matrice est symétrique, par conséquent seule la partie triangulaire supérieure ou inférieure de la matrice est conservée, en éliminant également la diagonale de la matrice, afin de supprimer les redondances. La fonction objectif peut donc s'écrire sous les deux formes équivalentes suivantes, qui sont respectivement une 5 fonction quadratique et une fonction valeur absolue : Min( > (SOCi - SOCi)2 E Min( > abs(SOCi - SOCi) iEN i 0) ; 20 l'état de charge SOC de la batterie i, i E {i, ..., N} : SOC est supérieure ou égale à 0 (SOCi > 0) ; la tension de fonctionnement Upack du système de stockage 100 : Upack est supérieure ou égale à 0 (Upack > O). 3059845 10 Les contraintes du problème d'optimisation, qui représentent le fonctionnement du système de stockage 100 de la figure 1, sont par exemples définies par les relations suivantes : Umin Xi * eps (C. 10) 15 cette contrainte signifiant que si la batterie i, i E {1, ..., N}, est connectée, elle doit délivrer ou absorber un courant minimal eps, par exemple fixé à 0.0001 A, afin d'éviter la symétrie de la solution qui se traduit en pratique par un contacteur fermé et un courant délivré ou absorbé égal à 0. Le problème d'optimisation ainsi construit est non-linéaire car la contrainte 20 (C.7) et la fonction objectif sont non-linéaires. Un tel problème est difficile à résoudre et il n'est pas certain d'obtenir une solution optimale. Avantageusement, un problème linéaire est construit à partir du problème non-linéaire afin de simplifier la résolution et d'assurer que la solution trouvée soit optimale, sans perte de performance. 3059845 12 Pour effectuer cette linéarisation, des variables de décisions supplémentaires sont introduites : - une variable de linéarisation Li du produit des variables Ui et Xi, Li étant supérieure ou égale à 0 (Li > 0) ; 5 - une variable de linéarisation Hi du produit des variables Upack et Xi, Hi étant supérieure ou égale à 0 (Hi > 0) ; - deux variables de linéarisation et Yi de la fonction valeur absolue, YiT étant la partie positive et Yilétant la partie négative de la fonction valeur absolue, Yii+ et Yi étant supérieures ou égales à 0 (Yi-JE, Yii > 0) ; 10 - Uij qui vaut 1 si la partie positive de la valeur absolue est prise en compte, ou 0 si la partie négative de la valeur absolue est prise en compte (Uii E {0,1}) ; - une erreur en tension Dvi des batteries connectées qui ne doit pas dépasser l'erreur maximale autorisée paramétrée par l'utilisateur (en cas d'erreur de mesure du capteur de tension). 15 La fonction objectif (Obj.1) peut donc être linéarisée de la manière suivante : Min( (Yi-JE +Yij (Obj.1') iEN j SOCmax. Les contraintes (C.12) et (C.13) permettent d'obliger l'une des parties (positive ou 3059845 13 négative) de la valeur absolue à être égale à zéro. Ainsi, en combinant avec la contrainte (C.11), on obtient : Yi-JE = SOCi - SOCi si SOCi > SOCi ; et Yii = abs(SOCi - SOCi) si SOCi > SOCi Les contraintes (C.1) à (C.10) décrites précédemment restent valides à 5 l'exception de la contrainte non-linéaire (C.7) qui est linéarisée de la manière suivante : Hi - (Li + Dvi) = 0 (C. 14) -Deltav Umax. Upack - Hi Umax. 15 Des paramètres caractéristiques des batteries 120 telles que la tension de circuit ouvert OCVi, la résistance interne Rsi, le courant de charge maximal IchMaxi et le courant de décharge maximal IdisMaxi sont utilisés dans les contraintes (C.3) et (C.4) du problème d'optimisation. Avantageusement, on procède à une étape 3059845 14 préalable de caractérisation d'au moins une batterie parmi les N batteries 120 pour obtenir les variations de leurs paramètres caractéristiques en fonction de l'état de charge. Les batteries 120 étant identiques, leurs paramètres caractéristiques sont également sensiblement identiques. Il est donc possible de ne caractériser qu'une 5 seule batterie 120 afin de réduire la complexité d'intégration de l'invention. La figure 3A représente, en fonction de l'état de charge de la batterie SOC, les variations des tensions de circuit ouvert OCVch et OCVdis respectivement pour le mode de charge et pour le mode de décharge. La figure 3B représente, en fonction 10 de l'état de charge de la batterie SOC, les variations des résistances internes Rch et Rdis respectivement pour le mode de charge et pour le mode de décharge. La figure 3C représente, en fonction de l'état de charge de la batterie SOC, les variations des courants maximaux Ichmax et Idismax respectivement pour le mode de charge et pour le mode de décharge. Les tensions de circuit ouvert OCVch et OCVdis, les 15 résistances internes Rch et Rdis, et les courants maximaux Ichmax et Idismax sont exprimées en pourcentage de leurs valeurs maximales respectives. L'invention peut également s'appliquer à des batteries ayant un vieillissement différent, ce qui se traduit par des valeurs de résistances internes différentes d'une batterie à l'autre. Ces valeurs de résistances internes peuvent en effet être calculées dynamiquement 20 au moyen des informations de courant et de tension de la période précédente. La résolution du problème d'optimisation définit précédemment, par exemple au moyen d'un algorithme de recherche arborescente, également appelé « Branchand-Cut », permet d'obtenir les cartographies des états de conduction des 25 contacteurs 160. Après l'étape de calcul des cartographies, on procède à une étape (c) de détermination de l'intervalle de courant auquel appartient le courant de fonctionnement I req du système de stockage 100, c'est-à-dire le courant de charge 30 fourni au système 100 dans le mode de charge, ou le courant de décharge fourni par le système 100 dans le mode de décharge. Ensuite, la cartographie correspondant à l'intervalle de courant déterminé au cours de l'étape (c) est appliquée aux éléments de connexion au cours d'une étape 3059845 15 (d). En fonction de l'intervalle dans lequel se trouve le courant de fonctionnement, les contacteurs 160 ne sont pas placés dans les mêmes états de conduction. Grâce à l'invention, le système de stockage 100 répond aux besoins en courant tout en faisant en sorte que les états de charges des batteries 120 convergent vers un 5 même état de charge. Le procédé de gestion assure l'homogénéité des états de charge des batteries 120 au cours du temps. Les états de charge des batteries 120 variant en fonction des sollicitations appliquées au système de stockage 100, les étapes (a) d'obtention des états de 10 charges et (b) de calcul des cartographies sont avantageusement répétées régulièrement pour maintenir à jour les cartographies par rapport aux variations des états de charge des batteries 120. Une période d'actualisation Tupdate peut donc être définie. Cette période est par exemple paramétrée par l'utilisateur. 15 L'étape (c) de détermination de l'intervalle de courant peut être répétée à chaque fois que le courant de fonctionnement du système de stockage 100 est modifié. Alternativement, l'étape (c) peut être répétée périodiquement, avec une période appelée « pas temps » ou « pas de contrôle » qui correspond à la fréquence d'échantillonnage du courant de sollicitation, cette période étant inférieure à la 20 période d'actualisation Tupdate des étapes (a) et (b). La figure 4A montre un exemple de profil de sollicitation en courant P appliqué au système de stockage 100 et à un système de stockage de l'art antérieur 100' dans lequel une mise à niveau des états de charge est effectuée par une source 25 d'énergie externe en cas de nécessité. Pendant cette mise à niveau, le système de stockage 100' est indisponible et ne peut donc pas répondre à des sollicitations, ce qui n'est pas le cas avec le procédé selon l'invention. Le système de stockage de l'art antérieur 100' comporte des batteries 120 identiques à celles du système de stockage 100. Le profil de courant s'étend sur une durée de 24h avec un pas temps 30 d'une minute. La figure 4B montre les états de charges initiaux SOCinit, c'est-à-dire avant la mise en oeuvre du procédé d'équilibrage, des différentes batteries 120 des systèmes de stockage 100, 100'. L'abscisse de la figure 4B correspond à l'indice Idx des batteries 120. En l'occurrence, le nombre N de batteries est égale à 32. 3059845 16 Pour l'application du profil de la figure 4A au système de stockage 100, le procédé d'équilibrage est paramétré avec les valeurs suivantes : - la période d'actualisation Tupdate est égale à 15 mn ; - l'erreur maximale en tension Delta, est égale à 5V ; 5 - le nombre minimum de batteries 120 connectées Min_On_Battery au bus de distribution 140 est égale à la moitié des batteries, soit 16 sur 32 ; - le nombre d'intervalles de courant est égal à 5 ; - les états de charge minimal SOCmin et maximal SOCmax sont fixés respectivement à 10% et 100%. 10 Les résultats de l'application du profil de la figure 4A sont illustrés sur les figures 5, 6A, 6B, 7A et 7B. La figure 5 montre les variations au cours du temps T des états de charge 15 minimaux SOCminioo et SOCminioo, respectivement pour le système 100 et pour le système 100', et les variations au cours du temps T des états de charge maximaux SOCmaxioo et SOCmaxioo, respectivement pour le système 100 et pour le système 100'. Il est à noter que la différence entre l'état de charge minimal et l'état de charge maximal est bien inférieure pour le système de stockage 100 dans lequel le procédé 20 d'équilibrage est mis en oeuvre. A chaque instant, un état de charge moyen de chaque système 100, 100' peut être obtenu en calculant la moyenne des états de charge de chaque batterie du système 100, 100'. Grâce à l'invention, la différence entre l'état de charge moyen du 25 système 100 et l'état de charge moyen du système 100' est positive au cours du temps, ce qui représente un gain d'autonomie pour le système de stockage 100. A titre d'exemple, à la fin du profil de sollicitation de la figure 4A, l'état de charge moyen des batteries 120 du système de stockage 100 est d'environ 63,7%, 30 alors que l'état de charge moyen des batteries 120 du système de stockage de l'art antérieur 100' est seulement d'environ 61%, ce qui représente un gain d'autonomie d'environ 2,6% en faveur du système de stockage 100. De plus, sur toute la durée du 3059845 17 profil de sollicitation, le gain d'autonomie présente une valeur moyenne d'environ 2,8% et une valeur maximale d'environ 5,8%. La figure 6A montre les variations au cours du temps T des courants 5 minimaux Iminioo et Iminioo, fournis (mode de décharge) ou absorbés (mode de charge) respectivement par les batteries du système 100 et par les batteries du système 100', et les variations au cours du temps T des courants maximaux Imaxioo et Imaxioo, fournis ou absorbés respectivement par le système 100 et par le système 100'. La figure 6B montre les variations au cours du temps T de la différence 10 Delta_lioo entre le courant Imaxioo et le courant Iminioo, et les variations au cours du temps T de la différence Delta_1100, entre le courant Imaxioo, et le courant Iminiou. Ces figures permettent de constater que les courants fournis ou absorbés par les batteries du système 100 sont quasiment identiques, ce qui implique un vieillissement homogène des batteries. Au contraire, les différences entre les 15 courants fournis ou absorbés par les batteries du système 100' étant plus importantes, le vieillissement des batteries est moins homogène. Les figures 7A et 7B montrent les variations au cours du temps T du nombre de batteries 120 connectées Nlc respectivement pour le système 100 et pour le 20 système 100'. On constate que le nombre de batteries 120 connectées varie davantage dans le cas système de stockage 100 (entre 16 et 26 batteries connectées) que dans le cas système de stockage 100' (constamment 26 ou 27 batteries connectées). Cela traduit bien l'optimisation qui est faite en modifiant les états de conduction des contacteurs 160. 25 Les différentes étapes du procédé selon l'invention peuvent être mises en oeuvre par un dispositif de gestion (non représenté) équipant le système de stockage 100 et comportant un processeur. Le procédé selon l'invention est par exemple mis en oeuvre par le biais d'une application. Une application désigne un programme 30 d'ordinateur comprenant des instructions exécutables par le processeur du dispositif de gestion. L'échange d'information entre le système de stockage 100 et l'application se fait par exemple via une communication TCP (Transmission Control Protocol en anglais) à travers une structure de données de type XML. 3059845 18 La figure 8 représente schématiquement un deuxième exemple de contexte de mise en oeuvre du procédé selon l'invention. La figure 8 montre un système de stockage d'énergie électrique 800 comportant toujours N batteries d 120 montés en parallèle. Par contre, chaque batterie 120 est connectée électriquement au bus de 5 distribution 140 par l'intermédiaire d'une résistance variable 170. Les valeurs que peuvent prendre ces résistances 170 forment un intervalle continu. La fonction objectif du problème d'optimisation reste inchangée. Par contre, le problème d'optimisation subit les modifications suivantes : 10 - une variable de décision Rvi correspondant à la valeur de la résistance variable 170 associée à la batterie i est ajoutée ; - la contrainte (C.3) devient : Ui = OCVi + (Rsi + Rvi) * li ; - la variable de décision Xi est supprimée et est remplacé par Xi = 1 dans les contraintes (C.5), (C.6) et (C.7) ; 15 - la contrainte (C.10) est supprimée. L'intérêt d'utiliser des résistances variables à la place des contacteurs est par exemple d'intégrer cette résistance au sein de la batterie afin de permettre aux états de charge des batteries de converger plus rapidement en connectant plus de 20 batteries. Les résistances variables peuvent également permettre de dissiper de la chaleur. Naturellement, l'invention n'est pas limitée aux modes de mise en oeuvre décrits en référence aux figures et des variantes pourraient être envisagées sans 25 sortir du cadre de l'invention. REVENDICATIONS1. Procédé d'équilibrage d'états de charge (SOC) d'une pluralité de dispositifs de stockage d'énergie électrique appartenant à un système de stockage d'énergie électrique (100, 800), les dispositifs de stockage d'énergie (120) étant connectés à un bus de distribution (140) par l'intermédiaire d'éléments de connexion (160, 170), le système de stockage présentant un courant de fonctionnement traversant le bus de distribution, le courant de fonctionnement appartenant à une plage de courant s'étendant entre un courant de charge maximal et un courant de décharge maximal, la plage de courant étant divisée en au moins deux intervalles de courant, un premier intervalle s'étendant entre le courant de charge maximal et un courant nul, et un deuxième intervalle s'étendant entre un courant nul et le courant de décharge maximal, chaque élément de connexion présentant au moins deux états de conduction électrique, le procédé comportant les étapes suivantes : - (a) obtenir l'état de charge des dispositifs de stockage par des moyens de mesures de l'état de charge ; - (b) calculer, à partir des états de charge obtenus au cours de l'étape (a), une cartographie des états de conduction des éléments de connexion associée à chaque intervalle de courant ; - (c) déterminer l'intervalle de courant auquel appartient le courant de fonctionnement du système de stockage ; - (d) appliquer aux éléments de connexion la cartographie associée à l'intervalle de courant déterminé au cours de l'étape (c) ; les étapes (a) et (b) étant répétées périodiquement pour mettre à jour les 25 cartographies des états de conduction des éléments de connexion. 2. Procédé selon la revendication 1 dans lequel l'étape (b) de calcul des cartographies comporte les opérations suivantes : - construire un problème d'optimisation défini par des variables de décision, 30 des contraintes, et une fonction objectif visant à minimiser l'écart entre les états de charge des dispositifs de stockage (120) ; - résoudre le problème d'optimisation. 3059845 20 3. Procédé selon la revendication 2 dans lequel les dispositifs de stockage (120) sont connectés en parallèle au bus de distribution par l'intermédiaire de contacteurs (160), le problème d'optimisation étant non-linéaire. 5 4. Procédé selon la revendication 3 dans lequel le problème d'optimisation est linéarisé. 5. Procédé selon la revendication 2 dans lequel les dispositifs de stockage (120) sont connectés en parallèle au bus de distribution (140) par l'intermédiaire de lo résistances variables (170), le problème d'optimisation étant non-linéaire. 6. Procédé selon l'une quelconque des revendications 2 à 5 comportant une étape préalable de caractérisation d'au moins un dispositif de stockage (120) pour déterminer, en fonction de l'état de charge du dispositif de stockage, les variations de 15 chacun des paramètres caractéristiques du dispositif de stockage suivants : la tension de circuit ouvert, la résistance interne, le courant de charge maximal et le courant de décharge maximum. 7. Procédé selon l'une quelconque des revendications 1 à 6 dans lequel le 20 premier intervalle de courant et/ou le deuxième intervalle de courant sont subdivisés en un nombre prédéterminé d'intervalles de courant. 8. Procédé selon la revendication 7 dans lequel le nombre prédéterminé d'intervalles de courant est paramétré par un utilisateur. 9. Produit programme d'ordinateur comprenant des instructions pour la mise en oeuvre d'un procédé selon l'une quelconque des revendications 1 à 8, lorsque ces instructions sont exécutées par un processeur. 30