Abstract

In general, battery packs are monitored by the battery management system (BMS) to ensure the efficiency and reliability of the energy storage system. SOC and SOH represent the battery’s energy and lifetime, respectively. They are the core aspects of the battery BMS. The traditional method assumes that the SOC is determined by the integral of the current input and output from the battery over time, which is an open-loop-based approach and often accompanies by poor estimation accuracy and the accumulation of sensor errors. The contribution of this work is to establish a new equivalent circuit model based on the lithium battery external characteristic, and the battery parameters are identified by considering the influence of capacity fade, voltage rebound, and internal capacitance-resistance performance. The correlation between the ohmic internal resistance and real capacity is obtained by degradation test. Then, the dual extended Kalman filter (DEKF) is used to perform real-time prediction of the lithium battery state. And through the simulation analysis and experiments, the feasibility and precision of the estimation method are well proved.

1. Introduction

Lithium battery has been widely used in the energy storage field due to its high energy density, long cycle life, high voltage, and outstanding security [1]. In general, in order to ensure the efficiency and reliability of the energy storage system, battery packs are monitored by the battery management system (BMS). And the state estimation of the battery is the essential function of the BMS. The accurate online state estimation can contribute to monitor the battery’s state of charge (SOC) and state of health (SOH) clearly so that the reasonable charge-discharge can be performed to keep its high-efficiency working conditions and long life [2].

SOC and SOH represent the battery’s energy and lifetime, respectively. They are the core aspects of the battery BMS. The traditional method assumes that the SOC is determined by the integral of the current input and output from the battery over time. This is an open-loop-based approach, the result of which is often accompanied by poor estimation accuracy and the accumulation of sensor errors [3]. Moreover, in the past, the battery’s capacity variations were not to be taken into account during the lifetime; thus, the estimation method had great deficiencies [4]. In fact, it comes with irreversible chemical reactions and physical changes in the working process of the battery, which make SOC and SOH tightly coupled with each other. Improper battery state estimation method can lead to premature damage and deterioration of the battery [5]. As the impedance and capacity of the battery change with time, the maximum power and energy that can be delivered by the battery will be reduced. This technical problem can be solved by simultaneously estimating the SOC and the SOH of the battery [6].

Traditional estimation methods, including the internal resistance method, coulomb counting (CC) method, open-circuit voltage method, and load charge method, have some limitations in practical application due to lack of abilities to correct errors and largely depending on high precision of sensors and low environment influences [7]. Nowadays, with the development of modern control theory, there are some new algorithms which are constantly being used in battery state estimation. Plett [8] proposed the sigma-point Kalman filtering (SPKF) method as applied to BMS algorithms through the analysis of batteries’ nonlinear characteristics. And a better result was obtained by comparison with the traditional Kalman filtering method. Fleischer et al. [9] presented a state-of-available-power prediction method based on the framework of the adaptive neuro-fuzzy inference system. The method achieved the accurate power prediction at low-temperature environment. Lee et al. [10] had researched the linear prediction error method and the neural network method in the battery state estimation and then found the latter had better precision in short-term prediction. Chen et al. [11] employed the genetic algorithm to estimate the battery model parameters based on analysis of the resistance-capacitance circuit model. The precision and robustness of results are much improved, but its real-time performance was limited.

The contribution of this paper is to establish a new equivalent circuit model based on the lithium battery external characteristic, and the battery parameters are identified by considering the influence of capacity fade, voltage rebound, and internal capacitance-resistance performance. The correlation between the ohmic internal resistance and real capacity is obtained by degradation test. Then, the dual extended Kalman filter (DEKF) is used to perform real-time prediction of the lithium battery state. And through the simulation analysis and experiments, the feasibility and precision of the estimation method are well proved.

The rest of the paper is organized as follows. An equivalent circuit model is established, and the characteristic parameters of the battery are identified in Section 2. In Section 3, the method of battery state estimation is proposed, and the computer-aided modeling analysis is carried out. Section 4 provides detailed experimental procedures and results. Section 5 presents some conclusions.

2. Battery Model and Parameter Identification

2.1. Battery Model

The battery performance is closely related to its internal parameters. And the change rule of characteristic parameters can be analyzed by establishing the reasonable equivalent model [12]. In order to achieve the real-time state analysis, not only the working characteristics, including voltage rebound effect, hysteresis, and electromotive force-SOC (EMF-SOC), need to be taken into consideration in the model but also the complex electrochemistry mechanism should be avoided as much as possible. Therefore, mainly focusing on the analysis of external characteristics, a new circuit containing a controlled source is established as the equivalent circuit model of the lithium battery. The model is based on Thevenin’s theorem and consists of an equivalent voltage source and equivalent impedance, as shown in Figure 1.

2.1.1. The Equivalent Voltage Source

The voltage source is related with the SOC, SOH, and working conditions of the battery. Therefore, the equivalent voltage source consists of electromotive force and hysteresis voltage . As the voltage-controlled voltage source (VCVS), is controlled by voltage . And is the capacity of the battery, which is related with the SOH. And the current flowing through is equal to the working current . Therefore, the value of can be used to represent the battery SOC:where is the working time and is the remaining capacity.

Hysteresis voltage is also a voltage-controlled voltage source (VCVS), which is controlled by voltage . The current flowing through inductance is . Hysteresis voltage is influenced by previous working current, in which the characteristic can be described by a circuit containing the inductance. The value of is determined by charging and discharging conditions:

Therefore, the equivalent voltage source is

2.1.2. The Equivalent Impedance

The capacitance and resistance performance of the battery are described by an equivalent impedance submodel, which consists of the ohmic internal resistance and capacitive-resistance network [13]. is caused by materials and the specific structure of the battery, which will increase with the decay of SOH. And the capacitive-resistance network is used to describe the voltage rebound effect caused by polarization capacitance and polarization resistance [14]. The rebound performance of the model is closer to the real value of the lithium battery with higher the order of the network. Meanwhile, the complexity and computation time of the model will increase with the orders. Considering the accuracy of the model and the computational efficiency, the third-order is much reasonable. Then, according to zero-state response of the network, the polarization voltage can be shown as

Therefore, according to Thevenin's theorem, the complete equivalent circuit model of the lithium battery can be expressed as

Through the equivalent circuit model, the relations of internal parameters and external working characteristic, including EMF-SOC, voltage hysteresis, rebound effect, and impedance characteristic of the lithium battery, can be described clearly.

2.2. Parameter Identification

At the constant temperature of 25°C, a pulse discharge test for the max load 18650 lithium battery is performed to identify the model characteristic parameters. The battery capacity is 2900mAh, and the current is 2.9 A. The test current is shown in Figure 2(a). SOC and voltage during the test are shown in Figure 2(b).

Based on the test data of the balance electromotive force and corresponding SOC, the characteristic curve can be obtained by the polynomial segmentation fitting method, as shown in Figure 3.

The corresponding mathematical expression iswhere is the coefficient matrix, and

The networks are in the zero-state respond during the voltage rebound period. Therefore, according to equation (5), the output voltage can be shown aswhere and are, respectively, the rebound voltages in the first stage and second stage and can be expressed aswhere , , , and .

According to test data and equations (8)∼(10), the voltage rebound characteristic curve can be obtained, as shown in Figure 4. And the parameters of the fitting curve are and . Therefore, the characteristic parameters of the battery can be identified, as shown in Table 1.

3. Battery State Estimation

3.1. Battery State

Battery state contains the state of charge (SOC) and the state of health (SOH). SOC is the remain charge that can be released from the battery, which is related with the initial charge state, actual capacity, and previous current [15]:where is the initial charge state and is the coefficient of the current ratio.

SOH is evaluated based on the battery capacity. According to battery failure standard, the battery is invalid when its actual capacity drops to 80% of rated capacity [16]. Therefore, SOH can be defined aswhere ,, and are, respectively, the rated capacity, actual capacity, and invalid capacity.

The lithium ion in the electrolyte is hindered by the solid electrolyte interface (SEI) adhering to the electrode of the lithium battery. With continuous electrochemical reaction, SEI becomes thicker and thicker, which means hindrance to the lithium ion gets more serious. This can be seen from the increased ohmic internal resistance in the macroscopic view [17]. Therefore, the battery capacity gradually decreases, and ohmic internal resistance increases during the battery degradation [18]. In the normal working condition, battery degradation is a long-term and slow process. In this work, the changes of the ohmic internal resistance and capacity are obtained by accelerating the growth of solid electrolyte membranes in the high-temperature working environment. The degradation curve is shown in Figure 5.

The corresponding mathematical expression is

According to the result analysis, SOC and SOH are closely related. The estimation accuracy of SOC can be improved by considering the influence of capacity degradation.

3.2. Estimation Principle of DEKF

Extended Kalman filter (EKF) is an improved method based on the traditional Kalman filter that can be used for the estimation of nonlinear state variables [19]. And when two sets of state variables need to be estimated, dual extended Kalman filtering (DEKF) can be used [20]. In this work, the battery SOC is estimated, which is related with the polarization resistances. The state variable is defined as . Besides, SOH can be used to calculate the actual capacity to improve the estimation accuracy of SOC. Although SOH cannot be measured directly, it can be obtained by the estimation of the ohmic internal resistance due to equation (13). So, another state variable is defined as . Besides, an observable variable related with state variables is needed in Kalman filter [21]. According to the battery model, the output voltage can be regarded as the observational variable for it can be measured by a sensor directly. The measuring equation at time k is shown aswhere is the characteristic equation and is the measurement noise.

The battery state estimation based on the DEKF method is shown in Figure 6.

According to the result of the last cycle, the battery capacity at time k can be updated. Then, based on the equivalent circuit model, state equations are established, and the preliminary state predictions and are calculated bywhere is the process variable, is the process noise, and is the small perturbation of state . And is the preliminary estimation of state variables.

The estimated values obtained by state equations cannot avoid the accumulation of process errors. Besides, the initial errors that may exist is hard to be recognized and corrected [22]. Therefore, an observable variable needs to be established to correct the previous error. This is the key to the algorithm. In the DEKF method, both of the two state variables need to be corrected by the observational variable:

The error corrections depend on the difference value of observational variable and the Kalman gain . is obtained by the voltage sensor, and is calculated by the measuring equation. The Kalman gain is calculated bywhere are the covariances of measurement noise that rely on the performance of the sensor. Jacobian matrixes of state variables are shown in equation (18). are the covariance matrixes of the state, which are used to reflect the change of state variables, as shown in equation (19).where are the covariances of the process noise.

3.3. Simulation Model and Analysis

The real cycle test is time consuming. And it is more efficient for the algorithm analysis to simulate the battery working conditions with the help of computer simulation technique [23]. Therefore, according to the principle of the battery equivalent circuit, a simulation model is established by using mathematical logic modules and circuit element modules in Matlab Simulink, as shown in Figure 7. In order to ensure the reliability of simulation analysis and reflect the real working characteristics of the battery, the parameters of the simulation model are set up corresponding with Table 1.

The SOC estimation is performed by using the DEKF method in the 1C intermittent discharge working condition as shown in Figure 8. The measurement noise and process noise are simulated by Gaussian white noises. According to the sensor performance, the noise covariances are set as and . To consider the ability of error correction, four situations are simulated, where the initial errors are, respectively, set as 0%, 5%, 10%, and 15%, and the step size is 0.5s. The simulation results are shown in Figure 9.

In Figure 9, the initial errors of DEKF1∼ DEKF4 are, respectively, 0%, 5%, 10%, and 15%. It is seen that, with the discharge process, in the DEKF1 situation, which has no initial error, the estimation is still near the real value. And in terms of DEKF2∼ DEKF4 that contain different initial errors, their estimations are gradually close to real value. Therefore, the results show that DEKF has good convergence and stability and can correct the initial errors in the state estimation of the battery.

4. Experiment

4.1. Experimental Setup

To verify the feasibility and accuracy of the DEKF method in the real SOC estimation, the experiments are performed by the BMS circuit. The circuit structure consists of a lithium battery pack, the battery monitoring circuit (BMC), the battery control unit (BCU), and the upper computer, as shown in Figure 10(a). In the experiment, the battery pack consists of 12 MAX LOAD 18650 cells connected in series. The LTC6804-1 chip is used to monitor the battery data. BCU is a Linduino microcontroller. The upper computer and BCU are communicated by the serial peripheral interface (SPI), as shown in Figure 10(b). The load is a resistor, and the current is controlled by a sliding rheostat.

The SOC estimation experiments are performed by, respectively, using the traditional CC method and DEKF method in the working conditions of the open-circuit state, constant-current discharge state, and intermittent discharge state.

4.2. Results and Discussion
4.2.1. Open-Circuit State

When the initial charge is 92%, the SOC estimation is performed in the open-circuit state. The measuring voltage nearly remains constant due to the zero working current, as shown in Figure 11. When initial error is 0, the SOC estimated by both methods is close to 92% with a little fluctuate, as shown in Figure 12(a). Because the current is zero, CC-SOC and REAL-SOC have the same value at this time, and the two curves are basically coincident. And when 2% initial error exists, the SOC estimated by the DEKF method can converge to real value, but the estimation still contains initial error by using the CC method as shown in Figure 12(b).

4.2.2. Constant-Current Discharge State

The measuring voltage constantly decreases in 1.5 C constant-current discharge state, as shown in Figure 13. It is seen that the gradient is great in the high and low SOC stages, while the gradient is much less in the middle stages, which reflects the characteristic of the lithium battery. When initial error is 0, the SOC estimated by the DEKF method is near the real value, while the value by the CC method gradually deviated from the real value, as shown in Figure 14(a). And when 10% initial error is introduced, the estimation of the CC method still contains nearly whole error, while the error is reduced obviously by the DEKF method, as shown in Figure 14(b).

4.2.3. Intermittent Discharge State

The change of measuring voltage in 1C intermittent discharge state is shown in Figure 15. It is seen that the voltage is in a decreasing tendency, but some rebounds and step changes occur in the intermittent phase due to the rebound characteristic of the battery.

When initial errors are 0 and 10%, the estimation results of SOC are shown in Figure 16. It is seen that the SOC can keep track with real value both in cutoff and discharge stage, while the CC method is serious influenced by errors.

4.2.4. Error Analysis

In the open-circuit state, the accuracy of the CC method completely depends on initial error because of zero working current, while DEKF can correct error within a short time. In working conditions, even if the initial error is 0, the maximum errors of CC, respectively, reach 8% and 12%, while the DEKF method is less than 2%. In terms of situation that 10% initial error is introduced, and the error of the CC method is still very big in all of the experiment processes. But, the DEKF method has good convergence, and the initial error can gradually be decreased by correction of the state equation and observational equation.

The average errors of two methods in different conditions are shown in Table 2. It is seen that whatever the initial errors exist or not, average errors of DEKF are less than 2.7%, which is also much less than the errors of the CC method. Besides, the DEKF has good robustness and can well deal with the influence of initial errors, while the CC method much relies on the stable initial value and the precision of sensors. Therefore, it is proved that the DEKF performs well in the SOC estimation of the lithium battery.

5. Conclusions

(1)A new equivalent circuit model that can describe the working characteristics of the lithium battery is established based on Thevenin’s theorem. The internal parameters, including capacity, the ohmic internal resistance, polarization resistance, and capacitance, are identified by fitting characteristic curve. The relation of equivalent internal resistance and capacity is obtained by battery life experiments.(2)The DEKF principle using in battery state estimation is analyzed based on the relation of SOC and SOH. A simulation model is established to analyze the estimation ability with the influence of different initial errors.(3)The SOC estimation experiments are performed by, respectively, using the traditional CC method and DEKF method in three different working conditions. By comparison with the estimation results of the CC method, the results of DEKF are still near the real values, even if the errors of sensors exist. When 10% initial errors are introduced, the DEKF can correct the errors, and the estimation can fast converge to real values. The average errors are less than 3% in all kinds of working conditions, which verifies the feasibility of the algorithm.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors would like to acknowledge the support of the National Key R&D Program of China (Grant no. 2016YFC0802900) and the Fundamental Research Funds for the Central Universities.