System and Method for Recursively Estimating Battery Cell Total Capacity

Abstract

A system and method for determining an estimated battery cell total capacity indicative of a total capacity of a battery cell is provided. The method includes receiving a first battery cell state-of-charge estimate at a first time and receiving a second battery cell state-of-charge estimate at a second time subsequent to the first time, measuring an integrated battery cell current value indicative of the integrated battery cell current between the first time and the second time, updating at least one recursive parameter based on the first battery cell state-of-charge estimate, the second battery cell state-of-charge estimate, and the integrated battery cell current value, determining the estimated battery cell total capacity based on at least one recursive parameter, and storing a value corresponding to the estimated battery cell total capacity in a memory.

Description

RELATED APPLICATIONS

This applications claims priority to U.S. Provisional Application Ser. No. 61/283,329 filed Dec. 2, 2009, the disclosure of which is hereby incorporated by reference.

FIELD OF THE APPLICATION

This application relates generally to battery cells that can be used in electric vehicles, and more particularly to a method and system for recursively estimating battery cell capacity.

BACKGROUND

Battery cells are used in a wide variety of electronic and electrical devices. In these applications it is often desirable to be able to ascertain the present “health” of a battery cell. One measure of battery cell health is its total capacity, which is the maximum amount of charge that the cell can hold, and is usually measured in ampere hours (Ah) or milliampere hours (mAh). As a battery cell ages, its total capacity tends to decrease. Therefore, it is desirable to be able to estimate the present value of battery cell total capacity.

A number of algorithms have been proposed for the purpose of battery cell total capacity estimation. Very often, there is a tradeoff between the complexity of the algorithm and the operational constraints imposed by the algorithm and the accuracy of the results that it produces. Very accurate results tend to be produced only by very complex algorithms (which require extensive processor resources to execute) or those that impose operational constraints on how the battery is used. Accurate results at a lower complexity and with fewer constraints are desirable.

Accordingly, the inventor herein has recognized a need for an improved method for more efficiently determining an estimated battery cell total capacity. The embodiments disclosed herein impose fewer operational constraints, and requires very little processing to provide accurate estimates of cell total capacity.

BRIEF DESCRIPTION OF THE EMBODIMENTS

A method for determining an estimated battery cell total capacity indicative of a total capacity of a battery cell in accordance with an exemplary embodiment is provided. The method includes receiving a first battery cell state-of-charge estimate at a first time, receiving a second battery cell state-of-charge estimate at a second time subsequent to the first time, and measuring an integrated battery cell current between the first and second times. The method further includes computing updating at least one recursive parameter based on the first battery cell state-or-charge estimate, the second battery cell state-of-charge estimate, and the integrated battery cell current. The method further includes determining an to estimated battery cell total capacity based on at least one recursive parameter and storing a value corresponding to the estimated battery cell total capacity in a memory.

A system for determining an estimated battery cell total capacity indicative of a total capacity of a battery cell in accordance with an exemplary embodiment is provided. The system includes a first subsystem configured to receive a first state-of-charge estimate at a first time, and a second subsystem configured to receive a second state-of-charge estimate at a second time subsequent to the first time. The system further includes a current sensor electrically coupled to the battery cell. The current sensor is configured to generate a signal indicative of an integrated battery cell current between the first and second times. The system further includes a computer comprising the first subsystem, the second subsystem, and a subsystem configured to receive the integrated current signal. The computer is further configured to update at least one recursive parameter based on the first state-of-charge estimate, the second state-of-charge estimate, and the integrated current signal. The computer is further configured to determine the estimated battery cell total capacity based on at least one recursive parameter and to store a value corresponding to the estimated battery cell total capacity in a memory.

A computer readable medium having computer-executable instructions for performing a method for determining an estimated battery cell total capacity indicative of a total capacity of a battery cell in accordance with an exemplary embodiment is provided. The method includes receiving a first state-of-charge estimate at a first time, receiving a second state-of-charge estimate at a second time subsequent to the first time, measuring an integrated battery cell current between the first and second time, updating at least one recursive parameter based on the first state-of-charge estimate, the second state-of-charge estimate, and the integrated battery cell current, and determining the estimated battery cell total capacity based on at least one recursive parameter. The method further includes storing a value corresponding to the estimated battery cell total capacity in a memory.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a schematic of a system for determining an estimated battery cell total capacity and available energy in accordance with an exemplary embodiment;

FIG. 2 illustrates the different objectives behind “least squares” regression and “total least squares” regression;

FIG. 3 illustrates several important distances between a data point and a regression line in the “total least squares” method;

FIG. 4 is a block diagram of an exemplary computational algorithm utilized by the system of FIG. 1; and

FIG. 5 is a result using an exemplary embodiment.

DETAILED DESCRIPTION

The total capacity of a battery cell is a value, usually expressed in ampere hours (Ah) or milliampere hours (mAh), that indicates the maximum electrical charge that the battery cell is capable of holding. New battery cells are manufactured with certain nominal total capacities, but as the cells age, their capacities generally decrease. Therefore, being able to accurately estimate the total capacity of a battery cell is important to being able to determine the health of that battery cell.

The state-of-charge (SOC) of a battery cell is a value between 0% and 100% that indicates the relative level of charge presently held by the battery cell. A state-of-charge of 100% corresponds to a “full” cell, while a state-of-charge of 0% corresponds to an “empty”' cell. State-of-charge is sometimes referred to as “residual capacity” and is not to be confused with the battery cell total capacity. However, the two are related by the equation

$\begin{array}{cc}z\left(t_2\right)=z\left(t_1\right)+\frac{1}{C}{\int }_{t_1}^{t_2}\eta i\left(t\right)t& \left(1\right)\end{array}$

where z(t₂) is the battery cell SOC at time t₂, z(t₁) is the battery cell SOC at time t₁, C is the battery cell total capacity, i(t) is the battery cell current at time t, in amperes, and η is a unitless efficiency factor, which may take on different values depending on whether the current is positive or negative. SOC itself is unitless. Note that in the convention used herein, discharge current is assumed to have a negative sign and charge current is assumed to have a positive sign. The opposite convention may be used with very simple changes to the equations, and the disclosed embodiment applies to that case as well.

Equation (1) is the mathematical basis for most capacity estimating methods. The terms can be rearranged to get:

$\begin{array}{cc}C=\frac{1}{z\left(t_2\right)-z\left(t_1\right)}{\int }_{t_1}^{t_2}\eta i\left(t\right)t\mathrm{or}{\int }_{t_1}^{t_2}\eta i\left(t\right)t=C\left(z\left(t_2\right)-z\left(t_1\right)\right)& \left(2\right)\end{array}$

Using a regression technique, one may compute estimates of C, especially using the second for, where the obvious linear structure of y=Cx becomes apparent. One only needs values for “x” and “y”.

The problem with using standard (least squares) linear regression techniques is that both the integrated current value y and the difference between state-of-charge values x have sensor noise or estimation noise associated with them. The least squares linear regression problem is a solution to the equation (y+Δy)=Cx: that is, there is noise assumed on the measurements y, but not on the independent variable x. However, Equation (2) is implicitly of the form (y+Δy)=C(x+Δx) since both the integrated current and SOC estimates have noise. That is, because estimates of SOC are generally imperfect, there will be noise on the x variable, and using standard least squares linear regression results in an inaccurate and biased estimate of battery cell total capacity.

The usual approach to counteract this problem is to try to ensure that the SOC estimates are as accurate as possible and then use standard least squares estimation anyway. For example, U.S. Pat. No. 6,892,148 puts constraints on how the capacity is estimated. It forces the cell current to be zero before the test begins (so that the cell is in an equilibrium state and the first SOC estimate is accurate) and similarly forces the cell current to be zero after the test ends (again, so that the cell is in an equilibrium state and the second SOC estimate is accurate). This procedure eliminates to a large extent (but not completely) the error in the x variable, and makes the regression reasonably accurate. A second example is U.S. Pat. No. 5,789,924, where battery current is periodically interrupted and a test performed to estimate the cell resistance so that the ongoing SOC estimates may be corrected for ohmic voltage losses. Again, the basic idea is to try to make x as noise-free as possible.

Both of these methods require a very structured and invasive setting for determining battery cell total capacity. The battery cell current must be controlled by the algorithm. Furthermore, neither one correctly handles the residual noise in x (while they minimize the noise, they never totally eliminate it). What is needed is a method for determining battery cell total capacity in a non-invasive setting where the method does not impose constraints on the battery cell current, and where noise in the SOC estimates is correctly accounted for. The system and method disclosed herein may be used to determine an accurate estimated battery cell total capacity in a computationally efficient manner, without imposing constraints on the battery cell current.

An overview of the method for estimating a battery cell total capacity will now be described. Referring to FIG. 1, the computer 28 receives a first state-of-charge estimate corresponding to a cell state-of-charge at a first time, and receives a second state-of-charge estimate corresponding to the cell state-of-charge at a second time. These estimates may be received from a subsystem within the computer configured to produce such estimates, or from some external source. The system 10 measures an integrated battery cell current between the first and second times. The computer 28 updates at least one recursive parameter based on the first state-of-charge estimate, the second state-of-charge estimate, and the integrated battery cell current. Further, the computer 28 determines an estimated battery cell total capacity based on at least one recursive parameter and stores a value corresponding to the estimated battery cell total capacity in the memory 45.

To further understand the deficiency of the methods used in the state of the art, consider FIG. 2. FIG. 2(a) illustrates the standard assumptions used when using least-squares linear regression to fit a straight line to data. Sample data points are drawn as squares, and error 9 bars indicating one standard deviation of measurement error are drawn as well. The x component of the data is assumed to have no measurement error, so the line is fit to minimize the sum of squares of all vertical distances between measured points and the regression line itself.

FIG. 2(b) illustrates the correct way to use regression to fit a straight line to data when both the x and y components have measurement error. This procedure is sometimes called “total least squares.” Notice that there are error bars associated with both the horizontal and vertical components of the measurement point. The regression line is the best fit when it minimizes the sum of squares of perpendicular distances between measured points and the regression line itself, not simply the vertical distances between measured points and the regression line.

FIG. 3 illustrates a representative example of perpendicular distance by denoting the distance between a measurement point and the regression line by R. Note that the slope of the regression line, which is what is desired to be computed, is C=Δy/Δx. Furthermore, if θ is defined as a tan(C), then R=Δy cos(θ). Using the trigonometric identity that

$\cos \left(\theta \right)=\frac{1}{\sqrt{1+{\tan }^2\left(\theta \right)}},$

and desiring to minimize the sum of squares of all perpendicular distances R_i from points indexed by i to the regression line, the total cost J must be minimized, where

$\begin{array}{cc}J=\sum _iR_i^2=\frac{1}{1+C^2}\sum _i{\left(y_i-Cx_i\right)}^2& \left(3\right)\end{array}$

To minimize J, 0 is set to ∂J/∂C and it is found that

$\begin{array}{cc}0=\frac{1}{1+C^2}\sum _i-2\left(y_i-Cx_i\right)x_i-\frac{2C}{{\left(1+C^2\right)}^2}\sum _i{\left(y_i-{\mathrm{Cx}}_i\right)}^2& \left(4\right)\end{array}$

which may be written as a quadratic equation

$0=s_2C^2+s_1C+s_0$ $\mathrm{where}$ $s_2=\sum _ix_iy_i,s_1=\sum _ix_i^2-y_i^2,\mathrm{and}s_0=-s_2.$

The slope of the line may then be solved by computing the solution to a quadratic equation

$\begin{array}{cc}C=\frac{-s_1+\sqrt{s_1^2+4s_2^2}}{2s_2}& \left(5\right)\end{array}$

Note that since s₀=−s₂ there is no need to compute both terms independently.

Generalizing this to the case with unequal errors Δy and Δx is straightforward: it merely involves pre-scaling they values by multiplying by some constant K so that the errors are equal, computing a preliminary value of the slope C′, and computing a final value of the slope as C=C′K. (A similar and equivalent procedure pre-scales the x values instead). Furthermore, generalizing to the case where different data points are given different confidence levels is straightforward. Furthermore, the math lends itself very naturally to a recursive algorithm for estimating the slope of the total least squares regression line, as will be described in a number of exemplary embodiments, below.

In one exemplary embodiment, the computer receives a first estimate of the battery cell SOC, z(t₁), at a first time, t₁. Furthermore, it receives a second estimate of the battery cell SOC, z(t₂), at a second time t₂. It computes x=z(t₂)−z(t₁). It measures an integrated battery cell current

$C=\frac{-s_1+\sqrt{s_1^2+4s_2^2}}{2s_2}$

It updates a first recursive parameter s₁=s₁+x²−y² and updates a second recursive parameter s₂=s₂+xy. Finally, it computes an estimate of battery cell total capacity as

$y={\int }_{t_1}^{t_2}\eta i\left(t\right)t.$

and stores that value in a memory. (s₁ and s₂ may be initialized to values indicative of a battery cell with nominal total capacity.) This exemplary embodiment is illustrated in FIG. 4.

In another exemplary embodiment, the first and second recursive parameters are instead computed as s₁=as₁+x²−y² and s₂=βs₂+xy, where 0<α<1 and 0<β<1.

This has the effect of placing more emphasis on the most recent data samples and placing less emphasis on the oldest data samples, allowing the estimate of battery cell total capacity to adapt more quickly to a changing total capacity.

In another exemplary embodiment, a weighted battery cell current is computed as

$y=K{\int }_{t_1}^{t_2}\eta i\left(t\right)t,$

where K is chosen to make the errors in x and y of similar size. The recursive parameters are calculated as before, using this new definition for y, and an initial battery cell total capacity estimate is calculated as

$C=\frac{-s_1+\sqrt{s_1^2+4s_2^2}}{2s_2}.$

The battery cell total capacity estimate is then computed as C=C′/K.

In yet another exemplary embodiment, individual sequential estimates of battery cell total capacity are computed and denoted as C_k⁽ⁱ⁾. The sequence is indexed by variable k The corresponding battery cell total capacity estimate C_k is then computed as

$\begin{array}{cc}{C}_k=\sum _{j=0}^{n_b}b_jC_{k-j}^{\left(i\right)}+\sum _{j=1}^{n_a}a_jC_{k-j}& \left(6\right)\end{array}$

where a_j, j=1L and b_j, j=1L n_b, are constants designed to implement a smoothing filter. FIG. 5 shows results using this exemplary embodiment. The line 70, referred to as the “intermediate estimate,” shows initial capacity estimates C_k^(j) over time for a lengthy battery test involving multiple full discharges from a fully charged state. The line 80, referred to as the “filtered estimate,” shows battery cell total capacity estimates C_k as the output of Eq. (6). In this particular example, an initial incorrect capacity estimate of 9.0 Ah is corrected to a value of 9.15 Ah in about 10 hours of driving. Since capacity is expected to change very slowly (perhaps 20 percent in a decade of use), this speed is more than sufficient.

The system and method for determining an estimated battery cell total capacity provides a substantial advantage over other systems and methods. In particular, the system and method provide a technical effect of accurately determining an estimated battery cell total capacity that is computationally efficient to compute, and does not impose operational constraints on how the battery cell is used when computing the estimate of battery cell total capacity.

The above-described methods can be embodied in the form of computer program code containing instructions embodied in tangible media, such as floppy diskettes, CD ROMs, hard drives, or any other computer-readable storage medium, wherein, when the computer program code is loaded into and executed by a computer, the computer becomes an apparatus for practicing the invention. The above-described methods can also be embodied in the form of computer program code, for example, whether stored in a storage medium, loaded into and/or executed by a computer, or transmitted over some transmission medium, loaded into and/or executed by a computer, or transmitted over some transmission medium, such as over electrical wiring or cabling, through fiber optics, or via electromagnetic radiation, wherein, when the computer program code is loaded into an executed by a computer, the computer becomes an apparatus for practicing the methods. When implemented on a general-purpose microprocessor, the computer program code segments configure the microprocessor to create specific logic circuits.

While the invention is described with reference to exemplary embodiments, it will be understood by those skilled in the art that various changes may be made and equivalent elements may be substituted for elements thereof without departing from the scope of the invention. In addition, many modifications may be made to the teachings of the invention to adapt to a particular situation without departing from the scope thereof. Therefore, it is intended that the invention not be limited the embodiments disclosed herein, but that the invention includes all embodiments falling with the scope of the appended claims. Moreover, the use of the terms first, second, etc. does not denote any order of importance, but rather the terms first, second, etc. are used to distinguish one element from another.

Claims

1. A method for determining an estimated battery cell total capacity indicative of a total capacity of a battery cell, comprising:

receiving a first battery cell state-of-charge estimate at a first time and receiving a second battery cell state-of-charge estimate at a second time subsequent to the first time;

measuring an integrated battery cell current value indicative of the integrated battery cell current between the first time and the second time;

updating at least one recursive parameter based on the first battery cell state-of-charge estimate, the second battery cell state-of-charge estimate, and the integrated battery cell current value;

determining the estimated battery cell total capacity based on at least one recursive parameter; and

storing a value corresponding to the estimated battery cell total capacity in a memory.

2. The method of claim 1, wherein updating at least one recursive parameter based on the first battery cell state-of-charge estimate, the second battery cell state-of-charge estimate, and the integrated battery cell current value comprises:

computing a state-of-charge difference by subtracting the first battery cell state-of-charge estimate from the second battery cell state-of-charge estimate;

updating a first recursive parameter by multiplying its prior value by a first scaling factor between 0 and 1 (inclusive);

further updating the first recursive parameter by adding the square of the state-of-charge difference and by subtracting the square of the integrated battery cell current value;

updating a second recursive parameter by multiplying its prior value by a second scaling factor between 0 and 1 (inclusive); and

further updating the second recursive parameter by adding the product of the state-of-charge difference and the integrated battery cell current value.

3. The method of claim 1, wherein determining the estimated battery cell total capacity based on at least one recursive parameter comprises:

receiving a first recursive parameter and a second recursive parameter; and

solving a quadratic equation involving the first recursive parameter and the second recursive parameter to determine the estimated battery cell total capacity.

4. The method of claim 1, wherein determining the estimated battery cell total capacity based on at least one recursive parameter comprises:

receiving a first recursive parameter and a second recursive parameter;

solving a quadratic equation involving the first recursive parameter and the second recursive parameter to determine an initial battery cell total capacity estimate;

determining a filtered battery cell total capacity estimate based on at least one of an initial battery cell total capacity estimate corresponding to a present or prior time and a filtered battery cell total capacity estimate corresponding to a prior time; and

determining an estimated battery cell total capacity based on the filtered battery cell total capacity estimate.

5. The method of claim 1 where measuring an integrated battery cell current value indicative of the integrated battery cell current between the first time and the second time comprises:

measuring an initial integrated battery cell current value indicative of the integrated battery cell current between the first time and the second time; and

computing the integrated battery cell current by multiplying the initial integrated battery cell current value by a scaling constant.

6. A system for determining an estimated battery cell total capacity indicative of a total capacity of a battery cell, comprising:

a first subsystem configured to receive a first state-of-charge estimate at a first time;

a second subsystem configured to receive a second state-of-charge estimate at a second time, subsequent to the first time;

a current sensor electrically coupled to the battery cell, the current sensor configured to generate a signal indicative of integrated battery cell current between the first and second times;

a computer configured to receive the first state-of-charge estimate, the second state-of-charge estimate, and the integrated battery cell current, the computer is further configured to update at least one adaptive parameter based on the first state-of-charge estimate, the second state-of-charge estimate, and the integrated battery cell current, and to determine the estimated battery cell total capacity based on at least one adaptive parameter, the computer further configured to store a value corresponding to the estimated battery cell total capacity in a memory.

7. A computer readable medium having computer-executable instructions for performing a method for determining an estimated battery cell total capacity indicative of a total capacity of a battery cell, the method comprising:

receiving a first battery cell state-of-charge estimate at a first time and receiving a second battery cell state-of-charge estimate at a second time subsequent to the first time;

measuring an integrated battery cell current value indicative of the integrated battery cell current between the first time and the second time;

updating at least one recursive parameter based on the first battery cell state-of-charge estimate, the second battery cell state-of-charge estimate, and the integrated battery cell current value;

determining the estimated battery cell total capacity based on at least one recursive parameter; and

storing a value corresponding to the estimated battery cell total capacity in a memory.

8. The method of claim 1, wherein the step of measuring an integrated battery cell current value indicative of the integrated battery cell current between the first time and the second time uses the relationship y = ∫ t 1 t 2  η   i  ( t )    t.

9. The method of claim 1, wherein the step of updating at least one recursive parameter based on the first battery cell state-of-charge estimate, the second battery cell state-of-charge estimate, and the integrated battery cell current value uses the relationship s1=s1+(z(t2)−z(t1))2−y2.

10. The method of claim 2, wherein the step of updating the second recursive parameter uses the relationship s2=s2+(z(t2)−z(t1))y.

11. The method of claim 3, wherein the step of computing an estimated battery cell total capacity is governed by the equation C = - s 1 + s 1 2 + 4  s 2 2 2  s 2.

12. The method of claim 3, wherein the first recursive parameter is governed by the equation s1=as1+x2−y2, and the second recursive parameter is governed by the equation s2=βs2+xy, where 0<α<1 and 0<β<1.