(57)【要約】【課題】雨量から指定された地点の水位を適切に予測する学習モデルを構築する。【解決手段】水位予測地点を選択する予測地点選択部と、前記水位予測地点、地形情報、および雨量観測地点データを基に、前記水位予測地点の水位を予測する学習モデルの学習用の雨量観測値データを選択する雨量観測地点選択処理部と、を有することを特徴とする。【選択図】図１４
2020002658.tif 000002
【特許請求の範囲】
【請求項１】
水位予測地点を選択する予測地点選択部と、
前記水位予測地点、地形情報、および雨量観測地点データを基に、前記水位予測地点の水位を予測する学習モデルの学習用の雨量観測値データを選択する雨量観測地点選択処理部と、
を有することを特徴とする情報処理装置。
【請求項２】
前記雨量観測地点選択処理部は、
前記水位予測地点より標高が高い前記雨量観測地点データであって、かつ、前記水位予測地点の位置情報と雨量観測地点データの位置情報が接続する線分が前記地形情報に含まれる分水嶺情報と交差しない前記雨量観測地点データの雨量観測値データを選択することを特徴とする請求項１に記載の情報処理装置。
【請求項３】
選択された前記雨量観測値データと、前記水位予測地点の水位観測値データを基に、前記水位予測地点の水位を予測する水位予測モデルを学習するモデル学習処理部とを有することを特徴とする請求項１および２記載の情報処理装置。
【請求項４】
前記水位予測モデルに、前記雨量観測値データを入力して、前記水位予測地点の水位を予測する水位予測処理部を有する請求項３記載の情報処理装置。
【請求項５】
水位予測地点を選択し、
前記水位予測地点、地形情報、および雨量観測地点データを基に、前記水位予測地点の水位を予測する学習モデルの学習用の雨量観測値データを選択し、
選択された前記雨量観測値データと、前記水位予測地点の水位観測値データを基に、前記水位予測地点の水位を予測する水位予測モデルを学習し、
前記水位予測モデルに、前記雨量観測値データを入力して、前記水位予測地点の水位を予測する
処理をコンピュータが実行することを特徴とする河川水位予測方法。
【請求項６】
水位予測地点を選択し、
前記水位予測地点、地形情報、および雨量観測地点データを基に、前記水位予測地点の水位を予測する学習モデルの学習用の雨量観測値データを選択し、
選択された前記雨量観測値データと、前記水位予測地点の水位観測値データを基に、前記水位予測地点の水位を予測する水位予測モデルを学習し、
前記水位予測モデルに、前記雨量観測値データを入力して、前記水位予測地点の水位を予測する
処理をコンピュータに実行させることを特徴とする河川水位予測プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、情報処理装置、河川水位予測方法及び河川水位予測プログラムに関する。
【背景技術】
【０００２】
河川の水位を予測する技術が従来から知られている。例えば、機械学習の手法により予め学習したモデル（機械学習モデル）を用いて、予測される雨量と、過去の雨量と、過去の水位とから、当該河川の水位を予測する技術が知られている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００７−２３６０８号公報
【特許文献２】特開２００７−６３８４９号公報
【特許文献３】特開２００７−２０５００１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ここで、従来の機械学習モデルを用いた水位の予測では、水位を予測する時点よりも前の過去の水位を示すデータが必要であった。このため、例えば、水位観測地点に設置されている水位観測機器の不具合等によって過去の水位を示すデータに欠損があった場合には、水位の予測を行うことができなかった。
【０００５】
他方で、過去の水位を用いずに機械学習モデルにより水位を予測する場合、多数の雨量観測地点で予測及び観測される雨量を用いて、当該機械学習モデルを学習させる必要がある。しかしながら、例えば、水位を予測する地点と関連性が低い雨量観測地点で予測及び観測される雨量を用いた場合、水位の予測精度が低下することがある。
【０００６】
開示の技術は、雨量から指定された地点の水位を適切に予測する学習モデルを構築することを目的とする。
【課題を解決するための手段】
【０００７】
開示の技術は、水位予測地点を選択する予測地点選択部と、前記水位予測地点、地形情報、および雨量観測地点データを基に、前記水位予測地点の水位を予測する学習モデルの学習用の雨量観測値データを選択する雨量観測地点選択処理部と、を有することを特徴とする。
【発明の効果】
【０００８】
雨量から指定された地点の水位を適切に予測する学習モデルを構築することができる。
【図面の簡単な説明】
【０００９】
【図１】本実施形態に係る水位予測装置の全体構成の一例を示す図である。
【図２】本実施形態に係る水位予測装置のハードウェア構成の一例を示す図である。
【図３】水位観測地点データの一例を示す図である。
【図４】標高情報データの一例を示す図である。
【図５】分水嶺情報データの一例を示す図である。
【図６】雨量観測地点データの一例を示す図である。
【図７】雨量観測値データの一例を示す図である。
【図８】水位観測値データの一例を示す図である。
【図９】本実施形態に係る雨量観測地点選択処理部の機能構成の一例を示す図である。
【図１０】本実施形態に係るモデル学習処理部の機能構成の一例を示す図である。
【図１１】本実施形態に係る水位予測処理部の機能構成の一例を示す図である。
【図１２】本実施形態に係る雨量観測地点の選択処理の一例を示すフローチャートである。
【図１３】流域の特定の一例を説明するための図である。
【図１４】雨量観測地点の選択の一例を説明するための図である。
【図１５】計算時間と予測精度との関係の一例を説明するための図である。
【図１６】本実施形態に係る水位予測モデルの学習処理の一例を示すフローチャートである。
【図１７】本実施形態に係る河川水位の予測処理の一例を示すフローチャートである。
【図１８】予測精度の一例を説明するための図である。
【発明を実施するための形態】
【００１０】
以下、本発明の実施の形態（以降、「本実施形態」とも表す。）について、図面を参照しながら説明する。本実施形態では、機械学習の手法により、雨量（予測される雨量及び過去の雨量）から、将来の或る時点（以降、「予測時点」とも表す。）の水位を予測する水位予測装置１０について説明する。本実施形態に係る水位予測装置１０には、機械学習の手法によりモデルを学習する「モデル学習時」と、学習済みのモデルを用いて雨量から水位を予測する「水位予測時」とがある。以降では、水位の予測に用いられるモデルを「水位予測モデル」とも表す。
【００１１】
まず、本実施形態に係る水位予測装置１０の全体構成について、図１を参照しながら説明する。図１は、本実施形態に係る水位予測装置１０の全体構成の一例を示す図である。
【００１２】
図１に示すように、本実施形態に係る水位予測装置１０は、雨量観測地点選択処理部１１０と、モデル学習処理部１２０と、水位予測処理部１３０とを有する。
【００１３】
また、本実施形態に係る水位予測装置１０は、水位観測地点記憶部２１０と、標高情報記憶部２２０と、分水嶺情報記憶部２３０と、雨量観測地点記憶部２４０と、雨量観測値記憶部２５０と、水位観測値記憶部２６０とを有する。
【００１４】
水位観測地点記憶部２１０は、河川の水位を観測する水位観測地点に関するデータ（以降、「水位観測地点データ」とも表す。）を記憶する。なお、水位観測地点とは、例えば、河川の水位を測るための水位観測機器が設置等されている観測所が位置する地点のことである。
【００１５】
標高情報記憶部２２０は、地図上に予め設定された各地点の標高に関するデータ（以降、「標高情報データ」とも表す。）を記憶する。
【００１６】
分水嶺情報記憶部２３０は、地図上で分水嶺を表す地点に関するデータ（以降、「分水嶺情報データ」とも表す。）を記憶する。
【００１７】
雨量観測地点記憶部２４０は、雨量を観測する雨量観測地点に関するデータ（以降、「雨量観測地点データ」とも表す。）を記憶する。なお、雨量観測地点とは、例えば、雨量を測るための雨量観測機器が設置等されている観測所が位置する地点のことである。
【００１８】
雨量観測値記憶部２５０は、各雨量観測地点でそれぞれ観測された雨量観測値に関するデータ（以降、「雨量観測値データ」とも表す。）を記憶する。
【００１９】
水位観測値記憶部２６０は、各水位観測地点でそれぞれ観測された水位観測値に関するデータ（以降、「水位観測値データ」とも表す。）を記憶する。
【００２０】
なお、水位観測地点データ、標高情報データ、分水嶺情報データ、雨量観測地点データ、雨量観測値データ及び水位観測値データの詳細については後述する。
【００２１】
雨量観測地点選択処理部１１０は、モデル学習時において、水位を予測する水位観測地点（以降、「水位予測地点」とも表す。）を選択した上で、水位予測地点の流域を特定する。そして、雨量観測地点選択処理部１１０は、特定した流域との位置関係に基づいて、複数の雨量観測地点のうち、水位予測地点と関連性が高い雨量観測地点を選択する処理を行う。このとき、雨量観測地点選択処理部１１０は、水位観測地点記憶部２１０と、標高情報記憶部２２０と、分水嶺情報記憶部２３０とを参照することで、水位予測地点の選択と、流域の特定とを行う。また、雨量観測地点選択処理部１１０は、雨量観測地点記憶部２４０を参照することで、雨量観測地点の選択を行う。なお、水位予測地点の流域とは、降雨があった場合に、水位予測地点に雨水が流入する領域のことである。
【００２２】
ここで、水位予測地点と関連性が高い雨量観測地点とは、当該雨量観測地点で観測される雨量の観測値と、水位予測地点における水位の観測値との間に相関関係（例えば、相関係数が所定の値以上の強さの相関関係）がある雨量観測地点のことである。
【００２３】
以降では、雨量観測地点選択処理部１１０により選択された雨量観測地点を「選択雨量観測地点」とも表す。
【００２４】
モデル学習処理部１２０は、モデル学習時において、選択雨量観測地点の雨量観測値と、水位予測地点の水位観測値とを用いて、水位予測モデルを学習する処理を行う。このとき、モデル学習処理部１２０は、雨量観測値記憶部２５０と、水位観測値記憶部２６０とを参照することで、水位予測モデルの学習を行う。なお、水位予測モデルの学習は、既知の最適化手法により当該水位予測モデルのパラメータが更新されることで行われる。
【００２５】
水位予測処理部１３０は、水位予測時において、学習済みの水位予測モデルを用いて、選択雨量観測地点の雨量観測値及び雨量予測値から水位予測値を算出する。ここで、水位予測値とは、水位予測地点の予測時点における水位を示す値のことである。また、雨量予測値とは、将来の各日時における選択雨量観測地点の雨量の予測値である。なお、雨量予測値としては、例えば、気象庁等により提供される雨量の予測値が用いられる。
【００２６】
以上のように、本実施形態に係る水位予測装置１０は、水位予測地点の流域との位置関係に基づいて、複数の雨量観測地点のうち、当該水位予測地点と関連性が高い雨量観測地点を選択する。そして、本実施形態に係る水位予測装置１０は、選択した雨量観測地点の雨量観測値と、水位予測地点の水位観測値とを用いて、水位予測モデルを学習する。このように、水位観測地点と関連性の高い雨量観測地点の雨量観測値を水位予測モデルの学習に用いることで、水位予測時における予測精度を向上させると共に、計算時間を削減することができる。
【００２７】
したがって、本実施形態に係る水位予測装置１０は、水位予測時において、雨量観測値及び雨量予測値から水位予測値を高い精度で、かつ、迅速に予測することができるようになる。
【００２８】
次に、本実施形態に係る水位予測装置１０のハードウェア構成について、図２を参照しながら説明する。図２は、本実施形態に係る水位予測装置１０のハードウェア構成の一例を示す図である。
【００２９】
図２に示すように、本実施形態に係る水位予測装置１０は、入力装置１１と、表示装置１２と、外部Ｉ／Ｆ１３と、通信Ｉ／Ｆ１４と、ＲＯＭ（Read  Only  Memory）１５と、ＲＡＭ（Random  Access  Memory）１６と、プロセッサ１７と、補助記憶装置１８とを有する。これら各ハードウェアは、バスＢを介して通信可能に接続されている。
【００３０】
入力装置１１は、例えばキーボードやマウス、タッチパネル等を含み、ユーザが各種操作を行うのに用いられる。表示装置１２は、水位予測装置１０の処理結果（例えば、水位予測値等）を表示する。なお、水位予測装置１０は、入力装置１１及び表示装置１２のうちの少なくとも一方を有していなくても良い。
【００３１】
外部Ｉ／Ｆ１３は、外部装置とのインタフェースである。外部装置には、記録媒体１３ａ等がある。水位予測装置１０は、外部Ｉ／Ｆ１３を介して、記録媒体１３ａの読み取りや書き込み等を行うことができる。記録媒体１３ａには、例えば、ＳＤメモリカード（SD  memory  card）やＵＳＢ（Universal  Serial  Bus）メモリ、ＣＤ（Compact  Disk）、ＤＶＤ（Digital  Versatile  Disk）等がある。
【００３２】
通信Ｉ／Ｆ１４は、水位予測装置１０が通信ネットワークに接続するためのインタフェースである。ＲＯＭ１５は、電源を切ってもデータを保持することができる不揮発性の半導体メモリである。ＲＡＭ１６は、プログラムやデータを一時保持する揮発性の半導体メモリである。
【００３３】
プロセッサ１７は、例えばＣＰＵ（Central  Processing  Unit）やＧＰＵ（Graphics  Processing  Unit）等の演算装置である。プロセッサ１７は、ＲＯＭ１５や補助記憶装置１８等の記憶装置からプログラムやデータをＲＡＭ１６上に読み出し、当該プログラムやデータに基づく処理を実行することで、水位予測装置１０全体の制御や機能を実現する。なお、水位予測装置１０が有する各処理部（雨量観測地点選択処理部１１０、モデル学習処理部１２０及び水位予測処理部１３０）は、例えば、当該水位予測装置１０にインストールされた１以上のプログラムがプロセッサ１７に実行させる処理により実現される。
【００３４】
補助記憶装置１８は、例えばＨＤＤ（Hard  Disk  Drive）やＳＳＤ（Solid  State  Drive）等であり、プログラムやデータを格納している不揮発性の記憶装置である。補助記憶装置１８には、例えば、水位予測装置１０全体を制御する基本ソフトウェアであるＯＳ（Operating  System）、ＯＳ上で各種機能を提供するアプリケーションソフトウェア等がある。また、補助記憶装置１８には、雨量観測地点選択処理部１１０やモデル学習処理部１２０、水位予測処理部１３０等を実現する１以上のプログラム等が格納されている。なお、水位予測装置１０が有する各記憶部（例えば、水位観測地点記憶部２１０、標高情報記憶部２２０、分水嶺情報記憶部２３０、雨量観測地点記憶部２４０、雨量観測値記憶部２５０及び水位観測値記憶部２６０）は、補助記憶装置１８を用いて実現可能である。
【００３５】
本実施形態に係る水位予測装置１０は、図２に示すハードウェア構成を有することにより、本実施形態に係る各種処理を実現することができる。
【００３６】
次に、水位観測地点記憶部２１０に記憶されている水位観測地点データについて、図３を参照しながら説明する。図３は、水位観測地点データの一例を示す図である。
【００３７】
図３に示すように、水位観測地点記憶部２１０には、１以上の水位観測地点データが記憶されている。また、各水位観測地点データには、データ項目として、水位観測地点名称と、緯度と、経度とが含まれる。
【００３８】
データ項目「水位観測地点名称」には、水位観測地点の名称（例えば、水位観測地点を識別する識別情報等）が設定される。データ項目「緯度」には、水位観測地点の緯度が設定される。データ項目「経度」には、水位観測地点の経度が設定される。例えば、水位観測地点名称が「水位観測地点１」の水位観測地点データには、緯度「Ｘ１１」と経度「Ｙ１１」とが含まれている。
【００３９】
このように、水位観測地点記憶部２１０には、水位観測地点の位置座標（すなわち、緯度及び経度の組）が含まれる水位観測地点データが記憶されている。以降では、水位観測地点の総数をＭとして、各水位観測地点をｍで表すものとする。すなわち、水位観測地点記憶部２１０には、Ｍ個の水位観測地点データが記憶されているものとする。
【００４０】
次に、標高情報記憶部２２０に記憶されている標高情報データについて、図４を参照しながら説明する。図４は、標高情報データの一例を示す図である。
【００４１】
図４に示すように、標高情報記憶部２２０には、１以上の標高情報データが記憶されている。また、各標高情報データには、データ項目として、地点名称と、緯度と、経度と、標高とが含まれる。
【００４２】
データ項目「地点名称」には、地図上に予め設定された地点の名称（例えば、地点を識別する識別情報等）が設定される。なお、このような地点は、例えば、所定の間隔毎に地図上に設定された格子点等である。
【００４３】
データ項目「緯度」には、当該地点の緯度が設定される。データ項目「経度」には、当該地点の経度が設定される。データ項目「標高」には、当該地点の標高が設定される。例えば、地点名称が「地点１」の標高情報データには、緯度「Ｘ２１」と、経度「Ｙ２１」と、標高「Ｚ２１」とが含まれている。
【００４４】
このように、標高情報記憶部２２０に記憶されている標高情報データには、地図上の各地点の位置座標と標高とが含まれる標高情報データが記憶されている。
【００４５】
次に、分水嶺情報記憶部２３０に記憶されている分水嶺情報データについて、図４を参照しながら説明する。図５は、分水嶺情報データの一例を示す図である。
【００４６】
図５に示すように、分水嶺情報記憶部２３０には、１以上の分水嶺情報データが記憶されている。また、各分水嶺情報データには、データ項目として、分水嶺名称と、接続関係とが含まれる。
【００４７】
データ項目「分水嶺名称」には、分水嶺の名称（例えば、分水嶺を識別する識別情報等）が設定される。データ項目「接続関係」には、分水嶺に含まれる地点の接続関係を表す１以上の組が設定される。また、各組は、当該組に含まれる地点同士が線分で接続されることを表す。これにより、各組がそれぞれ表す線分により分水嶺が表される。
【００４８】
例えば、分水嶺名称が「分水嶺１」の分水嶺情報データには、接続関係「（地点１，地点４），（地点４，地点５），・・・，（地点８，地点１１）」が含まれる。これは、分水嶺１は、地点１と地点４とを接続する線分と、地点４と地点５とを接続する線分と、・・・、地点８と地点１１とを接続する線分とで表されることを示している。
【００４９】
このように、分水嶺情報記憶部２３０に記憶されている分水嶺情報データには、分水嶺を表す１以上の線分の接続関係が含まれる分水嶺情報データが記憶されている。
【００５０】
次に、雨量観測地点記憶部２４０に記憶されている雨量観測地点データについて、図６を参照しながら説明する。図６は、雨量観測地点データの一例を示す図である。
【００５１】
図６に示すように、雨量観測地点記憶部２４０には、１以上の雨量観測地点データが記憶されている。また、各雨量観測地点データには、データ項目として、雨量観測地点名称と、緯度と、経度とが含まれる。
【００５２】
データ項目「雨量観測地点名称」には、雨量観測地点の名称（例えば、雨量観測地点を識別する識別情報等）が設定される。データ項目「緯度」には、雨量観測地点の緯度が設定される。データ項目「経度」には、雨量観測地点の経度が設定される。例えば、雨量観測地点名称が「雨量観測地点１」の雨量観測地点データには、緯度「Ｘ３１」と経度「Ｙ３１」とが含まれている。
【００５３】
このように、雨量観測地点記憶部２４０には、雨量観測地点の位置座標が含まれる雨量観測地点データが記憶されている。以降では、雨量観測地点の総数をＩとして、各雨量観測地点をｉで表すものとする。すなわち、雨量観測地点記憶部２４０には、Ｉ個の雨量観測地点データが記憶されているものとする。
【００５４】
次に、雨量観測値記憶部２５０に記憶されている雨量観測値データについて、図７を参照しながら説明する。図７は、雨量観測値データの一例を示す図である。
【００５５】
図７に示すように、雨量観測値記憶部２５０には、雨量観測地点毎に、１以上の雨量観測値データが記憶されている。また、各雨量観測値データには、データ項目として、雨量観測地点名称と、日時と、雨量観測値とが含まれる。
【００５６】
データ項目「雨量観測地点名称」には、雨量観測地点の名称が設定される。データ項目「日時」には、当該雨量観測地点で雨量が観測された日時（例えば、年月日及び時刻）が設定される。データ項目「雨量観測値」には、当該雨量観測地点で観測された雨量の観測値が設定される。
【００５７】
例えば、雨量観測地点名称が「雨量観測地点１」の日時「２０１８／４／１  １０：１０」における雨量観測値データには、雨量観測値「０．５」が含まれる。同様に、例えば、雨量観測地点名称が「雨量観測地点１」の日時「２０１８／４／１  １０：２０」における雨量観測値データには、雨量観測値「０．５」が含まれる。なお、雨量観測値の単位は、例えば、「ｍｍ」である。
【００５８】
このように、雨量観測値記憶部２５０には、雨量観測地点毎に、当該雨量観測地点で観測された雨量観測値の時系列データである雨量観測値データが記憶されている。なお、雨量観測値記憶部２５０には、所定の時間毎（例えば、１０分毎や１時間毎等）に、例えば、数か月分、数年分又は数十年分の雨量観測値データが記憶されているものとする。
【００５９】
次に、水位観測値記憶部２６０に記憶されている水位観測値データについて、図８を参照しながら説明する。図８は、水位観測値データの一例を示す図である。
【００６０】
図８に示すように、水位観測値記憶部２６０には、水位観測地点毎に、１以上の水位観測値データが記憶されている。また、各水位観測値データには、データ項目として、水位観測地点名称と、日時と、水位観測値とが含まれる。
【００６１】
データ項目「水位観測地点名称」には、水位観測地点の名称が設定される。データ項目「日時」には、当該水位観測地点で水位が観測された日時が設定される。データ項目「水位観測値」には、当該水位観測地点で観測された水位の観測値が設定される。
【００６２】
例えば、水位観測地点名称が「水位観測地点１」の日時「２０１８／４／１  １０：１０」における水位観測値データには、水位観測値「２．６８」が含まれる。同様に、例えば、水位観測地点名称が「水位観測地点１」の日時「２０１８／４／１  １０：２０」の水位観測値データには、水位観測値「２．６９」が含まれる。なお、水位観測値の単位は「ｍ」である。
【００６３】
このように、水位観測値記憶部２６０には、水位観測地点毎に、当該水位観測地点で観測された水位観測値の時系列データである水位観測値データが記憶されている。なお、水位観測値記憶部２６０には、所定の時間毎（例えば、１０分毎や１時間毎等）に、例えば、数か月分、数年分又は数十年分の水位観測値データが記憶されているものとする。
【００６４】
次に、本実施形態に係る雨量観測地点選択処理部１１０の機能構成について、図９を参照しながら説明する。図９は、本実施形態に係る雨量観測地点選択処理部１１０の機能構成の一例を示す図である。
【００６５】
図９に示すように、本実施形態に係る雨量観測地点選択処理部１１０には、予測地点選択部１１１と、流域特定部１１２と、雨量観測地点選択部１１３とが含まれる。
【００６６】
予測地点選択部１１１は、水位観測地点記憶部２１０に記憶されているＭ個の水位観測地点データがそれぞれ示すＭ個の水位観測地点のうち、水位を予測する水位観測地点（すなわち、水位予測地点）を選択する。
【００６７】
流域特定部１１２は、標高情報記憶部２２０に記憶されている標高情報データと、分水嶺情報記憶部２３０に記憶されている分水嶺情報データとを用いて、水位予測地点の流域を特定する。流域特定部１１２は、例えば、水位予測地点よりも標高が高い地点で、かつ、当該地点と水位予測地点とを接続する線分が分水嶺と交差しない地点の集合を流域と特定する。
【００６８】
雨量観測地点選択部１１３は、特定された流域との位置関係に基づいて、雨量観測地点記憶部２４０に記憶されているＩ個の雨量観測地点データがそれぞれ示すＩ個の雨量観測地点のうち、水位予測地点と関連性が高い雨量観測地点（選択雨量観測地点）を選択する。
【００６９】
雨量観測地点選択部１１３は、例えば、水位予測地点の流域内に含まれる雨量観測地点（第１の雨量観測地点）と、当該流域内に含まれない雨量観測地点のうち、所定の条件を満たす雨量観測地点（第２の雨量観測地点）とを選択雨量観測地点に選択する。
【００７０】
次に、本実施形態に係るモデル学習処理部１２０の機能構成について、図１０を参照しながら説明する。図１０は、本実施形態に係るモデル学習処理部１２０の機能構成の一例を示す図である。
【００７１】
図１０に示すように、本実施形態に係るモデル学習処理部１２０には、取得部１２１と、学習データ作成部１２２と、出力値算出部１２３と、パラメータ更新部１２４とが含まれる。
【００７２】
取得部１２１は、雨量観測値記憶部２５０に記憶されている雨量観測値データのうち、選択雨量観測地点の雨量観測値データを取得する。また、取得部１２１は、水位観測値記憶部２６０に記憶されている水位観測値データのうち、予測水位地点の水位観測値データを取得する。
【００７３】
学習データ作成部１２２は、取得部１２１により取得された雨量観測値データ及び水位観測値データを用いて、水位予測モデルの学習に用いる学習データを作成する。
【００７４】
出力値算出部１２３は、学習データ作成部１２２により作成された学習データに含まれる雨量観測値を用いて、学習済みでない水位予測モデルにより、予測水位地点の過去の各日時における水位を示す出力値を算出する。なお、学習済みでない水位予測モデルとは、未学習又は学習途中の水位予測モデルのことである。
【００７５】
パラメータ更新部１２４は、出力値算出部１２３により算出された出力値と、学習データに含まれる水位観測値との同一日時（又は互いに対応する日時）における誤差に基づいて、水位予測モデルのパラメータを更新する。このとき、パラメータ更新部１２４は、既知の最適化手法により、当該誤差が最小となるように、水位予測モデルのパラメータを更新する。これにより、水位予測モデルが学習される。
【００７６】
次に、本実施形態に係る水位予測処理部１３０の機能構成について、図１１を参照しながら説明する。図１１は、本実施形態に係る水位予測処理部１３０の機能構成の一例を示す図である。
【００７７】
図１１に示すように、本実施形態に係る水位予測処理部１３０には、取得部１３１と、水位予測部１３２とが含まれる。
【００７８】
取得部１３１は、雨量観測値記憶部２５０に記憶されている雨量観測値データのうち、選択雨量観測地点の雨量観測値データを取得する。また、取得部１３１は、選択雨量観測地点の予測時点の１つ前の日時までの雨量予測値を取得する。
【００７９】
水位予測部１３２は、雨量観測値データが示す雨量観測値と、雨量予測値とを用いて、学習済みの水位予測モデルにより、予測水位地点の予測時点における水位を示す水位予測値を算出する。
【００８０】
次に、雨量観測地点選択処理部１１０が選択雨量観測地点を選択する処理（雨量観測地点の選択処理）について、図１２を参照しながら説明する。図１２は、本実施形態に係る雨量観測地点の選択処理の一例を示すフローチャートである。
【００８１】
まず、予測地点選択部１１１は、水位観測地点記憶部２１０に記憶されているＭ個の水位観測地点データがそれぞれ示すＭ個の水位観測地点のうち、水位を予測する水位観測地点（すなわち、水位予測地点）を選択する（ステップＳ１２０１）。なお、予測地点選択部１１１は、例えば、ユーザにより指定された水位観測地点を選択しても良いし、予め決められた水位観測地点を選択しても良い。
【００８２】
次に、流域特定部１１２は、標高情報記憶部２２０に記憶されている標高情報データと、分水嶺情報記憶部２３０に記憶されている分水嶺情報データとを用いて、水位予測地点の流域を特定する（ステップＳ１２０２）。ここで、流域特定部１１２は、例えば、水位予測地点よりも標高が高い地点で、かつ、当該地点と水位予測地点とを接続する線分が分水嶺と交差しない地点の集合を流域と特定する。
【００８３】
より具体的には、流域特定部１１２は、標高情報記憶部２２０に記憶されている標高情報データを参照して、水位予測地点の標高を特定した上で、水位予測地点よりも標高の高い地点を特定する。また、流域特定部１１２は、分水嶺情報記憶部２３０に記憶されている分水嶺情報データを参照して、分水嶺を表す各線分を特定する。そして、流域特定部１１２は、特定した地点のうち、当該地点と水位予測地点とを接続する線分が、分水嶺を表す線分と交差しない地点の集合を流域と特定する。
【００８４】
例えば、図１３に示すように、水位予測地点Ｑよりも標高が高い或る地点をＶ１及びＶ２とした場合、水位予測地点Ｑと地点Ｖ１と接続する線分Ｌ１は分水嶺と交差しない。一方で、水位観測地点Ｑと地点Ｖ２とを接続する線分Ｌ２は分水嶺と交差する。このため、地点Ｖ１は領域と特定される一方で、地点Ｖ２は流域と特定されない。水位予測地点Ｑの流域Ｄは、このような地点Ｖ１の集合で特定される。
【００８５】
なお、本実施形態では、標高情報データと分水嶺情報データとを用いて流域を特定する場合を説明したが、これに限られない。例えば、水位予測地点の流域を示す流域情報データが存在する場合は、流域特定部１１２は、当該流域情報データから流域を特定しても良い。また、稜線又は尾根が分水嶺である場合が多いため、流域特定部１１２は、分水嶺情報データの代わりに、例えば、稜線を示す稜線情報データ又は尾根を示す尾根情報データ等から分水嶺を表す各線分を特定することで、流域を特定しても良い。このように、標高情報や分水嶺情報、流域情報、稜線情報、尾根情報等の種々の地形情報を示すデータから、分水嶺や流域を特定しても良い。
【００８６】
次に、雨量観測地点選択部１１３は、雨量観測地点記憶部２４０に記憶されているＩ個の雨量観測地点データがそれぞれ示すＩ個の雨量観測地点のうち、特定された流域内に含まれる雨量観測地点（第１の雨量観測地点）を選択する（ステップＳ１２０３）。これは、水位予測地点の流域内に含まれる雨量観測地点は、当該水位予測地点との関連性が高いためである。
【００８７】
次に、雨量観測地点選択部１１３は、雨量観測地点記憶部２４０に記憶されているＩ個の雨量観測地点データがそれぞれ示すＩ個の雨量観測地点に対して、当該雨量観測地点が含まれ、半径ｒの円形領域を設定する（ステップＳ１２０４）。例えば、雨量観測地点選択部１１３は、Ｉ個の雨量観測地点に対して、当該雨量観測地点を中心として、半径ｒの円形領域を設定すれば良い。なお、ｒとしては任意の値を設定することができるが、例えば、１ｋｍ等と設定すれば良い。円形領域中央の雨量観測地点で観測された雨量と同等の雨量がこの円形領域内で降雨している、と仮定出来る。
【００８８】
次に、雨量観測地点選択部１１３は、特定した流域内に含まれない雨量観測地点の中から、当該雨量観測地点に設定された円形領域と流域とが重なる領域（この領域を「重畳領域」と表す。）の面積が所定の閾値以上となる雨量観測地点を抽出する（ステップＳ１２０５）。なお、閾値としては任意の値を設定することができるが、例えば、πｒ２／２乃至πｒ２／４等の範囲内の値を設定すれば良い。
【００８９】
次に、雨量観測地点選択部１１３は、抽出した雨量観測地点のうち、当該雨量観測地点の重畳領域から第１の雨量観測地点に設定された円形領域を除いた領域（この領域を「差分領域」と表す。）の面積が閾値以上である雨量観測地点（第２の雨量観測地点）を選択する（ステップＳ１２０６）。
【００９０】
例えば、図１４に示すように、雨量観測地点ｉ１〜ｉ８が存在するものとする。この場合、雨量観測地点ｉ１〜ｉ５は流域に含まれるため、これら雨量観測地点ｉ１〜ｉ５が第１の雨量観測地点となる。一方で、雨量観測地点ｉ６〜ｉ８は流域に含まないため、重畳領域が閾値以上かつ差分領域が閾値以上である場合に、第２の雨量観測地点となる。例えば、雨量観測地点ｉ６は差分領域が存在しないため、重畳領域の面積が閾値以上であれば第２の雨量観測地点となる。また、雨量観測地点ｉ７は重畳領域の面積と差分領域の面積との両方が閾値以上であれば（すなわち、差分領域が閾値以上であれば）第２の雨量観測地点となる。また、雨量観測地点ｉ８は、仮に重畳領域の面積が閾値以上であった場合であっても、差分領域の面積が０となるため、第２の雨量観測地点とならない。
【００９１】
このように、流域に含まれない雨量観測地点のうち、所定の条件（すなわち、流域との重畳領域が閾値以上で、かつ、第１の雨量観測地点の円形領域との差分領域が閾値以上）を満たす雨量観測地点を第２の雨量観測地点に選択する。これは、水位予測地点の流域内に含まれない雨量観測地点であっても、このような条件を満たす雨量観測地点は、当該水位予測地点との関連性が高いためである。
【００９２】
以上のように、本実施形態に係る水位予測装置１０は、Ｉ個の雨量観測地点の中から、水位予測地点と関連性の高い選択雨量観測地点（第１の雨量観測地点及び第２の雨量観測地点）を選択する。これにより、本実施形態に係る水位予測装置１０は、後述する水位予測モデルの学習処理と河川水位の予測処理との計算時間を削減させることができると共に、河川水位の予測処理における予測精度を向上させることができるようになる。以降では、選択雨量観測地点の総数をＪ、各選択雨量観測地点をｊで表すものとする。
【００９３】
なお、上述した雨量観測地点の選択処理では、ステップＳ１２０４で各雨量観測地点に対して円形領域を設定したが、これに限られない。例えば、各雨量観測地点に対して、楕円領域や矩形領域、多角形領域等が設定されても良い。
【００９４】
また、図１５に示すように、選択雨量観測地点ｊの総数Ｊが多い程水位予測モデルの学習処理と河川水位の予測処理との計算時間が増大する一方で、Ｊが或る一定数を超えると誤差（予測誤差）の向上が認められないことがある。そこで、選択雨量観測地点ｊの総数Ｊが所定の値（又は所定の値以下）となるようにしても良い。このような所定の値は流域の面積等によって異なるが、例えば、図１５に示すように、Ｊ＝５やＪ＝６等とすれば良い。なお、予測誤差が低い程、予測精度が高いことを意味する。
【００９５】
また、例えば、第１の雨量観測地点の円形領域が流域を所定の割合以上覆っているような場合には、第２の雨量観測地点を選択しなくても良い。この場合、第１の雨量観測地点で観測される雨量のみでも高い予測精度が期待できるためである。
【００９６】
次に、モデル学習処理部１２０が水位予測モデルを学習する処理（水位予測モデルの学習処理）について、図１６を参照しながら説明する。図１６は、本実施形態に係る水位予測モデルの学習処理の一例を示すフローチャートである。
【００９７】
まず、取得部１２１は、雨量観測値記憶部２５０に記憶されている雨量観測値データのうち、各選択雨量観測地点ｊの雨量観測値データを取得する（ステップＳ１６０１）。すなわち、取得部１２１は、Ｊ個の選択雨量観測地点ｊの各日時における雨量観測値データを雨量観測値記憶部２５０からそれぞれ取得する。ここで、各日時をｔで表して、日時ｔ＝１〜Ｔまでの雨量観測値データが雨量観測値記憶部２５０に記憶されているとした場合、取得部１２１により取得された雨量観測値データが示す雨量観測値はｘｊ，ｔ（ｊ＝１，・・・，Ｊ，ｔ＝１，・・・，Ｔ）と表すことができる。
【００９８】
次に、取得部１２１は、水位観測値記憶部２６０に記憶されている水位観測値データのうち、予測水位地点の水位観測値データを取得する（ステップＳ１６０２）。このとき、取得部１２１は、上記のステップＳ１６０１で取得した雨量観測値データの各日時に対応する日時の水位観測値データを取得する。すなわち、取得部１２１は、予測水位地点ｍ´の各日時ｔ＝１〜Ｔまでの水位観測値データを水位観測値記憶部２６０から取得する。ここで、この場合、取得部１２１により取得された水位観測値データが示す水位観測値はｙｍ´，ｔ（ｔ＝１，・・・，Ｔ）と表すことができる。以降では、この水位観測値ｙｍ´，ｔを「ｙｔ」と表す。
【００９９】
次に、学習データ作成部１２２は、取得した雨量観測値データが示す雨量観測値ｘｊ，ｔと、水位観測値データが示す水位観測値はｙｔとを用いて、水位予測モデルの学習に用いる学習データを作成する（ステップＳ１６０３）。出力値算出部１２３は、例えば、予め設定されたＫを用いて、日時ｔ毎に、学習データ｛ｘｊ，ｔ−Ｋ，・・・，ｘｊ，ｔ−１，ｙｔ｝を作成する。
【０１００】
ここで、Ｋは、或る日時ｔの水位を予測する場合に、この日時ｔよりも過去の雨量観測値ｘｊ，ｔを用いる範囲を決定するハイパーパラメータである。このような範囲としては、例えば、数時間、数十時間、数日又は数十日等が挙げられる。例えば、日時ｔが１０分毎であり、或る日時ｔから１０時間前までの範囲の過去の雨量観測値ｘｊ，ｔを用いる場合、Ｋ＝６０となる。なお、学習データに含まれるｙｔが、学習済みでない水位予測モデルの日時ｔにおける出力値に対する正解データとなる。
【０１０１】
次に、出力値算出部１２３は、日時ｔ毎に、学習データに含まれる雨量観測値ｘｊ，ｔ−Ｋ，・・・，ｘｊ，ｔ−１を用いて、以下の式（１）に示す水位予測モデルにより、予測水位地点ｍ´の出力値を算出する（ステップＳ１６０３）。
【０１０２】
【数１】
2020002658.tif 000003
【０１０３】
ここで、αｊは選択雨量観測地点ｊで観測される雨量観測値に対する水位予測地点ｍ´における水位の感度を示すパラメータ、Ｄｊは選択雨量観測地点ｊの時間遅れを示すパラメータ、λｊは選択雨量観測地点ｊの時定数を示すパラメータである。これらのパラメータαｊ、Ｄｊ及びλｊが学習対象のパラメータである。なお、これらのパラメータαｊ、Ｄｊ及びλｊは適切な初期値に初期化されている。
【０１０４】
更に、φは以下の式（２）で表される減衰関数である。
【０１０５】
【数２】
2020002658.tif 000004
【０１０６】
なお、感度は、例えば、選択雨量観測地点ｊで単位雨量が観測された場合に、水位予測地点ｍ´でどの程度水位が上昇するかを表す。また、時間遅れＤｊは、選択雨量観測地点ｊで観測された雨量が、水位予測地点ｍ´に流れ込むまでに要する時間を表す。また、時定数λｊは、上記の式（２）の減衰関数に用いられ、選択雨量観測地点ｊで観測された雨量の水位上昇への影響がどの程度早く減衰していくかを表す。
【０１０７】
上記の式（１）に示すように、本実施形態では、雨量に対する水位応答曲線として減衰関数を用いたモデルを水位予測モデルとしている。
【０１０８】
次に、パラメータ更新部１２４は、上記のステップＳ１６０３で算出された各日時ｔの出力値と、当該日時ｔに対応する学習データに含まれる正解データｙｔとを用いて、以下の式（３）により誤差を算出する。そして、パラメータ更新部１２４は、既知の最適化手法により、当該誤差が最小となるように、水位予測モデルのパラメータを更新する（ステップＳ１６０４）。
【０１０９】
【数３】
2020002658.tif 000005
【０１１０】
これにより、水位予測地点ｍ´の水位を予測するための水位予測モデルが学習される。
【０１１１】
以上のように、本実施形態に係る水位予測装置１０は、Ｊ（＜Ｉ）個の選択雨量観測地点の雨量観測値と、水位予測地点の水位観測値とを用いて、水位予測モデルを学習する。このように、水位予測モデルの学習に用いる雨量観測値を、水位予測地点と関連性の高い雨量観測地点の雨量観測値とすることで、本実施形態に係る水位予測装置１０は、水位予測モデルの学習に要する計算時間の削減と、予測精度の向上とを図ることができる。
【０１１２】
なお、例えば、水位予測モデルの過学習（又は過適合）を防止するため、Ｊ個の選択雨量観測地点の中から、学習に用いる選択雨量観測地点を更に選択しても良い。言い換えれば、上記の式（１）において、任意の個数のαｊをαｊ＝０としても良い。αｊ＝０とするｊは、例えば、交差検定等を用いた水位予測モデルの評価結果から決定されれば良い。
【０１１３】
また、上述した水位予測モデルの学習処理では或る水位予測地点ｍ´の水位予測モデルを学習する場合について説明した。複数の水位予測地点ｍ´の水位予測モデルをそれぞれ学習する場合、例えば、図１２のステップＳ１２０１で選択した水位予測地点ｍ´毎に、図１２のステップＳ１２０２〜ステップＳ１２０６の処理と、水位予測モデルの学習処理とを繰り返し実行すれば良い。
【０１１４】
次に、学習済みの水位予測モデルを用いて、水位予測処理部１３０が水位予測地点の水位を予測する処理（河川水位の予測処理）について、図１７を参照しながら説明する。図１７は、本実施形態に係る河川水位の予測処理の一例を示すフローチャートである。以降では、予測時点Ｔ´、現在を示す日時をＳと表す。
【０１１５】
まず、取得部１３１は、雨量観測値記憶部２５０に記憶されている雨量観測値データのうち、各選択雨量観測地点ｊの雨量観測値ｘｊ，Ｔ´−Ｋ，ｘｊ，Ｔ´−Ｋ−１，・・・，ｘｊ，Ｓがそれぞれ含まれる雨量観測値データを取得する（ステップＳ１７０１）。なお、Ｋの値は、モデル学習時に用いた値と同様である。なお、雨量観測値記憶部２５０に記憶されている雨量観測値データを水位予測処理部１３０に入力する際に、雨量観測地点選択処理部１１０に雨量観測値データを選択させても良い。選択処理は学習時と同様である。
【０１１６】
次に、取得部１３１は、各選択雨量観測地点ｊの現在から予測時点の１つ前の日時までの雨量予測値ｘｊ，Ｓ＋１，・・・，ｘｊ，Ｔ´−１を取得する（ステップＳ１７０２）。なお、このような雨量予測値ｘｊ，Ｓ＋１，・・・，ｘｊ，Ｔ´−１として、取得部１３１は、例えば、気象庁等により提供される雨量の予測値を取得すれば良い。
【０１１７】
次に、水位予測部１３２は、取得された雨量観測値ｘｊ，Ｔ´−Ｋ，ｘｊ，Ｔ´−Ｋ−１，・・・，ｘｊ，Ｓと、雨量予測値ｘｊ，Ｓ＋１，・・・，ｘｊ，Ｔ´−１とを用いて、以下の式（４）に示す水位予測モデルにより、水位予測値を算出する（ステップＳ１７０３）。
【０１１８】
【数４】
2020002658.tif 000006
【０１１９】
ここで、αｊ、Ｄｊ及びλｊは学習済みの水位予測モデルのパラメータである。これにより、水位予測地点ｍ´の予測時点Ｔ´における水位を示す値（水位予測値）が得られる。
【０１２０】
以上のように、本実施形態に係る水位予測装置１０は、Ｊ（＜Ｉ）個の選択雨量観測地点の雨量（雨量観測値及び雨量予測値）を用いて、学習済みの水位予測モデルにより水位予測値を算出する。このように、水位予測モデルの学習に用いる雨量観測値を、水位予測地点と関連性の高い雨量観測地点の雨量観測値とすることで、本実施形態に係る水位予測装置１０は、予測精度の向上を図ることができる。
【０１２１】
また、本実施形態に係る水位予測装置１０は、水位観測値を用いることなく、雨量（雨量観測値及び雨量予測値）から水位を予測することができる。これにより、例えば、水位観測地点に設置されている水位観測機器の不具合等によって過去の水位を示すデータに欠損があった場合であっても、水位を予測することができるようになる。
【０１２２】
ここで、本実施形態に係る水位予測装置１０により予測した水位予測値の予測精度について、図１８を参照しながら説明する。図１８は、予測精度の一例を説明するための図である。
【０１２３】
図１８は、或る水位観測地点の各日時における水位観測値と水位予測値とを示している。また、図１８の雨量観測値は、当該水位観測地点と関連性が高い雨量観測地点で観測された雨量を示している。図１８に示すように、本実施形態に係る水位予測装置１０によれば、各日時における水位観測値と水位予測値との誤差が少なく、高い精度で予測できていることがわかる。
【０１２４】
以上、本発明の実施の形態について詳述したが、本発明は斯かる特定の実施形態に限定されるものではなく、特許請求の範囲に記載された本発明の要旨の範囲内において、種々の変形や変更が可能である。
【符号の説明】
【０１２５】
１０        水位予測装置
１１０      雨量観測地点選択処理部
１１１      予測地点選択部
１１２      流域特定部
１１３      雨量観測地点選択部
１２０      モデル学習処理部
１２１      取得部
１２２      学習データ作成部
１２３      出力値算出部
１２４      パラメータ更新部
１３０      水位予測処理部
１３１      取得部
１３２      水位予測部
２１０      水位観測地点記憶部
２２０      標高情報記憶部
２３０      分水嶺情報記憶部
２４０      雨量観測地点記憶部
２５０      雨量観測値記憶部
２６０      水位観測値記憶部
【図１】
2020002658.tif 000007
【図２】
2020002658.tif 000008
【図３】
2020002658.tif 000009
【図４】
2020002658.tif 000010
【図５】
2020002658.tif 000011
【図６】
2020002658.tif 000012
【図７】
2020002658.tif 000013
【図８】
2020002658.tif 000014
【図９】
2020002658.tif 000015
【図１０】
2020002658.tif 000016
【図１１】
2020002658.tif 000017
【図１２】
2020002658.tif 000018
【図１３】
2020002658.tif 000019
【図１４】
2020002658.tif 000020
【図１５】
2020002658.tif 000021
【図１６】
2020002658.tif 000022
【図１７】
2020002658.tif 000023
【図１８】
2020002658.tif 000024