所定状態の液体流の気泡から発せられる音響波形における、周波数に対する音圧のプロファイルが極大となる周波数を基準周波数に設定する基準周波数設定ステップと、任意状態の液体流の気泡から発せられる音響波形を測定する音響波形測定ステップと、前記測定された音響波形における周波数に対する音圧のプロファイルが極大となる周波数と前記基準周波数との周波数差を算出する周波数差算出ステップと、前記周波数差に基づいて、前記任意状態の液体流の水位を算出する水位算出ステップと、を含むことを特徴とする水位測定方法。前記液体流が、下水管渠内または雨水管渠内を流れる水流であることを特徴とする請求項1に記載の水位測定方法。所定状態の液体流の気泡から発せられる音響波形における、周波数に対する音圧のプロファイルにおいて音圧が極大となる周波数を基準周波数に設定する基準周波数設定手段と、任意状態の液体流の気泡から発せられる音響波形を測定する音響波形測定手段と、前記音響波形測定手段によって測定された音響波形における周波数に対する音圧のプロファイルにおいて音圧が極大となる所定周波数と、前記基準周波数との周波数差を算出する周波数差算出手段と、前記周波数差算出手段によって算出された前記周波数差に基づいて、前記任意状態の液体流の水位を算出する流速算出手段と、を含むことを特徴とする水位測定システム。前記液体流との接触によって気泡を発生させる気泡発生手段をさらに備えることを特徴とする請求項3に記載の水位測定システム。前記液体流との接触によって気泡を発生させるとともに衝突音を発する気泡発生手段をさらに備えることを特徴とする請求項3に記載の水位測定システム。前記気泡発生手段が、下水管渠または雨水管渠の内壁に、水位の増減方向に沿って複数設置されることを特徴とする請求項3~5のいずれか1項に記載の水位測定システム。前記液体流が、下水管渠内または雨水管渠内を流れる水流であることを特徴とする請求項3~6のいずれか1項に記載の水位測定システム。しかしながら、菅渠内の雨水などの流水といった液体流における物理量を測定する水位計などの各種測定計器には、測定値などの各種情報データを管理室の端末に送信するための通信手段や、それらを駆動するための電源が必要になる。さらに、これらの測定計器は、多数存在するマンホールごとに設置されているため、測定手段としての測定計器におけるメンテナンスの手間やメンテナンスコスト、さらに設置コストが高くなっていた。本発明は、上記に鑑みてなされたものであって、その目的は、液体流の水位を測定するための測定手段の設置数を低減することができ、測定手段の設置やメンテナンスに要するコストの低減を図ることができる水位測定方法および水位測定システムを提供することにある。上述した課題を解決し、上記目的を達成するために、本発明者は鋭意検討を行った。以下にその概要を説明する。まず、本発明者は、上述した従来技術における流速計や水位計の設置数を低減させる方法について検討を行った。そして、本発明者は、従来の流速計や水位計のように水流に直接接触させて測定する測定計器を用いるよりも、非接触で流速や水位を計測することができれば、測定計器の設置数を低減できることを想起した。そこで、本発明者が水流の流速や水位の変動に応じて変化する物理量について重ねて検討を行ったところ、水流から発せられる音響波形に基づいて水流の流速や水位を測定することができれば、菅渠内の長い距離において、非接触で流速や水位を測定できることを想起するに至った。すなわち、本発明者の知見によれば、液体の流れ液体流の音は、液体流中に生じる気泡により発生する。また、気泡の発生や破裂による音響波形に対して、フーリエ変換解析FT解析、具体的には例えば高速フーリエ変換解析FFT解析を行うと、数100Hz近傍の低周波域に音圧のピークが存在することも知られている非特許文献1参照。本発明者は、これらの知見に基づいて種々実験を行い、液体中の気泡群を検知する場合に、液体流中に気泡が発生する条件において、液体流が増加して水位が上昇すると液体流中の気泡群の密度が低下し、反対に液体流が減少して水位が低下すると気泡群の密度が上昇するという相関を見出した。さらに、本発明者は、実験およびこれに基づいた検討を行って、図1に示すような水位とピーク音圧における周波数との相関を見出した。すなわち、本発明者は、液体流中の気泡群の密度が低下すると、音圧Lpが極大となる周波数が、所定周波数f0から所定周波数f1に低下すること、反対に、液体流中の気泡群の密度が上昇すると、音圧Lpが極大となる周波数が、所定周波数f1から所定周波数f0に上昇することを見出した。これらのことから、本発明者は、液体流の水位が所定の水位から増加した場合には、液体流中の気泡の密度が低下して音圧が極大になる周波数が低下し、反対に液体流の水位が所定の水位から低下した場合には、液体流中の気泡の密度が上昇して音圧が極大となる周波数が上昇するという物理現象を、新たに見出すに至った。このような液体流と音圧が極大になる周波数との相関に関する物理現象に基づいて、本発明者は、所定状態の液体流において音圧が極大となる周波数を基準周波数としてあらかじめ測定し設定した後、菅渠内などを流れる液体流の通常の状態である任意状態において、液体流から発生する音響波形を測定して、この音響波形から周波数に対する音圧を解析し、音圧が極大となる周波数と基準周波数との差に応じた液体流の水位の相関に基づいて、液体流の水位を算出できることを想起するに至った。本発明は、以上の鋭意検討に基づいて案出されたものである。したがって、本発明に係る水位測定方法は、所定状態の液体流の気泡から発せられる音響波形における、周波数に対する音圧のプロファイルが極大となる周波数を基準周波数に設定する基準周波数設定ステップと、任意状態の液体流の気泡から発せられる音響波形を測定する音響波形測定ステップと、測定された音響波形における周波数に対する音圧のプロファイルが極大となる周波数と基準周波数との周波数差を算出する周波数差算出ステップと、周波数差に基づいて、任意状態の液体流の水位を算出する水位算出ステップと、を含むことを特徴とする。本発明に係る水位測定方法は、上記の発明において、液体流が、下水管渠内または雨水管渠内を流れる水流であることを特徴とする。本発明に係る水位測定システムは、所定状態の液体流の気泡から発せられる音響波形における、周波数に対する音圧のプロファイルにおいて音圧が極大となる周波数を基準周波数に設定する基準周波数設定手段と、任意状態の液体流の気泡から発せられる音響波形を測定する音響波形測定手段と、音響波形測定手段によって測定された音響波形における周波数に対する音圧のプロファイルにおいて音圧が極大となる所定周波数と、基準周波数との周波数差を算出する周波数差算出手段と、周波数差算出手段によって算出された周波数差に基づいて、任意状態の液体流の水位を算出する流速算出手段と、を含むことを特徴とする。本発明に係る水位測定システムは、上記の発明において、液体流との接触によって気泡を発生させる気泡発生手段をさらに備えることを特徴とする。また、本発明に係る水位測定システムは、上記の発明において、液体流との接触によって気泡を発生させるとともに衝突音を発する気泡発生手段をさらに備えることを特徴とする。この構成において、気泡発生手段は、下水管または雨水管などの菅渠の内壁に、水位の増減方向に沿って複数設置されることを特徴とする。本発明に係る水位測定システムは、上記の発明において、液体流が、下水管渠内または雨水管渠内を流れる水流であることを特徴とする。本発明に係る水位測定方法および水位測定システムによれば、測定された音波に基づいて液体流の水位を算出していることにより、従来の接触型の測定手段である測定計器の場合に比して、非接触型の測定手段の設置場所および設置数を低減することができるので、測定手段の設置やメンテナンスの低コスト化を図ることが可能となる。以下、本発明の一実施形態について図面を参照しつつ説明する。なお、以下の一実施形態の全図においては、同一または対応する部分には同一の符号を付す。また、本発明は以下に説明する一実施形態によって限定されるものではない。水位測定システムまず、本発明の一実施形態による水位測定システムについて説明する。図2は、この一実施形態による水位測定システムの全体構成を示す。図2に示すように、この一実施形態による水位測定システムにおいては、処理場管理室2内の情報処理装置21と、マンホール3内の耐水耐衝撃用の保護ケース31に格納された音響信号変換手段としての音響信号変換装置32の伝送装置32aとが、ネットワーク網1を通じて通信可能に接続されている。ネットワーク網1は、インターネットやイントラネット、または携帯端末通信網などの通信網である。また、音響信号変換装置32は、菅渠4内に設置された音響波形測定手段としてのマイクロホン41に接続されている。さらに、菅渠4内には、必要に応じて、気泡を発生させるための気泡発生装置42が設置される。処理場管理室2に備えられた情報処理装置21は、情報処理部21a、情報記録部21b、およびインターフェース21cを備える。情報処理部21aは、主として、情報データの入出力によって種々の演算を行う中央演算処理装置CPUとROMやRAMなどの記憶部とから構成されている。インターフェース21cは、情報処理部21aをネットワーク網1に接続して通信するための通信中継手段である。そして、マンホール3に備えられた伝送装置32aから送信される音響波形のデータや音響信号などの各種情報データは、ネットワーク網1を通じて情報処理装置21に供給され、インターフェース21cを介して情報処理部21aおよび情報記録部21bに供給される。情報記録部21bには、各種計算処理ソフトウェア、音響解析ソフトウェア、水位計算ソフトウェア、および流速計算ソフトウェアなどのプログラムがインストールされている。そして、情報記録部21bと情報処理部21aとの間において情報データの入出力を行うとともに、情報処理部21aがそれぞれのプログラムに従った処理を実行する。この情報処理部21aは、実行されるプログラムに応じて、種々の処理を行う処理手段である。また、情報記録部21bには、各種データベースがデータテーブルとして格納されている。この各種データベースとしては、菅渠4における管の勾配、材質、マニング粗度係数、および口径などの、いわゆる管の寸法に関する種々の寸法データベース、基準周波数からの周波数差に基づいた水位の加減値の水位計算用データベース、後述する流水43の水位を算出するための周波数差や、周波数と水流の水位との相関データベース、または気泡発生装置42ごとの固有の周波数のデータベースなどを挙げることができる。また、菅渠4内のマイクロホン41は、例えば指向性が高いパラボラ集音機などから構成される。このマイクロホン41においては、集音した音に基づいて電気信号が生成される。マイクロホン41により生成された電気信号は、音響信号変換装置32に供給される。音響信号変換装置32は、マイクロホン41から供給された電気信号を音響信号に変換し、この音響信号を、伝送装置32aを介しネットワーク網1を通じて処理場管理室2の情報処理装置21に供給する。菅渠4内に設置される気泡発生装置42は、液体流としての流水43に接触することで気泡を発生する気泡発生手段であり、菅渠4内において気泡が発生しにくい部分などに設置される。なお、気泡発生装置42の詳細については後述する。水位測定方法次に、以上のように構成された菅渠内の水流の水位測定システムによる水位測定方法について説明する。図3は、この一実施形態による水位測定方法を示すフローチャートである。図3に示すように、この一実施形態による水位測定方法においては、まず、ステップST1において、流速の計測対象となる流水43が存在する菅渠4内において基準周波数を設定する初期設定を行う。すなわち、測定対象となる菅渠4の内部において、従来公知の水位計を用いて菅渠4内の所定状態の流水43の水位を測定するとともに、そのときの音響波形を測定する。その後、測定された音響波形に対して、基準周波数設定手段としての、例えば携帯用の情報処理端末図示せずや情報処理装置21を用いてFFT解析が行われる。このFFT解析によって、図1に示すような周波数ごとの音圧値のデータの集合、すなわち周波数に対する音圧のプロファイルが導出される。そして、この周波数に対する音圧のプロファイルの導出結果に基づいて、音圧値が極大となる周波数を基準周波数とする。なお、この菅渠4において設定された所定の水位と、この所定の水位に対応する基準周波数の情報データは、図2に示す情報処理装置21の情報記録部21bにデータテーブルとして格納される。このようにして測定対象となる流水43が存在する菅渠4に関する初期設定が行われた後、図3に示すステップST2に移行する。ステップST2においては、菅渠4内の流水43から発せられる音響波形の測定を開始する。すなわち、図2に示すように、音響波形測定手段としてのマイクロホン41によって菅渠4内の音を集音する。マイクロホン41は、集音した音を電気信号に変換して音響信号変換装置32に供給する。音響信号変換装置32は、マイクロホン41からの電気信号を増幅して音響信号として伝送装置32aに供給する。伝送装置32aは、ネットワーク網1を介して音響信号を処理場管理室2の情報処理装置21に供給する。情報処理装置21においては、供給された音響信号を音響波形のデータとして情報記録部21bに格納する。その後、ステップST3に移行する。ステップST3においては、情報処理装置21に供給されて情報記録部21bに格納された音響波形のデータに対して、情報処理部21aによりFFT解析が行われる。すなわち、測定された音響波形のデータは時間軸に沿った音の変化のデータであることから、この音響波形のデータに対して所定時間ごと、例えば5分間ごとに区切って高速フーリエ変換FFTを行うことによって、音響波形から周波数に対する音圧のプロファイルが導出される。この導出された周波数に対する音圧のプロファイルのデータは情報記録部21bに格納される。その後、ステップST4に移行する。ステップST4においては、測定に基づいて算出された周波数と基準周波数との周波数差が算出される。すなわち、周波数差算出手段としての情報処理装置21において、情報記録部21bに記録された測定された周波数と基準周波数との周波数差が算出される。その後、ステップST5に移行する。ステップST5においては、水位算出手段としての情報処理装置21において、算出された周波数差に基づいて菅渠4の内部の流水43の水位が算出される。すなわち、情報処理装置21において、情報処理部21aが、算出した周波数差に基づいて情報記録部21bに格納されたデータテーブルから周波数差に応じた水位の加減値を抽出し、基準周波数に対応した流水43の水位に対して、抽出した加減値を加算する演算を行う。これにより、菅渠4内の流水43の水位が算出される。その後、図2に示す水位測定システムによる菅渠4内の音響波形の測定が継続して行われる。すなわち、ステップST1において初期設定がされた後は、ステップST2~ST5が順次繰り返し実行され、菅渠4内においてマイクロホン41により音響波形が測定されつつ、情報処理装置21において流水43の水位が継続して算出される。他の配管からの流入および段差での水流の落水また、上述した菅渠4内においては、段差が設けられる地点も存在する。図4Aは、この段差地点を示す模式図である。具体的には、菅渠4内においては、マンホール3などの様々な所から流水43が流入する段差地点や段差により流水43の落水が生じる段差地点が存在する。このような段差地点においても流水43の水位を計測する。以下に、この菅渠4内の段差地点における水位測定方法について説明する。図4Aに示すように、マンホール3内の他の配管44から菅渠4内に流水43が流入したり菅渠4内において流水43の落水が生じたりする地点図4A中、点線丸囲み部においては、流水43は常時落水している。このとき流水43は、落水することによって音を発する。ここで、本発明者の知見によれば、菅渠4内の段差地点における流水43は、落下する流水43の流量に応じて音圧が変化する。すなわち、図4Bに示すように、段差地点における流水43の流量が増加すると、所定周波数f2における音圧Lp0′が上昇して音圧Lp1′となり、反対に段差地点における流水43の流量が減少すると、所定周波数f2における音圧Lp0′が下降して音圧Lp2′になるという現象が生じる。また、所定周波数f2近傍においても、その音圧Lpが変化する。図5は、この一実施形態における段差部分での水位を測定する方法を示すフローチャートである。図5に示すように、まず、ステップST1において、従来公知の流量計を用いて菅渠4内の段差地点における流水43の流量を測定するとともに、そのときの音響波形を測定する。その後、測定された音響波形に対して、基準音圧設定手段としての、例えば携帯用の情報処理端末図示せずや情報処理装置21を用いてFFT解析が行われる。これにより、段差地点における流水43の発する音響波形に関する、周波数に対する音圧のプロファイルが導出される。そして、情報記録部21bに、例えば音圧値が極大となる周波数のデータを所定周波数f2として格納するとともに、この所定周波数f2における音圧Lpのデータを基準音圧Lp0′のデータとして格納する。これにより、情報処理装置21において所定周波数f2における基準音圧Lp0′が設定される。その後、ステップST2に移行する。続いて、ステップST2において、流水43が流れている任意状態において、マイクロホン41により段差を落下する流水43が発する音響波形を測定する。測定された音響波形は、マイクロホン41により電気信号とされ音響信号変換装置32に供給される。この電気信号は音響信号変換装置32により音響信号に変換された後、この音響信号は音響波形のデータとして、伝送装置32aおよびネットワーク網1を通じて情報処理装置21に供給される。その後、ステップST3に移行する。ステップST3においては、情報処理装置21に供給された音響信号の音響波形に対して、情報処理部21aによりFFT解析が実行されて、図4Bに示すような周波数に対する音圧のプロファイルが導出される。次に、図5に示すステップST4において、測定された音響波形から導出された周波数に対する音圧Lpのプロファイルと、ステップST1において設定された所定周波数f2における基準周波数Lp0′とに基づいて、音圧差算出手段としての情報処理部21aにより所定周波数f2における音圧差が算出される。その後、ステップST5において、情報処理部21aによって算出された音圧差に基づいて、情報記録部21bに格納された基準流量からの加減値のデータベースから、音圧差に対応する流量の加減値のデータが抽出される。そして、情報処理部21aによって、抽出された流量の加減値のデータと基準流量のデータとから、測定された音圧に対応する流量が算出される。このようにして、情報処理部21aにより、段差地点において生じる音の音響波形から、流水43の流量が算出される。一方、配管44の高さ、菅渠4の勾配のデータ、材質抵抗のデータおよび口径のデータなどの菅渠4における既知の各種データが、データテーブルとして情報記録部21bに格納されている。そして、情報処理部21aにより、下記の1式に従って、菅渠4における勾配、材質および口径のデータと、算出された流量とに基づいて、流水43の流速や水位が算出される。なお、1式において、V:平均流速、R:径深水理学的平均水深、I:動水水面勾配、n:マニング粗度係数菅渠の内周面の材質に固有の係数、A:通水断面積、S:潤辺流れの横断面で水に接している区間の長さである。ステップST5が終了すると、菅渠4内における段差地点においてもマイクロホン41による音響波形の測定が継続して行われる。すなわち、図5に示すステップST1において初期設定がされた後は、ステップST2~ST5が順次繰り返し実行され、マイクロホン41による音響波形の測定が継続して実行されるとともに、情報処理装置21による菅渠4内の段差地点における流水43の水位の算出が実行される。なお、他の配管44からの流水43の流入や段差地点などが複数存在する場合、測定される音響波形は、それぞれの他の配管44からの流入によって発生する音の音響波形と、それぞれの段差における流水43の落水によって発生する音の音響波形とが混合した音響波形となる。そこで、個々の他の配管44ごと、および個々の段差ごとにそれぞれ固有の所定周波数f2を測定して設定し、それぞれ基準周波数Lp0′を測定して設定する。これにより、合成された音響波形から他の配管44からの流入ごと、および段差ごとに、流水43の流量を個別に算出することができる。気泡発生装置次に、この一実施形態による気泡発生装置の詳細について説明する。図6A、図6Bおよび図6Cはそれぞれ、この一実施形態による気泡発生装置42の平面図、側面図および正面図である。図6A~図6Cに示すように、この一実施形態による気泡発生装置42は、例えば防錆加工処理された銅板や鉄板などの基台420上に、同様に防錆加工処理された突起部421が複数並べて設けられている。突起部421は、楕円球を長軸に沿って2等分割した形状において、その長軸に沿った一端部分に切り欠けられた窪み部421aが形成され、例えば蝋付けなどによって基台420上に接着されている。また、基台420には、例えば皿ねじを用いて菅渠4の壁面に設置するための、複数のねじ穴422が形成されている。このような構成によって、気泡発生装置42の強度が高くされて耐腐食性が向上されている。ここで、気泡発生装置42の寸法の一例を挙げると、図6Aに示すように、長さL0が100mm、幅W0が70mm、突起部421の長手方向に沿った長さl2が40mm、突起部421の長手方向に沿ったねじ穴422の形成領域の長さl1l3が30mmである。また、図6Bに示すように、基台420の板厚d0は例えば2mmであり、図6Cに示すように突起部421の高さh0は例えば10mm、幅w0は例えば40mmである。そして、以上のように構成された気泡発生装置42は、図6Bに示すように、突起部421の窪み部421aの形成位置に対して反対側から水流が到達するように設置される。これにより、気泡発生装置42に流水43が到達すると、流水43が窪み部421aに巻き込まれて窪み部421aの内側に渦流が生じ流水43中に気泡が発生する。これにより、菅渠4内において気泡が発生しにくい場所に気泡発生装置42を設置することによって、流水43が到達した時に意図的に気泡を発生させることが可能になる。また、気泡発生装置42を、流水43が窪み部421aに対して突起部421の上部分から巻き込むようにして設置することにより、流水43に伴って流れる砂、泥、または夾雑物などを、特に突起部421や窪み部421aに堆積させないようにすることができる。以上のように構成された気泡発生装置42によって発生される気泡群43aは、固有の周波数の音を発する。気泡発生装置42によって生じる気泡群43aから発せられる固有の周波数は、気泡発生装置42における突起部421や窪み部421aの数、形状、または寸法を変更することによって、異なる周波数にすることができる。そこで、複数の気泡発生装置42を用いて、流水43の水位を測定する方法について以下に説明する。ここで、複数の気泡発生装置42は、突起部421や窪み部421aの数、形状、または寸法を変更して、生成される気泡群から発する音の固有の周波数が相互に異なるように製造される。図7は、相互に異なる固有の周波数の音を発する気泡群43aを発生するように製造された複数の気泡発生装置42a~42dを、菅渠4内に設置した状態を示す。図7に示すように、菅渠4の内壁面に、例えば4台の気泡発生装置42a~42dが水位の増減方向である縦方向に沿って配置されている。具体的には、菅渠4の底部H0を基準として、気泡発生装置42aが高さH1、気泡発生装置42bが高さH2、気泡発生装置42cが高さH3、気泡発生装置42dが高さH4の位置に配置される。なお、これらの気泡発生装置42a~42dの設置位置、設置高さ、およびマイクロホン41からの距離などのデータは、情報処理装置21の情報記録部21bにデータベースとして格納される。そして、菅渠4内に雨水などの水が流れ込み、流水43が例えば高さH1まで増加すると、気泡発生装置42aに流水43が接触して、気泡群43aが発生し音を発する。この気泡群43aの発する音は、図2に示すマイクロホン41によって音響波形の一部として測定される。マイクロホン41によって測定された音響波形は、図2に示す音響信号変換装置32によって音響信号に変換された後、伝送装置32aを介して情報処理装置21に供給されてFFT解析が行われる。これにより、情報処理装置21において、測定された音響波形における周波数に対する音圧のプロファイルが導出される。この周波数に対する音圧のプロファイルにおいては、所定周波数としての気泡発生装置42aの固有の周波数において音圧が極大となる。そして、情報処理部21aによって、情報記録部21bに格納されたデータベースから、気泡発生装置42aの固有の周波数および設置高さH1のデータが抽出され、菅渠4内の流水43の水位が高さH1であることが算出される。なお、算出された径深の値、および情報記録部21bに格納された菅渠4のマニング粗度係数や動水勾配のデータに基づいて、必要に応じて、情報処理部21aによって上述した1式に基づいた演算が行われて流水43の流速が算出される。その後、流水43が増水したりして、流水43の水位が高さH2、高さH3、または高さH4となった場合においても同様にして水位が算出される。すなわち、マイクロホン41が、気泡発生装置42b~42dによって流水43中に生じる気泡群43aの発する音の音響波形を測定するとともに、この音響波形に対してFFT処理を行うことにより、情報処理装置21が菅渠4内の流水43の水位をそれぞれ算出する。以上のようにして流水43の水位が計測されて、さらに流速も計測されると、これらの計測値と菅渠4の管径や形状のデータとに基づいて、情報処理装置21によって菅渠4内の流水43の流量も併せて算出することが可能となる。気泡発生装置の他の実施例1次に、気泡発生装置の他の実施例1について説明する。図8は、この一実施形態の他の実施例1における気泡発生装置を示す。図8に示すように、この一実施形態の他の実施例1による気泡発生装置42eは、円柱状または円筒状の部材の外周曲面に複数の突起部が設けられている。また、気泡発生装置42eは、円柱状または円筒状の円断面のほぼ中心を軸として回転可能に構成されている。すなわち、気泡発生装置42eは、いわゆる羽根からなる突起部を有する水車から構成されている。この気泡発生装置42eは、その突起部の羽面に流水43が接触することにより軸を中心に回転する。気泡発生装置42eの突起部の羽面が流水43と衝突することによって、気泡が発生するとともに衝突音が連続して発生する。そこで、この突起部と流水43との衝突音の音響波形をマイクロホン41により測定することによって、流水43の流れに起因した気泡発生装置42eの回転により生じる音響波形を測定することができる。また、流水43の流速の変化に関連して、気泡発生装置42eの単位時間あたりの回転数が変化する。そこで、気泡発生装置42eの突起部と流水43との接触によって連続して生じる接触音の発生時間間隔を測定する。一方、気泡発生装置42eの寸法や突起部の数に基づいて、その回転数と流水43の流速との相関関係が情報処理装置21の情報記録部21bにデータテーブルとして必要に応じて抽出可能に格納されている。そのため、マイクロホン41により、気泡発生装置42eによって生じる接触音の音響波形を測定して発生時間間隔を計測するとともに、情報処理部21aによって、情報記録部21bに格納されたデータテーブルから情報データを抽出して、演算を行うことにより、情報処理装置21において流水43の流速を算出することが可能となる。また、流水43の水位の上昇に応じて、気泡発生装置42eの突起部421によってかき乱される水量が変化して、気泡発生装置42eから発せられる音響波形の特性も水量に応じて変化する。そこで、気泡発生装置42eから発せられる音響波形と流水43の水量との相関をあらかじめ測定し、水量と水位と音響波形との相関を、情報処理装置21の情報記録部21bのデータベースに格納しておく。これにより、気泡発生装置42eから発せられる音響波形から、流水43の水位を測定することが可能となる。気泡発生装置の他の実施例2次に、気泡発生装置の他の実施例2について説明する。図9は、この一実施形態の他の実施例2における気泡発生装置を示す。図9に示すように、この一実施形態の他の例2による気泡発生装置42fは、中空の細い導管から構成されている。この気泡発生装置42fの導管の一端は流水43中に浸漬されるとともに、他端が大気圧の外気に解放されている。そして、気泡発生装置42fの導管の一端を流水43の下流側に向けることで、他端から大気圧によって導管内の空気が押圧されて、流水43中に空気が放出され、気泡群を発生させることができる。上述したように、気泡発生装置42fによって流水43中に気泡群を発生させることで、マイクロホン41によって気泡群の発する音響波形を測定することができる。そのため、この実施例2においても、マイクロホン41により測定された音響波形のデータに基づいて、情報処理装置21において流水43の流速や水位を算出することが可能となる。以上説明した本発明の一実施形態によれば、菅渠4内に流れる流水43の気泡から発する音響波形を測定し、周波数に対する音圧のプロファイルを導出して、音圧が極大となる周波数の変化に基づいて流水43の水位を算出していることにより、菅渠4内における流水43の水位の計測を非接触で行うことが可能になるので、従来のように接触型の流速計をマンホール3ごとに設ける必要がなくなり、水位測定における機材の設置数を大幅に低減することができる。そのため、機材の設置に要するコストやメンテナンスコストを低減できる。また、従来、水位計がマンホール3ごとに設けられていたことから、水位計から情報処理装置21に送信する情報データ量が極めて多量であったのに対し、この一実施形態においては、測定装置の設置台数を大幅に低減できるので、通信量通信トラフィックも大幅に低減することができる。以上、本発明の一実施形態について具体的に説明したが、本発明は、上述の一実施形態に限定されるものではなく、本発明の技術的思想に基づく各種の変形が可能である。例えば、上述の一実施形態において挙げた数値はあくまでも例に過ぎず、必要に応じてこれと異なる数値を用いてもよい。上述の一実施形態においては、FFT解析を、音響波形における時間を5分間ごとに区切って行っているが、必ずしも5分間ごとに限定されるものではなく、菅渠4における種々の状況や状態に応じて種々の時間に設定することが可能である。また、FFT解析は、対象時間を区切ることなく連続的に行うことも可能である。上述の一実施形態においては、流量の測定時に採用する所定周波数f2として、周波数に対する音圧のプロファイルにおいて音圧が極大のピークになる周波数を採用しているが、必ずしもこれに限定されるものではない。所定周波数f2は、水流の流量との相関が明確となる周波数帯であれば、周波数に対する音圧のプロファイルにおいて音圧が極大になる周波数以外の周波数を採用することも可能である。上述の一実施形態においては、マンホール3の内部に設置される音響信号変換装置32と、菅渠4の内部に設置されるマイクロホン41とをケーブルなどの有線によって接続しているが、必ずしも有線で接続する方法に限定されるものではない。音響信号変換装置32とマイクロホン41とは、必要に応じて無線通信可能な構成を採用することも可能である。