JP3839361B2 - 雨水流出係数予測方法、雨水流入量予測方法、雨水流出係数予測プログラムおよび雨水流入量予測プログラム - Google Patents

Description

【０００１】
【発明の属する技術分野】
本発明は、降雨時の雨水排水を目的とする下水道ポンプ場等に設置されているポンプ施設のポンプ運転制御等のために必要となる雨水の流入量予測技術に関し、特に、降雨時に下水幹線に流入する排水区域の雨水流出係数の予測方法、予測プログラム、および下水処理場または下水ポンプ場または雨水ポンプ場への雨水流入量の予測方法、予測プログラムに関する。
【０００２】
【従来の技術】
下水処理施設は、汚水の処理だけでなく、雨水に起因する災害を防止する役割を果たしており、都市衛生の安全および環境の保持観点からも重要な施設である。
【０００３】
通常、雨は下水幹線（管渠）を経由して、ポンプ場や処理場の雨水排水ポンプから海や河川へ放出される。このため、ポンプ場や処理場の雨水排水ポンプは、ポンプ場や処理場への雨水流入量に応じて運転台数が決定されることが重要である。
【０００４】
近年、住宅地の密集化や舗装道路の普及に伴い、雨水が下水管へ流入する割合（雨水流出係数）が多くなってきている。さらには、ポンプの排水処理能力に限界があるため、貯留管や調節池とよばれる一時的な雨水貯留施設を設置することが多くなってきている。このような雨水貯留施設を利用する際には、雨水貯留施設へのゲートの開度を適切に制御することが重要である。
【０００５】
従来、降雨時に、排水区域から下水幹線に流入し、ポンプ場等へ流入する雨水量を演算するモデルには、下記のようなものがあった。
（１）地上雨量計やレーダー雨量計で計測された降雨量をもとに、ＲＲＬ（Ｒｏａｄ Ｒｅｓｅａｒｃｈ Ｌａｂｏｒａｔｏｒｙ）法または修正ＲＲＬ法で懸案地点の流出量を求める。
（２）地表面の雨水の挙動や管渠内の挙動を数式で厳密に表わし、与えられた境界条件、初期条件のもとで懸案地点の流出量を求める。
（３）降雨量や下水管渠内水位や流量計測値の時系列データと現在までの計測された流出量時系列データを入力として将来の流出量を出力とするニューラルネットワークにより、懸案地点の流出量を求める。
【０００６】
ＲＲＬ法とは、英国道路研究所で開発された雨水流入量を算定する方法で、概説すると以下のような方法である。
まず、対象領域の管渠の長さ、直径、勾配等の水力学的特徴を記入した管渠図を作成する。管渠図全体を一つの流域（単一流域）とみて、流域の最下流点を流量算定点として選択する。そして、開水路等における流量速度を算出し、地点ｐまでの雨水流達時間が流量計算時間間隔と等しくなるような等到達時間曲線を作成する。さらに、等到達時間曲線で区分される面積を時間域別面積Ａ_i［ｍ²］として算出し、時間面積図を作成する。
【０００７】
次に、流量計算時刻ｉでの流域に降る降雨強度Ｉ_i［ｍｍ／ｓ］から降雨量曲線を作成する。ここで、流出の割合を表す流出係数Ｃは、領域の土地の利用状態から求められる固定値で与える。作成した降雨量曲線と時間面積図から雨水流入量Ｐ_i［ｍ³／ｓ］を仮に下記の式１で算出する。
【０００８】
【数１】

【０００９】
一般に、降った雨はすぐには流出せず、一度管渠に貯えられてから流出するため、仮に求めた雨水流入量Ｐ_iを雨水貯留量と流出量との関係式（式２、式３）に当てはめて流出量Ｑを算定する。
Ｓ（ｔ）＝ＫＱ（ｔ）ⁿ ・・（式２）
Ｋ，ｎ：流域により定まる定数
ｄＳ（ｔ）／ｄｔ＝Ｐ（ｔ）−Ｑ（ｔ） ・・（式３）
Ｓ：貯留量［ｍ³］，Ｐ：仮想流入量［ｍ³／ｓ］，Ｑ：流出量［ｍ³／ｓ］
上記式２は、いわゆる運動方程式であり、式３は、いわゆる連続方程式である。通常は式３を差分化した下記の式４を利用して流出量の算出を行なう。
Ｑ_i＋（２／Δｔ）Ｓ_i＝（Ｐ_i＋Ｐ_i-1）−Ｑ_i-1＋（２／Δｔ）Ｓ_i-1 ・・（式４）
Δｔ：刻み時間［ｓ］
このようにして算出した流出量を流出算定地点ｐの雨水流入量Ｑ_iとして最終的に決定する。
【００１０】
以上のようなＲＲＬ法に、管内貯留量の計算方法ならびに有効降雨の算定方法に改良を加えた修正ＲＲＬ法という方法もある。
これらの中で、上記（１）のモデルでは降雨量に対して、どれだけの量が流出するかを表わすパラメータとして流出係数がある。この値は時間的変動を考慮せず、設計時に良く用いられる合理式の場合も含め、一義的（一定値）に与えられる。
【００１１】
一方、降雨現象の時間的変動特性を考慮した事例として「流出係数推定方法及び流出係数予測システム：特開平５−２６３７６７号公報」がある。この手法の特徴は、流出係数を一定値として捉えず、降雨開始直後から時々刻々と収束値（総括流出係数）に向けて変動する時間の関数として捉えている点にある。
【００１２】
【発明が解決しようとする課題】
従来のモデルに用いられている流出係数は使用する者にとって分かりやすいパラメータである。その値は総括流出係数も含め観測値や学術的な推定により決定される。しかも、実績値による調整要素がなく、降雨や地表面の状態による流出特性の変動を考慮していない。
【００１３】
そのため、予測流入量と実績値との予測誤差を大きくする原因となっていたという問題点があった。
本発明は、降水量のうち下水道に流入する水量の比率を示す雨水流出係数、および施設に流入する雨水流入量を、実測値との予測誤差が極力小さくなるように予測することが可能な雨水流出係数予測方法、雨水流入量予測方法、雨水流出係数予測プログラムおよび雨水流入量予測プログラムを提供することを目的とする。
【００１４】
【課題を解決するための手段】
本発明は、上記課題を解決するため、下記のような構成を採用した。
下水処理場または下水ポンプ場または雨水ポンプ場において、降雨時の流入量を予測するために、降雨計と上流幹線内に複数の水位計または流量計を設置し、降雨情報とそれから得られる降雨強度と幹線内の水位情報または該地点の直接計測する管渠内流量により得られる排水区域から幹線への雨水流入量を用いて、該排水区域の流出係数を演算し、上記時系列情報を用いて該排水区域の雨水流出係数を予測する手段を持つことを特徴とする。
【００１５】
また、上記で得られた雨水流出係数予測値を用いて、該排水区域から管渠へ流入する雨水流入量を予測し、その結果に基づいて対象施設に流入する雨水流入量を予測する手段を持つことを特徴とする。
【００１６】
すなわち、本発明の一態様によれば、本発明の雨水流出係数予測方法は、降水量のうち下水道に流入する水量の比率を示す雨水流出係数を予測するシステムにおいて実行される雨水流出係数予測方法であって、まず、降雨量を計測し、上記計測した降雨量に基づいて降雨強度を算出し、上記計測した降雨量および上記算出した現時点までの降雨強度を時系列に格納し、次に、下水幹線内の水位である管渠内水位を計測し、上記計測した管渠内水位に基づいて上記下水幹線内への雨水流入量を算出し、上記計測した管渠内水位および上記算出した現時点までの雨水流入量を時系列に格納し、そして、上記算出した降雨強度および上記算出した雨水流入量に基づいて、雨水流出係数を算出し、上記算出した現時点までの雨水流出係数を時系列に格納し、さらに、上記格納された現時点までの雨水流出係数と降雨強度と雨水流入量とに基づいて、将来の雨水流出係数を予測することを特徴とする。
【００１７】
また、本発明の雨水流出係数予測方法は、上記将来の雨水流出係数の予測には、カルマンフィルターを適用することが望ましい。
また、本発明の一態様によれば、本発明の雨水流入量予測方法は、対象施設に流入する雨水流入量を予測するシステムにおいて実行される雨水流入量予測方法であって、まず、上述の雨水流出係数予測方法において算出した降雨強度計測値と、上述の雨水流出係数予測方法において予測した雨水流出係数予測値とに基づいて、下水幹線への雨水流入量を予測し、次に、上記予測した下水幹線への雨水流入量に基づいて、上記対象施設に流入する雨水流入量を予測することを特徴とする。
【００１８】
また、本発明の雨水流入量予測方法は、さらに、上記下水幹線内の流量を計測し、そして、上記対象施設に流入する雨水流入量の予測において、上記予測した下水幹線への雨水流入量に加え、上記計測した下水幹線内の流量に基づいて行われることが望ましい。
【００１９】
また、本発明の雨水流出係数予測プログラムは、上述の雨水流出係数予測方法を実行させるための雨水流出係数予測プログラムである。
また、本発明の雨水流入量予測プログラムは、上述の雨水流入量予測方法を実行させるための雨水流入量予測プログラムである。
【００２０】
【発明の実施の形態】
以下、本発明の実施の形態を、図１乃至図３を参照しながら詳細に説明する。
図１は、本発明を適用した下水道プラントの概念図を示す図である。
【００２１】
図１において、排水区域１０１と排水区域１０２とから集められた雨水は、それぞれマンホール１０３およびマンホール１０４の地点で下水幹線（管渠）１０５に流入される。そして、下水幹線１０５へ流下された雨水は、その後、ポンプ場１０６からポンプ１０７で排水されることを表わしている。
【００２２】
図２は、本発明の第１の実施の形態に係る雨水流出係数予測システムの機能ブロック図である。
本発明に係る雨水流出係数予測方法について、図１および図２を用いて説明する。
【００２３】
排水区域１０１に降った雨の降雨量は、排水区域１０１周辺に設置された降雨量計測手段（雨量計１０８）２０１で計測された後、マンホール１０３に集水され、下水幹線１０５に流出する。
【００２４】
降雨量計測手段２０１で計測された排水区域１０１の降雨量は、降雨強度演算手段２０２によって降雨強度（降水量）に変換された後、降雨量計測手段２０１によって計測された降雨量とともに時系列データとして、降雨量、降雨強度記憶手段２０３に記憶される。
【００２５】
他方、マンホール１０３の下流に設置された管渠内水位計測手段（水位計１０９）２０４によって下水幹線内の水位を計測し、管渠への雨水流入量演算手段２０５によって下記のマニング流速公式を用いて下水幹線内の流量を算出する。
Ｖ＝（１／ｎ）Ｒ^2/3Ｉ^1/2 ・・（式５）
Ｑ₁＝Ｖ・Ａ ・・（式６）
Ｖ：流速、ｎ：粗度係数（下水管表面を水が流れるときの抵抗）、Ｒ：径深（水路形状および水位により変化する値）、Ｉ：管渠勾配、Ｑ₁：管渠内流量、Ａ：通水断面積
雨水流入量は、管渠内水位計測手段（水位計１０９）２０４の地点で計測された流量がマンホール１０３からの流入量そのものであるから、上記式６のＱ₁と同じである（この例では、マンホール１０３以外からの雨水の流れはないと想定）。そして、得られた管渠内流量および雨水流入量を水位とともに時系列データとして、水位、雨水流入量、管渠内流量記憶手段２０６に記憶する。
【００２６】
そして、雨水流出係数演算手段２０７は、降雨強度演算手段２０２によって求めた降雨強度と、管渠への雨水流入量演算手段２０５によって求めたマンホール１０３からの雨水流入量（式６のＱ₁に等しい）と、排水区域面積とを用いて、下記の式７により雨水流出係数（Ｃ）を求め、流出係数記憶手段２０８に記憶する。この場合、排水区域１０１からのマンホール１０３までの流達時間はほぼ一定（ｔｄ）とする。
Ｃ（ｔ）＝３６０×Ｑｍ₁（ｔ）／（Ｒ（ｔ−ｔｄ）・Ａ） ・・（式７）Ｃ（ｔ）：時刻ｔの雨水流入係数
Ｑｍ₁（ｔ）：時刻ｔの雨水流入量［ｍ³／ｓ］
ｔｄ：流達時間
Ｒ（t−ｔｄ）：時刻ｔからｔｄ前の降雨強度［ｍｍ／ｈ］
ただし、降雨開始からの降雨量［ｍｍ］が、部分排水区の地表流出が生じる降雨量限界値［ｍｍ］を越えるまでは、Ｒ（t−ｔｄ）の値を０（ゼロ）とする。また、演算間隔内で越えた場合、残りの降雨が流出に寄与するものとする。
Ａ：排水区域１０１の面積［ｈａ］
次に、排水区域１０２に降る雨について考える。
【００２７】
排水区域１０２に降った雨の降雨量は、排水区域１０２周辺に設置された降雨量計測手段（雨量計１１０）２０１で計測された後、マンホール１０４に集水され、下水幹線１０５に流出する。
【００２８】
降雨量計測手段２０１で計測された排水区域１０２の降雨量も、排水区域１０２の降雨量と同様に、降雨強度演算手段２０２によって降雨強度に変換された後、降雨量計測手段２０１によって計測された降雨量とともに時系列データとして、降雨量、降雨強度記憶手段２０３に記憶される。
【００２９】
また、マンホール１０４の下流に設置された管渠内水位計測手段（水位計１１１）２０４によって下水幹線内の水位を計測し、管渠への雨水流入量演算手段２０５によって上述のマニング流速公式（式５、式６）を用いて管渠内流量Ｑ₂（ｔ）を求める。
【００３０】
そして、下記の式８から上記流量Ｑ₁を用いて、マンホール１０４から流入する雨水流入量Ｑｍ₂（ｔ）を求める。
Ｑｍ₂（ｔ）＝Ｑ₂（ｔ）−Ｑ₁（ｔ−Δｔ） ・・（式８）
Δｔ：マンホール１０３とマンホール１０４の間の流達時間
次に、各部分排水区１１２乃至１１５を流達時間（例えば５分間隔）で区分し、その面積と降雨強度時系列データとを用いて、平均雨水流出係数（Ｃ）を下記の式９により求める。
Ｃ（ｔ）＝３６０×Ｑｍ₂（ｔ）／Σ（Ｒ_k（ｔ−ｔｄ_k）・Ａ_k）・・（式９）
Ｑｍ₂（ｔ）：時刻ｔの雨水流入量［ｍ³／ｓ］
Ｒ_k（ｔ−ｔｄ_k）：部分排水区ｋの時刻ｔからｔｄ_k前の降雨強度［ｍｍ／ｈ］ただし、降雨開始からの降雨量［ｍｍ］が、部分排水区ｋの地表流出が生じる降雨量限界値［ｍｍ］を越えるまでは、Ｒ_k（t−ｔｄ_k）の値を０（ゼロ）とする。また、演算間隔内で越えた場合、残りの降雨が流出に寄与するものとする。
Ａ_k：部分排水区ｋの面積［ｈａ］
ここで、Σはｋについての加算、すなわちすべての部分排水区に対して加算することを示す。図１に示した排水区域１０２では、部分排水区１１２乃至部分排水区１１５までに対する加算を示す。
【００３１】
このような演算を繰返すことにより、多数マンホールから雨水が流入しても同様な方法で雨水流出係数を求めることができる。このとき、▲１▼すべての流入マンホールの流出係数を必要としない場合、必要とするマンホールの上流側に水位計を設置し、管渠内流量を算出し、式８のＱ₁とすること、▲２▼管渠内水位計の代わりに直接流量を計測して演算に用いること、▲３▼地上雨量計の代わりに、レーダー雨量計の計測値を用いること、を行っても良い。
【００３２】
上述の式７または式９で得られた雨水流出係数には計測誤差や排水区域の降雨の非一様性あるいは地表貯留・浸透・流下変動に起因するノイズを含んでいる。
そこで、雨水流出係数予測手段２０９は、降雨量、降雨強度記憶手段２０３、水位、雨水流入量管渠内流量記憶手段２０６、流出係数記憶手段２０８に記憶された、現在までの雨水流出係数実績値や降雨強度、雨水流入量の時系列情報を入力とするカルマンフィルター手法を用いてノイズを除去し、将来の雨水流出係数を予測する。このとき、実績降雨強度に加えて、降雨強度予測値を用いて予測精度を向上させることもできる。
【００３３】
なお、カルマンフィルターの代わりにニューラルネットワークや回帰式を用いても良い。
図３は、本発明の第２の実施の形態に係る雨水流入量予測システムの機能ブロック図である。
【００３４】
本発明に係る雨水流入量予測方法について、図１および図３を用いて説明する。
管渠への雨水流入量予測手段３０３は、上述の第１の実施の形態で予測した雨水流出係数予測値３０１と、上述の第１の実施の形態で求めた降雨強度値３０２とを用いて、各排水区域（排水区域１０１および１０２）から下水幹線１０５への雨水流入量を下記の式（式１０、式１１）により予測する。
マンホール１の場合
Ｑｍ₁（ｔ）＝ Ｃ（ｔ）・Ｒ（ｔ−ｔｄ）・Ａ／３６０ ・・（式１０）
マンホール２の場合
Ｑｍ₂（ｔ）＝ Ｃ（ｔ）・Σ（Ｒ_k（t−ｔｄ_k）・Ａ_k）／３６０・・（式１１）
だだし、Σはｋについての加算
なお、降雨強度として計測値（降雨強度演算手段２０２による演算値）だけでなく予測値も用いることにより、将来の雨水流入量の予測ができる。また、地上雨量だけでなく、レーダー雨量を用いてもよい。
【００３５】
次に、ポンプ場（または処理場または雨水ポンプ場）１０６への雨水流入量の予測を行なう。予測演算には、すべての流入マンホール（マンホール１０３、１０４）に対して、ポンプ場１０６までの流達時間を考慮した下記の式１２を用いてポンプ場流入量予測手段３０４が実行する。
【００３６】
このとき、現在時刻までの雨水流入量は、管渠内流量計測手段３０５で計測され、管渠内流量記憶手段３０６に記憶された値を用い、現在時刻よりの予測値は式１０、式１１の値を用いる。
Ｑｐ（ｔ）＝ΣＱｍ_i（ｔ−ｔｄ_i） ・・（式１２）
Ｑｐ（ｔ）：時刻ｔにおけるポンプ場流入量［ｍ³／ｓ］
Ｑｍ_i（ｔ−ｔｄ_i）：流入マンホールｉの時刻ｔからｔｄ_i前の雨水流入量［ｍ³／ｓ］
ｔｄ_i：流入マンホールｉからポンプ場までの流達時間
Σ：ｉに対する加算
なお、雨水流入量の予測には、以下の方法を用いてもよい。
（１）上述の雨水流入量予測情報に加えて、実測した管渠内水位や流量、ポンプ場水位やポンプ運転状態をも入力とし、管渠内の流れを不定流運動方程式と連続の式を用いて水理解析により、ポンプ場への流入量を予測する。
（２）上述の雨水流入量予測情報に加えて、実測した管渠内水位や流量、ポンプ場水位の時系列データをも入力とするニーラルネットワークにより、ポンプ場への流入量を予測する。
【００３７】
これらの場合にも、管渠内流量計測手段３０５によって計測され管渠内流量記憶手段３０６に時系列に記憶されている管渠内流量値を用いる。
以上、本発明の実施の形態を、図面を参照しながら説明してきたが、本発明が適用される雨水流出係数予測システムおよび雨水流入量予測システムを制御するプログラムは、そのプログラムコードを記録したＲＯＭやＲＡＭのメモリ、外部記録装置、可搬記録媒体を、雨水流出係数予測システムおよび雨水流入量予測システムに供給し、その雨水流出係数予測システムおよび雨水流入量予測システムのコンピュータがプログラムコードを読み出し実行することによっても、達成されることは言うまでもない。
【００３８】
この場合、記録媒体から読み出されたプログラムコード自体が本発明の新規な機能を実現することになり、そのプログラムコードを記録した可搬記録媒体等は本発明を構成することになる。
【００３９】
プログラムコードを供給するための可搬記録媒体としては、例えば、フレキシブルディスク、ハードディスク、光ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＣＤ−Ｒ、ＤＶＤ−ＲＯＭ、ＤＶＤ−ＲＡＭ、磁気テープ、不揮発性のメモリーカード、ＲＯＭカード、電子メールやパソコン通信等のネットワーク接続装置（言い換えれば、通信回線）を介して記録した種々の記録媒体などを用いることができる。
【００４０】
また、コンピュータがメモリ上に読み出したプログラムコードを実行することによって、前述した実施の形態の機能が実現される他、そのプログラムコードの指示に基づき、コンピュータ上で稼動しているＯＳなどが実際の処理の一部または全部を行ない、その処理によっても前述した実施の形態の機能が実現される。
【００４１】
さらに、可搬型記録媒体から読み出されたプログラムコードが、コンピュータに挿入された機能拡張ボードやコンピュータに接続された機能拡張ユニットに備わるメモリに書き込まれた後、そのプログラムコードの指示に基づき、その機能拡張ボードや機能拡張ユニットに備わるＣＰＵなどが実際の処理の一部または全部を行ない、その処理によっても前述した実施の形態の機能が実現され得る。
【００４２】
すなわち、本発明は、以上に述べた実施の形態に限定されるものではなく、本発明の要旨を逸脱しない範囲内で種々の構成または形状を取ることが出来る。
【００４３】
【発明の効果】
以上説明してきたように、本発明によれば、排水区域の雨水流出係数を一定値に設定するのはなく、降雨強度と排水区域からの雨水流入量から、時々刻々演算するので、正確な雨水流出係数の予測と雨水流入量予測を行なうことができ、ポンプ場への流入量の予測精度を向上させることができる。これによりポンプの最適運転が可能となり、さらに、ポンプ省エネルギー運転、浸水回避、雨天汚濁負荷流出削減、運転員の負荷軽減が可能となる。
【図面の簡単な説明】
【図１】本発明を適用した下水道プラントの概念図を示す図である。
【図２】本発明の第１の実施の形態に係る雨水流出係数予測システムの機能ブロック図である。
【図３】本発明の第２の実施の形態に係る雨水流入量予測システムの機能ブロック図である。
【符号の説明】
１０１、１０２ 排水区域
１０３、１０４ マンホール
１０５ 下水幹線（管渠）
１０６ ポンプ場
１０７ ポンプ
１０８ 雨量計
１０９ 水位計
１１０ 雨量計
１１１ 水位計
１１２、１１３、１１４、１１５ 部分排水区
２０１ 降雨量計測手段
２０２ 降雨強度演算手段
２０３ 降雨量、降雨強度記憶手段
２０４ 管渠内水位計測手段
２０５ 管渠への雨水流入量演算手段
２０６ 水位、雨水流入量管渠内流量記憶手段
２０７ 雨水流出係数演算手段
２０８ 流出係数記憶手段
２０９ 雨水流出係数予測手段
３０１ 雨水流出係数予測値
３０２ 降雨強度計測値
３０３ 管渠への雨水流入量予測手段
３０４ ポンプ場流入量予測手段
３０５ 管渠内流量計測手段
３０６ 管渠内流量記憶手段

Claims (6)

1. 降水量のうち下水道に流入する水量の比率を示す雨水流出係数を予測するシステムにおいて実行される雨水流出係数予測方法であって、
   降雨量を計測し、前記計測した降雨量に基づいて降雨強度を算出し、前記計測した降雨量および前記算出した現時点までの降雨強度を時系列に格納し、
   下水幹線内の水位である管渠内水位を計測し、前記計測した管渠内水位に基づいて前記下水幹線内への雨水流入量を算出し、前記計測した管渠内水位および前記算出した現時点までの雨水流入量を時系列に格納し、
   前記算出した降雨強度および前記算出した雨水流入量に基づいて、雨水流出係数を算出し、前記算出した現時点までの雨水流出係数を時系列に格納し、
   前記格納された現時点までの雨水流出係数と降雨強度と雨水流入量とに基づいて、将来の雨水流出係数を予測することを特徴とする雨水流出係数予測方法。
2. 前記将来の雨水流出係数の予測に、カルマンフィルターを適用することを特徴とする請求項１に記載の雨水流出係数予測方法。
3. 対象施設に流入する雨水流入量を予測するシステムにおいて実行される雨水流入量予測方法であって、
   請求項１または２に記載の雨水流出係数予測方法において算出した降雨強度計測値と、請求項１または２に記載の雨水流出係数予測方法において予測した雨水流出係数予測値とに基づいて、下水幹線への雨水流入量を予測し、
   前記予測した下水幹線への雨水流入量に基づいて、前記対象施設に流入する雨水流入量を予測することを特徴とする雨水流入量予測方法。
4. さらに、前記下水幹線内の流量を計測し、
   前記対象施設に流入する雨水流入量の予測は、前記予測した下水幹線への雨水流入量に加え、前記計測した下水幹線内の流量に基づいて行われることを特徴とする請求項３に記載の雨水流入量予測方法。
5. コンピュータに請求項１または２に記載の雨水流出係数予測方法を実行させるための雨水流出係数予測プログラム。
6. コンピュータに請求項３または４に記載の雨水流入量予測方法を実行させるための雨水流入量予測プログラム。

Images