JP2002054460A - Combustion vibration predicting device - Google Patents

Combustion vibration predicting device

Info

Publication number
JP2002054460A
JP2002054460A JP2001022338A JP2001022338A JP2002054460A JP 2002054460 A JP2002054460 A JP 2002054460A JP 2001022338 A JP2001022338 A JP 2001022338A JP 2001022338 A JP2001022338 A JP 2001022338A JP 2002054460 A JP2002054460 A JP 2002054460A
Authority
JP
Japan
Prior art keywords
internal pressure
pressure fluctuation
data
nox
input
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.)
Granted
Application number
JP2001022338A
Other languages
Japanese (ja)
Other versions
JP3999467B2 (en
Inventor
Masumi Nomura
真澄 野村
Katsunori Tanaka
克則 田中
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.)
Mitsubishi Heavy Industries Ltd
Original Assignee
Mitsubishi Heavy Industries 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 Mitsubishi Heavy Industries Ltd filed Critical Mitsubishi Heavy Industries Ltd
Priority to JP2001022338A priority Critical patent/JP3999467B2/en
Priority to US10/057,955 priority patent/US6922612B2/en
Priority to CA002369539A priority patent/CA2369539C/en
Priority to DE60210157T priority patent/DE60210157T2/en
Priority to EP02002276A priority patent/EP1229226B1/en
Publication of JP2002054460A publication Critical patent/JP2002054460A/en
Application granted granted Critical
Publication of JP3999467B2 publication Critical patent/JP3999467B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • Y02T50/677

Abstract

PROBLEM TO BE SOLVED: To facilitate regulation of a combustion control system and early detect the occurrence of combustion vibration during operation by predicting combustion vibration generated by the combustor of a gas turbine by a numeral formula model. SOLUTION: A combustion vibration predicting device comprises an input means 4 to input plant data, weather data, and the limit value of an internal pressure fluctuation; an internal pressure fluctuation characteristics grasping means 1 to effect numeral-formula-model of the internal pressure fluctuation of the combustor from inputted plant data and weather data; a combustion vibration region estimating means 2 to apply a limit value of an internal pressure fluctuation to a numeral formula model determined by the internal pressure fluctuation characteristics grasping means 1; and an output means 5 to output a combustion vibration region estimated result by a combustion vibration region estimating means 2.

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は、燃焼振動予測装置
に関し、特に、事業用や航空機用のガスタービンの燃焼
器で発生する燃焼振動を予測する燃焼振動予測装置に関
するものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a combustion vibration predicting apparatus, and more particularly to a combustion vibration predicting apparatus for predicting combustion vibration generated in a combustor of a commercial or aircraft gas turbine.

【0002】[0002]

【従来の技術】制御系の調整に関しては、従来よりカル
マンフィルタやニューラルネットワークなどによるさま
ざまな取り組みがなされてきた。しかし、ガスタービン
の燃焼振動に関しては、現象が複雑で、現地試運転調整
は、依然として運転員の経験とノウハウに基づき調整さ
れているのが現状である。
2. Description of the Related Art Various approaches have been made to adjust a control system by using a Kalman filter or a neural network. However, the phenomena related to the combustion oscillation of the gas turbine are complicated, and the field trial operation adjustment is still adjusted based on the experience and know-how of the operator.

【0003】運転中の燃焼振動監視に関しては、燃焼器
に設置された圧力センサによりデータを収集し、そのデ
ータより燃焼振動の周波数分析をすることにより、燃焼
振動が異常であることを早期に捉え、健全性を重視して
燃焼器の運転状態を監視する技術が知られている。
[0003] Regarding the monitoring of combustion vibration during operation, data is collected by a pressure sensor installed in the combustor, and the frequency of combustion vibration is analyzed based on the data, so that it is possible to grasp early that combustion vibration is abnormal. There is known a technique for monitoring the operating state of a combustor with an emphasis on soundness.

【0004】図19は従来例として、特開平11−32
4725号公報に示された燃焼異常監視装置(ガスター
ビン異常監視装置)を示している。この監視装置は、ガ
スタービン燃焼器に設置された圧力センサ100と、圧
力センサ100からの検出信号をデジタルデータに変換
して受けるA/D変換器101と、そのデジタルデータ
を周波数成分に分解して解析する周波数解析装置102
と、ガスタービン燃焼器出力及びその燃料供給量で規定
されるパラメータに基づいて監視対象の周波数成分に関
する基準データを可変設定する判定条件設定部103
と、その基準データに基づいて周波数成分の解析データ
から燃焼振動現象に起因する周波数成分を抽出し、その
周波数成分の振幅値と監視対象の周波数成分に関するデ
ータの正常時の振幅値との比較によりガスタービンの燃
焼振動状態を判定する判定処理部104と、この判定結
果に関するデータを表示する結果表示部105とを有し
ている。
FIG. 19 shows a conventional example as disclosed in JP-A-11-32.
1 shows a combustion abnormality monitoring device (gas turbine abnormality monitoring device) disclosed in Japanese Patent No. 4725. This monitoring device includes a pressure sensor 100 installed in a gas turbine combustor, an A / D converter 101 that converts a detection signal from the pressure sensor 100 into digital data and receives it, and decomposes the digital data into frequency components. Analysis device 102 for analyzing
And a determination condition setting unit 103 for variably setting reference data on a frequency component to be monitored based on a parameter defined by a gas turbine combustor output and a fuel supply amount thereof.
And, based on the reference data, extract a frequency component caused by the combustion oscillation phenomenon from the analysis data of the frequency component, and compare the amplitude value of the frequency component with the normal amplitude value of the data on the monitored frequency component. It has a judgment processing unit 104 for judging the combustion oscillation state of the gas turbine, and a result display unit 105 for displaying data relating to this judgment result.

【0005】上述の監視装置では、圧力センサ100よ
りのデータにより燃焼振動の周波数分析をし、周波数毎
の正常時の振動の振幅値と実際値との比較によって燃焼
振動が異常であるか否かを判定することができる。
In the above-described monitoring device, the frequency of the combustion vibration is analyzed based on the data from the pressure sensor 100, and whether the combustion vibration is abnormal is determined by comparing the amplitude value of the normal vibration for each frequency with the actual value. Can be determined.

【0006】[0006]

【発明が解決しようとする課題】しかし、上述のような
従来の監視装置では、燃焼振動の正常時の振幅値を周波
数毎に事前に入力しておく必要があり、この正常時の値
は、実質的にはアラームなどの管理値となり、燃焼振動
を早期に検知することができない。
However, in the conventional monitoring device as described above, it is necessary to input the amplitude value of the combustion oscillation at the normal time for each frequency in advance. In effect, it becomes a management value such as an alarm, and it is not possible to detect combustion vibration at an early stage.

【0007】本発明は、上述の如き問題点を解消するた
めになされたもので、ガスタービンの燃焼器で発生する
燃焼振動を数式モデルにより予測することにより燃焼制
御系調整を容易化すると共に、運転中においては燃焼振
動発生を早期に検知し、設備破損の回避と、設備利用率
の向上、安全性の向上を図ることができる燃焼振動予測
装置を提供することを目的としている。
SUMMARY OF THE INVENTION The present invention has been made to solve the above-described problems, and facilitates adjustment of a combustion control system by predicting a combustion oscillation generated in a combustor of a gas turbine by a mathematical model. It is an object of the present invention to provide a combustion vibration prediction device that can detect the occurrence of combustion vibration at an early stage during operation, avoid damage to equipment, improve equipment utilization, and improve safety.

【0008】[0008]

【課題を解決するための手段】上述の目的を達成するた
めに、請求項1の発明による燃焼振動予測装置は、プラ
ントデータ、気象データから内圧変動を説明する数式モ
デルを構築し、構築した数式モデルに基づいて燃焼振動
の発生し易い領域、発生しにくい領域を求めて出力する
ものである。
In order to achieve the above object, a combustion vibration predicting apparatus according to the first aspect of the present invention constructs a mathematical model for explaining internal pressure fluctuation from plant data and weather data, Based on the model, a region where combustion vibration is likely to occur and a region where combustion vibration is unlikely to occur are obtained and output.

【0009】この構成によれば、プラントデータ、気象
データから構築された内圧変動を説明する数式モデルに
基づいて燃焼振動の発生し易い領域、発生しにくい領域
を求めることが行われ、その結果が出力される。
According to this configuration, a region where combustion oscillation is likely to occur and a region where combustion vibration is unlikely to occur are obtained based on a mathematical model for explaining internal pressure fluctuations constructed from plant data and weather data. Is output.

【0010】また、請求項2の発明による燃焼振動予測
装置は、プラントデータや気象データ、内圧変動の制限
値などを入力する入力手段と、入力されたプラントデー
タ、気象データから燃焼器の内圧変動を数式モデル化す
る内圧変動特性把握手段と、前記内圧変動特性把握手段
により求められた数式モデルに内圧変動の制限値を適用
して燃焼振動の発生し易い領域を求める燃焼振動領域推
定手段と、前記燃焼振動領域推定手段による燃焼振動領
域推定結果を出力する出力手段とを有しているものであ
る。
A combustion vibration prediction apparatus according to a second aspect of the present invention includes an input means for inputting plant data, weather data, a limit value of internal pressure fluctuation, etc., and an internal pressure fluctuation of a combustor based on the input plant data and weather data. Internal pressure fluctuation characteristic grasping means for modeling a mathematical model, and combustion vibration region estimating means for applying a limit value of the internal pressure fluctuation to the mathematical model obtained by the internal pressure fluctuation characteristic grasping means to obtain a region where combustion vibration is likely to occur, Output means for outputting a result of the combustion vibration region estimation by the combustion vibration region estimation means.

【0011】この構成によれば、内圧変動特性把握手段
が入力手段により入力されたプラントデータ、気象デー
タから燃焼器の内圧変動を数式モデル化し、燃焼振動領
域推定手段が数式モデルに内圧変動の制限値を適用して
燃焼振動の発生し易い領域を求め、その燃焼振動領域推
定結果が出力手段から出力される。
According to this configuration, the internal pressure fluctuation characteristic grasping means forms a mathematical model of the internal pressure fluctuation of the combustor from the plant data and weather data inputted by the input means, and the combustion vibration region estimating means restricts the internal pressure fluctuation to the mathematical model. A region where combustion oscillation is likely to occur is obtained by applying the value, and the combustion oscillation region estimation result is output from the output means.

【0012】また、請求項3の発明による燃焼振動予測
装置は、さらに、前記入力手段により入力されたプラン
トデータや気象データなどを時系列に格納するデータベ
ースを有し、前記内圧変動特性把握手段は前記データベ
ースよりデータを取得して燃焼器の内圧変動を数式モデ
ル化するものである。
Further, the combustion vibration predicting apparatus according to the third aspect of the present invention further includes a database for storing plant data, weather data, and the like input by the input means in a time series. Data is obtained from the database to formulate the internal pressure fluctuation of the combustor as a mathematical model.

【0013】この構成によれば、入力手段により入力さ
れたプラントデータや気象データなどがデータベースに
時系列に格納され、内圧変動特性把握手段はそのデータ
ベースよりデータを取得して燃焼器の内圧変動を数式モ
デル化する。
According to this configuration, the plant data, weather data, and the like input by the input means are stored in the database in chronological order, and the internal pressure fluctuation characteristic grasping means obtains the data from the database and detects the internal pressure fluctuation of the combustor. Make a mathematical model.

【0014】また、請求項4の発明による燃焼振動予測
装置は、プラントデータや気象データなどを入力する入
力手段と、入力されたプラントデータや気象データより
燃焼器の内圧変動を推定する内圧変動推定手段と、前記
内圧変動推定手段により推定された内圧変動推定結果を
出力する出力手段とを有しているものである。
According to a fourth aspect of the present invention, there is provided a combustion vibration predicting apparatus for inputting plant data, weather data and the like, and estimating an internal pressure fluctuation of a combustor from the input plant data and weather data. Means, and output means for outputting the internal pressure fluctuation estimation result estimated by the internal pressure fluctuation estimating means.

【0015】この構成によれば、入力手段によって入力
されたプラントデータや気象データより内圧変動推定手
段が燃焼器の内圧変動を推定し、推定された内圧変動推
定結果が出力手段から出力される。
According to this configuration, the internal pressure fluctuation estimating means estimates the internal pressure fluctuation of the combustor from the plant data and the weather data input by the input means, and the estimated internal pressure fluctuation estimation result is output from the output means.

【0016】また、請求項5の発明による燃焼振動予測
装置は、さらに前記入力手段により入力されたプラント
データや気象データなどを時系列に格納するデータベー
スを有し、前記内圧変動推定手段は、データベースに格
納されている最新時刻のデータにより内圧変動の予想値
を推定するものである。
The apparatus for predicting combustion vibration according to the invention of claim 5 further includes a database for storing plant data, weather data, and the like input by the input means in a time series, wherein the internal pressure fluctuation estimating means includes a database. Is used to estimate the expected value of the internal pressure fluctuation based on the latest time data stored in.

【0017】この構成によれば、入力手段により入力さ
れたプラントデータや気象データなどがデータベースに
時系列に格納され、内圧変動推定手段は、そのデータベ
ースに格納されている最新時刻のデータにより内圧変動
の予想値を推定する。
According to this configuration, the plant data, weather data, and the like input by the input means are stored in the database in chronological order, and the internal pressure fluctuation estimating means uses the latest time data stored in the database to determine the internal pressure fluctuation. Estimate the expected value of.

【0018】また、請求項6の発明による燃焼振動予測
装置は、プラントデータ、気象データから内圧変動、N
Ox排出量を説明する数式モデルを構築し、構築した数
式モデルに基づいて燃焼振動の発生しにくい領域、発生
し易い領域を求めて出力するものである。
Further, the combustion vibration predicting device according to the invention of claim 6 is a device for predicting internal pressure fluctuation, N based on plant data and weather data.
A mathematical model describing the amount of Ox emission is constructed, and based on the constructed mathematical model, a region where combustion oscillation is unlikely to occur and a region where combustion vibration is likely to occur are obtained and output.

【0019】この構成によれば、プラントデータ、気象
データから構築された内圧変動、NOx排出量を説明す
る数式モデルに基づいて燃焼振動の発生しにくい領域、
発生し易い領域を求めることがおこなわれ、その結果が
出力される。
According to this configuration, based on a mathematical model describing the internal pressure fluctuation and the NOx emission amount constructed from the plant data and the weather data, an area where combustion oscillation is unlikely to occur,
A region that is likely to occur is determined, and the result is output.

【0020】また、請求項7の発明による燃焼振動予測
装置は、プラントデータや気象データ、内圧変動の制限
値やNOxの規制値などを入力する入力手段と、入力さ
れたプラントデータ、気象データから燃焼器の内圧変動
を数式モデル化する内圧変動特性把握手段と、入力され
たプラントデータ、気象データからNOx排出量を数式
モデル化するNOx排出量特性把握手段と、前記内圧変
動特性把握手段により求められた数式モデルに内圧変動
の制限値を適用し、前記NOx排出量特性把握手段によ
り求められた数式モデルにNOxの規制値を適用して、
NOx排出量が規制値以下でかつ燃焼振動の発生しにく
い領域を求める安全領域推定手段と、前記安全領域推定
手段による安全領域推定結果を出力する出力手段と、を
有しているものである。
Further, the combustion vibration predicting apparatus according to the present invention comprises an input means for inputting plant data, weather data, a limit value of internal pressure fluctuation, a NOx control value, and the like, and an input means for inputting plant data and weather data. The internal pressure fluctuation characteristic grasping means for mathematically modeling the internal pressure fluctuation of the combustor, the NOx emission characteristic grasping means for mathematically modeling the NOx emission from the input plant data and weather data, and the internal pressure fluctuation characteristic grasping means are obtained. Applying the limit value of the internal pressure fluctuation to the obtained mathematical model, and applying the regulation value of NOx to the mathematical model obtained by the NOx emission characteristic grasping means,
The safety region estimating means for obtaining a region where the NOx emission amount is equal to or less than the regulation value and in which combustion oscillation is unlikely to occur, and an output means for outputting a safe region estimation result by the safety region estimating device.

【0021】この構成によれば、内圧変動特性把握手段
が入力手段により入力されたプラントデータ、気象デー
タから燃焼器の内圧変動を数式モデル化し、NOx排出
量特性把握手段が入力手段により入力されたプラントデ
ータ、気象データからNOx排出量を数式モデル化し、
安全領域推定手段が数式モデルに内圧変動の制限値とN
Oxの規制値を適用して燃焼振動の発生しにくい領域を
求め、その安全領域推定結果が出力手段から出力され
る。
According to this configuration, the internal pressure fluctuation characteristic grasping means forms a mathematical model of the internal pressure fluctuation of the combustor from the plant data and the weather data inputted by the input means, and the NOx emission characteristic grasping means is inputted by the input means. Formulate NOx emissions from plant data and weather data into mathematical models,
The safety region estimating means adds the internal pressure fluctuation limit value and N to the mathematical model.
A region where combustion oscillation is unlikely to occur is obtained by applying the regulation value of Ox, and the safety region estimation result is output from the output means.

【0022】また、請求項8の発明による燃焼振動予測
装置は、プラントデータ、気象データから内圧変動、N
Ox、CO排出量を説明する数式モデルを構築し、構築
した数式モデルに基づいて燃焼振動の発生しにくい領
域、発生し易い領域を求めて出力するものである。
Further, the combustion vibration predicting apparatus according to the invention of claim 8 is a device for predicting internal pressure fluctuation, N
A mathematical model for describing the amounts of Ox and CO emissions is constructed, and based on the constructed mathematical model, a region in which combustion oscillation is unlikely to occur and a region in which combustion vibration is likely to occur are obtained and output.

【0023】この構成によれば、プラントデータ、気象
データから構築された内圧変動、NOx排出量、CO排
出量を説明する数式モデルに基づいて燃焼振動の発生し
にくい領域、発生し易い領域を求めることがおこなわ
れ、その結果が出力される。
According to this configuration, a region in which combustion oscillation is unlikely to occur and a region in which combustion vibration is likely to occur are obtained based on a mathematical model for explaining internal pressure fluctuations, NOx emissions, and CO emissions constructed from plant data and weather data. Is performed and the result is output.

【0024】また、請求項9の発明による燃焼振動予測
装置は、プラントデータや気象データ、内圧変動の制限
値やNOx、COの規制値などを入力する入力手段と、
入力されたプラントデータ、気象データから燃焼器の内
圧変動を数式モデル化する内圧変動特性把握手段と、入
力されたプラントデータ、気象データからNOx排出量
を数式モデル化するNOx排出量特性把握手段と、入力
されたプラントデータ、気象データからCO排出量を数
式モデル化するCO排出量特性把握手段と、前記内圧変
動特性把握手段により求められた数式モデルに内圧変動
の制限値を適用し、前記NOx排出量特性把握手段によ
り求められた数式モデルにNOxの規制値を適用し、前
記CO排出量特性把握手段により求められた数式モデル
にCOの規制値を適用して、NOxおよびCO排出量が
規制値以下でかつ燃焼振動の発生しにくい領域を求める
安全領域推定手段と、前記安全領域推定手段による安全
領域推定結果を出力する出力手段と、を有しているもの
である。
Further, the combustion vibration predicting apparatus according to the ninth aspect of the present invention is characterized in that: input means for inputting plant data, weather data, a limit value of internal pressure fluctuation, a NOx, CO regulation value, and the like;
Internal pressure fluctuation characteristic grasping means for mathematically modeling the internal pressure fluctuation of the combustor from the input plant data and weather data, and NOx emission characteristic grasping means for mathematically modeling the NOx emission from the input plant data and weather data A CO emission characteristic grasping means for mathematically modeling the CO emissions from the input plant data and weather data; and a limit value of the internal pressure fluctuation applied to the mathematical model obtained by the internal pressure fluctuation characteristic grasping means. The NOx regulation value is applied to the mathematical expression model obtained by the emission characteristic grasping means, and the CO regulation value is applied to the mathematical expression model obtained by the CO emission characteristic grasping means, whereby NOx and CO emissions are regulated. A safety area estimating means for obtaining an area which is less than or equal to the value and in which combustion oscillation is unlikely to occur, and a safety area estimating result obtained by the safety area estimating means. And output means for, those having a.

【0025】この構成によれば、内圧変動特性把握手段
が入力手段により入力されたプラントデータ、気象デー
タから燃焼器の内圧変動を数式モデル化し、NOx排出
量特性把握手段が入力手段により入力されたプラントデ
ータ、気象データからNOx排出量を数式モデル化し、
CO排出量特性把握手段が入力手段により入力されたプ
ラントデータ、気象データからCO排出量を数式モデル
化し、安全領域推定手段が、数式モデルに内圧変動の制
限値とNOx、COの規制値を適用して、NOxおよび
CO排出量が規制値以下でかつ燃焼振動の発生しにくい
領域を求め、その安全領域推定結果が出力手段から出力
される。
According to this structure, the internal pressure fluctuation characteristic grasping means forms a mathematical model of the internal pressure fluctuation of the combustor from the plant data and weather data inputted by the input means, and the NOx emission characteristic grasping means is inputted by the input means. Formulate NOx emissions from plant data and weather data into mathematical models,
The CO emission characteristic grasping means forms a mathematical model of the CO emission from the plant data and the weather data input by the input means, and the safety region estimating means applies the limit value of the internal pressure fluctuation and the NOx and CO regulation values to the mathematical model. Then, an area where the NOx and CO emissions are less than the regulation value and in which combustion oscillation is less likely to be generated is obtained, and the safety area estimation result is output from the output means.

【0026】また、請求項10の発明による燃焼振動予
測装置は、プラントデータや気象データ、内圧変動の制
限値やNOx、COの規制値などを入力する入力手段
と、入力されたプラントデータ、気象データから数式モ
デル化に用いるデータを選択するフォーカス設定手段
と、選択されたプラントデータ、気象データから燃焼器
の内圧変動を数式モデル化する内圧変動特性把握手段
と、選択されたプラントデータ、気象データからNOx
およびCO排出量を数式モデル化する排出量特性把握手
段と、前記内圧変動特性把握手段により求められた数式
モデルに内圧変動の制限値を適用し、前記排出量特性把
握手段により求められた数式モデルにNOxおよびCO
の規制値を適用して、NOxおよびCO排出量が規制値
以下でかつ燃焼振動の発生しにくい領域を求める安全領
域推定手段と、前記安全領域推定手段による安全領域推
定結果を出力する出力手段と、を有しているものであ
る。
The combustion vibration predicting apparatus according to the tenth aspect of the present invention includes an input means for inputting plant data, weather data, a limit value of internal pressure fluctuation, a NOx, CO control value, and the like; Focus setting means for selecting data to be used for mathematical modeling from data, means for grasping internal pressure fluctuation characteristics for mathematically modeling the internal pressure fluctuation of a combustor from selected plant data and weather data, and selected plant data and weather data From NOx
And an emission characteristic grasping means for mathematically modeling the CO emission, and applying a limit value of the internal pressure fluctuation to the mathematical model obtained by the internal pressure fluctuation characteristic grasping means, and a mathematical expression model obtained by the emission characteristic grasping means. NOx and CO
A safety area estimating means for obtaining an area in which NOx and CO emissions are less than or equal to the regulatory value and in which combustion oscillation is unlikely to occur, and an output means for outputting a safety area estimation result by the safety area estimating means. .

【0027】この構成によれば、内圧変動特性把握手段
が、入力手段により入力され、かつフォーカス設定手段
により選択されたプラントデータ、気象データから燃焼
器の内圧変動を数式モデル化し、排出量特性把握手段
が、フォーカス設定手段により選択されたプラントデー
タ、気象データからNOxおよびCO排出量を数式モデ
ル化し、安全領域推定手段が、数式モデルに内圧変動の
制限値とNOx、COの規制値を適用して、NOxおよ
びCO排出量が規制値以下でかつ燃燃焼振動の発生しに
くい領域を求め、その安全領域推定結果が出力手段から
出力される。
According to this configuration, the internal pressure fluctuation characteristic grasping means forms a mathematical model of the internal pressure fluctuation of the combustor from the plant data and the weather data input by the input means and selected by the focus setting means, and grasps the emission characteristic. Means for mathematically modeling NOx and CO emissions from the plant data and weather data selected by the focus setting means, and the safety region estimating means applying the internal pressure fluctuation limit value and the NOx and CO regulation values to the mathematical model. Thus, an area in which the NOx and CO emissions are less than the regulation value and in which fuel combustion oscillation is unlikely to occur is obtained, and the safety area estimation result is output from the output means.

【0028】また、請求項11の発明による燃焼振動予
測装置は、さらに、前記フォーカス設定手段が、前記入
力手段により指定された領域、または設定モードに基づ
いて、前記入力手段から入力されたプラントデータ、気
象データを選択するものである。
[0028] In the combustion vibration predicting apparatus according to the eleventh aspect, the focus setting means may further include a plant data input from the input means based on an area specified by the input means or a setting mode. , To select weather data.

【0029】この構成によれば、入力手段により指定さ
れた領域、または設定モードに基づいて、フォーカス設
定手段は、入力手段から入力されたプラントデータ、気
象データの選択をおこなう。
According to this configuration, the focus setting means selects the plant data and the weather data input from the input means based on the area specified by the input means or the setting mode.

【0030】また、請求項12の発明による燃焼振動予
測装置は、プラントデータや気象データ、内圧変動の制
限値やNOx、COの規制値などを入力する入力手段
と、入力されたプラントデータ、気象データから数式モ
デル化に用いるデータを選択するフォーカス決定手段
と、選択されたプラントデータ、気象データから燃焼器
の内圧変動を数式モデル化する内圧変動特性把握手段
と、選択されたプラントデータ、気象データからNOx
およびCO排出量を数式モデル化する排出量特性把握手
段と、前記内圧変動特性把握手段により求められた数式
モデルに内圧変動の制限値を適用し、前記排出量特性把
握手段により求められた数式モデルにNOxおよびCO
の規制値を適用して、NOxおよびCO排出量が規制値
以下でかつ燃焼振動の発生しにくい領域を求める安全領
域推定手段と、前記安全領域推定手段による安全領域推
定結果を用いて、次に計測するべきポイントを求める調
整案生成手段と、前記安全領域推定手段による安全領域
推定結果および前記調整案生成手段による計測するべき
ポイントを出力する出力手段と、を有しているものであ
る。
A combustion vibration predicting apparatus according to a twelfth aspect of the present invention provides an input means for inputting plant data, weather data, internal pressure fluctuation limit values, NOx and CO control values, and the like. Focus determination means for selecting data to be used for mathematical modeling from data; internal pressure fluctuation characteristic grasping means for mathematically modeling internal pressure fluctuations of a combustor from selected plant data and weather data; selected plant data and weather data From NOx
And an emission characteristic grasping means for mathematically modeling the CO emission, and applying a limit value of the internal pressure fluctuation to the mathematical model obtained by the internal pressure fluctuation characteristic grasping means, and a mathematical expression model obtained by the emission characteristic grasping means. NOx and CO
The safety region estimating means for obtaining a region in which NOx and CO emissions are equal to or less than the regulation values and in which combustion oscillation is unlikely to occur, and the safety region estimation result obtained by the safety region estimating device, It has adjustment plan generating means for obtaining points to be measured, and output means for outputting the safety area estimation result by the safety area estimating means and the points to be measured by the adjustment plan generating means.

【0031】この構成によれば、内圧変動特性把握手段
が、入力手段により入力され、かつフォーカス決定手段
により選択されたプラントデータ、気象データから燃焼
器の内圧変動を数式モデル化し、排出量特性把握手段
が、フォーカス決定手段により選択されたプラントデー
タ、気象データからNOxおよびCO排出量を数式モデ
ル化し、安全領域推定手段が、数式モデルに内圧変動の
制限値とNOx、COの規制値を適用して、NOxおよ
びCO排出量が規制値以下でかつ燃焼振動の発生しにく
い領域を求め、調整案生成手段が、安全領域推定手段に
よる安全領域推定結果を用いて、次に計測するべきポイ
ントを求め、また、安全領域推定手段による安全領域推
定結果と調整案生成手段による次に計測すべきポイント
が出力手段から出力される。
According to this configuration, the internal pressure fluctuation characteristic grasping means forms a mathematical model of the internal pressure fluctuation of the combustor from the plant data and the weather data input by the input means and selected by the focus determining means, and grasps the emission characteristic. Means for mathematically modeling NOx and CO emissions from the plant data and weather data selected by the focus determining means, and the safety region estimating means applies the internal pressure fluctuation limit value and the NOx and CO regulated values to the mathematical model. Then, an area in which the NOx and CO emissions are less than or equal to the regulation value and in which combustion oscillation is unlikely is obtained, and the adjustment plan generation means obtains a point to be measured next using the safety area estimation result by the safety area estimation means. Also, the safety area estimation result by the safety area estimation means and the next point to be measured by the adjustment plan generation means are output from the output means. It is.

【0032】また、請求項13の発明による燃焼振動予
測装置は、さらに、前記フォーカス決定手段が、前回の
フォーカス決定により選択されたプラントデータ、気象
データに基づいて得られた数式モデルに基づいて、次の
フォーカスを決定するものである。
In the combustion vibration predicting apparatus according to the thirteenth aspect of the present invention, the focus determination means may further include a mathematical model obtained based on plant data and weather data selected by the previous focus determination. The next focus is determined.

【0033】この構成によれば、フォーカス決定手段
は、前回のフォーカス決定により選択されたプラントデ
ータ、気象データに基づいて得られた数式モデルに基づ
いて、次のフォーカスの決定をおこなう。
According to this configuration, the focus determining means determines the next focus based on the mathematical model obtained based on the plant data and weather data selected by the previous focus determination.

【0034】[0034]

【発明の実施の形態】以下に添付の図を参照して、本発
明の燃焼振動予測装置の実施の形態について詳細に説明
する。
BRIEF DESCRIPTION OF THE DRAWINGS FIG. 1 is a perspective view of a combustion vibration predicting apparatus according to an embodiment of the present invention.

【0035】実施の形態1 図1は本発明による燃焼振動予測装置の実施の形態1の
構成を示している。図1において、符号10は燃焼振動
予測装置の全体を示している。燃焼振動予測装置10
は、燃焼器の内圧変動をモデル化する内圧変動特性把握
手段1と、燃焼振動の発生し易い領域を求める燃焼振動
領域推定手段2と、プラントデータや気象データを時系
列に格納するデータベース3と、プラントデータや気象
データ、内圧変動の制限値などを入力する入力手段4
と、燃焼振動領域推定結果を出力する出力手段5とを有
し、入力手段4には燃焼振動予測対象のプラント30が
接続されている。
Embodiment 1 FIG. 1 shows the configuration of Embodiment 1 of a combustion vibration prediction device according to the present invention. In FIG. 1, reference numeral 10 indicates the entire combustion vibration prediction device. Combustion vibration prediction device 10
Are internal pressure fluctuation characteristic grasping means 1 for modeling the internal pressure fluctuation of the combustor, combustion vibration area estimating means 2 for obtaining an area where combustion vibration is likely to occur, and a database 3 for storing plant data and weather data in time series. Input means 4 for inputting plant data, weather data, internal pressure fluctuation limit values, and the like.
And an output unit 5 for outputting a combustion vibration region estimation result. The input vibration unit 4 is connected to a plant 30 whose combustion vibration is to be predicted.

【0036】内圧変動特性把握手段1は、データベース
3に格納されたデータを用いて、内圧変動を説明する数
式モデルを構築する。例えば、燃焼器数をnl 、モデル
化すべき周波数帯数をn2 とすると、次式(1)のよう
な重回帰モデルで内圧変動をモデル化する。 Yij=aij,0+aij,1×X11+aij,2×X12+aij,3×X21+aij,4×X22 …(1)
The internal pressure fluctuation characteristic grasping means 1 uses the data stored in the database 3 to construct a mathematical model for explaining the internal pressure fluctuation. For example, assuming that the number of combustors is n 1 and the number of frequency bands to be modeled is n 2 , the internal pressure fluctuation is modeled by a multiple regression model such as the following equation (1). Y ij = a ij, 0 + a ij, 1 × X 11 + a ij, 2 × X 12 + a ij, 3 × X 21 + a ij, 4 × X 22 (1)

【0037】ここで、 Yij:第i燃焼器の第j周波数帯の内圧変動値(i=1
…n1 、j=1…n2) X11:操作量1の値 X12:操作量2の値 X21:操作できない状態量1の値 X22:操作できない状態量2の値 aij,0、aij,1、aij,2、aij,3、aij,4:係数パラメ
ータ である。
Here, Y ij : the internal pressure fluctuation value of the j-th frequency band of the i-th combustor (i = 1
.. N 1 , j = 1... N 2 ) X 11 : Value of manipulated variable 1 X 12 : Value of manipulated variable 2 X 21 : Value of state variable 1 that cannot be operated X 22 : Value of state variable 2 that cannot be operated a ij, 0 , a ij, 1 , a ij, 2 , a ij, 3 , a ij, 4 are coefficient parameters.

【0038】内圧変動特性把握手段1は、データベース
3に時刻ごとに整理して格納された内圧変動値、操作
量、操作できない状態量を用いて(1)式の係数パラメ
ータa ij,0、aij,1、aij,2、aij,3、aij,4を求め、
これを燃焼振動領域推定手段2に送信する。係数パラメ
ータの解法には、例えば最小二乗法が用いられる。
The internal pressure fluctuation characteristic grasping means 1 includes a database
3. Internal pressure fluctuation values organized and stored by time of day, operation
Using the quantity and the state quantity that cannot be operated
Data a ij, 0, Aij, 1, Aij, 2, Aij, 3, Aij, 4,
This is transmitted to the combustion oscillation region estimation means 2. Coefficient parameter
For example, a least squares method is used for solving the data.

【0039】ここで云う内圧変動値とは、プラント30
に設置された圧力センサ(内圧センサ)31から得られ
たデータをA/D変換し、周波数解析した結果をn2
の周波数帯に区切り、夫々の周波数帯においてある時間
内に得られた最大振幅値である。なお、上記では、説明
の都合上、操作量を2変数、操作できない状態量を2変
数としてモデル式を記述しているが、特に2変数に限る
ものではない。
The internal pressure fluctuation value referred to here is the plant 30
A / D conversion is performed on the data obtained from the pressure sensor (internal pressure sensor) 31 installed in the apparatus, and the result of the frequency analysis is divided into n 2 frequency bands, and the maximum obtained within a certain time in each frequency band is obtained. It is an amplitude value. In the above description, for convenience of explanation, the operation amount is two variables and the state quantity that cannot be operated is two variables. However, the model expression is not limited to two variables.

【0040】燃焼振動領域推定手段2は、内圧変動特性
把握手段1により求められた数式モデルを用いて燃焼振
動の発生し易い領域を求める。
The combustion vibration region estimating means 2 uses the mathematical model obtained by the internal pressure fluctuation characteristic grasping means 1 to find a region where combustion vibration is likely to occur.

【0041】例えば、操作量1、操作量2、操作できな
い状態量1、操作できない状態量2が、それぞれ
X'11、X'12、X'21、X'22の時の第i燃焼器の第j周
波数帯の内圧変動予測値Y'ijは次式(2)で求める。 Y'ij=aij,0+aij,1×X'11+aij,2×X'12+aij,3×X'21+aij,4×X '22 …(2)
[0041] For example, the operation amount 1, the operation amount 2, the state quantity 1 that can not be operated, the state quantity 2 that can not be operated, each X '11, X' 12, X '21, X' i-th combustor when the 22 The internal pressure fluctuation predicted value Y ′ ij of the j-th frequency band is obtained by the following equation (2). Y 'ij = a ij, 0 + a ij, 1 × X' 11 + a ij, 2 × X '12 + a ij, 3 × X' 21 + a ij, 4 × X '22 ... (2)

【0042】ここで、aij,0、aij,1、aij,2
ij,3、aij,4は、内圧変動特性把握手段1から送られ
た係数パラメータである。
Where a ij, 0 , a ij, 1 , a ij, 2 ,
a ij, 3 and a ij, 4 are coefficient parameters sent from the internal pressure fluctuation characteristic grasping means 1.

【0043】第i燃焼器の第j周波数帯の内圧変動に
は、燃焼器や周囲の設備の構造面から制限値が設けられ
ている。入力手段4から送信された第i燃焼器の第j周
波数帯の内圧変動の制限値をZijとすると、 Zij=aij,0+aij,1×X'11+aij,2×X'12+aij,3×X'21+aij,4×X' 22 …(3) となるX'11、X'12、X'21、X'22が存在することにな
る。
The internal pressure fluctuation in the j-th frequency band of the i-th combustor
Limits are set due to the structural aspects of the combustor and surrounding equipment.
ing. The j-th round of the i-th combustor transmitted from the input means 4
The limit value of the internal pressure fluctuation in the wavenumber band is ZijThen, Zij= Aij, 0+ Aij, 1× X '11+ Aij, 2× X '12+ Aij, 3× X 'twenty one+ Aij, 4× X ' twenty two … (3) X '11, X '12, X 'twenty one, X 'twenty twoWill exist
You.

【0044】今、仮に、入力手段4において操作できな
い状態量1及び操作できない状態量2の値が入力され、
これら入力値が燃焼振動領域推定手段2に送信されたと
すると、(3)式のうち、X'11、X'12以外は定数とな
り、(3)式を満たす(X'1 1、X'12)を容易に求める
ことができる。入力手段4から送信されたαk(k=1
…n3)なるゲインにより、 αk×Zij=aij,0+aij,1×X'11+aij,2×X'12+aij,3×X'21+aij, 4 ×X'22 …(4) として(X'11、X'12)を求めれば、各燃焼器の各周波
数帯毎にn3本の線を求めることができる。図2は、こ
れを示したものである。ここで、係数パラメータaij,2
が正であれば、直線の上側が燃焼振動の発生し易い領
域、下側が発生しにくい領域となる。逆に、係数パラメ
ータaij,2が負であれば、直線の下側が燃焼振動の発生
し易い領域、上側が発生しにくい領域となる。
Now, suppose that the value of the state quantity 1 that cannot be operated and the value of the state quantity 2 that cannot be operated are input by the input means 4.
When these input values and sent to the combustion vibration region estimating unit 2, (3) out of the equation, X '11, X' 12 than is a constant, satisfying the expression (3) (X '1 1, X' 12 ) Can be easily obtained. Αk (k = 1) transmitted from the input means 4
... The n 3) becomes the gain, αk × Z ij = a ij , 0 + a ij, 1 × X '11 + a ij, 2 × X' 12 + a ij, 3 × X '21 + a ij, 4 × X' 22 ... If (X ′ 11 , X ′ 12 ) is obtained as (4), n 3 lines can be obtained for each frequency band of each combustor. FIG. 2 illustrates this. Here, the coefficient parameter a ij, 2
Is positive, the upper side of the straight line is a region where combustion vibration is likely to occur, and the lower side is a region where combustion vibration is unlikely to occur. Conversely, if the coefficient parameter a ij, 2 is negative, the lower side of the straight line is a region where combustion vibration is likely to occur, and the upper side is a region where combustion vibration is less likely to occur.

【0045】燃焼振動領域推定手段2は、入力手段4か
ら送信される第i燃焼器の第j周波数帯の内圧変動の制
限値Zij(i=1…n1、j=1…n2)、ゲインαk
(k=1…n3)及び特定の2つを除く変数の値と、内
圧変動特性把握手段1から送信される係数パラメータa
ij,0、aij,1、aij,2、aij,3、aij,4(i=1…
1、j=1…n2)から、全ての燃焼器の全ての周波数
帯について上記の直線を求め、線形計画法の手順に基づ
き、最終的に燃焼振動の発生し易い領域、発生しにくい
領域を求めて出力手段5に送信する。
The combustion oscillation region estimating means 2 transmits the limit value Z ij (i = 1... N 1 , j = 1... N 2 ) of the internal pressure fluctuation in the j-th frequency band of the i-th combustor transmitted from the input means 4. , Gain αk
(K = 1... N 3 ) and the values of the variables except for the specific two, and the coefficient parameter a transmitted from the internal pressure fluctuation characteristic grasping means 1.
ij, 0 , a ij, 1 , a ij, 2 , a ij, 3 , a ij, 4 (i = 1...
n 1 , j = 1... n 2 ), the above-mentioned straight lines are obtained for all the frequency bands of all the combustors. The area is obtained and transmitted to the output means 5.

【0046】データベース3には、内圧変動値、操作
量、操作できない状態量が時刻ごとに整理されて時系列
に格納されており、入力手段4からこれらのデータが送
信されると、このデータベース3内部に追加記憶され
る。
The internal pressure fluctuation value, the operation amount, and the inoperable state amount are stored in the database 3 in a time-series manner arranged for each time, and when these data are transmitted from the input means 4, the database 3 It is additionally stored internally.

【0047】入力手段4は、燃焼振動予測装置10外部
のプラント30から送信されるプラントデータ、気象デ
ータを受信し、データベース3へ送信する。プラントデ
ータ、気象データとは、上述の内圧変動値、操作量、操
作できない状態量を含むものである。また、本入力手段
4に備え付けられたキーボードやタッチスクリーンなど
の装置から、第i燃焼器の第j周波数帯の内圧変動の制
限値Zij(i=1…n 1、j=1…n2)、ゲインαk
(k=1…n3)及び特定の2つを除く変数の値を入力
し、燃焼振動領域推定手段2へ送信する。
The input means 4 is provided outside the combustion vibration prediction device 10.
Plant data and weather data transmitted from the plant 30
The data is received and transmitted to the database 3. Plant de
Data and weather data are the above-mentioned internal pressure fluctuation values, manipulated variables,
Includes state quantities that cannot be produced. This input means
Keyboard and touch screen provided for 4
Control of the internal pressure fluctuation in the j-th frequency band of the i-th combustor
Limit value Zij(I = 1 ... n 1, J = 1 ... nTwo), Gain αk
(K = 1 ... nThree) And the values of the variables except for two specific ones
Then, it is transmitted to the combustion oscillation region estimation means 2.

【0048】出力手段5は燃焼振動領域推定手段2より
送信された推定結果を出力する。図3は横軸をX11、縦
軸をX12として燃焼振動領域を出力した出力例を示して
いる。この例では、ゲインαkごとに等高線のように燃
焼振動領域を表現しており、中央部が燃焼振動の発生し
にくい領域、周辺部ほど発生し易い領域である。出力先
は出力手段に備え付けられたCRTなどの表示装置や印
刷装置などである。
The output means 5 outputs the estimation result transmitted from the combustion oscillation area estimation means 2. FIG. 3 shows an output example in which the combustion oscillation region is output with the horizontal axis being X 11 and the vertical axis being X 12 . In this example, the combustion oscillation region is expressed as a contour line for each gain αk, and the central portion is a region where combustion oscillation is less likely to occur, and the peripheral portion is a region where the combustion oscillation is more likely to occur. The output destination is a display device such as a CRT or a printing device provided in the output means.

【0049】以上に説明したように、燃焼振動予測装置
10によれば、ガスタービンの燃焼器で発生する燃焼振
動を数式モデルにより予測することが行われ、これに基
づき、燃焼制御系調整を容易とすることが可能となり、
設備破損の回避と設備利用率の向上、安全性の向上が可
能となる。
As described above, according to the combustion vibration predicting apparatus 10, the combustion vibration generated in the combustor of the gas turbine is predicted by the mathematical model, and the combustion control system can be easily adjusted based on the prediction. It becomes possible,
It is possible to avoid equipment damage, improve the equipment utilization rate, and improve safety.

【0050】なお、上述の実施の形態では、プラントデ
ータ、気象データをプラント30から入力する構成とし
ているが、これらデータは入力手段4に備え付けられた
キーボードなどから直接手入力しても構わない。また、
モデル構造として、線形の一次式として記述している
が、2次以上の高次モデルとしてもよい。また、プラン
ト30から入力された操作量や操作できない状態量を用
いたモデル式として記述しているが、物理特性などに基
づいて変換した値を用いてもよい。
In the above embodiment, the plant data and the weather data are inputted from the plant 30. However, these data may be directly inputted manually from a keyboard provided in the input means 4. Also,
Although the model structure is described as a linear linear expression, it may be a second-order or higher-order model. In addition, although described as a model formula using the operation amount or the state amount that cannot be operated input from the plant 30, a value converted based on physical characteristics or the like may be used.

【0051】実施の形態2 図4は本発明による燃焼振動予測装置の実施の形態2の
構成を示している。図4において、符号20は燃焼振動
予測装置の全体を示している。燃焼振動予測装置20
は、燃焼器の内圧変動を推定する内圧変動推定手段12
と、プラントデータや気象データを時系列に格納するデ
ータベース13と、プラントデータや気象データなどを
入力する入力手段14と、内圧変動推定結果を出力する
出力手段15とを有し、入力手段14には燃焼振動予測
対象のプラント30が接続されている。
Embodiment 2 FIG. 4 shows the configuration of Embodiment 2 of a combustion vibration prediction device according to the present invention. In FIG. 4, reference numeral 20 indicates the entire combustion vibration prediction device. Combustion vibration prediction device 20
Are internal pressure fluctuation estimating means 12 for estimating the internal pressure fluctuation of the combustor.
And a database 13 for storing plant data and weather data in time series, an input means 14 for inputting plant data and weather data, and the like, and an output means 15 for outputting an internal pressure fluctuation estimation result. Is connected to a plant 30 whose combustion vibration is to be predicted.

【0052】内圧変動推定手段12は、データベース1
3に格納されている最新時刻の内圧変動値、操作量、操
作できない状態量を用いて内圧変動の予想値を推定し、
内圧変動予想値をデータベース13に送信する。例え
ば、燃焼器数をnl 、周波数帯数をn2 とすると、次式
(5)のような重回帰モデルで内圧変動予測値を推定す
る。 Y'ij=aij,0+aij,1×X11+aij,2×X12+aij,3×X21+aij,4×X22 …(5)
The internal pressure fluctuation estimating means 12 stores the database 1
Estimating the expected value of the internal pressure fluctuation using the internal pressure fluctuation value, the operation amount, and the inoperable state amount at the latest time stored in 3;
The internal pressure fluctuation predicted value is transmitted to the database 13. For example, assuming that the number of combustors is n l and the number of frequency bands is n 2 , the internal pressure fluctuation predicted value is estimated using a multiple regression model as shown in the following equation (5). Y ′ ij = a ij, 0 + a ij, 1 × X 11 + a ij, 2 × X 12 + a ij, 3 × X 21 + a ij, 4 × X 22 (5)

【0053】ここで、 Y'ij:第i燃焼器の第j周波数帯の内圧変動予測値
(i=1…n1、j=1…n2) X11:操作量1の値 X12:操作量2の値 X21:操作できない状態量1の値 X22:操作できない状態量2の値 aij,0、aij,1、aij,2、aij,3、aij,4:係数パラメ
ータ である。
Here, Y ′ ij is a predicted value of internal pressure fluctuation in the j-th frequency band of the i-th combustor (i = 1... N 1 , j = 1... N 2 ). X 11 : Value of manipulated variable 1 X 12 : Value of manipulated variable 2 X 21 : Value of state variable 1 that cannot be operated X 22 : Value of state variable 2 that cannot be operated a ij, 0 , a ij, 1 , a ij, 2 , a ij, 3 , a ij, 4 : This is the coefficient parameter.

【0054】係数パラメータは事前解析により求められ
ており、内圧変動推定手段12内部に記憶している。こ
こで云う内圧変動値も、プラント30に設置された圧力
センサ(内圧センサ)31から得られたデータをA/D
変換し、周波数解析した結果をn2 個の周波数帯に区切
り、夫々の周波数帯においてある時間内に得られた最大
振幅値である。
The coefficient parameter is obtained by a preliminary analysis and is stored in the internal pressure fluctuation estimating means 12. The internal pressure fluctuation value referred to here is also obtained by A / D conversion of data obtained from a pressure sensor (internal pressure sensor) 31 installed in the plant 30.
The result of the conversion and frequency analysis is divided into n 2 frequency bands, and is the maximum amplitude value obtained within a certain time in each frequency band.

【0055】なお、上記では説明の都合上、操作量を2
変数、操作できない状態量を2変数としてモデル式を記
述しているが、この場合も、特に2変数に限るものでは
ない。また、モデル構造として線形の一次式として記述
しているが、2次以上の高次モデルやニューラルネット
ワークなどの非線形モデルとしてもよい。また、プラン
ト30から入力された操作量や操作できない状態量を用
いたモデル式として記述しているが、質量収支などの法
則に基づいて変換した値を用いてもよい。
In the above description, for convenience of explanation, the operation amount is set to 2
Although the model formula is described with variables and unoperable state quantities as two variables, this is not particularly limited to two variables. Further, the model structure is described as a linear linear expression, but may be a second-order or higher-order model or a nonlinear model such as a neural network. Further, the description is given as a model formula using the operation amount or the state amount that cannot be operated input from the plant 30, but a value converted based on a law such as a mass balance may be used.

【0056】データベース13には、内圧変動値、操作
量、操作できない状態量と、内圧変動予想値が時刻ごと
に整理されて時系列に格納されており、入力手段14や
内圧変動推定手段12からデータが送信されると、デー
タベース13内部に追加記憶される。
In the database 13, the internal pressure fluctuation value, the operation amount, the inoperable state amount and the estimated internal pressure fluctuation value are arranged in time sequence and stored in a time series, and are stored in the input means 14 and the internal pressure fluctuation estimating means 12. When the data is transmitted, it is additionally stored in the database 13.

【0057】入力手段14は、燃焼振動予測装置20外
部のプラント30から送信されるプラントデータ、気象
データを受信し、データベース13へ送信する。プラン
トデータ、気象データとは、上述の内圧変動値、操作
量、操作できない状態量を含むものである。
The input means 14 receives plant data and weather data transmitted from the plant 30 outside the combustion vibration prediction device 20 and transmits them to the database 13. The plant data and the weather data include the above-mentioned internal pressure fluctuation value, operation amount, and state amount that cannot be operated.

【0058】出力手段15はデータベース13内部に格
納されたデータを出力する。図5は横軸を時刻、縦軸を
第i燃焼器の第j周波数帯の内圧変動実測値Yijとその
予想値Y'ijとして出力した例を示している。出力先は
出力手段に備え付けられたCRTなどの表示装置や印刷
装置などである。
The output means 15 outputs data stored in the database 13. FIG. 5 shows an example in which the horizontal axis represents time and the vertical axis represents the measured internal pressure fluctuation value Y ij and its predicted value Y ′ ij in the j-th frequency band of the i-th combustor. The output destination is a display device such as a CRT or a printing device provided in the output means.

【0059】上述したように、この実施の形態によれ
ば、内圧変動推定値と実測値を同時に出力することがで
き、ガスタービン燃焼器の内圧変動が予定通りのレベル
なのかどうかを判断できるようになり、燃焼振動発生を
早期に検知することが可能で、設備破損の回避と設備利
用率の向上、安全性の向上が可能となる。
As described above, according to this embodiment, the estimated value of the internal pressure fluctuation and the actually measured value can be output simultaneously, and it can be determined whether the internal pressure fluctuation of the gas turbine combustor is at the expected level. Thus, it is possible to detect the occurrence of combustion vibration at an early stage, thereby avoiding equipment damage, improving the equipment utilization rate, and improving safety.

【0060】実施の形態3 図6は本発明による燃焼振動予測装置の実施の形態3の
構成を示している。図6において、符号40は燃焼振動
予測装置の全体を示している。燃焼振動予測装置40
は、燃焼器の内圧変動を数式モデル化する内圧変動特性
把握手段41と、NOx排出量が規制値以下でかつ燃焼
振動の発生しにくい領域を求める安全領域推定手段42
と、プラントデータや気象データを時系列に格納するデ
ータベース43と、プラントデータや気象データ、内圧
変動の制限値やNOxの規制値などを入力する入力手段
44と、安全領域推定結果を出力する出力手段45と、
NOx排出量を数式モデル化するNOx排出量特性把握
手段46とを有し、入力手段44には燃焼振動予測対象
のプラント50が接続されている。なお、図6において
符号49は、各手段を相互に接続する接続部を示す(他
の図においても同じ)。
Third Embodiment FIG. 6 shows the configuration of a third embodiment of the combustion vibration prediction device according to the present invention. In FIG. 6, reference numeral 40 indicates the entire combustion vibration prediction device. Combustion vibration prediction device 40
Are internal pressure fluctuation characteristic grasping means 41 for modeling the internal pressure fluctuation of the combustor as a mathematical model, and a safe area estimating means 42 for finding a region where the NOx emission amount is equal to or less than a regulation value and in which combustion oscillation is hardly generated.
And a database 43 for storing plant data and weather data in time series, input means 44 for inputting plant data and weather data, a limit value of internal pressure fluctuation, a NOx limit value, and the like, and an output for outputting a safety region estimation result. Means 45;
A NOx emission characteristic grasping means 46 for mathematically modeling the NOx emission amount is provided. The input means 44 is connected to a plant 50 whose combustion vibration is to be predicted. In FIG. 6, reference numeral 49 denotes a connecting portion for connecting each means to each other (the same applies to other drawings).

【0061】内圧変動特性把握手段41は、データベー
ス43に格納されたデータを用いて、内圧変動を説明す
る数式モデルを構築する。例えば、燃焼器数をnl 、モ
デル化すべき周波数帯数をn2 とすると、前記(1)式
のような重回帰モデルで内圧変動をモデル化する。
The internal pressure fluctuation characteristic grasping means 41 uses the data stored in the database 43 to construct a mathematical model for explaining the internal pressure fluctuation. For example, assuming that the number of combustors is n l and the number of frequency bands to be modeled is n 2 , the internal pressure fluctuation is modeled by a multiple regression model as in the above equation (1).

【0062】内圧変動特性把握手段41は、データベー
ス43に時刻ごとに整理して格納された内圧変動値、操
作量、操作できない状態量を用いて前記(1)式の係数
パラメータaij,0、aij,1、aij,2、aij,3、aij,4
求め、これを安全領域推定手段42に送信する。係数パ
ラメータの解法には、例えば最小二乗法が用いられる。
The internal pressure fluctuation characteristic grasping means 41 uses the internal pressure fluctuation value, the operation amount, and the inoperable state amount arranged and stored for each time in the database 43, and uses the coefficient parameters a ij, 0 of the above equation (1), a ij, 1 , a ij, 2 , a ij, 3 , and a ij, 4 are obtained and transmitted to the safe region estimation means 42. For solving the coefficient parameter, for example, the least squares method is used.

【0063】ここで云う内圧変動値とは、プラント50
に設置された図示しない圧力センサ(内圧センサ)から
得られたデータをA/D変換し、周波数解析した結果を
2個の周波数帯に区切り、夫々の周波数帯においてあ
る時間内に得られた最大振幅値である。なお、上記で
は、説明の都合上、操作量を2変数、操作できない状態
量を2変数としてモデル式を記述しているが、特に2変
数に限るものではない。
The internal pressure fluctuation value referred to here is the plant 50
A / D conversion is performed on data obtained from a pressure sensor (internal pressure sensor), not shown, which is installed in the system, and the result of frequency analysis is divided into n 2 frequency bands, and each frequency band is obtained within a certain time. This is the maximum amplitude value. In the above description, for convenience of explanation, the operation amount is two variables and the state quantity that cannot be operated is two variables. However, the model expression is not limited to two variables.

【0064】NOx排出量特性把握手段46は、データ
ベース43に格納されたデータを用いて、NOx排出量
を説明する数式モデルを構築する。例えば、つぎの
(6)式であらわされるような重回帰モデルでNOx排
出量をモデル化する。 E=b0+b1×X11+b2×X12+b3×X21+b4×X22 …(6)
The NOx emission characteristic grasping means 46 uses the data stored in the database 43 to construct a mathematical model for explaining the NOx emission. For example, the NOx emission amount is modeled by a multiple regression model represented by the following equation (6). E = b 0 + b 1 × X 11 + b 2 × X 12 + b 3 × X 21 + b 4 × X 22 (6)

【0065】ここで、 E:NOx排出量 X11:操作量1の値 X12:操作量2の値 X21:操作できない状態量1の値 X22:操作できない状態量2の値 b0、b1、b2、b3、b4:係数パラメータ である。Here, E: NOx emission amount X 11 : value of operation amount 1 X 12 : value of operation amount 2 X 21 : value of state amount 1 that cannot be operated X 22 : value of state amount 2 that cannot be operated b 0 , b 1 , b 2 , b 3 , b 4 : coefficient parameters

【0066】NOx排出量特性把握手段46は、データ
ベース43に時刻ごとに整理して格納されたNOx排出
量、操作量、操作できない状態量を用いて前記(6)式
の係数パラメータb0、b1、b2、b3、b4を求め、こ
れを安全領域推定手段42に送信する。係数パラメータ
の解法には、例えば最小二乗法が用いられる。
The NOx emission characteristic grasping means 46 uses the NOx emission amount, the operation amount, and the unoperable state amount arranged and stored for each time in the database 43 to obtain the coefficient parameters b 0 , b of the above equation (6). 1 , b 2 , b 3 , b 4 are obtained and transmitted to the safe area estimation means 42. For solving the coefficient parameter, for example, the least squares method is used.

【0067】なお、上記では、説明の都合上、操作量を
2変数、操作できない状態量を2変数としてモデル式を
記述しているが、特に2変数に限るものではない。
In the above description, for convenience of explanation, the operation amount is two variables, and the state quantity that cannot be operated is two variables. However, the model expression is not limited to the two variables.

【0068】安全領域推定手段42は、内圧変動特性把
握手段41とNOx排出量特性把握手段46により求め
られた数式モデルを用いてNOx排出量が規制値以下で
かつ燃焼振動の発生しにくい領域を求める。
The safety region estimating means 42 uses the mathematical model obtained by the internal pressure fluctuation characteristic grasping means 41 and the NOx emission characteristic grasping means 46 to determine an area where the NOx emission is equal to or less than the regulation value and in which combustion oscillation is unlikely to occur. Ask.

【0069】例えば、操作量1、操作量2、操作できな
い状態量1、操作できない状態量2が、それぞれ
X'11、X'12、X'21、X'22の時の第i燃焼器の第j周
波数帯の内圧変動予測値Y'ijは前記(2)で求める。
このとき、前記(2)式のaij,0、aij,1、aij,2、a
ij,3、aij,4は、内圧変動特性把握手段41から送られ
た係数パラメータである。
[0069] For example, the operation amount 1, the operation amount 2, the state quantity 1 that can not be operated, the state quantity 2 that can not be operated, each X '11, X' 12, X '21, X' i-th combustor when the 22 The internal pressure fluctuation predicted value Y'ij of the j-th frequency band is obtained by the above (2).
At this time, a ij, 0 , a ij, 1 , a ij, 2 , a ij, 0
ij, 3 and a ij, 4 are coefficient parameters sent from the internal pressure fluctuation characteristic grasping means 41.

【0070】第i燃焼器の第j周波数帯の内圧変動に
は、燃焼器や周囲の設備の構造面から制限値が設けられ
ている。入力手段44から送信された第i燃焼器の第j
周波数帯の内圧変動の制限値をZijとすると、前記
(3)式を満たすX'11、X'12、X'21、X'22が存在す
ることになる。
The internal pressure fluctuation in the j-th frequency band of the i-th combustor has a limit value in terms of the structure of the combustor and surrounding equipment. The j-th combustion chamber transmitted from the input means 44
Assuming that the limit value of the internal pressure fluctuation in the frequency band is Z ij , there are X ′ 11 , X ′ 12 , X ′ 21 , and X ′ 22 that satisfy the above expression (3).

【0071】今、仮に、入力手段44において操作でき
ない状態量1および操作できない状態量2の値が入力さ
れ、これら入力値が安全領域推定手段42に送信された
とすると、前記(3)式のうち、X'11、X'12以外は定
数となり、前記(3)式を満たす(X'11、X'12)を容
易に求めることができる。入力手段44から送信された
αk(k=1…n3)なるゲインを用いて、前記(4)
式より(X'11、X'12)を求めれば、各燃焼器の各周波
数帯毎にn3本の線を求めることができる。図7は、こ
れを示したものである。ここで、係数パラメータaij,2
が正であれば、直線の上側が燃焼振動の発生し易い領
域、下側が発生しにくい領域となる。逆に、係数パラメ
ータaij,2が負であれば、直線の下側が燃焼振動の発生
し易い領域、上側が発生しにくい領域となる。
Now, suppose that the values of the inoperable state quantity 1 and the inoperable state quantity 2 are input to the input means 44, and that these input values are transmitted to the safe area estimating means 42. , X ′ 11 , and X ′ 12 are constants, and (X ′ 11 , X ′ 12 ) satisfying the expression (3) can be easily obtained. Using the gain αk (k = 1... N 3 ) transmitted from the input means 44,
If (X ′ 11 , X ′ 12 ) is obtained from the equation, n 3 lines can be obtained for each frequency band of each combustor. FIG. 7 illustrates this. Here, the coefficient parameter a ij, 2
Is positive, the upper side of the straight line is a region where combustion vibration is likely to occur, and the lower side is a region where combustion vibration is unlikely to occur. Conversely, if the coefficient parameter a ij, 2 is negative, the lower side of the straight line is a region where combustion vibration is likely to occur, and the upper side is a region where combustion vibration is less likely to occur.

【0072】安全領域推定手段42は、入力手段44か
ら送信される第i燃焼器の第j周波数帯の内圧変動の制
限値Zij(i=1…n1、j=1…n2)、ゲインαk
(k=1…n3)および特定の2つを除く変数の値と、
内圧変動特性把握手段41から送信される係数パラメー
タaij,0、aij,1、aij,2、aij,3、aij,4(i=1…
1、j=1…n2)から、全ての燃焼器の全ての周波数
帯について上記の直線を求め、線形計画法の手順に基づ
き、最終的に燃焼振動の発生しにくい領域、発生し易い
領域を求める。図8は、これを示したものである。
The safety region estimating means 42 receives the limit value Z ij (i = 1... N 1 , j = 1... N 2 ) of the internal pressure fluctuation in the j-th frequency band of the i-th combustor transmitted from the input means 44. Gain αk
(K = 1... N 3 ) and the values of the variables except for two specific ones,
The coefficient parameters a ij, 0 , a ij, 1 , a ij, 2 , a ij, 3 , a ij, 4 (i = 1...) Transmitted from the internal pressure fluctuation characteristic grasping means 41.
n 1 , j = 1... n 2 ), the above-mentioned straight lines are obtained for all the frequency bands of all the combustors. Find the area. FIG. 8 illustrates this.

【0073】また、安全領域推定手段42は、上述した
処理と同様の処理により、入力手段44から送信される
NOx排出量の規制値F、ゲインβk(k=1…n3
および特定の2つを除く変数の値と、NOx排出量特性
把握手段46から送信される係数パラメータb0、b1
2、b3、b4から、NOxの発生しにくい領域、発生
し易い領域を求める。図9は、これを示したものであ
る。ここで、係数パラメータb2が正であれば、直線の
上側がNOxの発生し易い領域、下側が発生しにくい領
域となる。逆に、係数パラメータb2が負であれば、直
線の下側がNOxの発生し易い領域、上側が発生しにく
い領域となる。
The safe area estimating means 42 performs the same processing as the above-described processing to control the NOx emission amount F transmitted from the input means 44 and the gain βk (k = 1... N 3 ).
And the values of the variables except for the specific two, and the coefficient parameters b 0 , b 1 ,
From b 2 , b 3 , and b 4 , an area where NOx is hardly generated and an area where NOx is easily generated are obtained. FIG. 9 illustrates this. Here, if positive coefficient parameter b 2, upper straight lines is likely to occur region, hardly lower is generated region of NOx. Conversely, if the coefficient parameter b 2 is negative, the lower straight line is likely to occur region, the upper occurs hardly region of NOx.

【0074】そして、安全領域推定手段42は、燃焼振
動の発生しにくい領域、発生し易い領域と、NOxの発
生しにくい領域、発生し易い領域とについて、線形計画
法の手順に基づき、NOx排出量が規制値以下でかつ燃
焼振動の発生しにくい領域を求め、出力手段45に送信
する。図10は、これを示したものである。
Then, the safety region estimating means 42 determines NOx emission for the region where combustion oscillation is unlikely to occur and the region where combustion oscillation is likely to occur, and the region where NOx is unlikely to occur and the region where NOx is likely to occur based on the linear programming procedure. An area in which the amount is equal to or less than the regulation value and in which combustion oscillation is unlikely to occur is obtained and transmitted to the output means 45. FIG. 10 illustrates this.

【0075】データベース43には、内圧変動値、NO
x排出量、操作量、操作できない状態量が時刻ごとに整
理されて時系列に格納されており、入力手段44からこ
れらのデータが送信されると、このデータベース43内
部に追加記憶される。
In the database 43, the internal pressure fluctuation value, NO
The x-discharge amount, the operation amount, and the state amount that cannot be operated are arranged in time series and stored in chronological order. When these data are transmitted from the input means 44, they are additionally stored in the database 43.

【0076】入力手段44は、燃焼振動予測装置40外
部のプラント50から送信されるプラントデータ、気象
データを受信し、データベース43へ送信する。プラン
トデータ、気象データとは、上述の内圧変動値、NOx
排出量、操作量、操作できない状態量を含むものであ
る。また、本入力手段44に備え付けられたキーボード
やタッチスクリーンなどの装置から、第i燃焼器の第j
周波数帯の内圧変動の制限値Zij(i=1…n1、j=
1…n2)、NOx排出量の規制値F、ゲインαk、β
k(k=1…n3)および特定の2つを除く変数の値を
入力し、安全領域推定手段42へ送信する。
The input means 44 receives plant data and weather data transmitted from the plant 50 outside the combustion vibration prediction device 40 and transmits them to the database 43. The plant data and weather data are the above-mentioned internal pressure fluctuation value, NOx
This includes the amount of discharge, the amount of operation, and the amount of state that cannot be operated. In addition, a device such as a keyboard or a touch screen provided to the input means 44 may be used to determine whether the j-th combustor
Limit value Z ij (i = 1... N 1 , j =
1 ... n 2 ), NOx emission regulation value F, gain αk, β
k (k = 1... n 3 ) and the values of the variables except for two specific ones are input and transmitted to the safe area estimation means 42.

【0077】例えば、図11に示す構成のガスタービン
51において、プラントデータとは内圧変動値およびN
Ox排出量に加えて、吸気温度、吸気圧力、吸気流量、
圧縮機出口温度、圧縮機出口圧力、燃料流量、燃料温
度、燃料圧力、排ガス温度、入口案内翼角度、燃焼器バ
イパス弁開度、燃料流量制御弁開度などがある。燃料流
量、燃料圧力、燃料流量制御弁は、燃焼の主となるメイ
ン火炎用と、メイン火炎保炎のためのパイロット火炎用
がある。また、気象データには大気温度、大気圧力、湿
度などがある。前記重回帰モデルに用いる操作量および
操作できない状態量は、これらの中から選定されてい
る。なお、図11において、符号53は圧縮機、符号5
4はタービン、符号55は燃焼器、符号56はメイン燃
料流量制御弁、符号57はパイロット燃料流量制御弁、
符号58は燃焼器バイパス弁、符号59は入口案内翼で
ある。
For example, in the gas turbine 51 having the configuration shown in FIG.
In addition to Ox emissions, intake temperature, intake pressure, intake flow rate,
There are compressor outlet temperature, compressor outlet pressure, fuel flow, fuel temperature, fuel pressure, exhaust gas temperature, inlet guide vane angle, combustor bypass valve opening, fuel flow control valve opening, and the like. The fuel flow rate, fuel pressure, and fuel flow rate control valve include a main flame for main combustion and a pilot flame for main flame holding. The weather data includes atmospheric temperature, atmospheric pressure, humidity, and the like. The manipulated variable and the inoperable state variable used in the multiple regression model are selected from these. In FIG. 11, reference numeral 53 denotes a compressor, and reference numeral 5 denotes a compressor.
4 is a turbine, 55 is a combustor, 56 is a main fuel flow control valve, 57 is a pilot fuel flow control valve,
Reference numeral 58 denotes a combustor bypass valve, and reference numeral 59 denotes an inlet guide vane.

【0078】出力手段45は安全領域推定手段42より
送信された推定結果を出力する。図10は横軸をX11
縦軸をX12として安全領域を出力した出力例を示してい
る。この例では、ゲインαk、βkごとに等高線のよう
に安全領域を表現しており、中央部が燃焼振動の発生し
にくい領域、周辺部ほど発生し易い領域である。出力先
は出力手段45に備え付けられたCRTなどの表示装置
や印刷装置などである。
The output means 45 outputs the estimation result transmitted from the safe area estimation means 42. In FIG. 10, the horizontal axis is X 11 ,
The vertical axis X 12 shows an output example of outputting the safe area. In this example, the safety region is expressed as a contour line for each of the gains αk and βk, and the central portion is a region in which combustion vibration is less likely to occur, and the peripheral portion is a region in which combustion vibration is more likely to occur. The output destination is a display device such as a CRT or a printing device provided in the output unit 45.

【0079】以上に説明したように、燃焼振動予測装置
40によれば、ガスタービンの燃焼器で発生する燃焼振
動とNOx排出量を数式モデルにより予測することがお
こなわれ、これに基づき、燃焼制御系調整を容易とする
ことが可能となり、設備破損の回避と設備利用率の向
上、安全性の向上が可能となる。したがって、燃焼振動
予測装置40を用いることによって、たとえば従来、熟
練した調整員の経験に基づいて調整していたガスタービ
ンの燃焼制御パラメータについて、NOx排出量が規制
値以下で、かつ燃焼振動の発生しにくい安全領域が提示
されるので、たとえばガスタービンの現地調整期間の短
縮が可能であり、また熟練者でなくても現地調整を容易
におこなうことができる。ここで、図12に、ガスター
ビンの燃焼制御系調整に適用可能な出力結果の一例を具
体的に示す。
As described above, according to the combustion vibration predicting device 40, the combustion vibration and the NOx emission generated in the combustor of the gas turbine are predicted by the mathematical model, and the combustion control is performed based on this. System adjustment can be facilitated, and equipment damage can be avoided, the equipment utilization rate can be improved, and safety can be improved. Therefore, by using the combustion vibration prediction device 40, for example, regarding the combustion control parameters of the gas turbine that have been conventionally adjusted based on the experience of a skilled adjuster, the NOx emission is equal to or less than the regulation value and the generation of combustion vibration Since the safety area that is difficult to perform is presented, for example, the on-site adjustment period of the gas turbine can be shortened, and the on-site adjustment can be easily performed even by a non-expert. Here, FIG. 12 specifically shows an example of an output result applicable to the adjustment of the combustion control system of the gas turbine.

【0080】なお、上述の実施の形態では、プラントデ
ータ、気象データをプラント50から入力する構成とし
ているが、これらデータは入力手段44に備え付けられ
たキーボードなどから直接手入力しても構わない。ま
た、モデル構造として、線形の一次式として記述してい
るが、2次以上の高次モデルとしてもよいし、あるいは
操作量や操作できない状態量を乗除した相互項を加えた
モデルとしてもよい。また、プラント50から入力され
た操作量や操作できない状態量を用いたモデル式として
記述しているが、物理特性などに基づいて変換した値を
用いてもよい。
In the above embodiment, the plant data and the weather data are inputted from the plant 50. However, these data may be directly inputted manually from a keyboard provided in the input means 44. In addition, the model structure is described as a linear linear expression, but may be a second-order or higher-order model, or may be a model to which a mutual term obtained by multiplying or dividing the operation amount or the state amount that cannot be operated is added. In addition, although described as a model formula using the operation amount or the state amount that cannot be operated input from the plant 50, a value converted based on physical characteristics or the like may be used.

【0081】実施の形態4図13は本発明による燃焼振
動予測装置の実施の形態4の構成を示している。図13
において、符号60は燃焼振動予測装置の全体を示して
いる。燃焼振動予測装置60は、燃焼器の内圧変動を数
式モデル化する内圧変動特性把握手段41と、NOx排
出量およびCO排出量が規制値以下でかつ燃焼振動の発
生しにくい領域を求める安全領域推定手段62と、プラ
ントデータや気象データを時系列に格納するデータベー
ス63と、プラントデータや気象データ、内圧変動の制
限値、NOxの規制値、COの規制値などを入力する入
力手段64と、安全領域推定結果を出力する出力手段6
5と、NOx排出量を数式モデル化するNOx排出量特
性把握手段46と、CO排出量を数式モデル化するCO
排出量特性把握手段67とを有し、入力手段64には燃
焼振動予測対象のプラント50が接続されている。
Fourth Embodiment FIG. 13 shows the configuration of a fourth embodiment of the combustion vibration prediction device according to the present invention. FIG.
In the figure, reference numeral 60 indicates the entire combustion vibration prediction device. The combustion vibration prediction device 60 includes an internal pressure fluctuation characteristic grasping unit 41 that mathematically models the internal pressure fluctuation of the combustor, and a safe area estimation for obtaining a region where the NOx emission amount and the CO emission amount are equal to or less than the regulation values and in which the combustion oscillation hardly occurs. Means 62, a database 63 for storing plant data and weather data in time series, input means 64 for inputting plant data and weather data, internal pressure fluctuation limit values, NOx control values, CO control values, and the like. Output means 6 for outputting a region estimation result
5, NOx emission characteristic grasping means 46 for mathematically modeling NOx emissions, and CO for mathematically modeling CO emissions.
It has emission characteristic grasping means 67, and the plant 50 of which combustion vibration is to be predicted is connected to the input means 64.

【0082】内圧変動特性把握手段41およびNOx排
出量特性把握手段46は、ともに実施の形態3と同じで
ある。以下には、重複する説明を省略し、実施の形態3
と異なる点についてのみ説明する。
The internal pressure fluctuation characteristic grasping means 41 and the NOx emission characteristic grasping means 46 are the same as those in the third embodiment. In the following, redundant description is omitted, and the third embodiment is omitted.
Only the points different from the above will be described.

【0083】CO排出量特性把握手段67は、データベ
ース63に格納されたデータを用いて、CO排出量を説
明する数式モデルを構築する。例えば、つぎの(7)式
であらわされるような重回帰モデルでCO排出量をモデ
ル化する。 G=c0+c1×X11+c2×X12+c3×X21+c4×X22 …(7)
The CO emission characteristic grasping means 67 uses the data stored in the database 63 to construct a mathematical model for explaining the CO emission. For example, the CO emission amount is modeled by a multiple regression model represented by the following equation (7). G = c 0 + c 1 × X 11 + c 2 × X 12 + c 3 × X 21 + c 4 × X 22 (7)

【0084】ここで、 G:CO排出量 X11:操作量1の値 X12:操作量2の値 X21:操作できない状態量1の値 X22:操作できない状態量2の値 c0、c1、c2、c3、c4:係数パラメータ である。Here, G: CO emission amount X 11 : Value of manipulated variable 1 X 12 : Value of manipulated variable 2 X 21 : Value of inoperable state variable 1 X 22 : Value of inoperable state variable 2 c 0 , c 1 , c 2 , c 3 , c 4 : coefficient parameters

【0085】CO排出量特性把握手段67は、データベ
ース63に時刻ごとに整理して格納されたCO排出量、
操作量、操作できない状態量を用いて前記(7)式の係
数パラメータc0、c1、c2、c3、c4を求め、これを
安全領域推定手段62に送信する。係数パラメータの解
法には、例えば最小二乗法が用いられる。
The CO emission characteristic grasping means 67 stores the CO emissions stored in the database 63 for each time,
The coefficient parameters c 0 , c 1 , c 2 , c 3 , and c 4 of the above equation (7) are obtained using the operation amount and the state amount that cannot be operated, and are transmitted to the safe area estimation means 62. For solving the coefficient parameter, for example, the least squares method is used.

【0086】なお、上記では、説明の都合上、操作量を
2変数、操作できない状態量を2変数としてモデル式を
記述しているが、特に2変数に限るものではない。
In the above description, for convenience of explanation, the operation amount is two variables and the state quantity that cannot be operated is two variables. However, the model expression is not limited to two variables.

【0087】安全領域推定手段62は、内圧変動特性把
握手段41とNOx排出量特性把握手段46とCO排出
量特性把握手段67により求められた数式モデルを用い
てNOx排出量およびCO排出量がともに規制値以下で
かつ燃焼振動の発生しにくい領域を求める。
The safety region estimating means 62 uses the mathematical model obtained by the internal pressure fluctuation characteristic grasping means 41, the NOx emission characteristic grasping means 46, and the CO emission characteristic grasping means 67 to determine both the NOx emission amount and the CO emission amount. A region that is equal to or less than the regulation value and in which combustion vibration is unlikely to occur is determined.

【0088】燃焼振動の発生しにくい領域、発生し易い
領域と、NOxの発生しにくい領域、発生し易い領域の
求め方は、いずれも実施の形態3と同じである。したが
って、以下にはCOの発生しにくい領域、発生し易い領
域の求め方について説明する。
The method of finding the region where combustion vibration is less likely to occur, the region where combustion vibration is more likely to occur, the region where NOx is less likely to be generated and the region where NOx is more likely to occur is the same as in the third embodiment. Therefore, a method of obtaining an area where CO is less likely to be generated and an area where CO is more likely to be generated will be described below.

【0089】安全領域推定手段62は、燃焼振動やNO
xの場合と同様の処理により、入力手段64から送信さ
れるCO排出量の規制値H、ゲインγk(k=1…
3)および特定の2つを除く変数の値と、CO排出量
特性把握手段67から送信される係数パラメータc0
1、c2、c3、c4から、COの発生しにくい領域、発
生し易い領域を求める。
The safety region estimating means 62 determines whether combustion oscillation or NO
By the same processing as in the case of x, the regulated value H of the CO emission amount transmitted from the input means 64 and the gain γk (k = 1.
n 3 ) and the values of the variables except for the specific two, and the coefficient parameter c 0 transmitted from the CO emission characteristic grasping means 67,
From c 1 , c 2 , c 3 , and c 4 , an area where CO is less likely to be generated and an area where CO is more likely to be generated are obtained.

【0090】そして、安全領域推定手段62は、燃焼振
動の発生しにくい領域、発生し易い領域と、NOxの発
生しにくい領域、発生し易い領域と、さらにCOの発生
しにくい領域、発生し易い領域とについて、線形計画法
の手順に基づき、NOx排出量およびCO排出量がとも
に規制値以下でかつ燃焼振動の発生しにくい領域を求
め、出力手段65に送信する。図14は、これを示した
ものである。図14において、係数パラメータc2が正
であれば、直線の上側がCOの発生し易い領域、下側が
発生しにくい領域となる。逆に、係数パラメータc2
負であれば、直線の下側がCO発生し易い領域、上側が
発生しにくい領域となる。
The safety region estimating means 62 calculates a region where combustion oscillation is hardly generated, a region where combustion vibration is easily generated, a region where NOx is hardly generated, a region where NOx is easily generated, a region where CO is hardly generated, and a region where COx is hardly generated. With respect to the region, a region in which both the NOx emission amount and the CO emission amount are equal to or less than the regulation value and in which combustion oscillation is hardly generated is obtained based on the procedure of the linear programming, and transmitted to the output means 65. FIG. 14 illustrates this. 14, if positive coefficient parameter c 2, upper straight lines is likely to occur region, hardly lower is generated region of CO. Conversely, if the coefficient parameter c 2 is negative, the lower side of the straight line is an area where CO is easily generated, and the upper side is an area where CO is hardly generated.

【0091】データベース63には、内圧変動値、NO
x排出量、CO排出量、操作量、操作できない状態量が
時刻ごとに整理されて時系列に格納されており、入力手
段64からこれらのデータが送信されると、このデータ
ベース63内部に追加記憶される。
The database 63 stores the internal pressure fluctuation value, NO
The x-emission amount, the CO-emission amount, the operation amount, and the inoperable state amount are arranged in time series and stored in chronological order, and when these data are transmitted from the input means 64, they are additionally stored in the database 63. Is done.

【0092】入力手段64は、燃焼振動予測装置60外
部のプラント50から送信されるプラントデータ、気象
データを受信し、データベース63へ送信する。また、
本入力手段64に備え付けられたキーボードやタッチス
クリーンなどの装置から、第i燃焼器の第j周波数帯の
内圧変動の制限値Zij(i=1…n1、j=1…n2)、
NOx排出量の規制値F、CO排出量の規制値H、ゲイ
ンαk、βk、γk(k=1…n3)および特定の2つ
を除く変数の値を入力し、安全領域推定手段62へ送信
する。プラントデータ、気象データとは、上述の内圧変
動値、NOx排出量、CO排出量、操作量、操作できな
い状態量を含むものであり、例えば、実施の形態3にお
いて図11に関連して説明した各種データにCO排出量
を加えたものである。
The input means 64 receives plant data and weather data transmitted from the plant 50 outside the combustion vibration prediction device 60 and transmits them to the database 63. Also,
From a device such as a keyboard or a touch screen provided in the input means 64, the limit value Z ij (i = 1... N 1 , j = 1... N 2 ) of the internal pressure fluctuation in the j-th frequency band of the i-th combustor,
The control value F of the NOx emission amount, the control value H of the CO emission amount, the gains αk, βk, γk (k = 1... N 3 ) and the values of the variables except for two specific values are input to the safety region estimation means 62. Send. The plant data and the weather data include the above-described internal pressure fluctuation value, NOx emission amount, CO emission amount, operation amount, and state amount that cannot be operated, and are described, for example, in the third embodiment with reference to FIG. This is the sum of various data and the amount of CO emissions.

【0093】出力手段65は安全領域推定手段62より
送信された推定結果を出力する。図14は横軸をX11
縦軸をX12として安全領域を出力した出力例を示してい
る。この例では、ゲインαk、βk、γkごとに等高線
のように安全領域を表現しており、中央部が燃焼振動の
発生しにくい領域、周辺部ほど発生し易い領域である。
出力先は出力手段65に備え付けられたCRTなどの表
示装置や印刷装置などである。
The output means 65 outputs the estimation result transmitted from the safe area estimation means 62. In FIG. 14, the horizontal axis is X 11 ,
The vertical axis X 12 shows an output example of outputting the safe area. In this example, the safety region is expressed as a contour line for each of the gains αk, βk, and γk, and the central portion is a region in which combustion vibration is less likely to occur, and the peripheral region is more likely to occur.
The output destination is a display device such as a CRT or a printing device provided in the output unit 65.

【0094】以上に説明したように、燃焼振動予測装置
60によれば、ガスタービンの燃焼器で発生する燃焼振
動とNOx排出量とCO排出量を数式モデルにより予測
することがおこなわれ、これに基づき、燃焼制御系調整
を容易とすることが可能となり、設備破損の回避と設備
利用率の向上、安全性の向上が可能となる。したがっ
て、燃焼振動予測装置60を用いることによって、たと
えば従来、熟練した調整員の経験に基づいて調整してい
たガスタービンの燃焼制御パラメータについて、NOx
排出量およびCO排出量が規制値以下で、かつ燃焼振動
の発生しにくい安全領域が提示されるので、たとえばガ
スタービンの現地調整期間の短縮が可能であり、また熟
練者でなくても現地調整を容易におこなうことができ
る。
As described above, according to the combustion vibration predicting device 60, the combustion vibration, NOx emissions, and CO emissions generated in the combustor of the gas turbine are predicted by a mathematical model. Based on this, it becomes possible to easily adjust the combustion control system, thereby avoiding equipment damage, improving the equipment utilization rate, and improving safety. Therefore, by using the combustion vibration predicting device 60, for example, the NOx control of the combustion control parameters of the gas turbine, which has been conventionally adjusted based on the experience of a skilled adjuster,
Since a safety region where the emission amount and the CO emission amount are equal to or less than the regulation values and in which combustion oscillation is hardly generated is presented, for example, it is possible to shorten the on-site adjustment period of the gas turbine, and it is possible to perform on-site adjustment even for a non-expert. Can be easily performed.

【0095】なお、上述の実施の形態では、プラントデ
ータ、気象データをプラント50から入力する構成とし
ているが、これらデータは入力手段64に備え付けられ
たキーボードなどから直接手入力しても構わない。ま
た、モデル構造として、線形の一次式として記述してい
るが、2次以上の高次モデルとしてもよいし、あるいは
操作量や操作できない状態量を乗除した相互項を加えた
モデルとしてもよい。また、プラント50から入力され
た操作量や操作できない状態量を用いたモデル式として
記述しているが、物理特性などに基づいて変換した値を
用いてもよい。
In the above embodiment, the plant data and the weather data are inputted from the plant 50. However, these data may be directly inputted manually from a keyboard provided in the input means 64. In addition, the model structure is described as a linear linear expression, but may be a second-order or higher-order model, or may be a model to which a mutual term obtained by multiplying or dividing the operation amount or the state amount that cannot be operated is added. In addition, although described as a model formula using the operation amount or the state amount that cannot be operated input from the plant 50, a value converted based on physical characteristics or the like may be used.

【0096】実施の形態5 図15は本発明による燃焼振動予測装置の実施の形態5
の構成を示している。図15において、符号70は燃焼
振動予測装置の全体を示している。燃焼振動予測装置7
0は、燃焼器の内圧変動を数式モデル化する内圧変動特
性把握手段71と、NOx排出量およびCO排出量が規
制値以下でかつ燃焼振動の発生しにくい領域を求める安
全領域推定手段72と、プラントデータや気象データを
時系列に格納するデータベース63と、プラントデータ
や気象データ、内圧変動の制限値、NOxの規制値、C
Oの規制値などを入力する入力手段74と、安全領域推
定結果を出力する出力手段65と、NOx排出量および
CO排出量を数式モデル化する排出量特性把握手段76
と、数式モデル化に用いるデータを選択するフォーカス
設定手段78とを有し、入力手段74には燃焼振動予測
対象のプラント50が接続されている。
Fifth Embodiment FIG. 15 shows a combustion vibration predicting apparatus according to a fifth embodiment of the present invention.
Is shown. In FIG. 15, reference numeral 70 indicates the entire combustion vibration prediction device. Combustion vibration prediction device 7
0 is an internal pressure fluctuation characteristic grasping means 71 for modeling the internal pressure fluctuation of the combustor as a mathematical model; a safety region estimating means 72 for finding a region where the NOx emission amount and the CO emission amount are equal to or less than the regulation values and in which combustion oscillation is unlikely to occur; A database 63 for storing plant data and weather data in a time series; plant data and weather data; internal pressure fluctuation limit values; NOx regulation values;
An input unit 74 for inputting a regulation value of O, an output unit 65 for outputting a safety region estimation result, and an emission characteristic grasping unit 76 for numerically modeling the NOx emission amount and the CO emission amount
And a focus setting unit 78 for selecting data used for mathematical modeling. The input unit 74 is connected to the plant 50 whose combustion vibration is to be predicted.

【0097】内圧変動特性把握手段71は、データベー
ス63のしかるべきエリアに記された選択結果に基づい
て、数式モデル化に用いるデータベースの選択をおこな
う機能を有する。内圧変動特性把握手段71のその他の
構成および機能等は実施の形態3の内圧変動特性把握手
段41と同じであるため、重複する説明を省略する。
The internal pressure fluctuation characteristic grasping means 71 has a function of selecting a database to be used for modeling a mathematical expression based on a selection result described in an appropriate area of the database 63. Other configurations, functions, and the like of the internal pressure fluctuation characteristic grasping means 71 are the same as those of the internal pressure fluctuation characteristic grasping means 41 of the third embodiment.

【0098】排出量特性把握手段76は、データベース
63のしかるべきエリアに記された選択結果に基づい
て、数式モデル化に用いるデータベースの選択をおこな
う機能を有する。排出量特性把握手段76のその他の構
成および機能等は、実施の形態3のNOx排出量特性把
握手段46および実施の形態4のCO排出量特性把握手
段67の各構成および機能等を兼ね備えたものと同じで
あるため、重複する説明を省略する。
The emission characteristic grasping means 76 has a function of selecting a database to be used for mathematical modeling based on the selection result described in an appropriate area of the database 63. Other configurations, functions, and the like of the emission characteristic grasping means 76 have the respective configurations and functions of the NOx emission characteristic grasping means 46 of the third embodiment and the CO emission characteristic grasping means 67 of the fourth embodiment. Since it is the same as that of FIG.

【0099】安全領域推定手段72は、内圧変動特性把
握手段71と排出量特性把握手段76により求められた
数式モデルを用いてNOx排出量およびCO排出量がと
もに規制値以下でかつ燃焼振動の発生しにくい領域を求
める。
The safety region estimating means 72 uses the mathematical model obtained by the internal pressure fluctuation characteristic grasping means 71 and the emission characteristic grasping means 76 to determine that both the NOx emission amount and the CO emission amount are equal to or less than the regulation value and that the occurrence of combustion oscillation. Find areas that are difficult to do.

【0100】燃焼振動の発生しにくい領域、発生し易い
領域と、NOxの発生しにくい領域、発生し易い領域
と、COの発生しにくい領域、発生し易い領域のそれぞ
れの求め方は、いずれも実施の形態4と同じであるた
め、重複する説明を省略する。
The method of obtaining the region where combustion oscillation is unlikely to occur, the region where combustion oscillation is likely to occur, the region where NOx is unlikely to occur, the region where it is likely to generate CO, the region where CO is unlikely to occur, and the region where COX is likely to occur is any of: Since it is the same as the fourth embodiment, a duplicate description will be omitted.

【0101】フォーカス設定手段78は、データベース
63に格納されたデータの中から、入力手段74により
入力されたフォーカス設定情報に対応するデータを選択
し、その選択結果をデータベース63のしかるべきメモ
リに記す。ここで、フォーカス設定情報とは、特に限定
しないが、たとえば各変数の上限値および下限値や、選
択範囲とする中心点とその中心点からの最大距離などの
情報である。
The focus setting means 78 selects data corresponding to the focus setting information input by the input means 74 from the data stored in the database 63, and writes the selection result in an appropriate memory of the database 63. . Here, the focus setting information is not particularly limited, but is, for example, information such as an upper limit value and a lower limit value of each variable, a center point to be selected, and a maximum distance from the center point.

【0102】入力手段74は、備えつけのキーボードや
タッチスクリーンなどの装置から、数式モデル化に用い
るデータを選択するためのフォーカス設定情報を入力
し、フォーカス設定手段78に送信する。入力手段74
のその他の構成および機能等は、キーボードやタッチス
クリーンなどの装置から入力された各種データ等の送信
先が安全領域推定手段72(実施の形態4では安全領域
推定手段62)であることを除いて、実施の形態4の入
力手段64と同じであるため、重複する説明を省略す
る。また、プラントデータや気象データについても実施
の形態4と同じである。
The input means 74 inputs focus setting information for selecting data used for mathematical expression modeling from a device such as a built-in keyboard or touch screen, and transmits the information to the focus setting means 78. Input means 74
Other configurations, functions, and the like of, except that the transmission destination of various data and the like input from a device such as a keyboard or a touch screen is the safety region estimation unit 72 (the safety region estimation unit 62 in the fourth embodiment). Since it is the same as the input means 64 of the fourth embodiment, a duplicate description will be omitted. Further, plant data and weather data are the same as in the fourth embodiment.

【0103】出力手段65は、安全領域の推定結果の供
給元が安全領域推定手段72(実施の形態4では安全領
域推定手段62)であることを除いて、実施の形態4と
同じであるため、重複する説明を省略する。図16は横
軸をX11、縦軸をX12として安全領域を出力した出力例
を示している。この例では、ゲインαk、βk、γkご
とに等高線のように安全領域を表現しており、中央部が
燃焼振動の発生しにくい領域、周辺部ほど発生し易い領
域である。
The output means 65 is the same as that of the fourth embodiment except that the supply source of the safety area estimation result is the safety area estimation means 72 (the safety area estimation means 62 in the fourth embodiment). , Overlapping description will be omitted. FIG. 16 shows an output example in which the safety area is output with the horizontal axis being X 11 and the vertical axis being X 12 . In this example, the safety region is expressed as a contour line for each of the gains αk, βk, and γk, and the central portion is a region in which combustion vibration is less likely to occur, and the peripheral region is more likely to occur.

【0104】以上に説明したように、燃焼振動予測装置
70によれば、ガスタービンの燃焼器で発生する燃焼振
動とNOx排出量とCO排出量を数式モデルにより予測
し、燃焼振動の発生しにくい安全領域をマクロ的に広く
求めたり、ミクロ的に高い精度で求めることがおこなわ
れる。したがって、燃焼振動特性の把握が容易となる。
また、燃焼振動特性を把握することによって、燃焼制御
系調整を容易におこなうことが可能となり、設備破損の
回避と設備利用率の向上、安全性の向上が可能となる。
たとえば、燃焼振動予測装置70を用いることによっ
て、従来、熟練した調整員の経験に基づいて調整してい
たガスタービンの燃焼制御パラメータについて、NOx
排出量およびCO排出量が規制値以下で、かつ燃焼振動
の発生しにくい安全領域が提示されるので、たとえばガ
スタービンの現地調整期間の短縮が可能であり、また熟
練者でなくても現地調整を容易におこなうことができ
る。
As described above, according to the combustion vibration predicting apparatus 70, the combustion vibration, NOx emission and CO emission generated in the combustor of the gas turbine are predicted by the mathematical model, and the combustion vibration is hardly generated. It is necessary to obtain a safe area widely on a macro scale or with high precision on a micro scale. Therefore, it is easy to grasp the combustion vibration characteristics.
Further, by grasping the combustion oscillation characteristics, it is possible to easily adjust the combustion control system, and it is possible to avoid equipment damage, improve the equipment utilization rate, and improve safety.
For example, by using the combustion vibration prediction device 70, the NOx control of the combustion control parameters of the gas turbine conventionally adjusted based on the experience of a skilled adjuster
Since a safety region where the emission amount and the CO emission amount are equal to or less than the regulation values and in which combustion oscillation is hardly generated is presented, for example, it is possible to shorten the on-site adjustment period of the gas turbine, and it is possible to perform on-site adjustment even for a non-expert. Can be easily performed.

【0105】なお、上述の実施の形態では、プラントデ
ータ、気象データをプラント50から入力する構成とし
ているが、これらデータは入力手段74に備え付けられ
たキーボードなどから直接手入力しても構わない。ま
た、モデル構造として、線形の一次式として記述してい
るが、2次以上の高次モデルとしてもよいし、あるいは
操作量や操作できない状態量を乗除した相互項を加えた
モデルとしてもよい。また、プラント50から入力され
た操作量や操作できない状態量を用いたモデル式として
記述しているが、物理特性などに基づいて変換した値を
用いてもよい。
In the above-described embodiment, the plant data and the weather data are inputted from the plant 50. However, these data may be directly inputted manually from a keyboard provided in the input means 74. In addition, the model structure is described as a linear linear expression, but may be a second-order or higher-order model, or may be a model to which a mutual term obtained by multiplying or dividing the operation amount or the state amount that cannot be operated is added. In addition, although described as a model formula using the operation amount or the state amount that cannot be operated input from the plant 50, a value converted based on physical characteristics or the like may be used.

【0106】実施の形態6図17は本発明による燃焼振
動予測装置の実施の形態6の構成を示している。図17
において、符号80は燃焼振動予測装置の全体を示して
いる。燃焼振動予測装置80は、燃焼器の内圧変動を数
式モデル化する内圧変動特性把握手段71と、NOx排
出量およびCO排出量が規制値以下でかつ燃焼振動の発
生しにくい領域を求める安全領域推定手段82と、プラ
ントデータや気象データを時系列に格納するデータベー
ス63と、プラントデータや気象データ、内圧変動の制
限値、NOxの規制値、COの規制値などを入力する入
力手段84と、安全領域推定結果および計測すべきポイ
ントを出力する出力手段85と、NOx排出量およびC
O排出量を数式モデル化する排出量特性把握手段76
と、数式モデル化に用いるデータを選択するフォーカス
決定手段88と、安全領域推定結果を用いてつぎに計測
するべきポイントを求める調整案生成手段89とを有
し、入力手段84には燃焼振動予測対象のプラント50
が接続されている。
Sixth Embodiment FIG. 17 shows the configuration of a sixth embodiment of the combustion vibration prediction device according to the present invention. FIG.
In the figure, reference numeral 80 indicates the entire combustion vibration prediction device. The combustion vibration predicting device 80 includes an internal pressure fluctuation characteristic grasping unit 71 that mathematically models the internal pressure fluctuation of the combustor, and a safe area estimation for finding a region where the NOx emission amount and the CO emission amount are equal to or less than the regulation values and in which the combustion oscillation hardly occurs. Means 82, a database 63 for storing plant data and weather data in time series, input means 84 for inputting plant data and weather data, internal pressure fluctuation limit values, NOx control values, CO control values, and the like. Output means 85 for outputting the region estimation result and the points to be measured, NOx emission amount and C
Emission characteristics grasping means 76 for mathematically modeling O emissions
A focus determining means 88 for selecting data used for mathematical modeling, and an adjustment plan generating means 89 for obtaining a point to be measured next using the safety region estimation result. Target plant 50
Is connected.

【0107】内圧変動特性把握手段71および排出量特
性把握手段76は、ともに実施の形態5と同じであるた
め、重複する説明を省略する。
Since the internal pressure fluctuation characteristic grasping means 71 and the discharge amount characteristic grasping means 76 are the same as those in the fifth embodiment, the duplicate description will be omitted.

【0108】安全領域推定手段82は、線形計画法の手
順に基づいて、NOx排出量およびCO排出量がともに
規制値以下でかつ燃焼振動発生レベルの最も小さい予想
最適点を求め、それをフォーカス決定手段88および出
力手段85に送信する。安全領域推定手段82のその他
の構成および機能等は、実施の形態5の安全領域推定手
段72と同じであるため、重複する説明を省略する。
The safety region estimating means 82 obtains a predicted optimum point where both the NOx emission amount and the CO emission amount are equal to or less than the regulation value and the combustion oscillation occurrence level is the smallest, based on the procedure of the linear programming, and determines the focus point. To the means 88 and the output means 85. Other configurations, functions, and the like of the safe area estimating means 82 are the same as those of the safe area estimating means 72 of the fifth embodiment, and thus redundant description will be omitted.

【0109】燃焼振動の発生しにくい領域、発生し易い
領域と、NOxの発生しにくい領域、発生し易い領域
と、COの発生しにくい領域、発生し易い領域のそれぞ
れの求め方は、いずれも実施の形態4と同じであるた
め、重複する説明を省略する。
The method of obtaining the region in which combustion vibration is less likely to occur, the region in which combustion vibration is more likely to occur, the region in which NOx is less likely to be generated, the region in which COx is more likely to occur, the region in which CO is less likely to be generated, and the region in which COx is more likely to be determined, Since it is the same as the fourth embodiment, a duplicate description will be omitted.

【0110】フォーカス決定手段88は、調整初期の段
階においてはデータベース63に格納されたデータの中
から、入力手段84により入力された初期フォーカス設
定情報に対応するデータを選択し、その選択結果をデー
タベース63のしかるべきメモリに記す。ここで、フォ
ーカス設定情報とは、特に限定しないが、たとえば各変
数の上限値および下限値や、フォーカスの中心点座標な
どの情報である。
The focus determining means 88 selects data corresponding to the initial focus setting information input by the input means 84 from the data stored in the database 63 at the initial stage of the adjustment, and stores the result of the selection in the database. Write in 63 appropriate memories. Here, the focus setting information is, for example, but not limited to, information such as an upper limit value and a lower limit value of each variable, and coordinates of a center point of focus.

【0111】フォーカス決定手段88は、初期フォーカ
ス設定情報に対して、安全領域推定手段82により求め
られた予想最適点に基づく変更をおこない、新たなフォ
ーカス設定情報とする。また、フォーカス決定手段88
は、現在のフォーカス設定情報に対して、安全領域推定
手段82により求められた予想最適点に基づく変更をお
こない、新たなフォーカス設定情報とする。そして、フ
ォーカス決定手段88は、データベース63に格納され
たデータの中から、その変更後のフォーカス設定情報に
対応するデータを選択し、その選択結果をデータベース
63のしかるべきメモリに記す。予想最適点が現在のフ
ォーカスの外側に位置する場合には、予想最適点の方向
にフォーカスが移動することになる。また、予想最適点
がフォーカスの内側に位置していても、その位置がフォ
ーカスの周辺部である場合には、移動量を小さくしてフ
ォーカスの位置を設定し直す。図18はフォーカスの移
動の様子を一例として示したものである。この例ではフ
ォーカスの一部が重なりながら移動するようになってい
るが、移動量はこれに限るものではない。
The focus determining means 88 changes the initial focus setting information based on the predicted optimum point obtained by the safe area estimating means 82 to obtain new focus setting information. Further, the focus determining means 88
Makes a change to the current focus setting information based on the predicted optimum point obtained by the safe area estimation means 82, and sets the change as new focus setting information. Then, the focus determining means 88 selects data corresponding to the changed focus setting information from the data stored in the database 63, and records the selection result in an appropriate memory of the database 63. When the predicted optimum point is located outside the current focus, the focus moves in the direction of the predicted optimum point. Further, even if the predicted optimum point is located inside the focus, if the position is located at the periphery of the focus, the movement amount is reduced and the focus position is set again. FIG. 18 shows an example of the movement of the focus. In this example, the focus moves while partly overlapping, but the movement amount is not limited to this.

【0112】調整案生成手段89は、フォーカス決定手
段88により新たに設定されたフォーカス内で、たとえ
ばガスタービン実機のデータ計測がまばらなポイントを
探索し、そのまばらなポイントを追加計測するポイント
として出力手段85に送信する。図18にはこの様子が
示されている。たとえば現在のフォーカスが図18中の
フォーカス1のエリアに設定されていて、既に×印で示
すポイントのデータが計測されており、かつ予想最適点
が図18中の右上方向であると仮定する。この場合、フ
ォーカスが予想最適点方向(図18中、右上方向)に移
動し、フォーカス1よりも図18中の右上方向に位置す
るエリアが新たにフォーカス2として決定される。その
際、調整案生成手段89は、新たにフォーカス2内に含
まれことになった、たとえば図18中に△印で示すポイ
ントについて、データ計測をおこなうことを提案する。
The adjustment plan generating means 89 searches the focus newly set by the focus determining means 88, for example, for a point at which data measurement of the actual gas turbine is sparse, and outputs the sparse point as a point for additional measurement. Transmit to the means 85. FIG. 18 shows this state. For example, it is assumed that the current focus is set in the area of focus 1 in FIG. 18, the data of the point indicated by the mark “x” has already been measured, and the predicted optimum point is the upper right direction in FIG. In this case, the focus moves in the direction of the expected optimum point (upper right direction in FIG. 18), and an area located in the upper right direction in FIG. At this time, the adjustment plan generation means 89 proposes to perform data measurement on a point newly included in the focus 2, for example, a point indicated by a triangle in FIG. 18.

【0113】入力手段84は、備えつけのキーボードや
タッチスクリーンなどの装置から、数式モデル化に用い
るデータを選択するための初期フォーカス設定情報を入
力し、フォーカス決定手段88に送信する。入力手段8
4のその他の構成および機能等は、キーボードやタッチ
スクリーンなどの装置から入力された各種データ等の送
信先が安全領域推定手段82(実施の形態4では安全領
域推定手段62)であることを除いて、実施の形態4の
入力手段64と同じであるため、重複する説明を省略す
る。また、プラントデータや気象データについても実施
の形態4と同じである。
The input means 84 inputs initial focus setting information for selecting data to be used for mathematical modeling from a device such as a built-in keyboard or touch screen, and transmits the information to the focus determining means 88. Input means 8
Other configurations, functions, and the like of the fourth embodiment except that the transmission destination of various data and the like input from a device such as a keyboard and a touch screen is the safety region estimation unit 82 (the safety region estimation unit 62 in the fourth embodiment). Since the input unit 64 is the same as the input unit 64 according to the fourth embodiment, a duplicate description will be omitted. Further, plant data and weather data are the same as in the fourth embodiment.

【0114】出力手段85は、安全領域推定手段82よ
り送信された推定結果、フォーカス決定手段88により
設定されたフォーカスエリア、調整案生成手段89より
送信された計測ポイントを出力する。出力先は出力手段
85に備え付けられたCRTなどの表示装置や印刷装置
などである。図18は横軸をX11、縦軸をX12として安
全領域を出力した出力例を示している。この例では、ゲ
インαk、βk、γkごとに等高線のように安全領域を
表現しており、中央部が燃焼振動の発生しにくい領域、
周辺部ほど発生し易い領域である。図18において、フ
ォーカス1、フォーカス2およびフォーカス3で示す矩
形領域は、フォーカスが順次移動していく様子を示した
ものである。また、×印はフォーカス1でのデータ計測
済みのポイント、△印および○印はそれぞれフォーカス
2およびフォーカス3でのデータ追加計測ポイント案を
示す。
The output unit 85 outputs the estimation result transmitted from the safe area estimation unit 82, the focus area set by the focus determination unit 88, and the measurement points transmitted from the adjustment plan generation unit 89. The output destination is a display device such as a CRT or a printing device provided in the output unit 85. FIG. 18 shows an output example in which the horizontal axis is X 11 and the vertical axis is X 12 to output the safety region. In this example, the safety region is expressed as a contour line for each of the gains αk, βk, and γk, and the central portion is a region where combustion vibration is unlikely to occur,
This is an area that is more likely to occur in the peripheral area. In FIG. 18, rectangular areas indicated by focus 1, focus 2 and focus 3 show how the focus moves sequentially. In addition, crosses indicate points at which data has been measured at focus 1 and triangles and circles indicate data addition measurement points at focus 2 and focus 3, respectively.

【0115】以上に説明したように、燃焼振動予測装置
80によれば、ガスタービンの燃焼器で発生する燃焼振
動とNOx排出量とCO排出量を数式モデルにより予測
し、燃焼振動の発生しにくい安全領域を高い精度で求め
るとともに、より安全な運転ポイントを探索するために
新たな計測ポイントを提示することがおこなわれる。し
たがって、NOx排出量とCO排出量の規制値以下で、
かつ燃焼振動が最も発生しにくいポイント、すなわち最
適運転ポイントを求めることが可能となる。これによっ
て、燃焼制御系の調整が容易となり、従来よりも短い期
間で調整することが可能となる。また、設備破損の回避
と設備利用率の向上、安全性の向上が可能となる。たと
えば、燃焼振動予測装置80を用いることによって、従
来、熟練した調整員の経験に基づいて調整していたガス
タービンの燃焼制御パラメータについて、NOx排出量
およびCO排出量が規制値以下で、かつ燃焼振動が最も
発生しにくい運転ポイントが提示されるので、たとえば
ガスタービンの現地調整期間の短縮が可能であり、また
熟練者でなくても現地調整を容易におこなうことができ
る。
As described above, according to the combustion vibration predicting device 80, the combustion vibration, NOx emissions, and CO emissions generated in the combustor of the gas turbine are predicted by the mathematical model, and the combustion vibration is hardly generated. A safe area is obtained with high accuracy, and a new measurement point is presented to search for a safer driving point. Therefore, when the NOx emissions and CO emissions are below the regulation values,
In addition, it is possible to obtain a point where combustion vibration is least likely to occur, that is, an optimum operation point. This facilitates adjustment of the combustion control system, and allows adjustment in a shorter period than before. In addition, it is possible to avoid equipment damage, improve the equipment utilization rate, and improve safety. For example, by using the combustion vibration prediction device 80, regarding the combustion control parameters of the gas turbine, which have been conventionally adjusted based on the experience of a skilled adjuster, the NOx emission and the CO emission are equal to or less than the regulation values and the combustion is controlled. Since the operating point where vibration is least likely to occur is presented, for example, the on-site adjustment period of the gas turbine can be shortened, and on-site adjustment can be easily performed even by a non-expert.

【0116】なお、上述の実施の形態では、プラントデ
ータ、気象データをプラント50から入力する構成とし
ているが、これらデータは入力手段84に備え付けられ
たキーボードなどから直接手入力しても構わない。ま
た、モデル構造として、線形の一次式として記述してい
るが、2次以上の高次モデルとしてもよいし、あるいは
操作量や操作できない状態量を乗除した相互項を加えた
モデルとしてもよい。また、プラント50から入力され
た操作量や操作できない状態量を用いたモデル式として
記述しているが、物理特性などに基づいて変換した値を
用いてもよい。
In the above embodiment, the plant data and the weather data are inputted from the plant 50. However, these data may be directly inputted manually from a keyboard provided in the input means 84. In addition, the model structure is described as a linear linear expression, but may be a second-order or higher-order model, or may be a model to which a mutual term obtained by multiplying or dividing the operation amount or the state amount that cannot be operated is added. In addition, although described as a model formula using the operation amount or the state amount that cannot be operated input from the plant 50, a value converted based on physical characteristics or the like may be used.

【0117】[0117]

【発明の効果】以上の説明から理解される如く、本発明
による燃焼振動予測装置(請求項1)によれば、プラン
トデータ、気象データから構築された内圧変動を説明す
る数式モデルに基づいて燃焼振動の発生し易い領域、発
生しにくい領域を求めることが行われ、その結果が出力
されるから、燃焼制御系調整を容易とすることが可能と
なり、設備破損の回避と設備利用率の向上、安全性の向
上が可能となる。
As will be understood from the above description, according to the combustion vibration predicting apparatus (claim 1) of the present invention, combustion is performed based on a mathematical model describing internal pressure fluctuations constructed from plant data and weather data. An area where vibration is likely to occur and an area where vibration is unlikely to occur are obtained, and the results are output, so that it becomes possible to easily adjust the combustion control system, avoid equipment damage and improve equipment utilization, Safety can be improved.

【0118】また、本発明による燃焼振動予測装置(請
求項2)によれば、内圧変動特性把握手段が入力手段よ
り入力されたプラントデータ、気象データから燃焼器の
内圧変動を数式モデル化し、燃焼振動領域推定手段が数
式モデルに内圧変動の制限値を適用して燃焼振動の発生
し易い領域を求め、その燃焼振動領域推定結果が出力手
段に出力されるから、燃焼制御系調整を容易とすること
が可能となり、設備破損の回避と設備利用率の向上、安
全性の向上が可能となる。
Further, according to the combustion vibration predicting apparatus of the present invention, the internal pressure fluctuation characteristic grasping means forms a mathematical model of the internal pressure fluctuation of the combustor from the plant data and weather data inputted from the input means, and performs combustion. The vibration region estimating means applies the limit value of the internal pressure fluctuation to the mathematical model to obtain a region where combustion vibration is likely to occur, and the combustion vibration region estimation result is output to the output means, thereby facilitating adjustment of the combustion control system. This makes it possible to avoid equipment damage, increase the equipment utilization rate, and improve safety.

【0119】また、本発明による燃焼振動予測装置(請
求項3)によれば、入力手段により入力されたプラント
データや気象データなどがデータベースに時系列に格納
され、内圧変動特性把握手段はそのデータベースよりデ
ータを取得して燃焼器の内圧変動を数式モデル化するか
ら、より的確に、燃焼振動の発生し易い領域を求めるこ
とが行われる。
Further, according to the combustion vibration predicting apparatus of the present invention (claim 3), the plant data, weather data, and the like input by the input means are stored in a database in chronological order, and the internal pressure fluctuation characteristic grasping means is stored in the database. Since more data is acquired and the fluctuation of the internal pressure of the combustor is modeled as a mathematical model, a region where combustion oscillation is likely to occur is more accurately determined.

【0120】また、本発明による燃焼振動予測装置(請
求項4)によれば、入力手段によって入力されたプラン
トデータや気象データより内圧変動推定手段が燃焼器の
内圧変動を推定し、推定された内圧変動推定結果が出力
手段に出力されるから、ガスタービン燃焼器の内圧変動
が予定通りのレベルなのかどうかを判断できるようにな
り、燃焼振動発生を早期に検知することが可能で、設備
破損の回避と設備利用率の向上、安全性の向上が可能と
なる。
Further, according to the combustion vibration predicting device of the present invention (claim 4), the internal pressure fluctuation estimating means estimates the internal pressure fluctuation of the combustor from the plant data and the weather data input by the input means, and estimates the fluctuation. Since the internal pressure fluctuation estimation result is output to the output means, it is possible to determine whether the internal pressure fluctuation of the gas turbine combustor is at the expected level, it is possible to detect the occurrence of combustion vibration at an early stage, and the equipment is damaged. It is possible to avoid problems, improve the capacity factor, and improve safety.

【0121】また、本発明による、燃焼振動予測装置
(請求項5)によれば、入力手段により入力されたプラ
ントデータや気象データなどがデータベースに時系列に
格納され、内圧変動推定手段は、そのデータベースに格
納されている最新時刻のデータにより内圧変動の予想値
を推定するから、内圧変動の予想が的確に行われる。
Further, according to the combustion vibration predicting device (claim 5) of the present invention, the plant data, weather data, and the like input by the input means are stored in a database in chronological order. Since the expected value of the internal pressure fluctuation is estimated from the latest time data stored in the database, the prediction of the internal pressure fluctuation is accurately performed.

【0122】また、本発明による燃焼振動予測装置(請
求項6)によれば、プラントデータ、気象データから構
築された内圧変動、NOx排出量を説明する数式モデル
に基づいて燃焼振動の発生しにくい領域、発生し易い領
域を求めることがおこなわれ、その結果が出力されるか
ら、燃焼制御系調整を容易とすることが可能となり、設
備破損の回避と設備利用率の向上、安全性の向上が可能
となる。
Further, according to the combustion vibration predicting apparatus of the present invention (claim 6), it is difficult to generate combustion vibration based on a mathematical model describing the internal pressure fluctuation and NOx emission amount constructed from plant data and weather data. The region and the region that is likely to occur are determined, and the result is output, so that it is possible to easily adjust the combustion control system, and to avoid equipment damage, improve equipment utilization, and improve safety. It becomes possible.

【0123】また、本発明による燃焼振動予測装置(請
求項7)によれば、内圧変動特性把握手段が入力手段に
より入力されたプラントデータ、気象データから燃焼器
の内圧変動を数式モデル化し、NOx排出量特性把握手
段が入力手段により入力されたプラントデータ、気象デ
ータからNOx排出量を数式モデル化し、安全領域推定
手段が数式モデルに内圧変動の制限値とNOxの規制値
を適用して、NOx排出量が規制値以下でかつ燃焼振動
の発生しにくい領域を求め、その安全領域推定結果が出
力手段から出力されるから、燃焼制御系調整を容易とす
ることが可能となり、設備破損の回避と設備利用率の向
上、安全性の向上が可能となる。
Further, according to the combustion vibration predicting apparatus of the present invention, the internal pressure fluctuation characteristic grasping means forms a mathematical model of the internal pressure fluctuation of the combustor from the plant data and the weather data inputted by the input means, and NOx The emission characteristic grasping means forms a mathematical model of the NOx emission from the plant data and the weather data inputted by the input means, and the safety region estimating means applies the internal pressure fluctuation limit value and the NOx regulation value to the mathematical model to obtain NOx. An area in which the amount of emission is less than the regulation value and in which combustion oscillation is unlikely to occur is obtained, and the safety area estimation result is output from the output means, so that it becomes possible to easily adjust the combustion control system, thereby avoiding equipment damage. It is possible to improve equipment utilization and safety.

【0124】また、本発明による燃焼振動予測装置(請
求項8)によれば、プラントデータ、気象データから構
築された内圧変動、NOx排出量、CO排出量を説明す
る数式モデルに基づいて燃焼振動の発生しにくい領域、
発生し易い領域を求めることがおこなわれ、その結果が
出力されるから、燃焼制御系調整を容易とすることが可
能となり、設備破損の回避と設備利用率の向上、安全性
の向上が可能となる。
Further, according to the combustion vibration predicting apparatus of the present invention (claim 8), the combustion vibration is estimated based on a mathematical model describing the internal pressure fluctuation, NOx emission and CO emission constructed from plant data and weather data. Area where the occurrence of
Since the area where the occurrence is likely to occur is obtained and the result is output, it is possible to easily adjust the combustion control system, and it is possible to avoid equipment damage, improve the equipment utilization rate, and improve safety. Become.

【0125】また、本発明による燃焼振動予測装置(請
求項9)によれば、内圧変動特性把握手段が入力手段に
より入力されたプラントデータ、気象データから燃焼器
の内圧変動を数式モデル化し、NOx排出量特性把握手
段が入力手段により入力されたプラントデータ、気象デ
ータからNOx排出量を数式モデル化し、CO排出量特
性把握手段が入力手段により入力されたプラントデー
タ、気象データからCO排出量を数式モデル化し、安全
領域推定手段が、数式モデルに内圧変動の制限値とNO
x、COの規制値を適用して、NOxおよびCO排出量
が規制値以下でかつ燃焼振動の発生しにくい領域を求
め、その安全領域推定結果が出力手段から出力されるか
ら、燃焼制御系調整を容易とすることが可能となり、設
備破損の回避と設備利用率の向上、安全性の向上が可能
となる。
Further, according to the combustion vibration predicting apparatus of the present invention, the internal pressure fluctuation characteristic grasping means forms a mathematical model of the internal pressure fluctuation of the combustor from the plant data and the weather data inputted by the input means, and NOx The emission characteristic grasping means formulates the NOx emissions from the plant data and the weather data input by the input means, and the CO emission characteristic grasping means formulates the CO emissions from the plant data and the weather data input by the input means. The safety area estimating means calculates the limit value of the internal pressure fluctuation and the NO
By applying the regulated values of x and CO, a region where the NOx and CO emissions are less than the regulated value and in which combustion oscillation is unlikely to be generated is obtained, and the safety region estimation result is output from the output means. It is possible to avoid equipment damage, improve the equipment utilization rate, and improve safety.

【0126】また、本発明による燃焼振動予測装置(請
求項10)によれば、内圧変動特性把握手段が、入力手
段により入力され、かつフォーカス設定手段により選択
されたプラントデータ、気象データから燃焼器の内圧変
動を数式モデル化し、排出量特性把握手段が、フォーカ
ス設定手段により選択されたプラントデータ、気象デー
タからNOxおよびCO排出量を数式モデル化し、安全
領域推定手段が、数式モデルに内圧変動の制限値とNO
x、COの規制値を適用して、NOxおよびCO排出量
が規制値以下でかつ燃燃焼振動の発生しにくい領域を求
め、その安全領域推定結果が出力手段から出力されるか
ら、燃焼制御系調整を容易とすることが可能となり、設
備破損の回避と設備利用率の向上、安全性の向上が可能
となる。
Further, according to the combustion vibration predicting device of the present invention, the internal pressure fluctuation characteristic grasping means is provided with a combustor based on plant data and weather data inputted by the input means and selected by the focus setting means. Of the internal pressure fluctuation of the internal pressure fluctuation, the emission characteristic grasping means mathematically models the NOx and CO emissions from the plant data and the weather data selected by the focus setting means, and the safety region estimating means converts the internal pressure fluctuation into a mathematical model. Limit value and NO
Since the control values of x and CO are applied to obtain a region where the NOx and CO emissions are less than the control values and in which the occurrence of fuel-combustion oscillation is difficult, and the safety region estimation result is output from the output means, the combustion control system Adjustment can be facilitated, and equipment damage can be avoided, the equipment utilization rate can be improved, and safety can be improved.

【0127】また、本発明による燃焼振動予測装置(請
求項11)によれば、入力手段により指定された領域、
または設定モードに基づいて、フォーカス設定手段は、
入力手段から入力されたプラントデータ、気象データの
選択をおこなうため、燃焼振動の発生しにくい安全領域
をマクロ的に広く求めたり、ミクロ的に高い精度で求め
ることが可能となる。
Further, according to the combustion vibration predicting apparatus of the present invention (claim 11), the area specified by the input means,
Or, based on the setting mode, the focus setting means
Since the plant data and the weather data input from the input means are selected, it is possible to macro-widely obtain a safe area where combustion vibration is unlikely to occur, or to obtain microscopically high-precision.

【0128】また、本発明による燃焼振動予測装置(請
求項12)によれば、内圧変動特性把握手段が、入力手
段により入力され、かつフォーカス決定手段により選択
されたプラントデータ、気象データから燃焼器の内圧変
動を数式モデル化し、排出量特性把握手段が、フォーカ
ス決定手段により選択されたプラントデータ、気象デー
タからNOxおよびCO排出量を数式モデル化し、安全
領域推定手段が、数式モデルに内圧変動の制限値とNO
x、COの規制値を適用して、NOxおよびCO排出量
が規制値以下でかつ燃焼振動の発生しにくい領域を求
め、調整案生成手段が、安全領域推定手段による安全領
域推定結果を用いて、次に計測するべきポイントを求
め、また、安全領域推定手段による安全領域推定結果及
び調整案生成手段による計測すべきポイントが出力手段
から出力されるから、燃焼制御系調整を容易とすること
が可能となり、設備破損の回避と設備利用率の向上、安
全性の向上が可能となる。
Further, according to the combustion vibration predicting device of the present invention, the internal pressure fluctuation characteristic grasping means is provided with a combustor based on plant data and weather data inputted by the input means and selected by the focus determining means. Formula of the internal pressure fluctuation of the NOx and CO emissions from the plant data and the weather data selected by the focus determining means, and the safety region estimating means converts the internal pressure fluctuation into a mathematical model. Limit value and NO
Applying the regulated values of x and CO to obtain a region where NOx and CO emissions are less than or equal to the regulated values and in which combustion oscillation is unlikely to occur, and the adjustment plan generating means uses the safety area estimation result by the safety area estimating means. Then, the point to be measured next is obtained, and the safety area estimation result by the safety area estimation means and the point to be measured by the adjustment plan generation means are output from the output means, so that the combustion control system can be easily adjusted. As a result, it is possible to avoid equipment damage, improve the equipment utilization rate, and improve safety.

【0129】また、本発明による燃焼振動予測装置(請
求項13)によれば、フォーカス決定手段は、前回のフ
ォーカス決定により選択されたプラントデータ、気象デ
ータに基づいて得られた数式モデルに基づいて、次のフ
ォーカスの決定をおこなうため、最適運転ポイントの探
索が容易となる。
According to the combustion vibration predicting apparatus of the present invention (claim 13), the focus determination means is based on a mathematical model obtained based on the plant data and weather data selected by the previous focus determination. Since the next focus is determined, it is easy to search for the optimum driving point.

【図面の簡単な説明】[Brief description of the drawings]

【図1】本発明の実施の形態1に関する燃焼振動予測装
置の構成図である。
FIG. 1 is a configuration diagram of a combustion vibration prediction device according to a first embodiment of the present invention.

【図2】本発明の実施の形態1に関する燃焼振動予測装
置における燃焼振動領域の推定法に関する原理図であ
る。
FIG. 2 is a principle diagram regarding a method of estimating a combustion oscillation region in the combustion oscillation prediction device according to the first embodiment of the present invention.

【図3】本発明の実施の形態1に関する燃焼振動予測装
置における燃焼振動領域の出力例を示すグラフである。
FIG. 3 is a graph showing an output example of a combustion vibration region in the combustion vibration prediction device according to the first embodiment of the present invention.

【図4】本発明の実施の形態2に関する燃焼振動予測装
置の構成図である。
FIG. 4 is a configuration diagram of a combustion vibration prediction device according to a second embodiment of the present invention.

【図5】本発明の実施の形態2に関する燃焼振動予測装
置における出力例を示すグラフである。
FIG. 5 is a graph showing an output example in the combustion vibration prediction device according to the second embodiment of the present invention.

【図6】本発明の実施の形態3に関する燃焼振動予測装
置の構成図である。
FIG. 6 is a configuration diagram of a combustion vibration prediction device according to a third embodiment of the present invention.

【図7】本発明の実施の形態3に関する燃焼振動予測装
置における安全領域の推定法に関する原理図である。
FIG. 7 is a principle diagram relating to a method of estimating a safety region in a combustion vibration prediction device according to Embodiment 3 of the present invention.

【図8】本発明の実施の形態3に関する燃焼振動予測装
置における安全領域の推定法に関する原理図である。
FIG. 8 is a principle diagram relating to a method for estimating a safety region in a combustion vibration prediction device according to Embodiment 3 of the present invention.

【図9】本発明の実施の形態3に関する燃焼振動予測装
置におけるNOx規制レベルの推定法に関する原理図で
ある。
FIG. 9 is a principle diagram regarding a method of estimating a NOx regulation level in a combustion vibration prediction device according to a third embodiment of the present invention.

【図10】本発明の実施の形態3に関する燃焼振動予測
装置における安全領域の出力例を示すグラフである。
FIG. 10 is a graph showing an output example of a safety region in the combustion vibration prediction device according to the third embodiment of the present invention.

【図11】本発明を適用可能なガスタービンの系統構成
の一例を模式的に示す構成図である。
FIG. 11 is a configuration diagram schematically illustrating an example of a system configuration of a gas turbine to which the present invention can be applied.

【図12】本発明の実施の形態3に関する燃焼振動予測
装置の具体的な適用例を示すグラフである。
FIG. 12 is a graph showing a specific application example of the combustion vibration prediction device according to the third embodiment of the present invention.

【図13】本発明の実施の形態4に関する燃焼振動予測
装置の構成図である。
FIG. 13 is a configuration diagram of a combustion vibration prediction device according to a fourth embodiment of the present invention.

【図14】本発明の実施の形態4に関する燃焼振動予測
装置における安全領域の出力例を示すグラフである。
FIG. 14 is a graph showing an output example of a safety region in the combustion vibration prediction device according to the fourth embodiment of the present invention.

【図15】本発明の実施の形態5に関する燃焼振動予測
装置の構成図である。
FIG. 15 is a configuration diagram of a combustion vibration prediction device according to a fifth embodiment of the present invention.

【図16】本発明の実施の形態5に関する燃焼振動予測
装置における安全領域の出力例を示すグラフである。
FIG. 16 is a graph showing an output example of a safety region in the combustion vibration prediction device according to the fifth embodiment of the present invention.

【図17】本発明の実施の形態6に関する燃焼振動予測
装置の構成図である。
FIG. 17 is a configuration diagram of a combustion vibration prediction device according to a sixth embodiment of the present invention.

【図18】本発明の実施の形態6に関する燃焼振動予測
装置における安全領域の出力例を示すグラフである。
FIG. 18 is a graph showing an output example of a safety region in the combustion vibration prediction device according to the sixth embodiment of the present invention.

【図19】従来例を示す構成図である。FIG. 19 is a configuration diagram showing a conventional example.

【符号の説明】[Explanation of symbols]

10 燃焼振動予測装置 1 内圧変動特性把握手段 2 燃焼振動領域推定手段 3 データベース 4 入力手段 5 出力手段 20 燃焼振動予測装置 12 内圧変動推定手段 13 データベース 14 入力手段 15 出力手段 30 プラント 40 燃焼振動予測装置 41 内圧変動特性把握手段 42 安全領域推定手段 43 データベース 44 入力手段 45 出力手段 46 NOx排出量特性把握手段 50 プラント 60 燃焼振動予測装置 62 安全領域推定手段 63 データベース 64 入力手段 65 出力手段 67 CO排出量特性把握手段 70 燃焼振動予測装置 71 内圧変動特性把握手段 72 安全領域推定手段 74 入力手段 76 排出量特性把握手段 78 フォーカス設定手段 80 燃焼振動予測装置 82 安全領域推定手段 84 入力手段 85 出力手段 88 フォーカス決定手段 89 調整案生成手段 DESCRIPTION OF SYMBOLS 10 Combustion vibration prediction apparatus 1 Internal pressure fluctuation characteristic grasping means 2 Combustion vibration area estimation means 3 Database 4 Input means 5 Output means 20 Combustion vibration prediction apparatus 12 Internal pressure fluctuation estimation means 13 Database 14 Input means 15 Output means 30 Plant 40 Combustion vibration prediction apparatus 41 Internal pressure fluctuation characteristic grasping means 42 Safety area estimating means 43 Database 44 Input means 45 Output means 46 NOx emission characteristic grasping means 50 Plant 60 Combustion vibration prediction device 62 Safety area estimating means 63 Database 64 Input means 65 Output means 67 CO emission Characteristic grasping means 70 Combustion vibration predicting device 71 Internal pressure fluctuation characteristic grasping means 72 Safety area estimating means 74 Input means 76 Emission characteristic grasping means 78 Focus setting means 80 Combustion vibration predicting device 82 Safety area estimating means 84 Input means 85 Output Means 88 focus determination unit 89 proposed adjustment generating unit

Claims (13)

【特許請求の範囲】[Claims] 【請求項1】 プラントデータ、気象データから内圧変
動を説明する数式モデルを構築し、構築した数式モデル
に基づいて燃焼振動の発生し易い領域、発生しにくい領
域を求めて出力することを特徴とする燃焼振動予測装
置。
1. A formula model for explaining internal pressure fluctuation is constructed from plant data and meteorological data, and a region where combustion oscillation is likely to occur and a region where combustion vibration is unlikely to occur are obtained and output based on the constructed formula model. Combustion vibration prediction device.
【請求項2】 プラントデータや気象データ、内圧変動
の制限値などを入力する入力手段と、 入力されたプラントデータ、気象データから燃焼器の内
圧変動を数式モデル化する内圧変動特性把握手段と、 前記内圧変動特性把握手段により求められた数式モデル
に内圧変動の制限値を適用して燃焼振動の発生し易い領
域を求める燃焼振動領域推定手段と、 前記燃焼振動領域推定手段による燃焼振動領域推定結果
を出力する出力手段と、 を有していることを特徴とする燃焼振動予測装置。
2. Input means for inputting plant data, weather data, a limit value of internal pressure fluctuation, etc., internal pressure fluctuation characteristic grasping means for mathematically modeling the internal pressure fluctuation of a combustor from the input plant data and weather data, A combustion oscillation region estimating unit that obtains a region where combustion oscillation is likely to occur by applying a limit value of the internal pressure variation to the mathematical model obtained by the internal pressure variation characteristic grasping unit; A combustion vibration predicting device, comprising: output means for outputting the following.
【請求項3】 前記入力手段により入力されたプラント
データや気象データなどを時系列に格納するデータベー
スを有し、前記内圧変動特性把握手段は前記データベー
スよりデータを取得して燃焼器の内圧変動を数式モデル
化することを特徴とする請求項2に記載の燃焼振動予測
装置。
3. A data base for storing plant data, weather data, and the like input by the input means in time series, wherein the internal pressure fluctuation characteristic grasping means obtains data from the database and detects the internal pressure fluctuation of the combustor. The combustion vibration predicting apparatus according to claim 2, wherein the apparatus is modeled by a mathematical model.
【請求項4】 プラントデータや気象データなどを入力
する入力手段と、 入力されたプラントデータや気象データより燃焼器の内
圧変動を推定する内圧変動推定手段と、 前記内圧変動推定手段により推定された内圧変動推定結
果を出力する出力手段と、 を有していることを特徴とする燃焼振動予測装置。
4. An input means for inputting plant data or weather data, an internal pressure fluctuation estimating means for estimating an internal pressure fluctuation of a combustor from the input plant data or weather data, and an internal pressure fluctuation estimating means. A combustion vibration prediction device, comprising: output means for outputting an internal pressure fluctuation estimation result.
【請求項5】 前記入力手段により入力されたプラント
データや気象データなどを時系列に格納するデータベー
スを有し、前記内圧変動推定手段は、データベースに格
納されている最新時刻のデータにより内圧変動の予想値
を推定することを特徴とする請求項4に記載の燃焼振動
予測装置。
5. A database for storing plant data, weather data, and the like input by the input unit in a time series, wherein the internal pressure fluctuation estimating unit detects the internal pressure fluctuation based on the latest time data stored in the database. The combustion vibration prediction device according to claim 4, wherein the prediction value is estimated.
【請求項6】 プラントデータ、気象データから内圧変
動、NOx排出量を説明する数式モデルを構築し、構築
した数式モデルに基づいて燃焼振動の発生しにくい領
域、発生し易い領域を求めて出力することを特徴とする
燃焼振動予測装置。
6. A mathematical model for explaining internal pressure fluctuations and NOx emissions is constructed from plant data and weather data, and a region in which combustion oscillation is less likely to occur and a region in which combustion vibration is more likely are determined and output based on the constructed mathematical model. A combustion vibration prediction device, characterized in that:
【請求項7】 プラントデータや気象データ、内圧変動
の制限値やNOxの規制値などを入力する入力手段と、 入力されたプラントデータ、気象データから燃焼器の内
圧変動を数式モデル化する内圧変動特性把握手段と、 入力されたプラントデータ、気象データからNOx排出
量を数式モデル化するNOx排出量特性把握手段と、 前記内圧変動特性把握手段により求められた数式モデル
に内圧変動の制限値を適用し、前記NOx排出量特性把
握手段により求められた数式モデルにNOxの規制値を
適用して、NOx排出量が規制値以下でかつ燃焼振動の
発生しにくい領域を求める安全領域推定手段と、 前記安全領域推定手段による安全領域推定結果を出力す
る出力手段と、 を有していることを特徴とする燃焼振動予測装置。
7. An input means for inputting plant data, weather data, a limit value of internal pressure fluctuation, a NOx control value, etc., and an internal pressure fluctuation for formulating an internal pressure fluctuation of a combustor from the input plant data and weather data. Characteristic grasping means; NOx emission characteristic grasping means for mathematically modeling NOx emissions from the input plant data and weather data; and internal pressure fluctuation limit value applied to the mathematical expression model obtained by the internal pressure fluctuation characteristic grasping means. A safety region estimating means for applying a NOx regulation value to the mathematical model obtained by the NOx emission characteristic grasping means to find a region where the NOx emission amount is equal to or less than the regulation value and in which combustion vibration is unlikely to occur; Output means for outputting a safety area estimation result obtained by the safety area estimation means.
【請求項8】 プラントデータ、気象データから内圧変
動、NOx、CO排出量を説明する数式モデルを構築
し、構築した数式モデルに基づいて燃焼振動の発生しに
くい領域、発生し易い領域を求めて出力することを特徴
とする燃焼振動予測装置。
8. A mathematical model for explaining internal pressure fluctuation, NOx, and CO emissions is constructed from plant data and weather data, and a region in which combustion oscillation is unlikely to occur and a region in which combustion vibration is likely to occur are determined based on the constructed mathematical model. A combustion vibration prediction device characterized by outputting.
【請求項9】 プラントデータや気象データ、内圧変動
の制限値やNOx、COの規制値などを入力する入力手
段と、 入力されたプラントデータ、気象データから燃焼器の内
圧変動を数式モデル化する内圧変動特性把握手段と、 入力されたプラントデータ、気象データからNOx排出
量を数式モデル化するNOx排出量特性把握手段と、 入力されたプラントデータ、気象データからCO排出量
を数式モデル化するCO排出量特性把握手段と、 前記内圧変動特性把握手段により求められた数式モデル
に内圧変動の制限値を適用し、前記NOx排出量特性把
握手段により求められた数式モデルにNOxの規制値を
適用し、前記CO排出量特性把握手段により求められた
数式モデルにCOの規制値を適用して、NOx及びCO
排出量が規制値以下でかつ燃焼振動の発生しにくい領域
を求める安全領域推定手段と、 前記安全領域推定手段による安全領域推定結果を出力す
る出力手段と、 を有していることを特徴とする燃焼振動予測装置。
9. An input means for inputting plant data, weather data, internal pressure fluctuation limit values, NOx, CO regulation values, and the like, and formulating an internal pressure fluctuation of a combustor from input plant data and weather data. Internal pressure fluctuation characteristic grasping means, NOx emission characteristic grasping means for mathematically modeling NOx emissions from input plant data and weather data, and CO for mathematically modeling CO emissions from input plant data and weather data Emission amount characteristic grasping means, a limit value of internal pressure fluctuation is applied to the mathematical model obtained by the internal pressure fluctuation characteristic grasping means, and a NOx regulation value is applied to the mathematical expression model obtained by the NOx emission characteristic grasping means. By applying the regulated value of CO to the mathematical model obtained by the CO emission characteristic grasping means, NOx and CO
A safety area estimating means for obtaining an area in which the emission amount is equal to or less than a regulation value and in which combustion oscillation is unlikely to occur, and an output means for outputting a safety area estimation result by the safety area estimating means. Combustion vibration prediction device.
【請求項10】 プラントデータや気象データ、内圧変
動の制限値やNOx、COの規制値などを入力する入力
手段と、 入力されたプラントデータ、気象データから数式モデル
化に用いるデータを選択するフォーカス設定手段と、 選択されたプラントデータ、気象データから燃焼器の内
圧変動を数式モデル化する内圧変動特性把握手段と、 選択されたプラントデータ、気象データからNOx及び
CO排出量を数式モデル化する排出量特性把握手段と、 前記内圧変動特性把握手段により求められた数式モデル
に内圧変動の制限値を適用し、前記排出量特性把握手段
により求められた数式モデルにNOx及びCOの規制値
を適用して、NOx及びCO排出量が規制値以下でかつ
燃焼振動の発生しにくい領域を求める安全領域推定手段
と、 前記安全領域推定手段による安全領域推定結果を出力す
る出力手段と、 を有していることを特徴とする燃焼振動予測装置。
10. Input means for inputting plant data, weather data, internal pressure fluctuation limit values, NOx, CO control values, etc., and a focus for selecting data to be used for mathematical modeling from the input plant data and weather data. Setting means; internal pressure fluctuation characteristic grasping means for mathematically modeling the internal pressure fluctuation of the combustor from the selected plant data and weather data; and emissions for mathematically modeling the NOx and CO emissions from the selected plant data and weather data. Volume characteristic grasping means, applying a limit value of internal pressure fluctuation to the mathematical model obtained by the internal pressure fluctuation characteristic grasping means, and applying NOx and CO regulation values to the mathematical expression model obtained by the emission characteristic grasping means. Safety area estimating means for finding an area where NOx and CO emissions are less than or equal to a regulation value and in which combustion oscillation is unlikely to occur; Combustion oscillation prediction apparatus characterized in that it has an output means for outputting the safety region estimation result by the region estimation means.
【請求項11】 前記フォーカス設定手段は、前記入力
手段により指定された領域、または設定モードに基づい
て、前記入力手段から入力されたプラントデータ、気象
データを選択することを特徴とする請求項10に記載の
燃焼振動予測装置。
11. The apparatus according to claim 10, wherein said focus setting means selects plant data and weather data input from said input means based on an area designated by said input means or a setting mode. Combustion vibration prediction apparatus according to any one of the above.
【請求項12】 プラントデータや気象データ、内圧変
動の制限値やNOx、COの規制値などを入力する入力
手段と、 入力されたプラントデータ、気象データから数式モデル
化に用いるデータを選択するフォーカス決定手段と、 選択されたプラントデータ、気象データから燃焼器の内
圧変動を数式モデル化する内圧変動特性把握手段と、 選択されたプラントデータ、気象データからNOx及び
CO排出量を数式モデル化する排出量特性把握手段と、 前記内圧変動特性把握手段により求められた数式モデル
に内圧変動の制限値を適用し、前記排出量特性把握手段
により求められた数式モデルにNOx及びCOの規制値
を適用して、NOx及びCO排出量が規制値以下でかつ
燃焼振動の発生しにくい領域を求める安全領域推定手段
と、 前記安全領域推定手段による安全領域推定結果を用い
て、次に計測するべきポイントを求める調整案生成手段
と、 前記安全領域推定手段による安全領域推定結果及び前記
調整案生成手段による計測するべきポイントを出力する
出力手段と、 を有していることを特徴とする燃焼振動予測装置。
12. Input means for inputting plant data, weather data, internal pressure fluctuation limit values, NOx, CO control values, etc., and a focus for selecting data to be used in mathematical modeling from the input plant data and weather data. Determining means; means for grasping the internal pressure fluctuation characteristics of the combustor from the selected plant data and weather data as a mathematical model; and emissions for modeling the NOx and CO emissions from the selected plant data and the weather data. Volume characteristic grasping means, applying a limit value of internal pressure fluctuation to the mathematical model obtained by the internal pressure fluctuation characteristic grasping means, and applying NOx and CO regulation values to the mathematical expression model obtained by the emission characteristic grasping means. Safety area estimating means for finding an area where NOx and CO emissions are less than or equal to a regulation value and in which combustion oscillation is unlikely to occur; Adjustment plan generating means for obtaining a point to be measured next using the safety area estimation result by the area estimating means, and a safety area estimation result by the safety area estimating means and a point to be measured by the adjustment plan generating means are output. A combustion vibration prediction device, comprising: output means.
【請求項13】 前記フォーカス決定手段は、前回のフ
ォーカス決定により選択されたプラントデータ、気象デ
ータに基づいて得られた数式モデルに基づいて、次のフ
ォーカスを決定することを特徴とする請求項12に記載
の燃焼振動予測装置。
13. The apparatus according to claim 12, wherein said focus determining means determines a next focus based on a mathematical model obtained based on plant data and weather data selected by a previous focus determination. Combustion vibration prediction apparatus according to any one of the above.
JP2001022338A 2000-06-01 2001-01-30 Combustion vibration prediction device Expired - Lifetime JP3999467B2 (en)

Priority Applications (5)

Application Number Priority Date Filing Date Title
JP2001022338A JP3999467B2 (en) 2000-06-01 2001-01-30 Combustion vibration prediction device
US10/057,955 US6922612B2 (en) 2001-01-30 2002-01-29 Combustion vibration estimating apparatus, plant and gas turbine plant
CA002369539A CA2369539C (en) 2001-01-30 2002-01-29 Combustion vibration estimating apparatus, plant and gas turbine plant
DE60210157T DE60210157T2 (en) 2001-01-30 2002-01-30 Device for estimating combustion chamber pressure fluctuations, plant and gas turbine plant
EP02002276A EP1229226B1 (en) 2001-01-30 2002-01-30 Combustion vibration estimating apparatus, plant and gas turbine plant

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2000164251 2000-06-01
JP2000-164251 2000-06-01
JP2001022338A JP3999467B2 (en) 2000-06-01 2001-01-30 Combustion vibration prediction device

Publications (2)

Publication Number Publication Date
JP2002054460A true JP2002054460A (en) 2002-02-20
JP3999467B2 JP3999467B2 (en) 2007-10-31

Family

ID=26593150

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2001022338A Expired - Lifetime JP3999467B2 (en) 2000-06-01 2001-01-30 Combustion vibration prediction device

Country Status (1)

Country Link
JP (1) JP3999467B2 (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006511750A (en) * 2002-12-23 2006-04-06 ヌオーヴォ ピニォーネ ホールディング ソシエタ ペル アチオニ A system that evaluates and controls the concentration of pollutant gases at the gas turbine outlet.
JP2010249133A (en) * 2009-04-10 2010-11-04 General Electric Co <Ge> Corrected megawatt backup curve methodology
JP2012233480A (en) * 2011-05-05 2012-11-29 Alstom Technology Ltd Method for protecting gas turbine engine against high dynamical process value and gas turbine engine for conducting the method
CN103807873A (en) * 2014-02-25 2014-05-21 北京华清燃气轮机与煤气化联合循环工程技术有限公司 Active control method for burning instability of combustion gas turbine combustor
JP2014169701A (en) * 2014-05-20 2014-09-18 Mitsubishi Heavy Ind Ltd Method of controlling gas turbine and method of setting fuel-air ratio of gas turbine
JP2017508095A (en) * 2014-01-22 2017-03-23 シーメンス アクティエンゲゼルシャフト Gas turbine adjustment method

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101833835B1 (en) * 2016-08-26 2018-03-02 한국중부발전(주) Gas turbine operation margin prediction system using weather fore cast

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006511750A (en) * 2002-12-23 2006-04-06 ヌオーヴォ ピニォーネ ホールディング ソシエタ ペル アチオニ A system that evaluates and controls the concentration of pollutant gases at the gas turbine outlet.
JP2010249133A (en) * 2009-04-10 2010-11-04 General Electric Co <Ge> Corrected megawatt backup curve methodology
US9115648B2 (en) 2009-04-10 2015-08-25 General Electric Company Corrected megawatt backup curve methodology
JP2012233480A (en) * 2011-05-05 2012-11-29 Alstom Technology Ltd Method for protecting gas turbine engine against high dynamical process value and gas turbine engine for conducting the method
JP2017508095A (en) * 2014-01-22 2017-03-23 シーメンス アクティエンゲゼルシャフト Gas turbine adjustment method
US10309320B2 (en) 2014-01-22 2019-06-04 Siemens Aktiengesellschaft Method for regulating a gas turbine
CN103807873A (en) * 2014-02-25 2014-05-21 北京华清燃气轮机与煤气化联合循环工程技术有限公司 Active control method for burning instability of combustion gas turbine combustor
JP2014169701A (en) * 2014-05-20 2014-09-18 Mitsubishi Heavy Ind Ltd Method of controlling gas turbine and method of setting fuel-air ratio of gas turbine

Also Published As

Publication number Publication date
JP3999467B2 (en) 2007-10-31

Similar Documents

Publication Publication Date Title
US20200184131A1 (en) A method for prediction of key performance parameter of an aero-engine transition state acceleration process based on space reconstruction
US6922612B2 (en) Combustion vibration estimating apparatus, plant and gas turbine plant
KR920003499B1 (en) Applied numerical control system
US5559690A (en) Residual activation neural network
US7058617B1 (en) Method and apparatus for training a system model with gain constraints
US20170255177A1 (en) Control Parameter Automatic-Adjustment Apparatus, Control Parameter Automatic-Adjustment Method, and Control Parameter Automatic-Adjustment Apparatus Network
US7188019B2 (en) Gas turbine control apparatus, gas turbine system and gas turbine control method
US6278899B1 (en) Method for on-line optimization of a plant
US8214062B2 (en) Plant control system and thermal power generation plant control system
CN111290271B (en) Control parameter optimization system and operation control optimization device provided with same
KR101256038B1 (en) Gas turbine control method and controller
US20160138470A1 (en) Application of probabilistic control in liquid-fueled gas turbine tuning, related control systems, computer program products and methods
KR0149864B1 (en) Computer aided turning of turbine controls
CN101482062A (en) Methods and systems for providing real-time comparision with an alternate control strategy for a turbine
US20140005909A1 (en) Real time linearization of a component-level gas turbine engine model for model-based control
CN101802728A (en) Gaseous emission is carried out system and method based on the virtual sensing of experience set
WO2005054968A1 (en) Intelligent system for detection of process status, process fault and preventive maintenance
US6985781B2 (en) Residual activation neural network
RU2675965C2 (en) Method for determining emissions mode
JP2012141862A (en) Control device for plant and control device for fire power generation plant
CN109992880A (en) A kind of gas turbine performance adaptation method, equipment and system
JP2017110644A (en) Application of combined probabilistic control in gas turbine tuning for power output-emissions parameters with scaling factor, and related control systems, computer program products and methods
JP2002054460A (en) Combustion vibration predicting device
CN106050435B (en) Regulation and control system, computer program product and method for a gas turbine
JP2017125495A (en) Machine-specific probabilistic control in gas turbine tuning for power output-emissions parameters, related control systems, computer program products and methods

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20041217

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20060519

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20060613

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20060809

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20061114

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20070109

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20070501

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20070531

A911 Transfer to examiner for re-examination before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20070706

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20070724

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20070809

R151 Written notification of patent or utility model registration

Ref document number: 3999467

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100817

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100817

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110817

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110817

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120817

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130817

Year of fee payment: 6

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

EXPY Cancellation because of completion of term