Abstract

Capacitance aging of DC link capacitors in voltage source inverters (VSIs) is a common fault which can lead to instability of the DC voltage. In such a failure state, although the VSI can still work, its performance gradually deteriorates, resulting in a shorter service life of the equipment. Here, an online monitoring and fault diagnosis method for capacitance aging based on evidence reasoning (ER) rule is presented. Features from the DC link voltage data with different levels of capacitance aging are extracted, and data features are generated as pieces of diagnostic evidence, which are then combined according to the ER rule. Finally, capacitance aging fault levels were estimated using the combined results. This method has better diagnostic performance compared to a backpropagation (BP) neural network approach and can be used to flexibly define the relative weighting of each evidence parameter depending on the application. This approach can therefore be widely used for fault diagnosis of an array of different devices.

1. Introduction

Voltage source inverters (VSIs) are generally reliable, accurate, and effective and have become the most common inverter type used in industry [1]. A capacitor is a key component of VSIs, as they supply the input current to the inverter, compensate for the difference between the power requirement of the inverter and the output power of the rectifier, suppress the current harmonics, absorb feedback energy, meet instantaneous power demand, and stabilize voltage fluctuations [2]. There are three main types of capacitor: electrolytic, ceramic, and film capacitors. Of these, the electrolytic capacitor is most commonly used in VSIs. Electrolytic capacitors generally have a much shorter life expectancy compared with the service life of the inverter. Capacitor faults generally include short circuits, electric leakages, and capacitance attenuation. As inverter operation proceeds, evaporation of electrolyte will cause gradual capacitance attenuation. And although the inverter can still function as capacitance aging increases, overcurrent and under voltage situations increasingly arise when the inverter runs under load, increasing potential safety hazards [3]. This fault is classified as an incipient fault and does not generally trigger the protection system. However, if left undetected and untreated, this situation can significantly reduce the service life of the equipment. It is therefore increasingly important to diagnose capacitance aging in order to evaluate and carry out preventative maintenance in advance.

Historically, equivalent series resistance (ESR) and capacitance C are two important parameters used to evaluate capacitance aging. A fault situation is defined if the ESR is more than the double of its initial value. When capacitance attenuation exceeds 20% of the initial value, the capacitor is classified as failed [4]. The majority of methods used to detect capacitor faults are based on monitoring the ESR value, including offline and online modes. Offline detection methods are gradually becoming obsolete as they do not allow dynamic management [5]. On the other hand, for online detection methods, the most direct and effective method is to estimate the ESR value using the correlation between current and voltage values [6, 7]. The ratio of capacitor voltage to supply voltage can also be used to calculate the ESR by reference resistor [8]. However, the accuracy of the technique depends largely on the difference of the reference resistors. The influence of temperature on ESR has also been widely investigated [9, 10]. Improvements in the measuring terminal have been used to obtain more accurate ESR values [11]. In addition, a built-in self-test (BIST) method was used to monitor the attenuation of capacitance C [12]. However, there are few research studies which consider the capacitor fault state jointly determined by ESR value and capacitance C. Oukaour et al. [13] monitored both parameters, and the aging index was estimated using the least squares method, which is proved to be effective.

However, all above methods are based on the principle of using formulae to fit the relationship between the acquisition signal source and the capacitor ESR value and capacitance C which is not ideal. The monitoring parameters are various and complex. Moreover, it is difficult to identify incipient faults and predict serious faults using this approach. At present, data-driven fault monitoring and identification methods are increasingly used due to their broad applicability and high diagnostic accuracy [14, 15]. Data-driven methods using artificial neural networks (ANNs) have been used to estimate capacitance capacity [16]; however, the accuracy of ANN method is limited by the quantity of training samples. Extended neural network (ENN) algorithms and chaotic synchronization detection methods have also been used to identify capacitor faults [17]. An online track circuit compensation capacitor fault diagnosis method based on K- fault diagnosis method was reported [18]. Unfortunately, these methods only consider the identification of serious faults that have occurred and do not carry out further subdivision on the degree of capacitor faults which would allow much earlier intervention and fault correction.

Here, we propose a fault diagnosis method for capacitance aging in DC link capacitors used in inverters which is based on the evidence reasoning (ER) rule [19]. Based on the principle of data-driven decision making, this method does not assume any relationship between capacitance and other parameters [20, 21], but only generates types of diagnostic evidence and combines them by measuring the DC link voltage value under different capacitance attenuation degrees and then estimates the fault level of capacitance aging based on the combined results. The advantage of using the ER rule method is that it has the ability to deal with the variable information provided by incipient fault situations and it can also define the relative reliability of different information sources. Moreover, for multiple types of evidence, the relative weighting or importance between different types of evidences can be allocated. This weighted information can then be considered comprehensively to provide more correct data for final decision making. Different fault thresholds can be set according to expert user experience in different practical applications. The method also has the advantages of using fewer parameters but providing higher monitoring sensitivity.

2. Simulation Model and Experiment

2.1. A Simulation Model for a Marine Electric Propulsion System

In marine electric propulsion systems, the propeller is driven by an electric motor. This propulsion method abandons the traditional way in which the shaft system and propeller are driven by the prime mover and is efficient and energy saving and allows excellent maneuverability and flexible arrangements [22, 23]. The most common electric propulsion system is generally powered by diesel generator sets, which provide electric power to propulsion motors through switchboard and frequency converters and which then drives the propellers to generate thrust. The ship “Jiangsu ferry 3011” was used in this work as a case study and is a ferry ship which is in use and which runs across the Yangtze River between Zhenjiang and Yangzhou in Jiangsu province in China (Figure 1). The ship is 90 meters long and 15.4 meters wide and has a tonnage of 720 tons and is the first civil ship adopting an electric propulsion system (Figure 2) using DC power in China.

A DC power system is used on board, which means the electric power is generated from an alternating current (AC) generator and transformed into DC power by a rectifier to feed the DC bus. Power is then transmitted to inverters to drive propulsion motors and other auxiliary loads. The basic parameters of this marine power system are shown in Table 1.

The overall life cycle of electric propulsion ships is characterized by working in highly variable environments with multiple types of interference, strong maneuvering requirements, and complex river conditions for 24 hours per day, always accompanied by potential incipient faults in the electrical system. The inverter is the core component of the ship’s electric propulsion system. The failure of the DC link capacitor is one of the most common faults of the inverter [24]. Due to the different operating environments and different types of electric propulsion systems, the degree and characteristics of capacitor attenuation are irregular and uncertain. In addition, with current monitoring technology, it is difficult to track, collect, and extract data which can identify fault features. The accumulated data are often incomplete and invalid, which frustrates careful research-based and data-driven solutions. However, by building a computer model of this system, we can simulate the life cycle characteristics of the ship’s electric propulsion system and track the fault characteristics of the capacitor with time, laying the foundations for a new and rigorous research-based approach to incipient fault diagnosis which can be developed for use on board these electric ferries.

Based on the main parameters in Table 1, we developed a simulation model of this marine electric propulsion system using a Matlab/Simulink environment [25] as shown in Figure 3.

In the simulation model, the Jiangsu ferry 3011’s electric propulsion system consists of generator, rectifier, inverter, and propulsion motors. The DC link capacitor is made up of two identical capacitors in series, and N is the midpoint between the two capacitors. The propulsion system consists of two propulsion motors of the same type. Tm1 and Tm2 are the load torques of the two propulsion motors.

2.2. Fault Simulation

In this paper, the most frequent working condition was selected for research, and the load torque was set at 2228 Nm. In the DC link, three different levels of capacitance attenuation were set in the simulation model. The initial capacitance C1 is 7500e − 6 F, and 90∼100% of the initial capacitance is classified as normal. With the gradual deterioration of capacitance, if the capacitance is set to C2 (80∼90% of initial capacitance), it is defined as an incipient fault. If the capacitance value is set to C3 (less than 80% of initial capacitance), it is no longer an incipient fault, but a serious fault. At the same time, the data sampling frequency was set to 50 kHz to collect the instantaneous DC link voltage VDC for 24 seconds, including three levels of capacitance with 1.2 million data points. Among them, the data of C1 are generated in the1∼8th second, the data of C2 are generated in the 8∼16th second, and the data of C3 are generated in the 16∼24th second. In order to make the diagnosis process closer to the actual situation, the measured VDC is added with noise, and the awgn function in Matlab was used to add 30 dB white Gaussian noise to VDC as shown in Figure 4.

In this figure, The fluctuation range of the upper and lower limits of the DC voltage under the capacitance C1 is from −0.79% to +0.49%; for capacitance C2, the fluctuation range is from −0.87% to +0.56%, and for capacitance C3, the fluctuation range is from −1.09% to +0.65%. So, the instantaneous value of the DC bus voltage under different fault levels has no obvious feature distinction, and it is difficult to effectively diagnose using the general threshold method.

We further propose an inference model based on an ER rule to estimate the fault level of capacitance attenuation by collecting and analyzing instantaneous voltage data VDC from the DC link. The definition of fault levels and their thresholds can be set according to different capacitor types and different application systems. This method does not need to make assumptions about the relationship between ESR values, capacitance values, and voltage values from the DC link as it directly learns the fault characteristics and can then be used to diagnose a predetermined fault level using a data-driven approach.

3. Fault Diagnosis Based on Evidence Reasoning Rule

3.1. Outline of ER Rule

The evidence reasoning (ER) rule has been proposed by Yang et al. [26, 27], and the evidence can be defined here as follows:

Here, define as a frame of discernment; it contains mutually exclusive and collectively exhaustive hypotheses. and represent all the subsets of the power set . is an element of evidence ej which indicates the degree of ej approaching (not an empty set) is , when , is the focal element of ej.

The concepts of reliability rj and weight of ej are introduced in the ER rule. rj is the ability of the information source, which characterizes credibility of this ej for the fault diagnosis system. On the other hand, for , it is different from rj, which can be defined by the person who uses this diagnostic method according to subjective conditions such as different equipment or different environments. represents the importance degree of different types of evidence. The weighted belief distribution with reliability factor and weight factor is defined as follows:where represents the degree of ej approaching after combining rj and , which is defined as follows:where , are normalization factors; it makes , wherein .

If the two pieces of evidence are independent of each other, they can be combined with the ER rule as follows:where represents the degree of how the two pieces of evidence jointly support , and the fusion process includes the bounded sum and the orthogonal sum. The bounded sum is reflected in , which means the sum of the two sets of evidence e1 and e2 independently approaches , and it includes the interaction between e1 and e2. For example, when e1 is considered to be completely reliable , then , which means that e2 has no effect on the combined results. Conversely, e1 has no effect on the combined results. In addition, the orthogonal sum indicates the degree of all evidence conjointly approaching to. This method can also combine multiple pieces of evidence without order restriction while satisfying the exchange law and the combination law [28, 29].

3.2. Diagnostic Process

Figure 5 outlines the fault diagnosis process using the ER rule.

First, the instantaneous DC voltage value VDC is obtained as shown in Section 2.2 so that the peak-to-peak voltage value VPP can be calculated from VDC. Here, 0.02 seconds, including 1000 data points, is taken as a window length to obtain the average value of VPP of 400 data windows named as VP(t). VR(t) is obtained by calculating the root mean square value of the VDC with same window length, and then by normalization. and , where , , , and are, respectively, the minimum and maximum values of input characteristic signals VP(t) and VR(t). The capacitor aging fault level is then defined as C(t), C(t) = [1, 2, 3]. And VP(t), VR(t), and C(t) are denoted together as T = {[VP(t), VR(t), C(t)] | t = 1, 2, …, TS}, TS = 1200, where [VP(t), VR(t), C(t)] is a sample vector.

Secondly, the data [VP(t), VR(t), C(t)] in T are converted to the form of similarity degree to the reference value. The reference value set of the capacitor fault level is D = {Dn | n = 1, …, N}, where N is the number of reference values of the capacitor fault level. The input reference values set of the voltage signal Vi is Ai = { | j = 1, …, Ji}, i = P, R, where Ji is the number of reference values of Vi. The similarity distribution of Vi(t) to reference values is calculated as follows:where represents the similarity of Vi(t) to the reference value .

Similarly, the calculation process for the similarity distribution of the fault level C(t) to the reference value Dn is as follows:where denotes the similarity of C(t) to Dn

The parameter form similar to can be obtained after completing the calculation of all the sample sets in T, where denotes the comprehensive similarity of the input value in sample pair (Vi(t), C(t)) matching the reference value , Dn. Table 2 shows the cast statistics result of all sample pairs in T, and .

According to the casting statistics shown in Table 2, the degree of belief of ej can be calculated as follows:where . Then, the evidence corresponding to the reference value can be defined as follows:

Therefore, an evidence matrix, as shown in Table 3 can be constructed to describe the relationship between Vi(t) and C(t).

Now we define reliability rj and weight of evidence ej, after obtaining the evidence of the input Vi. ri defined in this paper is determined by the Spearman rank correlation coefficient for ranked data [30]. The Spearman rank correlation coefficient expresses the closeness of the relationship between the two sets of variables and it is suitable for describing the reliability between Vi(t) and C(t) in this case. The calculation procedures are as follows:where is the difference between Vi(t) and C(t) and is the reliability of the evidence.

For all of the input values Vi(t), they must be within the interval formed by the two adjacent reference values, and evidence and corresponding to the two reference values will be activated. The final evidence of the input values Vi(t) can be weighted by and with the calculated reliability .

The initial evidence weight is set as  = ri, and e1 and e2 are combined by ER rule. The fusion process is as follows:

After the fusion process is completed, the combined result is obtained. Furthermore, the fault level of the capacitor can then be estimated based on the combined results.

4. Analysis of Simulation Examples

4.1. Training Model

In the experiment, 80% of the total sample set T (sample size TS = 1200) was randomly selected as the training sample TR = {[VP(t), VR(t), C(t)] | t = 1, 2, …, TSR}, sample size TSR = 960, and the remaining 20% was selected as the test sample TE = {[VP(t), VR(t), C(t)]| t = 1, 2, …, TSE}, sample size TSE = 240. Based on experience, the reference value set for the initial parameter VP(t) is A1 = {6.7877, 7.8030, 8.1517, 8.4575, 8.8254, 9.1441, 9.5400, 9.9043, 10.3445, 10.7758}, including 10 reference values in total. The reference value set for the initial parameter VR(t) is A2 = {0, 0.0679, 0.1245, 0.1733, 0.2177, 0.2759, 0.3375, 0.4224, 0.5122, 0.6417}, including 10 reference values in total. The reference value set D of C(t) is D = {1, 2, 3}, including 3 reference values in total. The training sample set TR is processed through the steps described in Section 4 to obtain casting result tables of sample pairs (VP(t), C(t)), which are presented in Tables 4 and 5.

After obtaining the casting result tables, the evidence matrix table of Vi is obtained by normalization of likelihoods as shown in Tables 6 and 7.

After obtaining the evidence matrix tables, the reliability factors r1 = 0.8897 and r2 = 0.4117 of the evidence parameters e1 and e2 are calculated with the sample dataset T. The initial fault diagnosis result on the sample training set TR can be obtained by using the ER rule, and the results are presented in the form of a confusion matrix as shown in Table 8.

It can be seen from the initial diagnostic results that the accuracy is not satisfactory because the initial set of reference values Ai, reliability factor rj, and weight are not the optimal ER rule parameters, so the ideal diagnostic accuracy cannot be achieved. Therefore, it is necessary to establish a parameter optimization model for the ER rule process.

4.2. Parameter Optimization

The mean square error between the estimated value and the real value of the capacitor fault level is set as the objective function, and the set of optimization parameters P is determined at the same time. P = {, | i = P, R; j = 2,..., Ji − 1}, represents the weight of the evidence, and the other parameters are set to , , , .

The objective function is optimized by using the fmincon function in Matlab, and the optimal parameter set P is obtained after training. At the same time, the optimized casting result tables, optimized evidence matrix table, and optimized capacitance aging fault diagnosis results can be obtained.

4.3. Analysis of Diagnostic Results

As can be seen from Tables 413, the optimized diagnostic results are more accurate. Especially for capacitor fault levels 2 and 3, the judgment is more accurate, and misjudgment is greatly reduced. In order to further verify the reliability of the ER fault diagnosis model, test set TE was used for diagnosis. In the meanwhile, a classical three-layer BP neural network approach was selected for comparison. The learning rate was set as Ir = 0.05, and the maximum number of iterations was set to 100. The diagnostic results of TE are shown in Table 14 and Figures 6 and 7.

The difference between the true and estimated levels of capacitance at different level capacitance attenuation in test set is as shown in Figure 7. In this figure, the ordinate value of 1 indicates that the capacitance is normal, the value 2 indicates that the capacitance is in incipient fault, and the value 3 indicates that the capacitance is in serious fault. On the time axis of the abscissa, the true level of the capacitance is normal for 0 to 8 seconds, an incipient fault for 8 to 16, and a serious fault for 16 to 24.

Figure 8 shows the estimated levels of the test set calculated by the trained ER and BP neural network, respectively.

From the results presented, the diagnostic accuracy of the BP neural network approach is lower than that obtained using the ER rule model, which is mainly seen in the diagnosis of capacitor fault level 2. This fault level 2 represents an incipient fault. For the case of C2, the change of DC bus voltage is not as obvious as case of C3, which is a very small change. Therefore, it is difficult to distinguish the change from C1 to C2 in the specific diagnosis process. In this case, there are two reasons for the low accuracy of BP neural network. Firstly, BP neural network is sensitive to initial weights and bias, and these two values are randomly selected in each BP training, so every time, the training tends to converge to different local minimum, resulting in low diagnostic accuracy in some diagnosis experiments. Secondly, the choice of structure of BP neural network has not a unified and complete theoretical guidance, normally only selected by experience. The selected structure may not the most suitable, and there may be a better network structure, but it takes more time to select and adjust. However, the above problem does not exist in ER rule; its initial reliability and weights are derived from sample data, and their calculating method has clear physical meanings. The evidences generated method is based on data statistics and likelihood normalization method. It is strictly based on the principle of sample statistical reasoning, thus reducing the randomness and uncertainty in process of diagnosis.

In the process of actual dynamic monitoring, fault trend prediction is of great importance to marine electric propulsion systems. Nowadays, what the marine electric propulsion system operators need urgently is to be able to predict the state of the equipment before faults occur to allow the faults to be corrected prior to more serious damage or economic loss. Therefore, the ER rule method is a flexible and applicable approach which serves as an early warning and management assistance system for fault diagnosis in marine electric propulsion systems.

5. Conclusions

In this paper, an online monitoring and fault diagnosis method for capacitance aging of DC link capacitors of VSI based on ER rule is proposed. Diagnostic evidence is generated from information sources, and then reliability factors and weighting of each parameter are determined by evaluating the correlation between fault features and fault levels. The ER rule was then used to combine the evidence, and each fault level was estimated based on the combined results. Finally, we defined the objective function to optimize the ER rule parameters to make the diagnostic results more accurate. Compared with most classical data-driven method, BP neural networks, the ER rule method also has higher diagnostic accuracy and applicability in the field of incipient fault diagnosis. Unlike other types of fault diagnosis methods, this data-driven method has the advantage that it does not need to assume any mapping relationships between fault features and fault models. In addition, this method can be extended not only to serious fault diagnosis but also incipient fault diagnosis due to its ability to predefine different weighting to each evidence parameter depending on the application field. When using this method on ships, it is flexible enough to allow changes to the parameters in the ER rule to meet the needs of the different electric propulsion systems in use on board diverse types of vessel. The Jiangsu ferry 3011 is a new ship. However, the aging process of capacitor will take years. Therefore, the actual aging data cannot be obtained for research. The method described in this paper will be a promising diagnostic method to detect incipient fault for Jiangsu ferry 3011; further analysis and discussion will be carried out after collecting overall life cycle’s physical data.

Data Availability

The data used to support the findings of this study were supplied by [The 711th Research Institute of China Shipbuilding Industry Corporation] under license and so cannot be made freely avaliable. Requests for access to these data should be made to [Yelan He, email:heyl99@163.com].

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Linhao Liao and Haibo Gao contributed equally to modeling, calculation, and writing. Yelan He, Xiaobin Xu, Zhiguo Lin, Yajie Chen, and Fubing You gave useful advice and suggestions and contributed to the discussion and preparation of the manuscript.

Acknowledgments

The authors gratefully acknowledge the financial support provided by the NSFC-Zhejiang Joint Fund for the Integration of Industrialization and Informatization (Project no. U1709215).