(57)【特許請求の範囲】
【請求項１】
圧縮機（１）を作動させる方法（１００）であって、
−  前記圧縮機のそれぞれの吸込み部分または吐出し部分におけるそれぞれの複数の測定値から得られた複数の測定データ（Ｐs、Ｐd、Ｔs、Ｔd、ｈs、ｈd）を取得するステップ（１０５）と、
−  前記圧縮機（１）によって圧縮されたガスの分子量（Ｍw）を計算することにより、前記測定データ（Ｐs、Ｐd、Ｔs、Ｔd、ｈs、ｈd）の相合を確認するステップ（１２０）と、
−  前記測定データ（Ｐs、Ｐd、Ｔs、Ｔd、ｈs、ｈd）の第１の測定値に障害がある場合には、前記分子量（Ｍw）の最近の有効な値および前記測定データ（Ｐs、Ｐd、Ｔs、Ｔd、ｈs、ｈd）の有効な測定値に基づいて前記第１の測定値を推定値で置換するステップ（１３０）と、
−  前記推定値および前記測定データ（Ｐs、Ｐd、Ｔs、Ｔd、ｈs、ｈd）の有効な測定値に基づいて、サージ防止マップ（３００、４００）上の推定の動作点（３０２、４０２）を求めるステップと、
を含む、方法（１００）。
【請求項２】
前記置換するステップ（１３０）が、所定の安全時間期間（ｔ1）中に遂行される請求項１に記載の方法（１００）。
【請求項３】
前記測定データ（Ｐs、Ｐd、Ｔs、Ｔd、ｈs、ｈd）の第２の測定値に障害がある場合、または前記安全時間期間（ｔ1）の最後に、
−  前記第１および第２の測定値の最大値および／または最小値に基づいて、前記第１および第２の測定値をそれぞれの最悪の場合の値で置換するさらなるステップ（１４０）と、
−  前記最悪の場合の値および前記測定データ（Ｐs、Ｐd、Ｔs、Ｔd、ｈs、ｈd）の前記有効な測定値に基づいて、前記サージ防止マップ（３００、４００）上の最悪の場合のポイント（３０３、４０３）を求めるステップと、
をさらに含む、請求項２に記載の方法（１００）。
【請求項４】
前記測定データ（Ｐs、Ｐd、Ｔs、Ｔd、ｈs、ｈd）の前記相合を確認する前記ステップ（１２０）において、前記計算された分子量（Ｍw）が許容値の区間と比較される、請求項１乃至３のいずれかに記載の方法（１００）。
【請求項５】
前記サージ防止マップ（３００、４００）が無次元のサージ防止マップである、請求項１乃至４のいずれかに記載の方法（１００）。
【請求項６】
前記第１の測定値が、前記サージ防止マップのタイプと、前記圧縮機の流量エレメント（ＦＥ）の位置とに依拠するものである、請求項１乃至５のいずれかに記載の方法（１００）。
【請求項７】
前記第１の測定値が、
−  吸込み（Ｐs）における圧力、
−  吐出し（Ｐd）における圧力、
−  吸込みの流量エレメントにおける圧力降下（ｈs）または吐出しの流量エレメントにおける圧力降下（ｈd）、
−  吸込みの温度（Ｔs）、
−  吐出しの温度（Ｔd）、
のうちの１つである、請求項１乃至６のいずれかに記載の方法（１００）。
【請求項８】
１つまたは複数のデジタルコンピュータ上で実行されたとき、請求項１乃至７のいずれか１項記載の方法のステップを実行するのに適切であるソフトウェアコードの部分を含む、コンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、１つまたは複数の測定信号に障害があるとき、サージ防止コントローラに、サージ防止弁を開くことによる介入をさせることなく、圧縮機を作動させ続けると同時に、複数のフォールバック方策によって適切なレベルの保護をもたらす方法に関するものである。
【背景技術】
【０００２】
サージ防止コントローラは、不変の圧縮機マップにおける圧縮機の動作点の位置を識別するために、複数のセンサおよび伝送器による複数のフィールド測定値を必要とする。たとえば伝送器とコントローラの間の通信損失といった、必要な測定に障害がある場合には、動作点の位置は評価されない。このとき、圧縮機を安全に作動させるために、一般に、最悪の場合の手法が用いられる。この手法では、障害のある測定値が、動作点をサージラインの方へ、できる限り安全にシフトすることを可能にする値で置換される。たとえば吸込みに流量エレメントを含んでいる圧縮機の設備では、
−  吐出し圧力値の損失の場合には、後者は、その最大の可能な値で置換され、
−  流量エレメントの差圧値（ｈ）の損失の場合には、そのような差圧の最低の可能な値（すなわちゼロ値）が選択される。
【０００３】
いずれにせよ、この最悪の場合の手法は、実際の動作条件からすれば必要でないときさえサージ防止弁を開く傾向があり、通常、処理の有用性が失われる。
【０００４】
したがって、圧縮機を安全に作動させるのと同時に、公知の従来技術の上記の不都合を防止することを可能にする改良された方法を提供することが望ましいであろう。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】米国特許第４８２５３８０号明細書
【発明の概要】
【０００６】
第１の実施形態によれば、本発明は、圧縮機を作動させる方法を提供することによってそのような目的を達成するものであり、この方法は、
−  圧縮機のそれぞれの吸込み部分または吐出し部分におけるそれぞれの複数の測定値から得られた複数の測定データを取得するステップと、
−  圧縮機によって圧縮されたガスの分子量を計算することにより、測定データの相合を確認するステップと、
−  前記測定データの第１の測定値に障害がある場合には、前記分子量の最近の有効な値および前記測定データの有効な測定値に基づいて前記第１の測定値を推定値で置換するステップと、
−  前記推定値および前記測定データの有効な測定値に基づいて、サージ防止マップ上の推定の動作点を求めるステップとを含む。
【０００７】
本発明の別の態様によれば、前記第１の測定値を推定値で置換する前記ステップは、所定の安全時間期間中に遂行される。
【０００８】
本発明のさらなる態様によれば、この方法は、前記測定データの第２の測定値に障害がある場合、または安全時間期間の最後に、
−  前記第１および第２の測定値の最大値および／または最小値に基づいて、前記第１および第２の測定値をそれぞれの最悪の場合の値で置換するさらなるステップと、
−  前記最悪の場合の値および前記測定データの有効な測定値に基づいて、サージ防止マップ上の最悪の場合のポイントを求めるステップとを含む。
【０００９】
そのような方法では、無次元の解析によって与えられた圧縮機の挙動モデルを考慮して、１つの障害のある測定値が、残りの複数の健全な測定データを用いることによって計算される。マップ上で、測定された動作点を推定の動作点で置換すると、動作点の位置決めに関する不連続が防止され、したがって、サージ防止制御の不要な介在および処理の混乱が防止される。
【００１０】
本発明の他の目的の特徴および利益が、本発明の実施形態の以下の説明を以下の図面とともに理解することから明らかになるであろう。
【図面の簡単な説明】
【００１１】
【図１】本発明による、圧縮機を作動させる方法の全般的なブロック図である。
【図２】図１の方法の部分的ブロック図である。
【図３ａ】本発明の方法によって作動され得る圧縮機の第１の例の概略図である。
【図３ｂ】図３ａの圧縮機のサージ防止マップの図である。
【図４】図３ａの圧縮機に関して、図１の方法によって管理され得る障害状態に対応する図３ｂのサージ防止マップの図である。
【図５】図３ａの圧縮機に関して、図１の方法によって管理され得る、図４のものとは別の障害状態に対応する図３ｂのサージ防止マップの図である。
【図６】図３ａの圧縮機に関して、図１の方法によって管理され得る、図４、図５のものとは別の障害状態に対応する図３ｂのサージ防止マップの図である。
【図７ａ】本発明の方法によって作動され得る圧縮機の第２の例の概略図である。
【図７ｂ】図７ａの圧縮機のサージ防止マップの図である。
【図８】図７ａの圧縮機に関して、図１の方法によって管理され得る、障害状態に対応する図７ｂのサージ防止マップの図である。
【図９】図７ａの圧縮機に関して、図１の方法によって管理され得る、図８のものとは別の障害状態に対応する図７ｂのサージ防止マップの図である。
【図１０】図７ａの圧縮機に関して、図１の方法によって管理され得る、図８、図９のものとは別の障害状態に対応する図７ｂのサージ防止マップの図である。
【図１１】図７ａの圧縮機に関して、図１の方法によって管理され得る、図８〜図１０のものとは別の障害状態に対応する図７ｂのサージ防止マップの図である。
【図１２】図７ａの圧縮機に関して、図１の方法によって管理され得る、図８〜図１１のものとは別の障害状態に対応する図７ｂのサージ防止マップの図である。
【発明を実施するための形態】
【００１２】
図１と、図３ａおよび図７ａの例の概略図とを参照して、本発明による、遠心圧縮機１を作動させる方法が、全体として１００で示されている。方法１００は、サージ防止マップ上の動作点を求めるのに用いられる測定値の妥当性を検証することによって圧縮機１を作動させるものである。１つまたは複数の測定値が欠落している場合、フォールバック方策がもたらされる。方法１００の最後に、測定値または計算値である複数の値が、サージ防止マップ上の動作点を計算するのに利用可能になる。
【００１３】
この方法は、たとえばＰＬＣシステムといった、圧縮機１に関連した制御ユニットによって繰り返し実行される。方法１００の２つの連続した実行の間の時間間隔は、一般に制御ユニット（ＰＬＣ）の走査時間に対応するものである。
【００１４】
方法１００は、遠心圧縮機１の吸込みおよび吐出しに接続されたそれぞれの複数の測定器から複数の測定データを取得する予備のステップ１０５を含む。測定データは、
−  吸込み圧力Ｐsと、
−  吐出し圧力Ｐdと、
−  吸込み温度Ｔsと、
−  吐出し温度Ｔdと、
−  それぞれ吸込みまたは吐出しにおける流量エレメントＦＥ上の差圧ｈs＝ｄＰsまたはｈd＝ｄＰdとを含む。
【００１５】
上記のデータは、サージ防止マップ上の圧縮機１の動作点を求めるために通常用いられるものである。
【００１６】
方法１００に用いられるサージ防止マップは、無次元のサージ防止マップである。様々なタイプのサージ防止マップが用いられ得る。流量エレメントＦＥが圧縮機１の吸込み側に配置されている場合、ｈs／Ｐs（横座標）対Ｐd／Ｐs（縦座標）のマップ３００が用いられる（図３ｂ、図４〜図６）。無次元のマップ３００が用いられるとき、マップ上の動作点位置を識別するのに、ｈs、ＰsおよびＰdの３つの測定値が必要とされる。完全な無次元の解析には、以下でより詳細に説明されるように、吸込みのガス温度Ｔsおよび吐出しのガス温度Ｔdの測定値が必要である。流量エレメントＦＥが圧縮機１の吐出し側に配置されている場合、ｈs／Ｐs対Ｐd／Ｐsマップ４００が用いられる（図７ｂ、図８〜図１０）。しかしながら、後者の場合、ｈs＝ｄＰsは有効でなく、当技術分野で公知の次式を用いて計算する必要がある。
【００１７】
ｈs＝ｈd・（Ｐd／Ｐs）・（Ｔs／Ｔd）・（Ｚs／Ｚd）  （Ａ）
マップ４００上の動作点の位置を識別するために式Ａを適用するには、ｈd、Ｐs、Ｐd、Ｔs、Ｔdの５つの測定値の組が必要である。
【００１８】
あるいは、流量エレメントＦＥが吸込みに配置されていようと吐出しに配置されていようと、圧縮比Ｐd／Ｐsの代わりに、低減されたヘッドｈrが、横軸上のｈs／Ｐsとともに縦軸上にマッピングされ得る。後者のマップが用いられるとき、マップ上の動作点の位置を識別するには、ｈrを計算することにより、ｈs、Ｐs、Ｐd、Ｔs、Ｔdの５つの測定値が必要である。
【００１９】
方法１００は、予備のステップ１０５の後、圧縮機１の吸込みおよび吐出しに接続されている複数の測定器の間で測定器の障害を検出する第１の有効なステップ１１０を含む。
【００２０】
第１のステップ１１０の間に測定器の障害が検出されなければ、方法１００は、複数の測定データの相合を確認する第２の有効なステップ１２０を始める。第２のステップ１２０は、流量エレメントｈsまたはｈdの差圧における温度の測定データＴs、Ｔd、圧力の測定データＰs、Ｐdと、以下で、吸込み状態における分子量とガス圧縮性Ｚsの間の比Ｍw／Ｚsの計算に関して説明される（図２に表される）プロシージャ２００とに基づいて、圧縮機１によって圧縮されたガスの分子量Ｍwを計算する第１のサブステップ１２１を含む。
【００２１】
プロシージャ２００は、プロシージャ２００の先の実行で計算された値を用いて比Ｍw／Ｚsの第１の値を設定する初期化動作２０１を含む。プロシージャ２００が初めて実行されるので、そのような値が利用可能でない場合には、吸込み状態における分子量Ｍwおよびガス圧縮性Ｚsの設計条件値が用いられる。反復プロシージャ２００には、初期化動作２０１の後に、後続の操作２１１〜２２０を連続的に遂行するサイクル２１０が含まれる。
【００２２】
反復サイクル２１０の１回目の動作２１１中に、当技術分野で公知の次式によって吸込み密度γsが計算され、
γs＝Ｐs／（Ｒ・Ｔs）・（Ｍw／Ｚs）i-1  （Ｂ）
（Ｍw／Ｚs）i-1は、反復サイクル２１０の以前の反復で計算された、または反復サイクル２１０が初めて実行される以前に初期化動作２０１で計算されたＭw／Ｚsの値である。
【００２３】
反復サイクル２１０の２回目の動作２１２中に、当技術分野で公知の次式によって体積流量Ｑvsが計算され、
Ｑvs＝ｋFE  ｓｑｒｔ（ｈs・１００／γs）  （Ｃ）
ｋFEは流量エレメントＦＥの定数であり、「ｓｑｒｔ」は平方根関数である。流量エレメントＦＥが圧縮機１の吐出し側に配置されていてマップ４００が用いられる場合、ｈsは直接測定されず、式Ａを用いて計算され得る。
【００２４】
反復サイクル２１０の３回目の動作２１３中に、当技術分野で公知の次式によってインペラの先端の速度ｕ1が計算され、
ｕ1＝Ｎ・Ｄ・π／６０  （Ｄ）
Ｎはインペラの回転速度であり、Ｄはインペラの直径である。
【００２５】
反復サイクル２１０の４回目の動作２１４中に、当技術分野で公知の次式によって流量の無次元の係数φ1が計算される。
【００２６】
φ1＝４・Ｑvs／（π・Ｄ2・ｕ1）  （Ｅ）
反復サイクル２１０の５回目の動作２１５中に、当技術分野で公知の次式によって吸込みにおける音速ａsが計算され、
ａs＝ｓｑｒｔ（ｋv・ＲＴs／（Ｍw／Ｚs）i-1）  （Ｆ）
ｋvは断熱指数である。
【００２７】
反復サイクル２１０の６回目の動作２１６中に、吸込みにおけるマッハ数Ｍ1が、インペラの先端の速度ｕ1と吸込みにおける音速ａsの間の比として計算される。
【００２８】
反復サイクル２１０の７回目の動作２１７中に、φ1として既知の無次元のデータ配列からの補間およびマッハ数Ｍ1によって、ヘッドの無次元の係数τとポリトロープ効率ｅｔａｐの間の積が導出される。
【００２９】
反復サイクル２１０の８回目の動作２１８中に、当技術分野で公知の次式によってポリトロープヘッドＨpcが計算される。
【００３０】
Ｈpc＝τ・ｅｔａｐ・ｕ12  （Ｇ）
反復サイクル２１０の９回目の動作２１９中に、当技術分野で公知の次式によってポリトロープ指数ｘが計算される。
【００３１】
ｘ＝ｌｎ（Ｔd／Ｔs）／ｌｎ（Ｐd／Ｐs）  （Ｈ）
反復サイクル２１０の１０回目の最後の動作２２０中に、当技術分野で公知の次式によって比Ｍw／Ｚsの値が更新される。
【００３２】
（Ｍw／Ｚs）i＝ＲＴs・（（Ｐd／Ｐs）x−１）／（Ｈpc・ｘ）  （Ｉ）
第２のステップ１２０の第２のサブステップ１２２では、Ｍw／Ｚsの計算された値が、最小値と最大値の間で定義された許容値の区間と比較される。Ｍw／Ｚsの計算された値がそのような区間の外にあると、第２のステップ１２０の後続の第３のサブステップ１２３で警報が発生される。第２のサブステップ１２２中に遂行される比較検査により、遠心圧縮機１の吸込みおよび吐出しにおける複数の測定器によって遂行された複数の測定値Ｐs、Ｐd、Ｔs、Ｔd、ｈsまたはｈdを確認することが可能になる。これは、具体的には、オペレータが、較正されていない測定器を起動中に識別するのを支援するのに用いることができる。
【００３３】
第１の有効なステップ１１０の間に測定器の障害が検出されると、方法１００は、複数の測定器が障害状態にあるかどうか検出する第３のステップ１１３を始める。第３のステップ１１３中に遂行された検査がマイナスであって、測定器の障害が１つしか検出されなかった場合には、方法１００は、所定の安全時間期間ｔ1にわたって、分子量の最近の有効な値および他の有効な測定データの値に基づいて欠落データ（Ｐs、Ｐd、Ｔs、Ｔd、ｈsまたはｈdのうちの１つ）を推定値で置換するフォールバックステップ１３０を継続する。
【００３４】
方法１００は、フォールバックステップ１３０に入る以前に、安全時間期間ｔ1が第４のステップ１１４および第５のステップ１１５を含むかどうか識別するために、それぞれ、フォールバックステップ１３０が進行中であるか、また安全時間期間ｔ1が経過したか、検査する。第４のステップ１１４の間および第５のステップ１１５の間に遂行された検査のうちの１つがマイナスであって、フォールバックステップ１３０がまだ進行中でないか、または安全時間期間ｔ1がまだ経過していなければ、フォールバックステップ１３０が遂行される。
【００３５】
第４のステップ１１４中に遂行された検査がマイナスであれば、方法１００は、フォールバックステップ１３０の第１のサブステップ１３１を継続し、タイマが、安全時間期間ｔ1を測定し始める。第４のステップ１１４中に遂行された検査がプラスであって、フォールバックステップ１３０が既に進行中であれば、第５のステップ１１５が遂行される。第５のステップ１１５中に遂行されたマイナスの検査および第１のサブステップ１３１の後、すなわちフォールバックステップ１３０が進行中であって、安全時間期間ｔ1がまだ終了していなければ、方法１００は、フォールバックステップ１３０の第２のサブステップ１３２を継続し、欠落データの推定値が求められる。フォールバックステップ１３０は、第２のサブステップ１３２の後に、特に圧縮機１のオペレータに対して、測定器のうちの１つが障害状態にあり、関連したフォールバックステップ１３０が遂行されていることを伝えるために警報を発生する第３のサブステップ１３３を含む。
【００３６】
フォールバックステップ１３０の第２のサブステップ１３２中に遂行される動作は、測定器のどれが障害状態にあり、したがって、どの測定データが欠けているかということに依拠する。すべての場合に、フォールバックステップ１３０の第２のサブステップ１３２中に、Ｍw／Ｚsの最近の有効な適値、すなわち、測定器障害が生じる直前の第２のステップ１２０の第１のサブステップ１２１で計算された値が用いられる。
【００３７】
すべての場合に、任意選択で、安全性をさらに改善するために、フォールバックステップ１３０の第２のサブステップ１３２中に、サージ防止マップ３００、４００においてサージ防止マージンが増加される。
【００３８】
本発明の第１の実施形態（図３ａ、図３ｂ、図４〜図６）では、圧縮機１は、吸込み側に流量エレメントＦＥを含み、横座標の変数としてｈs／Ｐsがマッピングされ、縦座標の変数としてＰd／Ｐsがマッピングされている無次元のマップ３００が用いられる。正常状態では、マップ３００上の測定された動作点３０１を求めるには、流量エレメントＦＥからの差圧ｈsと、吸込みの圧力センサからのＰsおよび吐出しの圧力センサからのＰdとの測定値で十分である。障害状態では、ｈs、ＰsまたはＰdの測定値のうちの１つが欠けていると、測定された動作点３０１の割出しが妨げられ、フォールバック推定を遂行する必要がある。フォールバック推定の間、以下で明白になるように、吸込みにおける温度の値Ｔsおよび吐出しにおける温度の値Ｔdが必要とされる。
【００３９】
本発明の第１の実施形態では、障害状態の測定器が流量エレメントＦＥであれば、差圧ｈsは、フォールバックステップ１３０の第２のサブステップ１３２において、連続的に遂行される以下の動作によって推定される。
−  ポリトロープ指数ｘが式Ｈを用いて計算され、
−  ポリトロープヘッドＨpcが、最近の有効な適値Ｍw／Ｚsならびに既知のＴs、Ｐd／Ｐsおよびｘを用いて式Ｉから計算され、
−  Ｈpcが既知であり、ｕ1が式Ｄを用いて計算されるので、ポリトロープヘッドの無次元の係数τとポリトロープ効率ｅｔａｐの間の積が式Ｇから計算され、
−  音速ａsが、式Ｆおよび最近の有効な適値Ｍw／Ｚsを用いて計算され、
−  マッハ数Ｍ1がｕ1とａsの間の比として計算され、
−  積τ・ｅｔａｐが既知であるため、流量の無次元の係数φ1が、サイクル２１０の７回目の動作２１７で用いられるのと同一の無次元のデータ配列からの補間によって導出され、
−  体積流量Ｑvsが式Ｅから計算され、
−  吸込み密度γsが式Ｂによって計算され、
−  Ｑvs、ｋおよびγsが既知であるため、差圧ｈsが式Ｃから計算される。
【００４０】
図４を参照すると、測定値ＰsおよびＰdならびにｈsの推定に基づいて、マップ３００で、測定された動作点３０１が推定の動作点３０２によって置換される。ｈsを求めるのに用いられる計算および補間における誤差のマージンを考慮して、推定の動作点３０２は、測定された動作点３０１を含む円形の領域に入る。普通には、そのような領域は、ＳＬＬの右側の安全領域にあるか、または最悪シナリオの手法で計算された動作点よりも安全領域に少なくともより接近しているはずである。公知の方法で用いられる最悪のシナリオでは、測定された動作点３０１は、ｈs＝０という仮定に基づいて、マップ３００においてマップ３００の縦軸上で、最悪の場合のポイント３０３で置換される。したがって、最悪の場合のポイント３０３は、常にＳＬＬの左側にあり、サージ防止弁を完全に開けさせる。
【００４１】
本発明の第１の実施形態では、障害状態の下で測定器が吸込みの圧力センサである場合、吸込圧力Ｐsは、フォールバックステップ１３０の第２のサブステップ１３２において、繰り返し遂行される以下の動作によって推定される。
−  最初に、Ｐsが、障害状態に到達する前に吸込圧力センサによって測定された最近の有効な適値として定義され、
−  吸込み密度γsが、最近の有効な適値ＰsおよびＭw／Ｚsと既知のＴsとを用いて式Ｂによって計算され、
−  体積流量Ｑvsが式Ｃによって計算され、
−  流量の無次元の係数φ1が式Ｅによって計算され、
−  音速ａsが式Ｆを用いて計算され、
−  マッハ数Ｍ1がｕ1とａsの間の比として計算され、
−  ヘッドの無次元の係数τとポリトロープ効率ｅｔａｐの間の積が、マッハ数Ｍ1および上記で計算されたφ1の値を用いて、無次元のデータ配列からの補間によって導出され、
−  ポリトロープヘッドＨpcが式Ｉによって計算され、
−  ポリトロープ指数ｘが、Ｍw／Ｚsの最近の有効な適値を用いて、当技術分野で公知の次式で計算され、
ｘ＝Ｒ・（Ｔd−Ｔs）／（Ｍw／Ｚs）／Ｈpc  （Ｌ）
−  最後に、ｘ、Ｐd、ＴsおよびＴdが既知であるため、Ｐsの新規の値が式Ｈから計算される。
【００４２】
図５を参照すると、測定値ｈsおよびＰdならびにＰsの推定に基づいて、マップ３００で、測定された動作点３０１が推定の動作点３０２によって置換される。Ｐsを求めるのに用いられる計算および補間における誤差のマージンを考慮して、推定の動作点３０２は、測定された動作点３０１を含む円形の領域に入る。普通には、そのような領域は、ＳＬＬの右側の安全領域にあるか、または最悪シナリオの手法で計算された動作点よりも安全領域に少なくともより接近しているはずである。公知の方法で用いられる最悪のシナリオでは、測定された動作点３０１は、Ｐd／Ｐs＝Ｐd／Ｐs,minおよびｈs／Ｐs＝ｈs／Ｐs,maxという仮定に基づいて、マップ３００において、最悪の場合のポイント３０３で置換され、Ｐs,minおよびＰs,maxは、それぞれ、吸込みにおける圧力の可能な最小値および最大値である。一般に、最悪の場合のポイント３０３は、この場合もＳＬＬの左側にあり、サージ防止弁を開けさせる。
【００４３】
本発明の第１の実施形態では、障害状態の測定器が吐出しの圧力センサであれば、吐出し圧力Ｐdは、フォールバックステップ１３０の第２のサブステップ１３２において、以下の動作によって推定される。
−  吸込み密度γsが式Ｂによって計算され、
−  体積流量Ｑvsが式Ｃによって計算され、
−  流量の無次元の係数φ1が式Ｅによって計算され、
−  音速ａsが、最近の有効な適値Ｍw／Ｚsを用いて式Ｆによって計算され、
−  マッハ数Ｍ1がｕ1とａsの間の比として計算され、
−  ヘッドの無次元の係数τとポリトロープ効率ｅｔａｐの間の積が、マッハ数Ｍ1および上記で計算されたφ1の値を用いて、無次元のデータ配列からの補間によって導出され、
−  ポリトロープヘッドＨpcが式Ｇから計算され、
−  ポリトロープ指数ｘが、最近の有効な適値Ｍw／Ｚsを用いて式Ｌによって計算され、
−  ｘ、Ｐs、ＴsおよびＴdが既知であるため、Ｐdが式Ｈから計算される。
【００４４】
図６を参照すると、測定値ｈsおよびＰsならびにＰdの推定に基づいて、マップ３００で、測定された動作点３０１が推定の動作点３０２によって置換される。変数としてマップ３００の縦軸上にのみ現われるＰdを求めるのに用いられる計算および補間における誤差のマージンを考慮して、推定の動作点３０２は、測定された動作点３０１を含む細長い垂直な領域に入る。普通には、そのような領域は、ＳＬＬの右側の安全領域にあるか、または最悪シナリオの手法で計算された動作点よりも安全領域に少なくともより接近しているはずである。公知の方法で用いられる最悪のシナリオでは、測定された動作点３０１は、Ｐd／Ｐs＝Ｐd,max／Ｐsという仮定に基づいて、マップ３００において、最悪の場合のポイント３０３で置換され、Ｐd,maxは、吐出しにおける圧力の最大の可能な値である。一般に、最悪の場合のポイント３０３は、この場合もＳＬＬの左側にあり、サージ防止弁を開けさせる。
【００４５】
本発明の第２の実施形態（図７ａ、図７ｂ、図８〜図１２）では、圧縮機１は吐出し側の流量エレメントＦＥを含み、横座標の変数としてｈs／Ｐsがマッピングされ、縦座標の変数としてＰd／Ｐsがマッピングされている無次元のマップ４００が用いられる。差圧ｈsを測定値から入手することができないので、関連した値は式Ａによって計算される。正常状態では、マップ４００上の測定された動作点４０１を求めるには、流量エレメントＦＥからの差圧ｈdと、吸込みの圧力センサからのＰsおよび吐出しの圧力センサからのＰdと、吸込みの温度センサからのＴsおよび吐出しの温度センサからのＴdとの測定値が必要である。障害状態では、ｈd、Ｐs、Ｐd、ＴsまたはＴdの測定値のうちの１つが欠けていると、測定された動作点４０１の割出しが妨げられ、フォールバック推定を遂行する必要がある。フォールバックステップ１３０の第２のサブステップ１３２中に遂行される動作は、本発明の第１の実施形態を参照しながら上記で説明されたものに類似であり、したがって詳細には報告されない。結果は添付の図８〜図１２に示されている。
【００４６】
図８〜図１２を参照すると、欠落データの推定および依然として有効な他の測定データに基づいて、マップ４００で、測定された動作点４０１が推定の動作点４０２によって置換される。欠落データを推定するのに用いられる計算および補間の誤差のマージンを考慮して、推定の動作点４０２は、円形の領域（図８〜図１０の、ｈd、ＰsまたはＰdが推定される場合）または測定された動作点４０１を含む細長い水平の領域（図１１、図１２の、ＴsまたはＴdが推定される場合）に入る。普通には、そのような領域は、ＳＬＬの右側の安全領域にあるか、または最悪のシナリオの手法で計算された動作点よりも安全領域に少なくともより接近しているはずである。公知の方法で用いられる最悪のシナリオでは、測定された動作点４０１は、マップ４００において、欠落データが、関連した可能な最大値または最小値に等しいと仮定することによって求められた最悪の場合のポイント４０３で置換され、最大値または最小値のうちどちらでも、ケースバイケースで最悪の状態を決定する。一般に、最悪の場合のポイント４０３はＳＬＬの左側にあり、サージ防止弁を開けさせる。
【００４７】
本発明の様々な実施形態（図示せず）によれば、他の無次元のマップを用いることができ、たとえば、流量エレメントＦＥが圧縮機１の吸込み側に配置されている場合には、ｈr対ｈs／Ｐsのマップが用いられ得る。しかしながら、すべての場合に、測定された動作点は、無次元のマップにおいて、本発明の第１の実施形態を参照しながら上記で説明されたものに類似の動作によって求められる推定の動作点で置換される。結果は、すべての場合において、添付の図４〜図６および図８〜図１２で図示されたものと同一であるかまたは類似であり、すなわち、推定の動作点は、ＳＬＬの右側の安全領域にあるか、または最悪のシナリオの手法で計算された動作点よりも安全領域に少なくともより接近しており、サージング防止制御システムの不必要な介在が防止され、その結果、サージ防止弁が不必要に開くのが防止される。
【００４８】
第３のステップ１１３中に遂行された検査がプラスであって、複数の測定器の障害が検出された場合、または、第５のステップ１１５中に遂行された検査がマイナスであって、測定器の障害が１つしか検出されないが、安全時間期間ｔ1が経過している場合には、方法１００は、さらなる置換の最悪の場合のステップ１４０で、測定器の障害のために欠けている２つ以上の測定値の最大値および／または最小値に基づいて、無次元のマップ３００、４００において、測定された動作点３０１、４０１または最悪の場合のポイント３０３、４０３を、推定の動作点３０２、４０２で置換する。たとえば、第１および第２の実施形態では、最悪の場合のポイント３０３、４０３は、ケースバイケースで、上記で定義され、添付の図４〜図６および図８〜図１２に表されたものである。最悪の場合のステップ１４０中に、特にステップ１４０が遂行されている圧縮機１のオペレータに対して伝えるために、警報が発生される。
【００４９】
第２の測定器がもはや信頼できないとき、すなわち圧縮機の挙動モデルに基づく推定がもはや不可能なとき、または第１の測定器の障害が安全時間ｔ1よりも長期にわたって残存するとき（許容できると見なされる）、最悪の場合のステップ１４０を実行すると、フォールバックステップ１３０に関して、より高度な安全性が確実になる。
【符号の説明】
【００５０】
１  圧縮機
３００  マップ
３０１  動作点
３０２  動作点
３０３  ポイント
４００  マップ
４０１  動作点
４０２  動作点
４０３  ポイント
ＦＥ  流量エレメント
【図１】
0006310930.tif 000002
【図２】
0006310930.tif 000003
【図３ａ】
0006310930.tif 000004
【図３ｂ】
0006310930.tif 000005
【図４】
0006310930.tif 000006
【図５】
0006310930.tif 000007
【図６】
0006310930.tif 000008
【図７ａ】
0006310930.tif 000009
【図７ｂ】
0006310930.tif 000010
【図８】
0006310930.tif 000011
【図９】
0006310930.tif 000012
【図１０】
0006310930.tif 000013
【図１１】
0006310930.tif 000014
【図１２】
0006310930.tif 000015