Abstract
Microcharge induction has recently been applied as a dust detection method. However, in complex environments, the detection device can be seriously polluted by noise. To improve the quality of the measured signal, the characteristics of both the signal and the noise should be analyzed so as to determine an effective noise removal method. Traditional removal methods mostly deal with specific noise signals, and it is difficult to consider the correlation of measured signals between adjacent time periods. To overcome this shortcoming, we describe a method in which wavelet decomposition is applied to the measured signal to obtain sub-band components in different frequency ranges. A time-lapse Pearson method is then used to analyze the correlation of the sub-band components and the noise signal. This allows the sub-band component of the measurement signal that has the strongest correlation with the noise to be determined. Based on multifractal detrended fluctuation analysis combined with empirical mode decomposition, the similarity between the signal sub-band components and the noise sub-band components is analyzed and three indices are employed to determine the multifractal characteristics of the sub-band components. The consistency between noise components and signal components is obtained and the main signal components are verified. Finally, the sub-band components are used to reconstruct the signal, giving the noise-free measured (microcharge induction) signal. The filtered signal presents smoother, multifractal features.
1. Introduction
One of the notable properties of dust is its electrification during movement. There has been considerable research on the chargeability of dust and how to measure the electrostatic signal [1–6], and it has been widely used in the detection of gas flow velocity in the field of pneumatic conveying. Yan et al. [7–9] conducted extensive research on electrostatic induction for measuring gas-solid two-phase flows and proposed various types of measurement electrodes and measurement methods. In recent years, soft measurement and machine learning [10] have been applied to measure the parameters of gas-solid two-phase flows by the induced charge method. Zheng et al. [11] studied the electrification of sand dust and the spatial variation of the particle concentration along the flow direction by assuming a certain charge-to-mass ratio. Xu et al. [12, 13] proposed a new type of electrostatic sensor system for measuring the pulverized coal speed and relative mass flow, which determines the speed of pulverized coal from the autocorrelation of the output signal of the electrostatic sensor array. The relative concentration and mass flow rate are then obtained from the root mean square value output by two ring electrodes. The present authors [6] have conducted research on the theory of mine dust belt motors, charging modes, charging measurement models, and measuring devices and have identified nonlinear characteristics in the dust electrostatic signal that reflect the multiscale information contained in the signal. Because the measurement environment of the detection device is complex and changeable, the measured dust electrostatic signal always contains noise. The noise that exists when the detection device is working is also called measurement noise (in the remainder of this paper, the electrostatic signal of dust is referred to as the measurement signal). This noise seriously affects the measurement accuracy and must be removed.
The processing of electrostatic signals obtained by electrostatic sensors is mostly achieved by adding anti-aliasing or low-pass filters to the hardware circuit. There are relatively few signal processing methods for electrostatic signals. Wang et al. [14] compared the frequency spectra of the measured signals from external-surface electrostatic sensors and embedded electrostatic sensors and used a harmonic wavelet transform to decompose the measured signal. The frequency spectrum of inner flush-mounted electrostatic sensors was found to contain a critical frequency for decoupling the induced charge signal from the transferred charge signal, thus improving the cross-correlation speed of the signal. Kamel and Yan [15] proposed a hybrid removal method that combines the cut-off frequency method for removing the noise component outside the signal bandwidth with a median filter for smoothing the signal, thus minimizing the effect of noise in the signal. Wang [16] proposed a radial vibration measurement method based on electrostatic sensing and Hilbert–Huang transform (HHT) signal processing, which provides a simple and economical method for rotating-shaft radial vibration measurements. Zhang [17] used the HHT to study the characteristics of different electrodes’ electrostatic signals and extracted the features that characterize the changes in solid concentration. Wang et al. [18] proposed an adaptive decomposition method based on the discrete wavelet transform (DWT), which effectively reduces the average relative standard deviation of the solid velocity. The abovementioned studies can be combined with modern signal processing methods to achieve electrostatic signal processing, but the basic characteristics and frequency range of noise signals are preset, and so adaptive noise recognition and processing cannot be performed. To remove noise, it is necessary to separate the signal from the noise, or to identify the obvious differences in their respective characteristics. The multilayer wavelet decomposition transform is an effective signal-noise separation method in both the time domain and frequency domain. Therefore, many studies [19–23] have used the multilayer wavelet decomposition transform to remove noise. Wavelet transforms can decompose the noise and the signal into different frequency bands, but they do not have the ability to adaptively identify which sub-band is the noise.
The fractal theory proposed by Mandelbrot [24] provides the possibility of analyzing the similarity and correlation between the measured signal and the noise signal. In particular, multifractals [25, 26] enable a detailed description of the local characteristics of the signal. Combining the theory of multifractals and detrended fluctuation analysis, Kantelhardt et al. [27] developed multifractal detrended fluctuation analysis (MFDFA). MFDFA has since been widely applied to stock market analysis [28], temperature time series [29], seismic wave signals [30], vibration tomographic diagnosis [31], image processing [32], and voice signal analysis [33]. MFDFA can effectively describe the nonlinear measurement signal, especially the multifractal characteristics of the time series, but the analysis of the time series signal requires a detrending process, which causes pseudo-fluctuation errors to appear. There are two main reasons for this: one is that the sequence is over- or underfitted due to the uncertainty of the order of the fitted polynomial function; the other is that MFDFA uses a uniform sequence for data segmentation, resulting in sequence segmentation points that are not continuous. To solve these problems, Li et al. [34] proposed an MFDFA algorithm based on empirical mode decomposition (EMD) and template movement. This algorithm was used to analyze the multifractal characteristics of harmonic signals, resulting in a new method for determining the characteristics of harmonic signals. This method can be used to analyze the similarity between the noise and signal sub-bands given by wavelet analysis, allowing the signal component and noise component to be identified. The noise component can then be filtered out and the signal component reconstructed, thus realizing signal denoising.
In order to separate noise signal and measurement signal, this article proposes an asymmetric two-electrode measuring device, most signal components of one electrode are electrostatic signals, and the other contains most noise signals. A wavelet transform is used to decompose the electrostatic signals and noise signals, and a time-lapse Pearson correlation analysis method is developed to analyze the correlation between the signal wavelet and the noise sub-bands to confirm which signal components have the most correlation with noise. The result of time-lapse Pearson determines that is a noise component, but it is impossible to determine which other components are noise and which are signals. So the EMD-MFDFA algorithm is then used to analyze the similarity between the signal sub-band and the noise sub-band and determine the wavelet sub-bands in which the electrostatic signal is located, and the results show that , , , and are the signal components. Finally, the noise is filtered and the denoised electrostatic signal is reconstructed with , , , and .
2. Analysis Method Combining Multifractals and Empirical Mode Decomposition
2.1. Generalized Multifractal Detrending Algorithm [27]
(1)Determine the “contour” of as where .(2)Divide the contour into equal-length nonoverlapping segments with a length of . There may be a segment with a length of less than s at the end of the contour. To eliminate the influence of this part, the same process is repeated from the opposite end to obtain data segments.(3)Calculate the local trend of each of the segments through the least-squares fitting of the sequence, and then determine the variance as where each segment is denoted by , , and where . is the fitting polynomial of the data segment . Linear, quadratic, cubic, or higher-order polynomials can be used in the fitting process (generally called DFA1, DFA2, DFA3, and so on). As the time series is detrended by subtracting the polynomial fitting from the contour, the ability of the different-order DFAs to eliminate the trend of the series is different. In m-order MFDFA, the m-order trend (or, equivalently, the ()-order trend in the original sequence) in the contour is eliminated.(4)Take the average of all line segments to obtain the q-order wave function where is a q-order function of . As s increases, has a power relationship with s: where is the Hurst index, which is related to the multifractal quality index according to(5)According to the Legendre transformation, the relationship between the singular index and the multifractal spectrum and is
2.2. Empirical Mode Decomposition Algorithm [35, 36]
EMD enables the smoothing of nonstationary signals, with fluctuations and changing trends on different characteristic time scales identified by decomposing the signal. The decomposition result is composed of multiple intrinsic mode functions (IMFs) and residuals , where is a monotonic function that is the average trend of the signal.
The EMD process can be described as follows:
2.3. Improved Method Based on EMD-MFDFA
Analyzing the third step of the MFDFA algorithm, it can be seen that the order of the fitting polynomial cannot be determined at the initial stage of signal analysis. is obtained by curve fitting. If the order is very high, it will take a long time to determine the order, and the determined polynomial may not be optimal. Improper fitting of the volatility trend will introduce erroneous volatility and affect the analysis results. In the EMD result, reflects the average trend characteristics of the signal. Using instead of in the third step of the MFDFA constitutes a new EMD-MFDFA algorithm.
Third step: the local variance of each of the segments iswhere each segment is denoted by , , andwhere .
The first, second, and fourth steps of the EMD-MFDFA algorithm are the same as in the MFDFA algorithm. The following research compares the effects of the MFDFA algorithm and the EMD-MFDFA algorithm in multiple analysis applications.
2.4. Time-Lapse Pearson Analysis Method
The Pearson correlation coefficient describes the correlation between two variables at the same time and does not consider the time lag of the change between variables. This paper proposes a Pearson correlation coefficient method with time-lapse correlation analysis, which enables the time-lapse correlation between two variables to be analyzed.
The Pearson correlation coefficient between samples X and Y of a total of N samples is given bywhere , , , and .
The range of Pearson’s correlation coefficient r is [–1, 1]; |r| represents the degree of correlation between two variables, with |r| closer to 1 indicating that the correlation between the two variables is higher, and the relationship between them is closer.
The time-lapse Pearson correlation between and is defined aswhere , , , and .
The time-lapse Pearson correlation matrix is
3. Signal Acquisition
The study in [15] describes a number of treatments for weak signal-detection circuit design and noise shielding and finds that the measurement signal is still inevitably subjected to interference by surrounding noise signals. Obviously, in addition to the necessary hardware protection to reduce noise interference, noise removal is required. First, it is necessary to understand the morphological characteristics of the noise signal obtained by the detection circuit, and to distinguish the noise components in the measurement signal. The noise in the measurement signal can then be effectively removed.
3.1. Asymmetric Two-Electrode Measuring Device
To remove noise, an asymmetric two-electrode measurement unit is developed. One electrode is the noise collection terminal, which mainly collects noise signals in the detection environment, and the other electrode is the signal collection terminal, which mainly collects noise-containing microelectric signals. According to the positive correlation between the induction signal and the size of the electrode [6], it can be assumed that the induction signal obtained by the noise-collecting electrode is much smaller than the induction signal measured by the measuring electrode and that noise is the main component. The noise-collecting electrode is much smaller than the measuring electrode, having a length of 20–30 mm and a diameter of 1–2 mm. The measuring electrode has a length of 400 mm and a diameter of 10 mm. The structure of the measuring device is shown in Figure 1.

3.2. Acquisition of Measurement Signal
Figures 2(a) and 2(b) show the dust microelectric signal collected by the measuring device. It can be seen that the signal induced on the measuring electrode by the movement of the dust is much larger than the signal collected by the noise-collecting electrode. In the measurement signal recorded by the noise electrode, the effective induction signal is almost completely submerged in the noise.

(a)

(b)
4. Characteristic Analysis of Measurement Signal and Measurement Noise Signal
Figure 2 illustrates that the measurement signal exhibits nonlinear and random characteristics, as does the noise in the measurement environment. To effectively filter the noise from the measurement signal, the measurement noise needs to be identified. First, wavelet decomposition is applied to the signal recorded by the measuring electrode to obtain the frequency components in different frequency ranges. Time-lapse Pearson analysis is then performed on the different frequency components and noise signals to identify those components that have the strongest correlation with the noise. Using the MFDFA algorithm and the improved EMD-MFDFA algorithm, the similarity between the different frequency bands of the signal and the measurement noise is then analyzed, and the frequency components closest to the collected noise are removed. This analysis makes it possible to determine the characteristics of the noise that affect the measurement signal, and to obtain the noise signal as the basis for noise filtering. The architecture is shown in Figure 3.

4.1. Wavelet Analysis of Measurement Signals
4.1.1. db24 Wavelet Decomposition
The application of the db24 wavelet with five-layer decomposition to the measured signal shown in Figure 2(a) gives the signal waveforms of each sub-band, as shown in Figure 4.

The frequency range of each sub-band is listed in Table 1. Analysis shows that each decomposed sub-band contains corresponding frequency components, and the signal amplitudes of the high-frequency signal parts , and are small, but not negligible.
4.2. Time-Lapse Pearson Correlation Analysis of the Measured Signal Sub-Band and the Noise Signal
4.2.1. Analysis of Time-Lapse Pearson Correlation between the Measured Signal Sub-Band and the Noise Signal
Figure 5 shows the result of direct analysis of the correlation between the measured signal sub-band and the noise signal for about 10s (5000 data points). The focus is whether the change in noise occurs first, or whether the detection signal leads.

From Figure 5, the following points can be identified: (1) the time-lapse correlation coefficient between component and the noise is greater than 0.6, indicating a strong correlation between the two; (2) the two signals appear to change simultaneously (only 1-point difference, 0.002 s), and the frequency range of is 29.75–59.5 Hz.
4.2.2. Wavelet Sub-Band Decomposition of Noise Signal
By analyzing the time-lapse Pearson correlation between the wavelet sub-band of the measured signal and the measured noise, it is clear that there is a strong correlation between and the noise, although the two are not completely correlated. The noise contains a strong 50 Hz component, which masks the correlation between other components and the noise. First, the same method as for measurement signal analysis is used to decompose the noise signal into its sub-bands; the wavelet sub-band ranges of the noise are the same as in Table 1.
The time-lapse Pearson correlation between the measured signal wavelet sub-bands and the noise wavelet sub-bands is presented in Table 2. The following correlations can be identified:(1)– and – are highly correlated and exhibit the characteristics of simultaneous changes(2)– indicates that, given the strong correlation between and the overall noise signal, component can be considered as a noise signal(3)–, –, and – exhibit a certain correlation, but it is impossible to directly determine which is the noise component and which is the signal
4.3. Multifractal Characteristic Analysis of Signal and Noise
The MFDFA algorithm is used to study the multifractal characteristics of the measured signal and the measured noise, which is judged as follows.
Relationship between q and : for , the Hurst exponent has the following properties [37]: when = 0.5, the signal sequence has no correlation; when >0.5, the signal sequence has a positive long-range correlation, indicating that the time series is positively correlated; when <0.5, the time series is negatively correlated. If q is not related to , the signal is monofractal; if q is related to , the signal is multifractal.
Relationship between q and : if is a straight line, then the signal is simplex; if q and are nonlinear, then the signal is multifractal.
Relationship between and : if is a constant, the signal is a single-shaped signal; if the curve between and has a single-peak bell shape, the signal is multifractal.
MFDFA analysis was applied to the measured signal and the measured noise signal, and the results are shown in Figure 6. Both and indicate that the measured signal and measured noise do not have multifractal characteristics.

The multifractal features of the signal are submerged by noise pollution. Next, by analyzing the multifractal features of the measured signal and measured noise on each wavelet sub-band, it is possible to determine which sub-bands have noise signals as their main components and which are mainly measured signals.
4.3.1. Analysis of the Detrending Ability of MFDFA and EMD-MFDFA Algorithms
The curve of the wavelet sub-band d component of the signal was obtained using MFDFA and EMD-MFDFA. The detrending polynomial of the MFDFA algorithm is linear, and the EMD-MFDFA detrending term is the residual of EMD. The results are shown in Figure 7. From the value of in Figure 7, we find that the absolute value of the result obtained by EMD-MFDFA and its numerical range are smaller than those obtained by MFDFA. This shows that EMD-MFDFA has stronger detrending ability than MFDFA.

4.3.2. Analysis of MFDFA and EMD-MFDFA Characteristics: Wavelet Sub-Bands of Measurement Signals and Wavelet Sub-Bands of Noise Signals
(1) q-h(q) Characteristic Analysis. The wavelet sub-bands of the measured signal and the measured noise obtained from MFDFA and EMD-MFDFA are now analyzed. The results are shown in Figure 8. For , the Hurst exponent of each wavelet sub-band component of the measured signal and noise and its long-range correlation are presented in Tables 3 and 4. Analysis of the MFDFA and EMD-MFDFA algorithms shows that EMD-MFDFA can effectively eliminate the influence of the trend item, so that the trend characteristic of the affected component is eliminated. This highlights the multifractal characteristics of the component. By comparing the values of signal components and in Figure 8, it can be seen that EMD-MFDFA highlights the multifractal characteristics of the measured signal components. EMD-MFDFA gives smaller k values for acoustic components , , and than MFDFA, which leads to the fact that the measured signal does not have multifractal characteristics.

(a)

(b)
(2) Characteristic Analysis. Next, MFDFA and EMD-MFDFA were used to obtain the wavelet sub-bands of the measured signal and the measured noise, and the linearity R between and was calculated. The results are presented in Tables 5 and 6. Comparing the relationship between R and , it is apparent that values of greater than 0.96 indicate multifractal characteristics. EMD-MFDFA effectively reduces the trend characteristics of components and highlights their fractal characteristics.
(3) Characteristic Analysis. MFDFA and EMD-MFDFA were then used to obtain the wavelet sub-band signals of the measured signal and the measured noise, and the waveforms were calculated. The waveforms are shown in Figure 9, and the fractal feature judgments are presented in Tables 7 and 8.

(a)

(b)
(4). Determination of Signal Wavelet Sub-Band. The correlation, average error, and standard deviation are now used to evaluate the differences among the EMD-MFDFA analysis results of each sub-band using , , and . The results are presented in Tables 9–11. Combining the results of Figures 8–10, it can be seen that the multifractal features of , , , and have very high consistency. These components are mainly concentrated at low frequencies, so , , , and can be considered as real signal components of the measured signal.

(a)

(b)
5. Signal Noise Filtering
Combining the time-lapse Pearson analysis results and the analysis results using the EMD-MFDFA algorithm, it can be seen that component is highly correlated in the wavelet decomposition of the measured signal and the measured noise. This can be considered as the main noise component, and the results of multifractal analysis show that , , , and are the main real signal components. The wavelet decomposition components and of the two channel acquisition signals exhibit disorder. This indicates that the signal is polluted by disordered noise.
The , , , and components can be used to reconstruct a new signal, regarded as an electrostatic signal after noise filtering. The filtered waveform is shown in Figure 11. It can be seen that the filtered signal is smoother than the original signal, as the interference of power frequency noise and disordered noise has been eliminated, and provides a basis for further signal analysis. Even if the measured noise is completely submerged by noise, the real charge-induced signal can be obtained by filtering, as shown in Figure 12.

(a)

(b)

The filtered signal was analyzed using EMD-MFDFA, and the curve was found to exhibit multifractal characteristics (see Figure 13). This indicates that the filtered signal has multifractal characteristics.

6. Conclusion
(1)First, the measurement model was improved, and a signal acquisition method using an asymmetric two-electrode structure was proposed. In this method, large-sized electrodes are used to obtain the induced signals and small-sized electrodes are used to record the noise.(2)A method of combining EMD-MFDFA with time-lapse Pearson correlation analysis has been proposed. Based on the db24 wavelet decomposition method, the measured signal and measured noise were decomposed. The results of the time-lapse Pearson correlation analysis showed that the component of the induced signal wavelet decomposition had a strong overall correlation with the noise and exhibited the characteristics of simultaneous changes, without any delay or lag.(3)Further, based on the wavelet decomposition of the noise signal, time-lapse Pearson correlation analysis found that –, –, –, –, –, and – exhibited some correlation. However, it was impossible to directly determine which components were noise and which were part of the signal using this method.(4)Multifractal feature analysis of the measured signal found no obvious multifractal features. MFDFA and EMD-MFDFA analysis of each component after wavelet decomposition did, however, find obvious multifractal features in several components.(5)Comparing the multifractal results of the measured signal wavelet sub-band and the measured noise wavelet sub-band, there is obvious consistency between –, –, –, and –. It was verified that the components of the measured signal were mainly concentrated at low frequencies, and , , , and were considered to be the main components of the real signal.(6)The signal was reconstructed using components , , , and to obtain the noise-removed induction signal, which presents smoother characteristics. The filtered signal also presents obvious multifractal characteristics.In future work, the denoised electrostatic induction signal will be used for the measurement of dust concentration. Simultaneously, by combining Ensemble Empirical Mode Decomposition [36] (EEMD), Variational Mode Decomposition [38] (VMD) with stronger ability to extract trends, MFDFA’s de-trending capabilities will be further improved.
Data Availability
The data used to support the findings of this study are included within the supplementary information file.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the National Natural Science Foundation of China Program (no. 62073198) and the Major Research Development Program of Shandong Province of China (no. 2016GSF117009).
Supplementary Materials
“Measured_noise_data.docx” includes the measured noise used in this paper. “Measured_signal_data.docx” includes the measured signal used in this paper. (Supplementary Materials)