1. Introduction
Energy storage systems are very popular in the industry due to their increased demand in a number of applications, from smart micro grids to electrified transportation systems [
1]. Among these energy storage systems, batteries receive the most attention because of their wide range of power and energy densities and low cost as compared to other systems such as flywheels, ultra-capacitors, or hydrogen storage [
2]. There are different types of batteries available in the market like Li Ion, lead acid, and nickel cadmium, etc. Among lead acid batteries, Flooded Lead Acid (FLA) batteries are a versatile and dependable choice in numerous applications, including backup power systems, automotives, renewable energy storage, and aviation [
3]. Their enduring appeal lies in their cost-effectiveness, reliability, and adaptability, making them a reliable energy storage solution across a broad spectrum of industries.
The State of Charge (SOC) estimation of a battery is one of the fundamental requirements of a battery management system. SOC estimation becomes a critical concern in extending battery lifetime, preventing progressive unalterable damage, and enhancing battery operation performance; therefore, a proper analysis of the internal behavior of Sealed Flooded Lead Acid (SFLA) batteries for SOC estimation is required, as this has not been conducted before. A battery management system (BMS) is an application which uses a suitable equivalent model for the battery to protect it from over-voltages and over-currents and estimates the other required states of the battery like the SOC, State of Health (SOH), and State of Function (SOF), etc. [
4,
5,
6]. As the Open-Circuit Voltage (OCV) is computed through the SOC input and SOC–OCV relationship, then combining it with the battery polarization voltage can be used to model the battery dynamics. The SOC–OCV relationship is very important in the modeling of the battery equivalent circuit, since equivalent circuit models make the complex electrochemical processes occurring within a battery simpler by representing them with electrical components like resistors, capacitors, and voltage sources. The SOC–OCV relationship is key in determining the voltage source within this model. By obtaining an accurate SOC–OCV relationship, the equivalent circuit model can effectively mimic how the battery voltage changes during charging and discharging cycles. The relation of SOC–OCV is not simply a one-to-one mapping, but it exhibits the hysteresis effect (the OCV is different at the same SOC value depending upon the history of charge and discharge). Normally, the hysteresis effect is not very effective in most Li-Ion batteries like lithium manganese oxide and lithium cobalt, etc. [
4]. This is because the gap between the SOC–OCV charge and discharge curves is very small in most of the region. But the hysteresis effect is very prominent in lead acid batteries, especially in SFLA batteries [
7]. The low nominal cell voltage of lead acid batteries, changes in charging rates, and the effects of aging contribute to the prominence of the hysteresis effect. Understanding and addressing this effect is crucial for accurate SOC estimation and applicable battery management in applications using lead acid batteries. Therefore, accurately modeling the hysteresis phenomenon in FLA batteries is essential for achieving accurate SOC estimation results.
To model the hysteresis effect of batteries, different approaches have already been published. For those batteries in which the hysteresis effect is not very dominant, an average of the charge and discharge curves is used as the SOC–OCV relationship. This type of approximation will create a large SOC estimation error for those batteries which have a significant hysteresis phenomenon. Sirinivasan et al. proposed an empirical approach to model hysteretic curves based on the boundary curves of hysteresis in nickel hydroxide electrode, and the process worked well in predicting how the battery discharged, but it was not as accurate when predicting how it charged [
8]. Verbrugge et al. formulated a first-order differential equation to construct voltage behavior for NiMH batteries [
9]. Roscher et al. modeled the hysteresis of a battery by using a parallelogram [
10]. Thele et al. formulated second-order polynomial equations to model the shapes of the hysteresis of NiMH batteries [
11]. Similarly, the Takacs model and Preisach model of hysteresis were used in batteries by Windarko et al. and Tang et al., respectively [
12,
13]. All these hysteresis models are based on the fact that the hysteresis loop always touches the SOC–OCV charge curve whenever it is applied from the SOC–OCV discharge curve and vice versa. The hysteresis phenomenon in maintenance-free SFLA batteries has not been investigated before; therefore, the behavior of the OCV requires both proper investigation and modeling according to its behavior. The proposed research is unique in the sense that it fully investigates the OCV behavior and proposes a novel hysteresis model for SFLA batteries.
An online estimation of the parameters of an equivalent model of battery is very important for making the battery model more and more accurate. In this research, an Auto Regression for external Input (ARX) model for the battery combined with the Recursive Least Squares (RLS) method is used for the parameters’ estimations. The ARX model is a type of linear system model used to describe the relationship between the input and output of a dynamic system. It is often used for modeling and control applications. In the context of battery modeling, the ARX model represents the behavior of a battery as a linear combination of past inputs and outputs. The model assumes that the current output is a linear function of previous inputs and outputs, as well as external input signals. Various methods can be used for parameter estimation, and one of the common techniques is the RLS method. RLS is an adaptive filter algorithm used for estimating the parameters of a linear model. It is particularly useful when dealing with time-varying systems and noisy data. RLS updates the parameter estimates iteratively as new data become available. RLS minimizes the sum of squared errors between the observed output and the predicted output of the model. The algorithm recursively updates the parameter estimation by giving more weight to recent data and less weight to older data. The OCV information from the SOC estimation algorithm is sent to the ARX model, which estimates the values of the parameters by the RLS method and provides the parameters back to the SOC estimation algorithm. So, both the parameter estimation algorithm and the SOC estimation algorithm work in parallel and exchange information at every iteration [
14].
For the SOC estimation of a battery, different estimation algorithms are proposed, such as Artificial Neural Networks (ANNs) that control complex neural structures for accurate predictions, ampere-hour counting which integrates the current flow to estimate the SOC, Kalman filters that balance mathematical models and measurements for real-time SOC tracking, and fuzzy logics, which provide a robust, rule-based approach to SOC estimation, etc. Among them, Kalman filters have received the most attention because of their simplicity and accuracy. Kalman filters are usually used for linear systems. For non-linear systems, the Extended Kalman Filter (EKF) was developed, which uses Tylor Series Expansion to linearize any system across a certain point. As batteries are also very non-linear in their behavior, EKFs are widely adopted by most researchers for their SOC estimation. Recently, the Unscented Kalman Filter (UKF) was developed to further improve the estimation accuracy when estimating the state of a non-linear system, surpassing the capabilities of the EKF. While the EKF relies on linear approximations for non-linear systems, which can introduce errors, the UKF takes a more robust approach. It employs a deterministic sampling technique that captures a set of representative points, or sigma points, from the probability distribution of the state variables of the system [
15,
16,
17,
18,
19,
20,
21].
In a Kalman filter, process and measurement noise covariances play a vital role in the calculation of filter gain. Unfortunately, there is not any rule of thumb available for finding the values of these noises for any given system. Mostly, these noises are selected based on a trial-and-error method. In practical systems, especially in electric vehicles (EVs), sudden acceleration and deceleration can cause rapid shifts in the dynamics of the battery system. Maintaining constant noise values may result in higher estimation errors, highlighting the need for adaptive approaches that account for these dynamic changes in the system. Similarly, when constant values are used for the covariance of process and measurement noises in a Kalman filter, the Kalman gain, which determines the weight given to the measurement data, becomes constrained to a specific range. This constraint restricts the filter adaptability to varying noise conditions, limiting its ability to respond effectively to unexpected changes in the environment or system dynamics. In such cases, the filter performance may suffer when faced with real-world uncertainties, as it cannot freely adjust the Kalman gain to optimize the state estimation. Therefore, an adaptive system is required which can update the values of the process and measurement noise covariances based on the system error generated in each iteration, giving independence to the Kalman filter to work freely and more accurately in the SOC estimation process [
22,
23].
In this research, firstly, the most important hysteresis effect of the FLA battery is discussed in detail and, depending upon the results and analysis, a novel hysteresis modeling method is proposed and tested. Secondly, an AUKF method is proposed to optimize the results for a better SOC estimation for FLA batteries. Experimental results are presented to validate the proposed hysteresis model. Similarly, a comparison of simple (traditional) UKF and adaptive UKF is presented, which proves the superiority of adaptive UKF.
3. Unscented Kalman Filter
The Unscented Kalman Filter (UKF) offers a great improvement over all other filters for addressing the issue in SOC estimation. The UKF takes the approach of a minimal set of carefully selected sample points instead of using linearization. When these points are propagated through the actual nonlinear system, they capture the posterior mean and covariance with a high accuracy and use unscented transformation instead of solving a Jacobian matrix for non-linear systems, which makes them stabler and better in terms of performance. The more nonlinear a system is, the better the performance of the UKF. The UKF also consists of prediction, observation, and measurement stages. The discrete-time state equations for a non-linear system can be given as:
where
xk represents the state parameters,
f and
h are non-linear system functions,
uk is the input,
wk is the process noise vector,
vk is the measurement noise vector, and
yk is the observed measurement signal.
The Adaptive UKF algorithm can be summarized as:
where
and
are the initial state and its covariance matrix.
λ is known as a scaling factor. α determines the spread of sigma points. Its value can change between 10−4 and 1. κ is usually set to ‘0’. ‘L’ represents the number of state vectors.
- 3.
Prediction Transformation:
Propagate the sigma points through the input function and calculate the predicted state.
Calculate the propagated covariance as:
where
Qw is the covariance value of the process noise and:
β is usually equal to ‘2’ for Gaussian distributions. and are the mean and covariance weights, respectively.
- 4.
Observation Transformation:
Propagate the sigma points though the output function and calculate the predicted measurement.
Calculate the estimated covariance as:
where
Rv is the covariance value of the measurement noise.
Cross covariance is calculated as:
Calculate the Kalman gain as:
Update the state and covariance values as:
- 6.
Adjustment of Qw and Rv in AUKF:
In AUKF, the process noise covariance
Qw and measurement noise covariance
Rv do not have a constant value like in traditional UKF; instead, they are adaptively calculated in each iteration based on the covariance matching technique [
24]. The basic idea of covariance matching is to make the residuals consistent with their theoretical covariances.
Fk represents the covariance approximation of the residual voltage, which can be calculated as:
where
e is the residual voltage of the battery at time step
j and
Lq represents the window size for the covariance matching. The small window size of the covariance matching could lead the filter to diverge, whereas an adaptive filter can lose the ability of adaptation by using a large window size. Therefore, an optimal window size should be used. Based on the covariance matching technique,
Qw and
Rv can be updated as:
So, the adaptive UKF changes the values of Qw and Rv depending upon the residual voltage error in such a way to obtain covariance matching, which can further lead to minimizing the estimation error.
The state equations for a battery model can be shown as:
Cn represents the nominal capacity of the battery.
Ik and
Vk are the current and voltage of the battery.
OCVf represents the final Open-Circuit Voltage after considering the diffusion and hysteresis effect. The charge current is negative in the above equations.
wk and
vk are the process and measurement noise vectors, respectively. The overall working procedure of AUKF is shown in
Figure 8.
4. Experimental Verification
For experimentation, a 90 Ah, 12 V SFLA battery from the Sebang Company, Seoul, Republic of Korea is tested for its SOC–OCV results, where it is placed in a temperature-controlled chamber at 25 degrees Celsius [
25,
26,
27]. The OCV values were recorded during both charging and discharging phases after periods of 3 min and 3 h of rest, as illustrated in
Figure 2.
4.1. Validation of Hysteresis Modeling
To verify the hysteresis modeling at different SOC regions, a special profile was created in which SOC went to a certain value, then charged up to a certain value, and then went back again. The coulomb counting method is used as a reference for the comparison of SOC estimation.
Figure 9 shows the SOC estimation results of the proposed hysteresis modeling and fixed slopes method. The simple UKF is used for estimating the SOCs of both hysteresis models, and the values of the different parameters required in UKF are given in
Table 1.
In the proposed method, the slopes always change depending on the SOC value when the battery changes its state from discharging to charging. Normally, in previous hysteresis modeling, only one hysteresis model is used for all the regions of SOC. Therefore, a comparison of the fixed slope method and the proposed method is given.
The SOC estimation results clearly show that the proposed method of hysteresis modeling is very accurate, as the error is just about 3%. When fixed slopes are used, as the hysteresis voltage increases, the SOC estimation error also starts increasing, which clearly shows that the fixed hysteresis model is not suitable for an SFLA battery. The SOC error with the fixed slope method is about 6%.
4.2. Validation of Adaptive Unscented Kalman Filter
Generally, ampere-hour counting is a simple and intuitive method that integrates the current over time, but it suffers from cumulative errors due to current sensor drift and initial SOC uncertainty. Traditional UKF is a more advanced method that uses a non-linear battery model and exploits the statistical relationships between SOC, voltage, and current. However, the UKF assumes constant values for the process and measurement noises, which may not be valid in dynamic conditions. AUKF is an improved version of UKF that adapts the noise covariances dynamically based on the voltage error in each iteration, thus enhancing the accuracy and robustness of SOC estimation.
To validate the estimation accuracy of AUKF, a random charge/discharge profile is applied on the battery with a 100% nominal capacity, as shown in
Figure 10.
The applied profile is designed in such a way that the SOC of the battery moves in between 100% and 50% and has both charging and discharging currents, which is more practical for finding out the estimation accuracy of the proposed algorithm.
A variable resting time is also introduced in the profile to make it more dynamic and practical. The measured terminal voltage and SOC calculated by the ampere-counting method are also shown in
Figure 10.
The basic differences in both simple UKF and adaptive UKF are the values of the measurement and process noises. In simple UKF, their values are always taken as constant, while in adaptive UKF, those values are estimated based on the voltage error in each iteration, which makes AUKF more dynamic and practical in nature.
The estimation window
Lq in AUKF was taken as 20 time intervals.
Figure 11 shows a comparison of the simple and adaptive UKF algorithms. For simple UKF, the same parameters shown in
Table 1 are used.
As can be observed in
Figure 11 and
Figure 12, adaptive UKF is fast in terms of its convergence time. It converges in about 0.07 h to the right value, while the simple UKF takes about 0.2 h.
In traditional UKF, the Kalman gain is, in fact, limited because of the constant values of Rv and Qw, while in AUKF, Kalman gain can freely change depending upon the residual voltage error in the battery system, which makes AUKF fast in terms of convergence time.
Similarly, after converging to the right value SOC error, AUKF does not fluctuate very fast compared to simple UKF, whose SOC error has a high fluctuation. The AUKF takes an Lq number of previous samples of residual voltage error into account for updating Rv and Qw, which makes AUKF immune to abrupt current changes happening in a battery, which can cause voltage errors for a while as well.
Traditional UKF only takes current voltage error into account and it changes the state depending upon that error, which makes the state estimated by traditional UKF more subject to fluctuation. The overall error by simple UKF is about 5%, while by using adaptive UKF, the error is reduced to about 2%. In
Figure 13, the variations in the
Qw and
Rv are shown. Both values of these noises are updated adaptively in each iteration based on the residual voltage error.