Identification of Vibration Signal for Residual Pressure Utilization Hydraulic Unit Using MRFO-BP Neural Network

Combined with wavelet threshold denoising and Ensemble Empirical Mode Decomposition (EEMD) decomposition, an identication method based on Manta Ray Foraging Optimization-BP (MRFO-BP) neural network for vibration signals of residual pressure utilization hydraulic units is proposed to distinguish the vibration signal of each unit. e feature vectors of vibration signals are extracted by wavelet denoising and EEMD decomposition. e weights and thresholds in BP neural network are optimized by the MRFO algorithm. e feature vectors are input into the optimized BP neural network to realize the identication and classication of vibration signals. Compared with Particle Swarm Optimization-BP (PSO-BP) neural network, Bat Algorithm-BP (BA-BP) neural network, and BP neural network, the results show that the identication rate of each measuring point from the MRFO-BP neural network is greatly improved. e average identication rate of other measuring points is 98.514%, except the identication rate of the generator, which is 85.389%. erefore, the MRFO-BP neural network has better stability and higher identication accuracy and can identify and classify vibration signals more accurately. e conclusions can provide theoretical basis for vibration signals identication of residual pressure utilization hydraulic unit. When the vibration signal of each unit cannot be clearly distinguished, the vibration signals of the units are identied by the method proposed in this paper. According to the obtained results, a feasible classication method can be provided for the vibration signals belonging to dierent units.


Introduction
e residual pressure utilization hydraulic unit is used to recover the energy wasted in the process of throttling and pressure drop and is used for residual pressure power generation, pipe network optimization, e cient pump transformation, etc. Energy saving and e cient use are realized from the two aspects of open source and reducing expenditure. Vibration phenomenon is inevitable in the operation of residual pressure utilization hydraulic unit. Most faults of the unit are re ected in the vibration signal, so the vibration signal of the unit can directly re ect whether the unit runs safely and stably [1,2]. Comprehensive, accurate, and e ective identi cation of vibration signals is the premise of unit condition monitoring, and it is necessary to judge the stable operation of residual pressure utilization hydraulic unit [3].
In the identi cation of the vibration signal of the unit, the feature vector extraction of the signal is an indispensable part. In recent years, many scholars have proposed di erent methods in signal feature vector extraction and recognition combined with optimization algorithms. In order to promote the e ect of variational mode decomposition (VMD), Li et al. [4] proposed a feature extraction method based on GA-VMD and center frequency for the di culty of bearing fault feature extraction. In order to improve the ability of SlEn to distinguish di erent types of signals and solve the problem of two threshold parameters selection, Li et al. [5] proposed a new time series complexity indicator on the basis of SlEn by introducing fractional calculus and combining particle swarm optimization (PSO), named PSO fractional SlEn (PSO-FrSlEn). Aiming at the feature extraction and identi cation algorithm of AE signal for fatigue crack of pipeline structure, Wei et al. [6] proposed a fatigue crack identification method based on PNN combined with EMD feature extraction. In order to obtain short-term accurate power generation prediction, Hou et al. [7] used the Grey Wolf Optimization Algorithm (GWO) to optimize the initial weights and initial thresholds of the BP network. In order to predict the rock mass parameters of the TBM tunnel, Liu B. et al. [8] proposed a hybrid algorithm (SA-BPNN), which integrated the backpropagation neural network (BPNN) with simulated annealing (SA). Lee and Cheng [9] proposed a motor fault detection method based on wavelet transform (WT) and improved PSO-BP neural network for motor fault detection. Deepti Deshwal et al. [10] used hybrid robust feature extraction techniques for spoken language identification (LID) system and used forward backpropagation neural network (FFBPNN) for language identification or classification. Houssein et al. [11] proposed a new hybrid ECG arrhythmia classification method combining Manta Ray Foraging Optimization algorithm (MRFO) and Support Vector Machine (SVM). Yousri et al. [12] proposed a recombination and optimization method of the electrical photovoltaic array based on the Artificial Ecosystem-based Optimization (AEO) to solve the problem of harvesting the maximum power from a partially shaded photovoltaic array. Abdolrasol et al. [13] enhanced the neural network by optimizing the algorithm and obtained the best structural network pattern by manipulating its tuned parameters or training parameters. is paper includes some results of PSO, GA, ABC, and BSA optimization techniques to improve the performance of neural network. Xiong et al. [14] proposed a novel application based on Supply-Demandbased Optimization (SDO) to extract accurate and reliable parameters for photovoltaic models. In order to study the optimal control and operation of the power distribution system, Elattar et al. [15] applied the MRFO algorithm to IEEE 33-bus, 69-bus, and practical distribution network of 84-bus in the Taiwan Power Company (TPC). e results show that the MRFO algorithm is more effective and robust than other optimization techniques. Almagboul et al. [16] proposed a novel partially connected hybrid analog-digital receive beamformer based on the Atom Search Optimization algorithm (ASO), which is used to reduce the peak sidelobe level (SLL) and steer the nulls in the desired directions. Zhang and Li [17] used Particle Swarm Optimization (PSO) to optimize the BP neural network and used the trained neural network to predict the fault arc. Ismail et al. [18] coupled the PSO and BP algorithms and applied the algorithm to predict the load-deformation behavior of axially loaded piles, and the results showed that the proposed algorithms are more accurate than others. Shadman Abid et al. [19] solved the problem of RDG planning optimization by using Artificial Hummingbird Algorithm (AHA). arwat et al. [20] proposed Bat Algorithm (BA) to optimize the parameters of SVM, and the results of the BA-SVM algorithm are compared with those of the grid search, which is a conventional method of searching parameter values, and two well-known optimization algorithms (GA and PSO), and the results showed that the proposed model can find the optimal values of the SVM parameters.
In view of the nonlinear, nonstationary, and strong noise characteristics of the vibration signal of the unit, wavelet threshold denoising is used to denoise the vibration signal [21], which can suppress the noise interference and protect the useful information from being lost to the greatest extent. Vamsi et al. [22] used vibration data to monitor the damage of pultruded samples, and wavelet transform was used to process the vibration signals of healthy samples and damaged samples.
e Empirical Mode Decomposition (EMD) method does not need to set the wavelet basis function and the number of decomposition layers, so it is widely used in the field of processing nonlinear and nonstationary signals [23]. However, the EMD method has problems such as insufficient endpoint effect and obvious modal aliasing. Wu and Huang [24] proposed an Ensemble Empirical Mode Decomposition (EEMD). Adding Gaussian white noise to the signal to be measured can effectively reduce the problems of the EMD. Sun [25] researched the fault signal extraction and identification of vibration motor rolling bearings, and the results proved that EEMD is obviously stronger than EMD. e vibration signals are decomposed into Intrinsic Mode Function (IMF) components by EEMD. He et al. [26] used EEMD method to extract features of vibration signals of hydropower units, which verified the feasibility of this method for feature extraction of vibration signals. e advantages of manta ray foraging optimization algorithm (MRFO) are that it does not need to provide an accurate mathematical optimization model, and the optimization method of the objective function and constraint function is more relaxed. At the same time, it has strong global search ability and is not restricted by functions and continuous factors. When running the MRFO algorithm, it has the characteristics of strong optimization ability, fast convergence, and few parameters [27]. For the problems of unreasonable selection of initial weights and thresholds, poor network convergence, and prediction accuracy in BP neural network [7], MRFO algorithm can be used to optimize BP neural network. erefore, combined with wavelet threshold denoising and EEMD [28,29] feature extraction method, an analysis method named MRFO-BP neural network, which is based on the MRFO algorithm and BP neural network, is proposed to extract and identify the vibration signal features of residual pressure utilization hydraulic unit. e vibration signal is denoised by the wavelet threshold denoising method, and then the denoising vibration signal is decomposed by the EEMD method. e IMF components are extracted as the feature vectors and input into the optimized BP neural network for vibration signals identification. e MRFO-BP neural network is compared with PSO-BP neural network, BA-BP neural network, and BP neural network. e results show that the identification rate of MRFO-BP neural network is higher, and it can more effectively distinguish the vibration signals of different units in the same position.
is method can provide a theoretical basis for the subsequent identification of the vibration signals of residual pressure utilization hydraulic unit.  Figure 1. ere are five similar units, and the capacity is 200 kW. e runner is HL360A, and the speed is 600 r/min. ree units (1#, 3#, and 4#) are selected as the vibration test objects. e shaft of the 1# unit is shattered and just repaired. e 4# unit is running well. Compared with the 4# unit, the noise of the 3# unit is too big. e acceleration sensors are set in five parts of each unit as the measuring points, and the vertical and horizontal directions are arranged in coordination, as shown in Figures 2-4. e five measuring points are the generator, the inlet, the middle of the spiral casing, the tail of the spiral casing, and the outlet. e 10 measuring points are numbered from No. 1 to No. 10, respectively. e guide vanes opening of the units are selected as 40 mm.
DASP-V10 vibration data acquisition instrument is used to collect vibration signals. Ten models of INV9828 piezoelectric acceleration sensors are selected. e sampling frequency of vibration signal is 1024 Hz. e basic parameters of the acceleration sensors are shown in Table 1.
e experimental data of three residual pressure utilization hydraulic units are divided into five categories according to the different positions of measuring points. e specific classification is shown in Table 2. e three cases of each category are classified as 1# unit number 1, 3# unit number 2, and 4# unit number 3.

Feature Extraction Method.
e vibration frequency distribution of Gaussian white noise is uniform, which can effectively reduce the problem of insufficient endpoint effect and suppress the phenomenon of modal aliasing [30]. Based on the original EMD method, the EEMD method is proposed, which is realized by adding white noise with equal amplitude to the original signal for many times, and using EMD to decompose the signal after adding white noise [31,32]. In order to reduce the influence of white noise on the signal to be measured, the average value of the IMF components obtained by EMD decomposition is used as the final decomposition result. e EEMD method is selected to extract the features of the vibration signals of the residual pressure utilization hydraulic unit. Firstly, the wavelet threshold denoising method is applied to denoise the original signal, and then the denoised vibration signals are decomposed by the EEMD method, and the IMF components are extracted as the feature vectors to input the optimized BP neural network to identify the vibration signals. e feature vectors extraction process of the vibration signals is shown in Figure 4.

Vibration Signal Processing and Feature Extraction.
e wavelet threshold denoising method is used to denoise the vibration signals, and the "db1" wavelet function is selected, with three decomposition layers. Taking the vibration signals at the inlet of the three units as an example, the comparison waveform of the vibration signals before and after denoising is shown in Figure 5, with the abscissa as the sampling points and the ordinate as the amplitude (/mm). It can be seen from Figure 5 that the original vibration signals are disturbed by external noise, and feature extraction of the denoised signals can improve the accuracy of signal identification. e denoised vibration signals are decomposed by EEMD, and the IMF components are obtained as shown in Figure 6.
In order to unify the dimensions of the feature vectors, the feature vectors of different dimensions are substituted into the MRFO-BP neural network. After many times of training, IMF7-IMF12 are finally selected as the feature vector X � [X 1 , X 2 , X 3 , X 4 , X 5 , X 6 ] of each group of vibration signals.

Manta Ray Foraging
Optimization. e manta ray foraging optimization algorithm is a new intelligent bionic swarm algorithm proposed by Zhao et al. [33] in 2020. Manta rays are large marine creatures. ey usually appear in groups when foraging, and they feed on plankton with uneven or irregular distribution. Figure 7 depicts a foraging manta ray. e manta ray foraging optimization algorithm is a bionic optimization algorithm, which is inspired by the intelligent foraging strategy of manta ray. Manta rays usually appear in groups when foraging, and they have evolved a set of intelligent foraging strategies. e MRFO algorithm is obtained by mathematically modeling of the three predation strategies of manta rays. e MRFO algorithm searches for the optimal solution in the space by imitating chain foraging, spiral foraging, and flipping foraging behaviors of manta rays and achieves global optimization [34]. e first is chain foraging, in which the manta ray groups will line up one after the other. is arrangement is conducive to the former manta ray missing food, and the subsequent manta ray to make up for the former's mistake. e second is spiral foraging. When the manta rays detect that the food concentration is very high, the manta ray groups will gather together and connect their tails with their heads in a spiral shape, thus forming a spiral apex, and the filtered water will move up to the surface of the water. is will pull plankton into their open mouths. e third is flip foraging, which is an uncommon foraging behavior, but an extremely effective strategy. When the manta rays detect food, the manta rays will perform a series of random, frequent backflips centered on the food.

Chain Foraging.
During foraging, manta rays can observe the position of plankton and swim toward it, with the higher concentrations of plankton indicating the better location. Assuming that the position of high concentrations plankton indicates the optimal position, manta rays are connected end to end in a feeding chain. e first manta ray swims toward the optimal position, and the other manta rays swim not only toward the optimal position, but also toward Shock and Vibration the position of the previous manta ray. So, during the iteration process, each manta ray updates its position according to the current global optimal position and the previous individual position. e individual position update is given by the following equation:     where t and d represent the number of iterations and dimension, respectively. x d i (t) represents the individual position of the i generation. r represents a random number in the range of [0, 1]. x d best (t) represents the best position. α represents weight coefficient.

Spiral Foraging.
When the plankton is found in the deep water, the manta ray will move towards the food in a spiral form, and the individual manta ray will also follow the previous individual to update the position. e position update equation is as follows: where β represents the spiral weight coefficient. T represents the maximum number of iterations. r 1 represents a random In the experiment of measuring point 6, channel 11 is selected for 1# unit, channel 9 is selected for 3# unit and 4# unit, and the other parameters are completely the same.   (3) indicates that the MRFO algorithm searches near the optimal solution in the n-dimensional space and has good local search ability. e individual manta ray uses the food position as the reference position to search randomly, which can improve the global search ability of the algorithm and avoid falling into the local optimization. e expression is as follows: where x d rand (t) represents a random position randomly generated by the individual in the search space. Ub d and Lb d represent the upper and lower limits of the value, respectively.

Flip Foraging.
In the flipping foraging behavior, the manta rays flip to a new position with the food position as the fulcrum. e new position is any position in the search domain between the current position and the symmetrical position of the best position it has found so far. In the search space, as the distance between the individual position and the best position found so far decreases, all individuals gradually approach the optimal solution. e expression is as follows: where S represents the flip factor of the manta rays flip range. e size of S determines the individual flip distance of the manta rays.

Optimization of BP Neural Network by MRFO Algorithm.
e BP neural network [35] is a feed-forward neural network algorithm that conducts backward transmission training with error. It consists of input layer, hidden layer, and output layer, including feed-forward calculation process and error backpropagation process [36]. BP neural network is introduced to identify and classify the feature vectors obtained by EEMD decomposition.
Aiming at the unavoidable problems of BP neural network, such as easily falling into local minimum, slow convergence speed, and weak generalization ability, the MRFO algorithm is used to optimize BP neural network. e random initial population individuals within the boundary range are used as neural network parameters to assign initial values to the weights and thresholds of the network. en, the BP neural network is used as the fitness function, and the calculated error is used as the fitness value. e population individuals and fitness values are continuously updated through the MRFO algorithm. After many iterations, the optimal population individual and the global optimal fitness value are obtained. e optimal population individuals are substituted into the BP neural network to update the weights and thresholds, and the identification results of the residual pressure utilization hydraulic unit vibration signals are calculated. e identification process of MRFO-BP neural network is shown in Figure 8. e specific description is as follows: Step 1. Parameter initialization. Set the population size N, the maximum iteration number T, and the number of nodes in each layer of BP neural network; Step 2. Population initialization. Randomly initial population individuals, using the BP neural network as the fitness function to calculate the individual fitness value; Step 3. Generate a random number R 1 ; Step 4. If R 1 ≥ 0.5, update the position according to Equation (1); otherwise, generate a random number R 2 ; Step 5. If R 2 < t/T, update the position according to Equation (3); otherwise, update the position according to Equation (6); Step 6. Calculate the updated fitness value, and then update the position according to Equation (7); Step 7. Update the population and the global optimal position; Step 8. Judge whether the algorithm meets the termination condition. If so, output the optimal model of the BP neural network to identify the vibration signals. Otherwise, return to Step 3.

Training MRFO-BP Neural Network.
e parameters of MRFO-BP neural network are initialized: the initial population number is 30, the maximum number of iterations is 200, the number of nodes in the input layer is 6, the number of nodes in the hidden layer is 10, and the number of nodes in the output layer is 3.
Six groups of signals at the inlets of three units, each with 1000 IMF7-IMF12 components, are selected as the feature vectors of the vibration signals of residual pressure utilization hydraulic units. e feature vector sets are input into the MRFO-BP neural network for training and prediction. For each group of signals, the first 700 pieces of data are selected as training data, and the remaining 300 pieces of data are selected as the test data. e data are normalized to form a sample set. e training data set and initialization parameters are input into MRFO-BP neural network. After the BP neural network training and MRFO algorithm optimization, the optimal population and global optimal fitness value are obtained. e historical optimal fitness value curve of the inlet in the target space is shown in Figure 9.
It can be seen from Figure 9 that the initial fitness value at the inlet of the unit is 9084.238, which is maintained until the 12th iteration. e changes of the fitness value in the intermediate iteration process are as follows: the 13th to 6 Shock and Vibration It can be seen from Figure 10 that the initial fitness value at the generator of the unit is 8400.116, which is maintained until the 6th iteration. e fitness value is 8400.061 from the 7th to the 117th iteration. At the 118th iteration, the optimal fitness value is 8400.034.
It can be seen from Figure 11 that the initial fitness value in the middle of the spiral casing of the unit is 8818.357, which is maintained until the second iteration.
e changes of the fitness value in the intermediate iteration process are as follows: the 3rd to 42nd iterations are 8401.040, the 43rd to 49th iterations are 8400.524, the  Shock and Vibration 50th to 73rd iterations are 8400.406, the 74th to 101st iterations are 8400.167, the 102nd iteration is 8400.165, and the 103rd to 149th iterations are 8400.146. At the 150th iteration, the optimal fitness value is 8400.071. It can be seen from Figure 12 that the initial fitness value in the tail of the spiral casing of the unit is 8400.506, which is maintained until the 7th iteration. e changes of the fitness value in the intermediate iteration process are as follows: the 8th to 27th iterations are 8400.380, the 28th to 67th iterations are 8400.351, and the 68th to 72nd iterations are 8400.250. At the 73rd iteration, the optimal fitness value is 8400.072.
It can be seen from Figure 13 that the initial fitness value at the outlet of the unit is 8475. 034    to 8th iterations are 8400.684, the 9th to 26th iterations are 8400.309, the 27th to 124th iterations are 8400.160, the 125th to 175th iterations are 8400.151, and the 176th to 178th iterations are 8400.106. At the 179th iteration, the optimal fitness value is 8400.072.

MRFO-BP Neural Network
Identification. e weights and thresholds of the BP neural network are optimized by using the optimal population. en, 1000 pieces of data for each of the 6 groups of signals at the inlet of the three residual pressure utilization hydraulic units are divided into  the first 700 data as training data and the remaining 300 data as test data, which are, respectively, input into the BP neural network. After training and prediction, the identification rate of vibration signals is obtained through comparing the prediction results with each group of classification labels. In the MRFO-BP neural network, the identification rate of the three units at the inlet is 99.167%, and the image of the identification rate is shown in Figure 14.
In Figure 14, at the inlet, the 600 predicted values of the 1# unit completely coincide with the marker 1, the 600 predicted values of the 4# unit completely coincide with the marker 3, and the probability that the 600 predicted values of the 3# unit do not coincide with the marker 2 is 0.833%.

Shock and Vibration
Based on the MRFO-BP neural network, the vibration signals of the four measuring points of the three units are trained and predicted, respectively. e obtained identification rates are shown in Table 3, and the identification rate images are shown in Figure 15. e four measuring points are the generator, the middle of the spiral casing, the tail of the spiral casing, and the outlet.
In Figure 15, at the generator, the 600 predicted values of the 1# unit completely coincide with marker 1, the 600 predicted values of the 4# unit completely coincide with marker 3, and the probability that the 600 predicted values of the 3# unit do not coincide with marker 2 is 14.611%. In the middle of the spiral casing, the 600 predicted values of the 1# unit completely coincide with marker 1, the 600 predicted values of the 4# unit completely coincide with marker 3, and the probability that the 600 predicted values of the 3# unit do not coincide with marker 2 is 1.167%. In the tail of the spiral casing, the 600 predicted values of the 3# unit completely coincide with marker 2, the 600 predicted values of the 4# unit completely coincide with marker 3, and the probability that the 600 predicted values of the 1# unit do not coincide with marker 1 is 3.944%. At the outlet, the 600 predicted values of the 1# unit completely coincide with marker 1, the 600 predicted values of the 3# unit completely coincide with marker 2, and the 600 predicted values of the 4# unit completely coincide with marker 3. erefore, except for the generator, the average probability that the predicted values of other measuring points do not coincide with the classification markers is 1.981%.

Algorithm Comparison
In order to further verify the superiority of the MRFO-BP neural network, it is compared with PSO-BP neural network, BA-BP neural network, and BP neural network, respectively.

BP Neural Network Identification.
e same data set is directly input into the BP neural network. At the same time, it is initialized according to the parameter settings of the MRFO-BP neural network. e identification rates of the three units at each measuring point are shown in Table 4, and the identification rate images are shown in Figure 16.
In Figure 16, at the generator, the 600 predicted values of the 4# unit completely coincide with marker 3, and the probability that the 600 predicted values of the 1# unit and the 3# unit, respectively, do not coincide with marker     1 and marker 2 is 31.722%. At the inlet, the 600 predicted values of the 1# unit completely coincide with marker 1, the 600 predicted values of the 4# unit completely coincide with marker 3, and the probability that the 600 predicted values of the 3# unit do not coincide with marker 2 is 9.722%. In the middle of the spiral casing, the 600 predicted values of the 1# unit completely coincide with marker 1, the 600 predicted values of the 4# unit completely coincide with marker 3, and the probability that the 600 predicted values of the 3# unit do not coincide with marker 2 is 12.278%. In the tail of the spiral casing, the 600 predicted values of the 3# unit completely coincide with marker 2, the 600 predicted values of the 4# unit completely coincide with marker 3, and the probability that the 600 predicted values of the 1# unit do not coincide with marker 1 is 11.778%. At the outlet, the 600 predicted values of the 3# unit completely coincide with marker 2, the 600 predicted values of the 4# unit completely coincide with marker 3, and the probability that the 600 predicted values of the 1# unit do not coincide with marker 1 is 11.389%. erefore, except for the generator, the average probability that the predicted values of other measuring points do not coincide with the classification markers is 11.292%.

PSO-BP Neural Network Identification.
e PSO-BP neural network uses particle swarm algorithm (PSO) to optimize the BP neural network. e same data set is input into the PSO-BP neural network. At the same time, it is initialized according to the parameter settings of the MRFO-BP neural network. rough optimization, training, and prediction, the obtained prediction results are compared with each group of category labels. e obtained vibration signal identification rates of the three units at each measuring point are shown in Table 5, and the identification rate images are shown in Figure 17.
In Figure 17, at the generator, the 600 predicted values of the 1# unit completely coincide with marker 1, the 600 predicted values of the 4# unit completely coincide with marker 3, and the probability that the 600 predicted values of the 3# unit do not coincide with marker 2 is 25.556%. At the inlet, the 600 predicted values of the 1# unit completely coincide with marker 1, the 600 predicted values of the 4# unit completely coincide with marker 3, and the probability that the 600 predicted values of the 3# unit do not coincide with marker 2 is 5.333%. In the middle of the spiral casing, the 600 predicted values of the 4# unit completely coincide with marker 3, and the probability that the 600 predicted values of the 1# unit and the 3# unit,   respectively, do not coincide with marker 1 and marker 2 is 7.611%. In the tail of the spiral casing, the 600 predicted values of the 3# unit completely coincide with marker 2, the 600 predicted values of the 4# unit completely coincide with marker 3, and the probability that the 600 predicted values of the 1# unit do not coincide with marker 1 is 9.500%. At the outlet, the 600 predicted values of the 3# unit completely coincide with marker 2, the 600 predicted values of the 4# unit completely coincide with marker 3, and the probability that the 600 predicted values of the 1# unit do not coincide with marker 1 is 6.111%. erefore, except for the generator, the average probability that the predicted values of other measuring points do not coincide with the classification marks is 7.139%.

BA-BP Neural Network
Identification. e BA-BP neural network uses Bat algorithm (BA) to optimize the BP neural network. e same data set is input into the BA-BP neural network. At the same time, it is initialized according to the parameter settings of the MRFO-BP neural network. rough optimization, training, and prediction, the obtained prediction results are compared with each group of category labels. e obtained vibration signal identification rates of the three units at each measuring point are shown in Table 6, and the identification rate images are shown in Figure 18.
In Figure 18, at the generator, the 600 predicted values of the 1# unit completely coincide with marker 1, the 600 predicted values of the 4# unit completely coincide with marker 3, and the probability that the 600 predicted values of the 3# unit do not coincide with marker 2 is 22.444%. At the inlet, the 600 predicted values of the 1# unit completely coincide with marker 1, the 600 predicted values of the 4# unit completely coincide with marker 3, and the probability that the 600 predicted values of the 3# unit do not coincide with marker 2 is 5.000%. In the middle of the spiral casing, the 600 predicted values of the 1# unit completely coincide with marker 1, the 600 predicted values of the 4# unit completely coincide with marker 3, and the probability that the 600 predicted values of the 3# unit do not coincide with     Table 7.
According to the comparison in Table 7, although the identification rate of the vibration signals of the three units at the generator is lower than that of other measuring points, the identification rate of each measuring point of MRFO-BP neural network is greatly improved, which is much higher than that of the PSO-BP neural network, the BA-BP neural network, and the BP neural network. e analysis shows that the MRFO-BP neural network has obvious advantages in the

Shock and Vibration 29
identification of vibration signals of the residual pressure utilization hydraulic units.

Comparative Experiment
For the comparative experiment, 1# and 3# units with a guide vane opening of 30 mm and 4# unit with a guide vane opening of 33 mm were selected. Except for the different data of the unit used, the classification, processing, training, and identification of the data are the same as those of the unit with a guide vane opening of 40 mm.

Training MRFO-BP Neural Network.
In the comparative experiment, the wavelet threshold denoising method is used to denoise the vibration signal, and EEMD method is used to extract the features of the denoised signal. IMF7-IMF12 are selected as the feature vector X � [X 1 , X 2 , X 3 , X 4 , X 5 , X 6 ] of each group of vibration signals. e parameters of MRFO-BP neural network are initialized: the initial population number is 30, the maximum number of iterations is 200, the number of nodes in the input layer is 6, the number of nodes in the hidden layer is 10, and the number of nodes in the output layer is 3.
Six groups of signals from three units at each measuring point, each with 1000 IMF7-IMF12 components, are selected as the feature vectors of the vibration signals of residual pressure utilization hydraulic units. For each group of signals, the first 700 pieces of data are selected as training data, and the remaining 300 pieces of data are selected as the test data. e data are normalized, and the training data set and initialization parameters are input into the MRFO-BP neural network. After the BP neural network training and MRFO algorithm optimization, the optimal population and global optimal fitness value are obtained. e historical optimal fitness value curve of each measuring point in the target space is shown in Figures 19-23. e five measuring points are the generator, the inlet, the middle of the spiral casing, the tail of the spiral casing, and the outlet.
It can be seen from Figure 19 that the initial fitness value at the generator of the unit is 9462.632, which is maintained until the 9th iteration. e changes of the fitness value in the intermediate iteration process are as follows: the 10th to 76th iterations are 8400.241, the 77th to 80th iterations are 8400.217, and the 81st to 114th iterations are 8400.105. At the 115th iteration, the optimal fitness value is 8400.085.
It can be seen from Figure 20 that the initial fitness value at the inlet of the unit is 10353.948, which is maintained until the 5th iteration.
e changes of the fitness value in the intermediate iteration process are as follows: the 6th to 11th iterations are 10062.538, the 12th to 23rd iterations are 9925.532, the 24th to 98th iterations are 8961.488, and the 99th to 104th iterations are 8895.809. At the 180th iteration, the optimal fitness value is 8591.741.
It can be seen from Figure 21 that the initial fitness value in the middle of the spiral casing of the unit is 8400.455, which is maintained until the second iteration. e changes   of the fitness value in the intermediate iteration process are as follows: the 3rd to 5th iterations are 8400.453, the 6th iteration is 8400.276, the 7th to 43rd iterations are 8400.182, and the 44th to 146th iterations are 8400.149. At the 147th iteration, the optimal fitness value is 8400.142. It can be seen from Figure 22 that the initial fitness value in the tail of the spiral casing of the unit is 8400.375, which is maintained until the 4th iteration. e fitness value is 8400.228 from the 5th to the 40th iteration. At the 41st iteration, the optimal fitness value is 8400.078.
It can be seen from Figure 23 that the initial fitness value at the outlet of the unit is 8475.034, which is maintained until the 5th iteration. e changes of the fitness value in the intermediate iteration process are as follows: the 6th to 9th iterations are 8400.106, the 10th to 147th iterations are 8400.025, and the 148th to 191st iterations are 8400.0244. At the 192nd iteration, the optimal fitness value is 8400.0238.

MRFO-BP Neural Network Identification.
In the comparative experiment, the weights and thresholds of the BP neural network are optimized by using the optimal population. en, the training data sets and test data sets of the three units at each measuring point are, respectively, input into the BP neural network. After training and prediction, the identification rates of vibration signals are obtained through comparing the prediction results with each group of classification labels. In the MRFO-BP neural network, the identification rates of the three units at each measuring point are shown in Table 8, and the identification rate images are shown in Figure 24.
In Figure 24, at the generator, the 600 predicted values of the 3# unit completely coincide with marker 2, the 600 predicted values of the 4# unit completely coincide with marker 3, and the probability that the 600 predicted values of the 1# unit do not coincide with marker 1 is 14.278%. At the inlet, the 600 predicted values of the 3# unit completely coincide with marker 2, the 600 predicted values of the 4# unit completely coincide with marker 3, and the probability that the 600 predicted values of the 1# unit do not coincide with marker 1 is 1.5%. In the middle of the spiral casing, the 600 predicted values of the 1# unit completely coincide with marker 1, the 600 predicted values of the 4# unit completely coincide with marker 3, and the probability that the 600 predicted values of the 3# unit do not coincide with marker 2 is 3.167%. In the tail of the spiral casing, the 600 predicted values of the 1# unit completely coincide with marker 1, the 600 predicted values of the 4# unit completely coincide with marker 3, and the probability that the 600 predicted values of the 3# unit do not coincide with marker 2 is 4.222%. At the outlet, the 600 predicted values of the 3# unit completely coincide with marker 2, the 600 predicted values of the 4# unit completely coincide with marker 3, and the probability that the 600 predicted values of the 1# unit do not coincide with marker 1 is 1.389%. erefore, except for the generator, the average probability that the predicted values of other measuring points do not coincide with the classification markers is 2.570%.

Algorithm Comparison.
e same data sets are input into BP neural network, PSO-BP neural network, and BA-BP neural network, respectively. At the same time, it is initialized according to the parameter settings of the MRFO-BP neural network.
rough optimization, training, and prediction, the obtained prediction results are compared with each group of category labels. In each algorithm, the vibration signal identification rate of three units at each measuring point is obtained. e PSO-BP neural network, BA-BP neural network, and BP neural network are compared with MRFO-BP neural network. e vibration signal identification rates of three units in each measuring point of four neural networks are shown in Table 9.
In the comparative experiment, it can be obtained from Table 9 that although the identification rate of the vibration signals of the three units at the generator is lower than that of other measuring points, the identification rate of each measuring point of MRFO-BP neural network is greatly improved, which is much higher than that of the PSO-BP neural network, the BA-BP neural network, and the BP neural network. e analysis shows that the MRFO-BP neural network has obvious advantages in the identification of vibration signals of the residual pressure utilization hydraulic units.

Conclusions
Combined with wavelet threshold denoising and EEMD decomposition, an analysis method based on MRFO-BP neural network is proposed for feature extraction and identification of vibration signals of residual pressure utilization hydraulic units. After the optimization, training, and identification of MRFO-BP neural network, combined with the results of the comparative experiment, the original experimental conclusions are as follows: Wavelet threshold denoising and EEMD decomposition are performed on the vibration signals of the residual pressure utilization hydraulic unit collected on-site, and the feature vectors are extracted and constructed. e feature vector is input into the MRFO-BP neural network recognition model, and the identification rates of the three units at different measuring points are obtained.
e average identification rate of other measuring points is 98.514%, except 85.389% at the generator. e identification rates of MRFO-BP neural network are compared with PSO-BP neural network, BA-BP neural network, and BP neural network. e identification rate of the MRFO-BP neural network at different measuring points: the identification rate of the generator is 85.389%, and the average identification rate of other measuring points is 98.514%. e identification rate of the PSO-BP neural network at different measuring points: the identification rate of the generator is 74.444%, and the average identification rate of other measuring points is 92.861%. e identification rate of the BA-BP neural network at different measuring points: the identification rate of the generator is 77.556%, and the average identification rate of other measuring points is 93.903%. e identification rate of the BP neural network at different measuring points: the identification rate of the generator is 68.278%, and the average identification rate of other measuring points is 88.708%. It can be concluded that the MRFO-BP neural network has higher identification accuracy, except that the lifting amplitude of generator is 10.945% (PSO-BP), 7.833% (BA-BP), and 17.111% (BP neural network), and the average lifting amplitude of other measuring points is 5.653% (PSO-BP) and 4.611%, respectively. erefore, the identification accuracy of vibration signals of residual pressure utilization hydraulic units obtained by MRFO-BP neural network is significantly improved.
Data Availability e data are derived from field experiments and involves confidentiality so are not shared.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.