FR 3137457 A1 20240105 FR 2206730 A 20220704 Titre : Procédé de détection d'un risque de défaillance par déséquilibre d' un dispositif de stockage d'énergie Procédé de détection d'un risque de défaillance par déséquilibre d'un dispositif de stockage d'énergie (1) comprenant un ensemble d’étages (2) reliés électriquement entre eux en série et constitués de cellules électrochimiques (3) reliées électriquement entre elles en parallèle, caractérisé en ce qu'il comprend : - une étape (E2) de détermination d'une première fonction (f1) caractérisant un fonctionnement correct d'au moins un étage, - une étape (E3) de détermination d'une deuxième fonction (f2) caractérisant le fonctionnement d'un étage présentant la plus faible tension à ses bornes parmi l'ensemble des étages du dispositif de stockage d'énergie (1), - une étape (E4) de calcul d'une différence entre ladite première fonction (f1) et ladite deuxième fonction (f2), puis - une étape (E5) de comparaison de ladite différence à un seuil, ou plusieurs seuils. Figure à publier avec l’abrégé : Figure 1 Procédé de détection d'un risque de défaillance par déséquilibre d'un dispositif de stockage d'énergie comprenant un ensemble d’étages de cellules électrochimiques Domaine technique de l'invention L’invention concerne le domaine de la surveillance des dispositifs de stockage d'énergie comprenant un ensemble d’étages de cellules électrochimiques, les étages étant reliés électriquement en série, chaque étage comprenant une à plusieurs cellules électrochimiques en parallèle, notamment des cellules de type lithium-ion. Plus précisément, l’invention concerne un procédé de détection d'un risque de défaillance par déséquilibre entre les étages en série d'un tel dispositif de stockage d'énergie. L'invention porte aussi sur un équipement de surveillance configuré pour mettre en œuvre un tel procédé de détection. Etat de la technique antérieure Certains dispositifs de stockage d'énergie comprennent un ensemble d’étages de cellules électrochimiques, notamment de type lithium-ion, reliés électriquement en série afin d’obtenir une tension cible souhaitée, chaque étage comprenant une ou plusieurs cellules électrochimiques reliées électriquement en parallèle afin d’obtenir une capacité cible souhaitée. Pour différentes raisons, un déséquilibre d’état de charge couramment dénommé SOC (de l’anglicisme "State-of-Charge") peut apparaître entre les étages en série. Ce déséquilibre est couramment dénommé par l'anglicisme "cell imbalance". Parmi les raisons possibles, on peut citer des problèmes de dispersion d'état de charge lors de l’assemblage des cellules, des problèmes de dispersions d’autodécharge, de capacité, ou encore de résistance entre les cellules, pouvant eux-mêmes être la conséquence de problèmes de dispersion de fabrication des cellules, ou encore des problèmes de dispersion des conditions d’usage en opération entrainant des cinétiques de vieillissement différentes. Une fois observés, ces déséquilibres sont le plus souvent corrigés par un système électronique d’équilibrage. Toutefois, il arrive que le déséquilibre soit tel qu’il ne puisse pas être compensé. Ce déséquilibre peut alors conduire à ce qu'un étage atteigne de manière prématurée sa capacité de charge maximale, respectivement sa décharge maximale, avant les autres étages en série. Si la charge de l'étage est poursuivie après qu'il ait atteint sa capacité de charge maximale, respectivement sa capacité de décharge maximale, il peut résulter une surcharge, respectivement une sous-décharge. Ces états peuvent induire un échauffement indésirable de l’étage en série concerné, provoquer un emballement thermique, voire même un incendie de l'ensemble du dispositif de stockage d'énergie. Pour détecter un risque de déséquilibre d'un étage électrochimique d'un dispositif de stockage d'énergie, la méthode la plus répandue est basée sur l'observation de la tension aux bornes de chaque étage en série lors d'une charge ou d'une décharge complète du système de stockage électrochimique. L’étage de cellules électrochimiques présentant un risque de déséquilibre, ou déjà déséquilibré, possède en règle de générale une tension à ses bornes qui est significativement différente de la tension aux bornes des autres étages et peut ainsi être identifié. Toutefois, cette méthode présente des inconvénients. Les écarts de tension observés sont eux-mêmes fonction des conditions d’usage, à savoir de la température et du courant de charge et de décharge. Ces écarts de tension observés sont également fonction des états de charge et des états de santé (couramment dénommé SOH, de l’anglicisme "State-of-Health") considérés au moment de leur observation. Pour finir, des écarts d'état de charge ne se traduisent pas obligatoirement par des écarts en tension, notamment dans le cas des batteries Lithium Fer Phosphate (LFP), c’est-à-dire à base de phosphate de fer à l’électrode positive, qui présentent une valeur très stable de tension sur une large plage de fonctionnement, dit autrement sur une large plage d'état de charge. Cette méthode apparait ainsi en pratique difficile à calibrer pour éviter des fausses alertes, voire donc insuffisante dans certaines configurations. D'autre part, les écarts de tension entre étages sont parfois trop faibles ou ne deviennent suffisamment importants que très tardivement dans les conditions d’usage, susceptibles alors d’entrainer un problème de sécurité une fois détecté. Ainsi, lorsqu'un risque de déséquilibre est détecté par ce biais, il est généralement nécessaire d'interrompre en urgence l'utilisation du dispositif de stockage d'énergie ce qui perturbe fortement les différents équipements qui lui sont reliés. Les procédés de détection connus de l'état de la technique ne permettent donc pas une gestion simple et sereine de la maintenance des dispositifs de stockage d'énergie. On connait également des indicateurs d'état de santé fournissant un indicateur de l'état de vieillissement d'un dispositif de stockage d'énergie ou d'un étage composant le dispositif de stockage d'énergie. De tels indicateurs sont complexes à calculer et ne permettent pas de détecter un risque de défaillance par déséquilibre d'au moins un étage de cellules électrochimiques composant le dispositif de stockage d'énergie. Parallèlement, avec la généralisation des équipements embarquant une unité de stockage d'énergie, notamment les véhicules automobiles embarquant une batterie lithium-ion, il existe une quantité de plus en plus importante d'unités de stockage d'énergie, dites de seconde vie, qui peuvent être utilisées pour le stockage stationnaire d'énergie, notamment pour stocker de l'énergie électrique produite par une source de production d'énergie intermittente (par exemple de l'énergie solaire ou éolienne) en vue de restituer cette énergie de manière progressive. Ces différentes unités de stockage d'énergie sont rassemblées et connectées ensemble électriquement de manière à former un dispositif de stockage d'énergie de plus grande capacité. Comme les unités de stockage d'énergie qui composent de tels dispositifs de stockage d'énergie peuvent présenter des niveaux d'usure ou d'ancienneté différents, le risque d'observer un déséquilibre entre les unités de stockage d'énergie est particulièrement important. Présentation de l'invention Le but de l’invention est de fournir un procédé de détection d'un risque de défaillance par déséquilibre d'un dispositif de stockage d'énergie comprenant un ensemble d’étages de cellules électrochimiques, les étages étant reliés électriquement en série, le procédé de détection remédiant aux inconvénients ci-dessus et améliorant les procédés de détection connus de l’art antérieur. Plus précisément, un premier objet de l'invention est un procédé permettant de détecter de manière précoce un risque de défaillance par déséquilibre d'un dispositif de stockage d'énergie. Un deuxième objet de l’invention est un procédé de détection d'un risque de défaillance par déséquilibre qui puisse être mis en œuvre lors de charges et/ou de décharges partielles du dispositif de stockage d'énergie. Reprise des revendications L'invention se rapporte à un procédé de détection d'un risque de défaillance par déséquilibre d'un dispositif de stockage d'énergie comprenant un ensemble d’étages reliés électriquement entre eux en série et constitués de cellules électrochimiques reliées électriquement entre elles en parallèle, le procédé de détection comprenant : - une étape de détermination d'une première fonction caractérisant un fonctionnement correct d'au moins un étage, la première fonction définissant une relation entre d'une part une grandeur relative à une quantité de charges circulant dans l'au moins un étage et d'autre part un temps écoulé au cours d’une charge ou d’une décharge de l’au moins un étage, - une étape de détermination d'une deuxième fonction caractérisant le fonctionnement d'un étage présentant la plus faible tension à ses bornes parmi l'ensemble des étages du dispositif de stockage d'énergie, la deuxième fonction définissant une relation entre d'une part ladite grandeur relative à une quantité de charges circulant dans l’étage présentant la plus faible tension à ses bornes parmi l'ensemble des étages du dispositif de stockage d'énergie et d'autre part un temps écoulé au cours d’une charge ou d’une décharge de l’étage présentant la plus faible tension à ses bornes, puis - une étape de calcul d'une différence entre ladite première fonction et ladite deuxième fonction, puis - une étape de comparaison de ladite différence à un seuil. La première fonction peut définir une relation entre d'une part une moyenne sur l'ensemble des étages du dispositif de stockage d'énergie d'une grandeur relative à une quantité de charges circulant dans chaque étage et d'autre part un temps écoulé au cours d’une charge ou d’une décharge de l'ensemble des étages du dispositif de stockage d'énergie. Ladite grandeur relative à une quantité de charges circulant dans un étage peut être une capacité incrémentale de cet étage. L'étape de calcul d'une différence entre ladite première fonction et ladite deuxième fonction peut comprendre une sous-étape de calcul d'une valeur intégrale d'une différence entre la première fonction et la deuxième fonction. L'étape de calcul d'une différence entre ladite première fonction et ladite deuxième fonction peut comprendre: - une sous-étape d'estimation d'un extremum atteint par la première fonction, - une sous-étape d'estimation d'un extremum atteint par la deuxième fonction, puis - une sous-étape de calcul d'une différence entre l'extremum atteint par la première fonction et l'extremum atteint par la deuxième fonction. L'étape de calcul d'une différence entre ladite première fonction et ladite deuxième fonction peut comprendre: - une sous-étape de détection d'un extremum atteint par la première fonction, - une sous étape d'estimation d'une durée de charge ou de décharge au terme de laquelle la première fonction atteint son extremum, - une sous-étape de détection d'un extremum atteint par la deuxième fonction, - une sous-étape d'estimation d'une durée de charge ou de décharge au terme de laquelle la deuxième fonction atteint son extremum, - une sous-étape de calcul d'une différence entre la durée de charge ou de décharge au terme de laquelle la première fonction atteint son extremum et la durée de charge ou de décharge au terme de laquelle la deuxième fonction atteint son extremum. L'étape de calcul d'une différence entre ladite première fonction et ladite deuxième fonction peut comprendre : - une sous-étape de calcul d'un courant de charge moyen du dispositif de stockage d'énergie entre un instant où la première fonction atteint son extremum et un instant où la deuxième fonction atteint son extremum, puis - une sous-étape de calcul d'un déséquilibre d'état de charge par la formule : D_SOC = D_T x I_moy / Q, où : D_SOC désigne le déséquilibre d'état de charge, D_T désigne la différence entre la durée de charge ou de décharge au terme de laquelle la première fonction atteint un extremum et la durée de charge ou de décharge au terme de laquelle la deuxième fonction atteint un extremum correspondant, I_moy désigne le courant de charge moyen, et Q désigne la capacité totale restante de l’étage de cellules considéré, puis - une étape de comparaison dudit déséquilibre d'état de charge à un seuil. L'étape de calcul d'une différence entre ladite première fonction et ladite deuxième fonction peut comprendre : - une sous-étape d'estimation d'un premier extremum atteint par la première fonction, - une sous-étape d'estimation d'au moins un deuxième extremum atteint par la première fonction, - une sous étape d'estimation d'une durée de charge ou de décharge au terme de laquelle la première fonction atteint son premier extremum, - une sous-étape d'estimation d'une durée de charge ou de décharge au terme de laquelle la première fonction atteint son deuxième extremum, - une sous-étape d'estimation d'un premier extremum atteint par la deuxième fonction, - une sous-étape d'estimation d'au moins un deuxième extremum atteint par la deuxième fonction, - une sous-étape d'estimation d'une durée de charge ou de décharge au terme de laquelle la deuxième fonction atteint son premier extremum, - une sous-étape d'estimation d'une durée de charge ou de décharge au terme de laquelle la deuxième fonction atteint son deuxième extremum, puis : - une sous-étape de calcul d'une différence entre le premier extremum de la première fonction et le premier extremum de la deuxième fonction, et/ou - une sous-étape de calcul d'une différence entre le deuxième extremum de la première fonction et le deuxième extremum de la deuxième fonction, et/ou - une sous-étape de calcul d'une différence entre la durée de charge ou de décharge au terme de laquelle la première fonction atteint son premier extremum et la durée de charge ou de décharge au terme de laquelle la deuxième fonction atteint son premier extremum, et/ou - une sous-étape de calcul d'une différence entre la durée de charge ou de décharge au terme de laquelle la première fonction atteint son deuxième extremum et la durée de charge ou de décharge au terme de laquelle la deuxième fonction atteint son deuxième extremum. Ladite première fonction et/ou ladite deuxième fonction peuvent être déterminées : - soit lors d'une phase de charge ou de décharge du dispositif de stockage d'énergie selon un régime lent, notamment un régime inférieur ou égal à C/5, - soit lors d'une phase de charge ou de décharge du dispositif de stockage d'énergie selon un régime rapide, notamment un régime strictement supérieur à C/5, l'étape de détermination de la première fonction et/ou de la deuxième fonction comprenant alors une sous-étape de filtrage de la grandeur relative à une quantité de charges circulant dans un étage. L'étape de comparaison de ladite différence à un seuil peut comprendre : - une sous-étape de comparaison de ladite différence à un premier seuil et à un deuxième seuil, le deuxième seuil étant strictement supérieur au premier seuil, puis - une sous-étape de mémorisation d'un premier témoin d'alerte indiquant un risque modéré, si ladite différence est supérieure ou égale au premier seuil et strictement inférieur au deuxième seuil, et - une sous-étape de mémorisation d'un deuxième témoin d'alerte indiquant un risque élevé, si ladite différence est supérieure ou égale au deuxième seuil. Le premier seuil peut être déterminé en fonction d'une dispersion constatée de ladite différence, et le deuxième seuil peut être déterminé en fonction d'une surcharge admissible par au moins un étage de cellules électrochimiques du dispositif de stockage d'énergie. L'invention se rapporte également à un équipement de surveillance d'un dispositif de stockage d'énergie comprenant un ensemble d’étages électrochimiques reliés électriquement en série, l'équipement de surveillance comprenant des moyens matériels et logiciels configurés pour mettre en œuvre le procédé de détection d'un risque de défaillance par déséquilibre du dispositif de stockage d'énergie tel que défini précédemment. Présentation des figures Ces objets, caractéristiques et avantages de la présente invention seront exposés en détail dans la description suivante d’un mode de réalisation particulier fait à titre non-limitatif en relation avec les figures jointes parmi lesquelles : La est une vue schématique d'un dispositif de stockage d'énergie auquel est relié un équipement de surveillance selon un mode de réalisation de l'invention. La est un synoptique d'un procédé de détection d'un risque de défaillance par déséquilibre d'un dispositif du stockage d'énergie selon un mode de réalisation de l'invention. La est un graphique représentant la capacité incrémentale de différents étages du dispositif de stockage d'énergie en fonction du temps écoulé au cours d'une charge du dispositif du stockage d'énergie. La est un graphique représentant la capacité incrémentale des différents étages du dispositif de stockage d'énergie en fonction de la tension aux bornes de ces étages au cours d'une charge du dispositif du stockage d'énergie. Description détaillée La illustre schématiquement un dispositif de stockage d'énergie 1 comprenant un ensemble d’étages 2 de cellules électrochimiques reliés électriquement entre eux. Les étages 2 sont reliés électriquement en série. Chaque étage 2 peut comprendre une ou plusieurs cellules 3 électrochimiques, également dénommées "accumulateurs" ou "piles rechargeables", reliées électriquement entre elles en série et/ou en parallèle. Chaque cellule 3 comprend une électrode positive, ou cathode, et une électrode négative, ou anode. Les cathodes des différentes cellules 3 sont reliées directement ou indirectement à une borne positive d'un étage 2. De même, les anodes des différentes cellules 3 sont reliées directement ou indirectement à une borne négative d’un étage 2. Les bornes positives et négatives de chaque étage 2 sont reliées respectivement, de manière directe ou indirecte, à une borne positive et négative du dispositif de stockage d'énergie 1. Les différents étages peuvent être assemblés de manière amovible dans le dispositif de stockage d'énergie, de manière à pouvoir être retirés et/ou remplacés. Selon le mode de réalisation illustré sur la , le dispositif de stockage d'énergie 1 comprend quatre étages 2 reliés électriquement en série. Chaque étage 2 comprend six cellules 3 reliées électriquement en parallèle. En variante, le nombre d’étages 2 et/ou de cellules 3 pourrait être différent. Avantageusement, tous les étages 2 comprennent un nombre et un agencement des cellules 3 identiques. Ainsi, ils peuvent comprendre des modes de fonctionnement théoriques sensiblement identiques, notamment des tensions à leurs bornes et une capacité qui sont comparables. Les étages 2 et/ou les cellules 3 composant le dispositif de stockage d'énergie 1 peuvent éventuellement être respectivement des étages et/ou des cellules dits de seconde vie, c’est-à-dire des étages et/ou des cellules issus d'un procédé de re-fabrication après avoir été intégrés au sein d'un premier système. Par exemple, le dispositif de stockage d'énergie 1 peut être composé d'un ensemble de batteries de véhicules automobiles électriques ou hybrides. Ces batteries peuvent avoir servi à stocker de l'énergie pour la propulsion du véhicule au cours d'une première vie, puis avoir été démontées en vue d'une seconde vie lorsque le véhicule était usagé. Le dispositif de stockage d'énergie 1 peut être destiné à stocker de l'énergie électrique produite par une source de production d'énergie intermittente (par exemple de l'énergie solaire ou éolienne). Les cellules 3 composant le dispositif de stockage d'énergie 1 sont de préférence des cellules de type lithium-ion. Dans de telles cellules, des ions lithium peuvent être échangés réversiblement entre l'électrode positive et l'électrode négative. Toutes les cellules 3 d'un même dispositif de stockage d'énergie 1 ont de préférence la même composition chimique. L'électrode négative peut comprendre un matériau à base de graphite (LixC6) ou à base de titanate de lithium (LTO). L'électrode positive peut être à base d’un des matériaux suivants : - Lithium Fer Phosphate (LFP), - Lithium Nickel Manganese Cobalt Oxide (NMC), - Lithium Cobalt Oxide (LCO), - Lithium Nickel Cobalt Aluminium Oxide (NCA), - un mélange de Lithium Cobalt Oxide et de Lithium Nickel Cobalt Aluminium Oxide (Blend LCO-NCA). En variante, les cellules 3 composant le dispositif de stockage d'énergie 1 pourraient être de type sodium-ion. En tout état de cause, les différentes cellules 3 et les étages 2 qui comprennent les cellules 3 sont destinés à fonctionner de manière équilibrée. Le déséquilibre d'un étage 2 peut entraîner des pertes de performances, voire un emballement thermique de cet étage et donc à une défaillance du dispositif de stockage d'énergie 1. Le dispositif de stockage d'énergie 1 comprend également un système de contrôle électronique 4, couramment dénommé BMS (acronyme de l'anglais "Battery Management System"), qui est configuré pour contrôler l'état et/ou le fonctionnement du dispositif de stockage d'énergie 1. Le système de contrôle électronique 4 peut être configuré pour contrôler chaque cellule 3 individuellement ou bien un ensemble de cellules 3 reliées entre elles sous forme d'un étage 2. En particulier, selon le mode de réalisation présenté, le système de contrôle électronique 4 est configuré pour déterminer et/ou mesurer les données suivantes : - une tension moyenne U_moy, égale à la moyenne des tension aux bornes des différents étages 2 ; - une tension minimale U_min, égale à la tension aux bornes de l’étage 2 présentant la tension la plus faible parmi l'ensemble des étages 2; - une tension maximale U_max, égale à la tension aux bornes de l’étage 2 présentant la tension la plus haute parmi l'ensemble des étages 2; - un courant électrique I de charge ou de décharge traversant le dispositif de stockage d'énergie 1. Avantageusement, une large majorité des batteries ou unités de stockage d'énergie produites ou en service à travers le monde comprennent un système de contrôle électronique 4 qui est déjà configuré pour fournir ces données. Il n'est donc pas nécessaire de modifier les systèmes de contrôle électroniques 4 existants pour mettre en œuvre l'invention. En remarque, comme les différents étages 2 sont assemblés en série, le courant électrique traversant le dispositif de stockage d'énergie 1 est égal au courant électrique traversant chacun des étages 2. En outre, le système de contrôle électronique 4 peut également être configuré pour fournir d'autres données parmi lesquelles la tension aux bornes de chaque étage du dispositif de stockage d'énergie 1, l'état de charge du dispositif de stockage d'énergie 1 (couramment dénommé SOC), l'état de santé du dispositif de stockage d'énergie 1 (couramment dénommé SOH), etc... Le système de contrôle électronique 4 est relié via un réseau d'échange de données à un équipement de surveillance 5 selon un mode de réalisation de l'invention. L'équipement de surveillance 5 comprend notamment une mémoire 6, un microprocesseur 7, une interface d'entrée/sortie 8 configurée pour recevoir des données issues du système de contrôle électronique 4 et configurée pour communiquer avec une interface homme-machine 9, par exemple un ordinateur équipé d'un écran. La mémoire 6 est un support d'enregistrement de données comprenant des codes d'instruction qui, lorsqu'ils sont exécutés par le microprocesseur 7, conduisent celui-ci à mettre en œuvre un procédé de détection d'un risque de défaillance par déséquilibre du dispositif de stockage d'énergie 1, selon un mode de réalisation de l'invention. L'équipement de surveillance 5 peut être relié au système de contrôle électronique 4 via un réseau d'échange de données tel qu'Internet. En variante, l'équipement de surveillance 5 peut être intégré à un boîtier relié au système de contrôle électronique 4 par une liaison filaire directe, voire même être intégré au système de contrôle électronique 4. On décrit à présent un premier mode de réalisation d'un procédé de détection d'un risque de défaillance par déséquilibre du dispositif de stockage d'énergie 1 selon l'invention. Le procédé repose sur des données calculées ou mesurées, par le système de contrôle électronique 4, lors d'une phase de charge ou de décharge de l'unité de stockage d'énergie 1. Avantageusement, le procédé ne requiert pas une charge ou une décharge complète du dispositif de stockage d'énergie 1. Au contraire, seule une charge ou une décharge partielle suffisent pour la mise en œuvre du procédé. Par exemple, le procédé peut être mis en œuvre au cours d'une charge ou d'une décharge dans laquelle l'état de charge du dispositif de stockage d'énergie 1 varie entre 25% et 75% de sa capacité de charge totale. Le procédé de détermination peut être décomposé en cinq étapes E1, E2, E3, E4, E5 représentées schématiquement sur la . Dans une première étape E1, le système de contrôle électronique 4 transmet à l'équipement de surveillance 5 les valeurs des grandeurs suivantes : - la tension U_min de l’étage 2 présentant la plus basse tension, - la tension U_max de l’étage 2 présentant la plus haute tension, - la tension moyenne U_moy aux bornes des différents étages, - le courant électrique I circulant dans le dispositif de stockage d'énergie. Ces valeurs peuvent être par exemple transmises sous forme de séries temporelles, de manière périodique et/ou à l'issue de chaque phase de charge ou de décharge du dispositif de stockage d'énergie 1. Dans une deuxième étape E2, on détermine une première fonction f1, dite fonction de référence, caractérisant un fonctionnement correct d'au moins un étage 2. Par "fonctionnement correct", on comprend un fonctionnement normal ou nominal d'au moins un étage 2, c’est-à-dire le fonctionnement d'un étage non défaillant. Selon le premier mode de réalisation, la première fonction f1 est égale à une fonction moyenne calculée sur la base de l'ensemble des étages du dispositif de stockage d'énergie 1. Ce premier mode de réalisation repose donc sur l'hypothèse selon laquelle la moyenne de tous les étages est représentative d'un fonctionnement correct. La définition de la première fonction sur la base d'une moyenne de tous les étages du dispositif de stockage permet de rendre le procédé de détection plus robuste, et notamment de conserver une détection efficace même lorsqu'un des étages présente une tension anormalement haute à ses bornes. On peut éventuellement convenir que ce mode de réalisation ne peut être mis en œuvre que pour un dispositif de stockage d'énergie comprenant un nombre d’étages suffisant, afin que la moyenne calculée sur l'ensemble des étages traduise bien, selon les lois de la statistique, un fonctionnement correct. Alternativement, et comme nous le verrons par la suite, d'autres méthodes permettant de déterminer la fonction de référence peuvent être proposées. D'une manière générale, la première fonction f1 est une fonction mathématique, représentable sur un graphique tel que le graphique de la , et qui peut être définie par un ensemble de points. La première fonction définit une relation entre d'une part une grandeur relative à une quantité de charges circulant dans au moins un étage (représenté en ordonnées sur la ), et d'autre part un temps écoulé au cours d’une charge ou d’une décharge de l’au moins un étage (représenté en abscisse sur la , et exprimé par exemple en heures). En l'espèce, les fonctions représentées sur la sont représentatives d'une charge du dispositif de stockage d'énergie : l'état de charge augmente lorsqu'on progresse sur l'axe des abscisses. Ces fonctions peuvent donc être déterminées au cours d'une phase de charge du dispositif de stockage d'énergie. En variante, ces fonctions peuvent aussi être calculées lors de décharges du dispositif de stockage d'énergie en inversant l'orientation de l'axe des abscisses. Dans l'hypothèse où ladite grandeur relative à une quantité de charges circulant dans au moins un étage change de signe dans lors des phases de décharge comparativement aux phases de charge, on utilise avantageusement une valeur absolue de cette grandeur. Dans l'hypothèse où la première fonction f1 est égale à une fonction moyenne calculée sur la base de l'ensemble des étages, ledit au moins un étage correspond à l'ensemble des étages du dispositif de stockage d'énergie 1. Selon un mode de réalisation préféré, la grandeur relative à une quantité de charges circulant dans au moins un étage est égale à une capacité incrémentale (dQ/dU, exprimé par exemple en Ampères-heure par volt) de l'au moins un étage. La capacité incrémentale d'un étage est définie par un rapport d'un différentiel de quantité de charges dQ de cet étage sur un différentiel de tension dU aux bornes de cet étage. La première fonction f1 peut être déterminée de la manière suivante : tout d'abord on calcule une première fonction intermédiaire définissant une relation entre le courant électrique I circulant dans le dispositif de stockage d'énergie et le temps écoulé au cours d'une période de charge ou d'une décharge du dispositif de stockage d'énergie. Ensuite, on calcule une deuxième fonction intermédiaire définissant une relation entre une quantité de charges Q circulant dans chaque étage et le temps écoulé en intégrant la première fonction intermédiaire sur la période de charge ou de décharge considérée. Cette deuxième fonction intermédiaire est combinée avec une troisième fonction intermédiaire établissant une relation entre la tension moyenne U_moy et le temps écoulé. On peut ainsi calculer une quatrième fonction intermédiaire définissant une relation entre la quantité de charges Q circulant dans chaque étage et la tension moyenne U_moy. Ensuite on calcule une cinquième fonction intermédiaire en dérivant la quatrième fonction intermédiaire relativement à la tension moyenne U_moy. La cinquième fonction intermédiaire est donc une fonction du type dQ/dU_moy = f(U_moy). Enfin, cette cinquième fonction intermédiaire est combinée avec la troisième fonction intermédiaire établissant une relation entre la tension moyenne U_moy et le temps écoulé t de manière à obtenir la première fonction f1. La première fonction est donc une fonction du type dQ/dU_moy = f(t). Alternativement, pour déterminer la première fonction f1, on peut déterminer pour chaque étage la fonction définissant la relation entre la capacité incrémentale de cet étage et le temps écoulé. Ensuite, on peut effectuer une moyenne arithmétique des fonctions déterminées pour chaque étage. Cette méthode permet une détection plus précise mais requiert davantage de ressources de calcul car les calculs sont répétés pour chaque étage du dispositif de stockage d'énergie. De plus, cette méthode requiert que le système de contrôle électronique 4 fournisse la tension aux bornes de chaque étage du dispositif de stockage d'énergie. Selon une variante de réalisation de la deuxième étape E2, l'au moins un étage dont le fonctionnement est correct pourrait être défini comme l’étage 2 dont la tension à ses bornes est la plus proche de la tension moyenne de l'ensemble des étages 2 du dispositif de stockage d'énergie 1. Selon d'autres variantes de réalisation de la deuxième étape E2, la première fonction pourrait être définie différemment, par exemple au moyen d'une fonction théorique ou encore en identifiant par tout moyen un ou plusieurs étages du dispositif de stockage d'énergie 1 qui fonctionne correctement et en déterminant la relation entre la capacité incrémentale circulant dans ce ou ces étages et le temps écoulé au cours d'une charge ou d'une décharge. Finalement, à l'issue de la deuxième étape E2 on obtient une première fonction f1, représentative d'un fonctionnement normal d'un ou plusieurs étages. Cette première fonction peut être déterminée selon plusieurs méthodes différentes mais qui possèdent le point commun de définir une relation entre d'une part une grandeur relative à une quantité de charges circulant dans un étage et d'autre part un temps écoulé au cours d'une charge ou d'une décharge du dispositif de stockage d'énergie. Cette première fonction est donc une fonction de référence et sert de base de comparaison pour déterminer si un étage particulier présente un risque de défaillance par déséquilibre. Dans une troisième étape E3, on détermine une deuxième fonction f2 destinée à être comparée avec la première fonction précédemment définie. De même que pour la première fonction f1, la deuxième fonction f2 est une fonction mathématique, représentable sur un graphique tel que le graphique de la , et qui peut être définie par un ensemble de points. La troisième étape E3 peut être exécutée avant ou après la deuxième étape E2 ou encore en parallèle de la deuxième étape E2. La deuxième fonction f2 définit une relation entre d'une part la grandeur relative à une quantité de charges circulant dans l’étage présentant la plus faible tension à ses bornes parmi l'ensemble des étages du dispositif de stockage d'énergie et d'autre part un temps écoulé au cours d’une charge ou d’une décharge de cet étage. En particulier, selon le mode de réalisation préféré, la grandeur relative à la quantité de charges circulant dans l’étage présentant la plus faible tension à ses bornes est égale à la capacité incrémentale de cet étage. La méthode de détermination de la deuxième fonction f2 peut être analogue à la méthode de détermination de la première fonction f1. La deuxième fonction f2 peut être déterminée de la manière suivante : tout d'abord on calcule une première fonction intermédiaire définissant une relation entre le courant électrique I circulant dans le dispositif de stockage d'énergie et le temps écoulé au cours d'une période de charge ou d'une décharge du dispositif de stockage d'énergie. Ensuite, on calcule une deuxième fonction intermédiaire définissant une relation entre une quantité de charges Q circulant dans l’étage présentant la plus faible tension à ses bornes et le temps écoulé en intégrant la première fonction intermédiaire sur la période de charge ou de décharge. Cette deuxième fonction intermédiaire est combinée avec une troisième fonction intermédiaire établissant une relation entre la tension minimale U_min et le temps écoulé. On peut ainsi calculer une quatrième fonction intermédiaire définissant une relation entre la quantité de charges Q circulant dans l’étage présentant la plus faible tension à ses bornes et la tension minimale U_min. Ensuite on calcule une cinquième fonction intermédiaire en dérivant la quatrième fonction intermédiaire relativement à la tension minimale U_min. La cinquième fonction intermédiaire est donc une fonction du type dQ/dU_min = f(U_min). Enfin, cette cinquième fonction intermédiaire est combinée avec la troisième fonction intermédiaire établissant une relation entre la tension moyenne U_min et le temps écoulé de manière à obtenir la deuxième fonction f2. La deuxième fonction est donc une fonction du type dQ/dU_min = f(t). D'une manière générale, le graphique sur lequel la deuxième fonction est représentable est identique au graphique sur lequel la première fonction est représentable. Autrement dit, la forme de la première fonction est la même que la forme de la deuxième fonction de manière à permettre une comparaison de ces deux fonctions. En remarque, la quantité de charges circulant dans chaque étage est calculée de manière plus précise lors d'une phase de charge ou de décharge du dispositif de stockage d'énergie selon un régime lent, notamment un régime inférieur ou égal à C/5, c’est-à-dire avec un courant de charge permettant de recharger complètement le dispositif de stockage d'énergie en au moins cinq heures. Alternativement, la quantité de charges circulant dans chaque étage peut aussi être calculée lors d'une phase de charge ou de décharge du dispositif de stockage d'énergie selon un régime plus rapide, notamment un régime strictement supérieur à C/5. Dans ce cas, les étapes E2 et E3 de détermination de la première fonction et/ou de la deuxième fonction comprenant avantageusement une sous-étape de filtrage de la grandeur relative à une quantité de charges circulant dans un étage. La première fonction f1 et la deuxième fonction f2 sont définies sur une période de charge ou de décharge donnée qui peut correspondre à une charge ou d'une décharge partielle du dispositif de stockage d'énergie. Ainsi, la période sur laquelle la première fonction f1 et la deuxième fonction f2 sont définies peut être restreinte comparativement à la durée totale nécessaire pour charger ou décharger complètement le dispositif de stockage d'énergie, avec un même régime de charge, à partir d'un état de charge de 0% dans le cas d'une charge ou respectivement à partir d'un état de charge de 100% dans le cas d'une décharge. Par exemple, la période sur laquelle la première fonction f1 et la deuxième fonction f2 sont définies peut être inférieure ou égale à 75%, voire inférieure ou égale à 50%, voire même inférieure ou égale à 25% de la durée totale. De préférence, la période sur laquelle la première fonction f1 et la deuxième fonction f2 sont définies est suffisante pour identifier au moins un extremum de la première fonction et au moins un extremum de la deuxième fonction, voire au moins deux extremums de la première fonction et au moins deux extremums de la deuxième fonction, voire même trois extremums de la première fonction et trois extremums de la deuxième fonction. Comme cela est bien visible sur la , lorsque la période sur laquelle la première fonction f1 et la deuxième fonction f2 sont définies est suffisamment importante, la première fonction f1 et la deuxième fonction f2 possède une allure spécifique : chacune de ces deux fonctions atteint successivement deux valeurs maximales, référencées VM11 et VM12 pour la première fonction f1 et VM21 et VM22 pour la deuxième fonction f2. Les valeurs maximales VM11, VM12, VM21 et VM22 sont atteintes respectivement au terme d'une durée T11, T12, T21 et T22. Lorsque la première fonction f1 et la deuxième fonction f2 sont établies au cours d'une charge du dispositif de stockage d'énergie, la première valeur maximale VM11 atteinte par la première fonction est généralement inférieure ou égale à la deuxième valeur maximale VM12 atteinte par la première fonction. De même, la première valeur maximale VM21 de la deuxième fonction est généralement inférieure ou égale à la deuxième valeur maximale VM22 de la deuxième fonction. Par ailleurs, on observe un décalage entre la première fonction f1 et la deuxième fonction f2. En particulier, les valeurs maximales VM11 et VM12 de la première fonction sont strictement supérieures respectivement aux valeurs maximales VM21 et VM22 de la deuxième fonction. De plus, les durées T11 et T12 sont strictement inférieures respectivement aux durées T21 et T22. Entre ses deux valeurs maximales VM11 et VM12, la première fonction f1 atteint une valeur minimale VM13 au terme d'une durée T13. De même, entre ses deux valeurs maximales VM21 et VM22, la deuxième fonction f2 atteint une valeur minimale VM23 au terme d'une durée T23. On note que la valeur minimale VM13 est strictement inférieure à la valeur maximale VM23 et que la durée T13 est strictement inférieure à la durée T23. Les deux valeurs maximales VM11 et VM12 et la valeur minimale VM13 constituent trois extremums de la première fonction. De même, les deux valeurs maximales VM21 et VM22 et la valeur minimale VM23 constituent trois extremums de la deuxième fonction. De plus, on observe aussi qu'il existe des instants où la première fonction est strictement supérieure à la deuxième fonction et d'autres instants où la première fonction est strictement inférieure à la deuxième fonction. Dans une quatrième étape E4, on calcule une différence entre ladite première fonction f1 et ladite deuxième fonction f2. Il existe plusieurs manières de quantifier une telle différence. Selon un premier mode de réalisation, la quatrième étape E4 comprend une sous-étape de calcul d'une valeur intégrale d'une différence entre la première fonction f1 et la deuxième fonction f2. Ce calcul intégral peut notamment être réalisé sur l'ensemble de la période sur laquelle la première fonction et la deuxième fonction sont définies. Ce calcul revient donc sensiblement à calculer l'aire définie entre la première fonction et la deuxième fonction. Avantageusement, le calcul intégral peut être basé sur une valeur absolue de la différence entre la première fonction f1 et la deuxième fonction f2. Ainsi, le calcul intégral des périodes où la première fonction est strictement supérieure à la deuxième fonction s'additionne au calcul intégral des périodes où la première fonction est strictement inférieure à la deuxième fonction. On parvient ainsi à bien mettre en évidence les différences entre la première fonction et la deuxième fonction, ce qui permet d'améliorer la sensibilité de la détection. Un avantage à déterminer ladite différence sur la base d'un calcul intégral est que cette méthode peut être mise en œuvre sur n'importe quelle période de charge ou de décharge, y compris une période de charge ou de décharge partielle dans laquelle les fonctions f1 et/ou f2 n'atteignent pas tous leurs extremums. Ainsi, le procédé de détection permet de détecter un déséquilibre même lorsque le dispositif de stockage d'énergie subit des cycles de charge et de décharge incomplets. Un autre avantage à déterminer ladite différence sur la base d'un calcul intégral est que cette méthode permet de détecter une différence même lorsque les valeurs maximales VM11 et VM12 de la première fonction sont sensiblement égales aux valeurs maximales VM21 et VM22 de la deuxième fonction (on a alors uniquement un décalage temporel entre les deux fonctions f1 et f2). De même, cette méthode permet de détecter une différence même lorsque les durées T11 et T12 sont sensiblement égales aux durées T21 et T22 de la deuxième fonction (on a alors uniquement un décalage en amplitude entre les deux fonctions f1 et f2). Selon un autre mode de réalisation de la quatrième étape E4, la différence entre la première fonction f1 et la deuxième fonction f2 peut être calculée par une différence d'amplitude entre ces deux fonctions. Dans ce cas, la quatrième étape comprend une sous-étape d'estimation d'au moins un extremum (notamment les valeurs VM11, VM12 ou VM13) atteinte par la première fonction, et une sous-étape d'estimation d'au moins un extremum atteint par la deuxième fonction (respectivement les valeurs VM21, VM22 ou VM23). Ensuite, la quatrième étape E4 comprend une sous-étape de calcul d'une différence entre l'extremum atteint par la première fonction et l'extremum atteint par la deuxième fonction. Cette différence est donc égale à VM11 – VM21, ou à VM12 – VM22, ou encore à VM13 – VM23. La différence peut aussi être égale à un résultat calculé en fonction des trois différences à VM11 – VM21, VM12 – VM22, et VM13 – VM23, ou encore égale un résultat calculé en fonction de deux différences parmi ces trois différences. En raison de la similitude d'allure de la première fonction et de la deuxième fonction, il convient de comparer respectivement les extremums VM11, VM12 et VM13 avec les extremums VM21, VM22 et VM23. Le calcul d'une différence entre des extremums non correspondants conduirait à un résultat aberrant qui peut être filtré pour ne pas faire reposer le procédé sur une telle comparaison. La mise en œuvre de la quatrième étape E4 par comparaison des valeurs maximales ou minimales atteintes par la première fonction et par la deuxième fonction présente l'avantage d'être économe en calcul et simple à mettre en œuvre. Selon un autre mode de réalisation de la quatrième étape E4, la différence entre la première fonction f1 et la deuxième fonction f2 peut être calculée par une différence temporelle entre les deux fonctions f1 et f2. Dans ce cas, la quatrième étape comprend une sous-étape de détection d'un extremum VM11, VM12, VM13 atteint par la première fonction f1 suivi d'une sous étape d'estimation de la durée T11, T12, T13 de charge ou de décharge au terme de laquelle la première fonction atteint son extremum. De même pour la deuxième fonction : la quatrième étape comprend une sous-étape de détection d'un extremum VM21, VM22, VM23 atteint par la deuxième fonction, puis une sous étape d'estimation de la durée T21 T22, T23 de charge ou de décharge au terme de laquelle la deuxième fonction atteint son extremum. Ensuite, on calcule au moins une des différences T11 –T21, et/ou T12 – T22, et/ou T13 – T23. La différence peut ainsi être égale à T11 – T21, ou à T12 – T22, ou encore à T13 – T23. La différence peut aussi être égale à un résultat calculé en fonction des trois différences à T11 – T21, T12 – T22, et T13 – T23, ou encore égale un résultat calculé en fonction de deux différences parmi ces trois différences. Cette variante de réalisation présente l'avantage de détecter un déséquilibre entre les étages qui se traduit principalement par un décalage temporel entre les deux fonctions. On parvient ainsi à détecter un déséquilibre dans l'inertie de charge et de décharge des différents étages. Selon encore une autre variante de réalisation de la quatrième étape, ladite différence entre la première fonction et la deuxième fonction peut être égale à tout résultat calculé en fonction de tout ou partie des six différences VM11 – VM21, VM12 – VM22, VM13 – VM23, T11 – T21, T12 – T22, et T13 – T23. Dans une cinquième étape E5, on compare la différence calculée lors de la quatrième étape avec un seuil. Ensuite, si la différence est strictement supérieure audit seuil, un témoin peut être mémorisée dans la mémoire 6 de l'équipement de surveillance 5. Ce témoin peut être lu par l'interface homme-machine 9. Puis, l'interface homme-machine 9 peut générer un message d'alerte indiquant qu'un étage du dispositif de stockage d'énergie présente un risque de défaillance par déséquilibre. De manière avantageuse, la comparaison de fonctions définissant une relation entre une grandeur relative à une quantité de charges circulant dans un étage et un temps écoulé au cours d'une charge ou d'une décharge du dispositif de stockage d'énergie permet de détecter de manière très anticipée une dérive annonciatrice d’un risque d’emballement thermique. Il a été ainsi observé des dispositifs de stockage d'énergie 1 dont la simple observation de la tension aux bornes des différents étages ne permettait pas d'identifier une quelconque anomalie plusieurs mois avant qu'une défaillance se produise. En revanche, la mise en œuvre du procédé selon l'invention sur ce dispositif de stockage d'énergie permet d'identifier un risque de défaillance par déséquilibre plusieurs mois avant que celui-ci ne se produise. De plus, le procédé de détection permet en général d'identifier l’étage du dispositif de stockage d’énergie responsable de cette anomalie. L’étage en question peut alors facilement être retiré ou remplacé lors d'une opération de maintenance. La illustre de manière comparative deux fonction f1' et f2' caractérisant respectivement un fonctionnement correct d'au moins un étage et le fonctionnement d'un étage présentant la plus faible tension à ses bornes parmi l'ensemble des étages du dispositif de stockage d'énergie. Les deux fonctions f1' et f2' établissent une relation entre une capacité incrémentale dQ/dU (en ordonnée) et une tension U aux bornes de l'étage concerné. Les fonctions f1' et f2' sont construites sur la base du même dispositif de stockage d'énergie que précédemment décrit et lors de la même phase de charge que celle ayant permis de calculer les fonctions f1 et f2 représentées sur la . On observe que la différence entre la première fonction f1' et f2' sur la est nettement moins perceptible que la différence entre les fonctions f1 et f2 sur la . Notamment les fonctions f1 et f2 sont décalées l'une de l'autre à la fois en abscisse et en ordonnée alors que les fonctions f1' et f2' sont décalées uniquement en ordonnées. Comparativement à une méthode basée sur des fonctions établissant une relation entre la capacité incrémentale et une tension, la méthode qui vient d'être décrite, basée sur des fonctions établissant une relation entre la capacité incrémentale et un temps écoulé de charge ou de décharge, permet de mieux mettre en évidence une différence de comportement entre les différents étages et donc de détecter de manière plus précoce un risque de défaillance. Selon un perfectionnement de l'invention, la cinquième étape E5 peut comprendre : - une étape E51 de comparaison de ladite différence à un premier seuil et à un deuxième seuil, le deuxième seuil étant strictement supérieur au premier seuil, puis - une étape E52 de mémorisation d'un premier témoin d'alerte indiquant un risque modéré, si ladite différence est supérieure ou égale au premier seuil et strictement inférieure au deuxième seuil, et - une étape E53 de mémorisation d'un deuxième témoin d'alerte indiquant un risque élevé, si ladite différence est supérieure ou égale au deuxième seuil. Le premier témoin et le deuxième témoin sont destinés à être enregistrés dans la mémoire 6 de l'équipement de surveillance 5. Ces témoins peuvent ensuite être consultés par l'interface homme-machine 9 afin de produire un message d'alerte adapté à la situation. Avantageusement, le premier seuil est déterminé en fonction d'une dispersion normale de ladite grandeur relative à une quantité de charges circulant dans un étage. On peut par exemple déterminer ce premier seuil de manière expérimentale en observant des étages fonctionnant correctement dans un dispositif de stockage d'énergie. Le premier seuil peut être ainsi défini comme égal ou légèrement supérieur à la plus grande différence (telle que calculée lors de l'étape E4) observée sur une période suffisamment longue, avec un dispositif de stockage d'énergie dont tous les étages fonctionnent correctement. Le premier seuil peut aussi correspondre à un seuil au-delà duquel le système de contrôle électronique 4 ne peut plus parvenir à compenser des déséquilibres entre les différents étages du dispositif de stockage d'énergie. Le deuxième seuil peut être déterminé en fonction d'une surcharge admissible par au moins un étage de cellules électrochimiques du dispositif de stockage d'énergie. La surcharge admissible désigne le pourcentage de charges qu'est capable de supporter un étage avant une dégradation irréversible. Autrement dit, si cette différence dépasse la surcharge admissible, alors un emballement thermique se produira de manière certaine. On peut donc avantageusement définir le deuxième seuil comme une fraction de la surcharge admissible. Selon un autre perfectionnement de l'invention, la quatrième étape E4 peut comprendre une sous-étape de calcul d'un courant de charge moyen I_moy du dispositif de stockage d'énergie entre l'instant où la première fonction f1 atteint un des extremums VM11, VM12, VM13 et l'instant où la deuxième fonction f2 atteint un des extremums VM21, VM22, VM23 correspondant. De préférence, le courant de charge I (dont la valeur est fournie par le système de contrôle électronique 4) est sensiblement constant entre ces deux instants et le courant de charge moyen I_moy est égal à cette valeur. En variante, le courant de charge I peut subir certaines variations entre ces deux instants et dans ce cas une moyenne temporelle peut être calculée. Ensuite, la quatrième étape E4 peut comprendre une sous-étape de calcul d'une variation d'état de charge D_SOC au moyen d'une multiplication de la différence entre la durée T11, T12, T13 de charge ou de décharge au terme de laquelle la première fonction atteint un de ses extremums et la durée T21 T22, T23 de charge ou de décharge au terme de laquelle la deuxième fonction atteint l'extremum VM21, VM22, VM23 correspondant avec le courant de charge moyen I_moy. Autrement dit, le déséquilibre d'état de charge D_SOC peut être calculé par la formule suivante : D_SOC = D_T x I_moy / Q où Q désigne la capacité totale restante de l’étage de cellules concerné, c’est-à-dire la capacité totale de l'étage considéré lors de la mise en œuvre du procédé, et où D_T est égal au résultat de la soustraction T11-T21 ou T12-T22 ou T13-T23 . Ainsi, le déséquilibre d'état de charge D_SOC est égal à un rapport entre deux quantités de charges et peut être exprimé en pourcentage. Ensuite, lors d’une cinquième étape E5, on peut comparer le déséquilibre d’état de charge D_SOC à un seuil, ou à plusieurs seuils pour quantifier le niveau de criticité de l’alerte. En particulier, le déséquilibre d’état de charge D_SOC peut être comparé à un premier seuil et à un deuxième seuil comme expliqué précédemment. Le premier seuil peut être par exemple compris entre 5% et 20%. Le second seuil d’alerte, strictement supérieur au premier seuil d'alerte, peut être par exemple compris entre 20% et 40%. Finalement, grâce à l'invention, on dispose d'un procédé de détection précoce d'un risque de défaillance par déséquilibre d'un étage d'un dispositif de stockage d'énergie qui peut être mis en œuvre lors de charges et/ou de décharges partielles. Comparativement aux procédés connus, ce procédé permet de détecter un déséquilibre de manière précoce, ce qui permet une meilleure maintenance du dispositif de stockage d'énergie. En particulier, le procédé permet de détecter des déséquilibres jusqu'alors difficiles à détecter, notamment des déséquilibres causés par un écart d'état de charge (SOC) sur des chimies ne présentant pas de relation notable entre tension et état de charge (SOC) tel qu'une chimie Li-ion LFP, ou encore causé par un écart d'état de santé (SOH) des différents étages du dispositif de stockage d'énergie. Procédé de détection d'un risque de défaillance par déséquilibre d'un dispositif de stockage d'énergie (1) comprenant un ensemble d’étages (2) reliés électriquement entre eux en série et constitués de cellules électrochimiques (3) reliées électriquement entre elles en parallèle, caractérisé en ce qu'il comprend : - une étape (E2) de détermination d'une première fonction (f1) caractérisant un fonctionnement correct d'au moins un étage, la première fonction définissant une relation entre d'une part une grandeur relative à une quantité de charges circulant dans l'au moins un étage et d'autre part un temps écoulé au cours d’une charge ou d’une décharge de l’au moins un étage, - une étape (E3) de détermination d'une deuxième fonction (f2) caractérisant le fonctionnement d'un étage présentant la plus faible tension à ses bornes parmi l'ensemble des étages du dispositif de stockage d'énergie (1), la deuxième fonction définissant une relation entre d'une part ladite grandeur relative à une quantité de charges circulant dans l’étage présentant la plus faible tension à ses bornes parmi l'ensemble des étages du dispositif de stockage d'énergie et d'autre part un temps écoulé au cours d’une charge ou d’une décharge de l’étage présentant la plus faible tension à ses bornes, puis - une étape (E4) de calcul d'une différence entre ladite première fonction (f1) et ladite deuxième fonction (f2), puis - une étape (E5) de comparaison de ladite différence à un seuil. Procédé de détection selon la revendication précédente, caractérisé en ce que la première fonction (f1) définit une relation entre d'une part une moyenne sur l'ensemble des étages du dispositif de stockage d'énergie d'une grandeur relative à une quantité de charges circulant dans chaque étage et d'autre part un temps écoulé au cours d’une charge ou d’une décharge de l'ensemble des étages du dispositif de stockage d'énergie. Procédé de détection selon l'une des revendications précédentes, caractérisé en ce que ladite grandeur relative à une quantité de charges circulant dans un étage est une capacité incrémentale de cet étage. Procédé de détection selon l'une des revendications précédentes, caractérisé en ce que l'étape (E4) de calcul d'une différence entre ladite première fonction et ladite deuxième fonction comprend une sous-étape de calcul d'une valeur intégrale d'une différence entre la première fonction (f1) et la deuxième fonction (f2). Procédé de détection selon l'une des revendications précédentes, caractérisé en ce que l'étape (E4) de calcul d'une différence entre ladite première fonction (f1) et ladite deuxième fonction (f2) comprend : - une sous-étape d'estimation d'un extremum (VM11, VM12, VM13) atteint par la première fonction, - une sous-étape d'estimation d'un extremum (VM21, VM22, VM23) atteint par la deuxième fonction, puis - une sous-étape de calcul d'une différence entre l'extremum (VM11, VM12, VM13) atteint par la première fonction et l'extremum (VM21, VM22, VM23) atteint par la deuxième fonction. Procédé de détection selon l'une des revendications précédentes, caractérisé en ce que l'étape (E4) de calcul d'une différence entre ladite première fonction (f1) et ladite deuxième fonction (f2) comprend : - une sous-étape de détection d'un extremum (VM11, VM12, VM13) atteint par la première fonction (f1), - une sous étape d'estimation d'une durée (T11, T12, T13) de charge ou de décharge au terme de laquelle la première fonction atteint son extremum (VM11, VM12, VM13), - une sous-étape de détection d'un extremum (VM21, VM22, VM23) atteint par la deuxième fonction (f2), - une sous-étape d'estimation d'une durée (T21 T22, T23) de charge ou de décharge au terme de laquelle la deuxième fonction atteint son extremum (VM21, VM22, VM23), - une sous-étape de calcul d'une différence entre la durée (T11, T12, T13) de charge ou de décharge au terme de laquelle la première fonction atteint son extremum et la durée (T21 T22, T23) de charge ou de décharge au terme de laquelle la deuxième fonction atteint son extremum (VM21, VM22, VM23). Procédé de détection selon la revendication précédente, caractérisé en ce que l'étape (E4) de calcul d'une différence entre ladite première fonction (f1) et ladite deuxième fonction (f2) comprend : - une sous-étape de calcul d'un courant de charge moyen du dispositif de stockage d'énergie entre un instant où la première fonction atteint son extremum (VM11, VM12, VM13) et un instant où la deuxième fonction atteint son extremum (VM21, VM22, VM23), puis - une sous-étape de calcul d'un déséquilibre d'état de charge par la formule : D_SOC = D_T x I_moy / Q, où : D_SOC désigne le déséquilibre d'état de charge, D_T désigne la différence entre la durée (T11, T12, T13) de charge ou de décharge au terme de laquelle la première fonction atteint un extremum et la durée (T21 T22, T23) de charge ou de décharge au terme de laquelle la deuxième fonction atteint un extremum (VM21, VM22, VM23) correspondant, I_moy désigne le courant de charge moyen, et Q désigne la capacité totale restante de l’étage de cellules considéré, puis - une étape (E5) de comparaison dudit déséquilibre d'état de charge à un seuil. Procédé de détection selon l'une des revendications précédentes, caractérisé en ce que l'étape (E4) de calcul d'une différence entre ladite première fonction (f1) et ladite deuxième fonction (f2) comprend : - une sous-étape d'estimation d'un premier extremum (VM11) atteint par la première fonction, - une sous-étape d'estimation d'au moins un deuxième extremum (VM12, VM13) atteint par la première fonction, - une sous étape d'estimation d'une durée (T11) de charge ou de décharge au terme de laquelle la première fonction atteint son premier extremum (VM11), - une sous-étape d'estimation d'une durée (T12, T13) de charge ou de décharge au terme de laquelle la première fonction atteint son deuxième extremum (VM12, VM13), - une sous-étape d'estimation d'un premier extremum (VM21) atteint par la deuxième fonction, - une sous-étape d'estimation d'au moins un deuxième extremum (VM22, VM23) atteint par la deuxième fonction, - une sous-étape d'estimation d'une durée (T21) de charge ou de décharge au terme de laquelle la deuxième fonction atteint son premier extremum (VM21), - une sous-étape d'estimation d'une durée (T22, T23) de charge ou de décharge au terme de laquelle la deuxième fonction atteint son deuxième extremum (VM22, VM23), puis : - une sous-étape de calcul d'une différence entre le premier extremum (VM11) de la première fonction et le premier extremum (VM21) de la deuxième fonction, et/ou - une sous-étape de calcul d'une différence entre le deuxième extremum (VM12, VM13) de la première fonction et le deuxième extremum (VM22, VM23) de la deuxième fonction, et/ou - une sous-étape de calcul d'une différence entre la durée (T11) de charge ou de décharge au terme de laquelle la première fonction atteint son premier extremum (VM11) et la durée (T21) de charge ou de décharge au terme de laquelle la deuxième fonction atteint son premier extremum (VM21), et/ou - une sous-étape de calcul d'une différence entre la durée (T12, T13) de charge ou de décharge au terme de laquelle la première fonction atteint son deuxième extremum (VM12, VM13) et la durée (T22, T23) de charge ou de décharge au terme de laquelle la deuxième fonction atteint son deuxième extremum (VM22, VM23). Procédé de détection selon l’une des revendications précédentes, caractérisé en ce que ladite première fonction (f1) et/ou ladite deuxième fonction (f2) sont déterminées : - soit lors d'une phase de charge ou de décharge du dispositif de stockage d'énergie selon un régime lent, notamment avec un courant de charge permettant de recharger complètement le dispositif de stockage d'énergie en au moins cinq heures, - soit lors d'une phase de charge ou de décharge du dispositif de stockage d'énergie selon un régime rapide, notamment avec un courant de charge permettant de recharger complètement le dispositif de stockage d'énergie en moins de cinq heures, l'étape (E2, E3) de détermination de la première fonction et/ou de la deuxième fonction comprenant alors une sous-étape de filtrage de la grandeur relative à une quantité de charges circulant dans un étage. Procédé de détection selon l’une des revendications précédentes, caractérisé en ce que l'étape (E5) de comparaison de ladite différence à un seuil comprend : - une sous-étape (E51) de comparaison de ladite différence à un premier seuil et à un deuxième seuil, le deuxième seuil étant strictement supérieur au premier seuil, puis - une sous-étape (E52) de mémorisation d'un premier témoin d'alerte indiquant un risque modéré, si ladite différence est supérieure ou égale au premier seuil et strictement inférieur au deuxième seuil, et - une sous-étape (E53) de mémorisation d'un deuxième témoin d'alerte indiquant un risque élevé, si ladite différence est supérieure ou égale au deuxième seuil. Procédé de détection selon la revendication précédente, caractérisé en ce que le premier seuil est déterminé en fonction d'une dispersion constatée de ladite différence, et en ce que le deuxième seuil est déterminé en fonction d'une surcharge admissible par au moins un étage de cellules électrochimiques du dispositif de stockage d'énergie. Equipement de surveillance (5) d'un dispositif de stockage d'énergie (1) comprenant un ensemble d’étages (2) électrochimiques reliés électriquement en série, caractérisé en ce qu'il comprend des moyens matériels (6, 7, 8) et logiciels configurés pour mettre en œuvre le procédé de détection d'un risque de défaillance par déséquilibre du dispositif de stockage d'énergie selon l'une des revendications précédentes.