(57)【要約】【課題】本発明は、水道管における漏水位置を検出する処理において、センサと演算装置との間で伝送される情報量を低減することを目的とする。【解決手段】漏水検知システムは、水道管１４における２つの検知点のそれぞれに配置されたセンサ装置（１２−１，１２−２）から振動波形データを取得し、各振動波形データに対する演算処理に基づいて漏水位置を特定するものである。センサ装置１２−１および１２−２は、予め定められた距離Ｓを隔てて水道管１４の検知点に配置されている。各センサ装置が備えるセンサは音響振動を検出し、各センサ装置はセンサの検出値に基づいて振動波形データを生成する。各センサ装置は、コンセントレータ１０に振動波形データを送信する。コンセントレータ１０は、各振動波形データに基づいて測定データを生成し、演算装置２０に送信する。演算装置２０は、測定データに基づいて水道管１４の漏水位置を求める。【選択図】図１
2019124613.tif 000002
【特許請求の範囲】
【請求項１】
水道管における２つの検知点のそれぞれに配置されたセンサから、振動波形情報を取得する波形取得部と、
各前記センサから取得された振動波形情報についての相関処理であって、一方の波形をシフト時間だけ時間軸上でシフトした場合における当該シフト時間に対する相互相関値を求める相関処理を実行する相関処理部と、
前記シフト時間と前記相互相関値との関係に基づく測定情報であって、漏水位置を求める元となる測定情報を求める測定部と、を備え、
各前記振動波形情報は、
時間軸上に配列された複数のサンプル値によって時間波形を表す情報であり、
前記相関処理は、
所定サンプル数おきのサンプル値に対して相関演算を施す処理である、
ことを特徴とする漏水検出装置。
【請求項２】
水道管における２つの検知点のそれぞれに配置されたセンサから、振動波形情報を取得する波形取得部と、
各前記センサから取得された振動波形情報についての相関処理であって、一方の波形をシフト時間だけ時間軸上でシフトした場合における当該シフト時間に対する相互相関値を求める相関処理を実行する相関処理部と、
前記シフト時間と前記相互相関値との関係に基づく測定情報を送信する送信部と、を備え、
前記送信部は、
前記測定情報に基づいて漏水位置を求める漏水位置検出装置に前記測定情報を送信する、
ことを特徴とする漏水検出装置。
【請求項３】
請求項１または請求項２に記載の漏水検出装置において、
前記測定情報は、
前記シフト時間と前記相互相関値との関係における複数のピークについての、前記シフト時間、前記相互相関値のピーク値、および際立ち値であって、前記複数のピークを母集団としたときのピーク値の際立ちを示す際立ち値を含む、
ことを特徴とする漏水検出装置。  
【請求項４】
水道管における２つの検知点のそれぞれに配置されたセンサから、振動波形情報を取得する波形取得部と、
各前記センサから取得された振動波形情報についての相関処理であって、一方の波形をシフト時間だけ時間軸上でシフトした場合における当該シフト時間に対する相互相関値を求める相関処理を実行する相関処理部と、
前記シフト時間と前記相互相関値との関係に基づく測定情報であって、漏水位置を求める元となる測定情報を求める測定部と、を備え、
前記測定情報は、
前記シフト時間と前記相互相関値との関係における複数のピークについての、前記シフト時間、前記相互相関値のピーク値、および際立ち値であって、前記複数のピークを母集団としたときのピーク値の際立ちを示す際立ち値を含む、
ことを特徴とする漏水検出装置。
【請求項５】
請求項３または請求項４に記載の漏水検出装置において、
前記際立ち値は、
複数の前記シフト時間に対して求められた複数の前記相互相関値を母集団とする偏差値であることを特徴とする漏水検出装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、漏水検出装置に関し、特に、複数のセンサを用いて漏水位置を検出する装置に関する。
【背景技術】
【０００２】
地中に埋められた水道管からの漏水が問題となっている。漏水を検出する方法としては、水道管に設けられた水栓において作業員が振動音を耳で確認するものがあった。この場合、作業員は、複数箇所の水栓で振動音を確認し、経験や勘に基づいて漏水が生じている位置を特定する必要があるため、作業員の労力が大きいものとなっていた。
【０００３】
このような問題点を解決する手段として、以下の特許文献１および２に記載されているものがある。これらの特許文献には、配管に２つのセンサを設け、各センサによって検出された振動データの相互相関値（一方の振動波形と、他方の振動波形を時間軸上で所定時間だけシフトしたものとが似ている度合いを表す値）に基づいて、配管からの流体の漏れを検出する技術が記載されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開平０８−２２６８６５号公報
【特許文献２】特表２００３−５０２６７８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
特許文献１および２に記載されているように、２つのセンサによって検出された振動データの相互相関値を求めることで漏れが生じている位置を特定する場合、相互相関処理の対象となる情報量が多い。そのため、演算装置が各センサから振動データを取得して相互相関処理を実行する場合、各センサから演算装置に送信される情報量が多くなり、迅速な演算の妨げになってしまう。
【０００６】
本発明は、水道管における漏水位置を検出する処理において、センサと演算装置との間で伝送される情報量を低減することを目的とする。
【課題を解決するための手段】
【０００７】
本発明は、水道管における２つの検知点のそれぞれに配置されたセンサから、振動波形情報を取得する波形取得部と、各前記センサから取得された振動波形情報についての相関処理であって、一方の波形をシフト時間だけ時間軸上でシフトした場合における当該シフト時間に対する相互相関値を求める相関処理を実行する相関処理部と、前記シフト時間と前記相互相関値との関係に基づく測定情報であって、漏水位置を求める元となる測定情報を求める測定部と、を備え、各前記振動波形情報は、時間軸上に配列された複数のサンプル値によって時間波形を表す情報であり、前記相関処理は、所定サンプル数おきのサンプル値に対して相関演算を施す処理である、ことを特徴とする。
【０００８】
また、本発明は、水道管における２つの検知点のそれぞれに配置されたセンサから、振動波形情報を取得する波形取得部と、各前記センサから取得された振動波形情報についての相関処理であって、一方の波形をシフト時間だけ時間軸上でシフトした場合における当該シフト時間に対する相互相関値を求める相関処理を実行する相関処理部と、前記シフト時間と前記相互相関値との関係に基づく測定情報を送信する送信部と、を備え、前記送信部は、前記測定情報に基づいて漏水位置を求める漏水位置検出装置に前記測定情報を送信する、ことを特徴とする。
【０００９】
望ましくは、前記測定情報は、前記シフト時間と前記相互相関値との関係における複数のピークについての、前記シフト時間、前記相互相関値のピーク値、および際立ち値であって、前記複数のピークを母集団としたときのピーク値の際立ちを示す際立ち値を含む。
【００１０】
また、本発明は、水道管における２つの検知点のそれぞれに配置されたセンサから、振動波形情報を取得する波形取得部と、各前記センサから取得された振動波形情報についての相関処理であって、一方の波形をシフト時間だけ時間軸上でシフトした場合における当該シフト時間に対する相互相関値を求める相関処理を実行する相関処理部と、前記シフト時間と前記相互相関値との関係に基づく測定情報であって、漏水位置を求める元となる測定情報を求める測定部と、を備え、前記測定情報は、前記シフト時間と前記相互相関値との関係における複数のピークについての、前記シフト時間、前記相互相関値のピーク値、および際立ち値であって、前記複数のピークを母集団としたときのピーク値の際立ちを示す際立ち値を含む、ことを特徴とする。
【００１１】
望ましくは、前記際立ち値は、複数の前記シフト時間に対して求められた複数の前記相互相関値を母集団とする偏差値である。
【発明の効果】
【００１２】
本発明によれば、水道管における漏水位置を検出する処理において、センサと演算装置との間で伝送される情報量を低減することができる。
【図面の簡単な説明】
【００１３】
【図１】本発明の実施形態に係る漏水検出システムを示す図である。
【図２】コンセントレータおよび各センサ装置が実行する処理を示す図である。
【図３】各センサ装置による振動波形データが示す振動波形を示す図である。
【図４】相互相関値を求める演算を概念的に示す図である。
【図５】シフト時間に対する相互相関値の特性の例を示す図である。
【図６】相互相関値を求める演算を概念的に示す図である。
【図７】コンセントレータのハードウエアの例を示す図である。
【発明を実施するための形態】
【００１４】
図１には本発明の実施形態に係る漏水検出システムが示されている。漏水検出システムは、水道管１４における２つの検知点のそれぞれに配置されたセンサ装置（１２−１，１２−２）から振動波形データを取得し、各振動波形データに対する演算処理に基づいて漏水位置を特定するものである。図１では、センサ装置１２−１から漏水箇所１６までの距離ｄｗによって漏水位置が表されている。
【００１５】
漏水検出システムは、コンセントレータ１０、センサ装置１２−１およびセンサ装置１２−２を備える（以下の説明では、センサ装置１２−１および１２−２の両者に対して符号「１２」を用いる。）。センサ装置１２−１および１２−２は、予め定められた距離Ｓを隔てて水道管１４の検知点に配置されている。各検知点は、例えば、地中の水道管１４から地上に引き出された給水栓や空気抜き用バルブにある。各センサ装置１２は、地中において水道管１４に接触して固定されてもよい。
【００１６】
各センサ装置１２はセンサを備えている。センサには、例えば、圧電素子やコンデンサマイク等、音や機械振動等の音響振動を検出する素子が用いられる。各センサ装置１２が備えるセンサは音響振動を検出し、各センサ装置１２はセンサの検出値に基づいて振動波形データを生成する。振動波形データは、時間軸に対して振動の振れ幅を対応付けたものである。
【００１７】
漏水検出装置としてのコンセントレータ１０は、各センサ装置１２との間で無線通信を行う。各センサ装置１２は、コンセントレータ１０に振動波形データを送信する。コンセントレータ１０は、各センサ装置１２から送信された振動波形データに基づいて、漏水位置を求める元となる測定データを生成し、インターネット等の通信網１８を介して演算装置２０に送信する。演算装置２０は、漏水位置検出装置としての機能を有し、測定データに基づいて水道管１４の漏水位置を求める。演算装置２０は、自らのメモリに漏水位置を記憶し、または、自らのディスプレイに漏水位置を表示する。
【００１８】
図２には、コンセントレータ１０、センサ装置１２−１および１２−２が実行する処理を示すシーケンスチャートが示されている。コンセントレータ１０は、各センサ装置１２に測定タイミング信号を送信する（Ｓ１０１）。測定タイミング信号は、振動波形データを取得するタイミングを各センサ装置１２に伝える信号である。
【００１９】
各センサ装置１２は、測定タイミング信号に従うタイミングで、水道管１４の検知点で音響振動を検出し振動波形データを生成する（Ｓ１０２，Ｓ１０３）。各センサ装置１２は振動波形データを記憶し、振動波形データの生成および記憶が完了した旨を示す完了通知信号をコンセントレータ１０に送信する（Ｓ１０４）。
【００２０】
コンセントレータ１０は、完了通知信号を受信すると、各センサ装置１２にデータ要求信号を送信する（Ｓ１０５，Ｓ１０７）。各センサ装置１２は、データ要求信号に応答して、コンセントレータ１０に振動波形データを送信する（Ｓ１０６，Ｓ１０８）。コンセントレータ１０は、各センサ装置１２から送信された振動波形データに基づいて、以下に説明する測定データを生成し（Ｓ１０９）、通信網１８を介して演算装置２０に送信する。
【００２１】
測定データは、演算装置２０が漏水位置を求める元となるデータであり、センサ装置１２−１から取得した振動波形データが示す振動波形と、センサ装置１２−２から取得した振動波形データが示す振動波形をシフト時間だけシフトした時間波形との相互相関値に基づくデータである。
【００２２】
以下、測定データを生成する処理について説明する。図３には、センサ装置１２−１による振動波形データが示す振動波形ｆ＝ｙ（ｔ）と、センサ装置１２−２による振動波形データが示す振動波形ｇ＝ｘ（ｔ）が示されている。また、振動波形ｇ＝ｘ（ｔ）を時間軸上でシフト時間τだけ過去の方向にシフトした時間波形ｇ＝ｘ（ｔ＋τ）および振動波形ｇ＝ｘ（ｔ）を時間軸上でシフト時間τｍだけ過去の方向にシフトした時間波形ｇ＝ｘ（ｔ＋τｍ）が併せて示されている。振動波形ｆ＝ｙ（ｔ）および振動波形ｙ＝ｘ（ｔ＋τ）の相互相関値Ｍ（τ）は、次の（数１）および（数２）のように表される。
【００２３】
【数１】
2019124613.tif 000003
【００２４】
【数２】
2019124613.tif 000004
【００２５】
ただし、符号「−」が上に付されたｘおよびｙは、それぞれ、ｘ（ｔ）およびｙ（ｔ）の時間範囲Ｔｓ≦ｔ≦Ｔｅにおける時間平均値である。相互相関値Ｍ（τ）は、時間波形ｘ（ｔ＋τ）と時間波形ｙ（ｔ）とが近似している度合いを示す。したがって、時間波形ｘ（ｔ＋τ）のピークが現れる時間と時間波形（ｔ）のピークが現れる時間とが一致するときのシフト時間τに相互相関値Ｍ（τ）は極大値となる。
【００２６】
図３に示されている例では、シフト時間τを増加させながら相互相関値Ｍ（τ）を順次求めた場合、シフト時間がτ＝τｍのときに、時間波形ｆ＝ｙ（ｔ）のピークが現れる時間と、時間波形ｇ＝ｘ（ｔ＋τ）のピークが現れる時間とが一致し、相互相関値Ｍ（τ）が極大値となる。
【００２７】
本実施形態におけるコンセントレータ１０が取得する振動波形データは、ｘ（ｔ）およびｙ（ｔ）が時間軸上で離散化されたサンプル値群によって構成され、ｘ（ｔ）およびｙ（ｔ）は、それぞれ（数３）および（数４）のように表される。ただし、δ（ｔ）はｔ＝０で１となり、その他の時間で０となる関数である。また、整数ｉは時間変数としての整数であり、ｘｉは時間ｔｉにおけるサンプル値である。Ｎｅは、振動波形データのサンプル値の数である。
【００２８】
【数３】
2019124613.tif 000005
【００２９】
【数４】
2019124613.tif 000006
【００３０】
コンセントレータ１０は、（数５）〜（数７）に基づいて時間波形ｘ（ｔ）と時間波形ｙ（ｔ）との相互相関値Ｍ（ｔｊ）を求める。（数５）〜（数７）はｘ（ｔ）およびｙ（ｔ）を離散化した場合の（数１）および（数２）に相当し、（数１）および（数２）と同様の技術的意義を有する。ただし、シフト時間に対応するシフト量ｊはｊｓ〜ｊｅの整数であり、−ｎ〜ｎの範囲内の整数である。整数ｉ＝ｋｓ〜ｋｅは相互相関演算が行われる範囲を示す整数である。加算合計範囲の記載がないシグマ記号は、ｉ＝ｋｓ〜ｋｅに亘る加算合計値を示す。
【００３１】
【数５】
2019124613.tif 000007
【００３２】
【数６】
2019124613.tif 000008
【００３３】
【数７】
2019124613.tif 000009
【００３４】
相互相関値Ｍ（ｔｊ）は、時間波形ｙ（ｔ）と、時間波形ｘ（ｔ）をシフト時間ｊ・Δだけ過去に向かう方向にシフトした時間波形ｘ（ｔ＋ｊ・Δ）とが近似している度合いを示す。ただし、Δはｘ（ｔ）およびｙ（ｔ）のサンプル時間間隔である。
【００３５】
図４（ａ）には、シフト量ｊが−５０００である場合の相互相関値Ｍ（ｔ−５０００）を求める演算が概念的に示されており、図４（ｂ）には、シフト量ｊが−４９９９である場合の相互相関値Ｍ（ｔ−４９９９）を求める演算が概念的に示されている。これらの図の例では、ｎ＝５０００，ｋｓ＝１，ｋｅ＝３００００である。ただし、ｋｅ＝３００００については図面に現れていない。
【００３６】
すなわち、図４（ａ）にはｎ＝５０００，ｊ＝−５０００である場合に、ｘｉ＋ｎ＋ｊ（ｉ＋ｎ＋ｊ＝１，２，３・・・・３００００）に対し、ｙｉ＋ｎ（ｉ＋ｎ＝５００１，５００２，５００３・・・・３５０００）がそれぞれ（数５）における演算の対象となり、相互相関値Ｍ（ｔ−５０００）が求められることが、各ｘｉを表す長方形と各ｙｉを表す長方形とを結ぶ直線によって示されている。図４（ｂ）にはｎ＝５０００，ｊ＝−４９９９である場合に、ｘｉ＋ｎ＋ｊ（ｉ＝２，３，４・・・・３０００１）に対し、ｙｉ＋ｎ（ｉ＋ｎ＝５００１，５００２，５００３・・・・３５０００）がそれぞれ（数５）における演算の対象となり、相互相関値Ｍ（ｔ−４９９９）が求められることが、各ｘｉを表す長方形と各ｙｉを表す長方形とを結ぶ直線によって示されている。
【００３７】
コンセントレータ１０は、ｊ＝ｊｓ〜ｊｅの相互相関値Ｍ（ｔｊ）について大きいものから順にＮ個の極大値を探索する。そして、Ｎ個の極大値のそれぞれに対応するシフト時間（極大点シフト時間）、およびＮ個の極大値のそれぞれの偏差値を求める。偏差値は、相互相関値Ｍ（ｊｓ）〜Ｍ（ｊｅ）を母集団とする偏差値であり、この母集団において値が際立つ度合いを示す。
【００３８】
図５には、シフト時間に対する相互相関値の特性の例が示されている。横軸はシフト時間を示し、縦軸は相互相関値を示す。（表１）の右側の欄には、図５で例示された相互相関値について、大きいものから順にＮ＝４個の極大値を探索し、これら４個の極大値に対するシフト時間、相互相関値の極大値および偏差値を求めた例が示されている。
【００３９】
【表１】
2019124613.tif 000010
【００４０】
なお、左半分の欄には、漏水が生じていない場合におけるシフト時間、相互相関値の極大値および偏差値が示されている。
【００４１】
コンセントレータ１０は、相互相関値が極大値となるときのシフト時間、相互相関値の極大値、および、その極大値に対応する偏差値の組を１組の相互相関極大点データとして、Ｎ組の相互相関極大点データを含む測定データを通信網１８を介して演算装置２０に送信する。
【００４２】
演算装置２０は、測定データに含まれるＮ組の相互相関極大点データのうち、相互相関値の極大値が所定の閾値を超えるものがないときは漏水が発生していないと判断する。一方、測定データに含まれるＮ組の相互相関極大点データのうち、相互相関値の極大値が所定の閾値を超えるものがあるときは漏水が発生したものと判断し、次のような処理を実行する。
【００４３】
演算装置２０は、Ｎ組の相互相関極大点データのうち、偏差値が最大である１組の相互相関極大点データにおけるシフト時間を測定遅延時間ｔｍとして求める。そして、（数８）に基づいてセンサ装置１２−１の位置を基準とした漏水位置ｄｗを求める。（数８）のｖは音響振動の伝搬速度を示す定数であり、予め実験等によって求められる。この漏水位置ｄｗは、センサ装置１２−１からセンサ装置１２−２に向かう方向に距離ｄｗだけ離れた位置を示す。
【００４４】
【数８】
2019124613.tif 000011
【００４５】
（数８）は、距離ｄｗと漏水位置からセンサ装置１２−２までの距離Ｌとの間に、Ｓ＝ｄｗ＋Ｌ、およびｔｍ・ｖ＝ｄｗ―Ｌの関係が成立することから導かれる数式である。  測定遅延時間ｔｍは負の値を取り得るため、（数８）で求められる漏水位置ｄｗもまた負の値を取り得る。漏水位置ｄｗが負である場合、漏水位置ｄｗはセンサ装置１２−１からセンサ装置１２−２に向かう方向とは反対の方向に距離｜ｄｗ｜だけ離れた位置を示す。
【００４６】
演算装置２０は、漏水位置ｄｗを自らのメモリに記憶する。また、ディスプレイなどの表示装置を備える場合には、演算装置２０は漏水位置を表示してもよい。
【００４７】
このように、本実施形態に係る漏水検出システムでは、コンセントレータ１０が各センサ装置１２から振動波形データを取得して測定データを求める。コンセントレータ１０は測定データを通信網１８を介して演算装置２０に送信し、演算装置２０は測定データに基づいて漏水が発生しているか否かを判定する。演算装置２０は、漏水が発生していると判定したときは、漏水位置ｄｗを求める。コンセントレータ１０から演算装置２０に送信される情報は、相互相関値を求めるための振動波形データではなく情報量が削減されたＮ組の相互相関極大点データである。これによって、コンセントレータ１０から演算装置２０に伝送される情報量が低減される。
【００４８】
なお、コンセントレータ１０は、（数５）〜（数７）に示される演算における整数ｉ＝ｋｓ〜ｋｅに関する加算合計を、総ての時間変数ｉ＝ｋｓ〜ｋｅについて行うのではなく、２以上の整数をｐとして、ｐ−１おきに時間変数を増加させて加算合計を行ってもよい。すなわち、コンセントレータ１０は、ｐ（ｉ−ｋｓ）＋ｋｓを時間変数として加算合計を行ってもよい。このとき（数５）〜（数７）は、それぞれ、（数９）〜（数１１）のように変形される。（数９）〜（数１１）は、それぞれ、（数５）〜（数７）におけるｉを、ｐ（ｉ−ｋｓ）＋ｋｓに置き換えたものである。
【００４９】
【数９】
2019124613.tif 000012
【００５０】
【数１０】
2019124613.tif 000013
【００５１】
【数１１】
2019124613.tif 000014
【００５２】
加算合計範囲の上限を表す整数ｋｄは、ｐ（ｉ−ｋｓ）＋ｋｓがｋｅを超えない最大値Ｎｍａｘとなるときにおけるｉの値であり、ｋｄ＝（Ｎｍａｘ−ｋｓ）／ｐ＋ｋｓである。  図６（ａ）および（ｂ）には、（数９）〜（数１１）で表される時間波形ｘ（ｔ）およびｙ（ｔ）が概念的に示されている。図６（ａ）には、シフト量ｊが−５０００である場合の相互相関値Ｍ（ｔ−５０００）を求める演算が概念的に示されており、図６（ｂ）には、シフト量ｊが−４９９９である場合の相互相関値Ｍ（ｔ−４９９９）を求める演算が概念的に示されている。これらの図の例では、ｎ＝５０００，ｐ＝２，ｋｓ＝１，ｋｅ＝３００００である。ただし、ｋｅ＝３００００については図面に現れていない。
【００５３】
すなわち、図６（ａ）にはｎ＝５０００，ｊ＝−５０００である場合に、ｘｐ（ｉ−ｋｓ）＋ｋｓ＋ｎ＋ｊ（ｐ（ｉ−ｋｓ）＋ｋｓ＋ｎ＋ｊ＝１，３，５・・・・２９９９９）に対し、ｙｐ（ｉ−ｋｓ）＋ｋｓ＋ｎ（ｐ（ｉ−ｋｓ）＋ｋｓ＋ｎ＝５００１，５００３，５００５・・・・３４９９９）がそれぞれ（数９）における演算の対象となり、相互相関値Ｍ（ｔ−５０００）が求められることが示されている。図６（ｂ）にはｎ＝５０００，ｊ＝−４９９９である場合に、ｘｐ（ｉ−ｋｓ）＋ｋｓ＋ｎ＋ｊ（ｐ（ｉ−ｋｓ）＋ｋｓ＋ｎ＋ｊ＝２，４，６・・・・３００００）に対し、ｙｐ（ｉ−ｋｓ）＋ｋｓ＋ｎ（ｐ（ｉ−ｋｓ）＋ｋｓ＋ｎ＝５００１，５００３，５００５・・・・３４９９９）がそれぞれ（数９）における演算の対象となり、相互相関値Ｍ（ｔ−４９９９）が求められることが示されている。
【００５４】
このように、時間変数ｉを間引いて増加させたｘｉの各値（ｉをｐ−１おきに増加させたｘｉの各値）に対して、ｙ（ｔ）のサンプル値ｙｉ＋ｊを相互相関演算の相手とすることで、時間変数ｉを１ずつ増加させた場合に比べて、相互相関値Ｍ（ｔｊ）を求める際の演算量が少なくなる。
【００５５】
なお、上記では、相互相関極大点データに、相互相関値が極大値となるときのシフト時間、相互相関値の極大値、および、その極大値に対応する偏差値が含まれる例について説明した。このように、相互相関値の極大値が母集団から際立っている程度を示す際立ち値として偏差値が用いられている。際立ち値には、偏差値とは異なる値が用いられてもよい。例えば、相互相関値Ｍ（ｊｓ）〜Ｍ（ｊｅ）を母集団として、相互相関値Ｍ（ｊｓ）〜Ｍ（ｊｅ）の大きさの傾向を表す傾向値を求め、相互相関値の極大値から傾向値を減算した値が偏差値の代わりに用いられてもよい。傾向値としては、例えば、相互相関値Ｍ（ｊｓ）〜Ｍ（ｊｅ）の平均値、中央値、メジアン等がある。
【００５６】
図７には、コンセントレータ１０のハードウエアの例が示されている。コンセントレータ１０は、制御部２２、無線部３０、および通信インターフェース３２を備える。無線部３０は、センサ装置から無線送信された信号を受信し、制御部２２に出力する。また、無線部３０は、制御部２２から出力された信号を各センサ装置に無線送信する。通信インターフェース３２は、通信網１８から送信された信号を受信し、制御部２２に出力する。また、通信インターフェース３２は、制御部２２から出力された信号を通信網１８に送信する。
【００５７】
制御部２２は、予め記憶されたプログラムに従って動作するプロセッサを含んでもよい。この場合、制御部２２は、プログラムを実行することで通信処理部２４、相関処理部２６および測定部２８を構成する。これらの構成要素は、ディジタル回路によって個別に構成されてもよい。
【００５８】
通信処理部２４は無線部３０を介して、図２に示された測定タイミング信号の送信（Ｓ１０１）、完了通知信号の受信（Ｓ１０４）、データ要求信号の送信（Ｓ１０５，Ｓ１０７）、振動波形データの受信（Ｓ１０６，Ｓ１０８）等を行う。
【００５９】
相関処理部２６は、波形取得部としての通信処理部２４が、無線部３０を介して受信した各振動波形データに対し、上述の相互相関処理を実行する。測定部２８は、上述の測定データを求める処理を実行する。
【００６０】
このように漏水検出装置としてのコンセントレータ１０は、水道管における２つの検知点のそれぞれに配置されたセンサから、振動波形情報を取得する波形取得部（通信処理部２４）と、各センサから取得された振動波形データ（振動波形情報）についての相関処理であって、一方の波形をシフト時間だけ時間軸上でシフトした場合における当該シフト時間に対する相互相関値を求める相関処理を実行する相関処理部２６と、シフト時間と相互相関値との関係に基づく測定データ（測定情報）であって、漏水位置を求める元となる測定データを求める測定部２８とを備える。
【００６１】
一つの実施形態における振動波形データは、時間軸上に配列された複数のサンプル値によって時間波形を表す情報である。相互相関処理は、所定サンプル数おきのサンプル値に対して相関演算を施す処理であってもよい。
【００６２】
また、コンセントレータ１０は、測定データを送信する送信部としての通信処理部２４を備え、この通信処理部２４は、無線部３０を介して、測定データに基づいて漏水位置を求める漏水位置検出装置としての演算装置２０に測定データを送信する。
【００６３】
一つの実施形態における測定データは、複数組の相互相関極大点データを含み、各相互相関極大点データは、シフト時間と相互相関値との関係におけるピークについての、シフト時間、相互相関値のピーク値、および際立ち値を含む。
【００６４】
コンセントレータ１０によれば、水道管における漏水位置を検出する処理において、センサと演算装置２０との間で伝送される情報が低減される。
【００６５】
上記では、水道管１４における２つの検知点のそれぞれにセンサ装置が配置された実施形態が示された。センサ装置は、水道管１４における３つ以上の検知点のそれぞれに設けられてもよい。この場合、コンセントレータ１０および演算装置２０は、複数のセンサ装置から選択される２つのセンサ装置から取得される振動波形データに基づいて、その２つのセンサ装置の位置によって定義される漏水位置を求める。
【符号の説明】
【００６６】
１０  コンセントレータ、１２−１，１２−２  センサ装置、１４  水道管、１８  通信網、２０  演算装置、２２  制御部、２４  通信処理部、２６  相関処理部、２８  測定部、３０  無線部、３２  通信インターフェース。
【図１】
2019124613.tif 000015
【図２】
2019124613.tif 000016
【図３】
2019124613.tif 000017
【図４】
2019124613.tif 000018
【図５】
2019124613.tif 000019
【図６】
2019124613.tif 000020
【図７】
2019124613.tif 000021