FR 3129735 A1 20230602 FR 2112656 A 20211129 Procédé de diagnostic d’un système électrochimique comprenant : a. collecter (101) des données de mesures horodatées, b. déduire (102) pour chaque instant: - un état de charge de référence (Q t ) du système, - un écart d’état de charge (∆Q t,max ) du premier élément (C t,1 ) présentant la tension la plus élevée, - un écart d’état de charge (∆Q t,min ) du deuxième élément (C t,2 ) présentant la tension la plus faible, c. grouper (103) les jeux de données (Q t , ∆Q t,max , ∆Q t,min , C t,1 , C t,2 ,) par plages de valeurs (T) de l’état de charge de référence (Q t ), d. attribuer (104) à chaque paire (C T,1 , ∆Q T,max ; C T,2 , ∆Q T,min ) de données constituées de l’identifiant de l’élément l’électrochimique (C t,1 , C t,2 ) et de l’écart d’état de charge correspondant (∆Q t,max , ∆Q t,min ) : - une fonction de diagnostic parmi des fonctions préétablies. Figure de l’abrégé : Figure 1 Identification d’anomalie dans un système électrochimique La présente divulgation relève du domaine des dispositifs électrochimiques de stockage d’énergie, et plus particulièrement du diagnostic de batteries, ou accumulateurs. Il est connu de surveiller le comportement et les performances des batteries au fil du temps par exemple en mesurant certaines propriétés d’une batterie pour en déduire son état. L’un des paramètres de performance d’un dispositif électrochimique est une variable appelée « état de santé », ou SOH pour « State Of Health ». Ce paramètre, sans unité et exprimé généralement en pourcentage, se définit comme le rapport de la capacité électrique réelle C rea de l’élément testé sur la capacité électrique nominale C nom (ou la capacité en début de vie) du même élément. L’équation ci-après représente ce paramètre : Un autre paramètre utile à la qualification d’un état d’un dispositif électrochimique est « l’état de charge », ou SOC pour « State Of Charge ». Ce paramètre est aussi sans unité et exprimé généralement en pourcentage. Le SOC se définit comme le rapport de la charge réelle disponible Q t de l’élément testé sur la charge maximale Q max,t du même élément dans son état de santé actuel. L’équation ci-après représente ce paramètre : Il est possible de simuler le fonctionnement de systèmes et éléments électrochimiques en construisant des modèles théoriques. Dans ce cas, la pertinence des résultats virtuels dépend de la qualité de la modélisation. Il est aussi possible de mesurer des propriétés d’un élément électrochimique (réel donc) relatives à son état de santé. Généralement, cela nécessite de lui imposer des sollicitations particulières (cycles de caractérisation). Cela nécessite des équipements de contrôle particuliers et du temps. Les méthodes connues tendent même à dégrader les éléments testés et leur état de santé que l’on souhaite justement évaluer. De manière générale, un élément électrochimique soumis à de tels essais n’est plus opérationnel le temps des essais. Autrement dit, cela est réalisé au cours d’une période de maintenance, ou du moins hors d’un état de fonctionnement normal. Quelques méthodes peuvent être mises en œuvre pendant une phase opérationnelle du dispositif mais seulement dans certaines conditions de fonctionnement, par exemple lorsque la charge et/ou la décharge est complète. D’autres impliquent que des opérateurs entrent des données autres que des données de mesure du dispositif. En pratique, il existe des systèmes de supervision et de contrôle des systèmes électrochimiques, souvent appelés BMS pour « Battery Management System ». Un BMS reçoit en entrées des paramètres et mesures lors du fonctionnement nominal d’un système électrochimique et fournit en sortie des données de diagnostic, par exemple des signaux d’alerte, et quelques données de mesure à destination de niveaux de contrôle supérieurs, souvent appelées « plateforme de supervision ». Typiquement, un BMS est localement associé à un unique système électrochimique et les données de diagnostic issues de plusieurs BMS sont collectés à distance de manière centralisée sur la plateforme de supervision. Cela signifie que, en pratique, au niveau de la plateforme de supervision en aval des BMS, seule une partie des données est disponible, souvent avec retard (pas en temps réel), et non l’intégralité des mesures effectuées et traitées localement par chaque BMS. Par exemple, même si un BMS collecte des mesures à un niveau de détail élevé, distinctement pour chaque élément composant un système électrochimique, les données fournies en sortie du BMS concernent le système dans son ensemble plutôt que chaque élément considéré en tant que tel. Il serait théoriquement possible de modifier les BMS existants pour implémenter de nouveaux procédés de diagnostic. Cependant, cela impliquerait un grand nombre d’interventions directement au niveau de chaque système électrochimique. Il serait théoriquement possible de centraliser à distance et en temps réel l’intégralité des données mesurées localement aux niveaux des systèmes électrochimiques pour implémenter de nouveaux procédés de diagnostic au niveau de la plateforme de supervision. Cependant, le volume de données à transmettre puis à traiter au niveau de la plateforme de supervision impliquerait un surdimensionnement des infrastructures existantes, notamment en termes de réseau et de ressources de calcul. Il existe donc un besoin en diagnostic plus précis que ceux existants et qui ne nécessiterait pas plus de données de mesure que celles disponibles en sortie des BMS existants. En outre, des diagnostics indépendants de ceux fournis par les BMS permettent d’apporter une complémentarité dans la gestion des systèmes. Résumé La présente divulgation vient améliorer la situation. Il est proposé un procédé de diagnostic d’un système électrochimique réel, mis en œuvre par des moyens informatiques, ledit système étant composé d’une pluralité d’éléments électrochimiques électriquement connectés en série. Le procédé comprend : a. collecter, pour une pluralité d’instants, des données de mesures du système, lesdites données étant horodatées et incluant, pour chaque instant, au moins : - une tension et un courant aux bornes du système, - la tension maximale aux bornes d’un premier élément électrochimique parmi les tensions aux bornes de chacun des éléments électrochimiques, associée à un identifiant dudit premier élément électrochimique, - la tension minimale aux bornes d’un deuxième élément électrochimique parmi les tensions aux bornes de chacun des éléments électrochimiques, associée à un identifiant dudit deuxième élément électrochimique, b. déduire, à partir desdites données de mesures, pour chaque instant : - un état de charge de référence du système, - un premier écart d’état de charge du premier élément par rapport à l’état de charge de référence du système, - un deuxième écart d’état de charge du deuxième élément par rapport à l’état de charge de référence du système, de manière à obtenir, pour chaque instant, un jeu de données associé à un couple d’éléments d’intérêt formé par les premier et deuxième éléments, c. grouper les jeux de données par plages de valeurs de l’état de charge de référence de chaque jeu de données, de manière à retenir un jeu de données représentatif par plages de valeurs de l’état de charge de référence, d. attribuer à chaque paire de données constituées de l’identifiant de l’élément l’électrochimique et de l’écart d’état de charge correspondant : - une absence de fonction de diagnostic, ou l’une au moins des fonctions de diagnostic suivantes : - une fonction d’identification d’au moins un élément électrochimique présentant une différence de capacité résiduelle par rapport à une capacité résiduelle nominale des éléments électrochimiques du système dans la plage de valeur de l’état de charge de référence ; - une fonction de quantification de ladite ou lesdites différences de capacité résiduelle ; - une fonction de quantification du déséquilibre maximal entre les éléments chimiques du système. Selon un autre aspect, il est proposé un programme informatique comportant des instructions pour la mise en œuvre de tout ou partie d’un procédé tel que défini dans les présentes lorsque ce programme est exécuté par un processeur. Selon un autre aspect, il est proposé un support d’enregistrement non transitoire, lisible par un ordinateur, sur lequel est enregistré un tel programme. Un tel procédé permet d’effectuer un diagnostic différencié et précis de deux phénomènes physiques nuisibles au bon fonctionnement des systèmes électrochimiques à partir de données de mesure usuelles. Un tel diagnostic permet de mettre en œuvre des mesures correctives ciblées, de manière fiable, pour améliorer les performances ultérieures des systèmes, leur durabilité et la sécurité des équipements et opérateurs voisins. Par « diagnostic différencié », on entend ici qu’il devient possible de détecter et identifier distinctement les phénomènes l’un de l’autre et aussi qu’il est possible d’identifier la source d’une anomalie à l’échelle précise d’un ou plusieurs éléments constituant le système électrochimique plutôt qu’à l’échelle globale et imprécise du système dans son ensemble. Par « données de mesure usuelles », on entend ici des données qui sont classiquement disponibles en fonctionnement opérationnel (non seulement pendant des phases de diagnostic, de maintenance ou expérimentales), y compris en aval des BMS. Autrement dit, lesdites données de mesure peuvent être obtenues sans qu’il soit nécessaire d’imposer un régime de fonctionnement spécifique, ni d’ajouter des équipements de mesure supplémentaire aux systèmes opérationnels existants. Il s’agit par exemple du courant traversant le système et de tensions aux bornes des éléments du système. Les caractéristiques exposées dans les paragraphes suivants peuvent, optionnellement, être mises en œuvre, indépendamment les unes des autres ou en combinaison les unes avec les autres : Le groupement des jeux de données comprend, pour chaque plage de valeurs d’état de charge, le calcul : - d’une première valeur représentative des premiers écarts d’états de charge déduites des données de mesure, ladite première valeur représentative correspondant à une moyenne non pondérée, une moyenne pondéré ou des quantiles, tels qu’une médiane, des premiers écarts d’états de charge, - d’une deuxième valeur représentative des deuxièmes écarts d’états de charge déduites des données de mesure, ladite deuxième valeur représentative correspondant à une moyenne non pondérée, une moyenne pondéré ou une médiane des deuxième écarts d’états de charge, lesdites première et deuxièmes valeur représentative calculées pour chaque plage de valeurs d’état de charge étant celles retenues pour l’attribution de fonctions de diagnostic ou d’absence de fonction de diagnostic. Le nombre d’occurrences des identifiant des éléments électrochimiques parmi les données de mesures collectées est un paramètre d’attribution des fonctions de diagnostique ou d’attribution d’absence de fonction de diagnostic. Chaque paire de données constituées de l’identifiant de l’élément l’électrochimique et de l’écart d’état de charge correspondant peut en outre se voir attribuée la fonction suivante : - une fonction d’identification des éléments impliquées dans le déséquilibre maximal entre les éléments chimiques du système. Le procédé comprend en outre une étape de diagnostic : e. générer l’un au moins des signaux d’anomalie suivant en fonction des paires de données constituée de l’identifiant de l’élément l’électrochimique et de l’écart d’état de charge correspondant : - une identification d’au moins un élément électrochimique présentant une différence de capacité résiduelle par rapport à une capacité résiduelle nominale des éléments électrochimiques du système dans la plage de valeur de l’état de charge de référence ; - une quantification de ladite ou lesdites différences de capacité résiduelle ; - une quantification du déséquilibre maximal entre les éléments chimiques du système. Le procédé comprend en outre une étape de déclenchement d’au moins une mesure corrective en réponse à l’un au moins desdits signaux d’anomalie. Les paires de données constituées de l’identifiant de l’élément l’électrochimique et de l’écart d’état de charge correspondant auxquelles une fonction de quantification de ladite ou lesdites différences de capacité résiduelle est attribuée sont utilisées conjointement dans la construction d’un modèle de régression linéaire, puis une quantification d’une anomalie de capacité résiduelle est mise en œuvre en fonction de la valeur de pente dudit modèle de régression linéaire. Les paires de données constituées de l’identifiant de l’élément l’électrochimique et de l’écart d’état de charge correspondant auxquelles une fonction de quantification du déséquilibre maximal entre les éléments chimiques du système est attribuée sont utilisées conjointement pour calculer une moyenne, pondérée ou non, ou des quantiles, tels qu’une médiane, des valeurs des premiers et deuxièmes écarts d’états de charge. Les jeux de données sont obtenus pour une pluralité d’instants répartis sur une pluralité de cycles de charge et/ou décharge, complets ou partiels, dudit système électrochimique. D’autres caractéristiques, détails et avantages apparaîtront à la lecture de la description détaillée ci-après, et à l’analyse des dessins annexés, sur lesquels : Fig. 1 montre une architecture de gestion d’une pluralité de systèmes électrochimiques. Fig. 2 montre les tensions de plusieurs éléments d’un même système électrochimique au cours d’une décharge du système. Fig. 3 montre les extremums de différence de SOC par rapport à la moyenne de SOC du système parmi les éléments du système de la figure précédente en fonction du SOC du système. Fig. 4 montre les tensions de plusieurs éléments d’un même système électrochimique au cours d’une décharge du système. Fig. 5 montre les extremums de différence de SOC par rapport à la moyenne de SOC du système parmi les éléments du système de la figure précédente en fonction du SOC du système. Fig. 6 montre les tensions de plusieurs éléments d’un même système électrochimique au cours d’une décharge du système. Fig. 7 montre les extremums de différence de SOC par rapport à la moyenne de SOC du système parmi les éléments du système de la figure précédente en fonction du SOC du système. Fig. 8 montre les tensions de plusieurs éléments d’un même système électrochimique au cours d’une décharge du système. Fig. 9 montre les extremums de différence de SOC par rapport à la moyenne de SOC du système parmi les éléments du système de la figure précédente en fonction du SOC du système. Fig. 10 montre un diagramme de mise en œuvre d’un procédé. Les présentes solutions techniques peuvent trouver à s’appliquer sur tous types de batteries ou accumulateurs et notamment aux batteries de type Li-ion stationnaires (pour les réseaux ou le résidentiel par exemple) ou embarquées (véhicules par exemple). En conclusion, les modes de réalisation du procédé décrit ci-avant permettent de ne pas nécessiter de données de mesure propres à chacun des éléments composant un système, mais au contraire de nécessiter seulement des données de mesure plutôt usuelles, en particulier relatives aux valeurs extrêmes plutôt qu’à chaque valeur de chaque élément. Ainsi, et contrairement à de nombreuses méthodes existantes, il n’est pas nécessaire d’effectuer des calculs lourds et un diagnostic absolu pour chaque élément du système pour ensuite faire une comparaison des éléments entre eux. Le procédé décrit ici permet de quantifier une capacité résiduelle en termes relatifs (par rapport à une capacité résiduelle de référence), ce qui suffit dans la plupart des cas à identifier une anomalie. En outre, selon les méthodes connues, l’identification (et non la quantification) d’un élément ayant une capacité résiduelle plus faible/forte que la moyenne consiste à faire une analyse de la fréquence d’apparition des identifiants correspondant aux éléments en tension minimale ou maximale au cours du temps. Cela implique le postulat sous-jacent que, sur une période donnée, l’élément le plus faible sera celui le plus souvent associé à la tension minimum ou maximum de l’assemblage. Or, selon cette approche le résultat est fonction du temps pendant lequel le système se trouve dans chaque état de charge. En effet, en pratique, la cellule associée à la tension minimale n’est pas forcément la même à mi-charge et à pleine charge par exemple. En conséquence, dans une approche basée sur la fréquence d’apparition, les zones d’état de charge où le système passe le plus de temps vont être surpondérées. Au contraire, le procédé proposé ici, en se basant sur une approche dans le plan de l’état de charge plutôt que strictement temporelle, permet une identification et quantification des capacités résiduelles faibles/fortes par tranche d’état de charge : cela est plus pertinent pour aborder un problème lié par nature à la capacité des éléments et donc à leur état de charge. En particulier, les biais de pondération générés par les approches temporelles usuelles sont ici supprimés. Quant à la quantification du déséquilibre en état de charge, il s’agit d’un diagnostic assez rare dans l’art antérieur. En pratique, le circuit d’équilibrage (la plupart du temps un équilibrage passif par dissipation d’énergie dans des résistances) est simplement activé sur critère de seuil de tension. Dans ce cas, aucune quantification préalable du déséquilibre n’est effectuée. Une approche plus avancée repose sur l’exécution de l’équilibrage à la suite d’une quantification du besoin d’équilibrage de chaque cellule. Ainsi, une estimation de l’écart de SOC entre les éléments est effectuée pour chacune des cellules afin d’appliquer un équilibrage adapté à chacune. Une telle approche repose, par exemple, sur un calcul du SOC au niveau d’un élément (une cellule). Cette approche est impossible lorsque des données ne sont pas disponibles pour chaque élément. En outre, les méthodes connues jusque maintenant se basent souvent sur un diagnostic ponctuel à un instant donné. Au contraire, le procédé proposé repose quant à lui sur une quantité de données limitées mais portant sur une durée (un historique), par exemple de plusieurs jours, ce qui permet de tirer profit d’une plus grande variété de situations réelles et d’états du système. Pour autant, le procédé proposé peut tout à fait être mis en œuvre « en temps réel », par exemple sur une fenêtre glissante de temps se terminant à l’instant présent. Enfin, les approches traditionnelles traitent les problématiques de déséquilibre et de capacité résiduelle plus faible de manière indépendante, ce qui peut porter à confusion. En effet, regardé à un instant donné, un écart de tension d’un élément par rapport à la moyenne peut tout aussi bien être expliqué par un élément faible que par un déséquilibre. Le procédé proposé, en étant global et en exploitant un historique de données large permet d’obtenir un diagnostic plus précis et de meilleure qualité. La demanderesse a appliqué le procédé décrit sur un système « saint » d’une part (présentant seulement des déséquilibres) et sur un autre système défaillant (éléments faibles et déséquilibres). Le procédé a été appliqué lors de diverses périodes de temps et dans des états de fonctionnement des systèmes différents les uns des autres sans que cela ne nuise à la bonne répétabilité des résultats. L’analyse du système saint a permis la quantification des déséquilibres hauts et bas avec une répétabilité de l’estimation de l’ordre de ± 0,25% SOC (en termes absolus) sur des fenêtres temporelles indépendantes les unes des autres. L’analyse du système défaillant a permis la quantification des déséquilibres hauts et bas avec une répétabilité de l’estimation de l’ordre de ± 0,5% SOC (en termes absolus) ainsi qu’une quantification de la capacité résiduelle de l’élément faible avec une variabilité de l’ordre de ± 20% (en termes relatifs), ordre de grandeur suffisant pour la détection d’un élément faible. La présente divulgation ne se limite pas aux exemples de procédés, programmes informatiques et support d’enregistrement non transitoire lisible par un ordinateur décrits ci-avant, seulement à titre d’exemple, mais elle englobe toutes les variantes que pourra envisager l’homme de l’art dans le cadre de la protection recherchée. Procédé de diagnostic d’un système électrochimique réel, mis en œuvre par des moyens informatiques, ledit système étant composé d’une pluralité d’éléments électrochimiques électriquement connectés en série, le procédé comprenant : a. collecter (101), pour une pluralité d’instants (t), des données de mesures du système, lesdites données étant horodatées et incluant, pour chaque instant (t), au moins : - une tension (U t ) et un courant (I t ) aux bornes du système, - la tension maximale (U t,max ) aux bornes d’un premier élément électrochimique (C t,1 ) parmi les tensions aux bornes de chacun des éléments électrochimiques, associée à un identifiant dudit premier élément électrochimique (C t,1 ), - la tension minimale (U t,min ) aux bornes d’un deuxième élément électrochimique (C t,2 ) parmi les tensions aux bornes de chacun des éléments électrochimiques, associée à un identifiant dudit deuxième élément électrochimique (C t,2 ), b. déduire (102), à partir desdites données de mesures, pour chaque instant (t) : - un état de charge de référence (Q t ) du système, - un premier écart d’état de charge (∆Q t,max ) du premier élément (C t,1 ) par rapport à l’état de charge de référence (Q t ) du système, - un deuxième écart d’état de charge (∆Q t,min ) du deuxième élément (C t,2 ) par rapport à l’état de charge de référence (Q t ) du système, de manière à obtenir, pour chaque instant (t), un jeu de données (Q t , ∆Q t,max , ∆Q t,min , C t,1 , C t,2 ,) associé à un couple d’éléments d’intérêt formé par les premier et deuxième éléments (C t,1 , C t,2 ), c. grouper (103) les jeux de données (Q t , ∆Q t,max , ∆Q t,min , C t,1 , C t,2 ,) par plages de valeurs (T) de l’état de charge de référence (Q t ) de chaque jeu de données, de manière à retenir un jeu de données représentatif (T, ∆Q T ,max , ∆Q T ,min , C T ,1 , C T ,2 ,) par plages de valeurs (T) de l’état de charge de référence, d. attribuer (104) à chaque paire (C T,1 , ∆Q T,max ; C T,2 , ∆Q T,min ) de données constituées de l’identifiant de l’élément l’électrochimique (C t,1 , C t,2 ) et de l’écart d’état de charge correspondant (∆Q t,max , ∆Q t,min ) : - une absence de fonction de diagnostic, ou l’une au moins des fonctions de diagnostic suivantes : - une fonction d’identification d’au moins un élément électrochimique présentant une différence de capacité résiduelle par rapport à une capacité résiduelle nominale des éléments électrochimiques du système dans la plage de valeur (T) de l’état de charge de référence ; - une fonction de quantification de ladite ou lesdites différences de capacité résiduelle ; - une fonction de quantification du déséquilibre maximal entre les éléments chimiques du système. Procédé selon la revendication 1, dans lequel le groupement des jeux de données comprend, pour chaque plage de valeurs (T) d’état de charge, le calcul : - d’une première valeur représentative (∆Q T,max ) des premiers écarts d’états de charge (∆Q t,max ) déduites des données de mesure, ladite première valeur représentative correspondant à une moyenne non pondérée, une moyenne pondéré ou des quantiles, tels qu’une médiane, des premiers écarts d’états de charge (∆Q t,max ), - d’une deuxième valeur représentative (∆Q T,min ) des deuxièmes écarts d’états de charge (∆Q t,min ) déduites des données de mesure, ladite deuxième valeur représentative correspondant à une moyenne non pondérée, une moyenne pondéré ou une médiane des deuxième écarts d’états de charge (∆Q t,m in ), lesdites première et deuxièmes valeur représentative (∆Q T,max , ∆Q T,min ) calculées pour chaque plage de valeurs (T) d’état de charge étant celles retenues pour l’attribution de fonctions de diagnostic ou d’absence de fonction de diagnostic. Procédé selon l’une des revendications précédentes, dans lequel le nombre d’occurrences des identifiant (C t,1 , C t,2 ,) des éléments électrochimiques parmi les données de mesures collectées est un paramètre d’attribution des fonctions de diagnostique ou d’attribution d’absence de fonction de diagnostic. Procédé selon l’une des revendications précédentes, dans lequel chaque paire (C T,1 , ∆Q T,max ; C T,2 , ∆Q T,min ) de données constituées de l’identifiant de l’élément l’électrochimique (C t,1 , C t,2 ) et de l’écart d’état de charge correspondant (∆Q t,max , ∆Q t,min ) peut en outre se voir attribuée la fonction suivante : - une fonction d’identification des éléments impliquées dans le déséquilibre maximal entre les éléments chimiques du système. Procédé selon l’une des revendications précédentes comprenant en outre une étape de diagnostic : e. générer l’un au moins des signaux d’anomalie suivant en fonction des paires (C T,1 , ∆Q T,max ; C T,2 , ∆Q T,min ) de données constituée de l’identifiant de l’élément l’électrochimique (C t,1 , C t,2 ) et de l’écart d’état de charge correspondant (∆Q t,max , ∆Q t,min ) : - une identification d’au moins un élément électrochimique présentant une différence de capacité résiduelle par rapport à une capacité résiduelle nominale des éléments électrochimiques du système dans la plage de valeur (T) de l’état de charge de référence ; - une quantification de ladite ou lesdites différences de capacité résiduelle ; - une quantification du déséquilibre maximal entre les éléments chimiques du système. Procédé selon la revendication 5 comprenant en outre une étape de déclenchement d’au moins une mesure corrective en réponse à l’un au moins desdits signaux d’anomalie. Procédé selon l’une des revendications précédentes, dans lequel les paires (C T,1 , ∆Q T,max ; C T,2 , ∆Q T,min ) de données constituées de l’identifiant de l’élément l’électrochimique (C t,1 , C t,2 ) et de l’écart d’état de charge correspondant (∆Q t,max , ∆Q t,min ) auxquelles une fonction de quantification de ladite ou lesdites différences de capacité résiduelle est attribuée sont utilisées conjointement dans la construction d’un modèle de régression linéaire, puis une quantification d’une anomalie de capacité résiduelle est mise en œuvre en fonction de la valeur de pente dudit modèle de régression linéaire. Procédé selon l’une des revendications précédentes, dans lequel les paires (C T,1 , ∆Q T,max ; C T,2 , ∆Q T,min ) de données constituées de l’identifiant de l’élément l’électrochimique (C t,1 , C t,2 ) et de l’écart d’état de charge correspondant (∆Q t,max , ∆Q t,min ) auxquelles une fonction de quantification du déséquilibre maximal entre les éléments chimiques du système est attribuée sont utilisées conjointement pour calculer une moyenne, pondérée ou non, ou des quantiles, tels qu’une médiane, des valeurs des premiers et deuxièmes écarts d’états de charge (∆Q t,max , ∆Q t,min ). Procédé selon l’une des revendications précédentes, dans lequel lesdits jeux de données sont obtenus pour une pluralité d’instants (t) répartis sur une pluralité de cycles de charge et/ou décharge, complets ou partiels, dudit système électrochimique. Programme informatique comportant des instructions pour la mise en œuvre du procédé selon l’une des revendications précédentes lorsque ce programme est exécuté par un processeur. Support d’enregistrement non transitoire lisible par un ordinateur sur lequel est enregistré un programme pour la mise en œuvre du procédé selon l’une des revendications 1 à 9 lorsque ce programme est exécuté par un processeur.