JP2010150583A - Method for measuring layer thickness distribution of charged material in blast furnace, and apparatus for measuring layer thickness distribution using the same - Google Patents

Abstract

<P>PROBLEM TO BE SOLVED: To provide a method for measuring layer thickness distribution of charged materials in a blast furnace, with which the layer thickness distribution along the furnace radial direction from the furnace center to the furnace wall of the charged materials can be precisely obtained in a short period of time, and an apparatus for measuring the layer thickness distribution of charged materials in the blast furnace using the method. <P>SOLUTION: The method includes: a measuring step; a shaping line deriving step; a correcting step; and a layer thickness distribution deriving step. In the shaping line deriving step, a continuous function F1 for regulating an estimated shaping line 50 is set in advance. This continuous function F1 is a plurality of functions continued in the furnace radial direction, and regulated with a first linear function f11, a first curving function f21, a second linear function f12, and a second curving function f22, lined up from the furnace wall toward the furnace center side. The first primary function f11 is a function where a line regulated by the function is positioned in the vicinity of the furnace wall and is a straight line in the horizontal direction. <P>COPYRIGHT: (C)2010,JPO&INPIT

Description

  The present invention relates to a method for measuring the layer thickness distribution along the furnace radial direction from the furnace center to the furnace wall of the charge such as ore and coke charged in the blast furnace, and the layer thickness distribution using this method. It relates to a measuring device.

  In order to stably operate the blast furnace, along the radial direction of the furnace from the furnace center to the furnace wall of the layer (ore layer and coke layer) formed by the charge charged and deposited in the furnace Adjustment of the layer thickness distribution is important.

  As a method for adjusting the layer thickness distribution, the mover armor is operated in the bell-arm type blast furnace, or the turning chute is operated in the bell-less type blast furnace, and the charging position of the charges such as ore and coke is changed in the furnace. It is adjusted by. In this adjustment method, the layer thickness distribution of the charge is accurately grasped, and based on this grasped layer thickness distribution, the layer thickness distribution of the layer formed by the next charge becomes a more suitable distribution. In addition, operation of the charging device such as the movable armor and the turning chute is performed.

  In order to grasp the layer thickness distribution of the charge, the shape of the layer top surface of the charge in the corresponding batch (unit of charging of the charge into the furnace), that is, the surface profile and the previous one Use the surface profile of the batch and the total volume of the charge in that batch. Specifically, the volume of the portion between the surface profile of the corresponding batch and the surface profile of the previous batch is the same as the volume of the entire charge of the corresponding batch. Set the distance between the top and bottom surface profiles with the surface profile down. The layer thickness distribution is estimated (measured) by imitating the vertical width between the two upper and lower surface profiles set in this way as the layer thickness.

  In order to obtain the above surface profile, for example, a method of directly measuring the shape of the upper surface of the charge layer with a microwave profile meter 12 as shown in FIG. 2 is known. The microwave profile meter 12 measures the height of the surface profile (depth from the microwave irradiation position) at a constant pitch along the furnace radius at a predetermined position by utilizing the reflection of the microwave.

  Specifically, each measured value obtained by measuring along the furnace radius at the predetermined position on the upper surface of the charge layer in a predetermined batch, the furnace central axis as the Y axis, and the furnace radius at the predetermined position. Is plotted as the X axis. By connecting the plotted points, the top layer shape of the charge along the radius of the furnace at the predetermined position can be obtained. By rotating this shape once so that the furnace center axis (Y axis) is the center of rotation, the layer top surface shape of the charge of the batch, that is, the surface profile is obtained.

As another method, the relationship between the charging conditions and the surface profile is determined from the past performance data, and various charging conditions (for example, armature stroke, tilt angle of the chute, charging weight of charging material, and souding) There is known a method using a simulation model that calculates a surface profile using a level or the like as a calculation condition and estimates a layer thickness distribution based on the surface profile (see Patent Document 1).
JP 2001-323306 A

  However, when the surface profile is obtained using the microwave profile meter 12, an error occurs between the surface profile obtained from the measured value and the actual surface profile. In other words, each measurement value obtained by the measurement is plotted as described above, and the line obtained by connecting the plotted points often has a severe uneven shape different from the actual surface profile. . This is because a partial protrusion or the like that occurs when the charge is charged into the furnace is present in the measured layer top surface shape, and this protrusion or the like is included in the measured value as noise. Therefore, even if the top surface shape of the charge along the furnace radius at the predetermined position obtained by the plot is rotated around the furnace center axis, the surface profile is different from the actual surface profile. There are many.

  As described above, in the case of estimating the lowering amount of the lower surface profile so that the volume of the space sandwiched between the upper and lower surface profiles is the same as the volume of the entire charge of the corresponding batch, When the error generated on the wall side is large, the calculated interval between the profiles is greatly different, and an accurate layer thickness distribution cannot be calculated.

  In order to eliminate the error due to the partial protrusion or the like, it is conceivable that several tens of batches of the surface profile are measured, and the surface profile is calculated using a value obtained by averaging the measured values of several tens of times.

  However, in such a method, it is necessary to calculate the surface layer and calculate the latest layer thickness after obtaining the average value of every tens of batches measured in the past. When there is a change, the change cannot be grasped quickly and the operating conditions cannot be changed in response to the change.

  On the other hand, since the simulation model is a collection of various calculation conditions, in order to adjust the parameters of these calculation conditions for each operation condition, the relationship between many measurement profiles and each operation condition must be analyzed. It takes a very long time.

  Even if the parameters thus obtained are used, when the surface profile is calculated using the simulation model, only the same result can be calculated under the same operating conditions. However, the surface profile of the charge formed in an actual blast furnace may not be the same surface profile in spite of the same operating conditions, and may have a different surface profile for each batch. Many.

  As described above, when the simulation model is used, it takes a very long time to adjust the parameter of the calculation condition that matches the operation condition even if the operation condition is partially different. Moreover, the problem that it cannot respond to the change of the layer thickness distribution under the same operation condition arises.

  Therefore, in view of the above problems, the present invention provides a layer thickness distribution of a blast furnace charge that can accurately obtain a layer thickness distribution along the furnace radial direction from the furnace center to the furnace wall of the charge. It is an object of the present invention to provide a measuring method and an apparatus for measuring a layer thickness distribution of a blast furnace charge using this method.

  The inventors of the present invention have made extensive studies to solve the above problems, and as a result, have found the following.

  Normally, when the charge is charged into the blast furnace during the operation of the blast furnace, it is necessary to adjust the surface profile to maintain stable operation of the blast furnace. The charge is charged while changing the dropping position along the center, and is charged (center charge) so that the charged charge rises in the center of the furnace. When analyzing the data of many surface profiles such as past performance data in such an operating state, the surface profile at the furnace radius at a predetermined position has a curved line in any surface profile in the operating blast furnace. When simulated by the (estimated shape line), a horizontal or substantially horizontal straight line portion near the furnace wall and a straight down slope that is located on the furnace center side of the straight line portion toward the furnace center. And a curved portion that bulges upward in the furnace center.

  Therefore, based on this discovery, the inventors focused on an estimated shape line having a straight portion and a curved portion as a part thereof, and used the method for measuring the layer thickness distribution of the blast furnace charge having the following configuration and this method. Created a measuring device.

  The method for measuring the thickness of a blast furnace charge according to the present invention involves repeating the batch for charging the charge into the blast furnace n times (where n is a natural number arbitrarily selected within the range of n> 1). A method for measuring a layer thickness distribution along a furnace radial direction from a furnace center to a furnace wall in a layer formed by a charge of an n-th batch among charges stacked in a blast furnace, A measurement step of measuring the layer upper surface shape of the charge at a furnace radius at a predetermined position in each batch of the (n-1) th batch and the nth batch, and an estimated shape that simulates the layer upper surface shape as a line having a curve A shape line deriving step for deriving a line for each batch of the (n−1) -th batch and the n-th batch based on the measured value of each batch obtained in the measuring step, and deriving in the shape line deriving step 2 estimated shape lines in the furnace When the estimated surface profiles are obtained by rotating the shaft so that the axis is the center of rotation, the volume of the space sandwiched between the two estimated surface profiles is equal to the total volume of the charge in the nth batch. A correction step of correcting at least one of the two estimated shape lines so that the estimated shape line of the n-1st batch is positioned below the estimated shape line of the n-th batch; The furnace wall from the furnace center in the charge of the nth batch based on the vertical spacing at each position along the radial direction of the furnace in each estimated shape line of the n-1st batch and the nth batch A layer thickness distribution deriving step for deriving a layer thickness distribution along the furnace radial direction up to, and in the shape line deriving step, a first continuous function defining the estimated shape line is set in advance, The continuous function 1 is a plurality of functions that are continuous in the furnace radial direction at the predetermined position, and the coefficient of each function and the range in the furnace radial direction at the predetermined position are undetermined, and adjacent functions are continuous at the boundary. As defined by a first linear function, a first curve function, a second linear function, and a second curve function arranged in order from the furnace wall to the furnace center side, A line defined by the function is located in the vicinity of the furnace wall and becomes a horizontal or substantially horizontal straight line, and the second linear function is a line defined by the function from the first linear function. Is a function which is located on the furnace center side and becomes a straight line having a constant downward gradient toward the furnace center, and the first curve function is such that a line defined by the function is the first linear function and the first function. Between the first linear function and the first function A function that is a curve that smoothly connects a furnace center side end portion of a substantially horizontal straight line defined by a linear function and a furnace wall side end portion of a straight line of a downward slope defined by the second linear function, The second curve function is a function in which a line defined by the function is located at the center of the furnace and becomes a curve that bulges upward. In each batch, the measured value of the batch obtained in the measurement step The estimated shape line for each batch of the n-1st batch and the nth batch is derived by calculating the coefficient in each function and the range in the furnace radial direction of the predetermined position on a batch basis. It is characterized by doing. Here, in the present invention, the line means a finite length line having both ends.

  According to such a configuration, the layer upper surface shape (hereinafter also simply referred to as “surface profile”) along the furnace radius at the predetermined position in two consecutive batches (the n−1th batch and the nth batch) is measured. By simply doing, it is possible to measure the layer thickness distribution of the charge in the batch in which the charge was charged (the nth batch). Therefore, it is not necessary to measure the surface profile of several tens batches and obtain the average of the measured values as in the past, and the batch in which the top layer charge is charged in a short time (time required for two consecutive batches). The layer thickness distribution of the (nth batch) can be measured. Moreover, in order to define an estimated shape line that simulates the surface profile, the radial direction of the blast furnace (furnace radial direction) is obtained by using the first continuous function in which the portion near the furnace wall is a straight line (first linear function). ), It is possible to derive an estimated shape line that conforms to the actual surface profile even in the operation state in which the charge is charged to the furnace wall side while changing the fall position along this line. Can be measured.

  Specifically, a region (so-called terrace) that is horizontal or substantially horizontal in the vicinity of the furnace wall in the surface profile by charging the furnace wall side while changing the dropping position along the furnace radial direction. Part) is formed. Since the top surface shape of the terrace portion is a horizontal or substantially horizontal linear shape at the furnace radius at the predetermined position, in the first continuous function in which a plurality of functions are connected, the first linear function is provided in the section closest to the furnace wall. Is used, the shape of the terrace portion is imitated with high accuracy. As described above, the portion near the furnace wall in the estimated shape line simulates the surface profile with high accuracy, so that the layer thickness distribution can be measured with high accuracy.

  In other words, since the estimated surface profile is obtained by rotating the estimated shape line around the furnace center axis as the center of rotation, the moving distance of the part where the error occurs when the estimated shape line is rotated toward the furnace wall side becomes longer. As a result, errors on the furnace wall side have a greater effect on the estimated surface profile. Therefore, in the correction step, the estimated shape line in which an error occurs on the furnace wall side as described above is used, and the volume of the space sandwiched between the two upper and lower estimated surface profiles is the total of the charge of the nth batch. When correction is made so that the estimated shape line of the (n-1) th batch is positioned below the estimated shape line of the nth batch so as to be the same as the volume, the influence of the error appears greatly, and the layer thickness distribution The measurement accuracy decreases. Therefore, in the estimated shape line derived as described above, the portion in the vicinity of the furnace wall conforms to the shape of the terrace portion of the actual surface profile, that is, has a shape with less error, so that two consecutive batches (the first batch) The shape of each estimated surface profile of (n-1th batch and nth batch) and the distance between these two estimated surface profiles can be accurately derived.

  Further, the first continuous function can accurately simulate the surface profile not only in the terrace portion of the surface profile but also in the portion closer to the furnace center than the terrace portion at the furnace radius at a predetermined position. Specifically, a portion of the surface profile at the furnace radius at a predetermined position located on the furnace center side from the terrace portion and having a constant downward gradient toward the furnace center is a second linear function that is a straight line of the downward gradient. The second curve function, which is a curve that bulges upward in the center of the furnace and bulges upward from the center charging, can be simulated with high accuracy. Therefore, by using the first continuous function, it is possible to obtain an estimated shape line corresponding to the actual surface profile in the entire range from the furnace wall to the furnace center at the furnace radius at a predetermined position.

  Further, since the estimated shape line is used to derive the estimated surface profile, the layer thickness distribution can be easily and accurately calculated based on these measured values even if noise is included in the measured values obtained by measuring the surface profile. Can be derived.

  Specifically, an estimated shape line in the furnace radius at the predetermined position is defined, and a first continuous function in which a plurality of functions are connected is set in advance, and the coefficient of each function and the predetermined position of the predetermined position are set based on the plurality of measured values. By calculating the range of the furnace radius, an estimated shape line as a line having a curve is derived regardless of the magnitude and amount of noise included in each measurement value. The surface profile at the furnace radius at the predetermined position is simulated by this estimated shape line, so that even if there is a partial protrusion or the like (noise) that happens to occur in the measured top surface shape of the layer, A highly accurate surface profile can be obtained without being affected by protrusions or the like.

  Further, in the most recent batch (nth batch) and the previous batch (n-1th batch), the surface profile at the furnace radius at the predetermined position is actually measured, and the measured value (actually measured value). ) Is calculated in advance, the coefficient of each function that defines the first continuous function set in advance and the range of the furnace radius at the predetermined position are derived, so that an accurate estimated shape line corresponding to the surface profile is derived. Thus, the measurement of the layer thickness distribution corresponding to the change in the layer thickness distribution under the same operating condition of the blast furnace is performed with high accuracy.

  In addition, in order to solve the above-described problem, the method for measuring the layer thickness of the blast furnace charge according to the present invention charges the charge to the furnace wall side while changing the drop position along the furnace radial direction for each batch. In the blast furnace operation where the load is charged to the furnace wall side while keeping the drop position in the furnace radial direction at a fixed position, charging to the furnace wall side is performed while changing the drop position in the furnace radial direction. The continuous function used when deriving the estimated shape line may be changed depending on whether or not an object is inserted. That is, in the shape line derivation step, a second continuous function having an undetermined coefficient that defines the estimated shape line is further set, and the second continuous function includes a plurality of continuous functions in the furnace radial direction of the predetermined position. The function coefficient and the range of the predetermined position in the furnace radial direction are undetermined, and adjacent functions are continuous at the boundary from the furnace wall to the furnace center side along the furnace radius. Are defined by a third curve function, a third linear function, and a fourth curve function, which are arranged in order, and the line defined by the function is located near the furnace wall, The fourth curve function is a function that becomes a curve that bulges upward when a line defined by the function is located in the center of the furnace, and the fourth curve function is a function that becomes a curve that bulges upward. The linear function is defined by the function Is a function that is located between the third curve function and the fourth curve function and forms a straight line with a constant downward slope toward the furnace center, and the n-1st batch and the nth In each batch with the batch, when charging the charging material into the blast furnace while changing the dropping position along the furnace radial direction when charging the charging material to the furnace wall side, the measurement is performed. Deriving the estimated shape line of the batch by calculating the coefficient in each function defining the first continuous function based on the measured value of the batch obtained in the step and the range in the furnace radial direction of the predetermined position, When charging the charging material into the blast furnace without changing the dropping position in the furnace radial direction when charging the charging material to the furnace wall side, the measured value of the batch obtained in the measuring step And a coefficient in each function defining the second continuous function based on Preferably, to derive the estimated shape line of the batch by calculating the furnace radial range of positions.

  By being configured in this way, even in a batch in which a charge is charged to the furnace wall side without changing the dropping position in the furnace radial direction, the second continuous function is based on the measurement value obtained in the measurement step. It is possible to obtain an estimated shape line having a shape conforming to the surface profile of the corresponding batch by deriving the coefficient of each function that defines the above and the range in the furnace radius of the predetermined position.

  Specifically, by charging without changing the fall position in the furnace radial direction, the charge is shaped like a mountain with the fall position at the top at the furnace radius at the predetermined position. The terrace is difficult to form. Therefore, the furnace wall in this surface profile is obtained by using the second continuous function which is a shape in which the portion in the vicinity of the furnace wall is raised like a mountain (curve bulging upward (third curve function)). It is possible to accurately simulate the shape of a nearby part (a shape that rises like a mountain). That is, when the terrace portion is not formed, the first continuous function in which the vicinity of the furnace wall is a straight line is not used, but the second continuous function in which the vicinity of the furnace wall is a curve bulging upward is used. The shape of the surface profile can be simulated with high accuracy. Therefore, even in a batch in which the charge is charged to the furnace wall side without changing the fall position in the furnace radial direction, the surface profile in the vicinity of the furnace wall approximated a surface profile having a mountain-like area. An estimated shape line can be derived, whereby the layer thickness distribution can be accurately measured.

  In addition, the third curve function, the first curve function, the first curve function, the second curve function, and the continuous function defined by the second curve function (the first continuous function) are used. 3 using the linear function of 3 and the continuous function defined by the fourth curve function (second continuous function), the number of functions defining the continuous function is small, and the calculation is simplified and facilitated. Compared with the case where the estimated shape line is derived using the first continuous function regardless of whether or not the charge is charged to the furnace wall side while changing the fall position in the furnace radial direction, it is easier and shorter. Thus, it is possible to obtain an estimated shape line having a shape conforming to the actual surface profile.

  The first continuous function or the second continuous function is preferably a function represented on an xy plane in which a furnace center axis is a y-axis and a furnace radius at the predetermined position is an x-axis. By doing in this way, each function which prescribes | regulates a continuous function is represented in the form of y = f (x), and handling of each said function becomes easy.

The first curve function, the second curve function, the third curve function, or the fourth curve function may be expressed as x at an angle formed by a tangent to a curve defined by the curve function and the x axis. It is preferable that the rate of change in angle with respect to the axial position is a constant function. Specifically, the first curve function, the second curve function, the third curve function, or the fourth curve function is a function represented by the following expression (1) or (2): It is preferable that
y = (− 1 / a) · log | cos (ax + b) | + c (1)
y = αx ² + βx + γ (2)
Here, a, b, c, α, β, and γ are coefficients.

  Thus, since the first curve function, the second curve function, the third curve function, or the fourth curve function is a simple function and has few coefficients, the coefficient of each function based on the measurement value It is also easier to calculate the range in the x-axis direction, that is, the range of the furnace radius at the predetermined position. Therefore, it is easy to calculate the coefficient of each function and the range in the x-axis direction based on the measured value of the surface profile, and it is easy to derive the estimated shape line.

  Further, in the shape line derivation step, coefficients in all functions that define the first continuous function or the second continuous function using the steepest descent method from the measurement values obtained in the measurement step, and the x-axis direction By simultaneously obtaining the ranges in the above, it is possible to easily obtain the coefficients of all the functions that define the continuous function and the ranges in the x-axis direction easily from the measurement values obtained in the measurement step.

  In order to solve the above-mentioned problem, the apparatus for measuring the layer thickness of a blast furnace charge according to the present invention is configured to perform batch charging of the charge into the blast furnace n times (where n is in the range of n> 1). Layer thickness distribution along the furnace radial direction from the furnace center to the furnace wall in the layer formed by the charge of the nth batch among the charge stacked in the blast furnace by repeating. In accordance with the measurement value obtained by measuring the top surface shape of the charge at the furnace radius at a predetermined position in each of the (n-1) th batch and the nth batch. Layer thickness distribution deriving means for deriving the layer thickness distribution along the furnace radial direction from the furnace center to the furnace wall in the nth batch charge, and the layer thickness distribution derived by this layer thickness distribution deriving means to the outside Output means for outputting, the layer thickness distribution deriving means is Defines an estimated shape line that simulates the shape of the upper surface of the layer as a line having a curve, a plurality of functions are coefficients in each function and the range in the furnace radial direction of the predetermined position is undetermined, and adjacent functions are continuous at the boundary And a function storage unit that stores in advance a first continuous function defined by being connected along the furnace radial direction of the predetermined position, and a charge charged into the blast furnace in the nth batch A volume storage unit that stores the value of the entire volume, and the function storage unit that is stored in the function storage unit for each batch based on measured values of each batch of the (n-1) th batch and the nth batch. Deriving estimated shape lines for each batch of the (n-1) th batch and the nth batch by calculating a coefficient in each function defining a continuous function of 1 and a range in the furnace radial direction of the predetermined position Shi An estimated shape line deriving unit that stores the values of these two derived estimated shape lines, and two estimated shape lines whose values are stored in the estimated shape line deriving unit When the estimated surface profiles are respectively obtained by rotating the same, the volume of the space sandwiched between the two estimated surface profiles is equal to the volume of the entire charge of the nth batch stored in the volume storage unit. A position correction unit that corrects at least one of the two estimated shape lines so that the estimated shape line of the (n-1) th batch is positioned below the estimated shape line of the nth batch, In the charge of the nth batch based on the vertical spacing at each position along the furnace radial direction in each estimated shape line of the n-1st batch and the nth batch corrected in this way. Furnace wall from the furnace center A layer thickness distribution deriving unit for deriving a layer thickness distribution along the furnace radial direction and transmitting a value of the derived layer thickness distribution to the output means, and stored in advance in the function storage unit The first continuous function is a first linear function, a first curve function, a second linear function, and a second linear function that are arranged in order from the furnace wall toward the furnace center along the furnace radius at the predetermined position. The first linear function is a function in which a line defined by the function is located near the furnace wall and becomes a horizontal or substantially horizontal straight line, and the second linear function is The line defined by the function is a function that is located closer to the furnace center than the first linear function and becomes a straight line with a constant downward slope toward the furnace center, and the first curve function is the function A line defined by the first linear function and the second linear function The furnace center side end of the horizontal or substantially horizontal straight line defined by the first linear function and the furnace wall side end of the downwardly straight line defined by the second primary relation. The second curve function is a function that becomes a curve in which a line defined by the function is located in the center of the furnace and bulges upward.

  According to such a configuration, the furnace radius at a predetermined position in the surface profile of the last two consecutive batches (the n-1st batch and the nth batch) is measured by a measuring means that is practically used in a blast furnace that is already in operation. By obtaining the measured values at a plurality of locations along the line, it is possible to accurately measure the layer thickness distribution of the charge in the batch in which the charge in the top layer is charged (the nth batch). And since the site | part of the furnace wall vicinity of the estimated shape line which the 1st continuous function previously stored in the function memory | storage part prescribes | regulates is a linear shape (1st linear function), a fall position along a furnace radial direction It is possible to derive an estimated shape line in accordance with the surface profile having the terrace portion even in the operation state in which the charge is charged to the furnace wall side while changing the thickness, thereby accurately measuring the layer thickness distribution. Can do.

  In addition, by using the estimated shape line, the layer thickness distribution can be easily and accurately derived even if noise is included in the measurement value, while the derived layer thickness distribution is output by the output means. This makes it possible for an operator operating the blast furnace to quickly and accurately grasp the layer thickness distribution of the batch in which the top layer charge is charged.

  In addition, the apparatus for measuring the layer thickness distribution of a blast furnace charge according to the present invention is a method for determining whether or not the charge is charged to the furnace wall side while changing the radial drop position in the n-th batch. At least one of an input unit for inputting state information or a receiving unit to which the charging state information is transmitted from the outside, and a state storage unit for storing the input or transmitted charging state information; The function storage unit defines the estimated shape line, and the predetermined function is such that a plurality of functions have a coefficient in each function and a range in the furnace radial direction of the predetermined position and adjacent functions are continuous at the boundary. A second continuous function that continues along the furnace radial direction of the position is further stored, and the second continuous function stored in the function storage unit is stored in the furnace wall along the furnace radius at the predetermined position. The third lined up in order toward the center It is defined by a curve function, a third linear function, and a fourth curve function, and the third curve function is a curve in which a line defined by the function is located near the furnace wall and bulges upward. The fourth curve function is a function in which a line defined by the function is located in the furnace center and becomes a curve that bulges upward, and the third linear function is determined by the function The defined line is a function that is located between the third curve function and the fourth curve function and becomes a straight line having a constant downward slope toward the furnace center. In the estimated shape line deriving unit, Based on the charging state information stored in the state storage unit, in each batch of the (n-1) th batch and the nth batch, when charging the charging material to the furnace wall side Charge the charge into the blast furnace while changing the drop position along the furnace radial direction. When calculating the coefficient in each function defining the first continuous function stored in the function storage unit based on the measurement value of the batch and the range in the furnace radial direction of the predetermined position, When charging the charging material into the blast furnace without changing the dropping position in the furnace radial direction when charging the charging material to the furnace wall side, the function storage is performed based on the measured value of the batch. By calculating the coefficient in each function defining the second continuous function stored in the section and the range in the furnace radial direction of the predetermined position, for each batch of the n-1th batch and the nth batch The estimated shape line may be derived, and the values of the two derived estimated shape lines may be stored.

  According to such a configuration, whether or not the charge to the furnace wall side is charged while changing the dropping position of the furnace radius is input by the operator operating the blast furnace or transmitted from the control unit of the blast furnace or the like. Thus, the continuous function used in the estimated shape line deriving unit is switched so that an estimated shape line having a shape corresponding to the surface profile can be derived in a shorter time.

  In the blast furnace charge layer thickness distribution measuring device, the top surface shape of the charge at the furnace radius at the predetermined position is measured, and the measured value obtained by this measurement is derived from the estimated shape line deriving means. Measurement means for transmitting to may be further provided.

  As described above, according to the present invention, the method for measuring the layer thickness distribution of the blast furnace charge capable of accurately obtaining the layer thickness distribution along the furnace radial direction from the furnace center to the furnace wall of the charge, And the apparatus for measuring the layer thickness distribution of the blast furnace charge using this method can be provided.

  Hereinafter, a first embodiment of the present invention will be described with reference to the accompanying drawings.

  The blast furnace charge layer thickness distribution measuring device (hereinafter also simply referred to as “measuring device”) is a batch for charging coke, ore, or other charge into the blast furnace (charging the blast furnace into the blast furnace). Among the charges stacked in the blast furnace by repeating the charging unit) N times (where N is a natural number and 1 <N) from the furnace center in the layer formed by the charge of the Nth batch. It is an apparatus that measures the layer thickness distribution along the furnace radial direction to the wall. This layer thickness distribution is used to maintain the blast furnace operation in a stable manner by maintaining the deposited layer of the charge in an appropriate shape.

  Specifically, as shown in FIG. 1, the measuring apparatus includes a measuring unit 12, a layer thickness distribution deriving unit 20, and an output unit 14.

  The measuring unit 12 measures a surface profile at a furnace radius at a predetermined position, and transmits (transmits) a measurement value obtained by this measurement to the layer thickness distribution deriving unit 20 by an output signal. A so-called microwave profile meter is used. As shown in FIG. 2, the measuring means 12 has a measuring rod 12a inserted from the outside of the furnace so as to penetrate the furnace wall, and a measuring portion 12b is provided at the tip thereof. The measuring portion 12b at the tip of the measuring rod 12a reciprocates along the furnace radius at a predetermined position of the furnace, and the charge layer is based on the reflected wave of the microwave irradiated from the tip (measuring portion) 12b. The height of the upper surface (profile depth) is measured. The measuring means 12 is provided in the measuring device 10 in the present embodiment. For example, when a microwave profile meter is already provided in the blast furnace, this is measured by the measuring means 12 of the measuring device 10. It may be used as

  The layer thickness distribution deriving unit 20 is configured to charge the Nth batch based on the measurement values obtained by measuring the surface profile in the measuring unit 12 in each of the (N-1) th batch and the Nth batch. Of the layer thickness in the furnace radial direction (layer thickness distribution), a function storage unit 21, a volume storage unit 22, an estimated shape line deriving unit 23, a position correcting unit 24, a layer thickness distribution deriving unit 25, Have

  The function storage unit 21 stores the first continuous function F1. The first continuous function F1 is for defining an estimated shape line 50 for imitating the surface profile at the furnace radius at a predetermined position, and this estimated shape line 50 is defined as a line having a curve. In the present embodiment, the above-mentioned line that defines the estimated shape line 50 refers to a line having a straight line portion and a curved portion and having a finite length having both ends.

  The first continuous function F1 is obtained as a result of analyzing a number of surface profile data such as past actual data, and the actual data includes a drop in the furnace radial direction in the vicinity of the furnace wall. The charge is charged while changing the position, and in the center of the furnace, the premise is the operation state in which the charge is centrally charged so that the charged charge is elevated. Specifically, in any blast furnace in the above operating state, if the surface profile at the furnace radius at a predetermined position is simulated by a line having a curve (estimated shape line) in any surface profile, , A horizontal or substantially horizontal straight part near the furnace wall, an inclined part located on the furnace center side of the straight part and having a constant downward slope toward the furnace center, and bulging upward in the furnace center It is a function set by paying attention to the fact that a curved portion tends to appear.

  The first continuous function F1 set in this way is a function represented on the XY plane with the furnace center axis as the Y axis and the furnace radial direction at a predetermined position as the X axis, and this first continuous function F1. As shown in FIG. 3, the estimated shape line 50 defined by the straight line portion 52, the inclined portion 54, the connection curve portion 53 that smoothly connects the straight portion 52 and the inclined portion 54, and the central curve portion 55. Have Here, the straight line portion 52 is a portion that is located in the vicinity of the furnace wall of the line (estimated shape line) 50 defined by the first continuous function F1, and is a horizontal or substantially horizontal straight line, and is inclined. The part 54 is a part that is located closer to the furnace center than the straight part 52 among the lines defined by the first continuous function F1 and is a straight line having a constant downward slope toward the furnace center, and is a connection curve. The portion 53 is located between the straight portion 52 and the inclined portion 54 of the first continuous function F1, and the furnace center side end portion (left end portion in FIG. 3) of the straight portion 52 and the inclined portion 54. Is the portion that forms a curve that smoothly connects the furnace wall side end portion (the right end portion in FIG. 3), and the central curve portion 55 is located in the furnace center portion of the first continuous function F1. This is a part that forms a curve that bulges upward. The gradient and length of the straight portion 52, the gradient and length of the inclined portion 54, the shape of the connection curve portion 53, the shape of the central curve portion 55, and the like are not determined because they are different for each surface profile. 12 is used to calculate the coefficient of the first continuous function F1 based on a plurality of measurement values obtained by measuring the surface profile, so that the gradient and length of the linear portion 52 and the gradient and length of the inclined portion 54 are calculated. Then, the shape of the connection curve portion 53, the shape of the center curve portion 55, and the like are respectively determined, and the estimated shape line 50 that approximates the actual surface profile can be obtained.

  Specifically, the first continuous function F1 is defined by a plurality of functions that are continuous along the X-axis direction. The plurality of functions are set so that the coefficient and the range in the X-axis direction are undetermined and are continuous at the boundary between adjacent functions. In the present embodiment, the first linear function f11, the first curve function f21, the second linear function f12, and the second curve function f22 are arranged in order from the furnace wall toward the furnace center along the X axis. It is prescribed by. Each function excluding the second curve function f22, that is, the first linear function f11, the first curve function f21, and the second linear function f12 is differentiable at each of the boundary positions r1 and r2. Functions are connected.

  The first linear function f11 is a function in which a line defined by the function f11 is a straight line corresponding to the straight line portion 52, and the first curve function f21 is a function in which the line defined by the function f21 is a connection curve portion. 53, a second linear function f12 is a function in which a line defined by the function f12 is a straight line corresponding to the inclined portion 54, and a second curve function f22 is This is a function in which the line defined by the function f22 becomes a curve corresponding to the central curve portion 55.

If each function f11, f12, f21, f22 is shown more concretely, the first linear function f11 is
y = e ₁ x + f ₁ (hereinafter, also simply referred to as “Expression (10)”)
And the second linear function f12 is
y = e ₂ x + f ₂ (hereinafter also simply referred to as “Expression (12)”)
It is represented by The first curve function f21 is
y = (− 1 / a) · log | cos (ax + b) | + c (hereinafter also simply referred to as “expression (11)”)
And the second curve function f22 is
y = αx ² + βx + γ (hereinafter also simply referred to as “Expression (13)”)
It is represented by Note that a, b, c, e ₁ , e ₂ , f ₁ , f ₂ , α, β, and γ are coefficients.

  In the following, in the furnace radius (X axis) at the predetermined position, the boundary position between the straight portion 52 and the connection curve portion 53 is r1, the boundary position between the connection curve portion 53 and the inclination portion 54 is r2, and the inclination portion 54. And a boundary position between the central curve portion 55 and r3.

  The volume storage unit 22 is a part for storing a value of volume (volume information) of the entire charged material charged into the blast furnace in each batch transmitted by an output signal from the control unit of the blast furnace.

The estimated shape line deriving unit 23 obtains the first continuous function F1 stored in the function storage unit 21, and performs the operation based on the measurement values of the batches of the (N-1) th batch and the Nth batch. Each of the coefficients of the continuous function F1 of 1 is calculated, and estimated shape lines 50 _n-1 and 50 _n (see FIG. 5B) for each batch of the (N-1) th batch and the Nth batch are derived, These are portions that store the values of the two estimated shape lines 50 _n−1 and 50 _n derived. The calculation of the coefficient of the first continuous function F1 in the present embodiment refers to the coefficients of the functions f11, f12, f21, and f22 that define the first continuous function F1 and the functions f11, f12, f21, and f22. It means obtaining the range in the X-axis direction.

The position correction unit 24 obtains the volume value of the entire charge of the Nth batch stored in the volume storage unit 22 from the volume storage unit 22 and, based on this volume, the N-1th batch. estimated shape line 50 _n-1 and a portion for correcting the distance between the estimated shape line 50 _n of the n-th batch. Specifically, the position correction unit 24 rotates the two estimated shape lines 50 _n-1 and 50 _n stored in the estimated shape line deriving unit 23 so that the furnace center axis is the center of rotation, and thereby the estimated surface profile. Are estimated below the estimated shape line 50 _n at intervals such that the volume of the space sandwiched between both estimated surface profiles is equal to the volume of the entire charge acquired from the volume storage unit 22. The estimated shape line 50 _n-1 is pulled down (corrected) so that the shape line 50 _n-1 is positioned. In the present embodiment, the position correcting unit 24, by pulling the estimated shape line _{50 n-1} is configured to adjust the distance between the two estimated shape line _{50 n-1, 50 n,} to The present invention is not limited, and it is sufficient that the interval is adjusted by correcting at least one of the two estimated shape lines 50 _n−1 and 50 _n .

The layer thickness distribution deriving unit 25 derives a layer thickness distribution along the furnace radial direction from the furnace center to the furnace wall in the charge of the Nth batch, and outputs the value of the derived layer thickness distribution by an output signal. The layer thickness distribution is a part that is transmitted (transmitted) to the output unit 14, and the layer thickness distribution is estimated shape lines 50 _n−1 , 50 of the N−1th batch and the Nth batch corrected by the position correction unit 24. It is derived based on the vertical interval at each position along the furnace radial direction at _n .

  The output means 14 is for receiving the output signal transmitted from the layer thickness distribution deriving unit 25 and outputting it to the outside. In the present embodiment, a CRT display is used as the output means 14. However, the present invention is not limited to this, and other display means such as an FPD may be used, or a printing means or the like may be used. A combination of these may also be used.

  The measuring apparatus 10 according to the present embodiment has the above-described configuration. Next, the operation of the measuring apparatus 10 will be described with reference to FIG.

  The surface profile of the (N-1) th batch at a predetermined furnace radius is measured by the measuring means 12, and the measurement value obtained by this measurement is transmitted as an output signal to the layer thickness distribution deriving means 20 (estimated shape line deriving section 23). (Step 1). In this embodiment, measurement is performed toward the furnace center at 10 cm intervals from the furnace wall along the furnace radius (X axis) at a predetermined position. At this time, since the furnace radius of the blast furnace is 3.675 m, the measured values (depth data) of the top surface of 38 layers at equal intervals on the X axis are obtained (FIG. 5A and FIG. 5). b)).

  The estimated shape line deriving unit 23 that has received the output signal applies the first continuous function F1 (see FIG. 3) stored in the function storage unit 21 to the measurement value obtained by the output signal, so that A coefficient of the first continuous function F1 is derived (step 2). At this time, the coefficients and the ranges in the X-axis direction in all the functions f11, f12, f21, and f22 that define the first continuous function F1 are calculated as follows based on the measured values.

The value of the coefficient (parameter) of the first continuous function F1 is determined so that the first continuous function F1 best matches the measured value. Specifically, the furnace radius position at an arbitrary measurement point is X _i , the actual surface profile measurement value at the furnace radius position is Y _i, and the first continuous function F1 at the furnace radius X is Y = F1 ( X)
Σ {Y _i −F 1 (X _i )} ²
Of the first continuous function F1 (coefficients e ₁ and f _{1 in} the above equation (10), coefficients a, b and c in the above equation (11), and the above equation (12) E ₂ , f ₂ , and the coefficients α, β, γ of the above equation (13) and the range of the function in the X-axis direction (boundary position r 1 between the straight line portion 52 and the connection curve portion 53, the connection curve portion) 53 and the boundary position r3 between the inclined portion 54 and the central curve portion 55) are simultaneously obtained by the steepest descent method which is one of the nonlinear optimization methods. Thus, by using the steepest descent method, the coefficients in all the functions f11, f12, f21, and f22 and the range in the X-axis direction in which each function is defined can be easily obtained.

The estimated shape line deriving unit 23 substitutes the coefficients of the functions f11, f12, f21, and f22 thus obtained and the ranges in the X-axis direction, that is, the first continuous function for which the coefficients are determined. F1 is stored as the estimated shape line 50 _n-1 of the (N-1) th batch.

In these steps 1 and 2, for example, the estimated shape line 50 _n-1 of the (N-1) th batch has already been derived and stored in the estimated shape line deriving unit 23, for example, when measuring the previous layer thickness distribution. If omitted, it is omitted. That is, in this case, the measurement of the layer thickness distribution starts from step 3 below.

  The surface profile of the Nth batch is measured by the measuring means 12 in the same manner as the N-1th batch, and the measured value obtained by this measurement is the layer thickness distribution deriving means 20 (the estimated shape line deriving section 23). (Step 3).

The estimated shape line deriving unit 23 that has received the output signal, like the N-1th batch, has a first continuous function F1 stored in the function storage unit 21 in the measurement value obtained by this output signal. By applying (see FIG. 3), the coefficient of the first continuous function F1 is derived, and the derived first continuous function F1 is stored as the estimated shape line 50 _n of the Nth batch (step) 4).

Thus estimated shape line of the N-1 th batch 50 _n-1 and the estimated shape line 50 _n of the N-th batch is derived, and stored respectively in the estimated shape-guiding portion 23, the position correcting unit 24 Acquires the value of the entire volume of the Nth batch of charge from the volume storage unit 22 and stores two estimated shape lines 50 _n−1 , 50 whose values are stored in the estimated shape line deriving unit 23. The interval of _n is corrected (step 5). Specifically, when the estimated surface profiles are obtained by rotating the two estimated shape lines 50 _n-1 and 50 _{n so} that the furnace center axis (Y axis) is the rotation center, both estimated surfaces are obtained. Nth below the estimated shape line 50 _n of the Nth batch at an interval such that the volume of the space sandwiched between the profiles is equal to the volume of the entire charge of the Nth batch acquired from the volume storage unit 22. The interval between the two estimated shape lines 50 _n-1 and 50 _n is corrected so that the estimated shape line 50 _{n-1 of the} first batch is located. In the present embodiment, the interval between the two estimated shape lines 50 _n-1 and 50 _n is corrected by lowering the estimated shape line 50 _n-1 of the (N-1) th batch as described above (FIG. 5).

When the position correction unit 24 corrects the interval between the two estimated shape lines 50 _n-1 and 50 _n in this way, each estimated shape of the N-1th batch and the Nth batch is then used. Based on the vertical distance at each position along the X axis in the lines 50 _n-1 and 50 _n , the layer thickness distribution deriving unit 25 performs the X axis (from the furnace center to the furnace wall in the charge of the Nth batch) A layer thickness distribution along the furnace radius direction is derived, and the value of the layer thickness distribution (layer thickness distribution information) is transmitted to the output means 14 as an output signal (step 6).

  The output unit 14 that has received this output signal displays the layer thickness distribution in the Nth batch as layer thickness distribution information (step 7). Based on this display, the operator who operates the blast furnace can accurately and accurately grasp the layer thickness distribution in the furnace radial direction of the charge of the Nth batch (uppermost layer).

In this way, the layer thickness distribution information is displayed by the output means 14, and the measurement of the layer thickness distribution of the Nth batch is completed. When the layer thickness distribution of the (N + 1) th batch is continuously measured, the estimated shape line 50 _n of the Nth batch is already derived and stored in the estimated shape line deriving unit 23, and thus the process starts from step 3 ( Step 8).

  According to the measuring apparatus 10 as described above, in the surface profile of the last two consecutive batches (the (N-1) th batch and the Nth batch) by the measuring means 12 or the like that has been put into practical use in an already operating blast furnace. By obtaining measurement values at multiple locations along the furnace radius at a predetermined position, accurate measurement of the layer thickness distribution of the charge in the batch in which the top layer charge is charged (Nth batch) is possible. It becomes possible. Therefore, it is not necessary to measure the surface profile of several tens batches and obtain the average of the measured values as in the past, and the batch in which the top layer charge is charged in a short time (time required for two consecutive batches). The layer thickness distribution of the (Nth batch) can be measured.

  In addition, in order to define the estimated shape line 50 that simulates the surface profile, by using the first continuous function F1 in which the portion near the furnace wall is a straight line (first linear function f11), the radial direction of the blast furnace It is possible to derive the estimated shape line 50 in accordance with the actual surface profile even in the operation state in which the charge is charged to the furnace wall side while changing the drop position, and this makes it possible to accurately measure the layer thickness distribution. it can.

  Specifically, the region of the surface profile that is horizontal or substantially horizontal in the vicinity of the furnace wall by charging the furnace wall while changing the dropping position in the furnace radial direction (X-axis direction) A so-called terrace portion is formed (see FIG. 2). This is because when the charge is charged so that it falls to a certain position in the furnace radial direction, the shape of the layer top surface of the charge along the furnace radial direction rises like a mountain with this position at the top. If the charge is charged while changing the dropping position along the furnace radial direction, the charge is distributed almost uniformly in the furnace radial direction. This is because becomes substantially flat. Since the top surface shape of the terrace portion is a horizontal or substantially horizontal straight line shape at the furnace radius at a predetermined position, the first primary function F1 in the section closest to the furnace wall in the first continuous function F1 in which a plurality of functions are connected. By using the function f11, the shape of the terrace portion is imitated with high accuracy. As described above, the portion 52 in the vicinity of the furnace wall in the estimated shape line 50 simulates the surface profile with high accuracy, so that the layer thickness distribution can be measured with high accuracy.

That is, since the estimated shape line 50 is rotated around the furnace center axis to obtain the estimated surface profile, the moving distance of the site where the error occurs when the estimated shape line is rotated toward the furnace wall side becomes longer, As a result, errors on the furnace wall side have a greater effect on the estimated surface profile. Therefore, in the position correction unit 24, the estimated shape lines 50 _n-1 and 50 _{n in} which errors occur on the furnace wall side as described above are used, and the volume of the space sandwiched between the two upper and lower estimated surface profiles is the first. It is corrected so that the estimated shape line 50 _n-1 of the (N-1) th batch is positioned below the estimated shape line 50 _n of the Nth batch so as to be equal to the volume of the entire charge of the Nth batch. Then, the influence of the error appears greatly, and the measurement accuracy of the layer thickness distribution decreases. Accordingly, in the derivation of the estimated shape line 50, the portion in the vicinity of the furnace wall conforms to the shape of the terrace portion of the actual surface profile, that is, has a shape with less error, so that two consecutive batches (N-1th time) It is possible to accurately derive the shape of each estimated surface profile of the batch and the Nth batch) and the distance between these two estimated surface profiles.

  Further, the first continuous function F1 can simulate the surface profile with high accuracy not only in the terrace portion of the surface profile but also in the portion closer to the furnace center than the terrace portion at the furnace radius at a predetermined position. Specifically, a second linear function f12, which is a straight line of a downward gradient, is a portion of the surface profile at a predetermined position of the furnace radius that is located closer to the furnace center than the terrace portion and has a constant downward gradient toward the furnace center. The second curve function f22, which is a curve that bulges upward in the center of the furnace and that bulges up at the center, can be simulated with high accuracy. Therefore, by using the first continuous function F1, it is possible to obtain an estimated shape line 50 that conforms to the actual surface profile in the entire range from the furnace wall to the furnace center at the furnace radius at a predetermined position.

  Furthermore, since the estimated shape line 50 is used to derive the estimated surface profile, the layer thickness distribution can be easily and accurately calculated based on these measured values even if noise is included in the measured values obtained by measuring the surface profile. Can be derived.

  Specifically, an estimated shape line 50 at a predetermined furnace radius is defined, and a first continuous function F1 in which a plurality of functions f11, f12, f21, and f22 is connected is set in advance, and each function f11 is based on a plurality of measured values. , F12, f21, f22, and the range of the furnace radius at a predetermined position, the estimated shape line 50 as a line having a curve is derived regardless of the magnitude and amount of noise included in each measurement value. Is done. The surface profile at the furnace radius at a predetermined position is simulated by the estimated shape line 50, so that even if there is a partial protrusion or the like (noise) that occurs in the measured layer upper surface shape, A highly accurate surface profile can be obtained without being affected by protrusions or the like.

  Furthermore, in the latest batch (Nth batch) and the previous batch (N-1th batch), the surface profile at the furnace radius at a predetermined position was actually measured, and the measured value was measured. By calculating the coefficient of each function f11, f12, f21, f22 defining the first continuous function F1 set in advance based on (actually measured value) and the range of the furnace radius at a predetermined position, the surface profile is obtained. The estimated shape line 50 with good accuracy is derived, whereby the measurement of the layer thickness distribution corresponding to the change in the layer thickness distribution under the same operating condition of the blast furnace is performed with high accuracy.

  In the blast furnace in operation, the surface profile of two continuous batches (C2 batch (N-1st batch) and O1 batch (Nth batch)) is measured by the measuring means 12, and the measurement obtained by this measurement When the layer thickness distribution is estimated by the layer thickness distribution measuring method (estimated shape line defined by the first continuous function F1) according to the first embodiment based on the values, and the layer thickness distribution measuring method Comparison is made with the case where the layer thickness distribution is measured without using (by the method using the unprocessed measurement value), and the results are shown in FIGS. In this example, the method using the layer thickness distribution measuring method according to the first embodiment is “present method”, and the method not using is “unprocessed”.

  FIG. 6A shows a surface profile obtained by performing C2 batch and O1 batch one by one, and plotting the measurement values obtained for each batch unprocessed on XY coordinates. FIG. 6B shows a surface profile obtained by plotting on the XY coordinates values derived using the present method from the same measurement values used in FIG. FIG. 6C shows a surface profile obtained by performing 30 batches each of the C2 batch and the O1 batch, and plotting an average of the measurement values obtained for each batch on the XY coordinates. The average of 30 batches is a batch at a time when the operating conditions such as the charged amount of the charge and the armor stroke are constant and the furnace condition hardly changes. Therefore, it is considered that the obtained surface profile is less affected by noise. FIG. 7 shows the layer thickness of the O1 batch (Nth batch) shown in FIGS. 6 (a) to 6 (c).

  As can be seen from FIG. 7, the C2 batch and the O1 bat are each subjected to one batch, and the layer thickness of the surface profile (FIG. 6 (a)) using the measurement values obtained for each batch unprocessed is 30 batches. The average thickness of the surface profile (FIG. 6C) is greatly deviated at a position near the furnace radius of 1.775 m. On the other hand, the surface profile (FIG. 6B) using this method approximates the layer thickness of the surface profile averaged over 30 batches (FIG. 6C). Moreover, FIG. 8 shows the error of the layer thickness when the measurement value obtained for each batch is processed by this method with respect to the layer thickness of the surface profile averaged over 30 batches, and the error of the layer thickness when not processed. It is the table | surface which showed. From this table, it can be seen that the error is smaller when this method is processed than when it is not processed.

  From the above, by using the above-described layer thickness distribution measurement method, an accurate layer thickness distribution that could not be obtained unless several tens of batches were measured can be accurately obtained by measuring the surface profile of two consecutive batches. It can be confirmed that it is obtained.

  Next, a second embodiment of the present invention will be described with reference to FIG. 9. The same reference numerals are used for the same components as those in the first embodiment, and a detailed description thereof will be omitted. Only different components will be described in detail. explain.

  The layer thickness distribution deriving unit 20 of the measuring apparatus 110 according to the present embodiment further includes a receiving unit 30 and a state storage unit 32.

  The receiving unit 30 is connected to the control unit of the blast furnace and the like, to the furnace wall side while changing the fall position in the radial direction of the furnace when charging a new charge into the blast furnace transmitted from the control unit or the like. This is a part that receives, as an output signal, information as to whether or not charging has been performed (charging state information). The layer thickness distribution deriving means 20 is not limited to a configuration in which the receiving state information received from the outside is received by the receiving unit 30 as an output signal. May be configured to be input from the input unit. For example, the input unit includes a keyboard and a touch panel.

  The state storage unit 32 is a part that stores the charging state information received by the receiving unit 30. Further, when the input unit is provided in the layer thickness distribution deriving means 20, the state storage unit 32 can also store the charging state information input from the input unit.

  The function storage unit 121 stores two continuous functions F1 and F2 whose coefficients are undetermined. Specifically, in addition to the first continuous function F1 defined by a series of four functions f11, f12, f21, and f22 similar to the first embodiment, a series of three functions f13, f23, and f24. The prescribed second continuous function F2 is stored.

  As shown in FIG. 10, the second continuous function F2 is a function represented on the XY plane with the furnace center axis as the Y axis and the furnace radial direction at a predetermined position as the X axis. 51, an inclined portion 54, and a central curve portion 55. Here, the furnace wall curve portion 51 is a portion that is located in the vicinity of the furnace wall of the line defined by the second continuous function F2 and forms a curve that bulges upward, and the center curve portion 55 is The portion of the second continuous function F2 that is located in the center of the furnace and forms a curved portion that bulges upward, and the inclined portion 54 is the furnace of the line defined by the second continuous function F2. It is a part which is located between the wall curve part 51 and the center curve part 55, and becomes a straight line with a certain downward gradient toward the furnace center. And like 1st Embodiment, since the shape of the furnace wall curve part 51, the gradient and length of the inclination part 54, the shape of the center curve part 55, etc. differ for every surface profile, it is not decided, but by the measurement means 12 Each coefficient is determined by calculating the coefficient of the second continuous function F2 based on a plurality of measurement values obtained by measuring the surface profile, and an estimated shape line 150 approximating the actual surface profile can be obtained. .

  Specifically, the second continuous function F2 is defined by a plurality of functions that are continuous along the X-axis direction. The second continuous function F2 is defined by a third curve function f23, a third linear function f13, and a fourth curve function f24 that are arranged in order from the furnace wall toward the furnace center along the X axis. Yes. The third linear function f13 and the third curve function f23 are connected so as to be differentiable at the boundary position r4.

  The third curve function f23 is a function in which the line defined by the function f23 becomes a curve corresponding to the furnace wall curve portion 51, and the third linear function f13 is an inclination of the line defined by the function f13. The fourth curve function f24 is a function in which a line defined by the function f24 becomes a curve corresponding to the central curve portion 55.

If each function f13, f23, f24 is shown more concretely, the third linear function f13 is
y = e ₁₁ x + f ₁₁ (hereinafter, also simply referred to as “formula (110)”)
It is represented by The third curve function f23 is
y = (− 1 / a ₁ ) · log | cos (a ₁ x + b ₁ ) | + c ₁ (hereinafter, also simply referred to as “expression (111)”)
And the fourth curve function f24 is
y = α ₁ x ² + β ₁ x + γ ₁ (hereinafter, also simply referred to as “formula (113)”)
It is represented by Note that a ₁ , b ₁ , c ₁ , e ₁₁ , f ₁₁ , α ₁ , β ₁ , and γ ₁ are coefficients.

  Hereinafter, in the furnace radius (X axis) at the predetermined position, the boundary position between the furnace wall curved portion 51 and the inclined portion 54 is r4, and the boundary position between the second inclined portion and the central curved portion is r5.

  The estimated shape line deriving unit 123 is charged in the charging state information stored in the state storage unit 32 (information on whether or not the charging material has been charged to the furnace wall side while changing the dropping position in the furnace radial direction). Based on the continuous function (first continuous function) used to derive the estimated shape line 50 (or 150) so that the estimated shape line 50 (or 150) having a shape corresponding to the measured surface profile can be derived in a shorter time. F1 and the second continuous function F2) are switched, and the coefficient of one of the selected first continuous function F1 or second continuous function F2 is calculated. Then, the estimated shape line deriving unit 123 stores the continuous function F1 or F2 into which the calculated coefficient is substituted.

  The measuring apparatus 110 according to the present embodiment has the above-described configuration. Next, the operation of the measuring apparatus 110 will be described with reference to FIG.

  First, when a new charge is charged into the blast furnace in the (N-1) th batch, the charge state information (changes the fall position in the furnace radial direction) transmitted as an output signal from the control unit of this blast furnace. However, the receiving unit 30 of the layer thickness distribution deriving means 20 receives the information on whether or not the charging material has been charged to the furnace wall side, and the received charging state information is stored in the state storage unit 32. (Step 0).

  On the other hand, the surface profile at a predetermined furnace radius of the charge of the (N-1) th batch is measured by the measuring means 12, and the measured value obtained by this measurement is the layer thickness distribution deriving means 20 (the estimated shape line deriving section 123). ) By an output signal (step 1). In the present embodiment, when the charge is charged to the furnace wall side while changing the dropping position along the furnace radius, the top surfaces of 38 layers at equal intervals on the X axis as in the first embodiment. When the charge is charged to the furnace wall side without changing the dropping position of the furnace radius, 26 points are equally spaced on the X axis. The height of the upper surface of the layer is measured (see FIG. 12).

  Note that step 0 and step 1 may be performed in order or in parallel. Further, in step 1, the number of measurement points on the X-axis does not need to be changed depending on whether or not the charge is charged to the furnace wall side while changing the dropping position of the furnace radius, but is derived. A larger layer thickness distribution is preferable because the layer thickness distribution becomes more accurate.

  The estimated shape line deriving unit 123 that has received the output signal from the measuring means 12 acquires the charging state information of the (N-1) th batch stored in the state storage unit 32, and based on this information, the output signal By applying the first continuous function F1 (see FIG. 3) or the second continuous function F2 (see FIG. 10) stored in the function storage unit 121 to the obtained measurement value, the coefficient of the continuous function F1 or F2 Is derived (step 2a).

  Specifically, the estimated shape line deriving unit 123 inserts the charging state information of the (N-1) th batch stored in the state storage unit 32 into the furnace wall side while changing the dropping position along the furnace radial direction. In the case of information on the loading of an object, the coefficient of the first continuous function F1 is calculated by the steepest descent method as in the first embodiment, based on the measurement value obtained by the measurement of the measurement means 12. The first continuous function F1 to which this coefficient is substituted is stored as the estimated shape line 50 of the (N-1) th batch. On the other hand, the estimated shape line deriving unit 123 is information in which the charging state information stored in the state storage unit 32 is charged into the furnace wall without changing the dropping position in the furnace radial direction. In this case, the coefficient of the second continuous function F2 is calculated by the steepest descent method based on the measurement value obtained by the measuring means 12, and the second continuous function F2 to which this coefficient is substituted is the (N-1) th. Stored as the estimated shape line 150 of the second batch. In the present embodiment, the charged material is charged while changing the dropping position in the furnace radial direction in the (N-1) th batch (C2 batch), and the falling in the furnace radial direction in the Nth batch (O1 batch). The charge is charged without changing the position (see FIG. 12).

  Next, the receiving unit 30 of the layer thickness distribution deriving means 20 receives the charging state information transmitted as an output signal from the control unit of the blast furnace in the Nth batch, and the received charging state information is the state storage unit 32. (Step 3a). On the other hand, the surface profile of the Nth batch is measured by the measuring unit 12, and the measurement value obtained by this measurement is transmitted to the layer thickness distribution deriving unit 20 (estimated shape line deriving unit 123) (step 3).

  The estimated shape line deriving unit 123 that has received the output signal from the measuring unit 12 is based on the charging state information of the Nth batch stored in the state storage unit 32 as in the case of the N-1th batch. By applying the first continuous function F1 (see FIG. 3) or the second continuous function F2 (see FIG. 10) stored in the function storage unit 22a to the measurement value obtained by the above, the continuous function F1 or F2 The coefficient of is derived. Then, the estimated shape line deriving unit 123 assigns the derived coefficient to the first continuous function F1 or the second continuous function F2, and stores this as the estimated shape line 50 or 150 of the Nth batch (step) 4a).

As described above, when the estimated shape line 50 _n-1 of the (N _-1) th batch and the estimated shape line 150 _n of the Nth batch are derived and stored in the estimated shape line deriving unit 123, respectively, Similar to the second batch (step 2a in FIG. 11), the position correction unit 24 acquires the value of the volume of the entire charge of the Nth batch from the volume storage unit 22, and sends the value to the estimated shape line deriving unit 123. The interval between the two estimated shape lines 50 _n−1 and 150 _n is stored (step 5: refer to FIG. 12).

Next, the layer thickness distribution deriving unit 25 is based on the vertical distance at each position along the X axis in each of the estimated shape lines 50 _n−1 and 150 _n of the N−1th batch and the Nth batch. A layer thickness distribution along the X axis from the furnace center to the furnace wall in the Nth batch charge is derived, and the value of this layer thickness distribution is transmitted to the output means 14 by an output signal (step 6). The output means 14 that has received this output signal displays the layer thickness distribution in the Nth batch as layer thickness distribution information (step 7), and the measurement of the layer thickness distribution is completed. As in the first embodiment, when the layer thickness distribution of the (N + 1) th batch is measured continuously, the process starts from step 3a (step 8).

  The continuous function F1 or F2 used when the estimated shape line deriving unit 123 derives the estimated shape line 50 or 150 based on the charging state information acquired from the blast furnace or the like and stored in the state storage unit 32 as described above. The second continuous function F2 based on the measured value obtained by the measuring means 12 even in the batch in which the charge is charged to the furnace wall side without changing the dropping position in the furnace radial direction. It is possible to obtain an estimated shape line 150 having a shape conforming to the surface profile of the corresponding batch by deriving the coefficient of each function that defines the range and the range in the furnace radius at a predetermined position.

  Specifically, by charging without changing the fall position in the furnace radial direction, the charge is shaped like a mountain with the drop position at the top at the furnace radius (X axis) at a predetermined position. Therefore, it is difficult to form a terrace portion. Therefore, by using the second continuous function F2 which is a shape in which the portion in the vicinity of the furnace wall is raised like a mountain (curve bulging upward (third curve function f23)), It is possible to accurately simulate the shape of the portion in the vicinity of the furnace wall (the shape that rises like a mountain). That is, when the terrace portion is not formed, the first continuous function F1 in which the vicinity of the furnace wall is a straight line is not used, but the second continuous function F2 in which the vicinity of the furnace wall is a curve bulging upward is used. Thus, the shape of the surface profile can be simulated with high accuracy. Therefore, even in a batch in which the charge is charged to the furnace wall side without changing the fall position in the furnace radial direction, the surface profile in the vicinity of the furnace wall approximated a surface profile having a mountain-like area. The estimated shape line 150 can be derived, whereby the layer thickness distribution can be accurately measured.

  Further, the third linear function f1 is used rather than the continuous function (first continuous function) F1 defined by the first linear function f11, the first curve function f21, the second linear function f12, and the second curve function f22. The number of functions that define the continuous function is smaller when the continuous function (second continuous function) F2 defined by the curve function f23, the third linear function f13, and the fourth curve function f24 is used. Is simplified, and the estimated shape line 50 is derived using the first continuous function F1 regardless of whether or not the charge is charged to the furnace wall side while changing the fall position in the furnace radial direction. Compared to the case (measurement apparatus 100 of the first embodiment), it is possible to obtain the estimated shape line 150 having a shape conforming to the actual surface profile easily and in a short time.

  Note that the measuring devices 10 and 110 according to the present invention are not limited to the first and second embodiments described above, and various changes can be made without departing from the scope of the present invention.

  For example, in the present embodiment, the first continuous function F1 is defined by sequentially connecting the above equations (10) to (13) from the furnace wall to the furnace center side, but it is necessary to be limited to this. There is no. Specifically, in the first continuous function F1, each function defining the connection curve portion 53 and the center curve portion 55 is expressed in the x-axis direction at the angle formed by the tangent line of the curve defined by the function and the x-axis. Any function may be used as long as the angle change rate with respect to the position is constant. For example, in the first continuous function F1, the function that defines the connection curve portion 53 and the function that defines the center curve portion 55 are opposite, that is, the function that defines the connection curve portion 53 is the equation (13). The function that defines 55 may be equation (11), and the function that defines the connection curve portion 53 and the function that defines the center curve portion 55 are the same, that is, the function that defines the connection curve portion and the center curve. Both the function defining the part may be the expression (11) or the expression (13). Also in the second continuous function F2, as in the first continuous function F1, each function that defines the furnace wall curve portion 51 and the central curve portion 55 is the tangent to the curve defined by the function and x Any function may be used as long as the angle change rate with respect to the position in the x-axis direction at the angle formed by the axis is constant.

  Further, the number of functions defining the continuous function is not limited to three or four, and may be five or more.

  In addition, when a blast furnace operation in which the central charging is not performed is performed in the charging of the charged material into the blast furnace, the estimated shape line without the central curve portion 55 (the second curved function or the fourth curved function). Can be applied to the same method as described above to obtain the layer thickness distribution. That is, when the center charging is not performed, the continuous function defining the estimated shape line is the first linear function f11 and the first curve function that are arranged in order from the furnace wall toward the furnace center along the X axis. It is defined by f21 and the second linear function f12, or is defined by the third curve function f23 and the third linear function f13.

It is a block diagram of the blast furnace charge layer thickness distribution measuring apparatus concerning a 1st embodiment. It is the schematic of the measurement means of the layer thickness distribution measuring apparatus of the blast furnace charge which concerns on the same embodiment. It is a figure which shows the 1st continuous function in the layer thickness distribution measuring apparatus of the blast furnace charge which concerns on the same embodiment, and the several function which prescribes | regulates this 1st continuous function. It is a figure which shows the flow of operation | movement in the layer thickness distribution measuring apparatus of the blast furnace charge which concerns on the same embodiment. It is a figure which shows the correction | amendment in the position correction part of the layer thickness distribution measuring apparatus of the blast furnace charge which concerns on the same embodiment. In the measurement of the layer thickness distribution in the layer thickness distribution measuring device for the blast furnace charge according to the embodiment, (a) shows the surface profile of the actual measurement value, and (b) shows the layer thickness distribution estimation method according to the embodiment. (C) is a figure which shows the average surface profile of 30 batches. It is a figure which shows the layer thickness based on each profile in Fig.6 (a) thru | or FIG.6 (c) in the measurement of the layer thickness distribution in the layer thickness distribution measuring apparatus of the blast furnace charge which concerns on the same embodiment. It is a figure which shows the error of the layer thickness in FIG. It is a block diagram of the layer thickness distribution measuring apparatus of the blast furnace charge which concerns on 2nd Embodiment. It is a figure which shows the 2nd continuous function and the several function which prescribes | regulates this 1st continuous function in the layer thickness distribution measuring apparatus of the blast furnace charge which concerns on the same embodiment. It is a figure which shows the flow of operation | movement in the layer thickness distribution measuring apparatus of the blast furnace charge which concerns on the same embodiment. It is a figure which shows the correction | amendment in the position correction part of the layer thickness distribution measuring apparatus of the blast furnace charge which concerns on the same embodiment.

Explanation of symbols

10 Measuring device 12 Measuring means (microwave profile meter)
14 Output unit 20 Layer thickness distribution deriving unit 21 Function storage unit 22 Volume storage unit 23 Estimated shape line deriving unit 24 Position correcting unit 25 Layer thickness distribution deriving unit 50 Estimated shape line 50 _n-1 N-1th batch estimated shape Line 50 _n Estimated shape line of Nth batch 51 Furnace wall curve portion 52 Straight line portion 53 Connection curve portion 54 Inclination portion 55 Central curve portion F1 First continuous function f11 First linear function f12 Second linear function f21 1 curve function f22 2nd curve function

Claims (8)

Of the charges stacked in the blast furnace by repeating the batch for charging the charges into the blast furnace n times (where n is a natural number arbitrarily selected within the range of n> 1), the nth time A method of measuring a layer thickness distribution along a furnace radial direction from a furnace center to a furnace wall in a layer formed by a batch charge,
A measurement step of measuring the top surface shape of the charge in the furnace radius at a predetermined position in each batch of the (n-1) th batch and the nth batch;
An estimated shape line that simulates the upper surface shape of the layer as a line having a curve is derived for each batch of the (n−1) -th batch and the n-th batch based on the measurement value of each batch obtained in the measurement step. A shape line deriving step,
When the estimated surface profiles are obtained by rotating the two estimated shape lines derived in this shape line deriving step so that the furnace center axis is the center of rotation, the volume of the space sandwiched between both estimated surface profiles is the first. The two pieces are arranged such that the estimated shape line of the n-1st batch is positioned below the estimated shape line of the nth batch at an interval equal to the volume of the entire charge of the nth batch. A correction step for correcting at least one of the estimated shape lines;
The furnace in the charge of the nth batch based on the vertical spacing at each position along the furnace radial direction in each estimated shape line of the n-1st batch and the nth batch corrected in this way A layer thickness distribution deriving step for deriving a layer thickness distribution along the furnace radial direction from the center to the furnace wall;
In the shape line derivation step, a first continuous function that defines the estimated shape line is set in advance, and the first continuous function is a plurality of functions continuous in the furnace radial direction of the predetermined position, A first linear function in which the coefficient of each function and the range of the predetermined position in the furnace radial direction are undetermined and adjacent functions are arranged in order from the furnace wall to the furnace center side so that adjacent functions are continuous at the boundary, Defined by a curve function, a second linear function, and a second curve function;
The first linear function is a function in which a line defined by the function is positioned in the vicinity of the furnace wall and becomes a horizontal or substantially horizontal straight line, and the second linear function is defined by the function. The line is a function that is positioned closer to the furnace center than the first linear function and becomes a straight line having a constant downward slope toward the furnace center, and the first curve function is a line defined by the function. It is located between the first linear function and the second linear function, and is defined by a substantially horizontal straight furnace center side end defined by the first linear function and the second linear function. The second curve function is a function that smoothly connects a straight line with a straight line with a downward slope, and a line defined by the function is located at the center of the furnace and expands upward. It is a function that becomes a curve to be issued,
In each batch, by calculating the coefficient in each function and the range in the furnace radial direction of the predetermined position for each batch based on the measurement value of the batch obtained in the measurement step, the n-1st batch and A method for measuring a layer thickness distribution of a blast furnace charge, wherein an estimated shape line for each batch with the n-th batch is derived. In the method for measuring the layer thickness distribution of the blast furnace charge according to claim 1,
In the shape line derivation step, a second continuous function having an undetermined coefficient that defines the estimated shape line is further set, and the second continuous function is a plurality of functions that are continuous in the furnace radial direction at the predetermined position. The coefficient of each function and the range in the furnace radial direction of the predetermined position are undetermined, and adjacent functions continue from the furnace wall toward the furnace center along the furnace radius at the predetermined position so that the adjacent functions are continuous at the boundary. Defined by a third curve function, a third linear function, and a fourth curve function,
The third curve function is a function in which a line defined by the function is located in the vicinity of the furnace wall and bulges upward, and the fourth curve function is a line defined by the function. Is a function that is located in the center of the furnace and forms a curve that bulges upward, and the third linear function is such that the line defined by the function is the third curve function and the fourth curve function. Is a function that is a straight line with a constant downward slope toward the furnace center,
In each batch of the (n-1) th batch and the nth batch, the charge is placed in the blast furnace while changing the drop position along the furnace radial direction when the charge is charged to the furnace wall side. In the case of charging, the coefficient in each function defining the first continuous function and the range in the furnace radial direction of the predetermined position are calculated based on the measurement value of the batch obtained in the measurement step. Deriving the estimated shape line of the batch,
When charging the charging material into the blast furnace without changing the dropping position in the furnace radial direction when charging the charging material to the furnace wall side, the measured value of the batch obtained in the measuring step The estimated shape line of the batch is derived by calculating a coefficient in each function that defines the second continuous function and a range in the furnace radial direction of the predetermined position based on Layer thickness distribution measurement method. In the method for measuring the layer thickness distribution of the blast furnace charge according to claim 1 or 2,
The first continuous function or the second continuous function is a function represented on an xy plane having a furnace center axis as a y-axis and a furnace radius at the predetermined position as an x-axis.
The first curve function, the second curve function, the third curve function, or the fourth curve function is an x-axis direction at an angle between a tangent line of the curve defined by the curve function and the x-axis. A method for measuring a layer thickness distribution of a blast furnace charge, wherein an angle change rate with respect to a position of the blast furnace is a constant function. In the method for measuring the layer thickness distribution of the blast furnace charge according to claim 3,
The first curve function, the second curve function, the third curve function, or the fourth curve function is a function represented by the following expression (1) or (2): A method for measuring the layer thickness distribution of the blast furnace charge.
y = (− 1 / a) · log | cos (ax + b) | + c (1)
y = αx ² + βx + γ (2)
Here, a, b, c, α, β, and γ are coefficients. In the method for measuring the layer thickness distribution of the blast furnace charge according to claim 4,
In the shape line derivation step, coefficients in all the functions that define the first continuous function or the second continuous function using the steepest descent method from the measurement values obtained in the measurement step, and ranges in the x-axis direction The layer thickness distribution measuring method of the blast furnace charge characterized by calculating | requiring simultaneously. Of the charges stacked in the blast furnace by repeating the batch for charging the charges into the blast furnace n times (where n is a natural number arbitrarily selected within the range of n> 1), the nth time An apparatus for measuring a layer thickness distribution along a furnace radial direction from a furnace center to a furnace wall in a layer formed by a batch charge,
In each of the n-th batch and the n-th batch, in the charge of the n-th batch based on the measurement value obtained by measuring the top surface shape of the charge in the furnace radius at a predetermined position. A layer thickness distribution deriving means for deriving a layer thickness distribution along the furnace radial direction from the furnace center to the furnace wall;
Output means for outputting the layer thickness distribution derived by the layer thickness distribution deriving means to the outside,
The layer thickness distribution deriving means defines an estimated shape line that simulates the shape of the upper surface of the layer as a line having a curve, and a plurality of functions are coefficients in each function and a range in the furnace radial direction of the predetermined position are not determined and adjacent to each other A function storage unit that stores in advance a first continuous function defined by being continuous along the furnace radial direction of the predetermined position so that the functions to be performed are continuous at the boundary;
A volume storage unit for storing a volume value of the entire charge charged into the blast furnace in the n-th batch;
A coefficient in each function that defines the first continuous function stored in the function storage unit for each batch based on the measured values of each batch of the (n-1) th batch and the nth batch; By calculating a range in the furnace radial direction of the predetermined position, an estimated shape line for each batch of the (n−1) -th batch and the n-th batch is derived, and two of these derived estimated shape lines are derived. An estimated shape line deriving unit for storing each value;
When two estimated shape lines whose values are stored in this estimated shape line deriving section are rotated so that the furnace center axis is the center of rotation, and the estimated surface profiles are obtained, they are sandwiched between both estimated surface profiles. The n-1st batch below the estimated shape line of the nth batch at an interval such that the volume of the space is equal to the volume of the entire charge of the nth batch stored in the volume storage unit. A position correction unit that corrects at least one of the two estimated shape lines so that the estimated shape line is positioned;
In the charge of the nth batch based on the vertical spacing at each position along the furnace radial direction in each estimated shape line of the n-1st batch and the nth batch corrected in this way. A layer thickness distribution deriving unit for deriving a layer thickness distribution along the furnace radial direction from the furnace center to the furnace wall, and transmitting a value of the derived layer thickness distribution to the output means;
The first continuous function stored in advance in the function storage unit is a first linear function, a first curve function, a first curve function arranged in order from the furnace wall toward the furnace center along the furnace radius at the predetermined position. Defined by a linear function of 2 and a second curve function,
The first linear function is a function in which a line defined by the function is positioned in the vicinity of the furnace wall and becomes a horizontal or substantially horizontal straight line, and the second linear function is defined by the function. The line is a function that is positioned closer to the furnace center than the first linear function and becomes a straight line having a constant downward slope toward the furnace center, and the first curve function is a line defined by the function. Located between the first linear function and the second linear function, a horizontal or substantially horizontal straight furnace center side end defined by the first linear function and the second linear relation The second curve function is a function in which the line defined by the function is located in the furnace center, Of the blast furnace charge, which is a function that forms a curved line Thickness distribution measurement device. In the blast furnace charge layer thickness distribution measuring apparatus according to claim 6,
In the n-th batch, an input unit for inputting charging state information indicating whether charging has been performed on the furnace wall side while changing the radial drop position or the charging state information is transmitted from the outside. At least one of the receiving units,
A state storage unit for storing the input or transmitted charging state information,
The function storage unit defines the estimated shape line, and the predetermined function is such that a plurality of functions have a coefficient in each function and a range in the furnace radial direction of the predetermined position and adjacent functions are continuous at the boundary. A second continuous function that continues along the furnace radial direction of the position is further stored,
The second continuous function stored in the function storage unit includes a third curved function, a third linear function, and a third linear function arranged in order from the furnace wall toward the furnace center along the furnace radius at the predetermined position. 4 is defined by the curve function,
The third curve function is a function in which a line defined by the function is located in the vicinity of the furnace wall and bulges upward, and the fourth curve function is a line defined by the function. Is a function that is located in the center of the furnace and forms a curve that bulges upward, and the third linear function is such that the line defined by the function is the third curve function and the fourth curve function. Is a function that is a straight line with a constant downward slope toward the furnace center,
In the estimated shape line deriving unit, the furnace wall of the charged material in each batch of the (n-1) th batch and the nth batch based on the charging state information stored in the state storage unit In the case of charging the charged material into the blast furnace while changing the drop position along the furnace radial direction at the time of charging to the side, it is stored in the function storage unit based on the measured value of the batch. Calculating a coefficient in each function that defines the first continuous function and a range in the furnace radial direction of the predetermined position;
When charging the charging material into the blast furnace without changing the dropping position in the furnace radial direction when charging the charging material to the furnace wall side, the function is based on the measured value of the batch. By calculating a coefficient in each function defining the second continuous function stored in the storage unit and a range in the furnace radial direction of the predetermined position, the batch of the n-1st batch and the nth batch An apparatus for measuring a layer thickness distribution of a blast furnace charge, wherein each estimated shape line is derived and values of the two derived estimated shape lines are stored. In the blast furnace charge layer thickness distribution measuring apparatus according to claim 6 or 7,
Blast furnace charging characterized by further comprising measuring means for measuring the shape of the top surface of the charge at the furnace radius at the predetermined position and transmitting the measured value obtained by this measurement to the estimated shape line deriving means Equipment thickness distribution measuring device.

Images