(57)【特許請求の範囲】
【請求項１】
配水管路網を構成する管路の接続構成を示す情報に基づいて、前記配水管路網のグラフ構造を表す接続行列を生成する接続行列生成部と、
前記接続行列に基づいて、前記配水管路網を構成する管路を前記配水管路網において開管路を構成する木となる管路と、前記木とともに閉管路を構成する補木となる管路とに分類し、前記補木となる管路を抽出する補木抽出部と、
前記抽出された補木となる管路の両端の節点を水圧センサの配置位置として選択する選択部と、
を備えるセンサ配置位置選択装置。
【請求項２】
前記補木抽出部は、水圧センサの配置位置に選択される節点に関して予め設定された目標条件を満たすように、前記補木となる管路を抽出する、
請求項１に記載のセンサ配置位置選択装置。
【請求項３】
前記目標条件は、各節点における水圧センサの設置可否を示す条件であり、前記補木抽出部は、水圧センサの設置が可能な節点を持つ管路を補木として抽出する、
請求項２に記載のセンサ配置位置選択装置。
【請求項４】
前記木となる管路で構成される開管路のグラフ構造に基づいて、前記開管路を幹と枝とに分類し、前記枝なる管路が存在する場合、前記枝となる管路を前記幹となる管路に集約する集約部を備える、請求項１から３のいずれか一項に記載のセンサ配置位置選択装置。
【請求項５】
前記集約部は、口径の小さい管路を他の管路に集約する、
請求項４に記載のセンサ配置位置選択装置。
【請求項６】
前記集約部は、前記配置位置選択部によって選択される配置位置が既存の水圧センサの配置位置となる程度に前記管路の集約を行う、
請求項４又は５に記載のセンサ配置位置選択装置。
【請求項７】
前記補木抽出部は、前記補木となる管路の組み合わせを複数取得し、前記複数の組み合わせのうち、前記補木となる管路の数が最も少ない組み合わせを選択して前記補木となる管路を抽出する、
請求項１から６のいずれか一項に記載のセンサ配置位置選択装置。
【請求項８】
前記選択部は、前記補木となる管路の両端の節点から既存の水圧センサの配置位置に近い配置位置を選択する、
請求項１から７のいずれか一項に記載のセンサ配置位置選択装置。
【請求項９】
前記配水管路網における節点の接続関係を含む配水管路網情報を記憶する配水管路網情報記憶部と、
前記配水管路網情報記憶部により記憶された配水管路網情報よりも新しい情報が存在することを示す情報を取得したときに、前記配水管路網情報を更新すべき旨の情報を通知する更新通知部と、
を備える、請求項１から８のいずれか一項に記載のセンサ配置位置選択装置。
【請求項１０】
配水管路網を構成する管路の接続構成を示す情報に基づいて、前記配水管路網のグラフ構造を表す接続行列を生成する接続行列生成部と、
前記接続行列に基づいて、前記配水管路網を構成する管路を前記配水管路網において開管路を構成する木となる管路と、前記木とともに閉管路を構成する補木となる管路とに分類し、前記補木となる管路を抽出する補木抽出部と、
前記抽出された補木となる管路の両端の節点を水圧センサの配置位置として選択する選択部と、
を備えるセンサ配置位置選択装置として、
コンピュータを機能させるためのコンピュータプログラム。
【請求項１１】
請求項１から９のいずれか一項に記載されたセンサ配置位置選択装置によって選択された配置位置に設置された水圧センサによって計測された、配水管路網において補木となる管路の水圧を示す情報と、配水管路網の各節点を通る開管路の始点から終点までの差圧と、前記配水管路網内の漏水量の総量と、各節点への流入流量と、各節点での水の使用量と、に基づいて各節点における漏水係数を決定変数とする最適化問題を解くことで、各節点における漏水係数を決定する漏水係数最適化部と、
前記決定された漏水係数と、前記配水管路網における節点間の接続関係を反映した流体の運動方程式に基づいて節点ごとの漏水量を推定する漏水量推定部と、
を備える漏水量推定装置。
【請求項１２】
前記漏水量推定部により推定された漏水量に基づく画像を表示する表示部を備える、
請求項１１に記載の漏水量推定装置。
【請求項１３】
前記表示部は、前記漏水量推定部により推定された漏水量に応じた色彩で、前記節点を表すシンボルを描画した画像を表示する、
請求項１１又は１２に記載の漏水量推定装置。
【請求項１４】
管路毎の埋設条件を含む埋設条件情報を記憶する埋設条件情報記憶部と、
前記漏水量推定部により推定された節点毎の漏水量および前記埋設条件情報に基づいて、管路を更新する際の優先順位を算出する更新優先順位算出部と、
を備える請求項１１から１３のいずれか一項に記載の漏水量推定装置。
【請求項１５】
一部又は全部の節点における水の使用量を示す情報に水道メータにより取得された検針データを用いる、
請求項１１から１４のいずれか一項に記載の漏水量推定装置。
【請求項１６】
請求項１から９のいずれか一項に記載されたセンサ配置位置選択装置によって選択された配置位置に設置された水圧センサによって計測された、配水管路網において補木となる管路の水圧を示す情報と、配水管路網の各節点を通る開管路の始点から終点までの差圧と、前記配水管路網内の漏水量の総量と、各節点への流入流量と、各節点での水の使用量と、に基づいて各節点における漏水係数を決定変数とする最適化問題を解くことで、各節点における漏水係数を決定する漏水係数最適化部と、
前記決定された漏水係数と、前記配水管路網における節点間の接続関係を反映した流体の運動方程式に基づいて節点ごとの漏水量を推定する漏水量推定部と、
を備える漏水量推定装置として、
コンピュータを機能させるためのコンピュータプログラム。
【請求項１７】
配水管路網を構成する管路の接続構成を示す情報に基づいて、前記配水管路網のグラフ構造を表す接続行列を生成する接続行列生成部と、
前記接続行列に基づいて、前記配水管路網を構成する管路を前記配水管路網において開管路を構成する木となる管路と、前記木とともに閉管路を構成する補木となる管路とに分類し、前記補木となる管路を抽出する補木抽出部と、
前記抽出された補木となる管路の両端の節点を水圧センサの配置位置として選択する選択部と、
前記選択部によって選択された配置位置に設置された水圧センサによって計測された、配水管路網において補木となる管路の水圧を示す情報と、配水管路網の各節点を通る開管路の始点から終点までの差圧と、前記配水管路網内の漏水量の総量と、各節点への流入流量と、各節点での水の使用量と、に基づいて各節点における漏水係数を決定変数とする最適化問題を解くことで、各節点における漏水係数を決定する漏水係数最適化部と、
前記決定された漏水係数と、前記配水管路網における節点間の接続関係を反映した流体の運動方程式に基づいて節点ごとの漏水量を推定する漏水量推定部と、
を備える漏水診断システム。
【請求項１８】
配水管路網を構成する管路の接続構成を示す情報に基づいて、前記配水管路網のグラフ構造を表す接続行列を生成する接続行列生成ステップと、
前記接続行列に基づいて、前記配水管路網を構成する管路を前記配水管路網において開管路を構成する木となる管路と、前記木とともに閉管路を構成する補木となる管路とに分類し、前記補木となる管路を抽出する補木抽出ステップと、
前記抽出された補木となる管路の両端の節点を水圧センサの配置位置として選択する選択ステップと、
前記選択ステップにおいて選択された配置位置に設置された水圧センサによって計測された、配水管路網において補木となる管路の水圧を示す情報と、配水管路網の各節点を通る開管路の始点から終点までの差圧と、前記配水管路網内の漏水量の総量と、各節点への流入流量と、各節点での水の使用量と、に基づいて各節点における漏水係数を決定変数とする最適化問題を解くことで、各節点における漏水係数を決定する漏水係数最適化ステップと、
前記決定された漏水係数と、前記配水管路網における節点間の接続関係を反映した流体の運動方程式に基づいて節点ごとの漏水量を推定する漏水量推定ステップと、
を有する漏水診断方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、センサ配置位置選択装置、漏水量推定装置、漏水診断システム、漏水診断方法及びコンピュータプログラムに関する。
【背景技術】
【０００２】
一般に、配水管路網における漏水の調査には、漏水の有無を調査する一次調査と、漏水箇所を特定する二次調査とがある。一次調査は調査員によって定期的に行われる調査であり、調査員は音聴棒等を用いて配水管路網における漏水の有無を調査する。二次調査は、一次調査の結果、漏水が発生している可能性が高いと判断された箇所について行われる調査であり、相関式漏水探査機を用いて漏水箇所が特定される。一次調査は、対象区域について均等に行われ、どの地域を重点的に調査するかについては考慮されていないのが現状である。
【０００３】
一方で、環境問題に対する意識の高まりを背景に、水道スマートメータの導入が検討されている。水道スマートメータは、需要家の使用水量を計測することによって、配水を効率的に行うことを目的として各戸に設置されるものである。水道スマートメータが普及すれば、水圧センサと組み合わせることによって、配水管路網の漏水診断を行う漏水診断システムを構成することも可能になると考えられる。
【０００４】
しかしながら、現状では漏水診断に必要な水圧センサを効率的に配置する基準がないため、水圧センサを過剰に配置してしまう可能性があり、管理コストの増大を招く可能性が懸念される。そのため、漏水診断を行うための水圧センサを効率的に配置する基準の策定が望まれている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１３−１７８２０７号公報
【特許文献２】特開２００６−２８５３８９号公報
【特許文献３】特開２０１０−４８０５８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
本発明が解決しようとする課題は、水圧センサの効率的な配置位置を選択することができるセンサ配置位置選択装置、漏水量推定装置、漏水診断システム、漏水診断方法及びコンピュータプログラムを提供することである。
【課題を解決するための手段】
【０００７】
実施形態のセンサ配置位置選択装置は、接続行列生成部と、補木抽出部と、選択部とを持つ。接続行列生成部は、配水管路網を構成する管路の接続構成を示す情報に基づいて、前記配水管路網のグラフ構造を表す接続行列を生成する。補木抽出部は、前記接続行列に基づいて、前記配水管路網を構成する管路を前記配水管路網において開管路を構成する木となる管路と、前記木とともに閉管路を構成する補木となる管路とに分類し、前記補木となる管路を抽出する。選択部は、前記抽出された補木となる管路の両端の節点を水圧センサの配置位置として選択する。
【図面の簡単な説明】
【０００８】
【図１】第１の実施形態に係る配水システムの概要とセンサ配置位置選択装置１の構成を示す図。
【図２】センサ配置位置選択装置１の機能構成を示す機能ブロック図。
【図３】実施形態の配水管路網情報記憶部１１に格納される情報の一例を示す図。
【図４】配水管路網のモデルの一例を示す図。
【図５】配水管路網情報の一例を示す図。
【図６】配水管路網ＰＮが１本の管路とみなされる一例を示す図。
【図７】変形例のセンサ配置位置選択装置１ａの機能構成を示す機能ブロック図。
【図８】グラフ構造を集約する第１の方法の概要を示す概略図。
【図９】グラフ構造を集約する第２の方法の概要を示す概略図。
【図１０】異なる配置位置が選択される場合の具体例を示す図。
【図１１】第２の実施形態の漏水量推定装置２の機能構成を示す機能ブロック図。
【図１２】ある時刻の仮想管路における有効水圧差と各管路の差圧との関係の具体例を示す。
【図１３】制約付き非線形最適化問題の定式化の概略を示すブロック図。
【図１４】漏水量の推定のシミュレーション結果の一例を示す図。
【図１５】出力部２２により漏水量が表示される画面の一例を示す図。
【発明を実施するための形態】
【０００９】
以下、実施形態のセンサ配置位置選択装置、漏水量推定装置、漏水診断システム、漏水診断方法及びコンピュータプログラムを、図面を参照して説明する。
【００１０】
（第１の実施形態）
図１は、第１の実施形態に係る配水システムの概要を示す概略図である。この配水システムでは、配水池７０に蓄えられた水（浄水）が、ポンプやバルブ等によって家庭や事業所等を含む配水管路網ＰＮに供給される。配水管路網ＰＮに流入する水の量は、流量センサ８０によって検出される。
【００１１】
配水管路網ＰＮには、図中、白丸で表される複数の節点が便宜的に設定されている。図中、複数の節点にはスマートメータが取り付けられている。スマートメータは、例えば、水が提供される家庭や事業所等に取り付けられ、家庭や事業所等における水の使用量を３０分ないし１時間程度ごとに検出する。各節点から下に向く矢印は、需要家による水の使用や漏水等によって水が流出していることを表している。
【００１２】
センサ配置位置選択装置は、このような配水管路網の各節点から、節点ごとの漏水量を推定することが可能となるような水圧センサの配置位置を選択する。
【００１３】
図２は、第１の実施形態のセンサ配置位置選択装置１の機能構成を示す機能ブロック図である。
センサ配置位置選択装置１は、バスで接続されたＣＰＵ（Central Processing Unit）やメモリや補助記憶装置などを備え、センサ配置位置選択プログラムを実行する。センサ配置位置選択装置１は、センサ配置位置選択プログラムの実行によって配水管路網情報記憶部１１、入力部１２、出力部１３、接続行列生成部１４、補木抽出部１５及び配置位置選択部１６を備える装置として機能する。なお、センサ配置位置選択装置１の各機能の全て又は一部は、ＡＳＩＣ（Application Specific Integrated Circuit）やＰＬＤ（Programmable Logic Device）やＦＰＧＡ（Field Programmable Gate Array）等のハードウェアを用いて実現されてもよい。センサ配置位置選択プログラムは、コンピュータ読み取り可能な記録媒体に記録されてもよい。コンピュータ読み取り可能な記録媒体とは、例えばフレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置である。センサ配置位置選択プログラムは、電気通信回線を介して送信されてもよい。
【００１４】
配水管路網情報記憶部１１は、磁気ハードディスク装置や半導体記憶装置などの記憶装置を用いて構成される。
【００１５】
入力部１２は、例えば、キーボード、マウス、タッチパネル、スイッチ等の入力デバイスを用いて構成されユーザの操作の入力を受け付ける。
【００１６】
出力部１３（表示部）は、例えば、ＬＣＤ（Liquid  Crystal  Display）や有機ＥＬ（Electroluminescence）等の表示装置および表示制御部、プリンタ、スピーカ等を含んで構成され、選択された水圧センサの配置位置に関する情報を出力する。
【００１７】
接続行列生成部１４は、配水管路網情報記憶部１１に格納された情報に基づいて、配水管路網にグラフ理論を適用するための接続行列を生成する。接続行列とは、配水管路網のグラフ構造（接続構成）を示す行列である。
【００１８】
補木抽出部１５は、接続行列生成部１４により生成された接続行列に基づいて、配水管路網を構成する管路から補木となる管路を抽出する。補木とは、グラフ構造を構成する枝の種別の１つである。グラフ構造を構成する枝は木又は補木に分類され、閉回路を作らない枝集合のことを木（Tree）といい、木に含まれない枝集合を補木（Cotree）という。具体的には、補木抽出部１５は、配水管路網を構成する管路を、木となる管路と補木となる管路とに分類し、補木に分類された管路を抽出する。
【００１９】
配置位置選択部１６は、補木抽出部１５によって抽出された補木となる管路の両端の節点を水圧センサの配置位置として選択し、選択された配置位置を示す情報を出力部１３に出力する。
【００２０】
以下で説明する管路網モデルに基づけば、各節点における漏水量の推定に必要な水圧センサは、上記の補木となる管路の両端の節点に設置されればよいことが分かる。
【００２１】
［管路網モデルの構築］
配水管路網を流れる非圧縮流体の運動方程式は、式（１）で表される。式（１）において、ｉ、ｊは節点の番号である。ｖｉｊは管路ｉｊ（節点ｉと節点ｊを結ぶ管路）における水の流速である。Ｌｉｊは管路ｉｊの長さ［ｍ］である。ρは水密度［ｋｇ／ｍ３］である。Ｈｉは節点ｉの標高である。Ｄｉｊは管路ｉｊの口径［ｍ］である。λｉｊは管路ｉｊの管路摩擦抵抗である。
【００２２】
【数１】
0006334417.tif 000002
【００２３】
配水管路網情報記憶部１１には、節点間の接続関係の他、管路の長さ、口径、および管路摩擦抵抗と、節点の標高などの情報が格納されている。図３は、実施形態の配水管路網情報記憶部１１に格納される情報の一例を示す図である。図示するように、配水管路網情報記憶部１１には、節点数や管路数の他、節点毎の有効水頭［ｍ］や種別、設置標高［ｍ］などの情報が、節点番号に対応付けられて記述されている。また、配水管路網情報記憶部１１には、管路毎の始点および終点の節点番号（すなわち節点間の接続関係）、管長（長さ）、管路摩擦係数などの情報が、管路番号に対応付けられて記述されている。
【００２４】
式（１）を行列表現で表すと、式（２）となる。
【００２５】
【数２】
0006334417.tif 000003
【００２６】
式（２）において「．」はベクトルを成分毎に乗算したベクトルを表している。また、｜（ベクトル）｜は、ベクトルの各成分を絶対値に変換したベクトルを表している。式（２）における各行列やベクトルの成分は、式（３）〜（７）で表される。
【００２７】
【数３】
0006334417.tif 000004
【００２８】
【数４】
0006334417.tif 000005
【００２９】
【数５】
0006334417.tif 000006
【００３０】
【数６】
0006334417.tif 000007
【００３１】
【数７】
0006334417.tif 000008
【００３２】
また、各節点における質量保存式は、式（８）で表される。
【００３３】
【数８】
0006334417.tif 000009
【００３４】
式（８）のＡは接続行列である。Ｓは区間断面積［ｍ２］を対角がゼロの正方行列で表したものである。ｑは各節点ｉの水の使用量を要素とするベクトルである。ｌは各節点ｉにおける漏水量を要素とするベクトルである。行列Ｓの要素Ｓｉｊは、管路ｉｊの口径［ｍ］から、式（９）により導出される。
【００３５】
【数９】
0006334417.tif 000010
【００３６】
［接続行列の生成］
接続行列Ａとは、有向グラフに基づいて生成される行列であり、節点の数をｐ、管路の数をｋとした場合にｐ×ｋ行列として生成される。接続行列Ａの各要素Ａｐｋは、下記の定義によって決定される。
Ａｐｋ＝―１  ：管路ｋにおいて節点ｐが始点の場合
Ａｐｋ＝  １  ：管路ｋにおいて節点ｐが終点の場合
Ａｐｋ＝  ０  ：管路ｋにおいて節点ｐが始点でも終点でも無い場合
【００３７】
以下、配水管路網のモデルとして図４を参照し、質量保存式について説明する。図４は、図１とは異なる配水管路網ＰＮの一例を示す図である。図４における丸は節点を表し、各丸を結ぶ線は管路を表す。丸内の数字は節点番号を示す。また、各管路に付記された括弧内の数字は管路番号を示す。図５は、図４に示す配水管路網ＰＮに対応して配水管路網情報記憶部１１に格納される情報を例示した図である。
【００３８】
図４に示す配水管路網ＰＮの場合、接続行列Ａは、次の式（１０）で表される。
【００３９】
【数１０】
0006334417.tif 000011
【００４０】
前述したように、接続行列Ａの各行は節点に対応し、各列は管路に対応する。図４及び図５に示すモデルを式（８）に適用すると、次の式（１１）となる。式（１１）におけるｌｉは、節点毎の漏水量であり、未知数である。
【００４１】
【数１１】
0006334417.tif 000012
【００４２】
ここで、有効水圧差であるｄＰと各ノードの圧力値Ｐｉは、Ｐ１をゼロとすると式（１２）の関係にある。
【００４３】
【数１２】
0006334417.tif 000013
【００４４】
そして、式（２）、（１２）から式（１３）が得られる。
【００４５】
【数１３】
0006334417.tif 000014
【００４６】
式（１３）に基づけば、各管路を流れる水の流速を算出することができる。例えば、ニュートンラプソン法などを用いることで、式（１３）の非線形連立常微分方程式を解くことができる。式（１３）の非線形連立常微分方程式を解く手法には、その他の手法が用いられてもよい。
【００４７】
ここで、Ａ*は、接続行列Ａから、任意の節点に対応する行を取り除いて得られる既約接続行列である。以下では、接続行列Ａの１行目を取り除くこととする。式（１０）に示される接続行列Ａから得られる既約接続行列Ａ*は、式（１４）で表される。そして、行列Ａ*ｔは、既約接続行列Ａ*を、補木が前列にくるように入れ替えた場合の、後列のみの部分行列である。
【００４８】
【数１４】
0006334417.tif 000015
【００４９】
以下、行列Ａ*ｔについて説明する。グラフ構造が閉回路（ループ）を形成している場
合、ループを一巡すると圧力損失がゼロとなる法則を得ることができる。すなわち、ある始点となる節点からループを辿って始点に戻ると、ｄＰ＝０である。これを式で表現すると、式（１５）となる。
【００５０】
【数１５】
0006334417.tif 000016
【００５１】
式（１５）中、Ｂは基本閉路行列であり、各行ｉがグラフの基本閉路Ｆｉに、各列ｋがグラフの各枝（管路）に対応するｉ×ｋ行列である。基本閉路行列Ｂの各要素Ｂｉｋは、下記の定義によって決定される。
Ｂｉｋ＝―１  ：閉路Ｆｉが枝ｋを負の向きに含んでいる場合
Ｂｉｋ＝  １  ：閉路Ｆｉが枝ｋを正の向きに含んでいる場合
Ｂｉｋ＝  ０  ：閉路Ｆｉが枝ｋを含まない場合
【００５２】
図３および図４に示すモデルの場合、基本閉路行列Ｂは式（１６）で表される。
【００５３】
【数１６】
0006334417.tif 000017
【００５４】
このように、任意のひとつの補木と木の集合により一意に決定される閉路を基本閉路と呼ぶ。グラフの基本閉路の数は、補木の数と等しい。基本閉路行列Ｂは、既約接続行列Ａ*と基本閉路行列Ｂの間にある式（１７）の関係により求めることができる。式中、添え字のｃは補木を、ｔは木を表している。
【００５５】
【数１７】
0006334417.tif 000018
【００５６】
補木は、例えば深さ優先探索によって選択される。具体的には、管路網モデル構築部４０は、管路番号１の始点節点からスタートし、管路番号の終点を探索する。状態変数算出部４２は、この探索を管路番号Ｍまで繰り返し、その始点および終点をスタックに追加していき、スタック内で重複する節点番号が３回以上現れた時点の管路番号を補木とする。
【００５７】
図３および図４に示すモデルの場合、基本閉路行列Ｂは、補木である管路（３）、（５）に対応する列が前列になるように、列を入れ替えることにより、式（１８）に変形できる。基本閉路の定義から、補木に関する成分は単位行列となる。
【００５８】
【数１８】
0006334417.tif 000019
【００５９】
この結果を受けて、同じように列を入れ替えることで、式（１９）に示すように、行列Ａ*ｃ、Ａ*ｔが得られる。以上説明した理論に基づけば、上記の管路網モデルに応じた演算を行うことにより、各管路を流れる水の流速を算出することができる。
【００６０】
【数１９】
0006334417.tif 000020
【００６１】
［漏水量の推定］
このように、接続行列Ａの成分を木と補木の成分に分離することにより、漏水量の推定に係る演算を補木と木とで独立させて行うことが可能となる。支配方程式である式（２）、（３）及び（９）の３式をそれぞれグラフの木及び補木の部分に分けて書き下すと式（２０）、（２１）、（２２）及び（２３）となる。
【００６２】
【数２０】
0006334417.tif 000021
【００６３】
【数２１】
0006334417.tif 000022
【００６４】
【数２２】
0006334417.tif 000023
【００６５】
［水圧センサの配置位置の選択］
上記のように、グラフ理論を応用すれば、配水管路網を構成する管路を補木と木とに分類することができる。その結果、配水管路網の支配方程式も補木及び木の部分に分けて考えることが可能となる。すなわち、補木となる管路に関して圧力損失ｄＰｃが計測できれば、配水管路網は、木となる管路で構成される１本の管路とみなすことが可能となる。１本の管路であれば各節点での漏水量は全圧力損失に直接影響することになるため、１つの系として扱うことができる。そのため、補木となる管路の両端に水圧センサが配置されれば、全ての節点における漏水量を推定することが可能となる。配置位置選択部１６は、このような理由から、補木として抽出された管路の両端の節点を水圧センサの配置位置として選択する。
【００６６】
図６は、配水管路網ＰＮが１本の管路とみなされる一例を示す図である。図６の左図は、配水管路網ＰＮの４つのループに対して、黒四角で示される４つの管路が補木として抽出されたことを示している。この場合、配水管路網ＰＮは、図６の右図のような一本の管路と等価となる。そして、黒丸の節点が水圧センサの配置位置として選択される。
【００６７】
例えば、節点数がＮ、管路数がＭ、ループ数がＦであるグラフ構造を持つ配水管路網において、各節点での漏水量を算出するためには、従来、全ての節点において水圧センサが必要であった。すなわち、Ｎ個の水圧センサが必要であった。これに対して、第１の実施形態の配置位置選択装置１は、水圧センサが配置される位置として２×Ｆ−α＋β個の節点を選択する。このような、配置位置の選択によって、漏水診断に必要となる水圧センサの数を削減することが可能となる。αは補木となる管路において重複する節点の数を表し、βは、ループを構成する構成する経路上にある節点の数を表す。具体的には、図６の例の場合、節点１００がαに対応する節点である。また、節点１０１がβに対応する節点である。
【００６８】
このように構成された第１の実施形態のセンサ配置位置選択装置１は、配水管路網２の接続構成をグラフ構造として表す接続行列に基づいて、配水管路網２を構成する管路から、ループを構成する補木となる管路を抽出する。そして、センサ配置位置選択装置１は、補木として抽出された管路の両端を、水圧センサの設置位置として選択する。このような機能を備えることにより、センサ配置位置選択装置１は、複数の管路の節点から水圧センサの効率的な配置位置を選択することが可能となる。
【００６９】
以下、第１の実施形態のセンサ配置位置選択装置１の変形例について説明する。
【００７０】
センサ配置位置選択装置１は、管路網の構成に基づいて、配水管路網２のグラフ構造を集約するグラフ構造集約部１７を備えてもよい。図７は、変形例のセンサ配置位置選択装置１ａの機能構成を示す機能ブロック図である。グラフ構造集約部１７は、配水管路網情報記憶部１１に記憶された情報に基づいて、補木抽出部１５と同様の方法で、各管路を木と補木とに分類する。グラフ構造集約部１７は、木となる管路のグラフ構造について、以下に示すような集約を行う。
【００７１】
図８は、グラフ構造を集約する第１の方法の概要を示す概略図である。
図８は、グラフ構造を集約する第１の方法として、木となる管路のグラフ構造に基づいて管路を枝と幹とに分類し、枝となった管路を幹となった管路に集約する方法を示す。図８の左図は、集約前の管路の例を示す図である。図８の右図は、集約後の管路の例を示す図である。図８における各節点の大きさは、各節点における水の使用量を表している。この場合、例えば、グラフ構造集約部１７は、左図の４つの管路を右図の１つの管路に集約してもよい。グラフ構造集約部１７は、このような集約の結果を、配水管路網情報記憶部１１に記憶された情報に反映させる。
【００７２】
また、上記の管路の集約は、配置位置選択部１６の選択結果が、既存の水圧センサの配置位置となる程度に行われてもよい。このような集約が行われることによって、既存の水圧センサを用いた場合であっても、（漏水診断の解像度は荒くなるが）節点ごとの漏水量の推定を行うことができる。
【００７３】
図９は、グラフ構造を集約する第２の方法の概要を示す概略図である。
図９は、グラフ構造を集約する第２の方法として、口径の小さい管路を集約する方法を示す。図８と同様に、図９の左図は集約前の管路の例を示し、右図は集約後の管路の例を示す。この場合、グラフ構造集約部１７は、左図の破線に示された口径の小さい管路をそれ以外の管路に集約してもよい。この場合、グラフ構造集約部１７は、破線に示された管路は他の管路に集約するのではなく、単純に無視してもよい。グラフ構造集約部１７は、このような集約の結果を、配水管路網情報記憶部１１に記憶された情報に反映させる。
【００７４】
また、センサ配置位置選択装置１ａによって選択される配置位置は、抽出される補木の組み合わせによって異なる。図１０は、異なる配置位置が選択される場合の具体例を示す図である。図１０において、破線で示された管路が補木として抽出された管路である。図１０の左図と右図とは、それぞれ異なる組み合わせで補木が抽出されている。このとき、左図における水圧センサの数は７個、右図の水圧センサの数は８個となる。よって、左図の組み合わせを選択する方が、水圧センサの数を削減するという点で有利である。
【００７５】
なお、上述したグラフ構造の集約や補木の組み合わせの選択は、水圧センサの数を削減することを可能にする一方で、水圧センサによって取得される情報量が減少することになり、漏水診断の解像度が低下することになる。ここでいう漏水診断の解像度とは、漏水量を推定することが可能な単位の細かさを意味する。すなわち、水圧センサの数を減らすことと漏水診断の解像度とはトレードオフの関係にある。そのため、水圧センサの数の削減は、両者が適切なバランスとなるように行われるとよい。
【００７６】
また、補木抽出部１５は、水圧センサの配置位置に選択される節点に関して予め設定された目標条件を満たすように補木となる管路を抽出するように構成されてもよい。例えば、目標条件には、水圧センサの設置可否が設定されてもよい。この場合、補木抽出部１５は、補木となる管路が水圧センサの設置が可能な節点を持つ場合に当該管路を抽出する。
【００７７】
また、このような目標条件は、選択される水圧センサの配置位置について設定されてもよい。例えば、選択される水圧センサの配置位置が既存の水圧センサの配置位置となるように目標条件が設定されてもよい。この場合、配置位置選択部１６は、抽出された補木となる管路の両端の節点に近い位置に設置された既存の水圧センサの位置を配置位置として選択してもよい。
【００７８】
センサ配置位置選択装置は、ユーザに対して、配水管路網情報記憶部１１に記憶される配水管路網情報の更新を促すように構成されてもよい。例えば、センサ配置位置選択装置は、現在記憶されている配水管路網情報よりも新しい配水管路網情報が存在する場合、新しい配水管路網情報の存在を示す情報を外部から取得し、ユーザに更新を促す通知を行う更新通知部を備えてもよい。
【００７９】
（第２の実施形態）
第２の実施形態では、第１の実施形態のセンサ配置位置選択装置１によって選択された配置位置に設置された水圧センサと、配水管路網ＰＮへの水の供給量を計測する流量計と、各管路の節点に配置されたスマートメータと、を用いて各節点における漏水を診断する漏水量推定装置２について説明する。
【００８０】
図１１は、第２の実施形態の漏水量推定装置２の機能構成を示す機能ブロック図である。
漏水量推定装置２は、バスで接続されたＣＰＵやメモリや補助記憶装置などを備え、漏水診断プログラムを実行する。漏水量推定装置２は、漏水診断プログラムの実行によって入力部２１、出力部２２、計測データ取得部２３、管路網モデル構築部２４、漏水係数最適化部２５及び漏水量推定部２６を備える装置として機能する。なお、漏水量推定装置２の各機能の全て又は一部は、ＡＳＩＣやＰＬＤやＦＰＧＡ等のハードウェアを用いて実現されてもよい。漏水診断プログラムは、コンピュータ読み取り可能な記録媒体に記録されてもよい。コンピュータ読み取り可能な記録媒体とは、例えばフレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置である。漏水診断プログラムは、電気通信回線を介して送信されてもよい。
【００８１】
入力部２１は、例えば、キーボード、マウス、タッチパネル、スイッチ等の入力デバイスを用いて構成されユーザの操作の入力を受け付ける。
【００８２】
出力部２２は、例えば、ＬＣＤ（Liquid  Crystal  Display）や有機ＥＬ（Electroluminescence）等の表示装置および表示制御部、プリンタ、スピーカ等を含んで構成され、選択された水圧センサの配置位置に関する情報を出力する。
【００８３】
計測データ取得部２３は、配水管路網に設置された各種計測装置から、計測によって取得された計測データを取得する。
【００８４】
管路網モデル構築部２４は、第１の実施形態のセンサ配置位置選択装置１から配水管路網情報を取得する。管路網モデル構築部２４は、取得した配水管路網情報に基づいて、配水管路網の管路網モデルを構築する。管路網モデルの構築については、上述したとおりである。管路網モデル構築部２４は、構築した管路網モデルに関する情報を漏水係数最適化部２２に出力する。
【００８５】
漏水係数最適化部２５は、計測データ取得部２３によって取得された計測データに基づいて、各節点の漏水量の推定に必要となる漏水係数を決定する。漏水係数最適化部２５は、決定した漏水係数を漏水量推定部２６に出力する。
【００８６】
漏水量推定部２６は、漏水係数最適化部２５によって決定された漏水係数に基づいて、各節点における漏水量を算出する。漏水量推定部２６は、算出された各節点における漏水量を、漏水量の推定値として出力部２２に出力する。
【００８７】
［漏水係数の決定］
第１の実施形態でも説明したように、補木として抽出された管路の両端の水圧が計測できれば、Ｎ個の節点を持つ配水管路網２は、Ｎ−１の節点から水が引き抜かれる１本の管路（以下、「仮想管路」という。）と等価となる。そして、ある時刻の仮想管路における有効水圧差ｄＰと、各管路の差圧ｄＰｊとの関係は図１２の具体例のように表される。図１２に示されるように、１本の仮想管路であれば、各節点での漏水量と水圧差の関係が一意に決まるため、配水管路網の各節点における漏水量の推定は、制約付き非線形最適化問題として定式化することができる。以下、このような制約付き非線形最適化問題の定式化の方法について、２つの具体例を示す。
【００８８】
［第１の方法］
第１の方法は、全漏水量が各節点での漏水量の総和として表されることを制約とし、各時刻において計測された計測データに基づいて、仮想管路の始点から終点までの差圧が実測値に漸近するように漏水係数ｋｉを最適化する方法である。漏水係数ｋｉは、各節点における漏水量の推定値を決定する係数であり、水圧の１．１５乗に比例すると仮定すれば式（２４）で表される。
【００８９】
【数２３】
0006334417.tif 000024
【００９０】
この場合、第１の方法では、各節点における漏水量の推定は式（２５）のように定式化される。
【００９１】
【数２４】
0006334417.tif 000025
【００９２】
なお、仮想管路の始点から終点までの差圧ｄＰ（ｔ）及び漏水量の総量Ｌｔｏｔａｌは計測可能であり既知である。漏水量の総量Ｌｔｏｔａｌは式（２６）で表される。
【００９３】
【数２５】
0006334417.tif 000026
【００９４】
上記の最適化問題を解くことにより漏水係数ｋｉが求められる。漏水係数ｋｉが決まれば、式（２７）の非線形連立方程式を収束計算で求解することでＶｔが求められる。
【００９５】
【数２６】
0006334417.tif 000027
【００９６】
［第２の方法］
第２の方法は、第１の方法における全漏水量に関する制約条件を一部緩和して漏水係数ｋｉを最適化する方法である。この場合、第２の方法では、各節点における漏水量の推定は、例えば式（２８）のように定式化する。
【００９７】
【数２７】
0006334417.tif 000028
【００９８】
ここでは、最適化すべき目的関数をｆ１及びｆ２の和としたが、ｆ１及びｆ２の積など他の式で定式化されてもよい。
【００９９】
上述した定式化の方法をブロック図で表した場合、図１３のように表される。また、上記のように定式化された最適化問題を求解する手法には任意の手法が用いられてもよい。例えば、次の式（２９）に示されるニュートン法による探索によって、最適化問題の求解が行われてもよい。式（２９）において、Ｆは目的関数の行列表現であり、Ｊはヤコビ行列である。
【０１００】
【数２８】
0006334417.tif 000029
【０１０１】
このように構成された第２の実施形態の漏水量推定装置２は、第１の実施形態のセンサ配置位置選択装置１によって選択された節点に配置された水圧センサによって計測された水圧と、配水管路網の各節点を通る開管路の始点から終点までの差圧と、前記配水管路網内の漏水量の総量と、各節点への流入流量と、各節点での水の使用量と、に基づいて各節点における漏水係数を決定変数とする最適化問題を解くことで、各節点における漏水係数を決定する漏水係数最適化部と、前記決定された漏水係数と、前記配水管路網における節点間の接続関係を反映した流体の運動方程式に基づいて節点ごとの漏水量を推定する漏水量推定部とを持つことにより、より少ない水圧センサで各節点の漏水量を推定することが可能となる。
【０１０２】
図１４は、漏水量の推定のシミュレーション結果の一例を示す図である。
図１４における左図は、シミュレーション結果を示す図である。また、図１４の右図は、シミュレーションで想定された環境における実測値を示す図である。いずれの図においても、横軸は各節点を表し、縦軸は各節点における漏水量を表す。図からも明らかなように、上記方法によって各節点における漏水量が精度よく推定できることが分かる。このように、実施形態の漏水量推定装置２は、漏水量が圧力に依存して変化する場合であっても、各節点の漏水量を推定することが可能となる。
【０１０３】
以下、第２の実施形態の漏水量推定装置２の変形例について説明する。
【０１０４】
漏水量推定装置２は、このように推定された節点毎の漏水量を、例えば図１５に示す画面（画像）で、出力部２２に表示させてもよい。図１５は、実施形態の出力部２２により表示される画面の一例を示す図である。出力部２２の表示制御部は、例えば、節点毎の漏水量を棒グラフ等の態様で表示装置に表示させる。また、出力部２２の表示制御部は、漏水量の大小によって節点を示すシンボルの色を変えて表示装置に表示させてもよい。
【０１０５】
また、配水管路網にスマートメータが取り付けられていない節点が存在する場合、漏水量推定装置２は、当該節点における水の使用量を示す情報として水道メータにより取得された検針データを用いてもよい。
【０１０６】
漏水量推定装置２は、管路が埋設されてからの年数や材質、土壌などの条件に応じて、漏水に基づく管路の更新の優先順位を算出するように構成されてもよい。この場合、漏水量推定装置２は、管路が埋設されてからの年数や材質、土壌などの条件を含む埋設条件を示す埋設条件情報を記憶する埋設条件情報記憶部と、漏水量推定部１１によって推定された漏水量、及び埋設条件に基づいて管路の更新の優先順位を算出する更新優先順位算出部とを備える。
【０１０７】
以上説明した少なくともひとつの実施形態によれば、配水管路網を構成する管路の接続構成を示す情報に基づいて、配水管路網のグラフ構造を表す接続行列を生成する接続行列生成部と、接続行列に基づいて、配水管路網を構成する管路を配水管路網において開管路を構成する木となる管路と、木とともに閉管路を構成する補木となる管路とに分類し、補木となる管路を抽出する補木抽出部と、抽出された補木となる管路の両端の節点を水圧センサの配置位置として選択する選択部と、を持つことにより、水圧センサの効率的な配置位置を選択することができる。
【０１０８】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【０１０９】
１、１ａ…センサ配置位置選択装置，１１…配水管路網情報記憶部，１２…入力部，１３…出力部，１４…接続行列生成部，１５…補木抽出部，１６…配置位置選択部，１７…グラフ構造集約部，２…漏水量推定装置，２１…入力部，２２…出力部，２３…計測データ取得部，２４…管路網モデル構築部，２５…漏水係数最適化部，２６…漏水量推定部
【図１】
0006334417.tif 000030
【図２】
0006334417.tif 000031
【図３】
0006334417.tif 000032
【図４】
0006334417.tif 000033
【図５】
0006334417.tif 000034
【図６】
0006334417.tif 000035
【図７】
0006334417.tif 000036
【図８】
0006334417.tif 000037
【図９】
0006334417.tif 000038
【図１０】
0006334417.tif 000039
【図１１】
0006334417.tif 000040
【図１２】
0006334417.tif 000041
【図１３】
0006334417.tif 000042
【図１４】
0006334417.tif 000043
【図１５】
0006334417.tif 000044