JPS6394006A - Starting device for power generation plant - Google Patents

Starting device for power generation plant

Info

Publication number
JPS6394006A
JPS6394006A JP23901586A JP23901586A JPS6394006A JP S6394006 A JPS6394006 A JP S6394006A JP 23901586 A JP23901586 A JP 23901586A JP 23901586 A JP23901586 A JP 23901586A JP S6394006 A JPS6394006 A JP S6394006A
Authority
JP
Japan
Prior art keywords
schedule
startup
plant
time
starting
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.)
Pending
Application number
JP23901586A
Other languages
Japanese (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 JP23901586A priority Critical patent/JPS6394006A/en
Publication of JPS6394006A publication Critical patent/JPS6394006A/en
Pending legal-status Critical Current

Links

Landscapes

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

Abstract

PURPOSE:To make it possible shorten the required time for starting the captioned plant by optimizing the starting schedule through the repetitive calculation method under the estimation of the starting characteristic obtained in a dynamic characteristic model, thus permitting the operation restriction condition to be satisfied. CONSTITUTION:A starting schedule preparing means 1000 quantitatively calculates the starting characteristic in a dynamic characteristic model 100 according to the standard schedule in the starting mode selected in response to the stopping time of a power generation plant 2000. In order to increase the estimation precision, an initial value estimating means 200 estimates the initial value based on the difference between the standard cooling characteristic and the actual cooling state through the application of the fuzzy estimation 500. A schedule optimizing means 400 optimizes the starting schedule through the repetitive calculation method, evaluating the starting characteristic, in order to permit the starting in the shortest time and satisfy the operation restriction condition by using the initial schedule of a fundamental schedule preparing means 300. The starting control 3000 is thus performed. Therefore, the plant can be safely started in the shortest time.

Description

【発明の詳細な説明】 〔産業上の利用分野〕 本発明は発電プラントの起動装置に係り、特に高精度起
動特性予測に好適な初期値予測方法、及び起動スケジュ
ール最適化計算の収束性向上に好適な初期スケジュール
の作成方法に関する。
[Detailed Description of the Invention] [Field of Industrial Application] The present invention relates to a startup device for a power generation plant, and particularly to an initial value prediction method suitable for predicting startup characteristics with high accuracy, and to improving convergence of startup schedule optimization calculations. This invention relates to a method for creating a suitable initial schedule.

〔従来の技術〕[Conventional technology]

火力発電プラントの起動に関する従来の方法は、起動前
の停止時間や機鼎の温度状態に応じて、ボイラへの初期
投入燃料量、主蒸気の昇温及び昇圧の時間関数、タービ
ンの昇速及び負荷上昇の時間関数を起動スケジュールと
して決定し、この起動スケジュールをプラントの各系統
に設けられた制御系で実行するという方法が採られてい
た。
The conventional method for starting up a thermal power plant is based on the initial amount of fuel input to the boiler, the time function of main steam temperature and pressure rise, the turbine speed increase and A method has been adopted in which a time function of the load increase is determined as a startup schedule, and this startup schedule is executed by a control system installed in each system of the plant.

二の最も代表的な方法は、Elactrical Wo
rld。
The second most typical method is Electrical Wo
rld.

Vofl・、165.&6の論文“Thermal 5
tressInfluance Starting、 
Loading of BoilersTurbine
s”で述べられている。この方法は、プラントの限られ
た部分の初期状態によって一義的に起動スケジュールを
決定する方法である。即ち、ボイラ蒸気圧力、ボイラ出
口蒸気温度、蒸気タービンケーシング温度の初期値に応
じて、蒸気タービンの昇速率、初負荷量、速度保持並び
に負荷保持による蒸気タービンの暖機時間及び負荷変化
率を決定する方法である。この方法によると、運転制限
要因である蒸気タービンの熱応力を管理する上で重要な
ボイラ発生蒸気の昇温特性を起動前に予測できないため
、その不確定性を起動スケジュールに余裕を持たせるこ
とにより吸収している。
Vofl・, 165. &6 paper “Thermal 5
tressInfluence Starting,
Loading of BoilersTurbine
This method determines the startup schedule uniquely based on the initial conditions of a limited portion of the plant.In other words, the startup schedule is determined based on the initial conditions of a limited portion of the plant. This method determines the steam turbine speed increase rate, initial load amount, speed maintenance, and load maintenance time and load change rate of the steam turbine according to the initial values. Since the temperature rise characteristics of the boiler-generated steam, which is important for managing the thermal stress of the turbine, cannot be predicted before startup, this uncertainty is absorbed by allowing some leeway in the startup schedule.

そのため、作成される起動スケジュールは必要以上に長
くなりがちであった。
Therefore, the created startup schedule tends to be longer than necessary.

また、別の従来方法としては、U S P3,446.
224及びU S P4,228,359が知られてい
る。これらは蒸気タービンに発生する熱応力をオンライ
ンリアルタイムで監視しながら蒸気タービンの急速起動
を図るものであるが、前記従来方法と同様にボイラの起
動方法に関しては何ら言及していない。
Another conventional method is US P3,446.
No. 224 and US P 4,228,359 are known. These methods attempt to quickly start up a steam turbine while monitoring the thermal stress generated in the steam turbine online in real time, but, like the conventional methods described above, there is no mention of a method for starting the boiler.

ボイラの起動時間の短縮を目的とした従来方法としては
、特開昭59−157402が知られている。この方法
はボイラに発生する熱応力をオンラインリアルタイムで
監視しながらボイラ発生蒸気の急速昇温を図るものであ
る。しかし、この方式は蒸気タービンの起動に関しては
何ら言及していない。
As a conventional method aimed at shortening the start-up time of a boiler, Japanese Patent Application Laid-Open No. 157402/1984 is known. This method aims to rapidly raise the temperature of the steam generated by the boiler while monitoring the thermal stress generated in the boiler online in real time. However, this method makes no mention of starting the steam turbine.

プラント全体の起動時間は、ボイラと蒸気タービンの協
調により短縮が可能であるが、以上述べた従来の方法は
、何れもボイラもしくは蒸気タービンの片方のみに着目
した急速起動方法であり、この様な個別の方法を組合わ
せたとしてもプラント全体の起動時間が最短となる保証
は何も無い。
The startup time for the entire plant can be shortened by coordinating the boiler and steam turbine, but the conventional methods described above are all rapid startup methods that focus on only one side of the boiler or steam turbine. Even if individual methods are combined, there is no guarantee that the start-up time of the entire plant will be the shortest.

何故ならば、ボイラと蒸気タービンは相互干渉が極めて
強く、個々の最適化が必ずしも全体の最適化とならない
からである。
This is because the boiler and the steam turbine have extremely strong mutual interference, and individual optimization does not necessarily result in overall optimization.

さらに、上記従来方法では、ボイラ点火直前のプラント
初期温度状層の実測値に基づいて起動スケジュールを作
成するため、ボイラ点大前の任意一時刻にて起動スケジ
ュールを作成し、且つ中央給・電指令所(以下、中給と
呼び)から指定された起動完了時刻に正確に起動を完了
(以下、このことを定刻起動と呼ぶ)できる起動スケジ
ュールを作成することができなかった。
Furthermore, in the conventional method described above, the startup schedule is created based on the actual measured value of the initial temperature layer of the plant immediately before the boiler ignition, so the startup schedule is created at an arbitrary time before the boiler ignition, and the central supply/power supply It was not possible to create a startup schedule that would allow the startup to be completed accurately at the startup completion time specified by the command center (hereinafter referred to as "Nakagyo") (hereinafter referred to as "scheduled startup").

このように、従来方法では起動時間が長くなりがちであ
り、それに伴って起動時のエネルギ損失(以下、このこ
とを起動損失と呼ぶ)が大きくなるといったプラント運
用上の経済性の問題と、起動スケジュールを実行する上
での安全性の問題があった。
In this way, the conventional method tends to take a long time to start up, and as a result, energy loss during startup (hereinafter referred to as startup loss) increases, which is an economical problem in plant operation. There was a safety issue in executing the schedule.

〔発明が解決しようとする問題点〕[Problem that the invention seeks to solve]

本発明の目的は、前記従来方法における起動時プラント
初期状態の予測上の問題点を解決するために、運転員が
経験的に持っているプラントの冷却特性に関する知識を
活用し、初期状態の予測精度を向上させることと、繰返
し計算による起動スケジュール最適化の良好な収束性を
得るために、運転員が経験的に持っているプラント初期
状態と適切な起動スケジュールに関する定性的な知識を
活用して初期スケジュールの適正化を可能とする発電プ
ラント起動装置を提供するにある。
An object of the present invention is to solve the problem of predicting the initial state of the plant at startup in the conventional method by utilizing the knowledge of the cooling characteristics of the plant that the operator has experientially and predicting the initial state. In order to improve the accuracy and obtain good convergence of start-up schedule optimization through iterative calculation, operators' qualitative knowledge about the initial plant state and appropriate start-up schedule is used from experience. An object of the present invention is to provide a power generation plant startup device that enables optimization of an initial schedule.

〔問題点を解決するための手段〕[Means for solving problems]

上記目的は、プラントの起動特性を定量的に計算するた
めの動特性モデルと、この動特性モデルで使用するボイ
ラ点大時のプラント初期温度状態を予測するための初期
値予測手段と、最短時間での起動を実現するために上記
動特性モデルで得られる起動特性を評価しながら繰返し
計算手法により起動スケジュールを最適化するスケジュ
ール最適化手段と、上記繰返し計算用の初期スケジュー
ルを決定するための基本スケジュール作成手段を用いる
ことにより、達成できる。
The above objectives are to develop a dynamic characteristic model for quantitatively calculating the start-up characteristics of a plant, an initial value prediction means for predicting the initial temperature state of the plant when the boiler temperature is high, and a means for predicting the initial temperature in the shortest time. A schedule optimization means for optimizing the startup schedule using an iterative calculation method while evaluating the startup characteristics obtained by the above dynamic characteristic model in order to achieve startup, and a basic method for determining the initial schedule for the above-mentioned iterative calculation. This can be achieved by using a schedule creation means.

上記初期値予測手段及び基本スケジュール作成手段では
、ファジィ推論を適用して、熟練運転員の思考方法と似
た方法を用いることにより、それぞれ高精度起動特性予
測に必要な高精度初期値予測及び繰返し計算による良好
な収束性が得られる適切な初期スケジュールの作成が不
能となる。
The above-mentioned initial value prediction means and basic schedule creation means apply fuzzy inference and use a method similar to the thinking method of a skilled operator, thereby achieving high-precision initial value prediction and repetition necessary for high-precision starting characteristic prediction, respectively. It becomes impossible to create an appropriate initial schedule that provides good convergence through calculations.

〔作用〕[Effect]

ル辷従ってプラントを起動した場合の起動特性を定量的
に計算することができる。そのため、起動スケジュール
と起動特性の関係を、実際にプラントを起動する前に知
ることが可能であるとともに、プラント状態値が運転制
限条件を満足するか否かを事前に確認できるため、安全
性の高い起動スケジュールを作成できる。
Therefore, the startup characteristics when the plant is started can be quantitatively calculated. Therefore, it is possible to know the relationship between the startup schedule and startup characteristics before actually starting the plant, and it is also possible to confirm in advance whether the plant status values satisfy the operation limit conditions, which improves safety. You can create high startup schedules.

初期値予測手段は、熟練運転員の思考方法と似た方法で
初期値を予測するために、ファジィ推論を適用し、プラ
ントの標準的な冷却特性と実際の冷却状態との定性的な
差異により初期値を予測する方法を採っている。これに
よって、高精度の初期値予測が可能であり、その結果、
動特性モデルによる高精度な起動特性の予測が可能であ
る・スケジュール最適化手段は、熟練運転員の思考方法
と似た方法でスケジュールを最適化するために、ファジ
ィ推論を適用し、上記動特性モデルで得られた起動特性
を定性的に評価し、スケジュールを修正し、再度動特性
モデルを用いて起動特性を計算し、その結果を定性的に
評価し、スケジュールを修正するといった繰返し計算手
法により最適化を行なっている。このスケジュールと起
動特性の定性的因果関係に基づいてスケジュールの修正
量を決定するところが熟練運転員の思考方法と似ている
。これにより、最適スケジュール求解のための収束性が
極めて良好となるため、詳細な計算式を用いた大規模な
動特性モデルを用いることが可能で、起動特性を高精度
で予測することができる。
The initial value prediction method applies fuzzy reasoning to predict the initial value in a way similar to the thinking method of a skilled operator, and uses qualitative differences between the standard cooling characteristics of the plant and the actual cooling state. A method is used to predict the initial value. This enables highly accurate initial value prediction, and as a result,
It is possible to predict startup characteristics with high accuracy using a dynamic characteristics model.The schedule optimization means applies fuzzy inference to optimize the schedule using a method similar to the thinking method of a skilled operator. Using an iterative calculation method, we qualitatively evaluate the startup characteristics obtained by the model, modify the schedule, calculate the startup characteristics again using the dynamic characteristics model, qualitatively evaluate the results, and modify the schedule. Optimization is being performed. The way a skilled operator thinks is similar to the way a skilled operator thinks, in that the amount of schedule modification is determined based on the qualitative cause-and-effect relationship between the schedule and startup characteristics. As a result, the convergence for finding the optimal schedule is extremely good, so it is possible to use a large-scale dynamic characteristic model using detailed calculation formulas, and it is possible to predict startup characteristics with high accuracy.

基本スケジュール作成手段は、前記初期値予測手段で得
られた初期値に応じて、上記スケジュール最適化が用い
る初期スケジュールを作成するためのもので、予測した
初期値と標準的な初期値との定性的な差異により標準的
なスケジュールを修正することにより適切な初期スケジ
ュールが作成される。これにより、スケジュール作成に
関する熟練運転員の思考方法を真似ることができ、上記
スケジュール最適化の収束性を向上1きる。
The basic schedule creation means is for creating an initial schedule to be used in the schedule optimization according to the initial value obtained by the initial value prediction means, and is used to qualitatively compare the predicted initial value and the standard initial value. An appropriate initial schedule is created by modifying the standard schedule due to the differences. As a result, it is possible to imitate the thinking method of a skilled operator regarding schedule creation, and the convergence of the above-mentioned schedule optimization can be improved.

以上述べた手段により、プラントの運転制限条件を満足
し、起動所要時間を最小にし、中給から指定された時刻
通りに起動を完了できる最適起動スケジュールを作成す
ることが可能となる。
By means of the means described above, it is possible to create an optimal startup schedule that satisfies the plant operation restriction conditions, minimizes the time required for startup, and completes startup at the time designated by the intermediate supply.

(実施例〕 以下、本発明の一実施例を説明する。(Example〕 An embodiment of the present invention will be described below.

第1111は本発明を適用した発電プラント起動装置の
全体構成を示すものである0本装置は、ボイラ、蒸気タ
ービン及び発電機から成るプラント2000において、
ボイラ点大から目標負荷(−般的には中給から指令され
る負荷レベル)に到達するまでの全起動過程における機
器の操作時期及び制御目標を決定するための起動スケジ
ュール作成機能1oooと、作成された起動スケジュー
ルに従ってプラントを自動的に起動してゆくための起動
制御機能3000と、プラントの運転員6000からの
リクエストに応じてプラント起動に関する必要な情報を
表示装置115000に表示したり、スケジュール作成
機能1000内の情報を動特性モデル100.初期値予
測機能200.基本スケジュール作成機能300.スケ
ジュール最適化機能400.ファジィ推論機能500.
知識ベース600から成る。SOmべ−x6ooは、更
に、初期値予測ルール610.基本スケジュール作成ル
ール640.スケジュール最適化ルール670から成る
No. 1111 shows the overall configuration of a power plant starting device to which the present invention is applied. This device is used in a plant 2000 consisting of a boiler, a steam turbine, and a generator.
Startup schedule creation function 1ooo for determining equipment operation timing and control targets during the entire startup process from the boiler point high until reaching the target load (generally the load level commanded from the intermediate supply). A startup control function 3000 automatically starts the plant according to a set startup schedule, and displays necessary information regarding plant startup on a display device 115000 in response to a request from a plant operator 6000, and creates a schedule. The information in the function 1000 is converted into a dynamic characteristic model 100. Initial value prediction function 200. Basic schedule creation function 300. Schedule optimization function 400. Fuzzy inference function 500.
It consists of 600 knowledge bases. SOmbe-x6oo further includes an initial value prediction rule 610. Basic schedule creation rules 640. It consists of schedule optimization rules 670.

上記起動スケジュール作成機能1000を詳細に説明す
る前に、各機能の目的について説明しておく。
Before explaining the startup schedule creation function 1000 in detail, the purpose of each function will be explained.

動特性モデル100は、初期値予測機能200にて予測
されたボイラ点火時のプラントの状態を初期値210と
して、基本スケジュール作成機能300あるいはスケジ
ュール最適化機能400から与えられる起動スケジュー
ルに従ってプラントを起動した場合の起動特性140を
定量的に計算するためのものである。
The dynamic characteristic model 100 starts the plant according to the startup schedule given by the basic schedule creation function 300 or the schedule optimization function 400, with the plant state at the time of boiler ignition predicted by the initial value prediction function 200 as an initial value 210. This is for quantitatively calculating the starting characteristics 140 for the case.

初期値予測機能200は、起動前の任意時刻にト、艇い
て、ボイラ点火時のプラント状態を予測し、ミれを動特
性モデル100及び基本スケジュール′作成機能300
に設定するためのものである。
The initial value prediction function 200 predicts the plant state at the time of boiler ignition at an arbitrary time before startup, and generates a dynamic characteristic model 100 and a basic schedule 'creation function 300.
This is for setting.

基本スケジュール作成機能300は、スケジュール最適
化機能400における最適値求解計算で良好な収束性を
得るために初期スケジュールとしての基本スケジュール
310を決定し、動特性モデル100に設定するための
ものである。
The basic schedule creation function 300 is for determining a basic schedule 310 as an initial schedule and setting it in the dynamic characteristic model 100 in order to obtain good convergence in the optimum value calculation in the schedule optimization function 400.

スケジュール最適化機能400は、動特性モデル100
を用いて起動特性を予測し、その結果に応じてスケジュ
ールを修正し、再び動特性モデル100に起動スケジュ
ール410を設定し、起動特性を予測するといった繰返
し計算手法で最適スケジュール420を決定するもので
ある。
The schedule optimization function 400 uses the dynamic characteristic model 100
The optimal schedule 420 is determined by an iterative calculation method in which the startup characteristics are predicted using the Dynamic Characteristics Model 100, the schedule is corrected according to the results, the startup schedule 410 is set in the dynamic characteristics model 100 again, and the startup characteristics are predicted. be.

ファジィ推論機能500は、前記、初期値予測機能20
0.基本スケジュール作成機能300及びスケジュール
最適化機能400に対して作用し、それぞれの処理にお
いて、熟練運転員の思考方法を模擬することにより、起
動特性の予測精度の向上と最適スケジュール求解の高速
化を図るもので基本スケジュール作成ルール640.ス
ケジュール最適化ルール670を準備し、これを知識ベ
ースとしている。
The fuzzy inference function 500 includes the initial value prediction function 20 described above.
0. It acts on the basic schedule creation function 300 and the schedule optimization function 400, and in each process, by simulating the thinking method of a skilled operator, it aims to improve the prediction accuracy of startup characteristics and speed up the solution of the optimal schedule. Basic schedule creation rules 640. A schedule optimization rule 670 is prepared and used as a knowledge base.

ここで、スケジュール最適化機能400による起動スケ
ジュール最適化の基本的考え方を第2図を用いて説明し
ておく。
Here, the basic concept of startup schedule optimization by the schedule optimization function 400 will be explained using FIG. 2.

第2図の中で破線はスケジュール最適化前即ち基本スジ
ュールに従ってプラントを起動した場合のタービン応力
、起動パターン及び起動時間を示している。また、実線
はスケジュール最適化後のそれぞれを示している0本図
では、中給から併入時刻(発電機を電力系統に接続する
時刻)が指定された場合を示すものであるが、起動完了
時刻を指定された場合でも本発明は基本的に同じ方法ス
ケジュール最適化が可能である、最適化前の起動スケジ
ュールに従ってプラントを起動した場合、図のように、
タービン応力は起動の前半で制限値に対して大きな余裕
があり、後半では余裕が小さくなり部分的に制限値以上
の応力が発生している。
In FIG. 2, the broken lines indicate the turbine stress, startup pattern, and startup time before schedule optimization, that is, when the plant is started according to the basic schedule. In addition, the solid lines indicate the results after schedule optimization. In this figure, the connection time (the time to connect the generator to the power grid) is specified by the central supply, but the start-up is completed. Even if the time is specified, the present invention can basically optimize the schedule using the same method.If the plant is started according to the startup schedule before optimization, as shown in the figure,
Turbine stress has a large margin with respect to the limit value in the first half of startup, and in the second half, the margin becomes smaller and stress exceeding the limit value is generated in some parts.

この様な起動特性が前記動特性モデル100で予測され
ると、スケジュール最適化機能400ではスケジュール
最適化ルール670を用いたファジィ推論機能500を
動作させ、起動の前半では点火時刻を遅らせたり、ター
ビンの速度保持時間及び負荷保持時間を短縮するなどし
て起動時間の短縮を図っている。また、起動の後半では
負荷保持時間を延長することによりタービン応力の緩和
を図っている。このように1本発明を適用した起動装置
を用いると、運転制限要因であるタービン応力を制限値
以下に抑え、かつ最短時間で中給から指定された時刻通
りに起動できることができる。
When such startup characteristics are predicted by the dynamic characteristics model 100, the schedule optimization function 400 operates the fuzzy inference function 500 using the schedule optimization rules 670, and in the first half of startup, the ignition time is delayed or the turbine The startup time is reduced by shortening the speed holding time and load holding time. Additionally, in the latter half of startup, the load holding time is extended to alleviate turbine stress. As described above, by using the starting device to which the present invention is applied, it is possible to suppress the turbine stress, which is a factor limiting operation, to a limit value or less, and to start the turbine at the designated time from intermediate supply in the shortest possible time.

以下、上記に概説した起動スケジュール作成機能100
0について、各構成機能毎に詳細を説明する。
Below, 100 startup schedule creation functions outlined above
0 will be explained in detail for each constituent function.

(1)初期値予測機能200 動特性モデル100を用いて起動特性を高精度で予測す
るためには、動特性モデルで用いるプラント初期値(ボ
イラ点火時における値で、主に温度状fi)を高精度に
予測する必要がある。しかし、−現時点での計算値に基
づいて冷却特性を解析時に、′、ンノ、′) l!−求゛めることは、停止時の操作内容や現在プラン
トが置かれた環境などに影響されるため困難である。
(1) Initial value prediction function 200 In order to predict startup characteristics with high accuracy using the dynamic characteristics model 100, the plant initial values (values at the time of boiler ignition, mainly temperature condition fi) used in the dynamic characteristics model are It is necessary to predict with high accuracy. However, when analyzing the cooling characteristics based on the current calculated values, ′, nno, ′) l! -Determining this is difficult because it is affected by the operation details at the time of shutdown and the current environment of the plant.

ところが、豊富な運転経験をもつ運転員は、現時点での
温度状態が標準的な状態からどの程度ずれているかによ
って、将来のずれをかなり的確に予測できる0本初期値
予測機能200は、この点に着目したもので、標準冷却
特性で予測した値を、上記現時点でのずれに応じて、ど
の程度修正するかをルール化した運転員の知識を利用す
るものである。以下、本初期値予測機能200に関して
具体的に説明する。
However, operators with extensive operating experience can fairly accurately predict future deviations based on how much the current temperature condition deviates from the standard condition. This method uses the operator's knowledge to create rules for how much the values predicted by the standard cooling characteristics should be corrected in response to the current deviation. The present initial value prediction function 200 will be specifically explained below.

前回の解列時刻(発電機を電力系統から切り離した時刻
)tprと中給からの併入指定時刻tptとの差を停止
時間Δtriとする。この停止時間に応じて予め規定さ
れた4つの起動モード毎に準備された標準スケジュール
の中から該当する起動モードを決定し、起動スケジュー
ルを選択する。ここでいう起動スケジュールとは、ボイ
ラ点火がら起動完了までの起動パターンを規定する次の
パラメータで定義する。
The difference between the previous disconnection time (the time when the generator was disconnected from the power system) tpr and the specified time tpt for parallel connection from the intermediate supply is defined as the stop time Δtri. A corresponding startup mode is determined from standard schedules prepared for each of four startup modes defined in advance according to this stop time, and a startup schedule is selected. The startup schedule here is defined by the following parameters that define the startup pattern from boiler ignition to completion of startup.

(a)ボイラ起動所要時間(点火から通気まで)(b)
第1速度保持時間 (1000rp+mにおいて)(C
)第2速度保持時間 (2800rpmにおいて)(d
)第3速度保持時間 (3600rp■において)(a
)第1負荷保持時間 (初負荷レベルにおいて)(f)
第2負荷保持時間 (20%負荷レベルにおいて) (g)第3負荷保持時間 (40%負荷レベルにおいて
) また、タービン昇速率(rpm/分)及び負荷変化率(
%/分)も起動モードに対応して一義的に決定される。
(a) Boiler startup time (from ignition to ventilation) (b)
1st speed holding time (at 1000 rpm+m) (C
) 2nd speed holding time (at 2800 rpm) (d
) 3rd speed holding time (at 3600 rpm) (a
) 1st load holding time (at initial load level) (f)
Second load holding time (at 20% load level) (g) Third load holding time (at 40% load level) Also, turbine speed increase rate (rpm/min) and load change rate (
%/min) is also uniquely determined in accordance with the startup mode.

上記スケジュールパラメータが決まると、併入時刻から
の逆算により、ボイラ点火時刻traが決定される0次
に、起動モード毎に準備した標準初期値(ボイラ点火時
相当)の中から、今回選択された起動モードにより該当
する初期値を選択する。ここで使用する初期値は次のプ
ラント状態に関するものとする。
Once the above schedule parameters are determined, the boiler ignition time tra is determined by back calculation from the joining time. Select the appropriate initial value depending on the startup mode. The initial values used here relate to the following plant conditions.

(a)  ドラム温度 (b)過熱器出口蒸気温度 (C)再熱器出口蒸気温度 (d)主蒸気管メタル温度 (θ)再熱蒸気管メタル温度 (f)氷壁入口内部流体温度 (g)節炭器出口内部流体温度 (h)節炭器入口内部流体温度 (i)高圧タービン(HPT)第1段後メタル温度(j
)中圧タービン(IPT)第1段後メタル温度(k) 
 ドラム圧力 一方、停止中のタービン及びボイラの標準冷却特性を用
いて、上記状態量に関する現在即ち起動スケジュール作
成時の標準的な値を求める。これを現在標準値推定機能
と呼ぶことにする。なお、上記標準冷却特性を次式で定
義する。
(a) Drum temperature (b) Superheater outlet steam temperature (C) Reheater outlet steam temperature (d) Main steam pipe metal temperature (θ) Reheat steam pipe metal temperature (f) Ice wall inlet internal fluid temperature (g) Economizer outlet internal fluid temperature (h) Economizer inlet internal fluid temperature (i) High pressure turbine (HPT) 1st stage post metal temperature (j
) Intermediate pressure turbine (IPT) first stage after metal temperature (k)
Drum Pressure On the other hand, standard values for the above-mentioned state quantities at present, that is, at the time of creation of the start-up schedule, are determined using the standard cooling characteristics of the turbine and boiler while they are stopped. This will be called the current standard value estimation function. Note that the above standard cooling characteristic is defined by the following equation.

−互肌 T巳TC T工(TPFT−T^)θ    +T^   ・・・
(1)ここで、TPFT ニブラント停止時の標準温度
(℃) T^  :大気温度(”C) Tlsf!:解列後経過時間(分) 丁目Tc:冷却時定数(分) 上記ドラム圧力に関しても同様である。
- Mutual skin T TC T engineering (TPFT-T^) θ +T^ ・・・
(1) Here, TPFT Standard temperature when nibrant is stopped (℃) T^: Atmospheric temperature (''C) Tlsf!: Elapsed time after line separation (minutes) Chome Tc: Cooling time constant (minutes) Regarding the drum pressure mentioned above. The same is true.

また、同様にしてボイラ点火時の初期状態を予測する。In addition, the initial state at the time of boiler ignition is similarly predicted.

この予測を高精度に行なうために、現在状態の実測値2
30(第1図参照)と(1)式による標準値との差を考
慮したファジィ推論により上記予測値を修正する。
In order to perform this prediction with high accuracy, the actual measured value of the current state 2
The predicted value is corrected by fuzzy inference taking into consideration the difference between the standard value obtained by Equation (1) and the standard value obtained by Equation (1).

以下、この修正方法を説明する。This correction method will be explained below.

いま、現在値偏差E(1)、E(2)、・・・・・・E
(11)は、それぞれ、現時点(to )で実測したド
ラム温度X1(toL過熱器出口蒸気温度X5(to)
、・・・・・・ドラム圧力Xxx(to)と(1)式か
ら得られたその標準値Xzs(t ode Xts(t
 (1)、・++++ Xtns(t o)との差を標
準値で正規化したものであり、次式で定義する。
Now, the current value deviation E(1), E(2),...E
(11) is the drum temperature X1 (toL superheater outlet steam temperature X5 (to)) actually measured at the present moment (to), respectively.
,...Drum pressure Xxx(to) and its standard value Xzs(tode Xts(t) obtained from equation (1)
(1), +++++ This is the difference from Xtns(t o) normalized by a standard value, and is defined by the following formula.

第3図は、上記現在値偏差の大きさを定性的に評価する
ためのメンバーシップ関数である0図中の形を規定する
定数であり、PB、PM、PS。
FIG. 3 shows constants that define the shape in Figure 0, which is a membership function for qualitatively evaluating the magnitude of the current value deviation, PB, PM, and PS.

ZO,NS、NM、NBは偏差Eの大きさを定性的に評
価するためにメンバーシップ関数に与えた名称であり、
それぞれ下記の意味をもつ。
ZO, NS, NM, and NB are the names given to the membership functions to qualitatively evaluate the magnitude of the deviation E,
Each has the following meaning.

P B : Po5itive BigP M : P
o5itive MediumP S : Po5it
ive SmallZ O: Zar。
PB: Po5itive BigP M: P
o5itive MediumPS: Po5it
ive SmallZ O: Zar.

N S : Nagatlve SmallNM : 
Negative Mediu@N B  : Neg
ative Bigまた、図の縦軸はメンバーシップ値
である。このメンバーシップ関数を用いて、前記11個
の状態量に関する現在値偏差E(i)(i=1〜11)
を定性的に評価する。
NS: Nagatlve Small NM:
Negative Mediu@N B: Neg
active Big Also, the vertical axis of the figure is the membership value. Using this membership function, the current value deviation E(i) (i=1 to 11) regarding the 11 state quantities is calculated.
Evaluate qualitatively.

第4図は、現在値偏差がボイラ点火時初期に及ぼす影響
度を定性的に整理したものであり、これに従って作成し
た初期値予測ルールの一例を第5E(3)の定性的関係
において、ボイラ点火時過熱器出口蒸気温度の予測値(
(1)式で予測した値)の修正量DtIa(3)を定性
的に決定するためのものである6例えば、ルールNa3
2の場合はI F (E(2) is NS and 
E(3) is NM)THEN (DTIQ(2)i
s NS)という意味である。
Figure 4 qualitatively organizes the degree of influence that the current value deviation has on the initial stage of boiler ignition. Predicted value of superheater outlet steam temperature at ignition (
For example, rule Na3 is for qualitatively determining the correction amount DtIa(3) of
In the case of 2, I F (E(2) is NS and
E(3) is NM)THEN(DTIQ(2)i
s NS).

第6図は、定性的に決定された初期値修正量を定量的な
値に変換するためのメンバーシップ関数である0図中の
EIMB(i)(i=1〜7)はメンバーシップ関数の
形を規定する定数である・PB。
Figure 6 shows the membership function for converting the qualitatively determined initial value correction amount into a quantitative value.・PB is a constant that defines the shape.

PM、PS、 zO,NS、NM、NBjt修正量の大
きさを定性的に表わすためにメンバーシップ関数に与え
た名称であり、第5図の中で使用している名称に対応し
ている。また、図の縦軸はメンバーシップ値である。適
用されたルールにより、初期値修正量が定性的にどのメ
ンバーシップ関数に属するか決定される。成る初期値に
対する修正貴重平均値をもって実際の定量的修正量51
0(第1図参照)とする、これにより各初期値を予測し
たことになる。
PM, PS, zO, NS, NM, NBjt These are the names given to the membership functions to qualitatively represent the magnitude of the correction amount, and correspond to the names used in FIG. Moreover, the vertical axis of the figure is the membership value. The applied rule qualitatively determines which membership function the initial value modification amount belongs to. The actual quantitative correction amount 51 is calculated based on the corrected valuable average value for the initial value.
0 (see FIG. 1), thereby predicting each initial value.

(2)基本スケジュール作成機能300基本スケジュー
ルとはスケジュール最適化の収束計算の初期スケジュー
ルであり、良い収束性を得るには、できるだけ最適値の
近傍に設定するのが望ましい、豊富な運転経験をもつ運
転員は、ボイラ点火時の初期状態に応じて起動スケジュ
ールをかなり的確に決定することができる0本基本スケ
ジュール作成機能300は、この点に着目したもので、
初期状態の標準値からのずれに応じて予め準備した標準
スケジュールのパラメータをどの程度修正するかをルー
ル化した運転員の知識を利用するものである。以下、本
基本スケジュール作成機能300について具体的に説明
する。
(2) Basic schedule creation function 300 The basic schedule is the initial schedule for the convergence calculation of schedule optimization, and to obtain good convergence, it is desirable to set it as close to the optimal value as possible. The 0-line basic schedule creation function 300, which allows the operator to fairly accurately determine the startup schedule according to the initial state at the time of boiler ignition, focuses on this point.
It utilizes the knowledge of operators who have established rules for how much to modify the parameters of a standard schedule prepared in advance in response to deviations from standard values in the initial state. The basic schedule creation function 300 will be specifically explained below.

起動モードに対応して選択された前記起動スケジュール
は標準スケジュールの中から選ばれたちの′であり、必
ずしも今回の起動条件にマツチした表 工゛モのではない、そこで、前記標準初期値と上記力、
′4.。
The startup schedule selected in accordance with the startup mode is one selected from standard schedules, and is not necessarily the one that matches the current startup conditions. Therefore, the standard initial value and the above Power,
'4. .

1法で予測された初期値との差を考慮してスケジュール
を修正するのが本機能である。この修正量は上記予測値
の差に応じてファジィ推論により決定する。
This function corrects the schedule by taking into account the difference from the initial value predicted by the first method. This amount of correction is determined by fuzzy inference according to the difference between the predicted values.

以下、その修正量決定方法について説明する。The method for determining the amount of correction will be explained below.

いま、点火時予測偏差Ep(1)、Ep(2)、・・・
・・・Ep(11)は、それぞれ、初期値予測機能20
0で予測したボイラ点火時のドラム温度Xt(txoL
過熱器出口蒸気温度Xz(tto) *・・・・・・ド
ラム圧力Xzt(two)とその標準値Xt(ttaL
 X5(tXo)。
Now, the predicted deviation at the time of ignition Ep(1), Ep(2),...
... Ep (11) is the initial value prediction function 20, respectively.
The drum temperature Xt (txoL) at the time of boiler ignition predicted at
Superheater outlet steam temperature Xz (tto) *...Drum pressure Xzt (two) and its standard value Xt (ttaL
X5 (tXo).

・・・・・・Xzi(tto)との差を標準値で正規化
したものであり、次式で定義する。
. . . The difference from Xzi(tto) is normalized by a standard value, and is defined by the following formula.

第7図は、上記点火時予測偏差を定性的に評価するため
のメンバーシップ関数である0図中のEpxn(i)(
i=1〜7)はメンバーシップ関数の形を規定する定数
であり、PR,PM、PS。
FIG. 7 shows the membership function Epxn(i)(
i=1 to 7) are constants that define the form of membership functions, PR, PM, and PS.

ZO,NS、NM、NBは偏差Epの大きさを定性的に
評価するためにメンバーシップ関数に与えた名称であり
、意味するとろころは初期値予測機能200におけるも
のと同じである。また、図の縦軸はメンバーシップ値で
ある。
ZO, NS, NM, and NB are names given to membership functions in order to qualitatively evaluate the magnitude of the deviation Ep, and their meanings are the same as those in the initial value prediction function 200. Moreover, the vertical axis of the figure is the membership value.

第8図は、予測したボイラ点火時の初期値に応じて適切
な基本スケジュールを作成するための知識を整理したも
のであり、点火時予測偏差に応じて、どのスケジュール
パラメータを修正するのが効果的であるかを整理したも
のである。これに従つて作成した基本スケジュール作成
ルールの一例を第9図に示す0本図は、ボイラ点火時ド
ラム温度を偏差Hp(1)からボイラ起動所要時間修正
量Dp(1)、第1速度保持時間修正量Dp(2)、・
・・・・・第3負荷保持時間修正量DP(7)を定性的
に決定するためのものである0例えば、ルール&6の場
合は。
Figure 8 summarizes the knowledge needed to create an appropriate basic schedule according to the predicted initial value at the time of boiler ignition, and which schedule parameters are effective to modify depending on the predicted deviation at the time of ignition. This is an organized list of the targets. An example of the basic schedule creation rule created in accordance with this is shown in Figure 9. This figure shows how the drum temperature at boiler ignition is deviated from the deviation Hp (1), the boiler start time required correction amount Dp (1), and the first speed maintained. Time correction amount Dp(2),・
...0 for qualitatively determining the third load holding time correction amount DP(7) For example, in the case of rule &6.

I F(Ep(1)ls PM) THI N(Dp(1) is PM and Dp(
2) is P 5and Dp(3) 1s P S
 and Dp(4) is P 5and Dp(5
) is P S and Dp(6) is P 5
and Dp(7) is P S) という意味である。
I F(Ep(1)ls PM) THI N(Dp(1) is PM and Dp(
2) is P 5and Dp(3) 1s P S
and Dp(4) is P 5and Dp(5
) is P S and Dp(6) is P 5
and Dp(7) is P S).

第10図は、定性的に決定されたスケジュールパラメー
タ修正量を定量的な値に変換するためのメンバーシップ
関数である0図中のDpMa(i)(i=1〜7)はメ
ンバーシップ関数の形を規定する定数である。PB、P
M、PS、 zO,NS。
Figure 10 shows the membership function for converting the qualitatively determined schedule parameter modification amount into a quantitative value. It is a constant that defines the shape. P.B., P.
M, PS, zO, NS.

NM、NBは修正量の大きさを定性的に表わすためにメ
ンバーシップ関数に与えた名称であり、第7図で使用し
ている名称に対応している。また、図の縦軸はメンバー
シップ値である。適用されたルールにより、スケジュー
ルパラメータの修正量が定性的にどのメンバーシップ関
数に属するかが決定される。成るスケジュールパラメー
タに対する修正値が複数のルールによる複数のメンバー
シップで規定された場合は1、各メンバーシップ値に応
じた加重平均値をもって実際の定量的修正量520(第
1図参照)とする、これにより基本スケジュールが作成
されたことになる。
NM and NB are names given to membership functions to qualitatively represent the magnitude of the correction amount, and correspond to the names used in FIG. Moreover, the vertical axis of the figure is the membership value. The applied rule qualitatively determines which membership function the modification amount of the schedule parameter belongs to. If the modification value for the schedule parameter consisting of 1 is defined by multiple memberships according to multiple rules, the weighted average value corresponding to each membership value is set as the actual quantitative modification amount 520 (see Figure 1). This means that a basic schedule has been created.

(3)スケジュール最適化機能400 上記基本スケジュール作成機能300により作成された
基本スケジュールを基にして、最適スケジュール即ちボ
イラ点火から起動完了までの全起動過程において運転制
限条件を満足し、かつ起動所要時間を最小とする前記ス
ケジュールパラメータを決定するのが本スケジュール最
適化機能400である。第11図に本機能の全体処理手
順を示す。
(3) Schedule optimization function 400 Based on the basic schedule created by the above-mentioned basic schedule creation function 300, the optimal schedule is one that satisfies the operation restriction conditions in the entire startup process from boiler ignition to completion of startup, and the required startup time. This schedule optimization function 400 determines the schedule parameter that minimizes . FIG. 11 shows the overall processing procedure of this function.

本機能では、運転制限条件であるタービン熱応力の起動
特性を評価するためにプラントの動特性モデル100を
用いる。起動時の照応カバターンとスケジュールパラメ
ータは大きな相関があり、熱応力の制限値に対する余裕
(以下、マージンと呼ぶ)が小さい程短時間の起動が可
能となる。しかし、コンプレックス法などを利用した従
来の制約条件付非線形最適化アルゴリズムを用いると。
This function uses the plant dynamic characteristic model 100 in order to evaluate the starting characteristics of the turbine thermal stress, which is an operation limiting condition. There is a strong correlation between the anaphoric cover turn at startup and the schedule parameters, and the smaller the margin (hereinafter referred to as margin) for the thermal stress limit value, the shorter the startup time. However, when using conventional constrained nonlinear optimization algorithms such as the complex method.

最適解(最適スケジュール)を得るのに、本例のように
7変数がパラメータとなる場合、少なくとも100回程
度の繰返し計算(動特性モデルによる起動特性の計算)
が必要となり、現実的でない。
In order to obtain the optimal solution (optimal schedule), if seven variables are used as parameters as in this example, it is necessary to repeat calculations at least 100 times (calculation of startup characteristics using a dynamic characteristics model).
is required, which is not realistic.

そこで、ファジィ推論を応用した最適化アルゴリズムと
することにより収束性の大幅な改警を図る。
Therefore, we aim to significantly improve convergence by using an optimization algorithm that applies fuzzy inference.

プラント運転員は、起動前に熱応力特性の予測値を与え
られた場合、そのマージンに応じてどのパラメータをど
の程度短縮できるか経験的に知っている。この経験的か
つ定性的な知識を活用して最適化のためのスケジュール
修正量を決定する。
Plant operators know from experience which parameters can be shortened by how much, depending on the margin, given predicted values of thermal stress characteristics before startup. This empirical and qualitative knowledge is utilized to determine the amount of schedule modification for optimization.

具体的には、まず第11図に示すように、基本スケジュ
ール作成機能300から与えられたスケジュールパラメ
ータpt(i=1〜7)で定まる基本スケジュールに従
ってプラントを起動した場合のタービン熱応力特性を動
特性モデル100を用いて予測する。ここで、動特性モ
デル100は第12図に示すように、スケジュールパラ
メータが与えられた場合、ボイラ点火指令、タービン速
度及び負荷の目標値を計算するためのスケジュール計算
機能110と、ボイラの起動特性を計算するためのボイ
ラモデル120と、ボイラから発生する蒸気条件を受け
てタービンの熱応力を計算するためのタービンモデル1
30から構成されている。
Specifically, as shown in FIG. 11, first, the turbine thermal stress characteristics are controlled when the plant is started according to the basic schedule determined by the schedule parameters pt (i=1 to 7) given from the basic schedule creation function 300. Prediction is made using the characteristic model 100. Here, as shown in FIG. 12, the dynamic characteristic model 100 includes a schedule calculation function 110 for calculating target values of a boiler ignition command, turbine speed, and load when schedule parameters are given, and a boiler starting characteristic. A boiler model 120 for calculating , and a turbine model 1 for calculating thermal stress of the turbine in response to steam conditions generated from the boiler.
It consists of 30.

ここで計算するタービン熱応力は、高圧タービンのロー
タ表面応力、同ボア応力、中圧タービンのロータ表面応
力、同ボア応力の4個所であり、いずれもタービン起動
時に着目すべき重要な運転制限要因である。第11図に
示す起動特性評価機能140は、起動過程を7つの区間
に分割し、各区間毎の最小応力マージンma  (j=
1〜7)を求める0本図ではm4を区間jにおける高圧
ロータ表面応力と中圧ロータ表面応力のうち最小マージ
ンMs(j)と高圧ロータボア応力と中圧ロータボア応
力のうち最小マージンMa(j)の両者の意味で示して
いる0本起動特性評価機能140は、次に続く応力マー
ジン評価機能150で熱溶カバターンの特徴を抽出する
ためのものである。
The turbine thermal stress calculated here is from four areas: rotor surface stress and bore stress of the high-pressure turbine, and rotor surface stress and bore stress of the intermediate-pressure turbine, all of which are important operational limiting factors that should be noted when starting the turbine. It is. The starting characteristic evaluation function 140 shown in FIG. 11 divides the starting process into seven sections, and calculates the minimum stress margin ma (j=
1 to 7) In this figure, m4 is the minimum margin Ms(j) between high pressure rotor surface stress and medium pressure rotor surface stress in section j, and the minimum margin Ma(j) between high pressure rotor bore stress and medium pressure rotor bore stress. The zero starting characteristic evaluation function 140 shown in both senses is for extracting the characteristics of the hot melt cover turn in the subsequent stress margin evaluation function 150.

応力マージン評価機能150では、第13図に示すメン
バーシップ関数を用いて、ロータ表面応カマージンMs
(j)及びロータボア応カマージンMa(j)を定性的
に評価することにより、熱溶カバターンの特徴を抽出す
る1図中のMMB (i )(i、=1〜6)はメンバ
ーシップ関数の形を規定する’NM、NBは応力マージ
ンMs(j)及びMB(j )の大きさを定性的に評価
するためにメンバーシップ関数に与えた名称である。ま
た、図の縦軸はメンバーシップ値である。
The stress margin evaluation function 150 uses the membership function shown in FIG.
MMB (i) (i, = 1 to 6) in Figure 1 is a form of membership function. NM and NB are the names given to the membership functions in order to qualitatively evaluate the magnitudes of the stress margins Ms(j) and MB(j). Moreover, the vertical axis of the figure is the membership value.

スケジュール最適化ルール670は、「熱溶カバターン
がどの様であればどのスケジュールパラメータをどの程
度修正すれば良い」といった断片的知識である。第14
図は、動特性モデル100を用いて予測したタービン熱
溶カバターンに応じてスケジュールを修正するための知
識を整理したものである0本図は、起動過程の各区間に
おける最小マージンに応じて、どのスケジュールパラメ
ータを修正するのが効果的であるかを整理したものであ
る。ここで、Ms(I L Ms(2)、 ・・−・M
s(7)及びMa(1)、Ma(2)、−・”MB(7
)は、それぞれ第1.第2.・・・・・・第7区間にお
けるロータ表面最小応力マージン及びロータボア最小応
力マージンである。
The schedule optimization rule 670 is piecemeal knowledge such as ``What schedule parameters should be modified and how much should the hot melt cover pattern be?''. 14th
The figure summarizes the knowledge for modifying the schedule according to the turbine hot melt cover turn predicted using the dynamic characteristic model 100. This is a summary of whether it is effective to modify schedule parameters. Here, Ms(I L Ms(2), . . . M
s(7) and Ma(1), Ma(2), -・”MB(7
) are the first. Second. . . . The minimum stress margin on the rotor surface and the minimum stress margin on the rotor bore in the seventh section.

6第15図は、上記考え方に従って作成したスケ゛1°
?ニー)LtR適化ルールの一例を示すものである。
6 Figure 15 shows the scale 1 degree created according to the above idea.
? 2) shows an example of the LtR optimization rule.

本図は、ロータ表面応力に関して、第5区間最小応力マ
ージンMs(5)と第6区間最小応力マージンMs(6
)からスケジュールパラメータの修正量としての第3速
度保持時間修正量を量DPT(4)1第1負荷保持時間
修正量DPT(5)及び第2負荷保持時間修正量DPT
(6)を定性的に決定するためのものである0例えば、
ルールNQ54の場合は、I F (Ms(5) is
 P B and Ms(6) is PM)THEN
 (DPT(4) is NM and Dpt(5)
isNM  and  DPT(6)is  N Sと
いう意味である。
This figure shows the 5th section minimum stress margin Ms(5) and the 6th section minimum stress margin Ms(6) regarding the rotor surface stress.
), calculate the third speed holding time correction amount as the schedule parameter correction amount DPT(4)1 the first load holding time correction amount DPT(5) and the second load holding time correction amount DPT
(6) is for qualitatively determining 0. For example,
In the case of rule NQ54, I F (Ms(5) is
P B and Ms(6) is PM) THEN
(DPT(4) is NM and Dpt(5)
isNM and DPT(6)isNS.

第16図は、定性的に決定されたスケジュールパラメー
タ修正量を定量に変換するためのメンバーシップ関数で
ある0図中のDPTMB (i )  (i =1〜7
)はメンバーシップ関数の形を規定する定数である。P
B、PM、PS、20.NS、NM。
FIG. 16 shows DPTMB (i) (i = 1 to 7) in FIG.
) is a constant that defines the form of the membership function. P
B, PM, PS, 20. N.S., N.M.

NBは修正量の大きさを定性的に表わすためにメンバー
シップ関数に与えた名称である。また、図の縦軸はメン
バーシップ値である。
NB is a name given to the membership function to qualitatively represent the magnitude of the correction amount. Moreover, the vertical axis of the figure is the membership value.

第11図に示す修正パラメータ選択機能160は、応力
マージン評価機能150により抽出された応カバターン
の特徴とスケジュール最適化ルール670とを照合させ
ることにより修正パラメータを選択し、適用されたルー
ルによりパラメータの修正量が定性的に第16図のどの
メンバーシップ関数に属するかが決定される。
The modified parameter selection function 160 shown in FIG. 11 selects modified parameters by comparing the characteristics of the response pattern extracted by the stress margin evaluation function 150 with the schedule optimization rules 670, and selects the modified parameters by comparing the characteristics of the response pattern extracted by the stress margin evaluation function 150 with the schedule optimization rules 670. It is qualitatively determined to which membership function in FIG. 16 the modification amount belongs.

成るスケジュールパラメータに対する修正量が複数のル
ールによる複数のメンバーシップで規定された場合は、
各メンバーシップ値に応じた加重平均値をもって実際の
定量的修正量530(第1図参照)とする、これを行な
うのが第11図に示すスケジュール修正量決定機能17
0である。
If the amount of modification to a schedule parameter that consists of is specified by multiple memberships in multiple rules,
The weighted average value corresponding to each membership value is used as the actual quantitative modification amount 530 (see FIG. 1).The schedule modification amount determination function 17 shown in FIG. 11 performs this.
It is 0.

以上でスケジュールパラメータの修正量が決定され、修
正されたスケジュールに従って、プラントを起動した場
合の起動特性を予測するために、再び動特性モデル10
0を動作させる0以上を繰返すことにより最適起動スケ
ジュール420(第1図参照)を求めることができる。
The amount of modification of the schedule parameters is determined above, and the dynamic characteristic model 1 is used again to predict the startup characteristics when the plant is started according to the modified schedule.
The optimal activation schedule 420 (see FIG. 1) can be obtained by repeating 0 or more operations in which 0 is activated.

なお、第11図中の収束判定機能180は、上記繰返し
計算に一1袢いて作成される起動スケジュールの最適性
を判・/゛ 2定するためのものである。また、その判定基準は、全
起動過程においてタービン熱応力が制限値以下となる起
動スケジュールのうち起動所要時間が前回までのものと
比較した時間短縮率が所定値以下となった場合とし、そ
の中で最も起動所要時間が小さくなる起動スケジュール
を最適スケジュール420(第1図参照)とする。
Note that the convergence determination function 180 in FIG. 11 is for determining the optimality of the startup schedule created by adding the above-described iterative calculations. In addition, the criterion for this is that among the startup schedules in which the turbine thermal stress is below the limit value during the entire startup process, the time reduction rate compared to the previous startup schedule is below a predetermined value; The startup schedule that requires the shortest startup time is defined as the optimal schedule 420 (see FIG. 1).

ここで決定された最適スケジュール420は、第1図に
示したように、ユーザインターフェイス4000を介し
て表示装置115000に表示されると共に、起動制御
機能3000に設定される。起動制御機能3000に設
定された最適スケジュール420は運転員6000から
の実行指令430を受けて実行される。起動制御機能3
000は。
The optimal schedule 420 determined here is displayed on the display device 115000 via the user interface 4000 and set in the activation control function 3000, as shown in FIG. The optimal schedule 420 set in the startup control function 3000 is executed upon receiving an execution command 430 from the operator 6000. Startup control function 3
000 is.

これを実行するためにプラント2000からのプロセス
入力3010を受けて、制御出力3010をプラントに
与える。
To do this, it receives process input 3010 from plant 2000 and provides control output 3010 to the plant.

運転員6000は、知識ベース600の追加。Operator 6000 adds knowledge base 600.

変更、削除が必要となった場合は、ユーザインタ第17
図は、ファジィ推論を応用したスケジュール最適化機能
400におけるスケジュール最適化過程を示すものであ
る0図中の番号は繰返し計算の第1.3,5.200回
目作成された起動スケジュールとそのときの起動特性を
示す、ここで、第1回目は基本スケジュールに対応し、
200回目最適スケジュールに対応する。本図かられか
るように、5回目でほぼ最適値に近いスケジュールが得
られており1本アルゴリズムによる収束性は極めて良好
である。また、作成された起動スケジュールは併入指定
時刻(本図では解列時刻からの経過時間で示してあり4
80分(8時間))通りに起動できることを示している
。タービンの熱応力マージンが大きな起動前半では速度
保持や負荷保持を省略することにより起動時間の短縮す
る方向でスケジュールが修正されており、熱応力マージ
ンが負となる起動後半では負荷保持を延長することによ
り熱応力の緩和を図っている。このよう通りの起動が可
能となる。
If changes or deletions are necessary, please contact User Interface No. 17.
The figure shows the schedule optimization process in the schedule optimization function 400 that applies fuzzy inference. The numbers in the figure indicate the startup schedules created at the 1st, 3rd, 5th, and 200th iterations of repeated calculations, and the schedule optimization process at that time. Showing the startup characteristics, here, the first time corresponds to the basic schedule,
Corresponds to the 200th optimal schedule. As can be seen from this figure, a schedule close to the optimal value was obtained in the fifth time, and the convergence of the single algorithm is extremely good. In addition, the created startup schedule is displayed at the specified merging time (in this figure, it is shown as the elapsed time from the time of train separation).
80 minutes (8 hours)). The schedule has been revised to shorten the startup time by omitting speed holding and load holding in the first half of startup, when the thermal stress margin of the turbine is large, and by extending load holding in the second half of startup, when the thermal stress margin is negative. This is aimed at alleviating thermal stress. It is possible to start up in this way.

スケジュール誤差監視機能700について次に説明する
。スケジュール誤差監視機能700は次の2つの機能で
スケジュール誤差の発生を監視する。
The schedule error monitoring function 700 will be described next. The schedule error monitoring function 700 monitors the occurrence of schedule errors using the following two functions.

1つは、実際にスケジュール誤差が発生したことを監視
する誤差検知機能であり、別の1は、実際に誤差が発生
する前に誤差の発生を予測する誤差予測機能である。
One is an error detection function that monitors whether a schedule error has actually occurred, and the other is an error prediction function that predicts the occurrence of an error before it actually occurs.

誤差検知機能は、起動前に作成された起動スケジュール
パラメータで規定される起動パターンが実際の起動中に
正確に遂行されているか否かを監視するものである。し
たがって、この起動パターンのずれを監視するものであ
る。
The error detection function monitors whether a startup pattern defined by startup schedule parameters created before startup is accurately executed during actual startup. Therefore, this deviation in the activation pattern is monitored.

一方、誤差予測機能は、起動パターンのずれとなって現
われる前に、ずれの原因となるプラント状態量の挙動を
監視しておき、起動スケジュール作成時に動特性モデル
100で予測したプラントのである。
On the other hand, the error prediction function monitors the behavior of the plant state quantity that causes the deviation before it appears as a deviation in the starting pattern, and predicts the behavior of the plant state quantity using the dynamic characteristic model 100 when creating the starting schedule.

上記2つの機能でスケジュール誤差を監視し、誤差発生
を検知あるいは予測できた場合は、その旨710を関連
するプラント状態量と共にユーザインターフェイス40
00を介し1表示袋[5000に表示する。これにより
、運転員600oは速かに対応処置をすることが可能と
なる。
If schedule errors are monitored using the above two functions and the occurrence of an error is detected or predicted, a message to that effect 710 is sent to the user interface 40 along with related plant state quantities.
00 through 1 display bag [Display in 5000. This allows the operator 600o to take prompt action.

(発明の変形例) 以上述べた本発明の実施例では、スケジュールパラメー
タとしてボイラ起動所要時間とタービン速度保持時間及
び負荷保持時間に着目しているが、必ずしもこれらに限
る必要はなく、タービン速度変化率や負荷変化率あるい
はボイラ昇温速度や昇圧速度などプラントの起動パター
ンを代表するパラメータであれば本発明は基本原理を変
えることな〈実施できることは明らかである。
(Modified example of the invention) In the embodiment of the present invention described above, the schedule parameters are focused on the boiler startup time, the turbine speed holding time, and the load holding time, but it is not necessarily limited to these, and the turbine speed change It is clear that the present invention can be implemented without changing the basic principle as long as the parameters are representative of the plant start-up pattern, such as rate, load change rate, boiler temperature increase rate, pressure increase rate, etc.

また、本発明の実施例では、運転制限条件としてタービ
ン熱応力に着目しているが、必ずしもこれだけに限る必
要はなく、タービン入口蒸気温度及びその変化速度、あ
るいはタービンケーシング温度などタービン熱応力を間
接的に推定できる状態量、あるいは、タービン伸び差、
ボイラ蒸気圧力上昇速度、ボイラ燃焼ガス温度など運転
上重要となる制限要因はプラントにより異なるから、必
要に応じて考慮すれば良く、本発明は本質を変えること
な〈実施できることは明らかである。
Further, in the embodiments of the present invention, attention is focused on turbine thermal stress as an operation limiting condition, but it is not necessarily limited to this. The state quantity that can be estimated visually or the turbine expansion difference,
Important limiting factors for operation, such as boiler steam pressure increase rate and boiler combustion gas temperature, vary depending on the plant, so they can be taken into consideration as necessary, and it is clear that the present invention can be practiced without changing its essence.

また、本発明の実施例では中央給電指令所から指定され
る併入時刻を正確にまもった起動スケジュールを作成す
る方法を例に説明したが、本発明を実施するにあたり、
必ずしも併入時刻が指定される必要はなく、プラントに
応じてボイラ点火時刻、タービン通気時刻、目標負荷到
達時刻などであっても本発明の本質が変らないのは明ら
かである。
In addition, in the embodiment of the present invention, a method for creating a startup schedule that accurately adheres to the joining time specified by the central power dispatch center was explained as an example, but in implementing the present invention,
It is clear that the essence of the present invention does not change even if the joining time does not necessarily have to be designated, and the boiler ignition time, turbine ventilation time, target load arrival time, etc. can be specified depending on the plant.

また1本実施例でプラント初期値として着目している状
態は、本発明を実施するプラントの機器構成及び計測位
置に応じて適宜選定すべきもので必ずしも本実施例と同
じものを使う必要のないこ時刻からの時間関数として表
わしているが、ボイラ消火時刻などプラント冷却過程に
移行する時刻を基準とした時間関数で表わしたとしても
、本発明の実施にあたり本質が変わるものではない。
In addition, the conditions focused on as the initial plant values in this example should be selected appropriately depending on the equipment configuration and measurement position of the plant in which the present invention is implemented, and it is not necessarily necessary to use the same values as in this example. Although it is expressed as a time function from time, even if it is expressed as a time function based on the time when the plant cools down, such as the boiler extinguishing time, the essence of the present invention does not change.

また1本実施例のスケジュール最適化機能における収束
判定方法を必ずしも採用する必要はなく、繰返し計算を
所定回数だけ実施し、その中で運転制限条件を満足し、
起動所要時間が最小となる起動スケジュールを最適値と
見做す方法、あるいは運転制限条件を満足する所定個の
起動スケジュール候補が得られるまで繰返し計算を実施
し、その中で起動所要時間が最小となるものを最適値と
見做す方法を採用しても、本発明の本質が変わらないの
は勿論である。
In addition, it is not necessary to necessarily adopt the convergence determination method in the schedule optimization function of this embodiment, but it is necessary to perform repeated calculations a predetermined number of times and satisfy the operation restriction conditions,
A method that considers the startup schedule that minimizes the startup time as the optimal value, or performs repeated calculations until a predetermined number of startup schedule candidates that satisfy the operation restriction conditions are obtained, and then calculates the startup schedule that takes the minimum startup time among them. Of course, the essence of the present invention does not change even if a method is adopted in which the optimum value is regarded as the optimum value.

さらに、本実施例で用いているメンバーシップ関数は全
て三角型としたが、必ずしもこの形に拘わる必要はなく
、プラントの特性及び運転員の知識に応じて、2次曲線
や指数曲線を採用しても、゛本発明の本質が変わるもの
ではない、また、メン/、、 、  シップ関数の形だ
けでなく、その数も任意に設定しても、本発明の本質が
変わるものではない。
Furthermore, although all membership functions used in this example are triangular, they do not necessarily have to be of this shape; quadratic or exponential curves may be adopted depending on the characteristics of the plant and the knowledge of the operator. However, the essence of the present invention does not change even if not only the form of the men/ship functions but also the number of them is arbitrarily set, the essence of the present invention does not change.

〔発明の効果〕〔Effect of the invention〕

本発明によれば、プラントの起動所要時間を従来方式と
比較して約30%短縮することができるため、負荷需要
の変動に伴い発電プラントの頻繁な起動停止が必要とな
る電力系統の安定かつ経済的な運用が可能となる。また
、これにより運転員の負担が大幅に軽減される。また、
起動所要時間の短縮に伴って起動時のエネルギー損失も
約15%低減できるため発電プラントの運用コストを大
幅に低減できる。
According to the present invention, the time required to start up a plant can be reduced by about 30% compared to conventional systems, thereby stabilizing power systems that require frequent starting and stopping of power plants due to fluctuations in load demand. Economical operation becomes possible. This also significantly reduces the burden on the operator. Also,
As the start-up time is reduced, energy loss during start-up can be reduced by about 15%, resulting in a significant reduction in power plant operating costs.

また、本発明によれば起動特性の高精度予測が可能とな
り、運転制限条件および指定時刻を忠実に守った起動が
できるため、プラント運用上の安全性を向上し、電力系
統への正確な電力供給が可能となる。
In addition, according to the present invention, it is possible to predict startup characteristics with high accuracy, and startup can be performed while faithfully observing operating limit conditions and specified times, improving safety in plant operation and providing accurate power to the power system. supply becomes possible.

さらに、本発明によれば、起動スケジュール最適化計算
における計算所要時間を短縮できるため、演算性態の高
い高価な計算機を使用する必要がなく、安価な起動装置
とすることが可能となる。
Further, according to the present invention, the calculation time required for the startup schedule optimization calculation can be shortened, so there is no need to use an expensive computer with high computational performance, and an inexpensive startup device can be provided.

【図面の簡単な説明】[Brief explanation of the drawing]

第1図は発電プラント起動装置の全体構成を示す、第2
図は起動スケジュール最適化の基本的考え方を示す、第
3図は現在値偏差評価用メンバーシップ関数を示す、第
4図は現在値偏差がボイラ点火時初期値へ及ぼす影響度
を示す、第5図は初期予測ルールの一例を示す、第6図
は初期値修正電変換用メンバーシップ関数を示す、第7
図は点火時予測偏差評価用メンバーシップ関数を示す。 第8図は点火時予測偏差と修正対象スケジュールパラメ
ータの関係を示す、第9図は基本スケジュール作成のル
ールの一例を示す、第10図はスケジュールパラメータ
修正用メンバーシップ関数を示す、第11図はスケジュ
ール処理機能における全体処理手順を示す、第12図は
動特性モデルを示す、第13図は応力マージン評価用メ
ンバーシップ関数を示す、第14図は応力マージンと修
正対象スケジュールパラメータの関係を示す、第15図
はスケジュール最適化ルールの一例を示す。 第16図はスケジュールパラメータ修正用メンバ寮30 橘5 の ネηDiall*−rv 二E(Z>、ε(3) ff
らDTIC(Zル&1ジす3ル−JVUF (a(Z)
Is  NS −ε(j)isNM)   TMEN(
DT+4(Z〕Is  NS)寥i 第90 ルールQ吏aj、(ppo、ab*sl)@フgnA 
optD ts ps) 第10 口 不13図
Figure 1 shows the overall configuration of the power plant starting device;
Figure 3 shows the basic concept of startup schedule optimization, Figure 3 shows the membership function for evaluating the current value deviation, Figure 4 shows the degree of influence of the current value deviation on the initial value at the time of boiler ignition, and Figure 5. The figure shows an example of the initial prediction rule. Figure 6 shows the initial value correction electric conversion membership function.
The figure shows the membership function for evaluating the predicted deviation at the time of ignition. Fig. 8 shows the relationship between the predicted ignition deviation and the schedule parameters to be corrected, Fig. 9 shows an example of rules for creating a basic schedule, Fig. 10 shows membership functions for correcting schedule parameters, and Fig. 11 12 shows the dynamic characteristic model; FIG. 13 shows the membership function for stress margin evaluation; FIG. 14 shows the relationship between the stress margin and schedule parameters to be modified. FIG. 15 shows an example of a schedule optimization rule. Figure 16 shows member dormitory 30 Tachibana 5's ηDial*-rv 2E(Z>, ε(3) ff) for modifying schedule parameters.
raDTIC (Zle & 1jis3ru - JVUF (a(Z)
Is NS −ε(j)isNM) TMEN(
DT + 4 (Z] Is NS) 寥i 90th rule Q 吏aj, (ppo, ab*sl) @fugnA
optD ts ps) No. 10 Kuchifu 13

Claims (1)

【特許請求の範囲】 1、燃料の燃焼により発生する熱により蒸気を発生させ
るためのボイラと、発生した蒸気の熱エネルギを機械エ
ネルギに変換するための蒸気タービンと、変換された該
機械エネルギを電気エネルギに変換するための発電機と
で構成される発電プラントにおいて、該発電プラントの
起動前初期温度状態に応じて機器の操作時間及び制御目
標値を決定するための起動スケジュール作成手段を有し
、該初期温度状態を予測するために、プラント停止後の
プラント冷却特性を時間関数で表わした標準冷却特性関
数と、該起動スケジュールを作成する時刻(現在時刻)
におけるプラント温度状態を実測するための現在状態実
測手段と、該標準冷却特性関数を用いて現在時刻のプラ
ント温度状態を推定するための現在状態推定手段と、該
標準冷却特性関数を用いて起動開始時刻のプラント温度
状態を予測するための第1の初期値予測手段と、該現在
状態推定手段にて得られた推定現在値と該現在状態実測
手段にて得られた実測現在値との差である現在値推定誤
差を求めるための現在値推定誤差計算手段と、該第1の
初期値予測手段で得られた予測値を該現在値推定誤差に
応じて修正するための第2の初期値予測手段を有し、該
起動スケジュールを作成するためにプラントの停止時間
の長さに応じて予め定義された複数の起動モードの中か
ら実際のプラント停止時間により該当する起動モードを
選択する起動モード選択手段と、選択された該起動モー
ドに対応して標準スケジュールを選択する標準スケジュ
ール選択手段と、選択された該標準スケジュールに対応
した標準初期温度と実際の初期温度状態との差を求める
ための温度偏差計算手段と、該温度偏差に応じて選択さ
れた該標準スケジュールを修正し、新たなスケジュール
を作成する基本スケジュール作成手段を有することを特
徴とする発電プラント起動装置。 2、特許請求の範囲第1項記載の第2の該初期値予測手
段において、該現在値推定誤差の大きさを定性的に評価
するための現在値推定誤差評価手段と、定性的に評価さ
れた該現在値推定誤差との関係において予め準備した修
正ルールを適用して該第1の初期値予測手段にて得られ
た予測値の修正量を決定するための予測値修正手段を有
し、該基本スケジュールを作成するために、該温度偏差
の大きさを定性的に評価するための温度偏差評価手段と
、定性的に評価された該温度偏差との関係において予め
準備した修正ルールを適用して該選択された標準スケジ
ュールの修正量を決定するための標準スケジュール修正
手段を有することを特徴とする発電プラント起動装置。
[Claims] 1. A boiler for generating steam using heat generated by combustion of fuel, a steam turbine for converting the thermal energy of the generated steam into mechanical energy, and a steam turbine for converting the thermal energy of the generated steam into mechanical energy. A power generation plant comprising a generator for converting into electrical energy, and a power generation plant comprising a startup schedule creation means for determining equipment operation time and control target values according to an initial temperature state before startup of the power generation plant. , in order to predict the initial temperature state, a standard cooling characteristic function that expresses the plant cooling characteristics after the plant is stopped as a time function, and the time (current time) at which the startup schedule is created.
current state measuring means for actually measuring the plant temperature state at the current time; current state estimating means for estimating the plant temperature state at the current time using the standard cooling characteristic function; and start-up starting using the standard cooling characteristic function. The difference between the estimated current value obtained by the first initial value prediction means for predicting the plant temperature state at the time and the current state estimation means and the actual measured current value obtained by the current state actual measurement means. A current value estimation error calculation means for determining a certain current value estimation error, and a second initial value prediction for correcting the predicted value obtained by the first initial value prediction means according to the current value estimation error. a startup mode selection means for selecting a corresponding startup mode according to the actual plant shutdown time from among a plurality of startup modes predefined according to the length of the plant shutdown time in order to create the startup schedule; a standard schedule selection means for selecting a standard schedule corresponding to the selected startup mode; and a temperature for determining a difference between a standard initial temperature corresponding to the selected standard schedule and an actual initial temperature state. A power generation plant starting device characterized by having a deviation calculation means and a basic schedule creation means for modifying the selected standard schedule according to the temperature deviation and creating a new schedule. 2. The second initial value prediction means according to claim 1, further comprising a current value estimation error evaluation means for qualitatively evaluating the magnitude of the current value estimation error; a predicted value correction means for determining a correction amount of the predicted value obtained by the first initial value prediction means by applying a correction rule prepared in advance in relation to the current value estimation error; In order to create the basic schedule, a temperature deviation evaluation means for qualitatively evaluating the magnitude of the temperature deviation and a correction rule prepared in advance in relation to the qualitatively evaluated temperature deviation are applied. 1. A power generation plant starting device comprising standard schedule modification means for determining the amount of modification of the selected standard schedule.
JP23901586A 1986-10-09 1986-10-09 Starting device for power generation plant Pending JPS6394006A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP23901586A JPS6394006A (en) 1986-10-09 1986-10-09 Starting device for power generation plant

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP23901586A JPS6394006A (en) 1986-10-09 1986-10-09 Starting device for power generation plant

Publications (1)

Publication Number Publication Date
JPS6394006A true JPS6394006A (en) 1988-04-25

Family

ID=17038623

Family Applications (1)

Application Number Title Priority Date Filing Date
JP23901586A Pending JPS6394006A (en) 1986-10-09 1986-10-09 Starting device for power generation plant

Country Status (1)

Country Link
JP (1) JPS6394006A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02181001A (en) * 1988-12-29 1990-07-13 Hitachi Ltd Starting operation support system for thermal power plant

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02181001A (en) * 1988-12-29 1990-07-13 Hitachi Ltd Starting operation support system for thermal power plant

Similar Documents

Publication Publication Date Title
US6917838B2 (en) Open-loop and closed-loop control method, and a control device for starting up and shutting down a process component of a technical process
US4868754A (en) Method of starting thermal power plant
JPH0521241B2 (en)
RU2213997C2 (en) Method for evaluating temperature of blades in steam turbine
Arastou et al. Modeling and parameter estimation of a steam power plant including condenser back-pressure uncertainty using operational data
Prasad et al. A novel performance monitoring strategy for economical thermal power plant operation
JPH08339204A (en) Autonomous adaptive optimization control system for thermal power station
JP2003150233A (en) Method and device for evaluating performance deteriorating state for plant
JP3666035B2 (en) Thermal power plant autonomous adaptive control system
JPH09152903A (en) Method and device for generating plant operation paln
Holmes et al. Fuzzy damage-mitigating control of a fossil power plant
JPS6394006A (en) Starting device for power generation plant
JPH11242503A (en) Plant operation control support system
JP2000056805A (en) Predictive control unit
JPS6396401A (en) Power plant starter
JPS6394008A (en) Starting device for power generation plant
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
Zaccaria et al. Adaptive control of micro gas turbine for engine degradation compensation
JPS6394005A (en) Starting device for power generation plant
JPS6394007A (en) Starting device for power generation plant
JPS6394003A (en) Starting device for power generation plant
JPS6394004A (en) Starting device for power generation plant
JPH0573886B2 (en)
JPS6394009A (en) Starting device for power generation plant
JP2678216B2 (en) Thermal power plant startup operation support device