(57)【特許請求の範囲】
【請求項１】
可燃物燃焼時の廃熱を回収して蒸気を発生させる廃熱回収型の焼却炉について、少なくとも可燃物の炉内への投入操作量を含む１つ以上の操作量を要素とする操作ベクトルに基づいて、１つ以上の制御量を要素とする制御ベクトルを予測演算する第１の予測演算手段を備えた焼却炉の制御装置において、
前記第１の予測演算手段による制御ベクトルの第１の予測値に基づいて、前記操作ベクトルに所定サンプル時間後までに変動がある場合の、前記操作ベクトルの変動分に起因した前記制御ベクトルの第１の変動予測値を演算する第１の変動予測値演算手段と、
少なくとも炉内に投入された可燃物の投入量に対応した予測補正量を要素とする予測補正ベクトルに基づいて、前記制御ベクトルを予測演算する第２の予測演算手段と、
前記第２の予測演算手段による制御ベクトルの第２の予測値に基づいて、前記操作ベクトルに前記所定サンプル時間後までに変動がない場合の前記制御ベクトルの第２の変動予測値を演算する第２の変動予測値演算手段と、
前記制御ベクトルの実測値、前記第１の変動予測値、及び前記第２の変動予測値に基づいて、前記制御ベクトルの予測値を演算する予測演算手段とを具備してなり、
前記第１の予測演算手段が、次式（１）で表されるモデルを用いるものであり、前記第１の変動予測値演算手段が、前記第１の変動予測値ｙ１（ｋ＋ｉ）の導出に、次式（２）及び（３）を用いるものであり、前記第２の予測演算手段が、次式（４）で表されるモデルを用いるものであり、前記第２の変動予測値演算手段が、前記第２の変動予測値ｙ２（ｋ＋ｉ）の導出に、次式（５）及び（６）を用いるものであり、前記予測演算手段が、前記制御量の予測値ｙｐ（ｋ＋ｉ）の導出に、次式（７）を用いるものであることを特徴とする焼却炉の制御装置。
ｙＭ１（ｋ＋ｉ）＝ｙＭ１（ｋ）
＋Ａｉ（ｚ−１）ΔｙＭ１（ｋ）＋Ｂｉ（ｚ−１）Δｕ１（ｋ＋ｉ） （１）
ｙ１（ｋ＋ｉ）＝ｙＭ１（ｋ＋ｉ）−ｙＭ１（ｋ） （２）
ΔｙＭ１（ｋ−ｊ）＝０、Δｕ１（ｋ−ｊ−１）＝０ （３）
ｙＭ２（ｋ＋ｉ）＝ｙＭ２（ｋ）
＋Ｇｉ（ｚ−１）ΔｙＭ２（ｋ）＋Ｈｉ（ｚ−１）Δｕ２（ｋ＋ｉ） （４）
ｙ２（ｋ＋ｉ）＝ｙＭ２（ｋ＋ｉ）−ｙＭ２（ｋ） （５）
Δｕ２（ｋ＋ｊ）＝０ （６）
ｙｐ（ｋ＋ｉ）＝ｙ（ｋ）＋ｙ１（ｋ＋ｉ）＋ｙ１（ｋ＋ｉ） （７）
ただし、ｋは現時点、ｕ1 は前記操作ベクトル、ｙＭ１は操作量ｕ1 を用いた制御ベクトルの第１の予測値、Ａｉ（ｚ−１）、Ｂｉ（ｚ−１）は遅延演算子ｚ−１に関する多項式行列、ｕ２は予測補正ベクトル、ｙＭ２は予測補正ベクトルｕ２を用いた制御ベクトルの第２の予測値、ｙは制御ベクトルの実測値、Ｇｉ（ｚ−１）、Ｈｉ（ｚ−１）は遅延演算子ｚ−１に関する多項式行列、ｉ＝０、１、２、…、Δは前時刻に対する変量の偏差ベクトル（差分演算子）である。
【請求項２】
前記制御ベクトルが、ボイラドラム圧力と蒸気流量を要素とするベクトルであって、前記操作ベクトルが、炉内に可燃物を投入するための給塵ロータの回転速度と蒸気弁開度を要素とするベクトルであって、前記予測補正ベクトルが、炉内酸素濃度と蒸気弁開度を要素とするベクトルである請求項１に記載の焼却炉の制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、焼却炉の制御装置に係り、詳しくは、可燃物燃焼時の廃熱を回収して蒸気を発生させる廃熱回収型の焼却炉において、ボイラドラム圧力、流量などを要素とする制御ベクトルを、給塵ロータ回転速度、蒸気弁開度などを操作して制御する焼却炉の制御装置に関するものである。
【背景技術】
【０００２】
近年、地球環境問題が深刻化する中、未利用エネルギーの有効活用が叫ばれており、流動床式焼却炉のような都市ごみ焼却炉においても、ごみなどの可燃物を燃焼させた際に生じる廃熱を用いてボイラで蒸気を発生させエネルギーを回収する熱回収型の焼却炉の開発研究が精力的に行われている。このような熱回収型の焼却炉では、ボイラから発生する蒸気を如何に安定して回収するかが操業上重要である。
過去、蒸気系においては、ＰＩＤ補償器によって、蒸気弁や給塵ロータを操作し、ボイラドラム圧力や蒸気流量を制御する構成を用いることが多かった。しかしながら、これら複数の操作量と制御量との間には干渉が存在し、また、一部にむだ時間も存在することから、これらの系に対しては、多変数モデル予測制御を用いるのが有効であり、本発明者はこれまで多変数モデル予測制御を用いた制御システムを構築してきた（特許文献１等を参照）。
例えば、ボイラドラム圧力、蒸気流量等の複数の制御量を要素とする制御ベクトルの挙動を予測するモデルを作成するために、実炉にて同定実験を行い、部分空間法により同定したモデルを式変形すると、ｋを現サンプル時点、ｚ−１を遅延演算子としたとき、給塵ロータ回転速度、蒸気弁開度などの複数の操作量を要素とする操作ベクトルと上記制御ベクトルとの関係は、次式（８）により表すことができる。
ｙＭ１（ｋ＋ｉ）＝ｙＭ１（ｋ）
＋Ａｉ（ｚ−１）ΔｙＭ１（ｋ）＋Ｂｉ（ｚ−１）Δｕ１（ｋ＋ｉ） （８）
ただし、ｋは現時点、ｕ1 は操作ベクトル、ｙＭ１は操作ベクトルｕ1 を用いた制御ベクトルの予測値（第１の予測値）、Ａｉ（ｚ−１）、Ｂｉ（ｚ−１）はｚ−１に関する多項式行列、ｉ＝０、１、２、…、Δは前時刻に対する変量の偏差ベクトル、すなわち差分演算子である。
ここで、給塵ロータの回転速度と、蒸気弁開度を操作量としたときのボイラドラム圧力の時間的推移について、上記モデル式（８）による計算値と実測値とを比較した比較例を図５に示す。
図５から理解される通り、上記モデル式（８）は、ボイラドラム圧力、ひいては制御ベクトルのおおよそ挙動を表現することができる。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開平１０−３１８５０３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ところで，多種多様な都市ごみを処理する都市ごみ焼却炉や，これらの都市ごみなどを乾燥させて固形化したごみ固形化燃料，いわゆるＲＤＦ(Refuse Derived fuel) を処理する処理炉などの焼却炉においては，前記操作ベクトルを変化させていないにもかかわらず，ごみ質の変動などに伴って突発的に給塵量がばらつく場合があり，その影響を考慮する必要がある。しかしながら，一般的な１入力１出力の制御系はもちろん，前記したモデル式（８）に基づく多変数モデル予測制御系においても，実際の給塵量のばらつきが反映されていないため，制御量の挙動の表現精度に限界があった。
本発明は，このような従来の技術における課題を解決するために，焼却炉の制御装置を改良し，少なくとも可燃物の炉内への投入操作量を含む１つ以上の操作量を要素とする操作ベクトルに応じた第１のモデルを用いて求めた制御ベクトルの第１の予測値に基づいて，前記操作ベクトルに所定サンプル時間後までに変動がある場合の，前記操作ベクトルの変動分に起因した前記制御ベクトルの第１の変動予測値を求め，少なくとも炉内に投入された可燃物の投入量に対応した予測補正量を要素とする予測補正ベクトルに応じた第２のモデルを用いて求めた前記制御ベクトルの第２の予測値に基づいて，前記操作ベクトルに前記所定サンプル時間後までに変動がない場合の前記制御ベクトルの第２の変動予測値を求め，前記制御ベクトルの実測値，前記第１の変動予測値，及び前記第２の変動予測値に基づいて，前記制御ベクトルの予測値を求め，前記制御ベクトルの予測値に基づいて前記制御を行うことにより，実際の給塵量にばらつきがある場合でも，精度の良い安定した制御を行うことができる焼却炉の制御装置を提供するものである。
【課題を解決するための手段】
【０００５】
前記目的を達成するために、請求項１に係る発明は、可燃物燃焼時の廃熱を回収して蒸気を発生させる廃熱回収型の焼却炉について、少なくとも可燃物の炉内への投入操作量を含む１つ以上の操作量を要素とする操作ベクトルに基づいて、１つ以上の制御量を要素とする制御ベクトルを予測演算する第１の予測演算手段を備えた焼却炉の制御装置において、前記第１の予測演算手段による制御ベクトルの第１の予測値に基づいて、前記操作ベクトルに所定サンプル時間後までに変動がある場合の、前記操作ベクトルの変動分に起因した前記制御ベクトルの第１の変動予測値を演算する第１の変動予測値演算手段と、少なくとも炉内に投入された可燃物の投入量に対応した予測補正量を要素とする予測補正ベクトルに基づいて、前記制御ベクトルを予測演算する第２の予測演算手段と、前記第２の予測演算手段による制御ベクトルの第２の予測値に基づいて、前記操作ベクトルに前記所定サンプル時間後までに変動がない場合の前記制御ベクトルの第２の変動予測値を演算する第２の変動予測値演算手段と、前記制御ベクトルの実測値、前記第１の変動予測値、及び前記第２の変動予測値に基づいて、前記制御ベクトルの予測値を演算する予測演算手段とを具備してなり、前記第１の予測演算手段が、次式（１）で表されるモデルを用いるものであり、前記第１の変動予測値演算手段が、前記第１の変動予測値ｙ１（ｋ＋ｉ）の導出に、次式（２）及び（３）を用いるものであり、前記第２の予測演算手段が、次式（４）で表されるモデルを用いるものであり、前記第２の変動予測値演算手段が、前記第２の変動予測値ｙ２（ｋ＋ｉ）の導出に、次式（５）及び（６）を用いるものであり、前記予測演算手段が、前記制御量の予測値ｙｐ（ｋ＋ｉ）の導出に、次式（７）を用いるものであることを特徴とする焼却炉の制御装置である。
ｙＭ１（ｋ＋ｉ）＝ｙＭ１（ｋ）
＋Ａｉ（ｚ−１）ΔｙＭ１（ｋ）＋Ｂｉ（ｚ−１）Δｕ１（ｋ＋ｉ） （１）
ｙ１（ｋ＋ｉ）＝ｙＭ１（ｋ＋ｉ）−ｙＭ１（ｋ） （２）
ΔｙＭ１（ｋ−ｊ）＝０、Δｕ１（ｋ−ｊ−１）＝０ （３）
ｙＭ２（ｋ＋ｉ）＝ｙＭ２（ｋ）
＋Ｇｉ（ｚ−１）ΔｙＭ２（ｋ）＋Ｈｉ（ｚ−１）Δｕ２（ｋ＋ｉ） （４）
ｙ２（ｋ＋ｉ）＝ｙＭ２（ｋ＋ｉ）−ｙＭ２（ｋ） （５）
Δｕ２（ｋ＋ｊ）＝０ （６）
ｙｐ（ｋ＋ｉ）＝ｙ（ｋ）＋ｙ１（ｋ＋ｉ）＋ｙ１（ｋ＋ｉ） （７）
ただし、ｋは現時点、ｕ1 は前記操作ベクトル、ｙＭ１は操作量ｕ1 を用いた制御ベクトルの第１の予測値、Ａｉ（ｚ−１）、Ｂｉ（ｚ−１）は遅延演算子ｚ−１に関する多項式行列、ｕ２は予測補正ベクトル、ｙＭ２は予測補正ベクトルｕ２を用いた制御ベクトルの第２の予測値、ｙは制御ベクトルの実測値、Ｇｉ（ｚ−１）、Ｈｉ（ｚ−１）は遅延演算子ｚ−１に関する多項式行列、ｉ＝０、１、２、…、Δは前時刻に対する変量の偏差ベクトル（差分演算子）である。
【０００６】
また，請求項２に係る発明は，前記請求項１に記載の焼却炉の制御装置において，前記制御ベクトルが，ボイラドラム圧力と蒸気流量を要素とするベクトルであって，前記操作ベクトルが，炉内に可燃物を投入するための給塵ロータの回転速度と蒸気弁開度を要素とするベクトルであって，前記予測補正ベクトルが，炉内酸素濃度と蒸気弁開度を要素とするベクトルであることをその要旨とする。
前記請求項１又は２のいずれか１項に記載の焼却炉の制御装置では，少なくとも可燃物の炉内への投入操作量を含む１つ以上の操作量を要素とする操作ベクトルに応じた第１の予測演算手段を用いて求めた制御ベクトルの第１の予測値に基づいて，前記操作ベクトルに所定サンプル時間後までに変動がある場合の，前記操作ベクトルの変動分に起因した前記制御ベクトルの第１の変動予測値を演算し，少なくとも炉内に投入された可燃物の投入量に対応した予測補正量を要素とする予測補正ベクトルに応じた第２の予測演算手段を用いて求めた前記制御ベクトルの第２の予測値に基づいて，前記操作ベクトルに前記所定サンプル時間後までに変動がない場合の前記制御ベクトルの第２の変動予測値を演算し，前記制御ベクトルの実測値，前記第１の変動予測値，及び前記第２の変動予測値に基づいて，予測演算手段が，前記制御ベクトルの予測値を演算し，前記制御ベクトルの予測値に基づいて前記制御を行うことにより，実際の給塵量にばらつきがある場合でも，精度の良い安定した制御を行うことが可能な焼却炉の制御装置を提供することができる。
【発明の効果】
【０００７】
前記請求項１又は２のいずれか１項に記載の焼却炉の制御装置では，少なくとも可燃物の炉内への投入操作量を含む１つ以上の操作量を要素とする操作ベクトルに応じた第１の予測演算手段を用いて求めた制御ベクトルの第１の予測値に基づいて，前記操作ベクトルに所定サンプル時間後までに変動がある場合の，前記操作ベクトルの変動分に起因した前記制御ベクトルの第１の変動予測値を演算し，少なくとも炉内に投入された可燃物の投入量に対応した予測補正量を要素とする予測補正ベクトルに応じた第２の予測演算手段を用いて求めた前記制御ベクトルの第２の予測値に基づいて，前記操作ベクトルに前記所定サンプル時間後までに変動がない場合の前記制御ベクトルの第２の変動予測値を演算し，前記制御ベクトルの実測値，前記第１の変動予測値，及び前記第２の変動予測値に基づいて，予測演算手段が，前記制御ベクトルの予測値を演算し，前記制御ベクトルの予測値に基づいて前記制御を行うことにより，実際の給塵量にばらつきがある場合でも，精度の良い安定した制御を行うことが可能な焼却炉の制御装置を提供することができる。
【図面の簡単な説明】
【０００８】
【図１】本発明の実施の形態に係る焼却炉の制御方法が適用可能な制御装置の概略構成を示す図。
【図２】前記焼却炉の制御装置が備えるモデル予測制御器の概略構成を示す機能ブロック図。
【図３】予測補正ベクトルを用いた制御量のモデルの精度を説明するためのタイムチャート。
【図４】前記制御装置の制御結果を説明するためのタイムチャート。
【図５】従来の制御量のモデルの精度を説明するためのタイムチャート。
【発明を実施するための形態】
【０００９】
以下，添付図面を参照して，本発明の実施の形態につき説明し，本発明の理解に供する。尚，以下の実施の形態は，本発明の具体的な一例であって，本発明の技術的範囲を限定する性格のものではない。ここに，図１は本発明の実施の形態に係る焼却炉の制御方法が適用可能な焼却炉の制御装置の概略構成を示す図である。
本発明の実施の形態に係る焼却炉の制御方法は，可燃物燃焼時の廃熱を利用して発生させた蒸気から熱エネルギーを回収する焼却炉について，少なくとも可燃物の炉内への投入操作量を含む複数の操作量を要素とする操作ベクトルに基づいて，複数の制御量を要素とする制御ベクトルを予測演算する第１のモデルを用いて制御を行う点で従来のものと同様である。
一方，本発明の実施の形態に係る焼却炉の制御方法は，前記第１のモデルによる前記制御ベクトルの第１の予測値に基づいて，前記操作ベクトルに所定サンプル時間後までに変動がある場合の，前記操作ベクトルの変動分に起因した前記制御ベクトルの第１の変動予測値を演算し，少なくとも炉内に投入された可燃物の投入量に対応した予測補正量を要素とする予測補正ベクトルに基づいて，前記制御ベクトルを予測演算する第２のモデルによる，前記制御ベクトルの第２の予測値に基づいて，前記操作ベクトルに前記所定サンプル時間後までに変動がない場合の前記制御ベクトルの第２の変動予測値を演算し，前記制御ベクトルの実測値，前記第１の変動予測値，及び前記第２の変動予測値に基づいて，前記制御ベクトルの予測値を演算し，前記制御ベクトルの予測値に基づいて前記制御を行う点で従来のものと異なる。
【００１０】
本発明の実施の形態に係る焼却炉の制御方法は，例えば図１に示すような焼却炉の制御装置によって具体化される。
図１において，ごみ等の可燃物は，給塵ロータ１により炉内に投入され，砂層２内に落下する。砂層２には，一次空気が複数部位から送り込まれており，それによって砂が流動するとともに，前記投入された可燃物が砂層２内で一次燃焼し，発生したガスはフリーボード３に流れてゆく。フリーボード３には，２次空気が送り込まれ，ＣＯ等の未燃ガスが完全燃焼させられ，排ガスは排ガス出口４に送られる。その際，廃熱ボイラ５では，排ガスの熱エネルギーによって飽和蒸気が発生する。廃熱ボイラ５で発生した発生蒸気は，図示しないタービン等のエネルギー変換手段へと送られる。その蒸気の流量は蒸気弁６により調整され，該蒸気のボイラドラム圧力，流量はそれぞれ圧力計７，流量計８により測定される。前記圧力計７，流量計８の出力値は，コンピュータ等の演算手段を用いたモデル予測制御器１０に供給される。
前記モデル予測制御器１０は，例えば給塵ロータ１の回転速度，及び蒸気弁６の弁開度を操作し，ボイラドラム圧力，及び蒸気流量を予測制御するものであって，図２に示す如く，少なくとも可燃物の炉内への投入操作量を含む複数の操作量を要素とする操作ベクトルｕ1 に基づいて，複数の制御量を要素とする制御ベクトルを予測演算する第１の予測演算手段１０１と，前記第１の予測演算手段１０１による制御ベクトルの第１の予測値ｙM1に基づいて，前記操作ベクトルｕ1 に所定サンプル時間後までに変動がある場合の，前記操作ベクトルｕ1 の変動分Δｕ1 に起因した制御ベクトルｙp の第１の変動予測値ｙ1 を演算する第１の変動予測値演算手段１０２と，少なくとも炉内に投入された可燃物の投入量に対応した炉内酸素濃度（予測補正量）を要素とする予測補正ベクトルｕ2 に基づいて，前記制御ベクトルを予測演算する第２の予測演算手段１０３と，前記第２の予測演算手段１０３による前記制御ベクトルの第２の予測値ｙM2に基づいて，前記操作ベクトルｕ1 に前記所定サンプル時間後までに変動がない場合の制御ベクトルの第２の変動予測値ｙ2 を演算する第２の変動予測値演算手段１０４と，前記制御ベクトルの実測値ｙ，前記第１の変動予測値ｙ1 ，及び前記第２の変動予測値ｙ2 に基づいて，前記制御ベクトルの予測値ｙp を演算する予測演算手段１０５とを具備する。
【００１１】
前記モデル予測制御器１０において，第１の予測演算手段１０１で用いられるモデル（第１のモデルに相当）が，ボイラドラム圧力，及び蒸気流量を要素とする制御ベクトルの第１の予測値ｙM1を予測するための上式（８）（請求項１における式（１）と同じ式）である。ただし，上式（８）では，炉内への実際の給塵量（投入量）が反映されない。これは，平均的にみたときには給塵ロータ１の回転速度は実際のごみ投入量に比例するものの，ごみ投入量を短周期的にみたときには，前記給塵ロータ１の回転速度が一定であっても，極端な場合「ごみ切れ」や「どか落ち」などの変動が起こりうるからであり，上式（８）に従った制御ベクトルの予測精度には限界がある。
前記モデル予測制御器１０では，これを補うために，炉内酸素濃度を測定する炉内酸素計９が接続されている。炉内酸素濃度は，炉内に実際に投入された可燃物の投入量と相関（対応）があって，炉内に投入された可燃物の実際の投入量が変動したときに，制御ベクトルに用いられるボイラドラム圧力や蒸気流量と比べて，その値の変動が検出されるまでに必要な時間が短い。このため，前記「ごみ切れ」や「どか落ち」などのごみ投入量の短周期的な変動がある場合でも，該炉内酸素計９により測定された炉内酸素濃度（予測補正量の一例）を要素とする予測補正ベクトルｕ2 に基づいて，制御ベクトルを予測すれば，制御ベクトルの予測精度を向上させることが可能である。
上式（８）を同定する際に，給塵ロータ１の回転速度の代わりに，炉内酸素計９の出力値を用いると，次式（９）（請求項１における式（４）と同じ式）であらわされるモデル（第２のモデル）を得ることができる。前記第２の予測演算手段１０３は，次式（９）で表されるモデルを用いて制御量の演算を行うものである。
ｙM2（ｋ＋ｉ）＝ｙM2（ｋ）
＋Ｇi （ｚ-1）ΔｙM2（ｋ）＋Ｈi （ｚ-1）Δｕ2 （ｋ＋ｉ）  （９）
ただし，ｕ2 は予測補正ベクトル，ｙM2は予測補正ベクトルｕ2 を用いた制御ベクトルの第２の予測値，Ｇi （ｚ-1），Ｈi （ｚ-1）は遅延演算子ｚ-1に関する多項式行列である。
ここで，炉内酸素計９から出力された炉内酸度濃度の他，蒸気弁開度を予測補正ベクトルｕ2 の要素としたときのボイラドラム圧力の時間推移について，前記モデル式（９）による計算値と実測値とを比較した比較例を図３に示す。
図３から理解される通り，前記モデル式（９）は，ボイラドラム圧力，ひいては制御ベクトルを精度良く表現することができる。
【００１２】
ところで、モデル予測制御アルゴリズムにおける予測式の構造は、一般に次式（１０）（請求項１における式（７）と同じ式）により表現することができる。
ｙｐ（ｋ＋ｉ）＝ｙ（ｋ）＋ｙ１（ｋ＋ｉ）＋ｙ２（ｋ＋ｉ） （１０）
ただし、ｙｐは現時点からｉサンプル後の制御ベクトルの予測値、ｙは制御ベクトルの実測値、ｙ１は現時点からｉサンプル後までに操作ベクトルｕ１に変動があった場合に、操作ベクトルｕ１の変動分に起因する制御ベクトルの予測変動値（第１の変動予測値）、ｙ２は現時点からｉサンプル後までに操作量ｕ１に変動がない場合の制御ベクトルの予測変動値（第２の変動予測値）である。
そこで、上式（８）を用いてｙ１を計算し、上式（９）を用いてｙ２を計算することにより、制御量の予測値ｙｐを求めることを考える。つまり、ｉ＝０、１、２、…に対し、次式（１１）及び（１２）（請求項１における式（２）、（３）と同じ式）によりｙ１を与え、
ｙ１（ｋ＋ｉ）＝ｙＭ１（ｋ＋ｉ）−ｙＭ１（ｋ） （１１）
ΔｙＭ１（ｋ−ｊ）＝０、Δｕ１（ｋ−ｊ−１）＝０ （１２）
次式（１３）及び（１４）（請求項１における式（５）、（６）と同じ式）によりｙ２を与え、
ｙ２（ｋ＋ｉ）＝ｙＭ２（ｋ＋ｉ）−ｙＭ２（ｋ） （１３）
Δｕ２（ｋ＋ｊ）＝０ （１４）
上式（１０）の制御ベクトルの予測値ｙｐを求める。
上式（８）、（１１）、及び（１２）を用いて、第１の変動予測値ｙ１を求めるのが、第１の変動予測演算手段１０２であり、上式（９）、（１３）、及び（１４）を用いて、第２の変動予測値ｙ２を求めるのが、第２の変動予測演算手段１０４である。
前記第１の変動予測演算手段１０２及び第２の変動予測演算手段によりそれぞれ求められた変動予測値ｙ１、ｙ２は、予測演算手段１０５に供給される。
前記予測演算手段１０５では、供給された制御ベクトルの実測値ｙ、前記第１の変動予測値ｙ１、及び前記第２の変動予測値ｙ２に基づき、上式（１０）に従って、前記制御ベクトルの予測値ｙｐが演算される。
そして、上式（１０）について、一般的なモデル予測制御アルゴリズムと同様に、目標軌道ｙＲと予測値ｙｐが、未来の一定区間においてできるだけ一致するような偏差入力Δｕ１（ｋ）が、最小自乗法や２次計画法を用いて求められる。
このようにして求めた偏差入力Δｕ１（ｋ）に従って、前記予測モデル制御器１０は、前記操作ベクトルｕ１を操作し制御を行う。
【００１３】
現在操業中の実炉において，適用実験を行った結果を図４に示す。ここで，図４（ａ）は制御ベクトルの要素の一つであるボイラドラム圧力の目標値と実測値に関するタイムチャートであり，図４（ｂ）は制御ベクトルの他の要素である蒸気流量の目標値と実測値に関するタイムチャートであり，図４（ｃ）は給塵ロータの回転速度と蒸気弁開度の操作設定値に関するタイムチャートであり，図４（ｄ）は炉内酸素計出力値のタイムチャートである。
図４の適用実験では，制御周期は３秒に設定され，予測ホライズンは１００ステップ，操作ホライズンは１ステップである。
図４から理解される通り，ボイラドラム圧力も蒸気流量も目標値に安定に且つ精度良く制御されている。
このように，本実施の形態に係る焼却炉の制御方法及びその装置によれば，少なくとも可燃物の炉内への投入操作量を含む複数の操作量を要素とする操作ベクトルに応じた制御ベクトルの第１の予測値に基づいて，前記操作ベクトルに所定サンプル時間後までに変動がある場合の，前記操作ベクトルの変動分に起因した制御ベクトルの第１の変動予測値を演算し，少なくとも炉内に投入された可燃物の投入量に対応した予測補正量を要素とする予測補正ベクトルに応じた制御ベクトルの第２の予測値に基づいて，前記操作ベクトルに前記所定サンプル時間後までに変動がない場合の制御ベクトルの第２の変動予測値を演算し，制御ベクトルの実測値，前記第１の変動予測値，及び前記第２の変動予測値に基づいて，前記制御ベクトルの予測値を演算し，前記制御ベクトルの予測値に基づいて前記制御を行うことにより，実際の給塵量にばらつきがある場合でも，精度の良い安定した制御を行うことが可能な焼却炉の制御方法及びその装置を提供することができる。
【００１４】
尚，前記実施の形態では，流動床炉に本発明を適用したが，これに限られるものではなく，ストーカ炉や，ガス化溶融炉のような他の炉にも適用することは可能である。
さらに，本発明は，いわゆる都市ごみを処理する炉だけでなく，ごみを乾燥させて固形化したごみ固形化燃料ＲＤＦや石炭などの燃料も含む可燃物全般を処理する炉に適用することも可能である。
また，前記実施の形態では，操作ベクトルとして，ロータ回転速度と蒸気弁開度を用い，制御ベクトルとして，ボイラドラム圧力と蒸気流量を用い，予測補正ベクトルとして，炉内酸素濃度と蒸気弁開度を用いたが，これに限られるものではない。
前記予測補正ベクトルのうち，制御ベクトルの予測精度を向上させるために用いる，前記炉内酸素濃度に代わる予測補正量の選択条件としては，次の２つの条件（Ａ），及び（Ｂ）が挙げられる。
（Ａ）  炉内に投入された実際のごみ投入量を反映する（相関）があること
（Ｂ）  ごみ投入量の変動に対する影響が，制御ベクトルよりもはやく現れること
前記条件（Ａ），（Ｂ）を満たしていれば，「ごみ切れ」や「どか落ち」など給塵ロータ１の回転速度などにあらわれないごみ投入量の短周期的な変動があった場合でも，制御ベクトルの挙動を精度よく予測することが可能であり，炉内酸素濃度に限らず，予測補正量として，例えば炉内温度，炉内圧力，ロータ油圧，排ガス流量等を用いることが可能である。また，予測補正量に，投入されるごみや炉内火炎の画像評価値（ごみの投入量に対応）を用いてもよい。もちろん，実際に投入されるごみ投入量が測定できる場合には，ごみ投入量の測定値そのものを第２のモデルの入力（予測補正量）とすればよい。
【００１５】
さらに，前記制御ベクトルには，前記ボイラドラム圧力，蒸気流量の他，蒸気温度，砂層温度（流動床炉の場合），炉内温度（フリーボード部温度や，炉頂温度などを含む），炉内酸素濃度などを用いてもよいし，ボイラドラム圧力，蒸気流量も含めて，それらの組み合わせを用いることも可能である。但し，制御ベクトルの要素に，予測補正量に用いることができるような検出時間が他の制御量とくらべて相対的に短い，例えば炉内酸素濃度や炉内温度などを用いた場合には，これらを予測補正量として用いることはできず，これらよりもごみ投入量の変動に対する影響がはやく検出される変量を少なくとも一つは予測補正ベクトルに用いる必要がある。例えば画像やマイクロ波センサなどを用いて，炉内への実際の可燃物の投入量を検出した場合，炉内酸素濃度や炉内温度よりも，ごみ投入量の変動に対する影響がその検出値にはやく現れるので，炉内酸素濃度や炉内温度を制御ベクトルの要素に用いた場合には，画像やマイクロ波センサなどにより検出された炉内への実際の可燃物の投入量を予測補正量に用いる。
さらに，前記操作ベクトルも，給塵ロータ１の回転速度や蒸気弁開度に限られるものではなく，給塵ロータ１の代わりにベルトコンベアやスクリューコンベアを用いた場合には，ベルトコンベア速度やスクリューコンベア速度を用いることが可能であるし，ごみを後方から押して給塵するプッシャーが備えられている場合には，プッシャー速度を用いることも可能であり，ごみ投入量を調節する操作量であれば，他の操作量も同様に用いることができる。
また，前記実施の形態では，第２のモデルに上式（９）を用い，第２の予測演算手段１０３も上式（９）を用いて制御ベクトルの予測を行ったが，第２のモデルは，所定サンプル時間後までに変動がない場合の制御ベクトルを予測するものであるから，数式モデルである必要はなく，熱収支・物質収支などに基づく数値計算プログラムを用いたシミュレータであってもよいし，数量変化に対するオペレータの知識（ノウハウ）を蓄積した知識ベースであってもよい。
また，前記実施の形態では，２入力２出力の多変数モデル予測制御系について本発明を適用したが，本発明の構成は制御系の入力数と出力数に関係なく，それ以外の多変数モデル予測制御系に本発明を適用することも可能であるし，１入力１出力の制御系に本発明を適用することも可能である。
【符号の説明】
【００１６】
１…給塵ロータ
５…ボイラ
６…蒸気弁
７…圧力計
８…流量計
９…炉内酸素計
１０…モデル予測制御器
１０１…第１の予測演算手段
１０２…第１の変動予測演算手段
１０３…第２の予測演算手段
１０４…第２の変動予測演算手段
１０５…予測演算手段
【図１】
0004942788.tif 000002
【図２】
0004942788.tif 000003
【図３】
0004942788.tif 000004
【図４】
0004942788.tif 000005
【図５】
0004942788.tif 000006