1. Introduction
The distillation column is characterized by a large scale, continuity and multiple variables. It is widely used in chemical production and separation processes [
1,
2]. Many control parameters must be regulated in order to maintain the balance of the operation process [
3,
4,
5,
6,
7]. When the operating parameters change or malfunction occurs, the distillation column loses its original balance and enters an abnormal operation state [
8], which threatens the safety of the chemical production process if it is not solved in time. Therefore, it is important to develop quick fault diagnosis method of fluid flow in columns to reduce the probability of accidents [
9].
At present, the monitoring methods of chemical equipment mainly include ray scanning and ultrasound techniques [
10,
11]. However, these methods are not the ideal choices for column operation monitoring because of their high cost of operation and restrictions on application conditions. Recently, the rapid development of processing technology has promoted the application of acoustic signal in various fields [
12,
13,
14]. As a new technology, acoustic detection technology has been successful in the detection of chemical dynamic equipment [
15,
16].
Compared to the obvious mechanical acoustics of dynamic equipment, as static equipment, the distillation column has less obvious acoustics. The acoustics of fluid flows in sieve tray distillation columns are affected by numerous factors, such as the gas–liquid two phase interaction, and fluid interaction with the column wall [
17,
18]. These interactions are the main source of acoustics from the operation of the distillation column. Different operating states correspond to various acoustic characteristics. The feature database can be constructed by extracting acoustic features in different flow states. Based on the database, the unknown states can be recognized by a machine learning classification algorithm.
In this work, the microphone array was installed on the outer wall of the sieve tray distillation column to collect the acoustic signals of different flow states. Then, the acoustic signals were preprocessed [
19] and the characteristic parameters of the acoustic signals were extracted [
20]. Next, the acoustic feature database containing the flow states label was constructed. Finally, the acoustic method based on the k-NN algorithm was used to identify the fluid flow states. The acoustic method could judge the operation states, which provided a strong guarantee for the safe operation of the distillation column.
This paper is divided into six sections.
Section 2 introduces the initial state and text process of the experiment and test equipment.
Section 3 presents the feature extraction method of an acoustic signal.
Section 4 describes the theory of the k-NN algorithm and the selection process of the
k value.
Section 5 shows the experimental process and analyzes the experimental data.
Section 6 offers a synthesis of the objective and importance of this work.
2. Experiment and Test Device
The diameter of the distillation column is 50 mm, the number of trays is 9, the system is 15% ethanol- water solution, and the total reflux operation is carried out. The condensing capacity is sufficient for the system to achieve the apparent flooding operation by adjusting the amount of heating steam. The experimental device and acoustic measuring element are shown in
Figure 1. The geometrical parameters of the sieve tray distillation columns and operational parameters are shown in
Table 1.
Acoustic signal acquisition process: Place the recording pen vertically on the outer wall of the column insulation layer. The hand-held acoustic meter should be close to the two-layer tray. When the distillation column is fully refluxed to the state of flooding, recording begins; the same position is recorded and saved many times. The acoustic signal sampling frequency is 8000, and the sample quantity includes 45 normal data and 55 abnormal data.
3. Feature Extraction of Acoustic Signal
3.1. Double-Domain Feature Extraction Method
The original acoustic signal has many characteristics. Dual-domain feature extraction is used to select suitable features for research and classification. First, the energy ratio of the acoustic signals is judged based on the time domain characteristics. Then, the sixth-order linear prediction coefficient (LPC) of the acoustic signal is obtained based on the frequency domain characteristics.
3.1.1. Analysis of Sectional Energy Ratio
This paper studied the energy ratio of an acoustic signal produced by different attenuation speeds. Since the acoustic signal is a slow-changing short-time stable signal, the energy ratio at different time periods should be different. In order to improve the total energy ratio which could not be better used for comparison, a segmented energy ratio was adopted [
21]. First, the acoustic signal of the column was divided into three parts and the energy of each part was calculated. Then, the energy ratio could be obtained by comparing the total energy of the acoustic signal with the total energy of the column. The expression is as follows:
where
is the total energy of the input signal. The three sections of energy are represented by
,
and
, respectively. The energy proportion of each segment is
,
and
.
3.1.2. LPC
LPC is a parameter model based on acoustic synthesis. The purpose of linear prediction is to predict the present or future sample values. Therefore, the current sampling value of the input signal is calculated by the relationship between the past sampling value, the current sampling value and the future sampling value [
22]. In this paper, the relationship between the two was closer so that the error was also reduced to a minimum. Finally, a unique set of sixth-order linear prediction coefficients is obtained [
23].
For the sampled acoustic signal
at a given time, the linear combination of the preceding
samples can be used for prediction. The expression is as follows:
where
,
, …
is the linear prediction coefficient of the previous
samples, and
is the prediction error. Z transforms
to obtain
.
The transfer function of the linear predictive synthetic filter system can be written as
3.2. Analysis and Comparison of Characteristic Parameters
The acoustic signals of different parts of the column were extracted; for example, the acoustic signals of the normal operation of the upper layer of the tray tower, the overflow of the upper layer and the overflow of the lower part of the tray tower were extracted, respectively. The dual-domain feature extraction method was used to judge the energy ratio of the acoustic signal from the perspective of the time domain, and then to process the feature of the acoustic signal according to the sixth-order LPC based on the frequency domain characteristics. This method can realize the feature extraction of all acoustic signals through nine feature vector parameters. The acoustic signals in different states can be distinguished according to these eigenvector parameters.
4. k-NN Algorithm
4.1. Basic Process of the k-NN Algorithm
The k-NN algorithm is a non-parametric regression classification algorithm, which means that the category of any sample can be judged by its adjacent k samples. First, the algorithm selects a feature as the similarity measure and specifies the distance rule. Next, the k adjacent samples are found by calculating the similarity between the test sample and the training set sample to form a neighborhood of the samples to be tested. Finally, the classification of the test samples is determined by the category proportion of the most training samples in the neighborhood.
In this experiment, k-NN was realized through three steps. The first step was to calculate the distance d between the acoustic signal to be measured and the acoustic signal of each sample. Specifically, the center of the acoustic signal to be measured was used to find the
k neighboring points around it. The Euclidean distance between two sampling points was calculated, as shown in the equation below.
where
are the coordinates of its reference point. In addition, the signals are arranged in order of distance from small to large.
The second step is to arrange the signals in order of distance from the smallest to the largest. Then, the k points closest to the acoustic signal to be measured are selected as the reference points.
In the last step, the number of classes to which the
k reference points belong is compared. According to the principle of maximum winning, the test samples are classified into the winning category. The adjacent points are divided into two types: normal acoustic signal and fault acoustic signal, as indicated in Equation (6).
is the k-NN including
y, and
C is the label. When the label of
y is
j, the return value is 1 and the value of
is “true”. Otherwise,
[
24].
4.2. The Value of k
The k is the proximity number. The category of the acoustic signal to be measured is determined by selecting the k adjacent points.
Selection of the k value is very important. If the value of k is very small, the noise will have a great impact on the prediction results. The reduction in the k value is easy to overfit, which increases the complexity of the classification model. If the value of k is large, it is equivalent to using the acoustic signal in the larger neighborhood to predict the target point, which will increase the error.
The common method is to use the test set to estimate the classifier error rate from k = 1, repeat the process and increase k by 1 each time. The k value that produces the minimum error rate should be selected. Generally, the value of k does not exceed 20 and the upper limit is the root of N. With the increase in the data set, the value of k also increases.
4.3. K-Fold Cross-Validation Method
K-fold cross-validation was used to randomly divide data set A into k packets, one of which was used as a test set, and the remaining (k − 1) packets were used as a training set.
The specific operation process of the K-fold cross-validation method is as follows.
Data set D is split by k packets, and the process of splitting is random.
One of the k packets is used as the test set, and the remaining (k − 1) packets are used as the training set. The data set is divided into feature (train_ x), training set label (train_ y) and test set feature (test_ x). The k-NN model is trained by train_ x and train_ y, and then the trained model is used to predict train_ x.
The average of k times classification rate is the real classification rate of the model (hypothesis function).
Figure 2 shows the prediction probability corresponding to different
k values. We found that when the
k value was 5, the maximum prediction probability was 88.2353%; therefore, we chose the
k value as 5.
5. Experimental Data Analysis
5.1. Experimental Process
First, the acoustic signals obtained from the distillation column of Hebei University of Technology were used to establish a characteristic acoustic parameter database and verify the acoustic recognition method. The following steps were performed, as shown in
Figure 3.
For acoustic signals, 100 acoustic signals in two operating states (normal and flooding) were recorded, respectively.
MATLAB was used to convert all analog signals into digital signals to achieve A/D conversion.
The characteristic information of the acoustic signal was extracted from the time domain and frequency domain.
The characteristic parameter data of 100 acoustic signals were divided into a training set and test set. The ratio of training set to test set was 7:3. Each group included the characteristic value of the signal and the corresponding operation state. The training set was used as the signal feature database.
The k-NN was constructed by using the signal characteristic database of three flow states. Then, the characteristic parameters of the test set were input into the algorithm to predict the operation state of the distillation column. The precision of k-NN was calculated according to the matching degree between the predicted state and the actual state.
5.2. Experimental Result
Three kinds of acoustic signals of different parts and different operating states of the column were extracted and divided into two angles. First, from the perspective of the time domain, the energy ratio of the acoustic signal was used to judge. Then, from the perspective of the frequency domain, the acoustic signal was processed according to the sixth-order LPC. The characteristic vector parameters of the distillation column’s acoustic signal under different operating conditions are shown in
Table 2. The waveform of the original acoustic signal and the spectrum of linear prediction are shown in
Figure 4.
Figure 4a shows the waveform of the original acoustic signal and the linear prediction spectrum of the upper layer’s normal operation;
Figure 4b shows the waveform of the original acoustic signal and the linear prediction spectrum of the flooding of the upper layer;
Figure 4c shows the waveform of the original acoustic signal and the linear prediction spectrum of the flooding of the lower layer.
Table 2 shows that the sectional energy ratio of the acoustic of ‘Flooding of the upper layer’ increased first and then decreased, and its attenuation degree was convex. There was no large change in the section energy ratio of ‘normal operation of the upper layer’, as the section energy ratio represents the attenuation speed of the acoustic signal. The attenuation speed of the acoustic signal of the ‘normal operation of the upper layer‘ was very slow, while the section energy ratio of the ‘Flooding of the lower part’ accounted for 55.26% in the first section. Therefore, the two types of acoustic can be distinguished according to the segmented energy ratio.
The double-domain feature extraction method was analyzed taking the distillation column as an example. The nine feature vector parameters were applied to the feature extraction of the column’s acoustic signal.
5.3. Analysis of Test and Simulation Results
5.3.1. Testing of Simulation Results
Ten independent experiments were conducted on 70 types of acoustics, and the experimental results were the same. As shown in
Table 3, among the 70 groups of acoustic signals to be measured, five groups of abnormal signals were wrongly judged as normal signals, three groups of normal signals were wrongly judged as abnormal signals, and the prediction results of the remaining 62 groups were correct, so the accuracy of the simulation experiment was 88.57%.
5.3.2. Error Analysis of Test Simulation Experiment Results
Three typical acoustic signals were selected from the test data set for a comparative analysis. ‘Signal a’ and ‘Signal b’ were signals when the upper layer was seriously flooded, and ‘Signal c’ was the test signal when the column was in normal operation at the same position. The time–frequency domain characteristics comparison diagram of three test signals is shown in
Figure 5.
Figure 5a is composed of three subgraphs showing the waveform, linear prediction and spectrum of the ‘Signal a’ original acoustic signal from left to right;
Figure 5b shows the waveform, linear prediction and spectrum diagram of the ‘Signal b’ original acoustic signal; and
Figure 5c shows the waveform, linear prediction and spectrum of the ‘Signal c’ original acoustic signal. The time domain characteristic parameter used in the experiment was the sectional energy ratio, and the frequency domain characteristic parameter was the sixth-order LPC. The characteristic parameters of the three test signals are shown in
Table 4. The actual experimental result shows that ‘Signal b’ was judged as flooding, which is correct; however, ‘Signal a’ was incorrectly judged to be in a normal operation state.
In
Table 4, the sectional energy ratio of ’Signal a’ decreases in turn, and the gap is very small, which is highly similar to the ’Signal c’ sectional energy ratio, while the sectional energy ratio of ’Signal b’ is larger at both ends and smaller in the middle, which is the first reason for miscalculation. In
Figure 5, the linear prediction spectrum of ’Signal a’ has two peaks (4.5, 1.60) (17.4, 2.69), while the linear prediction spectrum of ’Signal b’ has two peaks (3.1, 2.31) (17.1, 2.43), and the peak value of ’Signal c’ is (15.2, 1.30), which is the second reason for miscalculation. The spectrum shows that the components of ’Signal a’ were concentrated in the low frequency and high frequency parts, while ’Signal b’ was concentrated in a small number of high frequency parts, which is the third reason for miscalculation.
6. Conclusions
This paper proposed a fault detection methodology for abnormal distillation column operations using acoustic signals. The k-NN classification algorithm based on the sectional energy ratio and LPC fusion feature vector proved to have high recognition accuracy (88.57%). In addition, the feasibility of acoustic signals monitoring in static equipment operation was demonstrated.
Based on the acoustic signals of a chemical column under different operating conditions, this study mainly monitored and classified the plate tower under normal and flooding operating conditions and demonstrated the feasibility of monitoring the operating conditions of the chemical column based on acoustic signals. In addition to the flooding mentioned in this paper, the chemical column also had many abnormal operating conditions, such as liquid leakage, downcomer blockage, plate hole abnormality, entrainment of mist, etc. The accurate classification of different types of chemical column under various abnormal operating conditions will be the focus of our next research work.
In the future, acoustic signal detection technology can be widely used in the chemical equipment of operation states identification and fault diagnosis, which is of great importance for safe chemical production.
Author Contributions
Conceptualization, G.-Y.W. and H.-H.W.; methodology, G.-Y.W.; software, Z.-H.Y. and Y.Z.; validation, Z.-H.Y. and Z.-X.Z.; formal analysis, Z.-H.Y.; investigation, G.-Y.W. and Z.-H.Y.; resources, Z.-X.Z. and H.-H.W.; data curation, Z.-X.Z. and Y.Z.; writing—original draft preparation, G.-Y.W.; writing—review and editing, Z.-H.Y.; visualization, Z.-H.Y.; supervision, G.-Y.W., H.-H.W. and B.-J.G.; project administration, G.-Y.W. and H.-H.W.; funding acquisition, H.-H.W. All authors have read and agreed to the published version of the manuscript.
Funding
This research was funded by the National Natural Science Foundation of China (21878066/22008050); the Natural Science Foundation of Hebei Province (B2019202167/B2020202066); the Natural science foundation of Tianjin (20JCQNJC00430); the Scientific Research and Innovation of Tianjin (2021YJSS290).
Data Availability Statement
Not applicable.
Acknowledgments
We thank for financial support from the National Natural Science Foundation of China (21878066/22008050); the Natural Science Foundation of Hebei Province (B2019202167/B2020202066); the Natural science foundation of Tianjin (20JCQNJC00430); the Scientific Research and Innovation of Tianjin (2021YJSS290).
Conflicts of Interest
The authors declare no conflict of interest.
Abbreviations
E | the total energy of the input signal |
Ei | Section i energy |
ki | Energy ratio of section i |
S(n) | sampled acoustic signal |
ap | the linear prediction coefficient of the previous p samples |
u(n) | prediction error |
S(z) | Z transform s(n) to obtain S(z) |
H(z) | The transfer function of linear predictive synthetic filter system |
References
- Finberg, E.A.; May, T.L.; Shiflett, M.B. Multicomponent Refrigerant Separation Using Extractive Distillation with Ionic Liquids. Ind. Eng. Chem. Res. 2022, 61, 9795–9812. [Google Scholar] [CrossRef]
- Jiao, Y.; Yin, K.; Liu, T.; Meng, F.; Li, X.; Zhong, L.; Zhu, Z.; Cui, P.; Wang, Y. Process design and mechanism analysis of reactive distillation coupled with extractive distillation to produce an environmentally friendly gasoline additive. J. Clean. Prod. 2022, 369, 133290. [Google Scholar] [CrossRef]
- Rumbo Morales, J.Y.; López López, G.; Alvarado Martínez, V.M.; Sorcia Vázquez, F.d.J.; Brizuela Mendoza, J.A.; Martínez García, M. Parametric study and control of a pressure swing adsorption process to separate the water-ethanol mixture under disturbances. Sep. Purif. Technol. 2020, 236, 116214. [Google Scholar] [CrossRef]
- Rumbo Morales, J.Y.; Brizuela Mendoza, J.A.; Ortiz Torres, G.; Sorcia Vázquez, F.d.J.; Rojas, A.C.; Pérez Vidal, A.F. Fault-Tolerant Control implemented to Hammerstein–Wiener model: Application to Bio-ethanol dehydration. Fuel 2022, 308, 121836. [Google Scholar] [CrossRef]
- Wang, H.H.; Wang, Z.B.; Zhou, Q.; Liang, J.; Yin, Y.; Su, W.Y.; Wang, G.Y. Optimization and Sliding Mode Control of Dividing-Wall Column. Ind. Eng. Chem. Res. 2020, 59, 20102–20111. [Google Scholar] [CrossRef]
- Torres Cantero, C.A.; Pérez Zúñiga, R.; Martínez García, M.; Ramos Cabral, S.; Calixto-Rodriguez, M.; Valdez Martínez, J.S.; Mena Enriquez, M.G.; Pérez Estrada, A.J.; Ortiz Torres, G.; Sorcia Vázquez, F.d.J.; et al. Design and Control Applied to an Extractive Distillation Column with Salt for the Production of Bioethanol. Processes 2022, 10, 1792. [Google Scholar] [CrossRef]
- Švandová, Z.; Markoš, J.; Jelemenský, Ľ. Impact of mass transfer coefficient correlations on prediction of reactive distillation column behaviour. Chem. Eng. J. 2008, 140, 381–390. [Google Scholar] [CrossRef]
- Taqvi, S.A.; Tufa, L.D.; Zabiri, H.; Maulud, A.S.; Uddin, F. Multiple Fault Diagnosis in Distillation Column Using Multikernel Support Vector Machine. Ind. Eng. Chem. Res. 2018, 57, 14689–14706. [Google Scholar] [CrossRef]
- Akula, A.; Goel, S.; Ghosh, R. Condition Monitoring Saves Money and Prevents Failures. Aiche J. 2017, 113, 58–62. [Google Scholar]
- Wibisono, W.; Sugiharto, S. Learning gamma-ray scanning technique through an educational simulation rig. J. Phys. Conf. Ser. 2021, 2104, 012001. [Google Scholar] [CrossRef]
- Du, Y.; Zhou, S.; Jing, X.; Peng, Y.; Wu, H.; Kwok, N. Damage detection techniques for wind turbine blades: A review. Mech. Syst. Signal Pr. 2020, 141, 106445. [Google Scholar] [CrossRef]
- Huang, Z.L.; Peng Zhang, P.; Yang, Y.; Sun, J.Y.; Wang, J.D.; Yang, Y.R.; Huang, C.H.; Huang, Z.X. Acoustic emission detection of slag performance in coal gasifiers. Chin. J. Chem. Eng. 2020, 28, 2200–2206. [Google Scholar] [CrossRef]
- Salvati, D.; Drioli, C.; Ferrin, G.; Foresti, G.L. Acoustic Source Localization from Multirotor UAVs. IEEE Trans. Ind. Electron. 2020, 67, 8618–8628. [Google Scholar] [CrossRef]
- De Roo, L.; Vergeynst, L.; De Baerdemaeker, N.; Steppe, K. Acoustic Emissions to Measure Drought-Induced Cavitation in Plants. Appl. Sci. 2016, 6, 71. [Google Scholar] [CrossRef] [Green Version]
- Zhu, Y.; Li, G.P.; Tang, S.N.; Wang, R.; Su, H.; Wang, C. Acoustic signal-based fault detection of hydraulic piston pump using a particle swarm optimization enhancement CNN. Appl. Acoust. 2022, 192, 108718. [Google Scholar] [CrossRef]
- Gong, C.A.; Lee, H.C.; Chuang, Y.C.; Li, T.H.; Su, C.S.; Huang, L.H.; Hsu, C.W.; Hwang, Y.S.; Lee, J.D.; Chang, C.H. Design and Implementation of Acoustic Sensing System for Online Early Fault Detection in Industrial Fans. J. Sens. 2018, 2018, 4105208. [Google Scholar] [CrossRef]
- Syeda, S.R.; Afacan, A.; Chuang, K.T. A Fundamental Model for Prediction of Sieve Tray Efficiency. Chem. Eng. Res. Des. 2007, 85, 269–277. [Google Scholar] [CrossRef]
- Zarei, A.; Hosseini, S.H.; Rahimi, R. CFD and experimental studies of liquid weeping in the circular sieve tray columns. Chem. Eng. Res. Des. 2013, 91, 2333–2345. [Google Scholar] [CrossRef]
- Kwon, H.; Oh, K.C.; Choi, Y.; Chung, Y.G.; Kim, J. Development and application of machine learning-based prediction model for distillation column. Int. J. Intell. Syst. 2021, 36, 1970–1997. [Google Scholar] [CrossRef]
- Xie, H.; Pierce, L.E.; Ulaby, F.T. SAR speckle reduction using wavelet denoising and Markov random field modeling. IEEE Trans. Geosci. Remote Sens. 2002, 40, 2196–2212. [Google Scholar] [CrossRef] [Green Version]
- Wu, S.; Zhao, C. A Study of Perceptually-Based Features for Objective Speech Quality Evaluation. Acta Electron. Sin. 1999, 27, 92–94. [Google Scholar]
- Wang, R.; Chai, P. Research on Speech Cepstral Features. Comput. Eng. 2003, 29, 31–33. [Google Scholar]
- Makhoul, J. Linear prediction: A tutorial review. Proc. IEEE 1975, 63, 561–580. [Google Scholar] [CrossRef]
- Mitchell, H.B.; Schaefer, P.A. A “soft” K-nearest neighbor voting scheme. Int. J. Intell. Syst. 2001, 16, 459–468. [Google Scholar] [CrossRef]
| Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).