JP3033466B2 - 高炉の操業方法 - Google Patents
高炉の操業方法Info
- Publication number
- JP3033466B2 JP3033466B2 JP7106498A JP10649895A JP3033466B2 JP 3033466 B2 JP3033466 B2 JP 3033466B2 JP 7106498 A JP7106498 A JP 7106498A JP 10649895 A JP10649895 A JP 10649895A JP 3033466 B2 JP3033466 B2 JP 3033466B2
- Authority
- JP
- Japan
- Prior art keywords
- furnace
- distribution
- model
- state
- blast furnace
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
Landscapes
- Manufacture Of Iron (AREA)
Description
【０００１】
【産業上の利用分野】本発明は、高炉の円周方向の分布
を含む３次元方向の非定常状態をモデル化した３次元非
定常モデルと炉頂からの装入物の分布モデルとを用いて
非定常状態における高炉の操業をシミュレートし、シミ
ュレート結果が実炉における計測値と一致するようにモ
デルの一部のプロセス定数を修正した後、修正したモデ
ルを用いて、溶銑温度、溶銑中Ｓｉ濃度、炉内の圧力分
布、燃料比の時間的変化等、操業結果の良否に影響を及
ぼす因子の非定常状態が目標の操業結果が得られる非定
常状態となるような送風条件及び装入物分布を推定して
実炉の非定常状態における操業条件とする高炉の操業方
法に関する。
を含む３次元方向の非定常状態をモデル化した３次元非
定常モデルと炉頂からの装入物の分布モデルとを用いて
非定常状態における高炉の操業をシミュレートし、シミ
ュレート結果が実炉における計測値と一致するようにモ
デルの一部のプロセス定数を修正した後、修正したモデ
ルを用いて、溶銑温度、溶銑中Ｓｉ濃度、炉内の圧力分
布、燃料比の時間的変化等、操業結果の良否に影響を及
ぼす因子の非定常状態が目標の操業結果が得られる非定
常状態となるような送風条件及び装入物分布を推定して
実炉の非定常状態における操業条件とする高炉の操業方
法に関する。
【０００２】
【従来の技術】高炉は、炉内において、気固液３層が共
存し、互いに反応、相変化、及び熱交換を行いながら移
動する複雑な反応容器である。このような高炉の炉内現
象を、基礎実験、解体調査等の知見を基に、移動速度論
的に整理し、組み合わせることにより、高炉を近似的に
数学モデル化することが可能である。数学モデルとして
は１次元定常モデル（特公昭６１−６１２２号、特公昭
６１−６０１２４号、特公昭６３−２４０４４号、特公
平１−１２８０５号公報）、２次元定常モデル（特公平
６−２８８６号公報）等が開発されてきた。これらのモ
デルでは、高炉の３次元的な挙動を高さ方向のみ取り扱
い（１次元モデル）、又は円周方向の分布を均一とし、
さらに時間的にも状態の変化がないとした定常状態のみ
の解析に限定し（２次元モデル）、モデルを簡略化する
ことで炉内状態値を算出する計算量の減少を図って計算
機の応答性を上げている。
存し、互いに反応、相変化、及び熱交換を行いながら移
動する複雑な反応容器である。このような高炉の炉内現
象を、基礎実験、解体調査等の知見を基に、移動速度論
的に整理し、組み合わせることにより、高炉を近似的に
数学モデル化することが可能である。数学モデルとして
は１次元定常モデル（特公昭６１−６１２２号、特公昭
６１−６０１２４号、特公昭６３−２４０４４号、特公
平１−１２８０５号公報）、２次元定常モデル（特公平
６−２８８６号公報）等が開発されてきた。これらのモ
デルでは、高炉の３次元的な挙動を高さ方向のみ取り扱
い（１次元モデル）、又は円周方向の分布を均一とし、
さらに時間的にも状態の変化がないとした定常状態のみ
の解析に限定し（２次元モデル）、モデルを簡略化する
ことで炉内状態値を算出する計算量の減少を図って計算
機の応答性を上げている。
【０００３】
【発明が解決しようとする課題】しかし、高炉では、炉
下部の円周方向に設けられている複数の羽口より熱風を
送風するため、炉壁，充填層等の固体温度は円周方向に
分布し、高温域における反応挙動（ＳｉＯガスのような
気体の固体からの蒸発反応等）に差を生じさせると考え
られるが、上述のような１次元モデルでは高さ方向のみ
の分布を考慮し、また２次元モデルでは周方向の分布を
均一とし、定常状態のみの解析を行って計算を簡略化し
ているので、炉内状態の解析精度が低い。従って、この
ようなモデルによる炉内状態の解析結果を用いた操業シ
ミュレーションで高炉の操業条件を決定した場合、高炉
のスタートアップ時等のように、炉内状態の時間経過に
ともなう変化が大きい非定常状態における高炉の的確な
操業条件を推定できないという問題がある。
下部の円周方向に設けられている複数の羽口より熱風を
送風するため、炉壁，充填層等の固体温度は円周方向に
分布し、高温域における反応挙動（ＳｉＯガスのような
気体の固体からの蒸発反応等）に差を生じさせると考え
られるが、上述のような１次元モデルでは高さ方向のみ
の分布を考慮し、また２次元モデルでは周方向の分布を
均一とし、定常状態のみの解析を行って計算を簡略化し
ているので、炉内状態の解析精度が低い。従って、この
ようなモデルによる炉内状態の解析結果を用いた操業シ
ミュレーションで高炉の操業条件を決定した場合、高炉
のスタートアップ時等のように、炉内状態の時間経過に
ともなう変化が大きい非定常状態における高炉の的確な
操業条件を推定できないという問題がある。
【０００４】本発明はこのような問題点を解決するため
になされたものであって、高炉の高さ方向、半径方向の
みならず、周方向を含む３次元方向の非定常状態を記述
する３次元非定常モデルを用いて高炉の操業をシミュレ
ートし、実炉での計測値と一致するように３次元非定常
モデルのプロセス定数及び装入物分布モデルによって算
出される炉頂での鉄鉱石及びコークスのそれぞれの堆積
層厚さ等の分布状態を修正し、修正後のモデルにより高
炉の操業を高精度にシミュレートすることにより、実炉
の操業条件が的確に決定できる高炉の操業方法の提供を
目的とする。
になされたものであって、高炉の高さ方向、半径方向の
みならず、周方向を含む３次元方向の非定常状態を記述
する３次元非定常モデルを用いて高炉の操業をシミュレ
ートし、実炉での計測値と一致するように３次元非定常
モデルのプロセス定数及び装入物分布モデルによって算
出される炉頂での鉄鉱石及びコークスのそれぞれの堆積
層厚さ等の分布状態を修正し、修正後のモデルにより高
炉の操業を高精度にシミュレートすることにより、実炉
の操業条件が的確に決定できる高炉の操業方法の提供を
目的とする。
【０００５】
【課題を解決するための手段】本発明の高炉の操業方法
は、その壁面の周方向に複数個設けられている羽口から
それぞれ炉内に吹き込まれた熱風が羽口から炉芯にかけ
て堆積しているコークスを通過して炉頂に上昇し、炉頂
から鉄鉱石とともに装入されているコークスをガス化し
て鉄鉱石を還元するとともに、還元された鉄鉱石を溶融
し、溶銑が炉芯のコークス層を通って滴下し、炉床に湯
溜りを形成すべくなしてある高炉の、時間経過に伴って
変化する非定常状態を気相・液相・固相それぞれの物質
収支、運動量収支、及びエネルギー収支に関連付けて表
した状態モデルを用いて推測して非定常状態における実
炉の操業条件を決定する高炉の操業方法において、炉頂
における鉄鉱石及びコークスの分布状態を表した装入物
分布モデルにより鉄鉱石及びコークスの炉頂における粒
子径及び堆積層厚さの分布状態を算出し、炉内の高さ方
向、径方向、及び周方向における炉内の非定常状態を気
相・液相・固相それぞれの物質収支、運動量収支、及び
エネルギー収支に関する微分方程式群により表した３次
元状態モデルに、時間経過に伴ってその値が変化し得る
還元反応速度定数、炉心の空隙率、固気間熱伝達係数及
び気液接触面積を含むプロセス定数を与え、前記装入物
分布モデルにより算出された前記分布状態を炉頂におけ
る境界条件として、また、羽口への送風条件を羽口にお
ける境界条件として与えて高炉の操業をシミュレート
し、銑鉄生産速度、炉頂ガス組成分布、炉頂温度分布、
炉内圧力分布、溶銑温度及び、溶銑中Ｓｉ濃度等を含
む、高炉の操業結果の良否に影響を及ぼす要因の将来の
時間変化状態を推定し、該要因の時間変化状態を実炉に
て測定し、推定値と実測値とが一致するように、３次元
状態モデルの前記プロセス定数と、前記装入物分布モデ
ルにより算出された前記分布状態とを修正した後、前記
要因の時間変化状態が、目的の操業結果が得られる時間
変化状態となるように、修正後の３次元状態モデルによ
る高炉の操業シミュレーションを基にして装入物分布と
羽口への送風条件とを推定し、推定した装入物分布及び
羽口への送風条件を実炉の操業条件とすることを特徴と
する。
は、その壁面の周方向に複数個設けられている羽口から
それぞれ炉内に吹き込まれた熱風が羽口から炉芯にかけ
て堆積しているコークスを通過して炉頂に上昇し、炉頂
から鉄鉱石とともに装入されているコークスをガス化し
て鉄鉱石を還元するとともに、還元された鉄鉱石を溶融
し、溶銑が炉芯のコークス層を通って滴下し、炉床に湯
溜りを形成すべくなしてある高炉の、時間経過に伴って
変化する非定常状態を気相・液相・固相それぞれの物質
収支、運動量収支、及びエネルギー収支に関連付けて表
した状態モデルを用いて推測して非定常状態における実
炉の操業条件を決定する高炉の操業方法において、炉頂
における鉄鉱石及びコークスの分布状態を表した装入物
分布モデルにより鉄鉱石及びコークスの炉頂における粒
子径及び堆積層厚さの分布状態を算出し、炉内の高さ方
向、径方向、及び周方向における炉内の非定常状態を気
相・液相・固相それぞれの物質収支、運動量収支、及び
エネルギー収支に関する微分方程式群により表した３次
元状態モデルに、時間経過に伴ってその値が変化し得る
還元反応速度定数、炉心の空隙率、固気間熱伝達係数及
び気液接触面積を含むプロセス定数を与え、前記装入物
分布モデルにより算出された前記分布状態を炉頂におけ
る境界条件として、また、羽口への送風条件を羽口にお
ける境界条件として与えて高炉の操業をシミュレート
し、銑鉄生産速度、炉頂ガス組成分布、炉頂温度分布、
炉内圧力分布、溶銑温度及び、溶銑中Ｓｉ濃度等を含
む、高炉の操業結果の良否に影響を及ぼす要因の将来の
時間変化状態を推定し、該要因の時間変化状態を実炉に
て測定し、推定値と実測値とが一致するように、３次元
状態モデルの前記プロセス定数と、前記装入物分布モデ
ルにより算出された前記分布状態とを修正した後、前記
要因の時間変化状態が、目的の操業結果が得られる時間
変化状態となるように、修正後の３次元状態モデルによ
る高炉の操業シミュレーションを基にして装入物分布と
羽口への送風条件とを推定し、推定した装入物分布及び
羽口への送風条件を実炉の操業条件とすることを特徴と
する。
【０００６】
【作用】図１は本発明の高炉の操業方法（以下、本発明
方法という）の手順の概略を示すフローチャート、図２
は本発明方法を実施する高炉の模式図である。炉内の気
相・液相・固相の状態を記述する物質収支・運動量収支
・エネルギー収支の数式からなり、高炉の周方向を含む
３次元方向の非定常状態を推定する３次元非定常モデル
（後に詳述する）に、鉄鉱石の還元反応速度定数・炉芯
の空隙率・固気間熱伝達係数・気液接触面積等、時間経
過に伴ってその値が変化し得るプロセス定数と、鉱石／
コークスの落下高さ，落下速度，鉱石／コークスの装入
具である大ベルの角度，装入物と大ベル面との摩擦係数
等を導入した所定の計算式からなる装入物分布モデルに
より算出された、炉頂における鉄鉱石・コークスのそれ
ぞれの粒子径、堆積層厚さ等の分布状態と、解析対象の
実炉における羽口からの送風条件とを与えて解析対象の
高炉の操業をシミュレートし、高炉の操業結果の良否に
影響を及ぼす要因である、銑鉄生産速度・炉頂のガス組
成分布・炉頂温度分布・炉内の圧力分布・溶銑温度・溶
銑中Ｓｉ濃度を推定する（Ｓ１〜Ｓ３）。一方、これら
の銑鉄生産速度・炉頂のガス組成分布・炉頂温度分布・
炉内の圧力分布・溶銑温度・溶銑中Ｓｉ濃度を、ガス分
析計、垂直ゾンデ、圧力計、水平ゾンデ、温度計等の計
測端により実炉において計測する（Ｓ４）。
方法という）の手順の概略を示すフローチャート、図２
は本発明方法を実施する高炉の模式図である。炉内の気
相・液相・固相の状態を記述する物質収支・運動量収支
・エネルギー収支の数式からなり、高炉の周方向を含む
３次元方向の非定常状態を推定する３次元非定常モデル
（後に詳述する）に、鉄鉱石の還元反応速度定数・炉芯
の空隙率・固気間熱伝達係数・気液接触面積等、時間経
過に伴ってその値が変化し得るプロセス定数と、鉱石／
コークスの落下高さ，落下速度，鉱石／コークスの装入
具である大ベルの角度，装入物と大ベル面との摩擦係数
等を導入した所定の計算式からなる装入物分布モデルに
より算出された、炉頂における鉄鉱石・コークスのそれ
ぞれの粒子径、堆積層厚さ等の分布状態と、解析対象の
実炉における羽口からの送風条件とを与えて解析対象の
高炉の操業をシミュレートし、高炉の操業結果の良否に
影響を及ぼす要因である、銑鉄生産速度・炉頂のガス組
成分布・炉頂温度分布・炉内の圧力分布・溶銑温度・溶
銑中Ｓｉ濃度を推定する（Ｓ１〜Ｓ３）。一方、これら
の銑鉄生産速度・炉頂のガス組成分布・炉頂温度分布・
炉内の圧力分布・溶銑温度・溶銑中Ｓｉ濃度を、ガス分
析計、垂直ゾンデ、圧力計、水平ゾンデ、温度計等の計
測端により実炉において計測する（Ｓ４）。
【０００７】３次元非定常モデルにおいて、時間の経過
とともに変化する流体の流れ状況、固体の接触状況等に
よりその値が変化し得るプロセス定数、例えば、鉄鉱石
の還元反応速度定数、炉芯の空隙率、固気間熱伝達係
数、気液接触面積等と、装入物分布モデルにより算出し
た堆積層厚さとを、前述の両モデルによる推定値と実炉
での実測値とが一致するように修正する（Ｓ５）。その
後、修正した３次元非定常モデルを用いて、操業結果に
影響を及ぼす要因である溶銑温度、溶銑中Ｓｉ濃度、炉
内の圧力分布、コークスの燃料比の時間的変化等の非定
常状態が目標の操業結果が得られる非定常状態となる装
入物分布及び羽口への送風条件を推測し（Ｓ６）、これ
らの条件を実炉の操業条件とする（Ｓ７）。
とともに変化する流体の流れ状況、固体の接触状況等に
よりその値が変化し得るプロセス定数、例えば、鉄鉱石
の還元反応速度定数、炉芯の空隙率、固気間熱伝達係
数、気液接触面積等と、装入物分布モデルにより算出し
た堆積層厚さとを、前述の両モデルによる推定値と実炉
での実測値とが一致するように修正する（Ｓ５）。その
後、修正した３次元非定常モデルを用いて、操業結果に
影響を及ぼす要因である溶銑温度、溶銑中Ｓｉ濃度、炉
内の圧力分布、コークスの燃料比の時間的変化等の非定
常状態が目標の操業結果が得られる非定常状態となる装
入物分布及び羽口への送風条件を推測し（Ｓ６）、これ
らの条件を実炉の操業条件とする（Ｓ７）。
【０００８】
【実施例】図３は本発明方法に用いる高炉３次元非定常
モデルの構成概要を示す模式図である。羽口から炉芯に
かけて逆Ｖ字状に堆積しているコークスのうち、羽口に
近いレースウェイのコークスが羽口からの送風によって
燃焼・運動し、羽口から送風された熱風とコークスとの
反応熱によりガス温度が上昇する。ガスはコークス層の
表面層を通って炉頂に上昇し、炉頂から交互に装入され
ている鉱石及びコークスに付着している水分を蒸発させ
て鉱石・コークスを乾燥、除々に加熱し、さらにコーク
スをＨ2 ＯとＣＯ2 とによりガス化し、ガス化により発
生するＨ2 とＣＯとによって鉱石中の鉄と結合している
酸素の大半が除かれる（間接還元）。間接還元された鉱
石は熱と荷重により半溶融状になり、隣り合った粒子同
志が溶着し、逆Ｖ字状に堆積しているコークス層の表面
に融着帯が形成される。一旦溶着した鉱石はより高温の
環境下でさらに還元され（直接還元）、溶銑は逆Ｖ字状
に堆積しているコークスの表面層（滴下帯）を浸炭反応
を起こしながら、炉芯コークスを通過して炉底に滴下
し、湯溜りを形成する。
モデルの構成概要を示す模式図である。羽口から炉芯に
かけて逆Ｖ字状に堆積しているコークスのうち、羽口に
近いレースウェイのコークスが羽口からの送風によって
燃焼・運動し、羽口から送風された熱風とコークスとの
反応熱によりガス温度が上昇する。ガスはコークス層の
表面層を通って炉頂に上昇し、炉頂から交互に装入され
ている鉱石及びコークスに付着している水分を蒸発させ
て鉱石・コークスを乾燥、除々に加熱し、さらにコーク
スをＨ2 ＯとＣＯ2 とによりガス化し、ガス化により発
生するＨ2 とＣＯとによって鉱石中の鉄と結合している
酸素の大半が除かれる（間接還元）。間接還元された鉱
石は熱と荷重により半溶融状になり、隣り合った粒子同
志が溶着し、逆Ｖ字状に堆積しているコークス層の表面
に融着帯が形成される。一旦溶着した鉱石はより高温の
環境下でさらに還元され（直接還元）、溶銑は逆Ｖ字状
に堆積しているコークスの表面層（滴下帯）を浸炭反応
を起こしながら、炉芯コークスを通過して炉底に滴下
し、湯溜りを形成する。
【０００９】次に、本発明方法に用いる３次元非定常モ
デルを記述する支配方程式である気液固３相それぞれの
物質・運動量・エネルギーの収支式について具体的に説
明する。なお、式中の添え字、ｇはガスを、ｓは固体
を、またｌは液体を表す。 (a) 気相の物質収支 気相の成分ｋの質量収支は式(1) で表される。
デルを記述する支配方程式である気液固３相それぞれの
物質・運動量・エネルギーの収支式について具体的に説
明する。なお、式中の添え字、ｇはガスを、ｓは固体
を、またｌは液体を表す。 (a) 気相の物質収支 気相の成分ｋの質量収支は式(1) で表される。
【００１０】
【数１】
【００１１】さらに、式(1) を成分ｋについて加算し、
気相に関する連続の式として式(2)を得る。
気相に関する連続の式として式(2)を得る。
【００１２】
【数２】
【００１３】(b) 気相のエネルギー収支 気相側のエネルギー収支は式(3) で表される。
【００１４】
【数３】
【００１５】(c) 気相の運動量収支 気体が充填粒子層を通過する抵抗が充填層における気体
の圧力損失をモデル化したERGUN 式に従うと仮定する
と、気相における運動量収支は式(4) で表される。
の圧力損失をモデル化したERGUN 式に従うと仮定する
と、気相における運動量収支は式(4) で表される。
【００１６】
【数４】
【００１７】(d) 固相の物質収支 固体粒子として、焼結鉱とコークスとを考える。焼結鉱
及びコークス粒子内の成分ｋの物質収支は式(5) で表さ
れる。但し、ｎは焼結鉱及びコークスを表す。
及びコークス粒子内の成分ｋの物質収支は式(5) で表さ
れる。但し、ｎは焼結鉱及びコークスを表す。
【００１８】
【数５】
【００１９】成分ｋについて式(5) を加算して式(6) を
得る。
得る。
【００２０】
【数６】
【００２１】さらに、式(6) を成分ｎについて加算し、
固相に関する連続の式として式(7)を得る。
固相に関する連続の式として式(7)を得る。
【００２２】
【数７】
【００２３】(e) 固相のエネルギー収支 コークスと焼結鉱とに温度差がないと仮定した場合、固
相のエネルギー収支は式(8) で表される。
相のエネルギー収支は式(8) で表される。
【００２４】
【数８】
【００２５】(f) 固相の運動量収支 固相の運動方程式として、本実施例ではKinematic モデ
ルを用いる。Kinematic モデルの支配方程式は、重力方
向をｙ方向、流速ベクトルの成分を（Ｕs ，Ｖ s ，
Ｗs ）とした場合、式(9) で表される。
ルを用いる。Kinematic モデルの支配方程式は、重力方
向をｙ方向、流速ベクトルの成分を（Ｕs ，Ｖ s ，
Ｗs ）とした場合、式(9) で表される。
【００２６】
【数９】
【００２７】Ｖs については、式(9) を固相の連続の式
である式(7) に代入して求める。 (g) 粒子体積に関する保存式 ガス化や、溶銑への浸炭反応等の固体粒子の消滅反応を
界面反応として記述すれば、粒子１個の体積をＶpnとし
た場合、粒子体積変化はラグランシュ微分を用いて式(1
0)のように表される。
である式(7) に代入して求める。 (g) 粒子体積に関する保存式 ガス化や、溶銑への浸炭反応等の固体粒子の消滅反応を
界面反応として記述すれば、粒子１個の体積をＶpnとし
た場合、粒子体積変化はラグランシュ微分を用いて式(1
0)のように表される。
【００２８】
【数１０】
【００２９】(h) 液相の物質収支 液相の成分ｋの物質収支は式(11)で表される。
【００３０】
【数１１】
【００３１】さらに、成分ｋについて式(11)を加算する
と式(12)が得られる。
と式(12)が得られる。
【００３２】
【数１２】
【００３３】(i) 液相のエネルギー収支 固体及び気体との熱交換を考慮した液相のエネルギー収
支は式(13)で表される。なお、液相は分散相であるか
ら、熱拡散は考慮しないものとする。
支は式(13)で表される。なお、液相は分散相であるか
ら、熱拡散は考慮しないものとする。
【００３４】
【数１３】
【００３５】(j) 液相の運動量収支 充填層内を滴下する流体は垂直滴下を仮定し、運動量収
支として滴下速度を与える。
支として滴下速度を与える。
【００３６】図４は以上のような式(1) 〜(13)で表され
る３次元非定常モデルを解析して炉内の３次元非定常状
態を推定する手順を示すフローチャートである。モデル
は完全自立型であり、装入物の炉頂での条件、羽口への
送風条件、及び炉壁での熱的条件を与えると、実炉運転
と基本的に同じ高炉の操業をシミュレートし、銑鉄生産
速度、及び３次元方向での炉頂のガス組成分布、炉頂温
度分布、炉内の圧力分布、溶銑温度、溶銑中Ｓｉ濃度等
を算出する。
る３次元非定常モデルを解析して炉内の３次元非定常状
態を推定する手順を示すフローチャートである。モデル
は完全自立型であり、装入物の炉頂での条件、羽口への
送風条件、及び炉壁での熱的条件を与えると、実炉運転
と基本的に同じ高炉の操業をシミュレートし、銑鉄生産
速度、及び３次元方向での炉頂のガス組成分布、炉頂温
度分布、炉内の圧力分布、溶銑温度、溶銑中Ｓｉ濃度等
を算出する。
【００３７】即ち、上述の式(1) 〜(13)を適切な境界条
件と初期条件との下に連成し、時間進展することによ
り、各格子点の３次元方向の非定常状態が算出される。
このとき、格子系にはスタガード格子を採用し、任意形
状領域の処理は、境界適合格子による差分法を用いた。
解析アルゴリズムは流体解析で一般的な流速と圧力との
同時緩和法であるSOLA法を改良して用いた。図４に示し
た順に、ガス流れモデル、固体流れモデル、ガス相温度
Ｔg ・固相温度Ｔs ・液相温度Ｔl の連立解、ガス相の
質量分率ｘgk、液相の質量分率ｘlk、固相の質量分率ｘ
sk、及び粒子径を求める。以上を格子点毎に繰り返し、
３次元方向の炉内状態が算出される。
件と初期条件との下に連成し、時間進展することによ
り、各格子点の３次元方向の非定常状態が算出される。
このとき、格子系にはスタガード格子を採用し、任意形
状領域の処理は、境界適合格子による差分法を用いた。
解析アルゴリズムは流体解析で一般的な流速と圧力との
同時緩和法であるSOLA法を改良して用いた。図４に示し
た順に、ガス流れモデル、固体流れモデル、ガス相温度
Ｔg ・固相温度Ｔs ・液相温度Ｔl の連立解、ガス相の
質量分率ｘgk、液相の質量分率ｘlk、固相の質量分率ｘ
sk、及び粒子径を求める。以上を格子点毎に繰り返し、
３次元方向の炉内状態が算出される。
【００３８】図５(a) は羽口近傍を３次元解析した場合
の周方向の各断面におけるガス流れと温度分布との状態
を示す図であり、図５(b) は羽口近傍の液相及び固相に
おけるガス流れと温度分布との解析結果を２次元解析と
比較して示した図である。図から明らかなように、３次
元解析では羽口近傍の固相における周方向の温度分布が
解析されている。従って、高温域における反応挙動（Ｓ
ｉＯガスのような気体の固体からの蒸発反応）の解析精
度が高いと考えられる。
の周方向の各断面におけるガス流れと温度分布との状態
を示す図であり、図５(b) は羽口近傍の液相及び固相に
おけるガス流れと温度分布との解析結果を２次元解析と
比較して示した図である。図から明らかなように、３次
元解析では羽口近傍の固相における周方向の温度分布が
解析されている。従って、高温域における反応挙動（Ｓ
ｉＯガスのような気体の固体からの蒸発反応）の解析精
度が高いと考えられる。
【００３９】また、図６は本発明方法に用いた３次元非
定常モデルにより高炉のスタートアップ時における液相
及びガス相での非定常な流れ及び温度を解析した結果を
示す図である。このように、急激な変化を伴う非定常な
炉内挙動についても安定に計算することができ、実炉で
の計測値（温度の推移、炉頂ガス組成等）との一致も良
いことが分かった。なお、167 時間以上では炉内の状態
変数はほとんど変化しなくなるので、定常解が得られ
る。
定常モデルにより高炉のスタートアップ時における液相
及びガス相での非定常な流れ及び温度を解析した結果を
示す図である。このように、急激な変化を伴う非定常な
炉内挙動についても安定に計算することができ、実炉で
の計測値（温度の推移、炉頂ガス組成等）との一致も良
いことが分かった。なお、167 時間以上では炉内の状態
変数はほとんど変化しなくなるので、定常解が得られ
る。
【００４０】
【発明の効果】以上のように、本発明方法は、高炉の高
さ方向、半径方向のみならず、周方向を含む３次元方向
の非定常状態を記述する３次元非定常モデルを用いて高
炉の操業をシミュレートし、実炉での計測値と一致する
ように３次元非定常モデルのプロセス定数及び装入物分
布モデルによって算出される炉頂での鉄鉱石及びコーク
スのそれぞれの堆積層厚さ等の分布状態を修正し、修正
後のモデルにより高炉の操業を高精度にシミュレートす
るので、実炉の操業条件が的確に決定できるという優れ
た効果を奏する。
さ方向、半径方向のみならず、周方向を含む３次元方向
の非定常状態を記述する３次元非定常モデルを用いて高
炉の操業をシミュレートし、実炉での計測値と一致する
ように３次元非定常モデルのプロセス定数及び装入物分
布モデルによって算出される炉頂での鉄鉱石及びコーク
スのそれぞれの堆積層厚さ等の分布状態を修正し、修正
後のモデルにより高炉の操業を高精度にシミュレートす
るので、実炉の操業条件が的確に決定できるという優れ
た効果を奏する。
【図１】本発明方法の手順を示すフローチャートであ
る。
る。
【図２】本発明方法を実施する高炉の模式図である。
【図３】本発明方法に用いる３次元非定常モデルの構成
概要を示す模式図である。
概要を示す模式図である。
【図４】本発明方法に用いる３次元非定常モデルの解析
手順を示すフローチャートである。
手順を示すフローチャートである。
【図５】２次元解析と３次元解析とによる解析結果の比
較を示す図である。
較を示す図である。
【図６】本発明方法に用いた３次元非定常モデルによる
高炉のスタートアップ時の計算例（気体及び液体の流れ
と温度）を示す図である。
高炉のスタートアップ時の計算例（気体及び液体の流れ
と温度）を示す図である。
Claims (1)
- 【請求項１】 その壁面の周方向に複数個設けられてい
る羽口からそれぞれ炉内に吹き込まれた熱風が羽口から
炉芯にかけて堆積しているコークスを通過して炉頂に上
昇し、炉頂から鉄鉱石とともに装入されているコークス
をガス化して鉄鉱石を還元するとともに、還元された鉄
鉱石を溶融し、溶銑が炉芯のコークス層を通って滴下
し、炉床に湯溜りを形成すべくなしてある高炉の、 時間経過に伴って変化する非定常状態を気相・液相・固
相それぞれの物質収支、運動量収支、及びエネルギー収
支に関連付けて表した状態モデルを用いて推測して非定
常状態における実炉の操業条件を決定する高炉の操業方
法において、 炉頂における鉄鉱石及びコークスの分布状態を表した装
入物分布モデルにより鉄鉱石及びコークスの炉頂におけ
る粒子径及び堆積層厚さの分布状態を算出し、 炉内の高さ方向、径方向、及び周方向における炉内の非
定常状態を気相・液相・固相それぞれの物質収支、運動
量収支、及びエネルギー収支に関する微分方程式群によ
り表した３次元状態モデルに、 時間経過に伴ってその値が変化し得る還元反応速度定
数、炉心の空隙率、固気間熱伝達係数及び気液接触面積
を含むプロセス定数を与え、前記装入物分布モデルによ
り算出された前記分布状態を炉頂における境界条件とし
て、また、羽口への送風条件を羽口における境界条件と
して与えて高炉の操業をシミュレートし、銑鉄生産速
度、炉頂ガス組成分布、炉頂温度分布、炉内圧力分布、
溶銑温度及び、溶銑中Ｓｉ濃度等を含む、高炉の操業結
果の良否に影響を及ぼす要因の将来の時間変化状態を推
定し、 該要因の時間変化状態を実炉にて測定し、 推定値と実測値とが一致するように、３次元状態モデル
の前記プロセス定数と、前記装入物分布モデルにより算
出された前記分布状態とを修正した後、 前記要因の時間変化状態が、目的の操業結果が得られる
時間変化状態となるように、修正後の３次元状態モデル
による高炉の操業シミュレーションを基にして装入物分
布と羽口への送風条件とを推定し、 推定した装入物分布及び羽口への送風条件を実炉の操業
条件とすることを特徴とする高炉の操業方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7106498A JP3033466B2 (ja) | 1995-04-28 | 1995-04-28 | 高炉の操業方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7106498A JP3033466B2 (ja) | 1995-04-28 | 1995-04-28 | 高炉の操業方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH08295910A JPH08295910A (ja) | 1996-11-12 |
| JP3033466B2 true JP3033466B2 (ja) | 2000-04-17 |
Family
ID=14435107
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP7106498A Expired - Lifetime JP3033466B2 (ja) | 1995-04-28 | 1995-04-28 | 高炉の操業方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3033466B2 (ja) |
Families Citing this family (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AT407993B (de) | 1999-03-03 | 2001-07-25 | Voest Alpine Ind Anlagen | Verfahren zur optimierung von auslegung und betrieb eines reduktionsverfahrens |
| JP4968388B2 (ja) * | 2010-02-05 | 2012-07-04 | Ｊｆｅスチール株式会社 | 流体系の温度推定方法、流体系における物質成分の濃度および温度推定方法、流体系の温度分布推定方法、流体系の温度モニタリング方法、溶融金属設備の溶融金属温度制御方法、溶融金属設備の濃度および温度推定方法、溶融亜鉛めっき鋼板ならびに流体系の温度推定装置 |
| CN102667881B (zh) | 2010-03-30 | 2013-11-27 | 新日铁住金系统集成株式会社 | 信息处理装置、信息处理方法以及程序 |
| JP4984000B1 (ja) * | 2010-09-30 | 2012-07-25 | Ｊｆｅスチール株式会社 | 流体系の温度推定方法、流体系の温度分布推定方法、流体系の温度分布モニタリング方法、温度推定装置、溶融亜鉛めっきポット内の溶融亜鉛温度制御方法、溶融亜鉛めっき鋼板、およびタンディッシュ内の溶鋼温度制御方法 |
| JP5583615B2 (ja) * | 2011-02-23 | 2014-09-03 | 新日鐵住金株式会社 | 高炉操業特徴量演算装置、高炉操業特徴量演算方法及びプログラム |
| JP5583614B2 (ja) * | 2011-02-23 | 2014-09-03 | 新日鐵住金株式会社 | 高炉操業状態制御システム、操業制御装置、高炉操業状態制御方法およびプログラム |
| JP5862470B2 (ja) * | 2012-06-13 | 2016-02-16 | 新日鐵住金株式会社 | 高炉休風方法 |
| CN102703627B (zh) * | 2012-06-19 | 2013-10-30 | 首钢总公司 | 一种预测高炉煤气利用率的方法 |
| JP6248550B2 (ja) * | 2013-11-01 | 2017-12-20 | 新日鐵住金株式会社 | 高炉操業条件の決定方法 |
| KR101482469B1 (ko) * | 2013-12-24 | 2015-01-13 | 주식회사 포스코 | 고로 공정 가상 현실 교육 시스템 |
| JP6311659B2 (ja) * | 2015-06-22 | 2018-04-18 | Ｊｆｅスチール株式会社 | 高炉内の層厚分布の推定方法、高炉の操業方法、および高炉内の層厚分布の推定装置 |
| JP6493447B2 (ja) * | 2016-08-02 | 2019-04-03 | Ｊｆｅスチール株式会社 | 溶銑温度予測方法、溶銑温度予測装置、高炉の操業方法、操業ガイダンス装置、溶銑温度制御方法、及び溶銑温度制御装置 |
| JP6531782B2 (ja) * | 2016-08-02 | 2019-06-19 | Ｊｆｅスチール株式会社 | 溶銑温度予測方法、溶銑温度予測装置、高炉の操業方法、操業ガイダンス装置、溶銑温度制御方法、及び溶銑温度制御装置 |
| JP7207118B2 (ja) * | 2019-04-11 | 2023-01-18 | 日本製鉄株式会社 | 装入方法決定方法、装入方法決定装置及び装入方法決定プログラム |
| JP7335500B2 (ja) * | 2019-09-17 | 2023-08-30 | 日本製鉄株式会社 | パラメータ決定装置、パラメータ決定プログラム及びパラメータ決定方法 |
| JP7464033B2 (ja) * | 2021-01-13 | 2024-04-09 | Ｊｆｅスチール株式会社 | 高炉融着層の通気性推定方法、高炉融着層の通気性推定装置、及び高炉の操業方法 |
| CN114058749B (zh) * | 2021-11-29 | 2023-02-10 | 山西晋煤集团技术研究院有限责任公司 | 一种无烟块煤在高炉中的劣化研究方法 |
-
1995
- 1995-04-28 JP JP7106498A patent/JP3033466B2/ja not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JPH08295910A (ja) | 1996-11-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3033466B2 (ja) | 高炉の操業方法 | |
| Yang et al. | Numerical modeling of reaction and flow characteristics in a blast furnace with consideration of layered burden | |
| Kitaev et al. | Heat exchange in shaft furnaces | |
| Fu et al. | CFD modeling of multiphase reacting flow in blast furnace shaft with layered burden | |
| Austin et al. | Analysis of actual blast furnace operations and evaluation of static liquid holdup effects by the four fluid model | |
| CN106957935A (zh) | 一种高炉内部软熔带形状的软测量方法 | |
| Wen et al. | Fluidized‐bed heat transfer: A generalized dense‐phase correlation | |
| JP5583614B2 (ja) | 高炉操業状態制御システム、操業制御装置、高炉操業状態制御方法およびプログラム | |
| Pintowantoro et al. | Numerical analysis of static holdup of fine particles in blast furnace | |
| Zhou et al. | Computational fluid dynamics analysis of 3D hot metal flow characteristics in a blast furnace hearth | |
| CN116306375A (zh) | 高炉生产运行状态的数字孪生方法 | |
| Jindal et al. | A reduced order thermo-chemical model for blast furnace for real time simulation | |
| Mondal et al. | Computational Study of Gas‐Flow and Temperature Dynamics at Blast Furnace Charging | |
| Fu | Numerical simulation of ironmaking blast furnace shaft | |
| JP2727563B2 (ja) | 高炉操業方法 | |
| JP5583615B2 (ja) | 高炉操業特徴量演算装置、高炉操業特徴量演算方法及びプログラム | |
| Panjkovic et al. | Computational Fluid Dynamics Modelling of Iron Flow and Heat Transfer in the Iron Blast Furnace Hearth | |
| Pomerleau et al. | A first principle simulator of an iron oxide pellet induration furnace–an application to optimal tuning | |
| RU2025495C1 (ru) | Способ контроля теплообмена в доменной печи | |
| Nisioka et al. | Development of Mathematical Models for Blast Furnaces | |
| JPS6055561B2 (ja) | 高炉の操業方法 | |
| WO1997012064A9 (en) | A method for the determination of the gas flux distribution in a blast furnace | |
| Dash et al. | Sinter Strand Control; A Simulation Investment | |
| Yang et al. | Numerical modelling of temperature distribution and pressure drop through the layered burden loading in a blast furnace | |
| Hekkala et al. | Mathematical model of heat and mass transfer in the steel belt sintering process |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080218 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090218 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100218 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100218 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110218 Year of fee payment: 11 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120218 Year of fee payment: 12 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120218 Year of fee payment: 12 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130218 Year of fee payment: 13 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130218 Year of fee payment: 13 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140218 Year of fee payment: 14 |
|
| EXPY | Cancellation because of completion of term |