1. Introduction
Because of its high energy density, extended cycle life, and environmental friendliness, lithium-ion batteries find extensive application in a variety of fields, including aerospace, special equipment, and microgrid energy storage [
1,
2]. Lithium-ion batteries will progressively transition to an unstable state as a result of an internal chemical reaction with prolonged operation. Lithium-ion batteries have the potential to explode and destroy equipment if they are used in an unstable state. Lithium-ion battery health must, therefore, be evaluated in a precise, secure, and trustworthy manner [
3].
In general, the micro-health parameters of lithium-ion batteries refer to microscopic indicators that can reflect the performance of the active substances and electrolytes inside the battery. Changes in these parameters are usually related to the chemical reaction of the battery, the aging of the material, and changes in the internal structure of the battery [
4]. By monitoring changes in micro-health parameters, the internal health state of the battery can be determined. Current studies based on the health status assessment of lithium-ion batteries are mainly divided into the following three types: data-driven model, physical model, and hybrid model. The physical method of the lithium-ion battery health assessment is to determine the health status of the battery by analyzing its electrochemical reaction, internal structural changes, and electrical performance parameters. Amuta et al. proposed a method for evaluating the health status of lithium-ion batteries based on voltage integration [
5]. The method calculates the health status of a similar battery by integrating the voltage at a given ambient temperature and constant current charge. Su et al. proposed a method for evaluating the health status of lithium-ion batteries based on the distribution of relaxation time (DRT) [
6]. Asiedu–Asante proposed a method for monitoring the health status of lithium-ion batteries via a frequency–domain reflectometry (FDR), which assesses battery aging by analyzing the correlation between the impedance measured by FDR and the equivalent series resistance (ESR) measured by EIS [
7]. Wang et al. proposed a variable separation algorithm based on first- and second-order quasi-Newton methods to estimate and predict the health state of a battery by reducing the dimension and increasing the convergence rate [
8]. However, the accuracy of physics-based techniques is limited because of intricate electrochemical reactions that occur inside lithium-ion batteries [
9].
A data-driven approach enables accurate assessment of lithium-ion battery health through historical data and machine learning. Singh et al. proposed a deep machine learning prediction technique that utilizes encoders and decoders to extract important features and achieve high-precision Li-ion battery charge state prediction through time series analysis [
10]. Teixeira et al. proposed a model for estimating the health state curve of lithium-ion batteries using gated cycle unit (GRU) neural networks, which provides a high-precision estimation of the state charge curve for smartphone battery exchange applications with low computational complexity and cost advantages [
11]. Alwabli proposed a method for battery health analysis using logistic regression and convolutional neural networks, which further improved the performance indicators of convolutional neural networks through particle swarm optimization, a process that is significantly superior to other models [
12]. Lin et al. suggested using feature optimization and convolutional neural networks to estimate the health state of lithium-ion batteries [
13]. The method extracts a variety of electrical, thermodynamic and electrochemical characteristics; combines principal component analysis and convolutional neural network; and optimizes the feature dimension to achieve good experimental results. While most data-driven models can yield relatively accurate estimates, these methods are often black-box models and lack interpretability.
The hybrid approach to lithium-ion battery health status assessment is to improve the accuracy and reliability of the assessment by combining two or more assessment techniques and feature extraction methods for modeling. Wen et al. proposed an N-CatBoost hybrid framework for accurate estimation of lithium-ion battery health status and its uncertainty, demonstrating superior accuracy and interpretability to other machine learning algorithms [
14]. Wang et al. proposed an MFE–GRU TCA hybrid model to accurately predict the health status, which showed root–mean–square errors (RMSE) of less than 0.832% and 0.614% on the NASA and CALCE datasets, respectively [
15]. Yang et al. proposed an evidential inference rule health assessment method for lithium-ion batteries based on dynamic reference values [
16]. This method improves the accuracy and robustness of the assessment through the dynamic adjustment and whale optimization algorithm (WOA), and the experimental results show that it has good generalization ability. Yuan et al. proposed a hybrid neural network based on variational mode decomposition and CNN–Transformer for efficient prediction of lithium-ion battery health, which performs well on CALCE datasets and has good generalization ability on NASA datasets [
17]. However, based on variational mode decomposition and CNN–Transformer, it is a modeling method that combines the two models. This fusion method will not only increase the complexity of the algorithm, but also has higher requirements for the fusion strategy [
18]. The DO-BRB-I model proposed in this paper can directly integrate expert knowledge and experience into the model, and has good interpretability. At the same time, because the BRB model adopts evidential reasoning (ER) for inference, its structure is simple, transparent, and easy to understand. This approach avoids the black-box characteristics of complex deep learning models and provides a more transparent decision-making basis when dealing with uncertainty and nonlinear problems of complex systems. The hybrid model is usually a balance of the physical model and data-driven model, which has better accuracy and interpretability [
19].
The belief rule base (BRB) is a hybrid model based on data and knowledge [
20]. It is a model proposed by Yang et al. in 2006 through the Dempster–Shafer (D–S) theory and IF–THEN rule [
21]. A nonlinear modeling technique called BRB can represent several types of uncertain information, including ignorance and randomness [
22]. In addition, because of its modeling method that combines expert knowledge and IF–THEN rules, BRB has strong causal reasoning ability and good interpretability [
23]. Therefore, in recent years, many researchers have used BRB to evaluate the health status of lithium-ion batteries. For example, Han et al. proposed a lithium-ion battery health assessment model based on interpretability belief rule base (BRB-I), which addresses uncertainty and interpretability deficiencies and improves optimization results through WOA [
3]. Zhao et al. proposed a method that combines an approximate belief rule base with a hidden Markov model to estimate lithium-ion battery capacity in orbit, using historical data and expert knowledge [
24]. The model is validated with satellite battery performance data and has proved effective for in-orbit capacity estimation.
However, in these studies, parameters such as belief degree were directly divided by experts. While expert knowledge can provide broad and scientific guidance in one context, there may be limitations to the applicability of such guidance to more specific data sets. In addition, the problem that the randomness of the optimization algorithm in BRB will destroy the interpretability of BRB should be further studied. Therefore, this paper proposed a lithium-ion battery health assessment method based on the double optimization belief rule base with interpretability (DO-BRB-I). In this method, the belief degree is doubly optimized, and four interpretability constraint strategies are proposed in the optimization process, which further enhances the interpretability of the model. The main contributions of this paper are as follows:
- (1)
A belief degree optimization method for Gaussian membership function with Bayesian updating (GMF-B) is proposed. This method can optimize the belief degree according to the data distribution while maintaining the original expert knowledge.
- (2)
To improve the accuracy of the model, the projection covariance matrix adaptive evolution strategy (P-CMA-ES) is used to further optimize the other parameters, such as the belief degree, to form a double optimization.
- (3)
In view of the randomness of the optimization algorithm, four interpretability constraint strategies are proposed to constrain the interpretability based on the characteristics of lithium-ion batteries.
The structure of this paper is as follows: Problem formulation and construction of DO-BRB-I model are presented in
Section 2. The optimization process of the model is presented in
Section 3.
Section 4 provides the reasoning process of the model.
Section 5 contains a case study. Finally, the conclusion is given in
Section 6.
2. Problem Formulation and DO-BRB-I Model Construction
This part first describes the problems existing in the health status assessment of lithium-ion batteries based on BRB. Then, the DO-BRB-I model is constructed.
2.1. Problem Formulation
This study aims to solve the following three problems in the health status assessment model of lithium-ion batteries based on BRB:
Problem 1: How to obtain a more reasonable belief distribution by combining the data distribution while maintaining the original expert knowledge.
Although experts can provide extensive and scientific guidance on the classification of belief degrees, this guidance is often based on the accumulation of experience and theory [
25]. However, in the application of specific data sets, this guidance may have certain limitations. In particular, expert recommendations are often directed at general trends or phenomena, while in actual data, the true characteristics of the data may not be fully reflected due to environmental factors, differences in data distribution, or other complexities. Therefore, a new GMF-B method is proposed in this paper. This method can combine the reference values of the data set and the antecedent attributes, and then obtain a more reasonable belief distribution through Gaussian membership function and Bayesian updating. The process of the method can be described as follows:
where
represents the reference value of the preceding attribute,
represents the initial reference value determined by expert knowledge, and
represents the input data.
represents a more reasonable belief degree after optimization, and
represents a function of the GMF-B method.
Problem 2: How to get a more realistic belief distribution and then further enhance the accuracy of the model.
After solving problem 1, according to the data distribution, the confidence degree divided by experts is adjusted to obtain a more reasonable belief distribution. The initial belief degree of these optimizations provides a good basis for subsequent optimizations. The belief degree is optimized using P-CMA-ES to increase the model’s accuracy even further. At the same time, in this step, rule weights and attribute weights, two important parameters that affect the model results, are optimized for the first time. The whole optimization process can be described as follows:
where
represents the set of parameters in the optimization process, which includes belief degree
, rule weight
, and attribute weight
.
represents the input data,
represents the inference function, and
represents other parameters in the optimization process.
Problem 3: How to restrain the behavior that destroys interpretability in optimization process.
The P-CMA-ES adopted in problem 2 is a global optimization algorithm with randomness. During optimization, this randomness may destroy the interpretability of the initial parameters [
20]. Therefore, it is necessary to propose a series of constraint strategies to limit the behavior that destroys interpretability. These strategies are described as follows:
where
represents the number of constraint strategies.
After adding constraints to the optimization algorithm, Equation (3) should be updated as follows:
2.2. Construction of DO-BRB-I Model
Lithium-ion battery health status assessment model based on DO-BRB-I supplemented the original BRB model with reasonable adjustment to the initial expert setting belief degree, and added interpretability constraint strategy in the optimization process. It is the
k-th rule is described as follows:
where
refers to the indicators of lithium-ion battery health status assessment,
represents the reference value set of the assessment indicators,
is the
health status assessment results of DO-BRB-I,
represents the belief degree corresponding to each result under the belief rule
,
represents the rule weight of the
i-th belief rule, and
represents the number of belief rules.
represents the attribute weight of the
i-th assessment indicator, and
represents the number of assessment indicators.
stands for
interpretable constraint strategies.
The constructed DO-BRB-I model for lithium-ion battery health status assessment is shown in
Figure 1. The whole framework can be divided into the following steps:
Step 1: Create the basic BRB model using the expert knowledge and input data.
Step 2: The GMF-B method constructed in this paper is used to optimize the belief degree constructed by the initial expert knowledge, and a more reasonable belief distribution is obtained.
Step 3: To further improve the progress of the model, the P-CMA-ES optimization algorithm is used to further optimize other parameters such as the belief degree. However, to prevent the interpretability of expert knowledge from being destroyed in the optimization process, four interpretability constraint strategies are added.
Step 4: According to the optimal parameters, the evidential reasoning (ER) algorithm is used to deduce the final lithium-ion battery health status results.
3. Optimization Process of DO-BRB-I Model
In this section, the initial optimization of the belief degree based on the data distribution is first performed in
Section 3.1. The interpretability constraint strategy is then described in
Section 3.2. The belief degree and other parameters are optimized twice in
Section 3.3.
3.1. The First Optimization Based on GMF-B
Initial belief degrees provided by experts are correct in the general direction, but may not be completely accurate for specific data sets and application scenarios, such as specific battery models. However, a good initial belief degree is an important guide for subsequent optimization. Therefore, it is very important to obtain a reasonable initial belief degree. The Gauss membership function provides a statistical and probabilistic representation of the membership degree of elements to fuzzy sets in fuzzy systems [
26]. This membership degree reflects the influence of the data distribution on belief degree to some extent, but the initial expert knowledge reflects the expert’s guidance and consideration in the overall direction. Therefore, the initial expert knowledge must be integrated when optimizing the belief degree. Bayesian updating is a method that uses new evidence to update prior probabilities, and it is very effective in many statistical and machine learning tasks. Therefore, a new GMF-B method is proposed in this section, which can be combined with the initial belief distribution and data distribution to optimize and obtain more reasonable belief degree, to further improve the accuracy and reliability of the model. The specific steps of this method are as follows:
Step 1: Preliminary preparation
Firstly, reference values of the preceding attributes and results of expert knowledge are and , respectively, which will be used to calculate the subsequent membership degree. Then, according to the initial belief degree defined by experts, a prior probability distribution matrix is generated, denoted as .
Step 2: Calculate the membership of the reference value
The membership degree for the reference value of the preceding attribute and the membership degree for the result reference value are calculated, respectively, according to the calculation formula of the Gaussian membership degree function, as shown in the following formulas:
Equation (7) is used to calculate the membership degree between the reference value of the preceding attribute and the input data. Equation (8) is used to calculate the membership degree between the reference value of the result and the input data. stands for standard deviation based on the data distribution.
Step 3: Calculate the comprehensive membership degree
The comprehensive membership corresponding to the result can be obtained by multiplying the membership of the reference value of the preceding attribute in the current rule with the membership of the result. The formula is as follows:
where
represents the comprehensive membership degree corresponding to the
i-th result in the current rule,
represent the membership degree of
antecedent attribute reference values, and
represents the membership degree of the
i-th result.
Step 4: Perform Bayesian updates
The above three steps calculate the comprehensive membership degree of the result according to the input data and reference values. This membership can be used as the data’s influence on belief degree. Therefore, the initial expert knowledge is used as the prior probability distribution, and the above membership degree is used as new evidence to update the prior probability by Bayesian method. This process can be expressed as follows:
where
is the likelihood function expressed as the sum of the comprehensive membership degree.
is the prior probability of the current combination and
is the posterior probability of the current combination.
Step 5: Normalized operation
Repeat the first four steps to calculate the posterior probability of each set of input data and each combination, and carry out normalization to get the final posterior probability distribution matrix. This posterior probability distribution matrix is the belief distribution of the expert knowledge optimized by the data distribution.
3.2. Interpretability Constraint Strategy
After completing the preliminary optimization of belief degree in
Section 3.1, a more reasonable belief distribution is obtained. To further improve the accuracy of the model, it is beneficial to further optimize the belief degree and other parameters. However, the P-CMA-ES selected in this paper is a random global optimization, and, if it is not constrained, the interpretability of the original expert knowledge will be seriously damaged. Therefore, four interpretability constraint strategies are proposed in this section to constrain the second optimization process of the model according to the actual lithium-ion battery health status background. The following is a detailed introduction to these four strategies:
Strategy 1. Ensure that activated rules participate in optimization and reasoning
The interpretive BRB model incorporates every potential state combination for evaluating the health status of lithium-ion batteries. However, due to the possibility of inadequate observational data, not all rules may be activated by the input data [
25]. Experts create a comprehensive rule base from an international standpoint, and data gathered from a particular lithium battery type may not activate all the rules, leaving some of them dormant and excluded from the BRB reasoning process. To keep valuable expert original information, certain non-activation rules should not be included in training. However, earlier research treated each parameter as an individual in a population and utilized a global optimization approach for model optimization. In global optimization, these individuals go through an evolutionary process, but for non-activation rules in BRB, this approach is not reasonable. Therefore, it is necessary to consider a reasonable way to identify non-active rules and retain their original relevant parameters as shown by the following:
Activation weights can be used to identify rules that are not activated. The activation of the
k-th rule can be expressed as follows:
where
indicates the number of input data, and
represents the active weight of the current rule for all input data. If all activation weights are 0, the current rule is never activated. Therefore, for this rule, the parameters in the
should be preserved.
Strategy 2. Set the optimization range of the parameter
Expert knowledge provides important guidance in the assessment of the health status of lithium-ion batteries [
20]. However, when it comes to a specific battery type, it needs to be optimized according to the actual situation to ensure the accuracy of the model. To maintain the integrity of expert knowledge in the optimization process, it is necessary to set a reasonable optimization range for each parameter. This approach not only makes use of expert experience, but can also be effectively adjusted in the actual application of a specific battery model, thus achieving a balanced optimization result. This process can be described as follows:
where
and
represent the upper and lower bounds of the optimization range, respectively.
Strategy 3. Ensure reasonable belief distribution
Due to the randomness of the optimization algorithm, some rules that do not conform to the health status of lithium ions may be generated to blindly obtain higher accuracy [
25].
In
Figure 2, several belief distributions that may be presented during the optimization process are shown, and, H1 to H4, respectively, represent several health states of lithium-ion batteries. The several scenarios shown in
Figure 2b are unreasonable because the health of a lithium-ion battery cannot be both good and bad at the same time [
3]. Therefore, a reasonable belief distribution should be either monotonic or convex. Therefore, the distribution of belief should be strictly constrained in the optimization process to make it conform to a reasonable distribution. The constraints on belief can be expressed as follows:
Strategy 4. Punish behavior that goes beyond limits.
Although the boundary constraint is added in the optimization process, the step length setting of the optimization algorithm may not be suitable for all parameters and iteration processes. In the process of optimization algorithm exploration, there may still be a phenomenon that exceeds the boundary. Although this parameter may obtain higher precision, its transgression behavior has destroyed the original interpretability, so it needs to be punished for this behavior. The principle of punishment is to add an extra value to this parameter so that it will be eliminated in the subsequent screening process. This value is calculated as follows:
where
represents the initial population generated by the optimization algorithm.
3.3. The Second Optimization Based on P-CMA-ES
To improve the accuracy of the model, this part uses P-CMA-ES [
27] to optimize the attribute weight, rule weight, and GMF-B adjusted belief degree of DO-BRB-I. This method can dynamically adjust the search strategy to better adapt to the characteristics of the current optimization problem. However, the original P-CMA-ES algorithm takes the form of random scattering points and iterates continuously to find the optimal parameter [
20]. This type of random global optimization may result in optimized parameters that deviate excessively from the original parameters, erasing the original expert knowledge and making the model difficult to comprehend. Thus, the incorporation of interpretability constraint techniques improves the original P-CMA-ES method. The enhanced P-CMA-ES, the optimization method that emerged as a result, is best explained as follows:
The DO-BRB-I modeling accuracy is represented by the mean square error (MSE). In this section, the rule weights, attribute weights, and belief degrees are optimized parameters, so the MSE can be expressed as follows:
where
is the training data volume,
is the actual output value and
is the predicted output value of the system. Therefore, the optimization objective is to minimize the MSE through iterative refinement. The specific process is as follows:
Step 1 (Initial operation): Provide the initial parameters as follows:
where
is the optimized parameter set.
Step 2 (Sampling operation): The initial population can be determined by the following:
where
represents the
i-th solution of the (
g+1)-th generation,
represents the mean of the population,
represents the step size,
represents the normal distribution, and
represents the covariance matrix of the
g-th generation.
Step 3 (Constraint operation): In this step, the four constraint strategies introduced by Equations (12)–(16) are added to the optimization process of the model. These four constraint strategies not only retain the parameters of rules that have never been activated, but also impose detailed constraints on them. Ensuring that expert knowledge is not compromised.
Step 4 (Projection operation): The solutions produced by sampling operations may not satisfy the constraints, thus necessitating projection operations to ensure adherence to the constraints:
A hyperplane can be represented as , where represents the number of equality constraint variables in solution , represents the number of equality constraints in solution , and represents a parameter vector.
Step 5 (Selection operation): Update the mean value by performing selection operations using the following formula:
where
represents the weight coefficient of the
i-th equation.
represents the
i-th solution in the (
g+1)-th generation.
represents the subpopulation size.
Step 6 (Updating operation): Update the covariance matrix through adaptive operations and determine the population search range and direction. The calculation process is illustrated in the following formula:
where
and
represent the learning rates.
represents the evolutionary path of the covariance, and 0 is the initial evolutionary path.
To sum up, the overall optimization process of DO-BRB-I is shown in
Figure 3.
4. Reasoning Process of DO-BRB-I
The ER is an inference rule for evidence composition. The method is mainly used to solve uncertainty and fuzziness problems, and it uses the concepts of Dempster–Shafer theory and Bayesian reasoning. The essence of the ER algorithm is to integrate information from multiple independent pieces of evidence to form a comprehensive and reliable conclusion or prediction [
28]. This fusion improves the reliability and accuracy of decision-making and achieves good results, so ER is often used in BRB reasoning. The DO-BRB-I model developed in this study is derived from the fundamental BRB model. The ER has been selected as the inference engine for this model. Its detailed process is described as follows:
Step 1: Calculate the degree of matching, which indicates the flexibility of the rules, between the input sample information and belief rules. The matching degree of the
k-th rule for the
i-th input is calculated as follows:
where
represents the degree of matching,
represents the reference value of the prerequisite attribute, and
represents the input data. This membership function ensures that at least one rule can be activated for each input of data.
Step 2: Determine the activation weight through the following calculation:
where
represents the attribute weight for the
kth index.
Step 3: Calculate the ultimate belief degree by applying rule inference with the ER analysis algorithm using the following formula:
where
represents the belief degree from the final belief distribution.
Step 4: Calculate the expected utility value. The final belief distribution result and utility conversion formula are as follows:
where
is the actual input vector,
is the utility of
, and
is the expected utility at the end.
is a set composed of belief distributions. This methodology, which combines the IF–THEN rule-based and utility-based methods, enhances the reliability of initial information and facilitates logical adjustments to belief structures during the inference process.