(57)【要約】【課題】高炉数学モデルの解析精度を向上できる還元速度定数の補正係数の決定方法を提供する。【解決手段】高炉の炉頂で採取した炉頂ガスから求められるＣＯガス利用率と、第１候補値を乗じた第１還元速度定数、第２候補値を乗じた第２還元速度定数、及び第３候補値を乗じた第３還元速度定数が用いられる高炉数学モデルの解析結果から求められる炉頂におけるＣＯガス利用率と、の差である炉頂ηＣＯ差を求め、第２所定高さに設置される水平ゾンデで採取した炉内ガスから求められるＣＯガス利用率と、解析結果から求められる第２所定高さにおけるＣＯガス利用率と、の差である第２ηＣＯ差を求め、炉頂ηＣＯ差が所定範囲αの範囲内になるとともに第２ηＣＯ差が所定範囲γの範囲内になる第１候補値、第２候補値、及び第３候補値を、それぞれ、第１補正係数、第２補正係数、及び第３補正係数として決定する、ことを特徴とする補正係数の決定方法。【選択図】図１
2020176320.tif 000002
【特許請求の範囲】
【請求項１】
物質収支、運動量収支、及びエネルギー収支を考慮した高炉数学モデルに用いられる還元速度定数を補正する補正係数の決定方法であって、
前記還元速度定数は、Ｆｅ２Ｏ３をＦｅ３Ｏ４に還元する第１還元反応の第１還元速度定数、Ｆｅ３Ｏ４をＦｅＯに還元する第２還元反応の第２還元速度定数、及びＦｅＯをＦｅに還元する第３還元反応の第３還元速度定数であり、
前記補正係数は、前記第１還元速度定数の基準値に乗じられる第１補正係数、前記第２還元速度定数の基準値に乗じられる第２補正係数、及び前記第３還元速度定数の基準値に乗じられる第３補正係数であり、
高炉の炉頂で採取した炉頂ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記第１補正係数の候補である第１候補値を前記第１還元速度定数の基準値に乗じた前記第１還元速度定数、前記第２補正係数の候補である第２候補値を前記第２還元速度定数の基準値に乗じた前記第２還元速度定数、及び前記第３補正係数の候補である第３候補値を前記第３還元速度定数の基準値に乗じた前記第３還元速度定数が用いられる前記高炉数学モデルの解析結果から求められる、前記炉頂におけるＣＯガス利用率ηＣＯの計算値と、の差である炉頂ηＣＯ差を求め、
前記第２還元反応が進行する前記高炉の第２所定高さに設置される水平ゾンデで採取した炉内ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記解析結果から求められる前記第２所定高さにおけるＣＯガス利用率ηＣＯの計算値と、の差である第２ηＣＯ差を求め、
前記炉頂ηＣＯ差が所定範囲αの範囲内になるとともに前記第２ηＣＯ差が所定範囲γの範囲内になる前記第１候補値、前記第２候補値、及び前記第３候補値を、それぞれ、前記第１補正係数、前記第２補正係数、及び前記第３補正係数として決定する、
ことを特徴とする補正係数の決定方法。
【請求項２】
前記第１還元反応が進行する前記高炉の第１所定高さに設置される水平ゾンデで採取した炉内ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記解析結果から求められる前記第１所定高さにおけるＣＯガス利用率ηＣＯの計算値と、の差である第１ηＣＯ差を求めることをさらに含み、
前記炉頂ηＣＯ差が前記所定範囲αの範囲内になり、前記第１ηＣＯ差が所定範囲βの範囲内になり、前記第２ηＣＯ差が前記所定範囲γの範囲内になる前記第１候補値、前記第２候補値、及び前記第３候補値を、それぞれ、前記第１補正係数、前記第２補正係数、及び前記第３補正係数として決定することを特徴とする請求項１に記載の補正係数の決定方法。
【請求項３】
前記第３還元反応が進行する前記高炉の第３所定高さに設置される水平ゾンデで採取した炉内ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記解析結果から求められる前記第３所定高さにおけるＣＯガス利用率ηＣＯの計算値と、の差である第３ηＣＯ差を求めることをさらに含み、
前記炉頂ηＣＯ差が前記所定範囲αの範囲内になり、前記第２ηＣＯ差が前記所定範囲γの範囲内になり、前記第３ηＣＯ差が所定範囲δの範囲内になる前記第１候補値、前記第２候補値、及び前記第３候補値を、それぞれ、前記第１補正係数、前記第２補正係数、及び前記第３補正係数として決定することを特徴とする請求項１に記載の補正係数の決定方法。
【請求項４】
前記第１還元反応が進行する前記高炉の第１所定高さに設置される水平ゾンデで採取した炉内ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記解析結果から求められる前記第１所定高さにおけるＣＯガス利用率ηＣＯの計算値と、の差である第１ηＣＯ差を求め、
前記第３還元反応が進行する前記高炉の第３所定高さに設置される水平ゾンデで採取した炉内ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記解析結果から求められる前記第３所定高さにおけるＣＯガス利用率ηＣＯの計算値と、の差である第３ηＣＯ差を求めること、をさらに含み、
前記炉頂ηＣＯ差が前記所定範囲αの範囲内になり、前記第１ηＣＯ差が所定範囲βの範囲内になり、前記第２ηＣＯ差が前記所定範囲γの範囲内になり、前記第３ηＣＯ差が所定範囲δの範囲内になる前記第１候補値、前記第２候補値、及び前記第３候補値を、それぞれ、前記第１補正係数、前記第２補正係数、及び前記第３補正係数として決定することを特徴とする請求項１に記載の補正係数の決定方法。
【請求項５】
前記炉頂ηＣＯ差が前記所定範囲αの範囲外となる場合、前記第１〜第３候補値を変更し、
前記第１ηＣＯ差が前記所定範囲βの範囲外となる場合、前記第１候補値を変更し、
前記第２ηＣＯ差が前記所定範囲γの範囲外となる場合、前記第２候補値を変更し、
前記第３ηＣＯ差が前記所定範囲δの範囲外となる場合、前記第３候補値を変更する、
ことをさらに含むことを特徴とする請求項４に記載の補正係数の決定方法。
【請求項６】
物質収支、運動量収支、及びエネルギー収支を考慮した高炉数学モデルに用いられる還元速度定数を補正する補正係数の決定方法であって、
前記還元速度定数は、Ｆｅ２Ｏ３をＦｅ３Ｏ４に還元する第１還元反応の第１還元速度定数、Ｆｅ３Ｏ４をＦｅＯに還元する第２還元反応の第２還元速度定数、及びＦｅＯをＦｅに還元する第３還元反応の第３還元速度定数であり、
前記補正係数は、前記第１還元速度定数の基準値に乗じられる第１補正係数、前記第２還元速度定数の基準値に乗じられる第２補正係数、及び前記第３還元速度定数の基準値に乗じられる第３補正係数であり、
高炉の炉頂で採取した炉頂ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記第１補正係数の候補である第１候補値を前記第１還元速度定数の基準値に乗じた前記第１還元速度定数、前記第２補正係数の候補である第２候補値を前記第２還元速度定数の基準値に乗じた前記第２還元速度定数、及び前記第３補正係数の候補である第３候補値を前記第３還元速度定数の基準値に乗じた前記第３還元速度定数が用いられる前記高炉数学モデルの解析結果から求められる、前記炉頂におけるＣＯガス利用率ηＣＯの計算値と、の差である炉頂ηＣＯ差を求め、
前記第３還元反応が進行する前記高炉の第３所定高さに設置される水平ゾンデで採取した炉内ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記解析結果から求められる前記第３所定高さにおけるＣＯガス利用率ηＣＯの計算値と、の差である第３ηＣＯ差を求め、
前記炉頂ηＣＯ差が所定範囲αの範囲内になるとともに前記第３ηＣＯ差が所定範囲δの範囲内になる前記第１候補値、前記第２候補値、及び前記第３候補値を、それぞれ、前記第１補正係数、前記第２補正係数、及び前記第３補正係数として決定する、
ことを特徴とする補正係数の決定方法。
【請求項７】
前記第１還元反応が進行する前記高炉の第１所定高さに設置される水平ゾンデで採取した炉内ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記解析結果から求められる前記第１所定高さにおけるＣＯガス利用率ηＣＯの計算値と、の差である第１ηＣＯ差を求めることをさらに含み、
前記炉頂ηＣＯ差が前記所定範囲αの範囲内になり、前記第１ηＣＯ差が所定範囲βの範囲内になり、前記第３ηＣＯ差が前記所定範囲δの範囲内になる前記第１候補値、前記第２候補値、及び前記第３候補値を、それぞれ、前記第１補正係数、前記第２補正係数、及び前記第３補正係数として決定することを特徴とする請求項６に記載の補正係数の決定方法。
【請求項８】
前記第１所定高さに設置される水平ゾンデで炉半径方向における異なる複数の位置で前記炉内ガスを採取して、当該複数の位置での前記ＣＯガス利用率ηＣＯの実測値をそれぞれ求め、
前記高炉数学モデルにおいて前記第１所定高さの炉半径方向に並ぶ複数の計算格子を設定して、当該複数の計算格子での前記ＣＯガス利用率ηＣＯの計算値をそれぞれ求め、
前記第１ηＣＯ差を規定する前記実測値として、前記第１所定高さにおける前記複数のＣＯガス利用率ηＣＯの実測値のうちの最大値と最小値を加算した加算値を用いるとともに、前記第１ηＣＯ差を規定する前記計算値として、前記第１所定高さにおける前記複数のＣＯガス利用率ηＣＯの計算値のうちの最大値と最小値を加算した加算値を用いることを特徴とする請求項２，４，５，７の何れか一項に記載の補正係数の決定方法。
【請求項９】
前記第２所定高さに設置される水平ゾンデで炉半径方向における異なる複数の位置で前記炉内ガスを採取して、当該複数の位置での前記ＣＯガス利用率ηＣＯの実測値をそれぞれ求め、
前記高炉数学モデルにおいて前記第２所定高さの炉半径方向に並ぶ複数の計算格子を設定して、当該複数の計算格子での前記ＣＯガス利用率ηＣＯの計算値をそれぞれ求め、
前記第２ηＣＯ差を規定する前記実測値として、前記第２所定高さにおける前記複数のＣＯガス利用率ηＣＯの実測値のうちの最大値と最小値を加算した加算値を用いるとともに、前記第２ηＣＯ差を規定する前記計算値として、前記第２所定高さにおける前記複数のＣＯガス利用率ηＣＯの計算値のうちの最大値と最小値を加算した加算値を用いることを特徴とする請求項１〜５の何れか一項に記載の補正係数の決定方法。
【請求項１０】
前記第３所定高さに設置される水平ゾンデで炉半径方向における異なる複数の位置で前記炉内ガスを採取して、当該複数の位置での前記ＣＯガス利用率ηＣＯの実測値をそれぞれ求め、
前記高炉数学モデルにおいて前記第３所定高さの炉半径方向に並ぶ複数の計算格子を設定して、当該複数の計算格子での前記ＣＯガス利用率ηＣＯの計算値をそれぞれ求め、
前記第３ηＣＯ差を規定する前記実測値として、前記第３所定高さにおける前記複数のＣＯガス利用率ηＣＯの実測値のうちの最大値と最小値を加算した加算値を用いるとともに、前記第３ηＣＯ差を規定する前記計算値として、前記第３所定高さにおける前記複数のＣＯガス利用率ηＣＯの計算値のうちの最大値と最小値を加算した加算値を用いることを特徴とする請求項３〜７の何れか一項に記載の補正係数の決定方法。
【請求項１１】
物質収支、運動量収支、及びエネルギー収支を考慮した高炉数学モデルに用いられる還元速度定数を補正する補正係数の決定方法であって、
前記還元速度定数は、Ｆｅ２Ｏ３をＦｅＯに還元する第１−２還元反応の第１−２還元速度定数、及びＦｅＯをＦｅに還元する第３還元反応の第３還元速度定数であり、
前記補正係数は、前記第１−２還元速度定数の基準値に乗じられる第１−２補正係数、及び前記第３還元速度定数の基準値に乗じられる第３補正係数であり、
高炉の炉頂で採取した炉頂ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記第１−２補正係数の候補である第１−２候補値を前記第１−２還元速度定数の基準値に乗じた前記第１−２還元速度定数、及び前記第３補正係数の候補である第３候補値を前記第３還元速度定数の基準値に乗じた前記第３還元速度定数が用いられる前記高炉数学モデルの解析結果から求められる、前記炉頂におけるＣＯガス利用率ηＣＯの計算値と、の差である炉頂ηＣＯ差を求め、
前記第１−２還元反応が進行する前記高炉の第１−２所定高さに設置される水平ゾンデで採取した炉内ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記解析結果から求められる前記第１−２所定高さにおけるＣＯガス利用率ηＣＯの計算値と、の差である第１−２ηＣＯ差を求め、
前記炉頂ηＣＯ差が所定範囲αの範囲内になるとともに前記第１−２ηＣＯ差が所定範囲βまたは所定範囲γの範囲内になる前記第１−２候補値及び前記第３候補値を、それぞれ、前記第１−２補正係数及び前記第３補正係数として決定する、
ことを特徴とする補正係数の決定方法。
【請求項１２】
物質収支、運動量収支、及びエネルギー収支を考慮した高炉数学モデルに用いられる還元速度定数を補正する補正係数の決定方法であって、
前記還元速度定数は、Ｆｅ２Ｏ３をＦｅＯに還元する第１−２還元反応の第１−２還元速度定数、及びＦｅＯをＦｅに還元する第３還元反応の第３還元速度定数であり、
前記補正係数は、前記第１−２還元速度定数の基準値に乗じられる第１−２補正係数、及び前記第３還元速度定数の基準値に乗じられる第３補正係数であり、
高炉の炉頂で採取した炉頂ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記第１−２補正係数の候補である第１−２候補値を前記第１−２還元速度定数の基準値に乗じた前記第１−２還元速度定数、及び前記第３補正係数の候補である第３候補値を前記第３還元速度定数の基準値に乗じた前記第３還元速度定数が用いられる前記高炉数学モデルの解析結果から求められる、前記炉頂におけるＣＯガス利用率ηＣＯの計算値と、の差である炉頂ηＣＯ差を求め、
前記第３還元反応が進行する前記高炉の第３所定高さに設置される水平ゾンデで採取した炉内ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記解析結果から求められる前記第３所定高さにおけるＣＯガス利用率ηＣＯの計算値と、の差である第３ηＣＯ差を求め、
前記炉頂ηＣＯ差が所定範囲αの範囲内になるとともに前記第３ηＣＯ差が所定範囲δの範囲内になる前記第１−２候補値及び前記第３候補値を、それぞれ、前記第１−２補正係数及び前記第３補正係数として決定する、
ことを特徴とする補正係数の決定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、還元速度定数を補正するための補正係数の決定方法に関する。
【背景技術】
【０００２】
近年、地球温暖化の観点から、ＣＯ２排出量の削減が喫緊の課題となっている。高炉からのＣＯ２排出量を削減するには、高炉で使用するコークスや微粉炭等、炭素を多く含む還元材（以下、「炭素還元材」ともいう）の削減、すなわち溶銑１ｔあたりの炭素消費量（＝炭素消費原単位）の削減が不可避である。
【０００３】
そこで、近年、高炉の炉頂から排出されるガスからＣＯ２およびＨ２Ｏを除去したガス（以下、「調整ＢＦＧ」ともいう）を、炉下部に設けられる通常羽口や、シャフト部に設けられるシャフト羽口から吹き込み、これらのガスに含まれるＣＯやＨ２を炭素還元材の代替として使用することが検討されている（特許文献１）。
【０００４】
このような新たな操業の設計を行うには、高炉数学モデルによる数値シミュレーション（以下、「解析」ともいう）を行うことが必要不可欠である。高炉数学モデルとしては、高炉の内部領域を複数の小領域（以下、「計算格子」ともいう）に分割し、各計算格子において、物質収支、運動量収支及びエネルギー収支の数式を解くことにより、高炉内の温度分布、ガス濃度分布、鉄鉱石の還元率分等の炉内状態を解析する数学モデルがある（非特許文献１）。特許文献２では、高炉数学モデルにより炉内状態を解析するとともに、解析した炉内状態と実測した炉内状態とに基づき、的確な操業条件を決定している。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１６−１１３６７７号公報
【特許文献２】特開平８−２９５９１０号公報
【特許文献３】特公昭６３−２５０４３号公報
【特許文献４】特許第５４２７０８４号公報
【非特許文献】
【０００６】
【非特許文献１】Ｋｏｕｊｉ  Ｔａｋａｔａｎｉ，Ｔａｋａｎｏｂｕ  Ｉｎａｄａ，Ｙｕｔａｋａ  Ｕｊｉｓａｗａ，「Ｔｈｒｅｅ−ｄｉｍｅｎｓｉｏｎａｌ Ｄｙｎａｍｉｃ Ｓｉｍｕｌａｔｏｒ ｆｏｒ Ｂｌａｓｔ Ｆｕｒｎａｃｅ」、ＩＳＩＪ Ｉｎｔｅｒｎａｔｉｏｎａｌ、一般社団法人日本鉄鋼協会、１９９９年、Ｖｏｌ．３９、Ｎｏ．１、ｐ１５−２２
【非特許文献２】原行明、土屋勝、近藤真一、「酸化鉄ペレットの還元時における粒子内温度」、鉄と鋼、Ｖｏｌ．６０（１９７４）、Ｎｏ．９、ｐｐ．１２６１−１２７０
【非特許文献３】原行明、坂輪光弘、近藤真一、「鉄鉱石還元用シャフト炉の数学的モデル」、鉄と鋼、Ｖｏｌ．６２（１９７６）、Ｎｏ．３、ｐｐ．３１５−３２３
【発明の概要】
【発明が解決しようとする課題】
【０００７】
高炉数学モデルにより炉内状態を解析するには、酸化鉄を鉄に還元するまでの各還元反応の還元速度定数が必要となる。具体的には、Ｆｅ２Ｏ３をＦｅ３Ｏ４に還元する還元反応の還元速度定数、Ｆｅ３Ｏ４をＦｅＯに還元する還元反応の還元速度定数、及び、ＦｅＯをＦｅに還元する還元反応の還元速度定数が必要となる。これらの還元速度定数は、一般的に、ラボで酸化鉄試料の還元試験を行い、各還元段階（Ｆｅ２Ｏ３→Ｆｅ３Ｏ４，Ｆｅ３Ｏ４→ＦｅＯ，ＦｅＯ→Ｆｅ）における、時間に対する還元率の変化から求められる。
【０００８】
しかしながら、ラボでの還元試験により得られた各還元速度定数を、そのまま高炉数学モデルで使用した場合、高炉の炉内状態を正しく再現できないことがある。例えば、炉頂におけるＣＯガス利用率（ＣＯ２ガス濃度を、ＣＯ２ガス濃度とＣＯガス濃度の総和で除算した値。以下、「ηＣＯ」と表記する。）に関し、高炉数学モデルの解析結果から求められる値（以下、「計算値」ともいう）は、高炉の炉頂で採取したガスから求められる値（以下、「実測値」ともいう）よりも大きくなることが多い。このため、従来の高炉数学モデルでは、炉頂におけるＣＯガス利用率、すなわち炉頂ηＣＯの計算値と実測値とが一致するように、酸化鉄の各還元段階の還元速度定数それぞれに０〜１の補正係数を乗じて、還元速度定数を補正している。
【０００９】
上述したように、酸化鉄の還元反応は、Ｆｅ２Ｏ３→Ｆｅ３Ｏ４、Ｆｅ３Ｏ４→ＦｅＯ、ＦｅＯ→Ｆｅというように、３段階で進行するため、還元速度定数は３つあり、還元速度定数に乗じる補正係数も３つある。これに対して、炉頂ηＣＯは１つしかない。このため、３つの補正係数を一意的に決めることができなかった。そこで、従来は、炉頂ηＣＯの計算値と実測値とが一致するように、（１）３つの補正係数を同じ値にして３つの補正係数を同時に調整する方法や、（２）２つの補正係数は固定値（例えば、１．０）としておき、残りの１つの補正係数を調整する方法等が用いられてきた。しかしながら、このような方法で決定した補正係数により還元速度定数を補正しても、高炉数学モデルによる解析精度は十分でないことがあった。特に、シャフト部から還元ガス（例えば、調整ＢＦＧ）を吹き込む高炉を対象に解析を行った場合には、解析結果（炉内状態）が実際の高炉の炉内状態と大きく異なることがあった。
【００１０】
ここで、特許文献３には、消耗型垂直ゾンデで計測した計測データ（炉高方向における複数の箇所で計測した、温度、ガス組成、圧力）と、装入物の層頂部に配置される固定ゾンデによって推定される炉内状態の推定結果を比較し、両者が一致しない場合には、炉内状態推定モデル中の熱交換係数や反応速度定数を修正して、消耗型垂直ゾンデの計測結果と一致させることが記載されている。特許文献３では、装入物の荷下がりに伴って繰り出される消耗型垂直ゾンデの繰り出し長さに基づいて、計測データの炉高方向における位置を特定しているが、装入物は炉内を垂直に降下（荷下がり）するものではない。従って、消耗型垂直ゾンデにより取得した計測データの炉高方向における位置は、実際の高炉の炉高方向における位置と一致するとは限らない。このため、消耗型垂直ゾンデの計測結果と炉内状態の推定結果とを一致するように、炉内状態推定モデル中の熱交換係数や反応速度定数を修正したとしても、解析精度を十分に向上できないおそれがある。
【００１１】
また、特許文献４では、垂直水平ゾンデにより、高炉内部のＣＯガス濃度やＣＯ２ガス濃度を測定することが開示されている。しかしながら、垂直水平ゾンデも、特許文献３の消耗型垂直ゾンデと同様に、装入物の荷下がりに伴って降下するゾンデである。このため、垂直水平ゾンデで取得した計測データの炉高方向における位置は、実際の高炉の炉高方向における位置と一致するとは限らない。また、特許文献４では、測定したＣＯガス濃度及びＣＯ２ガス濃度から、Ｃ＋ＣＯ２→２ＣＯの反応式で示されるガス化反応速度定数Ｋｓを求めているだけであり、還元速度定数を補正する補正係数を求めることについて一切着目されていない。
【００１２】
本発明は、高炉数学モデルの解析精度を向上することができる還元速度定数の補正係数の決定方法を提供することを目的とする。
【課題を解決するための手段】
【００１３】
上記課題を解決するために、本発明に係る補正係数の決定方法は、（１）物質収支、運動量収支、及びエネルギー収支を考慮した高炉数学モデルに用いられる還元速度定数を補正する補正係数の決定方法であって、前記還元速度定数は、Ｆｅ２Ｏ３をＦｅ３Ｏ４に還元する第１還元反応の第１還元速度定数、Ｆｅ３Ｏ４をＦｅＯに還元する第２還元反応の第２還元速度定数、及びＦｅＯをＦｅに還元する第３還元反応の第３還元速度定数であり、前記補正係数は、前記第１還元速度定数の基準値に乗じられる第１補正係数、前記第２還元速度定数の基準値に乗じられる第２補正係数、及び前記第３還元速度定数の基準値に乗じられる第３補正係数であり、高炉の炉頂で採取した炉頂ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記第１補正係数の候補である第１候補値を前記第１還元速度定数の基準値に乗じた前記第１還元速度定数、前記第２補正係数の候補である第２候補値を前記第２還元速度定数の基準値に乗じた前記第２還元速度定数、及び前記第３補正係数の候補である第３候補値を前記第３還元速度定数の基準値に乗じた前記第３還元速度定数が用いられる前記高炉数学モデルの解析結果から求められる、前記炉頂におけるＣＯガス利用率ηＣＯの計算値と、の差である炉頂ηＣＯ差を求め、前記第２還元反応が進行する前記高炉の第２所定高さに設置される水平ゾンデで採取した炉内ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記解析結果から求められる前記第２所定高さにおけるＣＯガス利用率ηＣＯの計算値と、の差である第２ηＣＯ差を求め、前記炉頂ηＣＯ差が所定範囲αの範囲内になるとともに前記第２ηＣＯ差が所定範囲γの範囲内になる前記第１候補値、前記第２候補値、及び前記第３候補値を、それぞれ、前記第１補正係数、前記第２補正係数、及び前記第３補正係数として決定する、ことを特徴とする。
【００１４】
（２）前記第１還元反応が進行する前記高炉の第１所定高さに設置される水平ゾンデで採取した炉内ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記解析結果から求められる前記第１所定高さにおけるＣＯガス利用率ηＣＯの計算値と、の差である第１ηＣＯ差を求めることをさらに含み、前記炉頂ηＣＯ差が前記所定範囲αの範囲内になり、前記第１ηＣＯ差が所定範囲βの範囲内になり、前記第２ηＣＯ差が前記所定範囲γの範囲内になる前記第１候補値、前記第２候補値、及び前記第３候補値を、それぞれ、前記第１補正係数、前記第２補正係数、及び前記第３補正係数として決定することを特徴とする（１）に記載の補正係数の決定方法。
【００１５】
（３）前記第３還元反応が進行する前記高炉の第３所定高さに設置される水平ゾンデで採取した炉内ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記解析結果から求められる前記第３所定高さにおけるＣＯガス利用率ηＣＯの計算値と、の差である第３ηＣＯ差を求めることをさらに含み、前記炉頂ηＣＯ差が前記所定範囲αの範囲内になり、前記第２ηＣＯ差が前記所定範囲γの範囲内になり、前記第３ηＣＯ差が所定範囲δの範囲内になる前記第１候補値、前記第２候補値、及び前記第３候補値を、それぞれ、前記第１補正係数、前記第２補正係数、及び前記第３補正係数として決定することを特徴とする（１）に記載の補正係数の決定方法。
【００１６】
（４）前記第１還元反応が進行する前記高炉の第１所定高さに設置される水平ゾンデで採取した炉内ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記解析結果から求められる前記第１所定高さにおけるＣＯガス利用率ηＣＯの計算値と、の差である第１ηＣＯ差を求め、前記第３還元反応が進行する前記高炉の第３所定高さに設置される水平ゾンデで採取した炉内ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記解析結果から求められる前記第３所定高さにおけるＣＯガス利用率ηＣＯの計算値と、の差である第３ηＣＯ差を求めること、をさらに含み、前記炉頂ηＣＯ差が前記所定範囲αの範囲内になり、前記第１ηＣＯ差が所定範囲βの範囲内になり、前記第２ηＣＯ差が前記所定範囲γの範囲内になり、前記第３ηＣＯ差が所定範囲δの範囲内になる前記第１候補値、前記第２候補値、及び前記第３候補値を、それぞれ、前記第１補正係数、前記第２補正係数、及び前記第３補正係数として決定することを特徴とする（１）に記載の補正係数の決定方法。
【００１７】
（５）前記炉頂ηＣＯ差が前記所定範囲αの範囲外となる場合、前記第１〜第３候補値を変更し、前記第１ηＣＯ差が前記所定範囲βの範囲外となる場合、前記第１候補値を変更し、前記第２ηＣＯ差が前記所定範囲γの範囲外となる場合、前記第２候補値を変更し、前記第３ηＣＯ差が前記所定範囲δの範囲外となる場合、前記第３候補値を変更する、ことをさらに含むことを特徴とする（４）に記載の補正係数の決定方法。
【００１８】
（６）物質収支、運動量収支、及びエネルギー収支を考慮した高炉数学モデルに用いられる還元速度定数を補正する補正係数の決定方法であって、前記還元速度定数は、Ｆｅ２Ｏ３をＦｅ３Ｏ４に還元する第１還元反応の第１還元速度定数、Ｆｅ３Ｏ４をＦｅＯに還元する第２還元反応の第２還元速度定数、及びＦｅＯをＦｅに還元する第３還元反応の第３還元速度定数であり、前記補正係数は、前記第１還元速度定数の基準値に乗じられる第１補正係数、前記第２還元速度定数の基準値に乗じられる第２補正係数、及び前記第３還元速度定数の基準値に乗じられる第３補正係数であり、高炉の炉頂で採取した炉頂ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記第１補正係数の候補である第１候補値を前記第１還元速度定数の基準値に乗じた前記第１還元速度定数、前記第２補正係数の候補である第２候補値を前記第２還元速度定数の基準値に乗じた前記第２還元速度定数、及び前記第３補正係数の候補である第３候補値を前記第３還元速度定数の基準値に乗じた前記第３還元速度定数が用いられる前記高炉数学モデルの解析結果から求められる、前記炉頂におけるＣＯガス利用率ηＣＯの計算値と、の差である炉頂ηＣＯ差を求め、前記第３還元反応が進行する前記高炉の第３所定高さに設置される水平ゾンデで採取した炉内ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記解析結果から求められる前記第３所定高さにおけるＣＯガス利用率ηＣＯの計算値と、の差である第３ηＣＯ差を求め、前記炉頂ηＣＯ差が所定範囲αの範囲内になるとともに前記第３ηＣＯ差が所定範囲δの範囲内になる前記第１候補値、前記第２候補値、及び前記第３候補値を、それぞれ、前記第１補正係数、前記第２補正係数、及び前記第３補正係数として決定する、ことを特徴とする補正係数の決定方法。
【００１９】
（７）前記第１還元反応が進行する前記高炉の第１所定高さに設置される水平ゾンデで採取した炉内ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記解析結果から求められる前記第１所定高さにおけるＣＯガス利用率ηＣＯの計算値と、の差である第１ηＣＯ差を求めることをさらに含み、前記炉頂ηＣＯ差が前記所定範囲αの範囲内になり、前記第１ηＣＯ差が所定範囲βの範囲内になり、前記第３ηＣＯ差が前記所定範囲δの範囲内になる前記第１候補値、前記第２候補値、及び前記第３候補値を、それぞれ、前記第１補正係数、前記第２補正係数、及び前記第３補正係数として決定することを特徴とする（６）に記載の補正係数の決定方法。
【００２０】
（８）前記第１所定高さに設置される水平ゾンデで炉半径方向における異なる複数の位置で前記炉内ガスを採取して、当該複数の位置での前記ＣＯガス利用率ηＣＯの実測値をそれぞれ求め、前記高炉数学モデルにおいて前記第１所定高さの炉半径方向に並ぶ複数の計算格子を設定して、当該複数の計算格子での前記ＣＯガス利用率ηＣＯの計算値をそれぞれ求め、前記第１ηＣＯ差を規定する前記実測値として、前記第１所定高さにおける前記複数のＣＯガス利用率ηＣＯの実測値のうちの最大値と最小値を加算した加算値を用いるとともに、前記第１ηＣＯ差を規定する前記計算値として、前記第１所定高さにおける前記複数のＣＯガス利用率ηＣＯの計算値のうちの最大値と最小値を加算した加算値を用いることを特徴とする（２），（４），（５），（７）の何れか一項に記載の補正係数の決定方法。
【００２１】
（９）前記第２所定高さに設置される水平ゾンデで炉半径方向における異なる複数の位置で前記炉内ガスを採取して、当該複数の位置での前記ＣＯガス利用率ηＣＯの実測値をそれぞれ求め、前記高炉数学モデルにおいて前記第２所定高さの炉半径方向に並ぶ複数の計算格子を設定して、当該複数の計算格子での前記ＣＯガス利用率ηＣＯの計算値をそれぞれ求め、前記第２ηＣＯ差を規定する前記実測値として、前記第２所定高さにおける前記複数のＣＯガス利用率ηＣＯの実測値のうちの最大値と最小値を加算した加算値を用いるとともに、前記第２ηＣＯ差を規定する前記計算値として、前記第２所定高さにおける前記複数のＣＯガス利用率ηＣＯの計算値のうちの最大値と最小値を加算した加算値を用いることを特徴とする（１）〜（５）の何れか一項に記載の補正係数の決定方法。
【００２２】
（１０）前記第３所定高さに設置される水平ゾンデで炉半径方向における異なる複数の位置で前記炉内ガスを採取して、当該複数の位置での前記ＣＯガス利用率ηＣＯの実測値をそれぞれ求め、前記高炉数学モデルにおいて前記第３所定高さの炉半径方向に並ぶ複数の計算格子を設定して、当該複数の計算格子での前記ＣＯガス利用率ηＣＯの計算値をそれぞれ求め、前記第３ηＣＯ差を規定する前記実測値として、前記第３所定高さにおける前記複数のＣＯガス利用率ηＣＯの実測値のうちの最大値と最小値を加算した加算値を用いるとともに、前記第３ηＣＯ差を規定する前記計算値として、前記第３所定高さにおける前記複数のＣＯガス利用率ηＣＯの計算値のうちの最大値と最小値を加算した加算値を用いることを特徴とする（３）〜（７）の何れか一項に記載の補正係数の決定方法。
【００２３】
（１１）物質収支、運動量収支、及びエネルギー収支を考慮した高炉数学モデルに用いられる還元速度定数を補正する補正係数の決定方法であって、前記還元速度定数は、Ｆｅ２Ｏ３をＦｅＯに還元する第１−２還元反応の第１−２還元速度定数、及びＦｅＯをＦｅに還元する第３還元反応の第３還元速度定数であり、前記補正係数は、前記第１−２還元速度定数の基準値に乗じられる第１−２補正係数、及び前記第３還元速度定数の基準値に乗じられる第３補正係数であり、高炉の炉頂で採取した炉頂ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記第１−２補正係数の候補である第１−２候補値を前記第１−２還元速度定数の基準値に乗じた前記第１−２還元速度定数、及び前記第３補正係数の候補である第３候補値を前記第３還元速度定数の基準値に乗じた前記第３還元速度定数が用いられる前記高炉数学モデルの解析結果から求められる、前記炉頂におけるＣＯガス利用率ηＣＯの計算値と、の差である炉頂ηＣＯ差を求め、前記第１−２還元反応が進行する前記高炉の第１−２所定高さに設置される水平ゾンデで採取した炉内ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記解析結果から求められる前記第１−２所定高さにおけるＣＯガス利用率ηＣＯの計算値と、の差である第１−２ηＣＯ差を求め、前記炉頂ηＣＯ差が所定範囲αの範囲内になるとともに前記第１−２ηＣＯ差が所定範囲βまたは所定範囲γの範囲内になる前記第１−２候補値及び前記第３候補値を、それぞれ、前記第１−２補正係数及び前記第３補正係数として決定する、ことを特徴とする補正係数の決定方法。
【００２４】
（１２）物質収支、運動量収支、及びエネルギー収支を考慮した高炉数学モデルに用いられる還元速度定数を補正する補正係数の決定方法であって、前記還元速度定数は、Ｆｅ２Ｏ３をＦｅＯに還元する第１−２還元反応の第１−２還元速度定数、及びＦｅＯをＦｅに還元する第３還元反応の第３還元速度定数であり、前記補正係数は、前記第１−２還元速度定数の基準値に乗じられる第１−２補正係数、及び前記第３還元速度定数の基準値に乗じられる第３補正係数であり、高炉の炉頂で採取した炉頂ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記第１−２補正係数の候補である第１−２候補値を前記第１−２還元速度定数の基準値に乗じた前記第１−２還元速度定数、及び前記第３補正係数の候補である第３候補値を前記第３還元速度定数の基準値に乗じた前記第３還元速度定数が用いられる前記高炉数学モデルの解析結果から求められる、前記炉頂におけるＣＯガス利用率ηＣＯの計算値と、の差である炉頂ηＣＯ差を求め、前記第３還元反応が進行する前記高炉の第３所定高さに設置される水平ゾンデで採取した炉内ガスから求められるＣＯガス利用率ηＣＯの実測値と、前記解析結果から求められる前記第３所定高さにおけるＣＯガス利用率ηＣＯの計算値と、の差である第３ηＣＯ差を求め、前記炉頂ηＣＯ差が所定範囲αの範囲内になるとともに前記第３ηＣＯ差が所定範囲δの範囲内になる前記第１−２候補値及び前記第３候補値を、それぞれ、前記第１−２補正係数及び前記第３補正係数として決定する、ことを特徴とする補正係数の決定方法。
【発明の効果】
【００２５】
本発明によれば、高炉数学モデルの解析精度を向上することができる還元速度定数の補正係数の決定方法を提供することができる。
【図面の簡単な説明】
【００２６】
【図１】補正係数を決定する手順を示すフローチャートである（第１実施形態）。
【図２】補正係数を決定する手順の一部を示すフローチャートである。
【図３】補正係数を決定する手順の一部を示すフローチャートである。
【図４】補正係数を決定する手順を示すフローチャートである（第２実施形態）。
【図５】補正係数を決定する手順を示すフローチャートである（第３実施形態）。
【図６】ＣＯガス利用率の分布を示す図である。
【図７】ＣＯガス利用率の分布を示す図である（調整ＢＦＧ吹込み操業）。
【発明を実施するための形態】
【００２７】
本発明の一実施形態について説明する。
【００２８】
（第１実施形態）
本実施形態は、酸化鉄の還元反応における還元速度定数の補正係数を決定する方法である。本実施形態で決定する補正係数は、酸化鉄（Ｆｅ２Ｏ３）を鉄（Ｆｅ）に還元するまでの各還元反応の還元速度定数を補正する係数である。補正係数により補正される還元速度定数は、高炉数学モデルに用いられる定数であり、Ｆｅ２Ｏ３をＦｅ３Ｏ４に還元する第１還元反応の第１還元速度定数と、Ｆｅ３Ｏ４をＦｅＯに還元する第２還元反応の第２還元速度定数と、ＦｅＯをＦｅに還元する第３還元反応の第３還元速度定数である。つまり、本実施形態では、第１還元速度定数に乗じる第１補正係数と、第２還元速度定数に乗じる第２補正係数と、第３還元速度定数に乗じる第３補正係数を決定する。
【００２９】
第１〜第３還元速度定数が用いられる高炉数学モデルは、物質収支、運動量収支及びエネルギー収支を考慮した高炉数学モデルである。物質収支、運動量収支及びエネルギー収支を考慮した高炉数学モデルは、高炉の内部領域を分割して複数の小領域（計算格子）に区画し、各計算格子において、物質収支、運動量収支及びエネルギー収支の数式を解く数学モデル（以下、単に「高炉数学モデル」ともいう）である。高炉数学モデルにおいて、高炉の内部領域は、２次元方向（炉半径方向及び炉高方向）に分割されて小領域を形成してもよく、３次元方向（炉半径方向、炉高方向及び炉周方向）に分割されて小領域を形成していてもよい。具体的な高炉数学モデルは、例えば、非特許文献１に記載されているため詳細な説明は省略する。
【００３０】
高炉数学モデルにおいて、第１〜第３還元速度定数は、各還元速度定数に対応する還元反応の反応速度を求めるために用いられる。反応速度は、物質収支、運動量収支及びエネルギー収支の数式で用いられる。第１〜第３還元速度定数を用いて各還元反応の反応速度を求めるには、例えば、非特許文献２や非特許文献３に記載される未反応核３界面モデルが用いられる。未反応核３界面モデルは、原料である鉄鉱石が還元される際に、未反応核であるＦｅ２Ｏ３層と、Ｆｅ２Ｏ３層の周囲に存在するＦｅ３Ｏ４層と、Ｆｅ３Ｏ４層の周囲に存在するＦｅＯ層と、ＦｅＯ層の周囲に存在するＦｅ層と、Ｆｅ層の周囲に存在するガス境膜が形成されると仮定した上で、第１〜第３還元反応の各反応速度を算出するモデルである。なお、未反応核３界面モデルを利用して各反応速度を求める方法は、非特許文献２や非特許文献３に記載されているため、詳細な説明は省略する。未反応核３界面モデルで求められた反応速度は、物質収支、運動量収支及びエネルギー収支の数式に用いられる。
【００３１】
未反応核３界面モデルにおいて、第１還元速度定数には、Ｈ２による第１還元反応の還元速度定数と、ＣＯによる第１還元反応の還元速度定数とが含まれるが、これらの還元速度定数の第１補正係数としては、同じ値を用いることができる。また、第２還元速度定数には、Ｈ２による第２還元反応の還元速度定数と、ＣＯによる第２還元反応の還元速度定数とが含まれるが、これらの還元速度定数の第２補正係数としては、同じ値を用いることができる。また、第３還元速度定数には、Ｈ２による第３還元反応の還元速度定数と、ＣＯによる第３還元反応の還元速度定数とが含まれるが、これらの還元速度定数の第３補正係数としては、同じ値を用いることができる。
【００３２】
図１を用いて、本実施形態の補正係数の決定方法について具体的に説明する。図１は、補正係数を決定する手順を示すフローチャートである。図１に示すように、本実施形態の補正係数の決定方法は、ステップブロックＡ〜Ｇを含む。
【００３３】
まず、ステップブロックＡの処理について説明する。ステップブロックＡでは、ステップＳ１０１〜ステップＳ１０３の処理を行う。
【００３４】
ステップＳ１０１の処理では、高炉で採取したガスを用いて、炉頂部におけるＣＯガス利用率（以下、「炉頂ηＣＯ」ともいう）の実測値と、第１還元反応が進行する高炉の第１所定高さにおけるＣＯガス利用率（以下、「第１ηＣＯ」ともいう）の実測値と、第２還元反応が進行する高炉の第２所定高さにおけるＣＯガス利用率（以下、「第２ηＣＯ」ともいう）の実測値と、第３還元反応が進行する高炉の第３所定高さにおけるＣＯガス利用率（以下、「第３ηＣＯ」ともいう）の実測値を求める。
【００３５】
ＣＯガス利用率ηＣＯは、下記式（１）で表される高炉内の還元状況を示す指標である。ここで、ＣＯガス利用率ηＣＯの実測値とは、下記式（１）中のＣＯ２ガス濃度ＣＣＯ２及びＣＯガス濃度ＣＣＯとして、実測した数値を用いたＣＯガス利用率を指す。
【００３６】
2020176320.tif 000003
【００３７】
上記式（１）中、ηＣＯはＣＯガス利用率［％］を示し、ＣＣＯ２はＣＯ２ガス濃度［％］を指し、ＣＣＯはＣＯガス濃度［％］を指す。
【００３８】
炉頂ηＣＯの実測値は、高炉の炉頂で採取した炉頂ガス中のＣＯ２ガス濃度とＣＯガス濃度を、上記式（１）に代入することにより求めることができる。ここで、高炉の炉頂とは、上昇管、下降管及び除塵機を含む排ガス回収設備が存在する領域を指し、炉頂ガスは、例えば除塵機出口において採取される。
【００３９】
炉頂ガスの採取方法やガス濃度（ＣＯ２ガス濃度とＣＯガス濃度）の測定方法は、任意であり、公知の方法を用いることができる。例えば、炉頂ガス中のガス濃度は、ガスクロマトグラフィーなどのガス分析計により測定することができる。
【００４０】
第１ηＣＯの実測値は、第１還元反応が進行する第１所定高さに設置される水平ゾンデで炉内ガスを採取し、採取した炉内ガス中のＣＯ２ガス濃度及びＣＯガス濃度を、上記式（１）に代入することで求めることができる。ここで、第１還元反応が進行する第１所定高さとは、炉内温度が約５００℃〜７００℃となる高さを指す。例えば、シャフト部の上部に水平ゾンデを設置すれば、水平ゾンデを炉内温度が５００℃〜７００℃となる高さに設置することができる。なお、シャフト部とは、炉下部に向かって炉周が大きくなる部位を指し、シャフト部の上部とは、シャフト部を炉高方向に等間隔に３分割したときに最も上方に位置する部位を指す。通常羽口からストックラインまでの距離を高炉全体の高さＨ［ｍ］とし、通常羽口から水平ゾンデまでの距離を第１所定高さ位置ｈ［ｍ］とするとき、第１所定高さは無次元表記ｈ／Ｈ［−］で例えば０．７〜０．９である。
【００４１】
炉内ガスの採取に用いられる水平ゾンデは、炉壁から炉中心に向けて水平に設置されるゾンデである。本実施形態で用いる水平ゾンデには、炉内ガスを採取するサンプリング孔が形成されており、水平ゾンデが設置される高さに位置する炉内ガスを採取することができる。なお、水平ゾンデには、サンプリング孔の他にも、炉内温度測定するセンサーが設けられていてもよい。
【００４２】
水平ゾンデには、先端にサンプリング孔があり、通常、径方向にスライドさせることにより、炉半径方向における異なる複数の位置で炉内ガスを採取する。これにより、採取した位置毎に第１ηＣＯの実測値を求めてもよい。すなわち、炉半径方向における第１ηＣＯの分布が得られ、当該分布から、第１ηＣＯの平均値、最大値及び最小値、最大値と最小値との差（減算値）、並びに、最大値と最小値との合計（加算値）などを求めることができる。なお、採取した炉内ガス中のガス濃度の測定方法は、任意であり、公知の方法を用いることができる。
【００４３】
第２ηＣＯの実測値は、第２還元反応が進行する第２所定高さに設置される水平ゾンデにより炉内ガスを採取し、採取した炉内ガス中のＣＯ２ガス濃度及びＣＯガス濃度を上記式（１）に代入することで得ることができる。
【００４４】
第２ηＣＯの実測値を求める方法は、水平ゾンデを設置する高さを除き、第１ηＣＯの実測値を求める方法と同じである。水平ゾンデを設置する高さは、第２還元反応が進行する第２所定高さである。第２還元反応が進行する第２所定高さとは、炉内温度が約７００℃〜１０００℃となる高さを指す。例えば、シャフト部の中部に水平ゾンデを設置すれば、水平ゾンデを炉内温度が７００℃〜１０００℃となる高さに設置することができる。なお、シャフト部の中部とは、シャフト部を炉高方向に等間隔に３分割したときに中間に位置する部位を指す。通常羽口からストックラインまでの距離を高炉全体の高さＨ［ｍ］とし、通常羽口から水平ゾンデまでの距離を第２所定高さ位置ｈ［ｍ］とするとき、第２所定高さは無次元表記ｈ／Ｈ［−］で例えば０．４〜０．７である。
【００４５】
第３ηＣＯの実測値は、第３還元反応が進行する第３所定高さに設置される水平ゾンデにより炉内ガスを採取し、採取した炉内ガス中のＣＯ２ガス濃度及びＣＯガス濃度を上記式（１）に代入することで得ることができる。
【００４６】
第３ηＣＯの実測値を求める方法は、水平ゾンデを設置する高さを除き、第１ηＣＯの実測値を求める方法と同じである。水平ゾンデを設置する高さは、第３還元反応が進行する所定高さである。第３還元反応が進行する第３所定高さとは、炉内温度が９００℃〜１２００℃となる高さを指す。例えば、シャフト部の下部に水平ゾンデを設置すれば、水平ゾンデを炉内温度が９００℃〜１２００℃となる高さに設置することができる。なお、シャフト部の下部とは、シャフト部を炉高方向に等間隔に３分割したときに最も下方に位置する部位を指す。通常羽口からストックラインまでの距離を高炉全体の高さＨ［ｍ］とし、通常羽口から水平ゾンデまでの距離を第３所定高さ位置ｈ［ｍ］とするとき、第３所定高さは無次元表記ｈ／Ｈ［−］で例えば０．２〜０．４である。
【００４７】
なお、ステップＳ１０１の処理において、炉頂ηＣＯ及び第１〜第３ηＣＯの実測値は、同じ高炉で採取されたガスを用いて求められる。
【００４８】
以上説明したステップＳ１０１の処理を終了すると、ステップＳ１０２の処理に進む。
【００４９】
ステップＳ１０２の処理では、高炉数学モデルの解析結果から炉頂ηＣＯ及び第１〜第３ηＣＯの計算値を求める。ここで、計算値とは、上記式（１）中のＣＯ２ガス濃度ＣＣＯ２及びＣＯガス濃度ＣＣＯとして、高炉数学モデルの解析結果（ＣＯ２ガス濃度及びＣＯ２ガス濃度）を用いたＣＯガス利用率である。
【００５０】
ステップＳ１０２の処理において、高炉数学モデルによる解析では、第１〜第３補正係数の候補として任意の値（以下、「候補値」ともいう）が設定され、この候補値により第１〜第３還元速度定数（基準値）を補正した上で行われる。なお、以下の説明では、第１補正係数として設定した候補値を第１候補値と呼び、第２補正係数として設定した候補値を第２候補値と呼び、第３補正係数として設定した候補値を第３候補値と呼ぶ。第１〜第３候補値には、例えば０よりも大きく１以下の値が用いられる。第１〜第３候補値は、互いに異なる値であってもよく、互いに同じ値であってもよい。
【００５１】
高炉数学モデルによる解析は、ステップＳ１０１の処理でガス濃度を測定した高炉を解析対象として行われる。ガス濃度を測定した高炉を解析対象とするには、高炉数学モデルに当該高炉の操業諸元を入力条件として設定する必要がある。これらの入力条件を用いて、物質収支、運動量収支及びエネルギー収支の数式を解くことで、各計算格子におけるＣＯ２ガス濃度及びＣＯガス濃度が求められる。なお、高炉数学モデルによる解析では、ＣＯ２ガス濃度及びＣＯガス濃度の他にも、温度、鉄鉱石の還元率等が求められてもよい。また、物質収支、運動量収支及びエネルギー収支の数式を解く方法は、特に限定されるものではなく、公知の方法（例えば、特開２０１２−１７２２２２号公報に記載の方法）を用いることができる。
【００５２】
炉頂ηＣＯの計算値は、炉頂に位置する計算格子のＣＯ２ガス濃度とＣＯガス濃度を上記式（１）に代入することで求めることができる。また、第１〜第３ηＣＯの計算値は、第１〜第３所定高さに位置する計算格子のＣＯ２ガス濃度とＣＯガス濃度をそれぞれ上記式（１）に代入することで求めることができる。
【００５３】
以上説明したステップＳ１０２の処理を終了すると、ステップＳ１０３の処理に進む。
【００５４】
ステップＳ１０３の処理では、炉頂ηＣＯの実測値と計算値との差である炉頂ηＣＯ差、第１ηＣＯの実測値と計算値との差である第１ηＣＯ差、第２ηＣＯの実測値と計算値との差である第２ηＣＯ差、及び、第３ηＣＯの実測値と計算値との差である第３ηＣＯ差を求める。
【００５５】
第１〜第３ηＣＯ差は、１つの実測値と１つの計算値の差であってもよいが、炉半径方向のサンプリング位置毎に求めた複数の実測値のうちの２つ以上の実測値の合計値と、炉半径方向の計算格子毎に求めた複数の計算値のうちの２つ以上の計算値の合計値の差を用いてもよい。また、第１〜第３ηＣＯ差は、それぞれ、炉半径方向のサンプリング位置毎に求めた複数の実測値の平均値と、炉半径方向の計算格子毎に求めた複数の計算値の平均値との差であってもよい。第１〜第３ηＣＯ差は、それぞれ、炉半径方向のサンプリング位置毎に求めた複数の実測値のうちの最大値から最小値を減算した減算値と、炉半径方向の計算格子毎に求めた複数の計算値のうちの最大値から最小値を減算した減算値との差であってもよい。高炉の状況を正確に反映させる観点からは、第１〜第３ηＣＯ差は、それぞれ、炉半径方向のサンプリング位置毎に求めた複数の実測値のうちの最大値と最小値を加算した加算値と、炉半径方向の計算格子毎に求めた複数の計算値のうちの最大値と最小値を加算した加算値との差であることが好ましい。
【００５６】
炉頂ηＣＯ差及び第１〜第３ηＣＯ差は、実測値から計算値を減じた値でもよく、計算値から実測値を減じた値でもよい。本実施形態では、実測値から計算値を減じた値を用いた形態について説明する。
【００５７】
以上説明した、ステップＳ１０１〜ステップＳ１０３の処理により、ステップブロックＡの処理は終了する。なお、ステップブロックＡの処理において、ステップＳ１０１の処理とステップＳ１０２の処理の順序は限定されない。例えば、ステップＳ１０２の処理の後に、ステップＳ１０１の処理を行ってもよい。ステップブロックＡの処理を終了すると、ステップブロックＢの処理に進む。
【００５８】
ステップブロックＢでは、炉頂ηＣＯ差が、所定範囲α内にあるか否かを判断し、その判断結果に応じた処理を行う。ステップブロックＢは、ステップＳ１０４〜ステップＳ１０７の処理を含む。
【００５９】
ステップＳ１０４の処理では、炉頂ηＣＯ差が、炉頂ηＣＯ差に関する所定範囲αの上限値よりも大きいか否かを判断する。所定範囲αの上限値は、例えば、０．１％以上１％以下の範囲内で設定できる。ステップＳ１０４の処理において、炉頂ηＣＯ差が、所定範囲αの上限値よりも大きいと判断された場合（ステップＳ１０４のＹＥＳ）には、ステップＳ１０５の処理に進む。
【００６０】
ステップＳ１０５の処理では、ステップＳ１０２の処理において（高炉数学モデルでの計算において）第１〜第３還元速度定数に乗じた第１〜第３候補値を所定量増加する。ここで、本実施形態では、炉頂ηＣＯ差は、炉頂ηＣＯの実測値から計算値を減じることで求めている。このため、炉頂ηＣＯ差が所定範囲αの上限値よりも大きいということは、実測値を求めた高炉（つまり、実際に稼働している高炉（以下、「実炉」ともいう））での還元反応が、計算値を求めた高炉（つまり、高炉数学モデルで解析した高炉）での還元反応よりも進んでいることを意味する。このため、ステップＳ１０５の処理では、第１〜第３候補値を所定量増加する。
【００６１】
第１〜第３候補値を増加する所定量は、任意の量を設定できる。例えば、所定量は、０．０１以上０．１以下の範囲内で設定できる。ステップＳ１０５の処理を終了すると、ステップブロックＣの処理に進む。
【００６２】
一方、ステップＳ１０４の処理において、炉頂ηＣＯ差が、所定範囲αの上限値よりも大きくないと判断された場合（ステップＳ１０４のＮＯ）には、ステップＳ１０６の処理を行う。ステップＳ１０６の処理では、炉頂ηＣＯ差が、所定範囲αの下限値よりも小さいか否かを判断する。所定範囲αの下限値は、例えば、−１％以上−０．１％以下の範囲内で設定できる。ステップＳ１０６の処理において、炉頂ηＣＯ差が、所定範囲αの下限値よりも小さいと判断された場合（ステップＳ１０６のＹＥＳ）には、ステップＳ１０７の処理を行う。
【００６３】
ステップＳ１０７の処理では、ステップＳ１０２の処理において（高炉数学モデルでの計算において）第１〜第３還元速度定数に乗じた第１〜第３候補値を所定量減少する。本実施形態では、炉頂ηＣＯ差は、炉頂ηＣＯの実測値から計算値を減じることで求めているため、炉頂ηＣＯ差が所定範囲αの下限値よりも小さいということは、実炉での還元反応が、高炉数学モデルで解析した高炉の還元反応よりも遅れていることを意味する。このため、ステップＳ１０７の処理では、第１〜第３候補値を所定量減少する。
【００６４】
第１〜第３候補値を減少する所定量は、所定量減少させた第１〜第３候補値が０以下にならない範囲内において、任意の量を設定できる。例えば、所定量は、０．０１以上０．１以下の範囲内で設定できる。ステップＳ１０７の処理で減少する所定量は、ステップＳ１０５の処理で増加する所定量と同一であってもよく、異なっていてもよい。ステップＳ１０７の処理を終了すると、ステップブロックＣの処理に進む。
【００６５】
一方、ステップＳ１０６の処理において、炉頂ηＣＯ差が、所定範囲αの下限値よりも小さくないと判断された場合（ステップＳ１０６のＮＯ）には、ステップブロックＣの処理に進む。
【００６６】
ここで、ステップブロックＢにおいて、ステップＳ１０４とステップＳ１０６の処理の順序は限定されるものではなく、ステップＳ１０６の処理を行った後、ステップＳ１０４の処理を行ってもよい。
【００６７】
ステップブロックＣでは、第１ηＣＯ差が、第１ηＣＯ差に関する所定範囲βの範囲内にあるか否かを判断し、その判断結果に応じた処理を行う。ステップブロックＣは、ステップＳ１０８〜ステップＳ１１１の処理を含む。
【００６８】
ステップＳ１０８の処理では、第１ηＣＯ差が、所定範囲βの上限値よりも大きいか否かを判断する。所定範囲βの上限値は、例えば、０．１％以上１％以下の範囲内で設定できる。ステップＳ１０８の処理において、第１ηＣＯ差が、所定範囲βの上限値よりも大きいと判断された場合（ステップＳ１０８のＹＥＳ）には、ステップＳ１０９の処理を行う。
【００６９】
ステップＳ１０９の処理では、第１候補値を所定量増加する。ステップＳ１０９の処理において、所定量増加する第１候補値は、ステップブロックＢの処理が行われた第１候補値を指す。つまり、ステップＳ１０５やステップＳ１０７の処理において、第１候補値を所定量増減している場合には、所定量増減された第１候補値を所定量増加する。第１候補値を増加する所定量は、任意の量を設定できる。例えば、所定量は、０．０１以上０．１以下の範囲内で設定できる。ステップＳ１０９の処理を終了すると、ステップブロックＤの処理に進む。
【００７０】
ここで、第１ηＣＯは、第１還元反応（Ｆｅ２Ｏ３→Ｆｅ３Ｏ４）が進行している高さにおけるＣＯガス利用率であるため、第１ηＣＯの実測値と計算値との差である第１ηＣＯ差は、第１還元反応の進行状況の差を示す指標となる。また、本実施形態において、第１ηＣＯ差は、第１ηＣＯの実測値から計算値を減じることで求めているため、第１ηＣＯ差が所定範囲βの上限値よりも大きいということは、実炉での第１還元反応が、高炉数学モデルで解析した高炉での第１還元反応よりも進んでいることを意味する。このため、ステップＳ１０９の処理では、第１候補値を所定量増加する。
【００７１】
一方、ステップＳ１０８の処理において、第１ηＣＯ差が、所定範囲βの上限値よりも大きくないと判断された場合（ステップＳ１０８のＮＯ）には、ステップＳ１１０の処理を行う。ステップＳ１１０の処理では、第１ηＣＯ差が、所定範囲βの下限値よりも小さいか否かを判断する。所定範囲βの下限値は、例えば、−０．１％以上−０．０１％以下の範囲内で設定できる。ステップＳ１１０の処理において、第１ηＣＯ差が、所定範囲βの下限値よりも小さいと判断された場合（ステップＳ１１０のＹＥＳ）には、ステップＳ１１１の処理を行う。
【００７２】
ステップＳ１１１の処理では、第１候補値を所定量減少する。上述した通り、第１ηＣＯ差は、第１還元反応の進行状況の差を示す指標となり、第１ηＣＯの実測値から計算値を減じることで求めている。このため、第１ηＣＯ差が所定範囲βの下限値よりも小さいということは、実炉での第１還元反応が、高炉数学モデルで解析した高炉での第１還元反応よりも遅れていることを意味する。このため、ステップＳ１１１の処理では、第１還元反応に対応する第１候補値を所定量減少する。
【００７３】
ステップＳ１１１の処理において、所定量を減少する第１候補値は、ステップブロックＢの処理が行われた第１候補値を指す。第１候補値を減少する所定量は、所定量減少させた第１候補値が０以下にならない範囲内において、任意の量を設定できる。例えば、所定量は、０．０１以上０．１以下の範囲内の量を設定できる。ステップＳ１１１の処理で減少する所定量は、ステップＳ１０９の処理で増加する所定量と同一であってもよく、異なっていてもよい。ステップＳ１１１の処理を終了すると、ステップブロックＤの処理に進む。
【００７４】
一方、ステップＳ１１０の処理において、第１ηＣＯ差が、所定範囲βの下限値よりも小さくないと判断された場合（ステップＳ１１０のＮＯ）には、ステップブロックＤの処理に進む。
【００７５】
なお、ステップブロックＣにおいて、ステップＳ１０８とステップＳ１１０の処理の順序は限定されるものではなく、ステップＳ１１０の処理を行った後、ステップＳ１０８の処理を行ってもよい。
【００７６】
ステップブロックＤでは、第２ηＣＯ差が、第２ηＣＯ差に関する所定範囲γの範囲内にあるか否かを判断し、その判断結果に応じた処理を行う。ステップブロックＤは、ステップＳ１１２〜ステップＳ１１５の処理を含む。
【００７７】
ステップＳ１１２の処理では、第２ηＣＯ差が、所定範囲γの上限値よりも大きいか否かを判断する。ここで、所定範囲γの上限値は、例えば、０．１％以上１％以下の範囲内で設定できる。ステップＳ１１２の処理において、第２ηＣＯ差が、所定範囲γの上限値よりも大きいと判断された場合（ステップＳ１１２のＹＥＳ）には、ステップＳ１１３の処理を行う。
【００７８】
ステップＳ１１３の処理では、第２候補値を所定量増加する。ステップＳ１１３の処理において、所定量を増加する第２候補値は、ステップブロックＢの処理が行われた第２候補値を指す。第２候補値を増加する所定量は、任意の量を設定できる。例えば、所定量は、０．０１以上０．１以下の範囲内で設定できる。ステップＳ１１３の処理を終了すると、ステップブロックＥの処理に進む。
【００７９】
ここで、第２ηＣＯは、第２還元反応（Ｆｅ３Ｏ４→ＦｅＯ）が進行している高さにおけるＣＯガス利用率であり、第２ηＣＯ差は、第２還元反応の進行状況の差を示す指標となる。また、本実施形態では、第２ηＣＯ差は、第２ηＣＯの実測値から計算値を減じることで求めているため、第２ηＣＯ差が所定範囲γの上限値よりも大きいということは、実炉での第２還元反応が、高炉数学モデルで解析した高炉での第２還元反応よりも進んでいることを意味する。このため、ステップＳ１１３の処理では、第２候補値を所定量増加する。
【００８０】
一方、ステップＳ１１２の処理において、第２ηＣＯ差が、所定範囲γの上限値よりも大きくないと判断された場合（ステップＳ１１２のＮＯ）には、ステップＳ１１４の処理を行う。ステップＳ１１４の処理では、第２ηＣＯ差が、所定範囲γの下限値よりも小さいか否かを判断する。所定範囲γの下限値は、例えば、−１％以上−０．１％以下の範囲内で設定できる。ステップＳ１１４の処理において、第２ηＣＯ差が、所定範囲γの下限値よりも小さいと判断された場合（ステップＳ１１４のＹＥＳ）には、ステップＳ１１５の処理を行う。
【００８１】
ステップＳ１１５の処理では、第２候補値を所定量減少する。上述したとおり、第２ηＣＯ差は、第２還元反応の進行状況の差を示す指標となり、第２ηＣＯの実測値から計算値を減じることで求めている。このため、第２ηＣＯ差が所定範囲γの下限値よりも小さいということは、実炉での第２還元反応が、高炉数学モデルで解析した高炉での第２還元反応よりも遅れていることを意味する。このため、ステップＳ１１５の処理では、第２候補値を所定量減少する。
【００８２】
ステップＳ１１５の処理において、所定量減少する第２候補値は、ステップブロックＢの処理が行われた第２候補値を指す。第２候補値を減少する所定量は、所定量減少させた第２候補値が０以下にならない範囲内において、任意の量を設定できる。例えば、所定量は、０．０１以上０．１以下の範囲内の量を設定できる。ステップＳ１１５の処理で減少する所定量は、ステップＳ１１３の処理で増加する所定量と同一であってもよく、異なっていてもよい。ステップＳ１１５の処理を終了すると、ステップブロックＥの処理に進む。
【００８３】
一方、ステップＳ１１４の処理において、第２ηＣＯ差が、所定範囲γの下限値よりも小さくないと判断された場合（ステップＳ１１４のＮＯ）には、ステップブロックＥの処理に進む。
【００８４】
なお、ステップブロックＤにおいて、ステップＳ１１２とステップＳ１１４の処理の順序は限定されるものではなく、ステップＳ１１４の処理を行った後、ステップＳ１１２の処理を行ってもよい。
【００８５】
ステップブロックＥでは、第３ηＣＯ差が、第３ηＣＯ差に関する所定範囲δの範囲内にあるか否かを判断し、その判断結果に応じた処理を行う。ステップブロックＥは、ステップＳ１１６〜ステップＳ１１９の処理を含む。
【００８６】
ステップＳ１１６の処理では、第３ηＣＯ差が、所定範囲δの上限値よりも大きいか否かを判断する。ここで、所定範囲δの上限値は、例えば、０．１％以上１％以下の範囲内で設定できる。ステップＳ１１６の処理において、第３ηＣＯ差が、所定範囲δの上限値よりも大きいと判断された場合（ステップＳ１１６のＹＥＳ）には、ステップＳ１１７の処理を行う。
【００８７】
ステップＳ１１７の処理では、第３候補値を所定量増加する。ステップＳ１１７の処理において、所定量増加する第３候補値は、ステップブロックＢの処理が行われた第３候補値を指す。第３候補値を増加する所定量は、任意の量を設定できる。例えば、所定量は、０．０１以上０．１以下の範囲内の量を設定できる。ステップＳ１１７の処理を終了すると、ステップブロックＦの処理に進む。
【００８８】
ここで、第３ηＣＯは、第３還元反応（ＦｅＯ→Ｆｅ）が進行している高さにおけるＣＯガス利用率であるため、第３ηＣＯの実測値と計算値との差である第３ηＣＯ差は、第３還元反応の進行状況の差を示す指標となる。また、本実施形態において、第３ηＣＯ差は、第３ηＣＯの実測値から計算値を減じた値であるため、第３ηＣＯ差が所定範囲δの上限値よりも大きいということは、実炉での第３還元反応が、高炉数学モデルで解析した高炉での第３還元反応よりも進んでいることを意味する。このため、ステップＳ１１７の処理では、第３候補値を所定量増加する。
【００８９】
一方、ステップＳ１１６の処理において、第３ηＣＯ差が、所定範囲δの上限値よりも大きくないと判断された場合（ステップＳ１１６のＮＯ）には、ステップＳ１１８の処理に進む。ステップＳ１１８の処理では、第３ηＣＯ差が、所定範囲δの下限値よりも小さいか否かを判断する。所定範囲δの下限値は、−１％以上−０．１％以下の範囲内で設定できる。ステップＳ１１８の処理において、第３ηＣＯ差が、所定範囲δの下限値よりも小さいと判断された場合（ステップＳ１１８のＹＥＳ）には、ステップＳ１１９の処理を行う。
【００９０】
ステップＳ１１９の処理では、第３候補値を所定量減少する。上述したとおり、第３ηＣＯ差は、第３還元反応の進行状況の差を示し、第３ηＣＯの実測値から計算値を減じて求めている。このため、第３ηＣＯ差が所定範囲δの下限値よりも小さいということは、実炉での第３還元反応が、高炉数学モデルで解析した高炉での第３還元反応よりも遅れていることを意味する。このため、ステップＳ１１９の処理では、第３候補値を所定量減少する。
【００９１】
ステップＳ１１９の処理において、所定量減少する第３候補値は、ステップブロックＢの処理が行われた第３候補値を指す。第３候補値を減少する所定量は、所定量減少させた第３候補値が０以下とならない範囲内において、任意の量を設定できる。例えば、所定量は、０．０１以上０．１以下の範囲内の量を設定できる。ステップＳ１１９の処理で減少する所定量は、ステップＳ１１７の処理で増加する所定量と同一であってもよく、異なっていてもよい。ステップＳ１１９の処理を終了すると、ステップブロックＦの処理に進む。
【００９２】
一方、ステップＳ１１８の処理において、第３ηＣＯ差が、所定範囲δの下限値よりも小さくないと判断された場合（ステップＳ１１８のＮＯ）には、ステップブロックＦの処理に進む。
【００９３】
なお、ステップブロックＥにおいて、ステップＳ１１６とステップＳ１１８の処理の順序は限定されるものではなく、ステップＳ１１８の処理を行った後、ステップＳ１１６の処理を行ってもよい。
【００９４】
また、ステップブロックＢ〜Ｅの順序は任意であり、必ずしも図１に示す順序で行う必要はない。また、ステップブロックＢ〜Ｅにおいて、第１〜第３候補値を増減する所定量は、互いに同じ量であってもよく、異なる量であってもよい。また、ステップブロックＢ〜Ｅにおいて、所定範囲α〜δは、互いに同一の範囲であってもよく、異なる範囲であってもよい。
【００９５】
次に、ステップブロックＦについて説明する。ステップブロックＦでは、ステップブロックＢ〜Ｅの処理が行われた第１〜第３候補値を使用し、炉頂ηＣＯ及び第１〜第３ηＣＯの計算値を求める。そして、求めた炉頂ηＣＯ及び第１〜第３ηＣＯの計算値を使用し、炉頂ηＣＯ差及び第１〜第３ηＣＯ差を求める。ステップブロックＦは、ステップＳ１２０〜Ｓ１２１の処理を含む。
【００９６】
ステップＳ１２０の処理では、ステップブロックＢ〜Ｅの処理が行われた第１〜第３候補値を使用して炉頂ηＣＯ及び第１〜第３ηＣＯの計算値を求める。具体的には、ステップブロックＢ〜Ｅの処理において増減された第１〜第３候補値を、それぞれ第１〜第３還元速度定数（基準値）に乗じ、これらの還元速度定数を用いて高炉数学モデルによる解析を行う。そして、解析して得られたガス濃度（ＣＯ２濃度及びＣＯ濃度）から炉頂ηＣＯ及び第１〜第３ηＣＯの計算値を求める。炉頂ηＣＯ及び第１〜第３ηＣＯの計算値の求め方は、ステップＳ１０２と同じであるため、詳細な説明は省略する。ステップＳ１２０の処理を終了すると、ステップＳ１２１の処理に進む。
【００９７】
ステップＳ１２１の処理では、ステップＳ１２０の処理で求めた炉頂ηＣＯ及び第１〜第３ηＣＯの計算値を用いて、炉頂ηＣＯ差及び第１〜第３ηＣＯ差を求める。ステップＳ１２１の処理において、炉頂ηＣＯ及び第１〜第３ηＣＯの実測値には、ステップＳ１０１の処理で求めた実測値を用いることができる。ステップＳ１２１の処理を終了すると、ステップブロックＧの処理に進む。
【００９８】
ステップブロックＧでは、ステップブロックＦで求めた炉頂ηＣＯ差及び第１〜第３ηＣＯ差が、それぞれ対応する所定範囲α〜δ内にあるか否かを判断し、その判断結果に応じて所定の処理を行う。ステップブロックＧは、ステップＳ１２２〜Ｓ１２５の処理を含む。
【００９９】
ステップＳ１２２の処理では、炉頂ηＣＯ差が所定範囲αの範囲内にあるか否かを判断する。所定範囲αは、ステップブロックＢの処理で説明した上限値及び下限値を用いて設定できる。炉頂ηＣＯ差が所定範囲αの範囲内にある場合（ステップＳ１２２のＹＥＳ）、ステップＳ１２３の処理に進む。一方、炉頂ηＣＯ差が所定範囲αの範囲内にない場合（ステップＳ１２２のＮＯ）、ステップＳ１０４の処理に戻る。
【０１００】
ステップＳ１２３の処理では、第１ηＣＯ差が所定範囲βの範囲内にあるか否かを判断する。所定範囲βは、ステップブロックＣの処理で説明した上限値及び下限値を用いて設定できる。第１ηＣＯ差が所定範囲βの範囲内にある場合（ステップＳ１２３のＹＥＳ）、ステップＳ１２４の処理に進む。一方、第１ηＣＯ差が所定範囲βの範囲内にない場合（ステップＳ１２３のＮＯ）、ステップＳ１０４の処理に戻る。
【０１０１】
ステップＳ１２４の処理では、第２ηＣＯ差が所定範囲γの範囲内にあるか否かを判断する。所定範囲γは、ステップブロックＤの処理で説明した上限値及び下限値を用いて設定できる。第２ηＣＯ差が所定範囲γの範囲内にある場合（ステップＳ１２４のＹＥＳ）、ステップＳ１２５の処理に進む。一方、第２ηＣＯ差が所定範囲γの範囲内にない場合（ステップＳ１２４のＮＯ）、ステップＳ１０４の処理に戻る。
【０１０２】
ステップＳ１２５の処理では、第３ηＣＯ差が所定範囲δの範囲内にあるか否かを判断する。所定範囲δは、ステップブロックＥの処理で説明した上限値及び下限値を用いて設定できる。第３ηＣＯ差が所定範囲δの範囲内にある場合（ステップＳ１２５のＹＥＳ）、ステップブロックＢ〜Ｅの処理が行われた第１候補値（ステップＳ１２０で用いた第１候補値）を、第１補正係数として決定し、ステップブロックＢ〜Ｅの処理が行われた第２候補値（ステップＳ１２０で用いた第２候補値）を、第２補正係数として決定し、ステップブロックＢ〜Ｅの処理が行われた第３候補値（ステップＳ１２０で用いた第３候補値）を、第３補正係数として決定する。
【０１０３】
一方、ステップＳ１２５の処理において、第３ηＣＯ差が所定範囲δの範囲内にない場合（ステップＳ１２５のＮＯ）、ステップＳ１０４の処理に戻る。
【０１０４】
なお、ステップブロックＧにおいて、ステップＳ１２２〜Ｓ１２５の順序は任意であり、必ずしも図１に示す順序で行う必要はない。
【０１０５】
以上説明したステップＳ１０１〜Ｓ１２５の処理により、第１〜第３還元速度定数に乗じられる第１〜第３補正係数を決定することができる。
【０１０６】
本実施形態の補正係数の決定方法では、炉頂ηＣＯの計算値と実測値の差に基づいて、第１〜第３候補値を調整（探索）するだけでなく、第１〜第３還元反応が進行する位置毎に求めたηＣＯ（第１〜第３ηＣＯ）の計算値と実測値の差に基づいて、第１〜第３還元反応に対応する候補値をそれぞれ調整（探索）している。このため、本実施形態で決定される第１〜第３補正係数は、実炉における第１〜第３還元反応の状況が正確に反映されやすい。従って、本実施形態の補正係数の決定方法によれば、高炉数学モデルの解析精度を向上することができる。
【０１０７】
ここで、本実施形態では、第１候補値〜第３候補値を増加するステップＳ１０５の処理の後、図２に示すステップブロックＢ１の処理を行ってもよい。ステップブロックＢ１の処理では、第１候補値〜第３候補値を増加した後の炉頂ηＣＯの計算値と炉頂ηＣＯの実測値との差を求め、その差に応じて所定の処理を行う。ステップブロックＢ１の処理は、図２に示すように、ステップＳ１０５−１〜ステップＳ１０５−６の処理を含む。
【０１０８】
ステップＳ１０５−１の処理では、ステップＳ１０５の処理で増加した第１〜第３候補値を使用して炉頂ηＣＯの計算値を求める。具体的には、ステップＳ１０５の処理で増加した第１〜第３候補値を、それぞれ対応する還元速度定数に乗じ、この還元速度定数を用いた高炉数学モデルにより炉頂ηＣＯの計算値を求める。なお、炉頂ηＣＯの計算値を求める方法は、ステップＳ１０２と同じであるため、詳細な説明は省略する。ステップＳ１０５−１の処理を終了すると、ステップＳ１０５−２に進む。
【０１０９】
ステップＳ１０５−２の処理では、第１〜第３候補値を増加した後の炉頂ηＣＯの計算値と、炉頂ηＣＯの実測値との差を求める。この差は、炉頂ηＣＯの実測値から、第１〜第３候補値を増加した後の計算値（つまり、ステップＳ１０５−１で求めた炉頂ηＣＯの計算値）を減じた値でもよく、増加した後の計算値から実測値を減じた値でもよい。本実施形態では、第１〜第３候補値を増加した後の計算値から実測値を減じた値を用いた形態について説明する。ステップＳ１０５−２の処理を終了すると、ステップＳ１０５−３に進む。
【０１１０】
ステップＳ１０５−３の処理では、ステップＳ１０５−２の処理で求めた差が０よりも大きいか否かを判断する。ステップＳ１０５−２の処理で求めた差が０よりも大きいと判断された場合（ステップＳ１０５−３のＹＥＳ）、ステップＳ１０５−４の処理に進む。ステップＳ１０５−４の処理では、ステップＳ１０５の処理で所定量増加した第１候補値〜第３候補値を、所定量減少する。ステップＳ１０５−４の処理で減少する所定量は、ステップＳ１０５の処理で増加した所定量よりも小さい量とし、ステップＳ１０５−４の処理で所定量減少された第１候補値〜第３候補値が、ステップＳ１０２で用いた第１候補値〜第３候補値より大きくなるようにする。例えば、所定量は、０．０１以上０．１以下の範囲内の量を設定できる。ステップＳ１０５−４の処理を終了すると、ステップＳ１０８の処理に進む。
【０１１１】
一方、ステップＳ１０５−３の処理において、ステップＳ１０５−２の処理で求めた差が０よりも大きくないと判断された場合（ステップＳ１０５−３のＮＯ）、ステップＳ１０５−５の処理に進む。ステップＳ１０５−５の処理では、ステップＳ１０５−２の処理で求めた差が０よりも小さいか否かを判断する。ステップＳ１０５−２の処理で求めた差が０よりも小さいと判断された場合（ステップＳ１０５−５のＹＥＳ）、ステップＳ１０５−６の処理に進む。また、ステップＳ１０５−２の処理で求めた差が０よりも小さくないと判断された場合（ステップＳ１０５−５のＮＯ）、ステップＳ１０８の処理に進む。
【０１１２】
ステップＳ１０５−６の処理では、ステップＳ１０５の処理で所定量増加した第１〜第３候補値を、さらに所定量増加する。ステップＳ１０５−６の処理で増加する所定量は、任意の量を設定でき、例えば、ステップＳ１０５の処理で増加した所定量よりも小さい量とすることができる。ステップＳ１０５−６の処理を終了すると、ステップＳ１０８の処理に進む。
【０１１３】
以上説明したステップブロックＢ１の処理を行うことで、ステップブロックＢ１の処理を行わない場合と比較して、早期に第１〜第３補正係数を決定することができる。
【０１１４】
また、本実施形態では、第１〜第３候補値を所定量減少するステップＳ１０７の処理の後、図３に示すステップブロックＢ２の処理を行ってもよい。ステップブロックＢ２の処理では、第１〜第３候補値を減少した後の炉頂ηＣＯの計算値と、炉頂ηＣＯの実測値との差を求め、その差に応じて所定の処理を行う。ステップＢ２の処理は、ステップＳ１０７−１〜ステップＳ１０７−６の処理を含む。
【０１１５】
ステップＳ１０７−１の処理では、ステップＳ１０７の処理で減少した第１〜第３候補値を使用して炉頂ηＣＯの計算値を求める。具体的には、ステップＳ１０７の処理で所定量減少した第１〜第３候補値を、それぞれ対応する還元速度定数に乗じ、この還元速度定数を用いた高炉数学モデルにより、炉頂ηＣＯの計算値を求める。なお、炉頂ηＣＯの計算値を求める方法は、ステップＳ１０２と同じであるため、詳細な説明は省略する。ステップＳ１０７−１の処理を終了すると、ステップＳ１０７−２に進む。
【０１１６】
ステップＳ１０７−２の処理では、第１候補値〜第３候補値を減少した後の炉頂ηＣＯの計算値と炉頂ηＣＯの実測値との差を求める。この差は、炉頂ηＣＯの実測値から、第１〜第３候補値を減少した後の計算値（つまり、ステップＳ１０７−１で求めた炉頂ηＣＯの計算値）を減じた値でもよく、減少した後の計算値から実測値を減じた値でもよい。本実施形態では、実測値から、減少した後の計算値を減じた値を用いた形態について説明する。ステップＳ１０７−２の処理を終了すると、ステップＳ１０７−３に進む。
【０１１７】
ステップＳ１０７−３の処理では、ステップＳ１０７−２の処理で求めた差が０よりも大きいか否かを判断する。ステップＳ１０７−２の処理で求めた差が０よりも大きいと判断された場合（ステップＳ１０７−３のＹＥＳ）、ステップＳ１０７−４の処理に進む。ステップＳ１０７−４の処理では、ステップＳ１０７の処理で所定量減少した第１候補値〜第３候補値を、所定量増加する。ステップＳ１０７−４の処理で増加する所定量は、ステップＳ１０７の処理で減少した所定量よりも小さい量とし、ステップＳ１０７−４の処理で所定量増加された第１〜第３候補値が、ステップＳ１０２で用いた第１候補値〜第３候補値より小さくなるようにする。ステップＳ１０７−４の処理を終了すると、ステップＳ１０８の処理に進む。
【０１１８】
一方、ステップＳ１０７−３の処理において、ステップＳ１０７−２の処理で求めた差が０よりも大きくないと判断された場合（ステップＳ１０７−３のＮＯ）、ステップＳ１０７−５の処理に進む。ステップＳ１０７−５の処理では、ステップＳ１０７−２の処理で求めた差が０よりも小さいか否かを判断する。ステップＳ１０７−２の処理で求めた差が０よりも小さいと判断された場合（ステップＳ１０７−５のＹＥＳ）、ステップＳ１０７−６の処理に進む。また、ステップＳ１０７−２の処理で求めた差が０よりも小さくないと判断された場合（ステップＳ１０７−５のＮＯ）、ステップＳ１０８の処理に進む。
【０１１９】
ステップＳ１０７−６の処理では、ステップＳ１０７の処理で所定量減少した第１候補値〜第３候補値を、さらに所定量減少する。ステップＳ１０７−６の処理で減少する所定量は、所定量減少させた第１〜第３候補値が０以下にならない範囲内において、任意の量を設定できる。例えば、所定量は、０．０１以上０．１以下の範囲内の量を設定できる。ステップＳ１０７−６の処理を終了すると、ステップＳ１０８の処理に進む。
【０１２０】
以上説明したステップブロックＢ２の処理を行うことで、ステップブロックＢ２の処理を行わない場合と比較して、早期に第１〜第３補正係数を決定することができる。
【０１２１】
本実施形態では、炉頂ηＣＯ差に加えて、第１〜第３ηＣＯ差を求めているが、第１〜第３ηＣＯ差については、必ずしも全て求める必要はなく、少なくとも、炉頂ηＣＯ差に加えて、第２ηＣＯ差と第３ηＣＯ差のいずれか一方の差を求めればよい。第２ηＣＯや第３ηＣＯの実測値は、炉頂から炉高方向に十分離間した位置において測定され、炉頂ηＣＯが示す還元状況を補完する情報を有する。よって、炉頂ηＣＯ差に加えて、第２ηＣＯ差と第３ηＣＯ差のいずれか一方の差を求めれば、実炉で進行している還元反応が正確に反映されやすく、高炉数学モデルの解析精度を向上できる第１〜第３補正係数を決定することができる。
【０１２２】
一方、炉頂ηＣＯ差に加えて、第１ηＣＯ差のみを用いて第１〜第３補正係数を決定した場合、高炉数学モデルの解析精度は向上しにくい。第１ηＣＯ差を求めるために使用される第１ηＣＯの実測値は、第１還元反応が進行する第１所定高さ（例えば、シャフト部の上部）に設置された水平ゾンデで採取した炉内ガスから求められるが、シャフト部の上部では十分に還元反応が進行していないことがある。このため、第１ηＣＯの実測値には、実炉で進行している還元反応が正確に反映されないことがある。従って、第１〜第３の補正係数を決定するために用いる指標として、炉頂ηＣＯ差と第１ηＣＯ差を用いるだけでは、高炉数学モデルの解析精度が向上しにくい。
【０１２３】
以下、炉頂ηＣＯ及び第１〜第３ηＣＯ差の４つの差のうち、３つ以下の差を用いた実施形態について説明する。
【０１２４】
（第２実施形態）
本実施形態は、炉頂ηＣＯ差及び第１〜第３ηＣＯ差のうち、３つの差を用いた形態である。本実施形態で用いる３つの差には、少なくとも、炉頂ηＣＯ差に加え、第２ηＣＯ差と第３ηＣＯ差のいずれか一方の差が含まれていればよい。本実施形態の具体的な一例について、図４を用いて説明する。図４では、３つの差として、炉頂ηＣＯ差及び第１〜第２ηＣＯ差を用いている。
【０１２５】
図４に示すように、本実施形態は、ステップブロックＨ〜Ｍを含む。まず、ステップブロックＨについて説明する。ステップブロックＨでは、炉頂ηＣＯ差及び第１〜第２ηＣＯ差を求める。炉頂ηＣＯ差及び第１〜第２ηＣＯ差の求め方は、第１実施形態のステップブロックＡで説明しているため、詳細な説明は省略する。ここで、ステップブロックＨにおけるステップＳ２０１〜ステップＳ２０３は、ステップブロックＡにおけるステップＳ１０１〜Ｓ１０３にそれぞれ対応する。ステップブロックＨを終了すると、ステップブロックＩに進む。
【０１２６】
ステップブロックＩでは、炉頂ηＣＯ差が、所定範囲αの範囲内にあるか否かを判断し、その判断結果に応じた処理を行う。ステップブロックＩは、ステップＳ２０４〜ステップＳ２０７の処理を含む。
【０１２７】
ステップＳ２０４の処理では、炉頂ηＣＯ差が所定範囲αの上限値よりも大きいか否かを判断する。所定範囲αの上限値は、図１のステップＳ１０４で説明した、所定範囲αの上限値を用いることができる。ステップＳ２０４の処理において、炉頂ηＣＯ差が、所定範囲αの上限値よりも大きいと判断された場合（ステップＳ２０４のＹＥＳ）には、ステップＳ２０５の処理を行う。
【０１２８】
ステップＳ２０５の処理では、ステップブロックＨの処理でηＣＯ（実測値及び計算値）を求めていない還元反応（つまり、その還元反応が進行する高さでηＣＯを求めていない還元反応）に対応する候補値を所定量増加する。本実施形態では、ステップブロックＨの処理において、第３還元反応が進行する高さでηＣＯ（実測値及び計算値）を求めていないため、第３還元反応に対応する第３候補値を所定量増加する。所定量は、図１のステップＳ１０５の処理で説明した所定量を用いることができる。ステップＳ２０５の処理を終了すると、ステップブロックＪに進む。
【０１２９】
一方、ステップＳ２０４の処理において、炉頂ηＣＯ差が、所定範囲αの上限値よりも大きくないと判断された場合（ステップＳ２０４のＮＯ）には、ステップＳ２０６の処理を行う。ステップＳ２０６の処理では、炉頂ηＣＯ差が所定範囲αの下限値よりも小さいか否かを判断する。所定範囲αの下限値は、図１のステップＳ１０６で説明した所定範囲αの下限値を用いることができる。ステップＳ２０６の処理において、炉頂ηＣＯ差が、所定範囲αの下限値よりも小さいと判断された場合（ステップＳ２０６のＹＥＳ）には、ステップＳ２０７の処理を行う。
【０１３０】
ステップＳ２０７の処理では、ステップブロックＨの処理でηＣＯ（実測値及び計算値）を求めていない還元反応に対応する候補値（本実施形態では、第３候補値）を所定量減少する。所定量は、図１のステップＳ１０７の処理で説明した所定量を用いることができる。ステップＳ２０７の処理を終了すると、ステップブロックＪに進む。
【０１３１】
一方、ステップＳ２０６の処理において、炉頂ηＣＯ差が、所定範囲αの下限値よりも小さくないと判断された場合（ステップＳ２０６のＮＯ）には、ステップブロックＪの処理に進む。
【０１３２】
なお、ステップブロックＩにおいて、ステップＳ２０４とステップＳ２０６の処理の順序は限定されるものではなく、ステップＳ２０６の処理を行った後、ステップＳ２０４の処理を行ってもよい。
【０１３３】
ステップブロックＪでは、第１ηＣＯ差が、所定範囲βの範囲内にあるか否かを判断し、その判断結果に応じた処理を行う。ステップブロックＪは、第１実施形態のステップブロックＣと同じ処理であるため、詳細な説明は省略する。ここで、ステップブロックＪにおけるステップＳ２０８〜ステップＳ２１１は、ステップブロックＣにおけるステップＳ１０８〜Ｓ１１１にそれぞれ対応する。ステップブロックＪの処理を終了すると、ステップブロックＫに進む。
【０１３４】
ステップブロックＫでは、第２ηＣＯ差が、所定範囲γの範囲内にあるか否かを判断し、その判断結果に応じた処理を行う。ステップブロックＫは、第１実施形態のステップブロックＤと同じ処理であるため、詳細な説明は省略する。ここで、ステップブロックＫにおけるステップＳ２１２〜ステップＳ２１５は、ステップブロックＤにおけるステップＳ１１２〜Ｓ１１５にそれぞれ対応する。ステップブロックＫの処理を終了すると、ステップブロックＬに進む。
【０１３５】
なお、ステップブロックＩ〜Ｋの順序は任意であり、必ずしも図４に示す順序で行う必要はない。
【０１３６】
次に、ステップブロックＬについて説明する。ステップブロックＬでは、ステップブロックＩ〜Ｋの処理が行われた第１〜第３候補値を使用し、炉頂ηＣＯ及び第１〜第２ηＣＯの計算値を求める。そして、求めた炉頂ηＣＯ及び第１〜第２ηＣＯの計算値を使用し、炉頂ηＣＯ差及び第１〜第２ηＣＯ差を求める。炉頂ηＣＯ差及び第１〜第２ηＣＯ差を求め方は、第１実施形態のステップブロックＦと同じであるため、詳細な説明は省略する。ここで、ステップブロックＬにおけるステップＳ２１６及びステップＳ２１７は、ステップブロックＦにおけるステップＳ１２０及びステップＳ１２１にそれぞれ対応する。ステップブロックＬの処理を終了すると、ステップブロックＭに進む。
【０１３７】
ステップブロックＭでは、ステップブロックＬで求めた炉頂ηＣＯ差及び第１〜第２ηＣＯ差が、対応する所定範囲α，β，γ内にあるか否かを判断し、その判断結果に応じて所定の処理を行う。ステップブロックＭに含まれるステップＳ２１８〜Ｓ２２０は、それぞれ、第１実施形態のステップブロックＧに含まれるステップＳ１２２〜Ｓ１２４と同じ処理であるため、詳細な説明は省略する。ステップブロックＭの処理により、炉頂ηＣＯ差及び第１〜第２ηＣＯ差のそれぞれが対応する所定範囲α，β，γ内にあると判断された場合、ステップブロックＩ〜Ｋの処理が行われた第１〜第３候補値を、それぞれ、第１〜第３補正係数として決定する。
【０１３８】
上述したステップＳ２０１〜Ｓ２２０の処理により第１〜第３還元速度定数に乗じられる第１〜第３補正係数を決定することができる。
【０１３９】
本実施形態の補正係数の決定方法では、炉頂におけるηＣＯの計算値と実測値の差に基づいて第３候補値を調整するだけでなく、還元反応が進行する位置におけるηＣＯの計算値と実測値の差に基づいて第１〜第２候補値を調整している。そして、本実施形態では、第２候補値の調整にあたり、実炉で進行している還元反応が正確に反映されやすい第２ηＣＯの実測値が用いられている。このため、本実施形態で決定される第１〜第３候補値は、実炉における還元反応の状況が正確に反映されやすい。従って、本実施形態の補正係数の決定方法によれば、高炉数学モデルの解析精度を向上することができる。
【０１４０】
以上説明した本実施形態では、ステップブロックＩの処理において、ηＣＯを求めていない還元反応に対応する候補値（本実施形態では第３候補値）を調整している。しかしながら、ステップブロックＩの処理において、第１〜第３候補値全てを増減してもよい。このような場合であっても、高炉数学モデルの解析精度を向上することができる。
【０１４１】
また、本実施形態では、ステップブロックＪの処理において、第１候補値のみを調整しているが、第１〜第３候補値すべてを調整してもよい。このような場合であっても、高炉数学モデルの解析精度を向上することができる。
【０１４２】
また、本実施形態では、ステップブロックＫの処理において、第２候補値のみを調整しているが、第１〜第３候補値すべてを調整してもよい。このような場合であっても、高炉数学モデルの解析精度を向上することができる。なお、ステップブロックＫとＪの両方の処理において、第１〜第３候補値すべてを調整してもよい。
【０１４３】
また、本実施形態において、ηＣＯを求めていない還元反応に対応する候補値（本実施形態では第３候補値）は、その還元反応の前後の還元段階の還元反応が進行する高さで求めたηＣＯの実測値と計算値に基づいて調整してもよい。具体的には、第３還元反応に対応する第３候補値は、第２還元反応（第３還元反応の前の還元段階の還元反応）が進行する高さにおけるηＣＯ（第２ηＣＯ）の実測値と計算値の差に基づいて調整されてもよい。言い換えれば、第３候補値は、ステップブロックＫの処理において、第２候補値とともに調整されてもよい。このような場合であっても、高炉数学モデルの解析精度を向上することができる。
【０１４４】
また、本実施形態では、炉頂ηＣＯ差及び第１〜第２ηＣＯ差に着目して、第１〜第３補正係数を決定しているが、これに限るものではない。具体的には、炉頂ηＣＯ差及び第２〜第３ηＣＯ差に着目して、第１〜第３補正係数を決定することができる。この場合には、炉頂ηＣＯ差が所定範囲αの範囲内になり、第２ηＣＯ差が所定範囲γの範囲内になり、かつ第３ηＣＯ差が所定範囲δの範囲内になるときの第１〜第３候補値を第１〜第３補正係数として決定することができる。また、炉頂ηＣＯ差及び第１，第３ηＣＯ差に着目して、第１〜第３補正係数を決定することができる。この場合には、炉頂ηＣＯ差が所定範囲αの範囲内になり、第１ηＣＯ差が所定範囲βの範囲内になり、かつ第３ηＣＯ差が所定範囲δの範囲内になるときの第１〜第３候補値を第１〜第３補正係数として決定することができる。
【０１４５】
（第３実施形態）
次に、本発明の第３実施形態について説明する。本実施形態は、炉頂ηＣＯ及び第１〜第３ηＣＯ差のうち、２つの差を用いた実施形態である。本実施形態で用いる２つの差には、少なくとも、炉頂ηＣＯ差と、第２ηＣＯ差と第３ηＣＯ差のいずれか一方の差が含まれていればよい。本実施形態の具体的な一例について、図５を用いて説明する。図５では、２つの差として、炉頂ηＣＯ差と第２ηＣＯ差を用いている。
【０１４６】
本実施形態は、ステップブロックＮ〜Ｒを含む。まず、ステップブロックＮについて説明する。ステップブロックＮでは、炉頂ηＣＯ差及び第２ηＣＯ差を求める。炉頂ηＣＯ差及び第２ηＣＯ差の求め方は、第１実施形態のステップブロックＡで説明しているため、詳細な説明は省略する。ここで、ステップブロックＮにおけるステップＳ３０１〜ステップＳ３０３は、ステップブロックＡにおけるステップＳ１０１〜ステップＳ１０３にそれぞれ対応する。ステップブロックＮを終了するとステップブロックＯに進む。
【０１４７】
ステップブロックＯでは、炉頂ηＣＯ差が、所定範囲αの範囲内にあるか否かを判断し、その判断結果に応じた処理を行う。ステップブロックＯは、ステップＳ３０４〜ステップＳ３０７の処理を含む。
【０１４８】
ステップＳ３０４の処理では、炉頂ηＣＯ差が所定範囲αの上限値よりも大きいか否かを判断する。所定範囲αの上限値は、図１のステップＳ１０４で説明した所定範囲αの上限値を用いることができる。ステップＳ３０４の処理において、炉頂ηＣＯ差が、所定範囲αの上限値よりも大きいと判断された場合（ステップＳ３０４のＹＥＳ）には、ステップＳ３０５の処理を行う。
【０１４９】
ステップＳ３０５の処理では、ステップブロックＮの処理でηＣＯ（実測値及び計算値）を求めていない還元反応（つまり、その還元反応が進行する高さでηＣＯを求めていない還元反応）に対応する候補値を所定量増加する。本実施形態では、ステップブロックＮの処理において、第１及び第３還元反応が進行する高さにおいてηＣＯ（実測値及び計算値）を求めていないため、第１還元反応及び第３還元反応に対応する第１候補値及び第３候補値を所定量増加する。所定量は、図１のステップＳ１０５の処理で説明した所定量を用いることができる。ステップＳ３０５の処理を終了すると、ステップブロックＰに進む。
【０１５０】
一方、ステップＳ３０４の処理において、炉頂ηＣＯ差が、所定範囲αの上限値よりも大きくないと判断された場合（ステップＳ３０４のＮＯ）には、ステップＳ３０６の処理を行う。ステップＳ３０６の処理では、炉頂ηＣＯ差が所定範囲αの下限値よりも小さいか否かを判断する。所定範囲αの下限値は、図１のステップＳ１０６で説明した所定範囲αの下限値を用いることができる。ステップＳ３０６の処理において、炉頂ηＣＯ差が、所定範囲αの下限値よりも小さいと判断された場合（ステップＳ３０６のＹＥＳ）には、ステップＳ３０７の処理を行う。
【０１５１】
ステップＳ３０７の処理では、ステップブロックＮの処理でηＣＯ（実測値及び計算値）を求めていない還元反応に対応する候補値（本実施形態では、第１及び第３候補値）を所定量減少する。所定量は、図１のステップＳ１０７の処理で説明した所定量を用いることができる。ステップＳ３０７の処理を終了すると、ステップブロックＰに進む。
【０１５２】
一方、ステップＳ３０６の処理において、炉頂ηＣＯ差が、所定範囲αの下限値よりも小さくないと判断された場合（ステップＳ３０６のＮＯ）には、ステップブロックＰの処理に進む。
【０１５３】
なお、ステップブロックＯにおいて、ステップＳ３０４とステップＳ３０６の処理の順序は限定されるものではなく、ステップＳ３０６の処理を行った後、ステップＳ３０４の処理を行ってもよい。
【０１５４】
ステップブロックＰでは、第２ηＣＯ差が、所定範囲γの範囲内にあるか否かを判断し、その判断結果に応じた処理を行う。ステップブロックＰは、第１実施形態のステップブロックＤと同じ処理であるため、詳細な説明は省略する。ここで、ステップブロックＰにおけるステップＳ３０８〜ステップＳ３１１は、ステップブロックＤにおけるステップＳ１１２〜ステップＳ１１５にそれぞれ対応する。ステップブロックＰの処理を終了すると、ステップブロックＱに進む。
【０１５５】
なお、ステップブロックＯ〜Ｐの順序は任意であり、必ずしも図５に示す順序で行う必要はない。
【０１５６】
次に、ステップブロックＱについて説明する。ステップブロックＱでは、ステップブロックＯ〜Ｐの処理が行われた第１〜第３候補値を使用し、炉頂ηＣＯ及び第２ηＣＯの計算値を求める。そして、求めた炉頂ηＣＯ及び第２ηＣＯの計算値を使用し、炉頂ηＣＯ差及び第２ηＣＯ差を求める。炉頂ηＣＯ差及び第２ηＣＯ差を求め方は、第１実施形態のステップブロックＦと同じであるため、詳細な説明は省略する。ここで、ステップブロックＱにおけるステップＳ３１２及びステップＳ３１３は、ステップブロックＦにおけるステップＳ１２０〜ステップＳ１２１にそれぞれ対応する。ステップブロックＱの処理を終了すると、ステップブロックＲに進む。
【０１５７】
ステップブロックＲでは、ステップブロックＱで求めた炉頂ηＣＯ差及び第２ηＣＯ差について、対応する所定範囲α，γ内にあるか否かを判断し、その判断結果に応じて所定の処理を行う。ステップブロックＲに含まれるステップＳ３１４及びＳ３１５は、第１実施形態のステップブロックＧに含まれるステップＳ１２２及びＳ１２４と同じ処理であるため、詳細な説明は省略する。ステップブロックＲの処理により、炉頂ηＣＯ差及び第２ηＣＯ差のそれぞれが対応する所定範囲α，γ内にあると判断された場合、ステップブロックＯ〜Ｐの処理が行われた第１〜第３候補値を、それぞれ、第１〜第３補正係数として決定する。
【０１５８】
以上説明したステップＳ３０１〜Ｓ３１５の処理により、第１〜第３還元速度定数に乗じられる第１〜第３補正係数を決定することができる。
【０１５９】
本実施形態の補正係数の決定方法では、炉頂におけるηＣＯの計算値と実測値の差に基づいて第１及び第３候補値を調整するだけでなく、還元反応が進行する位置におけるηＣＯの計算値と実測値の差に基づいて第２候補値を調整している。そして、本実施形態では、第２候補値の調整にあたり、実炉で進行している還元反応が正確に反映されやすい第２ηＣＯの実測値が用いられている。このため、本実施形態で決定される第１〜第３補正係数は、実炉における還元反応の状況が正確に反映されやすい。従って、本実施形態の補正係数の決定方法によれば、高炉数学モデルの解析精度を向上することができる。
【０１６０】
なお、本実施形態では、ステップブロックＯの処理において、ηＣＯ（実測値及び計算値）を求めていない還元反応に対応する候補値（本実施形態では第１候補値及び第３候補値）のみを所定量増加しているが、第１〜第３候補値全てを調整してもよい。この場合であっても、高炉数学モデルの解析精度を向上することができる。
【０１６１】
また、本実施形態では、ステップブロックＰの処理において、第２候補値のみを調整しているが、第１〜第３候補値全てを調整してもよい。この場合であっても、高炉数学モデルの解析精度を向上することができる。
【０１６２】
また、本実施形態では、炉頂ηＣＯ差及び第２ηＣＯ差に着目して、第１〜第３補正係数を決定しているが、これに限るものではない。具体的には、炉頂ηＣＯ差及び第３ηＣＯ差に着目して、第１〜第３補正係数を決定することができる。この場合には、炉頂ηＣＯ差が所定範囲αの範囲内になり、かつ第３ηＣＯ差が所定範囲δの範囲内になるときの第１〜第３候補値を第１〜第３補正係数として決定することができる。
【０１６３】
以上、３つの実施形態を説明したが、本発明は、同一段階の還元反応が進行する高さ位置に複数の水平ゾンデが設置されている場合にも適用できる。例えば、第２還元反応が進行するシャフト部中部の第２Ａ所定高さと第２Ｂ所定高さとに１本ずつ、計２本の水平ゾンデが設置されている場合、各高さ位置について第２ＡηＣＯ差及び第２ＢηＣＯ差を求め、その平均値を第２ηＣＯ差として、上述の各処理を行えばよい。このことは、第１還元反応が進行するシャフト部上部及び第３還元反応が進行するシャフト部下部に、複数の水平ゾンデが設置されている場合も同様である。水平ゾンデの本数は基本的には多いほど好ましく、各段階の還元反応が進行する高さ位置に分散して設置されていることが好ましい。
【０１６４】
また、以上の説明では未反応核３界面モデルによる解析を例に挙げたが、本発明は、未反応核２界面モデルを用いる場合にも適用できる。未反応核２界面モデルにおいては、例えば、Ｆｅ２Ｏ３からＦｅＯまでの還元が１段階で進むものとして扱う。Ｆｅ２Ｏ３→ＦｅＯの還元反応を第１−２還元反応とし、このときの還元速度定数を第１−２還元速度定数とするとともに、ＦｅＯ→Ｆｅの還元反応を第３還元反応とし、このときの還元速度定数を第３還元速度定数とすると、求めるべき補正係数は２つになる。したがって、炉頂ηＣＯに加えて水平ゾンデで測定されたηＣＯ（第１−２ηＣＯまたは第３ηＣＯ）を用いることにより、第１−２補正係数及び第３補正係数を決定して、高炉数学モデルの解析精度を向上することができる。ここで、第１−２補正係数は、第１−２還元速度定数（基準値）に乗じられる値であり、第３補正係数は、第３還元速度定数（基準値）に乗じられる値である。
【０１６５】
具体的には、炉頂ηＣＯ差及び第１−２ηＣＯ差を求め、炉頂ηＣＯ差が所定範囲αの範囲内になるとともに第１−２ηＣＯ差が所定範囲β（またはγ）の範囲内になる第１−２候補値及び第３候補値を、それぞれ、第１−２補正係数及び第３補正係数として決定することができる。また、炉頂ηＣＯ差及び第３ηＣＯ差を求め、炉頂ηＣＯ差が所定範囲αの範囲内になるとともに第３ηＣＯ差が所定範囲δの範囲内になる第１−２候補値及び第３候補値を、それぞれ、第１−２補正係数及び第３補正係数として決定することができる。
【実施例】
【０１６６】
次に、実施例を挙げて本発明をより具体的に説明する。ただし、本発明はこれらの実施例のみに限定されるものではない。
【０１６７】
試験高炉を用いて、下記表１に示す操業諸元で、所定期間操業を行った。この試験高炉のシャフト部上部、シャフト部中部、およびシャフト部下部には、それぞれ水平ゾンデが設置されていた。
【０１６８】
【表１】
2020176320.tif 000004
【０１６９】
所定期間が経過後、試験高炉の炉頂から排出される排ガスを炉頂において採取し、排ガス中のＣＯ２濃度及びＣＯ濃度を測定した。測定したＣＯ２濃度及びＣＯ濃度を上記式（１）に代入し、炉頂ηＣＯの実測値を求めた。また、試験高炉に設置される水平ゾンデを用いて、炉半径方向における異なる複数の位置で炉内ガスを採取した。採取した炉内ガス毎に、ＣＯ２濃度及びＣＯ濃度を測定し、測定したＣＯ２濃度及びＣＯ濃度を上記式（１）に代入して第１〜第３ηＣＯの実測値を求めた。
【０１７０】
試験高炉を解析対象として高炉数学モデルを実施した。高炉数学モデルには、第１〜第３候補値をそれぞれ乗じた第１〜第３還元速度定数を用いた。また、高炉数学モデルでは、試験高炉の操業諸元を入力条件として設定した。高炉数学モデルには、非特許文献１に記載の高炉数学モデルを使用し、各計算格子において、物質収支、運動量収支及びエネルギー収支の数式を解くことにより、各計算格子におけるＣＯ２濃度及びＣＯ濃度を求めた。
【０１７１】
炉頂部、シャフト部上部、シャフト部中部、およびシャフト部下部に位置する各計算格子の解析結果（ＣＯ２濃度及びＣＯ濃度）を、それぞれ上記式（１）に代入し、炉頂ηＣＯ及び第１〜第３ηＣＯの計算値を求めた。
【０１７２】
炉頂ηＣＯ及び第１〜第３ηＣＯのそれぞれについて、計算値と実測値の差を求め、炉頂ηＣＯ差及び第１〜第３ηＣＯ差を求めた。なお、第１〜第３ηＣＯ差については、実測値のうちの最大値と最小値を加算した加算値と、計算値のうちの最大値と最小値を加算した加算値との差とした。
【０１７３】
炉頂ηＣＯ差及び第１〜第３ηＣＯ差のそれぞれが所定範囲α〜所定範囲δになるまで、図１に示すステップブロックＢ〜Ｇの処理を繰り返し、第１〜第３補正係数を決定した。
【０１７４】
決定した第１〜第３補正係数をそれぞれ乗じた第１〜第３還元速度定数を使用して、上述した高炉数学モデルを実施し、高炉内におけるＣＯガス利用率の分布及び温度分布を得た。また、試験高炉に設置される水平ゾンデを用いて、炉内ガスを採取してＣＯガス利用率の分布を得た。ＣＯガス利用率の分布を図６に示す。なお、図６において、縦軸は、炉底部に設けられる通常羽口からの炉高方向における距離［ｍ］を示し、横軸はＣＯガス利用率［％］を示す。
【０１７５】
図６に示すように、ＣＯガス利用率の分布に関し、計算結果（高炉数学モデルで解析した結果）と測定結果（▲のプロット）は良く一致していた。
【０１７６】
次に、試験高炉において、炉頂から排出される排ガス（ＢＦＧ）からＨ２Ｏ及びＣＯ２を除去したガス（調整ＢＦＧ）をシャフト部に設けられるシャフト羽口から吹き込むとともに、コークス炉から排出されるコークス炉ガス（ＣＯＧ）を通常羽口に吹き込む操業（以下「調整ＢＦＧ吹込み操業」ともいう）に変更した。具体的な操業諸元を、下記表２に示す。所定期間が経過後、水平ゾンデを用いて、炉内ガスを採取し、ＣＯガス利用率の分布を得た。
【０１７７】
【表２】
2020176320.tif 000005
【０１７８】
一方、調整ＢＦＧ吹込み操業に変更した試験高炉を解析対象として、高炉数学モデルを実施し、定常状態に達するまで計算を行った。高炉数学モデルには、ベース操業に合わせて決定された第１〜第３補正係数をそれぞれ乗じた第１〜第３還元速度定数を用いた。また、高炉数学モデルの入力条件には、調整ＢＦＧ吹込み操業に変更した試験高炉の条件（表２）を用いた。定常状態に達した後、ＣＯガス利用率の分布を得た。
【０１７９】
ＣＯガス利用率の分布を図７に示す。なお、図７において、縦軸は、通常羽口からの炉高方向における距離［ｍ］を示し、横軸はＣＯガス利用率［％］を示す。
【０１８０】
図７に示すように、ＣＯガス利用率の分布に関し、計算結果（高炉数学モデルで解析した結果）と測定結果（▲のプロット）は良く一致していた。
【図１】
2020176320.tif 000006
【図２】
2020176320.tif 000007
【図３】
2020176320.tif 000008
【図４】
2020176320.tif 000009
【図５】
2020176320.tif 000010
【図６】
2020176320.tif 000011
【図７】
2020176320.tif 000012