JP2009097048A - Method for estimating distribution of layer thickness of charged material in blast furnace and instrument using the method - Google Patents

Abstract

<P>PROBLEM TO BE SOLVED: To provide a method for estimating the distribution of a layer thickness of charged materials in a blast furnace, with which the distribution of the layer thickness of the charged materials along the furnace radius direction from the center of the furnace to a furnace side wall can be obtained with high precision in a short period of time, and an instrument for measuring the distribution of the layer thickness of the charged materials in the blast furnace, using the method. <P>SOLUTION: A continuous function (f) having non-fixed factors is set, which regulates an estimated shape line pL for simulating the layer upper surface shape of the charged material in the furnace radius at the prescribed position, the layer upper surface is measured with a measuring means in the (n-1)th batch and the (n)th batch, respectively, and the respective estimated shape line at the (n-1)th batch and the (n)th batch are derived by calculating the factor in the continuous function (f) on the basis of the measured value obtained by the above measurement. <P>COPYRIGHT: (C)2009,JPO&INPIT

Description

  The present invention relates to a method for estimating a layer thickness distribution along a furnace radial direction from a furnace center to a furnace side wall of charges such as ore and coke charged in a blast furnace, and an apparatus using this method.

  In blast furnace operation, ore (iron ore, etc.), which is an iron source, and coke are alternately charged into the furnace from the top of the furnace, and air or pulverized coal is supplied from the tuyere and the coke, pulverized coal, etc. Burns. A cohesive zone is formed in the furnace as the ore is heated and reduced by the high temperature reducing gas generated by this combustion and rising from below in the furnace. A dripping zone in which the ore is further temperature-reduced and reduced is formed at the lower part of the cohesive zone, and becomes molten iron at the hearth and is discharged from the tap.

  In order to stably perform such blast furnace operation, it is very important to control the shape and formation position of the cohesive zone. Therefore, the layer thickness distribution along the furnace radial direction from the furnace center to the furnace side wall of the ore layer and coke layer stacked in the furnace is adjusted so that the shape and formation position of the cohesive zone are appropriate. Is called. The ore layer and the coke layer are formed by alternately charging ore and coke into the furnace from the top of the blast furnace furnace, and are alternately laminated.

  As a method for adjusting the layer thickness distribution, the mover armor is operated in the bell-arm type blast furnace, and the swivel chute is operated in the bell-less type blast furnace to change the charging position of the charges such as ore and coke into the furnace. It is adjusted by doing. In this adjustment method, the layer thickness distribution of the charged material is grasped, and based on the grasped layer thickness distribution, the next arm thickness distribution becomes a more suitable layer thickness distribution. The input device is operated.

  In order to grasp the layer thickness distribution of the charge, the shape of the top layer of the charge in the corresponding batch (unit of charging of the charge into the furnace), that is, the profile and the previous batch And the total charge of the batch in question. Specifically, the profile of the previous batch is set so that the volume of the portion sandwiched between the profile of the corresponding batch and the 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 profiles down. The layer thickness distribution is estimated by imitating the upper and lower width between the two upper and lower profiles set in this way as the layer thickness.

  Thus, in order to estimate the layer thickness distribution, it is necessary to obtain a charge profile. In order to obtain this profile, conventionally, a method of directly measuring the shape of the upper surface of the charge layer using a profile meter is known (see Patent Document 1).

  Specifically, the height of the upper surface of the charge layer at the furnace radius at a predetermined position is measured by a profile meter that uses microwaves or the like. When the measurement of the top surface of the charge in a given batch is completed, each measured value is plotted on the Y axis and the furnace radius is plotted on the X axis. By connecting the plotted points, the top layer shape of the charge along the radius at the predetermined furnace radius position can be obtained. By rotating this shape once so that the center axis of the furnace becomes the center of rotation, the shape of the layer upper surface of the charge of the batch, that is, the profile is obtained.

As another method, the relationship between the charging conditions and the profile is determined from past performance data, and various charging conditions (such as armature stroke, tilt angle of the chute, charging weight of charging material, and souding level) are determined. A method using a simulation model for estimating a profile as a calculation condition is known (see Patent Document 2).
Japanese Patent Laid-Open No. 04-24404 JP 2001-323306 A

  However, when the profile is obtained using the profile meter, an error occurs between the profile obtained from the measurement value and the actual profile. That is, there are many irregularities such as partial protrusions on the upper surface (surface) of the charge formed due to various causes at the time of charging only on the measured furnace radius. Therefore, even if the layer 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 profile may be different from the actual profile. Many.

  As described above, in the case of estimating the lowering amount of the lower profile so that the volume sandwiched between the two upper and lower profiles is the same as the volume of the charge, the error generated on the furnace wall side is particularly large. If it 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 to measure profiles of several tens of batches and calculate a profile 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 profile after calculating the average value of the measured values for every several tens of batches measured in the past, so that the latest layer thickness must be calculated. When changes occur, these changes cannot be grasped quickly. Therefore, there arises a problem that the operating conditions cannot be changed quickly 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 profile is estimated using the simulation model, only the same result can be calculated under the same operating conditions. However, the profile of the charge formed in the actual blast furnace is not always the same profile in spite of the same operating conditions, and is often different for each batch.

  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 has arisen.

  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 side wall of the charge. It is an object of the present invention to provide an estimation method and an apparatus for measuring a layer thickness distribution of a blast furnace charge using this method.

  Therefore, in order to solve the above-described problem, the method for estimating the layer thickness of the blast furnace charge according to the present invention performs N batches of charging the charge into the blast furnace (where N is a natural number and N> 1). Of the charges stacked in the blast furnace by repeating, the furnace center in the layer formed by the n-th charge (where n is a natural number arbitrarily selected within the range of 1 <n ≦ N) batch Is a method for estimating the layer thickness distribution along the furnace radial direction from the furnace side wall to the furnace side wall, and the coefficient undetermined continuous that defines the estimated shape line to simulate the shape of the top surface of the charge at the furnace radius at a predetermined position A function is set, and in each of the (n−1) -th batch and the n-th batch, the upper surface shape of the layer is measured by the measuring means, and the coefficient of the continuous function is calculated based on the measured value obtained by the measurement. By calculating, the n-1st batch and the nth time When the estimated profile lines are derived and the estimated profiles are obtained by rotating the derived two estimated profile lines so that the furnace center axis is the center of rotation, they are sandwiched between the two estimated profiles. The estimated shape line of the n-1st batch is positioned below the estimated shape line of the nth batch at an interval such that the volume of the space is equal to the charge volume of the nth batch. At least one of the two estimated shape lines is corrected, and the vertical distance at each position along the furnace radial direction in each of the estimated shape lines of the (n-1) th batch and the nth batch arranged in this way. And the layer thickness distribution along the furnace radial direction from the furnace center to the furnace side wall in the charge of the nth batch is estimated based on the calculation result.

  According to such a configuration, only by measuring the top surface shape of two consecutive batches (the (n-1) th batch and the nth batch), the furnace side wall from the furnace center in the desired charge of the nth batch is measured. The layer thickness distribution along the furnace radial direction is estimated. Therefore, there is no need to measure the top surface shape of several tens of batches and obtain the average of the measured values as in the prior art, and the layer thickness in the charge of the desired batch in a short time (the time required for two consecutive batches). Distribution is estimated.

  In addition, since a continuous function with an unknown coefficient that prescribes an estimated shape line to simulate the shape of the top surface of the charge at the furnace radius at a predetermined position is set in advance, A smooth estimated shape line is derived without being affected by a partial protrusion or the like that happens to occur when entering, and a more accurate profile can be obtained. Therefore, the interval between the two profiles is determined so that the volume of the space sandwiched between the two profiles of the (n-1) th batch and the nth batch is the same as the volume of the charge of the nth batch. When at least one of the two estimated shape lines is corrected, the vertical distance at each position along the furnace radial direction of both estimated shape lines is closer to the actual charge layer, so the layer thickness distribution is more accurate. Is obtained.

  Furthermore, in the batch where the layer thickness distribution is to be obtained and the previous batch, the layer top surface shape is measured, and a coefficient of a preset function is calculated based on the measured value. The estimated shape line with high accuracy is derived. As a result, the layer thickness distribution corresponding to the change of the layer thickness distribution under the same operating condition of the blast furnace is accurately estimated.

  In the blast furnace charge layer thickness estimation method according to the present invention, the continuous function includes a plurality of types of section functions defined for a plurality of sections sectioned along the furnace radius at the predetermined position. The plurality of interval functions are set such that the coefficient and the range in the furnace radial direction of the predetermined position are undecided and are continuous at the boundary between adjacent interval functions, and the n-1th batch In each of the n-th batch, the coefficients in all the interval functions and the ranges in the furnace radial direction of the predetermined positions are calculated based on the measured values, so that the n-1 th batch and the n-th batch are calculated. The configuration may be such that each estimated shape line of the batch is derived.

  According to this configuration, the continuous function is divided into a plurality of sections, and a section function that defines a line segment conforming to the shape of the upper layer surface is set for each of the plurality of sections. Therefore, it is possible to obtain an estimated shape line that more closely approximates the actual charge profile. As a result, the layer thickness distribution in the desired batch charge can be estimated more accurately.

  The N may be N = n and N> 2.

  According to such a configuration, it is only necessary to measure the layer upper surface shape in the uppermost layer charge in the blast furnace, that is, the layer upper surface shape exposed at the upper portion, so that measurement is facilitated. In addition, since the layer thickness distribution in the latest charge is obtained, it becomes easier to set operation conditions such as the charge position of the charge in the next batch.

  The interval function may be a function represented on an xy plane in which the furnace center axis is the y axis and the furnace radius at the predetermined position is the x axis.

  According to such a configuration, the plurality of interval functions are each expressed in the form y = f (x), and the handling of the functions becomes easy. Therefore, it is easy to calculate the coefficient in each interval function based on the measurement value of the layer upper surface shape and the range in the furnace radial direction of the predetermined position, and it is easy to derive the estimated shape line.

  Further, the continuous function is defined by three interval functions, and these three interval functions are first intervals in which the line segment defined by the corresponding interval function is a curve in order from the furnace center toward the furnace side wall. It may be a function, a linear function in which the line segment is a straight line, or a second interval function in which the line segment is a curve.

  According to such a configuration, an estimated shape line having a shape corresponding to the shape of the top surface of the charge can be obtained by connecting three line segments respectively defined by the three section functions in series. Since the estimated shape line can be defined with such a small number of interval functions, it is easy to calculate the coefficient in each interval function based on the measured value of the layer upper surface shape and the range in the furnace radial direction of the predetermined position, and derive the estimated shape line. Becomes easier. Moreover, since the derived estimated shape line is a shape that matches the shape of the measured layer upper surface, the layer thickness distribution of the charge can be estimated with high accuracy.

  The first interval function and the second interval function are expressed by a function having a constant angle change rate with respect to the x axis expressed by the following equation (1), or by the following equation (2). It may be a quadratic function.

y = −1 / a · log | cos (αx + b) | + c (1)
y = αx ² + βx + γ (2)
Here, a, b, c, α, β, and γ are coefficients.

  According to such a configuration, the first interval function and the second interval function are simple functions and have a small number of coefficients. Therefore, the coefficients in the interval functions based on the measurement values and the furnace radial direction of the predetermined position The range can be calculated more easily.

  Further, in each of the (n-1) th batch and the nth batch, the steepest descent method may be used from the measured values to simultaneously obtain the coefficients in all the interval functions and the ranges in the x-axis direction.

  According to such a configuration, the coefficients of all interval functions and the ranges in the x-axis direction in each batch can be easily obtained from the measured values at the same time.

  Further, in order to solve the above problems, the apparatus for estimating the layer thickness of the blast furnace charge according to the present invention measures the top surface shape of the charge in the blast furnace for each batch in which the charge is charged into the blast furnace. Possible measurement means, and among the charges stacked in the blast furnace by repeating the batch N times (where N is a natural number and 1 <N), the measurement means makes the n-1st batch The charging of the nth batch (where n is a natural number arbitrarily selected within the range of 1 <n ≦ N) based on the measured values obtained by measuring the shape of the upper surface of the two consecutive batches with the nth batch The layer thickness distribution estimating means for estimating the layer thickness distribution along the furnace radial direction from the furnace center to the furnace side wall in the layer formed by the object, and the value of the layer thickness distribution estimated by the layer thickness measuring means is output to the outside. An apparatus for estimating a layer thickness distribution of a blast furnace charge comprising an output means, wherein the layer thickness distribution is estimated The stage sets a continuous function with an undetermined coefficient that prescribes an estimated shape line for imitating the shape of the layer top surface of the charge at a furnace radius at a predetermined position, and the n-1th batch and the nth batch In each of the above, the upper surface shape of the layer is measured by the measuring means, and the coefficient of the continuous function is calculated based on the measurement value obtained by the measurement, so that the n−1th batch and the nth batch When each estimated shape line is derived and the estimated profiles are obtained by rotating these two estimated shape lines so that the furnace center axis is the center of rotation, the volume of the space sandwiched between both estimated profiles Is the same as the charge volume of the nth batch, and the two of the two pieces are positioned so that the estimated shape line of the (n-1) th batch is positioned below the estimated shape line of the nth batch. There are few estimated shape lines Both are corrected, and the vertical spacing at each position along the radial direction of the furnace is calculated for each estimated shape line of the n-1st batch and the nth batch arranged in this way. The value is transmitted to the output means as a value of the layer thickness distribution along the furnace radial direction from the furnace center to the furnace side wall in the charge of the nth batch, and the output means is transmitted from the layer thickness measuring means. The value of the layer thickness distribution is output to the outside.

  By adopting such a configuration, similarly to the above, it is possible to measure the top surface shape of two consecutive batches (the (n-1) th batch and the nth batch) from the center of the furnace in the charge of the nth batch. The layer thickness distribution along the furnace radial direction to the side wall is estimated, and the layer thickness distribution along the furnace radial direction from the furnace center to the furnace side wall in a predetermined batch of charge is estimated in a short time.

  In addition, since a continuous function with an indefinite coefficient that prescribes an estimated shape line for imitating the shape of the top surface of the charge at the furnace radius at a predetermined position is set in advance, the charge of the charge is added to the measured top surface shape of the layer. A smooth estimated shape line is derived without being affected by a partial protrusion or the like that happens to occur when entering, and a more accurate profile can be obtained. Therefore, the interval between the two profiles is determined so that the volume of the space sandwiched between the two profiles of the (n-1) th batch and the nth batch is the same as the volume of the charge of the nth batch. When at least one of the two estimated shape lines is corrected, the vertical distance at each position along the furnace radial direction of both estimated shape lines is closer to the actual charge layer, so the layer thickness distribution is more accurate. Is obtained.

  In addition, in the batch where the layer thickness distribution is to be obtained and the previous batch, the layer upper surface shape is measured, and a coefficient of a preset function is calculated based on the measured value. The estimated shape line with high accuracy is derived. As a result, the layer thickness distribution corresponding to the change in the layer thickness distribution under the same operating condition of the blast furnace is accurately estimated.

  Further, the value of the obtained layer thickness distribution is output by the output means, so that an operator operating the furnace can accurately grasp the layer thickness distribution of the desired batch charge.

  The N may be N = n and N> 2, and the measuring unit may be a profile meter capable of measuring the layer upper surface shape of the topmost charge.

  According to such a configuration, it is only necessary to measure the layer upper surface shape in the uppermost layer charge in the blast furnace, so that measurement is facilitated. In addition, since various profile meters have already been developed, it is possible to obtain a reliable measurement means at a low cost.

  As described above, according to the present invention, the method for estimating the layer thickness 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 side wall of the charge, and An apparatus for measuring the layer thickness distribution of the blast furnace charge using this method can be provided.

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

  First, based on FIG. 1, the blast furnace 10 used in this embodiment and the charging method of the charge to the said blast furnace 10 are demonstrated easily.

  The blast furnace 10 is a so-called bell armor type blast furnace that includes a movable armor 11b and that uses a bell charging method as a charging method. Specifically, the blast furnace 10 includes a central charging chute 12 at the top of the blast furnace 10 and a bell type charging device 11 in which a bell 11a and a moveable armor 11b are combined. The central charging chute 12 is a bowl-shaped member extending toward the furnace center in order to concentrate coke into the central portion (furnace center) of the blast furnace 10. The bell type charging device 11 is moved up and down, and the plate 11 reciprocates along the radial direction of the furnace with the bell 11a for charging the charged material (coke and ore) toward the vicinity of the furnace side wall in the furnace. This is a device for throwing the charge into the furnace peripheral edge using the moveable armor 11b that changes the drop position of the charge from the bell 11a.

  The charge, in the present embodiment, coke and ore are charged into the furnace from the top of the blast furnace 10 configured as described above, so that the coke layer and the ore layer are alternately stacked. In this embodiment, a unit of charging coke or ore into the blast furnace 10 is a batch. And 1 coke is comprised mainly by 2 coke batches for forming a coke layer, and 2 ore batches for forming an ore layer. The first coke batch in one charge is C1, the next coke batch is C2, the next ore batch is O1, and the next ore batch is O2. In addition, when forming an ore layer, in order to charge a small amount of coke intensively to the center of the furnace, each center charging to the center of the blast furnace 10 using the center charging dedicated chute performed alternately with the ore batch. The coke batch is CC1 and CC2.

  Specifically, the coke layer is formed by using a bell-type charging device (bell 11a and moveable armor 11b) 11 and performing C1 batch and C2 batch. Specifically, in the C1 batch, the moveable armor 11b moves backward to the furnace side wall position, and a coke dropping point is set in the vicinity of the furnace side wall. In this state, the bell 11a is lowered, so that coke is charged and a coke layer c1 is formed. Next, in the C2 batch, the moveable armor advances to the furnace center side, and a drop point is set from the furnace center. In this state, coke is charged to form a coke layer c2.

  Next, an ore layer is formed. At that time, the central charging coke is also charged into the furnace center. After the coke layers c1 and c2 are formed, the central charging chute 12 and the bell charging device 11 are alternately used. That is, each of the batches (the charge of the charge) in the order of the central charging coke batch CC1 to the furnace center, the ore batch O1 near the furnace side wall, the central charging coke batch CC2 to the furnace center, and the ore batch O2 to the furnace side wall. Charging) is performed. As described above, the centrally charged coke batch at the furnace center and the ore batch near the furnace side wall are alternately performed, so that the centrally charged coke layers cc1 and cc2 are divided by the ore layers o1 and o2. The coke layers c1 and c2 and the centrally charged coke layers cc1 and cc2 are continuous.

  Next, the blast furnace charge layer thickness distribution estimating apparatus used for the blast furnace 10 will be described.

  As shown in FIG. 2, a layer thickness distribution estimation device (hereinafter simply referred to as “estimation device”) 20 is a furnace center of charges such as ore and coke charged and stacked in a blast furnace 10. Is a device for estimating the layer thickness distribution along the furnace radial direction from the furnace side wall to the furnace side wall. Specifically, the estimation device 20 includes a measurement unit 21, a layer thickness distribution estimation unit 22, and an output unit 23.

  As the measuring means 21, a profile meter is used in the present embodiment. The profile meter 21 is a device for measuring the layer upper surface shape (profile) of the charge in the blast furnace 10 and can be measured for each batch in which the charge is charged into the blast furnace 10. As shown in FIG. 3, the profile meter 21 has a measuring rod 211 inserted from the outside of the furnace so as to penetrate the furnace side wall. The measuring part 212 at the tip of the measuring rod 211 reciprocates along the radial direction at a predetermined position of the furnace, and based on the reflected wave of the microwave irradiated from the tip (measuring part) 212, Height is measured.

  The layer thickness distribution estimation means 22 is composed of two computers, a process computer 22a and a layer thickness calculation computer 22b. Note that the layer thickness distribution estimation means 22 does not need to be configured by two computers as in this embodiment, and may be a single computer or may be configured by three or more.

  The process computer 22a controls the profile meter 21 and the bell system charging device 11 and the like, and measures the measured value from the profile meter 21 and the amount of charged material into the blast furnace 10 (volume of charged material). Is transmitted to the layer thickness calculation computer 22b. The process computer 22a is connected to an operation condition input device 24 for inputting operation conditions of the blast furnace 10 and the like.

  The layer thickness calculation computer 22b is a layer thickness distribution of the uppermost layer among the layers stacked in the blast furnace 10 based on information from the process computer 22a (measured value of the profile meter 21, the amount of charged material, etc.). It is a computer that estimates The details of the layer thickness distribution estimation method performed by the layer thickness calculation computer 22b will be described later.

  The output means 23 is for displaying the calculation result in the layer thickness calculation computer 22b, and a CRT display is used in this embodiment. However, the present invention is not limited to this, and may be an FPD, a printer, or the like, or a combination thereof.

  Next, a method for estimating the layer thickness distribution performed by the layer thickness calculation computer 22b will be described with reference to FIGS.

  The layer thickness calculation computer 22b repeats the batch N times (where N is a natural number and 2 <N), and among the charges stacked in the blast furnace 10, the N-1th time measured by the profile meter 21. The layer thickness distribution of the charge of the Nth batch is estimated based on the measured value of the layer upper surface shape of the batch and the Nth batch. The estimated layer thickness distribution is a layer thickness distribution along the furnace radial direction from the furnace center to the furnace side wall in the layer formed by the charge of the Nth batch.

  Specifically, the layer thickness calculation computer 22b estimates the layer thickness distribution in the charge of the Nth batch (hereinafter also simply referred to as “corresponding batch”) as follows.

"Setting a continuous function with an unknown coefficient that defines an estimated shape line"
In the layer thickness calculation computer 22b, a continuous function f with an undetermined coefficient defining the estimated shape line pL at the furnace radius at a predetermined position (position measured by the profile meter 21) is preset and stored. The estimated shape line pL is a curve for imitating the top layer shape of the uppermost layer at the furnace radius at the predetermined position. This estimated shape line pL is obtained as a result of analyzing the layer upper surface shape of the charge in the blast furnace 10 under various operating conditions.

  Specifically, the continuous function f is defined by a plurality of types of section functions f1, f2,... That are defined for a plurality of sections partitioned along the furnace radius at the predetermined position. . The plurality of interval functions f1, f2,... Are set such that the coefficient and the range of the predetermined position in the furnace radial direction are undecided and are continuous at the boundary between adjacent interval functions.

  Specifically, in the present embodiment, this continuous function is defined by three interval functions f1, f2, and f3. These three interval functions f1, f2, and f3 are, in order from the furnace center toward the furnace side wall, a first interval function f1 in which the line segment defined by the corresponding interval function becomes a curve, and the line segment becomes a straight line. It comprises a linear function f2 and a second interval function f3 in which the line segment is a curve. Note that the section in which the first section function f1 is defined is the L1 section, the section defined by the linear function f2 is the L2 section, and the section in which the second section function f3 is defined is the L3 section. Further, in the furnace radius at the predetermined position, the boundary position between the L1 section and the L2 section is r1, and the boundary position between the L2 section and the L3 section is r2.

More specifically, regarding the interval function, the first interval function f1 is represented by y = αx ² + βx + γ (hereinafter also simply referred to as “Expression (10)”), and the linear function f2 is expressed as y = dx + e ( Hereinafter, the second interval function f3 is simply expressed as “expression (11)”, and y = −1 / a · log | cos (αx + b) | + c (hereinafter simply referred to as “expression (12)”). Designated). Here, each section function is a function represented on the xy plane where the furnace center axis is the y axis and the furnace radius at the predetermined position is the x axis. Further, a, b, c, d, e, α, β, and γ are coefficients.

  Each section function is set in this way because the profile of the uppermost layer under various operating conditions in the blast furnace 10 according to the present embodiment is analyzed, and as a result, the section charged by the bell 11a is loaded. This is because the rate of change of the charge height relative to the furnace radius is often constant in the vicinity of the entry fall point. In addition, in the L2 section, it is a part that corresponds to the foot of the mountain due to the charge accumulated in the L3 section, and therefore it is often a straight slope. In the L1 section, the centrally charged coke that is centrally charged. This is because the layers cc1 and cc2 are often deposited on the parabola. Here, the inclination angle of the profile of the L2 section is called a representative inclination angle. This representative inclination angle accumulates in the L3 section when it continues to be charged with the bell 11a and becomes a mountain. However, when the slope of the mountain (bottom) becomes steeper than the representative inclination angle, The angle is such that the collapse is subsided when the tilted surface is broken and becomes the representative tilt angle again.

"Derivation of estimated shape line by calculating coefficients and range of each interval function"
In each of the (N-1) th batch (hereinafter, also simply referred to as “previous batch”) and the Nth batch, two consecutive batches, the previous batch and the corresponding batch, are measured based on the measurement values of the profile meter 21. Each estimated shape line pL of the batch is derived.

  Specifically, in each batch (each of the previous batch and the corresponding batch), the coefficients in all the interval functions f1, f2, and f3 and the ranges in the x-axis direction are simultaneously determined based on the measurement values of the profile meter 21. Calculated. At this time, using the measurement value of the profile meter 21, the coefficients and the ranges in the x-axis direction in all the interval functions are obtained simultaneously by the steepest descent method which is a kind of optimization method. In this way, the estimated shape lines pL1 and pL2 in the corresponding batch and the previous batch are derived.

"Correction of estimated shape line of previous batch"
The two estimated shape lines pL1 and pL2 thus derived are rotated so that the furnace center axis is the center of rotation, thereby obtaining estimated profiles. The estimation of the previous batch is made below the estimated shape line pL2 of the corresponding batch at an interval such that the volume of the space sandwiched between the two estimated profiles thus obtained becomes equal to the total volume of the charge of the corresponding batch. The estimated shape line pL1 of the previous batch is corrected so that the shape line pL1 is positioned (the estimated shape line pL1 is pulled downward in FIG. 5).

"Calculation of layer thickness at each position along the furnace radial direction"
The vertical intervals at the respective positions along the furnace radial direction in the estimated shape lines pL1, pL2 between the previous batch and the corresponding batch corrected in this way are calculated. Based on this calculated value, the layer thickness of the corresponding batch at each position is calculated, and the layer thickness distribution along the furnace radial direction is estimated.

  As described above, the layer thickness distribution is estimated in the layer thickness calculation computer 22b. As described above, the continuous function f whose coefficient is undetermined is set in advance as the estimated shape line pL, and the measured value of the profile meter 21 is obtained. Based on the determination of the coefficients, the laminar pressure distribution is accurately estimated. That is, even if there is a partial protrusion, etc. that happened when charging the measured profile, the smooth estimated shape line pL is derived without being affected by this partial protrusion, etc. An accurate profile can be obtained. Therefore, the interval between the two estimated profiles is such that the volume of the space sandwiched between the two estimated profiles of the N-1st batch and the Nth batch is the same as the volume of the charge of the Nth batch. And the estimated shape line pL1 of the (N-1) th batch is corrected, the vertical distance at each position along the furnace radial direction of both estimated shape lines pL1, pL2 is close to the actual charge layer. Therefore, an accurate layer thickness distribution can be obtained.

  Further, the continuous function f that defines the estimated shape line pL is divided into a plurality of sections, and section functions f1, f2,... That define line segments that conform to the shape of the upper surface of the layer are set for each of the plurality of sections. Thus, it is possible to obtain an estimated shape line pL that more closely approximates the actual charge profile. Therefore, the layer thickness distribution in the charge of the Nth batch is estimated with higher accuracy.

  In addition, in this embodiment, since the continuous function is defined by three interval functions, it is easy to calculate the coefficient and the range in the x-axis direction in each interval function based on the measured value of the profile, and the estimated shape line pL Derivation becomes easy. In addition, since the derived estimated shape line pL is a shape conforming to the profile of the charge stacked in the blast furnace according to the present embodiment, the layer thickness distribution of the charge is accurately estimated. .

  Moreover, since the coefficient in each interval function and the range in the x-axis direction are calculated based on the measured value (actual value) of the profile, the estimation corresponding to the change in the furnace condition in the blast furnace under the same operating conditions The shape line pL can be calculated.

  Next, the operation of the estimation device 20 will be described using the operation when estimating the layer thickness distribution of the O1 batch along the furnace radial direction with reference to FIG. The layer thickness distribution of the O1 batch along the furnace radius refers to the distribution of the layer thickness of the O1 batch at each position from the furnace center to the furnace side wall along the predetermined path radius. Is the vertical distance (interval) from the profile of the C2 batch (coke layer c2) to the profile of the O1 batch (ore layer o1) at a certain position on the furnace radius.

  First, whether or not the measurement value by the profile meter 21 for the profile (layer upper surface shape) of the C2 batch exists (stores) in the storage area of the layer thickness calculation computer 22b or the process computer 22a. It is determined whether or not (step S1 in FIG. 6).

  If the measured value of the profile of the C2 batch does not exist, the layer thickness calculation computer 22b gives the profile of the C2 batch profile, that is, the top layer shape of the uppermost layer (coke layer c2) to the process computer 22a. A command is issued to obtain the measured value. In response to this command, the process computer 22a controls the profile meter 21 to measure the profile of the C2 batch and obtain the measured value. At this time, the profile meter 21 measures the distance from the furnace center to the furnace side wall at intervals of 20 cm along the furnace radius at a predetermined position. In the present embodiment, since the furnace radius is 5 m, the measurement values of the heights of the 26 layer upper surfaces are obtained at equal intervals in the furnace radius (see FIGS. 4 and 5). When the measurement value of the profile of the C2 batch exists, the measurement of the profile of the O1 batch described later is started.

  Next, the layer thickness calculation computer 22b acquires the measurement value of the profile of the C2 batch from the process computer 22a (step S2 in FIG. 6).

  Then, the layer thickness calculation computer 22b calculates the estimated shape line pL of the C2 batch based on the measurement value (step S3 in FIG. 6).

  At this time, the coefficients and the ranges in the x-axis direction of all the section functions of the estimated shape line pL are calculated as follows based on the measured values of the profile of the C2 batch.

  In the vicinity of the boundary position r2 between the L2 section and the L3 section, the actual profile inclination angle is equal to the representative inclination angle (see FIG. 4). Therefore, it is assumed that the estimated shape line pL has not only continuity but also differentiability at the boundary position r2. In this way, when the coefficients a, b, and c in the above equation (12) and the coefficient d in the above equation (11) are determined, the coefficient e and the L2 interval in the above equation (11) and L3 The boundary position r2 of the section is determined dependently. Further, since the profile is continuous even at the boundary position r1 between the L1 section and the L2 section, the continuity of the function can be assumed, and the coefficients a, b, c in the above expression (12), and the above expression (11) ), The coefficients α and β in the above equation (10), and the boundary position r1 between the L1 and L2 intervals, the coefficient γ in the above equation (10) can be determined in a dependent manner. In this way, the coefficients a, b, c in the above equation (12), the coefficient d in the above equation (11), the coefficients α, β in the above equation (10), and the boundary position between the L1 and L2 intervals. r1 is a decision variable, the coefficient e in (11) above, the boundary position r2 between the L2 and L3 sections, and the coefficient γ in the above equation (10) are dependent variables.

Next, a furnace radius position at an arbitrary measurement point is set to x _i , a measured value of an actual profile at the furnace radius position is set to y _i, and a continuous function f defining an estimated shape line pL at the furnace radius x is y = f (X)
Σ {y _i −f (x _i )} ²
Of the continuous function that defines the estimated shape line (coefficients a, b, c in the above equation (12), coefficient d in the above equation (11), and the above equation (10). Coefficients α, β) and a section where the function is defined (boundary position r1 between the L1 section and the L2 section) are simultaneously obtained by the steepest descent method which is one of nonlinear optimization techniques. By determining the parameters in this way, the coefficient e in the above equation (11), the coefficient γ in the above equation (10), and the boundary position r2 between the L2 and L3 intervals, which are the dependent variables, can be determined. . That is, the range in the x-axis direction of all parameters and each interval function can be determined by the steepest descent method.

  The estimated shape line pL1 of the C2 batch is calculated by substituting each coefficient thus obtained and the range in the x-axis direction into the above equations (10) to (12).

  Next, the O1 batch is performed, and thereafter, the layer thickness calculation computer 22b instructs the process computer 22a to obtain the measurement value by the profile meter 21 for the profile of the O1 batch. By this command, the process computer 22a controls the profile meter 21 to measure the profile of the O1 batch and obtain the measured value.

  Then, the layer thickness calculation computer 22b acquires the measured value for the profile of the O1 batch from the process computer 22a (step S4 in FIG. 6), and similarly to the estimated shape line of the O1 batch based on the acquired measured value. pL2 is calculated (step S5 in FIG. 6).

  The layer thickness calculation computer 22b acquires the volume of the entire charged material (ore) charged in the O1 batch from the process computer 22a (step S6 in FIG. 6).

  The layer thickness calculation computer 22b rotates the two estimated shape lines pL1 and pL2 of the C2 batch and the O1 batch derived (calculated) based on the measurement value of the profile meter 21 so that the furnace center axis is the rotation center. The estimated profiles of the C2 batch and the O1 batch are calculated respectively. The C2 batch below the estimated shape line pL2 of the O1 batch at an interval such that the volume of the space sandwiched between the two estimated profiles is equal to the total volume of the O1 batch charge obtained from the process computer 22a. The estimated shape line of the C2 batch is corrected so that the estimated shape line pL1 is positioned (see FIG. 5) (step S7 in FIG. 6).

  In this way, by calculating the vertical spacing at each position in the x-axis direction of the pair of estimated shape lines pL2, pL1 arranged vertically, the layer thickness at each position is calculated (estimated). The layer thickness distribution is calculated (estimated) by arranging the layer thicknesses at the respective positions calculated in this way along the x-axis direction (step S8 in FIG. 6).

  The layer thickness calculation computer 22b displays the layer thickness distribution thus calculated on the CRT (output means) 23 (step S9 in FIG. 6).

  The layer thickness calculation computer 22b determines whether to continuously estimate the layer thickness distribution of the next batch (O2 batch in the present embodiment) (step S10 in FIG. 6). When estimating the layer thickness distribution, the above series of steps is repeated. In this case, since the estimated shape line pL2 of the O1 batch has already been calculated, the layer thickness calculation computer 22b starts from the step of obtaining the measured value of the profile of the O2 batch and calculating the estimated shape line of the O2 batch. That's fine. If the layer thickness distribution of the next batch is not estimated, the process ends.

  The operating conditions of the blast furnace according to the above embodiment were made the same, and the layer thickness distribution of a predetermined batch was estimated at two different dates (date 1 and date 2). At that time, the profile of two consecutive batches at each date and time is measured with a profile meter, and based on the measured value obtained by this measurement, the layer thickness distribution estimation method according to the above embodiment (the estimated shape defined by the continuous function) A method using a line) and a case where the layer thickness distribution is estimated without using the layer thickness distribution estimation method (using the unprocessed measurement value), The results are shown in FIG. In FIG. 7, the method using the layer thickness distribution estimation method according to the above embodiment is “this method”, and the method not used is “unprocessed”.

  The estimation results of the layer thickness distribution of the predetermined batch on the date 1 are “present method 1” and “unprocessed 1” in FIG. The measurement results of the layer thickness distribution of a predetermined batch at the date 2 are “present method 2” and “unprocessed 2” in FIG. In this way, the profile shape changes even if the operating conditions are the same. That is, it can be seen that, unlike the conventional simulation, it is possible to estimate the layer thickness distribution that accurately reflects the change in the layer thickness distribution due to the change in the furnace condition. In addition, since the change in the layer thickness in the furnace radial direction is smoother in the one using this method than in the untreated one, it is affected by unevenness such as partial protrusions that exist only in the measured profile. It can be seen that the layer thickness distribution was estimated with high accuracy without being affected or with little influence.

  Next, the profile of two continuous batches (C2 batch and O1 batch) in the blast furnace according to the above embodiment is measured with a profile meter, and the estimated profile of each batch is calculated from the measured value, and the volume between both estimated profiles is calculated. Was corrected to lower the estimated shape line of the previous batch (C2 batch) or the profile in the furnace radial direction at a predetermined position so that the total volume of the charge of the corresponding batch (O1 batch) becomes equal.

  FIG. 8A to FIG. 8C show the results. FIG. 8A shows a profile in the furnace radial direction at a predetermined position where the measurement values are used unprocessed, and FIG. 8B shows the above embodiment from the same batch measurement values as in FIG. 8A. FIG. 8C shows an estimated shape line (profile in the furnace radial direction at a predetermined position) calculated by the layer thickness distribution estimation method, and FIG. 8C shows 30 averages of the C2 and O1 batches, and the averaged furnace radius at the predetermined position. The profile in is shown. The profile shown in FIG. 8 (c) obtained by averaging 30 batches is a profile at a time when the operating conditions such as the charging amount and the armor stroke are constant and the situation in the blast furnace is hardly changed. It is considered that there is almost no noise generated only in one batch (unevenness of only one batch profile).

  From these FIGS. 8 (a) to 8 (c), by using the layer thickness distribution estimation method according to the above-described embodiment, a profile (estimated shape) having almost no noise or the like can be obtained only by measuring two consecutive batch profiles. It can be seen that a profile with a small measurement error is obtained (see FIG. 8B and FIG. 8C).

  Fig.9 (a) is a figure which shows the layer thickness of O1 batch shown by Fig.8 (a) thru | or FIG.8 (c). From this figure, the measured layer thickness distribution of each batch (each measured value of the profile of two consecutive batches) was used unprocessed. On the other hand, there is a large deviation at the furnace radii of 1.2m and 3.4m. On the other hand, it can be seen that the profile (estimated shape line) using the layer thickness distribution estimation method according to the above embodiment is close to the layer thickness of the profile obtained by averaging 30 batches.

  FIG. 9B shows the case where the estimated shape line is calculated and the layer thickness distribution is estimated with respect to the layer thickness of the profile obtained by averaging 30 batches (when the layer thickness distribution estimation method according to the above embodiment is used). And the layer thickness error when the layer thickness distribution is estimated using the measured values as they are (when the measured values are used unprocessed). As can be seen from this figure, it can be seen that the layer thickness distribution using the layer thickness distribution estimation method according to the above embodiment has fewer errors. Thus, the thickness distribution can be estimated (obtained) with high accuracy by measuring the profile of each of the upper and lower batches (two consecutive batches).

  Next, the section functions of the L1 section and the L3 section in the continuous function f are interchanged, that is, the section function of the L1 section is set as the second section function f3, and the section function of the L3 section is set as the first section function f1. , Using the estimation device 20 with the other being the same, the profile of two continuous batches (C2 batch and O1 batch) in the blast furnace according to the above-described embodiment is measured with a profile meter, and each batch is determined from the measured values. The estimated profile was calculated, and correction was made to lower the estimated shape line of the previous batch or the profile in the furnace radial direction at a predetermined position so that the volume between both estimated profiles is equal to the volume of the entire charge of the corresponding batch.

  FIG. 10A to FIG. 10C show the results. 10 (a) is the same as FIG. 8 (a), and FIG. 10 (c) is the same as FIG. 8 (c). FIG. 10B shows an estimation calculated by the layer thickness distribution estimation method in which the interval functions of the L1 interval and the L3 interval of the continuous function f are interchanged in the above embodiment from the same batch measurement values as in FIG. A shape line is shown.

  From these FIG. 10 (a) to FIG. 10 (c), even by using the layer thickness distribution estimation method in which the section functions of the L1 section and the L3 section of the continuous function f in the above embodiment are replaced, It can be seen that a profile with a small measurement error can be obtained just by measuring the profile (see FIGS. 10B and 10C).

  Fig.11 (a) is a figure which shows the layer thickness of O1 batch shown by Fig.10 (a) thru | or FIG.10 (c). From this figure, the thickness of the profile in which the measured values for determining the layer thickness distribution of one batch were used untreated is the position of the furnace radius of 1.2 m and 3.4 m with respect to the layer thickness of the profile obtained by averaging 30 batches. It is greatly shifted by. On the other hand, it can be seen that the profile using the layer thickness distribution estimation method in which the interval functions of the L1 interval and L3 interval of the continuous function f are interchanged is close to the layer thickness of the profile obtained by averaging 30 batches.

  Further, FIG. 11B shows the layer thickness distribution using the measured thickness and the error of the layer thickness when the estimated shape line is calculated and the layer thickness distribution is estimated with respect to the layer thickness of the profile averaged over 30 batches. It shows the error of the layer thickness when. As can be seen from this figure, even in the layer thickness distribution using the layer thickness distribution estimation method in which the interval functions of the L1 interval and the L3 interval of the continuous function f are interchanged in the above embodiment, there are few errors in the same manner as described above. I understand. As described above, by using the layer thickness distribution estimation method in which the interval functions of the L1 interval and the L3 interval are interchanged in the above embodiment, it is possible to measure the accuracy of each of the upper and lower batches (two consecutive batches). The layer thickness distribution can be estimated well.

  It should be noted that the layer thickness distribution estimation method and apparatus of the present invention are not limited to the above-described embodiment, and various modifications can be made without departing from the scope of the present invention.

  For example, in this embodiment, the continuous function that defines the estimated shape line is defined by three interval functions, but is not limited to this, and may be two interval functions. The interval function described above may be used. Further, the interval function need not be limited to the above formulas (10) to (12). That is, it is only necessary that a plurality of types of interval functions are connected to form a curve that can accurately simulate the top surface shape (profile) of the uppermost layer of the furnace radius at a predetermined position. By doing so, it is possible to accurately estimate the layer thickness distribution of the corresponding batch only by measuring the profile of two consecutive batches as described above.

  Accordingly, it is not necessary to limit to the blast furnace of the central charging method as in the present embodiment, and various charging methods can be obtained simply by measuring the profile of two consecutive batches by appropriately selecting the number of sections and the section function. It is possible to accurately estimate the layer thickness distribution in a blast furnace (for example, a bell-less system in which a charge is charged by a chute).

  Moreover, in this embodiment, although an estimated shape line is prescribed | regulated by said Formula (10) thru | or (12) continuing in order, it is not necessary to be limited to this order. That is, in the above L1 to L3 sections, if the L2 section is the formula (11), the L1 section may be the formula (10) or (12), and the L3 section may be the formula (10) or (12). That's fine. Therefore, the expression (10), the expression (11), and the expression (10) may be sequentially connected in order from the L1 section, or the expression (12), the expression (11), and the expression (10) may be sequentially connected in this order.

  In this embodiment, the charge to the blast furnace is coke and ore. However, the charge is coke or ore in which auxiliary materials (limestone, silica, converter slag, etc.) are mixed. May be. However, in the case where a large amount of the auxiliary material of a predetermined amount or more is charged only with the auxiliary material, the layer thickness estimation method according to the above embodiment is applied to the auxiliary material layer formed by the auxiliary material. be able to.

It is a partial expanded block diagram of the blast furnace top part concerning this embodiment. It is a block diagram of the layer thickness distribution estimation apparatus which concerns on the same embodiment. It is the schematic of the profile meter of the layer thickness distribution estimation apparatus which concerns on the same embodiment. It is a figure which shows the estimated shape line and interval function in the layer thickness distribution estimation method which concerns on the embodiment. It is a figure which shows correction | amendment of the estimated shape line of the last batch in the layer thickness distribution estimation method which concerns on the embodiment. It is a figure which shows the flow of operation | movement in the layer thickness distribution estimation apparatus which concerns on the same embodiment. It is a figure which shows the comparison with the layer thickness computed by the layer thickness distribution estimation method which concerns on this embodiment on the same operation conditions, and the layer thickness computed by un-processing in the blast furnace concerning the embodiment. In the estimation of the layer thickness distribution in the blast furnace according to the embodiment, (a) shows the profile of the actual measurement value, (b) shows the profile obtained by the layer thickness distribution estimation method according to the embodiment, (c) FIG. 3 is a diagram showing an average profile of 30 batches. In the estimation of the layer pressure distribution in the blast furnace according to the embodiment, (a) shows the layer thickness in each profile in FIG. 8 (a) to FIG. 8 (c), and (b) in FIG. 9 (a). It is a figure which shows the difference | error of layer thickness. In the estimation of the layer thickness distribution in the blast furnace according to another embodiment, (a) shows the profile of the actual measurement value, (b) shows the profile obtained by the layer thickness distribution estimation method according to the embodiment, and (c) FIG. 3 is a diagram showing an average profile of 30 batches. In the estimation of the layer pressure distribution in the blast furnace according to the embodiment, (a) shows the layer thickness in each profile in FIG. 10 (a) to FIG. 10 (c), and (b) in FIG. 11 (a). It is a figure which shows the difference | error of layer thickness.

Explanation of symbols

10 Blast furnace 21 Profile meter (measuring means)
f Continuous function pL Estimated shape line

Claims (9)

Of the charges stacked in the blast furnace by repeating the batch of charging the charges into the blast furnace N times (where N is a natural number and N> 1), the nth time (where n is A natural number arbitrarily selected within the range of 1 <n ≦ N) A method for estimating the layer thickness distribution along the furnace radial direction from the furnace center to the furnace side wall in the layer formed by the batch charge,
Set a continuous function with an unknown coefficient that defines an estimated shape line to simulate the shape of the top surface of the charge at the furnace radius at a predetermined position,
In each of the (n−1) -th batch and the n-th batch, the upper surface shape of the layer is measured by a measuring unit, and the coefficient of the continuous function is calculated based on the measured value obtained by the measurement. Deriving each estimated shape line of the (n-1) th batch and the nth batch,
When the estimated profiles are obtained by rotating these two estimated shape lines so that the furnace center axis is the center of rotation, the volume of the space sandwiched between the two estimated profiles is the charge of the nth batch. At least one of the two estimated shape lines is corrected 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 object volume. And
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 arranged in this way is calculated,
A method for estimating a layer thickness distribution of a blast furnace charge, comprising estimating a layer thickness distribution along a furnace radial direction from a furnace center to a furnace side wall in the charge of the nth batch based on the calculation result. In the method for estimating the layer thickness distribution of the blast furnace charge according to claim 1,
The continuous function is defined by connecting a plurality of types of section functions defined for a plurality of sections partitioned along the furnace radius at the predetermined position,
The plurality of interval functions are set such that the coefficient and the range of the predetermined position in the furnace radial direction are undetermined and are continuous at the boundary between adjacent interval functions,
In each of the (n−1) th batch and the nth batch, the coefficient in all the interval functions and the range in the furnace radial direction of the predetermined position are calculated based on the measured value, thereby the n−th batch. A method for estimating a layer thickness distribution of a blast furnace charge, wherein each estimated shape line of the first batch and the n-th batch is derived. In the method for estimating the layer thickness distribution of the blast furnace charge according to claim 1 or 2,
The method for estimating the layer thickness distribution of a blast furnace charge, wherein N is N = n and N> 2. In the method for estimating the layer thickness distribution of the blast furnace charge according to claim 2 or 3,
The method of estimating a layer thickness distribution of a blast furnace charge, wherein the interval function is a function represented on an xy plane in which the furnace center axis is the y axis and the furnace radius at the predetermined position is the x axis. In the method for estimating the layer thickness distribution of the blast furnace charge according to claim 4,
The continuous function is defined by three interval functions,
These three section functions are, in order from the furnace center toward the furnace side wall, a first section function in which the line segment defined by the corresponding section function is a curve, a linear function in which the line segment is a straight line, and the line segment Is a second interval function in which the curve is a curve. In the method for estimating the layer thickness distribution of the blast furnace charge according to claim 5,
The first interval function and the second interval function are a function having a constant angle change rate with respect to the x-axis expressed by the following equation (1), or a secondary expressed by the following equation (2). A method for estimating the layer thickness distribution of a blast furnace charge, which is a function.
y = −1 / a · log | cos (αx + b) | + c (1)
y = αx ² + βx + γ (2)
Here, a, b, c, α, β, and γ are coefficients. In the method for estimating the layer thickness distribution of the blast furnace charge according to any one of claims 4 to 6,
A blast furnace characterized in that, in each of the (n-1) -th batch and the n-th batch, a coefficient in all the interval functions and a range in the x-axis direction are obtained simultaneously from the measured value using the steepest descent method. Method for estimating the layer thickness distribution of charges. Measuring means capable of measuring the top surface shape of the charge in the blast furnace for each batch in which the charge is charged into the blast furnace;
Of the charges stacked in the blast furnace by repeating the batch N times (where N is a natural number and 1 <N), the measurement means uses the n-1st batch and the nth batch. Based on the measurement values obtained by measuring the top surface shape of the two consecutive batches, the nth (where n is a natural number arbitrarily selected within the range of 1 <n ≦ N) in the layer formed by the charge of the batch A layer thickness distribution estimating means for estimating a layer thickness distribution along the furnace radial direction from the furnace center to the furnace side wall;
An output means for outputting the value of the layer thickness distribution estimated by the layer thickness measuring means to the outside, and a layer thickness distribution estimating apparatus for a blast furnace charge,
The layer thickness distribution estimating means sets a continuous function with an unknown coefficient defining an estimated shape line for imitating the shape of the layer upper surface of the charge at a furnace radius at a predetermined position,
In each of the (n−1) -th batch and the n-th batch, the upper surface shape of the layer is measured by a measuring unit, and the coefficient of the continuous function is calculated based on the measured value obtained by the measurement. Deriving each estimated shape line of the (n-1) th batch and the nth batch,
When the estimated profiles are obtained by rotating these two estimated shape lines so that the furnace center axis is the center of rotation, the volume of the space sandwiched between the two estimated profiles is the charge of the nth batch. At least one of the two estimated shape lines is corrected 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 object volume. And
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 arranged in this way is calculated,
This calculated value is transmitted to the output means as the value of the layer thickness distribution along the furnace radial direction from the furnace center to the furnace side wall in the charge of the nth batch,
The apparatus for estimating a layer thickness distribution of a blast furnace charge, wherein the output unit outputs the value of the layer thickness distribution transmitted from the layer thickness measuring unit to the outside. In the apparatus for estimating the layer thickness distribution of the blast furnace charge according to claim 9,
N is N = n and N> 2.
The apparatus for estimating a layer thickness of a blast furnace charge, wherein the measuring means is a profile meter capable of measuring the shape of the upper surface of the charge of the uppermost layer.

Images