JPH06348686A - 配管系における圧力脈動の低減のためのオリフィス選定方法 - Google Patents

Abstract

(57)【要約】 【目的】往復圧縮機などの加圧源の圧力変動による配管
系の圧力脈動を低減するためのオリフィス選定方法に関
し、ニューラルネットワークのモデルの一つであるボル
ツマンマシンを用い、配管系におけるオリフィスの種類
及び位置の最適の組み合わせを自動的に選定することを
目的とする。 【構成】オリフィスの種類とオリフィスを挿入可能な位
置とを示すユニットマトリックスＵＭを定義し、ユニッ
トマトリックスＵＭの各行及び各列に１であるユニット
Ｕｉｋが１つずつあるときに最小となる第１の関数と、
脈動圧力値と基準値との差を評価する第２の関数との加
重加算によってネットワークのエネルギー関数を定義
し、模擬焼き鈍しの方法によって収束したときのユニッ
トＵｉｋの状態によって、オリフィスの種類と位置とを
決定する。

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は、往復圧縮機などの加圧
源の圧力変動による配管系の圧力脈動を低減するための
オリフィス選定方法に関する。

【０００２】

【従来の技術】往復圧縮機から吐出されるガスは間欠的
であるため、往復圧縮機に接続された配管内又は機器内
に圧力脈動が発生する。圧力脈動が大きい場合には、圧
縮機の性能低下をもたらしたり配管を破損させたりする
おそれがあるので、圧力脈動を所定の基準値以下に抑え
るように対策を講じる必要がある。そのような基準値を
定めた規格として、例えば米国石油協会規格（ＡＰＩ規
格）がある。

【０００３】配管系の圧力脈動を低減する方法として、
適当な種類（縮径比）のオリフィスを配管系の適所に挿
入することが知られている。オリフィスを挿入すると圧
力脈動は減少するが、しかし、オリフィスを挿入するこ
とによりエネルギーロスが生じて効率が低下するので、
挿入するオリフィスの個数はなるべく少なく且つその縮
径比はなるべく１に近いものが好ましい。

【０００４】したがって、圧力脈動を基準値以下に抑え
且つ効率の低下を最小限に抑えるために、基準値を超え
ずにできるだけ基準値に近い圧力脈動となるようなオリ
フィスの種類と位置との最適な組み合わせを選定する必
要がある。

【０００５】従来において、オリフィスの種類と位置を
選定するために、配管系における圧力脈動をシミュレー
ションする圧力脈動解析装置（又は圧力脈動解析プログ
ラム「ＰＵＬＳＡＳ」）が用いられている（特開平４−
６２４４６号公報、「往復圧縮機配管系の圧力脈動解
析」神戸製鋼技報Ｖｏｌ．３７Ｎｏ．１、「往復圧縮機
の圧力脈動解析」同Ｖｏｌ．４１Ｎｏ．１）。

【０００６】つまり、まず、圧力脈動解析装置を用いて
配管系の圧力脈動解析を行うことによって、配管系の位
置及び周波数に対する脈動圧力値を求める。基準値を超
える脈動圧力値がある場合に、操作者の勘と経験に基づ
いて適当な位置にオリフィスを模擬的に挿入し、又は位
置を変更し、再度圧力脈動解析を行う。従来において
は、これらの操作を試行錯誤で繰り返すことによって、
オリフィスの種類と位置とについての最適と思われる組
み合わせを見つけ出していた。

【０００７】

【発明が解決しようとする課題】しかし、上述した従来
の選定方法では、どんな種類のオリフィスを配管系のど
こに配置するか、という重要な点について、操作者の勘
と経験に頼っているため、操作者に対して高度な熟練と
工学的センスが要求され、且つ多くの時間を要している
とともに、操作者の個人差によっても結果が異なってく
るという問題があった。

【０００８】本発明は、ニューラルネットワークのモデ
ルの一つであるボルツマンマシンを用い、配管系におけ
るオリフィスの種類及び位置の最適の組み合わせを自動
的に選定するための方法を提供することを目的とする。

【０００９】

【課題を解決するための手段】請求項１の発明に係る方
法は、上述の課題を解決するため、配管系における圧力
脈動の低減のために前記配管系内に挿入するオリフィス
の種類と位置とをボルツマンマシンを用いて選定する方
法であって、オリフィスの種類とオリフィスを挿入可能
な位置とを示すユニットマトリックスを定義し、前記ユ
ニットマトリックスの各行及び各列に１であるユニット
が１つずつあるときに最小となる第１の関数と、脈動圧
力値と基準値との差を評価する第２の関数との加重加算
によってネットワークのエネルギー関数を定義し、模擬
焼き鈍しの方法によって収束したときのユニットの状態
によって、オリフィスの種類と位置とを決定する方法で
ある。

【００１０】請求項２の発明に係る方法は、オリフィス
の種類とオリフィスを挿入可能な位置及び脈動圧力値を
評価する要素とを示すユニットマトリックスを定義し、
ネットワークのエネルギー関数Ｅを次式のように定義
し、 Ｅ＝（Ａ・Ｅ１＋Ｂ・Ｅ２）／（ｎ＋１） 但し、Ｅ１：ユニットマトリックスの拘束条件を与える
関数 Ｅ２：脈動圧力値と基準値との差を評価する関数 Ａ：Ｅ１の重み Ｂ：Ｅ２の重み ｎ：脈動圧力値が基準値を超えた要素の個数 模擬焼き鈍しの方法によって収束したときのユニットの
状態によって、オリフィスの種類と位置とを決定する方
法である。

【００１１】請求項３の発明に係る方法は、圧力脈動解
析を行って基準値を超える要素を求める第１ステップ、
オリフィスの種類とオリフィスを挿入可能な位置及び脈
動圧力値を評価する要素とを示すユニットマトリックス
を定義する第２ステップ、エネルギー関数に基づいてネ
ットワークのエネルギーを計算する第３ステップ、ネッ
トワークのエネルギー状態が収束したか否かを判断し、
収束した場合には処理を終了し、収束していない場合に
は次の処理に移行する第４ステップ、各ユニット間の結
合の重み、各ユニットのしきい値、及び各ユニットへの
入力の総和を計算する第５ステップ、入力の総和に基づ
いてボルツマンの確率を求める第６ステップ、確率に応
じて各ユニットの状態を決定する第７ステップ、圧力脈
動解析を行って脈動圧力を計算する第８ステップ、第８
ステップで計算された脈動圧力が基準値を超えているか
否かを判断する第９ステップ、第９ステップで基準値を
超えていないと判断された場合に、第３ステップ以降の
処理を繰り返す第１０ステップ、及び、第９ステップで
基準値を超えていると判断された場合に、第５ステップ
で計算された結合の重み及びしきい値を１回前の値に戻
し、第６ステップ以降の処理を繰り返す第１１ステッ
プ、を含んでなる方法である。

【００１２】請求項４の発明に係る方法は、前記第７ス
テップで決定された各ユニットの状態が、前回に決定し
たユニットの状態と同じであるか否かを判断し、同じで
あれば第３ステップ以降の処理を繰り返す方法である。

【００１３】請求項５の発明に係る方法は、複数の往復
圧縮機によって運転する場合に、これらの往復圧縮機の
異なる組み合わせに対して請求項１乃至請求項４のいず
れかに記載の方法により最適のオリフィスを選定し、且
つ選定したオリフィスを可変可能に取り付けておき、往
復圧縮機の運転時の組み合わせに応じて前記オリフィス
を最適の状態に可変する方法である。

【００１４】

【作用】ボルツマンマシンは、原則としてネットワーク
のエネルギーを下げるように動作するので、最適化問題
の変数をニューロン（ユニット）の状態に、目的関数を
ネットワークのエネルギーに対応付けることができれ
ば、これらのメカニズムによって組み合わせ最適化問題
が解けることが知られている。

【００１５】本発明においては、ユニットマトリックス
の各行及び各列に１であるユニットが１つずつあるとき
に最小となる第１の関数と、脈動圧力値と基準値との差
を評価する第２の関数との加重加算によってネットワー
クの目的関数（エネルギー関数）を定義することによっ
て、オリフィスの種類及び位置の組み合わせを選定する
という最適化問題の解が得られる。

【００１６】オリフィスの種類及び位置は、収束したと
きのユニットマトリックスにおけるユニットの状態によ
り決定される。

【００１７】

【実施例】図７は本発明の方法に用いるニューロン（ユ
ニット）Ｕｉのモデルを示す図である。なお、ユニット
Ｕの添字であるｉは、使用される箇所に応じて、ｊ，ｉ
ｋ，ｊｌなどのように異ならせている。

【００１８】図７において、ユニットＵｉへの入力の総
和Ｉｉは、次の（１）式で示される。

【００１９】

【数１】

【００２０】但し、Ｗｉｊ：ユニットＵｉとユニットＵ
ｊとの結合の重み Ｓｉ ：外部からの入力 θｉ ：ユニットＵｉのしきい値 ユニットＵｉの状態（これも「Ｕｉ」で示す）は、入力
の総和Ｉｉに応じて変化する。ボルツマンマシンでは、
ニューラルネットワークにおけるユニットの状態の変化
が確率的に行われる。すなわち、入力の総和Ｉｉに基づ
いて、次の（２）式に示す確率でユニットＵｉの新しい
状態が「１」に設定される。

【００２１】

【数２】

【００２２】但し、Ｔ：ネットワークの温度（＞０） ここで、ネットワークにおけるある状態ａのエネルギ関
数Ｅ（ａ）を次の（３）式のように定義する。

【００２３】

【数３】

【００２４】ここで、Ｕｉ（ａ）は、状態ａにおけるユ
ニットＵｉの値である。上述の（２）式に示すような確
率でユニットＵｉの状態を変化させると、Ｅ（ａ）は最
小値に近づくことが知られている。

【００２５】さて、次に、本実施例のニューラルネット
ワークで用いるユニットマトリックスＵＭ及びエネルギ
ー関数Ｅについて説明する。図１は本発明に係るユニッ
トマトリックスＵＭの例を示す図、図４は配管系ＰＬの
モデルの一例を示す図である。

【００２６】図４に示す配管系ＰＬでは、２台の圧縮機
ＣＡ，ＣＢによって圧縮されたガスが、配管ＰＰ、スナ
ッバＳＮＡ，ＳＮＢ、連結要素Ｅ９、クーラーＣＬ、ス
クラバＳＲなどを経由して開放端から出力される。但
し、圧縮機ＣＡと配管ＰＰとの接続部分は閉じられた状
態である。

【００２７】この配管系ＰＬにおける圧力脈動を解析す
るために、先に述べた公知の圧力脈動解析装置（ＰＵＬ
ＳＡＳ）が用いられる。圧力脈動解析装置は、有限要素
法を適用して配管内における流体の圧力脈動を解析し、
配管系ＰＬの各位置（要素）における基本波成分及び高
調波成分の圧力脈動の数値解を得るものである。

【００２８】圧力脈動解析装置によると、例えば図５に
示すように、各周波数成分毎に、横軸に沿った位置（要
素）に対する脈動圧力値が示される。破線は基準値（Ａ
ＰＩ規格値）ＶＳＴである。また、図６に示すように、
横軸を周波数、縦軸を圧力として、基本波及びｎ次高調
波の脈動圧力値が実線の棒グラフで示され、基準値ＶＳ
Ｔが破線で示される。

【００２９】図１において、ユニットマトリックスＵＭ
は、オリフィスの種類Ｏ１，Ｏ２，…Ｏｋが行方向に、
オリフィスの位置（流量節点）Ｑ１，Ｑ２，…Ｑｍ、及
び脈動圧力値が基準値を超える要素Ｐ１，Ｐ２，…Ｐｎ
が列方向に、それぞれ配置されて構成されている。

【００３０】オリフィスの種類として、例えば、縮径比
が１／２、１／３、１／４などのものをＯ１，Ｏ２，
…，Ｏｋ−１というように対応付け、縮径比が１／１で
あるダミーオリフィスをＯｋに対応付け、合計ｋ種類の
オリフィスを用いることとしている。しかし、各種類の
オリフィスをｍ個の各節点（流量節点）Ｑ１，Ｑ２，…
Ｑｍに対応付けることを可能とするため、各種類のオリ
フィスをｍ個ずつ準備し、合計（ｋ×ｍ）個のオリフィ
スを配置して（ｋ×ｍ）の列を設けている。なお、ダミ
ーオリフィスは、節点Ｑにオリフィスを挿入しないとき
に使用するためのものであり、これをｍ個準備すること
によって、実際に挿入するオリフィスが１個以下である
場合にも対応可能なようになっている。

【００３１】オリフィスの位置として、配管系ＰＬ上に
とったｍ個の節点Ｑ１，Ｑ２，…Ｑｍが準備されてい
る。要素Ｐ１，Ｐ２，…Ｐｎは、圧力脈動解析を行って
脈動圧力値が基準値を超えるｌ個の要素Ｐ１，Ｐ２，…
Ｐｌを得た場合に、実際にはこれらの全部の要素を用い
るのではなく、各次数の振動モードを考慮することによ
ってそれよりも少ない要素Ｐ１，Ｐ２，…Ｐｎ（ｎ＜
ｌ）としたものである。

【００３２】このようにして、オリフィスを挿入する位
置（行）と種類（列）を示すユニットマトリックスＵＭ
ができ上がる。なお、ユニットマトリックスＵＭを正方
行列とするために以下の操作を行う。 （１）ｎ＋ｍ≧ｍ×ｋの場合には、（ｎ＋ｍ−ｍ×ｋ）
個のＯｋを追加する。 （２）ｎ＋ｍ＜ｍ×ｋの場合には、〔ｍ×ｋ−（ｎ＋
ｍ）〕個の流量節点Ｑｍ＋１，Ｑｍ＋２…を追加する。

【００３３】ユニットマトリックスＵＭにおいて、ユニ
ットＵｉの値が「１」である場合に、それに該当する位
置Ｑ１，Ｑ２，…Ｑｍに該当する種類Ｏ１，Ｏ２，…，
Ｏｋのオリフィスを挿入することを示す。

【００３４】このようなユニットマトリックスＵＭに対
して、各要素の脈動圧力値と基準値との差の２乗が最小
になるときにネットワークのエネルギが最小となるよう
に、ユニットＵｉの結合を決める。

【００３５】まず、次の２つの関数Ｅ１，Ｅ２を定義す
る。 （１）ユニットマトリックスＵＭの拘束条件を与える関
数Ｅ１を次の（４）式に示すように定める。関数Ｅ１
は、ユニットマトリックスＵＭの各行各列に「１」であ
るユニットＵｉがそれぞれ１つずつあるときに最小値が
０となる関数である。

【００３６】

【数４】

【００３７】（２）要素Ｐｉでの脈動圧力値と基準値と
の差を評価する関数Ｅ２を次の（５）式に示すように定
める。

【００３８】

【数５】

【００３９】但し、係数Ｃｉは次の（６）式によって示
される。 Ｃｉ＝（Ａｉ−ｐｉ）² ……（６） ここで、Ａｉ及びｐｉは要素Ｐｉでの基準値および脈動
圧力値である。また、節点に関する係数Ｃｋは「０」と
する。

【００４０】次に、エネルギー関数（目的関数）Ｅを、
次の（７）式に示すように、上述した２つの関数Ｅ１，
Ｅ２の重み付けの和として定義する。 Ｅ＝（Ａ・Ｅ１＋Ｂ・Ｅ２）／（ｎ＋１） ……（７） ここで、ｎは脈動圧力値が基準値を超えた要素の個数で
ある。

【００４１】さらに、γ＝Ａ／Ｂ、δ＝Ｂ／Ｔ０（Ｔ０
は初期温度）とすると、（７）式は次の（８）式のよう
に変形できる。 Ｅ＝Ｂ（γＥ１＋Ｅ２）／（ｎ＋１） ……（８） この（８）式を展開すると、次の（９）式のようにな
る。

【００４２】

【数６】

【００４３】ここで、δｉｊ、δｋｌはクロネッカのデ
ルタである。（３）式で外部からの入力Ｓｉを「０」と
したときのエネルギ関数Ｅ（ａ）の式と（９）式とを比
較すると、ネットワークにおけるユニットＵｉｋとユニ
ットＵｊｌとの結合の重みＷｉｋｊｌ、及びユニットＵ
ｉｋのしきい値θｉｋは、次の（１０）式及び（１１）
式のようになる。

【００４４】

【数７】

【００４５】これらの（１０）式及び（１１）式を
（１）式に代入することにより、ユニットＵｉｋへの入
力の総和Ｉｉｋが求められる。この入力の総和Ｉｉｋに
基づいて、（２）式に示す確率によってユニットＵｉｋ
の新しい状態を設定していく。

【００４６】このとき、模擬焼き鈍しの方法を適用す
る。つまり、ネットワークのエネルギーを最小値付近に
到達させるために、ネットワークの状態変化を行わせな
がら、ネットワークの温度Ｔを次の（１２）式に示すよ
うに除々に下げていく。

【００４７】Ｔ＝Ｔ０ｅｘｐ（Ｄ・Ｌ） ……（１２） ここで、Ｔ０は先に述べたように初期温度（＝Ｂ／δ）
であり、Ｄは負の定数（本実施例では「−０．０２」と
した）であり、Ｌは状態変化の繰り返し数である。

【００４８】次に、図２及び図３に示すフローチャート
に基づいてオリフィス選定のための処理又は操作につい
て説明する。まず、収束計算の繰り返し回数を表す変数
Ｌを「０」に初期化し（＃１）、圧力脈動解析装置によ
り圧力脈動解析を行って基準値を超える要素Ｐ１，Ｐ
２，…Ｐｎを求める（＃２）。

【００４９】オリフィスの種類Ｏ１，Ｏ２，…，Ｏｋ及
び取り付け可能な位置Ｑ１，Ｑ２，…Ｑｍを入力し（＃
３）、ユニットマトリックスＵＭの各ユニットＵｉｋの
初期値を入力して「０」に初期化し（＃４）、変数Ｌに
１を加算する（＃５）。ネットワークの温度Ｔを（１
２）式に基づいて求め（＃６）、γを自動算出する（＃
７）。

【００５０】ここで、γは２つの関数Ｅ１，Ｅ２の重み
の比（Ａ／Ｂ）であり、各δに対して最適な値がプログ
ラムによって自動的に決定される。すなわち、収束計算
を行っていく中で、関数Ｅ１が「０」にならない場合に
は、関数Ｅ１に関する重みＡが小さ過ぎると見なされ、
γがγ＝ａ・Ｈにしたがって増加するように処理され、
これによって関数Ｅ１の重みＡを増加させる。なお、Ｈ
はＥ１が「０」でない回数である。また、ａは２つの関
数Ｅ１，Ｅ２の比に関連する定数であり、２つの関数Ｅ
１，Ｅ２のバランスが保持されるように、例えば、０．
００１〜０．００４程度の値が用いられる。

【００５１】エネルギー関数Ｅの計算を行い（＃８）、
収束したか否かの判断を行う（＃９）。本実施例では、
全部のユニットＵｉｋの状態が連続１０回にわたって何
の変化もない場合には、ネットワークがその状態に収束
したものと判断し、処理を終了する。

【００５２】収束していない場合には、結合の重みＷｉ
ｋｊｌ及びしきい値θｉｋを（１０）式及び（１１）式
に基づいて計算し、（１）式に基づいて入力の総和Ｉｉ
ｋを求める（＃１１）。

【００５３】次に、入力の総和Ｉｉｋから上述の（２）
式に基づいて確率Ｇｉｋを求め（＃１２）、求めた確率
Ｇｉｋの値をその都度発生させた乱数Ｒと比較する（＃
１３）。

【００５４】確率Ｇｉｋが乱数Ｒよりも大きいときには
（＃１３でイエス）、当該ユニットＵｉｋの状態を
「１」とし（＃１４）、そうでない場合には当該ユニッ
トＵｉｋの状態を「０」とする（＃１５）。

【００５５】ユニットマトリックスＵＭの状態、つまり
オリフィスの組み合わせが前回と異なっている場合には
（＃１６でイエス）、圧力脈動解析を行って脈動圧力値
を計算し（＃１７）、計算した脈動圧力値が基準値を超
えていない場合には（＃１８でノー）、ステップ＃１１
で計算された結合の重みＷｉｋｊｌ及びしきい値θｉｋ
をレジスタＷ，θにそれぞれ保持し（＃１９）、ステッ
プ＃５以降を繰り返す。

【００５６】脈動圧力値が基準値を超えている場合には
（＃１８でイエス）、結合の重みＷｉｋｊｌ及びしきい
値θｉｋを、前回に計算してレジスタＷ，θに保持した
値とし、且つ入力の総和Ｉｉを算出し（＃２０）、ステ
ップ＃１２以降を繰り返す。つまり、この場合には、結
合の重みＷｉｋｊｌ及びしきい値θｉｋを更新すること
なく、新たに発生した乱数Ｒを用いてユニットＵｉｋの
状態を決定し、これによって新たなオリフィスの組み合
わせを選出する。

【００５７】前回と同一のオリフィスの組み合わせが選
出されている場合には（＃１６でイエス）、ステップ＃
５以降を繰り返す。ステップ＃５以降が繰り返されるこ
とにより、ネットワークの温度Ｔが徐々に低下していき
ながら、エネルギー関数Ｅが最小値に近づいていき、ユ
ニットマトリックスＵＭの各行各列のユニットＵｉｋが
それぞれ１つだけ「１」となり、且つ脈動圧力値と基準
値との差が最小となる。脈動圧力値が基準値を超える要
素が含まれていた場合には、ステップ＃１２以降が繰り
返されることによってそのようなオリフィスの組み合わ
せが自動的に取り除かれる。つまり、基準値を超える要
素がなくなるまで、次々と発生された乱数Ｒとの比較に
よってユニットＵｉｋの状態が決定され、ユニットマト
リックスＵＭの状態が変化する。

【００５８】このように、フローチャートの内容が繰り
返して実行されることにより、脈動圧力値が基準値を超
えない範囲で基準値に近づくとともに、予め設定した各
節点Ｑ１，Ｑ２，…Ｑｍに対して挿入すべきオリフィス
の種類Ｏ１，Ｏ２，…，Ｏｋが自動的に決定される。

【００５９】したがって、収束した後のユニットマトリ
ックスＵＭにより示される位置及び種類のオリフィスを
選定し、それらを配管系ＰＬに実際に挿入することによ
って、脈動圧力値を基準値以下に抑え且つ効率の低下を
最小限に抑えることができる。

【００６０】図６はオリフィスを取り付ける前後におけ
る脈動圧力値の状態の一例を示す図である。これらの図
で、図６（ａ）はある要素Ｐａにおけるオリフィス取り
付け前の状態、図６（ｂ）はある要素Ｐａにおけるオリ
フィス取り付け後の状態、図６（ｃ）はある要素Ｐｂに
おけるオリフィス取り付け前の状態、図６（ｄ）はある
要素Ｐｂにおけるオリフィス取り付け後の状態を示す。

【００６１】図６（ｂ）（ｄ）に示されるように、基準
値ＶＳＴを越えることなく基準値ＶＳＴに近い解を選定
していることが分かる。ところで、上述した方法によ
り、配管系ＰＬにおける種々の運転状態におけるオリフ
ィスの最適な位置と種類を選定し、且つ選定したオリフ
ィスを可変可能に取り付けておき、実際の運転状態に応
じてオリフィスを最適の状態に可変調整することが可能
である。

【００６２】例えば、図４に示すような圧縮機ＣＡと配
管ＰＰとの接続部分が閉じられた運転状態の他に、圧縮
機ＣＡと配管ＰＰとの接続部分が開かれた運転状態、さ
らには圧縮機ＣＢと配管ＰＰとの接続部分が閉じられた
運転状態などにおいて、それぞれ最適のオリフィスの位
置と種類を選定しておき、それぞれ選定された状態を実
現することが可能なように可変型のオリフィスを取り付
けておくことによって、実際の種々の運転状態に合わせ
てオリフィスを可変調整し、脈動圧力値を基準値以下に
抑え且つ効率の低下を最小限に抑えることができる。

【００６３】上述の実施例においては、脈動圧力値をＡ
ＰＩ規格値で定められた基準値ＶＳＴ以下に抑えるよう
にした例について説明したが、基準値ＶＳＴとして他の
適当な値を設定することが可能である。また、脈動圧力
値を基準値以下に抑えることなく、単に基準値にできる
だけ近づけたい場合にも適用することができる。

【００６４】上述の実施例において、ボルツマンマシン
は、例えば高速の汎用コンピュータ上でシミュレーショ
ンを行うことにより、また専用のニューロチップを用い
ることにより、それぞれ実現することができる。ユニッ
トマトリックスＵＭの行及び列の数は、配管系ＰＬの状
態、オリフィスの種類、節点又は要素のとり方などに応
じて種々変更することができる。フローチャートの内容
又は順序は本発明の主旨に沿って変更することが可能で
ある。

【００６５】

【発明の効果】本発明によると、配管系におけるオリフ
ィスの種類及び位置の最適の組み合わせを自動的に選定
することができ、高度の熟練度や工学的センスを有して
いなくても少ない時間で最適のオリフィスを選定するこ
とができる。

【図面の簡単な説明】

【図１】本発明に係るユニットマトリックスの例を示す
図である。

【図２】オリフィス選定のための処理又は操作を示すフ
ローチャートである。

【図３】オリフィス選定のための処理又は操作を示すフ
ローチャートである。

【図４】配管系のモデルの一例を示す図である。

【図５】オリフィスを取り付ける前におけるある周波数
成分の脈動圧力値の状態の解析結果を示す図である。

【図６】オリフィスを取り付ける前後における脈動圧力
値の状態の解析結果の一例を示す図である。

【図７】本発明の方法に用いるニューロンのモデルを示
す図である。

【符号の説明】

ＵＭ ユニットマトリックス Ｏ１，Ｏ２，Ｏｋ 種類（オリフィスの種類） Ｑ１，Ｑ２，Ｑｍ 節点（オリフィスの位置） Ｐ１，Ｐ２，Ｐｎ 要素（脈動圧力値を評価する要素）

Claims (5)

【特許請求の範囲】
1. 【請求項１】配管系における圧力脈動の低減のために前
   記配管系内に挿入するオリフィスの種類と位置とをボル
   ツマンマシンを用いて選定する方法であって、 オリフィスの種類とオリフィスを挿入可能な位置とを示
   すユニットマトリックスを定義し、 前記ユニットマトリックスの各行及び各列に１であるユ
   ニットが１つずつあるときに最小となる第１の関数と、
   脈動圧力値と基準値との差を評価する第２の関数との加
   重加算によってネットワークのエネルギー関数を定義
   し、 模擬焼き鈍しの方法によって収束したときのユニットの
   状態によって、オリフィスの種類と位置とを決定するこ
   とを特徴とする配管系における圧力脈動の低減のための
   オリフィス選定方法。
2. 【請求項２】配管系における圧力脈動の低減のために前
   記配管系内に挿入するオリフィスの種類と位置とをボル
   ツマンマシンを用いて選定する方法であって、 オリフィスの種類とオリフィスを挿入可能な位置及び脈
   動圧力値を評価する要素とを示すユニットマトリックス
   を定義し、 ネットワークのエネルギー関数Ｅを次式のように定義
   し、 Ｅ＝（Ａ・Ｅ１＋Ｂ・Ｅ２）／（ｎ＋１） 但し、Ｅ１：ユニットマトリックスの拘束条件を与える
   関数 Ｅ２：脈動圧力値と基準値との差を評価する関数 Ａ：Ｅ１の重み Ｂ：Ｅ２の重み ｎ：脈動圧力値が基準値を超えた要素の個数 模擬焼き鈍しの方法によって収束したときのユニットの
   状態によって、オリフィスの種類と位置とを決定するこ
   とを特徴とする配管系における圧力脈動の低減のための
   オリフィス選定方法。
3. 【請求項３】配管系における圧力脈動の低減のために前
   記配管系内に挿入するオリフィスの種類と位置とをボル
   ツマンマシンを用いて選定する方法であって、 圧力脈動解析を行って基準値を超える要素を求める第１
   ステップ、 オリフィスの種類とオリフィスを挿入可能な位置及び脈
   動圧力値を評価する要素とを示すユニットマトリックス
   を定義する第２ステップ、 エネルギー関数に基づいてネットワークのエネルギーを
   計算する第３ステップ、 ネットワークのエネルギー状態が収束したか否かを判断
   し、収束した場合には処理を終了し、収束していない場
   合には次の処理に移行する第４ステップ、 各ユニット間の結合の重み、各ユニットのしきい値、及
   び各ユニットへの入力の総和を計算する第５ステップ、 入力の総和に基づいてボルツマンの確率を求める第６ス
   テップ、 確率に応じて各ユニットの状態を決定する第７ステッ
   プ、 圧力脈動解析を行って脈動圧力を計算する第８ステッ
   プ、 第８ステップで計算された脈動圧力が基準値を超えてい
   るか否かを判断する第９ステップ、 第９ステップで基準値を超えていないと判断された場合
   に、第３ステップ以降の処理を繰り返す第１０ステッ
   プ、 及び、 第９ステップで基準値を超えていると判断された場合
   に、第５ステップで計算された結合の重み及びしきい値
   を１回前の値に戻し、第６ステップ以降の処理を繰り返
   す第１１ステップ、 を含んでなることを特徴とする配管系における圧力脈動
   の低減のためのオリフィス選定方法。
4. 【請求項４】請求項３記載のオリフィス選定方法におい
   て、 前記第７ステップで決定された各ユニットの状態が、前
   回に決定したユニットの状態と同じであるか否かを判断
   し、同じであれば第３ステップ以降の処理を繰り返すこ
   とを特徴とする配管系における圧力脈動の低減のための
   オリフィス選定方法。
5. 【請求項５】複数の往復圧縮機によって運転する場合
   に、これらの往復圧縮機の異なる組み合わせに対して請
   求項１乃至請求項４のいずれかに記載の方法により最適
   のオリフィスを選定し、且つ選定したオリフィスを可変
   可能に取り付けておき、往復圧縮機の運転時の組み合わ
   せに応じて前記オリフィスを最適の状態に可変すること
   を特徴とする配管系における圧力脈動の低減のためのオ
   リフィス選定方法。