(57)【特許請求の範囲】
【請求項１】
無端帯状になるように連結されて周回する複数のパレットのうち、原料給鉱部に位置するパレットに順次焼結原料を給鉱し、当該パレットが移動する過程で当該焼結原料の焼結を行って焼結鉱を生産し、当該パレットが焼結鉱排鉱部に到達すると、順次生産した焼結鉱を排鉱するドワイト・ロイド式の焼結機の操業状態を監視する焼結プロセス操業監視装置であって、
前記パレット上に形成されている焼結層内の測定温度データであって、高さ方向及びパレット幅方向のそれぞれに間隔を有して前記パレット上に配置された複数の温度測定部が前記焼結機の原料給鉱部から焼結鉱排鉱部まで移動する間に当該複数の温度測定部によって測定された時系列の測定温度データを収集する測定温度データ収集手段と、
前記測定温度データ収集手段により収集された測定温度データを用いて、前記焼結層を焼結機の機長方向に、原料給鉱部から焼結鉱排鉱部まで区間を初期原料帯、水分凝縮帯、乾燥帯、赤熱帯、及び焼結完了帯の５つのゾーンに分けるゾーン解析手段と、
前記ゾーン毎に個別に設定される関数であって、前記焼結層内の機長方向の温度の変化を定量的に表現する関数のパラメータとして、前記複数の温度測定部が配置された高さ位置のそれぞれについて、個別に与えられるパラメータを前記複数の温度測定部によって測定された測定温度データに基づいてそれぞれ同定する測定温度関数のパラメータ同定手段と、
前記測定温度関数のパラメータ同定手段により同定されたパラメータが与えられた前記複数の関数上の点であって、前記ゾーン内での位置が相互に対応する点に基づく外挿又は内挿を行って、前記焼結層内の温度測定部が配置された高さ位置以外の指定された高さ位置の温度データである導出対象温度データを導出する導出対象温度データ導出手段と、
前記焼結層内の機長方向の温度の変化を定量的に表現する関数のパラメータを、前記導出対象温度データに基づいて同定する導出対象温度関数のパラメータ同定手段と、
前記焼結層内の温度の分布の情報を表示装置に表示する表示手段と、を有し、
前記測定温度関数のパラメータ同定手段及び前記導出対象温度関数のパラメータ同定手段で同定するパラメータは、
前記初期原料帯では同定するパラメータがなく、
前記水分凝縮帯では空気流入温度、水分気化速度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、
前記乾燥帯では空気流入温度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、
前記赤熱帯では空気流入温度、燃料燃焼速度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、
前記焼結完了帯では空気流入温度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、
前記同定されたパラメータが与えられた導出対象温度関数によって、前記焼結層内の高さ方向における温度分布が得られるようにしたことを特徴とする焼結プロセス操業監視装置。
【請求項２】
前記パレット幅方向で異なる位置の前記測定温度データ又は前記導出対象温度データであって、ゾーン内で相互に対応する位置の前記測定温度データ又は前記導出対象温度データに基づく外挿又は内挿を行って、前記焼結層内のパレット幅方向における温度の変化を表す関数を導出するパレット幅方向関数導出手段を有し、
前記焼結層内のパレット幅方向における温度の変化を表す関数によって、前記焼結層内のパレット幅方向における温度分布が得られるようにしたことを特徴とする請求項１に記載の焼結プロセス操業監視装置。
【請求項３】
前記指定された高さ位置は、前記焼結層の床部に略対応する位置を含み、
前記５つのゾーンの少なくとも１つのゾーンには、赤熱帯が含まれており、
前記導出対象温度関数のパラメータ同定手段により同定されたパラメータが与えられた関数は、前記赤熱帯における前記焼結層の床部における機長方向の温度の変化を定量的に表現する関数を含み、
前記パレット幅方向関数導出手段は、前記パレット幅方向で異なる位置の前記導出対象温度データであって、前記赤熱帯の終点に対応する位置の前記導出対象温度データに基づく外挿又は内挿を行って、ＢＴＰ（Burn Through Point）のパレット幅方向の分布を導出することを特徴とする請求項２に記載の焼結プロセス操業監視装置。
【請求項４】
前記焼結層を構成する焼結鉱粒子が球形粒子であるとして取り扱い、
前記関数は、球形粒子の熱収支を表す時間微分方程式であることを特徴とする請求項１〜３の何れか１項に記載の焼結プロセス操業監視装置。
【請求項５】
前記表示手段は、焼結層内の温度が同値である等値線を用いて、前記焼結層内の温度の分布の情報を表示することを特徴とする請求項１〜４の何れか１項に記載の焼結プロセス操業監視装置。
【請求項６】
無端帯状になるように連結されて周回する複数のパレットのうち、原料給鉱部に位置するパレットに順次焼結原料を給鉱し、当該パレットが移動する過程で当該焼結原料の焼結を行って焼結鉱を生産し、当該パレットが焼結鉱排鉱部に到達すると、順次生産した焼結鉱を排鉱するドワイト・ロイド式の焼結機の操業状態を監視する焼結プロセス操業監視方法であって、
前記パレット上に形成されている焼結層内の測定温度データであって、高さ方向及びパレット幅方向のそれぞれに間隔を有して前記パレット上に配置された複数の温度測定部が前記焼結機の原料給鉱部から焼結鉱排鉱部まで移動する間に当該複数の温度測定部によって測定された時系列の測定温度データを収集する測定温度データ収集工程と、
前記測定温度データ収集工程により収集された測定温度データを用いて、前記焼結層を焼結機の機長方向に、原料給鉱部から焼結鉱排鉱部まで区間を初期原料帯、水分凝縮帯、乾燥帯、赤熱帯、及び焼結完了帯の５つのゾーンに分けるゾーン解析工程と、
前記ゾーン毎に個別に設定される関数であって、前記焼結層内の機長方向の温度の変化を定量的に表現する関数のパラメータとして、前記複数の温度測定部が配置された高さ位置のそれぞれについて、個別に与えられるパラメータを前記複数の温度測定部によって測定された測定温度データに基づいてそれぞれ同定する測定温度関数のパラメータ同定工程と、
前記測定温度関数のパラメータ同定工程により同定されたパラメータが与えられた前記複数の関数上の点であって、前記ゾーン内での位置が相互に対応する点に基づく外挿又は内挿を行って、前記焼結層内の温度測定部が配置された高さ位置以外の指定された高さ位置の温度データである導出対象温度データを導出する導出対象温度データ導出工程と、
前記焼結層内の機長方向の温度の変化を定量的に表現する関数のパラメータを、前記導出対象温度データに基づいて同定する導出対象温度関数のパラメータ同定工程と、
前記焼結層内の温度の分布の情報を表示装置に表示する表示工程と、を有し、
前記測定温度関数のパラメータ同定工程及び前記導出対象温度関数のパラメータ同定工程で同定するパラメータは、
前記初期原料帯では同定するパラメータがなく、
前記水分凝縮帯では空気流入温度、水分気化速度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、
前記乾燥帯では空気流入温度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、
前記赤熱帯では空気流入温度、燃料燃焼速度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、
前記焼結完了帯では空気流入温度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、
前記同定されたパラメータが与えられた導出対象温度関数によって、前記焼結層内の高さ方向における温度分布が得られるようにしたことを特徴とする焼結プロセス操業監視方法。
【請求項７】
無端帯状になるように連結されて周回する複数のパレットのうち、原料給鉱部に位置するパレットに順次焼結原料を給鉱し、当該パレットが移動する過程で当該焼結原料の焼結を行って焼結鉱を生産し、当該パレットが焼結鉱排鉱部に到達すると、順次生産した焼結鉱を排鉱するドワイト・ロイド式の焼結機の操業状態を監視するためのプログラムであって、
前記パレット上に形成されている焼結層内の測定温度データであって、高さ方向及びパレット幅方向のそれぞれに間隔を有して前記パレット上に配置された複数の温度測定部が前記焼結機の原料給鉱部から焼結鉱排鉱部まで移動する間に当該複数の温度測定部によって測定された時系列の測定温度データを収集する測定温度データ収集処理と、
前記測定温度データ収集処理により収集された測定温度データを用いて、前記焼結層を焼結機の機長方向に、原料給鉱部から焼結鉱排鉱部まで区間を初期原料帯、水分凝縮帯、乾燥帯、赤熱帯、及び焼結完了帯の５つのゾーンに分けるゾーン解析処理と、
前記ゾーン毎に個別に設定される関数であって、前記焼結層内の機長方向の温度の変化を定量的に表現する関数のパラメータとして、前記複数の温度測定部が配置された高さ位置のそれぞれについて、個別に与えられるパラメータを前記複数の温度測定部によって測定された測定温度データに基づいてそれぞれ同定する測定温度関数のパラメータ同定処理と、
前記測定温度関数のパラメータ同定処理により同定されたパラメータが与えられた前記複数の関数上の点であって、前記ゾーン内での位置が相互に対応する点に基づく外挿又は内挿を行って、前記焼結層内の温度測定部が配置された高さ位置以外の指定された高さ位置の温度データである導出対象温度データを導出する導出対象温度データ導出処理と、
前記焼結層内の機長方向の温度の変化を定量的に表現する関数のパラメータを、前記導出対象温度データに基づいて同定する導出対象温度関数のパラメータ同定処理と、
前記焼結層内の温度の分布の情報を表示装置に表示する表示処理と、をコンピュータに実行させ、
前記測定温度関数のパラメータ同定処理及び前記導出対象温度関数のパラメータ同定処理で同定するパラメータは、
前記初期原料帯では同定するパラメータがなく、
前記水分凝縮帯では空気流入温度、水分気化速度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、
前記乾燥帯では空気流入温度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、
前記赤熱帯では空気流入温度、燃料燃焼速度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、
前記焼結完了帯では空気流入温度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、
前記同定されたパラメータが与えられた導出対象温度関数によって、前記焼結層内の高さ方向における温度分布が得られるようにしたことを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、焼結プロセス操業監視装置、焼結プロセス操業監視方法、及びプログラムに関し、特に、高炉に供給する焼結鉱を生成する焼結機における焼結層の状態を解析するために用いて好適なものである。
【背景技術】
【０００２】
高炉に供給する焼結鉱を生産する焼結機では、銑鉄の主原料である微粉鉄鉱石を燃料となるコークス等とともに焼き固め、その粒度・強度・化学成分を調整して焼結鉱を製造する。図１に示したようなドワイト・ロイド式の焼結機による焼結鉱の生産方法の概要を説明すると、まず、微粉鉄鉱石を主原料とする焼結原料を、原料供給ホッパから、無端帯状に連結され移動する複数のパレット上へ移動方向上流側で給鉱して焼結原料層を形成し、当該焼結原料層の上部を、点火炉によって着火する。そして、パレットが移動しながら複数のパレットの下方に配置された複数のウインドボックス（風箱）を介して、焼結原料が積載されたパレットを下方から吸気する。これにより、複数のパレットに積載された焼結原料は、焼結機の機長方向に搬送されながら、焼結原料層の上方から下方に向けて焼結反応が進行する。そして、焼結原料層の上部から下部までの焼成が完了すると、得られた焼結鉱が各パレットから順次排鉱される。本願では、焼結反応の進行に際して焼結原料層から焼結鉱層へと変化する材料層を、ともに焼結層と記す。
【０００３】
このような焼結機による操業には、以下の（１）〜（３）の事項が求められる。
（１）微粉鉄鉱石の配合割合等の原料条件や操業条件の変動を極力少なくする。
（２）原料条件や操業条件の変動に対し、系全体をうまく制御することで、高炉で還元・溶融するための一定の粒度と強度をもつ焼結鉱を連続的に安定して生産する。
（３）鉄鋼需給の変動や製鉄所の方針等、種々の制約条件のもとで、コストの低減・歩留まりの向上等の操業目標値を最大にするように操業する。
【０００４】
また、生産物としての焼結鉱の品質指標には、冷間強度指数（ＳＩ：Shutter Index）、還元粉化性指数（ＲＤＩ：Reduction Disintegration Index）、被還元性指数（ＲＩ：Reducibility Index）があり、焼結鉱には、以下の（１）〜（３）を満足するような品質が求められる。
（１）製造後の搬送時や高炉装入時に粉化しないよう十分に大きい冷間強度を焼結鉱が有すること。
（２）高炉内での通気性を確保するため、著しい還元粉化を起こさないような十分に小さい還元粉化性指数を焼結鉱が有すること。
（３）高炉装入物として良好な被還元性を保証できるような十分に大きい被還元性指数を焼結鉱が有すること。
【０００５】
微粉鉄鉱石がコークスの燃焼によって適切に溶融・塊成化しているか否かの監視指標として、（焼結層底面における）赤熱帯終点位置を焼結完了点（ＢＴＰ：Burn Through Point）と定義している。一般的に、このＢＴＰが焼結機のどこに位置されているかを操業管理する。前述したような焼結鉱の品質を満足し、且つ高生産性（高焼結鉱生産量［ｔｏｎ／ｈ］）を実現することを操業目標として、基本的にはＢＴＰがなるべく焼結機の排鉱部側に安定的に存在する範囲でパレット速度を最大化するような操業を行う。
そのために、従来は、ウインドボックスにおける排ガス温度が最高となる点がＢＴＰと対応するという仮定のもと、この排ガス温度から、ＢＴＰを予測することが行われている。
【０００６】
特許文献１では、排ガス温度計で測定された排ガス温度に基づいて、機長方向における温度が最高となる位置を焼結完了点ＢＴＰとして求める。また、排ガス温度計で測定された排ガス温度に基づいて、温度上昇開始点から焼結終了点までの間において所定温度になる位置を所定温度位置ＢＲＰとして求める。そして、焼結完了点ＢＴＰと、所定温度位置ＢＲＰとの相関を算出しておき、この相関を利用して、現時点で求められた所定温度位置ＢＲＰから焼結完了点ＢＴＰを求める。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００６−３０７２５９号公報
【特許文献２】特開２０１０−２４３４４３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
しかしながら、前述したように、従来の技術では、ウインドボックスにおける排ガス温度が最高となる点がＢＴＰと対応するという仮定のもと、この排ガス温度から、ＢＴＰを予測しているに過ぎない。ウインドボックスの寸法は、機長方向で約５［ｍ］、幅方向で約５［ｍ］あり、このような大きな空間を通過してくる排ガス温度の測定の結果からでは、焼結層における焼成の進行のばらつきを正確に把握することができない。このように従来の技術では、ＢＴＰの管理精度に限界がある。
【０００９】
そこで、本出願人らは、特許文献２において、焼結機の焼結層の温度を直接的に測定することができる温度測定装置（以下、この温度測定装置を必要に応じて「パレット温度計」と称する）を提案した。さらに、本出願人らは、特願２０１１−１８４０１７号において、パレット温度計により測定された温度データから、焼結層内の物理現象を定量的に表現するパラメータである特徴量を、焼結層のゾーン毎に同定する技術を提案した。
【００１０】
この特願２０１１−１８４０１７号に記載の技術によれば、焼結層の底部に相当する位置にパレット温度計の測温体を配置することによって、ＢＴＰを直接的に求めることができる。
しかしながら、特願２０１１−１８４０１７号に記載の技術では、パレット温度計の測温体が配置されている箇所の温度データしか得られないため、焼結層内の所望の位置での温度分布を求めることができない。したがって、例えば、焼結層の底部に相当する位置にパレット温度計の測温体を配置しなければ、ＢＴＰを求めることが困難である。
【００１１】
本発明は、このような問題点に鑑みてなされたものであり、焼結機における焼結層内の所望の領域での温度分布をオンラインで正確に同定して、焼結機の操業を監視できるようにすることを目的とする。
【課題を解決するための手段】
【００１２】
本発明の焼結プロセス操業監視装置は、無端帯状になるように連結されて周回する複数のパレットのうち、原料給鉱部に位置するパレットに順次焼結原料を給鉱し、当該パレットが移動する過程で当該焼結原料の焼結を行って焼結鉱を生産し、当該パレットが焼結鉱排鉱部に到達すると、順次生産した焼結鉱を排鉱するドワイト・ロイド式の焼結機の操業状態を監視する焼結プロセス操業監視装置であって、前記パレット上に形成されている焼結層内の測定温度データであって、高さ方向及びパレット幅方向のそれぞれに間隔を有して前記パレット上に配置された複数の温度測定部が前記焼結機の原料給鉱部から焼結鉱排鉱部まで移動する間に当該複数の温度測定部によって測定された時系列の測定温度データを収集する測定温度データ収集手段と、前記測定温度データ収集手段により収集された測定温度データを用いて、前記焼結層を焼結機の機長方向に、原料給鉱部から焼結鉱排鉱部まで区間を初期原料帯、水分凝縮帯、乾燥帯、赤熱帯、及び焼結完了帯の５つのゾーンに分けるゾーン解析手段と、前記ゾーン毎に個別に設定される関数であって、前記焼結層内の機長方向の温度の変化を定量的に表現する関数のパラメータとして、前記複数の温度測定部が配置された高さ位置のそれぞれについて、個別に与えられるパラメータを前記複数の温度測定部によって測定された測定温度データに基づいてそれぞれ同定する測定温度関数のパラメータ同定手段と、前記測定温度関数のパラメータ同定手段により同定されたパラメータが与えられた前記複数の関数上の点であって、前記ゾーン内での位置が相互に対応する点に基づく外挿又は内挿を行って、前記焼結層内の温度測定部が配置された高さ位置以外の指定された高さ位置の温度データである導出対象温度データを導出する導出対象温度データ導出手段と、前記焼結層内の機長方向の温度の変化を定量的に表現する関数のパラメータを、前記導出対象温度データに基づいて同定する導出対象温度関数のパラメータ同定手段と、前記焼結層内の温度の分布の情報を表示装置に表示する表示手段と、を有し、前記測定温度関数のパラメータ同定手段及び前記導出対象温度関数のパラメータ同定手段で同定するパラメータは、前記初期原料帯では同定するパラメータがなく、前記水分凝縮帯では空気流入温度、水分気化速度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、前記乾燥帯では空気流入温度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、前記赤熱帯では空気流入温度、燃料燃焼速度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、前記焼結完了帯では空気流入温度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、
前記同定されたパラメータが与えられた導出対象温度関数によって、前記焼結層内の高さ方向における温度分布が得られるようにしたことを特徴とする。
また、本発明の焼結プロセス操業監視装置の他の特徴とするところは、前記パレット幅方向で異なる位置の前記測定温度データ又は前記導出対象温度データであって、ゾーン内で相互に対応する位置の前記測定温度データ又は前記導出対象温度データに基づく外挿又は内挿を行って、前記焼結層内のパレット幅方向における温度の変化を表す関数を導出するパレット幅方向関数導出手段を有し、前記焼結層内のパレット幅方向における温度の変化を表す関数によって、前記焼結層内のパレット幅方向における温度分布が得られるようにした点にある。さらに、前記指定された高さ位置は、前記焼結層の床部に略対応する位置を含み、前記５つのゾーンの少なくとも１つのゾーンには、赤熱帯が含まれており、前記導出対象温度関数のパラメータ同定手段により同定されたパラメータが与えられた関数は、前記赤熱帯における前記焼結層の床部における機長方向の温度の変化を定量的に表現する関数を含み、前記パレット幅方向関数導出手段は、前記パレット幅方向で異なる位置の前記導出対象温度データであって、前記赤熱帯の終点に対応する位置の前記導出対象温度データに基づく外挿又は内挿を行って、ＢＴＰ（Burn Through Point）のパレット幅方向の分布を導出する点にある。
また、本発明の焼結プロセス操業監視装置の他の特徴とするところは、前記焼結層を構成する焼結鉱粒子が球形粒子であるとして取り扱い、前記関数は、球形粒子の熱収支を表す時間微分方程式である点にある。
また、本発明の焼結プロセス操業監視装置の他の特徴とするところは、前記表示手段は、焼結層内の温度が同値である等値線を用いて、前記焼結層内の温度の分布の情報を表示する点にある。
本発明の焼結プロセス操業監視方法は、無端帯状になるように連結されて周回する複数のパレットのうち、原料給鉱部に位置するパレットに順次焼結原料を給鉱し、当該パレットが移動する過程で当該焼結原料の焼結を行って焼結鉱を生産し、当該パレットが焼結鉱排鉱部に到達すると、順次生産した焼結鉱を排鉱するドワイト・ロイド式の焼結機の操業状態を監視する焼結プロセス操業監視方法であって、前記パレット上に形成されている焼結層内の測定温度データであって、高さ方向及びパレット幅方向のそれぞれに間隔を有して前記パレット上に配置された複数の温度測定部が前記焼結機の原料給鉱部から焼結鉱排鉱部まで移動する間に当該複数の温度測定部によって測定された時系列の測定温度データを収集する測定温度データ収集工程と、前記測定温度データ収集工程により収集された測定温度データを用いて、前記焼結層を焼結機の機長方向に、原料給鉱部から焼結鉱排鉱部まで区間を初期原料帯、水分凝縮帯、乾燥帯、赤熱帯、及び焼結完了帯の５つのゾーンに分けるゾーン解析工程と、前記ゾーン毎に個別に設定される関数であって、前記焼結層内の機長方向の温度の変化を定量的に表現する関数のパラメータとして、前記複数の温度測定部が配置された高さ位置のそれぞれについて、個別に与えられるパラメータを前記複数の温度測定部によって測定された測定温度データに基づいてそれぞれ同定する測定温度関数のパラメータ同定工程と、前記測定温度関数のパラメータ同定工程により同定されたパラメータが与えられた前記複数の関数上の点であって、前記ゾーン内での位置が相互に対応する点に基づく外挿又は内挿を行って、前記焼結層内の温度測定部が配置された高さ位置以外の指定された高さ位置の温度データである導出対象温度データを導出する導出対象温度データ導出工程と、前記焼結層内の機長方向の温度の変化を定量的に表現する関数のパラメータを、前記導出対象温度データに基づいて同定する導出対象温度関数のパラメータ同定工程と、前記焼結層内の温度の分布の情報を表示装置に表示する表示工程と、を有し、前記測定温度関数のパラメータ同定工程及び前記導出対象温度関数のパラメータ同定工程で同定するパラメータは、前記初期原料帯では同定するパラメータがなく、前記水分凝縮帯では空気流入温度、水分気化速度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、前記乾燥帯では空気流入温度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、前記赤熱帯では空気流入温度、燃料燃焼速度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、前記焼結完了帯では空気流入温度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、前記同定されたパラメータが与えられた導出対象温度関数によって、前記焼結層内の高さ方向における温度分布が得られるようにしたことを特徴とする。
本発明のプログラムは、無端帯状になるように連結されて周回する複数のパレットのうち、原料給鉱部に位置するパレットに順次焼結原料を給鉱し、当該パレットが移動する過程で当該焼結原料の焼結を行って焼結鉱を生産し、当該パレットが焼結鉱排鉱部に到達すると、順次生産した焼結鉱を排鉱するドワイト・ロイド式の焼結機の操業状態を監視するためのプログラムであって、前記パレット上に形成されている焼結層内の測定温度データであって、高さ方向及びパレット幅方向のそれぞれに間隔を有して前記パレット上に配置された複数の温度測定部が前記焼結機の原料給鉱部から焼結鉱排鉱部まで移動する間に当該複数の温度測定部によって測定された時系列の測定温度データを収集する測定温度データ収集処理と、前記測定温度データ収集処理により収集された測定温度データを用いて、前記焼結層を焼結機の機長方向に、原料給鉱部から焼結鉱排鉱部まで区間を初期原料帯、水分凝縮帯、乾燥帯、赤熱帯、及び焼結完了帯の５つのゾーンに分けるゾーン解析処理と、前記ゾーン毎に個別に設定される関数であって、前記焼結層内の機長方向の温度の変化を定量的に表現する関数のパラメータとして、前記複数の温度測定部が配置された高さ位置のそれぞれについて、個別に与えられるパラメータを前記複数の温度測定部によって測定された測定温度データに基づいてそれぞれ同定する測定温度関数のパラメータ同定処理と、前記測定温度関数のパラメータ同定処理により同定されたパラメータが与えられた前記複数の関数上の点であって、前記ゾーン内での位置が相互に対応する点に基づく外挿又は内挿を行って、前記焼結層内の温度測定部が配置された高さ位置以外の指定された高さ位置の温度データである導出対象温度データを導出する導出対象温度データ導出処理と、前記焼結層内の機長方向の温度の変化を定量的に表現する関数のパラメータを、前記導出対象温度データに基づいて同定する導出対象温度関数のパラメータ同定処理と、前記焼結層内の温度の分布の情報を表示装置に表示する表示処理と、をコンピュータに実行させ、前記測定温度関数のパラメータ同定処理及び前記導出対象温度関数のパラメータ同定処理で同定するパラメータは、前記初期原料帯では同定するパラメータがなく、前記水分凝縮帯では空気流入温度、水分気化速度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、前記乾燥帯では空気流入温度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、前記赤熱帯では空気流入温度、燃料燃焼速度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、前記焼結完了帯では空気流入温度、対流熱伝達率、時定数、焼結鉱粒子の形状係数、及び空気流速であり、前記同定されたパラメータが与えられた導出対象温度関数によって、前記焼結層内の高さ方向における温度分布が得られるようにしたことを特徴とする。
【発明の効果】
【００１３】
本発明によれば、ゾーン毎に個別に設定される、焼結層内の機長方向の温度の変化を定量的に表現する関数のパラメータを、異なる高さ位置で測定された測定温度データに基づいて、当該高さ位置のそれぞれで同定する。そして、同定した複数の関数上の、ゾーン内の位置が相互に対応する点に基づく外挿又は内挿を行って、焼結層内の指定された高さ位置の温度データである導出対象温度データを導出する。そして、導出した導出対象温度データに基づいて、焼結層内の機長方向の温度の変化を定量的に表現する関数のパラメータであって、当該指定された高さ位置におけるパラメータを同定する。したがって、焼結機における焼結層内の高さ方向の温度分布をオンラインで正確に同定して、焼結機の操業を監視することができる。
また、本発明によれば、ゾーン内で相互に対応する位置の測定温度データ又は導出対象温度データに基づく外挿又は内挿を行って、焼結層内のパレット幅方向の温度の変化を表す関数を導出するようにした。したがって、焼結機における焼結層内のパレット幅方向の温度分布をオンラインで正確に同定して、焼結機の操業を監視することができる。
以上により、焼結機における焼結層内の機長方向、高さ方向、及びパレット幅方向の３次元の温度分布をオンラインで正確に同定して、例えば焼結機の操業指標として重要なＢＴＰ（Burn Through Point）のパレット幅方向の分布を導出して、焼結機の操業を監視することができる。
【図面の簡単な説明】
【００１４】
【図１】本発明の実施形態を示し、焼結機の概略構成の一例を示す図である。
【図２】本発明の実施形態を示し、パレット温度計が設置されたパレットの概略構成の一例を示す図である。
【図３】本発明の実施形態を示し、パレット温度計の概略構成の一例を示す図である。
【図４】本発明の実施形態を示し、焼結プロセス操業監視装置の機能的な構成の一例を示すブロック図である。
【図５】本発明の実施形態を示し、焼結層を区画するゾーンの一例を示す図である。
【図６】本発明の実施形態を示し、焼結層のある高さ位置における測定温度データの一例を概念的に示す図である。
【図７】本発明の実施形態を示し、焼結層内の焼結鉱粒子充填層を概念的に示す図である。
【図８】本発明の実施形態を示し、焼結層のゾーンIV（赤熱帯）の領域と、測定温度データと、導出対象温度データの一例を示す図である。
【図９】本発明の実施形態を示し、測定温度データに基づく関数から導出対象温度データに基づく関数を求める様子の一例を概念的に示す図である。
【図１０】本発明の実施形態を示し、あるゾーンにおける相互に対応する位置でのパレット幅方向の温度の変化を表す関数の一例を示す図である。
【図１１】本発明の実施形態を示し、焼結層のゾーンIV（赤熱帯）の概略の一例を３次元的に表現して示した図である。
【図１２】本発明の実施形態を示し、焼結プロセス操業監視装置の処理動作の一例を説明するフローチャートである。
【発明を実施するための形態】
【００１５】
本出願人は、特許文献２において、焼結層内の温度を、連続して安定的に直接的に測定することができるパレット温度計を提案した。このパレット温度計を用いることにより、焼結層内の温度分布を初めて実測することができた。本発明者らは、このように従来では得られなかった「焼結層内の温度分布の実測値」を鋭意検討することにより、後述するようにして、焼結層内の物理現象を定量的に表現するパラメータを焼結状態の特徴量としてオンラインで正確に同定することができることを見出した。すなわち、従来では、ウインドボックスにおける排ガス温度から、焼結層内の温度を間接的に推測するに留まっているに過ぎず（焼結層内の温度を経験に基づいて求めているに過ぎず）、実プロセスの焼結層内の温度データ実測値に基づく焼結層の状態の理論的な解明は十分に行われていなかった。これに対し、本発明者らは、焼結層内の温度分布の実測値から、焼結層の状態の理論的な解明を行い、その結果に基づいて、焼結層内の物理現象を定量的に表現するパラメータを表現するモデルを構築し、構築したモデルを用いて、焼結層内の物理現象を定量的に表現するパラメータを、焼結状態の特徴量として同定し、これを用いて、焼結機の操業をオンラインで監視できることを見出した。
【００１６】
更に、本発明者らは、このようにして焼結機による焼結プロセスの操業状態の特徴量を同定した結果から、焼結層内の所望の領域での２次元又は３次元の温度分布を導出することができることを見出した。これにより、例えば、焼結層全体の温度分布を等値線等で表現したり、焼結層の底部にパレット温度計の測温体を配置しなくてもＢＴＰを求めたりすることができる。そして、これらの情報を表示し、さらにこれらの表示を逐次更新することにより、焼結機の操業をオンラインで監視することができる。
【００１７】
以下に、図面を参照しながら、本発明の一実施形態について説明する。
＜焼結機＞
まず、焼結機の概略について説明する。
図１は、焼結機の概略構成の一例を示す図である。図１は、焼結機１００をその側方から見た図である。また、図１では、説明の都合上、必要な部分のみを簡略化して示している。
図１において、焼結機１００は、複数のパレット１１０からなるパレット群と、原料供給ホッパ１２０と、点火炉１３０と、ウインドボックス１４０とを備えて構成される。
【００１８】
パレット群は、焼結原料を積載して搬送する複数のパレット１１０を焼結機１００の長手方向に移動可能に連結して構成される。図１に示すように、各パレット１１０は、先端と後端とが連接されて全体として無端帯状になっており、駆動ローラ１０２ａ、１０２ｂの回転にしたがってエンドレスに焼結機１００内を周回する。図１では、駆動ローラ１０２ａ、１０２ｂが時計回りに回転することによって、パレット群は時計回りに周回している。パレット群を構成するパレット１１０のうち少なくとも１つのパレット１１０に、パレット１１０上に積載された焼結原料の温度を測定する複数のパレット温度計１５０ａ〜１５０ｆ（図２を参照）が設けられる。
【００１９】
原料供給ホッパ１２０は、焼結鉱の原料となる焼結原料を供給する。原料供給ホッパ１２０は、点火炉１３０の位置を基準としてパレット群の回転方向上流側に設けられる。原料供給ホッパ１２０から供給された焼結原料（燃料を含む）は、原料供給ホッパ１２０の下方を通過するパレット１１０の収容空間に収容され、パレット群の進行方向へ搬送される。
点火炉１３０は、パレット１１０上に積載された焼結原料に着火する点火部である。点火炉１３０と対向するパレット１１０上に積載された焼結原料に着火して、焼結原料の焼結過程を開始させる。
【００２０】
ウインドボックス１４０は、焼結原料が積載されたパレット１１０を下方から吸気するための減圧空間を構成する。ウインドボックス１４０は、パレット１１０の下方に設けられた空間であり、ウインドボックス１４０の上部は、パレット１１０の底面下部と連通し、ウインドボックス１４０の下部は、吸気管１４５に連通している。ウインドボックス１４０は、吸気管１４５を介してブロア（図示せず）と接続されており、ブロアによって吸気／排気されることでウインドボックス１４０の内部空間が減圧される。
【００２１】
パレット群の周回中、まず、各パレット１１０に原料供給ホッパ１２０から焼結原料が供給され、焼結原料の搬送が開始される。次いで、パレット１１０の積載する焼結原料の上部が点火炉１３０により着火されると共に、ウインドボックス１４０及び吸気管１４５を介してブロアにより吸気される。これにより、パレット群が移動する間に焼結層の上側表面から下方に燃焼帯を進行させ、焼結鉱を連続的に生産する。生産された焼結鉱は、排鉱部にパレット１１０が到達すると、パレット１１０上から排出される。そして、焼結鉱が排出されたパレット１１０は、原料供給ホッパ１２０まで周回して搬送される。
【００２２】
＜パレット温度計＞
次に、パレット温度計の概略について説明する。尚、パレット温度計の詳細については、特許文献２に記載されているので、ここでは、その概略のみを説明する。
図２は、パレット温度計が設置されたパレットの概略構成の一例を示す図である。
図２において、パレット１１０は、複数のグレートバー１１２と、サイドウォール１１４と、車輪１１６と、６個のパレット温度計１５０ａ〜１５０ｆとを備えて構成される。
【００２３】
グレートバー１１２は、パレット１１０の底面を構成する部材である。複数のグレートバー１１２のうち、パレット幅方向（ｘ軸方向）について所定の間隔で６個のグレートバー１１２の代わりにパレット温度計１５０ａ〜１５０ｆが配置される。尚、パレット温度計１５０ａ〜１５０ｆはいずれも同じ構成を有し、以下では「パレット温度計１５０」とも記す。
サイドウォール１１４は、複数のグレートバー１１２によって構成された底面の、パレット幅方向（ｘ軸方向）における両側に設けられた壁板である。サイドウォール１１４は、パレット１１０上に積載された焼結原料や生成された焼結鉱がパレット１１０の進行方向に対して側方向（ｘ軸方向）に落下するのを防止する。
【００２４】
車輪１１６は、パレット１１０を焼結機１００の機長方向（ｙ軸方向）に移動させるための部材である。図２に示す例では、車輪１１６は、１つのパレット１１０に対して各サイドウォール１１４の外側面にそれぞれ２つずつ設けられる。
パレット温度計１５０は、パレット１１０に積載された焼結層内の温度を測定する装置である。ここで、焼結層とは、パレット１１０上に積載された焼結原料の層と、生成された焼結鉱の層との何れの層も含むものとする。パレット温度計１５０は、前述したように、パレット１１０の底面を構成するために配列されたグレートバー１１２と共に配列される。本実施形態では、ｙ軸正方向側のグレートバー列において、パレット幅方向（ｘ軸方向）に６個のパレット温度計１５０ａ〜１５０ｆを設けるようにしている。具体的に、本実施形態では、ｙ軸正方向側のグレートバー列における、サイドウォール１１４に直近の位置に、パレット温度計１５０ａ、１５０ｆが配置されるようにしている。また、相互に隣り合うパレット温度計１５０の間隔は略同じになるように、各パレット温度計１５０ａ〜１５０ｆが配置されるようにしている。
【００２５】
尚、図２では、パレット温度計１５０を６個設けた場合を例に挙げて示しているが、１つのパレット１１０に設けるパレット温度計１５０の数は、６個に限定されない。パレット幅方向（ｘ軸方向）の所望の位置における高さ方向（ｚ軸方向）の温度分布を求める際には、１つのパレット１１０に設けるパレット温度計１５０の数は１個でもよい。ただし、３次元の温度分布を求める際には、パレット幅方向（ｘ軸方向）の複数の位置での温度のデータが必要になるので、１つのパレット１１０においてパレット幅方向（ｘ軸方向）に設けるパレット温度計１５０の数は２個以上必要になる。１つのパレット１１０に複数のパレット温度計１５０を設ける場合、それら複数のパレット温度計１５０の配置は、焼結層内のパレット幅方向（ｘ軸方向）の複数の箇所での温度が測定可能な配置であれば、どのような配置であってもよい。そして、１つのパレット１１０内において複数のパレット温度計１５０を配置するとき、複数のパレット温度計１５０を図２で例示するにようにパレット幅方向（ｘ軸方向）に一列に配置してもよいし、機長方向（ｙ軸方向）に対して複数の列を構成した配置であってもよい。また、本実施形態では、１つのパレット１１０にパレット温度計１５０を設ける場合を例に挙げて説明するが、パレット温度計１５０は、全てのパレット１１０に設けても、一部（複数）のパレット１１０に設けてもよい。
【００２６】
＜パレット温度計の構成＞
図３は、パレット温度計１５０の概略構成の一例を示す図である。具体的に図３は、パレット温度計１５０をｙ軸方向及びｚ軸方向の双方に沿うようにして切ったときの断面図である。
パレット温度計１５０は、高さ方向（ｚ軸方向）で間隔を隔てて配置された２つの温度測定部１７０ａ、１７０ｂを保持する筐体（スタンド）１６０と、パレット１１０上に積載された焼結層内の温度を測定する温度測定部１７０とを備えて構成される。
【００２７】
筐体１６０は、グレートバー１１２の基台部と共にパレット１１０の底面を構成する基台部１６１と、基台部１６１からパレット１１０の収容空間に向かってｚ軸正方向に突出する突出部１６２と、グレートバー１１２と共に筐体１６０をパレット１１０の底面に配列するための脚部１６３ａ、１６３ｂとを備えて構成される。筐体１６０は、例えば耐熱鋳物であって、各部が一体に形成される。
突出部１６２の内部には、温度測定部１７０ａ、１７０ｂが挿通される空間が形成される。すなわち、突出部１６２の内部には、温度測定部１７０ａが挿入される水平孔１６４ａと、温度測定部１７０ｂが挿入される水平孔１６４ｂと、水平孔１６４ａ、１６４ｂと突出部１６２の底面とを連結する垂直孔１６５とが形成される。
【００２８】
水平孔１６４ａ、１６４ｂに挿入された温度測定部１７０ａ、１７０ｂの先端部は、パレット１１０の進行方向に向かって筐体１６０から突出している。
また、水平孔１６４ｂ（の中心）は、基台部１６１の上面からｚ軸正方向に向かって３０［ｍｍ］の位置に形成され、水平孔１６４ａ（の中心）は、基台部１６１の上面からｚ軸正方向に向かって１８０［ｍｍ］の位置（水平孔１６４ｂ（の中心）からｚ軸正方向に向かって１５０［ｍｍ］の位置）に形成されている。水平孔１６４ａ、水平孔１６４ｂのｚ軸方向における位置は、測温したい焼結層の位置を規定する。すなわち、これらの水平孔１６４ａ、１６４ｂに設けられる温度測定部１７０ａ、１７０ｂによって、パレット１１０に積載された焼結層内の温度を高さ方向（ｚ軸方向）に複数箇所で測定することができる。
【００２９】
尚、ここでは、高さ方向（ｚ軸方向）に間隔を有して２つの水平孔１６４ａ、１６４ｂ（温度測定部１７０ａ、１７０ｂ）を、１つのパレット温度計１５０に設けるようにした場合を例に挙げて示した。しかしながら、１つのパレット温度計１５０に設ける温度測定部１７０の数は２つに限定されない。ある高さにおける、パレット幅方向（ｘ軸方向）及び機長方向（ｙ軸方向）の温度分布を求める場合には、１つのパレット温度計１５０に設ける温度測定部１７０の数は１つであってもよい。ただし、３次元の温度分布を求める際には、高さ方向（ｚ軸方向）の複数の位置での温度のデータが必要になるので、１つのパレット温度計１５０に設ける温度測定部１７０の数は高さ方向（ｚ軸方向）に２つ以上必要になる。１つのパレット温度計１５０に複数の温度測定部１７０を設ける場合、それら複数の温度測定部１７０の配置は、焼結層内の高さ方向（ｚ軸方向）の複数の箇所での温度が測定可能な配置であれば、どのような配置であってもよい。
【００３０】
温度測定部１７０ａ、１７０ｂは、それぞれ同じ構成を有し、焼結層内の温度を測定する測温体と、測温体を包囲する保護管とを備える。測温体は、例えばシース熱電対である。測温体の後端側には、測定結果を出力するためのリード線１７８が連結されている。保護管は、測温体の損傷を防止するためのものであり、焼結鉱との接触及び衝突がある環境下で、望ましくは１ヶ月程度以上の長期の連続使用に耐え得る高温耐摩耗性のある材質から形成され、例えば鉄系合金を用いて形成される。また、保護管の厚みは、測温体の感温部のある先端部分が、他の部分よりも薄くなっている。保護管による測温体の感度への影響を低減するためである。
【００３１】
水平孔１６４ａ、１６４ｂの他端側は、ｚ軸方向に延びる垂直孔１６５に連結されている。垂直孔１６５の一端側は、基台部１６１の底面において外部と連通している。垂直孔１６５は、水平孔１６４ａ、１６４ｂ及び外部と連通し、水平孔１６４ａ、１６４ｂに挿入された温度測定部１７０のリード線１７８をパレット温度計１５０（筐体１６０）の外部へ引き出すための孔である。尚、パレット温度計１５０で測定された温度データは、焼結機の外部に、データ伝送部（図示せず）により取り出される。
【００３２】
また、筐体１６０には、水平孔１６４ａ、１６４ｂの一部と交わるようにｘ軸方向に延びる貫通孔１６６ａ、１６６ｂが形成される。貫通孔１６６ａ、１６６ｂには、水平孔１６４ａ、１６４ｂに挿入された温度測定部１７０ａ、１７０ｂを固定する固定部材１８０ａ、１８０ｂが挿入されている。温度測定部１７０ａ、１７０ｂの、筐体１６０から突出する部分には、焼結原料の荷重から多大な荷重がかかる。このとき、温度測定部１７０ａ、１７０ｂが筐体１６０から脱落しないように、筐体１６０と、温度測定部１７０ａ、１７０ｂとが固定部材１８０ａ、１８０ｂによって固定されるようにしている。
【００３３】
＜焼結プロセス操業監視装置＞
図４は、焼結プロセス操業監視装置の機能的な構成の一例を示すブロック図である。焼結プロセス操業監視装置４００は、例えば、ＣＰＵ、ＲＯＭ、ＲＡＭ、ＨＤＤ、及び各種のインターフェースを備えたコンピュータ装置を用いることにより実現することができる。
【００３４】
図４に示すように、焼結プロセス操業監視装置４００には、パレット温度計１５０で測定された焼結層内の温度データが、パレット温度計１５０から直接又は間接的に入力される。また、焼結プロセス操業監視装置４００には、プロセスコンピュータ４５０から、焼結機１００の操業データが入力される。また、焼結プロセス操業監視装置４００は、入力装置４６０からユーザによる指示を受け付ける。また、焼結プロセス操業監視装置４００は、後述するようにして計算された「焼結機１００における焼結層内の温度分布」等を表示装置４７０に表示する。以下、焼結プロセス操業監視装置４００の各機能について説明する。
【００３５】
［温度データ収集部４０１］
温度データ収集部４０１は、パレット温度計１５０で測定された焼結層内の温度の時系列データが上記のデータ伝送部経由で入力され、ＨＤＤ等に記憶する。ここで、本実施形態では、温度データ収集部４０１は、どのパレット温度計１５０のどの温度測定部１７０で測定された温度データであるのかが識別できるように、各温度データを記憶する。例えば、パレット温度計１５０毎、温度測定部１７０毎に温度データを記憶するデータベースに各温度データを記憶することができる。また、本実施形態では、温度データ収集部４０１は、パレット温度計１５０ａ〜１５０ｆが、原料給鉱部（端）から焼結鉱排鉱部（端）まで移動したときに得られた温度データを１単位としてＨＤＤ等に記憶する。ここでは、点火炉１３０の真下の位置を原料給鉱部とする。以下の説明では、パレット温度計１５０で測定された焼結層内の温度データを必要に応じて「測定温度データ」と称する。本実施形態では、６個のパレット温度計１５０ａ〜１５０ｆの夫々が２つの温度測定部１７０ａ、１７０ｂを有しているので、１２個の測定温度データ群が得られる。
【００３６】
［ゾーン解析部４０２］
ゾーン解析部４０２は、パレット温度計１５０ａ〜１５０ｆで測定された測定温度データに基づいて、焼結層の各ゾーンの始点温度θp,is［℃］、終点温度θp,ie［℃］、ゾーン通過時間ｔi［ｓｅｃ］、平均パレット速度ｖp［ｍ／ｓｅｃ］、及びゾーン長ｌi［ｍ］等を求める。ここで、ｉは、ゾーンを識別するものであり、１〜５の整数となる。
ここで、ゾーン通過時間ｔiとは、パレット温度計１５０ａ〜１５０ｆ（の測温体）が各ゾーンを通過するのに要した時間をいう。平均パレット速度ｖpとは、パレット温度計１５０ａ〜１５０ｆ（の測温体）が原料給鉱部から焼結鉱排鉱部まで移動している間のパレット１１０の平均速度をいう。ゾーン長ｌiは、パレット温度計１５０ａ〜１５０ｆの測温体が配置されている高さにおける、各ゾーンの機長方向の長さをいう。
尚、以後の説明において、「パレット温度計１５０ａ〜１５０ｆで測定された測定温度データ」をまとめて「パレット温度計１５０で測定された測定温度データ」として説明するが、本実施形態では、パレット幅方向（ｘ軸方向）に複数配置するパレット温度計１５０ａ〜１５０ｆにおいて、それぞれ個別に同じ処理を実施する。
【００３７】
本実施形態では、ゾーン解析部４０２は、焼結機１００で焼結されている焼結層を５つのゾーン（ｉ＝１〜５）に分けるようにしている。以下に、これらのゾーンについて説明する。
図５は、焼結層を区画するゾーンの一例を示す図である。
図５において、焼結層は、焼結機の機長方向に、原料給鉱部から焼結鉱排鉱部までが５つのゾーンI〜Vにより区画される。ゾーンIは、初期原料帯であり、焼結層に何も反応が起こっていない領域である。ゾーンIIは、水分凝縮帯であり、焼結層に含まれる水分が凝縮して気化する領域である。ゾーンIIIは、乾燥帯であり、水分が気化した焼結層が乾燥する領域である。ゾーンIVは、赤熱帯であり、乾燥後の焼結層内で燃料が着火し、焼結原料が焼結する領域である。ゾーンVは、焼結完了帯であり、燃料が全て燃え尽きて焼結反応が完了後に空冷され、焼結鉱が生成されている領域である。ｉ＝１〜５は、それぞれゾーンI〜Vに対応するものである。言い換えれば、あるパレット上の焼結層は焼結反応の進行につれて、当初は全体がゾーンI（初期原料帯）に属し、その後ゾーンII（水分凝縮帯）、III（乾燥帯）、IV（赤熱帯）が上部から下部へ移動していき、焼結の完了によって全体がゾーンV（焼結完了帯）で占められる。図５に示すように、赤熱帯ゾーンIVの底部の終点（焼結完了帯ゾーンVの始点）が焼結完了点（ＢＴＰ）となる。
【００３８】
図６は、焼結層のある高さ位置における測定温度データの一例を概念的に示す図である。
〔始点温度θp,is、終点温度p,ie、ゾーン通過時間ｔi〕
（ゾーンI）
まず、パレット温度計１５０の測定結果から、ゾーンI（初期原料帯）に相当する温度の波形（温度と時間との関係）を確認することができなかった。そこで、本実施形態では、ゾーンIを省略して取り扱うこととする。すなわち、パレット温度計１５０の測定温度データθp(ｔa)の初期値（時間ｔa＝０における測定温度θp(０)）を、ゾーンII（水分凝縮帯）の始点温度θp,2sとする。すなわち、ゾーン解析部４０２は、ゾーンIの始点温度θp,1sと、ゾーンIの終点温度θp,1eとを以下の（１）式のように設定する。
θp,1s＝θp,1e＝θp(０)  ・・・（１）
【００３９】
また、ゾーン解析部４０２は、ゾーンIのゾーン通過時間ｔ1を以下の（２）式のように設定する。
ｔ1＝０  ・・・（２）
尚、以下に示す時間ｔaは、このゾーン通過時間ｔ1＝０と設定した時間をｔa＝０として基準とした時間である。
【００４０】
（ゾーンII）
前述したように、ゾーン解析部４０２は、ゾーンIの終点温度θp,1eを、ゾーンII（水分凝縮帯）の始点温度θp,2sとして設定する。すなわち、ゾーン解析部４０２は、ゾーンIIの始点温度θp,2sを以下の（３）式のように設定する。
θp,2s＝θp,1e  ・・・（３）
また、ゾーン解析部４０２は、パレット温度計１５０の測定温度データθp(ｔa)の最小値θp,minを、ゾーンIIの終点温度θp,2eとして設定する。すなわち、ゾーン解析部４０２は、ゾーンIIの終点温度θp,2eを以下の（４）式のように設定する。
θp,2e＝θp,min＝θp(ｔa,min)  ・・・（４）
ここで、時間ｔa,min［ｓｅｃ］は、パレット温度計１５０の測定温度データθp(ｔa)が最小値θp,minとなったときの時間である。
さらに、ゾーン解析部４０２は、ゾーンIIのゾーン通過時間ｔ2を、以下の（５）式のように設定する。
ｔ2＝ｔa,min  ・・・（５）
【００４１】
（ゾーンIII）
ゾーン解析部４０２は、ゾーンIIの終点温度θp,2eを、ゾーンIII（乾燥帯）の始点温度θp,3sとして設定する。すなわち、ゾーン解析部４０２は、ゾーンIIIの始点温度θp,3sを以下の（６）式のように設定する。
θp,3s＝θp,2e  ・・・（６）
また、ゾーン解析部４０２は、パレット温度計１５０の測定温度データθp(ｔa)が、予め定められている焼結原料の着火温度θc,f［℃］になったときの温度を、ゾーンIIIの終点温度θp,3eとして設定する。すなわち、ゾーン解析部４０２は、ゾーンIIIの終点温度θp,3eを以下の（７）式のように設定する。
θp,3e＝θc,f＝θp(ｔa,cf)  ・・・（７）
ここで、時間ｔa,cf［ｓｅｃ］は、パレット温度計１５０の測定温度データθp(ｔa)が、焼結原料の着火温度θc,fとなったときの時間である。
さらに、ゾーン解析部４０２は、ゾーンIIIのゾーン通過時間ｔ3を、以下の（８）式に従って計算し設定する。
ｔ3＝ｔa,cf−ｔa,min  ・・・（８）
ただし、ｔa,min＜ｔa,cf＜ｔa,maxの関係が成立するものとする（ｔa,max［ｓｅｃ］は、後述するようにパレット温度計１５０の測定温度データθp（ｔa）が最大値θp,maxとなったときの時間である）。
【００４２】
（ゾーンIV）
ゾーン解析部４０２は、ゾーンIIIの終点温度θp,3eを、ゾーンIV（赤熱帯）の始点温度θp,4sとして設定する。すなわち、ゾーン解析部４０２は、ゾーンIVの始点温度θp,4sを以下の（９）式のように設定する。
θp,4s＝θp,3e  ・・・（９）
また、ゾーン解析部４０２は、パレット温度計１５０の測定温度データθp(ｔa)の最大値θp,maxを、ゾーンIVの終点温度θp,4eとして設定する。すなわち、ゾーン解析部４０２は、ゾーンIVの終点温度θp,4eを以下の（１０）式のように設定する。
θp,4e＝θp,max＝θp(ｔa,max)  ・・・（１０）
ここで、時間ｔa,max［ｓｅｃ］は、パレット温度計１５０の測定温度データθp(ｔa)が最大値θp,maxとなったときの時間である。
さらに、ゾーン解析部４０２は、ゾーンIVのゾーン通過時間ｔ4を、以下の（１１）式に従って計算し設定する。
ｔ4＝ｔa,max−ｔa,cf  ・・・（１１）
【００４３】
（ゾーンV）
後述するように、本実施形態では、ゾーンＩ〜IV（初期原料帯、水分凝縮帯、乾燥帯、及び赤熱帯）の温度分布を求めるので、ゾーンV（焼結完了帯）については必ずしも求める必要はない。しかしながら、本実施形態は、ゾーンV（焼結完了帯）の温度分布を求めることを排除するものではない。また、平均パレット速度ｖp等を求める際には、ゾーンVの情報（ｔa,5等）を用いることができる。したがって、ゾーンVの温度分布を求めるための情報の導出についても説明する。
【００４４】
ゾーン解析部４０２は、ゾーンIVの終点温度θp,4eを、ゾーンVの始点温度θp,5sとして設定する。すなわち、ゾーン解析部４０２は、ゾーンVの始点温度θp,5sを以下の（１２）式のように設定する。
θp,5s＝θp,4e  ・・・（１２）
また、ゾーン解析部４０２は、パレット温度計１５０の測定温度データθp(ｔa)の極小値θ´p,minを、ゾーンVの終点温度θp,5eとして設定する。すなわち、ゾーン解析部４０２は、ゾーンVの終点温度θp,5eを以下の（１３）式のように設定する。
θp,5e＝θ´p,min＝θp(ｔ´a,min)  ・・・（１３）
ただし、ｔa,max＜ｔ´a,min＜ｔa,eの関係が成立するものとする。
ここで、時間ｔ´a,min［ｓｅｃ］は、図６に例示するパレット温度計１５０の測定温度データθp(ｔa)が極小値θ´p,minとなったときの時間である。
さらに、ゾーン解析部４０２は、ゾーンVのゾーン通過時間ｔ5を、以下の（１４）式に従って計算し設定する。
ｔ5＝ｔ´a,min−ｔa,max  ・・・（１４）
【００４５】
〔平均パレット速度ｖp〕
ゾーン解析部４０２は、以下の（１５）式に従って、平均パレット速度ｖpを計算する。
ｖp＝ｌa,5／ｔa,5  ・・・（１５）
ここで、ｌa,5［ｍ］は、パレット温度計１５０（の測温体）が原料給鉱部（端）から焼結鉱排鉱部（端）まで移動する間の機長方向における長さであり、予め設定されるものである。また、ｔa,5［ｓｅｃ］は、パレット温度計１５０（の測温体）が前記機長方向の長さｌa,5を移動している間の時間であり、各ゾーンのゾーン通過時間ｔiの和として、以下の（１６）式で求められる。
ｔa,5＝ｔ1＋ｔ2＋ｔ3＋ｔ4＋ｔ5  ・・・（１６）
後述するように本実施形態では、ゾーンＩ〜IV（初期原料帯、水分凝縮帯、乾燥帯、及び赤熱帯）の温度分布を求めることにより、ＢＴＰを導出する。よって、平均パレット速度ｖpを以下の（１５´）式、（１６´）式に従って計算するようにしてもよい。
ｖp＝ｌa,4／ｔa,4  ・・・（１５´）
ｔa,4＝ｔ1＋ｔ2＋ｔ3＋ｔ4  ・・・（１６´）
【００４６】
〔ゾーン長ｌi〕
ゾーン解析部４０２は、以下の（１７）式に従って、各ゾーンのゾーン長ｌiを計算する。
ｌi＝ｔi×ｖp  ・・・（１７）
【００４７】
尚、本実施形態では、ゾーン解析部４０２は、ゾーンIのゾーン位置ｌa,1を０（ゼロ）に設定する。ここで、ゾーンｉのゾーン位置ｌa,iは、原点（点火炉１３０から焼結原料が点火される位置（原料給鉱部））からゾーンｉの終点までの機長方向における距離を示すものである。
具体的にゾーン解析部４０２は、ソーンｉのゾーン位置ｌa,i［ｍ］を、以下の（１８）式で計算する。
ｌa,i＝Σｌm  （Σはｍ＝１からｉ（１≦ｉ≦５）までのｌmの総和）・・・（１８）
ゾーン解析部４０２は、以上の処理を、全てのパレット温度計１５０ａ〜１５０ｆについて個別に行う。
【００４８】
［規定データ設定部４０３］
規定データ設定部４０３は、ユーザによる入力装置４６０の操作に基づいて、各種の定数と、必要に応じてこれらの定数の温度依存関数とを設定（ＨＤＤ等に記憶）する。これらの定数や温度依存関数は、公知の文献に示されているものや、過去の知見から得られたものを採用することができる。
【００４９】
規定データ設定部４０３は、焼結鉱粒子（焼結層を構成する粒子）の密度ρp［ｋｇ／ｍ3］を設定する。本実施形態では、焼結鉱粒子の密度ρpとして、ゾーンI〜Vで共通の値を用いるようにしているが、ゾーン毎に個別の値を用いても構わない。
【００５０】
また、規定データ設定部４０３は、焼結鉱粒子の直径Ｄp,i［ｍ］を設定する。本実施形態では、焼結鉱粒子の直径Ｄp,iとして、ゾーンI〜V毎に値を設定するようにしているが、ゾーンI〜Vで共通の値を用いても構わない。
図７は、焼結層内の焼結鉱の粒子充填層を概念的に示す図である。図７に示すように、本実施形態では、解析に際し、同じゾーン内では、焼結鉱粒子８０２のそれぞれが同じ直径Ｄp,iを有する球であるとしているが、焼結鉱粒子それぞれを表現する立体形状として他の多面体や曲面体を用いても構わない。尚、図７において、上から下に向かう矢印８０１は、空気の流れを示している。
【００５１】
また、規定データ設定部４０３は、燃料（粉コークス）のモル重量Ｗc［ｋｇ／ｍｏｌ］、密度ρc［ｋｇ／ｍ3］及び燃焼熱量Ｑc［Ｊ／ｍｏｌ］をそれぞれ設定する。
また、規定データ設定部４０３は、焼結原料の着火温度θc,f［℃］を設定する。本実施形態では、焼結原料の着火温度θc,fは、（７）式に示したように、ゾーンIIとゾーンIIIとの境界を規定するための温度として用いる（図６も参照）。尚、ゾーンIIとゾーンIIIとの境界を規定するための温度には、必要に応じて別の温度を用いても構わない。
【００５２】
また、規定データ設定部４０３は、水のモル重量Ｗw［ｋｇ／ｍｏｌ］及び蒸発潜熱Ｑw［Ｊ／ｍｏｌ］をそれぞれ設定する。
【００５３】
さらに、本実施形態では、規定データ設定部４０３は、焼結鉱粒子の比熱Ｃp［ｋＪ／（ｋｇ・Ｋ）］、水の密度ρw［ｋｇ／ｍ3］、空気の密度ρg［ｋｇ／ｍ3］、定圧比熱Ｃpg［ｋＪ／（ｋｇ・Ｋ）］、粘性率μg［Ｐａ・ｓ］及び熱伝導率λg［Ｗ／（ｍ・Ｋ）］のそれぞれを、温度θ［℃］に依存する特性を有するものとして設定する。例えば、これらのそれぞれを、温度θ［℃］を独立変数とする低次の線形多項式として設定する。これらの線形多項式は、例えば、公知の文献の物性値のデータ表の値を近似関数化することにより得られる。
尚、これらのそれぞれの温度依存特性を非線形多項式で近似関数化して用いても構わないし、公知の文献の物性値のデータ表をそのまま用いたデータテーブル参照方式によって用いても構わない。また、これらのそれぞれのいくつか又は全ての温度依存特性は小さいものとして、ある値で代表させて用いても構わない。
また、規定データ設定部４０３は、例えば、ユーザによる入力装置４６０の操作の内容に基づいて、機長、すなわち原料給鉱部（端）から焼結鉱排鉱部（端）までの機長方向における長さla,5［ｍ］を設定する。
さらに、規定データ設定部４０３は、例えば、ユーザによる入力装置４６０の操作の内容に基づいて、基台部１６１の上面の位置と、温度測定部１７０ａ、１７０ｂ（の中心）のｚ軸方向における位置を設定する。尚、本実施形態では、温度測定部１７０ｂ（の中心）から温度測定部１７０ａ（の中心）までのｚ軸方向における距離が１５０［ｍｍ］であるとする。
【００５４】
［初期値設定部４０４］
初期値設定部４０４は、ユーザによる入力装置４６０の操作に基づいて、各ゾーンI〜Vについて、空気流入温度θgi,i［℃］、形状係数φs,i［−］の初期値θgi,i,0、φs,i,0を設定（ＨＤＤ等に記憶）する。これらの値としては、過去の知見から得られた値を採用する。ここで、空気流入温度θgi,iは、各ゾーンI〜Vの焼結層内（のパレット温度計１５０の測温体が配置されている位置）に流入する空気の温度である。また、形状係数φs,iは、焼結鉱粒子の直径Ｄp,iに掛ける係数であり０〜１の間の値を有する。前述したように、本発明者らは、焼結層内の機長方向における温度分布を解析した結果と実測した結果とを比較した結果、焼結鉱粒子の伝熱に寄与する表面積が変化しないまま加熱・冷却が進行する場合と、変化しながら加熱・冷却が進行する場合とがあることを見出した。このような傾向は焼結完了帯で顕著となった。このことから、本発明者らは、焼結鉱粒子が融着を起こすこと等により、焼結鉱粒子の形状や焼結層内の粒子充填構造が変化し、図７に示したそれぞれの焼結鉱粒子の表面積の全てが伝熱に寄与しないと仮定した。そこで、焼結鉱粒子の表面積のうち、伝熱に寄与する面積を導出すべく、本発明者らは、本実施形態において、形状係数φs,iを導入した。
【００５５】
［操業データ収集部４０５］
操業データ収集部４０５は、プロセスコンピュータ４５０から、圧力損失（ブロア吸引圧）ΔＰ［Ｐａ］と、焼結層の層厚Ｌ［ｍ］と、初期水分比Ｋw0［−］と、初期燃料比Ｋc0［−］が入力され、ＨＤＤ等に記憶する。ここで、水分比Ｋwとは、焼結機１００の前処理の造粒過程（造粒機、ミキサー）で投入する水分の焼結原料に対する重量比［ｋｇ-ｗａｔｅｒ／ｋｇ-焼結原料］である。また、燃料比Ｋcとは、焼結原料に対する燃料（粉コークス）の重量比［ｋｇ-ｃｏｋｅ／ｋｇ-焼結原料］である。
【００５６】
［モデル格納部４０６］
本発明者らは、パレット温度計１５０の測定温度データから、焼結層内の機長方向における温度が指数関数的に変化するという知見を得た。また、伝熱工学における粒子体の熱収支を表す時間微分方程式の解析解の一般解には指数関数の項が含まれることが知られている。これらのことを考慮して、本発明者らは、伝熱工学における粒子体の熱収支を表す時間微分方程式の解析解により、焼結層内の機長方向における温度の変化を求めることを見出した。また、焼結層を構成する焼結鉱粒子を粒子体とみなした場合、粉体工学（粒子充填層工学）における知見から、本発明者らは、粒子充填層における流動抵抗、すなわち圧力損失ΔＰ［Ｐａ］を焼結層の層厚Ｌ［ｍ］で除した圧力変化率ΔＰ／Ｌ［Ｐａ／ｍ］を、焼結層内（パレット温度計１５０の測温体が配置されている位置）に流入する空気流速ｖg,i［ｍ／ｓｅｃ］に関する多項式関数ｆ1で表現することを見出した。
【００５７】
〔粉体工学における粒子充填層内の圧力変化率と流速に関する関係式〕
粉体工学における粒子充填層内の圧力変化率と流速に関する関係は、以下の（１９）式で表される。
【００５８】
【数１】
0005729251.tif 000002
【００５９】
（１９）式において、カンマの後のｉは、前述したようにゾーンを識別するものである（ｉ＝１〜５はそれぞれゾーンI〜Vに対応する）。ｖg,iは、ゾーンｉの焼結層内（パレット温度計１５０の測温体が配置されている位置）に流入する空気の流速（空気流速）［ｍ／ｓｅｃ］である。ここで、Ｃkは圧力変化率ΔＰ／Ｌ［Ｐａ／ｍ］を表現する多項式関数の各項の係数であり、ｋは多項式関数を構成する空気流速ｖg［ｍ／ｓｅｃ］の次数で、ｎは多項式関数の最高次数で、焼結鉱の粒子充填層の各ゾーンを模擬した粒子充填層の流動抵抗実験結果や公知の文献に記載のデータ表の値を多項式関数近似することによって、予め決定しておく。本実施形態では、（１９）式で例示した如く、圧力変化率ΔＰ／Ｌ［Ｐａ／ｍ］を空気流速ｖg,iの線形多項式関数として表現し、例えばNewton法による収束計算法を用いて、空気流速ｖg,iをゾーン毎に計算するようにしている。次数nが低次であれば、解析的に空気流速ｖg,iを求解してもよい。また、（１９）式は圧力変化率ΔＰ／Ｌ［Ｐａ／ｍ］を空気流速ｖg,iの線形多項式関数として表現した例を示しているが、（１９）式のような線形多項式関数でなく、圧力変化率ΔＰ／Ｌ［Ｐａ／ｍ］を空気流速ｖg,iの非線形多項式関数で表現しても構わない。この場合は、例えば二分法のような収束計算法を用いて、空気流速ｖg,iをゾーン毎に計算することができる。或いは、粒子充填層内の圧力変化率と流速の関係式として一般に知られているErgunの式に前記定義の形状係数φs,iを付加して用いた式で空気流速ｖg,iをゾーン毎に計算しても構わない。尚、Ergunの式は、例えば「Ergun, S. : Fluid Flow through Packed Columns, Chem. Eng. Prog., Vol.48(1952), pp.89-94」に記載されている。
【００６０】
〔伝熱工学に関する式；焼結鉱粒子の熱収支を表す時間微分方程式〕
本実施形態では、焼結鉱粒子と空気の対流熱伝達における焼結鉱粒子の熱収支を表す時間微分方程式を、以下の（２０）式で表すものとする。
【００６１】
【数２】
0005729251.tif 000003
【００６２】
（２０）式においても、カンマの後のｉは、ゾーンを識別するものである（ｉ＝１〜５はそれぞれゾーンI〜Vに対応する）。また、θp,i(ｔ)は、各ゾーンｉの始点を通過した時間を基準（すなわち、t＝０）とする時間ｔにおけるゾーンｉの焼結層内（のパレット温度計１５０の測温体が配置されている位置）の温度［℃］である。また、ｈi［Ｗ／（ｍ2・Ｋ）］は、焼結鉱粒子まわりの対流熱伝達率であり、以下の（２１）式で表される。
【００６３】
【数３】
0005729251.tif 000004
【００６４】
（２１）式に示すＮｕiは、焼結鉱粒子まわりの平均ヌッセルト数［−］である。また、（２２）式に示すＮｕiは、例えば、焼結鉱粒子を球形と仮定したときの、焼結鉱球形粒子まわりの流体による対流熱伝達特性を、固体球まわりの平均ヌッセルト数として表現したものを例示したものである。一般に、伝熱工学における対流熱伝達について、固体球まわりの流体による対流熱伝達特性、すなわち固体球まわりの平均ヌッセルト数は、固体球まわりの流体の流速ｖを固体球径Ｄp等で無次元化表現するレイノルズ数Ｒｅ［−］と、流体の物性値である粘性率μ、定圧比熱Ｃp、熱伝導率λを総括して無次元化表現するプラントル数Ｐｒ［−］との関数で表現できることが知られている。
そこで、本発明者らは、焼結層を構成する焼結鉱粒子まわりの対流熱伝達特性も同様にして、焼結層内（パレット温度計１５０の測温体が配置されている位置）に流入する空気の流速ｖg,i［ｍ／ｓｅｃ］に関するレイノルズ数Ｒｅiと当該空気の物性値に関するプラントル数Ｐｒiの関数ｆ2で表現できることを見出し、焼結鉱の粒子充填層の各ゾーンを模擬した粒子充填層の対流熱伝達特性実験結果や公知の文献に記載のデータ表の値から当該関数の係数を同定し、予め決定しておく。本実施形態では、予め係数を決定しておいた（２２）式を用いてＮｕiを決定し、その値を（２１）式に代入することによって、ｈiの値を求めるようにしている。
尚、（２２）式には、Ｎｕiの必要に応じて任意の関数形を用いることが可能である。或いは、伝熱工学における固体球まわりの流体による対流熱伝達特性、すなわち固体球まわりの平均ヌッセルト数の特性式として一般に知られているRanz-Marchallの式や、当該式中の係数の値を焼結鉱の粒子充填層の各ゾーンを模擬した粒子充填層の対流熱伝達特性実験結果や公知の文献に記載のデータ表の値から必要に応じて同定した式を用いても構わない。尚、Ranz-Marchallの式は、例えば「Ranz, W. E. and Marsha11, W. R. : Evaporation from Drops I＆II, Chem. Eng. Prog., Vol.48(1952), pp.141-146, 173-180」に記載されている。
また、Ｖp,iは、ゾーンｉにおける焼結鉱粒子の体積［ｍ3］であり、Ａp,iは、ゾーンｉにおける焼結鉱粒子の伝熱表面積［ｍ2］である。本実施形態では、焼結鉱粒子を球形粒子とし、前記定義の形状係数φs,iを用いて伝熱表面積を定義して用いることとする。このとき、焼結鉱粒子の体積Ｖp,i及び伝熱表面積Ａp,iは、それぞれ以下の（２３）式、（２４）式で表され、焼結鉱球形粒子まわりの空気のレイノルズ数Ｒｅi及びプラントル数Ｐｒiは、（Ａ）式、（Ｂ）式で表される。
【００６５】
【数４】
0005729251.tif 000005
【００６６】
Ｒｅi＝（ｖg,i・Ｄp,i・ρg,i）／μg,i  ・・・（Ａ）
Ｐｒi＝（μg,i・Ｃp,i）／λg,i  ・・・（Ｂ）
【００６７】
また、Ｒwは、水分気化速度［ｍｏｌ／ｓｅｃ］であり、Ｒcは、燃料燃焼速度［ｍｏｌ／ｓｅｃ］である。
本実施形態では、（２０）式の微分方程式を、各ゾーンの始点温度を直前のゾーンの終点温度（最終到達温度）とする境界条件で解くことにより、各ゾーンｉの焼結層内の温度θp,i(ｔ)を解析解としてそれぞれ求めることとする。この解析解において、以下の（２５）式に示す時定数Ｔi［ｓｅｃ］を定義する。また、（２５）式に（２４）式を代入して、形状係数φs,iについて解くと、以下の（２６）式に示す形状係数φs,iを求めることができる。
【００６８】
【数５】
0005729251.tif 000006
【００６９】
〔水分凝縮帯（ゾーンII）における焼結鉱粒子の水分の質量収支を表す時間微分方程式〕
本実施形態では、焼結鉱粒子の水分の質量収支を表す時間微分方程式を、以下の（２７）式で表すものとする。
【００７０】
【数６】
0005729251.tif 000007
【００７１】
ここで、Ｋw(ｔ)［−］は、ゾーンIIの始点を通過した時間を基準とする時間ｔにおける水分比である。水分比の定義は前述した通りである。水分凝縮帯であるゾーンIIの終点（ｔ＝ｔ2）では水分は全て蒸発しているので、Ｋw(ｔ2)＝０となる。この条件のもとで、例えば水分気化速度ＲwがゾーンIIにおいて一定であると仮定して、（２７）式を水分比Ｋw(ｔ)について解くことにより解析解、すなわちゾーンIIにおける水分比の時間変化を求めることができる。この場合、解析解の一般解は１次関数である。或いは、水分気化速度ＲwがゾーンII内で一定であると仮定せずに任意の関数形を設定して（２７）式を水分比Ｋw(ｔ)について解いても構わない。尚、（２７）式のこの他の変数は前述した通りである。
【００７２】
〔赤熱帯（ゾーンIV）における焼結鉱粒子の燃料の質量収支を表す時間微分方程式〕
本実施形態では、焼結鉱粒子の燃料の質量収支を表す時間微分方程式を、以下の（２８）式で表すものとする。
【００７３】
【数７】
0005729251.tif 000008
【００７４】
ここで、Ｋc(ｔ)［−］は、ゾーンIVの始点を通過した時間を基準とする時間ｔにおける燃料比である。燃料比の定義は前述した通りである。赤熱帯であるゾーンIVの終点（ｔ＝ｔ4）では燃料は全て燃焼しているので、Ｋc(ｔ4)＝０となる。この条件のもとで、例えば燃料燃焼速度ＲcがゾーンIVにおいて一定であると仮定して、式（２８）を燃料比Ｋc(ｔ)について解くことにより解析解、すなわちゾーンIVにおける燃料比の時間変化を求めることができる。この場合、解析解の一般解は１次関数である。或いは、燃料燃焼速度ＲcがゾーンIV内で一定であると仮定せずに任意の関数形を設定して（２８）式を燃料比Ｋc(ｔ)について解いても構わない。尚、（２８）式のこの他の変数は前述した通りである。
【００７５】
尚、モデル格納部４０６には、前述した（１９）式〜（２８）式そのものの他に、（１９）式〜（２８）式に基づく各ゾーンI〜V毎の式、例えば各ゾーンI〜Vにおける（２０）式の微分方程式を前記記載の境界条件でθp,i(ｔ)について解いた式、Ｋw(ｔ2)＝０として（２７）式の微分方程式をＫw(ｔ)について解いた式、Ｋc(ｔ4)＝０として（２８）式の微分方程式をＫc(ｔ)について解いた式等を記憶しておいてもよい。
【００７６】
［オンライン計算処理部４０７］
［［概要］］
オンライン計算処理部４０７で行われる処理の概要の一例を説明する。
まず、オンライン計算処理部４０７は、ゾーン解析部４０２で解析された結果と、操業データ収集部４０５で収集されたデータと、規定データ設定部４０３及び初期値設定部４０４で設定されたデータとを取得する。オンライン計算処理部４０７は、取得したデータを、モデル格納部４０６に記憶されたモデルに適用して、「温度測定部１７０ａ、１７０ｂが配置されている位置での測定温度データに基づく関数（モデル式）θp,iA(ｔiA)、θp,iB(ｔiB)を定めるパラメータ」を、各ゾーンII〜IVのそれぞれにおいて同定する。
【００７７】
次に、オンライン計算処理部４０７は、同定したパラメータを使って、「『温度測定部１７０ａが配置されている位置での測定温度データに基づく関数θp,iA(ｔiA)』と『温度測定部１７０ｂが配置されている位置での測定温度データに基づく関数θp,iB(ｔiB)』」を、例えば、（２０）式で例示する時間微分方程式の解として、各ゾーンII〜IVのそれぞれにおいて求める。
【００７８】
次に、オンライン計算処理部４０７は、「『温度測定部１７０ａが配置されている位置での測定温度データに基づく関数θp,iA(ｔiA)上の点』と、当該点とゾーン内における位置が相互に対応する『温度測定部１７０ｂが配置されている位置での測定温度データに基づく関数θp,iB(ｔiB)上の点』」を、各ゾーンII〜IVのそれぞれにおいて抽出する。例えば、時間ｔの最大値と最小値を含み、且つ、それぞれの時間隔が略等間隔になる２０個の時間における温度を抽出する。そして、オンライン計算処理部４０７は、抽出した「温度と時間」の値から外挿又は内挿を行って、「測定温度データが得られていない『焼結層の高さ方向（ｚ軸方向）の位置』における温度データ」を、各ゾーンII〜IVのそれぞれにおいて導出する。以下の説明では、測定温度データが得られていない「焼結層の高さ方向（ｚ軸方向）の位置」における温度データを、必要に応じて、「導出対象温度データ」と称する。そして、オンライン計算処理部４０７は、当該位置における「導出対象温度データに基づく関数θp,iC(ｔiC)」を定めるパラメータを導出する。そして、オンライン計算処理部４０７は、導出したパラメータを使って、「測定温度データが得られていない『焼結層の高さ方向（ｚ軸方向）の位置』における導出対象温度データに基づく関数θp,iC(ｔiC)」を、各ゾーンII〜IVのそれぞれにおいて求める。
【００７９】
ここで、オンライン計算処理部４０７は、IIゾーンにおけるパラメータとして、空気流入温度θgi,2、水分気化速度Ｒw、対流熱伝達率ｈ2、時定数Ｔ2、形状係数φs,2及び空気流速ｖg,2を同定する。
また、オンライン計算処理部４０７は、IIIゾーンにおけるパラメータとして、空気流入温度θgi,3、対流熱伝達率ｈ3、時定数Ｔ3、形状係数φs,3及び空気流速ｖg,3を同定する。
また、オンライン計算処理部４０７は、IVゾーンにおけるパラメータとして、空気流入温度θgi,4、燃料燃焼速度Ｒc、対流熱伝達率ｈ4、時定数Ｔ4、形状係数φs,4及び空気流速ｖg,4を同定する。
尚、本実施形態では、オンライン計算処理部４０７は、Vゾーンについては、前述した処理を行わないが、Vゾーンについても処理する場合、Vゾーンにおけるパラメータとして、空気流入温度θgi,5、対流熱伝達率ｈ5、時定数Ｔ5、形状係数φs,5及び空気流速ｖg,5を同定する。
【００８０】
また、オンライン計算処理部４０７は、同一の高さ位置にある温度データ（測定温度データ、又は導出対象温度データ）であって、ゾーン内での位置が相互に対応する温度データから、（ｍ−１）次多項式関数を導出して内挿又は外挿を行い、焼結層のパレット幅方向（ｘ軸方向）の位置における温度を求める。ここで、ｍは、１つのパレット１１０において、パレット幅方向（ｘ軸方向）に配置されているパレット温度計１５０の数であり、本実施形態では６になる。
オンライン計算処理部４０７は、パレット温度計１５０ａ〜１５０ｆの温度測定部１７０ａ、１７０ｂで測定された測定温度データのそれぞれについて以上の処理を個別に行う。
以上の処理を、パレット温度計１５０を具備したパレット１１０が焼結機１００を１スイープする毎に逐次行うことにより、焼結層の任意の位置における温度（焼結層全体の３次元温度分布）を逐次得ることができる。
【００８１】
［［パラメータの同定］］
次に、前述したゾーンII〜IVにおける温度データ（測定温度データ及び導出対象温度データ）に基づく関数θp,2A(ｔ2A)、θp,2B(ｔ2B)、θp,2C(ｔ2C)、θp,3A(ｔ3A)、θp,3B(ｔ3B)、θp,3C(ｔ3C)、θp,4A(ｔ4A)、θp,4B(ｔ4B)、θp,4C(ｔ4C)を定めるパラメータを同定する処理の一例を説明する。尚、ここではゾーンVにおける温度データに基づく関数θp,5A(ｔ5A)、θp,5B(ｔ5B)、θp,5C(ｔ5C)を定めるパラメータを同定する処理の一例も説明する。前述したように本実施形態では、ゾーンIの始点温度θp,1sと、ゾーンIの終点温度θp,1eとをゾーンII（水分凝縮帯）の始点温度θp,2sとし（（１）式を参照）、且つ、ゾーンIのゾーン通過時間ｔ1を０（ゼロ）としている（（２）式を参照）。したがって、ゾーンIにおける温度θp,1(ｔ)を定めるパラメータについては同定しない。
【００８２】
（ゾーンII）
オンライン計算処理部４０７は、規定データ設定部４０３で設定された「空気の密度ρg,2、空気の粘性率μg,2及び焼結鉱粒子の粒子径Ｄp,2」と、初期値設定部４０４で設定された「形状係数の初期値φs,2,0」と、操業データ収集部４０５で得られた「圧力損失ΔＰ及び焼結層の層厚Ｌ」とにより（１９）式からゾーンIIにおける空気流速ｖg,2を計算する。
【００８３】
次に、オンライン計算処理部４０７は、初期値設定部４０４で設定された「初期水分比Ｋw0」と、ゾーン解析部４０２で得られた「ゾーンIIのゾーン通過時間ｔ2」と、規定データ設定部４０３で設定された「水の密度ρw、及び水のモル重量Ｗw」と、（２３）式で計算される「ゾーンIIにおける焼結鉱粒子の体積Ｖp,2」とにより（２７）式に基づく式から水分気化速度Ｒwを計算する。
【００８４】
次に、オンライン計算処理部４０７は、前述したようにして計算した「ゾーンIIにおける空気流速ｖg,2」と、規定データ設定部４０３で設定された「空気の密度ρg,2、粘性率μg,2、熱伝導率λg,2及び定圧比熱Ｃp,2、及び焼結鉱粒子の直径Ｄp,2」とにより（２２）式からＮｕ2を求める。そして、オンライン計算処理部４０７は、（２１）式から対流熱伝達率ｈ2を計算する。
また、オンライン計算処理部４０７は、初期値設定部４０４で設定された「形状係数の初期値φs,2,0」と、規定データ設定部４０３で設定された「ゾーンIIにおける焼結鉱粒子の直径Ｄp,2」とにより（２４）式からゾーンIIにおける焼結鉱粒子の伝熱表面積Ａp,2を求める。
【００８５】
そして、オンライン計算処理部４０７は、ゾーン解析部４０２で得られた「ゾーンIIの終点温度θp,2e」と、規定データ設定部４０３で設定された「水分蒸発潜熱Ｑw」と、前述したようにして計算された「水分気化速度Ｒw」、「対流熱伝達率ｈ2」及び「ゾーンIIにおける焼結鉱粒子の伝熱表面積Ａp,2」とにより（２０）式、（２７）式に基づく式からゾーンIIの空気流入温度θgi,2の中間計算値θgi,2,tmpを求める。
【００８６】
次に、オンライン計算処理部４０７は、（２５）式より、ゾーンIIにおける時定数Ｔ2を計算する。そして、オンライン計算処理部４０７は、（２６）式より、形状係数φs,2の中間計算値φs,2,tmpを求める。そして、オンライン計算処理部４０７は、形状係数φs,2の中間計算値φs,2,tmpと、それまでに設定されている形状係数φs,2（最初は初期値φs,2,0）との差の絶対値が収束判定基準値ＥＰＳ以下であるか否かを判定する。
【００８７】
この判定の結果、形状係数φs,2の中間計算値φs,2,tmpと、それまでに設定されている形状係数φs,2との差の絶対値が収束判定基準値ＥＰＳ以下でない場合、オンライン計算処理部４０７は、それまでに設定されている形状係数φs,2を、中間計算値φs,2,tmpに置き換える（最初は初期値φs,2,0を中間計算値φs,2,tmpに置き換える）。そして、オンライン計算処理部４０７は、この中間計算値φs,2,tmpに置き換えた形状係数φs,2を用いて、前述した「空気流速ｖg,2」、「水分気化速度Ｒw」、「空気流入温度θgi,2の中間計算値θgi,2,tmp」、「形状係数φs,2の中間計算値φs,2,tmp」を再度求める。オンライン計算処理部４０７は、このような処理を、形状係数φs,2の中間計算値φs,2,tmpと、それまでに設定されている形状係数φs,2との差の絶対値が収束判定基準値ＥＰＳ以下になるまで繰り返し行う。
【００８８】
そして、形状係数φs,2の中間計算値φs,2,tmpと、それまでに設定されている形状係数φs,2との差の絶対値が収束判定基準値ＥＰＳ以下になると、オンライン計算処理部４０７は、空気流入温度θgi,2の中間計算値θgi,2,tmpと、それまでに設定されている空気流入温度θgi,2（最初は初期値θgi,2,0）との差の絶対値が収束判定基準値ＥＰＳ以下であるか否かを判定する。この判定の結果、空気流入温度θgi,2の中間計算値θgi,2,tmpと、それまでに設定されている空気流入温度θgi,2との差の絶対値が収束判定基準値ＥＰＳ以下でない場合、オンライン計算処理部４０７は、それまでに設定されている空気流入温度θgi,2を、中間計算値θgi,2,tmpに置き換える（最初は初期値θgi,2,0を中間計算値θgi,2,tmpに置き換える）。そして、オンライン計算処理部４０７は、この中間計算値θgi,2,tmpを用いて、前述した「空気流速ｖg,2」、「水分気化速度Ｒw」、「空気流入温度θgi,2の中間計算値θgi,2,tmp」、「形状係数φs,2の中間計算値φs,2,tmp」の計算を繰り返し行う。
【００８９】
そして、ゾーンIIの空気流入温度θgi,2の中間計算値θgi,2,tmpと、それまでに設定されているゾーンIIの空気流入温度θgi,2（最初は初期値θgi,2,0）との差の絶対値が収束判定基準値ＥＰＳ以下になると、オンライン計算処理部４０７は、そのときに得られた「空気流速ｖg,2」、「水分気化速度Ｒw」、「対流熱伝達率ｈ2」、「空気流入温度θgi,2」、「時定数Ｔ2」、「形状係数φs,2」を最終計算値として確定し、これらを同定パラメータとしてパラメータ蓄積部４０８に記憶する。
【００９０】
(ゾーンIII)
オンライン計算処理部４０７は、規定データ設定部４０３で設定された「空気の密度ρg,3、空気の粘性率μg,3及び焼結鉱粒子の粒子径Ｄp,3」と、初期値設定部４０４で設定された「形状係数の初期値φs,3,0」と、操業データ収集部４０５で得られた「圧力損失ΔＰと焼結層の層厚Ｌ」とにより（１９）式からゾーンIIIにおける空気流速ｖg,3を計算する。
【００９１】
前述したように本実施形態では、乾燥帯であるゾーンIIIの空気流入温度θgi,3は、赤熱帯であるゾーンIVの空気流入温度θgi,4と等しいと仮定している。そこで、オンライン計算処理部４０７は、ゾーンIIIの空気流入温度θgi,3として、ゾーンIVの空気流入温度θgi,4を設定する。ゾーンIVの空気流入温度θgi,4を求め方については後述する。尚、空気流入温度θgi,3には、他の知見等に基づき、空気流入温度θgi,4以外の別の温度を設定しても構わない。
【００９２】
次に、オンライン計算処理部４０７は、ゾーン解析部４０２で得られた「ゾーンIIIの始点温度θp,3s及び終点温度θp,3e」と、初期値設定部４０４で設定された「空気流入温度θgi,3,0」と、ゾーン解析部４０２で得られた「ゾーンIIIのゾーン通過時間ｔ3」とを用いてゾーンIIIにおける時定数Ｔ3を求める。また、オンライン計算処理部４０７は、「焼結鉱粒子の比熱Ｃp,3」を規定データ設定部４０３から読み取る。
【００９３】
また、オンライン計算処理部４０７は、前述したようにして計算した「ゾーンIIIにおける空気流速ｖg,3」と、規定データ設定部４０３で設定された「空気の密度ρg,3、粘性率μg,3、熱伝導率λg,3及び定圧比熱Ｃpg,3、及び焼結鉱粒子の直径Ｄp,3」とにより（２２）式からＮｕ3を求める。そして、オンライン計算処理部４０７は、（２１）式から対流熱伝達率ｈ3を計算する。
【００９４】
そして、オンライン計算処理部４０７は、（２３）式〜（２６）式より、形状係数φs,3の中間計算値φs,3,tmpを求める。そして、オンライン計算処理部４０７は、形状係数φs,3の中間計算値φs,3,tmpと、それまでに設定されている形状係数φs,3（最初は初期値φs,3,0）との差の絶対値が収束判定基準値ＥＰＳ以下であるか否かを判定する。
この判定の結果、形状係数φs,3の中間計算値φs,3,tmpと、それまでに設定されている形状係数φs,3との差の絶対値が収束判定基準値ＥＰＳ以下でない場合、オンライン計算処理部４０７は、それまでに設定されている形状係数φs,3を、中間計算値φs,3,tmpに置き換える（最初は初期値φs,3,0を中間計算値φs,3,tmpに置き換える）。そして、オンライン計算処理部４０７は、この中間計算値φs,3,tmpに置き換えた形状係数φs,3を用いて、前述した「空気流速ｖg,3」、「形状係数φs,3の中間計算値φs,3,tmp」を再度求める。オンライン計算処理部４０７は、このような処理を、形状係数φs,3の中間計算値φs,3,tmpと、それまでに設定されている形状係数φs,3との差の絶対値が収束判定基準値ＥＰＳ以下になるまで繰り返し行う。
【００９５】
そして、形状係数φs,3の中間計算値φs,3,tmpと、それまでに設定されている形状係数φs,3との差の絶対値が収束判定基準値ＥＰＳ以下になると、オンライン計算処理部４０７は、そのときに得られた「空気流速ｖg,3」、「空気流入温度θgi,3」、「時定数Ｔ3」、「対流熱伝達率ｈ3」、「形状係数φs,3」を最終計算値として確定し、これらを同定パラメータとしてパラメータ蓄積部４０８に記憶する。
【００９６】
(ゾーンIV)
オンライン計算処理部４０７は、規定データ設定部４０３で設定された「空気の密度ρg,4、空気の粘性率μg,4及び焼結鉱粒子の粒子径Ｄp,4」と、初期値設定部４０４で設定された「形状係数の初期値φs,4,0」と、操業データ収集部４０５で得られた「圧力損失ΔＰと焼結層の層厚Ｌ」とにより（１９）式からゾーンIVにおける空気流速ｖg,4を計算する。
【００９７】
次に、オンライン計算処理部４０７は、初期値設定部４０４で設定された「初期燃料比Ｋc0」と、ゾーン解析部４０２で得られた「ゾーンIVのゾーン通過時間ｔ4」と、規定データ設定部４０３で設定された「燃料の密度ρc、及び燃料のモル重量Ｗc」と、（２３）式で計算される「ゾーンIVにおける焼結鉱粒子の体積Ｖp,4」とにより（２８）式に基づく式から燃料燃焼速度Ｒcを計算する。
【００９８】
次に、オンライン計算処理部４０７は、前述したようにして計算した「ゾーンIVにおける空気流速ｖg,4」と、規定データ設定部４０３で設定された「空気の密度ρg,4、粘性率μg,4、熱伝導率λg,4及び定圧比熱Ｃpg,4、及び焼結鉱粒子の直径Ｄp,4」とにより（２２）式からＮｕ4を求める。そして、オンライン計算処理部４０７は、（２１）式から対流熱伝達率ｈ4を計算する。
また、オンライン計算処理部４０７は、初期値設定部４０４で設定された「形状係数の初期値φs,4,0」と、規定データ設定部４０３で設定された「ゾーンIVにおける焼結鉱粒子の直径Ｄp,4」とにより（２４）式からゾーンIVにおける焼結鉱粒子の伝熱表面積Ａp,4を求める。
【００９９】
そして、オンライン計算処理部４０７は、ゾーン解析部４０２で得られた「ゾーンIVの終点温度θp,4e」と、規定データ設定部４０３で設定された「燃料の燃焼熱量Ｑc」と、前述したようにして計算された「燃料燃焼速度Ｒc」、「対流熱伝達率ｈ4」及び「ゾーンIVにおける焼結鉱粒子の伝熱表面積Ａp,4」とにより（２０）式、（２８）式に基づく式からゾーンIVの空気流入温度θgi,4の中間計算値θgi,4,tmpを求める。
【０１００】
次に、オンライン計算処理部４０７は、（２５）式より、ゾーンIVにおける時定数Ｔ4を計算する。そして、オンライン計算処理部４０７は、（２６）式より、形状係数φs,4の中間計算値φs,4,tmpを求める。そして、オンライン計算処理部４０７は、形状係数φs,4の中間計算値φs,4,tmpと、それまでに設定されている形状係数φs,4（最初は初期値φs,4,0）との差の絶対値が収束判定基準値ＥＰＳ以下であるか否かを判定する。
【０１０１】
この判定の結果、形状係数φs,4の中間計算値φs,4,tmpと、それまでに設定されている形状係数φs,4との差の絶対値が収束判定基準値ＥＰＳ以下でない場合、オンライン計算処理部４０７は、それまでに設定されている形状係数φs,4を、中間計算値φs,4,tmpに置き換える（最初は初期値φs,4,0を中間計算値φs,4,tmpに置き換える）。そして、オンライン計算処理部４０７は、この中間計算値φs,4,tmpに置き換えた形状係数φs,4を用いて、前述した「空気流速ｖg,4」、「燃料燃焼速度Ｒc」、「空気流入温度θgi,4の中間計算値θgi,4,tmp」、「形状係数φs,4の中間計算値φs,4,tmp」を再度求める。オンライン計算処理部４０７は、このような処理を、形状係数φs,4の中間計算値φs,4,tmpと、それまでに設定されている形状係数φs,4との差の絶対値が収束判定基準値ＥＰＳ以下になるまで繰り返し行う。
【０１０２】
そして、形状係数φs,4の中間計算値φs,4,tmpと、それまでに設定されている形状係数φs,4との差の絶対値が収束判定基準値ＥＰＳ以下になると、オンライン計算処理部４０７は、空気流入温度θgi,4の中間計算値θgi,4,tmpと、それまでに設定されている空気流入温度θgi,4（最初は初期値θgi,4,0）との差の絶対値が収束判定基準値ＥＰＳ以下であるか否かを判定する。この判定の結果、空気流入温度θgi,4の中間計算値θgi,4,tmpと、それまでに設定されている空気流入温度θgi,4との差の絶対値が収束判定基準値ＥＰＳ以下でない場合、オンライン計算処理部４０７は、それまでに設定されている空気流入温度θgi,4を、中間計算値θgi,4,tmpに置き換える（最初は初期値θgi,4,0を中間計算値θgi,4,tmpに置き換える）。そして、オンライン計算処理部４０７は、中間計算値θgi,4,tmpに置き換えた空気流入温度θgi,4を用いて、前述した「空気流速ｖg,4」、「燃料燃焼速度Ｒc」、「空気流入温度θgi,4の中間計算値θgi,4,tmp」、「形状係数φs,4の中間計算値φs,4,tmp」の計算を繰り返し行う。
【０１０３】
そして、ゾーンIVの空気流入温度θgi,4の中間計算値θgi,4,tmpと、それまでに設定されているゾーンIVの空気流入温度θgi,4（最初は初期値θgi,4,0）との差の絶対値が収束判定基準値ＥＰＳ以下になると、オンライン計算処理部４０７は、そのときに得られた「空気流速ｖg,4」、「燃料燃焼速度Ｒc」、「対流熱伝達率ｈ4」、「空気流入温度θgi,4」、「時定数Ｔ4」、「形状係数φs,4」を最終計算値として確定し、これらを同定パラメータとしてパラメータ蓄積部４０８に記憶する。
【０１０４】
（ゾーンV）
オンライン計算処理部４０７は、規定データ設定部４０３で設定された「空気の密度ρg,5、空気の粘性率μg,5及び焼結鉱粒子の粒子径Ｄp,5」と、初期値設定部４０４で設定された「形状係数の初期値φs,5,0」と、操業データ収集部４０５で得られた「圧力損失ΔＰと焼結層の層厚Ｌ」とにより（１９）式からゾーンVの空気流速ｖg,5を計算する。
【０１０５】
前述したように本実施形態では、焼結完了帯であるゾーンVの空気流入温度θgi,5は、水分凝縮帯であるゾーンIIの空気流入温度θgi,2と等しいと仮定している。そこで、オンライン計算処理部４０７は、空気流入温度θgi,5として、ゾーンIIの空気流入温度θgi,2を設定する。ゾーンIIの空気流入温度θgi,2の求め方は前述した通りである。尚、空気流入温度θgi,5には、他の知見等に基づき、空気流入温度θgi,2以外の別の温度を設定しても構わない。
【０１０６】
次に、オンライン計算処理部４０７は、ゾーン解析部４０２で得られた「ゾーンVの始点温度θp,5s及び終点温度θp,5eと、初期値設定部４０４で設定された「空気流入温度θgi,5,0」とを用いてゾーンVの時定数Ｔ5を求める。また、オンライン計算処理部４０７は、「焼結鉱粒子の比熱Ｃp,5」を規定データ設定部４０３から読み取る。
【０１０７】
また、オンライン計算処理部４０７は、前述したようにして計算した「ゾーンVの空気流速ｖg,5」と、規定データ設定部４０３で設定された「空気の密度ρg,5、粘性率μg,5、熱伝導率λg,5及び定圧比熱Ｃpg,5、及び焼結鉱粒子の直径Ｄp,5」とにより（２２）式からＮｕ5を求める。そして、オンライン計算処理部４０７は、（２１）式から対流熱伝達率ｈ5を計算する。
【０１０８】
そして、オンライン計算処理部４０７は、（２３）式〜（２６）式より、形状係数φs,5の中間計算値φs,5,tmpを求める。そして、オンライン計算処理部４０７は、形状係数φs,5の中間計算値φs,5,tmpと、それまでに設定されている形状係数φs,5（最初は初期値φs,5,0）との差の絶対値が収束判定基準値ＥＰＳ以下であるか否かを判定する。
この判定の結果、形状係数φs,5の中間計算値φs,5,tmpと、それまでに設定されている形状係数φs,5との差の絶対値が収束判定基準値ＥＰＳ以下でない場合、オンライン計算処理部４０７は、それまでに設定されている形状係数φs,5を、中間計算値φs,5,tmpに置き換える（最初は初期値φs,5,0を中間計算値φs,5,tmpに置き換える）。そして、オンライン計算処理部４０７は、この中間計算値φs,5,tmpに置き換えた形状係数φs,5を用いて、前述した「空気流速ｖg,5」、「形状係数φs,5の中間計算値φs,5,tmp」を再度求める。オンライン計算処理部４０７は、このような処理を、形状係数φs,5の中間計算値φs,5,tmpと、それまでに設定されている形状係数φs,5との差の絶対値が収束判定基準値ＥＰＳ以下になるまで繰り返し行う。
【０１０９】
そして、形状係数φs,5の中間計算値φs,5,tmpと、それまでに設定されている形状係数φs,5との差の絶対値が収束判定基準値ＥＰＳ以下になると、オンライン計算処理部４０７は、そのときに得られた「空気流速ｖg,5」、「空気流入温度θgi,5」、「時定数Ｔ5」、「対流熱伝達率ｈ5」、「形状係数φs,5」を最終計算値として確定し、これらを同定パラメータとしてパラメータ蓄積部４０８に記憶する。
【０１１０】
［［測定温度データに基づく関数の導出］］
次に、オンライン計算処理部４０７は、以上のようにして同定された「各ゾーンII〜IV(V)におけるパラメータ」をそれぞれ（２０）式で例示する時間微分方程式の解として得られる関数に適用して、各ゾーンII〜IV(V)における「『温度測定部１７０ａが配置されている位置での測定温度データに基づく関数θp,iA(ｔiA)』と『温度測定部１７０ｂが配置されている位置での測定温度データに基づく関数θp,iB(ｔiB)』」を求める。ここで、ｉは、既に説明したようにゾーンを識別するものであり、１〜５の整数となる。
【０１１１】
［［処理の流れ］］
以上の［［パラメータの同定］］と［［測定温度データに基づく関数の導出］］は、例えば、以下の順序で行われる。
（１）焼結機１００において、焼結原料を積載して搬送する無端帯状のパレット群の中から、パレット温度計１５０を装備したパレット１１０を１つ指定する。
（２）（１）で指定したパレット１１０に配置された複数のパレット温度計１５０ａ〜ｆの中から、パレット温度計１５０を１つ指定する。
（３）（２）で指定したパレット温度計１５０の中から、一方の温度測定部１７０を指定する。
（４）（３）で指定したパレット温度計１５０の温度測定部１７０（例えば温度測定部１７０ａ）で得られた測定温度データについての、ゾーン解析部４０２で解析された結果を使って、パラメータを同定し、測定温度データに基づく関数（例えばθp,iA(ｔiA)）を導出する。
（５）（３）で指定したパレット温度計１５０の他方の温度測定部１７０（例えば温度測定部１７０ｂ）で得られた測定温度データについての、ゾーン解析部４０２で解析された結果を使って、パラメータを同定し、測定温度データに基づく関数（例えばθp,iB(ｔiB)）を導出する。
尚、パラメータの同定に際し、導出する対象の関数から明らかに外れている測定温度データについては採用しないようにする。
【０１１２】
（６）（１）で指定したパレット１１０に配置された複数のパレット温度計１５０ａ〜ｆの中から、（２）で指定したパレット温度計１５０以外の、未指定のパレット温度計１５０を１つ指定する。
（７）（６）で指定したパレット温度計１５０に対して、（３）〜（５）の処理を行う。
（８）（１）で指定したパレット１１０に配置された複数のパレット温度計１５０ａ〜ｆの中から、未指定のパレット温度計１５０がなくなるまで（６）、（７）の処理を繰り返し行う。
本実施形態では、（１）で指定するパレット温度計１５０を装備したパレット１１０において、パレット温度計１５０の数が６個であるので、（３）〜（５）の処理が６回行われることになる。
また、焼結機１００において、焼結原料を積載して搬送する無端帯状のパレット群において、パレット温度計１５０を装備したパレット１１０が複数ある場合には、その数（回数）だけ（１）〜（８）の処理を繰り返し行う。
【０１１３】
例えば、ゾーンIV（赤熱帯）における「『温度測定部１７０ａが配置されている位置での測定温度データに基づく関数θp,iA(ｔiA)』と『温度測定部１７０ｂが配置されている位置での測定温度データに基づく関数θp,iB(ｔiB)』」は、それぞれ（−ｔiA／Ｔ4A）および（−ｔiB／Ｔ4B）を変数とする指数関数で表される。
【０１１４】
当該関数それぞれにおいて、添字Ａは、温度測定部１７０ａが配置されている位置のものであることを示し、添字Ｂは、温度測定部１７０ｂが配置されている位置のものであることを示す。また、時間ｔ4A、ｔ4Bは、それぞれ、パレット温度計１５０の温度測定部１７０ａ、１７０ｂがゾーンIVに入ったときの時間を０（ゼロ）としたときの時間である。このことは、他のゾーンII〜IV(V）の温度θp,iA(ｔiA)、θp,iB(ｔiB)についても同じである。
【０１１５】
［［導出対象温度データを定めるためのパラメータの同定］］
オンライン計算処理部４０７は、各ゾーンII〜IV(V)における「温度測定部１７０ａ、１７０ｂが配置されている位置での測定温度データに基づく関数θp,iA(ｔiA)、θp,iB(ｔiB)」のうち、同一のパレット温度計１５０から得られたものを選択する。
次に、オンライン計算処理部４０７は、温度測定部１７０ａ、１７０ｂの位置の情報を、規定データ設定部４０３から読み出す。例えば、オンライン計算処理部４０７は、基台部１６１の上面から温度測定部１７０ｂ（の中心）までのｚ軸方向における距離（前述したように本実施形態では３０［ｍｍ］）と、温度測定部１７０ｂ（の中心）から温度測定部１７０ａ（の中心）までのｚ軸方向における距離（本実施形態では１５０［ｍｍ］）とを、温度測定部１７０ａ、１７０ｂの位置の情報の一例として、規定データ設定部４０３から読み出す。
【０１１６】
次に、オンライン計算処理部４０７は、「温度測定部１７０ａが配置されている位置での測定温度データに基づく関数θp,iA(ｔiA)上の点」と、当該点とゾーン内での位置が相互に対応する「温度測定部１７０ｂが配置されている位置での測定温度データに基づく関数θp,iB(ｔiB)上の点」との組を複数抽出する。前述したように本実施形態では、各ゾーンII〜IV(V)における時間ｔの最大値と最小値を含み、且つ、それぞれの時間隔が略等間隔になる２０個の時間における温度を、それぞれの関数θp,iA(ｔiA)、θp,iB(ｔiB)から抽出する。オンライン計算処理部４０７は、この処理をゾーンII〜IV(V)に個別に行う。
【０１１７】
次に、オンライン計算処理部４０７は、例えば、抽出した点の組の値に対して外挿又は内挿を行って、ｚ軸方向の所定の位置における導出対象温度データを求める。オンライン計算処理部４０７は、この処理をゾーンII〜IV(V)毎に個別に行う。
次に、オンライン計算処理部４０７は、導出対象温度データを使って、前記ｚ軸方向の指定された位置における導出対象温度データに基づく関数θp,iC(ｔiC)を定めるためのパラメータを同定する。このパラメータの同定は、例えば、前述した［［パラメータの同定］］で説明したのと同じ方法で行うことができる。オンライン計算処理部４０７は、この処理をｚ軸方向に任意の所定の位置における導出対象温度データに対して実施し、ｚ軸方向に任意の所定の位置における関数θp,iC(ｔiC)を定めるためのパラメータを同定することができる。そして、オンライン計算処理部４０７は、この処理をゾーンII〜IV(V)毎に個別に行う。尚、パラメータの同定に際し、導出する対象の関数から明らかに外れている導出対象温度データについては採用しないようにする。
【０１１８】
［［導出対象温度データに基づく関数の導出］］
次に、オンライン計算処理部４０７は、以上のようにして同定された「各ゾーンII〜IV(V)におけるパラメータ」をそれぞれ（２０）式で例示する時間微分方程式の解として得られる関数に適用して、「ｚ軸方向の指定された位置における導出対象温度データに基づく関数θp,iC(ｔiC)」を、各ゾーンII〜IV(V)においてそれぞれ求める。
【０１１９】
例えば、ゾーンIVにおける「ｚ軸方向の指定された位置における導出対象温度データに基づく関数θp,4C(ｔ4C)」は、（−ｔ4C／Ｔ4C）を変数とする指数関数で表される。
【０１２０】
オンライン計算処理部４０７は、以上の処理を、全てのパレット温度計１５０ａ〜１５０ｆについて実行する。
【０１２１】
［［焼結層の高さ方向における所望の位置での温度の導出方法の一例］］
図８は、焼結層のゾーンIV（赤熱帯）の領域と、測定温度データと、導出対象温度データの一例を示す図である。
図８に示す例では、温度測定部１７０ａにより測定されたゾーンIVの測定温度データに基づく関数９０１と、温度測定部１７０ｂにより測定されたゾーンIVの測定温度データに基づく関数９０２とが得られたものとする。また、図８に示す例では、基台部１６１の上面（焼結層の床部）から温度測定部１７０ｂ（の中心）までのｚ軸方向における距離をＬBCと表記し、温度測定部１７０ｂ（の中心）から温度測定部１７０ａ（の中心）までのｚ軸方向における距離をＬABと表記している。
【０１２２】
また、図８に示す例では、温度測定部１７０ａがゾーンIVに入った地点をＡ１で表記し、温度測定部１７０ａがゾーンIVから出た地点をＡ２で表記している。同様に、温度測定部１７０ｂがゾーンIVに入った地点をＢ１で表記し、温度測定部１７０ｂがゾーンIVから出た地点をＢ２で表記している。
【０１２３】
図８の点Ａ１と点Ａ２とを結ぶ直線と、点Ｂ１と点Ｂ２とを結ぶ直線と、で囲まれる領域（図８の斜線で示す領域）がゾーンIVに対応する。尚、図８の横軸は時間であるので、当該時間と、平均パレット速度ｖpとを掛け合わせることにより、当該時間を、原料給鉱部からの距離に換算することができる。また、全てのパレット温度計１５０ａ〜１５０ｆについて図８の斜線で示す領域を求め、それらの領域の関数を基にパレット幅方向に内挿又は外挿を行うことで、ゾーンIVの３次元空間を得ることができる。このことは、その他のゾーンについても同じである。
【０１２４】
図９は、測定温度データに基づく関数９０１、９０２から導出対象温度データに基づく関数９０３を求める様子の一例を概念的に示す図である。
図９において、例えば、測定温度データに基づく関数９０１上の点９０１ｅは、温度測定部１７０ａがゾーンIVから出た地点における測定温度データによるものである。また、測定温度データに基づく関数９０２上の点９０２ｅは、温度測定部１７０ｂがゾーンIVから出た地点における測定温度データによるものである。よって、これらの点９０１ｅ、９０２ｅは、ゾーンIVにおける位置が相互に対応する点となる。
【０１２５】
また、本実施形態では、基台部１６１の上面（焼結層の床部）から温度測定部１７０ｂ（の中心）までのｚ軸方向における距離をＬBCは３０［ｍｍ］であり、温度測定部１７０ｂ（の中心）から温度測定部１７０ａ（の中心）までのｚ軸方向における距離をＬABは１５０［ｍｍ］である。よって、以下の（２９）式が成り立つ。
ＬAB：ＬBC＝５：１  ・・・（２９）
また、図５に示したように、各ゾーンＩ〜IVは、下部であるほど（焼結層の底部に近い位置であるほど）機長方向の長さが長くなることが知られている。
【０１２６】
以上のことから、点９０１ｅ、９０２ｅを結ぶ直線を、点９０２ｅの方向に外挿し、点９０１ｅ、９０２ｅを結ぶ直線の長さの１／５の長さになる点９０３ｅを導出対象温度データ９０３ｅとする。このような処理を、例えば、時間ｔ4A、ｔ4Bの最大値と最小値を含み、且つ、それぞれの時間隔が略等間隔になる２０個の時間について行い、２０個の導出対象温度データを求め、それらから、導出対象温度データに基づく関数９０３を求める。尚、点９０１ｅ、９０２ｅの間の高さの点を導出対象温度データとする際には、点９０１ｅ、９０２ｅに基づく内挿を行うようにする。さらに、求めたい導出対象温度データ点９０３ｅの高さが点９０１ｅよりも上にある場合には、点９０１ｅ、９０２ｅを結ぶ直線を、点９０１ｅの方向に外挿を行うようにする。
【０１２７】
図８に示す例では、基台部１６１の上面（焼結層の床部、すなわちｚ＝０）において、IVゾーンの導出対象温度データに基づく関数９０３を求めた例を示す。基台部１６１の上面（焼結層の床部）におけるゾーンIVの開始地点をＣ１で表記し、終了地点をＣ２で表記する。この基台部１６１の上面（焼結層の床部）におけるゾーンIVの終了地点Ｃ２がＢＴＰに対応する。
尚、曲線９０３を求めるに際し、図９に示す例では、説明のため、指定した一つのパレット温度計１５０の測定温度データ９０１、９０２及び導出対象温度データ９０３について、これらのゾーンIVの開始地点の温度θp,4sA、θp,4sB、θp,4sCの点が一致する例を示しているが、サンプリング誤差や測定ノイズによる誤差等の外乱により、これらの点が一致しない場合には、例えば、これらの点が一致するように、θp,4sA、θp,4sB、θp,4sC値を変更したりすることができる。
オンライン計算処理部４０７は、以上の処理を、全てのパレット温度計１５０ａ〜１５０ｆについて個別に行う。
【０１２８】
［［焼結層のパレット幅方向における所望の位置での温度の導出］］
前述したように、本実施形態では、指定した一つのパレット温度計１５０において、焼結層の高さ方向における所望の位置での温度を、焼結層内の物理現象を定量的に表現するパラメータを表現するモデルに基づいて求めるようにした。本実施形態では、これに加えて、焼結層のパレット幅方向における所望の位置での温度を、焼結層の同一の高さ位置における温度データ（測定温度データ、導出対象温度データ）に基づき、空間的に補間して求めるようにする。
【０１２９】
具体的に説明すると、本実施形態では、オンライン計算処理部４０７は、まず初めに、各ゾーンII〜IV(V)における指定する機長方向のｙ座標位置について、同一の焼結層高さ位置での温度データとして、焼結層高さ位置が相互に対応する温度データ（測定温度データ、導出対象温度データ）を、前記同定したパラメータに基づく関数θp,ijA(ｔiA)、θp,ijB(ｔiB)、θp,ijC(ｔiC)と前記指定した機長方向ｙ座標位置に対応する時間ｔiA、ｔiB、ｔiCの情報から抽出又は生成する。
ここで、ｉは前述したようにゾーンを識別するものであり（ｉ＝１〜５はそれぞれゾーンI〜Vに対応する）、ｊはパレット幅方向（ｘ軸方向）に配置されているパレット温度計１５０を識別する１〜ｍまでの整数値である。本実施形態ではｍ＝６であり、ｊ＝１〜６がそれぞれパレット温度計１５０ａ〜ｆに対応する。
次に、オンライン計算処理部４０７は、前記抽出したｍ個の温度データから、各ゾーンII〜IV(V)における指定する機長方向のｙ座標位置について、同一の焼結層高さ位置における焼結層のパレット幅方向（ｘ軸方向）の温度分布を表現する（ｍ−１）次の多項式関数を導出する。このとき、（ｍ−１）次の多項式関数の係数はｍ個の温度データから代数的に求めることができるし、必要に応じて最小二乗法等を組み合わせて用いてもよい。
【０１３０】
図１０は、あるゾーンにおけるある機長方向位置での同一の焼結層高さ位置のパレット幅方向の温度の変化を表す関数の一例を示す図である。
図１０において、温度データ１１０１〜１１０７は、同一の焼結層高さ位置の温度データであり、曲線１１００が、それらから導出された多項式関数を示している。
尚、本実施形態では、サイドウォール１１４に直近の位置にパレット温度計１５０ａ、１５０ｆを配置しているので、サイドウォール１１４に直近の領域における焼結層の温度も直接的に且つ正確に測定することができる。
そして、本実施形態では、前述した各ゾーンII〜IV(V)において指定する機長方向位置について、同一の焼結層高さ位置での温度データの抽出又は生成にあたっては、各ゾーンII〜IV(V)における時間ｔの最大値と最小値を含み、且つ、それぞれの時間隔が略等間隔になる時間における温度データを抽出又は生成するようにして、それぞれ同一の焼結層高さ位置における焼結層のパレット幅方向（ｘ軸方向）の温度分布を表現する（ｍ−１）次の多項式関数を導出することにより、これら導出した（ｍ−１）次の多項式関数を用いて、最終的に、焼結層のパレット幅方向（ｘ座標）・機長方向（ｙ座標）・高さ方向（ｚ座標）における所望の位置での焼結層の温度を得る。
【０１３１】
［［ＢＴＰの導出］］
以上のようにして、焼結層のパレット幅方向（ｘ座標）・機長方向（ｙ座標）・高さ方向（ｚ座標）における所望の位置での焼結層の温度が得られるので、オンライン計算処理部４０７は、焼結層内の３次元温度分布を得ることができる。この他に、本実施形態では、オンライン計算処理部４０７は、焼結機の操業指標として重要なＢＴＰを導出するようにしている。
そのために、本実施形態では、オンライン計算処理部４０７は、一つのパレット温度計１５０で得られた基台部１６１の上面（焼結層の床部、すなわちｚ＝０）におけるIIゾーン〜IVゾーンの導出対象温度データに基づく関数（θp,2C(ｔ2C)、θp,3C(ｔ3C)、θp,4C(ｔ4C)）をゾーンの境界で繋ぐ。次に、オンライン計算処理部４０７は、当該繋いだ関数から、ゾーンIVのゾーン通過時間ｔ4C及びゾーン境界時間ｔa,4Cを求め（図６や（１１）式を参照）、ゾーンIVのゾーン通過時間ｔ4C又はゾーン境界時間ｔa,4Cと平均パレット速度ｖpとを掛け合わせることによりＢＴＰを導出する。このとき、平均パレット速度ｖpは、パレット温度計１５０が焼結鉱排鉱部に到達後に（１５）式及び（１６）式で算出して用いてもよいし、パレット温度計１５０がIVゾーンから出た時点で（１５´）式及び（１６´）式で算出して用いてもよい。
オンライン計算処理部４０７は、このようなＢＴＰの導出を、パレット幅方向（ｘ軸方向）に配置されているｍ個のパレット温度計１５０毎に実施し、それぞれで得られたパレット幅方向（ｘ軸方向）のｍ個のＢＴＰを用いて、パレット幅方向（ｘ軸方向）のＢＴＰ分布を表現する（ｍ−１）次の多項式関数を導出する。これにより、本実施形態では、ＢＴＰのパレット幅方向（ｘ軸方向）の分布を導出することができる。
【０１３２】
本実施形態では、オンライン計算処理部４０７は、パレット温度計１５０が原料給鉱部から焼結鉱排鉱部まで移動する度に（１スイープする度に）、「焼結層の各ゾーンI〜Vの領域」と、「焼結層の任意の位置における温度（焼結層全体の３次元温度分布）」と、「ＢＴＰのパレット幅方向（ｘ軸方向）分布」とを繰り返し導出する。
【０１３３】
［オンライン操業監視表示部４０９］
オンライン操業監視表示部４０９は、オンライン計算処理部４０７で得られた「『焼結層の各ゾーンI〜Vの領域』と、『焼結層の任意の位置における温度（焼結層全体の３次元温度分布）』と、『ＢＴＰのパレット幅方向（ｘ軸方向）分布』とに基づいて、それらを表示する表示データを生成して表示装置４７０に表示する。
【０１３４】
図１１は、焼結層のゾーンIVの概略の一例を３次元的に表現して示した図である。図１１において、パレット幅方向（ｘ軸方向）の側面と底面とを斜線で示した領域がゾーンIV（赤熱帯）となる。また、図１１では図示を省略するが、この領域内における温度の分布を、例えば、温度が同値である任意の等値線の表示や、所定の温度範囲毎の色分け表示等を行うことによって、表示することができる。尚、温度の分布を表示する方法は、公知の技術で実現できるので、その詳細な説明を省略する。等値線については、例えば、特許第３８１４１４３号公報に記載の技術を用いて探索することができる。
【０１３５】
また、図１１において、ＢＴＰの表記の横に付されている矢印が指している破線１２０１がＢＴＰのパレット幅方向（ｘ軸方向）分布となる。これにより、オペレータは、パレット幅方向（ｘ軸方向）にＢＴＰがどのように分布しているのかを把握することができる。尚、図１１では、ＢＴＰのパレット幅方向（ｘ軸方向）分布を表現する破線１２０１を直線で示しているが、オンライン計算処理部４０７におけるＢＴＰのパレット幅方向（ｘ軸方向）分布の導出結果によっては、当該破線１２０１は曲線になる。このことは、後述の[変形例]で述べるように、ゾーンIVの側面を構成する直線１２０２ａ〜１２０２ｅについても同じである。
【０１３６】
本実施形態では、オンライン計算処理部４０７は、パレット温度計１５０が原料給鉱部から焼結鉱排鉱部まで移動する度に（１スイープする度に）、「焼結層の各ゾーンI〜Vの領域」と、「焼結層の任意の位置における温度（焼結層全体の３次元温度分布）」と、「ＢＴＰのパレット幅方向（ｘ軸方向）分布」とを得る。よって、オンライン操業監視表示部４０９は、これらの情報が得られる度に（パレット温度計１５０が１スイープする度に）、表示の更新を行うことができる。また、このとき、パレット温度計１５０を装備するパレット１１０は、無端帯状に連結されて周回する複数のパレットの中で１つに限定される必要はなく、複数設置し、各パレットのスイープを追いながら、前記表示の更新を行うことができる。
【０１３７】
尚、ここでは、ゾーンIVについてのみ、領域と温度の分布との表示を行う場合を例に挙げて説明した。しかしながら、その他のゾーンについても、ゾーンIVと同様に、領域と温度の分布とを表示するようにしてもよい。また、ゾーンI〜Vのうちの全部又は一部の複数のゾーンについて、それらのゾーンの領域と温度の分布とを組み合わせて同時に（例えば同じ画面上に）表示するようにしてもよい。
【０１３８】
［動作フローチャート］
次に、図１２のフローチャートを参照しながら、焼結プロセス操業監視装置４００の処理動作の一例を説明する。
まず、ステップＳ１において、事前設定が行われる。具体的に説明すると、規定データ設定部４０３は、ユーザによる入力装置４６０の操作に基づいて、各種の定数と必要に応じてその温度依存特性関数とを設定（ＨＤＤ等に記憶）する。
次に、ステップＳ２において、初期値設定が行われる。具体的に説明すると、初期値設定部４０４は、ユーザによる入力装置４６０の操作に基づいて、各ゾーンI〜Vについて、空気流入温度θgi,0［℃］、形状係数φs,i,0［−］を設定（ＨＤＤ等に記憶）する。
【０１３９】
次に、温度データ収集部４０１は、パレット温度計１５０が、原料給鉱部から焼結鉱排鉱部まで移動したときに得られたパレット温度計１５０の測定温度データθp（ｔa）（１スイープ分のデータ）を取得するまで待機する。パレット温度計１５０の測定温度データθp（ｔa）が得られると（ステップＳ３）、ステップＳ４に進む。ステップＳ４に進むと、操業データ収集部４０５は、プロセスコンピュータ４５０から入力されてＨＤＤ等に記憶しておいた操業データとして、圧力損失（ブロア吸引圧）ΔＰ［Ｐａ］と、焼結層の層厚Ｌ［ｍ］と、初期水分比Ｋw0［−］と、初期燃料比Ｋc0［−］とを読み出す。
【０１４０】
次に、ステップＳ５において、ゾーン解析部４０２は、パレット温度計１５０の測定温度データθp（ｔa）に基づいて、各ゾーンの境界の温度（ゾーン境界温度）θp（０）、θp（ｔa,min）、θp（ｔa,cf）、θp（ｔa,max）θp（ｔ´a,max）と、ゾーン境界時間ｔa,1（＝０）、ｔa,2（＝ｔa,min）、ｔa,3（＝ｔa,cf）、ｔa,4（＝ｔa,max）、ｔa,5（＝ｔ´a,min）とを読み出す。
次に、ステップＳ６において、ゾーン解析部４０２は、パレット温度計１５０の測定温度データθp（ｔa）に基づいて、各ゾーンの始点温度θp,is（（１）式、（３）式、（６）式、（９）式、（１２）式）・終点温度θp,ie（（１）式、（４）式、（７）式、（１０）式、（１３）式）と、ゾーン通過時間ｔi（（２）式、（５）式、（８）式、（１１）式、（１４）式）とを設定する。
【０１４１】
次に、ステップＳ７において、ゾーン解析部４０２は、パレット温度計１５０の測定温度データθp（ｔa）に基づいて、（１５）式の計算を行って、平均パレット速度ｖpを設定する。また、ゾーン解析部４０２は、ゾーン通過時間ｔiと平均パレット速度ｖpとに基づいて、（１７）式の計算を行って、各ゾーンのゾーン長ｌiを設定する。また、ゾーン解析部４０２は、各ゾーンのゾーン長ｌiに基づいて（１８）式の計算を行って、各ゾーンのゾーン位置ｌa,iを設定する。
【０１４２】
次に、ステップＳ８において、オンライン計算処理部４０７によって、上述したようにII〜IVゾーンのオンライン計算処理が行われ、「各ゾーンI〜IVの測定温度データに基づく関数（θp,2A(ｔ2A)、θp,3A(ｔ3A)、θp,4A(ｔ4A)、θp,2B(ｔ2B)、θp,3B(ｔ3B)、θp,4B(ｔ4B)）と、導出対象温度データに基づく関数（θp,2C(ｔ2C)、θp,3C(ｔ3C)、θp,4C(ｔ4C)）」を得る。
【０１４３】
次に、ステップＳ９において、オンライン計算処理部４０７は、ステップＳ８で得られた「各ゾーンI〜IVの測定温度データに基づく関数（θp,2A(ｔ2A)、θp,3A(ｔ3A)、θp,4A(ｔ4A)、θp,2B(ｔ2B)、θp,3B(ｔ3B)、θp,4B(ｔ4B)）と、導出対象温度データに基づく関数（θp,2C(ｔ2C)、θp,3C(ｔ3C)、θp,4C(ｔ4C)）」に基づいて、各ゾーンＩ〜IVにおける焼結層内の３次元温度分布と、ＢＴＰのパレット幅方向（ｘ軸方向）分布と、を含む操業監視情報を生成する。
【０１４４】
次に、ステップＳ１０において、オンライン操業監視表示部４０９は、操業監視情報に基づく表示データを生成して表示装置４７０に表示する（図１１を参照）。
最後に、ステップＳ１１において、温度データ収集部４０１は、焼結機１００の操業監視を終了するか否かを判定する。この判定の結果、焼結機１００の操業監視を終了しない場合には、ステップＳ３に戻り、次のスイープ分のデータに基づく焼結機１００の操業監視を行う。一方、焼結機１００の操業監視を終了する場合には、図１２のフローチャートによる処理を終了する。
【０１４５】
［まとめ］
以上のように本実施形態では、パレット温度計１５０が、原料給鉱部から焼結鉱排鉱部まで移動したときに得られたパレット温度計１５０の温度測定部１７０ａ、１７０ｂで測定された測定温度データをオンラインで取得する。そして、測定温度データから、焼結層の温度を定量的に表現するパラメータを導出し、当該パラメータから、測定温度データに基づく関数θp,iA(ｔiA)、θp,iB(ｔiB)を求める。そして、求めた２つの関数上の、ゾーン内での位置が相互に対応する点に基づいて外挿又は内挿を行って、所望の高さ方向（ｚ軸方向）の位置における導出対象温度データを求める。そして、導出対象温度データから、焼結層の温度を定量的に表現するパラメータを導出し、当該パラメータから、導出対象温度データに基づく関数θp,iC(ｔiC)を求める。このような導出対象温度データに基づく関数θp,iC(ｔiC)の導出を、パレット１１０のパレット幅方向（ｘ軸方向）に並べられた全てのパレット温度計１５０ａ〜１５０ｆについてゾーン毎に行う。したがって、焼結層内のパレット幅方向（ｘ座標）・機長方向（ｙ座標）・高さ方向（ｚ座標）における所望の位置での焼結層の温度、すなわち焼結層内の３次元温度分布をオンラインで正確に導出することができ、導出した焼結層内の３次元温度分布に基づき、温度の等値線や所定の温度範囲毎の色分け表示等を行うことによって、焼結機の操業をオンラインで正確に監視することができる。
【０１４６】
また、本実施形態では、ゾーンIVの終点の位置におけるパレット幅方向（ｘ軸方向）の焼結層の温度の分布を求めることにより、ＢＴＰのパレット幅方向の分布をオンラインで正確に導出することができ、導出したＢＴＰのパレット幅方向分布に基づき、ＢＴＰの表示等を行うことによって、焼結機の操業をオンラインで正確に監視することができる。
【０１４７】
［変形例］
本実施形態では、１つのパレット温度計１５０に、２つの温度測定部１７０ａ、１７０ｂを設けた場合を例に挙げて説明した。前述したように１つのパレット温度計１５０に設ける温度測定部１７０の数は２つに限定されない。１つのパレット温度計１５０に、高さ方向（ｚ軸方向）にｎ（ｎは３以上の整数）個の温度測定部１７０を設ける場合には、次のようにして導出対象温度データを求めることができる。例えば、ｎ個の温度測定部１７０で測定された測定温度データに基づくｎ個の曲線上の相互に対応する点に基づいて外挿又は内挿を行って、所望の高さ方向の位置における導出対象温度データを求める。すなわち、ｎ個の温度測定部１７０で測定された測定温度データに基づくｎ個の曲線上の相互に対応する点から（ｎ−１）次多項式関数で表される曲線を求める。そして、求めた曲線上の点であって、ｎ個の温度測定部１７０の高さ方向の位置と、所望の高さ方向の位置とにより定まる点を導出対象温度データとして求める。この場合、例えば、図８に示す斜線の領域を区画する直線（点Ａ１、Ｂ１、Ｃ１を結ぶ曲線、点Ａ２、Ｂ２、Ｃ２を結ぶ直線）は曲線で表される。
【０１４８】
また、本実施形態では、各ゾーンを一つの区分として取り扱い、前記焼結層内の物理現象を定量的に表現するパラメータである焼結状態の特徴量を同定する一例を示したが、必要に応じてゾーンI〜Vのそれぞれを機長方向で２つ以上の区分に分割し、各区分について、各区分の始点及び終点の温度、物性値及び（１９）式〜（２８）式を用い、前記と同様の手順によって、各ゾーン内の各区分毎に前記特徴量を同定しても構わない。
【０１４９】
また、本実施形態を実現するためには、焼結層内の機長方向における温度分布の実測データを取得することは必須であるが、当該温度分布を測定する装置は、必ずしもパレット温度計１５０である必要はない。パレット温度計１５０以外の温度計で当該温度分布を測定することができれば、当該温度計を用いるようにしてもよい。
【０１５０】
尚、以上説明した本発明の実施形態は、コンピュータがプログラムを実行することによって実現することができる。また、プログラムをコンピュータに供給するための手段、例えばかかるプログラムを記録したＣＤ−ＲＯＭ等のコンピュータ読み取り可能な記録媒体、又はかかるプログラムを伝送する伝送媒体も本発明の実施の形態として適用することができる。また、前記プログラムを記録したコンピュータ読み取り可能な記録媒体などのプログラムプロダクトも本発明の実施の形態として適用することができる。前記のプログラム、コンピュータ読み取り可能な記録媒体、伝送媒体及びプログラムプロダクトは、本発明の範疇に含まれる。
また、以上説明した本発明の実施形態は、何れも本発明を実施するにあたっての具体化の例を示したものに過ぎず、これらによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその技術思想、またはその主要な特徴から逸脱することなく、様々な形で実施することができる。
【０１５１】
［本発明と実施形態との対応関係］
＜請求項１、４、５＞
温度測定部は、例えば、温度測定部１７０ａ、１７０ｂを用いることにより実現される。
測定温度データ収集手段は、例えば、温度データ収集部４０１を用いることにより実現される。
ゾーン解析手段は、例えば、ゾーン解析部４０２を用いることにより実現される。
測定温度関数のパラメータ同定手段は、例えば、オンライン計算処理部４０７が、各ゾーンII〜IVにおける「『温度測定部１７０ａが配置されている位置での測定温度データに基づく関数θp,iA(ｔiA)』と『温度測定部１７０ｂが配置されている位置での測定温度データに基づく関数θp,iB(ｔiB)』」を定めるパラメータを同定することにより実現される。また、このパラメータが、測定温度関数のパラメータ同定手段により同定されるパラメータの一例である。
導出対象温度データ導出手段は、例えば、オンライン計算処理部４０７が、温度測定部１７０ａが配置されている位置での測定温度データに基づく関数θp,iA(ｔiA)と、温度測定部１７０ｂが配置されている位置での測定温度データに基づく関数θp,iB(ｔiB)とから、ゾーンにおける位置が相互に対応する点（例えば図９の点９０１ａ、９０２ａ）を抽出し、抽出した点に基づく外挿又は内挿を行って、オペレータ等によって指定された高さ位置における点（例えば図９の点９０３ａ）を導出することにより実現される。
導出対象温度関数のパラメータ同定手段は、例えば、オンライン計算処理部４０７が、各ゾーンII〜IVにおける「導出対象温度データに基づく関数θp,iC(ｔiC)を定めるパラメータを同定することにより実現される。また、このパラメータが、導出対象温度関数のパラメータ同定手段により同定されるパラメータの一例である。
表示手段は、例えば、オンライン操業監視表示部４０９を用いることにより実現される。
焼結層内の高さ方向における温度分布は、例えば、温度測定部１７０ａが配置されている位置での測定温度データに基づく関数θp,iA(ｔiA)と、温度測定部１７０ｂが配置されている位置での測定温度データに基づく関数θp,iB(ｔiB)と、導出対象温度データに基づく関数θp,iC(ｔiC)とから、同値の温度の等値線を探索することにより実現される。
＜請求項２＞
パレット幅方向関数導出手段は、例えば、オンライン計算処理部４０７が、同一の高さ位置の温度データ（測定温度データ、導出対象温度データ）から得られた多項式関数（例えば図１０の曲線１１００等）を導出することにより実現される。また、この多項式関数が、焼結層内のパレット幅方向における温度の変化を表す関数の一例である。
焼結層内のパレット幅方向における温度分布は、例えば、多項式関数（例えば図１０の曲線１１００等）で表現することが可能であるし、更に、例えば、焼結層内の同一高さにおいて、機長方向に複数地点のパレット幅方向の温度分布を示す前記多項式関数で構成される温度分布から、同値の温度の等値線を探索することにより実現される。
＜請求項３＞
赤熱帯における前記焼結層の床部における機長方向の温度の変化を定量的に表現する関数は、例えば、図８等に示す導出対象温度データに基づく関数９０３により実現される。
赤熱帯の終点（ＢＴＰ）に対応する位置の導出対象温度データは、例えば、導出対象温度データ９０３ａにより実現される。そして、パレット幅方向の異なる位置で導出対象温度データ９０３ａを求めることにより（例えば、図１０に示す温度データ１１０１〜１１０７を求めることにより）、パレット幅方向で異なる位置の導出対象温度データの一例が得られる。したがって、赤熱帯の終点（ＢＴＰ）に対応する位置で、図１０に示すような温度データ１１０１〜１１０７を求めることにより、パレット幅方向で異なる位置の導出対象温度データであって、赤熱帯の終点（ＢＴＰ）に対応するパレット幅方向の異なる位置の導出対象温度データの一例が得られる。
【符号の説明】
【０１５２】
１００    焼結機
１０２    駆動ローラ
１１０    パレット
１２０    原料供給ホッパ
１３０    点火炉
１４０    ウインドボックス
１４５    吸気管
１５０ａ〜１５０ｆ    パレット温度計
１７０ａ、１７０ｂ    温度測定部
４００    焼結プロセス操業監視装置
４０１    温度データ収集部
４０２    ゾーン解析部
４０３    規定データ設定部
４０４    初期値設定部
４０５    操業データ収集部
４０６    モデル格納部
４０７    オンライン計算処理部
４０８    パラメータ蓄積部
４０９    オンライン操業監視表示部
４５０    プロセスコンピュータ
４６０    入力装置
４７０    表示装置
【図１】
0005729251.tif 000009
【図２】
0005729251.tif 000010
【図３】
0005729251.tif 000011
【図４】
0005729251.tif 000012
【図５】
0005729251.tif 000013
【図６】
0005729251.tif 000014
【図７】
0005729251.tif 000015
【図８】
0005729251.tif 000016
【図９】
0005729251.tif 000017
【図１０】
0005729251.tif 000018
【図１１】
0005729251.tif 000019
【図１２】
0005729251.tif 000020