(57)【要約】【課題】圧縮機の性能を適切に予測する圧縮機の性能予測装置等を提供する。【解決手段】性能予測装置３は、圧縮機の実測データを取得する実測データ取得手段３１と、試験ガス物性値補正式が格納される試験ガス物性値補正式データベース３２と、圧縮機の試験パラメータを算出する試験パラメータ算出部３３ａと、対象となる試験ガスに含まれるガスの種類及び混合比に基づいて、試験ガス物性値補正式データベース３２から試験ガス物性値補正式を選択し、選択した前記試験ガス物性値補正式を用いて試験パラメータを補正する試験パラメータ補正部３３ｂと、を備える。【選択図】図２
2016205237.tif 000002
【特許請求の範囲】
【請求項１】
相似試験の試験対象である圧縮機によって、複数種類のガスを含む試験ガスを圧縮しているときに、前記圧縮機の流量、吸込温度、吸込圧力、吐出温度、及び吐出圧力の実測データを取得する実測データ取得手段と、
前記相似試験で実際に用いられる試験ガスの圧縮係数及び定容比熱を含む試験ガス物性値と、前記ガスの混合比が異なる複数の試験ガスの前記試験ガス物性値と、の関係を示す試験ガス物性値補正式が、前記ガスの種類及び混合比に対応付けて格納される試験ガス物性値補正式データベースと、
前記実測データ取得手段によって取得される前記実測データに基づいて、前記圧縮機の性能を示すポリトロープヘッド及びポリトロープ効率を含む試験パラメータを算出する試験パラメータ算出手段と、
前記圧縮機の性能を予測する際、対象となる試験ガスに含まれる前記ガスの種類及び混合比に基づいて、前記試験ガス物性値補正式データベースから前記試験ガス物性値補正式を選択し、選択した前記試験ガス物性値補正式を用いて前記試験パラメータを補正する試験パラメータ補正手段と、を備えること
を特徴とする圧縮機の性能予測装置。
【請求項２】
前記試験ガス物性値補正式データベースには、前記試験ガス物性値補正式が、前記ガスの種類、混合比、及び分子量に対応付けて格納され、
前記試験パラメータ補正手段は、前記試験ガス物性値補正式データベースに格納されている試験ガスの分子量と、前記対象となる試験ガスの分子量と、の大小関係に基づき、前記対象となる試験ガスに対応する試験ガス物性値補正式の係数を線形補間によって求めること
を特徴とする請求項１に記載の圧縮機の性能予測装置。
【請求項３】
前記相似試験の試験設備とは異なる現地で前記圧縮機が稼動する際の現地稼動条件と、前記試験パラメータ算出手段によって算出される前記試験パラメータと、に基づいて、前記現地における前記圧縮機の性能を示す吐出圧力及び動力を含む現地性能パラメータを算出する現地性能パラメータ算出手段を備えること
を特徴とする請求項１又は請求項２に記載の圧縮機の性能予測装置。
【請求項４】
前記試験パラメータ補正手段による補正後の前記試験パラメータが所定範囲に含まれ、かつ、前記現地性能パラメータ算出手段によって算出される前記現地性能パラメータが所定範囲に含まれる場合、前記圧縮機は性能に関する要件を満たしていると判定する合否判定手段を備えること
を特徴とする請求項３に記載の圧縮機の性能予測装置。
【請求項５】
複数種類のガスを含むとともに前記現地で前記圧縮機によって圧縮されることが想定される現地ガスの圧縮係数及び定容比熱を含む現地ガス物性値と、前記ガスの混合比が異なる複数の現地ガスの前記現地ガス物性値と、の関係を示す現地ガス物性値補正式が、前記ガスの種類及び混合比に対応付けて格納される現地ガス物性値補正式データベースと、
前記現地で前記圧縮機によって実際に圧縮される現地ガスに含まれる前記ガスの種類及び混合比に基づいて、前記現地ガス物性値補正式データベースから前記現地ガス物性値補正式を選択し、選択した前記現地ガス物性値補正式を用いて前記現地性能パラメータを補正する現地性能パラメータ補正手段と、を備えること
を特徴とする請求項３に記載の圧縮機の性能予測装置。
【請求項６】
前記現地ガス物性値補正式データベースには、前記現地ガス物性値補正式が、前記ガスの種類、混合比、及び分子量に対応付けて格納され、
前記現地性能パラメータ補正手段は、前記現地ガス物性値補正式データベースに格納されている現地ガスの分子量と、前記現地で前記圧縮機によって実際に圧縮される現地ガスの分子量と、の大小関係に基づき、当該現地ガスに対応する現地ガス物性値補正式の係数を線形補間によって求めること
を特徴とする請求項５に記載の圧縮機の性能予測装置。
【請求項７】
前記試験パラメータ補正手段による補正後の前記試験パラメータが所定範囲に含まれ、かつ、前記現地性能パラメータ補正手段による補正後の前記現地性能パラメータが所定範囲に含まれる場合、前記圧縮機は性能に関する要件を満たしていると判定する合否判定手段を備えること
を特徴とする請求項５に記載の圧縮機の性能予測装置。
【請求項８】
相似試験の試験対象である圧縮機の性能を示すポリトロープヘッド及びポリトロープ効率を含む試験パラメータを取得する試験パラメータ取得手段と、
前記相似試験の試験設備とは異なる現地で前記圧縮機が使用される際の現地稼動条件と、前記試験パラメータ取得手段によって取得される前記試験パラメータと、に基づいて、前記現地における前記圧縮機の性能を示す吐出圧力及び動力を含む現地性能パラメータを算出する現地性能パラメータ算出手段と、
複数種類のガスを含むとともに前記現地で前記圧縮機によって圧縮されることが想定される現地ガスの圧縮係数及び定容比熱を含む現地ガス物性値と、前記ガスの混合比が異なる複数の現地ガスの前記現地ガス物性値と、の関係を示す現地ガス物性値補正式が、前記ガスの種類及び混合比に対応付けて格納される現地ガス物性値補正式データベースと、
前記現地で前記圧縮機によって実際に圧縮される現地ガスに含まれる前記ガスの種類及び混合比に基づいて、前記現地ガス物性値補正式データベースから前記現地ガス物性値補正式を選択し、選択した前記現地ガス物性値補正式を用いて前記現地性能パラメータを補正する現地性能パラメータ補正手段と、を備えること
を特徴とする圧縮機の性能予測装置。
【請求項９】
前記現地ガス物性値補正式データベースには、前記現地ガス物性値補正式が、前記ガスの種類、混合比、及び分子量に対応付けて格納され、
前記現地性能パラメータ補正手段は、前記現地ガス物性値補正式データベースに格納されている現地ガスの分子量と、前記現地で前記圧縮機によって実際に圧縮される現地ガスの分子量と、の大小関係に基づき、当該現地ガスに対応する現地ガス物性値補正式の係数を線形補間によって求めること
を特徴とする請求項８に記載の圧縮機の性能予測装置。
【請求項１０】
前記現地性能パラメータ補正手段による補正後の前記現地性能パラメータが所定範囲に含まれる場合、前記圧縮機は性能に関する要件を満たしていると判定する合否判定手段を備えること
を特徴とする請求項８又は請求項９に記載の圧縮機の性能予測装置。
【請求項１１】
相似試験の試験対象である圧縮機によって、複数種類のガスを含む試験ガスを圧縮しているときに、前記圧縮機の流量、吸込温度、吸込圧力、吐出温度、及び吐出圧力の実測データを取得する実測データ取得ステップと、
前記実測データ取得ステップで取得される前記実測データに基づいて、前記圧縮機の性能を示すポリトロープヘッド及びポリトロープ効率を含む試験パラメータを算出する試験パラメータ算出ステップと、
前記相似試験で実際に用いられる試験ガスの圧縮係数及び定容比熱を含む試験ガス物性値と、前記ガスの混合比が異なる複数の試験ガスの前記試験ガス物性値と、の関係を示す試験ガス物性値補正式が前記ガスの種類及び混合比に対応付けて格納される試験ガス物性値補正式データベースから、前記圧縮機の性能を予測する際の対象となる試験ガスに含まれる前記ガスの種類及び混合比に基づいて前記試験ガス物性値補正式を選択し、選択した前記試験ガス物性値補正式を用いて前記試験パラメータを補正する試験パラメータ補正ステップと、を含むこと
を特徴とする圧縮機の性能予測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、圧縮機の性能予測装置及び性能予測方法に関する。
【背景技術】
【０００２】
圧縮機は、化学プラントや機械設備等において広く用いられている。圧縮機をユーザに提供する前には、例えば、米国機械学会（American Society of Mechanical Engineers：ＡＳＭＥ）によるパフォーマンス・テスト・コード１０（ＰＴＣ１０）に準拠した相似試験を行い、圧縮機が、ユーザから指定された満足すべき性能などの要件を満たしているか否かが試される。前記した「相似試験」とは、試験設備において圧縮機を実際に稼動し、その効率等が満足すべき範囲に入っているか否かを確認するための試験である。このような相似試験に関する技術として、例えば、以下に示すものが知られている。
【０００３】
すなわち、特許文献１には、「４０ｇ／ｇｍｏｌから１５０ｇ／ｇｍｏｌの間の分子量、７００未満の地球温暖化係数（ＧＷＰ）、及び１から１．５の間のガス比熱比を有する試験ガス」を用いて、圧縮機の相似試験を行うことが記載されている。なお、「試験ガス」とは、圧縮機の相似試験で用いられるガスである。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１２−１３７０８７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
特許文献１に記載の相似試験システムでは、ユーザから指定された現地（例えば、化学プラント）で圧縮機を運転する際に使用する現地ガス組成に対してＰＴＣ１０に基づいて圧縮機メーカが選択した試験ガスを用いて圧縮機を稼動し、圧縮機の吸込側の温度・圧力、吐出側の温度・圧力に基づいて試験パラメータが計算される。そして、前記した試験パラメータと、これに対応する指定パラメータと、を比較することで、圧縮機が相似試験に合格したか否かが判定される。
【０００６】
ところで、圧縮機の相似試験では、既存の計算手段を用いて試験ガスの物性値（例えば、圧縮係数）が計算されることが多いが、相似試験で使用する試験ガスの種類は多様であり、また、複数種類のガスを混合した試験ガスを用いる場合もある。したがって、既存の計算手段を用いて計算した試験ガス物性値が、相似試験での吸込温度、吸込圧力、吐出温度、及び吐出圧力の条件下において実測値と良好に一致するとは限らない。
仮に、試験ガス物性値の計算値と実測値との誤差が大きい場合には、圧縮機が相似試験に合格したか否かの判定に用いる試験パラメータの計算値が実際の値と良好に一致せず、圧縮機が現地に据付けられて稼動した際、満足すべき性能が得られない可能性がある。
【０００７】
そこで、本発明は、圧縮機の性能を適切に予測する性能予測装置等を提供することを課題とする。
【課題を解決するための手段】
【０００８】
前記課題を解決するために、本発明は、相似試験の試験対象である圧縮機によって、複数種類のガスを含む試験ガスを圧縮しているときに、前記圧縮機の流量、吸込温度、吸込圧力、吐出温度、及び吐出圧力の実測データを取得する実測データ取得手段と、前記相似試験で実際に用いられる試験ガスの圧縮係数及び定容比熱を含む試験ガス物性値と、前記ガスの混合比が異なる複数の試験ガスの前記試験ガス物性値と、の関係を示す試験ガス物性値補正式が、前記ガスの種類及び混合比に対応付けて格納される試験ガス物性値補正式データベースと、前記実測データ取得手段によって取得される前記実測データに基づいて、前記圧縮機の性能を示すポリトロープヘッド及びポリトロープ効率を含む試験パラメータを算出する試験パラメータ算出手段と、前記圧縮機の性能を予測する際、対象となる試験ガスに含まれる前記ガスの種類及び混合比に基づいて、前記試験ガス物性値補正式データベースから前記試験ガス物性値補正式を選択し、選択した前記試験ガス物性値補正式を用いて前記試験パラメータを補正する試験パラメータ補正手段と、を備えることを特徴とする。
【発明の効果】
【０００９】
発明によれば、圧縮機の性能を適切に予測する性能予測装置等を提供できる。
【図面の簡単な説明】
【００１０】
【図１】本発明の第１実施形態に係る性能予測装置によって性能が予測される圧縮機の試験設備の構成図である。
【図２】圧縮機の性能予測装置の機能ブロック図である。
【図３】（ａ）は相似試験での圧縮係数Ｚtの計算値と、相似試験での圧縮係数Ｚtの実測値との関係を示す説明図であり、（ｂ）は相似試験での定容比熱Ｃｖtの計算値と、相似試験での定容比熱Ｃｖtの実測値との関係を示す説明図である。
【図４】試験ガス物性値補正式データベースに格納される情報を示す説明図である。
【図５】性能予測装置が実行する処理を示すフローチャートである。
【図６】本発明の第２実施形態に係る性能予測装置の機能ブロック図である。
【図７】現地ガス物性値補正式データベースに格納される情報を示す説明図である。
【図８】性能予測装置が実行する処理を示すフローチャートである。
【図９】本発明の第３実施形態に係る性能予測装置の機能ブロック図である。
【図１０】性能予測装置が実行する処理を示すフローチャートである。
【発明を実施するための形態】
【００１１】
≪第１実施形態≫
以下では、圧縮機２（図１参照）の相似試験を行うための試験設備１（図１参照）について説明した後、本実施形態に係る性能予測装置３（図２参照）について詳細に説明する。
【００１２】
＜試験設備の構成＞
図１は、第１実施形態に係る性能予測装置３（図２参照）によって性能が予測される圧縮機２の試験設備１の構成図である。
圧縮機２は、例えば、一軸多段型の遠心圧縮機であり、図１に示す駆動軸２ａと、この駆動軸２ａと一体で回転するロータ（図示せず）と、ロータに固定される動翼（図示せず）と、ロータ及び動翼を収容するケーシング（図示せず）と、を備えている。そして、ケーシングと、回転する動翼と、の間を試験ガスが通流する過程で、動翼によって試験ガスにエネルギを与えて、試験ガスを昇圧するようになっている。
【００１３】
ここで、「試験ガス」とは、圧縮機２の相似試験で用いられるガスである。「試験ガス」には、相似試験で圧縮機２によって実際に圧縮するガスの他、後記する性能予測装置３（図２参照）によって圧縮機２の性能計算を行う際、圧縮機２によって圧縮すると仮定されるガスも含まれる。
また、「相似試験」とは、現地（例えば、化学プラント）で圧縮機２を実際に使用する前に、この圧縮機２が、ユーザから指定された満足すべき性能を有しているか否かを確認するための試験である。
【００１４】
ちなみに、圧縮機２が実際に使用される現地では、圧縮機２によって圧縮したガスが下流側の機器（図示せず）に供給されるが、相似試験の試験設備１では、圧縮したガスを自身の吸込側に戻すように圧縮機２が設置されている。
【００１５】
図１に示す試験設備１は、後記する性能予測装置３（図２参照）によって圧縮機２の現地稼動条件での性能を予測するために、試験ガスを使用して圧縮機２を実際に稼動させて、少なくとも圧縮機２の吸込側と吐出側における試験ガスの圧力・温度等を検出するための設備である。なお、「現地稼動条件」とは、圧縮機２が実際に現地（例えば、化学プラント）で使用される際の回転速度、吸込位置での温度・圧力・流量などの稼動条件である。
試験設備１は、ガス供給源１１と、ガス供給弁１２と、ガスパージ弁１３と、ガスリザーブ容器１４と、熱交換器１５と、吸込絞り弁１６と、モータ１７と、変速機１８と、流量センサ１９ａと、吸込温度センサ１９ｂと、吸込圧力センサ１９ｃと、吐出温度センサ１９ｄと、吐出圧力センサ１９ｅと、を備えている。図１に示すように、圧縮機２と、熱交換器１５と、吸込絞り弁１６と、流量センサ１９ａと、は環状に順次接続されている。
【００１６】
ガス供給源１１は、相似試験で用いる試験ガスの供給源であり、配管ｐ１，ｐ３（一部）を介して、圧縮機２の吸込側に接続されている。試験ガスとして、例えば、窒素、二酸化炭素、ヘリウム、フロン、メタン、エタン、プロパンのうち一つを用いてもよいし、これらのうち複数のガスを所定の比率で混合してもよい。ガス供給弁１２は、ガス供給源１１からのガスの供給／遮断を切り替える弁であり、配管ｐ１に設置されている。
【００１７】
ガスパージ弁１３は、圧縮機２において圧縮される試験ガスの濃度を調整する弁であり、配管ｐ２に設置されている。ガスリザーブ容器１４は、ガスパージ弁１３が開かれることで、配管ｐ１（一部），ｐ２を介して分流したガスを貯留する容器である。
熱交換器１５は、圧縮機２から吐出される高温のガスを、冷却水などの冷媒との熱交換によって冷却するものである。吸込絞り弁１６は、圧縮機２の吸込側に向かうガスの流量を調整する弁である。モータ１７は、圧縮機２に軸動力を付与する動力源である。変速機１８は、モータ１７の動力を所定の変速比で駆動軸２ａに伝達するものである。
【００１８】
流量センサ１９ａは、ノズル１９１ａにおけるガスの差圧に基づいて、ガスの流量を測定するセンサである。
吸込温度センサ１９ｂは、圧縮機２の吸込温度を検出するセンサである。吸込圧力センサ１９ｃは、圧縮機２の吸込圧力を検出するセンサである。吸込温度センサ１９ｂ及び吸込圧力センサ１９ｃは、圧縮機２の吸込口付近に設置されている。
吐出温度センサ１９ｄは、圧縮機２の吐出温度を検出するセンサである。吐出圧力センサ１９ｅは、圧縮機２の吐出圧力を検出するセンサである。吐出温度センサ１９ｄ及び吐出圧力センサ１９ｅは、圧縮機２の吐出口付近に設置されている。
【００１９】
流量センサ１９ａ、吸込温度センサ１９ｂ、吸込圧力センサ１９ｃ、吐出温度センサ１９ｄ、及び吐出圧力センサ１９ｅの検出値は、次に説明する性能予測装置３（図２参照）に出力される。
【００２０】
＜性能予測装置の構成＞
図２は、圧縮機２の性能予測装置３の機能ブロック図である。
性能予測装置３は、各センサ１９ａ〜１９ｅの検出値を入力値として使用した性能計算を行うことで、圧縮機２の性能を予測する装置である。性能予測装置３は、図示はしないが、ＣＰＵ（Central Processing Unit）、ＲＯＭ（Read Only Memory）、ＲＡＭ（Random Access Memory）、各種インタフェース等の電子回路を含んで構成される。そして、ＲＯＭに記憶されたプログラムをＲＡＭに展開し、ＣＰＵが処理を実行するようになっている。
【００２１】
図２に示すように、性能予測装置３は、実測データ取得手段３１と、試験ガス物性値補正式データベース３２と、演算処理手段３３と、表示制御手段３４と、を備えている。
実測データ取得手段３１は、流量センサ１９ａ、吸込温度センサ１９ｂ、吸込圧力センサ１９ｃ、吐出温度センサ１９ｄ、及び吐出圧力センサ１９ｅの検出値（実測データ）を、例えば、所定周期で取得する機能を有している。つまり、実測データ取得手段３１は、相似試験の試験対象である圧縮機２（図１参照）の稼動中、試験ガスを圧縮する圧縮機２の実測データを取得するようになっている。
【００２２】
試験ガス物性値補正式データベース３２には、相似試験で実際に用いられる試験ガスの試験ガス物性値（圧縮係数、定容比熱）と、相似試験で実際には用いられないものの、ガスの混合比が異なる複数の試験ガスの試験ガス物性値に関して、事前に別途実施するガス物性値測定実験で得られたそれぞれのガス物性値の実測値と、ガス物性値をガス混合比などから計算する既存の計算手段による計算値と、の関係を示す試験ガス物性値補正式が格納されている。
【００２３】
例えば、圧縮機２を使用する顧客から、所定の混合比で２種類のガスＧ１t，Ｇ２tを含む試験ガスを用いた場合について、圧縮機２の性能が知りたいという要望が寄せられることが多い。なお、圧縮機２の構成が同一であっても、試験ガスの組成（含まれるガスの種類・混合比）が異なれば、試験ガスの圧縮係数や定容比熱も異なる値になり、結果的に圧縮機２の効率等も異なる値になる。
【００２４】
顧客から指定されるたびに、その都度、試験ガスを作成して圧縮機２の相似試験を行うことも考えられるが、多大な時間・コストを要する。したがって、本実施形態では、将来、顧客から指定される可能性が高いガスの組合せ（例えば、ガスＧ１t，Ｇ２t）を想定し、その種類・混合比・分子量に対応付けて試験ガス物性値補正式をデータベース化しておくという方法をとっている。なお、ガスＧ１t，Ｇ２tの下付きの‘ｔ’は、（現地の仕様ではなく）圧縮機２の相似試験に関するものであることを示している。
【００２５】
以下では、２種類のガスＧ１t，Ｇ２tを混合した試験ガス（一方が０％、他方が１００％の場合を含む）を例に、試験ガス物性値補正式データベース３２に格納されている情報について説明する。
【００２６】
図３（ａ）は、相似試験での圧縮係数Ｚtの計算値と、相似試験での圧縮係数Ｚtの実測値との関係を示す説明図である。図３（ａ）の縦軸は、試験設備１（図１参照）とは別に、事前に別途実施するガス物性値測定実験装置（図示せず）において、２種類のガスＧ１t，Ｇ２tを所定の比率で混合してなる試験ガスをチャンバ（図示せず）に圧送し、その温度Ｔt_cor・圧力Ｐt_corを変化させたときの密度ρt_corの実測値を使用して、以下の（数式１）から算出される圧縮係数の実測値Ｚt_corである。なお、（数式１）に示すＲt［Ｊ／ｋｇ・Ｋ］は試験ガスの気体定数である。
【数１】
2016205237.tif 000003
【００２７】
例えば、事前に別途実施するガス物性値測定実験装置（図示せず）において、試験ガスＭｉｘ１t（ガスＧ１t：１００％、ガスＧ２t：０％）をチャンバに圧送し、その温度Ｔt_cor・圧力Ｐt_corを変化させることで、図３（ａ）に示す５つのｑ１の圧縮係数Ｚt_corが実測されたとする。なお、Ｚt_corの下付きの‘cor’は、「補正に使用する実測値」という意味である。
【００２８】
その後、相似試験で実際に使用する試験ガス（例えば、ガスＧ１t：３０％、ガスＧ２t：７０％）について、各センサ１９ｂ〜１９ｅ（図１参照）の検出値に対応する温度・圧力等が性能予測装置３に入力される。さらに、以下の（数式２）に基づき、その圧縮係数Ｚt［−］が、圧縮機２の相似試験で実測された温度・圧力に応じて算出される。この圧縮係数Ｚtの計算値が、図３（ａ）の横軸である。
なお、（数式２）に示すＰtは、相似試験実施中にセンサ１９ｃ，１９ｅ（図１参照）で検出された試験ガスの圧力であり、Ｔt［Ｋ］は相似試験実施中にセンサ１９ｂ，１９ｄ（図１参照）で検出された試験ガスの温度である。ここでは、センサ１９ｃ，１９ｅで検出された異なる圧力を平均した値をＰｔとし、かつ、センサ１９ｂ，１９ｄで検出された異なる温度を平均した値をＴｔとして使用する場合を例にとって説明する。なお、（数式２）に示すρt［ｋｇ／ｍ3］は試験ガスを対象に既存の計算手段で算出した密度であり、Ｒt［Ｊ／ｋｇ・Ｋ］は試験ガスの気体定数である。
【００２９】
【数２】
2016205237.tif 000004
【００３０】
性能予測装置３は、例えば、最小二乗法に基づいて５つの点ｑ１を直線近似し、図３（ａ）に示す直線Ａ１を表す関数を保有している。同様にして、性能予測装置３は、試験ガスＭｉｘ２t（ガスＧ１t：８０％、ガスＧ２t：２０％）、試験ガスＭｉｘ３t（ガスＧ１t：５０％、ガスＧ２t：５０％）、試験ガスＭｉｘ４t（ガスＧ１t：２０％、ガスＧ２t：８０％）、及び試験ガスＭｉｘ５t（ガスＧ１t：０％、ガスＧ２t：１００％）のそれぞれについて、圧縮係数の補正式を保有している。つまり、性能予測装置３は、図３（ａ）に示す直線Ａ２〜Ａ５を表す関数を保有している。これらの情報は、試験ガス物性値補正式データベース３２（図２参照）に格納される。
【００３１】
ちなみに、傾きが１の直線近似となり、かつ、軸切片が０であれば、圧縮係数Ｚt（計算値）と、圧縮係数Ｚt_cor（実測値）とが等しくなる（図３（ａ）の破線：直線Ｂを参照）。
【００３２】
図３（ｂ）は、相似試験での定容比熱Ｃｖtの計算値と、相似試験での定容比熱Ｃｖtの実測値との関係を示す説明図である。図３（ｂ）の縦軸は、前記した試験ガスＭｉｘ１t〜Ｍｉｘ５tのそれぞれについて、事前に別途実施するガス物性値測定実験装置（図示せず）において、チャンバ（図示せず）における温度・圧力を変化させたときの試験ガスの定容比熱Ｃｖt_cor（実測値）である。図３（ｂ）に示すように、チャンバ内の温度・圧力を変化させることで５つの点ｒ１において定容比熱Ｃｖt_corが検出されたとする。
図３（ｂ）の横軸は、相似試験で実際に使用する試験ガスについて、５つの点ｒ１のひとつひとつに対応する温度・圧力等に基づき、周知の方法で算出される定容比熱Ｃｖt（計算値）である。
【００３３】
性能予測装置３は、例えば、最小二乗法に基づいて５つの点ｒ１を直線近似し、図３（ｂ）に示す直線Ｃ１を表す関数を保有している。同様にして、性能予測装置３は、相似試験で実際に使用する試験ガスの定容比熱Ｃｖt（計算値）から、前記した試験ガスＭｉｘ２t〜Ｍｉｘ５tの定容比熱Ｃｖt_cor（実測値）を導くための関数を保有している。つまり、性能予測装置３は、図３（ｂ）に示す直線Ｃ２〜Ｃ５を表す関数を保有している。これらの情報は、圧縮機２の相似試験に先立って、試験ガス物性値補正式データベース３２（図２参照）に予め記憶されている。
【００３４】
ちなみに、傾きが１の直線近似となり、かつ、軸切片が０であれば、定容比熱Ｃｖt（計算値）と、定容比熱Ｃｖt_cor（実測値）とが等しくなる（図３（ｂ）の破線：直線Ｄを参照）。
【００３５】
図４は、試験ガス物性値補正式データベース３２に格納される情報を示す説明図である。図４に示すように、試験ガスＭｉｘ１t〜Ｍｉｘ５tの圧縮係数の補正式と、定容比熱の補正式と、がガスＧ１t，Ｇ２tの種類と、その混合比（モル分率）と、試験ガスの分子量と、に対応付けて、試験ガス物性値補正式データベース３２（図２参照）に記憶されている。
【００３６】
例えば、図４に示すガスＧ１t：１００％、ガスＧ２t：０％の混合比の試験ガスＭｉｘ１tについて、その圧縮係数の補正式は、関数：Ｚt_cor＝Ａｚ１t×Ｚt＋Ｂｚ１tであり、図３（ａ）に示す直線Ａ１に対応している。
また、例えば、図４に示すガスＧ３t：５０％、ガスＧ２t：５０％の混合比の試験ガスＭｉｘ３tについて、その定容比熱の補正式は、関数：Ｃｖt_cor＝Ａｃｖ３t×Ｃｖt＋Ｂｃｖ３tであり、図３（ｂ）に示す直線Ｃ３に対応している。
【００３７】
例えば、図４に示す試験ガスの分子量が小さくなるにつれて（例えば、Ｍｗ_Ｍｉｘ１_t＞Ｍｗ_Ｍｉｘ2_t＞…＞Ｍｗ_Ｍｉｘ5_tとすると）、図３（ａ）に示す直線Ａ１から直線Ａ５に向けて傾きが徐々に小さくなり、また、切片も徐々に小さくなる。つまり、試験ガスの分子量が連続的に変化すると、それに伴って、圧縮係数Ｚt，Ｚt_corの関係を与える直線の傾き・切片も連続的に変化する。このように試験ガスのガス組成成分が同じであれば、分子量の変化に対して傾きと切片が連続的に変化する。したがって、実際にこのデータベース作成時にガス物性値を計測していない混合比のガスが試験ガスとして使用された場合でも、後記するように、試験ガスの分子量に基づく線形補間によって、圧縮係数Ｚt，Ｚt_corの関係を与える直線の傾き及び切片を導くことができる。
なお、定容比熱についても同様のことがいえる（図３（ｂ）参照）。
【００３８】
また、試験ガスに含まれるガスの種類が異なっていれば（例えば、図示しないガスＧ３t，Ｇ４tを混合してなる試験ガスについては）、試験ガス物性値補正式データベース３２における別の記憶領域に各情報が格納される。
【００３９】
再び、図２に戻って説明を続ける。演算処理手段３３は、圧縮機２（図１参照）の性能を表す性能パラメータに関する演算処理を行うものであり、試験パラメータ算出部３３ａと、試験パラメータ補正部３３ｂと、現地性能パラメータ算出部３３ｃと、合否判定部３３ｄと、を備えている。
試験パラメータ算出部３３ａは、実測データ取得手段３１によって取得される実測データに基づいて、圧縮機２の試験パラメータを算出する機能を有している。ここで、「試験パラメータ」とは、圧縮機２の性能の評価基準となる状態量であり、本実施形態では相似試験における圧縮機２のポリトロープヘッドと、ポリトロープ効率とを示している。
【００４０】
前記した「ポリトロープヘッド」とは、圧縮機２における真の圧縮過程に代えて、周知のポリトロープ圧縮を仮定することで近似的に求められる圧力ヘッドである。また、「ポリトロープ効率」とは、ポリトロープ圧縮の仮定に基づく有効仕事に対して、実際に要する比仕事が占める割合である。
【００４１】
試験パラメータ補正部３３ｂは、顧客等から指定される試験ガスに含まれるガスＧ１t,Ｇ２tの種類・混合比（図２に示す「試験ガス情報」）と、試験ガス物性値補正式データベース３２に格納されている情報と、に基づいて、圧縮機２の試験パラメータを補正する機能を有している。
【００４２】
現地性能パラメータ算出部３３ｃは、相似試験の試験設備１とは異なる現地で圧縮機２が稼動する際の現地稼動条件と、試験パラメータ算出部３３ａによって算出される試験パラメータと、に基づいて、圧縮機２の現地性能パラメータを算出する機能を有している。ここで、「現地性能パラメータ」とは、圧縮機２の性能の評価基準となる状態量であり、本実施形態では現地での圧縮機２の吐出圧力と、圧縮機２の稼動に要する動力とを示している。
【００４３】
合否判定部３３ｄは、試験パラメータ補正部３３ｂによる補正後の試験パラメータと、現地性能パラメータ算出部３３ｃによって算出される現地性能パラメータと、に基づいて、圧縮機２が、性能に関する所定の要件を満たしているか否かを判定する機能を有している。
なお、試験パラメータ算出部３３ａ、試験パラメータ補正部３３ｂ、現地性能パラメータ算出部３３ｃ、及び合否判定部３３ｄの処理については後記する。
【００４４】
表示制御部３４は、演算処理手段３３の処理結果を画像として表示装置４（例えば、ディスプレイ）に表示させる機能を有している。
【００４５】
＜性能予測装置の動作＞
図５は、性能予測装置３が実行する処理を示すフローチャートである。
ステップＳ１０１において性能予測装置３は、試験設備１において圧縮機２が実際に稼動しているときに、実測データ取得手段３１によって、各センサ１９ａ〜１９ｅから実測データを取得する（実測データ取得ステップ）。
【００４６】
ステップＳ１０２において性能予測装置３は、ステップＳ１０１で取得した実測データに基づき、試験パラメータ算出部３３ａによって、圧縮機２の試験パラメータを算出する（試験パラメータ算出ステップ）。
まず、性能予測装置３は、以下に示す（数式３）を用いて、相似試験における圧縮機２のポリトロープヘッドＨpol_t［Ｊ／ｋｇ］を算出する。なお、（数式３）に示すｎt［−］は、相似試験における圧縮機２のポリトロープ指数であり、ｆt［−］は、相似試験における圧縮機２のポリトロープファクタである。
また、Ｐd_t［Ｐａ］は、吐出圧力センサ１９ｅ（図１参照）によって検出される吐出圧力であり、Ｐi_t［Ｐａ］は、吸込圧力センサ１９ｃによって検出される吸込圧力である。ｖd_t［ｍ3／ｋｇ］は吐出ガス比容積であり、ｖi_t［ｍ3／ｋｇ］は吸込ガス比容積である。吐出ガス比容積ｖd_t及び吸込ガス比容積ｖi_tは、各センサ１９ａ〜１９ｅ（図１参照）の検出値に基づき、ガス物性値計算ソフト等を用いて周知の方法で算出される。
【００４７】
【数３】
2016205237.tif 000005
【００４８】
なお、（数式３）に示すポリトロープ指数ｎtは、以下の（数式４）に基づいて算出される。
【００４９】
【数４】
2016205237.tif 000006
【００５０】
また、（数式３）に示すポリトロープファクタｆtは、以下の（数式５）に基づいて算出される。なお、（数式５）に示すｈd_t’［Ｊ／ｋｇ］は、圧縮機２において等エンタルピ変化を仮定した場合における吐出ガスのエンタルピであり、ｈi_t［Ｊ／ｋｇ］は、吸込ガスのエンタルピである。ｖd_t’［ｍ3／ｋｇ］は、等エンタルピ変化を仮定した場合における吐出ガス比容積である。
【００５１】
【数５】
2016205237.tif 000007
【００５２】
前記したように、圧縮機２の性能を予測する際、相似試験で実際に使用される試験ガス（例えば、ガスＧ１t：３０％、ガスＧ２t：７０％）と、対象となる試験ガス（例えば、図４に示す試験ガスＭｉｘ３）と、は異なっている場合。つまり、両者は、圧縮係数・定容比熱を含むガス物性値が異なっている。
したがって、（数式３）に基づいて算出されるポリトロープヘッドＨpolと、求めようとする真のポリトロープヘッドと、の間には誤差がある。この誤差をゼロに近づけるために、本実施形態では、試験ガス物性値補正式データベース３２に格納されている情報に基づいて、ポリトロープヘッドＨpolを含む試験パラメータを補正するようにしている。
【００５３】
図５のステップＳ１０３において性能予測装置３は、試験ガス物性値補正式データベース３２から試験ガス物性値補正式を選択する。例えば、顧客からの要望に基づく試験ガスは、ガスＧ１t，Ｇ２tが所定の比率で混合されたものであり、その分子量Ｍｗ_tが、図４に示す分子量Ｍｗ_Ｍｉｘ3_tに等しかったとする。この場合、性能予測装置３は、試験ガスＭｉｘ３に対応する圧縮係数の補正式（Ｚt_cor＝Ａｚ３t×Ｚt＋Ｂｚ３t：図４参照）を試験ガス物性値補正式データベース３２から取得する。
【００５４】
また、実際の相似試験で使用する試験ガスの分子量Ｍｗ_tが、試験ガス物性値補正式データベース３２に格納されている各分子量のいずれにも等しくないこともある。例えば、試験ガスの分子量Ｍｗ_tが、図４に示すＭｉｘ１tの分子量Ｍｗ_Ｍｉｘ1_tよりも大きく、かつ、Ｍｉｘ２tの分子量Ｍｗ_Ｍｉｘ2_t未満であったとする。この場合、性能予測装置３は、以下の（数式６）・（数式７）に基づいて圧縮係数の補正式の係数Ａｚt，Ｂｚtを求める。
【００５５】
【数６】
2016205237.tif 000008
【００５６】
【数７】
2016205237.tif 000009
【００５７】
このように性能予測装置３は、ステップＳ１０３において圧縮係数の補正式で表される直線の傾きＡｚtと切片Ｂｚtとを、それぞれ、各試験ガスの分子量に基づいて算出する。つまり、性能予測装置３は、試験ガス物性値補正式データベース３２に格納されている各試験ガスの分子量（Ｍｗ_Ｍｉｘ1_t，Ｍｗ_Ｍｉｘ2_t）と、対象となる試験ガスの分子量Ｍｗ_tと、の大小関係に基づき、圧縮係数の補正式の係数Ａｚt，Ｂｚtを線形補間（比例計算）によって求めるようにしている。これによって、試験ガス物性値補正式データベース３２に格納すべき補正式の個数（図４では、５個）が比較的少なくても、対象となる試験ガスの圧縮係数Ｚt_corを適切に算出できる。
【００５８】
同様にして性能予測装置３は、試験ガス物性値補正式データベース３２に格納されている各試験ガスの分子量と、対象となる試験ガスの分子量と、の大小関係に基づき、定容比熱の補正式の係数Ａｃｖt，Ｂｃｖtを線形補間によって求め、さらに、補正後の定容比熱Ｃｖt_corを算出する。
【００５９】
以下では、試験ガス物性値補正式データベース３２に格納されている情報を直接又は間接に用いて状態量を算出する場合、「補正計算に基づく…」と記載する。
【００６０】
図５のステップＳ１０４において性能予測装置３は、試験パラメータ補正部３３ｂによって、試験パラメータを補正する（試験パラメータ補正ステップ）。具体的に説明すると、性能予測装置３は、以下の（数式８）を用いて、補正計算に基づくポリトロープヘッドＨpol_t_cor［Ｊ／ｋｇ］を算出する。
なお、（数式８）に記載したＨpol_t［Ｊ／ｋｇ］は（数式３）に基づく補正前のポリトロープヘッドであり、κt［−］は試験ガスの比熱比である。Ｚt[−]は補正前の圧縮係数であり、Ｚt_cor［−］は補正後の圧縮係数である。Ｒt［Ｊ／ｋｇ・Ｋ］は試験ガスの気体定数であり、Ｔi_t［Ｋ］は吸込温度センサ１９ｂによって検出される吸込温度である。Ａｚt，Ｂｚtは、試験ガス物性値補正式データベース３２に格納されている情報と、前記した（数式６）・（数式７）と、に基づく圧縮係数の補正式の係数である。
【００６１】
【数８】
2016205237.tif 000010
【００６２】
（数式８）の最上段の右辺の分母・分子は、それぞれ、試験ガスを理想気体として扱った場合の圧縮係数（分母はＺt、分子はＺt_cot）を含む断熱ヘッドを乗算した形になっている。これによって、（数式８）の次段に示すように試験ガスを実在気体として扱って補正する場合に比べて式を簡単化でき、さらに最下段に示すように、試験ガス物性値補正式データベース３２に格納されている情報（補正式の係数や分子量）に基づいて、補正後のポリトロープヘッドＨpol_t_corを算出できる。
【００６３】
なお、（数式８）に示す試験ガスの比熱比κt［−］は、以下の（数式９）に基づいて求められる。ここで、Ｃｐt［Ｊ／ｋｇ・Ｋ］は、相似試験における試験ガスの定容比熱であり、Ｃｖt［Ｊ／ｋｇ・Ｋ］は、相似試験における試験ガスの定圧比熱である。
【００６４】
【数９】
2016205237.tif 000011
【００６５】
さらに、性能予測装置３は、以下の（数式１０）を用いて、補正計算に基づく圧縮機２の理論ヘッドＨth_t_corを算出する。なお、理論ヘッドとは、圧縮機２の有効仕事を示す圧力ヘッドである。
また、（数式１０）に示すＣｖt[−]は補正前の定容比熱であり、Ｃｖt_cor［−］は補正後の定容比熱である。Ｔd_t［Ｋ］は、吐出温度センサ１９ｄによって検出される吐出温度であり、Ａｃｖt，Ｂｃｖtは、試験ガス物性値補正式データベース３２の情報に基づく定容比熱の補正式の係数である。
【００６６】
【数１０】
2016205237.tif 000012
【００６７】
次に、性能予測装置３は、前記した（数式８）、（数式１０）の算出結果を、以下の（数式１１）に代入し、補正計算に基づくポリトロープ効率ηpol_t_corを求める。このようにして性能予測装置３は、図５のステップＳ１０４において、補正計算に基づくポリトロープヘッドＨpol_t_cor（数式８）と、ポリトロープ効率ηpol_t_cor（数式１１）と、を含む「試験パラメータ」を算出する。
【００６８】
【数１１】
2016205237.tif 000013
【００６９】
図５のステップＳ１０５において性能予測装置３は、現地性能パラメータ算出部３３ｃによって、現地性能パラメータ（現地での圧縮機２の吐出圧力、動力）を求める。
例えば、性能予測装置３は、以下に示す一連の収束計算によって、補正計算に基づく現地稼動条件での圧縮機２の吐出圧力Ｐd_sp［Ｐａ］を求める。なお、下付きのｓｐは、現地稼動条件に基づくことを示しており、現地稼動条件での吐出圧力Ｐd_sp［Ｐａ］は、通常、相似試験における吐出圧力センサ１９ｅ（図１参照）の検出値とは異なる値になる。
【００７０】
まず、性能予測装置３は、以下の（数式１２）に基づいて、現地稼動条件における圧縮機２の吐出側におけるエンタルピｈd_sp［Ｊ／ｋｇ］を求める。なお、ｈi_sp［Ｊ／ｋｇ］は、現地稼動条件における圧縮機２の吸入側のエンタルピである。Ｈpol_t［Ｊ／ｋｇ］は、相似試験における圧縮機２のポリトロープヘッドであり、前記した（数式３）に基づいて求められる。ηpol_t[−]は、相似試験における圧縮機２のポリトロープ効率であり、ポリトロープヘッドＨpol_tに基づいて周知の方法で求められる。
【００７１】
【数１２】
2016205237.tif 000014
【００７２】
次に、性能予測装置３は、（数式１２）で算出したエンタルピｈd_spを一定とする等エンタルピ条件において、所定の吐出圧力Ｐd_sp_as［Ｐａ］を仮定し、以下の（数式１３）を用いて、暫定のポリトロープヘッドＨpol_t_as［Ｊ／ｋｇ］を算出する。
なお、（数式１３）に示すｎsp［−］は、現地稼動条件におけるポリトロープ指数であり、（数式４）と同様の方法で算出される。ｆt［−］はポリトロープファクタであり、（数式５）に基づいて算出される。Ｐd_sp［Ｋ］、ｖd_sp［ｍ3／ｋｇ］は、現地稼動条件における圧縮機２の吐出圧力と、比容積であり、Ｐi_sp［Ｋ］、ｖi_sp［ｍ3／ｋｇ］は、現地稼動条件における圧縮機２の吸込圧力と、比容積である。
【００７３】
【数１３】
2016205237.tif 000015
【００７４】
性能予測装置３は、暫定のポリトロープヘッドＨpol_t_as［Ｊ／ｋｇ］が、相似試験に基づくポリトロープヘッドＨpol_t［Ｊ／ｋｇ］よりも小さければ、暫定の吐出圧力Ｐd_sp_as［Ｐａ］を前回よりも大きくし、エンタルピｈd_spを一定として現地稼動条件における吐出温度Ｔd_sp［Ｋ］、比容積ｖd_sp等を再び算出する。そして、性能予測装置３は、暫定のポリトロープヘッドＨpol_t_asが、相似試験におけるポリトロープヘッドＨpol_tに一致するまで、（数式１２）・（数式１３）に基づく計算を繰り返す。
このようにして、性能予測装置３は、補正計算に基づく現地稼動条件での圧縮機２の吐出圧力Ｐd_sp_corを算出する。
【００７５】
また、性能予測装置３は、補正計算に基づく現地稼動条件での圧縮機２の動力Ｐｗsp_cor［Ｗ］を求めるために、それに先立って以下の（数式１４）に基づき、補正計算に基づく現地稼動条件での圧縮機２の吸込質量流量Ｇi_sp_cor［ｋｇ／ｓ］を算出する。
なお、（数式１４）に示すＱi_sp［ｍ3／ｓ］は、ユーザから仕様として与えられる現地稼動条件での圧縮機２の吸込体積流量である。Ｚsp［−］は、現地稼動条件での圧縮係数の計算値であり、（数式２）と同様の方法で求められる。Ｔi_sp［Ｋ］は、現地稼動条件での圧縮機２の吸込温度であり、Ｐi_sp［Ｋ］は、現地稼動条件での圧縮機２の吸込圧力である。
【００７６】
【数１４】
2016205237.tif 000016
【００７７】
そして、性能予測装置３は、以下の（数式１５）を用いて、補正計算に基づく現地稼動条件での圧縮機２の動力Ｐｗsp_cor［Ｗ］を算出する。なお、（数式１５）に示すκsp［−］は、現地ガスの比熱比であり、（数式９）と同様の方法で求められる。Ｃｖsp［Ｊ／ｋｇ・Ｋ］は、試験ガスの定容比熱であり、ユーザから予め提示される現地ガスの組成に基づいて、仕様として与えられる現地稼動条件に応じて既存の計算手段で算出される。
【００７８】
【数１５】
2016205237.tif 000017
【００７９】
このようにして性能予測装置３は、図５のステップＳ１０５において、現地稼動条件での吐出圧力Ｐd_sp_cor（収束計算）と、動力Ｐｗsp_cor（数式１５）を含む「現地性能パラメータ」を算出する。
【００８０】
図５のステップＳ１０６において性能予測装置３は、合否判定部３３ｄによって、圧縮機２の性能に関する合否判定処理を行う。例えば、性能予測装置３は、補正計算に基づくポリトロープヘッドＨpol_t_corが所定の要求値の１００％以上かつ１０５％未満であり、かつ、現地稼動条件での動力Ｐｗsp_corが所定の要求値の１０７％以下である場合、圧縮機２は性能に関する要件を満たしていると判定する。
また、性能予測装置３は、補正計算に基づくポリトロープヘッドＨpol_t_corが前記した範囲外であるか、又は、現地稼動条件での動力Ｐｗsp_corが前記した範囲外である場合、圧縮機２は性能に関する要件を満たしていないと判定する。
【００８１】
なお、補正計算に基づくポリトロープ効率ηpol_t_corや現地稼動条件での吐出圧力Ｐd_sp_corを合否判定の基準に加えてもよい。
【００８２】
図５のステップＳ１０７において性能予測装置３は、例えば、（数式１）〜（数式１５）の算出結果と、ステップＳ１０６の合否判定処理の結果と、を表示制御手段３４によって表示装置４に表示させる。これによって性能予測装置３の管理者は、圧縮機２の性能に関する情報を把握し、合否の判定結果を考慮して所定の対応をとることができる。例えば、圧縮機２のポリトロープヘッドが不足している場合には、ガスが通流する圧縮機２の動翼やケーシング内部の表面（図示せず）を研磨したり、圧縮機２の稼動方法を変更したりして、圧縮機２の相似試験が再び行われる。
【００８３】
＜効果＞
本実施形態によれば、試験ガスの物性値に関する情報を試験ガス物性値補正式データベース３２に予め格納しておくことで、試験ガス物性値補正式に基づく圧縮係数Ｚt_cor・定容比熱Ｃｖt_corを用いて、試験パラメータを補正できる。したがって、顧客から試験ガスが指定されるたびに、その試験ガスを大量に作成して圧縮機２の相似試験を行う必要がなくなるとともに、圧縮機２の試験パラメータを精度よく算出できる。
【００８４】
また、本実施形態では、前記した（数式６）・（数式７）に基づく線形補間によって、ガス物性値補正式の係数を算出するようにしている。したがって、２種類のガスＧ１t，Ｇ２tの混合比に対応付けて、例えば、５つの試験ガス物性値補正式（図４参照）を用意しておけば、所望の試験ガスに関する各係数Ａｚt，Ｂｚt，Ａｃｖt，Ｂｃｖtを線形補間に基づいて算出できる。
【００８５】
≪第２実施形態≫
第２実施形態に係る性能予測装置３Ａ（図６参照）は、現地ガス物性値補正式データベース３５と、現地性能パラメータ補正部３３ｅと、を備える点が第１実施形態と異なっている。また、第２実施形態は、合否判定部３３ｆ（図６参照）の処理内容が、第１実施形態で説明した合否判定部３３ｄ（図２参照）とは異なっている。なお、その他の構成については第１実施形態（図２参照）と同様である。したがって、第１実施形態と異なる部分について説明し、重複する部分については説明を省略する。
【００８６】
＜性能予測装置の構成＞
図６は、第２実施形態に係る性能予測装置３Ａの機能ブロック図である。
図６に示すように、性能予測装置３Ａは、実測データ取得手段３１と、試験ガス物性値補正式データベース３２と、現地ガス物性値補正式データベース３５と、演算処理手段３３Ａと、表示制御手段３４と、を備えている。
【００８７】
現地ガス物性値補正式データベース３５には、現地で圧縮機２によって圧縮されることが想定される現地ガスの現地ガス物性値（圧縮係数、定容比熱）と、ガスの混合比が異なる複数の現地ガスの現地ガス物性値と、の関係を示す現地ガス物性値補正式が格納されている。
【００８８】
前記した「現地ガス」とは、試験設備１（図１参照）とは異なる現地（例えば、化学プラント）で、圧縮機２によって実際に圧縮されるガスである。なお、圧縮機２の納入先が把握できていれば、性能予測装置３Ａの管理者が、現地ガスに含まれる主要なガスを想定することは可能である。
本実施形態では、事前のガス物性値測定実験に基づいて現地ガスのガス物性値を予めデータベース化しておき、その後に顧客から通知される実際の現地ガスの組成に基づいて、現地性能パラメータ（吐出圧力、動力）を補正するようにしている。
【００８９】
図７は、現地ガス物性値補正式データベース３５に格納される情報を示す説明図である。本実施形態では、現地ガスとして、２種類のガスＧａsp，Ｇｂsp（一方が１００％、他方が０％の場合を含む）が混合されたものを用いる場合を例に説明する。
図７に示すように、現地ガスＭｉｘ１sp〜Ｍｉｘ５spの圧縮係数の補正式、及び定容比熱の補正式が、各現地ガスに含まれるガスＧａsp，Ｇｂspの種類・混合比、及び現地ガスの分子量に対応付けて、現地ガス物性値補正式データベース３５に記憶されている。
【００９０】
なお、圧縮係数の補正式の導出方法や、定容比熱の補正式の導出方法については、第１実施形態と同様である。例えば、Ｚsp_cor＝Ａｚ２sp×Ｚsp＋Ｂｚ２spという関数の各係数Ａｚ２sp，Ｂｚ２spは、ガスＧａspが８０％、ガスＧｂspが２０％の混合比の（想定される）現地ガスの圧縮係数（実測値）と、ガスＧａsp，Ｇｂspを所定の比率で含む現地ガスの圧縮係数（計算値）と、で特定される各点を、最小二乗法等に基づいて直線近似することで求められる。また、定容比熱Ｃｖspについても同様である。
【００９１】
ちなみに、大量の試験ガスを用いる相似試験とは異なり、図７に示す各補正式を導くための事前のガス物性値測定実験は、比較的少量の現地ガスがあれば足りる。したがって、組成の複雑な現地ガスを用いて相時試験を行うことは困難であっても、各補正式を導くための事前実験については比較的容易に行うことができる。
【００９２】
図６に示す演算処理手段３３Ａは、試験パラメータ算出部３３ａと、試験パラメータ補正部３３ｂと、現地性能パラメータ算出部３３ｃと、現地性能パラメータ補正部３３ｅと、合否判定部３３ｆと、を備えている。
【００９３】
現地性能パラメータ補正部３３ｅは、現地で圧縮機２によって実際に圧縮される現地ガスの組成（現地ガスに含まれるガスの種類・混合比）に基づいて、現地性能パラメータを補正する機能を有している。この現地性能パラメータは、第１実施形態で説明したように、圧縮機２の吐出圧力Ｐd_sp_corと、圧縮機２の稼動に要する動力Ｐｗsp_corとを示している。
【００９４】
合否判定部３３ｆは、試験パラメータ補正部３３ｂによる補正後の試験パラメータと、現地性能パラメータ補正部３３ｅによる補正後の現地性能パラメータと、に基づいて、圧縮機２が性能に関する所定の要件を満たしているか否かを判定する機能を有している。なお、現地性能パラメータ補正部３３ｅ及び合否判定部３３ｆが実行する処理については後記する。
【００９５】
＜性能予測装置の処理＞
図８は、性能予測装置３Ａが実行する処理を示すフローチャートである。
ステップＳ２０１〜Ｓ２０４については、第１実施形態で説明したステップＳ１０１〜Ｓ１０４（図５参照）と同様であるから、説明を省略する。なお、現地では圧縮機２の上流側・下流側に所定の機器（図示せず）が設置されており、ＰＴＣ１０に準拠した吸込側の温度・圧力計測位置までの配管の直線部長さなどの試験条件を実現できないため（吐出側についても同様）、現地稼動条件で圧縮機２のポリトロープヘッドＨpol_t_cor（数式８）及びポリトロープ効率ηpol_t_cor（数式１１）でユーザ要求仕様に対して性能の合否を正当に判定することは困難である。したがって、本実施形態でも、試験ガスを用いた相似試験の結果に基づいて、第１実施形態と同様の方法でポリトロープヘッドＨpol_t_cor及びポリトロープ効率ηpol_t_corが算出される。
【００９６】
図８のステップＳ２０５において性能予測装置３Ａは、現地性能パラメータ算出部３３ｃによって、圧縮機２の現地性能パラメータを算出する。まず、性能予測装置３Ａは、以下の（数式１６）に基づき、現地稼動条件におけるポリトロープ指数ｎspを算出する。なお、ηpol_t［−］は、相似試験におけるポリトロープ効率であり、κsp［−］は、現地ガスの比熱比である。
【００９７】
【数１６】
2016205237.tif 000018
【００９８】
（数式１６）に示すポリトロープ効率ηpol_tは、例えば、以下の（数式１７）に基づいて算出される。なお、（数式１７）に示すポリトロープ指数ｎtは、第１実施形態で説明した（数式４）に基づいて算出され、ポリトロープファクタｆtは（数式５）に基づいて算出される。
また、（数式１７）に示す吐出ガスエンタルピｈd_t［Ｊ／ｋｇ］、吸込ガスエンタルピｈi_t［Ｊ／ｋｇ］、吐出圧力Ｐd_t［ｍ3／ｋｇ］は、吸込圧力Ｐi_t［ｍ3／ｋｇ］、吐出ガス比容積ｖd_t［ｍ3／ｋｇ］、吸込ガス比容積ｖi_t［ｍ3／ｋｇ］についても、相似試験の結果に基づいて周知の方法で求めれられる。
【００９９】
【数１７】
2016205237.tif 000019
【０１００】
また、性能予測装置３Ａは、以下の（数式１８）に基づいて、現地ガスの比熱比κspを算出する。なお、（数式１８）に示すＣｐsp［Ｊ／ｋｇ・Ｋ］は、現地ガスの定圧比熱であり、Ｃｖsp［Ｊ／ｋｇ・Ｋ］は、現地ガスの定容比熱である。
【０１０１】
【数１８】
2016205237.tif 000020
【０１０２】
また、性能予測装置３Ａは、以下の（数式１９）を用いて、補正計算に基づくポリトロープ指数ｎsp_cor［−］を算出する。なお、（数式１９）に示すηpol_t_cor［−］は、第１実施形態の（数式１１）で説明したように、補正計算に基づくポリトロープ効率である。
【０１０３】
【数１９】
2016205237.tif 000021
【０１０４】
次に、性能予測装置３Ａは、以下の（数式２０）を用いて、現地ガスを理想気体として扱った場合における、現地稼動条件での圧縮機２の吐出圧力Ｐd_sp_id［Ｐａ］を求める。
なお、（数式２０）に示すＰi_sp［Ｐａ］は、現地の仕様に基づいて与えられる圧縮機２の吸込圧力である。ポリトロープヘッドＨpol_t［Ｊ／ｋｇ］は、第１実施形態の（数式３）で説明したように、相似試験における圧縮機２のポリトロープヘッドである。
また、（数式２０）に示すＺsp［−］は、想定される現地ガスの圧縮係数であり、第１実施形態で説明した圧縮係数Ｚt（計算値）と同様に、事前に計算で求められる。また、Ｒsp［Ｊ／ｋｇ・Ｋ］は現地ガスの気体定数であり、Ｔi_sp［Ｋ］は現地の仕様に基づいて与えられる圧縮機２の吸込温度である。
【０１０５】
【数２０】
2016205237.tif 000022
【０１０６】
（数式２０）に示す圧縮係数Ｚsp［−］は、例えば、２種類のガスＧａsp、Ｇｂspが現地ガスに含まれるとの事前の想定に基づき、このガスＧａsp、Ｇｂspを所定の比率で混合してなる現地ガスについて、周知の方法で求められる圧縮係数である。つまり、事前の想定に基づいたガスＧａsp、Ｇｂspの所定の比率と、顧客から通知される現地ガスの組成とが異なる場合、（数式２０）に基づいて求められる吐出圧力Ｐd_sp_id［Ｐａ］と、実際に圧縮機２によって現地ガスを圧縮する際の吐出圧力と、の間には誤差がある。
したがって、本実施形態では、顧客から通知される現地ガスの組成（図６に示す「現地ガス情報」）と、現地ガス物性値補正式データベース３５に格納されている情報と、に基づいて、現地性能パラメータ（吐出圧力、動力）を補正するようにしている。
なお、顧客から通知される現地ガスの組成（現地ガスに含まれるガスの種類・混合比）は、管理者によって性能予測装置３Ａに入力される。
【０１０７】
図８のステップＳ２０６において性能予測装置３Ａは、現地性能パラメータ補正部３３ｅによって、現地ガス物性値補正式データベース３５から現地ガス物性値補正式を選択する。
例えば、顧客から通知される現地ガスの分子量Ｍｗ_spが、現地ガス物性値補正式データベース３５に格納されている分子量Ｍｗ_Ｍｉｘ3_spに等しければ、性能予測装置３Ａは、現地ガスＭｉｘ３spに対応する圧縮係数の補正式（Ｚsp_cor＝Ａｚ３sp×Ｚsp＋Ｂｚ３sp）を取得する。
【０１０８】
また、顧客からの通知される現地ガスの分子量Ｍｗ_spが、現地ガス物性値補正式データベース３５に格納されている各分子量のいずれにも等しくないこともある。例えば、試験ガス分子量Ｍｗ_spが、図７に示すＭｉｘ１spの分子量Ｍｗ_Ｍｉｘ1_spよりも大きく、かつ、Ｍｉｘ２spの分子量Ｍｗ_Ｍｉｘ2_sp未満であったとする。この場合、性能予測装置３Ａは、以下の（数式２１）・（数式２２）に基づいて圧縮係数の補正式の係数Ａｚsp，Ｂｚspを求める。
【０１０９】
【数２１】
2016205237.tif 000023
【数２２】
2016205237.tif 000024
【０１１０】
このように、性能予測装置３Ａは、現地ガス物性値補正式データベース３５に格納されている現地ガスの分子量と、現地で圧縮機２によって実際に圧縮される現地ガスの分子量と、の大小関係に基づき、圧縮係数の補正式の係数Ａzsp，Ｂzspを線形補間で求めるようにしている。同様に、定容比熱の補正式の係数Ａｃｖsp，Ｂｃｖspも線形補間で求められる。
【０１１１】
図８のステップＳ２０７において性能予測装置３Ａは、現地性能パラメータ補正部３３ｅによって、現地性能パラメータ（吐出圧力、動力）を補正する。まず、性能予測装置３Ａは、以下に示す（数式２３）に基づき、現地ガスを理想気体として扱った場合における、現地稼動条件での圧縮機２の吐出圧力Ｐd_sp_id_cor［Ｐａ］を算出する。
なお、（数式２３）に示すＰi_sp［Ｐａ］は、現地の仕様に基づいて与えられる吸込圧力である。Ｈpol_t_cor[Ｊ／ｋｇ]は、補正計算に基づく相似試験でのポリトロープヘッドであり、前記した（数式８）に基づいて求められる。ｎsp_cor［−］は、補正計算に基づく現地稼動条件でのポリトロープ指数であり、前記した（数式１９）に基づいて求められる。
また、（数式２３）に示すＺsp_cop［−]は、現地ガスの圧縮係数の補正値であり、Ｚsp［−]は、現地ガスに含まれるとの事前の想定に基づき、このガスＧａsp，Ｇｂspを所定の比率で混合してなる現地ガスについて、周知の方法で求められる圧縮係数である。また、Ａｚsp，Ｂｚspは、現地ガス物性値補正式データベース３５の情報に基づく圧縮係数の補正式の係数である。
【０１１２】
【数２３】
2016205237.tif 000025
【０１１３】
そして、性能予測装置３Ａは、以下の（数式２４）に基づき、補正計算に基づく現地稼動条件での圧縮機２の吐出圧力Ｐd_sp_cor［Ｐａ］を求める。なお、（数式２４）に示すＰd_sp［Ｐａ］は、現地稼動条件における圧縮機２の吸込圧力であり、第１実施形態と同様に、所定の収束計算（（数式１２）、（数式１３））に基づいて求められる。
【０１１４】
【数２４】
2016205237.tif 000026
【０１１５】
また、性能予測装置３Ａは、補正計算に基づく現地稼動条件での圧縮機２の動力Ｐｗsp_cor［Ｗ］を求めるために、それに先立って、以下の（数式２５）を用いて、現地稼動条件における圧縮機２の吸込質量流量Ｇi_sp_cor［ｋｇ／ｓ］を算出する。
なお、（数式２５）に示すＱi_sp［ｍ3／ｓ］は、仕様として与えられる現地稼動条件での圧縮機２の吸込体積流量である。圧縮係数Ｚsp［−］と、その補正式の係数Ａｚsp，Ｂｚspと、現地ガスの気体定数Ｒsp［Ｊ／Ｋ・ｋｇ］と、現地稼動条件での圧縮機２の吸込温度Ｔi_sp［Ｋ］と、吸込圧力Ｐi_sp［Ｐａ］と、については、前記したとおりである。
【０１１６】
【数２５】
2016205237.tif 000027
【０１１７】
さらに、性能予測装置３Ａは、以下の（数式２６）を用いて、補正計算に基づく現地稼動条件での圧縮機２の動力Ｐｗsp_cor［Ｗ］を算出する。なお、（数式２６）に示すκsp［−］は、現地ガスの比熱比であり、（数式９）と同様の方法で算出される。Ｃｖsp_cor［Ｊ／ｋｇ・Ｋ］は、現地稼動条件における現地ガスの定容比熱の補正値であり、現地ガス物性値補正式データベース３５の情報に基づいて求められる。係数Ａｃｖsp、Ｂｃｖspについては、前記した圧縮係数に関する係数Ａzsp、Ｂzsp（（数式２１）、（数式２２））と同様の方法で、線形補間によって求めることができる。
【０１１８】
【数２６】
2016205237.tif 000028
【０１１９】
このようにして性能予測装置３Ａは、図８のステップＳ２０７において、補正計算に基づく現地稼動条件での圧縮機２の吐出圧力Ｐd_sp_corと、動力Ｐｗsp_corと、を含む現地性能パラメータを算出する。
図８のステップＳ２０８において性能予測装置３Ａは、合否判定部３３ｆによって、圧縮機２の性能に関する合否判定処理を行う。つまり、性能予測装置３Ａは、補正後の試験パラメータ（ポリトロープヘッドＨpol_t_cor、ポリトロープ効率ηpol_t_cor）が所定範囲に含まれ、かつ、補正後の現地性能パラメータ（吐出圧力Ｐd_sp_cor、動力Ｐｗsp_cor）が所定範囲に含まれる場合、圧縮機２は性能に関する要件を満たしていると判定する。
図８のステップＳ２０９において性能予測装置３Ａは、例えば、ステップＳ２０１〜Ｓ２０８の一連の処理結果を、表示制御手段３４によって表示装置４に表示させる。
【０１２０】
＜効果＞
本実施形態によれば、現地ガスの物性値に関する情報を現地ガス物性値補正式データベース３５に予め格納しておくことで、現地ガス物性値補正式に基づく補正後の圧縮係数Ｚsp_cor・定容比熱Ｃｖsp_corを用いて、現地性能パラメータの補正値を算出できる。したがって、第１実施形態よりも圧縮機２の性能に関する合否を適切に判定できる。
【０１２１】
≪第３実施形態≫
第３実施形態は、事前に相似試験に相当する何らかの性能を評価する試験を実施した結果から試験パラメータが得られているときに実施する。すなわち、第３実施形態に係る性能予測装置３Ｂ（図９参照）は、第２実施形態で説明した実測データ取得手段３１（図６参照）に代えて、試験パラメータ取得手段３６（図９参照）を備え、さらに、第２実施形態で説明した構成（図６参照）から、試験ガス物性値補正式データベース３２、試験パラメータ算出部３３ａ、及び試験パラメータ補正部３３ｂを省略した構成になっている。なお、その他の各構成ついては、第２実施形態と同様である。したがって、第２実施形態と異なる部分について説明し、重複する部分については説明を省略する。
【０１２２】
＜性能予測装置の構成＞
図９は、第３実施形態に係る性能予測装置３Ｂの機能ブロック図である。
図９に示すように、性能予測装置３Ｂは、試験パラメータ取得手段３６と、現地ガス物性値補正式データベース３５と、演算処理手段３３Ｂと、表示制御手段３４と、を備えている。
【０１２３】
試験パラメータ取得手段３６は、圧縮機２のポリトロープヘッド及びポリトロープ効率を含む試験パラメータを取得する機能を有している。例えば、別のコンピュータ（図示せず）において、補正計算に基づくポリトロープヘッドＨpol_t_cor（数式８）及びポリトロープ効率ηpol_t_cor（数式１１）を含む試験パラメータを、相似試験の結果に基づいて算出し、このコンピュータから性能予測装置３Ｂに試験パラメータが入力するようにしてもよい。また、例えば、管理者によるキーボード（図示せず）の操作によって、性能予測装置３に試験パラメータの数値を入力するようにしてもよい。
【０１２４】
演算処理手段３３Ｂは、現地性能パラメータ算出部３３ｃと、現地性能パラメータ補正部３３ｅと、合否判定部３３ｆと、を備えている。
現地性能パラメータ算出部３３ｃは、現地で圧縮機２が使用される際の現地稼動条件と、試験パラメータ取得手段３６によって取得される試験パラメータと、に基づいて、圧縮機２の現地性能パラメータ（吐出圧力、動力）を算出する機能を有している。
現地性能パラメータ算出部３３ｃ及び合否判定部３３ｆについては、第２実施形態と同様であるから、説明を省略する。
【０１２５】
＜性能予測装置の処理＞
図１０は、性能予測装置３Ｂが実行する処理を示すフローチャートである。
ステップＳ３０１において性能予測装置３Ｂは、試験パラメータ取得手段３６によって、ポリトロープヘッド及びポリトロープ効率を含む試験パラメータを取得する。前記したように、別のコンピュータ（図示せず）から試験パラメータを取得するようにしてもよいし、管理者の操作によって試験パラメータを入力するようにしてもよい。
【０１２６】
ステップＳ３０２において性能予測装置３Ｂは、現地性能パラメータ算出部３３ｃによって、補正前の吐出圧力を含む現地性能パラメータを算出する。つまり、性能予測装置３Ｂは、ステップＳ３０１で取得した試験パラメータと、管理者によって入力される圧縮機２の現地稼動条件と、に基づいて、圧縮機２の現地性能パラメータを算出する。なお、ステップＳ３０２の処理は、第２実施形態で説明したステップＳ２０５（図８参照）と同様であるから、詳細な説明を省略する。
【０１２７】
また、ステップＳ３０３，Ｓ３０４の処理は、第２実施形態で説明したステップＳ２０６，Ｓ２０７（図８参照）と同様であるから、説明を省略する。
ステップＳ３０５において性能予測装置３Ｂは、ステップＳ３０４で求めた補正後の現地性能パラメータが所定範囲に含まれる場合、合否判定部３３ｆによって、圧縮機２は性能に関する要件を満たしていると判定する。
ステップＳ３０６において性能予測装置３Ｂは、例えば、ステップＳ３０１〜Ｓ３０５の一連の処理結果を、表示制御手段３４によって表示装置４に表示させる。
【０１２８】
＜効果＞
本実施形態によれば、試験パラメータ取得手段３６によって取得される試験パラメータ等に基づいて圧縮機２の現地性能パラメータを算出し、さらに、現地ガス物性値補正式データベース３５に格納されている情報に基づいて現地性能パラメータを補正できる。したがって、圧縮機２の性能に関する合否を容易かつ適切に判定できる。
【０１２９】
≪変形例≫
以上、本発明に係る性能予測装置３，３Ａ，３Ｂについて説明したが、本発明はこれらの記載に限定されるものではなく、種々の変更を行うことができる。
例えば、第１実施形態では、ガスＧ１t，Ｇ２tの混合比が異なる５つの試験ガスＭｉｘ１t〜Ｍｉ５tに対応して、試験ガス物性値補正式データベース３２（図４参照）に各情報が格納される場合について説明したが、これに限らない。すなわち、ガスＧ１t，Ｇ２tの混合比が異なる試験ガスは、４つ以下であってもよいし、また、６つ以上であってもよい。第２、第３実施形態で説明した現地ガス物性値補正式データベース３５（図７、図９参照）についても同様のことがいえる。
【０１３０】
また、各実施形態では、試験ガスに２種類のガスＧ１t，Ｇ２tが含まれる場合について説明したが、試験ガスに含まれるガスの種類は３種類以上であってもよい。この場合でも、（数式６）・（数式７）と同様にして、線形補間によって試験ガス物性値を補正できる。
なお、現地ガスについても同様である。
【０１３１】
また、第１実施形態では、対象とする試験ガスの分子量Ｍｗが、試験ガス物性値補正式データベース３２に格納されている各分子量のいずれにも等しくない場合、（数式６）・（数式７）を用いた線形補間によって係数Ａｚt，Ｂｚtを求める構成について説明したが、これに限らない。すなわち、試験ガス物性値補正式データベース３２に格納されている試験ガスＭｉｘ１t〜Ｍｉ５tのうち、その分子量が、対象とする試験ガスの分子量に最も近いものを選択し、選択した試験ガスに対応する補正式に基づいて圧縮係数Ｚtを算出するようにしてもよい。なお、試験ガスの定容比熱Ｃｖt、現地ガスの圧縮係数Ｚsp、及び、現地ガスの定容比熱Ｃｖspについても同様のことがいえる。
【０１３２】
また、第１実施形態では、試験パラメータ補正部３３ｂや現地性能パラメータ算出部３３ｃの処理結果に基づき、合否判定部３３ｄによって圧縮機２の性能に関する合否を判定する場合について説明したが、これに限らない。すなわち、合否判定部３３ｄを省略し、試験パラメータ補正部３３ｂや現地性能パラメータ算出部３３ｃの処理結果を表示装置４に表示させるようにしてもよい。この場合でも、性能予測装置３の管理者は、試験パラメータ・現地性能パラメータの数値を把握し、その数値に基づいてその後の対応を考えることができる。なお、第２、第３実施形態についても同様のことがいえる。
また、第１実施形態から合否判定部３３ｄ及び現地性能パラメータ算出部３３ｃを省略し、試験パラメータ補正部３３ｂの処理結果を表示装置４に表示するようにしてもよい。
【０１３３】
また、第２、第３実施形態では、現地性能パラメータ算出部３３ｃによって、現地での圧縮機２の吐出圧力（補正前）を算出する場合について説明したが、これに限らない。例えば、現地性能パラメータ算出部３３ｃによって、圧縮機２の吐出圧力（補正前）と、動力（補正前）と、の両方を算出するようにしてもよい。
【０１３４】
また、各実施形態では、「試験ガス物性値」として、試験ガスの圧縮係数及び定容比熱を用いる場合について説明したが、例えば、試験ガスのマッハ数を「試験ガス物性値」に含めてもよい（現地ガスについても同様）。
また、各実施形態では、「試験パラメータ」として、圧縮機２のポリトロープヘッド及びポリトロープ効率を算出する場合について説明したが、例えば、圧縮機２の理論ヘッドを「試験パラメータ」に含めてもよい。
また、各実施形態では、「現地性能パラメータ」として、現地での圧縮機２の吐出圧力及び動力を算出する場合について説明したが、例えば、圧縮機２の周速マッハ数（圧縮機２の回転速度／マッハ数）を「現地性能パラメータ」に含めてもよい。
【０１３５】
また、各実施形態では、試験ガス物性値補正式として、直線を表す一次関数を用いる場合について説明したが（図３（ａ）、（ｂ）参照）、曲線を表す所定の関数であってもよい。なお、現地ガス物性値補正式についても、同様のことがいえる。
また、各実施形態では、圧縮機２が一軸多段型の遠心圧縮機である場合について説明したが、これに限らない。すなわち、圧縮機２は、斜流圧縮機や軸流圧縮機であってもよいし、また、単段型の圧縮機であってもよい。
【０１３６】
また、各実施形態は本発明を分かりやすく説明するために詳細に記載したものであり、必ずしも説明した全ての構成を備えるものに限定されない。
また、前記した各構成、機能、処理部、処理手段等は、それらの一部又は全部を、例えば集積回路で設計する等によりハードウェアで実現してもよい。また、機構や構成は説明上必要と考えられるものを示しており、製品上必ずしも全ての機構や構成を示しているとは限らない。
【符号の説明】
【０１３７】
１  試験設備
２  圧縮機
３，３Ａ，３Ｂ  性能予測装置
３１  実測データ取得手段
３２  試験ガス物性値補正式データベース
３３，３３Ａ，３３Ｂ  演算処理手段
３３ａ  試験パラメータ算出部（試験パラメータ算出手段）
３３ｂ  試験パラメータ補正部（試験パラメータ補正手段）
３３ｃ  現地性能パラメータ算出部（現地性能パラメータ算出手段）
３３ｄ，３３ｆ  合否判定部（合否判定手段）
３３ｅ  現地性能パラメータ補正部（現地性能パラメータ補正手段）
３４  表示制御手段
３５  現地ガス物性値補正式データベース
３６  試験パラメータ取得手段
４  表示装置
【図１】
2016205237.tif 000029
【図２】
2016205237.tif 000030
【図３】
2016205237.tif 000031
【図４】
2016205237.tif 000032
【図５】
2016205237.tif 000033
【図６】
2016205237.tif 000034
【図７】
2016205237.tif 000035
【図８】
2016205237.tif 000036
【図９】
2016205237.tif 000037
【図１０】
2016205237.tif 000038