Hydraulic-Pump Fault-Diagnosis Method Based on Mean Spectrogram Bar Graph of Voiceprint and ResNet-50 Model Transfer

<jats:p>The vibration signal of a pump is often used for analysis in the study of hydraulic-pump fault diagnosis methods. In this study, for the analysis, sound signals were used, which can be used to acquire data in a non-contact manner to expand the use scenarios of hydraulic-pump fault-diagnosis methods. First, the original data are denoised using complete ensemble empirical mode decomposition with adaptive noise and the minimum redundancy maximum relevance algorithm. Second, the noise-reduced data are plotted as mean spectrogram bar graphs, and the datasets are divided. Third, the training set graphs are input into the ResNet-50 network to train the base model for fault diagnosis. Fourth, all the layers of the base model are frozen, except for the fully connected and softmax layers, and the support set graphs are used to train the base model through transfer learning. Finally, a fault diagnosis model is obtained. The model is tested using data from two test pumps, resulting in accuracies of 86.1% and 90.8% and providing evidence for the effectiveness of the proposed method for diagnosing faults in hydraulic plunger pumps.</jats:p>

section_content: The hydraulic pump is the power component of the hydraulic system, and the plunger pump has advantages such as a high rated pressure, compact structure, and high volumetric efficiency compared to other types of pumps. Therefore, they are widely used in heavy machinery, national defense equipment, and other fields. Plunger pumps typically operate under high-speed and heavy-load conditions, resulting in a shorter service life than other types of pumps. To prevent the occurrence of serious safety incidents, reduce hydraulic system maintenance costs, and shorten the maintenance duration, it is necessary to conduct research on fault diagnosis methods for plunger pumps. 

Currently, vibration signals are commonly used for the fault diagnosis of rotating mechanical components such as bearings, gearboxes, and plunger pumps, and good research results have been obtained [1] [2] [3] [4]. However, the collection of vibration signals is contactbased, which is inconvenient or impossible to achieve in some practical applications. In contrast, sound signals, which are generated by vibrations, not only contain rich equipment status information but can also be acquired in a non-contact-based manner, with less influence from spatial restrictions. The use of sound signals for fault diagnosis will significantly expand the application scope of fault diagnosis technologies. In addition, the use of a single sound-level meter to sample sound signals can aid in diagnosing faults in multiple components of equipment, while the collection of vibration signals often requires multiple sensors. Therefore, in this study, sound signals were analyzed, and a fault diagnosis was conducted for a constant-pressure variable-displacement axial plunger pump. 

Traditional fault diagnosis techniques usually require complex signal processing tasks, such as noise reduction, filtering, and feature extraction, which require technical expertise and knowledge for implementation. However, with the development of deep-learning technology and the advent of the big data era, feature-based approaches that rely on manual selection will gradually be replaced by deep-learning algorithms. Deep-learningbased methods can adaptively extract features from equipment status data. These features not only contain well-known equipment information but also potential information that can help in identifying and locating faults [5]. Consequently, an increasing number of researchers have been interested in applying deep-learning methods in the field of fault diagnosis. In 2017, Qi Y. M. et al. proposed a fault diagnosis method based on a stacked sparse autoencoder and applied it to the diagnosis of rotating machinery faults [6]. In the same year, Verstraete D. et al. generated images from the time-frequency information of the raw data and inputted them into a deep convolutional neural network (CNN) for fault diagnosis. This method was tested using a publicly available rolling-bearing dataset [7]. In 2019, Mao W. T. et al. proposed a fault diagnosis method based on a generative adversarial network (GAN), considering the data imbalance problem in fault diagnosis. The proposed method was demonstrated to address the data-imbalance fault-diagnosis problem in bearing datasets [8]. In the same year, Peng D. D. et al. proposed a one-dimensional deep CNN based on a new residual block and tested it on a rolling-bearing dataset, thus achieving a good diagnostic performance [9]. In 2020, Li X. Q. et al. proposed an enhanced selective ensemble deep learning method and used a beetle antennae search algorithm to improve the performance of a rolling-bearing fault diagnosis algorithm. The experimental results showed that this method is more accurate and robust than baseline models and other ensemble learning methods [10]. In 2021, Li T. F. et al. proposed an effective intelligent fault diagnosis method based on a multireceptive field-graph convolutional network. This method has the advantages of strong data relationship mining and feature representation effects. By converting data samples into weighted graphs and extracting and fusing features from multiple receptive fields, this method solves the limitations of existing graph convolutional networks in terms of weights and receptive fields [11]. In 2022, Tang S. N. et al. used a Bayesian optimization algorithm to automatically select the hyperparameters of CNN models for the intelligent fault diagnosis of hydraulic axial plunger pumps. The results indicated that the proposed method provides excellent performance [12]. 

ResNet-50, which is used in this study, is a deep learning structure proposed by He K. M. et al., and along with it are other models with different network depths, such as ResNet-18, ResNet-34, ResNet-101, and ResNet-152 [13]. These models are composed of residual blocks and are therefore also known as deep residual networks (DRNs). Normalized initialization and intermediate normalization can largely solve the problems of vanishing gradients and exploding gradients caused by the increasing network depth. The DRN is designed to solve this degradation problem. The degradation problem is not overfitting, which refers to the phenomenon in which, as the network depth increases, the training accuracy quickly decreases after reaching saturation; in contrast, overfitting results in extremely high accuracy in the training set but low accuracy in the testing set. 

However, deep learning has significant limitations, which can be summarized as follows. Firstly, the effectiveness of deep learning largely depends on large-scale samples, and the lack of prior knowledge may lead to results that deviate from human senses or expert knowledge. Secondly, deep learning is essentially mapping-i. e., a feature relationship between the input and output-and it lacks causal reasoning. Thirdly, deep learning lacks interpretability, as it is an end-to-end model that contains numerous neurons and parameters, and people cannot explain the meaning of each parameter clearly, which is also one of the greatest shortcomings of deep learning [14]. 

With the development of deep learning, the serious problem of dependence on large sample sizes has become apparent. Researchers have attempted to use models trained for a specific problem to solve a different but related problem. This learning method is known as transfer learning. Compared with training a separate model for the second problem, this approach can reduce the data and time of the required training. In 1993, Pratt L. proposed a discriminability-based transfer algorithm that can achieve transfer between neural networks and described in detail how to use transfer learning to leverage existing neural networks [15]. In 2020, Chen Z. Y. et al. proposed a transferable CNN to improve target tasks and applied this network to the fault diagnosis of rotating machinery [16]. In 2021, Deng Y. F. et al. proposed a method based on a double-layer attention-based GAN (DA-GAN) to address the problem of local transfer in the field of mechanical fault diagnosis. The DA-GAN method constructs two attention matrices that guide the model to focus on or ignore certain data components before domain adaptation [17]. In 2022, Wang Z. J. et al. proposed a fault diagnosis model called the subdomain adaptation transfer learning network, which not only extracts transferable features but also constructs target subdomains for each category using pseudo-label learning. It also reduces the bias of marginal and conditional distributions and adaptively adjusts the contribution of each network layer using dynamic weight terms [18]. 

The aim of this study is to investigate the constant-pressure variable-displacement axial plunger pump, and the noise reduction of the pump's sound signals is achieved through complete ensemble empirical mode decomposition with adaptive noise (CEEMDAN) and maximum relevance minimum redundancy (mRMR) algorithms. The denoised signals are then plotted as mean spectrogram bar graphs. Subsequently, the ResNet-50 model is used to extract features and classify the graphs generated in the previous step, thus completing the training of the source domain model for the fault diagnosis of a single pump. Finally, through transfer learning, the final diagnostic model is generalized and is able to diagnose faults in multiple pumps with minimal training. 

The first part of this article is the introduction, which introduces the research content, research background, and research purpose. The second part is the theoretical background, which details the relevant information and mathematical derivation of the proposed algorithm, including CEEMDAN, the generation of the mean spectrogram bar graph, and the transfer learning algorithm based on the ResNet-50 model. The third part comprises the test and processes, wherein the pump fault simulation test and the dataset generated from it are introduced. The data are then input into the diagnostic algorithm, and the results of each step are presented and analyzed. The fourth part comprises the conclusions, which summarize this article and highlight the advantages of the proposed method. Finally, an outlook for future research on this topic is presented.
section_name: Introduction

section_name: Theoretical Background

section_content: In 1998, Huang N. E. et al. proposed an empirical mode decomposition (EMD) algorithm, which can adaptively decompose signals into a finite number of intrinsic mode functions (IMFs) and is suitable for analyzing non-linear and non-stationary signals [19]. In 2009, Wu Z. et al. proposed an ensemble empirical mode decomposition (EEMD) algorithm to solve the endpoint effect and mode mixing problems in EMD by adding auxiliary white noise [20]. In 2011, Torres M. E. et al. proposed the CEEMDAN algorithm based on EEMD, which eliminates the noise residue caused by adding auxiliary white noise [21]. The EMD serves as the foundation for the EEMDAN algorithm; therefore, it is necessary to introduce the EMD before discussing the CEEMDAN.
section_name: CEEMDAN

section_content: The basic concept of the EMD algorithm is to decompose a signal sequence until the decomposition-stopping conditions are satisfied. The conditions are as follows: (1) The number of extreme points in the current component is less than two and (2) the remainder of the current component cannot be decomposed further. A flowchart of the EMD algorithm is presented in Figure 1.
section_name: EMD

section_content: The basic concept of the EMD algorithm is to decompose a signal sequence until the decomposition-stopping conditions are satisfied. The conditions are as follows: (1) The number of extreme points in the current component is less than two and (2) the remainder of the current component cannot be decomposed further. A flowchart of the EMD algorithm is presented in Figure 1. The specific calculation steps of the EMD algorithm are as follows:

Step 1: All the local maximum and minimum points of the original signal sequence ( ) x t are found, and spline functions are used to fit the upper envelope 1 ( ) l t  and lower The specific calculation steps of the EMD algorithm are as follows:

Step 1: All the local maximum and minimum points of the original signal sequence x(t) are found, and spline functions are used to fit the upper envelope l 1 (t) and lower envelope l ˜1(t). The average line l 1 (t) of the upper and lower envelopes are then calculated:

Step 2: The average line l 1 (t) is subtracted from the original signal sequence x(t) to obtain the new sequence h 1 (t):

Typically, h 1 (t) is not considered the first IMF component. Therefore, it is necessary to repeat the above steps k times until the average line l 1k (t) tends to zero and obtain the first IMF component c 1 (t):

Step 3: c 1 (t) is subtracted from the original sequence x(t), which results in the residual sequence r 1 (t):

The second IMF component c 2 (t) is obtained by repeating the first and second steps for r 1 (t). This process is repeated until the n-th residual sequence r n (t), which cannot be decomposed, is obtained as:

At this point, the EMD decomposition process is completed, and r n (t) is called the residue. The original signal x(t) can be reconstructed by the various IMF components and residue:

Since the emergence of the EMD algorithm, it has been widely applied in several nonlinear fields. However, it also suffers from the problems of endpoint effects and mode mixing. Methods such as the mirror method [22], the extrema extension method [23], the parallel extension method [24], and the boundary local characteristic scale extension method [25] have been developed to address the endpoint effects. To address the problem of mode mixing, Wu Z. et al. proposed the EEMD based on EMD, which solves the aforementioned problem by introducing Gaussian white noise, thus making the decomposition results more stable and reliable. However, because of the repeated addition of Gaussian white noise, the reconstructed signal contains more residual noise, which increases the number of IMF components in the decomposition results. To address this problem, Torres M. E. improved the EEMD method and proposed the CEEMDAN algorithm.
section_name: EMD

section_content: As in the case of the EEMD, the CEEMDAN adds white noise to the original signal x(t). However, in contrast to EEMD, after the first addition of the complete Gaussian white noise, CEEMDAN adds the IMF components of the white noise obtained using EMD. A flowchart of the CEEMDAN algorithm is presented in Figure 2. 

The specific steps of the CEEMDAN algorithm are as follows:

Step 1: The calculation parameters are initialized, including the noise standard deviation ε i the and number of EMD executions l. 

Step 2: Different Gaussian white noises n i (t) are added to the original signal x(t) to obtain a new signal x i (t):

where i = 1, 2,. . ., l. 

Step 3: The EMD is executed l times on x i (t), the average of l IMF components is calculated, and the first IMF component of the CEEMDAN c 1 (t) is obtained:

Figure 2. Flowchart of the CEEMDAN algorithm. 

The specific steps of the CEEMDAN algorithm are as follows:

Step 1: The calculation parameters are initialized, including the noise standard deviation i  the and number of EMD executions l. 

Step 2: Different Gaussian white noises ( ) i n t are added to the original signal ( ) Step 4: EMD is executed on l Gaussian white noises n i (t). The first IMF components are added to r 1 (t), and the l resulting sequences are subjected to EMD. The average of the first components is the second IMF component c 2 (t):

where E 1 (•) denotes the first IMF component of the sequence in parentheses. Subtracting c 2 (t) from r 1 (t) results in the second residue r 2 (t):

Step 5:

Step 4 is repeated to calculate the remaining IMF components c k (t) and residue r k (t):

where k = 3, 4,. . ., K. The K-th residue cannot be further decomposed. The original signal x(t) is reconstructed from the obtained IMF components and residues:

The two important parameters in the CEEMDAN algorithm are the number of iterations l and the noise standard deviation ε i. The number of executions l is generally set to a large value. Torres M. E. et al. set l as 500 in their study [21]. The standard deviation of the added noise ε i can be used to select the signal-to-noise ratio (SNR) at each step. Wu Z. et al. suggested that adding small-amplitude noise signals during the decomposition process is beneficial [20]. As a result, ε i is usually set to a small value, such as 0. 02 [21]. In this study, all ε i are set as the same value.
section_name: CEEMDAN

section_content: The spectrogram was invented by researchers at Bell Laboratories in 1941. This is an image that displays sounds in three dimensions. A spectrogram typically presents the intensity of sound using colors, with time along the horizontal axis and frequency along the vertical axis. The spectrogram provides a visual representation of how the different frequency components of a sound signal change over time, thus making it easy to observe the frequency content of the signal at different moments. In 1945, Kesta L. G. et al. completed spectrogram matching using visual observation and first proposed the concept of a voiceprint. Spectrograms can be classified into two types: Wideband and narrowband spectrograms. 

The wideband spectrogram is based on the windowed Fourier Transform technique. It divides a long signal into several time segments of equal length, performs a short-time Fourier transform (STFT) on each segment, and then concatenates the multiple frames of the spectra obtained by the STFT in the time dimension with overlap to form an image. Typically, the wideband spectrogram uses a wider frequency range and shorter time windows with shorter intervals between adjacent windows. Therefore, it allows for the observation of the transient properties of the signal in the time domain, thus presenting a more distinct visualization of high-frequency components and transient signals and providing a high temporal resolution. 

Narrowband spectrograms are similar to wideband spectrograms, but usually select a relatively narrow frequency band with a high-frequency resolution, thus making them more suitable for analyzing the frequency domain details of audio signals. Considering that the state information of the equipment is more clearly reflected in the frequency domain, in this study, narrowband spectrograms are used for the analysis.
section_name: Spectrogram and Mean Spectrogram Bar Graph

section_content: Generating a spectrogram involves several steps:

Step 1: Pre-emphasis is applied to the input audio signals. As the main analysis object of a spectrogram is human speech, a first-order high-pass filter is typically used to improve the SNR of the high-frequency part of the signal. This is related to the characteristics of speech itself and the weighting method used in the sound sampling process. In the fault diagnosis process, this step can be omitted to reflect the real-state information of the equipment. 

Step 2: The pre-emphasized signals are divided into a series of overlapping frames and windowing is applied to each data frame. The frame length affects the temporal resolution of the spectrogram. Frame overlap can smooth the transition between frames and avoid signal leakage caused by window boundaries; it is usually half the frame length. Windowing can reduce the Gibbs phenomenon caused by frame division and is beneficial for accurately extracting the frequency characteristics of the signals in the subsequent Fourier transform. Commonly used window functions include the rectangular, Hamming, and Hann windows. In this study, the Hann window is used, and the window function is represented as follows:

where i is the index of the data points within the frame and N is the length of the frame. 

Step 3: The STFT is performed on each segment to observe the frequency domain of the signals. The frequency-domain information X(m, k) of the m-th time window at frequency band k is expressed as follows:

where x(n) is the original signal, 0 ≤ m ≤ M -1, M is the number of frames, and 0

Step 4: The square of the amplitude of the STFT value is obtained and logarithmically processed. The purpose of logarithmic processing is to amplify low-amplitude components to observe signals masked by low-amplitude noise, the unit of which is dB. The processed signal Y(m, k) is expressed as follows:

Step 5: The processed signal Y(m, k) of each frame is concatenated to obtain the N 2 × M spectrogram matrix SP:

The matrix is plotted as an image, which can be presented in grayscale or pseudocolor to indicate the magnitude of the power. In this article, grayscale images are used. The elements in the matrix SP are mapped to the range of 0-255, and a grayscale image is drawn, where 0 represents black and 255 represents white. 

A traditional spectrogram is designed for human speech signals and used for speaker recognition and speech recognition. As the object of this study is mechanical equipment, it is necessary to improve the traditional spectrogram to satisfy the requirements of fault diagnosis.
section_name: Generating Spectrogram

section_content: Human voice signals have a higher energy in the low-frequency region; therefore, A-weighting is typically used to enhance the high-frequency region during the sampling process. The diagnostic object of this study is a hydraulic plunger pump that rotates periodically during operation. To obtain more realistic sound signals, a sound-level meter with Z-weighting is selected for the test process, which ensures that the sampled signals decay as little as possible over a wide frequency range. Therefore, when generating spectrograms for equipment fault diagnosis, the pre-emphasis processing of the sound signals can be omitted. 

Steps 2-5 of the traditional spectrogram generation are retained to generate the SP matrix, which can be used to draw a mean spectrogram bar graph. To comprehensively monitor the health status of the plunger pump, the sample used to draw a spectrogram should contain at least one rotation period, and the sample length l s should satisfy the following requirements:

where f s is the sampling frequency and n is the rotational frequency. 

Based on the analysis of the sound signals of the plunger pump, its changes in the time domain are relatively small, while the frequency domain can better reflect the fault information. Therefore, after obtaining the matrix SP, the mean value of each row element in the matrix is calculated to obtain the F order vector SSP:

The elements in the matrix SSP are mapped to a range of 0-255, and a grayscale image is drawn. This processing not only reduces the amount of data input into the deep learning model and speeds up the model calculation but also helps to eliminate occasional interference when l s is long. Figure 3 presents the grayscale spectrogram and mean spectrogram bar graphs drawn using the sound data of the plunger pump. The mean spectrogram bar graph is stretched horizontally for better observation. 

As shown in Figure 3, owing to the second mapping, the improved spectrogram has a stronger contrast, which can effectively improve the recognizability of the fault characteristics, reduce the learning difficulty, and improve the accuracy of the fault diagnosis.
section_name: Mean Spectrogram Bar Graph

section_content: Deep learning, which originates from artificial neural networks, is a machine learning technique that uses a multilayer perceptron with multiple hidden layers to automatically extract abstract features from massive amounts of data, thereby establishing direct relationships between these features and the target outputs. However, deep learning often requires a large number of samples to achieve optimal results. In contrast, transfer learning allows the use of models trained for specific tasks for solving related but different problems. This article proposes a deep transfer learning method that combines deep learning with model-based transfer learning. By integrating these two methods, deep transfer can achieve stronger feature extraction capabilities and end-to-end diagnostic functions with better universality. Moreover, deep transfer learning can be used to efficiently extract features from the data, thereby addressing the problem of insufficient data for fault diagnosis. As shown in Figure 3, owing to the second mapping, the improved spectrogram has a stronger contrast, which can effectively improve the recognizability of the fault characteristics, reduce the learning difficulty, and improve the accuracy of the fault diagnosis.
section_name: Deep Transfer Learning Algorithm Based on ResNet-50

section_content: Deep learning, which originates from artificial neural networks, is a machine learning technique that uses a multilayer perceptron with multiple hidden layers to automatically extract abstract features from massive amounts of data, thereby establishing direct relationships between these features and the target outputs. However, deep learning often requires a large number of samples to achieve optimal results. In contrast, transfer learning allows the use of models trained for specific tasks for solving related but different problems. This article proposes a deep transfer learning method that combines deep learning with model-based transfer learning. By integrating these two methods, deep transfer can achieve stronger feature extraction capabilities and end-to-end diagnostic functions with better universality. Moreover, deep transfer learning can be used to efficiently extract features from the data, thereby addressing the problem of insufficient data for fault diagnosis.
section_name: Deep Transfer Learning Algorithm Based on ResNet-50

section_content: ResNet-50 achieves interlayer connections through shortcuts, by adding the input to the output after convolution across layers, which fully trains the deeper-level network and improves the accuracy significantly with increasing depth, thus solving the model degradation problem. This structure is called a residual block. In this study, the convolutional residual block (CRB) and identity residual block (IRB) are used, as shown in Figure 4.
section_name: ResNet-50 Deep Learning Model

section_content: ResNet-50 achieves interlayer connections through shortcuts, by adding the input to the output after convolution across layers, which fully trains the deeper-level network and improves the accuracy significantly with increasing depth, thus solving the model degradation problem. This structure is called a residual block. In this study, the convolutional residual block (CRB) and identity residual block (IRB) are used, as shown in Figure 4. The ResNet-50 model used in this study contains 49 convolution layers and 1 fully connected layer, and the activation function used is the rectified linear unit (ReLU). The model structure is presented in Figure 5. The ResNet-50 model used in this study contains 49 convolution layers and 1 fully connected layer, and the activation function used is the rectified linear unit (ReLU). The model structure is presented in Figure 5. The ResNet-50 model used in this study contains 49 convolution layers and 1 fully connected layer, and the activation function used is the rectified linear unit (ReLU). The model structure is presented in Figure 5. Let the input sample be x, and by computing the convolution of the convolution kernel inp k and x, the feature vector inp f can be obtained:
section_name: ResNet-50 Deep Learning Model

section_content: where r  represents the ReLU activation function and inp k and inp b are trainable pa- rameters. Let the input sample be x, and by computing the convolution of the convolution kernel k inp and x, the feature vector f inp can be obtained:

where σ r represents the ReLU activation function and k inp and b inp are trainable parameters. Batch normalization is performed on f inp to improve the model training speed. Let f inp be an m × l matrix and f ij be the element in the i-th row and j-th column. The element n inp ij in the normalized result n inp is represented as follows:

where e is a small number added to prevent the mean squared error from becoming equal to zero. n inp ij is the element in the i-th row and j-th column. γ and β are trainable parameters. At this point, the distribution of f inp is adjusted to a standard normal distribution, which allows the input data to fall into the sensitive region of the activation function and prevents the occurrence of the vanishing gradients problem. However, this can result in a decrease in the expressive power of the network, and the depth of the network loses its meaning. Therefore, it is necessary to perform an inverse operation on the transformed n inp to enable the model to learn the optimal values of γ and β. 

The normalized result is split into several non-overlapping one-dimensional segments, the maximum value of each segment element is returned, and maximum pooling is performed. This can reduce the feature dimensions and number of trainable parameters. The output of the max pooling layer is denoted as P inp, which is a q × l matrix, where p inp k,j is the element in the k-th row and j-th column of P inp :

where s represents the length of the non-overlapping segments. Subsequently, the abstract features of P inp are extracted using multiple residual units. The residual unit calculates the sum of the residual function and input feature, where the residual function is a non-linear mapping F r (X) consisting of three convolution layers:

where

3 denotes the trainable parameter set for the t-th residual unit. In this study, the number of residual units is 16 when t = 1 and x 0 = P inp. 

In IRB, the shortcut connection performs element-wise addition between two features of equal dimensions, which is an identity mapping. The output Y t of the t-th residual unit can be expressed as follows:

where Y t denotes a multidimensional matrix. The IRB does not introduce additional parameters or computational complexity into the model, which has significant advantages in practical applications. 

In the CRB, as the output channels of the convolution layers are modified, the dimensions of the features are not equal when they are added; therefore, dimension matching is required to be performed in the shortcut connection (implemented through a 1 × 1 convolution). This increases the network parameters and improves its performance. This process can be represented as follows:

where K t 4 and B t 4 are parameters to be trained. After extracting the device fault feature matrix Y 16 through a 16-layer residual network, Global Average Pooling (GAP) [26] is performed, followed by a fully connected layer. The GAP calculates the average value of each element in the feature matrix of every dimension of Y 16, thus obtaining a feature vector g out with a length equal to the dimension of Y 16. Compared with directly connecting to a fully connected layer, this approach reduces the network parameters and prevents overfitting. 

Then, a fully connected layer is connected after the GAP layer, the features are mapped to the label space of the samples, and fc out is output as follows:

where θ fc = {w fc, b fc denotes the parameter set of the fully connected layer to be trained. 

Finally, the softmax function is used to calculate the probability distribution of fc out in the label space, and the probability of belonging to the s-th fault type is given as:

where θ soft = {w soft, b soft denotes the set of parameters to be trained for the softmax layer. In this study, the total number of types S = 8. Thus, the device fault diagnosis task based on deep learning is completed.
section_name: 7×7 convolution C

section_content: There are two basic concepts in transfer learning: Domains and tasks. The existing knowledge is called the source domain, which includes data knowledge and model knowledge. The source domain generates training samples containing rich annotation information. The new knowledge expected to be obtained is called the target domain, which generates test samples, usually without annotations or with only a small amount of annotation information. A task refers to a specific problem that is required to be solved, and a model is required to learn from a large amount of data. Transfer learning can transfer the annotated data and knowledge structure of the source domain to the target domain to complete the task of the target domain. Typically, the tasks of the source and target domains are identical; however, the data distributions are different. Pan et al. referred to this type of transfer learning as transductive transfer learning [27]. 

In model-based transfer learning, it is assumed that the source domain has a large amount of labeled data D s = x s i, y s i ns i=1, where x s i represents the i-th sample in the source domain, and y s i is its label. ns is the number of samples in the source domain. It is also assumed that the target domain has a small amount of labeled data D t = x t j, y t j nt j=1 wherein x t j represents the j-th sample in the target domain, y t j is its label, and nt is the number of samples in the target domain, where nt ns. In model transfer, an attempt is made to transfer the abundant structural knowledge learned from the source domain data to the target domain and obtain a more accurate prediction model using only a small amount of labeled data from the target domain. 

When transferring the model in this study, all the layer parameters, except for the fully connected layer, and softmax are frozen, which retains the knowledge of deep learning for the feature extraction of the device states. Transfer learning training is performed using the labeled data in the target domain. Parameter training comprises the use of a stochastic gradient descent with a momentum optimizer, which adds momentum to the stochastic gradient descent, and its iterative update formula for the parameters is as follows:

where m t represents the momentum, β is the decay coefficient, typically, β = 0. 9 [28], g t is the gradient of the objective function with respect to the current parameters, θ t is the optimized parameter, and α is the initial learning rate, which is obtained from optimization algorithms. 

The cost function J MSE uses the mean squared error, which is expressed as:

where y t j is the predicted label for sample x t j in the j-th sample of the target domain. L2 regularization is applied to the cost function to prevent overfitting.
section_name: Model Transfer Learning Based on ResNet-50

section_name: Test and Calculation

section_content: In the hydraulic plunger pump fault simulation test, the normal components of the pump are replaced with artificial faulty components to simulate seven types of faults: Shoe wear, port plate wear, loose shoe, plunger wear, and support bearing faults (outer ring, inner ring, and rolling element faults). The test also distinguishes the degree of some faults, which is, however, not introduced in this study. The faulty components used in the tests are presented in Figure 6. 

In the hydraulic plunger pump fault simulation test, the normal components of the pump are replaced with artificial faulty components to simulate seven types of faults: Shoe wear, port plate wear, loose shoe, plunger wear, and support bearing faults (outer ring, inner ring, and rolling element faults). The test also distinguishes the degree of some faults, which is, however, not introduced in this study. The faulty components used in the tests are presented in Figure 6. To produce worn shoe, plunger, and port plate components, we used 80-grit sandpaper to polish the surface of the components, and the weight of the former two decreased by 0. 4 g while the latter decreased by 1 g. To manufacture the loose-shoe component, we employed a tool to stretch the shoe and plunger, resulting in a 0. 3 mm increase in the distance between the two. The three types of fault-bearing components are obtained using electrical discharge machining. For the faulty inner and outer rings of the bearings, a through groove with a width of 1 mm and a depth of 1 mm is machined on the rolling track. For the faulty rolling elements of the bearings, a small pit with a diameter of 1 mm and a depth of 1 mm is machined on a rolling element. 

The fault simulation test bench of the hydraulic plunger pump is designed based on a constant-pressure variable-displacement axial plunger pump (see (2) in Figure 7 ), which uses a proportional relief valve (see (2) in Figure 4 ) to control the pressure of the test system and adjusts the opening size of the throttle valve (see (2) in Figure 7 ) to change the load pressure and flow rate. A hydraulic schematic of the test bench is presented in Figure 7. 

During the test, various types of fault signals are generated by replacing normal components with faulty components to simulate faulty pumps. In addition to the flow, pressure, and temperature sensors indicated in Figure 7, a sound level meter is also placed 0. 5 m directly above the pump to collect sound signals during the pump operation. In the test, the vibration signals from the pump are also sampled; however, this is not discussed in detail in this article. A photograph of the test bench is presented in Figure 8, and the main components and their performance indicators are listed in Table 1. 

track. For the faulty rolling elements of the bearings, a small pit with a diameter of 1 mm and a depth of 1 mm is machined on a rolling element. 

The fault simulation test bench of the hydraulic plunger pump is designed based on a constant-pressure variable-displacement axial plunger pump (see (2) in Figure 7 ), which uses a proportional relief valve (see (2) in Figure 4 ) to control the pressure of the test system and adjusts the opening size of the throttle valve (see (2) in Figure 7 ) to change the load pressure and flow rate. A hydraulic schematic of the test bench is presented in Figure 7. During the test, various types of fault signals are generated by replacing normal components with faulty components to simulate faulty pumps. In addition to the flow, pressure, and temperature sensors indicated in Figure 7, a sound level meter is also placed 0. 5 m directly above the pump to collect sound signals during the pump operation. In the test, the vibration signals from the pump are also sampled; however, this is not discussed in detail in this article. A photograph of the test bench is presented in Figure 8, and the main components and their performance indicators are listed in Table 1. The signals obtained using the various sensors are transmitted to the data acquisition card, and the data are displayed and stored in real time using the developed LabVIEW program. The front panel of the LabVIEW program is presented in Figure 9. During the test, the oil temperature is maintained at 35-40 • C by controlling the air-cooled unit to reduce the influence of changes in oil viscosity on the test. The tests are conducted at different pressures (5, 10, and 15 MPa) and flow rates (3, 6, and 9 L/min). Three different sampling frequencies (20, 30, and 40 kHz) are used to continuously sample the pump 10 times under each state, with each sampling duration being 5 s. Six pumps with the same model are individually installed on the bench, and seven faulty components are replaced to collect eight types of state data for the six pumps.
section_name: Fault Simulation Test of the Hydraulic Plunger Pump

section_content: The data-collection conditions for the sound signal dataset used in this study are as follows: A sampling frequency of 40 kHz, system pressure of 10 MPa, and flow rate of 9 L/min. The dataset includes data from six pumps, each with eight state signals (normal N, wear of shoe WS, wear of port plate WPP, loose shoe LS, wear of plunger WP, fault of bearing outer ring FOR, fault of bearing inner ring FIR, and fault of bearing rolling element FRE). Figure 10 presents the corresponding time-domain and frequency-domain graphs of the pump under these eight states. The bearings in the hydraulic pump are deep groove ball bearings, with model number 6205. The pitch diameter is 39. 04 mm, and there are 9 rolling elements with a diameter of 7. 94 mm. The rated speed of the driving motor is 1440 r/min, and the actual measured speed is 1488 r/min, which is equivalent to 24. 8 Hz. The characteristic frequency coefficients for the three types of bearing faults are 5. 4152 (FIR), 3. 5848 (FOR), and 2. 3567 (FRE). The corresponding characteristic frequencies are 134. 3 Hz, 88. 9 Hz, and 58. 4 Hz. The hydraulic pump has a total of seven plungers, which results in a pressure shock at 173. 6 Hz. We enlarged the frequency domain plot to display the range of 0-700 Hz, as shown in Figure 11. The bearings in the hydraulic pump are deep groove ball bearings, with model number 6205. The pitch mm, and there are 9 rolling elements with a diameter of 7. 94 mm. The rated speed of the driving motor is 1440 r/min, and the actual measured speed is 1488 r/min, which is equivalent to 24. 8 Hz. The characteristic frequency coefficients for the three types of bearing faults are 5. 4152 (FIR), 3. 5848 (FOR), and 2. 3567 (FRE). The corresponding characteristic frequencies are 134. 3 Hz, 88. 9 Hz, and 58. 4 Hz. The hydraulic pump has a total of seven plungers, which results in a pressure shock at 173. 6 Hz. We enlarged the frequency domain plot to display the range of 0-700 Hz, as shown in Figure 11. In Figure 11, peak values at the reference frequency of the pressure shock and its harmonics can be clearly observed. The fifth harmonic of the fault frequency of three types of bearing faults can be observed. 

The most commonly used method for measuring sound signals is the sound pressure level (dB). The measurement for the signals obtained in this test is that of sound pressure in Pa. The conversion formula between the two is:

where SPL represents the sound pressure level, P is the sound pressure, and ref P is the reference sound pressure. As the minimum sound pressure perceived by the human ear is 20 µPa, the reference sound pressure in air is generally considered 20 µPa. As shown in Figure 10, the sound pressure amplitude of the time domain in the normal state is less than that in the other states, and there are slight differences in the frequency domain graphs corresponding to different states; however, it is difficult to recognize their corresponding state categories. In transfer learning, the dataset is divided into four parts: Training, validation, support, and query sets. The training and validation sets are obtained from the source domain, In Figure 11, peak values at the reference frequency of the pressure shock and its harmonics can be clearly observed. The fifth harmonic of the fault frequency of three of faults can be observed. 

The most commonly used method for measuring sound signals is the sound pressure level (dB). The measurement for the signals obtained in this test is that of sound pressure in Pa. The conversion formula between the two is:

where SPL represents the sound pressure level, P is the sound pressure, and P re f is the reference sound pressure. As the minimum sound pressure perceived by the human ear is 20 µPa, the reference sound pressure in air is generally considered 20 µPa. As shown in Figure 10, the sound pressure amplitude of the time domain in the normal state is less than that in the other states, and there are slight differences in the frequency domain graphs corresponding to states; it is difficult to recognize their corresponding state categories. In transfer learning, the dataset is divided into four parts: Training, validation, support, and query sets. The training and validation sets are obtained from the source domain, and the samples are labeled; the support set comprises the training samples in the target domain, and the samples are labeled; and the query set comprises the test samples in the target domain, and the samples are unlabeled. Each data file collected through the test contains 5 s of sound signals, with 200 k data points, and each file is divided into 47 samples through overlapping partitioning. The dataset divisions used in this study are listed in Table 2. The data obtained from Pump_1-Pump_4 are used as the source domain data, while the data obtained from Pump_5 and Pump_6 are used as the target domain data. Considering that in practical engineering applications, normal data are usually abundant while failure data are scarce, all the normal samples in the target domain are selected as the support set.
section_name: Test Dataset

section_name: Fault Diagnosis Calculation Process

section_content: Before denoising, the original data are zero-centered. CEEMDAN is then used to decompose the sound signal into multiple modal components and selected components for reconstruction using the mRMR algorithm to complete the denoising process. For example, we demonstrated the denoising process of a sound signal in the normal state of Pump_1. The other signals are processed similarly. Figure 12 presents a comparison of the modal components after the EEMD and CEEMDAN are performed. 

The EEMD produced 17 IMF components and 1 residual, while the CEEMDAN produced 13 IMF components and 1 residual. Compared with EEMD, there are fewer signal components after the CEEMDAN, which indicates that CEEMDAN can effectively avoid the interference of white noise that is added multiple times by the EEMD to the original signal. 

In the mRMR algorithm, the mutual information I( c j ; x) between each component c j (s) and the original signal x(t) is first calculated:

where p[ c j (s), x(t)] represents the probability of their joint occurrence; p[ c j (s)] and p[x(t)] represent their individual occurrence probabilities; and S and T are equal, both of which are the number of signal points. The higher the mutual information, the higher the degree of dependence between the two random variables. The features are then sorted according to their mutual information values. Next, the first feature is selected as the most relevant feature from a sorted list of features. Subsequently, the mutual information of the remaining features with the selected feature is calculated, and their correlation with the selected feature is computed. The scores are then totaled to obtain a comprehensive score mRMR j :

where k denotes the number of selected features. Based on this, the features are sorted and the feature with the highest score is selected as the next feature. This process is repeated until sufficient features are selected or all the features are selected [29]. After sorting all the IMFs using the mRMR algorithm, the correlation coefficients between the signals are reconstructed with different numbers of selected IMFs and the original signal. In this study, the Pearson correlation coefficient is used and is calculated as follows:

where cov(x, y) represents the covariance of x and y, and var(x) represents the variance of x. Subsequently, the calculated correlation coefficient is used to plot the curve, as shown in Figure 13. The EEMD produced 17 IMF components and 1 residual, while the CEEMDAN produced 13 IMF components and 1 residual. Compared with EEMD, there are fewer signal components after the CEEMDAN, which indicates that CEEMDAN can effectively avoid the interference of white noise that is added multiple times by the EEMD to the original signal. 

In the mRMR algorithm, the mutual information ( ; )

and the original signal x(t) is first calculated: Subsequently, the calculated correlation coefficient is used to plot the curve, as shown in Figure 13. In Figure 13, the horizontal axis represents the ordering number of the IMFs obtained using the mRMR algorithm. The closer to the left end, the higher the mRMR comprehensive score of the component, with 14 indicating a residual. To elaborate, the score of IMF7 is the highest, with a correlation coefficient of 0. 74 with the original signal; the score of IMF2 ranks second, and the signal reconstructed using IMF2 and IMF7 has a correlation coefficient of 0. 78 with the original signal. Thus, the signal reconstructed using all the components and residual has a correlation coefficient of 1 with the original signal. The slope of the curve is relatively steep in the first half of the curve, and the curve growth slows after the addition of IMF1, at the location of the red dot marker. Therefore, the first eight IMFs (7/2/8/6/4/5/9/1) are selected to reconstruct the signal. The frequency-domain graphs of the reconstructed and original signals are presented in Figure 14. In Figure 13, the horizontal axis represents the ordering number of the IMFs obtained using the mRMR algorithm. The closer to the left end, the higher the mRMR comprehensive score of the component, with 14 indicating a residual. To elaborate, the score of IMF7 is the highest, with a correlation coefficient of 0. 74 with the original signal; the score of IMF2 ranks second, and the signal reconstructed using IMF2 and IMF7 has a correlation coefficient of 0. 78 with the original signal. Thus, the signal reconstructed using all the components and residual has a correlation coefficient of 1 with the original signal. The slope of the curve is relatively steep in the first half of the curve, and the curve growth slows after the addition of IMF1, at the location of the red dot marker. Therefore, the first eight IMFs (7/2/8/6/4/5/9/1) are selected to reconstruct the signal. The frequency-domain graphs of the reconstructed and original signals are presented in Figure 14. From Figure 14, it can be observed that, compared with the original signal, the reconstructed signal retains the data in the 0-2000 Hz frequency band, and the two frequency bands with high amplitudes of 5000-9000 Hz have some a enuation, but their peaks are preserved. The other frequency bands are suppressed to lower levels.
section_name: CEEMDAN Noise Reduction

section_content: The reconstructed signal is linearly normalized, and the data are mapped to [0, 1], which is expressed as:

where max

x and min x are the maximum and minimum values of all the data, respectively. 

The normalized data are plo ed as a spectrogram. It is worth noting that converting a spectrogram generated by a false color to a black-and-white image will result in the loss of a considerable amount of information and produce the same file size as a grayscale spectrogram. In this study, grayscale images are directly generated using raw data, which can compress the file size, reduce the network computing power, and fully record the device state information. For example, a comparison of the three visualizations using the From Figure 14, it can be observed that, compared with the original signal, the reconstructed signal retains the data in the Hz frequency band, and the two frequency bands with high amplitudes of 5000-9000 attenuation, but are preserved. bands suppressed to lower levels.
section_name: Mean Spectrogram Bar Graphs Generation

section_content: The reconstructed signal is linearly normalized, and the data are mapped to [0, 1], which is expressed as:

where x max and x min are the maximum and minimum values of all the data, respectively. The normalized data are plotted as a spectrogram. It is worth noting that converting a spectrogram generated by a false color to a black-and-white image will result in the loss of a considerable amount of information and produce the same file size as a grayscale spectrogram. In this study, grayscale images are directly generated using raw data, which can compress the file size, reduce the network computing power, and fully record the device state information. For example, a comparison of the three visualizations using the normal signal of Pump_1 is presented in Figure 15. 

which is expressed as:

where max

x and min x are the maximum and minimum values of all the data, respectively. 

The normalized data are plo ed as a spectrogram. It is worth noting that converting a spectrogram generated by a false color to a black-and-white image will result in the loss of a considerable amount of information and produce the same file size as a grayscale spectrogram. In this study, grayscale images are directly generated using raw data, which can compress the file size, reduce the network computing power, and fully record the device state information. For example, a comparison of the three visualizations using the normal signal of Pump_1 is presented in Figure 15. In the field of deep learning, color images can provide more feature information than grayscale images and are therefore widely used. However, this requires a premise: The RGB three-color channels of a color image contain different feature information. The spectrogram drawn using pseudocolors uses a three-dimensional tensor to express a two- In the field of deep learning, color images can provide more feature information than grayscale images and are therefore widely used. However, this requires a premise: The RGB three-color channels of a color image contain different feature information. The spectrogram drawn using pseudocolors uses a three-dimensional tensor to express a twodimensional matrix, which means that the feature information contained in the pseudocolor and grayscale spectrograms is identical. Using color images as the input to the network will inevitably result in a multiplication of the computational load owing to the increase in the number of channels. Therefore, in this study, grayscale images are used instead of pseudocolor images as network inputs. Figure 16 presents the grayscale spectrograms corresponding to the data of Pump_1 in the eight different states considered. 

dimensional matrix, which means that the feature information contained in the pseudocolor and grayscale spectrograms is identical. Using color images as the input to the network will inevitably result in a multiplication of the computational load owing to the increase in the number of channels. Therefore, in this study, grayscale images are used instead of pseudocolor images as network inputs. Figure 16 From Figure 16, it can be observed that the grayscale spectrogram has more obvious horizontal stripes. This indicates that the frequency domain has a be er feature representation ability for the pump-state data. Therefore, in this study, grayscale spectrograms are compressed into mean frequency spectrum bar graphs as inputs to the network. Figure 17 presents the mean spectrogram bar graphs corresponding to the eight states of Pump_1. To be er observe the compressed images, they are stretched horizontally. From Figure 16, it can be observed that the grayscale spectrogram has more obvious horizontal stripes. This indicates that the frequency domain has a better feature representation ability for the pump-state data. Therefore, in this study, grayscale spectrograms are compressed into mean frequency spectrum bar graphs as inputs to the network. Figure 17 presents the mean spectrogram bar graphs corresponding to the eight states of Pump_1. To better observe the compressed images, they are stretched horizontally. From Figure 16, it can be observed that the grayscale spectrogram has more obvious horizontal stripes. This indicates that the frequency domain has a be er feature representation ability for the pump-state data. Therefore, in this study, grayscale spectrograms are compressed into mean frequency spectrum bar graphs as inputs to the network. Figure 17 presents the mean spectrogram bar graphs corresponding to the eight states of Pump_1. To be er observe the compressed images, they are stretched horizontally. From Figure 17, it can be observed that the mean spectrogram bar graphs corresponding to each state are significantly different from those of the other states. The mean spectrum bar graphs of the three types of bearing faults are relatively similar; however, some differences can still be observed. Compared with the grayscale spectrogram, in addition to compressing the file size, the mean spectrogram bar graph has a higher contrast and more obvious differences between the images of each state. These From Figure 17, it can be observed that the mean spectrogram bar graphs corresponding to each state are significantly different from those of the other states. The mean spectrum bar graphs of the three types of bearing faults are relatively similar; however, some differences can still be observed. Compared with the grayscale spectrogram, in addition to compressing the file size, the mean spectrogram bar graph has a higher contrast and more obvious differences between the images of each state. These significant differences can reduce the difficulty in learning the feature information and improve the accuracy of the diagnostic model.
section_name: Mean Spectrogram Bar Graphs Generation

section_content: Previously, the dataset used in many studies on fault-diagnosis equipment, such as bearings and hydraulic plunger pumps, is obtained using a single device. Therefore, the similarity between the datasets was extremely high and the diagnostic accuracy was often close to 100%. This study demonstrates this in the same manner by dividing the Pump _1 data into training, validation, and test sets in a 5:2:3 ratio; the corresponding sample sizes are divided in the ratio of 1880:752:1128. Via overlapping segmentation of the original data, 47 mean spectrogram bar graphs are generated for each data file. The mean spectrogram bar graphs generated from the training and validation sets are imported into the ResNet-50 model for training. The training process parameters are listed in Table 3. The numbers in the confusion matrix represent the number of samples, with the last row indicating the recall rate for each true class and the last column indicating the precision rate for each predicted class. The bo om-right corner presents the testing accuracy of the proposed model. It can be observed that the diagnostic model's accuracy reached 100% when tested with samples obtained from Pump_1, but the accuracy decreased to an unacceptable level when tested with other pump samples. The testing accuracy of the model for the samples obtained from Pump_2-Pump_6 is 22. 7%, 13. 4%, 12. 8%, 12. 3% and 13. 1%, respectively. The confusion matrices shown in Figure 20 illustrate the performance of the model when tested using samples from different pumps. The numbers in the confusion matrix represent the number of samples, with the last row indicating the recall rate for each true class and the last column indicating the precision rate for each predicted class. The bottom-right corner presents the testing accuracy of the proposed model. It can be observed that the diagnostic model's accuracy reached 100% when tested with samples obtained from Pump_1, but the accuracy decreased to an unacceptable level when tested with other pump samples. The testing accuracy of the model for the samples obtained from Pump_2-Pump_6 is 22. 7%, 13. 4%, 12. 8%, 12. 3% and 13. 1%, respectively. The confusion matrices shown in Figure 20 illustrate the performance of the model when tested using samples from different pumps. For practical applications, a well-trained diagnostic model should be able to diagnose the majority of pumps in the same model and have a satisfactory diagnostic effect. Therefore, it is necessary to conduct research on transfer learning to increase the generalization of the diagnostic model. 

The ResNet-50 diagnostic model is trained on a dataset that is divided as shown in Table 2. The parameters used in the first training process are the same as those in Table 3 For practical applications, a well-trained diagnostic model should be able to diagnose the majority of pumps in the same model and have a satisfactory diagnostic effect. Therefore, it is necessary to conduct research on transfer learning to increase the generalization of the diagnostic model. 

The ResNet-50 diagnostic model is trained on a dataset that is divided as shown in Table 2. The parameters used in the first training process are the same as those in Table 3, but the maximum number of epochs is adjusted to 40 to increase the training time of the model. The model training process curves are presented in Figure 21. 

Although the training accuracy quickly reached 100%, it can be observed that the validation accuracy can only reach a maximum of 65. 1%, at the location of the red dot marker and the minimum validation loss is 1. 32, at the location of the red dot marker. Because the selected dataset for the validation set comprises data on the seven types of faults of Pump_4, which is quite different from the data of Pump_1-Pump_3 used for training. Therefore, Bayesian optimization is used to optimize the initial learning rate of the model between [0. 00001, 0. 1], a logarithmic scale is used for the horizontal axis, and 30 attempts are made. The loss and accuracy corresponding to different initial learning rates during the optimization process are presented in Figure 22. 

the majority of pumps in the same model and have a satisfactory diagnostic effect. Therefore, it is necessary to conduct research on transfer learning to increase the generalization of the diagnostic model. 

The ResNet-50 diagnostic model is trained on a dataset that is divided as shown in Table 2. The parameters used in the first training process are the same as those in Table 3 Although the training accuracy quickly reached 100%, it can be observed that the validation accuracy can only reach a maximum of 65. 1%, at the location of the red dot marker and the minimum validation loss is 1. 32, at the location of the red dot marker. Because the selected dataset for the validation set comprises data on the seven types of faults of Pump_4, which is quite different from the data of Pump_1-Pump_3 used for training. Therefore, Bayesian optimization is used to optimize the initial learning rate of the model between [0. 00001,0. 1], a logarithmic scale is used for the horizontal axis, and 30 a empts are made. The loss and accuracy corresponding to different initial learning rates during the optimization process are presented in Figure 22. As shown in Figure 22, in the 30 a empts, the validation accuracy realized is the highest at 95. 2%, at the location of the red dot marker and the validation loss is 0. 13, at the location of the red dot marker when the initial learning rate is set as 0. 023. The model trained at this learning rate, which had the minimum validation loss, is exported as the base model for the transfer learning. The confusion matrices of the diagnostic results obtained for the data of Pump _5 and Pump _6 using the base model are presented in Figure 23. As shown in Figure 22, in the 30 attempts, the validation accuracy realized is the highest at 95. 2%, at the location of the red dot marker and the validation loss is 0. 13, at the location of the red dot marker when the initial learning rate is set as 0. 023. The model trained at this learning rate, which had the minimum validation loss, is exported as the base model for the transfer learning. The confusion matrices of the diagnostic results obtained for the data of Pump _5 and Pump _6 using the base model are presented in Figure 23. 

As can be observed from During the training process, the highest validation accuracy is 99. 5%, at the location of the red dot marker and the minimum validation loss is 0. 28, at the location of the red dot marker. The model with the minimum validation loss is selected as the final diagnostic model and is used to diagnose the Pump_5 and Pump_6 data. The confusion matrices for the diagnostic results are presented in Figure 25. 

location of the red dot marker when the initial learning rate is set as 0. 023. The model trained at this learning rate, which had the minimum validation loss, is exported as the base model for the transfer learning. The confusion matrices of the diagnostic results obtained for the data of Pump _5 and Pump _6 using the base model are presented in Figure 23. During the training process, the highest validation accuracy is 99. 5%, at the location of the red dot marker and the minimum validation loss is 0. 28, at the location of the red dot marker. The model with the minimum validation loss is selected as the final diagnostic model and is used to diagnose the Pump_5 and Pump_6 data. The confusion matrices for the diagnostic results are presented in Figure 25. It can be observed that the transferred model can also provide better diagnosis results for the new data. The diagnostic accuracies of Pump_5 and Pump_6 are 86. 5% and 90. 8%, respectively. The diagnostic results are 38. 4% and 36. 8% higher than those of the baseline model, respectively. The diagnostic accuracies of the different models used in this study are listed in Table 4. It can be observed that the transferred model can also provide be er diagnosis results for the new data. The diagnostic accuracies of Pump_5 and Pump_6 are 86. 5% and 90. 8%, respectively. The diagnostic results are 38. 4% and 36. 8% higher than those of the baseline model, respectively. The diagnostic accuracies of the different models used in this study are listed in Table 4. In this paper, a transfer learning study based on SqueezeNet and GoogleLeNet is also conducted. Using the same dataset mentioned previously, the accuracy of the fault diagnosis with different transfer models is presented in Table 5. From Table 5, it can be observed that the hydraulic-pump fault-diagnosis model based on the transfer of ResNet-50 has higher diagnostic accuracy than the other two models. 

In addition, to observe whether the model in this study exhibits the overfi ing phenomenon, models trained with different numbers of epochs are used to diagnose the data obtained from Pump_6, and the result is presented in Figure 26. In this paper, a transfer learning study based on SqueezeNet and GoogleLeNet is also conducted. Using the same dataset mentioned previously, the accuracy of the fault diagnosis with different transfer models is presented in Table 5. From Table 5, it can be observed that the hydraulic-pump fault-diagnosis model based on the transfer of ResNet-50 has higher diagnostic accuracy than the other two models. 

In addition, to observe whether the model in this study exhibits the overfitting phenomenon, models trained with different numbers of epochs are used to diagnose the data obtained from Pump_6, and the result is presented in Figure 26. From Figure 26, it can be observed that, with an increase in the number of model epochs, the test accuracy realized using data obtained from Pump_6 increases significantly when the number of training epochs is increased from one to five. Subsequently, the growth becomes slower, and two plateau periods are observed: [4, 7] and [9, 26]. Finally, when the number of training epochs reaches 39, the test accuracy exceeds 90%, at the location of the red dot marker and remains unchanged thereafter. It can be considered that the model did not exhibit overfi ing owing to the increase in the number of training epochs. 

In addition, there is one point that needs to be recognized. The sampled signals contain environmental noise, the sound generated by various components of the hydraulic system, and specific frequency sounds caused by faults. Environmental noise From Figure 26, it can be observed that, with an increase in the number of model epochs, the test accuracy realized using data obtained from Pump_6 increases significantly when the number of training epochs is increased from one to five. Subsequently, the growth becomes slower, and two plateau periods are observed: [4, 7] and [9, 26]. Finally, when the number of training epochs reaches 39, the test accuracy exceeds 90%, at the location of the red dot marker and remains unchanged thereafter. It can be considered that the model did not exhibit overfitting owing to the increase in the number of training epochs. 

In addition, there is one point that needs to be recognized. The sampled signals contain environmental noise, the sound generated by various components of the hydraulic system, and specific frequency sounds caused by faults. Environmental noise can be effectively suppressed through CEEMDAN and mRMR. However, the sound generated by the operation of the hydraulic system cannot be eliminated by noise-reduction algorithms, and the sound signals will also change when the working pressure or spatial arrangement changes. It can be expected that this will have an adverse effect on the accuracy of fault diagnosis. Fortunately, this problem can be solved. By using re-sampled sound signals to conduct transfer learning, the fault diagnosis algorithm can be adapted to the changed hydraulic system.
section_name: Use of Deep Transfer Learning for Fault Location

section_content: This article first elaborates on the CEEMDAN and mRMR denoising algorithms, the method of drawing mean spectrogram bar graphs based on sound signals, and the transfer learning algorithm based on the ResNet-50 model used in the research process. The test settings and sampled data are then described, and the data calculation process and diagnostic results are displayed. The following conclusions are thus drawn:

1. 

Compared with the IMF components obtained via EEMD, the CEEMDAN can provide fewer IMF components and can effectively denoise the signal by combining it with the mRMR algorithm.
section_name: Conclusions

section_content: The mean spectrogram bar graph significantly compresses the input diagnostic model images while ensuring the diagnostic accuracy of the diagnostic model.
section_name: 2.

section_content: A diagnostic model trained using samples from a single pump cannot be used to diagnose the faults in other pumps. In this study, the diagnostic accuracy of the diagnostic model trained using data from Pump_1 for diagnosing other pumps is found to be less than 22. 7%.
section_name: 3.

section_content: The baseline diagnostic model trained directly with samples from multiple pumps has a higher diagnostic accuracy when applied with other pumps than the diagnostic model trained with samples obtained from a single pump. In this study, the baseline diagnostic model trained using samples from Pump_1-Pump_3 achieved diagnostic accuracies of 48. 1% and 58% for Pump_5 and Pump_6, respectively, which reflect increases of 35. 8% and 44. 9%, respectively, compared with the diagnostic accuracy of the model trained with samples obtained from a single pump. 5. 

After the transfer learning, the diagnostic model has a higher diagnostic accuracy than the baseline diagnostic model trained using samples from multiple pumps. The former model can achieve a diagnostic accuracy of 86. 1% and 90. 8% for Pump_5 and Pump_6, respectively, which is an increase of 38. 4% and 36. 8%, respectively, compared with the diagnostic accuracy of the latter model. 6. 

Compared with SqueezeNet and GoogleLeNet, the hydraulic-pump fault-diagnosis model based on the transfer of ResNet-50 has higher diagnostic accuracy. The diagnostic accuracies obtained for the Pump_5 validation data are 67. 7%, 77. 5%, and 86. 1% for the three models, respectively. The accuracies obtained for the Pump_6 data are 68. 4%, 84. 1%, and 90. 8%, respectively. 

Research on fault diagnosis based on sound signals has expanded the applications of fault diagnosis methods. In the future, signals such as sound, vibration, and pressure can be fused and used for the fault diagnosis of hydraulic equipment. As inputs to the model, the three channels of the image can be used separately to store different input information. The mean spectrogram bar graph presented in this study can significantly compress the input image file size required for fault diagnosis and can be promoted in appropriate fields. Transfer learning can be used to effectively solve the problem of insufficient data in fault diagnosis and other research fields and has great research value. Ensemble learning, which combines multiple deep-learning models, has been proven to have better generality than a single model [30]. Each model can be specialized for processing specific data subsets during
section_name: 4.

section_content: the training phase [31], which can aid in fully leveraging the performances of different models. Therefore, future research should be focused on this direction.

section_name: Institutional Review Board Statement: Not applicable.

section_content: Funding: This research was supported by the National Natural Science Foundation of China, China (Grant Nos. 52275067, 51875498, and 52275069 ) and the Natural Science Foundation of Hebei Province, China (Grant No. F2020203058 ).

section_content: Informed Consent Statement: Not applicable.
section_name: Institutional Review Board Statement: Not applicable.

section_name: Data Availability Statement:

section_content: No data sharing available.
section_name: Data Availability Statement:

section_content: The authors declare no conflict of interest.
section_name: Conflicts of Interest:

section_content: Author Contributions: Methodology, P. Z. ; Software, P. Z. ; Validation, P. Z. ; Formal analysis, P. Z. ; Investigation, S. Z. (Sheng Zhang); Resources, W. J. and S. Z. (Sheng Zhang); Data curation, P. Z. ; Writing-original draft, P. Z. and Y. Z. ; Writing-review & editing, W. J. ; Visualization, P. Z. ; Supervision, W. J. ; Funding acquisition, W. J., S. Z. (Shuqing Zhang) and S. L. All authors have read and agreed to the published version of the manuscript.