JPH0573886B2 - - Google Patents

Info

Publication number
JPH0573886B2
JPH0573886B2 JP23902086A JP23902086A JPH0573886B2 JP H0573886 B2 JPH0573886 B2 JP H0573886B2 JP 23902086 A JP23902086 A JP 23902086A JP 23902086 A JP23902086 A JP 23902086A JP H0573886 B2 JPH0573886 B2 JP H0573886B2
Authority
JP
Japan
Prior art keywords
schedule
startup
time
plant
function
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
Application number
JP23902086A
Other languages
English (en)
Other versions
JPS6394010A (ja
Inventor
Hiroshi Matsumoto
Yasuo Morooka
Daizo Iba
Seiitsu Nikawara
Isamu Sano
Kazuharu Aoyanagi
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Tokyo Electric Power Co Holdings Inc
Original Assignee
Tokyo Electric Power Co Inc
Hitachi Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tokyo Electric Power Co Inc, Hitachi Ltd filed Critical Tokyo Electric Power Co Inc
Priority to JP23902086A priority Critical patent/JPS6394010A/ja
Publication of JPS6394010A publication Critical patent/JPS6394010A/ja
Publication of JPH0573886B2 publication Critical patent/JPH0573886B2/ja
Granted legal-status Critical Current

Links

Landscapes

  • Control Of Steam Boilers And Waste-Gas Boilers (AREA)
  • Feedback Control In General (AREA)

Description

【発明の詳細な説明】
〔産業上の利用分野〕 本発明は発電プラントの起動装置に係り、特に
スケジユール誤差が所定以上となつた場合に、そ
の旨を関連するプラント状態量と共に運転員に知
らせるのに好適なスケジユール誤差監視方法に関
する。 〔従来の技術〕 火力発電プラントの起動に関する従来の方法
は、起動前の停止時間や機器の温度状態に応じ
て、ボイラへの初期投入燃料量、主蒸気の昇温及
び昇圧の時間関数、タービンの昇速及び負荷上昇
の時間関数を起動スケジユールとして決定し、こ
の起動スケジユールをプラントの各系統に設けら
れた制御系で実行するという方法が採られてい
た。 この最も代表的な方法は、Electrical World、
Vol.165、No.6の論文“Thermal Stress
Influence Starting、Loading of Boilers
Turbines”で述べられている。この方法は、プ
ラントの限られた部分の初期状態によつて一義的
に起動スケジユールを決定する方法である。即
ち、ボイラ蒸気圧力、ボイラ出口蒸気温度、蒸気
タービンケーシング温度の初期値に応じて、蒸気
タービンの昇速率、初負荷量、速度保持並びに負
荷保持による蒸気タービンの暖機時間及び負荷変
化率を決定する方法である。この方法によると、
運転制限要因である蒸気タービンの熱応力を管理
する上で重要なボイラ発生蒸気の昇温特性を起動
前に予測できないため、その不確定性を起動スケ
ジユールに余裕を持たせることにより吸収してい
る。そのため、作成される起動スケジユールは必
要以上に長くなりがちであつた。 また、別の従来方法としては、USP3446224及
びUSP4228359が知られている。これらは蒸気タ
ービンに発生する熱応力をオンラインリアルタイ
ムで監視しながら蒸気タービンの急速起動を図る
ものであるが、前記従来方法と同様にボイラの起
動方法に関しては何ら言及していない。 ボイラの起動時間の短縮を目的とした従来方法
としては、特開昭59−157402が知られている。こ
の方法はボイラに発生する熱応力をオンラインリ
アルタイムで監視しながらボイラ発生蒸気の急速
昇温を図るものである。しかし、この方式は蒸気
タービンの起動に関しては何ら言及していない。 プラント全体の起動時間は、ボイラと蒸気ター
ビンの協調により短縮が可能であるが、以上述べ
た従来方法は、何れもボイラもしくは蒸気タービ
ンの片方のみに着目した急速起動方法であり、こ
の様な個別の方法を組合わせたとしてもプラント
全体の起動時間が最短となる保証は何も無い。何
故ならば、ボイラと蒸気タービンは相互干渉が極
めて強く、個々の最適化が必ずしも全体の最適化
とならないからである。 さらに、上記従来方法では、ボイラ点火直前の
プラント初期温度状態の実測値に基づいて起動ス
ケジユールを作成するため、ボイラ点火前の任意
時刻にて起動スケジユールを作成し、且つ中央給
電指令所(以下、中給と呼び)から指定された起
動完了時刻に正確に起動を完了(以下、このこと
を定刻起動と呼ぶ)できる起動スケジユールを作
成することができなかつた。 このように、従来方法では起動時間が長くなり
がちであり、それに伴つて起動時のエネルギ損失
(以下、このことを起動損失と呼ぶ)が大きくな
るといつたプラント運用上の経済性の問題と、起
動スケジユールを実行する上での安全性の問題が
あつた。 さらに、従来方法では起動前に作成した起動ス
ケジユールが実際に遂行される起動過程におい
て、発生するスケジユール誤差を監視する機能を
有していない。そのため、スケジユール誤差発生
要因が異常監視対象の状態量にまで影響が波及し
た時始めて運転員はスケジユール誤差の発生に気
がつき、対応が遅れがちであつた。そのため、ス
ケジユール誤差が大きくなりがちであり、中給か
ら指定される起動時刻を守れない場合があつた。 〔発明が解決しようとする問題点〕 本発明の目的は、前記従来方法にない起動中の
スケジユール監視機能を設けることにより、スケ
ジユール誤差発生の早期発見を可能にし、中給か
ら指定される起動時刻をまもるための早期対策を
運転員がとり易くする発電プラント起動装置を提
供することにある。 〔問題点を解決するための手段〕 上記目的は、プラントの起動特性を定量的に計
算するための動特性モデルと、この動特性モデル
で使用するボイラ点火時のプラント初期温度状態
を予測するための初期値予測手段と、最短時間で
の起動を実現するために上記動特性モデルで得ら
れる起動特性を評価しながら繰返し計算手法によ
り起動スケジユールを最適化するスケジユール最
適化手段と、上記繰返し計算用の初期スケジユー
ルを決定するための基本スケジユール作成手段
と、上記スケジユール最適化手段により作成され
た起動スケジユールと実際の起動過程で遂行され
る起動スケジユールとの誤差を監視するスケジユ
ール誤差監視手段を用いることにより、達成され
る。 〔作 用〕 上記動特性モデルは、仮定した起動スケジユー
ルに従つてプラントを起動した場合の起動特性を
定量的に計算することができる。そのため、起動
スケジユールと起動特性の関係を、実際にプラン
トを起動する前に知ることが可能であるととも
に、プラント状態値が運転制限条件を満足するか
否かを事前に確認できるため、安全性の高い起動
スケジユールを作成できる。 初期値予測手段は、熟練運転員の思想方法と似
た方法で初期値を予測するために、フアジイ推論
を適用し、プラントの標準的な冷却特性と実際の
冷却状態との定性的な差異により初期値を予測す
る方法を採つている。これによつて、高精度の初
期値予測が可能であり、その結果、動特性モデル
による高精度な起動特性の予測が可能である。 スケジユール最適化手段は、熟練運転員の思考
方法と似た方法でスケジユールを最適化するため
に、フアジイ推論を適用し、上記動特性モデルで
得られた起動特性を定性的に評価し、スケジユー
ルを修正し、再度動特性モデルを用いて起動特性
を計算し、その結果を定性的に評価し、スケジユ
ールを修正するといつた繰返し計算方法により最
適化を行なつている。このスケジユールと起動特
性の定性的因果関係に基づいてスケジユールの修
正量を決定するところが熟練運転員の思考方法と
似ている。これにより、最適スケジユール求解の
ための収束性が極めて良好となるため、詳細な計
算式を用いた大規模な動特性モデルを用いること
が可能で、起動特性を高精度で予測することがで
きる。 基本スケジユール作成手段は、前記初期値予測
手段で得られた初期値に応じて、上記スケジユー
ル最適化が用いる初期スケジユールを作成するた
めのもので、予測した初期値と標準的な初期値と
の定性的な差異により標準的なスケジユールを修
正することにより適切な初期スケジユールが作成
される。これにより、スケジユール作成に関する
熟練運転員の思考方法を真似ることができ、上記
スケジユール最適化の収束性を向上できる。 以上述べた手段により、プラントの運転制限条
件を満足し、起動所要時間を最小にし、中給から
指定された時刻通りに起動を完了できる最適起動
スケジユールを作成することが可能となる。 また、スケジユール誤差監視手段では、起動中
に遂行されるスケジユールに誤差が発生すると、
その旨を関係するプラント状態量と共に運転員に
提示することが可能となる。 〔実施例〕 以下、本発明の一実施例を説明する。 第1図は本発明を適用した発電プラント起動装
置の全体構成を示すものである。本装置は、ボイ
ラ、蒸気タービン及び発電機から成るプラント2
000において、ボイラ点火から目標負荷(一般
的には中給から指令される負荷レベル)に到達す
るまでの全起動過程における機器の操作時期及び
制御目標を決定するための起動スケジユール作成
機能1000と、作成された起動スケジユールに
従つてプラントを自動的に起動してゆくための起
動制御機能3000と、プラントの運転員600
0からのリクエストに応じてプラント起動に関す
る必要な情報を表示装置5000に表示したり、
スケジユール作成機能1000内の情報を書換え
たりするためのユーザインターフエイス4000
から構成されている。 起動スケジユール作成機能1000は、更に、
動特性モデル100、初期値予測機能200、基
本スケジユール作成機能300、スケジユール最
適化機能400、スケジユール誤差監視機能70
0、フアジイ推論機能500、知識ベース600
から成る。知識ベース600は、更に、初期値予
測ルール610、基本スケジユール作成ルール6
40、スケジユール最適化ルール670から成
る。 上記起動スケジユール作成機能1000を詳細
に説明する前に、各機能の目的について説明して
おく。 動特性モデル100は、初期値予測機能200
にて予測されたボイラ点火時のプラントの状態を
初期値210として、基本スケジユール作成機能
300あるいはスケジユール最適化機能400か
ら与えられる起動スケジユールに従つてプラント
を起動した場合の起動特性140を定量的に計算
するためのものである。 初期値予測機能200は、起動前の任意時刻に
おいて、ボイラ点火時のプラント状態を予測し、
これを動特性モデル100及び基本スケジユール
作成機能300に設定するためのものである。 基本スケジユール作成機能300は、スケジユ
ール最適化機能400における最適値求解計算で
良好な収束性を得るために初期スケジユールとし
ての基本スケジユール310を決定し、動特性モ
デル100に設定するためのものである。 スケジユール最適化機能400は、動特性モデ
ル100を用いて起動特性を予測し、その結果に
応じてスケジユールを修正し、再び動特性モデル
100に起動スケジユール410を設定し、起動
特性を予測するといつた繰返し計算手法で最適ス
ケジユール420を決定するものである。 スケジユール誤差監視機能700は、起動前に
作成された最適スケジユールを実行中にスケジユ
ールの遂行状態即ちスケジユール誤差を監視する
ものであり、誤差が所定値以上となつた場合は、
その旨を関連するプラント状態量とともにユーザ
インターフエイス4000を介して表示装置50
00に表示する。 フアジイ推論機能500は、前記、初期値予測
機能200、基本スケジユール作成機能300及
びスケジユール最適化機能400に対して作用
し、それぞれの処理において、熟練運転員の思考
方法を模擬することにより、起動特性の予測精度
の向上と最適スケジユール求解の高速化を図るも
のである。そのために、上記3機能に対して、熟
練運転員の知識としてそれぞれ初期値予測ルール
610、基本スケジユール作成ルール640、ス
ケジユール最適化ルール670を準備し、これを
知識ベースとしている。 ここで、スケジユール最適化機能400による
起動スケジユール最適化の基本的考え方を第2図
を用いて説明しておく。 第2図の中で破線はスケジユール最適化前即ち
基本スケジユールに従つてプラントを起動した場
合のタービン応力、起動パターン及び起動時間を
示している。また、実線はスケジユール最適化後
のそれぞれを示している。本図では、中給から併
入時刻(発電機を電力系統に接続する時刻)が指
定された場合を示すものであるが、起動完了時刻
を指定された場合でも本発明は基本的に同じ方法
スケジユール最適化が可能である、最適化前の起
動スケジユールに従つてプラントを起動した場
合、図のように、タービン応力は起動の前半で制
限値に対して大きな余裕があり、後半では余裕が
小さくなり部分的に制限値以上の応力が発生して
いる。この様な起動特性が前記動特性モデル10
0で予測されると、スケジユール最適化機能40
0ではスケジユール最適化ルール670を用いた
フアジイ推論機能500を動作させ、起動の前半
では点火時刻を遅らせたり、タービンの速度保持
時間及び負荷保持時間を短縮するなどして起動時
間の短縮を図つている。また、起動の後半では負
荷保持時間を延長することによりタービン応力の
緩和を図つている。このように、本発明を適用し
た起動装置を用いると、運転制限要因であるター
ビン応力を制限値以下に抑え、かつ最短時間で中
給から指定された時刻通りに起動できることがで
きる。 以下、上記に概説した起動スケジユール作成機
能1000について、各構成機能毎に詳細を説明
する。 (1) 初期値予測機能200 動特性モデル100を用いて起動特性を高精度
で予測するためには、動特性モデルで用いるプラ
ント初期値(ボイラ点火時における値で、主に温
度状態)を高精度に予測する必要がある。しか
し、現時点での計算値に基づいて冷却特性を解析
時に求めることは、停止時の操作内容や現在プラ
ントが置かれた環境などに影響されるため困難で
ある。ところが、豊富な運転経験をもつ運転員
は、現時点での温度状態が標準的な状態からどの
程度ずれているかによつて、将来のずれをかなり
的確に予測できる。本初期値予測機能200は、
この点に着目したもので、標準冷却特性で予測し
た値を、上記現時点でのずれに応じて、どの程度
修正するかをルール化した運転員の知識を利用す
るものである。以下、本初期値予測機能200に
関して具体的に説明する。 前回の解列時刻(発電機を電力系統から切り離
した時刻)tPFと中給からの併入指定時刻tPIとの
差を停止時間ΔtTSとする。この停止時間に応じて
予め規定された4つの起動モード毎に準備された
標準スケジユールの中から該当する起動モードを
決定し、起動スケジユールを選択する。ここでい
う起動スケジユールとは、ボイラ点火から起動完
了までの起動パターンを規定する次のパラメータ
で定義する。 (a) ボイラ起動所要時間(点火から通気まで) (b) 第1速度保持時間(1000rpmにおいて) (c) 第2速度保持時間(2800rpmにおいて) (d) 第3速度保持時間(3600rpmにおいて) (e) 第1負荷保持時間(初負荷レベルにおいて) (f) 第2負荷保持時間(20%負荷レベルにおい
て) (g) 第3負荷保持時間(40%負荷レベルにおい
て) また、タービン昇速率(rpm/分)及び負荷変
化率(%/分)も起動モードに対応して一義的に
決定される。上記スケジユールパラメータが決ま
ると、併入時刻からの逆算により、ボイラ点火時
刻tIGが決定される。次に、起動モード毎に準備
した標準初基値(ボイラ点火時相当)の中から、
今回選択された起動モードにより該当する初期値
を選択する。ここで使用する初期値は次のプラン
ト状態に関するものとする。 (a) ドラム温度 (b) 過熱器出口蒸気温度 (c) 再熱器出口蒸気温度 (d) 主蒸気管メタル温度 (e) 再熱蒸気管メタル温度 (f) 水壁入口内部流体温度 (g) 節炭器出口内部流体温度 (h) 節炭器入口内部流体温度 (i) 高圧タービン(HPT)第1段後メタル温
度 (j) 中圧タービン(IPT)第1段後メタル温度 (k) ドラム圧力 一方、停止中のタービン及びボイラの標準冷却
特性を用いて、上記状態量に関する現在即ち起動
スケジユール作成時の標準的な値を求める。これ
を現在標準値推定機能と呼ぶことにする。なお、
上記標準冷却特性を次式で定義する。 T=(TPFT−TA)e-TIME/TETC+TA …(1) ここで、TPFT:プラント停止時の標準温度
(℃) TA :大気温度(℃) TIME:解列後経過時間(分) TETC:冷却時定数(分) 上記ドラム圧力に関しても同様である。 また、同様にしてボイラ点火時の初期状態を予
測する。この予測を高精度に行なうために、現在
状態の実測値230(第1図参照)と(1)式による標準
値との差を考慮したフアジイ推論により上記予測
値を修正する。 以下、この修正方法を説明する。 いま、現在値偏差E(1),E(2),……E(11)は、そ
れぞれ、現時点(t0)で実測したドラム温度X1
(t0)、過熱器出口蒸気温度X2(t0),……ドラム圧
力X11(t0)と(1)式から得られたその標準値X1S
(t0),X2S(t0),……X11S(t0)との差を標準値で
正規化したものであり、次式で定義する。 E(i)=Xi(t0)−XiS(t0)/XiS(t0)…(2) 第3図は、上記現在値偏差の大きさを定性的に
評価するためのメンバーシツプ関数である。図中
のESMB(i)(i=1〜7)はメンバーシツプ関数の
形を規定する定数であり、PB,PM,PS,ZO,
NS,NM,NBは偏差Eの大きさを定性的に評
価するためにメンバーシツプ関数に与えた名称で
あり、それぞれ下記の意味をもつ。 PB:Positive Big PM:Positive Medium PS:Positive Small ZO:Zero NS:Negative Small NM:Negative Medium NB:Negative Big また、図の縦軸はメンバーシツプ値である。こ
のメンバーシツプ関数を用いて、前記11個の状態
量に関する現在値偏差E(i)(i=1〜11)を定性
的に評価する。 第4図は、現在値偏差がボイラ点火時初期に及
ぼす影響度を定性的に整理したものであり、これ
に従つて作成した初期値予測ルールの一例を第5
図に示す。本図は、過熱器出口蒸気温度の現在値
偏差E(2)と再熱器出口蒸気温度の現在値偏差E(3)
の定性的関係において、ボイラ点火時過熱器出口
蒸気温度の予測値((1)式で予測した値)の修正量
DTIG(3)を定性的に決定するためのものである。例
えば、ルールNo.32の場合は IF(E(2) is NS and E(3) is NM)THEN
(DTIG(2) is NS) という意味である。 第6図は、定性的に決定された初期値修正量を
定量的な値に変換するためのメンバーシツプ関数
である。図中のEIMB(i)(i=1〜7)はメンバー
シツプ関数の形を規定する定数である。PB,
PM,PS,ZO,NS,NM,NBは修正量の大き
さを定性的に表わすためにメンバーシツプ関数に
与えた名称であり、第5図の中で使用している名
称に対応している。また、図の縦軸はメンバーシ
ツプ値である。適用されたルールにより、初期値
修正量が定性的にどのメンバーシツプ関数に属す
るか決定される。或る初期値に対する修正量が複
数のルールによる複数のメンバーシツプで規定さ
れた場合は、各メンバーシツプ値に応じた加重平
均値をもつて実際の定量的修正量510(第1図
参照)とする。これにより各初期値を予測したこ
とになる。 (2) 基本スケジユール作成機能300 基本スケジユールとはスケジユール最適化の収
束計算の初期スケジユールであり、良い収束性を
得るには、できるだけ最適値の近傍に設定するの
が望ましい。豊富な運転経験をもつ運転員は、ボ
イラ点火時の初期状態に応じて起動スケジユール
をかなり的確に決定することができる。本基本ス
ケジユール作成機能300は、この点に着目した
もので、初期状態の標準値からのずれに応じて予
め準備した標準スケジユールのパラメータをどの
程度修正するかをルール化した運転員の知識を利
用するものである。以下、本基本スケジユール作
成機能300について具体的に説明する。 起動モードに対応して選択された前記起動スケ
ジユールは標準スケジユールの中から選ばれたも
のであり、必ずしも今回の起動条件にマツチした
ものではない。そこで、前記標準初期値と上記方
法で予測された初期値との差を考慮してスケジユ
ールを修正するのが本機能である。この修正量は
上記予測値の差に応じてフアジイ推論により決定
する。 以下、その修正量決定方法について説明する。 いま、点火時予測偏差EP(1),EP(2),……EP(11)
は、それぞれ、初期値予測機能200で予測した
ボイラ点火時のドラム温度X1(tIG),過熱器出口
蒸気温度X2(tIG),……ドラム圧力X11(tIG)とそ
の標準値X1(tIG),X2(tIG),……X11(tIG)との差
を標準値で正規化したものであり、次式で定義す
る。 EP(i)=XiS(tIG)−Xi(tIG)/XiS(tIG) …(3) 第7図は、上記点火時予測偏差を定性的に評価
するためのメンバーシツプ関数である。図中の
EPMB(i)(i=1〜7)はメンバーシツプ関数の形
を規定する定数であり、PB,PM,PS,ZO,
NS,NM,NBは偏差EPの大きさを定性的に評
価するためにメンバーシツプ関数に与えた名称で
あり、意味するとろころは初期値予測機能200
におけるものと同じである。また、図の縦軸はメ
ンバーシツプ値である。 第8図は、予測したボイラ点火時の初期値に応
じて適切な基本スケジユールを作成するための知
識を整理したものであり、点火時予測偏差に応じ
て、どのスケジユールパラメータを修正するのが
効果的であるかを整理したものである。これに従
つて作成した基本スケジユール作成ルールの一例
を第9図に示す。本図は、ボイラ点火時ドラム温
度を偏差EP(1)からボイラ起動所要時間修正量DP
(1)、第1速度保持時間修正量DP(2),……第3負
荷保持時間修正量DP(7)を定性的に決定するため
のものである。例えば、ルールNo.6の場合は、 IF(EP(1) is PM) THIN(DP(1) is PM and DP(2) is PS and DP(3) is PS and DP(4) is PS and DP(5) is PS and DP(6) is PS and DP(7) is PS) という意味である。 第10図は、定性的に決定されたスケジユール
パラメータ修正量を定量的な値に変換するための
メンバーシツプ関数である。図中のDPMB(i)(i
=1〜7)はメンバーシツプ関数の形を規定する
定数である。PB,PM,PS,ZO,NS,NM,
NBは修正量の大きさを定性的に表わすためにメ
ンバーシツプ関数に与えた名称であり、第7図で
使用している名称に対応している。また、図の縦
軸はメンバーシツプ値である。適用されたルール
により、スケジユールパラメータの修正量が定性
的にどのメンバーシツプ関数に属するかが決定さ
れる。或るスケジユールパラメータに対する修正
値が複数のルールによる複数のメンバーシツプで
規定された場合は、、各メンバーシツプ値に応じ
た加重平均値をもつて実際の定量的修正量520
(第1図参照)とする。これにより基本スケジユ
ールが作成されたことになる。 (3) スケジユール最適化機能400 上記基本スケジユール作成機能300により作
成された基本スケジユールを基にして、最適スケ
ジユール即ちボイラ点火から起動完了までの全起
動過程において運転制限条件を満足し、かつ起動
所要時間を最小とする前記スケジユールパラメー
タを決定するのが本スケジユール最適化機能40
0である。第11図に本機能の全体処理手順を示
す。 本機能では、運転制限条件であるタービン熱応
力の起動特性を評価するためにプラントの動特性
モデル100を用いる。起動時の熱応力パターン
とスケジユールパラメータは大きな相関があり、
熱応力の制限値に対する余裕(以下、マージンと
呼ぶ)が小さい程短時間の起動が可能となる。し
かし、コンプレツクス法などを利用した従来の制
御条件付非線形最適化アルゴリズムを用いると、
最適解(最適スケジユール)を得るのに、本例の
ように7変数がパラメータとなる場合、少なくと
も100回程度の繰返し計算(動特性モデルによる
起動特性の計算)が必要となり、現実的でない。
そこで、フアジイ推論を応用した最適化アルゴリ
ズムとすることにより収束性の大幅な改善を図
る。 プラント運転員は、起動前に熱応力特性の予測
値を与えられた場合、そのマージンに応じてどの
パラメータをどの程度短縮できるか経験的に知つ
ている。この経験的かつ定性的な知識を活用して
最適化のためのスケジユール修正量を決定する。
具体的には、まず第11図に示すように、基本ス
ケジユール作成機能300から与えられたスケジ
ユールパラメータpi(i=1〜7)で定まる基本
スケジユールに従つてプラントを起動した場合の
タービン熱応力特性を動特性モデル100を用い
て予測する。ここで、動特性モデル100は第1
2図に示すように、スケジユールパラメータが与
えられた場合、ボイラ点火指令、タービン速度及
び負荷の目標値を計算するためのスケジユール計
算機能110と、ボイラの起動特性を計算するた
めのボイラモデル120と、ボイラから発生する
蒸気条件を受けてタービンの熱応力を計算するた
めのタービンモデル130から構成されている。
ここで計算するタービン熱応力は、高圧タービン
のロータ表面応力、同ボア応力、中圧タービンの
ロータ表面応力、同ボア応力の4個所であり、い
ずれもタービン起動時に着目すべき重要な運転制
限要因である。第11図に示す起動特性評価機能
140は、起動過程を7つの区間に分割し、各区
間毎の最小応力マージンmj(j=1〜7)を求め
る。本図ではmjを区間jにおける高圧ロータ表
面応力と中圧ロータ表面応力のうち最小マージン
MS(j)と高圧ロータボア応力と中圧ロータボア応
力のうち最小マージンMB(j)の両者の意味で示し
ている。本起動特性評価機能140は、次に続く
応力マージン評価機能150で熱応力パターンの
特徴を抽出するためのものである。 応力マージン評価機能150では、第13図に
示すメンバーシツプ関数を用いて、ロータ表面応
力マージンMS(j)及びロータボア応力マージンMB
(j)を定性的に評価することにより、熱応力パター
ンの特徴を抽出する。図中のMMB(i)(i=1〜
6)はメンバーシツプ関数の形を規定する定数で
あり、PB,PM,PS,ZO,NS,NM,NBは応
力マージンMS(j)及びMB(j)の大きさを定性的に評
価するためにメンバーシツプ関数に与えた名称で
ある。また、図の縦軸はメンバーシツプ値であ
る。 スケジユール最適化ルール670は、「熱応力
パターンがどの様であればどのスケジユールパラ
メータをどの程度修正すれば良い」といつた断片
的知識である。第14図は、動特性モデル100
を用いて予測したタービン熱応力パターンに応じ
てスケジユールを修正するための知識を整理した
ものである。本図は、起動過程の各区間における
最小マージンに応じて、どのスケジユールパラメ
ータを修正するのが効果的であるかを整理したも
のである。ここで、MS(1),MS(2),……MS(7)及び
MB(1),MB(2),……MB(7)は、それぞれ第1、第
2、……第7区間におけるロータ表面最小応力マ
ージン及びロータボア最小応力マージンである。 第15図は、上記考え方に従つて作成したスケ
ジユール最適化ルールの一例を示すものである。
本図は、ロータ表面応力に関して、第5区間最小
応力マージンMS(5)と第6区間最小応力マージン
MS(6)からスケジユールパラメータの修正量とし
ての第3速度保持時間修正量を量DPT(4)、第1負
荷保持時間修正量DPT(5)及び第2負荷保持時間修
正量DPT(6)を定性的に決定するためのものであ
る。例えば、ルールNo.54の場合は、 IF(MS(5) is PB and MS(6) is PM) THEN(DPT(4) is NM and DPT(5) is NM and DPT(6) is NS という意味である。 第16図は、定性的に決定されたスケジユール
パラメータ修正量を定量に変換するためのメンバ
ーシツプ関数である。図中のDPTMB(i)(i=1〜
7)はメンバーシツプ関数の形を規定する定数で
ある。PB,PM,PS,ZO,NS,NM,NBは修
正量の大きさを定性的に表わすためにメンバーシ
ツプ関数に与えた名称である。また、図の縦軸は
メンバーシツプ値である。 第11図に示す修正パラメータ選択機能160
は、応力マージン評価機能150により抽出され
た応力パターンの特徴とスケジユール最適化ルー
ル670とを照合させることにより修正パラメー
タを選択し、適用されたルールによりパラメータ
の修正量が定性的に第16図のどのメンバーシツ
プ関数に属するかが決定される。 或るスケジユールパラメータに対する修正量が
複数のルールによる複数のメンバーシツプで規定
された場合は、各メンバーシツプ値に応じた加重
平均値をもつて実際の定量的修正量530(第1
図参照)とする。これを行なうのが第11図に示
すスケジユール修正量決定機能170である。 以上でスケジユールパラメータの修正量が決定
され、修正されたスケジユールに従つて、プラン
トを起動した場合の起動特性を予測するために、
再び動特性モデル100を動作させる。以上を繰
返すことにより最適起動スケジユール420(第
1図参照)を求めることができる。なお、第11
図中の収束判定機能180は、上記繰返し計算に
おいて作成される起動スケジユールの最適性を判
定するためのものである。また、その判定基準
は、全起動過程においてタービン熱応力が制限値
以下となる起動スケジユールのうち起動所要時間
が前回までのものと比較した時間短縮率が所定値
以下となつた場合とし、その中で最も起動所要時
間が小さくなる起動スケジユールを最適スケジユ
ール420(第1図参照)とする。 ここで決定された最適スケジユール420は、
第1図に示したように、ユーザインターフエイス
4000を介して表示装置5000に表示される
と共に、起動制御機能3000に設定される。起
動制御機能3000に設定された最適スケジユー
ル420は運転員6000からの実行指令430
を受けて実行される。起動制御機能3000は、
これを実行するためにプラント2000からのプ
ロセス入力3010を受けて、制御出力3010
をプラントに与える。 運転員6000は、知識ベース600の追加、
変更、削除が必要となつた場合は、ユーザインタ
ーフエイスを介して知識ベース管理情報690に
より、これらを実行する。 第17図は、フアジイ推論を応用したスケジユ
ール最適化機能400におけるスケジユール最適
化過程を示すものである。図中の番号は繰返し計
算の第1、3、5、20回目に作成された起動スケ
ジユールとそのときの起動特性を示す。ここで、
第1回目は基本スケジユールに対応し、20回目は
最適スケジユールに対応する。本図からわかるよ
うに、5回目でほぼ最適値に近いスケジユールが
得られており、本アルゴリズムによる収束性は極
めて良好である。また、作成された起動スケジユ
ールは併入指定時刻(本図では解列時刻からの経
過時間で示してあり480分(8時間)通りに起動
できることを示している。タービンの熱応力マー
ジンが大きな起動前半では速度保持や負荷保持を
省略することにより起動時間の短縮する方向でス
ケジユールが修正されており、熱応力マージンが
負となる起動後半では負荷保持を延長することに
より熱応力の緩和を図つている。このように、本
スケジユール最適化機能400によると、運転制
限条件を満足し、最短起動時間で指定時刻通りの
起動が可能となる。 スケジユール誤差監視機能700について次に
説明する。スケジユール誤差監視機能700は次
の2つの機能でスケジユール誤差の発生を監視す
る。 1つは、実際にスケジユール誤差が発生したこ
とを監視する誤差検知機能であり、別の1は、実
際に誤差が発生する前に誤差の発生を予測する誤
差予測機能である。 誤差検知機能は、起動前に作成された起動スケ
ジユールパラメータで規定される起動パターンが
実際の起動中に正確に遂行されているか否かを監
視するものである。したがつて、この起動パター
ンのずれを監視するものである。 一方、誤差予測機能は、起動パターンのずれと
なつて現われる前に、ずれの原因となるプラント
状態量の挙動を監視しておき、起動スケジユール
作成時に動特性モデル100で予測したプラント
状態量の挙動と比較し、その差が所定値以上とな
つたとき、スケジユール誤差の発生を予言するも
のである。 上記2つの機能でスケジユール誤差を監視し、
誤差発生を検知あるいは予測できた場合は、その
旨710を関連するプラント状態量と共にユーザ
インターフエイス4000を介し、表示装置50
00に表示する。これにより、運転員6000は
速かに対応処置をすることが可能となる。 (発明の変形例) 以上述べた本発明の実施例では、スケジユール
パラメータとしてボイラ起動所要時間とタービン
速度保持時間及び負荷保持時間に着目している
が、必ずしもこれらに限る必要はなく、タービン
速度変化率や負荷変化率あるいはボイラ昇温速度
や昇圧速度などプラントの起動パターンを代表す
るパラメータであれば本発明は基本原理を変える
ことなく実施できることは明らかである。 また、本発明の実施例では、運転制限条件とし
てタービン熱応力に着目しているが、必ずしもこ
れだけに限る必要はなく、タービン入口蒸気温度
及びその変化速度、あるいはタービンケーシング
温度などタービン熱応力を間接的に推定できる状
態量、あるいは、タービン伸び差、ボイラ蒸気圧
力上昇速度、ボイラ燃焼ガス温度など運転上重要
となる制限要因はプラントにより異なるから、必
要に応じて考慮すれば良く、本発明は本質を変え
ることなく実施できることは明らかである。 また、本発明の実施例では中央給電指令所から
指定される併入時刻を正確にまもつた起動スケジ
ユールを作成する方法を例に説明したが、本発明
を実施するにあたり、必ずしも併入時刻が指定さ
れる必要はなく、プラントに応じてボイラ点火時
刻、タービン通気時刻、目標負荷到達時刻などで
あつても本発明の本質が変らないのは明らかであ
る。 また、本実施例でプラント初期値として着目し
ている状態は、本発明を実施するプラントの機器
構成及び計測位置に応じて適宜選定すべきもので
必ずしも本実施例と同じものを使う必要のないこ
とは勿論であり、これにより本発明の本質が変わ
るものではない。 また、本実施例では、プラント冷却特性を解列
時刻からの時間関数として表わしているが、ボイ
ラ消火時刻などプラント冷却過程に移行する時刻
を基準とした時間関数で表わしたとしても、本発
明の実施にあたり本質が変わるものではない。 また、本実施例のスケジユール最適化機能にお
ける収束判定方法を必ずしも採用する必要はな
く、繰返し計算を所定回数だけ実施し、その中で
運転制限条件を満足し、起動所要時間が最小とな
る起動スケジユールを最適値と見做す方法、ある
いは運転制限条件を満足する所定個の起動スケジ
ユール候補が得られるまで繰返し計算を実施し、
その中で起動所要時間が最小となるものを最適値
と見做す方法を採用しても、本発明の本質が変わ
らないのは勿論である。 さらに、本実施例で用いているメンバーシツプ
関数は全て三角型としたが、必ずしもこの形に拘
わる必要はなく、プラントの特性及び運転員の知
識に応じて、2次曲線や指数曲線を採用しても、
本発明の本質が変わるものではない。また、メン
バーシツプ関数の形だけでなく、その数も任意に
設定しても、本発明の本質が変わるものではな
い。 〔発明の効果〕 本発明によれば、プラントの起動所要時間を従
来方式と比較して約30%短縮することができるた
め、負荷需要の変動に伴い発電プラントの頻繁な
起動停止が必要となる電力系統の安定かつ経済的
な運用が可能となる。また、これにより運転員の
負担が大幅に軽減される。また、起動所要時間の
短縮に伴つて起動時のエネルギー損失も約15%低
減できるため発電プラントの運用コストを大幅に
低減できる。さらに、運転制限条件および指定時
刻を忠実に守つた起動ができるため、プラント運
用上の安全性を向上し、電力系統への正確な電力
供給が可能となる。
【図面の簡単な説明】
第1図は発電プラント起動装置の全体構成を示
す。第2図は起動スケジユール最適化の基本的考
え方を示す。第3図は現在値偏差評価用メンバー
シツプ関数を示す。第4図は現在値偏差がボイラ
点火時初期値へ及ぼす影響度を示す。第5図は初
期予測ルールの一例を示す。第6図は初期値修正
量変換用メンバーシツプ関数を示す。第7図は点
火時予測偏差評価用メンバーシツプ関数を示す。
第8図は点火時予測偏差と修正対象スケジユール
パラメータの関係を示す。第9図は基本スケジユ
ール作成のルールの一例を示す。第10図はスケ
ジユールパラメータ修正用メンバーシツプ関数を
示す。第11図はスケジユール処理機能における
全体処理手順を示す。第12図は動特性モデルを
示す。第13図は応力マージン評価用メンバーシ
ツプ関数を示す。第14図は応力マージンと修正
対象スケジユールパラメータの関係を示す。第1
5図はスケジユール最適化ルールの一例を示す。
第16図はスケジユールパラメータ修正用メンバ
ーシツプ関数を示す。第17図はスケジユール最
適化過程を示す。

Claims (1)

  1. 【特許請求の範囲】 1 燃料の燃焼により発生する熱により蒸気を発
    生させるためのボイラと、発生した蒸気の熱エネ
    ルギを機械エネルギに変換するための蒸気タービ
    ンと、交換された該機械エネルギを電気エネルギ
    に変換するための発電機とで構成される発電プラ
    ントにおいて、該発電プラントの起動開始前に起
    動スケジユールとしての機器の操作時期及び制御
    目標を決定するための起動前スケジユール作成手
    段と、起動を開始した後に実際に遂行される起動
    スケジユールと該起動前スケジユールとの差を監
    視するスケジユール誤差監視手段と、該スケジユ
    ール誤差が所定値以上となつた場合にその旨を運
    転員に知らせるための表示手段を有することを特
    徴とする発電プラント起動装置。 2 特許請求の範囲第1項記載の該スケジユール
    誤差監視手段において、該スケジユール誤差が所
    定値以上となつた場合、その旨を関連する状態量
    と共に運転員に知らせるための表示手段を有する
    ことを特徴とする発電プラント起動装置。 3 特許請求の範囲第2項記載の該スケジユール
    誤差監視手段において、関連する該状態量とし
    て、起動前スケジユール作成時に予測したプラン
    ト起動特性とそれに対応した実際の起動時のプラ
    ント状態量とすることを特徴とする発電プラント
    起動装置。
JP23902086A 1986-10-09 1986-10-09 発電プラント起動装置 Granted JPS6394010A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP23902086A JPS6394010A (ja) 1986-10-09 1986-10-09 発電プラント起動装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP23902086A JPS6394010A (ja) 1986-10-09 1986-10-09 発電プラント起動装置

Publications (2)

Publication Number Publication Date
JPS6394010A JPS6394010A (ja) 1988-04-25
JPH0573886B2 true JPH0573886B2 (ja) 1993-10-15

Family

ID=17038695

Family Applications (1)

Application Number Title Priority Date Filing Date
JP23902086A Granted JPS6394010A (ja) 1986-10-09 1986-10-09 発電プラント起動装置

Country Status (1)

Country Link
JP (1) JPS6394010A (ja)

Also Published As

Publication number Publication date
JPS6394010A (ja) 1988-04-25

Similar Documents

Publication Publication Date Title
US4868754A (en) Method of starting thermal power plant
RU2213997C2 (ru) Способ оценки температуры лопаток в паровой турбине
CN113991647B (zh) 一种面向频率响应容量规划的电力系统随机生产模拟方法
JP2002155708A (ja) 発電プラントのガイダンス提供システムおよび方法
JP3666035B2 (ja) 火力発電プラント自律適応制御システム
US11287126B2 (en) System and method for controlling operation of boiler
JP3333674B2 (ja) プラント運転計画作成方法および装置
JPH08339204A (ja) 火力発電プラント自律適応最適化制御システム
Prasad et al. A novel performance monitoring strategy for economical thermal power plant operation
Arastou et al. Modeling and parameter estimation of a steam power plant including condenser back-pressure uncertainty using operational data
JP3666036B2 (ja) 火力発電プラント起動制御システム及び起動制御方法
Holmes et al. Fuzzy damage-mitigating control of a fossil power plant
JPH0573885B2 (ja)
JPH0573886B2 (ja)
JPS6396401A (ja) 発電プラント起動装置
Matsumoto et al. Turbine control system based on prediction of rotor thermal stress
Matsumoto et al. An operation support expert system based on on-line dynamics simulation and fuzzy reasoning for startup schedule optimization in fossil power plants
Aslanidou et al. Towards an integrated approach for micro gas turbine fleet monitoring, control and diagnostics
JPS6394006A (ja) 発電プラント起動装置
JPS6394007A (ja) 発電プラント起動装置
JPS6394005A (ja) 発電プラント起動装置
JP2678216B2 (ja) 火力発電プラント起動時運転支援装置
JPS6394004A (ja) 発電プラント起動装置
JPS6394009A (ja) 発電プラント起動装置
JPS6394003A (ja) 発電プラント起動装置