JP2009167751A - Inflow forecasting system, inflow forecasting method and program - Google Patents
Inflow forecasting system, inflow forecasting method and program Download PDFInfo
- Publication number
- JP2009167751A JP2009167751A JP2008009568A JP2008009568A JP2009167751A JP 2009167751 A JP2009167751 A JP 2009167751A JP 2008009568 A JP2008009568 A JP 2008009568A JP 2008009568 A JP2008009568 A JP 2008009568A JP 2009167751 A JP2009167751 A JP 2009167751A
- Authority
- JP
- Japan
- Prior art keywords
- amount
- period
- inflow
- value
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000013277 forecasting method Methods 0.000 title description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 41
- 238000000611 regression analysis Methods 0.000 claims abstract description 23
- 239000005413 snowmelt Substances 0.000 claims description 85
- 238000001556 precipitation Methods 0.000 claims description 69
- 238000002844 melting Methods 0.000 claims description 35
- 230000008018 melting Effects 0.000 claims description 35
- 238000000034 method Methods 0.000 claims description 18
- 238000009825 accumulation Methods 0.000 claims description 13
- 238000013179 statistical model Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 239000000155 melt Substances 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
Description
本発明は、貯水施設への水の流入量を予測するシステム、方法及びプログラムに関する。 The present invention relates to a system, method, and program for predicting the amount of water flowing into a water storage facility.
ダム等の貯水施設における水位管理や河川の防災対策などのために、貯水施設への水の流入量や河川における水の流量(以下、これらを流入量という。)を予測することが行われている。例えば、特許文献1では、予報雨量の予報誤差の傾向を解析し、その結果を用いて予報雨量を修正し、現在の流量、現在の雨量、及び修正予報雨量を用いて将来の流量を予測している。
しかしながら、特許文献1では河川等の流入に影響を与える多くの地点の降雨量を観察し、降雨が流入量の計測地点に変化を生じさせる時間も考慮に入れて、過去の累積降雨量を説明変数として流入量を推計しており、複雑な推計が行われている。
However,
本発明は、このような背景を鑑みてなされたものであり、貯水施設への水の流入量を容易に予測することのできる、流入量予測システム、流入量予測方法及びプログラムを提供することを目的とする。 The present invention has been made in view of such a background, and provides an inflow amount prediction system, an inflow amount prediction method, and a program capable of easily predicting the inflow amount of water into a water storage facility. Objective.
上記課題を解決するための本発明のうち請求項1に記載の発明は、貯水施設への水の流入量を予測するシステムであって、所定期間ごとに、前記流入量の実績値及び降雨量の実績値を記憶するデータベースと、ある第1期間における前記降雨量、及び、前記第1期間から所定期前の第2期間における前記流入量と所定の定数である均衡流入量との差を説明変数とし、前記第2期間から前記第1期間までの前記流入量の増加量を目的変数とする回帰モデルを記憶するモデル記憶部と、前記データベースに記憶されている前記流入量の実績値及び前記降雨量の実績値、並びに前記モデル記憶部に記憶されている前記回帰モデルに基づいて回帰分析を行い、前記説明変数のそれぞれについての回帰係数を推計する流入量モデル推計部と、予測対象の期間である予測期間における前記降雨量の予測値を取得する予測降雨量取得部と、前記予測期間から前記所定期前の期間における前記流入量の実績値を前記データベースから取得し、前記降雨量の予測値、前記流入量の実績値、及び前記各回帰係数を前記回帰モデルに適用して、前記予測期間における前記増加量の予測値を算出する増加量予測部と、前記流入量の実績値に前記増加量の予測値を加算して、前記予測期間における前記流入量の予測値を算出する流入量予測部と、を備えることとする。
The invention according to
前記所定期間とは、例えば、1時間や1日などの単位期間である。前記所定期は、例えば、系列相関がない場合には1期とすることができる。また、前記所定期は、系列相関がある場合には、例えば、Prais-Winstein変換やコクラン・オーカット法により変数変換した期間の数とすることができる。 The predetermined period is a unit period such as one hour or one day. For example, the predetermined period may be one period when there is no series correlation. Further, when there is a series correlation, the predetermined period can be the number of periods subjected to variable conversion by, for example, Prais-Winstein conversion or the Cochran-Ocut method.
本発明の流入量予測システムによれば、流入量の実績値に基づいて、流入量を予測することができるので、降雨量のみに基づいて流入量を予測する場合に比べ、より容易かつ精度よく流入量の予測を行うことができる。 According to the inflow amount prediction system of the present invention, since the inflow amount can be predicted based on the actual value of the inflow amount, it is easier and more accurate than the case where the inflow amount is predicted based only on the rainfall amount. The amount of inflow can be predicted.
また、本発明のうち請求項2に記載の発明は、請求項1に記載の流入量予測システムであって、前記データベースは、前記所定期間ごとに、融雪量の実績値をさらに記憶し、前記モデル記憶部に記憶される前記回帰モデルには、説明変数として、さらに前記第1期間における前記融雪量が含まれ、前記予測期間における前記融雪量の予測値を取得する予測融雪量取得部を備え、前記流入量モデル推計部は、前記データベースに記憶されている前記流入量の実績値、前記降雨量の実績値、及び前記融雪量の実績値、並びに前記モデル記憶部に記憶されている前記回帰モデルに基づいて回帰分析を行い、前記回帰係数を推計することとする。
Moreover, invention of
この場合、融雪量の予測値を考慮して貯水施設への水の流入量を予測することができるので、降水量の予測のみに基づいて流入量を予測する場合に比べて、降水がなかった場合にも融雪により流入量が増加しうることを盛り込んで予測を行うことができる。したがって、より精度の高い予測を行うことができる。 In this case, the amount of water flowing into the water storage facility can be predicted in consideration of the predicted value of the amount of snowmelt, so there was no precipitation compared to the case where the amount of inflow was predicted based solely on the prediction of precipitation. Even in such a case, it is possible to predict that the inflow amount may increase due to melting snow. Therefore, more accurate prediction can be performed.
また、本発明のうち請求項3に記載の発明は、請求項2に記載の流入量予測システムであって、前記データベースはさらに、気温の実績値及び積雪量の実績値を記憶し、前記モデル記憶部は、前記回帰モデルである流入量モデルに加えて、前記第1期間における気温から融雪の始まる気温を示す定数を引いた値又はゼロの大きい方の値に、前記第2期間における前記積雪量を乗じた値と、前記第1期間における前記降雨量とを説明変数とし、前記第1期間における前記融雪量を目的変数とする回帰モデルである融雪量モデルを記憶し、前記データベースに記憶されている前記気温の実績値、前記積雪量の実績値及び前記降雨量の実績値、並びに前記モデル記憶部に記憶されている前記融雪量モデルに基づいて回帰分析を行い、前記融雪量モデルに係る前記説明変数のそれぞれについての回帰係数を推計する融雪量モデル推計部と、前記予測期間における前記気温の予測値を取得する予測気温取得部と、を備え、前記予測融雪量取得部は、前記予測期間の前記所定期前の期間における前記積雪量の実績値を前記データベースから取得し、前記気温の予測値、前記取得した積雪量の実績値、前記降雨量の予測値、及び前記融雪量モデルに係る前記各説明変数を、前記融雪量モデルに適用して、前記予測期間における前記融雪量の予測値を算出することとする。
Moreover, invention of
ある気温以上では気温が上昇するほど融雪量は多くなり、また降雨によっても雪は溶かされるので、気温、積雪量及び降雨量を説明変数として、融雪量を回帰させることができる。一般的に気温や降雨量は気象予報の手法を用いて容易に予測することが可能であり、これらの容易に取得可能な気温や降雨量の予測値と、積雪量の実績値とに基づいて、融雪量を予測することができるので、融雪量そのものの予測が困難である場合にも対応することができる。 Above a certain temperature, the amount of snow melt increases as the temperature rises, and the snow melts due to rain. Therefore, the amount of snow melt can be regressed using the temperature, the amount of snow and the amount of rainfall as explanatory variables. In general, temperature and rainfall can be easily predicted using a weather forecasting method. Based on these easily obtainable predicted values of temperature and rainfall, and actual values of snow cover. Since the snow melting amount can be predicted, it is possible to cope with the case where it is difficult to predict the snow melting amount itself.
また、本発明のうち請求項4に記載の発明は、請求項3に記載の流入量予測システムであって、前記データベースはさらに、降雪量の実績値を記憶し、前記融雪モデル推計部は、前記積雪量モデルに含まれる前記積雪量の説明変数に、式
前記第1期間における積雪量=max(前記第1期間における前記降雪量+前記第2期間における前記積雪量−前記第1期間における前記融雪量,0)
を代入したモデルを回帰分析することにより、前記融雪量モデルに係る前記説明変数のそれぞれについての前記回帰係数を推計することとする。
Moreover, invention of
The regression coefficient for each of the explanatory variables related to the snowmelt amount model is estimated by performing regression analysis on the model in which is substituted.
また、本発明のうち請求項5に記載の発明は、請求項1に記載の流入量予測システムであって、前記データベースには、前記降雨量の実績値に代えて、降水量の実績値が記憶されており、前記予測期間における前記降水量の予測値を取得する予測降水量取得部と、前記予測期間における前記気温の予測値を取得する予測気温取得部と、を備え、前記予測降雨量取得部は、前記気温の予測値が所定の気温を超えている場合に前記降水量を前記降雨量とし、前記気温の予測値が前記所定の気温以下である場合には前記降雨量を0とすることとする。
Moreover, invention of
また、本発明うち請求項6に記載の発明は、請求項5に記載の流入量予測システムであって、前記データベースは、前記降雨量の実績値に代えて、降水量の実績値を記憶し、前記データベースはさらに、気温の実績値及び降雪量の実績値を記憶し、前記気温、及び前記降水量を説明変数とし、前記降雪量を目的変数とした回帰モデルである降雪モデル、並びに前記データベースに記憶されている前記気温の実績値、前記降水量の実績値、及び前記降雪量の実績値に基づいて回帰分析を行い、前記説明変数のそれぞれについての回帰係数を推計する降雪モデル推計部と、前記データベースに記憶されている前記気温の実績値のそれぞれについて、前記気温がδ以下の場合には、前記降雪モデルに前記推計した回帰係数と前記気温の実績値及び前記降水量の実績値とを前記降雪モデルに適用して算出した値を推計降雪量とし、前記気温がδより高い場合には0を前記推計降雪量としたときに、前記降雪量の実績値と前記推計降雪量との誤差が最も小さくなる前記δを前記所定の気温として決定する降雪気温推計部と、を備えることとする。
Further, the invention according to
また、本発明のうち請求項7に記載の発明は、貯水施設への水の流入量を予測するシステムであって、所定期間ごとに、前記流入量の実績値、降雨量の実績値、気温の実績値、及び積雪量の実績値を記憶するデータベースと、ある第1期間における前記降雨量、前記第1期間における融雪量、及び前記第1期間から所定期前の第2期間における前記流入量と所定の定数である均衡流入量との差を説明変数とし、前記第2期間から前記第1期間までの前記流入量の増加量を目的変数とする回帰モデルである流入量モデル、及び、前記第1期間における気温から融雪の始まる気温を示す定数を引いた値又はゼロの大きい方の値に、前記第2期間における前記積雪量を乗じた値と、前記第1期間における前記降雨量とを説明変数とし、前記第1期間における前記融雪量を目的変数とする回帰モデルである融雪量モデルを記憶するモデル記憶部と、前記データベースに記憶されている前記気温の実績値、前記積雪量の実績値及び前記降雨量の実績値、並びに前記モデル記憶部に記憶されている前記融雪量モデルに基づいて回帰分析を行い、前記融雪量モデルに係る前記説明変数のそれぞれについての回帰係数を推計する融雪量モデル推計部と、前記データベースに記憶されている前記気温の実績値、前記積雪量の実績値、前記降雨量の実績値、及び前記推計した回帰係数を前記融雪量モデルに適用して、前記所定期間ごとの前記融雪量の理論値を算出し、前記データベースに記憶されている前記流入量の実績値、及び前記降雨量の実績値、前記算出した融雪量の理論値、並びに前記モデル記憶部に記憶されている前記流入量モデルに基づいて回帰分析を行い、前記流入量モデルに係る前記説明変数のそれぞれについての回帰係数を推計する流入量モデル推計部と、予測対象の期間である予測期間における前記降雨量の予測値を取得する予測降雨量取得部と、前記予測期間における前記気温の予測値を取得する予測気温取得部と、前記予測期間の前記所定期前の期間における前記積雪量の実績値を前記データベースから取得し、前記気温の予測値、前記取得した積雪量の実績値、前記降雨量の予測値、及び前記融雪量モデルに係る前記各説明変数を、前記融雪量モデルに適用して、前記予測期間における前記融雪量の予測値を算出する予測融雪量取得部と、前記予測期間から前記所定期前の期間における前記流入量の実績値を前記データベースから取得し、前記降雨量の予測値、前記融雪量の予測値、前記流入量の実績値、及び前記各回帰係数を前記流入量モデルに適用して、前記予測期間における前記増加量の予測値を算出する増加量予測部と、前記流入量の実績値に前記増加量の予測値を加算して、前記予測期間における前記流入量の予測値を算出する流入量予測部と、を備えることとする。
Further, the invention according to
この場合、気温、降雨量、及び積雪量を説明変数とする融雪量モデルを用いて融雪量の理論値を算出し、その融雪量を考慮して、貯水施設への水の流入量を予測することができる。したがって、融雪量の測定値などを気象庁や民間気象会社などから取得することができない場合でも、一般的に取得が容易な気温や降雨量の予測値と積雪量の実績値とに基づいて融雪量の予測値を算出することができる。これにより、融雪量の予測値が容易に取得できない場合にも、降水がなかった場合にも融雪により流入量が増加しうることを盛り込んで、貯水施設への水の流入量の予測を行うことが可能となり、精度の高い流入量の予測を実現できる。 In this case, the theoretical value of the amount of snowmelt is calculated using a snowmelt amount model with temperature, rainfall, and snow cover as explanatory variables, and the amount of water flowing into the water storage facility is predicted in consideration of the amount of snowmelt. be able to. Therefore, even if measured values of snowmelt cannot be obtained from the Japan Meteorological Agency or private weather companies, etc., the amount of snowmelt is generally based on the predicted values of temperature and rainfall that are easy to obtain and the actual value of snow cover. The predicted value can be calculated. This makes it possible to predict the inflow of water into the water storage facility, taking into account that the inflow can increase due to snow melting, even if the predicted value of the amount of snow melt cannot be easily obtained or when there is no precipitation. Therefore, it is possible to predict the inflow amount with high accuracy.
その他本願が開示する課題やその解決方法については、発明の実施形態の欄及び図面により明らかにされる。 Other problems and solutions to be disclosed by the present application will be made clear by the embodiments of the invention and the drawings.
本発明によれば、貯水施設への水の流入量を容易に予測することができる。 According to the present invention, it is possible to easily predict the amount of water flowing into a water storage facility.
以下、本発明の一実施形態に係る流入量予測システム10について説明する。
本実施形態の流入量予測システム10は、ダムなどの貯水施設に河川から流入する水量(以下、流入量という。)を予測する。流入量予測システム10は、過去の降水量、降雪量、積雪量、融雪量、流入量を、後述する統計モデルを用いて回帰分析することにより回帰係数を推計し、気温及び降水量の予測値(例えば、気象予報による予報値を採用することができる。)と、推計した回帰係数とを統計モデルに適用して流入量を予測する。なお、流入量の予測は一日単位で行うものとする。
Hereinafter, an inflow
The inflow
図1〜3は、流入量の変化を示すグラフである。降水や融雪などがない場合にも、例えば山林などからの滲出によって、所定の流入量は存在する。したがって、降水や融雪がないと、流入量は所定の均衡値(以下、均衡流入量という。)に逓減していく(図1)。これに対し、降水があると、それに応じて流水量は一時的に増加するが、降水が止むとともに、再度均衡流入量に向けて逓減を始める(図2)。一方、気温が上昇すると融雪が発生し、それに応じて流水量も増加するが、融雪は降水に比べて流入量に与える影響の変化が緩やかである(図3)。これは、例えば冬季から春季に向けての時期などにおいて、平均気温が上昇している場合に、継続的に融雪が発生するような場合である。 1 to 3 are graphs showing changes in the inflow amount. Even when there is no precipitation or snowmelt, there is a predetermined inflow due to, for example, oozing from a forest. Therefore, if there is no precipitation or snowmelt, the inflow will gradually decrease to a predetermined equilibrium value (hereinafter referred to as equilibrium inflow) (Fig. 1). On the other hand, when there is precipitation, the amount of water flowing temporarily increases accordingly, but the precipitation stops and begins to decrease again toward the equilibrium inflow (Fig. 2). On the other hand, when the temperature rises, snow melts and the amount of water flowing increases accordingly, but the effect of snow melt on the inflow is more gradual than that of precipitation (Fig. 3). This is a case where, for example, snow melting continuously occurs when the average temperature rises during the period from winter to spring.
そこで、本実施形態では、降水量、融雪量及び均衡流入量に着目し、流入量を算出する回帰モデルを用いて流入量を予測する。本実施形態の流入量予測システム10は、過去の気温や降水量などの気象データの実績値と後述する回帰モデルとに基づいてパラメータを推計し、推計したパラメータと、例えば気象予報などにより求められる気温の予測値(以下、予測気温という。)及び降水量の予測値(以下、予測降水量という。)とを回帰モデルに適用して流入量の予測値を算出する。以下、詳細について説明する。
Therefore, in this embodiment, paying attention to precipitation, snowmelt, and balanced inflow, the inflow is predicted using a regression model that calculates the inflow. The
==システム構成==
本実施形態では、流入量予測システム10は、1台のコンピュータにより構成されることを想定している。図4は、本実施形態の流入量予測システム10のハードウェア構成を示す図である。流入量予測システム10は、CPU11、メモリ12、記憶装置13、入力装置14及び出力装置15を備える。記憶装置13は、各種のプログラムやデータを記憶する、例えば、ハードディスクドライブやフラッシュメモリ、CD−ROMドライブなどである。CPU11は、記憶装置13に記憶されているプログラムをメモリ12に読み出して実行することにより各種の機能を実現する。入力装置14は、ユーザからデータの入力を受け付ける、例えばキーボードやマウス、タッチパネル、マイクロフォンなどである。出力装置15は、データを出力する、例えばディスプレイやプリンタ、スピーカなどである。
== System configuration ==
In the present embodiment, it is assumed that the inflow
図5は、本実施形態の流入量予測システム10のソフトウェア構成を示す図である。本実施形態の流入量予測システム10は、降雪気温推計部21、融雪量モデル推計部22、流入量モデル推計部23、予測気温取得部24、予測降水量取得部25、流入量予測部26、モデル記憶部31、パラメータ記憶部32、及び気象実績データベース40を備えている。
FIG. 5 is a diagram illustrating a software configuration of the inflow
なお、降雪気温推計部21、融雪量モデル推計部22、流入量モデル推計部23、予測気温取得部24、予測降水量取得部25及び流入量予測部26は、流入量予測システム10のCPU11が、記憶装置13に記憶されているプログラムをメモリ12に読み出して実行することにより実現される。
The snowfall
気象実績データベース40には、気象の各種実績値を含む情報(以下、気象実績情報という。)の履歴が記憶される。図6は、気象実績データベース40に記憶される気象実績情報の構成例を示す図である。同図に示すように、気象実績情報には、日付に対応付けて、気温、降水量、降雪量、積雪量、融雪量、流入量が含まれている。気温は、一日の平均気温である。降水量、降雪量、融雪量は一日の降水量、降雪量、融雪量の累積値である。積雪量は、その日に観測された積雪量である。流入量は、一日にダムなどの貯水施設に流入した水量の累計値である。気温、降水量、降雪量及び積雪量は、例えば、気象庁や民間気象会社などが提供するデータである。流入量は、貯水施設において測定した測定値である。流入量は、例えば、河川を管理する自治体などが提供する河川の流量の測定値としてもよい。融雪量は、気象庁や民間気象会社、測量会社などが測定したものであってもよいし、後述するモデルにより計算した値を実績値として記録するようにしてもよい。
The
モデル記憶部31には、気象実績情報に基づく各種の統計モデルが記憶され、パラメータ記憶部32には、モデル記憶部31に記憶されているモデルに適用される回帰係数や定数などのパラメータが記憶される。モデル記憶部31には、降水量の統計モデル(以下、降水量モデル1という。)、降雪量の統計モデル(以下、降雪量モデル2という。)、積雪量の統計モデル(以下、積雪量モデル3という。)、融雪量の統計モデル(以下、融雪量モデル4という。)、及び前日からの流入量の増加量に係る統計モデル(以下、流入量モデル5という。)が記憶される。パラメータ記憶部32には、降雨が降雪に変わる気温(以下、降雪気温という。)δ、均衡流入量μ1、融雪が始まる気温μ2、回帰係数α1〜α3、β1、β2がパラメータ記憶部32に記憶される。なお、均衡流入量μ1及び融雪が始まる気温μ2については、所定の定数として、予めパラメータ記憶部32に記憶されているものとする。
Various statistical models based on weather performance information are stored in the
なお以下の説明において、日付tにおける気温、降水量、降雪量、積雪量、融雪量及び流入量をそれぞれ、Tt、Pt、St、Dt、Mt及びFtとする。前日からの流入量の増加量をΔFtとする。気温Ttがδより高い場合の降水量をP1t、気温Ttがδ以下である場合の降水量をP2tとする。 In the following description, the temperature at the date t, rainfall, snowfall, snowfall, snowmelt and inflow respectively, and T t, P t, S t , D t, M t , and F t. Let ΔFt be the increase in the inflow from the previous day. The precipitation when the temperature Tt is higher than δ is P1 t , and the precipitation when the temperature T t is δ or less is P2 t .
降水量モデル1は、降水量Ptが、降雪気温δを境に、P1tまたはP2tとなることを示すモデルであり、次式で表される。
降雪量モデル2は、気温Ttがδ以下である場合の降水量P2tを説明変数とし、降雪量Stを目的変数とした回帰モデルであり、次式で表される。
積雪量モデル3は、積雪量Dtが、前日までの積雪量Dt−1に当日の降雪量Stを加え、そこから融雪量Mtを引いたものに一致するという関係を示すモデルであり、次式で表される。
融雪量モデル4は、気温Ttから融雪が始まる気温μ2を減じた値に積雪量Dtを乗じた値と、P1tとを説明変数とし、融雪量Mtを目的変数とした回帰モデルであり、次式で表される。
融雪量モデル4において、気温Ttがμ2よりも低ければ第1項は0になり、前日までの積雪量Dt−1が0であれば融雪量Mtは0になる。
流入量モデル5は、均衡流入量μ1から前日の流入量Ft−1を減じた値と、当日の気温Ttがδより高い場合の降水量P1tと、当日の融雪量Mtとを説明変数とし、流入量の増加量ΔFtを目的変数とした回帰モデルであり、次式で表される。
The
The snow
In the
降雪気温推計部21は、降水量モデル1及び気象実績情報に基づいて降雪気温δを推計し、推計した降雪気温δをパラメータ記憶部32に登録する。なお、降雪気温推計部21によるδの推計処理の詳細については後述する。
The snowfall
また、降雪気温推計部21は、各日付tについて、日付tに対応する気象実績情報を気象実績データベース40から読み出し、読み出した気象実績情報の降水量及び気温と、上記推計した降雪気温δを前記降水量モデル1に適用してP2tを算出する。降雪気温推計部21は、気象実績情報、P2t、及び降雪量モデル2に基づいて回帰分析を行い、回帰変数γを推計する。降雪気温推計部21は、推計した回帰変数γをパラメータ記憶部32に登録する。
Moreover, the snowfall
融雪量モデル推計部22は、積雪量モデル3を融雪量モデル4に代入した式
を回帰分析して、回帰係数α2及びα3を推計する。融雪量モデル推計部22は、推計した回帰係数α2及びα3をパラメータ記憶部32に登録する。
The snow melting amount
And regression coefficients α 2 and α 3 are estimated. The snowmelt amount
流入量モデル推計部23は、流入量モデル5及び気象実績情報に基づいて回帰係数α1、β1、及びβ2を推計する。具体的には、流入量モデル推計部23は、各日付tについて、降水量モデル1、δ及び日付tに対応する気象実績情報の降水量に基づいてP1tを算出し、各日付tの気象実績情報と、対応するP1tを用いて、流入量モデル5を回帰分析し、回帰係数α1、β1、及びβ2を推計する。気象実績情報に融雪量がない場合は,日付t−1の積雪量と日付tのその他の気象情報を用いて,融雪量モデル4を用いて算出したMtを用いて,流入量モデル5を回帰分析し、回帰係数α1、β1、及びβ2を推計することもできる。流入量モデル推計部23は、推計した回帰係数α1、β1、及びβ2をパラメータ記憶部32に登録する。
The inflow amount
予測気温取得部24は、気温の予測値(以下、予測気温という。)を取得する。予測気温取得部24は、例えば、ユーザから予測気温の入力を受け付けてもよいし、気象庁や民間気象会社のコンピュータにアクセスして予測気温を取得するようにしてもよい。また、予測気温取得部24は、気象実績情報に基づいて気温の予測を行うようにしてもよい。この場合、例えば、一般的な気温の予測に用いられる統計モデルをモデル記憶部31に記憶しておき、予測気温取得部24がその統計モデルと気象実績情報とに基づいて回帰分析を行ってパラメータを推計し、推計したパラメータと気象実績情報とを統計モデルに適用して予測気温を算出することができる。
The predicted
予測降水量取得部25は、降水量の予測値(以下、予測降水量という。)を取得する。予測降水量取得部25は、例えば、ユーザから予測降水量の入力を受け付けてもよいし、気象庁や民間気象会社のコンピュータにアクセスして予測降水量を取得するようにしてもよい。また、予測降水量取得部25は、予測気温取得部24と同様に、気象実績情報に基づいて降水量の予測を行うようにしてもよい。
The predicted
流入量予測部26は、予測気温取得部24が取得した予測気温、予測降水量取得部25が取得した予測降水量、降雪気温推計部21が推計したパラメータ、気象実績情報、及び流入量モデル5を用いて、流入量の予測値(以下、予測流入量という。)を算出する。なお、予測流入量の算出処理の詳細については後述する。
The inflow
以下、本実施形態の流入量予測システムにおける処理の詳細について説明する。 Hereinafter, details of processing in the inflow amount prediction system of the present embodiment will be described.
==推計処理==
まず、降雪気温推計部21による降雪気温δの推計処理について説明する。図7は、降雪気温推計部21による降雪気温δの推計処理の流れを示す図である。
== Estimation processing ==
First, the process for estimating the snowfall temperature δ by the snowfall
まず、降雪気温推計部21は、気象実績データベース40から降雪量が0より大きい気象実績情報を取得する(S501)。これにより、図8に示すように、気象実績データベース40に記憶されている気象実績情報のうち、降雪量が0より大きいもののみが抽出される。
次に、降雪気温推計部21は、抽出した気象実績情報に基づき、モデル
降雪量=a×気温+b×降水量
を回帰分析して、回帰係数a及びbを推計する(S502)。
降雪気温推計部21は、気温がδ以下の場合には「a×気象実績情報の気温+b×気象実績情報の降水量」、気温がδより高い場合には「0」を推計降雪量として、推計降雪量と気象実績情報の降雪量との差を2乗した値が最小になるδを算出する(S503)。
First, the snowfall
Next, the snowfall
The snowfall
降雪気温推計部21は、例えば、所定範囲の気温を所定ステップごとに増加させたδについて、図9に示すように、各気象実績情報について、気温がδより高ければ、上記回帰係数aを気温に乗じた値と、上記か域係数bを降水量に乗じた値とを合計して推計降雪量として算出し、推計降雪量と気象実績情報の降雪量との差を2乗した値を誤差の2乗として算出していき、誤差の2乗が最も小さくなったものをδとして決定することができる。なお、上記誤差の2乗が最も小さくなるように上記δを決定する処理については、一般的な統計手法を利用することが可能である。
For example, as shown in FIG. 9, the snowfall
降雪気温推計部21は、上記のようにして決定したδを、パラメータ記憶部32に登録する(S504)。
以上のようにして融雪気温δが決定される。
The snowfall
The snow melting temperature δ is determined as described above.
==流入量の予測処理==
次に、流入量の予測について説明する。図10は、本実施形態の流入量予測システム10における流入量予測処理の流れを示す図である。なお、図10の処理では、日付tが予測対象となる日付であるものとしている。
== Inflow prediction process ==
Next, prediction of the inflow amount will be described. FIG. 10 is a diagram illustrating a flow of the inflow amount prediction process in the inflow
予測気温取得部24は、予測気温Ttを取得する(S601)。予測降水量取得部25は、予測降水量Ptを取得する(S602)。上述したように、予測気温取得部24及び予測降水量取得部25は、例えば、ユーザからデータの入力を受け付けたり、気象庁や民間気象会社のコンピュータにアクセスしてデータを取得することにより、TtやPtを取得することができる。
The predicted
流入量予測部26は、パラメータ記憶部32からδを読み出し(S603)、降水量モデル1にδ、Tt及びPtを適用してP1tを算出する(S604)。すなわち、予測気温Ttがδより大きければP1t=Ptとなり、Ttがδ以下であればP1t=0となる。
The inflow
流入量予測部26は、パラメータ記憶部32から、α1〜α3、β1、β2、μ1、μ2を読み出す(S605)。流入量予測部26は、気象実績データベース40から、前日t−1に対応する気象実績情報を読み出し、読み出した気象実績情報の積雪量をDt−1とし(S606)、読み出した気象実績情報の流入量をFt−1とする(S607)。
The inflow
流入量予測部26は、融雪量モデル4に、α2、Tt、μ2、Dt−1、α3、P1tを代入して、融雪量の予測値Mtを算出し(S608)、流入量モデル5に、α1、μ1、Ft−1、β1、P1t、β2、Mtを代入して流入増加量の予測値ΔFtを算出する(S609)。流入量予測部26は、Ft−1にΔFtを加算して、予測流入量Ftを算出する(S610)。
The inflow
以上のようにして、本実施形態の流入量予測システム10によれば、予測気温及び予測降水量と、気象実績情報とに基づいて、降水量及び融雪量を考慮した流入量の予測を行うことができる。降水がなかい場合にも融雪により流入量は増加するため、融雪量を考慮して流入量の予測を行うことで、予測の精度を向上することができる。
As described above, according to the inflow
また、本実施形態の流入量予測システム10によれば、融雪量は、降水量及び気温から算出することができる。降水量及び気温の予測は、気象予報の手法として様々な手法が存在し、容易に入手可能である。したがって、融雪量の予測が困難である場合でも、容易に取得可能な降水量や気温の予測値に基づいて融雪量の予測を行うことで、融雪量を考慮した流入量の予測を容易に実現することができる。
Moreover, according to the inflow
また、上記流入量モデル5では、第1項を均衡流入量と流入量との差であるものとして、均衡流入量を考慮に入れているので、単に流入量を説明変数とする場合に比べ、より精度の高い流入量予測を行うことができる。
Also, in the
なお、本実施形態の流入量予測システムでは、ダムなどの貯水施設に河川から流入する水量についての予測を行うものとしたが、河川を流れる水量の予測を行うシステムにも容易に適用することができる。この場合、河川の上流域における気温や降水量の予報値及び実績値などを取得及び記録するものとする。 In the inflow amount prediction system of the present embodiment, the amount of water flowing into the water storage facility such as a dam is predicted from the river. However, the inflow amount prediction system can be easily applied to a system that predicts the amount of water flowing through the river. it can. In this case, the forecast value and the actual value of the temperature and precipitation in the upstream area of the river shall be acquired and recorded.
また、本実施形態では、均衡流入量μ1及び融雪が始まる気温μ2については、予めパラメータ記憶部32に記憶されているものとしたが、これに限らず、過去の気象実績情報に基づいて当てはまりのよい値を推計するようにしてもよい。
In the present embodiment, the equilibrium inflow amount μ 1 and the temperature μ 2 at which snow melting starts are stored in the
また、貯水施設が降雪や積雪の少ない地域に存在する場合などには、流入量モデル5は融雪量Mを考慮しないようにしてもよい。この場合流入量モデル5は、次式で表される。
Further, when the water storage facility exists in an area where there is little snowfall or snowfall, the
また、本実施形態の各回帰モデルについて、誤差項に系列相関があるような場合には、降雪気温推計部21がPrais-Winstein変換やコクラン・オーカット法によりパラメータを推計するようにしてもよい。この場合、例えば、コクラン・オーカット法により流入量モデル5のパラメータを推計する場合、残差をεtとし、説明変数Xt*=Xt−ρXt−1として(すなわち、Ft*=Ft−ρFt−1、P1t*=P1t−ρP1t−1、Mt*=Mt−ρMt−1などとして)、流入量モデル5を式のようにすることができる。
上記式(5”)を変形すると次式となる。
In addition, for each regression model of the present embodiment, when there is a series correlation in the error term, the snow
When the above equation (5 ″) is modified, the following equation is obtained.
上記のようにして、残差εtに系列相関がなくなるn期前まで遡って式変形を行った流入量モデル5を求めることができる、その他のモデルについても同様に、誤差項に系列相関がなくなるようにパラメータを推計することができる。上記のようにして残差εtの系列相関がなくなるn期遅れのモデルをモデル記憶部31に記憶しておき、n期遅れのモデルを用いて流入量を予測することにより、残差項に系列相関がある場合にも、適切なパラメータの推計を行うことが可能となり、適切なパラメータを用いて精度よく流入量を予測することができる。
As described above, it is possible to obtain the
以上、本実施形態について説明したが、上記実施形態は本発明の理解を容易にするためのものであり、本発明を限定して解釈するためのものではない。本発明は、その趣旨を逸脱することなく、変更、改良され得ると共に、本発明にはその等価物も含まれる。 Although the present embodiment has been described above, the above embodiment is intended to facilitate understanding of the present invention and is not intended to limit the present invention. The present invention can be changed and improved without departing from the gist thereof, and the present invention includes equivalents thereof.
10 流入量予測システム
11 CPU
12 メモリ
13 記憶装置
14 入力装置
15 出力装置
21 降雪気温推計部
22 融雪量モデル推計部
23 流入量モデル推計部
31 モデル記憶部
32 パラメータ記憶部
40 気象実績データベース
10
DESCRIPTION OF
Claims (9)
所定期間ごとに、前記流入量の実績値、及び降雨量の実績値を記憶するデータベースと、
ある第1期間における前記降雨量、及び前記第1期間から所定期前の第2期間における前記流入量と所定の定数である均衡流入量との差を説明変数とし、前記第2期間から前記第1期間までの前記流入量の増加量を目的変数とする回帰モデルを記憶するモデル記憶部と、
前記データベースに記憶されている前記流入量の実績値及び前記降雨量の実績値、並びに前記モデル記憶部に記憶されている前記回帰モデルに基づいて回帰分析を行い、前記説明変数のそれぞれについての回帰係数を推計する流入量モデル推計部と、
予測対象の期間である予測期間における前記降雨量の予測値を取得する予測降雨量取得部と、
前記予測期間から前記所定期前の期間における前記流入量の実績値を前記データベースから取得し、前記降雨量の予測値、前記流入量の実績値、及び前記各回帰係数を前記回帰モデルに適用して、前記予測期間における前記増加量の予測値を算出する増加量予測部と、
前記流入量の実績値に前記増加量の予測値を加算して、前記予測期間における前記流入量の予測値を算出する流入量予測部と、
を備えることを特徴とする流入量予測システム。 A system for predicting the amount of water flowing into a water storage facility,
A database for storing the actual value of the inflow and the actual value of the rainfall for each predetermined period;
The difference between the rainfall amount in a certain first period and the inflow amount in the second period before the first period from the first period and the equilibrium inflow amount that is a predetermined constant is used as an explanatory variable, and the second period from the second period to the second period. A model storage unit for storing a regression model whose objective variable is an increase amount of the inflow until one period;
Based on the actual value of the inflow and the actual value of the rainfall stored in the database, and the regression model stored in the model storage unit, regression analysis is performed for each of the explanatory variables. An inflow model estimation unit for estimating a coefficient;
A predicted rainfall amount acquisition unit that acquires a predicted value of the rainfall amount in a prediction period that is a prediction target period;
The actual value of the inflow from the prediction period to the period before the predetermined period is acquired from the database, and the predicted value of the rainfall, the actual value of the inflow, and the regression coefficients are applied to the regression model. An increase amount prediction unit for calculating a predicted value of the increase amount in the prediction period;
An inflow amount prediction unit that calculates the predicted value of the inflow amount in the prediction period by adding the predicted value of the increase amount to the actual value of the inflow amount;
An inflow rate prediction system comprising:
前記データベースは、前記所定期間ごとに、融雪量の実績値をさらに記憶し、
前記モデル記憶部に記憶される前記回帰モデルには、説明変数として、さらに前記第1期間における前記融雪量が含まれ、
前記予測期間における前記融雪量の予測値を取得する予測融雪量取得部を備え、
前記流入量モデル推計部は、前記データベースに記憶されている前記流入量の実績値、前記降雨量の実績値、及び前記融雪量の実績値、並びに前記モデル記憶部に記憶されている前記回帰モデルに基づいて回帰分析を行い、前記回帰係数を推計すること、
を特徴とする流入量予測システム。 The inflow amount prediction system according to claim 1,
The database further stores the actual value of the amount of snow melt for each predetermined period,
The regression model stored in the model storage unit further includes the amount of snow melt in the first period as an explanatory variable,
A predicted snowmelt amount acquisition unit for acquiring a predicted value of the snowmelt amount in the prediction period;
The inflow amount model estimation unit includes the actual value of the inflow amount stored in the database, the actual value of the rainfall amount, the actual value of the snowmelt amount, and the regression model stored in the model storage unit. Performing a regression analysis based on the above and estimating the regression coefficient,
An inflow prediction system characterized by
前記データベースはさらに、気温の実績値及び積雪量の実績値を記憶し、
前記モデル記憶部は、前記回帰モデルである流入量モデルに加えて、前記第1期間における気温から融雪の始まる気温を示す定数を引いた値又はゼロの大きい方の値に、前記第2期間における前記積雪量を乗じた値と、前記第1期間における前記降雨量とを説明変数とし、前記第1期間における前記融雪量を目的変数とする回帰モデルである融雪量モデルを記憶し、
前記データベースに記憶されている前記気温の実績値、前記積雪量の実績値及び前記降雨量の実績値、並びに前記モデル記憶部に記憶されている前記融雪量モデルに基づいて回帰分析を行い、前記融雪量モデルに係る前記説明変数のそれぞれについての回帰係数を推計する融雪量モデル推計部と、
前記予測期間における前記気温の予測値を取得する予測気温取得部と、
を備え、
前記予測融雪量取得部は、前記予測期間の前記所定期前の期間における前記積雪量の実績値を前記データベースから取得し、前記気温の予測値、前記取得した積雪量の実績値、前記降雨量の予測値、及び前記融雪量モデルに係る前記各説明変数を、前記融雪量モデルに適用して、前記予測期間における前記融雪量の予測値を算出すること、
を特徴とする流入量予測システム。 The inflow amount prediction system according to claim 2,
The database further stores an actual temperature value and an actual snow amount value,
In addition to the inflow model that is the regression model, the model storage unit subtracts a constant indicating the temperature at which snow melting starts from the temperature in the first period or a larger value of zero in the second period. Storing a snowmelt amount model which is a regression model having the value obtained by multiplying the amount of snow and the amount of rainfall in the first period as explanatory variables, and the amount of snowmelt in the first period as an objective variable;
Performing regression analysis based on the actual temperature value stored in the database, the actual snow accumulation value, the actual rainfall value, and the snowmelt model stored in the model storage unit, A snowmelt amount model estimation unit for estimating a regression coefficient for each of the explanatory variables related to the snowmelt amount model;
A predicted temperature acquisition unit that acquires a predicted value of the temperature in the prediction period;
With
The predicted snowmelt amount acquisition unit acquires the actual value of the snow cover amount in the period before the predetermined period of the prediction period from the database, the predicted value of the temperature, the actual value of the acquired snow cover amount, the rainfall amount Applying the predicted value of each and the explanatory variables related to the snowmelt amount model to the snowmelt amount model to calculate the predicted value of the snowmelt amount in the prediction period;
An inflow prediction system characterized by
前記データベースはさらに、降雪量の実績値を記憶し、
前記融雪モデル推計部は、前記積雪量モデルに含まれる前記積雪量の説明変数に、式
前記第1期間における積雪量=max(前記第1期間における前記降雪量+前記第2期間における前記積雪量−前記第1期間における前記融雪量,0)
を代入したモデルを回帰分析することにより、前記融雪量モデルに係る前記説明変数のそれぞれについての前記回帰係数を推計すること、
を特徴とする流入量予測システム。 The inflow amount prediction system according to claim 3,
The database further stores the actual amount of snowfall,
The snow melting model estimation unit uses an expression of the snow accumulation amount included in the snow accumulation model as an expression variable: snow accumulation amount in the first period = max (the snowfall amount in the first period + the snow accumulation amount in the second period) -The amount of snow melt in the first period, 0)
Estimating the regression coefficient for each of the explanatory variables related to the snowmelt amount model by performing regression analysis of the model into which
An inflow prediction system characterized by
前記データベースには、前記降雨量の実績値に代えて、降水量の実績値が記憶されており、
前記予測期間における前記降水量の予測値を取得する予測降水量取得部と、
前記予測期間における前記気温の予測値を取得する予測気温取得部と、
を備え、
前記予測降雨量取得部は、前記気温の予測値が所定の気温を超えている場合に前記降水量を前記降雨量とし、前記気温の予測値が前記所定の気温以下である場合には前記降雨量を0とすること、
を特徴とする流入量予測システム。 The inflow amount prediction system according to claim 1,
In the database, instead of the actual rainfall value, the actual rainfall value is stored,
A predicted precipitation acquisition unit for acquiring a predicted value of the precipitation in the prediction period;
A predicted temperature acquisition unit that acquires a predicted value of the temperature in the prediction period;
With
The predicted rainfall amount acquisition unit sets the precipitation amount as the rainfall when the predicted value of the temperature exceeds a predetermined temperature, and the rainfall when the predicted value of the temperature is equal to or lower than the predetermined temperature. Set the amount to 0,
An inflow prediction system characterized by
前記データベースはさらに、気温の実績値及び降雪量の実績値を記憶し、
前記気温、及び前記降水量を説明変数とし、前記降雪量を目的変数とした回帰モデルである降雪モデル、並びに前記データベースに記憶されている前記気温の実績値、前記降水量の実績値、及び前記降雪量の実績値に基づいて回帰分析を行い、前記説明変数のそれぞれについての回帰係数を推計する降雪モデル推計部と、
前記データベースに記憶されている前記気温の実績値のそれぞれについて、前記気温がδ以下の場合には、前記推計した回帰係数と前記気温の実績値及び前記降水量の実績値とを前記降雪モデルに適用して算出した値を推計降雪量とし、前記気温がδより高い場合には0を前記推計降雪量としたときに、前記降雪量の実績値と前記推計降雪量との誤差が最も小さくなる前記δを前記所定の気温として決定する降雪気温推計部と、
を備えることを特徴とする流入量予測システム。 The inflow amount prediction system according to claim 5,
The database further stores an actual temperature value and an actual snowfall value,
A snowfall model that is a regression model with the temperature and the precipitation as explanatory variables, and the snowfall as an objective variable, and the actual value of the temperature stored in the database, the actual value of the precipitation, and the A snowfall model estimation unit that performs regression analysis based on the actual amount of snowfall, and estimates a regression coefficient for each of the explanatory variables;
For each of the actual temperature values stored in the database, when the temperature is δ or less, the estimated regression coefficient, the actual temperature value, and the actual precipitation value are stored in the snowfall model. When the estimated snowfall amount is a value calculated by application, and 0 is the estimated snowfall amount when the temperature is higher than δ, the error between the actual value of the snowfall amount and the estimated snowfall amount is the smallest. A snow temperature estimating unit for determining δ as the predetermined temperature;
An inflow rate prediction system comprising:
所定期間ごとに、前記流入量の実績値、降雨量の実績値、気温の実績値、及び積雪量の実績値を記憶するデータベースと、
ある第1期間における前記降雨量、前記第1期間における融雪量、及び前記第1期間から所定期前の第2期間における前記流入量と所定の定数である均衡流入量との差を説明変数とし、前記第2期間から前記第1期間までの前記流入量の増加量を目的変数とする回帰モデルである流入量モデル、及び、前記第1期間における気温から融雪の始まる気温を示す定数を引いた値又はゼロの大きい方の値に、前記第2期間における前記積雪量を乗じた値と、前記第1期間における前記降雨量とを説明変数とし、前記第1期間における前記融雪量を目的変数とする回帰モデルである融雪量モデルを記憶するモデル記憶部と、
前記データベースに記憶されている前記気温の実績値、前記積雪量の実績値及び前記降雨量の実績値、並びに前記モデル記憶部に記憶されている前記融雪量モデルに基づいて回帰分析を行い、前記融雪量モデルに係る前記説明変数のそれぞれについての回帰係数を推計する融雪量モデル推計部と、
前記データベースに記憶されている前記気温の実績値、前記積雪量の実績値、前記降雨量の実績値、及び前記推計した回帰係数を前記融雪量モデルに適用して、前記所定期間ごとの前記融雪量の理論値を算出し、前記データベースに記憶されている前記流入量の実績値、及び前記降雨量の実績値、前記算出した融雪量の理論値、並びに前記モデル記憶部に記憶されている前記流入量モデルに基づいて回帰分析を行い、前記流入量モデルに係る前記説明変数のそれぞれについての回帰係数を推計する流入量モデル推計部と、
予測対象の期間である予測期間における前記降雨量の予測値を取得する予測降雨量取得部と、
前記予測期間における前記気温の予測値を取得する予測気温取得部と、
前記予測期間の前記所定期前の期間における前記積雪量の実績値を前記データベースから取得し、前記気温の予測値、前記取得した積雪量の実績値、前記降雨量の予測値、及び前記融雪量モデルに係る前記各説明変数を、前記融雪量モデルに適用して、前記予測期間における前記融雪量の予測値を算出する予測融雪量取得部と、
前記予測期間から前記所定期前の期間における前記流入量の実績値を前記データベースから取得し、前記降雨量の予測値、前記融雪量の予測値、前記流入量の実績値、及び前記各回帰係数を前記流入量モデルに適用して、前記予測期間における前記増加量の予測値を算出する増加量予測部と、
前記流入量の実績値に前記増加量の予測値を加算して、前記予測期間における前記流入量の予測値を算出する流入量予測部と、
を備えることを特徴とする流入量予測システム。 A system for predicting the amount of water flowing into a water storage facility,
For each predetermined period, a database that stores the actual value of the inflow, the actual value of the rainfall, the actual value of the temperature, and the actual value of the amount of snow;
The difference between the rainfall amount in a certain first period, the snowmelt amount in the first period, and the inflow amount in the second period before the predetermined period from the first period and the equilibrium inflow amount that is a predetermined constant is used as an explanatory variable. , An inflow model that is a regression model with an increase in the inflow from the second period to the first period as an objective variable, and a constant indicating the temperature at which snow melting starts is subtracted from the temperature in the first period. A value obtained by multiplying the larger value or zero by the amount of snow in the second period and the amount of rainfall in the first period are explanatory variables, and the amount of snow melt in the first period is an objective variable. A model storage unit that stores a snowmelt amount model that is a regression model
Performing regression analysis based on the actual temperature value stored in the database, the actual snow accumulation value, the actual rainfall value, and the snowmelt model stored in the model storage unit, A snowmelt amount model estimation unit for estimating a regression coefficient for each of the explanatory variables related to the snowmelt amount model;
Applying the actual temperature value, the actual snow accumulation value, the actual rainfall value, and the estimated regression coefficient stored in the database to the snow melting model, the snow melting for each predetermined period. The theoretical value of the amount is calculated, the actual value of the inflow amount stored in the database, the actual value of the rainfall amount, the theoretical value of the calculated snow melting amount, and the model storage unit An inflow model estimation unit that performs a regression analysis based on the inflow model, and estimates a regression coefficient for each of the explanatory variables related to the inflow model;
A predicted rainfall amount acquisition unit that acquires a predicted value of the rainfall amount in a prediction period that is a prediction target period;
A predicted temperature acquisition unit that acquires a predicted value of the temperature in the prediction period;
The actual value of the snow cover amount in the period before the predetermined period of the prediction period is acquired from the database, the predicted value of the temperature, the acquired actual value of the snow cover amount, the predicted value of the rainfall amount, and the snow melting amount Applying each of the explanatory variables related to the model to the snowmelt amount model, and calculating a predicted value of the snowmelt amount in the prediction period;
The actual value of the inflow amount in the period from the prediction period to the predetermined period is acquired from the database, the predicted value of the rainfall amount, the predicted value of the snowmelt amount, the actual value of the inflow amount, and each regression coefficient Applying to the inflow amount model, an increase amount prediction unit that calculates a predicted value of the increase amount in the prediction period;
An inflow amount prediction unit that calculates the predicted value of the inflow amount in the prediction period by adding the predicted value of the increase amount to the actual value of the inflow amount;
An inflow rate prediction system comprising:
コンピュータが、
所定期間ごとに、前記流入量の実績値、及び降雨量の実績値をメモリに記憶し、
ある第1期間における前記降雨量、及び前記第1期間から所定期前の第2期間における前記流入量と所定の定数である均衡流入量との差を説明変数とし、前記第2期間から前記第1期間までの前記流入量の増加量を目的変数とする回帰モデルを前記メモリに記憶し、
前記メモリに記憶されている前記流入量の実績値及び前記降雨量の実績値、並びに前記モデル記憶部に記憶されている前記回帰モデルに基づいて回帰分析を行い、前記説明変数のそれぞれについての回帰係数を推計し、
予測対象の期間である予測期間における前記降雨量の予測値を取得し、
前記予測期間から前記所定期前の期間における前記流入量の実績値を前記メモリから取得し、前記降雨量の予測値、前記流入量の実績値、及び前記各回帰係数を前記回帰モデルに適用して、前記予測期間における前記増加量の予測値を算出し、
前記流入量の実績値に前記増加量の予測値を加算して、前記予測期間における前記流入量の予測値を算出すること、
を特徴とする流入量予測方法。 A method for predicting the amount of water flowing into a water storage facility,
Computer
For each predetermined period, the actual value of the inflow and the actual value of the rainfall are stored in a memory,
The difference between the rainfall amount in a certain first period and the inflow amount in the second period before the first period from the first period and the equilibrium inflow amount that is a predetermined constant is used as an explanatory variable, and the second period from the second period to the second period. Storing in the memory a regression model whose objective variable is the amount of increase in inflow up to one period;
The regression analysis is performed based on the actual value of the inflow amount and the actual value of the rainfall amount stored in the memory, and the regression model stored in the model storage unit, and the regression for each of the explanatory variables. Estimate the coefficient,
Get the predicted value of the rainfall in the forecast period, which is the forecast target period,
The actual value of the inflow from the prediction period to the period before the predetermined period is acquired from the memory, and the predicted value of the rainfall, the actual value of the inflow, and the regression coefficients are applied to the regression model. Calculating a predicted value of the increase amount in the prediction period,
Adding the predicted value of the increase amount to the actual value of the inflow amount to calculate the predicted value of the inflow amount in the prediction period;
An inflow rate prediction method characterized by
コンピュータに、
所定期間ごとに、前記流入量の実績値、及び降雨量の実績値をメモリに記憶するステップと、
ある第1期間における前記降雨量及び前記第1期間から所定期前の第2期間における前記流入量と所定の定数である均衡流入量との差を説明変数とし、前記第2期間から前記第1期間までの前記流入量の増加量を目的変数とする回帰モデルを前記メモリに記憶するステップと、
前記メモリに記憶されている前記流入量の実績値、及び前記降雨量の実績値、並びに前記モデル記憶部に記憶されている前記回帰モデルに基づいて回帰分析を行い、前記説明変数のそれぞれについての回帰係数を推計するステップと、
予測対象の期間である予測期間における前記降雨量の予測値を取得するステップと、
前記予測期間から前記所定期前の期間における前記流入量の実績値を前記メモリから取得し、前記降雨量の予測値、前記融雪量の予測値、前記流入量の実績値、及び前記各回帰係数を前記回帰モデルに適用して、前記予測期間における前記増加量の予測値を算出するステップと、
前記流入量の実績値に前記増加量の予測値を加算して、前記予測期間における前記流入量の予測値を算出するステップと、
を実行させるためのプログラム。 A program for predicting the amount of water flowing into a water storage facility,
On the computer,
Storing the actual value of the inflow and the actual value of the rainfall in a memory for each predetermined period;
The difference between the rainfall amount in a certain first period and the inflow amount in the second period before the predetermined period from the first period and the equilibrium inflow amount that is a predetermined constant is used as an explanatory variable, and the first period to the first period Storing, in the memory, a regression model whose objective variable is an increase in the amount of inflow until a period;
The regression analysis is performed based on the actual value of the inflow amount stored in the memory, the actual value of the rainfall amount, and the regression model stored in the model storage unit. Estimating a regression coefficient;
Obtaining a predicted value of the rainfall amount in a prediction period that is a period to be predicted;
The actual value of the inflow amount in the period before the predetermined period from the prediction period is acquired from the memory, the predicted value of the rainfall amount, the predicted value of the snowmelt amount, the actual value of the inflow amount, and each regression coefficient Applying to the regression model to calculate a predicted value of the increase in the prediction period;
Adding the predicted value of the increase amount to the actual value of the inflow amount to calculate the predicted value of the inflow amount in the prediction period;
A program for running
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008009568A JP4827861B2 (en) | 2008-01-18 | 2008-01-18 | Inflow prediction system, inflow prediction method and program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008009568A JP4827861B2 (en) | 2008-01-18 | 2008-01-18 | Inflow prediction system, inflow prediction method and program |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2009167751A true JP2009167751A (en) | 2009-07-30 |
JP4827861B2 JP4827861B2 (en) | 2011-11-30 |
Family
ID=40969236
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2008009568A Expired - Fee Related JP4827861B2 (en) | 2008-01-18 | 2008-01-18 | Inflow prediction system, inflow prediction method and program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4827861B2 (en) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011065379A (en) * | 2009-09-16 | 2011-03-31 | Chugoku Electric Power Co Inc:The | Operation support system, operation support method, and program for water storage facility |
JP2012012907A (en) * | 2010-07-05 | 2012-01-19 | Chubu Electric Power Co Inc | Dam measurement data evaluation method and dam measurement data evaluation system |
JP5014522B1 (en) * | 2011-06-28 | 2012-08-29 | 中国電力株式会社 | Inflow amount prediction apparatus, inflow amount prediction method, and program |
JP2013079530A (en) * | 2011-10-04 | 2013-05-02 | Chugoku Electric Power Co Inc:The | Inflow prediction device, inflow prediction method, and program |
WO2014178116A1 (en) * | 2013-04-30 | 2014-11-06 | 中国電力株式会社 | Flow quantity forecasting device and flow quantity forecasting system |
JP2015049177A (en) * | 2013-09-03 | 2015-03-16 | 中国電力株式会社 | Inflow prediction device, inflow prediction method, water level prediction device and program |
JP2015146065A (en) * | 2014-01-31 | 2015-08-13 | 中国電力株式会社 | Forecasting system and forecasting method |
JP2015525809A (en) * | 2012-07-03 | 2015-09-07 | スリーエム イノベイティブ プロパティズ カンパニー | Method for manufacturing structured hybrid adhesive article |
JP2016130428A (en) * | 2015-01-14 | 2016-07-21 | 中国電力株式会社 | System and method for predicting inflow water volume into reservoir |
JP2019219978A (en) * | 2018-06-21 | 2019-12-26 | 中国電力株式会社 | Power output prediction device, power output prediction system and power output prediction method |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08219828A (en) * | 1995-02-08 | 1996-08-30 | Hitachi Ltd | Method and device for estimating inflow to dam |
JPH11323885A (en) * | 1998-05-11 | 1999-11-26 | Mitsubishi Heavy Ind Ltd | River water level preestimating system |
JP2007205001A (en) * | 2006-02-01 | 2007-08-16 | Fuji Electric Systems Co Ltd | Discharge forecasting apparatus |
-
2008
- 2008-01-18 JP JP2008009568A patent/JP4827861B2/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08219828A (en) * | 1995-02-08 | 1996-08-30 | Hitachi Ltd | Method and device for estimating inflow to dam |
JPH11323885A (en) * | 1998-05-11 | 1999-11-26 | Mitsubishi Heavy Ind Ltd | River water level preestimating system |
JP2007205001A (en) * | 2006-02-01 | 2007-08-16 | Fuji Electric Systems Co Ltd | Discharge forecasting apparatus |
Non-Patent Citations (3)
Title |
---|
山田富士宏, 他8名: "逓減時定数を用いた降雨後における発電用ダム流入量予測", 電気学会電力技術研究会資料, vol. PE-08巻, 116-120号, JPN6011005685, 2008, JP, pages 1 - 6, ISSN: 0002011650 * |
市山健二: "三国川ダムにおける融雪出水予測システムの構築", ダム技術, vol. 第258巻, JPN6011005687, 2008, JP, pages 52 - 58, ISSN: 0002011652 * |
玉腰尚久: "オンラインによる出水予測システムの開発 飛騨川出水予測システムの開発", 中部電力株式会社技術開発ニュース, vol. 第79巻, JPN6011005686, 1999, JP, pages 13 - 14, ISSN: 0002011651 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011065379A (en) * | 2009-09-16 | 2011-03-31 | Chugoku Electric Power Co Inc:The | Operation support system, operation support method, and program for water storage facility |
JP2012012907A (en) * | 2010-07-05 | 2012-01-19 | Chubu Electric Power Co Inc | Dam measurement data evaluation method and dam measurement data evaluation system |
JP5014522B1 (en) * | 2011-06-28 | 2012-08-29 | 中国電力株式会社 | Inflow amount prediction apparatus, inflow amount prediction method, and program |
WO2013001600A1 (en) * | 2011-06-28 | 2013-01-03 | 中国電力株式会社 | Inflow prediction device, and inflow prediction method and program |
JP2013079530A (en) * | 2011-10-04 | 2013-05-02 | Chugoku Electric Power Co Inc:The | Inflow prediction device, inflow prediction method, and program |
JP2015525809A (en) * | 2012-07-03 | 2015-09-07 | スリーエム イノベイティブ プロパティズ カンパニー | Method for manufacturing structured hybrid adhesive article |
WO2014178116A1 (en) * | 2013-04-30 | 2014-11-06 | 中国電力株式会社 | Flow quantity forecasting device and flow quantity forecasting system |
JPWO2014178116A1 (en) * | 2013-04-30 | 2017-02-23 | 中国電力株式会社 | Flow prediction device and flow prediction system |
JP2015049177A (en) * | 2013-09-03 | 2015-03-16 | 中国電力株式会社 | Inflow prediction device, inflow prediction method, water level prediction device and program |
JP2015146065A (en) * | 2014-01-31 | 2015-08-13 | 中国電力株式会社 | Forecasting system and forecasting method |
JP2016130428A (en) * | 2015-01-14 | 2016-07-21 | 中国電力株式会社 | System and method for predicting inflow water volume into reservoir |
JP2019219978A (en) * | 2018-06-21 | 2019-12-26 | 中国電力株式会社 | Power output prediction device, power output prediction system and power output prediction method |
JP7139711B2 (en) | 2018-06-21 | 2022-09-21 | 中国電力株式会社 | Power generation output prediction device, power generation output prediction system, power generation output prediction method |
Also Published As
Publication number | Publication date |
---|---|
JP4827861B2 (en) | 2011-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4827861B2 (en) | Inflow prediction system, inflow prediction method and program | |
JP5823459B2 (en) | Inflow amount prediction device, inflow amount prediction method, water level prediction device, and program | |
JP5014213B2 (en) | Water storage facility operation support system, operation support method and program | |
Arguez et al. | The definition of the standard WMO climate normal: The key to deriving alternative climate normals | |
Li et al. | Application of Bayesian model averaging in modeling long-term wind speed distributions | |
JP6880841B2 (en) | Photovoltaic power generation output prediction device considering snow cover | |
JP6794878B2 (en) | Photovoltaic power generation output prediction device considering the remaining snow | |
JP2008241433A (en) | Observation data assimilation method | |
KR20200080496A (en) | System and method for predicting watershed runoff information | |
JP2006092058A (en) | Flow rate estimation device | |
JP5034021B2 (en) | Wave estimation method, system and program | |
CN109299521B (en) | Method and device for determining runoff total control rate | |
JP5014522B1 (en) | Inflow amount prediction apparatus, inflow amount prediction method, and program | |
Bishop et al. | Bayesian model averaging’s problematic treatment of extreme weather and a paradigm shift that fixes it | |
JP4413202B2 (en) | Dam basin snow amount estimation device, dam inflow water amount estimation device and program | |
JP7067306B2 (en) | Snow cover prediction device, snow cover prediction method | |
KR20180060286A (en) | Calibration method and system to improve accuracy of short-term precipitation prediction | |
JP6191704B2 (en) | Inflow amount prediction apparatus, inflow amount prediction method, and inflow amount prediction program | |
WO2019207622A1 (en) | Power demand prediction device, power demand prediction method, and program therefor | |
JP5383767B2 (en) | Inflow amount prediction apparatus, inflow amount prediction method, and program | |
JP5425985B2 (en) | Water storage facility operation support system, operation support method and program | |
JP2018080509A (en) | Water level estimating program, water level estimating method, and water level estimating device | |
JP2008184782A (en) | Rainwater inflow estimating device and rainwater inflow estimating method | |
JP2011180025A (en) | Action prediction device, method and program | |
JP2006045892A (en) | Apparatus for predicting amount of water supply |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20090703 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20100621 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20110208 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20110906 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20110913 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140922 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4827861 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140922 Year of fee payment: 3 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |