JP2002194405A - Method and device for monitoring operation of blast furnace, and computer-readable recording medium - Google Patents

Abstract

PROBLEM TO BE SOLVED: To correctly monitor the operational condition of a blast furnace, and to predict any operational abnormality such as channeling. SOLUTION: The distribution condition of the measurement data collected by a plurality of sensors installed on a blast furnace facility is arranged on the two-dimensional plane or three-dimensional space reflecting the installation positions of each sensor on the blast furnace facility, the pattern formed by the measurement data using isopleths or the like, and the pattern or the characteristics information of the pattern is operated by the image processing. The operational condition of the blast furnace is monitored by comparing the obtained pattern or characteristics information of the pattern with the preset pattern or characteristics information of the pattern, the pattern or characteristics information of the pattern is updated corresponding to the temporal transition of the measurement data, and any operational abnormality of the blast furnace such as channeling is predicted by comparing the temporal transition with the transition condition of the preset pattern or characteristics information of the pattern.

Description

DETAILED DESCRIPTION OF THE INVENTION

[0001]

BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a method for monitoring the operating state of a blast furnace during operation and predicting an abnormal operation such as a blow-through, an apparatus used for implementing the method, and a computer-readable recording medium.

[0002]

2. Description of the Related Art Conventionally, a method for monitoring and estimating abnormal operation of a blast furnace has been disclosed in JP-A-5-156328.
And Japanese Patent Application Laid-Open No. 11-140520. Each of these monitoring and prediction methods collects measurement data from each sensor without reflecting the installation position information of each sensor on the blast furnace equipment, and uses a reference value or a simple physical model set in advance. The operation state is monitored and the operation abnormality is predicted by comparing with the threshold value.

[0003]

However, the blast furnace process to which the present invention is directed is to be treated as a distributed constant process having dynamic characteristics. Therefore, measurement data of a plurality of various sensors installed with distribution on the blast furnace equipment are collected independently of each other,
It should not be evaluated, but should be collected and evaluated in relation to the installation location on the blast furnace facility where each sensor is installed.

[0004] In the conventional system, the installation positions of such sensors are not collected and evaluated in association with the measurement data. As a result, there is a problem that the accuracy of monitoring and prediction of the operation state of the blast furnace is low. Was.

The present invention has been made in view of the above circumstances, and has as its object to enable accurate monitoring of the operation state of a blast furnace and prediction of an operation abnormality.

[0006]

The operation monitoring method for blast furnace operation according to the present invention uses a two-dimensional plane or a two-dimensional plane reflecting the installation position of each sensor by measuring the measurement target amount from a plurality of sensors installed in the blast furnace. A method of monitoring the operating state of a blast furnace by arranging in a three-dimensional space, expressing the distribution state of each measurement data as a figure or characteristic information of the figure formed by the three-dimensional space, and evaluating these, The method is characterized in that an arbitrary contour having the same measurement target quantity is calculated on a plane or a three-dimensional space, and a figure formed by the contours or characteristic information of the figure is evaluated.

[0007] Another operation monitoring method in blast furnace operation according to the present invention is a method of monitoring measurement data of an object to be measured from a plurality of sensors installed in a blast furnace in a two-dimensional plane or a three-dimensional space reflecting the installation position of each sensor. A method of monitoring the operating state of a blast furnace by arranging and expressing the distribution state of each measurement data as a graphic or characteristic information of the graphic formed by the two or three-dimensional planes. Calculate any isolines with the same measurement target quantity in space, and the number, position, area, center of gravity, aspect ratio of the figures formed by the isolines, the maximum or minimum value in the figures, The feature is that at least one feature information of the average value and the variance is calculated and evaluated by image processing.

Another feature of the operation monitoring method in the blast furnace operation according to the present invention is that a figure formed by isolines or feature information of the figure is updated in accordance with the time transition of each measurement data. The point is to compare these temporal transitions with transition conditions set in advance.

Another feature of the operation monitoring method in blast furnace operation according to the present invention is that the number, position and area of figures formed by isolines corresponding to the time transition of each measurement data. , The center of gravity, the aspect ratio of the figure, the maximum value or the minimum value in the figure, the average value, the temporal transition of at least one characteristic information of the variance is calculated by image processing, and these temporal transitions are set in advance with transition conditions. The point is to compare.

Another feature of the operation monitoring method in the blast furnace operation according to the present invention is that the procedure for calculating the isolines is performed on data arranged in an uneven positional relationship on a two-dimensional plane. Select a quadrilateral element in which one of the internal angles does not exceed 180 degrees, set the average value of the data of the four vertices at the intersection of the diagonal lines, and search for the isoline using the triangular element having this intersection as the vertex. It is a technique of drawing.

An operation monitoring apparatus for blast furnace operation according to the present invention includes a data collection unit that collects measurement data measured by various sensors installed on a plurality of blast furnace facilities, and a distribution state of the collected measurement data of each sensor. A contour line calculation unit for arranging in a two-dimensional plane or three-dimensional space reflecting the installation position on the blast furnace equipment and calculating an arbitrary contour line having the same measurement target quantity, and a figure formed by the contour lines Or a graphic feature information calculation unit for calculating graphic feature information by image processing;
It is characterized in that it comprises an operation monitoring unit for comparing the graphic obtained by the graphic characteristic information calculation unit or the characteristic information of the graphic with the preset graphic or characteristic information of the graphic and monitoring the operation.

Another feature of the operation monitoring apparatus in the blast furnace operation of the present invention is that the operation of the contour line calculating section and the graphic characteristic information calculating section is repeated in accordance with the time transition of each measurement data, An operation for predicting an operation abnormality such as a blow-by by comparing a graphic characteristic information transition calculating unit for calculating these temporal transitions and a transition condition of a graphic or graphic characteristic information in which the graphic characteristic information transition information is set in advance. And a prediction unit.

Further, another feature of the operation monitoring apparatus in the blast furnace operation of the present invention is that the contour line calculating section comprises:
For data arranged in a non-uniform positional relationship on a two-dimensional plane, a quadrangle element whose inner angle does not exceed 180 degrees is selected, and the intersection of diagonal lines is set to the average value of the data of four vertices. This is a point having a function of searching for and drawing isolines using a triangular element having the intersection as a vertex.

Another feature of the operation monitoring apparatus for blast furnace operation according to the present invention is that the apparatus has an output unit for visualizing the transition of the figure and the characteristic information of the figure.

The computer-readable recording medium of the present invention is characterized in that a program for causing a computer to execute the processing procedure of the operation monitoring method in the blast furnace operation is recorded.

Further, another computer-readable recording medium of the present invention is characterized in that a program for causing a computer to function as each unit of the operation monitoring device in the blast furnace operation is recorded.

[0017]

BRIEF DESCRIPTION OF THE DRAWINGS FIG. 1 is a block diagram showing a blast furnace operation monitoring method and apparatus according to an embodiment of the present invention;

FIG. 1 is a block diagram showing a configuration of an operation monitoring device in blast furnace operation of the present embodiment. In FIG. 1, a plurality of various sensors for measuring a stave temperature, a shaft pressure, and the like are provided on a blast furnace facility 1. FIG.
In the example, stave temperature and shaft pressure are taken as examples to show the positions of multiple sensors installed on the outer surface of the blast furnace equipment, but the same applies when multiple sensors are installed inside the blast furnace equipment. The arrangement of the various sensors on the blast furnace equipment may be at irregular intervals.

A plurality of various sensors 2 on the blast furnace equipment 1 measure physical quantities such as temperature, pressure, flow rate, particle size, density, and composition. In the following, the sensor for measuring the temperature is shown in FIG.
An example in which a plurality of blast furnaces are arranged on the outer surface of the blast furnace as shown in FIG.

The installation position information of each of the plurality of temperature sensors arranged on the outer surface of the blast furnace is represented by three-dimensional spatial coordinates (x (i), y
(i), z (i)), where i = 1, 2, 3,..., N (N: the number of temperature sensors).

In the data collection device 3, measurement data output from a plurality of temperature sensors arranged on the blast furnace equipment is sampled and collected at a preset sampling period Δt. The sampling period Δt can be arbitrarily set at a time interval of several ms or more corresponding to the processing capacity of the data collection device 3 and the processing capability of the data processing device 4 and the time intervals required for operation monitoring and operation prediction. The temperature data collected by the data collection device 3 is sent to the data processing device 4 in real time.

In the contour line calculating section 5, the temperature data input from the data collection device 3 is arranged on a two-dimensional plane or a three-dimensional space reflecting information on the positions of the sensors on the blast furnace equipment. An equivalent isovalue is calculated, and a figure formed by the isovalue is calculated.

Hereinafter, an example of a method of calculating the contour line in the contour line calculating section 5 will be described. FIG. 2 shows that the contour line calculator 5
The figure shows an example of a figure formed by defining a two-dimensional plane having an r-axis in the furnace circumferential direction of the blast furnace and an h-axis in the furnace height direction, and forming contour lines on the two-dimensional plane. In FIG. 2, ●
The mark indicates that a plurality of temperature sensor installation positions arranged on the outer shape of the blast furnace are arranged by transforming their three-dimensional spatial coordinates (x (i), y (i), z (i)). .

In FIG. 2, the coordinate transformations are tuyere diameter, furnace belly diameter,
The measurement was performed using a geometric relationship such as calculating a projection onto the two-dimensional plane from a furnace bottom diameter, a shaft angle, a Bosch angle (morning glory angle), and the like. The method in the present invention defines 2
The dimensional plane need not be limited to a square plane as shown in FIG. 2, and a partially fan-shaped two-dimensional plane may be defined according to the shaft angle and the Bosch angle (morning glory angle).

FIG. 2 defines and uses a two-dimensional plane having an r-axis in the furnace circumferential direction and an h-axis in the furnace height direction of the blast furnace for the purpose of explanation. The same applies to the case of arranging in a three-dimensional space according to spatial coordinates.

By disposing the corresponding measurement data at the point indicated by a black mark indicating the temperature sensor installation position on the two-dimensional plane in FIG. 2, the distribution state of the temperature data at a certain time t can be expressed. At this time, the intervals between the marks may be spatially unequal, and need not be spatially equal, by an isoline search method described later.

Based on the temperature data arranged at the position of the mark, the temperature data in the mutual space of the mark is spatially interpolated to search for an isoline. Here, the isoline is obtained by connecting points indicating the same value from the spatially distributed temperature data by a line.

In order to search for isolines with respect to temperature data distributed at spatially unequal positions, a method using a triangular element composed of temperature sensor installation points is reliable.
There are enormous degrees of freedom in combining triangle elements in space. Further, when the number of measurement points is small in the space, there is a problem that the shape of the obtained isoline differs depending on the selection of the triangular element.

Therefore, as a method of reducing the degree of freedom of element selection and facilitating the selection, and reducing the error of the isoline shape due to the element selection, a method using a triangular element using a quadrilateral element as an apex is used. The "contour search method" is exemplified.

This method will be described with reference to FIG. 2 of FIG.
Temperature sensor installation position on the dimensional plane ● For all the points indicated by the mark, each point is associated in advance so that one of the inner angles is constituted by a square element not exceeding 180 degrees. According to the element selection condition for the square element, the degree of freedom of element selection is reduced, and element selection can be facilitated. In the case of a blast furnace facility, the coordinates of each sensor position are known, so that the association may be made once, or the association may be made automatically using an automatic search algorithm as a combination problem.

In FIG. 3, any quadrangular element in which one of the interior angles does not exceed 180 degrees, ie, each vertex P1, P2, P3,
Temperature sensor measurement data at point P4 are T1, T2, T
Here is an example that is 3, T4. The temperature at the intersection of the diagonal lines of the square elements, that is, the temperature at the point Pm indicated by the circle in FIG. 3 is defined as Tm. Tm is an average value calculated from T1, T2, T3, and T4, and is defined as, for example, an arithmetic average. Tm = (T1 + T2 + T3 + T4) / 4 (1)

Next, four triangular elements having a vertex at the intersection Pm on the diagonal are defined inside the quadrilateral element, and the temperature data on the side of each triangular element is obtained by calculating the vertices of the vertices constituting both ends of the side. It shall be obtained by interpolation with temperature data. In interpolation, primary interpolation method, etc.
Any method may be used.

Suppose that the value of the contour line to be searched for is T, and the following equation T1 <T <T4 (2) T1 <T <T2 (3) It is assumed that there is a relationship.

In the example shown in FIG. 3, T
Always exists on the straight line connecting P1 and P4, and always P1
Exists as a temperature data point interpolated on a straight line connecting Pm and Pm or a straight line connecting Pm and P4. Here, if it is assumed that T1 <T <Tm (4), a temperature data point of T exists on a straight line connecting P1 and Pm. The points of these temperature data T are indicated by a mark. Similarly, from the condition of equation (3), T always exists as a temperature data point interpolated on a straight line connecting P1 and P2, which is indicated by a triangle. Temperature obtained from above
By connecting the points of T with straight lines, it is possible to search for a contour line of the temperature T in the square element of interest.

In the above example, if Tm <T <T4... (5) instead of equation (4), and T2 <T <T3. Are as shown by the squares, and the isolines connecting these with straight lines can be shown by broken lines.

By repeating the above processing for all the quadrangular elements in the space, the search and drawing of the isolines in the space are completed. As illustrated in FIG. 2, the temperature data forms a certain figure in a two-dimensional plane by the obtained isolines. In particular, isolines that form closed curves form certain characteristic figures. In FIG. 2, the isolines of a certain temperature T are indicated by solid lines, and the figure enclosed by the closed curve is indicated by hatching. Dashed lines are other temperature isolines.

As described above, for data arranged in a spatially unequal positional relationship, a quadrilateral element whose one of the internal angles does not exceed 180 degrees is selected, and the data of the four vertices is set at the intersection of the diagonal lines. The method of setting the average value and searching for and drawing isolines using a triangular element having this intersection point as a vertex is different from the method of searching for isolines using only triangular elements,
While reducing the degree of freedom of element selection and making selection easier,
Since a triangular element having an average value of each vertex of a quadrilateral element as a vertex is used, this is an effective method that can reduce a search error of an isoline depending on element selection. Since triangular elements are used in the final stage of the search, it goes without saying that there is no problem that the contour to be searched intersects with another contour in the middle or that the contour is interrupted in the middle.

The search method is not limited to a two-dimensional plane, but is a feasible and effective method for a three-dimensional space figure composed of rectangular plane elements.

In the present invention, it is not necessary to limit the method of searching for contour lines, and it is possible to draw contour lines using other methods or triangular elements on a two-dimensional plane or three-dimensional space. I do not care.

The graphic characteristic information calculating section 6 performs image processing on the graphic calculated by the iso-line calculating section 5 and obtains the graphic and characteristic information of the graphic, that is, the number, position, area, center of gravity, aspect ratio of the graphic. , The maximum value or the minimum value, the average value, and the variance in the figure.

The operation monitoring unit 7 can monitor the operation of the blast furnace by comparing the figure calculated by the figure characteristic information calculation unit 6 and the characteristic information of the figure with the predetermined figure and the characteristic information of the figure. is there. 4 to 7,
The operation monitoring method will be described using temperature data as an example.

FIGS. 4 to 7 show examples of graphic feature information calculation performed by the graphic feature information calculation unit 6 by performing image processing, which are arranged in accordance with the transition of time.
FIGS. 4 to 7 show gray levels such that a graphic area in an isoline at a higher temperature is whiter and a graphic area in an isoline at a lower temperature is blacker. Further, in FIGS. 4 to 7, a figure formed by isolines at a certain high temperature is calculated by image processing, the figure is surrounded by lines, and labeling is performed.

FIG. 4 shows a case where the operation of the blast furnace is stable, and a figure formed by isolines at a certain high temperature caused by the root of the cohesive zone widely exists over the entire furnace periphery at the lower part of the blast furnace. It is in the state that it is.

FIG. 5 shows a state in which a certain time has elapsed from the state of FIG. 4, and a figure formed by isothermal lines at a certain high temperature caused by the root of the cohesive zone due to operational disturbance is shown in the furnace circumferential direction. At a certain point, the state is expanding in the furnace height direction.

FIG. 6 shows the state of the figure formed by the contour lines at a certain high temperature when the time further elapses from the state of FIG. This shows a state in which the position of the blast furnace has moved upward and has reached the center position in the height direction of the blast furnace equipment.

FIG. 7 shows a state in which more time has passed since the state of FIG. 6, and a figure formed by isolines at a certain high temperature shows that most of the area has passed to the upper part of the blast furnace equipment.
FIG. 4 shows a state in which the area of the remaining figure exists at a substantially central position in the furnace height direction and at a position about three-quarters in the furnace circumferential direction, and the operation has not returned to the stable operation state as shown in FIG. This indicates an abnormal state.

FIGS. 6 and 7 show a situation in which a figure formed by isolines due to a certain high temperature is flowing to the upper part of the blast furnace equipment.
That is, it indicates an operation abnormality, and in this case, a so-called blow-by phenomenon.

At this time, the operation monitoring unit 7 can monitor the operation by comparing the figure calculated by the figure characteristic information calculation unit 6 and the characteristic information of the figure with the preset figure and the characteristic information of the figure. It is.

Further, the graphic feature information transition calculating section 8 calculates the temporal transition of the graphic calculated by the graphic characteristic information calculating section 6 and the characteristic information of the graphic. The operation prediction unit 9 compares the time transition of the graphic and the characteristic information of the graphic calculated by the graphic characteristic information transition calculation unit 8 with the predetermined transition condition of the graphic and the characteristic information of the graphic to predict the operation state. It is possible.

Referring to FIG. 8, an operation prediction method using temperature data as an example will be described. FIG. 8 shows the position of the center of gravity of the figure formed by isolines at a certain high temperature illustrated in FIGS.
The process of transition from FIG. 4 to FIG. 7 is illustrated by taking the position of the center of gravity on the vertical axis and time on the horizontal axis.

In the operation prediction section 9, an upper limit management value of the position of the center of gravity is set in advance in order to predict the operation state. In the stable operation state as shown in FIG. 4, the center of gravity G (t) of the target graphic is smaller than the set upper limit management value Gu. From the position of the center of gravity G (t) and its time rate of change dG (t), the operating state after a certain time Δt, that is, the position of the center of gravity G (t + Δt) after Δt is given by G (t + Δt) = G (t ) + DG (t) · Δt (7), and if G (t + Δt) <Gu (8), it is possible to predict that the stable operation state will continue after Δt. It is.

On the other hand, in an operation fluctuation state as shown in FIG. 6, the center of gravity G (t) of the target graphic is set to the set upper limit management value.
Although smaller than Gu, the center of gravity position G (t) and its time rate of change dG
From (t), when the operation state after a certain time Δt, that is, the gravity center position G (t + Δt) after Δt is predicted in the same manner as in the equation (7), G (t + Δt)> Gu (9) ), It is possible to predict that an operation abnormality such as blow-by will occur in the operation state after Δt.

In FIG. 8, as the characteristic information of the figure formed by the isolines at a certain high temperature, the value of the center of gravity G (t) of the figure obtained by image processing and its time rate of change dG (t) are shown as an example. It has been shown that the method according to the present invention makes it possible to predict an operation abnormality.However, in addition to the position of the center of gravity, the above-described graphic characteristic information obtained by image processing, a method for evaluating the time change rate, It is also effective to use a method of evaluating by combining some pieces of characteristic information, a method of using not only an upper limit management value but also a lower limit management value, or a method of evaluating by combining an upper limit management value and a lower limit management value.

Further, in the recording unit 10, the calculation result in the graphic feature information transition calculating unit 8 is recorded as a file in a text format or the like, and is made into a database. In recording the transition of the graphic and the graphic feature information, the calculation result can be recorded as a moving image file in an AVI format or the like.
At this time, in implementing the method of monitoring the operation of the blast furnace according to the present invention, it is possible to efficiently record and create a database by removing redundant moving image information using various data compression methods as necessary. is there. In the method of the present invention, it is not necessary to limit the data compression method.

It is also possible to input the information recorded by the recording unit 10 into a file and evaluate the operating state of the blast furnace off-line.

The output unit 11 outputs the transition of the graphic and the characteristic information of the graphic, the operation monitoring result, and the operation prediction result on a screen by a monitor or the like.

In the present embodiment, the method of the present invention has been described using stave temperature data as an example. However, the method of the present invention does not need to be limited to stave temperature data. It goes without saying that the data is also effective.

The data processing apparatus 4 according to the embodiment described above
Is a computer CPU or MPU, RAM, RO
M and the like, and can be realized by operating a program recorded in a RAM or a ROM. Therefore, the present invention can be realized by recording a program that causes a computer to perform the above function on a recording medium and causing the computer to read the program. As a recording medium, a CD-ROM, a DVD, a floppy disk, a hard disk, a magnetic tape, a magneto-optical tape, a nonvolatile memory card, or the like can be used.

The functions of the above-described embodiments can be realized not only when the computer executes the supplied program, but also when the program code is executed by the operating system (OS) or other application software running on the computer. Needless to say, the program code is also included in the embodiment of the present invention when the functions of the above-described embodiment are realized in cooperation with the above.

[0060]

The method of the present invention described in detail above evaluates measurement data of a plurality of various sensors arranged in a blast furnace facility with a spatial distribution by reflecting the installation position information of each sensor. And evaluates the spatial distribution information of the measurement data as graphic information on a two-dimensional plane or a three-dimensional space reflecting the installation position information of each sensor. By calculating and evaluating by image processing, it becomes possible to accurately monitor the operation state of the blast furnace and predict the operation abnormality.

[Brief description of the drawings]

FIG. 1 is a block diagram illustrating a configuration of an operation monitoring device according to an embodiment.

FIG. 2 is a diagram showing a figure formed by searching for isolines from measurement data of a plurality of sensors distributed at spatially unequal positions and forming certain isolines.

FIG. 3 is a diagram illustrating an example of a method in which a contour line calculation unit searches for a contour line from measurement data of a plurality of sensors distributed at spatially unequal positions. FIG. 11 is a diagram for explaining a “contour search method using elements”.

FIG. 4 is a diagram illustrating an example of calculation of graphic feature information in a process from a stable operation to an abnormal operation (bypass) in a graphic feature information calculating unit.

FIG. 5 is a diagram showing an example of graphic characteristic information calculation in a process from an operation stability to an operation abnormality (bypass) in the graphic characteristic information calculation unit.

FIG. 6 is a diagram illustrating an example of graphic characteristic information calculation in a process from stable operation to abnormal operation (bypass) in the graphic characteristic information calculation unit.

FIG. 7 is a diagram showing an example of graphic characteristic information calculation in a process from operation stability to operation abnormality (bypass) in the graphic characteristic information calculation unit.

FIG. 8 is a diagram illustrating an operation prediction method using a temporal transition of graphic feature information in an operation prediction unit.

[Explanation of symbols]

 REFERENCE SIGNS LIST 1 blast furnace equipment 2 multiple sensors on blast furnace equipment 3 data collection device 4 data processing device 5 isoline calculation unit 6 graphic feature information calculation unit 7 operation monitoring unit 8 graphic feature information transition calculation unit 9 operation prediction unit 10 recording unit 11 Output section

Claims (11)

[Claims] 1. A measurement data of an amount to be measured from a plurality of sensors installed in a blast furnace is arranged in a two-dimensional plane or a three-dimensional space reflecting an installation position of each sensor, and a distribution state of each measurement data is determined. A method of monitoring the operating state of a blast furnace by expressing the information as a graphic or characteristic information of the graphic formed by these, and evaluating the graphic information, the arbitrary amount having the same value in the two-dimensional plane or the three-dimensional space. An operation monitoring method in a blast furnace operation, comprising calculating an isoline and evaluating a graphic or characteristic information of the graphic formed by the isoline. 2. A method for arranging measurement data of measurement quantities from a plurality of sensors installed in a blast furnace in a two-dimensional plane or a three-dimensional space reflecting an installation position of each sensor. A method for monitoring the operating state of a blast furnace by expressing the information as a graphic or characteristic information of the graphic formed by the evaluation and evaluating the information, wherein the arbitrary amount of the measurement target is the same on the two-dimensional plane or the three-dimensional space. Calculate isolines, and calculate at least one feature information of the number, position, area, center of gravity, aspect ratio of the figure, maximum or minimum value in the figure, average value, and variance of the figure formed by the isolines. An operation monitoring method in blast furnace operation, wherein the method is calculated and evaluated by image processing. 3. A method of updating a figure formed by isolines or feature information of the figure corresponding to a temporal transition of each measurement data, and comparing these temporal transitions with a preset transition condition. The operation monitoring method for blast furnace operation according to claim 1. 4. The number, position, area, center of gravity, aspect ratio of the figure, maximum or minimum value in the figure, average value, the number of figures formed by isolines corresponding to the temporal transition of each measurement data, 3. The operation monitoring method in blast furnace operation according to claim 2, wherein a temporal transition of at least one characteristic information of the variance is calculated by image processing, and the temporal transition is compared with a preset transition condition. 5. A method of calculating an isoline, comprising selecting, for data arranged in an uneven positional relationship on a two-dimensional plane, a quadrilateral element whose one of the internal angles does not exceed 180 degrees, 5. A method of setting an average value of data of four vertices at an intersection of, and searching for and drawing isolines using a triangle element having the intersection as a vertex. 2. The operation monitoring method for blast furnace operation according to item 1. 6. A data collection unit for collecting measurement data measured by a plurality of sensors installed on a plurality of blast furnace facilities, and a distribution state of the collected measurement data reflecting an installation position of each sensor on the blast furnace facility. A contour line calculating unit for arranging the contour lines in a two-dimensional plane or a three-dimensional space, and calculating an arbitrary contour line having the same measurement target amount, and a graphic or feature information of the graphic formed by the contour lines by image processing A graphic feature information calculating unit to be calculated; and an operation monitoring unit that monitors the operation by comparing the graphic or graphic characteristic information obtained by the graphic characteristic information calculating unit with the preset graphic or graphic characteristic information. An operation monitoring device for blast furnace operation, characterized in that: 7. The calculation of the contour line calculation unit and the graphic feature information calculation unit is repeated in accordance with the temporal transition of each measurement data,
A graphic feature information transition calculating unit that calculates these temporal transitions, and an operation that predicts an operating abnormality such as a blow-by by comparing a transition condition of the graphic or graphic characteristic information in which the graphic characteristic information transition information is set in advance. The operation monitoring device for blast furnace operation according to claim 6, further comprising a prediction unit. 8. The isoline calculating unit selects a quadrilateral element whose one of the internal angles does not exceed 180 degrees with respect to data arranged in a non-uniform positional relationship on a two-dimensional plane, 8. The blast furnace according to claim 6, further comprising a function of setting an average value of data of four vertices at the intersection, and searching for and drawing isolines using a triangle element having the intersection as a vertex. Operation monitoring device in operation. 9. An output unit for visualizing the transition of the figure and characteristic information of the figure.
The operation monitoring device in blast furnace operation according to any one of the above. 10. A computer-readable recording medium on which a program for causing a computer to execute the processing procedure of the operation monitoring method according to claim 1 is recorded. 11. A computer-readable recording medium on which a program for causing a computer to function as each unit of the operation monitoring device according to claim 6 is recorded.