(57)【要約】【課題】高炉の炉内の状態をオンラインで高精度に推定するとともに、推定した炉内の状態から現在または将来における操業の異常度を推定することができる操業異常推定方法および操業異常推定装置を提供すること。【解決手段】操業異常推定方法は、炉内に装入物を装入して銑鉄を製造する高炉の操業中の炉内ガスのガス成分の実績に基づいて、炉内における装入物の装入状態を推定し、推定した装入物の装入状態に基づいて、操業中の高炉の融着帯の状態を推定する炉内状態推定ステップと、炉内状態推定ステップで推定した融着帯の形状から、統計的解析によって、高炉の現在または将来における操業の異常度を推定する異常度推定ステップと、を含んでいる。【選択図】図１０
2018003044.tif 000002
【特許請求の範囲】
【請求項１】
炉内に装入物を装入して銑鉄を製造する高炉の操業中の炉内ガスのガス成分の実績に基づいて、前記炉内における前記装入物の装入状態を推定し、推定した前記装入物の装入状態に基づいて、操業中の前記高炉の融着帯の状態を推定する炉内状態推定ステップと、
前記炉内状態推定ステップで推定した前記融着帯の形状から、統計的解析によって、前記高炉の現在または将来における操業の異常度を推定する異常度推定ステップと、
を含むことを特徴とする操業異常推定方法。
【請求項２】
前記異常度推定ステップは、前記融着帯の中心部、炉壁部の層厚、前記径方向における前記融着帯の層厚の平均値、前記融着帯の高さ、前記融着帯の傾斜のうちの少なくとも１つを説明変数とし、前記高炉の現在または将来における通気抵抗指数を目的関数として、回帰分析によって、前記高炉の現在または将来における操業の異常度を推定することを特徴とする請求項１に記載の操業異常推定方法。
【請求項３】
前記炉内状態推定ステップは、前記ガス成分の実績から算出されるガス利用率分布を用いて前記装入物の装入状態を推定することを特徴とする請求項１または請求項２に記載の操業異常推定方法。
【請求項４】
前記炉内状態推定ステップは、前記装入物の装入状態として、前記装入物である鉱石およびコークスの前記径方向に沿った層厚分布を推定することを特徴とする請求項１から請求項３のいずれか一項に記載の操業異常推定方法。
【請求項５】
前記炉内状態推定ステップは、主成分分析および重回帰分析、あるいは部分的最小二乗法を用いて、前記ガス成分の実績データから前記装入物の装入状態を推定することを特徴とする請求項１から請求項４のいずれか一項に記載の操業異常推定方法。
【請求項６】
炉内に装入物を装入して銑鉄を製造する高炉の操業中の炉内ガスのガス成分の実績に基づいて、前記炉内における前記装入物の装入状態を推定し、推定した前記装入物の装入状態に基づいて、操業中の前記高炉の融着帯の状態を推定する炉内状態推定部と、
前記炉内状態推定部が推定した前記融着帯の状態の形状から、統計的解析によって、前記高炉の現在または将来における操業の異常度を推定する異常度推定部と、
を備えることを特徴とする操業異常推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、操業異常推定方法および操業異常推定装置に関する。
【背景技術】
【０００２】
高炉操業の高度化に伴い、鉱石およびコークスからなる装入物の炉内における装入状態（例えば層厚分布）を把握することは非常に重要となっている。装入物の層厚分布を推定する従来技術としては、例えば特許文献１や特許文献２で示される方法がある。
【０００３】
特許文献１は、高炉の炉壁に設置された圧力データや温度データの時間的・空間的推移から融着帯の根部、すなわち壁面に接する部分の融着帯の形状を推定する方法に関するものである。特許文献２は、高炉の炉頂における温度分布とガス流速分布の変化から高炉装入物の層厚分布を推定する方法に関するものである。
【０００４】
また、高炉の炉内状態を推定する従来技術として、例えば非特許文献１で示される方法がある。非特許文献１の推定方法は、高炉を軸対称として２次元のメッシュに分割し、与えられた条件のもとに、固体（鉱石、コークス）の流れ、液体（溶銑）の流れ、ガスの流れ、伝熱および化学反応を計算し、炉内の状態を推定するものである。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００３−１９３１２０号公報
【特許文献２】特開平９−２８７００８号公報
【非特許文献】
【０００６】
【非特許文献１】佐藤健ら、「高炉操業シミュレータの開発と溶銑シリコン低減への適用」、川崎製鉄技報、1997年、Vol.29、p.30-36
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、特許文献１の推定方法は、炉体表面の温度や圧力の計測値から推定を行うため、炉内の形状までは推定することができないという問題がある。また、特許文献２の推定方法は、高炉の炉頂における温度分布とガス流速分布の変化から装入物の層厚分布を推定する方法であるため、当該方法で推定した層厚分布を、例えば非特許文献１で示されているような高炉物理モデルの入力として与えたとしても、炉内の状態を高精度に推定することができないという問題がある。
【０００８】
非特許文献１の推定方法は、オフライン利用のために開発されたものであり、例えば新しい操業形態を実施するとき等の事前検討として利用される。しかし、入力データは必ずしもオンラインで計測または推定できていないため、非特許文献１の技術をオンラインで利用することができないという問題がある。
【０００９】
また、高炉操業にあたっては、炉内の状態をオンラインで高精度に推定するということ以外にも、例えば現在または将来における高炉の操業中に発生する異常を推定（予測）することも望まれているが、前記した特許文献１，２および非特許文献１では、いずれもこのような異常推定を行うことはできない。
【００１０】
本発明は、上記課題に鑑みてなされたものであって、高炉の炉内の状態をオンラインで高精度に推定するとともに、推定した炉内の状態から現在または将来における操業の異常度を推定することができる操業異常推定方法および操業異常推定装置を提供することを目的とする。
【課題を解決するための手段】
【００１１】
上述した課題を解決し、目的を達成するために、本発明に係る操業異常推定方法は、炉内に装入物を装入して銑鉄を製造する高炉の操業中の炉内ガスのガス成分の実績に基づいて、前記炉内における前記装入物の装入状態を推定し、推定した前記装入物の装入状態に基づいて、操業中の前記高炉の融着帯の状態を推定する炉内状態推定ステップと、前記炉内状態推定ステップで推定した前記融着帯の形状から、統計的解析によって、前記高炉の現在または将来における操業の異常度を推定する異常度推定ステップと、を含むことを特徴とする。
【００１２】
また、本発明に係る操業異常推定方法は、上記発明において、前記異常度推定ステップは、前記融着帯の中心部、炉壁部の層厚、前記径方向における前記融着帯の層厚の平均値、前記融着帯の高さ、前記融着帯の傾斜のうちの少なくとも１つを説明変数とし、前記高炉の現在または将来における通気抵抗指数を目的関数として、回帰分析によって、前記高炉の現在または将来における操業の異常度を推定することを特徴とする。
【００１３】
また、本発明に係る操業異常推定方法は、上記発明において、前記炉内状態推定ステップは、前記ガス成分の実績から算出されるガス利用率分布を用いて前記装入物の装入状態を推定することを特徴とする。
【００１４】
また、本発明に係る操業異常推定方法は、上記発明において、前記炉内状態推定ステップは、前記装入物の装入状態として、前記装入物である鉱石およびコークスの前記径方向に沿った層厚分布を推定することを特徴とする。
【００１５】
また、本発明に係る操業異常推定方法は、上記発明において、前記炉内状態推定ステップは、主成分分析および重回帰分析、あるいは部分的最小二乗法を用いて、前記ガス成分の実績データから前記装入物の装入状態を推定することを特徴とする。
【００１６】
上述した課題を解決し、目的を達成するために、本発明に係る操業異常推定装置は、炉内に装入物を装入して銑鉄を製造する高炉の操業中の炉内ガスのガス成分の実績に基づいて、前記炉内における前記装入物の装入状態を推定し、推定した前記装入物の装入状態に基づいて、操業中の前記高炉の融着帯の状態を推定する炉内状態推定部と、前記炉内状態推定部が推定した前記融着帯の状態の形状から、統計的解析によって、前記高炉の現在または将来における操業の異常度を推定する異常度推定部と、を備えることを特徴とする。
【発明の効果】
【００１７】
本発明によれば、高炉の炉内の状態をオンラインで高精度に推定することができるとともに、推定した炉内の状態から現在または将来における操業の異常度を的確かつ早期に推定することができる。
【図面の簡単な説明】
【００１８】
【図１】図１は、本発明の実施形態に係る操業異常推定装置の構成を模式的に示す図である。
【図２】図２は、装入物の層厚分布を説明するための説明図である。
【図３】図３は、ガス利用率分布を説明するための説明図である。
【図４】図４は、高炉物理モデルによるガス利用率分布の計算方法を説明するための説明図である。
【図５】図５は、入力推定逆モデルの構築手順を示すフローチャートである。
【図６】図６は、入力初期値を説明するための説明図である。
【図７】図７は、入力セットを説明するための説明図である。
【図８】図８は、高炉物理モデルによる炉内状態の推定方法を説明するための説明図である。
【図９】図９は、操業異常度としての通気抵抗指数の計算方法を説明するための説明図である。
【図１０】図１０は、本発明の実施形態に係る操業異常推定装置による操業異常推定方法の処理手順を示すフローチャートである。
【図１１】図１１は、本発明の実施形態に係る操業異常推定装置において、計算結果に係るデータの表示方法を示すフローチャートである。
【図１２】図１２は、本発明の実施形態に係る操業異常推定装置による炉内状態の推定結果の一例を示す図である。
【図１３】図１３は、本発明の実施形態に係る操業異常推定装置による炉内状態の推定結果の一例を示す図である。
【図１４】図１４は、本発明の実施形態に係る操業異常推定装置による炉内状態の推定結果の一例を示す図である。
【図１５】図１５は、本発明の実施形態に係る操業異常推定装置による操業異常度の推定結果の一例を示す図である。
【発明を実施するための形態】
【００１９】
以下、本発明に係る操業異常推定方法および操業異常推定装置の実施形態について、図面を参照しながら説明する。なお、本発明は以下の実施形態に限定されるものではない。また、以下の実施形態における構成要素には、当業者が置換可能かつ容易なもの、あるいは実質的に同一のものが含まれる。
【００２０】
[操業異常推定装置]
本実施形態に係る操業異常推定装置の構成について、図１〜図９を参照しながら説明する。操業異常推定装置１は、高炉１００の炉内の状態（以下、「炉内状態」という）をオンラインで推定するとともに、推定した現在の炉内状態から現在または将来における操業の異常度を推定するものである。操業異常推定装置１は、図１に示すように、ガスサンプラー２と、炉内状態推定部３と、判定処理部４と、異常度推定部５と、入力部８と、表示部９と、記憶部１１と、を備えている。
【００２１】
高炉１００は、巨大な向流移動層反応炉であり、図１に示すように、炉頂部に炉頂バンカー１０１を、炉下部に羽口１０２を備えている。炉頂バンカー１０１は、高炉１００の炉頂部から内部（炉内）に装入物１１０を装入する。炉内に装入される装入物１１０は、酸化鉄を主体とする鉄原料である鉱石（例えば鉄鉱石、焼結鉱、塊鉱石等）およびコークスである。
【００２２】
炉頂バンカー１０１は、高炉１００の炉頂部から炉下部に向かって鉱石層１１１とコークス層１１２とが交互に積層するように、装入物１１０を炉内に順次装入する。これにより、図１に示すように、鉱石層１１１とコークス層１１２との多層構造を有する炉内混合物１１３が形成される。なお、この炉内混合物１１３は、同図に示したような鉱石層１１１とコークス層１１２とからなる多層構造の他に、塊状帯、融着帯および滴下帯等を含んでいる。
【００２３】
高炉１００は、羽口１０２から吹き込まれた高温ガスを炉内ガスとして炉下部側から炉頂部側に向かって流通させることにより、炉頂バンカー１０１によって炉内に装入された各装入物１１０（鉱石およびコークス）のうち、コークスを炉内ガスによって燃焼させる。そして、高炉１００は、コークスの燃焼によって生成された一酸化炭素を含む還元性ガスを利用して、炉内の各装入物１１０のうちの鉱石に含まれる酸化鉄を還元し、銑鉄１１５を製造する。この銑鉄１１５は、図１に示すように、溶融した状態の溶銑となって高炉１００の炉下部に流れ、その後、高炉１００の外部に送出される。
【００２４】
高炉１００には、図１に示すように、高炉１００の炉体温度を計測する温度計６と、高炉１００内のガス圧力を計測する圧力計７とが設けられている。温度計６は、高炉１００の壁部に配置されており、操業中の高炉１００の炉本体の温度である炉体温度を連続的または断続的に測定する。
【００２５】
圧力計７は、操業中の高炉１００におけるガス圧力を連続的または断続的に測定する。本実施形態における圧力計７は、例えば操業中の高炉１００のガス圧力として、高炉１００の炉内に供給される高温ガス（熱風）の圧力、高炉１００の炉内ガスの圧力等を測定する。
【００２６】
入力部８は、キーボードやマウス等の入力デバイスによって構成され、作業者（オペレータ）の入力操作に従って各種情報を操業異常推定装置１に入力可能となっている。
【００２７】
表示部９は、液晶ディスプレイ等の表示デバイスによって構成され、各種データを表示する。この表示部９によって表示されるデータとしては、例えば操業中の高炉１００の炉内状態の推定データ、高炉１００の現在または将来における操業の異常度（以下、「操業異常度」という）の推定データ、高炉１００の各種測定データ、高炉１００の操業実績データ等が挙げられる。
【００２８】
記憶部１１は、例えば不揮発性のメモリ等を有する記憶装置である。記憶部１１には、後記する入力推定逆モデル１１ａ、高炉モデル１１ｂおよび高炉１００の操業に係る実績データ１１ｃが保存される。なお、入力推定逆モデル１１ａ、高炉モデル１１ｂおよび高炉１００の操業に係る実績データ１１ｃは、ハードディスク等に一旦保存し、これを読み込んでもよい。
【００２９】
ガスサンプラー２は、操業管理対象である高炉１００の炉内を流れる高温ガス（以下、「炉内ガス」という）を計測する。ガスサンプラー２は、図１に示すように、高さ方向における高炉１００の上部に配置されており、例えば予め定められた時間間隔で炉内ガスのサンプリングを行う。なお、炉内ガスのサンプリングは、予め定められた時間間隔だけではなく、必要に応じて適宜実施してもよい。また、ガスサンプラー２は、細長い筒状のゾンデ２ａを有している。ガスサンプラー２のゾンデ２ａは、操業中の高炉１００の炉壁から中心部に向けてバッチ的に装入される。
【００３０】
ガスサンプラー２は、具体的には、炉内ガスの温度や、炉内ガスの成分比を計測する。炉内ガスの成分には、一酸化炭素ＣＯ、二酸化炭素ＣＯ２、水素Ｈ２、窒素Ｎ２等が含まれる。また、ガスサンプラー２は、炉内ガスの成分比として、炉内ガス中の各ガス成分の濃度［体積％］、具体的には、一酸化炭素のガス成分濃度［ＣＯ］、二酸化炭素のガス成分濃度［ＣＯ２］、水素のガス成分濃度［Ｈ２］、および窒素のガス成分濃度［Ｎ２］等を測定する。
【００３１】
ガスサンプラー２は、高炉１００の炉内の径方向に沿った複数の位置において炉内ガスをサンプリングし、炉内ガスの成分比を計測する。つまり、本実施形態におけるガスサンプラー２は、高炉１００の炉内の径方向に沿ったガス成分の分布を測定するガス測定部として機能する。また、ガスサンプラー２は、高炉１００の炉内の径方向に沿ったガス成分の分布を計測するガス測定ステップを実行する。
【００３２】
炉内状態推定部３は、操業中の高炉１００の炉内状態を推定する炉内状態推定ステップを実行する。炉内状態推定部３は、まず操業中の高炉１００の炉内における装入物１１０の装入状態を推定する。装入物１１０の装入状態とは、具体的には鉱石層１１１とコークス層１１２との層厚分布（以下、単に「層厚分布」という）のことを示している。そして、炉内状態推定部３は、層厚分布の推定結果に基づいて、操業中の高炉１００の炉内状態を推定する。なお、炉内状態推定部３は、推定処理に先立ち、推定処理に必要なパラメータを設定するための情報（高炉１００の諸元等）等をデータベースから読み込む。以下、層厚分布の推定処理と、炉内状態の推定処理の詳細について説明する。
【００３３】
（層厚分布の推定）
炉内状態推定部３は、ガスサンプラー２によって計測されたガス成分の分布から、予め定められた対応関係に基づいて、炉内の径方向に沿った鉱石層１１１とコークス層１１２との層厚分布を推定する。
【００３４】
ここで、層厚分布とは、図２に示すように、高炉１００の径方向に沿った層厚比の分布のことを示している。また、同図において、横軸は高炉１００における径方向位置を、層厚比は全体に対する鉱石層１１１の厚さの割合をそれぞれ示している。また、中心軸線ＣＰとは、径方向位置の原点（中心）のことを示しており、径方向位置の座標値が大きくなるに従って高炉１００の炉壁に近づく。
【００３５】
層厚比は、具体的には図２に示すように、鉱石層１１１の厚さＬＯとコークス層１１２の厚さＬＣとを合わせた層厚（ＬＯ＋ＬＣ）に対する鉱石層１１１の厚さＬＯの比であり、下記式（１）で表すことができる。
【００３６】
層厚比＝ＬＯ／（ＬＯ＋ＬＣ）  ・・・（１）
【００３７】
炉内状態推定部３は、具体的には、下記式（２）で表されるガス利用率に基づいて、層厚分布を推定する。ガス利用率は、下記式（２）に示すように、炉内ガス中の一酸化炭素の濃度［ＣＯ］と二酸化炭素の濃度［ＣＯ２］の合計濃度に対する二酸化炭素の濃度［ＣＯ２］の比のことを示している。
【００３８】
ガス利用率＝［ＣＯ２］／（［ＣＯ２］＋［ＣＯ］）  ・・・（２）
【００３９】
ガス利用率は、炉内で計測される指標のうちで安定的に計測可能な指標であり、非常に重要である。ガス利用率は、ガスサンプラー２によって計測されたガス成分比から計算され、高炉１００の炉内での還元反応の進行を推定するために用いられる。図３は、ガス利用率分布の一例を示す図である。同図において、横軸は高炉１００の径方向位置を示し、縦軸はガス利用率を示している。なお、ガス利用率分布は、具体的には同図に示すように、高炉１００の径方向に沿ったガス利用率の分布のことを示している。
【００４０】
炉内状態推定部３が層厚分布を推定する際に用いる所定の対応関係は、高炉物理モデル（以下、「高炉モデル」という）に基づいて定められている。高炉モデルは、少なくとも層厚分布を含む入力データおよび高炉１００の諸元に基づいて操業中の高炉１００の炉内をシミュレートし、ガス利用率分布を推定する数学モデルである。高炉モデルとしては、例えば、非特許文献１に記載されたモデルが用いられる。
【００４１】
本実施形態における前記所定の対応関係は、以下に説明するように、ガス利用率分布の実績が高炉モデルによる計算結果と一致するように予め定められている。本実施形態では、高炉モデルに基づいて構築された、ガス利用率分布から層厚分布を推定する入力推定逆モデルが用いられる。
【００４２】
炉内状態推定部３は、入力推定逆モデルによってガス利用率分布の実績から層厚分布を推定する。つまり、炉内状態推定部３は、ガス利用率分布の実績が高炉モデルによる計算結果と一致するように、高炉モデルの入力である層厚分布を推定する。そして、炉内状態推定部３は、推定した層厚分布のデータを、その都度、記憶部１１に保存する。
【００４３】
以下、前記した入力推定逆モデルの構築方法について、図４を参照しながら説明する。同図には、高炉モデルのモデル入力である層厚分布（図４（ａ）参照）、およびモデル出力であるガス利用率分布（図４（ｂ）参照）が示されている。まず、径方向に沿った層厚分布ｘから計算される径方向に沿ったガス利用率分布ｙへの順方向の推定を考える。
【００４４】
この場合、ガス利用率分布ｙと層厚分布ｘとの関係は、下記式（３）によって表すことができる。ここで、下記式（３）において、層厚分布ｘの各要素ｘｉ（ｉ＝１，２，…，ｍ）は、径方向位置Ｐｉ（ｉ＝１，２，…，ｍ）における層厚比である。また、ガス利用率分布ｙの各要素ｙｉ（ｉ＝１，２，…，ｍ）は、径方向位置Ｐｉ（ｉ＝１，２，…，ｍ）におけるガス利用率である。また、重回帰係数ａは、入力の層厚分布ｘのうち径方向のある位置Ｐｉにおいて変化があったときに、出力であるガス利用率分布ｙに対してその変化のあった径方向位置Ｐｉを中心にどれだけ影響があったかという影響係数のような意味を有している。
【００４５】
【数１】
2018003044.tif 000003
【００４６】
従って、出力ｙから入力ｘを逆推定するためには、影響行列Ａの逆行列Ａ-1を求めればよいことになる。しかし、対象としているのは分布系であるため、隣接関係が大きいほど相関が強く、逆行列を求めることができない。そこで、本実施形態では、主成分分析を適用して入力推定逆モデルを構築する。主成分分析により出力ｙの次元を適当な次元に圧縮することが可能となり、出力ｙから入力ｘへの逆影響係数ｂを求めることが可能となる。求められた逆影響係数ｂは、入力推定逆モデルとして記憶部１１に保存される。
【００４７】
以下、前記した入力推定逆モデルの構築手順について図５〜図７を参照しながら説明する。なお、図５で示した処理フローは、例えば高炉モデルの計算処理を実行可能なコンピュータ等の演算装置においてオフラインで実行される。
【００４８】
まず、入力初期値の設定が実行される（ステップＳ１）。図６に示す入力初期値の層厚分布（以下、「デフォルト分布」という）Ｓｘｄは、層厚分布のデフォルト値である。デフォルト分布Ｓｘｄは、例えば、高炉１００における層厚分布の目標値である。入力初期値の設定は、例えば作業者の手入力によってなされても、あるいは記憶装置に保存されたデータファイルの読み込みによってなされてもよい。
【００４９】
続いて、入力セットが作成される（ステップＳ２）。入力セットＸｓは、図７に示すように、デフォルト分布Ｓｘｄと、デフォルト分布Ｓｘｄに対して各々の径方向位置Ｐｉの層厚比を所定値だけ増減させた層厚分布Ｓｘ１，Ｓｘ２，…，Ｓｘ２ｍとを合わせた層厚分布のセットのことを示している。
【００５０】
図７において、層厚分布Ｓｘ１は、位置Ｐ１の層厚比をデフォルト分布Ｓｘｄに対して所定値だけ増加側（鉱石層１１１の厚さＬＯの比率が増加する側）に振った層厚分布である。また、層厚分布Ｓｘ２は、位置Ｐ１の層厚比をデフォルト分布Ｓｘｄに対して減少側（鉱石層１１１の厚さＬＯの比率が減少する側）に所定値だけ振った層厚分布である。同様にして、位置Ｐ２，Ｐ３，…，Ｐｍについて、デフォルト分布Ｓｘｄに対して層厚比を増減させた層厚分布Ｓｘ３，Ｓｘ４，…，Ｓｘ２ｍが作成される。
【００５１】
続いて、出力セットが作成される（ステップＳ３）。出力セットＹｓは、入力セットＸｓの各層厚分布Ｓｘｄ，Ｓｘ１，…，Ｓｘ２ｍに対する高炉モデルの計算結果のセットのことを示している。
【００５２】
ここで、入力推定逆モデルを構築する演算装置には、高炉モデルの計算処理プログラムおよび各パラメータが記憶されている。計算処理プログラムおよび各パラメータは、高炉１００の諸元に基づいて予め定められており、高炉１００の諸元には、例えばコークス比、送風量、送風温度等が含まれている。演算装置は、１つの層厚分布Ｓｘが与えられると、高炉モデルの計算処理を行い、入力された層厚分布Ｓｘに対するガス利用率分布の計算結果Ｒｙを出力する。例えば、デフォルト分布Ｓｘｄに対しては計算結果Ｒｙｄが出力され、層厚分布Ｓｘｊ（ｊ＝１，２，…，２ｍ）に対して計算結果Ｒｙｊ（ｊ＝１，２，…，２ｍ）が出力される。
【００５３】
続いて、入出力セット間の影響係数が算出され、入力推定逆モデルが構築される（ステップＳ４）。本ステップでは、下記式（４）の逆影響行列Ｂの各要素である逆影響係数ｂが算出される。
【００５４】
ｘ＝Ｂｙ  ・・・（４）
【００５５】
逆影響係数ｂの算出に際しては、主成分分析が適用されてガス利用率分布の次元が圧縮される。また、重回帰によって、出力セットＹｓから入力セットＸｓを推定する推定精度が最適となるように逆影響係数ｂが決定される。本実施形態では、以上のような処理を経て決定された逆影響係数ｂによって入力推定逆モデルが構築される。そして、構築された入力推定逆モデルは、操業異常推定装置１の記憶部１１に保存される。入力推定逆モデルが構築および保存されると、本処理フローは終了する。
【００５６】
（炉内状態の推定）
炉内状態推定部３は、前記したように推定した高炉１００の炉内における層厚分布に基づいて、操業中の高炉１００の炉内状態を推定する。炉内状態推定部３は、炉内状態として、操業中の高炉１００の炉内に形成される融着帯の状態（例えば形状および温度分布等）を推定する。炉内状態推定部３は、具体的には図８に示すように、高炉１００の炉内の層厚分布を入力変数（入力セットＸｓ）として用いることにより、高炉モデルによって、炉内推定データ１２３を算出する。
【００５７】
ここで、炉内推定データ１２３は、操業中の高炉１００の炉内状態、具体的には炉内の装入物１１０等による炉内混合物１１３（図１参照）の形状および温度分布等の状態を等温線によって示した推定データである。そして、図８に示すように、この炉内推定データ１２３において、等温線の１２００℃〜１４００℃の温度域が融着帯推定データ１２５である。炉内状態推定部３は、得られた炉内推定データ１２３を、その都度、記憶部１１に保存する。
【００５８】
判定処理部４は、炉内状態推定部３による推定処理の入力変数が異常値であるか否かを判定する異常値判定ステップを実行する。本実施形態における判定処理部４は、推定処理の入力変数として炉内状態推定部３が算出したガス成分の実績データ、具体的には、操業中の高炉１００の炉内におけるガス利用率分布について、異常値判定処理を行う。
【００５９】
判定処理部４は、異常値判定処理の結果、ガス利用率分布が異常値ではないと判定した場合、その旨を炉内状態推定部３に通知するとともに、炉内状態推定部３に対して融着帯の状態の推定処理の実行を許可する。一方、判定処理部４は、異常値判定処理の結果、ガス利用率分布が異常値であると判定した場合、その旨を炉内状態推定部３に通知するとともに、炉内状態推定部３に対して融着帯の状態の推定処理の実行を禁止する。
【００６０】
異常度推定部５は、高炉１００の現在または将来における操業異常度を推定する異常度推定ステップを実行する。異常度推定部５は、炉内状態推定部３で推定した融着帯の状態（形状および温度分布等）を特徴付ける特徴量を用いて、統計的解析によって、高炉１００の現在または将来における操業の異常度を推定する。
【００６１】
異常度推定部５は、具体的には、融着帯の中心部の層厚（以下、「融着帯の中心部厚」という）、融着帯の炉壁部の層厚（以下、「融着帯の炉壁部厚」という）、高炉１００の径方向における融着帯の層厚の平均値（以下、「融着帯の平均厚」という）、融着帯の高さ、融着帯の傾斜のうちの少なくとも１つを説明変数とし、高炉１００の現在または将来における通気抵抗指数を目的関数として、回帰分析によって、高炉１００の現在または将来における操業異常度を推定する。すなわち、異常度推定部５は、融着帯推定データ１２５において、高炉１００の炉内の通気に関係する特徴量を選択し、回帰分析を行うことにより、操業異常度を示す通気抵抗指数を推定する。なおその際、主成分分析や部分的最小二乗法（ＰＬＳ：Partial Least Squares）によって次元圧縮を行ってもよい。
【００６２】
ここで、通気抵抗指数とは、高炉１００の炉内の通気抵抗を指数化したものであり、高炉操業の安定性を示す指標の一つである。また、融着帯の中心部厚とは、例えば図９の符号Ａで示した位置の層厚であり、当該位置の層厚が薄いほど、中心部の通気が良くなる。また、融着帯の炉壁部厚とは、例えば同図の符号Ｂで示した位置の層厚であり、当該位置の層厚が薄いほど、壁側の通気が良くなる。また、融着帯の平均厚についても、値が小さいほど、炉内の通気が良くなる。
【００６３】
融着帯の高さは、具体的には、高さの最小値、高さの最大値および高さの平均値のことであり、例えば図９の符号Ｃで示した位置が融着帯の高さの最小値である。融着帯の傾斜は、具体的には、融着帯の最大傾斜、最小傾斜および平均傾斜のことであり、例えば同図の符号Ｄで示した位置である。
【００６４】
異常度推定部５は、前記した特徴量のうち、例えば融着帯の平均厚、融着帯の中心部厚、融着帯の炉壁部厚、融着帯の傾斜、融着帯の高さの最小値を説明変数とした場合、下記式（５）により目的関数である通気抵抗指数を推定する。なお、下記式（５）におけるａ１〜ａ５は、重回帰係数である。
【００６５】
通気抵抗指数＝（ａ１×融着帯の平均厚）＋（ａ２×融着帯の中心部厚）＋（ａ３×融着帯の炉壁部厚）＋（ａ４×融着帯の傾斜）＋（ａ５×融着帯の高さの最小値）  ・・・（５）
【００６６】
[操業異常推定方法]
以下、本実施形態に係る操業異常推定装置１による操業異常推定方法について、図１０を参照しながら説明する。同図に示した処理フローは、高炉１００の操業中にオンラインで実行される。
【００６７】
まず、操業異常推定装置１は、高炉モデル計算に必要な高炉１００の諸元を含むパラメータ設定等を行う（ステップＳ１１）。
【００６８】
続いて、操業異常推定装置１は、データ収集のタイミングであるか否かを判定する（ステップＳ１２）。すなわち、本実施形態における操業異常推定装置１は、予め定められた時間間隔でガスサンプラー２に炉内ガスのサンプリング（後記するステップＳ１３）を実行させる。操業異常推定装置１は、炉内ガスのサンプリングタイミングが到来している場合（ステップＳ１２でＹｅｓ）、ステップＳ１３の処理に進む。一方、操業異常推定装置１は、炉内ガスのサンプリングタイミングが到来していない場合（ステップＳ１２でＮｏ）、ステップＳ１２の判定を繰り返す。
【００６９】
続いて、操業異常推定装置１は、データの収集および変数の定義を行う（ステップＳ１３）。この場合、操業異常推定装置１は、ガスサンプラー２に対して炉内ガスのサンプリングを指令する。これを受けて、ガスサンプラー２は高炉１００の中心軸線ＣＰに向けてゾンデ２ａを装入し、径方向の各位置Ｐｉにおいて炉内ガスをサンプリングするとともに、炉内ガスの温度を計測する。
【００７０】
また、操業異常推定装置１は、炉内状態推定部３によって変数定義を行う。炉内状態推定部３は、ガスサンプラー２によって収集された最新の実績データを取得し、当該実績データから層厚分布ｘの推定に必要な変数を得る計算、例えばガスサンプラー２で測定された一酸化炭素の濃度［ＣＯ］および二酸化炭素の濃度［ＣＯ２］からガス利用率を計算する。また、炉内状態推定部３は、炉内の圧力、炉体の温度等の操業実績データを取得する。
【００７１】
続いて、操業異常推定装置１は、判定処理部４によって異常値判定の実行を行う（ステップＳ１４）。判定処理部４は、ステップＳ１３で収集された実績データや計算された変数を炉内状態推定部３から受け取り、これらの値が異常値であるか否かを判定し、その判定結果を炉内状態推定部３に出力する。その際、判定処理部４は、例えば実績データや変数が予め定められた許容範囲から外れている場合、異常値であると判定する。
【００７２】
判定処理部４によって異常値であると判定された場合（ステップＳ１４でＹｅｓ）、ステップＳ１９の処理に進む。一方、判定処理部４によって異常値ではないと判定された場合（ステップＳ１４でＮｏ）、ステップＳ１５の処理に進む。
【００７３】
続いて、炉内状態推定部３は、モデル入力の推定を行う（ステップＳ１５）。モデル入力の推定において、炉内状態推定部３は、ステップＳ１３で収集された実績データおよび計算された変数を用いて、入力推定逆モデルによって層厚分布を推定する。その際、炉内状態推定部３は、具体的には上記式（４）のガス利用率分布ｙの各要素にステップＳ１３で計算したそれぞれの径方向位置Ｐｉのガス利用率分布の値を代入することにより、層厚分布ｘを算出する。
【００７４】
続いて、炉内状態推定部３は、高炉モデルの計算を行う（ステップＳ１６）。炉内状態推定部３は、具体的には、高炉１００の炉内の層厚分布を入力変数として用いて、高炉モデルによって高炉１００の炉内推定データ１２３（図８参照）を算出する。
【００７５】
続いて、異常度推定部５は、操業異常度の計算を行う（ステップＳ１７）。異常度推定部５は、ステップＳ１６で算出された炉内推定データ１２３において、融着帯の中心部厚、融着帯の炉壁部厚、融着帯の平均厚、融着帯の高さ、融着帯の傾斜のうちの少なくとも１つを説明変数とし、高炉１００の現在または将来における通気抵抗指数を目的関数として、回帰分析によって高炉１００の操業異常度を推定する。
【００７６】
続いて、炉内状態推定部３および異常度推定部５は、計算結果の保存を行う（ステップＳ１８）。本ステップにおいて、炉内状態推定部３は、収集した実績データ、およびステップＳ１５で算出した層厚分布ｘを記憶部１１に保存する。また、異常度推定部５は、ステップＳ１７で計算した操業異常度を記憶部１１に保存する。
【００７７】
続いて、操業異常推定装置１は、処理を終了するか否かを判定する（ステップＳ１９）。操業異常推定装置１は、例えば作業者によって処理を終了する指示がなされている場合等、本処理フローの終了条件が成立している場合（ステップＳ１９でＹｅｓ）、本処理フローを終了する。一方、操業異常推定装置１は、前記した本処理フローの終了条件が成立していない場合（ステップＳ１９でＮｏ）、ステップＳ１２に戻り、ステップＳ１２以降を繰り返す。
【００７８】
（計算結果の表示方法）
以下、本実施形態に係る操業異常推定装置１により、前記した処理フロー（図１０参照）を通じて取得された高炉１００の炉内状態や操業異常度等の計算結果を表示する方法について、図１１を参照しながら説明する。
【００７９】
まず、操業異常推定装置１は、計算結果に係るデータを表示するための設定を行う（ステップＳ２１）。この場合、操業異常推定装置１は、作業者等の入力操作に応じて、表示部９に表示する計算結果に係るデータの表示設定を行う。
【００８０】
本ステップにおいて、操業異常推定装置１は、表示部９に表示する計算結果に係るデータの表示内容、表示期間等を設定する。例えば、操業異常推定装置１は、操業中の高炉１００の炉内状態を示すオンライン推定結果を、図８で示したような融着帯推定データ１２５を含む炉内推定データ１２３として表示するか、あるいは融着帯推定データ１２５のみを表示するか等を決定することにより、表示内容の設定を行う。また、操業異常推定装置１は、高炉１００の炉内における操業異常度の推定結果を示すデータ（以下、「操業異常度推定データ」という）とともに、炉内推定データ１２３や融着帯推定データ１２５を表示するか等を決定することにより、表示内容の設定を行う。
【００８１】
続いて、操業異常推定装置１は、ステップＳ２１における表示設定に基づいて、計算結果に係るデータを表示する（ステップＳ２２）。この場合、操業異常推定装置１は、設定した表示内容に該当する計算結果に係るデータを記憶部１１から読み出し、表示部９に出力するとともに、その表示期間を指示する。これを受けて、表示部９は、指示された表示期間中、計算結果に係るデータを表示する。
【００８２】
以上のような処理を行う操業異常推定装置１によれば、高炉１００の炉内状態をオンラインで高精度に推定することができるとともに、推定した炉内状態から現在または将来における操業異常度を的確かつ早期に推定することができる。
【実施例】
【００８３】
（炉内状態の推定結果の一例）
図１２〜図１４は、本実施形態に係る操業異常推定装置１によって取得した炉内状態の推定結果の一例を示す図である。
【００８４】
これらの図において、径方向位置とは、高炉１００の内径と、炉内の中心軸線ＣＰ（図１参照）からの方向の位置との比によって表現される無次元の位置（無次元化径方向位置）のことを示している。また、実線Ｌ１，Ｌ３，Ｌ５は、ガスサンプラー２による測定データに基づいて算出されたガス利用率分布（炉内ガスのガス成分の実績データ）から、入力推定逆モデルによって算出された層厚分布の推定データを示している。また、実線Ｌ２，Ｌ４，Ｌ６は、入力推定逆モデルによって推定された層厚分布を入力変数として用い、高炉モデルによって算出されたガス利用率分布の推定データを示している。また、破線Ｌ１１，Ｌ１２，Ｌ１３は、前記したようにガスサンプラー２による測定データに基づいて算出されたガス利用率分布の実績データを示している。
【００８５】
図１２〜図１４に示すように、炉内状態のオンライン推定結果である融着帯推定データ１２５を高炉モデルによって算出する際の入力変数、すなわちガス利用率分布の推定データ（実線Ｌ２，Ｌ４，Ｌ６参照）は、測定データに基づくガス利用率分布の実績データ（破線Ｌ１１，Ｌ１２，Ｌ１３）とほぼ同様となっている。このことから、本実施形態に係る操業異常推定装置１を用いることにより、ガス利用率分布をその実績データとほぼ同様となるように推定できるということがわかる。すなわち、これらの図で示した融着帯推定データ１２５は、操業中の高炉１００の炉内における融着帯の形状等の炉内状態をリアルタイムかつ高精度に示している。
【００８６】
また、図１２〜図１４に示すような融着帯推定データ１２５を表示部９に表示することにより、高炉１００の操業中における融着帯の状態を把握することができる。例えば、図１２に示す融着帯推定データ１２５を参照することにより、高炉１００の径方向の中心部（同図の矢印部分）において融着帯の位置が高いことを確認することができる。また、図１３に示す融着帯推定データ１２５を参照することにより、高炉１００の炉壁近傍の部分（同図の矢印部分）において融着帯が下方に垂れ込む形状をなしていることを確認することができる。そして、図１４に示す融着帯推定データ１２５を参照することにより、高炉１００の炉内において融着帯の層厚の一部（同図の矢印部分）が分厚くなっていることを確認することができる。
【００８７】
（操業異常度の推定結果の一例）
図１５は、本実施形態に係る操業異常推定装置１によって、現在の融着帯の状態から、現在および将来の操業異常度を推定した結果の一例を示す図である。同図は、具体的には操業異常推定装置１によって推定した現在から２．５時間後までの３０分ごとの通気抵抗指数（推定値）と、同じ時間における実際の通気抵抗指数（実績値）との相関を解析したものであり、縦軸が時間、横軸が通気抵抗指数の推定値と通気抵抗指数の実績値との間の相関の大きさを示している。
【００８８】
図１５に示すように、本実施形態に係る操業異常推定装置１によって推定した通気抵抗指数は、１時間後の実際の通気抵抗指数と最も相関があることがわかる。これにより、本実施形態に係る操業異常推定装置１を用いることにより、現在の融着帯の状態（融着帯推定データ１２５）から、１時間後の高炉１００の操業異常度を的確に推定可能であることが確認された。
【００８９】
以上、本発明に係る操業異常推定装置および操業異常推定方法について、発明を実施するための形態および実施例により具体的に説明したが、本発明の趣旨はこれらの記載に限定されるものではなく、特許請求の範囲の記載に基づいて広く解釈されなければならない。また、これらの記載に基づいて種々変更、改変等したものも本発明の趣旨に含まれることはいうまでもない。
【００９０】
例えば、本実施形態に係る操業異常推定装置１は、炉内状態推定ステップにおいて、主成分分析および重回帰分析を用いて層厚分布を推定していたが、主成分分析および重回帰分析に代えて、部分的最小二乗法を用いて層厚分布を推定してもよい。この方法によってもほぼ同様の結果を得ることができる。
【符号の説明】
【００９１】
１  操業異常推定装置
２  ガスサンプラー
２ａ  ゾンデ
３  炉内状態推定部
４  判定処理部
５  異常度推定部
６  温度計
７  圧力計
８  入力部
９  表示部
１１  記憶部
１１ａ  入力推定逆モデル
１１ｂ  高炉モデル
１１ｃ  実績データ
１００  高炉
１０１  炉頂バンカー
１０２  羽口
１１０  装入物
１１１  鉱石層
１１２  コークス層
１１３  炉内混合物
１１５  銑鉄
１２３  炉内推定データ
１２５  融着帯推定データ
ＣＰ  中心軸線
Ｓxd  デフォルト分布
ｘ  層厚分布
Ｘｓ  入力セット
ｙ  ガス利用率分布
Ｙｓ  出力セット
【図１】
2018003044.tif 000004
【図２】
2018003044.tif 000005
【図３】
2018003044.tif 000006
【図４】
2018003044.tif 000007
【図５】
2018003044.tif 000008
【図６】
2018003044.tif 000009
【図７】
2018003044.tif 000010
【図８】
2018003044.tif 000011
【図９】
2018003044.tif 000012
【図１０】
2018003044.tif 000013
【図１１】
2018003044.tif 000014
【図１２】
2018003044.tif 000015
【図１３】
2018003044.tif 000016
【図１４】
2018003044.tif 000017
【図１５】
2018003044.tif 000018