Method and device for detecting battery cell states and battery cell parameters

Abstract

Device (1) and method for detecting battery cell states, BZZ, and/or battery cell parameters, BZP, of at least one battery cell (BZ), comprising a dual Kalman filter (2) which includes a state estimator (2A) for estimating battery cell states, BZZ, and a parameter estimator (2B) for estimating battery cell parameters, BZP, and comprising a determination unit (3) which is suitable for determining noise components (n, v) of the state estimator (2A) and of the parameter estimator (2B) on the basis of a stored characteristic parameter behaviour of the battery cell (BZ), wherein the battery cell states, BZZ, and the battery cell parameters, BZP, can be adapted automatically to a specified battery model (BM) of the battery cell (BZ) by means of the dual Kalman filter (2) on the basis of the noise components (n, v) determined by the determination unit.

Claims

1. A method for detecting battery cell states and battery cell parameters of at least one battery cell, comprising the steps of: (a) determining first noise components of a state estimator of a dual Kalman filter and determining second noise components of a parameter estimator of the dual Kalman filter on the basis of a characteristic parameter behaviour of the battery cell in relation to at least one battery cell state in dependence upon a change in the battery cell state over time or in dependence upon a read-out parameter behaviour in relation to at least one battery cell parameter; and (b) adapting the battery cell states and the battery cell parameters to a pre-determined battery model of the battery cell by means of the dual Kalman filter on the basis of the determined first and second noise components, wherein the determined first and second noise components include a process noise and a measurement noise, wherein the process noise includes a process noise of the battery cell parameters and/or a process noise of the battery cell states, wherein the process noise of the battery cell parameters is determined in dependence upon its read-out characteristic parameter behaviour in relation to a battery cell state and in dependence upon a change in this battery cell state over time, and wherein the process noise of the battery cell states is determined in dependence upon its read-out parameter behaviour in relation to at least one battery cell parameter; (c) outputting the adapted battery cell states and battery cell parameters to a control unit; and (d) operating the control unit to control a load connected to the battery cell and/or a current source connected to the battery cell in response to the adapted battery cell states and battery cell parameters.

2. The method as claimed in claim 1, wherein: the battery cell states include a state of charge of the battery cell and/or dynamic battery cell states including a diffusion voltage of the battery cell, and the battery cell parameters include one or more of an internal resistance of the battery cell, a rated capacitance of the battery cell, resistive dynamic components including a diffusion resistance of the battery cell, and capacitive dynamic components including a diffusion capacitance of the battery cell.

3. The method as claimed in claim 1, wherein a state vector of the state estimator of the dual Kalman filter includes the battery cell states of the battery cell and a state vector of the parameter estimator of the dual Kalman filter includes the battery cell parameters of the battery cell.

4. The method as claimed in claim 1, further comprising controlling a loading of the battery cell on the basis of the detected battery cell states and the detected battery cell parameters.

5. The method as claimed in a claim 1, further comprising determining the characteristic parameter behaviour of the battery cell on the basis of measurement variables which are detected by means of sensors and include one or more of a terminal current, a terminal voltage and a temperature of the battery cell, wherein the characteristic parameter behaviour of the battery cell indicates, for each of the battery cell parameters its average value and/or its variance in relation to each of the battery cell states and/or on the basis of measurement variables of the battery cell which are detected by means of sensors.

6. The method as claimed in claim 1, wherein said process noise of the battery cell parameters is determined on the basis of measurement variables detected by means of sensors.

7. The method as claimed in claim 1, further comprising, for each of measurement variables detected by means of sensors, calculating a measurement variable noise of the measurement variable on the basis of an average value and/or a variance of a noise behaviour of the corresponding measurement variable sensor in relation to the measurement variable which is read out from a data store, wherein the calculated measurement variable noises of all the detected measurement variables are added in order to determine the measurement noise.

8. The method as claimed in claim 1, further comprising monitoring a stability of the dual Kalman filter and the adaptation of the respective battery cell states by means of the state estimator of the dual Kalman filter, and suspending the adaptation of the respective battery cell parameters by means of the parameter estimator of the dual Kalman filter as soon as a threat of instability of the dual Kalman filter is recognised, wherein the stability of the dual Kalman filter is monitored on the basis of a covariance behaviour of the battery cell states and a covariance behaviour of the battery cell parameters.

9. The method as claimed in a claim 8, further comprising comparing the covariance behaviour of the battery cell state and/or the covariance behaviour of the battery cell parameter of the dual Kalman filter with an associated stored desired covariance behaviour, wherein the adaptation of a battery cell state by means of the state estimator of the dual Kalman filter is suspended as soon as its covariance behaviour exceeds the associated stored desired covariance behaviour of the respective battery cell state and wherein the adaptation of a battery cell parameter by means of the parameter estimator of the dual Kalman filter is suspended as soon as its covariance behaviour exceeds the associated stored desired covariance behaviour of the respective battery cell parameter wherein the associated stored desired covariance behaviour includes a funnel function which decreases exponentially over time and of which the coefficients are configured.

10. A device for detecting battery cell states and/or battery cell parameters of at least one battery cell, comprising: a dual Kalman filter, in communication with the control unit, which includes a state estimator for estimating battery cell states and a parameter estimator for estimating battery cell parameters; and a determination unit which is arranged to determine noise components of the state estimator and of the parameter estimator on the basis of a stored characteristic parameter behaviour of the battery cell in relation to at least one battery cell state in dependence upon a change in the battery cell state over time or in dependence upon a read-out parameter behaviour in relation to at least one battery cell parameter, wherein the dual Kalman filter is arranged to adapt the batter cell states and the battery cell parameters automatically to a pre-determined battery model of the battery cell on the basis of the noise components determined by the determination unit, wherein the determined noise components include a process noise and a measurement noise, wherein the process noise includes a process noise of the battery cell parameters and/or a process noise of the battery cell states, wherein the process noise of the battery cell parameters is determined in dependence upon its read-out characteristic parameter behaviour in relation to the battery cell states and in dependence upon a change in the battery cell state over time, wherein the process noise of the battery cell states is determined in dependence upon its read-out parameter behaviour in relation to at least one battery cell parameter, and wherein the adapted battery cell states and the battery cell parameters of the battery cell are output to the control unit which is operated to control a load connected to the battery cell and/or a current source connected to the battery cell in response to the adapted battery cell states and the battery cell parameters.

11. The device as claimed in claim 10, wherein the state estimator of the dual Kalman filter is formed by a first Kalman filter and the parameter estimator of the dual Kalman filter is formed by a second Kalman filter.

12. The device as claimed in claim 11, wherein the two Kalman filters of the dual Kalman filter each include a linear Kalman filter, an extended Kalman filter, an unscented Kalman filter, a square-root unscented Kalman filter or a central-difference Kalman filter.

13. The device as claimed in claim 10, wherein the load is an electric motor.

Description

(1) Possible embodiments of the various aspects of the invention will be explained in greater detail hereinafter with reference to the enclosed figures.

(2) In the figures:

(3) FIG. 1 shows an electrical equivalent circuit for a battery cell, as can be used in the case of a conventional battery management system;

(4) FIG. 2 shows a simple flow diagram to illustrate one possible exemplified embodiment of a method in accordance with the invention for detecting battery cell states and battery cell parameters of at least one battery cell according to one aspect of the present invention;

(5) FIG. 3 shows a block diagram to illustrate an exemplified embodiment of a device in accordance with the invention for detecting battery cell states and/or battery cell parameters of a battery cell according to a further aspect of the invention;

(6) FIG. 4 shows a block diagram to illustrate an exemplified embodiment of a Kalman filter which can be used in the case of the method in accordance with the invention and in the case of the device in accordance with the invention;

(7) FIG. 5 shows a further block diagram to illustrate the structure of a dual Kalman filter, as can be used in the case of the method in accordance with the invention and the device in accordance with the invention;

(8) FIG. 6 shows an example of a battery model of a battery cell which can be used in the case of the method in accordance with the invention;

(9) FIGS. 7A, 7B show by way of example one possible noise behaviour of voltage and current sensors in the case of a battery management system in accordance with the invention;

(10) FIGS. 8A, 8B show by way of example a noise behaviour of a battery model used in the case of the method in accordance with the invention;

(11) FIGS. 9A to 9D show diagrams to illustrate a temperature and state of charge dependency of parameters of the battery model;

(12) FIGS. 10 to 10D show diagrams to illustrate a normalised discharge capacity in the case of various conventional lithium-ion batteries;

(13) FIGS. 11A to 11D show diagrams to illustrate a variance of the change in the battery cell parameters and a state change in relation to the state of charge of a battery cell;

(14) FIGS. 12A to 12D show diagrams to illustrate the variance of the change in a battery cell capacitance in relation to a terminal current in the case of a battery cell;

(15) FIG. 13 shows an image to illustrate a theoretical progression of a covariance;

(16) FIG. 14 shows by way of example the progression of a diffusion voltage within a funnel function, as can be used in the case of the method in accordance with the invention and the device in accordance with the invention for stability monitoring.

(17) As can be seen in FIG. 2, the method in accordance with the invention for detecting battery cell states BZZ and battery cell parameters BZP comprises essentially two steps.

(18) In a first step S1, noise components n, v of a state estimator of a dual Kalman filter and noise components n, v of a parameter estimator of the dual Kalman filter are determined on the basis of a characteristic parameter behaviour of a battery cell in relation to at least one battery cell characteristic variable BZK of the battery cell BZ.

(19) In a further step S2, battery cell states BZZ and battery cell parameters BZP are adapted to a specified battery model BM of the battery cell BZ by means of the dual Kalman filter on the basis of the noise components n, v determined in step S1.

(20) In the case of one possible embodiment, the noise components n, v determined in step S1 include a process noise v and a measurement noise n. The determined process noise v can include a process noise v.sub.BZZ of the battery cell states BZZ and/or a process noise v.sub.BZP of the battery cell parameters BZP.

(21) The dual Kalman filter DKF used in the case of the method in accordance with the invention as shown in FIG. 1 has two Kalman filters KF coupled together, wherein a first Kalman filter is used as a state estimator and a second Kalman filter is used as a parameter estimator. A state vector of the state estimator of the dual Kalman filter DKF includes the battery cell states x.sub.BZZ of the battery cell BZ. A state vector of the parameter estimator of the dual Kalman filter DKF includes the battery cell parameters x.sub.BZP of the battery cell BZ. In the case of one possible embodiment, the characteristic parameter behaviour of the battery cell BZ as used in step S1 is read out from a data store. In the case of one possible embodiment, the characteristic parameter behaviour of the battery cell BZ includes, for each battery cell parameter BZP, its average value μ.sub.BZP and/or its variance σ.sub.BZP in relation to each battery cell state BZZ of the battery cell BZ. In the case of one embodiment, the stored characteristic parameter behaviour of the battery cell BZ can be determined on the basis of measurement variables MG detected by means of sensors. These measurement variables MG include e.g. a terminal current I.sub.KL, a terminal voltage U.sub.KL and/or a temperature of the battery cell BZ.

(22) In the case of one possible embodiment, in step S1 a process noise v.sub.BZP of a battery cell parameter BZP is determined in dependence upon its read-out characteristic parameter behaviour in relation to a battery cell state BZZ and/or on the basis of measurement variables MG detected by means of sensors in dependence upon a change in the battery cell state BZZ over time and/or on the basis of measurement variables detected by means of sensors.

(23) In the case of one possible embodiment, the battery cell states BZZ of the battery cell BZ include a state of charge SOC of the battery cell and/or dynamic battery cell states BZZ, in particular a diffusion voltage U.sub.d of the battery cell BZ. The state of charge SOC of the battery cell BZ is between 0 and 1 or between 0 and 100%. The diffusion voltage U.sub.d is the particular voltage in a battery model BM, which is illustrated e.g. in FIG. 6, which decreases over the RC member and thereby demonstrates a dynamic behaviour. The battery cell states BZZ are the particular physical variables of the battery model BM which also occur as states in differential equations.

(24) In one possible embodiment, the battery cell parameters BZP of the battery cell BZ include an internal resistance R.sub.i of the battery cell, a rated capacitance C.sub.n of the battery cell Bz and resistive dynamic components, in particular a diffusion resistance R.sub.d of the battery cell and/or capacitive dynamic components, in particular a diffusion capacitance C.sub.D of the battery cell BZ. The battery cell parameters BZP include the particular physical variables of the battery model BM which occur as a parameter or value in the differential equations. The rated capacitance C.sub.n forms a battery cell parameter BZP and indicates the present capacitance of the battery cell BZ based upon the present influence variables. The diffusion resistance R.sub.d likewise represents a battery cell parameter BZP and indicates the resistive component of the dynamic RC component of the battery model BM used, which is illustrated by way of example in FIG. 6. The diffusion capacitance C.sub.D of the battery cell BZ likewise forms a battery cell parameter BZP and represents the capacitive component of the dynamic RC component within the battery model shown in FIG. 6.

(25) In one possible embodiment, the battery cell characteristic variables BZK of the battery cell BZ include a loading characteristic variable of the battery cell BZ, a temperature T of the battery cell BZ and/or an ageing characteristic variable of the battery cell BZ. An example of a loading variable of the battery cell BZ which is used as the battery cell characteristic variable BZK is the so-called C-rate. The C-rate is the loading multiple of the rated capacitance and describes the loading of the battery cell BZ. For example, in the case of a capacitance of the battery cell BZ of 1 Ah, 1C denotes a loading of 1 A. For example, cycle numbers or characteristic variables in terms of the ageing of the battery cell can be used as ageing characteristic variables of the battery cell BZ.

(26) In the case of one possible embodiment, the noise components n, v used in step S1 include a process noise v and/or a measurement noise n. In one possible embodiment, the process noise can include a process noise v.sub.BZZ of the various battery cell states BZZ and/or a process noise v.sub.BZP of the various battery cell parameters BZP. In the case of one possible embodiment, the process noise includes a process noise of a state of charge SOL of the battery cell BZ and/or a process noise of a diffusion voltage U.sub.d of the battery cell BZ. Furthermore, the process noise can also include a process noise v.sub.BZP of various battery cell parameters BZP, e.g. an internal resistance of the battery cell BZ, a rated capacitance C.sub.N of the battery cell BZ, a diffusion resistance R.sub.d of the battery cell BZ and/or a diffusion capacitance C.sub.D of the battery cell BZ.

(27) The process noise v.sub.BZP of the battery cell parameter BZP is determined in dependence upon its read-out characteristic parameter behaviour in relation to a battery cell state BZZ. Alternatively, the process noise can also be determined on the basis of measurement variables MG detected by means of sensors, and in dependence upon a change in a battery cell state BZZ over time. The process noise v.sub.BZZ of a battery cell state BZZ is determined in dependence upon its read-out parameter behaviour in relation to at least one battery cell parameter BZP. In the case of one possible embodiment of the method in accordance with the invention, for each measurement variable MG detected by means of sensors, its measurement variable noise is calculated on the basis of an average value μ.sub.MG and/or a variance σ.sub.MS of a noise behaviour of the corresponding measurement variable sensor in relation to the measurement variable MG. The calculated measurement variable noise of all measurement variables can be added in order to determine the measurement noise n.

(28) FIG. 3 shows a block diagram to illustrate an exemplified embodiment of a device 1 in accordance with the invention for detecting battery cell states BZZ and/or battery cell parameters BZP of at least one battery cell BZ. As can be seen in FIG. 3, the detection device 1 in accordance with the invention has a dual Kalman filter 2 in the illustrated exemplified embodiment. The dual Kalman filter 2 consists of two coupled Kalman filters. The first Kalman filter 2A of the dual Kalman filter 2 is used as a state estimator ZS for estimating battery cell states BZZ. The second Kalman filter 2B of the dual Kalman filter 2 is used as a parameter estimator PS for estimating battery cell parameters BZP of the battery cell BZ. Furthermore, the detection device 1 has a determination unit 3 which is suitable for determining noise components n, v of the state estimator 2A and of the parameter estimator 2B on the basis of a stored characteristic parameter behaviour of the battery cell BZ. The characteristic parameter behaviour of the battery cell BZ can be stored e.g. in a data store 4. The battery cell states BZZ and the battery cell parameters BZP can be adapted automatically to a specified battery model BM of the battery cell BZ by means of the dual Kalman filter 2 on the basis of the noise components n, v determined by the determination unit 3. The state estimator ZS of the dual Kalman filter 2 is formed by a first Kalman filter 2A and the parameter estimator PS of the dual Kalman filter 2 is formed by a second Kalman filter 2B. In the case of one possible embodiment of the detection device 1 in accordance with the invention, the two Kalman filters 2A, 2B of the dual Kalman filter 2 are each formed by a linear Kalman filter LKF, an extended Kalman filter EKF, an unscented Kalman filter UKF, a square-root unscented Kalman filter SR-UKF or by a central-difference Kalman filter CD-KF. The Kalman filter 2 provides estimated values for system state variables x,
{dot over (x)}=A.Math.x+B.Math.u+v
y=C.Math.x+n
where x is the state vector,
A represents the system matrix,
B represents the input matrix,
u is the deterministic input vector or control vector,
v represents a Gaussian-distributed process noise,
y forms an output value which is produced from the output matrix
C and the Gaussian-distributed measurement noise n.

(29) In the case of one possible embodiment, the terminal current I.sub.KL and/or the terminal voltage U.sub.KL of the battery cell BZ can be used as the deterministic input variable or control vector u.

(30) FIG. 4 shows the structure of a linear discrete Kalman filter 2. The linear discrete Kalman filter KF illustrated in FIG. 4 can be used as a state estimator 2A and/or as a parameter estimator 2B of the dual Kalman filter 2. The upper part of the Kalman filter 2 illustrated in FIG. 4 depicts the battery model BM of the battery cell BZ, whereas the lower part of the Kalman filter KF as an adaptive filter minimises errors in the states, in particular in the battery cell states BZZ, by means of Kalman amplification. x represents the state vector of the battery cell BZ and y represents the output vector of the battery cell BZ. Located at the input of the Kalman filter KF is a deterministic control vector u, such as e.g. the measured terminal voltage U.sub.KL and/or the measured terminal current I.sub.KL of the battery cell BZ.

(31) For non-linear systems, discretization via the Laplace transform is not possible. The output form is a state space representation which consists of non-linear first-order differential equations:
{dot over (x)}=f(x,u)+v
y=h(x,u)+n
wherein v is the process noise, n represents the measurement noise,
x is the state vector of the battery cell,
x is the output vector of the battery cell,
u is the deterministic input vector, and
f, h represent non-linear functions.

(32) In the case of one possible embodiment, the Kalman filter KF used as a state estimator ZS or parameter estimator PS is formed by a linear Kalman filter. In the case of a further possible embodiment, extended Kalman filters EKF can also be used. In this case, the discrete system matrix A is approximated from a Taylor series. The continuous system matrix A and the output matrix C are approximated with the aid of a first-order approximation (Jacobi matrix) from the non-linear differential equations.

(33) In the case of a further possible alternative embodiment, the two Kalman filters KF of the dual Kalman filter 2 are formed by unscented Kalman filters UKF. The unscented Kalman filter UKF belongs to the group of sigma-point Kalman filters. In the case of this type of filter, the problem of the EKF in terms of a linearisation and the precision limited thereby is solved by deduction-free approaches. In this case, the system equations are not approximated via Taylor series, but instead are calculated by function evaluations at various points, the so-called sigma points. The sigma points are selected as inputs of the non-linear function such that they precisely detect the average value and the covariance of the sought-after state.

(34) In the case of a further possible embodiment of the detection device 1 in accordance with the invention, the two Kalman filters 2A, 2B of the dual Kalman filter 2 are formed by square-root unscented Kalman filters SR-UKF. A disadvantage when using UKF filters is that of calculating a necessary square root of the sigma points S.sub.k in order to determine the state covariance P.sub.k=S.sub.kSTk. For example, a Cholesky factorisation which, however, requires significant computational effort can be used for the sigma points. The sigma points S.sub.k are carried directly with a square-root unscented Kalman filter, without a re-factorisation being required at every time segment. The Kalman filter thereby becomes more stable because the state covariance assumes positive-semi-defined values and thus has a similarly recursive behaviour as the linear Kalman filter. For this purpose, in the case of one possible embodiment the square-root unscented Kalman filter SR-UKF can employ various linearly algebraic techniques, e.g. QR decomposition, the updating of the Cholesky factor or the method of least squares.

(35) In the case of a further possible embodiment of the detection device 1 in accordance with the invention, the two Kalman filters 2A, 2B of the dual Kalman filter 2 can each be formed by a central-difference Kalman filter CD-KF. A central-difference Kalman filter CD-KF is based upon an interpolation formula by Stirling and is similar in terms of its circuitry-wise structure to the unscented Kalman filter UKF.

(36) FIG. 5 shows a structure of a dual Kalman filter 2 which can be used in the detection device 1 in accordance with the invention. The dual Kalman filter KF includes a state estimator 2A and a parameter estimator 2B. Both the state estimator 2A and the parameter estimator 2B can be formed by a Kalman filter 2, wherein each Kalman filter 2 includes a prediction part PT and a correction part KT. The state estimator 2A is used for estimating battery cell states BZZ of the battery cell BZ. The parameter estimator 2B of the dual Kalman filter 2 is used for estimating battery cell parameters BZP of the battery cell BZ. The battery cell states BZZ and the battery cell parameters BZP are adapted automatically to a specified battery model BM of the battery cell BZ by means of the dual Kalman filter 2 on the basis of the determined noise components n, v. These noise components n, v are determined by the determination unit 3 illustrated in FIG. 3.

(37) In one possible embodiment, the battery cell states BZZ adapted by the state estimator 2A of the dual Kalman filter 2 and the battery cell parameters BZP of the battery cell BZ which are adapted by the parameter estimator 2B of the dual Kalman filter 2 are output to a control unit 5, as illustrated in FIG. 3. In one possible embodiment, this control unit 5 is used for controlling a load, e.g. an electric motor, connected to the battery cell BZ. Furthermore, in an alternative embodiment, the control unit 5 can be provided for controlling a current source connected to the battery source BZ. As illustrated in FIG. 3, in the case of one possible embodiment, the control unit 5 provides control signals CRTL-RL for an electric load 6 connected to the battery cell BZ. This load 6 can be e.g. an electric motor of an electric vehicle. Furthermore, in the exemplified embodiment shown in FIG. 3, the control unit 5 provides a control signal CTRL-SQ for a current source 7 connected to the battery cell BZ. The current source 7 is connected to a first battery cell terminal 8 of the battery cell BZ. The load 6 is connected to a second battery cell terminal 9 of the battery cell BZ.

(38) In the case of the exemplified embodiment illustrated in FIG. 3, various sensors are provided on the battery cell BZ. In the case of one possible embodiment, these sensors include a voltage sensor 10, a current sensor 11 and a temperature sensor 12. The voltage sensor 10 measures a terminal voltage U.sub.KL between the two battery cell terminals 8, 9. In the case of the illustrated exemplified embodiment, the current sensor 11 is used for measuring a terminal current I.sub.KL of the battery cell BZ. Furthermore, the temperature sensor 12 is used for measuring a temperature T of the battery cell BZ. The terminal voltage U.sub.KL, the terminal current I.sub.KL and the temperature T of the battery cell BZ represent measurement variables MG detected by means of sensors. These measurement variables MG are supplied as deterministic input measurement variables to the dual Kalman filter 2, as can be seen in the block diagram shown in FIG. 3. The measurement variables MG, i.e. the terminal voltage U.sub.KL, the terminal current I.sub.KL and/or the temperature T of the battery cell BZ are also supplied to the determination unit 3. The determination unit 3 determines the noise components n, v of the state estimator 2A and of the parameter estimator 2B on the basis of a stored characteristic parameter behaviour of the battery cell BZ. In the case of one possible embodiment, the data store 4 stores, for each battery cell parameter BZP, its average value μ.sub.BZP, it variance σ.sub.BZP in relation to each battery cell state BZZ and the measurement variables MG of the battery cell BZ which are detected by means of sensors. The process noise v.sub.BZP of the battery cell parameters BZP is determined by the determination unit 3 in dependence upon its read-out characteristic parameter behaviour in relation to a battery cell state BZZ on the basis of the measurement variables MG detected by means of sensors. Furthermore, a process noise v.sub.BZZ of the battery cell state BZZ is determined in dependence upon its read-out parameter behaviour by means of the determination unit 3 in relation to at least one battery cell parameter BZP. For each measurement variable MG detected by means of sensors, the determination unit 3 calculates its measurement variable noise n on the basis of an average value μ.sub.MG and/or a variance σ.sub.MS of a noise behaviour of the corresponding measurement variable sensor 10, 11, 12 in relation to the measurement variable MG which are read out from the data store 4. The calculated measurement variable noise of the various measurement variables MG can be added in order to determine the measurement noise n by means of the determination unit 3. The determined measurement noise n is supplied to the dual Kalman filter 2 as shown in FIG. 3. The dual Kalman filter 2 comprises two Kalman filters 2A, 2B which can be adjusted in parallel to each other. The two Kalman filters 2A, 2B of the dual Kalman filter 2 can have the same input vector u.sub.k−1 and the same measurement variables y.sub.k. In the respective prediction part PT of the two Kalman filters 2A, 2B, the previously calculated states and parameters form the input variables of the current time step. The two Kalman filters 2A, 2B within the dual Kalman filter 2 can have the same filter type or can be formed by different filter types. The state vector x of the dual Kalman filter 2 includes the state vector of the state estimator 2A and the state vector of the parameter estimator 2B. The state vector x.sub.BZZ of the state estimator 2A includes the battery cell states BZZ of the battery cell BZ. The state vector of the parameter estimator 2B of the dual Kalman filter 2 includes the battery cell parameters x.sub.BZP of the battery cell BZ. Since the supplied measurement variables MG of the state estimator 2A and of the parameter estimator 2B are the same, only the additional dependency of the parameter vector in the non-linear output equation changes in the continuous state space representation of the parameter estimator 2B:
y=g(x,u,w)+n

(39) In this embodiment, the adaptation of the battery cell parameter BZP is effected primarily via the parameter process noise.

(40) In one possible embodiment of the method in accordance with the invention, stability of the dual Kalman filter 2 is monitored. An adaptation of the respective battery cell states BZZ by means of the state estimator 2A of the dual Kalman filter 2, and an adaptation of the respective battery cell parameters BZP by means of the parameter estimator 2B of the dual Kalman filter 2 is suspended as soon a threat of instability of the dual Kalman filter 2 is recognised. In the case of one possible embodiment, the stability of the dual Kalman filter 2 is monitored on the basis of a covariance behaviour P.sub.BZZ of the battery cell states BZZ and on the basis of a covariance behaviour P.sub.BZP of the battery cell parameters BZP. In one possible embodiment, the covariance behaviour P.sub.BZZ of the battery cell state BZZ and/or the covariance behaviour P.sub.BZP of a battery cell parameter BZP of the dual Kalman filter 2 is compared with an associated stored desired covariance. An adaptation of a battery cell state BZZ by means of the state estimator 2A of the dual Kalman filter 2 is suspended as soon as its covariance behaviour P.sub.BZZ exceeds the stored associated desired covariance behaviour of the respective battery cell state BZZ. Furthermore, the adaptation of a battery cell parameter BZP by means of the parameter estimator 2B of the dual Kalman filter 2 is suspended as soon as its covariance behaviour P.sub.BZP exceeds the stored associated desired covariance behaviour of the respective battery cell parameter BZP. In the case of one possible embodiment, the stored desired covariance behaviour includes a funnel function which decreases exponentially over time and of which the coefficients are configured.

(41) In the case of the exemplified embodiment illustrated in FIG. 3, the detection device 1 includes a stability monitoring unit 13 which is suitable for suspending the adaptation of a battery cell state BZZ by means of the state estimator 2A of the dual Kalman filter 2 and/or suspending the adaptation of a battery cell parameter BZP by means of the parameter estimator 2B of the dual Kalman filter 2 as soon as the corresponding covariance behaviour exceeds a stored associated desired covariance behaviour. In the case of one possible embodiment, the desired covariance behaviour is stored in a data store 14 of the detection device 1. In the case of one possible implementation, the stored desired covariance behaviour includes a funnel function which decreases exponentially over time, as illustrated in FIGS. 13, 14.

(42) FIG. 6 shows by way of example a battery model BM which can be used in the case of the method in accordance with the invention and the detection device 1 in accordance with the invention.

(43) The resistance R.sub.DC1s is a combined resistance which combines the real cell internal resistance R.sub.i of the battery cell BZ and all dynamic procedures occurring within 1 sec, e.g. charge transfer procedures. The resistance R.sub.D represents the diffusion resistance of the battery cell BZ. The capacitance C.sub.D represents the diffusion capacitance of the battery cell BZ. Furthermore, FIG. 6 illustrates the terminal voltage U.sub.KL of the battery cell BZ and the terminal current I.sub.KL of the battery cell BZ. Furthermore, the battery model BM shows a diffusion voltage U.sub.D which is applied across the diffusion resistance R.sub.D and the diffusion capacitance C.sub.D. U.sub.0(SOC) is the state of charge-dependent open-circuit voltage (OCV) of the battery cell.

(44) The state of charge SoC is defined as the integrated terminal current over the loading time period and the deduction of the state of charge SoC is the terminal current I.sub.KL(t) in relation to the capacitance of the battery cell C.sub.n.

(45) $\begin{array}{cc}S\stackrel{.}{O}C=\frac{1}{C_n}{I}_{KL}\left(t\right)& \left(1\right)\end{array}$

(46) The diffusion voltage U.sub.D can be calculated from a first-order differential equation:

(47) $\begin{array}{cc}{\stackrel{.}{U}}_D=-\frac{1}{R_D.Math.C_D}.Math.U_D\left(t\right)+\frac{1}{C_D}{I}_{KL}\left(t\right)& \left(2\right)\end{array}$

(48) The state space representation can now be made up therefrom as follows:

(49) $\begin{array}{cc}\left[\begin{array}{c}{\stackrel{.}{U}}_D\\ S\stackrel{.}{O}C\end{array}\right]=\left[\begin{array}{cc}\frac{1}{R_D{C}_D}& 0\\ 0& 0\end{array}\right]\left[\begin{array}{c}{U}_D\\ \mathrm{SOC}\end{array}\right]+\left[\begin{array}{c}\frac{1}{C_D}\\ \frac{1}{C_n}\end{array}\right]I_{KL}\left(t\right)& \left(3\right)\end{array}$

(50) The terminal voltage U.sub.KL is calculated from the internal voltage drops of the battery cell BZ, caused by the terminal current I.sub.KL as:
U.sub.KL(t)=U.sub.0(SOC)+.sub.RDc,1sI.sub.KL(t)+U.sub.D(t)  (4)

(51) In this simple exemplified embodiment, the state vector x(t) therefore contains:
[x.sub.1,x.sub.2].sup.T=[U.sub.D(t), SOC(t)].sup.T  (5)

(52) An essential variable, which can influence the behaviour of the Kalman filter 2, is the noise, i.e. the process noise v and the measurement noise n. The measurement noise n takes into consideration the errors in the measurement value and the effects thereof on the model adaptation traced therefrom. The process noise v takes into consideration the simplification and error tolerance of the battery model BM. Hereinunder, the respective noise term for the dual Kalman filter 2 is derived in order to carry out rapid and efficient adaptation in a possible real application.

(53) The consideration of the measurement noise n is carried out in the output equation:
y (t)=C.sub.x(.sub.t)+n   (6)

(54) When examining the output equation of the battery model BM used, various measurement noise components can be identified, specifically the measurement noise of the cell voltage measurement, the measurement noise of the current sensor and a noise value in relation to the diffusion voltage. From the variance σ of the respective terms, the measurement noise n can be derived:
U.sub.KL=R U.sub.0(SOC)+R.sub.DC,1s(I.sub.KL(t)+σ.sub.i.sub.kl)+(U.sub.D(t)+σ.sub.U.sub.D) +σ.sub.U.sub.KL
n=R.sub.DC,1sσ.sub.I.sub.KLσ.sub.U.sub.Dσ.sub.U.sub.KL   (7)

(55) The weighting of the three noise components can differ strongly from one to another.

(56) FIGS. 7A, 7B show by way of example the noise behaviour of voltage and current sensors of the battery management system in accordance with the invention. FIG. 7A shows the noise behaviour of a voltage sensor of the battery management system, e.g. of the voltage sensor 10 shown in FIG. 3. FIG. 7A shows the distribution of measurements with respect to the terminal voltage U.sub.KL. For example, the average value μis 4.01 V and the variance is 6.546e-0.5.

(57) FIG. 7B shows by way of example the behaviour of a current sensor within the battery management system, e.g. of the current sensor 11 shown in FIG. 3. It is possible to see the number of measurements distributed over different terminal currents I.sub.KLThe average value μ is e.g. μ=−1.14 A and the variance σ is 1.389e-01.

(58) FIGS. 8A, 8B show a noise behaviour of the battery model BM. FIG. 8A shows an error between measurement and simulation over time t. FIG. 8A shows an absolute error in the battery modellings in relation to a true cell voltage of the battery cell BZ.

(59) FIG. 8B shows a histogram of the deviation. FIG. 8B thus shows an associated noise behaviour. FIG. 8B shows e.g. the number of measurements for various voltages U.

(60) Since the internal resistance R.sub.DCIs illustrated in the battery model BM shown in FIG. 6 behaves in terms of noise in exactly the same way as a parameter, it can be compared together with the parameter noise. When examining the process noise v1 of the first state U.sub.D, it can be derived as follows:

(61) $\begin{array}{cc}{\stackrel{.}{U}}_D=-\frac{1}{R_D{C}_D}\left(U_D\left(t\right)+{\sigma }_{U_D}\right)+\frac{1}{C_D}\left(I_{KL}\left(t\right)+{\sigma }_{I_{\mathrm{KL}}}\right)v_1=-\frac{1}{R_D{C}_D}{\sigma }_{U_D}+\frac{1}{C_D}{\sigma }_{I_{KL}}& \left(8\right)\end{array}$

(62) For examining the second process noise component v2 of the second battery cell state SOC, the derivation produces the following:

(63) $\begin{array}{cc}S\stackrel{.}{O}C=\frac{1}{C_n}\left(I_{KL}\left(l\right)+{\sigma }_{I_{KL}}\right)v_2=\frac{1}{C_n}{\sigma }_{I_{KL}}& \left(9\right)\end{array}$

(64) When modelling the process noise v, the third state R.sub.dc1s can also be observed because its behaviour is the same as the parameters. The process noise of the battery cell parameters BZP establishes how greatly the individual battery cell parameters BZP must change. In the case of conventional arrangements, the parameter estimator 2B is used only in order to estimate fixed parameters or battery cell parameters BZP which no longer change. A changing battery cell parameter BZP can result in instability of the dual Kalman filter 2.

(65) The process noise component v.sub.BZP occurs in the state space representation of the battery cell parameters BZP. Within the scope of basic considerations with regard to the Kalman filter 2, this noise value means nothing other than a specific probability of change in the corresponding battery cell parameter BZP. If a battery cell parameter BZP can change very greatly over wide ranges, this corresponds to a high variance σ of the relevant battery cell parameter BZP, whereas, conversely, a smaller change corresponds to a very small variance.

(66) FIGS. 9A to 9D and 10A to 10D illustrate a dynamic dependency of the individual battery cell parameters BZP upon the temperature T, the state of charge SOC and the loading. If, from this, the variance σ of the corresponding battery cell parameter BZP is deduced in relation to the respective influence variable or battery cell characteristic variable BZK, a noise term of the battery cell parameter BZP can be determined. In this manner, it is possible e.g. in the case of the third battery cell parameter R.sub.dc1s and the two battery cell parameters R.sub.D, C.sub.D to observe the state of charge dependency.

(67) FIG. 9A shows the dependency of the internal resistance R.sub.dc1s upon the state of charge SOC as the battery cell state BZZ for various temperatures T.

(68) FIG. 9B shows the dependency of a time constant τ.sub.CD (τ=R.sub.D×C.sub.D) upon the state of charge SOC for various temperatures.

(69) FIG. 9C shows the dependency of the diffusion resistance R.sub.D upon the state of charge SOC for various temperatures T. Furthermore, FIG. 9D shows the dependency of the diffusion capacitance C.sub.D upon the state of charge SOC for various temperatures T.

(70) The capacitance C.sub.n likewise forms a battery cell parameter BZP. FIGS. 10A to 10D shows the normalised discharge capacity for various typical lithium-ion cell chemistries. FIG. 10A shows a normalised discharge capacity for an NMC battery.

(71) FIG. 10B shows a normalised discharge capacity for an LFP battery for various temperatures T.

(72) FIG. 10C shows the normalised discharge capacity in relation to the C-rate for various temperatures T. FIG. 10D shows the normalised discharge capacity over the C-rate for various temperatures T for an LTO battery. With the battery cell parameter C.sub.n, its change in relation to the terminal current I.sub.KL or with regard to the C-rate is observed in FIGS. 10A to 10D. It can assumed that the temperature T changes only relatively slowly in comparison with the other variables and therefore constitutes only a fraction of the noise covariance.

(73) FIGS. 11A to 11D illustrate a variance of the change in the parameters and a change in state in relation to the state of charge SOC. FIG. 11A shows the variance of the internal resistance R.sub.dc,1s in terms of the state of charge SOC. FIG. 11A shows the variance of a time constant in terms of the state of charge SOC. FIG. 11C shows the variance of the diffusion resistance R.sub.D in terms of the state of charge SOC. FIG. 11D shows the variance of a diffusion capacitance C.sub.D in terms of the state of charge SOC.

(74) Since the variance, which is dependent upon the state of charge SOC, is present, the difference in the state of charge SOC for the cycle period t.sub.2−t.sub.2, i.e. (SOC(t.sub.2)−SOC(t.sub.1)) can be used:
v.sub.3=((SOC(t.sub.2)−SOC(t.sub.1))σR.sub.DC,1s
r.sub.2=((SOC(t.sub.2)−SOC(t.sub.1))σR.sub.D
r.sub.3=((SOC(t.sub.2)−SOC(t.sub.1))σC.sub.D
V.sub.3.sub.k=(X.sub.2.sub.k−x.sub.2.sub.k-1)σR.sub.DC,1s
r.sub.2.sub.k=(x.sub.2.sub.k−x.sub.2.sub.k-1)σ.sub.R.sub.D
r.sub.3.sub.k=(x.sub.2.sub.k−x.sub.2.sub.k-1)σ.sub.C.sub.D    (10)

(75) As already shown in FIGS. 10A to 10D, the battery capacitance is dependent on various factors, in particular upon the terminal current, the temperature, the discharge depth, aging and recovery. On the basis of the various dependencies, it is possible to implement an adaptive method on the basis of modelling of the process noise component of the cell capacitance, the adaptation effort being minimised with said method.

(76) In the case of this modelling, the dependency of the dynamically changing terminal current I.sub.KL can be taken into consideration, wherein it is assumed that the temperature T changes comparatively slowly and therefore makes a subordinated contribution to the variance σ. Since at lower temperatures T the capacitance of the battery cell BZ drops considerably, the cell behaviour at a temperature T=0° C. is used for the variance examination.

(77) FIGS. 12A to 12D illustrate various histograms for various typical lithium-ion cell chemistries.

(78) FIGS. 12A to 12D show a variance of the change of cell capacitance with respect to the terminal current I.sub.KL in four characteristic lithium-ion cell chemistries.

(79) FIG. 12A shows a capacitance variance for an NMC battery. FIG. 12B shows a capacitance variance in the case of an LFP battery. FIG. 12C shows a capacitance variance in the case of an NCA battery and FIG. 12D shows a capacitance variance in the case of an LTO battery.

(80) Since the variance is present in dependence upon the terminal current I.sub.KL as a multiple of the original cell capacitance after the production C0, the difference in the terminal current can be used for the cycle period:

(81) $\begin{array}{cc}{r}_1=\left(\frac{I_{\mathrm{KL}}\left(t_2\right)}{C_0}-\frac{I_{\mathrm{KL}}\left(t_1\right)}{C_0}\right){\sigma }_{C_n}& \left(11\right)\end{array}$

(82) The discrete noise modelling produces the following:

(83) $\begin{array}{cc}{r}_{1_k}=\left(\frac{U_k}{C_0}-\frac{U_{k-1}}{C_0}\right){\sigma }_{C_n}& \left(12\right)\end{array}$

(84) FIGS. 10A to 10D show that the two cell chemistries NCA and NMC hardly differ and the LFP and LTO chemistry includes strong fluctuations. An LFP battery cell displays the slightest noise behaviour. This means that the cell capacitance of the battery cell BZ changes much less when using this cell chemistry than in the case of other cell chemistries. LTO battery cell chemistry lies therebetween. For the adaptation of the filter to the respective application the cell chemistry can be selected and the noise behaviour thereof can be considered accordingly.

(85) In the case of a preferred embodiment of the device 1 in accordance with the invention and of the method in accordance with the invention, stability monitoring of the dual Kalman filter 2 additionally takes place. This preferably takes place on the basis of a covariance behaviour of the battery cell states BZZ and of a covariance behaviour of the battery cell parameters BZP. This stability analysis is based on the fact that, in the case of the threat of instability, the covariance of the corresponding battery cell state BZZ or battery cell parameter BZP increases although this would actually have to decrease when the filter excitation is present. This happens long before the respective value of the battery cell state BZZ or battery cell parameter BZP assumes a non-plausible value.

(86) In order to use the behaviour of the covariance of each battery cell state BZZ and of each battery cell parameter BZP to deduce whether an erroneous adaptation or an error event is possibly present, which could lead to instability of the dual Kalman filter 2, the correct behaviour of the covariance is examined. In so doing, the diagonal entries of a covariance matrix of the battery cell states BZZ and of the battery cell parameters BZP are preferably used since these represent the autocorrelation of the estimation error of each battery cell state BZZ or battery cell parameter BZP. The individual entries represent a measurement of the quality of the present estimation. The smaller the entries or diagonal entries of the covariance matrix, the more precise the present estimation of the dual Kalman filter 2. In order to evaluate the entries, a reference variable is used. Its absolute value does not represent a statement relating to the correct mode of function of the dual Kalman filter 2. On the contrary, the change in this reference variable in comparison with the preceding values is a measure of the correct mode of operation, i.e. if the covariance falls, the error in the estimation becomes smaller, whereas, if the covariance rises, it is highly probable that the error in the estimation will also increase.

(87) In order to examine or analyse the general behaviour of the covariance, its progression is preferably firstly examined based upon a linear Kalman filter. For this purpose, constant values can be assumed for the matrices A, B, C which are variable over time and for the noise terms. The covariance P also does not depend upon the inputs of the dual Kalman filter 2, whereby the progression thereof can be performed independently of the input variables.

(88) FIG. 13 illustrates the progression of the two diagonal values of the covariance matrix P. FIG. 13 shows theoretical progression of the covariance of the two battery cell states U.sub.d, SOC, i.e. of the diffusion voltage and of the state of charge of the battery cell BZ. FIG. 13 shows the progression over time of a covariance for these two battery cell states BZZ.

(89) As can be seen in FIG. 13, the two covariance values decrease exponentially and approximate asymptotically an end value. This end value depends upon the noise components n, v.

(90) An increase in the measurement noise n results in a slight increase in the limit value, i.e. a flat progression. An increase in the process noise v results in a direct increase in the limit value and in a slight change in the progression.

(91) Based upon these examinations, for each battery cell state BZZ and for each battery cell parameter BZP, a funnel which is adapted to the covariance behaviour is defined in exponential form. If dynamics are present in the system, the respective covariance P must be located within the specified funnel. As soon as the covariance P departs from the funnel or exceeds the value of the funnel function in the event of sufficient excitation, a threat of instability of the dual Kalman filter 2 can be detected. Then, the adaptation of the respective battery cell state BZZ or battery cell parameter BZP is preferably suspended by the dual Kalman filter 2.

(92) The progression of the covariance P within the funnel function TF is illustrated e.g. in FIG. 14. FIG. 14 shows the progression of the diffusion voltage U.sub.D as the battery cell state BZZ within or below the funnel function TF(t). It is necessary to ensure that there is a sufficient distance between the progression of the funnel function and the correct covariance progression. In the case of a preferred embodiment, the funnel function TF is a funnel function which decreases exponentially and has configurable coefficients a, b, c.

(93) In the case of one possible embodiment, the funnel function TF(t) is in the following form:
TF(t)=a.Math.e.sup.b.Math.t+c,
wherein the coefficients a, b, determine the pitch of the curve and the coefficient c determines the final limit value of the function, wherein b<0, a>0 and c>0. Since the dual Kalman filter 2 is a non-linear dynamic system in which even the individual noise terms depend upon specific influence factors, the form of the funnel or the funnel function TF is adapted preferably on the basis of the real reaction of the individual battery cell states BZZ and battery cell parameters BZP.

(94) The method in accordance with the invention and the device 1 in accordance with the invention render it possible to estimate battery cell parameters BZP of the battery cell BZ in a dynamic manner with the battery cell states BZZ. Furthermore, the method in accordance with the invention allows the present battery capacitance to be estimated and to also be taken into consideration in the calculation of the state of charge SOC. Furthermore, the method in accordance with the invention allows the state of charge SOC to be estimated on the basis of the present battery capacitance and not, as was previously incorrectly the case, via the rated capacitance C0 which is indicated in the data sheet of the battery cell. Furthermore, the method in accordance with the invention and the device 1 in accordance with the invention render it possible to specify the noise behaviour on the basis of the characteristic parameter influences of the respective cell chemistry. In this case, it does not have to be readjusted by filter-tuning. Furthermore, it is possible also to take into consideration only the respective cell chemistry, wherein the measurement does not have to be performed separately for different battery cells.

(95) Furthermore, the stability monitoring method in accordance with the invention renders it possible to recognise a threat of instability of each battery cell state BZZ and/or each battery cell parameter BZP early. Therefore, it is possible to suspend the estimation of the respective battery cell state BZZ and battery cell parameter BZP until the detected instability is rectified. This can ensure a stable operation of the dual Kalman filter 2. Furthermore, it is possible also to use the dual Kalman filter 2 as an adaptive filter in battery management systems BMS.

(96) In the case of the method in accordance with the invention, the filter-tuning of the dual Kalman filter 2 is effected with the aid of measurable noise components. This allows a rapid adaptation of the dual Kalman filter 2 to different new battery systems or battery cells BZ. The respective influence variables can be controlled in an exact manner.

(97) The method in accordance with the invention allows the state of charge of a battery cell BZ or of a battery to be determined with very high precision. If this battery cell BZ is used e.g. for driving an electric motor of a vehicle, the range prediction for the vehicle is significantly improved thereby.

(98) The detection method in accordance with the invention operates robustly and reliably. It can be flexibly adapted for different systems of battery cell types. The method in accordance with the invention can be implemented with the aid of a dual Kalman filter 2 with a relatively small amount of technical outlay. In the case of one possible embodiment variant, the detection device 1 in accordance with the invention can be integrated in the housing of a battery which includes one or a plurality of battery cells BZ, in order to detect battery cell states BZZ and/or battery cell parameters BZP of the battery cell BZ. The detected battery cell states BZZ and/or detected battery cell parameters BZP can be communicated to a control unit. In dependence upon the detected battery cell parameters BZP and detected battery cell states BZZ, this control unit can activate a current source for charging the battery cell and/or an electric load operated by the battery cell. The detection method in accordance with the invention can be carried out in real time. The battery cell states BZZ and/or battery cell parameters BZP detected in real time can be evaluated in real time and used for closed-loop and open-loop processes e.g. in an electric motor or a photovoltaic installation. In a preferred embodiment, the detection device 1 in accordance with the invention has integrated stability monitoring and is automatically deactivated when there is a threat of instability.

(99) According to a further aspect, the invention provides a battery management system BMS for a battery which consists of one or a plurality of battery cells, having a detection device 1 in accordance with the invention for detection of battery cell states BZZ and/or battery cell parameters BZP.

(100) According to a further aspect, the invention provides an electric vehicle with such a battery management system and one or a plurality of battery cells BZ.

(101) According to a further aspect, the invention provides a photovoltaic installation with rechargeable battery cells BZ and a battery management system BMS, said installation having a detection device 1 in accordance with the invention.

(102) The method in accordance with the invention and the device 1 in accordance with the invention for detection of battery cell states BZZ and/or battery cell parameters BZP can be used for multiple purposes and can be used with any rechargeable battery cells BZ or energy stores.