La présente invention concerne un procédé pour estimer l'état de santé d'une cellule électrochimique de stockage d'énergie électrique. Elle s'applique notamment, mais pas exclusivement, aux véhicules automobiles électriques ou hybrides.

Dans le contexte actuel de consensus autour du réchauffement climatique, la diminution des émissions de dioxyde de carbone (CO2) est un défi majeur auquel sont confrontés les constructeurs automobiles, les normes étant toujours plus exigeantes en la matière.

Outre l'amélioration constante des rendements des moteurs thermiques classiques, qui s'accompagne d'une baisse des émissions de CO2, les véhicules électriques (« EV » d'après la terminologie anglo-saxonne « Electric Vehicle ») et les véhicules hybrides thermique-électrique (« HEV » d'après la terminologie anglo-saxonne « Hybrid Electric Vehicle ») sont aujourd'hui considérés comme la solution la plus prometteuse pour diminuer les émissions de CO2.

Différentes technologies de stockage de l'énergie électrique ont été testées dans les dernières années afin de répondre aux besoins des EV. Il apparaît aujourd'hui que les batteries à cellules lithium-ion (Li-ion) sont celles qui permettent d'obtenir le meilleur compromis entre la densité de puissance, qui favorise les performances en termes d'accélération notamment, et la densité d'énergie, qui favorise l'autonomie.

Malheureusement, la densité de puissance et la densité d'énergie diminuent tout au long de la vie de la batterie, notamment sous l'effet des variations de température. Ainsi, après un délai d'utilisation suffisamment long, un EV peut présenter des performances dégradées en termes d'autonomie et/ou de puissance. Il convient de contrôler cette dégradation, afin de maintenir des niveaux de prestation et de sécurité suffisants.

Afin de quantifier cette dégradation, un indicateur appelé « Etat de Santé » (« SOH » d'après la terminologie anglo-saxonne « State Of Health ») a été défini, qui est le ratio de la capacité courante d'une batterie par rapport à sa capacité initiale en début de vie. Il est connu, notamment du brevet US6653817, d'estimer le SOH à partir de l'estimation par impédancemétrie de la résistance interne de la batterie (ou « DCR » d'après la terminologie anglo-saxonne « Direct Current Résistance »). La résistance interne d'une batterie caractérisant la variation de la tension U aux bornes de la batterie pour une certaine variation de l'intensité I du courant traversant la batterie, ce brevet propose de piloter finement la variation de l'intensité I traversant la batterie et de mesurer la variation de la tension U aux bornes de la batterie. Un inconvénient majeur de cette solution est que si la mesure de tension n'est pas précise, notamment en cas de bruits transitoires, l'estimation de la résistance peut être impactée et l'estimation du SOH peut être imprécise. Il s'agit de l'un des problèmes que la présente invention se propose de résoudre.

L'invention a notamment pour but de résoudre les inconvénients précités, notamment d'éviter les problèmes de mesure, ceci en linéarisant la mesure de la résistance interne. A cet effet, l'invention a pour objet un procédé pour estimer l'état de santé d'une cellule électrochimique de stockage d'énergie électrique. Il comporte une étape d'application à la cellule d'au moins un pic de courant d'intensité Ι _γ , le pic de courant traversant la cellule. Il comporte également une étape de mesure de la variation, en fonction du temps t s'écoulant après l'application du pic de courant, de la tension U aux bornes de la cellule. Il comporte également une étape de calcul d'au moins un coefficient a ₇ et d'au moins un coefficient U _0J tels que la fonction t→ x [t + U _{0 Ii} est une approximation linéaire de la variation de la tension U en fonction de t pour t > C , où C > 0 .

Dans un mode de réalisation, les coefficients a ₇ et U _0J peuvent être calculés une première fois en début de vie de la cellule. Le procédé peut alors comporter une étape de calcul du ratio de la valeur couramment calculée de a par rapport à sa valeur calculée en début de vie, une augmentation du coefficient a au-delà d'une valeur prédéterminée pouvant indiquer une incapacité de la cellule à délivrer un courant dans ses gammes de puissance les plus élevées. Le procédé peut également comporter une étape de calcul du ratio de la valeur couramment calculée de U ₀ par rapport à sa valeur calculée en début de vie, une augmentation du coefficient U _{0 Ii} au-delà d'une valeur prédéterminée pouvant indiquer une incapacité de la cellule à délivrer un courant dans ses gammes de puissance les plus élevées.

Dans un autre mode de réalisation, des tables ou des graphes de vieillissement peuvent être remplies au préalable. Le procédé peut alors comporter une étape de comparaison de la valeur couramment calculée de a avec des valeurs contenues dans une table ou un graphe de vieillissement associant des niveaux de vieillissement à des valeurs de a , de manière à déduire le niveau de vieillissement de la cellule. Le procédé peut également comporter une étape de comparaison de la valeur couramment calculée de U _{0 Ii} avec des valeurs contenues dans une table ou un graphe de vieillissement associant des niveaux de vieillissement à des valeurs de U _{0 Ii} , de manière à déduire le niveau de vieillissement de la cellule.

Dans un mode de réalisation préférentiel, une pluralité de pics de courant d'intensités (l„) _n≥2 peuvent être appliqués à la cellule

( / _j , I ₂ , I ₃ , 1 ₄ , I ₅ )■ Le procédé peut alors comporter une étape de calcul d'un coefficient β tel que la fonction I→ β χ ΐ est une approximation linéaire de la variation de ^Δ ^ ^ ^en fonction de / .

Dans un autre mode de réalisation, le coefficient β peut être calculé une première fois en début de vie de la cellule. Le procédé peut alors comporter une étape de calcul du ratio de la valeur couramment calculée de β par rapport à sa valeur calculée en début de vie, une augmentation du coefficient β au-delà d'une valeur prédéterminée pouvant indiquer une incapacité de la cellule à délivrer un courant dans ses gammes de puissance les plus élevées.

Dans un autre mode de réalisation préférentiel, des tables ou des graphes de vieillissement peuvent être remplies au préalable. Le procédé peut alors comporter une étape de comparaison de la valeur couramment calculée de β avec des valeurs contenues dans une table ou un graphe de vieillissement associant des niveaux de vieillissement à des valeurs de β , de manière à déduire le niveau de vieillissement de la cellule. Dans un mode de réalisation préférentiel, le procédé peut comporter une étape de calcul d'un coefficient y tel que la fonction

I→y x I + OCV , où OCV est la tension en circuit ouvert de la cellule

(« Open Circuit Voltage » selon la terminologie anglo-saxonne), est une approximation linéaire de la variation de U _{0 I} en fonction de / .

Dans un mode de réalisation, le coefficient y peut être calculé une première fois en début de vie de la cellule. Le procédé peut alors comporter une étape de calcul du ratio de la valeur couramment calculée de y par rapport à sa valeur calculée en début de vie, une augmentation du coefficient y au-delà d'une valeur prédéterminée pouvant indiquer une incapacité de la cellule à délivrer un courant dans ses gammes de puissance les plus élevées.

Dans un mode de réalisation préférentiel, des tables ou des graphes de vieillissement peuvent être remplies au préalable. Le procédé peut alors comporter une étape de comparaison de la valeur couramment calculée de y avec des valeurs contenues dans une table ou un graphe de vieillissement associant des niveaux de vieillissement à des valeurs de y , de manière à déduire le niveau de vieillissement de la cellule.

La présente invention a encore pour principal avantage qu'elle ne nécessite que la mise à jour logicielle des dispositifs actuels d'estimation de l'état de santé d'une batterie.

D'autres caractéristiques et avantages de l'invention apparaîtront à l'aide de la description qui suit faite en regard de dessins annexés qui représentent :

- la figure 1 , par un graphe, un exemple de variation en fonction du temps t après application d'un pic de courant d'intensité de la tension U aux bornes d'une cellule Li-ion ;

- la figure 2, par un graphe, la variation de U en fonction de y[t ;

- la figure 3, par un graphe, différents exemples de variation de U en fonction de y[t pour différents pics de courant ; - la figure 4, par un graphe, la variation de ^Δ ¾ en fonction de l'intensité I ;

la figure 5, par un graphe, la variation en fonction de coefficient selon l'invention ;

- les figures 6, 7, 8 et 9, des graphes de vieillissement selon l'invention.

La présente invention propose de calculer différents coefficients caractéristiques de l'état de santé d'une cellule Li-ion.

Comme illustré dans l'exemple de la figure 1 , la tension U aux bornes d'une cellule Li-ion ayant un niveau de charge de l'ordre de 40% peut, après qu'un pic de courant d'intensité l _x = ! (ampères) ait traversé cette cellule à un instant t=0, évoluer de t=0 à t=60 secondes conformément au graphe illustré par cette figure. A partir de cette simple courbe représentant la variation de U en fonction du temps t , il est possible selon l'invention d'extraire divers coefficients directement reliés à la capacité de la cellule à fonctionner à forte densité de puissance, cette capacité étant désignée « SOH _P » dans la suite de la présente demande, par opposition à sa capacité à fonctionner à forte densité d'énergie qui sera désignée « SOH _E ».

Dans l'exemple de la figure 1 , la tension aux bornes de la cellule décroît de 3,86V à t=0s jusqu'à 3,71V à t=60s de manière clairement non linéaire. La valeur initiale de tension à 3,86V correspondant à la tension aux bornes de la cellule en circuit ouvert, c'est-à-dire lorsque la cellule n'est pas traversée par un courant, cette tension étant couramment désignée par l'acronyme anglo-saxon « OCV » signifiant « Open Circuit Voltage ».

La figure 2 illustre, en fonction de y[t et non plus en fonction de t, l'évolution de la même tension U après que le pic de courant d'intensité ait traversé la cellule. On constate que, à partir de t = 1 , la tension décroît linéairement en fonction de y[t . Et en effet, une méthode de linéarisation montre que la variation de U en fonction de t peut être approchée, pour > 1 , par la droite d'équation U = a x t + U ₀ de pente a = -0,011 et d'ordonnée à l'origine U ₀ = 3,7936 , ceci avec un coefficient de corrélation

R ² = 0,9984 .

Pour un courant I donné, la pente a ₇ donne une indication sur le

SOHp de la cellule, notamment sur sa capacité à fonctionner sur des temps longs supérieurs à 1 seconde à un courant donné.

Pour un courant I donné, l'ordonnée à l'origine U _{0 I} donne également une indication sur le SOH _P de la cellule, en particulier :

• plus U _{0 I} est important par rapport à l'OCV, moins la cellule est apte à fonctionner à haute puissance de l'ordre de 100 à 1000 watts.; · une ordonnée à l'origine U _{0 I} important par rapport à l'OCV peut aussi indiquer que la cellule a un problème à des temps courts, c'est- à-dire à hautes fréquences, comme des problèmes de connectique ou de soudure.

Optionnellement, il est possible de mettre en œuvre une stratégie de diagnostic de l'état de santé à plusieurs courants selon la présente invention, comme illustré par les figure 3, 4 et 5 qui suivent, cette stratégie permettant un diagnostic plus affiné.

En effet, la figure 3 illustre en fonction de Vt l'évolution de la tension U aux bornes de la même cellule ayant un niveau de charge de l'ordre de 40%, non seulement après que le pic de courant d'intensité 7^ = 74 A ait traversé la cellule, mais également après que 4 autres pics d'intensités respectives 7 ₂ = 37 A, 7 ₃ = 18,5A, 7 ₄ = 7,4A et 7 ₅ = 3,7A aient traversé cette même cellule. Là encore, quel que soit le pic considéré, on constate que, à partir de Vt = 1 , la tension décroît linéairement en fonction de y[t . La même méthode de linéarisation montre que les variations des tensions en fonction de Vt peuvent être approchées, pour t≥ 1 :

• après le pic d'intensité , par la droite d'équation U = a x t + U ₀ de pente a = -0,0058 et d'ordonnée à l'origine

U ₀ = 3,8268 , ceci avec un coefficient de corrélation R ² = 0,9981 ; après le pic d'intensité h, par la droite d'équation U = a x t +U _0I} de pente a = -0,003 et d'ordonnée à l'origine

U ₀ = 3,8424, ceci avec un coefficient de corrélation R ² = 0,9989 ; après le pic d'intensité l ₄ , par la droite d'équation U = a _l4 yft + U _0l4 de pente a _h =-0,0012 et d'ordonnée à l'origine

U _0l4 = 3,8515 , ceci avec un coefficient de corrélation R ² = 0,9938 ; après le pic d'intensité l ₅ , par la droite d'équation U = a x t +U _0Is de pente a = -0,0007 et d'ordonnée à l'origine

U _0I = 3,8546 , ceci avec un coefficient de corrélation R ² = 0,9652.

La figure 4 illustre, d'une part, l'évolution de ^Δ ^ ^ ^en fonction de l'intensité I, ceci pour les valeurs I _x =74 A, I ₂ =37 A, I ₃ =18,5 A, I ₄ =7,4A et I ₅ =3,7A de la figure 3. On constate que ^Δ ^ ^ décroît linéairement avec I. Et en effet, une méthode de linéarisation montre que la variation de ^Δ ^ ^ en fonction de / peut être approchée par la droite d'équation ^Δ ^ ^ = β ^{x 1} de pente β =-0,000151 et d'ordonnée nulle à l'origine, ceci avec un coefficient de corrélation R ² = 0,998.

La pente β révèle la sensibilité de la cellule au courant sur des constantes de temps importantes, c'est-à-dire des phénomènes de diffusion. La pente β indique la capacité de la cellule à fonctionner à des courants importants : plus la valeur absolue de la pente β augmente, plus la cellule est sensible à l'utilisation de forts courants.

La figure 5 illustre, d'autre part, l'évolution de l'ordonnée à l'origine U _0I des courbes illustrées par la figure 3, en fonction de l'intensité I, ceci pour les valeurs I =74 A, I ₂ =37 A, I ₃ =18,5 A, I ₄ =7,4 A et I ₅ =3,7 A de la figure 3. On constate que U _0I décroît linéairement avec I. Et en effet, une méthode de linéarisation montre que la variation de U _0I en fonction de / peut être approchée par la droite d'équation U _0I =γ ^χ ΐ + OCV de pente γ =-0,000868 et d'ordonnée à l'origine OCT = 3,86, ceci avec un coefficient de corrélation R ² = 1 . La pente y correspond quant à elle à la résistance interne de la cellule, qui vaut donc 0,868 milliohms (ητιΩ) à ce niveau de charge de 40% et qui correspond à une impédance mesurée à 3,5 hertz (Hz).

Puisqu'elle correspond à la résistance interne de la cellule, la pente y donne donc également des indications sur le SOH _P de la cellule : plus la pente γ augmente, plus le SOH _P de la cellule est dégradé. De manière plus précise que l'ordonnée à l'origine U _{0 I} en fonction de l'intensité

I, la pente γ estimée à partir des U _{0 I} donne une indication sur la capacité de la cellule à fonctionner à des temps courts, c'est-à-dire à hautes fréquences.

Elle donne une indication sur la résistance à hautes fréquences de la cellule qui, si elle est élevée, peut s'expliquer par un problème de connectique ou de vieillissement important de la cellule.

Dans ce dernier cas, les coefficients a ₇ et β devraient indiquer un vieillissement. Donc si les coefficients a ₇ et β sont acceptables et que le coefficient γ ne l'est pas, on peut déduire que le problème à hautes fréquences est dû à un problème de connectique.

Une fois les coefficients α ₇ , β , U _{0 I} et y calculés conformément à la présente invention, on peut les exploiter de différentes manières.

Comme explicité précédemment, une première manière est de les utiliser à des fins de diagnostic de la cellule, afin notamment d'estimer sa capacité à fonctionner en puissance, c'est-à-dire d'estimer son SOH _P , voire même de diagnostiquer un défaut de connectique. Par exemple, pour α ₇ , β et y , on peut calculer un ratio entre la valeur calculée à l'instant courant et la valeur calculée initialement, à savoir a _{I B0L} , β _ΒΟι et /^ respectivement où l'abbréviation anglo-saxonne « BOL » signifie « Beginning Of Life ». On peut ainsi observer l'évolution relative des ratios ^a y , _R et y : si un

/ ^α /,ΒΟ / P BOL / Y BOL coefficient à un moment donné augmente trop par rapport à sa valeur initiale, tandis que les autres coefficients, eux, montrent une évolution attendue, alors la cellule présente sans doute un défaut de connectique. Il est également possible d'observer des ratios du type y _R ou y . Dans l'exemple illustré par les figures, au cours de la vie de la cellule, le ratio

varie entre 4,59 et 5,78. Cependant, un défaut de connectique de 0,2ηηΩ fait passer cette variation entre 5,78 et 6,814 tandis qu'un défaut de I mQ fait passer cette variation entre 10,56 et 1 1 ,17. En estimant au préalable ces différentes valeurs, il est possible de détecter des problèmes de connectiques en début de vie et lors de la vie de la cellule.

Une autre manière de les exploiter est d'estimer le SOH _E de la cellule en utilisant des tables ou des graphes de vieillissement, comme les graphes illustrés aux figures 6, 7, 8 et 9 pour la pente a ₇ , la pente β , l'ordonnée à l'origine U _{0 I} et la pente y respectivement. Ces graphes associent, à des valeurs desdits coefficients, des valeurs dans une unité arbitraire (u.a.), ces valeurs dans une unité arbitraire caractérisant le vieillissement d'une cellule. Par exemple, une valeur 0 en abscisse caractérise le début de la vie de la cellule et une valeur 9 caractérise la fin de la vie de la cellule. Ces graphes sont remplis préalablement à l'utilisation de la cellule, lors de campagnes d'étude du processus de vieillissement des cellules tout au long de leur vie. Ainsi, d'après la figure 6, la pente ₁ croît de sensiblement 0,0062 en début de vie de la cellule jusqu'à sensiblement 0,0076 en fin de vie de la cellule. D'après la figure 7, la pente β croît de sensiblement 0,000163 en début de vie de la cellule jusqu'à sensiblement

0,000194 en fin de vie de la cellule. D'après la figure 8, l'ordonnée à l'origine U _{0 I} croît de sensiblement 3,829 en début de vie de la cellule jusqu'à sensiblement 3,848 en fin de vie de la cellule. Enfin, d'après la figure 9, la pente y décroît d'abord de sensiblement 0,00082 à sensiblement 0,00077 en début de vie de la cellule, avant de croître jusqu'à sensiblement 0,001 12 en fin de vie de la cellule.

Une autre manière d'exploiter les coefficients ₁ , β , U _{0 I} et y calculés conformément à la présente invention est de comparer une cellule à une autre dans un module ou un pack comportant plusieurs cellules, ou encore de manière à comparer l'évolution de divers types de cellules dans le cas de cellules basées sur des chimies différentes ou n'étant pas issues du même fournisseur. L'invention décrite précédemment a encore pour principal avantage que, ne nécessitant qu'une mise à jour logicielle des dispositifs actuels d'estimation de l'état de santé, elle a un coût de déploiement très faible.