(57)【要約】      （修正有）【課題】エンジンの制御器に組み込み可能であり且つオンラインでリアルタイムにエンジンの計測不能項目を計測可能項目から推定し、その推定結果を用いて適切なエンジン制御を実現するガスタービンエンジンの性能推定システムを提供する。【解決手段】エンジンのある作動点で最適に設計されたカルマンフィルタ・ゲインＫと、実機エンジンを忠実にモデル化した非線形ダイナミックシミュレーション・モデルとから成る一定ゲイン拡張カルマンフィルタ(CGEKF)を利用して、観測変数と観測変数推定値との間に差が生じる場合はモデルの状態変数の一部分を成すチューニングパラメータをその差が最小となるように調整することにより、その非線形ダイナミックシミュレーション・モデルを常に実機に忠実なモデルとする。【選択図】図１
2009068359.tif 000002
【特許請求の範囲】
【請求項１】
2009068359.tif 000025
【請求項２】
2009068359.tif 000026
【請求項３】
2009068359.tif 000027
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ガスタービンエンジンの性能推定システム、特にエンジンの制御器に組み込み可能であり且つオンラインでリアルタイムにエンジンの計測不能項目を計測可能項目から推定し、その推定結果を用いて適切なエンジン制御を実現するガスタービンエンジンの性能推定システムに関するものである。
【背景技術】
【０００２】
航空用ガスタービンエンジンは航空機に搭載されて使用される性質上、運転中に発生する推力などの性能に関するパラメータの多くは直接計測することが難しい。また、エンジンの内部の温度・圧力といったエンジンの状態量をあらわすパラメータも安全上その他の理由により計測できない箇所が多い。そのため、通常、エンジンの制御においては本来制御を行いたいが計測することができない推力等そのものを制御するのではなく、エンジンの回転数等の計測可能項目を目安に間接的に制御を行ったり、監視計測などにおいて重要な部位の温度・圧力は計測できないが、他の計測できる箇所の温度・圧力等から両者の相関等により推測したりすることが一般的である。このような方法では、経年劣化や故障・不具合等による性能の変化があり、その結果、回転数等の計測可能項目と推力等の計測できないパラメータとの相関関係がそれまでの関係から変化したとしても、制御の方法を適切に調整することはできない上に、重要部位の温度・圧力の推定値が実際の計測値からズレることにもなる。
上記のように経年劣化、損傷等によるエンジンの性能変化があった場合に、エンジン制御系を適切に調整し最適な制御を行うためには、常にエンジンの性能状態を監視し、性能変化を推定する仕組みが必要となる。これまでガスパスアナリシスと呼ばれるエンジンのガス通過要素の性能変化を同定するコンディションモニタリングの手法や、カルマンフィルタを利用したエンジンの性能推定の方法等が実用化されている。
しかしこれらの手法はいずれも線形理論に基づいた手法であり、航空用エンジンのような非線形システムには適用できないため、通常、線形区分法によりあらかじめシステム・マトリクス等の必要データを計算しデータテーブルを作成しておき、実行時にはこのデータテーブルを内挿して使用する手法を採用している。このデータテーブルは、制御変数、状態変数、計測変数の数、更には飛行高度・飛行マッハ数などの異なるエンジン作動点の近似点数によって増大し、航空用エンジンの作動範囲を網羅しようとすると莫大な点数のデータが必要となり実用的でない。また、線形システム理論に基づいたカルマンフィルタの欠点を補うため、実機エンジンに忠実な非線形ダイナミックシミュレーション・モデル等を利用し、運転条件が変わる都度カルマンフィルタ・ゲインを再計算する拡張カルマンフィルタの手法が知られている(例えば、特許文献１を参照。)。この拡張カルマンフィルタ手法を用いた場合、莫大な点数のデ一タを用意しておく必要はないが、計算負荷の高いカルマンフィルタ・ゲインの計算を都度行う必要があり計算能力の高い計算機器が必要となるが、信頼性が最優先される航空エンジン組み込み機器等においてはこのような計算能力の高い機器を用いるのは現実的でない。
【０００３】
【特許文献１】特開２００５−２４８９４６号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
上述した通り、拡張カルマンフィルタを利用して非計測又は非感知パラメータを推定することによりエンジンの状態を推定する従来の制御システムは、計算負荷の高いカルマンフィルタ・ゲインの計算を都度行う必要があるため、演算処理能力の高いコンピュータが必要となる。
しかしながら、信頼性が最優先される航空エンジン組み込み機器においては、このような計算能力の高い機器を用いるのは現実的でないという問題がある。
そこで、本発明は、上記実情に鑑み創案されたものであって、エンジンの制御器に組み込み可能であり且つオンラインでリアルタイムにエンジンの計測不能項目を計測可能項目から推定し、その推定結果を用いて適切なエンジン制御を実現するガスタービンエンジンの性能推定システムを提供することを目的とする。
【課題を解決するための手段】
【０００５】
上記目的を達成するために請求項１に記載のガスタービンエンジンの性能推定システムでは、
2009068359.tif 000003
上記ガスタービンエンジンの性能推定システムは、観測変数(ym)と観測変数の推定値
2009068359.tif 000004
との差が最小となるようにカルマンフィルタ・ゲイン(Ｋ)をエンジンの動作中にリアルタイムに更新するのではなく、あるエンジンの作動点でその差が最小となるように設計された一定カルマンフィルタ・ゲインを、実機エンジンをモデル化した非線形ダイナミックシミュレーション・モデルに組み合わせ、そして観測変数とその推定値との間に差が生じる場合はその差が最小となるようにモデルの状態変数の一部分を成すチューニングパラメータによって非線形ダイナミックシミュレーション・モデルを逐次チューニングすることによって常に実機に忠実なモデルとする。従って、上記非線形ダイナミックシミュレーション・モデルを必要に応じて逐次チューニングすることによって、エンジンの動作中に計算負荷の高いカルマンフィルタ・ゲインを再計算・再設計する必要がなくなり、その結果、計算機に対する負荷が好適に軽減され上記性能システムをエンジンの制御器に組み込むことが可能となる。また、非線形ダイナミックシミュレーション・モデルを逐次チューニングすることによって、計測可能項目と計測不能項目との間の相関関係を変化させる経年劣化や損傷等の影響が十分に反映された高精度な非線形ダイナミックシミュレーション・モデルとすることが可能となる。従って、この高精度な非線形ダイナミックシミュレーション・モデルから実機では計測不能な項目（観測不能変数）を高精度に推定することによって上記の目的を達成することが出来るようになる。つまり、エンジンのある作動点で最適に設計されたカルマンフィルタ・ゲインは、他の全作動点に対しても有効であり、そして観測変数と観測変数の推定値との間に差が生じる場合はモデルの状態変数の一部分を成すチューニングパラメータを使用してその差が最小となるように非線形ダイナミックシミュレーション・モデルを調整することにより、常に実機に忠実なモデルとする。その結果、上記性能推定システムは、計算負荷の高いカルマンフィルタ・ゲインの更新を都度行う必要がなくなり、計算機に対する負荷を好適に軽減すると共にエンジンの制御器に組み込むことが可能となる。
【０００６】
請求項２に記載のガスタービンエンジンの性能推定システムでは、
2009068359.tif 000005
上記ガスタービンタービンエンジンの性能推定システムでは、上記構成とすることにより、オンラインでリアルタイムにエンジンの計測不能項目を計測可能項目から好適に推定することが出来るようになる。
【０００７】
請求項３に記載のガスタービンエンジンの性能推定システムでは、
2009068359.tif 000006
上記ガスタービンエンジンの性能推定システムでは、上記構成とすることにより、これまで実機においては不可能であった推力や燃料消費率等の観測不能変数の推定値を入力または制御目標としてフィードバック制御器に組み込むことが可能となり、その結果、観測することができないパラメータの変化をも考慮した高精度なエンジンの制御が可能となる。
【発明の効果】
【０００８】
ジェットエンジンを含む航空用機器に使用される電子機器類は本質的に安全性・確実性が求められ、それゆえ実績が十分ある電子デバイス等を使用して設計されることが多く、またシステムの信頼性を確保するために回路自体はなるべくシンプルで部品点数を少なくするように設計されることが多い。それゆえ、計算処理能力の高いプロセッサ等を使うことが難しく、航空機・ジェットエンジンに組み込み可能な実用的なエンジン性能推定システムを開発することは困難であった。本発明による一定ゲイン拡張カルマンフィルタを利用したエンジン性能推定システムでは、推定に必要なカルマンフィルタ・ゲインの計算を予めあるエンジン作動点で行った結果を用いることにより、エンジン動作中にカルマンフィルタ・ゲインを再計算・再設計する必要がなく、計算処理負荷が非常に少ないエンジン牲能推定システムを実現できる。
また、計測項目についても非計測項目として他の計測パラメータから推定し、両者を比較することによりセンサ計測値異常などを検知することが可能となる。
【発明を実施するための最良の形態】
【０００９】
以下、図に示す実施の形態により本発明をさらに詳細に説明する。
【００１０】
本発明の実施形態は、ジェットエンジン・ガスタービン等のシステムにおいて、製造上許容されるエンジン毎の性能の差や、エンジン使用による経年変化や損傷による性能劣化を、エンジンの非線形ダイナミックシミュレーション・モデルと推定フィルタにより常に推定・検出・監視するエンジン性能推定システムおよび方法を提供する。更に本発明によるシステムおよび方法は、観測不能変数をも目標値としたエンジン制御システムを提供する。
【００１１】
図１は、本発明に係るガスタービンエンジンの性能推定システム１００を示す概要説明図である。エンジンは制御変数・環境変数ｕによって動作状態を変更することができる。このときエンジンのセンサの計測値を観測変数ymとし、制御変数・環境変数ｕおよび観測変数ymは一定ゲイン拡張カルマンフィルタ(CGEKF)に入力される。観測変数ymと観測変数推定値
2009068359.tif 000007
の差は、エンジンのある１又は複数の動作状態において最小になるように設計されたカルマンフィルタ・ゲインＫによって増幅され、一定ゲイン拡張カルマンフィルタ(CGEKF)内の非線形ダイナミックシミュレーション・モデルに入力される。制御変数・環境変数ｕもまた入力される。非線形ダイナミックシミュレーション・モデル内では、エンジンの状態変数とチューニング変数（パラメータ）とから成る拡張状態変数
2009068359.tif 000008
、観測変数推定値
2009068359.tif 000009
および観測不能変数推定値
2009068359.tif 000010
がそれぞれ計算される。一連の計算が繰り返されることにより、観測変数ymと観測変数推定値
2009068359.tif 000011
の差が最小となり、観測不能変数の最尤推定値が得られる。
【００１２】
図２は、従来の観測変数を制御目標値としたエンジン制御方法を示す説明図である。従来のエンジン制御では本来制御したい変数が観測不能であることが多く、他の相関のある観測変数ymをフィードバックしてエンジンを制御している。
【００１３】
図３は、本発明の実施形態による観測不能変数を目標値としたエンジン制御方法を示す説明図である。制御目標値に観測不能変数をも与えることが可能であり、エンジンの観測変数ymと制御変数uを一定ゲイン拡張カルマンフィルタ(CGEKF)に入力し、得られた観測不能変数推定値
2009068359.tif 000012
をフィードバックしてエンジンを制御することができる。
【００１４】
一定ゲイン拡張カルマンフィルタ（ＣＧＥＫＦ）
非線形ジェットエンジンのダイナミクスは、
【数（１）】
2009068359.tif 000013
で表せる。ここで、X：状態変数ベクトル、u：制御変数ベクトル、ym：計測可能変数ベクトル、yu：計測不能変数ベクトル、v：システムノイズベクトル、ｗ：計測ノイズベクトル、G：システムノイズ伝達関数、f(),gm(),gu()：関数ベクトルである。数(1)を任意作動点近傍で線形化すると、
【数（２）】
2009068359.tif 000014
となる。ここで、A,B,C,D：システム行列、Cu,Du：計測不能変数観測行列である。また、システムノイズ、計測ノイズの平均値および共分散を、
【数（３）】
2009068359.tif 000015
とする。状態変数、計測可能変数、計測不能変数の最尤推定値
2009068359.tif 000016
は、
【数（４）】
2009068359.tif 000017
と表せる。ここで、
【数（５）】
2009068359.tif 000018
であり、Kは線形カルマンフィルタ・ゲインである。これを非線形システムに拡張すると、数(4)に対応して、最尤推定値
2009068359.tif 000019
は、
【数（６）】
2009068359.tif 000020
となり、カルマンフィルタ・ゲイン
2009068359.tif 000021
は作動点の非線形関数である。作動点が移動する場合、数(5)を繰り返し解かなければならず、制御計算機への負荷が過剰になるため、代表的作動点において求められたカルマンフィルタ・ゲイン
2009068359.tif 000022
を全作動領域に適応することにする。これを一定ゲイン拡張カルマンフィルタ(CGEKF)と呼ぶ。
【００１５】
１軸ターボジェットエンジンのカルマンフィルタゲインの導出
１軸ターボジェットエンジンのエンジン状態変数Xeは、
Xe=(N,m3,u3,m4,u4,m5,u5)T    ・・・・（数(7))
ただし、N：ロータ回転速度、m3：コンプレッサ出口ボリュームの蓄積質量、u3：コンプレッサ出口ボリュームの蓄積内部エネルギー、m4：燃焼器ボリュームの蓄積質量、u4：燃焼器ボリュームの蓄積内部エネルギー、m5：ノズルボリュームの蓄積質量、u5：ノズルボリュームの蓄積内部エネルギー。
エンジン要素特性変化を同定するため、チューニング変数Xcは、
Xc=(Gc,ηc,Gt,ηt,An)T    ・・・・(数(8))
とする。ただし、Gc：コンプレッサ流量パラメータ、ηc：コンプレッサ効率パラメータ、Gt：タービン流量パラメータ、ηt：タービン効率パラメータ、An：ノズル面積パラメータ。カルマンフィルタは状態変数を推定するものであるため、
【数（９）】
2009068359.tif 000023
により、Xcを人為的な状態変数として組み込む。エンジン状態変数Xeとチューニング状態変数Xcを結合して、拡張状態変数Xを、
X=(XeT,XcT)T
=(N,m3,u3,m4,u4,m5,u5,Gc,ηc,Gt,ηt,An)T       ・・・・(数(10))
と定義する。CGEKFは、計測値ymと計測値の推定値
2009068359.tif 000024
に差がある時、それを０にするようにXeとXcとを変化させる。例えば、経年性能劣化や異物吸い込み等でエンジン要素特性が変化すると、実機エンジンとモデルエンジンに差異が生じ、計測値とモデルエンジンによる推定値に差がでることになるが、CGEKFはその差をなくす最も合理的なXcの変化を算出する。
制御変数uは、ここでは一変数のみで、
u=Wf         ・・・・(数(11))
ただし、Wf：燃料流量。計測可能変数ymは任意にとることができるが、ここでは、
ym=(N,P3,T3,P6,T6)T         ・・・・(数(12))
とする。ただし、P3：コンプレッサ出口圧力、T3：コンプレッサ出口温度、P6：タービン出口圧力、T6：タービン出口温度。計測不能変数yuは、制御やモニタリングに有用な変数として、
yu=(F,TIT,SFC,ηc,ηt,・・・)T        ・・・・(数(13))
とする。ただし、F：推力、TIT：タービン入口温度、SFC：燃料消費率、ηc：コンプレッサ効率、ηt：タービン効率。
各変数の次元は、次のようにまとめられる。
(a)エンジン状態変数 Xe=7
(b)同定したいエンジン変数 Xc=5
(c)計測可能変数 ym=5
(d)制御変数 u=1
カルマンフィルタが存在する必要条件は、モデルが可観測(Observable)であることである。ここで考えている１軸ターボジェットエンジンの場合、
(Xcの次元）≦(ymの次元)         ・・・・(数(14))
が必要条件となる。すなわち、同定したい変数の数は、計測点数に等しいか、少なくなければならない。上の例では、両者とも５であるので、この必要条件は満足している。
【００１６】
エンジンの任意の作動点におけるカルマンフィルタ・ゲインは次の手順で計算される。
（１）図１のCGEKFモデルを作成し、代表的作動点における拡張システム行列A,B,C,Dを導出する。
（２）数値計算上の問題を避けるため、拡張状態変数、制御変数、計測可能変数、計測不能変数のスケーリング値を設定し、拡張システム行列A,B,C,Dを規準化する。
（３）システムノイズ共分散行列Q、計測ノイズ共分散行列Rおよびシステムノイズ伝達行列Gを仮定する。Rは計測信号の統計的性質であり、設定することはできるが、Q,Gを明確に設定することは困難である。従って、Q,Gはカルマンフィルタ設計時のチューニングのための自由パラメータとする。ノイズ特性、応答速度等を勘案して設定する。
（４）行列(A,C)でシステムの可観測性（Observability)を確認し、カルマンフィルタ・ゲインKを求める。
【産業上の利用可能性】
【００１７】
本発明のガスタービンエンジンの性能推定システムは、航空用ジェットエンジン、ロケットエンジンのターボポンプ、および産業用・発電用ガスタービンシステムの性能推定、センサ・アクチュエータ等の健全性モニタ・異常診断に対し好適に適用することが可能である。
【図面の簡単な説明】
【００１８】
【図１】本発明に係るガスタービンエンジンの性能推定システムを示す概要説明図である。
【図２】従来の観測変数を目標値としたエンジン制御方法を示す説明図である。
【図３】本発明の実施形態による観測不能変数を目標値としたエンジン制御方法を示す説明図である。
【符号の説明】
【００１９】
１００   ガスタービンエンジンの性能推定システム
【図１】
2009068359.tif 000028
【図２】
2009068359.tif 000029
【図３】
2009068359.tif 000030