JP6776881B2 - 粉砕プラントにおける処理装置、方法、およびプログラム - Google Patents
粉砕プラントにおける処理装置、方法、およびプログラム Download PDFInfo
- Publication number
- JP6776881B2 JP6776881B2 JP2016249067A JP2016249067A JP6776881B2 JP 6776881 B2 JP6776881 B2 JP 6776881B2 JP 2016249067 A JP2016249067 A JP 2016249067A JP 2016249067 A JP2016249067 A JP 2016249067A JP 6776881 B2 JP6776881 B2 JP 6776881B2
- Authority
- JP
- Japan
- Prior art keywords
- mill
- inlet temperature
- value
- mill inlet
- temperature
- 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.)
- Active
Links
Images
Landscapes
- Crushing And Grinding (AREA)
- Disintegrating Or Milling (AREA)
- Feedback Control In General (AREA)
Description
本発明は、粉砕プラントにおける処理装置、方法、およびプログラムに関し、特に、粉砕プラントに用いて好適なものである。
微粉炭やセメント等を製造するための粉砕プラントとの一例として、高炉への微粉炭の吹込み（ＰＣＩ；Pulverized Coal Injection）を行うために石炭を粉砕するＰＣＩプラントがある。
ＰＣＩプラントでは、まず、燃料ガス（燃焼ガス）と燃焼エアとを熱ガス発生装置に供給し、熱ガス発生装置において、熱風を排ガスとして発生させる。排ガスは、原料の粉砕を行う粉砕機の内部に供給される。粉砕機で粉砕された原料（粉体）は、排ガスと共にバグフィルター（濾布（繊布や不織布））に供給され、バグフィルターで捕集される。
負圧式・排ガス循環系のＰＣＩプラントでは、熱風（排ガス）は、循環ファンで昇圧されて循環ガスとして再び熱ガス発生装置に供給される。一方、１パス方式のＰＣＩプラントでは、熱風（排ガス）は、循環されずにそのまま煙突から大気中に放散される。
ＰＣＩプラントでは、まず、燃料ガス（燃焼ガス）と燃焼エアとを熱ガス発生装置に供給し、熱ガス発生装置において、熱風を排ガスとして発生させる。排ガスは、原料の粉砕を行う粉砕機の内部に供給される。粉砕機で粉砕された原料（粉体）は、排ガスと共にバグフィルター（濾布（繊布や不織布））に供給され、バグフィルターで捕集される。
負圧式・排ガス循環系のＰＣＩプラントでは、熱風（排ガス）は、循環ファンで昇圧されて循環ガスとして再び熱ガス発生装置に供給される。一方、１パス方式のＰＣＩプラントでは、熱風（排ガス）は、循環されずにそのまま煙突から大気中に放散される。
何れのＰＣＩプラントにおいても、粉砕機の入口側の配管内のガスの温度であるミル入口温度と、粉砕機の出側の配管内の微粉炭の温度であるミル出口温度を管理して操業が行われる。ミル入口温度が高くなり過ぎると、配管を繋ぐ伸縮管が燃焼したり、粉砕機に供給された原料（石炭）が発火したりする虞がある。そこで、ミル入口温度には上限値が設定される。
一方、原料（石炭）を粉砕してできる微粉炭内に含まれる水分を一定に保つために、ミル出口温度は、一定に保たれることが望ましい。そこで、特許文献１に記載されているように、ミル出口温度を一定に保つためのミル出口温度制御が行われる。ミル出口温度制御とは、ミル出口温度を目標温度に追従させるために熱ガス発生装置を構成するバーナーの負荷を操作する制御である。尚、特許文献１では、ヒートアップが終了して給炭を開始するとミル出口温度が目標値に低下するまで、操業条件に応じた一定の流量のバーナー負荷を熱ガス発生装置に与え、ミル出口温度が目標値に低下した後にミル出口温度制御を行う。
以上のようにミル入口温度には上限値があり、且つ、ミル出口温度が一定になるように熱ガス発生装置が操作されるので、ミル入口温度が高くなる場合には、給炭量を下げることになる。給炭量を下げると微粉炭の生産量が低下する。そのため、ミル入口温度が、その上限値を上回らない範囲で給炭量を自動または手動で制御する。ミル入口温度が一定になるように給炭量を制御しても、給炭量が正確に設定値の通りにならないこと等に起因して、ＰＣＩプラントには常に外乱が入る。その結果、ミル入口温度にはばらつきが生じる。ミル入口温度を制御する目的は、ミル入口温度が上限値を上回らないようにすることであるので、一般に、ミル入口温度の目標値をこのようなばらつきを見越した値に設定する手法が採られる。このとき、ミル入口温度のばらつきを小さく見積もりすぎればミル入口温度の目標値は必要以上に高くなり、ミル入口温度が上限値を上回る危険性が高まる。逆に、ばらつきを大きく見積もりすぎればミル入口温度の目標値は必要以上に低くなり、微粉炭の生産量を十分に上げることができない。
そこで、ミル入口温度のばらつきを予測し、このばらつきに応じてミル入口温度の目標値を設定する方法が考えられる。プラントにおける目標値を設定する方法として特許文献２、３に記載の技術がある。
特許文献２には、プラントの制御系を模擬する模擬制御系（モデル）を設け、この模擬制御系に参考軌道に基づき参考目標値を供給して予測応答を求め、求めた予測応答を評価した結果に基づいて参考軌道を修正して最適目標軌道を求め、求めた最適目標軌道に基づき最適目標値を決定することが開示されている。
特許文献２には、プラントの制御系を模擬する模擬制御系（モデル）を設け、この模擬制御系に参考軌道に基づき参考目標値を供給して予測応答を求め、求めた予測応答を評価した結果に基づいて参考軌道を修正して最適目標軌道を求め、求めた最適目標軌道に基づき最適目標値を決定することが開示されている。
特許文献３には、プラントの運転状態ごとに計測信号のばらつき量を求め、求めた計測信号のばらつき量とプラントの運転上の制約値とに基づいてプラントの運転効率が大きくなる制御目標値を求めることが開示されている。
しかしながら、特許文献２に記載の技術では、プラントの制御系を模擬する模擬制御系（モデル）を用いる。前述したように、ミル入口温度のばらつきには外乱に起因するものがある。したがって、特許文献２に記載の技術を用いても、外乱をモデルに表すことが容易ではないため、ミル入口温度のばらつきを正確に予測することが容易ではない。
また、特許文献３に記載の技術では、運転状態ごとの計測信号は定常値を基準として変動することを前提とする。ミル入口温度の定常値は給炭量、原料水分、外気温と強い相関があるが、これらは時々刻々と変化する。したがって、特許文献２に記載の技術を用いても、定常値を正確に定めることが容易ではないため、ミル入口温度のばらつきを正確に予測することが容易ではない。
以上のように従来の技術では、ミル入口温度のばらつきを正確に予測することが容易ではないため、ミル入口温度が上限値を可及的に上回らないようにすることと、原料の供給量を可及的に多くすることを両立させることが容易ではなかった。
以上のように従来の技術では、ミル入口温度のばらつきを正確に予測することが容易ではないため、ミル入口温度が上限値を可及的に上回らないようにすることと、原料の供給量を可及的に多くすることを両立させることが容易ではなかった。
本発明は、以上のような問題点に鑑みてなされたものであり、ミル入口温度が上限値を可及的に上回らないようにすることと、原料の供給量を可及的に多くすることを両立させることを目的とする。
本発明の粉砕プラントにおける処理装置は、熱風を排ガスとして発生する熱風発生装置と、原料を粉砕し、粉砕後の原料を、前記排ガスの流れに乗せて外部に放出する粉砕機と、前記粉砕機に原料を供給する原料供給装置と、前記粉砕機から前記排ガスの流れに乗って放出された粉砕後の原料を捕集する捕集機と、を有し、前記粉砕機の出口側の所定の位置における温度であるミル出口温度の測定値と目標値との偏差に応じて前記熱風発生装置を操作する制御が行われる粉砕プラントに対する処理を行う、粉砕プラントにおける処理装置であって、前記ミル出口温度のばらつき、または、前記ミル出口温度の測定値と目標値との偏差のばらつきをミル出口温度ばらつきとして導出するミル出口温度ばらつき導出手段と、前記粉砕機の入口側の所定の位置における温度であるミル入口温度の上限値と、前記ミル出口温度ばらつき導出手段により導出された前記ミル出口温度ばらつきを前記ミル入口温度のばらつきに換算した値と、を用いて、前記ミル入口温度の設定値を導出するミル入口温度設定値導出手段と、を有し、前記ミル入口温度設定値導出手段により導出された前記ミル入口温度の設定値は、前記ミル入口温度の測定値が前記ミル入口温度の設定値を上回らないように、前記原料供給装置から前記粉砕機に供給される前記原料の単位時間当たりの供給量を導出するために用いられることを特徴とする。
本発明の粉砕プラントにおける処理方法は、熱風を排ガスとして発生する熱風発生装置と、原料を粉砕し、粉砕後の原料を、前記排ガスの流れに乗せて外部に放出する粉砕機と、前記粉砕機に原料を供給する原料供給装置と、前記粉砕機から前記排ガスの流れに乗って放出された粉砕後の原料を捕集する捕集機と、を有し、前記粉砕機の出口側の所定の位置における温度であるミル出口温度の測定値と目標値との偏差に応じて前記熱風発生装置を操作する制御が行われる粉砕プラントに対する処理を行う、粉砕プラントにおける処理方法であって、前記ミル出口温度のばらつき、または、前記ミル出口温度の測定値と目標値との偏差のばらつきをミル出口温度ばらつきとして導出するミル出口温度ばらつき導出工程と、前記粉砕機の入口側の所定の位置における温度であるミル入口温度の上限値と、前記ミル出口温度ばらつき導出工程により導出された前記ミル出口温度ばらつきを前記ミル入口温度のばらつきに換算した値と、を用いて、前記ミル入口温度の設定値を導出するミル入口温度設定値導出工程と、を有し、前記ミル入口温度設定値導出工程により導出された前記ミル入口温度の設定値は、前記ミル入口温度の測定値が前記ミル入口温度の設定値を上回らないように、前記原料供給装置から前記粉砕機に供給される前記原料の単位時間当たりの供給量を導出するために用いられることを特徴とする。
本発明のプログラムは、前記粉砕プラントにおける処理装置の各手段としてコンピュータを機能させることを特徴とする。
本発明によれば、ミル入口温度が上限値を可及的に上回らないようにすることと、原料の供給量を可及的に多くすることを両立させることができる。
以下、図面を参照しながら、本発明の一実施形態を説明する。本実施形態では、粉砕プラントが、負圧式・排ガス循環系のＰＣＩプラントである場合を例に挙げて説明する。尚、後述するように、本実施形態の適用範囲は、負圧式・排ガス循環系のＰＣＩプラントに限定されない。
（負圧式・排ガス循環系のＰＣＩプラントの構成）
図１は、負圧式・排ガス循環系のＰＣＩプラントの構成の一例を示す図である。図１において、各構成要素を繋ぐ実線は配管を示し、破線は信号の伝達経路を示す。また、矢印線は、配管内のガスや石炭の進行方向を示す。尚、負圧式・排ガス循環系のＰＣＩプラントの構成は、例えば、特許文献１に記載の技術等の公知の技術で実現できるので、ここでは、各構成について簡単に説明し、詳細な説明を省略する。
図１は、負圧式・排ガス循環系のＰＣＩプラントの構成の一例を示す図である。図１において、各構成要素を繋ぐ実線は配管を示し、破線は信号の伝達経路を示す。また、矢印線は、配管内のガスや石炭の進行方向を示す。尚、負圧式・排ガス循環系のＰＣＩプラントの構成は、例えば、特許文献１に記載の技術等の公知の技術で実現できるので、ここでは、各構成について簡単に説明し、詳細な説明を省略する。
図１において、熱ガス発生装置（ＨＧＧ）１０１は、バーナーを有し、燃料ガスおよび燃焼エア（空気）をバーナーへの入力として、バーナーの空燃比を制御し、排ガス（熱風）を発生させる。本実施形態では、燃料ガスとしてＢＦＧ（Blast Furnace Gas）を利用する。尚、燃焼エアは、燃焼エアファン１０２により熱ガス発生装置１０１に送り込まれる。
バンカー１０３は、原料である石炭を貯蔵する。
給炭機１０４は、チェーンコンベアを有し、バンカー１０３内に貯蔵されている石炭をチェーンコンベアにより切り出してミル１０５に投入する。
ミル１０５は、給炭機１０４から投入された石炭を粉砕する粉砕機である。ミル１０５の入側の位置における圧力が負圧に保たれるようにすることにより、ミル１０５の内部の圧力は負圧に保たれる。ミル１０５は、例えば、ロールミル１０５ａと粉砕テーブル１０５ｂとを有する。ミル１０５の上部から投入された石炭をロールミル１０５ａと粉砕テーブル１０５ｂとの間に供給する。回転している粉砕テーブル１０５ｂに対してロールミル１０５ａを押し付けながら回転させることにより、石炭は押し潰されて粉砕される。粉砕された石炭は、熱ガス発生装置１０１から供給された排ガスの流れにのって、ミル１０５の上部に供給され、分級機で分級された後、外部に放出される。
給炭機１０４は、チェーンコンベアを有し、バンカー１０３内に貯蔵されている石炭をチェーンコンベアにより切り出してミル１０５に投入する。
ミル１０５は、給炭機１０４から投入された石炭を粉砕する粉砕機である。ミル１０５の入側の位置における圧力が負圧に保たれるようにすることにより、ミル１０５の内部の圧力は負圧に保たれる。ミル１０５は、例えば、ロールミル１０５ａと粉砕テーブル１０５ｂとを有する。ミル１０５の上部から投入された石炭をロールミル１０５ａと粉砕テーブル１０５ｂとの間に供給する。回転している粉砕テーブル１０５ｂに対してロールミル１０５ａを押し付けながら回転させることにより、石炭は押し潰されて粉砕される。粉砕された石炭は、熱ガス発生装置１０１から供給された排ガスの流れにのって、ミル１０５の上部に供給され、分級機で分級された後、外部に放出される。
この際、シールエアファン１０６からミル１０５の内部（粉砕テーブル１０５ｂの軸受部）の隙間にシールエアを供給することにより、その隙間から外部に放出されようとする微粉炭を、熱ガス発生装置１０１から供給された排ガスの流れに押し戻す。ミル１０５の内部の圧力がシールエアの圧力未満になるように、シールエアの流速が定められる。このように、シールエアは、粉砕テーブル１０５ｂの軸受部に微粉炭が進入し、その結果として、粉砕テーブル１０５ｂの軸受部の潤滑不良が起こることと、粉砕テーブル１０５ｂの軸受部から外部に放出されることとを防止するためのものである。
ミル入口温度計１２０は、ミル１０５の入口側（熱ガス発生装置１０１とミル１０５との間）の所定の位置における（配管内の（ガスの））温度であるミル入口温度を測定する。
バグフィルター１０７は、ミル１０５から放出された微粉炭を、濾布を用いて捕集する濾過式の捕集機である。ミル１０５と同様に、バグフィルター１０７の内部の圧力も負圧に保たれている。微粉炭以外の異物がバグフィルター１０７で捕集されることがある。異物除去装置１０８は、この異物を除去するためのものである。このように異物除去装置１０８で異物が除去された後、リザーバタンク１０９に微粉炭が貯蔵される。リザーバタンク１０９に貯蔵された微粉炭は、高炉の羽口から高炉の内部に吹き込まれる（微粉炭吹き込みが行われる）。
バグフィルター１０７は、ミル１０５から放出された微粉炭を、濾布を用いて捕集する濾過式の捕集機である。ミル１０５と同様に、バグフィルター１０７の内部の圧力も負圧に保たれている。微粉炭以外の異物がバグフィルター１０７で捕集されることがある。異物除去装置１０８は、この異物を除去するためのものである。このように異物除去装置１０８で異物が除去された後、リザーバタンク１０９に微粉炭が貯蔵される。リザーバタンク１０９に貯蔵された微粉炭は、高炉の羽口から高炉の内部に吹き込まれる（微粉炭吹き込みが行われる）。
ミル出口温度計１１０は、ミル１０５の出口側（ミル１０５とバグフィルター１０７との間）の所定の位置における（配管内の（微粉炭の））温度であるミル出口温度を測定する。
ベンチュリ管１１１は、バグフィルター１０７を通過した排ガスの流量を測定する。
ダンパー１１２は、バグフィルター１０７を通過した排ガスの流量を調整する。
循環ファン１１３は、ダンパー１１２を通過した排ガスを熱ガス発生装置１０１に循環させることができるように、排ガスを昇圧する。
循環ファン１１３により昇圧された排ガスの一部は、煙突１１４を介して大気中に放出される。放散系圧力調整弁１１５は、このようにして大気中に放出される排ガスの圧力を調整するためのものである。
ベンチュリ管１１１は、バグフィルター１０７を通過した排ガスの流量を測定する。
ダンパー１１２は、バグフィルター１０７を通過した排ガスの流量を調整する。
循環ファン１１３は、ダンパー１１２を通過した排ガスを熱ガス発生装置１０１に循環させることができるように、排ガスを昇圧する。
循環ファン１１３により昇圧された排ガスの一部は、煙突１１４を介して大気中に放出される。放散系圧力調整弁１１５は、このようにして大気中に放出される排ガスの圧力を調整するためのものである。
循環系圧力調整弁１１６は、循環ファン１１３により昇圧された排ガスのうち、煙突１１４を介して大気中に放出されずに熱ガス発生装置１０１に循環させる排ガスの圧力を調整するためのものである。このようにして、熱ガス発生装置１０１で発生した排ガスは、循環ガスとして再び熱ガス発生装置１０１に供給され、熱ガス発生装置１０１、ミル１０５、バグフィルター１０７、ベンチュリ管１１１、ダンパー１１２、循環ファン１１３、循環系圧力調整弁１１６、熱ガス発生装置１０１の経路を循環する。
本実施形態では、大気中の空気（エア）を、希釈エアとして、負圧式・排ガス循環系のＰＣＩプラントに供給する。オリフィス流量計１１７は、この希釈エアの流量を測定する。エア流量調整弁１１８は、負圧式・排ガス循環系のＰＣＩプラントに供給される希釈エアの流量を調整するためのものである。希釈エアファン１１９は、エア流量調整弁１１８で流量が調整された希釈エアを昇圧し、希釈エアを熱ガス発生装置１０１の入側の配管に押し込む。これにより、循環ガスの酸素濃度を調整することができる。
粉砕制御装置２００は、ＰＩＤ制御を行うことにより、ミル出口温度計１１０で測定されたミル出口温度の目標値に対する偏差が０（ゼロ）になるようなバーナー負荷を導出して熱ガス発生装置１０１に出力する。また、粉砕制御装置２００は、ミル入口温度計１２０で測定されたミル入口温度が可及的に上限値を上回らない給炭量を導出して給炭機１０４に出力する。
（粉砕制御装置２００の機能構成）
図２は、粉砕制御装置２００の機能構成の一例を示す図である。粉砕制御装置２００は、例えば、ＣＰＵ、ＲＯＭ、ＲＡＭ、ＨＤＤ、および各種のインターフェースを備える情報処理装置（ＰＣ）や、プログラマブルロジックコントローラ（ＰＬＣ）や、専用のハードウェアを用いることにより実現される。
図２は、粉砕制御装置２００の機能構成の一例を示す図である。粉砕制御装置２００は、例えば、ＣＰＵ、ＲＯＭ、ＲＡＭ、ＨＤＤ、および各種のインターフェースを備える情報処理装置（ＰＣ）や、プログラマブルロジックコントローラ（ＰＬＣ）や、専用のハードウェアを用いることにより実現される。
＜ミル出口温度目標値記憶部２０１＞
ミル出口温度目標値記憶部２０１は、ミル出口温度の目標値［℃］を記憶する。ミル出口温度の目標値は、微粉炭の乾燥度等に応じて、オペレータにより設定される。
＜ミル出口温度測定値取得部２０２＞
ミル出口温度測定値取得部２０２は、ミル出口温度計１１０で測定されたミル出口温度［℃］を、例えば一定周期で取得する。一定周期としては、例えば、１［分］を採用することができる。
ミル出口温度目標値記憶部２０１は、ミル出口温度の目標値［℃］を記憶する。ミル出口温度の目標値は、微粉炭の乾燥度等に応じて、オペレータにより設定される。
＜ミル出口温度測定値取得部２０２＞
ミル出口温度測定値取得部２０２は、ミル出口温度計１１０で測定されたミル出口温度［℃］を、例えば一定周期で取得する。一定周期としては、例えば、１［分］を採用することができる。
＜ミル出口温度偏差導出部２０３＞
ミル出口温度偏差導出部２０３は、ミル出口温度目標値記憶部２０１に記憶されたミル出口温度の目標値から、ミル出口温度測定値取得部２０２で取得されたミル出口温度の測定値を減算して、ミル出口温度の測定値の目標値に対する偏差［℃］を導出する。尚、以下の説明では、ミル出口温度の測定値の目標値に対する偏差を必要に応じてミル出口温度偏差と称する。
ミル出口温度偏差導出部２０３は、ミル出口温度目標値記憶部２０１に記憶されたミル出口温度の目標値から、ミル出口温度測定値取得部２０２で取得されたミル出口温度の測定値を減算して、ミル出口温度の測定値の目標値に対する偏差［℃］を導出する。尚、以下の説明では、ミル出口温度の測定値の目標値に対する偏差を必要に応じてミル出口温度偏差と称する。
＜ＰＩＤ制御部２０４＞
ＰＩＤ制御部２０４は、ミル出口温度偏差導出部２０３で導出されたミル出口温度偏差を入力として、比例動作、積分動作、および微分動作を行い、操作量としてバーナー負荷を導出して熱ガス発生装置１０１に出力することを繰り返して、ミル出口温度の測定値を目標値に近づける制御（すなわちＰＩＤ制御）を行う。
ＰＩＤ制御部２０４は、ミル出口温度偏差導出部２０３で導出されたミル出口温度偏差を入力として、比例動作、積分動作、および微分動作を行い、操作量としてバーナー負荷を導出して熱ガス発生装置１０１に出力することを繰り返して、ミル出口温度の測定値を目標値に近づける制御（すなわちＰＩＤ制御）を行う。
＜ミル出口温度記憶部２０５＞
ミル出口温度記憶部２０５は、ミル出口温度測定値取得部２０２で取得されたミル出口温度の測定値であって、一定時間分のミル出口温度の測定値を記憶する。すなわち、ミル出口温度記憶部２０５は、ミル出口温度測定値取得部２０２により新たなミル出口温度の測定値が取得されると、ミル出口温度記憶部２０５に記憶されているミル出口温度の測定値のうち最古に記憶したミル出口温度の測定値を破棄し、当該新たなミル出口温度の測定値を記憶する。一定時間としては、例えば１０００［分］を採用することができる。
ミル出口温度記憶部２０５は、ミル出口温度測定値取得部２０２で取得されたミル出口温度の測定値であって、一定時間分のミル出口温度の測定値を記憶する。すなわち、ミル出口温度記憶部２０５は、ミル出口温度測定値取得部２０２により新たなミル出口温度の測定値が取得されると、ミル出口温度記憶部２０５に記憶されているミル出口温度の測定値のうち最古に記憶したミル出口温度の測定値を破棄し、当該新たなミル出口温度の測定値を記憶する。一定時間としては、例えば１０００［分］を採用することができる。
＜ミル出口温度ばらつき導出部２０６＞
ミル出口温度ばらつき導出部２０６は、ミル出口温度記憶部２０５に記憶されている一定時間分のミル出口温度の測定値を用いて、ミル出口温度のばらつき［℃］を導出する。本実施形態では、ばらつきとして標準偏差を用いる場合を例に挙げて説明する。尚、ミル出口温度ばらつき導出部２０６は、ミル出口温度の測定値の標準偏差を導出する際に、必ずしもミル出口温度記憶部２０５に記憶されている一定時間分のミル出口温度の測定値の全てを用いる必要はない。
ミル出口温度ばらつき導出部２０６は、ミル出口温度記憶部２０５に記憶されている一定時間分のミル出口温度の測定値を用いて、ミル出口温度のばらつき［℃］を導出する。本実施形態では、ばらつきとして標準偏差を用いる場合を例に挙げて説明する。尚、ミル出口温度ばらつき導出部２０６は、ミル出口温度の測定値の標準偏差を導出する際に、必ずしもミル出口温度記憶部２０５に記憶されている一定時間分のミル出口温度の測定値の全てを用いる必要はない。
＜ミル入口温度ばらつき導出部２０７＞
ミル入口温度ばらつき導出部２０７は、ミル出口温度ばらつき導出部２０６で導出されたミル出口温度の測定値のばらつきに、ミル出口温度をミル入口温度に換算するための係数を乗算した値をミル入口温度のばらつきの推定値［℃］として導出する。
本実施形態では、ＰＩＤ制御部２０４のゲインを用いて、ミル出口温度をミル入口温度に換算するための係数を定める。
ミル入口温度ばらつき導出部２０７は、ミル出口温度ばらつき導出部２０６で導出されたミル出口温度の測定値のばらつきに、ミル出口温度をミル入口温度に換算するための係数を乗算した値をミル入口温度のばらつきの推定値［℃］として導出する。
本実施形態では、ＰＩＤ制御部２０４のゲインを用いて、ミル出口温度をミル入口温度に換算するための係数を定める。
石炭の粉砕が安定して行われている場合、ミル出口温度の目標値は一定であり、ミル出口温度の測定値は、ＰＩＤ制御部２０４による制御によって目標値に近い値になる。一方、ミル入口温度は、バーナー負荷の変動によりばらつきが生じる。このバーナー負荷の変動は、ＰＩＤ制御部２０４によるミル出口温度の制御の結果として現れる。したがって、ＰＩＤ制御部２０４のゲインが大きいほど、バーナー負荷を大きく変動させることになり、その結果、ミル入口温度のばらつきも大きくなる。このことから、本実施形態では、ＰＩＤ制御部２０４のゲインを基準にして、ミル出口温度をミル入口温度に換算するための係数を定める。ミル出口温度をミル入口温度に換算するための係数として、ＰＩＤ制御部２０４のゲインそのものを用いてもよいし、ＰＩＤ制御部２０４のゲインに調整値を加算または減算した値を用いてもよい。また、ＰＩＤ制御部２０４のゲインと積分時間とを用いて、ミル出口温度をミル入口温度に換算するための係数を定めてもよい。
また、前述したように本実施形態では、ばらつきとして標準偏差を用いる。
したがって、本実施形態では、ミル入口温度ばらつき導出部２０７は、以下の（１）式により、ミル入口温度の標準偏差の推定値を導出する。
ミル入口温度の標準偏差の推定値＝ミル出口温度の測定値の標準偏差×係数 ・・・（１）
したがって、本実施形態では、ミル入口温度ばらつき導出部２０７は、以下の（１）式により、ミル入口温度の標準偏差の推定値を導出する。
ミル入口温度の標準偏差の推定値＝ミル出口温度の測定値の標準偏差×係数 ・・・（１）
＜ミル入口温度上限値記憶部２０８＞
ミル入口温度上限値記憶部２０８は、ミル入口温度の上限値［℃］を記憶する。ミル入口温度の上限値としては、配管を繋ぐ伸縮管の燃焼やミル１０５に供給された原料（石炭）の発火などが生じない範囲の温度が、オペレータにより設定される。ミル入口温度の上限値は、例えば、３００［℃］〜４００［℃］の範囲で設定される。
ミル入口温度上限値記憶部２０８は、ミル入口温度の上限値［℃］を記憶する。ミル入口温度の上限値としては、配管を繋ぐ伸縮管の燃焼やミル１０５に供給された原料（石炭）の発火などが生じない範囲の温度が、オペレータにより設定される。ミル入口温度の上限値は、例えば、３００［℃］〜４００［℃］の範囲で設定される。
＜外気温取得部２０９＞
外気温取得部２０９は、外気温の現在値［℃］を取得する。外気温は、大気中の温度を測定する温度計で測定された測定値であっても推定値であってもよい。例えば、希釈エアの温度から所定の温度を増減した値を外気温の推定値として用いることができる。
外気温取得部２０９は、外気温の現在値［℃］を取得する。外気温は、大気中の温度を測定する温度計で測定された測定値であっても推定値であってもよい。例えば、希釈エアの温度から所定の温度を増減した値を外気温の推定値として用いることができる。
＜石炭水分取得部２１０＞
石炭水分取得部２１０は、粉砕中の石炭の（単位質量当たりの）水分量［質量％］の現在値を取得する。以下の説明では、石炭の水分量を必要に応じて石炭水分と称する。石炭水分は、測定値であっても推定値であってもよい。例えば、特開２０１１−１３３４５０号公報に記載のようにして石炭水分の測定値を得ることができる。また、例えば、特開平９−４８３９号公報に記載のようにして石炭水分の推定値を得ることができる。このように石炭水分を得る方法は公知の技術で実現できるので、ここでは、その詳細な説明を省略する。また、この他に、特願２０１４−１８９２５２号明細書に記載のように、燃料ガス（ＢＦＧ）の流量、燃焼エアの流量、希釈エアの流量、煙突１１４から大気中に放出されるガス（放散ガス）の流量、外気温、燃料ガスの温度、燃焼エアの温度、希釈エアの温度、および石炭水分を含む物理量を変数として有する計算式であって、負圧式・排ガス循環系のＰＣＩプラントにおける熱の収支のバランスをとる計算を行うための複数の計算式の計算を行うことにより、石炭水分の推定値を得ることができる。
石炭水分取得部２１０は、粉砕中の石炭の（単位質量当たりの）水分量［質量％］の現在値を取得する。以下の説明では、石炭の水分量を必要に応じて石炭水分と称する。石炭水分は、測定値であっても推定値であってもよい。例えば、特開２０１１−１３３４５０号公報に記載のようにして石炭水分の測定値を得ることができる。また、例えば、特開平９−４８３９号公報に記載のようにして石炭水分の推定値を得ることができる。このように石炭水分を得る方法は公知の技術で実現できるので、ここでは、その詳細な説明を省略する。また、この他に、特願２０１４−１８９２５２号明細書に記載のように、燃料ガス（ＢＦＧ）の流量、燃焼エアの流量、希釈エアの流量、煙突１１４から大気中に放出されるガス（放散ガス）の流量、外気温、燃料ガスの温度、燃焼エアの温度、希釈エアの温度、および石炭水分を含む物理量を変数として有する計算式であって、負圧式・排ガス循環系のＰＣＩプラントにおける熱の収支のバランスをとる計算を行うための複数の計算式の計算を行うことにより、石炭水分の推定値を得ることができる。
＜給炭量取得部２１１＞
給炭量取得部２１１は、給炭量［ｔ／ｈ］の現在値を取得する。給炭量は、給炭機１０４からミル１０５に供給される石炭の単位時間当たり（ここでは１時間当たり）の重量である。
本実施形態では、オペレータが、給炭機１０４に対して給炭量を手動で設定できるものとする。オペレータは、例えば、操業計画に基づく生産量を達成できるかどうかと、設備の損傷が発生しないかどうかとを指標として、給炭量を決定し、決定した給炭量に応じて、例えば、給炭機１０４に備わるチェーンコンベアのモータの回転数を設定する。このようにして設定された回転数でチェーンコンベアのモータが動作することにより、所望の給炭量で石炭がミル１０５に供給される。
給炭量取得部２１１は、給炭量［ｔ／ｈ］の現在値を取得する。給炭量は、給炭機１０４からミル１０５に供給される石炭の単位時間当たり（ここでは１時間当たり）の重量である。
本実施形態では、オペレータが、給炭機１０４に対して給炭量を手動で設定できるものとする。オペレータは、例えば、操業計画に基づく生産量を達成できるかどうかと、設備の損傷が発生しないかどうかとを指標として、給炭量を決定し、決定した給炭量に応じて、例えば、給炭機１０４に備わるチェーンコンベアのモータの回転数を設定する。このようにして設定された回転数でチェーンコンベアのモータが動作することにより、所望の給炭量で石炭がミル１０５に供給される。
＜最大給炭可能量導出部２１２＞
最大給炭可能量導出部２１２は、ミル入口温度ばらつき導出部２０７により導出されたミル入口温度のばらつきと、ミル入口温度上限値記憶部２０８により記憶されたミル入口温度の上限値とに基づいて、ミル入口温度の設定値［℃］を導出する。ミル入口温度の設定値は、ミル入口温度の測定値がミル入口温度の上限値を上回らないにようするためのミル入口温度の目標値である。
最大給炭可能量導出部２１２は、ミル入口温度ばらつき導出部２０７により導出されたミル入口温度のばらつきと、ミル入口温度上限値記憶部２０８により記憶されたミル入口温度の上限値とに基づいて、ミル入口温度の設定値［℃］を導出する。ミル入口温度の設定値は、ミル入口温度の測定値がミル入口温度の上限値を上回らないにようするためのミル入口温度の目標値である。
前述したように本実施形態では、ばらつきとして標準偏差を用いる。また、本実施形態では、ミル入口温度は、定常値に対し正規分布に従ってばらつくものと仮定する。そこで、本実施形態では、ミル入口温度の測定値のうち、約９９［％］が上限値以下になるような目標値をミル入口温度の設定値とする。以上のことから本実施形態では、最大給炭可能量導出部２１２は、以下の（２）式により、ミル入口温度の設定値を導出する。
ミル入口温度の設定値＝ミル入口温度の上限値−３×ミル入口温度の標準偏差の推定値 ・・・（２）
ミル入口温度の設定値＝ミル入口温度の上限値−３×ミル入口温度の標準偏差の推定値 ・・・（２）
次に、最大給炭可能量導出部２１２は、ミル入口温度の設定値と、外気温取得部２０９で取得された外気温と、石炭水分取得部２１０で取得された石炭水分と、を入力として、最大給炭可能量［ｔ／ｈ］を導出する。最大給炭可能量とは、ミル入口温度の測定値がミル入口温度の設定値を上回らない範囲で最大となる定常時の給炭量をいう。
本実施形態では、以下の（３）式で表されるミル入口温度の予測式に基づいて最大可能給炭量の計算式を定める。
ミル入口温度の予測値＝α×給炭量×石炭水分＋β×給炭量×石炭水分×外気温＋γ×給炭量×外気温＋ε ・・・（３）
ミル入口温度の予測値＝α×給炭量×石炭水分＋β×給炭量×石炭水分×外気温＋γ×給炭量×外気温＋ε ・・・（３）
（３）式は、ミル入口温度の予測値［℃］は、給炭量×水分と、給炭量×水分×外気温と、給炭量×外気温に影響していると考えられることに基づいて得られるものである。すなわち、ミル１０５の前後の熱バランスに基づいている。ミル入口温度予測係数（水分潜熱）α［℃／ｔ・ｈ］、ミル入口温度予測係数（水分顕熱）β［ｋｃａｌ／（Ｎｍ3・℃）］、ミル入口温度予測係数（石炭顕熱）γ［℃／％・ｔ・ｈ］は、それぞれ、これらの係数である。また、ミル入口温度予測係数（定数補正）ε［ｋｃａｌ／（Ｎｍ3・℃）］は補正項である。
ミル入口温度予測係数α、β、γ、εは、例えば、過去の操業データ（例えば６か月分の操業データ）を用いて重回帰分析を行うことにより導出される値である。したがって、ミル入口温度予測係数α、β、γ、εは、オペレータにより予め最大給炭可能量導出部２１２に設定される。
（３）式のミル入口温度の予測値をミル入口温度の設定値に、給炭量を最大給炭可能量にそれぞれ置き換えることにより、以下の（４）式が得られる。
最大給炭可能量＝（ミル入口温度の設定値−ε）／｛（α×石炭水分）＋β×石炭水分×外気温＋γ×外気温｝ ・・・（４）
本実施形態では、最大給炭可能量導出部２１２は、（４）式の計算を行う事により、最大給炭可能量を導出する。
最大給炭可能量＝（ミル入口温度の設定値−ε）／｛（α×石炭水分）＋β×石炭水分×外気温＋γ×外気温｝ ・・・（４）
本実施形態では、最大給炭可能量導出部２１２は、（４）式の計算を行う事により、最大給炭可能量を導出する。
＜給炭量設定部２１３＞
給炭量設定部２１３は、最大給炭可能量導出部２１２により導出された最大可能給炭量と、給炭量取得部２１１により取得された給炭量（の現在値）とを比較し、比較した結果に基づいて、給炭量の設定を変更するか否かを判定し、給炭量を手動で変更する場合には給炭量の設定の変更を指示する情報を出力する。
給炭量設定部２１３は、最大給炭可能量導出部２１２により導出された最大可能給炭量と、給炭量取得部２１１により取得された給炭量（の現在値）とを比較し、比較した結果に基づいて、給炭量の設定を変更するか否かを判定し、給炭量を手動で変更する場合には給炭量の設定の変更を指示する情報を出力する。
本実施形態では、給炭量設定部２１３は、最大給炭可能量が給炭量（の現在値）を上回る状態が一定時間継続した場合、（最新の）最大給炭可能量の小数点以下を切り捨てた値を給炭量の変更値として導出する。また、給炭量設定部２１３は、最大給炭可能量が給炭量（の現在値）を下回る状態が一定時間継続した場合にも、（最新の）最大給炭可能量の小数点以下を切り捨てた値を給炭量の変更値として導出する。ここで、本実施形態では、最大給炭可能量が給炭量（の現在値）を上回る状態が継続しているか否かを判定するための前記一定時間を、最大給炭可能量が給炭量（の現在値）を下回る状態が継続しているか否かを判定するための前記一定時間よりも長くする。例えば、前者の一定時間として３００［分］を、後者の一定時間といて１２０［分］をそれぞれ採用することができる。
本実施形態では、これら以外の場合には、給炭量（の現在値）を変更しない。
本実施形態では、これら以外の場合には、給炭量（の現在値）を変更しない。
給炭量設定部２１３は、以上のようにして導出した変更値で石炭が給炭機１０４からミル１０５に供給されるように給炭量の設定の変更を指示する情報を出力する。粉砕制御装置２００が自動的に給炭量を変更する場合、給炭量設定部２１３は、導出した変更値に対応する回転数で給炭機１０４に備わるチェーンコンベアを動作させることを指示する制御信号を、給炭機１０４に備わるチェーンコンベア、または、給炭機１０４に備わるチェーンコンベアの駆動装置に送信する。また、オペレータが手動で給炭量を変更する場合、給炭量設定部２１３は、導出した変更値を示す情報を表示する。オペレータは、この表示に基づいて、変更値に対応する回転数で給炭機１０４に備わるチェーンコンベアを動作させるための操作を行う。
（動作フローチャート）
次に、図３のフローチャートを参照しながら、給炭量を変更するための処理を行う際の粉砕制御装置２００の動作の一例を説明する。
まず、ステップＳ３０１において、ミル出口温度測定値取得部２０２は、ミル出口温度の測定値を取得する。前述したように、ミル出口温度の測定値は、例えば１［分］周期で繰り返し取得される。この場合、図３のフローチャートは、１［分］周期で繰り返し行われる。
次に、図３のフローチャートを参照しながら、給炭量を変更するための処理を行う際の粉砕制御装置２００の動作の一例を説明する。
まず、ステップＳ３０１において、ミル出口温度測定値取得部２０２は、ミル出口温度の測定値を取得する。前述したように、ミル出口温度の測定値は、例えば１［分］周期で繰り返し取得される。この場合、図３のフローチャートは、１［分］周期で繰り返し行われる。
次に、ステップＳ３０２において、ミル出口温度記憶部２０５は、ステップＳ３０１で取得された（最新の）ミル出口温度の測定値を記憶すると共に、既に記憶しているミル出口温度の測定値のうち最古に記憶したミル出口温度の測定値を破棄する。これにより、一定時間分のミル出口温度の測定値が記憶される。
次に、ステップＳ３０３において、ミル出口温度ばらつき導出部２０６は、ミル出口温度記憶部２０５に記憶されている一定時間分のミル出口温度の測定値を用いて、ミル出口温度の測定値の標準偏差を導出する。
次に、ステップＳ３０４において、ミル入口温度ばらつき導出部２０７は、（１）式の計算を行うことにより、ミル入口温度の標準偏差の推定値を導出する。
次に、ステップＳ３０４において、ミル入口温度ばらつき導出部２０７は、（１）式の計算を行うことにより、ミル入口温度の標準偏差の推定値を導出する。
次に、ステップＳ３０５において、最大給炭可能量導出部２１２は、ステップＳ３０４で導出されたミル入口温度の標準偏差の推定値と、ミル入口温度上限値記憶部２０８により記憶されたミル入口温度の上限値とに基づいて、（２）式の計算を行うことにより、ミル入口温度の設定値を導出する。
次に、ステップＳ３０６において、外気温取得部２０９は、外気温の現在値を取得する。
次に、ステップＳ３０７において、石炭水分取得部２１０は、石炭水分の現在値を取得する。
次に、ステップＳ３０８において、給炭量取得部２１１は、給炭量の現在値を取得する。
次に、ステップＳ３０７において、石炭水分取得部２１０は、石炭水分の現在値を取得する。
次に、ステップＳ３０８において、給炭量取得部２１１は、給炭量の現在値を取得する。
次に、ステップＳ３０９において、最大給炭可能量導出部２１２は、ステップＳ３０５で導出されたミル入口温度の設定値と、ステップＳ３０６で取得された外気温と、ステップＳ３０７で取得された石炭水分と、ステップＳ３０８で取得された給炭量とに基づいて、（４）式の計算を行うことにより、最大給炭可能量を導出する。
次に、ステップＳ３１０において、給炭量設定部２１３は、ステップＳ３０９で導出された最大給炭可能量が、ステップＳ３０８で取得された給炭量を上回る状態が一定時間継続したか否かを判定する。この判定の結果、最大給炭可能量が給炭量を上回る状態が一定時間継続した場合には、ステップＳ３１１に進む。ステップＳ３１１に進むと、給炭量設定部２１３は、ステップＳ３０９で導出された最大給炭可能量の小数点以下を切り捨てた値を給炭量の変更値として導出し、導出した変更値で石炭が給炭機１０４からミル１０５に供給されるように給炭量の設定の変更を指示する情報を出力する。そして、後述するステップＳ３１３に進む。
ステップＳ３１０において、最大給炭可能量が給炭量を上回る状態が一定時間継続していないと判定された場合には、ステップＳ３１２に進む。ステップＳ３１２に進むと、給炭量設定部２１３は、ステップＳ３０９で導出された最大給炭可能量が、ステップＳ３０８で取得された給炭量を下回る状態が一定時間継続したか否かを判定する。この判定の結果、最大給炭可能量が給炭量を下回る状態が一定時間継続した場合には、ステップＳ３１１に進む。ステップＳ３１１に進むと、前述したように、給炭量設定部２１３は、ステップＳ３０９で導出された最大給炭可能量の小数点以下を切り捨てた値を給炭量の変更値として導出し、導出した変更値で石炭が給炭機１０４からミル１０５に供給されるように給炭量の設定の変更を指示する情報を出力する。そして、後述するステップＳ３１３に進む。
一方、最大給炭可能量が給炭量を下回る状態が一定時間継続していない場合には、ステップＳ３１１の処理を行わずにステップＳ３１３に進む。ステップＳ３１３に進むと、粉砕制御装置２００は、粉砕を終了するか否かを判定する。この判定は、例えば、負圧式・排ガス循環系のＰＣＩプラントの操業を管理する上位のコンピュータから送信される情報に基づいて行うことができる。
この判定の結果、粉砕を終了する場合には、図３のフローチャートによる処理を終了する。一方、粉砕を終了しない場合には、ステップＳ３０１に戻り、粉砕を終了すると判定するまで、ステップＳ３０１〜Ｓ３１３を繰り返し行う。
この判定の結果、粉砕を終了する場合には、図３のフローチャートによる処理を終了する。一方、粉砕を終了しない場合には、ステップＳ３０１に戻り、粉砕を終了すると判定するまで、ステップＳ３０１〜Ｓ３１３を繰り返し行う。
（実施例）
次に、実施例を説明する。
図４は、ミル入口温度の測定値と時間との関係の一例を示す図である。
図４に示すように、一般にミル入口温度の測定値には、１０［分］〜１５［分］程度の短周期の変動の他に、給炭量の設定値の変更や石炭水分の変化などによる緩やかな定常値の変化が生じる。
次に、実施例を説明する。
図４は、ミル入口温度の測定値と時間との関係の一例を示す図である。
図４に示すように、一般にミル入口温度の測定値には、１０［分］〜１５［分］程度の短周期の変動の他に、給炭量の設定値の変更や石炭水分の変化などによる緩やかな定常値の変化が生じる。
図５は、ミル出口温度の測定値と時間との関係の一例を示す図である。ここでは、ミル出口温度の目標値が一定である間は、図５に示すように、一般にミル出口温度の測定値には、短周期の変動のみが生じる（ミル入口温度の測定値に生じる緩やかな定常値の変化が生じない）。
図６は、ミル入口温度の標準偏差と時間との関係の一例を示す図である。
図６において、グラフ６０１は、図５に示したミル出口温度の測定値から、（１）式の計算を行うことにより得られるミル入口温度の標準偏差の推定値を示すグラフである。グラフ６０１は発明例である。一方、グラフ６０２は、図４に示したミル入口温度の測定値から得られるミル入口温度の標準偏差を示すグラフである。グラフ６０２は比較例である。
図６において、グラフ６０１は、図５に示したミル出口温度の測定値から、（１）式の計算を行うことにより得られるミル入口温度の標準偏差の推定値を示すグラフである。グラフ６０１は発明例である。一方、グラフ６０２は、図４に示したミル入口温度の測定値から得られるミル入口温度の標準偏差を示すグラフである。グラフ６０２は比較例である。
尚、本実施例では、何れの場合も（１）式の係数として同じ値（実際の操業で想定される値）を使用した。また、何れの場合も、直近の１日分の実績（１４４０サンプル）から標準偏差を計算した。
図４に示したようにミル入口温度の測定値には緩やかな定常値の変化がある。したがって、グラフ６０２に示すように、ミル入口温度の測定値からミル入口温度の標準偏差を導出すると、ミル入口温度の測定値に緩やかな定常値の変化があることにより、ミル入口温度の標準偏差（ばらつき量）が大きく見積もられる。これに対し、グラフ６０１に示すように、図５に示したミル出口温度の測定値からミル入口温度の標準偏差を導出すると、グラフ６０２に比べ、ミル入口温度の標準偏差（ばらつき量）が過大に見積まれることはない。前述したように、図５に示したミル出口温度の測定値には、短周期の変動のみが生じる（ミル入口温度の測定値に生じる緩やかな定常値の変化が生じない）からである。
本実施例では、以上のようなミル入口温度の測定値、ミル出口温度の測定値、およびミル入口温度の標準偏差を用いて、コンピュータシミュレーションを行った。
＜実施例１＞
実施例１では、以下の３つの場合について、どの程度の割合でミル入口温度が、その上限値を上回るのかを比較した。
（Ａ） ミル出口温度の測定値に基づいてミル入口温度の設定値を導出した場合：発明例
（Ｂ） ミル入口温度の測定値に基づいてミル入口温度の設定値を導出した場合：比較例１
（Ｃ） ミル入口温度の設定値を固定値とした場合：比較例２
図４に示したようにミル入口温度の測定値には緩やかな定常値の変化がある。したがって、グラフ６０２に示すように、ミル入口温度の測定値からミル入口温度の標準偏差を導出すると、ミル入口温度の測定値に緩やかな定常値の変化があることにより、ミル入口温度の標準偏差（ばらつき量）が大きく見積もられる。これに対し、グラフ６０１に示すように、図５に示したミル出口温度の測定値からミル入口温度の標準偏差を導出すると、グラフ６０２に比べ、ミル入口温度の標準偏差（ばらつき量）が過大に見積まれることはない。前述したように、図５に示したミル出口温度の測定値には、短周期の変動のみが生じる（ミル入口温度の測定値に生じる緩やかな定常値の変化が生じない）からである。
本実施例では、以上のようなミル入口温度の測定値、ミル出口温度の測定値、およびミル入口温度の標準偏差を用いて、コンピュータシミュレーションを行った。
＜実施例１＞
実施例１では、以下の３つの場合について、どの程度の割合でミル入口温度が、その上限値を上回るのかを比較した。
（Ａ） ミル出口温度の測定値に基づいてミル入口温度の設定値を導出した場合：発明例
（Ｂ） ミル入口温度の測定値に基づいてミル入口温度の設定値を導出した場合：比較例１
（Ｃ） ミル入口温度の設定値を固定値とした場合：比較例２
図７は、ミル入口温度の設定値と時間との関係の一例を示す図である。
図７において、グラフ７０１は、図６に示したグラフ６０１のミル入口温度の標準偏差の推定値から（２）式の計算を行うことにより得られるミル入口温度の設定値を示すグラフである。グラフ７０１は発明例（前記（Ａ）の場合）である。一方、グラフ７０２は、図６に示したグラフ６０２のミル入口温度の測定値の標準偏差を（２）式の右辺のミル入口温度の標準偏差の推定値として与えることにより得られるミル入口温度の設定値を示すグラフである。グラフ７０２は比較例１（前記（Ｂ）の場合）である。本実施例では、何れの場合も（２）式のミル入口温度の上限値を３５０［℃］とした。
図７において、グラフ７０１は、図６に示したグラフ６０１のミル入口温度の標準偏差の推定値から（２）式の計算を行うことにより得られるミル入口温度の設定値を示すグラフである。グラフ７０１は発明例（前記（Ａ）の場合）である。一方、グラフ７０２は、図６に示したグラフ６０２のミル入口温度の測定値の標準偏差を（２）式の右辺のミル入口温度の標準偏差の推定値として与えることにより得られるミル入口温度の設定値を示すグラフである。グラフ７０２は比較例１（前記（Ｂ）の場合）である。本実施例では、何れの場合も（２）式のミル入口温度の上限値を３５０［℃］とした。
グラフ７０２に示すように、ミル入口温度の測定値の標準偏差からミル入口温度の設定値を導出すると、ミル入口温度の測定値に緩やかな定常値の変化があることにより（図４を参照）、ミル入口温度の設定値が小さく見積もられる。これに対し、グラフ７０１に示すように、ミル入口温度の標準偏差の推定値からミル入口温度の設定値を導出すると、グラフ７０２に比べ、ミル入口温度の設定値が過小に見積まれることはない。
図８は、実施例１における給炭量と時間との関係の一例を示す図である。
図８において、グラフ８０１は、図７に示したグラフ７０１のミル入口温度の設定値から（４）式の計算を行うことにより最大給炭可能量を導出し、導出した最大給炭可能量からステップＳ３１０〜Ｓ３１２の処理を行うことにより設定される給炭量を示すグラフである。グラフ８０１は発明例（前記（Ａ）の場合）である。グラフ８０２は、図７に示したグラフ７０２のミル入口温度の設定値から（４）式の計算を行うことにより最大給炭可能量を導出し、導出した最大給炭可能量からステップＳ３１０〜Ｓ３１２の処理を行うことにより設定される給炭量を示すグラフである。グラフ８０２は比較例１（前記（Ｂ）の場合）である。グラフ８０３は、ミル入口温度の設定値を固定値として（４）式の計算を行うことにより最大給炭可能量を導出し、導出した最大給炭可能量からステップＳ３１０〜Ｓ３１２の処理を行うことにより設定される給炭量を示すグラフである。グラフ８０３は比較例２（前記（Ｃ）の場合）である。比較例２（前記（Ｃ））では、前記固定値として、グラフ８０１に示す給炭量の平均値と同等の給炭量が得られるような値を用いた。
また、何れの場合でも、ステップＳ３１０における一定時間を３００［分］とし、ステップＳ３１２における一定時間を１２０［分］とした。
図８において、グラフ８０１は、図７に示したグラフ７０１のミル入口温度の設定値から（４）式の計算を行うことにより最大給炭可能量を導出し、導出した最大給炭可能量からステップＳ３１０〜Ｓ３１２の処理を行うことにより設定される給炭量を示すグラフである。グラフ８０１は発明例（前記（Ａ）の場合）である。グラフ８０２は、図７に示したグラフ７０２のミル入口温度の設定値から（４）式の計算を行うことにより最大給炭可能量を導出し、導出した最大給炭可能量からステップＳ３１０〜Ｓ３１２の処理を行うことにより設定される給炭量を示すグラフである。グラフ８０２は比較例１（前記（Ｂ）の場合）である。グラフ８０３は、ミル入口温度の設定値を固定値として（４）式の計算を行うことにより最大給炭可能量を導出し、導出した最大給炭可能量からステップＳ３１０〜Ｓ３１２の処理を行うことにより設定される給炭量を示すグラフである。グラフ８０３は比較例２（前記（Ｃ）の場合）である。比較例２（前記（Ｃ））では、前記固定値として、グラフ８０１に示す給炭量の平均値と同等の給炭量が得られるような値を用いた。
また、何れの場合でも、ステップＳ３１０における一定時間を３００［分］とし、ステップＳ３１２における一定時間を１２０［分］とした。
図９は、発明例におけるミル入口温度と時間との関係を示す図である。
図９に示すグラフ９０１は、ミル入口温度のシミュレーション結果を示すグラフである。このときのミル入口温度は、グラフ８０１に示す給炭量から、（３）式の計算式に基づいて算出した。グラフ９０１は発明例である。
図１０は、比較例１におけるミル入口温度と時間との関係を示す図である。
図１０に示すグラフ１００１は、ミル入口温度のシミュレーション結果を示すグラフである。このときのミル入口温度は、グラフ８０２に示す給炭量から、（３）式の計算式に基づいて算出した。グラフ１００１は比較例１である。
図１１は、比較例２におけるミル入口温度と時間との関係を示す図である。
図１１に示すグラフ１１０１は、ミル入口温度のシミュレーション結果を示すグラフである。このときのミル入口温度は、グラフ８０３に示す給炭量から、（３）式の計算式に基づいて算出した。グラフ１１０１は比較例２である。
図９に示すグラフ９０１は、ミル入口温度のシミュレーション結果を示すグラフである。このときのミル入口温度は、グラフ８０１に示す給炭量から、（３）式の計算式に基づいて算出した。グラフ９０１は発明例である。
図１０は、比較例１におけるミル入口温度と時間との関係を示す図である。
図１０に示すグラフ１００１は、ミル入口温度のシミュレーション結果を示すグラフである。このときのミル入口温度は、グラフ８０２に示す給炭量から、（３）式の計算式に基づいて算出した。グラフ１００１は比較例１である。
図１１は、比較例２におけるミル入口温度と時間との関係を示す図である。
図１１に示すグラフ１１０１は、ミル入口温度のシミュレーション結果を示すグラフである。このときのミル入口温度は、グラフ８０３に示す給炭量から、（３）式の計算式に基づいて算出した。グラフ１１０１は比較例２である。
図１２は、ミル入口温度突破頻度を示す図である。
ミル入口温度突破頻度は、図９、図１０、図１１に示すグラフ９０１、１００１、１１０１の値のうち、ミル入口温度の上限値を上回る値の割合［％］である。図１２に示すようにミル入口温度突破頻度は、グラフ９０１、１００１、１１０１のそれぞれについて個別に導出される。
ミル入口温度突破頻度は、図９、図１０、図１１に示すグラフ９０１、１００１、１１０１の値のうち、ミル入口温度の上限値を上回る値の割合［％］である。図１２に示すようにミル入口温度突破頻度は、グラフ９０１、１００１、１１０１のそれぞれについて個別に導出される。
図１２に示すように、高い平均給炭量を確保する場合、ミル入口温度の設定値を固定値とするよりも（比較例２よりも）、ミル入口温度の標準偏差の推定値からミル入口温度の設定値を導出した方が（発明例の方が）、設備保全性能が高く、ミル入口温度突破頻度は、ミル入口温度の設定値を固定値とする場合の１割程度になる。
＜実施例２＞
実施例２では、以下の３つの場合について、平均給炭量を比較した。
（Ａ） ミル出口温度の測定値に基づいてミル入口温度の設定値を導出した場合：発明例
（Ｂ） ミル入口温度の測定値に基づいてミル入口温度の設定値を導出した場合：比較例１
（Ｄ） ミル入口温度の設定値を固定値とした場合：比較例３
実施例２では、以下の３つの場合について、平均給炭量を比較した。
（Ａ） ミル出口温度の測定値に基づいてミル入口温度の設定値を導出した場合：発明例
（Ｂ） ミル入口温度の測定値に基づいてミル入口温度の設定値を導出した場合：比較例１
（Ｄ） ミル入口温度の設定値を固定値とした場合：比較例３
図１３は、実施例２における給炭量と時間との関係の一例を示す図である。
図１３において、グラフ１３０１、１３０２は、それぞれ図８に示したグラフ８０１、８０２と同じである。
一方、グラフ１３０３は、ミル入口温度の設定値を固定値とした場合のグラフである。グラフ１３０３は比較例３（前記（Ｄ）の場合）である。
図１３において、グラフ１３０１、１３０２は、それぞれ図８に示したグラフ８０１、８０２と同じである。
一方、グラフ１３０３は、ミル入口温度の設定値を固定値とした場合のグラフである。グラフ１３０３は比較例３（前記（Ｄ）の場合）である。
ここで、比較例３（前記（Ｄ））では、前記固定値として、図１２に示したミル入口温度突破頻度が発明例と同じ値（０．０４［％］）になるような値を用いた。
図１４は、比較例３におけるミル入口温度と時間との関係を示す図である。
図１４に示すグラフ１４０１は、グラフ１３０３に示す給炭量から（３）式の計算を行うことにより得られるミル入口温度の予測値のグラフである。グラフ１４０１は比較例３である。尚、発明例および比較例１におけるミル入口温度と時間との関係は、それぞれ、図９に示したグラフ９０１、図１０に示したグラフ１００１と同じである。
図１４は、比較例３におけるミル入口温度と時間との関係を示す図である。
図１４に示すグラフ１４０１は、グラフ１３０３に示す給炭量から（３）式の計算を行うことにより得られるミル入口温度の予測値のグラフである。グラフ１４０１は比較例３である。尚、発明例および比較例１におけるミル入口温度と時間との関係は、それぞれ、図９に示したグラフ９０１、図１０に示したグラフ１００１と同じである。
図１５は、平均給炭量を示す図である。
図１５に示すように、ミル入口温度の標準偏差の推定値からミル入口温度の設定値を導出すると（すなわち発明例では）、ミル入口温度の設定値を固定値とする場合（すなわち比較例３の場合）、ミル入口温度の測定値の標準偏差からミル入口温度の設定値を導出する場合（すなわち比較例１の場合）に対し、それぞれ５［％］、７［％］高い生産量を確保することができる。
以上の実施例１、２より、発明例では、比較例と比較して、設備保全性能の向上と平均給炭量の上昇との両立を図ることができることが分かる。
図１５に示すように、ミル入口温度の標準偏差の推定値からミル入口温度の設定値を導出すると（すなわち発明例では）、ミル入口温度の設定値を固定値とする場合（すなわち比較例３の場合）、ミル入口温度の測定値の標準偏差からミル入口温度の設定値を導出する場合（すなわち比較例１の場合）に対し、それぞれ５［％］、７［％］高い生産量を確保することができる。
以上の実施例１、２より、発明例では、比較例と比較して、設備保全性能の向上と平均給炭量の上昇との両立を図ることができることが分かる。
ミル入口温度の測定値の標準偏差を導出する場合（すなわち比較例１の場合）、図４に示したようなミル入口温度の測定値の定常値の変化が、ミル入口温度の標準偏差（ばらつき量）に反映される。このため、図７に示したように、無用にミル入口温度を下げた運転になる。その結果として、図１５に示すように、平均給炭量が小さくなる。一方、ミル出口温度の測定値からミル入口温度の標準偏差を推定する場合（すなわち発明例の場合）には、このような定常値の変化の影響を受けることなくミル入口温度の標準偏差（ばらつき量）を評価できる。このため、比較例１に比べ発明例では、ミル入口温度を上限値以下に維持することと給炭量を大きくすることとの両立を図ることができる。
（まとめ）
以上のように本実施形態では、一定時間分のミル出口温度の測定値の標準偏差にＰＩＤ制御のゲインを乗算してミル入口温度の標準偏差の推定値を導出し、ミル入口温度の上限値とミル入口温度の標準偏差の推定値とに基づいてミル入口温度の設定値を導出する。そして、このミル入口温度の設定値と、外気温と、石炭水分とに基づいて、最大給炭可能量を導出し、給炭量が最大給炭可能量を上回るまたは下回る状態が一定時間継続した場合に、給炭量が最大給炭可能量に近づくように給炭量の設定値を変更する。このように、目標値になるようにフィードバック制御されるミル出口温度の測定値を用いて、ミル入口温度の標準偏差の推定値を導出するので、給炭量の設定値や石炭水分の変化等の影響を大きく受けることなく、ミル入口温度の標準偏差を導出することができる。したがって、ミル入口温度のばらつきを正確に予測し、ミル入口温度が上限値を可及的に上回らないようにすることと、原料の供給量を可及的に多くすることを両立させることができる。
以上のように本実施形態では、一定時間分のミル出口温度の測定値の標準偏差にＰＩＤ制御のゲインを乗算してミル入口温度の標準偏差の推定値を導出し、ミル入口温度の上限値とミル入口温度の標準偏差の推定値とに基づいてミル入口温度の設定値を導出する。そして、このミル入口温度の設定値と、外気温と、石炭水分とに基づいて、最大給炭可能量を導出し、給炭量が最大給炭可能量を上回るまたは下回る状態が一定時間継続した場合に、給炭量が最大給炭可能量に近づくように給炭量の設定値を変更する。このように、目標値になるようにフィードバック制御されるミル出口温度の測定値を用いて、ミル入口温度の標準偏差の推定値を導出するので、給炭量の設定値や石炭水分の変化等の影響を大きく受けることなく、ミル入口温度の標準偏差を導出することができる。したがって、ミル入口温度のばらつきを正確に予測し、ミル入口温度が上限値を可及的に上回らないようにすることと、原料の供給量を可及的に多くすることを両立させることができる。
（変形例）
＜変形例１＞
本実施形態では、負圧式・排ガス循環系のＰＣＩプラントに粉砕制御装置２００を適用する場合を例に挙げて説明した。しかしながら、負圧式・排ガス循環系のＰＣＩプラント以外の負圧式・排ガス循環系の粉砕プラントにも粉砕制御装置２００を適用することができる。例えば、セメントを製造するための負圧式・排ガス循環系の粉砕プラントにも粉砕制御装置２００を適用することができる。
＜変形例１＞
本実施形態では、負圧式・排ガス循環系のＰＣＩプラントに粉砕制御装置２００を適用する場合を例に挙げて説明した。しかしながら、負圧式・排ガス循環系のＰＣＩプラント以外の負圧式・排ガス循環系の粉砕プラントにも粉砕制御装置２００を適用することができる。例えば、セメントを製造するための負圧式・排ガス循環系の粉砕プラントにも粉砕制御装置２００を適用することができる。
＜変形例２＞
本実施形態では、粉砕プラントが負圧式・排ガス循環系の粉砕プラントである場合を例に挙げて説明した。しかしながら、粉砕プラントは１パス方式の粉砕プラントであってもよい。
図１６は、１パス方式のＰＣＩプラントの構成の一例を示す図である。図１６において、図１に示した負圧式・排ガス循環系のＰＣＩプラントの構成と同じ部分については、図１に付した符号と同一の符号を付す。また、図１６でも、図１と同様に、各構成要素を繋ぐ実線は配管を示し、破線は信号の伝達経路を示す。また、矢印線は、配管内のガスや石炭の進行方向を示す。
本実施形態では、粉砕プラントが負圧式・排ガス循環系の粉砕プラントである場合を例に挙げて説明した。しかしながら、粉砕プラントは１パス方式の粉砕プラントであってもよい。
図１６は、１パス方式のＰＣＩプラントの構成の一例を示す図である。図１６において、図１に示した負圧式・排ガス循環系のＰＣＩプラントの構成と同じ部分については、図１に付した符号と同一の符号を付す。また、図１６でも、図１と同様に、各構成要素を繋ぐ実線は配管を示し、破線は信号の伝達経路を示す。また、矢印線は、配管内のガスや石炭の進行方向を示す。
図１６において、熱ガス発生装置（ＨＧＧ）１０１は、バーナーの空燃比を制御し、排ガス（熱風）を発生させる。燃焼エアファン１０２は、燃焼エアを熱ガス発生装置１０１に送り込むために昇圧するファンである。
ダンパー１２１は、熱ガス発生装置１０１から排出される排ガスの流量を調整する。
ミル入口温度計１２０は、ミル１０５の入口側（熱ガス発生装置１０１とダンパー１２１との間）の所定の位置における（配管内の（ガスの））温度であるミル入口温度を測定する。
ダンパー１２１は、熱ガス発生装置１０１から排出される排ガスの流量を調整する。
ミル入口温度計１２０は、ミル１０５の入口側（熱ガス発生装置１０１とダンパー１２１との間）の所定の位置における（配管内の（ガスの））温度であるミル入口温度を測定する。
供給ファン１２２は、ダンパー１２１を通過した排ガスをミル１０５に供給させるために排ガスを昇圧するファンである。
バンカー１０３は、原料である石炭を貯蔵する。
給炭機１０４は、チェーンコンベアを有し、バンカー１０３内に貯蔵されている石炭をチェーンコンベアにより切り出してミル１０５に投入する。
バンカー１０３は、原料である石炭を貯蔵する。
給炭機１０４は、チェーンコンベアを有し、バンカー１０３内に貯蔵されている石炭をチェーンコンベアにより切り出してミル１０５に投入する。
ミル１０５は、給炭機１０４から投入された石炭を粉砕する粉砕機である。ミル１０５は、例えば、ロールミル１０５ａと粉砕テーブル１０５ｂとを有する。
シールエアファン１０６は、ミル１０５の内部（粉砕テーブル１０５ｂの軸受部）の隙間にシールエアを供給する。
シールエアファン１０６は、ミル１０５の内部（粉砕テーブル１０５ｂの軸受部）の隙間にシールエアを供給する。
バグフィルター１０７は、ミル１０５から放出された微粉炭を、濾布を用いて捕集する濾過式の捕集機である。
異物除去装置１０８は、この異物を除去するためのものである。このように異物除去装置１０８で異物が除去された後、リザーバタンク１０９に微粉炭が貯蔵される。リザーバタンク１０９に貯蔵された微粉炭は、高炉の羽口から高炉の内部に吹き込まれる（微粉炭吹き込みが行われる）。
異物除去装置１０８は、この異物を除去するためのものである。このように異物除去装置１０８で異物が除去された後、リザーバタンク１０９に微粉炭が貯蔵される。リザーバタンク１０９に貯蔵された微粉炭は、高炉の羽口から高炉の内部に吹き込まれる（微粉炭吹き込みが行われる）。
負圧式・排ガス循環系のＰＣＩプラントでは、ミル１０５とバグフィルター１０７の内部の圧力は負圧に保たれているが、１パス方式のＰＣＩプラントでは、ミル１０５とバグフィルター１０７の内部の圧力は大気圧である。
バグフィルター１０７を通過した排ガスは、煙突１２３を介して大気中に放出される。
尚、１パス方式のＰＣＩプラントにおいても粉砕制御装置２００の構成および処理は本実施形態で説明したものと同じもので実現できる。
バグフィルター１０７を通過した排ガスは、煙突１２３を介して大気中に放出される。
尚、１パス方式のＰＣＩプラントにおいても粉砕制御装置２００の構成および処理は本実施形態で説明したものと同じもので実現できる。
＜変形例３＞
本実施形態では、ミル出口温度制御と給炭量制御とを１つの装置（粉砕制御装置２００）で行う場合を例に挙げて説明した。しかしながら、必ずしもこのようにする必要はない。例えば、ミル出口温度制御と給炭量制御とを別々の装置で行うようにしてもよい。このようにする場合、例えば、粉砕制御装置２００が、ミル出口温度制御を行う機能（ミル出口温度目標値記憶部２０１、ミル出口温度測定値取得部２０２、ミル出口温度偏差導出部２０３、およびＰＩＤ制御部２０４の機能）を有さず、この機能を粉砕制御装置２００と通信可能な別の装置で実現してもよい。
本実施形態では、ミル出口温度制御と給炭量制御とを１つの装置（粉砕制御装置２００）で行う場合を例に挙げて説明した。しかしながら、必ずしもこのようにする必要はない。例えば、ミル出口温度制御と給炭量制御とを別々の装置で行うようにしてもよい。このようにする場合、例えば、粉砕制御装置２００が、ミル出口温度制御を行う機能（ミル出口温度目標値記憶部２０１、ミル出口温度測定値取得部２０２、ミル出口温度偏差導出部２０３、およびＰＩＤ制御部２０４の機能）を有さず、この機能を粉砕制御装置２００と通信可能な別の装置で実現してもよい。
＜変形例４＞
本実施形態では、ミル出口温度のばらつきとして標準偏差を導出する場合を例に挙げて説明した。しかしながら、一般的な統計的手法で利用されるばらつきを表す指標を用いてれば、ミル出口温度のばらつきは標準偏差に限定されない。この場合、ミル入口温度のばらつきの推定値も、ミル出口温度のばらつきとして選択されたばらつきを表す指標と同じ指標の推定値になる。
本実施形態では、ミル出口温度のばらつきとして標準偏差を導出する場合を例に挙げて説明した。しかしながら、一般的な統計的手法で利用されるばらつきを表す指標を用いてれば、ミル出口温度のばらつきは標準偏差に限定されない。この場合、ミル入口温度のばらつきの推定値も、ミル出口温度のばらつきとして選択されたばらつきを表す指標と同じ指標の推定値になる。
＜変形例５＞
本実施形態では、ミル出口温度をミル入口温度に換算するための係数が、ＰＩＤ制御のゲインを基準にして設定される係数である場合を例に挙げて説明した。しかしながら、ミル出口温度をミル入口温度に換算するための係数は、ＰＩＤ制御のゲインを基準にして設定される係数に限定されない。例えば、ミル出口温度制御をＰＩＤ制御以外のフィードバック制御で行う場合には、ミル出口温度をミル入口温度に換算するための係数として、当該フィードバック制御を行う制御器における制御ゲインを基準して設定される係数を用いてもよい。また、ミル出口温度とミル入口温度との関係を事前に調査し、調査した結果から係数を設定してもよい。
本実施形態では、ミル出口温度をミル入口温度に換算するための係数が、ＰＩＤ制御のゲインを基準にして設定される係数である場合を例に挙げて説明した。しかしながら、ミル出口温度をミル入口温度に換算するための係数は、ＰＩＤ制御のゲインを基準にして設定される係数に限定されない。例えば、ミル出口温度制御をＰＩＤ制御以外のフィードバック制御で行う場合には、ミル出口温度をミル入口温度に換算するための係数として、当該フィードバック制御を行う制御器における制御ゲインを基準して設定される係数を用いてもよい。また、ミル出口温度とミル入口温度との関係を事前に調査し、調査した結果から係数を設定してもよい。
＜変形例６＞
本実施形態では、ミル入口温度の上限値から、ミル入口温度の標準偏差の推定値に「３」を乗算した値（３σ）を減算した値をミル入口温度の設定値として導出する場合を例に挙げて説明した。しかしながら、ミル入口温度の設定値は、必ずしもこのようにして導出する必要はない。例えば、前述した「３」の代わりに「２」を用いてもよい。また、ミル入口温度は、定常値に対し正規分布以外の分布に従ってばらつくものと仮定してミル入口温度の設定値を導出してもよい。すなわち、ミル入口温度の上限値と、ミル入口温度のばらつきの推定値とに基づいて、ミル入口温度の測定値の所定の割合が上限値以下になるようにするためのミル入口温度の目標値をミル入口温度の設定値として導出していればよい。
本実施形態では、ミル入口温度の上限値から、ミル入口温度の標準偏差の推定値に「３」を乗算した値（３σ）を減算した値をミル入口温度の設定値として導出する場合を例に挙げて説明した。しかしながら、ミル入口温度の設定値は、必ずしもこのようにして導出する必要はない。例えば、前述した「３」の代わりに「２」を用いてもよい。また、ミル入口温度は、定常値に対し正規分布以外の分布に従ってばらつくものと仮定してミル入口温度の設定値を導出してもよい。すなわち、ミル入口温度の上限値と、ミル入口温度のばらつきの推定値とに基づいて、ミル入口温度の測定値の所定の割合が上限値以下になるようにするためのミル入口温度の目標値をミル入口温度の設定値として導出していればよい。
＜変形例７＞
本実施形態では、ミル出口温度の測定値の一定時間分（例えば１０００［分］）のばらつき（標準偏差）を導出する場合を例に挙げて説明した。しかしながら、必ずしもこのようにする必要はない。例えば、ミル出口温度偏差の一定時間分のばらつき（標準偏差）を導出してもよい。このようにすれば、例えば、ミル出口温度の目標値に変更があっても、当該変更前後におけるミル出口温度のばらつきをそのまま比較することができる。
本実施形態では、ミル出口温度の測定値の一定時間分（例えば１０００［分］）のばらつき（標準偏差）を導出する場合を例に挙げて説明した。しかしながら、必ずしもこのようにする必要はない。例えば、ミル出口温度偏差の一定時間分のばらつき（標準偏差）を導出してもよい。このようにすれば、例えば、ミル出口温度の目標値に変更があっても、当該変更前後におけるミル出口温度のばらつきをそのまま比較することができる。
＜変形例８＞
本実施形態では、給炭量取得部２１１は、オペレータが手動で設定した給炭量を取得する場合を例に挙げて説明した。しかしながら、必ずしもこのようにする必要はない。例えば、給炭量取得部２１１は、上位のコンピュータから給炭量を受信してもよい。
本実施形態では、給炭量取得部２１１は、オペレータが手動で設定した給炭量を取得する場合を例に挙げて説明した。しかしながら、必ずしもこのようにする必要はない。例えば、給炭量取得部２１１は、上位のコンピュータから給炭量を受信してもよい。
＜変形例９＞
本実施形態では、（１）式、（２）式でそれぞれ係数を乗算する場合を例に挙げて説明した。しかしながら、（１）式および（２）式の係数を乗算した値を１つの係数として用いてミル入口温度の標準偏差の推定値を導出してもよい。例えば、以下の（５）式により、ミル入口温度の標準偏差の推定値を導出してもよい。
ミル入口温度の標準偏差の推定値＝ミル入口温度の上限値−係数×ミル出口温度の測定値の標準偏差 ・・・（５）
（５）式における係数は、例えば、（１）式に示す係数と（２）式に示す「３」とを乗算した値である。
このようにした場合、ミル入口温度ばらつき導出部２０７は不要になる。また、（２）式の右辺の第２項の値（３×ミル入口温度の標準偏差の推定値）をミル入口温度ばらつき導出部２０７で導出してもよい。
本実施形態では、（１）式、（２）式でそれぞれ係数を乗算する場合を例に挙げて説明した。しかしながら、（１）式および（２）式の係数を乗算した値を１つの係数として用いてミル入口温度の標準偏差の推定値を導出してもよい。例えば、以下の（５）式により、ミル入口温度の標準偏差の推定値を導出してもよい。
ミル入口温度の標準偏差の推定値＝ミル入口温度の上限値−係数×ミル出口温度の測定値の標準偏差 ・・・（５）
（５）式における係数は、例えば、（１）式に示す係数と（２）式に示す「３」とを乗算した値である。
このようにした場合、ミル入口温度ばらつき導出部２０７は不要になる。また、（２）式の右辺の第２項の値（３×ミル入口温度の標準偏差の推定値）をミル入口温度ばらつき導出部２０７で導出してもよい。
＜変形例１０＞
本実施形態では、給炭量設定部２１３は、最大給炭可能量が給炭量（の現在値）を上回る状態または下回る状態が一定時間継続した場合に、（最新の）最大給炭可能量の小数点以下を切り捨てた値を給炭量の変更値として導出する場合を例に挙げて説明した。しかしながら、最大給炭可能量に基づいて給炭量の変更値を導出するようにしていれば、給炭量の変更値を導出する方法は、このような方法に限定されない。例えば、給炭量が最大給炭可能量に近づくように給炭量の設定値を変更する手法は、前述した手法以外に種々の手法がある。また、給炭量が最大給炭可能量に近づくようにせず、給炭量が最大給炭可能量よりも所定の値だけ下回る値に近づくようにしてもよい。
また、ミル入口温度の設定値を表示し、オペレータが、ミル入口温度の設定値に基づいて給炭量を手動で設定してもよい。
本実施形態では、給炭量設定部２１３は、最大給炭可能量が給炭量（の現在値）を上回る状態または下回る状態が一定時間継続した場合に、（最新の）最大給炭可能量の小数点以下を切り捨てた値を給炭量の変更値として導出する場合を例に挙げて説明した。しかしながら、最大給炭可能量に基づいて給炭量の変更値を導出するようにしていれば、給炭量の変更値を導出する方法は、このような方法に限定されない。例えば、給炭量が最大給炭可能量に近づくように給炭量の設定値を変更する手法は、前述した手法以外に種々の手法がある。また、給炭量が最大給炭可能量に近づくようにせず、給炭量が最大給炭可能量よりも所定の値だけ下回る値に近づくようにしてもよい。
また、ミル入口温度の設定値を表示し、オペレータが、ミル入口温度の設定値に基づいて給炭量を手動で設定してもよい。
＜その他の変形例＞
尚、以上説明した本発明の実施形態は、コンピュータがプログラムを実行することによって実現することができる。また、前記プログラムを記録したコンピュータ読み取り可能な記録媒体及び前記プログラム等のコンピュータプログラムプロダクトも本発明の実施形態として適用することができる。記録媒体としては、例えば、フレキシブルディスク、ハードディスク、光ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、磁気テープ、不揮発性のメモリカード、ＲＯＭ等を用いることができる。
また、以上説明した本発明の実施形態は、何れも本発明を実施するにあたっての具体化の例を示したものに過ぎず、これらによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその技術思想、またはその主要な特徴から逸脱することなく、様々な形で実施することができる。
尚、以上説明した本発明の実施形態は、コンピュータがプログラムを実行することによって実現することができる。また、前記プログラムを記録したコンピュータ読み取り可能な記録媒体及び前記プログラム等のコンピュータプログラムプロダクトも本発明の実施形態として適用することができる。記録媒体としては、例えば、フレキシブルディスク、ハードディスク、光ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、磁気テープ、不揮発性のメモリカード、ＲＯＭ等を用いることができる。
また、以上説明した本発明の実施形態は、何れも本発明を実施するにあたっての具体化の例を示したものに過ぎず、これらによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその技術思想、またはその主要な特徴から逸脱することなく、様々な形で実施することができる。
（請求項との関係）
熱風発生装置は、例えば、熱ガス発生装置１０１により実現される。
粉砕機は、例えば、ミル１０５により実現される。
捕集機は、例えば、バグフィルター１０７により実現される。
処理装置は、例えば、粉砕制御装置２００により実現される。
ミル出口温度ばらつき導出手段は、例えば、ミル出口温度ばらつき導出部２０６（ステップＳ３０３の処理）により実現される。
ミル入口温度設定値導出手段は、例えば、最大給炭可能量導出部２１２（ステップＳ３０５の処理）により実現される。
ミル入口温度ばらつき導出手段は、例えば、ミル入口温度ばらつき導出部２０７（ステップＳ３０４の処理）により実現される。
原料供給量導出手段は、例えば、給炭量設定部２１３（ステップＳ３１１の処理）により実現される。
最大原料供給可能量導出手段は、例えば、最大給炭可能量導出部２１２（ステップＳ３０９の処理）により実現される。最大原料供給可能量は、例えば、最大給炭可能量に対応する。
熱風発生装置は、例えば、熱ガス発生装置１０１により実現される。
粉砕機は、例えば、ミル１０５により実現される。
捕集機は、例えば、バグフィルター１０７により実現される。
処理装置は、例えば、粉砕制御装置２００により実現される。
ミル出口温度ばらつき導出手段は、例えば、ミル出口温度ばらつき導出部２０６（ステップＳ３０３の処理）により実現される。
ミル入口温度設定値導出手段は、例えば、最大給炭可能量導出部２１２（ステップＳ３０５の処理）により実現される。
ミル入口温度ばらつき導出手段は、例えば、ミル入口温度ばらつき導出部２０７（ステップＳ３０４の処理）により実現される。
原料供給量導出手段は、例えば、給炭量設定部２１３（ステップＳ３１１の処理）により実現される。
最大原料供給可能量導出手段は、例えば、最大給炭可能量導出部２１２（ステップＳ３０９の処理）により実現される。最大原料供給可能量は、例えば、最大給炭可能量に対応する。
１０１：熱ガス発生装置、１０３：バンカー、１０４：給炭機、１０５：ミル、１１０：ミル出口温度計、２００：粉砕制御装置、２０１：ミル出口温度目標値記憶部、２０２：ミル出口温度測定値取得部、２０３：ミル出口温度偏差導出部、２０４：ＰＩＤ制御部、２０５：ミル出口温度記憶部、２０６：ミル出口温度ばらつき導出部、２０７：ミル入口温度ばらつき導出部、２０８：ミル入口温度上限値記憶部、２０９：外気温取得部、２１０：石炭水分取得部、２１１：給炭量取得部、２１２：最大給炭可能量導出部、２１３：給炭量設定部
Claims (8)
- 熱風を排ガスとして発生する熱風発生装置と、
原料を粉砕し、粉砕後の原料を、前記排ガスの流れに乗せて外部に放出する粉砕機と、
前記粉砕機に原料を供給する原料供給装置と、
前記粉砕機から前記排ガスの流れに乗って放出された粉砕後の原料を捕集する捕集機と、を有し、
前記粉砕機の出口側の所定の位置における温度であるミル出口温度の測定値と目標値との偏差に応じて前記熱風発生装置を操作する制御が行われる粉砕プラントに対する処理を行う、粉砕プラントにおける処理装置であって、
前記ミル出口温度のばらつき、または、前記ミル出口温度の測定値と目標値との偏差のばらつきをミル出口温度ばらつきとして導出するミル出口温度ばらつき導出手段と、
前記粉砕機の入口側の所定の位置における温度であるミル入口温度の上限値と、前記ミル出口温度ばらつき導出手段により導出された前記ミル出口温度ばらつきを前記ミル入口温度のばらつきに換算した値と、を用いて、前記ミル入口温度の設定値を導出するミル入口温度設定値導出手段と、を有し、
前記ミル入口温度設定値導出手段により導出された前記ミル入口温度の設定値は、前記ミル入口温度の測定値が前記ミル入口温度の設定値を上回らないように、前記原料供給装置から前記粉砕機に供給される前記原料の単位時間当たりの供給量を導出するために用いられることを特徴とする粉砕プラントにおける処理装置。 - 前記ミル出口温度ばらつき導出手段により導出された前記ミル出口温度ばらつきと、前記ミル出口温度を前記ミル入口温度に換算する係数と、を用いて、前記ミル入口温度のばらつきの推定値を、前記ミル出口温度ばらつきを前記ミル入口温度のばらつきに換算した値として導出するミル入口温度ばらつき導出手段をさらに有することを特徴とする請求項１に記載の粉砕プラントにおける処理装置。
- 前記ミル入口温度設定値導出手段は、前記ミル入口温度ばらつき導出手段により導出された前記ミル入口温度のばらつきの推定値と、前記ミル入口温度の上限値と、を用いて、前記ミル入口温度の測定値の所定の割合が当該上限値以下になるようにするための前記ミル入口温度の目標値を前記ミル入口温度の設定値として導出することを特徴とする請求項２に記載の粉砕プラントにおける処理装置。
- 前記ミル入口温度を前記ミル出口温度に換算する係数は、前記制御を行う制御器における制御ゲインを用いて導出される係数であることを特徴とする請求項２または３に記載の粉砕プラントにおける処理装置。
- 前記原料供給装置から前記粉砕機に供給される前記原料の単位時間当たりの供給量として、前記ミル入口温度の測定値が、前記ミル入口温度設定値導出手段により導出された前記ミル入口温度の設定値を上回らないようになる供給量を導出する原料供給量導出手段をさらに有することを特徴とする請求項１〜４の何れか１項に記載の粉砕プラントにおける処理装置。
- 前記ミル入口温度設定値導出手段により導出された前記ミル入口温度の設定値と、外気温と、前記原料に含まれる単位質量当たりの水分量と、を用いて、前記ミル入口温度の測定値が前記ミル入口温度の設定値を上回らない範囲で前記原料供給装置から前記粉砕機に供給できる原料の単位時間当たりの最大供給量である最大原料供給可能量を導出する最大原料供給可能量導出手段をさらに有し、
前記原料供給量導出手段は、前記最大原料供給可能量導出手段により導出された前記最大原料供給可能量に基づいて、前記供給量を導出することを特徴とする請求項５に記載の粉砕プラントにおける処理装置。 - 熱風を排ガスとして発生する熱風発生装置と、
原料を粉砕し、粉砕後の原料を、前記排ガスの流れに乗せて外部に放出する粉砕機と、
前記粉砕機に原料を供給する原料供給装置と、
前記粉砕機から前記排ガスの流れに乗って放出された粉砕後の原料を捕集する捕集機と、を有し、
前記粉砕機の出口側の所定の位置における温度であるミル出口温度の測定値と目標値との偏差に応じて前記熱風発生装置を操作する制御が行われる粉砕プラントに対する処理を行う、粉砕プラントにおける処理方法であって、
前記ミル出口温度のばらつき、または、前記ミル出口温度の測定値と目標値との偏差のばらつきをミル出口温度ばらつきとして導出するミル出口温度ばらつき導出工程と、
前記粉砕機の入口側の所定の位置における温度であるミル入口温度の上限値と、前記ミル出口温度ばらつき導出工程により導出された前記ミル出口温度ばらつきを前記ミル入口温度のばらつきに換算した値と、を用いて、前記ミル入口温度の設定値を導出するミル入口温度設定値導出工程と、を有し、
前記ミル入口温度設定値導出工程により導出された前記ミル入口温度の設定値は、前記ミル入口温度の測定値が前記ミル入口温度の設定値を上回らないように、前記原料供給装置から前記粉砕機に供給される前記原料の単位時間当たりの供給量を導出するために用いられることを特徴とする粉砕プラントにおける処理方法。 - 請求項１〜６の何れか１項に記載の粉砕プラントにおける処理装置の各手段としてコンピュータを機能させることを特徴とするプログラム。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2016249067A JP6776881B2 (ja) | 2016-12-22 | 2016-12-22 | 粉砕プラントにおける処理装置、方法、およびプログラム |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2016249067A JP6776881B2 (ja) | 2016-12-22 | 2016-12-22 | 粉砕プラントにおける処理装置、方法、およびプログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2018103067A JP2018103067A (ja) | 2018-07-05 |
| JP6776881B2 true JP6776881B2 (ja) | 2020-10-28 |
Family
ID=62784650
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2016249067A Active JP6776881B2 (ja) | 2016-12-22 | 2016-12-22 | 粉砕プラントにおける処理装置、方法、およびプログラム |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP6776881B2 (ja) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6798407B2 (ja) * | 2017-04-20 | 2020-12-09 | Ｊｆｅスチール株式会社 | 微粉炭の製造方法 |
| CN110124842B (zh) * | 2019-05-15 | 2020-11-17 | 广东电网有限责任公司 | 一种锅炉机组及其磨煤机出口温度控制方法与系统 |
| JP7311571B2 (ja) * | 2019-07-01 | 2023-07-19 | ホソカワミクロン株式会社 | 粉体処理システムの診断装置、診断方法、及びコンピュータプログラム |
| CN115400865B (zh) * | 2021-05-28 | 2023-11-14 | 北京京能电力股份有限公司 | 一种降低中速磨煤机最小出力的控制方法 |
| CN113898968B (zh) * | 2021-09-09 | 2023-06-27 | 华电电力科学研究院有限公司 | 一种防止直吹式制粉系统堵磨的平衡给煤控制方法 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4177951A (en) * | 1978-06-28 | 1979-12-11 | Combustion Engineering Inc. | Pulverizer air flow and temperature control |
| JPS60122059A (ja) * | 1983-12-07 | 1985-06-29 | 株式会社日立製作所 | 石炭ミル温度監視方式 |
| JP4638708B2 (ja) * | 2004-10-06 | 2011-02-23 | バブコック日立株式会社 | ミル装置ならびにそれを備えた石炭焚ボイラ設備 |
| JP5949414B2 (ja) * | 2012-10-05 | 2016-07-06 | 新日鐵住金株式会社 | 粉砕プラント排ガス制御装置、粉砕プラント排ガス制御方法、及びコンピュータプログラム |
| JP5962476B2 (ja) * | 2012-12-07 | 2016-08-03 | 新日鐵住金株式会社 | 粉砕プラント温度制御装置、粉砕プラント温度制御方法、及びコンピュータプログラム |
-
2016
- 2016-12-22 JP JP2016249067A patent/JP6776881B2/ja active Active
Also Published As
| Publication number | Publication date |
|---|---|
| JP2018103067A (ja) | 2018-07-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6776881B2 (ja) | 粉砕プラントにおける処理装置、方法、およびプログラム | |
| US20130030573A1 (en) | Computer-based method and device for automatically providing control parameters for a plurality of coal mills supplying coal powder to a plant | |
| CN101334666B (zh) | 双进双出钢球磨煤机直吹式制粉系统优化控制方法 | |
| EP2777813B1 (en) | Pulverizer monitoring | |
| CN110848733B (zh) | 一种基于煤质在线监测的燃烧优化方法 | |
| JP5962476B2 (ja) | 粉砕プラント温度制御装置、粉砕プラント温度制御方法、及びコンピュータプログラム | |
| CN110124842A (zh) | 一种锅炉机组及其磨煤机出口温度控制方法与系统 | |
| JP5949414B2 (ja) | 粉砕プラント排ガス制御装置、粉砕プラント排ガス制御方法、及びコンピュータプログラム | |
| CN102639937B (zh) | 用于监测和控制电厂的系统和关联的方法 | |
| CN106423526B (zh) | 钢球磨煤机风系统智能控制方法及模块 | |
| JP2017136540A (ja) | 固体燃料粉砕装置およびその制御方法 | |
| CN101407382A (zh) | 干法水泥余热利用中窑头余热锅炉取风工艺的控制方法及其装置 | |
| CN102046759B (zh) | 生产粉煤的方法 | |
| JP6702110B2 (ja) | 循環系粉砕プラントにおける粉砕機の入口温度の予測方法、装置及びプログラム、並びに原料供給量の演算方法、装置及びプログラム | |
| JP6492478B2 (ja) | 粉砕プラントにおける給炭量の演算方法、装置及びプログラム | |
| JP6801438B2 (ja) | 粉砕プラント酸素濃度制御装置、粉砕プラント酸素濃度制御方法、およびプログラム | |
| TW201934927A (zh) | 用於操作燃燒室之系統及方法 | |
| JP6729349B2 (ja) | 粉砕プラント酸素濃度制御装置、粉砕プラント酸素濃度制御方法、およびプログラム | |
| JP6702109B2 (ja) | 粉砕プラントにおけるプロセス制御方法、粉砕プラントにおけるプロセス制御装置、およびプログラム | |
| EP0943864B1 (en) | Combustion control method for refuse incinerator | |
| JP6347100B2 (ja) | 排ガス循環系粉砕プラントのミル出口温度制御方法、装置及びプログラム | |
| Dahl-Soerensen et al. | Pulverized fuel control using biased flow measurements | |
| Kurek et al. | Detection of malfunctions and abnormal working conditions of a coal mill | |
| Zagorodniy et al. | Expert diagnostics of slagging of heating surfaces of coal-fired boilers |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20190805 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20200422 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20200602 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20200722 |
|
| 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: 20200908 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20200921 |
|
| R151 | Written notification of patent or utility model registration |
Ref document number: 6776881 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |