JP3488623B2 - 漏水位置検出方法および漏水位置検出装置 - Google Patents
漏水位置検出方法および漏水位置検出装置Info
- Publication number
- JP3488623B2 JP3488623B2 JP05200098A JP5200098A JP3488623B2 JP 3488623 B2 JP3488623 B2 JP 3488623B2 JP 05200098 A JP05200098 A JP 05200098A JP 5200098 A JP5200098 A JP 5200098A JP 3488623 B2 JP3488623 B2 JP 3488623B2
- Authority
- JP
- Japan
- Prior art keywords
- leak
- intensity
- ground
- vector
- water
- 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.)
- Expired - Fee Related
Links
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
Landscapes
- Measuring Volume Flow (AREA)
- Examining Or Testing Airtightness (AREA)
Description
【０００１】
【発明の属する技術分野】本発明は、地中に埋設されて
いる水道管等の漏水位置を地表において推定検出する漏
水位置検出方法および漏水位置検出装置に関する。
いる水道管等の漏水位置を地表において推定検出する漏
水位置検出方法および漏水位置検出装置に関する。
【０００２】
【従来の技術】老朽化等によって地中に埋設されている
水道管が破損すると、この破損箇所からいわるゆる漏水
が発生する。水道管からの漏水によって地中に浸透する
水の量は莫大なものであるといわれている。最近のよう
に水不足が重大な社会問題となつているときには、水道
管からの漏水をできるだけ少なくすることが水資源確保
のためにも重要である。
水道管が破損すると、この破損箇所からいわるゆる漏水
が発生する。水道管からの漏水によって地中に浸透する
水の量は莫大なものであるといわれている。最近のよう
に水不足が重大な社会問題となつているときには、水道
管からの漏水をできるだけ少なくすることが水資源確保
のためにも重要である。
【０００３】ところで、地中に埋設されている水道管の
漏水位置を検知する方法としては、従来から幾つか考え
られている。
漏水位置を検知する方法としては、従来から幾つか考え
られている。
【０００４】代表的な例は人間の聴覚による方法であ
る。この方法では、水道管の地上に露出している部分、
たとえば弁の設けられている場所において、熟練作業者
が水道管を通して伝わる音を聴くことにより、その水道
管に漏水音（振動）が発生しているか否かを検知する。
この検知作業は、人間の聴覚および経験によるよるた
め、車の振動等による外部雑音の少ない深夜等に行われ
るのが一般的である。そして、水道管から漏水音が発生
していると検知されたときには、漏水音の最も大きい位
置を把握するために、作業者が地上を移動しながらマイ
ク（音聴センサ）や音聴棒を用いて音源を探査する。漏
水音の最も大きい位置を把握した時点で、その地面を掘
って水道管を露出させ、現実に漏水しているか否かを調
べ、必要な修復を行う。
る。この方法では、水道管の地上に露出している部分、
たとえば弁の設けられている場所において、熟練作業者
が水道管を通して伝わる音を聴くことにより、その水道
管に漏水音（振動）が発生しているか否かを検知する。
この検知作業は、人間の聴覚および経験によるよるた
め、車の振動等による外部雑音の少ない深夜等に行われ
るのが一般的である。そして、水道管から漏水音が発生
していると検知されたときには、漏水音の最も大きい位
置を把握するために、作業者が地上を移動しながらマイ
ク（音聴センサ）や音聴棒を用いて音源を探査する。漏
水音の最も大きい位置を把握した時点で、その地面を掘
って水道管を露出させ、現実に漏水しているか否かを調
べ、必要な修復を行う。
【０００５】しかし、この検出方法では、時間的な制約
があるばかりか、作業者の感覚を主体にして行われるた
め、作業能率が非常に悪く、しかも漏水箇所を精度良く
見つけ出し難いという問題があった。
があるばかりか、作業者の感覚を主体にして行われるた
め、作業能率が非常に悪く、しかも漏水箇所を精度良く
見つけ出し難いという問題があった。
【０００６】一方、最近では、音響信号処理技術を応用
し、地表面に伝わる漏水エネルギの流れ方向を検出する
ことで、人間の感覚に頼ることなく、地中に存在する漏
水箇所を地上から検出する方法が提案されている。
し、地表面に伝わる漏水エネルギの流れ方向を検出する
ことで、人間の感覚に頼ることなく、地中に存在する漏
水箇所を地上から検出する方法が提案されている。
【０００７】この方法では、図１８に示すように、地中
に埋設されている水道管１の一部分に漏水位置Ｅからの
漏水音を検知する漏水音検知センサ２を取付けるととも
に地表面３の２箇所に所定間隔ｄだけ離間させて地表振
動検知センサ４，５を配置し、この地表振動検知センサ
４，５を水道管１の配設経路に沿って移動させ、このと
きに信号処理装置６で地表振動検知センサ４，５の出力
信号から漏水音検知センサ２で得られた信号成分に関わ
りのある信号成分を相互相関処理で抽出し、抽出された
信号成分から地表面に伝わる漏水エネルギのインテンシ
ティを演算し、このインテンシティの大きさおよび極性
から漏水位置Ｅを示すベクトルを得るようにしている。
に埋設されている水道管１の一部分に漏水位置Ｅからの
漏水音を検知する漏水音検知センサ２を取付けるととも
に地表面３の２箇所に所定間隔ｄだけ離間させて地表振
動検知センサ４，５を配置し、この地表振動検知センサ
４，５を水道管１の配設経路に沿って移動させ、このと
きに信号処理装置６で地表振動検知センサ４，５の出力
信号から漏水音検知センサ２で得られた信号成分に関わ
りのある信号成分を相互相関処理で抽出し、抽出された
信号成分から地表面に伝わる漏水エネルギのインテンシ
ティを演算し、このインテンシティの大きさおよび極性
から漏水位置Ｅを示すベクトルを得るようにしている。
【０００８】すなわち、この方法では、漏水位置Ｅから
放射された漏水音が地中を伝播して地表振動検知センサ
４，５に入射するとき、地表振動検知センサ４，５の位
置の違いに起因して位相差が生じることに着目し、この
位相差をインテンシティとして捉え、このインテンシテ
ィの大きさおよび極性から漏水位置Ｅを示すベクトル
（大きさおよび方向）を得るようにしている。
放射された漏水音が地中を伝播して地表振動検知センサ
４，５に入射するとき、地表振動検知センサ４，５の位
置の違いに起因して位相差が生じることに着目し、この
位相差をインテンシティとして捉え、このインテンシテ
ィの大きさおよび極性から漏水位置Ｅを示すベクトル
（大きさおよび方向）を得るようにしている。
【０００９】しかしながら、このような音響信号処理技
術を応用した検出方法にあっても次のような問題あっ
た。
術を応用した検出方法にあっても次のような問題あっ
た。
【００１０】すなわち、図１８に示す方法で用いるイン
テンシティは、音響分野でいうアクティブインテンシテ
ィであり（以下、これを用いる方法をＡＩ法と称す）、
音の進行方向に垂直な単位面積を単位時間に通過する音
のエネルギで定義される。これは音圧（振動）と粒子速
度の時間平均の量であり、粒子速度項は２つの地表振動
検知センサ４，５で検出した音圧（振動）の差を積分す
ることで得られる。この手法は直接積分法と呼ばれ、空
間域での音源探査では広く用いられている。
テンシティは、音響分野でいうアクティブインテンシテ
ィであり（以下、これを用いる方法をＡＩ法と称す）、
音の進行方向に垂直な単位面積を単位時間に通過する音
のエネルギで定義される。これは音圧（振動）と粒子速
度の時間平均の量であり、粒子速度項は２つの地表振動
検知センサ４，５で検出した音圧（振動）の差を積分す
ることで得られる。この手法は直接積分法と呼ばれ、空
間域での音源探査では広く用いられている。
【００１１】このＡＩ法を用いた装置で得られるインテ
ンシティＩａ（ｘ，ｙ）は、図１８に示される条件を用
いて示すと、（１）式で評価される。なお、漏水音は点
音源から放射され、地中の伝播経路では減衰はあるもの
の、反射はないものとする。
ンシティＩａ（ｘ，ｙ）は、図１８に示される条件を用
いて示すと、（１）式で評価される。なお、漏水音は点
音源から放射され、地中の伝播経路では減衰はあるもの
の、反射はないものとする。
【数１】
【００１２】この（１）式から判るように、ＡＩ法にお
けるインテンシティＩａ（ｘ，ｙ）は、漏水位置Ｅから
地表振動検知センサ４，５までの伝播距離の差とセンサ
設置誤差αとによって大きく変化する。このうちセンサ
設置誤差αとは、地表面と各地表振動センサとの接触状
態から生じる漏水振動検出時間遅れのことである。
けるインテンシティＩａ（ｘ，ｙ）は、漏水位置Ｅから
地表振動検知センサ４，５までの伝播距離の差とセンサ
設置誤差αとによって大きく変化する。このうちセンサ
設置誤差αとは、地表面と各地表振動センサとの接触状
態から生じる漏水振動検出時間遅れのことである。
【００１３】今、センサ設置誤差αが零であると仮定す
ると、２つの地表振動検知センサ４，５を、その間隔ｄ
を保って水道管１に沿って移動させたときの測定位置ｘ
（ここでは漏水位置Ｅから２つの地表振動検知センサ
４，５の中間位置までの水平距離とする）に対するイン
テンシティＩａ（ｘ，ｙ）の大きさ（振動振幅）および
極性は、図１９（ｂ）に示すように、漏水位置Ｅに近付
くにしたがってインテンシティＩａ（ｘ，ｙ）が次第に
大きくなり、さらに近付くと一転して小さくなり、漏水
位置Ｅの真上では零となる。そして、２つの地表振動検
知センサ４，５をさらに移動させると、インテンシティ
Ｉａ（ｘ，ｙ）の極性が反転し、その大きさが次第に増
加する。このときに得られる漏水位置特定ベクトルは図
１９（ａ）に示すように変化する。
ると、２つの地表振動検知センサ４，５を、その間隔ｄ
を保って水道管１に沿って移動させたときの測定位置ｘ
（ここでは漏水位置Ｅから２つの地表振動検知センサ
４，５の中間位置までの水平距離とする）に対するイン
テンシティＩａ（ｘ，ｙ）の大きさ（振動振幅）および
極性は、図１９（ｂ）に示すように、漏水位置Ｅに近付
くにしたがってインテンシティＩａ（ｘ，ｙ）が次第に
大きくなり、さらに近付くと一転して小さくなり、漏水
位置Ｅの真上では零となる。そして、２つの地表振動検
知センサ４，５をさらに移動させると、インテンシティ
Ｉａ（ｘ，ｙ）の極性が反転し、その大きさが次第に増
加する。このときに得られる漏水位置特定ベクトルは図
１９（ａ）に示すように変化する。
【００１４】したがって、２つの地表振動検知センサ
４，５を水道管１に沿って移動させ、インテンシティＩ
ａ（ｘ，ｙ）の極性が＋（−）から−（＋）に転じる零
点位置、つまり漏水位置特定ベクトルが零の位置を求め
れば、その位置が漏水位置Ｅの真上（あるいは真上の
横）であることから、漏水位置Ｅを知るとができる。な
お、漏水位置Ｅから離れた場所では、漏水位置特定ベク
トルが反転するが、振幅が小さくなることから検知精度
には影響を与えない。
４，５を水道管１に沿って移動させ、インテンシティＩ
ａ（ｘ，ｙ）の極性が＋（−）から−（＋）に転じる零
点位置、つまり漏水位置特定ベクトルが零の位置を求め
れば、その位置が漏水位置Ｅの真上（あるいは真上の
横）であることから、漏水位置Ｅを知るとができる。な
お、漏水位置Ｅから離れた場所では、漏水位置特定ベク
トルが反転するが、振幅が小さくなることから検知精度
には影響を与えない。
【００１５】このように、図１８に示されるＡＩ法で
は、条件が理想的に満たされた場合でも、そのインテン
シティＩａ（ｘ，ｙ）が漏水位置Ｅの真上で零となって
しまうため、漏水位置特定ベクトルも零となる。漏水位
置近傍では、本来漏水から受けるエネルギが大きいにも
拘らず、この有効な情報をＡＩ法では漏水位置特定ベク
トルに反映できないという問題があった。
は、条件が理想的に満たされた場合でも、そのインテン
シティＩａ（ｘ，ｙ）が漏水位置Ｅの真上で零となって
しまうため、漏水位置特定ベクトルも零となる。漏水位
置近傍では、本来漏水から受けるエネルギが大きいにも
拘らず、この有効な情報をＡＩ法では漏水位置特定ベク
トルに反映できないという問題があった。
【００１６】なお、通常の空間域における音源探査で利
用するインテンシティも上述したアクティブインテンシ
ティの原理を使っている。しかし、この場合には漏水位
置の真上で見られるような上述した問題は生じない。こ
れは、音源探査の場合には空間に音源信号検出用マイク
を互いに直交するｘｙｚ軸上に複数設置することがで
き、１つの軸上のインテンシティ射影成分がたとえ零に
なつても、他の軸上の射影成分が大きさを持つことか
ら、最終的にインテンシティは３軸のベクトルの合成と
なり、エネルギ量を少なく算出することがないことによ
る。
用するインテンシティも上述したアクティブインテンシ
ティの原理を使っている。しかし、この場合には漏水位
置の真上で見られるような上述した問題は生じない。こ
れは、音源探査の場合には空間に音源信号検出用マイク
を互いに直交するｘｙｚ軸上に複数設置することがで
き、１つの軸上のインテンシティ射影成分がたとえ零に
なつても、他の軸上の射影成分が大きさを持つことか
ら、最終的にインテンシティは３軸のベクトルの合成と
なり、エネルギ量を少なく算出することがないことによ
る。
【００１７】漏水位置検出の場合においても、もしも地
中にセンサを埋設できれば、深さ方向に伝わるエネルギ
量を検出できるので、漏水位置真上の地表面インテンシ
ティが零になっても、深さ（鉛直）方向のインテンシテ
ィが最大となり、通常の空間域における音源探査と同様
の探査機能を発揮させることができる。しかし、漏水位
置検出では現実問題として地下にセンサを埋設すること
はできないので、地下方向のインテンシティを漏水位置
特定ベクトルに加えることができず、漏水エネルギを最
大限に利用することはできない。
中にセンサを埋設できれば、深さ方向に伝わるエネルギ
量を検出できるので、漏水位置真上の地表面インテンシ
ティが零になっても、深さ（鉛直）方向のインテンシテ
ィが最大となり、通常の空間域における音源探査と同様
の探査機能を発揮させることができる。しかし、漏水位
置検出では現実問題として地下にセンサを埋設すること
はできないので、地下方向のインテンシティを漏水位置
特定ベクトルに加えることができず、漏水エネルギを最
大限に利用することはできない。
【００１８】
【発明が解決しようとする課題】上述の如く、音響信号
処理技術を応用した漏水位置検出方法は、人間の聴覚に
頼る方法に比べて自由度が大きい。しかし、アクティブ
インテンシティ法を用いた従来の漏水位置検出方法で
は、漏水エネルギの有効利用を図ることができないた
め、この面から精度の高い漏水位置検出ができないとい
う問題があった。
処理技術を応用した漏水位置検出方法は、人間の聴覚に
頼る方法に比べて自由度が大きい。しかし、アクティブ
インテンシティ法を用いた従来の漏水位置検出方法で
は、漏水エネルギの有効利用を図ることができないた
め、この面から精度の高い漏水位置検出ができないとい
う問題があった。
【００１９】そこで本発明は、上述した従来のアクティ
ブインテンシティ法が持つ不具合を解消でき、装置構成
も単純化でき、検出精度の向上に寄与できる漏水位置検
出方法および漏水位置検出装置を提供することを目的と
している。
ブインテンシティ法が持つ不具合を解消でき、装置構成
も単純化でき、検出精度の向上に寄与できる漏水位置検
出方法および漏水位置検出装置を提供することを目的と
している。
【００２０】
【課題を解決するための手段】上記目的を達成するため
に、請求項１に係る漏水位置検出方法では、地中に埋設
されている配管の漏水位置を推定検出するに当り、上記
配管の一部分に漏水音検知センサを取付けるとともに少
なくとも一対の地表振動検知センサを地表面に移動可能
に配置し、これら地表振動検知センサの出力信号から前
記漏水音検知センサで得られた信号とは関わりのない雑
音信号を相互相関処理でそれぞれ除去し、この処理を経
た対をなす信号を用いて地表面に伝わる漏水エネルギの
アクティブインテンシティを演算し、さらに前記相互相
関処理を経た一方の信号と時間積分処理された他方の信
号とを用いて周波数重み付けのあるリアクティブインテ
ンシテイを演算し、漏水の特定位置を示すベクトルの大
きさには前記リアクティブインテンシティの大きさを用
い、上記ベクトルの方向には前記アクティブインテンシ
ティの極性を用いて漏水位置特定ベクトルを求めるよう
にしている。
に、請求項１に係る漏水位置検出方法では、地中に埋設
されている配管の漏水位置を推定検出するに当り、上記
配管の一部分に漏水音検知センサを取付けるとともに少
なくとも一対の地表振動検知センサを地表面に移動可能
に配置し、これら地表振動検知センサの出力信号から前
記漏水音検知センサで得られた信号とは関わりのない雑
音信号を相互相関処理でそれぞれ除去し、この処理を経
た対をなす信号を用いて地表面に伝わる漏水エネルギの
アクティブインテンシティを演算し、さらに前記相互相
関処理を経た一方の信号と時間積分処理された他方の信
号とを用いて周波数重み付けのあるリアクティブインテ
ンシテイを演算し、漏水の特定位置を示すベクトルの大
きさには前記リアクティブインテンシティの大きさを用
い、上記ベクトルの方向には前記アクティブインテンシ
ティの極性を用いて漏水位置特定ベクトルを求めるよう
にしている。
【００２１】また、上記目的を達成するために、請求項
２に係る漏水位置検出方法では、地中に埋設されている
配管の漏水位置を推定検出するに当り、上記配管の一部
分に漏水音検知センサを取付けるとともに少なくとも一
対の地表振動検知センサを地表面に移動可能に配置し、
これら地表振動検知センサの出力信号から前記漏水音検
知センサで得られた信号とは関わりのない雑音信号を相
互相関処理でそれぞれ除去し、この処理を経た対をなす
信号を用いて地表面に伝わる漏水エネルギのアクティブ
インテンシティを演算し、さらに前記相互相関処理を経
た一方の信号と時間積分処理または時間微分処理された
他方の信号とを用いて最終的に周波数重み付けの除かれ
た真のリアクティブインテンシテイを演算し、漏水の特
定位置を示すベクトルの大きさには前記真のリアクティ
ブインテンシティの大きさを用い、上記ベクトルの方向
には前記アクティブインテンシティの極性を用いて漏水
位置特定ベクトルを求めるようにしている。
２に係る漏水位置検出方法では、地中に埋設されている
配管の漏水位置を推定検出するに当り、上記配管の一部
分に漏水音検知センサを取付けるとともに少なくとも一
対の地表振動検知センサを地表面に移動可能に配置し、
これら地表振動検知センサの出力信号から前記漏水音検
知センサで得られた信号とは関わりのない雑音信号を相
互相関処理でそれぞれ除去し、この処理を経た対をなす
信号を用いて地表面に伝わる漏水エネルギのアクティブ
インテンシティを演算し、さらに前記相互相関処理を経
た一方の信号と時間積分処理または時間微分処理された
他方の信号とを用いて最終的に周波数重み付けの除かれ
た真のリアクティブインテンシテイを演算し、漏水の特
定位置を示すベクトルの大きさには前記真のリアクティ
ブインテンシティの大きさを用い、上記ベクトルの方向
には前記アクティブインテンシティの極性を用いて漏水
位置特定ベクトルを求めるようにしている。
【００２２】また、上記目的を達成するために、請求項
５に係る漏水位置検出装置では、地中に埋設されている
配管の一部分に取付けられる漏水音検知センサと、地表
面に移動可能に配置される少なくとも一対の地表振動検
知センサと、これら地表振動検知センサの出力信号から
前記漏水音検知センサで得られた信号とは関わりのない
雑音信号を相互相関処理でそれぞれ除去する手段と、こ
の手段による処理を経た対をなす信号を用いて地表面に
伝わる漏水エネルギのアクティブインテンシティを演算
する手段と、前記相互相関処理を経た一方の信号と時間
積分処理された他方の信号とを用いて周波数重み付けの
あるリアクティブインテンシティを演算する手段と、漏
水の特定位置を示すベクトルの大きさには前記リアクテ
ィブインテンシティの大きさを用い、上記ベクトルの方
向には前記アクティブインテンシティの極性を用いて漏
水位置特定ベクトルを求める手段とを備えている。
５に係る漏水位置検出装置では、地中に埋設されている
配管の一部分に取付けられる漏水音検知センサと、地表
面に移動可能に配置される少なくとも一対の地表振動検
知センサと、これら地表振動検知センサの出力信号から
前記漏水音検知センサで得られた信号とは関わりのない
雑音信号を相互相関処理でそれぞれ除去する手段と、こ
の手段による処理を経た対をなす信号を用いて地表面に
伝わる漏水エネルギのアクティブインテンシティを演算
する手段と、前記相互相関処理を経た一方の信号と時間
積分処理された他方の信号とを用いて周波数重み付けの
あるリアクティブインテンシティを演算する手段と、漏
水の特定位置を示すベクトルの大きさには前記リアクテ
ィブインテンシティの大きさを用い、上記ベクトルの方
向には前記アクティブインテンシティの極性を用いて漏
水位置特定ベクトルを求める手段とを備えている。
【００２３】また、上記目的を達成するために、請求項
６に係る漏水位置検出装置では、地中に埋設されている
配管の一部分に取付けられる漏水音検知センサと、地表
面に移動可能に配置される少なくとも一対の地表振動検
知センサと、これら地表振動検知センサの出力信号から
前記漏水音検知センサで得られた信号とは関わりのない
雑音信号を相互相関処理でそれぞれ除去する手段と、こ
の手段による処理を経た対をなす信号を用いて地表面に
伝わる漏水エネルギのアクティブインテンシティを演算
する手段と、前記相互相関処理を経た一方の信号と時間
積分処理または時間微分処理された他方の信号とを用い
て最終的に周波数重み付けの除かれた真のリアクティブ
インテンシティを演算する手段と、漏水の特定位置を示
すベクトルの大きさには前記真のリアクティブインテン
シティの大きさを用い、上記ベクトルの方向には前記ア
クティブインテンシティの極性を用いて漏水位置特定ベ
クトルを求める手段とを備えている。
６に係る漏水位置検出装置では、地中に埋設されている
配管の一部分に取付けられる漏水音検知センサと、地表
面に移動可能に配置される少なくとも一対の地表振動検
知センサと、これら地表振動検知センサの出力信号から
前記漏水音検知センサで得られた信号とは関わりのない
雑音信号を相互相関処理でそれぞれ除去する手段と、こ
の手段による処理を経た対をなす信号を用いて地表面に
伝わる漏水エネルギのアクティブインテンシティを演算
する手段と、前記相互相関処理を経た一方の信号と時間
積分処理または時間微分処理された他方の信号とを用い
て最終的に周波数重み付けの除かれた真のリアクティブ
インテンシティを演算する手段と、漏水の特定位置を示
すベクトルの大きさには前記真のリアクティブインテン
シティの大きさを用い、上記ベクトルの方向には前記ア
クティブインテンシティの極性を用いて漏水位置特定ベ
クトルを求める手段とを備えている。
【００２４】なお、請求項１，２に係る漏水位置検出方
法および請求項５，６に係る漏水位置検出装置におい
て、漏水音に含まれている周波数成分に対応させて間隔
を異ならせた実質的に複数対の地表振動検知センサを移
動可能に配置し、対象とする周波数成分毎に前記漏水位
置特定ベクトルをそれぞれ求めることが好ましい。この
場合、複数の周波数成分に対応させて配置される実質的
に複数対の地表振動センサは、対象とする周波数が低い
ほど間隔が広く設定される。そして、各周波数成分毎に
求められたインテンシティを合成して漏水位置特定ベク
トルを求めるようにしてもよい。
法および請求項５，６に係る漏水位置検出装置におい
て、漏水音に含まれている周波数成分に対応させて間隔
を異ならせた実質的に複数対の地表振動検知センサを移
動可能に配置し、対象とする周波数成分毎に前記漏水位
置特定ベクトルをそれぞれ求めることが好ましい。この
場合、複数の周波数成分に対応させて配置される実質的
に複数対の地表振動センサは、対象とする周波数が低い
ほど間隔が広く設定される。そして、各周波数成分毎に
求められたインテンシティを合成して漏水位置特定ベク
トルを求めるようにしてもよい。
【００２５】請求項１，２に係る漏水位置検出方法およ
び請求項５，６に係る漏水位置検出装置では、相互相関
処理を経た対をなす信号を用いて地表面に伝わる漏水エ
ネルギのアクティブインテンシティと、音圧（振動）の
変化のみでエネルギの伝播はないリアクティブインテン
シティ（電気系や機械系での無効パワーに相当する）と
を求め、漏水の特定位置を示すベクトルの大きさ（振動
振幅）にはリアクティブインテンシティの大きさ（絶対
値）を用い、上記ベクトルの方向にはアクティブインテ
ンシティの極性を用いて漏水位置特定ベクトルを求めて
いる。
び請求項５，６に係る漏水位置検出装置では、相互相関
処理を経た対をなす信号を用いて地表面に伝わる漏水エ
ネルギのアクティブインテンシティと、音圧（振動）の
変化のみでエネルギの伝播はないリアクティブインテン
シティ（電気系や機械系での無効パワーに相当する）と
を求め、漏水の特定位置を示すベクトルの大きさ（振動
振幅）にはリアクティブインテンシティの大きさ（絶対
値）を用い、上記ベクトルの方向にはアクティブインテ
ンシティの極性を用いて漏水位置特定ベクトルを求めて
いる。
【００２６】このようにアクティブインテンシティとリ
アクティブインテンシティとを用いて漏水位置特定ベク
トルを算出する方法を以後、ＷＬＩ法（Ｗａｔｅｒ−Ｌ
ｅａｋ−Ｉｎｔｅｎｓｉｔｙ）と呼ぶことにする。な
お、リアクティブインテンシティだけを用いて漏水位置
特定ベクトルを算出する方法も考えられる。これをＲＩ
法と呼ぶことにする。ＷＬＩ法は、ＡＩ法とＲＩ法とを
組み合わせたものといえる。
アクティブインテンシティとを用いて漏水位置特定ベク
トルを算出する方法を以後、ＷＬＩ法（Ｗａｔｅｒ−Ｌ
ｅａｋ−Ｉｎｔｅｎｓｉｔｙ）と呼ぶことにする。な
お、リアクティブインテンシティだけを用いて漏水位置
特定ベクトルを算出する方法も考えられる。これをＲＩ
法と呼ぶことにする。ＷＬＩ法は、ＡＩ法とＲＩ法とを
組み合わせたものといえる。
【００２７】ＷＬＩ法で最終的に求められるインテンシ
ティＩ（ｘ，ｙ）は、（１）式で用いている変数を用い
て示すと、（２）式のようになる。
ティＩ（ｘ，ｙ）は、（１）式で用いている変数を用い
て示すと、（２）式のようになる。
【００２８】
【数２】
【００２９】この（２）式は、後述するように、本来、
漏水エネルギが最も大きい位置である漏水位置真上にお
ける漏水位置特定ベクトルの大きさを最大にできること
を意味している。したがって、ＷＬＩ法によれば、より
精度の高い漏水位置特定が可能となる。
漏水エネルギが最も大きい位置である漏水位置真上にお
ける漏水位置特定ベクトルの大きさを最大にできること
を意味している。したがって、ＷＬＩ法によれば、より
精度の高い漏水位置特定が可能となる。
【００３０】また、特に請求項１，５に係る漏水位置検
出方法および漏水位置検出装置では、リアクティブイン
テンシティの演算に、相互相関処理を経た一方の信号と
時間積分処理された他方の信号とを用いており、請求項
２，６に係る漏水位置検出方法および漏水位置検出装置
では、相互相関処理を経た一方の信号と時間積分処理ま
たは時間微分処理された他方の信号とを用いて最終的に
周波数重み付けの除かれた真のリアクティブインテンシ
ティを演算するようにしている。
出方法および漏水位置検出装置では、リアクティブイン
テンシティの演算に、相互相関処理を経た一方の信号と
時間積分処理された他方の信号とを用いており、請求項
２，６に係る漏水位置検出方法および漏水位置検出装置
では、相互相関処理を経た一方の信号と時間積分処理ま
たは時間微分処理された他方の信号とを用いて最終的に
周波数重み付けの除かれた真のリアクティブインテンシ
ティを演算するようにしている。
【００３１】このように、リアクティブインテンシティ
の演算過程に時間積分処理といった簡単なソフトウェア
で実現できる手法を取り入れているので、演算系にバン
ドパスフィルタ等を用いる場合に比べて演算系の単純化
を図ることができる。
の演算過程に時間積分処理といった簡単なソフトウェア
で実現できる手法を取り入れているので、演算系にバン
ドパスフィルタ等を用いる場合に比べて演算系の単純化
を図ることができる。
【００３２】すなわち、アクティブインテンシティＩａ
（ｘ，ｙ）は、（１）式で示すとおりであり、これに対
して真のリアクティブインテンシティＩｒ（ｘ，ｙ）は
（３）式で示される。
（ｘ，ｙ）は、（１）式で示すとおりであり、これに対
して真のリアクティブインテンシティＩｒ（ｘ，ｙ）は
（３）式で示される。
【００３３】
【数３】
【００３４】実際にリアクティブインテンシティＩｒ
（ｘ，ｙ）を算出するには、相互相関処理を経た２つの
信号のうちの一方の信号の位相を９０度ずらす必要があ
る。漏水音は単一周波数ではないので、対象とする周波
数全域に亘って位相を９０度ずらす必要がある。これを
実現する手段として、複数のバンドパスフィルタを用い
るとともに、それぞれのバンドの中心周波数の位相を９
０度ずらすことが考えられるが、この手段では数多くの
フィルタを必要とするので、装置を構成したときに大掛
かりになるばかりか、演算結果に誤差が入り込み易い。
しかし、相互相関処理を経た一方の信号と時間積分処理
された他方の信号とを用いてリアクティブインテンシテ
ィＩｒ（ｘ，ｙ）を算出する方法であると、上述した問
題を解決でき、ＷＬＩ法の利点を最大限に発揮させた状
態で、なおかつ漏水位置特定精度を向上させることがで
きる。
（ｘ，ｙ）を算出するには、相互相関処理を経た２つの
信号のうちの一方の信号の位相を９０度ずらす必要があ
る。漏水音は単一周波数ではないので、対象とする周波
数全域に亘って位相を９０度ずらす必要がある。これを
実現する手段として、複数のバンドパスフィルタを用い
るとともに、それぞれのバンドの中心周波数の位相を９
０度ずらすことが考えられるが、この手段では数多くの
フィルタを必要とするので、装置を構成したときに大掛
かりになるばかりか、演算結果に誤差が入り込み易い。
しかし、相互相関処理を経た一方の信号と時間積分処理
された他方の信号とを用いてリアクティブインテンシテ
ィＩｒ（ｘ，ｙ）を算出する方法であると、上述した問
題を解決でき、ＷＬＩ法の利点を最大限に発揮させた状
態で、なおかつ漏水位置特定精度を向上させることがで
きる。
【００３５】なお、相互相関処理を経た一方の信号と時
間積分処理された他方の信号とを用いて算出されたリア
クティブインテンシティＩｒ（ｘ，ｙ）は、アクティブ
インテンシティＩａ（ｘ，ｙ）に対して９０度位相がず
れた真のリアクティブインテンシティの振幅に周波数の
重み付け、具体的には１／（２πｆ）の重みを付けたも
のとなるが、地中を伝播する音波は周波数が高くなるほ
ど距離減衰作用が大きく、実際にはセンサまで到達しな
いことが多いので、重み付けの影響は極めて少ない。勿
論、請求項２および６に係る漏水位置検出方法および漏
水位置検出装置のように、最終的に重み付けを除く処理
を施してもよい。
間積分処理された他方の信号とを用いて算出されたリア
クティブインテンシティＩｒ（ｘ，ｙ）は、アクティブ
インテンシティＩａ（ｘ，ｙ）に対して９０度位相がず
れた真のリアクティブインテンシティの振幅に周波数の
重み付け、具体的には１／（２πｆ）の重みを付けたも
のとなるが、地中を伝播する音波は周波数が高くなるほ
ど距離減衰作用が大きく、実際にはセンサまで到達しな
いことが多いので、重み付けの影響は極めて少ない。勿
論、請求項２および６に係る漏水位置検出方法および漏
水位置検出装置のように、最終的に重み付けを除く処理
を施してもよい。
【００３６】
【発明の実施の形態】以下、図面を参照しながら発明の
実施形態を説明する。
実施形態を説明する。
【００３７】図１には本発明の一実施形態に係る漏水位
置検出装置のブロック構成図が示されている。
置検出装置のブロック構成図が示されている。
【００３８】この装置は、大きく分けて、地中に埋設さ
れている水道管１１の一部分、たとえば地上に露出して
いる弁部分等に取付けられる漏水音検知センサ１２と、
漏水音信号に含まれる複数の周波数帯域を検出対象にし
て地表面１３に移動可能に配置され、対象とする周波数
帯域に対応させて間隔をそれぞれ異ならせた実質的に複
数対の地表振動検知センサ１４〜１８と、これら地表振
動検知センサ１４〜１８の出力信号から漏水音検知セン
サ１２で得られた信号とは関わりのない雑音信号を相互
相関処理でそれぞれ除去する相互相関処理部１９と、こ
の処理部１９を経た各対をなす信号を用いて地表面１３
に伝わる漏水エネルギのうちの上記信号を出力した対を
なす地表振動検知センサが対象にしている周波数成分を
それぞれ抽出処理する信号処理部２０と、アクティブイ
ンテンシティ算出部２１とリアクティブインテンシティ
算出部２２とを備え、信号処理部２０で抽出された対を
なす各周波数帯域信号を導入してアクティブインテンシ
ティを算出するとともにリアクティブインテンシティを
算出し、漏水の特定位置を示すベクトルの大きさ（振動
振幅）にはリアクティブインテンシティの大きさ（絶対
値）を用い、上記ベクトルの方向にはアクティブインテ
ンシティの極性を用いて各周波数帯域毎に漏水位置特定
ベクトルを求めるベクトル演算部２３と、信号処理部２
０とリアクティブインテンシティ算出部２２との間に設
けられて信号処理部２０で抽出された対をなす各周波数
帯域信号のうちの基準となる信号については積分処理を
施してリアクティブインテンシティ算出部２２に導入す
る時間積分部２４と、ベクトル演算部２３によって演算
された各周波数帯域毎のインテンシティを合成演算する
演算部２５と、この演算部２５の出力に基づいて漏水位
置特定ベクトルを表示する表示部２６とで構成されてい
る。
れている水道管１１の一部分、たとえば地上に露出して
いる弁部分等に取付けられる漏水音検知センサ１２と、
漏水音信号に含まれる複数の周波数帯域を検出対象にし
て地表面１３に移動可能に配置され、対象とする周波数
帯域に対応させて間隔をそれぞれ異ならせた実質的に複
数対の地表振動検知センサ１４〜１８と、これら地表振
動検知センサ１４〜１８の出力信号から漏水音検知セン
サ１２で得られた信号とは関わりのない雑音信号を相互
相関処理でそれぞれ除去する相互相関処理部１９と、こ
の処理部１９を経た各対をなす信号を用いて地表面１３
に伝わる漏水エネルギのうちの上記信号を出力した対を
なす地表振動検知センサが対象にしている周波数成分を
それぞれ抽出処理する信号処理部２０と、アクティブイ
ンテンシティ算出部２１とリアクティブインテンシティ
算出部２２とを備え、信号処理部２０で抽出された対を
なす各周波数帯域信号を導入してアクティブインテンシ
ティを算出するとともにリアクティブインテンシティを
算出し、漏水の特定位置を示すベクトルの大きさ（振動
振幅）にはリアクティブインテンシティの大きさ（絶対
値）を用い、上記ベクトルの方向にはアクティブインテ
ンシティの極性を用いて各周波数帯域毎に漏水位置特定
ベクトルを求めるベクトル演算部２３と、信号処理部２
０とリアクティブインテンシティ算出部２２との間に設
けられて信号処理部２０で抽出された対をなす各周波数
帯域信号のうちの基準となる信号については積分処理を
施してリアクティブインテンシティ算出部２２に導入す
る時間積分部２４と、ベクトル演算部２３によって演算
された各周波数帯域毎のインテンシティを合成演算する
演算部２５と、この演算部２５の出力に基づいて漏水位
置特定ベクトルを表示する表示部２６とで構成されてい
る。
【００３９】ここで、地表振動センサ１４〜１８は、対
をなすもの同士が所定の間隔を設けて直線上に配置され
ている。すなわち、この例では図２に示すように、地表
振動検知センサ１４を共通（基準）センサとして用い、
この地表振動検知センサ１４から距離Ｘ1 ，Ｘ2 ，
Ｘ3 ，Ｘ4 離れた位置に地表振動検知センサ１５〜１８
が配置されている。この配置から判るように、地表振動
検知センサ１４と１５、地表振動検知センサ１４と１
６、地表振動検知センサ１４と１７、地表振動検知セン
サ１４と１８とがそれぞれ対をなしている。
をなすもの同士が所定の間隔を設けて直線上に配置され
ている。すなわち、この例では図２に示すように、地表
振動検知センサ１４を共通（基準）センサとして用い、
この地表振動検知センサ１４から距離Ｘ1 ，Ｘ2 ，
Ｘ3 ，Ｘ4 離れた位置に地表振動検知センサ１５〜１８
が配置されている。この配置から判るように、地表振動
検知センサ１４と１５、地表振動検知センサ１４と１
６、地表振動検知センサ１４と１７、地表振動検知セン
サ１４と１８とがそれぞれ対をなしている。
【００４０】さらに詳しく説明すると、地表振動検知セ
ンサ１４と１５とは、たとえば1kHz成分を対象にしてＸ
1 ＝20cmに設定されており、地表振動検知センサ１４と
１６とはたとえば600Hz 成分を対象にしてＸ2 ＝30cmに
設定されており、地表振動検知センサ１４と１７とはた
とえば400Hz 成分を対象にしてＸ3 ＝40cmに設定されて
おり、地表振動検知センサ１４と１８とはたとえば200H
z 成分を対象にしてＸ4 ＝50cmに設定されている。
ンサ１４と１５とは、たとえば1kHz成分を対象にしてＸ
1 ＝20cmに設定されており、地表振動検知センサ１４と
１６とはたとえば600Hz 成分を対象にしてＸ2 ＝30cmに
設定されており、地表振動検知センサ１４と１７とはた
とえば400Hz 成分を対象にしてＸ3 ＝40cmに設定されて
おり、地表振動検知センサ１４と１８とはたとえば200H
z 成分を対象にしてＸ4 ＝50cmに設定されている。
【００４１】これら地表振動検知センサ１４〜１８は、
支持部材２７によって共通に支持されており、水道管１
１の配設経路に沿って移動可能に地表面３上に載置され
る。なお、対をなして配置される地表振動検知センサの
設け方は、上記例に限られるものではなく、たとえば図
３に示すように、各周波数成分専用のものを２個ずつ用
意し、これら対をなすもの同士で間隔の小さいものほど
中央部に位置するように配置してもよい。また、地表振
動検知センサの移動には、動力車を用いることもできる
が、作業者の手で行うようにしてもよい。また、一対の
地表振動検知センサだけ用意し、両センサの間隔を作業
者の手で順次変更する方式も採用できる。
支持部材２７によって共通に支持されており、水道管１
１の配設経路に沿って移動可能に地表面３上に載置され
る。なお、対をなして配置される地表振動検知センサの
設け方は、上記例に限られるものではなく、たとえば図
３に示すように、各周波数成分専用のものを２個ずつ用
意し、これら対をなすもの同士で間隔の小さいものほど
中央部に位置するように配置してもよい。また、地表振
動検知センサの移動には、動力車を用いることもできる
が、作業者の手で行うようにしてもよい。また、一対の
地表振動検知センサだけ用意し、両センサの間隔を作業
者の手で順次変更する方式も採用できる。
【００４２】ここで、相互相関処理部１９、演算部２
０、ベクトル演算部２３、時間積分部２４、演算部２５
における、たとえば地表振動検知センサ１４と１５との
対、地表振動検知センサ１４と１６との対の処理系を代
表して取り出すと、これらの処理系は図４に示すように
構成されている。
０、ベクトル演算部２３、時間積分部２４、演算部２５
における、たとえば地表振動検知センサ１４と１５との
対、地表振動検知センサ１４と１６との対の処理系を代
表して取り出すと、これらの処理系は図４に示すように
構成されている。
【００４３】この構成から判るように、演算部２０を構
成しているデジタルフィルタによって抽出された対象と
する周波数帯域の対をなす信号を、一方においてはアク
ティブインテンシティ算出部２１に導入してアクティブ
インテンシティを算出させ、他方においては基準となる
地表振動検知センサ１４の信号のみを時間積分部２４、
具体的には積分器に通し、この積分器を通った信号と残
りの信号とをリアクティブインテンシティ算出部２２に
導入してリアクティブインテンシティを算出させ、さら
にベクトル演算部２３において漏水の特定位置を示すベ
クトルの大きさには前記リアクティブインテンシティの
大きさを用い、上記ベクトルの方向には前記アクティブ
インテンシティの極性を用いる漏水位置特定ベクトルを
演算させるようにしている。
成しているデジタルフィルタによって抽出された対象と
する周波数帯域の対をなす信号を、一方においてはアク
ティブインテンシティ算出部２１に導入してアクティブ
インテンシティを算出させ、他方においては基準となる
地表振動検知センサ１４の信号のみを時間積分部２４、
具体的には積分器に通し、この積分器を通った信号と残
りの信号とをリアクティブインテンシティ算出部２２に
導入してリアクティブインテンシティを算出させ、さら
にベクトル演算部２３において漏水の特定位置を示すベ
クトルの大きさには前記リアクティブインテンシティの
大きさを用い、上記ベクトルの方向には前記アクティブ
インテンシティの極性を用いる漏水位置特定ベクトルを
演算させるようにしている。
【００４４】このような構成であると、次のようにして
漏水位置特定ベクトルが算出される。なお、ここでは説
明を簡単にするために、対をなす地表振動検知センサ１
４，１５に注目して説明する。
漏水位置特定ベクトルが算出される。なお、ここでは説
明を簡単にするために、対をなす地表振動検知センサ１
４，１５に注目して説明する。
【００４５】先に説明したように、地表振動検知センサ
１４，１５によって検出された信号のうち相互相関処理
部１９によって相関処理された信号は演算部２０に導入
され、この演算部２０によって対象とする周波数帯域の
信号成分のみが抽出される。抽出された一対の対象とす
る周波数帯域信号は、一方においてはアクティブインテ
ンシティ算出部２１に導入される。このアクティブイン
テンシティ算出部２１で算出されるアクティブインテン
シティＩａは、（１）式によって表現され、測定位置に
よって図５（ｂ）に示すように変化する。
１４，１５によって検出された信号のうち相互相関処理
部１９によって相関処理された信号は演算部２０に導入
され、この演算部２０によって対象とする周波数帯域の
信号成分のみが抽出される。抽出された一対の対象とす
る周波数帯域信号は、一方においてはアクティブインテ
ンシティ算出部２１に導入される。このアクティブイン
テンシティ算出部２１で算出されるアクティブインテン
シティＩａは、（１）式によって表現され、測定位置に
よって図５（ｂ）に示すように変化する。
【００４６】抽出された一対の対象とする周波数帯域信
号は、他方においては時間積分部２４を介して（一方の
信号のみ）リアクティブインテンシティ算出部２２に導
入される。このりアクティブインテンシティ算出部２１
で算出されるリアクティブインテンシティＩｒは、理想
的な条件においては（３）式によって表現され、測定位
置によって図５（ｃ）に示すように変化する。
号は、他方においては時間積分部２４を介して（一方の
信号のみ）リアクティブインテンシティ算出部２２に導
入される。このりアクティブインテンシティ算出部２１
で算出されるリアクティブインテンシティＩｒは、理想
的な条件においては（３）式によって表現され、測定位
置によって図５（ｃ）に示すように変化する。
【００４７】ベクトル演算部２３は、漏水の特定位置を
示すベクトルの大きさにはリアクティブインテンシティ
Ｉｒの大きさを用い、ベクトルの方向にはアクティブイ
ンテンシティＩａの極性を用いてインテンシティＩを演
算する。このインテンシティＩは、（２）式によって表
現され、測定位置によって図５（ｄ）に示すように変化
する。
示すベクトルの大きさにはリアクティブインテンシティ
Ｉｒの大きさを用い、ベクトルの方向にはアクティブイ
ンテンシティＩａの極性を用いてインテンシティＩを演
算する。このインテンシティＩは、（２）式によって表
現され、測定位置によって図５（ｄ）に示すように変化
する。
【００４８】そして、最終的に演算部２５を介して表示
部２６に表示される漏水位置特定ベクトルは、図５
（ａ）に示すように、漏水位置に近づくほど大きくな
り、真上で最大となり、これを境にして極性が反転して
小さくなっていく。
部２６に表示される漏水位置特定ベクトルは、図５
（ａ）に示すように、漏水位置に近づくほど大きくな
り、真上で最大となり、これを境にして極性が反転して
小さくなっていく。
【００４９】したがって、この例に係る装置によれば、
漏水エネルギが最も大きい位置である漏水位置真上にお
ける漏水位置特定ベクトルの大きさを最大にすることが
でき、効率よく精度の高い漏水位置推定が可能となる。
なお、この例に係る装置では漏水位置から離れた位置で
もインテンシティが反転しないように振幅項に絶対値を
付けている。
漏水エネルギが最も大きい位置である漏水位置真上にお
ける漏水位置特定ベクトルの大きさを最大にすることが
でき、効率よく精度の高い漏水位置推定が可能となる。
なお、この例に係る装置では漏水位置から離れた位置で
もインテンシティが反転しないように振幅項に絶対値を
付けている。
【００５０】ところで、この例に係る装置ではリアクテ
ィブインテンシティＩｒの演算に当たって、基準となる
信号については時間積分部２４において積分された信号
を用いている。このため、算出されたリアクティブイン
テンシティＩｒ′は、（４）式に示すように、振幅に−
（１／ω）、つまり−（１／２πｆ）の重み付けがなさ
れたものとなる。
ィブインテンシティＩｒの演算に当たって、基準となる
信号については時間積分部２４において積分された信号
を用いている。このため、算出されたリアクティブイン
テンシティＩｒ′は、（４）式に示すように、振幅に−
（１／ω）、つまり−（１／２πｆ）の重み付けがなさ
れたものとなる。
【００５１】
【数４】
【００５２】しかし、このように、（２πｆ）で除算す
る重み付けが行われても、重み付けの影響はほとんど表
れない。
る重み付けが行われても、重み付けの影響はほとんど表
れない。
【００５３】すなわち、（２πｆ）で除算することは、
換言すると、高周波数の影響を小さくし、低周波数を強
調することになる。地表に伝わる漏水信号が低周波数域
から高周波数域まで同じレベルの場合には、低周波数の
信号成分で算出されたインテンシティの数値の方が高周
波数よりも支配的となり、真のインテンシティに近づけ
難くなる。しかし、実際には、地中にある音源から地表
に音が伝わる間に距離減衰が働き、高周波数域ほど減衰
が大きい。
換言すると、高周波数の影響を小さくし、低周波数を強
調することになる。地表に伝わる漏水信号が低周波数域
から高周波数域まで同じレベルの場合には、低周波数の
信号成分で算出されたインテンシティの数値の方が高周
波数よりも支配的となり、真のインテンシティに近づけ
難くなる。しかし、実際には、地中にある音源から地表
に音が伝わる間に距離減衰が働き、高周波数域ほど減衰
が大きい。
【００５４】図６（ａ），（ｂ）には、図６（ｃ）に示
す配置で実際に地中に加振源を埋設して漏水信号相当の
加速度レベルで加振し、その加振源信号と地表面に伝わ
る信号との間でのコヒーレンス（相関度）を調べた結果
が示されており、また図７には地表でのコヒーレンスの
分布を調べた結果が示されている。これらから判るよう
に、漏水位置近傍であっても、高周波数の信号ほど地表
に伝わり難い。したがって、地表で検出される漏水信号
はもともと高周波成分が小さいので、積分処理の影響を
受け難いことになる。
す配置で実際に地中に加振源を埋設して漏水信号相当の
加速度レベルで加振し、その加振源信号と地表面に伝わ
る信号との間でのコヒーレンス（相関度）を調べた結果
が示されており、また図７には地表でのコヒーレンスの
分布を調べた結果が示されている。これらから判るよう
に、漏水位置近傍であっても、高周波数の信号ほど地表
に伝わり難い。したがって、地表で検出される漏水信号
はもともと高周波成分が小さいので、積分処理の影響を
受け難いことになる。
【００５５】発明者らは、さらに「積分処理によるイン
テンシティに及ぼす影響」をシミュレーションによって
考察した。
テンシティに及ぼす影響」をシミュレーションによって
考察した。
【００５６】(1) センサ間隔一定（周波数分割なし）に
おけるシミュレーション結果 これは間隔が一定に保たれた一対の地表振動検知センサ
を用いる場合である。この条件で、真のリアクティブイ
ンテンシティＩｒを用いた場合、ＷＬＩ法で得られるイ
ンテンシティＩ（ｘ，ｙ）は（５）式のようになる。
おけるシミュレーション結果 これは間隔が一定に保たれた一対の地表振動検知センサ
を用いる場合である。この条件で、真のリアクティブイ
ンテンシティＩｒを用いた場合、ＷＬＩ法で得られるイ
ンテンシティＩ（ｘ，ｙ）は（５）式のようになる。
【００５７】
【数５】
【００５８】これに対して、時間積分処理を導入して算
出されたリアクティブインテンシティを用いた場合、Ｗ
ＬＩ法で得られるインテンシティＩ（ｘ，ｙ）は、
（６）式のようになる。
出されたリアクティブインテンシティを用いた場合、Ｗ
ＬＩ法で得られるインテンシティＩ（ｘ，ｙ）は、
（６）式のようになる。
【００５９】
【数６】
【００６０】図８、図９には配管の埋設深度等の条件を
同じくし、センサ間隔のみを０．２ｍ、０．４ｍと異な
らせて、（５）式に基づいて算出された漏水位置特定ベ
クトル（図８（ｂ）、図９（ｂ））と、（６）式に基づ
いて算出された漏水位置特定ベクトル（図８（ａ）、図
９（ａ））とが示されている。
同じくし、センサ間隔のみを０．２ｍ、０．４ｍと異な
らせて、（５）式に基づいて算出された漏水位置特定ベ
クトル（図８（ｂ）、図９（ｂ））と、（６）式に基づ
いて算出された漏水位置特定ベクトル（図８（ａ）、図
９（ａ））とが示されている。
【００６１】これらの図から判るように、リアクティブ
インテンシティの算出過程に積分処理を導入し、この結
果として（１／２πｆ）の重み付けが行われても、２０
０Ｈｚから１ｋＨｚのランダム音を対象にしたインテン
シティの結果にはほとんど変わりがない。これは、高周
波域では距離減衰が大きいため、低周波域に比べて地表
に伝わる漏水インテンシティそのものが小さく、（１／
２πｆ）の重み付けの影響が現れにくいことによる。
インテンシティの算出過程に積分処理を導入し、この結
果として（１／２πｆ）の重み付けが行われても、２０
０Ｈｚから１ｋＨｚのランダム音を対象にしたインテン
シティの結果にはほとんど変わりがない。これは、高周
波域では距離減衰が大きいため、低周波域に比べて地表
に伝わる漏水インテンシティそのものが小さく、（１／
２πｆ）の重み付けの影響が現れにくいことによる。
【００６２】したがって、リアクティブインテンシティ
の算出過程に積分処理を導入しても真のリアクティブイ
ンテンシティを用いた場合とほとんど同じ結果が得ら
れ、漏水位置特定精度を低下させることはないといえ
る。
の算出過程に積分処理を導入しても真のリアクティブイ
ンテンシティを用いた場合とほとんど同じ結果が得ら
れ、漏水位置特定精度を低下させることはないといえ
る。
【００６３】(2) 周波数帯域毎にセンサ間隔を変えた条
件におけるシミュレーション結果 ここでは、間隔の異なる２対の地表振動検知センサを用
いる場合について説明する。このように対象とする周波
数帯域に最適な間隔の地表振動検知センサを用いる理由
は、地表面と地表振動検知センサとの接触状態から生じ
る検出時間遅れ（位相遅れ）による漏水位置特定精度の
低下を抑制するためで、この点については後述する。
件におけるシミュレーション結果 ここでは、間隔の異なる２対の地表振動検知センサを用
いる場合について説明する。このように対象とする周波
数帯域に最適な間隔の地表振動検知センサを用いる理由
は、地表面と地表振動検知センサとの接触状態から生じ
る検出時間遅れ（位相遅れ）による漏水位置特定精度の
低下を抑制するためで、この点については後述する。
【００６４】対象とする周波数域を低周波数帯域と高周
波数帯域とに分け、各周波数帯域に適した間隔を持つ実
質的に２対の地表振動検知センサ用いる条件で、真のリ
アクティブインテンシティＩｒを用いた場合、ＷＬＩ法
で算出されるインテンシティＩ（ｘ，ｙ）は（７）式の
ようになる。
波数帯域とに分け、各周波数帯域に適した間隔を持つ実
質的に２対の地表振動検知センサ用いる条件で、真のリ
アクティブインテンシティＩｒを用いた場合、ＷＬＩ法
で算出されるインテンシティＩ（ｘ，ｙ）は（７）式の
ようになる。
【００６５】
【数７】
【００６６】これに対して、時間積分処理を導入してリ
アクティブインテンシティを算出した場合、ＷＬＩ法で
得られるインテンシティＩ（ｘ，ｙ）は（８）式のよう
になる。
アクティブインテンシティを算出した場合、ＷＬＩ法で
得られるインテンシティＩ（ｘ，ｙ）は（８）式のよう
になる。
【００６７】
【数８】
【００６８】図１０、図１１には配管の埋設深度、セン
サ間隔（低域０．４ｍ、高域０．２ｍ）等の条件を同じ
くし、センサ設置誤差（位相差）のみを０度、３０度と
異ならせて、（７）式に基づいて算出された漏水位置特
定ベクトル（図１０（ｂ）、図１１（ｂ））と、（８）
式に基づいて算出された漏水位置特定ベクトル（図１０
（ａ）、図１１（ａ））とが示されている。
サ間隔（低域０．４ｍ、高域０．２ｍ）等の条件を同じ
くし、センサ設置誤差（位相差）のみを０度、３０度と
異ならせて、（７）式に基づいて算出された漏水位置特
定ベクトル（図１０（ｂ）、図１１（ｂ））と、（８）
式に基づいて算出された漏水位置特定ベクトル（図１０
（ａ）、図１１（ａ））とが示されている。
【００６９】これらの図から判るように、センサ設置誤
差が０度、３０度の場合でも、またセンサ間隔を対象と
する周波数帯域毎に変化させた場合においても、リアク
ティブインテンシティの算出過程に積分処理を導入した
ことによる漏水位置特定精度劣化は生じていない。した
がって、センサ設置誤差対策として、複数対の地表振動
検知センサを用いる場合であってもリアクティブインテ
ンシティの算出過程に積分処理を導入することによって
真のリアクティブインテンシティを用いた場合とほとん
ど同程度の漏水位置特定精度を確保できるといえる。
差が０度、３０度の場合でも、またセンサ間隔を対象と
する周波数帯域毎に変化させた場合においても、リアク
ティブインテンシティの算出過程に積分処理を導入した
ことによる漏水位置特定精度劣化は生じていない。した
がって、センサ設置誤差対策として、複数対の地表振動
検知センサを用いる場合であってもリアクティブインテ
ンシティの算出過程に積分処理を導入することによって
真のリアクティブインテンシティを用いた場合とほとん
ど同程度の漏水位置特定精度を確保できるといえる。
【００７０】このように、リアクティブインテンシティ
の算出過程に積分処理を導入しても真のリアクティブイ
ンテンシティに近づけることができるので、信号処理系
の複雑化をなく漏水位置特定精度を向上させることがで
きる。
の算出過程に積分処理を導入しても真のリアクティブイ
ンテンシティに近づけることができるので、信号処理系
の複雑化をなく漏水位置特定精度を向上させることがで
きる。
【００７１】次に、リアクティブインテンシティの算出
過程に積分処理を導入したことの有効性を検証した模擬
漏水実験および実際の漏水現場試験の結果を説明する。
過程に積分処理を導入したことの有効性を検証した模擬
漏水実験および実際の漏水現場試験の結果を説明する。
【００７２】（１）模擬漏水実験
鉛管に漏水個所となる直径２mmの穴を開け、この鉛管を
地中１ｍに埋設した。メータますに漏水音検知センサを
取り付けた。図１２に示すように、地表面に低周波数帯
域用および高周波数帯域用の３個１組の地表振動検知セ
ンサを配置し、鉛管から漏水させている状態で鉛管に沿
って４０cm間隔で図１に示した装置と同様の手法で漏水
位置特定ベクトルを測定した。その結果、図１３に示す
ように、漏水位置の真上で漏水位置特定ベクトルが最大
となり、リアクティブインテンシティの算出過程に積分
処理を導入したことの有効性が確認された。
地中１ｍに埋設した。メータますに漏水音検知センサを
取り付けた。図１２に示すように、地表面に低周波数帯
域用および高周波数帯域用の３個１組の地表振動検知セ
ンサを配置し、鉛管から漏水させている状態で鉛管に沿
って４０cm間隔で図１に示した装置と同様の手法で漏水
位置特定ベクトルを測定した。その結果、図１３に示す
ように、漏水位置の真上で漏水位置特定ベクトルが最大
となり、リアクティブインテンシティの算出過程に積分
処理を導入したことの有効性が確認された。
【００７３】（２）漏水現場試験
地中１ｍに埋設された漏水配管を用いて実際の漏水位置
特定を実施した。図１４に示すように、メータますに漏
水音検知センサを取り付けるとともに地表面に低周波数
帯域用および高周波数帯域用の３個１組の地表振動検知
センサを配置し、配管に沿って４０cm間隔で図１に示し
た装置と同様の手法で漏水位置特定ベクトルを測定し
た。その結果、図１４に示すように実際の漏水位置周囲
５０cm内で漏水位置特定ベクトルの向きが反転し、リア
クティブインテンシティの算出過程に積分処理を導入し
たことの有効性が確認された。
特定を実施した。図１４に示すように、メータますに漏
水音検知センサを取り付けるとともに地表面に低周波数
帯域用および高周波数帯域用の３個１組の地表振動検知
センサを配置し、配管に沿って４０cm間隔で図１に示し
た装置と同様の手法で漏水位置特定ベクトルを測定し
た。その結果、図１４に示すように実際の漏水位置周囲
５０cm内で漏水位置特定ベクトルの向きが反転し、リア
クティブインテンシティの算出過程に積分処理を導入し
たことの有効性が確認された。
【００７４】なお、上述した例では、リアクティブイン
テンシティの算出過程に時間積分処理を導入し、振動振
幅に周波数重み付けのあるリアクティブインテンシティ
を算出しているが、これに限られるものではなく、相互
相関処理を経た一方の信号と時間積分処理または時間微
分処理された他方の信号とを用いて最終的に重み付けの
逆数を掛けることによって、周波数重み付けの除かれた
真のリアクティブインテンシテイを演算するようにして
もよい。
テンシティの算出過程に時間積分処理を導入し、振動振
幅に周波数重み付けのあるリアクティブインテンシティ
を算出しているが、これに限られるものではなく、相互
相関処理を経た一方の信号と時間積分処理または時間微
分処理された他方の信号とを用いて最終的に重み付けの
逆数を掛けることによって、周波数重み付けの除かれた
真のリアクティブインテンシテイを演算するようにして
もよい。
【００７５】次に、対象とする周波数帯域に最適な間隔
の地表振動検知センサを用いなけれならない理由につい
て説明する。
の地表振動検知センサを用いなけれならない理由につい
て説明する。
【００７６】図１５から図１７には地表振動検知センサ
の間隔を２０cm、５０cmに設定し、これらの間隔で２０
０Ｈｚ、４００Ｈｚ、６００Ｈｚ、８００Ｈｚ、１ｋＨ
ｚの漏水信号をＷＬＩ法で測定したときのインテンシテ
ィＩの変化が示されている。インテンシティＩの極性が
＋（−）から−（＋）に急激に変わる場所が漏水箇所の
真上に当たる。なお、これらの図で破線はＡＩ法で測定
されたインテンシティを示している。
の間隔を２０cm、５０cmに設定し、これらの間隔で２０
０Ｈｚ、４００Ｈｚ、６００Ｈｚ、８００Ｈｚ、１ｋＨ
ｚの漏水信号をＷＬＩ法で測定したときのインテンシテ
ィＩの変化が示されている。インテンシティＩの極性が
＋（−）から−（＋）に急激に変わる場所が漏水箇所の
真上に当たる。なお、これらの図で破線はＡＩ法で測定
されたインテンシティを示している。
【００７７】これらの図において、たとえばセンサ間隔
５０cmを例にとると、周波数が２００Ｈｚのとには、漏
水箇所の真上だけにおいて極性転換が行われるが、周波
数が高くなると、漏水箇所の真上以外の場所においても
極性転換が行われる。前述の如く、ＷＬＩ法では極性転
換を有効に活用しているので、漏水箇所の真上以外の複
数の場所において極性転換が行われると漏水位置特定精
度が低下する。したがって、高い漏水位置特定精度を確
保するには、対象とする周波数が２００Ｈｚの場合には
センサ間隔をたとえば５０cmにする必要がある。同様
に、対象とする周波数が１ｋＨｚの場合には、図１７に
示されるように、センサ間隔が５０cmの場合より、セン
サ間隔が２０cmのほうが有利である。
５０cmを例にとると、周波数が２００Ｈｚのとには、漏
水箇所の真上だけにおいて極性転換が行われるが、周波
数が高くなると、漏水箇所の真上以外の場所においても
極性転換が行われる。前述の如く、ＷＬＩ法では極性転
換を有効に活用しているので、漏水箇所の真上以外の複
数の場所において極性転換が行われると漏水位置特定精
度が低下する。したがって、高い漏水位置特定精度を確
保するには、対象とする周波数が２００Ｈｚの場合には
センサ間隔をたとえば５０cmにする必要がある。同様
に、対象とする周波数が１ｋＨｚの場合には、図１７に
示されるように、センサ間隔が５０cmの場合より、セン
サ間隔が２０cmのほうが有利である。
【００７８】このように、同じセンサ間隔でも周波数に
よってインテンシティの変化特性に差が出るのは、地表
面と地表振動検知センサとの接触状態から生じるセンサ
設置誤差α、つまり検出時間遅れが周波数によって異な
ることによる。
よってインテンシティの変化特性に差が出るのは、地表
面と地表振動検知センサとの接触状態から生じるセンサ
設置誤差α、つまり検出時間遅れが周波数によって異な
ることによる。
【００７９】図１５から図１７に示す特性から判るよう
に、低周波数帯域では地表振動検知センサの間隔ｄを大
きくした方がセンサ設置誤差αの許容範囲が拡大する。
一方、高周波数帯域では波長が短いことから、間隔ｄを
小さくした方がよいことが判る。
に、低周波数帯域では地表振動検知センサの間隔ｄを大
きくした方がセンサ設置誤差αの許容範囲が拡大する。
一方、高周波数帯域では波長が短いことから、間隔ｄを
小さくした方がよいことが判る。
【００８０】
【発明の効果】以上説明したように、本発明によれば、
人間の聴覚に頼ることなく、また雑音の多い昼間におい
ても容易に、しかも装置構成の複雑化を招くことなく地
中内配管に存在する漏水位置を精度よく特定することが
できる。
人間の聴覚に頼ることなく、また雑音の多い昼間におい
ても容易に、しかも装置構成の複雑化を招くことなく地
中内配管に存在する漏水位置を精度よく特定することが
できる。
【図１】本発明の一実施形態に係る漏水位置検出装置の
ブロック構成図
ブロック構成図
【図２】同装置における地表振動検知センサの配置例を
説明するための図
説明するための図
【図３】地表振動検知センサの異なる配置例を説明する
ための図
ための図
【図４】同装置の信号処理系を局部的に取り出して示す
図
図
【図５】同装置によって算出された漏水位置特定ベクト
ルの特徴を説明するための図
ルの特徴を説明するための図
【図６】地中に埋設された加振源の信号と地表面に伝わ
る信号との間でのコヒーレンスを調べた実験結果を示す
図
る信号との間でのコヒーレンスを調べた実験結果を示す
図
【図７】同実験において地表のコヒーレンス分布を調べ
た結果を示す図
た結果を示す図
【図８】リアクティブインテンシティの演算に積分処理
を導入したことによる影響を調べた結果を説明するため
の図
を導入したことによる影響を調べた結果を説明するため
の図
【図９】リアクティブインテンシティの演算に積分処理
を導入したことによる影響を調べた結果を説明するため
の図
を導入したことによる影響を調べた結果を説明するため
の図
【図１０】リアクティブインテンシティの演算に積分処
理を導入したことによる影響を調べた結果を説明するた
めの図
理を導入したことによる影響を調べた結果を説明するた
めの図
【図１１】リアクティブインテンシティの演算に積分処
理を導入したことによる影響を調べた結果を説明するた
めの図
理を導入したことによる影響を調べた結果を説明するた
めの図
【図１２】模擬漏水実験で採用したセンサ位置関係を説
明するための図
明するための図
【図１３】模擬漏水実験の結果を説明するための図
【図１４】漏水現場試験の状況と観測された漏水位置特
定ベクトルとの関係を示す図
定ベクトルとの関係を示す図
【図１５】ＷＬＩ法におけるセンサ間隔と周波数とイン
テンシティとの関係の一例を示す図
テンシティとの関係の一例を示す図
【図１６】ＷＬＩ法におけるセンサ間隔と周波数とイン
テンシティとの関係の一例を示す図
テンシティとの関係の一例を示す図
【図１７】ＷＬＩ法におけるセンサ間隔と周波数とイン
テンシティとの関係の一例を示す図
テンシティとの関係の一例を示す図
【図１８】ＡＩ法を用いた従来の漏水位置検出装置の概
略構成図
略構成図
【図１９】従来の装置で算出した漏水位置特定ベクトル
を説明するための図
を説明するための図
１１…水道管（配管）
１２…漏水音検知センサ
１３…地表面
１４〜１８，１４ａ，１４ｂ，１５ａ，１５ｂ，１６
ａ，１６ｂ，１７ａ，１７ｂ…地表振動検知センサ １９…相互相関処理部 ２０…信号処理部 ２１…アクティブインテンシティ算出部 ２２…リアクティブインテンシティ算出部 ２３…ベクトル演算部 ２４…時間積分部 ２５…演算部 ２６…漏水位置特定ベクトル表示装置 Ｅ…漏水位置
ａ，１６ｂ，１７ａ，１７ｂ…地表振動検知センサ １９…相互相関処理部 ２０…信号処理部 ２１…アクティブインテンシティ算出部 ２２…リアクティブインテンシティ算出部 ２３…ベクトル演算部 ２４…時間積分部 ２５…演算部 ２６…漏水位置特定ベクトル表示装置 Ｅ…漏水位置
フロントページの続き
(72)発明者 田崎 拓也
東京都新宿区西新宿２丁目８番１号 東
京都水道局内
(72)発明者 江波戸 明彦
神奈川県川崎市幸区小向東芝町１番地
株式会社東芝研究開発センター内
(72)発明者 佐藤 義之
東京都府中市東芝町１番地 株式会社東
芝府中工場内
(72)発明者 環 省二郎
東京都港区芝浦１丁目１番１号 株式会
社東芝本社事務所内
(56)参考文献 特開 平５−256726（ＪＰ，Ａ)
特開 平６−273223（ＪＰ，Ａ)
(58)調査した分野(Int.Cl.7，ＤＢ名)
G01M 3/24
G01F 1/00
Claims (8)
- 【請求項１】地中に埋設されている配管の漏水位置を推
定検出するに当り、上記配管の一部分に漏水音検知セン
サを取付けるとともに少なくとも一対の地表振動検知セ
ンサを地表面に移動可能に配置し、これら地表振動検知
センサの出力信号から前記漏水音検知センサで得られた
信号とは関わりのない雑音信号を相互相関処理でそれぞ
れ除去し、この処理を経た対をなす信号を用いて地表面
に伝わる漏水エネルギのアクティブインテンシティを演
算し、さらに前記相互相関処理を経た一方の信号と時間
積分処理された他方の信号とを用いて周波数重み付けの
あるリアクティブインテンシテイを演算し、漏水の特定
位置を示すベクトルの大きさには前記リアクティブイン
テンシティの大きさを用い、上記ベクトルの方向には前
記アクティブインテンシティの極性を用いて漏水位置特
定ベクトルを求めるようにしたことを特徴とする漏水位
置検出方法。 - 【請求項２】地中に埋設されている配管の漏水位置を推
定検出するに当り、上記配管の一部分に漏水音検知セン
サを取付けるとともに少なくとも一対の地表振動検知セ
ンサを地表面に移動可能に配置し、これら地表振動検知
センサの出力信号から前記漏水音検知センサで得られた
信号とは関わりのない雑音信号を相互相関処理でそれぞ
れ除去し、この処理を経た対をなす信号を用いて地表面
に伝わる漏水エネルギのアクティブインテンシティを演
算し、さらに前記相互相関処理を経た一方の信号と時間
積分処理または時間微分処理された他方の信号とを用い
て最終的に周波数重み付けの除かれた真のリアクティブ
インテンシテイを演算し、漏水の特定位置を示すベクト
ルの大きさには前記真のリアクティブインテンシティの
大きさを用い、上記ベクトルの方向には前記アクティブ
インテンシティの極性を用いて漏水位置特定ベクトルを
求めるようにしたことを特徴とする漏水位置検出方法。 - 【請求項３】漏水音に含まれている周波数成分に対応さ
せて間隔を異ならせた実質的に複数対の地表振動検知セ
ンサを移動可能に配置し、周波数成分毎に前記漏水位置
特定ベクトルをそれぞれ求めるようにしたことを特徴と
する請求項１または２に記載の漏水位置検出方法。 - 【請求項４】前記複数の周波数成分に対応させて配置さ
れる実質的に複数対の地表振動センサは、対象とする周
波数が低いほど間隔が広く設定されることを特徴とする
請求項３に記載の漏水位置検出方法。 - 【請求項５】地中に埋設されている配管の漏水位置を推
定検出するための装置であって、 前記配管の一部分に取付けられる漏水音検知センサと、 地表面に移動可能に配置される少なくとも一対の地表振
動検知センサと、 これら地表振動検知センサの出力信号から前記漏水音検
知センサで得られた信号とは関わりのない雑音信号を相
互相関処理でそれぞれ除去する手段と、 この手段による処理を経た対をなす信号を用いて地表面
に伝わる漏水エネルギのアクティブインテンシティを演
算する手段と、 前記相互相関処理を経た一方の信号と時間積分処理され
た他方の信号とを用いて周波数重み付けのあるリアクテ
ィブインテンシティを演算する手段と、 漏水の特定位置を示すベクトルの大きさには前記リアク
ティブインテンシティの大きさを用い、上記ベクトルの
方向には前記アクティブインテンシティの極性を用いて
漏水位置特定ベクトルを求める手段とを具備してなるこ
とを特徴とする漏水位置検出装置。 - 【請求項６】地中に埋設されている配管の漏水位置を推
定検出するための装置であって、 前記配管の一部分に取付けられる漏水音検知センサと、 地表面に移動可能に配置される少なくとも一対の地表振
動検知センサと、 これら地表振動検知センサの出力信号から前記漏水音検
知センサで得られた信号とは関わりのない雑音信号を相
互相関処理でそれぞれ除去する手段と、 この手段による処理を経た対をなす信号を用いて地表面
に伝わる漏水エネルギのアクティブインテンシティを演
算する手段と、 前記相互相関処理を経た一方の信号と時間積分処理また
は時間微分処理された他方の信号とを用いて最終的に周
波数重み付けの除かれた真のリアクティブインテンシテ
ィを演算する手段と、 漏水の特定位置を示すベクトルの大きさには前記真のリ
アクティブインテンシティの大きさを用い、上記ベクト
ルの方向には前記アクティブインテンシティの極性を用
いて漏水位置特定ベクトルを求める手段とを具備してな
ることを特徴とする漏水位置検出装置。 - 【請求項７】地表面に移動可能に配置される地表振動検
知センサは漏水音に含まれている周波数成分に対応させ
て間隔が異なるように配置された実質的に複数対設けら
れており、周波数成分毎に前記漏水位置特定ベクトルを
それぞれ求めるようにしていることを特徴とする請求項
５または６に記載の漏水位置検出装置。 - 【請求項８】前記複数の周波数成分に対応させて配置さ
れる実質的に複数対の地表振動センサは、対象とする周
波数が低いほど間隔が広く設定されることを特徴とする
請求項７に記載の漏水位置検出装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP05200098A JP3488623B2 (ja) | 1998-03-04 | 1998-03-04 | 漏水位置検出方法および漏水位置検出装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP05200098A JP3488623B2 (ja) | 1998-03-04 | 1998-03-04 | 漏水位置検出方法および漏水位置検出装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH11248591A JPH11248591A (ja) | 1999-09-17 |
| JP3488623B2 true JP3488623B2 (ja) | 2004-01-19 |
Family
ID=12902580
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP05200098A Expired - Fee Related JP3488623B2 (ja) | 1998-03-04 | 1998-03-04 | 漏水位置検出方法および漏水位置検出装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3488623B2 (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2014157539A1 (ja) * | 2013-03-29 | 2014-10-02 | 日本電気株式会社 | 欠陥分析装置、欠陥分析方法及びプログラム |
| CN111811635A (zh) * | 2020-07-07 | 2020-10-23 | 美钻石油钻采系统(上海)有限公司 | 水下目标监测装置 |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2001031308A1 (de) * | 1999-10-26 | 2001-05-03 | Peter Martinek | Verfahren und messsonde zur durchführung von messungen in wasserversorgungssystemen |
| JP5846015B2 (ja) | 2012-03-30 | 2016-01-20 | 日本電気株式会社 | 漏洩検知方法、漏水検知方法、漏洩検知装置および漏水検知装置 |
| CN108488637A (zh) * | 2018-03-22 | 2018-09-04 | 哈尔滨工业大学 | 一种基于自相关和近似熵原理的供水管道泄漏辨识方法 |
| KR102038689B1 (ko) | 2018-06-14 | 2019-10-30 | 한국원자력연구원 | 거리차-주파수 분석을 이용한 배관의 누설 감지장치 및 방법 |
-
1998
- 1998-03-04 JP JP05200098A patent/JP3488623B2/ja not_active Expired - Fee Related
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2014157539A1 (ja) * | 2013-03-29 | 2014-10-02 | 日本電気株式会社 | 欠陥分析装置、欠陥分析方法及びプログラム |
| CN111811635A (zh) * | 2020-07-07 | 2020-10-23 | 美钻石油钻采系统(上海)有限公司 | 水下目标监测装置 |
| CN111811635B (zh) * | 2020-07-07 | 2021-06-15 | 美钻石油钻采系统(上海)有限公司 | 水下目标监测装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH11248591A (ja) | 1999-09-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Pal et al. | Detecting & locating leaks in water distribution polyethylene pipes | |
| EP0552044B1 (en) | Method and apparatus for detecting the position of an abnormal site of a buried pipe | |
| US6453247B1 (en) | PC multimedia-based leak detection system for water transmission and distribution pipes | |
| JPH11201859A (ja) | 周波数帯域分割による配管漏洩検知方法 | |
| GB2421311A (en) | Assessing the size of a leak in a pipeline by detecting leak noise and pressure | |
| US5127267A (en) | Acoustic method for locating concealed pipe | |
| GB2318640A (en) | Underground leak location | |
| US10458878B2 (en) | Position determination device, leak detection system, position determination method, and computer-readable recording medium | |
| JP3488623B2 (ja) | 漏水位置検出方法および漏水位置検出装置 | |
| JP2013044612A (ja) | 埋設配管破損位置の検出方法及び装置 | |
| Yoon et al. | Improvement of cross-correlation technique for leak detection of a buried pipe in a tonal noisy environment | |
| JP3032185B2 (ja) | 配管の漏洩検知システム | |
| JP3488579B2 (ja) | 漏水位置検出方法および漏水位置検出装置 | |
| JPH10281921A (ja) | ガス配管の漏洩監視方法および装置 | |
| JP2004125628A (ja) | 配管の漏洩位置検知方法および装置 | |
| JP4306409B2 (ja) | 配管の漏洩位置検知方法および装置 | |
| JPH0894481A (ja) | 埋設管のガス漏洩検知方法 | |
| JPH11210999A (ja) | 相関法による管路系の漏洩位置特定方法 | |
| JP6789042B2 (ja) | 漏水場所の特定方法 | |
| Glentis et al. | Using Generalized Cross-Correlation estimators for leak signal velocity estimation and spectral region of operation selection | |
| WO2014157539A1 (ja) | 欠陥分析装置、欠陥分析方法及びプログラム | |
| RU2681552C1 (ru) | Способ обнаружения несанкционированных врезок в трубопровод | |
| Brennan et al. | Some recent research results on the use of acoustic methods to detect water leaks in buried plastic water pipes | |
| JP3262469B2 (ja) | 漏水検出装置 | |
| KR20110060632A (ko) | 산업설비의 결함감지를 위한 음향진동 복합 센싱유닛 및 이를 이용한 고압배관의 결함진단시스템 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081031 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081031 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091031 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101031 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111031 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111031 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121031 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131031 Year of fee payment: 10 |
|
| LAPS | Cancellation because of no payment of annual fees |