JP2000078749A - 発電プラントの最適化制御方法及び最適化制御装置 - Google Patents

Abstract

(57)【要約】 【課題】 時々刻々と変化する蒸気、電気の需要変化に
対応してボイラ、タービンを常に最適な状態で運転でき
るよう発電プラントを効率的に制御する最適化制御方法
及び装置を提供する。 【解決手段】 並列化された複数の発電要素を有する発
電プラントにおいて、該複数の発電要素を最適負荷配分
により制御する最適化制御方法であって、上記複数の発
電要素の各々に対する最適負荷配分の演算に必要なパラ
メータのうち、直接測定できない発電要素の特性や状態
に関わるパラメータを、発電要素の物理モデルと入力プ
ロセス状態量とに基づき演算により推定し、得られた各
発電要素のパラメータを用いて、発電プラント全体のコ
ストを最小にする各発電要素に対する最適な負荷配分を
決定して各発電要素を制御することを特徴とする発電プ
ラントの最適化制御方法。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は、電気や蒸気等のエ
ネルギを工場等に供給する発電プラントを効率的に制御
する最適化制御方法及び最適化制御装置に関する。本発
明は特に、発電プラントが複数のボイラやタービンを並
列運転するものである場合に、並列のボイラ間やタービ
ン間での負荷配分を最適化するための有効な手段を提供
する。

【０００２】

【従来の技術】図５は、石油化学工場等で採用されてい
る発電プラントの一例の構成図であり、重油や燃料ガス
等を燃料とするボイラ１，２と、蒸気及び電気を発生す
るタービン３，４と、発電専用の復水タービン５、エネ
ルギを電気に変換する発電機６，７，８と、ボイラ主蒸
気母管２０と、工場２３に蒸気及び電気を供給する蒸気
母管２１、電気母線２２とを有している。

【０００３】この発電プラントの一部にはＰＩＤ制御装
置が設けてある。すなわち、ボイラ主蒸気圧力を制御す
るＰＩＤ制御装置１３、ボイラ２の燃料流量を制御する
ＰＩＤ制御装置１４、蒸気圧力を制御するＰＩＤ制御装
置１５、タービン４の蒸気流量を制御するＰＩＤ制御装
置１６、発電機８の発電量を制御するＰＩＤ制御装置１
７が設けてあり、これらによってエネルギの需給バラン
スを保つようにしている。なお、図を分りやすくするた
め、ボイラやタービンに対するＰＩＤ制御装置は一部の
もののみが示されている。

【０００４】また、ボイラの負荷配分は、運転員が燃料
流量を制御するＰＩＤ制御装置１４の設定値を手動で調
整することによってなされ、タービンの負荷配分は蒸気
流量を制御するＰＩＤ制御装置１６の設定値を手動で調
整することによってなされ、電力会社からの受電量の調
整は発電量を制御するＰＩＤ制御装置１７の設定値を手
動で調整することによってなされている。このような手
動調整の際に、プラント全体としての経済性を考慮した
最適な運転指針を定周期で運転員に提示することが試み
られている（例えば、ＳｔｅａｍＢａｌａｎｃｅ Ｏｐ
ｔｉｍｉｚａｔｉｏｎ ｉｎ Ｃｈｅｍｉｃａｌ Ｐｌ
ａｎｔ，ｉｎ Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ＩＦＡ
Ｃ ＳｙｍｐｏｓｉｕｍＡＤＣＨＥＭ '９４（１９９
４））。これは、プラントの線形モデルをもとに最適計
算を実施するが、その結果を基にした運転調整は運転員
に委ねられていた。

【０００５】

【発明が解決しようとする課題】しかしながら、上記制
御システムでは次の３つの問題が発生した。第一には、
発電プラント全体としての経済性を考慮し、運転員が時
々刻々と変化する蒸気、電気の需要変化に応じボイラ、
タービンを常に最適な状態で運転すべく頻繁に調整する
ことは精神的疲労が大きく、きわめて困難であり、ラフ
な運転とならざるを得なかった。第二には、発電プラン
ト全体の経済性を評価するモデルは、過去の実績データ
をもとに得られた線形モデルであったため、プラントの
非線形特性を考慮できず運転条件が変化すると、真の最
適運転条件から乖離することがあった。第三には、ボイ
ラ１，２及びタービン３，４，５は、ボイラ主蒸気母管
２０、蒸気母管２１を介して密接に関連しあっているこ
と、発電機６，７，８はタービン３，４，５にそれぞれ
依存すること、また受電量と発電量が電気母線を介して
関連していることから、発電プラントのそれぞれの制御
システムが相互に干渉するため発電プラント全体として
電気、蒸気の需給バランスの変動を短時間のうちに吸収
することはきわめて困難であった。

【０００６】

【課題を解決するための手段】上記課題を解決するた
め、本発明は、ボイラやタービンなどの並列化された複
数の発電要素からなる発電プラントにおいて、複数の発
電要素の各々に対する最適負荷配分を自動的に演算して
制御を行なうようにするものである。この場合、最適負
荷配分演算では、タービン発電機の損失パラメータや熱
交換器の汚れ係数などのように、プラントからは直接測
定できない発電要素の特性や状態に関わるパラメータが
必要となる。本発明はこれらの直接測定できないパラメ
ータを、発電要素の物理モデルと入力プロセス状態量と
に基づき演算により推定し、得られた各発電要素のパラ
メータを用いて、発電プラント全体のコストを最小にす
る各発電要素に対する最適な負荷配分を決定する。この
決定された最適な負荷配分と入力プロセス状態量とを用
いて、多変数モデル予測制御を行なうことにより、各発
電要素を最適に制御する。本発明においては、時々刻々
と変化する蒸気、電気需要に応じ運転員の能力に左右さ
れることなく、発電プラント全体としての運転効率の向
上が可能となる。

【０００７】

【実施例】以下、本発明の一実施例について図面を参照
して説明する。図１は本発明の一実施例を示す最適化制
御装置の構成図である。図において、最適化制御装置１
００は、機能上大別して最適化システム１１０及び多変
数制御システム１２０からなる。最適化システム１１０
は、データ入力部１１１、パラメータ推定演算部１１
２、物理モデル１１３、最適化演算部１１４、最適化演
算結果出力部１１５により構成されている。多変数制御
システム１２０は、データ入力部１２１、最適化演算結
果入力部１２２、多変数制御演算部１２３、制御演算結
果出力部１２４により構成されている。

【０００８】データ入力部１１１は、発電プラントの状
態量をリアルタイムに入力するためのもので、入力され
た状態量はパラメータ推定演算部１１２及び最適化演算
部１１４に供給される。パラメータ推定演算部１１２
は、運転条件の変化により変化する発電プラントの物理
モデル１１３のパラメータを推定するもので推定された
パラメータは最適化演算部１１４に供給される。物理モ
デル１１３は、発電プラントのボイラ、タービン等の物
質収支、熱収支等から構成されるモデルでありパラメー
タ推定演算部１１２及び最適化演算部１１４にて利用さ
れる。

【０００９】最適化演算部１１４は、データ入力部１１
１及びパラメータ推定演算部１１２より供給されるデー
タと発電プラントの物理モデル１１３をもとに、発電プ
ラント全体の運転コストが最小となる最適化演算を実施
するものであり、最適化演算結果は最適化演算結果出力
部１１５に供給される。最適化演算結果出力部１１５
は、多変数制御システム１２０における最適化演算結果
入力部１２２に供給される。

【００１０】データ入力部１２１は、発電プラントの状
態量をリアルタイムに入力するためのもので、入力され
た状態量は多変数制御演算部１２３に供給される。最適
化演算結果入力部１２２は、最適化演算結果出力部１１
５よりデータを受け取り、多変数制御演算部１２３に供
給するためのものである。多変数制御演算部１２３は、
データ入力部１２１及び最適化演算結果入力部１２２よ
り得たデータを用い、多変数制御演算を行い、演算結果
を制御演算結果出力部１２４に供給するものである。制
御演算結果出力部１２４は、制御演算結果を下位制御装
置２００に供給するためのものである。

【００１１】図２は、図５で説明した石油化学工場等で
採用されている発電プラントの一例に図１の最適化制御
装置を適用した場合の実施例を示す図である。以下、図
１及び図２を参照して詳細に説明する。図２に示すよう
な構成によると、最適化システム１１０においては、蒸
気、電気の需給バランスを満たし、しかも発電プラント
全体として経済的に最適なボイラ１，２の負荷配分及び
タービン３，４，５の負荷配分が算出される。

【００１２】すなわち、図１の最適化システム１１０に
おける最適化演算部１１４では、数理計画法により発電
プラントの状態量、パラメータ推定演算部１１２で推定
されたパラメータ、及び発電プラントの物理モデル１１
３をもとに発電プラント全体の運転コストを最小とする
ための最適化演算を行う。一般的に発電プラント全体の
物理モデルは大規模な非線形方程式となるため、大規模
非線形問題を比較的効率的に解くことが知られている逐
次二次計画法を適用したが、該数理計画法以外の手法で
も適用可能である。

【００１３】発電プラントの物理モデル１１３は、ボイ
ラ、タービン等の各構成機器について物質収支、熱収支
等からなる方程式で構成されている。例えばボイラの物
理モデルの場合は、熱交換器の配置や数に応じて異なる
モデルが用意される。しかし、熱交換器自体は汎用的な
モデルで表現し、能力の相違はパラメータで与える。ま
たタービンの場合は、タービンの段数によってモデルは
異なるが、一つの段は汎用的なモデルによって共通化さ
れている。一例として、タービンの物理モデルについて
説明する。

【００１４】タービンの各段に関する物理モデルは、タ
ービン通過蒸気が過熱状態の時、次に示す式（１）〜
（７）で記述される。 Ｐｏｗ＝η＊Ｆ＊（Ｈｉｎ−Ｈｉｓｅｎｔ）＋β （１） η ＝ｆ１（Ｆ） （２） Ｈｉｓｅｎｔ＝ｆ２（Ｐｏｕｔ．Ｔｉｓｅｎｔ） （３） Ｓｉｎ＝Ｓｏｕｔ （４） Ｓｉｎ＝ｆ３（Ｐｉｎ，Ｔｉｎ） （５） Ｓｏｕｔ＝ｆ３（Ｐｏｕｔ，Ｔｉｓｅｎｔ） （６） Ｈｉｎ＝ｆ４（Ｐｉｎ，Ｔｉｎ） （７）

【００１５】ここで、Ｐｏｗはタービン発電機出力、η
はタービン効率、Ｆはタービン段を通過する蒸気流量、
Ｈｉｎはタービン段入口蒸気のエンタルピ、Ｈｉｓｅｎ
ｔはタービン段出口において等エントロピ変化したと仮
定した時のエンタルピ、βは発電機損失パラメータ、Ｐ
ｉｎ、Ｐｏｕｔはそれぞれタービン段入口、出口蒸気の
圧力、Ｔｉｎはタービン段入口蒸気の温度、Ｔｉｓｅｎ
ｔはタービン段出口において等エントロピ変化したと仮
定した時の蒸気温度、Ｓｉｎ、Ｓｏｕｔはそれぞれター
ビン入口、出口における蒸気のエントロピを表す。

【００１６】タービン効率ηは、タービンの設計データ
あるいはプロセス状態量の実績データを解析することに
よって得られ、蒸気流量Ｆの非線形関数で表現される。
エンタルピＨｉｎ、Ｈｉｓｅｎｔ及びエントロピＳｉ
ｎ、Ｓｏｕｔはいずれも蒸気の温度、圧力の非線形の関
数で表現される。発電機損失パラメータβはタービン入
口、出口における蒸気のエンタルピ落差を電気に変換す
る際の効率を示すパラメータである。なお、ｆ１、ｆ
２、ｆ３、ｆ４はいずれも非線形関数を表す。

【００１７】パラメータ推定演算部１１２においては、
ボイラやタービンごとにリアルタイムに入力されるプロ
セスの状態量と物理モデル１１３をもとにパラメータを
推定する。該パラメータは、例えばプロセス状態量をも
とに算出される式（１）で示されるタービン発電機の損
失パラメータβやプラントの運転状況により変化する熱
交換器の汚れ係数などである。

【００１８】このように推定された各ボイラやタービン
に関するパラメータは最適化演算部１１４にて最適化演
算に使用されるため、最適化演算結果はリアルタイムに
変化するプロセス状態を反映した結果となっている。最
適化演算結果出力部１１５は、最適化演算部１１４にて
算出された最適なボイラ１，２の負荷及びタービン３，
４，５の負荷を多変数制御システム１２０に出力する。

【００１９】一方、多変数制御システム１２０は、発電
プラントの動特性に基づき最適化システム１１０にて算
出された経済的に最適な運転条件を短時間に実現すると
ともに、蒸気、電気の需給バランスの変動を短時間に吸
収するように制御する。すなわち、最適化演算結果入力
部１２２にて取り込まれる最適なボイラ１，２の負荷及
びタービン３，４，５の負荷は多変数制御演算部１２３
に取り込まれ、制御演算の設定値として設定される。多
変数制御演算部１２３では、データ入力部１２１より供
給されるプロセス状態量を制御変数、外乱変数、操作変
数として取り込み、最適化演算結果入力部１２２に取り
込まれる最適化演算結果を設定値として多変数制御演算
を行う。

【００２０】発電プラントにおけるボイラ、タービン等
の制御は相互干渉性があり、常に電気、蒸気の需要変動
による外乱の影響を受けることを考慮し、相互干渉のあ
る多変数システムにおける設定値変更、外乱除去の性能
に優れている多変数モデル予測制御を適用した。多変数
制御演算結果は、多変数制御演算結果出力部１２４を介
して、下位制御装置２００のＰＩＤ制御装置１４，１
６，１７の設定値として供給される。

【００２１】図３に、前述した式（１）〜（７）で示さ
れるタービンの物理モデルを用いてパラメータ推定を行
なう場合の、パラメータ推定演算処理の概略フローを示
す。図３において、パラメータ推定演算部１１２は、タ
ービンのパラメータ推定を行なう際、物理モデル１１３
からパラメータ推定に必要な式（１）〜（７）等を取り
込むとともに、データ入力部１１１から、リアルタイム
で変化するプロセスデータＦ，Ｐｉｎ，Ｐｏｕｔ，Ｔｉ
ｎ，Ｐｏｗ等を入力する。

【００２２】パラメータ推定演算部１１２は、入力した
プロセスデータについて式（１）〜（７）等からなる連
立方程式を解いてパラメータ値を決定する。パラメータ
推定演算部１１２は、発電プラント内のボイラ、タービ
ン等の各発電要素についてパラメータ推定を行ない、そ
れぞれの推定パラメータ値を最適化演算部１１４に出力
する。最適化演算部１１４は、パラメータ推定演算部１
１２が推定したパラメータ値とデータ入力部１１１から
入力したプロセスデータ、及び物理モデル１１３から取
り込んだ発電プラントの物理モデルを用いて、発電プラ
ント全体としての運転コストが最小になるような最適化
演算を行なう。最適化演算に用いられるコスト計算式
は、ボイラ、タービン、ボイラの給水予熱系、蒸気配管
系などの、発電プラント内の各部のモデルを等式制約と
して表現し、運転コストを最小化するような評価関数を
用いた最適化問題に定式化して解く。例えば、ある発電
プラントの場合には、１２，０００もの膨大な等式制約
を組み合わせたものとなる。

【００２３】図１に示される最適化システム１１０及び
多変数制御システム１２０の各機能は、それぞれプログ
ラムによって実現されるが、市場で入手可能な一部の汎
用パッケージソフトを利用して構築することができる。
以上のように、本発明により最適化システム１１０及び
多変数制御システム１２０からなる最適化制御装置を構
成することにより、制御性能が向上し、しかも発電プラ
ント全体としての最適化を図ることができた。図４は、
その一例として従来の手動調整の途中で本発明による最
適化制御を適用した場合の受電量の推移を示したもので
あるが、最適化制御投入後は格段に制御性能が向上して
いることが見られる。以上、本発明を図２に示す発電プ
ラントの一例に適用した場合について説明したが、本発
明は他の任意の発電プラントにも適用できることは勿論
である。

【００２４】

【発明の効果】以上のように本発明によれば、時々刻々
と変化する蒸気、電気需要に応じ運転員の能力に左右さ
れることなく、発電プラント全体としての運転効率の向
上を図ることができる。

【図面の簡単な説明】

【図１】本発明の１実施例を示す最適化制御装置の構成
図である。

【図２】発電プラントに本発明の最適化制御装置を適用
した実施例システムの構成図である。

【図３】本発明によるパラメータ推定演算処理の概略フ
ローである。

【図４】本発明による最適化制御の性能を示す受電量の
推移グラフである。

【図５】発電プラントの一例の構成図である。

【符号の説明】 １００ 最適化制御装置 １１０ 最適化システム １１１ データ入力部 １１２ パラメータ推定演算部 １１３ 物理モデル １１４ 最適化演算部 １１５ 最適化演算結果出力部 １２０ 多変数制御システム １２１ データ入力部 １２２ 最適化演算結果出力部 １２３ 多変数制御演算部 １２４ 制御演算結果出力部 ２００ 下位制御装置

───────────────────────────────────────────────────── フロントページの続き (72)発明者 中川 誠司 岡山県倉敷市潮通３丁目10番地 三菱化学 株式会社水島事業所内 Ｆターム(参考） 5G066 AA03 AA05 AE09 5H004 GA14 GA16 GB04 HA02 HA03 HA14 HA16 HB02 HB03 HB14 JA03 KA71 KB02 KB04 KB06 KC06 KC08 KC28 LA05 LA15 LA18 5H590 AA02 EA16 EB25 HA06 HA12 HA15 JA02 JA12 JA13 JA14

Claims (4)

【特許請求の範囲】
1. 【請求項１】 並列化された複数の発電要素を有する発
   電プラントにおいて、該複数の発電要素を最適負荷配分
   により制御する最適化制御方法であって、 上記複数の発電要素の各々に対する最適負荷配分の演算
   に必要なパラメータのうち、直接測定できない発電要素
   の特性や状態に関わるパラメータを、発電要素の物理モ
   デルと入力プロセス状態量とに基づき演算により推定
   し、得られた各発電要素のパラメータを用いて、発電プ
   ラント全体のコストを最小にする各発電要素に対する最
   適な負荷配分を決定して各発電要素を制御することを特
   徴とする発電プラントの最適化制御方法。
2. 【請求項２】 請求項１において、並列化された複数の
   発電要素には、複数のボイラと複数のタービンが含まれ
   ることを特徴とする発電プラントの最適化制御方法。
3. 【請求項３】 並列化された複数の発電要素からなる発
   電プラントの最適化制御装置であって、 上記複数の発電要素の各々に対する最適負荷配分の演算
   に必要なパラメータのうち、直接測定できない発電要素
   の特性や状態に関わるパラメータを、発電要素の物理モ
   デルと入力プロセス状態量とに基づき演算により推定
   し、得られた各発電要素のパラメータを用いて、発電プ
   ラント全体のコストを最小にする各発電要素に対する最
   適な負荷配分を決定する最適化システムと、 上記最適化システムで決定された各発電要素に対する最
   適な負荷配分と入力プロセス状態量とを用いて、多変数
   モデル予測制御により各発電要素に対する制御量を求
   め、設定を行なう多変数制御システムと、を備えている
   ことを特徴とする発電プラントの最適化制御装置。
4. 【請求項４】 請求項３において、並列化された複数の
   発電要素には、複数のボイラと複数のタービンが含まれ
   ることを特徴とする発電プラントの最適化制御装置。