JP4361582B2 - ガスタービンの性能診断方法及び性能診断システム - Google Patents

Description

本発明はガスタービンの性能の診断方法及び診断システムに関する。

ガスタービンの性能は、燃料流量などの操作条件だけでなく、吸気温度などの大気条件によっても大きく変化するため、単純にそのまま絶対値では評価できない。そこで対応策として、実機の運転条件を性能計算プログラムに入力し、計算された性能の計算値すなわち、あるべき予測値に対して、実測値を比較して評価するのが一般的である。

実際には、ガスタービンの性能計算プログラムの計算値が実測値と完全に一致することはないため、このずれをなくすように、計算の仕方に何らかの補正が加えられる。これは一般にチューニング、モデル調整などと呼ばれる。その方法として、例えば特許文献１では、性能モデルの計算値が実測値と一致するようにモデルの公称値を調整し、この調整結果に基づいてモデル計算値と計測値を比較し、欠陥の有無を診断する方法を提唱している。

特開２００１−１７４３６６号公報

しかし、特許文献１をはじめとする従来の調整方法は、以下に述べるように複雑になりがちで、多量の分析が必要なため、多大な時間と労力を要するという課題があった。例えば、或る計測項目について、実測値と計算値が一致するように、設備の或る項目の公称値を調整するとする。この場合、目的の計測項目は計算値と実測値が一致するようになるが、ガスタービン内部の現象の因果関係によって、他の計測項目の計算値も変化し、実測値と一致しなくなることがある。このような副作用を避けるためには、さらに別の項目の公称値を調整する必要がある。しかし、これによってまた別の計測項目の計算値が変化し、実測値と一致しなくなるおそれがある。このような公称値の調整の連鎖をどのように扱うか、影響の低い項目を切り捨てて収束させるにはどのような基準で扱えばよいかは明確になっていない。

従って、これに対応するには、実際の現象の複雑な因果関係とその影響範囲を精査し、どの計測値のずれに対して公称値のどの項目を調整するか、そのときの他の計測値への影響はどの程度か、これら全体の複雑な相互関係を適切に把握し、適切な調整手順を定めることが必要になる。これには、複雑で多量の事前分析と、時間と労力が必要になる。

本発明の目的は、複雑な手順によらず、数少ない因子の調整によって、ガスタービンの性能診断を高精度化できるようにしたガスタービンの性能診断方法と性能診断システムを提供することにある。

本発明は、実機で計測された吸気温度と圧縮機圧力比と燃料流量の入力を受けて、ガスタービンの発電出力と排気温度を計算して出力する性能計算回路と、該性能計算回路から出力された発電出力と排気温度の実機計測値からの偏差に基づいて性能を評価する評価回路を備えたガスタービン性能診断システムであって、吸気温度と圧縮機圧力比と燃料流量の実機計測値に基づいて、燃料流量の値の補正値を計算し、かかる計算で得られた補正燃料流量を前記性能計算回路への入力として設定する補正回路を備えたことを特徴とする。

また、本発明は、ガスタービンの性能計算回路に実機で計測された吸気温度と圧縮機圧力比と燃料流量を入力して発電出力と排気温度を計算し、該計算で得られた発電出力と排気温度が実機での計測値に対して有する偏差に基づいて性能を評価するガスタービンの性能診断方法であって、吸気温度と圧縮機圧力比と燃料流量の実機計測値に基づいて、予め定められた関数関係に従って、燃料流量の値を補正し、かかる計算で得られた補正燃料流量を前記性能計算回路に入力することを特徴とする。

本発明の診断システムには、吸気温度と圧縮機圧力比と燃料流量の実機計測値についての一定期間の時系列データまたは運転条件を複数に変えたデータ（これらを以下、複数データセットと呼ぶ）を入力して計算出力された発電出力と排気温度が実機計測値に対して有する偏差が、予め定められた許容基準内になるように、前記補正回路の関数関係を設定する演算回路（以下、補正関数設定回路と呼ぶ）を設けることができる。

また補正回路では、吸気温度と圧縮機圧力比と燃料流量の実機計測値に基づいて、予め定められた関数関係に従って、燃料流量の値を補正する計算を実行し、実機で計測された燃料流量の値を、補正された燃料流量の値で差替えて性能計算回路へ入力する。

本発明により、複雑な手順によらず、数少ない因子の調整によって、ガスタービンの性能診断を高精度化することができた。本発明では、実機の計測値と性能計算値とのずれに関わる諸要因に対して、数多くの情報を使った複雑な手順によらず、３項目の入力情報を用いて１つの入力情報を補正するだけで、誤差に関係する諸要因を包括的に網羅して効率的にシンプルに調整できる。また、この補正のために必要な情報も５項目の計測値情報で済む。

我々は、ガスタービンの性能計算において、計算値と実測値のずれを引き起こす要因を分析し、有力な要因の１つとして実測値の熱収支の不均衡に着目した。そして、この不均衡に影響する因子を分析し、影響が大きい重要因子を摘出した。これに基づき、この熱収支の不均衡によって生じる誤差を解消するモデル調整手段を構築し、本発明を見出した。以下、本発明の内容を詳細に説明する。
（従来方式の分析とモデル立案）
本発明と対比するため、はじめに、従来行ってきたガスタービン性能診断システムについて説明する。

図６は、従前の典型的なガスタービン診断システムの構成図である。以下、これを従来方式と呼ぶ。ここでは後述のように課題を分析するため、モデル調整しない場合を示している。この診断システム６０は、実機計測値情報１（以下、実測値と略するときがある）のうち性能計算入力用の情報２を性能計算回路７に入力する入力手段５と、入力された前記情報２に基づいてガスタービンの性能の予測値を計算する性能計算回路７と、該計算の結果出力された計算結果情報８と前記実機計測値情報１のうち評価用の情報３との偏差に基づいて、実測された性能の状態を評価する評価回路９と、該評価回路９から出力された評価結果情報を格納する格納手段１０と、該格納手段１０に記憶された評価結果情報または前記評価回路９から出力された前記評価結果情報を診断情報１３として外部出力する出力手段１１からなる。

この従来方式の診断システム６０の動作原理は以下のようなものである。

入力となる実機計測値情報１は、ガスタービンの実機の運転条件や運転状態の計測値であり、運転条件には、ガスタービン本体の燃料流量などの操作条件の他に、ガスタービン性能に影響する大気温度・圧力・湿度などの環境条件が含まれる。これら環境条件の計測値は実機が設置された現地の計測値でなく、最寄の気象台のデータなど、実機データに代用できる別のデータのこともある。

実機計測値情報１は入力用の情報２と評価用の情報３の２通りの用いられ方がある。入力用の情報２の代表的なものとしては、吸気温度と圧縮機圧力比と燃料流量があり、これらは性能計算回路７に入力する運転条件の情報として用いられる。評価用の情報３の代表的なものとしては、発電出力と排気温度があり、これらは性能計算回路７で計算された発電出力と排気温度と比較して実機の状態を評価するのに用いられる。

入力手段５は、このような実機計測値情報１のうち入力用の情報２を、実機から、あるいは実機に接続したデータサーバから、有線または無線の通信回線を通じて電子データとして受け取り、性能計算回路７に入力する。具体的には、入力用の情報２を通信回路から情報受信装置によって受信し、受信した情報を情報変換手段によって性能計算回路７への入力書式に変換して入力する。実機計測値情報１からの入力用の情報２の分離は、入力手段５で受信される前を含めて、性能計算回路７に入力されるまで任意の段階で実行される。

また、入力手段５は、同じようにして、実機計測値情報１のうち評価用の情報３を、実機から、あるいは実機に接続したデータサーバから、有線または無線の通信回線を通じて電子データとして受け取り、評価回路９に入力する。具体的には、評価用の情報３を通信回路から情報受信装置によって受信し、受信した情報を情報変換手段によって性能計算回路９への入力書式に変換して入力する。実機計測値情報１からの評価用の情報３の分離は、入力手段５で受信される前を含めて、評価回路９に入力されるまでの任意の段階で実行される。

性能計算回路７は、前記入力手段５から性能計算入力用の実機計測値の情報２の入力を受け、ガスタービンの予想される性能を計算し、計算結果情報８を出力する。計算結果情報８には典型的には発電出力と排気温度が含まれる。この性能計算回路７は、熱力学的・空力的なサイクル計算、あるいはニューラルネットのような人工知能的なアルゴリズムなど、公知の計算方法をコンピュータプログラムまたは専用の計算回路として実装したものである。

なお、本発明において回路と記すものは、実装手段としては、コンピュータプログラムまたは演算チップ、専用回路など、電子的手段によって実現された公知の計算手段を指すものとする。また、これらは１つの回路単体で動作する場合だけでなく、複数の回路によって１つのプログラムや回路を構成する部品であってもよい。例えば、前記性能計算回路７と後に述べる評価回路９は、機能構成を明確に説明するためにそれぞれ別の回路名で区別しているが、別々の異なるハードウェアである必要はなく、実際には１つの診断プログラムを構成するサブモジュールとして実装されるのが一般的である。

前記性能計算回路７によって計算された計算結果情報８は、電子情報として出力され、評価回路９に入力される。

評価回路９では、計算結果情報８の発電出力と排気温度の値が、前記評価用の実機計測値の情報３との偏差に基づいて評価される。例えば、発電出力または排気温度についての、計算結果情報８と評価用実測値の情報３の値の偏差が、予め設定したしきい値を超過すれば、警告を発報するための情報が生成される。あるいは、次に述べる格納手段１０に格納された過去の蓄積結果を読み込んで、その時系列変化から異常の有無の判定処理や、異常の程度の数量化処理、性能の変化量の計算などを実行することもある。本発明では、これら評価回路９の諸処理を総称して比較評価、またその結果情報を総称して比較結果情報と呼ぶことがある。

前記評価回路９での比較結果情報は格納手段１０に格納される。格納される情報は典型的には実機計測値情報１と計算結果情報８である。格納手段１０はハードディスクや光ディスク、メモリなど公知の情報格納手段である。一般には評価回路９が実行されるごとに、これらの比較結果情報が順次、格納手段１０に蓄積される。

前記評価回路９での比較結果情報、あるいは前記格納手段１０に格納された情報は、診断情報１３として出力手段１１によって外部出力される。該出力手段１１は診断システム６０から外部への情報出力手段のことであり、具体的には、モニタなどの表示画面、あるいはファイル転送プロトコルのようなファイル転送手段、データファイルを外部出力するファイル出力手段、ハードディスクなどの情報記録手段、プリンタのような電子媒体から別媒体への出力手段、などである。

このような従来方式の有効性と精度、課題を分析するために、実際にガスタービンの性能を評価した事例を図７に示す。ここで誤差率は、正の場合は実測値よりも計算値が大きいことを、負の場合は実測値よりも計算値が小さいことを示す。図７の従来事例では、発電出力と排気温度について、実測値と計算値で概ね最大２％の誤差がある。そして、図７の領域Ａでは、発電出力と排気温度の両方について実測値よりも計算値の方が大きく、逆に領域Ｂでは、発電出力と排気温度の両方について実測値よりも計算値の方が小さい。

これはガスタービンの入出熱収支、すなわち、ガスタービン設備への入力側にあたる吸込空気と燃料の熱量和（入熱）と、出力側にあたる発電出力と排ガスと損失の熱量和（出熱）が合っていないことを示す。すなわち、実測値を基準にするならば、計算結果は領域Ａでは入熱に対して出熱が過剰であり、領域Ｂでは出熱が不足している。しかし、本ケースの性能計算回路７の計算モデルは入出熱収支が均衡するように計算が実行されている。このことは、入出力の熱収支が均衡していないのは、性能計算回路７の計算モデルではなく、実測値の側であることを示している。これが我々の着目点である。

このような実測値の入出熱収支の偏差の原因としては種々考えられるが、その１つに排気温度がある。排気温度はガスタービンのタービン出口側で周方向に１０数点が計測され、一般にはその中央値または平均値が代表値として性能評価やプラントの運転制御に用いられる。しかし、実際には、タービン出口から排ガスダクトへの温度降下や、温度計測断面での周方向と直径方向の温度分布、流れの旋回など複雑な要因があるため、前記の代表値は入出熱収支上の排気温度の真値とは一致しない。

しかし、これらの複雑な諸要因を全て考慮して排気温度の真値を推定するのは極めて困難である。そこで、これらの熱収支と誤差の問題を解決し、高精度に性能を評価できる、簡便で効率的な方策が必要である。

上述の排気温度についての複雑な現象を考慮しながらも、実用に使える簡便な解決方法を見出すため、我々は、まず熱力学的考察により熱収支の偏差に支配的に影響する重要因子を摘出し、そして、これに基づくモデル調整方法を立案した。

はじめに、熱収支の偏差に影響する重要因子の摘出については、以下のように、ガスタービンの動力サイクルにそって関連する因子を整理し、さらに単純化することにより、重要因子を絞り込んだ。

すなわち、まず計測された入出熱収支の偏差ΔＨについて、ガスタービンの動作原理から、圧縮仕事に影響する圧縮機効率ηｃと、タービン仕事に影響するタービン効率ηｔと、入熱量の大半を占める投入燃料熱量ＱＦの３要素の関数として式（１）のようにモデル化した。ここでｆは関数を示す。以下においても同様である。

ΔＨ ＝ ｆ（ηｃ， ηｔ， ＱＦ） …式（１）
式（１）の右辺の圧縮機効率とタービン効率は、熱力学的には各出入口の温度と圧力から計算できるが、単純化して出口温度を省略し、またタービンの圧力比は圧縮機の圧力比で代用し、式（２）、式（３）の関係を仮定した。ここで、Ｔｃｉは吸気温度、Ｔｔｉはタービン入口温度、Φｃは圧縮機圧力比を示す。

ηｃ ＝ ｆ（Ｔｃｉ， Φｃ） …式（２）
ηｔ ＝ ｆ（Ｔｔｉ， Φｃ） …式（３）
式（３）の右辺のタービン入口温度Ｔｔｉは、熱収支から式（４）のように考えた。ここでＴｃｄは圧縮機出口温度を示す。

Ｔｔｉ ＝ ｆ（Ｔｃｄ， ＱＦ） …式（４）
式（４）の右辺の圧縮機出口温度Ｔｃｄは、熱力学的に式（５）のように示される。

Ｔｃｄ ＝ ｆ（Ｔｃｉ， ηｃ） …式（５）
以上の式（１）〜式（５）の関係を整理して、式（６）を得た。

ΔＨ ＝ ｆ（Ｔｃｉ， Φｃ， ＱＦ） …（６）
このようにして、入出熱収支の偏差に支配的に影響する重要因子として、大きくは吸気温度と、圧縮機圧力比、燃料熱量の３つを摘出した。

次に、こうにして摘出した重要因子である吸気温度と、圧縮機圧力比、燃料熱量を用いて熱収支偏差を補正する方法を検討した。偏差の不均衡を補正するために、具体的にどのような因子を補正するのが適切か検討し、燃料流量を補正することとした。

この根拠は以下の通りである。すなわち、熱収支の不均衡を補正するためには、モデルの計算部分の内部か、入出力のいずれかを補正する必要がある。このうち、計算部分の内部を調整・改造するのは、モデルは熱収支が取れていることが重要な前提であり、これを崩すことになるため不適切である。また、出力側にあたる発電出力と排気温度は、性能解析で性能の劣化や異常を診断する根拠として用いられるため、これらを補正するのも不適切である。従って、入力側にあたる燃料流量か、吸気流量を補正する必要がある。このうち熱量的には燃料流量が圧倒的に大きく影響も大きいことから、補正項目として燃料流量を選定した。

このように補正項目を燃料流量にしたことに基づき、熱収支偏差を補正するための補正式を次のようにした。式（８）では関数を単純化して線形仮定し、また式（７）、式（８）ともに、計算しないと得られない燃料熱量の代わりに、直接計測できる燃料流量で代用した。

Ｇｆ＿ｃｏｒ ＝ γ×Ｇｆ …式（７）
γ ＝ ｆ（Ｔｃｉ，Φｃ，Ｇｆ）
＝ （ａ１×Ｔｃｉ ＋ ａ２×Φｃ ＋ ａ３×Ｇｆ）＋ｂ …式（８）
ここで、Ｇｆは燃料流量、Ｇｆ＿ｃｏｒは補正された燃料流量、γは補正係数、ａｉ（ｉ＝１，２，３）は係数、ｂは定数項を示す。
（発明システムの構成）
図１に以上の分析に基づく本発明のガスタービン診断システム４の構成図を示す。本システムは、従来方式の図６の構成要素に加え、性能計算用に入力する実機計測値情報のうち燃料流量の値を前述の式（７）、式（８）の関係によって補正する補正回路６と、該補正回路６の内容を設定する補正関数設定回路１２を備えている。

補正回路６は、式（７）、式（８）の関数が実装された計算回路である。この回路では、これらの計算式に従って、実測値の燃料流量を吸気温度と圧縮機圧力比と燃料流量の関数によって補正する計算を実行し、性能計算入力用データのうち、もとの燃料流量の値を補正された燃料流量（補正燃料流量）の値で差替えて性能計算回路７に入力する。

ここで、式（７）、式（８）の具体的内容、すなわち、式（８）の関数ｆの内容、あるいは係数ａ１〜ａ３、定数項ｂの値は、補正関数設定回路１２によって予め設定されている。補正関数設定回路１２は、図２に示すように、ある手順で取得された複数のデータ（以下、複数データセットと呼ぶ）２０に基づいて、前記補正回路６内の補正関数（式（７）（８））の未定定数を具体的に決定する計算を実行する。

これに用いられる複数データセット２０とは、実機の吸気温度、圧縮機圧力比、燃料流量、発電出力、排気温度を含む計測値についての、特定期間の時系列のデータ、または運転条件を複数に変えたデータを指す。以下、特定の１条件または１時点での実機の吸気温度、圧縮機圧力比、燃料流量、発電出力、排気温度を含む計測値一式を１つのデータセット、複数条件または複数時点の場合をまとめて複数データセットと呼んで区別する。また、吸気温度と圧縮機圧力比と燃料流量を入力用の情報２、発電出力と排気温度を評価用の情報３と呼んで適宜区別する。

なお、前記式（８）などに示す関数ｆは、数式化されていなくとも、右辺のカッコ内の項目を入力に含み、左辺の項目を出力に含むアルゴリズムを含めて指すものである。したがって、前記補正関数設定回路１２における補正関数の設定とは、式（８）の関数のｆを表す数式と、これに含まれる係数などの定数の決定に限らず、広義には、明示的に陽に数式化されていないアルゴリズムであっても、その計算の実行に必要な設定値の決定を含むものである。すなわち、吸気温度と圧縮機圧力比と燃料流量を入力すると、補正された燃料流量が出力されるような計算が実行されるために必要な設定全般をさす。
（補正関数設定回路の処理フロー）
補正関数設定回路１２の実行内容を詳しく説明するために、図３にその処理フローを示す。処理フローは、大きく分けて、前記複数データセット２０の各々のデータセットに対して、性能計算による発電出力と排気温度の計算値が実測値に対して有する偏差が、予め定めた許容基準内になるような燃料流量の修正値（修正燃料流量）と修正係数（後述）を収束計算により求める前半部分（工程２２〜２８）と、こうして求められた全データセット（複数データセット２０の全て）の修正係数とこれに対応する吸気温度、圧縮機圧力比、燃料流量のデータから、式（８）の係数と定数項をフィッティングして決定する後半部分（工程２９）からなる。以下、各工程を順に説明する。

性能計算工程２２では、後述の工程２８で次のデータセットに進む繰返しループにおける、その時点の（複数データセット２０のうちの）単一データセットを入力として、ガスタービンの性能計算が実施される。この性能計算の内容は既述した性能計算回路７（図６）で実行される計算のことである。

前記性能計算工程２２で計算された性能計算結果、すなわち発電出力と排気温度の計算値は、当該データセットに含まれる実測値すなわち評価用の情報２と比較され、計算値と実測値の偏差が予め定められた許容基準内にあるかどうかが、判定工程２３によって判定される。

この判定結果が偽であれば、続く修正工程２４によって燃料流量の値が修正され、このように修正された燃料流量の値で当該データセットの燃料の値をもとの値と差替えて、再度性能計算が実行される。この燃料流量の修正とこれを反映した性能計算は、前記判定工程２３の判定結果が真になるまで繰り返される。

上述した前記判定工程２３における、予め定められた許容基準とは、例えば、発電出力と排気温度のそれぞれの実測値と計算値の、偏差の割合ないし偏差の絶対値が、予め定められた上下限値内にあるかどうかである。あるいはさらに、これらの偏差の割合の和が予め定められた上下限値内にあるかどうかも判定するとなおよい。この判定工程２３の基準に関しては、これらに限らず、発電出力と排気温度のそれぞれの実測値と計算値の偏差の大小について、予め定めた基準によって許容範囲かどうか判定するものであればよい。

前記修正工程２４における燃料流量の修正方法の一例としては、発電出力と排気温度の計算値の、実測値に対する偏差をエネルギー換算し（ここで排気温度については排ガスのエンタルピーとして計算する）、これらの偏差の総和ΔＨ＿ｅｒと燃料入熱量ＱＦを使って、次の式（９）のように修正係数γを計算し、これを式（７）に入力することで、燃料流量の修正値を決めることができる。ここでｋは繰返し計算の回数や精度を調整するための係数である。

γ ＝ １ − ｋ×ΔＨ＿ｅｒ／ＱＦ …式（９）
なお、ここでγは「修正」係数として説明したが、これはこの係数が、前述の式（７）、式（８）で述べた「補正」係数と意味が異なるためである。すなわち、前述の式（７）、式（８）で述べた「補正」係数は、予め決められた関数によって燃料流量の「補正値」を計算するのに使われる「補正」係数であった。これに対して、ここで述べる「修正」係数とは、上述したような収束計算によって計算値と実測値がよく一致するように修正しながら求められた燃料流量の「修正値」に関する係数である。そこで両者を区別するために、このように呼び分けている。

以上の燃料流量の修正方法は一例であり、この他の方法であっても、前記判定工程２３の判定基準を満たす値を探索するために燃料流量の修正量を生成するものであれば、公知の最適化計算や収束計算の方法でよい。あるいは、前記判定工程２３のように個別のデータセットに対して判定するのでなく、複数データセットに対して計算された発電出力と排気温度が実機実測値に対して有する偏差の平均値や最大値など、データセット全体での偏差を示す指標値が、予め定められた許容範囲内に収まるかどうかという許容基準で、前記式（７）、式（８）の関数の内容を設定すなわち、前述したように関数の数式や未定定数を決定、あるいはアルゴリズムを動作させるための設定値を決定する最適化ないし収束計算であってもよい。

このような燃料流量の修正とこれを反映した性能計算によって前記判定工程２３の判定結果が真になれば、上述のようにして計算された燃料流量の修正値と修正係数γの値が、出力工程２５によって出力され、データ蓄積工程２６に送られる。この際、各データセットの吸気温度、圧縮機圧力比、元の燃料流量の値が合わせて一緒に送られ、データ蓄積工程２６に蓄積される。

ここまでの工程２２〜２５までの処理は、複数データセット２０の全データセットに対して終了するまで、データセットの繰返し終了判定工程２７と、次のデータセットに進む工程２８によって繰り返しが制御される。

この全データセットに対する繰返しが終了すると、後半の工程２９では、以上の繰返しで前記データ蓄積工程２６に蓄積された、全データセットに対する前記補正係数γとこれに対応する吸気温度、圧縮機圧力比、燃料流量のデータに基づいて、フィッティングにより式（８）の関数が特定される。フィッティングの具体的方法は、最も単純には式（８）の第二式の線形関数に対する最小二乗法によって同式の係数と定数項が決定される。しかし、関数の数式は線形に限ることはなく、修正係数γが吸気温度、圧縮機圧力比、燃料流量の関数として表されたものであればよい。また、フィッティングの手法も線形最小二乗法に限らず、公知の最適化や回帰分析の手法で、これらの関数の未定定数を、計算値と実測値の偏差が小さくなるように決定するものであればよい。

またここでは、関数の数式説明をわかりやすくするため式（８）の左辺は修正係数としたが、このように修正係数を使わず、左辺を修正燃料流量にして、直接燃料流量を計算する形式にしてもよいことはいうまでもない。この場合の関数の設定手順は、修正係数を左辺として上で述べたのと同様の手順で実行される。
（補正関数設定回路の構成図）
図４に、以上に説明した処理フローを実現するための補正関数設定回路１２の構成図を示す。補正関数設定回路１２は、前記複数データセット２０を用いて性能計算回路７を起動させて前記燃料流量の修正値または修正係数を計算する逆算回路３１と、このように逆算回路３１で計算された前記燃料流量の修正値または修正係数を、各データセットの吸気温度、圧縮機圧力比、元の燃料流量の値と一緒に順次一時記憶して蓄積するデータ保持手段３２と、該データ保持手段に蓄積されたこれらの情報に基づいて、前記補正関数の設定情報２１を同定して出力する補正関数同定回路３３からなる。なお、ここで「回路」と記したものは、図６の診断システムの説明において既述したように、コンピュータプログラムまたは専用チップなどの電子的な演算回路を示す。以下、それぞれを順に説明する。

逆算回路３１には、少なくとも以下の３要素が実装されている。第１に、上述したように性能計算回路７を起動して前記燃料流量の修正値または修正係数を決定する計算を実行するために、ガスタービンの性能を計算する前記性能計算工程２２（性能計算回路７）を呼び出して実行させる機能と、計算結果を判定する前記判定工程２３と、判定結果が偽の場合に前記燃料流量を修正する前記修正工程２４からなる収束計算ループが実装されている。第２に、この収束計算によって計算された燃料流量の修正値または修正係数を、対応する各データセットの吸気温度、圧縮機圧力比、元の燃料流量の値とともに出力する前記出力工程２５が実装されている。第３に、この収束計算を複数データセット２０の全データセットに対して繰返し処理するループ制御回路（工程２７、２８）が実装されている。

性能計算回路７は、前記性能計算工程２２の性能計算を実行する実装回路であり、実行される計算内容は性能診断について説明した図６で述べたように、吸気温度、圧縮機圧力比、燃料流量を含む入力情報に基づいて、性能を計算し、発電出力と排気温度を含む計算結果情報を出力する。

データ保持手段３２は、前記データ蓄積工程２６を実現する手段であり、前記逆算回路３１によって計算された全データセットに対する前記補正係数γの値が、対応する各データセットの吸気温度、圧縮機圧力比、元の燃料流量の値とともに格納され蓄積される。本データ保持手段３２は、典型的にはコンピュータのメモリなどの一時記憶領域として実現されるが、これに限らず、情報を電子的に格納する媒体であればよい。

補正関数同定回路３３は、補正関数の設定情報を同定する工程２９の演算手順を実装した回路である。式（８）の関数の未定定数についての具体的な設定値が工程２９の演算により決定され、補正関数設定情報２１として出力される。
（設定済み補正回路を用いた性能診断の実施方法）
本発明の性能診断システムでは、以上で述べたようにして補正関数設定回路１２で決定された補正関数設定情報２１（図４）が、前述の補正回路６（図１）に設定されて利用される。その診断の処理フローを図５に示す。

処理フローは、実機の吸気温度、圧縮機圧力比、燃料流量、発電出力、排気温度の実測値（以下、実機データセット４０と呼ぶ）を入力とし、前述の図３、図４で述べたように設定された補正回路６を使って燃料流量の補正値を計算する補正計算工程４１と、該補正計算工程４１で計算された燃料流量の前記補正値を入力に使ってガスタービンの性能を計算する性能計算工程２２と、該性能計算工程の計算結果を評価する評価処理工程４３と、該評価工程４３から出力された情報を蓄積するデータ蓄積工程４４と、前記評価処理工程４３からの出力情報または前記データ蓄積工程４４に蓄積された情報を外部出力する結果出力工程４５と、これら工程４１〜工程４５の処理を実機データセットに対して繰り返し処理するための終了判定工程４６と、次のデータセットに進む工程４７と、以上のデータセット全体に対する評価結果を出力する結果出力工程４８からなる。以下、各工程について説明する。

補正計算工程４１では、実機データセットのうち、前述の図３、図４で説明したように設定された補正回路６、すなわち、前述の式（７）（８）の関数の未定定数の値が決定された状態になった燃料流量の補正回路６に対して、実機の吸気温度、圧縮機圧力比、燃料流量の実測値が入力され、該燃料流量の補正値が計算される。

続く性能計算工程２２では、実機データセットの性能計算用の入力データ（吸気温度、圧縮機圧力比、燃料流量）のうち、燃料流量の値が前記補正計算工程４１で補正された値によって差替えられて入力される。具体的には、この入力に基づいて既述の性能計算回路７が実行されてガスタービンの性能が計算され、発電出力と排気温度を含む項目が計算される。なお、この計算結果は前述の図１の計算結果情報８となる。

このように性能計算工程２２で計算された発電出力と排気温度の計算値が、評価処理工程４３でそれぞれの実測値、すなわち、前記実機データセット４０のうちの、評価用の発電出力と排気温度の実測値と、比較評価される。この評価の処理は具体的には、既述の評価回路９（図１、図６）によって実行される。

前記評価処理工程４３で評価された結果の情報、あるいは評価に供された情報すなわち、実機データセット４０の各項目の情報と、性能計算工程２２の計算結果情報８（図１）は、データ蓄積工程４４に蓄積される。該データ蓄積工程４４は、具体的には既述の格納手段１０（図１、図６）と同様のハードウェアで実現される。

こうして前記評価処理工程４３で評価された結果情報、または前記データ蓄積工程４４に蓄積された情報は、続く結果出力工程４５によって外部出力される。外部出力は具体的には既述の出力手段１１（図１、図６）によって実行される。

以上に述べたような、補正計算工程４１ならびに性能計算工程２２の計算と、これに基づく計算値と実測値の評価処理(評価処理工程４３)、データ蓄積(データ蓄積工程４４)、結果出力（結果出力工程４５）は、実機データセット４０の一定期間のデータまたは全体について繰り返される。この繰返しは、終了判定工程４６で真すなわち終了と判定されるまで、次のデータの処理に進む工程４７が繰り返されることによって制御される。ここでの終了判定は、前記実機データセット４０の全てのデータに対して終了しているか、あるいは、診断を中止する指令がユーザ等によって入力されたか、などの診断のための計算を続行するかどうかの判定である。

こうして、終了判定されるまで繰り返された結果、前記データ蓄積工程４４に蓄積された評価結果情報は、結果出力工程４８で外部出力される。外部出力は具体的には既述の出力手段１１（図１、図６）によって実行される。

なお、結果の出力について、上記結果出力工程４５と４８の組み合わせは一例である。両者の違いは大まかには出力するタイミングの違いであり、扱っているデータの量が実機データセットの全体か、一部分のデータセットかということを別にすると、出力内容は同様であり、出力するタイミングと内容はこれらのいずれにあってもよい。例えば、バッチ処理によって特定期間中の全ての実機データセットに対して繰返し処理する場合、結果出力は結果出力工程４８のタイミングで実施される。また、例えば、一日に一回というように定期的に自動起動して、前回起動時との間のデータセットに対して繰返し処理する場合も同様である。逆に、例えば実機計測システムとオンライン接続されており、リアルタイムで評価する場合は、結果出力は、特に警報などは結果出力工程４５のタイミングで出力される。
（解析事例）
発明方式による性能評価精度の改善効果を検証するため、発明のモデル調整方式を実機の性能評価に適用して精度を評価した。表１に比較評価したモデル調整方式を示す。従来方式の例として一定値の補正係数を乗じて補正する単純な方式（定率型）を示した。また、発明方式については、影響度の低い因子を除いて、吸気温度だけで補正するようにした方式（簡略型）も示した。発明方式で全ての因子を考慮したケースは、簡略型と区別するため、ここでは完全型と呼ぶ。実機の運転データを対象に各方式でモデル調整し、この結果を使って図５の手順で性能計算した場合の、発電出力と排気温度の計算値と実測値の誤差の最大値（絶対値）、平均値（絶対値）、標準偏差を評価した。

この際、定率型のモデル調整は、具体的には前述の式（７）の補正係数γを評価期間全体で誤差が最も小さくなるような一定値に固定して実施した。簡略型のモデル調整は、前述の式（８）で圧縮機圧力比と燃料流量の項を除外したものをモデル式とし、このほかは完全型と同様の手順で実施した。完全型のモデル調整は図３で述べた手順で実施した。

結果の誤差率の推移を図７〜図９に、また、誤差の統計値を表２から表４に示す。結果をみると、誤差率は定率型では、図７に示すように発電出力と排気温度ともに最大２％程度を超えたのに対して、簡略型では図８に示すように概ね１．５％以内に縮小し、完全型では図９に示すようにほぼ１％以内に縮小した。従来の定率型に対して、発明方式の簡略型は、表２および表３に示すように排気温度の誤差は殆ど改善されなかったが、発電出力の誤差を平均値と標準偏差で２／３程度に低減できた。簡略型に対して完全型ではさらに、表３と表４に示すように、発電出力と排気温度ともに誤差は最大値がほぼ半減し、平均値と標準偏差もさらに２／３程度に低減した。

以上のように、簡略型でも従来型よりも精度を向上できるが、完全型ではさらに改善できることがわかった。簡略型と完全型の違いは、発電出力の誤差改善の程度の差だけでなく、排気温度の誤差を低減できるかどうかという差として現れた。これは完全型では、熱収支の偏差を補正するように摘出した因子を全て使っていることによる効果といえる。

なお、簡略型には吸気温度だけで燃料流量を補正した例を示したが、吸気温度と圧縮機圧力比、または吸気温度と燃料流量で、燃料流量を補正するようにすれば、より精度が向上し、完全型に近づくものと推定される。

（発明方式の代替案との比較）
本発明は以上で述べたように、入出熱収支の偏差を補正するために、燃料流量を補正することを提案している。これに対して、入出熱収支の偏差の主原因は排気温度の計測値と理論値のずれにあると考えられることから、排気温度を補正するアプローチもありうる。それにもかかわらず、燃料流量を補正するのは以下に述べる実用上の長所があるためである。

すなわち、排気温度を補正するとした場合、先の図６の従来方式で実施すると、発電出力と排気温度の誤差のうち、排気温度は補正されるが、発電出力は補正されずに誤差がそのまま残る（図７の解析事例の場合、発電出力の誤差の値はそのまま変化しない）。この発電出力の誤差を補正するために考慮すべき要因としては、入出熱収支以外にもガスタービン性能計算回路での性能特性の実機とのずれなど様々な要因がある。これを補正するには、ガスタービン性能特性に基づく発電出力の補正回路を新たに別に設ける必要があり、このために、ガスタービン性能計算モデルにおける各要素機器の性能特性マップの実機とのずれなど、複雑で多岐にわたる要因を検討する必要がある。さらに、発電出力を補正すると、排気温度にも影響が出るため、発電出力の補正と、排気温度の補正の間を適切に関係づけてモデル化し、両補正回路に組み合わせる必要がある。これらは全て、補正という目的に対しては複雑すぎ、システムが複雑化し、運用するコストと手間が増加する要因となる。

これに対して、本発明の燃料流量の補正では、補正された燃料流量が性能計算回路に入力されることにより、計算される排気温度だけでなく発電出力も変化する。このため、単一の補正回路６（式（７）（８））を補正関数設定回路１２で設定するだけで、発電出力と排気温度の両方の誤差が小さくなるようにモデル調整できる。発明の方式では、入出熱量の偏差やガスタービン性能特性の実機とのずれなど多岐にわたる誤差要因に対して、このように燃料流量という１点によって補正できるように企図されている。このため、性能診断の中心部分である性能計算回路７の本体の改造や調整を最小限に抑え、診断システムの調整や運用が容易かつ効率的に実施できる。

本発明では、上述したような実機と計算のずれに関わる複雑な諸要因の調整を、吸気温度と圧縮機圧力比と燃料流量という３つ入力項目により燃料流量を補正する（補正回路６）だけで実現できる。またこの補正回路６を設定するための補正関数設定回路１２も、上記３項目に発電出力と排気温度を加えた５項目のデータだけで実行できる。上述した複雑な誤差要因に対して、このように数少ない項目のデータで調整できるのは、前述したようにガスタービンサイクル全体について包括的に熱力学的に分析し、さらに単純化により重要因子を抽出する工夫をしたこと（式（１）〜（６））によるものである。

本発明のガスタービン性能診断システムの構成図。 補正関数設定回路の機能を示す入出力図。 補正関数設定回路の処理フロー図。 補正関数設定回路の構成図。 設定済み補正回路を用いた性能診断の処理フロー図。 従来の診断システムの構成図。 従来方式による性能評価事例（定率型）を示す図。 発明方式による性能評価事例（簡略型）を示す図。 発明方式による性能評価事例（完全型）を示す図。

符号の説明

１…実機計測値情報、４…ガスタービン診断システム、６…補正回路、７…性能計算回路、８…計算結果情報、１２…補正関数設定回路、３１…逆算回路、３３…補正関数同定回路、４１…補正計算工程。

Claims (9)

1. 実機で計測された吸気温度と圧縮機圧力比と燃料流量の入力を受けて、ガスタービンの発電出力と排気温度を計算して出力する性能計算回路と、該性能計算回路から出力された発電出力と排気温度の実機計測値からの偏差に基づいて性能を評価する評価回路を備えたガスタービン性能診断システムであって、
   前記吸気温度と圧縮機圧力比と燃料流量の実機計測値に基づいて、該燃料流量の値の補正値を計算し、かかる計算で得られた補正燃料流量を前記性能計算回路への入力として設定する補正回路を備えたことを特徴とするガスタービンの性能診断システム。
2. 前記補正回路が、吸気温度と圧縮機圧力比と燃料流量の実機計測値に基づいて、予め定められた関数関係に従って、該燃料流量の値を補正する計算を実行し、実機で計測された燃料流量の値を、補正された燃料流量の値で差替えて前記性能計算回路へ入力するものであることを特徴とする請求項１に記載のガスタービンの性能診断システム。
3. 前記性能計算回路に、実機での吸気温度と圧縮機圧力比と燃料流量の計測値についての一定期間の時系列データまたは運転条件を複数に変えたデータよりなる入力用複数データセットを入力して計算出力された発電出力と排気温度の、実機計測値に対する偏差が、予め定められた許容基準内になるように、前記補正回路の前記関数関係を設定する補正関数設定回路を備えたことを特徴とする請求項２に記載のガスタービンの性能診断システム。
4. 前記補正関数設定回路が、
   前記入力用複数データセットの各データセットのデータを前記性能計算回路入力して計算出力された発電出力と排気温度の値の実測値に対する偏差が、予め定められた上下限範囲内になるように、各データセットの燃料流量のデータを修正する逆算回路と、
   該逆算回路によって修正された燃料流量のデータを、当該データセットの吸気温度と圧縮機圧力比と燃料流量のデータとともに保持するデータ保持手段と、
   該データ保持手段に保持された複数のデータセットの修正燃料流量と吸気温度と圧縮機圧力比と燃料流量のデータに基づいて、修正燃料流量を吸気温度と圧縮機圧力比と燃料流量の関数としてフィッティングする補正関数同定回路を具備することを特徴とする請求項３に記載のガスタービンの性能診断システム。
5. 実機で計測された吸気温度と圧縮機圧力比と燃料流量の入力を受けて、ガスタービンの発電出力と排気温度を計算して出力する性能計算回路と、該性能計算回路から出力された発電出力と排気温度の実機計測値からの偏差に基づいて性能を評価する評価回路を備えたガスタービン性能診断システムであって、
   前記吸気温度の実機計測値、吸気温度と圧縮機圧力比の実機計測値、または吸気温度と燃料流量の実機計測値に基づいて、該燃料流量の値の補正値を計算し、かかる計算で得られた補正燃料流量を前記性能計算回路への入力として設定する補正回路を備えたことを特徴とするガスタービンの性能診断システム。
6. 前記補正回路が、吸気温度の実機計測値、または吸気温度と圧縮機圧力比の実機計測値、または吸気温度と燃料流量の実機計測値に基づいて、予め定められた関数関係に従って、該燃料流量の値を補正する計算を実行し、かかる計算で得られた補正燃料流量を前記性能計算回路への入力として設定するものであることを特徴とする請求項５に記載のガスタービンの性能診断システム。
7. ガスタービンの性能計算回路に実機で計測された吸気温度と圧縮機圧力比と燃料流量を入力して発電出力と排気温度を計算し、該計算で得られた発電出力と排気温度が実機での計測値に対して有する偏差に基づいて性能を評価するガスタービンの性能診断方法であって、
   前記吸気温度と圧縮機圧力比と燃料流量の実機計測値に基づいて、予め定められた関数関係に従って、該燃料流量の値を補正し、かかる計算で得られた補正燃料流量を前記性能計算回路に入力することを特徴とするガスタービンの性能診断方法。
8. 燃料流量を補正するための予め定められた前記関数関係が、前記性能計算回路に、実機の吸気温度と圧縮機圧力比と燃料流量についての一定期間の時系列データまたは運転条件を複数に変えたデータよりなる入力用複数データセットを入力して計算出力された発電出力と排気温度が、実機計測値に対して有する偏差が、予め定められた許容基準内になるように予め設定された関数関係であって、吸気温度と圧縮機圧力比と燃料流量を入力として、該燃料流量の値の補正値を出力するものであることを特徴とする請求項７に記載のガスタービンの性能診断方法。
9. ガスタービンの性能計算回路に実機で計測された吸気温度と圧縮機圧力比と燃料流量を入力して発電出力と排気温度を計算し、該計算で得られた発電出力と排気温度が実機での計測値に対して有する偏差に基づいて性能を評価するガスタービンの性能診断方法であって、
   前記吸気温度の実機計測値、または吸気温度と圧縮機圧力比の実機計測値、または吸気温度と燃料流量の実機計測値に基づいて、予め定められた関数関係に従って、該燃料流量の値を補正し、かかる計算で得られた補正燃料流量を前記性能計算回路に入力することを特徴とするガスタービンの性能診断方法。

Images