GB2610366A

GB2610366A - "Standardized temperature" -based online estimation method for state of health of battery within wide temperature range

Info

Publication number: GB2610366A
Application number: GB2218907.0A
Authority: GB
Inventors: Wang Limei; Lu Dong; Zhao Xiuliang; Cai Yingfeng; Lian Yubo; Ling Heping; Qiao Sibing; Pan Chaofeng; Chen Long
Original assignee: Jiangsu University
Current assignee: Jiangsu University
Priority date: 2020-05-28
Filing date: 2021-01-05
Publication date: 2023-03-01
Anticipated expiration: 2041-01-05
Also published as: WO2021238228A1; CN111693881B; CN111693881A; US20220146591A1; GB202218907D0; GB2610366B; US11442113B2

Abstract

A "standardized temperature"-based online estimation method for the state of health of a battery within a wide temperature range. The method mainly comprises solving a battery IC curve, establishing a quantitative relationship between temperature sensitive feature points and temperatures of a standard battery, performing standardized transformation on the IC curve at different temperatures, and establishing a BOX-COX-transformation-based relationship between capacity sensitive feature points and capacities, wherein the quantitative relationship between the temperature sensitive feature points and the temperatures of the standard battery is mainly obtained by calculating the difference between feature point voltage values at other temperatures and a feature point voltage value at a standard temperature, and obtaining a correspondence between temperatures and feature point voltage offsets by using an Arrhenius function; the relationship between the capacity sensitive feature points and the capacities is established on the basis of BOX-COX transformation; parameter λ in formula (I) is calculated by using a maximum likelihood function, and a transformed feature point height y is thus obtained; and linear fitting is then performed on the SOH of the battery and the transformed feature point height y to obtain a fitted curve, so as to determine the SOH of the battery.

Description

"STANDARDIZED TEMPERATURE"-BASED ONLINE ESTIMATION METHOD FOR STATE OF HEALTH OF BATTERY WITHIN WIDE TEMPERATURE RANGE

TECHNICAL FIELD

The present invention belongs to the electric vehicle technology, specifically relates to the State of Health estimation of a battery.

BACKGROUND

Accurate State of Health estimation is of great significance for improving battery safety and extending battery life. However, due to the complicated working conditions and various influencing factors in the actual conditions, it is difficult to estimate the State of Health accurately and effectively.

State of Health is generally estimated by a maximum available capacity method or an internal resistance method. At present, battery State of Health is mostly obtained by establishing models, such as an electrochemical model, an empirical model and an equivalent circuit model. Though the calculation results of the electrochemical model are extremely accurate, the model parameters of the electrochemical model are too numerous to obtain and the calculation process is very complicated. The empirical model requires a large amount of test data and takes a long time, resulting in a poor versatility. The equivalent circuit model equates a battery by some resistors and capacitors with a series-parallel sti-ucture. Although the State of Health is monitored in real time, the accuracy is not high, and multiple optimization algorithms are usually required.

An incremental capacity curve is obtained by deriving a charging curve of a battery. The incremental capacity curve directly reflects the multi-stage lithium ions intercalation process when the battery is charged. It is a non-destructive method to detect the internal reaction of the battery. By studying the relationship between the inflection points of the Incremental Capacity curve and the State of Health at a standard temperature, a State of Health estimation method by the Incremental Capacity curve is proposed. However, the internal chemical reaction of a battery is directly affected by the temperature, so many basic battery parameters, such as the capacity and internal resistance, will change with the temperature. Meanwhile, the Incremental Capacity curve will also shift with the temperature. The accuracy of the State of Health estimation method by the Incremental Capacity curve varies greatly at different temperatures. Therefore, a State of Health estimation method in a wide temperature range needs to be explored.

SUMMARY

To solve the above problems, based on "standardized temperature", an online State of Health estimation method in a wide temperature range is proposed. It mainly includes calculating Incremental Capacity curve of a battery, establishing a quantitative relationship between the voltage shift of the temperature-sensitive feature point and the temperature of a standard battery, standardized transformation of Incremental Capacity curves at different temperatures and establishing a quantitative relationship between the transformed height of the capacity-sensitive feature point and the State of Health of the battery based on a BOX-COX transformation.

The mentioned Incremental Capacity curve is often calculated by a conventional numerical differentiation method, and by a polynomial fitting and differentiation method. It can also be calculated by referring to the invention patent CN 109632138 A. The present invention only provides examples, the specific calculated method is not limited.

The method of establishing the quantitative relationship between the voltage shift of the temperature-sensitive feature point and the temperature of a standard battery is as follows: The internal resistance of a battery increases as the temperature decreases, resulting in the increase of kinetic loss. The increase of internal resistance is mainly manifested as the right shift of the Incremental Capacity curve and the voltage shift of inflection points in the Incremental capacity curve, as shown in FIG. 1. Therefore, the second inflection point in the Incremental Capacity curve is selected as the temperature-sensitive feature point to establish the quantitative relationship between the voltage shift of the temperature-sensitive feature point and the temperature. Firstly, two standard batteries, new batteries, are selected to stand for 2 hours at each fixed ambient temperature (-5°C, 0°C, 5°C, 10°C, 15°C, 20°C, 25°C, 30°C, 35°C, 40°C, 45°C, 50°C, 55°C) to ensure the consistent temperature inside and outside the standard batteries. Then, the standard batteries are charged with 0.1 C, and the Incremental Capacity curves at different temperatures are extracted from the charging curves. Subsequently, the voltages of the temperature-sensitive feature point at each temperature are recorded. To describe the quantitative relationship between voltage shift of the temperature-sensitive feature and the temperature, taking the standard temperature (25°C) as a reference, the quantitative relationship between the voltage shift of the temperature-sensitive feature point and the temperature is obtained by subtracting voltage of the temperature-sensitive feature point at the standard temperature from the voltage of the temperature-sensitive feature point at other temperatures. And an Arrhenius fitting function is utilized to fit the quantitative relationship: b, y = a expk-) + (1) Where a, b, c are the fitting parameters, T is the temperature, and y represents the voltage shift of temperature-sensitive feature point.

The above standardized transformation of Incremental Capacity curves at different temperatures is as follows: The Incremental Capacity curves at different temperatures are standardized to ensure the accurate estimation of State of Health in a wide range temperature. According to the quantitative relationship between the voltage shift of the temperature-sensitive feature point and the temperature by Arrhenius fitting function, a charging 0-V curve at each high temperature (upper 30 °C) is shifted with the corresponding voltage shift (bringing the temperature into the Arrhenius fitting equation) to get a normalized charging Q-17 curve. Then, a normalized Incremental Capacity curve by standardized temperature is obtained through the conventional numerical differentiation method.

Based on a BOX-COX transformation, the quantitative relationship between the transformed height of the capacity-sensitive feature point and the State of Health of the battery is established as follows: As the battery ages, battery's cathode and anode active materials and recyclable lithium ions are gradually lost. They are mainly reflected in the height decrease of the Incremental Capacity curve, as shown in FIG. 1. Therefore, the height of the second inflection point is used as the capacity-sensitive feature point.

In this method, the Box-Cox transformation is used to increase the linearity between the height of the capacity-sensitive feature point and the State of Health. A linear regression equation is expressed as: = X +6 (2) Where, Y is the dependent variable, X is the independent variable, ,a is the coefficient matrix, and E represents the fitting error.

The Box-Cox transformation is expressed as: (3) Where, y in the right side of the equation is the original variable, and the subscript k corresponding toy represents the k-th variable. A is the conversion parameter. Yk in the left side of the equation is the k-th conversion variable.

According to the above formula, the inverse transformation of)! is: (I + 2.y(2))1/2 2 # 0 exp(y,(A1) A = 0 (4) A maximum likelihood function is used to calculate the optimal A. Assuming that E is independent and obeys normal distribution, and y conforms to 31-(X16,a21) , X is the independent variable matrix, ig is the coefficient matrix, C2 is the variance, I is the identity matrix, it is the number of samples. The maximum likelihood function is expressed as: L(6, , X) = -L72.1n2n--721n 0-2 --1 62 EY -x,e1T -x 131+ (A -1): 2 k 1 = (5) In the above formula (5), Yin is the dependent variable after conversion The estimated parameters and 02 can be expressed as: /1(2)= (XTX) 11{TY2 02(2) = -1 [Y(2) -XfliT [PA) -XII] Substituting formula (6) into (5) and taking logarithmic: L(2) = -11(ln 27r + ln c2 + + (2 - in y, 2 k =1 The optimal A is obtained by maximizing formula (7).

According to the linear fitting formula (2), the quantitative relationship between transformed height of capacity-sensitive feature point and the State of Health is established. The benefits of the invention: I. According to standardized temperature, the temperature range of the State of Health estimation method by the Incremental Capacity curve is effectively broadened. The accuracy of the State of Health estimation by the Incremental Capacity curve in a wide temperature range is also improved by "standardized temperature" transformation.

2. The Box-Cox transformation is introduced to reduce the influence of random errors on the estimation result, increase the linearity between the height of the capacity-sensitive feature point and the State of Health, and improves the stability of the State of Health estimation by the Incremental Capacity curve.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG shows the Incremental Capacity curves at different State of Health and temperature FIG 2 gives the procedure of the proposed State of Health estimation method FIG 3 shows the fitting curve of the voltage shifts with temperatures FIG. 4 -FIG 9 shows the normalized Incremental Capacity curve at 30°C, 35°C, 40°C, 45°C, 50°C and 55°C respectively.

FIG. 10 shows the fitting results of State of Health and the transformed height of the (6) (7) capacity-sensitive feature point.

DETAILED DESCRIPTION OF THE EMBODIMENTS

It will be further explained below about the figures.

FIG. 1 gives the Incremental Capacity curves at different State of Health and temperature. When the temperature decreases, the Incremental Capacity curve has a right shift, and the height of inflection point decreases as the battery ages. But the height of the inflection point is also affected by the voltage, so the voltage shift caused by temperature is excluded by temperature standardized transformation.

FIG. 2 gives the procedure of the proposed State of Health estimation in the invention. It mainly includes two parts, an offline calibration stage and an online estimation stage The State of Health of the selected five batteries are shown in Table. I.

Table 1

Battery No. 1 No. 2 No. 3 No. 4 No. 5 State of Health 0.98 0.88 0.82 0.70 0.56 In the offl ne calibration stage, the standard batteries, new batteries, are stand for 2 hours at each fixed ambient temperature (-5°C, 0°C, 5°C, 10°C, 15°C, 20°C, 25°C, 30°C, 35°C, 40° C, 45°C, 50°C, 55°C) to ensure the consistent temperature inside and outside the standard batteries. Then, the standard batteries are charged with 0.1 C, and the Incremental Capacity curves at different temperatures are extracted from the charging curves by numerical calculation. The quantitative relationship between the voltage shift of the temperature-sensitive feature point and the temperature is obtained by subtracting the voltage of the temperature-sensitive feature point at the standard temperature from the voltage of the temperature-sensitive feature point at other temperatures. Then, the Arrhenius fitting function is utilized to fit the quantitative relationship. The fitting results are shown in FIG. 3, and the parameters a = 3,66328E-]2, b = 6300.19841, c = -0.00621.

According to the fitting results, the charging curve at each high temperature (upper 30 °C) is shifted with the corresponding voltage shift to get the normalized charging 0-V curve. Then, the normalized Incremental Capacity curve at standardized temperature is obtained by the conventional numerical differentiation method, by the polynomial fitting and differentiation method or by referring to the invention patent CN 109632138 A. FIG. 4 -FIG. 9 are shows the normalized Incremental Capacity curves under 30°C, 35°C, 40°C, 45°C, 50°C, 55°C, respectively.

After that, the height of the capacity-sensitive feature point is extracted and transformed by the Box-Cox method.

The Box-Cox transformation is expressed as: (A)

-2=0 (8)

Where, y in the right side of the equation is the original variable, and the subscript k corresponding toy represents the k-th variable. A is the calculated conversion parameter. in the left side of the equation is the k-th conversion variable.

According to the above formula, the inverse transformation of y is: f(1 +23,,,(11" A#0 (9) The maximum likelihood function is used to calculate the optimal A. Assuming that E is independent and obeys normal distribution, and y conforms top: (Xle'crf) , X is the independent variable matrix, P is the coefficient matrix, 02 is the variance, identity matrix, the maximum likelihood function is expressed as: is the L(13,o-2Y,X) n In 27r 2 2 -lexp(y,(t)) A= 0 1 Y") Xfl Y(A) -Xfl +(A-1)1 lnyfr (10) --CT k=1 Where, ti is the total number of samples. The estimated parameters le and c can be expressed as: ifl(A)= (X' X)-' XT Y2 a2 (A)= -1 Y") -xfir yuo-xfre 77 (11) Substituting formula (11) into (10) and taking logarithmic: 1,(2)= 2g +In a' -F1)+(A -I)Llny, k-1 (12) Taking its maximum value to get A = 3.3, and bringing it into formula (8) to get the transformed height of the capacity-sensitive feature point The relationship between the transformed height of the capacity-sensitive feature point and the State of Health is established, as shown in FIG. 10. The coefficient matrix P = [0.52055, 1.24712E-14]T.

In the online estimation stage, the charging voltage, current and temperature data of a measured battery are firstly recorded when the measured battery is charged, and then the Incremental Capacity curve is obtained. When the voltage of the temperature-sensitive feature point is reached, the normalized Incremental Capacity curve by standardized temperature is obtained according to the Arrhenius fitting function calibrated in the offline calibration stage. Then the Box-Cox transformation is performed on the height of the capacity-sensitive feature point of the measured battery, and the conversion parameter A has been calculated by the offline calibration stage. Finally, the transformed data is brought into the offline calibrated linear regression equation to estimate the State of Health of the measured battery. The details above are only specific descriptions of the feasible implementations of the present invention. They are not used to limit the scope of protection of the present invention. Any equivalent methods or modifications created by the above technology should be included in the protection scope of the present invention.

Claims

CLAIMSWhat is claimed is: I. An on-line State of Health estimation method for a battery in a wide temperature range based on "standardized temperature", characterized by comprising the following steps: SI, calculating the Incremental Capacity curve of a battery; S2. establishing a quantitative relationship between a voltage shift of the temperature-sensitive feature point and the temperature of a standard battery; S3. standardized transformation of Incremental Capacity curves at different temperatures; S4. establishing a quantitative relationship between a transformed height of a capacity-sensitive feature point and a State of Health of the battery based on a BOX-COX transformation.
2. The on-line State of Health estimation method of the battery in the wide temperature range based on "standardized temperature" according to claim 1, characterized in that the implementation of S2 comprises the selection of the temperature-sensitive feature point: internal resistance of the battery increases as the temperature decreases, resulting in the increase of kinetic loss; the increase of internal resistance is mainly manifested as the right shift of the Incremental Capacity curve and the voltage shift of the inflection points in the Incremental capacity curve; therefore, the second inflection point in the Incremental Capacity curve is selected as the temperature-sensitive feature point to establish the quantitative relationship between the voltage shift of the temperature-sensitive feature point and the temperature.
3. The on-line State of Health estimation method of the battery in the wide temperature range based on "standardized temperature" according to claim 2, characterized in that based on the temperature-sensitive feature point selected in claim 2, the specifics of establishing the quantitative relationship between the voltage shift of the temperature-sensitive feature point and the temperature in S2 is as followed: firstly, two standard batteries, new batteries, are selected to stand for 2 hours at each fixed ambient temperature (5°C, 0°C, 5°C, 10°C, 15°C, 20°C, 25°C, 30°C, 35°C, 40°C, 45°C, 50°C, 55°C) to ensure the consistent temperature inside and outside the standard batteries; the standard batteries are then charged with 0.1 C, and the Incremental Capacity curves at different temperatures are extracted from the charging curves; subsequently, the voltages of the temperature-sensitive feature point at each temperature are recorded; to quantitatively describe the relationship between voltage shift of the temperature-sensitive feature point and the temperature, taking the standard temperature (25 ° C) as a reference, the quantitative relationship between the voltage shift of the temperature-sensitive feature point and the temperature is obtained by subtracting the voltage of the temperature-sensitive feature point at the standard temperature from the voltage of the temperature-sensitive feature point at other temperatures; finally, an Arrhenius fitting function is utilized to fit the quantitative relationship: y=aexp(-)+c (1) wherein a, h, c are the fitting parameters, Tis the temperature, andy represents the voltage shift of temperature-sensitive feature point.
4. The on-line State of Health estimation method of the battery in the wide temperature range based on "standardized temperature" according to claim 3, characterized in that standardized transformation of the Incremental Capacity curves at different temperatures in S3 is as followed: according to the quantitative relationship between the voltage shift of the temperature-sensitive feature point and the temperature by the Arrhenius fitting function, a charging O-V curve at each high temperature is shifted with a corresponding voltage shift to get a normalized charging 0-V curve; then, a normalized incremental Capacity curve by the standardized temperature is obtained through a numerical differentiation method.
5. The on-line State of Health estimation method of the battery in the wide temperature range based on "standardized temperature" according to claim 4, characterized in that the charging O-V curve at high temperature refers to charging O-V curves with the temperature greater than 30 ° C; the corresponding voltage shift of the charging 047 curve is achieved by bringing the temperature into the Arrhenius fitting equation.
6. The on-line State of Health estimation method of the battery in the wide temperature range based on "standardized temperature" according to claim 1, characterized in that the selection of the capacity-sensitive feature point in S4 is as followed: as the battery ages, battery's cathode and anode active materials and recyclable lithium ions are gradually lost, it is mainly reflected in the height decrease of the Incremental Capacity curve, therefore, the height of the second inflection point is used as the capacity-sensitive feature point.
7. The on-line State of Health estimation method of the battery in the wide temperature range based on "standardized temperature" according to claim 6, characterized in that based on the capacity-sensitive feature point, the implementation of step 4 comprises: the Box-Cox transformation is used to increase the linearity between the height of the capacity-sensitive feature point and the State of Health; a linear regression equation is expressed as: Y = X if+ c (2) wherein Yis a dependent variable, Xis an independent variable, /I is a coefficient matrix, and E represents a fitting error; Box-Cox transformation is represented as: YA2--1 2#0AIn y, A=0 (3) wherein y in the right side of the equation is an original variable, and the subscript k y(it) corresponding toy represents a k-th variable, A is a conversion parameter, ./k in the left side of the equation is a k-th conversion variable.
8. The on-line State of Health estimation method of the battery in the wide temperature range based on "standardized temperature" according to claim 7, characterized in that the implementation of step 4 also includes: an inverse transformation of y is: (1+ Ay:2)u A, # 0 exp(y;21) 2=0 (4) a maximum likelihood function is used to calculate an optimal A assuming that E is = independent and obeys a normal distribution, and y conforms to Y (X Acr), X is an independent variable matrix, fl is the coefficient matrix, ci is a variance, I is an identity matrix, 71 is a number of samples, the maximum likelihood function is expressed as: L(13,0-211' , X) = in 27r -111 In o-2 2 2 1 -X f31-F (A -1)E In y, (5) Y'A'-xpl in the formula (5), is a dependent variable after conversion, estimated parameters fl and ci can be expressed as: ma,)=(xlx)'x' cr2 (2) Ly( A) x /if Ly(A) xpi substituting formula (6) into (5) and taking logarithmic: L(2) = -1172 (ln 27-c In c2 + 1) + (2 -1)± In y k 1 optimal A is obtained by maximizing formula (7); according to the linear regression equation (2), the quantitative relationship between the transformed height of the capacity-sensitive feature point and the State of Health is then established. (6) (7)