JP3761337B2 - 発電プラントの最適化制御方法及び最適化制御装置 - Google Patents
発電プラントの最適化制御方法及び最適化制御装置 Download PDFInfo
- Publication number
- JP3761337B2 JP3761337B2 JP25910098A JP25910098A JP3761337B2 JP 3761337 B2 JP3761337 B2 JP 3761337B2 JP 25910098 A JP25910098 A JP 25910098A JP 25910098 A JP25910098 A JP 25910098A JP 3761337 B2 JP3761337 B2 JP 3761337B2
- Authority
- JP
- Japan
- Prior art keywords
- power generation
- optimization
- power plant
- generation elements
- optimization control
- 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.)
- Expired - Lifetime
Links
Images
Landscapes
- Control Of Eletrric Generators (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Feedback Control In General (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、電気や蒸気等のエネルギを工場等に供給する発電プラントを効率的に制御する最適化制御方法及び最適化制御装置に関する。本発明は特に、発電プラントが複数のボイラやタービンを並列運転するものである場合に、並列のボイラ間やタービン間での負荷配分を最適化するための有効な手段を提供する。
【０００２】
【従来の技術】
図５は、石油化学工場等で採用されている発電プラントの一例の構成図であり、重油や燃料ガス等を燃料とするボイラ１，２と、蒸気及び電気を発生するタービン３，４と、発電専用の復水タービン５、エネルギを電気に変換する発電機６，７，８と、ボイラ主蒸気母管２０と、工場２３に蒸気及び電気を供給する蒸気母管２１、電気母線２２とを有している。
【０００３】
この発電プラントの一部にはＰＩＤ制御装置が設けてある。すなわち、ボイラ主蒸気圧力を制御するＰＩＤ制御装置１３、ボイラ２の燃料流量を制御するＰＩＤ制御装置１４、蒸気圧力を制御するＰＩＤ制御装置１５、タービン４の蒸気流量を制御するＰＩＤ制御装置１６、発電機８の発電量を制御するＰＩＤ制御装置１７が設けてあり、これらによってエネルギの需給バランスを保つようにしている。なお、図を分りやすくするため、ボイラやタービンに対するＰＩＤ制御装置は一部のもののみが示されている。
【０００４】
また、ボイラの負荷配分は、運転員が燃料流量を制御するＰＩＤ制御装置１４の設定値を手動で調整することによってなされ、タービンの負荷配分は蒸気流量を制御するＰＩＤ制御装置１６の設定値を手動で調整することによってなされ、電力会社からの受電量の調整は発電量を制御するＰＩＤ制御装置１７の設定値を手動で調整することによってなされている。
このような手動調整の際に、プラント全体としての経済性を考慮した最適な運転指針を定周期で運転員に提示することが試みられている（例えば、ＳｔｅａｍＢａｌａｎｃｅ Ｏｐｔｉｍｉｚａｔｉｏｎ ｉｎ Ｃｈｅｍｉｃａｌ Ｐｌａｎｔ，ｉｎ Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ＩＦＡＣ ＳｙｍｐｏｓｉｕｍＡＤＣＨＥＭ '９４（１９９４））。これは、プラントの線形モデルをもとに最適計算を実施するが、その結果を基にした運転調整は運転員に委ねられていた。
【０００５】
【発明が解決しようとする課題】
しかしながら、上記制御システムでは次の３つの問題が発生した。
第一には、発電プラント全体としての経済性を考慮し、運転員が時々刻々と変化する蒸気、電気の需要変化に応じボイラ、タービンを常に最適な状態で運転すべく頻繁に調整することは精神的疲労が大きく、きわめて困難であり、ラフな運転とならざるを得なかった。
第二には、発電プラント全体の経済性を評価するモデルは、過去の実績データをもとに得られた線形モデルであったため、プラントの非線形特性を考慮できず運転条件が変化すると、真の最適運転条件から乖離することがあった。
第三には、ボイラ１，２及びタービン３，４，５は、ボイラ主蒸気母管２０、蒸気母管２１を介して密接に関連しあっていること、発電機６，７，８はタービン３，４，５にそれぞれ依存すること、また受電量と発電量が電気母線を介して関連していることから、発電プラントのそれぞれの制御システムが相互に干渉するため発電プラント全体として電気、蒸気の需給バランスの変動を短時間のうちに吸収することはきわめて困難であった。
【０００６】
【課題を解決するための手段】
上記課題を解決するため、本発明は、ボイラやタービンなどの並列化された複数の発電要素からなる発電プラントにおいて、複数の発電要素の各々に対する最適負荷配分を自動的に演算して制御を行なうようにするものである。この場合、最適負荷配分演算では、タービン発電機の損失パラメータや熱交換器の汚れ係数などのように、プラントからは直接測定できない発電要素の特性や状態に関わるパラメータが必要となる。本発明はこれらの直接測定できないパラメータを、発電要素の物理モデルと入力プロセス状態量とに基づき演算により推定し、得られた各発電要素のパラメータと物理モデルを用いて最適化演算を行い、その結果に基づいて発電プラント全体のコストを最小にする各発電要素に対する最適な負荷配分を決定する。この決定された最適な負荷配分と入力プロセス状態量とを用いて、多変数モデル予測制御を行なうことにより、各発電要素を最適に制御する。
本発明においては、時々刻々と変化する蒸気、電気需要に応じ運転員の能力に左右されることなく、発電プラント全体としての運転効率の向上が可能となる。
【０００７】
【実施例】
以下、本発明の一実施例について図面を参照して説明する。
図１は本発明の一実施例を示す最適化制御装置の構成図である。図において、最適化制御装置１００は、機能上大別して最適化システム１１０及び多変数制御システム１２０からなる。
最適化システム１１０は、データ入力部１１１、パラメータ推定演算部１１２、物理モデル１１３、最適化演算部１１４、最適化演算結果出力部１１５により構成されている。
多変数制御システム１２０は、データ入力部１２１、最適化演算結果入力部１２２、多変数制御演算部１２３、制御演算結果出力部１２４により構成されている。
【０００８】
データ入力部１１１は、発電プラントの状態量をリアルタイムに入力するためのもので、入力された状態量はパラメータ推定演算部１１２及び最適化演算部１１４に供給される。
パラメータ推定演算部１１２は、運転条件の変化により変化する発電プラントの物理モデル１１３のパラメータを推定するもので推定されたパラメータは最適化演算部１１４に供給される。物理モデル１１３は、発電プラントのボイラ、タービン等の物質収支、熱収支等から構成されるモデルでありパラメータ推定演算部１１２及び最適化演算部１１４にて利用される。
【０００９】
最適化演算部１１４は、データ入力部１１１及びパラメータ推定演算部１１２より供給されるデータと発電プラントの物理モデル１１３をもとに、発電プラント全体の運転コストが最小となる最適化演算を実施するものであり、最適化演算結果は最適化演算結果出力部１１５に供給される。
最適化演算結果出力部１１５は、多変数制御システム１２０における最適化演算結果入力部１２２に供給される。
【００１０】
データ入力部１２１は、発電プラントの状態量をリアルタイムに入力するためのもので、入力された状態量は多変数制御演算部１２３に供給される。
最適化演算結果入力部１２２は、最適化演算結果出力部１１５よりデータを受け取り、多変数制御演算部１２３に供給するためのものである。
多変数制御演算部１２３は、データ入力部１２１及び最適化演算結果入力部１２２より得たデータを用い、多変数制御演算を行い、演算結果を制御演算結果出力部１２４に供給するものである。
制御演算結果出力部１２４は、制御演算結果を下位制御装置２００に供給するためのものである。
【００１１】
図２は、図５で説明した石油化学工場等で採用されている発電プラントの一例に図１の最適化制御装置を適用した場合の実施例を示す図である。
以下、図１及び図２を参照して詳細に説明する。
図２に示すような構成によると、最適化システム１１０においては、蒸気、電気の需給バランスを満たし、しかも発電プラント全体として経済的に最適なボイラ１，２の負荷配分及びタービン３，４，５の負荷配分が算出される。
【００１２】
すなわち、図１の最適化システム１１０における最適化演算部１１４では、数理計画法により発電プラントの状態量、パラメータ推定演算部１１２で推定されたパラメータ、及び発電プラントの物理モデル１１３をもとに発電プラント全体の運転コストを最小とするための最適化演算を行う。一般的に発電プラント全体の物理モデルは大規模な非線形方程式となるため、大規模非線形問題を比較的効率的に解くことが知られている逐次二次計画法を適用したが、該数理計画法以外の手法でも適用可能である。
【００１３】
発電プラントの物理モデル１１３は、ボイラ、タービン等の各構成機器について物質収支、熱収支等からなる方程式で構成されている。例えばボイラの物理モデルの場合は、熱交換器の配置や数に応じて異なるモデルが用意される。しかし、熱交換器自体は汎用的なモデルで表現し、能力の相違はパラメータで与える。またタービンの場合は、タービンの段数によってモデルは異なるが、一つの段は汎用的なモデルによって共通化されている。一例として、タービンの物理モデルについて説明する。
【００１４】
タービンの各段に関する物理モデルは、タービン通過蒸気が過熱状態の時、次に示す式（１）〜（７）で記述される。
Ｐｏｗ＝η＊Ｆ＊（Ｈｉｎ−Ｈｉｓｅｎｔ）＋β （１）
η ＝ｆ１（Ｆ） （２）
Ｈｉｓｅｎｔ＝ｆ２（Ｐｏｕｔ．Ｔｉｓｅｎｔ） （３）
Ｓｉｎ＝Ｓｏｕｔ （４）
Ｓｉｎ＝ｆ３（Ｐｉｎ，Ｔｉｎ） （５）
Ｓｏｕｔ＝ｆ３（Ｐｏｕｔ，Ｔｉｓｅｎｔ） （６）
Ｈｉｎ＝ｆ４（Ｐｉｎ，Ｔｉｎ） （７）
【００１５】
ここで、Ｐｏｗはタービン発電機出力、ηはタービン効率、Ｆはタービン段を通過する蒸気流量、Ｈｉｎはタービン段入口蒸気のエンタルピ、Ｈｉｓｅｎｔはタービン段出口において等エントロピ変化したと仮定した時のエンタルピ、βは発電機損失パラメータ、Ｐｉｎ、Ｐｏｕｔはそれぞれタービン段入口、出口蒸気の圧力、Ｔｉｎはタービン段入口蒸気の温度、Ｔｉｓｅｎｔはタービン段出口において等エントロピ変化したと仮定した時の蒸気温度、Ｓｉｎ、Ｓｏｕｔはそれぞれタービン入口、出口における蒸気のエントロピを表す。
【００１６】
タービン効率ηは、タービンの設計データあるいはプロセス状態量の実績データを解析することによって得られ、蒸気流量Ｆの非線形関数で表現される。エンタルピＨｉｎ、Ｈｉｓｅｎｔ及びエントロピＳｉｎ、Ｓｏｕｔはいずれも蒸気の温度、圧力の非線形の関数で表現される。
発電機損失パラメータβはタービン入口、出口における蒸気のエンタルピ落差を電気に変換する際の効率を示すパラメータである。なお、ｆ１、ｆ２、ｆ３、ｆ４はいずれも非線形関数を表す。
【００１７】
パラメータ推定演算部１１２においては、ボイラやタービンごとにリアルタイムに入力されるプロセスの状態量と物理モデル１１３をもとにパラメータを推定する。該パラメータは、例えばプロセス状態量をもとに算出される式（１）で示されるタービン発電機の損失パラメータβやプラントの運転状況により変化する熱交換器の汚れ係数などである。
【００１８】
このように推定された各ボイラやタービンに関するパラメータは最適化演算部１１４にて最適化演算に使用されるため、最適化演算結果はリアルタイムに変化するプロセス状態を反映した結果となっている。
最適化演算結果出力部１１５は、最適化演算部１１４にて算出された最適なボイラ１，２の負荷及びタービン３，４，５の負荷を多変数制御システム１２０に出力する。
【００１９】
一方、多変数制御システム１２０は、発電プラントの動特性に基づき最適化システム１１０にて算出された経済的に最適な運転条件を短時間に実現するとともに、蒸気、電気の需給バランスの変動を短時間に吸収するように制御する。
すなわち、最適化演算結果入力部１２２にて取り込まれる最適なボイラ１，２の負荷及びタービン３，４，５の負荷は多変数制御演算部１２３に取り込まれ、制御演算の設定値として設定される。
多変数制御演算部１２３では、データ入力部１２１より供給されるプロセス状態量を制御変数、外乱変数、操作変数として取り込み、最適化演算結果入力部１２２に取り込まれる最適化演算結果を設定値として多変数制御演算を行う。
【００２０】
発電プラントにおけるボイラ、タービン等の制御は相互干渉性があり、常に電気、蒸気の需要変動による外乱の影響を受けることを考慮し、相互干渉のある多変数システムにおける設定値変更、外乱除去の性能に優れている多変数モデル予測制御を適用した。
多変数制御演算結果は、多変数制御演算結果出力部１２４を介して、下位制御装置２００のＰＩＤ制御装置１４，１６，１７の設定値として供給される。
【００２１】
図３に、前述した式（１）〜（７）で示されるタービンの物理モデルを用いてパラメータ推定を行なう場合の、パラメータ推定演算処理の概略フローを示す。
図３において、パラメータ推定演算部１１２は、タービンのパラメータ推定を行なう際、物理モデル１１３からパラメータ推定に必要な式（１）〜（７）等を取り込むとともに、データ入力部１１１から、リアルタイムで変化するプロセスデータＦ，Ｐｉｎ，Ｐｏｕｔ，Ｔｉｎ，Ｐｏｗ等を入力する。
【００２２】
パラメータ推定演算部１１２は、入力したプロセスデータについて式（１）〜（７）等からなる連立方程式を解いてパラメータ値を決定する。
パラメータ推定演算部１１２は、発電プラント内のボイラ、タービン等の各発電要素についてパラメータ推定を行ない、それぞれの推定パラメータ値を最適化演算部１１４に出力する。
最適化演算部１１４は、パラメータ推定演算部１１２が推定したパラメータ値とデータ入力部１１１から入力したプロセスデータ、及び物理モデル１１３から取り込んだ発電プラントの物理モデルを用いて、発電プラント全体としての運転コストが最小になるような最適化演算を行なう。最適化演算に用いられるコスト計算式は、ボイラ、タービン、ボイラの給水予熱系、蒸気配管系などの、発電プラント内の各部のモデルを等式制約として表現し、運転コストを最小化するような評価関数を用いた最適化問題に定式化して解く。例えば、ある発電プラントの場合には、１２，０００もの膨大な等式制約を組み合わせたものとなる。
【００２３】
図１に示される最適化システム１１０及び多変数制御システム１２０の各機能は、それぞれプログラムによって実現されるが、市場で入手可能な一部の汎用パッケージソフトを利用して構築することができる。
以上のように、本発明により最適化システム１１０及び多変数制御システム１２０からなる最適化制御装置を構成することにより、制御性能が向上し、しかも発電プラント全体としての最適化を図ることができた。
図４は、その一例として従来の手動調整の途中で本発明による最適化制御を適用した場合の受電量の推移を示したものであるが、最適化制御投入後は格段に制御性能が向上していることが見られる。
以上、本発明を図２に示す発電プラントの一例に適用した場合について説明したが、本発明は他の任意の発電プラントにも適用できることは勿論である。
【００２４】
【発明の効果】
以上のように本発明によれば、時々刻々と変化する蒸気、電気需要に応じ運転員の能力に左右されることなく、発電プラント全体としての運転効率の向上を図ることができる。
【図面の簡単な説明】
【図１】本発明の１実施例を示す最適化制御装置の構成図である。
【図２】発電プラントに本発明の最適化制御装置を適用した実施例システムの構成図である。
【図３】本発明によるパラメータ推定演算処理の概略フローである。
【図４】本発明による最適化制御の性能を示す受電量の推移グラフである。
【図５】発電プラントの一例の構成図である。
【符号の説明】
１００ 最適化制御装置
１１０ 最適化システム
１１１ データ入力部
１１２ パラメータ推定演算部
１１３ 物理モデル
１１４ 最適化演算部
１１５ 最適化演算結果出力部
１２０ 多変数制御システム
１２１ データ入力部
１２２ 最適化演算結果出力部
１２３ 多変数制御演算部
１２４ 制御演算結果出力部
２００ 下位制御装置
【発明の属する技術分野】
本発明は、電気や蒸気等のエネルギを工場等に供給する発電プラントを効率的に制御する最適化制御方法及び最適化制御装置に関する。本発明は特に、発電プラントが複数のボイラやタービンを並列運転するものである場合に、並列のボイラ間やタービン間での負荷配分を最適化するための有効な手段を提供する。
【０００２】
【従来の技術】
図５は、石油化学工場等で採用されている発電プラントの一例の構成図であり、重油や燃料ガス等を燃料とするボイラ１，２と、蒸気及び電気を発生するタービン３，４と、発電専用の復水タービン５、エネルギを電気に変換する発電機６，７，８と、ボイラ主蒸気母管２０と、工場２３に蒸気及び電気を供給する蒸気母管２１、電気母線２２とを有している。
【０００３】
この発電プラントの一部にはＰＩＤ制御装置が設けてある。すなわち、ボイラ主蒸気圧力を制御するＰＩＤ制御装置１３、ボイラ２の燃料流量を制御するＰＩＤ制御装置１４、蒸気圧力を制御するＰＩＤ制御装置１５、タービン４の蒸気流量を制御するＰＩＤ制御装置１６、発電機８の発電量を制御するＰＩＤ制御装置１７が設けてあり、これらによってエネルギの需給バランスを保つようにしている。なお、図を分りやすくするため、ボイラやタービンに対するＰＩＤ制御装置は一部のもののみが示されている。
【０００４】
また、ボイラの負荷配分は、運転員が燃料流量を制御するＰＩＤ制御装置１４の設定値を手動で調整することによってなされ、タービンの負荷配分は蒸気流量を制御するＰＩＤ制御装置１６の設定値を手動で調整することによってなされ、電力会社からの受電量の調整は発電量を制御するＰＩＤ制御装置１７の設定値を手動で調整することによってなされている。
このような手動調整の際に、プラント全体としての経済性を考慮した最適な運転指針を定周期で運転員に提示することが試みられている（例えば、ＳｔｅａｍＢａｌａｎｃｅ Ｏｐｔｉｍｉｚａｔｉｏｎ ｉｎ Ｃｈｅｍｉｃａｌ Ｐｌａｎｔ，ｉｎ Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ＩＦＡＣ ＳｙｍｐｏｓｉｕｍＡＤＣＨＥＭ '９４（１９９４））。これは、プラントの線形モデルをもとに最適計算を実施するが、その結果を基にした運転調整は運転員に委ねられていた。
【０００５】
【発明が解決しようとする課題】
しかしながら、上記制御システムでは次の３つの問題が発生した。
第一には、発電プラント全体としての経済性を考慮し、運転員が時々刻々と変化する蒸気、電気の需要変化に応じボイラ、タービンを常に最適な状態で運転すべく頻繁に調整することは精神的疲労が大きく、きわめて困難であり、ラフな運転とならざるを得なかった。
第二には、発電プラント全体の経済性を評価するモデルは、過去の実績データをもとに得られた線形モデルであったため、プラントの非線形特性を考慮できず運転条件が変化すると、真の最適運転条件から乖離することがあった。
第三には、ボイラ１，２及びタービン３，４，５は、ボイラ主蒸気母管２０、蒸気母管２１を介して密接に関連しあっていること、発電機６，７，８はタービン３，４，５にそれぞれ依存すること、また受電量と発電量が電気母線を介して関連していることから、発電プラントのそれぞれの制御システムが相互に干渉するため発電プラント全体として電気、蒸気の需給バランスの変動を短時間のうちに吸収することはきわめて困難であった。
【０００６】
【課題を解決するための手段】
上記課題を解決するため、本発明は、ボイラやタービンなどの並列化された複数の発電要素からなる発電プラントにおいて、複数の発電要素の各々に対する最適負荷配分を自動的に演算して制御を行なうようにするものである。この場合、最適負荷配分演算では、タービン発電機の損失パラメータや熱交換器の汚れ係数などのように、プラントからは直接測定できない発電要素の特性や状態に関わるパラメータが必要となる。本発明はこれらの直接測定できないパラメータを、発電要素の物理モデルと入力プロセス状態量とに基づき演算により推定し、得られた各発電要素のパラメータと物理モデルを用いて最適化演算を行い、その結果に基づいて発電プラント全体のコストを最小にする各発電要素に対する最適な負荷配分を決定する。この決定された最適な負荷配分と入力プロセス状態量とを用いて、多変数モデル予測制御を行なうことにより、各発電要素を最適に制御する。
本発明においては、時々刻々と変化する蒸気、電気需要に応じ運転員の能力に左右されることなく、発電プラント全体としての運転効率の向上が可能となる。
【０００７】
【実施例】
以下、本発明の一実施例について図面を参照して説明する。
図１は本発明の一実施例を示す最適化制御装置の構成図である。図において、最適化制御装置１００は、機能上大別して最適化システム１１０及び多変数制御システム１２０からなる。
最適化システム１１０は、データ入力部１１１、パラメータ推定演算部１１２、物理モデル１１３、最適化演算部１１４、最適化演算結果出力部１１５により構成されている。
多変数制御システム１２０は、データ入力部１２１、最適化演算結果入力部１２２、多変数制御演算部１２３、制御演算結果出力部１２４により構成されている。
【０００８】
データ入力部１１１は、発電プラントの状態量をリアルタイムに入力するためのもので、入力された状態量はパラメータ推定演算部１１２及び最適化演算部１１４に供給される。
パラメータ推定演算部１１２は、運転条件の変化により変化する発電プラントの物理モデル１１３のパラメータを推定するもので推定されたパラメータは最適化演算部１１４に供給される。物理モデル１１３は、発電プラントのボイラ、タービン等の物質収支、熱収支等から構成されるモデルでありパラメータ推定演算部１１２及び最適化演算部１１４にて利用される。
【０００９】
最適化演算部１１４は、データ入力部１１１及びパラメータ推定演算部１１２より供給されるデータと発電プラントの物理モデル１１３をもとに、発電プラント全体の運転コストが最小となる最適化演算を実施するものであり、最適化演算結果は最適化演算結果出力部１１５に供給される。
最適化演算結果出力部１１５は、多変数制御システム１２０における最適化演算結果入力部１２２に供給される。
【００１０】
データ入力部１２１は、発電プラントの状態量をリアルタイムに入力するためのもので、入力された状態量は多変数制御演算部１２３に供給される。
最適化演算結果入力部１２２は、最適化演算結果出力部１１５よりデータを受け取り、多変数制御演算部１２３に供給するためのものである。
多変数制御演算部１２３は、データ入力部１２１及び最適化演算結果入力部１２２より得たデータを用い、多変数制御演算を行い、演算結果を制御演算結果出力部１２４に供給するものである。
制御演算結果出力部１２４は、制御演算結果を下位制御装置２００に供給するためのものである。
【００１１】
図２は、図５で説明した石油化学工場等で採用されている発電プラントの一例に図１の最適化制御装置を適用した場合の実施例を示す図である。
以下、図１及び図２を参照して詳細に説明する。
図２に示すような構成によると、最適化システム１１０においては、蒸気、電気の需給バランスを満たし、しかも発電プラント全体として経済的に最適なボイラ１，２の負荷配分及びタービン３，４，５の負荷配分が算出される。
【００１２】
すなわち、図１の最適化システム１１０における最適化演算部１１４では、数理計画法により発電プラントの状態量、パラメータ推定演算部１１２で推定されたパラメータ、及び発電プラントの物理モデル１１３をもとに発電プラント全体の運転コストを最小とするための最適化演算を行う。一般的に発電プラント全体の物理モデルは大規模な非線形方程式となるため、大規模非線形問題を比較的効率的に解くことが知られている逐次二次計画法を適用したが、該数理計画法以外の手法でも適用可能である。
【００１３】
発電プラントの物理モデル１１３は、ボイラ、タービン等の各構成機器について物質収支、熱収支等からなる方程式で構成されている。例えばボイラの物理モデルの場合は、熱交換器の配置や数に応じて異なるモデルが用意される。しかし、熱交換器自体は汎用的なモデルで表現し、能力の相違はパラメータで与える。またタービンの場合は、タービンの段数によってモデルは異なるが、一つの段は汎用的なモデルによって共通化されている。一例として、タービンの物理モデルについて説明する。
【００１４】
タービンの各段に関する物理モデルは、タービン通過蒸気が過熱状態の時、次に示す式（１）〜（７）で記述される。
Ｐｏｗ＝η＊Ｆ＊（Ｈｉｎ−Ｈｉｓｅｎｔ）＋β （１）
η ＝ｆ１（Ｆ） （２）
Ｈｉｓｅｎｔ＝ｆ２（Ｐｏｕｔ．Ｔｉｓｅｎｔ） （３）
Ｓｉｎ＝Ｓｏｕｔ （４）
Ｓｉｎ＝ｆ３（Ｐｉｎ，Ｔｉｎ） （５）
Ｓｏｕｔ＝ｆ３（Ｐｏｕｔ，Ｔｉｓｅｎｔ） （６）
Ｈｉｎ＝ｆ４（Ｐｉｎ，Ｔｉｎ） （７）
【００１５】
ここで、Ｐｏｗはタービン発電機出力、ηはタービン効率、Ｆはタービン段を通過する蒸気流量、Ｈｉｎはタービン段入口蒸気のエンタルピ、Ｈｉｓｅｎｔはタービン段出口において等エントロピ変化したと仮定した時のエンタルピ、βは発電機損失パラメータ、Ｐｉｎ、Ｐｏｕｔはそれぞれタービン段入口、出口蒸気の圧力、Ｔｉｎはタービン段入口蒸気の温度、Ｔｉｓｅｎｔはタービン段出口において等エントロピ変化したと仮定した時の蒸気温度、Ｓｉｎ、Ｓｏｕｔはそれぞれタービン入口、出口における蒸気のエントロピを表す。
【００１６】
タービン効率ηは、タービンの設計データあるいはプロセス状態量の実績データを解析することによって得られ、蒸気流量Ｆの非線形関数で表現される。エンタルピＨｉｎ、Ｈｉｓｅｎｔ及びエントロピＳｉｎ、Ｓｏｕｔはいずれも蒸気の温度、圧力の非線形の関数で表現される。
発電機損失パラメータβはタービン入口、出口における蒸気のエンタルピ落差を電気に変換する際の効率を示すパラメータである。なお、ｆ１、ｆ２、ｆ３、ｆ４はいずれも非線形関数を表す。
【００１７】
パラメータ推定演算部１１２においては、ボイラやタービンごとにリアルタイムに入力されるプロセスの状態量と物理モデル１１３をもとにパラメータを推定する。該パラメータは、例えばプロセス状態量をもとに算出される式（１）で示されるタービン発電機の損失パラメータβやプラントの運転状況により変化する熱交換器の汚れ係数などである。
【００１８】
このように推定された各ボイラやタービンに関するパラメータは最適化演算部１１４にて最適化演算に使用されるため、最適化演算結果はリアルタイムに変化するプロセス状態を反映した結果となっている。
最適化演算結果出力部１１５は、最適化演算部１１４にて算出された最適なボイラ１，２の負荷及びタービン３，４，５の負荷を多変数制御システム１２０に出力する。
【００１９】
一方、多変数制御システム１２０は、発電プラントの動特性に基づき最適化システム１１０にて算出された経済的に最適な運転条件を短時間に実現するとともに、蒸気、電気の需給バランスの変動を短時間に吸収するように制御する。
すなわち、最適化演算結果入力部１２２にて取り込まれる最適なボイラ１，２の負荷及びタービン３，４，５の負荷は多変数制御演算部１２３に取り込まれ、制御演算の設定値として設定される。
多変数制御演算部１２３では、データ入力部１２１より供給されるプロセス状態量を制御変数、外乱変数、操作変数として取り込み、最適化演算結果入力部１２２に取り込まれる最適化演算結果を設定値として多変数制御演算を行う。
【００２０】
発電プラントにおけるボイラ、タービン等の制御は相互干渉性があり、常に電気、蒸気の需要変動による外乱の影響を受けることを考慮し、相互干渉のある多変数システムにおける設定値変更、外乱除去の性能に優れている多変数モデル予測制御を適用した。
多変数制御演算結果は、多変数制御演算結果出力部１２４を介して、下位制御装置２００のＰＩＤ制御装置１４，１６，１７の設定値として供給される。
【００２１】
図３に、前述した式（１）〜（７）で示されるタービンの物理モデルを用いてパラメータ推定を行なう場合の、パラメータ推定演算処理の概略フローを示す。
図３において、パラメータ推定演算部１１２は、タービンのパラメータ推定を行なう際、物理モデル１１３からパラメータ推定に必要な式（１）〜（７）等を取り込むとともに、データ入力部１１１から、リアルタイムで変化するプロセスデータＦ，Ｐｉｎ，Ｐｏｕｔ，Ｔｉｎ，Ｐｏｗ等を入力する。
【００２２】
パラメータ推定演算部１１２は、入力したプロセスデータについて式（１）〜（７）等からなる連立方程式を解いてパラメータ値を決定する。
パラメータ推定演算部１１２は、発電プラント内のボイラ、タービン等の各発電要素についてパラメータ推定を行ない、それぞれの推定パラメータ値を最適化演算部１１４に出力する。
最適化演算部１１４は、パラメータ推定演算部１１２が推定したパラメータ値とデータ入力部１１１から入力したプロセスデータ、及び物理モデル１１３から取り込んだ発電プラントの物理モデルを用いて、発電プラント全体としての運転コストが最小になるような最適化演算を行なう。最適化演算に用いられるコスト計算式は、ボイラ、タービン、ボイラの給水予熱系、蒸気配管系などの、発電プラント内の各部のモデルを等式制約として表現し、運転コストを最小化するような評価関数を用いた最適化問題に定式化して解く。例えば、ある発電プラントの場合には、１２，０００もの膨大な等式制約を組み合わせたものとなる。
【００２３】
図１に示される最適化システム１１０及び多変数制御システム１２０の各機能は、それぞれプログラムによって実現されるが、市場で入手可能な一部の汎用パッケージソフトを利用して構築することができる。
以上のように、本発明により最適化システム１１０及び多変数制御システム１２０からなる最適化制御装置を構成することにより、制御性能が向上し、しかも発電プラント全体としての最適化を図ることができた。
図４は、その一例として従来の手動調整の途中で本発明による最適化制御を適用した場合の受電量の推移を示したものであるが、最適化制御投入後は格段に制御性能が向上していることが見られる。
以上、本発明を図２に示す発電プラントの一例に適用した場合について説明したが、本発明は他の任意の発電プラントにも適用できることは勿論である。
【００２４】
【発明の効果】
以上のように本発明によれば、時々刻々と変化する蒸気、電気需要に応じ運転員の能力に左右されることなく、発電プラント全体としての運転効率の向上を図ることができる。
【図面の簡単な説明】
【図１】本発明の１実施例を示す最適化制御装置の構成図である。
【図２】発電プラントに本発明の最適化制御装置を適用した実施例システムの構成図である。
【図３】本発明によるパラメータ推定演算処理の概略フローである。
【図４】本発明による最適化制御の性能を示す受電量の推移グラフである。
【図５】発電プラントの一例の構成図である。
【符号の説明】
１００ 最適化制御装置
１１０ 最適化システム
１１１ データ入力部
１１２ パラメータ推定演算部
１１３ 物理モデル
１１４ 最適化演算部
１１５ 最適化演算結果出力部
１２０ 多変数制御システム
１２１ データ入力部
１２２ 最適化演算結果出力部
１２３ 多変数制御演算部
１２４ 制御演算結果出力部
２００ 下位制御装置
Claims (4)
- 並列化された複数の発電要素を有する発電プラントにおいて、該複数の発電要素を最適負荷配分により制御する最適化制御方法であって、
上記複数の発電要素の各々に対する最適負荷配分の演算に必要なパラメータのうち、直接測定できない発電要素の特性や状態に関わるパラメータを、発電要素の物理モデルと入力プロセス状態量とに基づき演算により推定し、得られた各発電要素のパラメータと物理モデルを用いて最適化演算を行い、その結果に基づいて発電プラント全体のコストを最小にする各発電要素に対する最適な負荷配分を決定して各発電要素を制御することを特徴とする発電プラントの最適化制御方法。 - 請求項１において、並列化された複数の発電要素には、複数のボイラと複数のタービンが含まれることを特徴とする発電プラントの最適化制御方法。
- 並列化された複数の発電要素からなる発電プラントの最適化制御装置であって、
上記複数の発電要素の各々に対する最適負荷配分の演算に必要なパラメータのうち、直接測定できない発電要素の特性や状態に関わるパラメータを、発電要素の物理モデルと入力プロセス状態量とに基づき演算により推定し、得られた各発電要素のパラメータと物理モデルを用いて最適化演算を行い、その結果に基づいて発電プラント全体のコストを最小にする各発電要素に対する最適な負荷配分を決定する最適化システムと、
上記最適化システムで決定された各発電要素に対する最適な負荷配分と入力プロセス状態量とを用いて、多変数モデル予測制御により各発電要素に対する制御量を求め、設定を行なう多変数制御システムと、
を備えていることを特徴とする発電プラントの最適化制御方法。 - 請求項３において、並列化された複数の発電要素には、複数のボイラと複数のタービンが含まれることを特徴とする発電プラントの最適化制御装置。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP25910098A JP3761337B2 (ja) | 1998-08-31 | 1998-08-31 | 発電プラントの最適化制御方法及び最適化制御装置 |
| PCT/JP1999/004703 WO2004084371A1 (ja) | 1998-08-31 | 1999-08-31 | 発電プラントの最適化制御方法及び最適化制御装置 |
| US09/763,861 US6681155B1 (en) | 1998-08-31 | 1999-08-31 | Optimizing control method and optimizing control system for power plant |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP25910098A JP3761337B2 (ja) | 1998-08-31 | 1998-08-31 | 発電プラントの最適化制御方法及び最適化制御装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2000078749A JP2000078749A (ja) | 2000-03-14 |
| JP3761337B2 true JP3761337B2 (ja) | 2006-03-29 |
Family
ID=17329331
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP25910098A Expired - Lifetime JP3761337B2 (ja) | 1998-08-31 | 1998-08-31 | 発電プラントの最適化制御方法及び最適化制御装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3761337B2 (ja) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20030095881A (ko) * | 2002-06-15 | 2003-12-24 | 삼성정밀화학 주식회사 | Ｃｏｄ 예측제어를 통한 ｅｃｈ 공정 폐수처리 최적화방법 |
| JP4043879B2 (ja) | 2002-07-24 | 2008-02-06 | 三菱電機株式会社 | プラント最適運転制御装置 |
| JP4664842B2 (ja) * | 2006-03-20 | 2011-04-06 | 株式会社東芝 | エネルギープラントの最適運用システムと方法、およびプログラム |
| JP4932294B2 (ja) * | 2006-03-23 | 2012-05-16 | 新日本製鐵株式会社 | 製造仕様決定支援システム、製造仕様決定支援方法、コンピュータプログラム、及びコンピュータ読み取り可能な記録媒体 |
| JP5077831B2 (ja) * | 2008-09-25 | 2012-11-21 | 横河電機株式会社 | プラント制御システムおよびプラント制御方法 |
| JP5649596B2 (ja) * | 2012-01-27 | 2015-01-07 | 株式会社日立製作所 | プロセスデータ推定方法 |
| US9477214B2 (en) * | 2013-03-07 | 2016-10-25 | General Electric Company | Plant control systems and methods |
| CN103293961B (zh) * | 2013-06-07 | 2015-03-25 | 天津大学 | 一种基于需求响应控制的能效电厂建模和集成方法 |
| US9476584B2 (en) | 2013-12-12 | 2016-10-25 | General Electric Company | Controlling boiler drum level |
| CN111221253B (zh) * | 2020-03-11 | 2021-01-26 | 国网江苏省电力有限公司宿迁供电分公司 | 适用于三相并网逆变器的鲁棒模型预测控制方法 |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0814241B2 (ja) * | 1988-08-10 | 1996-02-14 | 三菱重工業株式会社 | タービン制御装置 |
| JPH0286905A (ja) * | 1988-09-22 | 1990-03-27 | Mitsubishi Heavy Ind Ltd | タービン制御装置 |
| JP2561265Y2 (ja) * | 1989-07-03 | 1998-01-28 | 三菱重工業株式会社 | 自家発電プラント運転支援装置 |
| JPH04266101A (ja) * | 1991-02-20 | 1992-09-22 | Toshiba Corp | 多変数モデル予測制御装置 |
| JP3061450B2 (ja) * | 1991-07-30 | 2000-07-10 | 株式会社東芝 | モデル予測制御装置 |
| JP2823371B2 (ja) * | 1991-02-26 | 1998-11-11 | 株式会社東芝 | 電力最適化管理装置 |
| JPH06236202A (ja) * | 1993-02-10 | 1994-08-23 | Hitachi Ltd | プラントの運転方法及び装置 |
| JPH07238803A (ja) * | 1994-02-28 | 1995-09-12 | Mitsubishi Heavy Ind Ltd | 発電プラント群に適用される運転制御装置 |
| JPH07281727A (ja) * | 1994-04-05 | 1995-10-27 | Mitsubishi Electric Corp | 熱交換器内部状態のシミュレーション方法 |
| JP3215279B2 (ja) * | 1995-03-15 | 2001-10-02 | 三菱電機株式会社 | エネルギー設備の設備運用方法およびその装置 |
| JPH09179604A (ja) * | 1995-09-13 | 1997-07-11 | Toshiba Corp | プラントの運転制御システム及び方法 |
-
1998
- 1998-08-31 JP JP25910098A patent/JP3761337B2/ja not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JP2000078749A (ja) | 2000-03-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6681155B1 (en) | Optimizing control method and optimizing control system for power plant | |
| JP3758862B2 (ja) | 発電プラントの最適運転制御方法及び最適運転制御装置 | |
| Velasco-Garcia et al. | Utility systems operation: Optimisation-based decision making | |
| CA2595739C (en) | Steam temperature control using integrated function block | |
| CA2914954C (en) | Model-based combined cycle power plant load control | |
| US9841185B2 (en) | Steam temperature control using model-based temperature balancing | |
| JP3761337B2 (ja) | 発電プラントの最適化制御方法及び最適化制御装置 | |
| CN102374518B (zh) | 使用动态矩阵控制的蒸汽温度控制 | |
| US10240775B2 (en) | Multi-objective steam temperature control | |
| CA2943023C (en) | Feedforward control with intermittent re-initialization based on estimated state information | |
| JP2012117523A (ja) | 蒸気タービンの装荷方法及び装荷システム | |
| US4577280A (en) | Control system for fluid flow distribution | |
| Bindlish | Power scheduling and real-time optimization of industrial cogeneration plants | |
| Zotică et al. | Optimal operation and control of heat to power cycles: A new perspective from a systematic plantwide control approach | |
| Kocaarslan et al. | An adaptive control application in a large thermal combined power plant | |
| Kordestani et al. | Predictive control of large steam turbines | |
| JP2000097001A (ja) | タービンの最適運転制御方法及び最適運転制御装置 | |
| CZ290525B6 (cs) | Způsob diagnózy a prognózy provozního chování zařízení turbíny a zařízení k provádění tohoto způsobu | |
| Kaya et al. | Methods of energy efficient control and optimization for combined-cycle cogeneration | |
| Śladewski et al. | Close loop optimisation of large CHP based on approximation model | |
| Pablos de la Fuente et al. | Modeling on-site combined heat and power systems coupled to main process operation | |
| El-Guindy et al. | Centralized multivariable feedback control of steam drums in combined cycle power plants | |
| Hristov | MODELING OF THE ENERGY EFFICIENCY OF CHP STEAM TURBINE INSTALLATIONS | |
| Majanne | Simulation as a tool for process design and disturbance analysis in industrial power plants | |
| Srovnal et al. | Distributed control systems of pipeline networks using multi-agent technologies |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20040120 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20040319 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050426 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050624 |
|
| 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: 20051213 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20060110 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090120 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100120 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110120 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110120 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120120 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130120 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130120 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140120 Year of fee payment: 8 |
|
| EXPY | Cancellation because of completion of term |