Research Article  Open Access
Yifan Jian, Xianguo Qing, Yang Zhao, Liang He, Xiao Qi, "Application of ModelBased Deep Learning Algorithm in Fault Diagnosis of Coal Mills", Mathematical Problems in Engineering, vol. 2020, Article ID 3753274, 14 pages, 2020. https://doi.org/10.1155/2020/3753274
Application of ModelBased Deep Learning Algorithm in Fault Diagnosis of Coal Mills
Abstract
The coal mill is one of the important auxiliary engines in the coalfired power station. Its operation status is directly related to the safe and steady operation of the units. In this paper, a modelbased deep learning algorithm for fault diagnosis is proposed to effectively detect the operation state of coal mills. Based on the system mechanism model of coal mills, massive fault data are obtained by analyzing and simulating the different types of faults. Then, stacked autoencoders (SAEs) are established by combining the said data with the deep learning algorithm. The SAE model is trained by the fault data, which provide it with the learning and identification capability of the characteristics of faults. According to the simulation results, the accuracy of fault diagnosis of coal mills based on SAE is high at 98.97%. Finally, the proposed SAEs can well detect the fault in coal mills and generate the warnings in advance.
1. Introduction
The coal mill is one of the important auxiliary equipment of coalfired units, and its operating status is directly related to the safe and stable operation of the units. When a fault occurs in the coal mill, the fuel supply of the boiler cannot be guaranteed which creates the mismatch between boiler energy output and the turbine power output. Under this situation, a quick load rejection operation will occur, which directly leads to fire extinguishing in the furnace. The fault in the coal mill will cause large economic loss to power generation enterprises and decrease the safety and stability of the power system. Therefore, it is of great necessity to guarantee the normal operation through effective fault warning and diagnosing of coal mill.
Agrawal et al. [1] divided the fault diagnosis methods into three categories: model, signal, and historical operation databased fault diagnosis methods. Modelbased fault diagnosis methods need to establish the mathematical model of the coal mill. Odgaard and Mataji [2] used a simplified energy balance equation to monitor and diagnose abnormal energy flow in the coal mill. Andersen et al. [3] designed a Kalman filter to estimate the moisture in the coal that enters and exists a coal mill to determine whether the energy in the coal mill is in normal condition. Based on the multisegment model of coal mills established by Wei et al. [4], Guo et al. [5] realized the monitoring of the state of coal mills by identifying the abnormal variation in the model parameters. Modelbased fault diagnosis methods analyze the mathematical model of the actual object for fault diagnosis, and thus, the physical meaning is clear. However, establishing the exact model in practical application is difficult [6–10]. Thus, the operability of these methods is poor.
Signalbased fault diagnosis systems are widely used to evaluate the health of mechanical equipment. Many signals of high frequency change during the operation of the coal mill, such as current of coal mill, outlet primary air flow of coal mill, and differential pressure of primary air. Su et al. [11] designed a system that records the vibration signals of the coal mill and shifts them to energy amplitudes by use of wavelet analysis. Whether the coal mill is in coal interruption or coal choking, other fault operations can be determined by analyzing the relationship between the vibration signals and the amount of coal in the mill. Kisić et al. [12] proposed a method to detect the wear degree of grinding roller and analyzed the multivariate control chart on the frequency spectrum to find the appropriate time to replace the worn parts. Collura et al. [13] utilized model identification and signal processing techniques to develop a coal mill performance monitoring tool based on realtime detection of the fineness of pulverized coal. Compared with the modelbased fault diagnosis method, the signalbased fault diagnosis method does not need to establish complex object model. Only through the analysis of collected data can a fault in the system be found. However, these methods often need to install a large number of sensors to collect signal, thereby resulting in high implementation and maintenance costs.
Fault diagnosis based on historical operation data is mainly done by analyzing the differences between the normal operation data and fault operation data to determine the health status of coal mill. Han and Jiang [14] proposed a fault diagnosis method based on fuzzy decision clustering and used a singlelayer neural network to realize three kinds of fault identification of coal mill. Qin et al. [15] utilized the abnormal operation data of the coal mill to establish an expert system to determine the operating status of the coal mill by comparing the trend of the model output with the expert system. A databased fault analysis method is a datadriven approach, and even researchers who are unfamiliar with the system can use relevant algorithms for analysis. However, fault types and fault data in the mass historical data of the thermal power units are incomplete and a datadriven method requires analyzing a large amount of fault data. Thus, selecting the fault data from the vast amount of historical data one by one is difficult [16–18].
The modelbased fault diagnosis method needs to establish an accurate model of the coal mill in order to obtain good fault diagnosis results. However, the coal mill is a complex object with multiparameter coupling. It is difficult to establish an accurate mathematical model. The premise of applying signalbased fault diagnosis methods is to be able to measure the monitored parameters. Therefore, it is necessary to install a large number of new sensors on the shell of the coal mill. However, when the coal mill was initially constructed, it usually did not consider reserving the mechanical interface for new sensors. So, it is not easy to install new sensors on the shell of the coal mill. The fault diagnosis method based on historical operation data firstly needs to obtain a large amount of fault operation data of the coal mill. However, the fault data of the coal mill is usually mixed with the normal operation data, which is difficult to classify and identify. Based on the above analysis, the existing methods are difficult to achieve good application results for the fault diagnosis of coal mills. Although the above three types of traditional methods have shortcomings, combining their advantages can find a simpler and more effective method to solve the fault diagnosis of coal mills. The basic idea is to obtain fault simulation data based on a simplified model and use big data analysis for fault identification. In recent years, the rapid development of deep learning algorithms has provided the possibility of big data analysis. Guo et al. [19] constructed an adaptive convolution neural network, which greatly improves accuracy of fault diagnosis of motor bearing. Duan et al. [20] used a deep learning algorithm to study the missing traffic data to implement the interpolation of missing data. In this study, a modelbased datadriven fault diagnosis method is proposed to obtain a fault diagnosis method with simple operation, low cost, and high accuracy. First, on the basis of the simplified coal mill model, massive fault data of the coal mill are obtained by analyzing the fault principle and simulating the fault operation status of the coal mill. This method solves the difficulty in obtaining a large amount of fault data manually from the massive data. Then, stacked autoencoders (SAEs) with multilayer neural networks are established on the basis of the theory of deep learning algorithm. The numerous fault data obtained by the steps above are used in the training of the networks to fully motivate the nonlinear characteristics of deep neural networks. Thus, the built network can accurately learn the essential characteristics of all kinds of faults and then achieve the early warning and diagnosis of the fault in coal mills. The above method can greatly improve the fault diagnosis accuracy of the coal mill, and at the same time, it can also provide fault warning to the operator, which is of great significance for ensuring the safe operation of the power plant and ensuring the safety of the equipment.
The rest of the paper is organized as follows. Section 2 introduces the working principle of the coal mill and its nonlinear dynamic model. Section 3 analyzes the mechanism of two typical coal mill faults and obtains a large number of fault data by simulation experiments. Section 4 introduces the working principle of the SAEs and makes certain improvements to the model. Section 5 is the simulation analysis that aims to verify the effectiveness of the proposed method in the fault diagnosis of the coal mill. Section 6 elaborates the conclusions of the study.
2. Brief Introduction of the Coal Mill System
2.1. Working Principle
MPStype mediumspeed coal mill [21] is a rollertype coal mill designed and manufactured by Babcock, Germany. Such mills are characterized by smooth output, low energy consumption, and long maintenance period. In this study, MPS180HPII mediumspeed coal mill is used in the analysis. The maximum output is 44.496 t/h, and the fineness of coal powder R90 is 22% (Figure 1). R90 indicates the probability that coal powders cannot pass through a sieve with a pore size of 90 μm.
The raw coal falls into the coal mill through the coal dropping pipe and is milled into coal powder under the squeezing effect of two milling parts (grinding disks and rollers) [4]. The primary air enters the coal mill through the annulus around the grinding disk to dry the coal powders and bring them into the coarse coal separator for separation. The qualified fine coal powders are blown into the boiler for combustion while large ones return into the coal for subsequent milling.
2.2. Mathematical Model of the Coal Mill
The operation of the coal mill involves the mass balance of coal and the energy balance of the entire coal mill. Establishing an effective dynamic mathematical model of coal mills is an important prerequisite for the state monitoring of coal mills. Zeng et al. [22–24] established an MPS mediumspeed mill model (equation (1)), which includes three inputs and three outputs based on the mass and energy balance of the primary air and coal moisture in the mill. The proposed method is based on this model, and the symbolic description of the model is shown in Table 1:

In Equation (1), U_{L}, U_{H}, and are the control quantities of the model; W_{air}, T_{out}, and W_{pf} are the output quantities of the model; and K_{i} and T_{i} (i = 1, 2, …, 15, j = 1, 2) are the model parameters to be identified and of which the values are shown in Table 2.

3. ModelBased Coal Mill Fault Simulation
The fault types and fault data in the vast amount of historical data of the thermal power units are incomplete, and selecting the fault data one by one from the massive historical data is difficult. Therefore, effectively obtaining a large number of fault data is the key to solve the fault diagnosis of coal mill. The simulation results in [22] showed that the mathematical model of MPStype medium speed coal mill presents high precision. In the current study, the coal mill model is used in the analysis and two typical coal mill faults (coal interruption and coal choking) are simulated by analyzing the fault mechanism of coal mill. The simulation experiments obtain a large number of fault data, which can effectively solve the difficulty in obtaining fault data manually from the massive data.
First, a control scheme is designed for the coal mill model. The purpose is to ensure that the simulation experiments are conducted under the closed loop regulation, such that the fault data obtained by the simulation experiments can be significantly close to the real operation status of coal mill. The control scheme is shown in Figure 2. The entire control scheme consists of three controlled, three control, and four state variables as presented in Table 3. The control circuits are composed of three singleloop proportionintegraldifferential (PID) controllers. Specifically, PID1, where the setting parameters, respectively, are K_{p} = 1, K_{i} = 0.05, K_{d} = 0, controls the outlet temperature of coal mill by adjusting the valve position of cold air. PID2, where the setting parameters, respectively, are K_{p} = 2, K_{i} = 0.5, K_{d} = 0, controls the inlet primary air flow of coal mill by adjusting the valve position of hot air, and PID3, where the setting parameters, respectively, are K_{p} = 0.1, K_{i} = 0.1, K_{d} = 0, controls the outlet pulverized coal flow of coal mill by adjusting the inlet coal flow of coal mill. After designing the control scheme, the fault operation status of the coal mill can be simulated by adjusting the corresponding controllers.

3.1. Fault Simulation of Coal Interruption
When an obstruction exists in the coal dropping pipe or a fault occurs in the coal feeder, the amount of coal into the coal mill will reduce directly and coal interruption will occur when the case is serious, thereby endangering the stability of the boiler combustion. The process of simulating coal interruption is as follows. When the coal mill is in stable operation, a negative slope signal is superimposed on PID3, such that the mass of coal entering the coal mill is gradually reduced to 0. The data generated in this process can be considered the coal interruption samples. A large number of coal interruption samples can be obtained by adjusting the set value to run the coal mill in other operation status and repeating the steps above to record fault data.
To verify the effectiveness of the simulation experiments of coal interruption, the variables that change significantly during the period of coal interruption are selected and their varying curves are drawn. Figure 3 shows the result of an arbitrary selection of experimental data. Figure 3(a) shows that coal interruption decreases the mass of coal entering the coal mill, which then decreases the outlet pulverized coal flow of the coal mill. Meanwhile, the heat consumption of the inlet primary air flow of coal mill through the coal mill reduces, thereby resulting in an upward trend in the outlet temperature of coal mill; the valve position of cold air is then rapidly opened, thereby making the outlet temperature of coal mill fall (Figure 3(b)). The reduction in the mass of coal stored in the coal mill results in the reduction of the current of coal mill and the differential pressure of primary air. The trend is consistent with that of the curves described in Figures 3(c) and 3(d). The ramp signal is removed after 75 s, and the variables are returned to the original set value under the control of the controller.
The research object in [19] is a MPStype medium speed coal mill in a power plant in Hainan, China. The current study obtains sets of fault data of coal interruption by looking for the historical operation data of the coal mill and draws the varying curves of key variables as shown in Figure 4. According to the accident analysis, the coal interruption fault occurs because of the malfunctioning of the coal feeder; as a result, the actual supply of coal gradually reduces to 0 (Figure 4(a)). Figure 4 shows that, when coal interruption fault occurs, the outlet temperature of the coal mill rises (Figure 4(b)); however, the current of coal mill (Figure 4(c)) and the differential pressure of primary air (Figure 4(d)) decrease. This changing trend is similar to that of the key variables in the simulation of coal interruption. Therefore, the simulation of coal interruption in this study is reasonable. Accordingly, the data in the rectangular frame in Figure 3 can be recorded as fault samples.
(a)
(b)
(c)
(d)
3.2. Fault Simulation of Coal Choking
Coal choking may be caused by too little inlet primary air flow of coal mill, excessive coal feed, or too much moisture in raw coal. The process of simulating coal choking is as follows. A positive step signal is superimposed on PID3 to make the mass of coal in coal mill quickly reach the upper limit, and the data are recorded as the coal choking samples. Similar to previous simulation experiments, a large number of coal choking samples can be obtained by adjusting the set value to run the coal mill in other operation status and repeating the steps above to record fault data. The upper limit of the mass of coal stored in the coal mill is set as 60 kg.
The varying curves of the variables are drawn, and Figure 5 shows the result of an arbitrary selection of experimental data. As shown in the figure, the sudden increase in the set value of the outlet pulverized coal flow of coal mill causes the mass of raw coal in the coal mill to rise continuously (Figure 5(a)) and increases the resistance along the way, which results in the increase in differential pressure of primary air (Figure 5(c)). At the same time, the work load of coal mill increases accordingly, such that the current of the coal mill increases as well (Figure 5(b)). When the thickness of the raw coal reaches a certain degree, the grinding efficiency drops significantly, which then reduces the current of the coal mill. The outlet pulverized coal flow of coal mill is reduced to 0 until the mass of raw coal in the coal mill reaches the upper limit, at which time the primary air pipe is blocked and the pulverized coal cannot be blown out (Figure 5(d)). The analysis above shows that the simulation experiment results are consistent with the fault characteristics of coal choking; therefore, the data in the rectangular frame can be used as fault samples.
4. Stacked Autoencoders
The fault diagnosis method based on historical operation data is a datadriven approach, which aims to obtain the nonlinear mapping relationship between the data and fault features. When sufficient data are available for learning, the deep neural networks can theoretically approximate any nonlinear function. This section describes a deep neural network called SAE, which is stacked by autoencoders (AEs), for fault diagnosis of coal mill and proposes two ways to improve the network performance.
4.1. Fundamentals of Autoencoder
An AE neural network can be considered a threelayer neural network. This network applies unsupervised learning algorithm to train and adjust the network weight and ultimately sets the network output to be equal to the network input. A typical example is shown in Figure 6, where {x_{1}, x_{2}, …, x_{n}; x_{i} ∈ R^{n}} can be treated as a set of unlabeled raw data and represents the network output. The circles with b are called bias units and correspond to the intercept term.
The transfer process of raw data from the input layer to the hidden layer is called encoding, and the transfer process from the hidden layer to the output layer is called decoding, which can be described bywhere S(·) represents sigmoid function, W_{1} represents the weight matrix between the input and hidden layers, W_{2} represents the weight matrix between the hidden and output layers, and b_{1} and b_{2} represent the bias.
According to the concepts mentioned in this section, AE tries to learn a function h_{w,b}(x)≈x. In other words, AE is trained to learn an approximate function such that the network output is similar to the network input. In fact, by putting constraints into AE, such as limiting the number of nodes in hidden layer, AE can obtain the lowdimensional feature of the data by compressing the highdimensional input data. In Section 3, a large number of fault data have been obtained by fault simulation experiments. The remaining parts focus on the establishment of a suitable AE network to find the relationship between data and fault characteristics by effective learning of the complex fault data.
The backpropagation algorithm is used for AE training. A training set of m training samples is assumed. The network can be trained using batch gradient descent. For a single training example , the cost function can be defined as
For a training set of m samples, the overall cost function iswhere represents weight decay coefficient that controls the relative importance of the two terms in equation (5). represents the synaptic weight between the ith neuron in layer l and jth neuron in layer l + 1. represents the number of layers in AE. In other words, can represent the output layer of the network, and represents the number of the total neurons in layer l. The first term in the definition of is an average sumofsquares error term. The second term is a weight decay term that can decrease the magnitude of the weights and prevent overfitting.
The weight W and bias b are updated with gradient descent as follows:where represents the learning rate. The partial derivatives in the equations above are derived as follows:wherewhere represents the activation of unit j in layer l and represents the error term of layer , given bywhere represents the input weighted sum of unit i in layer l and represents partial deflection of sigmoid function. The error term of the output layer n_{l} is given bywhere represents the activation of unit i of layer l, represents the activation of unit i in the output layer, and represents the input weighted sum of unit i in the output layer.
Repeating the above equations can make the output of AE equal to the input of AE by minimizing the overall cost function (equation (5)).
4.2. Improvement of AE
As mentioned in Section 4.1, limiting the number of nodes in hidden layer is conducive to helping AE learn the relationship between the input data and fault features, because reducing the number of neurons can simplify the structure of the hidden layer and reduce the dimension of the input data.
Restricting the number of neurons can reduce the dimension of the data, but the network can learn few features in the hidden layer. On the basis of guaranteeing the diversity of the features in the hidden layer, a method called sparse constraint is introduced in this study to improve AE. The main idea is not to reduce the number of neurons but to consider restrictions to limit the activities of the neurons and thus reduce the dimension of the input data. Accordingly, the original overall cost function (equation (5)) should be modified to introduce an additional penalty factor, given bywherewhere represents the sparsity penalty term, controls the weight of the sparsity penalty term, represents the average activation of unit j in hidden layer, represents a sparsity parameter, and s represents the number of units in one hidden layer.
The penalty term has the following property: if = , then ; the value increases monotonically with the difference between and . Therefore, the activations of hidden units are sufficiently small when is set close to zero.
Random noise is introduced into the input data to make the network learn rich information and thus prevent the AE from learning only the equivalent representation of the original data. The main idea is to set a small number of nodes in the input layer to zero at a small probability. However, the probability of introducing random noise should be appropriate; otherwise, the noise may cause irreversible damage to the input data.
SAEs are deep neural networks consisting of multiple layers of the improved AEs in which the output of each layer is wired to the input of the next layer. The SAE model is connected with a Softmax classifier to complete the construction of the deep neural network (Figure 7). The SAE model can identify the fault in the coal mill by learning the labeled data obtained from the fault simulation experiments of the coal mill.
5. Fault Diagnosis Based on SAE
5.1. Data Preprocessing and Health State Definition
In accordance with the fault simulation method described in Section 3, the simulation experiments are conducted repeatedly, and then 5000 sets of experimental data are obtained including three kinds of data samples; namely, coal mill operates in coal broken condition (coal interruption), fullofcoal condition (coal choking), and normal condition (normal operation). To facilitate the training of the SAE model, the three different operation conditions of the coal mill are labeled. The definition is shown in Table 4. The 5000 sets of data are randomly divided into training data and test data (Table 5). At the same time, two experiments are conducted to validate the effects of the proposed fault diagnosis method. The test samples of two experiments are the same.


5.2. Establishment of the SAE Model
The optimum SAE model has important effects on the accuracy rate of fault identification. In [25, 26], the unsupervised learning effect of SAE is reported to be affected by parameters of model, such as the number of nodes in the input and hidden layers, sparse parameter, and the number of times of network training. The experimental data of Experiment 1 are used as training samples, and relevant experiments are conducted to determine the optimum parameters of the SAE model. The evaluation index is the reconstruction error of the first layer of the SAE model, which is calculated by equation (4), and the experimental results are shown in Figure 8. According to the analysis in Section 3, the significantly changed variables during the fault period of the coal mill include differential pressure of primary air, outlet temperature of coal mill, and current of coal mill. Therefore, the three variables are used as the input nodes of the SAE model.
(a)
(b)
(c)
(d)
The SAE model can learn much information when the number of nodes in the input layer is large. However, the number of nodes in the input layer cannot be increased indefinitely because of the computational complexity. Figure 8(a) shows that, when the number of nodes in the input layer increases from 40 to 100, the reconstruction error of the network decreases continuously. If the number of nodes in the input layer increases further, then the reconstruction error will remain unchanged.
The number of nodes in the hidden layer determines the degree to which the model compresses the input data. The degree of compression is high when the number of nodes in the hidden layer is small. An experiment is employed using the first layer of SAE. In the experiment, the input size is set to 120 on the basis of the experiment above to determine the appropriate hidden layer parameters by analyzing the influence on the reconstruction results. As shown in Figure 8(b), when the number of hidden layer nodes is less than the input layer nodes, the reconstruction error fluctuates in a small range. This result indicates that the original data can obtain better compression when the number of hidden layer nodes is small, and this situation is conducive for the model to learn data characteristics. However, when the number of nodes in the hidden layer exceeds the number of nodes in the input layer, the reconstruction error increases rapidly, and the training effect is poor. The reason is that the sparse parameter is set to 0 at this time, and the activities of neurons in the hidden layer cannot be limited, which then leads to the poor compression effect of the SAE model on original data. Combining the constraints of complexity of network structure and computational efficiency, the number of hidden layers is set to three, and the number of nodes in each layer is 100, 50, and 25.
Sparse constraint is introduced to improve the capability of the SAE model to compress input data. Figure 8(c) shows that, when the value of is between 0.05 and 0.15, the reconstruction error of the network continues to decrease, showing that the inhibitory effect on neurons is appropriate. With the increase in the value of , the inhibitory effect on neurons is excessive, and the reconstruction error increases rapidly.
Random noise is introduced to the input data to prevent the SAE model from learning only the equivalent representation of the original data. As shown in Figure 8(d), when the probability of introducing noise is in the range of 0 to 0.1, reconstruction error decreases with the increase in noise. However, with the increase in noise, the reconstruction error increases rapidly because excessive noise causes nondestructive damage to the raw data.
In combination with the analysis above, the key parameters of the SAE model are shown in Table 6. The allocation of the nodes in the input layer is shown in Table 7. In the table, 1–40 nodes are sampled values of the differential pressure of primary air in four seconds (the sampling time is set to 0.1 s), 41–80 nodes are sampled values of outlet temperature of coal mill, and 81–120 nodes are sampled values of current of coal mill.


5.3. Validation of the Proposed Method
After determining the parameters of SAE, the network is trained by the training data obtained from simulation experiments, and then the test data are sent into the network to test the result of fault identification of the SAE model, as shown in Figure 9. The simulation results of Experiment 1 show that all samples of false diagnosis in the 194 groups of test samples are from the normal operation data. Two sets of normal operation samples are mistakenly diagnosed as coal interruption, and seven sets of normal samples are mistakenly diagnosed as coal choking. The accuracy rate of fault identification is 95.4%. To analyze the experimental results, the data curves corresponding to the misdiagnosed test samples are plotted and shown in Figure 10. The two sets of misdiagnosed samples contained in the rectangular box in Figure 9(a) correspond to data contained in the rectangular box in Figure 10. The change trend of data during this period indicates that differential pressure of primary air decreases, outlet temperature of coal mill rises, and current of coal mill decreases. These trends are consistent with the characteristics of coal mill when the coal interruption fault occurs, and these characteristics have been described in Section 3.1. Therefore, when normal operation samples are insufficient, SAE cannot fully study the differences between the two types of data and wrongly diagnoses the normal operation data as coal interruption fault.
(a)
(b)
Similarly, from the data contained in the ellipse box in Figure 10, the change trend of the data in this period is found to be consistent with the characteristics of the coal choking, differential pressure of primary air rises, outlet temperature of coal mill decreases, and current of coal mill decreases. Thus, SAE can mistakenly diagnose the normal operation data as coal choking fault. To improve the accuracy of fault diagnosis of SAE, the training samples of normal operation are increased to 1500 groups, and an experiment (Experiment 2) is conducted again. Figure 9(b) shows that, although two sets of misdiagnosed samples are still present, the accuracy of fault diagnosis of SAE has been improved to 98.7%. Therefore, if the training samples continue to increase, then the accuracy of fault diagnosis of SAE will theoretically be close to 100%.
To illustrate the effectiveness of the method that improves the performance of SAE, the accuracy rate of fault diagnosis of SAE before and after the algorithm improvement is compared. The comparison results are shown in Table 8. Training and test samples are from the data in Experiment 1. Table 8 shows that, when no improvement is implemented in SAE, the fault recognition rate is 84.02%. When sparse constraint is introduced in SAE, the fault recognition rate is increased to 85.05%. After introducing random noise into input data, the fault identification rate increases to 89.18%. When two improved methods are introduced into SAE, the network fault recognition rate further increases to 95.4%. The above analysis shows that the two improved methods proposed in this study can improve the fault recognition capability of SAE.

The coal mill is characterized by a large delay system. Detecting changes in outlet pulverized coal flow of coal mill to find the operation fault in coal mill often cannot establish early warning. Through realtime monitoring of differential pressure of primary air, outlet temperature of coal mill, and current of coal mill, and the three kinds of fast changing signals, the trained SAE can find the operation fault in the coal mill in advance. Coal interruption is taken as an example in this study. As shown in Figure 11, coal interruption fault occurs by artificial simulation. As a result, the outlet pulverized coal flow of coal mill reduces to 0 in 110 s, while the output of the SAE jumps from the normal operation state to coal interruption fault in 75 s. The network has advanced 35 s to predict the fault in the coal mill. With the adjustment of PID, the outlet pulverized coal flow of the coal mill rises gradually and goes back to the safety limit in 160 s. At this point, the output of the SAE returns to normal. Thus, the proposed method based on the deep learning algorithm can play an important role in the fault diagnosis of the coal mill.
(a)
(b)
6. Conclusions
In this study, a deep learning algorithm based on a datadriven model is proposed for fault diagnosis of coal mills. On the basis of the mechanism model of coal mills, the fault operation of coal mills is simulated and numerous fault data are obtained. Thus, the difficulty in obtaining the fault data using traditional methods is addressed. The performance of SAE is improved by introducing sparse constraints and random noise in the input layer. At the same time, the accuracy of fault diagnosis of coal mills is effectively improved, thereby enabling the possible prediction of fault in coal mills. The method proposed in the paper greatly improves the accuracy of the fault diagnosis of coal mills, which is of great significance for ensuring the safe operation of power plants. In addition, the proposed method is easy to generalize. Complex mechanical equipment in other industrial fields can use this method for fault diagnosis. The method can reduce the use of sensors for fault diagnosis of large equipment and the investment of human resources, which is essential to improve the economy and safety of the industry.
It should be noted that the paper does not consider online training for SAE. The main reason is that the online training of deep neural networks will take a lot of time, which puts stricter requirements on the performance of computers and optimization algorithms. So, we adopt a simplified method which chooses to directly use the offline training model to achieve fault diagnosis. It greatly saves the cost of calculation. However, using the online training model can continuously optimize the accuracy of the model [27–31] as the data accumulate. Therefore, how to realize the online training of the model will be the focus of our followup research.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
References
 V. Agrawal, B. K. Panigrahi, and P. M. V. Subbarao, “Review of control and fault diagnosis methods applied to coal mills,” Journal of Process Control, vol. 32, pp. 138–153, 2015. View at: Publisher Site  Google Scholar
 P. F. Odgaard and B. Mataji, “Observerbased fault detection and moisture estimating in coal mills,” Control Engineering Practice, vol. 16, no. 8, pp. 909–921, 2008. View at: Publisher Site  Google Scholar
 P. Andersen, J. D. Bendtsen, T. S. Pedersen, and B. Mataji, “Coal moisture estimation in power plant mills,” in Proceedings of the 2009 17th Mediterranean Conference on Control and Automation, pp. 1066–1071, Thessalonik, Greece, 2009. View at: Google Scholar
 J.L. Wei, J. Wang, and Q. H. Wu, “Development of a multisegment coal mill model using an evolutionary computation technique,” IEEE Transactions on Energy Conversion, vol. 22, no. 3, pp. 718–727, 2007. View at: Publisher Site  Google Scholar
 S. Guo, J. Wang, J. Wei, and P. Zachariades, “A new modelbased approach for power plant tubeball mill condition monitoring and fault detection,” Energy Conversion and Management, vol. 80, pp. 10–19, 2014. View at: Publisher Site  Google Scholar
 Y. Wang, H. R. Karimi, H. Shen, Z. Fang, and M. Liu, “Fuzzymodelbased sliding mode control of nonlinear descriptor systems,” IEEE Transactions on Cybernetics, vol. 49, no. 9, pp. 3409–3419, 2019. View at: Publisher Site  Google Scholar
 Y. Wang, X. Yang, and H. Yan, “Reliable fuzzy tracking control of nearspace hypersonic vehicle using aperiodic measurement information,” IEEE Transactions on Industrial Electronics, vol. 66, no. 12, pp. 9439–9447, 2019. View at: Publisher Site  Google Scholar
 E. A. Garcia, Z. Han, and P. M. Frank, “FDI based on parameter and output estimation: an integrated approach,” in Proceedings of the European Control Conference, Karlsruhe, Germany, 1999. View at: Google Scholar
 T. Hofling and R. Isermann, “Adaptive parity equations and advanced parameter estimation for fault detection and diagnosis,” in Proceedings of the IFAC World Congress, pp. 55–60, San Francisco, CA, USA, July 1996. View at: Google Scholar
 P. Nomikos and J. F. MacGregor, “Monitoring batch processes using multiway principal component analysis,” AICHE Journal, vol. 40, no. 8, pp. 1361–1375, 1994. View at: Publisher Site  Google Scholar
 Z.G. Su, P.H. Wang, X.J. Yu, and Z.Z. Lv, “Experimental investigation of vibration signal of an industrial tubular ball mill: monitoring and diagnosing,” Minerals Engineering, vol. 21, no. 10, pp. 699–710, 2008. View at: Publisher Site  Google Scholar
 E. Kisić, V. Petrović, S. Vujnović, Ž. Đurović, and M. Ivezić, “Analysis of the condition of coal grinding mills in thermal power plants based on the T^{2} multivariate control chart applied on acoustic measurements,” Facta UniversitatisSeries: Automatic Control and Robotics, vol. 11, no. 2, pp. 141–151, 2012. View at: Google Scholar
 S. Collura, D. Possanzini, D. Pestonesi, and M. Gualerci, Coal Mill Performances Optimization through Noninvasive Online Coal Fineness Monitoring, Powergen, Vienna, Austria, 2013.
 X. Han and X. Jiang, “Fault diagnosis of pulverizing system based on fuzzy decisionmaking fusion method,” in Fuzzy Information and Engineering, vol. 2, pp. 1045–1056, Springer Berlin Heidelberg, Berlin, Germany, 2009. View at: Google Scholar
 W. Qin, W. Yan, and J. Xu, “Application of fault diagnosis expert system in grinding process,” in Proceedings of the Automation and Logistics (ICAL), pp. 290–295, Hong Kong, China, 2010. View at: Google Scholar
 Y. Wang, H. R. Karimi, H.K. Lam, and H. Shen, “An improved result on exponential stabilization of sampleddata fuzzy systems,” IEEE Transactions on Fuzzy Systems, vol. 26, no. 6, pp. 3875–3883, 2018. View at: Publisher Site  Google Scholar
 R. Isermann and P. Balle, “Trends in the application of model based fault detection and diagnosis of technical processes,” in Proceedings of the IFAC World Congress, pp. 55–60, San Francisco, CA, USA, July 1996. View at: Google Scholar
 P. M. Frank, “Fault diagnosis in dynamic systems using analytical and knowledgebased redundancy,” Automatica, vol. 26, no. 3, pp. 459–474, 1990. View at: Publisher Site  Google Scholar
 X. Guo, L. Chen, and C. Shen, “Hierarchical adaptive deep convolution neural network and its application to bearing fault diagnosis,” Measurement, vol. 93, pp. 490–502, 2016. View at: Publisher Site  Google Scholar
 Y. Duan, Y. Lv, Y. L. Liu et al., “An efficient realization of deep learning for traffic data imputation,” Transportation Research Part C: Emerging Technologies, vol. 72, pp. 168–181, 2016. View at: Google Scholar
 Y. Gao, D. Zeng, J. Liu, and Y. Jian, “Optimization control of a pulverizing system on the basis of the estimation of the outlet coal powder flow of a coal mill,” Control Engineering Practice, vol. 63, pp. 69–80, 2017. View at: Publisher Site  Google Scholar
 D.L. Zeng, Y. Hu, S. Gao, and J.Z. Liu, “Modelling and control of pulverizing system considering coal moisture,” Energy, vol. 80, pp. 55–63, 2015. View at: Publisher Site  Google Scholar
 D. Simon, “Kalman filtering with state constraints: a survey of linear and nonlinear algorithms,” IET Control Theory & Applications, vol. 4, no. 8, pp. 1303–1318, 2010. View at: Publisher Site  Google Scholar
 V. Rashtchi, E. Rahimpour, and S. Fazli, “Genetic algorithm application to detect broken rotor bar in three phase squirrel cage induction motors,” International Review of Electrical Engineering, vol. 6, no. 5, pp. 2286–2292, 2011. View at: Google Scholar
 H. Pan and X. W. Wang, “Survey on collaborative filtering recommendation algorithm based on extreme learning machine stacked denoising autoencodes,” Application Research of Computers, vol. 33, no. 8, 2016. View at: Google Scholar
 Y. Jian, X. Qing, H. Liang, Y. Zhao, X. Qi, and M. Du, “Fault diagnosis of motor bearing based on deep learning,” Advances in Mechanical Engineering, vol. 11, no. 9, pp. 1–9, 2019. View at: Publisher Site  Google Scholar
 A. Das and F. L. Lewis, “Distributed adaptive control for synchronization of unknown nonlinear networked systems,” Automatica, vol. 46, no. 12, pp. 2014–2021, 2010. View at: Publisher Site  Google Scholar
 S. ElFerik, H. A. Hashim, and F. L. Lewis, “Neuroadaptive distributed control with prescribed performance for the synchronization of unknown nonlinear networked systems,” IEEE Transactions on Systems,Man, and Cybernetics: Systems, vol. 48, no. 12, pp. 2135–2144, 2017. View at: Publisher Site  Google Scholar
 S. ElFerik, A. Qureshi, and F. L. Lewis, “Neuroadaptive cooperative tracking control of unknown higherorder affine nonlinear systems,” Automatica, vol. 50, no. 3, pp. 798–808, 2014. View at: Publisher Site  Google Scholar
 H. A. Hashim, S. ElFerik, and F. L. Lewis, “Adaptive synchronisation of unknown nonlinear networked systems with prescribed performance,” International Journal of Systems Science, vol. 48, no. 4, pp. 885–898, 2017. View at: Publisher Site  Google Scholar
 F. L. Lewis, H. Zhang, K. HengsterMovric, and A. Das, Cooperative Control of MultiAgent Systems: Optimal and Adaptive Design Approaches, Springer Science & Business Media, Berlin, Germany, 2013.
Copyright
Copyright © 2020 Yifan Jian et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.