1. Introduction
With the increasing popularity of electric/hybrid transportation vehicles, the storage and utilization of energy is an important topic, especially regarding secondary batteries. Among them, the lithium-ion battery has emerged as a leading mobile energy storage solution for its remarkable characteristics, including high energy density, elevated operating voltage, eco-friendliness, and absence of memory effect [
1]. Therefore, more attention and research interest are focused on lithium-ion batteries [
2]. However, Nickel–cadmium batteries have also been widely used in direct-current power supply, distribution substation, and auxiliary power supply of rail transit vehicles for the advantages of low cost, high safety, wide temperature range, high reliability, and high mechanical strength [
3,
4]. Especially in rail transit, vehicle safety is always the first priority; however, there is still a lack of a Battery Management System (BMS) due to the inability to accurately estimate the state-of-charge (SOC) of the battery, because SOC estimation is the basis of BMS system development. Accurate SOC estimation can effectively prevent battery overcharging and over discharging, improve the reliability and safety, and extend battery life [
5,
6,
7]. With an increasingly urgent demand for intelligent and automated rail transit, it is necessary to carry out SOC estimation for rail transit vehicle nickel–cadmium batteries.
The existing SOC estimation methods can be divided into two categories: electrical testing and intelligent algorithm. The former includes the discharging experiment method, ampere-hour integration method, open-circuit voltage (OCV) method, internal resistance method, etc. [
8,
9]. The latter includes fuzzy logic, neural networks [
10], and the Kalman filter method [
11]. Among them, the discharging experiment method takes a long time and is an offline estimation. The main disadvantage of ampere-hour integration method is that it relies too much on an accurate SOC initial value, and the error will gradually increase in the integration process [
8]. The open-circuit voltage method assumes that there is a specific functional mapping relationship between the OCV and the SOC of the battery. However, the OCV measurement of the battery usually requires a long static time, which is not conducive to online estimation. Usually, the OCV method is chosen to combine with other estimation methods to estimate the SOC. The internal resistance method obtains the functional mapping relationship between the internal resistance and SOC by measuring the internal resistance of the battery. However, the measurement of the internal resistance of the battery itself has a large error, and the internal resistance is also affected by environmental variables, such as temperature, so the application scope of this method is relatively narrow [
9]. Fuzzy logic algorithms and neural network algorithms have high accuracy in SOC estimation, but they require a large amount of computation and rely on a large number of sample data [
10]. The Kalman filter (KF) has the advantages of high precision and strong anti-noise interference [
11]. In summary, KF is currently one of the most effective methods for battery SOC estimation and has received the most attention and applications. However, it relies on accurate battery models.
The existing battery models can be mainly divided into three categories according to different research focuses: the electrochemical model [
12], thermal model (including electrothermal coupling model) [
13], and performance model [
14]. The electrochemical model is a battery model based on the electrochemical theory of ion diffusion and polarization effect, and the partial differential equation was used to describe the changes of electrode surface concentration and electrode overpotential inside the battery. However, it is difficult to obtain certain parameters, such as solid phase concentration and solid phase potential. Moreover, the solution of partial differential equations, which are mostly used in the design and development of high-performance batteries, is complex. Based on the conservation of energy, the thermal model is a mathematical model established based on the three thermal processes of heat generation, heat conduction, and heat loss. However, it is difficult to obtain the parameters of reversible heat and irreversible heat in the thermal model, so it is rarely used in SOC estimation, and is instead often used in the design of battery thermal management system and the study of single-battery heat production. The performance model is widely used in the study of external characteristics of batteries. The two mainstream performance models are the black box model and the equivalent circuit model [
15]. Black-box models, including neural networks, fuzzy logic, and support vector machines, rely on a large number of sample data and have a large amount of computation. BMS system chips are difficult to meet their computational requirements. The equivalent circuit model simulates the external characteristics of a battery based on the experimental data of battery charging and discharging, using electrical components such as capacitors, resistors, and constant voltage sources. It has the advantages of a simple model, the easy identification of parameters, low computational complexity, and high accuracy, and is widely used in SOC estimation. It should be emphasized that the existing battery equivalent circuit model is mainly for lithium batteries, while the main reactions and side reactions and working characteristics of nickel–cadmium batteries are significantly different from those of lithium batteries, so there is a lack of mature models for reference at present.
Therefore, by analyzing the basic characteristics of nickel–cadmium batteries, an equivalent circuit model of nickel–cadmium battery was established, the model parameters were identified, and the accuracy of the model was verified by experimental data and simulation. Then, based on the constructed equivalent circuit model, a SOC estimation method of Ni-Cd batteries based on adaptive unscented Kalman filter (AUKF) is proposed, and the results are compared with unscented Kalman filter (UKF) under different operating conditions.
3. AUKF Estimation Algorithm
In
Section 2, the equivalent circuit model of the nickel–cadmium battery was established. Based on the model and KF algorithm, the SOC of the nickel–cadmium battery can be estimated. The standard KF algorithm is a pure time-domain filter, which has the advantages of a simple structure and good robustness. However, it is only suitable for the modeling analysis of linear systems. The nickel–cadmium battery has polarization effect and diffusion effect in its internal and nonlinear external performance, so it is necessary to use the improved Kalman filter algorithm. At present, the commonly improved forms are the extended Kalman filter (EKF) and unscented Kalman filter (UKF) [
23]. When the EKF algorithm performs Taylor expansion on the nonlinear part, it will produce a Taylor truncation error and may make the filter difficult to converge and increase the error. Moreover, the EKF algorithm needs to calculate the Jacobian matrix once in each iteration process, which increases the amount of calculations. The UKF algorithm [
24,
25,
26] can also be used for the state estimation of nonlinear systems. Its core method is unscented transformation (UT) [
27,
28], in which sampling points are used to approximate the distribution of nonlinear state variables, without the linearization of nonlinear functions and Jacobian matrix solution. Compared to the EKF algorithm, the UKF algorithm has a simpler principle, higher accuracy, less computation, and wider application.
3.1. UKF Algorithm
The discrete state space model of linear system is generally composed of state equation and observation equation, and the general expression is shown in Equation (12).
where
xk is the state quantity at time
k;
xk+1 is the state quantity at the next time;
uk is the system input quantity;
yk is the system observation value at time
k; and
A,
B,
C, and
D are the state transition matrix, input matrix, output matrix, and feedforward matrix of the system, respectively;
wk and
vk refer to process noise and observation noise, respectively, which are generally set as Gaussian white noise, as shown in Equation (3).
where
Qk(
Rk) is the covariance matrix of process noise (observation noise).
The UT in the UKF algorithm [
27,
28] is used to sample near the state points of nonlinear functions with reference to certain sampling strategies. The sampling points are collectively referred to as the sigma point set. The sigma point set and the state variables of the original state distribution have the same mean and variance. The sigma point set is substituted into the original nonlinear function relationship for nonlinear transfer, and the nonlinear function value point set is calculated. The mean and covariance of the transformed point set are further calculated. The UKF algorithm can approximate the mean and covariance to the third-order Taylor series. For the n-dimensional state variable,
x,
, and P are used to represent the mean and variance of the state variable, respectively. Firstly, 2
n + 1 sampling points are obtained by using the symmetrical sampling strategy, and then the corresponding weights of each sampling point are calculated. Finally, the mean and covariance of the output variable are obtained. The specific steps are as follows:
(1) Calculate 2
n + 1 sampling points according to Equation (14).
where
n refers to the dimension of the state variable,
xi refers to sigma point
i,
refers to the number
i column of the root mean square of the
matrix,
refers to the scale adjustment factor, and the distance between the sampling point and the mean
x can be adjusted by changing
so as to adjust the accuracy of the algorithm.
(2) Calculate the weight corresponding to each sampling point.
In Equation (15), is the mean weight; represents covariance weight; is the dispersion degree factor, which is used to describe the dispersion of sampling points around the mean and meets the relationship ; k is the auxiliary scale factor, usually taken as 0 or ; and is the prior distribution factor, which is used to suppress the error caused by higher-order terms. For Gaussian distribution, the best result can be obtained when is taken.
(3) Calculate the mean
and covariance
of the output variable.
Based on UT, the detailed steps of UKF algorithm can be obtained as follows:
(1) Initialize the calculation of state variable
and covariance
.
(2) Combined with UT, calculate the sigma point sets using the symmetrical sampling strategy.
(3) Substitute the sigma point sets into the state equation to complete the backward propagation of the sigma point sets.
(4) According to the weight corresponding to the sigma point sets, complete the prediction of the state quantity and the calculation of the covariance matrix.
In Equation (21), is the average value of process noise.
(5) Substitute the sigma point sets into the observation equation to calculate the system observation value.
(6) According to the weight corresponding to sigma point set, complete the prediction of observation and the calculation of covariance and cross covariance matrix.
In Equation (23), is the average value of observation noise.
(7) Obtainthe Kalman gain matrix from the covariance matrix.
(8) Update the system state matrix and its covariance matrix
3.2. AUKF Algorithm
When using UKF algorithm to deal with battery SOC estimation, it is necessary to obtain the prior system noise and observation noise information, which is set as a constant value by the standard UKF algorithm. However, there will be a certain deviation between the set constant value and the actual noise value in each iteration process, and the error will gradually accumulate with the increase in iteration times, which will reduce the estimation accuracy of the UKF algorithm, and may even lead to filtering divergence. Therefore, in order to further improve the accuracy of SOC estimation, the Sage–Husa adaptive filtering algorithm [
29,
30,
31,
32] is selected to correct the covariance of process noise and observation noise in real time. The main steps of the standard Sage-Husa adaptive filtering algorithm for noise adaptation are as follows:
(1) Calculate the estimated average value of process noise.
In Equation (26), is the forgetting coefficient, and is the forgetting factor, which generally meets the relationship .
(2) Calculate the covariance estimated value of process noise.
In Equation (27), represents the difference between the true value of the observation value and the estimated value of the observation value at time , which is the observation residual, also known as the innovation sequence of the observation value.
(3) Calculate the estimated average value of the observation noise
(4) Calculate the estimated covariance value of the system observation noise
Under some special conditions (such as large current fluctuations), it is difficult to ensure that the process noise covariance matrix
and the observation noise covariance matrix
maintain non-negative definiteness. If
and
are negative definite, the covariance matrix
of the state variable will be non-semi-positive definite, resulting in the failure of the sigma sampling point calculation in the UKF algorithm and the failure of the program [
33]. Therefore,
and
are improved as follows:
In Equation (30), means to take the main diagonal elements of the original matrix to form a diagonal matrix, and Equations (26)–(30) constitute an improved Sage–Husa adaptive filtering algorithm. When combined with the UKF algorithm, it can effectively ensure the semi-positive definiteness of the state variable covariance matrix so that the UKF algorithm can continue to iterate and improve the stability and persistence of filtering.
The improved Sage–Husa adaptive filtering algorithm is used to adaptively modify the noise covariance matrix in the UKF algorithm, and
,
,
, and
in the UKF algorithm are replaced by
,
,
, and
, respectively, to form the AUKF algorithm. The specific flow of the AUKF algorithm is shown in
Figure 12.
4. Simulation and Experimental Verification
Before using the UKF algorithm and the AUKF algorithm to estimate the SOC, the state equation and observation equation of the model need to be discretized. Based on the equivalent circuit model shown in
Figure 5, the battery SOC,
branch voltage
and
branch voltage
are selected as the three state variables of the system, the charging and discharging current is selected as the system input value, and the battery terminal voltage is selected as the system observation value, and the battery state equation is derived, as shown in Equation (31).
The battery observation equation is shown in Equation (32).
The above equation is discretized, and the discretized equation of state is obtained, as shown in Equation (33).
The discretized output equation is shown in Equation (34).
The nonlinear form of the model can be further written as
The above nonlinear function is substituted into the AUKF algorithm flow, and the SOC estimation of nickel–cadmium battery is completed according to Equations (17)–(30).
4.1. SOC Estimation Verification under HPPC Test
Due to the typical voltage plateau period of nickel–cadmium battery, the voltage change during the plateau period is small, so a small voltage error may produce a large SOC error. In addition, in the actual operation of the battery, the initial value of the SOC is usually unknown, and it will spend more time to obtain the initial value of the SOC by standing for a long time, so the SOC estimation algorithm needs to be robust and can accept an initial value error within a certain range.
The initial values of covariance matrix, process noise, and observation noise are set as follows:
Two parameters, MAE and RMSE, are selected to comprehensively evaluate the performance of the two algorithms. The comparison results are shown in
Table 4.
HPPC test steps are shown in
Section 2.2.1. Before the experiment, the battery was in a fully charged and static state for a long time, and the polarization effect in the battery was not obvious. It can be considered that the initial value of RC network terminal voltage was 0, and that the initial value of the SOC was 1. When the initial value of the SOC was set to 0.7, there was an initial error of 30% with the actual initial value. The robustness and accuracy of the algorithm were tested.
Figure 13 shows the estimation results of two algorithms for SOC estimation under HPPC conditions.
Based on the analysis in
Figure 13 and
Table 4, the UKF algorithm and AUKF algorithm can still converge to the real value of the SOC when there is a large error in the initial value of the SOC, which has good robustness. As can be seen from
Figure 13b,d, the noise error of UKF algorithm will gradually accumulate with the increase in iteration times, which shows that the error increases at the end of the working condition, and the maximum error reaches 2.93%. The AUKF algorithm updates the noise adaptively in each iteration process, and the error is smaller and more stable. Under the three initial values, MAE and RMSE are stable within 0.5%, and the accuracy is higher. Simulation results show that the AUKF algorithm has higher accuracy and robustness in estimating the SOC of the nickel–cadmium battery.
4.2. Verification of SOC Estimation under Constant Power Discharging Cycle
In the actual use of nickel–cadmium batteries, the discharging current is usually dynamic and in a state of recycling. Therefore, the constant power discharging cycle is further selected to verify the model. The discharging process is set at a constant power of 66.5 W, the charging process is set at 0.2 C constant current, the charging and discharging capacity is 20 Ah, and the cycle repeats five times. The voltage and current changes under this working condition are shown in
Figure 14.
Under this working condition, the current changes continuously during the discharging process, including both the discharging and charging processes, which tests the dynamic following ability of the algorithm. Similarly, the SOC initial value error was set to 30% to further evaluate the robustness of the algorithm under this working condition, as shown in
Figure 15.
The SOC estimation results of the two algorithms under constant power cycle are shown in
Table 5.
Based on
Figure 15 and
Table 5, when the battery changes from discharging to charging, the UKF estimation error becomes larger, while the AUKF can maintain a smaller error throughout the cycle. The MAE and RMSE of the AUKF algorithm in the convergence stage are always kept within 0.8%, which has good accuracy. The simulation results show that the AUKF algorithm can still maintain good dynamic following ability under cyclic conditions, and has high accuracy. Hence, the AUKF is more suitable for the SOC estimation of nickel–cadmium batteries.
5. Conclusions
In this paper, the open-circuit voltage characteristics and polarization voltage characteristics of nickel–cadmium battery were studied. It was found that the intermittent charge–discharge method used to determine the OCV-SOC curve of a nickel–cadmium battery was more applicable than the micro-current method. It showed that the nickel–cadmium battery had typical voltage hysteresis characteristics, as well as the characteristics of ohmic polarization, electrochemical polarization, and concentration polarization during the whole process. All these characteristics can provide a reference for the modeling of nickel–cadmium batteries.
Then, considering the characteristics of open-circuit voltage and polarized voltage and the complexity of the model, the equivalent circuit model of nickel–cadmium battery was constructed, in which the voltage hysteresis, ohmic polarization, electrochemical polarization, and concentration polarization were considered. Based on the HPPC experimental data, the average absolute error of the constructed model was only 4.05 mV, the root mean square error was only 7.40 mV, the maximum absolute error was 49.79 mV, and the relative error was less than 1% and the maximum relative error was less than 5% most of the time.
The SOC estimation method of nickel–cadmium battery based on the established dynamic model and the AUKF algorithm was proposed. According to the estimated results from the measured data, it can be seen that, although the method is insensitive to the initial value of SOC, it can converge to the real values even if the initial value has a 30% error. The average absolute error and root mean square error of the final SOC estimation results are less than 1%.