JPH01216099A - 遠心圧縮機の制御方法 - Google Patents

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 ［産業上の利用分野］ 本発明は、酸素製造プラントや各種プラントに　　′お
いて用いられる原料空気圧縮機、工場空気源用圧縮機、
化学プラント用ガス圧縮機等の遠心圧縮機の容量、圧力
を制御するための制御方法に関するものである。

［従来の技術］ 一般に、酸素製造プラントや各種プラントにおける遠心
圧縮機としては多段構成のものが用いられている。この
ような多段遠心圧縮機において、第１２図に示すように
、遠心圧縮機１は、駆動機２からの回転を増速する動力
伝達歯車３によって駆動される第１段目圧縮機４．第２
段目圧縮機５゜第３段目圧縮機６および第４段目圧縮機
７をそなえるとともに、圧縮機４，５間に中間冷却器８
を、圧縮機５，６間に中間冷却器９を、また圧縮機６゜
７間に中間冷却器１０をそなえて構成されている。

なお、圧縮機４と５および圧縮機６と７はそれぞれ同一
軸端にオーバハングされている。

このような遠心圧縮機１においては、第１段目の圧縮機
４に吸い込まれた空気は、各圧縮機５〜７および中間冷
却器８〜１０によって順次圧縮および冷却され、第４段
目の圧縮機７からプロセスに送出されるようになってい
る。

そして、各段の圧縮機４〜７の入口側には、角度可変式
の入口ガイドベーン（ＧＶ）１１〜１４が設けられ、こ
れらの入口ガイドベーン１１〜１４の角度を調整するこ
とにより、各圧縮機４〜７に流入する空気容量を調整で
きるようになっている。

また、各段の圧縮機４〜７の出口側には、ディフューザ
ベーン（ＤＶ）１５〜１８が設けられており。

これらのディフューザベーン１５〜１８の角度を調整す
ることによって、各圧縮機４〜７から流出する空気容量
を調整できるようになっている。

これらの入口ガイドベーン１１〜１４およびディフュー
ザベーン１５〜１８の角度は、それぞれ駆動装置１９に
よって任意の値に調整される。

さらに、この遠心圧縮機１全体もしくは各段の圧縮機４
〜７の運転状態、例えば、空気流量、温度、圧力等の運
転状態量は、それぞれ流量センサ２０、温度センサ２１
．圧力センサ２２などの検出手段によって検出されるほ
か１回転数センサ２４により遠心圧縮機１の回転数が検
出されるとともに、湿度センサ２７により遠心圧縮機１
に流入する空気の湿度が検出されるようになっている。

そして、各センサ２０〜２２，２４．２７と駆動装置１
９との間には、制御装置２３が設けられている。

上述のような遠心圧縮機を、各種の運転条件に応じて常
に所定の空気容量（流量）を最適運転効率で得られるよ
うに制御すべく、従来、第１３図に示すような制御手段
が開示されている（特開昭５５−６０６９２号公報）、
この制御手段においては、遠心圧縮機１全体もしくは各
段の圧縮機４〜７における空気流量、温度、圧力等によ
って表した各種の運転状態に対して最適運転状態を実現
するための操作量として、各段の入口ガイドベーン１１
〜１４およびディフューザベーン１５〜１８の角度の最
適な組合せ値を、予め制御袋＠２３内の記憶部にプログ
ラムして記憶させておく（ステップＳｏ）。

そして、第１３図に示すように、制御装置２３は、セン
サ２０〜２２から状態検出値を受けると。

その検出値から現在の多段遠心圧縮機の運転状態を演算
して監視する（ステップ８１）とともに、その運転状態
に対応し、予め決定しておいた操作プログラム（ベーン
角度の組合せ値）を照合し最適運転効率を実現できるよ
うな操作量を求めて、この操作量を駆動装置１９に出力
する（ステップＳ２）。

この後、センサ２０〜２２からの状態検出値をもとに演
算される運転状態が予め求められた最適運転効率状態で
あるか否か判定される（ステップＳ３）、最適運転効率
状態であると判定されると、その時点で制御を終了し選
択された操作プログラムでのベーン角度を維持する一方
、その運転状態が最適状態ではないと判定されると、選
択された操作プログラムによる制御出力（ベーン角度）
を修正して出力しなおしくステップＳ４）、状態検出値
から効率が向上したか否かを判定する（ステップＳ５）
。

このようにして、各段の運転状態をフィードバックせし
め、予めプログラムされた各段の入口ガイドベーン１１
〜１４およびディフューザベーン１５〜１８の角度の組
合せが最適か否かを監視し。

経年変化や性能変化等の運転状態の変化に対応して常時
最適運転効率が得られるように、操作プログラムにおけ
るベーン角度の組合せを修正する制御が行なわれる（ス
テップＳ６）。

ところで、多段遠心圧縮機を、各種の運転条件に応じ常
に所定の空気容量を最適運転効率で得られるように制御
するために、上述のごとく入口ガイドベーン１１〜１４
およびディフューザベーン１５〜１８の角度を制御する
手段のほかに、従来、第１４図に示すような制御手段も
開示されている（特開昭５６−６６４９０号公報）、こ
の制御手段は、入口ガイドベーン１１〜１４およびディ
フューザベーン１５〜１８の角度を制御する代わりに、
入口ガイドベーン１１（第１２図参照）の角度と駆動機
２（第１２図参照）の回転数とを制御するもので、第１
４図に示すように、入口ガイドベーン１１のみの角度を
調整する駆動袋［１９Ａと、駆動機２の回転数を制御す
るための駆動機制御袋［２５とがそなえられている。そ
して、遠心圧縮機１全体もしくは各段の圧縮機４〜７に
おける空気流量、温度、圧力等によって表した各種の運
転状態に対して最適運転状態を実現するための操作量と
して、入口ガイドベーン１１の角度と駆動機２の回転数
との最適な組合せ値の操作量（吸込流量に対応して、最
適な運転効率を得るための入口ガイドベーン１１の角度
と駆動機２の回転数とを与えたもの）を、予め、制御装
置２３Ａ内の記憶部２６にプログラムして記憶させてお
く。

そして、第１４図に示すように、制御装置２３Ａは、セ
ンサ２０〜２２．２４から状態検出値を受けると、その
検出値から現在の多段遠心圧縮機の運転状態（運転条件
や吸込流量など）を演算して監視する（ステップＴｌ）
とともに、その運転状態に対応し、特に運転条件の変化
によって吸込流量が変化した場合、記憶部２６における
操作量に基づいて、最適運転効率を実現できるような操
作量が求められ、その操作量が駆動装置１９Ａおよび駆
動機制御装置２５に出力される。

即ち、現時点における駆動機２および入口ガイドベーン
１１の操作状態が、記憶部２６における操作量と比較さ
れ（ステップＴ２）、その比較結果から最適運転効率状
態であるか否かが判定される（ステップＴ３）。

このとき、最適運転効率状態であると判定されれば、そ
の時点で制御を終了し操作状態（ベーン角度および回転
数）を維持する一方、最適運転効率状態ではないと判定
されると、入口ガイドベーン１１の角度および駆動機２
の回転数が、操作量に基づいて修正された後、得られた
操作量が駆動装置１９Ａおよび駆動機制御装置２５に出
力される（ステップＴ４）。

また、従来、上述した制御手段では、いずれもステップ
８３．Ｔ３において、最適運転状態であるか否かの判定
を行なうが、その際、各種のセンサ２０〜２２，２４．
２７による検出結果に基づき運転効率を演算する。この
演算に際しては、空気のガス定数Ｒを必要とするが、ガ
ス定数Ｒは。

第１５図に示すように、温度（ｔｏ〜ｔｚ）と湿度Ｈと
の関数になっているので、検出された温度および湿度に
基づきｉ正されたガス定数Ｒを用いて効率演算を行なっ
ている。

［発明が解決しようとする課題］ しかしながら、上述のような従来の遠心圧縮機の制御手
段では、入口ガイドベーン１１〜１４の角度とディフュ
ーザベーン１５〜１８の角度との組合せ、または、入口
ガイドベーン１１の角度と駆動機２の回転数との組合せ
を予め設定しておかなければならず、このような組合せ
設定のために。

遠心圧縮機の実際の稼動前もしくは稼動初期において相
当の試運転期間が必要であり、特に最適な組合せを設定
するには膨大なデータを要し、現実には設定不可能と考
えられる。また、温度、圧力等の検出可能な運転状態以
外の要因の変動によって、最適操作量（最適組合せ）が
変化した場合には。

それに対応することができない。

また、圧縮機の流量特性、効率特性、サージング発生条
件などは、すべて流入する空気の温度。

湿度の影響を受けるほか、実際には、効率特性は温度、
湿度以外の項目の影響も受けるので、効率特性について
のみ空気の温度、湿度を考慮しても、最高運転効率を保
証することはできない。

本発明は、これらの課題の解決をはかろうとするもので
、上述のような試運転期間や膨大なデータ等を必要とす
ることなく、運転状層の変化に対して自動的に追従でき
るようにして、常に最適な運転状態を維持できるように
した、遠心圧縮機の制御方法を提供することを目的とす
る。

［課題を解決するための手段］ 上記目的を達成するため、本発明の遠心圧縮機の制御方
法は、■入口ガイドベーンおよびディフューザベーンの
角度の関数でありこれらの角度に応じた遠心圧縮機の動
作特性を与え且つ所定のパラメータにより形状を決定さ
れる特性関数を予め設定するとともに、■上記遠心圧縮
機の複数種の運転状態に応じた上記特性関数の形状を決
定しうる複数のパラメータを、上記複数種の運転状態に
それぞれ対応して予め設定しておいてから、■上記遠心
圧縮機の運転状態を実測し、■実測された実運転状態が
上記複数種の運転状態と異なる場合には、■上記実運転
状態と上記複数種の運転状態との差異に応じ上記複数の
パラメータに重み付けをして上記実運転状態に対応した
推定パラメータを合成し、■この推定パラメータにより
上記実運転状態に応じた形状の上記特性関数を決定した
後、■決定された上記特性関数に基づいて上記遠心圧縮
機が所定の運転状態となるように上記の入口ガイドベー
ンおよびディフューザベーンのベーン角度操作量を決定
して上記遠心圧縮機を制御することを特徴としている。

［作　　　用］ 上述した本発明の遠心圧縮機の制御方法では、遠心圧縮
機における流量特性、効率特性、サージング発生条件等
の各種動作特性が、特性関数にて与えられ、この特性関
数の運転状態（温度、湿度。

圧力など）に応じた形状が、物理的に意味をもつパラメ
ータ（運転状態の関数）で定義される。そして、遠心圧
縮機が、運転状態の変化により、予め設定される複数種
の運転状態以外の運転状態となった場合には、上記複数
種の運転状態に対応して設定された複数のパラメータを
用い、特性関数の性質を保存しつつ運転状態のずれの量
に応じて、現在の運転状態（実運転状態）に応じた推定
パラメータが合成・補間される。

これにより、運転状態の変化が遠心圧縮機の特性に反映
され、従来では、温度、湿度など物理的なモデルが作成
できる場合に限ってしかその変化の影響を考慮できなか
ったが、本発明では、それら以外の要因、例えば物理モ
デルが明らかでない要因、理論値からのずれが著しい項
目などをも、特性関数の物理的特性を維持した状態で考
慮して、未知の状態における特性関数が推定され、推定
された特性関数に基づいて遠心圧縮機が所定の運転状態
となるように上記の入口ガイドベーンおよびディフュー
ザベーンのベーン角度操作量が決定され、遠心圧縮機が
制御される。

［発明の実施例］ 以下、図面により本発明の一実施例としての遠心圧縮機
の制御方法に、ついて説明すると、第１図はそのフロー
チャートであるが、まず１本実施例による方法を説明す
るに先立ち、第１１図により本実施例の方法が適用され
る遠心圧縮機の構成について説明しておく。なお、第１
１図において、第１２図中の符号と同一のものはほぼ同
様の部分を示しているので、その説明は省略する。ただ
し。

第１１図における本実施例の制御装置が適用される多段
遠心圧縮機では、圧縮機４〜７が全て同一軸上に配置さ
れ動力伝達歯車３が省略されている点が、第１２図にお
ける多段遠心圧縮機と異なっているが、第１２図に示す
多段遠心圧縮機にも本発明の方法を適用することはでき
る。

第１１図に示すように、本実施例における多段遠心圧縮
機では、入口ガイドベーン（ＧＶ）１１〜１４はそれぞ
れ入口ガイドベーン駆動装置１９ａ〜１９ｄにより駆動
されるとともに、ディフューザベーン（ＤＶ）１５〜１
８はそれぞれディフューザベーン駆動装ｆｌ　１９　ｅ
〜１９ｈにより駆動されるようになっている。また、セ
ンサとしては、従来と同様に、流量センサ２０．温度セ
ンサ２１．圧力センサ２２２回転数センサ２４．湿度セ
ンサ２７が設けられている。

そして、センサ２０〜２２，２４．２７からの検出信号
は、すべて制御装置２８へ入力されるようになっており
、この制御装置２８は、遠心圧縮機１における流量を制
御すべく、入口ガイドベーン１１〜１４およびディフュ
ーザベーン１５〜１８の角度をそれぞれ調整するよう各
検出信号を受けて各駆動装置ｉｉ　１９　ａ〜１９ｈへ
適当な制御信号を演算して出力できるようになっている
。

次に、上述のような遠心圧縮機に本実施例の制御方法を
適用する場合について、第１〜１０図を用いて説明する
。

まず、制御を開始する前に、制御装置２８の記憶部には
、入口ガイドベーン１１〜１４およびディフューザベー
ン１５〜１８の角度（α、β）により決まるベーン角度
平面（αβ平面）上において、ベーン角度（α、β）の
関数でありこれらのベーン角度（α、β）に応じた流量
Ｑおよび効率ηを与える流量特性関数Ｑ＝ｆ、（α、β
）および効率特性関数η＝　ｆｚ（α、β）がそれぞれ
予め数学的にモデル化されて設定されるとともに、サー
ジング領域（サージング発生条件；第８図中の斜線で示
す領域）を規定するサージングラインＳＬ（第８図参照
）もＧ＝１３（α、β）＝Ｏとして予め数学的にモデル
化されて設定されている。なお、流量特性関数ｆ１およ
び効率特性関数ｆ２は、それぞれ第８図の点線および実
線で示すようになっている。

特性関数ｆ１．　ｆ、、　ｆ、は、前述の通すベーン角
度（α、β）の関数であるが、その形状は、それぞれ所
定のパラメータｐ、、　ｐ２．　ｐ、によって決定でき
るため、 と表わすことができる０本実施例では、以下、効率特性
関数ｆ２を中心にして説明するが、流量特性関数ｆ□や
サージング発生条件ｆ３も同様に取り扱われる。

ここで、効率特性関数ｆ２の形状を、例えば、第２図に
示すように１等効率線を楕円として近似し、効率特性関
数ｆ２をα、βの２次式で表す。

この場合、パラメータＰ２は、第２図に示すように・α
ａｔ　βｏｓ　Ｃ１９σ２ｔ　θ、Ｑの６つの物理的に
意味のある量【後述の（２）式参照；ｐｚ＝（α。。

β。、σ１．σ２．θ、Ｃ）Ｔ）であり、これらによっ
て効率特性関数ｆ２の位置、大きさ、形状などが表され
る。

さらに、（１）式により示した各特性に影響を与える項
目は、すべてパラメータＰユ　ｐ２．ｐ：ｌに含まれる
ことになるので、遠心圧縮機１の運転状態Ｃに応じて、 と表わせる。ここで、運転状態Ｃは、温度センサ２１に
より検出される温度Ｔ、圧カセンサ２２により検出され
る圧力Ｐ、湿度センサ２７により検出される湿度Ｈなど
であり、Ｃ＝（Ｔ、Ｐ、Ｈ，・・・）Ｔと表わせる。

そして、本実施例では、複数種の既知の運転状態Ｃ１［
ｉ＝１．２．・・・；温度、圧力、湿度などの上下限値
、平均値あるいは特徴的な値等の点を用゛、）る；具体
的な例については第５，６図にて詳述する］について、
各運転状態Ｃ１に応じた各特性関数ｆ□＊　Ｌ＊　ｆ、
の形状を決定するパラメータＰ８゜Ｐ、、Ｐ、を、各運
転状態Ｃ１ごとに（２）式に基づきｐｘ（ｃｉＬ　ｐｚ
（ｃｉＬ　Ｐｉ（Ｃｉ）（ｉ＝Ｌ　Ｌ・・・）として、
予め実測により設定しておき、制御を開始する前に、制
御装置２８の記憶部に（１）式により示す関数とともに
記憶させておく。

このようにして特性関数ｆ１．ｆ２．ｆ、および運転状
態ＣＩに応じた設定パラメータＰ□（Ｃ１）。

ｐｇ（ｃｉＬ　Ｐｉ（Ｃｉ）　（ｉ：ｉ、　２．・・・
）を予め与えておいてから、本実施例における制御方法
では、第１図に示すように、ステップＡ２〜Ａ７により
本発明の特徴とする制御方法を実施した後、ステップ八
８〜Ａ１７による山登り法と呼ばれる制御方法を実施し
て、制御装置２８により遠心圧縮機１の流量を制御して
いる。なお、後述するように。

ステップＡ１８〜Ａ２０はサージング解消方法を実施す
るものであるとともに、ステップＡ２１〜Ａ２３はサー
ジング回避方法を実施するものである。

さて、遠心圧縮機１の流量を制御する際には、まず、目
標流量ＱＦＸＮを設定して（ステップＡｌ）、遠心圧縮
機ｌの流量の制御が開始される。この後、まず、流量セ
ンサ２０．温度センサ２１．圧力センサ２２１回転数セ
ンサ２４．湿度、センサ２７により現在の実運転状態Ｃ
を実測しくステップＡ２）、この実運転状態Ｃが、予め
設定されている既知の運転状態Ｃ１（ｉ＝１．２．・・
・）に等しいか否かを判定する（ステップＡ３）。

判定の結果、Ｃ＝Ｃ，であれば、その運転状態に応じた
パラメータｐ、、ｐ、、ｐ３は、前述の通り設定パラメ
ータｐｔ（ｃｉＬ　Ｐｚ（Ｃｉ）−Ｐｉ（Ｃｔ）（ｉ＝
１．２．・・・）として予め設定されているので、これ
らの設定パラメータの中から適切なものをそれぞれ選択
し、運転状態Ｃに応じたパラメータＰ１．Ｐ、、Ｐ、を
決定する（ステップＡ４）。

一方、ステップＡ３における判定の結果、Ｃ≠Ｃ１であ
れば、実運転状態Ｃは、既知運転状態Ｃ１と異なる未知
の運転状態であるので、特性関数ｆ□ｔ　ｆｚｔ　ｆ２
の形状を決定するためのパラメータｐ、、ｐ、、ｐ、も
未知である。そこで、以下のうようにして、未知の実運
転状態Ｃに対応するパラメータを推定する。

まず、未知のパラメータＰ、（Ｃ）、Ｐ、（Ｃ）、Ｐ。

（Ｃ）を推定するために、実測された実運転状態Ｃと既
知の運転状態Ｃ１との差異に応じて、適当な設定パラメ
ータ〔例えば、実運転状態Ｃに近い既知の運転状態Ｃ１
についての設定パラメータ〕を、ｐｘ（ｃｉＬ　Ｐａ（
ＣｉＬ　Ｐ３（Ｃ１）（ｉ＝＝１．２１・・・）からそ
れぞれ選択する（ステップＡ５）。ここでは、例えば、
２つの既知運転状態Ｃ１，Ｃ２についての設定パラメー
タｐ、（ｃ、）、ｐ、（ｃｚ）が選択されたものとして
、パラメータＰ！（Ｃ）を推定し効率特性関数ｆ２を決
定する場合を説明する。

即ち、設定パラメータＰａ（Ｃｔ）＊　ｐｚ（ｃａ）の
選択後、これらの設定パラメータと運転状態Ｃ，Ｃ，。

Ｃ２とに基づいて、各設定パラメータｐ、（ｃｍ）。

ｐｉ（ｃａ）に重み付けをして、実運転状態Ｃに対応し
た推定パラメータＩ’ａ（Ｃ）が合成される（ステップ
Ａ６）、なお、他のパラメーター工（Ｃ）、ｐａ（ｃ）
についての推定パラメーターｚ　（ｃ　）　−Ｐ　３　
（ｃ　）も同様にして合成される。

例えば、未知のパラメータＰｓ（ｃ）は、次式（３）で
示すような推定パラメータｆ）、（Ｃ）として合成され
る。

１）、（Ｃ）＝ｇｔ（Ｃ）・ｐｚ（ｃｘ）＋ｇｔ（Ｃ）
・Ｐ　２　（Ｃ２）・・・（３） ただし、 この（３）式は、選択された設定パラメータＰ。

（ｃ　ｚ　）　、ｐ　ｚ　（ｃ　２）に、基準状態から
のずれに応じた重みｇ、（Ｃ）およびｇｉ（Ｃ）をそれ
ぞれ付けて補間して、推定パラメータＩ）、（Ｃ）を合
成している。

このようにして、ステップＡ４により選択されたパラメ
ータｐｔ（ｃｉＬ　ｔ’５（ｃｉＬ　Ｐａ（Ｃｉ）、も
しくは、ステップＡ５により合成された推定パラメータ
ー１（ＣＬｆ□（Ｃ）、Ｉ’３（Ｃ）を、（１）式にお
ける各式に代入することにより、実運転状態Ｃに応じた
特性関数Ｑ、η、Ｇが決定される（ステップＡ７）。

ところで、（３）式により合成して得られた推定パラメ
ータ］）２（Ｃ）に基づき形状を決定される具体的な効
率特性関数の例を、第３図（ａ）〜（ｄ）により説明す
る。

選択された設定パラメータＰ　２　（Ｃｉ　）　−Ｐ　
Ｚ　（Ｃ２）によりそれぞれ決定される効率特性関数η
＝ｆ２（α、β、ｐ、（ｃｍ））、η＝ｆ２（α、β−
ｐ２（ｃｚ））の形状が、それぞれ第３図（ａ）、（ｂ
）に示すようなものである場合、一般的な関数の線形近
似（ｆ。

（α、β、　ＰＺ（Ｃ））＝（ｆ！（α、β、　Ｐａ（
ＣＩ））＋ｆｚ（α、β、Ｐｉ（ｃ、）））／２；第４
図参照〕を行なうと、第３図（ｃ）に示すような形状を
有する効率特性関数が得られる。一方、（３）式による
パラメータ補間を行なって推定パラメータｐｚ＜ｃ＞を
求めてから、この推定パラメータｐｇ（ｃ）を（１）式
に代入すると、第３図（ｄ）に示すような形状の効率特
性関数が得られることになる。

第３図（Ｑ）と第３図（ｄ）とに示す効率特性関数を比
較して明らかなように、線形近似による第３図（Ｃ）の
効率特性関数は、第３図（ａ）、（ｂ）に示すものの中
間状態を近似したものとはいえず、効率特性関数の物理
的特性を維持していないが、パラメータ補間による第３
図（ｄ）の効率特性関数は、効率特性関数の物理的特性
を十分に維持しており。

第３図（ａ）、（ｂ）に示す特性関数の中間状態をよく
近似している。

次に、さらに第５〜７図により、既知の運転状態ＣＩ（
ｉ＝１．２．・・・）の設定例や、第３図および（３）
式により示したものとは異なるパラメータ補間手段の例
を第１図に対応させて具体的に説明する。

例えば、運転状態Ｃの量として、遠心圧縮機１の入側温
度Ｔと気圧Ｐとの２個を考える。つまり、Ｃ＝（Ｔ、Ｐ
）”とする。

ここで説明する例においても、前述した場合と同様に、
オフライン状態で予め何通りかの既知の運転状態Ｃｔ（
ｉ＝１．２．・・・）について、特性関数ｆ１．　ｆ２
．　ｆ、の形状を決定するｐｘ（ｃｉＬ　Ｐ。

（Ｃｔ）、　Ｐｉ（Ｃｉ）　（ｉ　＝Ｌ２ｔ・・・）を
実測により設定しておくが、既知の運転状態Ｃｉ（ｚ　
＝　１１２　ｅ・・・）の選び方は、例えば、第５図に
示すように、運転状態として考えられる最低温度Ｔしと
最高温度Ｔ。

どの間、および、最低気圧ｐＬと最高気圧ＰＨとの間を
それぞれ何等分（第５図ではｎ等分）かして得られる格
子点を、運転状ｊｉｔｃｔとして選ぶ。

そして、実際に遠心圧縮機１を制御する際には、ステッ
プＡ２で実運転状態Ｃ＝（Ｔ、Ｐ）Ｔを実測し、ステッ
プＡ３でＣ＃Ｃ１と判定された場合、ステップＡ５にお
いて、予め設定されたパラメータｐｉ（ｃｉ）ｔ　Ｐａ
（ＣｘＬ　Ｐｌ（Ｃ１）　（ｉ＝　１ｔ　２１・・・）
のうちどれを用いて未知のパラメータＰｌ（Ｃ）−ｐ、
（ｃ）、ｐ、（ｃ）を推定するかを決定する。例えば、
第５図７中の斜線部分に実運転状態Ｃがある場合には、
第５図に示す既知運転状態ＣＡ＝（Ｔ工。

Ｐ　１）ｅ　Ｃａ＝　（Ｔｔ　＊　Ｐ　ａＬ　Ｃｃ＝　
（’ｒ　ｚ　ｅ　Ｐ　ｘ）　ｅ　Ｃ。

：（Ｔ、、Ｐｉ）の４点を選択する。

ステップ八６では、ステップＡ５にて選択した各既知運
転状態に対応した設定パラメータＰ２（ＣＡ）＝　Ｐ　
ｚ（Ｔｔｅ　Ｐ□）ｅ　Ｐ　１（Ｃａ）＝　Ｐ　ｚ（Ｔ
ｔｅ　Ｐ　ｚ）ｔＰ　１（Ｃｃ）＝　Ｐａ（Ｔｔｅ　Ｐ
　ｔＬ　Ｐ　ａ（Ｃｏ）＝　Ｐ　ｚ（Ｔ２ｅ　Ｐ　ａ）
から、実運転状態Ｃ＝（Ｔ、Ｐ）”における未知のパラ
メータｐ＊（ｃ）を合成により推定する。なお。

ここでは、パラメータｐｇ（ｃ）を推定する場合につい
てのみ説明するが、他の未知のパラメータＰ工（Ｃ）、
Ｐｌ（Ｃ）も同様にして推定される。

合成の手段としては、前述した（３）式を拡張してこの
場合に適用すると１次式（４）により推定パラメータ？
、（ｃ）が合成される。

・・・（４） この（４）式において、（Ｔ−Ｔ、）／（Ｔ１−Ｔ、）
と（’ｒ−’ｒ、）／（’ｒｔ−’ｒ□）とは、それぞ
れ既知点Ｔ２およびＴ２から未知点Ｔへの温度差の比率
となっているとともに、（ｐ−ｐ　＊　）　／　（ｐ　
ｔ−Ｐ　２　）と（ｐ−ｐｔ）／（ｐｉ−ｐｉ）とは、
それぞれ既知点Ｐ２およ・　びＰｌから未知点Ｐへの気
圧差の比率となっておす、これらの比率によって選択さ
れた設定パラメータの重み付けがなされ、推定パラメー
タｐｚ＜ｃ＞が補間・合成される。

（４）式により合成して得られた推定パラメータｆ’、
（Ｃ）に基づき形状を決定される具体的な効率特性関数
の例を、第７図（ａ）〜（ｉ）に示す。第７図（ａ）〜
（ｄ）はそれぞれ選択された設定パラメータＰｉ（ＣＡ
）−Ｐｘ（Ｃａ）−Ｐａ（Ｃｃ）ｅ　Ｐｇ（Ｃｏ）によ
り決定される効率特性関数ｆ２の形状を示しており、第
７図（ｅ）〜（ｉ）は、それぞれ、第６図に示す未知の
運転状態Ｃ＋：ｙＣｒｐＣａｐＣｏ＊Ｃ工について、設
定パラメータＰｇ（ＣＡ）ｔ　ｐｉ（ｃａ）ｅ　Ｐｉ（
Ｃｃ）＊Ｐｇ（Ｃｏ）により（４）式に基づいて推定し
たパラメータＰｇ（Ｃｉ：）ｓ　Ｐａ（Ｃｒ）ｓ　Ｐｇ
（Ｃａ）＊　Ｐｇ（ＣｕＬＰ　ｚ　（Ｃｘ）によって決
定される効率特性関数ｆ２の形状を示している。

なお、（４）式では、入側温度Ｔと気圧Ｐとを考慮して
いるが、これに湿度Ｈなどの他の計測量を加えて、実運
転状態Ｃ＝（Ｔ、Ｐ、Ｈ，・・・）としてもよく、この
場合には、推定パラメータＰｇ（Ｃ）は次式（５）のよ
うになる。

Ｔ、−Ｔユ　Ｐｌ−Ｐ、　　Ｈｌ−Ｈ。

Ｔｉ−Ｔｔ　　Ｐｉ−Ｐｉ　　Ｈｚ−Ｈｚ・・・（５） 上述のようにして、各特性関数ｆ１．　ｆ、、　ｆ３の
形状を実運転状態Ｃに応じて決定するパラメータＰユ（
Ｃ）、Ｐｉ（Ｃ）、Ｐｌ（Ｃ）が、ステップＡ４もしく
はＡ６において選択もしくは推定により得られ、ステッ
プＡ７で（１）式に基づいて実運転状態Ｃに応じた特性
関数Ｑ、η、Ｇが決定された後には、山登り法によるベ
ーン角度制御を実施する。

即ち、まず、入口ガイドベーン１１〜１４およびディフ
ューザベーン１５〜１−８の現在角度（α。

β）を測定する（ステップＡ８）、そして、測定された
ベーン角度（α、β）から、この現在角度が予め設定さ
れたサージング領域内にあるか否かを判定しくステップ
Ａ９）、サージング領域内にある場合［Ｇ（α、β）＜
Ｏのとき］には後述するステップＡ１８に移る。一方、
サージング領域内にない場合には次のステップＡＩＯへ
移る。このとき、上記ベーン角度（α、β）における現
在流量Ｑ　ｐｉｏｗは。

ステップＡ７で決定された流量特性関数ｆ１により予測
する。

現在のベーン角度（α、β）がサージング領域内にない
場合には、まず、目標流量ＱＦＩＮへ向かう第１のベー
ン角度操作量ベクトル（Δα１．Δβ１）を次のように
して求める。つまり、流量特性関数ｆ１により決定され
る現在のベーン角度（α、β）における等流量曲線Ｑ（
α、β）＝ＱＮｏｖについて、この等流量曲線の法線方
向（ａｆ□／ａα、ａｆｉ／ａβ）で目標流量ＱＦＩＮ
側へ向かう方向にベーン角度（α、β）を駆動操作すれ
ば、最短操作経路が得られるので、この法線方向ベクト
ルであり目標流量ＱＦＩＮへ向かうものを第１のベーン
角度操作量ベクトル（Δα□、Δβ、）として求める。

ここで、ベーン角度操作量ベクトル（Δα１．Δβ、）
に基づいて修正される流量をΔＱとすると。

Ａ　Ｑ　＝（ａ　ｆ　−／　ａ　ａ　）　Ａ　ａ　ｘ　
＋　（ａ　ｆ　ｘ　／ａ　ｆ’　）　Ａ　ｆＪ　１・・
・・（６） と表され、ベーン角度操作量Δα４．Δβ１は、目標流
量ＱＦＩＮと現在流量Ｑ　Ｎ０ＩＩ＋との差に比例する
ので、上記法線方向も考慮すると、 ・・・・（７） として、第１のベーン角度操作量ベクトル（Δα、。

Δβ、）を求めることができる。なお、Ｋ、はチューニ
ングパラメータを示す。

このようにして、ステップＡＩＯにおいて第１のベーン
角度操作量ベクトル（Δα１．Δβ１）が求められた後
、ＩＩｔ在のベーン角度（α、β）が第８図に示すサー
ジングラインＳＬの近傍にあるか否かを判定する（ステ
ップＡ１１）、そして、サージングラインＳＬの近傍に
ある場合［Ｇ（α、β）＜ｅａのとき；ただしｅＱはＯ
に近い正の値］には後述するステップＡ２１に移る一方
、サージングラインＳＬの近傍にない場合には次のステ
ップＡ１２に移る。

ステップＡ１．２においては、目標流量ＱＦ工Ｎに最短
経路で近づくように求められた第１のベーン角度操作量
ベクトル（Δα、、Δβ□）に対して、その流量を変化
させることなく、効率が増加するように第２のベーン角
度操作量ベクトル（Δα２゜Δβ２）を求める。つまり
、前述したステップＡｌｌでの流量制御の場合と同様に
、効率特性関数ｆ２により決定される現在のベーン角度
（α、β）における等効率面纏り（α、β）；ηＮＯＶ
について、この等効率曲線の法線方向（ａｆ、／ａα、
ａｆ、／ａβ）で且つ改善効率Δηが正となる方向にベ
ーン角度（α、β）を駆動操作すれば、効率ηは改善さ
れる。

ここで、第２のベーン角度操作量ベクトル（Δα２゜Δ
β２）に基づいて改善される効率Δηは、Δη＝（ａｆ
、／ａα）Ａ　ｒｘ、＋Ｃａｆ、／ａβ）Ａ　ｆＪ２・
・・・（８） と表される。このとき、効率改善の操作によって流量が
変化するのを防止するために、ベーン角度操作量Δα２
．Δβ２による流量変化ΔＱをゼロとすべく、（６）式
より。

（ａｆ、／ａａ）ｈａ、＋Ｃａｆ１／ａβ）Δβ２＝０
・・・・（９） を満足するようにベーン角度操作量Δα２．Δβ２を決
定する。即ち、上記等流量曲線の接線方向ベクトルで改
善効率Δη〉Ｏとなるものを、（８）。

（９）式から求めると。

・・・・（ｌＯ） として、第２のベーン角度操作量ベクトル（Δα２゜Δ
β２）を求めることができる。なお、Ｋ２はチューニン
グパラメータを示す。

そして、第１のベーン角度操作量ベクトル（Δα□、Δ
β、）および第２のベーン角度操作量ベクトル（Δα２
．Δβ２）を加算することにより、流量を目標流量ＱＦ
ＩＮに近づけながら効率を高めることができる第３のベ
ーン角度操作量ベクトル（Δα３．Δβ３）を求めベー
ン角度操作量ベクトル（Δα、Δβ）として出力する（
ステップＡ１３）。

つまり、下記（１１）式に基づく演算・を行なう。

この後、求められた各ベーン角度操作量Δα。

Δβが、入口ガイドベーン１１〜１４およびディフュー
ザベーン１５〜１８の最大許容操作量を越える場合、つ
まり１回で駆動することができる角度を越えている場合
があるので、ベーン角度操作量ΔαまたはΔβを、第９
図に示すようなリミッタ関数ＦＬ（Ｘ）に代入しておき
、上記最大許容操作量を越える場合にはベーン角度操作
量ΔαまたはΔβを縮小するようにする（ステップＡ１
４）。

さらに、このステップＡ１４までに求められたベーン角
度操作量ベクトル（Δα、Δβ）に基づいてベーン角度
（α、β）を駆動操作した場合に、１回の操作で変更さ
れた新たなベーン角度（α＋Δα。

β＋Δβ）がサージング領域内に入ると判断されたとき
は、このベーン角度操作量ベクトル（Δα。

Δβ）の大きさを縮小して出力する（ステップＡ１５）
。

このようにして得られたベーン角度操作量ベクトル（Δ
α、Δβ）に基づいて制御装置ｉ！２８から駆動装置１
９ａ〜１９ｈへ制御信号を出力し、入口ガイドベーン１
１〜１４およびディフューザベーン１５〜１８を駆動制
御する（ステップＡ１６）。

そして１以上のベーン駆動制御により変更された流量が
目標流量Ｑ　ＦＩＮの近傍に到達したか否かを判断しく
ステップＡ１７）、到達していない場合には、再びステ
ップＡ２に戻り、実運転状態Ｃに応じた特性関数の決定
および出登り法によるベーン角度制御を続行する一方、
目標流量ＱＦＩＮの近傍に到達したと判断された場合に
はその時点で流量制御を終了する。ここで、ステップＡ
１７における判断方法としては次のようなものがある。

（ｉ）上記ベーン駆動制御により得られた現在のベーン
角度（α、β）に基づいて、予め設定された流量特性関
数ｆ１から現在流量を予測し、この予測値Ｑと目標流量
ＱＦ工Ｎとを比較して、両者の差が一定値と８以内、即
ち。

Ｉ　Ｑ　−ＱＦｚｓｌ≦ｃ８ となったならば、目標流量ＱＦＸＮの近傍に到達したと
判断する。

（ｉｆ）ベーン角度操作量Δα、Δβが、ｌΔα１≦６
かつ１Δβ１≦ｉ４となった場合（それ以上操作されな
くなった場合）に、目標流量ＱＦＸＨの近傍に到達した
と判断する。

（ｉｎ）ベーン角度操作量Δα、Δβが。

ｈａ　　＋τ了Ｐ≦ｆｓ となった場合に、目標流量Ｑ　ＦＩＮの近傍に到達した
と判断する。

以上のようにして、ステップ八８〜Ａ１７による山登り
法により、ベーン角度の始点（α。、β、）から目標流
量ＱＦＩ関の近傍に到達するとともに高効率となるよう
に、ベーン角度制御が行なわれる。

さて、次にステップＡ１８〜Ａ２０によるサージング解
消方法およびステップＡ２１−Ａ２３にサージング回避
方法について説明する。まず、ステップＡ９において現
在のベーン角度（α。、β。）がサージング領域内にあ
ると［Ｇ（α、β）＜Ｏのとき］判定された場合には、
第１０図（ａ）に示すように、サージングを解消するた
めの最短操作経路は、当然、サージラインＳ　Ｌ［ｇ−
Ｇ（α、β）＝０］の法線方向（ａｇ／ａα、ａｇ／ａ
β）”ｒサージング領域外へ向かう方向であるので、こ
の法線方向ベクトルでありサージング領域外へ向かうも
のをサージング解消用ベーン角度操作量ベクトル（Δα
４．Δβ４）として求める。ここで、操作量は境界線Ｇ
（α、β）＝０と現在のサージング関数値Ｇ（α。、β
。）との差に比例するので。

・・・・（１２） として、サージング解消用ベーン角度操作量ベクトル（
Δα、、Δβ４）を求め（ステップＡ１８）、これをベ
ーン角度操作量ベクトル（Δα、Δβ）として出力する
（ステップＡ１９）。

そして、ステップＡ１４の場合と同様に、求められたベ
ーン角度操作量Δα、Δβが、入口ガイドベーン１１〜
１４およびディフューザベーン１５〜１８の最大許容操
作量を越える場合には、ベーン角度操作量ΔαまたはΔ
βを縮小してから、駆動装置１９ａ〜１９ｈへ制御信号
を出力して。

入口ガイドベーン１１〜１４およびディフューザベーン
１５〜１８を駆動制御する（ステップＡ２０）、この後
、再びステップＡ８へ戻りステップＡ９により現在のサ
ージング領域内にあるか否かを判定し、ベーン角度（α
、β）がサージング領域外となるまで、上述したステッ
プＡ１８〜Ａ２０の手順を繰り返す。

また、ステップＡｌｌにおいて、現在のベーン角度（α
、β）がサージングラインＳＬの近傍にあると［Ｇ（α
、β）＜ｅａのとき；ただしｅＱは０に近い正の値］判
定された場合には、まず、ステップＡＩＯにより求めら
れた第１のベーン角度操作量ベクトル（Δα１．Δβ１
）が、第１０図（ｂ）に示すように、サージング領域内
に向かっているか否かを判定しくステップＡ２１）、向
かっていない場合には、第２のベーン角度操作量ベクト
ルを演算することなくステップＡ１２およびＡ１３をと
ばして、第１のベーン角度操作量ベクトル（Δα１゜Δ
β１）をそのままベーン角度操作量ベクトル（Δα、Δ
β）として出力しくステップＡ２３）ステップＡ１４へ
移ることにより、効率の改善を考慮することなく逸早く
サージングラインＳＬの近傍から回避することができる
。

一方、ステップＡ２１において、第１のベーン角度操作
量ベクトル（Δα４．Δβ１）が、第１０図（ｂ）に示
すように、サージング領域内に向かっていると判定され
た場合には、次のようにして第１のベーン角度操作量ベ
クトル（Δα４．Δβ□）からサージング領域へ向かう
成分を除去することにより、サージングラインに沿うサ
ージング回避用ベーン角度操作量ベクトル（Δα６．Δ
β、）を求めて、サージングを回避する。ここで、サー
ジング領域内に向かうサージングラインＳＬの法線方向
単位ベクトルをｎＱとすると。

・・・・（１３） であり、この（１３）式を用いて第１のベーン角度操作
量ベクトル（Δα、、Δβ１）から次式（１４）に示す
ようなサージング回避用ベーン角度操作量ベクトル（Δ
α６．Δβ、）が求められる。

・・・・（１４） であり、Ｋ、はチューニングパラメータである。

このようにして、第１のベーン角度操作量ベクトル（Δ
α□、Δβｉ）からサージングラインＳＬに沿う方向の
サージング回避用ベーン角度操作量ベクトル（Δα５．
Δβ、）を求め（ステップＡ２２）、このベーン角度操
作量ベクトル（Δα５．Δβ、）をベーン角度操作量ベ
クトル（Δα、Δβ）として出力しくステップＡ２３）
ステップＡ１４へ移ることにより、効率の改善を考慮す
ることなく、且つ、サージング領域に入るのを回避しな
がら、遠心圧縮機１の流量制御が行なわれるのである。

なお、ステップＡ１４以降の手順は前述の通りであるの
で、その説明は省略する。

上述のようにして、サージングを解消あるいは回避しな
がら、山登り法により、ステップＡ２〜Ａ７により決定
されたベーン角度操作量の特性関数ｆ、、ｆ、、ｆ、に
従って、目標流量ＱＦＩＮへの流量制御が行なわれるの
である。

以上のように、本実施例の方法によれば、遠心圧縮機１
における流量特性Ｑ、効率特性η、サージング発生条件
Ｇの各種動作特性が、それぞれ特性関数ｆ１．　ｆ、、
　ｆ、で与えられ、実運転状態Ｃ（温度、湿度、圧力な
ど）に応じた形状の特性関数ｆ１．　ｆ、、　ｆ、が、
物理的に意味をもつパラメータｐｍ（ｃ）、ｐｇ（ｃ）
、ｐａ（ｃ）を補間により推定することで決定される（
ステップＡ５．Ａ６）。

即ち、遠心圧縮機１が、運転状態Ｃの変化により、予め
設定された複数種の運転状態ＣＩ（ｉ＝１゜２、・・・
）以外の運転状態となると、複数種の運転状態Ｃ１に対
応する設定パラメータＰ工（Ｃ１）、ｐ。

（Ｃｉ）　ｅ　Ｐ　ｚ　（Ｃｉ）を用いて、（３）〜（
５）式により各特性関数ｆユ＃　ｆｌｌｌ　ｆａの性質
を保存しつつ運転状態Ｃのずれの量に応じて、実運転状
ｓＣに対応する推定パラメータｆ’１（Ｃ）、かｚ（ｃ
）、ｐａ（ｃ）が補間・合成され、これを（１）式に代
入することにより実運転状態Ｃに応じた形状の特性関数
ｆｉｌｒ、、ｆ、が決定される。

これにより、運転状態Ｃの変化が遠心圧縮機１の特性に
反映され、従来では、温度、湿度など物理的なモデルが
作成できる場合に限ってしかその変化の影響を考慮でき
なかったが、本実施例では、それら以外の要因、例えば
物理モデルが明らかでない要因、理論値からのずれが著
しい項目などをも、特性関数ｆユＩ　ｆａｔ　ｆａの物
理的特性を維持した状態で考慮し、未知の運転状態Ｃに
おける特性関数ｆｔ＝　ｆｚ−ｆａが推定され、それに
基づいて遠心圧縮機１が所定の流量ＱＦＩＮとなるよう
に入口ガイドベーン１１〜１４およびディフューザベー
ン１５〜１８のベーン角度操作量が決定され、出登り法
によって常に最適な運転状態を維持しながら遠心圧縮機
１が制御される。

従って、（１）式に対応する基本的な特性関数ｆｉｔｆ
、、　ｆ、と複数種の運転状態Ｃ１に対応する設定パラ
メータｐｘ（ｃｉＬ　Ｐｇ（Ｃｔ）＊　Ｐａ（Ｃｉ）と
を予め設定しておくことで、従来のように相当の試運転
期間や膨大なデータを必要とすることなく、運転状態Ｃ
の変化に追従することができるので、最適な運転状態を
維持できて、常に最高運転効率が保証されるのである。

なお、上記実施例の制御方法を実施するに際しては、全
段の入口ガイドベーン１１〜１４およびディフューザベ
ーン１５〜１８を各段ごとに上記制御方法に別々に制御
してもよいし、全段の入口ガイドベーン１１〜１４およ
びディフューザベーン１５〜１８を同一角度で駆動制御
するようにしてもよい、また、入口ガイドベーン１１〜
１４およびディフューザベーン１５〜１８の角度を、各
段の遠心圧縮機４〜７の運転流量が設計流量に対して同
−比の相似運転流量となる１組の無次元入口ガイドベー
ン角度および無次元ディフューザベーン角度としてそれ
ぞれ表し、各段の入口ガイドベーン１１〜１４およびデ
ィフューザベーン１５〜１８の角度を１組のものとみな
して、ベーン角度制御を行なうようにしてもよい。

また、上記実施例では、本発明の方法を多段遠心圧縮機
に適用した場合を示しているが、単段の遠心圧縮機にも
同様にして適用することができるほか、本実施例の方法
は流量制御だけでなく圧力制御にも適用できる。

［発明の効果］ 以上詳述したように、本発明の遠心圧縮機の制御方法に
よれば、運転状態が予め設定された複数種の運転状態以
外の運転状態となると、上記複数種の運転状態に対応す
る設定パラメータを用いて、運転状態のずれの量に応じ
、実運転状態に対応する推定パラメータを補間・合成し
、この推定パラメータにより形状を決定された特性関数
に基づいて遠心圧縮機の制御を行なうように構成したの
で。

運転開始時に特性関数と複数種の運転状態に対応する設
定パラメータとを与えてしまえば、その後は、自動的に
運転状態のの変化に追従することができ、最適な運転状
態を維持でき、常に最高運転効率を保証できる効果があ
る。

【図面の簡単な説明】

第１〜１１図は本発明の一実施例としての遠心圧縮機の
制御方法を示すもので、第１図はそのフローチャート、
第２図はその効率特性関数についてのパラメータを説明
するためのグラフ、第３図（ａ）〜（ｄ）は本発明の方
法の作用を説明するための数値計算例を示すグラフ、第
４図は一般的な線形補間を説明するためのグラフ、第５
，６図は複数種の運転状態の具体的な設定例を示す図、
第７図（ａ）〜（ｉ）は本発明の方法の作用を説明する
ための数値計算例を示すグラフ、第８図は予め設定され
る流量特性関数、効率特性関数およびサージングライン
のモデル例を示すグラフ、第９図はそのリミッタ関数の
例を示すグラフ、第１０図（ａ）はサージング解消方法
を説明するためのグラフ。 第１０図（ｂ）はサージング回避方法を説明するための
グラフ、第１１図は本実施例の方法を適用する遠心圧縮
機を示すブロック図であり、第１２図は従来の一般的な
多段遠心圧縮機を示すブロック図、第１３図は従来の遠
心圧縮機の制御手段を説明するためのフローチャート、
第１４図は他の従来の遠心圧縮機の制御手段を説明する
ためのフローチャート、第１５図はガス定数と湿度およ
び温度との関係を示すグラフである。 図において、１−遠心圧縮機、２−駆動機、４〜７−圧
縮機、８〜１〇−中間冷却器、１１〜１４−人口ガイド
ベーン、１５〜１８・・−ディフューザベーン、１９８
〜１９ｄ−人口ガイドベーン駆動装置、１９ｆ〜１９ｈ
−ディフューザベーン駆動装置、２〇−流量センサ、２
１・−温度センサ、２２−圧力センサ、２４−回転数セ
ンサ、２７−・・湿度センサ、２８−制御装置。 特許出願人　株式会社　神戸製鋼所

Claims (1)

【特許請求の範囲】
1. 入口側および出口側にそれぞれ角度可変式の入口ガイド
   ベーンおよびディフューザベーンを有する遠心圧縮機に
   ついて、上記の入口ガイドベーンおよびディフューザベ
   ーンのベーン角度を操作して上記遠心圧縮機を制御する
   方法において、上記の入口ガイドベーンおよびディフュ
   ーザベーンの角度の関数でありこれらの角度に応じた上
   記遠心圧縮機の動作特性を与え且つ所定のパラメータに
   より形状を決定される特性関数を予め設定するとともに
   、上記遠心圧縮機の複数種の運転状態に応じた上記特性
   関数の形状を決定しうる複数のパラメータを、上記複数
   種の運転状態にそれぞれ対応して予め設定しておいてか
   ら、上記遠心圧縮機の運転状態を実測し、実測された実
   運転状態が上記複数種の運転状態と異なる場合には、上
   記実運転状態と上記複数種の運転状態との差異に応じ上
   記複数のパラメータに重み付けをして上記実運転状態に
   対応した推定パラメータを合成し、この推定パラメータ
   により上記実運転状態に応じた形状の上記特性関数を決
   定した後、決定された上記特性関数に基づいて上記遠心
   圧縮機が所定の運転状態となるように上記の入口ガイド
   ベーンおよびディフューザベーンのベーン角度操作量を
   決定して上記遠心圧縮機を制御することを特徴とする遠
   心圧縮機の制御方法。