Method and device for monitoring blast furnace state based on multi-modes fusion

11898212 ยท 2024-02-13

Assignee

Inventors

Cpc classification

International classification

Abstract

A method and device for blast furnace state monitoring based on multi-modes fusion. In the process of blast furnace state monitoring, the concept of sub-mode is introduced, and the method for pre-acquisition and multi-modes weighted fusion by sub-mode is proposed to monitor blast furnace state. After preprocessing the main parameter data of the blast furnace, several sub-modes of the blast furnace and their corresponding blast furnace state indication variable data are obtained by means of a mean shift clustering algorithm. Taking the real-time parameter data of the blast furnace as inputs, the Euclidean distance between the real-time parameter data of the blast furnace and each sub-mode is calculated, and the sub-modes are weighted and fused according to the distance, and the fusion result is the real-time blast furnace state, so as to realize the real-time state monitoring of the blast furnace in a smooth operation.

Claims

1. A method for monitoring blast furnace state based on multi-modes fusion, comprising: step (1) pre-calculation of sub-modes of blast furnace, comprising: step (1.1) acquiring blast furnace historical parameter variable data and constructing a data set; step (1.2) filling in missing values in the data set by a moving average method; step (1.3) detecting abnormal values in the data set by a boxplot method and directly eliminating the detected abnormal values; step (1.4) obtaining blast furnace state indication variable data and preprocessing, comprising filling the missing values and eliminating the abnormal values; step (1.5) calculating a grey relation degree value between each blast furnace parameter variable data and blast furnace state indication variable data, and selecting top N blast furnace parameter variables with the highest correlation as characteristic variables; wherein a calculation formula of the grey relation degree is as follows: ( x 0 , x i ) = 1 n .Math. p = 1 n i ( t p ) where i ( t p ) = min + max .Math. "\[LeftBracketingBar]" x 0 ( t p ) - x i ( t p ) .Math. "\[RightBracketingBar]" + max
.sub.min=custom character(custom character|x.sub.0(t.sub.p)x.sub.i(t.sub.p)|)
.sub.max=custom character(custom character|x.sub.0(t.sub.p)x.sub.i(t.sub.p)|) where x.sub.0 is the blast furnace state indication variable data of n samples, x; is each blast furnace parameter variable of the n samples, i=1, 2, . . . , m, m is a total number of blast furnace parameter variables, is a grey relation coefficient between the two variables, t.sub.p is a p-th moment, p=1, 2, . . . , n, n is a total number of samples contained in each parameter variable, is a resolution coefficient, .sub.min is a two level minimum difference, .sub.max is a two level maximum difference, x.sub.0(t.sub.p) is a blast furnace state indication variable value at t.sub.p and x.sub.i(t.sub.p) is a blast furnace parameter value at t.sub.p; step (1.6) selecting characteristic variable data calculated in step (1.5) from the blast furnace historical parameter variable data set to form a historical characteristic variable data set; step (1.7) applying a mean shift clustering algorithm in the historical characteristic variable data set to obtain a plurality of cluster centers; and step (1.8) calculating the Euclidean distance between the samples in the historical characteristic variable data set and each cluster center, selecting a sample point with a smallest Euclidean distance from each cluster center as the sub-modes of the blast furnace, respectively, and constructing a sub-mode characteristic variable data set by all selected sample point data, storing the blast furnace state indication variable data at a corresponding time of each sub-mode, and constructing a sub-mode indication variable data set; step (2) blast furnace sub-mode fusion, comprising: step (2.1) obtaining real-time parameter variable data of the blast furnace, and selecting the data corresponding to the N characteristic variables calculated in step (1.5) as input variables; step (2.2) calculating the Euclidean distance between input variable data and the characteristic variable data of each sub-mode, and constructing an Euclidean distance matrix D=(d.sub.1, d.sub.2, . . . , d.sub.u, . . . , d.sub.q), where d.sub.u is the Euclidean distance between the input variable data and the characteristic variable data of a u-th sub-mode, and u=1, 2, . . . , q, q is a number of types of the sub-modes obtained in step (1.8); step (2.3) calculating a weight on the basis of an exponential function and a sub-mode contribution rate, comprising: step (2.3.1) scaling the Euclidean distance calculated in step (2.2) based on the exponential function and taking a reciprocal, wherein a calculation formula is as follows: d u_t = 1 e r d u where d.sub.u_t is a transformed Euclidean distance between the input variable data and the characteristic variable data of the u-th sub-mode, and r is a scaling coefficient; step (2.3.2) arranging the sub-modes in an order according to corresponding d.sub.u_t from big to small, and selecting first L sub-modes whose total contribution rate is greater than a set threshold; wherein the calculation formula of the contribution rate con.sub.u of each sub-mode is as follows: c o n u = d u_t .Math. u = 1 q d u_t step (2.3.3) normalizing the d.sub.u_t of the selected sub-mode to obtain a corresponding weight w.sub.u, the corresponding weight of the unselected sub-mode being 0, and obtaining a weight matrix W of the sub-modes; step (3) monitoring of blast furnace state, comprising: reading the indication variable data corresponding to each sub-mode from the sub-mode indication variable data set, and performing weighted summation; wherein the calculation formula is as follows: y z = .Math. u = 1 q w u y z u where .sub.z is an estimated value of a z-th blast furnace state indication variable, z=1, 2, . . . , v, v is a number of types of the blast furnace state indication variables, and y.sub.z.sup.u is a z-th blast furnace state indication variable value corresponding to the u-th sub-mode; and obtaining an estimated value of the corresponding real-time state indication variable of the blast furnace after sub-mode fusion, which realizes real-time state monitoring of the blast furnace.

2. The method according to claim 1, wherein the calculation formula of filling missing values in step (1.2) is as follows: x ( t ) = 1 2 k + 1 .Math. j = - k k x ( t + j ) where t is a time at which there is a missing value, and k is a number of times of moving averages, and x(t) is a parameter value at t.

3. The method according to claim 1, wherein, in step (1.3), the formulas for calculating data distribution inner limits L.sub.1, L.sub.2 and data distribution outer limits L.sub.3, L.sub.4 by the boxplot method are as follows:
L.sub.1=Q.sub.3+1.5(Q.sub.3Q.sub.1)
L.sub.2=Q.sub.31.5(Q.sub.3Q.sub.1)
L.sub.3=Q.sub.3+3(Q.sub.3Q.sub.1) L.sub.4=Q.sub.33(Q.sub.3Q.sub.1) where Q.sub.1 is a lower quartile of data distribution, and Q.sub.3 is an upper quartile of data distribution.

4. The method according to claim 1, wherein in the step (1.7), said applying a mean shift clustering algorithm in the historical characteristic variable data set to obtain a plurality of cluster centers comprises: step (1.7.1) randomly selecting a sample from the historical characteristic variable data set of the blast furnace as a center point c to construct a sliding window; step (1.7.2) calculating an average distance from other sample points to the center point in the sliding window, and taking the average distance as an offset M, wherein the calculation formula is as follows: M = 1 w .Math. y ( t p ) S h [ y ( t p ) - c ] where w is a number of sample points in the sliding window, y is the historical characteristic variable data of the blast furnace, and S.sub.h is a set of all sample points in the sliding window; step (1.7.3) updating a position of the center point according to the offset M; step (1.7.4) repeating steps (1.7.1) to (1.7.3) until the position of the center point does not change; step (1.7.5) when the distance between the center point of a current cluster and a center point of another cluster is less than a set threshold, merging the two clusters into one cluster, and when the distance between the center point of a current cluster and a center point of another cluster is greater or equal to than a set threshold, increasing the number of center points by 1; and step (1.7.6) repeating the steps (1.7.1) to (1.7.5) until all the sample points are traversed, and obtaining a cluster center set of the historical characteristic variable data set of the blast furnace.

5. A device for monitoring blast furnace state based on multi-modes fusion, comprising a memory and one or more processors, wherein executable codes are stored in the memory, and wherein when the processor executes the executable codes, the device is configured to implement the blast furnace state monitoring method based on multi-modes fusion according to claim 1.

6. A non-transitory computer readable storage medium on which a program is stored, wherein the program, when executed by a processor, implements the blast furnace state monitoring method based on multi-modes fusion according to claim 1.

Description

BRIEF DESCRIPTION OF DRAWINGS

(1) FIG. 1 is a flowchart of a method for monitoring blast furnace state based on multi-modes fusion in the embodiment of the present disclosure.

(2) FIG. 2 is an explanatory diagram of the mean shift clustering method in the embodiment of the present disclosure.

(3) FIG. 3 is an estimation result diagram of silicon content in molten iron when the method in the embodiment of the present disclosure is applied to the data of a blast furnace.

(4) FIG. 4 is a structural diagram of a device for monitoring blast furnace state based on multi-modes fusion in the embodiment of the present disclosure.

DESCRIPTION OF EMBODIMENTS

(5) The present disclosure will be further explained with reference to the drawings and specific examples.

(6) FIG. 1 shows the overall flow of the blast furnace state monitoring method based on multi-modes fusion in the embodiment of the present disclosure. There are many variables indicating the state of a blast furnace, including silicon content in molten iron, sulfur content in molten iron, gas utilization rate, etc. This application takes silicon content in molten iron as an example for further explanation. Preprocessing the data of the blast furnace historical parameter variable data set includes filling missing values in the data set and eliminating abnormal values; taking the silicon content in molten iron, an blast furnace state indication variable, as the monitoring target, the characteristic variables of various parameters obtained are selected to form a historical characteristic variable data set; the mean shift clustering algorithm is applied in the historical characteristic variable data set to obtain the sub-modes of the stable state of the blast furnace; according to the Euclidean distance between the real-time parameter variable data and each sub-mode, sub-mode weighted fusion is carried out to obtain the real-time blast furnace state.

(7) The embodiment of the present disclosure relates to a blast furnace state monitoring method based on multi-modes fusion, and each step is specifically described as follows:

(8) (1) The concept of sub-mode is introduced in the process of blast furnace state monitoring and the sub-modes of the blast furnace are obtained; the specific steps include:

(9) (1.1) Obtaining historical state parameter variable data and historical operation parameter variable data of a blast furnace; here, a total of 720 groups of blast furnace parameter variable data in a certain year and a month in a 2650 m.sup.3 blast furnace database in China are selected here. Specifically:

(10) The historical parameter variable data of the blast furnace include top pressure (kPa), total pressure difference (kPa), permeability index (m gas property (min.Math.kPa)), top temperature southeast ( C.), top temperature southwest ( C.), top temperature northeast ( C.), top temperature northwest ( C.), actual wind speed (m/s), bosh gas index (m bosh (min.Math.m.sup.2)), drag coefficient, cold wind flow (10,000 m wind flow), oxygen enrichment flow (m, oxygen enrichment), oxygen enrichment rate (%), blast momentum (kg.Math.m/s), cold wind pressure (Pa), blast humidity (g/m.sup.3), pulverized coal injection (ton), theoretical burning temperature ( C.).

(11) (1.2) The missing values of the blast furnace historical parameter variable data are filled. During the operation of the industrial blast furnace, due to the failure or maintenance of sensors, the data in the database is often missing, and the real-time information of the blast furnace is incomplete, which leads to the decrease of the accuracy of its state monitoring. Therefore, it is necessary to fill in the missing data. The moving average method is adopted to fill the missing values in the data set, and the missing value x at time t is calculated by the following formula:

(12) x ( t ) = 1 2 k + 1 .Math. j = - k k x ( t + j )

(13) where k is the number of moments of moving averages, which is set to 3 here.

(14) (1.3) The abnormal values of blast furnace historical parameter variable data are eliminated. In the process of industrial blast furnace operation, due to the violent fluctuation in the blast furnace or the blast furnace failure, the observed values of some variables will obviously deviate from the normal range of the variable data in the database, so the data that cannot represent the state of the blast furnace in a smooth operation will be eliminated directly. A boxplot method is used to detect abnormal values in blast furnace data. The boxplot method uses five statistical values of data to describe data distribution, including maximum, minimum, median, upper quartile Q.sub.3 and lower quartile Q.sub.1. According to these statistical values, two inner limits L.sub.1, L.sub.2 and two outer limits L.sub.3, L.sub.4 of data distribution are calculated. Their calculation formulas are as follows:
L.sub.1=Q.sub.3+1.5(Q.sub.3Q.sub.1)
L.sub.2=Q.sub.31.5(Q.sub.3Q.sub.1)
L.sub.3=Q.sub.3+3(Q.sub.3Q.sub.1)
L.sub.4=Q.sub.33(Q.sub.3Q.sub.1)

(15) When a data is between the inner limit and the outer limit, it is called a mild abnormal value; when a data is outside the outer limit, it is called an extreme abnormal value. Here, both kinds of abnormal values are eliminated.

(16) (1.4) Taking the silicon content (%) in molten iron, an indication variable of direct state of the blast furnace, as the monitoring target of blast furnace state; reading the corresponding data and preprocessing, including filling the missing values in the data set and eliminating the abnormal values; the specific operation method is the same as that described in steps (1.2) and (1.3).

(17) (1.5) The correlation degree between the blast furnace parameter variables read in step (1.1) and the silicon content in molten iron are calculated by using the grey relation analysis method. The state indication variable of the blast furnace, the silicon content in molten iron, is a reference variable, and each parameter variable of the blast furnace is a comparison variable. Then, the grey relation coefficient between the reference variable x.sub.0 containing n samples and the comparison variable x; containing n samples is calculated as follows:

(18) 0 i ( t p ) = min + max | X 0 ( t p ) - x i ( t p ) | + max

(19) where i=1, 2, . . . , m, m is the total number of blast furnace parameter variables; p=1, 2, . . . , n, n is the total number of samples contained in each parameter; is the resolution coefficient, taking 0.5 here; .sub.min is a two level minimum difference, .sub.max is a two level maximum difference. Their calculation formulas are as follows:
.sub.min=min.sub.i(min.sub.p|x.sub.0(t.sub.p)x.sub.i(t.sub.p)|)
.sub.max=max.sub.i(max.sub.p|x.sub.0(t.sub.p)x.sub.i(t)|)

(20) The average value of the grey relation coefficient of two variables is the grey relation degree between the two, and the calculation formula is as follows:

(21) ( x 0 , x i ) = 1 n .Math. p = 1 n i ( t p )

(22) m kinds of blast furnace parameter variables are rearranged according to their corresponding grey relation degree values from big to small, and the first N kinds of blast furnace parameter variables with the highest relation degree with the silicon content in molten iron, which is an indication variable of blast furnace state, are selected as the final characteristic variables, where Nis set to 7. Some results calculated in this example are shown in the following table:

(23) TABLE-US-00001 Parameters of blast Grey relation furnace degree Top temperature 0.7569 northeast Top temperature 0.7539 northwest Top temperature 0.7486 southeast Drag coefficient 0.7428 Top temperature 0.7416 southwest Gas permeability 0.7260 Oxygen enrichment rate 0.7236 Cold wind pressure 0.7209 Oxygen enrichment 0.7190 flow rate Pulverized coal 0.7186 injectoin . . . . . .

(24) Therefore, seven characteristic variables are selected: top temperature northeast, top temperature northwest, top temperature southeast, drag coefficient, top temperature southwest, gas permeability and oxygen enrichment rate.

(25) (1.6) In the historical characteristic variable data set, a number of blast furnace sub-modes under the condition of smooth operation are obtained by applying the mean shift clustering algorithm. The specific steps include:

(26) FIG. 2 is an explanatory diagram of the mean shift clustering method in the embodiment of the present disclosure.

(27) (1.6.1) Randomly selecting a sample from the historical characteristic variable data set of the blast furnace as a center point c to construct a sliding window with a radius R, wherein the radius R is determined according to the distribution of samples.

(28) (1.6.2) Calculating an average distance from other sample points to the center point in the sliding window, and taking the average distance as an offset M, wherein the calculation formula is as follows:

(29) M = 1 w .Math. y ( t p ) S h [ y ( t p ) - c ]

(30) where w is a number of sample points in the sliding window, y is the historical characteristic variable data of the blast furnace, and S.sub.h is a set of all sample points in the sliding window.

(31) (1.6.3) Updating a position of the center point according to the offset M; wherein the calculation formula is as follows:
c.sup.j+1=M.sup.l+c.sup.l
where l is the number of iterations;

(32) (1.6.4) Repeating steps (1.6.1) to (1.6.3) until the position of the center point does not change.

(33) (1.6.5) If the distance between the center point of a current cluster and a center point of another cluster is less than a set threshold, merging the two clusters into one cluster, otherwise, increasing the number of center points by 1.

(34) (1.6.6) Repeating steps (1.6.1) to (1.6.5) until all the sample points are traversed, and obtaining a cluster center set of the historical characteristic variable data set of the blast furnace.

(35) (1.7) In this embodiment, a total of 24 sub-modalities are obtained, the Euclidean distance between the samples in the historical characteristic variable data set and each cluster center is calculated, and the sample points with the smallest Euclidean distance from each cluster center are selected as the sub-modes of the blast furnace, and all the selected sample point data are used to construct a sub-mode characteristic variable data set; at the same time, the data of blast furnace state indication variables at the corresponding time of each sub-mode are stored, and a sub-mode indication variable data set is constructed.

(36) (2) In the process of blast furnace state monitoring, a sub-mode weighted fusion strategy is introduced and the weights of sub-modes are calculated. The specific steps include:

(37) (2.1) Reading the real-time parameter variable data of the blast furnace, and selecting the seven characteristic variables selected in step (1.5), including: top temperature northeast, top temperature northwest, top temperature southeast, drag coefficient, top temperature southwest, gas permeability and oxygen enrichment rate, and forming the input variable data y.sub.c of the model.

(38) (2.2) Calculating the Euclidean distance between the input variable data and the characteristic variable data of each sub-mode, and constructing an Euclidean distance matrix D=(d.sub.1, d.sub.2, . . . , d.sub.u, . . . , d.sub.q), where d.sub.u is the Euclidean distance between the input variable data and the characteristic variable data of the u-th sub-mode, u=1, 2, . . . , q, q is the number of sub-modes obtained in step (1.7), and in this embodiment, q is 24.

(39) (2.3) Calculating the weight of each sub-mode based on the Euclidean distance between the input variable data and each sub-mode; the specific steps are as follows:

(40) (2.3.1) Scaling the Euclidean distance calculated in step (2.2) based on the exponential function and taking a reciprocal, wherein a calculation formula is as follows:

(41) d u - t = 1 e r d u where d.sub.u_t is a transformed Euclidean distance between the input variable data and the characteristic variable data of the u-th sub-mode, and r is a scaling coefficient;

(42) (2.3.2) Arranging the sub-modes in an order according to corresponding d.sub.u_t from big to small, and selecting first L sub-modes whose total contribution rate is greater than a set threshold; wherein the calculation formula of the contribution rate con.sub.u of each sub-mode is as follows:

(43) c o n u = d u_t .Math. u = 1 q d u_t

(44) (2.3.3) Normalizing the d.sub.u_t of the selected sub-mode to obtain a corresponding weight w.sub.u, and the corresponding weight of the unselected sub-mode being 0, thus obtaining a weight matrix W of the sub-modes;

(45) w u = d u_t .Math. u U L d u_t

(46) where U.sub.L is a set of selected L sub-modes.

(47) Part of the results of sub-mode weights calculated in this embodiment is shown in the following table:

(48) TABLE-US-00002 Moment Moment Moment Moment Moment Moment Moment Moment Moment Moment 1 2 3 4 5 6 7 8 9 10 . . . Sub-mode 0 0 0 0 0 0 0 0 0 0 . . . 1 Sub-mode 0 0 0 0 0 0 0 0 0 0 . . . 2 Sub-mode 0 0 0.407 0.128 0 0 0 0 0 0 . . . 3 Sub-mode 0 0 0 0 0 0 0 0 0 0 . . . 4 Sub-mode 0 0.624 0.539 0.872 1 1 0.118 0.129 0.782 0 . . . 5 Sub-mode 0 0 0 0 0 0 0 0 0 0 . . . 6 Sub-mode 0 0 0 0 0 0 0 0 0 0 . . . 7 Sub-mode 0 0 0 0 0 0 0 0 0 0 . . . 8 Sub-mode 0.063 0 0 0 0 0 0.334 0.383 0.219 0.356 . . . 9 Sub-mode 0 0 0 0 0 0 0 0 0 0 . . . 10 Sub-mode 0.071 0 0 0 0 0 0.183 0.169 0 0.201 . . . 11 Sub-mode 0 0 0 0 0 0 0 0 0 0 . . . 12 Sub-mode 0.868 0.271 0 0 0 0 0 0 0 0 . . . 13 Sub-mode 0 0 0 0 0 0 0 0 0 0 . . . 14 Sub-mode 0 0 0 0 0 0 0 0 0 0 . . . 15 Sub-mode 0 0 0 0 0 0 0 0 0 0 . . . 16 Sub-mode 0 0 0 0 0 0 0.365 0.318 0 0.443 . . . 17 Sub-mode 0 0 0 0 0 0 0 0 0 0 . . . 18 Sub-mode 0 0 0 0 0 0 0 0 0 0 . . . 19 Sub-mode 0 0 0 0 0 0 0 0 0 0 . . . 20 Sub-mode 0 0 0.054 0 0 0 0 0 0 0 . . . 21 Sub-mode 0 0 0 0 0 0 0 0 0 0 . . . 22 Sub-mode 0 0 0 0 0 0 0 0 0 0 . . . 23 Sub-mode 0 0.106 0 0 0 0 0 0 0 0 . . . 24

(49) (3) The state monitoring of the blast furnace is realized based on the sub-mode weight calculation results, and the specific steps include:

(50) After the real-time parameter variable data of the blast furnace is input, the weight matrix is calculated according to step (2), and the silicon content data in molten iron corresponding to each sub-mode is read from the sub-mode silicon content data set, and the weighted sum is carried out. The calculation formula is as follows:

(51) y Si ^ = .Math. u = 1 q w u y S i u where .sub.Si is the estimated value of silicon content in molten iron, and y.sub.Si.sup.u is the value of silicon content in molten iron corresponding to the u-th sub-mode. Therefore, the estimated value of silicon content in blast furnace molten iron after sub-mode fusion is obtained.

(52) In this embodiment, 720 sets of blast furnace parameter variable data from the blast furnace database in another month of the same year are selected to test the model.

(53) Two statistical indicators, Mean Absolute Percentage Error (MAPE) and Root Mean Square Error (RMSE), are used to evaluate the estimation effect of the model. Their calculation formulas are as follows:

(54) RMSE = 1 n .Math. p = 1 n [ s ( t p ) - s ( t p ) _ ] 2

(55) MAPE = 1 n .Math. p = 1 n .Math. "\[LeftBracketingBar]" s ( t p ) _ - s ( t p ) s ( t p ) .Math. "\[RightBracketingBar]" 100 %

(56) where s(t.sub.p) and s(t.sub.p) are the real and estimated values of the data, respectively.

(57) As shown in FIG. 3, because there are too many data samples, only the test results of the first 200 samples are selected for display. The dotted line shows the estimated value of silicon content in molten iron based on multi-modes fusion, and the solid line shows the actual value of silicon content in molten iron. As shown in the figure, the change trends of the solid line and the dotted line are roughly the same, and the fitting effect is ideal, which proves the effectiveness and applicability of the method provided by the present disclosure. In addition, RMSE and MAPE of the estimated silicon content based on multi-modes fusion are 0.0803 and 14.99% respectively, which more intuitively proves the effectiveness and applicability of the proposed method.

(58) Corresponding to the embodiment of the blast furnace state monitoring method based on multi-modes fusion, the present disclosure also provides an embodiment of the device for monitoring blast furnace state based on multi-modes fusion.

(59) Referring to FIG. 4, a device for monitoring blast furnace state based on multi-modes fusion provided by an embodiment of the present disclosure includes a memory and one or more processors. The executable code is stored in the memory, and when the processors execute the code, the above blast furnace condition monitoring method based on multi-modes fusion is implemented.

(60) The embodiment of the device for monitoring blast furnace state based on multi-modes fusion can be applied to any equipment with data processing capability, and any equipment with data processing capability can be equipment or devices such as computers. The embodiment of the device can be realized by software, or by hardware or a combination of hardware and software. Take the software implementation as an example, as a logical device, it is formed by reading the corresponding computer program instructions in the nonvolatile memory into the memory through the processor of any device with data processing capability. From the aspect of hardware, as shown in FIG. 4, it is a hardware structure diagram of any equipment with data processing capability in the device monitoring blast furnace state based on multi-modes fusion according to the present disclosure. Except the processor, memory, network interface and non-volatile memory shown in FIG. 4, any equipment with data processing capability in the embodiment of the present disclosure can usually include other hardware according to the actual functions of the equipment with data processing capability, which will not be repeated here.

(61) The realization process of the functions and functions of each unit in the above-mentioned device is detailed in the realization process of the corresponding steps in the above-mentioned method, which will not be repeated here.

(62) An embodiment of the present disclosure also provides a computer readable storage medium on which a program is stored. When the program is executed by a processor, it implements the blast furnace state monitoring based on multi-modes fusion in the above embodiment.

(63) It should be noted that the above embodiments are only used to illustrate, rather than to limit the technical solution of the present disclosure. Although the present disclosure has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that the technical solution of the present disclosure can be modified or equivalently replaced without departing from the spirit and scope of the technical solution of the present disclosure, which should be covered by the claims of the present disclosure. The above specific embodiments have explained the technical solution and beneficial effects of the present disclosure in detail. It should be understood that the above embodiments are only the most preferred embodiments of the present disclosure, and are not intended to limit the present disclosure. Any modification, supplement, equivalent substitution and the like made within the scope of the principles of the present disclosure shall be included in the scope of protection of the present disclosure.