JPH09969A - 微粉炭燃焼ボイラ用ミルの制御方法および装置 - Google Patents

Abstract

(57)【要約】 【目的】 負荷スケジュールどおりの運用を可能とし、
最適タイミングでミル投入できる微粉炭燃焼用ミルの制
御装置を提供する。 【構成】 各時刻における負荷を入力する負荷入力部１
と、負荷スケジュールに基づきミルからの出炭を開始す
る給炭機自動時刻を算出するミル投入時刻演算部２と、
起動前のミル内残炭量とその粒度分布、ミルへの給炭量
の時間変化、給炭の粒径、破砕性等のデータをミルプロ
セスデータ部３から入力してこれらデータに基づきミル
への給炭開始から出炭可能となるまでの出炭待ち時間予
測値を算出するミルシミュレーションモデル部４と、給
炭機自動時刻とミル出炭待ち時間に基づきミルへの給炭
開始時刻を算出する演算部５と、給炭開始時刻と給炭機
自動時刻に基づきミルを制御する制御部６とを備えた微
粉炭燃焼ボイラ用ミルの制御装置。

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は、微粉炭燃焼ボイラ用ミ
ルの制御方法および装置に係り、特に石炭粉砕用ミルを
使用する微粉炭燃焼ボイラを備えた発電プラントの制御
装置であって、あらかじめ決められた負荷スケジュール
どおりに運用するのに好適な微粉炭燃焼ボイラ用ミルの
制御方法および装置に関する。

【０００２】

【従来の技術】まず、石炭粉砕用ミルに石炭を供給する
給炭機起動について（ミル投入）説明する。ミルのウォ
ーミングが完了し、ミルの出口空気温度が設定温度（例
えば６５℃）以上となった段階でミルを起動し、引き続
き給炭機を起動して給炭を開始する。その際いきなり目
標給炭量まで給炭せず、図１のｂに示すように最低給炭
量まで給炭後、一定時間待って給炭機を自動投入して
（図１のａ）目標給炭量まで上げていく。これは給炭を
開始しても石炭はミル内部で循環し、分級機を通り抜け
ることができるほどの微粉となって出炭が安定するまで
に時間遅れが生じ、この間目標給炭まで給炭を増加させ
るとミル内の保有炭量が多くなって粉砕できなくなる
（いわゆるハラボテ現象）ためである。

【０００３】次に、ミルの出炭待ち時間について説明す
る。図１に示したようにこの最低給炭量で時間待ちをす
る時間を出炭待ち時間という。この時間は石炭の水分や
粉砕性（ＨＧＩ；ハード・グルーブ・インデックス）に
よって大きく変化し（１〜５分）予測は困難であった。
そこで従来は使用する炭種の中で、粉砕性の最も悪い炭
種に合わせた一律の待ち時間（５分程度）としていた。

【０００４】次に、ミルに対する給炭開始について説明
する。図２に示すようにボイラ負荷を変化する場合、ボ
イラ負荷により必要なミル台数は決まるので、ミルを投
入しなければならない負荷に到達すると対象となるミル
に給炭を開始（図２のｂ）していた。

【０００５】

【発明が解決しようとする課題】上記従来技術はミル起
動に対し炭種銘柄によって異なる出炭待ち時間を考慮し
ていなかったので下記に示す問題点があった。すなわ
ち、（１）負荷によって決定されるミル投入タイミング
で給炭を開始（図２のｂ）するので、実際の出炭は遅れ
（図２のａ）、真に必要なタイミングで投入されたこと
にならない（ａ−ｂの時間だけ遅れる）。したがってボ
イラの蒸気圧力・温度偏差が大きくなる。（２）炭種銘
柄に無関係な出炭待ち時間とするため、余裕を含んだ大
きな設定となる（例えば５分）。このため、ボイラの高
負荷変化率制御に対応できない。例えばボイラ負荷５０
％から１００％負荷へ５％／分での負荷上昇を想定した
場合、約１０分で負荷上昇が完了してしまうので、出炭
待ち時間を５分もとると、ミル２台の投入も困難とな
り、その分負荷定値となり、負荷スケジュール遅延とな
る。

【０００６】本発明の目的は、上記従来技術の問題点を
解決し、負荷スケジュールどおりの運用を可能とし、か
つ最適タイミングミル投入を行うことができる微粉炭燃
焼ボイラ用ミルの制御方法および装置を提供することで
ある。

【０００７】

【課題を解決するための手段】上記目的を達成するため
本願で特許請求される発明は以下のとおりである。 （１）負荷スケジュールに合わせて微粉炭燃焼ボイラ用
ミルを制御する方法において、各時刻における目標負荷
スケジュールを入力する工程と、入力された負荷スケジ
ュールに従ってミルから出炭を開始する給炭機自動時刻
を算出する工程と、起動前のミル内残炭量とその粒度分
布、ミルへの給炭量の時間変化、給炭の粒度分布および
粉砕性に基づき算出したミル内保有炭量、ミル出炭量お
よび粒度分布の値が規定範囲で安定したときに起動完了
と判断してミル起動後の出炭待ち時間を算出する工程
と、上記給炭機自動時刻と出炭待ち時間とに基づき給炭
機、ダンパを含むミルの制御を行う工程とを有すること
を特徴とする微粉炭燃焼ボイラ用ミルの制御方法。

【０００８】（２）負荷スケジュールに合わせて微粉炭
燃焼ボイラ用ミルを制御する装置において、各時刻にお
ける負荷スケジュールを入力する負荷入力部と、上記負
荷スケジュールに基づきミルからの出炭を開始する給炭
機自動時刻を算出するミル投入時刻演算部と、起動前の
ミル内残炭量とその粒度分布、ミルに対する給炭量の時
間変化、給炭の粒度、粉砕性などのデータをミルシミュ
レーションモデルに入力するミルプロセス入力部と、上
記入力データに基づきミルに給炭を開始してから出炭可
能となるまでの出炭待ち時間予測値を算出するミルシミ
ュレーションモデル部と、上記給炭機自動時刻とミル出
炭待ち時間に基づきミルへの給炭開始時刻を算出する給
炭開始時刻演算部と、該演算部からの給炭開始時刻と前
記給炭機自動時刻に基づきミルを制御する制御部とを備
えたことを特徴とする微粉炭燃焼ボイラ用ミルの制御装
置。

【０００９】

【作用】従来はミルを投入しなければならない負荷に到
達すると給炭を開始（図２のｂ）していたが、これに代
わりこのポイントを給炭自動時刻（図２のａ）とする。
そして出炭待ち時間を予測し、 給炭開始時刻（図２のｂ）＝給炭自動時刻（図２のａ）
−出炭待ち時間 から、給炭開始時刻を計算する。そしてこのスケジュー
ルに従って各時刻に到達したらミルを投入すること（給
炭開始／自動）により、ミル投入の必要な負荷において
目標給炭でのミル投入が可能となるので、スケジュール
の遅延および出炭不足による蒸気温度・圧力の低下がな
い。

【００１０】出炭待ち時間の予測は下記にて行う。 （１）ミル起動前のミル内残炭量、当該残炭の粒度分布
を初期条件として設定する。 （２）ミルへの給炭の粒度分布、粉砕性（ハード・グル
ーブ指数）をパラメータとして設定する。またミル起動
時のミルへの給炭量の時間変化を境界条件として仮定す
る（ミル起動操作法に応じて与える）。 （３）ミル内の粉砕、分級、混合、蓄積の各課題におけ
る微粉炭流量、粒度分布の変化を上記（１）、（２）で
与えた初期条件、パラメータ、境界条件を用いて解く
（時間応答を求める）。 （４）上記（３）で求めたミル内保有炭量、ミル出炭
量、粒度分布の値が規定範囲で安定した段階でミル起動
完了（ミルへの給炭量をプラント全体の制御の観点から
逐次に操作できる状態となること）と判断し、当該状態
への所要時間をもって出炭待ち時間とする。

【００１１】

【実施例】本発明の一実施例を図３により説明する。図
３において１は入力部であり、オペレータにより負荷ス
ケジュール（時刻と負荷）を入力する入力表示装置であ
る。２はミル投入時刻演算部であり、ミル給炭自動時刻
を演算する。５はミルへの給炭開始時刻演算部であり、
ミルへの給炭開始時刻を演算する。６はミル制御部であ
り、ミル投入時刻演算部２および給炭開始時刻演算部５
の指令を受け、ミル関係のダンパ、給炭機等を制御する
部分である。７は石炭粉砕ミルであり、ミルモータコン
タクタや給炭機コンタクタなどのミル本体の電気接点を
含む。３はミルプロセス入力部であり、ミルモデル４の
演算に必要なミルプロセスデータの入力装置である。４
はミルモデルであり、ミル出炭待ち時間を演算する部で
ある。このシミュレーションモデルは、粒度分布をわず
か４つの変数で模擬可能で、低計算量で高精度を得られ
る点に特徴を有する。なお、本発明に適用可能な石炭粉
砕機の動特性シミュレーションモデルの演算法の内容を
以下に述べる。演算に使用する記号、略号などは次のと
おりである。 １．１ 添字 下記の添字は特定の部位、機構、関数、時点を示す。

【００１２】ｒ ：時刻ｔ＝τ△ｔ ｂ ：混合機構 ｉｂ：ミル入口（給炭） ｉｊ：第ｊ番目の分級機構入口 ｉｐ：粉砕機構入口 ｊ ：第ｊ番目の分級機構 ｏｂ：混合機構出口 ｏｐ：粉砕機構出口 ｏｊ：第ｊ番目の分級機構による循環流 ｓ ：粉砕分布定数 τｊ：第ｊ番目の分級機構による通過流 １．２ 記号 β_k ：ｋ次キュムラント μ ：平均 ξ ：対数粒径（ノミナルな値） σ ：標準偏差 ｂ（ξ） ：粉砕速度定数 ｃ（ξ） ：部分分級効率 ｆ（ξ） ：一般に確率密度関数 ｇ（ξ） ：粒度分布確率密度関数 ｇ（ξ｜η） ：条件付確率密度関数 Ξ ：粒子の対数粒径（確率変数） ｒ ：循環比（ｒ≡分級機捕集流量）／分級
機入口流量） ｓ（ξ） ：粉砕分布定数 ｔ ：時刻 ｖ_k （λ，ρ）：（λ，ρ）に規格化したｋ次モーメン
ト Ｅ｛・｝ ：期待値 Ｅ｛・｜・｝ ：条件付期待値 Ｖａｒ｛・｝ ：分散 Ｐｒ｛・｝ ：確率 Ｑ ：重量流量（確率変数） Θ ：分級過程である粒子が通過（Θ＝
１）、または捕集（Θ＝０）捕集の事象を示すインジケ
ータ（確率変数） １．３ 語句 確率変数 ：「ゆらぎ」や情報の不足のため、確率的
な変数値の評価を要す（値を確定的に知り得ない；雑
音、粉体流量、石炭性状が典型的）変数。したがって、
確率Ｐｒ｛ｘ＜Ｘ≦ｘ＋△ｘ｝（確率変数Ｘの値がｘの
近傍となる確率）をｆ_x （ｘ）△ｘ（ｆ_x （ｘ）：Ｘの
確率密度関数）で与え、必要な統計量を計算する。

【００１３】キュムラント：特性関数（確率密度関数の
フーリエ変換）の対数をローラン展開した際、第ｋ次の
展開係数をｋ次キュムラントと呼ぶ。応用上、「モーメ
ントと相互変換可能（簡単な計算）」、「重畳積分（粉
砕機構での粒度分布変化）が加算に帰着」等、便利な性
質がある。

【００１４】Ｓｋｅｗｎｅｓｓ：歪度と訳され、分布の
左右非対称性の程度を示す（＝０で左右対称）。 Ｅｘｃｅｓｓ ：光度と訳され、分布の中心部、裾の
部分の割合を示す（正規分布では＝３）。 条件付確率密度 ：２つの事象の因果関係を考えるとき
一方の情報を得たときの他方の確率密度を示す。具体的
には粉砕機構入口で粒径がη＜Ξ_ip≦η＋ｄηの粒子が
粉砕後にξ＜Ξ_op≦ξ＋ｄξになる確率はΞ_opのΞ_ipに
対する確率密度関数（ｃｏｎｄｉｔｉｏｎａｌ ｐｒｏ
ｂａｂｉｌｉｔｙ ｄｅｎｓｉｔｙ ｆｕｎｃｔｉｏｎ
ｏｆ Ξ_op ｇｉｖｅｎ ｂｙΞ_ip）を用いてｇ
_op|ip（ξ｜η）ｄξと表わされる。 （１）粉砕機内の現象 （１．１）粒度分布の表記 断面を微小時間に通過する粒子中、粒径ξ以下なる質量
割合により粒度分布が定義可能で、その密度関数をｇ
（ξ）と表記し、適宜に場所を示す添字を付加する。サ
ンプルされた静止状態の質量粒度分布密度ｆ（ξ）との
関係は質量流量Ｑを用いて次式となる。

【００１５】

【数１】 ｇ（ξ）≡Ｅ｛Ｑ｜（ξ，ξ＋ｄξ）｝ｆ（ξ）／Ｅ｛Ｑ｝ （１） （１．２）粉砕機構 粉砕前後の諸量にそれぞれ添字ｉｐ、ｏｐを与えると、
粒度分布について次の関係がある。

【００１６】

【数２】

【００１７】ここに、粒径ξを対数軸にとると条件付確
率密度ｇｏｐ ｉｐは、Ｌ．Ａｕｓｔｉｎらの解明した
粉砕分布定数（Ｐｏｗｅｒ Ｔｅｃｈｎｏｌｏｇｙ，Ｖ
ｏｌ．２９，ｐｐ．２６３−２７５（１９８１），同、
Ｖｏｌ．３３，ｐｐ．１１３−１２５（１９８２），
同、Ｖｏｌ．３３，ｐｐ．１２７−１３４（１９８２）
に記載）と一致し、これをｓとする。

【００１８】

【数３】 ｇ_op|ip（ξ｜η）＝ｓ（ξ−η） （３） 質量流量については、粉砕機構内で蓄積はないと仮定し
て次式を得る。

【００１９】

【数４】 Ｅ｛Ｑ_op｝＝Ｅ｛Ｑ_ip｝ （４） （１．３）分級機構 第ｊ番目の分級機構について、各「粒子の通過」は互い
に独立事象であって、Θ_j をインジケータとすれば、実
験により解明されている分級効率ｃ_j （ξ）（廻ほか、
粉体工学誌、Ｖｏｌ．２５，ｐｐ．４３０−４３６（１
９８８）に記載）と次の関係がある。

【００２０】

【数５】 Ｐ_r ｛Θ_j ＝０｜（ξ，ξ＋ｄξ）｝＝ｃ_j （ξ） （５）

【数６】 Ｐ_r ｛Θ_j ＝１｜（ξ，ξ＋ｄξ）｝＝１−ｃ_j （ξ） （６） 分級機入口炭、循環炭、および通過炭に係る諸量にそれ
ぞれ添字ｉｊ、ｒｊ、ｏｊを与えると、ベイズ定理によ
り粒度分布密度を表わす次式を得る。

【００２１】

【数７】 ｇ_rj（ξ）＝ｃ_j （ξ）ｇ_ij（ξ）／ｒ_j （７）

【００２２】

【数８】 ｇ_oj（ξ）＝〔１−ｃ_j （ξ）〕ｇ_ij（ξ）／（１−ｒ_j ） （８） ここに、

【００２３】

【数９】

【００２４】分級機構周辺の流量は次のとおり求められ
る。

【００２５】

【数１０】 Ｅ｛Ｑ_rj｝＝Ｅ｛Θ_j Ｑ_ij｝＝ｒ_j Ｅ｛Ｑ_ij｝ （１０）

【００２６】

【数１１】 Ｅ｛Ｑ_oj｝＝Ｅ｛（１−Θ_j ）Ｑ_ij｝＝（１−ｒ_j ）Ｅ｛Ｑ_ij｝ （１１） （１．４）混合機構 分級機構（ｊ＝０、……、ｎ）からの循環炭と、原料炭
（添字ｉｂ）とを混合して流出炭（添字ｏｂ）となす機
構を考える。ここで混合機構保有炭Ｇ_b とＱ_obとの間に
次の関係を仮定する。

【００２７】

【数１２】 Ｅ｛Ｑ_ob｝＝Ｅ｛Ｐ｝Ｅ｛Ｇ_b ｝ （１２） Ｐは粒径と独立とし、この仮定を正当化するため、混合
機構と続く粉砕機構の間に仮想的な分級機構（ｊ＝０）
を設けて前出のＬ．Ａｕｓｔｉｎらが解明したξに依存
する粉砕速度定数を考慮する。

【００２８】ここで、前出の第（１）式と第（１２）式
に着目し、混合により粒径は変化しないと考えて、
（ξ，ξ＋ｄξ）に属する粒子のマスバランス式が得ら
れる。

【００２９】

【数１３】

【００３０】（２）モデルの数学的記述 （２．１）分布密度のパラメタライズ Ξが分布密度ｇ（ξ）に従うときλ、ρで規準化（アフ
イン変換）したモーメントを考える。

【００３１】

【数１４】 ｖ_k （λ，ρ）＝Ｅ｛〔（Ξ−λ）／ρ〕^k ｝ （１４） このとき、キユムラントβ_k （λ，ρ）が対応して求め
られる。本モデルでは、分布密度を次の４パラメータで
整理する。

【００３２】

【数１５】 μ＝ｖ₁ （０，１），σ＝〔ｖ₂ （０，１）〕^1/2 （１５）

【００３３】

【数１６】 Ｓｋｅｗｎｅｓｓ：β₃ （μ，σ）＝ｖ₃ （μ，σ） （１６）

【００３４】

【数１７】 Ｅｘｃｅｓｓ：β₄ （μ，σ）＝ｖ₄ （μ，σ）^-3 （１７） これらより一意にエッジワース展開係数α_k が求まり、
分布密度を具体的に表示できる。

【００３５】

【数１８】 ｇ（ξ）＝Σ_k α_k ｐ（ξ；μ，σ）ｈ_k （〔ξ−μ〕／σ） （１８） ここにｐ（ξ；μ，σ）はガウス分布、ｈ_k はｋ次のエ
ルミート多項式である。 （２．２）粉砕機構 前出の第（３）式を第（２）式に代入すると重畳積分で
あって、キユムラントの和に帰着し、以下を得る。

【００３６】

【数１９】 ｖ_1op （μ_{ip ,}σ_ip）＝ｖ_1s（μ_{ip ,}σ_ip） （１９）

【００３７】

【数２０】 ｖ_2op （μ_{ip ,}σ_ip）＝１＋ｖ_2s（μ_{ip ,}σ_ip） （２０）

【００３８】

【数２１】 ｖ_3op （μ_{ip ,}σ_ip）＝ｖ_3ip （μ_{ip ,}σ_ip）＋ｖ_3s（μ_{ip ,}σ_ip） （２１）

【００３９】

【数２２】 ｖ_4op （μ_{ip ,}σ_ip） ＝ｖ_4ip （μ_{ip ,}σ_ip）＋ｖ_4s（μ_{ip ,}σ_ip）＋６ｖ_2s（μ_{ip ,}σ_ip） （２２） ここに、添字ｓは粉砕分布定数ｓを、それ以外は各粒度
分布密度ｇを指す。さらに、μ_op、σ_op、β₃ （μ_{op ,}
σ_op）、β₄ （μ_{op ,}σ_op）は前出の第（１５）式〜第
（１７）式および次式を用いて計算できる。

【００４０】

【数２３】

【００４１】（２．３）分級機構 ｃ_j （ξ）は、適当なτ_mj、λ_mj、ρ_mjを用いて近似で
きる。

【００４２】

【数２４】 ｃ_j （ξ）≒Σ_m τ_mjｐ（ξ；λ_{mj ,}ρ_mj） （２４） ｇ_ij（ξ）は第（１８）式の形式であり、第（７）式お
よび第（１０）式より循環炭の諸量が具体的に求められ
る。

【００４３】

【数２５】 Ｅ｛Ｑ_rj｝ｇ_rj（ξ） ＝Ｅ｛Ｑ_ij｝Σ_m ｐ（ξ；λ_{mrj ,} ρ_mrj ） ・Σ_k α_rjmkｈ_k （〔ξ−λ_mrj 〕／ρ_mrj ） （２５） ここに次の関係がある。

【００４４】

【数２６】 α_rjmk＝α_0rjmk τ_mj ・ｅｘｐ｛−（μ_ij−λ_mj）² ／〔２（σ_ij ² ＋ρ_mj ² ）〕｝ （２６）

【００４５】

【数２７】 λ_mrj ＝（ρ_mj ² μ_ij＋σ_ij ² λ_mj）／（σ_ij ² ＋ρ_mj ² ） （２７）

【００４６】

【数２８】 ρ_mrj ＝ρ_mjσ_ij／（σ_ij ² ＋ρ_mj ² ）^1/2 （２８） またα_0rjmk は、次式にエルミート多項式の加法定理を
適用し、係数を整理して得られる。

【００４７】

【数２９】 Σ_k α_ijk ｈ_k （〔ρ_mrj ／σ_ij〕〔ξ−λ_mrj 〕／ρ_mrj ＋〔λ_mrj −μ_ij〕／σ_ij） （２９） 前出の第（２５）式は分布密度の重みつき混合であり、
添字ｍについてα_{0rjm k} からｖ_k （λ_{mrj ,} ρ_mrj ）が
一意に求まり、同一λ、ρのｖ_k は重みつき加算が可能
だから、結局前出の第（２３）式、第（１５）式〜第
（１７）式を用いて、μ_rj、σ_rj、β₃ （μ
_{rj ,}σ_rj）、β₄ （μ_{rj ,}σ_rj）が計算できる。添字ｏ
ｊの通過炭についても同様の議論である。 （２．４）混合機構 適当に選んだλ_b 、ρ_b で規準化すると、前出の第（１
３）式よりｖ_kob についての微分方程式を得る。

【００４８】

【数３０】

【００４９】第（２３）式、第（１５）式〜第（１７）
式を適用すれば、一般にμ、σ、β ₃ （μ，σ）、β₄
（μ，σ）とｖ_k （λ，ρ）の相互変換が可能だから、
（３．２）、（３．３）の結論を代入して（３０）を解
くことができる。このときＰａｄｅ近似の採用で安定な
数値計算が可能となった。したがって、初期条件として
Ｅ｛Ｑ_ob｝、ｖ_kob （λ_b ，ρ_b ）を与え、起動操作手
順に応じてミルへの給炭量とその粒度分布に対応するＥ
｛Ｑ_ib｝、ｖ_kib（λ_{ib ,}ρ_ib）を与え、以上の諸式を
逐次解くことにより、ミル出口出炭量とその粒度分布で
あるＥ｛Ｑ_oj｝、ｖ_koj （λ_{ib ,}ρ_ib）および保有炭量
Ｅ｛Ｇｂ｝を求めることができる。

【００５０】以上述べた演算を行うミルモデル４を使っ
て図３に示すようにミルの制御を行う。入力部１からオ
ペレータが時刻と目標負荷をキーイン（入力）し、投入
時刻演算部２ではキーインされた負荷スケジュールから
給炭器自動時刻を演算する。一方、ミルプロセス入力部
３により入力されたミルのプロセスデータからミルモデ
ル４によりミル出炭待ち時間を前述した演算式を使って
計算する。次いで給炭開始時刻演算部５では給炭自動時
刻とミル出炭待ち時間からミル給炭開始時刻を計算す
る。計算された給炭開始時刻および給炭器自動時刻によ
りミル制御部６では必要なダンパや給炭器等を含めてミ
ル７を制御する。

【００５１】ミル動特性モデルは前述のとおりミル内の
物理現象を忠実に模擬しているので、起動停止のような
大幅な状態変化に際しても高精度である。本動特性モデ
ルの動作例を図４を用いて説明する。給炭器を自動にし
て（図４のａ）から目標給炭量まで給炭量を増加させる
が、この時の出炭もミル内での循環のため遅れることと
なり実出炭が目標給炭量となるまでの遅れを出炭遅れ時
間という。

【００５２】この出炭遅れ時間を前記と同じ方法で予測
し、目標給炭量到達時刻（図４のｃ）をミル投入ポイン
トとし、給炭自動時刻を逆算し、スケジューリングする
ことにより、さらにボイラの蒸気温度・圧力に対する安
定化の効果がある。

【００５３】

【発明の効果】本発明によればミル出炭時間を正確に予
測し、ボイラ負荷に対し必要なポイント（タイミング）
でミルを投入できるので、プラントの負荷スケジュール
に対し遅延なく運用できるとともに、ボイラの蒸気温度
・圧力の要求に対し不足なく石炭を供給でき、プラント
の良好な運転を行うことができる。

【図面の簡単な説明】

【図１】給炭量と出炭量の関係を示す説明図。

【図２】ボイラ負荷変化に対し、給炭開始と給炭自動の
遅延を説明する図。

【図３】本発明の実施例全体構成図。

【図４】他の実施例の説明図。

【符号の説明】

１…負荷スケジュール入力部、２…ミル投入時刻演算
部、３…ミルプロセス入力部、４…ミルモデル、５…給
炭時刻演算部、６…ミル制御部、７…ミル。

フロントページの続き (72)発明者 井山 望 東京都中央区銀座６丁目15番１号 電源開 発株式会社内 (72)発明者 山本 晃二 広島県呉市宝町６番９号 バブコック日立 株式会社呉工場内 (72)発明者 深山 幸穂 広島県呉市宝町３番36号 バブコック日立 株式会社呉研究所内 (72)発明者 津村 俊一 広島県呉市宝町６番９号 バブコック日立 株式会社呉工場内 (72)発明者 川瀬 隆世 広島県呉市宝町６番９号 バブコック日立 株式会社呉工場内

Claims (2)

【特許請求の範囲】
1. 【請求項１】 負荷スケジュールに合わせて微粉炭燃焼
   ボイラ用ミルを制御する方法において、各時刻における
   目標負荷スケジュールを入力する工程と、入力された負
   荷スケジュールに従ってミルから出炭を開始する給炭機
   自動時刻を算出する工程と、起動前のミル内残炭量とそ
   の粒度分布、ミルへの給炭量の時間変化、給炭の粒度分
   布および粉砕性に基づき算出したミル内保有炭量、ミル
   出炭量および粒度分布の値が規定範囲で安定したときに
   起動完了と判断してミル起動後の出炭待ち時間を算出す
   る工程と、上記給炭機自動時刻と出炭待ち時間とに基づ
   き給炭機、ダンパを含むミルの制御を行う工程とを有す
   ることを特徴とする微粉炭燃焼ボイラ用ミルの制御方
   法。
2. 【請求項２】 負荷スケジュールに合わせて微粉炭燃焼
   ボイラ用ミルを制御する装置において、各時刻における
   負荷スケジュールを入力する負荷入力部と、上記負荷ス
   ケジュールに基づきミルからの出炭を開始する給炭機自
   動時刻を算出するミル投入時刻演算部と、起動前のミル
   内残炭量とその粒度分布、ミルに対する給炭量の時間変
   化、給炭の粒度、粉砕性などのデータをミルシミュレー
   ションモデルに入力するミルプロセス入力部と、上記入
   力データに基づきミルに給炭を開始してから出炭可能と
   なるまでの出炭待ち時間予測値を算出するミルシミュレ
   ーションモデル部と、上記給炭機自動時刻とミル出炭待
   ち時間に基づきミルへの給炭開始時刻を算出する給炭開
   始時刻演算部と、該演算部からの給炭開始時刻と前記給
   炭機自動時刻に基づきミルを制御する制御部とを備えた
   ことを特徴とする微粉炭燃焼ボイラ用ミルの制御装置。