(57)【要約】【課題】本開示は、河川、ダム又は下水等の水位を予測するために、ニューラルネットワークを利用するにあたり、入力データを適切に抽出することにより、水位を十分に精度高く予測し、計算量を低減することを目的とする。【解決手段】本開示は、水位予測地点及びある時刻の水位の学習に必要な、水位予測地点及び当該時刻の水位との相関が高い雨量及び水位の訓練データを抽出する訓練データ抽出ステップＥと、ニューラルネットワークを用いて、水位予測地点及び当該時刻の水位との相関が高い抽出された雨量及び水位の訓練データを学習する訓練データ学習ステップＣ１〜Ｃ４と、を順に備えることを特徴とする水位予測方法である。【選択図】図１
2019095240.tif 000002
【特許請求の範囲】
【請求項１】
水位予測地点及びある時刻の水位の学習に必要な、前記水位予測地点及び当該時刻の水位との相関が高い雨量及び水位の訓練データを抽出する訓練データ抽出ステップと、
ニューラルネットワークを用いて、前記水位予測地点及び当該時刻の水位との相関が高い抽出された雨量及び水位の訓練データを学習する訓練データ学習ステップと、
を順に備えることを特徴とする水位予測方法。
【請求項２】
前記訓練データ抽出ステップでは、前記水位予測地点の水位との相関が高い水位の訓練データとして、前記水位予測地点の水位の訓練データを抽出し、前記水位予測地点の水位との相関が高い雨量の訓練データとして、前記水位予測地点の周辺の地形に基づいて、前記水位予測地点の上流域の雨量の訓練データを抽出する
ことを特徴とする、請求項１に記載の水位予測方法。
【請求項３】
前記訓練データ抽出ステップでは、前記水位予測地点の水位との相関が高い水位の訓練データとして、前記水位予測地点の水位の訓練データと、前記水位予測地点より上流地点の水位の訓練データと、を抽出し、前記水位予測地点の水位との相関が高い雨量の訓練データとして、前記水位予測地点の周辺の地形及び前記上流地点の周辺の地形に基づいて、前記上流地点の上流域を除く前記水位予測地点の上流域の雨量の訓練データを抽出する
ことを特徴とする、請求項１に記載の水位予測方法。
【請求項４】
前記訓練データ抽出ステップでは、学習対象として水位ピーク期間の水位の訓練データを抽出する際に、抽出される水位の訓練データの個数が削減されすぎず、水位不変動期間の水位の訓練データが抽出されないように、水位ピーク判別閾値を設定する
ことを特徴とする、請求項１から３のいずれかに記載の水位予測方法。
【請求項５】
前記訓練データ抽出ステップでは、学習対象として水位ピーク期間の水位の訓練データを抽出する際に、前記水位ピーク期間の全体が抽出され、水位不変動期間の水位の訓練データが抽出されないように、水位ピーク時刻のまわりの水位ピーク抽出幅を設定する
ことを特徴とする、請求項１から４のいずれかに記載の水位予測方法。
【請求項６】
前記訓練データ抽出ステップでは、ある時刻の水位との相関が高い雨量の訓練データとして、当該時刻より以前の所定雨量期間の雨量の訓練データを抽出する際に、雨量の訓練データと水位の訓練データとの相関が最適化されるような、雨量の訓練データと水位の訓練データとの間のシフト時間を、前記所定雨量期間として設定する
ことを特徴とする、請求項１から５のいずれかに記載の水位予測方法。
【請求項７】
前記訓練データ抽出ステップでは、ある時刻の雨量の訓練データとして、当該時刻より以前の所定累加期間の雨量の訓練データの累加結果を採用する際に、前記シフト時間を施された雨量の訓練データと水位の訓練データとの相関が飽和するような、各時刻の雨量の訓練データにおける雨量の訓練データの累加期間を、前記所定累加期間として設定する
ことを特徴とする、請求項６に記載の水位予測方法。
【請求項８】
前記訓練データ抽出ステップでは、ある時刻の水位との相関が高い水位の訓練データとして、当該時刻より過去の所定水位期間の水位の訓練データを抽出する際に、水位の訓練データの上昇傾向又は下降傾向が持続する時間を、前記所定水位期間として設定する
ことを特徴とする、請求項１から７のいずれかに記載の水位予測方法。
【請求項９】
前記水位予測地点及びある時刻の水位の予測に必要な、前記水位予測地点及び当該時刻の水位との相関が高い雨量及び水位の評価データを抽出する評価データ抽出ステップと、
前記ニューラルネットワークを用いて、前記水位予測地点及び当該時刻の水位との相関が高い抽出された雨量及び水位の評価データに基づいて、前記水位予測地点及び当該時刻の水位を予測する水位予測ステップと、
を前記訓練データ学習ステップの後に順に備える
ことを特徴とする、請求項１から８のいずれかに記載の水位予測方法。
【請求項１０】
請求項１から９のいずれかに記載の水位予測方法をコンピュータに実行させるための水位予測プログラム。
【請求項１１】
請求項１０に記載の水位予測プログラムをインストールされた水位予測装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、河川、ダム又は下水等の水位を予測する技術に関する。
【背景技術】
【０００２】
水害の予報又は警報を発令するために、河川、ダム又は下水等の水位を予測する技術が存在する（例えば、非特許文献１〜８を参照。）。水位予測方法として、貯留関数法、タンクモデル又は分布型モデル等の流出モデルを利用する方法と、時系列解析又は機械学習等の統計的手法を利用する方法と、が存在する。機械学習の中でも、ニューラルネットワークを利用する方法では、水位をリアルタイムにある程度精度高く予測することができる。
【０００３】
ニューラルネットワークを利用する方法では、学習時には、雨量及び水位の時系列の訓練データを１次元のアレイ状に全結合レイヤに入力し、雨量及び水位の時系列の訓練データを学習する。そして、評価時には、雨量及び水位の時系列の評価データを１次元のアレイ状に全結合レイヤに入力し、水位予測地点及び水位予測時刻の水位を予測する。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】阿部清明、外３名、「ニューラルネットワークによる流出解析手法（日流量）に関する研究」、土木学会論文集、２０００年８月、Ｎｏ．６５６／ＩＩ−５２、ｐｐ．１−１３
【非特許文献２】平岡透、外２名、「降雨量からの河川水位予測の一手法」、土木学会論文集Ｂ、２０１０年３月、Ｖｏｌ．６６、Ｎｏ．１、ｐｐ．９３−９８
【非特許文献３】高崎忠勝、外２名、「ニューラルネットワークによる都市中小河川感潮域の水位推定」、土木学会論文集Ｂ１（水工学）、２０１１年、Ｖｏｌ．６７、Ｎｏ．４、ｐｐ．Ｉ＿１６０３−Ｉ＿１６０８
【非特許文献４】一言正之、外２名、「ニューラルネットワークを用いた洪水予測システムの開発」、こうえいフォーラム、２０１２年３月、Ｖｏｌ．２０
【非特許文献５】関基、外６名、「河川特性を反映したニューラルネットワーク洪水予測の精度向上」、河川技術論文集、２０１３年６月、Ｖｏｌ．１９
【非特許文献６】椿涼太、外１名、「危機管理の実務に供する洪水予測技術」、［ｏｎｌｉｎｅ］、２０１３年、河川技術に関するシンポジウム、［２０１７年１１月１日検索］、インターネット＜ＵＲＬ：ｈｔｔｐ：／／ｃｏｍｍｉｔｔｅｅｓ．ｊｓｃｅ．ｏｒ．ｊｐ／ｈｙｄｒａｕｌｉｃ０１／ｓｙｓｔｅｍ／ｆｉｌｅｓ／２０１３＿ＯＳ１．ｐｄｆ＞
【非特許文献７】一言正之、外２名、「深層学習を用いた河川水位予測手法の開発」、土木学会論文集Ｂ１（水工学）、２０１６年、Ｖｏｌ．７２、Ｎｏ．４、ｐｐ．Ｉ＿１８７−Ｉ＿１９２
【非特許文献８】一言正之、外１名、「深層ニューラルネットワークと分布型モデルを組み合わせたハイブリッド河川水位予測手法」、土木学会論文集Ｂ１（水工学）、２０１７年、Ｖｏｌ．７３、Ｎｏ．１、ｐｐ．２２−３３
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ニューラルネットワークを利用する方法では、水位予測地点及び水位予測時刻の水位と相関が高いデータを、全結合レイヤに入力するのみならず、水位予測地点及び水位予測時刻の水位と相関が低いデータを、全結合レイヤに入力している。よって、水位を十分に精度高く予測することができず、計算量を低減することができなかった。
【０００６】
そこで、前記課題を解決するために、本開示は、河川、ダム又は下水等の水位を予測するために、ニューラルネットワークを利用するにあたり、入力データを適切に抽出することにより、水位を十分に精度高く予測し、計算量を低減することを目的とする。
【課題を解決するための手段】
【０００７】
前記課題を解決するために、水位予測地点及びある時刻の水位の学習に必要な、水位予測地点及び当該時刻の水位との相関が高い雨量及び水位の訓練データを抽出する。
【０００８】
具体的には、本開示は、水位予測地点及びある時刻の水位の学習に必要な、前記水位予測地点及び当該時刻の水位との相関が高い雨量及び水位の訓練データを抽出する訓練データ抽出ステップと、ニューラルネットワークを用いて、前記水位予測地点及び当該時刻の水位との相関が高い抽出された雨量及び水位の訓練データを学習する訓練データ学習ステップと、を順に備えることを特徴とする水位予測方法である。
【０００９】
この構成によれば、雨量及び水位の訓練データを適切に抽出することにより、雨量及び水位の訓練データを適切に学習し、訓練時の計算量を低減することができる。
【００１０】
また、本開示は、前記訓練データ抽出ステップでは、前記水位予測地点の水位との相関が高い水位の訓練データとして、前記水位予測地点の水位の訓練データを抽出し、前記水位予測地点の水位との相関が高い雨量の訓練データとして、前記水位予測地点の周辺の地形に基づいて、前記水位予測地点の上流域の雨量の訓練データを抽出することを特徴とする水位予測方法である。
【００１１】
この構成によれば、上記の上流域より他の領域の雨量の訓練データを、学習に利用しない（データをゼロクリア等）ことにより、訓練時の計算量を低減することができる。
【００１２】
また、本開示は、前記訓練データ抽出ステップでは、前記水位予測地点の水位との相関が高い水位の訓練データとして、前記水位予測地点の水位の訓練データと、前記水位予測地点より上流地点の水位の訓練データと、を抽出し、前記水位予測地点の水位との相関が高い雨量の訓練データとして、前記水位予測地点の周辺の地形及び前記上流地点の周辺の地形に基づいて、前記上流地点の上流域を除く前記水位予測地点の上流域の雨量の訓練データを抽出することを特徴とする水位予測方法である。
【００１３】
この構成によれば、上記の上流地点の上流域の雨量の訓練データを、上記の上流地点の水位の訓練データに置換することにより、訓練時の計算量を低減することができる。
【００１４】
また、本開示は、前記訓練データ抽出ステップでは、学習対象として水位ピーク期間の水位の訓練データを抽出する際に、抽出される水位の訓練データの個数が削減されすぎず、水位不変動期間の水位の訓練データが抽出されないように、水位ピーク判別閾値を設定することを特徴とする水位予測方法である。
【００１５】
この構成によれば、十分なデータ量の水位変動期間（洪水時等）のみを学習するため、水位の訓練データを適切に学習し、訓練時の計算量を低減することができる。
【００１６】
また、本開示は、前記訓練データ抽出ステップでは、学習対象として水位ピーク期間の水位の訓練データを抽出する際に、前記水位ピーク期間の全体が抽出され、水位不変動期間の水位の訓練データが抽出されないように、水位ピーク時刻のまわりの水位ピーク抽出幅を設定することを特徴とする水位予測方法である。
【００１７】
この構成によれば、上昇開始から下降終了までの水位変動期間（洪水時等）のみを学習するため、水位の訓練データを適切に学習し、訓練時の計算量を低減することができる。
【００１８】
また、本開示は、前記訓練データ抽出ステップでは、ある時刻の水位との相関が高い雨量の訓練データとして、当該時刻より以前の所定雨量期間の雨量の訓練データを抽出する際に、雨量の訓練データと水位の訓練データとの相関が最適化されるような、雨量の訓練データと水位の訓練データとの間のシフト時間を、前記所定雨量期間として設定することを特徴とする水位予測方法である。
【００１９】
降雨地点から水位予測地点までの距離及び傾斜に応じて、降雨が降雨地点から水位予測地点へと到達するために必要である時間が決まる。この構成によれば、降雨が降雨地点から水位予測地点へと到達するために必要である時間分だけ、つまり、学習すべき時刻の水位と相関が高い期間分だけ、雨量の訓練データを抽出することにより、雨量の訓練データを適切に学習し、訓練時の計算量を低減することができる。
【００２０】
また、本開示は、前記訓練データ抽出ステップでは、ある時刻の雨量の訓練データとして、当該時刻より以前の所定累加期間の雨量の訓練データの累加結果を採用する際に、前記シフト時間を施された雨量の訓練データと水位の訓練データとの相関が飽和するような、各時刻の雨量の訓練データにおける雨量の訓練データの累加期間を、前記所定累加期間として設定することを特徴とする水位予測方法である。
【００２１】
降雨地点での土壌内又は貯水池等への貯水量に応じて、どのくらいの時間分だけ降雨が続けば降雨が土壌内又は貯水池等から流出するかが決まる。この構成によれば、どのくらいの時間分だけ降雨が続けば降雨が土壌内又は貯水池等から流出するかに応じて、雨量の訓練データを適切に学習し、訓練時の計算量を低減することができる。
【００２２】
また、本開示は、前記訓練データ抽出ステップでは、ある時刻の水位との相関が高い水位の訓練データとして、当該時刻より過去の所定水位期間の水位の訓練データを抽出する際に、水位の訓練データの上昇傾向又は下降傾向が持続する時間を、前記所定水位期間として設定することを特徴とする水位予測方法である。
【００２３】
降雨地点から水位予測地点までの距離及び傾斜に応じて、水位が上昇又は下降の傾向にある時間が決まる。この構成によれば、水位が上昇又は下降の傾向にある時間分だけ、つまり、学習すべき時刻の水位と相関が高い期間分だけ、水位の訓練データを抽出することにより、水位の訓練データを適切に学習し、訓練時の計算量を低減することができる。
【００２４】
また、本開示は、前記水位予測地点及びある時刻の水位の予測に必要な、前記水位予測地点及び当該時刻の水位との相関が高い雨量及び水位の評価データを抽出する評価データ抽出ステップと、前記ニューラルネットワークを用いて、前記水位予測地点及び当該時刻の水位との相関が高い抽出された雨量及び水位の評価データに基づいて、前記水位予測地点及び当該時刻の水位を予測する水位予測ステップと、を前記訓練データ学習ステップの後に順に備えることを特徴とする水位予測方法である。
【００２５】
この構成によれば、雨量及び水位の評価データを適切に抽出することにより、水位を十分に精度高く予測し、評価時の計算量を低減することができる。
【００２６】
また、本開示は、以上に記載の水位予測方法をコンピュータに実行させるための水位予測プログラムである。
【００２７】
この構成によれば、雨量及び水位の訓練データを適切に抽出することにより、雨量及び水位の訓練データを適切に学習し、訓練時の計算量を低減することができる。
【００２８】
また、本開示は、以上に記載の水位予測プログラムをインストールされた水位予測装置である。
【００２９】
この構成によれば、雨量及び水位の訓練データを適切に抽出することにより、雨量及び水位の訓練データを適切に学習し、訓練時の計算量を低減することができる。
【発明の効果】
【００３０】
このように、本開示は、河川、ダム又は下水等の水位を予測するために、ニューラルネットワークを利用するにあたり、入力データを適切に抽出することにより、水位を十分に精度高く予測し、計算量を低減することができる。
【図面の簡単な説明】
【００３１】
【図１】本開示の水位予測方法の処理手順を示す図である。
【図２】本開示の訓練データ抽出の処理手順を示す図である。
【図３】本開示の訓練データ抽出の処理手順を示す図である。
【図４】本開示の水位予測地点の上流域抽出の処理内容を示す図である。
【図５】本開示の上流域が広大な場合のデータ抽出の処理内容を示す図である。
【図６】本開示の水位ピーク期間抽出の処理内容を示す図である。
【図７】本開示の過去雨量期間設定の処理内容を示す図である。
【図８】本開示の雨量累加期間設定の処理内容を示す図である。
【図９】本開示の雨量累加期間設定の処理内容を示す図である。
【図１０】本開示の評価データ抽出の処理手順を示す図である。
【図１１】本開示の水位予測方法の予測結果を示す図である。
【発明を実施するための形態】
【００３２】
添付の図面を参照して本開示の実施形態を説明する。以下に説明する実施形態は本開示の実施の例であり、本開示は以下の実施形態に制限されるものではない。
【００３３】
（本開示の水位予測方法の処理手順）
本開示の水位予測方法の処理手順を図１に示す。本開示の水位予測方法をコンピュータに実行させるための水位予測プログラムを、コンピュータにインストールすることにより、本開示の水位予測方法を実行することができる。
【００３４】
まず、訓練時について説明する。訓練・評価データ抽出レイヤＥにおいて、雨量及び水位の訓練データＲｔ（ｒ、ｔ）、Ｌｔ（ｒ、ｔ）を入力する。そして、水位予測地点ｒｐ及びある時刻ｔｐの水位Ｌｔ（ｒｐ、ｔｐ）の学習に必要な、水位予測地点ｒｐ及び当該時刻ｔｐの水位Ｌｔ（ｒｐ、ｔｐ）との相関が高い雨量及び水位の訓練データＲｔ（ｒ、ｔ）、Ｌｔ（ｒ、ｔ）を抽出する。訓練データ抽出については、図２から図９までを用いて詳述する。
【００３５】
全結合入力レイヤＣ１、全結合隠れレイヤＣ２、活性関数レイヤＣ３及び全結合出力レイヤＣ４において、訓練・評価データ抽出レイヤＥで抽出された雨量及び水位の訓練データＲｔ（ｒ、ｔ）、Ｌｔ（ｒ、ｔ）に対して、全結合レイヤの処理を実行する。ここで、訓練・評価データ抽出レイヤＥで抽出された雨量及び水位の訓練データＲｔ（ｒ、ｔ）、Ｌｔ（ｒ、ｔ）は、時刻同期をとったうえで、全結合入力レイヤＣ１に入力される。そして、水位予測地点ｒｐ及び当該時刻ｔｐの水位Ｌｔ（ｒｐ、ｔｐ）を予測し、全結合入力レイヤＣ１及び全結合隠れレイヤＣ２の重み付け係数に対して、予測誤差を最適化するように誤差逆伝搬を実行する。さらに、以上に記載の処理が繰り返される。
【００３６】
このように、雨量及び水位の訓練データを適切に抽出することにより、雨量及び水位の訓練データを適切に学習し、訓練時の計算量を低減することができる。
【００３７】
次に、評価時について説明する。訓練・評価データ抽出レイヤＥにおいて、雨量及び水位の評価データＲｅ（ｒ、ｔ）、Ｌｅ（ｒ、ｔ）を入力する。そして、水位予測地点ｒｐ及びある時刻ｔｐの水位Ｌｐ（ｒｐ、ｔｐ）の予測に必要な、水位予測地点ｒｐ及び当該時刻ｔｐの水位Ｌｐ（ｒｐ、ｔｐ）との相関が高い雨量及び水位の評価データＲｅ（ｒ、ｔ）、Ｌｅ（ｒ、ｔ）を抽出する。評価データ抽出については、図１０を用いて詳述する。
【００３８】
全結合入力レイヤＣ１、全結合隠れレイヤＣ２、活性関数レイヤＣ３及び全結合出力レイヤＣ４において、訓練・評価データ抽出レイヤＥで抽出された雨量及び水位の評価データＲｅ（ｒ、ｔ）、Ｌｅ（ｒ、ｔ）に対して、全結合レイヤの処理を実行する。ここで、訓練・評価データ抽出レイヤＥで抽出された雨量及び水位の評価データＲｅ（ｒ、ｔ）、Ｌｅ（ｒ、ｔ）は、時刻同期をとったうえで、全結合入力レイヤＣ１に入力される。そして、水位予測地点ｒｐ及び当該時刻ｔｐの水位Ｌｐ（ｒｐ、ｔｐ）を予測する。
【００３９】
このように、雨量及び水位の評価データを適切に抽出することにより、水位を十分に精度高く予測し、評価時の計算量を低減することができる。
【００４０】
（本開示の訓練データ抽出の処理手順）
本開示の訓練データ抽出の処理手順を図２及び図３に示す。ステップＥ１及びステップＥ２については、いずれか一方のステップの処理を実行すればよい。
【００４１】
ステップＥ１では、水位予測地点ｒｐの水位との相関が高い水位の訓練データＬｔ（ｒ、ｔ）として、水位予測地点ｒｐの水位の訓練データを抽出する。そして、水位予測地点ｒｐの水位との相関が高い雨量の訓練データＲｔ（ｒ、ｔ）として、水位予測地点ｒｐの周辺の地形に基づいて、水位予測地点ｒｐの上流域ｒｐｕの雨量の訓練データを抽出する。
【００４２】
本開示の水位予測地点の上流域抽出の処理内容を図４に示す。水位予測地点ｒｐの水位との相関が高い雨量の訓練データを抽出する際に、水位予測地点ｒｐの上流域ｒｐｕは、水位予測地点ｒｐの周辺の地形に基づいて設定される。例えば、地上での測量情報又は衛星による地形情報に基づいて、水位予測地点ｒｐの周辺の分水境界を探索する。そして、水位予測地点ｒｐの周辺の分水境界に基づいて、水位予測地点ｒｐの上流域ｒｐｕを設定する。
【００４３】
具体的には、ステップＥで入力した水位の訓練データＬｔ（ｒ、ｔ）のうち、ｒ＝ｒｐを満たす水位の訓練データを抽出する。そして、ステップＥで入力した雨量の訓練データＲｔ（ｒ、ｔ）のうち、ｒ＝ｒｐｕを満たす雨量の訓練データを抽出する。
【００４４】
このように、上流域ｒｐｕより他の領域の雨量の訓練データを、学習に利用しない（データをゼロクリア等）ことにより、訓練時の計算量を低減することができる。
【００４５】
ステップＥ２では、水位予測地点ｒｐの水位との相関が高い水位の訓練データＬｔ（ｒ、ｔ）として、水位予測地点ｒｐの水位の訓練データと、水位予測地点ｒｐより上流地点ｒｕの水位の訓練データと、を抽出する。そして、水位予測地点ｒｐの水位との相関が高い雨量の訓練データＲｔ（ｒ、ｔ）として、水位予測地点ｒｐの周辺の地形及び上流地点ｒｕの周辺の地形に基づいて、上流地点ｒｕの上流域ｒｕｕを除く水位予測地点ｒｐの上流域ｒｐｕの雨量の訓練データを抽出する。
【００４６】
本開示の上流域が広大な場合のデータ抽出の処理内容を図５に示す。水位予測地点ｒｐの水位との相関が高い雨量の訓練データを抽出する際に、上流地点ｒｕの上流域ｒｕｕは、上流地点ｒｕの周辺の地形に基づいて設定される。例えば、地上での測量情報又は衛星による地形情報に基づいて、上流地点ｒｕの周辺の分水境界を探索する。そして、上流地点ｒｕの周辺の分水境界に基づいて、上流地点ｒｕの上流域ｒｕｕを設定する。なお、水位予測地点ｒｐの上流域ｒｐｕは、図４で示したように設定される。
【００４７】
具体的には、ステップＥで入力した水位の訓練データＬｔ（ｒ、ｔ）のうち、ｒ＝ｒｐ、ｒｕを満たす水位の訓練データを抽出する。そして、ステップＥで入力した雨量の訓練データＲｔ（ｒ、ｔ）のうち、ｒ＝ｒｐｕ、ｒ≠ｒｕｕを満たす雨量の訓練データを抽出する。
【００４８】
このように、上流地点ｒｕの上流域ｒｕｕの雨量の訓練データを、上流地点ｒｕの水位の訓練データに置換することにより、訓練時の計算量を低減することができる。
【００４９】
ステップＥ３では、ステップＥ１又はＥ２で抽出した水位の訓練データＬｔ（ｒ、ｔ）のうち、学習対象として水位ピーク期間の水位の訓練データを抽出する。
【００５０】
ここで、水位ピーク判別閾値を大きく設定しすぎると、十分なデータ量の水位変動期間（洪水時等）を学習することができず、一方で、水位ピーク判別閾値を小さく設定しすぎると、水位不変動期間（平常時等）を抽出するため予測精度を劣化させてしまう。そして、水位ピーク抽出幅を大きく設定しすぎると、水位不変動期間（平常時等）を抽出するため予測精度を劣化させてしまい、一方で、水位ピーク抽出幅を小さく設定しすぎると、上昇開始から下降終了までの水位変動期間（洪水時等）を学習することができない。
【００５１】
そこで、抽出される水位の訓練データの個数が削減されすぎず、水位不変動期間（平常時等）の水位の訓練データが抽出されないように、水位ピーク判別閾値を設定する。そして、水位ピーク期間の全体が抽出され、水位不変動期間（平常時等）の水位の訓練データが抽出されないように、水位ピーク時刻のまわりの水位ピーク抽出幅を設定する。
【００５２】
本開示の水位ピーク期間抽出の処理内容を図６に示す。図６では、水位ピーク判別閾値は、約１．５ｍに設定されているため、抽出される水位の訓練データの個数が削減されすぎず、水位ピーク時刻から離れた水位不変動期間（平常時等）の水位の訓練データが抽出されない。そして、水位ピーク抽出幅は、水位ピーク時刻の前後の約１日又は約２日に設定されているため、水位ピーク期間の全体が抽出され、水位ピーク時刻から離れた水位不変動期間（平常時等）の水位の訓練データが抽出されない。
【００５３】
ここで、水位ピーク抽出幅は、水位ピーク時刻の前後で、対称であってもよく非対称であってもよい。そして、水位ピーク抽出幅は、過去水位期間を満たす水位の訓練データ（ステップＥ８及びステップＥ９で詳述）を抽出可能な程度に十分に広いことが望ましい。
【００５４】
このように、水位ピーク判別閾値を適切に設定することにより、十分なデータ量の水位変動期間（洪水時等）のみを学習するため、水位の訓練データを適切に学習し、訓練時の計算量を低減することができる。そして、水位ピーク抽出幅を適切に設定することにより、上昇開始から下降終了までの水位変動期間（洪水時等）のみを学習するため、水位の訓練データを適切に学習し、訓練時の計算量を低減することができる。
【００５５】
ステップＥ４及びステップＥ５では、ステップＥ１又はＥ２で抽出した雨量の訓練データＲｔ（ｒ、ｔ）のうち、ある時刻ｔｐの水位との相関が高い雨量の訓練データとして、当該時刻ｔｐより以前の所定の過去雨量期間の雨量の訓練データを抽出する。
【００５６】
ここで、学習すべき時刻ｔｐの水位と相関が高い期間は、降雨が降雨地点から水位予測地点ｒｐへと到達するために必要である時間に等しい。そして、降雨が降雨地点から水位予測地点ｒｐへと到達するために必要である時間は、雨量ピーク時刻から水位ピーク時刻までに等しく、降雨地点から水位予測地点ｒｐまでの距離及び傾斜に応じて決まる。
【００５７】
そこで、ステップＥ４では、ステップＥ１又はＥ２で抽出した雨量の訓練データＲｔ（ｒ、ｔ）及びステップＥ３で抽出した水位の訓練データＬｔ（ｒ、ｔ）について、Ｒｔ（ｒ、ｔ）とＬｔ（ｒ、ｔ）との相関が最適化されるような、Ｒｔ（ｒ、ｔ）とＬｔ（ｒ、ｔ）との間のシフト時間ｔｓを、上記の過去雨量期間として設定する。そして、ステップＥ５では、ステップＥ１又はＥ２で抽出した雨量の訓練データＲｔ（ｒ、ｔ）のうち、過去雨量期間（ｔｐ−ｔｓ≦ｔ≦ｔｐ）を満たすデータを抽出する。
【００５８】
ここで、Ｒｔ（ｒ、ｔ）とＬｔ（ｒ、ｔ）との相関を計算するにあたり、Ｒｔ（ｒ、ｔ）については、上昇期間及び下降期間をともに含めることが望ましく、Ｌｔ（ｒ、ｔ）については、上昇期間のみ含めて下降期間を含めないことが望ましい。なぜならば、Ｒｔ（ｒ、ｔ）においては、ピークに向けて徐々に上昇し、ピークの後には急激に下降するが、Ｌｔ（ｒ、ｔ）においては、ピークに向けて急激に上昇し、ピークの後には徐々に下降する。そして、Ｌｔ（ｒ、ｔ）についても下降期間を含めれば、Ｒｔ（ｒ、ｔ）とＬｔ（ｒ、ｔ）との相関を計算しにくく、シフト時間ｔｓつまり過去雨量期間が大きくなりすぎる。
【００５９】
本開示の過去雨量期間設定の処理内容を図７に示す。図７の左欄では、Ｒｔ（ｒ、ｔ）とＬｔ（ｒ、ｔ）との相関についての、シフト時間ｔｓ依存性を示す。あるシフト時間ｔｓにおいて、Ｒｔ（ｒ、ｔ）とＬｔ（ｒ、ｔ）との相関が最適化されている。
【００６０】
Ｒｔ（ｒ、ｔ）とＬｔ（ｒ、ｔ）との相関が最適化されるシフト時間ｔｓは、図１１に示す多摩川、鈴鹿川、山国川、佐波川及び渡川において、それぞれ、３時間、５時間、５時間、６時間及び１１時間であり、河川毎に異なる値となった。
【００６１】
図７の右欄では、評価時の予測誤差（訓練時の予測誤差ではない）についての、過去雨量期間依存性を示す。実線は、予測誤差の標準誤差を示し、破線は、予測誤差の平均の絶対値を示し、一点鎖線は、両者の加算値を示す。全結合入力レイヤＣ１及び全結合隠れレイヤＣ２の重み付け係数の学習前初期値は、ランダムに設定されるため、予測誤差にばらつきがあるから、図７の右欄のように予測誤差の標準誤差や平均を考えたのである。
【００６２】
過去雨量期間をＲｔ（ｒ、ｔ）とＬｔ（ｒ、ｔ）との相関が最適化されるシフト時間ｔｓに設定したときに、図１１に示す多摩川、鈴鹿川、山国川、佐波川及び渡川において、評価時の予測誤差がほぼ最適化されることが、図７の左欄と右欄との比較で分かった。
【００６３】
このように、降雨地点から水位予測地点ｒｐまでの距離及び傾斜に応じて決まる、降雨が降雨地点から水位予測地点ｒｐへと到達するために必要である時間分だけ、つまり、学習すべき時刻ｔｐの水位と相関が高い期間分だけ、雨量の訓練データを抽出することにより、雨量の訓練データを適切に学習し、訓練時の計算量を低減することができる。
【００６４】
ステップＥ６及びステップＥ７では、ある時刻ｔｐの雨量の訓練データとして、当該時刻ｔｐより以前の所定の雨量累加期間の雨量の訓練データの累加結果を採用する。
【００６５】
ここで、降雨地点での土壌内又は貯水池等への貯水量に応じて、どのくらいの時間分だけ降雨が続けば降雨が土壌内又は貯水池等から流出するかが決まる。そして、短時間のみの降雨は、土壌内又は貯水池等に貯水されるため、水位予測地点ｒｐの水位に反映されにくいが、長時間に渡る降雨は、土壌内又は貯水池等から流出するため、水位予測地点ｒｐの水位に反映されやすい。そこで、ある時刻ｔｐの雨量の訓練データとして、当該時刻ｔｐの雨量の訓練データそのものを採用するのではなく、当該時刻ｔｐより以前の上記の雨量累加期間の雨量の訓練データの累加結果を採用するのである。
【００６６】
つまり、ステップＥ６では、ステップＥ５で抽出した雨量の訓練データＲｔ（ｒ、ｔ）について、数１のように変換を行う。ここで、ｔａは、雨量累加期間である。
【数１】
2019095240.tif 000003
そして、ステップＥ７では、ステップＥ６で変換した雨量の訓練データＲｔ（ｒ、ｔ）及びステップＥ３で抽出した水位の訓練データＬｔ（ｒ、ｔ）について、Ｒｔ（ｒ、ｔ＋ｔｓ）とＬｔ（ｒ、ｔ）との相関が飽和するような、雨量累加期間ｔａを設定する。ここで、ｔｓは、上記のシフト時間である。つまり、Ｒｔ（ｒ、ｔ＋ｔｓ）とＬｔ（ｒ、ｔ）との相関が飽和するまで、ステップＥ６及びステップＥ７を繰り返すのである。
【００６７】
本開示の雨量累加期間設定の処理内容を図８に示す。図８では、Ｒｔ（ｒ、ｔ＋ｔｓ）とＬｔ（ｒ、ｔ）との相関についての、雨量累加期間ｔａ依存性を示す。雨量累加期間ｔａが０から増加するにつれて、Ｒｔ（ｒ、ｔ＋ｔｓ）とＬｔ（ｒ、ｔ）との相関は急激に増加する。これは、短時間のみの降雨は、水位予測地点ｒｐの水位に反映されにくいが、長時間に渡る降雨は、水位予測地点ｒｐの水位に反映されやすいからである。
【００６８】
しかし、雨量累加期間ｔａがさらに増加するにつれて、Ｒｔ（ｒ、ｔ＋ｔｓ）とＬｔ（ｒ、ｔ）との相関は徐々にしか増加しなくなる。これは、Ｒｔ（ｒ、ｔ＋ｔｓ）において、時刻が近く相関が高い雨量の訓練データを次々と累加しているからである。
【００６９】
そして、Ｒｔ（ｒ、ｔ＋ｔｓ）において、時刻が近く相関が高い雨量の訓練データを次々と累加しても、図９に示すように水位の予測誤差を増加させるだけである。そこで、Ｒｔ（ｒ、ｔ＋ｔｓ）とＬｔ（ｒ、ｔ）との相関が飽和するような、雨量累加期間ｔａを設定するのである。ここで、「Ｒｔ（ｒ、ｔ＋ｔｓ）とＬｔ（ｒ、ｔ）との相関が飽和する」とは、図８の具体例では、雨量累加期間ｔａの増加に対する相関係数の増加割合が、雨量累加期間ｔａが０近傍であるときと比べて半減することを意味している。
【００７０】
本開示の雨量累加期間設定の処理内容を図９にも示す。図９の左欄では、Ｒｔ（ｒ、ｔ＋ｔｓ）とＬｔ（ｒ、ｔ）との相関についての、雨量累加期間ｔａ依存性を示す。ある雨量累加期間ｔａにおいて、Ｒｔ（ｒ、ｔ＋ｔｓ）とＬｔ（ｒ、ｔ）との相関が飽和している。
【００７１】
Ｒｔ（ｒ、ｔ＋ｔｓ）とＬｔ（ｒ、ｔ）との相関が飽和する雨量累加期間ｔａは、図１１に示す多摩川、鈴鹿川、山国川、佐波川及び渡川において、それぞれ、３時間、３時間、７時間、４時間及び１０時間であり、河川毎に異なる値となった。
【００７２】
図９の右欄では、評価時の予測誤差（訓練時の予測誤差ではない）についての、雨量累加期間ｔａ依存性を示す。実線は、予測誤差の標準誤差を示し、破線は、予測誤差の平均の絶対値を示し、一点鎖線は、両者の加算値を示す。全結合入力レイヤＣ１及び全結合隠れレイヤＣ２の重み付け係数の学習前初期値は、ランダムに設定されるため、予測誤差にばらつきがあるから、図９の右欄のように予測誤差の標準誤差や平均を考えたのである。
【００７３】
雨量累加期間ｔａをＲｔ（ｒ、ｔ＋ｔｓ）とＬｔ（ｒ、ｔ）との相関が飽和する相関飽和時間に設定したときに、図１１に示す多摩川、鈴鹿川、山国川、佐波川及び渡川において、評価時の予測誤差がほぼ最適化されることが、図９の左欄と右欄との比較で分かった。
【００７４】
このように、降雨地点での土壌内又は貯水池等への貯水量に応じて決まる、どのくらいの時間分だけ降雨が続けば降雨が土壌内又は貯水池等から流出するかに応じて、雨量の訓練データを適切に学習し、訓練時の計算量を低減することができる。
【００７５】
ステップＥ８及びステップＥ９では、ステップＥ３で抽出した水位の訓練データＬｔ（ｒ、ｔ）のうち、ある時刻ｔｐの水位との相関が高い水位の訓練データとして、当該時刻ｔｐより過去の所定の過去水位期間の水位の訓練データを抽出する。
【００７６】
ここで、学習すべき時刻ｔｐの水位と相関が高い期間は、水位が上昇又は下降の傾向にある時間に等しい。そして、水位が上昇又は下降の傾向にある時間は、降雨地点から水位予測地点ｒｐまでの距離及び傾斜に応じて決まる。
【００７７】
そこで、ステップＥ８では、ステップＥ３で抽出した水位の訓練データＬｔ（ｒ、ｔ）について、Ｌｔ（ｒ、ｔ）の上昇傾向又は下降傾向が持続する時間ｔｃを、上記の過去水位期間として設定する。例えば、図６によれば、期間ｔｃを約１日程度に設定すればよい。そして、ステップＥ９では、ステップＥ３で抽出した水位の訓練データＬｔ（ｒ、ｔ）のうち、過去水位期間（ｔｐ−ｔｃ≦ｔ＜ｔｐ）を満たすデータを抽出する。
【００７８】
このように、降雨地点から水位予測地点ｒｐまでの距離及び傾斜に応じて決まる、水位が上昇又は下降の傾向にある時間分だけ、つまり、学習すべき時刻ｔｐの水位と相関が高い期間分だけ、水位の訓練データを抽出することにより、水位の訓練データを適切に学習し、訓練時の計算量を低減することができる。
【００７９】
なお、雨量の訓練データＲｔ（ｒ、ｔ）及び水位の訓練データＬｔ（ｒ、ｔ）を、時空間内で適切に絞り込むことができるならば、ステップＥ１からＥ９までの順序は、図２及び図３に示したものに限定されない。ただし、ステップＥ６及びＥ７は、ステップＥ４及びＥ５を前提とするため、ステップＥ４及びＥ５の後に実行すればよい。
【００８０】
（本開示の評価データ抽出の処理手順）
本開示の評価データ抽出の処理手順を図１０に示す。ステップＥ１１及びステップＥ１２については、いずれか一方のステップの処理を実行すればよい。本開示の評価データ抽出の処理手順は、本開示の訓練データ抽出の処理手順と、ほぼ同様である。
【００８１】
ステップＥ１１では、水位予測地点ｒｐの水位との相関が高い水位の評価データＬｅ（ｒ、ｔ）として、水位予測地点ｒｐの水位の評価データを抽出する。そして、水位予測地点ｒｐの水位との相関が高い雨量の評価データＲｅ（ｒ、ｔ）として、水位予測地点ｒｐの周辺の地形に基づいて、水位予測地点ｒｐの上流域ｒｐｕの雨量の評価データを抽出する。
【００８２】
すると、上流域ｒｐｕより他の領域の雨量の評価データを、予測に利用しない（データをゼロクリア等）ことにより、評価時の計算量を低減することができる。
【００８３】
ステップＥ１２では、水位予測地点ｒｐの水位との相関が高い水位の評価データＬｅ（ｒ、ｔ）として、水位予測地点ｒｐの水位の評価データと、水位予測地点ｒｐより上流地点ｒｕの水位の評価データと、を抽出する。そして、水位予測地点ｒｐの水位との相関が高い雨量の評価データＲｅ（ｒ、ｔ）として、水位予測地点ｒｐの周辺の地形及び上流地点ｒｕの周辺の地形に基づいて、上流地点ｒｕの上流域ｒｕｕを除く水位予測地点ｒｐの上流域ｒｐｕの雨量の評価データを抽出する。
【００８４】
すると、上流地点ｒｕの上流域ｒｕｕの雨量の評価データを、上流地点ｒｕの水位の評価データに置換することにより、評価時の計算量を低減することができる。
【００８５】
ステップＥ１３では、ステップＥ１１又はＥ１２で抽出した雨量の評価データＲｅ（ｒ、ｔ）のうち、過去雨量期間（ｔｐ−ｔｓ≦ｔ≦ｔｐ）を満たすデータを抽出する。ここで、過去雨量期間は、ステップＥ５及びステップＥ１３において同様である。そして、過去雨量期間（ｔｐ−ｔｓ≦ｔ≦ｔｐ）を満たす雨量の評価データは、現在時刻までの雨量の実測データを含んでいるが、水位をさらに精度高く予測するためには、将来時刻での雨量の予測データを含んでもよい。
【００８６】
すると、降雨地点から水位予測地点ｒｐまでの距離及び傾斜に応じて決まる、降雨が降雨地点から水位予測地点ｒｐへと到達するために必要である時間分だけ、つまり、予測すべき時刻ｔｐの水位と相関が高い期間分だけ、雨量の評価データを抽出することにより、水位を十分に精度高く予測し、評価時の計算量を低減することができる。
【００８７】
ステップＥ１４では、ステップＥ１３で抽出した雨量の評価データＲｅ（ｒ、ｔ）について、数２のように変換を行う。ここで、ｔａは、雨量累加期間である。そして、雨量累加期間に渡る雨量の評価データは、現在時刻までの雨量の実測データを含んでいるが、水位をさらに精度高く予測するためには、将来時刻での雨量の予測データを含んでもよい。
【数２】
2019095240.tif 000004
【００８８】
すると、降雨地点での土壌内又は貯水池等への貯水量に応じて決まる、どのくらいの時間分だけ降雨が続けば降雨が土壌内又は貯水池等から流出するかに応じて、水位を十分に精度高く予測し、評価時の計算量を低減することができる。
【００８９】
ステップＥ１５では、ステップＥ１１又はステップＥ１２で抽出した水位の評価データＬｅ（ｒ、ｔ）のうち、過去水位期間（ｔｐ−ｔｃ≦ｔ＜ｔｐ）を満たすデータを抽出する。ここで、過去水位期間は、ステップＥ９及びステップＥ１５において同様である。そして、過去水位期間（ｔｐ−ｔｃ≦ｔ＜ｔｐ）を満たす水位の評価データは、現在時刻までの水位の実測データを含んでいるが、水位をさらに精度高く予測するためには、将来時刻での水位の予測データを含んでもよい。
【００９０】
すると、降雨地点から水位予測地点ｒｐまでの距離及び傾斜に応じて決まる、水位が上昇又は下降の傾向にある時間分だけ、つまり、予測すべき時刻ｔｐの水位と相関が高い期間分だけ、水位の評価データを抽出することにより、水位を十分に精度高く予測し、評価時の計算量を低減することができる。
【００９１】
なお、雨量の評価データＲｅ（ｒ、ｔ）及び水位の評価データＬｅ（ｒ、ｔ）を、時空間内で適切に絞り込むことができるならば、ステップＥ１１からＥ１５までの順序は、図１０に示したものに限定されない。ただし、ステップＥ１４は、ステップＥ１３を前提とするため、ステップＥ１３の後に実行すればよい。
【００９２】
（本開示の水位予測方法の予測結果）
本開示の水位予測方法の予測結果を図１１に示す。図１１では、以下の（１）〜（４）の視点に基づいて、多摩川、鈴鹿川、佐波川、渡川及び山国川を選定している。
（１）河川の特徴（長さ及び傾斜等）が異なる、（２）雨量データの欠損が少ない、（３）水位に反映されにくい積雪が少ない、（４）水位予測を困難にするダムがない又は少ない。
【００９３】
各河川についての訓練期間及び評価期間は、以下の通りである。
多摩川：訓練期間＝２００８〜２０１３、２０１５年、評価期間＝２０１４年
鈴鹿川：訓練期間＝２００８〜２０１５年、評価期間＝２０１６年
佐波川：訓練期間＝２００８〜２０１５年、評価期間＝２０１６年
渡川：訓練期間＝２００９〜２０１５年、評価期間＝２０１６年
山国川：訓練期間＝２００８〜２０１５年、評価期間＝２０１６年
【００９４】
図１１では、上記の評価期間のうち、水位が最も高くなった期間について、水位の予測及び実測の結果を示している。１時間毎に、１〜３時間後だけ、水位が予測されている。各河川によらず、増水及び減水の傾向がとらえられていることが分かる。
【産業上の利用可能性】
【００９５】
本開示の水位予測方法、水位予測プログラム及び水位予測装置は、水害の予報又は警報を発令するために、河川、ダム又は下水等の水位を予測することができる。
【符号の説明】
【００９６】
Ｅ：訓練・評価データ抽出レイヤ
Ｃ１：全結合入力レイヤ
Ｃ２：全結合隠れレイヤ
Ｃ３：活性関数レイヤ
Ｃ４：全結合出力レイヤ
【図１】
2019095240.tif 000005
【図２】
2019095240.tif 000006
【図３】
2019095240.tif 000007
【図４】
2019095240.tif 000008
【図５】
2019095240.tif 000009
【図６】
2019095240.tif 000010
【図７】
2019095240.tif 000011
【図８】
2019095240.tif 000012
【図９】
2019095240.tif 000013
【図１０】
2019095240.tif 000014
【図１１】
2019095240.tif 000015