ECT Image Recognition of Pipe Plugging Flow Patterns Based on Broad Learning System in Mining Filling

,e process of mining filling, when the slurry is transported to the goaf by the filling pipeline, is very important to find the location and size of the caking in the filling pipeline in time for the safe and stable operation of the mine filling pipeline. It is an important research work to detect different flow patterns after two-dimensional section reconstruction in closed filling pipeline based on ECT (electrical capacitance tomography) visualization method. Slurry flow in pipeline is regarded as a two-phase flow, and the multishape distribution was reconstructed into images by ECT and intelligently recognized by broad learning system (BLS) algorithm. BLS is a feedforward neural network with few optimization parameters and fast training speed. In this paper, three features of two-phase sample images, the number of regional blocks, the roundness of regional blocks, and barycenter of regional blocks, are combined with network structure of BLS to recognize different flow patterns. ,rough the simulation, the recognition accuracy of two-phase fillback image is more than 99%.,is conclusion indicates the effectiveness of BLS to predict different twophase flow patterns; it also provides a new solution for the pattern recognition of the flow pattern in the mining filling pipeline.


Introduction
In the filling pipeline from the surface to the goaf, due to the unstable slurry flow pattern, there is a large density difference between the upper and lower parts of the pipe. e transport velocity is too low to make the coarse slurry deposit at the bottom of the pipe, which can easily cause serious blockage in the junction section between the horizontal pipeline and the vertical (or inclined) pipeline [1,2]. ECT visualization and BLS image recognition technology are helpful for rapid detection and slurry location in pipe blockage.
As slurry materials and proportions, strength and other factors lead to different flow conditions [3]; ECT visualization of pipeline cross section focuses on the permittivity of the medium in backfill pipeline. A possible option is to dynamically observe the internal conditions of the pipeline in the area where the blockage is most likely to occur a priori. ECT is a visualization technique for two-dimensional or three-dimensional image reconstruction of limited mediums or mixtures in a closed container [4]. ECT technology has the advantages of high accuracy and low cost. Its hardware mainly includes capacitive sensor, acquisition circuit, and computer. e capacitive sensor may be composed of 8 or 12 or 16 electrodes, in which each electrode travels around the closed container at an equal distance and is attached to the outer wall of backfill pipeline. When the ECT system is working, a set of independent capacitance data is collected through different excitation orders of electrodes to the computer for image reconstruction. From this set of independent capacitance data, the permittivity distribution of the slurry flow is inverted by image reconstruction algorithm. ECT technology has great potential for applications such as multiphase flow detection [5], combustion imaging [6], and fluidized bed solid particle monitoring [7]. In this paper, ECT technology is used to detect the insulation of mining fillings. e classic image reconstruction algorithms of ECT technology include linear back projection (LBP) [8], Landweber [9], and Tikhonov regularization [10]. ese algorithms have advantages and disadvantages. For example, using LBP algorithm has fast imaging speed but low accuracy, while the result of using the Landweber algorithm is that the image quality depends on parameters. Accuracy and efficiency have become the main obstacle for the field application of ECT technology, and it is also the main research direction of relevant scientific researchers. Wang Pai used least squares-support vector machines (LS-SVM) method to predict the capacitance error caused by the soft field, the fitness function of annealing particle swarm optimization (APSO) algorithm is constructed according to the capacitance error, and the image reconstruction with high precision is realized by searching for the best resolution [11]. Zhang proposed the ECT image reconstruction algorithm based on the compressed sensing theory and used the discrete Fourier transform basis to sparse the original image gray signal, and l 1 − norm regularization and primal-dual interior point method were used to reconstruct the image [12]. ECT is becoming more mature with the advent of many fast and accurate reconstruction algorithms. ECT algorithms include forward and inverse problem solving. e forward problem refers to the establishment of mathematical model based on the distribution of permittivity in sensitive field and boundary conditions of sensitive field so as to obtain the mutual capacitance between electrode pairs. e inverse problem is to reconstruct the distribution of permittivity in sensitive field by the sensitivity of the field and the measured capacitance data and take the image as the final representation. However, the unprocessed reconstructed images cannot meet the requirements of the site. In order to achieve intelligent and integrated operation flow, ECT technology still needs various postprocessing, among which flow pattern recognition is the basic requirement. AbdeelJ Roman extracted the features of the space average permittivity and the position of the center of mass and combined them and used neural network to automatically identify the liquid-vapor two-phase flow patterns of refrigerant in horizontal tubes. e final experiment showed that the classification results of the nine flow patterns reached an average accuracy of 98.1% [13]. By using long short-term memory (LSTM) method, which can overcome the gradient disappearance or escalation, Rafael Johansen provided a successful application scenario for the flow pattern identification of annular, plug, slug, stratified, and wavy flow regimes [14]. Image reconstruction of twophase flow and flow pattern identification as a link in the field assembly line, their training speed, and classification accuracy has important influence on the subsequent operation.
In this paper, the ECT system is used to reconstruct the two-dimensional cross-sectional image inside the mining filling pipeline, and the broad learning system (BLS) provides a new idea for the flow pattern recognition inside the filling pipeline [15].
BLS is a novel neural network framework, which has the advantages of fewer parameters to be optimized and shorter training time compared with deep structural neural network. e simple BLS network is suitable for processing small sample sets of low dimensions and establishing online learning model efficiently. BLS network structure is flexible, and it can be used in combination with other networks. For example, the features extracted by convolutional neural network (CNN) are combined with the output of BLS feature layer to enhance the ability of network prediction or classification under the premise of optimal training time. In this paper, the idea of combining features is adopted to realize flow pattern identification.
is paper is organized as follows: in Section 2, the basic principles of BLS theory are introduced. In Section 3, three artificial feature extraction methods are described. In Section 4, artificial feature extraction method and BLS are combined to enhance the ability of model recognition.
e reconstruction results are analyzed and evaluated through simulation in Section 5. Concluding remarks are presented in Section 6.

BLS Theory
e predecessor of BLS is random vector functional-link neural network [16] (RVFLNN), which essentially belongs to feedforward neural network (FNN). BLS network consists of mapping layer, enhancement layer, and output layer. A BLS network is shown in Figure 1. First, input data is mapped to "feature nodes" through sparse weight mapping, and these "feature nodes" form the mapping layer. Second, the output of the mapping layer serves as the input of the enhancement layer and is activated by the "enhancement node." Finally, the output of the mapping layer and the output of the enhancement layer are combined as the input vector of the output network, and the corresponding weight is calculated by the Moore−Penrose generalized inverse. e main network relations of each layer are as follows: where X is the input data; Y is known label data; ϕ(·) is the identity mapping function; ξ(·) is the tan sig activation function (tan sig(x) � (2/(1 + e − 2x )) − 1); W e , W h , and W m are weights of mapping layer, enhancement layer, and output layer, respectively; β e and β h are biases; Z and H are the outputs of the mapping layer and the enhancement layer, respectively. Set the input data to X m×n (m is the number of samples in the training set and n is the number of pixels). A column vector whose entries are all constants adds to the end of X m×n and the new matrix is written as X m×(n+1) , so equation (1) can be abbreviated as where W ee is a randomly generated (n + 1) × N 1 weight matrix (N 1 is the number of "feature nodes" in the mapping layer). It is a necessary processing to sparse W ee to effectively reduce the linear correlation of new "feature nodes." e output Z m×N 1 of the mapping layer is the input to the enhancement layer. Similarly, equation (1) can be abbreviated as 2 Advances in Civil Engineering where W hh is the randomly generated weight matrix of (N 1 + 1) × N 2 (N 2 is the number of "enhancement nodes" of the enhancement layer). e orthonormal basis of W hh was used to improve the feature extraction ability of the enhancement layer.
[Z|H] m×(N 1 +N 2 ) is obtained by combining the output features of the mapping layer and the enhancement layer. A(A � [Z|H] m×(N 1 +N 2 ) ) directly establishes the mapping equation (1) with label data Y. In order to solve the weight W m of the output layer, the optimization of equation (1) is established as where λ is a parameter close to the l 2 − norm and regularization method is adopted to prevent model overfitting.
To solve for weight W m is to solve for the Moore−Penrose generalized inverse.
e BLS network model is mainly composed of W1 of the mapping layer, W2 of the enhancement layer, and W3 of the output layer. Since there is no need for backpropagation during training, the training time is very short. In addition, BLS can provide incremental learning similar to migration learning when dealing with big data, including incremental enhancement node, incremental feature node, and incremental input data. Incremental learning avoids data retraining and greatly saves training time.

Feature Extraction
In this paper, the BLS network is used to identify the flow pattern inside the mining filling pipeline, and the characteristics of the filling body are added to the branches of the original BLS network to further improve the recognition accuracy. ese are three different functions and their extraction methods.
3.1. Number of Regional Blocks. Obtaining the number of area blocks is an important index to judge the pipeline blockage. Based on OpenCV library, the binary image of two-phase flow is used to obtain the number of blocks in the pipeline.
3.2. Circularity of Regional Block. Block shape complexity can be described by shape coefficient. For example, the circularity coefficient of stratified flow is small, while that of bubble flow is large.
where S represents the area of the block and C represents the perimeter of the block.
3.3. Barycenter of Regional Block. Barycenter reflects the location information of the region block. For example, pixels barycenter of stratified flow is near the bottom of the pipe, and pixels barycenter of circulation is generally near the center. Barycenter position of the block is calculated by means of moment based binary image of two-phase flow section [17].
where M 00 is the area of the white area of the image block.

Proposed Method
For the three features extracted from sample sets, the number of regional blocks, circularity, and the barycentric coordinates are represented by P 1 , P 2 , P 3 , respectively. After combining them, P � P 1 P 2 P 3 is obtained, where P 1 , P 2 , P 3 is the matrix of m × 1, m × q, and m × 2q, respectively,m is the number of sample sets, and the value of q is equal to the number of categories to be classified. e number of categories in this article is 6; then q � 6 is set.
calculated from the ith sample image, then the calculated value of the circular degree corresponding to the sample is placed in the front position of matrix P (i) 2 , and the remaining position is filled with 0; the  Advances in Civil Engineering barycentric coordinate is placed in the front position of matrix P (i) 3 , and the remaining position is filled with 0. As shown in Figure 2, the BLS mapping layer is augmented with "feature node," and the output of the mapping layer is combined with the feature matrixPto form a new matrix [Z|P]. After the normalization of P, the input of BLS's enhancement layer is also [Z|P]: where H * is the output of the enhancement layer; A * � [Z|P|H * ], representing the input to the output layer. Additional features are added to the BLS architecture to ensure the accuracy of identification; and because the amount of data for additional features is small, it hardly affects the training speed.

e Data Set.
e data set in this article simulates ECT to reconstruct a two-dimensional cross-sectional image of a mining filling pipeline, and the size of each sample image is 64 × 64 pixels.
ere are 6 types of flow patterns: single bubble, double bubbles, three bubbles, four bubbles, stratified, and annular flow. e data set includes training set, training set label, test set, and test set label. e training set has a total of 3000 images, including 500 images of different categories. e test set has a total of 500 images, with 100 images for each category. Part of the training set is shown in Figure 3.
After the pixels of the sample image are compressed to 32 × 32 pixels, and row by row into a 1 by 2 row vector 1 × 1024, the training set is a 3000 × 1024 matrix and the test set is a 500 × 1024 matrix. e training set label and the test set label are matrices of 3000 × 6 and 500 × 6, respectively. e elements of the label matrix are 0 and 1. e row corresponds to different samples and the column corresponds to different categories.

Feature Representation.
Based on the feature extraction method proposed in Section 4, three features are extracted from the reconstructed binary image by ECT. Figure 4 shows a sample of stratified flow. e image size is 64 × 64 pixels. e number of regional blocks is 1, the barycenter coordinate is (31,45), and the circularity is 0.5344. It can be seen from this sample that the barycenter of the stratified flow image is downward, and the circularity is small.

Simulation Result.
In order to verify the influence of additional features on BLS network flow pattern identification, this paper discusses three groups of schemes, namely, BLS network to train data, the combination of BLS networks with individual features to train data, and the combination of BLS networks with three features to train data.

BLS Network Training
Results. Adjust the number of feature nodes and enhancement nodes in BLS network, and the training accuracy, test accuracy, and training time are shown in Table 1.
As shown in the first row of Table 1, the numbers of feature nodes and enhancement nodes were set as 200 and 2000, respectively, and the training accuracy was 99.57% after the training set data was input into the BLS network, and the training time was 0.5787 s; the test accuracy is 66.96% after the test set data was input into BLS network.
In Table 1, the number of feature nodes was set to 200. With the increase of the number of enhanced nodes, the test accuracy is improved by a small margin. e number of enhancement nodes was fixed at 10000. With the increase of the number of feature nodes, the accuracy of the test set decreases by a small margin.
is indicates that the enhancement layer of BLS greatly influences the training results. Since the activation function of BLS enhancement layer is tan sig(·), it has the ability to fit nonlinear functions, which is also the core function of BLS network. e training accuracy of BLS network is close to 100%, but the test accuracy is low, which is a typical overfitting phenomenon.
ere are two possible reasons for the overfitting phenomenon in BLS network: (1) e mapping layer is a direct linear transformation system, whose extracted features are not sufficiently expressive or have few extracted features (2) ere are too many parameters in the enhancement layer, which overfits the training set, resulting in poor performance of the model on the test set

Combination of Features.
e main reasons for the poor classification results of two-phase flow by single BLS network are the limitations of network architecture and insufficient feature extraction ability. In this section, the additional extracted features, including P 1 (the number of regional blocks), P 2 (the circularity of the regional blocks), and P 3 (the barycenter coordinates of the regional blocks), are combined with the "feature nodes" of the BLS mapping layer, respectively. e classification results of the two-phase flow by this network are shown in Table 2.
As shown in the first row of Table 2, 200 feature nodes and 2000 enhancement nodes were set, respectively. When feature nodes of BLS and P 1 are combined, the accuracy rate of test set is 76.73% and the training time is 0.6187 s. When P 2 is combined with feature nodes, the accuracy of test set is 97.13% and the training time is 0.6754 s. When P 3 is combined with feature nodes, the accuracy of test set is 97.13% and the training time is 0.6754 s.
Combine the three additional features P 1 , P 2 , and P 3 into P � P 1 P 2 P 3 ; three additional features, P 1 , P 2 , and P 3 , are combined into P � P 1 P 2 P 3 , and P is combined with feature nodes of BLS mapping layer. e identification results of the two-phase flow data set are shown in Table 3.
As shown in the first row of  into BLS network; the test accuracy is 66.54% when test data was input into BLS network. e results of Table 2 show that P 1 improves BLS network to some extent, but the effect is not significant; P 2 and P 3 have significant improvement effect on BLS network and high accuracy of test set. e results of Table 3 show that the test accuracy is close to 100% when P 1 , P 2 , and P 3 jointly affect the BLS network. e training time of the BLS network

Summary and Prospect
In order to better analyze the distribution of blockage material in filling pipeline comprehensively, this paper proposes to adopt ECT method to observe the blockage area visually. According to the permittivity of the plugging material, six shapes of fluids (single-bubble, double-bubble, three-bubble, four-bubble, stratified, and annular flow) in the cross-sectional view are simulated. A fusion method of artificial features and BLS was proposed to recognize the flow pattern of blocked pipes in the gray image reconstructed by ECT. e training set and test set have 3000 and 500 sample images, respectively. In the simulation, the BLS network that can be quickly trained is used to identify the flow pattern of two-phase flow images. Due to the lack of feature extraction ability of BLS network, the simulation results show that the identification effect of test set is poor. In view of this situation, this paper used the method of nonneural network to extract three additional features, which are the number of regional blocks, the circularity of regional blocks, and the barycenter coordinates of regional blocks. ree different features are combined with feature nodes of BLS mapping layer to realize high accuracy test set recognition. In addition, the training time is almost unaffected because of the small amount of data of the additional features selected. Further work is to extract fewer and more prominent features from ECT reconstructed images, thereby improving the accuracy and speed of recognition.
Data Availability e data used in this paper were randomly generated by programming and simulation software, including the training set and the test set (the reconstructed images of the ECT system were taken as samples).

Conflicts of Interest
e authors declare that they have no conflicts of interest.