Random Forest Model of Flow Pattern Identiﬁcation in Scavenge Pipe Based on EEMD and Hilbert Transform

: Complex oil and gas two-phase ﬂow exists within an aero-engines bearing cavity scavenge pipe, prone to lubricated self-ignition and coking. Lubricant system designers must be able to accurately identify and understand the ﬂow state of the scavenge pipe. The prediction accuracy of previous models is insufﬁcient to meet the more demanding needs. This paper establishes a visualized ﬂow pattern identiﬁcation test system for the scavenge pipe, with a test temperature of up to 370 k, using a high-speed camera to photograph four ﬂow patterns, decomposing the pressure signals obtained from high-frequency dynamic pressure sensors using the ensemble empirical mode decomposition (EEMD) method, and then performing Hilbert transform, using the Hilbert spectrum to quantify the changes of amplitude and frequency with time, and establishing the energy and ﬂow pattern correspondence analysis. Then the energy percentage of IMFs is used as the input of feature values, and the random forest algorithm machine learning is used for predictive classiﬁcation. The experimental results show that the ﬂow pattern recognition rate established in this paper can reach 98%, which can identify the two-phase ﬂow pattern in the scavenge pipe more objectively and accurately.


Introduction
In modern aero engines, lubrication systems play a crucial role [1].Slip oils are not only used to lubricate bearings, gears, or seals, but they also carry away a lot of heat and have a cooling effect on the aero engine.However, high temperatures can also cause spontaneous combustion and coking of the lubricating oil, leading to engine reliability and safety problems.Usually, there is a vent pipe at the top of the bearing cavity and scavenge pipe at the bottom, which is connected to the oil return pump; the vent pipe discharges most of the sealed air to avoid high pressure in the bearing cavity and the scavenge pipe discharges most of the lubricating oil and a small amount of air to make the lubricating oil circulate.The aero-engine bearing cavity scavenge pipe has not only a large amount of high-temperature lubricating oil but also a small amount of high-pressure seal gas, and there is a complex oil and gas two-phase flow, which is more likely to occur the spontaneous combustion and coking of slip oil.The bearing cavity scavenging the pipe's flow state is vital in correctly designing the bearing cavity lubrication and significantly impacts the flow's pressure loss and heat transfer characteristics [2,3].It can be seen that the study of the flow pattern identification method of the scavenge pipe is of great significance to the design of bearing cavity lubrication.
Traditional two-phase flow pattern discrimination mainly uses flow pattern diagrams and transformation criterion relational equations, which have shortcomings such as low recognition accuracy and more restrictive conditions for use.The rapid development of electronic computer technology, gas-liquid two-phase flow measurement technology, nonlinear signal analysis theory, modern information processing technology, artificial neural Energies 2023, 16, 6084 2 of 15 network technology [4,5], and data fusion technology provides new ideas and methods for the flow pattern identification problem.Ong C.L. et al. proposed a new macro-microscale flow pattern map for application to both the lower end of macroscale flows and the upper end of microscale flows [6].Rafalko G. et al. used images directly for two-phase flow pattern identification [7,8], Sizikov V. S. et al. proposed the restoration of nonuniformly smeared images [9], VT Nguyen et al. used wavelets for research [10][11][12].Ji H used the EMD method [13,14] for signal decomposition by obtaining the signal of two-phase flow pressure fluctuation in the tube for analysis.However, the wavelet decomposition can only achieve fixed resolution decomposition of the signal, which may lose some meaning of the signal itself.In contrast, EMD decomposition [15] has orthogonality and completeness and can decompose the signal adaptively, but there may be a modal aliasing phenomenon and endpoint effect.Azzedine Dliou proposed a combination of EMD and the Choi-Williams time-frequency techniques for analyzing a noised ECG signal [16].Huang proposed the EEMD method [17] for the shortcomings of the EMD method, which suppresses the modal aliasing sensation by adding normal and spectrally uniformly distributed Gaussian white noise [18].The Hilbert transform based on the EMD method is a new method with adaptive time-frequency analysis [19].It can obtain a high resolution with good time-frequency aggregation, suitable for studying non-smooth and nonlinear signals.V. G. Kurbatskii et al. used Hilbert-Huang transform for signal analysis and prediction [20,21].Therefore, this paper adopts EEMD decomposition and Hilbert transform for signal processing, which can improve the decomposition accuracy.
Machine learning, the core of artificial intelligence, uses cutting-edge algorithms to find and learn the implicit laws in the data by analyzing and learning from large amounts of systematically stored data to form predictive ideas [22,23].It has been applied to various fields that require signal analysis, and B Yadav et al. used neural networks for prediction and classification [24,25].However, neural network models often require extended training and detailed data pre-processing.Random forest [26] models perform better for data containing different features.It does not require iterative adjustment of parameters to give good results or data scaling.Therefore, improving data quality with EEMD and assisting random forest model machine learning can effectively enhance the performance of flow pattern recognition models.
In this paper, a visualized flow pattern identification test system [27,28] was established with a test temperature of 370 k in the scavenge pipe, a high-speed camera was used to photograph the flow pattern, and a high-frequency dynamic pressure sensor was used to obtain the pressure fluctuation signals at both ends of the test section.EEMD is used to decompose the signal, and Hilbert transforms, and the Hilbert spectrum is used to analyze the changes of amplitude and frequency with time quantitatively.Then the energy share of IMFs is input as the feature value, and the random forest algorithm machine learning is used to predict the classification.Finally, a random forest model based on EEMD decomposition and Hilbert transform is constructed to identify the oil and gas two-phase flow patterns in the scavenge pipe.

EEMD
The empirical mode decomposition (EMD) method is suitable for analyzing and processing non-stationary and nonlinear signals, but there are problems and shortcomings in this method, mainly: (i) the existence of modal aliasing phenomenon with EMD decomposition; (ii) the end effect affects the decomposition effect.To suppress the mode-mixing phenomenon of EMD, Huang proposed an improved EMD method based on noise-assisted analysis, namely ensemble empirical mode decomposition (EEMD), by using EMD to perform statistics on the results of white noise decomposition [14].
The EEMD method is a multiple empirical mode decomposition with superimposed Gaussian white noise, which uses the statistical property of Gaussian white noise with uniform frequency distribution to change the polar characteristics of the signal by adding Energies 2023, 16, 6084 3 of 15 different white noise of equal amplitude each time, and then the overall average of the corresponding IMF obtained by multiple EMD to cancel the added white noise, to suppress the generation of modal aliasing effectively.The decomposition steps are as follows: (1) Set the overall average number of times M; (2) Add a white noise with standard normal distribution to the original signal to produce a new signal: x m (t) = x(t) + n m (t), where n m (t) denotes the mth addition of white noise sequence, x m (t) denotes the signal with the ith additional noise, m = 1, 2. ..M; (3) EMD decomposition is performed separately for the resulting noise-containing signals, and the respective sums are obtained in the form: where, (5) Using the principle that the statistical mean of uncorrelated sequences is zero, the above corresponding IMFs are subjected to a pooled averaging operation to obtain the final IMF after EEMD decomposition: where I i (t) is the nth IMF of the EEMD decomposition, m = 1, 2, . . .M, i = 1, 2, . . .n, and the residual component is the difference between the original data and j IMFs.
The final decomposition form of the signal is obtained:

Hilbert Transform
The Hilbert transform is an essential mathematical tool widely used to convert an actual signal into a complex signal for easy frequency domain analysis and processing.Firstly, the EEMD method decomposes the signal into the sum of several IMF components.Then the Hilbert transform is applied to each IMF component to obtain the instantaneous frequency and instantaneous amplitude, which results in the Hilbert spectrum of the signal, and the Hilbert spectrum represents the complete time-frequency distribution of the signal [13].
Hilbert transform for each IMF component I i (t): Constructing parsed signals: Energies 2023, 16, 6084 where a i (t) is the instantaneous amplitude function, φ i (t) is the instantaneous phase function; The instantaneous frequency of each IMF component is: As can be seen from the above equation, the instantaneous amplitude and instantaneous frequency of the IMF component are functions of time, with time t and instantaneous frequency ω i (t) as the independent variables, the amplitude of the signal can be expressed as a function H(ω, t) of t and ω, and Re is taken as the real part.The time-frequency distribution of this amplitude is called the Hilbert amplitude spectrum, or the Hilbert spectrum for short.
Defining the Hilbert marginal spectrum: T is the total length of the signal.H(ω, t) describes precisely the variation pattern of the amplitude of the signal over the whole frequency band with time and frequency while h(ω) reflecting the variation of the amplitude of the signal over the entire frequency band with frequency.

Random Forest
Random Forest is an integrated learning method that improves classification accuracy by combining multiple decision trees.In the Random Forest algorithm, for each decision tree construction, the data used are randomly selected from a certain percentage of the original data (with put-back sampling).This method can effectively reduce the overfitting problem and handle high-dimensional and large amounts of data (Figure 1).

Experimental Apparatus
This paper selected a horizontal pipe with a diameter of 10 mm as the scavenge pipe, and the visualization line was 200 mm long for the return flow simulation test.The mixing temperature of the two phases was set at 350 k, the liquid phase commutation speed was guaranteed to be constant for each group of tests, and the inlet flow rate was increased sequentially.The training process of the random forest classifier is as follows: Energies 2023, 16, 6084 5 of 15 (1) Randomly select a portion of data and features from the original data to construct a decision tree.(2) Repeat n times to construct n decision trees.
(3) For each new data point, pass it into each tree for classification and get the classification result for each tree.(4) Each new data point is classified as the most classified result of the random forest.
The final classification result is shown in Equation: where H(x) is the final voting result of the combined classification model, h i (x) is a single decision tree model, Y denotes the target variable, and I(•) is an indicator function.

Experimental Apparatus
This paper selected a horizontal pipe with a diameter of 10 mm as the scavenge pipe, and the visualization line was 200 mm long for the return flow simulation test.The mixing temperature of the two phases was set at 350 k, the liquid phase commutation speed was guaranteed to be constant for each group of tests, and the inlet flow rate was increased sequentially.
Where J L and J G denote the converted velocities of the liquid and gas phases, respectively, defined as: where: V L , V G is the volume flow rate of the liquid and gas phases, respectively; A is the cross-sectional area of the flow channel.
Figure 2 shows the scavenge pipe's flow pattern identification test system, the twophase flow pattern images in the horizontal tube, and the pressure fluctuation signals at both ends of the test tube obtained under different flow conditions.Ensure that the lubricant entering the pipeline is the target temperature.The lubricating oil in the tank is heated first to increase the base temperature.Then the lubricating oil enters the lubricating oil heater for heating, which has temperature feedback to ensure the outlet temperature is the target temperature.The liquid flow controller adjusts the flow rate entering the mixer.The compressed air from the air compressor enters the pressure stabilization tank.Then it enters the air heater with temperature control to ensure the target outlet temperature.The gas flow controller adjusts the gas flow rate into the mixer.The gas-liquid mixing device mixes the two-phase flow and enters the visualization pipeline, as shown in Figure 3.The temperature and pressure measurement points are located at both ends of the visualization pipe, with the pressure measurement point at the top and the temperature measurement point at the bottom, and the visualization section is made of high-strength quartz glass so that a high-speed camera can capture the flow pattern of the scavenge pipe.The high-speed camera used in the resolution of 1280 × 1024 has a shooting frame rate of up to 9300 fps.This test selected a frame rate of 3000 fps, which can shoot the flow image; using PCB Piezotronics model M105C02, high-frequency dynamic pressure sensor pressure range of 0~690 kPa, the maximum measurement frequency is 256 kHz, non-linearity ≤ 2%FS.The mixed fluid enters the gas-liquid separator at the end of the test section, and the separated slip oil flows back into the tank to reach the system's circulation.The parameters of the two-phase medium used for the test are shown in Table 1.
pipe.The high-speed camera used in the resolution of 1280 × 1024 has a shooting frame rate of up to 9300 fps.This test selected a frame rate of 3000 fps, which can shoot the flow image; using PCB Piezotronics model M105C02, high-frequency dynamic pressure sensor pressure range of 0~690 kPa, the maximum measurement frequency is 256 kHz, nonlinearity ≤2%FS.The mixed fluid enters the gas-liquid separator at the end of the test section, and the separated slip oil flows back into the tank to reach the system's circulation.The parameters of the two-phase medium used for the test are shown in Table 1.This study selected a horizontal pipe with a diameter of 10 mm as the scavenge pipe, the visualization line length was 200 mm, and the flow simulation test of the scavenge pipe was carried out.Set the two-phase mixing temperature to 350 k.This group of analysis images ensures that the liquid phase commutation velocity is 0.85 m/s constant, sequentially increasing the inlet flow rate and the gas phase commutation velocity change range of 0~2.12 m/s.The resulting flow pattern change filming results are shown in Table 2. To better display the shooting effect, the images were processed by simple binarization, and the graphics and contours were obtained by contrasting the prints and selecting a suitable threshold value, which made the photos show a black-and-white effect.

Decomposing the Pressure Signal of Oil Gas Two-Phase Flow in the Scavenge Pipe
The oil-gas two-phase flow in the scavenge pipe was tested to obtain the pressure fluctuation signals at both ends of its tube, and the fluctuation of the pressure signals showed relatively unstable characteristics from which information could not be directly obtained.Therefore, EEMD decomposition and Hilbert transformation of the obtained high-frequency pressure signal curves are carried out, which helps to carry out a more indepth analysis.Four typical flow types, namely, bubble flow, plug flow, slug flow, and annular flow, are selected for research, and a set of pressure fluctuation signals are selected for each.The IMFs, Hilbert spectra, and marginal spectra of the pressure fluctuation signals of the four flow types are shown in Figures 4-7. (a)

(a)
Ensure the liquid phase flow rate is constant and gradually increase the inlet air volume.When the inlet gas is small, the bubble flow is bubble-like; due to the influence of gravity, the bubbles are located in the upper part of the tube, and when the gas phase conversion speed reaches 0.127 m/s, the tiny bubbles merge into gas plugs, forming plug flow.The flow structure is also located in the upper part of the tube flow.With the further increase of the gas phase conversion velocity to 0.254 m/s, the gas phase forms a pop-up structure in the tube, called pop-up flow, which will be further elongated with the further increase of the inlet volume.As the inlet volume continues to increase, the gas phase velocity is high and thus forms a circular flow pattern, the liquid phase is pressed against the wall to form a circle, and the gas phase flows in the center of the tube.Still, the liquid film on the lower tube wall is thicker than the upper tube wall due to the effect of gravity during horizontal flow.

Decomposing the Pressure Signal of Oil Gas Two-Phase Flow in the Scavenge Pipe
The oil-gas two-phase flow in the scavenge pipe was tested to obtain the pressure fluctuation signals at both ends of its tube, and the fluctuation of the pressure signals showed relatively unstable characteristics from which information could not be directly obtained.Therefore, EEMD decomposition and Hilbert transformation of the obtained high-frequency pressure signal curves are carried out, which helps to carry out a more in-depth analysis.Four typical flow types, namely, bubble flow, plug flow, slug flow, and annular flow, are selected for research, and a set of pressure fluctuation signals are selected for each.The IMFs, Hilbert spectra, and marginal spectra of the pressure fluctuation signals of the four flow types are shown in Figures 4-7.
Comparing the original pressure fluctuation signal s1 collected from the four flow types, the original signal of vesicular flow is relatively stable with a small signal fluctuation amplitude, and with the increase of gas phase flow, the signal fluctuation amplitude increases more obviously, and it is nonlinear and non-stationary.The EEMD method decomposes the pressure signal into IMFs with different scales and multiple resolution levels.The original pressure signal is plotted in the first row of the decomposition diagram.The EMD method decomposition and sieving process are direct and adaptive, decomposing the signal based on its information to obtain the IMFs components, and each IMF component represents the real physical information contained within the signal.Summing the individual IMF components and the residual term yields the original signal, proving the completeness of the EMD decomposition.The EEMD method decomposition results in IMFs that satisfy the characteristic conditions, and these IMFs exhibit modulation characteristics in frequency or amplitude, further indicating that the EEMD method is well suited to handle nonlinear and non-stationary signals.Comparing the original pressure fluctuation signal s1 collected from the four flow types, the original signal of vesicular flow is relatively stable with a small signal fluctuation amplitude, and with the increase of gas phase flow, the signal fluctuation amplitude increases more obviously, and it is nonlinear and non-stationary.The EEMD method decomposes the pressure signal into IMFs with different scales and multiple resolution levels.The original pressure signal is plotted in the first row of the decomposition diagram.The EMD method decomposition and sieving process are direct and adaptive, decomposing the signal based on its information to obtain the IMFs components, and each IMF component represents the real physical information contained within the signal.Summing the individual IMF components and the residual term yields

Establishing Flow Pattern Recognition Model
To quantitatively analyze the transformation law of pressure fluctuation signal with flow pattern, it is necessary to select the eigenvalues to establish the widespread recognition model.Based on the division of IMF components from high to low frequency, the energy of each IMF component of the signal in the frequency band as a percentage of the total energy as a characteristic value to reflect the information of the flow pattern, the typical energy value k is introduced and defined as: Energies where E i is the energy of each IMF component of the frequency band (i = 1, . .., 12), the total energy E = E 1 + E 2 + . . .+ E 12 , and accordingly, k i is the ratio of the energy of each IMF component of the frequency band to the total energy, respectively.The following equation calculates the energy of the ith IMF component Twelve feature values are obtained, and the random forest integrated learning method is used for classification and identification.Classification using machine learning is more objective.The random forest classification method has better accuracy, can run effectively on large data sets, has good noise immunity, but can overfit in the case of relatively noisy data, can handle very high dimensional data, and does not use dimensionality reduction.
For easy identification of the records, the individual flow types are coded, as shown in Table 3.The 150 samples will be tested for training; from Figure 8a, its total recognition rate can reach 98%, and the recognition effect is excellent.As shown in Figure 8b, seeing the recognition classification effect through the confusion matrix can be more intuitive.Longitudinally, it is found that the two groups are incorrectly recognized as bullet-like flow, and the feature values obtained at the boundary of flow type conversion have certain similarities, thus producing specific false recognition.For easy identification of the records, the individual flow types are coded, as sh in Table 3.The 150 samples will be tested for training; from Figure 8a, its total recogn rate can reach 98%, and the recognition effect is excellent.As shown in Figure 8b, s the recognition classification effect through the confusion matrix can be more intu Longitudinally, it is found that the two groups are incorrectly recognized as bulle flow, and the feature values obtained at the boundary of flow type conversion have ce similarities, thus producing specific false recognition.The flow pattern in the process of its accurate flow pattern determination artifi has a particular subjectivity and lack of objective realism, and its error rate will be hi Based on the energy feature values obtained by EEMD decomposition and H transform, the flow pattern identification model constructed using a random forest m can quickly and accurately identify the two-phase flow pattern in the scavenge Machine learning is used to analyze and identify the flow pattern more objectively accurately.

Conclusions
This paper establishes a test system to visualize the oil and gas two-phase flow scavenge pipe and capture a variety of flow types: bubble flow, plug flow, slug flow annular flow by high-speed photography at 370 k; and obtain non-smooth and nonl pressure fluctuation signals by pressure sensors.
To deeply reveal the characteristics of pressure fluctuations and establish relationship with flow patterns, the EEMD was used to decompose the pre fluctuation signals, generate the corresponding IMFs, and perform the Hilbert trans followed by the Hilbert spectrum to quantitatively analyze the changes of amplitud frequency with time, and establish the correspondence between the changes of ampl and the changes of flow patterns in the tube, indicating that the extracted en characteristics can reflect the changes of flow patterns.
Random forest machine learning successfully established the scavenge pipe's o gas two-phase flow pattern identification model using the energy share as the fe value.The recognition rate of the model is close to 98%, and the accuracy of the mo high.However, the flow pattern recognition for the flow pattern transition boundar The flow pattern in the process of its accurate flow pattern determination artificially has a particular subjectivity and lack of objective realism, and its error rate will be higher.Based on the energy feature values obtained by EEMD decomposition and Hilbert transform, the flow pattern identification model constructed using a random forest model can quickly and accurately identify the two-phase flow pattern in the scavenge pipe.Machine learning is used to analyze and identify the flow pattern more objectively and accurately.

Conclusions
This paper establishes a test system to visualize the oil and gas two-phase flow in the scavenge pipe and capture a variety of flow types: bubble flow, plug flow, slug flow, and annular flow by high-speed photography at 370 k; and obtain non-smooth and nonlinear pressure fluctuation signals by pressure sensors.
To deeply reveal the characteristics of pressure fluctuations and establish the relationship with flow patterns, the EEMD was used to decompose the pressure fluctuation signals, generate the corresponding IMFs, and perform the Hilbert transform, followed by the Hilbert spectrum to quantitatively analyze the changes of amplitude and frequency with time, and establish the correspondence between the changes of amplitude and the changes of flow patterns in the tube, indicating that the extracted energy characteristics can reflect the changes of flow patterns.
Random forest machine learning successfully established the scavenge pipe's oil and gas two-phase flow pattern identification model using the energy share as the feature value.The recognition rate of the model is close to 98%, and the accuracy of the model is high.However, the flow pattern recognition for the flow pattern transition boundary has some limitations, and other recognition tools are needed to assist in better recognition.In conclusion, the present work will help identify the flow pattern more objectively and accurately and better understand the flow characteristics of the gas-liquid two-phase flow in the scavenge pipe.

Energies 2023 , 16 ( 4 ) 1 (
16, x FOR PEER REVIEW 5 of Each new data point is classified as the most classified result of the random forest.The final classification result is shown in Equation: the final voting result of the combined classification model, ( ) i h x is a single decision tree model, Y denotes the target variable, and ( ) I ⋅ is an indicator function.

Figure 3 .
Figure 3.The photo of the visualization pipe.

2 . 2 . 2 .
Decomposing the Pressure Signal of Oil Gas Two-Phase Flow in the Scavenge PipeThe oil-gas two-phase flow in the scavenge pipe was tested to obtain the pressure fluctuation signals at both ends of its tube, and the fluctuation of the pressure signals showed relatively unstable characteristics from which information could not be directly obtained.Therefore, EEMD decomposition and Hilbert transformation of the obtained high-frequency pressure signal curves are carried out, which helps to carry out a more indepth analysis.Four typical flow types, namely, bubble flow, plug flow, slug flow, and annular flow, are selected for research, and a set of pressure fluctuation signals are selected for each.The IMFs, Hilbert spectra, and marginal spectra of the pressure fluctuation signals of the four flow types are shown in Figures 4-7.Decomposing the Pressure Signal of Oil Gas Two-Phase Flow in the Scavenge Pipe The oil-gas two-phase flow in the scavenge pipe was tested to obtain the pressure fluctuation signals at both ends of its tube, and the fluctuation of the pressure signals showed relatively unstable characteristics from which information could not be directly obtained.Therefore, EEMD decomposition and Hilbert transformation of the obtained high-frequency pressure signal curves are carried out, which helps to carry out a more indepth analysis.Four typical flow types, namely, bubble flow, plug flow, slug flow, and annular flow, are selected for research, and a set of pressure fluctuation signals are selected for each.The IMFs, Hilbert spectra, and marginal spectra of the pressure fluctuation signals of the four flow types are shown in Figures 4-7.Decomposing the Pressure Signal of Oil Gas Two-Phase Flow in the Scavenge Pipe The oil-gas two-phase flow in the scavenge pipe was tested to obtain the pressure fluctuation signals at both ends of its tube, and the fluctuation of the pressure signals showed relatively unstable characteristics from which information could not be directly obtained.Therefore, EEMD decomposition and Hilbert transformation of the obtained high-frequency pressure signal curves are carried out, which helps to carry out a more indepth analysis.Four typical flow types, namely, bubble flow, plug flow, slug flow, and annular flow, are selected for research, and a set of pressure fluctuation signals are selected for each.The IMFs, Hilbert spectra, and marginal spectra of the pressure fluctuation signals of the four flow types are shown in Figures 4-7.(a)

2 .
Decomposing the Pressure Signal of Oil Gas Two-Phase Flow in the Scavenge PipeThe oil-gas two-phase flow in the scavenge pipe was tested to obtain the pressure fluctuation signals at both ends of its tube, and the fluctuation of the pressure signals showed relatively unstable characteristics from which information could not be directly obtained.Therefore, EEMD decomposition and Hilbert transformation of the obtained high-frequency pressure signal curves are carried out, which helps to carry out a more indepth analysis.Four typical flow types, namely, bubble flow, plug flow, slug flow, and annular flow, are selected for research, and a set of pressure fluctuation signals are selected for each.The IMFs, Hilbert spectra, and marginal spectra of the pressure fluctuation signals of the four flow types are shown in Figures 4-7.

Figure 4 .Figure 5 .
Figure 4. Bubble flow signal decomposition.(a) IMF components; (b) Hilbert spectrum; (c) Hilbert marginal spectrum.Figure 4. Bubble flow signal decomposition.(a) IMF components; (b) Hilbert spectrum; (c) Hilbert marginal spectrum.The 12 IMFs obtained from the decomposition are subjected to the Hilbert transform, and the signal changes in the way shown in Section 2.2 to obtain the Hilbert spectrum of the pressure signal.The visible two-dimensional contours of the fluctuations are generated from the Hilbert spectrum.The color bar next to each Hilbert spectrum from dark blue to yellow indicates the change of the amplitude from minimum to maximum, which shows the evolution of the amplitude of the pressure fluctuations with time and frequency, and the non-smoothness of the oscillations can be seen from the color shift of the amplitude color band.Comparing the Hilbert spectra of four flow types: bubble flow, plunger flow, elastic flow, and annular flow, shows that the amplitude of signal fluctuation increases sequentially with the increase of inlet volume.The non-smooth signal is caused by the radial and axial velocity variations of the medium in the sensing area of the pressure sensor and the medium distribution.Figures 4-7 show that the Hilbert spectrum is integrated over time

Figure 8 .
Figure 8.(a) Comparison of Prediction Results; (b) Confusion Matrix for Test Data.

Figure 8 .
Figure 8.(a) Comparison of Prediction Results; (b) Confusion Matrix for Test Data.
t) is the nth IMF obtained by decomposition after m additions of white

Table 1 .
Physical properties of the two-phase medium.
4. Results and Discussion4.1.Experimental Study on Oil Gas Two-Phase Flow in Scavenge Pipe