FR 3049352 A1 20170929 FR 1652537 A 20160324 PROCEDE DE DETERMI NATI ON DE L'ETAT DE SANTE D'UNE BATTERI E L'invention se rapporte à un procédé d'estimation de l'état de santé d'une batterie. L'algorithme de calcul du SOHE (de l'anglais State of Health Energy) ou de l'état de santé énergétique d'une batterie est connu en soi. Le SOHE est utilisé pour déterminer le niveau de vieillissement d'une batterie et une fois que ce paramètre atteint une valeur seuil, la batterie est jugée trop vieille et peut être retirée du circuit d'utilisation. Le SOHE est également utilisé pour d'autres calculs internes au BMS (de l'anglais Battery Management System), comme l'énergie disponible par exemple. Pour rappel, le BMS est un système électronique permettant de gérer et de contrôler les différents paramètres d'une batterie, comme par exemple son état de charge ou SOC (de l'anglais State Of Charge). Les calculs du SOHE sont habituellement assurés par le BMS lors de l'utilisation de la batterie. Or, il arrive bien souvent que ces calculs ne soient pas très fiables, avec parfois des écarts observés par rapport à la véritable valeur du SOHE, pouvant atteindre jusqu'à 10%, voire 20% de la valeur du SOHE. Afin de connaitre avec précision la valeur du SOHE, il est important de se donner les moyens de pouvoir sélectionner les valeurs du SOHE calculées, qui soient les plus proches possibles de la véritable valeur de ce paramètre, et ainsi d'éliminer les valeurs qui en sont éloignées. Il est également possible d'envisager de pondérer différemment les valeurs du SOHE retenues en fonction de l'erreur par rapport à la vraie valeur, mais il demeure cependant difficile de déterminer avec rigueur et précision, l'erreur de la valeur du SOHE calculée par le BMS par rapport à la vraie valeur. Une possibilité toutefois de connaitre cette erreur, consiste à effectuer un grand nombre d'expérimentations, puis de construire une table et/ou une 30 équation. Les procédés de détermination par calcul de l'état de santé d'une batterie existent et ont déjà fait l'objet de brevet. On peut par exemple citer la 3049352 2 demande US2015120225, dont le procédé présente les caractéristiques suivantes : - une première valeur de l'état de santé SOH1 (de l'anglais State Of Health) est calculée si la variation de l'état de charge SOC est 5 supérieure à un seuil. C'est un calcul classique du type soin. = f idt Capactt laie M'OC. Si SOH1 est supérieure à un seuil, alors une deuxième valeur de l'état de santé SOH2 est calculée, via - un calcul de capacité utilisant des couplets (Jidt, ASOC) mesurés (où i 10 est l'intensité du courant traversant la batterie). Il faut plusieurs couplets, mais rien n'est précisé concernant le nombre. Ce SOH2 est en fait un rapport entre la capacité courante et la capacité initiale, la capacité courante étant calculée par la méthode des moindres carrés avec les couplets (Jidt, ASOC). 15 Le procédé décrit dans ce document, présente l'inconvénient d'engendrer des erreurs de calcul, si le temps est grand, à cause du terme Jidt. De plus, il est difficile, voire impossible d'évaluer l'erreur commise sur la capacité calculée, et donc sur SOH2. Un procédé selon l'invention permet de déterminer de façon plus 20 rigoureuse, fiable et plus précise l'état de santé d'une batterie en introduisant notamment des facteurs correctifs permettant d'affiner la valeur du SOHE. L'invention a pour objet un procédé de détermination de l'état de santé énergétique d'une batterie. La principale caractéristique d'un procédé selon l'invention est qu'il 25 comprend les étapes suivantes, - une étape de calcul dudit état de santé à partir d'au moins un paramètre mesuré sur la batterie, - une étape de calcul de l'imprécision E réalisée sur le calcul dudit état de santé et résultant d'erreurs de mesure sur ledit au moins un paramètre, 3049352 3 - une étape de consolidation dudit état de santé consistant à incorporer le calcul de l'imprécision E au calcul de l'état de santé de la batterie. Le principe d'un tel procédé est d'introduire un facteur correctif 5 permettant de consolider la valeur calculée de l'état de santé de la batterie au moyen d'équations conventionnelles, ledit facteur correctif prenant en compte une imprécision E due notamment aux erreurs de mesures sur les paramètres de la batterie. Le terme « imprécision » est assimilable au terme « erreur », et la présence d'une telle imprécision a pour conséquence néfaste d'engendrer un 10 écart entre la valeur du SOHE obtenue par le calcul et la véritable valeur de celui-ci. De façon préférentielle, ledit au moins un paramètre est à choisir parmi une tension aux bornes de la batterie et un état de charge. Avantageusement, l'étape de calcul de l'état de santé est réalisée à 15 partir de l'équation, 1 Iftt I (T) d-cl SOHE =Qo x 100 If (OCV2)-f(OCV1)1 L'imprécision E résultant de ce calcul. Les différents paramètres intervenant dans cette équation sont les suivants, 20 Qo est la capacité à l'état neuf de la batterie OCV1 et OCV2 sont les tensions aux bornes de la batterie en circuit ouvert à l'instant t1 de début de la mesure et à l'instant t2 de fin de la mesure respectivement, 1(t) est l'intensité du courant traversant la batterie, f(OCV2) - SOC2 est 25 l'état de charge à l'instant t2, et f(OCV1) = SOC1 est l'état de charge à l'instant t1. De manière idéale, les SOC1 et SOC2 peuvent correspondre à des états relaxés de la batterie. Dans ces cas, le SOC peut être directement calculé à 3049352 4 partir de la tension cellule ou batterie et de la caractéristique SOC-OCV (SOC=f(OCV)). Dans l'hypothèse où il n'est pas possible d'avoir des cellules relaxées, pour lesquelles la détermination directe du SOC par la mesure de la tension 5 cellule n'est pas adaptée, le SOC peut être déterminé à l'aide d'estimateur, par exemple un filtrage de Kalman. Cependant, pour le calcul de l'erreur (ou d'imprécision) du SOHE, l'hypothèse d'un SOC déterminé à partir d'une mesure de tension relaxée est considérée en particulier. 10 L'imprécision sur le SOC due à d'autres causes, comme l'erreur de SOC calculé par un estimateur par exemple, peut être ajoutée par ailleurs dans le calcul d'erreur. De façon préférentielle, l'étape de consolidation comprend, - une étape de calcul des imprécisions Ei faites durant les calculs de 15 l'état de santé à M instants successifs ti respectivement, où 1M et où M est un entier, - une étape de mémorisation desdites imprécisions Ei, - une étape de calcul de l'état de santé de la batterie à partir de l'équation, M E SOHEagrégé = (1 i SOHEi i=1 Ej'L El 20 où SOHE; désigne l'état de santé à l'instant ti. Connaissant l'erreur faite sur le SOHE, il est possible alors de l'utiliser pour produire une valeur de SOHE agrégeant une ou plusieurs valeurs successives de SOHE avec un niveau de précision connu. Un exemple illustré par l'équation ci-dessus, consiste à pondérer 25 chaque valeur de SOHE calculée avec un poids établi à partir de l'erreur E associée, afin de calculer un SOHE agrégé. Préférentiellement, l'imprécision de calcul E résulte, 3049352 5 - d'une erreur sur la mesure de l'état de charge de la batterie, estimée par l'équation suivante : Af(OCVi) = -dSOC I x (a + dOCVocv, fi) dSOC Où - correspond à l'erreur en pourcentage d'état de charge par dOCV10CV/ unité de tension au point de tension en circuit ouvert à l'instant t; ; 5 où a représente l'erreur de mesure de la tension, et peut être connue soit expérimentalement, soit en consultant la fiche technique du constructeur du système de mesure (figure 1); où fl représente l'écart entre la tension relaxée et la tension au moment de la mesure à l'instant t,, 10 L'expression générale du facteur multiplicatif est (a+,0), mais ledit facteur peut aussi être égal, soit à a, ou soit à p. d'une erreur sur la mesure sur la mesure du courant traversant la batterie, estimée par l'équation, t,N 3,3=ft yd-rç'r'sY=NTsY i=1 Où Ts est la période d'échantillonnage du signal mesuré, N le nombre 15 de points échantillonnés utilisés pour l'intégration et y le biais du capteur de courant. De façon avantageuse, un procédé selon l'invention comprend une étape de calcul des imprécisions Ei faites durant les calculs de l'état de santé aux instants successifs ti au moyen de l'équation suivante, 20 A(SOHE)13^3 E = SOHE I = Af(Ocv1) f (ocv2) - f (ocvi) Af(Ocv2) f (ocv2) - f (ocvi) Où 3 est l'intégrale du courant mesuré entre les instants tl et t2. Cette équation, qui fait intervenir la somme de trois termes est la plus complète et 3049352 6 la plus précise. Mais elle ne pourrait faire intervenir qu'un seul des trois termes au choix, ou chacune des trois combinaisons possibles d'addition de deux desdits termes. Selon un autre mode de réalisation préféré d'un procédé selon 5 l'invention, ledit procédé comprend une étape de calcul des imprécisions Ei faites durant les calculs de l'état de santé aux instants successifs ti au moyen de l'équation suivante, E' = A(SOHE) SOHE (OCV1) 2 2 ( 3 ) + of - f(OCV1)) + g (OCV2 ) 2 f(OCV2) - f(OCV1) \ i ) i où 3 est l'intégrale du courant mesuré entre les instants t1 et t2. Cette équation, qui fait intervenir la somme de trois termes est la plus complète et 10 la plus précise. Mais elle ne pourrait faire intervenir qu'un seul des trois termes au choix, ou chacune des trois combinaisons possibles d'addition de deux desdits termes. Un procédé selon l'invention présente l'avantage de pouvoir déterminer un état de santé de la batterie qui soit parfaitement exploitable, 15 dans le sens où la valeur de cet état de santé est très voisine de la véritable valeur de cet état de santé. Un tel procédé apparait incontournable lors de la détection d'un état de santé d'une batterie destinée à être embarquée dans un véhicule, car la connaissance précise de cet état va notamment conditionner certains paramètres de roulage dudit véhicule, tels que par exemple, son 20 autonomie. On donne ci-après une description détaillée d'un procédé de détermination de l'état de santé d'une batterie selon l'invention, en se référant aux figures suivantes, - La figure 1 est un diagramme illustrant un exemple de distribution 25 de la mesure de tension, - La figure 2 est un diagramme illustrant un exemple d'un écart entre l'instant de mesure et la valeur après relaxation des cellules de la batterie, 3049352 7 - La figure 3 est un diagramme illustrant un exemple de variation dSOC/dOCV en fonction de l'état de charge SOC d'une cellule de la batterie, - La figure 4 est un logigramme simplifié illustrant les différentes 5 étapes du calcul juste de l'état de santé d'une batterie à partir d'un procédé selon l'invention. L'algorithme de calcul du SOHE (de l'anglais State of Health Energy ou état de santé énergie) d'une batterie est classique. Ce SOHE est utilisé pour déterminer le niveau de vieillissement d'une batterie. Une fois un seuil atteint, 10 la batterie peut être retirée du circuit d'utilisation. Ce SOHE est également utilisé pour d'autres calculs internes au BMS (de l'anglais Battery Management System), comme l'énergie disponible par exemple. Ce calcul du SOHE est d'ordinaire calculé par le BMS lors de l'utilisation 15 de la batterie. Plusieurs valeurs de SOHE sont alors calculées et peuvent être différentes pour des raisons qui seront indiquées dans la suite de la description. Le problème ici considéré est de pouvoir sélectionner les valeurs de SOHE qui soient le plus proches possibles de la valeur vraie et d'éliminer les 20 valeurs éloignées. Il est possible aussi d'envisager de pondérer différemment les valeurs de SOHE retenues en fonction de l'erreur par rapport à la vraie valeur. Il est, cependant, difficile de déterminer l'erreur de la valeur de SOHE calculée par le BMS par rapport à la vraie valeur. 25 Une possibilité est d'estimer cette erreur sur un grand nombre d'expérimentations, puis de construire une table ou une équation (modèle) L'idée ici est de calculer cette erreur par une équation simple fondée sur la théorie. De manière classique, le SOHE est défini par l'expression suivante : SOHE A -E0x100 [1] 3049352 8 Où E est l'énergie en [Wh] lorsque la batterie est complètement chargée, E0 est l'énergie en [Wh] lorsque la batterie est complètement chargée et à l'état neuf (début de vie). Cette définition est établie pour une température de 25°C. 5 L'inconvénient de cette définition est que l'énergie est difficilement calculable directement, à la différence de l'état de charge (SOC) qui peut être déterminé directement par la mesure de la tension cellule avec la caractéristique SOC-OCV (de l'anglais Open Circuit Voltage) fournie par le constructeur de cellule. 10 C'est pourquoi, dans la pratique, on calcule plutôt le SOHE de la manière suivante : SOHE A -Q0 x 100 Q Où Q est la capacité en [Ah] de la batterie, Q0 est la capacité à l'état neuf ou BOL (de l'anglais Begin Of Life). Cette définition est établie à une 15 température de 25°C. La capacité Q peut être calculée en ligne par le BMS par la relation suivante : If» 1(1 d-cl Q ,i ''` I SOC2 - SOCi I 3] Où 1(t) est le courant de la batterie, t1 est l'instant du début de la mesure, t2 est l'instant de la fin de la mesure, SOC1 est l'état de charge à 20 l'instant t1, SOC2 est l'état de charge à l'instant t2. De manière idéale, les SOC1 et SOC2 peuvent correspondre à des états relaxés de la batterie. Dans ces cas, le SOC peut être directement calculé à partir de la tension cellule ou batterie et de la caractéristique SOC-OCV (SOC=f(OCV)). 25 Dans ce cas, Q peut s'écrire tel que : Q 1 Iftt21(t) cl-c `'` If (OCV2) - f (OCVI) 2] [ 4] D'où : 3049352 9 1 Iftt: I(r) d-cl SOHE = Q0 I x 100 f (OCV2) - f (OCVID I Dans l'hypothèse où il n'est pas possible d'avoir des cellules relaxées, pour lesquelles la détermination directe du SOC par la mesure de la tension cellule n'est pas adaptée, le SOC peut être déterminé à l'aide d'estimateur, par exemple un filtrage de Kalman. 5 Cependant, pour le calcul de l'erreur (ou d'imprécision) de SOHE, l'hypothèse d'un SOC déterminé à partir d'une mesure de tension relaxée est considérée en particulier. L'erreur sur le SOC due à d'autres causes, comme l'erreur de SOC calculé par un estimateur par exemple, pourra être ajoutée par ailleurs dans le 10 calcul d'erreur. L'expression de l'erreur sur SOHE est obtenue de la manière suivante à partir de l'équation [5] 1 QN rSOHE a (soocHE21 d(SOHE) = c13-F Hafa(So0cvEi)d(f(OCV1))az af v) ci(f(ocv2)) Avec : .3- ft2 l(r) cit ti 15 Puis, d(SOHE)1 -1 1 f8] SOHE = cL3+f (ocv2) - f (ocvi)d(f (OCV1)) + f (ocv2) - f (ocvi) d(f (OCV2) ) De cette expression, on peut tirer 2 expressions pour l'erreur de SOHE. Une première expression de l'erreurs, dite erreur « pire cas » : A(SOHE)13,3 Af. (Ocv1) f (0 cv2) - f (004) Af. (0 cv2) f (0 cv2) - f (004) [9] £ SOHE I - + Une autre expression est également possible, dite erreur « norme (0CV1) )2 + ( (0CV2) 2 P (312 f (ocv2) - f (004) (0 0/ 2) - f (ocvi) A partir des équations [9] ou [10], il est possible de constater que : 20 1. les erreurs relatives s'accumulent, [5] [6] [7] A(SOHE) = SOHE 3049352 10 2. plus l'écart entre SOC1=f(OCV1) et SOC2=f(OCV2) devient grand, plus l'impact de l'erreur Af(OCV1) et de l'erreur Af(OCV2) surs devient faible. Causes d'erreur Af(OCV,) et calcul de Af(OCV,1 5 L'erreur Af(OCV,) correspond à l'erreur ASOC, qui est l'erreur sur le SOC au point i. Comme on considère ici le cas dans lequel le SOC est déterminé à partir de la mesure de tension, deux causes d'erreur distinctes sont combinées dans la caractéristiques SOC-OCV de la cellule : 10 - l'imprécision a de la mesure de la tension cellule, - l'écart f3 de tension entre l'instant de mesure et la valeur après relaxation. En se référant à la figure 1, si l'imprécision a de la mesure de la tension cellule a une distribution selon une loi normale, a peut alors être prise 15 comme l'erreur de mesure à ±xcy, où x=1 ou 2 ou 3 par exemple, 6 étant l'écart-type. La figure 2 illustre l'écart f3 de tension entre l'instant de mesure et la valeur après relaxation. Le courant est coupé ou devient nul à l'instant to après une décharge par exemple, la tension se relaxe pour atteindre sa tension 20 relaxée ou OCV (open circuit voltage). Si la mesure de tension intervient à l'instant tO+At, alors un écart f3 par rapport à l'OCV est introduit tel que : V = OCV + 8fl Où V est la tension mesurée à l'instant tO+At, OCV est la tension relaxée, f3 est l'écart en tension, 6 est défini tel que : +1 si le courant à to est un courant de charge [12] = -1 si le courant à to est un courant de décharge Ces erreurs a et f3 peuvent s'ajouter et se combinent ensuite dans la 25 caractéristique SOC-OCV de la cellule considérée. Il peut être utile de calculer la courbe (dSOC/dOCV)-OCV qui correspond à la dérivée de la courbe SOC-OCV pour tous les points d'OCV. Cette courbe en %/mV correspond à l'erreur en % de SOC par mV d'erreur sur la mesure d'OCV. 3049352 11 La figure 3 fournit un exemple de courbe (dSOC/dOCV)-OCV pour une cellule Li-ion. L'erreur en SOC ASOC,=3.f(OCV,) est calculée tel que : dSOC Af. (OCVi) = - x + /3)[ 13] dOCVlocv, Ce qui correspond à l'erreur en % de SOC par mV au point d'OCV, et 5 multipliée par l'erreur totale maximale (a+[3). La figure suivante illustre le calcul de l'équation [13]. Causes d'erreur A et calcul de l'erreur 4.3 La principale source d'erreur A3 sur l'intégrales du courant est l'offset 10 ou biais sur la mesure du courant. Selon la technologie du capteur de courant utilisé, cette erreur 3.73 peut avoir plusieurs raisons physiques. Par exemple pour un capteur à effet Hall utilisant un tore magnétique, le biais de mesure peut être classé en deux catégories : le biais électrique et le biais magnétique, les deux se combinant 15 pour former le biais total. L'erreur 3.73 est définie d'après la relation suivante : t, t, t, 3,3= rte + dt = ft1 I cit + f y cit Où y est le biais du capteur de courant et supposé constant au cours d'une mission. 20 Le calcul de 3.73 est approché de la manière suivante : t, AZ's= f y cit Ts y = N Ts y i=1 Où Ts est la période d'échantillonnage du signal mesuré, N le nombre de points échantillonnés utilisés pour l'intégration. Utilisation du calcul d'erreur E Connaissant l'erreur faite sur le SOHE, il est possible alors de l'utiliser 25 pour produire une valeur de SOHE agrégeant une ou plusieurs valeurs successives de SOHE avec un niveau de précision connu. [ 14] N [ 15] 3049352 12 Un exemple consiste à pondérer chaque valeur de SOHE calculée avec un poids établi à partir de l'erreurs associée, afin de calculer un SOHE agrégé donné par : M [1 6] Ei SOHEagrégé = (1 ,m SOHEi i=1 Ej Concrètement, un tableau de M colonnes de 2 lignes, soit Mx2 cases, est provisionné. La aère ligne recueille les valeurs de SOHE calculé comme à l'équation [5], la seconde la valeur d'imprécision associée donnée à l'équation [9] ou [10]. Un exemple est illustré par le Tableau 1 ci-après. SOHE SOHE1 ... SOHE; ... SOHEM imprécision El ... Si - - - EM Tableau 1 - Stockage des SOHE et de leur imprécision associée Ensuite, lorsque le tableau est plein, les valeurs de la première colonne sont retirées, les valeurs restantes sont décalées, puis la dernière colonne est remplie avec la nouvelle valeur de SOHE disponible et son imprécision 15 associée. C'est donc un stockage FIFO (first-in first-out) décrit ici, ce qui permet de faire un calcul de SOHEagrégé sur une fenêtre glissante de largeur M. Il est cependant possible de calculer SOHEagrégé avant que les M colonnes ne soient remplies. Il est possible également de traiter les M cases du tableau d'une manière différente. Par exemple, en retirant les ml cases 20 correspondant à des valeurs de SOHE les plus grandes et/ou les m2 cases correspondant à des valeurs de SOHE les plus faibles, avec m1+m2 REVENDICATIONS1. Procédé de détermination de l'état de santé énergétique d'une batterie, caractérisé en ce qu'il comprend les étapes suivantes, - une étape de calcul dudit état de santé à partir d'au moins un paramètre mesuré sur la batterie, - une étape de calcul de l'imprécision E réalisée sur le calcul dudit état de santé et résultant d'erreurs de mesure sur ledit au moins un paramètre, - une étape de consolidation dudit état de santé consistant à incorporer le calcul de l'imprécision E au calcul de l'état de santé de la batterie. 2. Procédé selon la revendication 1, caractérisé en ce que ledit au moins un paramètre est à choisir parmi une tension aux bornes de la batterie et un état de charge. 3. Procédé selon l'une quelconque des revendications 1 ou 2, caractérisé en ce que l'étape de calcul de l'état de santé est réalisée à partir de l'équation 1 Iftt I (T) d-cl SOHE =Qo x 100 If (OCV2)-f(OCV1)1 l'imprécision E résultant de ce calcul, Où Qo est la capacité à l'état neuf de la batterie, Où OCV1 et OCV2 sont les tensions aux bornes de la batterie en circuit ouvert à l'instant t1 de début de la mesure et à l'instant t2 de fin de la mesure respectivement, Où 1(T) est l'intensité du courant traversant la batterie, f(OCV2) = SOC2 est l'état de charge à l'instant t2, et f(OCV1) = SOC1 est l'état de charge à l'instant t1. 4. Procédé selon l'une quelconque des revendications 1 ou 2, caractérisé en ce que l'étape de consolidation comprend, 3049352 15 - une étape de calcul des imprécisions El faites durant les calculs de l'état de santé à M instants successifs ti respectivement, où 1M et où M est un entier, - une étape de mémorisation desdites imprécisions Ei, 5 - une étape de calcul de l'état de santé de la batterie à partir de l'équation M Ei SOHEagrégé = (1 SOHEi i=1 Lj=i El où SOHE; désigne l'état de santé à l'instant ti. 5. Procédé selon l'une quelconque des revendications 1 ou 2, caractérisé en ce que l'imprécision de calcul E résulte, 10 - d'une erreur sur la mesure de l'état de charge de la batterie estimée par l'équation, x (a + ) dSOC pf(OCVi) = -dOCV1 /j)ocv, dSOCI Où correspond à l'erreur en pourcentage d'état de charge par dOCV I 0CV/ unité de tension au point de tension en circuit ouvert à l'instant t; ; où a représente l'erreur de mesure de la tension ; 15 où /1 représente l'écart entre la tension relaxée et la tension au moment de la mesure à l'instant t,, - d'une erreur sur la mesure du courant traversant la batterie estimée par l'équation, t,N 3,3=ft yd-rçr'sY=NTsY i=1 Où Ts est la période d'échantillonnage du signal mesuré, N le 20 nombre de points échantillonnés utilisés pour l'intégration et y le biais du capteur de courant. 3049352 16 6. Procédé selon les revendications 4 et 5, caractérisé en ce qu'il comprend une étape de calcul des imprécisions Ei faites durant les calculs de l'état de santé aux instants successifs ti au moyen de l'équation suivante, 5 A(SOHE)13, 3 3 Af(OCV1) f (0 cv2) - f (ocvi) + Af(OCV2) f (0 cv2) - f (ocvi) E = SOHE I - 1 + où 3 est l'intégrale du courant mesuré entre les instants t1 et t2 7. Procédé selon les revendications 4 et 5, caractérisé en ce qu'il comprend une étape de calcul des imprécisions Ei faites durant les calculs de l'état de santé aux instants successifs ti au moyen de l'équation suivante, , A(SOHE) - ,\I (12 ( Af(OCV1) ) + 2 ( Af (OCV2) )2 E = SOHE 7:5 + f (OCV 2)- f (OCV 1) f (OCV 2)- f (OCV 1) .1 où 3 est l'intégrale du courant mesuré entre les instants t1 et t2 10