Next Article in Journal
Synthesizing Rolling Bearing Fault Samples in New Conditions: A Framework Based on a Modified CGAN
Next Article in Special Issue
High-Precision Real-Time Detection of Blast Furnace Stockline Based on High-Dimensional Spatial Characteristics
Previous Article in Journal
Predicting Activity Duration in Smart Sensing Environments Using Synthetic Data and Partial Least Squares Regression: The Case of Dementia Patients
Previous Article in Special Issue
Metal Surface Defect Detection Method Based on TE01 Mode Microwave
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Real-Time Detection and Short-Term Prediction of Blast Furnace Burden Level Based on Space-Time Fusion Features

School of Automation, Central South University, Changsha 410083, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(14), 5412; https://doi.org/10.3390/s22145412
Submission received: 23 June 2022 / Revised: 16 July 2022 / Accepted: 17 July 2022 / Published: 20 July 2022
(This article belongs to the Special Issue Intelligent Sensors for Industrial Process Monitoring)

Abstract

:
Real-time, continuous and accurate blast furnace burden level information is of great significance for controlling the charging process, ensuring a smooth operation of a blast furnace, reducing energy consumption and emissions and improving blast furnace output. However, the burden level information measured by conventional mechanical stock rods and radar probes exhibit problems of weak anti-interference ability, large fluctuations in accuracy, poor stability and discontinuity. Therefore, a space-time fusion prediction and detection method of burden level based on a long-term focus memory network (LFMN) and an efficient structure self-tuning RBF neural network (ESST-RBFNN) is proposed. First, the space dimensional features are extracted by the space regression model based on radar data. Then, the LFMN is designed to predict the burden level and extract the time dimensional features. Finally, the ESST-RBFNN based on a proposed fast eigenvector space clustering algorithm (ESC) is constructed to obtain reliable and continuous burden level information with high accuracy. Both the simulation results and industrial verification indicate that the proposed method can provide real-time and continuous burden level information in real-time, which has great practical value for industrial production.

1. Introduction

The iron and steel industry is a major energy consumer and carbon emitter. As the most critical equipment for iron and steel smelting, blast furnaces account for about 60%~70% of the comprehensive energy consumption and cost of iron and steel [1], and CO2 emissions account for more than 90% [2]. Blast furnace burden level is a key parameter of blast furnace production [3]. Real-time and continuous blast furnace burden level with high accuracy is of great significance for controlling the charging process, ensuring smooth operation of blast furnace [4], reducing energy consumption and emissions and improving the blast furnace output [5]. With the increasing requirements for refinement and intelligence of the smelting process, real-time continuous high-precision measurement of blast furnace burden levels is becoming increasingly important [6].
A blast furnace is an enclosed large countercurrent reactor [7], which has high temperatures, high pressure and much dust, and its environment is extremely harsh [8]. At present, the detection methods for the burden level of blast furnaces are mainly divided into contact methods and non-contact methods [9]. The mechanical stock rod (mechanical probe) is the main contact detection device, which is composed of a steel wire rope, a heavy hammer, a winch motor outside the furnace, a reducer, an encoder and a master controller for detection control [10]. This detection device has a simple structure and reliable mechanical transmission mechanism and can adapt to the complex and harsh environment in the blast furnace. Therefore, the burden level information obtained by a mechanical probe has high accuracy [11]. However, it can neither achieve continuous measurement to ensure the continuity of burden level information nor track the change in burden level in time when the blast furnace shows slip, inspect, or collapse conditions [12]. A radar probe is the most widely used non-contact detection device, which consists of an antenna transmitter, a receiver, and a display screen [13]. The electromagnetic wave signal is firstly transmitted to the burden surface, and then it is reflected by the burden surface and received by the receiver to form a radar echo signal to measure the burden level. Its working principle determines that the radar probe can obtain real-time and continuous burden level information. However, due to the high-speed gas flow and high concentration of dust in the furnace, it is easy to interfere with the echo signal of the radar probe, which distorts the measurement data and reduces the measurement accuracy and reliability [14]. Therefore, obtaining high-precision and real-time continuous measurement data of the burden level is challenging to solve.
Z.H. Ma et al. proposed a closed-loop speed control method, which realizes the anti-overturning control in the variable torque control by changing the anti-overturning torque, so as to make the mechanical probe drop at a uniform speed and ensure its stable operation. On the basis of the mechanical probe tracking process, the descending speed of the burden surface is used to judge whether an abnormal condition occurs [15]. Based on the mechanical probe data, J. Hinnel and H. Saxén established the relationship between the burden thickness and the detection data by using the neural network based on mechanical probe data and systematically analyzed the falling rate of the furnace burden [16]. Although this research has improved the accuracy of measuring the burden level with a mechanical probe, it still cannot overcome the problem of discontinuity in the measuring process. For the research of measuring burden level with a radar probe, in order to solve the problem of low accuracy and reliability of a radar probe, S. Schuster et al. used an integrated millimeter-wave design circuit and advanced signal processing technology to obtain the average burden level information of a blast furnace in real-time [17]. Dominik Zankl et al. proposed a BLASTDAR radar sensing system for real-time radar imaging of blast furnace burden surface [18]. X.Z. Chen et al. proposed a new forward-looking staring high-resolution imaging method for synthetic aperture radar (SAR) to further improve the definition and resolution of radar imaging [19]. The above methods have made a breakthrough in the real-time measurement of burden level with radar probe, especially for the radar imaging of the whole burden level, which is of great significance to blast furnace regulation and smelting. However, there is still a lack of efficient and feasible technical methods to overcome the defects of low measurement accuracy and large fluctuations of radar probes. In recent years, with the rapid development of deep learning technology, X.P. Liu et al. proposed a new codec structure composed of a convolutional neural network and long-term and a short-term memory network to classify distorted signals and perform regression in a learning manner to model long-term historical measurements [5]. Compared with the traditional manual denoising method, this method saves time and energy, but its accuracy still requires improvements. Q.D. Shi et al. combined high-temperature metallurgy, radar detection and image processing, and a new blast furnace surface deep learning detection method of a blast furnace smelting state visualization system for a burden surface based on energy weight is proposed, which realizes the visualization of blast furnace smelting state and digitization of burden surface information [20]. H. Wang et al. proposed a key point estimation method based on learning combined with a key point-based connected region noise reduction algorithm (KP-CRNR) to reconstruct the key points in the BSP image measured by the radar probe. This method improves the measurement accuracy of the radar probe from the perspective of the working principle of the sensor [21]. From the above research, the introduction of deep learning technology improves the accuracy of radar probes for burden level measurement to a certain extent. However, most of the existing studies ignore the important significance of a mechanical probe for burden level measurement [22]. It is difficult to fundamentally solve the current issue of low accuracy and large fluctuation of radar probe level measurement data just from the characteristics of radar data.
Therefore, the blast furnace smelting burden distribution mechanism is further analyzed, and it is discovered that the burden level has the features of slowly changing in the time dimension and rapidly changing in the space dimension. That is, the burden level fluctuates rapidly in the space dimension due to the rapid change in smelting conditions and environment in the furnace, while its periodic features, alternating frequency and other time dimensional features do not follow the rapid change in dimensional space features. Based on this, a prediction and detection method of burden level is proposed based on a long-term focus memory network and an efficient structure self-tuning RBF neural network. Figure 1 shows the research framework. Firstly, the space feature regression model of radar data based on mechanical probe data is established to extract the space dimensional features of blast furnace burden level. Then, the mechanical stock rod data are set as the standard long-term memory, and a long-term focus memory depth network structure is constructed. Combining the historical periodic radar data and the mechanical probe data, not only the burden level information of the next cycle is predicted with high precision but also the time dimensional features of which are accurately extracted. Finally, based on a proposed fast eigenvector space clustering algorithm, an efficient structure self-tuning RBF neural network is constructed. By integrating the advantages of mechanical probe data and radar data, the real-time continuous detection of burden level information with high precision is realized. Simulation results and industrial data verification demonstrate that the proposed method for blast furnace burden level prediction and continuity-accuracy detection has good real-time performance, high precision and obvious industrial practical value.
The remainder of the paper is organized as follows. In Section 2, the space dimensional feature extraction based on the space feature regression model of radar data is introduced. In Section 3, the long-term focus memory network is developed to realize the burden level prediction and tim-dimensional feature extraction. In Section 4, the efficient structure self-tuning RBF neural network based on a fast eigenvector space clustering algorithm is presented to realize real-time, continuous and accurate burden level detection. The simulation and industrial verification results are given in Section 5. Section 6 gives the discussion and the conclusions.

2. Space Dimensional Feature Extraction

2.1. Analysis of the Mechanical Probe Data and Radar Data

The distribution mechanism of a blast furnace is batch charging, which has the nature of periodicity. A distribution cycle of a blast furnace is 3 to 6 min. Generally, the first 1.5 min is the feeding period in a cycle, and the last 1.5 to 6 min is the idle period in a cycle. As shown in Figure 2, the change in blast furnace burden level has the periodic subsection nonlinear characteristic. At the same time, it can be seen that in a distribution cycle, the burden level decreases with the increase in burden surfaces in the feeding period and increases with the decrease in burden surfaces in the idle period. Analyzing the mechanical probe data and radar data indicates that the burden level data sampled by the mechanical probe is discrete, and the mechanical probe is randomly sampled once in each distribution cycle, while the burden level data sampled by the radar probe is continuous. The different sampling properties of the two kinds of probe data make it difficult for the two detection data to match accurately within an infinite time distribution, so it is difficult to integrate the advantages of the two detection data, which makes it very challenging to obtain the blast furnace burden level in real-time with high accuracy. In order to accurately match the data of the two probes and build a strong correlation between the two detection data, a sliding window with a width of N is used to construct a single mapping relationship between the data of the radar probe and the mechanical probe, and Figure 3 shows the schematic of the data processing process of the sliding window. In this way, the matching problem of the two infinite dimensional detection data is transformed into a corresponding problem of the two detection data, providing a data foundation for the subsequent algorithm in this paper. According to the previous research, when the window size is N = 10 , the accuracy of burden level detection is the highest [22].

2.2. Space Feature Regression Model of the Radar Data

The change in the burden level has the space features of material falling speed, acceleration, and jerk, etc. Since the radar data can fully reflect the space change in the burden level, based on the mechanical probe data, a piecewise space feature regression model is established for the radar data to obtain the space feature parameters of blast furnace burden level change. We choose the fifth order curve to fit it, with each coefficient order corresponding to some physical significance. The time sequence of the radar data in a single period is set as Q i = { ( t i 1 , l i 1 ) , ( t i 2 , l i 2 ) , ( t i m , l i m ) } . Piecewise nonlinear regression fitting of the radar data by the fifth-order polynomial is described as
l i = s d i + v d i t + a d i t 2 + j d i t 3 + h d i t 4 + k d i t 5 + ε i ,   t max i t t min i s u i + v u i t + a u i t 2 + j u i t 3 + h u i t 4 + k u i t 5 + ε i ,   t min i t t max i + 1
where i = 1 , 2 , , n presents the i th cycle; u and d present the feeding period and the idle period, respectively; s , v , a , j , h and k are the initial position, falling speed, acceleration, jerk, jerk changing rate and high-order term of the burden level, respectively; t max i and t max i + 1 are the beginning time and the ending time of the cycle; t min i is the dividing time between the feeding period and the idle period of this cycle; ε i is the measurement error.
To improve the accuracy and reliability of the data, combining the measurement data sequence of the mechanical probe t R i , l R i as the standard reference point into Equation (1), the space feature regression model of the radar data based on the mechanical probe data is described as
l i = s d i + v d i t + a d i t 2 + j d i t 3 + h d i t 4 + k d i t 5 + ε i + L R l R i ,   t ' max i t t min i + 1 s u i + v u i t + a u i t 2 + j u i t 3 + h u i t 4 + k u i t 5 + ε i + L R l R i ,   t min i + 1 t t ' max i + 1
where L R is the radar data closest to the detection time t = t i of the mechanical probe in the i th cycle. To ensure that the continuous regression curve of radar data segments is continuous at the periodic junction, the dividing points of the cycle t max i and t max i + 1 divided based on discrete data are adjusted to t ' max i and t ' max i + 1 . Setting l d i t and l u i t as the burden level value in the falling period and rising period of the i th cycle, for l d i t max i ' in the feeding period, l d i t max i ' = l u i 1 t max i + 1 ' , and for l u i t max i + 1 ' in the idle period, l u i t max i + 1 ' = l d i + 1 t max i ' . Figure 4 shows the original burden level measurement curve obtained based on the space feature regression model of radar data and the piecewise regression nonlinear (PRNN) curve after fine-tuning the periodic dividing points. The area between the red solid lines in the figure represents the i th cycle, and the area between the red dotted lines represents the idle period of the i th cycle and the feeding period of the i + 1 th cycle. Compared with the original radar data curve, the space feature regression curve of radar data is smoother and more coherent.

3. Time Dimensional Feature Extraction and Burden Level Prediction Based on Long-Term Focus Memory Network

The time dimensional change in blast furnace burden level corresponds to its space change features, and it has the time dimensional features of extreme time points such as speed, acceleration, jerk and jerk changing rate. Considering the complex and changeable environment in the blast furnace, the radar data are different in different environments. However, in the adjacent cycle, due to the extremely short sampling time interval and the small change in the detection environment, the feature changes in the radar detection data tend to be consistent in the time dimension. According to the distribution mechanism of the blast furnace, the distribution of the blast furnace does not change significantly in a short time; that is, the periodic feature in the time dimension does not change, so the slow variation feature of radar data reflects the features that the distribution mechanism does not change in a short time, and it can well predict the change trend of the burden level. Therefore, based on the slow variation feature of radar data in the time dimension, the LFMN deep network structure is proposed. By training the LFMN network, the radar data of the i + 1 th cycle is predicted. The biggest advantage of this deep network in burden level prediction is that it introduces the attention mechanism with the mechanical probe data as accurate and reliable, slowly varying burden level information, which improves the reliability and accuracy of burden level information predicted based on the slow variation features of radar data. According to the actual distribution mechanism of the blast furnace, we combine the radar data of the i 1 th and the i th cycle with the predicted radar data of the i + 1 th cycle, and the time dimensional features of the i + 1 can be extracted.
The LFMN network is composed of an input layer, a long-term focus memory layer consisting of multiple memory units, and an output layer. Its structure is shown in Figure 5. The input samples of the network input layer are set as follows,
x i ( k ) = t ' max i , t min i + 1 , s d i , v d i , a d i , j d i , γ i , T i , f i , l k N + 1 i , , l k i , t max i ' < t k i t min i + 1 t min i + 1 , t ' max i + 1 , s u i , v u i , a u i , j u i , γ i , T i , f i , l k N + 1 i , , l k i , t min + 1 < t k i t max i + 1 '
where i represents the cycle of the current input data, i q , and q represent the last one of the detected cycles. For the i th cycle, t max i ' , t min i + 1 is the feeding period, t min i + 1 , t max i + 1 ' is the idle period, γ i , T i , f i are the working condition parameters that have the greatest correlation with the blast furnace burden level, representing the opening of the material flow valve, the furnace top temperature and the furnace top air volume, respectively. l k N + 1 i , , l k i is the N radar data in the k th window of the i th cycle processed by the sliding window, and the sliding windows are shown in Figure 3. The output of the network is set as the predicted n radar data of the q + 1 th cycle l 1 q + 1 , , l n q + 1 .
For each memory unit in the long-term focus memory layer, it contains two state variables h and C , saving the short-term memory and long-term focus memory of the network. At time t , each memory unit has three inputs, which are the input sample x t of the time t , the short-term memory h t 1 of the time t 1 , and the long-term focus memory C t 1 . The internal structure of each long-term focus memory unit is shown in the right figure of Figure 4, including four control gates, which are the input gate, forgetting gate, guiding gate, and output gate. The long-term focus memory C is controlled by the input gate and the forgetting gate, which determines how many input states at the current time t are updated into the long-term focus memory.
The forgetting gate determines how much of the long-term focus memory C t 1 of the time t 1 is retained for the current time t . Therefore, the long-term focus memory C t of the time t can be described as
C t = f t C t 1 + i t C t f t = σ w f h t 1 , x i k + b f i t = σ w i h t 1 , x i k + b i C t = tanh w c h t 1 , x i k + b c
where f t , i t represent the forgetting gate and the input gate; w f , w i are the weight matrix of forgetting gate and input gate in network training; b f , b i are the bias of the forgetting gate and input gate; σ represents the s i g m o i d activation function; C t is the long-term memory of the current time t ; w c , b c represent the weight matrix and the bias of the current inputting long-term state.
The guiding gate is set to retain the focus information of the long-term focus memory unit. This is because the mechanical probe data contains the most important long-term changing information of the burden level, and it has the advantages of high accuracy and reliability. The model of the guiding gate is described as
C t ' = θ C y + 1 θ C t C y = j = 1 i f j + L R j l R j
where C t ' represents the long-term focus memory of the time t stored in the long-term focus memory layer; C y represents the long-term burden level focus information, which combines radar data and mechanical probe data; L R j is the radar data closest to the mechanical probe data in the j th cycle. l R j is the mechanical probe data in the j th cycle; f j represents the burden level prediction model based on radar data of the j th cycle, which can be expressed as
f j = r = 1 5 α j r exp t b j r / c j r 2 a j r b j r c j r = a j 1 a j 2 a j 3 a j 4 a j 5 b j 1 b j 2 b j 3 b j 4 b j 5 c j 1 c j 2 c j 3 c j 4 c j 5
where α j r , b j r , c j r , r = 1 , 2 , , 5 are the radar model’s fifth order Gaussian fitting coefficients of the j th cycle of the radar data.
The output gate is responsible for controlling the contribution of long-term focus memory to the network output value at the current time t . Its effect can be quantitatively described as
h t = o t tanh C t ' o t = σ w o h t 1 , x i k + b o
where o t represents the output gate, w o , b o are the weight matrix and the bias of output gate. The LFMN network can be trained by the error back propagation algorithm, and the predicted n burden level values of the i + 1 th cycle are
h i + 1 = l 1 i + 1 , , l n i + 1
Combining the radar data of the i 1 th and the i th cycle with the predicted radar data of the i + 1 th cycle, the time dimension features of the i + 1 can be extracted through the following Equation (9)
l k i = α l k i 1 + β l k i + ν l k i + 1 3 α + β + ν = 1       ,       k = 1 , 2 , , n
where k represents the predicted burden level information of the k th cycle, and α , β , ν are the data fusion weight of the i 1 , i , i + 1 th cycle, respectively. The weights are selected based on experience, which are set as 0.25 , 0.25 , 0.5 in this paper. Then, for the burden level sequence l k i ' of the i th cycle with time dimensional features of multiple cycles, the space feature regression model of radar data is used for regression to obtain the continuous function l i ' of the burden level change in the cycle, and the time dimensional features of blast furnace burden level of the i th cycle can be extracted, as shown below
c t i d = t i d v min , t i d v max , t i d a min , t i d a max , t i d j min , t i d j max c t i u = t i u v min , t i u v max , t i u a min , t i u a max , t i u j min , t i u j max
where t i d v min , t i d v max , t i d a min , t i d a max , t i d j min , t i d j max and t i u v min , t i u v max , t i u a min , t i u a max , t i u j min , t i u j max are the time points at which the extreme values of burden level speed, acceleration and jerk are taken in the feeding period and the idle period of the i th cycle.

4. ESST-RBFNN

Based on the slow change features in time dimension and fast change features in the space dimension of the burden level, integrating the advantages of the two detection data and overcoming the common shortcomings, the efficient structure self-tuning RBF neural network (ESST-RBFNN) is designed, and the real-time continuous detection of blast furnace burden level with high accuracy can be realized. The structure of the network will change with the change in input, which adapts to different blast furnace environments and have industrial universality. At the same time, an efficient eigenvector space clustering (ESC) algorithm is proposed, which improves the efficiency of training and ensures the real-time detection of the burden level. Firstly, the input samples of the algorithm are composed of the time dimensional features and space dimensional features of the burden level, the radar data processed by the sliding window and the key working condition parameters. Secondly, based on the proposed fast eigenvector space clustering algorithm, the clustering center and width of the input samples are determined. Thirdly, based on the RBF neural network, a dynamic self-tuning network structure is designed. According to the characteristics of input samples, the number of hidden nodes is automatically increased and clipped to realize that the network structure changes with the characteristics of input samples, and the relevant parameters of the network structure are adjusted and updated independently to realize the self-tuning of the structure and complete the learning and training process of the ESST-RBFNN algorithm. Finally, the trained network is used to output the burden level information at the next time, and the whole algorithm is realized.

4.1. Fast Eigenvector Space Clustering Algorithm (ESC)

The input sample consists of four parts, which are the key operating condition parameters, space dimensional features, time dimensional features of each cycle and the radar data processed by the sliding window in this cycle. Its specific form is shown as follows
x i ( k ) = t ' max i , t min i + 1 , c t i d , s d i , v d i , a d i , j d i , γ i , T i , f i , l k N + 1 i , , l k i ,                                                                                                                                           t max i ' < t k i t min i + 1 t min i + 1 , t ' max i + 1 , c t i u , s u i , v u i , a u i , j u i , γ i , T i , f i , l k N + 1 i , , l k i ,                                                                                                                                           t min + 1 < t k i t max i + 1 '
where x i k represents the input sample of the k th window in the i th cycle. It can be analyzed that radar data are the most critical data in the input samples. In a single cycle, burden level changes are usually up and down the standard burden level, which leads to the typical accumulation of radar data. In addition, the environment in the furnace is bad, meaning that the echo-based radar data fluctuate greatly, so the accuracy is not high, and the signal-to-noise ratio is low. To overcome the interference of the furnace environment, the fast eigenvector space clustering algorithm (ESC) is proposed to reveal the intrinsic correlation of the input data. The algorithm clusters n input sample datasets X ^ = x 1 , x 2 , x k , x n into K categories. The clustering center vector set is C = c 1 , c 2 , c K , and the input sample data set is centralized as X = ( x 1 , x 2 , , x n ) T , where x i d , i = 1 , 2 , , n . In the d -dimensional data space, the problem of making the clustering center vector describe the internal clustering characteristics of input samples as much as possible can be equivalent to the optimization problem, as shown in the following Equation (12)
min J K = i n x i 2 k = 1 K 1 n k X i , c K C x i c K T
where n K is the number of input samples in category K . Set the original indicator vector matrix of the category feature space as G K = ( g 1 , g 2 , , g K ) T ,where each row of the original indicator vector matrix of the category feature space has only one non-zero value, and each column can have multiple non-zero values, which indicates that each datum only belongs to the category represented by the column and does not belong to other categories. The matrix H K = ( h 1 , h 2 , , h K ) T is obtained by the rotation equivalent transformation of the matrix G K , where h K = ( 0 , , 0 , 1 , , 1 n K , 0 , , 0 ) T / n K , and the non-zero value in h K represents the corresponding sample belonging to the category K in the category feature space. After introducing the feature space indicator vector, the optimization objective function can be converted into the following Equation (13).
min J K = i n x i 2 T r ( H K X X T H K T )
It should be noted that among the K clustering categories in the sample space X , the mean center of the sample has an impact on the embodiment of the intrinsic clustering characteristics of the samples. In order to eliminate this impact, the input sample space is centered into X ˜ = ( x ˜ 1 , x ˜ 2 , , x ˜ n ) T , where x ˜ i = x i x ¯ is the centralized data matrix obtained by removing the average value of the given data set x ¯ = i n x i / n . The optimization objective function is rewritten as the following Equation (14)
min J K = i n x ˜ i + x ¯ 2 T r H K X ˜ + X ¯ X ˜ + X ˜ T H K T
To further simplify the optimization objective function, a K × K orthogonal transformation matrix T is constructed, and the last column of the matrix is set as
T K = ( n 1 / n , n 2 / n , , n K / n ) T
Then the transformation matrix is linearly transformed into the indicator vector matrix of the category feature space H K , and the transformed matrix is Q K ,
Q K = H K T T
The last column of the matrix Q K is given as
q k = ( n 1 / n h 1 + n 2 / n h 2 + + n K / n h K ) T             = 1 / n e T     ,             e T = ( 1 , , 1 , , 1 ) T
By analyzing the relevant properties of the indicator vector matrix H K in the feature space, it is found that its vector h K has the property of pairwise orthogonality, and it can be expressed in the following form,
h K T h L = δ K L δ K L           =       1         ,         K = L δ K L             =     0       ,         o t h e r w i s e
For the matrix Q K , the following Equation (19) holds
q K T q L = T K T h K T h L T L = δ K L
Set Q K 1 = ( q 1 , q 2 , , q K 1 ) , then using Equation (19) derives the following Equation (20)
Q K 1 T Q K 1 = I K 1
Substituting Equations (16) and (17) together, the rewritten matrix Q K = [ Q K 1 , q K ] into the optimization objective Equation (14) for a series of derivation, transformation and simplification, the following Equation (21) can be obtained,
J K = T r X ˜ X ˜ T T r Q K 1 T X ˜ X ˜ T Q K 1
Considering that T r X ˜ X ˜ T is a constant larger than zero and it is determined by the input sample space, so minimizing J K and maximizing J D = T r Q K 1 T X ˜ X ˜ T Q K 1 is equivalently.
Using K y F a n theorem, for a real symmetric matrix A = X ˜ X ˜ T with eigenvalues, λ 1 λ 2 λ n are its eigenvalues and ν 1 , ν 2 , , v n are its eigenvectors. Under the constraint of Q K 1 T Q K 1 = I K 1 , the maximum value of J D is Q K 1 = ν 1 , ν 2 , , ν K 1 U , where U is a K × K orthogonal matrix, and the maximum value is
m a x J D = λ 1 + λ 2 + + λ K 1
In combination with Equation (21), the upper and lower boundary of the optimization objective can be obtained after proper and reasonable expansion as
n X ˜ 2 ¯ k = 1 K λ k < J K < n X ˜ 2 ¯
where n X ˜ 2 ¯ is the total variance of samples X .
Because the indicator vector matrices H K 1 and Q K 1 of the category feature space have the following relationship
Q K 1 = H K 1 T T K 1
where T K 1 is a K 1 × K 1 orthogonal matrix. Then the indicator vector matrix K   H K 1 of the category feature space can be
H K 1 T = ν 1 , ν 2 , , v K 1
In this way, the optimization objective function is guaranteed to obtain the minimum value. Then, according to the indicator vector matrix H K 1 , for any one of the former K 1 clustering categories in which the input samples X are divided into K categories, the corresponding indicator vector of category feature space is H j = v j 1 , v j 2 , , v j n , which can be used to judge whether the n samples belong to the category j . The discriminant is as follows,
I C j = x i v j i > 0 i = 1 , 2 , , n         j = 1 , 2 , K 1    
where I C j represents the input sample set belonging to category j . Note that the distribution of the former K 1 categories of the input sample set X ^ can be determined by using Equation (26). For the category, it is common practice to classify the samples that do not belong to the former K 1 categories in the input sample set X ^ as the category K . The mathematical description is as the following
I C K = X ^ j = 1 K 1 I C j
For each clustering center vector, it can be calculated by the following Equation (28)
c j = x s I ( C j ) x s / m j = 1 , 2 , , K
where m is the elements number of the set I C j .

4.2. Efficient Structure Self-Tuning RBF Neural Network (ESST-RBFNN)

Due to the strong nonlinear correlation between the radar data and the mechanical probe data, in order to effectively integrate the advantages of the two-detection data and improve the burden level measurement accuracy, a network named ESST-RBFNN with a flexible and self-tuning structure is presented. The network structure is shown in Figure 6, and the mathematical model of the network is described as
g j x i k = exp x c j 2 σ j 2   ,   j = 1 , 2 , , K f x i k = y i k = j = 1 m w j g j x i k
where x i k is the k th fusion input sample of the i th cycle; y i k is the k th output burden level data of the i th cycle; g j x i k is gaussian basis excitation function of the j th hidden node; c j is the center of the j th basis function; σ j is the width of the j th basis function; K is the number of hidden nodes; w j is the output weight. The number of hidden nodes, the center and the width of the basis function can be determined by the clustering results of the ESC algorithm for the input samples set X ^ . Generally, the number of hidden nodes is the number of clustering classes, and the centers of the basis functions are the clustering centers. The width of the basis function σ j is determined by the following Equation (30)
σ j = 2 d max j 3               ( 1 j K ) G i = x i k I ( C j ) d max j = max x i k G i x i k c j 2
where d max j is the maximum distance from the clustering center to the sample point.
To determine the weight w j , the output matrix of the hidden layer H can be calculated according to Equation (29), which is a n × K matrix. The mechanical probe data corresponding to the input samples used for neural network training can be set as a matrix Y , and the network weight vector W can be obtained by solving the generalized inverse as the following
W = [ H T H ] 1 H T Y
According to the network structure shown in Figure 6, the network has the characteristic of dynamic structure self-tuning, which can automatically increase the number of hidden layer nodes, change the network structure and tune the network parameters itself with the increase in training samples. The main steps of ESST-RBFNN are as follows,
Step 1: According to the characteristics of input samples, set the value of category discrimination radius b as a positive real number.
Step 2: The algorithm begins. When the radar data of the first sliding window of the first cycle is sampled, the first input sample is constructed as x 1 1 according to Equation (11), and the first cluster is formed and K = 1 , the clustering center C = c 1 , c 1 = x 1 1 , and the first hidden node of the network hidden layer h 1 is generated.
Step 3: When the k th input sample data of the i th cycle is sampled, the number of samples n = k , k = 2 , 3 , . There have been D 1 clustering categories. The number of classes is K = D 1 . The clustering centers C = c 1 , c 2 , c D 1 , and the hidden nodes are h p , p = 1 , 2 , 3 , , K . The input sample data matrix is X = x i 1 , x i 2 , , x i n T . Set ς = ( x i s , y i s ) x i s d , y i s l , s = 1 , 2 , , k 1 as the network training set, where d and l are the training input sample set and its corresponding mechanical probe data set. If m i n x i k c j c j C 2 < b , keep the current clustering result and the network structure unchanged, and go to Step 7. Otherwise, dynamically adjust the network structure, increase the number of the hidden layer nodes by 1, go to Step 4, update the clustering result, and tune the network structure itself.
Step 4: Centralizing the input sample data matrix X into X ˜ = ( X ˜ 1 , X ˜ 2 , , X ˜ n ) T , and calculate the former D eigenvectors v 1 , , v D of X ˜ X ˜ T to obtain the feature space indicator vector of the category to which each sample belongs. All samples are reclassified by using Equations (26) and (27) to obtain a new set of categories I = I C 1 , I C 2 , , I C D .
Step 5: Use Equation (28) to calculate and update the clustering center vector set and combine Equation (30) to synchronously update the center and width of the basis function of the hidden nodes in the hidden layer of the network.
Step 6: Use Equation (31) to calculate and update the network weight vector to complete the establishment and learning of ESST-RBFNN.
Step 7: Use Equation (29), set x i k as the model input, and the current burden level information y i k can be obtained through real-time calculation. With the loop going on, the real-time burden level measurement can be realized with high precision.

5. The Simulation and Industrial Verification Results

For the training of the network model, the hardware environment included Intel Xeon (R) w-2133CPU, NVIDIA Geforce, and one RTX2080 graphics card; the language we chose is python, and the platform we used is pytorch 1.7.1 and the CUDA version was 10.2. All our simulation diagrams were drawn based on MATLAB. To verify the correctness of the proposed methods of prediction and measurement of blast furnace burden level, the actual detection data of mechanical probe and radar probe in the same area of a steel plant blast furnace were collected for simulation. The radar probe sampled every 10   s , and the number of samples is 76,426 rows, where 66,426 rows of the data were used for model training and another 10,000 rows were used for model validation. The width of the sliding window was preferred as N = 10 . The data format is 76,424 × 10. The sampling interval of mechanical stock rod is 3 6 min, and the number of samples is 3149.

5.1. Verification Results of the Blast Furnaces Burden Level Prediction and the Time Features Attraction Based on LFMN

LFMN is proposed to predict the blast furnace level in order to extract the time dimensional features of extreme time points such as speed, acceleration, jerk, and the change rate of jerk. For this purpose, it is required that the predicted burden level information can reflect the periodicity of burden level change and fit the real-time mechanical probe measurement data as much as possible, with a certain high accuracy. Figure 7 shows the trend of the prediction of burden under normal working conditions. It can be seen that the burden level data measured by the mechanical probe are mostly above and below the standard control stockline 1.5     m of the blast furnace, and the time of burden level data in a single cycle measured by the radar probe is kept within the 500     s left and right range. The reason for the radar data being much higher or lower than the mechanical probe data is that the radar uses microwaves to detect the burden level, and the microwaves are easily affected by the high-speed gas flow and dust in the blast furnace. However, due to the introduced guide layer mechanism and the slow variation features of the radar data, although the burden level predicted by LFMN has a certain deviation from the measured burden level data of the mechanical stock rod, it is highly consistent with the radar probe data in the change trend, especially in the fluctuation cycle of the data. This shows that the method of LFMN has high confidence in predicting blast furnace burden level information and extracting its time dimensional features for normal working conditions of the blast furnace.
Figure 8 shows that when the blast furnace is abnormal, the measured value of the mechanical probe usually deviates greatly from the standard control 1.5     m stockline of the burden, or even reaches the limit position 2.1     m , as shown in the figure, while the measured data of the radar probe often jumps greatly and frequently. From the prediction result of LFMN, under abnormal conditions, the accuracy of LFMN for burden level prediction decreases significantly for the reason that the network relies too much on radar data. In the case of an extremely bad blast furnace environment, the burden level changes periodically. The radar probe is affected by extreme environments, and the predicted burden level deviation is large. For example, at the time of 3600     s , the error between the predicted value of burden level and the measured value of the mechanical probe even exceeds 0.2     m . However, it can be found that due to abnormal furnace conditions, the blast furnace operator has significantly adjusted the blast furnace charging cycle into the period of 6030     s 7030     s . The charging cycle in this period is much larger than 6   min in this cycle. After this cycle, the operator artificially shortened the blast furnace charging cycle to maintain the burden level at the standard control stockline. The figure shows that the length of charging cycle is roughly in the range of 300     s 480     s after the time of 7210     s . At this time, the prediction method of LFMN can perfectly track the periodic change in the charging cycle and maintain considerable accuracy. In conclusion, it shows that the prediction method of LFMN proposed in this paper can not only effectively reflect the change trend of the burden level but also accurately predict both the periodicity and time change trend of the burden level under normal or abnormal conditions. In order to further verify the statistical significance and effectiveness of the proposed method of LFMN, a comparative statistical analysis on the consistency between the accuracy of the burden level predicted by the proposed method and the time period and the radar probe is conducted. The conclusions are shown in Table 1 below. From the statistical table, it can be found that the LFMN predicted error of the burden level is close to the radar probe measured error. The absolute error and relative error are 1.392 % and 0.06 , respectively. This shows that the method of LFMN can be used to directly predict the burden level of the blast furnace, which has a good reference value for the actual operation of the blast furnace. In addition, from the point of view of the coincidence degree of the blast furnace charging cycle, the method of LFMN reached the coincidence degree of 98.87 % , which shows that the proposed method has high reliability in extracting time dimensional features of the blast furnace burden level information.

5.2. Verification Results of the Blast Furnace Burden Level Detection Based on ESST-RBFNN

ESST-RBFNN is proposed to realize real-time, accurate, continuous, and reliable detection of burden level information. Figure 9 shows the comparison of different methods for measuring burden levels under normal working conditions. According to the discrete burden level data measured by the mechanical probe, the burden level fluctuates around the standard 1.5     m stockline, and the length of the charging cycle is basically fixed. Although the radar data are continuous in real-time, it deviates from the accurate and reliable mechanical probe data at different time points to different degrees, indicating that its accuracy still needs to be greatly improved. The blast furnace burden level data predicted by traditional RBFNN can fit the periodicity and the change trend of the radar data to a certain extent, but it has a strong inhibition on the fluctuation of the burden level, meaning that the fluctuation range of the measured burden level is within 0.2     m , which obviously does not conform to the actual blast furnace charging site. At the same time, it cannot accurately track the changes in the mechanical probe data, and the overall detection level is high. This is because the traditional RBFNN does not take the correlation between the data measured by the mechanical probe and the radar probe into account. It only considers the radar data but ignores the accuracy and reliability of the mechanical probe data, which makes the accuracy and reliability of the detected data low and does not conform to the real production process. The ESST-RBFNN extracts accurate time dimensional features and space dimensional features of the burden level through the method of LFMN and the space feature regression model based on radar data. It also deeply excavates and utilizes the internal relationship between the two detection data. The detected burden level data can not only perfectly fit the periodicity and change trend of radar probe data, but also track the changes in mechanical probe data in real-time with high precision. It greatly improves the accuracy and reliability of the detection data, which shows that the burden level data detected by this method not only has the real-time and continuous advantages of radar data but also has the accurate and reliable advantages of mechanical probe data. This shows its high practical value.
Figure 10 shows the comparison between the radar data, the detection data of the traditional RBFNN and the ESST-RBFNN under abnormal working conditions. Under this abnormal working condition, the mechanical probe data jumps up and down greatly. It can be seen from Figure 10 that at this time, the radar data can hardly track the mechanical probe data, and its accuracy is very low. This shows that under abnormal working conditions, due to a large amount of dust in the blast furnace, the reference value of the radar data is limited. The traditional RBFNN also has the disadvantage of obviously restraining the fluctuation of burden level, and its sensitivity is poor. At the same time, similar to the radar data, the detection data of the traditional RBFNN cannot track the mechanical probe data, so its reliability and accuracy are not high. The ESST-RBFNN proposed in this paper is not only perfectly consistent with the periodicity and change trend of the radar data but can also accurately track the mechanical probe data most of the time. Only in the case of extremely accidental anomalies, such as the three outliers shown in the figure, when the burden level suddenly jumps sharply and changes greatly, the tracking effect at the outlier position is poor. However, compared with the radar data, the reliability has been significantly improved.
Figure 11 and Figure 12 show the absolute error and relative error distributions of the burden level data measured by radar and the methods of traditional RBFNN and ESST-RBFNN compared with the measured data of mechanical stock rod. The area between the black dotted lines represents the main distribution area of the data measured by ESST-RBFNN. It can be seen from the figure that the absolute error and relative error of radar data fluctuate greatly. The absolute error is within the range of ± 0.5     m , and the relative error even reaches 40 % . The absolute error and relative error of the traditional RBFNN method are significantly improved compared with the radar probe; however, the conditions of burden level lower than 1.4     m are minimal during the actual operation of the blast furnace, and the traditional RBFNN method relies on the sample data too much, which leads to the error distribution points of the method concentrated in the upper area of the figure, that is, the measured value has an obvious larger error compared with the real value. Therefore, the maximum relative error of the traditional RBFNN method reaches up to 25 % . For the actual blast furnace charging process controlled by the standard 1.5     m stockline, the errors of the above two methods for burden level measurement are unacceptable in engineering, and the practical application value of which is thus low. The proposed ESST-RBFNN overcomes the problem that the traditional RBFNN is excessively dependent on the data samples by using the space-time features of burden level distribution. The absolute error and relative error range shown in the figure are within the range highlighted by the red line. The absolute error is within the range of ± 0.1     m and the relative error is within the range of ± 5 % . At the same time, considering the relative error of the mechanical probe data itself, on the premise of realizing the real-time continuous measurement of the blast furnace burden level, the measured error of the burden level detection method proposed in this paper is very close to the measured error of the mechanical probe, which satisfies the requirements of the blast furnace smelting process and has very high practicality.
Figure 13 shows the 45° line diagram of the measurement effect of the above three methods to reflect the accuracy of the three material level measurement methods more intuitively. It can be seen from the figure that the detection data obtained by the method of ESST-RBFNN are distributed in the red measured standard line, forming a vertebral body, whose accuracy is significantly high. However, compared with the detection methods of radar probe and traditional RBFNN, the method of ESST-RBFNN has the closest detection data to the standard line with the highest accuracy. When the standard line is within 1.5 2.0     m , the data of three measurement methods are more distributed. The radar data deviate far from the standard line with low accuracy. The traditional RBFNN measurement data also deviate relatively far from the standard line and have the inherent error of high detection level, and its accuracy is also low. The detection data of the ESST-RBFNN method are concentrated and closely distributed around the standard line, which can track the standard burden level well, with high accuracy and reliability. When the standard line is higher than 2.0     m , it indicates that the blast furnace has an abnormal working condition. From the figure, it indicates that the detection methods of a radar probe and traditional RBFNN are far from the standard line, and the accuracies of which are much lower. The method of ESST-RBFNN has little fluctuations and can still track the measured burden level with high precision.
To further verify the correctness and effectiveness of the proposed method, the 24-h operation data based on the 10,000 rows of validation data of the ESST-RBFNN method were collected to draw the operation effect diagram, as shown in Figure 14. The absolute error, relative error and hit rate of radar probe, and the methods of traditional RBFNN and ESST-RBFNN for up to 3 months are counted, as shown in Table 2. It can be seen from Figure 14 that based on the 10,000 rows of validation data within 24 h of operation, whether under normal or abnormal conditions, the detection method of ESST-RBFNN can well track the mechanical probe data with high accuracy and reliability. At the same time, it can well fit the periodicity and variation trend of the radar data and has good real-time continuity. From the 3-month statistical error analysis in Table 2, the MRE (Mean Relative Error) and RMSE (Root Mean Square Error) of the burden level detection method of ESST-RBFNN are nearly four times higher than those of the radar probe and the method of traditional RBFNN. The hit rate of samples with relative error lower than 2 % reaches up to 91.17 % . The hit rate of samples with relative error lower than 5 % even reaches up to 99.33 % . Compared with the radar probe and the method of traditional RBFNN, there is a significant qualitative improvement. The statistical analysis data further illustrate that the burden level detection method of ESST-RBFNN proposed in this paper can obtain continuous and accurate blast furnace burden level information with high precision.

6. Conclusions and Discussion

6.1. Conclusions

Real-time and continuous blast furnace burden level detection with high accuracy and high reliability is of great significance to realizing green and fine production of blast furnaces. The traditional mechanical probe detection method and radar detection method cannot obtain real-time accurate and continuous burden level information. The proposed method combines the advantages of high precision and reliability of mechanical probe data and the real-time continuous characteristics of radar data. Based on the slow variation features of burden level change in the time dimension and fast variation features in the space dimension, it can obtain real-time accurate and continuous burden level information with high precision and great reliability.
Based on the slow variation features of radar data, the LFMN network can predict the blast furnace burden level with relatively high accuracy and extract time dimensional features of the burden level information with great reliability. Combined with the proposed ESC algorithm, the ESST-RBFNN can realize a self-renewing dynamic structure with the ability of online training and provide a better approximation performance than the traditional RBFNN. The simulation results and industrial verification indicate that the proposed method can realize real-time and continuous blast furnace burden level detection with high accuracy, high stability, and great reliability. Therefore, the proposed method has a high industrial application value.

6.2. Discussion

The proposed method can track the burden level information well under normal working conditions and most abnormal working conditions, which obtains real-time and continuous burden level information with high precision and great reliability. However, under some extremely abnormal working conditions, when some extremely large fluctuations occur, the accuracy and reliability of burden level detection will decrease to some extent, where its practical value can be improved by further optimization. The main reason is that the LFMN network relies too much on the slow variation features of radar data, which means that when the burden level fluctuates significantly, the prediction of burden level based on the slow variation features of radar data will not track the real information of the burden level change trend, which needs to be further studied and improved in the future.

Author Contributions

Methodology, W.G.; Supervision, C.Y.; Writing—original draft, Y.C.; Writing—review & editing, Z.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the General Project of Hunan Natural Science Foundation of China under Grant 2022JJ30748, the National Natural Science Foundation for Young Scholars of China under Grant 61903325, the Foundation for Innovative Research Groups of the National Natural Science Foundation of China under Grant 61621062, the National Major Scientific Research Equipment of China under Grant 61927803.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data available on request due to privacy. The data presented in this study are available on request from the corresponding author. The data are not publicly available for the reason that it is real-time industrial operation data. In order to be responsible for the industry, it is difficult for us to disclose these data, but you can obtain them from the corresponding author by email.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Nomenclature

Abbreviations
PRNNpiecewise regression nonlinear
LFMNlong-term focus memory network
ESCeigenvector space clustering
RBFNNradial basis function neural network
ESSTefficient structure self-tuning
MREmean relative error
RMSEroot mean square error
Symbols
Q i time sequence of the radar data in a single period (-)
l i piecewise nonlinear regression fitting of radar data ( m )
s the initial position ( m )
v falling speed of the burden level ( m s 1 )
a falling acceleration of the burden level ( m s 2 )
j falling jerk of the burden level ( m s 3 )
X ^ n input sample datasets (-)
K Categories (-)
X Centralized input data set (-)
n K number of input sample in category K (-)
G K original indicator vector matrix of the category feature space (-)
h K corresponding sample belonging to the category K (-)
T transformation matrix (-)
h jerk changing rate of burden level ( m s 4 )
k high-order term of burden level ( m )
t max i beginning time of the i th cycle ( s )
t max i + 1 ending time of the i th cycle ( s )
t min i dividing time between feeding and idle period of i th cycle ( s )
ε i measurement error ( m )
l d i t burden level value in the falling period of the i th cycle ( m )
l u i t burden level value in the rising period of the i th cycle ( m )
x i ( k ) input samples of LFMN network (-)
γ i opening material flow valve ( m )
T i furnace top temperature ( ° C )
f i furnace top air volume ( c f m )
h state variables saving the short-term memory of LFMN network (-)
C state variables saving the long-term focus memory of LFMN network (-)
w f the weight matrix of forgetting gate in network training (-)
w i the weight matrix of input gate in network training (-)
b f the bias of the forgetting gate (-)
b i the bias of the input gate (-)
I C j input sample set belonging to category j (-)
y i k k th output burden level data of the i th cycle (-)
σ j width of the j th basis function (-)
d max j maximum distance from clustering center to the sample point (-)
ς network training set (-)
d training input sample set (-)
A real symmetric matrix (-)

References

  1. Mikhailova, U.V.; Kalugina, O.B.; Afanasyeva, M.V.; Konovalov, M.V. Development of Automated Control System of Blast-Furnace Melting Operation. In Proceedings of 2019 International Russian Automation Conference (RusAutoCon), Sochi, Russia, 8–14 September 2019; pp. 1–4. [Google Scholar]
  2. Gurin, I.A.; Spirin, N.A.; Lavrov, V.V. MES development for optimal distribution of fuel and energy resources in blast-furnace production. In Proceedings of 2017 International Conference on Industrial Engineering, Applications and Manufacturing (ICIEAM), St. Petersburg, Russia, 16–19 May 2017; pp. 1–6. [Google Scholar]
  3. Zhou, T.; Zhang, Z.-X.; Wu, Z.-X. Experimental Study on Low NOx Emission Using Blast Furnace Gas Reburning and Industrial Application in Stoker Boiler. In Proceedings of 2011 Second International Conference on Digital Manufacturing & Automation, Zhangjiajie, China, 5–7 August 2011; pp. 523–526. [Google Scholar]
  4. Gao, L.; Chen, H. Abnormal detection of blast furnace condition using PCA similarity and spectral clustering. In Proceedings of 2018 13th IEEE Conference on Industrial Electronics and Applications (ICIEA), Wuhan, China, 31 May–2 June 2018; pp. 2198–2203. [Google Scholar]
  5. Liu, X.; Liu, Y.; Zhang, M.; Chen, X.; Li, J. Improving stockline detection of Radar Sensor Array Systems in Blast Furnaces Using a Novel Encoder–Decoder Architecture. Sensors 2019, 19, 3470. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Zhang, R.X.; Cheng, Y.X.; Li, Y.; Zhou, D.D.; Cheng, S.S. Image-based flame detection and combustion analysis for blast furnace raceway. IEEE Trans. Instrum. Meas. 2019, 68, 1120–1131. [Google Scholar] [CrossRef]
  7. He, Z.K.; Yin, Y.H.; Zhang, S. Burden surface fitting and the design of simulation platform for the blast furnace burden distribution. In Proceedings of 2015 IEEE International Conference on Communication Problem-Solving (ICCP), Guilin, China, 16–18 October 2015; pp. 72–77. [Google Scholar]
  8. Ito, M.; Matsuzaki, S.; Kakiuchi, K.; Isobe, M.; Ishida, T.; Fujii, A.; Naito, M. Development of the visual information technique of blast furnace process data. In Proceedings of the 2004 IEEE International Conference on Control Applications, Taipei, Taiwan, 2–4 September 2004; pp. 884–889. [Google Scholar]
  9. Blencowe, M.P.; Buks, E. Quantum analysis of a linear dc SQUID mechanical displacement detector. Phys. Rev. B 2007, 76, 014511. [Google Scholar] [CrossRef] [Green Version]
  10. Fiorina, D. Novel Triple-GEM Mechanical Design for the CMS-ME0 Detector, Preliminary Performance and R&D Results. In Proceedings of 2020 IEEE Nuclear Science Symposium and Medical Imaging Conference (NSS/MIC), Boston, MA, USA, 31 October–7 November 2020; pp. 1–5. [Google Scholar]
  11. Li, Y.J.; Zhang, S.; Zhang, J.; Yin, Y.X.; Xiao, W.D.; Zhang, Z.Q. Data-driven multiobjective optimization for burden surface in blast furnace with feedback compensation. IEEE Trans. Ind. Inform. 2020, 16, 2233–2244. [Google Scholar] [CrossRef] [Green Version]
  12. Zhou, H.; Zhang, H.F.; Yang, C.J. Hybrid-Model-Based Intelligent Optimization of Ironmaking Process. IEEE Trans. Ind. Electron. 2020, 67, 2469–2479. [Google Scholar] [CrossRef]
  13. Du, Z.; Zhang, F.; Zhang, Z.; Yu, W. Radar Detector in Uncoordinated Communication Interference Plus Partially Homogeneous Clutter. IEEE Commun. Lett. 2021, 25, 1999–2003. [Google Scholar] [CrossRef]
  14. Coluccia, A.; Fascista, A.; Ricci, G. A KNN-Based Radar Detector for Coherent Targets in Non-Gaussian Noise. IEEE Signal Processing Lett. 2021, 28, 778–782. [Google Scholar] [CrossRef]
  15. Ma, Z.H. Control method and application of furnace top rod system based on S120 converter. Metall. Ind. Autom. 2015, 39, 68–71. [Google Scholar]
  16. Hinnel, J.; Saxén, H. Neural network model of burden layer formation dynamics in the blast furnace. ISIJ Int. 2001, 41, 142–150. [Google Scholar] [CrossRef] [Green Version]
  17. Schuster, S.; Zankl, D.; Scheiblhofer, S.; Feilmayr, C.; Reisinger, J.; Feger, R.; Stelzer, A.; Schmid, C. Massive MIMO Radar based Burden Surface Imaging: How mm-Wave Integrated Circuits Enable Optimization of Blast Furnace Charging. In Proceedings of 2020 IEEE MTT-S International Conference on Microwaves for Intelligent Mobility (ICMIM), Linz, Austria, 23 November 2020; pp. 1–4. [Google Scholar]
  18. Zankl, D.; Schuster, D.; Feger, R.; Stelzer, A.; Scheiblhofer, S.; Schmid, C.M.; Ossberger, G.; Stegfellner, L.; Lengauer, G.; Feilmayr, C.; et al. BLASTDAR—A Large Radar Sensor Array System for Blast Furnace Burden Surface Imaging. IEEE Sensors J. 2015, 15, 5893–5909. [Google Scholar] [CrossRef]
  19. Ni, Z.M.; Chen, X.Z. Development of Blast Furnace Charge Level Imaging Technology and Application. Metall. Autom. 2022, 46, 34–45. [Google Scholar]
  20. Shi, Q.D.; Wu, J.X. A Blast Furnace Burden Surface Deep learning Detection System Based on Radar Spectrum Restructured by Entropy Weight. IEEE Sensors J. 2021, 21, 7928–7939. [Google Scholar] [CrossRef]
  21. Wang, H.; Li, W.; Zhang, T.; Li, J.; Chen, X. Learning-Based Key Points Estimation Method for Burden Surface Profile Detection in Blast Furnace. IEEE Sensors J. 2022, 22, 9589–9597. [Google Scholar] [CrossRef]
  22. Chen, Z.; Jiang, Z.H.; Yang, C.H.; Gui, W.H. Detection of Blast Furnace Stockline Based on a Spatial–Temporal Characteristic Cooperative Method. IEEE Trans. Instrum. Meas. 2021, 70, 1–13. [Google Scholar] [CrossRef]
Figure 1. Research framework of the blast furnace burden level prediction and detection.
Figure 1. Research framework of the blast furnace burden level prediction and detection.
Sensors 22 05412 g001
Figure 2. Detection data of radar and mechanical probe.
Figure 2. Detection data of radar and mechanical probe.
Sensors 22 05412 g002
Figure 3. Schematic of the data processing process of the sliding window.
Figure 3. Schematic of the data processing process of the sliding window.
Sensors 22 05412 g003
Figure 4. Space feature regression model curve of radar data.
Figure 4. Space feature regression model curve of radar data.
Sensors 22 05412 g004
Figure 5. Structure of the long-term focus memory network.
Figure 5. Structure of the long-term focus memory network.
Sensors 22 05412 g005
Figure 6. Network structure of ESST-RBFNN.
Figure 6. Network structure of ESST-RBFNN.
Sensors 22 05412 g006
Figure 7. Trend of burden level prediction under normal working conditions.
Figure 7. Trend of burden level prediction under normal working conditions.
Sensors 22 05412 g007
Figure 8. Trend of burden level prediction under abnormal working conditions.
Figure 8. Trend of burden level prediction under abnormal working conditions.
Sensors 22 05412 g008
Figure 9. Comparison chart of different methods for measuring blast furnace burden level under normal working conditions.
Figure 9. Comparison chart of different methods for measuring blast furnace burden level under normal working conditions.
Sensors 22 05412 g009
Figure 10. Comparison chart of different methods for measuring blast furnace burden level under abnormal working conditions.
Figure 10. Comparison chart of different methods for measuring blast furnace burden level under abnormal working conditions.
Sensors 22 05412 g010
Figure 11. Compassion diagram of absolute error.
Figure 11. Compassion diagram of absolute error.
Sensors 22 05412 g011
Figure 12. Compassion diagram of relative error.
Figure 12. Compassion diagram of relative error.
Sensors 22 05412 g012
Figure 13. The 45° line diagram of the measurement effect of the three methods.
Figure 13. The 45° line diagram of the measurement effect of the three methods.
Sensors 22 05412 g013
Figure 14. Operation effect diagram of the detection method of ESST-RBFNN based on the validation data.
Figure 14. Operation effect diagram of the detection method of ESST-RBFNN based on the validation data.
Sensors 22 05412 g014
Table 1. Statistical prediction error comparison.
Table 1. Statistical prediction error comparison.
MethodStatistical Indices
MRERMSECycle-Fit
Radar Probe7.182%0.17699.53%
LFMN8.574%0.23698.87%
Table 2. Statistical error comparison.
Table 2. Statistical error comparison.
MethodStatistical Indices
MRERMSEError-2%Error-5%
Radar8.573%0.184713.97%32.53%
RBFNN7.291%0.182528.91%66.50%
Proposed2.361%0.048091.17%99.33%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, Y.; Chen, Z.; Gui, W.; Yang, C. Real-Time Detection and Short-Term Prediction of Blast Furnace Burden Level Based on Space-Time Fusion Features. Sensors 2022, 22, 5412. https://doi.org/10.3390/s22145412

AMA Style

Chen Y, Chen Z, Gui W, Yang C. Real-Time Detection and Short-Term Prediction of Blast Furnace Burden Level Based on Space-Time Fusion Features. Sensors. 2022; 22(14):5412. https://doi.org/10.3390/s22145412

Chicago/Turabian Style

Chen, Yanli, Zhipeng Chen, Weihua Gui, and Chunhua Yang. 2022. "Real-Time Detection and Short-Term Prediction of Blast Furnace Burden Level Based on Space-Time Fusion Features" Sensors 22, no. 14: 5412. https://doi.org/10.3390/s22145412

APA Style

Chen, Y., Chen, Z., Gui, W., & Yang, C. (2022). Real-Time Detection and Short-Term Prediction of Blast Furnace Burden Level Based on Space-Time Fusion Features. Sensors, 22(14), 5412. https://doi.org/10.3390/s22145412

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop