(57)【要約】【課題】距離差−周波数分析を用いた配管の漏洩感知装置および方法が提供される。【解決手段】距離差−周波数分析を用いた配管の漏洩感知装置は、配管の長さ方向に沿って予め設定された間隔で設置され、配管から伝達される振動信号を感知する振動感知センサーと、振動信号を用いて相互相関関数を算出し、相互相関関数に対する時間遅延−周波数分析を行って算出された相互相関関数の周波数別時間遅延情報を算出し、相互相関関数の周波数別時間遅延情報のそれぞれの周波数成分に対してそれぞれ周波数別伝播速度情報を用いた補正変換を行って距離差−周波数分析図を算出し、距離差−周波数分析図から前記配管の漏洩成分を感知する制御部とを含む。【選択図】図１
2019219377.tif 000002
【特許請求の範囲】
【請求項１】
配管の長さ方向に沿って予め設定された間隔で設置され、前記配管から伝達される振動信号を感知する振動感知センサーと、
前記振動信号を用いて相互相関関数を算出し、前記相互相関関数に対する時間遅延−周波数分析を行って相互相関関数の周波数別時間遅延情報を算出し、前記相互相関関数の周波数別時間遅延情報のそれぞれの周波数成分に対してそれぞれ周波数別伝播速度情報を用いた補正変換を行って距離差−周波数分析図を算出し、前記距離差−周波数分析図から前記配管の漏洩成分を感知する制御部と
を含む、距離差−周波数分析を用いた配管の漏洩感知装置。
【請求項２】
前記制御部は、前記配管の漏洩成分を感知して漏洩位置表示信号を発生する、請求項１に記載の距離差−周波数分析を用いた配管の漏洩感知装置。
【請求項３】
前記漏洩位置表示信号に基づいて前記配管の漏洩位置を表示する表示部をさらに含む、請求項２に記載の距離差−周波数分析を用いた配管の漏洩感知装置。
【請求項４】
配管の長さ方向両側にそれぞれ設置された一対の振動感知センサーから振動信号（ｘ１（ｔ）、ｘ２（ｔ））を計測する段階と、
前記振動信号を用いて相互相関関数Ｃ１２（τ）を算出する段階と、
前記相互相関関数に対する時間遅延−周波数分析を行って、算出された相互相関関数の周波数別時間遅延情報Ｃ１２（τ、ｆ）を算出する段階と、
前記相互相関関数の周波数別時間遅延情報Ｃ１２（τ、ｆ）のそれぞれの周波数成分に対してそれぞれ周波数別伝播速度情報ｃ（ｆ）を用いた補正変換δ＝ｃτを行って距離差−周波数分析図Ｄ１２（δ、ｆ）を算出する段階と、
前記距離差−周波数分析図Ｄ１２（δ、ｆ）から前記配管の漏洩成分を感知する段階と
を含む、距離差−周波数分析を用いた配管の漏洩感知方法。
【請求項５】
前記相互相関関数Ｃ１２（τ）は下記の式によって算出する、請求項４に記載の距離差−周波数分析を用いた配管の漏洩感知方法。
［数１］
Ｃ１２（τ）＝Ｅ［ｘ１（ｔ）ｘ２（ｔ＋τ）］
（ここで、Ｅ［］は期待値、τは時間遅延を示す）
【請求項６】
前記周波数別伝播速度情報ｃ（ｆ）は下記の式によって算出する、請求項４に記載の距離差−周波数分析を用いた配管の漏洩感知方法。
【数２】
2019219377.tif 000004
（ここで、ｃｆは流体の音速、Ｂｆは流体の体積弾性係数、ｒは配管の半径、ｈは配管の厚さ、Ｅは配管材質のヤング率、ρは配管の密度、ｆは周波数を示す）
【請求項７】
前記時間遅延−周波数分析は、相互相関関数に対するウィグナーヴィル分布または短時間フーリエ変換の方法を用いる、請求項４に記載の距離差−周波数分析を用いた配管の漏洩感知方法。
【請求項８】
前記配管の漏洩成分に基づいて前記配管の漏洩位置を表示する段階をさらに含み、
前記配管の漏洩位置を表示する段階で垂直線エッジ検出方法を用いる、請求項４に記載の距離差−周波数分析を用いた配管の漏洩感知方法。
【発明の詳細な説明】
【技術分野】
【０００１】
距離差−周波数分析を用いた配管の漏洩感知装置および方法が提供される。
【背景技術】
【０００２】
発電所を含むプラントの場合、初期建設時に設置された配管の老朽化によって配管の漏洩が発生することがある。配管の漏洩探知のために主に使用される方法は、相互相関関数（ｃｒｏｓｓ−ｃｏｒｒｅｌａｔｉｏｎ  ｆｕｎｃｔｉｏｎ）を用いた方法である。発電所のようにポンプ、電動機などが持続的に運転される状況では機械運転による雑音成分が相互相関関数に影響を与える。これにより、一般的な相互相関関数技法を使用する場合、漏洩の真偽判別が難しいことがある。
【０００３】
プラントの機械運転による周辺雑音環境でクロススペクトル（ｃｒｏｓｓ−ｓｐｅｃｔｒｕｍ）の位相情報傾きを用いて漏洩を推定する技術が開発された。しかし、このような方法は漏洩が一箇所である場合に効果的に適用され、複数の漏洩が発生する場合には適用が難しいことがある。したがって、機械運転雑音が存在する環境で複数の漏洩が存在する場合にも適用可能な漏洩の真偽を判別する技術が要求されている。
【０００４】
関連先行文献として、日本公開特許第１９９９−２０１８５８号（特許文献１）は「相関法を用いる管路系の振動測定方法」を開示し、日本登録特許第３２９５６８２号（特許文献２）は「配管漏洩位置特定方法における信号処理方法」を開示し、日本登録特許第３６３９１６０号（特許文献３）は「異常箇所検出装置」を開示する。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開平１１−２０１８５８号公報
【特許文献２】特許第３２９５６８２号公報
【特許文献３】特許第３６３９１６０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
本発明の一実施形態は、配管に機械運転雑音が伝達される環境で複数の漏洩が存在する場合にも距離差−周波数分析を用いて漏洩の真偽判別に対する信頼度を向上させることができる距離差−周波数分析を用いた配管の漏洩感知装置および方法を提供するためのものである。
【課題を解決するための手段】
【０００７】
距離差−周波数分析を用いた配管の漏洩感知装置は、配管の長さ方向に沿って予め設定された間隔で設置され、前記配管から伝達される振動信号を感知する振動感知センサーと、振動信号を用いて相互相関関数を算出し、相互相関関数に対する時間遅延−周波数分析を行って相互相関関数の周波数別時間遅延情報を算出し、相互相関関数の周波数別時間遅延情報のそれぞれの周波数成分に対してそれぞれ周波数別伝播速度情報を用いた補正変換を行って距離差−周波数分析図を算出し、距離差−周波数分析図から前記配管の漏洩成分を感知する制御部とを含む。
【発明の効果】
【０００８】
本発明の一実施形態によれば、配管の漏洩が複数で発生する場合にも距離差−周波数分析を用いて配管の漏洩信号に対する特性を視覚化して漏洩真偽判別および漏洩推定をより正確に行うことができる。
【０００９】
本発明の一実施形態によれば、配管に機械運転雑音が伝達される環境でも、漏洩成分は垂直線形態で表現され、回転体の機械運転雑音は水平線形態で表現されるので、より効率的に漏洩を推定することができる。
【図面の簡単な説明】
【００１０】
【図１】本発明の一実施形態による距離差−周波数分析を用いた配管漏洩感知装置を示した図である。
【図２】本発明の一実施形態による距離差−周波数分析を用いた配管漏洩感知方法を示したフローチャートである。
【図３】本発明の一実施形態による距離差−周波数分析を用いた配管漏洩成分感知過程を表示部に表示したフローチャートである。
【図４】本発明の一実施形態による配管に伝達される回転体雑音信号がある場合に距離差−周波数分析を用いた配管漏洩成分感知過程を表示部に表示したフローチャートである。
【図５】本発明の一実施形態による配管の漏洩が複数である場合に距離差−周波数分析を用いた配管漏洩成分感知過程を表示部に表示したフローチャートである。
【図６】本発明の一実施形態による配管に伝達される回転体雑音信号がある状況で配管の漏洩が複数である場合に距離差−周波数分析を用いた配管漏洩成分感知過程を表示部に表示したフローチャートである。
【発明を実施するための形態】
【００１１】
添付の図面を参照して本発明の実施形態について本発明の属する技術分野における通常の知識を有する者が容易に実施できるように詳しく説明する。本発明は様々な異なる形態に実現でき、ここで説明する実施形態に限定されない。図面において本発明を明確に説明するために説明上不必要な部分は省略し、明細書全体にわたって同一または類似の構成要素については同一図面符号を使用した。また、広く知られている公知技術の場合はその具体的な説明は省略する。
【００１２】
明細書全体で、ある部分がある構成要素を「含む」という時、これは特に反対になる記載がない限り、他の構成要素を除くのではなく他の構成要素をさらに含むことができるのを意味する。
【００１３】
以下、本発明の一実施形態による配管漏洩感知装置および方法について説明する。
【００１４】
図１は、本発明の一実施形態による距離差−周波数分析を用いた配管漏洩感知装置を示した図である。図１を参照すれば、本発明の一実施形態による距離差−周波数分析を用いた配管の漏洩感知装置は、振動感知センサー１０、制御部３０、そして表示部４０を含み、距離差−周波数分析を用いた配管１００の漏洩成分を感知して、機械運転雑音が存在する環境で複数の漏洩が存在する場合にも、従来の技術による方法に比べて漏洩真偽判別の信頼度が向上できる。
【００１５】
まず、振動感知センサー１０は、配管１００の長さ方向に沿って予め設定された間隔をおいて設置されて配管１００から伝達される振動信号を感知する。振動感知センサー１０は、配管１００で漏洩が発生すると予想される漏洩位置（Ｐ）の両側に位置するように装着され得る。振動感知センサー１０は、配管１００の長さ方向に沿って両端にそれぞれ設置され得る。振動感知センサー１０は、配管１００の漏洩位置（Ｐ）および回転体１１０から伝播される信号が混合された振動信号（ｘ１（ｔ）、ｘ２（ｔ））を感知することができる。この時、漏洩位置（Ｐ）の把握のために超音波信号またはマイクロフォンを用いて測定された音響信号が適用され得るが、これに限定されない。
【００１６】
制御部３０は、振動感知センサー１０から計測された振動信号を用いて相互相関関数（ｃｒｏｓｓ−ｃｏｒｒｅｌａｔｉｏｎ  ｆｕｎｃｔｉｏｎ）を算出する。制御部３０は、算出された相互相関関数に対する時間遅延−周波数分析（ｔｉｍｅ  ｄｅｌａｙ−ｆｒｅｑｕｅｎｃｙ  ａｎａｌｙｓｉｓ）を行って算出された相互相関関数の周波数別時間遅延情報を算出する。そして、制御部３０は、算出された相互相関関数の周波数別時間遅延情報のそれぞれの周波数成分に対してそれぞれ周波数別伝播速度情報を用いた補正変換を行って距離差−周波数分析図を算出する。制御部３０は、算出された距離差−周波数分析図から配管１００の漏洩成分を感知して漏洩位置表示信号を発生する。制御部３０は配管１００の漏洩位置（Ｐ）感知と処理に関する全般的な制御動作を行うマイクロプロセッサを含むことができる。制御部３０は、一般的なコンピュータと連結されて関連情報を共有することができる。制御部３０は、別途のメモリ部を備えることができる。メモリ部は、制御部３０で演算処理された情報と配管１００の漏洩位置（Ｐ）感知と制御のために関連するデータを保存することができる。
【００１７】
表示部４０は、制御部３０から入力される漏洩位置表示信号に基づいて配管１００の漏洩位置（Ｐ）を表示する。表示部４０は、配管１００の漏洩位置（Ｐ）に該当する異常部分の有無および特定の信号を表示することができる。表示部４０は距離差−周波数領域で視覚的に配管１００の漏洩成分を垂直線で表すことができる。そして、回転体１１０の雑音は水平線で表されるので、従来の技術による方法に比べて容易に漏洩が判別され得る。
【００１８】
即ち、表示部４０を通じて配管１００の漏洩位置（Ｐ）を検査者が目で見て容易に判断することができる。
【００１９】
図２は本発明の一実施形態による距離差−周波数分析を用いた配管漏洩感知方法を示したフローチャートであり、図３は本発明の一実施形態による距離差−周波数分析を用いた配管漏洩成分感知過程を表示部４０に表示したフローチャートであり、図４は本発明の一実施形態による配管１００に伝達される回転体１１０の雑音信号がある場合に距離差−周波数分析を用いた配管漏洩成分感知過程を表示部４０に表示したフローチャートである。
【００２０】
図１〜図４を参照して本発明の一実施形態による距離差−周波数分析を用いた配管漏洩感知方法を説明する。
【００２１】
まず、制御部３０は、配管１００の漏洩位置（Ｐ）に対する真偽を判別するために配管１００の両端に設置された振動感知センサー１０から計測された配管１００の振動信号を計測する（Ｓ２１０）。そして、制御部３０は、計測された配管１００の振動信号を用いて相互相関関数を算出する（Ｓ２２０）。制御部３０は、配管１００の両端に設置された振動感知センサー１０を用いて振動信号（ｘ１（ｔ）、ｘ２（ｔ））を計測し、計測された振動信号を用いて相互相関関数Ｃ１２（τ）を算出することができる。振動感知センサー１０は、図１に示したように、予想される漏洩位置（Ｐ）の両側に位置するように装着され、漏洩位置（Ｐ）から伝播される漏洩信号および回転体１１０から伝播される周辺雑音信号（ｎ１（ｔ）、ｎ２（ｔ））を感知する。
【００２２】
第１振動感知センサー１０ａには漏洩位置（Ｐ）から第１漏洩信号（ｓ１（ｔ））が到達し得る。また、回転体１１０から第１機械雑音信号（ｎ１（ｔ））が到達し得る。そして、第２振動感知センサー１０ｂには漏洩位置（Ｐ）から第２漏洩信号（ｓ２（ｔ））が到達し得る。また、回転体１１０から第２機械雑音信号（ｎ２（ｔ））が到達し得る。第１振動感知センサー１０ａには第１漏洩信号（ｓ１（ｔ））と第１機械雑音信号（ｎ１（ｔ））が混合された第１振動信号（ｘ１（ｔ））が到達し得る。そして、第２振動感知センサー１０ｂには第２漏洩信号（ｓ２（ｔ））と第２機械雑音信号（ｎ２（ｔ））が混合された第２振動信号（ｘ２（ｔ））が到達し得る。第１振動信号（ｘ１（ｔ））と第２振動信号（ｘ２（ｔ））は数式１で表すことができる。
［数１］
ｘ１（ｔ）＝ｓ１（ｔ）＋ｎ１（ｔ）、ｘ２（ｔ）＝ｓ２（ｔ）＋ｎ２（ｔ）
【００２３】
振動信号（ｘ１（ｔ）、ｘ２（ｔ））を通じて相互相関関数を求めることによって時間遅延（τｄ）を求めることができ、時間遅延（τｄ）を通じて配管１００の漏洩位置（Ｐ）を推定することができる。この時、数式２が適用され得る。
［数２］
ｄ１−ｄ２＝ｃτｄ
ｄ１＝（Ｄ＋ｃτｄ）／２、ｄ２＝（Ｄ−ｃτｄ）／２、
ここで、ｄ１は漏洩位置（Ｐ）から第１振動感知センサーまでの距離、ｄ２は漏洩位置（Ｐ）から第２振動感知センサーまでの距離、Ｄは第１振動感知センサーと第２振動感知センサーの間の距離、ｃは振動信号の配管１００での伝播速度である。
【００２４】
前述のように、振動信号（ｘ１（ｔ）、ｘ２（ｔ））を通じて計算される相互相関関数は数式３のように表すことができる。
［数３］
Ｃ１２（τ）＝Ｅ［ｘ１（ｔ）ｘ２（ｔ＋τ）］
ここで、Ｅ［］は期待値、τは時間遅延を示す。
【００２５】
次いで、制御部３０は相互相関関数に対する時間遅延−周波数分析を行って、相互相関関数の周波数別時間遅延情報Ｃ１２（τ、ｆ）を生成する（Ｓ２３０）。この時、時間遅延−周波数分析過程では相互相関関数に対するウィグナーヴィル分布（Ｗｉｇｎｅｒ  Ｖｉｌｌｅ  Ｄｉｓｔｒｉｂｕｔｉｏｎ）または短時間フーリエ変換（ＳＴＦＴ；Ｓｈｏｒｔ  Ｔｉｍｅ  Ｆｏｕｒｉｅｒ  Ｔｒａｎｓｆｏｒｍ）などの方法が使用され得る。実際実験を行った結果、時間遅延−周波数分析段階では周波数が高いほど垂直成分の撓みがより大きく現れる。
【００２６】
一方、制御部３０は相互相関関数の周波数別時間遅延情報を示すＣ１２（τ、ｆ）のそれぞれの周波数成分に対して、それぞれ周波数別伝播速度ｃ（ｆ）情報を用いた補正変換δ＝ｃτを行って、距離差−周波数分析図Ｄ１２（δ、ｆ）を算出する（Ｓ２４０）。
【００２７】
周波数別伝播速度は数式４のように表すことができる。
【数４】
2019219377.tif 000003
ここで、（ｃｆ：流体の音速、Ｂｆ：流体の体積弾性係数、ｒ：配管の半径、ｈ：配管の厚さ、Ｅ：配管材質のヤング率（Ｙｏｕｎｇ’ｓｍｏｄｕｌｕｓ）、ρ：配管の密度、ｆ：周波数）を示す。
【００２８】
この時、配管１００の漏洩成分は広帯域周波数成分が同時に発生する特徴を有するので、Ｄ１２（δ、ｆ）で配管１００の漏洩位置（Ｐ）に該当するδ（距離差）に対して広帯域周波数成分が存在し、これはＤ１２（δ、ｆ）で垂直線の形態で表示部４０に表される。配管１００の漏洩成分は垂直線の形態で表示部４０に表される反面、持続的に発生する機械運転などによる雑音成分は水平線形態で表示部４０に表される。したがって、表示部４０に表される垂直成分を効果的に感知することができる方法を適用して配管１００の漏洩成分が容易に検出され得る（Ｓ２５０）。例えば、表示部４０に表される垂直成分を効果的に感知することができる方法は、垂直線エッジ検出（ｅｄｇｅ  ｄｅｔｅｃｔｉｏｎ）方法がある。周波数別伝播速度情報を用いた補正変換で距離差−周波数分析図を算出し、算出された距離差−周波数分析図で垂直線エッジ検出方法を適用することによって、配管１００の漏洩位置（Ｐ）に該当する距離差に対して存在する漏洩成分が容易且つ効率的に感知され得る。表示部４０に表された垂直成分が撓みがなく垂直な場合は配管１００に漏洩が発生した位置（Ｐ）が判断され得る。但し、表示部４０に表された垂直成分が垂直な場合でなければ、例えば垂直成分が大きく撓んだ状態や斜線形態を維持する場合は配管１００に漏洩が発生していないと判断され得る。
【００２９】
図４は、本発明の一実施形態による配管に伝達される回転体の雑音信号がある場合に距離差−周波数分析を用いた配管漏洩成分感知過程を表示部４０に表示したフローチャートである。図４を参照すれば、本発明の一実施形態による距離差−周波数分析方法を適用する場合、距離差−周波数領域で視覚的に漏洩成分は垂直線で表される反面、回転体１１０雑音は水平線で表されるので、従来の技術による方法に比べて容易に配管１００の漏洩位置（Ｐ）を判別することができる。
【００３０】
しかし、従来の技術による相互相関関数を用いた方法の場合、回転体１１０運転雑音が大きくなるほど配管１００の漏洩位置（Ｐ）の感知が難しい。
【００３１】
図５は本発明の一実施形態による配管の漏洩が複数である場合に距離差−周波数分析を用いた配管漏洩成分感知過程を表示部に表示したフローチャートであり、図６は本発明の一実施形態による配管に伝達される回転体の雑音信号がある状況で配管の漏洩が複数である場合に距離差−周波数分析を用いた配管漏洩成分感知過程を表示部に表示したフローチャートである。
【００３２】
図５と図６を参照すれば、配管１００での漏洩位置（Ｐ）を感知時に漏洩位置（Ｐ）が複数である場合であるか、または回転体１１０による雑音が大きい場合にも配管１００の漏洩位置（Ｐ）が垂直成分で正確に表示され得る。配管１００の漏洩成分は広帯域周波数成分が同時に発生するので、配管１００の漏洩位置（Ｐ）に該当する距離差に対して広帯域周波数成分が存在する。これは、周波数別伝播速度情報を用いた補正変換を行って算出された距離差−周波数分析図Ｄ１２（δ、ｆ）で垂直線の形態で表示部４０に表される。配管１００での漏洩位置（Ｐ）が複数である場合であるか、または回転体１１０による雑音が大きい場合にも配管１００の漏洩成分は垂直線の形態で明確に表示部４０に表されるので、配管１００の漏洩位置（Ｐ）が正確に診断され得る。
【００３３】
以上で本発明の好ましい実施形態について詳細に説明したが、本発明の権利範囲はこれに限定されるのではなく、次の請求範囲で定義している本発明の基本概念を用いた当業者の様々な変形および改良形態も本発明の権利範囲に属するものである。
【符号の説明】
【００３４】
１０    振動感知センサー
３０    制御部
４０    表示部
１００  配管
１１０  回転体
【図１】
2019219377.tif 000005
【図２】
2019219377.tif 000006
【図３】
2019219377.tif 000007
【図４】
2019219377.tif 000008
【図５】
2019219377.tif 000009
【図６】
2019219377.tif 000010