JP2018044719A - 循環系粉砕プラントにおける粉砕機の入口温度の予測方法、装置及びプログラム、並びに原料供給量の演算方法、装置及びプログラム - Google Patents
循環系粉砕プラントにおける粉砕機の入口温度の予測方法、装置及びプログラム、並びに原料供給量の演算方法、装置及びプログラム Download PDFInfo
- Publication number
- JP2018044719A JP2018044719A JP2016179904A JP2016179904A JP2018044719A JP 2018044719 A JP2018044719 A JP 2018044719A JP 2016179904 A JP2016179904 A JP 2016179904A JP 2016179904 A JP2016179904 A JP 2016179904A JP 2018044719 A JP2018044719 A JP 2018044719A
- Authority
- JP
- Japan
- Prior art keywords
- inlet temperature
- pulverizer
- raw material
- amount
- heat
- 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.)
- Granted
Links
Images
Landscapes
- Disintegrating Or Milling (AREA)
Abstract
Description
特に循環系粉砕プラントでは、加熱された排ガスがライン内を循環するため、排ガスを全て排気する1パス式の粉砕プラントに比べると、粉砕機の入口温度が高くなる傾向がある。また、排ガス循環によって排ガス排出量を抑制できる構造から粉砕量を大型化することがあり、やはり粉砕機の入口温度が高くなる傾向にある。そのため、循環系粉砕プラントにおいては、粉砕機の入口温度を高精度に予測し、運転条件を適正化することへの要求が1パス式の粉砕プラントよりも高いと言える。
また、循環系粉砕プラントの操業中においても、粉砕機の入口温度を高精度に予測できれば、プラント制御に有効活用することができる。
特許文献1では、循環系粉砕プラントにおいて熱・物質収支計算を行うことが開示されているが、粉砕機の入口温度を予測することは開示されていない。
また、循環系粉砕プラントにおける粉砕機の入口温度を、1パス式の粉砕プラントにおける粉砕機の入口温度の予測方法に倣って予測することが行われているが、実際の温度に比較して遥かに低い温度を計算してしまうことが多いという実情がある。
[1] 熱風を排ガスとして発生する熱ガス発生装置と、原料を粉砕し、粉砕後の原料を前記熱ガス発生装置で発生する排ガスの流れに乗せて外部に放出する粉砕機と、前記粉砕機から排ガスの流れに乗って放出された粉砕後の原料を捕集する捕集機とを備え、循環ファンにより、前記捕集機からの排ガスを、その一部を大気中に放散するとともに、前記熱ガス発生装置に循環させる循環系粉砕プラントにおいて、前記粉砕機の入口温度を予測する方法であって、
前記粉砕機の入口温度は、前記熱ガス発生装置が有するバーナーにより発生する熱量と、前記循環ファンによる仕事量と、大気中への排ガスの放散による排熱量とを補うものとする熱収支に基づいて、前記粉砕機の入口温度を計算することを特徴とする循環系粉砕プラントにおける粉砕機の入口温度の予測方法。
[2] 大気中への排ガスの放散による排熱量は、放散ガスを外気温から前記粉砕機の出口温度まで加熱するのに必要な熱量として表わすことを特徴とする[1]に記載の循環系粉砕プラントにおける粉砕機の入口温度の予測方法。
[3] 前記バーナーにより発生する熱量ΔQHGGは、燃料ガスカロリーが
一定であるとして、式(101)によりで与えられ、
ΔQHGG=燃料ガスカロリー×燃料ガス流量・・・(101)
前記循環ファンによる仕事量は、既知の値が与えられることを特徴とする[1]又は[2]に記載の循環系粉砕プラントにおける粉砕機の入口温度の予測方法。
[4] 前記粉砕機の入口温度TINを、ガスの比熱C、前記捕集機の出口排ガス流量Fbug、放散ガス流量Fexh、前記粉砕機の出口温度TOUT、前記バーナーにより発生する熱量ΔQHGG、前記循環ファンによる仕事量ΔQFAN、外気温T0を用いて、式(102)により、
C・(Fbug−Fexh)・(TIN―TOUT)=
ΔQHGG+ΔQFAN+C・Fexh・(TOUT−T0)・・・(102)
計算することを特徴とする[1]乃至[3]のいずれか一つに記載の循環系粉砕プラントにおける粉砕機の入口温度の予測方法。
[5] 前記ガスの比熱C、前記循環ファンによる仕事量ΔQFANは、既知の値が与えられ、
前記放散ガス流量Fexh、前記バーナーにより発生する熱量ΔQHGGは、運転条件として設定される、前記捕集機の出口排ガス流量Fbug、前記粉砕機の出口温度TOUT、及び前記外気温T0を含む条件を用いて、熱及び物質収支に基づいて計算されることを特徴とする[4]に記載の循環系粉砕プラントにおける粉砕機の入口温度の予測方法。
[6] 熱風を排ガスとして発生する熱ガス発生装置と、原料を粉砕し、粉砕後の原料を前記熱ガス発生装置で発生する排ガスの流れに乗せて外部に放出する粉砕機と、前記粉砕機から排ガスの流れに乗って放出された粉砕後の原料を捕集する捕集機とを備え、循環ファンにより、前記捕集機からの排ガスを、その一部を大気中に放散するとともに、前記熱ガス発生装置に循環させる循環系粉砕プラントにおいて、前記粉砕機の入口温度を予測する装置であって、
前記粉砕機の入口温度は、前記熱ガス発生装置が有するバーナーにより発生する熱量と、前記循環ファンによる仕事量と、大気中への排ガスの放散による排熱量とを補うものとする熱収支に基づいて、前記粉砕機の入口温度を計算する手段を備えたことを特徴とする循環系粉砕プラントにおける粉砕機の入口温度の予測装置。
[7] 熱風を排ガスとして発生する熱ガス発生装置と、原料を粉砕し、粉砕後の原料を前記熱ガス発生装置で発生する排ガスの流れに乗せて外部に放出する粉砕機と、前記粉砕機から排ガスの流れに乗って放出された粉砕後の原料を捕集する捕集機とを備え、循環ファンにより、前記捕集機からの排ガスを、その一部を大気中に放散するとともに、前記熱ガス発生装置に循環させる循環系粉砕プラントにおいて、前記粉砕機の入口温度を予測するためのプログラムであって、
前記粉砕機の入口温度は、前記熱ガス発生装置が有するバーナーにより発生する熱量と、前記循環ファンによる仕事量と、大気中への排ガスの放散による排熱量とを補うものとする熱収支に基づいて、前記粉砕機の入口温度を計算する処理をコンピュータに実行させるためのプログラム。
前記粉砕機の入口温度と原料供給量との複数の組み合わせに基づいて、前記粉砕機の入口温度が上限を超えない範囲で最大となる原料供給量を求めるステップとを有することを特徴とする循環系粉砕プラントにおける原料供給量の演算方法。
[9] [1]乃至[5]のいずれか一つに記載の循環系粉砕プラントにおける粉砕機の入口温度の予測方法により、前記粉砕機の入口温度を原料供給量を変更して繰り返し計算し、前記粉砕機の入口温度と原料供給量との複数の組み合わせを得る手段と、
前記粉砕機の入口温度と原料供給量との複数の組み合わせに基づいて、前記粉砕機の入口温度が上限を超えない範囲で最大となる原料供給量を求める手段とを備えたことを特徴とする循環系粉砕プラントにおける原料供給量の演算装置。
[10] [1]乃至[5]のいずれか一つに記載の循環系粉砕プラントにおける粉砕機の入口温度の予測方法により、前記粉砕機の入口温度を原料供給量を変更して繰り返し計算し、前記粉砕機の入口温度と原料供給量との複数の組み合わせを得る処理と、
前記粉砕機の入口温度と原料供給量との複数の組み合わせに基づいて、前記粉砕機の入口温度が上限を超えない範囲で最大となる原料供給量を求める処理とをコンピュータに実行させるためのプログラム。
(第1の実施形態)
図1は、高炉への微粉炭の吹込み(PCI:Pulverized Coal Injection)を行うために石炭を粉砕する、負圧式の循環系粉砕プラント(以下、単に循環系粉砕プラントと呼ぶ)の構成例を示す図である。図1において、各構成要素を繋ぐ実線は配管を示し、矢印線は配管内のガスや石炭の進行方向を示す。
燃焼エアファン102は、熱ガス発生装置101に燃焼エアを送り込む。
給炭機104は、チェーンコンベアを有し、バンカー103内に貯蔵されている石炭をチェーンコンベアにより切り出してミル105に投入する。
ミル105は、給炭機104から投入された石炭を粉砕する粉砕機である。ミル105は、例えばロールミル105aと粉砕テーブル105bとを有する。ミル105の上部から投入された石炭をロールミル105aと粉砕テーブル105bとの間に供給する。回転している粉砕テーブル105bに対してロールミル105aを押し付けながら回転させることにより、石炭は押し潰されて粉砕される。粉砕された石炭(以下、微粉炭と称する)は、熱ガス発生装置101から供給された排ガスの流れにのって、ミル105の上部に供給され、分級機で分級された後、外部に放出される。
異物除去装置108は、バグフィルター107で捕集された異物を除去する。
リザーブタンク109は、異物除去装置108で異物が除去された微粉炭を貯蔵する。リザーブタンク109に貯蔵された微粉炭は、高炉の羽口から高炉の内部に吹き込まれる(微粉炭吹き込みが行われる)。
循環ファン113は、ダンパー112を通過した排ガスを熱ガス発生装置101に循環させることができるように、排ガスを昇圧する。
放散塔(煙突)114は、循環ファン113により昇圧された排ガスの一部(放散ガス)を大気中に放散する。
放散系圧力調整弁115は、放散塔114から大気中に放散される排ガスの圧力を調整する。
エア流量調整弁118は、循環系粉砕プラントに供給されるエアの流量を調整する。
希釈エアファン119は、エア流量調整弁118で流量が調整された希釈エアを昇圧し、希釈エアを熱ガス発生装置101の入側の配管に押し込む。これにより、循環ガスの酸素濃度を大きくすることができる。
ミル入口温度は、熱ガス過熱防止のための管理指標として用いられる。ミル入口伸縮管の耐熱温度や石炭の発火温度の関係から、ミル入口温度には上限が設定される。
ミル出口温度は、石炭乾燥のための管理指標として用いられる。ミル出口温度が下がると石炭を十分に乾燥できなくなるため、ミル出口温度は一定にすることが好ましい。ミル出口温度を一定に保つため、ミル出口温度制御を行う。ミル出口温度制御とは、ミル出口温度を目標温度に追従させるためにバーナー負荷を操作するものを指し、ミル出口温度計で測定されたミル出口温度の偏差からフィードバック制御を行う。フィードバック制御としては、PID制御が一般的である。
既述したように、プラントの運転条件に応じて粉砕機の入口温度を高精度に予測できれば、設備仕様を適正化することが可能となる。本実施形態では、プラント設計の事前評価に有効活用できるように、循環系粉砕プラントにおけるミル入口温度を、オフラインの熱・物質収支モデルに基づいて計算する。
<熱・物質収支モデルに基づくガス流量及び熱量の計算>
まず、熱・物質収支モデルに基づくガス(エアを含む)流量及び熱量の計算を行う(ステップS1)。ミル入口温度の計算を行うには、プロセス全体として熱・物質収支を満足するガス流量、熱量が必要である。
そこで、以下に述べる[A]〜[D]のモデルを連立して、BFG流量Fbfg、燃焼エア流量Fcair、希釈エア流量Fpairを計算する。
次のパラメータが与えられたとする。
ミル出口温度(製品温度) :TOUT=90[℃]
給炭量 :[ton/hr]
バグフィルター出口排ガス流量 :Fbug[Nm3/hr]
バグフィルター出口のO2濃度 :10[%]
石炭水分 :[%] ※粉砕前の石炭に含まれる水分
製品水分 :1.5[%] ※粉砕・乾燥後の微粉炭の水分
外気温 :T0[℃]
バンカー侵入エア流量 :Fbair[Nm3/hr]
シールエア流量 :Fsair[Nm3/hr]
一方、バグフィルター出口排ガス流量Fbugは、例えば循環ファン113の吸引能力を上限とした適切な設定値が与えられる。また、ミル出口温度TOUT、バグフィルター出口のO2濃度は、プラントが爆発せず、かつ、プラント内を循環する水蒸気が液滴化して配管に付着・腐食しない範囲での最適な設定値が与えられる。
製品水分は、例えばバグフィルター107で捕集された微粉炭のサンプリングを行って、平均値を設定すればよい。一例としては経験的に1.5[%]が与えられる。
外気温T0は、例えば夏場及び冬場それぞれの温度を設定してもよいし、その中間の温度を設定してもよい。
オフライン計算では、バンカー侵入エア流量Fbair、シールエア流量Fsair、石炭水分の3つが不明量であるが、既知として値を与えた上で計算を行うことにする。例えば石炭水分は、原料サンプリングを行って、その平均値を設定すればよい。
全ての物質は初期温度が外気温T0であり、バーナーが排出する熱風によって全てのガス(水蒸気含む)、微粉炭の温度がミル出口温度TOUTになる。
また、系内に注入、或いは発生するガスは、全て放散塔114から放散され、系内の圧力は一定の微負圧に保たれている。
熱・物質収支モデルに基づくガス流量の計算方法を示す。
[A]放散ガス流量モデル
放散塔114からの放散ガス流量Fexhは、式(1)で表わされる。Fcvは石炭由来の蒸気流量である。
Fexh=Fbfg+Fcair+Fcv+Fsair+Fbair+Fpair・・・(1)
バグフィルター出口のO2濃度は、上記のガス成分によって、式(2)で表わされる。O2濃度目標値は例えば10[%]であるとする。
O2濃度=(Fpair+Fbair+Fsair+Fcair×0.1)×21/Fexh・・・(2)
石炭由来の蒸気流量Fcvと、水蒸気として存在する水の単位時間あたりの重量WVを紐付けることができる。製品水分は例えば1.5[%]であるとする。
Fcv=WV×18/22.4
WV=給炭量×1000×{石炭水分/(100−石炭水分)−製品水分/(100−製品水分)}・・・(3)
[A]〜[C]のモデルに加えて、以下の熱収支式を連立して、BFG流量Fbfg、燃焼エア流量Fcair、希釈エア流量Fpairを計算することができる。
図3に、熱収支モデルのイメージを示す。図3に示すように、バーナーにより発生する熱量ΔQHGG、循環ファン113による仕事量(圧縮熱)ΔQFANについて、式(4)が成立する。
ΔQHGG+ΔQFAN
=ΔQ潜熱+ΔQ顕熱+ΔQ石炭+ΣΔQGAS(i)・・・(4)
ΔQGAS(i)=ガス比熱(i)[kcal/Nm3/℃]×ガス流量(i)×(製品温度−ガス注入温度(i))・・・(5)
具体的には、i=1〜5を考える。i=1はBFGガス、i=2はバーナー燃焼エア、i=3はシールエア、i=4はバンカー侵入エア、i=5は希釈エアであるとすれば、以下のように解釈することができる。
ガス比熱(1)はBFGの比熱、ガス比熱(2)〜(5)は空気の比熱とする。ガス流量は以下のようになる。
ガス流量(1)=Fbfg
ガス流量(2)=Fcair
ガス流量(3)=Fsair
ガス流量(4)=Fbair
ガス流量(5)=Fpair
なお、ガス注入温度(1)〜(5)は外気温T0とする。
ここで、
Fcair=Fbfg×理論空気量×過剰空気量
である。例えば、BFGガスの理論空気量=0.638、過剰空気量=1.1とする。
ΔQ石炭=石炭比熱[kcal/kg/℃]×給炭量×1000×(製品温度−原料温度)・・・(6)
なお、原料温度は外気温T0とする。
ΔQ顕熱=水比熱[kcal/kg/℃]×WM×(製品温度−原料温度)・・・(7)
ここで、石炭持ち込みの水の単位時間あたりの重量WMは、式(8)で表わされる。
WM=給炭量×1000×石炭水分/(100−石炭水分)・・・(8)
なお、原料温度は外気温T0とする。
ΔQ潜熱=水潜熱[kcal/kg]×WV・・・(9)
ΔQHGG=BFGカロリー×Fbfg・・・(10)
ΔQFAN=一定・・・(11)
具体的には、式(4)で記述される熱・物質収支計算と式(1)、式(2)で記述されるO2濃度に関する制約条件を同時に満たすFbfg、Fcair、Fpairを計算することになる。反復計算とは、以下の偏差の絶対値が一定値(例えば、ΔO2濃度が0.01%)以下になるまで、繰り返すことを指す。
ΔQ=ΔQHGG+ΔQFAN−{ΔQ潜熱+ΔQ顕熱+ΔQ石炭+ΣΔQGAS(i)}
ΔO2濃度=O2濃度目標値
−(Fpair+Fbair+Fsair+Fcair×0.1)×21/Fexh
まず、式(4)で希釈エア流量をゼロと仮定してBFG流量・燃焼エア流量を計算する。次に、式(1)に基づいて放散ガス流量Fexhを計算して、式(2)からO2濃度に関する制約条件を満足する希釈エアの流量Fpairを導出する。再び、導出した希釈エアの流量を式(4)に与えて、BFG流量・燃焼エア流量を計算する。これらの反復計算をΔO2濃度の絶対値が一定値(例えば、0.01%)以下になるまで繰り返す。
結果として、熱・物質収支を満足するBFG流量Fbfg、燃焼エア流量Fcair、希釈エア流量Fpairが計算される。
次に、ミル入口温度の計算を行う(ステップS2)。
[A]〜[D]のモデルで計算された熱・物質収支計算の結果に基づいて、式(12)により、ミル入口温度TINを計算する。
ここで、ガスの比熱Cは、既知として扱う。ミル入口温度TIN以外のパラメータは、設定値、或いは[A]〜[D]のモデルで計算された値であるので、式(12)から直ちに「ミル入口温度」を計算することができる。
C・(Fbug−Fexh)・(TIN―TOUT)=
ΔQHGG+ΔQFAN+C・Fexh・(TOUT−T0)・・・(12)
1パス式に倣った計算式は、式(13)で表わされる。Fcirは循環ガス流量、FINはミル入口流量である。
C・FIN・TIN=ΔQHGG+ΔQFAN+C・Fcir・TOUT+C・(Fpair+Fbfg+Fcair)・・・(13)
ただし、
FIN=Fbug−Fexh+Fpair+Fbfg+Fcair
Fcir=Fbug−Fexh
まず、1パス式の場合、次のようにミル入口温度を計算する。図5(a)には、1パス式の粉砕プラントの構造を単純化したモデルを示す。ミル出口温度TOUTに見合ったΔQHGGを計算してから、式(14)のように、単純な熱収支でミル入口温度TINを計算する。
ΔQHGG+ΔQFAN=C・FIN・(TIN−T0)・・・(14)
FIN=冷風流量+燃焼排ガス流量
FIN=Fcir+FHGG+Fpair・・・(15)
ΔQHGG+C・(Fpair+FHGG)・T0+C・Fcir・Tcir
=C.FIN・TIN・・・(16)
C・Fcir・(Tcir−TOUT)=ΔQFAN・・・(17)
これら式を整理すると、式(13)が得られる。
図6(a)には、実運転データにおけるBFG流量、希釈エア流量、給炭量設定変化率制限を示す。
図6(b)に、実際のミル入口温度と、発明法により計算したミル入口温度と、従来法により計算したミル入口温度とを示す。また、図7に、実際のミル入口温度に対する、発明法により計算したミル入口温度と、従来法により計算したミル入口温度との乖離度を示す。
図6(b)、図7に示すように、従来法により計算したミル入口温度は実際の温度に比較して遥かに低い温度を計算してしまい、経験的には、粉砕量が大きくなるほどに実際の温度と乖離し、粉砕量によっては100℃以上の誤差が生じてしまう。それに対して、発明法により計算したミル入口温度は実際の温度の略一致しており、実運転データを良く再現できていることがわかる。
入力部1001は、ミル入口温度の予測に必要な条件、具体的には、ミル出口温度TOUT、給炭量、バグフィルター出口排ガス流量Fbug、バグフィルター出口のO2濃度、石炭水分、製品水分、外気温T0、バンカー侵入エア流量Fbair、シールエア流量Fsairを入力する。
また、入力部1001は、各種の比熱や水潜熱、循環ファン113による仕事量、BFGカロリーを入力する。
これらの情報は、不図示の入力装置を介してオペレータから入力されるようにしてもよいし、ネットワークを介して外部機器から入力されるようにしてもよい。
ミル入口温度計算部1003は、熱・物質収支計算部1002での計算結果を用いて、ステップS2で述べたミル入口温度の計算を行う。
出力部1004は、ミル入口温度計算部1003で計算したミル入口温度を出力する。例えばミル入口温度を、不図示のディスプレイを表示したり、外部機器に送出したりする。
第1の実施形態では、プラント設計の事前評価に有効活用できるように、オフライン計算を行う例を説明した。
第2の実施形態では、循環系粉砕プラントの操業中にプラント制御を行うのに有効活用できるように、オンライン計算を行う例を説明する。具体的には、循環系粉砕プラントにおけるミル入口温度を、オンラインの熱・物質収支モデルに基づいて計算し、ミル入口温度が上限(例えば400℃)を超えない範囲で最大となる給炭量(最大給炭量)を求める。
<石炭水分の推定>
まず、石炭水分を推定する(ステップS11)。ミル入口温度の計算を高精度に行うには、実操業における石炭水分を推定することが求められる。
ここで、図9に示すように、循環系粉砕プラントにおいては、各種ガスの流量と、バグフィルター出口のO2濃度と、ミル入口温度及びミル出口温度とが計測され、石炭水分と、バンカー侵入エア流量Fbairと、シールエア流量Fsairとが未知量となる。
次のパラメータが与えられたとする。
ミル出口温度(製品温度) :実績値TOUT
給炭量 :設定値[ton/hr]
バグフィルター出口排ガス流量 :実績値Fbug[Nm3/hr]
バグフィルター出口のO2濃度 :実績値[%]
製品水分 :1.5[%] ※粉砕・乾燥後の微粉炭の水分
外気温 :実績値T0[℃]
放散ガス :実績値Fexh[Nm3/hr]
BFG流量 :実績値Fbfg[Nm3/hr]
燃焼エア流量 :実績値Fcair[Nm3/hr]
希釈エア流量 :実績値Fpair[Nm3/hr]
放散塔114からの放散ガス流量Fexhは、式(1)で表わされる。
また、バグフィルター出口のO2濃度[%]は、式(2)で表わされる。
式(1)、式(2)においては、石炭由来の蒸気流量Fcvと、バンカー侵入エア流量Fbairと、シールエア流量Fsairとが未知量となる。式(1)、式(2)の連立を解くと、石炭由来の蒸気流量Fcvと、バンカー侵入エア流量Fbairとシールエア流量Fsairの和とを計算することができる。
そして、石炭由来の蒸気流量Fcvが得られると、式(3)から石炭水分を計算することができる。結果として、石炭水分と、バンカー侵入エア流量Fbairとシールエア流量Fsairの和とを計算することができる。
次に、ステップS11で得られた石炭水分と、バンカー侵入エア流量Fbairとシールエア流量Fsairの和とを用いて、熱・物質収支モデルに基づくガス流量及び熱量の計算を行い(ステップS12)、ミル入口温度の計算を行う(ステップS13)。
ステップS12、S13は、第1の実施形態で説明したステップS1、S2と同様であり、ここではその詳細な説明は省略するが、ここでは、給炭量を複数回にわたって変更してミル入口温度の計算を繰り返し、ミル入口温度と給炭量との複数の組み合わせを得る。
パラメータをまとめると、以下のようになる。
・複数回変更するパラメータ
給炭量 :[ton/hr]
・固定パラメータ(設定値)
ミル出口温度(製品温度) :TOUT=90[℃]
バグフィルター出口排ガス流量 :Fbug[Nm3/hr]
バグフィルター出口のO2濃度 :10[%]
製品水分 :1.5[%]
・固定パラメータ(実績値)
外気温 :実績値T0[℃]
・固定パラメータ(推定値)
石炭水分[%]
バンカー侵入エア流量Fbairとシールエア流量Fsairの和[Nm3/hr]
次に、ミル入口温度と給炭量との複数の組み合わせに基づいて、ミル入口温度が上限(例えば400℃)を超えない範囲で最大となる給炭量(最大給炭量)を求める(ステップS14)。
[計算方法]
ステップS12、S13で得られたデータ集合を作成して、Y:ミル入口温度とX:給炭量との近似関数を得て、その近似関数からミル入口温度が400℃となる給炭量を計算する。
具体的には、まず給炭量の設定値「給炭量_k」を与えて、ステップS12、S13でミル入口温度を予測する。このミル入口温度の予測値を「TIN_k」と呼ぶ。なお、ステップS12、S13を簡単に表わすと、式(18)のようになる。式(18)において、φは給炭量を変数としてミル入口温度を計算するための関数を表わす。
TIN_k=φ(給炭量_k)・・・(18)
TIN_1=φ(給炭量_1)
TIN_2=φ(給炭量_2)
TIN_3=φ(給炭量_3)
TIN_4=φ(給炭量_4)
給炭量_1=最小給炭量(給炭機104のベルト速度の最小値で規定される量)、給炭量_4=最大給炭量(給炭機104のベルト速度の最大値で規定される量)として、
給炭量_2=(給炭量_4−給炭量_1)×1/3+給炭量_1
給炭量_3=(給炭量_4−給炭量_1)×2/3+給炭量_1
とする。
Y=A・X2+B・X+C・・・(19)
400=A・X2+B・X+C
を満足する給炭量X(>0)を計算する。
入力部1101は、石炭水分の推定に必要な条件、具体的には、ミル出口温度の実績値TOUT、給炭量、バグフィルター出口排ガス流量の実績値Fbug、バグフィルター出口のO2濃度の実績値、製品水分、外気温の実績値T0、放散ガスの実績値Fexh、BFG流量の実績値Fbfg、燃焼エア流量の実績値Fcair、希釈エア流量の実績値Fpairを入力する。
また、入力部1101は、ミル入口温度の予測に必要な条件、具体的には、複数回変更するパラメータである給炭量、固定パラメータ(設定値)である、製品温度(ミル出口温度)TOUT、バグフィルター出口排ガス流量Fbug、バグフィルター出口のO2濃度、製品水分を入力する。
また、入力部1101は、各種の比熱や水潜熱、循環ファン113による仕事量、BFGカロリーを入力する。
これらの情報は、不図示の入力装置を介してオペレータから入力されるようにしてもよいし、ネットワークを介して外部機器から入力されるようにしてもよい。
熱・物質収支計算部1103は、石炭水分計算部1102で推定された石炭水分、及び入力部1101で入力された条件に基づいて、ステップS12で述べた熱・物質収支モデルに基づくガス流量及び熱量の計算を行う。
ミル入口温度計算部1104は、熱・物質収支計算部1103での計算結果を用いて、ステップS13で述べたミル入口温度の計算を行う。
最大給炭量計算部1105は、熱・物質収支計算部1103及びミル入口温度計算部1104での繰り返し計算により得られたミル入口温度と給炭量との複数の組み合わせに基づいて、ミル入口温度が上限を超えない範囲で最大となる給炭量を求める。
出力部1106は、最大給炭量計算部1105で求めた最大給炭量を出力する。例えば最大給炭量を、不図示のディスプレイを表示する。また、最大給炭量を、給炭量を制御する制御装置に送出し、この最大給炭量に基づいて制御装置が給炭量制御を実行するように構成してもよい。
また、本発明は、本発明の機能を実現するソフトウェア(プログラム)を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータがプログラムを読み出して実行することによっても実現可能である。
Claims (10)
- 熱風を排ガスとして発生する熱ガス発生装置と、原料を粉砕し、粉砕後の原料を前記熱ガス発生装置で発生する排ガスの流れに乗せて外部に放出する粉砕機と、前記粉砕機から排ガスの流れに乗って放出された粉砕後の原料を捕集する捕集機とを備え、循環ファンにより、前記捕集機からの排ガスを、その一部を大気中に放散するとともに、前記熱ガス発生装置に循環させる循環系粉砕プラントにおいて、前記粉砕機の入口温度を予測する方法であって、
前記粉砕機の入口温度は、前記熱ガス発生装置が有するバーナーにより発生する熱量と、前記循環ファンによる仕事量と、大気中への排ガスの放散による排熱量とを補うものとする熱収支に基づいて、前記粉砕機の入口温度を計算することを特徴とする循環系粉砕プラントにおける粉砕機の入口温度の予測方法。 - 大気中への排ガスの放散による排熱量は、放散ガスを外気温から前記粉砕機の出口温度まで加熱するのに必要な熱量として表わすことを特徴とする請求項1に記載の循環系粉砕プラントにおける粉砕機の入口温度の予測方法。
- 前記バーナーにより発生する熱量ΔQHGGは、燃料ガスカロリーが一定であるとして、式(101)によりで与えられ、
ΔQHGG=燃料ガスカロリー×燃料ガス流量・・・(101)
前記循環ファンによる仕事量は、既知の値が与えられることを特徴とする請求項1又は2に記載の循環系粉砕プラントにおける粉砕機の入口温度の予測方法。 - 前記粉砕機の入口温度TINを、ガスの比熱C、前記捕集機の出口排ガス流量Fbug、放散ガス流量Fexh、前記粉砕機の出口温度TOUT、前記バーナーにより発生する熱量ΔQHGG、前記循環ファンによる仕事量ΔQFAN、外気温T0を用いて、式(102)により、
C・(Fbug−Fexh)・(TIN―TOUT)=
ΔQHGG+ΔQFAN+C・Fexh・(TOUT−T0)・・・(102)
計算することを特徴とする請求項1乃至3のいずれか1項に記載の循環系粉砕プラントにおける粉砕機の入口温度の予測方法。 - 前記ガスの比熱C、前記循環ファンによる仕事量ΔQFANは、既知の値が与えられ、
前記放散ガス流量Fexh、前記バーナーにより発生する熱量ΔQHGGは、運転条件として設定される、前記捕集機の出口排ガス流量Fbug、前記粉砕機の出口温度TOUT、及び前記外気温T0を含む条件を用いて、熱及び物質収支に基づいて計算されることを特徴とする請求項4に記載の循環系粉砕プラントにおける粉砕機の入口温度の予測方法。 - 熱風を排ガスとして発生する熱ガス発生装置と、原料を粉砕し、粉砕後の原料を前記熱ガス発生装置で発生する排ガスの流れに乗せて外部に放出する粉砕機と、前記粉砕機から排ガスの流れに乗って放出された粉砕後の原料を捕集する捕集機とを備え、循環ファンにより、前記捕集機からの排ガスを、その一部を大気中に放散するとともに、前記熱ガス発生装置に循環させる循環系粉砕プラントにおいて、前記粉砕機の入口温度を予測する装置であって、
前記粉砕機の入口温度は、前記熱ガス発生装置が有するバーナーにより発生する熱量と、前記循環ファンによる仕事量と、大気中への排ガスの放散による排熱量とを補うものとする熱収支に基づいて、前記粉砕機の入口温度を計算する手段を備えたことを特徴とする循環系粉砕プラントにおける粉砕機の入口温度の予測装置。 - 熱風を排ガスとして発生する熱ガス発生装置と、原料を粉砕し、粉砕後の原料を前記熱ガス発生装置で発生する排ガスの流れに乗せて外部に放出する粉砕機と、前記粉砕機から排ガスの流れに乗って放出された粉砕後の原料を捕集する捕集機とを備え、循環ファンにより、前記捕集機からの排ガスを、その一部を大気中に放散するとともに、前記熱ガス発生装置に循環させる循環系粉砕プラントにおいて、前記粉砕機の入口温度を予測するためのプログラムであって、
前記粉砕機の入口温度は、前記熱ガス発生装置が有するバーナーにより発生する熱量と、前記循環ファンによる仕事量と、大気中への排ガスの放散による排熱量とを補うものとする熱収支に基づいて、前記粉砕機の入口温度を計算する処理をコンピュータに実行させるためのプログラム。 - 請求項1乃至5のいずれか1項に記載の循環系粉砕プラントにおける粉砕機の入口温度の予測方法により、前記粉砕機の入口温度を原料供給量を変更して繰り返し計算し、前記粉砕機の入口温度と原料供給量との複数の組み合わせを得るステップと、
前記粉砕機の入口温度と原料供給量との複数の組み合わせに基づいて、前記粉砕機の入口温度が上限を超えない範囲で最大となる原料供給量を求めるステップとを有することを特徴とする循環系粉砕プラントにおける原料供給量の演算方法。 - 請求項1乃至5のいずれか1項に記載の循環系粉砕プラントにおける粉砕機の入口温度の予測方法により、前記粉砕機の入口温度を原料供給量を変更して繰り返し計算し、前記粉砕機の入口温度と原料供給量との複数の組み合わせを得る手段と、
前記粉砕機の入口温度と原料供給量との複数の組み合わせに基づいて、前記粉砕機の入口温度が上限を超えない範囲で最大となる原料供給量を求める手段とを備えたことを特徴とする循環系粉砕プラントにおける原料供給量の演算装置。 - 請求項1乃至5のいずれか1項に記載の循環系粉砕プラントにおける粉砕機の入口温度の予測方法により、前記粉砕機の入口温度を原料供給量を変更して繰り返し計算し、前記粉砕機の入口温度と原料供給量との複数の組み合わせを得る処理と、
前記粉砕機の入口温度と原料供給量との複数の組み合わせに基づいて、前記粉砕機の入口温度が上限を超えない範囲で最大となる原料供給量を求める処理とをコンピュータに実行させるためのプログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2016179904A JP6702110B2 (ja) | 2016-09-14 | 2016-09-14 | 循環系粉砕プラントにおける粉砕機の入口温度の予測方法、装置及びプログラム、並びに原料供給量の演算方法、装置及びプログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2016179904A JP6702110B2 (ja) | 2016-09-14 | 2016-09-14 | 循環系粉砕プラントにおける粉砕機の入口温度の予測方法、装置及びプログラム、並びに原料供給量の演算方法、装置及びプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2018044719A true JP2018044719A (ja) | 2018-03-22 |
JP6702110B2 JP6702110B2 (ja) | 2020-05-27 |
Family
ID=61694686
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2016179904A Active JP6702110B2 (ja) | 2016-09-14 | 2016-09-14 | 循環系粉砕プラントにおける粉砕機の入口温度の予測方法、装置及びプログラム、並びに原料供給量の演算方法、装置及びプログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6702110B2 (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110110462A (zh) * | 2019-05-15 | 2019-08-09 | 中国神华能源股份有限公司 | 热风始端温度的调控方法 |
CN116116554A (zh) * | 2023-02-10 | 2023-05-16 | 华能曲阜热电有限公司 | 一种冷却风管积粉自动吹扫系统及方法 |
-
2016
- 2016-09-14 JP JP2016179904A patent/JP6702110B2/ja active Active
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110110462A (zh) * | 2019-05-15 | 2019-08-09 | 中国神华能源股份有限公司 | 热风始端温度的调控方法 |
CN116116554A (zh) * | 2023-02-10 | 2023-05-16 | 华能曲阜热电有限公司 | 一种冷却风管积粉自动吹扫系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
JP6702110B2 (ja) | 2020-05-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yu et al. | Multistage relative humidity control strategy enhances energy and exergy efficiency of convective drying of carrot cubes | |
JP6776881B2 (ja) | 粉砕プラントにおける処理装置、方法、およびプログラム | |
JPH0127337B2 (ja) | ||
CN109577064B (zh) | 用于卫生纸机干燥部能耗与蒸发量预测的机理建模方法 | |
CN111365733B (zh) | 一种燃煤锅炉经济氧量动态寻优方法、系统及设备 | |
KR20130018697A (ko) | 석탄 분말을 플랜트에 공급하는 복수의 석탄 밀을 위한 제어 파라미터를 자동 제공하기 위한 컴퓨터-기반 방법 및 장치 | |
CN104919160B (zh) | 燃气轮机燃料供给方法和装置 | |
JP5065653B2 (ja) | 石炭ガス化炉運転制御方法、石炭ガス化炉運転制御装置および石炭ガス化炉運転制御プログラム | |
JP6702110B2 (ja) | 循環系粉砕プラントにおける粉砕機の入口温度の予測方法、装置及びプログラム、並びに原料供給量の演算方法、装置及びプログラム | |
JP6492478B2 (ja) | 粉砕プラントにおける給炭量の演算方法、装置及びプログラム | |
US20130122433A1 (en) | Smart Firing Control in a Rankine Cycle Power Plant | |
JP7427154B2 (ja) | 固体燃料ベースの燃焼プロセスのフィードバック制御を改善するための動的熱発生計算 | |
KR101804477B1 (ko) | 과열기 온도 제어 방법 | |
CN103033065B (zh) | 一种余热发电系统的控制方法及装置 | |
JP6801438B2 (ja) | 粉砕プラント酸素濃度制御装置、粉砕プラント酸素濃度制御方法、およびプログラム | |
JP6347100B2 (ja) | 排ガス循環系粉砕プラントのミル出口温度制御方法、装置及びプログラム | |
CN115345512A (zh) | 燃煤锅炉掺烧污泥耦合发电调控方法、系统、设备及介质 | |
JP3481352B2 (ja) | 微粉炭ボイラの制御装置 | |
JP7120167B2 (ja) | 微粉炭の製造装置及び製造方法 | |
JP2018103065A (ja) | 粉砕プラント酸素濃度制御装置、粉砕プラント酸素濃度制御方法、およびプログラム | |
JP6794948B2 (ja) | 微粉炭の製造装置及び製造方法 | |
JP6702109B2 (ja) | 粉砕プラントにおけるプロセス制御方法、粉砕プラントにおけるプロセス制御装置、およびプログラム | |
WO2008103131A1 (fr) | Procédé de réglage du mode de fonctionnement d'une installation d'extinction à sec du coke et dispositif de mise en oeuvre de ce procédé | |
CN112146887B (zh) | 一种燃煤机组最低稳燃负荷实时计算方法及系统 | |
JP3565616B2 (ja) | 微粉炭燃焼ボイラ用ミルの制御方法および装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20190415 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20200206 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20200218 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20200319 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20200407 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20200420 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 6702110 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |