(57)【要約】【課題】  水害のリスクを定量的にリアルタイムで評価できるようにする。【解決手段】  河道網作成装置１０及び流域探知装置１４により、各メッシュ化標高データの各メッシュについて、上流メッシュを示す上流メッシュリストと、各上流メッシュからの流下距離のリストを有する流域データ（１６）を作成する。最大到達時間決定装置２４が、各メッシュについて上流メッシュからの最大到達時間を算出する。リスク算定装置２８が、過去の降雨量データから水害リスクポテンシャルを算定する。リスク算定装置３２が現在の降雨からリスクを算定し、評価装置３４が、現在のリスクと水害リスクポテンシャルを比較して、危険度の有無を評価する。【選択図】  図１
2006127156.tif 000002
【特許請求の範囲】
【請求項１】
メッシュ化標高データの各メッシュについて、上流メッシュを示す上流メッシュリストと、各上流メッシュからの流下距離のリストを有する流域データ（１６）を生成する流域データ生成ステップ（１０，１２，１４）と、
各メッシュについて上流メッシュからの最大到達時間を算定する最大到達時間算定ステップ（２０，２４）と、
所望のメッシュについて、その上流メッシュにおける最大到達時間内での降雨量データから水害リスクを評価するリスク評価ステップ（２８，３０，３２）、
とを具備することを特徴とする水害リスク評価方法。
【請求項２】
当該流域データ生成ステップが、
メッシュ化標高データの各メッシュについて、隣接メッシュとの間の落水方向と落水距離を算定する河道網作成ステップ（１０）と、
当該各メッシュの落水方向から各メッシュに対する上流メッシュを探知する流域探知ステップ（１４）
とを具備することを特徴とする請求項１に記載の水害リスク評価方法。
【請求項３】
当該最大到達時間算定ステップが、
当該所望のメッシュに対する水源メッシュを抽出する水源メッシュ抽出ステップ（２０）と、
当該所望のメッシュについて当該水源メッシュからの最大到達時間を算定する到達時間算定ステップ（２４）
とを具備することを特徴とする請求項１又は２に記載の水害リスク評価方法。
【請求項４】
当該到達時間算定ステップは、当該所望のメッシュと当該所望のメッシュの各水源メッシュとの間の勾配と流下距離から到達時間を算定し、算定された到達時間の内の最大値を最大到達時間とすることを特徴とする請求項３に記載の水害リスク評価方法。
【請求項５】
当該最大到達時間算定ステップが、各メッシュについて、各上流メッシュからの到達時間を算定し、算定された到達時間の内の最大値を最大到達時間とすることを特徴とする請求項１又は２に記載の水害リスク評価方法。
【請求項６】
当該リスク評価ステップが、
所望のメッシュについて、その上流メッシュリストと過去の降雨量データとから水害リスクポテンシャルを算定するリスクポテンシャル算定ステップと、
所望のメッシュについて、その上流メッシュにおける当該最大到達時間内の降雨量から水害リスクを算定するリスク算定ステップと、
当該リスク算定ステップの算定値と、当該リスクポテンシャル算定ステップに
よる算定値とを比較し、警報の可否を評価する評価ステップ
とを具備することを特徴とする請求項１乃至５の何れか１項に記載の水害リスク評価方法。
【請求項７】
更に、当該リスク評価ステップの評価結果により危険と判断された場合に、関係機関に警報する警報ステップ（３６）を具備することを特徴とする請求項１乃至５の何れか１項に記載の水害リスク評価方法。
【請求項８】
メッシュ化標高データの各メッシュについて、上流メッシュを示す上流メッシュリストと、各上流メッシュからの流下距離のリストと、上流メッシュからの最大到達時間データとを有する流域データベースと、
当該流域データベースを参照し、所望のメッシュについて、その上流メッシュにおける最大到達時間内での降雨量データから水害リスクを評価するリスク評価手段（２８，３０，３２）、
とを具備することを特徴とする水害リスク評価システム。
【請求項９】
更に、メッシュ化標高データの各メッシュについて、上流メッシュを示す上流メッシュリストと、各上流メッシュからの流下距離のリストを有する流域データ（１６）を生成する流域データ生成手段（１０，１２，１４）と、
各メッシュについて当該最大到達時間を算定する最大到達時間算定手段（２０，２４）
とを具備することを特徴とする請求項８に記載の水害リスク評価システム。
【請求項１０】
当該流域データ生成手段が、
メッシュ化標高データの各メッシュについて、隣接メッシュとの間の落水方向と落水距離を算定する河道網作成手段（１０）と、
当該各メッシュの落水方向から各メッシュに対する上流メッシュを探知する流域探知手段（１４）
とを具備することを特徴とする請求項９に記載の水害リスク評価システム。
【請求項１１】
当該最大到達時間算定手段が、
当該所望のメッシュに対する水源メッシュを抽出する水源メッシュ抽出手段（２０）と、
当該所望のメッシュについて当該水源メッシュからの最大到達時間を算定する到達時間算定手段（２４）
とを具備することを特徴とする請求項９又は１０に記載の水害リスク評価システム。
【請求項１２】
当該到達時間算定手段が、当該所望のメッシュと当該所望のメッシュの各水源メッシュとの間の勾配と流下距離から到達時間を算定し、算定された到達時間の内の最大値を最大到達時間とすることを特徴とする請求項１１に記載の水害リスク評価システム。
【請求項１３】
当該最大到達時間算定手段が、各メッシュについて、各上流メッシュからの到達時間を算定し、算定された到達時間の内の最大値を最大到達時間とすることを特徴とする請求項９又は１０に記載の水害リスク評価システム。
【請求項１４】
当該リスク評価手段が、
所望のメッシュについて、その上流メッシュリストと過去の降雨量データとから水害リスクポテンシャルを算定するリスクポテンシャル算定手段（２８，３０）と、
所望のメッシュについて、その上流メッシュにおける当該最大到達時間内の降雨量から水害リスクを算定するリスク算定ステップ（３２）と、
当該リスク算定手段の算定値と、当該リスクポテンシャル算定ステップによる算定値とを比較し、警報の可否を評価する評価ステップ（３４）
とを具備することを特徴とする請求項８乃至１３の何れか１項に記載の水害リスク評価システム。
【請求項１５】
更に、当該リスク評価手段の評価結果に従い危険と判断された場合に、関係機関に警報する警報手段（３６）を具備することを特徴とする請求項８乃至１４の何れか１項に記載の水害リスク評価システム。
【請求項１６】
メッシュ化標高データの各メッシュについて、上流メッシュを示す上流メッシュリストと、各上流メッシュからの流下距離のリストと、上流メッシュからの最大到達時間データとを有することを特徴とする流域データベース。
【請求項１７】
更に、各メッシュの面積を示すデータを具備することを特徴とする請求項１６に記載の流域データベース。
【請求項１８】
各メッシュについて、水源メッシュであるか否かを示すデータを具備することを特徴とする請求項１６又は１７に記載の流域データベース。
【請求項１９】
メッシュ化標高データから水の流れを示す流域データベースを作成するコンピュータプログラムであって、
メッシュ化標高データの各メッシュについて、隣接メッシュとの間の落水方向と落水距離を算定する河道網作成機能（１０）と、
当該各メッシュの落水方向から各メッシュに対する上流メッシュを探知し、各メッシュについて、上流メッシュを示す上流メッシュリストと、各上流メッシュからの流下距離のリストを有する流域データ（１６）を生成する流域探知機能（１４）
とをコンピュータに実施させることを特徴とする流域データベース作成プログラム。
【請求項２０】
更に、各メッシュについて上流メッシュからの最大到達時間を算定する最大到達時間算定機能（２４）を具備することを特徴とする請求項１９に記載の流域データベース作成プログラム。
【請求項２１】
更に、  当該所望のメッシュに対する水源メッシュを抽出する水源メッシュ抽出機能（２０）を具備することを特徴とする請求項１９又は２０に記載の流域データベース作成プログラム。
【請求項２２】
当該最大到達時間算定機能が、
当該所望のメッシュに対する水源メッシュを抽出する水源メッシュ抽出機能（２０）と、
当該所望のメッシュについて当該水源メッシュからの最大到達時間を算定する到達時間算定機能（２４）
とを具備することを特徴とする請求項２０に記載の流域データベース作成プログラム。
【請求項２３】
当該到達時間算定機能は、当該所望のメッシュと当該所望のメッシュの各水源メッシュとの間の勾配と流下距離から到達時間を算定し、算定された到達時間の内の最大値を最大到達時間とすることを特徴とする請求項２２に記載の流域データベース作成プログラム。
【請求項２４】
当該最大到達時間算定機能が、各メッシュについて、各上流メッシュからの到達時間を算定し、算定された到達時間の内の最大値を最大到達時間とすることを特徴とする請求項２０に記載の流域データベース作成プログラム。
【請求項２５】
メッシュ化標高データの各メッシュについて、上流メッシュを示す上流メッシュリストと、各上流メッシュからの流下距離のリストを有する流域データ（１６）を生成する流域データ生成ステップ（１０，１２，１４）と、
所望のメッシュの上流メッシュでの降雨量データと当該流下距離とから当該所望のメッシュの水量を算定する算定ステップ
とを具備することを特徴とする水害シミュレーション方法。
【請求項２６】
当該流域データ生成ステップが、
メッシュ化標高データの各メッシュについて、隣接メッシュとの間の落水方向と落水距離を算定する河道網作成ステップ（１０）と、
当該各メッシュの落水方向から各メッシュに対する上流メッシュを探知する流域探知ステップ（１４）
とを具備することを特徴とする請求項２５に記載の水害シミュレーション方法。
【請求項２７】
更に、各メッシュについて当該最大到達時間を算定する最大到達時間算定ステップ（２４）を具備し、当該算定ステップは、当該最大到達時間内で、当該上流メッシュでの降雨量データを使用することを特徴とする請求項２５に記載の水害シミュレーション方法。
【請求項２８】
当該最大到達時間算定ステップが、
当該所望のメッシュに対する水源メッシュを抽出する水源メッシュ抽出ステップ（２０）と、
当該所望のメッシュについて当該水源メッシュからの最大到達時間を算定する到達時間算定ステップ（２４）
とを具備することを特徴とする請求項２７に記載の水害シミュレーション方法。
【請求項２９】
当該到達時間算定ステップが、当該所望のメッシュと当該所望のメッシュの各水源メッシュとの間の勾配と流下距離から到達時間を算定し、算定された到達時間の内の最大値を最大到達時間とすることを特徴とする請求項２８に記載の水害シミュレーション方法。
【請求項３０】
当該最大到達時間算定ステップが、各メッシュについて、各上流メッシュからの到達時間を算定し、算定された到達時間の内の最大値を最大到達時間とすることを特徴とする請求項２７に記載の水害シミュレーション方法。
【請求項３１】
当該降雨量データが、予測降雨量データであることを特徴とする請求項２５に記載の水害シミュレーション方法。
【請求項３２】
当該降雨量データが、現在から過去に遡る所定期間内の降雨量データであることを特徴とする請求項２５に記載の水害シミュレーション方法。
【請求項３３】
水害をシミュレートするコンピュータプログラムであって、
メッシュ化標高データの各メッシュについて、上流メッシュを示す上流メッシュリストと、各上流メッシュからの流下距離のリストを有する流域データ（１６）を生成する流域データ生成機能（１０，１２，１４）と、
所望のメッシュの上流メッシュでの降雨量データと当該流下距離とから当該所望のメッシュの水量を算定する算定機能
とをコンピュータに実施させること特徴とする水害シミュレーションプログラム。
【請求項３４】
当該流域データ生成機能が、
メッシュ化標高データの各メッシュについて、隣接メッシュとの間の落水方向と落水距離を算定する河道網作成機能（１０）と、
当該各メッシュの落水方向から各メッシュに対する上流メッシュを探知する流域探知機能（１４）
とを具備することを特徴とする請求項３３に記載の水害シミュレーションプログラム。
【請求項３５】
更に、各メッシュについて当該最大到達時間を算定する最大到達時間算定機能（２４）を具備し、当該算定機能は、当該最大到達時間内で、当該上流メッシュでの降雨量データを使用することを特徴とする請求項３３に記載の水害シミュレーションプログラム。
【請求項３６】
当該最大到達時間算定機能が、
当該所望のメッシュに対する水源メッシュを抽出する水源メッシュ抽出機能（２０）と、
当該所望のメッシュについて当該水源メッシュからの最大到達時間を算定する到達時間算定機能（２４）
とを具備することを特徴とする請求項３５に記載の水害シミュレーションプログラム。
【請求項３７】
当該到達時間算定機能が、当該所望のメッシュと当該所望のメッシュの各水源メッシュとの間の勾配と流下距離から到達時間を算定し、算定された到達時間の内の最大値を最大到達時間とすることを特徴とする請求項３６に記載の水害シミュレーションプログラム。
【請求項３８】
当該最大到達時間算定機能が、各メッシュについて、各上流メッシュからの到達時間を算定し、算定された到達時間の内の最大値を最大到達時間とすることを特徴とする請求項３５に記載の水害シミュレーションプログラム。
【請求項３９】
当該降雨量データが、予測降雨量データであることを特徴とする請求項３３に記載の水害シミュレーションプログラム。
【請求項４０】
当該降雨量データが、現在から過去に遡る所定期間内の降雨量データであることを特徴とする請求項３３に記載の水害シミュレーションプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、降雨による水害のリスクを定量的に評価する方法及びシステム、流域データーベースとその作成プログラム、並びに水害シミュレーション方法及びプログラムに関する。
【背景技術】
【０００２】
降雨、特に集中的な降雨による浸水又は洪水等の水害は、物的及び人的に多大な被害をもたらす。
【０００３】
従来は、各自治体又は関係機関が、定点観測される雨量又は近辺の河川の水位から経験的に洪水の危険を判断し、浸水又は洪水に対する警告を発している。
【０００４】
洪水等をシミュレーション又は評価するシステムが、特許文献１乃至６に記載されている。
【特許文献１】特開２００４−１９７５５４号公報
【特許文献２】特開２００３−２９６５１６号公報
【特許文献３】特開２００３−１４１６７０号公報
【特許文献４】特開２００３−０１４８６８号公報
【特許文献５】特開２００２−２５６５２５号公報
【特許文献６】特開２００２−１０７４６２号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかし、従来、しばしば警告の遅れが指摘されている。これは、降雨による水害のリスクを定量的に把握する手段が開発されていないからである。
【０００６】
本発明は、集中的な降雨による水害のリスクを定量的に評価する方法及びシステムを提示することを目的とする。
【０００７】
本発明はまた、集中的な降雨による水害の程度をシミュレートする水害シミュレーション方法及びプログラムを提示することを目的とする。
【０００８】
本発明は更に、このような水害リスク評価又は水害シミュレートを実現する流域データベースとその作成プログラムを提示することを目的とする。
【課題を解決するための手段】
【０００９】
本発明に係る水害リスク評価方法は、メッシュ化標高データの各メッシュについて、上流メッシュを示す上流メッシュリストと、各上流メッシュからの流下距離のリストを有する流域データを生成する流域データ生成ステップと、各メッシュについて上流メッシュからの最大到達時間を算定する最大到達時間算定ステップと、所望のメッシュについて、その上流メッシュにおける最大到達時間内での降雨量データから水害リスクを評価するリスク評価ステップとを具備することを特徴とする。
【００１０】
本発明に係る水害リスク評価システムは、メッシュ化標高データの各メッシュについて、上流メッシュを示す上流メッシュリストと、各上流メッシュからの流下距離のリストと、上流メッシュからの最大到達時間データとを有する流域データベースと、当該流域データベースを参照し、所望のメッシュについて、その上流メッシュにおける最大到達時間内での降雨量データから水害リスクを評価するリスク評価手段とを具備することを特徴とする。
【００１１】
本発明に係る流域データベースは、メッシュ化標高データの各メッシュについて、上流メッシュを示す上流メッシュリストと、各上流メッシュからの流下距離のリストと、上流メッシュからの最大到達時間データとを有することを特徴とする。
【００１２】
本発明に係る流域データベース作成プログラムは、メッシュ化標高データから水の流れを示す流域データベースを作成するコンピュータプログラムであって、メッシュ化標高データの各メッシュについて、隣接メッシュとの間の落水方向と落水距離を算定する河道網作成機能と、当該各メッシュの落水方向から各メッシュに対する上流メッシュを探知し、各メッシュについて、上流メッシュを示す上流メッシュリストと、各上流メッシュからの流下距離のリストを有する流域データを生成する流域探知機能とをコンピュータに実施させることを特徴とする。
【００１３】
本発明に係る水害シミュレーション方法は、メッシュ化標高データの各メッシュについて、上流メッシュを示す上流メッシュリストと、各上流メッシュからの流下距離のリストを有する流域データを生成する流域データ生成ステップと、所望のメッシュの上流メッシュでの降雨量データと当該流下距離とから当該所望のメッシュの水量を算定する算定ステップとを具備することを特徴とする。
【００１４】
本発明に係る水害シミュレーションプログラムは、水害をシミュレートするコンピュータプログラムであって、メッシュ化標高データの各メッシュについて、上流メッシュを示す上流メッシュリストと、各上流メッシュからの流下距離のリストを有する流域データを生成する流域データ生成機能と、所望のメッシュの上流メッシュでの降雨量データと当該流下距離とから当該所望のメッシュの水量を算定する算定機能とをコンピュータに実施させること特徴とする。
【発明の効果】
【００１５】
本発明によれば、上流域の降雨の影響を定量的に把握でき、より正確に警報を発することができるようになる。これにより、河川の近傍に居住する人々や、仕事、レジャー等で河川付近に滞在する人々が、降雨による河川水位や浸水に対する危険の程度を定量的に認識できるようになる。河川の急激な増水や氾濫による災害を軽減することが可能となる。
【発明を実施するための最良の形態】
【００１６】
以下、図面を参照して、本発明の実施例を詳細に説明する。
【００１７】
図１は、本発明の一実施例の概略構成ブロック図を示す。
【００１８】
河道網作成装置１０は、デジタル標高（ＤＥＭ）データから等緯度及び等経度で区切った図２に示すような複数のメッシュを作成する。例えば、２５０ｍメッシュのＤＥＭデータから１ｋｍ×１ｋｍ程度のメッシュ、即ちメッシュ化標高データを作成する。そして、各メッシュについて隣接メッシュ間の距離Ｄｉｓｔと落水方向Ｆｌｏｗを計算する。機械的な方法に拠らずに、河川位置とその近辺に適宜にメッシュを配置しても良い。
【００１９】
本実施例では、注目メッシュの東、西、南、北、東北、東南、西南及び西北に隣接する合計８個のメッシュを、注目メッシュに対する隣接メッシュであるとする。各メッシュは、東西方向と南北方向の２方向で特定され得るが、説明の便宜上、図３に示すように、本明細書では、注目する地域のメッシュに、西から東に、そして、北から南に一連の番号を付すこととする。図３に示す例では、２５個のメッシュがある。
【００２０】
図２に示すメッシュ分割例では、東西方向に隣接するメッシュとの間の距離はΔｘ、南北方向に隣接するメッシュとの間の距離はΔｙ、斜め方向（東北、東南、西南及び西北）に隣接するメッシュとの間の距離は、
（Δｘ２＋Δｙ２）１／２
である。各メッシュの中央又は平均の座標をメッシュの水平位置（緯度と経度）Ｐｏｓとし、各メッシュの中央又は平均の標高値をメッシュの標高値Ｈとする。各メッシュの面積Ｓは、緯度により異なる。面積値そのものを保持する代わりに、代表的なメッシュの面積値に対する補正係数を各メッシュの面積値の代わりに記憶しても良いことは明らかである。
【００２１】
各メッシュの落水方向Ｆｌｏｗは、次のように決定する。即ち、注目するメッシュに対する８個の隣接メッシュに対して、標高差とメッシュ間距離から勾配θを、下記式
θ＝（標高差）／（メッシュ間距離）      （１）
に従い、計算する。最も勾配の大きな隣接メッシュの方向を落水方向とする。落水方向Ｆｌｏｗには、落水方向の隣接メッシュを示す番号を入力する。８個の隣接メッシュには、注目メッシュの東北に位置する隣接メッシュの番号を１とし、時計方向に順番に番号を付ける。注目メッシュの北に位置する隣接メッシュの番号は８となる。どの隣接メッシュよりも注目メッシュの標高が低い場合、落水方向Ｆｌｏｗには０を代入する。例えば、窪地でこのような状況になりうる。流れ込んだ水は、何れは外部に流出すると見るのが自然である。従って、注目メッシュと隣接メッシュとの間で標高値の差が所定値以下の場合には、実質的に同じ標高値であるとして処理する。又は、このような窪地の全体を、平均的な標高値を持つ１つにメッシュに統合する。
【００２２】
図４は、河道網作成装置１０の処理のフローチャート例を示す。対象地域を図２に示すように、複数のメッシュに分割し（Ｓ１）、各メッシュの水平位置Ｐｏｓ、面積Ｓ及び標高Ｈを上述のように決定する（Ｓ２）。
【００２３】
また、各メッシュについて、落水方向Ｆｌｏｗと落水距離ＤｏｗｎＤｉｓｔを決定する（Ｓ３?Ｓ８）。即ち、注目メッシュを示すループ変数ｉを１で初期化する（Ｓ３）。注目メッシュｉに対して、８つの隣接メッシュに対する距離Ｄｉｓｔ（ｉ，ｍ）を上述のように計算し（Ｓ４）、隣接メッシュとの標高差を隣接メッシュとの間の距離で除算して勾配を算出する（Ｓ５）。最も勾配の大きい隣接メッシュを示す番号を落水方向Ｆｌｏｗ（ｉ）に代入し、この隣接メッシュとの間の距離を落水距離ＤｏｗｎＤｉｓｔ（ｉ）に代入する（Ｓ６）。ループ変数ｉをインクリメントし（Ｓ７）、ｉがメッシュ数Ｎを越えるまで（Ｓ８）、ステップＳ４〜Ｓ７を繰り返す。
【００２４】
落水方向を決定するより簡易な方法として、隣接メッシュの内で最も標高の低いものを、落水方向としても良い。更に別の方法として、注目メッシュを挟む両側の隣接メッシュの標高差を算出し、その標高差が最も大きい２つの隣接メッシュの内の低いほうを、落水方向としても良い。
【００２５】
このようにして、各メッシュｉについて、水平位置Ｐｏｓ（ｉ），面積Ｓ（ｉ），標高Ｈ（ｉ），落水方向Ｆｌｏｗ（ｉ）及び落水距離ＤｏｗｎＤｉｓｔ（ｉ）を具備する河道網データ１２が得られる。河道網データ１２は、ハードディスク等に記憶される。
【００２６】
図５は、標高値例とこの標高値例から算出される落水方向の例を示す。図５（ａ）が、標高値の例を示し、図５（ｂ）は、落水方向Ｆｌｏｗの値例を示す。図５（ｃ）は、図５（ｂ）に示す落水方向を矢印で示す。図５（ｄ）は、落水距離ＤｏｗｎＤｉｓｔを示す。図５（ｄ）に示す距離例では、東西方向のメッシュのサイズを１．１とし、南北方向のメッシュのサイズを０．９としている。
【００２７】
流域探知装置１４は、河道網データ１２から各メッシュの上流域を探知する。即ち。流域探知装置１４は、河道網データ１２の落水方向Ｆｌｏｗから各メッシュに流れ込む水を供給するメッシュ（上流メッシュ）を決定する。図６は、流域探知装置１４の処理のフローチャート例を示す。なお、上流メッシュリストＵｐＭｅｓｈ（ｉ）は、メッシュｉの上流メッシュのリストを保有し、流下距離リストＦｌｏｗＤｉｓｔ（ｉ）は、ＵｐＭｅｓｈ（ｉ）にリストされる各上流メッシュへの流下距離、即ち水の流れに沿った距離のリストを保有し、ＵｐＮｕｍ（ｉ）は、ＵｐＭｅｓｈ（ｉ）に収容される上流メッシュの数を保有する。この実施例では、上流メッシュに自メッシュを含める。
【００２８】
上流メッシュリストＵｐＭｅｓｈ（ｉ）と流下距離リストＦｌｏｗＤｉｓｔ（ｉ）の同じ順番位置に、それぞれ、互いに対応する上流メッシュの番号と、その上流メッシュまでの流下距離が格納される。例えば、上流メッシュリストＵｐＭｅｓｈ（ｉ）のｎ番目にリストされる上流メッシュまでの流下距離が、流下距離リストＦｌｏｗＤｉｓｔ（ｉ）のｎ番目にリストされる。
【００２９】
先ず、各メッシュについて、ＵｐＭｅｓｈ（ｉ）、ＦｌｏｗＤｉｓｔ（ｉ）及びＵｐＮｕｍ（ｉ）を初期化する（Ｓ２１）。本実施例では、上流メッシュに自メッシュを含めるので、各メッシュｉについて、クリアしたＵｐＭｅｓｈ（ｉ）に自メッシュを示すｉを追加し、クリアしたＦｌｏｗＤｉｓｔ（ｉ）に距離０を追加し、上流メッシュ数ＵｐＮｕｍ（ｉ）に１を代入する。ＵｐＭｅｓｈ（ｉ）とＦｌｏｗＤｉｓｔ（ｉ）を１つの構造体に含まれる要素として対で管理するのが便利である。
【００３０】
注目メッシュｉの８つの隣接メッシュの内で、注目メッシュｉに落水する隣接メッシュを探索し、その隣接メッシュの上流メッシュリストを、注目メッシュｉに追加する。そのために、先ず、注目メッシュを示すループ変数ｉを１で初期化し（Ｓ２１）、隣接メッシュを示す変数ｊを１で初期化する（Ｓ２３）。隣接メッシュｊの落水方向Ｆｌｏｗを調べて、注目メッシュｉに落水するか否かを調べる（Ｓ２４）。当然であるが、隣接メッシュｊが存在しない場合、ステップＳ２４では、隣接メッシュｊから注目メッシュｉへの落水が無いと判断される。隣接メッシュｊから注目メッシュｉに落水していない場合（Ｓ２４）、次の隣接メッシュを調べる為にステップＳ２８に移行する。
【００３１】
隣接メッシュｊから注目メッシュｉに落水している場合（Ｓ２４）、隣接メッシュｊの上流メッシュリストＵｐＭｅｓｈを注目メッシュｉの上流メッシュリストＵｐＭｅｓｈ（ｉ）に追加する（Ｓ２５）。そして、隣接メッシュｊの流下距離リストＦｌｏｗＤｉｓｔの各距離データに、隣接メッシュｊと注目メッシュｉとの間の距離を加算して、注目メッシュｉの流下距離リストＦｌｏｗＤｉｓｔ（ｉ）に追加する（Ｓ２６）。隣接メッシュｊの上流メッシュ数ＵｐＮｕｍを注目メッシュｉの上流メッシュ数ＵｐＮｕｍ（ｉ）を加算する（Ｓ２７）。
【００３２】
ステップＳ２４又はＳ２７の後、隣接メッシュを示すループ変数ｊをインクリメントする（Ｓ２８）。注目メッシュｉの全隣接メッシュを調べ終えるまで（Ｓ２９）、ステップＳ２４〜Ｓ２７を繰り返す。
【００３３】
注目メッシュｉの全隣接メッシュを調べ終えたら（Ｓ２９）、変数ｉをインクリメントする（Ｓ３０）。ｉがメッシュ数Ｎを越えるまで（Ｓ３１）、ステップＳ２３〜Ｓ３０の処理を繰り返す。
【００３４】
このような処理により、降雨により注目メッシュに流れ込む上流の流路上に位置するメッシュ、即ち上流メッシュのリストが、上流メッシュリストＵｐＭｅｓｈ（ｉ）に格納され、各上流メッシュからの流下距離のリストが流下距離リストＦｌｏｗＤｉｓｔ（ｉ）に格納される。各メッシュｉの上流メッシュリストＵｐＭｅｓｈ（ｉ）、流下距離リストＦｌｏｗＤｉｓｔ（ｉ）及び上流メッシュ数ＵｐＮｕｍ（ｉ）からなる流域データ１６が、ハードディスクに格納される。
【００３５】
上流メッシュリストＵｐＭｅｓｈ（ｉ）にリストされる上流メッシュの面積Ｓ（ｉ）の累計を計算することで、流域の面積、即ち、集水面積を決定できる。
【００３６】
図７は、図５（ｃ）に示す落水方向例に対する各メッシュの上流メッシュリストＵｐＭｅｓｈのリスト例を示す。主としてメッシュ１３，１８，２３に水が流れ込んでいることが分かる。
【００３７】
図６では、上流から下流に移動したが、下流から順に上流に遡っても良い。例えば、注目メッシュに対する８つの隣接メッシュの中で注目メッシュに落水する隣接メッシュを探索する処理を、探索された隣接メッシュを注目メッシュとする再帰的に実行することで、下流から上流に向かって上流メッシュを探索できる。
【００３８】
また、注目メッシュを順次、変更しながら、注目メッシュから出た水の落水方向をたどり、注目メッシュと途中の全メッシュのＵｐＭｅｓｈに注目メッシュのメッシュ番号を追加する。流下距離は別途、計算するか、下流に下りつつ、メッシュ間距離を累積していけばよい。
【００３９】
水源探知装置１８は、流域データ１６にリストされる上流メッシュの内から水源メッシュを探知し、水源メッシュには水源フラグＳｒｃ＿Ｆ（ｉ）に１をセットし、水源メッシュ以外のメッシュには、水源フラグＳｒｃ＿Ｆ（ｉ）に０をセットする。水源メッシュは、何れの隣接メッシュからも落水しないメッシュをいう。図８は、水源探知装置１８の処理のフローチャート例を示す。本実施例では、上流メッシュ数ＵｐＮｕｍが１であるメッシュは、水源メッシュである。
【００４０】
注目メッシュを示すループ変数ｉを１で初期化し（Ｓ４１）、上流メッシュ数ＵｐＮｕｍ（ｉ）が１かどうかを調べる（Ｓ４２）。ＵｐＮｕｍ（ｉ）が１である場合（Ｓ４２）、水源フラグＳｒｃ＿Ｆ（ｉ）に１をセットし（Ｓ４３）、ＵｐＮｕｍが１でない場合（Ｓ４２）、水源フラグＳｒｃ＿Ｆ（ｉ）に０をセットする（Ｓ４４）。変数ｉをインクリメントする（Ｓ４５）。ｉがメッシュ数Ｎを越えるまで（Ｓ４６）、ステップＳ４２〜Ｓ４５の処理を繰り返す。
【００４１】
なお、上流メッシュリストＵｐＭｅｓｈに自メッシュを含めない場合には、上流メッシュ数ＵｐＮｕｍが０であるメッシュが水源メッシュである。
【００４２】
次に、各メッシュに対して、水源メッシュからの最大到達時間を算出する。最大到達時間を算出することで、降雨に対する時間的な影響範囲を制限できる。最大到達時間よりも前の降雨による水は、既に注目メッシュを通過済みだからである。
【００４３】
注目メッシュから流下時間で最も遠い水源メッシュを特定する為に、水源メッシュ抽出装置２０が、上流メッシュリストＵｐＭｅｓｈから水源メッシュを抽出する。図９は、水源メッシュ抽出装置２０の処理のフローチャートを示す。
【００４４】
処理対象のメッシュを指定する変数ｉを初期化する（Ｓ５１）。水源メッシュ抽出装置２０は、流域データ１６から上流メッシュリストＵｐＭｅｓｈ（ｉ）を読み出し、そのコピーを水源メッシュリストＵｐＳｒｃＭｅｓｈ（ｉ）とし（Ｓ５２）、流域データ１６から流下距離リストＦｌｏｗＤｉｓｔ（ｉ）を読み出し、そのコピーを水源流下距離リストＵｐＳｒｃＤｉｓｔ（ｉ）とする（Ｓ５３）。
【００４５】
上流メッシュリストＵｐＭｅｓｈ（ｉ）にリストされている各メッシュ番号を指定する変数ｊを１で初期化し（Ｓ５４）、水源メッシュリストＵｐＳｒｃＭｅｓｈ（ｉ）からｊ番目のメッシュ番号を読み出し、変数ｋに収容する（Ｓ５５）。ｋがｉに等しければ（Ｓ５６）、それは自メッシュであるので、次のメッシュの処理（Ｓ５９）に移行する。自メッシュでない場合（Ｓ５６）、メッシュｋが水源メッシュか否か、即ち、水源フラグＳｒｃ＿Ｆ（ｋ）が１か否かを調べる（Ｓ５７）。メッシュｋが水源メッシュでなければ（Ｓ５７）、ＵｐＳｒｃＭｅｓｈ（ｉ）及びＵｐＳｒｃＤｉｓｔ（ｉ）からｊ番目のデータを削除又はＮＵＬＬとする（Ｓ５８）。ステップＳ５６，Ｓ５７又はＳ５８の後、ｊをインクリメントする（Ｓ５９）。水源メッシュリストＵｐＳｒｃＭｅｓｈ（ｉ）にリストされている全メッシュを調べるまで（Ｓ６０）、ステップＳ５５〜Ｓ５９を繰り返す。
【００４６】
これにより、水源メッシュリストＵｐＳｒｃＭｅｓｈ（ｉ）は、注目メッシュｉの上流に位置する水源メッシュの番号のリストを保持し、水源流下距離リストＵｐＳｒｃＤｉｓｔ（ｉ）は、水源メッシュリストＵｐＳｒｃＭｅｓｈ（ｉ）にリストされる水源メッシュまでの流下距離のリストを、ＵｐＳｒｃＭｅｓｈ（ｉ）と同じ順番で保持する。
【００４７】
メッシュｉの上流メッシュを全て調べたら（Ｓ６０）、ｉをインクリメントする（Ｓ６１）、対象地域の全メッシュを調べ終えるまで（Ｓ６２）、ステップＳ５２〜Ｓ６１を繰り返す。
【００４８】
水源メッシュ抽出装置２０は、各メッシュに対する水源メッシュリストＵｐＳｒｃＭｅｓｈ（ｉ）及び及び水源流下距離リストＵｐＳｒｃＤｉｓｔ（ｉ）を、水源リスト２２としてハードディスク２２に格納する。図１０は、図５（ｃ）に示す落水方向例に対する水源メッシュリストＵｐＳｒｃＭｅｓｈ（ｉ）のリスト例を示す。
【００４９】
次に、最大到達時間決定装置２４が、水源リスト２２の水源メッシュリストＵｐＳｒｃＭｅｓｈ（ｉ），水源流下距離リストＵｐＳｒｃＤｉｓｔ（ｉ）、並びに、河道網データ１２の水平位置Ｐｏｓ（ｉ）及び標高Ｈ（ｉ）から、各メッシュｉについての最大到達時間Ｔｍａｘ（ｉ）を算出及び決定する。図１１は、最大到達時間決定装置２４の処理のフローチャートを示す。
【００５０】
先ず、注目メッシュを示す変数ｉを１で初期化する（Ｓ７１）。注目メッシュｉの水源メッシュリストＵｐＳｒｃＭｅｓｈ（ｉ）にリストされた水源メッシュを示す変数ｊを１で初期化する（Ｓ７２）。
【００５１】
水源メッシュリストＵｐＳｒｃＭｅｓｈ（ｉ）からｊ番目の水源メッシュの番号を読みだし、変数ｋに代入する（Ｓ７３）。水源流下距離リストＵｐＳｒｃＤｉｓｔ（ｉ）からｊ番目の流下距離データを読みだす（Ｓ７４）。メッシュｉ，ｋの標高Ｈ（ｉ），Ｈ（ｋ）とメッシュｉ，ｋ間の直線距離から勾配θを算出する（Ｓ７５）。即ち、θは、
θ＝｜（Ｈ（ｉ）−Ｈ（ｋ））／（Ｐｏｓ（ｉ）−Ｐｏｓ（ｋ））｜    （２）
である。そして、メッシュｋからメッシュｉへの到達時間Ｔｉｍｅ（ｊ）を、下記式に従い、
Ｔｉｍｅ（ｊ）＝α（Ｌ／θｐ）ｑ＋β        （３）
計算する。但し、α、β、ｐ及びｑは、流域の土地利用及び流域面積に応じて設定される係数であり、Ｌは、水源流下距離リストＵｐＳｒｃＤｉｓｔ（ｉ）から読み出される水源メッシュｋから注目メッシュｉへの流下距離、即ち水の流れに沿った距離である。
【００５２】
α、β、ｐ及びｑの各係数値は、例えば、過去の降水観測所の降水量データと下流の水位又は流量データとの関係を、上述のようにして決定された流域データベースに適用することで、算出できる。例えば、土木研究所から、以下のような数値が報告されている。都市流域では、α＝２．４０×１０−４、β＝０、ｐ＝０．５、１＝０．７、自然流域では、α＝１．６７×１０−３、β＝０、ｐ＝０．５、１＝０．７である（「改訂新版  建設省河川砂防技術基準（案）  調査編」建設省河川局監修、（社）日本河川協会編、平成９年９月  山海堂発行）。
【００５３】
ｊをインクリメントする（Ｓ７７）。そして、注目メッシュｉの水源メッシュリストデータＵｐＳｒｃＭｅｓｈ（ｉ）にリストされた全水源メッシュで到達時間Ｔｉｍｅ（ｊ）を計算し終えるまで（Ｓ７８）、ステップＳ７３〜Ｓ７７を繰り返す。
【００５４】
注目メッシュｉの水源メッシュリストデータＵｐＳｒｃＭｅｓｈ（ｉ）にリストされた全水源メッシュで到達時間Ｔｉｍｅ（ｊ）を計算し終えると（Ｓ７８）、ステップＳ７６で計算されたＴｉｍｅ（ｊ）の内の最大値を最大到達時間Ｔｍａｘ（ｉ）とする（Ｓ７９）。
【００５５】
ｉをインクリメントする（Ｓ８０）。全てのメッシュについて最大到達時間Ｔｍａｘ（ｉ）を計算するまで（Ｓ８１）、ステップＳ７２〜Ｓ８０を繰り返す。
【００５６】
以上により、対象地域内の全メッシュについて、水源メッシュからの最大到達時間Ｔｍａｘを算出できる。算出された各メッシュの最大到達時間データ２６は、ハードディスク等に格納される。
【００５７】
最大到達時間Ｔｍａｘ（ｉ）を決定する別の方法として、全上流メッシュについて上式に従い到達時間Ｔｉｍｅを算定し、その内の最大値を最大到達時間Ｔｍａｘとしてもよい。また、簡易的には、直線距離で最も遠い水源メッシュについて、上式に従い到達時間を計算し、その到達時間を最大到達時間Ｔｍａｘ（ｉ）として扱っても良い。
【００５８】
リスク算定装置２８は、任意のメッシュｉについて、上流メッシュリストＵｐＭｅｓｈ（ｉ）と最大到達時間Ｔｍａｘ（ｉ）内の過去の降雨量データとから、メッシュｉとその上流域に最大到達時間内に降った雨量の総計を推計する。推計された総降雨量は、最悪の場合で一度に注目メッシュに流入する水量でもある。従って、この総雨量は、注目メッシュｉにおける水害の危険を定量的に示す指標として利用可能である。この意味で、本実施例では、リスク算定装置２８の算定結果をリスクポテンシャルと呼ぶ。具体的には、リスク算定装置２８は、注目メッシュｉに対するリスクポテンシャルＲＰ（ｉ）を、下記式
ＲＰ（ｉ）＝ΣｊΣｔＲａｉｎ（ｊ，−ｔ）Ｓ（ｊ）        （４）
に従い計算する。但し、Ｒａｉｎ（ｊ，−ｔ）は、上流メッシュｊに対するｔ時間前の降雨量を示す。Σｔは、ｔ＝０から最大到達時間Ｔｍａｘ（ｉ）までの累積加算を示し、Σｊは、上流メッシュｊの累積加算を示す。即ち、現在から最大到達時間前の時間内で、自メッシュを含む全上流メッシュｊに降った雨量の総和を計算する。
【００５９】
例えば、図７に示す例のメッシュ２３に対して、最大到達時間が６０分であるとする。メッシュ１３の上流メッシュは、メッシュ１，２，３，４，６，７，８，９，１１，１２，１３，１４である。メッシュ１３に対するリスクポテンシャルＲＰ（１３）は、上流メッシュ１，２，３，４，６，７，８，９，１１，１２，１３，１４について、注目時点から過去の６０分以内の降雨量を積算した結果となる。勿論、各上流メッシュについて、到達時間内の降雨量を考慮してもよいが、最大のリスクに備えるという意味で、本実施例では、最大到達時間内の上流域の降雨量が全て一時に注目メッシュに集まる状況を想定する。最大到達時間を越えた過去の降雨量は、既に注目メッシュを通過していると考えられるので、リスクとして評価しない。
【００６０】
このようにして計算された各メッシュｉのリスクポテンシャルＲＰ（ｉ）を水害リスクポテンシャルデータ３０としてハードディスク等に保存する。各メッシュｉのリスクポテンシャルＲＰ（ｉ）を対比して、リスクの高い順に順位付けしておく。実際の被害と被害の可能性を考慮して、リスクポテンシャルＲＰの警報に関する閾値を決定する。実際の被害と被害の可能性を考慮して、α、β、ｐ及びｑを修整しても良い。
【００６１】
流域の総雨量の代わりに、面積当たりの平均雨量を計算しても良い。しかし、上流域の面積が広いほど、より多くの雨が注目メッシュに集中するので、危険度が高い。総雨量を計算することで、このようなリスクを正しく評価できる。
【００６２】
リスク算定装置３２は、現在の降雨量データに従い、各メッシュｉのリスクを算定する。現在の降雨量データは、例えば、リアルタイムでレーダアメダスから入手可能である。リスク算定装置３２の処理内容は、リスク算定装置２８と全く同じである。入力する降雨量データが、過去のものか、現在のものかの相違しかない。即ち、リスク算定装置３２は、現在の降雨量データＲａｉｎ（ｉ，ｔ）に対して、リスクＲ（ｉ）を下記式
Ｒ（ｉ）＝ΣｊΣｔＲａｉｎ（ｊ，−ｔ）Ｓ（ｊ）        （５）
に従い算出する。
【００６３】
評価装置３４は、リスク算定装置３２で算定されたリスクＲ（ｉ）を水害リスクポテンシャル３０からのリスクポテンシャルＲＰ（ｉ）と対比して、危険度を評価及び判定する。リスクＲ（ｉ）がリスクポテンシャルＲＰ（ｉ）を上回るとき、評価装置３４は、関係部署又は地域に警報を送信するように警報装置３６に指示する。
【００６４】
また、将来の降雨量の予測値をリスク算定装置３２に入力することで、将来のリスクを定量的に評価できる。これにより、現在の降雨量データを使用する場合よりも更に早期に、危険度を関係機関又は地域住民に警告できる。
【００６５】
河道網作成装置１０、流域探知装置１４、水源メッシュ抽出装置２０、水源探知装置１８、最大到達時間決定装置２４及びリスク算定装置２８，３２の各機能の一部又は全部は、コンピュータ上のソフトウエアプログラムにより実現され得る。
【００６６】
リスクポテンシャルＲＰ又はリスクＲをその量に応じてランク分けてして地図上に重ねて表示することにより、河道に沿った水害リスクを容易に把握できるようになる。
【００６７】
本実施例の流域データベースを利用することで、予想される降雨に対する任意の地域での水量の時間変化をシミュレートすることができる。注目メッシュの各上流メッシュでの降雨は、注目メッシュまでの到達時間の遅れをもって、注目メッシュに到達するので、この遅れを考慮して、注目メッシュ上で累積すればよい。これにより、上流メッシュでの降雨の時間変化に対して、注目メッシュでの水量の時間変化をシミュレートできる。勿論、一定時間内の予想される降雨に対する、上流メッシュから注目メッシュまでの到達時間を無視することで、最大水量を見積もることも可能になる。
【００６８】
特定の説明用の実施例を参照して本発明を説明したが、特許請求の範囲に規定される本発明の技術的範囲を逸脱しないで、上述の実施例に種々の変更・修整を施しうることは、本発明の属する分野の技術者にとって自明であり、このような変更・修整も本発明の技術的範囲に含まれる。
【図面の簡単な説明】
【００６９】
【図１】本発明の一実施例の概略構成ブロック図である。
【図２】メッシュ標高データの構成例である。
【図３】メッシュの番号付け方法の説明図である。
【図４】河道網作成装置１０の処理のフローチャート例である。
【図５】標高値例とこの標高値例から算出される落水方向の例を示す。
【図６】流域探知装置１４の処理のフローチャート例を示す。
【図７】図５（ｃ）に示す落水方向例に対する各メッシュの上流メッシュリストＵｐＭｅｓｈのリスト例を示す。
【図８】水源探知装置１８の処理のフローチャート例を示す。
【図９】水源メッシュ抽出装置２０の処理のフローチャート例を示す。
【図１０】図５（ｃ）に示す落水方向例に対する水源メッシュリストＵｐＳｒｃＭｅｓｈ（ｉ）のリスト例を示す。
【図１１】最大到達時間決定装置２４の処理のフローチャート例を示す。
【符号の説明】
【００７０】
１０：河道網作成装置
１２：河道網データ
１４：流域探知装置
１６：流域データ
１８：水源探知装置
２０：水源メッシュ抽出装置
２２：水源リスト
２４：最大到達時間決定装置
２６：最大到達時間データ
２８：リスク算定装置
３０：水害リスクポテンシャルデータ
３２：リスク算定装置
３４：評価装置
３６：警報装置
Ｆｌｏｗ：落水方向
ＤｏｗｎＤｉｓｔ：落水距離
ＵｐＭｅｓｈ：上流メッシュリスト
ＦｌｏｗＤｉｓｔ：流下距離リスト
Ｓｒｃ＿Ｆ：水源フラグ
ＵｐＳｒｃＭｅｓｈ：上流水源メッシュリスト
ＵｐＳｒｃＤｉｓｔ：上流水源流下距離リスト
Ｔｍａｘ：最大到達時間
【図１】
2006127156.tif 000003
【図２】
2006127156.tif 000004
【図３】
2006127156.tif 000005
【図４】
2006127156.tif 000006
【図５】
2006127156.tif 000007
【図６】
2006127156.tif 000008
【図７】
2006127156.tif 000009
【図８】
2006127156.tif 000010
【図９】
2006127156.tif 000011
【図１０】
2006127156.tif 000012
【図１１】
2006127156.tif 000013