JP2017504753A - ガスタービンの動作パラメータの測定誤差の分析方法、及び、制御装置 - Google Patents

Abstract

ガスタービンの動作パラメータの測定誤差を検出する分析方法は、動作パラメータを、それぞれ予め定められた分析範囲を有する複数の関数の変数の少なくとも一部としたとき、当該動作パラメータの複数の測定値を取得するステップと、動作パラメータの測定値を、動作パラメータを含む変数を有する関数へ投入して、関数の結果を計算するステップと、関数の計算結果を、対応する関数の分析範囲と比較するステップと、計算結果が対応する関数の予め定められた分析範囲に入らない関数の変数の少なくとも一部である動作パラメータの１つもしくは複数の測定値を求めるステップと、求められた１つもしくは複数の測定値のなかで、計算結果が対応する関数の予め定められた分析範囲に入る関数の変数ではない動作パラメータの測定値が測定誤差の最尤度を有することを検出するステップとを含む。

Description

本発明は、ガスタービン分野、より具体的にはガスタービンの制御分野に関する。

関連技術
ガスタービンは、圧縮器セクション、燃焼器セクション及びタービンセクションを含む。通常はこれら３つのセクションの複数の位置に複数のセンサが配置されている。これらのセンサは、ガスタービンの制御装置がガスタービンの動作状態を監視してその動作パラメータを制御できるよう、動作圧、動作温度、圧縮器内のガス流量などの動作パラメータをこの制御装置に供給する。例えば、或るセンサは、タービンセクションの圧力、温度及び流量などの動作パラメータの値を検出するために、このセクションの燃料装置、空気装置又はガス装置に配置される必要がある。

センサが検出した測定値は、通常、ガスタービンの所定の位置での動作パラメータの実際値を反映している。しかし、センサに誤りが発生した場合、測定値は実際値から偏差する。この状況では、測定値は当該位置での動作状態を正確に反映しない。ガスタービンの制御装置が実際値を正確に反映しない測定値を受信すると、制御装置が検出したガスタービンの動作状態が不正確となり、このことはガスタービンの動作不安定、ひいてはバルブトリップを生じさせかねない。

概要
本発明の実施形態により、実際値から偏差した誤差値を求めるために、ガスタービンの動作パラメータの測定誤差を検出する分析方法が提供される。

当該分析方法は、動作パラメータを、それぞれ予め定められた分析範囲を有する複数の関数の変数の少なくとも一部としたとき、当該動作パラメータの複数の測定値を取得するステップと、当該動作パラメータの測定値を、当該動作パラメータを含む変数を有する関数へ投入して関数の結果を計算するステップと、関数の計算結果を対応する関数の分析範囲と比較するステップと、計算結果が対応する関数の予め定められた分析範囲に入らない関数の変数の少なくとも一部である動作パラメータの１つもしくは複数の測定値を求めるステップと、求められた１つもしくは複数の測定値のなかで、計算結果が対応する関数の予め定められた分析範囲に入る関数の変数ではない動作パラメータの測定値が測定誤差の最尤度を有することを検出するステップとを含む。

上述した関数の分析範囲は、当該関数の閾値範囲であってよく、又は、閾値範囲に許容可能な誤差範囲を加えたものであってもよい。「動作パラメータに関連する関数」とは、動作パラメータを含む変数を有する関数のことである。

一実施形態では、分析方法はさらに、測定誤差の最尤度を有すると検出された測定値の計算値を、当該測定値を含むいずれかの関数の分析範囲の中間値を用いて取得するステップと、当該測定値に代えて、補正値としての計算値を用いるステップとを含む。

別の実施形態では、分析方法はさらに、測定誤差の最尤度を有すると検出された測定値の計算値を、当該測定値を含む関数の分析範囲の中間値を用いて取得し、当該測定値を含む複数の関数（すなわち当該測定値を含む関数の変数）から計算値の平均を取得するステップと、当該測定値に代えて、補正値としての計算値を用いるステップとを含む。

別の実施形態では、ガスタービンの動作パラメータは、ガスタービン内の流体の温度、圧力、流量及び圧力差を含む。動作パラメータの複数の測定値は、ガスタービンのオイル装置、空気装置及びガス装置の動作パラメータの測定値を含む。

別の実施形態では、ガスタービン用コントローラが、ガスタービンの動作パラメータの測定誤差を検出するための分析方法を実行するようにプログラミングされる。当該分析方法は、動作パラメータを、それぞれ予め定められた分析範囲を有する複数の関数の変数の少なくとも一部としたとき、当該動作パラメータの複数の測定値を取得するステップと、当該動作パラメータの測定値を、当該動作パラメータを含む変数を有する関数へ投入して関数の結果を計算するステップと、関数の計算結果を対応する関数の分析範囲と比較するステップと、計算結果が対応する関数の予め定められた分析範囲に入らない関数の変数の少なくとも一部である動作パラメータの１つもしくは複数の測定値を求めるステップと、求められた１つもしくは複数の測定値のなかで、計算結果が対応する関数の予め定められた分析範囲に入る関数の変数ではない動作パラメータの測定値が測定誤差の最尤度を有することを検出するステップとを含む。

別の実施形態では、ガスタービンの燃料ガス装置が、燃料ガスをガスタービンへ輸送する燃料管と、燃料管内の燃料の物理的な動作パラメータの複数の値を測定する複数のセンサと、上述した分析方法を実行するようにプログラミングされたプログラマブルコントローラとを含む。

別の実施形態では、コントローラが、どのセンサで測定誤差が発生したかを報知する。

以下の詳細な説明から、開示内容がより良く理解されるであろう。ただし、詳細な説明は例示のためのものにすぎず、よって開示内容を限定しない。なお、本発明を、単なる典型例として、添付図を参照しながら説明する。各図は本発明の範囲を限定しない。
ガスタービンでの測定値の分析方法における分析プロセスを説明するための図である。 ガスタービンのガス供給装置の構造を示す概略図である。 図２の測定値の分析方法を説明するための図である。

参照番号
１０ 主管路、 １２ トリップ弁、 １４，１８，２１，２６，３１，３６ 圧力センサ、 １６，２８，３８ 温度センサ、 ２０ 第１の分岐管路、 ２２ 第１の制御弁、 ２４，３４ 流量センサ、 ２５，３５ 圧力差センサ、 ３０ 第２の分岐管路、 ３２ 第２の制御弁、 ３３ コントローラ

詳細な説明
本発明の技術的特徴、目的及び効果をより良く理解してもらうために、添付図を参照しながら、本発明の特定の幾つかの構成形態を説明する。図中、同じ要素には同じ参照番号を付してある。

本明細書で「典型的」とは、「実例、事例もしくは詳細例として用いられる」ことを表し、本明細書で「典型的」として説明される図又は構成形態は、より好ましいもしくはより有利な技術手段と解されるべきでない。

図示を簡単化するため、本発明に関連する要素のみを各図に典型的に示しており、各要素は模式的なものであって、製品において使用される実際の要素を表していないこともある。さらに、図示を簡単化しかつわかりやすくするため、幾つかの図にわたって同様の構造もしくは機能を有する要素がある場合、そのうち１つのみを典型的に表示又は標示している。

本明細書での「１つの」とは「唯一の」ケースのみならず「２つ以上」のケースも表しうる。本明細書での「第１の」「第２の」及び同様の記載は、相対的な重要度や序列を表すというよりもむしろ単に相互の区別のために用いている。

本明細書での“ｆ１”，“ｆ２”，“ｆ３”及び“ｆ１０”，“ｆ１１”，“ｆ１２”，“ｆ１３”，“ｆ１４”，“ｆ１５”，“ｆ１６”，“ｆ１７”，“ｆ１８”は、それぞれ異なる関数を表すために用いている。

図１を用いて、ガスタービンの測定値の分析方法における分析プロセスを説明する。ガスタービンの測定値は、１つの第１の測定値Ｓ１と１つの第２の測定値Ｓ２と１つの第３の測定値Ｓ３とを含み、これらの測定値は、ガスタービンの種々異なる組み込み位置に設けられた複数のセンサでの測定によって相互に別個に取得されたものである。また、これらの測定値間には関数関係が存在する。すなわち、例えば、第１の分析関数ｙ１＝ｆ１（Ｓ１，Ｓ２）は変数として第１の測定値Ｓ１及び第２の測定値Ｓ２を含み、第２の分析関数ｙ２＝ｆ２（Ｓ１，Ｓ３）は変数として第１の測定値Ｓ１及び第３の測定値Ｓ３を含み、第３の分析関数ｙ３＝ｆ３（Ｓ２，Ｓ３）は変数として第２の測定値Ｓ２及び第３の測定値Ｓ３を含む。

ガスタービンが通常動作を行い、各測定値が正確にガスタービンの実際の動作状態を反映し、その期間においてさらに第１の測定値Ｓ１及び第２の測定値Ｓ２及び第３の測定値Ｓ３の測定精度を考慮する必要がある場合、第１の分析関数の範囲を１つの第１の分析範囲Ｙ１と定義し、第２の分析関数の範囲を１つの第２の分析範囲Ｙ２と定義し、第３の分析関数の範囲を１つの第３の分析範囲Ｙ３と定義する。

ステップＳ１０では、第１の測定値Ｓ１及び第２の測定値Ｓ２及び第３の測定値Ｓ３が取得される。次に、プロセスはステップＳ２０へ移行する。

ステップＳ２０では、第１の測定値Ｓ１及び第２の測定値Ｓ２が第１の分析関数ｙ１＝ｆ１（Ｓ１，Ｓ２）へ代入されて計算により第１の分析値ｙ１が取得され、第１の測定値Ｓ１及び第３の測定値Ｓ３が第２の分析関数ｙ２＝ｆ２（Ｓ１，Ｓ３）へ代入されて計算により第２の分析値ｙ２が取得され、第２の測定値Ｓ２及び第３の測定値Ｓ３が第３の分析関数ｙ３＝ｆ３（Ｓ２，Ｓ３）へ代入されて計算により第３の分析値ｙ３が取得される。ついで、プロセスはステップＳ３０へ移行する。

ステップＳ３０では、第１の分析値ｙ１が第１の分析範囲Ｙ１と比較され、第２の分析値ｙ２が第２の分析範囲Ｙ２と比較され、第３の分析値ｙ３が第３の分析範囲Ｙ３と比較される。その後、プロセスは並列的にステップＳ４０，Ｓ４２，Ｓ４４へ移行する。

ステップＳ４０では、第１の分析値ｙ１が第１の分析範囲Ｙ１に入らず、かつ、第２の分析値ｙ２が第２の分析範囲Ｙ２に入らず、かつ、第３の分析値ｙ３が第３の分析範囲Ｙ３に入るか否かが判別される。判別結果がノーである場合、分析プロセスは終了する。判別結果がイエスである場合、第１の分析関数及び第２の分析関数の双方に含まれるが第３の分析関数には含まれない第１の測定値Ｓ１が測定誤差の最尤度を有することが検出される。つまり、最も発生確率が高いのは、第１の測定値がセンサの組み込み位置での実際値から偏差することである。こうして分析プロセスは終了する。

ステップＳ４２では、第１の分析値ｙ１が第１の分析範囲Ｙ１に入らず、かつ、第３の分析値ｙ３が第３の分析範囲Ｙ３に入らず、かつ、第２の分析値ｙ２が第２の分析範囲Ｙ２に入るか否かが判別される。判別結果がノーである場合、分析プロセスは終了する。判別結果がイエスである場合、第１の分析関数及び第３の分析関数の双方に共通の変数であって第２の分析関数の変数でない第２の測定値Ｓ２が測定誤差の最尤度を有することが検出される。つまり、最も発生確率が高いのは、第２の測定値がセンサの組み込み位置での動作パラメータの実際値から偏差することである。こうして分析プロセスは終了する。

ステップＳ４４では、第２の分析値ｙ２が第２の分析範囲Ｙ２に入らず、かつ、第３の分析値ｙ３が第３の分析範囲Ｙ３に入らず、かつ、第１の分析値ｙ１が第１の分析範囲Ｙ１に入るか否かが判別される。判別結果がノーである場合、分析プロセスは終了する。判別結果がイエスである場合、第２の分析関数及び第３の分析関数の双方の変数であって第１の分析関数の変数でない第３の測定値Ｓ３が測定誤差の最尤度を有することが検出される。つまり、最も発生確率が高いのは、第３の測定値がセンサの組み込み位置での動作パラメータの実際値から偏差することである。こうして分析プロセスは終了する。

ガスタービンにおける測定値の分析方法では、各測定値が対応する分析関数に代入されて計算により分析関数値が取得され、分析関数値がここでの分析関数に対応する分析範囲に入らない分析関数が探索される。実際値から偏差する確率が最も高い測定値をオンライン又はリアルタイムで発見及び分析するために、全ての関数に共通に含まれる測定値が探索され、これにより、ガスタービンの制御装置が誤った測定値を使用して、ガスタービン全体の動作が不安定となったり、最悪の場合にはバルブトリップが発生したりすることが回避される。

ガスタービンの測定値の分析方法の典型的な構成形態では、測定誤差が第１の測定値Ｓ１に発生したことが検出された場合、第２の測定値Ｓ２と第１の分析範囲の１つの値とが第１の分析関数へ代入され、第１の測定値の補正値として使用可能な１つの第１の計算値が取得される。当該第１の計算値は、センサの組み込み位置でのガスタービンの実際の動作状態（すなわち実際値）を反映可能である。好ましくは、第１の分析範囲の中間値が第１の分析関数に代入される。第２の測定値Ｓ２は正しいので、第１の分析関数においては第１の測定値Ｓ１のみが未知である。よって、第１の計算値を第１の分析関数から計算によって取得できる。さらに、上述した方式で複数の関数から取得可能な複数の計算値の平均を計算によって取得することもでき、測定誤差の発生した対応する測定値に代えてこの平均を用いることができる。

以下では、第１の測定値Ｓ１に測定誤差が発生したことが検出された場合であって、第１の測定値Ｓ１が第１の分析関数ｙ１＝ｆ１（Ｓ１，Ｓ２）と第２の分析関数ｙ２＝ｆ２（Ｓ１，Ｓ３）との変数である例を用いて説明する。まず、第１の分析範囲の中間値が第１の分析関数に代入される。第２の測定値Ｓ２は既知でありかつ正しいので、第１の計算値を計算によって取得できる。次に、第２の分析範囲の中間値が第２の分析関数に代入される。第３の測定値Ｓ３は既知でありかつ正しいので、第２の計算値も計算によって取得できる。第１の計算値と第２の計算値との平均値が計算され、この平均値は第１の測定値のセンサの組み込み位置での動作パラメータの実際の動作状態を反映可能である。当該平均値は、第１の測定値Ｓ１の補正値として使用される。

測定誤差が第２の測定値Ｓ２に発生したことが検出された場合、第３の測定値Ｓ３と第３の分析範囲の１つの値とが第３の分析関数へ代入され、１つの第２の計算値が計算によって取得される。この第２の計算値は第２の測定値の補正値として使用可能である。当該第２の計算値は、第２の測定値のセンサの組み込み位置での実際の動作状態を反映可能である。計算プロセスは第１の計算値の計算プロセスと同様であるので、これ以上は説明しない。

測定誤差が第３の測定値Ｓ３に発生したことが検出された場合、第２の測定値Ｓ２と第２の分析範囲の１つの値とが第２の分析関数へ代入され、１つの第３の計算値が計算によって取得される。この第３の計算値は第３の測定値の補正値として使用可能である。当該第３の計算値は、第３の測定値のセンサの組み込み位置での実際の動作状態を反映可能である。計算プロセスは第１の計算値の計算プロセスと同様であるので、これ以上は説明しない。

図２は、ガスタービンの測定値の分析方法を説明するための、ガスタービンのガス供給装置の構造の概略図である。ただし、ガスタービンの測定値の分析方法は、図２のガス供給装置には限定されず、ガスタービンの他の装置、例えば燃料装置、ガス装置又は空気装置などにも適用可能である。図２に示されているように、ガスタービンのガス供給装置は、１つの主管路１０と１つの第１の分岐管路２０と１つの第２の分岐管路３０とを含む。

１つのトリップ弁１２が第１の主管路１０に配置されている。第１の主管路１０では、ガスの移動経路に沿って、１つの圧力センサ１４と１つの温度センサ１６とがトリップ弁１２の上流に配置されている。トリップ弁１２の下流には、１つの圧力センサ１８が配置されている。圧力センサ１４は、この圧力センサ１４の組み込み位置での圧力を表す測定値Ｐ０を出力し、温度センサ１６はこの温度センサ１６の組み込み位置での温度を表す測定値Ｔ０を出力し、圧力センサ１８はこの圧力センサ１８の組み込み位置での圧力を表す測定値Ｐ１を出力する。

１つの第１の制御弁２２が第１の分岐管路２０に配置されている。第１の分岐管路２０では、ガスの移動経路に沿って、１つの流量センサ２４が第１の制御弁２２の下流に配置されている。流量センサ２４と第１の制御弁２２との間には１つの圧力センサ２１が配置されている。流量センサ２４の下流には１つの圧力センサ２６と１つの温度センサ２８とが配置されている。さらに、第１の制御弁２２には、その入口と出口との間のガス圧力差を測定可能な１つの圧力差センサ２５が配置されている。流量センサ２４はこの流量センサ２４の組み込み位置でのガス流量を表す測定値ｍ＿ＰＭを出力でき、圧力センサ２１はこの圧力センサ２１の組み込み位置での圧力を表す測定値Ｐ２ＰＭを出力でき、圧力センサ２６はこの圧力センサ２６の組み込み位置での圧力を表す測定値Ｐ３ＰＭを出力でき、温度センサ２８はこの温度センサ２８の組み込み位置での温度を表す測定値Ｔ３ＰＭを出力でき、圧力差センサ２５はこの圧力差センサ２５の組み込み位置での圧力差を表す測定値ｄｐ＿ＣＶＰＭを出力できる。各パラメータの添え字ＰＭは、第１の分岐管路の各センサの測定値を表す。

１つの第２の制御弁３２が第２の分岐管路３０に配置されている。第２の分岐管路３０では、ガスの移動経路に沿って、１つの流量センサ３４が第２の制御弁３２の下流に配置されている。流量センサ３４と第２の制御弁３２との間には１つの圧力センサ３１が配置されている。流量センサ３４の下流には１つの圧力センサ３６と１つの温度センサ３８とが配置されている。さらに、第２の制御弁３２には、その入口と出口との間のガス圧力差を測定可能な１つの圧力差センサ３５が配置されている。

流量センサ３４はこの流量センサ３４の組み込み位置でのガス流量を表す測定値ｍ＿ＰＧを出力でき、圧力センサ３１はこの圧力センサ３１の組み込み位置での圧力を表す測定値Ｐ２ＰＧを出力でき、圧力センサ３６はこの圧力センサ３６の組み込み位置での圧力を表す測定値Ｐ３ＰＧを出力でき、温度センサ３８はこの温度センサ３８の組み込み位置での温度を表す測定値Ｔ３ＰＧを出力でき、圧力差センサ３５はこの圧力差センサ３５の組み込み位置での圧力差を表す測定値ｄｐ＿ＣＶＰＧを出力できる。各パラメータの添え字ＰＧは、第２の分岐管路の各センサの測定値を表す。

第１の分析関数ｙ１＝ｆ１０（Ｐ０，Ｐ１，Ｔ０）は測定値Ｐ０，Ｐ１，Ｔ０に関連しており、すなわち、Ｐ０，Ｐ１，Ｔ０は第１の分析関数の変数である。第１の分析関数の特異的表現は流体力学を反映した既存の計算式から導出可能である。ガスタービンが通常動作を行い、測定値Ｐ０，Ｐ１，Ｔ０が実際の動作状態を正確に反映している場合、測定値Ｐ０，Ｐ１，Ｔ０の測定精度を考慮しつつこれらの測定値Ｐ０，Ｐ１，Ｔ０を第１の分析関数に代入する計算によって得られるｙ１は、第１の分析範囲Ｙ１すなわち第１の分析関数の範囲に入る。

第２の分析関数ｙ２＝ｆ１１（Ｐ１，Ｐ２ＰＭ，ｄｐ＿ＣＶＰＭ）は測定値Ｐ１，Ｐ２ＰＭ，ｄｐ＿ＣＶＰＭに関連しており、第２の分析関数の特異的表現は流体力学に関する既存の計算式から導出可能である。ガスタービンが通常動作を行い、測定値Ｐ１，Ｐ２ＰＭ，ｄｐ＿ＣＶＰＭが実際の動作状態を正確に反映している場合、測定値Ｐ１，Ｐ２ＰＭ，ｄｐ＿ＣＶＰＭの測定精度を考慮しつつこれらの測定値Ｐ１，Ｐ２ＰＭ，ｄｐ＿ＣＶＰＭを第２の分析関数に代入する計算によって得られるｙ２は、第２の分析範囲Ｙ２すなわち第２の分析関数の範囲に入る。

第３の分析関数ｙ３＝ｆ１２（Ｐ１，Ｐ２ＰＧ，ｄｐ＿ＣＶＰＧ）は測定値Ｐ１，Ｐ２ＰＧ，ｄｐ＿ＣＶＰＧに関連しており、第３の分析関数の特異的表現は流体力学に関する既存の計算式から導出可能である。ガスタービンが通常動作を行い、測定値Ｐ１，Ｐ２ＰＧ，ｄｐ＿ＣＶＰＧが実際の動作状態を正確に反映している場合、測定値Ｐ１，Ｐ２ＰＧ，ｄｐ＿ＣＶＰＧの測定精度を考慮しつつこれらの測定値Ｐ１，Ｐ２ＰＧ，ｄｐ＿ＣＶＰＧを第３の分析関数に代入する計算によって得られるｙ３は、第３の分析範囲Ｙ３すなわち第３の分析関数の範囲に入る。

第４の分析関数ｙ４＝ｆ１３（Ｐ２ＰＭ，Ｐ３ＰＭ，Ｔ３ＰＭ，ｍ＿ＰＭ）は測定値Ｐ２ＰＭ，Ｐ３ＰＭ，Ｔ３ＰＭ，ｍ＿ＰＭに関連しており、第４の分析関数の特異的表現は流体力学に関する既存の計算式から導出可能である。ガスタービンが通常動作を行い、測定値Ｐ２ＰＭ，Ｐ３ＰＭ，Ｔ３ＰＭ，ｍ＿ＰＭが実際の動作状態を正確に反映している場合、測定値Ｐ２ＰＭ，Ｐ３ＰＭ，Ｔ３ＰＭ，ｍ＿ＰＭの測定精度を考慮しつつこれらの測定値Ｐ２ＰＭ，Ｐ３ＰＭ，Ｔ３ＰＭ，ｍ＿ＰＭを第４の分析関数に代入する計算によって得られるｙ４は、第４の分析範囲Ｙ４すなわち第４の分析関数の範囲に入る。

第５の分析関数ｙ５＝ｆ１４（Ｐ２ＰＧ，Ｐ３ＰＧ，Ｔ３ＰＧ，ｍ＿ＰＧ）は測定値Ｐ２ＰＧ，Ｐ３ＰＧ，Ｔ３ＰＧ，ｍ＿ＰＧに関連しており、第５の分析関数の特異的表現は流体力学に関する既存の計算式から導出可能である。ガスタービンが通常動作を行い、測定値Ｐ２ＰＧ，Ｐ３ＰＧ，Ｔ３ＰＧ，ｍ＿ＰＧが実際の動作状態を正確に反映している場合、測定値Ｐ２ＰＧ，Ｐ３ＰＧ，Ｔ３ＰＧ，ｍ＿ＰＧの測定精度を考慮しつつこれらの測定値Ｐ２ＰＧ，Ｐ３ＰＧ，Ｔ３ＰＧ，ｍ＿ＰＧを第５の分析関数に代入する計算によって得られるｙ５は、第５の分析範囲Ｙ５すなわち第５の分析関数の範囲に入る。

第６の分析関数ｙ６＝ｆ１５（Ｔ０，Ｔ３ＰＭ，Ｔ３ＰＧ）は測定値Ｔ０，Ｔ３ＰＭ，Ｔ３ＰＧに関連しており、第６の分析関数の特異的表現は流体力学に関する既存の計算式から導出可能である。ガスタービンが通常動作を行い、測定値Ｔ０，Ｔ３ＰＭ，Ｔ３ＰＧが実際の動作状態を正確に反映している場合、測定値Ｔ０，Ｔ３ＰＭ，Ｔ３ＰＧの測定精度を考慮しつつこれらの測定値Ｔ０，Ｔ３ＰＭ，Ｔ３ＰＧを第６の分析関数に代入する計算によって得られるｙ６は、第６の分析範囲Ｙ６すなわち第６の分析関数の範囲に入る。

第７の分析関数ｙ７＝ｆ１６（Ｐ１，ｄｐ＿ＣＶＰＭ，Ｔ０，ｍ＿ＰＭ）は測定値Ｐ１，ｄｐ＿ＣＶＰＭ，Ｔ０，ｍ＿ＰＭに関連しており、第７の分析関数の特異的表現は流体力学に関する既存の計算式から導出可能である。ガスタービンが通常動作を行い、測定値Ｐ１，ｄｐ＿ＣＶＰＭ，Ｔ０，ｍ＿ＰＭが実際の動作状態を正確に反映している場合、測定値Ｐ１，ｄｐ＿ＣＶＰＭ，Ｔ０，ｍ＿ＰＭの測定精度を考慮しつつこれらの測定値Ｐ１，ｄｐ＿ＣＶＰＭ，Ｔ０，ｍ＿ＰＭを第７の分析関数に代入する計算によって得られるｙ７は、第７の分析範囲Ｙ７すなわち第７の分析関数の範囲に入る。

第８の分析関数ｙ８＝ｆ１７（Ｐ１，ｄｐ＿ＣＶＰＧ，Ｔ０，ｍ＿ＰＧ）は測定値Ｐ１，ｄｐ＿ＣＶＰＧ，Ｔ０，ｍ＿ＰＧに関連しており、第８の分析関数の特異的表現は流体力学に関する既存の計算式から導出可能である。ガスタービンが通常動作を行い、測定値Ｐ１，ｄｐ＿ＣＶＰＧ，Ｔ０，ｍ＿ＰＧが実際の動作状態を正確に反映している場合、測定値Ｐ１，ｄｐ＿ＣＶＰＧ，Ｔ０，ｍ＿ＰＧの測定精度を考慮しつつこれらの測定値Ｐ１，ｄｐ＿ＣＶＰＧ，Ｔ０，ｍ＿ＰＧを第８の分析関数に代入する計算によって得られるｙ８は、第８の分析範囲Ｙ８すなわち第８の分析関数の範囲に入る。

第９の分析関数ｙ９＝ｆ１８（Ｐ０，Ｐ１，Ｔ０，ｍ＿ＰＧ，ｍ＿ＰＭ）は測定値Ｐ０，Ｐ１，Ｔ０，ｍ＿ＰＧ，ｍ＿ＰＭに関連しており、第９の分析関数の特異的表現は流体力学に関する既存の計算式から導出可能である。ガスタービンが通常動作を行い、測定値Ｐ０，Ｐ１，Ｔ０，ｍ＿ＰＧ，ｍ＿ＰＭが実際の動作状態を正確に反映している場合、測定値Ｐ０，Ｐ１，Ｔ０，ｍ＿ＰＧ，ｍ＿ＰＭの測定精度を考慮しつつこれらの測定値Ｐ０，Ｐ１，Ｔ０，ｍ＿ＰＧ，ｍ＿ＰＭを第９の分析関数に代入する計算によって得られるｙ９は、第９の分析範囲Ｙ９すなわち第９の分析関数の範囲に入る。

図３を用いて、図２の測定値の測定誤差の分析方法を説明する。ここで、図３は、図３ａ，図３ｂ，図３ｃに分割されて示されている。図３に示されているように、ステップＳ１０で、測定値Ｐ０，Ｐ１，Ｐ２ＰＭ，Ｐ２ＰＧ，Ｐ３ＰＭ，Ｐ３ＰＧ，ｄｐ＿ＣＶＰＧ，ｄｐ＿ＣＶＰＭ，ｍ＿ＰＧ，ｍ＿ＰＭ，Ｔ０，Ｔ３ＰＭ，Ｔ３ＰＧが取得される。プロセスはステップＳ２０へ移行する。

ステップＳ２０では、取得された測定値Ｐ０，Ｐ１，Ｔ０が第１の分析関数ｙ１＝ｆ１０（Ｐ０，Ｐ１，Ｔ０）に代入され、計算により第１の分析値ｙ１が取得される。取得された測定値Ｐ１，Ｐ２ＰＭ，ｄｐ＿ＣＶＰＭは第２の分析関数ｙ２＝ｆ１１（Ｐ１，Ｐ２ＰＭ，ｄｐ＿ＣＶＰＭ）に代入され、計算により第２の分析値ｙ２が取得される。取得された測定値Ｐ１，Ｐ２ＰＧ，ｄｐ＿ＣＶＰＧは第３の分析関数ｙ３＝ｆ１２（Ｐ１，Ｐ２ＰＧ，ｄｐ＿ＣＶＰＧ）に代入され、計算により第３の分析値ｙ３が取得される。取得された測定値Ｐ２ＰＭ，Ｐ３ＰＭ，Ｔ３ＰＭ，ｍ＿ＰＭは第４の分析関数ｙ４＝ｆ１３（Ｐ２ＰＭ，Ｐ３ＰＭ，Ｔ３ＰＭ，ｍ＿ＰＭ）に代入され、計算により第４の分析値ｙ４が取得される。取得された測定値Ｐ２ＰＧ，Ｐ３ＰＧ，Ｔ３ＰＧ，ｍ＿ＰＧは第５の分析関数ｙ５＝ｆ１４（Ｐ２ＰＧ，Ｐ３ＰＧ，Ｔ３ＰＧ，ｍ＿ＰＧ）に代入され、計算により第５の分析値ｙ５が取得される。取得された測定値Ｔ０，Ｔ３ＰＭ，Ｔ３ＰＧは第６の分析関数ｙ６＝ｆ１５（Ｔ０，Ｔ３ＰＭ，Ｔ３ＰＧ）に代入され、計算により第６の分析値ｙ６が取得される。取得された測定値Ｐ１，ｄｐ＿ＣＶＰＭ，Ｔ０，ｍ＿ＰＭは第７の分析関数ｙ７＝ｆ１６（Ｐ１，ｄｐ＿ＣＶＰＭ，Ｔ０，ｍ＿ＰＭ）に代入され、計算により第７の分析値ｙ７が取得される。取得された測定値Ｐ１，ｄｐ＿ＣＶＰＧ，Ｔ０，ｍ＿ＰＧは第８の分析関数ｙ８＝ｆ１７（Ｐ１，ｄｐ＿ＣＶＰＧ，Ｔ０，ｍ＿ＰＧ）に代入され、計算により第８の分析値ｙ８が取得される。取得された測定値Ｐ０，Ｐ１，Ｔ０，ｍ＿ＰＧ，ｍ＿ＰＭは第９の分析関数ｙ９＝ｆ１８（Ｐ０，Ｐ１，Ｔ０，ｍ＿ＰＧ，ｍ＿ＰＭ）に代入され、計算により第９の分析値ｙ９が取得される。プロセスはステップＳ３０へ移行する。

ステップＳ３０では、第１の分析関数の関数値ｙ１が、第１の分析範囲Ｙ１と比較される。第１の分析関数の関数値ｙ１が第１の分析範囲Ｙ１に入らない場合、標識Ｆ１が添付される。そうでなく、第１の分析関数の関数値ｙ１が第１の分析範囲Ｙ１に入る場合には、標識Ｔ１が添付される。

第２の分析関数の関数値ｙ２は、第２の分析範囲Ｙ２と比較される。第２の分析関数の関数値ｙ２が第２の分析範囲Ｙ２に入らない場合、標識Ｆ２が添付される。そうでなく、第２の分析関数の関数値ｙ２が第２の分析範囲Ｙ２に入る場合には、標識Ｔ２が添付される。

第３の分析関数の関数値ｙ３は、第３の分析範囲Ｙ３と比較される。第３の分析関数の関数値ｙ３が第３の分析範囲Ｙ３に入らない場合、標識Ｆ３が添付される。そうでなく、第３の分析関数の関数値ｙ３が第３の分析範囲Ｙ３に入る場合には、標識Ｔ３が添付される。

第４の分析関数の関数値ｙ４は、第４の分析範囲Ｙ４と比較される。第４の分析関数の関数値ｙ４が第４の分析範囲Ｙ４に入らない場合、標識Ｆ４が添付される。そうでなく、第４の分析関数の関数値ｙ４が第４の分析範囲Ｙ４に入る場合には、標識Ｔ４が添付される。

第５の分析関数の関数値ｙ５は、第５の分析範囲Ｙ５と比較される。第５の分析関数の関数値ｙ５が第５の分析範囲Ｙ５に入らない場合、標識Ｆ５が添付される。そうでなく、第５の分析関数の関数値ｙ５が第５の分析範囲Ｙ５に入る場合には、標識Ｔ５が添付される。

第６の分析関数の関数値ｙ６は、第６の分析範囲Ｙ６と比較される。第６の分析関数の関数値ｙ６が第６の分析範囲Ｙ６に入らない場合、標識Ｆ６が添付される。そうでなく、第６の分析関数の関数値ｙ６が第６の分析範囲Ｙ６に入る場合には、標識Ｔ６が添付される。

第７の分析関数の関数値ｙ７は、第７の分析範囲Ｙ７と比較される。第７の分析関数の関数値ｙ７が第７の分析範囲Ｙ７に入らない場合、標識Ｆ７が添付される。そうでなく、第７の分析関数の関数値ｙ７が第７の分析範囲Ｙ７に入る場合には、標識Ｔ７が添付される。

第８の分析関数の関数値ｙ８は、第８の分析範囲Ｙ８と比較される。第８の分析関数の関数値ｙ８が第８の分析範囲Ｙ８に入らない場合、標識Ｆ８が添付される。そうでなく、第８の分析関数の関数値ｙ８が第８の分析範囲Ｙ８に入る場合には、標識Ｔ８が添付される。

第９の分析関数の関数値ｙ９は、第９の分析範囲Ｙ９と比較される。第９の分析関数の関数値ｙ９が第９の分析範囲Ｙ９に入らない場合、標識Ｆ９が添付される。そうでなく、第９の分析関数の関数値ｙ９が第９の分析範囲Ｙ９に入る場合には、標識Ｔ９が添付される。プロセスは、それぞれ並列的にステップＳ４０，Ｓ４１，Ｓ４２，Ｓ４３，Ｓ４４，Ｓ４５，Ｓ４６，Ｓ４７，Ｓ４８，Ｓ４９，Ｓ５０，Ｓ５１，Ｓ５２に移行する。

ステップＳ４０では、ステップＳ３０での比較結果がＦ１，Ｔ２，Ｔ３，Ｔ４，Ｔ５，Ｔ６，Ｔ７，Ｔ８，Ｆ９であるか否かが判別される。判別結果がイエスである場合、第１の分析関数及び第９の分析関数の双方に含まれる測定値Ｐ０が測定誤差の最尤度を有することが検出され、分析プロセスは終了する。判別結果がノーである場合も分析プロセスは終了する。

ステップＳ４１では、ステップＳ３０での比較結果がＦ１，Ｆ２，Ｆ３，Ｔ４，Ｔ５，Ｔ６，Ｆ７，Ｆ８，Ｆ９であるか否かが判別される。判別結果がイエスである場合、第１の分析関数及び第２の分析関数及び第３の分析関数及び第７の分析関数及び第８の分析関数及び第９の分析関数の全てに共通に含まれる測定値Ｐ１が測定誤差の最尤度を有することが検出され、分析プロセスは終了する。判別結果がノーである場合も分析プロセスは終了する。

ステップＳ４２では、ステップＳ３０での比較結果がＴ１，Ｆ２，Ｔ３，Ｆ４，Ｔ５，Ｔ６，Ｔ７，Ｔ８，Ｔ９であるか否かが判別される。判別結果がイエスである場合、第２の分析関数及び第４の分析関数の双方に含まれる測定値Ｐ２ＰＭが測定誤差の最尤度を有することが検出され、分析プロセスは終了する。判別結果がノーである場合も分析プロセスは終了する。

ステップＳ４３では、ステップＳ３０での比較結果がＴ１，Ｔ２，Ｆ３，Ｔ４，Ｆ５，Ｔ６，Ｔ７，Ｔ８，Ｔ９であるか否かが判別される。判別結果がイエスである場合、第３の分析関数及び第５の分析関数の双方に含まれる測定値Ｐ２ＰＧが測定誤差の最尤度を有することが検出され、分析プロセスは終了する。判別結果がノーである場合も分析プロセスは終了する。

ステップＳ４４では、ステップＳ３０での比較結果がＴ１，Ｔ２，Ｔ３，Ｔ４，Ｆ５，Ｔ６，Ｔ７，Ｔ８，Ｔ９であるか否かが判別される。判別結果がイエスである場合、測定値Ｐ３ＰＧが測定誤差の最尤度を有することが検出され、分析プロセスは終了する。判別結果がノーである場合も分析プロセスは終了する。

ステップＳ４５では、ステップＳ３０での比較結果がＴ１，Ｔ２，Ｔ３，Ｆ４，Ｔ５，Ｔ６，Ｔ７，Ｔ８，Ｔ９であるか否かが判別される。判別結果がイエスである場合、測定値Ｐ３ＰＭが測定誤差の最尤度を有することが検出され、分析プロセスは終了する。判別結果がノーである場合も分析プロセスは終了する。

ステップＳ４６では、ステップＳ３０での比較結果がＴ１，Ｆ２，Ｔ３，Ｔ４，Ｔ５，Ｔ６，Ｆ７，Ｔ８，Ｔ９であるか否かが判別される。判別結果がイエスである場合、第２の分析関数及び第７の分析関数の双方に含まれる測定値ｄｐ＿ＣＶＰＭが測定誤差の最尤度を有することが検出され、分析プロセスは終了する。判別結果がノーである場合も分析プロセスは終了する。

ステップＳ４７では、ステップＳ３０での比較結果がＴ１，Ｔ２，Ｆ３，Ｔ４，Ｔ５，Ｔ６，Ｔ７，Ｆ８，Ｔ９であるか否かが判別される。判別結果がイエスである場合、第３の分析関数及び第８の分析関数の双方に含まれる測定値ｄｐ＿ＣＶＰＧが測定誤差の最尤度を有することが検出され、分析プロセスは終了する。判別結果がノーである場合も分析プロセスは終了する。

ステップＳ４８では、ステップＳ３０での比較結果がＦ１，Ｔ２，Ｔ３，Ｔ４，Ｔ５，Ｆ６，Ｆ７，Ｆ８，Ｆ９であるか否かが判別される。判別結果がイエスである場合、第１の分析関数及び第６の分析関数及び第７の分析関数及び第８の分析関数及び第９の分析関数の全てに共通に含まれる測定値Ｔ０が測定誤差の最尤度を有することが検出され、分析プロセスは終了する。判別結果がノーである場合も分析プロセスは終了する。

ステップＳ４９では、ステップＳ３０での比較結果がＴ１，Ｔ２，Ｔ３，Ｆ４，Ｔ５，Ｆ６，Ｔ７，Ｔ８，Ｔ９であるか否かが判別される。判別結果がイエスである場合、第４の分析関数及び第６の分析関数の双方に含まれる測定値Ｔ３ＰＭが測定誤差の最尤度を有することが検出され、分析プロセスは終了する。判別結果がノーである場合も分析プロセスは終了する。

ステップＳ５０では、ステップＳ３０での比較結果がＴ１，Ｔ２，Ｔ３，Ｔ４，Ｆ５，Ｆ６，Ｔ７，Ｔ８，Ｔ９であるか否かが判別される。判別結果がイエスである場合、第５の分析関数及び第６の分析関数の双方に含まれる測定値Ｔ３ＰＧが測定誤差の最尤度を有することが検出され、分析プロセスは終了する。判別結果がノーである場合も分析プロセスは終了する。

ステップＳ５１では、ステップＳ３０での比較結果がＴ１，Ｔ２，Ｔ３，Ｆ４，Ｔ５，Ｔ６，Ｆ７，Ｔ８，Ｆ９であるか否かが判別される。判別結果がイエスである場合、第４の分析関数及び第７の分析関数及び第９の分析関数の全てに共通に含まれる測定値ｍ＿ＰＭが測定誤差の最尤度を有することが検出され、分析プロセスは終了する。判別結果がノーである場合も分析プロセスは終了する。

ステップＳ５２では、ステップＳ３０での比較結果がＴ１，Ｔ２，Ｔ３，Ｔ４，Ｆ５，Ｔ６，Ｔ７，Ｆ８，Ｆ９であるか否かが判別される。判別結果がイエスである場合、第５の分析関数及び第８の分析関数及び第９の分析関数の全てに共通に含まれる測定値ｍ＿ＰＧが測定誤差の最尤度を有することが検出され、分析プロセスは終了する。判別結果がノーである場合も分析プロセスは終了する。

測定誤差が測定値Ｐ０に発生したことが検出された場合、第１の分析範囲Ｙ１の中間値がｙ１＝ｆ１０（Ｐ０，Ｐ１，Ｔ０）へ代入される。ここで、式中の測定値Ｐ１，Ｔ０は既知であるので、Ｐ０に対応する１つの計算値が計算により取得される。当該計算値は、圧力センサ１４の組み込み位置での実際のガス圧を反映可能である。

測定誤差が測定値Ｐ１に発生したことが検出された場合、第２の分析範囲Ｙ２の中間値がｙ２＝ｆ１１（Ｐ１，Ｐ２ＰＭ，ｄｐ＿ＣＶＰＭ）へ代入される。ここで、式中の測定値Ｐ２ＰＭ，ｄｐ＿ＣＶＰＭは既知であるので、Ｐ１に対応する１つの計算値が計算により取得される。当該計算値は、圧力センサ１８の組み込み位置での実際のガス圧を反映可能である。

測定誤差が測定値Ｐ２ＰＭに発生したことが検出された場合、第２の分析範囲Ｙ２の中間値がｙ２＝ｆ１１（Ｐ１，Ｐ２ＰＭ，ｄｐ＿ＣＶＰＭ）へ代入される。ここで、式中の測定値Ｐ１，ｄｐ＿ＣＶＰＭは既知であるので、Ｐ２ＰＭに対応する１つの計算値が計算により取得される。当該計算値は、圧力センサ２１の組み込み位置での実際のガス圧を反映可能である。

測定誤差が測定値Ｐ２ＰＧに発生したことが検出された場合、第３の分析範囲Ｙ３の中間値がｙ３＝ｆ１２（Ｐ１，Ｐ２ＰＧ，ｄｐ＿ＣＶＰＧ）へ代入される。ここで、式中の測定値Ｐ１，ｄｐ＿ＣＶＰＧは既知であるので、Ｐ２ＰＧに対応する１つの計算値が計算により取得される。当該計算値は、圧力センサ３１の組み込み位置での実際のガス圧を反映可能である。

測定誤差が測定値Ｐ３ＰＭに発生したことが検出された場合、第４の分析範囲Ｙ４の中間値がｙ４＝ｆ１３（Ｐ２ＰＭ，Ｐ３ＰＭ，Ｔ３ＰＭ，ｍ＿ＰＭ）へ代入される。ここで、式中の測定値Ｐ２ＰＭ，Ｔ３ＰＭ，ｍ＿ＰＭは既知であるので、Ｐ３ＰＭに対応する１つの計算値が計算により取得される。当該計算値は、圧力センサ２６の組み込み位置での実際のガス圧を反映可能である。

測定誤差が測定値Ｐ３ＰＧに発生したことが検出された場合、第５の分析範囲Ｙ５の中間値がｙ５＝ｆ１４（Ｐ２ＰＧ，Ｐ３ＰＧ，Ｔ３ＰＧ，ｍ＿ＰＧ）へ代入される。ここで、式中の測定値Ｐ２ＰＧ，Ｔ３ＰＧ，ｍ＿ＰＧは既知であるので、Ｐ３ＰＧに対応する１つの計算値が計算により取得される。当該計算値は、圧力センサ３６の組み込み位置での実際のガス圧を反映可能である。

測定誤差が測定値ｄｐ＿ＣＶＰＭに発生したことが検出された場合、第２の分析範囲Ｙ２の中間値がｙ２＝ｆ１１（Ｐ１，Ｐ２ＰＭ，ｄｐ＿ＣＶＰＭ）へ代入される。ここで、式中の測定値Ｐ１，Ｐ２ＰＭは既知であるので、ｄｐ＿ＣＶＰＭに対応する１つの計算値が計算により取得される。当該計算値は、圧力差センサ２５の組み込み位置での実際のガス圧力差を反映可能である。

測定誤差が測定値ｄｐ＿ＣＶＰＧに発生したことが検出された場合、第３の分析範囲Ｙ３の中間値がｙ３＝ｆ１２（Ｐ１，Ｐ２ＰＧ，ｄｐ＿ＣＶＰＧ）へ代入される。ここで、式中の測定値Ｐ１，Ｐ２ＰＧは既知であるので、ｄｐ＿ＣＶＰＧに対応する１つの計算値が計算により取得される。当該計算値は、圧力差センサ３５の組み込み位置での実際のガス圧力差を反映可能である。

測定誤差が測定値Ｔ０に発生したことが検出された場合、第６の分析範囲Ｙ６の中間値がｙ６＝ｆ１５（Ｔ０，Ｔ３ＰＭ，Ｔ３ＰＧ）へ代入される。ここで、式中のＴ３ＰＭ，Ｔ３ＰＧは既知であるので、Ｔ０に対応する１つの計算値が計算により取得される。当該計算値は、温度センサ１６の組み込み位置での実際のガス温度を反映可能である。

測定誤差が測定値Ｔ３ＰＭに発生したことが検出された場合、第６の分析範囲Ｙ６の中間値がｙ６＝ｆ１５（Ｔ０，Ｔ３ＰＭ，Ｔ３ＰＧ）へ代入される。ここで、式中のＴ０，Ｔ３ＰＧは既知であるので、Ｔ３ＰＭに対応する１つの計算値が計算により取得される。当該計算値は、温度センサ２８の組み込み位置での実際のガス温度を反映可能である。

測定誤差が測定値Ｔ３ＰＧに発生したことが検出された場合、第６の分析範囲Ｙ６の中間値がｙ６＝ｆ１５（Ｔ０，Ｔ３ＰＭ，Ｔ３ＰＧ）へ代入される。ここで、式中のＴ０，Ｔ３ＰＭは既知であるので、Ｔ３ＰＧに対応する１つの計算値が計算により取得される。当該計算値は、温度センサ３８の組み込み位置での実際のガス温度を反映可能である。

測定誤差が測定値ｍ＿ＰＭに発生したことが検出された場合、第７の分析範囲Ｙ７の中間値がｙ７＝ｆ１６（Ｐ１，ｄｐ＿ＣＶＰＭ，Ｔ０，ｍ＿ＰＭ）へ代入される。ここで、式中のＰ１，ｄｐ＿ＣＶＰＭ，Ｔ０は既知であるので、ｍ＿ＰＭに対応する１つの計算値が計算により取得される。当該計算値は、流量センサ２４の組み込み位置での実際のガス流量を反映可能である。

測定誤差が測定値ｍ＿ＰＧに発生したことが検出された場合、第８の分析範囲Ｙ８の中間値がｙ８＝ｆ１７（Ｐ１，ｄｐ＿ＣＶＰＧ，Ｔ０，ｍ＿ＰＧ）へ代入される。ここで、式中のＰ１，ｄｐ＿ＣＶＰＧ，Ｔ０は既知であるので、ｍ＿ＰＧに対応する１つの計算値が計算により取得される。当該計算値は、流量センサ３４の組み込み位置での実際のガス流量を反映可能である。

分析方法を実行するためのプログラムはソフトウェアによって実現可能であり、ソフトウェアはプログラマブルコントローラ３３に記憶可能である。コントローラ３３は、圧力センサ１４，１８，２１及び温度センサ１６及び流量センサ２４及び圧力センサ２６などで検出された、ガスタービンの動作パラメータの測定値を受信する。次に、測定誤差の発生している測定値を検出するために、上述した方法の各ステップが実行される。また、コントローラ３３は、どのセンサで測定誤差が発生したかを表示する信号を送信することもできる。

本発明の別の特徴として、ガスタービンの燃料ガス装置が開示される。燃料ガス装置は、燃料ガスをガスタービンへ輸送する燃料管１０，２０，３０と、燃料管１０，２０，３０内の燃料の物理的な動作パラメータの複数の値を測定するセンサ１４，１６，１８，２１，２４，２５，２６，２８，３１，３４，３５，３６，３８と、プログラマブルコントローラ３３とを含む。コントローラ３３は、どのセンサで測定誤差が発生したかを報知することができる。

各実施形態を具体例に則して説明したが、各実施形態が唯一の独立した技術手段しか含まないわけではないことを理解されたい。こうした具体例に即した説明の方式は単に明瞭性のためのものであって、当業者はこうした具体例を全体として捉えるべきであり、各実施形態の技術手段を相互に適切に組み合わせて当業者に理解可能な他の実現形態を形成することもできる。

上述した一連の説明は本発明の可能な実施形態の単なる詳細な説明であって、本発明の権利範囲を限定する意図のものでない。各特徴と等価のあらゆる実現手段、又は、本発明の技術及び思想を逸脱することなく行われた各特徴の結合もしくは分割もしくは反復などの変更も、本発明の権利範囲に該当する。

Claims (8)

1. ガスタービンの動作パラメータ（Ｐ０，Ｐ１，Ｐ２ＰＭ，Ｐ２ＰＧ，Ｐ３ＰＭ，Ｐ３ＰＧ，ｄｐ＿ＣＶＰＧ，ｄｐ＿ＣＶＰＭ，ｍ＿ＰＧ，ｍ＿ＰＭ，Ｔ０，Ｔ３ＰＭ，Ｔ３ＰＧ）の測定誤差を検出するための分析方法であって、
   動作パラメータを、それぞれ予め定められた分析範囲（Ｙ１，Ｙ２…Ｙｍ）を有する複数の関数（ｙ１，ｙ２…ｙｍ、ここでｍ≧３）の変数の少なくとも一部としたとき、当該動作パラメータ（Ｓ１，Ｓ２…Ｓｎ，ここでｎ≧３）の複数の測定値を取得するステップと、
   前記動作パラメータの前記測定値を、前記動作パラメータを含む変数を有する関数へ投入して、関数の結果を計算するステップと、
   関数の計算結果を、対応する関数の分析範囲と比較するステップと、
   前記計算結果が対応する関数の予め定められた分析範囲に入らない関数の変数の少なくとも一部である動作パラメータの１つもしくは複数の測定値を求めるステップと、
   求められた前記１つもしくは複数の測定値のなかで、前記計算結果が対応する関数の予め定められた分析範囲に入る関数の変数ではない動作パラメータの測定値が測定誤差の最尤度を有することを検出するステップと
   を含む、分析方法。
2. 前記分析方法はさらに、
   測定誤差の最尤度を有すると検出された測定値の計算値を、当該測定値を含むいずれかの関数の分析範囲の中間値を用いて取得するステップと、
   当該測定値に代えて、補正値としての前記計算値を用いるステップと
   を含む、
   請求項１記載の分析方法。
3. 前記分析方法はさらに、
   測定誤差の最尤度を有すると検出された測定値を含む関数の計算値を、前記測定値を含む関数の分析範囲の中間値を用いて取得し、前記計算値の平均を取得するステップと、
   当該測定値に代えて、補正値としての前記計算値を用いるステップと
   を含む、
   請求項１記載の分析方法。
4. ガスタービンの動作パラメータは、ガスタービン内の流体の温度、圧力、流量及び圧力差を含む、
   請求項１記載の分析方法。
5. 前記動作パラメータ（Ｓ１，Ｓ２…Ｓｎ）の複数の測定値は、ガスタービンのオイル装置、空気装置及びガス装置の動作パラメータの測定値を含む、
   請求項４記載の方法。
6. ガスタービンの動作パラメータの測定誤差を検出するための分析方法を実行するようにプログラミングされた、ガスタービン用コントローラであって、特に請求項１から５までのいずれか１項記載の測定誤差を検出するための分析方法、すなわち、
   前記動作パラメータの前記測定値を、前記動作パラメータを含む変数を有する関数へ投入して、関数の結果を計算するステップと、
   前記関数の計算結果を、対応する関数の分析範囲と比較するステップと、
   前記計算結果が対応する関数の予め定められた分析範囲に入らない関数の変数の少なくとも一部である動作パラメータの１つもしくは複数の測定値を求めるステップと、
   求められた前記１つもしくは複数の測定値のなかで、前記計算結果が対応する関数の予め定められた分析範囲に入る関数の変数ではない動作パラメータの測定値が測定誤差の最尤度を有することを検出するステップと
   を含む分析方法
   を実行するようにプログラミングされた、ガスタービン用コントローラ。
7. ガスタービンの燃料ガス装置であって、
   燃料ガスをガスタービンへ輸送する燃料管（１０，２０，３０）と、
   前記燃料管（１０，２０，３０）内の燃料の物理的な動作パラメータの複数の値を測定する複数のセンサ（１４，１６，１８，２１，２４，２５，２６，２８，３１，３４，３５，３６，３８）と、
   請求項６に記載のプログラマブルコントローラと
   を含む、燃料ガス装置。
8. 前記燃料ガス装置はさらに、どのセンサで測定誤差が発生したかを報知するコントローラを含む、
   請求項７記載の燃料ガス装置。

Images