JPH039162B2 - - Google Patents
Info
- Publication number
- JPH039162B2 JPH039162B2 JP4215387A JP4215387A JPH039162B2 JP H039162 B2 JPH039162 B2 JP H039162B2 JP 4215387 A JP4215387 A JP 4215387A JP 4215387 A JP4215387 A JP 4215387A JP H039162 B2 JPH039162 B2 JP H039162B2
- Authority
- JP
- Japan
- Prior art keywords
- value
- prediction means
- furnace
- temperature
- evaluation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired
Links
- 238000011156 evaluation Methods 0.000 claims description 47
- 238000000034 method Methods 0.000 claims description 26
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 4
- 229910052799 carbon Inorganic materials 0.000 claims description 4
- 230000007423 decrease Effects 0.000 description 28
- 239000002184 metal Substances 0.000 description 24
- 229910052751 metal Inorganic materials 0.000 description 24
- 239000007789 gas Substances 0.000 description 17
- 230000008569 process Effects 0.000 description 14
- 238000005070 sampling Methods 0.000 description 14
- 230000009471 action Effects 0.000 description 9
- 238000010438 heat treatment Methods 0.000 description 8
- 230000002159 abnormal effect Effects 0.000 description 7
- 230000008859 change Effects 0.000 description 7
- 238000012545 processing Methods 0.000 description 7
- 238000004817 gas chromatography Methods 0.000 description 6
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 4
- 210000001015 abdomen Anatomy 0.000 description 4
- 238000009529 body temperature measurement Methods 0.000 description 4
- 238000009434 installation Methods 0.000 description 4
- 239000002994 raw material Substances 0.000 description 4
- 230000005856 abnormality Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N iron Substances [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 108010077544 Chromatin Proteins 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 229910002091 carbon monoxide Inorganic materials 0.000 description 2
- 210000003483 chromatin Anatomy 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000945 filler Substances 0.000 description 2
- 230000017525 heat dissipation Effects 0.000 description 2
- 229910052757 nitrogen Inorganic materials 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 238000009423 ventilation Methods 0.000 description 2
- HSRJKNPTNIJEKV-UHFFFAOYSA-N Guaifenesin Chemical compound COC1=CC=CC=C1OCC(O)CO HSRJKNPTNIJEKV-UHFFFAOYSA-N 0.000 description 1
- 240000001549 Ipomoea eriocarpa Species 0.000 description 1
- 235000005146 Ipomoea eriocarpa Nutrition 0.000 description 1
- 239000011449 brick Substances 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 239000011810 insulating material Substances 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 229910052745 lead Inorganic materials 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000000520 microinjection Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 238000010079 rubber tapping Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000003466 welding Methods 0.000 description 1
Landscapes
- Manufacture Of Iron (AREA)
- Blast Furnaces (AREA)
Description
【発明の詳細な説明】
(産業上の利用分野)
この発明は、高炉の安定な操業を行なうための
高炉炉熱低下予測方法に関する。DETAILED DESCRIPTION OF THE INVENTION (Industrial Field of Application) The present invention relates to a method for predicting heat drop in a blast furnace for stable operation of a blast furnace.
(従来の技術とその問題点)
高炉の安定操業の維持のためには、溶銑温度を
一定にすることが必要であることが従来より知ら
れている。このため、高炉操業者は常に高炉炉熱
変化を予測する必要性があつた。(Prior art and its problems) It has been known for a long time that in order to maintain stable operation of a blast furnace, it is necessary to keep the temperature of hot metal constant. For this reason, blast furnace operators have always needed to predict changes in blast furnace heat.
高炉炉熱変化において、特に温度低下によつて
溶銑が凝固し、高炉から流出しなくなる可能性が
あるため、温度低下の予測は極めて重要なものと
なる。 Prediction of temperature drop is extremely important in blast furnace furnace thermal changes, especially since hot metal may solidify due to temperature drop and may no longer flow out of the blast furnace.
高炉炉熱の予測方法としては、特開昭60−
39107に開示されたものがある。この方法は炉腹
部周辺装入物温度が溶銑温度と強い相関関係を持
つという見地から、予め第14図に示す様に高炉
1に設置したセンサ(炉腹ゾンデ)2により検出
される炉腹部周辺部温度と溶銑温度との関係を第
15図に示す如く直線回帰する。この直線式に基
づき、炉腹部周辺部温度から溶銑温度Tpigを予測
するのである。 As a method for predicting blast furnace furnace heat,
There is one disclosed in 39107. This method is based on the viewpoint that the charge temperature around the blast furnace belly has a strong correlation with the hot metal temperature. The relationship between the temperature of the hot metal and the temperature of the hot metal is linearly regressed as shown in FIG. Based on this linear equation, the hot metal temperature T pig is predicted from the temperature around the furnace belly.
しかしながら、この方法では炉内の内壁近傍の
温度を測定するために炉腹ゾンデ2を挿入する必
要があり、このため温度測定を間欠時点でしか行
なえず、溶銑温度予測精度も当然悪化してしまう
という問題点があつた。 However, with this method, it is necessary to insert the furnace probe 2 in order to measure the temperature near the inner wall of the furnace, which means that temperature measurements can only be carried out intermittently, which naturally deteriorates the accuracy of hot metal temperature prediction. There was a problem.
また、溶銑温度が同じ値でも、生産計画や原料
装入条件等の変化により、炉内温度が変化する場
合がある。したがつて第15図で示した炉壁温度
の絶対値に基づく直線式では、必ずしも正確な予
測ができないという問題点があつた。 Further, even if the hot metal temperature is the same, the furnace temperature may change due to changes in the production plan, raw material charging conditions, etc. Therefore, the linear equation based on the absolute value of the furnace wall temperature shown in FIG. 15 has the problem that accurate prediction cannot necessarily be made.
一方、従来より、高炉の還元状態の良否を示す
ソリユーシヨンロスカーボン量(以下「ソルロス
C量」と言う)の増減により、高炉炉熱温度の予
測が別の予測方法として行なわれている。ソルロ
スC量の増加は、以下に示すいわゆるソルロス反
応が促進することを示している。 On the other hand, another prediction method has conventionally been used to predict the blast furnace thermal temperature based on an increase or decrease in the amount of solution loss carbon (hereinafter referred to as "sol loss C amount"), which indicates the quality of the reduction state of the blast furnace. An increase in the amount of Solloss C indicates that the so-called Solloss reaction described below is promoted.
C+CO2→2CO
この反応は、吸熱反応であるため高炉炉熱が低
下することが予測できる。 C+CO 2 →2CO Since this reaction is an endothermic reaction, it can be predicted that the blast furnace heat will decrease.
ソルロスC量は、通常炉頂ガスの組成を分析す
るガスクロマトグラフイーの分析周期(3分程
度)毎に、灯頂ガス中のCO、CO2、N2等の割合
や送風条件や原料装入条件をもとに計算され、従
来は1時間毎のソルロスC量の平均値により炉熱
低下を管理していた。 The amount of Sol loss C is determined by determining the proportion of CO, CO 2 , N 2, etc. in the lamp top gas, ventilation conditions, and raw material charging every analysis cycle (about 3 minutes) of gas chromatography, which normally analyzes the composition of the furnace top gas. It is calculated based on the conditions, and conventionally, the decrease in furnace heat was managed by the average value of the Sol loss C amount every hour.
第16図a,bにおいて、同図aは、3分毎の
ソルロスC量l1、1時間毎のソルロスC量平均
値l2の経時変化を示し、同図bは溶銑温度の経
時変化を示すグラフである。同図において、時刻
17時に閾値ε1を越えているが昇熱アクシヨンをと
らず、その後のように溶銑温度は大幅低下して
いる。第17図は、閾値ε2を越えた時刻13時に昇
熱アクシヨンAを起した時の各々の経時変化を示
している。なお第16図と同様、図中l1が3分
毎の瞬時値、l2が1時間平均のソルロスC量を
示している。第16図、第17図を比較すること
により、昇熱アクシヨンAにより第16図bの
のような溶銑温度低下が、第17図で示すように
ある程度回避できているのがわかる。 In Figures 16a and 16b, Figure a shows the change over time in the Sol loss C amount l1 every 3 minutes and the average Sol loss C amount l2 every hour, and Figure b shows the change over time in the hot metal temperature. It is. In the figure, time
Although the threshold value ε 1 is exceeded at 17:00, no heating action is taken, and the hot metal temperature drops significantly after that. FIG. 17 shows the respective changes over time when the heating action A was initiated at time 13:00 when the threshold value ε 2 was exceeded. Note that, similarly to FIG. 16, l1 in the figure indicates an instantaneous value every 3 minutes, and l2 indicates an hourly average Sol loss C amount. By comparing FIG. 16 and FIG. 17, it can be seen that the temperature increase of the hot metal as shown in FIG. 16b can be avoided to some extent by the heating action A, as shown in FIG. 17.
しかしながら、1時間毎のソルロスC量の平均
値の予測では、急激なソルロスC量の増加があつ
た時、最悪の場合、ほぼ1時間程度も炉熱低下の
予測に遅れが生じてしまう問題点があつた。例え
ば第17図の場合にしても、予測遅れのため昇熱
アクシヨンAをとるのは溶銑温度が管理温度Tc
をある程度下まわつてからになつてしまつてい
る。そこで、この問題点を回避するため、3分程
度の間隔で測定したソルロスC量の瞬時値で炉熱
低下予測した場合、第16図、第17図のl1で
示したように個々のバラツキが大きく、ノイズ成
分が大きいためデータの持続性がない。したがつ
てソルロスC量の瞬時値では炉熱低下予測が不可
能に近い。 However, the problem with predicting the average value of Sol loss C amount every hour is that when there is a sudden increase in Sol loss C amount, in the worst case, there is a delay in predicting the decrease in furnace heat by approximately one hour. It was hot. For example, in the case of Fig. 17, the reason why the heating action A is taken due to the prediction delay is that the molten metal temperature is the control temperature T c
It has become a problem since it has fallen below a certain level. Therefore, in order to avoid this problem, if the furnace heat drop is predicted using the instantaneous value of the Sol loss C amount measured at intervals of about 3 minutes, individual variations will be reduced as shown by l1 in Figs. 16 and 17. The data is not sustainable because it is large and the noise component is large. Therefore, it is almost impossible to predict the decrease in furnace heat using the instantaneous value of the Sol loss C amount.
(発明の目的)
この発明の目的は、上記従来技術の問題点を解
消し、予測ができるたけ早く得られ、しかも溶銑
温度の低下を正確に予測することのできる高炉炉
熱低下予測方法を提供することである。(Object of the Invention) The object of the present invention is to provide a blast furnace furnace heat drop prediction method that can solve the problems of the prior art described above, obtain predictions as quickly as possible, and accurately predict the drop in hot metal temperature. It is to be.
(目的を達成するための手段)
上記目的を達成するため、この発明における高
炉炉熱低下予測方法は、高炉の所定箇所に内壁温
度計を設置し、該内壁温度計にて、所定時間間隔
ごとの内壁温度差を測定し、ある時刻における前
記内壁温度差の正の値を示す部分の合計値が閾値
を越えた時に所定期間評価点を与える第1の予測
手段と、ある時刻における前記内壁温度差の負の
値を示す部分の合計値が閾値を越えた時に所定期
間評価点を与える第2の予測手段と、ある時刻に
おける前記内壁温度差の正の値を示す部分の時間
幅の移動平均値の総和が閾値を越えた時に所定期
間評価点を与える第3の予測手段とのうちの少な
くとも1つ備え、かつ、ソリユーシヨンロスカー
ボン量を所定時間間隔ごとに求め、この求めた値
の所定時間幅における移動平均値が閾値を越えた
時に評価点を与える第4の予測手段と、炉頂ガス
成分中の窒素量を所定時間間隔ごとに求め、この
求めた値の所定時間幅における移動平均値が閾値
を下回つた時に評価点を与える第5の予測手段の
うちの少なくとも1つをさらに備え、前記第1〜
第3の予測手段のうちの少なくとも1つの評価点
と前記第4、第5の予測手段のうちの少なくとも
1つの評価点の総合評価に従い高炉炉熱低下予測
を行なうようにしている。(Means for Achieving the Object) In order to achieve the above object, the blast furnace furnace heat drop prediction method according to the present invention includes installing an inner wall thermometer at a predetermined location of the blast furnace, and measuring the temperature at predetermined time intervals with the inner wall thermometer. a first prediction means that measures the inner wall temperature difference of and gives an evaluation point for a predetermined period when the total value of the portion showing a positive value of the inner wall temperature difference at a certain time exceeds a threshold; a second prediction means for giving an evaluation point for a predetermined period when the total value of the portion showing a negative value of the difference exceeds a threshold; and a moving average of the time width of the portion showing a positive value of the inner wall temperature difference at a certain time. and a third prediction means that gives an evaluation point for a predetermined period when the sum of the values exceeds a threshold, and calculates the amount of solution loss carbon at each predetermined time interval, and calculates the amount of solution loss carbon at each predetermined time interval. A fourth prediction means that gives an evaluation point when the moving average value in a predetermined time width exceeds a threshold value, and a fourth prediction means that calculates the amount of nitrogen in the furnace top gas component at each predetermined time interval and moves the determined value in the predetermined time width. Further comprising at least one of fifth prediction means for giving an evaluation point when the average value falls below the threshold,
The blast furnace heat drop is predicted according to the comprehensive evaluation of at least one evaluation point of the third prediction means and at least one evaluation point of the fourth and fifth prediction means.
(実施例)
A 第1の炉熱低下理由
高炉の炉熱低下の一因として、以下に示すも
のが考えられる。(Example) A First Reason for Reduction in Furnace Heat The following factors can be considered as causes for the decrease in furnace heat in a blast furnace.
高炉羽口から吹き上げる溶銑温度及び溶銑量
調整のための高温空気(ガス流)は通常、炉内
中央部に吹き込んでいる。ところが原料装入条
件、装入物分布等の理由により、急にガス流が
炉内周辺部に多く流れる場合がある。その結
果、
FeO+C→Fe+CO
の吸熱反応が促進され、炉熱低下が起こる。 High-temperature air (gas flow) blown up from the blast furnace tuyeres to adjust the temperature and amount of hot metal is usually blown into the center of the furnace. However, due to reasons such as raw material charging conditions and charge distribution, a large amount of gas may suddenly flow to the periphery of the furnace. As a result, the endothermic reaction of FeO+C→Fe+CO is promoted and the furnace heat decreases.
ところで、ガス流が多量に炉内周辺部に流れ
ると、Na、K、Pb等の炉内付着物及び停滞層
が剥離し、壁落ちすることにより、その部分の
炉壁温度が急激に上昇する。この急激な温度上
昇を検知すれば炉熱低下が予測できる。 By the way, when a large amount of gas flows into the periphery of the furnace, the deposits and stagnation layer inside the furnace such as Na, K, and Pb peel off and fall off the wall, causing the temperature of the furnace wall in that area to rise rapidly. . If this rapid temperature rise is detected, a decrease in furnace heat can be predicted.
B 第2の炉熱低下理由
また、高炉の炉熱低下の一因として、以下に
示すものが考えられる。B. Second Reason for Decrease in Furnace Heat The following may be considered as a cause of the decrease in furnace heat in the blast furnace.
高炉内の荷下がり速度がAと同様の理由で上
がると、いわゆる生鉱下りにより高炉内の融着
帯レベルが下がり、炉熱低下が起こる。 When the unloading speed in the blast furnace increases for the same reason as A, the level of the cohesive zone in the blast furnace decreases due to so-called raw ore unloading, resulting in a decrease in furnace heat.
ところで、融着帯レベルが下がると、該当部
分における炉壁温度も急激に下降する。この急
激な温度下降を検知すれば炉熱低下が予測でき
る。 By the way, when the cohesive zone level decreases, the furnace wall temperature at the corresponding portion also decreases rapidly. If this rapid temperature drop is detected, a decrease in furnace heat can be predicted.
C 第3の炉熱低下理由
さらに、高炉の炉熱低下の一因として、以下
に示すものが考えられる。C. Third Reason for Decrease in Furnace Heat Furthermore, the following may be considered as a contributing factor to the decrease in furnace heat in the blast furnace.
高炉羽口から吹き上げる溶銑温度及び溶銑量
調整のための高温空気(ガス流)は通常、炉内
中央部に吹き込んでいる。ところが、A、Bと
同様の理由により、ガス流の一部が炉内周辺部
に流れる場合がある。この状態が長時間続く
と、高炉の炉壁からのガス流の熱放散が正常操
業時より多くなり、その結果、炉熱低下が起こ
る。 High-temperature air (gas flow) blown up from the blast furnace tuyeres to adjust the temperature and amount of hot metal is usually blown into the center of the furnace. However, for the same reasons as A and B, part of the gas flow may flow to the periphery of the furnace. If this state continues for a long time, heat dissipation from the gas flow from the blast furnace wall becomes greater than during normal operation, resulting in a decrease in furnace heat.
ところで、ガス流の一部が定常的に炉内周辺
部に流れると、炉壁温度が徐々に上昇する。こ
のような比較的長時間のゆるやかな温度上昇を
検知すれば炉熱低下が予測できる。 By the way, when a part of the gas flow steadily flows into the periphery of the furnace, the temperature of the furnace wall gradually increases. If such a gradual temperature rise over a relatively long period of time is detected, a decrease in furnace heat can be predicted.
D 第1〜第3の予測手段
第1図a,bは、各々この発明の一実施例で
用いられる内壁温度計の配置を示す側面断面
図、平面断面図である。内壁温度計3は同図a
に示すように、高炉1の高さ方向に7個(背部
3個、腹部2個、朝顔部2個)、同図bに示す
ように高炉1の周方向に4個設置する。つま
り、4方向7レベルで計28個の内壁温度計3を
設置する訳である。D First to Third Prediction Means FIGS. 1a and 1b are a side sectional view and a plan sectional view, respectively, showing the arrangement of an inner wall thermometer used in an embodiment of the present invention. Inner wall thermometer 3 is shown in Figure a.
As shown in Figure b, seven pieces are installed in the height direction of the blast furnace 1 (three on the back, two on the abdomen, and two on the morning glory part), and four pieces are installed in the circumferential direction of the blast furnace 1, as shown in Figure b. In other words, a total of 28 inner wall thermometers 3 are installed at 7 levels in 4 directions.
内壁温度計は例えば、本出願人による実開昭
57−81531、実公昭59−16816に開示されたもの
を用いてもよく、第2図は後者に開示された内
壁温度計(以下これを「FMセンサ」という。)
を示す概念図である。 For example, the inner wall thermometer is a
57-81531 and Utility Model Publication No. 59-16816 may be used, and Figure 2 shows the inner wall thermometer disclosed in the latter (hereinafter referred to as "FM sensor").
FIG.
同図において、4は2本の導線5が絶縁的に
平行して埋設され前方端側に感温部6を有する
シース型測温体であり、シース型測温体4は複
数本を、夫々の感温部6が長さ方向の異なる部
位に配置される様に平行配列されており、さら
にシース型ダミー棒7を感温部6の先端に接続
して、最先端を揃えている。シース型ダミー棒
7は2本の導線5が絶縁的に平行して埋設さ
れ、シース型測温体4と実質的に同一の熱伝導
性を有する。FMセンサ3はこのシース型測温
体4を絶縁材8で相互に非接触に保ち、シース
管9内に収納することにより形成される。 In the figure, reference numeral 4 denotes a sheath type thermometer having two conductive wires 5 buried in parallel insulatively and having a temperature sensing part 6 on the front end side. The temperature-sensing parts 6 are arranged in parallel so that they are arranged at different lengthwise positions, and a sheath-type dummy rod 7 is connected to the tips of the temperature-sensing parts 6 to align the leading ends. The sheath type dummy rod 7 has two conductive wires 5 buried in parallel insulating manner, and has substantially the same thermal conductivity as the sheath type temperature measuring element 4. The FM sensor 3 is formed by storing the sheath type temperature sensing element 4 in a sheath tube 9 while maintaining them in a non-contact manner with an insulating material 8.
第3図はFMセンサ3の設置説明図である。
同図において、10〜13は高炉の炉壁であ
り、10はレンガ、11はステーブ、12はス
タンプ、13は鉄皮である。FMセンサ3は同
図に示すように、パツキン14及び溶接部15
への溶接により、炉壁内部に設置されている。
なお、16は充填材であり、17はミクロ注入
口であり充填材16を注ぎ込む箇所である。 FIG. 3 is an explanatory diagram of the installation of the FM sensor 3.
In the figure, 10 to 13 are the walls of the blast furnace, 10 is a brick, 11 is a stave, 12 is a stamp, and 13 is an iron shell. As shown in the figure, the FM sensor 3 is connected to a packing 14 and a welded part 15.
It is installed inside the furnace wall by welding to.
Note that 16 is a filler, and 17 is a micro injection port into which the filler 16 is poured.
なお、ここで説明したFMセンサ3はその設
置及び構造上、炉壁の侵食と共にFMセンサ3
自体も侵食され、シース型測温体4が炉壁近傍
の炉内に露出する場合もあり、実際には「炉壁
温度」と共に「炉壁近傍の炉内温度」を測定し
ていることになる。以下、両者を含めた概念を
「炉壁温度」として述べる。FMセンサ3は上
述のように従来のシース熱電対等の温度計に比
べ、多数の測定点を有し、迅速な測温応答を満
足し、長期の連続的な温度測定が可能であり、
信頼性の向上、耐久性の向上、施工性の向上等
が計られている。 Furthermore, due to the installation and structure of the FM sensor 3 described here, the FM sensor 3 may deteriorate due to erosion of the furnace wall.
The sheath type temperature measuring element 4 may also be eroded and exposed inside the furnace near the furnace wall, meaning that it is actually measuring the ``furnace temperature near the furnace wall'' as well as the ``furnace wall temperature.'' Become. Hereinafter, a concept including both will be described as "furnace wall temperature." As mentioned above, the FM sensor 3 has a large number of measurement points compared to conventional thermometers such as sheathed thermocouples, satisfies rapid temperature measurement response, and is capable of long-term continuous temperature measurement.
Efforts are being made to improve reliability, durability, and workability.
各FMセンサ3は、第4図に示すように所定
サンプリング時間Δtごとに高炉1の内壁温度
を測定している。ここで、時刻jのi番目の
FMセンサ3の内壁温度をTj,iとし、時刻jの
1サンプリング時間Δt前の内壁温度をTj-1,iと
すると、Tj,iとTj-1,iとの内壁温度差(差分値)
ΔTj,iは、
ΔTj,i=Tj,i−Tj-1,i ……(1)
となる。この状態を第5図に示す。 Each FM sensor 3 measures the inner wall temperature of the blast furnace 1 at every predetermined sampling time Δt, as shown in FIG. Here, the i-th at time j
If the inner wall temperature of the FM sensor 3 is T j,i and the inner wall temperature one sampling time Δt before time j is T j-1,i , then the inner wall temperature difference between T j,i and T j-1,i is (difference value)
ΔT j,i becomes ΔT j,i =T j,i −T j-1,i (1). This state is shown in FIG.
この差分値ΔTj,iに、各FMセンサ3毎の高
さ、周方向等を考慮して重みwiを乗ずる。さら
に、差分値ΔTj,iが負のものに対しては、vi=
0、それ以外のものに対しては、vi=1を示す
正負パラメータviも乗じ、時刻jの補正差分値
(正の差分値)CTj,iを得る。 This difference value ΔT j,i is multiplied by a weight w i in consideration of the height, circumferential direction, etc. of each FM sensor 3. Furthermore, for negative difference values ΔT j,i , v i =
0 and other values, it is also multiplied by a positive/negative parameter v i indicating v i =1 to obtain a corrected difference value (positive difference value) CT j,i at time j.
CTj,i=wi・vi・ΔTj,i ……(2)
次に、補正差分値CTj,iの全FMセンサ3に対
する総和をとり、これをST1jとする。 CT j,i = w i · v i · ΔT j,i (2) Next, the sum of the corrected difference values CT j,i for all FM sensors 3 is calculated, and this is set as ST1 j .
ST1j=28
〓i=1
CTj,i ……(3)
そして次(4)式に従い、この差分値総和ST1jの
値が予め定められた閾値ε1より大きくなれば、急
激な温度上昇があつたとみなし所定期間評価点1
を与える。 ST1 j = 28 〓 i=1 CT j,i ...(3) Then, according to the following equation (4), if the value of this sum of difference values ST1 j becomes larger than the predetermined threshold ε 1 , there will be a rapid temperature rise. 1 evaluation point for a specified period of time
give.
ST1j≧ε1 ……(4)
以上がAの理由に基づく第1の予測手段であ
る。 ST1 j ≧ε 1 ...(4) The above is the first prediction means based on reason A.
Bの理由に基づく第2の予測手段は、以下に示
す通りである。 The second prediction means based on reason B is as shown below.
(2)式において、正負パラメータviは差分値ΔTj,i
が正のものに対しては、vi=0、それ以外のもの
に対しては、vi=1とし、次に、補正差分値CTj,i
の絶対値の全FMセンサ3に対する総和をとり、
これをST2jとする。 In equation (2), the positive and negative parameters v i are the difference values ΔT j,i
For positive values, set v i = 0; for other cases, set v i = 1, and then set the corrected difference value CT j,i
Take the sum of the absolute values of for all FM sensors 3,
Let this be ST2 j .
ST2j=28
〓i=1
|CTj,i| ……(3)′
そして次(4)′式に従い、(3)′式に基づく差分値総
和ST2jの値が予め定められた閾値ε2より大きく
なれば、生鉱下りによる急激な温度下降があつた
とみなし、所定期間備価点1を与える。 ST2 j = 28 〓 i=1 | CT j,i | ...(3)′ Then, according to the following equation (4)′, the value of the sum of difference values ST2 j based on equation (3)′ is set to a predetermined threshold ε If it is greater than 2 , it is assumed that there has been a sudden temperature drop due to raw ore descent, and a reserve point of 1 is given for a predetermined period of time.
ST2j≧ε2 ……(4)′
Cの理由に基づく第3の予測手段は、以下に示
す通りである。 ST2 j ≧ε 2 ...(4)' The third prediction means based on the reason C is as shown below.
(2)式の正負パラメータviは第1の予測手段と同
様、差分値ΔTj,iが負のものに対しては、vi=0、
それ以外のものに対してはvi=1とする。また時
刻jのkサンプリング時間前(すなわちΔt×k
時間前)の補正差分値をCTj-k,iとし、この補正差
分値の所定の時間幅nΔtの移動平均の時間jにお
ける値の全FMセンサ3に対する総和をとり、こ
れをST3jとする。 Similar to the first prediction means, the positive/negative parameter v i in equation (2) is v i =0, when the difference value ΔT j,i is negative.
For other cases, v i =1. Also, k sampling times before time j (i.e., Δt×k
CT jk,i is the corrected difference value (time before), and the sum of the moving average of this corrected difference value over a predetermined time width nΔt at time j for all FM sensors 3 is taken, and this is set as ST3 j .
ST3j
=28
〓i=1
{(1/(n+1))o
〓k=0
CTj-k,i} ……(3)″
そして次(4)″式に従い、この移動平均総和ST3j
の値が予め定められた閾値ε3より大きくなれば、
ゆるやかな温度上昇が長期間あつたとみなし、所
定期間評価点1を与える。ST3 j = 28 〓 i=1 {(1/(n+1)) o 〓 k=0 CT jk,i } ...(3)'' Then, according to the following formula (4), this moving average sum ST3 j
If the value of becomes larger than the predetermined threshold ε 3 , then
It is assumed that there has been a gradual temperature rise for a long period of time, and an evaluation point of 1 is given for the specified period.
ST3j≧ε3 ……(4)″
上記した第1〜第3の予測手段は、各々炉壁温
度差(差分値)により行なつているため、炉壁温
度の絶対値の上下によらず、正確な予測を行なう
ことができる。しかも、FMセンサ3はその施工
性の良さ及び測温応答性の良さから、高炉全周を
覆うように配置でき、連続的な内壁温度差が把握
できることで、さらに正確な予測を行なうことが
できる。 ST3 j ≧ε 3 ……(4)″ Since the first to third prediction methods described above are each performed based on the furnace wall temperature difference (difference value), they are independent of the upper and lower absolute values of the furnace wall temperature. , accurate predictions can be made.Furthermore, due to its ease of installation and good temperature measurement response, the FM sensor 3 can be placed to cover the entire circumference of the blast furnace, making it possible to grasp continuous inner wall temperature differences. , more accurate predictions can be made.
また、上記した第1〜第3の予測手段は、コン
ピユータにより実現が可能となる。第6図は第1
の予測手段の処理の流れを示すフローチヤートで
ある。同図において、ステツプS1で各FMセンサ
3の炉壁温度Tj,iをサンプリング時間Δt毎に測定
する。次に、ステツプS2において各FMセンサ3
の差分値を(1)式に基づき計算する。 Further, the first to third prediction means described above can be realized by a computer. Figure 6 is the first
3 is a flowchart showing the processing flow of the prediction means. In the figure, in step S1, the furnace wall temperature T j,i of each FM sensor 3 is measured at every sampling time Δt. Next, in step S2, each FM sensor 3
Calculate the difference value based on equation (1).
そして、ステツプS3において、(2)、(3)式に基
づく正の差分値総和ST1jを求める。さらに、ス
テツプS4において、この正の差分値総和ST1jと
予め定められた閾値ε1との比較を行い、(4)式を満
足すればステツプS5においてガス流の急激な炉
内周辺流化による炉熱低下が起こるであろうとみ
なし、所定期間評価点1を与える。一方、(4)式を
満足しない場合は、異常なしとみなしステツプ
S1に戻り、以下ステツプS1〜ステツプS4を繰り
返すことで炉熱低下評価を行なう。 Then, in step S3, the sum of positive difference values ST1 j is determined based on equations (2) and (3). Furthermore, in step S4, this positive difference value sum ST1 j is compared with a predetermined threshold value ε 1 , and if equation (4) is satisfied, in step S5, the gas flow is rapidly shifted to the periphery of the furnace. It is assumed that a decrease in furnace heat will occur, and an evaluation point of 1 is given for the predetermined period. On the other hand, if equation (4) is not satisfied, it is assumed that there is no abnormality and the step
Returning to S1, the furnace heat reduction evaluation is performed by repeating steps S1 to S4.
第7図は第2の予測手段の処理の流れを示すフ
ローチヤートである。同図において、ステツプ
S11で各FMセンサの炉壁温度Tj,iをサンプリング
時間Δt毎に測定する。次に、ステツプS12におい
て各FMセンサ3の差分値を(1)式に基づき計算す
る。 FIG. 7 is a flowchart showing the processing flow of the second prediction means. In the figure, the steps
In S11, the furnace wall temperature T j,i of each FM sensor is measured at every sampling time Δt. Next, in step S12, the difference value of each FM sensor 3 is calculated based on equation (1).
そして、ステツプS13において、(2)′、(3)′式に
基づく負の差分値の絶対値総和ST2jを求める。
さらに、ステツプS14において、この負の差分値
の絶対値総和ST2jと予め定められた閾値ε2との
比較を行い、(4)′式を満足すればステツプS15に
おいて荷下り速度が上つたことによる炉熱低下が
起こるであろうとみなし、所定期間評価点1を与
える。一方、(4)′式を満足しない場合は、異常な
しとみなしステツプS11に戻り、以下ステツプ
S11〜ステツプS14を繰り返すことで炉熱低下評
価を行なう。 Then, in step S13, the total absolute value ST2j of the negative difference values is determined based on equations (2)' and (3)'.
Furthermore, in step S14, the absolute value sum ST2 j of the negative difference values is compared with a predetermined threshold value ε 2 , and if equation (4)' is satisfied, it is determined that the unloading speed has increased in step S15. It is assumed that a decrease in the furnace heat will occur due to the above, and an evaluation point of 1 will be given for the predetermined period. On the other hand, if formula (4)' is not satisfied, it is assumed that there is no abnormality and the process returns to step S11.
Furnace heat drop evaluation is performed by repeating steps S11 to S14.
第8図は第3の予測手段の処理の流れを示すフ
ローチヤートである。同図において、ステツプS
21で各FMセンサ3の炉壁温度Tj,iをサンプリ
ング時間Δt毎に測定する。次に、ステツプS22に
おいて各FMセンサ3の差分値を(1)式に基づき計
算する。 FIG. 8 is a flowchart showing the processing flow of the third prediction means. In the figure, step S
At 21, the furnace wall temperature T j,i of each FM sensor 3 is measured at every sampling time Δt. Next, in step S22, the difference value of each FM sensor 3 is calculated based on equation (1).
そして、ステツプS23において(2)″、(3)″式に基
づく正の差分値の時間幅nΔtにおける移動平均総
和ST3jを求める。さらに、ステツプS24におい
て、この正の差分値移動平均総和ST3jと予め定
められた閾値ε3との比較を行い、(4)″式を満足す
ればステツプS25において炉体熱放散による炉熱
低下が起こるであろうとみなし、所定期間評価点
1を与える。一方、(4)″式を満足しない場合は、
異常なしとみなしステツプS21に戻り、以下ステ
ツプS21〜ステツプS24を繰り返すことで炉熱低
下評価を行なう。 Then, in step S23, the moving average sum ST3j of the positive difference values in the time width nΔt is determined based on equations (2)'' and (3)''. Furthermore, in step S24, this positive moving average sum ST3 j is compared with a predetermined threshold value ε 3 , and if the formula (4) is satisfied, then in step S25, the furnace heat decreases due to furnace heat dissipation. It is assumed that this will occur, and an evaluation point of 1 is given for the predetermined period.On the other hand, if formula (4) is not satisfied,
It is assumed that there is no abnormality, and the process returns to step S21, whereupon steps S21 to S24 are repeated to evaluate the decrease in furnace heat.
E 第4、第5の予測手段
ガスクロマトグラフイーによる炉頂ガス成分
分析、送風条件、原料装入条件などにより、ソ
ルロスC量(Kg/t−p)をサンプリング時間
Δtごとに算出する。ここで、時刻jにおける
ソルロスC量をxjとし、時刻jよりもkサンプ
リング時間前(すなわちΔt×k時間前)のソ
ルロスC量をxj-kとすると、現在の時刻jにお
ける所定時間幅nΔtの移動平均xMは、
xM=1/n+1o
〓k=0
xj-k ……(5)
で計算できる。E Fourth and Fifth Prediction Means Sol loss C amount (Kg/t-p) is calculated for each sampling time Δt based on top gas component analysis by gas chromatography, ventilation conditions, raw material charging conditions, etc. Here, if the amount of Sol loss C at time j is x j and the amount of Sol loss C at k sampling times before time j (that is, Δt×k hours before) is x jk , then The moving average x M can be calculated as x M =1/n+1 o 〓 k=0 x jk (5).
(5)式に基づくxMをサンプリング時間Δtごとに
計算し、下記(6)式により、xMが予め定めておい
た閾値εxi(i=1〜n)(εx1<εx2…<εxo)を越
え
た時の最大の閾値εxiにより評価点iを与え、評
価を行なう。 x M based on equation (5) is calculated for each sampling time Δt, and x M is set to a predetermined threshold ε xi (i=1 to n) (ε x1 <ε x2 ...< An evaluation point i is given based on the maximum threshold value ε xi when ε xo ) is exceeded, and evaluation is performed.
xM>εxi(i=1〜n) ……(6) 以上が第4の予測手段である。 x M > ε xi (i=1 to n) (6) The above is the fourth prediction means.
また、ガスクロマトグラフイーにより検出され
る炉頂ガス中の窒素量(%)(以下、「ガスクロ
N2量」と言う。)はソルロスC量と強い負の相関
があり、ソルロスC量の増加に代え、ガスクロ
N2量の減少により、高炉炉熱低下が予測できる。 In addition, the amount of nitrogen (%) in the furnace top gas detected by gas chromatography (hereinafter referred to as “gas chromatography”)
2 amounts of N. ) has a strong negative correlation with the amount of Sol loss C, and instead of increasing the amount of Sol loss C, gas chromatography
A decrease in blast furnace heat can be predicted by decreasing the amount of N2 .
このことにより、現在の時刻jにおけるガスク
ロN2量をyjとし、時刻jよりもkサンプリング
時間前(すなわちΔt×k時間前)のガスクロN2
量をyj-kとすると、現在の時刻jにおける所定時
間幅nΔtの移動平均yMは、
yM=1/n+1o
〓k=0
yj-k ……(7)
で計算できる。 As a result, the amount of gas chromatographic N 2 at the current time j is set as y j , and the amount of gas chromatic N 2 at k sampling times before time j (that is, Δt × k hours before)
When the quantity is y jk , the moving average y M of the predetermined time width nΔt at the current time j can be calculated as follows: y M =1/n+1 o 〓 k=0 y jk (7).
(7)式に基づくyMをサンプリング時間Δt毎に計
算し、下記(8)式により、yMが予め定めておいた
閾値εyj(j=1〜m)(εy1>εy2>…εyn)を下回
つた時の最小の閾値εyjにより評価点jを与え、
炉熱低下評価を行なう。 Calculate y M based on equation (7) for each sampling time Δt, and use equation (8) below to calculate y M to a predetermined threshold ε yj (j=1 to m) (ε y1 > ε y2 >... An evaluation point j is given by the minimum threshold ε yj when the value falls below ε yn ),
Evaluate furnace heat drop.
yM<εyj(j=1〜1m) ……(8) 以上が第5の予測手段である。 y M <ε yj (j=1 to 1 m) (8) The above is the fifth prediction means.
さらに、第4の予測手段であるソルロスC量の
移動平均を求めるに際し、ソルロスC量の瞬時値
が第9図aに示すようにノイズ等の原因で異常値
E1、E2を発生する場合がある。ここで、時刻j
のソルロスC量をxj、1サンプリング時間Δt前
のソルロスC量をxj-1とすると、ソルロスC量の
差分値の絶対値Δxjは
Δxj=|xj−xj-1| ……(9)
となる。このΔxjを閾値εzと同図bのように比較
することで異常値E1、E2を見つけだし、同図c
に示すように直前の測定値と置き換えることによ
り平滑化をはかる方法が考えられる。この方法を
適用することにより、より正確なソルロスC量の
移動平均が求まり、その結果、かなり精度の高い
予測が可能となる。 Furthermore, when calculating the moving average of the Sol loss C amount, which is the fourth prediction means, the instantaneous value of the Sol loss C amount is an abnormal value due to noise etc. as shown in Figure 9a.
E1 and E2 may occur. Here, time j
When the Sol loss C amount is x j and the Sol loss C amount one sampling time Δt before is x j-1 , the absolute value of the difference value of the Sol loss C amount Δx j is Δx j = | …(9) becomes. By comparing this Δx j with the threshold value ε z as shown in b in the same figure, abnormal values E1 and E2 are found, and
One possible method is to smooth the data by replacing it with the previous measured value, as shown in the figure below. By applying this method, a more accurate moving average of the Sol loss C amount can be obtained, and as a result, prediction with considerably high accuracy is possible.
このような異常値補正を含んだソルロスC量の
移動平均による炉熱低下予測方法はコンピユータ
を用いて実現することが可能である。第10図は
その処理の流れを示すフローチヤートである。同
図において、まず、ステツプS31において、n段
階にεx1<εx2…εxoの大きさで、閾値x1〜εxoを設定
する。そして、ステツプS32でソルロスC量の瞬
時値xjをサンプリング時間Δt毎に求める。そし
て、ステツプS33においてソルロスC量の差分値
の絶対値Δxjを求め、次にステツプS34において
差分値の絶対値Δxjが閾値εzと比較して大きい場
合、ステツプS35において、この瞬時値xjは異常
値となし、直前の測定値xj-1に置き換え、ステツ
プS36に移行する。一方、ステツプS34において
閾値εzより小さい場合は、瞬時値xjを変更するこ
となく、ステツプS36に移行する。ステツプS36
では時間幅nΔtの移動平均xMを求め、次のステツ
プS37において評価点iを0に初期設定する。 A method of predicting a decrease in furnace heat by using a moving average of the Sol loss C amount, which includes such an abnormal value correction, can be realized using a computer. FIG. 10 is a flowchart showing the flow of the process. In the figure, first, in step S31, threshold values x1 to ε xo are set in n stages with a magnitude of ε x1 <ε x2 . . . ε xo . Then, in step S32, the instantaneous value x j of the amount of Sol loss C is obtained for each sampling time Δt. Then, in step S33, the absolute value Δx j of the difference value of the Sol loss C amount is determined. Next, in step S34, if the absolute value Δx j of the difference value is larger than the threshold value ε z , in step S35, this instantaneous value x j is determined to be an abnormal value and replaced with the previous measured value x j-1 , and the process moves to step S36. On the other hand, if it is smaller than the threshold value ε z in step S34, the process moves to step S36 without changing the instantaneous value x j . Step S36
Then, the moving average x M of the time width nΔt is determined, and the evaluation point i is initially set to 0 in the next step S37.
そして、ステツプS38において、ソルロスC量
移動平均値xMと閾値εx1(i=0より)との比較が
行なわれ、xM≧εx1ならば、ステツプS39において
iの値を0→1と1増し、ステツプS40において
i=nと判定されるか、ステツプS38においてxM
<εx(i+1)と判定されるまで閾値εx(i+1)の値を段階的
に増加させながらステツプS38〜S40を繰返して
評価点iを算出し、ステツプS41に移行する。ま
たxM<εx1ならば、ステツプS39、S40は1度も実
行されず、評価点iは0としてステツプS41に移
行する。最後にステツプS41において、ステツプ
S38〜S40により求められた評価点iを出力する。 Then, in step S38, the moving average value x M of the Sol loss C amount is compared with the threshold value ε x1 (from i=0), and if x M ≥ε x1 , the value of i is changed from 0 to 1 in step S39. 1, and it is determined that i=n in step S40, or x M is determined in step S38.
The evaluation point i is calculated by repeating steps S38 to S40 while increasing the threshold value ε x( i +1) step by step until it is determined that <ε x(i+1) , and the process proceeds to step S41. If x M <ε x1 , steps S39 and S40 are never executed, the evaluation point i is set to 0, and the process moves to step S41. Finally, in step S41, the step
The evaluation point i obtained in S38 to S40 is output.
なお、当然のことながら、上記した異常値処理
のコンピユータへの適用は、ガスクロN2量の移
動平均値yMによる炉熱低下予測の場合において
も同様に実現できる。 Note that, as a matter of course, the above-described abnormal value processing can be applied to a computer in the same way in the case of predicting a decrease in furnace heat using the moving average value yM of the amount of gas chromatin N2 .
上述した第4、第5の予測手段はサンプリング
時間Δt毎の移動平均に基づいているため、予測
を早く得ることができ、しかも精度も十分確かな
ものといえる。 Since the fourth and fifth prediction means described above are based on a moving average for each sampling time Δt, predictions can be obtained quickly and the accuracy can be said to be sufficiently reliable.
F 総合予測手段
D、Eで述べた第1〜第5の予測手段の評価
点を用いることで、以下に述べるように総合予
測を行なう。F Comprehensive prediction means Comprehensive prediction is performed as described below by using the evaluation points of the first to fifth prediction means described in D and E.
第11図は、その処理の流れを示すフローチ
ヤートであり、以下、同図を参照しつつ説明す
る。ステツプS51で第1〜第3の予測手段の評
価点f1〜f3を求める。評価点f1〜f3は1度0→
1になれば、後に述べるデータホールド期間中
は、その値を保持する。したがつて、ホールド
時間hr(2時間程度)が時刻t1で設定されたと
すれば、時刻t1〜時刻t1+1hr(すなわちデータ
ホールド期間)内は評価点f1〜f3各々は、一旦
0→1に変化すれば、再び1→0に変化するこ
とはない。このデータホールド期間は、FMセ
ンサ合計値F(=f1+f2+f3)が、F=0→F>
0に変化した時に変動時刻tCを初期化すると共
に設定され、これ以降hrにより規定される時間
内が前述したようにデータホールド期間とな
る。 FIG. 11 is a flowchart showing the flow of the process, and the following description will be made with reference to the same figure. In step S51, evaluation points f 1 to f 3 of the first to third prediction means are determined. Evaluation points f 1 to f 3 are 0 once →
If it becomes 1, that value is held during a data hold period, which will be described later. Therefore , if the hold time h r (approximately 2 hours) is set at time t 1 , then each of the evaluation points f 1 to f 3 will be , once it changes from 0 to 1, it will never change from 1 to 0 again. During this data hold period, the FM sensor total value F (=f 1 + f 2 + f 3 ) is changed from F=0→F>
It is set at the same time as initializing the fluctuation time t C when it changes to 0, and from then on, the time specified by h r becomes the data hold period as described above.
次にステツプS52において、FMセンサ合計
値Fが0か否かを判別することにより、データ
ホールド時間hrが設定されているか否かを識別
し、F=0であれば、データホールド時間hrの
設定は無いので、ステツプS55に移る。一方、
F≠0であれば、データホールド時間hrが既に
設定されているので、ステツプS53で変動時刻
tcとホールド時間hrと比較することで、現在が
データホールド期間中か否かのチエツクを行な
い、tc≦hrならば、データホールド期間中であ
るので、FMセンサ合計値Fを初期化する必要
はないのでステツプS55へ移る。しかしなが
ら、ステツプS53においてtc>hrならば、デー
タホールド期間は終えたと判断し、ステツプ
S54においてFMセンサ合計値Fを“0”に初
期化する。 Next, in step S52, by determining whether the FM sensor total value F is 0 or not, it is determined whether the data hold time h r is set. If F = 0, the data hold time h r Since there is no setting, the process moves to step S55. on the other hand,
If F≠0, the data hold time h r has already been set, so the variable time is set in step S53.
By comparing t c with the hold time h r , it is checked whether the current data hold period is in progress. If t c ≦ h r , then the data hold period is in progress, so the FM sensor total value F is initialized. Since there is no need to convert, the process moves to step S55. However, if t c > hr in step S53, it is determined that the data hold period has ended, and step
In S54, the FM sensor total value F is initialized to "0".
このように、FMセンサ合計値Fはデータホ
ールド期間を考慮しながら求められる。これは
第1〜第3の予測手段が第4、第5の予測手段
に比べ先見性が高い(予測が早く得られる)た
め、将来の同一時点に対する予測結果を総合評
価するためには第1〜第3の予測手段の予測結
果を一定時間ホールドしておく必要があるから
である。 In this way, the FM sensor total value F is determined while taking into account the data hold period. This is because the first to third prediction means have higher foresight (predictions can be obtained faster) than the fourth and fifth prediction means, so in order to comprehensively evaluate the prediction results for the same point in the future, it is necessary to use the first to third prediction means. This is because it is necessary to hold the prediction result of the third prediction means for a certain period of time.
そして、ステツプS55において、第4の予測
手段による評価点iを算出し、さらにステツプ
S56において、第5の予測手段による評価点j
を算出する。次にステツプS57において総合評
価Cを次式に従い求める。 Then, in step S55, the evaluation score i by the fourth prediction means is calculated, and further step S55 is performed.
In S56, the evaluation score j by the fifth prediction means
Calculate. Next, in step S57, the overall evaluation C is determined according to the following formula.
C=wAF+wBi+wCj ……(9)
ここで、wA、wB、wCは第1〜第3、第4、
第5の予測手段の各々に対する重みである。こ
の総合評価CがステツプS58において吟味さ
れ、C=0ならば、全く炉熱低下の傾向なしと
みなし、ステツプS51に戻り、ステツプS51〜
ステツプ57により再び総合評価Cを求める。一
方、C>0ならばステツプS59において、総合
評価Cの値に応じて重要性を変えたアラームを
出力する。以降ステツプS51に戻り、継続して
炉熱低下予測が続けられる。 C=w A F+w B i+w C j ...(9) Here, w A , w B , w C are the first to third, fourth,
It is a weight for each of the fifth prediction means. This comprehensive evaluation C is examined in step S58, and if C=0, it is assumed that there is no tendency for the furnace heat to decrease at all, and the process returns to step S51.
In step 57, the overall evaluation C is determined again. On the other hand, if C>0, in step S59, an alarm whose importance is changed according to the value of the comprehensive evaluation C is output. Thereafter, the process returns to step S51, and the furnace heat drop prediction continues.
第12図、第13図は総合評価Cによる総合
予測の稼動例を示すグラフであり、aは溶銑温
度代表値、bはFMセンサ正の差分値総和(第
1の予測手段)、cはFMセンサ負の差分値絶
対値総和(第2の予測手段)、dはFMセンサ
正の差分値移動平均総和(第3の予測手段)、
eはガスクロN2量移動平均値(第4の予測手
段)、fはソルロスC量移動平均値(第5の予
測手段)、gは総合評価Cの経時変化を示して
いる。aの溶銑温度代表値は、溶銑鍋への出銑
時(1出銑約2時間)の約30分毎の実測温度の
平均値である。総合評価Cは(9)式において、
wA=wB=wC=1、(6)、(8)式においてn=m=
3とし、総合評価Cの変域を0〜9としてい
る。 Figures 12 and 13 are graphs showing operation examples of comprehensive prediction based on comprehensive evaluation C, where a is the representative value of hot metal temperature, b is the sum of positive difference values of the FM sensor (first prediction means), and c is FM d is the sum of absolute values of negative sensor differences (second prediction means); d is the moving average sum of positive difference values of the FM sensor (third prediction means);
e represents the moving average value of the gas chromatography N2 amount (fourth prediction means), f represents the moving average value of the Sol loss C amount (fifth prediction means), and g represents the change over time in the comprehensive evaluation C. The representative hot metal temperature value a is the average value of the temperatures actually measured approximately every 30 minutes during tapping into the hot metal ladle (one tap for approximately 2 hours). The overall evaluation C is given by equation (9),
w A = w B = w C = 1, n = m = in equations (6) and (8)
3, and the range of the comprehensive evaluation C is 0 to 9.
第12図において総合評価Cの最大値は9を
示し、重大アラームを出力する必要性があり、
実際に同図aに示すように管理温度Tcをか
なり下回つた溶銑温度代表値が計測されてい
る。一方、第13図において、総合評価Cの最
大値は3を示しており、この程度の値であれば
軽アラーム出力で充分であり、実際に同図a
に示す程度の炉熱低下しか起こつていない。 In Fig. 12, the maximum value of the comprehensive evaluation C is 9, and it is necessary to output a serious alarm.
In fact, as shown in Figure a, representative hot metal temperatures have been measured that were considerably lower than the control temperature Tc . On the other hand, in Fig. 13, the maximum value of the comprehensive evaluation C is 3, and if this value is enough, a light alarm output is sufficient, and in fact,
The furnace heat decreased only to the extent shown in the figure.
このように総合評価Cの値に応じて、アラー
ムの度合いを変えることで、昇熱アクシヨンの
程度を細かく変化させることができる。その結
果、必要十分な昇熱アクシヨンを選択すること
が実現できるようになり、炉熱低下を確実に防
止できるのは勿論、過度の昇熱アクシヨンによ
る不要な炉熱上昇を招くことがなく、安定かつ
経済的な高炉操業が可能になつた。 In this way, by changing the degree of alarm according to the value of the comprehensive evaluation C, the degree of heat raising action can be finely changed. As a result, it has become possible to select the necessary and sufficient heating action, which not only reliably prevents a drop in furnace heat, but also prevents unnecessary rises in furnace heat due to excessive heating actions, resulting in stable This also made it possible to operate a blast furnace economically.
G 補足
尚、この実施例における第1〜第3の予測手
段では内壁温度計にFMセンサを用いたが、通
常の測温センサ(例えばシース熱電対)でも寿
命の点で問題はあるものの代用可能である。ま
た、ステーブ温度計、レンガ埋め込み温度計を
用いてもその信頼性、測温応答性の低さから予
測精度は若干低下するものの、代用可能であ
る。G Supplement: In the first to third prediction means in this example, an FM sensor was used as the inner wall thermometer, but a normal temperature sensor (for example, a sheathed thermocouple) can be used as a substitute, although there are problems with the service life. It is. Furthermore, a stave thermometer or a brick-embedded thermometer can be used as a substitute, although the prediction accuracy will be slightly lower due to their reliability and low temperature response.
また、この実施例における第1〜第3の予測
手段では、FMセンサ3を7レベル4方向に28
個設置したが、高炉の特性により適当に設置す
れば良いのは勿論である。 In addition, in the first to third prediction means in this embodiment, the FM sensor 3 is arranged in 7 levels and 28 directions in 4 directions.
Of course, it is sufficient to install them appropriately depending on the characteristics of the blast furnace.
さらに、総合予測はFで述べたように第1〜
第5の予測手段全てを用いるのが望ましいが、
最低限第1〜第3の予測手段の少なくとも1
つ、第4、第5の予測手段の少なくとも1つを
用いることで、Fで述べた例とほぼ同様の効果
が期待できる。またデータホールド期間の設定
も第1〜第3の予測手段で各々独立に行なうこ
ともでき、閾値を第4、第5の予測手段と同様
に複数設ける等の変形例も考えられる。 Furthermore, as mentioned in F, the comprehensive prediction is
Although it is desirable to use all of the fifth prediction means,
At least one of the first to third prediction means
By using at least one of the first, fourth, and fifth prediction means, substantially the same effect as in the example described in F can be expected. Furthermore, the data hold period can also be set independently by the first to third prediction means, and modifications such as providing a plurality of threshold values as in the fourth and fifth prediction means are also conceivable.
(発明の効果)
以上説明したように、この発明によれば、予測
が早く得られ、しかも第1〜第5の予測手段の総
合判断によるため溶銑温度の低下をより正確に予
測し、必要に応じた昇熱アクシヨンをとることが
できる。(Effects of the Invention) As explained above, according to the present invention, a prediction can be obtained quickly, and since it is based on a comprehensive judgment of the first to fifth prediction means, a decrease in hot metal temperature can be predicted more accurately, and when necessary. The corresponding heating action can be taken.
第1図a,bは各々、この発明の一実施例に用
いられるFMセンサの高炉炉壁内の配置を示す側
面断面図、平面断面図、第2図、第3図は各々
FMセンサの概念図、設置説明図、第4図はFM
センサによる測定炉壁温度の経時変化を示すグラ
フ、第5図はFMセンサによる測定炉壁温度の差
分値の経時変化を示すグラフ、第6図は第1の予
測手段の処理の流れを示すフローチヤート、第7
図は第2の予測手段の処理の流れを示すフローチ
ヤート、第8図は第3の予測手段の処理の流れを
示すフローチヤート、第9図a,b,cは各々異
常値を含んだソルロスC量の瞬時値、ソルロスC
量の差分値の絶対値、異常値を取り除いたソルロ
スC量の瞬時値を示すグラフ、第10図は第4の
予測手段の処理の流れを示すフローチヤート、第
11図は総合予測手段の処理の流れを示すフロー
チヤート、第12図、第13図a,b,c,d,
e,f,gは各々、溶銑温度代表値、FMセンサ
正の差分値総和、負の差分値絶対値総和、正の差
分値移動平均総和、ガスクロN2量の移動平均値、
ソルロスC量の移動平均値、総合評価Cの経時変
化を示すグラフ、第14図は従来技術における炉
腹ゾンデの高炉内の配置を示す側面断面図、第1
5図は溶銑温度と炉腹部周辺部温度の相関を示す
グラフ、第16図はソルロスC量の1時間平均値
と溶銑温度の経時変化を時間的に対応させて示し
たグラフ、第17図は昇熱アクシヨンを起こした
時のソルロス量の1時間平均値と溶銑温度の経時
変化を時間的に対応させて示したグラフである。
1……高炉、3……FMセンサ。
1a and 1b are a side sectional view and a plan sectional view showing the arrangement of an FM sensor used in an embodiment of the present invention in the blast furnace wall, respectively; FIGS. 2 and 3 are sectional views, respectively.
Conceptual diagram of FM sensor, installation explanation diagram, Figure 4 is FM
A graph showing changes over time in the furnace wall temperature measured by the sensor, Figure 5 is a graph showing changes over time in the difference value of the furnace wall temperature measured by the FM sensor, and Figure 6 is a flowchart showing the processing flow of the first prediction means. Chaat, 7th
The figure is a flowchart showing the process flow of the second prediction means, FIG. 8 is a flowchart showing the process flow of the third prediction means, and FIGS. 9a, b, and c are Sol losses including abnormal values. Instantaneous value of C amount, Sol loss C
A graph showing the absolute value of the difference value of the quantity and the instantaneous value of the Sol loss C quantity after removing abnormal values. Fig. 10 is a flowchart showing the processing flow of the fourth prediction means. Fig. 11 is the processing of the comprehensive prediction means. Flowcharts showing the flow of Fig. 12, Fig. 13 a, b, c, d,
e, f, and g are the representative value of hot metal temperature, the total sum of positive difference values of the FM sensor, the sum of absolute negative difference values, the moving average sum of positive difference values, the moving average value of gas chromatin N2 amount,
Figure 14 is a graph showing the moving average value of the Sol loss C amount and changes over time in the comprehensive evaluation C.
Figure 5 is a graph showing the correlation between hot metal temperature and temperature around the furnace abdomen, Figure 16 is a graph showing the hourly average value of Sol loss C amount and temporal change in hot metal temperature, and Figure 17 is a graph showing the correlation over time. It is a graph showing the temporal correspondence between the one-hour average value of the Sol loss amount when the heating action occurs and the change over time in the hot metal temperature. 1... Blast furnace, 3... FM sensor.
Claims (1)
度差を測定し、 ある時刻における前記内壁温度差の正の値を示
す部分の合計値が閾値を越えた時に所定期間評価
点を与える第1の予測手段と、 ある時刻における前記内壁温度差の負の値を示
す部分の合計値が閾値を越えた時に所定期間評価
点を与える第2の予測手段と、 ある時刻における前記内壁温度差の正の値を示
す部分の所定時間幅の移動平均値の総和が閾値を
越えた時に所定期間評価点を与える第3の予測手
段とのうちの少なくとも1つ備え、 かつソリユーシヨンロスカーボン量を所定時間
間隔ごとに求め、この求めた値の所定時間幅にお
ける移動平均値が閾値を越えた時に評価点を与え
る第4の予測手段と、 炉頂ガス成分中の窒素量を所定時間間隔ごとに
求め、この求めた値の所定時間幅における移動平
均値が、閾値を下回つた時に評価点を与える第5
の予測手段とのうちの少なくとも1つをさらに備
え、 前記第1〜第3の予測手段のうちの少なくとも
1つの評価点と前記第4、第5の予測手段のうち
の少なくとも1つの評価点の総合評価に従い高炉
炉熱低下予測を行なう高炉炉熱低下予測方法。[Claims] 1. An inner wall thermometer is installed at a predetermined location in the blast furnace, and the inner wall temperature difference is measured at predetermined time intervals with the inner wall thermometer, and the positive value of the inner wall temperature difference at a certain time is determined. a first prediction means for giving an evaluation point for a predetermined period when the total value of the indicated portion exceeds a threshold; and a first prediction means for giving an evaluation point for a predetermined period when the total value of the portion indicating a negative value of the inner wall temperature difference at a certain time exceeds the threshold. a second prediction means for giving a point; and a third prediction means for giving an evaluation point for a predetermined period when the sum of moving average values for a predetermined time width of the portion showing a positive value of the inner wall temperature difference at a certain time exceeds a threshold value. and a fourth means for determining the amount of solution loss carbon at each predetermined time interval and giving an evaluation point when the moving average value of the determined value over a predetermined time width exceeds a threshold value. a prediction means;
further comprising at least one of the prediction means, and the evaluation score of at least one of the first to third prediction means and the evaluation score of at least one of the fourth and fifth prediction means. A blast furnace heat drop prediction method that predicts a blast furnace heat drop based on comprehensive evaluation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP4215387A JPS63210211A (en) | 1987-02-24 | 1987-02-24 | Prediction for lowering of furnace heat in blast furnace |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP4215387A JPS63210211A (en) | 1987-02-24 | 1987-02-24 | Prediction for lowering of furnace heat in blast furnace |
Publications (2)
Publication Number | Publication Date |
---|---|
JPS63210211A JPS63210211A (en) | 1988-08-31 |
JPH039162B2 true JPH039162B2 (en) | 1991-02-07 |
Family
ID=12627997
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP4215387A Granted JPS63210211A (en) | 1987-02-24 | 1987-02-24 | Prediction for lowering of furnace heat in blast furnace |
Country Status (1)
Country | Link |
---|---|
JP (1) | JPS63210211A (en) |
-
1987
- 1987-02-24 JP JP4215387A patent/JPS63210211A/en active Granted
Also Published As
Publication number | Publication date |
---|---|
JPS63210211A (en) | 1988-08-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JPH039162B2 (en) | ||
JPH0692609B2 (en) | Blast furnace furnace heat drop prediction method | |
JPH0587561B2 (en) | ||
JPH0635606B2 (en) | Blast furnace furnace heat drop prediction method | |
JPH0635607B2 (en) | Blast furnace furnace heat drop prediction method | |
JPS63210210A (en) | Prediction for lowering of furnace heat in blast furnace | |
JPS63210219A (en) | Prediction of furnace heat in blast furnace | |
JP3598824B2 (en) | Blast furnace operation method | |
JPS63206414A (en) | Method for predicting reduction of furnace heat of blast furnace | |
JPS63210221A (en) | Prediction of furnace heat in blast furnace | |
JPS63210215A (en) | Prediction of furnace heat in blast furnace | |
JPS63297517A (en) | Method for predicting lowering of furnace heat in blast furnace | |
JPS63210220A (en) | Prediction of furnace heat in blast furnace | |
JPS63210213A (en) | Prediction of furnace heat in blast furnace | |
JPS63206412A (en) | Method for predicting reduction of blast furnace heat | |
JPS63297516A (en) | Method for predicting lowering of furnace temperature in blast furnace | |
JPS63210217A (en) | Prediction of furnace heat in blast furnace | |
JPS63312910A (en) | Method for predicting lowering of furnace heat in blast furnace | |
JPH02294419A (en) | Method for controlling furnace heat in blast furnace | |
JPS596888B2 (en) | Method for detecting deposits on blast furnace walls | |
JPS63210214A (en) | Prediction of furnace heat in blast furnace | |
RU2044058C1 (en) | Method for control of erosion of blast-furnace well | |
JPS63312911A (en) | Method for predicting lowering of furnace heat in blast furnace | |
JPS63206413A (en) | Method for predicting reduction of blast furnace heat | |
JP3990560B2 (en) | Blast furnace bottom potential difference measuring apparatus and soot level evaluation method in blast furnace |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
LAPS | Cancellation because of no payment of annual fees |