Next Article in Journal / Special Issue
Classifying Image Stacks of Specular Silicon Wafer Back Surface Regions: Performance Comparison of CNNs and SVMs
Previous Article in Journal
Symmetry of Gait in Underweight, Normal and Overweight Children and Adolescents
Previous Article in Special Issue
Research on Modeling and Analysis of Generative Conversational System Based on Optimal Joint Structural and Linguistic Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Hybrid Multi-Objective Optimization Model for Vibration Tendency Prediction of Hydropower Generators

1
MOE Key Laboratory of Image Processing and Intelligence Control, School of Artificial Intelligence and Automation, Huazhong University of Science and Technology, Wuhan 430074, China
2
School of Hydropower and Information Engineering, Huazhong University of Science and Technology, Wuhan 430074, China
3
School of Mechanical Engineering and Electronic Information, China University of Geosciences, Wuhan 430074, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(9), 2055; https://doi.org/10.3390/s19092055
Submission received: 9 April 2019 / Revised: 28 April 2019 / Accepted: 29 April 2019 / Published: 2 May 2019

Abstract

:
The hydropower generator unit (HGU) is a vital piece of equipment for frequency and peaking modulation in the power grid. Its vibration signal contains a wealth of information and status characteristics. Therefore, it is important to predict the vibration tendency of HGUs using collected real-time data, and achieve predictive maintenance as well. In previous studies, most prediction methods have only focused on enhancing the stability or accuracy. However, it is insufficient to consider only one criterion (stability or accuracy) in vibration tendency prediction. In this paper, an intelligence vibration tendency prediction method is proposed to simultaneously achieve strong stability and high accuracy, where vibration signal preprocessing, feature selection and prediction methods are integrated in a multi-objective optimization framework. Firstly, raw sensor signals are decomposed into several modes by empirical wavelet transform (EWT). Subsequently, the refactored modes can be obtained by the sample entropy-based reconstruction strategy. Then, important input features are selected using the Gram-Schmidt orthogonal (GSO) process. Later, the refactored modes are predicted through kernel extreme learning machine (KELM). Finally, the parameters of GSO and KELM are synchronously optimized by the multi-objective salp swarm algorithm. A case study and analysis of the mixed-flow HGU data in China was conducted, and the results show that the proposed model performs better in terms of predicting stability and accuracy.

1. Introduction

The role of hydropower stations in the frequency control and peak regulation of the power grid has become increasingly important. Thus, the security and stable operation of the hydropower generator unit (HGU) is crucial to the stability of the power system [1,2]. In practical engineering, the vibration signal is a typical indicator for evaluating the health status and stability of the HGU [3]. Therefore, vibration tendency prediction of the HGU is essential in ensuring safe operation and conducting condition maintenance. It contributes to improving the comprehensive economic benefits of the station. Thus, to guarantee the security of the HGU and avoid economic losses at the station, it is of great importance to study the condition monitoring and tendency prediction of HGUs [4].
Essentially, the vibration tendency prediction of HGU is time series prediction according to historical values. The key aspects for vibration tendency prediction are signal processing and fitting regression [4]. The essential information that reflects the operating states under the environment of noise and electromagnetic interference can be obtained by signal processing [5,6,7]. Then, the vibration tendency of an HGU can be derived by generalized regression. Recently, many advanced methods for tendency forecasting have been applied in engineering, as summarized in Table 1 [4,8,9,10,11,12,13,14,15,16,17]. For instance, Pham et al. proposed a method to forecast the condition of the machine using the auto-regressive moving average (ARMA) [10]. Milovancevic et al. utilized an artificial neural network (ANN) to predict the gear transmission vibration of pellet mills [11]. Fu et al. combined aggregated ensemble empirical mode decomposition (EEMD) with support vector regression (SVR) to predict the vibration tendency of HGUs [12]. Javed et al. applied extreme learning machine (ELM) to forecast the vibration of a turbine engine [13]. In the above methods, ARMA is mainly applied to solve linear problems. However, it is not sufficient to deal with the issue of random fluctuations. Despite the fast calculation of ANN, its performance can be affected by the initial parameters. Moreover, SVR based on structural risk minimization tends to be time-consuming when dealing with constrained optimization issues [18]. Although ELM is widely used for prediction due to its faster convergence speed, considering only the empirical risk minimization principle, it cannot easily obtain the optimal model. In this paper, the kernel extreme learning machine (KELM) is employed as the prediction model, which subtly applies the generalized inverse and least squares to save the computing time [19].
To obtain better prediction results, time-series preprocessing and feature selection methods are also necessary. In term of data processing, short-time Fourier transform (STFT [20]), wavelet transform (WT [21]), empirical mode decomposition (EMD [22]) and its modified methods have been utilized to analyze the vibration data series. In recent years, Gilles proposed a new signal processing method, i.e., empirical wavelet transform (EWT) [23]. It combines the advantages of WT and EMD, which not only has a mathematical theory like WT, but can also decompose signals adaptively like EMD [24]. Hu et al. provided a hybrid model combining EWT with least square support vector machine for wind speed forecasting [25]. The results verified the predictive ability of the hybrid model. Nevertheless, the existing EWT-based forecasting models mainly accumulate all the modes of predicted values. However, there may be some false components among the modes, which will affect the model precision. In addition, it is more time-consuming to predict all the modes. The feature selection is helpful to select the optimal inputs for the prediction model, which improves model accuracy [18,26]. Some methods have been successfully applied in engineering such as partial autocorrelation function (PACF) [27], principal component analysis (PCA) [28] and the Gram–Schmidt orthogonal (GSO) process [26,29].
It has been proven that model optimization can enhance the prediction performance. In general, the objective of model optimization is to obtain the optimal parameters by minimizing the prediction errors [30,31]. For instance, ELM models are optimized by the grey wolf optimizer (GWO) [24] and gravitational search algorithm (GSA) [26]. However, the above models ignore the importance of stability and accuracy, which are equally necessary for evaluating the performance of the prediction model. Unfortunately, these two criteria are almost independent, and can be expressed as multi-objective problems (MOPs) [32]. Single objective optimization algorithms like GWO and GSA will not work. Thus, to handle MOPs, it is necessary to realize the two objectives of good stability and high accuracy simultaneously. In recent decades, driven by engineering practices like water resources scheduling, wind speed forecasting and power load prediction, MOPs have aroused great interest for researchers [33,34]. Pareto optimal solutions can be obtained by multi-objective optimization techniques. Then, they can be presented to the decision maker as a set of feasible optimal options in practice. The multi-objective salp swarm algorithm (MOSSA) is an effective and practical method, which has been widely used in solving MOPs [32,35]. Thus, the MOSSA is introduced to simultaneously realize good stability and high accuracy in this paper.
According to the above discussions, a hybrid model based on AEWT, GSO, KELM and MOSSA is proposed for forecasting vibration tendency. The main contributions are as follows: (1) To eliminate illusive components and reduce the computing time, a mode reconstruction strategy is presented with EWT and sample entropy (SE); (2) MOSSA is applied to find the optimal parameters of GSO and KELM, which can simultaneously obtain good stability and high accuracy of the forecast model; and (3) The hybrid model of AEWT-GSO-KELM in a MOSSA-based framework is introduced for vibration tendency forecasting of HGUs.
The rest of this paper is presented as follows. The framework of the proposed model is illustrated in detail in Section 2, including data preprocess, feature selection, multi-objective optimization, the forecasting and evaluation. Then, Section 3 describes the application results and analysis of the experimental cases. Finally, the conclusions are given in Section 4.

2. Proposed Predicting Model and Framework

To simultaneously obtain higher accuracy and more stable prediction results, a hybrid model is proposed for HGU vibration trend forecasting. Figure 1 shows the basic flow chart. The general procedure contains signal preprocess, features selection, multi-objective optimization, the best compromise solution selection, forecasting and evaluation. The specific procedures are described as follows:
Step 1: Use AEWT to decompose the vibration series into several modes and refactor all modes with the SE-based strategy.
Step 2: Apply GSO method to select the appreciated input features for KELMs.
Step 3: Optimize the parameters N of GSO, C and σ of KELM by MOSSA according to the two objective functions.
Step 4: Use the fuzzy evaluation method to select the compromise solution, i.e., the optimized parameters from the available Pareto front.
Step 5: Establish an optimized KELM prediction model with the optimal parameters of the model and calculate all the forecasting values.
Step 6: Evaluate the final results through several assessment indices quantitatively.

2.1. Data Preprocess

The vibration signal is influenced by electromagnetic and environmental factors, which may affect the forecasting result. Thus, it should be preprocessed before the establishment of the forecasting model.

2.1.1. Empirical Wavelet Transform (EWT)

EWT is an adaptive signal processing method [23]. In the EWT framework, the first main step is to divide the frequency domain of the signal according to the Fourier spectrum [36]. Another step is contracting the empirical wavelets, which can extract the AM-FM components with compact Fourier spectrum support [37].
The empirical wavelet function ψ ^ n ( ω ) the empirical scale function ϕ ^ n ( ω ) can be respectively determined according to Meyer’s wavelet [38] as follows:
ψ ^ n = { 1 i f ( 1 + γ ) ω n + 1 | ω | ( 1 γ ) ω n cos [ π 2 β ( 1 2 γ ω n + 1 ( | ω | ( 1 γ ) ω n + 1 ) ) ] i f ( 1 γ ) ω n + 1 | ω | ( 1 + γ ) ω n + 1 sin [ π 2 β ( 1 2 γ ω n ( | ω | ( 1 γ ) ω n ) ) ] i f ( 1 γ ) ω n + 1 | ω | ( 1 + γ ) ω n + 1 0 o t h e r w i s e ,
ϕ ^ n ( ω ) = { 1 i f | ω | ( 1 γ ) ω n cos [ π 2 β ( 1 2 γ ω n ( | ω | ( 1 γ ) ω n ) ) ] i f ( 1 γ ) ω n | ω | ( 1 + γ ) ω n 0 o t h e r w i s e ,
where γ min n ( ω n + 1 ω n ω n + 1 + ω n ) is a coefficient that satisfies the tight frame.
Based on the above, the EWT ( W f ε ( n , t ) ) can be constructed as classic WT. The approximation coefficients W f ε ( 0 , t ) and W f ε ( n , t ) are illustrated as follows, respectively.
W f ε ( 0 , t ) = f , ϕ 1 W f ε ( n , t ) = f , ψ n
Finally, the empirical modes fk are given by:
f 0 ( t ) = ϕ 1 ( t ) W x ( 0 , t ) f k ( t ) = ψ k ( t ) W x ( k , t )

2.1.2. AEWT Based on Sample Entropy Theory

The decomposition results of EWT include many modes, some of which are illusive components, and this consumes a lot of time and reduces the forecasting accuracy. To solve the above problems, a mode reconstruction strategy based on sample entropy (SE) is proposed.
According to Richman et al.’s work, the SE is utilized to indicate the randomness of the series [39]. And the SE of obtained modes { x ( i )   i = 1 , 2 , , j } can be calculated as:
P S E ( i ) = S a m p E n ( x ( i ) , m , r , N ) ,
where PSE represents the sample entropy, and m, r are the reconstruction dimension and threshold for calculating the SE, respectively.
With the idea of approximate aggregation, the reconstructed mode (RM) can be obtained with the SE values of the modes series PSE1, PSE2, …, PSEn.
R M = i i + j M o d e i ,
s . t . | P S E ( n ) P S E ( n i ) | P S E ( max ) P S E ( min ) n / 2 .
where n represents the number of modes.
Through the reconstruction strategy, the modes will be adaptively divided into groups, where the PE values are closer. The obtained RMs can reduce the illusive components and time consumption. Thus, RMs can be used as input for the proposed model to predict vibration tendency of HGU.

2.2. Feature Selection

GSO is a simple and effective method for feature selection, which sorts the correlation between input feature and output target in descending order by a criterion [40]. Suppose the input vector with k features of M patterns in a system is X k = [ x k 1 , x k 2 , , x k M ] T and Y = [ y 1 , y 2 , , y M ] T is the corresponding output target. To select the features most relevant to output, the angle between each input X k and output target Y can be applied as the criterion:
cos ( φ k ) = X k · Y X k Y , k = 1 , 2 , , N
where φ k denotes the angle between the k-th input X k and the output target Y, N represents the number of features. If the input is completely unrelated to output, the φ k is π / 2 . Conversely, φ k is 0.
In the iteration process, the features containing the maximum value of the mentioned criterion are chosen as the most relevant feature. For the next iteration, the other candidates and the output vector will be projected into null space of the selected feature. The procedure does not end until all candidate are ranked [41]. Finally, the first Nf features are chosen as the target features. Meanwhile, the number of selected features Nf affects the prediction accuracy and stability.

2.3. Kernel Extreme Learning Machine (KELM)

By introducing the kernel function into ELM, KELM gets the least-squares optimal solution, which performs better in generalization [42]. Meanwhile, compared to the traditional single output SVR, KELM can achieve multi-output, which reduces the training time.
The dataset ( x i , t i ) , x i = [ x i 1 , x i 2 x i n ] T R n and y i = [ y i 1 , y i 2 y i m ] T R m represent the input vector and the corresponding output vector, respectively. With the activation function including L hidden layer nodes, the output of ELM is calculated as follows [43]:
f L ( x i ) = j = 1 L β j h j ( x i ) = h ( x ) β = y i j = 1 , 2 , , N ,
where β is the output weight vector, and h ( x i ) represents the output of the hidden nodes concerning the inputs x i .
The kernel matrix of ELM is obtained by applying the Mercer condition [43]:
Ω E L M = H H T Ω E L M i , j = h ( x i ) · h ( x j ) = K ( X i , X j )
where H is the output matrix of the hidden layer:
Then, based on the ridge regression theory, the output function of KELM is formulated:
f ( x i ) = h ( x i ) H T ( I C + H H T ) 1 Y = [ K ( x , x 1 ) K ( x , x 2 ) K ( x , x N ) ] ( I C + Ω E L M ) ) 1 Y
The kernel width σ and penalty parameter C are both essential for predictive performance. Hence, the best combinations of such parameters should be optimized by advanced intelligent algorithms.

2.4. Multi-Objective Optimization

Multi-objective optimization involves multiple objectives, all of which should be optimized simultaneously [44]. In single objective optimization, it is easy to select the best one from a number of solutions through relational operators, while in MOPs, there is more than one criterion, which makes it impossible to compare the solutions. As described in [45,46], Pareto optimal dominance is the main operator that can compare different solutions in MOPs.

2.4.1. Multi-Objective Salp Swarm Algorithm (MOSSA)

MOSSA is an evolutionary computation method based on swarm intelligence. It is inspired by precisely positioning and searching of salps in complex environments [35,47]. With mathematical simulation of the salp chains, they are composed of leaders and followers. Assuming that the food source F is the target of swarm in the search space [32], the leader can update its position according to the food source that is defined as follows:
x j 1 = { F j + c 1 ( ( u b j l b j ) c 2 + l b j ) c 3 0 F j c 1 ( ( u b j l b j ) c 2 + l b j ) c 3 0 ,
where x j 1 and F j represent the position of leader salp and the food source in j-th dimension, respectively.
The coefficient c1 can balance exploitation and exploration:
c 1 = 2 e ( ( 4 l / L ) 2 ) ,
where L represents the maximum iteration, and l is the current iteration.
The coefficient parameters c2 and c3 are random numbers in [0,1].
Then the position of the followers is updated by utilizing Newton’s law of motion:
x j i = v f i n a l 2 v 0 t 2 + v 0 t i > 1 ,
v = ( x x 0 ) t ,
where x j i is the position of follower salp. v 0 and t are the initial speed and time, respectively.
Since the time is iteration in optimization, considering v 0 = 0 and the discrepancy between iterations is equal to 1, Equation (14) is written as follows.
x j i = ( x j i + x j i 1 ) / 2 i > 1
To explain the MOSSA more clearly and intuitively, the pseudo code is detailed in Algorithm 1.
Algorithm 1: MOSSA
  Parameters:
  iter_max—The maximum number of iterationslb—Variables lower bound
  obj_no—Number of objective functionsub—Variables upper bound
  dim—Number of decision variablesN—The salp population
  archive_maxsize—The maximum number of archive
  /* Set the parameters of MOSSA*/
  /* Initialize the salp population x i ( i = 1 , 2 , , n ) considering ub and lb*/
  while (iter < iter_max) do
  /* Calculate the objective values for each salp *//* Determine the non-dominated salps */
  /* Update repository with the obtained non-dominated salps */
  if (the repository becomes full) then
  /* Remove the repository resident by calling the repository maintenance procedure*/
  /* Add the non-dominated salp to the repository */
  end if
  /* choose the food source: F = Select Food (repository) */
  update c1
  for each salp (i = 1 to N) do
  if (i = = 1) then
  update the position of the leading salp
  else if ( i 1 ) then
  update the position of the follower salp
  end if
  end for
  /* Adjuse the salps considering the ub and lb */
  end while
  /* Return repository */

2.4.2. Optimization Strategy

The MOSSA is firstly employed to optimize the parameters of GSO and KELMs simultaneously. As illustrated in Figure 2, the parameters are encoded in the location vector in MOSSA. The first part represents the number of selected features. The remained part is the parameters C and σ , which have great influence on the performance of KELM.

2.4.3. The Fitness Function

In this paper, the effectiveness of the prediction model should be simultaneously evaluated by the prediction accuracy and stability. Therefore, the bias-variance indexes are utilized as the fitness functions of MOSSA [30,31]. The value of B i a s ( y ^ ) represents the precision, which indicates the difference between the true value and the predicted value.
B i a s ( y ^ ) = y E ( y ^ )
where y ^ is the predicted value the n-th datum, y and E ( y ^ ) denote the expectation of the observed value and forecasting value, respectively.
Similar to B i a s ( y ^ ) , the V a r ( y ^ ) indicates the stability of the model. It can represent the variability of the forecast results [48]. The V a r ( y ^ ) is defined as follows:
V a r ( y ^ ) = E ( y ^ E ( y ^ ) ) 2
Consequently, the fitness function focusing on high accuracy and stability is formulated as follows:
min { o b f 1 ( x ) = B i a s ( y ^ ) o b f 2 ( x ) = V a r ( y ^ )

2.5. Evaluation Criterion

To evaluate the performance of different models quantitatively, four evaluated indexes, such as root mean square error (RMSE), mean absolute error (MAE), mean absolute percentage error (MAPE) and coefficient of correlation (R), are employed [18,49].
R M S E = i = 1 N ( y i y ^ i ) 2 N
M A E = 1 N I = 1 N | y i y ^ i |
M A P E = 100 % × 1 N I = 1 N | y i y ^ i y i |
R = ( 1 / N ) i = 1 N ( y i y ¯ ) ( y ^ i y ^ ¯ ) ( 1 / N ) i = 1 N ( y i y ¯ ) 2 × ( 1 / N ) i = 1 N ( y ^ i y ^ ¯ ) 2
where N represents the sample size. y ^ i and y i are the predicted value and actual value, respectively.

3. Engineering Application and Analysis

To verify the effectiveness of the proposed model, the experiment was carried out in the Ertan Hydropower Station in China. The S8000 monitoring system was installed [4], which is made by Shenzhen Strongwish Co., Ltd. (Shenzhen, China). It contains a data acquisition and monitoring substation (NET8000) and central server (WEB8000). The NET8000 acquires the raw signal with sensors and conducts network communication. It can collect up to 4 key phase signals, 24 vibration signals, 12 static signals and 255 process signals obtained by MODBUS communication. The WEB8000 can store and manage the data from NET8000. The general structures of HGU is illustrated in Figure 3, where the location of sensors is shown in on the right. The BENTLY3300 sensors were mounted on the lower guide to obtain the swing data. The sensor is a kind of electric eddy current sensor of the 3300 series, whose output is 4–20 mA or 1–5 V. Its frequency response range is 0–10 kHz and the linear range is 2 mm, average sensitivity is 200 mV/mil [50].

3.1. Data Collection

The complicated structure and frequent switching of operating conditions may easily lead to different monitoring data with different time intervals. Hence, the monitoring data that meets the average time interval was chosen for analysis, which is in line with engineering practice [12]. Thus, the swing data in the X direction from 24 July 2011 to 27 July 2011 were obtained for the experiment. The detailed description of the monitoring data is presented in Table 2. As shown in Figure 4, there are 300 samples with a time interval of 10 min and each sample contains one point. It is obvious that the series has strong non-stationary characteristics. Among all experiments, 3/4 of the entire data were selected randomly as the training set, and the remaining 1/4 were the testing set.

3.2. Model Description

To prove the superiority of the proposed model, four types of models were applied and analyzed. The first type includes EEMD-KELM, EWT-KELM and AEWT-KELM, which consist of signal decomposing and model prediction. The second includes AEWT-PCA-KELM and AEWT-GSO-KELM, where PCA and GSO are used for feature selection. The third includes AEWT-GSO-SVR and AEWT-GSO-KELM, where SVR and KELM are utilized to predict the vibration trend. The final type includes AEWT-GSO-SSA-KELM, AEWT-GSO-MOPSO-KELM and AEWT-GSO-MOSSA-KELM.
The above experiments were carried out using the raw data shown in Figure 4. Among the above models, the parameters of the first three models were optimized by MOSSA, while the parameters in the fourth model were optimized by the MOPSO and MOSSA, respectively. The parameters of MOPSO and MOSSA are set as follows: the number of individual agents is 80, the maximum iteration is 100. In addition, for the GSO, the number of selected features is determined in the range [3,12] and the step size is 1.

3.3. Vibration Tendency Prediction of the HGU

3.3.1. Vibration Signal Decomposition and Mode Reconstruction

As shown in Figure 4, the vibration signal fluctuates severely with no clear regularity. Hence, EWT was used to decompose the signal series for decreasing the noise components. The results are shown in Figure 5, where the series is decomposed into several modes.
Due to the great impact of IMF complexity on prediction accuracy and computation efficiency of the model, a SE-based IMF reconstruction method was proposed for signal processing. The detailed SE values of the modes are presented in Table 3. It can be seen that the SE values are gradually increasing. This indicates that the complexities of the obtained modes are becoming more and more complex. Based on the proposed reconstruction method, the interval of SE H s max H S min n / 2 is set as 0.056. The constructed RMs after mode reconstruction are presented in Table 4. According to the proximity of the SE value, Mode 4, Mode 5 and Mode 6 were selected to construct the RM3. The RM1, RM2 and RM4 were generated by a similar principle. All the obtained RMs are presented in Figure 6, which show significantly different characteristics. More specifically, RM4 shows the strongest nonlinear and the highest frequency, while RM1 is the most stable along the time going. Since the RMs have obvious characteristics and relatively few components, they can be used as inputs of the prediction model to further improve the efficiency and accuracy.

3.3.2. The Selection of the Best Compromise Solution

To avoid randomness and contingency in the optimization, all the experiments were independently repeated 10 times under the same conditions. Therefore, the non-dominated solutions constitute the Pareto front as depicted in Figure 7. It can be seen that MOSSA can easily generate feasible solutions that cover the entire Pareto front. These solutions are distributed with two contradictory objectives. Thus, the fuzzy evaluation method was used to efficiently search for the best compromise solution from the Pareto-optimal front [32,51].
The fitness function of each solution can be obtained with the fuzzy membership function ui:
u i = f i , max f i f i , max f i , min ,
where f i , max and f i , min are the maximum and minimum fitness values, respectively. f i is the fitness value of the i-th non-dominated solution.
The standardized fuzzy evaluation value u can be defined as:
u = i = 1 m u i / m
where m represents the number of non-dominated solutions.
The maximum value u can be chosen as the best compromise solution. In addition, it should be emphasized that the fuzzy evaluation method is only utilized to select the compromise solution, but does not focus on the improvement of the method.

3.3.3. Analysis and Discussion

Based on RMs, the vibration tendency prediction results of the HGU can be obtained by using the above models. The inductors including RMSE, MAE, MAPE and the computing time are utilized to evaluate the performance of the models. Comparing the models further indicates that the integration of reconstruction, feature selection and multi-objective optimization are helpful to improve the prediction accuracy.
Driven by test data, the predicted vibration trend of the HGU, achieved by eight models is illustrated in Figure 8. The results are worthy of further analysis and discussion.
(1)
Compared with the other models of the first type, AEWT-KELM obtains better results than EEMD-KELM and EWT-KELM, which demonstrates that the SE-based reconstruction strategy can improve the forecasting precision.
(2)
Among the second type, AEWT-GSO-KELM performs better than AEWT-PCA-KELM, which indicates that GSO is more effective as a feature selection method in these experiments.
(3)
By comparing the first and second type models, i.e., AEWT-KELM vs. AEWT-GSO-KELM, we can conclude that the hybrid model integrating the GSO method enhances its performance in forecasting.
(4)
Compared with the SVR-based model, the KELM can realize a higher degree of forecasting, which suggests process that the ability of KELM is enhanced by the kernel function.
(5)
By comparing the final type models, it can be observed that the proposed model optimized by MOSSA performs better than that by MOPSO, which indicates MOSSA is superior to MOPSO in handling MOPs.
Overall, it is obvious that the prediction curve of the proposed model matches the actual curve best. In contrast, the other models do not fit well with the actual vibration trend. For the SVR-based model, some points of the output are almost divergent in Figure 8e.
The dispersion of the prediction errors is exhibited by box plots as shown in Figure 9, where the eight models are compared. The error index is the relative error between the predicted values and the real ones. The coordinates of the x-axis represent the compared model from EEMD-KELM to the proposed model as ordered in Figure 8. Compared to other models, it can be found that the errors of the proposed model are normally distributed to zero. The errors of the all eight models decrease on tiny scales. This means that the proposed model can achieve uniform precision at almost all points.
To quantify the prediction results, the evaluated indexes of the eight models are illustrated in Table 5. The first three rows represent the results of the first type models. The next four rows are the results of the second and third type models. The remained rows represent the results of the fourth type models, including the proposed model.
From Table 5, the RMSE of EWT-KELM is 1.207, which is lower than that of EEMD-KELM (1.236). The reason is that the EWT method can avoid the mode mixing of EEMD to some extent. Moreover, the RMSE, MAE, MAPE and R of AEWT-KELM are superior to those of EWT-KELM. This shows that the SE-based refactor strategy can obviously enhance the prediction performance and reduce the complexity of the modes. Next, the RMSE of AEWT-GSO-KELM is 0.823, which is lower than that of AEWT-PCA-KLEM (0.919) because GSO can capture the key feature of the series and simplifies the input of KELM. Compared with the SVR-based model, the KELM-based model has better performance, which indicates that KELM has great generalization ability and the over-fitting is reduced by adopting the kernel function. Finally, for MOSSA, the four indexes of the proposed model are 0.823, 0.65, 0.0068 and 0.913 respectively, while the single-objective optimization SSA obtains 0.939, 0.743, 0.0079 and 0.906. This indicates that the performance of the multi-objective optimization algorithm is better than that of single-objective optimization algorithms. Meanwhile, the performance of AEWT-GSO-MOSSA-KELM is superior to that of AEWT-GSO-MOPSO-KELM. For example, the RMSE and MAE of AEWT-GSO-MOSSA-KELM are 0.823 and 0.65, respectively, while those of AEWT-GSO-MOPSO-KELM are 0.841 and 0.731, respectively. This proves that the ability of MOSSA to solve MOP is better than that of MOPSO.
In addition, the correlation coefficient R of the real vibration tendency series and the predicted series by different models are exhibited in the fifth column in Table 5. The results show that the proposed model obtains the highest similarity and correlation with the real series among the different models. In summary, it can be concluded that the proposed model performs best in prediction compared with other models.
Apart from accuracy, efficiency is also a critical goal for model design. Both of these are essential for model evaluation. Thus, the computing time was used to evaluated the efficiency of different models. All the experiments were conducted on the platform with 16 GB of RAM and 3.0 GHz CPU in the MATLAB R2014a environment. The computing time of the nine models in vibration tendency prediction is illustrated in Table 4. As can be seen from the results, the proposed model achieved a better balance between accuracy and efficiency. For example, the time of EWT-GSO-KELM is 105.138 s, which is higher than 97.276 s of AEWT-GSO-KELM. Meanwhile, the RMSE of the former is 1.105, which indicates that the SE-based reconstruction method effectively reduces the computing time and improves the accuracy. In addition, the time of AEWT-GSO-SSA-KELM is 95.716 s, which is lower than that of AEWT-GSO-MOSSA-KELM and AEWT-GSO-MOPSO-KELM. However, its RMSE is 0.939, which is not a satisfactory result. In conclusion, the proposed model achieves a better balance between accuracy and efficiency.

4. Conclusions

In this paper, a new vibration tendency prediction model is proposed by integrating SE reconstruction, feature selection and multi-objective optimization. Firstly, raw monitoring data is decomposed into several modes by EWT. The RMs are obtained from all the modes according to the proposed SE reconstruction strategy. Subsequently, the GSO is used to select features for the KELM input for dimensionality reduction. Later, the AEWT-GSO-KELM model is established with the parameters of GSO and KELM optimized by MOSSA, thus simultaneously achieving high accuracy and stability. The proposed AEWT-GSO-MOSSA-KELM model is verified in different comparative experiments, which include the EEMD-based, PCA-based, SVR-based, SSA-based and MOPSO-based models. The results demonstrate that the proposed hybrid model has better comprehensive performance in prediction considering the two objectives of accuracy and stability.
The proposed model was only applied for analyzing data in the existing Ertan station in China. However, if other power stations or other parts of the HGU are to be promoted, some improvements should be made to enhance its applicability and robustness. For example, the fuzzy evaluation method can be improved to select the composited solution. The multi-prediction model will be constructed to predict the data series from different power stations. Therefore, the above-mentioned is also the focus of our subsequent research work.

Author Contributions

Conceptualization, Y.S.; Formal analysis, Y.S. and G.-N.C.; Methodology, K.-B.Z. and Y.S.; Project administration, J.-Y.Z.; Validation, Z.-Y.G.; Writing—Original draft, M.-F.G.; Writing—Review & editing, K.-B.Z.

Funding

This work is supported in part by the National High Technology Research and Development Program of China (No. 2013AA040705), the National Natural Science Foundation of China (No. 61873101) and the Fundamental Research Funds for the Central Universities (No. 2019kfyXJJS137).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zheng, Y.; Zhou, J.; Zhu, W.; Zhang, C.; Li, C.; Fu, W. Design of a multi-mode intelligent model predictive control strategy for hydroelectric generating unit. Neurocomputing 2016, 207, 287–299. [Google Scholar] [CrossRef]
  2. Xu, Y.; Zheng, Y.; Du, Y.; Yang, W.; Peng, X.; Li, C. Adaptive condition predictive-fuzzy PID optimal control of start-up process for pumped storage unit at low head area. Energy Convers. Manag. 2018, 177, 592–604. [Google Scholar] [CrossRef]
  3. Zhang, X.; Zhou, J.; Guo, J.; Zou, Q.; Huang, Z. Vibrant fault diagnosis for hydroelectric generator units with a new combination of rough sets and support vector machine. Expert Syst. Appl. 2012, 39, 2621–2628. [Google Scholar] [CrossRef]
  4. Fu, W.; Wang, K.; Li, C.; Li, X.; Li, Y.; Zhong, H. Vibration trend measurement for a hydropower generator based on optimal variational mode decomposition and an LSSVM improved with chaotic sine cosine algorithm optimization. Meas. Sci. Technol. 2018, 30, 015012. [Google Scholar] [CrossRef]
  5. Liu, J.; Hu, Y.; Wang, Y.; Wu, B.; Fan, J.; Hu, Z. An integrated multi-sensor fusion-based deep feature learning approach for rotating machinery diagnosis. Meas. Sci. Technol. 2018, 29, 055103. [Google Scholar] [CrossRef] [Green Version]
  6. Wang, Y.; Hu, Y.; Wu, B.; Liu, J. An integrated condition-monitoring method for a milling process using reduced decomposition features. Meas. Sci. Technol. 2017, 28, 085101. [Google Scholar]
  7. Hu, Z.; Wang, Y.; Ge, M.-F.; Liu, J. Data-driven Fault Diagnosis Method based on Compressed Sensing and Improved Multi-scale Network. IEEE Trans. Ind. Electron. 2019. [Google Scholar] [CrossRef]
  8. Xia, X.; Meng, Y.; Shi, B.; Qiu, M. Bootstrap Forecasting Method of Uncertainty for Rolling Bearing Vibration Performance Based on GM (1,1). J. Grey Syst. 2015, 27, 78–92. [Google Scholar]
  9. Dong, Y.; Zhang, Z.; Hong, W.-C. A Hybrid Seasonal Mechanism with a Chaotic Cuckoo Search Algorithm with a Support Vector Regression Model for Electric Load Forecasting. Energies 2018, 11, 1009. [Google Scholar] [CrossRef]
  10. Pham, H.T.; Yang, B.-S. Estimation and forecasting of machine health condition using ARMA/GARCH model. Mech. Syst. Signal Process. 2010, 24, 546–558. [Google Scholar] [CrossRef]
  11. Pavlovic, N.T.; Veg, A.; Troha, S.; Milovancevic, M.; Nikolic, V. Vibration prediction of pellet mills power transmission by artificial neural network. Assem. Autom. 2017, 37, 464–470. [Google Scholar]
  12. Fu, W.; Zhou, J.; Zhang, Y.; Zhu, W.; Xue, X.; Xu, Y. A state tendency measurement for a hydro-turbine generating unit based on aggregated EEMD and SVR. Meas. Sci. Technol. 2015, 26, 125008. [Google Scholar] [CrossRef]
  13. Javed, K.; Gouriveau, R.; Zerhouni, N.; Nectoux, P. Enabling Health Monitoring Approach Based on Vibration Data for Accurate Prognostics. IEEE Trans. Ind. Electron. 2015, 62, 647–656. [Google Scholar] [CrossRef] [Green Version]
  14. Fei, S.-W. Kurtosis forecasting of bearing vibration signal based on the hybrid model of empirical mode decomposition and RVM with artificial bee colony algorithm. Expert Syst. Appl. 2015, 42, 5011–5018. [Google Scholar] [CrossRef]
  15. Hong, W.-C.; Li, M.-W.; Geng, J.; Zhang, Y. Novel Chaotic Bat Algorithm for Forecasting Complex Motion of Floating Platforms. Appl. Math. Model. 2019, 72, 425–443. [Google Scholar] [CrossRef]
  16. Elsaid, A.; El Jamiy, F.; Higgins, J.; Wild, B.; Desell, T. Optimizing long short-term memory recurrent neural networks using ant colony optimization to predict turbine engine vibration. Appl. Soft Comput. 2018, 73, 969–991. [Google Scholar] [CrossRef]
  17. Fan, G.-F.; Peng, L.-L.; Hong, W.-C. Short term load forecasting based on phase space reconstruction algorithm and bi-square kernel regression model. Appl. Energy 2018, 224, 13–33. [Google Scholar] [CrossRef]
  18. Zhang, C.; Zhou, J.; Li, C.; Fu, W.; Peng, T. A compound structure of ELM based on feature selection and parameter optimization using hybrid backtracking search algorithm for wind speed forecasting. Energy Convers. Manag. 2017, 143, 360–376. [Google Scholar] [CrossRef]
  19. Chen, Z.; Wu, L.; Cheng, S.; Lin, P.; Wu, Y.; Lin, W. Intelligent fault diagnosis of photovoltaic arrays based on optimized kernel extreme learning machine and I-V characteristics. Appl. Energy 2017, 204, 912–931. [Google Scholar] [CrossRef]
  20. Vulli, S.; Dunne, J.; Potenza, R.; Richardson, D.; King, P. Time-frequency analysis of single-point engine-block vibration measurements for multiple excitation-event identification. J. Sound Vib. 2009, 321, 1129–1143. [Google Scholar] [CrossRef]
  21. Li, N.; Yang, J.; Zhou, R.; Wang, Q. Knock detection in spark ignition engines using a nonlinear wavelet transform of the engine cylinder head vibration signal. Meas. Sci. Technol. 2014, 25, 115002. [Google Scholar] [CrossRef]
  22. Huang, N.E.; Shen, Z.; Long, S.R.; Wu, M.C.; Shih, H.H.; Zheng, Q.; Yen, N.-C.; Tung, C.C.; Liu, H.H. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 1998, 454, 903–995. [Google Scholar] [CrossRef]
  23. Gilles, J. Empirical Wavelet Transform. IEEE Trans. Signal Process. 2013, 61, 3999–4010. [Google Scholar] [CrossRef]
  24. Liu, H.; Wu, H.; Li, Y. Smart wind speed forecasting using EWT decomposition, GWO evolutionary optimization, RELM learning and IEWT reconstruction. Energy Convers. Manag. 2018, 161, 266–283. [Google Scholar] [CrossRef]
  25. Hu, J.; Wang, J.; Ma, K. A hybrid technique for short-term wind speed prediction. Energy 2015, 81, 563–574. [Google Scholar] [CrossRef]
  26. Li, C.; Xiao, Z.; Xia, X.; Zou, W.; Zhang, C. A hybrid model based on synchronous optimisation for multi-step short-term wind speed forecasting. Appl. Energy 2018, 215, 131–144. [Google Scholar] [CrossRef]
  27. Ramsey, F.L. Characterization of the Partial Autocorrelation Function. Ann. Stat. 1974, 2, 1296–1301. [Google Scholar] [CrossRef]
  28. Polat, K.; Güneş, S. An expert system approach based on principal component analysis and adaptive neuro-fuzzy inference system to diagnosis of diabetes disease. Digit. Signal Process. 2007, 17, 702–710. [Google Scholar] [CrossRef]
  29. Abdoos, A.A. A new intelligent method based on combination of VMD and ELM for short term wind power forecasting. Neurocomputing 2016, 203, 111–120. [Google Scholar] [CrossRef]
  30. Xiao, L.; Shao, W.; Wang, C.; Zhang, K.; Lu, H. Research and application of a hybrid model based on multi-objective optimization for electrical load forecasting. Appl. Energy 2016, 180, 213–233. [Google Scholar] [CrossRef]
  31. Wang, J.; Heng, J.; Xiao, L.; Wang, C. Research and application of a combined model based on multi-objective optimization for multi-step ahead wind speed forecasting. Energy 2017, 125, 591–613. [Google Scholar] [CrossRef]
  32. Jiang, P.; Li, R.; Li, H. Multi-objective algorithm for the design of prediction intervals for wind power forecasting model. Appl. Math. Model. 2019, 67, 101–122. [Google Scholar] [CrossRef]
  33. Liu, Y.; Qin, H.; Mo, L.; Wang, Y.; Chen, D.; Pang, S.; Yin, X. Hierarchical Flood Operation Rules Optimization Using Multi-Objective Cultured Evolutionary Algorithm Based on Decomposition. Water Resour. Manag. 2019, 33, 337–354. [Google Scholar] [CrossRef]
  34. Galván, I.M.; Valls, J.M.; Cervantes, A.; Aler, R. Multi-objective evolutionary optimization of prediction intervals for solar energy forecasting with neural networks. Inf. Sci. 2017, 418, 363–382. [Google Scholar] [CrossRef]
  35. Mirjalili, S.; Gandomi, A.H.; Mirjalili, S.Z.; Saremi, S.; Faris, H.; Mirjalili, S.M. Salp Swarm Algorithm: A bio-inspired optimizer for engineering design problems. Adv. Eng. Softw. 2017, 114, 163–191. [Google Scholar] [CrossRef]
  36. Hu, Y.; Li, F.; Li, H.; Liu, C. An enhanced empirical wavelet transform for noisy and non-stationary signal processing. Digit. Signal Process. 2017, 60, 220–229. [Google Scholar] [CrossRef]
  37. Kedadouche, M.; Thomas, M.; Tahan, A. A comparative study between Empirical Wavelet Transforms and Empirical Mode Decomposition Methods: Application to bearing defect diagnosis. Mech. Syst. Signal Process. 2016, 81, 88–107. [Google Scholar] [CrossRef]
  38. Daubechies, I.; Heil, C. Ten Lectures on Wavelets. Comput. Phys. 1998, 6, 1671. [Google Scholar]
  39. Richman, J.S.; Lake, D.E.; Moorman, J.R. Sample Entropy. Methods Enzymol. 2004, 384, 172–184. [Google Scholar]
  40. Oussar, Y. Ranking a random feature for variable and feature selection. J. Mach. Learn. Res. 2003, 3, 1399–1414. [Google Scholar]
  41. Abdoos, A.A.; Mianaei, P.K.; Ghadikolaei, M.R. Combined VMD-SVM based feature selection method for classification of power quality events. Appl. Soft Comput. 2016, 38, 637–646. [Google Scholar] [CrossRef]
  42. Chen, H.-L.; Wang, G.; Ma, C.; Cai, Z.-N.; Liu, W.-B.; Wang, S.-J. An efficient hybrid kernel extreme learning machine approach for early diagnosis of Parkinson’s disease. Neurocomputing 2016, 184, 131–144. [Google Scholar] [CrossRef]
  43. Huang, G.-B.; Zhou, H.; Ding, X.; Zhang, R. Extreme Learning Machine for Regression and Multiclass Classification. IEEE Trans. Syst. Man Cybern. Part B 2012, 42, 513–529. [Google Scholar] [CrossRef] [Green Version]
  44. Lai, X.; Li, C.; Zhang, N.; Zhou, J. A multi-objective artificial sheep algorithm. Neural Comput. Appl. 2018, 1–35. [Google Scholar] [CrossRef]
  45. Zhang, C.; Peng, T.; Li, C.; Fu, W.; Xia, X.; Xue, X. Multiobjective Optimization of a Fractional-Order PID Controller for Pumped Turbine Governing System Using an Improved NSGA-III Algorithm under Multiworking Conditions. Complexity 2019, 2019, 1–18. [Google Scholar] [CrossRef]
  46. Zhou, J.; Xu, Y.; Zheng, Y.; Zhang, Y. Optimization of Guide Vane Closing Schemes of Pumped Storage Hydro Unit Using an Enhanced Multi-Objective Gravitational Search Algorithm. Energies 2017, 10, 911. [Google Scholar] [CrossRef]
  47. Abbassi, R.; Abbassi, A.; Heidari, A.A.; Mirjalili, S. An efficient salp swarm-inspired algorithm for parameters identification of photovoltaic cell models. Energy Convers. Manag. 2019, 179, 362–372. [Google Scholar] [CrossRef]
  48. Yang, X.-S. Bat Algorithm for Multi-objective Optimisation. Int. J. Bio-Inspired Comput. 2012, 3, 267–274. [Google Scholar] [CrossRef]
  49. Hyndman, R.J.; Koehler, A.B. Another look at measures of forecast accuracy. Int. J. Forecast. 2006, 22, 679–688. [Google Scholar] [CrossRef] [Green Version]
  50. BENTLY3300 Sensors. Available online: http://www.amy00.com/a/jszl/2017/63.html (accessed on 27 April 2019).
  51. Dhillon, J.; Parti, S.; Kothari, D. Stochastic economic emission load dispatch. Electr. Power Syst. Res. 1993, 26, 179–186. [Google Scholar] [CrossRef]
Figure 1. The framework of the proposed model.
Figure 1. The framework of the proposed model.
Sensors 19 02055 g001
Figure 2. Encoding strategy in the location vector in the multi-objective salp swarm algorithm (MOSSA).
Figure 2. Encoding strategy in the location vector in the multi-objective salp swarm algorithm (MOSSA).
Sensors 19 02055 g002
Figure 3. The structure of the hydropower generator unit (HGU) and the location of the sensors.
Figure 3. The structure of the hydropower generator unit (HGU) and the location of the sensors.
Sensors 19 02055 g003
Figure 4. The swing data of lower guide in X-direction.
Figure 4. The swing data of lower guide in X-direction.
Sensors 19 02055 g004
Figure 5. Modes of the vibration data with empirical wavelet transform (EWT).
Figure 5. Modes of the vibration data with empirical wavelet transform (EWT).
Sensors 19 02055 g005
Figure 6. RM obtained from EWT and SE-based reconstruction.
Figure 6. RM obtained from EWT and SE-based reconstruction.
Sensors 19 02055 g006
Figure 7. The best solution selection from Pareto-optimal sets of MOSSA.
Figure 7. The best solution selection from Pareto-optimal sets of MOSSA.
Sensors 19 02055 g007
Figure 8. Comparisons of vibration data and prediction results with different methods: (a) EEMD-KELM model, (b) EWT-KELM model, (c) AEWT-KELM model, (d) AEWT-PCA-KELM model, (e) AEWT-GSO-SVR model, (f) AEWT-GSO-SSA-KELM model, (g) AEWT-GSO-MOPSO-KELM model, and (h) AEWT-GSO-MOSSA-KELM model.
Figure 8. Comparisons of vibration data and prediction results with different methods: (a) EEMD-KELM model, (b) EWT-KELM model, (c) AEWT-KELM model, (d) AEWT-PCA-KELM model, (e) AEWT-GSO-SVR model, (f) AEWT-GSO-SSA-KELM model, (g) AEWT-GSO-MOPSO-KELM model, and (h) AEWT-GSO-MOSSA-KELM model.
Sensors 19 02055 g008aSensors 19 02055 g008bSensors 19 02055 g008c
Figure 9. The box plots of error distribution with different models.
Figure 9. The box plots of error distribution with different models.
Sensors 19 02055 g009
Table 1. Summary of literature on vibration tendency prediction.
Table 1. Summary of literature on vibration tendency prediction.
Forecasting ModelsMethodsData SetAuthorsReference
Statistical modelsAutoregressive moving average (ARMA) methodThe low methane compressorPham et al.[10]
Grey prediction methodRolling bearing vibrationXia et al.[8]
Artificial intelligence (AI) modelsArtificial neural network (ANN)The gear transmission vibration of pellet millsMilovancevic et al.[11]
Support Vector RegressionThe vibration trend of hydro-turbine generating unitFu et al.[12]
Extreme learning machinesThe vibration data of cutting tools and bearingJaved et al.[13]
Long short-term memory recurrent neural networksThe vibration of turbine engineEl Said et al.[16]
Hybrid modelsLS-SVR and chaotic sine cosine algorithm optimizationVibration trend of hydropower generatorFu et al.[4]
Empirical mode decomposition and relevance vector machineThe vibration signal of bearingsFei S.-W.[14]
Table 2. Detailed description of the monitoring data.
Table 2. Detailed description of the monitoring data.
SensorTimeTime IntervalThe Number of Samples
BENTLY330024–27 July 201110 min300
Table 3. The SE values of modes decomposed with EWT.
Table 3. The SE values of modes decomposed with EWT.
IndicatorMode 1Mode 2Mode 3Mode 4Mode 5Mode 6Mode 7Mode 8
SE0.00850.03810.08690.1290.1130.1680.2120.233
Table 4. Reconstructed RM derived by EWT according to SE.
Table 4. Reconstructed RM derived by EWT according to SE.
RMsModes ContainedHS
1Mode 1[0, 0.0085]
2Mode 2, Mode 3, [0.0308,0.0869]
3Mode 4, Mode 5, Mode 6[0.112, 0.168]
4Mode 7, Mode 8[0.177, 0.233]
Table 5. The performance of different models.
Table 5. The performance of different models.
ModelPrecision of Model PredictionComputing Time
RMSE (μm)MAE (μm)MAPE (%)RTime (s)
EEMD-KELM1.2361.0311.0930.827106.730
EWT-KELM1.2070.9971.0530.874105.138
AEWT-KELM1.1051.0271.1570.87998.876
AEWT-PCA-KELM0.9190.7980.8460.902100.263
AEWT-GSO-MOSSA-SVR0.8570.6410.6810.885109.114
AEWT-GSO-SSA-KELM0.9390.7430.7910.90695.716
AEWT-GSO-MOPSO-KELM0.8410.7040.7470.911109.987
AEWT-GSO-MOSSA-KELM0.8230.6500.6820.913102.920

Share and Cite

MDPI and ACS Style

Zhou, K.-B.; Zhang, J.-Y.; Shan, Y.; Ge, M.-F.; Ge, Z.-Y.; Cao, G.-N. A Hybrid Multi-Objective Optimization Model for Vibration Tendency Prediction of Hydropower Generators. Sensors 2019, 19, 2055. https://doi.org/10.3390/s19092055

AMA Style

Zhou K-B, Zhang J-Y, Shan Y, Ge M-F, Ge Z-Y, Cao G-N. A Hybrid Multi-Objective Optimization Model for Vibration Tendency Prediction of Hydropower Generators. Sensors. 2019; 19(9):2055. https://doi.org/10.3390/s19092055

Chicago/Turabian Style

Zhou, Kai-Bo, Jian-Yu Zhang, Yahui Shan, Ming-Feng Ge, Zi-Yue Ge, and Guan-Nan Cao. 2019. "A Hybrid Multi-Objective Optimization Model for Vibration Tendency Prediction of Hydropower Generators" Sensors 19, no. 9: 2055. https://doi.org/10.3390/s19092055

APA Style

Zhou, K. -B., Zhang, J. -Y., Shan, Y., Ge, M. -F., Ge, Z. -Y., & Cao, G. -N. (2019). A Hybrid Multi-Objective Optimization Model for Vibration Tendency Prediction of Hydropower Generators. Sensors, 19(9), 2055. https://doi.org/10.3390/s19092055

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop