WO2023218930A1

WO2023218930A1 - ガスタービン制御装置、ガスタービン制御方法、及びプログラム

Info

Publication number: WO2023218930A1
Application number: PCT/JP2023/016133
Authority: WO
Inventors: 一茂高木; 竜児竹中; 皓士郎福本; 修平戸田; 伸明林田
Original assignee: 三菱重工業株式会社
Priority date: 2022-05-09
Filing date: 2023-04-24
Publication date: 2023-11-16
Also published as: JP2023166083A

Abstract

ガスタービン制御装置は、燃焼器に流入する燃料の量を示す燃料流量を算出する燃料流量算出部と、圧縮機が吸入する空気の量を示す空気流量を算出する空気流量算出部と、燃料流量と、空気流量と、燃焼器に関する熱エネルギ収支に関する物理モデル式とに基づいてタービン入口温度を算出するタービン入口温度算出部と、ガスタービンのエネルギ収支に関する計測データに基づいてガスタービンのヒートバランスを満たす演算を行って演算結果を出力するヒートバランス演算部と、ガスタービンの状態を検出する状態検出部と、状態検出部が検出するガスタービンの状態が静定状態の場合にヒートバランス演算部が出力する演算結果に基づいて、燃料流量と空気流量のいずれか一方、または、両方を補正してタービン入口温度算出部に与える補正処理部と、を備える。

Description

ガスタービン制御装置、ガスタービン制御方法、及びプログラム

　本開示は、ガスタービン制御装置、ガスタービン制御方法、及びプログラムに関する。
　本願は、２０２２年５月９日に、日本に出願された特願２０２２－０７６８４８号に基づき優先権を主張し、その内容をここに援用する。

　ガスタービンでは、例えば、圧縮機で生成された圧縮空気を取り込んだ燃焼器において燃料を噴射して燃焼ガスを生成する。生成した燃焼ガスがタービンを駆動することにより動力が生成される。燃焼器に燃料を供給する燃料供給系統は、燃焼効率や燃焼安定性の観点から、複数の燃料供給系統に分けられることがある。複数の燃料供給系統の各々に対して燃料を配分する比率、すなわち燃料配分比は、タービン入口温度に基づいて算出することが可能である。

　ただし、タービン入口温度を安定して計測できる計器が存在しないため、ガスタービンに備えられたセンサが計測することにより得られる計測データに基づいてタービン入口温度を推定し、推定したタービン入口温度を用いて燃料配分比を算出することが一般的に行われる。そのため、推定するタービン入口温度の精度を担保することが重要となる。例えば、特許文献１には、燃焼器に流入する燃料流量や空気流量から燃焼後のタービン入口温度を推定し、推定したタービン入口温度をセンサによって計測したタービン出口温度、すなわち排ガス温度を用いて補正することにより、タービン入口温度の精度を高める技術が開示されている。

日本国特許第６９２０４５６号公報

　しかしながら、タービン入口温度と排ガス温度の関係であるタービン効率を一意に定めることは難しい。そのため、特許文献１に開示されている技術には、タービン入口温度を排ガス温度によって補正しても、タービン入口温度の精度を高めるように補正できない場合があるという課題がある。また、タービン効率を一意に定めることができたとしても、ガスタービン１０が長期間運転されることにより、経年劣化等によってガスタービン１０の性能が変化する。そのため、長期間にわたって、一意に定めたタービン効率を使い続けることができないため、安定してタービン入口温度の精度を維持することが難しくなるという課題がある。

　本開示は、上記課題を解決するためになされたものであって、タービン効率を用いずに、長期間にわたってタービン入口温度の精度を維持することができるガスタービン制御装置、ガスタービン制御方法、及びプログラムを提供することを目的とする。

　上記課題を解決するために、本開示に係るガスタービン制御装置は、ガスタービンの燃焼器に燃料を供給する燃料供給系統に関する計測データに基づいて前記燃焼器に流入する燃料の量を示す燃料流量を算出する燃料流量算出部と、前記ガスタービンの圧縮機に関する計測データに基づいて前記圧縮機が吸入する空気の量を示す空気流量を算出する空気流量算出部と、前記燃料流量と、前記空気流量と、前記燃焼器に関する熱エネルギ収支に関する物理モデル式とに基づいてタービン入口温度を算出するタービン入口温度算出部と、前記ガスタービンのエネルギ収支に関する計測データに基づいて前記ガスタービンのヒートバランスを満たす演算を行って演算結果を出力するヒートバランス演算部と、前記ガスタービンの状態を検出する状態検出部と、前記状態検出部が検出する前記ガスタービンの状態が静定状態の場合に前記ヒートバランス演算部が出力する前記演算結果に基づいて、前記燃料流量と前記空気流量のいずれか一方、または、両方を補正して前記タービン入口温度算出部に与える補正処理部と、を備える。

　本開示に係るガスタービン制御方法は、ガスタービンの燃焼器に燃料を供給する燃料供給系統に関する計測データに基づいて前記燃焼器に流入する燃料の量を示す燃料流量を算出するステップと、前記ガスタービンの圧縮機に関する計測データに基づいて前記圧縮機が吸入する空気の量を示す空気流量を算出するステップと、前記燃料流量と、前記空気流量と、前記燃焼器に関する熱エネルギ収支に関する物理モデル式とに基づいてタービン入口温度を算出するステップと、前記ガスタービンのエネルギ収支に関する計測データに基づいて前記ガスタービンのヒートバランスを満たす演算を行って演算結果を出力するステップと、前記ガスタービンの状態を検出するステップと、前記状態を検出するステップによって検出された前記ガスタービンの状態が静定状態の場合に前記ヒートバランスを満たす演算を行うステップにおいて出力された前記演算結果に基づいて、前記燃料流量と前記空気流量のいずれか一方、または、両方を補正して前記タービン入口温度を算出するステップに与えるステップと、を含む。

　本開示に係るプログラムは、コンピュータを、ガスタービンの燃焼器に燃料を供給する燃料供給系統に関する計測データに基づいて前記燃焼器に流入する燃料の量を示す燃料流量を算出する燃料流量算出手段、前記ガスタービンの圧縮機に関する計測データに基づいて前記圧縮機が吸入する空気の量を示す空気流量を算出する空気流量算出手段、前記燃料流量と、前記空気流量と、前記燃焼器に関する熱エネルギ収支に関する物理モデル式とに基づいてタービン入口温度を算出するタービン入口温度算出手段、前記ガスタービンのエネルギ収支に関する計測データに基づいて前記ガスタービンのヒートバランスを満たす演算を行って演算結果を出力するヒートバランス演算手段、前記ガスタービンの状態を検出する状態検出手段、前記状態検出手段が検出する前記ガスタービンの状態が静定状態の場合に前記ヒートバランス演算手段が出力する前記演算結果に基づいて、前記燃料流量と前記空気流量のいずれか一方、または、両方を補正して前記タービン入口温度算出手段に与える補正処理手段、として機能させるためのプログラムである。

　本開示のガスタービン制御装置、ガスタービン制御方法、及びプログラムによれば、タービン効率を用いずに、長期間にわたってタービン入口温度の精度を維持することができる。

本開示の第１の実施形態に係るガスタービン発電プラントの構成例を示す模式図である。本開示の第１の実施形態に係るガスタービン制御装置の構成例を示すブロック図である。本開示の第１の実施形態に係るガスタービン制御装置の動作例を示すフローチャート（その１）である。本開示の第１の実施形態に係るガスタービン制御装置の動作例を示すフローチャート（その２）である。本開示の第１の実施形態に係るガスタービン制御装置の動作例を示すフローチャート（その３）である。本開示の第１の実施形態の他の構成例に係るガスタービン制御装置の動作例を示すフローチャート（その１）である。本開示の第１の実施形態の他の構成例に係るガスタービン制御装置の動作例を示すフローチャート（その２）である。本開示の第２の実施形態に係るガスタービン発電プラントの構成例を示す模式図である。本開示の第２の実施形態に係るガスタービン制御装置の構成例を示すブロック図である。本開示の第２の実施形態に係るガスタービン制御装置の動作例を示すフローチャート（その１）である。本開示の第２の実施形態に係るガスタービン制御装置の動作例を示すフローチャート（その２）である。本開示の第２の実施形態に係るガスタービン制御装置の動作例を示すフローチャート（その３）である。本開示の第３の実施形態に係るガスタービン制御装置の構成例を示すブロック図である。本開示の第３の実施形態に係るガスタービン制御装置の動作例を示すフローチャート（その１）である。本開示の第３の実施形態に係るガスタービン制御装置の動作例を示すフローチャート（その２）である。少なくとも１つの実施形態に係るコンピュータの構成を示す概略ブロック図である。

　以下、本開示の実施形態に係るガスタービン制御装置、ガスタービン制御方法、及びプログラムについて、図１～図１６を参照して説明する。図１、図８は、それぞれ、本開示の第１及び第２の実施形態に係るガスタービン発電プラント１００，１００ａの構成例を示す模式図である。図２，図９，図１３は、それぞれ、本開示の第１から第３の実施形態に係るガスタービン制御装置２０，２０ａ，２０ｂの構成例を示すブロック図である。図３～図５は、本開示の第１の実施形態に係るガスタービン制御装置２０の動作例を示すフローチャートである。図６、図７は、本開示の第１の実施形態の他の構成例に係るガスタービン制御装置２０の動作例を示すフローチャートである。図１０～図１２は、本開示の第２の実施形態に係るガスタービン制御装置２０ａの動作例を示すフローチャートである。図１４，図１５は、本開示の第３の実施形態に係るガスタービン制御装置２０ｂの動作例を示すフローチャートである。図１６は、本開示の少なくとも１つの実施形態に係るコンピュータの構成を示す概略ブロック図である。なお、各図において同一または対応する構成には同一の符号を用いて説明を適宜省略する。
　なお、本開示に係る各実施形態においては、図２、図９等に示すごとく、全ての計算が単一のガスタービン制御装置内で行われるものとして説明するが、これはあくまで構成の一例にすぎない。すなわち、他の実施形態においては、例えば、燃料流量・タービン入口温度算出はガスタービン制御装置で行い、ヒートバランス計算や補正係数の算出は異なる計算機で行う態様であってもよい。その他、各計算処理と、当該計算処理を実行するハードウェアとの組み合わせは適宜変更可能である。

＜第１の実施形態＞
（第１の実施形態のガスタービン発電プラントの構成例）
　図１は、第１の実施形態のガスタービン発電プラント１００の全体構成図である。ガスタービン発電プラント１００は、ガスタービン１０、ガスタービン１０を制御するガスタービン制御装置２０、燃料を供給する燃料供給装置３０、燃料供給装置３０が供給する燃料をガスタービン１０に導く燃料供給系統４０、ガスタービン１０の動力で駆動されることにより発電を行う発電機５０、及び、様々な状態量を計測し、計測により検出する計測データをガスタービン制御装置２０に送信する各種のセンサ２５ａ～２５ｌを備える。

　ガスタービン１０は、圧縮機１１、燃焼器１２、及びタービン１３を備える。圧縮機１１は、空気を吸入し、吸入した空気を圧縮することにより、高圧の圧縮空気を生成する。圧縮機１１で生成された圧縮空気は下流側にある燃焼器１２及びタービン１３に対して供給される。圧縮機１１は、主軸Ａｓを中心軸として回転する圧縮機ロータ１１１と、圧縮機ロータ１１１を外周側から覆う圧縮機車室１１２と、圧縮機１１が吸入する空気量を調整するためのインレットガイドベーン１４（ＩＧＶ(Inlet Guide Vane)）とを備える。インレットガイドベーン１４は、点線の矢印で示すように通信回線などの制御回線によりガスタービン制御装置２０に接続されており、ガスタービン制御装置２０からの指令値を受けたインレットガイドベーン１４により圧縮機１１が吸入する空気量が調整される。

　タービン１３は、主軸Ａｓを中心軸として回転するタービンロータ１３１と、タービンロータ１３１を外周側から覆うタービン車室１３２とを備える。圧縮機ロータ１１１と、タービンロータ１３１とは、連結されており、主軸Ａｓを中心軸として一体的に回転する。燃焼器１２は、圧縮機１１によって生成された圧縮空気に対して、燃料供給系統４０から供給される燃料を噴射して燃焼させることにより、高圧力、高温度の燃焼ガスを生成する。燃焼器１２は、燃焼器内筒１２１と、燃焼器内筒１２１を覆う燃焼器車室１６とを備える。燃焼器内筒１２１は、燃料供給系統４０から供給される燃料を噴射して点火することにより、燃焼器車室１６内の燃料を燃焼させる。燃焼器車室１６は、圧縮機車室１１２と、タービン車室１３２とに連結されており、燃焼器車室１６で生成された燃焼ガスが、タービン車室１３２に供給されて、タービンロータ１３１が回転し、タービンロータ１３１の回転により圧縮機ロータ１１１が一体的に回転して空気を吸入する。

　燃料供給系統４０は、燃料供給装置３０が供給する燃料を、燃料配管を通じて所定の圧力及び流量で燃焼器１２に供給する。燃料供給系統４０は、図１において、実線、または、破線の矢印で示す燃料配管であって燃料供給装置３０とガスタービン１０の間を接続する燃料配管の他に、燃料配管に設けられる流量調整弁（以下、流調弁という）４１と、ノズル４２と、マニホールド配管４３とを備える。

　燃料供給系統４０は、燃料供給系統４０－１、燃料供給系統４０－２、及び燃料供給系統４０－３の３つの燃料供給系統を含んでいる。燃料供給装置３０は、供給する燃料を３つに分岐し、分岐した燃料の各々を燃料供給系統４０－１、燃料供給系統４０－２、及び燃料供給系統４０－３に供給する。例えば、燃料供給系統４０－１は、燃焼器１２のトップハット部に燃料を供給するためのトップハット燃料供給系統であり、燃料供給系統４０－２は燃焼器内筒１２１の中心部に燃料を供給するパイロット燃料供給系統であり、燃料供給系統４０－３は燃焼器内筒１２１の中心部を囲む部分に燃料を供給するメイン燃料供給系統である。

　燃料供給系統４０が３つの燃料供給系統４０－１～４０－３を含んでいることから、流調弁４１は、燃料供給系統４０－１～４０－３の各々に対応する流調弁４１－１，４１－２，４１－３を有する。流調弁４１－１，４１－２，４１－３の各々は、点線の矢印で示すように通信回線などの制御回線によりガスタービン制御装置２０に接続されている。流調弁４１－１，４１－２，４１－３の各々は、制御回線を通じてガスタービン制御装置２０から各々に対する弁開度を示す指令値を受けると、受けた指令値にしたがった開度で弁の開度を大きくしたり小さくしたりする。

　流調弁４１と同様に、ノズル４２は、燃料供給系統４０－１～４０－３の各々に対応するノズル４２－１，４２－２，４２－３を有する。マニホールド配管４３については、図１では、紙面の都合上、燃料供給系統４０－１に対応するマニホールド配管４３のみを示しているが、実際には、燃料供給系統４０－１～４０－３ごとに、１つずつマニホールド配管４３が設けられている。以下、燃料供給系統４０－１～４０－３の各々に対応するマニホールド配管４３を個別に説明する必要がある場合、符号の枝番号を対応させて、それぞれマニホールド配管４３－１～４３－３という。マニホールド配管４３－１～４３－３の各々は、燃焼器１２に接続されており、例えば、マニホールド配管４３－１は、燃焼器１２のトップハット部に燃料を供給し、マニホールド配管４３－２は、燃焼器内筒１２１の中心部に燃料を供給し、マニホールド配管４３－３は、燃焼器内筒１２１の中心部を囲む部分に燃料を供給する。

　ロータ１５は、圧縮機ロータ１１１に連結されており、圧縮機ロータ１１１及びタービンロータ１３１と共に、主軸Ａｓを中心軸として一体的に回転する。発電機５０は、ロータ１５を介して圧縮機ロータ１１１の一端に接続される。発電機５０は、ロータ１５の回転によって駆動されて発電を行って電力を生成する。

　次に、ガスタービン制御装置２０に接続する各種のセンサ２５ａ～２５ｌについて説明する。各種のセンサ２５ａ～２５ｌの各々は、点線の矢印で示すように通信回線などの制御回線によりガスタービン制御装置２０に接続されており、計測により検出する計測データを、制御回線を通じてガスタービン制御装置２０に送信する。

　流調弁前圧力センサ２５ａ、流調弁後圧力センサ２５ｂ、及び燃料温度センサ２５ｃは、以下に示すような、燃料供給系統４０に関する計測データを検出する。流調弁前圧力センサ２５ａは、燃料供給系統４０－１，４０－２，４０－３ごとに備えられる流調弁前圧力センサ２５ａ－１，２５ａ－２，２５ａ－３を有する。流調弁前圧力センサ２５ａ－１，２５ａ－２，２５ａ－３の各々は、各々に対応する流調弁４１－１，４１－２，４１－３に流入する燃料の圧力である流調弁前圧力Ｐ_１ＦＶ（１），Ｐ_１ＦＶ（２），Ｐ_１ＦＶ（３）を検出する。流調弁後圧力センサ２５ｂは、燃料供給系統４０－１，４０－２，４０－３ごとに備えられる流調弁後圧力センサ２５ｂ－１，２５ｂ－２，２５ｂ－３を有する。流調弁後圧力センサ２５ｂ－１，２５ｂ－２，２５ｂ－３の各々は、各々に対応する流調弁４１－１，４１－２，４１－３の各々から流出する燃料の圧力である流調弁後圧力Ｐ_２ＦＶ（１），Ｐ_２ＦＶ（２），Ｐ_２ＦＶ（３）を検出する。

　燃料温度センサ２５ｃは、燃料供給系統４０－１，４０－２，４０－３ごとに備えられる燃料温度センサ２５ｃ－１，２５ｃ－２，２５ｃ－３を有する。燃料温度センサ２５ｃ－１，２５ｃ－２，２５ｃ－３の各々は、各々に対応する燃料供給系統４０－１，４０－２，４０－３の燃料配管を流れる燃料の温度である燃料温度Ｔ_ｆ（１），Ｔ_ｆ（２），Ｔ_ｆ（３）を検出する。

　電力計センサ２５ｄ、インデックス差圧センサ２５ｅ、圧縮機入口圧力センサ２５ｆ、圧縮機入口温度センサ２５ｇ、圧縮機出口圧力センサ２５ｈ、圧縮機出口温度センサ２５ｉ、燃焼器車室温度センサ２５ｊ、排ガス圧力センサ２５ｋ、及び排ガス温度センサ２５ｌは、以下に示すような、ガスタービンに関する計測データを検出する。電力計センサ２５ｄは、発電機５０が発電する電力を計測し、計測した電力をガスタービン出力Ｇ_ｏｕｔとして検出する。

　インデックス差圧センサ２５ｅは、圧縮機インデックス差圧Ｐ_{ｉｎｄｅｘ}を検出する。ここで、圧縮機インデックス差圧Ｐ_{ｉｎｄｅｘ}とは、圧縮機１１のケーシング吸込口における圧力と、圧縮機１１の内部の翼付近との圧力差であり、圧縮機１１が吸入する空気流量の指標となる値である。圧縮機入口圧力センサ２５ｆは、圧縮機１１のケーシング吸込口における圧力である圧縮機入口圧力Ｐ_１Ｃを検出する。圧縮機入口温度センサ２５ｇは、圧縮機１１のケーシング吸込口における温度である圧縮機入口温度Ｔ_１Ｃを検出する。圧縮機出口圧力センサ２５ｈは、圧縮機１１の出口の圧力である圧縮機出口圧力Ｐ_２Ｃを検出する。圧縮機出口温度センサ２５ｉは、圧縮機１１の出口の温度である圧縮機出口温度Ｔ_２Ｃを検出する。燃焼器車室温度センサ２５ｊは、燃焼ガス温度、すなわち燃焼器車室１６内の空気の温度である燃焼器車室温度Ｔ_ＣＳを検出する。排ガス圧力センサ２５ｋは、タービン１３の出口の圧力である排ガス圧力Ｐ_２Ｔを検出する。排ガス温度センサ２５ｌは、タービン１３の出口の温度である排ガス温度Ｔ_２Ｔを検出する。

（第１の実施形態のガスタービン制御装置の構成例）
　ガスタービン制御装置２０は、図２に示すように、燃料流量算出部２０１、空気流量算出部２０２、ヒートバランス演算部２０３、補正処理部２０４、状態検出部２０５、タービン入口温度算出部２０６、燃料配分比算出部２０７、及び弁開度算出部２０８を備える。

　ガスタービン制御装置２０において、燃料流量算出部２０１は、予め定められる燃料流量算出関数ｆ_１（・）に対して、流調弁４１に対して指令値として与えられる弁開度Ｏ_１，Ｏ_２，Ｏ_３を示すデータの各々と、燃料供給系統４０－１，４０－２，４０－３の各々に関する計測データとを代入して単位時間当たりにガスタービン１０の燃焼器１２に流入する燃料の量を示す燃料流量Ｇ_ｆｔを算出する。ここで、燃料供給系統４０－１，４０－２，４０－３の各々に関する計測データとは、流調弁前圧力センサ２５ａが検出する流調弁前圧力Ｐ_１ＦＶ（１），Ｐ_１ＦＶ（２），Ｐ_１ＦＶ（３）、流調弁後圧力センサ２５ｂが検出する流調弁後圧力Ｐ_２ＦＶ（１），Ｐ_２ＦＶ（２），Ｐ_２ＦＶ（３）、及び燃料温度センサ２５ｃが検出する燃料温度Ｔ_ｆ（１），Ｔ_ｆ（２），Ｔ_ｆ（３）のデータである。これらのデータに基づいて、燃料流量算出部２０１は、次式（１），（２），（３）に示すように、燃料供給系統４０－１，４０－２，４０－３ごとに燃料流量算出関数ｆ_１（・）を用いた演算を行い、単位時間当たりに燃料供給系統４０－１，４０－２，４０－３の各々を通じて燃焼器１２に流入する燃料流量Ｇ_ｆｔ（１）と、燃料流量Ｇ_ｆｔ（２）と、燃料流量Ｇ_ｆｔ（３）とを算出する。

　燃料流量算出部２０１は、次式（４）に示すように、燃料供給系統４０－１に対応する燃料流量Ｇ_ｆｔ（１）と、燃料供給系統４０－２に対応する燃料流量Ｇ_ｆｔ（２）と、燃料供給系統４０－３に対応する燃料流量Ｇ_ｆｔ（３）とを合計して単位時間当たりに燃料供給系統４０を通じて燃焼器１２に流入する燃料流量Ｇ_ｆｔを算出する。

　空気流量算出部２０２は、ガスタービン１０の圧縮機１１に関する計測データ、すなわち、インデックス差圧センサ２５ｅが検出する圧縮機インデックス差圧Ｐ_{ｉｎｄｅｘ}と、圧縮機入口圧力センサ２５ｆが検出する圧縮機入口圧力Ｐ_１Ｃと、圧縮機入口温度センサ２５ｇが検出する圧縮機入口温度Ｔ_１Ｃとを、次式（５）に示すように、予め定められる空気流量算出関数ｆ_２（・）に代入して単位時間当たりにガスタービン１０の圧縮機１１が吸入する空気量を示す空気流量Ｇ_ａｔを算出する。

　ヒートバランス演算部２０３は、ガスタービン１０の全体のエネルギ収支、すなわち仕事の収支や熱の収支が一致することを示す計算式（以下、ガスタービンヒートバランス計算式という）に対して、計測した計測データを適用することにより、計測データが得られた際の燃料流量Ｇ_ｆｓと、空気流量Ｇ_ａｓとを算出する。ヒートバランス演算部２０３が、ガスタービンヒートバランス計算式に適用する計測データは、例えば、電力計センサ２５ｄが検出するガスタービン出力Ｇ_ｏｕｔ、排ガス圧力センサ２５ｋが検出する排ガス圧力Ｐ_２Ｔ、排ガス温度センサ２５ｌが検出する排ガス温度Ｔ_２Ｔ、圧縮機入口圧力センサ２５ｆが検出する圧縮機入口圧力Ｐ_１Ｃ、圧縮機入口温度センサ２５ｇが検出する圧縮機入口温度Ｔ_１Ｃ、圧縮機出口圧力センサ２５ｈが検出する圧縮機出口圧力Ｐ_２Ｃ、及び圧縮機出口温度センサ２５ｉが検出する圧縮機出口温度Ｔ_２Ｃである。

　状態検出部２０５は、ガスタービン制御装置２０に接続する各種のセンサ２５ａ～２５ｌが検出する各種の計測データによって表されるガスタービン発電プラント１００の状態量の時間的な変化量や、発電機５０が発電する電力などに基づいて、ガスタービン発電プラント１００の状態が、過渡状態であるか、静定状態であるかを検出する。状態検出部２０５は、例えば、状態量の時間的な変化量が、予め定める閾値未満である場合にガスタービン発電プラント１００の状態を静定状態として検出する。これに対して、状態検出部２０５は、時間的な変化量が、当該閾値以上である場合にガスタービン発電プラント１００の状態を過渡状態として検出する。

　補正処理部２０４は、燃料流量補正部２１０と、空気流量補正部２２０とを備える。燃料流量補正部２１０は、燃料側管理部２１１、燃料側記憶部２１２、燃料補正係数算出部２１３、及び燃料補正係数適用部２１４を備える。燃料側管理部２１１は、状態検出部２０５が検出する状態にしたがって、燃料流量算出部２０１が算出する燃料流量Ｇ_ｆｔと、ヒートバランス演算部２０３が算出する燃料流量Ｇ_ｆｓとを関連付けて燃料側記憶部２１２に書き込んで記憶させる。燃料側管理部２１１は、燃料流量Ｇ_ｆｔと燃料流量Ｇ_ｆｓの最新の組み合わせを燃料補正係数算出部２１３に出力する。

　燃料側記憶部２１２は、燃料流量Ｇ_ｆｔと、燃料流量Ｇ_ｆｓとの組み合わせを時系列順に記憶する。燃料補正係数算出部２１３は、最新の燃料流量Ｇ_ｆｓを、最新の燃料流量Ｇ_ｆｔで除算して燃料補正係数Ｋ_ｆ（＝Ｇ_ｆｓ／Ｇ_ｆｔ）を算出する。燃料補正係数適用部２１４は、例えば、乗算器であり、燃料流量算出部２０１が算出する燃料流量Ｇ_ｆｔに、燃料補正係数算出部２１３が算出する燃料補正係数Ｋ_ｆを乗算することにより燃料流量Ｇ_ｆｔを補正し、補正後の燃料流量Ｇ_ｆｔである補正燃料流量Ｇ_ｆｃｏｒ（＝Ｇ_ｆｔ×Ｋ_ｆ）を出力する。

　空気流量補正部２２０は、空気側管理部２２１、空気側記憶部２２２、空気補正係数算出部２２３、及び空気補正係数適用部２２４を備える。空気側管理部２２１は、状態検出部２０５が検出する状態にしたがって、空気流量算出部２０２が算出する空気流量Ｇ_ａｔと、ヒートバランス演算部２０３が算出する空気流量Ｇ_ａｓとを関連付けて空気側記憶部２２２に書き込んで記憶させる。空気側管理部２２１は、空気流量Ｇ_ａｔと空気流量Ｇ_ａｓの最新の組み合わせを空気補正係数算出部２２３に出力する。

　空気側記憶部２２２は、空気流量Ｇ_ａｔと、空気流量Ｇ_ａｓとの組み合わせを時系列順に記憶する。空気補正係数算出部２２３は、最新の空気流量Ｇ_ａｓを、最新の空気流量Ｇ_ａｔで除算して空気補正係数Ｋ_ａ（＝Ｇ_ａｓ／Ｇ_ａｔ）を算出する。空気補正係数適用部２２４は、例えば、乗算器であり、空気流量算出部２０２が算出する空気流量Ｇ_ａｔに、空気補正係数算出部２２３が算出する空気補正係数Ｋ_ａを乗算することにより空気流量Ｇ_ａｔを補正し、補正後の空気流量Ｇ_ａｔである補正空気流量Ｇ_ａｃｏｒ（＝Ｇ_ａｔ×Ｋ_ａ）を出力する。

　タービン入口温度算出部２０６は、補正燃料流量Ｇ_ｆｃｏｒと、補正空気流量Ｇ_ａｃｏｒと、燃焼器車室温度センサ２５ｊが検出する燃焼器車室温度Ｔ_ＣＳと、燃料温度センサ２５ｃ－１～２５ｃ－３が検出する燃料温度Ｔ_ｆ（１），Ｔ_ｆ（２），Ｔ_ｆ（３）の平均値である燃料温度Ｔ_ｆとを、次式（６）に示すように、予め定められる燃焼器１２周りの熱収支を解く式、すなわち、燃焼器１２に関する熱エネルギ収支に関する物理モデル式である関数ｆ_３（・）に代入してタービン入口温度Ｔ_１Ｔを算出する。

　ここで、上記の式（６）の関数ｆ_３（・）によって表される熱エネルギ収支に関する物理モデル式とは、例えば、非定常物理モデルを変形して得られる式である。非定常物理モデルを表す式とは、例えば、燃焼器１２に流入する熱エネルギと、燃焼器１２から流出する熱エネルギとが等しいことを表すモデル式であり、当該モデル式に含まれるタービン入口温度Ｔ_１Ｔを左辺とし、他の要素を右辺とした式に変形することにより熱エネルギ収支に関する物理モデル式が得られる。なお、非定常物理モデルにおける燃焼器１２に流入する熱エネルギは、燃料の熱エネルギ、空気の熱エネルギ、燃焼ガスの発熱エネルギの合計により表される。燃焼器１２から流出する熱エネルギはタービン１３の入口における熱エネルギにより表される。

　燃料配分比算出部２０７は、タービン入口温度算出部２０６が算出するタービン入口温度Ｔ_１Ｔに基づいて、予め定められるタービン入口温度Ｔ_１Ｔと燃料配分比の関係式により燃料供給系統４０－１，４０－２，４０－３の各々に対する燃料配分比Ｄ_１，Ｄ_２，Ｄ_３を算出する。燃料配分比算出部２０７が算出する３つの燃料配分比Ｄ_１，Ｄ_２，Ｄ_３の合計は１００％になる。

　弁開度算出部２０８は、燃料配分比算出部２０７が算出する燃料配分比Ｄ_１，Ｄ_２，Ｄ_３と、外部から与えられる燃料制御信号指令値（ＣＳＯ：Control Signal Output）とに基づいて、予め定められる弁開度算出式により、流調弁４１－１，４１－２，４１－３の各々に対する弁開度Ｏ_１，Ｏ_２，Ｏ_３を算出する。なお、燃料制御信号指令値（ＣＳＯ）とは、燃焼器１２に供給する燃料流量を制御する制御出力信号である。弁開度算出部２０８は、算出した弁開度Ｏ_１，Ｏ_２，Ｏ_３の各々を、各々に対応する流調弁４１－１，４１－２，４１－３に対する指令値として与える。これにより、流調弁４１－１は、弁開度Ｏ_１にしたがって開度を調整し、流調弁４１－２は、弁開度Ｏ_２にしたがって開度を調整し、流調弁４１－３は、弁開度Ｏ_３にしたがって開度を調整する。弁開度算出部２０８は、算出した弁開度Ｏ_１，Ｏ_２，Ｏ_３を燃料流量算出部２０１にフィードバックする。

（第１の実施形態のガスタービン制御装置の動作例）
　図３から図５を参照しつつガスタービン制御装置２０による処理について説明する。流調弁前圧力センサ２５ａ－１～２５ａ－３、流調弁後圧力センサ２５ｂ－１～２５ｂ－３、及び燃料温度センサ２５ｃ－１～２５ｃ－３は、検出周期ごとに計測を行い計測により検出した計測データを、制御回線を通じてガスタービン制御装置２０に送信する。

　燃料流量算出部２０１は、流調弁前圧力センサ２５ａ－１～２５ａ－３が送信する流調弁前圧力Ｐ_１ＦＶ（１），Ｐ_１ＦＶ（２），Ｐ_１ＦＶ（３）と、流調弁後圧力センサ２５ｂ－１～２５ｂ－３が送信する流調弁後圧力Ｐ_２ＦＶ（１），Ｐ_２ＦＶ（２），Ｐ_２ＦＶ（３）と、燃料温度センサ２５ｃ－１～２５ｃ－３が送信する燃料温度Ｔ_ｆ（１），Ｔ_ｆ（２），Ｔ_ｆ（３）とを受信して取り込む。燃料流量算出部２０１は、弁開度算出部２０８が前回の検出周期において算出して出力した弁開度Ｏ_１，Ｏ_２，Ｏ_３を取り込む（Ｓ１－１）。

　燃料流量算出部２０１は、ステップＳ１－１の処理において取り込んだデータを、上記したように、式（１），（２），（３）に代入して、燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｔ（２），Ｇ_ｆｔ（３）を算出する。燃料流量算出部２０１は、式（４）に示すように、燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｔ（２），Ｇ_ｆｔ（３）を合計して単位時間当たりにガスタービン１０の燃焼器１２に流入する燃料の量を示す燃料流量Ｇ_ｆｔを算出する。燃料流量算出部２０１は、算出した燃料流量Ｇ_ｆｔを燃料側管理部２１１と、燃料補正係数適用部２１４とに出力する（Ｓ２－１）。

　電力計センサ２５ｄ、排ガス圧力センサ２５ｋ、排ガス温度センサ２５ｌ、圧縮機入口圧力センサ２５ｆ、圧縮機入口温度センサ２５ｇ、圧縮機出口圧力センサ２５ｈ、及び圧縮機出口温度センサ２５ｉは、検出周期ごとに計測を行い計測により検出した計測データを、制御回線を通じてガスタービン制御装置２０に送信する。

　ヒートバランス演算部２０３は、電力計センサ２５ｄが送信するガスタービン出力Ｇ_ｏｕｔと、排ガス圧力センサ２５ｋが送信する排ガス圧力Ｐ_２Ｔと、排ガス温度センサ２５ｌが送信する排ガス温度Ｔ_２Ｔと、圧縮機入口圧力センサ２５ｆが送信する圧縮機入口圧力Ｐ_１Ｃと、圧縮機入口温度センサ２５ｇが送信する圧縮機入口温度Ｔ_１Ｃと、圧縮機出口圧力センサ２５ｈが送信する圧縮機出口圧力Ｐ_２Ｃと、圧縮機出口温度センサ２５ｉが送信する圧縮機出口温度Ｔ_２Ｃとを受信して取り込む（Ｓ１－２）。ヒートバランス演算部２０３は、ステップＳ１－２の処理において取り込んだデータと、ガスタービンヒートバランス計算式とに基づいて、燃料流量Ｇ_ｆｓと、空気流量Ｇ_ａｓとを算出する。ヒートバランス演算部２０３は、算出した燃料流量Ｇ_ｆｓを燃料側管理部２１１に出力し、算出した空気流量Ｇ_ａｓを空気側管理部２２１に出力する（Ｓ２－２）。

　インデックス差圧センサ２５ｅは、検出周期ごとに計測を行い計測により検出した計測データである圧縮機インデックス差圧Ｐ_{ｉｎｄｅｘ}のデータを、制御回線を通じてガスタービン制御装置２０に送信する。空気流量算出部２０２は、インデックス差圧センサ２５ｅが送信する圧縮機インデックス差圧Ｐ_{ｉｎｄｅｘ}と、圧縮機入口圧力センサ２５ｆが送信する圧縮機入口圧力Ｐ_１Ｃと、圧縮機入口温度センサ２５ｇが送信する圧縮機入口温度Ｔ_１Ｃとを受信して取り込む（Ｓ１－３）。空気流量算出部２０２は、ステップＳ１－３の処理において取り込んだデータを、上記したように、式（５）に代入して単位時間当たりにガスタービン１０の圧縮機１１が吸入する空気量を示す空気流量Ｇ_ａｔを算出する。空気流量算出部２０２は、算出した空気流量Ｇ_ａｔを空気側管理部２２１と、空気補正係数適用部２２４とに出力する（Ｓ２－３）。なお、図３において、ステップＳ１－１，Ｓ２－１の処理のセットと、ステップＳ１－２，Ｓ２－２の処理のセットと、ステップＳ１－３，Ｓ２－３の処理のセットとは、並列に行われる。

　補正処理部２０４は、燃料流量Ｇ_ｆｔと、燃料流量Ｇ_ｆｓと、空気流量Ｇ_ａｔと、空気流量Ｇ_ａｓとを受けると、図５に示す記憶部管理処理のサブルーチンを開始する（Ｓ３）。

　燃料流量補正部２１０において、燃料側管理部２１１は、燃料流量算出部２０１が出力する燃料流量Ｇ_ｆｔと、ヒートバランス演算部２０３が出力する燃料流量Ｇ_ｆｓとを取り込む。空気流量補正部２２０において、空気側管理部２２１は、空気流量算出部２０２が出力する空気流量Ｇ_ａｔと、ヒートバランス演算部２０３が出力する空気流量Ｇ_ａｓとを取り込む（Ｓ１００）。燃料側管理部２１１及び空気側管理部２２１の各々は、各々が行うステップＳ１００の処理を終了した時点で状態検出部２０５が検出するガスタービン発電プラント１００の状態を示すデータを状態検出部２０５から取得する（Ｓ１０１）。燃料側管理部２１１及び空気側管理部２２１の各々は、取得した状態を示すデータが、静定状態を示しているか否かを判定する（Ｓ１０２）。

　燃料側管理部２１１は、ステップＳ１０２の処理において、取得した状態を示すデータが、静定状態を示していると判定した場合（Ｓ１０２、Ｙｅｓ）、ステップＳ１００の処理で取り込んだ燃料流量Ｇ_ｆｔと、燃料流量Ｇ_ｆｓとを関連付けて燃料側記憶部２１２に書き込んで記憶させる。燃料側管理部２１１は、燃料側記憶部２１２に順次書き込む燃料流量Ｇ_ｆｔと燃料流量Ｇ_ｆｓの組み合わせが燃料側記憶部２１２において時系列順に記憶されるように、例えば、新たに書き込む組み合わせが、時系列順の先頭になるように燃料側記憶部２１２に書き込む（Ｓ１０３－１）。燃料側管理部２１１は、燃料側記憶部２１２に書き込んだ燃料流量Ｇ_ｆｔと燃料流量Ｇ_ｆｓの組み合わせを燃料補正係数算出部２１３に出力する（Ｓ１０４－１）。

　一方、燃料側管理部２１１は、ステップＳ１０２の処理において、取得した状態を示すデータが、静定状態を示していないと判定した場合（Ｓ１０２、Ｎｏ）、ステップＳ１００の処理で取り込んだ燃料流量Ｇ_ｆｔと、燃料流量Ｇ_ｆｓとを破棄し、燃料側記憶部２１２に記憶されている組み合わせのうち最新の燃料流量Ｇ_ｆｔと燃料流量Ｇ_ｆｓの組み合わせを読み出す。上記したように燃料側記憶部２１２が最新の組み合わせが先頭になるように時系列順に記憶している場合、燃料側管理部２１１は、燃料側記憶部２１２の先頭に記憶されている組み合わせを最新の燃料流量Ｇ_ｆｔと燃料流量Ｇ_ｆｓの組み合わせとして読み出す。燃料側管理部２１１は、読み出した最新の燃料流量Ｇ_ｆｔと燃料流量Ｇ_ｆｓの組み合わせを燃料補正係数算出部２１３に出力する（Ｓ１０５－１）。

　空気側管理部２２１が、ステップＳ１０２の処理において、取得した状態を示すデータが、静定状態を示していると判定した場合（Ｓ１０２、Ｙｅｓ）、ステップＳ１０３－２，Ｓ１０４－２の処理が行われる。一方、空気側管理部２２１が、ステップＳ１０２の処理において、取得した状態を示すデータが、静定状態を示していないと判定した場合（Ｓ１０２、Ｎｏ）、ステップＳ１０５－２の処理が行われる。ここで、ステップＳ１０３－２，Ｓ１０４－２，Ｓ１０５－２の処理の各々は、各々のステップ番号の枝番号を「－１」に置き換えたステップＳ１０３－１，Ｓ１０４－１，Ｓ１０５－１の処理において、燃料流量Ｇ_ｆｔを空気流量Ｇ_ａｔに読み替え、燃料流量Ｇ_ｆｓを空気流量Ｇ_ａｓに読み替え、燃料側管理部２１１を空気側管理部２２１に読み替え、燃料側記憶部２１２を空気側記憶部２２２に読み替え、燃料補正係数算出部２１３を空気補正係数算出部２２３に読み替えた処理である。

　これにより、燃料側管理部２１１は、状態検出部２０５が検出する状態が静定状態である場合、直近で取り込んだ燃料流量Ｇ_ｆｔと燃料流量Ｇ_ｆｓの組み合わせを、燃料補正係数算出部２１３に出力し、状態検出部２０５が検出する状態が静定状態でない場合、燃料側記憶部２１２に記憶されている最新の燃料流量Ｇ_ｆｔと燃料流量Ｇ_ｆｓの組み合わせを、燃料補正係数算出部２１３に出力することになる。言い換えると、燃料側管理部２１１は、状態検出部２０５が検出する状態がいずれの状態であったとしても、静定状態における最新の燃料流量Ｇ_ｆｔと燃料流量Ｇ_ｆｓの組み合わせを、燃料補正係数算出部２１３に出力することになる。同様に、空気側管理部２２１も、状態検出部２０５が検出する状態がいずれの状態であったとしても、静定状態における最新の空気流量Ｇ_ａｔと空気流量Ｇ_ａｓの組み合わせを、空気補正係数算出部２２３に出力することになる。ステップＳ１０４－１，Ｓ１０４－２，Ｓ１０５－１，Ｓ１０５－２の処理が終了すると、図５の記憶部管理処理のサブルーチンが終了して、図３の処理に戻る。なお、図５において、燃料側管理部２１１が行う処理と、空気側管理部２２１が行う処理とは、並列に行われる。

　図３のステップＳ４－１に示すように、燃料補正係数算出部２１３は、燃料側管理部２１１が出力する燃料流量Ｇ_ｆｔと燃料流量Ｇ_ｆｓの組み合わせを取り込む（Ｓ４－１）。燃料補正係数算出部２１３は、燃料流量Ｇ_ｆｓを燃料流量Ｇ_ｆｔで除算して燃料補正係数Ｋ_ｆを算出する。燃料補正係数算出部２１３は、算出した燃料補正係数Ｋ_ｆを燃料補正係数適用部２１４に出力する（Ｓ５－１）。燃料補正係数適用部２１４は、燃料流量算出部２０１が出力する燃料流量Ｇ_ｆｔと、燃料補正係数算出部２１３が出力する燃料補正係数Ｋ_ｆを取り込む。燃料補正係数適用部２１４は、燃料流量Ｇ_ｆｔに、燃料補正係数Ｋ_ｆを乗算して補正燃料流量Ｇ_ｆｃｏｒを算出する。燃料補正係数適用部２１４は、算出した補正燃料流量Ｇ_ｆｃｏｒをタービン入口温度算出部２０６に出力する（Ｓ６－１）。

　空気流量補正部２２０によって、図３のステップＳ４－２，Ｓ５－２，Ｓ６－２の処理が行われる。ここで、ステップＳ４－２，Ｓ５－２，Ｓ６－２の処理の各々は、各々のステップ番号の枝番号を「－１」に置き換えたＳ４－１，Ｓ５－１，Ｓ６－１の処理において、燃料流量Ｇ_ｆｔを空気流量Ｇ_ａｔに読み替え、燃料流量Ｇ_ｆｓを空気流量Ｇ_ａｓに読み替え、燃料補正係数算出部２１３を空気補正係数算出部２２３に読み替え、燃料補正係数適用部２１４を空気補正係数適用部２２４に読み替え、燃料流量算出部２０１を空気流量算出部２０２に読み替え、燃料補正係数Ｋ_ｆを空気補正係数Ｋ_ａに読み替え、補正燃料流量Ｇ_ｆｃｏｒを補正空気流量Ｇ_ａｃｏｒに読み替えた処理である。

　なお、図３において、ステップＳ４－１，Ｓ５－１，Ｓ６－１の処理のセットと、ステップＳ４－２，Ｓ５－２，Ｓ６－２の処理のセットとは、並列に行われる。図３に示す「Ａ」のマークは、処理が継続して行われることを示すマークであり、「Ａ」のマークが示すように、その後、図４の処理が行われる。

　図４に示すように、タービン入口温度算出部２０６は、燃料補正係数適用部２１４が出力する補正燃料流量Ｇ_ｆｃｏｒと、空気補正係数適用部２２４が出力する補正空気流量Ｇ_ａｃｏｒとを取り込む。タービン入口温度算出部２０６は、燃焼器車室温度センサ２５ｊが送信する燃焼器車室温度Ｔ_ＣＳと、燃料温度センサ２５ｃ－１～２５ｃ－３が送信する燃料温度Ｔ_ｆ（１），Ｔ_ｆ（２），Ｔ_ｆ（３）とを受信して取り込む。なお、燃焼器車室温度センサ２５ｊは、検出周期ごとに計測を行い計測により検出した計測データである燃焼器車室温度Ｔ_ＣＳのデータを、制御回線を通じてガスタービン制御装置２０に送信する。タービン入口温度算出部２０６は、取り込んだ燃料温度Ｔ_ｆ（１），Ｔ_ｆ（２），Ｔ_ｆ（３）の平均値である燃料温度Ｔ_ｆを算出する。タービン入口温度算出部２０６は、上記したように、補正燃料流量Ｇ_ｆｃｏｒと、補正空気流量Ｇ_ａｃｏｒと、燃焼器車室温度Ｔ_ＣＳと、燃料温度Ｔ_ｆとを式（６）に代入してタービン入口温度Ｔ_１Ｔを算出する。タービン入口温度算出部２０６は、算出したタービン入口温度Ｔ_１Ｔを燃料配分比算出部２０７に出力する（Ｓ２０）。

　燃料配分比算出部２０７は、タービン入口温度算出部２０６が出力するタービン入口温度Ｔ_１Ｔを取り込む。燃料配分比算出部２０７は、取り込んだタービン入口温度Ｔ_１Ｔに基づいて、タービン入口温度Ｔ_１Ｔと燃料配分比の関係式により、燃料供給系統４０－１，４０－２，４０－３の各々に対する燃料配分比Ｄ_１，Ｄ_２，Ｄ_３を算出する。燃料配分比算出部２０７は、算出した燃料配分比Ｄ_１，Ｄ_２，Ｄ_３を弁開度算出部２０８に出力する（Ｓ２１）。

　弁開度算出部２０８は、燃料配分比算出部２０７が出力する燃料配分比Ｄ_１，Ｄ_２，Ｄ_３と、外部から与えられる燃料制御信号指令値（ＣＳＯ）とを取り込む。弁開度算出部２０８は、取り込んだ燃料配分比Ｄ_１，Ｄ_２，Ｄ_３と、燃料制御信号指令値（ＣＳＯ）と基づいて、弁開度算出式により、流調弁４１－１，４１－２，４１－３の各々に対する弁開度Ｏ_１，Ｏ_２，Ｏ_３を算出する。弁開度算出部２０８は、算出した弁開度Ｏ_１，Ｏ_２，Ｏ_３の各々を、各々に対応する流調弁４１－１，４１－２，４１－３に出力する。弁開度算出部２０８は、算出した弁開度Ｏ_１，Ｏ_２，Ｏ_３を燃料流量算出部２０１に出力し（Ｓ２２）、処理が終了する。

　上記の図３及び図４に示す一連の処理が、１回の検出周期の処理として行われ、周期ごとに図３及び図４に示す一連の処理が行われることになる。

（第１の実施形態の構成における作用・効果）
　上記したように、第１の実施形態のガスタービン制御装置２０は、２通りの燃料流量Ｇ_ｆｔ，Ｇ_ｆｓを算出する。一方の燃料流量Ｇ_ｆｔは、燃料流量算出部２０１によって、弁開度Ｏ_１，Ｏ_２，Ｏ_３、流調弁前圧力Ｐ_１ＦＶ（１），Ｐ_１ＦＶ（２），Ｐ_１ＦＶ（３）と、流調弁後圧力Ｐ_２ＦＶ（１），Ｐ_２ＦＶ（２），Ｐ_２ＦＶ（３）と、燃料温度Ｔ_ｆ（１），Ｔ_ｆ（２）Ｔ_ｆ（３）というパラメータを用いて算出される。燃料流量算出部２０１が用いるパラメータは、応答が速いため、ガスタービン発電プラント１００の状態が過渡状態であっても、過渡状態における激しい変化に追随して、一定の精度を有する燃料流量Ｇ_ｆｔを算出することができる。ただし、燃料流量算出部２０１が燃料流量Ｇ_ｆｔの算出に用いる燃料流量算出関数ｆ_１（・）は、当該関数を構築する際に予め定める流調弁４１のＣｖ特性が、真の流調弁４１のＣｖ特性と相違しているといった諸元のズレによる誤差、燃料流量を算出するためにモデル化を行った際の誤差といった絶対的な誤差を含んでいる可能性がある。なお、Ｃｖ特性とは、流調弁４１においてどれだけ燃料を流すことができるのかということを示す特性である。このような絶対的な誤差を含んでいる可能性があるため、燃料流量算出関数ｆ_１（・）に代入するパラメータの精度を高くしたとしても、当該関数を用いて得られる燃料流量Ｇ_ｆｔついては、絶対的な精度が低い可能性がある。

　これに対して、他方の燃料流量Ｇ_ｆｓは、ヒートバランス演算部２０３によって、ガスタービン出力Ｇ_ｏｕｔと、排ガス圧力Ｐ_２Ｔと、排ガス温度Ｔ_２Ｔと、圧縮機入口圧力Ｐ_１Ｃと、圧縮機入口温度Ｔ_１Ｃと、圧縮機出口圧力Ｐ_２Ｃと、圧縮機出口温度Ｔ_２Ｃというパラメータを用いて算出される。ヒートバランス演算部２０３が、燃料流量Ｇ_ｆｓの算出に用いるガスタービンヒートバランス計算式は、ガスタービン１０の全体のエネルギ収支が一致することを示す計算式であるため、燃料流量算出部２０１が用いる燃料流量算出関数ｆ_１（・）に比べると、絶対的な精度の高い燃料流量Ｇ_ｆｓを算出することができる。ただし、ヒートバランス演算部２０３が用いるパラメータの中で、排ガス温度Ｔ_２Ｔなどのガスタービン１０から得られる温度の計測データは、応答遅れが大きい場合がある。また、発電機５０が発電する電力を示すガスタービン出力Ｇ_ｏｕｔは、発電機５０に接続している電力系統において需給バランスの崩れなどが発生するなどの外乱の影響を受ける場合がある。そのため、ヒートバランス演算部２０３が算出する燃料流量Ｇ_ｆｓは、ガスタービン発電プラント１００の状態が静定状態である場合には、燃料流量Ｇ_ｆｔよりも精度が高くなる一方で、過渡状態の場合には、燃料流量Ｇ_ｆｔよりも精度が低くなる可能性がある。

　補正処理部２０４の燃料流量補正部２１０は、上記の２通りの燃料流量Ｇ_ｆｔ，Ｇ_ｆｓを用いて、図３から図５を参照して説明した処理により、静定状態と過渡状態の両方の状態において、タービン入口温度算出部２０６に与える燃料流量の精度が低くならないようにするために、以下のような処理を行っているということがいえる。例えば、時刻ｔ１において、ヒートバランス演算部２０３が燃料流量Ｇ_ｆｓを算出し、燃料流量算出部２０１が燃料流量Ｇ_ｆｔを算出したとする。更に、燃料側管理部２１１が、ヒートバランス演算部２０３が出力する時刻ｔ１の燃料流量Ｇ_ｆｓを取り込んだ際に、状態検出部２０５が検出するガスタービン発電プラント１００の状態が静定状態であるとする。この場合、燃料側管理部２１１は、時刻ｔ１の燃料流量Ｇ_ｆｓ，Ｇ_ｆｔの組み合わせを燃料側記憶部２１２に記録すると共に、時刻ｔ１の燃料流量Ｇ_ｆｔ，Ｇ_ｆｓの組み合わせを燃料補正係数算出部２１３に出力する。ここで、燃料補正係数算出部２１３が算出する燃料補正係数Ｋ_ｆは、Ｋ_ｆ＝Ｇ_ｆｓ／Ｇ_ｆｔである。そのため、燃料補正係数適用部２１４は、Ｇ_ｆｔ×Ｋ_ｆ＝Ｇ_ｆｓという演算を行うことになり、燃料補正係数適用部２１４が算出する補正燃料流量Ｇ_ｆｃｏｒは、時刻ｔ１の燃料流量Ｇ_ｆｓに等しくなる。この場合の燃料流量Ｇ_ｆｓは、静定状態のものであるため、補正燃料流量Ｇ_ｆｃｏｒは、燃料流量Ｇ_ｆｔよりも高い精度になる。

　一方、燃料側管理部２１１は、ヒートバランス演算部２０３が出力する時刻ｔ１の燃料流量Ｇ_ｆｓを取り込んだ際に、状態検出部２０５が検出するガスタービン発電プラント１００の状態が静定状態でない場合、すなわち、過渡状態である場合、時刻ｔ１の燃料流量Ｇ_ｆｓ，Ｇ_ｆｔの組み合わせを燃料側記憶部２１２に記録しない。燃料側管理部２１１は、時刻ｔ１の燃料流量Ｇ_ｆｔ，Ｇ_ｆｓの組み合わせに替えて、燃料側記憶部２１２に書き込まれている最新の燃料流量Ｇ_ｆｔ，Ｇ_ｆｓの組み合わせを燃料補正係数算出部２１３に出力する。この最新の燃料流量Ｇ_ｆｔ，Ｇ_ｆｓの組み合わせが得られた時刻は、時刻ｔ１以前の時刻であって時刻ｔの直近の静定状態であった時刻である。この時刻を、時刻ｔ２とする。この場合、燃料補正係数算出部２１３が算出する燃料補正係数Ｋ_ｆは、Ｋ_ｆ＝（時刻ｔ２のＧ_ｆｓ）／（時刻ｔ２のＧ_ｆｔ）である。そのため、燃料補正係数適用部２１４は、（時刻ｔ１のＧ_ｆｔ）×Ｋ_ｆ＝（時刻ｔ１のＧ_ｆｔ）×（時刻ｔ２のＧ_ｆｓ）／（時刻ｔ２のＧ_ｆｔ）という演算を行うことになる。ここで、（時刻ｔ２のＧ_ｆｓ）／（時刻ｔ２のＧ_ｆｔ）という燃料補正係数Ｋ_ｆを算出する演算は、時刻ｔ１の直近において静定状態であった時刻である時刻ｔ２の燃料流量Ｇ_ｆｔ，Ｇ_ｆｓの組み合わせが用いられている演算である。静定状態における燃料流量Ｇ_ｆｓの精度は高いことから、燃料補正係数Ｋ_ｆの精度も高くなる。したがって、この精度の高い燃料補正係数Ｋ_ｆを、時刻ｔ１の燃料流量Ｇ_ｆｔに乗算することにより、時刻ｔ１の燃料流量Ｇ_ｆｔよりも精度の高い補正燃料流量Ｇ_ｆｃｏｒが得られることになる。

　このように、燃料側管理部２１１が、静定状態の場合のみ、燃料流量Ｇ_ｆｔ，Ｇ_ｆｓを関連付けて燃料側記憶部２１２に書き込んでおくことにより、過渡状態に遷移した場合であっても、燃料補正係数算出部２１３は、直近の静定状態において得られた精度の高い燃料流量Ｇ_ｆｓが得られた際の燃料流量Ｇ_ｆｔを用いて、精度の高い燃料補正係数Ｋ_ｆを算出することができる。また、上記したように、燃料補正係数算出部２１３は、静定状態においても精度の高い燃料補正係数Ｋ_ｆを算出することができる。したがって、燃料補正係数適用部２１４が算出する補正燃料流量Ｇ_ｆｃｏｒの精度は、静定状態と過渡状態のいずれの状態であっても、燃料流量算出部２０１が算出する燃料流量Ｇ_ｆｔよりも高くなる。

　ガスタービン制御装置２０では、空気流量についても、２通りの空気流量Ｇ_ａｔ，Ｇ_ａｓを算出するようにしている。空気流量Ｇ_ａｔ，Ｇ_ａｓの関係についても、燃料流量Ｇ_ｆｔ，Ｇ_ｆｓの関係と同様のことがいえるため、空気補正係数適用部２２４が算出する補正空気流量Ｇ_ａｃｏｒの精度は、静定状態と過渡状態のいずれの状態であっても、空気流量算出部２０２が算出する空気流量Ｇ_ａｔよりも高くなる。そのため、タービン入口温度算出部２０６は、静定状態と過渡状態のいずれの状態であっても、一意に定めることが難しいタービン効率を用いることなく、精度の高いタービン入口温度Ｔ_１Ｔを算出することが可能になる。

　ヒートバランス演算部２０３は、ガスタービンヒートバランス計算式を用いていることから、直近のガスタービン１０の性能を反映した燃料流量Ｇ_ｆｓと、空気流量Ｇ_ａｓとを算出していることになる。燃料流量補正部２１０は、直近のガスタービン１０の性能が反映された燃料流量Ｇ_ｆｓにより、燃料流量算出部２０１が算出する燃料流量Ｇ_ｆｔを補正する。同様に、空気流量補正部２２０は、直近のガスタービン１０の性能が反映された空気流量Ｇ_ａｓにより、空気流量算出部２０２が算出する空気流量Ｇ_ａｔを補正する。そのため、タービン入口温度算出部２０６は、直近のガスタービン１０の性能が反映されたタービン入口温度Ｔ_１Ｔを算出することができることになる。そうしてみると、補正処理部２０４は、経年劣化等によってガスタービン１０の性能が変化した場合であっても、その性能の変化に追随した補正を燃料流量Ｇ_ｆｔと、空気流量Ｇ_ａｔとに対して行うことを可能にしている。すなわち、補正処理部２０４による補正により、長期間にわたって精度が高く、かつ経年劣化等に対してロバスト性の高いタービン入口温度Ｔ_１Ｔを推定することが可能になる。

（燃料流量と空気流量の各々の平均値より各々の補正係数を算出する手法）
　上記の第１の実施形態のガスタービン制御装置２０では、燃料流量補正部２１０において、燃料側管理部２１１は、図５のステップＳ１０４－１の処理において、直前のステップＳ１００の処理において取り込んだ１組の燃料流量Ｇ_ｆｔと燃料流量Ｇ_ｆｓを燃料補正係数算出部２１３に出力し、ステップＳ１０５－１の処理において、燃料側記憶部２１２に記憶されている最新の１組の燃料流量Ｇ_ｆｔと燃料流量Ｇ_ｆｓを燃料補正係数算出部２１３に出力している。空気流量補正部２２０においてもステップＳ１０４－２，Ｓ１０５－２の処理において、同様に１組の空気流量Ｇ_ａｔと空気流量Ｇ_ａｓを空気補正係数算出部２２３に出力している。これに対して、例えば、以下に示す図６，図７に示すような処理を行うようにしてもよい。

　図６は、図５に示す記憶部管理処理のサブルーチンに置き換えて行われる記憶部管理処理のサブルーチンのフローチャートである。図７は、図３に示すフローチャートのステップＳ３以降の処理を示すフローチャートであり、ステップＳ４－１～Ｓ６－１及びステップＳ４－２～Ｓ６－２に替えて行われる処理を示している。図３において、ステップＳ２－１，Ｓ２－２，Ｓ２－３までの処理が行われた後、ステップＳ３の処理として、図６に示す記憶部管理処理のサブルーチンが開始され、図６のステップＳ１０２の処理までは、図５と同一の処理が行われる。

　燃料側管理部２１１は、ステップＳ１０２の処理において、取得した状態を示すデータが、静定状態を示していると判定した場合（Ｓ１０２、Ｙｅｓ）、図５のステップＳ１０３－１と同一の処理、すなわち、ステップＳ１００の処理で取り込んだ燃料流量Ｇ_ｆｔと、燃料流量Ｇ_ｆｓとを関連付けて燃料側記憶部２１２に書き込んで記憶させる処理を行う（Ｓ１０３－１）。燃料側管理部２１１は、燃料側記憶部２１２に記憶されている燃料流量Ｇ_ｆｔと燃料流量Ｇ_ｆｓの組み合わせの数が、予め定められる所定数以上であるか否かを判定する（Ｓ１０５－１）。燃料側管理部２１１は、燃料側記憶部２１２に記憶されている燃料流量Ｇ_ｆｔと燃料流量Ｇ_ｆｓの組み合わせの数が、予め定められる所定数以上であると判定した場合（Ｓ１０５－１、Ｙｅｓ）、燃料側記憶部２１２に記憶されている全ての燃料流量Ｇ_ｆｔと燃料流量Ｇ_ｆｓの組み合わせを読み出し、読み出した燃料流量Ｇ_ｆｔと燃料流量Ｇ_ｆｓの組み合わせを含む更新指示信号を燃料補正係数算出部２１３に出力する（Ｓ１０６－１）。一方、燃料側管理部２１１は、燃料側記憶部２１２に記憶されている燃料流量Ｇ_ｆｔと燃料流量Ｇ_ｆｓの組み合わせの数が、予め定められる所定数以上ではないと判定した場合（Ｓ１０５－１、Ｙｅｓ）、処理をステップＳ１０７に進める。

　空気流量補正部２２０によって、図６のステップＳ１０３－２，Ｓ１０５－２，Ｓ１０６－２の処理が行われる。ここで、ステップＳ１０３－２，Ｓ１０５－２，Ｓ１０６－２の処理の各々は、各々のステップ番号の枝番号を「－１」に置き換えたＳ１０３－１，Ｓ１０５－１，Ｓ１０６－１の処理において、燃料流量Ｇ_ｆｔを空気流量Ｇ_ａｔに読み替え、燃料流量Ｇ_ｆｓを空気流量Ｇ_ａｓに読み替え、燃料側管理部２１１を空気側管理部２２１に読み替え、燃料側記憶部２１２を空気側記憶部２２２に読み替え、燃料補正係数算出部２１３を空気補正係数算出部２２３に読み替えた処理である。

　燃料側管理部２１１と空気側管理部２２１の各々は、ステップＳ１０２の処理において、取得した状態を示すデータが、静定状態を示していないと判定した場合（Ｓ１０２、Ｎｏ）、各々がステップＳ１０７の処理を行う。ここで、ステップＳ１０７の処理は、燃料側管理部２１１によって行われる場合、燃料側管理部２１１が非更新指示信号を燃料補正係数算出部２１３に出力する処理であり、空気側管理部２２１によって行われる場合、空気側管理部２２１が非更新指示信号を空気補正係数算出部２２３に出力する処理である（Ｓ１０７）。

　ステップＳ１０６－１，Ｓ１０６－２，Ｓ１０７の処理が終了すると、記憶部管理処理のサブルーチンが終了して、図７に示すステップＳ１０－１，Ｓ１０－２以降の処理が開始される。なお、図６において、燃料側管理部２１１が行う処理と、空気側管理部２２１が行う処理とは、並列に行われる。

　燃料補正係数算出部２１３は、燃料側管理部２１１が出力する指示信号の種類を判定する（Ｓ１０－１）。燃料補正係数算出部２１３は、燃料側管理部２１１が出力する指示信号が、更新指示信号であると判定した場合（Ｓ１０－１，更新指示信号）、更新指示信号に含まれる複数の燃料流量Ｇ_ｆｔと燃料流量Ｇ_ｆｓの組み合わせを読み出して取り込む（Ｓ１１－１）。燃料補正係数算出部２１３は、複数の燃料流量Ｇ_ｆｔの平均値と、複数の燃料流量Ｇ_ｆｓの平均値とを算出する（Ｓ１２－１）。燃料補正係数算出部２１３は、算出した燃料流量Ｇ_ｆｓの平均値を、算出した燃料流量Ｇ_ｆｔの平均値で除算して燃料補正係数Ｋ_ｆを算出する。燃料補正係数算出部２１３は、算出した燃料補正係数Ｋ_ｆを内部の記憶領域に書き込んで記憶させると共に、算出した燃料補正係数Ｋ_ｆを燃料補正係数適用部２１４に出力する（Ｓ１３－１）。その後、ステップＳ１７－１の処理として、図３のステップＳ６－１の処理と同一の処理が行われる（Ｓ１７－１）。

　一方、燃料補正係数算出部２１３は、ステップＳ１０－１の処理において、燃料側管理部２１１が出力する指示信号が、非更新指示信号であると判定した場合（Ｓ１０－１，非更新指示信号）、内部の記憶領域に前回の燃料補正係数Ｋ_ｆが記憶されているか否かを判定する（Ｓ１４－１）。燃料補正係数算出部２１３は、内部の記憶領域に前回の燃料補正係数Ｋ_ｆが記憶されていると判定した場合（Ｓ１４－１、Ｙｅｓ）、内部の記憶領域から前回の燃料補正係数Ｋ_ｆを読み出し、読み出した燃料補正係数Ｋ_ｆを燃料補正係数適用部２１４に出力する（Ｓ１５－１）。その後、処理は、ステップＳ１７－１に進められる。燃料補正係数算出部２１３は、内部の記憶領域に前回の燃料補正係数Ｋ_ｆが記憶されていないと判定した場合（Ｓ１４－１、Ｎｏ）、Ｋ_ｆ＝１である燃料補正係数Ｋ_ｆを燃料補正係数適用部２１４に出力する（Ｓ１６－１）。その後、処理は、ステップＳ１７－１に進められる。なお、燃料補正係数Ｋ_ｆの値が「１」である場合、ステップＳ１７－１の処理において、燃料補正係数適用部２１４は、燃料流量Ｇ_ｆｔに対して「１」を乗算することになるため、実質的には、補正が行われないことになり、燃料補正係数適用部２１４は、燃料流量算出部２０１が出力する燃料流量Ｇ_ｆｔを補正燃料流量Ｇ_ｆｃｏｒとしてタービン入口温度算出部２０６に出力することになる。

　空気流量補正部２２０によって、図７のステップＳ１０－２～Ｓ１７－２の処理が行われる。ここで、ステップＳ１０－２～Ｓ１７－２の処理の各々は、各々のステップ番号の枝番号を「－１」に置き換えたＳ１０－１～Ｓ１７－１の処理において、燃料流量Ｇ_ｆｔを空気流量Ｇ_ａｔに読み替え、燃料流量Ｇ_ｆｓを空気流量Ｇ_ａｓに読み替え、燃料側管理部２１１を空気側管理部２２１に読み替え、燃料補正係数算出部２１３を空気補正係数算出部２２３に読み替え、燃料補正係数適用部２１４を空気補正係数適用部２２４に読み替え、燃料流量算出部２０１を空気流量算出部２０２に読み替え、燃料補正係数Ｋ_ｆを空気補正係数Ｋ_ａに読み替え、補正燃料流量Ｇ_ｆｃｏｒを補正空気流量Ｇ_ａｃｏｒに読み替えた処理である。

　なお、図７において、ステップＳ１０－１～Ｓ１７－１の処理のセットと、ステップＳ１０－２～Ｓ１７－２の処理のセットとは、並列に行われる。図７に示す「Ａ」のマークは、処理が継続して行われることを示すマークであり、「Ａ」のマークが示すように、その後、図４の処理が行われる。

　上記の図７の処理によって得られる燃料補正係数Ｋ_ｆは、複数の燃料流量Ｇ_ｆｔの平均値と、複数の燃料流量Ｇ_ｆｓの平均値とによって算出された値であることから、過去の静定状態における燃料流量Ｇ_ｆｔ，Ｇ_ｆｓの関係が反映された値になり、空気補正係数Ｋ_ａも、同様に、過去の静定状態における空気流量Ｇ_ａｔ，Ｇ_ａｓの関係が反映された値になる。そのため、上記した第１の実施形態のガスタービン制御装置２０のように１組の燃料流量Ｇ_ｆｔ，Ｇ_ｆｓと、１組の空気流量Ｇ_ａｔ，Ｇ_ａｓとを用いる場合よりも、信頼性の面で精度の高い燃料補正係数Ｋ_ｆ及び空気補正係数Ｋ_ａを得ることができる。精度の高い燃料補正係数Ｋ_ｆ及び空気補正係数Ｋ_ａを用いて補正を行うことにより、補正によって得られる補正燃料流量Ｇ_ｆｃｏｒ及び補正空気流量Ｇ_ａｃｏｒの精度を高めることが可能になる。したがって、精度の高い補正燃料流量Ｇ_ｆｃｏｒ及び補正空気流量Ｇ_ａｃｏｒに基づいてタービン入口温度算出部２０６が算出するタービン入口温度Ｔ_１Ｔの精度も向上する。

　図６のステップＳ１０５－１，１０５－２の処理において用いられる所定数として、信頼性を高める程度の組み合わせの個数を示す値が予め定められることになる。図７のステップＳ１６－１，Ｓ１６－２の処理、すなわち燃料補正係数Ｋ_ｆ＝１として出力する処理、及び空気補正係数Ｋ_ａ＝１として出力する処理は、燃料側記憶部２１２及び空気側記憶部２２２に記憶されている組み合わせの個数が、所定値未満である間に行われる処理であり、燃料流量Ｇ_ｆｔ，Ｇ_ｆｓ及び空気流量Ｇ_ａｔ，Ｇ_ａｓの各々の平均値を算出することができない間、燃料流量Ｇ_ｆｔと空気流量Ｇ_ａｔの各々に対して補正を行わないことを意図する処理である。図７のステップＳ１５－１，Ｓ１５－２の処理、すなわち前回の燃料補正係数Ｋ_ｆを補正に用いる処理、及び前回の空気補正係数Ｋ_ａを補正に用いる処理は、状態検出部２０５が検出した状態が、静定状態でない場合、すなわち過渡状態の場合に行われる処理であり、過渡状態の場合には、燃料補正係数Ｋ_ｆ及び空気補正係数Ｋ_ａの更新を行わないことを意図する処理である。

　なお、図６のステップＳ１０６－１，Ｓ１０６－２の処理として、以下のような処理を行うようにしてもよい。燃料側管理部２１１は、図６のステップＳ１０６－１の処理において、上記のように、燃料側記憶部２１２に記憶されている全ての燃料流量Ｇ_ｆｔ，Ｇ_ｆｓの組み合わせを読み出すのではなく、最新の組み合わせを先頭として時系列順において新しい組み合わせが順に含まれるように予め定められる所定数の組み合わせを読み出すようにしてもよい。同様に、空気側管理部２２１は、図６のステップＳ１０６－２の処理において、空気側記憶部２２２に記憶されている全ての空気流量Ｇ_ａｔ，Ｇ_ａｓの組み合わせを読み出すのではなく、最新の組み合わせを先頭として時系列順において新しい組み合わせが順に含まれるように予め定められる所定数の組み合わせを読み出すようにしてもよい。このようにすることで、移動平均値に基づく燃料補正係数Ｋ_ｆ及び空気補正係数Ｋ_ａを得ることができるので、燃料側記憶部２１２と空気側記憶部２２２に記憶されている全ての組み合わせを用いる場合よりも計算量を削減することができ、また、より直近のガスタービン１０の性能が反映された精度の高いタービン入口温度Ｔ_１Ｔを得ることが可能になる。

　また、燃料側記憶部２１２、及び空気側記憶部２２２に、所定数以上の組み合わせが記憶される段階になると、時系列順に記憶される個々の組み合わせに含まれる値の変化が小さくなっていると想定される。このような想定を前提として、燃料側管理部２１１は、図６のステップＳ１０６－１の処理において、最新の１組の燃料流量Ｇ_ｆｔ，Ｇ_ｆｓの組み合わせを読み出し、図７においてステップＳ１２－１の処理をスキップして、ステップＳ１３－１の処理において、燃料補正係数算出部２１３が、最新の燃料流量Ｇ_ｆｓを、最新の燃料流量Ｇ_ｆｔで除算して燃料補正係数Ｋ_ｆ（＝Ｇ_ｆｓ／Ｇ_ｆｔ）を算出するようにしてもよい。同様に、空気側管理部２２１は、図６のステップＳ１０６－２の処理において、最新の１組の空気流量Ｇ_ａｔ，Ｇ_ａｓの組み合わせを読み出し、図７においてステップＳ１２－２の処理をスキップして、ステップＳ１３－２の処理において、空気補正係数算出部２２３が、最新の空気流量Ｇ_ａｓを、最新の空気流量Ｇ_ａｔで除算して空気補正係数Ｋ_ａ（＝Ｇ_ａｓ／Ｇ_ａｔ）を算出するようにしてもよい。

（第１の実施形態におけるその他の構成例）
　上記の第１の実施形態では、補正処理部２０４は、燃料流量補正部２１０と、空気流量補正部２２０とを備えている。これに対して、補正処理部２０４は、燃料流量補正部２１０のみを備えるようにしてもよい。この場合、ヒートバランス演算部２０３は、燃料流量Ｇ_ｆｓのみを算出して燃料側管理部２１１に出力し、空気流量算出部２０２が算出する空気流量Ｇ_ａｔが、そのままタービン入口温度算出部２０６に与えられることになる。また、補正処理部２０４は、空気流量補正部２２０のみを備えるようにしてもよい。この場合、ヒートバランス演算部２０３は、空気流量Ｇ_ａｓのみを算出して空気側管理部２２１に出力し、燃料流量算出部２０１が算出する燃料流量Ｇ_ｆｔが、そのままタービン入口温度算出部２０６に与えられることになる。

＜第２の実施形態＞
　第１の実施形態では、ヒートバランス演算部２０３が算出する燃料流量Ｇ_ｆｓと、空気流量Ｇ_ａｓとを用いて、燃料流量算出部２０１が算出する燃料流量Ｇ_ｆｔと、空気流量算出部２０２が算出する空気流量Ｇ_ａｔとを補正するようにしていた。ところで、燃料補正係数Ｋ_ｆ及び空気補正係数Ｋ_ａは、ヒートバランス演算部２０３において参照している計測データのみに依存するわけではなく、ガスタービン発電プラント１００における様々な状態量に依存して変化するものと考えられる。第２の実施形態では、以下に説明するように、燃料補正係数Ｋ_ｆ及び空気補正係数Ｋ_ａを、このような状態量の変化に応じて可変にする構成を備える。

（第２の実施形態のガスタービン発電プラントの構成例）
　図８は、第２の実施形態のガスタービン発電プラント１００ａの全体構成図である。第１の実施形態に係るガスタービン発電プラント１００と、第２の実施形態に係るガスタービン発電プラント１００ａとの違いは、第２の実施形態に係るガスタービン発電プラント１００ａは、ガスタービン制御装置２０に替えてガスタービン制御装置２０ａを備えており、更に、燃焼器車室圧力センサ２５ｍを備える点である。燃焼器車室圧力センサ２５ｍは、点線の矢印で示すように通信回線などの制御回線によりガスタービン制御装置２０ａに接続されており、検出周期ごとに燃焼器車室１６内の圧力である燃焼器車室圧力Ｐ_ｏｕｔを計測により検出し、検出した燃焼器車室圧力Ｐ_ｏｕｔを、制御回線を通じてガスタービン制御装置２０ａに送信する。

（第２の実施形態のガスタービン制御装置の構成例）
　ガスタービン制御装置２０ａは、図９に示すように、燃料流量算出部２０１、空気流量算出部２０２、ヒートバランス演算部２０３、補正処理部２０４ａ、状態検出部２０５、タービン入口温度算出部２０６、燃料配分比算出部２０７、及び弁開度算出部２０８を備える。補正処理部２０４ａは、燃料流量補正部２１０ａと、空気流量補正部２２０ａとを備える。燃料流量補正部２１０ａは、燃料側管理部２１１ａ、燃料側記憶部２１２ａ、燃料補正係数算出部２１３ａ、燃料補正係数適用部２１４、及び燃料側回帰分析部２１５を備える。

　燃料流量補正部２１０ａにおいて、燃料側管理部２１１ａは、状態検出部２０５が検出する状態にしたがって、燃料流量算出部２０１が算出する燃料流量Ｇ_ｆｔと、ヒートバランス演算部２０３が算出する燃料流量Ｇ_ｆｓと、燃焼器車室温度センサ２５ｊが検出する燃焼器車室温度Ｔ_ＣＳと、燃焼器車室圧力センサ２５ｍが検出する燃焼器車室圧力Ｐ_ｏｕｔとを関連付けて燃料側記憶部２１２ａに書き込んで記憶させる。燃料側管理部２１１ａは、予め定められる所定数の燃料流量Ｇ_ｆｔ、燃料流量Ｇ_ｆｓ、燃焼器車室温度Ｔ_ＣＳ、及び燃焼器車室圧力Ｐ_ｏｕｔの組み合わせを燃料側回帰分析部２１５に出力する。

　燃料側記憶部２１２ａは、燃料流量Ｇ_ｆｔ、燃料流量Ｇ_ｆｓ、燃焼器車室温度Ｔ_ＣＳ、及び燃焼器車室圧力Ｐ_ｏｕｔの組み合わせを時系列順に記憶する。燃料側回帰分析部２１５は、所定数の燃料流量Ｇ_ｆｔ、燃料流量Ｇ_ｆｓ、燃焼器車室温度Ｔ_ＣＳ、及び燃焼器車室圧力Ｐ_ｏｕｔの組み合わせに基づいて、説明変数を燃焼器車室温度Ｔ_ＣＳ及び燃焼器車室圧力Ｐ_ｏｕｔとし、目的変数を燃料補正係数Ｋ_ｆ（＝燃料流量Ｇ_ｆｓ／燃料流量Ｇ_ｆｔ）とする所定の回帰分析を行い、当該所定の回帰分析に対応する回帰式の係数を算出する。燃料補正係数算出部２１３ａは、燃料側回帰分析部２１５が算出する回帰式の係数を適用した回帰式に、燃焼器車室温度Ｔ_ＣＳと、燃焼器車室圧力Ｐ_ｏｕｔとを代入して燃料補正係数Ｋ_ｆを算出する。

　空気流量補正部２２０ａは、空気側管理部２２１ａ、空気側記憶部２２２ａ、空気補正係数算出部２２３ａ、空気補正係数適用部２２４、及び空気側回帰分析部２２５を備える。空気側管理部２２１ａは、状態検出部２０５が検出する状態にしたがって、空気流量算出部２０２が算出する空気流量Ｇ_ａｔと、ヒートバランス演算部２０３が算出する空気流量Ｇ_ａｓと、圧縮機入口温度センサ２５ｇが検出する圧縮機入口温度Ｔ_１Ｃと、燃焼器車室圧力センサ２５ｍが検出する燃焼器車室圧力Ｐ_ｏｕｔとを関連付けて空気側記憶部２２２ａに書き込んで記憶させる。空気側管理部２２１ａは、予め定められる所定数の空気流量Ｇ_ａｔ、空気流量Ｇ_ａｓ、圧縮機入口温度Ｔ_１Ｃ、及び燃焼器車室圧力Ｐ_ｏｕｔの組み合わせを空気側回帰分析部２２５に出力する。

　空気側記憶部２２２ａは、空気流量Ｇ_ａｔ、空気流量Ｇ_ａｓ、圧縮機入口温度Ｔ_１Ｃ、及び燃焼器車室圧力Ｐ_ｏｕｔの組み合わせを時系列順に記憶する。空気側回帰分析部２２５は、所定数の空気流量Ｇ_ａｔ、空気流量Ｇ_ａｓ、圧縮機入口温度Ｔ_１Ｃ、及び燃焼器車室圧力Ｐ_ｏｕｔの組み合わせに基づいて、説明変数を圧縮機入口温度Ｔ_１Ｃ及び燃焼器車室圧力Ｐ_ｏｕｔとし、目的変数を空気補正係数Ｋ_ａ（＝空気流量Ｇ_ａｓ／空気流量Ｇ_ａｔ）とする所定の回帰分析を行い、当該所定の回帰分析に対応する回帰式の係数を算出する。空気補正係数算出部２２３ａは、空気側回帰分析部２２５が算出する回帰式の係数を適用した回帰式に、圧縮機入口温度Ｔ_１Ｃと、燃焼器車室圧力Ｐ_ｏｕｔとを代入して空気補正係数Ｋ_ａを算出する。

　ここで、燃料側回帰分析部２１５及び空気側回帰分析部２２５が行う所定の回帰分析とは、例えば、線形重回帰分析である。また、上記の所定数として、例えば、所定の回帰分析において精度の良い係数を算出するのに必要な組み合わせの個数を示す値が予め定められる。

（第２の実施形態のガスタービン制御装置の動作例）
　図１０から図１２を参照しつつガスタービン制御装置２０ａによる処理について説明する。図１０に示す処理と、図１０に示す処理の中で記憶部管理処理のサブルーチンとして実行される図１１に示す処理は、燃料側管理部２１１ａ、燃料側回帰分析部２１５、空気側管理部２２１ａ、空気側回帰分析部２２５によって行われる処理であり、以下、「回帰分析処理」という。図１２に示す処理は、燃料補正係数算出部２１３ａ、燃料補正係数適用部２１４、空気補正係数算出部２２３ａ、空気補正係数適用部２２４によって行われる処理であり、以下、「補正処理」という。

（回帰分析処理について）
　図１０において、ステップＳａ１－１，Ｓａ２－１は、図３のステップＳ１－１，Ｓ２－１と同一の処理が行われる。ステップＳａ１－２，Ｓａ２－２は、図３のステップＳ１－２，Ｓ２－２と同一の処理が行われる。ステップＳａ１－３，Ｓａ２－３は、図３のステップＳ１－３，Ｓ２－３と同一の処理が行われる。なお、ステップＳａ１－１，Ｓａ２－１の処理のセットと、ステップＳａ１－２，Ｓａ２－２の処理のセットと、ステップＳａ１－３，Ｓａ２－３の処理のセットとは、並列に行われる。ステップＳａ２－１，Ｓａ２－２，Ｓａ２－３の処理の後、ステップＳａ３の処理として、図１１に示す記憶部管理処理のサブルーチンが行われる。

　燃料側管理部２１１ａは、燃焼器車室温度センサ２５ｊが送信する燃焼器車室温度Ｔ_ＣＳと、燃焼器車室圧力センサ２５ｍが送信する燃焼器車室圧力Ｐ_ｏｕｔとを受信して取り込む。燃料側管理部２１１ａは、燃料流量算出部２０１が出力する燃料流量Ｇ_ｆｔと、ヒートバランス演算部２０３が出力する燃料流量Ｇ_ｆｓとを取り込む。空気側管理部２２１ａは、圧縮機入口温度センサ２５ｇが送信する圧縮機入口温度Ｔ_１Ｃと、燃焼器車室圧力センサ２５ｍが送信する燃焼器車室圧力Ｐ_ｏｕｔとを受信して取り込む。空気流量算出部２０２が出力する空気流量Ｇ_ａｔと、ヒートバランス演算部２０３が出力する空気流量Ｇ_ａｓとを取り込む（Ｓａ１００）。

　燃料側管理部２１１ａ及び空気側管理部２２１ａの各々は、各々が行うステップＳａ１００の処理を終了した時点で状態検出部２０５が検出するガスタービン発電プラント１００ａの状態を示すデータを状態検出部２０５から取得する（Ｓａ１０１）。燃料側管理部２１１ａ及び空気側管理部２２１ａは、取得した状態を示すデータが、静定状態を示しているか否かを判定する（Ｓａ１０２）。

　燃料側管理部２１１ａは、取得した状態を示すデータが、静定状態を示していると判定した場合（Ｓａ１０２、Ｙｅｓ）、ステップＳａ１００の処理で取り込んだ燃料流量Ｇ_ｆｔ、燃料流量Ｇ_ｆｓ、燃焼器車室温度Ｔ_ＣＳ、及び燃焼器車室圧力Ｐ_ｏｕｔを関連付けて燃料側記憶部２１２ａに書き込んで記憶させる。燃料側管理部２１１ａは、燃料側記憶部２１２ａに順次書き込む燃料流量Ｇ_ｆｔ、燃料流量Ｇ_ｆｓ、燃焼器車室温度Ｔ_ＣＳ、及び燃焼器車室圧力Ｐ_ｏｕｔの組み合わせが燃料側記憶部２１２ａにおいて時系列順に記憶されるように、例えば、新たに書き込む組み合わせが、時系列順の先頭になるように燃料側記憶部２１２ａに書き込む（Ｓａ１０３－１）。

　燃料側管理部２１１ａは、燃料側記憶部２１２ａに記憶されている組み合わせの数が、所定数以上であるか否かを判定する（Ｓａ１０４－１）。燃料側管理部２１１ａが、燃料側記憶部２１２ａに記憶されている組み合わせの数が、所定数以上でないと判定したとする（Ｓａ１０４－１、Ｎｏ）。この場合、燃料側回帰分析部２１５が、所定の回帰分析を行うのに十分なデータが存在せず、図１０のステップＳａ４－１の処理に進められないことから、処理を終了する。

　一方、燃料側管理部２１１ａは、燃料側記憶部２１２ａに記憶されている組み合わせの数が、所定数以上であると判定した場合（Ｓａ１０４－１、Ｙｅｓ）、最新の組み合わせを先頭として時系列順に所定数の燃料流量Ｇ_ｆｔ、燃料流量Ｇ_ｆｓ、燃焼器車室温度Ｔ_ＣＳ、及び燃焼器車室圧力Ｐ_ｏｕｔの組み合わせを読み出し、読み出した所定数の組み合わせを燃料側回帰分析部２１５に出力して（Ｓａ１０５－１）、記憶部管理処理のサブルーチンの処理を終了する。

　空気側管理部２２１が、ステップＳａ１０２の処理において、取得した状態を示すデータが、静定状態を示していると判定した場合（Ｓａ１０２、Ｙｅｓ）、ステップＳａ１０３－２，Ｓａ１０４－２，Ｓａ１０５－２の処理が行われる。ここで、ステップＳａ１０３－２，Ｓａ１０４－２，Ｓａ１０５－２の処理の各々は、各々のステップ番号の枝番号を「－１」に置き換えたステップＳａ１０３－１，Ｓａ１０４－１，Ｓａ１０５－１の処理において、燃料流量Ｇ_ｆｔを空気流量Ｇ_ａｔに読み替え、燃料流量Ｇ_ｆｓを空気流量Ｇ_ａｓに読み替え、燃焼器車室温度Ｔ_ＣＳを圧縮機入口温度Ｔ_１Ｃに読み替え、燃料側管理部２１１ａを空気側管理部２２１ａに読み替え、燃料側記憶部２１２ａを空気側記憶部２２２ａに読み替え、燃料側回帰分析部２１５を空気側回帰分析部２２５に読み替えた処理である。

　一方、ステップＳａ１０２の処理において、燃料側管理部２１１ａ及び空気側管理部２２１ａが、取得した状態を示すデータが、静定状態を示していないと判定したとする（Ｓａ１０２、Ｎｏ）。この場合、燃料側管理部２１１ａ及び空気側管理部２２１ａの各々は、ステップＳａ１００の処理で取り込んだデータが、過渡状態において得られたデータであるため、各々に対応する燃料側記憶部２１２ａ及び空気側記憶部２２２ａに記録する必要がなく、また、各々に対応する燃料側記憶部２１２ａ及び空気側記憶部２２２ａが記憶するデータの更新がないため、回帰式の係数の更新を行う必要もないことから、処理を終了する。

　ステップＳａ１０５－１，Ｓａ１０５－２の処理が終了すると、記憶部管理処理のサブルーチンが終了して、図１０の処理に戻る。なお、図１１において、燃料側管理部２１１ａが行う処理と、空気側管理部２２１ａが行う処理とは、並列に行われる。

　図１０に戻り、燃料側回帰分析部２１５は、燃料側管理部２１１ａが出力する所定数の燃料流量Ｇ_ｆｔ、燃料流量Ｇ_ｆｓ、燃焼器車室温度Ｔ_ＣＳ、及び燃焼器車室圧力Ｐ_ｏｕｔの組み合わせを取り込む（Ｓａ４－１）。燃料側回帰分析部２１５は、取り込んだ所定数の燃料流量Ｇ_ｆｔ、燃料流量Ｇ_ｆｓ、燃焼器車室温度Ｔ_ＣＳ、及び燃焼器車室圧力Ｐ_ｏｕｔの組み合わせに基づいて、所定の回帰分析を行い、当該所定の回帰分析に対応する回帰式の係数を算出する。燃料側回帰分析部２１５は、算出した係数を燃料補正係数算出部２１３ａに出力して（Ｓａ５－１）、処理を終了する。

　同様に、空気側回帰分析部２２５は、ステップＳａ４－２，Ｓａ５－２の処理を行う。ここで、ステップＳａ４－２，Ｓａ５－２の処理の各々は、各々のステップ番号の枝番号を「－１」に置き換えたステップＳａ４－１，Ｓａ５－１の処理において、燃料流量Ｇ_ｆｔを空気流量Ｇ_ａｔに読み替え、燃料流量Ｇ_ｆｓを空気流量Ｇ_ａｓに読み替え、燃焼器車室温度Ｔ_ＣＳを圧縮機入口温度Ｔ_１Ｃに読み替え、燃料側管理部２１１ａを空気側管理部２２１ａに読み替え、燃料側回帰分析部２１５を空気側回帰分析部２２５に読み替え、燃料補正係数算出部２１３ａを空気補正係数算出部２２３ａに読み替えた処理である。

　なお、図１０において、ステップＳａ４－１，Ｓａ５－１の処理のセットと、ステップＳａ４－２，Ｓａ５－２の処理のセットとは、並列に行われる。上記の図１０に示す処理が、１回の検出周期の処理として行われ、周期ごとに図１０に示す処理が行われることになる。

（補正処理について）
　図１２に示すように、燃料補正係数算出部２１３ａは、燃焼器車室温度センサ２５ｊが送信する燃焼器車室温度Ｔ_ＣＳと、燃焼器車室圧力センサ２５ｍが送信する燃焼器車室圧力Ｐ_ｏｕｔとを受信して取り込む（Ｓａ１０－１）。燃料補正係数算出部２１３ａは、燃料側回帰分析部２１５から回帰式の係数を受けているか否かを判定する（Ｓａ１１－１）。燃料補正係数算出部２１３ａは、回帰式の係数を受けていると判定した場合（Ｓａ１１－１、Ｙｅｓ）、受けた係数を取り込んで内部の記憶領域に書き込んで記憶させる。燃料補正係数算出部２１３ａは、取り込んだ係数を適用した回帰式に、燃焼器車室温度Ｔ_ＣＳと、燃焼器車室圧力Ｐ_ｏｕｔとを代入して燃料補正係数Ｋ_ｆを算出する（Ｓａ１２－１）。

　一方、燃料補正係数算出部２１３ａは、ステップＳａ１１－１の処理において、回帰式の係数を受けていないと判定した場合（Ｓａ１１－１、Ｎｏ）、内部の記憶領域に前回の係数が記憶されているか否かを判定する（Ｓａ１３－１）。燃料補正係数算出部２１３は、内部の記憶領域に前回の係数が記憶されていると判定した場合（Ｓａ１３－１、Ｙｅｓ）、内部の記憶領域から前回の係数を読み出し、読み出した係数を適用した回帰式に、燃焼器車室温度Ｔ_ＣＳと、燃焼器車室圧力Ｐ_ｏｕｔとを代入して燃料補正係数Ｋ_ｆを算出する（Ｓａ１４－１）。燃料補正係数算出部２１３は、内部の記憶領域に前回の係数が記憶されていないと判定した場合（Ｓａ１３－１、Ｎｏ）、Ｋ_ｆ＝１である燃料補正係数Ｋ_ｆを燃料補正係数適用部２１４に出力する（Ｓａ１５－１）。ステップＳａ１２－１，Ｓａ１４－１，Ｓａ１５－１の処理の後、ステップＳａ１６－１の処理として、図３のステップＳ６－１の処理と同一の処理が行われる（Ｓａ１６－１）。

　空気補正係数算出部２２３ａは、圧縮機入口温度センサ２５ｇが送信する圧縮機入口温度Ｔ_１Ｃと、燃焼器車室圧力センサ２５ｍが送信する燃焼器車室圧力Ｐ_ｏｕｔとを受信して取り込む（Ｓａ１０－２）。その後、空気補正係数算出部２２３ａ及び空気補正係数適用部２２４によってステップＳａ１１－２～Ｓａ１６－２の処理が行われる。ここで、ステップＳａ１１－２～Ｓａ１６－２の処理の各々は、各々のステップ番号の枝番号を「－１」に置き換えたステップＳａ１１－１～Ｓａ１６－１の処理において、燃料流量Ｇ_ｆｔを空気流量Ｇ_ａｔに読み替え、燃料流量Ｇ_ｆｓを空気流量Ｇ_ａｓに読み替え、燃焼器車室温度Ｔ_ＣＳを圧縮機入口温度Ｔ_１Ｃに読み替え、燃料補正係数算出部２１３ａを空気補正係数算出部２２３ａに読み替え、燃料補正係数適用部２１４を空気補正係数適用部２２４に読み替え、燃料補正係数Ｋ_ｆを空気補正係数Ｋ_ａに読み替え、燃料流量算出部２０１を空気流量算出部２０２に読み替え、補正燃料流量Ｇ_ｆｃｏｒを補正空気流量Ｇ_ａｃｏｒに読み替えた処理である。

　なお、ステップＳａ１０－１～Ｓａ１６－１の処理のセットと、Ｓａ１０－２～Ｓａ１６－２の処理のセットとは、並列に行われる。図１２に示す「Ａ」のマークは、処理が継続して行われることを示すマークであり、「Ａ」のマークが示すように、その後、第１の実施形態の図４の処理と同一の処理が行われる。上記の図１２及び図４に示す一連の処理が、１回の検出周期の処理として行われ、周期ごとに図１２及び図４に示す一連の処理が行われることになる。

（第２の実施形態の構成における作用・効果）
　上記したように、第２の実施形態の燃料流量補正部２１０ａは、状態検出部２０５が検出する状態が、静定状態である場合、燃料流量Ｇ_ｆｔ、燃料流量Ｇ_ｆｓ、燃焼器車室温度Ｔ_ＣＳ、及び燃焼器車室圧力Ｐ_ｏｕｔの組み合わせを燃料側記憶部２１２ａに記録し、燃料側記憶部２１２ａに記録されている所定数の組み合わせを用いて燃料補正係数Ｋ_ｆを更新する。一方、燃料流量補正部２１０ａは、状態検出部２０５が検出する状態が、静定状態でない場合、燃料側記憶部２１２ａへの記録を行わず、燃料補正係数Ｋ_ｆの更新も行わない。同様に、空気流量補正部２２０ａは、状態検出部２０５が検出する状態が、静定状態である場合、空気流量Ｇ_ａｔ、空気流量Ｇ_ａｓ、圧縮機入口温度Ｔ_１Ｃ、及び燃焼器車室圧力Ｐ_ｏｕｔの組み合わせを空気側記憶部２２２ａに記録し、空気側記憶部２２２ａに記録されている所定数の組み合わせを用いて空気補正係数Ｋ_ａを更新する。一方、空気流量補正部２２０ａは、状態検出部２０５が検出する状態が、静定状態でない場合、空気側記憶部２２２ａへの記録を行わず、空気補正係数Ｋ_ａの更新も行わない。そのため、第２の実施形態の構成は、第１の実施形態の構成と同様に、静定状態と過渡状態の両方の状態において、タービン入口温度算出部２０６に与える燃料流量と空気流量の精度が低くならないようにして、タービン入口温度算出部２０６が、一意に定めることが難しいタービン効率を用いることなく、長期間にわたってタービン入口温度Ｔ_１Ｔの精度を維持することを可能にしている。

　さらに、第２の実施形態では、燃料流量補正部２１０ａにおいて、燃料流量Ｇ_ｆｔ、燃料流量Ｇ_ｆｓ、燃焼器車室温度Ｔ_ＣＳ、及び燃焼器車室圧力Ｐ_ｏｕｔの所定数の組み合わせを用いて、説明変数を燃焼器車室温度Ｔ_ＣＳ及び燃焼器車室圧力Ｐ_ｏｕｔとし、目的変数を燃料補正係数Ｋ_ｆとする所定の回帰分析を行っている。そのため、燃料補正係数Ｋ_ｆを、燃焼器車室温度Ｔ_ＣＳや燃焼器車室圧力Ｐ_ｏｕｔといったガスタービン発電プラント１００ａの状態を示す状態量に応じて可変にすることができる。同様に、空気流量補正部２２０ａにおいて、空気流量Ｇ_ａｔ、空気流量Ｇ_ａｓ、圧縮機入口温度Ｔ_１Ｃ、及び燃焼器車室圧力Ｐ_ｏｕｔの所定数の組み合わせを用いて、説明変数を圧縮機入口温度Ｔ_１Ｃ及び燃焼器車室圧力Ｐ_ｏｕｔとし、目的変数を空気補正係数Ｋ_ａとする所定の回帰分析を行っている。そのため、空気補正係数Ｋ_ａを、圧縮機入口温度Ｔ_１Ｃや燃焼器車室圧力Ｐ_ｏｕｔといったガスタービン発電プラント１００ａの状態を示す状態量に応じて可変にすることができる。それにより、ガスタービン発電プラント１００ａの状態に応じた信頼性の面で精度の高い燃料補正係数Ｋ_ｆ及び空気補正係数Ｋ_ａを得ることができる。精度の高い燃料補正係数Ｋ_ｆ及び空気補正係数Ｋ_ａを用いて補正を行うことにより、補正によって得られる補正燃料流量Ｇ_ｆｃｏｒ及び補正空気流量Ｇ_ａｃｏｒの精度を高めることが可能になる。したがって、精度の高い補正燃料流量Ｇ_ｆｃｏｒ及び補正空気流量Ｇ_ａｃｏｒに基づいてタービン入口温度算出部２０６が算出するタービン入口温度Ｔ_１Ｔの精度も向上する。

　図１２のステップＳａ１５－１，Ｓａ１５－２の処理、すなわち燃料補正係数Ｋ_ｆ＝１として出力する処理、及び空気補正係数Ｋ_ａ＝１として出力する処理は、燃料側記憶部２１２ａ及び空気側記憶部２２２ａに記憶されている組み合わせの個数が、所定値未満である間に行われる処理であり、回帰式の係数を算出することができない間、燃料流量Ｇ_ｆｔと空気流量Ｇ_ａｔの各々に対して補正を行わないことを意図する処理である。また、ステップＳａ１４－１，Ｓａ１４－２の処理、すなわち前回の係数を回帰式に適用する用いる処理は、状態検出部２０５が検出した状態が、静定状態でない場合、すなわち過渡状態の場合に行われる処理であり、過渡状態の場合には、回帰式の係数の更新を行わないことを意図する処理である。

（第２の実施形態におけるその他の構成例）
　上記の第２の実施形態では、所定の回帰分析として、線形重回帰分析を適用しているが、所定の回帰分析として、線形重回帰分析以外の任意の回帰分析の手法を適用するようにしてもよい。例えば、ＳＶＲ（Support Vector Regression)やランダムフォレストなどの機械学習による回帰分析の手法を適用するようにしてもよい。

　上記の第２の実施形態では、燃料側回帰分析部２１５は、説明変数を燃焼器車室温度Ｔ_ＣＳ及び燃焼器車室圧力Ｐ_ｏｕｔとして所定の回帰分析を行うようにしている。また、空気側回帰分析部２２５は、説明変数を圧縮機入口温度Ｔ_１Ｃ及び燃焼器車室圧力Ｐ_ｏｕｔとして所定の回帰分析を行うようにしている。ただし、これらの説明変数は、一例として示したものであり、燃料側回帰分析部２１５は、燃焼器車室温度Ｔ_ＣＳ及び燃焼器車室圧力Ｐ_ｏｕｔのいずれか一方を説明変数としてもよいし、空気側回帰分析部２２５は、圧縮機入口温度Ｔ_１Ｃ及び燃焼器車室圧力Ｐ_ｏｕｔのいずれか一方を説明変数としてもよいし、ガスタービン発電プラント１００ａに備えられている任意のセンサ２５ａ～２５ｍが計測する計測データのいずれか１つ、または、複数を組み合わせたものを説明変数としてもよい。また、ガスタービン制御装置２０ａがガスタービン１０の制御のために出力する任意の指令値、例えば、上記した燃料制御信号指令値（ＣＳＯ）などを説明変数としてもよいし、計測データと、指令値とを組み合わせたものを説明変数としてもよい。

　上記の第２の実施形態では、図１１のステップＳａ１０５－１の処理において、燃料側管理部２１１ａは、燃料側記憶部２１２ａから所定数の燃料流量Ｇ_ｆｔ、燃料流量Ｇ_ｆｓ、燃焼器車室温度Ｔ_ＣＳ、及び燃焼器車室圧力Ｐ_ｏｕｔの組み合わせを読み出して燃料側回帰分析部２１５に出力するようにしている。これに対して、燃料側管理部２１１ａは、ステップＳａ１０５－１の処理において、燃料側記憶部２１２ａに記憶されている全ての組み合わせを読み出して燃料側回帰分析部２１５に出力するようにしてもよい。同様に、図１１のステップＳａ１０５－２の処理において、空気側管理部２２１ａは、空気側記憶部２２２ａから所定数の空気流量Ｇ_ａｔ、空気流量Ｇ_ａｓ、圧縮機入口温度Ｔ_１Ｃ、及び燃焼器車室圧力Ｐ_ｏｕｔの組み合わせを読み出して空気側回帰分析部２２５に出力するようにしている。これに対して、空気側管理部２２１ａは、ステップＳａ１０５－２の処理において、空気側記憶部２２２ａに記憶されている全ての組み合わせを読み出して空気側回帰分析部２２５に出力するようにしてもよい。

　上記の第２の実施形態では、補正処理部２０４ａは、燃料流量補正部２１０ａと、空気流量補正部２２０ａとを備えている。これに対して、補正処理部２０４ａは、燃料流量補正部２１０ａのみを備えるようにしてもよい。この場合、ヒートバランス演算部２０３は、燃料流量Ｇ_ｆｓのみを算出して出力し、空気流量算出部２０２が算出する空気流量Ｇ_ａｔが、タービン入口温度算出部２０６に与えられることになる。また、補正処理部２０４ａは、空気流量補正部２２０ａのみを備えるようにしてもよい。この場合、ヒートバランス演算部２０３は、空気流量Ｇ_ａｓのみを算出して出力し、燃料流量算出部２０１が算出する燃料流量Ｇ_ｆｔが、タービン入口温度算出部２０６に与えられることになる。また、補正処理部２０４ａが、燃料流量補正部２１０ａと、第１の実施形態の空気流量補正部２２０とを備えるようにしてもよいし、空気流量補正部２２０ａと、第１の実施形態の燃料流量補正部２１０とを備えるようにしてもよい。

＜第３の実施形態＞
　第１及び第２の実施形態では、燃料流量算出部２０１が算出する燃料流量Ｇ_ｆｔ及びヒートバランス演算部２０３が算出する燃料流量Ｇ_ｆｓは、３つの燃料供給系統４０－１，４０－２，４０－３の各々によって供給される燃料流量の合計値になっている。これに対して、第１及び第２の実施形態のガスタービン発電プラント１００，１００ａは、３つの燃料供給系統４０－１，４０－２，４０－３を備えていることから、燃料供給系統４０－１，４０－２，４０－３の各々に対する燃料流量Ｇ_ｆｔと、燃料流量Ｇ_ｆｓとを算出して、個々の燃料流量Ｇ_ｆｔに対して個別に補正を行う方が、より精度の高いタービン入口温度Ｔ_１Ｔを得ることができると考えられる。第３の実施形態のガスタービン制御装置２０ｂは、以下に説明するように、このような個別の補正を行う構成を備える。以下、説明の便宜上、第１の実施形態のガスタービン発電プラント１００において、ガスタービン制御装置２０を、第３の実施形態のガスタービン制御装置２０ｂに置き換えた構成をガスタービン発電プラント１００ｂという。

（第３の実施形態のガスタービン制御装置の構成例）
　図１３は、第３の実施形態のガスタービン制御装置２０ｂの構成を示すブロック図である。ガスタービン制御装置２０ｂは、燃料流量算出部２０１ｂ、空気流量算出部２０２、ヒートバランス演算部２０３ｂ、補正処理部２０４ｂ、状態検出部２０５、タービン入口温度算出部２０６、燃料配分比算出部２０７、弁開度算出部２０８、及び燃料配分量算出部２０９を備える。燃料流量算出部２０１ｂは、以下に示す点以外については、第１の実施形態の燃料流量算出部２０１と同一の構成を備える。第１の実施形態の燃料流量算出部２０１は、燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｔ（２），Ｇ_ｆｔ（３）を算出し、算出した燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｔ（２），Ｇ_ｆｔ（３）を合計して燃料流量Ｇ_ｆｔを算出している。これに対して、第３の実施形態の燃料流量算出部２０１ｂは、合計することなく、算出した燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｔ（２），Ｇ_ｆｔ（３）を出力する。ヒートバランス演算部２０３ｂは、以下に示す点以外については、第１の実施形態のヒートバランス演算部２０３と同一の構成を備える。第１の実施形態のヒートバランス演算部２０３は、燃料流量Ｇ_ｆｓと空気流量Ｇ_ａｓを算出して出力している。これに対して、第３の実施形態のヒートバランス演算部２０３ｂは、燃料流量Ｇ_ｆｓのみを算出して出力する。

　燃料配分量算出部２０９は、ヒートバランス演算部２０３ｂが算出する燃料流量Ｇ_ｆｓが、燃料配分比算出部２０７が算出する燃料配分比Ｄ_１，Ｄ_２，Ｄ_３の各々の比率になるように、３つの燃料流量Ｇ_ｆｓ（１），Ｇ_ｆｓ（２），Ｇ_ｆｓ（３）を算出する。燃料配分比Ｄ_１は、燃料供給系統４０－１に対応しているため、燃料配分比Ｄ_１に対応して算出する燃料流量Ｇ_ｆｓ（１）は、燃料供給系統４０－１に対応することになる。同様に、燃料供給系統４０－２と、燃料配分比Ｄ_２と、燃料流量Ｇ_ｆｓ（２）とが対応し、燃料供給系統４０－３と、燃料配分比Ｄ_３と、燃料流量Ｇ_ｆｓ（３）とが対応することになる。

　補正処理部２０４ｂは、燃料流量補正部２１０ｂを備える。燃料流量補正部２１０ｂは、燃料側管理部２１１ｂ、燃料側記憶部２１２ｂ、燃料補正係数算出部２１３－１，２１３－２，２１３－３、燃料補正係数適用部２１４－１，２１４－２，２１４－３、及び加算器２１６を備える。燃料側記憶部２１２ｂは、自らに予め設けられる記憶領域である燃料側記憶領域２１２ｂ－１，２１２ｂ－２，２１２ｂ－３を備える。

　燃料側管理部２１１ｂは、燃料供給系統４０－１に対応する燃料流量Ｇ_ｆｔ（１）と、燃料流量Ｇ_ｆｓ（１）とを関連付ける。燃料側管理部２１１ｂは、燃料供給系統４０－２に対応する燃料流量Ｇ_ｆｔ（２）と、燃料流量Ｇ_ｆｓ（２）とを関連付ける。燃料側管理部２１１ｂは、燃料供給系統４０－３に対応する燃料流量Ｇ_ｆｔ（３）と、燃料流量Ｇ_ｆｓ（３）とを関連付ける。燃料側管理部２１１ｂは、状態検出部２０５が検出する状態にしたがって、関連付けた燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｓ（１）と、関連付けた燃料流量Ｇ_ｆｔ（２），Ｇ_ｆｓ（２）と、関連付けた燃料流量Ｇ_ｆｔ（３），Ｇ_ｆｓ（３）とを、それぞれ燃料側記憶領域２１２ｂ－１，２１２ｂ－２，２１２ｂ－３に書き込んで記憶させる。燃料側管理部２１１ｂは、最新の燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｓ（１）の組み合わせと、最新の燃料流量Ｇ_ｆｔ（２），Ｇ_ｆｓ（２）の組み合わせと、最新の燃料流量Ｇ_ｆｔ（３），Ｇ_ｆｓ（３）の組み合わせとを、それぞれ燃料補正係数算出部２１３－１，２１３－２，２１３－３に出力する。

　燃料側記憶部２１２ｂにおいて、燃料側記憶領域２１２ｂ－１は、燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｓ（１）の組み合わせを時系列順に記憶する。燃料側記憶領域２１２ｂ－２は、燃料流量Ｇ_ｆｔ（２），Ｇ_ｆｓ（２）の組み合わせを時系列順に記憶する。燃料側記憶領域２１２ｂ－３は、燃料流量Ｇ_ｆｔ（３），Ｇ_ｆｓ（３）の組み合わせを時系列順に記憶する。燃料補正係数算出部２１３－１，２１３－２，２１３－３の各々は、第１の実施形態の燃料補正係数算出部２１３と同一の構成である。燃料補正係数適用部２１４－１，２１４－２，２１４－３の各々は、第１の実施形態の燃料補正係数適用部２１４と同一の構成である。加算器２１６は、燃料補正係数適用部２１４－１，２１４－２，２１４－３の各々が算出する補正燃料流量Ｇ_ｆｃｏｒ（１），Ｇ_ｆｃｏｒ（２），Ｇ_ｆｃｏｒ（３）を合計して補正燃料流量Ｇ_ｆｃｏｒを算出する。

（第３の実施形態のガスタービン制御装置の動作例）
　図１４，図１５を参照しつつガスタービン制御装置２０ｂによる処理について説明する。図１４において、ステップＳｂ１－１の処理は、図３のステップＳ１－１と同一の処理が燃料流量算出部２０１ｂによって行われる。ステップＳｂ１－２の処理は、図３のステップＳ１－２と同一の処理がヒートバランス演算部２０３ｂによって行われる。燃料配分量算出部２０９は、前回の検出周期において燃料配分比算出部２０７が算出して出力した燃料配分比Ｄ_１，Ｄ_２，Ｄ_３を取り込む（ステップＳｂ１－３）。

　燃料流量算出部２０１ｂは、ステップＳｂ１－１の処理において取り込んだデータを、式（１），（２），（３）に代入して、燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｔ（２），Ｇ_ｆｔ（３）を算出する。燃料流量算出部２０１ｂは、算出した燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｔ（２），Ｇ_ｆｔ（３）を、燃料側管理部２１１ｂに出力する。燃料流量算出部２０１ｂは、算出した燃料流量Ｇ_ｆｔ（１）を、燃料補正係数適用部２１４－１に出力し、燃料流量Ｇ_ｆｔ（２）を、燃料補正係数適用部２１４－２に出力し、燃料流量Ｇ_ｆｔ（３）を、燃料補正係数適用部２１４－３に出力する（Ｓｂ２－１）。

　ヒートバランス演算部２０３ｂは、ステップＳｂ１－２の処理において取り込んだデータと、ガスタービンヒートバランス計算式とに基づいて、燃料流量Ｇ_ｆｓを算出する。ヒートバランス演算部２０３ｂは、算出した燃料流量Ｇ_ｆｓを燃料配分量算出部２０９に出力する（Ｓｂ２－２）。

　燃料配分量算出部２０９は、ヒートバランス演算部２０３ｂが算出する燃料流量Ｇ_ｆｓが、ステップＳｂ１－３の処理において取り込んだ燃料配分比Ｄ_１，Ｄ_２，Ｄ_３の各々の比率になるように、３つの燃料流量Ｇ_ｆｓ（１），Ｇ_ｆｓ（２），Ｇ_ｆｓ（３）を算出する。例えば、燃料配分量算出部２０９は、Ｇ_ｆｓ（１）＝Ｇ_ｆｓ×Ｄ_１／（Ｄ_１＋Ｄ_２＋Ｄ_３）、Ｇ_ｆｓ（２）＝Ｇ_ｆｓ×Ｄ_２／（Ｄ_１＋Ｄ_２＋Ｄ_３）、Ｇ_ｆｓ（３）＝Ｇ_ｆｓ×Ｄ_３／（Ｄ_１＋Ｄ_２＋Ｄ_３）として算出する。燃料配分量算出部２０９は、算出した燃料流量Ｇ_ｆｓ（１），Ｇ_ｆｓ（２），Ｇ_ｆｓ（３）を燃料側管理部２１１ｂに出力する（Ｓｂ２－３）。

　なお、ステップＳｂ１－１，Ｓｂ１－２の処理のセットと、ステップＳｂ１－２，Ｓｂ２－２の処理のセットと、ステップＳｂ１－３の処理とは、並列に行われる。

　補正処理部２０４ｂの燃料流量補正部２１０ｂは、燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｔ（２），Ｇ_ｆｔ（３）と、燃料流量Ｇ_ｆｓ（１），Ｇ_ｆｓ（２），Ｇ_ｆｓ（３）とを受けると、図１５に示す記憶部管理処理のサブルーチンを開始する（Ｓｂ３）。

　燃料流量補正部２１０ｂにおいて、燃料側管理部２１１ｂは、燃料流量算出部２０１ｂが出力する燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｔ（２），Ｇ_ｆｔ（３）と、燃料配分量算出部２０９が出力する燃料流量Ｇ_ｆｓ（１），Ｇ_ｆｓ（２），Ｇ_ｆｓ（３）とを取り込む。燃料流量補正部２１０ｂは、取り込んだ燃料流量Ｇ_ｆｔ（１）と、燃料流量Ｇ_ｆｓ（１）とを関連付け、取り込んだ燃料流量Ｇ_ｆｔ（２）と、燃料流量Ｇ_ｆｓ（２）とを関連付け、取り込んだ燃料流量Ｇ_ｆｔ（３）と、燃料流量Ｇ_ｆｓ（３）とを関連付ける（Ｓｂ１００）。燃料側管理部２１１ｂは、ステップＳｂ１００の処理を終了した時点で状態検出部２０５が検出するガスタービン発電プラント１００ｂの状態を示すデータを状態検出部２０５から取得する（Ｓｂ１０１）。

　燃料側管理部２１１ｂは、取得した状態を示すデータが、静定状態を示しているか否かを判定する（Ｓｂ１０２）。燃料側管理部２１１ｂは、取得した状態を示すデータが、静定状態を示していると判定した場合（Ｓｂ１０２、Ｙｅｓ）、ステップＳｂ１０３－１，Ｓｂ１０３－２，Ｓｂ１０３－３の処理を行い、その後、ステップＳｂ１０４－１，Ｓｂ１０４－２，Ｓｂ１０４－３の処理を行う。

　すなわち、燃料側管理部２１１ｂは、関連付けた燃料流量Ｇ_ｆｔ（１）と燃料流量Ｇ_ｆｓ（１）の組み合わせを燃料側記憶領域２１２ｂ－１に書き込んで記憶させる。燃料側管理部２１１ｂは、燃料側記憶領域２１２ｂ－１に順次書き込む燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｓ（１）の組み合わせが燃料側記憶領域２１２ｂ－１において時系列順に記憶されるように、例えば、新たに書き込む組み合わせが、時系列順の先頭になるように燃料側記憶領域２１２ｂ－１に書き込む（Ｓｂ１０３－１）。燃料側管理部２１１ｂは、燃料側記憶領域２１２ｂ－１に書き込んだ燃料流量Ｇ_ｆｔ（１）と燃料流量Ｇ_ｆｓ（１）の組み合わせを燃料補正係数算出部２１３－１に出力する（Ｓｂ１０４－１）。

　燃料側管理部２１１ｂは、ステップＳｂ１０３－２，Ｓｂ１０４－２の処理として、ステップＳｂ１０３－１，Ｓｂ１０４－１の処理において、燃料流量Ｇ_ｆｔ（１）を燃料流量Ｇ_ｆｔ（２）に読み替え、燃料流量Ｇ_ｆｓ（１）を燃料流量Ｇ_ｆｓ（２）に読み替え、燃料側記憶領域２１２ｂ－１を燃料側記憶領域２１２ｂ－２に読み替え、燃料補正係数算出部２１３－１を燃料補正係数算出部２１３－２に読み替えた処理を行う。同様に、燃料側管理部２１１ｂは、ステップＳｂ１０３－３，Ｓｂ１０４－３の処理として、ステップＳｂ１０３－１，Ｓｂ１０４－１の処理において、燃料流量Ｇ_ｆｔ（１）を燃料流量Ｇ_ｆｔ（３）に読み替え、燃料流量Ｇ_ｆｓ（１）を燃料流量Ｇ_ｆｓ（３）に読み替え、燃料側記憶領域２１２ｂ－１を燃料側記憶領域２１２ｂ－３に読み替え、燃料補正係数算出部２１３－１を燃料補正係数算出部２１３－３に読み替えた処理を行う。

　燃料側管理部２１１ｂは、ステップＳｂ１０２の処理において、取得した状態を示すデータが、静定状態を示していないと判定した場合（Ｓｂ１０２、Ｎｏ）ステップＳｂ１００の処理で取り込んだ燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｔ（２），Ｇ_ｆｔ（３）と、燃料流量Ｇ_ｆｓ（１），Ｇ_ｆｓ（２），Ｇ_ｆｓ（３）とを破棄し、燃料側記憶領域２１２ｂ－１，２１２ｂ－２，２１２ｂ－３の各々において記憶されている組み合わせのうち最新の組み合わせを読み出す。上記したように燃料側記憶領域２１２ｂ－１，２１２ｂ－２，２１２ｂ－３が最新の組み合わせが先頭になるように時系列順に記憶している場合、燃料側管理部２１１ｂは、燃料側記憶領域２１２ｂ－１，２１２ｂ－２，２１２ｂ－３の各々の先頭に記憶されている組み合わせを最新の組み合わせとして読み出す。燃料側管理部２１１ｂは、燃料側記憶領域２１２ｂ－１から読み出した最新の燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｓ（１）の組み合わせを燃料補正係数算出部２１３－１に出力し、燃料側記憶領域２１２ｂ－２から読み出した最新の燃料流量Ｇ_ｆｔ（２），Ｇ_ｆｓ（２）の組み合わせを燃料補正係数算出部２１３－２に出力し、燃料側記憶領域２１２ｂ－３から読み出した最新の燃料流量Ｇ_ｆｔ（３），Ｇ_ｆｓ（３）の組み合わせを燃料補正係数算出部２１３－３に出力する（Ｓｂ１０５）。

　なお、図１５の処理において、燃料側管理部２１１ｂは、ステップＳｂ１０３－１，Ｓｂ１０４－１の処理のセットと、ステップＳｂ１０３－２，Ｓｂ１０４－２の処理のセットと、ステップＳｂ１０３－３，Ｓｂ１０４－３の処理のセットとをいずれの順番で行うようにしてもよいし、並列に行うようにしてもよい。また、燃料側管理部２１１ｂは、ステップＳｂ１０５の処理において、燃料側記憶領域２１２ｂ－１，２１２ｂ－２，２１２ｂ－３の各々からの最新の組み合わせを読み出して燃料補正係数算出部２１３－１，２１３－２，２１３－３の各々に出力する処理を、いずれの順番で行ってもよいし、並列に行うようにしてもよい。

　図１４に戻り、燃料補正係数算出部２１３－１は、燃料側管理部２１１ｂが出力する燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｓ（１）の組み合わせを取り込む（Ｓｂ４－１）。燃料補正係数算出部２１３－１は、燃料流量Ｇ_ｆｓ（１）を燃料流量Ｇ_ｆｔ（１）で除算して燃料補正係数Ｋ_ｆ（１）を算出する。燃料補正係数算出部２１３－１は、算出した燃料補正係数Ｋ_ｆ（１）を燃料補正係数適用部２１４－１に出力する（Ｓｂ５－１）。燃料補正係数適用部２１４－１は、燃料流量算出部２０１ｂが出力する燃料流量Ｇ_ｆｔ（１）と、燃料補正係数算出部２１３－１が出力する燃料補正係数Ｋ_ｆ（１）を取り込む。燃料補正係数適用部２１４－１は、燃料流量Ｇ_ｆｔ（１）に、燃料補正係数Ｋ_ｆ（１）を乗算して補正燃料流量Ｇ_ｆｃｏｒ（１）を算出する。燃料補正係数適用部２１４－１は、算出した補正燃料流量Ｇ_ｆｃｏｒ（１）を加算器２１６に出力する（Ｓｂ６－１）。

　同様に、ステップＳｂ４－２，Ｓｂ５－２，Ｓｂ６－２の処理と、ステップＳｂ４－３，Ｓｂ５－３，Ｓｂ６－３の処理とが行われる。ここで、ステップＳｂ４－２，Ｓｂ５－２，Ｓｂ６－２の処理の各々は、各々のステップ番号の枝番号を「－１」に置き換えたステップＳｂ４－１，Ｓｂ５－１，Ｓｂ６－１の処理において、燃料流量Ｇ_ｆｔ（１）を燃料流量Ｇ_ｆｔ（２）に読み替え、燃料流量Ｇ_ｆｓ（１）を燃料流量Ｇ_ｆｓ（２）に読み替え、燃料補正係数算出部２１３－１を燃料補正係数算出部２１３－２に読み替え、燃料補正係数Ｋ_ｆ（１）を燃料補正係数Ｋ_ｆ（２）に読み替え、燃料補正係数適用部２１４－１を燃料補正係数適用部２１４－２に読み替え、補正燃料流量Ｇ_ｆｃｏｒ（１）を補正燃料流量Ｇ_ｆｃｏｒ（２）に読み替えた処理である。また、ステップＳｂ４－３，Ｓｂ５－３，Ｓｂ６－３の処理の各々は、各々のステップ番号の枝番号を「－１」に置き換えたステップＳｂ４－１，Ｓｂ５－１，Ｓｂ６－１の処理において、燃料流量Ｇ_ｆｔ（１）を燃料流量Ｇ_ｆｔ（３）に読み替え、燃料流量Ｇ_ｆｓ（１）を燃料流量Ｇ_ｆｓ（３）に読み替え、燃料補正係数算出部２１３－１を燃料補正係数算出部２１３－３に読み替え、燃料補正係数Ｋ_ｆ（１）を燃料補正係数Ｋ_ｆ（３）に読み替え、燃料補正係数適用部２１４－１を燃料補正係数適用部２１４－３に読み替え、補正燃料流量Ｇ_ｆｃｏｒ（１）を補正燃料流量Ｇ_ｆｃｏｒ（３）に読み替えた処理である。

　加算器２１６は、燃料補正係数適用部２１４－１，２１４－２，２１４－３の各々が出力する補正燃料流量Ｇ_ｆｃｏｒ（１），Ｇ_ｆｃｏｒ（２），Ｇ_ｆｃｏｒ（３）を取り込み、取り込んだ補正燃料流量Ｇ_ｆｃｏｒ（１），Ｇ_ｆｃｏｒ（２），Ｇ_ｆｃｏｒ（３）を合計して補正燃料流量Ｇ_ｆｃｏｒを算出する。加算器２１６は、算出した補正燃料流量Ｇ_ｆｃｏｒをタービン入口温度算出部２０６に出力する（Ｓｂ７）。図１４に示す「Ａ」のマークは、処理が継続して行われることを示すマークであり、「Ａ」のマークが示すように、その後、第１の実施形態の図４の処理と同一の処理が行われる。上記の図１４及び図４に示す一連の処理が、１回の検出周期の処理として行われ、周期ごとに図１４及び図４に示す一連の処理が行われることになる。

（第３の実施形態における作用・効果）
　上記したように、第３の実施形態の燃料流量補正部２１０ｂは、状態検出部２０５が検出する状態が、静定状態と過渡状態のいずれの状態であっても、直近の静定状態の燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｓ（１）、燃料流量Ｇ_ｆｔ（２），Ｇ_ｆｓ（２）、及び燃料流量Ｇ_ｆｔ（３），Ｇ_ｆｓ（３）を用いて燃料補正係数Ｋ_ｆを算出する。そのため、第３の実施形態の構成は、第１の実施形態の構成と同様に、静定状態と過渡状態の両方の状態において、タービン入口温度算出部２０６に与える燃料流量の精度が低くならないようにして、タービン入口温度算出部２０６が、一意に定めることが難しいタービン効率を用いることなく、長期間にわたってタービン入口温度Ｔ_１Ｔの精度を維持することを可能にしている。

　さらに、第３の実施形態では、燃料配分量算出部２０９は、ヒートバランス演算部２０３ｂが算出する燃料流量Ｇ_ｆｓを、燃料配分比Ｄ_１，Ｄ_２，Ｄ_３にしたがって配分した燃料流量Ｇ_ｆｓ（１），Ｇ_ｆｓ（２），Ｇ_ｆｓ（３）を算出する。燃料補正係数算出部２１３－１，２１３－２，２１３－３の各々は、各々に対応する燃料流量Ｇ_ｆｓ（１），Ｇ_ｆｓ（２），Ｇ_ｆｓ（３）に基づいて、燃料補正係数Ｋ_ｆ（１），Ｋ_ｆ（２），Ｋ_ｆ（３）を算出する。燃料補正係数適用部２１４－１～２１４－３の各々は、燃料供給系統４０－１，４０－２，４０－３ごとの燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｔ（２），Ｇ_ｆｔ（３）の各々を、各々に対応する燃料補正係数Ｋ_ｆ（１），Ｋ_ｆ（２），Ｋ_ｆ（３）によって補正する。そのため、第３の実施形態では、燃料供給系統４０－１，４０－２，４０－３の各々に存在する固有の誤差、すなわち、静定状態の際に記録した燃料流量Ｇ_ｆｓ（１）と燃料流量Ｇ_ｆｔ（１）の誤差、燃料流量Ｇ_ｆｓ（２）と燃料流量Ｇ_ｆｔ（２）の誤差、及び燃料流量Ｇ_ｆｓ（３）と燃料流量Ｇ_ｆｔ（３）の誤差を、個々に補正することができる。したがって、第３の実施形態では、第１の実施形態の構成において補正処理部２０４が燃料流量補正部２１０のみを備えた場合の構成よりも、高い精度で燃料流量Ｇ_ｆｔの補正を行うことを可能とし、それによって、タービン入口温度算出部２０６が算出するタービン入口温度Ｔ_１Ｔの精度をより向上させることが可能になる。

（第３の実施形態におけるその他の構成例）
　上記の第３の実施形態では、燃料配分量算出部２０９は、燃料配分比算出部２０７が算出する燃料配分比Ｄ_１，Ｄ_２，Ｄ_３を用いて、ヒートバランス演算部２０３ｂが算出する燃料流量Ｇ_ｆｓから燃料流量Ｇ_ｆｓ（１），Ｇ_ｆｓ（２），Ｇ_ｆｓ（３）を算出するようにしている。これに対して、燃料配分量算出部２０９は、弁開度算出部２０８が算出する弁開度Ｏ_１，Ｏ_２，Ｏ_３を用いて、燃料流量Ｇ_ｆｓから燃料流量Ｇ_ｆｓ（１），Ｇ_ｆｓ（２），Ｇ_ｆｓ（３）を算出するようにしてもよい。

　また、燃料供給系統４０－１，４０－２，４０－３の各々において、燃料供給系統４０－１，４０－２，４０－３の各々を流れる燃料の流量を計測するセンサを備え、当該センサから得られる燃料供給系統４０－１，４０－２，４０－３の各々の燃料流量を用いて、燃料流量Ｇ_ｆｓから燃料流量Ｇ_ｆｓ（１），Ｇ_ｆｓ（２），Ｇ_ｆｓ（３）を算出するようにしてもよい。ただし、燃料供給系統４０－１，４０－２，４０－３の各々を流れる燃料の流量をセンサによって計測することにより得られる燃料流量のデータは、燃料供給系統４０－１，４０－２，４０－３の各々を流れる凡その比率を知るには十分な精度のデータであるが、計測における絶対的な誤差を含んだデータであるため、燃料流量算出部２０１ｂが算出する燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｔ（２），Ｇ_ｆｔ（３）と異なる値になり、また、当該データの合計値をタービン入口温度算出部２０６に与えたとしても、精度の良いタービン入口温度Ｔ_１Ｔが得られるデータではない。

　上記の第３の実施形態において、燃料流量補正部２１０ｂに対して、第１の実施形態において説明した「燃料流量と空気流量の各々の平均値より各々の補正係数を算出する手法」を適用するようにしてもよい。また、燃料流量補正部２１０ｂに対して、第２の実施形態の構成を適用するようにしてもよい。また、第３の実施形態において、ヒートバランス演算部２０３ｂを、第１の実施形態のヒートバランス演算部２０３に置き換えた上で、補正処理部２０４ｂが、第１の実施形態の空気流量補正部２２０を備えるようにしてもよいし、補正処理部２０４ｂが、第２の実施形態の空気流量補正部２２０ａを備えるようにしてもよい。

（第１から第３の実施形態におけるその他の構成例）
　上記の第１から第３の実施形態において、ガスタービン発電プラント１００，１００ａ，１００ｂは、燃料供給系統４０－１，４０－２，４０－３の各々に対応する燃料温度センサ２５ｃ－１，２５ｃ－２，２５ｃ－３を備えている。例えば、燃料供給系統４０－１，４０－２，４０－３の各々を流れる燃料の温度に大きな差がない場合には、燃料温度センサ２５ｃ－１，２５ｃ－２，２５ｃ－３のいずれか１つを備えるようにしてもよいし、燃料供給装置３０において３つに分岐される前の燃料の温度を測定するセンサを備えるようにしてもよい。例えば、燃料温度センサ２５ｃ－１のみを備える場合、燃料温度Ｔ_ｆ（２），Ｔ_ｆ（３）の各々を燃料温度Ｔ_ｆ（１）と同一値であるとみなして、上記した第１から第３の実施形態における処理が行われることになる。一例として、燃料流量算出部２０１，２０１ｂ、及びタービン入口温度算出部２０６には、燃料温度Ｔ_ｆ（１）のみが与えられることになり、燃料流量算出部２０１，２０１ｂは、式（２），式（３）の演算を行う場合、燃料温度Ｔ_ｆ（２），Ｔ_ｆ（３）に替えて燃料温度Ｔ_ｆ（１）を代入することになる。また、タービン入口温度算出部２０６は、式（６）に対して燃料温度Ｔ_ｆを代入する際に、燃料温度Ｔ_ｆ（１），Ｔ_ｆ（２），Ｔ_ｆ（３）の平均値を算出する必要はなく、燃料温度Ｔ_ｆ（１）を燃料温度Ｔ_ｆとして式（６）に代入することになる。

　上記の第１から第３の実施形態において、燃料流量算出部２０１，２０１ｂは、流調弁前圧力センサ２５ａ－１～２５ａ－３の各々が検出する流調弁前圧力Ｐ_１ＦＶ（１），Ｐ_１ＦＶ（２），Ｐ_１ＦＶ（３）と、流調弁後圧力センサ２５ｂ－１～２５ｂ－３の各々が検出する流調弁後圧力Ｐ_２ＦＶ（１），Ｐ_２ＦＶ（２），Ｐ_２ＦＶ（３）とに基づいて、燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｔ（２），Ｇ_ｆｔ（３）を算出する手法を採用している。ただし、当該手法は一例であり、燃料流量算出部２０１，２０１ｂは、燃料供給系統４０－１，４０－２，４０－３に関する計測データ、より詳細には、燃料供給系統４０－１，４０－２，４０－３の各々の２か所において計測することができる圧力に関する計測データがあれば、当該計測データが適用できるように変形した燃料流量算出関数ｆ_１（・）によって燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｔ（２），Ｇ_ｆｔ（３）を算出することができる。例えば、流調弁前圧力センサ２５ａ－１～２５ａ－３に替えて、ノズル４２－１，４２－２，４２－３に流入する燃料の圧力を計測するノズル前圧力センサと、流調弁後圧力センサ２５ｂ－１～２５ｂ－３に替えて、ノズル４２－１，４２－２，４２－３から流出する燃料の圧力を計測するノズル後圧力センサとを備えるようにしてもよい。この場合、燃料流量算出部２０１，２０１ｂは、ノズル前圧力センサが計測する圧力と、ノズル後圧力センサが計測する圧力とを取り込み、これら２つの圧力が適用できるように変形した燃料流量算出関数ｆ_１（・）と、弁開度Ｏ_１～Ｏ_３と、当該２つの圧力と、燃料温度Ｔ_ｆ（１）～Ｔ_ｆ（３）とによって燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｔ（２），Ｇ_ｆｔ（３）を算出することになる。また、例えば、流調弁前圧力センサ２５ａ－１～２５ａ－３に替えて、燃料供給系統４０－１，４０－２，４０－３の各々に対応するマニホールド配管４３－１，４３－２，４３－３内の燃料の圧力を計測する第１のセンサと、流調弁後圧力センサ２５ｂ－１～２５ｂ－３に替えて、マニホールド配管４３－１，４３－２，４３－３の各々から流出する燃料の圧力、すなわち、燃焼ノズルの出口圧力を計測する第２のセンサとを備えるようにしてもよい。なお、第２のセンサの一例は、第２の実施形態において示した燃焼器車室圧力Ｐ_ｏｕｔを計測する燃焼器車室圧力センサ２５ｍである。

　上記の第１から第３の実施形態では、ガスタービン制御装置２０，２０ａ，２０ｂが備える弁開度算出部２０８は、流調弁４１－１，４１－２，４１－３の各々に対応する弁開度Ｏ_１，Ｏ_２，Ｏ_３を算出するようにしている。これに対して、弁開度算出部２０８は、流調弁４１－１，４１－２，４１－３の各々に対応する３つのＣｖ値を算出するようにしてもよい。この場合、燃料流量算出部２０１，２０１ｂは、弁開度算出部２０８が算出する３つのＣｖ値を取り込み、Ｃｖ値が適用できるように変形した燃料流量算出関数ｆ_１（・）に弁開度Ｏ_１～Ｏ_３に替えて３つのＣｖ値を代入して燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｔ（２），Ｇ_ｆｔ（３）を算出することになる。

　上記の第１から第３の実施形態において、ヒートバランス演算部２０３，２０３ｂは、ガスタービン出力Ｇ_ｏｕｔと、排ガス圧力Ｐ_２Ｔと、排ガス温度Ｔ_２Ｔと、圧縮機入口圧力Ｐ_１Ｃと、圧縮機入口温度Ｔ_１Ｃと、圧縮機出口圧力Ｐ_２Ｃと、圧縮機出口温度Ｔ_２Ｃとを取り込むようにしているが、これら以外の計測データを取り込むようにしてもよい。

　上記の第１から第３の実施形態において、一例として、燃料供給系統４０が、３つの燃料供給系統４０－１～４０－３を備える場合を示しているが、２つの燃料供給系統４０－１，４０－２を備える構成や、４つ以上の燃料供給系統４０－１，４０－２，４０－３，…を備える構成などであってもよい。また、図１，図８では、一例として、ガスタービン１０が、１つの燃焼器１２を備える場合を示しているが、ガスタービン１０が、複数の燃焼器１２であって各々が圧縮機１１と、タービン１３とに連結する燃焼器１２を備え、複数の燃焼器１２の各々に対応する複数の燃焼器内筒１２１と、燃料供給系統４０－１，４０－２，４０－３，…とがそれぞれ備えられるように構成されていてもよい。また、１つの燃焼器１２において複数の燃焼器内筒１２１が備えられ、複数の燃焼器内筒１２１の各々に対応する燃料供給系統４０－１，４０－２，４０－３，…が備えられるように構成されていてもよい。このように、第１から第３の実施形態において示した燃料供給系統４０－１，４０－２，４０－３の個数や、燃焼器１２の個数や、燃焼器内筒１２１の個数を変える構成変更が行われる場合、当該構成変更に応じてガスタービン発電プラント１００，１００ａ，１００ｂにおいて備えられる流調弁４１などの部品の数が変わることになり、ガスタービン制御装置２０，２０ａ，２０ｂが備える機能部の構成も変わることになる。

　上記の第１から第３の実施形態において、ガスタービン制御装置２０，２０ａ，２０ｂは、図２，図９，図１３に示す構成に加えて、インレットガイドベーン１４に対して与えるインレットガイドベーン１４が吸入する空気量を調整する指令値を算出する構成を備えているものとする。

　上記の第１から第３の実施形態において、ガスタービン発電プラント１００，１００ａ，１００ｂは、流調弁４１－１，４１－２，４１－３を備えているが、流量調整弁に替えて圧力調整弁を適用するようにしてもよい。

　上記の第１から第３の実施形態において、ガスタービン発電プラント１００，１００ａ，１００ｂは、更にスチームタービンを備え、ロータ１５に当該スチームタービンが連結されたガスタービンコンバインドサイクル発電プラント（ＧＴＣＣ）として構成されていてもよい。

　上記の第１から第３の実施形態において、タービン入口温度算出部２０６が利用する関数ｆ_３（・）は、非定常物理モデルに基づいて構築された関数であり、非定常物理モデルにおける燃焼器１２に流入する熱エネルギは、上記したように、燃料の熱エネルギ、空気の熱エネルギ、燃焼ガスの発熱エネルギの合計により表される。ところで、空気の熱エネルギは、燃焼器車室１６に流入する空気に含まれる水蒸気による湿度の影響を考慮した比エンタルピーに応じて変動する。そのため、タービン入口温度算出部２０６は、関数ｆ_３（・）に代入するパラメータ、すなわち補正燃料流量Ｇ_ｆｃｏｒ、補正空気流量Ｇ_ａｃｏｒ、燃焼器車室温度Ｔ_ＣＳ、及び燃料温度Ｔ_ｆ（１），Ｔ_ｆ（２），Ｔ_ｆ（３）の平均値である燃料温度Ｔ_ｆのいずれかによって空気の熱エネルギを算出するのではなく、以下のようにして空気の熱エネルギを算出するようにしてもよい。すなわち、燃焼器車室１６に流入する空気の湿度を計測するセンサを新たに備え、タービン入口温度算出部２０６は、当該センサによって計測した燃焼器車室１６に流入する空気の湿度に基づいて比エンタルピーを算出し、算出した比エンタルピーを用いて空気の熱エネルギを算出するようにしてもよい。

　また、上記の第１から第３の実施形態の構成では、状態検出部２０５が、ガスタービン発電プラント１００の状態を検出する処理、図６のステップＳ１０５－１，Ｓ１０５－２の処理、及び図１１のステップＳａ１０４－１，Ｓａ１０４－２の処理において、等号付き不等号を用いた判定処理を行っている。しかしながら、本開示は、当該実施の形態に限られるものではなく、「以上であるか否か」という判定処理は一例に過ぎず、閾値の定め方に応じて、ぞれぞれ「超過するか否か」という判定処理に置き換えられてもよい。

　以上、本開示の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施の形態に限られるものではなく、本開示の要旨を逸脱しない範囲の設計等も含まれる。

＜コンピュータ構成＞
　上記の第１から第３の実施形態において、ガスタービン制御装置２０，２０ａ，２０ｂは、例えば、ＣＰＵ（Ｃｅｎｔｒａｌ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ）、ＲＡＭ（Ｒａｎｄｏｍ　Ａｃｃｅｓｓ　Ｍｅｍｏｒｙ）、ＲＯＭ（Ｒｅａｄ　Ｏｎｌｙ　Ｍｅｍｏｒｙ）、及びコンピュータ読み取り可能な記憶媒体等から構成されている。そして、ガスタービン制御装置２０，２０ａ，２０ｂが備える各種機能を実現するための一連の処理は、一例として、プログラムの形式で記憶媒体等に記憶されており、このプログラムをＣＰＵがＲＡＭ等に読み出して、情報の加工・演算処理を実行することにより、各種機能が実現される。なお、プログラムは、ＲＯＭやその他の記憶媒体に予めインストールしておく形態や、コンピュータ読み取り可能な記憶媒体に記憶された状態で提供される形態、有線又は無線による通信手段を介して配信される形態等が適用されてもよい。コンピュータ読み取り可能な記憶媒体とは、磁気ディスク、光磁気ディスク、ＣＤ－ＲＯＭ、ＤＶＤ－ＲＯＭ、半導体メモリ等である。

　より具体的には、上記の第１から第３の実施形態に係るガスタービン制御装置２０，２０ａ，２０ｂの各々は、図１６に示すコンピュータの構成によって実現される。コンピュータ９０は、プロセッサ９１、メインメモリ９２、ストレージ９３、及びインタフェース９４を備える。上記のガスタービン制御装置２０，２０ａ，２０ｂが備える燃料流量算出部２０１，２０１ｂ、ヒートバランス演算部２０３，２０３ｂ、空気流量算出部２０２、補正処理部２０４，２０４ａ，２０４ｂ、状態検出部２０５、タービン入口温度算出部２０６、燃料配分比算出部２０７、弁開度算出部２０８、及び燃料配分量算出部２０９の機能部は、コンピュータ９０に実装される。これらの機能部の動作は、プログラムの形式でストレージ９３に記憶される。プロセッサ９１は、プログラムをストレージ９３から読み出してメインメモリ９２に展開し、メインメモリ９２に展開されたプログラムにしたがって、第１から第３の実施形態において示した処理を実行する。また、プロセッサ９１は、プログラムにしたがって、例えば、燃料側記憶部２１２，２１２ａ，２１２ｂ、及び空気側記憶部２２２，２２２ａの領域をメインメモリ９２に確保する。また、プロセッサ９１は、プログラムにしたがって、ガスタービン制御装置２０，２０ａ，２０ｂが備える機能部に対して、各種のセンサ２５ａ～２０ｍ、インレットガイドベーン１４、及び流調弁４１との間でデータの送受信を行うインタフェース９４を割り当てる。

　プログラムは、コンピュータ９０に発揮させる機能の一部を実現するためのものであってもよい。例えば、プログラムは、ストレージ９３に既に記憶されている他のプログラムとの組み合わせ、または他の装置に実装された他のプログラムとの組み合わせによって機能を発揮させるものであってもよい。なお、他の実施形態においては、コンピュータは、上記構成に加えて、または上記構成に代えてＰＬＤ（Programmable Logic Device）などのカスタムＬＳＩ（Large Scale Integrated Circuit）を備えてもよい。ＰＬＤの例としては、ＰＡＬ(Programmable Array Logic)、ＧＡＬ(Generic Array Logic)、ＣＰＬＤ(Complex Programmable Logic Device)、ＦＰＧＡ（Field Programmable Gate Array）等が挙げられる。この場合、プロセッサによって実現される機能の一部または全部が当該集積回路によって実現されてよい。

　ストレージ９３の例としては、ＨＤＤ（Hard Disk Drive）、ＳＳＤ（Solid State Drive）、磁気ディスク、光磁気ディスク、ＣＤ－ＲＯＭ（Compact Disc Read Only Memory）、ＤＶＤ－ＲＯＭ（Digital Versatile Disc Read Only Memory）、半導体メモリ等が挙げられる。ストレージ９３は、コンピュータ９０のバスに直接接続された内部メディアであってもよいし、インタフェース９４または通信回線を介してコンピュータ９０に接続される外部メディアであってもよい。また、このプログラムが通信回線によってコンピュータ９０に配信される場合、配信を受けたコンピュータ９０が当該プログラムをメインメモリ９２に展開し、上記処理を実行してもよい。また、ストレージ９３は、一時的でない有形の記憶媒体である。

　上記の第１から第３の実施形態に係るガスタービン制御装置２０，２０ａ，２０ｂの各々は、一体のコンピュータとして実現されることに限定されない。例えば、第１及び第３の実施形態では、燃料流量算出部２０１，２０１ｂ、空気流量算出部２０２、燃料補正係数適用部２１４，２１４－１～２１４－３、空気補正係数適用部２２４、タービン入口温度算出部２０６、燃料配分比算出部２０７、及び弁開度算出部２０８の構成と、当該構成に含まれない機能部を含む構成とが、それぞれ異なるコンピュータに備えられており、２つのコンピュータが、通信回線によりデータの送受信を行うようにしてもよい。また、第２の実施形態では、燃料流量算出部２０１、空気流量算出部２０２、燃料補正係数算出部２１３ａ、燃料補正係数適用部２１４、空気補正係数算出部２２３ａ、空気補正係数適用部２２４、タービン入口温度算出部２０６、燃料配分比算出部２０７、及び弁開度算出部２０８の構成と、当該構成に含まれない機能部を含む構成とが、それぞれ異なるコンピュータに備えられており、２つのコンピュータが、通信回線によりデータの送受信を行うようにしてもよい。

＜付記＞
　各実施形態に記載のガスタービン制御装置（例えば、ガスタービン制御装置２０，２０ａ，２０ｂ）は、例えば以下のように把握される。

（１）第１の態様に係るガスタービン制御装置は、ガスタービン１０の燃焼器１２に燃料を供給する燃料供給系統４０に関する計測データ（例えば、弁開度Ｏ_１，Ｏ_２，Ｏ_３、流調弁前圧力Ｐ_１ＦＶ（１），Ｐ_１ＦＶ（２），Ｐ_１ＦＶ（３）、流調弁後圧力Ｐ_２ＦＶ（１），Ｐ_２ＦＶ（２），Ｐ_２ＦＶ（３）、燃料温度Ｔ_ｆ（１），Ｔ_ｆ（２），Ｔ_ｆ（３））に基づいて前記燃焼器に流入する燃料の量を示す燃料流量Ｇ_ｆｔを算出する燃料流量算出部（例えば、燃料流量算出部２０１，２０１ｂ）と、前記ガスタービンの圧縮機１１に関する計測データ（例えば、圧縮機インデックス差圧Ｐ_{ｉｎｄｅｘ}、圧縮機入口圧力Ｐ_１Ｃ、圧縮機入口温度Ｔ_１Ｃ）に基づいて前記圧縮機が吸入する空気の量を示す空気流量Ｇ_ａｔを算出する空気流量算出部２０２と、前記燃料流量Ｇ_ｆｔと、前記空気流量Ｇ_ａｔと、前記燃焼器に関する熱エネルギ収支に関する物理モデル式（例えば、燃料流量算出関数ｆ_１（・））とに基づいてタービン入口温度Ｔ_１Ｔを算出するタービン入口温度算出部２０６と、前記ガスタービンのエネルギ収支に関する計測データ（例えば、ガスタービン出力Ｇ_ｏｕｔ、排ガス圧力Ｐ_２Ｔ、排ガス温度Ｔ_２Ｔ、圧縮機入口圧力Ｐ_１Ｃ、圧縮機入口温度Ｔ_１Ｃ、圧縮機出口圧力Ｐ_２Ｃ、及び圧縮機出口温度Ｔ_２Ｃ）に基づいて前記ガスタービンのヒートバランスを満たす演算を行って演算結果を出力するヒートバランス演算部（例えば、ヒートバランス演算部２０３，２０３ｂ）と、前記ガスタービンの状態を検出する状態検出部２０５と、前記状態検出部が検出する前記ガスタービンの状態が静定状態の場合に前記ヒートバランス演算部が出力する前記演算結果に基づいて、前記燃料流量と前記空気流量のいずれか一方、または、両方を補正して前記タービン入口温度算出部に与える補正処理部（例えば、補正処理部２０４，２０４ａ，２０４ｂ）と、を備える。第１の態様、及び以下の各態様によれば、ヒートバランス演算部の演算結果に基づいて、燃料流量Ｇ_ｆｔと空気流量Ｇ_ａｔの補正が可能になるため、一意に定めることが難しいタービン効率を用いることなく補正を行うことができ、更に、ヒートバランス演算部は、繰り返し計測される計測データに基づいて演算を行うため、長期間にわたってタービン入口温度Ｔ_１Ｔの精度を維持することができる。

（２）第２の態様に係るガスタービン制御装置は、（１）のガスタービン制御装置であって、前記燃料流量算出部が算出する燃料流量を第１の燃料流量Ｇ_ｆｔとし、前記空気流量算出部が算出する空気流量を第１の空気流量Ｇ_ａｔとし、前記補正処理部は、前記状態検出部が検出する前記ガスタービンの状態が静定状態の場合に前記ヒートバランス演算部が出力する前記演算結果に含まれる第２の燃料流量Ｇ_ｆｓを取得し、取得した前記第２の燃料流量と、当該第２の燃料流量を取得した際の前記第１の燃料流量とに基づいて、前記燃料流量算出部が算出する前記第１の燃料流量を補正し、補正した前記第１の燃料流量を前記タービン入口温度算出部に与える燃料流量補正部（例えば、燃料流量補正部２１０，２１０ａ，２１０ｂ）と、前記状態検出部が検出する前記ガスタービンの状態が静定状態の場合に前記ヒートバランス演算部が出力する前記演算結果に含まれる第２の空気流量Ｇ_ａｓを取得し、取得した前記第２の空気流量と、当該第２の空気流量を取得した際の前記第１の空気流量とに基づいて、前記空気流量算出部が算出する前記第１の空気流量を補正し、補正した前記第１の空気流量を前記タービン入口温度算出部に与える空気流量補正部（例えば、空気流量補正部２２０，２２０ａ）とのいずれか一方の補正部、または、両方の補正部を備える。第２の態様によれば、ヒートバランス演算部の演算結果に含まれる精度の高い第２の燃料流量Ｇ_ｆｓによって第１の燃料流量Ｇ_ｆｔを補正することが可能になり、当該演算結果に含まれる精度の高い第２の空気流量Ｇ_ａｓによって第１の空気流量Ｇ_ａｔを補正することが可能になるので、精度の高いタービン入口温度Ｔ_１Ｔを算出することが可能になる。

（３）第３の態様に係るガスタービン制御装置は、（２）のガスタービン制御装置であって、前記燃料流量補正部は、前記状態検出部が検出する前記ガスタービンの状態が静定状態の場合に取得した中で最新の前記第２の燃料流量と、当該最新の第２の燃料流量を取得した際の前記第１の燃料流量とに基づいて燃料補正係数を算出する燃料補正係数算出部（例えば、燃料補正係数算出部２１３，２１３－１～２１３－３）と、前記燃料補正係数算出部が算出する前記燃料補正係数によって前記燃料流量算出部が算出する前記第１の燃料流量を補正する燃料補正係数適用部（例えば、燃料補正係数適用部２１４，２１４－１～２１４－３）と、を備える。第３の態様によれば、最新の第２の燃料流量Ｇ_ｆｓと、最新の第２の燃料流量Ｇ_ｆｓを取得した際の第１の燃料流量Ｇ_ｆｔとを用いて燃料補正係数Ｋ_ｆを算出し、算出した燃料補正係数Ｋ_ｆによって第１の燃料流量Ｇ_ｆｔを補正するため、精度の高いタービン入口温度Ｔ_１Ｔを算出することが可能になると共に、長期間にわたってタービン入口温度Ｔ_１Ｔの精度を維持することができる。

（４）第４の態様に係るガスタービン制御装置は、（２）のガスタービン制御装置であって、前記燃料流量補正部は、前記状態検出部が検出する前記ガスタービンの状態が静定状態の場合に取得した複数の前記第２の燃料流量と、当該複数の第２の燃料流量の各々を取得した際の前記第１の燃料流量とに基づいて燃料補正係数を算出する燃料補正係数算出部（例えば、燃料補正係数算出部２１３，２１３－１～２１３－３）と、前記燃料補正係数算出部が算出する前記燃料補正係数によって前記燃料流量算出部が算出する前記第１の燃料流量を補正する燃料補正係数適用部（例えば、燃料補正係数適用部２１４，２１４－１～２１４－３）と、を備える。第４の態様によれば、複数の第２の燃料流量Ｇ_ｆｓと、複数の第１の燃料流量Ｇ_ｆｔとに基づいて燃料補正係数Ｋ_ｆを算出し、算出した燃料補正係数Ｋ_ｆによって第１の燃料流量Ｇ_ｆｔを補正するため、信頼性の面で精度の高い燃料補正係数Ｋ_ｆが得られることになる。

（５）第５の態様に係るガスタービン制御装置は、（２）から（４）のいずれか１つのガスタービン制御装置であって、前記空気流量補正部は、前記状態検出部が検出する前記ガスタービンの状態が静定状態の場合に取得した中で最新の前記第２の空気流量と、当該最新の第２の空気流量を取得した際の前記第１の空気流量とに基づいて空気補正係数を算出する空気補正係数算出部２２３と、前記空気補正係数算出部が算出する前記空気補正係数によって前記空気流量算出部が算出する前記第１の空気流量を補正する空気補正係数適用部２２４と、を備える。第５の態様によれば、最新の第２の空気流量Ｇ_ａｓと、最新の第２の空気流量Ｇ_ａｓを取得した際の第１の空気流量Ｇ_ａｔとを用いて空気補正係数Ｋ_ａを算出し、算出した空気補正係数Ｋ_ａによって第１の空気流量Ｇ_ａｔを補正するため、精度の高いタービン入口温度Ｔ_１Ｔを算出することが可能になると共に、長期間にわたってタービン入口温度Ｔ_１Ｔの精度を維持することができる。

（６）第６の態様に係るガスタービン制御装置は、（２）から（４）のいずれか１つのガスタービン制御装置であって、前記空気流量補正部は、前記状態検出部が検出する前記ガスタービンの状態が静定状態の場合に取得した複数の前記第２の空気流量と、当該複数の第２の空気流量の各々を取得した際の前記第１の空気流量とに基づいて空気補正係数を算出する空気補正係数算出部２２３と、前記空気補正係数算出部が算出する前記空気補正係数によって前記空気流量算出部が算出する前記第１の空気流量を補正する空気補正係数適用部２２４と、を備える。第６の態様によれば、複数の第２の空気流量Ｇ_ａｓと、複数の第１の空気流量Ｇ_ａｔとに基づいて空気補正係数Ｋ_ａを算出し、算出した空気補正係数Ｋ_ａによって第１の空気流量Ｇ_ａｔを補正するため、信頼性の面で精度の高い空気補正係数Ｋ_ａが得られることになる。

（７）第７の態様に係るガスタービン制御装置は、（２），（５），（６）のいずれか１つのガスタービン制御装置であって、前記燃料流量補正部は、前記状態検出部が検出する前記ガスタービンの状態が静定状態の場合に取得した前記第２の燃料流量と、当該第２の燃料流量を取得した際の前記第１の燃料流量及び説明変数のデータ（例えば、燃焼器車室温度Ｔ_ＣＳ及び燃焼器車室圧力Ｐ_ｏｕｔ）とに基づいて所定の回帰分析（例えば、線形重回帰分析）を行って前記所定の回帰分析に対応する回帰式の係数を算出する燃料側回帰分析部２１５と、前記燃料側回帰分析部が算出する前記係数を適用した前記回帰式に対して、前記燃料流量算出部が前記第１の燃料流量を算出した際の前記説明変数のデータを代入して燃料補正係数を算出する燃料補正係数算出部（例えば、燃料補正係数算出部２１３ａ，２１３－１～２１３－３）と、前記燃料補正係数算出部が算出する前記燃料補正係数によって、前記回帰式に代入した前記説明変数のデータに対応する前記第１の燃料流量を補正する燃料補正係数適用部（例えば、燃料補正係数適用部２１４，２１４－１～２１４－３）と、を備える。第７の態様によれば、回帰分析を利用することにより、信頼性の面で精度の高い燃料補正係数Ｋ_ｆが得られることになる。

（８）第８の態様に係るガスタービン制御装置は、（２），（３），（４），（７）のいずれか１つのガスタービン制御装置であって、前記空気流量補正部は、前記状態検出部が検出する前記ガスタービンの状態が静定状態の場合に取得した前記第２の空気流量と、当該第２の空気流量を取得した際の前記第１の空気流量及び説明変数のデータとに基づいて所定の回帰分析を行って前記所定の回帰分析に対応する回帰式の係数を算出する空気側回帰分析部２２５と、前記空気側回帰分析部が算出する前記係数を適用した前記回帰式に対して、前記空気流量算出部が前記第１の空気流量を算出した際の前記説明変数のデータを代入して空気補正係数を算出する空気補正係数算出部２２３ａと、前記空気補正係数算出部が算出する前記空気補正係数によって、前記回帰式に代入した前記説明変数のデータに対応する前記第１の空気流量を補正する空気補正係数適用部２２４と、を備える。第８の態様によれば、回帰分析を利用することにより、信頼性の面で精度の高い空気補正係数Ｋ_ａが得られることになる。

（９）第９の態様に係るガスタービン制御装置は、（７）または（８）のガスタービン制御装置であって、前記説明変数のデータは、前記計測データと、前記ガスタービンの制御のための指令値のデータとを任意に組み合わせたデータである。第９の態様によれば、計測によって取得する計測データと、ガスタービン１０に与える指令値のデータとを反映した燃料補正係数Ｋ_ｆ、または、空気補正係数Ｋ_ａが得られることになる。

（１０）第１０の態様に係るガスタービン制御装置は、（２）から（９）のいずれか１つのガスタービン制御装置であって、前記補正処理部が、前記燃料流量補正部を備えている場合であって、かつ前記燃料供給系統が複数存在している場合、複数の前記燃料供給系統ごとの燃料配分比率（例えば、燃料配分比Ｄ_１，Ｄ_２，Ｄ_３）に基づいて、前記ヒートバランス演算部が出力する前記演算結果に含まれる前記第２の燃料流量から前記燃料供給系統ごとの前記第２の燃料流量を算出する燃料配分量算出部２０９を備え、前記燃料流量算出部２０１ｂは、前記燃料供給系統ごとの前記第１の燃料流量（例えば、燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｔ（２），Ｇ_ｆｔ（３））を算出し、前記燃料流量補正部２１０ｂは、前記状態検出部が検出する前記ガスタービンの状態が静定状態の場合に前記ヒートバランス演算部２０３ｂが出力する前記演算結果に含まれる第２の燃料流量に替えて、前記燃料配分量算出部が算出する前記燃料供給系統ごとの前記第２の燃料流量（例えば、燃料流量Ｇ_ｆｓ（１），Ｇ_ｆｔ（２），Ｇ_ｆｔ（３））を取得し、取得した前記第２の燃料流量の各々と、前記第２の燃料流量の各々を取得した際の前記燃料供給系統ごとの前記第１の燃料流量の各々とに基づいて、前記燃料流量算出部が算出する前記燃料供給系統ごとの前記第１の燃料流量を補正し、補正した前記燃料供給系統ごとの前記第１の燃料流量の合計値を前記タービン入口温度算出部に与える。第１０の態様によれば、燃料供給系統４０－１，４０－２，４０－３ごとの燃料流量Ｇ_ｆｔ（１），Ｇ_ｆｔ（２），Ｇ_ｆｔ（３）を個別に補正することが可能になる。

　上述した一態様によれば、タービン効率を用いずに、長期間にわたってタービン入口温度の精度を維持する。

２０　ガスタービン制御装置
２０１　燃料流量算出部
２０２　空気流量算出部
２０３　ヒートバランス演算部
２０４　補正処理部
２０５　状態検出部
２０６　タービン入口温度算出部
２０７　燃料配分比算出部
２０８　弁開度算出部
２１０　燃料流量補正部
２１１　燃料側管理部
２１２　燃料側記憶部
２１３　燃料補正係数算出部
２１４　燃料補正係数適用部
２２０　空気流量補正部
２２１　空気側管理部
２２２　空気側記憶部
２２３　空気補正係数算出部
２２４　空気補正係数適用部

Claims

　ガスタービンの燃焼器に燃料を供給する燃料供給系統に関する計測データに基づいて前記燃焼器に流入する燃料の量を示す燃料流量を算出する燃料流量算出部と、
　前記ガスタービンの圧縮機に関する計測データに基づいて前記圧縮機が吸入する空気の量を示す空気流量を算出する空気流量算出部と、
　前記燃料流量と、前記空気流量と、前記燃焼器に関する熱エネルギ収支に関する物理モデル式とに基づいてタービン入口温度を算出するタービン入口温度算出部と、
　前記ガスタービンのエネルギ収支に関する計測データに基づいて前記ガスタービンのヒートバランスを満たす演算を行って演算結果を出力するヒートバランス演算部と、
　前記ガスタービンの状態を検出する状態検出部と、
　前記状態検出部が検出する前記ガスタービンの状態が静定状態の場合に前記ヒートバランス演算部が出力する前記演算結果に基づいて、前記燃料流量と前記空気流量のいずれか一方、または、両方を補正して前記タービン入口温度算出部に与える補正処理部と、
　を備えるガスタービン制御装置。
　前記燃料流量算出部が算出する燃料流量を第１の燃料流量とし、
　前記空気流量算出部が算出する空気流量を第１の空気流量とし、
　前記補正処理部は、
　前記状態検出部が検出する前記ガスタービンの状態が静定状態の場合に前記ヒートバランス演算部が出力する前記演算結果に含まれる第２の燃料流量を取得し、取得した前記第２の燃料流量と、当該第２の燃料流量を取得した際の前記第１の燃料流量とに基づいて、前記燃料流量算出部が算出する前記第１の燃料流量を補正し、補正した前記第１の燃料流量を前記タービン入口温度算出部に与える燃料流量補正部と、
　前記状態検出部が検出する前記ガスタービンの状態が静定状態の場合に前記ヒートバランス演算部が出力する前記演算結果に含まれる第２の空気流量を取得し、取得した前記第２の空気流量と、当該第２の空気流量を取得した際の前記第１の空気流量とに基づいて、前記空気流量算出部が算出する前記第１の空気流量を補正し、補正した前記第１の空気流量を前記タービン入口温度算出部に与える空気流量補正部とのいずれか一方の補正部、または、両方の補正部
　を備える請求項１に記載のガスタービン制御装置。
　前記燃料流量補正部は、
　前記状態検出部が検出する前記ガスタービンの状態が静定状態の場合に取得した中で最新の前記第２の燃料流量と、当該最新の第２の燃料流量を取得した際の前記第１の燃料流量とに基づいて燃料補正係数を算出する燃料補正係数算出部と、
　前記燃料補正係数算出部が算出する前記燃料補正係数によって前記燃料流量算出部が算出する前記第１の燃料流量を補正する燃料補正係数適用部と、
　を備える請求項２に記載のガスタービン制御装置。
　前記燃料流量補正部は、
　前記状態検出部が検出する前記ガスタービンの状態が静定状態の場合に取得した複数の前記第２の燃料流量と、当該複数の第２の燃料流量の各々を取得した際の前記第１の燃料流量とに基づいて燃料補正係数を算出する燃料補正係数算出部と、
　前記燃料補正係数算出部が算出する前記燃料補正係数によって前記燃料流量算出部が算出する前記第１の燃料流量を補正する燃料補正係数適用部と、
　を備える請求項２に記載のガスタービン制御装置。
　前記空気流量補正部は、
　前記状態検出部が検出する前記ガスタービンの状態が静定状態の場合に取得した中で最新の前記第２の空気流量と、当該最新の第２の空気流量を取得した際の前記第１の空気流量とに基づいて空気補正係数を算出する空気補正係数算出部と、
　前記空気補正係数算出部が算出する前記空気補正係数によって前記空気流量算出部が算出する前記第１の空気流量を補正する空気補正係数適用部と、
　を備える請求項２に記載のガスタービン制御装置。
　前記空気流量補正部は、
　前記状態検出部が検出する前記ガスタービンの状態が静定状態の場合に取得した複数の前記第２の空気流量と、当該複数の第２の空気流量の各々を取得した際の前記第１の空気流量とに基づいて空気補正係数を算出する空気補正係数算出部と、
　前記空気補正係数算出部が算出する前記空気補正係数によって前記空気流量算出部が算出する前記第１の空気流量を補正する空気補正係数適用部と、
　を備える請求項２に記載のガスタービン制御装置。
　前記燃料流量補正部は、
　前記状態検出部が検出する前記ガスタービンの状態が静定状態の場合に取得した前記第２の燃料流量と、当該第２の燃料流量を取得した際の前記第１の燃料流量及び説明変数のデータとに基づいて所定の回帰分析を行って前記所定の回帰分析に対応する回帰式の係数を算出する燃料側回帰分析部と、
　前記燃料側回帰分析部が算出する前記係数を適用した前記回帰式に対して、前記燃料流量算出部が前記第１の燃料流量を算出した際の前記説明変数のデータを代入して燃料補正係数を算出する燃料補正係数算出部と、
　前記燃料補正係数算出部が算出する前記燃料補正係数によって、前記回帰式に代入した前記説明変数のデータに対応する前記第１の燃料流量を補正する燃料補正係数適用部と、
　を備える請求項２に記載のガスタービン制御装置。
　前記空気流量補正部は、
　前記状態検出部が検出する前記ガスタービンの状態が静定状態の場合に取得した前記第２の空気流量と、当該第２の空気流量を取得した際の前記第１の空気流量及び説明変数のデータとに基づいて所定の回帰分析を行って前記所定の回帰分析に対応する回帰式の係数を算出する空気側回帰分析部と、
　前記空気側回帰分析部が算出する前記係数を適用した前記回帰式に対して、前記空気流量算出部が前記第１の空気流量を算出した際の前記説明変数のデータを代入して空気補正係数を算出する空気補正係数算出部と、
　前記空気補正係数算出部が算出する前記空気補正係数によって、前記回帰式に代入した前記説明変数のデータに対応する前記第１の空気流量を補正する空気補正係数適用部と、
　を備える請求項２に記載のガスタービン制御装置。
　前記説明変数のデータは、前記計測データと、前記ガスタービンの制御のための指令値のデータとを任意に組み合わせたデータである、
　請求項７または請求項８に記載のガスタービン制御装置。
　前記補正処理部が、前記燃料流量補正部を備えている場合であって、かつ前記燃料供給系統が複数存在している場合、
　複数の前記燃料供給系統ごとの燃料配分比率に基づいて、前記ヒートバランス演算部が出力する前記演算結果に含まれる前記第２の燃料流量から前記燃料供給系統ごとの前記第２の燃料流量を算出する燃料配分量算出部を備え、
　前記燃料流量算出部は、
　前記燃料供給系統ごとの前記第１の燃料流量を算出し、
　前記燃料流量補正部は、
　前記状態検出部が検出する前記ガスタービンの状態が静定状態の場合に前記ヒートバランス演算部が出力する前記演算結果に含まれる第２の燃料流量に替えて、前記燃料配分量算出部が算出する前記燃料供給系統ごとの前記第２の燃料流量を取得し、取得した前記第２の燃料流量の各々と、前記第２の燃料流量の各々を取得した際の前記燃料供給系統ごとの前記第１の燃料流量の各々とに基づいて、前記燃料流量算出部が算出する前記燃料供給系統ごとの前記第１の燃料流量を補正し、補正した前記燃料供給系統ごとの前記第１の燃料流量の合計値を前記タービン入口温度算出部に与える、
　請求項２から請求項８のいずれか一項に記載のガスタービン制御装置。
　ガスタービンの燃焼器に燃料を供給する燃料供給系統に関する計測データに基づいて前記燃焼器に流入する燃料の量を示す燃料流量を算出するステップと、
　前記ガスタービンの圧縮機に関する計測データに基づいて前記圧縮機が吸入する空気の量を示す空気流量を算出するステップと、
　前記燃料流量と、前記空気流量と、前記燃焼器に関する熱エネルギ収支に関する物理モデル式とに基づいてタービン入口温度を算出するステップと、
　前記ガスタービンのエネルギ収支に関する計測データに基づいて前記ガスタービンのヒートバランスを満たす演算を行って演算結果を出力するステップと、
　前記ガスタービンの状態を検出するステップと、
　前記状態を検出するステップによって検出された前記ガスタービンの状態が静定状態の場合に前記ヒートバランスを満たす演算を行うステップにおいて出力された前記演算結果に基づいて、前記燃料流量と前記空気流量のいずれか一方、または、両方を補正して前記タービン入口温度を算出するステップに与えるステップと、
　を含むガスタービン制御方法。
　コンピュータを、
　ガスタービンの燃焼器に燃料を供給する燃料供給系統に関する計測データに基づいて前記燃焼器に流入する燃料の量を示す燃料流量を算出する燃料流量算出手段、
　前記ガスタービンの圧縮機に関する計測データに基づいて前記圧縮機が吸入する空気の量を示す空気流量を算出する空気流量算出手段、
　前記燃料流量と、前記空気流量と、前記燃焼器に関する熱エネルギ収支に関する物理モデル式とに基づいてタービン入口温度を算出するタービン入口温度算出手段、
　前記ガスタービンのエネルギ収支に関する計測データに基づいて前記ガスタービンのヒートバランスを満たす演算を行って演算結果を出力するヒートバランス演算手段、
　前記ガスタービンの状態を検出する状態検出手段、
　前記状態検出手段が検出する前記ガスタービンの状態が静定状態の場合に前記ヒートバランス演算手段が出力する前記演算結果に基づいて、前記燃料流量と前記空気流量のいずれか一方、または、両方を補正して前記タービン入口温度算出手段に与える補正処理手段、
　として機能させるためのプログラム。