JP5807107B1 - 分析データ作成方法、周波数フィルター作成方法、異常音発生位置の特定方法、分析データ作成装置、周波数フィルター作成装置および異常音発生位置の特定装置 - Google Patents
分析データ作成方法、周波数フィルター作成方法、異常音発生位置の特定方法、分析データ作成装置、周波数フィルター作成装置および異常音発生位置の特定装置 Download PDFInfo
- Publication number
- JP5807107B1 JP5807107B1 JP2014231581A JP2014231581A JP5807107B1 JP 5807107 B1 JP5807107 B1 JP 5807107B1 JP 2014231581 A JP2014231581 A JP 2014231581A JP 2014231581 A JP2014231581 A JP 2014231581A JP 5807107 B1 JP5807107 B1 JP 5807107B1
- Authority
- JP
- Japan
- Prior art keywords
- cross
- measurement
- correlation
- frequency
- unit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Examining Or Testing Airtightness (AREA)
Abstract
【課題】本発明の目的は、信頼性の高い分析データ作成方法、周波数フィルター作成方法、異常音発生位置の特定方法、分析データ作成装置、周波数フィルター作成装置および異常音発生位置の特定装置を提供することである。【解決手段】本発明にかかる分析データ作成方法は、少なくとも2個の振動センサ200により振動を計測する計測工程と、計測工程による計測値を周波数フィルターにより処理する周波数フィルター工程と、周波数フィルター工程により処理された計測値から異なる周波数帯域で複数の計測値を抽出する抽出工程と、抽出工程により抽出された複数の計測値から複数の相互相関関数を演算する相互相関関数演算工程と、複数の相互相関関数から所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択工程と、を含む。【選択図】図1
Description
本発明は、分析データ作成方法、周波数フィルター作成方法、異常音発生位置の特定方法、分析データ作成装置、周波数フィルター作成装置および異常音発生位置の特定装置に関する。
従来、配管の漏れの位置を決定する方法として、センサにより振動を検知し、検知された信号から相互相関関数を生成し、音響の伝搬速度を用いて異常音発生位置を特定する方法が存在する。
例えば、特許文献1(特開平8−226865号公報)には、掘削が容易でなく、雑音が多い環境において、侵入的でない方法で、邪魔な雑音源を排除して、正確に導管の漏れの位置を決定する方法について開示されている。
特許文献1(特開平8−226865号公報)記載の導管の漏れの位置を決定する方法は、a)第1時間差の生プロットを得るため、導管に沿って離間して配置された第1センサ対から得られる漏れ雑音データから相互相関関数を計算する段階と、b)第2時間差の生プロットを得るため、導管に沿って離間して配置された第2センサ対から得られる漏れ雑音データから相互相関関数を計算する段階と、c)各プロットの時間差ピークを得るために、時間差の各生プロットを平滑化する段階と、d)第1時間差ピーク、及び第1センサ対間の既知の間隔を用いることで、導管の漏れ雑音に対する伝搬速度を決定する段階と、e)伝搬速度、第2時間差ピーク、及び第2センサ対間の間隔を用いることに よって、漏れの位置を決定する段階とからなることを特徴とする導管の漏れの位置を決定するものである。
また、特許文献2(特表2003-502678号公報)には、相関に基づく技術を使って流体搬送管内での漏れを検知および位置特定する方法および装置が記載されている。
特許文献2(特表2003-502678号公報)記載の共通信号の検知と位置特定の方法においては、相関に基づく技術を使って2つの入力信号の共通信号を検知および位置特定する方法であって、周波数領域での入力信号の位相の分析により、少なくとも1個の フィルターを準備する工程と、前記の少なくとも1個のフィルターを使って周波数領域で入力信号をフィルター処理する工程と、フィルター処理した信号の相互相関を実行する工程とから成るものである。
このように、上水道またはガス管等の老朽化が進み、欠陥からの流体漏洩が問題となっている。そこで、特許文献1(特開平8−226865号公報)または特許文献2(特表2003−502678号公報)に記載の方法を適用して、欠陥位置を特定することが考えられる。
しかしながら、実際の漏洩現場においては、異常音の他に自動車の通行音等の外騒音、その他ノイズが入力されるという問題がある。
また、漏洩する流体量が少なく異常音が小さい場合、漏洩位置と振動センサ位置とが離れており、検出される異常音が小さいという問題がある。
そのため、特許文献2記載の方法では、充分なコヒーレンス性を示さない周波数を除去あるいは、阻止することで、相互相関関係のピークを増強することにより、異常音を検出しようとしている。
また、漏洩する流体量が少なく異常音が小さい場合、漏洩位置と振動センサ位置とが離れており、検出される異常音が小さいという問題がある。
そのため、特許文献2記載の方法では、充分なコヒーレンス性を示さない周波数を除去あるいは、阻止することで、相互相関関係のピークを増強することにより、異常音を検出しようとしている。
しかしながら、特許文献2記載の方法では、不完全なデジタル化の結果として極めて相関性の高いノイズが含まれてしまうため、自動帯域フィルターを併用して当該ノイズを削除している。
その削除の結果、合成樹脂管の微小な漏水のように、ノイズよりも小さな異常音が入力された場合に、必要な波形成分が除去されてしまうという問題が生じる。
その削除の結果、合成樹脂管の微小な漏水のように、ノイズよりも小さな異常音が入力された場合に、必要な波形成分が除去されてしまうという問題が生じる。
本発明の目的は、信頼性の高い分析データ作成方法、周波数フィルター作成方法、異常音発生位置の特定方法、分析データ作成装置、周波数フィルター作成装置および異常音発生位置の特定装置を提供することである。
(1)
一局面に従う分析データ作成方法は、少なくとも2個の振動センサにより振動を計測する計測工程と、計測工程による計測値を周波数フィルターにより処理する周波数フィルター工程と、周波数フィルター工程により処理された計測値から異なる周波数帯域で複数の計測値を抽出する抽出工程と、抽出工程により抽出された複数の計測値から複数の相互相関関数を演算する相互相関関数演算工程と、複数の相互相関関数から所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択工程と、を含むものである。
一局面に従う分析データ作成方法は、少なくとも2個の振動センサにより振動を計測する計測工程と、計測工程による計測値を周波数フィルターにより処理する周波数フィルター工程と、周波数フィルター工程により処理された計測値から異なる周波数帯域で複数の計測値を抽出する抽出工程と、抽出工程により抽出された複数の計測値から複数の相互相関関数を演算する相互相関関数演算工程と、複数の相互相関関数から所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択工程と、を含むものである。
この場合、振動センサにより計測された計測値に対して、異なる周波数帯域で、複数個の計測値を抽出し、抽出された複数の計測値から複数の相互相関関数を演算し、所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択することで、ノイズ成分を除去することができるため、信頼性の高いピークを有する分析データを作成することができる。
(2)
第2の発明に係る分析データ作成方法は、一局面に従う分析データ作成方法において、計測工程、周波数フィルター工程、抽出工程、相互相関関数演算工程、および相互相関関数選択工程のうち少なくとも1または複数を、複数回繰り返してもよい。
第2の発明に係る分析データ作成方法は、一局面に従う分析データ作成方法において、計測工程、周波数フィルター工程、抽出工程、相互相関関数演算工程、および相互相関関数選択工程のうち少なくとも1または複数を、複数回繰り返してもよい。
この場合、計測工程、周波数フィルター工程、抽出工程、相互相関関数演算工程、および相互相関関数選択工程の1または複数回繰り返すことにより、よりノイズ成分を除去することができるので、さらに信頼性の高い分析データを作成することができる。
(3)
第3の発明に係る分析データ作成方法は、一局面または第2の発明に係る分析データ作成方法において、周波数フィルター工程は、コヒーレンスフィルタからなり、抽出工程における異なる周波数帯域は、計測値のコヒーレンス値に基づき決定してもよい。
第3の発明に係る分析データ作成方法は、一局面または第2の発明に係る分析データ作成方法において、周波数フィルター工程は、コヒーレンスフィルタからなり、抽出工程における異なる周波数帯域は、計測値のコヒーレンス値に基づき決定してもよい。
この場合、計測値のコヒーレンス値に基づいて、コヒーレンスフィルタを形成するので、ノイズを有効に除去し、相関性の高い計測値(振動波形)のみを抽出することができる。特に、抽出工程においてコヒーレンス値に基づいて抽出作業を行うため、ノイズ除去を確実に行うことができる。
(4)
第4の発明に係る分析データ作成方法は、第3の発明に係る分析データ作成方法において、周波数フィルター工程は、コヒーレンス値を所定の周波数幅で周波数帯域を変化させつつ、周波数帯域内のコヒーレンス値の積分値が大きい順に、抽出工程の異なる周波数帯域を複数決定してもよい。
第4の発明に係る分析データ作成方法は、第3の発明に係る分析データ作成方法において、周波数フィルター工程は、コヒーレンス値を所定の周波数幅で周波数帯域を変化させつつ、周波数帯域内のコヒーレンス値の積分値が大きい順に、抽出工程の異なる周波数帯域を複数決定してもよい。
この場合、コヒーレンス値の積分値が大きい順に抽出工程を行うため、ノイズの影響が多い部分をフィルターで削除することができ、ノイズの影響の少ない部分からデータ抽出することができ、信頼性の高い分析データを作成することができる。
(5)
第5の発明に係る分析データ作成方法は、第3の発明に係る分析データ作成方法において、周波数フィルター工程は、コヒーレンス値の第n最大値(nは、整数)における周波数幅を変化させることにより、抽出工程の異なる周波数帯域を複数決定してもよい。
第5の発明に係る分析データ作成方法は、第3の発明に係る分析データ作成方法において、周波数フィルター工程は、コヒーレンス値の第n最大値(nは、整数)における周波数幅を変化させることにより、抽出工程の異なる周波数帯域を複数決定してもよい。
この場合、コヒーレンス値の第1最大値、第2最大値、第3最大値、第4最大値、第5最大値、から第n最大値までにおいて所定の周波数幅において抽出工程において異なる周波数帯域を複数決定してもよい。
なお、nが2の場合、第1最大値および第2最大値において所定の周波数帯域で抽出することができる。
その結果、コヒーレンス値が高いとは、信頼性が高い値であるため、信頼性の高い分析データを作成することができる。
なお、nが2の場合、第1最大値および第2最大値において所定の周波数帯域で抽出することができる。
その結果、コヒーレンス値が高いとは、信頼性が高い値であるため、信頼性の高い分析データを作成することができる。
(6)
他の局面に従う周波数フィルター作成方法は、少なくとも2個の振動センサにより振動を計測する計測工程と、計測工程による計測値を周波数フィルターにより処理する周波数フィルター工程と、周波数フィルター工程により処理された計測値から異なる周波数帯域で複数の計測値を抽出する抽出工程と、抽出工程により抽出された複数の計測値から複数の相互相関関数を演算する相互相関関数演算工程と、複数の相互相関関数から所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択工程と、相互相関関数選択工程において選択された頻度に応じて1または複数の周波数フィルターを決定するフィルター選択工程と、フィルター選択工程により決定された1または複数の周波数フィルターに基づき異なる周波数帯域が決定される周波数帯域決定工程を含むものである。
他の局面に従う周波数フィルター作成方法は、少なくとも2個の振動センサにより振動を計測する計測工程と、計測工程による計測値を周波数フィルターにより処理する周波数フィルター工程と、周波数フィルター工程により処理された計測値から異なる周波数帯域で複数の計測値を抽出する抽出工程と、抽出工程により抽出された複数の計測値から複数の相互相関関数を演算する相互相関関数演算工程と、複数の相互相関関数から所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択工程と、相互相関関数選択工程において選択された頻度に応じて1または複数の周波数フィルターを決定するフィルター選択工程と、フィルター選択工程により決定された1または複数の周波数フィルターに基づき異なる周波数帯域が決定される周波数帯域決定工程を含むものである。
この場合、振動センサにより計測された計測値に対して、異なる周波数帯域で、複数個の計測値を抽出し、抽出された複数の計測値から複数の相互相関関数を演算し、所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択することで、ノイズ成分を除去することができるため、信頼性の高いピークを有する周波数フィルターを作成することができる。
(7)
第7の発明に係る周波数フィルター作成方法は、他の局面に従う周波数フィルター作成方法において、計測工程、周波数フィルター工程、抽出工程、相互相関関数演算工程、相互相関関数選択工程、フィルター選択工程、および周波数帯域決定工程の少なくとも1または複数を繰り返してもよい。
第7の発明に係る周波数フィルター作成方法は、他の局面に従う周波数フィルター作成方法において、計測工程、周波数フィルター工程、抽出工程、相互相関関数演算工程、相互相関関数選択工程、フィルター選択工程、および周波数帯域決定工程の少なくとも1または複数を繰り返してもよい。
この場合、計測工程、周波数フィルター工程、抽出工程、相互相関関数演算工程、相互相関関数選択工程、フィルター選択工程、および周波数帯域決定工程の少なくとも1または複数を繰り返すので、よりノイズ成分を除去することができ、さらに信頼性の高いピークを有する周波数フィルターを作成することができる。
(8)
さらに他の発明に係る異常音発生位置の特定方法は、少なくとも2個の振動センサにより振動を計測する計測工程と、計測工程による計測値を周波数フィルターにより処理する周波数フィルター工程と、周波数フィルター工程により処理された計測値から異なる周波数帯域で複数個の計測値を抽出する抽出工程と、抽出工程により抽出された複数の計測値から複数の相互相関関数を演算する相互相関関数演算工程と、複数の相互相関関数の中から所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択工程と、相互相関関数演算工程により選択された相互相関関数に基づき、各振動センサに対する異常音の到達時間差を求め、音の発生位置を算出する算出工程と、を含むものである。
さらに他の発明に係る異常音発生位置の特定方法は、少なくとも2個の振動センサにより振動を計測する計測工程と、計測工程による計測値を周波数フィルターにより処理する周波数フィルター工程と、周波数フィルター工程により処理された計測値から異なる周波数帯域で複数個の計測値を抽出する抽出工程と、抽出工程により抽出された複数の計測値から複数の相互相関関数を演算する相互相関関数演算工程と、複数の相互相関関数の中から所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択工程と、相互相関関数演算工程により選択された相互相関関数に基づき、各振動センサに対する異常音の到達時間差を求め、音の発生位置を算出する算出工程と、を含むものである。
この場合、振動センサにより計測された計測値に対して、異なる周波数帯域で、複数個の計測値を抽出し、抽出された複数の計測値から複数の相互相関関数を演算し、所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択することで、ノイズ成分を除去することができるため、信頼性の高い音の発生位置を算出することができる。
(9)
第9の発明に係る異常音発位置の特定方法は、さらに他の発明に係る異常音発生位置の特定方法において、相互相関関数選択工程において選択された頻度に応じて、1または複数の周波数フィルターを決定するフィルター選択工程と、フィルター選択工程により決定された1または複数の周波数フィルターに基づき異なる周波数帯域が決定される周波数帯域決定工程と、を含んでもよい。
第9の発明に係る異常音発位置の特定方法は、さらに他の発明に係る異常音発生位置の特定方法において、相互相関関数選択工程において選択された頻度に応じて、1または複数の周波数フィルターを決定するフィルター選択工程と、フィルター選択工程により決定された1または複数の周波数フィルターに基づき異なる周波数帯域が決定される周波数帯域決定工程と、を含んでもよい。
この場合、振動センサにより計測された計測値に対して、異なる周波数帯域で、複数個の計測値を抽出し、抽出された複数の計測値から複数の相互相関関数を演算し、所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択することで、ノイズ成分を除去することができるため、信頼性の高い音の発生位置を算出することができる。
(10)
第10の発明に係る異常音発生位置の特定方法は、さらに他の発明または第9の発明に係る異常音発生位置の特定方法において、計測工程、周波数フィルター工程、抽出工程、相互相関関数演算工程、相互相関関数選択工程、および算出工程の少なくとも1または複数を繰り返してもよい。
第10の発明に係る異常音発生位置の特定方法は、さらに他の発明または第9の発明に係る異常音発生位置の特定方法において、計測工程、周波数フィルター工程、抽出工程、相互相関関数演算工程、相互相関関数選択工程、および算出工程の少なくとも1または複数を繰り返してもよい。
この場合、計測工程、周波数フィルター工程、抽出工程、相互相関関数演算工程、相互相関関数選択工程、および算出工程の少なくとも1または複数を繰り返すことができるので、よりノイズ成分を除去することができ、信頼性の高い音の発生位置を算出することができる。
(11)
さらに他の発明に係る異常音発生位置の特定方法は、少なくとも2個の振動センサにより振動を計測する計測工程と、計測工程による計測値を周波数フィルターにより処理する周波数フィルター工程と、周波数フィルター工程により処理された計測値を、異なる周波数帯域で複数個の計測値を抽出する抽出工程と、抽出工程により抽出された複数の計測値のから複数の相互相関関数を演算する相互相関関数演算工程と、複数の相互相関関数の中から所定の相関最大値または所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択工程と、相互相関関数選択工程で選択された相互相関関数の累積結果に基づき、各振動センサに対する異常音の到達時間差を求め、音の発生位置を算出する算出工程と、を含むものである。
さらに他の発明に係る異常音発生位置の特定方法は、少なくとも2個の振動センサにより振動を計測する計測工程と、計測工程による計測値を周波数フィルターにより処理する周波数フィルター工程と、周波数フィルター工程により処理された計測値を、異なる周波数帯域で複数個の計測値を抽出する抽出工程と、抽出工程により抽出された複数の計測値のから複数の相互相関関数を演算する相互相関関数演算工程と、複数の相互相関関数の中から所定の相関最大値または所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択工程と、相互相関関数選択工程で選択された相互相関関数の累積結果に基づき、各振動センサに対する異常音の到達時間差を求め、音の発生位置を算出する算出工程と、を含むものである。
この場合、振動センサにより計測された計測値に対して、異なる周波数帯域で、複数個の計測値を抽出し、抽出された複数の計測値から複数の相互相関関数を演算し、所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択し、累積した結果に基づき、ノイズ成分を除去することができるため、より信頼性の高い音の発生位置を算出することができる。
(12)
第12の発明に係る異常音発生位置の特定方法は、さらに他の発明に係る異常音発生位置の特定方法において、相互相関関数選択工程において選択された頻度に応じて、1または複数の周波数フィルターを決定するフィルター選択工程と、
フィルター選択工程により決定された1または複数の周波数フィルターに基づき異なる周波数帯域が決定される周波数帯域決定工程と、を含んでもよい。
第12の発明に係る異常音発生位置の特定方法は、さらに他の発明に係る異常音発生位置の特定方法において、相互相関関数選択工程において選択された頻度に応じて、1または複数の周波数フィルターを決定するフィルター選択工程と、
フィルター選択工程により決定された1または複数の周波数フィルターに基づき異なる周波数帯域が決定される周波数帯域決定工程と、を含んでもよい。
この場合、振動センサにより計測された計測値に対して、異なる周波数帯域で、複数個の計測値を抽出し、抽出された複数の計測値から複数の相互相関関数を演算し、所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択することで、ノイズ成分を除去することができるため、信頼性の高い音の発生位置を算出することができる。
(13)
第13の発明に係る異常音発生位置の特定方法は、さらに他の発明または第12の発明に係る異常音発生位置の特定方法において計測工程、周波数フィルター工程、抽出工程、相互相関関数演算工程、相互相関関数選択工程、および算出工程の少なくとも1または複数を繰り返してもよい。
第13の発明に係る異常音発生位置の特定方法は、さらに他の発明または第12の発明に係る異常音発生位置の特定方法において計測工程、周波数フィルター工程、抽出工程、相互相関関数演算工程、相互相関関数選択工程、および算出工程の少なくとも1または複数を繰り返してもよい。
この場合、計測工程、周波数フィルター工程、抽出工程、相互相関関数演算工程、相互相関関数選択工程、および算出工程の少なくとも1または複数を繰り返すことができるので、よりノイズ成分を除去することができ、信頼性の高い音の発生位置を算出することができる。
(14)
さらに他の局面に従う分析データ作成装置は、少なくとも2個の振動センサにより振動を計測する計測部と、計測部による計測値を周波数フィルターにより処理する周波数フィルター部と、周波数フィルター部により処理された計測値から異なる周波数帯域で複数の計測値を抽出する抽出部と、抽出部により抽出された複数の計測値から複数の相互相関関数を演算する相互相関関数演算部と、複数の相互相関関数から所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択部と、を含むものである。
さらに他の局面に従う分析データ作成装置は、少なくとも2個の振動センサにより振動を計測する計測部と、計測部による計測値を周波数フィルターにより処理する周波数フィルター部と、周波数フィルター部により処理された計測値から異なる周波数帯域で複数の計測値を抽出する抽出部と、抽出部により抽出された複数の計測値から複数の相互相関関数を演算する相互相関関数演算部と、複数の相互相関関数から所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択部と、を含むものである。
この場合、振動センサにより計測された計測値に対して、異なる周波数帯域で、複数個の計測値を抽出し、抽出された複数の計測値から複数の相互相関関数を演算し、所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択することで、ノイズ成分を除去することができるため、信頼性の高いピークを有する分析データを作成することができる。
(15)
さらに他の局面に従う周波数フィルター作成装置は、少なくとも2個の振動センサにより振動を計測する計測部と、計測部による計測値を周波数フィルターにより処理する周波数フィルター部と、周波数フィルター部により処理された計測値から異なる周波数帯域で複数の計測値を抽出する抽出部と、抽出部により抽出された複数の計測値から複数の相互相関関数を演算する相互相関関数演算部と、複数の相互相関関数から所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択部と、相互相関関数選択部において選択された頻度に応じて、1または複数の周波数フィルターを決定するフィルター選択部と、フィルター選択部により決定された1または複数の周波数フィルターに基づき異なる周波数帯域が決定される周波数帯域決定部と、を含むものである。
さらに他の局面に従う周波数フィルター作成装置は、少なくとも2個の振動センサにより振動を計測する計測部と、計測部による計測値を周波数フィルターにより処理する周波数フィルター部と、周波数フィルター部により処理された計測値から異なる周波数帯域で複数の計測値を抽出する抽出部と、抽出部により抽出された複数の計測値から複数の相互相関関数を演算する相互相関関数演算部と、複数の相互相関関数から所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択部と、相互相関関数選択部において選択された頻度に応じて、1または複数の周波数フィルターを決定するフィルター選択部と、フィルター選択部により決定された1または複数の周波数フィルターに基づき異なる周波数帯域が決定される周波数帯域決定部と、を含むものである。
この場合、振動センサにより計測された計測値に対して、異なる周波数帯域で、複数個の計測値を抽出し、抽出された複数の計測値から複数の相互相関関数を演算し、所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択することで、ノイズ成分を除去することができるため、信頼性の高いピークを有する周波数フィルターを作成することができる。
(16)
さらに他の局面に従う異常音発生位置の特定装置は、少なくとも2個の振動センサにより振動を計測する計測部と、計測部による計測値を周波数フィルターにより処理する周波数フィルター部と、周波数フィルター部により処理された計測値から異なる周波数帯域で複数個の計測値を抽出する抽出部と、抽出部により抽出された複数の計測値から複数の相互相関関数を演算する相互相関関数演算部と、複数の相互相関関数から所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択部と、相互相関関数演算部により選択された相互相関関数に基づき、各振動センサに対する異常音の到達時間差を求め、音の発生位置を算出する算出部と、を含むものである。
さらに他の局面に従う異常音発生位置の特定装置は、少なくとも2個の振動センサにより振動を計測する計測部と、計測部による計測値を周波数フィルターにより処理する周波数フィルター部と、周波数フィルター部により処理された計測値から異なる周波数帯域で複数個の計測値を抽出する抽出部と、抽出部により抽出された複数の計測値から複数の相互相関関数を演算する相互相関関数演算部と、複数の相互相関関数から所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択部と、相互相関関数演算部により選択された相互相関関数に基づき、各振動センサに対する異常音の到達時間差を求め、音の発生位置を算出する算出部と、を含むものである。
この場合、振動センサにより計測された計測値に対して、異なる周波数帯域で、複数個の計測値を抽出し、抽出された複数の計測値から複数の相互相関関数を演算し、所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択することで、ノイズ成分を除去することができるため、信頼性の高い音の発生位置を算出することができる。
以下、図面を参照しつつ、本発明の実施の形態について説明する。以下の説明では、同一の部品には同一の符号を付してある。それらの名称および機能も同じである。したがって、それらについての詳細な説明は繰り返さない。
<異常音発生位置の特定方法の状況説明>
図1は、異常音発生位置の特定方法の状況を説明するための模式図である。
図1は、異常音発生位置の特定方法の状況を説明するための模式図である。
図1に示すように、地中に管網110が設けられている。管網110には、一定間隔で、縦孔(マンホール)120が設けられている。本実施の形態においては、ポイントAおよびポイントBの間隔で縦孔120が設けられている。この場合、図1のポイントAおよびポイントBの縦孔120に、それぞれ振動センサ200を設ける。
<振動センサの説明>
図2は、振動センサを含む異常音の発生位置特定装置の一例を示す模式図であり、図3は図2の振動センサの特徴の一例を示す模式図である。
図2は、振動センサを含む異常音の発生位置特定装置の一例を示す模式図であり、図3は図2の振動センサの特徴の一例を示す模式図である。
図2に示すように、本実施の形態にかかる異常音の発生特定装置100は、演算装置300および少なくとも一対の振動センサ200を含む。一対の振動センサ200は、共振型の振動センサ200である。
図2の振動センサ200は、台座210、支柱220、薄膜電極230,240、リード線231,241、圧電素子250、錘260およびGPS装置262を含む。
演算装置300は、後述するフィルターを生成する生成部310、および演算部320からなる。
図2の振動センサ200は、台座210、支柱220、薄膜電極230,240、リード線231,241、圧電素子250、錘260およびGPS装置262を含む。
演算装置300は、後述するフィルターを生成する生成部310、および演算部320からなる。
図2に示すように、振動センサ200は、鉄製の台座210上に支柱220が固定される。支柱220の上端部に圧電素子250が設けられる。圧電素子250の一端部は、支柱220の上端部に片持ち支持されている。
圧電素子250の両面に銀ペーストを塗布して形成された上下一対の薄膜電極230,240が設けられる。支柱220および一対の薄膜電極230,240の間は、絶縁されている。
また、圧電素子250の他端部で、かつ薄膜電極230上に錘260が載置されている。
また、圧電素子250の他端部で、かつ薄膜電極230上に錘260が載置されている。
薄膜電極230には、リード線231が接続されており、薄膜電極240には、リード線241が接続されており、リード線231,241はそれぞれ演算装置300につながっている。
リード線231,241から出力される電位差を、コンピュータ等の処理装置により振動波形として出力する。
なお、本実施の形態においては、リード線231,241を用いることとしているが、これに限定されず、演算装置300との間で送受信可能な機能部を設けてもよい。
リード線231,241から出力される電位差を、コンピュータ等の処理装置により振動波形として出力する。
なお、本実施の形態においては、リード線231,241を用いることとしているが、これに限定されず、演算装置300との間で送受信可能な機能部を設けてもよい。
また、圧電素子250は、高分子圧電材料であるポリフッ化ビニリデンの延伸フィルム(PVDFフィルム)によって形成されている。
具体的なパラメータが、圧電材料の弾性Eと、断面二次モーメントJと、長さLと、幅bと、高さhとである場合、バネ定数kは、以下のように示される。
k=3EJ/L3(J=bh3/12)・・・(1)
圧電素子250と錘260とからなる系の共振周波数foは、以下のように示される。
fo=√(k/M)/2π・・・(2)
また、共振型の振動センサ200は、共振周波数foが、60Hz以上1000Hz未満の範囲内に少なくとも1個存在するように形成する。
本実施の形態に係る共振型の振動センサ200は、図3の実線に示すように、100Hzから500Hzまでの間に共振周波数foが4個存在するように形成している。その理由としては、管網110を伝わる異常音、特に漏水音は、可聴音が多く、中でも1000Hz未満が多いからである。
具体的には、図3に示すように、共振周波数foは、260Hz近辺のピークP1、310Hz近辺のピークP2、350Hz近辺のピークP3、480Hz近辺のピークP4の4個のピークを有する。
一方、図3の破線は、従来の振動センサの一例を示すものである。この場合、ピークを形成していない。
本実施の形態に係る共振型の振動センサ200は、図3の実線に示すように、100Hzから500Hzまでの間に共振周波数foが4個存在するように形成している。その理由としては、管網110を伝わる異常音、特に漏水音は、可聴音が多く、中でも1000Hz未満が多いからである。
具体的には、図3に示すように、共振周波数foは、260Hz近辺のピークP1、310Hz近辺のピークP2、350Hz近辺のピークP3、480Hz近辺のピークP4の4個のピークを有する。
一方、図3の破線は、従来の振動センサの一例を示すものである。この場合、ピークを形成していない。
また、本実施の形態に係る振動センサ200は、管網110のいずれの位置に設置されたかをGPS装置262により演算装置300に送信する。その結果、演算装置300は、振動センサ200の位置を認識することができる。例えば、演算装置300は、振動センサ200の埋設位置を認識できる。
<異常音発生位置の特定方法のフローチャート>
続いて、異常音発生位置の特定方法について具体例を示しつつ説明する。
続いて、異常音発生位置の特定方法について具体例を示しつつ説明する。
本実施の形態にかかる異常音発生位置の特定方法は、管網110の少なくとも2ヶ所(ポイントAおよびポイントB)に振動センサ200を設置し、管網110の欠陥等によって発生する異常音または振動を振動センサ200により検知する。
各振動センサ200に入力された波形のコヒーレンス関数を用いてフィルターを作成する。このフィルターを作成する際に、異常音発生位置の精度を高めるために後述する処理を行い、フィルターを作成し、当該フィルターを適用した後、相互相関関数から振動の伝達時間差Tdを求め、伝達時間差Tdと振動の伝搬速度Vとから異常音発生位置を特定する異常音発生位置の特定方法である。
各振動センサ200に入力された波形のコヒーレンス関数を用いてフィルターを作成する。このフィルターを作成する際に、異常音発生位置の精度を高めるために後述する処理を行い、フィルターを作成し、当該フィルターを適用した後、相互相関関数から振動の伝達時間差Tdを求め、伝達時間差Tdと振動の伝搬速度Vとから異常音発生位置を特定する異常音発生位置の特定方法である。
図1において、ポイントAの振動センサ200から距離Lの位置で流体の漏洩が発生したと仮定する。すなわち、距離Lの位置が異常音の発生位置(流体の漏洩位置)である。この場合、漏洩音は、ポイントBの振動センサ200に到達するまでにポイントAの振動センサ200の距離Lよりも距離にして距離Nだけ長い距離(L+N)を伝搬する。
したがって、ポイントAの振動センサ200およびポイントBの振動センサ200の距離をDと仮定した場合、漏洩音がポイントAの振動センサ200とポイントBの振動センサ200とに到着する伝達時間差Tdとすると、漏洩音の伝搬速度V、2つの振動センサ間の距離をDとして以下の式で求めることができる。
Td=N/V・・・(3)
また、
N=D−2L・・・(4)
で示すことができる。
また、
N=D−2L・・・(4)
で示すことができる。
式(4)を式(3)に代入することにより、
L=(D−V・Td)/2・・・(5)
と表すことができる。
以上のように距離Lを求めることができる。
L=(D−V・Td)/2・・・(5)
と表すことができる。
以上のように距離Lを求めることができる。
続いて、本実施の形態にかかる異常音発生位置の特定方法の具体例について説明する。図4は、本実施の形態にかかる異常音発生位置の特定方法の一例を示すフローチャートである。
また、図5および図6は、振動センサ200の周波数帯域における強度を示す模式図であり、図7は、振動センサ200のコヒーレンス関数の一例を示す模式図である。
また、図5および図6は、振動センサ200の周波数帯域における強度を示す模式図であり、図7は、振動センサ200のコヒーレンス関数の一例を示す模式図である。
まず、図4に示すように、管網110のポイントAの振動センサ200から漏洩音の波形を取得する(ステップS11)。同様に、管網110のポイントBの振動センサ200から漏洩音の波形を取得する(ステップS21)。
ここで、本発明に係る振動センサ200は、図3に示したように、共振周波数foが、60Hz以上1000Hz未満の間で4個あるため、流量の少ない漏れなどの振動波形を感度よく検出することができる。特に低周波数の帯域の振動波形の感度を高く維持することができる。すなわち、振動センサは、対象音の周波数が共振周波数と大きく異なる場合に、感度が極端に小さくなり、共振周波数に近づく程、感度を大きくすることができる。
なお、共振周波数foは、100Hz以上500Hz以下に少なくとも1個以上存在することが望ましい。
特に、振動センサの共振周波数foは、60Hz以上1000Hz未満の間において、好ましくは2個以上6個以下で、所定の間隔で有することが好ましいが、一定感度以下に落ち込んだ周波数帯域がある場合、当該周波数帯域に共振点を追加してもよい。
特に、振動センサの共振周波数foは、60Hz以上1000Hz未満の間において、好ましくは2個以上6個以下で、所定の間隔で有することが好ましいが、一定感度以下に落ち込んだ周波数帯域がある場合、当該周波数帯域に共振点を追加してもよい。
次いで、ポイントAの振動センサ200から取得した漏洩音の波形をフーリエ変換処理(ステップS12)し、フーリエスペクトルAを取得する(ステップS13)。同様に、ポイントBの振動センサ200から取得した漏洩音の波形をフーリエ変換処理(ステップS22)し、フーリエスペクトルBを取得する(ステップS23)。
<フーリエ変換処理(ステップS12およびステップS22)>
フーリエ変換処理に際しては、ポイントAの振動センサ200およびポイントBの振動センサ200で得られた波形から、同時刻を始点とする一定時間の波形を取り出して、その波形をフーリエ変換処理する。フーリエスペクトルをX(f)とすると、X(f)は、下記の(6)式のように、複素関数として表現される。
フーリエ変換処理に際しては、ポイントAの振動センサ200およびポイントBの振動センサ200で得られた波形から、同時刻を始点とする一定時間の波形を取り出して、その波形をフーリエ変換処理する。フーリエスペクトルをX(f)とすると、X(f)は、下記の(6)式のように、複素関数として表現される。
X(f)=∫-∞ ∞x(t)e-j2πftdt・・・(6)
式(6)は、下記の式(7)のように、実数部と虚数部とに分けて表現することが出来る。
X(f)=XR(f)+jXI(f)=|X(f)|ejθ(f)・・・(7)
式(7)において、|X(f)|は、フーリエスペクトルの振幅を表し、θ(f)は、位相を表す。フーリエスペクトルの振幅|X(f)|は、下記(8)式で求めることができる。
|X(f)|=√{XR(f)2+XI(f)2}・・・(8)
ここで、図5および図6で示すように、振動センサ200のそれぞれの周波数帯域に対する感度を取得することができる。ここで、図5および図6に示すように、振動センサ200により取得された波形は、周波数250Hzの近辺において高い値を示していることが判る。
<コヒーレンス関数処理(ステップS31)>
次いで、フーリエスペクトルAおよびフーリエスペクトルBに対してコヒーレンス関数処理を行う(ステップS31)。
次いで、フーリエスペクトルAおよびフーリエスペクトルBに対してコヒーレンス関数処理を行う(ステップS31)。
ここで、コヒーレンス関数とは、振動センサ200の信号に一定の関係が有るか否かを、0以上1以下の値で周波数毎に評価するための関数である。コヒーレンス値が1に近い程、2つの振動センサ200の関係性が大きいことを示し、2つの振動センサ200において共通の波形であると言える。
一方、コヒーレンス値が0に近い程、2つの振動センサ200のうち、いずれか一方に取得された波形である可能性が高く、その波形は、いわゆるノイズとして見なすことができるというものである。
一方、コヒーレンス値が0に近い程、2つの振動センサ200のうち、いずれか一方に取得された波形である可能性が高く、その波形は、いわゆるノイズとして見なすことができるというものである。
図8は、コヒーレンス関数処理の前段の一例を示すフローチャートであり、図9はコヒーレンス関数処理の後段の一例を示すフローチャートである。
図8に示すように、まず、一対の振動センサ200の計測値から上述したコヒーレンス関数を算出する(ステップS51)。ここで、A地点の所定時刻の波形をA−1とし、B地点の同時刻の波形をB−1とする。
次に、コヒーレンス値の高い周波数帯を所定の周波数幅で選択する(ステップS52)。
次に、コヒーレンス値の高い周波数帯を所定の周波数幅で選択する(ステップS52)。
具体的には、コヒーレンス値の高い周波数帯とは、所定の周波数幅におけるコヒーレンス値の積分値が大きい部分である。
すなわち、図7の所定の周波数幅におけるコヒーレンス値の面積が大きい部分を選択する。
本実施の形態においては、所定の周波数幅は、100Hzであるとし、面積が大きい順で、第1最大値、第2最大値、第3最大値の3か所の周波数帯を選択することとする。
すなわち、図7の所定の周波数幅におけるコヒーレンス値の面積が大きい部分を選択する。
本実施の形態においては、所定の周波数幅は、100Hzであるとし、面積が大きい順で、第1最大値、第2最大値、第3最大値の3か所の周波数帯を選択することとする。
本実施の形態におけるステップS52の処理は、以下に示す。図7に示したコヒーレンス値の積分値で大きい周波数帯は、650Hz以上750Hz以下(f1−1とする。)、750Hz以上850Hz以下(f1−2とする。)、300Hz以上400Hz以下(f1−3とする。)の3か所であるため、当該3か所を選択する。
なお、本実施の形態においては、3か所の周波数帯を選択することとしたが、これに限定されず、第P箇所(Pは任意の整数)の周波数帯を選択してもよい。
なお、本実施の形態においては、3か所の周波数帯を選択することとしたが、これに限定されず、第P箇所(Pは任意の整数)の周波数帯を選択してもよい。
ここで、選択された3か所の周波数帯を説明の上、650Hz以上750Hz以下の周波数帯を第1の周波数帯とし、750Hz以上850Hz以下の周波数帯を第2の周波数帯とし、300Hz以上400Hz以下の周波数帯を第3の周波数帯とする。
次に、一対の振動センサ200の元波形に対して選択された周波数帯の総当たりのバンドパスフィルターを通して波形を得る(ステップS53)。
ここで、具体的に、一対の振動センサ200の元波形に対して、第1周波数帯の周波数帯でバンドパスフィルターを通した波形を第1波形(A地点波形A’−1−1、B地点波形B‘−1−1)とする。
同様に、一対の振動センサ200の元波形に対して、第2周波数帯の周波数帯でバンドパスフィルターを通した波形を第2波形(A地点波形A’−1−2、B地点波形B‘−1−2)とし、第3周波数帯の周波数帯でバンドパスフィルターを通した波形を第3波形として作成する。
さらに、第1周波数帯および第2周波数帯のバンドパスフィルター、第2周波数帯および第3周波数帯のバンドパスフィルター、第1、第2、第3周波数帯のいずれも用いないバンドパスフィルター、第1、第2、第3周波数帯のすべてを用いたバンドパスフィルターを通した波形を第8波形(A地点波形A’−1−8、B地点波形B‘−1−8)として作成する。
ここで、具体的に、一対の振動センサ200の元波形に対して、第1周波数帯の周波数帯でバンドパスフィルターを通した波形を第1波形(A地点波形A’−1−1、B地点波形B‘−1−1)とする。
同様に、一対の振動センサ200の元波形に対して、第2周波数帯の周波数帯でバンドパスフィルターを通した波形を第2波形(A地点波形A’−1−2、B地点波形B‘−1−2)とし、第3周波数帯の周波数帯でバンドパスフィルターを通した波形を第3波形として作成する。
さらに、第1周波数帯および第2周波数帯のバンドパスフィルター、第2周波数帯および第3周波数帯のバンドパスフィルター、第1、第2、第3周波数帯のいずれも用いないバンドパスフィルター、第1、第2、第3周波数帯のすべてを用いたバンドパスフィルターを通した波形を第8波形(A地点波形A’−1−8、B地点波形B‘−1−8)として作成する。
次に、得られた波形に対して相互相関関数を算出する(ステップS54)。
具体的には、第1波形(A地点波形A’−1−1、B地点波形B‘−1−1)について第1相互相関関数(Td−1−1)を算出し、第2波形(A地点波形A’−1−2、B地点波形B‘−1−2)について第2相互相関関数(Td−1−2)を算出し、第3波形乃至第8波形(A地点波形A’−1−8、B地点波形B‘−1−8)について第3乃至第8相互相関関数(Td−1−8)を算出する。
具体的には、第1波形(A地点波形A’−1−1、B地点波形B‘−1−1)について第1相互相関関数(Td−1−1)を算出し、第2波形(A地点波形A’−1−2、B地点波形B‘−1−2)について第2相互相関関数(Td−1−2)を算出し、第3波形乃至第8波形(A地点波形A’−1−8、B地点波形B‘−1−8)について第3乃至第8相互相関関数(Td−1−8)を算出する。
続いて、得られた複数の相互相関関数において、ピークが明確に表れているものを1つ選択する(ステップS55)。
具体的には、ステップS54の処理において、第1相互相関関数から第8相互相関関数までの中より結果のピークが明確なものを選定する。
ここで、ピークが明確に表れているとは、8個の中からピークが最も高いものを選定する。
具体的には、ステップS54の処理において、第1相互相関関数から第8相互相関関数までの中より結果のピークが明確なものを選定する。
ここで、ピークが明確に表れているとは、8個の中からピークが最も高いものを選定する。
なお、本実施の形態においては、ピークが明確に表れているとは、8個の中からピークが最も高いものを選定することとしたが、これに限定されず、他の所定の条件を設定してもよい。
続いて、図9に示すように、ステップS51からステップS55の処理をn回数(nは任意の正の整数)繰り返した否かを判定する(ステップS55L)。n回数繰り返すまでステップS51からステップS55の処理を繰り返す。
続いて、n回数繰り返したステップS51からステップS55の処理から選択した相互相関関数の周波数帯の頻度を積算する(ステップS56)。
次いで、最大頻度の50%以上の頻度を有する周波数帯のみ通すフィルタ(バンドパスフィルター)を作成し、振動センサ200の波形に適用する(ステップS57)。
続いて、n回数繰り返したステップS51からステップS55の処理から選択した相互相関関数の周波数帯の頻度を積算する(ステップS56)。
次いで、最大頻度の50%以上の頻度を有する周波数帯のみ通すフィルタ(バンドパスフィルター)を作成し、振動センサ200の波形に適用する(ステップS57)。
具体的に、図10は、ステップS55Lの処理を32回(n=32)繰り返し、相互相関関数の周波数帯の頻度を積算した一例を示す模式図である。なお、当然のことながら、当該処理とは、図8および図9のステップS51からステップS57までの処理を意味する。
図10に示すように、10回以上(最大となった20回の半分)選択された周波数帯を選定する。
その結果、60Hz以上200Hz以下の範囲、300Hz以上410Hz以下の範囲、460Hz以上560Hz以下の範囲となる。
続いて、選択された周波数帯以外を排除するフィルターを作成し、ステップS57により作成されたフィルターを振動センサ200の元波形に適用する。
その結果、60Hz以上200Hz以下の範囲、300Hz以上410Hz以下の範囲、460Hz以上560Hz以下の範囲となる。
続いて、選択された周波数帯以外を排除するフィルターを作成し、ステップS57により作成されたフィルターを振動センサ200の元波形に適用する。
なお、ステップS57の処理において使用する元波形は、フィルター作成時に使用した一対の振動センサ200からの波形であってもよく、一対の振動センサ200から再度、ステップS11からステップS13およびステップS21からステップS23を行い、再取得した波形を用いてもよい。
以上の説明がステップS31におけるコヒーレンス関数処理である。
以上の説明がステップS31におけるコヒーレンス関数処理である。
次いで、図4に示すように、フーリエスペクトルAおよびフーリエスペクトルBを上述したフィルターに通し、フーリエ逆変換し(ステップS32)、相互相関関数(ステップS33)を求める。相互相関関数は、N乗(Nは1以外の自然数)され、これにより相互相関関数の波形のピーク点が保存された状態で、それ以外の点が圧縮される。
フーリエ逆変換することで、式(9)に示すような振幅情報を排除した波形X‘(t)を得ることができる。
x‘(t)=∫-∞ ∞X(f)e-j2πftdf・・・(9)
このx‘(t)を各振動センサ200の波形により算出し、これらの相互相関をとることによって伝達時間差Tdを求めて、位置算出処理(ステップS34)を行うことができる。
<第2の実施の形態>
次に、第2の実施の形態について説明する。第2の実施の形態においては、主に第1の実施の形態と異なる点について説明を行う。
図11は、図8のコヒーレンス関数処理の他の例を示すフローチャートである。
次に、第2の実施の形態について説明する。第2の実施の形態においては、主に第1の実施の形態と異なる点について説明を行う。
図11は、図8のコヒーレンス関数処理の他の例を示すフローチャートである。
第2の実施の形態においては、ステップS31(ステップS51からS57の処理)のコヒーレンス関数の処理において、ステップS54の処理の得られた波形に対して相互相関関数を算出した後、全ての相互相関関数の平均化を算出する(ステップS55a)。
次に、ステップS11からS13、ステップS21からS23および、ステップS51からステップS55aまでの処理を繰り返し(図9のステップS55Lと同様に)行い、さらに平均化する(ステップS56a)。
この場合、相互相関関数の処理を繰り返し行い、平均化するため、ノイズを除去することができる。
この場合、相互相関関数の処理を繰り返し行い、平均化するため、ノイズを除去することができる。
図12は、実施例の相互相関関数の一例を示す模式図であり、図13は、比較例の相互相関関数の一例を示す模式図である。
(実施例)
実施例として、振動センサ200のそれぞれの場所で振動波形を計測して、ロガーに保存した。ロガーのサンプリングレートを50kHzとし、60秒間の波形(1秒間の波形データを60個)を取得した。その波形について、コンピュータを用いて解析を行った。
実施例として、振動センサ200のそれぞれの場所で振動波形を計測して、ロガーに保存した。ロガーのサンプリングレートを50kHzとし、60秒間の波形(1秒間の波形データを60個)を取得した。その波形について、コンピュータを用いて解析を行った。
実施例においては、第1の実施の形態における分析データ作成方法、周波数フィルター作成方法、異常音発生位置の特定方法を用いて実際の相互相関関数を算出した。
(比較例)
比較例は、実施例とほぼ同じ処理で、異なる点は、本実施の形態に係るコヒーレンス関数処理を行わず、相互相関関数を算出した。
比較例は、実施例とほぼ同じ処理で、異なる点は、本実施の形態に係るコヒーレンス関数処理を行わず、相互相関関数を算出した。
図12に示すように、実施例においては、−0.05secにおいて明らかなピークが表れた。そのため、そのピーク位置が漏水位置であることが容易に認識できた。
一方、図13に示すように、比較例においては、ピークが複数表れている。具体的には、−0.05secおよび0.03secにおいてピークが表れた。
一方、図13に示すように、比較例においては、ピークが複数表れている。具体的には、−0.05secおよび0.03secにおいてピークが表れた。
以上のように、本実施の形態によると、振動センサ200により計測された計測値に対して、異なる周波数帯域で、コヒーレンス値の積分値が大きい順に抽出工程を行ない、抽出された複数の計測値から複数の相互相関関数を演算し、所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択することで、ノイズ成分を除去することができるため、信頼性の高いピークを有する分析データを作成することができる。
また、上記の処理を複数回繰り返すことにより、よりノイズ成分を除去することができるので、さらに信頼性の高い分析データを作成することができる。
また、上記の処理を複数回繰り返すことにより、よりノイズ成分を除去することができるので、さらに信頼性の高い分析データを作成することができる。
さらに、コヒーレンス値に基づいて、フィルタを形成するので、ノイズを有効に除去し、相関性の高い計測値(振動波形)のみを抽出することができる。
また、当該分析データにより信頼性の高いピークを有する周波数フィルターを作成することができる。
さらに、分析データおよび周波数フィルターを用いて、信頼性の高い異常音の発生位置を算出することができる。
さらに、分析データおよび周波数フィルターを用いて、信頼性の高い異常音の発生位置を算出することができる。
なお、上記異常音発生位置の特定方法については、各種の管網110に適用することができる。例えば、水道の配管からの漏水を検出する他、水道以外の各種配管内の漏水を検出する用途、または、工場内の薬液等の配管における薬液等の流体の漏洩を検出する用途などでも使用することができる。
本発明においては、振動センサ200による計測処理が、「計測工程」に相当し、振動センサ200が「計測部」に相当し、ステップS51の処理が「周波数フィルター工程」に相当し、ステップS52の処理が「抽出工程」に相当し、ステップS54の処理が「相互相関関数演算工程」に相当し、ステップS55の処理が「相互相関関数選択工程」に相当し、ステップS51からS55までの処理が、「分析データ作成方法」に相当し、ステップS56の処理が、「フィルター選択工程」に相当し、ステップS57の処理が、「周波数帯域決定工程」に相当し、ステップS51からS57までの処理が、「周波数フィルター作成方法」に相当し、ステップS11からS13、ステップS21からS23、ステップS31からS34、ステップS51からS57までの処理が、「異常音発生位置の特定方法」に相当し、ステップS34の処理が「算出工程」に相当し、演算装置300は、「周波数フィルター部、抽出部、相互相関関数演算部、相互相関関数選択部、フィルター選択部、周波数帯域決定部、算出部」に相当し、異常音の発生特定装置100が「異常音の発生位置の特定装置」に相当する。
本発明の好ましい一実施の形態は上記の通りであるが、本発明はそれだけに制限されない。本発明の精神と範囲から逸脱することのない様々な実施形態が他になされることは理解されよう。さらに、本実施形態において、本発明の構成による作用および効果を述べているが、これら作用および効果は、一例であり、本発明を限定するものではない。
110 管網
200 振動センサ
200 振動センサ
Claims (9)
- 少なくとも2個の振動センサにより振動を計測する計測工程と、
前記計測工程による計測値のコヒーレンス値に基づいて異なる周波数帯域を有する複数の周波数フィルターを作成するフィルター作成工程と、
前記フィルター作成工程で作成されたフィルターの総当たりの組合せを設定する設定工程と、
前記設定工程で設定された組合せのフィルターで計測値を処理する処理工程と、
前記処理工程で処理された処理値から複数の相互相関関数を演算する相互相関関数演算工程と、
前記複数の相互相関関数から所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択工程と、を含む分析データ作成方法。 - 前記周波数フィルター作成工程は、前記コヒーレンス値を所定の周波数幅で周波数帯域を変化させつつ、前記周波数帯域内の前記コヒーレンス値の積分値が大きい順に異なる周波数帯域を複数決定する、請求項1に記載の分析データ作成方法。
- 少なくとも2個の振動センサにより振動を計測する計測工程と、
前記計測工程による計測値のコヒーレンス値に基づいて異なる周波数帯域を有する複数の周波数フィルターを作成するフィルター作成工程と、
前記フィルター作成工程で作成されたフィルターの総当たりの組合せを設定する設定工程と、
前記設定工程で設定された組合せのフィルターで計測値を処理する処理工程と、
前記処理工程で処理された処理値から複数の相互相関関数を演算する相互相関関数演算工程と、
前記複数の相互相関関数から所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択工程と、
前記計測工程から前記相互相関関数選択工程はフィルター選択工程の前に繰り返し複数回行われるものとし、
前記相互相関関数選択工程において選択された頻度に応じて、1または複数の周波数フィルターを決定するフィルター選択工程と、
前記フィルター選択工程により決定された1または複数の周波数フィルターに基づき周波数帯域が決定される周波数帯域決定工程と、を含む、
周波数フィルター作成方法。 - 少なくとも2個の振動センサにより振動を計測する計測工程と、
前記計測工程による計測値のコヒーレンス値に基づいて異なる周波数帯域を有する複数の周波数フィルターを作成するフィルター作成工程と、
前記フィルター作成工程で作成されたフィルターの総当たりの組合せを設定する設定工程と、
前記設定工程で設定された組合せのフィルターで計測値を処理する処理工程と、
前記処理工程で処理された処理値から複数の相互相関関数を演算する相互相関関数演算工程と、
前記複数の相互相関関数から所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択工程と、
前記相互相関関数選択工程により選択された前記相互相関関数に基づき、各振動センサに対する異常音の到達時間差を求め、音の発生位置を算出する算出工程と、を含む、異常音発生位置の特定方法。 - 少なくとも2個の振動センサにより振動を計測する計測工程と、
前記計測工程による計測値のコヒーレンス値に基づいて異なる周波数帯域を有する複数の周波数フィルターを作成するフィルター作成工程と、
前記フィルター作成工程で作成されたフィルターの総当たりの組合せを設定する設定工程と、
前記設定工程で設定された組合せのフィルターで計測値を処理する処理工程と、
前記処理工程で処理された処理値から複数の相互相関関数を演算する相互相関関数演算工程と、
前記複数の相互相関関数から所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択工程と、
前記計測工程から前記相互相関関数選択工程はフィルター選択工程の前に繰り返し複数回行われるものとし、
前記相互相関関数選択工程において選択された頻度に応じて、1または複数の周波数フィルターを決定するフィルター選択工程と、
前記フィルター選択工程により決定された1または複数の周波数フィルターに基づき周波数帯域が決定される周波数帯域決定工程と、
前記周波数帯域決定工程により決定された周波数帯域を通過するフィルターにより計測値を処理し、処理された計測値から相互相関関数を演算し、演算された相互相関関数に基づき、各振動センサに対する異常音の到達時間差を求め、音の発生位置を算出する算出工程と、を含む、異常音発生位置の特定方法。 - 少なくとも2個の振動センサにより振動を計測する計測工程と、
前記計測工程による計測値のコヒーレンス値に基づいて異なる周波数帯域を有する複数の周波数フィルターを作成するフィルター作成工程と、
前記フィルター作成工程で作成されたフィルターの総当たりの組合せを設定する設定工程と、
前記設定工程で設定された組合せのフィルターで計測値を処理する処理工程と、
前記処理工程で処理された処理値から複数の相互相関関数を演算する相互相関関数演算工程と、
前記複数の相互相関関数から所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択工程と、
前記計測工程から前記相互相関関数選択工程は繰り返し複数回行われるものとし、選択された相互相関関数の累積結果に基づき、各振動センサに対する異常音の到達時間差を求め、音の発生位置を算出する算出工程と、を含む、異常音発生位置の特定方法。 - 少なくとも2個の振動センサにより振動を計測する計測部と、
前記計測部による計測値のコヒーレンス値に基づいて異なる周波数帯域を有する複数の周波数フィルターを作成するフィルター作成部と、
前記フィルター作成部で作成されたフィルターの総当たりの組合せを設定する設定部と、
前記設定部で設定された組合せのフィルターで計測値を処理する処理部と、
前記処理部で処理された処理値から複数の相互相関関数を演算する相互相関関数演算部と、
前記複数の相互相関関数から所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択部と、を含む分析データ作成装置。 - 少なくとも2個の振動センサにより振動を計測する計測部と、
前記計測部による計測値のコヒーレンス値に基づいて異なる周波数帯域を有する複数の周波数フィルターを作成するフィルター作成部と、
前記フィルター作成部で作成されたフィルターの総当たりの組合せを設定する設定部と、
前記設定部で設定された組合せのフィルターで計測値を処理する処理部と、
前記処理部で処理された処理値から複数の相互相関関数を演算する相互相関関数演算部と、
前記複数の相互相関関数から所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択部と、
前記計測部による計測から前記相互相関関数選択部による選択はフィルター選択部による決定の前に繰り返し複数回行われるものとし、
前記相互相関関数選択部において選択された頻度に応じて、1または複数の周波数フィルターを決定するフィルター選択部と、
前記フィルター選択部により決定された1または複数の周波数フィルターに基づき周波数帯域が決定される周波数帯域決定部と、を含む、周波数フィルター作成装置。 - 少なくとも2個の振動センサにより振動を計測する計測部と、
前記計測部による計測値のコヒーレンス値に基づいて異なる周波数帯域を有する複数の周波数フィルターを作成するフィルター作成部と、
前記フィルター作成部で作成されたフィルターの総当たりの組合せを設定する設定部と、
前記設定部で設定された組合せのフィルターで計測値を処理する処理部と、
前記処理部で処理された処理値から複数の相互相関関数を演算する相互相関関数演算部と、
前記複数の相互相関関数から所定の相関最大値および所定の波形の少なくとも一方を有する関数を選択する相互相関関数選択部と、
前記相互相関関数選択部により選択された前記相互相関関数に基づき、各振動センサに対する異常音の到達時間差を求め、音の発生位置を算出する算出部と、を含む、異常音発生位置の特定装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2014231581A JP5807107B1 (ja) | 2014-11-14 | 2014-11-14 | 分析データ作成方法、周波数フィルター作成方法、異常音発生位置の特定方法、分析データ作成装置、周波数フィルター作成装置および異常音発生位置の特定装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2014231581A JP5807107B1 (ja) | 2014-11-14 | 2014-11-14 | 分析データ作成方法、周波数フィルター作成方法、異常音発生位置の特定方法、分析データ作成装置、周波数フィルター作成装置および異常音発生位置の特定装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP5807107B1 true JP5807107B1 (ja) | 2015-11-10 |
JP2016095231A JP2016095231A (ja) | 2016-05-26 |
Family
ID=54545761
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2014231581A Active JP5807107B1 (ja) | 2014-11-14 | 2014-11-14 | 分析データ作成方法、周波数フィルター作成方法、異常音発生位置の特定方法、分析データ作成装置、周波数フィルター作成装置および異常音発生位置の特定装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5807107B1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115439990A (zh) * | 2022-08-29 | 2022-12-06 | 中国联合网络通信集团有限公司 | 人孔异常确定方法、装置、存储介质和人孔 |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019054433A1 (ja) * | 2017-09-14 | 2019-03-21 | 日本電気株式会社 | 診断装置、診断方法及びプログラム記憶媒体 |
KR102038689B1 (ko) * | 2018-06-14 | 2019-10-30 | 한국원자력연구원 | 거리차-주파수 분석을 이용한 배관의 누설 감지장치 및 방법 |
JP7304301B2 (ja) * | 2020-02-04 | 2023-07-06 | 株式会社日立製作所 | 音響診断方法、音響診断システム、及び音響診断プログラム |
-
2014
- 2014-11-14 JP JP2014231581A patent/JP5807107B1/ja active Active
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115439990A (zh) * | 2022-08-29 | 2022-12-06 | 中国联合网络通信集团有限公司 | 人孔异常确定方法、装置、存储介质和人孔 |
CN115439990B (zh) * | 2022-08-29 | 2023-09-19 | 中国联合网络通信集团有限公司 | 人孔异常确定方法、装置、存储介质和人孔 |
Also Published As
Publication number | Publication date |
---|---|
JP2016095231A (ja) | 2016-05-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2020261087B2 (en) | Method and system for detecting a structural anomaly in a pipeline network | |
Pal et al. | Detecting & locating leaks in water distribution polyethylene pipes | |
JP5807107B1 (ja) | 分析データ作成方法、周波数フィルター作成方法、異常音発生位置の特定方法、分析データ作成装置、周波数フィルター作成装置および異常音発生位置の特定装置 | |
Zan et al. | Event detection and localization in urban water distribution network | |
JP2005308511A5 (ja) | ||
GB2421311A (en) | Assessing the size of a leak in a pipeline by detecting leak noise and pressure | |
EP2326933B1 (en) | Improvements in and relating to apparatus for the airborne acoustic inspection of pipes | |
JP2008051776A (ja) | 漏水検出装置および漏水検出方法 | |
WO2015072130A1 (ja) | 漏洩判定システムおよび漏洩判定方法 | |
CN116577037B (zh) | 一种基于非均匀频率语谱图的风管泄漏信号检测方法 | |
JP6789042B2 (ja) | 漏水場所の特定方法 | |
WO2014050358A1 (ja) | 漏洩検査装置、漏洩検査方法、及び漏洩検査プログラム | |
JP2017083291A (ja) | 異常音の発生位置特定方法および異常音の発生位置特定装置 | |
JP6557576B2 (ja) | 異常音の発生位置特定方法および異常音の発生位置特定装置 | |
JP6316131B2 (ja) | 異常音発生位置の特定方法 | |
US10948376B2 (en) | Apparatus and method of detecting leak sound in plant equipment using time-frequency transformation | |
Bakhti et al. | Experimental validation of hybrid EMD-correlation acoustic digital leaks detector in water distribution network system | |
JP6408929B2 (ja) | 分析データ作成方法、漏水位置検知装置および漏水位置特定方法 | |
JP2018205192A (ja) | 漏水位置特定方法、漏洩有無判定方法 | |
KR20180110543A (ko) | 누출음의 시공간 특성을 이용한 미세누출 탐지 장치 및 이를 이용한 미세누출 탐지 방법 | |
KR102061616B1 (ko) | 비선형 스펙트럼 상관관계를 이용한 구조물 진단 방법 및 구조물 진단 시스템 | |
JP2017173034A (ja) | 管路の異常箇所特定方法 | |
JP3639160B2 (ja) | 異常箇所検出装置 | |
Shi et al. | Leak detection in virtually isolated pipe sections within a complex pipe system using a two-source-four-sensor transient testing configuration | |
JP2014219342A (ja) | 埋設管路の漏洩検出方法および装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20150812 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20150907 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 5807107 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |