Abstract
Owing to the degradation of an echelon-use lithium-ion battery (EULIB), the Ohmic internal resistance (OIR) and actual capacity (AE) have both changed greatly, and the state of energy (SOE) can more accurately represent the state of a EULIB than the state of charge (SOC) because of the working voltage. To improve the accuracy and adaptability of SOE estimation, in the paper, we study the energy state estimation of a EULIB. First, the four-order resistor–capacitance equivalent model of a EULIB is established, and an unscented transformation is introduced to further improve the estimation accuracy of the SOE. Second, a EULIB’s SOE is estimated based on adaptive unscented Kalman filter (AUKF), and the OIR and AE of a EULIB are estimated based on the AUKF. Third, a Takagi–Sugeno fuzzy model is introduced to optimize the OIR and AE of the EULIB, and the SOE estimation method is established based on an adaptive dual unscented Kalman filter (ADUKF). Through simulation experiments, verification, and comparison, energy decayed to 80%, 60%, and 40% of the rated energy, respectively, even with a large initial error; with the initial value of the SOE starting at 100%, 60%, or 20%, the estimated SOE can track the actual value. It can be seen that the method has a strong adaptive ability, and the estimation accuracy error is less than 1.0%, indicating that the algorithm has high accuracy. The method presented in this paper provides a new perspective for SOE estimation of EULIBs.
1 Introduction
An echelon-use lithium-ion battery (EULIB) refers to a power lithium battery with less than 80% capacity, which can be used as a backup power supply and on other occasions. Owing to the attenuation of EULIB performance, the working voltage consistency of a lithium battery is poor, so the use of the state of charge (SOC) cannot accurately represent the state of a lithium battery. Therefore, considering the influence of real-time working voltage, the state of energy (SOE) is used to represent the state of the lithium battery.
To ensure the safe, reliable, and efficient operation of a lithium-ion battery (LIB), the battery management system (BMS) plays an important role in monitoring the running process of a LIB and providing system status (; ). At the same time, the BMS accurately estimated key parameters of LIBs, such as estimating the SOC (; ), the state of health (SOH) (), SOE (), and the remaining useful life (RUL) () of LIBs. The SOE is an important parameter of a BMS, which is the ratio of the remaining available energy to the maximum available energy (; ; ), and is one of the most critical parameters in a BMS. To improve the performance of electric vehicle BMSs, a high-precision SOE estimation algorithm is needed (). SOE estimation methods for LIBs are generally divided into four categories (as shown in Figure 1): power integral, data-driven, model-based, and hybrid methods ().
FIGURE 1
A general SOE estimation method is the power integral method (detailed in ), in which a new indicator of the energetic reserve, SOE, is proposed to deal with modern BMSs. An improved scheme to estimate the SOE by mapping the relationship between discharge power, remaining energy, and SOE was proposed in . Although this method is advanced to some extent, it needs calibration and is time consuming and expensive.
With the advancements in machine learning methods, the data-driven method has been gaining popularity (). Generally, machine learning methods, for example, the genetic algorithm (), neural network (NN) (), deep neural networks (DNNs) (), and support vector regression (SVR) (), are widely used for SOE estimation. In , a method based on a feedforward NN is proposed for the joint estimation of the SOC, SOE, and power loss. The method has been verified by testing, and its prediction accuracy is high. A data-driven method is based on a long short-term memory (LSTM) to jointly estimate the SOC and SOE in . The accuracy and robustness of the method have been verified by the dynamic cycling condition. In , two different data-driven SOE estimation methods, that is, using the DNN and SVR, are compared. The SOE estimation results demonstrate the high accuracy of the DNN over SVR under the same dynamic operating conditions. proposed a method based on a radial basis function (RBF) neural network to estimate the energy recovered and the actual energy released when the battery current direction changes. The results show that this method can effectively improve SOE estimation accuracy.
The model-based method uses a resistor–capacitance equivalent circuit model (; ; ; ; ), which mainly includes particle filter and Kalman filter families. An unscented particle filter (PF) is used in to solve the non-linear and noise problems of the system. The experimental results show that the error of the SOE is less than 1.8%. In , to ensure a precise estimate of the SOE, dual adaptive particle filters are proposed based on the third-order resistor–capacitance equivalent model (TRCEM). Through the dynamic stress test and supplemental federal test, the feasibility and effectiveness of the method are verified. In , using the second-order resistor–capacitance equivalent model (SRCEM), the SOE is estimated based on an adaptive unscented Kalman filter (AUKF) algorithm. The simulation results indicate that the proposed method is adaptive, regardless of whether the initial SOE value is consistent with the real value, and shows high accuracy. A joint estimation method of the maximum available energy and SOE that adopted the Kalman filter algorithm is proposed in and Zhang, S.; Zhang, X (2021). The experimental results show that the method has good accuracy and robustness. In , a method based on a fifth-order simplex square radius cubature Kalman filter is developed to achieve SOE estimation accuracy and robustness. The error of SOE estimation is less than 3%. In , a multiple timescale dual-extended Kalman filter is used to simultaneously estimate the SOC, SOE, state of power (SOP), and SOH, and the experimental results show that the errors of the estimated SOC and SOE are all less than 1%.
It can be seen from the aforementioned literature that a single method has a limited effect on improving the accuracy of SOE estimation, and increasingly more scholars are trying hybrid methods. In , the LSTM combined with an AUKF method is proposed to estimate the SOC and SOE simultaneously. The test results show that the proposed method can effectively estimate the SOC and SOE of LIBs with high accuracy and low complexity. In , the SOE method based on a particle filter and an extended Kalman filter is presented. The experimental results show that the maximum error of SOE estimation is less than 3%. In , an approach for battery SOE prediction is proposed based on an adaptive square root unscented Kalman filter. Owing to the strong non-linear characteristics of LIBs, even with a 20% initial SOE error, the predicted SOE can converge to the actual value. The estimation error of LIB’s SOE is less than 2%. A method based on the extended Kalman filter algorithm and Markov chain model was advanced in to predict the future output power, as well as SOE. The results show that this method has high precision and good robustness.
Researchers have not only worked to improve the SOE estimation accuracy but have also studied the adaptability of SOE estimation. In , an online adaptive estimation method is proposed that can achieve a high-precision estimation of the SOC, SOH, and SOP. The effectiveness and accuracy of this method have been verified in a BMS. In , an adaptive SOE estimation method based on a cubature Kalman filter applied to series-connected LIB packs is proposed. Even if the initial error is large, in this method, it can be adjusted quickly and has high accuracy; in other words, the root mean square error is less than 2.2%. Both the single SOE (; Zhang, S.; Zhang, X.; 2021; ; ) and hybrid SOE estimation methods () enhance the adaptability and accuracy of the algorithm.
Owing to the performance attenuation of a EULIB, both the Ohmic internal resistance (OIR) and actual energy (AE) have different variations. Accurate estimation is needed to ensure accurate SOE estimation. At the same time, when an error occurs between the actual and initial energy, the estimation method should be able to quickly adjust, track the actual SOE, and have self-adaptability.
To improve the accuracy and adaptability of SOE estimation, in the paper, we study the energy state estimation of a EULIB. First, the four-order resistor–capacitance equivalent model of a EULIB is established, and an unscented transformation is introduced to further improve the estimation accuracy of SOE. Second, a EULIB’s SOE is estimated based on an AUKF, and the OIR and AE of a EULIB are estimated based on the AUKF. Third, a Takagi–Sugeno (TS) fuzzy model is introduced to optimize the OIR and AE of the EULIB, and the SOE estimation method is established based on an adaptive dual unscented Kalman filter (ADUKF). Finally, the accuracy and adaptability of the algorithm are verified and compared in simulation experiments.
The objective of the present study is to propose an SOE estimation method for EULIBs and verify the superiority of the ADUKF. The original contributions of this paper are as follows.
1) The accuracy of SOE estimation is improved, and a four-order resistor–capacitance equivalent model (FRCEM) of a EULIB is established.
2) The TS fuzzy model is introduced to optimize the OIR and AE of the EULIB and to improve the accuracy of SOE estimation.
3) The unscented transformation (UT) is introduced to establish the SOE estimation method based on the ADUKF, and the EULIB’s SOE is estimated based on the AUKF, as is its OIR and AE.
2 Estimation model of the state of energy
2.1 State of energy
Compared with the SOC, the estimation of the SOE increases the influence of the working voltage, which is of great significance in the safety evaluation of EULIBs. The SOE of a EULIB is defined as the ratio of the remaining energy to the maximum available energy and is a direct expression of the remaining available mileage of a EULIB, which can be derived as follows:where and are the SOE of EULIB at k and k + 1 time in a discrete state, respectively; is the energy of a EULIB; and are the current and working voltage of a EULIB at k time in a discrete state; and is the sampling period.
2.2 State of energy estimation model based on the FRCEM
Owing to the complex non-linear system of a EULIB, in order to simulate the characteristics of a EULIB more accurately, a higher-order battery equivalent model is required.
The FRCEM not only fully considers the influence of the polarization impedance of the EULIB on the performance, but also can reflect the state information of the EULIB in timely and accurate manner. The model has high accuracy and is easy to implement. However, the higher-order resistor–capacitance (RC) model is highly complex and has high requirements on the system, so it is difficult to apply. On the premise of considering accuracy, complexity, and practical value, the FRCEM is adopted in this article.
In Figure 2, is the open-circuit voltage (OCV) of the battery; is the working voltage of the EULIB; is the ohmic internal resistance of the EULIB; and are the polarization resistance; and are the polarization capacitance; , and are the voltages at both ends of capacitors , , and , respectively; is the charge/discharge current; and , , , and are time constants.
FIGURE 2
According to Figure 2, the discrete state equation of EULIB’s FRCEM is as follows:
According to Figure 2, the discrete observation equation of EULIB’s FRCEM is as follows:
As
So,where , , and are the estimated voltage values of , and at k time in a discrete state, respectively; , , and are the estimated voltage values of , and at k+1 time in a discrete state, respectively; are independent system noises; is the OCV of a EULIB corresponding to the SOE value of a EULIB at k time in a discrete state.
2.3 Model parameter identification
The working voltage , charge/discharge current , and the OCV of the EULIB are collected through the charge/discharge test. Model parameter identification based on the least square method is not repeated in this article because the method is described in detail in and . Both and adopt the off-line identification method (; ).
3 SOE estimation of the EULIB
In this paper, UT is introduced and the ADUKF algorithm is used to estimate the SOE and OIR and AE, respectively, which are optimized by the TS fuzzy model. The estimation accuracy and adaptability characteristics are compared and analyzed (; ).
3.1 UT
UT is the core technology of the unscented Kalman filter (UKF). The UT’s idea is to approximate the distribution of a probability density function through a set of carefully selected sample points and the corresponding weights of the sample points. A certain number of sigma points are selected from the prior distribution, according to a certain strategy, and non-linear transformation is performed on each sigma point to obtain the corresponding transformed sampling points. In addition, the posterior mean and variance are calculated by weighting these sample points. The sample points obtained by UT are neither linearized nor have lost their higher-order terms, so the accuracy of UT is higher than that of the extended Kalman filter algorithm ().
For any non-linear function , around the n-dimensional state variable with the mean value of and variance of , 2n +1 sigma points and corresponding weights can be selected by the following UT process:
In the aforementioned formula, is the column or row of the square root of . is the scale factor used to represent the range of sampling points around equilibrium value points, and is the weight of the corresponding sample points. In this paper, .
From , the sample point with non-linear transformation can be obtained:
Then, statistical characteristics of the output sampling points can be obtained by weighted approximation:
3.2 State of energy estimation based on the AUKF
From Eqs 04 and 05, the variable of a EULIB system is the SOE. The state and observation are as follows:where is the state variable OIR and AE, , is the system state variable of a EULIB; and are the input and observation variables of the system, which are the current and the working voltages of a EULIB, respectively; are the zero-mean Gaussian white noise, the error covariance matrices of are .
The AUKF algorithm is as follows:
Step 1:Initialize :
Step 2:Generating sigma points:
Step 3:Time update of :
Step 4:Status update of :The Kalman gain can be expressed as follows:The optimal estimation of the state variable can be expressed as follows:The optimal estimate of the covariance is as follows:
Step 5:Process noise covariance equation:
Step 6:Observation noise covariance equation:where , is the forgetting factor, 0< <1; and are the estimation and optimal estimation of the state variable at k time, respectively; and are the estimated value and actual observation value at k time, respectively; and are, respectively, the estimation and optimal estimation of the error covariance at k time. The subscript U in the formula represents the algorithm formula of the AUKF.According to equations 2, 3, and 23,where is the optimal value of the SOE based on the AUKF.
3.3 OIR and AE estimation based on the AUKF
The state and observation formulas of the system with the newly added state parameters are as follows:where is the noise on the input variable, and it is the zero-mean Gaussian white noise; is the noise on the output variable, and it is the zero-mean Gaussian white noise; the error covariance matrices of are , respectively; the state variable is estimated based on the AUKF algorithm, and the estimated values of the EULIB’s OIR and AE are calculated. In order to improve the accuracy, the error between the actual value and the estimated value of the working voltage is optimized.
The AUKF algorithm flow is as follows:
Step 1:Initialize :
Step 2:Generating sigma points:
Step 3:Time update of :
Step 4:Status update of :The Kalman gain is as follows:The optimal estimation of the state variable is as follows:The optimal estimate of the covariance is as follows:
Step 5:Process noise covariance equation:
Step 6:Observation of the noise covariance equation:where , and is the forgetting factor of , 0 < < 1 , .
3.4 OIR based on TS fuzzy optimization
In order to further improve the estimation accuracy of the SOE, this paper adopts the TS fuzzy model to the OIR and AE in the Kalman filter process ().
The OIR theory domain: [1 3].
The OIR fuzzy language variables: , , and .
The OIR membership function expression:
Definition:
The change rate of the OIR theory domain: [−0.3 0.6].
The change rate of OIR fuzzy language variables: , , and .
The change rate of OIR membership function expression:
According to Table 1, the optimal updated value of OIR is as follows:where is the number of fuzzy rules, is the conclusion equation from rule , and represents the membership degree of rule corresponding to this generalized input vector.
TABLE 1
Fuzzy optimization decision table for the OIR.
3.5 AE based on TS fuzzy optimization
The AE theory domain: [44 66].
The AE fuzzy language variables: , , and .
The AE membership function expression:
The change rate of the AE theory domain: [−30 50].
The change rate of AE fuzzy language variables: , , and .
The change rate of AE membership function expression:
According to Table 2, the optimal updated value of AE is as follows:where is the number of fuzzy rules, is the conclusion equation from rule , and represents the membership degree of rule corresponding to this generalized input vector.
TABLE 2
Fuzzy optimization decision table for AE.
3.6 SOE estimation based on the ADUKF
The SOE estimation flow chart of the EULIB based on the ADUKF is shown in Figure 3.
FIGURE 3
4 Simulation and discussion
4.1 Experiment
According to Figure 4 and Table 3, the EULIB is selected to conduct the charge–discharge experiment using test equipment (BTS20) at room temperature. In order to simulate the working conditions, the fully charged EULIB was discharged several times. Different discharge currents were used each time for simulation, verification, and analysis based on MATLAB R2022a.
FIGURE 4
TABLE 3
| Item | Parameter | Remark (A) |
|---|---|---|
| Nominal discharge capacity | 60 A h | 60 |
| Nominal voltage | 3.2 V | |
| Working voltage | 2.5 V–3.65 V | |
| Standard charging time | 3 h | 20 |
| Charging temperature | 25°C | |
| Discharging temperature | 25°C | |
| Rated power | 192 W |
Parameters of the battery.
In order to verify the adaptive characteristics of the ADUKF algorithm, a test experiment was carried out on a fully charged EULIB, starting from an SOE of 100% to ending at an SOE of 25%. In the paper, selected battery energy decays to 80%, 60%, and 40% separately of rated energy, and the initial SOE values were changed to 100%, 60%, and 20% separately, and the adaptive and error curves were observed and analyzed.
In the process of simulation and verification, the estimated value of the SOE was calculated based on the ADUKF algorithm, and the actual value was acquired by BTS20.
According to Eq. 27, the formula for error is as follows:
The SOE error of the ADUKF formula is as follows:where is the value acquired by a piece of test equipment.
4.2 Energy decays to 80% and the SOE starts at 100%
The simulation comparative verification curve when energy decays to 80% and the SOE starts at 100% is shown in Figure 5.
FIGURE 5
The SOE estimation curves of the EULIB are shown in Figure 5A. The bottom graph is the SOE error curve, and the top graph is the SOE adaptive estimation curve. As shown in Table 4 and Figure 5, the SOE error is 0%–0.84% based on the ADUKF with TS fuzzy optimization.
TABLE 4
| Estimation error | SOE error of ADUKF | |
|---|---|---|
| Decay to 80% | SOE = 100% | −0%–0.84% |
| SOE = 60% | −0.90%–0.96% | |
| SOE = 20% | −0.92%–0.97% | |
| Decay to 60% | SOE = 100% | −0%–0.88% |
| SOE = 60% | −0.91%–0.97% | |
| SOE = 20% | −0.94%–0.99% | |
| Decay to 40% | SOE = 100% | −0%–0.90% |
| SOE = 60% | −0.93%–0.99% | |
| SOE = 20% | −0.95%–1.00% | |
Estimation error of the EULIB’s SOE.
The observation noise covariance and process noise covariance curves of a EULIB are shown in Figure 5B. The bottom graph is the process noise covariance curve, and the top graph is the observation noise covariance curve.
According to the variation trend of observation and process noise covariance curves, the method is convergent.
4.3 Energy decays to 80% and SOE starts at 60%
The simulation comparative verification curve when energy decays to 80% and the SOE starts at 60% is shown in Figure 6.
FIGURE 6
The SOE estimation curves of EULIB are shown in Figure 6A. The bottom graph is SOE error curve, and the top graph is the SOE adaptive estimation curve. As shown in Table 4 and Figure 6, the SOE error is −0.90%–0.96% (maximum and minimum errors in the SOE error curve) based on the ADUKF with TS fuzzy optimization.
The observation noise covariance and process noise covariance curves of a EULIB are shown in Figure 6B. The bottom graph is the process noise covariance curve, and the top graph is observation noise covariance curve.
According to the variation trend of observation and process noise covariance curves, the method is convergent. The change in observation noise covariance is small, but the change in process noise covariance is large, compared with that when SOE is 100%. Since the initial value is different, the algorithm is adjusted and the actual value is tracked very high in the end.
4.4 Energy decays to 80% and the SOE starts at 20%
The simulation comparative verification curve when energy decays to 80% and the SOE starts at 20% is shown in Figure 7.
FIGURE 7
The SOE estimation curves of a EULIB are shown in Figure 7A. The bottom graph is the SOE error curve, and the top graph is the SOE adaptive estimation curve. As shown in Table 4 and Figure 7, the SOE error is −0.92%–0.97% (maximum and minimum errors in the SOE error curve) based on the ADUKF with TS fuzzy optimization.
The observation noise covariance and process noise covariance curves of a EULIB are shown in Figure 7B The bottom graph is process noise covariance curve, and the top graph is the observation noise covariance curve.
According to the variation trend of observation and process noise covariance curves, the method is convergent. The changes in observation and process noise covariance are both increased, compared with that when the SOE is 60%. Since the initial value is different, the algorithm is adjusted and the actual value is tracked very high in the end.
4.5 Energy decays to 60% and SOE starts at 100%
The simulation comparative verification curve when energy decays to 60% and the SOE starts at 100% is shown in Figure 8.
FIGURE 8
The SOE estimation curves of a EULIB are shown in Figure 8A. The bottom graph is the SOE error curve, and the top graph is the SOE adaptive estimation curve. As shown in Table 4 and Figure 8, the SOE error is 0%–0.88% (maximum and minimum errors in the SOE error curve) based on the ADUKF with TS fuzzy optimization.
The observation noise covariance and process noise covariance curves of a EULIB are shown in Figure 8B. The bottom graph is process noise covariance curve, and the top graph is the observation noise covariance curve.
According to the variation trend of observation and process noise covariance curves, the method is convergent. As the energy decays more, the algorithm needs stronger adjustment ability and the error is larger than that before the decays.
4.6 Energy decays to 60% and the SOE starts at 60%
The simulation comparative verification curve when energy decays to 60% and the SOE starts at 60% is shown in Figure 9.
FIGURE 9
The SOE estimation curves of a EULIB are shown in Figure 9A The bottom graph is the SOE error curve, and the top graph is the SOE adaptive estimation curve. As shown in Table 4 and Figure 9, the SOE error is −0.91%–0.97% (maximum and minimum errors in the SOE error curve) based on the ADUKF with TS fuzzy optimization.
The observation noise covariance and process noise covariance curves of a EULIB are shown in Figure 9B. The bottom graph is the process noise covariance curve, and the top graph is the observation noise covariance curve.
According to the variation trend of observation and process noise covariance curves, the method is convergent. With the intensification of energy attenuation, the algorithm needs stronger adjustment ability and the error is larger than that of energy attenuation to 80%. However, tracking the actual value is good.
4.7 Energy decays to 60% and the SOE starts at 20%
The simulation comparative verification curve when energy decays to 60% and the SOE starts at 20% is shown in Figure 10.
FIGURE 10
The SOE estimation curves of a EULIB are shown in Figure 10A The bottom graph is the SOE error curve, and the top graph is the SOE adaptive estimation curve. As shown in Table 4 and Figure 10, the SOE error is −0.94%–0.99% (maximum and minimum errors in the SOE error curve) based on the ADUKF with TS fuzzy optimization.
The observation noise covariance and process noise covariance curves of a EULIB are shown in Figure 10B The bottom graph is the process noise covariance curve, and the top graph is the observation noise covariance curve.
According to the variation trend of observation and process noise covariance curves, the method is convergent. With the intensification of energy attenuation, the algorithm needs stronger adjustment ability and the error is larger than that of energy attenuation to 80%. However, tracking the actual value is good.
4.8 Energy decays to 40% and the SOE starts at 100%
The simulation comparative verification curve when energy decays to 40% and the SOE starts at 100% is shown in Figure 11.
FIGURE 11
The SOE estimation curves of a EULIB are shown in Figure 11A. The bottom graph is the SOE error curve, and the top graph is the SOE adaptive estimation curve. As shown in Table 4 and Figure 11, the SOE error is 0%–0.90% (maximum and minimum errors in the SOE error curve) based on ADUKF with TS fuzzy optimization.
The observation noise covariance and process noise covariance curves of a EULIB are shown in Figure 11B The bottom graph is the process noise covariance curve, and the top graph is the observation noise covariance curve.
According to the variation trend of observation and process noise covariance curves, the method is convergent. As the energy decays more, the algorithm needs stronger adjustment ability and the error is larger than that before the decays.
4.9 Energy decays to 40% and the SOE starts at 60%
The simulation comparative verification curve when energy decays to 40% and the SOE starts at 60% is shown in Figure 12.
FIGURE 12
The SOE estimation curves of a EULIB are shown in Figure 12A, The bottom graph is the SOE error curve, and the top graph is the SOE adaptive estimation curve. As shown in Table 4 and Figure 12, the SOE error is −0.93%–0.99% (maximum and minimum errors in the SOE error curve) based on the ADUKF with TS fuzzy optimization.
The observation noise covariance and process noise covariance curves of a EULIB are shown in Figure 12B The bottom graph is the process noise covariance curve, and the top graph is the observation noise covariance curve.
According to the variation trend of observation and process noise covariance curves, the method is convergent. With the intensification of energy attenuation, the algorithm needs stronger adjustment ability and the error is larger than that of energy attenuation to 60%. However, tracking the actual value is good.
4.10 Energy decays to 40% and the SOE starts at 20%
The simulation comparative verification curve when energy decays to 40% and the SOE starts at 20% is shown in Figure 13.
FIGURE 13
The SOE estimation curves of a EULIB are shown in Figure 13A The bottom graph is the SOE error curve, and the top graph is the SOE adaptive estimation curve. As shown in Table 4 and Figure 13, the SOE error is −0.95%–1.00% (maximum and minimum errors in the SOE error curve) based on the ADUKF with TS fuzzy optimization.
The observation noise covariance and process noise covariance curves of a EULIB are shown in Figure 13B. The bottom graph is the process noise covariance curve, and the top graph is the observation noise covariance curve.
According to the variation trend of observation and process noise covariance curves, the method is convergent. With the intensification of energy attenuation, the algorithm needs stronger adjustment ability and the error is larger than that of energy attenuation to 60%. However, tracking the actual value is good.
4.11 Discussion
The simulation results show that the estimation accuracy error of the EULIB’s SOE is less than 1.00% based on a TS fuzzy-optimized ADUKF. With the continuous attenuation of energy, the algorithm is required to have increasingly more adaptive ability, and although the error also increases, the maximum error is less than 1%. As shown in Table 5, compared with the SOE estimation accuracy of 1.8% in , 2% in , 3% in and , 2.34% in , and 1% in , the method proposed in this paper has higher accuracy. Furthermore, the proposed method is adaptive, regardless of whether the initial SOE value is consistent with the actual value.
TABLE 5
| Reference | Accuracy of estimation (%) | Adaptability | Method |
|---|---|---|---|
| Method of this paper | 1 | Yes | TS fuzzy optimized ADUKF |
| 1.8 | No | Unscented particle filter | |
| 2 | Yes | Adaptive square root unscented Kalman filter | |
| 3 | No | Particle swarm optimization and the cubature Kalman filter | |
| 3 | Yes | Particle filter and the extended Kalman filter | |
| 2.34 | Yes | ADUKF | |
| 1 | No | Multi-time scale dual-extended Kalman filter |
Comparison of optimization methods.
5 Conclusion
In this paper, the SOE estimation method based on a TS fuzzy optimization ADUKF is established of a EULIB. To improve the accuracy of SOE estimation, the four-order resistor–capacitance equivalent model is used to estimate the SOE, OIR, and AE based on an adaptive dual unscented Kalman filter. To improve the accuracy of the estimation model, the OIR and AE have been optimized based on the TS fuzzy model. To enhance the adaptive ability, the process error and the observation error are estimated. Through simulation, energy decayed to 80%, 60%, and 40% of the rated energy; with the initial value of the SOE starting at 100%, 60%, or 20%, the estimated SOE can track the actual value. The estimation accuracy error of the EULIB’s SOE is less than 1.00% based on a TS fuzzy optimization ADUKF, and the algorithm has higher accuracy.
Also, positive results have been achieved in this paper, both in accuracy and adaptability. However, due to the high-order equivalent model, compared with the low-order model, the calculation amount is larger.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
EH: writing—original draft preparation, conceptualization, validation, and methodology. ZXW: project administration, funding acquisition, and resources. ZW: writing—review and editing and supervision. XQ: visualization, data curation, and software. GL: formal analysis and investigation. All authors have read and agreed to the published version of the manuscript.
Funding
This work was supported by the Science and Technology Major Project of Inner Mongolia Autonomous Region (2020ZD0014). This research was conducted within the University of Shandong Jiao Tong, supported by the charging discharging equipment BTS20 (5V/4*300A/WD).
Acknowledgments
The authors thank the testing team members for their support with the experiment.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AnF.JiangJ.ZhangW.ZhangC.FanX. (2022). State of energy estimation for lithium-ion battery pack via prediction in electric vehicle applications. IEEE Trans. Veh. Technol.71 (1), 184–195. 10.1109/tvt.2021.3125194
2
BaraiA.UddinK.WidanalageW. D.McGordonA.JenningsP. (2016). The effect of average cycling current on total energy of lithium-ion batteries for electric vehicles. J. Power Sources303, 81–85. 10.1016/j.jpowsour.2015.10.095
3
ChangJ.ChiM.ShenT. (2020). Model based state-of-energy estimation for LiFePO4 batteries using unscented particle filter. J. Power Electron.20 (2), 624–633. 10.1007/s43236-020-00051-5
4
ChenY.YangX.LuoD.WenR. (2021). Remaining available energy prediction for lithium-ion batteries considering electrothermal effect and energy conversion efficiency. J. Energy Storage40, 102728. 10.1016/j.est.2021.102728
5
FanT.-E.LiuS.-M.TangX.QuB. (2022). Simultaneously estimating two battery states by combining a long short-term memory network with an adaptive unscented Kalman filter. J. Energy Storage50, 104553. 10.1016/j.est.2022.104553
6
GaoT.LuW. (2021). Machine learning toward advanced energy storage devices and systems. iScience24 (1), 101936. 10.1016/j.isci.2020.101936
7
HeP.WangC.ZhaoW.WangW.WuG.ChangC. (2022). A SOE estimation method for lithium batteries considering available energy and recovered energy. Proc. Institution Mech. Eng. Part D J. Automob. Eng.237, 273–290. 10.1177/09544070211070441
8
Hossain LipuM. S.AnsariS.MiahM. S.MerajS. T.HasanK.ShihavuddinA. S. M.et al (2022). Deep learning enabled state of charge, state of health and remaining useful life estimation for smart battery management system: Methods, implementations, issues and prospects. J. Energy Storage55, 105752. 10.1016/j.est.2022.105752
9
Hossain LipuM. S.HannanM. A.HussainA.AyobA.SaadM. H. M.KarimT. F.et al (2020). Data-driven state of charge estimation of lithium-ion batteries: Algorithms, implementation factors, limitations and future trends. J. Clean. Prod.277, 124110. 10.1016/j.jclepro.2020.124110
10
HouE.QiaoX.LiuG. (2017). Estimation of power lithium-ion battery SOC based on fuzzy optimal decision. Chin. J. Power Sources41, 920–922.
11
HouE.QiaoX.LiuG. (2014). Modeling and simulation of power lithium-ion battery SOC. Chin. Comput. Simul.31, 193–196.
12
HouE.WangZ.QiaoX.LiuG. (2022). Remaining useful cycle life prediction of lithium-ion battery based on TS fuzzy model. Front. Energy Res.10. 10.3389/fenrg.2022.973487
13
HouE.XuY.QiaoX.LiuG.WangZ. (2022). Research on state of power estimation of echelon-use battery based on adaptive unscented kalman filter. Symmetry14 (5), 919. 10.3390/sym14050919
14
HouE.XuY.QiaoX.LiuG.WangZ. (2021). State of power estimation of echelon-use battery based on adaptive dual extended kalman filter. Energies14 (17), 5579. 10.3390/en14175579
15
HuJ. N.HuJ. J.LinH. B.LiX. P.JiangC. L.QiuX. H.et al (2014). State-of-charge estimation for battery management system using optimized support vector machine for regression. J. Power Sources269, 682–693. 10.1016/j.jpowsour.2014.07.016
16
HuX.LiS. E.YangY. (2016). Advanced machine learning approach for lithium-ion battery state estimation in electric vehicles. IEEE Trans. Transp. Electrification2 (2), 140–149. 10.1109/tte.2015.2512237
17
HuX.YuanH.ZouC.LiZ.ZhangL. (2018). Co-estimation of state of charge and state of health for lithium-ion batteries based on fractional-order calculus. IEEE Trans. Veh. Technol.67 (11), 10319–10329. 10.1109/tvt.2018.2865664
18
LaiX.HuangY.HanX.GuH.ZhengY. (2021). A novel method for state of energy estimation of lithium-ion batteries using particle filter and extended Kalman filter. J. Energy Storage43, 103269. 10.1016/j.est.2021.103269
19
LiK.WeiF.TsengK. J.SoongB.-H. (2018). A practical lithium-ion battery model for state of energy and voltage responses prediction incorporating temperature and ageing effects. IEEE Trans. Industrial Electron.65 (8), 6696–6708. 10.1109/tie.2017.2779411
20
LiX.XuJ.HongJ.TianJ.TianY. (2021). State of energy estimation for a series-connected lithium-ion battery pack based on an adaptive weighted strategy. Energy214, 118858. 10.1016/j.energy.2020.118858
21
MaD.GaoK.MuY.WeiZ.DuR. (2022). An adaptive tracking-extended kalman filter for SOC estimation of batteries with model uncertainty and sensor error. Energies15 (10), 3499. 10.3390/en15103499
22
MaL.HuC.ChengF. (2021). State of charge and state of energy estimation for lithium-ion batteries based on a long short-term memory neural network. J. Energy Storage37, 102440. 10.1016/j.est.2021.102440
23
MamadouK.LemaireE.DelailleA.RiuD.HingS. E.BultelY. (2012). Definition of a state-of-energy indicator (SoE) for electrochemical storage devices: Application for energetic availability forecasting. J. Electrochem. Soc.159 (8), A1298–A1307. 10.1149/2.075208jes
24
PurohitK.SrivastavaS.NookalaV.JoshiV.ShahP.SekharR.et al (2021). Soft sensors for state of charge, state of energy, and power loss in formula student electric vehicle. Appl. Syst. Innov.4 (4), 78. 10.3390/asi4040078
25
QiaoX.WangZ.HouE.LiuG.CaiY. (2022). Online estimation of open circuit voltage based on extended kalman filter with self-evaluation criterion. Energies15 (12), 4373. 10.3390/en15124373
26
QinD.LiJ.WangT.ZhangD. (2019). Modeling and simulating a battery for an electric vehicle based on modelica. Automot. Innov.2 (3), 169–177. 10.1007/s42154-019-000660
27
RahimifardS.AhmedR.HabibiS. (2021). Interacting multiple model strategy for electric vehicle batteries state of charge/health/power estimation. IEEE Access9, 109875–109888. 10.1109/access.2021.3102607
28
ShrivastavaP.SoonT. K.IdrisM. Y. I. B.MekhilefS.AdnanS. B. R. S. (2022). Model‐based state of X estimation of lithium‐ion battery for electric vehicle applications. Int. J. Energy Res.46, 10704–10723. 10.1002/er.7874
29
ShrivastavaP.Tey KokS.Yamani Bin IdrisM.MekhilefS.Bahari Ramadzan Syed AdnanS. (2021). Lithium-ion battery state of energy estimation using deep neural network and support vector regression. 10.1109/ECCE-Asia49820.2021.9479413
30
TranM. K.PanchalS.KhangT. D.PanchalK.FraserR.FowlerM. (2022). Concept review of a cloud-based smart battery management system for lithium-ion batteries: Feasibility, logistics, and functionality. Batter. (Basel)8 (2), 19. 10.3390/batteries8020019
31
WangQ.WangZ.ZhangL.LiuP.ZhouL. (2022). A battery capacity estimation framework combining hybrid deep neural network and regional capacity calculation based on real-world operating data. IEEE Trans. Industrial Electron., 1–11. 10.1109/tie.2022.3229350
32
WangS.ChenQ.WangK.WangZ.WangY. (2019). Parameter identification of lithium iron phosphate battery model for battery electric vehicle. IOP Conf. Ser. Mater. Sci. Eng.677 (3), 032048. 10.1088/1757-899x/677/3/032048
33
WangS.FanY.JinS.Takyi-AninakwaP.FernandezC. (2023). Improved anti-noise adaptive long short-term memory neural network modeling for the robust remaining useful life prediction of lithium-ion batteries. Reliab. Eng. Syst. Saf.230. 10.1016/j.ress.2022.108920
34
XuW.XuJ.LangJ.YanX. (2019). A multi-timescale estimator for lithium-ion battery state of charge and state of energy estimation using dual H infinity filter. IEEE Access7, 181229–181241. 10.1109/access.2019.2959396
35
YongF.FanB.ChenX.FenglingH. (2021). State-of-Charge and state-of-energy estimation for lithium-ion batteries using sliding-mode observers.
36
ZhaiG.LiuS.WangZ.ZhangW.MaZ. (2017). State of energy estimation of lithium titanate battery for rail transit application. Energy Procedia105, 3146–3151. 10.1016/j.egypro.2017.03.681
37
ZhangS.PengN.ZhangX. (2021). A variable multi-time-scale based dual estimation framework for state-of-energy and maximum available energy of lithium-ion battery. Int. J. Energy Res.46, 2876–2892. 10.1002/er.7350
38
ZhangS.ZhangX. (2021). Joint estimation method for maximum available energy and state-of-energy of lithium-ion battery under various temperatures. J. Power Sources506, 230132. 10.1016/j.jpowsour.2021.230132
39
ZhangZ.ZhangX.MaoS.HuaiR.YuZ. (2021). Estimation of state‐of‐energy for lithium batteries based on dual adaptive particle filters considering variable current and noise effects. Int. J. Energy Res.45 (11), 15921–15935. 10.1002/er.6823
40
ZhouX.ZhouS.LiQ.FanZ.huaY.YangS.et al (2021). Innovative method for state of energy estimation based on improved Cubature Kalman filter. J. Phys. Conf. Ser.1754 (1), 012069. 10.1088/1742-6596/1754/1/012069
Summary
Keywords
state of energy (SOE), echelon-use lithium-ion battery (EULIB), Takagi–Sugeno (TS) fuzzy optimization, four-order resistor–capacitance equivalent model (FRCEM), adaptive dual unscented Kalman filter (ADUKF)
Citation
Hou E, Wang Z, Wang Z, Qiao X and Liu G (2023) State of energy estimation of the echelon-use lithium-ion battery based on Takagi–Sugeno fuzzy optimization. Front. Energy Res. 11:1137358. doi: 10.3389/fenrg.2023.1137358
Received
04 January 2023
Accepted
06 February 2023
Published
10 March 2023
Volume
11 - 2023
Edited by
Shunli Wang, Southwest University of Science and Technology, China
Reviewed by
Sharifuddin Mondal, National Institute of Technology Patna, India
Lei Zhang, Beijing Institute of Technology, China
Updates
Copyright
© 2023 Hou, Wang, Wang, Qiao and Liu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhixue Wang, 215035@sdjtu.edu.cn
This article was submitted to Energy Storage, a section of the journal Frontiers in Energy Research
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.