(57)【特許請求の範囲】
【請求項１】
雨天時に複数の区域で発生する不明水を監視する不明水監視装置であって、
各区域の時刻毎の降雨量を蓄積する降雨量データを記憶する降雨量記憶部と、
前記複数の区域から下水処理場に流入する汚水の時刻毎の流入量を蓄積する流入量データを記憶する流入量記憶部と、
前記流入量記憶部に記憶されている晴天時の流入量と雨天時の流入量との差を降雨イベントとして定められる所定期間の雨天時の不明水の総量として求める総不明水量算出手段と、
前記降雨量記憶部に記憶されている前記降雨イベントの降雨量を加算し、各区域における前記降雨イベントにおける総降雨量を求める総降雨量算出手段と、
前記総不明水量算出手段で算出した不明水の総量と前記総降雨量算出手段で算出した各区域の総降雨量を用いて各区域の不明水量を特定する値を求め、当該値から該当する値を選択して不明水の発生する区域を推定する不明水発生区域推定手段と、
を備えることを特徴とする不明水監視装置。
【請求項２】
前記不明水発生区域推定手段によって第１期間の不明水の総量及び総降雨量とから求められた第１期間の値と、前記第１期間に続く第２期間の不明水の総量及び総降雨量とから求められた第２期間の値とを比較し、不明水の変化を推定する変化推定手段をさらに備えることを特徴とする請求項１に記載の不明水監視装置。
【請求項３】
前記不明水発生区域推定手段は、重回帰分析を利用して前記不明水の総量及び各区域の総降雨量を用いた複数の推定式を定め、当該複数の推定式から各区域の不明水量を特定する係数を求め、求めた複数の係数を比較して値の大きい係数を選択し、選択した係数に対応する区域を不明水が発生している区域と推定することを特徴とする請求項１又は２に記載の不明水監視装置。
【請求項４】
前記不明水発生区域推定手段は、前記不明水の総量及び各区域の総降雨量を用いてステップワイズ法で利用する変数を求め、求めた複数の変数から該当する変数を選択し、選択された変数に対応する区域で不明水が発生していると推定することを特徴とする請求項１又は２に記載の不明水監視装置。
【請求項５】
前記不明水発生区域推定手段で選択された値と降雨量の積を該当する区域の不明水の推定値として求める不明水発生量推定手段を備えることを特徴とする請求項１乃至４のいずれかに記載の不明水監視装置。
【請求項６】
前記降雨量記憶部に記憶されている降雨量データを読み出すとともに、所定量以上の降雨量かつ前回の雨天期間と次回の雨天期間との間隔が所定間隔以上である雨天期間を前記降雨イベントと特定する降雨イベント特定手段をさらに備えることを特徴とする請求項１乃至５のいずれかに記載の不明水監視装置。
【請求項７】
汚水管及び雨水管の工事の履歴として工事した管路の場所、工事の時期及び工事の概要を関連付ける工事データを記憶する工事履歴記憶部と、
前記不明水発生区域推定手段で不明水が発生していると推定された区域に関するデータ及び前記変化推定手段で推定された変化に関するデータに、前記工事データ記憶部に記憶される工事データを関連付けて工事による不明水改善の効果を推定する効果推定画面を生成する効果推定手段と、
前記効果推定画面を表示する表示処理手段と、
を備えることを特徴とする請求項１乃至６のいずれかに記載の不明水監視装置。
【請求項８】
前記効果推定手段は、前記効果推定画面に地図情報を組み合わせることを特徴とする請求項７に記載の不明水監視装置。
【請求項９】
雨天時に複数の区域で発生する不明水を監視する不明水監視方法であって、
各区域の時刻毎の降雨量を蓄積する降雨量データを降雨量記憶部に記憶するステップと、
前記複数の区域から下水処理場に流入する汚水の時刻毎の流入量を蓄積する流入量データを流入量記憶部に記憶するステップと、
前記流入量記憶部に記憶されている晴天時の流入量と雨天時の流入量との差を降雨イベントとして定められる所定期間の雨天時の不明水の総量として求めるステップと、
前記降雨量記憶部に記憶されている前記降雨イベントの降雨量を加算し、各区域における前記降雨イベントにおける総降雨量を求めるステップと、
算出された不明水の総量と各区域の総降雨量を用いて各区域の不明水量を特定する値を求め、当該値から該当する値を選択して不明水の発生する区域を推定するステップと、
を備えることを特徴とする不明水監視方法。
【請求項１０】
不明水の発生する区域の推定手段によって第１期間の不明水の総量及び総降雨量とから求められた第１期間の値と、前記第１期間に続く第２期間の不明水の総量及び総降雨量とから求められた第２期間の値とを比較し、不明水の変化を推定するステップをさらに備えることを特徴とする請求項９に記載の不明水監視方法。
【請求項１１】
不明水の発生する区域の推定では、重回帰分析を利用して前記不明水の総量及び各区域の総降雨量を用いた複数の推定式を定め、当該複数の推定式から各区域の不明水量を特定する係数を求め、求めた複数の係数を比較して値の大きい係数を選択し、選択した係数に対応する区域を不明水が発生している区域と推定することを特徴とする請求項９又は１０に記載の不明水監視方法。
【請求項１２】
不明水の発生する区域の推定では、前記不明水の総量及び各区域の総降雨量を用いてステップワイズ法で利用する変数を求め、求めた複数の変数から該当する変数を選択し、選択された変数に対応する区域で不明水が発生していると推定することを特徴とする請求項９又は１０に記載の不明水監視方法。
【請求項１３】
前記選択された値と降雨量の積を該当する区域の不明水の推定値として求めることを特徴とする請求項９乃至１２のいずれかに記載の不明水監視方法。
【請求項１４】
前記降雨量記憶部に記憶されている降雨量データを読み出すステップと、
所定量以上の降雨量かつ前回の雨天期間と次回の雨天期間との間隔が所定間隔以上である雨天期間を前記降雨イベントと特定するステップと
をさらに備えることを特徴とする請求項９乃至１３のいずれかに記載の不明水監視方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、生活排水や工場廃水等の下水処理に雨天時に流入する不明水の発生を監視する不明水監視装置及び不明水監視方法に関する。
【背景技術】
【０００２】
下水道には、合流式下水道と分流式下水道の２つの方式がある。合流式下水道は、排水と雨水が同一管渠を流下して、下水処理及び雨水排水処理を行なう方法である。また、分流式下水道は、汚水と雨水は別々の管渠を流下して、下水処理及び雨水排水処理を行なう方式である。近年では、分流式下水道の整備が進んでいる。
【０００３】
分流式下水道では、通常、汚水管渠には汚水のみが流入する。一方、管渠の損傷、管渠への雨水流入、管の誤接続等によって、分流式下水道には計画した汚水量よりも多くの汚水量が流入する場合がある。この通常の汚水量を超える汚水は不明水とよばれ、汚水の水質に影響を与え、下水処理における運用方法への影響や、処理コストの増大につながることから、下水処理場における大きな問題となっている。
【０００４】
分流式下水道で不明水の影響を低減するためには、不明水の発生地点を把握し、その地点を改修する必要がある。しかしながら、不明水の発生地点の把握には、管渠の埋設状況、下水処理場の運用データ又は調査用の流量計測データ等を収集・解析し、不明水の発生箇所の絞込みを行なう必要があり、これらの調査には多くの労力、時間、コストがかかる。したがって、不明水の発生箇所をできるだけ短時間かつ簡単な方法で絞り込む技術が求められている。
【０００５】
従来から、不明水の発生区域を推定する方法が様々ある（例えば、特許文献１及び非特許文献１参照）。例えば、特許文献１に記載の方法では、雨量データと不明水量データの時系列データ間のパターン解析（相関値）をもとに不明水発生分布を推定している。こｋの特許文献１に記載される推定方法では、雨量データと不明水量データの時系列データの間のパターン解析（相関値）をもとに不明水発生分布の推定を行なうものであり、各種データから不明水発生分布の推定を行なうことができる。しかしながら、推定にあたって長期間のデータや土木データが必要であり、パターン解析にあたって流達時間等を考慮した補正処理が必要になるため、程度の経験的なノウハウが必要である。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特許第３８５７６７０号公報
【非特許文献】
【０００７】
【非特許文献１】松島  修、「事例ベースモデリング技術を用いた雨天時浸入水発生領域の絞込みに関する技術マニュアルの概要」、「月間下水道」、Vol.31、No.2、p.58−64
【発明の概要】
【発明が解決しようとする課題】
【０００８】
上述したように、従来の不明水の発生区域の推定方法では、推定が困難であったり、経験的なノウハウが必要であった。
【０００９】
本発明は、上記課題に鑑み、下水処理場に流入する雨天時不明水の発生箇所や量を容易に把握することができる不明水監視装置及び不明水監視方法を提供することを目的とする。
【課題を解決するための手段】
【００１０】
上記の課題を解決するために、本発明は、雨天時に複数の区域で発生する不明水を監視する不明水監視装置であって、各区域の時刻毎の降雨量を蓄積する降雨量データを記憶する降雨量記憶部と、前記複数の区域から下水処理場に流入する汚水の時刻毎の流入量を蓄積する流入量データを記憶する流入量記憶部と、前記流入量記憶部に記憶されている晴天時の流入量と雨天時の流入量との差を降雨イベントとして定められる所定期間の雨天時の不明水の総量として求める総不明水量算出手段と、前記降雨量記憶部に記憶されている前記所定期間の降雨量を加算し、各区域における所定期間の総降雨量を求める総降雨量算出手段と、前記総不明水量算出手段で算出した不明水の総量と前記総降雨量算出手段で算出した各区域の総降雨量を用いて各区域の不明水量を特定する値を求め、当該値から該当する値を選択して不明水の発生する区域を推定する不明水発生区域推定手段とを備える。
【００１１】
また、他の形態に係る本発明は、雨天時に複数の区域で発生する不明水を監視する不明水監視方法であって、各区域の時刻毎の降雨量を蓄積する降雨量データを降雨量記憶部に記憶するステップと、前記複数の区域から下水処理場に流入する汚水の時刻毎の流入量を蓄積する流入量データを流入量記憶部に記憶するステップと、前記流入量記憶部に記憶されている晴天時の流入量と雨天時の流入量との差を降雨イベントとして定められる所定期間の雨天時の不明水の総量として求めるステップと、前記降雨量記憶部に記憶されている前記所定期間の降雨量を加算し、各区域における所定期間の総降雨量を求めるステップと、算出された不明水の総量と各区域の総降雨量を用いて各区域の不明水量を特定する値を求め、当該値から該当する値を選択して不明水の発生する区域を推定するステップとを備える。
【発明の効果】
【００１２】
本発明は、下水処理場に流入する雨天時不明水の発生箇所や量を容易に把握することができる。
【図面の簡単な説明】
【００１３】
【図１】本発明の実施形態に係る不明水監視装置を備える不明水監視システムの一例について説明する図である。
【図２】図１の不明水監視装置の一例について説明するブロック図である。
【図３】降雨量と不明水量の関係について説明する図である。
【図４】降雨イベントの決定方法について説明する図である。
【図５】不明水発生区域の推定に利用するデータの一例である。
【図６】重回帰分析を使用した不明水発生区域の推定について説明する図である。
【図７】不明水量が増加している場合の推定式の回帰係数の一例を示す図である。
【図８】不明水量が減少している場合の推定式の回帰係数の一例を示す図である。
【図９】ステップワイズ法を使用した不明水発生区域の推定について説明する図である。
【図１０】不明水監視装置で表示される画面の一例を示す図である。
【図１１】不明水監視装置で表示される工事に関する画面の一例を示す図である。
【図１２】不明水監視装置によって表示される推移に関する画面の一例を示す図である。
【図１３】図１の不明水監視装置における処理の一例を説明するフローチャートである。
【図１４】不明水監視システムの変形例について説明する図である。
【発明を実施するための形態】
【００１４】
以下に、図面を用いて本発明の実施形態に係る不明水監視装置について説明する。図１に示すように、本発明の第１の実施形態に係る不明水監視装置１０を有する不明水監視システム１は、汚水と雨水とが別々の管渠を流下する分流式下水道で使用される。この分流式下水道は、対象の地域Ａの家庭や事業所等から発生した汚水のみを汚水管２０を通して下水処理場３０に送水して処理する。下水処理場３０において適切に処理された処理水は、河川等に放流する。一方、雨水は汚水管２０とは別の管渠（図示せず）に流入して放流する。
【００１５】
この分流式下水道では、晴天時は、汚水のみが汚水管２０に流入する。一方、雨天時は、汚水とともに雨水が流入し雨天時不明水が発生することがある。この雨天時不明水の水量や発生箇所の把握に、自治体等が保有する降雨レーダ４０からの降雨量や気象情報配信事業者（図示せず）が配信する降雨量をもとにデータ処理装置５０で算出される降雨量データや、下水処理場へ流入する汚水量測定センサ６０で測定される汚水の流入量データを利用する。
【００１６】
図２に示すように、不明水監視装置１０は、中央処理装置（ＣＰＵ）１００、通信インタフェース（通信Ｉ／Ｆ）１１０、キーボードやマウス等の入力装置１２０、ディスプレイ等の出力装置１３０及び記憶装置１４０を備える情報処理装置である。この不明水監視装置１０の記憶装置１４０は、不明水の監視に利用する不明水監視プログラムＰを記憶しており、この不明水監視プログラムＰを読み出して実行することで、ＣＰＵ１００に降雨量記憶処理手段１００ａ、流入量記憶処理手段１００ｂ、降雨イベント特定手段１００ｃ、総降雨量算出手段１００ｄ、総不明水量算出手段１００ｅ、不明水発生区域推定手段１００ｆ、変化推定手段１００ｇ、不明水発生量推定手段１００ｈ、効果推定手段１００ｉ及び表示処理手段１００ｊが実装される。
【００１７】
降雨量記憶処理手段１００ａは、データ処理装置５０が送信した降雨量データを入力すると、入力した降雨量データを記憶装置１４０に記憶する。例えば、対象の地域Aは、複数の区域ａ1，ａ2，・・・ａnに区分されており、降雨量記憶処理手段１００ａは各区域に関する降雨量を降雨量データとして入力する。降雨量は時間によって変化するため、記憶装置１４０の降雨量データ１４０ａでは、異なる時間に取得された複数のデータが蓄積されている。
【００１８】
流入量記憶処理手段１００ｂは、汚水量測定センサ６０が測定した汚水の流入量データを入力すると、入力した汚水の流入量データを記憶装置１４０に記憶する。流入量も時間によって変化するため、記憶装置１４０の流入量データ１４０ｂでは、異なる時間に取得された複数のデータが蓄積されている。
【００１９】
降雨イベント特定手段１００ｃは、降雨量データを利用して対象の地域Ａの降雨イベントを特定する。
【００２０】
雨天時不明水は、図３に一例を示すように、降雨の影響により発生する不明水であるが、雨天時不明水の発生や水量は線形的な現象ではなく、土壌による浸透の影響や降雨状況による流達時間の違い等に影響される。例えば、降雨の発生後に不明水の発生する時間は一定しておらず、降雨量と雨天時不明水量の関係も異なる等、雨天時不明水は非線形性の強い現象である。しかしながら、降雨量の総量と不明水の総量との関係で捉えた場合、その関係はある程度線形的な現象として捉えることができる。そこで、不明水監視装置１０では、降雨イベントにおける区域の総降雨量と不明水の総量のデータを利用し、雨天時不明水を監視する。
【００２１】
不明水の監視の際、各降雨イベントが前回の降雨イベントの影響を受けることがないように、無降雨時間の間隔が一定時間空くようにすることが望ましい。また、このとき、降雨量が少ない場合には、誤差が生じやすく、不明水の発生区域や発生量を求めるデータとして採用することが望ましくない。したがって、降雨イベント特定手段１００ｃでは、記憶装置１４０に記憶されている降雨量データ１４０ａを読み出し、無降雨時間の間隔を利用して降雨イベントを特定する。
【００２２】
例えば、図４に示すように、降雨イベント特定手段１００ｃは、降雨量データ１４０ａとして記憶装置１４０に記憶される雨量をパルスとして表わし、前回のパルス群から次のパルス群までの間隔を無降雨時間と判断する。降雨イベント特定手段１００ｃは、無降雨時間の間隔が一定時間以上であれば、別々の降雨イベントとし、無降雨時間の間隔が一定時間未満であれば、同一の降雨イベントとする。このとき、降雨イベント特定手段１００ｃは、予め定められる所定量の雨量があった場合にのみ、処理対象の降雨イベントとして特定する。
【００２３】
総降雨量算出手段１００ｄは、降雨量データ１４０ａから対象の降雨イベントの期間における各区域ａ1〜ａnの降雨量を抽出し、区域毎に抽出した降雨量を加算して対象の降雨イベントにおける各区域ａ1〜ａnの総降雨量ｒ1〜ｒnを算出する。
【００２４】
総不明水量算出手段１００ｅは、流入量データ１４０ｂに含まれる雨天時（対象の降雨イベント）の汚水の流入量と晴天時の汚水の流入量との差を地域Ａで対象の降雨イベントにおいて発生する不明水の総量Ｑとして算出する。
【００２５】
不明水発生区域推定手段１００ｆは、総降雨量算出手段１００ｄで算出した総降雨量ｒ1〜ｒn及び総不明水量算出手段１００ｅで算出した不明水の総量Ｑを利用して、不明水監視システム１が対象とする地域Ａ内で不明水が発生している区域を推定する。例えば、不明水発生区域推定手段１００ｆは、不明水が発生している区域の推定の際に重回帰分析を利用することができる。具体的には、不明水発生区域推定手段１００ｆは、図５に示すように、複数の降雨イベントに対して総降雨量算出手段１００ｄで求められた各区域の総降雨量ｒiと総不明水量算出手段１００ｅで求められた不明水の総量Ｑを含むデータ（data1〜dataｍ）を利用して重回帰分析を行う。なお、重回帰分析を行なうにあたっては、対象とする地域Ａの区域数（ｎ）よりも多い数（ｍ＞ｎ）のデータが必要である。
【００２６】
このデータを利用して重回帰分析により得られる式は、下記の式（１）のような各区域の総雨量と不明水の総量との関係を表わす推定式である。
【００２７】
Ｑ＝α0＋α1・ｒ1＋α2・ｒ2＋α3・ｒ3＋・・・＋αn-1・ｒn-1＋αn・ｒn＋ε  ・・・（１）
Ｑ：降雨イベントにおける不明水の総量［ｍ3］
α0：切片
αi：ｉ番目の回帰係数（ｉ＝１，・・・，ｎ）
ｒi：降雨イベントにおける区域ａiの総降雨量［mm］（ｉ＝１，・・・，ｎ）
ε：誤差
不明水発生区域推定手段１００ｆは、全てのデータを利用して複数の推定式を求め、各推定式に適合するように、各区域ａ1〜ａnと対応する回帰係数α1〜αnを求める。また不明水発生区域推定手段１００ｆは、全ての回帰係数α1〜αnを求めると、求めた回帰係数を比較する。不明水発生区域推定手段１００ｆは、比較の結果、著しく大きい値の回帰係数αkと対応する区域ａkに降った雨が不明水量に大きく関係し、不明水が発生している区域であると推定する。
【００２８】
例えば、図６に示す例では、他の回帰係数と比較して、回帰係数α8が明らかに大きな値であるため、不明水発生区域推定手段１００ｆは、ａ8で不明水が発生していると推定する。
【００２９】
変化推定手段１００ｇは、降雨イベントのデータの蓄積によって、不明水発生区域推定手段１００ｆで求められた推定式を利用して、不明水の変化を推定する。具体的には、変化推定手段１００ｇは、不明水発生区域推定手段１００ｆによって第１期間の不明水の総量及び総降雨量とから求められた第１期間の推定式と、第１期間に続く第２期間の不明水の総量及び総降雨量とから求められた第２期間の推定式とを比較し、不明水の変化を推定する。
【００３０】
例えば、ある年（第１期間）の不明水の推定式（式（２））と翌年（第２期間）の不明水の推定式（式（３））があるとする。
【００３１】
Ｑb＝αb0＋αb1・ｒ1＋αb2・ｒ2＋αb3・ｒ3＋・・・＋αbn-1・ｒn-1＋αbn・ｒn＋ε  ・・・（２）
Ｑa＝αa0＋αa1・ｒ1＋αa2・ｒ2＋αa3・ｒ3＋・・・＋αan-1・ｒn-1＋αan・ｒn＋ε  ・・・（３）
このとき例えば、区域ａkに関し、ある年の回帰係数αakと比較して、翌年の回帰係数αbkが著しく増加しているとき、区域ａkにおいて不明水が発生していると推定できる。したがって、変化推定手段１００ｇは、管路の損傷の拡大、新たな管路の損傷の発生又は管路の誤接続の発生等、不明水の要因を推測することができる。
【００３２】
例えば、図７（ａ）に示すように、ある地域の２００９年の降雨イベントデータ（ｄａｔａ1〜ｄａｔａm）を利用して重回帰分析を行い回帰係数がαi1〜αi16のように求められたとする。２００９年の降雨イベントデータから求められた回帰係数αi1〜αi16によると、区域ａ８の回帰係数αi8が最も値が大きいため、区域ａ８で不明水が多く発生している可能性が高いと推測される。
【００３３】
また、図７（ｂ）は、同一の地域の２０１０年の降雨イベントデータ（ｄａｔａ1〜ｄａｔａn）を利用して重回帰分析を行った場合に求められた回帰係数αii1〜αii16の一例である。２００９年の降雨イベントデータから求めた回帰係数と２０１０年の降雨イベントデータから求めた回帰係数を比較すると、区域ａ8についての回帰係数が、１９．９６４［m3/mm］から２７．９５６［m3/mm］へと大きく増加している。区域ａ8は、２００９年にも不明水が発生している可能性が高いと推測されていた区域であるが、これにより、２０１０年にこの区域ａ8において、新たな不明水要因が発生したと推測される。したがって、この区域ａ8を優先的に管路の調査、補修又は更新等を行なう必要があると判断できる。
【００３４】
また逆に、区域ａkに関し、ある年の回帰係数αakと比較して、翌年の回帰係数αbkが著しく減少しているとき、区域ａkにおいて、管路の補修や更新等の工事が行なわれている場合、この工事が不明水の除去に効果があったと推測することができる。
【００３５】
例えば、図８（ａ）に示すように、ある地域の２００９年の降雨イベントデータ（ｄａｔａ1〜ｄａｔａmp）を利用して重回帰分析を行い回帰係数αi1〜αi16が求められたとする。２００９年の降雨イベントから求められた回帰係数αi1〜αi16によると、区域ａ8の回帰係数αi8が最も大きいため、区域ａ8で不明水が多く発生している可能性が高いと推測される。
【００３６】
また、図８（ｂ）は、同一の地域の２０１０年の降雨イベントデータ（ｄａｔａ1〜ｄａｔａq）を利用して重回帰分析を行った場合に求められた回帰係数αii1〜αii16の一例である。２００９年の降雨イベントデータから求めた回帰係数と２０１０年の降雨イベントから求めた回帰係数を比較すると、区域ａ8についての回帰係数が、２７．９５６［m3/mm］から１９．９６４［m3/mm］へと大きく減少している。これにより、区域ａ8では、２０１０年に管渠の補修等が行なわれ、不明水の要因が除去された結果、不明水が削減できたものと推測される。
【００３７】
すなわち、変化推定手段１００ｇは、不明水発生区域推定手段１００ｆによって第１期間の不明水の総量及び総降雨量とから求められた第１期間の推定式と、第１期間に続く第２期間の不明水の総量及び総降雨量とから求められた第２期間の推定式とを比較し、第１期間の推定式の係数が第２期間の推定式の係数より低い場合は、不明水が発生している可能性があると推定し、第１期間の推定式の係数が第２期間の推定式の係数より高い場合は、不明水の原因が解決した可能性が高いと推定する。
【００３８】
不明水発生量推定手段１００ｈは、不明水発生区域推定手段１００ｆで求めた値を利用して、各区域ａiにおける不明水の量を推定する。具体的には、式（１）の推定式において、回帰係数αiの単位は［m3/mm］であることから、不明水発生量推定手段１００ｈは、ある区域ａiにおいて、１mmの降雨量に対してαim3の不明水が発生する可能性があると推定することができる。したがって、不明水発生量推定手段１００ｈは、不明水発生区域推定手段１００ｆで選択された回帰係数αiと降雨量の積を発生する不明水の推定量と求める。
【００３９】
上述の説明では、不明水発生区域推定手段１００ｆは、重回帰分析を利用していたが、ステップワイズ法を利用してもよい。重回帰分析では、全ての総降雨量と不明水の総量との関係を表わす推定式（１）を用いていたが、不明水の発生に寄与しない区域もその推定式には含まれることになる。そのため、不明水発生区域推定手段１００ｆは、ステップワイズ法を利用してもよい。変数の数を取捨選択しながら、相関の高い変数のみを選択して目的とする変数を表現することが可能である。不明水発生区域推定手段１００ｆは、ステップワイズ分析で利用する各区域ａｉに対応する変数を求め、図９に示すように、求めた複数の変数から、最適な値を選択する。また、不明水発生区域推定手段１００ｆは、最適な値を選択すると、選択された変数に対応する区域ａkで不明水が発生していると推定する。
【００４０】
ステップワイズ法の場合、選択された変数で不明水発生区域を推定することが可能である。したがって、ステップワイズ法を利用した場合、不明水発生量推定手段１００ｈは、選択された変数の値を該当する区域での降雨量１mmあたりの不明水発生量と推定する。したがって、不明水発生量推定手段１００ｈは、不明水発生区域推定手段１００ｆで選択された変数と降雨量の積を発生する不明水の推定量と求める。
【００４１】
総降雨量のデータや不明水の総量のデータは、ある程度の誤差があるものと推測される。例えば、降雨量が少ない場合には、降雨量に対して発生する誤差が大きくなることが想定される。したがって、ある一定以上の降雨量があった降雨イベントを対象とすることにより、降雨量における誤差の影響を低減することができると考えられる。そのため、重回帰分析やステップワイズ分析により求められた値を精度の高い値とすることができる。
【００４２】
また、不明水発生区域推定手段１００ｆでは、重回帰分析で求めた結果とステップワイズ分析で求めた結果を蓄積し、その蓄積された結果から、最適な結果を求める手法を判断できたとき、重回帰分析又はステップワイズ分析から適切な手法を選択してその後の処理を実行するようにしてもよい。
【００４３】
上述したように、ある程度の周期（期間）で推定式を更新することにより、不明水発生の可能性の高い区域や不明水発生量の情報を経年的に捉えることで、不明水発生所の調査に有効な情報ができるものと考えられる。なお、上述した例で、推定式の更新の周期を１年としているが、対象とする区域数よりも多い降雨イベントデータを確保することができれば、必ずしも１年周期である必要はない。また、２００９年１月から２００９年１２月分、２００９年７月から２０１０年６月分、２０１０年１月から２０１０年１２月分等のように、同じ期間の降雨イベントデータを重複して利用して推定式を更新してもよい。
【００４４】
記憶装置１４０では、汚水管や雨水管の工事の履歴として、工事が行なわれた管路の場所、工事の時期及び工事の概要を少なくとも関連付ける工事データ１４０ｃを記憶している。
【００４５】
効果推定手段１００ｉは、不明水発生区域推定手段１００ｆで不明水が発生していると推定された区域に関するデータ及び変化推定手段１００ｇで推定された不明水の変化に関するデータに、記憶装置１４０で記憶されている工事データを関連付けた画面データを生成する。
【００４６】
また、表示処理手段１００ｊは、効果推定手段１００ｉで生成された画面データを出力装置１３０や、不明水監視装置１０と接続される外部の表示装置（図示せず）に表示する。例えば、表示処理手段１００ｊは、画面データから、まず、図１０に示すような推定した不明水に関する画面Ｄ１を表示する。この画面Ｄ１では、最新の不明水の推定情報に関する情報として、区域の番号と、値が大きい不明水推定量を表示する。図１０に示す例では、区域８と区域１５が、不明水が多く発生している可能性のある区域である。
【００４７】
図１０に示すように、画面Ｄ１には、工事に関する予定や履歴に関する画面Ｄ２（図１１）を表示させる工事情報表示ボタンｂ１と、不明水の推定結果の推移に関する画面Ｄ３（図１３）を表示させる不明水推定結果推移表示ボタンｂ２を有している。図１１及び図１２の例では、区域８では、２０１０年の推定結果で不明水推定量が減少しており、２０１０年に実施した工事がこの不明水の推定量の減少に効果をあげたと推測することができる。また、区域１５では、不明水推定量が徐々に増加しており、新たな不明水要因が発生していると推測することができ、この区域１５で優先的に管路の調査、補修、更新等を行なう必要があると判断することができる。
【００４８】
なお、図１０乃至図１２に示す画面Ｄ１〜Ｄ３は、数値やグラフを有する構成であるが、不明水発生区域や工事実施箇所については、地図情報と重ね合わせることも有効な標示方法である。したがって、地図情報と重ね合わせた表示画面を表示してもよい。
【００４９】
図１３に示すフローチャートを用いて、不明水監視装置１０における処理の流れについて説明する。ここでは、不明水発生区域推定手段１００ｆで重回帰分析を利用する場合の一例について説明する。
【００５０】
降雨量記憶処理手段１００ａは、各区域ａiの降雨量に関するデータを収集して記憶するとともに、流入量記憶処理手段１００ｂは、各区域ａiの汚水の流入量に関するデータを収集して記憶する（Ｓ１１）。
【００５１】
その後、降雨イベント特定手段１００ｃが記憶装置１４０の降雨量データ１４０ａを用いて降雨イベントを特定すると（Ｓ１２でＹＥＳ）、総降雨量算出手段１００ｄは、記憶装置１４０に記憶される降雨量データ１４０ａから該当する降雨イベントに関する降雨量を抽出するとともに、各区域ａiにおける総降雨量ｒiを算出し、総不明水量算出手段１００ｅは、記憶装置１４０に記憶される流入量データ１４０ｂから雨天時の流入量と晴天時の流入量を抽出して不明水の総量Ｑを算出する（Ｓ１３）。
【００５２】
必要な降雨イベントに関するデータを抽出すると（Ｓ１４でＹＥＳ）、不明水発生区域推定手段１００ｆは、推定式を生成し、各エリアにおける回帰係数を決定する（Ｓ１５）。不明水発生区域推定手段１００ｆは、全ての回帰係数を決定すると、各回帰係数を比較し（Ｓ１６）、回帰係数の高い区域を不明水発生区域として推定する（Ｓ１７）。
【００５３】
その後、不明水発生量推定手段１００ｈは、回帰係数を利用して不明水の発生量を推定する（Ｓ１８）。
【００５４】
上述したように、本発明の実施形態に係る不明水監視装置１０によれば、容易に不明水発生区域や特定の区域における不明水発生量を推定することができる。これにより、不明水の発生箇所を調査し、改修を効率的に行なうことが可能となる。
【００５５】
〈変形例〉
図１に示した不明水監視システム１では、降雨量の測定に降雨レーダ４０を利用していたが、図１４に示すように、各区域ａiに設置され、設置される各区域ａiの雨量を測定する地上雨量計７０を使用してもよい。図１４に示す例では、図１と同一の構成については、同一の符号を用いて説明を省略する。地上雨量計７０を利用する場合、不明水監視装置１０では、降雨レーダ４０からの降雨量データに代えて、地上雨量計７０を利用する他には、図１を用いて上述した場合と同様に処理される。
【符号の説明】
【００５６】
１…不明水監視システム
１０…不明水監視装置
１００…ＣＰＵ
１００ａ…降雨量記憶処理手段
１００ｂ…流入量記憶処理手段
１００ｃ…降雨イベント特定手段
１００ｄ…総降雨量算出手段
１００ｅ…総不明水量算出手段
１００ｆ…不明水発生区域推定手段
１００ｇ…変化推定手段
１００ｈ…不明水発生量推定手段
１００ｉ…効果推定手段
１００ｊ…表示処理手段
１１０…通信インタフェース
１２０…入力装置
１３０…出力装置
１４０…記憶装置（降雨量記憶部、流入量記憶部、工事履歴記憶部）
１４０ａ…降雨量データ
１４０ｂ…流入量データ
１４０ｃ…工事データ
２０…汚水管
３０…下水処理場
４０…降雨レーダ
５０…データ処理装置
６０…汚水量測定センサ
７０…地上雨量計
【図１】
0005574769.tif 000002
【図２】
0005574769.tif 000003
【図３】
0005574769.tif 000004
【図４】
0005574769.tif 000005
【図５】
0005574769.tif 000006
【図６】
0005574769.tif 000007
【図７】
0005574769.tif 000008
【図８】
0005574769.tif 000009
【図９】
0005574769.tif 000010
【図１０】
0005574769.tif 000011
【図１１】
0005574769.tif 000012
【図１２】
0005574769.tif 000013
【図１３】
0005574769.tif 000014
【図１４】
0005574769.tif 000015