Brevet FR3006450A1 — publié le 20141205

Procédé pour estimer l'état de santé d'une cellule électrochimique de stockage d'énergie électrique 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 5 électriques ou hybrides. Dans le contexte actuel de consensus autour du réchauffement climatique, la diminution des émissions de dioxyde de carbone (002) est un 10 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 002, les véhicules électriques (« EV » d'après la terminologie anglo-saxonne 15 « 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 002. Différentes technologies de stockage de l'énergie électrique ont 20 é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. 25 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, 30 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 Resistance »). 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 'H) 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. 15 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 20 stockage d'énergie électrique. Il comporte une étape d'application à la cellule d'au moins un pic de courant d'intensité 1, , 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 25 calcul d'au moins un coefficient al, et d'au moins un coefficient U0,11 tels que la fonction N/7 x .%/7 + U0,11 est une approximation linéaire de la variation de la tension U en fonction de N/7 pour -Nfi C, où C >0 O. Dans un mode de réalisation, les coefficients al, et U0,1, peuvent être calculés une première fois en début de vie de la cellule. Le procédé peut 30 alors comporter une étape de calcul du ratio de la valeur couramment calculée de al, par rapport à sa valeur calculée en début de vie, une augmentation du coefficient al, 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 35 étape de calcul du ratio de la valeur couramment calculée de U0,1, par rapport à sa valeur calculée en début de vie, une augmentation du coefficient U0,11 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 al, avec des valeurs contenues dans une table ou un graphe de vieillissement associant des niveaux de vieillissement à des valeurs de al, , to 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 U0,11 avec des valeurs contenues dans une table ou un graphe de vieillissement associant des niveaux de vieillissement à des valeurs de U0,11, de manière à déduire le niveau de vieillissement de la 15 cellule. Dans un mode de réalisation préférentiel, une pluralité de pics de courant d'intensités (/),2 peuvent être appliqués à la cellule ( , 12 , 13 , 14 , 15 ). Le procédé peut alors comporter une étape de calcul d'un coefficient ,6 tel que la fonction / x / est une approximation linéaire de 20 la variation de AU/,- en fonction de I. .Vt Dans un autre mode de réalisation, le coefficient ,6 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 ,6 par rapport à sa valeur calculée en début de vie, une augmentation du 25 coefficient ,6 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é 30 peut alors comporter une étape de comparaison de la valeur couramment calculée de ,6 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 7 tel que la fonction / -> 7 x / + OCV , où OCV est la tension en circuit ouvert de la cellule (« Open Circuit Voltage » selon la terminologie anglo-saxonne), est une 5 approximation linéaire de la variation de U0,1 en fonction de I. Dans un mode de réalisation, le coefficient 7 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 7 par rapport à sa valeur calculée en début de vie, une augmentation du coefficient o 7 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é 15 peut alors comporter une étape de comparaison de la valeur couramment calculée de 7 avec des valeurs contenues dans une table ou un graphe de vieillissement associant des niveaux de vieillissement à des valeurs de 7, de manière à déduire le niveau de vieillissement de la cellule. 20 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. 25 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 30 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 la figure 3, par un graphe, différents exemples de variation de U en fonction de -N/7 pour différents pics de courant ; la figure 4, par un graphe, la variation de Aty,- en fonction de l'intensité I ; la figure 5, par un graphe, la variation en fonction de / d'un 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 10 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é /, =74A (ampères) ait traversé cette cellule à un instant t=0, évoluer de t=0 à t=60 secondes conformément 15 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 « SOHp » dans la suite de la présente demande, par opposition à 20 sa capacité à fonctionner à forte densité d'énergie qui sera désignée « SOHE ». Dans l'exemple de la figure 1, la tension aux bornes de la cellule décroît de 3,86V à t=Os jusqu'à 3,71V à t=60s de manière clairement non linéaire. La valeur initiale de tension à 3,86V correspondant à la tension aux 25 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 ». 30 La figure 2 illustre, en fonction de N/7 et non plus en fonction de t, l'évolution de la même tension U après que le pic de courant d'intensité l ait traversé la cellule. On constate que, à partir de V7=1, la tension décroît linéairement en fonction de .N/7 . Et en effet, une méthode de linéarisation montre que la variation de U en fonction de V7 peut être approchée, pour 35 -'/7 1, par la droite d'équation U =cx1 x + um, de pente al, = -0,011 et d'ordonnée à l'origine U0,11 = 3,7936, ceci avec un coefficient de corrélation R2 = 0,9984. Pour un courant 1 donné, la pente al donne une indication sur le SOHp de la cellule, notamment sur sa capacité à fonctionner sur des temps 5 longs supérieurs à 1 seconde à un courant donné. Pour un courant 1 donné, l'ordonnée à l'origine U0,1 donne également une indication sur le SOHp de la cellule, en particulier : - plus U0,1 est important par rapport à l'OCV, moins la cellule est apte à fonctionner à haute puissance de l'ordre de 100 à 1000 watts.; 10 - une ordonnée à l'origine U0,1 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. 15 Optionnellement, il est possible de mettre en oeuvre 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é. 20 En effet, la figure 3 illustre en fonction de .,/7 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é 25 1 = 74 A ait traversé la cellule, mais également après que 4 autres pics d'intensités respectives 12 = 37 A, /3 = 18,5 A, /4 =7,4A et 15 = 3,7 A aient traversé cette même cellule. Là encore, quel que soit le pic considéré, on constate que, à partir de -%/7 = 1, la tension décroît linéairement en fonction de -%/7 . La même méthode de linéarisation montre que les variations des 30 tensions en fonction de ,/7 peuvent être approchées, pour V7 1 : - après le pic d'intensité 12, par la droite d'équation U = x,1 ±(10,1, de pente al, = -0,0058 et d'ordonnée à l'origine U0,12 = 3,8268, ceci avec un coefficient de corrélation R2 = 0,9981 ; - après le pic d'intensité 13, par la droite d'équation U = cx13 xI 7 +U0,1, de pente al, = -0,003 et d'ordonnée à l'origine U0,13 = 3,8424, ceci avec un coefficient de corrélation R2 = 0,9989 ; - après le pic d'intensité 14, par la droite d'équation U = cx14x -\17 +U0,14 de pente al, = -0,0012 et d'ordonnée à l'origine U0,14 = 3,8515, ceci avec un coefficient de corrélation R2 = 0,9938 ; - après le pic d'intensité 15, par la droite d'équation U = x +U0,1, de pente al, = -0,0007 et d'ordonnée à l'origine U0,15 = 3,8546, ceci avec un coefficient de corrélation R2 = 0,9652. La figure 4 illustre, d'une part, l'évolution de AU4 en fonction de l'intensité I, ceci pour les valeurs /, = 74 A, 12 = 37 A, /3 = 18,5A, /4 =7,4A et 15 = 3,7A de la figure 3. On constate que AU4 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 AU/,- = 6x / de pente fi -0,000151 et d'ordonnée nulle à l'origine, ceci avec un coefficient de corrélation R2 = 0,998. La pente fi 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 fi indique la capacité de la cellule à fonctionner à des courants importants : plus la valeur absolue de la pente fi 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 U0,1 des courbes illustrées par la figure 3, en fonction de l'intensité I, ceci pour les valeurs /, = 74 A, 12 = 37 A, /3 = 18,5A, /4 = 7,4 A et 15 = 3,7 A de la figure 3. On constate que U0,1 décroît linéairement avec I. Et en effet, une méthode de linéarisation montre que la variation de U0,1 en fonction de / peut être approchée par la droite d'équation U0,1 = x / + OCV de pente 7 = -0,000868 et d'ordonnée à l'origine OCV = 3,86, ceci avec un coefficient de corrélation R2 =1. La pente 7 correspond quant à elle à la résistance interne de la cellule, qui vaut donc 0,868 milliohms (mn) à 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 7 donne donc également des indications sur le SOHp de la cellule : plus la pente y augmente, plus le SOHp de la cellule est dégradé. De manière plus précise que l'ordonnée à l'origine U0,1 en fonction de l'intensité I, la pente y estimée à partir des U0,1 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 al et ,6 devraient indiquer 15 un vieillissement. Donc si les coefficients al et ,6 sont acceptables et que le coefficient y ne l'est pas, on peut déduire que le problème à hautes fréquences est dû à un problème de connectique. 20 Une fois les coefficients al , 6, U0,1 et 7 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 SOHp, voire 25 même de diagnostiquer un défaut de connectique. Par exemple, pour al , fi et 7 , on peut calculer un ratio entre la valeur calculée à l'instant courant et la valeur calculée initialement, à savoir cxJBQL, PBOL et 7 Bol, respectivement où l'abbréviation anglo-saxonne « BOL » signifie « Beginning Of Life ». On peut ainsi observer l'évolution relative des ratios°y /6/7 et : si un a ,BOL BOL 7 BOL 30 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 2/6 ou 7 Dans l'exemple illustré par les figures, au cours de la vie de la cellule, le ratio 2/6 varie entre 4,59 et 5,78. Cependant, un défaut de connectique de 0,2mS2 fait passer cette variation entre 5,78 et 6,814 tandis qu'un défaut de 1mS2 fait passer cette variation entre 10,56 et 11,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 SOHE 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 an la pente l'ordonnée à l'origine U0,1 et la pente 7 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 a, 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 fi 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 U0,1 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 7 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,00112 en fin de vie de la cellule. Une autre manière d'exploiter les coefficients an U0,1 et 7 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.5 REVENDICATIONS1. Procédé pour estimer l'état de santé d'une cellule électrochimique de stockage d'énergie électrique, le procédé comportant : - une étape d'application à la cellule d'au moins un pic de courant d'intensité 1, , le pic de courant traversant la cellule ; - 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 ; le procédé étant caractérisé en ce qu'il comporte une étape de calcul d'au moins un coefficient ail et d'au moins un coefficient U0,11 tels que la fonction %fi -> cx11 x +U0,11 est une approximation linéaire de la variation de la tension U en fonction de -'/7 pour .^/7 C, où C > o. 2. Procédé selon la revendication 1, caractérisé en ce que, les coefficients ail et U0,11 ayant été calculés une première fois en début de vie de la cellule, le procédé comporte : - 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 ail au-delà d'une valeur prédéterminée dudit ratio indiquant une incapacité de la cellule à délivrer un courant dans ses gammes de puissance les plus élevées, et/ou ; - une étape de calcul du ratio de la valeur couramment calculée de U0,11 par rapport à sa valeur calculée en début de vie, une augmentation du coefficient U0,11 au-delà d'une valeur prédéterminée dudit ratio indiquant une incapacité de la cellule à délivrer un courant dans ses gammes de puissance les plus élevées. 3. Procédé selon la revendication 1, caractérisé en ce que, des tables ou des graphes de vieillissement ayant été remplies au préalable, le procédé comporte : - une étape de comparaison de la valeur couramment calculée de ail avec des valeurs contenues dans une table ou un graphe de vieillissement associant des niveaux de vieillissement à des valeursde al, , de manière à déduire le niveau de vieillissement de la cellule, et/ou ; - une étape de comparaison de la valeur couramment calculée de U0,11 avec des valeurs contenues dans une table ou un graphe de vieillissement associant des niveaux de vieillissement à des valeurs de U0,11, de manière à déduire le niveau de vieillissement de la cellule. 4. Procédé selon la revendication 1, caractérisé en ce que, une pluralité de pics de courant d'intensités (/),2 étant appliqués à la cellule (11, 12 , 13 , 14 , 15 ), il comporte en outre une étape de calcul d'un coefficient fi tel que la fonction / -> fi x I est une approximation linéaire de la variation de Aty,- en fonction de I. -%/t 5. Procédé selon la revendication 4, caractérisé en ce que, le coefficient fi ayant été calculé une première fois en début de vie de la cellule, le procédé comporte une étape de calcul du ratio de la valeur couramment calculée de fi par rapport à sa valeur calculée en début de vie, une augmentation du coefficient fi au-delà d'une valeur prédéterminée dudit ratio indiquant une incapacité de la cellule à délivrer un courant dans ses gammes de puissance les plus élevées. 6. Procédé selon la revendication 4, caractérisé en ce que, des tables ou des graphes de vieillissement ayant été remplies au préalable, le procédé comporte, une étape de comparaison de la valeur couramment calculée de ,6 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. 7. Procédé selon la revendication 4, caractérisé en ce qu'il comporte en outre une étape de calcul d'un coefficient 7 tel que la fonction / 7 x / + OCV , où OCV est la tension en circuit ouvert de la cellule, est une approximation linéaire de la variation de U0,1 en fonction de I. 8. Procédé selon la revendication 7, caractérisé en ce que, le coefficient 7 ayant été calculé une première fois en début de vie de la cellule, leprocédé comporte une étape de calcul du ratio de la valeur couramment calculée de 7 par rapport à sa valeur calculée en début de vie, une augmentation du coefficient 7 au-delà d'une valeur prédéterminée dudit ratio indiquant une incapacité de la cellule à délivrer un courant dans ses gammes de puissance les plus élevées. 9. Procédé selon la revendication 8, caractérisé en ce que, des tables ou des graphes de vieillissement ayant été remplies au préalable, le procédé comporte, une étape de comparaison de la valeur couramment calculée de 7 avec des valeurs contenues dans une table ou un graphe de vieillissement associant des niveaux de vieillissement à des valeurs de 7 , de manière à déduire le niveau de vieillissement de la cellule.