JP2014234675A - Flow rate prediction device, flow rate prediction method, flow rate prediction program and flow rate prediction system - Google Patents

Flow rate prediction device, flow rate prediction method, flow rate prediction program and flow rate prediction system Download PDF

Info

Publication number
JP2014234675A
JP2014234675A JP2013117920A JP2013117920A JP2014234675A JP 2014234675 A JP2014234675 A JP 2014234675A JP 2013117920 A JP2013117920 A JP 2013117920A JP 2013117920 A JP2013117920 A JP 2013117920A JP 2014234675 A JP2014234675 A JP 2014234675A
Authority
JP
Japan
Prior art keywords
flow rate
rainfall intensity
flow
intensity value
value
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
Application number
JP2013117920A
Other languages
Japanese (ja)
Other versions
JP6222689B2 (en
Inventor
由紀夫 平岡
Yukio Hiraoka
由紀夫 平岡
眞哉 近藤
Masaya Kondo
眞哉 近藤
理 山中
Osamu Yamanaka
理 山中
直人 吉澤
Naoto Yoshizawa
直人 吉澤
孝典 石井
Takanori Ishii
孝典 石井
浩嗣 山本
Koji Yamamoto
浩嗣 山本
賢治 梅田
Kenji Umeda
賢治 梅田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toshiba Corp
Original Assignee
Toshiba Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Toshiba Corp filed Critical Toshiba Corp
Priority to JP2013117920A priority Critical patent/JP6222689B2/en
Publication of JP2014234675A publication Critical patent/JP2014234675A/en
Application granted granted Critical
Publication of JP6222689B2 publication Critical patent/JP6222689B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

PROBLEM TO BE SOLVED: To provide a flow rate prediction device, a flow rate prediction method, a flow rate prediction program and a flow rate prediction system capable of accurately predicting a flow rate of rainwater through a simple method with a reduced quantity of engineering work.SOLUTION: A flow rate prediction device comprises: a precipitation intensity acquisition section; a flow rate acquisition section; and a flow rate estimation section. The precipitation intensity acquisition section acquires precipitation intensity on a precipitation event at predetermined intervals in each of a plurality of partitioned regions partitioning an object region. The flow rate acquisition section acquires a flow rate of rainwater flowing out of the object region in synchronization with an acquisition timing of the precipitation intensity acquisition section. The flow rate estimation section estimates the flow rate of the rainwater flowing out of the object region at a future time on the basis of the precipitation intensity acquired by the precipitation intensity acquisition section and the flow rate acquired by the flow rate acquisition section.

Description

本発明の実施形態は、流量予測装置、流量予測方法、流量予測プログラム、および流量予測システムに関する。   Embodiments described herein relate generally to a flow prediction device, a flow prediction method, a flow prediction program, and a flow prediction system.

従来、降雨レーダによる降雨状況の観測結果を用いて、雨水処理場への雨水流入量を予測する技術が知られている(例えば、特許文献1参照)。   2. Description of the Related Art Conventionally, a technique for predicting the amount of rainwater inflow into a rainwater treatment plant using the observation result of a rainfall situation by a rain radar is known (see, for example, Patent Document 1).

特開2008−184783号公報JP 2008-184783 A

しかしながら、従来の雨水流入量の予測技術では、予測精度が十分ではなく、この予測精度を高めるためには、詳細なモデリングを行う等の多大なエンジニアリング作業を要していた。
そこで、本発明が解決しようとする課題は、エンジニアリング作業量を少なく抑えた簡便な方法によって雨水流量の予測を精度よく行うことができる、流量予測装置、流量予測方法、流量予測プログラム、および流量予測システムを提供することである。
However, the conventional rainwater inflow prediction technology does not have sufficient prediction accuracy, and enormous engineering work such as detailed modeling is required to improve the prediction accuracy.
Therefore, the problem to be solved by the present invention is that a flow rate prediction device, a flow rate prediction method, a flow rate prediction program, and a flow rate prediction that can accurately predict a rainwater flow rate by a simple method with a small amount of engineering work. Is to provide a system.

上記の課題を解決するため、本発明の実施形態である流量予測装置は、降雨強度値取得部と、流量値取得部と、流量値推定部とを備える。降雨強度値取得部は、対象領域を区分した複数の区分領域それぞれにおいて、降雨イベントにおける降雨強度値をあらかじめ一定時間ごとに取得する。流量値取得部は、前記対象領域から流出する雨水の流量値を、前記降雨強度値取得部による取得のタイミングに同期させて取得する。流量値推定部は、前記降雨強度値取得部が取得した降雨強度値と前記流量値取得部が取得した流量値とに基づいて、前記対象領域から流出する雨水の未来時刻における流量値を推定する。   In order to solve the above problems, a flow rate prediction apparatus according to an embodiment of the present invention includes a rainfall intensity value acquisition unit, a flow rate value acquisition unit, and a flow rate value estimation unit. The rainfall intensity value acquisition unit acquires the rainfall intensity value in the rain event for each predetermined time in advance in each of the plurality of divided areas obtained by dividing the target area. The flow value acquisition unit acquires the flow value of rainwater flowing out from the target area in synchronization with the acquisition timing by the rainfall intensity value acquisition unit. The flow rate estimation unit estimates a flow rate value at a future time of rainwater flowing out from the target region based on the rainfall intensity value acquired by the rainfall intensity value acquisition unit and the flow rate value acquired by the flow rate value acquisition unit. .

第1実施形態である流量予測装置を適用した雨水排水管理システムの全体構成を模式的に示す図である。It is a figure which shows typically the whole structure of the rainwater drainage management system to which the flow volume prediction apparatus which is 1st Embodiment is applied. 同実施形態において、流量予測装置の機能構成を示すブロック図である。In the embodiment, it is a block diagram which shows the function structure of a flow volume prediction apparatus. 同実施形態において、対象領域を区分した複数の区分領域それぞれに対応する降雨強度値と、当該対象領域から流出する雨水の流量値とを模式的に表した図である。In the same embodiment, it is the figure which represented typically the rainfall intensity value corresponding to each of the some division area which divided the object area | region, and the flow volume value of the rainwater which flows out out of the said object area | region. 同実施形態において、流量予測装置が取り込む降雨イベントデータの具体例を示すデータ構成図である。In the embodiment, it is a data block diagram which shows the specific example of the rain event data which a flow volume prediction apparatus takes in. 同実施形態において、流量予測装置の処理を示すフローチャートである。In the same embodiment, it is a flowchart which shows the process of a flow volume prediction apparatus. 同実施形態において、流量値推定部に記憶された係数のうち、k番目の区分領域(メッシュk)における各係数の値を示したグラフである。In the same embodiment, it is the graph which showed the value of each coefficient in the kth division field (mesh k) among the coefficients memorized by the flow rate estimation part. 同実施形態において、流量予測装置が生成した二次元データ上で、同一または同等の遅れ時間を対応付けることによって得た等到達時間曲線を表した図である。In the same embodiment, it is the figure showing the equal arrival time curve obtained by matching the same or equivalent delay time on the two-dimensional data which the flow volume prediction apparatus produced | generated.

以下、実施形態について、図面を参照して詳細に説明する。
[第1の実施形態]
図1は、第1実施形態である流量予測装置を適用した雨水排水管理システムの全体構成を模式的に示す図である。同図に示すように、雨水排水の対象となる地上の対象領域には、雨水管10(集水溝も含む。)が敷設されている。ただし、同図における雨水管10は、簡略化して表されている。雨水管10の排出口側には雨水排水施設80が設けられている。雨水排水施設80には、排水ポンプ20が雨水管10の排水口を接続して設置されている。よって、対象領域に降った雨は、その地表から雨水管10に流入または浸入して雨水排水施設80に集められる。雨水排水施設80の排水ポンプ20は、集水した雨水を、例えば河川や海に排出する。
Hereinafter, embodiments will be described in detail with reference to the drawings.
[First Embodiment]
FIG. 1 is a diagram schematically illustrating an overall configuration of a rainwater drainage management system to which a flow rate prediction apparatus according to the first embodiment is applied. As shown in the figure, a rainwater pipe 10 (including a water collecting groove) is laid in a target area on the ground that is a target of rainwater drainage. However, the rainwater pipe 10 in the figure is represented in a simplified manner. A rainwater drainage facility 80 is provided on the drain outlet side of the rainwater pipe 10. In the rainwater drainage facility 80, a drainage pump 20 is installed with a drainage port of the rainwater pipe 10 connected thereto. Therefore, the rain that has fallen on the target area flows into or enters the rainwater pipe 10 from the ground surface and is collected in the rainwater drainage facility 80. The drainage pump 20 of the rainwater drainage facility 80 discharges the collected rainwater to, for example, a river or the sea.

図1において、雨水排水管理システムは、降雨レーダ30と、信号処理装置40と、配信装置50と、流量予測装置60と、流量測定装置70とを含んでいる。   In FIG. 1, the stormwater drainage management system includes a rainfall radar 30, a signal processing device 40, a distribution device 50, a flow rate predicting device 60, and a flow rate measuring device 70.

降雨レーダ30は、少なくとも上記の対象領域を観測範囲に含めて降雨の状況を観測する降雨観測装置である。具体的に、降雨レーダ30は、空中に向けて電波を発信し、雨滴等から反射される電波を受信することによって降雨の状況を示す信号を生成する。降雨レーダ30は、生成した信号を信号処理装置40に供給する。   The rain radar 30 is a rain observation apparatus that observes the state of rainfall by including at least the above-mentioned target area in the observation range. Specifically, the rainfall radar 30 generates a signal indicating a rain condition by transmitting a radio wave toward the air and receiving a radio wave reflected from a raindrop or the like. The rain radar 30 supplies the generated signal to the signal processing device 40.

信号処理装置40は、降雨レーダ30が供給する信号を受信し、この信号を信号処理することによって、対象領域を区分した複数の区分領域(メッシュ)それぞれにおける降雨の強度を表す降雨強度値のデータ(降雨強度値データ)を生成する。ある時刻における降雨強度値は、当該時刻に降っている雨が1時間降り続いたときに得られる雨量の値である。降雨強度値の単位は、(mm/h)である。図1には、雨水管10が敷設された地上に、降雨レーダ30の観測範囲内にある対象領域を区分した複数の区分領域Mを図示してある。各区分領域は、例えば、62500(250×250)(m)の面積を有する正方形領域である。なお、区分領域の形状は矩形に限られない。例えば、区分領域は、三角形、六角形、円形等による領域であってもよい。信号処理装置40は、各区分領域における降雨強度値データを配信装置50に供給する。 The signal processing device 40 receives a signal supplied by the rain radar 30 and processes the signal to process rainfall intensity value data representing the intensity of rainfall in each of a plurality of segmented regions (mesh) that segment the target region. (Rainfall intensity value data) is generated. The rainfall intensity value at a certain time is a value of the amount of rainfall obtained when the rain falling at the time continues for one hour. The unit of the rainfall intensity value is (mm / h). FIG. 1 illustrates a plurality of divided areas M obtained by dividing a target area within the observation range of the rainfall radar 30 on the ground where the rainwater pipe 10 is laid. Each divided region is a square region having an area of 62500 (250 × 250) (m 2 ), for example. Note that the shape of the segmented region is not limited to a rectangle. For example, the divided area may be an area formed by a triangle, a hexagon, a circle, or the like. The signal processing device 40 supplies the rainfall intensity value data in each segment area to the distribution device 50.

配信装置50は、信号処理装置40が供給する複数の区分領域それぞれに対応する降雨強度値データを取り込み、各降雨強度値データを外部に配信するためのデータフォーマットに変換する。そして、配信装置50は、フォーマット変換後の各降雨強度値データを流量予測装置60に供給する。   The distribution device 50 takes in the rainfall intensity value data corresponding to each of the plurality of divided regions supplied by the signal processing device 40, and converts each rainfall intensity value data into a data format for distribution to the outside. Then, the distribution device 50 supplies each rainfall intensity value data after format conversion to the flow rate prediction device 60.

流量測定装置70は、雨水管10の排出口もしくは排出口付近、または雨水排水施設80の入水口付近に設けられ、雨水排水施設80に流入する雨水の流量を測定する装置である。流量測定装置70は、流量を表す流量値のデータ(流量値データ)を流量予測装置60に供給する。流量値は、単位時間(例えば、1秒)当たりに流れる水量である。流量値の単位は、(m/s)である。 The flow rate measuring device 70 is a device that is provided near the discharge port of the rainwater pipe 10 or near the discharge port, or near the water inlet of the rainwater drainage facility 80, and measures the flow rate of rainwater flowing into the rainwater drainage facility 80. The flow rate measuring device 70 supplies flow rate value data (flow rate value data) representing the flow rate to the flow rate predicting device 60. The flow rate value is the amount of water flowing per unit time (for example, 1 second). The unit of the flow value is (m 3 / s).

流量予測装置60は、一または複数の降雨イベントにおいて、配信装置50が供給する降雨強度値データと、流量測定装置70が供給する流量値データとを、一定時間ごと(例えば、1分または5分ごと)に取り込んで記憶する。言い換えると、流量予測装置60は、対象領域を区分した複数の区分領域それぞれに対応する降雨強度値データと、当該対象領域から流出する雨水の流量値データとを、一または複数の降雨イベントにおける一定時間ごとに取り込み、複数の区分領域分の降雨強度値データ(降雨強度値データセット)と流量値データとを関係付けて時系列に記憶する。これら記憶される時系列データを降雨イベントデータともいう。降雨イベントとは、対象地域の地表に雨が降り始める時点から降り止む時点までの一連の降雨の事象である。ただし、一旦雨が降り止んでから再び降り始めるまでの無降雨時間が所定時間(例えば、4時間)未満である場合は、その無降雨時間前後の降雨を一連の降雨とみなして降雨イベントを区別しない。   The flow rate predicting device 60 uses the rainfall intensity value data supplied from the distribution device 50 and the flow rate value data supplied from the flow measuring device 70 at regular intervals (for example, 1 minute or 5 minutes) in one or more rain events. Every time). In other words, the flow rate predicting device 60 determines the rainfall intensity value data corresponding to each of a plurality of divided areas obtained by dividing the target area and the flow rate value data of rainwater flowing out from the target area at a certain level in one or more rainfall events. The data is taken in every time, and rainfall intensity value data (rainfall intensity value data set) and flow value data for a plurality of divided areas are related and stored in time series. These stored time-series data are also called rain event data. The rain event is a series of rain events from the time when it starts to rain on the surface of the target area to the time when it stops. However, if the rain-free time from when the rain has stopped once until it starts to fall again is less than a predetermined time (for example, 4 hours), the rain before and after the rain-free time is regarded as a series of rain to distinguish the rain event. do not do.

流量予測装置60は、記憶された、一または複数の降雨イベントにおける、降雨強度値データセットと流量値データとを有する時系列データを用いて、流量測定装置70からの未来時刻における雨水の流量値を推定し、推定結果として推定流量値を出力する。流量予測装置60が出力する推定流量値は、例えば、雨水排水施設80における排水ポンプ20、その他の機器の運用状態の監視に利用される。   The flow rate prediction device 60 uses the stored time series data having the rainfall intensity value data set and the flow rate value data in one or a plurality of rainfall events, and the flow rate value of rainwater from the flow rate measurement device 70 at a future time. And an estimated flow rate value is output as an estimation result. The estimated flow rate value output from the flow rate predicting device 60 is used, for example, for monitoring the operation state of the drainage pump 20 and other devices in the rainwater drainage facility 80.

次に、第1実施形態である流量予測装置60の構成について説明する。
図2は、流量予測装置60の機能構成を示すブロック図である。同図に示すように、流量予測装置60は、降雨強度値取得部61と、流量値取得部62と、記憶部63と、流量値推定部64と、係数更新部65とを備える。
Next, the structure of the flow volume prediction apparatus 60 which is 1st Embodiment is demonstrated.
FIG. 2 is a block diagram showing a functional configuration of the flow rate predicting device 60. As shown in the figure, the flow rate prediction device 60 includes a rainfall intensity value acquisition unit 61, a flow rate value acquisition unit 62, a storage unit 63, a flow rate value estimation unit 64, and a coefficient update unit 65.

降雨強度値取得部61は、配信装置50が供給する降雨強度値データを取り込む。具体的に、降雨強度値取得部61は、対象領域を区分した複数の区分領域それぞれ(各メッシュ)において、一または複数の降雨イベントにおける降雨強度値データを一定時間ごとに取得する。つまり、降雨強度値取得部61は、降雨強度値データセットを時系列に取り込む。降雨強度値取得部61は、取り込んだ降雨強度値データセットを記憶部63に記憶させる。   The rainfall intensity value acquisition unit 61 captures rainfall intensity value data supplied by the distribution device 50. Specifically, the rainfall intensity value acquisition unit 61 acquires rainfall intensity value data for one or a plurality of rainfall events at regular intervals in each of a plurality of divided areas (each mesh) obtained by dividing the target area. That is, the rainfall intensity value acquisition unit 61 takes in the rainfall intensity value data set in time series. The rainfall intensity value acquisition unit 61 causes the storage unit 63 to store the captured rainfall intensity value data set.

流量値取得部62は、流量測定装置70が供給する流量値データを取り込む。具体的に、流量値取得部62は、流量値データを、降雨強度値取得部61による取得のタイミングに同期させて取得する。ただし、タイミングに同期させることは、同一の時刻とすることも含む。流量値取得部62は、取り込んだ流量値データを記憶部63に記憶させる。   The flow value acquisition unit 62 takes in the flow value data supplied by the flow measuring device 70. Specifically, the flow value acquisition unit 62 acquires the flow value data in synchronization with the acquisition timing by the rainfall intensity value acquisition unit 61. However, synchronizing with timing includes setting the same time. The flow value acquisition unit 62 stores the acquired flow value data in the storage unit 63.

記憶部63は、降雨強度値取得部61が供給する降雨強度値データセットと、流量値取得部62が供給する流量値データとを関係付けて、順次記憶する。つまり、記憶部63は、一または複数の降雨イベントにおいて一定時間ごとに取り込まれた、複数の区分領域それぞれに対応する降雨強度値データと、対象領域から流出する雨水の流量値データとを関係付けて時系列に記憶する。記憶部63は、例えば、磁気ハードディスク装置、または半導体記憶装置により実現される。   The storage unit 63 sequentially associates and stores the rainfall intensity value data set supplied by the rainfall intensity value acquisition unit 61 and the flow rate value data supplied by the flow rate value acquisition unit 62. That is, the storage unit 63 correlates the rainfall intensity value data corresponding to each of the plurality of segmented regions, which is captured at regular intervals in one or a plurality of rainfall events, and the flow rate value data of rainwater flowing out from the target region. And memorize them in time series. The storage unit 63 is realized by, for example, a magnetic hard disk device or a semiconductor storage device.

流量値推定部64は、記憶部63に記憶された、一または複数の降雨イベントにおける降雨強度値データセットおよび流量値データの時系列データ(つまり、過去の降雨イベントデータ)に基づいて、対象領域から流出する雨水の未来時刻における流量値を推定する。具体的に、流量値推定部64は、一または複数の降雨イベントにおける、降雨強度値データセットおよび流量値データの時系列データを記憶部63から読み込む。そして、流量値推定部64は、読み込んだ時系列データに基づいて、重回帰分析処理を実行することによって複数の区分領域における降雨強度値に対する流量値の予測式(流量推定式)の係数(回帰係数)を算出し、これら係数を内蔵する記憶部に記憶させる。なお、ここでは、係数は切片および誤差も含む。流量値推定部64は、記憶した係数を適用した予測式の計算処理を実行することによって、未来時刻における流量値(推定流量値)を推定する。流量値推定部64は、算出した推定流量値のデータ(推定流量値データ)を出力する。   The flow value estimation unit 64 is based on the rainfall intensity value data set and the time series data of the flow value data (that is, past rainfall event data) in one or a plurality of rainfall events stored in the storage unit 63. Estimate the flow rate of rainwater flowing out from the future time. Specifically, the flow value estimation unit 64 reads the rainfall intensity value data set and the time series data of the flow value data from one or more rain events from the storage unit 63. Then, the flow value estimation unit 64 executes a multiple regression analysis process based on the read time-series data, thereby calculating a coefficient (regression) of a prediction formula (flow estimation formula) of the flow value with respect to the rainfall intensity value in a plurality of divided regions. (Coefficients) are calculated and stored in a storage unit incorporating these coefficients. Here, the coefficient includes an intercept and an error. The flow rate estimation unit 64 estimates a flow rate value (estimated flow rate value) at a future time by executing a calculation process of a prediction formula to which the stored coefficient is applied. The flow rate estimation unit 64 outputs the calculated estimated flow rate data (estimated flow rate data).

係数更新部65は、流量値推定部64が算出した係数で、流量値推定部64の記憶部に記憶している係数を所定のタイミングで更新する。例えば、係数更新部65は、流量値推定部64が係数を算出するたびに、算出された係数で記憶部に記憶された係数を書き換える(更新する)。または、係数更新部65は、更新後一定時間(例えば、数時間または数日等)が経過してから更新する。または、係数更新部65は、内蔵した時計の計時結果に応じて、あらかじめ決められた日時を経過してから更新する。   The coefficient update unit 65 updates the coefficient stored in the storage unit of the flow value estimation unit 64 at a predetermined timing with the coefficient calculated by the flow value estimation unit 64. For example, every time the flow rate estimation unit 64 calculates a coefficient, the coefficient update unit 65 rewrites (updates) the coefficient stored in the storage unit with the calculated coefficient. Alternatively, the coefficient updating unit 65 updates after a certain time (for example, several hours or days) after the update. Alternatively, the coefficient updating unit 65 updates after a predetermined date and time has elapsed according to the time measurement result of the built-in clock.

次に、流量予測装置60が取り込む降雨イベントデータについて説明する。
図3は、対象領域を区分した複数の区分領域それぞれに対応する降雨強度値と、当該対象領域から流出する雨水の流量値とを模式的に表した図である。同図は、ある時刻における、対象領域を区分したn個(nは、2以上の整数)の区分領域に対応した降雨強度値r,r,・・・,rと、当該対象領域から流出する雨水の流量値Qとを表す。流量予測装置60は、一または複数の降雨イベントにおいて、同図に示す降雨イベントデータを一定時間ごとに取り込んで記憶する。
Next, the rain event data captured by the flow rate prediction device 60 will be described.
FIG. 3 is a diagram schematically showing the rainfall intensity value corresponding to each of a plurality of divided areas obtained by dividing the target area and the flow value of rainwater flowing out from the target area. The figure, at a certain time, (the n, 2 or more integer) n number obtained by dividing the target area rainfall intensity value corresponding to the segmented region of the r 1, r 2, · · ·, and r n, the target area Represents the flow rate value Q of rainwater flowing out from the water. The flow rate predicting device 60 captures and stores the rain event data shown in the figure at regular intervals in one or a plurality of rain events.

図4は、流量予測装置60が取り込む降雨イベントデータの具体例を示すデータ構成図である。同図によれば、流量予測装置60は、降雨強度値データセット(r,r,・・・,r)および流量値Qを、1分ごとに取り込んでいることを示している。 FIG. 4 is a data configuration diagram illustrating a specific example of the rain event data captured by the flow rate prediction device 60. According to the figure, the flow rate predicting device 60 indicates that the rain intensity value data set (r 1, r 2, ··· , r n) a and flow rate value Q, is taken every minute.

次に、流量予測装置60における流量値推定部64が実行する重回帰分析処理について説明する。重回帰分析処理においては、区分領域(メッシュ)の数量をn、対象とする時間範囲を現在時刻からm(単位は分であるとする。)前までの時間とすると、n×(m+1)よりも多いデータ数が必要である。例えば、20個の区分領域に区分される対象領域で、現在時刻から20分前までの降雨イベントデータを用いて重回帰分析処理を実行するためには、20×(20+1)=420個以上のデータが必要である。   Next, the multiple regression analysis process executed by the flow rate estimation unit 64 in the flow rate prediction device 60 will be described. In the multiple regression analysis process, if the number of segmented areas (mesh) is n, and the target time range is the time from the current time to m (the unit is minutes), n × (m + 1) A large number of data is required. For example, in order to execute the multiple regression analysis process using the rain event data from the current time to 20 minutes before the target area divided into 20 divided areas, 20 × (20 + 1) = 420 or more I need data.

流量値推定部64は、記憶部64に記憶された、時刻(t−m)から時刻tまでの降雨イベントデータを読み込んで重回帰分析処理を実行することにより、下記の式(1)に示す回帰式における係数ai,j、切片a、および誤差εを算出する。式(1)において、r(t−j)は、時刻(t−j)におけるi番目のメッシュに対応する降雨強度値である。ただし、(i=1,2,・・・,n)、(j=0,1,2,・・・,m)である。また、Q(ハット)(t+1)は、時刻(t+1)における推定流量値である。なお、“(ハット)”は、その直前の文字“Q”の上部に付される“^”のことである。 The flow value estimation unit 64 reads the rain event data from the time (t−m) to the time t stored in the storage unit 64 and executes the multiple regression analysis process, thereby expressing the following equation (1). A coefficient a i, j , intercept a 0 , and error ε in the regression equation are calculated. In Expression (1), r i (t−j) is a rainfall intensity value corresponding to the i-th mesh at time (t−j). However, (i = 1, 2,..., N) and (j = 0, 1, 2,..., M). Q (hat) (t + 1) is an estimated flow rate value at time (t + 1). Note that “(hat)” is “^” added to the upper part of the character “Q” immediately before it.

Figure 2014234675
Figure 2014234675

流量値推定部64は、算出した係数ai,j、切片a、および誤差εを適用した(1)式による回帰式に、時刻(t−m)から時刻tまでの降雨強度値データセットを代入して計算することにより、時刻(t+1)における対象領域からの推定流量値Q(ハット)(t+1)を得る。 The flow rate estimation unit 64 uses the calculated coefficient a i, j , intercept a 0 , and error ε as a regression equation based on the equation (1), and sets the rainfall intensity value data set from time (tm) to time t. By substituting and calculating, an estimated flow rate value Q (hat) (t + 1) from the target region at time (t + 1) is obtained.

次に、第1実施形態である流量予測装置60の動作について説明する。
図5は、流量予測装置60の処理を示すフローチャートである。
ステップS1において、流量予測装置60は、降雨イベントが開始されるまで待機している(ステップS1:NO)。具体的に、流量予測装置60は、無降雨時間が所定時間(例えば、4時間)以上経過したか否かを判定する。無降雨時間が所定時間以上経過して雨が降り始めた場合、流量予測装置60は、降雨イベントが開始したと判定する。降雨イベントが開始したと流量予測装置60が判定した場合(ステップS1:YES)、ステップS2の処理に移す。
Next, operation | movement of the flow volume prediction apparatus 60 which is 1st Embodiment is demonstrated.
FIG. 5 is a flowchart showing the process of the flow rate prediction device 60.
In step S1, the flow rate predicting device 60 stands by until a rainfall event is started (step S1: NO). Specifically, the flow rate predicting device 60 determines whether or not a non-rainy time has passed a predetermined time (for example, 4 hours) or more. When the rainless time has elapsed for a predetermined time or more and the rain starts, the flow rate prediction device 60 determines that the rain event has started. When the flow rate prediction device 60 determines that a rain event has started (step S1: YES), the process proceeds to step S2.

ステップS2において、降雨強度値取得部61は、配信装置50が供給する降雨強度値データセットを取り込み、この降雨強度値データセットを記憶部64に記憶させる。また、流量値取得部62は、流量測定装置70が供給する流量値データを、降雨強度値取得部61による取得のタイミングに同期させて取り込み、その流量値データを降雨強度値データセットに関係付けて記憶部64に記憶させる。   In step S <b> 2, the rainfall intensity value acquisition unit 61 takes in the rainfall intensity value data set supplied by the distribution device 50 and stores the rainfall intensity value data set in the storage unit 64. The flow value acquisition unit 62 captures the flow value data supplied from the flow measurement device 70 in synchronization with the acquisition timing by the rainfall intensity value acquisition unit 61, and associates the flow value data with the rainfall intensity value data set. To be stored in the storage unit 64.

次に、ステップS3において、流量予測装置60は、降雨イベントデータの取り込みを終了するか否かを判定する。例えば、流量予測装置60は、ある一定の降雨量以上の降雨となる一または複数の降雨イベントにおいて降雨イベントデータを取得したか否かを判定する。流量予測装置60は、降雨イベントデータの取り込みを終了しないと判定した場合(ステップS3:NO)、ステップS1の処理に戻す。一方、流量予測装置60は、降雨イベントデータの取り込みを終了すると判定した場合(ステップS3:YES)、ステップS4の処理に移す。   Next, in step S <b> 3, the flow rate prediction device 60 determines whether or not to end the rain event data capture. For example, the flow rate predicting device 60 determines whether or not the rain event data has been acquired in one or a plurality of rain events in which the amount of rainfall exceeds a certain amount of rainfall. If the flow rate predicting device 60 determines not to finish capturing the rain event data (step S3: NO), it returns to the process of step S1. On the other hand, when it determines with the flow volume prediction apparatus 60 complete | finishing taking in of rain event data (step S3: YES), it moves to the process of step S4.

ステップS4において、流量値推定部64は、一または複数の降雨イベントにおける、降雨強度値データセットおよび流量値データの時系列データを記憶部63から読み込む。次に、流量値推定部64は、読み込んだ時系列データに基づいて、重回帰分析処理を実行することによって複数の区分領域における降雨強度値に対する流量値の予測式の係数を算出し、これら係数を内蔵する記憶部に記憶させる。   In step S <b> 4, the flow value estimation unit 64 reads the rainfall intensity value data set and the time series data of the flow value data from one or more rain events from the storage unit 63. Next, the flow value estimation unit 64 calculates the coefficient of the prediction formula of the flow value with respect to the rainfall intensity value in the plurality of divided regions by executing the multiple regression analysis process based on the read time series data, and these coefficients Is stored in a built-in storage unit.

次に、ステップS5において、流量値推定部64は、記憶した係数を適用した予測式の計算処理を実行することによって、未来時刻における推定流量値を算出する。次に、流量値推定部64は、算出した推定流量値のデータである推定流量値データを出力する。   Next, in step S <b> 5, the flow rate estimation unit 64 calculates an estimated flow rate value at a future time by executing a prediction formula calculation process using the stored coefficient. Next, the flow value estimation unit 64 outputs estimated flow value data that is data of the calculated estimated flow value.

次に、ステップS6において、係数更新部65は、流量値推定部64が算出した係数で、流量値推定部64に記憶された係数を更新するか否かを判定する。例えば、係数更新部65は、流量値推定部64が係数を算出したら必ず更新すると判定する。または、係数更新部65は、前回の更新から一定時間が経過したか否かを判定する。ただし、未だ更新されたことがない場合は、必ず更新すると判定する。または、係数更新部65は、内蔵する時計により、あらかじめ決定された日時(例えば、6月1日午前0時0分)になったか否かを判定する。そして、係数更新部65は、更新すると判定した場合(ステップS6:YES)、ステップS7の処理に移し、更新しないと判定した場合(ステップS6:NO)、ステップS1の処理に戻す。   Next, in step S <b> 6, the coefficient update unit 65 determines whether or not to update the coefficient stored in the flow value estimation unit 64 with the coefficient calculated by the flow value estimation unit 64. For example, the coefficient updating unit 65 determines to update the flow rate estimation unit 64 whenever the flow rate estimation unit 64 calculates the coefficient. Alternatively, the coefficient updating unit 65 determines whether a certain time has elapsed since the previous update. However, if it has not been updated yet, it is determined that it will always be updated. Alternatively, the coefficient updating unit 65 determines whether or not a predetermined date and time (for example, 0:00 am on June 1) has been reached using a built-in clock. If the coefficient update unit 65 determines to update (step S6: YES), the coefficient update unit 65 proceeds to the process of step S7. If determined not to update (step S6: NO), the coefficient update unit 65 returns to the process of step S1.

ステップS7において、係数更新部65は、ステップS5の処理において流量値推定部64が算出した係数で、流量値推定部64に記憶している係数を更新する。次に、流量予測装置60は、ステップS1の処理に戻す。   In step S7, the coefficient updating unit 65 updates the coefficient stored in the flow value estimation unit 64 with the coefficient calculated by the flow value estimation unit 64 in the process of step S5. Next, the flow rate prediction device 60 returns to the process of step S1.

次に、流量予測装置60による、対象領域における等到達時間曲線の取得方法について説明する。
図6は、流量値推定部64に記憶された係数のうち、k番目の区分領域(メッシュk)における各係数の値を示したグラフである。同図において、横軸はメッシュkにおける係数、縦軸は係数の値である。なお、0≦l(エル)≦mである。同図に示すように、メッシュkにおける係数ak,0から係数ak,mまでのうち、係数ak,lが、係数の値として最大となる係数である。よって、メッシュkに対応する降雨強度値の取得時点から、対象領域からの雨水の流量時点までの間には、時間l(エル)分の遅れ時間があるといえる。
Next, the acquisition method of the equal arrival time curve in the target region by the flow rate prediction device 60 will be described.
FIG. 6 is a graph showing the value of each coefficient in the kth segmented area (mesh k) among the coefficients stored in the flow rate estimation unit 64. In the figure, the horizontal axis represents the coefficient in the mesh k, and the vertical axis represents the coefficient value. Note that 0 ≦ l (el) ≦ m. As shown in the figure, among the coefficients a k, 0 to the coefficients a k, m in the mesh k, the coefficient a k, l is the maximum coefficient value. Therefore, it can be said that there is a delay time of time l (el) between the time point when the rainfall intensity value corresponding to the mesh k is acquired and the time point when the rainwater flows from the target region.

流量予測装置60は、対象領域における全ての区分領域それぞれについて、係数の値が最大となる係数を検出し、検出した係数から遅れ時間を取得する。そして、流量予測装置60は、対象領域の全ての区分領域における遅れ時間を区分領域に対応させた二次元データを生成し、この二次元データにおいて同一または同等の遅れ時間を相互に対応付けることにより、等到達時間を得る。   The flow rate predicting device 60 detects the coefficient having the maximum coefficient value for each of the divided areas in the target area, and acquires the delay time from the detected coefficient. Then, the flow rate predicting device 60 generates two-dimensional data in which the delay times in all the divided areas of the target area are associated with the divided areas, and by associating the same or equivalent delay times with each other in the two-dimensional data, Get equal arrival time.

図7は、流量予測装置60が生成した二次元データ上で、同一または同等の遅れ時間を対応付けることによって得た等到達時間曲線を表した図である。同図では、視覚的に見やすくするよう、等到達時間を線で表したが、例えば、メッシュを色分けしたり、三次元方向の長さを変えたりして、等到達時間を表してもよい。
等到達時間を求めることにより、対象領域における流下特性を視覚的にとらえることができる。
FIG. 7 is a diagram showing an equal arrival time curve obtained by associating the same or equivalent delay times on the two-dimensional data generated by the flow rate prediction device 60. In the figure, the equal arrival time is represented by a line so as to be visually easy to see, but the equal arrival time may be represented by, for example, color-coding the mesh or changing the length in the three-dimensional direction.
By obtaining the equal arrival time, it is possible to visually grasp the flow characteristics in the target region.

第1実施形態である流量予測装置は、対象領域を区分した複数の区分領域それぞれにおいて、降雨イベントにおける降雨強度値をあらかじめ一定時間ごとに取得する降雨強度値取得部を備える。また、この流量予測装置は、対象領域から流出する雨水の流量値を、降雨強度値取得部による取得のタイミングに同期させて取得する流量値取得部を備える。また、この流量予測装置は、降雨強度値取得部が取得した降雨強度値と流量値取得部が取得した流量値とに基づいて、対象領域から流出する雨水の未来時刻における流量値を推定する流量値推定部を備える。   The flow rate prediction apparatus according to the first embodiment includes a rainfall intensity value acquisition unit that acquires, in advance, a rainfall intensity value in a rain event every predetermined time in each of a plurality of divided areas obtained by dividing the target area. In addition, the flow rate prediction apparatus includes a flow value acquisition unit that acquires the flow value of rainwater flowing out from the target region in synchronization with the acquisition timing of the rainfall intensity value acquisition unit. In addition, the flow rate predicting apparatus estimates a flow rate value at a future time of rainwater flowing out from the target area based on the rainfall intensity value acquired by the rainfall intensity value acquiring unit and the flow rate value acquired by the flow rate acquiring unit. A value estimation unit is provided.

第1実施形態である流量予測装置は、降雨イベントについて一定時間ごとに計測した、対象領域における降雨強度値および流量値に基づいて、未来時刻における当該対象領域からの雨水の流量を推定する。よって、第1実施形態である流量予測装置によれば、エンジニアリング作業量を少なく抑えた簡便な方法によって雨水流量の予測を行うことができる。   The flow rate predicting apparatus according to the first embodiment estimates the flow rate of rainwater from the target region at a future time based on the rainfall intensity value and the flow rate value in the target region measured at regular intervals for the rain event. Therefore, according to the flow volume prediction apparatus which is 1st Embodiment, the rainwater flow volume can be estimated with the simple method which restrained the engineering work amount small.

また、第1実施形態である流量予測装置における流量値推定部は、それぞれ取得した降雨強度値および流量値に基づいて、重回帰分析処理を実行することによって複数の区分領域における降雨強度値に対する流量値の予測式の係数を算出し、これら係数を適用した予測式を計算することによって未来時刻における流量値を推定する。   In addition, the flow rate estimation unit in the flow prediction device according to the first embodiment performs a multiple regression analysis process on the basis of the acquired rainfall intensity value and the flow rate value, and thereby the flow rate for the rainfall intensity value in the plurality of segmented regions. A coefficient of a value prediction formula is calculated, and a flow rate value at a future time is estimated by calculating a prediction formula to which these coefficients are applied.

これにより、第1実施形態である流量予測装置によれば、エンジニアリング作業量を少なく抑えた簡便な方法によって雨水流量の予測を精度よく行うことができる。   Thereby, according to the flow volume prediction apparatus which is 1st Embodiment, the rainwater flow volume can be accurately estimated by the simple method which restrained the engineering work amount small.

また、第1実施形態である流量予測装置は、流量値推定部が算出する係数を所定のタイミングで更新する係数更新部を備える。更新のタイミングは、常に最新の算出結果である係数を保つように制御されてもよいし、あらかじめ決められたときに更新されるよう制御されてもよい。   In addition, the flow rate prediction apparatus according to the first embodiment includes a coefficient update unit that updates a coefficient calculated by the flow value estimation unit at a predetermined timing. The update timing may be controlled so as to always keep the coefficient which is the latest calculation result, or may be controlled to be updated when determined in advance.

これにより、第1実施形態である流量予測装置は、適切な係数を用いることによって、雨水流量の予測を精度よく行うことができる。   Thereby, the flow volume prediction apparatus which is 1st Embodiment can predict a rainwater flow volume accurately by using a suitable coefficient.

[第2の実施形態]
第1実施形態における降雨強度値取得部61は、配信装置50から取り込んだ降雨強度値データの降雨強度値をそのままの値として記憶部63に記憶させた。すなわち、第1実施形態では、降雨レーダ30が生成した信号に基づき信号処理装置40が生成した区分領域ごとの降雨強度値データをそのまま利用していたため、地形変化に起因したノイズや信号の減衰による降雨強度値への影響が考慮されていなかった。第2実施形態は、降雨強度値取得部が、取得した降雨強度値データを、降雨強度に応じて区分した複数のレンジに対応させて正規化する例である。
[Second Embodiment]
The rainfall intensity value acquisition unit 61 in the first embodiment stores the rainfall intensity value of the rainfall intensity value data captured from the distribution device 50 in the storage unit 63 as it is. That is, in the first embodiment, since the rainfall intensity value data for each divided region generated by the signal processing device 40 based on the signal generated by the rain radar 30 is used as it is, it is caused by noise or signal attenuation due to terrain change. The effect on rainfall intensity was not considered. The second embodiment is an example in which the rainfall intensity value acquisition unit normalizes the acquired rainfall intensity value data in association with a plurality of ranges divided according to the rainfall intensity.

第2実施形態である流量予測装置の機能構成図は、第1実施形態である流量予測装置60の機能構成を示した図2を用いることとする。第2実施形態においては、第1実施形態と相違する部分について説明し、共通する部分についての説明を省略する。   As the functional configuration diagram of the flow rate prediction apparatus according to the second embodiment, FIG. 2 showing the functional configuration of the flow rate prediction apparatus 60 according to the first embodiment is used. In the second embodiment, parts that are different from the first embodiment will be described, and descriptions of common parts will be omitted.

第2実施形態における降雨強度値取得部61は、第1実施形態と同様に、配信装置50が供給する降雨強度値データを取り込む。そして、降雨強度値取得部61は、取得した降雨強度値データを、降雨強度に応じて区分した複数のレンジに対応させて正規化する。   The rainfall intensity value acquisition unit 61 in the second embodiment takes in the rainfall intensity value data supplied by the distribution device 50 as in the first embodiment. Then, the rainfall intensity value acquisition unit 61 normalizes the acquired rainfall intensity value data in correspondence with a plurality of ranges divided according to the rainfall intensity.

具体的に、降雨強度値取得部61は、複数のレンジを示すレンジデータを、内蔵する記憶部にあらかじめ記憶させておく。降雨強度値が200(mm/h)以下である場合の例として、降雨強度値取得部61の記憶部は、例えば、0(mm/h)、0(mm/h)より大きく1(mm/h)以下、1(mm/h)より大きく5(mm/h)以下、5(mm/h)より大きく10(mm/h)以下、10(mm/h)より大きく20(mm/h)以下、20(mm/h)より大きく50(mm/h)以下、50(mm/h)より大きく100(mm/h)以下、および100(mm/h)より大きく200(mm/h)以下を示すレンジデータをあらかじめ記憶する。   Specifically, the rainfall intensity value acquisition unit 61 stores range data indicating a plurality of ranges in advance in a built-in storage unit. As an example when the rainfall intensity value is 200 (mm / h) or less, the storage unit of the rainfall intensity value acquisition unit 61 is, for example, 0 (mm / h), greater than 0 (mm / h), and 1 (mm / h). h) or less, greater than 1 (mm / h), 5 (mm / h) or less, greater than 5 (mm / h), less than 10 (mm / h), greater than 10 (mm / h), 20 (mm / h) Hereinafter, greater than 20 (mm / h) and not greater than 50 (mm / h), greater than 50 (mm / h) and not greater than 100 (mm / h), and greater than 100 (mm / h) and not greater than 200 (mm / h) Is stored in advance.

降雨強度値取得部61は、取り込んだ降雨強度値データを記憶されたレンジデータにおけるレンジに応じて、該当するレンジの中央値に換算することによって正規化する。つまり、降雨強度値取得部61は、取り込んだ降雨強度値データを、下記の表1のように換算する。   The rainfall intensity value acquisition unit 61 normalizes the captured rainfall intensity value data by converting it into the median value of the corresponding range according to the range in the stored range data. That is, the rainfall intensity value acquisition unit 61 converts the captured rainfall intensity value data as shown in Table 1 below.

Figure 2014234675
Figure 2014234675

第2実施形態における降雨強度値取得部61は、取得した降雨強度値データを、降雨強度に応じて区分した複数のレンジに対応させて正規化することにより、簡易な構成により正規化前の降雨強度値データに含まれる誤差を吸収して、雨水流量の予測を精度よく行うことができる。   The rainfall intensity value acquisition unit 61 in the second embodiment normalizes the acquired rainfall intensity value data in correspondence with a plurality of ranges divided according to the rainfall intensity, so that the rainfall before normalization can be achieved with a simple configuration. The error contained in the intensity value data can be absorbed, and the rainwater flow rate can be predicted accurately.

[第3の実施形態]
第3実施形態は、降雨強度値取得部が、取得した降雨強度値データを、降雨強度に応じて区分した複数のレンジに対応させて正規化する例であり、第2の実施形態の変形例である。
[Third Embodiment]
The third embodiment is an example in which the rainfall intensity value acquisition unit normalizes the acquired rainfall intensity value data corresponding to a plurality of ranges divided according to the rainfall intensity, and is a modification of the second embodiment. It is.

第3実施形態である流量予測装置の機能構成図は、第1実施形態における流量予測装置60の機能構成を示した図2を用いることとする。第3実施形態においては、第2実施形態と相違する構成について説明し、第1実施形態および第2実施形態と共通する構成についての説明を省略する。   As the functional configuration diagram of the flow rate prediction apparatus according to the third embodiment, FIG. 2 showing the functional configuration of the flow rate prediction apparatus 60 according to the first embodiment is used. In the third embodiment, a configuration that is different from the second embodiment will be described, and a description of a configuration that is common to the first embodiment and the second embodiment will be omitted.

第3実施形態における降雨強度値取得部61は、第1実施形態と同様に、配信装置50が供給する降雨強度値データセットを時系列に取り込む。そして、降雨強度値取得部61は、取得した降雨強度値データセットを、降雨強度に応じて区分した複数のレンジに対応させて正規化する。   As in the first embodiment, the rainfall intensity value acquisition unit 61 in the third embodiment captures the rainfall intensity value data set supplied by the distribution device 50 in time series. Then, the rainfall intensity value acquisition unit 61 normalizes the acquired rainfall intensity value data set in association with a plurality of ranges divided according to the rainfall intensity.

具体的に、降雨強度値取得部61は、第2実施形態と同様に、複数のレンジを示すレンジデータを、内蔵する記憶部にあらかじめ記憶させておく。ここでのレンジデータの内容は、第2実施形態と同一の例とする。   Specifically, the rainfall intensity value acquisition unit 61 stores range data indicating a plurality of ranges in advance in a built-in storage unit, as in the second embodiment. The contents of the range data here are the same examples as in the second embodiment.

降雨強度値取得部61は、取り込んだ降雨強度値データセットにおいて、降雨強度値が最大となる降雨強度値データについては換算せず、その他の降雨強度値データを、記憶されたレンジデータにおけるレンジに応じて、該当するレンジの中央値に換算することによって正規化する。例えば、降雨強度値データセットに最大値として40(mm/h)の降雨強度値データが含まれている場合、降雨強度値取得部61は、取り込んだ降雨強度値データセットを、下記の表2のように換算する。   The rainfall intensity value acquisition unit 61 does not convert the rainfall intensity value data having the maximum rainfall intensity value in the captured rainfall intensity value data set, and converts other rainfall intensity value data into the range in the stored range data. Accordingly, normalization is performed by converting to the median of the corresponding range. For example, when the rainfall intensity value data set includes 40 (mm / h) of rainfall intensity value data as the maximum value, the rainfall intensity value acquisition unit 61 selects the captured rainfall intensity value data set as shown in Table 2 below. Convert as follows.

Figure 2014234675
Figure 2014234675

第3実施形態における降雨強度値取得部61は、取得した降雨強度値データセットを、降雨強度値が最大となる降雨強度値データを除き、降雨強度に応じて区分した複数のレンジに対応させて正規化することにより、簡易な構成により正規化前の降雨強度値データに含まれる誤差を吸収して、雨水流量の予測を精度よく行うことができる。   The rainfall intensity value acquisition unit 61 according to the third embodiment associates the acquired rainfall intensity value data set with a plurality of ranges divided according to the rainfall intensity, except for the rainfall intensity value data having the maximum rainfall intensity value. By normalizing, the error contained in the rainfall intensity value data before normalization can be absorbed with a simple configuration, and the rainwater flow rate can be predicted accurately.

[第4の実施形態]
第4実施形態は、降雨強度値取得部が、取得した降雨強度値データを、降雨強度に応じて区分した複数のレンジに対応させて正規化する例であり、第3の実施形態の変形例である。
[Fourth Embodiment]
The fourth embodiment is an example in which the rainfall intensity value acquisition unit normalizes the acquired rainfall intensity value data in association with a plurality of ranges divided according to the rainfall intensity, and is a modification of the third embodiment. It is.

第4実施形態である流量予測装置の機能構成図は、第1実施形態における流量予測装置60の機能構成を示した図2を用いることとする。第4実施形態においては、第3実施形態と相違する構成について説明し、第1実施形態から第3実施形態までと共通する構成についての説明を省略する。   As the functional configuration diagram of the flow rate prediction apparatus according to the fourth embodiment, FIG. 2 showing the functional configuration of the flow rate prediction apparatus 60 according to the first embodiment is used. In the fourth embodiment, the configuration different from the third embodiment will be described, and the description of the configuration common to the first to third embodiments will be omitted.

第4実施形態における降雨強度値取得部61は、第1実施形態と同様に、配信装置50が供給する降雨強度値データセットを時系列に取り込む。そして、降雨強度値取得部61は、取得した降雨強度値データセットを、降雨強度に応じて区分した複数のレンジに対応させて正規化する。   The rainfall intensity value acquisition unit 61 in the fourth embodiment takes in the rainfall intensity value data set supplied by the distribution device 50 in time series, as in the first embodiment. Then, the rainfall intensity value acquisition unit 61 normalizes the acquired rainfall intensity value data set in association with a plurality of ranges divided according to the rainfall intensity.

具体的に、降雨強度値取得部61は、第2実施形態と同様に、複数のレンジを示すレンジデータを、内蔵する記憶部にあらかじめ記憶させておく。ここでのレンジデータの内容は、第2実施形態と同一の例とする。   Specifically, the rainfall intensity value acquisition unit 61 stores range data indicating a plurality of ranges in advance in a built-in storage unit, as in the second embodiment. The contents of the range data here are the same examples as in the second embodiment.

降雨強度値取得部61は、取り込んだ降雨強度値データセットにおいて、降雨強度値が最大となる降雨強度値データ(最大降雨強度値データ)については換算しない。そして、降雨強度値取得部61は、その他の降雨強度値データを、最大降雨強度値データが含まれるレンジ内での最大降雨強度値の割合を基準に、最少レンジまでの間の各レンジにおいて一定の割合で増加または減少する割合に対応した値に換算することによって正規化する。例えば、降雨強度値データセットに最大値として75(mm/h)の降雨強度値データが含まれている場合、降雨強度値取得部61は、取り込んだ降雨強度値データセットを、下記の表3のように換算する。   The rainfall intensity value acquisition unit 61 does not convert the rainfall intensity value data (maximum rainfall intensity value data) having the maximum rainfall intensity value in the captured rainfall intensity value data set. Then, the rainfall intensity value acquisition unit 61 makes other rainfall intensity value data constant in each range up to the minimum range based on the ratio of the maximum rainfall intensity value in the range including the maximum rainfall intensity value data. Normalize by converting to a value corresponding to the rate of increase or decrease at a rate of. For example, when the rainfall intensity value data set includes the rainfall intensity value data of 75 (mm / h) as the maximum value, the rainfall intensity value acquisition unit 61 sets the captured rainfall intensity value data set as the following Table 3. Convert as follows.

Figure 2014234675
Figure 2014234675

第4実施形態における降雨強度値取得部61は、取得した降雨強度値データセットを、降雨強度値が最大となる降雨強度値データを除き、降雨強度に応じて区分した複数のレンジに対応させて正規化することにより、簡易な構成により正規化前の降雨強度値データに含まれる誤差を吸収して、雨水流量の予測を精度よく行うことができる。   The rainfall intensity value acquisition unit 61 according to the fourth embodiment associates the acquired rainfall intensity value data set with a plurality of ranges divided according to the rainfall intensity, except for the rainfall intensity value data having the maximum rainfall intensity value. By normalizing, the error contained in the rainfall intensity value data before normalization can be absorbed with a simple configuration, and the rainwater flow rate can be predicted accurately.

[第5の実施形態]
第5実施形態は、流量値推定部が、降雨強度に応じて区分した複数のレンジそれぞれについて予測式における係数を算出し、レンジごとに予測式を計算する例である。
[Fifth Embodiment]
The fifth embodiment is an example in which the flow value estimation unit calculates a coefficient in a prediction formula for each of a plurality of ranges divided according to rainfall intensity, and calculates the prediction formula for each range.

第5実施形態である流量予測装置の機能構成図は、第1実施形態である流量予測装置60の機能構成を示した図2を用いることとする。第5実施形態においては、第1実施形態と相違する部分について説明し、共通する部分についての説明を省略する。   The functional configuration diagram of the flow rate predicting apparatus according to the fifth embodiment uses FIG. 2 showing the functional configuration of the flow rate predicting apparatus 60 according to the first embodiment. In the fifth embodiment, portions that are different from the first embodiment will be described, and descriptions of common portions will be omitted.

流量値推定部64は、降雨強度に応じて区分した複数のレンジそれぞれに対応させて、一または複数の降雨イベントにおける、降雨強度値データセットおよび流量値データの時系列データを記憶部63から読み込む。降雨強度値が200(mm/h)以下である場合の例として、レンジは、例えば、0(mm/h)、0(mm/h)より大きく1(mm/h)以下、1(mm/h)より大きく5(mm/h)以下、5(mm/h)より大きく10(mm/h)以下、10(mm/h)より大きく20(mm/h)以下、20(mm/h)より大きく50(mm/h)以下、50(mm/h)より大きく100(mm/h)以下、および100(mm/h)より大きく200(mm/h)以下である。   The flow value estimation unit 64 reads the rainfall intensity value data set and the time series data of the flow value data from the storage unit 63 in one or a plurality of rain events in association with each of a plurality of ranges divided according to the rainfall intensity. . As an example when the rainfall intensity value is 200 (mm / h) or less, the range is, for example, 0 (mm / h), greater than 0 (mm / h) and 1 (mm / h) or less, 1 (mm / h) h) greater than 5 (mm / h) and greater than 5 (mm / h) greater than 10 (mm / h) and greater than 10 (mm / h) and not greater than 20 (mm / h) and 20 (mm / h) It is larger than 50 (mm / h), larger than 50 (mm / h) and smaller than 100 (mm / h), and larger than 100 (mm / h) and smaller than or equal to 200 (mm / h).

そして、レンジごとに、流量値推定部64は、読み込んだ時系列データに基づいて、重回帰分析処理を実行することによって複数の区分領域における降雨強度値に対する流量値の予測式の係数を算出し、これら係数を記憶する。なお、係数は、切片および誤差も含む。そして、レンジごとに、流量値推定部64は、記憶した係数を適用した予測式を計算することによって未来時刻における流量値を推定する。流量値推定部64は、算出した推定流量値の推定流量値データを出力する。すなわち、流量値推定部64は、降雨強度に応じて区分した複数のレンジそれぞれについて予測式における係数を算出し、レンジごとに予測式を計算することによって、未来時刻における流量値をレンジごとに推定する。   Then, for each range, the flow value estimation unit 64 calculates the coefficient of the prediction formula of the flow value with respect to the rainfall intensity values in the plurality of divided regions by executing the multiple regression analysis process based on the read time series data. These coefficients are stored. The coefficient includes an intercept and an error. For each range, the flow rate estimation unit 64 estimates a flow rate value at a future time by calculating a prediction formula to which the stored coefficient is applied. The flow value estimation unit 64 outputs estimated flow value data of the calculated estimated flow value. That is, the flow rate estimation unit 64 calculates a coefficient in a prediction formula for each of a plurality of ranges divided according to rainfall intensity, and estimates a flow rate value at a future time for each range by calculating a prediction formula for each range. To do.

これにより、対象領域への降雨の度合いによって対象領域からの雨水の流量の度合いが変わる場合、言い換えると、降雨の度合いに応じて予測式の係数が変わる場合にも、雨水流量の予測を精度よく行うことができる。   As a result, when the degree of rainwater flow from the target area changes depending on the degree of rainfall in the target area, in other words, even when the coefficient of the prediction formula changes according to the degree of rain, the prediction of the rainwater flow rate is accurate. It can be carried out.

[第6の実施形態]
対象領域を区分した複数の区分領域において、雨水排水施設80が設置された場所から離れた区分領域ほど、当該対象領域からの流量に寄与するまでには時間がかかる。そこで、第6実施形態は、流量値推定部が、降雨イベントにおける経過時間が長いほど、対象領域における流出箇所から離れた区分領域における降雨強度値に対する係数をも算出する例である。言い換えると、第6実施形態における流量値推定部は、降雨イベントにおける経過時間が短いときは、対象領域における流出箇所に近い区分領域における降雨強度値に対する係数を算出し、経過時間が長くなるほど、流出箇所側から区分領域を広げるように算出する係数を増やしていく。
[Sixth Embodiment]
In a plurality of divided areas obtained by dividing the target area, it takes time to contribute to the flow rate from the target area as the divided area is farther from the place where the rainwater drainage facility 80 is installed. Therefore, the sixth embodiment is an example in which the flow value estimation unit also calculates a coefficient for the rain intensity value in the segmented region farther from the outflow point in the target region as the elapsed time in the rain event is longer. In other words, when the elapsed time in the rain event is short, the flow value estimation unit in the sixth embodiment calculates a coefficient for the rain intensity value in the segmented area close to the outflow location in the target area. The coefficient to be calculated is increased so as to widen the segment area from the location side.

第6実施形態である流量予測装置の機能構成図は、第1実施形態である流量予測装置60の機能構成を示した図2を用いることとする。第6実施形態においては、第1実施形態と相違する部分について説明し、共通する部分についての説明を省略する。   As the functional configuration diagram of the flow rate prediction apparatus according to the sixth embodiment, FIG. 2 showing the functional configuration of the flow rate prediction apparatus 60 according to the first embodiment is used. In the sixth embodiment, portions that are different from the first embodiment will be described, and descriptions of common portions will be omitted.

流量値推定部64は、一または複数の降雨イベントにおける、降雨強度値データセットおよび流量値データの時系列データを記憶部63から読み込む。そして、流量値推定部64は、係数を算出する際に、降雨イベントにおける経過時間が長いほど、対象領域における流出箇所から離れた区分領域における降雨強度値に対する係数をも算出する。   The flow rate estimation unit 64 reads the rainfall intensity value data set and the time series data of the flow rate value data from one or more rainfall events from the storage unit 63. Then, when calculating the coefficient, the flow rate estimation unit 64 also calculates a coefficient for the rainfall intensity value in the segmented area farther from the outflow location in the target area as the elapsed time in the rain event is longer.

これにより、雨水流量の予測精度を落とすことなく、係数算出の負荷を軽くし、また処理時間を短縮させることができる。   Thereby, the load of coefficient calculation can be reduced and the processing time can be shortened without reducing the prediction accuracy of the rainwater flow rate.

なお、各実施形態である流量予測装置における降雨強度値取得部、流量値取得部、記憶部、流量値推定部、および係数更新部における少なくとも一つの構成を装置として、分離させて流量予測システムとして構成してもよい。   In addition, at least one configuration in the rainfall intensity value acquisition unit, the flow value acquisition unit, the storage unit, the flow value estimation unit, and the coefficient update unit in the flow rate prediction device according to each embodiment is separated as a device and used as a flow rate prediction system. It may be configured.

また、上述した実施形態である流量予測装置60の一部の機能をコンピュータで実現するようにしてもよい。この場合、その機能を実現するための流量予測プログラムをコンピュータ読み取り可能な記録媒体に記録し、この記録媒体に記録された流量予測プログラムをコンピュータシステムに読み込ませて、このコンピュータシステムが実行することによって実現してもよい。なお、このコンピュータシステムとは、オペレーティング・システム(Operating System;OS)や周辺装置のハードウェアを含むものである。また、コンピュータ読み取り可能な記録媒体とは、フレキシブルディスク、光磁気ディスク、光ディスク、メモリカード等の可搬型記録媒体、コンピュータシステムに備えられる磁気ハードディスクやソリッドステートドライブ等の記憶装置のことをいう。さらに、コンピュータ読み取り可能な記録媒体とは、インターネット等のコンピュータネットワーク、および電話回線や携帯電話網を介してプログラムを送信する場合の通信回線のように、短時間の間、動的にプログラムを保持するもの、さらには、その場合のサーバ装置やクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持するものを含んでもよい。また上記の流量予測プログラムは、前述した機能の一部を実現するためのものであってもよく、さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせにより実現するものであってもよい。   Moreover, you may make it implement | achieve a part of function of the flow volume prediction apparatus 60 which is embodiment mentioned above with a computer. In this case, a flow rate prediction program for realizing the function is recorded on a computer-readable recording medium, the flow rate prediction program recorded on the recording medium is read by the computer system, and the computer system executes the program. It may be realized. This computer system includes an operating system (OS) and hardware of peripheral devices. The computer-readable recording medium is a portable recording medium such as a flexible disk, a magneto-optical disk, an optical disk, or a memory card, and a storage device such as a magnetic hard disk or a solid state drive provided in the computer system. Furthermore, a computer-readable recording medium dynamically holds a program for a short time, such as a computer network such as the Internet, and a communication line when transmitting a program via a telephone line or a cellular phone network. In addition, a server that holds a program for a certain period of time, such as a volatile memory inside a computer system serving as a server device or client in that case, may be included. Further, the above flow rate prediction program may be for realizing a part of the above-described functions, and further, realizing the above-described functions by a combination with a program already recorded in the computer system. May be.

以上、詳述したとおり、各実施形態における流量予測装置によれば、対象領域において降雨イベントについて一定時間ごとに計測した降雨強度値および流量値に基づいて、未来時刻における当該対象領域からの雨水の流量を推定する流量推定部を有していることにより、エンジニアリング作業量を少なく抑えた簡便な方法によって、雨水流量の予測を精度よく行うことができる。   As described above in detail, according to the flow rate prediction device in each embodiment, rainwater from the target region at a future time based on the rainfall intensity value and the flow rate value measured at regular intervals for a rain event in the target region. By having the flow rate estimation unit for estimating the flow rate, the rainwater flow rate can be predicted with high accuracy by a simple method with a small amount of engineering work.

以上、本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。   As mentioned above, although some embodiment of this invention was described, these embodiment is shown as an example and is not intending limiting the range of invention. These embodiments can be implemented in various other forms, and various omissions, replacements, and changes can be made without departing from the spirit of the invention. These embodiments and their modifications are included in the scope and gist of the invention, and are also included in the invention described in the claims and the equivalents thereof.

10…雨水管、20…排水ポンプ、30…降雨レーダ、40…信号処理装置、50…配信装置、60…流量予測装置、61…降雨強度値取得部、62…流量値取得部、63…記憶部、64…流量値推定部、65…係数更新部、70…流量測定装置、80…雨水排水施設   DESCRIPTION OF SYMBOLS 10 ... Rainwater pipe, 20 ... Drain pump, 30 ... Rainfall radar, 40 ... Signal processing device, 50 ... Distribution device, 60 ... Flow rate prediction device, 61 ... Rainfall intensity value acquisition unit, 62 ... Flow rate value acquisition unit, 63 ... Storage unit 64 ... flow rate estimation unit, 65 ... coefficient updating unit, 70 ... flow rate measuring device, 80 ... rain drainage facility

Claims (9)

対象領域を区分した複数の区分領域それぞれにおいて、降雨イベントにおける降雨強度値をあらかじめ一定時間ごとに取得する降雨強度値取得部と、
前記対象領域から流出する雨水の流量値を、前記降雨強度値取得部による取得のタイミングに同期させて取得する流量値取得部と、
前記降雨強度値取得部が取得した降雨強度値と前記流量値取得部が取得した流量値とに基づいて、前記対象領域から流出する雨水の未来時刻における流量値を推定する流量値推定部と、
を備える流量予測装置。
In each of a plurality of divided areas that divide the target area, a rainfall intensity value acquisition unit that acquires a rainfall intensity value in a rain event every predetermined time in advance,
A flow value acquisition unit that acquires a flow value of rainwater flowing out of the target area in synchronization with an acquisition timing by the rainfall intensity value acquisition unit;
Based on the rainfall intensity value acquired by the rainfall intensity value acquisition unit and the flow rate value acquired by the flow rate acquisition unit, a flow value estimation unit that estimates a flow value at a future time of rainwater flowing out from the target area;
A flow rate prediction apparatus comprising:
前記流量値推定部は、重回帰分析処理を実行することによって前記複数の区分領域における降雨強度値に対する流量値の予測式の係数を算出し、前記係数を適用した前記予測式を計算することによって前記未来時刻における流量値を推定する、
請求項1記載の流量予測装置。
The flow value estimation unit calculates a coefficient of a prediction formula of a flow value with respect to a rainfall intensity value in the plurality of divided areas by executing a multiple regression analysis process, and calculates the prediction formula to which the coefficient is applied. Estimating a flow value at the future time;
The flow rate prediction apparatus according to claim 1.
前記降雨強度値取得部は、取得した降雨強度値を、降雨強度に応じて区分した複数のレンジに対応させて正規化する、
請求項1または請求項2いずれか記載の流量予測装置。
The rainfall intensity value acquisition unit normalizes the acquired rainfall intensity value corresponding to a plurality of ranges divided according to the rainfall intensity,
The flow rate prediction apparatus according to claim 1 or 2.
前記流量値推定部は、降雨強度に応じて区分した複数のレンジそれぞれについて前記予測式における係数を算出し、レンジごとに前記予測式を計算する、
請求項2または請求項3いずれか記載の流量予測装置。
The flow rate estimation unit calculates a coefficient in the prediction formula for each of a plurality of ranges divided according to rainfall intensity, and calculates the prediction formula for each range.
The flow volume prediction apparatus in any one of Claim 2 or Claim 3.
前記流量値推定部は、前記降雨イベントにおける経過時間が長いほど、前記対象領域における流出箇所から離れた区分領域における降雨強度値に対する係数をも算出する、
請求項2から請求項4いずれか一項記載の流量予測装置。
The flow rate estimation unit also calculates a coefficient for a rain intensity value in a segmented area away from an outflow location in the target area, as the elapsed time in the rain event is longer.
The flow volume prediction apparatus as described in any one of Claims 2-4.
前記流量値推定部が算出する係数を所定のタイミングで更新する係数更新部、
をさらに備える請求項2から請求項5いずれか一項記載の流量予測装置。
A coefficient updating unit that updates the coefficient calculated by the flow rate estimation unit at a predetermined timing;
The flow rate prediction apparatus according to claim 2, further comprising:
降雨強度値取得部が、対象領域を区分した複数の区分領域それぞれにおいて、降雨イベントにおける降雨強度値をあらかじめ一定時間ごとに取得する降雨強度値取得ステップと、
流量値取得部が、前記対象領域から流出する雨水の流量値を、前記降雨強度値取得ステップにおける取得のタイミングに同期させて取得する流量値取得ステップと、
前記降雨強度値取得ステップにおいて取得された降雨強度値と前記流量値取得ステップにおいて取得された流量値とに基づいて、流量値推定部が、前記対象領域から流出する雨水の未来時刻における流量値を推定する流量値推定ステップと、
を有する流量予測方法。
A rainfall intensity value acquisition step in which a rainfall intensity value acquisition unit acquires a rainfall intensity value in a rain event in advance every predetermined time in each of a plurality of divided areas into which the target area has been divided,
A flow value acquisition unit acquires a flow value of rainwater flowing out of the target area in synchronization with an acquisition timing in the rainfall intensity value acquisition step, and
Based on the rainfall intensity value acquired in the rainfall intensity value acquisition step and the flow rate value acquired in the flow value acquisition step, the flow value estimation unit calculates a flow value at a future time of rainwater flowing out from the target area. A flow value estimation step to be estimated;
A flow rate prediction method comprising:
コンピュータを、
対象領域を区分した複数の区分領域それぞれにおいて、降雨イベントにおける降雨強度値をあらかじめ一定時間ごとに取得する降雨強度値取得部と、
前記対象領域から流出する雨水の流量値を、前記降雨強度値取得部による取得のタイミングに同期させて取得する流量値取得部と、
前記降雨強度値取得部が取得した降雨強度値と前記流量値取得部が取得した流量値とに基づいて、前記対象領域から流出する雨水の未来時刻における流量値を推定する流量値推定部と、
として機能させるための流量予測プログラム。
Computer
In each of a plurality of divided areas that divide the target area, a rainfall intensity value acquisition unit that acquires a rainfall intensity value in a rain event every predetermined time in advance,
A flow value acquisition unit that acquires a flow value of rainwater flowing out of the target area in synchronization with an acquisition timing by the rainfall intensity value acquisition unit;
Based on the rainfall intensity value acquired by the rainfall intensity value acquisition unit and the flow rate value acquired by the flow rate acquisition unit, a flow value estimation unit that estimates a flow value at a future time of rainwater flowing out from the target area;
Flow prediction program to function as.
対象領域を区分した複数の区分領域それぞれにおいて、降雨イベントにおける降雨強度値をあらかじめ一定時間ごとに取得する降雨強度値取得装置と、
前記対象領域から流出する雨水の流量値を、前記降雨強度値取得装置による取得のタイミングに同期させて取得する流量値取得装置と、
前記降雨強度値取得装置が取得した降雨強度値と前記流量値取得装置が取得した流量値とに基づいて、前記対象領域から流出する雨水の未来時刻における流量値を推定する流量値推定装置と、
を備える流量予測システム。
In each of a plurality of divided areas into which the target area is divided, a rainfall intensity value acquisition device that acquires a rainfall intensity value in a rain event every predetermined time in advance,
A flow value acquisition device that acquires the flow value of rainwater flowing out of the target area in synchronization with the acquisition timing by the rainfall intensity value acquisition device;
A flow value estimation device that estimates a flow value at a future time of rainwater flowing out from the target area, based on a rainfall intensity value acquired by the rainfall intensity value acquisition device and a flow value acquired by the flow value acquisition device;
A flow rate prediction system comprising:
JP2013117920A 2013-06-04 2013-06-04 Flow prediction device, flow prediction method, flow prediction program, and flow prediction system Active JP6222689B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2013117920A JP6222689B2 (en) 2013-06-04 2013-06-04 Flow prediction device, flow prediction method, flow prediction program, and flow prediction system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2013117920A JP6222689B2 (en) 2013-06-04 2013-06-04 Flow prediction device, flow prediction method, flow prediction program, and flow prediction system

Publications (2)

Publication Number Publication Date
JP2014234675A true JP2014234675A (en) 2014-12-15
JP6222689B2 JP6222689B2 (en) 2017-11-01

Family

ID=52137584

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013117920A Active JP6222689B2 (en) 2013-06-04 2013-06-04 Flow prediction device, flow prediction method, flow prediction program, and flow prediction system

Country Status (1)

Country Link
JP (1) JP6222689B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116485034A (en) * 2023-05-15 2023-07-25 深圳市嘉闰州生态建工有限公司 Urban drainage prediction method and system

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05134056A (en) * 1991-11-13 1993-05-28 Toshiba Corp Estimating apparatus for inflow of rain water
JPH06102911A (en) * 1992-09-21 1994-04-15 Hitachi Ltd Method and device for predicting inflow water rate in drainage pump place
JPH0961545A (en) * 1995-08-21 1997-03-07 Toshiba Corp Device for predicting flowout of rainfall
JP2000345604A (en) * 1999-06-08 2000-12-12 Yaskawa Electric Corp Device for estimating inflow sewage quantity
JP2002285634A (en) * 2001-03-27 2002-10-03 Fuji Electric Co Ltd Device for estimating inflow of rainwater
WO2004038343A1 (en) * 2002-10-25 2004-05-06 Waste Water Control Aps Determining a flow in a wastewater system
JP2006184206A (en) * 2004-12-28 2006-07-13 Mitsui Kyodo Kensetsu Consultant Kk System and program for predicting distributed outflow
JP2006266072A (en) * 2004-12-14 2006-10-05 Toshiba Corp Rainwater drain support system and method, its control system and method
JP2006295473A (en) * 2005-04-08 2006-10-26 Nippon Telegr & Teleph Corp <Ntt> Apparatus, method and program for detecting electronic watermark
JP2007205001A (en) * 2006-02-01 2007-08-16 Fuji Electric Systems Co Ltd Discharge forecasting apparatus
JP2008106477A (en) * 2006-10-24 2008-05-08 Yamatake Corp Rainwater outflow estimating equipment and program
JP2010084493A (en) * 2008-10-03 2010-04-15 Yaskawa Electric Corp Drainage operation support device in combined pump station and operation method thereof

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05134056A (en) * 1991-11-13 1993-05-28 Toshiba Corp Estimating apparatus for inflow of rain water
JPH06102911A (en) * 1992-09-21 1994-04-15 Hitachi Ltd Method and device for predicting inflow water rate in drainage pump place
JPH0961545A (en) * 1995-08-21 1997-03-07 Toshiba Corp Device for predicting flowout of rainfall
JP2000345604A (en) * 1999-06-08 2000-12-12 Yaskawa Electric Corp Device for estimating inflow sewage quantity
JP2002285634A (en) * 2001-03-27 2002-10-03 Fuji Electric Co Ltd Device for estimating inflow of rainwater
WO2004038343A1 (en) * 2002-10-25 2004-05-06 Waste Water Control Aps Determining a flow in a wastewater system
JP2006266072A (en) * 2004-12-14 2006-10-05 Toshiba Corp Rainwater drain support system and method, its control system and method
JP2006184206A (en) * 2004-12-28 2006-07-13 Mitsui Kyodo Kensetsu Consultant Kk System and program for predicting distributed outflow
JP2006295473A (en) * 2005-04-08 2006-10-26 Nippon Telegr & Teleph Corp <Ntt> Apparatus, method and program for detecting electronic watermark
JP2007205001A (en) * 2006-02-01 2007-08-16 Fuji Electric Systems Co Ltd Discharge forecasting apparatus
JP2008106477A (en) * 2006-10-24 2008-05-08 Yamatake Corp Rainwater outflow estimating equipment and program
JP2010084493A (en) * 2008-10-03 2010-04-15 Yaskawa Electric Corp Drainage operation support device in combined pump station and operation method thereof

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116485034A (en) * 2023-05-15 2023-07-25 深圳市嘉闰州生态建工有限公司 Urban drainage prediction method and system
CN116485034B (en) * 2023-05-15 2023-09-08 深圳市嘉闰州生态建工有限公司 Urban drainage prediction method and system

Also Published As

Publication number Publication date
JP6222689B2 (en) 2017-11-01

Similar Documents

Publication Publication Date Title
CN101876546B (en) MEMS (Micro Electronic Mechanical System) gyro data processing method based on wavelet threshold de-noising and FAR (Finite Automaton Recognizable) model
KR101829560B1 (en) Method, system and computer device for predicting a capacity based on kalman filter
JP6189254B2 (en) Rainwater drainage pump control device, rainwater drainage pump control method, rainwater drainage pump control program, and parameter providing device
CN108416468B (en) Flood advanced early warning and forecasting method
CN107526880B (en) Method for calculating upstream river-entering pollution flux based on river water quality observation
JP6880841B2 (en) Photovoltaic power generation output prediction device considering snow cover
CN113051529B (en) Localized uniform weight particle filtering data assimilation method based on statistical observation
JP6222689B2 (en) Flow prediction device, flow prediction method, flow prediction program, and flow prediction system
CN105488352A (en) Long-term deflection monitoring data based concrete bridge stiffness reliability evaluation method
KR101456524B1 (en) Rainfall information system and rainfall information server using data from rain sensor of vehicle, rainfall meter and radar rainfall, with big data process
JP2018148743A (en) Photovoltaic power output prediction device considering residual snow
CN109902266B (en) Riverway flow calculation method based on Copula function
CN104407366A (en) Pseudo-range smooth processing method
WO2016098250A1 (en) Waveform estimation device and waveform estimation method
CN106153046B (en) Gyro random noise AR modeling method based on self-adaptive Kalman filtering
JP2014211330A (en) Target tracking device and target tracking method
JP7067306B2 (en) Snow cover prediction device, snow cover prediction method
CN110941032A (en) Typhoon forecasting method, device, equipment and computer-readable storage medium
WO2019207622A1 (en) Power demand prediction device, power demand prediction method, and program therefor
JP6116398B2 (en) Waveform estimation apparatus and waveform estimation method
JP6290072B2 (en) Unknown water generation area estimation device, unknown water generation area estimation method, and computer program
CN101509774A (en) ARMA time-series north-searching method based on optical fiber gyroscope
Yi et al. Cycle slip detection and correction of GPS carrier phase based on wavelet transform and neural network
JP5383767B2 (en) Inflow amount prediction apparatus, inflow amount prediction method, and program
CN103761442B (en) Forecasting method and device for flow parameters of micro areas

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20160308

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20161212

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20170110

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170313

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: 20170829

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20170911

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20170911

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20170927

R150 Certificate of patent or registration of utility model

Ref document number: 6222689

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150