JPS59231604A - 火力発電プラントの運転制御方法 - Google Patents

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 〔発明の利用分野〕 本発明は、火力発電プラントの運転制御方法に係り、特
に、プラントの肉厚部に発生する熱応力を許容値以下に
抑制しながら急速な起動等を可能とするに好適な火力発
電プラントの運転制御方法に関するものである。

〔発明の背景〕

最近の火力発電プラントは、原子力発電プラントとの調
和をとるため、中間負荷運転用として溝成されている。

斯かる火力発電プラントは、与えられるプラントの起動
及び運転の各パターンに基づいて得たプラント運転パラ
メータに従って蒸気発生装置及びタービンの運転を制御
する。斯かる火力発電プラントは、急速な起動停止、大
幅な負荷変化に追従することが要請されている。したが
って、蒸気発生装置、タービン等の圧力部材の肉厚部に
発生する熱応力を精度よく求め、その部分の部材の寿命
消費量を抑えて起動・停止等の運転制御をすることが重
要である。プラント起動時において、大きな熱応力を発
生する箇所には、蒸気発生装置の２次過熱器出力ヘッダ
管寄部と、タービン第１段ラビリンスパツキン部のロー
タ表面及びボア（中心孔）とがある。

これら部材の応力を求めることや、その応力計算のもと
となる部材の温度分布を実測することは甚だしく困難で
あり、後者については回転体であることから特に困難で
ある。加えて、蒸気状態が時々刻々と変化するため、発
生する熱応力を精度よく求めることは困難である。した
がって、起動前の蒸気条件に応じて起動スケジュールを
含む運転パラメータを決定していた。しかしながら、こ
の方法の場合は、起動過程において、蒸気条件が予定値
からずれることを予め見込んで余裕値を大きくとってい
るため、起動時間が必要以上に長くなりがちであった。

しかも、熱応力の発生のもととなる蒸気温度の制御は、
応答遅れが大きく、熱応力に基づき先行制御を行うこと
は困難であった。

上記のことを考慮して、最近、タービンルータに発生す
る熱応力を予測し、昇速率、負荷変化率等の運転パラメ
ータを修正する方法が提案されている。しかしながら、
斯かる方法は、蒸気発生装置側で発生する蒸気条件を前
提とし、かつ蒸気発生装置側の蒸気状態の制御と独立し
であるため、途中で昇ね停止するなどプラント相互協調
が図れず、起動時間は必要以上に長くなりがちであった
。

〔発明の目的〕

本発明の目的は、プラントの１回の起動及び運転に許容
される寿命消費量を有効に利用でき、かつプラントの肉
厚部に発生する熱応力を許容値以下に抑えながら起動時
間を最短にする火力発電プラントの運転制御方法を提供
することにある。

〔発明の概要〕

本発明は、上記目的を達成するため、主蒸気温度昇温率
、タービン昇速率、負荷変化率等の各プラント状純変化
率に門するプラント運転パラメータを仮に設定し、この
設定されたプラント運転パラメータを基に主蒸気温度の
状態量の変化率を予υ１１ｊシ、この予測値に基づいて
蒸気発生装置及びタービンの応力評価部の熱応力を予測
し、予め設定した１回の起動及び運転に許されるスナ命
消費付に対応して決足された熱応力許容値と比較し、当
該許容値との偏差の小さい方で上ｉピプラント状態度化
率が最大、つまシブラント起動を最短時間とする運転パ
ラメータを得て、これをもって蒸気発生装置及びタービ
ンの運転を制御することを特徴とする。

〔発明の実施例〕

以下、本発明の実施例を図面に基づいて説明する。

第１図は、本発明に係る火力発電プラントの運転制御方
法が適用される火力プラントの概念的な構成を示すブロ
ック図である。

図中符号１０は操作盤、２０はディジタル計算機、３０
は石炭ミルシステム、４０は蒸気発生装置（以下ボイラ
システムと称す）、５０はタービン発電機システムであ
る。

火力プラントのオペレータは、＠ｔｇ＋ａ２０を介して
与えられるプラントの谷部のデータ、図示しない中央給
電指令所のような上位の制御所から与えられるデータ等
を参考に必要な操作を操作盤１０上にて行う。計ｑ６機
２０は事前に与えられているプログラムに従って、プラ
ントの谷部のデータ、操作盤１０の操作に応じた信号等
から、プラントの各部に対して必要な制御信号を与える
。石炭ミルシステム３０は、コールバンカ３０１、給炭
機３０２、ミル３１０、プロワ３２１及び３２２、ダン
パ３２３及び３２４を含んで構成されている。

バンカ３０１、給炭機−３０２を介してミル３１０内に
供給された石炭は、ミル３１０内で粉砕される。このよ
うに微粉炭にされた石炭は、プロワ３２１及び３２２に
より吹込まれた空気により、ボイラシステム４０のバー
ナ４０７に搬送され、ここで燃焼する。

ボイラシステム４０は、給水ポンプ４０１、給水制御弁
４０２、エバポレータ４０３、−次ス−パヒータ４０４
、二次スーパヒータ４０５、ガス書循項プロワ４０６及
びバーナ４０７より植成され°Ｃいる。該給水ポンプ４
０１より供給δれた水は、エバポレータ４０３で蒸気に
され、−次、二次スーパヒータ４０４及び４０５で過熱
蒸気である主蒸気とされ°Ｃ後述するタービン発電機シ
ステム５０に供宿される。バーナ４０７の供伶する熱針
は、エバポレータで水を蒸気に度え、スーパヒータで蒸
気を過熱然気にするために使用さｈるとともに、一部は
排出ガスとともに煙突から梗気中に逃げる。煙突から排
出されるガスの一部はガス書循環ブロソ４０６によって
ボイラに戻される。

ボイラが供給する蒸気流量を制御する１ヒめ、計算機２
０の出力により制御弁４０２が操作される。

主蒸気温度を制御するためのボイラ関連データとして給
水量、二次スーパヒータ入口温度、蒸気流量、主蒸気温
度、主蒸気圧力及びガス再循環流量等がそれぞれ検出器
４１１，４１２，４１３゜４１４　、４１５及ヒ４１６
−Ｃ’検出すｉ、Ｍｔｎ観２０へ送られる。

タービン発電機システム５０は、タービン制御弁５０１
、高圧タービン５０２、中低圧タービン５０３、ｔｉ水
器５０４及びタービンロータに直結された発電機５０５
より構成される。計算機２０の出力により制御弁５０１
が操作され、その開度に応じた量の主蒸気がタービン５
０２，５０３に供給され、発電機５０５を回転させる。

タービンで仕事を終えた蒸気は狽水器５０４によシ水に
戻される。この水は前述した給水ポンプ４０１によりボ
イラに再度供給される。発電機出力が検出器５０９で検
出され計算機２０へ送られる。

計ｎ機２０は、操作盤１０からプラント運転に　　　′
関する独々の要求を与えられ、プラントから得られるデ
ータと事前に与えられているプログラムから、所定の制
御信号を出力し、プラント目標状態に制御する。

以上火力プラントの構成の概要を説明した。

次に上記火力発電プラントの起動時における熱応力の発
生機構について説明する。

火力発電プラントにおいて起動時における熱応力の発生
する部分は、既に述べたように、大別すると、タービン
の１段ラビリンスパツキン部と、ボイラシステム４０の
２次過熱器出ロヘッダ管寄部とである。そこで、まず該
管寄部の熱応力の算出過程を第２図を用いて説明し、次
いでタービンのラビリンスパツキン部の熱応力の算出過
程を鋲３図を用いて説明する。

第２図は、ボイラの２次過熱器４０５の出口ヘッダ管寄
部４０５Ａを示す断面図である。この部分は、外面は加
熱されず、内面のみが蒸気（内部流体）によって加熱さ
れる。しかして、半径方向の肉厚が大きいことから、内
外面の温度差が顕著に現われる。このため、特にノズル
コーナ部ＮＣには、大きな熱応力が生じることとなる。

この部分ＮＣに生ずる熱応力を推定するため、主蒸気流
量ＭＳＦ、主蒸気温度ＭＳＴ及び主蒸気圧力ＭＳＰをそ
れぞれ検出器４１３，４１４及び４１５で検出して用い
ることにする。

それでは、まず部材の温度分布について検討する。

２次過熱器出ロヘッダ管寄部４０５Ａを無限円筒とみな
すとすれば、主蒸気と管寄部部材との伝熱によるメタル
温度分布は、次式で与えられる。

ここに　Ｔ：時刻ｔ１円筒中心から半径ｒの点のメタル
温度 α：メタル温度拡散率 境界値条件として、主蒸気とメタル内面、メタル外面と
外部との間に次式が成立する。

ただし、ａ：円筒の内半径 ｂ＝円筒の外半径 Ｔｏ：時刻ｔの円筒内面（ｒ＝ａ）のメタル温度 Ｔｔ：時刻ｔの主蒸気温度 λ：メタル熱伝導率 ｈ：主蒸気とメタル間の伝達率 また、ｈは次式で与えられる。

ｈ＝０．０２３１＋１ＩＯ−８Ｐｒ”’・−・・−・・
（３）ａ ただし、Ｋ：流体（主蒸気）熱伝達率 几、：レイノズル数 Ｐ、ニブラントル数 次に、上記温度分布に基づく半径ｒの点の熱応力を極座
標表示で求めると、半径方向熱応力σ、（す、周方向熱
応力σ０（す、軸方向熱応力σ、（ｒ）は次式で与えら
れる。

Ｅα　　２ｂ ａ、＜ｒ）＝−１−／　　Ｔ（ｒ）ｒａｒ−Ｔ（ｒ））
　　　　　−（６１１−νｂ２ａ２ｍ ここに、Ｅ：ヤング率 α；線膨張率 シ：ボアソン比 しかして、熱応力が最も厳しいのは、既に述べたが、内
面のノズルコーナ部Ｎ　Ｃである。−収部の熱応力σθ
（ａ）＝σ１（ａ）に応力集中定数係数Ｃを乗じ、この
部分の熱応力をσとすると、σ＝Ｃ・σθ（ａ） −−（−Ｊ　Ｔ（すｒｄｒ−Ｔ（ａ））　　　　　”’
　　（７）ＣＥα　　２ｂ ｌ−ν　ｂ２−３２　　。

以上のように、該コーナ部の熱応力を求めることができ
る。

ところで、ボイラ起動時の熱応力予測値を求めるために
は、上記（７）式からも理解できるように、主蒸気温度
の予測を行う必要がある。そこで、主蒸気温度の予測の
手法について一例を挙げて睨明するＯ 主蒸気温度の予測値は、２次過熱器の動特性モデルを用
いて、ｎ・τ（ｎは整数、τは演昇周期ン先の主蒸気温
度ＭＳ　Ｔ　（ｔ　＋　ｎτ）を次式のＸ（τ）をｎ回
繰り返すことによって求める。

すなわち、時刻ｔにおける主蒸気温度をＸｌ、２次過熱
器部拐温度をχ２．２次過熱器入口蒸気ｖＡＩＫをｕｌ
　、２次過熱器蒸気鴎度をｕ２とすると、である。また
、 − 十 ここに、 ＣＰ：主蒸気定圧比熱 Ｆ８　：２次過熱器内部流体（主蒸気）流量ＦＩＩＲ：
　２次過熱器内部流体（主蒸気）定格流量 ｒｓ：２次過熱器内部流体（主蒸気）比重危ｖ　：２次
過熱器内部流木（主蒸気）容積Ｆｚ　ＩＦ：ゲイ２再循
環ガス流量 Ｆｇ＊ｙｎ：ボイラ再循環ガス定格流最Ｍｍ　：２次過
熱器メタル重量 Ｃ，：２次過熱器メタル比熱 Ａ　：２次過熱器伝熱面積 α１１：蒸気からメタルへのだ格状態での熱伝達率 α１■：メタルから蒸気への定格状態での熱伝達率 ところでｕｅ式は、 Ｘ（ｉ）＝Φ（ｉ−ｔ）ｘ（ｉ−ｔ）＋Ｈ（ｉ−ｔ）　
ｕ（ｉ−ｘ）　　・・−（１７１と表現できる。この時
、プロセスの観測過ａが、ｙ（ｉ）＝　ｃ（ｉ）ｘ（ｉ
）＋ｗ（ｉ）　　　　　　　　　　・−ＱＳｙ（す：０
１次元観測ベクトル Ｃ（ｉ）：ｍＸｎ観測行列 ｗ（ｉ）：ｍ次元観測ノイズベクトル であるとすると、信号Ｘ（ｉ）の最大推定値ｘ（ｉ）は
、カルマンフィルタ理論を適用すると、次式で求められ
る。すなわち、 Ｘ　（ｉ）＝　ｘ　（ｉ）十ｐ（ｉ）Ｃ’（ｉ）Ｗ−’
　（Ｖ（ｉ）−（Ｃ（ｉ）Ｘ（ｉ）＋　Ｖｌ／（ｉ））
）　−ｔｉｌである。

ここにＸはモデルの推定量であり、 Ｘ（ｉ）＝Φ（ｉ−１１Ｘ（ｉ−１）＋Ｈ（ｉ−１）Ｕ
（ｉ　−１）　　−（２１で与えられる。

よって、カルマンフィルタを適用することによ９００式
で与えられるＸ（ｉ）は、より鞘朋の高い推定値が求ま
ることになる。

したがって、０１．６１式で求めたｎτ先の主蒸気温度
予測値を（１）式の温就分布計算に用い、次いで計算さ
れた温度分布値を（方式に適用することにより１１τ時
間先の熱応力が求まることになる。もちろん、主蒸気温
ｖ−の予測値は、例えばプラント運転パラメータとして
設定されているプラント状態量変化率から一定時間後の
ものを推定してもよい。

第３図は、タービン発電機システム５０の高圧タービン
５０２の第１没後ラビリンスパツキン部５４１を示す断
面図である。この部分は、タービン側でタービン起動時
に最も熱応力の発生が大きい部分であることは既に述べ
た。ここを洩れる蒸気の温度、圧力及び洩れ流速が、タ
ービン起動過程においては大幅に変動するため、この附
近のロータは過酷な条件にさらされることになる。この
ｌヒめ、当該部分は、加熱又は冷却の繰り返しが急激と
なり、過大な熱応力が発生しやすくなるのである。この
熱応力を推定するため、主蒸気温度、主蒸気圧力、第１
没後蒸気温度Ｔｌ８Ｔ及び第１没後蒸気圧力をそれぞれ
検出器４１４，４１５゜５０５及び５０６で検出して用
いる。

まず、ロータの部材の温度分布を求める。

２次過熱器出ロヘッダ管寄部４０５Ａと同様に無限円柱
とみなすと、ロータの温度分布は前記（１）式で与えら
れる。

ただし、α：ロータＩの温度伝導率 Ｔ：時刻ｔ、０−タ軸がら半径ｒの点 のロータ内温度 次に、ロータの温度を求める。

仮にロータを半径方向に軸に平行に６分割し、かつ、蒸
気条件を前記囲及び１！ｊ式を用いて、τ時間先のロー
タの表面温１Ｔｆ（ｔ＋τ）、ロータボア温度Ｔｂ（ｔ
＋τ）を求めると、 で与えられる。

λＪ　Ｂｔ：第１没後蒸気の熱伝導率 Ｎ６：ヌーセルト数 さらに、上記Ｎ、は、 で与えられる。

ここに、δ：パッキン間隙 ’１’ｌ［ｌＴ：第１没後蒸気温度 である。

上ｄＣの温度分布に基づく熱応力予測値は、次のように
与えられる。すなわち、ロータ表面熱応力Ｏｆ及びロー
タボアの熱応力σｂは、 である。

ただし、’Ｉ　ＭＳ　Ｓロータ表面での体積平均温度’
１”ＭＢ：ロータボアでの体積平均温度以上説明してき
たように熱応力を詩、出することができる。

上記説明からも理解できるように、蒸気子側の正確さが
熱応力を算出するには重要な要素であるといえよう。

次に、以上で求まる熱応力予測計算を用いたプラント起
動方法について第４図を用いて説明する。

第４図は、前記火力発電プラントのプラント起動曲Ｍを
示す波形図である。Ｍ４図において、横軸は時間ｔを示
し、縦軸は各挿置の大きさを示す。

第４図において、符号ＭＳＴは主蒸気ｆＵ［（Ｃ１、Ｔ
Ｖはタービン速度（ＲＰＭ）、ＰＬは負荷（ＭＷ＞、ま
た時刻１．は点火、時刻ｔ２は通気、時刻ｔ３は併入、
時刻ｔ４は弁切換、の荷時間を示している。

１）　ボイラ点火（１＋）からタービン通気（ｔ２）直
前。

この期間は、タービン５０２には通気をしていないので
、タービン５０２に関係なく、ボイラ熱応力ｔＦｒ’８
値一杯にしてボイラの昇温・昇圧；Ｕｌ」岬を行うこと
ができる。

ＩＩ）　　タービン通気（ｔ２　）から併入（ｔ３　）
まで。

この期向は、タービン本体の諸々の制約、例えば危険速
度の振動の問題等、制約が多く、タービン５０２を中心
として起！ｖＩ運転することが望ましい。これは、ター
ビン熱応力を予測計算し、その許容１直内で、最大の昇
速率を選択し、短時間で昇速することを意味している。

したがって、仁の期間におけるボイラ制御は、タービン
熱応力許容値内で最大昇温・昇圧制御を行わざるをえな
い。

１１１）併入（ｉｇ）から加減弁切換完了（ｔ４　）ま
で。

この期間においては、ボイラ温度が大きく変動しやすい
ものの、他方タービン５０２としては負荷変化が比較的
小さいという傾向を呈することから、ボイラシステム４
０は熱応力許容値で最大の昇温・昇圧率を選択しこれら
選択値に基づいて昇温・昇圧制御を行わせる。このよう
な条件下において、タービン熱応力ＩｒＦ８値内に入る
初負荷Ｉｒ（、弁切換から負荷までの負荷上昇、弁切換
負荷と負荷上昇パターン制ｔｉｌｌを行うこととする。

１ｖ）加減弁切換児了（ｔ４　）から全負荷（ｔ５　）
まで。

この期間は、ボイラとしては定格然気条件に極力早く、
また、タービンとしては最短で全負荷をとるように制御
する必要があるのはいうまでもない。

したがって、ボイラ熱応力とタービンの熱応力の計算値
に対する余裕度の小さい方で主Ｔｏ１１ν１を行い、他
はこれに追従制御をさせる。

上述のことを具体的に説明すると、例えば、タービン熱
応力許容値が小さい場合、タービン熱応力許容値内での
最大の負荷五化率、昇温率、外圧率を選択し、これをも
って制御を行い、ボイラシステム、４０は他のプラント
状態変化率で制御せしめる。

なお、運転条件によっては、負荷又は蒸気条件を最短で
定格に制御したい場合がある。

＋ＹＪ者に対しては、タービン熱応力の許容値内での最
大負荷変化率を選択し、昇温率、昇圧率はその負荷に追
従させる。

一方、後者に対しては、ボイラ熱応力許容値内での最大
昇温外圧率を選択し、その蒸気条件下で、最大負荷変化
率を選択して制御することとする。

以上１）〜＋Ｖ　＞に説明してきたように火力発電プラ
ントの運転状態に対応して熱応力に余裕をもち、安全か
つ、最短時間での起動ができるように運転する。

上記説明をまとめてみると、要するに、本発明に係る運
転制御方法は、ボイラタービンの起動運転時の状況に応
じて主に応力許容値との偏差が小さい方で最大プラント
状態変化率を選択し、これをもってボイラ又はタービン
を制御し他はこれに追従制御される。

上述の運転制御方法を第５図に示すブロックを用いて説
明する。

プラント起動に当って、ステップ２００では、運転員１
が操作盤１０を介して計７ｖ機２０の記憶装置の運転パ
ラメータ設定エリアに、プラント起動パターン及び運転
パターン、ボイラ（２次過熱器ヘッダ管寄部）熱応力ｉ
ｉ’ｆ　８値、高圧タービン（ｍ１段ラビリンスパツキ
ン部）ローフ熱応力許容値等の運転パラメータを設定す
る。ステップ２０１では、該記憶装置の所定のエリア部
に設定されたプラント起動及び運転パターンに基づき、
タービン５０２の負荷変化率及び昇速率と、ボイラシス
テム４０の昇温率及び昇圧率等の運転パラメータの最大
値をｎＱ憶装置の他のエリアに仮に設定してステップ２
０２に移る。次に、ステップ２０２では、どの程度光の
熱応力を予測して制御すべきかを、すなわち予測時間を
決定してステップ２０３に移る。ここで、ステップ２０
２における予測時間は、応力評価部（例えば、該管寄部
４０５Ａ、該ラビリンスパツキン部５４１）’Ｊの熱伝
達量の大きさで、すなわちプラントの運転状態をもって
決定するものとする。ステップ２０３では、前記決定し
た予測時間に基づき、蒸気条件を、例えば前述の（１６
）、ｔ１９式を用いて予測計算を行いステップ２０４に
移る。ステップ２０４では、この蒸気栄件予測値に基づ
きボイラシステム４０の応力計画部（２次過熱器ヘッダ
管寄部）の昌度分布をｉｎして２０５に移行する。ステ
ップ２０５では、その計ｐ、結呆から、２次過熱器熱応
力予測値を求める（これは上記仮設定した変化率での熱
応力予測値となる）。次に２０６で、運転員ｌによシ計
ｎ磯２０の記憶装置の設定エリアに設定された熱応力許
容値との余裕値を求め、ステップ２０７に移る。次にス
テップ２０４乃至２０６と同様に、ステップ２０７乃至
２０９において、タービンロータに関しても熱応力を求
め、最終的に熱余裕値を求め４ｔ７！：ａ２ｏの記憶装
置の所定の設定エリアに余裕値を格納する。しかして、
ステップ２０６及び２０９において求め、かつ計ｎ＋Ａ
２０の記憶装置の所定の設定エリアに格納した２次過熱
器４０５Ａ及びタービンロータ５４１の６各の熱応力許
容値の余裕値をステップ２１０にて取り込み、余裕値の
より小さい方を、すなわち低値を選択してステップ２１
１に移る。ステップ２１１では、第４図の説明の如く、
どの期間の運転かを判定するとともに、前述ｉＶ）項の
期間の場合は負荷上昇最短又は蒸気条件の最短時間での
定格値までの昇温、外圧等の優先条件を評価する、。

しかして、その結呆をステップ２１２に与える。

ステップ２１２では、運転員１により設定された熱応力
許容値と比較し熱応力許容値以下で、かクプラント状態
変化率が最大となるプラント運転パラメータを選択し、
計１１機２０の記憶装置の所定の設定エリアに仮設定さ
れたプラント状態変化率の値を修正させる。さらに、ス
テップ２１２では、ボイラ蒸気温度制御機能２２０に、
昇温率及び昇圧率を与える。また、ステップ２１２では
、タービン昇速負荷制御機能２３０に、昇速率及び負菊
変化率を与える。このような、動作をステップ２１２で
行いステップ２１３に移る。ステップ２１３では、上記
プラント状態変化率設定値を１足する毎に目標値（起動
又は運転の完了）に到達したかどうかを判定し、目標値
に達しない場合は、ステップ２０１に移る。またステッ
プ２１３で目標値に達した場合は、当該処理を終了とす
る。

〔発明の効果〕

以上述べたように、本発明によれば、熱応力を許容値以
下に抑えながら、ボイラ・タービンの協調制御をできる
ので、安全かつ急速なプラント起動及び運転が可能とな
るという効果がある。

【図面の簡単な説明】

第１図は本発明に係る火力発電プラントの運転制御方法
が適用される火力発・亀プラントの構成を示すブロック
図、第２図はボイラシステムの応力評価部の熱応力発生
及びその算定を説明するために示す図、第３図はタービ
ンシステムの応力評価部の熱応力発生及びその算定を説
明するために示す図、第４図は火力発電プラントの起動
曲線を示す線図、第５図は本発明に係る火力発電プラン
トの運転制御方法の一例を示すフローチャートである。 １・・・運転員、１０・・・操作盤、２０・・・計算機
、３０・・・ミルシステム、４０・・・ボイラシステム
、５０・・・タービン発電機システム。 代理人　弁理士　鵜沼辰之 芽２　目 ２０ ＄３目　　、。 ＄　４　目

Claims (1)

【特許請求の範囲】
1. １、与えられたプラントの起動及び運転の各パターンに
   基づいて得たプラント運転パラメータに従って蒸気発生
   装置及びタービンの運転を制御する火力発電プラントの
   運転制御方法において、上記パターンに基づいてタービ
   ンの昇速率及び負荷変化率と、蒸気発生装置の主蒸気昇
   温率及び主蒸気昇圧率との各プラント状態変化率に関す
   るプラント運転パラメータを仮り設定し、該設定された
   プラント状態変化率を基に指定時間後の主蒸気の状態量
   の変化率を予測し、この予測値をもって蒸気発生装置及
   びタービンの応力評価部の各熱応力をそれぞれ予測し、
   該予測された画然応力を予め設定された１回の起動及び
   運転で許容される寿命消費址に応じて定めた各熱応力許
   容値と比較し、許容値との差の小さい方で上記の如く設
   定されたプラント状態変化率が最大となる運転パラメー
   タを得ると共に、目標制御状態となるまで上記行程を繰
   シ返し、かつ得られたプラント運転パラメータをもって
   蒸気発生装置及びタービンの運転を制御する火力発電プ
   ラントの運転制御方法。