(57)【特許請求の範囲】
【請求項１】
あらかじめ設定された回転速度で回転するロータと、固定羽根及び可動羽根の翼列とを有するタービンエンジンのコンプレッサにおける、航空力学的流体流出の不安定性の早期検知方法であって、
固定羽根と可動羽根との翼列の間でコンプレッサの円周上に配置された、角位置θの少なくとも１つのセンサによって与えられる信号の獲得（２０）を行い、この獲得は、あらかじめ設定された時間内に行われ、
獲得された信号を正規化し（３０）、信号の正規化は、信号の平均値によって構成される静的成分を取り除き、信号の変動によって得られた結果を分けることからなる、信号の適応フィルタリングによって行なわれ、
正規化された信号が、センサの前を連続する２つの羽根が通過する時間Ｔａより少なくとも長い時間Ｔｚ中、ゼロの値を通らないことからなる特性を有するような事象を、正規化された信号から検知し（４０）、
検知された事象全体を分析し（５０）、
事象の分析の結果に応じて、コンプレッサにおける航空力学的流体流出の不安定性が目前に迫っているかどうかについての決定を下す（６０）
ことを含むことを特徴とする、航空力学的流体流出の不安定性の早期検知方法。
【請求項２】
位置θのセンサによる、瞬間ｔにおける事象の検知結果が、時間Ｔｚ中にゼロを通らないという特性が確認されるか否かに応じて、論理値１または０をとる関数Ｅ（ｔ，θ）によって表わされることを特徴とする、請求項１に記載の航空力学的流体流出の不安定性の早期検知方法。
【請求項３】
事象の分析が、
事象の検知頻度の値を得る（５２）ために、センサの前を連続する２つの羽根が通過する時間Ｔａより長いあらかじめ設定された時間中に、センサによって検知された事象を数え上げ（５１）、
不安定性が目前に迫っているかどうかの少なくとも第一の指示を得るために、あらかじめ設定された第一の閾値と得られた検知頻度の値とを比較する（５３）ことを含むことを特徴とする、請求項１に記載の航空力学的流体流出の不安定性の早期検知方法。
【請求項４】
検知された事象の分析が、
コンプレッサの円周上に等配分された少なくとも２つのセンサによって与えられる正規化された信号から、検知された事象の空間−時間的位置を読取り（５５）、
時間／方位（ｔ，θ）の表現空間において、コンプレッサのロータの回転速度に等しい値とゼロの値との間に含まれる伝播速度に対応する事象の伝播方向を求め（５６）、
求められた各伝播方向について、検知された各事象に対して、求められた伝搬方向において検知された事象の空間的伝播の存在に関する蓋然性のレベルを割当て（５７）、
不安定性が目前に迫っているかどうかの少なくとも一つの指示を得るために、あらかじめ設定された第二の閾値と各蓋然性のレベルとを比較する（５８）ことを含み、
角速度ｖで伝播する、角度位置θで時間ｔにおいて検知された事象Ｅ０についての蓋然性のレベルＪ０（ｔ，θ，ｖ）が、角速度ｖに対応する伝播方向において、当該の事象Ｅ０より先に検知された事象の数の関数であることを特徴とする、請求項１に記載の航空力学的流体流出の不安定性の早期検知方法。
【請求項５】
検知された事象の分析が、さらに、
コンプレッサの円周上に等配分された少なくとも２つのセンサによって与えられる正規化された信号から、検知された事象の空間−時間的位置を読取り（５５）、
時間／方位（ｔ，θ）の表現空間において、コンプレッサのロータの回転速度に等しい値とゼロの値との間に含まれる伝播速度に対応する事象の伝播方向を求め（５６）、
求められた各伝播方向について、検知された各事象に対して、求められた伝搬方向において検知された事象の空間的伝播の存在に関する蓋然性のレベルを割当て（５７）、
不安定性が目前に迫っているかどうかの少なくとも第二の指示を得るために、あらかじめ設定された第二の閾値と、各蓋然性のレベルとを比較する（５８）ことを含み、
角速度ｖで伝播する、角度位置θで時間ｔにおいて検知された事象Ｅ０についての蓋然性のレベルＪ０（ｔ，θ，ｖ）が、角速度ｖに対応する伝播方向において、当該の事象Ｅ０より先に検知された事象の数の関数であること特徴とする、請求項３に記載の航空力学的流体流出の不安定性の早期検知方法。
【請求項６】
検知された事象の分析が、
コンプレッサの円周上に等配分された少なくとも２つのセンサによって与えられる正規化された信号から、検知された事象の空間−時間的位置を読取り（５５）、
時間／方位（ｔ，θ）の表現空間において、コンプレッサのロータの回転速度に等しい値とゼロの値との間に含まれる伝播速度に対応する事象の伝播方向を求め（５６）、
求められた各伝播方向について、検知された各事象に対して、求められた伝搬方向において検知された事象の空間的伝播の存在に関する蓋然性のレベルを割当て（５７）、
不安定性が目前に迫っているかどうかの少なくとも一つの指示を得るために、あらかじめ設定された第二の閾値と各蓋然性のレベルとを比較する（５８）ことを含み、
角速度ｖで伝播する、角度位置θで時間ｔにおいて検知された事象Ｅ０についての蓋然性のレベルＪ０（ｔ，θ，ｖ）が、角速度ｖで伝播する、事象Ｅ０より先に検知された事象について計算された前の蓋然性のレベルＪｋ（ｔ，θ，ｖ）（ここでｋは正の整数である）の関数であることを特徴とする、請求項１に記載の航空力学的流体流出の不安定性の早期検知方法。
【請求項７】
検知された事象の分析が、さらに、
コンプレッサの円周上に等配分された少なくとも２つのセンサによって与えられる正規化された信号から、検知された事象の空間−時間的位置を読取り（５５）、
時間／方位（ｔ，θ）の表現空間において、コンプレッサのロータの回転速度に等しい値とゼロの値との間に含まれる伝播速度に対応する事象の伝播方向を求め（５６）、
求められた各伝播方向について、検知された各事象に対して、求められた伝搬方向において検知された事象の空間的伝播の存在に関する蓋然性のレベルを割当て（５７）、
不安定性が目前に迫っているかどうかの少なくとも第二の指示を得るために、あらかじめ設定された第二の閾値と、各蓋然性のレベルとを比較する（５８）ことを含み、
角速度ｖで伝播する、角度位置θで時間ｔにおいて検知された事象Ｅ０についての蓋然性のレベルＪ０（ｔ，θ，ｖ）が、角速度ｖで伝播する、事象Ｅ０より先に検知された事象について計算された前の蓋然性のレベルＪｋ（ｔ，θ，ｖ）（ここでｋは正の整数である）の関数であること特徴とする、請求項３に記載の航空力学的流体流出の不安定性の早期検知方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、タービンエンジンのコンプレッサにおける航空力学的不安定性の早期検知方法に関するものである。本発明は、あらゆるタイプのタービンエンジン、特に航空機のタービンジェットエンジンに適用される。
【０００２】
【従来の技術】
特定の作動条件においては、タービンエンジンのコンプレッサは、エンジンのエレメントの寿命の単なる短縮から、事故を招く恐れのあるエンジンの完全な燃焼中断にいたるまでの損傷を引き起こすような、不安定な作動モードに移行する恐れがある。
【０００３】
回転剥離（ｄｅｃｏｌｌｅｍｅｎｔ  ｔｏｕｒｎａｎｔ）及びサージの名で知られている航空力学的不安定性が、従来は、エンジンの作動条件の変化が、結果的に、流体の流量、出力温度または圧力のようなコンプレッサのパラメータのいずれか１つの急激な変化を招くときに発生する。
【０００４】
サージは、コンプレッサを介した流体流出の主に軸方向の振動によって特徴づけられる。
【０００５】
回転剥離は、剥離セル（ｃｅｌｌｕｌｅｓ  ｄｅ  ｄｅｃｏｌｌｅｍｅｎｔ）と呼ばれる局所化された１つまたは複数のゾーンを生み出す現象であり、この剥離セルが、コンプレッサの羽根の翼列における、コンプレッサの回転速度の一部分において円周方向に伝播する。
【０００６】
剥離セルの性質に応じて、流動経路の高さ全体に広がった剥離（英語でｆｕｌｌ  ｓｐａｎ  ｓｔａｌｌ）と、部分的剥離（英語でｐａｒｔ  ｓｐａｎ  ｓｔａｌｌ）とが区別される。部分的剥離は、同時に限られた数の羽根にしか影響を与えず、半径方向においては、コンプレッサの流動経路の一部分に限定される、局所的擾乱を構成する複数の回転するセルで構成されることもある。
【０００７】
逆に、全体に広がった剥離は、正常なゾーンと剥離したゾーンとの間により大きい振幅の変化をともない、羽根の翼列の半径方向及び角度方向のより大きな領域に影響を与える。
【０００８】
一般に、剥離セルの伝播速度は、全体に広がった剥離に対応する速度より、部分的剥離の場合の方が大きい。
【０００９】
剥離セルの伝播速度は、部分的剥離についてはコンプレッサのロータの回転速度のおよそ７０から８０％であるのに対して、全体に広がった剥離についてはおよそ４０から５０％である。
【００１０】
それぞれ局所的擾乱とモード擾乱と呼ばれる、２つのタイプの航空力学的不安定性の前兆が区別され、それら２つのタイプの前兆に対して、有効な検知方法を確立することが必要である。
【００１１】
局所的擾乱は、一般に、不安定性がコンプレッサ全体に広がる前に、コンプレッサの数十または数百回の回転において見られる。しかしながら、それらは、検知するのが難しい。なぜなら、それらは限られた数の羽根と、コンプレッサの流動経路の一部分とに限定されるからである。それらの振幅は小さく、挙動は大きく偶然に左右される。局所的擾乱は、ほぼ瞬間的に現われ、コンプレッサの不可逆的な全体の不安定性を引き起こすことなく、コンプレッサの回転の一部分の間に伝播する。さらに、これら局所的擾乱は、不安定性ゾーンに近づくと、コンプレッサの特定の段が、他の段より早く航空力学的負荷の限界に達するという状況の特徴的なものであり、その結果、センサが、好ましくは、これら局所的擾乱を検知することができるように、この段の近傍に置かれることが必要となる。
【００１２】
モード擾乱は同様に、不安定性より早く、数十または数百回のコンプレッサの回転において現われる。それらは、コンプレッサが、流量に応じた圧力特性曲線の最大値に達すると、コンプレッサ全体の流体流出の変動に対する動的レスポンスに直接関連する。これらのモード擾乱は、このようにして、コンプレッサの種々の段が、いっしょに航空力学的負荷の限界に達するといった状況の特徴的なものである。モード擾乱は、一定の瞬間にコンプレッサの円周全体に影響を与え、したがって、それらはセンサ全体で観測されることから、局所的擾乱より相対的に検知しやすい。
【００１３】
航空力学的不安定性の問題を解決するために、不安定性の発生を防ぐことができる予防用機械的システムを設置する方法が知られている。これらの機械的システムは、たとえば、可変予備回転タービン翼列、可変ステータ、およびコンプレッサの内部に取付けられた排出バルブとすることができる。しかしながら、これらの予防用システムによって、大きな作動上の制約と性能の限界が余儀なくされる。
【００１４】
さらに、サージ状況を検知し、そこから脱することができる修復システムが存在する。
【００１５】
最近の民間及び軍事用エンジンの大半に使用されている不安定性の能動的制御技術は、不安定性サイクルが生じたときにそれを検知するというものである。この検知は、次に、この不安定性ゾーンから最良の条件へ抜け出すために、コンプレッサまたはエンジンのパラメータの補正操作を開始させる。
【００１６】
しかしながら、この技術は、エンジンを損傷し、そのエンジンが搭載されている航空機の墜落を招くという危険を有する。
【００１７】
十分な安全性を保証しながらもエンジンの性能を改良するために、コンプレッサの作動点を、不安定な作動ゾーンの近くのゾーンにできるだけ長時間位置させることができる、また別の不安定性の能動的制御技術を使用する方法が知られている。
【００１８】
この不安定性の能動的制御技術は、コンプレッサが不安定性サイクルに入る前に、発生する航空力学的不安定性の前兆を検知することからなる。これらの前兆とは、たとえば、コンプレッサが不安定性サイクルに入る前に、コンプレッサの数回転でセンサの信号上に現われる、圧力、速度、温度の変化とすることができる。この検知が十分早期に行なわれる場合には、続いて、不安定性が実際に生じる前に、コンプレッサの作動点または幾何学的構成の補正操作を実施することができる。この不安定性の能動的制御技術は、非常に小さいサージ幅とともに作動することができるという利点を有する。なぜなら、この技術は、リアルタイムで、これら不安定性の予防を行うからである。しかしながら、不安定性の前兆の検知の現行の方法は、実施が難しく、あまり性能が良くない。なぜなら、それらの方法は、一般に、非常にノイズの多い環境において測定される信号の振幅について適用される、検知の閾値の調整に関連するからである。
【００１９】
特許ＷＯ９６／３４２０７は、航空力学的不安定性の前兆を検知することができ、コンプレッサの回転周波数における周波数信号のエネルギーを測定し、さらに、このエネルギーを、加速のときに正常なエンジンで得られる信号を観察することによって調整される経験的な閾値と比較することからなる、方法について記載している。
【００２０】
しかしながら、この方法では、コンプレッサの回転周波数における信号のエネルギーだけが分析され、その結果、コンプレッサの回転周波数とは別の周波数において発生する恐れのある、不安定性のあらゆる前兆を検知することはできない。
【００２１】
このことは、局所的擾乱の伝播速度が、多くの場合、コンプレッサの回転速度より小さいことから、とりわけ大きな支障となる。
【００２２】
【発明が解決しようとする課題】
本発明の目的は、高性能で、効果的で、実施が容易で、かつ、伝播速度に関係なくあらゆる擾乱を検知することができる、タービンエンジンのコンプレッサにおける航空力学的不安定性の新しい早期検知方法を決定することにある。
【００２３】
【課題を解決するための手段】
とりわけ、本発明は、固定羽根及び可動羽根の翼列の間でコンプレッサの円周上にほぼ均等に配置された、圧力または流出速度の非定常センサーによって与えられる信号上に非常に早期に現われる、局所的擾乱の検知方法に関するものである。
【００２４】
そのために、本方法は、不安定性の前兆の特徴的事象を検知し、コンプレッサの中で航空力学的不安定性が目前に迫っているかどうかを判断することができる、決定基準を適用することからなる。
【００２５】
決定は、不安定性の前兆の特徴的事象の検知頻度に基づいた基準に照らして下すことができる。
【００２６】
また決定は、センサ全体において検知された事象の空間−時間位置と、検知された局所的擾乱の伝播の存在の蓋然性のレベルの評価とに基づく基準に照らして下すことができる。さらに、決定は、連続的にまたは同時に適用される、先の２つの基準に照らして下すことができる。
【００２７】
本発明によれば、あらかじめ設定された回転速度で回転するロータと、固定羽根列及び可動羽根の翼列とを有するタービンエンジンのコンプレッサにおける、航空力学的流体流出の不安定性の早期検知方法は、
固定羽根と可動羽根との翼列の間でコンプレッサの円周上に配置された、角位置θの少なくとも一つのセンサによって与えられる信号の獲得を行い、この獲得は、あらかじめ設定された時間内に行われ、
獲得された信号を正規化し、
正規化された信号が、センサの前を連続する２つの羽根が通過する時間Ｔａよりも少なくとも長い時間Ｔｚ中、ゼロの値を通らないことからなる特性を有するような事象を、正規化された信号から検知し、
検知された事象全体を分析し、
事象の分析結果に応じて、コンプレッサにおける航空力学的流体流出の不安定性が目前に迫っているかどうかについて決定を下すことを含むことによって特徴づけられる。
【００２８】
本発明の第一の実施形態においては、検知された事象の分析は、
事象の検知頻度の値を得るために、センサの前を連続する２つの羽根が通過する時間Ｔａより長いあらかじめ設定された時間中に、センサによって検知された事象を数え上げ、
不安定性が目前に迫っているかどうかの少なくとも第１の指示を得るために、あらかじめ設定された閾値と得られた検知頻度の値とを比較することからなる。
【００２９】
本発明の第二の実施形態においては、検知された事象の分析は、
コンプレッサの円周上に等配分された少なくとも２つのセンサによって与えられる正規化された信号から、検知された事象の空間−時間的位置を読取り、
時間／方位の表現空間（ｔ，θ）において、コンプレッサのロータの回転速度に等しい値とゼロの値との間に含まれる、伝播速度に対応する事象の伝播方向を求め、
求められた各伝播方向について、検知された各事象に対して、求められた方向において検知された事象の空間的伝播の存在に関する蓋然性のレベルを割当て、
不安定性が目前に迫っているかどうかの少なくとも一つの指示を得るために、あらかじめ設定された第二の閾値と各蓋然性のレベルとを比較することからなる。
【００３０】
本発明の第三の実施形態においては、検知された事象の分析は、
事象の検知頻度の値を得るために、センサの前を連続する２つの羽根が通過する時間Ｔａより長いあらかじめ設定された時間中に、センサによって検知された事象を数え上げ、
不安定性が目前に迫っているかどうかの少なくとも第一の指示を得るために、得られた検知頻度の値と、あらかじめ設定された第一の閾値とを比較し、
コンプレッサの円周上に等配分された少なくとも２つのセンサによって与えられる正規化された信号から、検知された事象の空間−時間的位置を読取り、
時間／方位の表現空間（ｔ，θ）において、コンプレッサのロータの回転速度に等しい値とゼロの値との間に含まれる、伝播速度に対応する事象の伝播方向を求め、
求められた各伝播方向について、検知された各事象に対して、求められた方向において検知された事象の空間的伝播の存在に関連する蓋然性のレベルを割当て、
不安定性が目前に迫っているかどうかについて少なくとも第二の指示を得るために、あらかじめ設定された第二の閾値と、各蓋然性のレベルとを比較することからなる。
【００３１】
添付の図面を参照して、限定的でなく例示的なものとして与えられる以下の説明によって、本発明の他の特徴及び利点が明らかになるだろう。
【００３２】
【発明の実施の形態】
図１に示されているように、センサ１０は、羽根の翼の前縁の最も近くに位置するコンプレッサのロータの羽根の翼列１１の上流側に置かれ、剥離セルを観察する確率を最適化するために、コンプレッサの円周上にほぼ均等に配分される。使用されるセンサは、本質的に非定常である。すなわち、圧力または速度の急激な変化を読取ることが可能である。以下の説明は、圧力センサを使用する実施形態に限定される。速度センサによって観察される挙動は、反対方向の変化であるが類似している。
【００３３】
図２には、コンプレッサが正常に、したがって安定作動するときに、ロータの羽根の前縁の上流側に置かれたセンサによって観察される、圧力信号Ｐと羽根の翼列１１とが示されている。
【００３４】
この信号は、圧力センサの場合には、コンプレッサの各段によって行なわれる圧力の増大に対応する、静的成分Ｐｓｔａｔｉｑｕｅを有する。この静的成分を中心とする変化はより複雑である。ここには図示されていないスペクトル分析は、数多くの周波数の存在、特に回転速度の高調波、さらに、羽根の通過周波数への主要な寄与を示すことができる。
【００３５】
このことは、以下のように解釈される。ただ一つの羽根が固定されている翼列を対象にする場合、圧力センサ（固定されている）によって観察される信号に対するその局所的寄与は、各回転時にのみ見ることができるだろう。１回転の時間をＴｒとすると、この寄与は、Ｔｒ秒ごとに、したがって、１／Ｔｒに等しい回転周波数Ｆｒにおいて現われる。Ｎａ個の羽根が固定されている翼列の最も一般的ケースにおいては、システムが完璧であれば、各羽根は連続的に同一の寄与を行う。このようにして、圧力センサによって観察される信号は、１回転する間に、同一の寄与、すなわちＮａ／Ｔｒに等しい羽根の通過周波数Ｆａ（１／Ｔａ）への寄与を、Ｎａ回読取る。
【００３６】
この寄与は、単に、羽根の前縁の通過が、圧力の局所的増大とその後の減少を引き起こすことを考慮することによって、特徴づけることができる。
【００３７】
圧力信号を正規化すると、コンプレッサの正常作動に対応する正規化された信号は、定期的にゼロの値を通る。すなわち、連続する２つの羽根の間の通過時間に対応する圧力信号は、少なくとも一回の符号の変化を有する。
【００３８】
信号の正規化は、信号の平均値によって構成される静的成分を取り除き、信号の変動によって得られた結果を分けることからなる、信号の適応フィルタリングによって行なわれる。
【００３９】
図３は、サージサイクルの発生前に、コンプレッサの数回転で各センサによって与えられる正規化された圧力信号の一例を示している。
【００４０】
この図では、サージの発生前に８０回の回転で現われる圧力のピークが観察される。
【００４１】
あらゆるセンサによって与えられる信号において見られるこれらピークは、コンプレッサの回転速度の一部分において伝播する局所的擾乱である。
【００４２】
サージに入る前に発生する局所的擾乱は、短時間で複数の羽根に伝播する部分的剥離に対応する。
【００４３】
したがって、これら擾乱を体系的に検知することができる方法が決定される。
【００４４】
図４には、（ａ）回転剥離セルの伝播によって影響を受ける羽根の翼列、（ｂ）回転剥離セルの発生時にセンサによって観察される圧力信号Ｐが示されている。
【００４５】
羽根間チャネルの中に剥離セルが出現すると、センサの信号に対する剥離されたチャネルの羽根の寄与が、大きく変更される。実際、剥離セルは、羽根間チャネルの局所的及び瞬間的な障害となって現われる。
【００４６】
圧力信号の変更は、同一センサの前の連続する２つの羽根の通過時間Ｔａより長い時間Ｔｚ中に、ゼロの値を通らないという点である。
【００４７】
ゼロを通過しないことは、正規化された圧力信号上の剥離セルの特徴を構成し、コンプレッサ内の航空力学的不安定性の前兆の属性として使用される。
【００４８】
したがって、図５から７に示されているように、本発明による方法は、各センサについて、獲得され（２０）、正規化された（３０）圧力信号に基づき、ゼロを通らないというこの特性が確認されたあらゆる事象を読取ることからなる。
【００４９】
角度位置θで時間ｔにおいて、コンプレッサの円周上に置かれたセンサによる事象検知４０の結果は、関数Ｅ（ｔ，θ）によって表され、ここにおいてＥ（ｔ，θ）は、ゼロを通らないという特性が確認された場合には、論理値１をとり、そうでない場合には値０をとる。
【００５０】
ピークの伝播速度は、コンプレッサの羽根の通過速度に対応する周波数と比べて低い周波数に対応するにもかかわらず、信号のゼロを通らない特性を評価するために、圧力信号のサンプリング周波数は、これらのピークの体系的検知を可能にするために、十分高くなければならない。すなわち、羽根の通過周波数の少なくとも２．５倍を越えなければならない。
【００５１】
次に本発明は、航空力学的不安定性が目前に迫っているかどうかを決定する（６０）ことができるように、検知された事象全体を分析する（５０）ことからなる。
【００５２】
種々のタイプの分析を行うことができる。
【００５３】
たとえば、図５に示されているように、事象の分析は、センサの前を連続する２つの羽根が通過する時間Ｔａより長い、あらかじめ設定された特定の時間中に、センサによって検知された事象を数え上げ（５１）、事象の検知頻度を決定する（５２）ことからなることもできる。こうして得られた検知頻度の値は、不安定性が目前に迫っていかどうかの少なくとも第一の指示を得るために、あらかじめ設定された第一の閾値と比較される（５３）。実際に、剥離セルの出現数が多くなればなるほど、航空力学的不安定性に到達する確率も大きくなる。こうして、決定（６０）は、航空力学的流体流出の安定性について下されることもできる。
【００５４】
この分析の利点は、その分析によって、ただ一つのセンサによって与えられるただ一つの圧力信号に基づき、流体流出の安定性について、第一の迅速な評価を得ることができる点である。しかしながら、この第一の評価は、センサ数がより多くなれば、その信頼度がより大きくなる。
【００５５】
より信頼度の高い評価を得ることができる別のタイプの分析は、複数のセンサが、コンプレッサの円周上に均等に配分されるように配置される場合に行うことができる。このタイプの分析の種々の段階が図６に示され、コンプレッサの円周上に均等に配分された４つのセンサによる実施形態が図８に示されている。この分析は、検知された事象の空間−時間的位置を読取る（５５）ことからなり、次にコンプレッサのロータの回転速度と０との間に含まれる速度の走査を行うことによって、さらに、検討された各速度値について、この速度に対応する伝播方向に沿って並べられる検知された事象を探索することによって、検知された各事象に対応する擾乱の伝播方向の探索（５６）が、すべてのセンサから与えられる信号に基づいて行なわれる。図８における例では、１０個の事象が、４つのセンサの角位置θ１、θ２、θ３、θ４について検知され、速度ｖ１、ｖ２、ｖ３に対応する３つの伝播方向が選択された。
【００５６】
このとき、コンプレッサのロータの回転速度の一部分における局所的擾乱の空間的伝播の存在に関する蓋然性（ｖｒａｉｓｅｍｂｌａｎｃｅ）のレベルＪは、求められた各伝播方向に割当てられる（５７）。時間及び空間において位置決定された各事象について、蓋然性のレベルＪは、時間ｔ、空間位置θ、そして伝播速度ｖによって異なる関数である。
【００５７】
時間ｔ、位置θ、速度ｖにおいて検知された事象についての、蓋然性のレベルＪ０（ｔ，θ，ｖ）の決定のために、種々の基準を利用することができる。しかしながら、それらの基準はすべて、速度ｖに対応する伝播方向における当該事象より先に検知された事象数の関数である。
【００５８】
速度ｖにおいて先に検知された事象数がゼロである場合には、すなわち、あらゆる前の値Ｅｋ（ここでｋは正の整数である）が０である場合、蓋然性のレベルはゼロにならなければならない。実際に、この方向においていかなる事象も検知されない場合には、速度ｖに対応する方向では、擾乱の伝播が存在するという可能性はない。
【００５９】
逆に、あらゆる前の値Ｅｋが１である場合には、蓋然性のレベルは最大でなければならない。図８では、検知された１０個の事象のうちの、２つの事象が伝播方向ｖ１にしたがって並べられ、また２つの事象が伝播方向ｖ２にしたがって並べられるの対して、７つの事象は、伝播方向ｖ３にしたがって並んでいることが分かる。したがって、速度ｖ１及びｖ２における伝播の蓋然性のレベルＪは、信頼度が低くなければならず、不安定性が存在しないことの決定を導き出すのに対して、速度ｖ３における伝播の蓋然性のレベルは大きく、不安定性の存在の決定を導き出さなければならない。
【００６０】
前の値は、先験的に無限の数であるので、前の値の有限の数だけを関与させる近似を行うことが必要になる。
【００６１】
近似は、たとえばｍ＝７である図９に示されているように、前の値のあらかじめ設定された一定数ｍしか考慮に入れないことによって、丸めるという手段によって行うことができる。この場合、蓋然性のレベルＪは、考慮された前の値Ｅｋの関数である。この場合、蓋然性のレベルＪは、検知された最後のｍ個の事象に分析を限定することによって、前の事象Ｅｋの値の和によって加重された当該の事象Ｅ０の値に等しくなるように選択することができる。
【００６２】
【数１】
0004155715.tif 000002
係数Ｅ０（ｔ，θ，ｖ）は、（ｔ，θ）においていかなる事象も検知されない場合には、Ｊに値ゼロを割当てることができる。図９に示されている例においては、Ｊ＝６である。
【００６３】
この近似は、実施が容易で迅速であるという利点を持っているが、前の数のうちの限定された数ｍしか組み入れず、検知された事象においていかなる連続性の区別も行わない。
【００６４】
前の値Ｅｋに対する無限の従属性を組み込む帰納的方法を利用し、しかも、たとえば図１０に示されているような蓋然性のレベルの値の有限数ｍのみを考慮に入れることによって、別のタイプの近似を行うことができる。この例においては、蓋然性のレベルＪ（ｔ，θ，ｖ）は、速度ｖに対応する伝播方向において、伝播する事象について計算される前の値Ｊｋの関数である。
【００６５】
【数２】
0004155715.tif 000003
図１０に示されている例においては、検知された最後の７つの事象に限定することによって（ｍ＝７）、
Ｊ７＝１、Ｊ６＝２、Ｊ５＝０、Ｊ４＝４、Ｊ３＝０、Ｊ２＝８、Ｊ１＝１６、Ｊ０（ｔ，θ，ｖ）＝３２となる。
【００６６】
この近似の利点は、検知された事象の前の値の無限の寄与を考慮に入れることができるという点であるが、この近似は、図９を参照して実施される方法よりもコストがかかる。なぜなら、各伝播速度について各反復において計算される蓋然性のレベルの過去の値を保存することが必要となるからである。さらに図９の例におけるように、この近似は、検知された事象の連続性についての区別を行うことができない。
【００６７】
蓋然性のレベルの決定の中に、検知された事象の連続性を考慮に入れるためには、単純な和の異なる関数を選択することが必要である。これらの関数は、連続する複数のセンサによる方向ｖにおける事象の伝播方向が、同じ数のセンサによって検知された断続的伝播より決定的となるように、選択されなければならない。
【００６８】
たとえば、λｋと記されるネグレクション係数（ｆａｃｔｅｕｒ  ｄ’ｏｕｂｌｉ）によって、帰納的方法の中で使用される前の値Ｊｋの和を加重することからなる関数を選択することができる。ここでλｋは、厳密に１より小さな正の数である。たとえば、
【００６９】
【数３】
0004155715.tif 000004
この数式において、ｍは事象の数であり、分析はこの数に限定される。
【００７０】
この加重は、断続的な伝播を、連続する伝播よりも小さくすることができる。
【００７１】
蓋然性のレベルを評価するために選択された方法に関係なく、次に、蓋然性のレベルは、コンプレッサの航空力学的安定性における第二の評価を得ることができるように、あらかじめ設定された第二の閾値と比較される（５８）。閾値は、上述の方法の各々について異なる。たとえば、閾値は、コンプレッサの円周の半分以上において連続的に検知される局所的擾乱の伝播が、サージの早期検知信号を始動させるように規定されることもできる。
【００７２】
一定の瞬間にＪの最大値が、選択された閾値に等しいまたはそれより大きくなる場合には、伝播の存在のための決定が下される（６０）。
【００７３】
図７に示されているように、上述の２つのタイプの分析はまた、決定を下す前に、コンプレッサの航空力学的安定性についての２つの指示を得るために、同時に行うことができる。
【図面の簡単な説明】
【図１】本発明による非定常圧力センサを備えたタービンエンジンのコンプレッサの一例を示す図である。
【図２】ロータの羽根の翼の前縁の上流側に置かれたセンサによって見られる圧力信号の一例を示す図である。
【図３】サージサイクルが発生する前に、コンプレッサの数回転数において種々のセンサによって与えられる正規化された圧力信号の一例を示す図である。
【図４ａ】回転剥離セルの伝播によって影響を与えられる羽根の翼列の一例を示す図である。
【図４ｂ】回転剥離セルの発生時にセンサによって見られる圧力信号の一例を示す図である。
【図５】本発明の種々の実施形態に応じた方法の主な段階を示すフローチャートである。
【図６】本発明の種々の実施形態に応じた方法の主な段階を示すフローチャートである。
【図７】本発明の種々の実施形態に応じた方法の主な段階を示すフローチャートである。
【図８】本発明による方法の適用例を示す図である。
【図９】本発明による方法の適用例を示す図である。
【図１０】本発明による方法の適用例を示す図である。
【符号の説明】
１０  センサ
１１  羽根の翼列
２０  信号の獲得
３０  信号の正規化
４０  事象の検知
４５  擾乱の伝播方向の探索
５０  検知された事象の分析
５１  事象の列挙
５２  検知頻度の決定
５３  第一の閾値との比較
５５  事象の空間−時間的位置読取り
５６  伝播方向の探索
５７  蓋然性のレベルの計算
５８  第二の閾値との比較
６０  決定
【図１】
0004155715.tif 000005
【図２】
0004155715.tif 000006
【図３】
0004155715.tif 000007
【図４ａ】
0004155715.tif 000008
【図４ｂ】
0004155715.tif 000009
【図５】
0004155715.tif 000010
【図６】
0004155715.tif 000011
【図７】
0004155715.tif 000012
【図８】
0004155715.tif 000013
【図９】
0004155715.tif 000014
【図１０】
0004155715.tif 000015