特許 6394269 高炉状態解析装置、高炉状態解析方法、およびプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6394269

(24)【登録日】2018年9月7日

(45)【発行日】2018年9月26日

(54)【発明の名称】高炉状態解析装置、高炉状態解析方法、およびプログラム

(51)【国際特許分類】

   C21B 5/00 20060101AFI20180913BHJP

【ＦＩ】

   C21B5/00 310

   C21B5/00 323

【請求項の数】10

【全頁数】32

(21)【出願番号】特願2014-211911(P2014-211911)

(22)【出願日】2014年10月16日

(65)【公開番号】特開2015-206107(P2015-206107A)

(43)【公開日】2015年11月19日

【審査請求日】2017年6月5日

(31)【優先権主張番号】特願2014-79426(P2014-79426)

(32)【優先日】2014年4月8日

(33)【優先権主張国】JP

(73)【特許権者】

【識別番号】000006655

【氏名又は名称】新日鐵住金株式会社

(74)【代理人】

【識別番号】100090273

【弁理士】

【氏名又は名称】國分 孝悦

(72)【発明者】

【氏名】西岡 浩樹

【審査官】 池ノ谷 秀行

(56)【参考文献】

【文献】 特開昭６０−２６２９１０（ＪＰ，Ａ）

【文献】 特開平０５−１２５４１４（ＪＰ，Ａ）

【文献】 特開２００２−１４６４１５（ＪＰ，Ａ）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｃ２１Ｂ ５／００

(57)【特許請求の範囲】

【請求項1】

高炉が有する複数の出銑口のうち、出銑に利用される複数の出銑口に対する出銑順の繰り返しパターンと、前記出銑順で隣り合う２つの出銑口のうち先行する出銑口と後行する出銑口との双方が１回の出銑において同時期に開口している時間であるラップ時間と、前記出銑に利用される複数の出銑口における開口径の初期値である初期出銑口径と、を設定する設定手段と、
前記出銑順の繰り返しパターンと前記ラップ時間とに基づき、次に開口する前記出銑口を開口するタイミングであるか否かを判定する開口タイミング判定手段と、
前記高炉の炉下部に単位時間当たりに滴下する溶銑および溶融スラグの量と、開口している出銑口から排出される溶銑および溶融スラグの排出速度とを境界条件として用いて、前記高炉の炉下部における溶銑および溶融スラグの速度分布を数値解析により導出する速度分布導出手段と、
前記高炉の炉下部に単位時間当たりに滴下する溶銑および溶融スラグの量から求められる溶銑および溶融スラグの体積消滅速度、並びに、前記速度分布導出手段により導出された前記高炉の炉下部における溶銑および溶融スラグの速度分布に基づいて、前記高炉の炉下部における溶融スラグと気体との界面の高さ位置であるスラグレベルを数値解析により導出するスラグレベル導出手段と、
前記スラグレベル導出手段により導出されたスラグレベルが、開口している前記出銑口の炉内側の所定の高さ位置になると、当該出銑口を閉塞するタイミングであると判定する閉塞タイミング判定手段と、を有し、
前記速度分布導出手段は、前記開口タイミング判定手段により判定されたタイミングで、前記次に開口する出銑口が前記初期出銑口径で開口され、その後、当該出銑口が時間の経過とともに拡大されるものとし、且つ、前記閉塞タイミング判定手段により判定されたタイミングで、前記開口している出銑口が閉塞されるものとして、前記高炉の炉下部における溶銑および溶融スラグの速度分布を導出することを特徴とする高炉状態解析装置。

【請求項2】

前記スラグレベル導出手段により導出された前記スラグレベルの状態が、安定的な操業状態を維持するための所定の条件を満たすか否かを判定する収束判定手段を更に有することを特徴とする請求項１に記載の高炉状態解析装置。

【請求項3】

前記出銑順の繰り返しパターンの一周期よりも長い所定の期間において前記スラグレベル導出手段により導出された前記スラグレベルの最大値である最大スラグレベルを導出する最大スラグレベル導出手段と、
前記出銑順の繰り返しパターンの一周期よりも長い所定の期間において前記スラグレベル導出手段により導出された前記スラグレベルの最大値と最小値との差であるスラグレベルの変動幅とを導出するスラグレベル変動幅導出手段と、を更に有し、
前記収束判定手段は、
前記最大スラグレベル導出手段により導出された最大スラグレベルが、安定的な操業状態を維持するための所定の上限値以下であるか否かを判定する最大スラグレベル判定手段と、
前記スラグレベル変動幅導出手段により導出されたスラグレベルの変動幅が、安定的な操業状態を維持するための所定の上限値以下であるか否かを判定するスラグレベル変動幅判定手段と、を更に有し、
前記最大スラグレベル判定手段による判定の結果と、前記スラグレベル変動幅判定手段による判定の結果とに基づいて、前記スラグレベル導出手段により導出された前記スラグレベルの状態が、安定的な操業状態を維持するための所定の条件を満たすか否かを判定することを特徴とする請求項２に記載の高炉状態解析装置。

【請求項4】

前記収束判定手段により、前記スラグレベル導出手段により導出された前記スラグレベルの状態が安定的な操業状態を維持するための所定の条件を満たさないと判定された場合、前記スラグレベル導出手段により導出される前記スラグレベルの状態が前記安定的な操業状態を維持するための所定の条件を満たす方向に変更されるように、前記ラップ時間と前記初期出銑口径の少なくとも何れか一方を変更する出銑方法変更手段を更に有し、
前記開口タイミング判定手段は、前記出銑方法変更手段によりラップ時間が変更された場合には、前記出銑順の繰り返しパターンと前記変更された後のラップ時間とに基づき、次に開口する前記出銑口を開口するタイミングであるか否かを判定し、
前記速度分布導出手段は、前記出銑方法変更手段により初期出銑口径が変更された場合には、前記開口タイミング判定手段により判定されたタイミングで、前記次に開口する出銑口が、前記変更された後の初期出銑口径で開口されるものとして、前記高炉の炉下部における溶銑および溶融スラグの速度分布を導出することを特徴とする請求項２または３に記載の高炉状態解析装置。

【請求項5】

前記速度分布導出手段と前記スラグレベル導出手段による数値解析の結果に基づいて、所定の期間における出銑回数を導出する出銑回数導出手段を更に有し、
前記設定手段は、前記所定の期間における出銑回数の上限値を更に設定し、
前記出銑方法変更手段は、前記出銑回数導出手段により導出された所定の期間における出銑回数が上限値を上回る場合には、前記出銑回数導出手段により導出される所定の期間における出銑回数が小さくなるように、前記ラップ時間と前記初期出銑口径の少なくとも何れか一方を変更することを特徴とする請求項４に記載の高炉状態解析装置。

【請求項6】

前記出銑方法変更手段は、オペレータによる入力操作に基づいて、前記ラップ時間と前記初期出銑口径の少なくとも何れか一方の変更値を入力し、入力した変更値を設定することを特徴とする請求項４または５に記載の高炉状態解析装置。

【請求項7】

前記速度分布導出手段と前記スラグレベル導出手段による数値解析の結果に基づいて、所定の期間における出銑回数を導出する出銑回数導出手段を更に有し、
前記出銑方法変更手段は、前記スラグレベル導出手段により導出されたスラグレベルと、前記出銑回数導出手段により導出された出銑回数と、前記設定手段により設定されたラップ時間と、に基づいて、前記ラップ時間と前記初期出銑口径の少なくとも何れか一方の変更値を自動的に設定することを特徴とする請求項４または５に記載の高炉状態解析装置。

【請求項8】

前記スラグレベル導出手段は、ＶＯＦ（Volume Of Fluid）法を用いて前記スラグレベルを導出することを特徴とする請求項１〜７の何れか１項に記載の高炉状態解析装置。

【請求項9】

高炉が有する複数の出銑口のうち、出銑に利用される複数の出銑口に対する出銑順の繰り返しパターンと、前記出銑順で隣り合う２つの出銑口のうち先行する出銑口と後行する出銑口との双方が１回の出銑において同時期に開口している時間であるラップ時間と、前記出銑に利用される複数の出銑口における開口径の初期値である初期出銑口径と、を設定する設定工程と、
前記出銑順の繰り返しパターンと前記ラップ時間とに基づき、次に開口する前記出銑口を開口するタイミングであるか否かを判定する開口タイミング判定工程と、
前記高炉の炉下部に単位時間当たりに滴下する溶銑および溶融スラグの量と、開口している出銑口から排出される溶銑および溶融スラグの排出速度とを境界条件として用いて、前記高炉の炉下部における溶銑および溶融スラグの速度分布を数値解析により導出する速度分布導出工程と、
前記高炉の炉下部に単位時間当たりに滴下する溶銑および溶融スラグの量から求められる溶銑および溶融スラグの体積消滅速度、並びに、前記速度分布導出工程により導出された前記高炉の炉下部における溶銑および溶融スラグの速度分布に基づいて、前記高炉の炉下部における溶融スラグと気体との界面の高さ位置であるスラグレベルを数値解析により導出するスラグレベル導出工程と、
前記スラグレベル導出工程により導出されたスラグレベルが、開口している前記出銑口の炉内側の所定の高さ位置になると、当該出銑口を閉塞するタイミングであると判定する閉塞タイミング判定工程と、を有し、
前記速度分布導出工程は、前記開口タイミング判定工程により判定されたタイミングで、前記次に開口する出銑口が前記初期出銑口径で開口され、その後、当該出銑口が時間の経過とともに拡大されるものとし、且つ、前記閉塞タイミング判定工程により判定されたタイミングで、前記開口している出銑口が閉塞されるものとして、前記高炉の炉下部における溶銑および溶融スラグの速度分布を導出することを特徴とする高炉状態解析方法。

【請求項10】

請求項１〜８の何れか１項に記載の高炉状態解析装置の各手段としてコンピュータを機能させることを特徴とするプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、高炉状態解析装置、高炉状態解析方法、およびプログラムに関し、特に、高炉の炉下部の状態を解析するために用いて好適なものである。

【背景技術】

【0002】

高炉内で生成される銑鉄を、高炉内に過剰に滞留させずに出銑することは、高炉の操業にとって重要である。したがって、炉内で生成される銑鉄と、出銑口から排出される銑鉄とを管理しなければならない。炉内のスラグレベル（溶融スラグと気体の界面の高さ位置）が所定のレベル以上に上昇すると、送風圧力の変動が大きくなり、安定的な操業状態が維持できなくなる。特に、炉内のスラグレベルが送風羽口のレベルまで到達すると、スラグが送風羽口を閉塞してしまい操業不能となる。また、炉内のスラグレベルの変動幅が大きい場合には、炉内のスラグレベルが送風羽口のレベルに到達していなくても、送風圧力が大きく変動するため、やはり安定的な操業状態の維持が困難になる。

【0003】

従来、炉内で生成される銑鉄と、出銑口から排出される銑鉄とを管理する方法として以下の方法が採られていた。すなわち、炉内で生成される銑鉄の量を、高炉に装入する鉄鉱石の量から導出する。また、出銑口から排出される銑鉄・スラグの量を、それぞれ、受銑容器（トーピードカー）・水砕のコンベアの秤量を行うことにより導出する。
しかしながら、このような方法では、秤量のタイミングのずれや、秤量誤差等により、正確な管理が困難であった。

【0004】

近年、高炉操業は、低コークス比操業へ移行しており、炉頂から炉内に装入する鉄鉱石の量がコークスの量に対して増加している。そのため、コークス粒子1個あたりの反応負荷の増大に起因する、炉下部におけるコークス粒径の低下や、還元反応により劣化したコークスから発生する粉の影響による、炉下部の通液性の悪化が懸念されている。加えて、出銑比も増加傾向にあり、溶銑やスラグが十分に排出されずに炉内に残留することにより、炉内のスラグレベルの上昇と、炉内のスラグレベルの変動の増大が生じ、安定的な高炉操業の妨げとなる。

【0005】

高炉の出銑比が増大すると、炉内に貯留する溶銑と溶融スラグを速やかに排出することができず、炉内に残留する溶銑と溶融スラグが増えると安定な高炉操業ができない。
これに対して、近年の高炉の大型化に伴う出銑量の増加に伴い、複数の出銑口から同時に出銑するラップ出銑の形態がとられることがある。炉内での溶銑および溶融スラグの生成速度が早くなってきており、出銑滓（高炉の出銑口から流出する溶銑と溶融スラグの混合物）を定量的に且つ正確に管理するシステムが望まれる。

【0006】

特許文献１には、高炉の炉底の外壁鉄皮に複数個の検出接点を取り付け、底盤レベルの検出接点と、それよりも高い位置の検出接点と、の電位差を測定することが開示されている。特許文献１では、この電位差の測定値が予め定めた設定値を超えた時に、溶銑と溶融スラグの生成速度を減少させるアクションと、溶銑と溶融スラグの排出速度を増加させるアクションとのいずれか一方または両方を実施する。

【0007】

また、非特許文献１には、出銑口内圧力損失予測モデルと炉床内流動解析モデルとで構成される高炉炉下部数学モデルにより、炉下部における界面（溶融スラグと気体の界面、および溶銑と溶融スラグの界面）の形状と、炉内に残っている溶銑・溶融スラグの量を推定することが開示されている。

【先行技術文献】

【特許文献】

【0008】

【特許文献1】特開２００３−３４８０９号公報

【非特許文献】

【0009】

【非特許文献1】西岡浩樹、外２名、「高炉炉床部における液流れの３次元数学モデル」、鉄と鋼、一般社団法人日本鉄鋼協会、２００６年、Vol.９２、No.１２、p.９６７−９７５

【非特許文献2】Kouji Takatani、外２名、「Three-dimensional Dynamic Simulator for Blast Furnace」、ISIJ International、一般社団法人日本鉄鋼協会、１９９９年、Vol.３９、No.１、p.１５−２２

【発明の概要】

【発明が解決しようとする課題】

【0010】

しかしながら、特許文献１に記載の発明では、測定される電位差が炉内の銑滓レベル（炉下部における溶銑・溶融スラグの表面の高さ位置）に応じて変化するメカニズムの詳細が不明である。炉内の耐火物の損耗の程度により、検知接点で検知する電位が変化すると考えられる。したがって、この電位の変化が、推定される炉内の銑滓レベルの推定精度に影響すると考えられる。また、特許文献１には、銑滓レベルが上昇した際の対応、すなわち、銑滓レベルを低減するための具体的な操業については記述されていない。

【0011】

一方、非特許文献１に記載の発明では、高炉炉下部数学モデルを用いることにより、炉内の銑滓レベルの推定精度を向上させることができる。しかしながら、非特許文献１に記載の発明では、同時に使用する出銑口の数が１つである場合を前提としている。したがって、複数の出銑口を同時に使用する出銑形態に対応することが容易ではない。

【0012】

そこで、本発明は、複数の出銑口から同時に出銑する出銑形態であっても、高炉の炉下部における溶銑・溶融スラグの状態を解析し、その解析の結果を出銑方法の画定に供することができるようにすることを目的とする。

【課題を解決するための手段】

【0013】

本発明の高炉状態解析装置は、高炉が有する複数の出銑口のうち、出銑に利用される複数の出銑口に対する出銑順の繰り返しパターンと、前記出銑順で隣り合う２つの出銑口のうち先行する出銑口と後行する出銑口との双方が１回の出銑において同時期に開口している時間であるラップ時間と、前記出銑に利用される複数の出銑口における開口径の初期値である初期出銑口径と、を設定する設定手段と、前記出銑順の繰り返しパターンと前記ラップ時間とに基づき、次に開口する前記出銑口を開口するタイミングであるか否かを判定する開口タイミング判定手段と、前記高炉の炉下部に単位時間当たりに滴下する溶銑および溶融スラグの量と、開口している出銑口から排出される溶銑および溶融スラグの排出速度とを境界条件として用いて、前記高炉の炉下部における溶銑および溶融スラグの速度分布を数値解析により導出する速度分布導出手段と、前記高炉の炉下部に単位時間当たりに滴下する溶銑および溶融スラグの量から求められる溶銑および溶融スラグの体積消滅速度、並びに、前記速度分布導出手段により導出された前記高炉の炉下部における溶銑および溶融スラグの速度分布に基づいて、前記高炉の炉下部における溶融スラグと気体との界面の高さ位置であるスラグレベルを数値解析により導出するスラグレベル導出手段と、前記スラグレベル導出手段により導出されたスラグレベルが、開口している前記出銑口の炉内側の所定の高さ位置になると、当該出銑口を閉塞するタイミングであると判定する閉塞タイミング判定手段と、を有し、前記速度分布導出手段は、前記開口タイミング判定手段により判定されたタイミングで、前記次に開口する出銑口が前記初期出銑口径で開口され、その後、当該出銑口が時間の経過とともに拡大されるものとし、且つ、前記閉塞タイミング判定手段により判定されたタイミングで、前記開口している出銑口が閉塞されるものとして、前記高炉の炉下部における溶銑および溶融スラグの速度分布を導出することを特徴とする。

【0014】

本発明の高炉状態解析方法は、高炉が有する複数の出銑口のうち、出銑に利用される複数の出銑口に対する出銑順の繰り返しパターンと、前記出銑順で隣り合う２つの出銑口のうち先行する出銑口と後行する出銑口との双方が１回の出銑において同時期に開口している時間であるラップ時間と、前記出銑に利用される複数の出銑口における開口径の初期値である初期出銑口径と、を設定する設定工程と、前記出銑順の繰り返しパターンと前記ラップ時間とに基づき、次に開口する前記出銑口を開口するタイミングであるか否かを判定する開口タイミング判定工程と、前記高炉の炉下部に単位時間当たりに滴下する溶銑および溶融スラグの量と、開口している出銑口から排出される溶銑および溶融スラグの排出速度とを境界条件として用いて、前記高炉の炉下部における溶銑および溶融スラグの速度分布を数値解析により導出する速度分布導出工程と、前記高炉の炉下部に単位時間当たりに滴下する溶銑および溶融スラグの量から求められる溶銑および溶融スラグの体積消滅速度、並びに、前記速度分布導出工程により導出された前記高炉の炉下部における溶銑および溶融スラグの速度分布に基づいて、前記高炉の炉下部における溶融スラグと気体との界面の高さ位置であるスラグレベルを数値解析により導出するスラグレベル導出工程と、前記スラグレベル導出工程により導出されたスラグレベルが、開口している前記出銑口の炉内側の所定の高さ位置になると、当該出銑口を閉塞するタイミングであると判定する閉塞タイミング判定工程と、を有し、前記速度分布導出工程は、前記開口タイミング判定工程により判定されたタイミングで、前記次に開口する出銑口が前記初期出銑口径で開口され、その後、当該出銑口が時間の経過とともに拡大されるものとし、且つ、前記閉塞タイミング判定工程により判定されたタイミングで、前記開口している出銑口が閉塞されるものとして、前記高炉の炉下部における溶銑および溶融スラグの速度分布を導出することを特徴とする。

【0015】

本発明のプログラムは、前記高炉状態解析装置の各手段としてコンピュータを機能させることを特徴とする。

【発明の効果】

【0016】

本発明によれば、出銑順の繰り返しパターンとラップ時間とに基づき、次に開口する出銑口を開口するタイミングであるか否かを判定すると共に、スラグレベルが開口している出銑口の炉内側の所定の高さ位置になると、当該出銑口を閉塞するタイミングであると判定するようにし、これらの判定の結果を反映して、高炉の炉下部における溶銑および溶融スラグの速度分布の導出とスラグレベルの導出とを行う。したがって、複数の出銑口から同時に出銑する出銑形態であっても、高炉の炉下部における溶銑・溶融スラグの状態を解析し、その解析の結果を出銑方法の画定に供することができる。

【図面の簡単な説明】

【0017】

【図1】高炉の炉下部の概略構成の一例を示す図である。

【図2】高炉に設けられている出銑口の一例を示す図である。

【図3】タッチ出銑、ラップ出銑、およびパラレル出銑の一例を説明する図である。

【図4】タッチ出銑の場合の炉下部の溶銑と溶融スラグの状態の一例を示す図である。

【図5】パラレル出銑の場合の炉下部の溶銑と溶融スラグの状態の一例を示す図である。

【図6】出銑を開始してからのスラグレベルの経時変化の一例を示す図である。

【図7】高炉状態解析装置の機能的な構成の一例を示す図である。

【図8】高炉状態解析装置の全体処理の一例を説明するフローチャートである。

【図9】ステップＳ８０２の処理の詳細の一例を説明するフローチャートである。

【図10】ステップＳ８０３の処理の詳細の一例を説明するフローチャートである。

【図11】ステップＳ８０４の処理の詳細の一例を説明するフローチャートである。

【発明を実施するための形態】

【0018】

以下、図面を参照しながら、本発明の一実施形態を説明する。尚、各図では、説明および表記の都合上、説明に必要な構成のみを必要に応じて簡略化して示す。
（高炉の操業の概要）
図１は、高炉の炉下部の概略構成の一例を示す図である。図１は、高炉の軸に沿って切った断面を示す。
高炉では、炉頂部から、原料である鉄鉱石とコークスを層状に装入すると共に、炉下部の送風羽口１１０から約１２００［℃］の高温の熱風を吹き込み、高温還元反応によって溶銑を製造する。溶銑と副生成物である溶融スラグは炉内を滴下し（図１に示す白抜きの矢印線を参照）、炉底に溜まって湯溜まりを形成する。尚、銑鉄の比重はスラグの比重よりも大きいことから、湯溜まりにおいては、溶銑の上に溶融スラグが位置する。

【0019】

この湯溜まりに向かって炉外から貫通孔（出銑口ＴＨ）を形成できるように、炉底部の出銑口形成位置には、マッド材と称される閉塞材（耐火物）が充填されている。このマッド材の部分を炉外から開口機で機械的に開口することにより出銑口が形成され、溶銑と溶融スラグとの混合液体が出銑口から炉外に排出される（図１の出銑口ＴＨから高炉の外側に向かう矢印線を参照）。出銑口ＴＨの大きさは、時間の経過とともに大きくなる。溶銑と溶融スラグとの混合液体の排出を継続させた後、開口した出銑口をマッド材で閉塞する。以上のようにして、１つの出銑口についての１回の出銑作業が終了する。

【0020】

（高炉の具体的な操業）
表１に、同時に使用する出銑口の最大数が１つである場合の高炉の操業方法の一例を示す。また、表２に、同時に使用する出銑口の最大数が２つである場合の高炉の操業方法の一例を示す。尚、表１および表２は、高炉の操業方法の概念を説明するためのものであり、大凡の表記としている。例えば、溶銑と溶融スラグは、実際には分離した相として炉内に存在するが、表１および表２では、両者が混合した相であるとした場合について表記する。

【0021】

【表1】

【0022】

【表2】

【0023】

図２は、高炉に設けられている出銑口の一例を示す図である。図２では、４つの出銑口１ＴＨ〜４ＴＨが設けられている場合を例に挙げて示す。また、図２は、高炉を、その高さ方向に沿って見た図である（図２に示す円は、高炉の炉壁（側壁）を示す）。
図３は、タッチ出銑、ラップ出銑、およびパラレル出銑の一例を説明する図である。図３に示すように、４つの出銑口がある場合、通常、そのうちの１つは出銑樋の補修のために使用されない。したがって、３つの出銑口が使用される。そこで、本実施形態では、図２に示す３つの出銑口１ＴＨ、２ＴＨ、３ＴＨをこの順で開口・閉塞して出銑作業を行う場合を例に挙げて示す。尚、図３において、黒塗りの横長の矩形で示す間、出銑（溶銑と溶融スラグとの混合液体の排出）が行われることを示す。

【0024】

表１および表２において、出銑口稼働本数は、同時に使用する出銑口の数を表す。銑滓バランスは、単位時間当たりに排出される銑滓（溶銑・溶融スラグ）の量と単位時間当たりに生成（炉下部に滴下）される銑滓の量との大小関係を表す。銑滓レベルは、出銑口を基準とした場合の銑滓相の表面（溶融スラグと気体との界面）の高さ位置を表す（図１を参照）。

【0025】

表２において、ラップ率（Ｒ）とは、或る出銑口における１回の出銑時間に対する、当該出銑口と当該出銑口の次に開口する出銑口の、１回の出銑において重複する出銑時間の割合をいう。ラップ率（Ｒ）の具体例を説明すると、図３（ｂ）および図３（ｃ）において、出銑口１ＴＨにおける１回の出銑時間Ｔ_Aに対する、出銑口１ＴＨ、２ＴＨの１回の出銑において重複する出銑時間Ｔ_Bの割合（＝Ｔ_B÷Ｔ_A）がラップ率（Ｒ）になる。尚、以下の説明では、当該出銑口と当該出銑口の次に開口する出銑口の、１回の出銑において重複する出銑時間（すなわち、１回の出銑においてこれら２つの出銑口が同時期に開口する時間）を必要に応じてラップ時間Ｔ_Bと称する。尚、図３（ａ）では、出銑口１ＴＨ・２ＴＨ、２ＴＨ・３ＴＨ、３ＴＨ・１ＴＨの１回の出銑において重複する出銑時間は０（ゼロ）であるので、ラップ時間Ｔ_Bは０（ゼロ）になる。

【0026】

表１において、「間欠出銑」とは、或る出銑口を開口した後、当該出銑口からガスが噴き出すと当該出銑口を閉塞し、その状態を一定時間保持してから、当該出銑口を再び開口する操業方法である。この場合には、当該出銑口が開口されている期間における単位時間当たりに排出される銑滓（溶銑・溶融スラグ）の量は、単位時間当たりに生成される銑滓（溶銑・溶融スラグ）の量を上回る。すなわち、時々刻々変化する当該出銑口が開口されている期間の銑滓の平均排出速度（＝（当該出銑口からの銑滓の排出量の積算値）÷（当該出銑口からの出銑時間））は、銑滓の平均生成速度を上回る。したがって、閉塞時には、鉄滓層相の表面の高さ位置は上昇する。

【0027】

「連続出銑」とは、或る出銑口を連続的に開口する操業方法である。この場合には、単位時間当たりに排出される銑滓（溶銑・溶融スラグ）の量は、単位時間当たりに生成される銑滓（溶銑・溶融スラグ）の量と等しくなる。ただし、このような操業は、一般には行われない。
「不能」とは、単位時間当たりに生成される銑滓（溶銑・溶融スラグ）の量が、単位時間当たりに排出される銑滓（溶銑・溶融スラグ）の量を上回るために、銑滓相の表面の高さ位置が常に上昇することを表す。この場合、同時に使用する出銑口の数が１つでは、生成した銑滓を排出できないので、同時に使用する出銑口の数を２つにする必要がある（表１の「不能」の欄の「２本出銑へ」を参照）。

【0028】

表２において、「溜め出し出銑」とは、或る出銑口を開口した後、当該出銑口からガスが噴き出すと当該出銑口を閉塞し、その状態を一定時間保持してから、当該出銑口の次の出銑口を開口することを、３つの出銑口１ＴＨ〜３ＴＨについてこの順で繰り返し行う操業方法である。この場合には、当該出銑口が開口されている期間における単位時間当たりに排出される銑滓（溶銑・溶融スラグ）の量は、単位時間当たりに生成される銑滓（溶銑・溶融スラグ）の量を上回る。したがって、閉塞時には、鉄滓層相の表面の高さ位置は上昇する。また、この場合には、同時に使用される出銑口の数は０（ゼロ）または１であるので、ラップ率（Ｒ）は０（ゼロ）になる。

【0029】

「タッチ出銑」とは、図３（ａ）に示すように、或る出銑口を開口した後、当該出銑口からガスが噴き出すと当該出銑口を閉塞し、すぐに、当該出銑口の次の出銑口を開口することを、３つの出銑口１ＴＨ〜３ＴＨについてこの順で繰り返し行う操業方法である。この場合には、単位時間当たりに排出される銑滓（溶銑・溶融スラグ）の量は、単位時間当たりに生成される銑滓（溶銑・溶融スラグ）の量と等しくなる。また、この場合には、同時に使用される出銑口の数は１であるので、前述したように、ラップ率（Ｒ）は０（ゼロ）になる。

【0030】

「ラップ出銑」とは、図３（ｂ）に示すように、或る出銑口を開口した後、当該出銑口からガスが噴き出すタイミングよりも前に、当該出銑口の次の出銑口を開口することを、３つの出銑口１ＴＨ〜３ＴＨについて順番に行う操業方法である。「ラップ出銑」では、ラップ率（Ｒ）が０（ゼロ）を上回り、１／２を下回るようにする。この場合、単位時間当たりに生成される銑滓（溶銑・溶融スラグ）の量は、当該出銑口から単位時間当たりに排出される銑滓の量を上回り、且つ、当該出銑口から単位時間当たりに排出される銑滓（溶銑・溶融スラグ）の量の２倍を下回る。

【0031】

「パラレル出銑」とは、図３（ｃ）に示すように、或る出銑口を開口した後、当該出銑口からガスが噴き出すタイミングよりも前に、当該出銑口の次の出銑口を開口することを、３つの出銑口１ＴＨ〜３ＴＨについてこの順で繰り返し行う操業方法である。この点は「ラップ出銑」と同じである。ただし、「パラレル出銑」では、ラップ率（Ｒ）が１／２になるようにする。したがって、常に２つの出銑口から出銑が行われる（表２の「パラレル出銑」の欄の「出銑口稼働本数」（＝２）を参照）。この場合、単位時間当たりに生成される銑滓（溶銑・溶融スラグ）の量は、当該出銑口から単位時間当たりに排出される銑滓（溶銑・溶融スラグ）の量の丁度２倍になる。

【0032】

「不能」とは、単位時間当たりに生成される銑滓の量が、当該出銑口から単位時間当たりに排出される銑滓の量の２倍を上回るために、銑滓相の表面の高さ位置が常に上昇することを表す。この場合、同時に使用する出銑口の数が２つでは、出銑作業を行えないので、同時に使用する出銑口の数を３つにする必要がある（表１の「不能」の欄の「３本出銑へ」を参照）。ただし、現在の高炉では、このような状態になることは想定され難いので、本実施形態でも、このような状態を想定しない場合を例に挙げて説明する。

【0033】

（経緯）
前述したように、本発明者等は、非特許文献１に記載の高炉炉下部数学モデルを開発した。しかし、非特許文献１に記載の高炉炉下部数学モデルは、１つの出銑口により出銑する形態に限定したものであり、複数の出銑口を同時に使用する出銑形態に対応できない。すなわち、非特許文献１には、或る１つの出銑口を開口した場合の炉下部の溶銑・溶融スラグの状態を解析することが記載されている。しかしながら、２つの出銑口を開口することや、２つの出銑口を開口した場合の炉下部の溶銑・溶融スラグの状態を解析することや、炉下部の溶銑・溶融スラグの状態に基づいて出銑方法を決定することについては記載されていない。

【0034】

表２に示すように、１つの出銑口により出銑する出銑形態には、タッチ出銑（連続出銑）または溜め出し出銑（間欠出銑）があり、出銑比がさほど高くない高炉では、この形態がとられてきた。一方、複数の出銑口から同時に出銑する出銑形態としては、ラップ出銑やパラレル出銑があり、これらの出銑形態は、高炉の高出銑比化に伴い必須の出銑形態となってきた。

【0035】

そこで、本実施形態では、非特許文献１に記載の高炉炉下部数学モデルを、ラップ出銑およびパラレル出銑のような、複数の出銑口から同時に出銑する出銑形態に適用できるように変更し、且つ、炉下部における最大スラグレベルの低減とスラグレベルの変動幅の低減とを達成する出銑方法を導出できるようにする。
図３に示すように、同時に使用する出銑口の数が２つ以上である場合には、ラップ時間を考慮する必要がある。そこで、本実施形態では、出銑口の開口径の初期値である初期出銑口径に加えて、ラップ時間を調整して、スラグレベル推定モデルによる計算を行うことにより、最大スラグレベルとスラグレベルの変動幅が、操業上要求される所定の条件を満たす場合の出銑方法を導出する。

【0036】

ここで、出銑方法は、複数の出銑口のそれぞれを開口するタイミング（すなわち、ラップ時間）と、当該複数の出銑口の初期出銑口径により定められるものとする。
また、最大スラグレベルとは、出銑パターンの一周期よりも長い所定の期間において、炉下部に溜まっている溶融スラグの表面（溶融スラグと気体との界面）の高さ位置が最も高くなる位置をいう。
スラグレベルの変動幅とは、出銑パターンの一周期よりも長い所定の期間において、炉下部に溜まっている溶融スラグの表面（溶融スラグと気体との界面）の高さ位置が最も高くなる位置と、炉下部に溜まっている溶融スラグの表面（溶融スラグと気体との界面）の高さ位置が最も低くなる位置との差をいう。
出銑パターンとは、出銑に利用される出銑口１ＴＨ〜３ＴＨに対する出銑順の繰り返しパターンをいい、本実施形態の例では、出銑口１ＴＨ→出銑口２ＴＨ→出銑口３ＴＨ→出銑口１ＴＨ→出銑口２ＴＨ→・・・のパターンをいう。また、出銑パターンの一周期とは、かかる出銑パターンの繰り返し周期をいう（図３に示す期間Ｔ_pを参照）。

【0037】

（スラグレベル推定モデル）
本実施形態のスラグレベル推定モデルについて説明する。尚、非特許文献１に記載されている内容については、その詳細な説明を省略する。また、ここでは説明を簡単にするために、温度に関わる部分の説明を省略する。ただし、各時間ステップにおいて、炉下部の溶銑・溶融スラグの温度分布を求め、温度に依存する変数（粘度等）を、当該温度分布に基づいて計算してもよい。

【0038】

まず、非特許文献１に示されているように、以下の（１）式に示す連続の式と（２）式〜（４）式に示す運動量保存式（ナビエ−ストークス方程式（Navier-Stokes equations））とを、例えば、有限体積法により差分化し、ＨＳＭＡＣ法により解を求める。求める解は、炉下部の各位置（各計算セル）、各時間ステップにおける流体の速度ベクトル（ｕ_r，ｕθ，ｕ_z）と圧力ｐである。

【0039】

【数1】

【0040】

ここで、ｒは、高炉内の炉下部の半径方向の座標（ｍ）、θは、高炉の軸を原点とする円周方向の角度（ｒａｄ）、ｚは、高炉の鉛直方向の座標（ｍ）、ｕ_rは、半径方向における流体の速度（ｍ／ｓ）、ｕ_zは、鉛直方向における流体の速度（ｍ／ｓ）、ｕθは、円周方向における流体の速度（ｍ／ｓ）、εは、空隙率（−）、Ｓは、流体と充填層との相互作用力（ｍ／ｓ²）、ρは、密度（ｋｇ／ｍ³）、ｐは、圧力（Ｐａ）、ｔは、時間（ｓ）、ｇは、重力加速度（ｍ／ｓ²）である。
生成項Ｓ（流体と充填層との相互作用力）は、例えば、以下の（５）式で表すことができる。

【0041】

【数2】

【0042】

ここで、μは、流体の粘度（Ｐａ・ｓ）、ｄ_pは、コークスの粒径（ｍ）、ｕは、流体の速度ベクトル（ｕ_r，ｕθ，ｕ_z）である。

【0043】

その際の境界条件として、以下の（Ｉ）および（II）を設定する。
（Ｉ）溶銑・溶融スラグの生成速度
炉下部に溜まっている溶融スラグの表面（溶融スラグと気体との界面）における溶銑・溶融スラグの滴下速度（すなわち、溶銑・溶融スラグの生成速度）を境界条件として設定する。溶銑・溶融スラグの生成速度は、溶銑と溶融スラグとで個別に設定される。
例えば、炉頂から旋回シュートを介して単位時間に炉内に装入される原料（焼結鉱、ペレット、塊鉱石等の鉄源を含む原料）の体積の、炉の半径方向の分布に基づいて、溶銑・溶融スラグの生成速度を設定することができる。溶銑・溶融スラグの生成速度は、例えば、非特許文献２等に記載されている公知の高炉数学モデルにより求めることができる。尚、溶銑・溶融スラグの生成速度に、炉の半径方向の分布を持たせず、溶銑・溶融スラグの生成速度の平均値を設定してもよい。尚、以下の説明では、溶銑・溶融スラグの生成速度を必要に応じて銑滓生成速度と称する。

【0044】

（II）溶銑・溶融スラグの排出速度
溶銑・溶融スラグの排出速度ｖ（ｍ³／ｓ）を境界条件として設定する。溶銑・溶融スラグの排出速度ｖは、溶銑と溶融スラグのトータルの排出速度である。本実施形態では、出銑口を、直径Ｄ（ｍ）、出銑口深度Ｌ（ｍ）、表面粗さｅ（ｍ）の円管であると仮定し、以下の（６）式により、溶銑・溶融スラグの排出速度ｖを設定する。

【0045】

【数3】

【0046】

（６）式において、ｆ_cは、摩擦係数（−）であり、例えば、以下の（７）式で表される。また、出銑口の直径Ｄ（ｍ）は、例えば、以下の（８）式で表される。また、ρは、密度（ｋｇ／ｍ³）であり、ΔＰは、出銑口の炉内側の圧力と炉外側の圧力（大気圧）との圧力差（Ｐａ）である。
本実施形態では、出銑口深度Ｌと表面粗さｅを一定であるとする。

【0047】

【数4】

【0048】

（７）式において、Ｒｅは、レイノルズ数（−）である。また、（８）式において、Ｄ₀（ｍ）は、初期出銑口径（出銑口の直径の初期値）であり、Ｔ（分）は、出銑口の出銑作業（開口作業）を開始してからの経過時間である。
その他、炉内の気体の圧力は、送風圧力により定められ、炉外の圧力は、大気圧であるものとする。

【0049】

以上の（１）式〜（８）式に基づいて、出銑口が開口した後の各時間ステップにおいて、炉下部の各位置（各計算セル）における流体の速度ベクトル（ｕ_r，ｕθ，ｕ_z）と圧力ｐとを計算する。
本実施形態では、ラップ時間Ｔ_Bと初期出銑口径Ｄ₀とを指定して、炉下部の各位置（各計算セル）における流体の速度ベクトル（ｕ_r，ｕθ，ｕ_z）と圧力ｐとを、各時間ステップにおいて計算する。

【0050】

このラップ時間Ｔ_Bにより、表２に示す「溜め出し出銑」、「タッチ出銑」、「ラップ出銑」および「パラレル出銑」の何れで操業されるのかが定められる。ラップ時間Ｔ_Bが負の値である場合には「溜め出し出銑」となり、ラップ時間Ｔ_Bが０（ゼロ）である場合には「タッチ出銑」となり、ラップ時間Ｔ_Bが正の値である場合には「ラップ出銑」または「パラレル出銑」になる。

【0051】

本実施形態では、各出銑口１ＴＨ〜３ＴＨを開口するタイミングを、例えば、以下のようにして定める。
すなわち、現在、出銑している出銑口の出銑開始時刻からの経過時間と、指定されたラップ時間Ｔ_Bとの和が、前回の当該出銑口の出銑時間と一致すると、次の出銑口を開口する。例えば、図３（ｃ）において、現在、出銑している出銑口が出銑口１ＴＨであり、その出銑開始時刻がｔ₀であるとする。この場合、出銑開始時刻ｔ₀からの経過時間とラップ時間Ｔ_Bとの和が、前回の出銑口１ＴＨの出銑時間Ｔ_A1になる時刻ｔ_sに、出銑口２ＴＨの開口を開始する。尚、図３（ｃ）では、出銑口１ＴＨの前回の出銑時間と今回の出銑時間とが共にＴ_Aである場合を例に挙げて示す。

【0052】

以上のようにして、出銑口が開口した後の各時間ステップにおいて、炉下部の各位置（各計算セル）における流体の速度ベクトル（ｕ_r，ｕθ，ｕ_z）と圧力ｐとを計算する。本実施形態では、流体の速度ベクトル（ｕ_r，ｕθ，ｕ_z）を、溶銑と溶融スラグとで区別されないもの（１つの流体の速度ベクトル）とする。

【0053】

出銑口が開口した後の各時間ステップにおける流体の速度ベクトル（ｕ_r，ｕθ，ｕ_z）が得られると、当該各時間ステップにおける流体の速度ベクトル（ｕ_r，ｕθ，ｕ_z）に基づいて、当該各時間ステップにおける流体の流速分布を計算することができる。さらに、当該各時間ステップにおける流体の流速分布に基づいて、溶融スラグと気体との界面の位置と、溶融スラグと溶銑との界面の位置とを、当該各時間ステップにおいて計算することができる。

【0054】

本実施形態では、ＶＯＦ（Volume Of Fluid）法を用いて、溶融スラグと気体との界面の位置と、溶融スラグと溶銑との界面の位置とを計算する。ＶＯＦ法は、公知の技術であるが、その概要を簡単に説明する。
ＶＯＦ法とは、或る流体ｉの体積存在率ｆ_i（−）を用いて表現した移流方程式を解いた結果に基づいて界面を移動させ、界面の形状を、或る流体とそれ以外の流体の境界線（または境界面）として間接的に追跡する方法である。体積存在率ｆ_iは０から１までの間の値をとる。

【0055】

或る流体ｉの或る計算セルにおける体積存在率ｆ_iが「１」である場合には、当該計算セルには当該流体ｉのみが存在することを示す。或る流体ｉの或る計算セルにおける体積存在率ｆ_iが「０（ゼロ）」である場合には、当該計算セルには当該流体ｉが存在しないことを示す。或る流体ｉの或る計算セルにおける体積存在率ｆ_iが「０（ゼロ）」を上回り「１」を下回る場合には、当該計算セルには当該流体ｉとそれ以外の流体とが存在することを示す。

【0056】

本実施形態では、変数ｉを、溶銑および溶融スラグを示す変数とする。具体的に変数ｉは「１」または「２」をとる。変数ｉが「１」であることは、溶銑であることを示し、変数ｉが「２」であることは、溶融スラグであることを示す。前述したように、溶銑の比重は溶融スラグの比重よりも重いため、溶銑の方が溶融スラグよりも下に位置する。
例えば、或る計算セルにおいて、溶銑の体積存在率ｆ₁が０．３であるとする。この場合、当該計算セルの最下部から、当該計算セルの高さの３０（％）の高さの領域まで溶銑が存在し、残りの領域に溶融スラグが存在することになる。すなわち、当該計算セルに溶銑と溶融スラグの界面が存在することになる。

【0057】

本実施形態では、各計算セルにおける、溶銑の体積存在率ｆ₁と溶融スラグの体積存在率ｆ₂とを、以下の（９）式の計算を行うことにより計算する。

【0058】

【数5】

【0059】

（９）式において、Ｒ_iは、体積消滅速度（１／ｓ）であり、ｔは、時間（ｓ）である。体積消滅速度Ｒ_iの境界条件としては、例えば、前述した溶銑・溶融スラグの生成速度を用いることができる。化学反応により生成・消滅する溶銑・溶融スラグを考慮して、体積消滅速度Ｒ_iを求めてもよい。
尚、本実施形態では、非特許文献１と同様に、（９）式を計算する際には、ドナー・アクセプター法を用いて離散化する。
また、前述したように、本実施形態では、流体の速度ベクトル（ｕ_r，ｕθ，ｕ_z）を、溶銑と溶融スラグとで区別されないもの（１つの流体の速度ベクトル）とする。すなわち、同一の計算セルにおいて、溶銑の速度ベクトルと溶融スラグの速度ベクトルとを個別に計算しない。

【0060】

以上のようにして計算された溶銑の体積存在率ｆ₁が「０（ゼロ）」になる高さ位置を、炉底から炉頂に向けて探索し、この高さ位置を、溶融スラグと溶銑との界面の位置とする。さらに、溶銑の体積存在率ｆ₁と溶融スラグの体積存在率ｆ₂との和が「０（ゼロ）」になる（すなわち、溶銑の体積存在率ｆ₁と溶融スラグの体積存在率ｆ₂とが共に「０（ゼロ）」になる）高さ位置を、炉底から炉頂に向けて探索し、この高さ位置を、溶融スラグと気体との界面の位置とする。

【0061】

以上のような探索を、炉の円周方向および半径方向の位置（計算セル）を異ならせて行うことにより、或る時間ステップにおける炉下部全体の「「溶融スラグと溶銑との界面の位置」と「溶融スラグと気体との界面の位置」」とが得られる。そして、各時間ステップにおいて以上の計算を行うことにより、各時間ステップにおける炉下部全体の「「溶融スラグと溶銑との界面の位置」と「溶融スラグと気体との界面の位置」」とが得られる。

【0062】

本実施形態では、このようにして得られる「溶融スラグと気体との界面の位置」に基づいて、各出銑口１ＴＨ〜３ＴＨを閉塞するタイミングを定める。
具体的には、溶融スラグと気体との界面の位置が、開口している出銑口の高さ位置であるか否かを判定する。この判定の結果、溶融スラグと気体との界面の位置が、開口している出銑口の高さ位置であると、当該出銑口を閉塞する。出銑口を閉塞する場合には、当該出銑口における溶銑・溶融スラグの排出速度ｖとして０（ゼロ）が設定される。

【0063】

ここで、操業上の安全性を確保する観点から、開口している出銑口の高さ位置は、当該出銑口の炉内側開口部の最上端の高さ位置であるのが好ましい。ただし、当該出銑口の炉内側開口部の最下端から最上端までの高さ位置を、開口している出銑口の高さ位置とすることができる。

【0064】

尚、本実施形態では、流体の速度ベクトル（ｕ_r，ｕθ，ｕ_z）の計算セルに関しては、体積存在率ｆ₁、ｆ₂の計算セルに対して半分ずれたセル（スタガード格子）を使用する。ただし、これらの計算セルを同じにしてもよい。また、炉下部の各位置における流体の速度ベクトル（ｕ_r，ｕθ，ｕ_z）と圧力ｐとを計算する時間ステップと、溶銑の体積存在率ｆ₁と溶融スラグの体積存在率ｆ₂を計算する時間ステップとを同じにする。

【0065】

そして、炉下部の各位置における流体の速度ベクトル（ｕ_r，ｕθ，ｕ_z）および圧力ｐを計算する際には、各計算セルにおける溶銑の体積存在率ｆ₁または溶融スラグの体積存在率ｆ₂に基づいて、各計算セルにおける溶銑と溶融スラグの存在割合を計算し、当該計算割合に基づいて、各計算セルにおける次の時間ステップの密度ρおよび粘度μを計算する。

【0066】

また、炉内の初期状態は、予め設定された状態であるものとする。例えば、前述した（Ｉ）溶銑・溶融スラグの生成速度に従って、溶銑と溶融スラグを炉下部に一定時間滴下させた状態とし、前述したＶＯＦ法を用いて、開口する出銑口に最も近い計算セルにおける溶銑と溶融スラグの存在割合を計算する。そして、当該計算した存在割合に基づいて、当該計算セルにおける密度ρを計算し、その値と初期出銑口径Ｄ₀を、前述した（６）式および（７）式に与えて、当該出銑口における溶銑・溶融スラグの排出速度ｖの初期値を境界条件として設定する。

【0067】

以降の各時間ステップにおいても、初期状態と同様に、開口している出銑口に最も近い計算セルにおける溶銑と溶融スラグの存在割合を計算する。そして、当該計算した存在割合に基づいて、当該計算セルにおける密度ρを計算し、その値と前述した（８）式に示す出銑口の直径Ｄを、前述した（６）式および（７）式に与えて、当該出銑口における溶銑・溶融スラグの排出速度ｖを境界条件として設定する。尚、（８）式に示すように、時間Ｔの経過と共に、出銑口の直径Ｄは拡大する。また、以下の説明では、溶銑・溶融スラグの排出速度を必要に応じて出銑滓速度と称する。

【0068】

以上のように本実施形態では、複数の出銑口１ＴＨ〜３ＴＨのそれぞれに対応して、（１）式〜（９）式を適用する。すなわち、それぞれの出銑口における出銑を開始するタイミングを異ならせ、それぞれの出銑口における溶銑・溶融スラグの排出速度ｖを導出する。このとき、複数の出銑口が開口している場合には、それら複数の出銑口における出銑滓速度を同じ時間ステップにおいて導出する。
そして、それぞれの出銑口における溶銑・溶融スラグの排出速度ｖを境界条件として導出された各計算セルにおける流体の速度ベクトル（ｕ_r，ｕθ，ｕ_z）に基づいて、溶銑・溶融スラグの体積存在率ｆ_iを導出する。この溶銑・溶融スラグの体積存在率ｆ_iは、出銑口の位置（隣り合う出銑口における、出銑口の中心と炉の軸とを結ぶ直線のなす角度）により影響を受ける。

【0069】

そこで、本実施形態では、それぞれの出銑口１ＴＨ〜３ＴＨに対して（６）式〜（９）式を適用可能とするために、それぞれの出銑口の情報（圧力差ΔＰ、出銑口の直径Ｄ、密度ρ、出銑口深度Ｌ、出銑口の表面粗さｅ、出銑開始からの経過時間、前回の出銑時間、開口・閉塞状況等）を個別に管理する機能を、スラグレベル推定モデルに追加した。

【0070】

（炉下部の状態の計算結果の概要）
本実施形態のスラグレベル推定モデルにより計算される炉下部の溶銑と溶融スラグの状態の概要の一例を説明する。
図４は、タッチ出銑の場合の炉下部の溶銑と溶融スラグの状態の一例を示す図である。図４は、本実施形態のスラグレベル推定モデルによる計算の結果から得られたものである。
溶銑の粘度は、溶融スラグの粘度に対し、ほぼ１００倍小さいので、粘度の小さな溶銑が優先して出銑口１ＴＨから排出される。その結果、図４に示す位置Ｄは、位置Ｃより低い。すなわち、溶銑と溶融スラグとの界面の高さは、出銑口１ＴＨから遠くなるほど、低くなる。一方、位置Ａは、位置Ｂより大きく高い位置にある。すなわち、溶融スラグと気体との界面の高さ（スラグレベル）は、出銑口１ＴＨから遠くなるほど、高くなる。そして、出銑末期には、溶融スラグと気体との界面の高さ（スラグレベル）は出銑口１ＴＨに向かって大きく傾斜し、この傾斜によって、出銑終了時に炉内に残留する溶融スラグの量が決まる。

【0071】

図５は、パラレル出銑の場合の炉下部の溶銑と溶融スラグの状態の一例を示す図である。図５も、図４と同様に、本実施形態のスラグレベル推定モデルによる計算の結果から得られたものである。
本実施形態のスラグレベル推定モデルの計算の結果から、本発明者らは、以下の知見を得た。
まず、炉下部の溶銑と溶融スラグの流動・排出の挙動には、複数の出銑口からの出銑滓速度と、出銑口の位置関係とが複雑に影響するという知見を得た。また、炉下部の溶銑との溶融スラグ体積存在率ｆ_iは、それぞれの出銑口からの出銑滓速度、コークスの粒径、溶融スラグの粘度等に依存し、出銑滓速度の増加およびコークスの粒径の低下に伴い、溶融スラグと気体との界面の高さ（スラグレベル）の傾きは大きくなるという知見を得た。

【0072】

図６は、出銑を開始してからのスラグレベルの経時変化の一例を示す図である。具体的に図６（Ａ）は、タッチ出銑の場合のスラグレベルの経時変化の一例を示す。また、図６（Ｂ）は、ラップ出銑の場合のスラグレベルの経時変化の一例を示す。また、図６（Ｃ）は、パラレル出銑の場合のスラグレベルの経時変化の一例を示す図である。図６も、図４および図５と同様に、本実施形態のスラグレベル推定モデルによる計算の結果から得られたものである。

【0073】

高炉の出銑口の開口は、出銑口の閉塞時に充填したマッド材に対して開口機を挿入することにより行われる。初期の出銑滓速度は小さいが、出銑を開始してからの時間の経過に伴うマッド材の侵食により、出銑口の直径は大きくなり、出銑滓速度は徐々に速くなる。
図６（Ａ）に示すように、タッチ出銑の場合、出銑初期の出銑滓速度ｖ₁は小さく、出銑滓速度ｖ₁は、銑滓生成速度Ｐよりも下回る。このため、溶銑・溶融スラグは炉内に貯留する（出銑時間Ｔ₁における出銑滓速度ｖ₁と銑滓生成速度Ｐを参照）。
その後、時間の経過に従い、出銑口の直径は徐々に拡大する。この出銑口の直径の拡大に伴い、出銑滓速度も上昇し、時刻Ｔ_Cで、出銑滓速度ｖ₁と銑滓生成速度Ｐはバランスし、スラグレベルは、最大スラグレベルＭ₁になる。時刻Ｔ_Cを過ぎると出銑滓速度ｖ₁は銑滓生成速度Ｐを上回るので、スラグレベルは、減少に転ずる。その後、スラグレベルは徐々に低下し、スラグレベルが出銑口の高さ位置に到達した時点で、ガスが噴き出し、当該出銑口からの出銑が終了となる。

【0074】

銑滓生成速度と出銑滓速度との差を小さくするために、初期出銑口径を拡大するとともに、マッド材の溶損速度を低下させ、出銑口の直径の拡大速度を低下させる方法もある。しかしながら、マッド材の溶損速度を低下させるためには、マッド材の改良や、マッド材の焼成時間の確保等、解決すべき多くの問題があり容易ではない。よって、マッド材の溶損速度の低下に依らない別の方法で、出銑初期と末期の出銑滓速度の差を低減する必要が望まれる。

【0075】

出銑初期と末期の出銑滓速度の差を低減する方策として、複数の出銑口から同時に出銑するラップ出銑およびパラレル出銑がある。
図６（Ｂ）に示すように、ラップ出銑では、出銑口１ＴＨからの出銑の末期と、出銑口２ＴＨからの出銑の初期とを一部ラップさせることにより、見かけ上、出銑初期の出銑滓速度をタッチ出銑の場合よりも向上させることができる。出銑初期の出銑滓速度Ｖ₁を向上させることにより、出銑初期の出銑滓速度Ｖ₁と銑滓生成速度Ｐとの差を、タッチ出銑の場合よりも小さくすることができ、炉内に貯留する溶銑・溶融スラグの量を少なくすることができる（図６（Ａ）および図６（Ｂ）の出銑時間Ｔ₁における出銑滓速度ｖ₁と銑滓生成速度Ｐとの差を参照）。尚、図６（Ｂ）において、ｎ回目、ｍ回目、ｋ回目は、それぞれ、出銑口１ＴＨ、２ＴＨ、３ＴＨごとに個別に数えた場合の出銑回数を示す。また、Ｖ_s1、Ｖ_s2、Ｖ_s3は、それぞれ、出銑口１ＴＨ、２ＴＨ、３ＴＨにおける出銑初期の出銑滓速度を示し、Ｖ_f1、Ｖ_f2、Ｖ_f3は、それぞれ、出銑口１ＴＨ、２ＴＨ、３ＴＨにおける出銑末期の出銑滓速度を示す。これらの符号の意味は、図６（Ｃ）でも同じである。

【0076】

図６（Ｃ）に示すように、パラレル出銑では、高炉全体の出銑滓速度は、出銑口１ＴＨからの出銑滓速度と出銑口２ＴＨからの出銑滓速度との合計である。したがって、タッチ出銑時の初期出銑口径よりも小さな初期出銑口径を用いることにより、パラレル出銑では、見かけ上、出銑初期の出銑滓速度を、ラップ出銑の場合よりも向上させ、出銑末期の出銑滓速度を、ラップ出銑の場合よりも低下させることができる。なぜなら、出銑初期の出銑滓速度ｖ₁は、出銑口２ＴＨの初期の出銑滓速度Ｖ_S2と、出銑口１ＴＨの中期の出銑滓速度Ｖ_m1との合計だからである。同様に、出銑末期の出銑滓速度ｖ₂も、ラップ出銑の場合よりも低下させることができる。
出銑滓速度の変化を低減することにより、出銑終了時に炉内に残る溶銑・溶融スラグの量が低下すると共に、出銑初期と末期の出銑滓速度の差が低下するため、最大スラグレベルの低下と、スラグレベルの変動幅の低減とが可能となる。

【0077】

以上のように、最大スラグレベルの低下と、スラグレベルの変動幅の低減との双方を実現するという観点から、溜め出し出銑、タッチ出銑、ラップ出銑、およびパラレル出銑のうち、最適な解は、パラレル出銑になる。しかしながら、例えば、所定期間（例えば１日当たり）の出銑回数が多すぎて現実的な出銑回数にならない場合には、パラレル出銑以外の出銑を選択する必要がある。また、最大スラグレベルやスラグレベルの変動幅がある程度の大きさであれば許容できる場合には、必ずしもパラレス出銑を選択する必要はない。

【0078】

そこで、本実施形態では、所定期間（例えば１日当たり）の出銑回数の指定と、最大スラグレベルの指定と、スラグレベルの変動幅の指定とを行うことができるようにする。また、例えば、これらの指定が必要ないような場合のために、初期出銑口径Ｄ₀の上限値の指定も行うことができるようにする。そして、これらの指定に基づいて、高炉の炉下部の状態を解析し、その解析の結果に基づいて、出銑方法（ラップ時間Ｔ_Bと初期出銑口径Ｄ₀）を決定する。

【0079】

ここで、高炉の大型化に伴う出銑量の増加に対し、複数の出銑口から同時に出銑するラップ出銑の形態を採用する場合、ラップ時間の増加に伴う出銑時間の低下、すなわち１日当たりの出銑回数の増加が問題となる。出銑回数の増加は、開口機による開口作業、マッドガンによる閉塞作業等の炉前作業頻度の増加はもとより、出銑口の閉塞に必要なマッド材の使用量の増加、すなわち、溶銑コストの上昇を引き起こす。このため、１日当たりの出銑回数は一定回数以下に抑えることが望ましい。
しかしながら、一般に、最大スラグレベルおよびスラグレベルの変動幅は、出銑回数の増加に伴い低下する。このため、炉前作業頻度、溶銑コスト、最大スラグレベルおよびスラグレベルの変動幅の全てを考慮した上で、初期出銑口径Ｄ₀およびラップ時間Ｔ_Bを決定する必要がある。
尚、１日あたりの出銑回数を低減するために、ラップ出銑時およびパラレル出銑時の初期出銑口直径を、タッチ出銑時の出銑滓速度の経時変化に基づき、タッチ出銑時よりも小さな値に設定するのが好ましい。
以上のような観点で、本実施形態の高炉状態解析装置７００を構成するに至った。

【0080】

（高炉状態解析装置７００の処理）
図７は、高炉状態解析装置７００の機能的な構成の一例を示す図である。高炉状態解析装置７００のハードウェアの構成は、例えば、ＣＰＵ、ＲＯＭ、ＲＡＭ、ＨＤＤ、および各種のインターフェースを用いることにより実現される。また、専用のハードウェアを構築して高炉状態解析装置７００を構成してもよい。
図８〜図１１のフローチャートを参照しながら、高炉状態解析装置７００の処理の一例を説明する。図８〜図１１のフローチャートは、例えば、ＣＰＵがＲＯＭに記憶されたプログラムを実行することにより実現される。

【0081】

図８は、高炉状態解析装置７００の全体処理の一例を説明するフローチャートである。前述したように、ここでは、図２に示す３つの出銑口１ＴＨ、２ＴＨ、３ＴＨをこの順で開口・閉塞して出銑作業を行う場合を例に挙げて示す。また、高炉の炉下部の１日における状態を解析し、その結果に基づいて出銑方法を決定する場合を例に挙げて説明する。また、前述したように、説明を簡単にするために、ここでは、温度に関わる部分の説明を省略する。

【0082】

ステップＳ８０１において、計算条件設定部７０１は、スラグレベル推定モデルによる計算に必要な情報である計算条件を設定する。
本実施形態では、計算条件設定部７０１は、高炉状態解析装置７００のユーザインターフェースの操作等に基づいて、以下の情報を計算条件として入力して設定する。

【0083】

・溶銑の１日当たりの生成量
・溶融スラグの１日当たりの生成量
・溶銑の単位時間当たりの生成量の半径方向の分布（炉内の半径方向の各位置における溶銑の単位時間当たりの生成量）
・溶融スラグの単位時間当たりの生成量の半径方向の分布（炉内の半径方向の各位置における溶融スラグの単位時間当たりの生成量）
・空隙率ε
・コークスの粒径ｄ_p
・溶銑・溶融スラグの密度、粘度
・送風圧力
・出銑口深度Ｌ
・出銑口の表面粗さｅ
・出銑口１ＴＨ、２ＴＨ、３ＴＨの位置
・出銑口１ＴＨ、２ＴＨ、３ＴＨの開口順
・１日当たりの出銑回数の指定の有無
・１日当たりの出銑回数の指定がある場合の指定回数
・最大スラグレベルの上限値の指定の有無
・最大スラグレベルの上限値の指定がある場合の当該上限値
・スラグレベルの変動幅の上限値の指定の有無
・スラグレベルの変動幅の上限値の指定がある場合の当該上限値
・初期出銑口径Ｄ₀の上限値の指定の有無
・初期出銑口径Ｄ₀の上限値の指定がある場合の当該上限値

【0084】

尚、これらの情報は、代表的なものであり、これらの情報以外にも種々の情報を計算条件に含めることができる。例えば、ラップ時間Ｔ_Bと初期出銑口径Ｄ₀等の初期値を計算条件に含めることができる。また、コークスの粒径および空隙率以外にも、炉下部の充填構造に関する情報が必要になる場合には、当該情報を計算条件に含めることができる。また、温度を考慮する場合には、生成される溶銑・溶融スラグの温度に関する情報を計算条件に含めることができる。

【0085】

次に、ステップＳ８０２において、炉下部状態解析部７０２は、ステップＳ８０１で設定された計算条件を用いて、炉下部状態解析処理を行う。炉下部状態解析処理は、前述したスラグレベル推定モデルを用いることにより、１日における炉下部の溶銑・溶融スラグの状態の変動を、ラップ時間Ｔ_Bと初期出銑口径Ｄ₀とを指定して計算する処理をいう。ステップＳ８０２の処理の詳細については、図９を参照しながら、後述する。

【0086】

次に、ステップＳ８０３において、収束判定部７０３は、ステップＳ８０２における炉下部状態解析処理の結果に基づいて、炉下部の状態が収束したか否かを判定する。ステップＳ８０３の判定の詳細については、図１０を参照しながら、後述する。
ステップＳ８０３の結果、炉下部の状態が収束していない場合には、ステップＳ８０４に進む。ステップＳ８０４に進むと、出銑方法変更部７０４は、ラップ時間Ｔ_Bまたは初期出銑口径Ｄ₀を変更する出銑方法変更処理を行う。そして、ステップＳ８０２に戻り、変更後のラップ時間Ｔ_Bまたは初期出銑口径Ｄ₀を用いて、炉下部状態解析処理を再び行う。ステップＳ８０４の処理の詳細については、図１１を参照しながら、後述する。

【0087】

一方、ステップＳ８０３の判定の結果、炉下部の状態が収束した場合には、ステップＳ８０５に進む。ステップＳ８０５に進むと、出力部７０５は、ステップＳ８０３で炉下部の状態が収束したと判定されたときの出銑方法（ラップ時間Ｔ_Bおよび初期出銑口径Ｄ₀）を示す情報を出力する。この他、出力部７０５は、ステップＳ８０３で炉下部の状態が収束したと判定されたときのステップＳ８０２における炉下部状態解析処理の結果を示す情報を出力してもよい。例えば、最大スラグレベルやスラグレベルの変動幅の時間変化を示す情報として、炉下部を模擬したグラフィックスや、グラフを出力してもよい。
そして、図８のフローチャートによる処理を終了する。

【0088】

次に、図９のフローチャートを参照しながら、図８のステップＳ８０２の処理（炉下部状態解析処理）の詳細の一例を説明する。
ステップＳ９０１において、炉下部状態解析部７０２は、出銑口を開口する条件を満足したか否かを判定する。前述したように、本実施形態では、現在、出銑している出銑口の出銑開始時刻ｔ₀からの経過時間と、指定されたラップ時間Ｔ_Bとの和が、前回の当該出銑口の出銑時間Ｔ_A1になる時刻ｔ_sになると、次の出銑口を開口する（図３（ｃ）を参照）。

【0089】

尚、最初にステップＳ９０１の判定を行う際には、炉下部における溶銑・溶融スラグの状態（速度分布、圧力分布、溶銑・溶融スラグの体積存在率ｆ₁、ｆ₂、界面の位置）が、予め設定された基準状態になっているものとし、且つ、出銑口１ＴＨを開口する条件を満足したと判定するものとする。

【0090】

この判定の結果、出銑口を開口する条件を満足していない場合には、ステップＳ９０２を省略して後述するステップＳ９０３に進む。
一方、出銑口を開口する条件を満足した場合には、ステップＳ９０２に進む。ステップＳ９０２に進むと、炉下部状態解析部７０２は、ステップＳ９０１で開口する条件を満足した出銑口を開口する処理を行う。ステップＳ９０２の具体例を説明すると、炉下部状態解析部７０２は、開口する条件を満足した出銑口に対して初期出銑口径Ｄ₀を与える。そして、ステップＳ９０３に進む。

【0091】

ステップＳ９０３に進むと、炉下部状態解析部７０２は、現在の時間ステップにおける流体の速度ベクトル（ｕ_r，ｕθ，ｕ_z）を導出し（（１）式〜（８）式を参照）、現在の時間ステップにおける炉下部における流体の速度分布を導出する。尚、このステップＳ９０３において、炉下部状態解析部７０２は、現在の時間ステップにおける炉下部における圧力分布も導出する。

【0092】

次に、ステップＳ９０４において、炉下部状態解析部７０２は、ステップＳ９０３で導出された速度分布に基づいて、溶銑の体積存在率ｆ₁と溶融スラグの体積存在率ｆ₂とを移流させ、次の時間ステップにおける溶銑の体積存在率ｆ₁と溶融スラグの体積存在率ｆ₂とを導出する（（９）式を参照）。これにより、溶融スラグと気体との界面（スラグレベル）と、溶融スラグと溶銑との界面が移動する。
次に、ステップＳ９０５において、炉下部状態解析部７０２は、ステップＳ９０４で導出された溶銑の体積存在率ｆ₁・溶融スラグの体積存在率ｆ₂と基づいて、次の時間ステップにおけるスラグレベル（溶融スラグと気体との界面の位置）を導出する。尚、このステップＳ９０５において、炉下部状態解析部７０２は、次の時間ステップにおける溶融スラグと溶銑との界面の位置も導出する。

【0093】

次に、ステップＳ９０６において、炉下部状態解析部７０２は、（８）式に基づき、次の時間ステップにおける当該出銑口の直径Ｄの計算を行うとともに、ステップＳ９０３で導出された炉下部における流体の流速分布・圧力分布と、ステップＳ９０４で導出された溶銑の体積存在率ｆ₁・溶融スラグの体積存在率ｆ₂とに基づいて、次の時間ステップにおける出銑口からの溶銑・溶融スラグの排出速度ｖを導出する（（６）式〜（８）式を参照）。尚、現在の時間ステップにおいて開口している出銑口が複数ある場合、それら複数の出銑口のそれぞれからの溶銑・溶融スラグの排出速度ｖが個別に導出される。また、溶銑・溶融スラグの排出速度ｖは、次の時間ステップにおけるステップＳ９０３の計算の境界条件として使用される。

【0094】

ステップＳ９０６の具体例を説明すると、炉下部状態解析部７０２は、開口している出銑口に最も近い計算セルにおける溶銑と溶融スラグの存在割合を導出する。次に、炉下部状態解析部７０２は、当該計算した存在割合に基づいて、当該計算セルにおける密度ρを導出する。また、炉下部状態解析部７０２は、ステップＳ９０３で導出された炉下部における圧力分布と大気圧から、圧力差ΔＰを導出する。そして、炉下部状態解析部７０２は、これらの値と、当該出銑口の直径Ｄと、計算条件（出銑口深度Ｌおよび出銑口の表面粗さｅ）とを（６）式および（７）式に与えて、当該出銑口からの溶銑・溶融スラグの排出速度ｖを導出する。

【0095】

ステップＳ９０７において、炉下部状態解析部７０２は、ステップＳ９０５で導出されたスラグレベルに基づいて、スラグレベル（溶融スラグと気体との界面の位置）が、開口している出銑口の高さ位置に到達したか否かを判定する。
この判定の結果、現在のスラグレベルが、開口している出銑口の高さ位置に到達していない場合には、ステップＳ９０８を省略して後述するステップＳ９０９に進む。

【0096】

一方、現在のスラグレベルが、開口している出銑口の高さ位置に到達している場合には、ステップＳ９０８に進む。ステップＳ９０８に進むと、炉下部状態解析部７０２は、当該出銑口を閉塞する処理を行う。ステップＳ９０８の具体例を説明すると、炉下部状態解析部７０２は、当該出銑口の直径Ｄを０（ゼロ）にする。すなわち、炉下部状態解析部７０２は、当該出銑口については、ステップＳ９０６で導出された溶銑・溶融スラグの排出速度ｖを０（ゼロ）にする。そして、ステップＳ９０９に進む。

【0097】

ステップＳ９０９に進むと、炉下部状態解析部７０２は、１日分の計算が終了したか否かを判定する。
ステップＳ９０９の具体例を説明すると、炉下部状態解析部７０２は、溶銑・溶融スラグの排出速度ｖと、出銑口の直径Ｄと、溶銑・スラグの比重と、開口している出銑口に最も近い計算セルにおける溶銑と溶融スラグの存在割合とに基づいて当該時間ステップにおける出銑量を導出する。そして、炉下部状態解析部７０２は、当該導出した当該時間ステップにおける出銑量を、当該時間ステップまでに出銑された出銑量に加えることにより、総出銑量を導出し、総出銑量が、計算条件に含まれる１日当たりの溶銑・溶融スラグの生成量以上になった場合に、１日分の計算が終了したと判定する。

【0098】

この判定の結果、１日分の計算が終了していない場合には、ステップＳ９１０に進む。ステップＳ９１０に進むと、炉下部状態解析部７０２は、時間ステップを１つ進める。そして、ステップＳ９０１に戻り、次の時間ステップにおける処理を行う。
一方、１日分の計算が終了した場合には、ステップＳ９１１に進む。ステップＳ９１１に進むと、炉下部状態解析部７０２は、１日当たりの出銑回数を導出する。ステップＳ９１１の具体例を説明すると、炉下部状態解析部７０２は、ステップＳ９０２で出銑口を開口した回数を出銑回数として導出する。

【0099】

尚、このステップＳ９１１において、炉下部状態解析部７０２は、１つの出銑口における１回当たりの出銑時間と、全ての出銑口における１日当たりの出銑時間（各出銑口における１日当たりの出銑時間の合計）も導出する。尚、１つの出銑口の１回当たりの出銑時間は、当該出銑口についてステップＳ９０２で開口したときの時間ステップから、当該出銑口についてステップＳ９０８で閉塞したときの時間ステップまでの時間である。

【0100】

さらに、炉下部状態解析部７０２は、最大スラグレベルとスラグレベルの変動幅を導出する。例えば、炉下部状態解析部７０２は、ステップＳ９０５で導出したスラグレベルの最大値（高さ位置が最も高い値）を最大スラグレベルとして導出する。また、炉下部状態解析部７０２は、ステップＳ９０５で導出したスラグレベルの最大値から最小値を引いた値をスラグレベルの変動幅として導出する。ただし、前述したように、１日ではなく、出銑パターンの１周期よりも長い所定の期間において、最大スラグレベルとスラグレベルの変動幅を導出してもよい。
そして、図９のフローチャートによる処理を終了する。

【0101】

次に、図１０のフローチャートを参照しながら、図８のステップＳ８０３の処理の詳細の一例を説明する。
ステップＳ１００１において、収束判定部７０３は、図８のステップＳ８０１で設定された計算条件において、最大スラグレベルの上限値が指定されているか否かを判定する。この判定の結果、最大スラグレベルの上限値が指定されていない場合には、後述するステップＳ１００５に進む。
一方、最大スラグレベルの上限値が指定されている場合には、ステップＳ１００２に進む。ステップＳ１００２に進むと、収束判定部７０３は、図９のステップＳ９１１で導出された最大スラグレベルが、最大スラグレベルの上限値以下であるか否かを判定する。この判定の結果、最大スラグレベルが、最大スラグレベルの上限値以下でない場合、収束判定部７０３は、炉下部の状態が収束していないと判定する。そして、図１０のフローチャートによる処理を終了する。

【0102】

一方、最大スラグレベルが、最大スラグレベルの上限値以下である場合には、ステップＳ１００３に進む。ステップＳ１００３に進むと、収束判定部７０３は、図８のステップＳ８０１で設定された計算条件において、スラグレベルの変動幅の上限値が指定されているか否かを判定する。
この判定の結果、スラグレベルの変動幅の上限値が指定されていない場合には、後述するステップＳ１００６に進む。
一方、スラグレベルの変動幅の上限値が指定されている場合には、ステップＳ１００４に進む。ステップＳ１００４に進むと、収束判定部７０３は、図９のステップＳ９１１で導出されたスラグレベルの変動幅が、スラグレベルの変動幅の上限値以下であるか否かを判定する。

【0103】

この判定の結果、スラグレベルの変動幅が、スラグレベルの変動幅の上限値以下でない場合、収束判定部７０３は、炉下部の状態が収束していないと判定する。そして、図１０のフローチャートによる処理を終了する。
一方、スラグレベルの変動幅が、スラグレベルの変動幅の上限値以下である場合、収束判定部７０３は、炉下部の状態が収束していると判定する。そして、図１０のフローチャートによる処理を終了する。

【0104】

前述したように、ステップＳ１００１の判定の結果、最大スラグレベルの上限値が指定されていない場合には、ステップＳ１００５に進む。ステップＳ１００５に進むと、収束判定部７０３は、最大スラグレベルが最小値であるか否かを判定する。
この判定の結果、最大スラグレベルが最小値である場合には、前述したステップＳ１００３に進む。
一方、最大スラグレベルが最小値でない場合、収束判定部７０３は、炉下部の状態が収束していないと判定する。そして、図１０のフローチャートによる処理を終了する。

【0105】

また、前述したように、ステップＳ１００３の判定の結果、スラグレベルの変動幅の上限値が指定されていないと判定された場合には、ステップＳ１００６に進む。ステップＳ１００６に進むと、収束判定部７０３は、スラグレベルの変動幅が最小値であるか否かを判定する。
この判定の結果、スラグレベルの変動幅が最小値である場合、収束判定部７０３は、炉下部の状態が収束していると判定する。そして、図１０のフローチャートによる処理を終了する。
一方、スラグレベルの変動幅が最小値でない場合、収束判定部７０３は、炉下部の状態が収束していないと判定する。そして、図１０のフローチャートによる処理を終了する。

【0106】

ステップＳ１００５（Ｓ１００６）の具体例を説明すると、収束判定部７０３は、最大スラグレベル（スラグレベルの変動幅）が、予め設定されている下限値以下になった場合に、最大スラグレベル（スラグレベルの変動幅）が最小値であると判定することができる。この他、収束判定部７０３は、山登り法やＧＡ（遺伝的アルゴリズム）等の公知の数理計画法における最適化探索手法を用いることにより、最大スラグレベル（スラグレベルの変動幅）の最小値を求めることもできる。

【0107】

次に、図１１のフローチャートを参照しながら、図８のステップＳ８０４の処理（出銑方法変更処理）の詳細の一例を説明する。
ステップＳ１１０１において、出銑方法変更部７０４は、図８のステップＳ８０１で設定された計算条件において、１日当たりの出銑回数が指定されており、且つ、図９のステップＳ９１１で導出された１日当たりの出銑回数が指定回数を上回るか否かを判定する。

【0108】

この判定の結果、１日当たりの出銑回数が指定されていない場合、または、１日当たりの出銑回数が指定回数を上回らない場合には、後述するステップＳ１１０６に進む。
一方、１日当たりの出銑回数が指定されており、且つ、１日当たりの出銑回数が、１日当たりの出銑回数の指定値を上回る場合には、ステップＳ１１０２に進む。

【0109】

ステップＳ１１０２に進むと、出銑方法変更部７０４は、前回のステップＳ８０４（図１１のフローチャート）の処理で、ラップ時間Ｔ_Bを変更したか否かを判定する。この判定の結果、前回のステップＳ８０４の処理で、ラップ時間Ｔ_Bを変更した場合には、ステップＳ１１０３に進む。ステップＳ１１０３に進むと、出銑方法変更部７０４は、初期出銑口径Ｄ₀を１（ｍｍ）縮小する。そして、図１１のフローチャートによる処理を終了する。

【0110】

一方、前回のステップＳ８０４の処理で、ラップ時間Ｔ_Bを変更していない場合には、ステップＳ１１０４に進む。ステップＳ１１０４に進むと、出銑方法変更部７０４は、現在のラップ時間Ｔ_Bが１（ｍｉｎ）以上であるか否かを判定する。この判定の結果、現在のラップ時間Ｔ_Bが１（ｍｉｎ）以上でない場合には、前述したステップＳ１１０３に進み、出銑方法変更部７０４は、初期出銑口径Ｄ₀を１（ｍｍ）縮小し、図１１のフローチャートによる処理を終了する。

【0111】

一方、現在のラップ時間Ｔ_Bが１（ｍｉｎ）以上である場合には、ステップＳ１１０５に進む。ステップＳ１１０５に進むと、出銑方法変更部７０４は、ラップ時間Ｔ_Bを１（ｍｉｎ）短縮する。そして、図１１のフローチャートによる処理を終了する。
前述したように、ステップＳ１１０１において、１日当たりの出銑回数が指定されていない場合、または、１日当たりの出銑回数が指定回数を上回らないと判定された場合には、ステップＳ１１０６に進む。

【0112】

ステップＳ１１０６に進むと、出銑方法変更部７０４は、図８のステップＳ８０１で設定された計算条件において、最大スラグレベルの上限値が指定されていないか、または、図９のステップＳ９１１で導出された最大スラグレベルが、最大スラグレベルの上限値を上回るか、を判定する。この判定の結果、最大スラグレベルの上限値が指定されており、且つ、最大スラグレベルが、最大スラグレベルの上限値を上回らない場合には、後述するステップＳ１１１２に進む。

【0113】

一方、最大スラグレベルの上限値が指定されていない場合、または、最大スラグレベルが、最大スラグレベルの上限値を上回る場合には、ステップＳ１１０７に進む。
ステップＳ１１０７に進むと、出銑方法変更部７０４は、前回のステップＳ８０４（図１１のフローチャート）の処理で、初期出銑口径Ｄ₀を変更したか否かを判定する。この判定の結果、前回のステップＳ８０４の処理で、初期出銑口径Ｄ₀を変更していない場合には、後述するステップＳ１１１０に進む。

【0114】

一方、前回のステップＳ８０４の処理で、初期出銑口径Ｄ₀を変更した場合には、ステップＳ１１０８に進む。ステップＳ１１０８に進むと、出銑方法変更部７０４は、現在のラップ時間Ｔ_Bが、ステップＳ９１１で導出された１つの出銑口における１回当たりの出銑時間の０．５倍を下回るか否かを判定する。この判定の結果、現在のラップ時間Ｔ_Bが、１つの出銑口における１回当たりの出銑時間の０．５倍を下回らない場合には、後述するステップＳ１１１０に進む。

【0115】

一方、現在のラップ時間Ｔ_Bが、１つの出銑口における１回当たりの出銑時間の０．５倍を下回る場合には、ステップＳ１１０９に進む。ステップＳ１１０９に進むと、出銑方法変更部７０４は、ラップ時間Ｔ_Bを１（ｍｉｎ）延長する。そして、図１１のフローチャートによる処理を終了する。
前述したように、ステップＳ１１０７において、前回のステップＳ８０４の処理で、初期出銑口径Ｄ₀を変更していないと判定された場合と、ステップＳ１１０８において、現在のラップ時間Ｔ_Bが、１つの出銑口における１回当たりの出銑時間の０．５倍を下回らないと判定された場合には、ステップＳ１１１０に進む。

【0116】

ステップＳ１１１０に進むと、出銑方法変更部７０４は、図８のステップＳ８０１で設定された計算条件において、初期出銑口径Ｄ₀の上限値が指定されており、且つ、現在の初期出銑口径Ｄ₀が、初期出銑口径Ｄ₀の上限値を下回るか否かを判定する。
この判定の結果、初期出銑口径Ｄ₀の上限値が指定されていない場合、または、現在の初期出銑口径Ｄ₀が、初期出銑口径Ｄ₀の上限値を下回らない場合には、図１１のフローチャートによる処理を終了する。
一方、初期出銑口径Ｄ₀の上限値を下回る場合には、ステップＳ１１１１に進む。ステップＳ１１１１に進むと、出銑方法変更部７０４は、初期出銑口径Ｄ₀を１（ｍｍ）拡大する。そして、図１１のフローチャートによる処理を終了する。

【0117】

前述したように、ステップＳ１１０６において、最大スラグレベルの上限値が指定されており、且つ、最大スラグレベルが、最大スラグレベルの上限値を上回らないと判定された場合には、ステップＳ１１１２に進む。ステップＳ１１１２に進むと、出銑方法変更部７０４は、図８のステップＳ８０１で設定された計算条件において、スラグレベルの変動幅の上限値が指定されていないか、または、図９のステップＳ９１１で導出されたスラグレベルの変動幅が、スラグレベルの変動幅の上限値を上回るか、を判定する。この判定の結果、スラグレベルの変動幅の上限値が指定されており、且つ、スラグレベルの変動幅が、スラグレベルの変動幅の上限値を上回らない場合には、図１１のフローチャートによる処理を終了する。
一方、スラグレベルの変動幅の上限値が指定されていない場合、または、スラグレベルの変動幅が、スラグレベルの変動幅の上限値を上回る場合には、前述したステップＳ１１０７に進む。

【0118】

（実施例）
次に、実施例を説明する。
［実施例１］
内容積が５７７５（ｍ³）の高炉に、コークス比が３００（ｋｇ／ｔ）、微粉炭比が１８０（ｋｇ／ｔ）、還元材比が４８０（ｋｇ／ｔ）、スラグ比が３１３（ｋｇ／ｔ）の条件で原料が装入され、且つ、図２に示す位置関係にある出銑口１ＴＨ〜３ＴＨからこの順でパラレル出銑を行うものとした場合であって、１日当たりの出銑量（ｔ／ｄ）、炉下部のコークスの粒径（ｍｍ）、炉下部の空隙率（％）を変更した場合の、最大スラグレベルおよびスラグレベルの変動幅を、前述した本実施形態のスラグレベル推定モデルを用いて計算した。この際、初期出銑口径Ｄ₀を６０（ｍｍ）として前述した（８）式で出銑口直径Ｄを計算した。また、出銑口深度Ｌを４（ｍ）とし、出銑口の表面粗さｅを１（ｍｍ）とした。
その結果を、表３に示す。出銑口基準の最大スラグレベル（出銑口の高さ位置を基準とした最大スラグレベル）は、０．７７（ｍ）〜１．４６（ｍ）であり、出銑口基準の羽口レベル（出銑口の高さ位置を基準とした送風羽口の高さ位置）である５．３ｍに対し十分に小さいものであり、安定した高炉操業が可能となることが分かる。

【0119】

【表3】

【0120】

（実施例２）
内容積が５７７５（ｍ³）の高炉において、コークス比が３００（ｋｇ／ｔ）、微粉炭比が（１８０ｋｇ／ｔ）、還元材比が４８０（ｋｇ／ｔ）、スラグ比が３１３（ｋｇ／ｔ）の条件で原料が装入され、且つ、図２に示す位置関係にある出銑口１ＴＨ〜３ＴＨからこの順でパラレル出銑を行うものとした場合であって、炉下部のコークスの粒径を３０（ｍｍ）、炉下部の空隙率を３０（％）、１日当たりの出銑回数の指定回数を１２（回）および１０（回）とした場合の、最大スラグレベルおよびスラグレベルの変動幅を、前述した本実施形態のスラグレベル推定モデルを用いて計算した。この際、出銑口基準の最大スラグレベルの上限値を１．７（ｍ）とした。また、初期出銑口径Ｄ₀を６０（ｍｍ）として前述した（８）式で出銑口直径Ｄを計算した。また、出銑口深度Ｌを４（ｍ）とし、出銑口の表面粗さｅを１（ｍｍ）とした。

【0121】

その結果を表４に示す。１日当たりの出銑回数の指定回数が１２（回）の場合には、初期出銑口径が７４（ｍｍ）、ラップ時間が１３（ｍｉｎ）となった。また、１日当たりの出銑回数の指定回数が１０（回）の場合には、初期出銑口径が７６（ｍｍ）、ラップ時間が６（ｍｉｎ）となった。何れの場合も、出銑口基準の最大スラグレベルが１．７ｍ以下となり、安定した高炉操業が可能となることが分かる。

【0122】

【表4】

【0123】

（まとめ）
前述したように従来は、炉内で生成される溶銑・溶融スラグの量を、高炉に装入される鉄鉱石の量から算出し、出銑口から排出される溶銑の量を、受銑容器（トーピードカー）の秤量により測定し、出銑口から排出される溶融スラグの量を、水砕のコンベアの秤量により測定していた。しかしながら、秤量のタイミングのずれや、秤量誤差等により正確な管理が困難であった。
また、非特許文献１に記載の技術では、同時に出銑することが可能な出銑口の数が１つである場合の炉下部の状態を計算することはできるが、複数の出銑口から同時に出銑することが可能な出銑形態で出銑を行う場合の炉下部の状態の計算はできない。

【0124】

これに対して本実施形態では、炉下部における溶銑・溶融スラグの流速の分布を各時間ステップにおいて導出し、導出した流速の分布に基づいて、溶銑・溶融スラグを移動し、移動後のスラグレベルを導出する。この際、ラップ時間Ｔ_Bと、出銑作業を行う出銑口１ＴＨ〜３ＴＨの順序とに基づいて各出銑口１ＴＨ〜３ＴＨを開口するタイミングを判定し、出銑口を開口するタイミングになると初期出銑口径Ｄ₀で当該出銑口が開口する状態にする。また、スラグレベルが出銑口の高さ位置に達した場合、当該出銑口が閉塞する状態にする。このようにして１日の操業におけるスラグレベルの変動を導出し、最大スラグレベル、スラグレベルの変動幅、初期出銑口径Ｄ₀とのうち、予め指定されたものが所定の条件を満足するようになるラップ時間Ｔ_Bと初期出銑口径Ｄ₀を探索する。

【0125】

したがって、複数の出銑口から同時に出銑する出銑形態であっても、炉下部の状態を高精度に推定することができる。したがって、炉下部におけるにおける最大スラグレベルとスラグレベルの変動幅を低減することが可能な出銑方法を導出することができる。

【0126】

（変形例）
［変形例１］
本実施形態では、高炉状態解析装置７００が、ラップ時間Ｔ_Bと初期出銑口径Ｄ₀の変更を、図１１に示す所定のアルゴリズムに基づいて自動的に行う場合を例に挙げて説明した（図１１を参照）。しかしながら、炉下部の状態が収束したか否かの判定の結果、炉下部の状態が収束していない場合、オペレータが、高炉状態解析装置７００のユーザインターフェースに対して入力操作を行うことにより、ラップ時間Ｔ_Bと初期出銑口径Ｄ₀の少なくとも何れか一方の変更値を入力してもよい。例えば、図１１のフローチャートによる手順をオペレータ自身の判断で行ってもよいが、その他の手順で変更値を定めてもよい。このことは、ラップ時間Ｔ_Bと初期出銑口径Ｄ₀の変更を高炉状態解析装置７００が自動的に行う場合も同様である。例えば、図１１のステップＳ１１０１の判定を、図８のステップＳ８０２とステップＳ８０３との間で行い、１日当たりの出銑回数が指定されており、且つ、図９のステップＳ９１１で導出された１日当たりの出銑回数が指定回数を上回る場合には、図１１のステップＳ１１０２に進み、そうでない場合であって、ステップＳ８０３で炉下部の状態が収束していないと判定された場合には、図１１のステップＳ１１０６に進むようにしてもよい。

【0127】

［変形例２］
本実施形態では、ラップ時間Ｔ_Bと初期出銑口径Ｄ₀の双方を変更し得る場合を例に挙げて説明した（図１１を参照）。しかしながら、ラップ時間Ｔ_Bと初期出銑口径Ｄ₀の一方を固定し、他方を変更してもよい。本変形例は、高炉状態解析装置７００が変更を自動的に行う場合と、オペレータが変更値の入力を行う場合の双方に適用される。
［変形例３］
本実施形態では、最大スラグレベルの上限値およびスラグレベルの変動幅の上限値を指定し、最大スラグレベルおよびスラグレベルの変動幅が、それぞれ上限値以下であるか否かを判定する場合を例に挙げて説明した（図１０のＳ１００１〜Ｓ１００４を参照）。しかしながら、最大スラグレベルの上限値およびスラグレベルの変動幅の上限値を指定せずに、最大スラグレベルおよびスラグレベルの変動幅が、それぞれ最小値である場合に、炉下部の状態が収束したと判定し、そうでない場合には、炉下部の状態が収束していないと判定してもよい。また、これとは逆に、最大スラグレベルおよびスラグレベルの変動幅が、それぞれ最小値であるか否かを判定せずに、最大スラグレベルおよびスラグレベルの変動幅が、それぞれ上限値以下であるか否かを判定してもよい。このように、スラグレベルの状態が操業上要求される所定の条件を満たすか否かを判定していれば、炉下部の状態が収束したか否かを判定する手法は、本実施形態で説明した手法に限定されない。

【0128】

［変形例４］
本実施形態では、初期出銑口径Ｄ₀、出銑口深度Ｌおよび出銑口の表面粗さｅが出銑口に関わらず同じ値である場合を例に挙げて説明した。しかしながら、これらの少なくとも１つを出銑口ごとに異ならせてもよい。また、初期出銑口径Ｄ₀（ｍ）を、出銑口の直径の初期値としたが、直径の代わりに半径を用いてもよい。このように、初期出銑口径は、出銑口の開口径の初期値であればよい。
［変形例５］
本実施形態では、出銑に利用する出銑口の数が３つである場合を例に挙げて説明した。しかしながら、出銑に利用する出銑口の数は２つ以上であれば、幾つであってもよい。出銑に利用する出銑口の数が２つである場合、出銑形態としてパラレル出銑を含めてもよいが、出銑形態をパラレル出銑にすると、２つの出銑口を常に開口することになる。したがって、出銑に利用する出銑口の数が２つである場合には、出銑形態が、「溜め出し出銑」、「タッチ出銑」および「ラップ出銑」の何れかになるように、ラップ時間Ｔ_Bと初期出銑口径Ｄ₀を探索するのが好ましい。また、出銑に利用する出銑口の数がＮ（Ｎは３以上の整数）である場合には、常に（Ｎ−１）の出銑口から出銑が行われる出銑形態を採用することができる。例えば、出銑に利用する出銑口の数が４つである場合には、常に３つの出銑口から出銑が行われる出銑形態を採用することができる。

【0129】

［変形例６］
本実施形態では、連続の式と運動量保存式を有限体積法により差分化する場合を例に挙げて説明した。しかしながら、炉下部における溶銑・溶融スラグの流速および圧力を数値解析により計算していれば、必ずしも有限体積法を用いる必要はない。例えば、有限差分法を用いてもよい。
［変形例７］
本実施形態では、ＶＯＦ法を用いて、体積存在率ｆ_iを導出し、溶融スラグと気体との界面の高さ位置と、溶融スラグと溶銑との界面の高さ位置とを導出する場合を例に挙げて説明した。しかしながら、界面の高さ位置を移流方程式に基づき数値解析により計算していれば、必ずしも、ＶＯＦ法を用いる必要はない。例えば、ＡＬＥ（Arbitrary Lagrangian and Eulerian）法を用いてもよい。

【0130】

［変形例８］
本実施形態では、１日当たりの出銑回数が指定回数を上回るか否かを判定する場合を例に挙げて説明した（図１１のステップＳ１１０１）。しかしながら、出銑回数を求める期間は、必ずしも１日である必要はなく、１日よりも短い期間でも長い期間でもよい。
［変形例９］
本実施形態では、３次元の解析を行う場合を例に挙げて説明した。しかしながら、例えば、高炉の高さ方向と、当該高炉の高さ方向に垂直な方向とにより定まる２次元の解析を行うようにしてもよい。
［変形例１０］
本実施形態では、炉下部の状態が収束したか否かを判定し、炉下部の状態が収束していない場合には、初期出銑口径Ｄ₀およびラップ時間Ｔ_Bの何れか一方を修正する場合を例に挙げて説明した。しかしながら、必ずしもこれらの処理を行わなくてもよい。例えば、或る初期出銑口径Ｄ₀およびラップ時間Ｔ_Bにおける最大スラグレベルおよびスラグレベルの変動幅を確認することを目的とする場合には、これらの処理を行わなくてもよい。

【0131】

［その他の変形例］
尚、以上説明した本発明の実施形態は、コンピュータがプログラムを実行することによって実現することができる。また、前記プログラムを記録したコンピュータ読み取り可能な記録媒体及び前記プログラム等のコンピュータプログラムプロダクトも本発明の実施形態として適用することができる。記録媒体としては、例えば、フレキシブルディスク、ハードディスク、光ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、磁気テープ、不揮発性のメモリカード、ＲＯＭ等を用いることができる。
また、以上説明した本発明の実施形態は、何れも本発明を実施するにあたっての具体化の例を示したものに過ぎず、これらによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその技術思想、またはその主要な特徴から逸脱することなく、様々な形で実施することができる。

【0132】

（請求項との関係）
［請求項１、８］
設定手段は、例えば、計算条件設定部７０１を用いることにより実現される（図８のステップＳ８０１、および、変形例５も参照）。
開口タイミング判定手段は、例えば、炉下部状態解析部７０２を用いることにより実現される（図９のステップＳ９０１も参照）。
速度分布導出手段は、例えば、炉下部状態解析部７０２を用いることにより実現される（図９のステップＳ９０２、Ｓ９０３、Ｓ９０６、Ｓ９０８、スラグレベル推定モデルの説明（（１）式〜（８）式）、および、変形例６も参照）。
スラグレベル導出手段は、例えば、炉下部状態解析部７０２を用いることにより実現される（図９のステップＳ９０５、Ｓ９０６、スラグレベル推定モデルの説明（（９）式）、および、変形例７も参照）。
閉塞タイミング判定手段は、例えば、炉下部状態解析部７０２を用いることにより実現される（図９のステップＳ９０７、および、スラグレベル推定モデルの説明も参照）。
［請求項２］
収束判定手段は、例えば、収束判定部７０３を用いることにより実現される（図８のステップＳ８０３、スラグレベル推定モデルの説明、および、変形例３、１０も参照）。
［請求項３］
最大スラグレベル導出手段およびスラグレベル変動幅導出手段は、例えば、炉下部状態解析部７０２を用いることにより実現される（図９のステップＳ９１１も参照）。
最大スラグレベル判定手段は、例えば、収束判定部７０３を用いることにより実現される（図１０のステップＳ１００２、Ｓ１００５も参照）。
スラグレベル変動幅判定手段は、例えば、収束判定部７０３を用いることにより実現される（図１０のステップＳ１００４、Ｓ１００６も参照）。
［請求項４、６、７］
出銑方法変更手段は、例えば、出銑方法変更部７０４を用いることにより実現される（図８のステップＳ８０４、図１１、および変形例１、２も参照）。
［請求項５］
出銑回数導出手段は、例えば、炉下部状態解析部７０２を用いることにより実現される（図９のステップＳ９１１、および、変形例８も参照）。
前記出銑方法変更手段は、例えば、出銑方法変更部７０４を用いることにより実現される（図８のステップＳ８０４、図１１のステップＳ１１０１〜１１０５、および変形例１も参照）。

【符号の説明】

【0133】

ＴＨ：出銑口、７００：高炉状態解析装置、７０１：計算条件設定部、７０２：炉下部状態解析部、７０３：収束判定部、７０４：出銑方法変更部、７０５：出力部

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】