JP2006002401A - Sewage inflow quantity predicting device - Google Patents

Sewage inflow quantity predicting device Download PDF

Info

Publication number
JP2006002401A
JP2006002401A JP2004178618A JP2004178618A JP2006002401A JP 2006002401 A JP2006002401 A JP 2006002401A JP 2004178618 A JP2004178618 A JP 2004178618A JP 2004178618 A JP2004178618 A JP 2004178618A JP 2006002401 A JP2006002401 A JP 2006002401A
Authority
JP
Japan
Prior art keywords
data
inflow
prediction
sewage
sewage inflow
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
JP2004178618A
Other languages
Japanese (ja)
Other versions
JP4326408B2 (en
Inventor
Tatsuya Iizaka
達也 飯坂
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.)
Fuji Electric Co Ltd
Original Assignee
Fuji Electric Systems Co Ltd
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 Fuji Electric Systems Co Ltd filed Critical Fuji Electric Systems Co Ltd
Priority to JP2004178618A priority Critical patent/JP4326408B2/en
Publication of JP2006002401A publication Critical patent/JP2006002401A/en
Application granted granted Critical
Publication of JP4326408B2 publication Critical patent/JP4326408B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

<P>PROBLEM TO BE SOLVED: To provide a sewage inflow quantity predicting device for improving accuracy in predicting a sewage inflow quantity, by accurately predicting a sewage inflow quantity and a rainwater inflow quantity by an uncomplicated sewage inflow quantity predicting model and a rainwater inflow quantity predicting model which utilize a general technology. <P>SOLUTION: This sewage inflow quantity predicting device forms sewage inflow quantity predicting data, by adding up formed sewage inflow quantity predicting data and rainwater inflow quantity predicting data, by forming the rainwater inflow quantity predicting data by inputting meteorological data to the rainwater inflow quantity predicting model, by forming the sewage inflow quantity predicting data by inputting calendar data to the sewage inflow quantity predicting model, by separately constructing the sewage inflow quantity predicting model and the rainwater inflow quantity predicting model, by using rainwater inflow quantity data and sewage inflow quantity data separately formed from sewage inflow quantity data being a sewage inflow quantity to a pump well. <P>COPYRIGHT: (C)2006,JPO&NCIPI

Description

本発明は、下水ポンプ場または終末処理場のポンプ井へ流入する下水の下水流入量を予測する下水流入量予測装置に関する。   The present invention relates to a sewage inflow prediction device for predicting a sewage inflow amount of sewage flowing into a pump well of a sewage pump station or a final treatment plant.

終末処理場において効率よく下水を処理するためには、1日の下水流入量の変動が少ないことが望ましい。このように下水流入量の変動を少なくするために、現状では下水道の中間または終末に多数配置されるポンプ井に流入する下水流入量を予測し、ポンプの起動停止を制御することで、終末処理場に流れ込む下水流入量の変動を小さくしている。   In order to efficiently treat sewage at the final treatment plant, it is desirable that the fluctuation of the daily sewage inflow is small. In order to reduce fluctuations in the amount of sewage inflow in this way, at present, the sewage inflow amount that flows into the pump wells that are placed in the middle or at the end of the sewage system is predicted, and the start / stop of the pump is controlled, so that The fluctuation of the amount of inflow of sewage flowing into the field is reduced.

従来技術による下水流入量を予測する方法として、主にRRL法(ロード・リサーチ・ラボラトリ法)が用いられている。また、近年ではRRL法の他にも下水流入量予測手法が開発されており、例えば特許文献1〜3が開示されている。   The RRL method (Road Research Laboratory Method) is mainly used as a method for predicting the amount of sewage inflow according to the prior art. In recent years, in addition to the RRL method, a sewage inflow prediction method has been developed. For example, Patent Documents 1 to 3 are disclosed.

特許文献1(特開平5−134715号公報、発明の名称「ニューラルネットワーク応用雨水流入量予測装置」)に記載された従来技術は、レーダー雨量計とニューラルネットワークとを用いることを特徴とした雨水流入量予測であり、ニューラルネットワークの学習機能により経年変化に対応でき、また、レーダー雨量計を利用することで、1地点の雨量計情報だけでなく、面的な雨量情報を使用できるようにしている。   The conventional technique described in Patent Document 1 (Japanese Patent Laid-Open No. Hei 5-134715, “Neural Network Applied Rainwater Inflow Prediction Device”) uses a radar rain gauge and a neural network. It is an amount prediction, can cope with aging by the learning function of the neural network, and by using a radar rain gauge, it can use not only the rain gauge information at one point but also the area rainfall information .

また、特許文献2(特開2000−345604号公報、発明の名称「流入下水量予測装置」)に記載された従来技術は、合流式下水流入量が予測できる手法である。つまり、汚水流入量と雨水流入量との区別無く予測するというものである。
一般的に汚水流入量は、生活排水によるものであるため、気象に関係なく時間により流入量が変化するという特徴がある。例えば、図10の通常時の汚水流入量−時間特性図で示すように、夕食時から就寝前までにかけて増加した汚水流入量が夜間は減少し、朝になるとまた汚水流入量が増大するというように、朝と夜とに流入ピークを持つというような特徴である。一方、雨水流入量は当然に降雨量により変化する。
特許文献2の従来技術では、これらのような特徴を持つ雨水流入量と汚水流入量との予測手法を異ならせており、汚水流入量は過去の平均流入量により算出し、雨水流入量は自己回帰モデルにより算出している。
Moreover, the prior art described in patent document 2 (Unexamined-Japanese-Patent No. 2000-345604, the name of an invention "inflow sewage amount prediction apparatus") is a method which can predict a combined sewage inflow amount. That is, the prediction is made without distinguishing between the sewage inflow and the rainwater inflow.
In general, the amount of inflow of sewage is due to domestic wastewater, so that the amount of inflow changes with time regardless of the weather. For example, as shown in the normal sewage inflow-time characteristic diagram of FIG. 10, the sewage inflow increased from dinner to before bedtime decreases at night, and the sewage inflow increases again in the morning. In addition, it is characterized by having an inflow peak in the morning and at night. On the other hand, the amount of rainwater inflow naturally varies depending on the amount of rainfall.
In the prior art of Patent Document 2, the method for predicting the amount of stormwater inflow and the amount of sewage inflow with these characteristics is different, and the amount of sewage inflow is calculated from the past average inflow. Calculated using a regression model.

また、特許文献3(特開2003−27567号公報、発明の名称「下水流入量予測装置および方法、サーバ装置」)に記載された従来技術は、TCBM(Topological Case-Based Modelling)による予測方式である。   The prior art described in Patent Document 3 (Japanese Patent Laid-Open No. 2003-27567, “Sewage Inflow Prediction Device and Method, Server Device”) is a prediction method based on TCBM (Topological Case-Based Modeling). is there.

特開平5−134715号公報Japanese Patent Laid-Open No. 5-134715 特開2000−345604号公報JP 2000-345604 A 特開2003−27567号公報JP 2003-27567 A

従来技術であるRRL法は、精度の高い予測を行なうためにはパラメータ設定のため膨大な計算が必要であるなど予測式構築に膨大な作業が必要であり、また経年変化に対応できないという問題点があった。つまり、人口、舗装率などの変化により予測式がずれてきた場合には、予測式修正のため再度膨大な作業を要するという問題がある。また、予測可能なのは分流式下水道の雨水流入量のみであり、汚水流入量や合流式下水道の流入量は予測できないという問題もあった。   The RRL method, which is a conventional technology, requires a huge amount of work to construct a prediction formula, such as a large amount of calculation for parameter setting in order to make a highly accurate prediction, and cannot cope with aging. was there. In other words, when the prediction formula is shifted due to changes in the population, the pavement rate, etc., there is a problem that a huge amount of work is required again to correct the prediction formula. In addition, only the rainwater inflow of the diversion sewerage can be predicted, and the inflow of sewage and the inflow of the combined sewerage cannot be predicted.

特許文献1に記載の従来技術は、ニューラルネットワークの学習機能により人口、舗装率が変化しても自動的に予測式が修正される利点があるが、雨水流入量しか予測できず、汚水流入量や合流式下水道の流入量は予測できないという問題があった。   The prior art described in Patent Document 1 has an advantage that the prediction formula is automatically corrected even if the population and the pavement rate change by the learning function of the neural network, but only the rainwater inflow can be predicted, and the sewage inflow In addition, there was a problem that the inflow of the combined sewerage system could not be predicted.

特許文献2に記載の従来技術は、人口、舗装率の変化に自動的に対応できる利点がある。しかし、汚水流入量は単純な過去平均値であるため、特異データがある場合には適切な汚水平均値が算出できない。
例えば、下水道管内の汚れ清掃を行なう場合、下水ポンプ場内の汚水を貯め、一気に送水して下水管内に付着した汚れを取る方法を採用している。このため図11の特異時の汚水流入量−時間特性図のように、清掃前では汚水流入量が無くなり、また清掃時に汚水流入量が尖塔状になり、通常の特性図と大きく異なる。このような下水管特有の清掃方法のため、特許文献2の発明による汚水予測方法では、適切な平均化を行なうことができないという問題があった。
また、雨水流入量は自己回帰式により予測しているため、雨という気象予報があったとしても予測値に反映させることができず、未来の気象状態を考慮することができないという問題もあった。
The prior art described in Patent Document 2 has an advantage that it can automatically cope with changes in population and pavement rate. However, since the sewage inflow is a simple past average value, an appropriate sewage average value cannot be calculated when there is singular data.
For example, when cleaning dirt in the sewer pipe, a method is adopted in which the dirty water in the sewage pumping station is stored, and the dirt adhering to the sewage pipe is removed at a stretch. For this reason, as shown in the sewage inflow amount-time characteristic chart in FIG. 11, the sewage inflow quantity disappears before cleaning, and the sewage inflow quantity becomes a spire shape at the time of cleaning, which is greatly different from the normal characteristic chart. Due to such a cleaning method peculiar to the sewage pipe, the sewage prediction method according to the invention of Patent Document 2 has a problem that appropriate averaging cannot be performed.
In addition, since the amount of rainwater inflow is predicted by an autoregressive equation, there is a problem that even if there is a weather forecast of rain, it cannot be reflected in the predicted value and the future weather conditions cannot be taken into account. .

特許文献3に記載の従来技術は、上記発明の問題点を解決しているが、TCBMという特殊な手法を用いているため、一般開発者は容易に適用することができないものであり、周知技術を活用して予測したいという要請があった。   The prior art described in Patent Document 3 solves the problems of the invention described above, but since a special technique called TCBM is used, general developers cannot easily apply it. There was a request to make predictions using.

まとめると、
(1)経年変化に対応した予測、
(2)合流式下水道の予測、
(3)特異データ(下水管清掃時など)の悪影響のない予測、
(4)雨量等気象情報を用いた予測、
(5)特殊な方法を用いない周知技術を活用しての予測、
という要件を全て満たすような下水流入量予測装置が求められていた。
Summary,
(1) Forecast corresponding to secular change,
(2) Forecast of combined sewerage,
(3) Prediction with no adverse effect on specific data (such as when cleaning sewer pipes)
(4) Prediction using weather information such as rainfall,
(5) Prediction using well-known technology that does not use special methods,
There was a need for a sewage inflow prediction device that would satisfy all the requirements.

本発明は、これら従来技術の問題点を解決するためになされたものであり、その目的は、一般技術を活用した複雑でない汚水流入量予測モデルと雨水流入量予測モデルとにより汚水流入量と雨水流入量とを正確に予測し、下水流入量の予測精度の向上を実現する下水流入量予測装置を提供することにある。   The present invention has been made to solve these problems of the prior art. The purpose of the present invention is to solve the problems of sewage inflow and rainwater by using an uncomplicated sewage inflow prediction model and a stormwater inflow prediction model using general techniques. An object of the present invention is to provide a sewage inflow amount prediction device that accurately predicts the inflow amount and improves the prediction accuracy of the sewage inflow amount.

本発明の請求項1に係る下水流入量予測装置は、
下水ポンプ場または終末処理場のポンプ施設へ流入する下水の流入量を予測する下水流入量予測装置において、
気象実績及び気象予報についての気象データを収集する気象データ収集手段と、
暦についての暦データを収集する暦データ収集手段と、
ポンプ井内の下水の水位についての水位データを収集する水位データ収集手段と、
ポンプ井の下水のポンプ揚水量(送水量)についての揚水量データを収集するポンプ揚水量データ収集手段と、
暦データ、気象データ、水位データおよびポンプ揚水量データを蓄積するデータ蓄積手段と、
暦データ、水位データおよび揚水量データに基づいてポンプ井への下水流入量である下水流入量データを算出する下水流入量データ算出手段と、
下水流入量データを分離して汚水流入量データと雨水流入量データとを算出する分離手段と、
暦データを入力因子とし、また、汚水流入量データを出力因子とするような汚水流入量予測モデルを構築する汚水流入量予測モデル構築手段と、
気象データを入力因子とし、また、雨水流入量データを出力因子とするような雨水流入量予測モデルを構築する雨水流入量予測モデル構築手段と、
予測に必要な暦データを汚水流入量予測モデルに入力して汚水流入量予測データを生成する汚水流入量予測手段と、
予測に必要な気象データを雨水流入量予測モデルに入力して雨水流入量予測データを生成する雨水流入量予測手段と、
生成した汚水流入量予測データと雨水流入量予測データとを合算して下水流入量予測データを生成する合算手段と、
を有することを特徴とする。
A sewage inflow prediction device according to claim 1 of the present invention is
In the sewage inflow prediction device that predicts the inflow of sewage flowing into the pump facility at the sewage pump station or the final treatment plant,
Meteorological data collection means for collecting meteorological data on weather results and weather forecasts;
Calendar data collection means for collecting calendar data about the calendar;
Water level data collection means for collecting water level data about the sewage water level in the pump well,
Pump pumping data collection means for collecting pumping volume data about pumping volume (pumped volume) of sewage of the pump well;
Data storage means for storing calendar data, meteorological data, water level data and pump pumping volume data;
Sewage inflow data calculation means for calculating sewage inflow data that is sewage inflow to the pump well based on the calendar data, the water level data, and the pumped water data;
Separation means for separating sewage inflow data and calculating sewage inflow data and rainwater inflow data;
A sewage inflow prediction model construction means for constructing a sewage inflow prediction model with calendar data as an input factor and sewage inflow data as an output factor;
A stormwater inflow prediction model construction means for constructing a stormwater inflow prediction model using weather data as an input factor and stormwater inflow data as an output factor;
A sewage inflow forecasting means for generating sewage inflow forecast data by inputting calendar data necessary for prediction into the sewage inflow forecast model,
A stormwater inflow forecasting means for generating stormwater inflow forecast data by inputting meteorological data necessary for the forecast into the stormwater inflow forecast model;
A summing means for generating the sewage inflow prediction data by adding the generated sewage inflow prediction data and the rainwater inflow prediction data;
It is characterized by having.

また、本発明の請求項2に係る下水流入量予測装置は、
請求項1に記載の下水流入量予測装置において、
前記汚水流入量予測モデルは、
特異データが除去された上で、暦データと汚水流入量データとが連関して登録され、暦データが入力されると対応する汚水流入量データが抽出されることで予測する事例ベース推論による予測モデルであることを特徴とする。
The sewage inflow prediction device according to claim 2 of the present invention is
In the sewage inflow prediction device according to claim 1,
The sewage inflow prediction model is
Prediction based on case-based reasoning, in which calendar data and sewage inflow data are linked and registered after singular data is removed, and the corresponding sewage inflow data is extracted when calendar data is input It is a model.

また、本発明の請求項3に係る下水流入量予測装置は、
請求項1または請求項2に記載の下水流入量予測装置において、
前記雨水流入量予測モデルはニューラルネットワークによる予測モデルであることを特徴とする。
A sewage inflow prediction device according to claim 3 of the present invention is
In the sewage inflow prediction device according to claim 1 or 2,
The rainwater inflow prediction model is a prediction model using a neural network.

また、本発明の請求項4に係る下水流入量予測装置は、
請求項3に記載の下水流入量予測装置において、
前記雨水流入量予測モデルはn時間単位で雨水流量予測データを算出する予測モデルであり、
前記雨水流入量予測手段は、現在の雨水流量データと雨水流量予測データとを、または、隣接する二つの雨水流量予測データを補間することで時間間隔がn時間よりも短い雨水流量予測データを生成する手段であることを特徴とする。
A sewage inflow prediction device according to claim 4 of the present invention is
In the sewage inflow prediction device according to claim 3,
The rainwater inflow prediction model is a prediction model for calculating rainwater flow rate prediction data in units of n hours,
The rainwater inflow prediction means generates rainwater flow prediction data whose time interval is shorter than n hours by interpolating the current rainwater flow data and rainwater flow prediction data or two adjacent rainwater flow prediction data. It is a means to do.

また、本発明の請求項5に係る下水流入量予測装置は、
請求項1〜請求項4の何れか一項に記載の下水流入量予測装置において、
下水流入量予測データを、過去の日時における下水流入量予測データと下水流入量データとから算出した誤差量より補正する補正手段を備えることを特徴とする。
A sewage inflow prediction device according to claim 5 of the present invention is
In the sewage inflow prediction device according to any one of claims 1 to 4,
The present invention is characterized by comprising correction means for correcting the sewage inflow prediction data from an error amount calculated from the sewage inflow prediction data and the sewage inflow data at a past date and time.

また、本発明の請求項6に係る下水流入量予測装置は、
請求項5に記載の下水流入量予測装置において、
前記補正手段は、過去の予測誤差を用いて未来の予測値を補正するときに、近い将来の予測値の補正量を大きく、遠い将来の予測値の補正量を小さくする手段であることを特徴とする。
A sewage inflow prediction device according to claim 6 of the present invention is
In the sewage inflow prediction device according to claim 5,
The correction unit is a unit that increases a correction amount of a near future prediction value and decreases a correction amount of a far future prediction value when correcting a future prediction value using a past prediction error. And

以上のような本発明によれば、一般技術を活用した複雑でない汚水流入量予測モデルと雨水流入量予測モデルとにより汚水流入量と雨水流入量とを正確に予測し、下水流入量の予測精度の向上を実現する下水流入量予測装置を提供することができる。   According to the present invention as described above, a sewage inflow amount and a stormwater inflow amount are accurately predicted by an uncomplicated sewage inflow amount prediction model and a stormwater inflow amount prediction model utilizing general techniques, and a prediction accuracy of a sewage inflow amount is obtained. It is possible to provide a sewage inflow amount prediction device that realizes improvement of the above.

以下、本発明を実施するための最良の形態の下水流入量予測装置について図を参照しつつ説明する。図1は、本形態の下水流入量予測装置の構成図である。
下水流入量予測装置1は、入力手段10、出力手段20、記録媒体読書手段30、データ発信手段40、データ収集手段50、データ蓄積手段60、下水流入量データ算出手段70、分離手段80、モデル構築手段90、予測手段100を備えている。
Hereinafter, the best mode sewage inflow prediction device for carrying out the present invention will be described with reference to the drawings. FIG. 1 is a configuration diagram of a sewage inflow prediction device according to this embodiment.
The sewage inflow prediction apparatus 1 includes an input means 10, an output means 20, a recording medium reading means 30, a data transmission means 40, a data collection means 50, a data storage means 60, a sewage inflow data calculation means 70, a separation means 80, a model. The construction means 90 and the prediction means 100 are provided.

入力手段10は、手入力を行うキーボードである。
出力手段20は、ディスプレイやプリンタなどであり、各種データ表示・印刷も行う。
記憶媒体読書手段30は、FD(Flexible Disc)、MO(Magnet Optical Disc)などの記憶媒体に対して情報の読み書きを行う装置である。
データ発信手段40は、ネットワーク2と接続するための手段であり、データ蓄積手段60からデータ発信手段40・ネットワーク2を経由して、外部へ情報の出力を行うことができるようになされている。
The input means 10 is a keyboard for performing manual input.
The output unit 20 is a display, a printer, or the like, and also performs various data display / printing.
The storage medium reading means 30 is a device that reads and writes information from and to a storage medium such as an FD (Flexible Disc) or an MO (Magnet Optical Disc).
The data transmission means 40 is a means for connecting to the network 2 and can output information from the data storage means 60 via the data transmission means 40 and the network 2 to the outside.

データ収集手段50は、各下水ポンプ場や気象情報会社などからネットワーク2を介して各種データを取得する手段であり、さらに気象データ収集手段51、暦データ収集手段52、水位データ収集手段53、ポンプ揚水量データ収集手段54を備える。
気象データ収集手段51は、気象実績および気象予報についての気象データを収集する手段である。ここでいう気象データとは、例えば天候(晴、雨、曇、雪)、雨量、気温、湿度、日照量等の全てまたは何れかを指している。
気象実績とは過去の記録であり、気象実績についての気象データは、各下水ポンプ場を経由せずに気象情報サービス会社からネットワーク2を介して配信される気象データや、また、各下水ポンプ場に独自で設置した雨量計・湿度計などの計測装置・センサから収集整理して各下水ポンプ場を経由して入力された気象データである。
気象予報とは未来の予報であり、気象予報についての気象データは、気象情報サービス会社から気象予報を受信して得たデータである。
このように気象データ収集手段51では、予測対象の下水ポンプ場の気象実績に係る気象データ、予測対象の下水ポンプ場の気象予報に係る気象データを収集する。
The data collection means 50 is a means for acquiring various data from each sewage pump station, a weather information company, etc. via the network 2, and further includes a weather data collection means 51, a calendar data collection means 52, a water level data collection means 53, a pump Pumped water data collection means 54 is provided.
The meteorological data collecting unit 51 is a unit that collects meteorological data on the weather results and the weather forecast. The weather data here refers to all or any of, for example, weather (sunny, rain, cloudy, snow), rainfall, temperature, humidity, and amount of sunlight.
The meteorological results are past records, and meteorological data about the meteorological results are met by the meteorological data distributed from the weather information service company via the network 2 without going through each sewage pump station, and each sewage pump station The weather data is collected and arranged from measuring devices and sensors such as rain gauges and hygrometers, which are installed independently, and input via each sewage pump station.
The weather forecast is a future forecast, and the weather data about the weather forecast is data obtained by receiving a weather forecast from a weather information service company.
Thus, the meteorological data collection means 51 collects meteorological data related to the weather results of the sewage pump station to be predicted and meteorological data related to the weather forecast of the sewage pump station to be predicted.

暦データ収集手段52は、暦に関するデータを収集する手段である。
暦データは年月日時分などの日時データ、その日時の特徴を表す平日・休日・曜日・季節(春夏秋冬)などの特徴データである。
なお、暦データは入力手段10のキーボードなどの装置にて手入力しても良いし、計算により算出しても良い。
The calendar data collecting means 52 is a means for collecting data related to the calendar.
The calendar data is date / time data such as year / month / day / hour / minute and the like, and feature data such as weekdays / holidays / day of the week / season (spring / summer / autumn / winter) representing the characteristics of the date / time.
The calendar data may be manually input by a device such as a keyboard of the input unit 10 or may be calculated by calculation.

水位データ収集手段53は、水位データを収集する手段である。
水位データは、下水ポンプ場にあるポンプ井内の下水の水位を表すデータである。
水位データはネットワーク2を介して複数下水ポンプ場のポンプ井にそれぞれ設置された多数の水位計から所定期間毎または常時入力される。従って水位データには下水ポンプ場を識別する識別データも付加される。
The water level data collecting means 53 is means for collecting water level data.
The water level data is data representing the sewage water level in the pump well in the sewage pump station.
Water level data is input via a network 2 from a plurality of water level gauges installed in pump wells of a plurality of sewage pump stations at predetermined intervals or constantly. Accordingly, identification data for identifying the sewage pump station is also added to the water level data.

ポンプ揚水量データ収集手段54は、ポンプ揚水量データを収集する手段である。
ポンプ揚水量データは、下水ポンプ場にあるポンプによるポンプ井からの揚水量を表すデータである。ポンプ揚水量データはネットワーク2を介して複数下水ポンプ場のポンプに設けられた流量計から所定期間毎または常時入力される。この場合も下水ポンプ場を識別する識別データも付加される。
これら気象データ、暦データ、水位データおよび揚水量データはデータ蓄積手段60に保存し、日々更新してこれらデータを蓄積していく。また、必要時に何時でも取り出せるようにする。このデータ蓄積手段60では計測日時・曜日・季節が特定できる暦データを主キーとして気象データ、水位データおよび揚水量データが関連付けられて登録されてデータベース化されている。
The pump yield data collecting means 54 is means for collecting pump yield data.
Pump yield data is data representing the yield from the pump well by the pumps in the sewage pump station. Pump pumping amount data is inputted via a network 2 from a flow meter provided in a pump of a plurality of sewage pump stations every predetermined period or constantly. In this case, identification data for identifying the sewage pump station is also added.
These meteorological data, calendar data, water level data, and pumped-up amount data are stored in the data storage means 60 and updated daily to store these data. Also, it can be taken out whenever necessary. In this data storage means 60, weather data, water level data, and pumped water amount data are associated and registered in a database using calendar data that can specify the measurement date, day of the week, and season as a main key.

下水流入量データ算出手段70は、暦データから検索して各時間の水位データおよび揚水量データをデータ蓄積手段60から読み出し、これらデータに基づいて下水ポンプ場のポンプ井への下水流入量である下水流入量データを算出する。   The sewage inflow data calculating means 70 retrieves the water level data and the pumped amount data for each hour from the calendar data by retrieving from the calendar data, and is the sewage inflow quantity to the pump well of the sewage pump station based on these data. Calculate sewage inflow data.

下水ポンプ場のポンプ井では、例えば、図2のポンプ井の説明図で示すように、通常は水位計が設置されているが、流入量計は設置されていることが少なく、水位からポンプ井内の下水流入量を算出する。
ポンプ井内のある時点の下水量をM1、同じくある時点より一定時間前の下水量をM2、一定時間にポンプ井から揚水して送水したポンプ揚水量をP1とすると、下水流入量は以下のようになる。
In the pump well of the sewage pump station, for example, as shown in the explanatory diagram of the pump well in FIG. 2, a water level gauge is usually installed, but an inflow meter is rarely installed, Calculate the amount of sewage inflow.
If the amount of sewage at a certain point in the pump well is M1, the amount of sewage before a certain time is M2, and the amount of pumped water pumped from the pump well at a certain time is P1, the amount of sewage inflow is as follows: become.

Figure 2006002401
Figure 2006002401

ここに下水量M1はある時点におけるポンプ井の水位h1から算出した量、下水量M2はある時点からさらに一定時間前のポンプ井の水位h2から算出した量である。この下水量M1,M2の算出式は、ポンプ井の水位を入力因子とする回帰式が用いられることが多い。最も簡単な例では単純比例関係であり、算出式は、M1=a×h1,M2=a×h2となる。
なお、ポンプ井によっては複雑な算出式が必要になる場合もあるが、何れの場合でもある時点での水位を入力因子とする算出式である。
Here, the sewage amount M1 is an amount calculated from the water level h1 of the pump well at a certain time, and the sewage amount M2 is an amount calculated from the water level h2 of the pump well at a certain time before a certain time. As a calculation formula for the sewage amounts M1 and M2, a regression formula using the water level of the pump well as an input factor is often used. The simplest example is a simple proportional relationship, and the calculation formula is M1 = a × h1, M2 = a × h2.
Depending on the pump well, a complicated calculation formula may be required. In any case, the calculation formula uses the water level at a certain point in time as an input factor.

なお、下水流入量の中には、上流にある下水ポンプ場からの揚水量(送水量)も含まれており、これは人為的な操作により変化するものであるため、次式のように上流ポンプ揚水量(送水量)を除いた流入量としても良い。上流ポンプ揚水量P2は、一定時間で上流側にあるポンプ井から送り出された量である。   The amount of sewage inflow includes the amount of pumped water from the upstream sewage pumping station (the amount of water delivered), which changes due to human operation. It is good also as an inflow amount except pump pumping amount (water supply amount). The upstream pump pumping amount P2 is an amount sent from the pump well on the upstream side in a certain time.

Figure 2006002401
Figure 2006002401

この場合、図2における他の流入量M3が下水流入量として算出される。
もちろん、ポンプ井への流入量を計測する流量計がある場合には、数式1,2を省略し、流量計からの流入量データを用いても良い。この流入量データ算出手段70によりある時点における下水流入量データが算出される。データ蓄積手段60では暦データを主キーとして気象データ、水位データ、揚水量データおよび下水流入量データが関連付けられて登録される。
In this case, the other inflow amount M3 in FIG. 2 is calculated as the sewage inflow amount.
Of course, if there is a flow meter that measures the inflow to the pump well, the equations 1 and 2 may be omitted and the inflow data from the flow meter may be used. The inflow data calculation means 70 calculates sewage inflow data at a certain time. In the data storage means 60, weather data, water level data, pumping amount data, and sewage inflow amount data are associated and registered with calendar data as a main key.

分離手段80は、下水流入量を雨水流入量と汚水流入量とに分離するため、下水流入量データを用いて雨水流入量データと汚水流入量データとを算出する手段である。合流式下水道において必要な処理であり、下水流入量のうち、雨水による雨水流入量と汚水による汚水流入量とを分離する。もちろん雨水流入量と汚水流入量とを合計するともとの下水流入量になる。
これは、晴れている時間の下水流入量は、雨水流入量がないため、汚水流入量にほぼ一致するものと考え、このような晴れている時間の汚水流入量を時間ごと(0時、1時、2時、・・・、23時)に多数日にわたり取得し、さらに同時刻における多数の日の平均を取ることで、時間別の平均汚水流量を算出する。これらを汚水流入量データとする。そして、時間ごとの下水流入量データから対応する時間の汚水流入量データを差し引けば、雨水流量データとなる。晴れの日ではほぼ0となるが、雨の日では雨水流量データが算出される。そして、これら汚水流入量データおよび雨水流入量データをデータ蓄積手段60に登録する。データ蓄積手段60では暦データを主キーとして気象データ、水位データ、揚水量データ、下水流入量データ、汚水流入量データおよび雨水流入量データが関連付けられて登録される。
Separating means 80 is means for calculating rainwater inflow data and sewage inflow data using sewage inflow data in order to separate the sewage inflow into stormwater inflow and sewage inflow. This is a necessary treatment in the combined sewer system, and the stormwater inflow due to rainwater and the sewage inflow due to sewage are separated from the sewage inflow. Of course, the sum of rainwater inflow and sewage inflow yields the original sewage inflow.
This is because the amount of sewage inflow during sunny hours is almost the same as the amount of sewage inflow because there is no inflow of rainwater. Hour, 2 o'clock,..., 23:00) over a number of days, and by taking the average of a number of days at the same time, the average sewage flow rate for each hour is calculated. These are used as sewage inflow data. Then, by subtracting the sewage inflow data for the corresponding time from the sewage inflow data for each hour, it becomes rainwater flow data. Although it is almost 0 on a clear day, rainwater flow rate data is calculated on a rainy day. Then, the sewage inflow data and the rainwater inflow data are registered in the data storage means 60. In the data storage means 60, weather data, water level data, pumped water amount data, sewage inflow data, sewage inflow data and rainwater inflow data are associated and registered with the calendar data as the main key.

モデル構築手段90は、モデルを構築する手段である。図3はモデル構築手段の構造図である。モデル構築手段90は、詳しくは図3で示すように、汚水流入量予測モデル構築手段91、雨水流入量予測モデル構築手段92を備える。なお、雨水流入量予測モデル構築手段91、汚水流入量予測モデル構築手段92で用いられる学習用のデータはデータ蓄積手段60に登録された過去の実績値である。
汚水流入量予測モデル構築手段91は、暦データを入力因子とし、また、汚水流入量データを出力因子とするような汚水流入量予測モデルを構築する手段である。構築された汚水流入量予測モデルはデータ蓄積手段60により蓄積保存される。
雨水流入量予測モデル構築手段92は、気象データを入力因子とし、また、雨水流入量データを出力因子として雨水流入量予測モデルの構築する手段である。構築された雨水流入量予測モデルはデータ蓄積手段60により蓄積保存される。
このような予測モデル構築手段90は、一定期間(例えば一ヶ月ごと)ごとに起動され、これら予測モデルを更新するようにする。なお、雨水流入量予測モデル・汚水流入量予測モデルとしては、RRL法や、ニューラルネットワーク、事例ベース推論、などによる予測するモデルを採用することができるが、最も好ましい予測モデルの詳細については後述される。
The model construction means 90 is a means for constructing a model. FIG. 3 is a structural diagram of the model construction means. As shown in detail in FIG. 3, the model construction unit 90 includes a sewage inflow prediction model construction unit 91 and a stormwater inflow prediction model construction unit 92. Note that the learning data used in the rainwater inflow prediction model construction unit 91 and the sewage inflow prediction model construction unit 92 are past performance values registered in the data storage unit 60.
The sewage inflow prediction model construction means 91 is a means for constructing a sewage inflow prediction model using calendar data as an input factor and sewage inflow data as an output factor. The constructed sewage inflow prediction model is stored and stored by the data storage means 60.
The rainwater inflow prediction model construction means 92 is means for constructing a rainwater inflow prediction model using weather data as an input factor and rainwater inflow data as an output factor. The constructed rainwater inflow prediction model is stored and stored by the data storage means 60.
Such a prediction model construction unit 90 is activated every predetermined period (for example, every month) and updates these prediction models. In addition, as a rainwater inflow prediction model / sewage inflow prediction model, a prediction model using RRL method, neural network, case-based reasoning, or the like can be adopted. Details of the most preferable prediction model will be described later. The

続いて予測手段100では、これら雨水流入量予測モデルおよび汚水流入量予測モデルを用いて予測を行う。図4は予測手段の構造図である。予測手段100は、詳しくは図4で示すように、汚水流入量予測手段101、雨水流入量予測手段102、合算手段103を備えている。
汚水流入量予測手段101は、予測対象時の暦データを汚水流入量予測モデルに入力して汚水流入量予測データを生成する手段である。汚水流入量は時間に依存しており、暦データを入力すれば、汚水流入量予測データを得られる。
雨水流入量予測手段102は、予測対象時の予報である気象データを雨水流入量予測モデルに入力して雨水流入量予測データを生成する手段である。気象データは未来の気象データとなり、予測対象日時の気象データを入力することにより、雨水流入量予測データが得られる。
Subsequently, the prediction unit 100 performs prediction using the rainwater inflow prediction model and the sewage inflow prediction model. FIG. 4 is a structural diagram of the prediction means. As shown in detail in FIG. 4, the prediction unit 100 includes a sewage inflow amount prediction unit 101, a rainwater inflow amount prediction unit 102, and a summation unit 103.
The sewage inflow amount prediction means 101 is a means for generating sewage inflow amount prediction data by inputting the calendar data at the time of prediction into the sewage inflow amount prediction model. The amount of sewage inflow depends on time, and if calendar data is input, sewage inflow prediction data can be obtained.
The stormwater inflow prediction means 102 is means for generating stormwater inflow prediction data by inputting meteorological data, which is a forecast at the time of prediction, into the stormwater inflow prediction model. The meteorological data becomes future meteorological data, and the rainwater inflow prediction data can be obtained by inputting the meteorological data of the forecast date and time.

合算手段102は、生成した汚水流入量予測データと雨水流入量予測データとを合算して下水流入量予測データを生成する手段である。
生成された下水流入量予測データは勿論のこと、雨水流入量予測データと汚水流入量予測データもデータ蓄積手段60により保存される。
このような予測手段100は、一定期間(例えば1時間ごと)ごとに起動され、下水流入量を予測する。
The summing unit 102 is a unit that generates the sewage inflow prediction data by adding the generated sewage inflow prediction data and the rainwater inflow prediction data.
In addition to the generated sewage inflow prediction data, the rainwater inflow prediction data and the sewage inflow prediction data are also stored by the data storage means 60.
Such a prediction unit 100 is activated every certain period (for example, every hour) and predicts the amount of sewage inflow.

このようにして得られた下水流入量予測データは、データ蓄積手段60のデータ発信手段50・ネットワーク2を介して、例えば下水ポンプ場にあるローカル制御用コントローラへ情報を発信することもできる。   The sewage inflow prediction data obtained in this way can be transmitted to the local control controller in the sewage pump station, for example, via the data transmission means 50 of the data storage means 60 and the network 2.

続いて、汚水流入量予測モデルの具体例について説明する。このモデルは特に事例ベース推論によるモデルを採用することが好ましい。出力因子は学習時では汚水流入量データであったが、予測時は汚水流入量予測データであり、相違している。
汚水流入量は、日々の変動が少ないことが知られており、過去の平均をシステムに設定し、それを汚水流入量予測とすることが多い。しかし、人間が設定する方法では、季節が変わるごと、経年変化を起こすごとに再設定しなければならない問題がある。自動化する方法では、先に述べたように特異データを除去すれば適切な平均値を得ることができる。
本発明の事例ベース推論による汚水流入量データ予測モデルでは、以下3点の条件より推論するモデルである。
(1)気象(晴天時)
(2)前後日の汚水流量形状(相関係数で一定値以上)
(3)暦情報(曜日、期間など)
Then, the specific example of a sewage inflow amount prediction model is demonstrated. This model is preferably a model based on case-based reasoning. The output factor was sewage inflow data at the time of learning, but it was sewage inflow prediction data at the time of prediction.
It is known that the amount of sewage inflow varies little every day, and it is often the case that a past average is set in the system and used as a sewage inflow prediction. However, the method set by humans has a problem that it has to be reset every time the season changes and every time the change occurs. In the method of automating, an appropriate average value can be obtained by removing singular data as described above.
The sewage inflow data prediction model based on case-based reasoning of the present invention is a model inferred from the following three conditions.
(1) Weather (when sunny)
(2) Sewage flow shape of the day before and after (more than a certain value in correlation coefficient)
(3) Calendar information (day of week, period, etc.)

(1)気象(晴天時)について
前述のように、下水流入量データは汚水流入量データと雨水流入量データの合算値である。しかしながら、下水への雨水流入量は雨天時にしかない流量であり、晴天時には雨水流入量は0になる。よって、事例ベース推論による予測汚水流入量予測モデルの学習では、過去一定期間内のデータから、晴天時の汚水流入量のみを抽出し、それを汚水流入量データとする。
(2)前後日の汚水流量形状(相関係数で一定値以上)について
ただし、これら抽出条件では特異データを排除することができない。そこで、次に通常データのみ抽出する。特異データとは、主に下水管内清掃によるものであり、一時的に汚水流入量が無くなったり多くなったりする特徴がある。またこのようなデータは毎日行なわれるものではなく、数ヶ月に数回程度の頻度でしかない。よって、特異データは前後日のデータと比較して汚水流入量の形状が大きく異なる特徴がある。逆をいうと通常データは前後日の汚水流入量データが同じ形状をしている。この特徴を利用して、汚水流入量データが前後日と同じ形状のものを抽出する。
ここで、抽出条件は、前後日の形状と比較して相関係数が一定値以上とすることで簡単に抽出できる。
なお、相関係数とは、2つのデータの相関度合いを測る統計指標である。相関係数は−1≦r≦+1の値をとり、r>0のとき正相関,r<0のとき負相関、r=0のときに無相関であることを意味する。
前後日の汚水流入量の形状を判定するときには、例えば、毎時データを用いるときには、対象日の汚水流入量をx1,x2,・・・,x24(mm3/h)とし、前日の汚水流入量をy1,y2,...,y24(mm3/h)とすると、次式にて相関係数が算出される。
(1) About the weather (at the time of fine weather) As mentioned above, the sewage inflow data is the sum of the sewage inflow data and the rainwater inflow data. However, the amount of rainwater inflow into the sewage is a flow rate that can only be obtained in rainy weather, and the amount of rainwater inflow is 0 in fine weather. Therefore, in learning of the predicted sewage inflow prediction model based on case-based reasoning, only the sewage inflow amount during clear weather is extracted from the data within a certain past period and used as sewage inflow amount data.
(2) Concerning the shape of sewage flow before and after the day (more than a certain value in the correlation coefficient) However, these extraction conditions cannot exclude singular data. Therefore, only normal data is extracted next. The peculiar data is mainly due to cleaning in the sewer pipe, and has a characteristic that the amount of inflow of sewage temporarily disappears or increases. In addition, such data is not performed every day, but only once every few months. Therefore, the peculiar data has a feature that the shape of the sewage inflow is greatly different from the data of the previous and next days. In other words, normal data has the same shape as the sewage inflow data on the previous and next days. Using this feature, the sewage inflow data having the same shape as the previous and next days is extracted.
Here, the extraction condition can be easily extracted by setting the correlation coefficient to a certain value or more as compared with the shape of the previous and next days.
The correlation coefficient is a statistical index for measuring the degree of correlation between two data. The correlation coefficient takes a value of −1 ≦ r ≦ + 1, meaning that positive correlation is obtained when r> 0, negative correlation is given when r <0, and no correlation is given when r = 0.
When determining the shape of the amount of sewage inflow on the preceding and following days, for example, when using hourly data, the amount of sewage inflow on the target day is x1, x2, ..., x24 (mm 3 / h), and the amount of sewage inflow on the previous day. Y1, y2,. . . , Y24 (mm 3 / h), the correlation coefficient is calculated by the following equation.

Figure 2006002401
Figure 2006002401

この相関係数rが1に近いとき(例えば0.8以上)に同じ形状と判断し、それ以下のときに、対象日は特異データであると判断する。なお,ここでは前日データのみ比較したが、対象日と対象日の翌日とのデータを比較しても良い。   When the correlation coefficient r is close to 1 (for example, 0.8 or more), it is determined that the shape is the same, and when the correlation coefficient r is less than that, it is determined that the target date is singular data. Although only the previous day data is compared here, data of the target day and the next day of the target day may be compared.

(3)暦情報(曜日、期間など)
さらに、汚水流入量は、曜日により量や立ち上がり時間が若干異なる性質がある。よって、予測対象日と同じ曜日を抽出する。こうして条件に適合した汚水流入量データは、複数日抽出されることになるが、日々のばらつきをなくすため、平均化する。
なお春夏秋冬を分けて平均化すると、さらに予測精度が高まる。
このように曜日や季節(春夏秋冬)により異なる場合もあるため、季節別・曜日別の汚水流入量データとする。
(3) Calendar information (day of week, period, etc.)
Furthermore, the amount of sewage inflow has the property that the amount and rise time differ slightly depending on the day of the week. Therefore, the same day of the week as the prediction target day is extracted. In this way, sewage inflow data that meets the conditions is extracted for a plurality of days, but is averaged to eliminate daily variations.
If the average of spring, summer, autumn and winter is divided, the prediction accuracy will be further increased.
In this way, since it may vary depending on the day of the week and the season (spring, summer, autumn and winter), the sewage inflow data is classified by season and day of the week.

このような汚水流入量予測モデルでは、事例ベース推論にて、特異データを自動的に識別して除去し、通常時の汚水流入量データのみで汚水予測するため、入力される暦データの月日・時間・曜日・季節に対応した汚水流入量データが出力され、従来技術での単なる平均値よりも大幅に予測精度が向上し、上記問題点が解決できる利点がある。   In such a sewage inflow prediction model, singular data is automatically identified and removed by case-based reasoning, and sewage prediction is performed only with normal sewage inflow data. -Sewage inflow data corresponding to time, day of the week, and season is output, and the prediction accuracy is greatly improved over the mere average value in the prior art, and there is an advantage that the above problems can be solved.

以下事例ベース推論による予測汚水流入量予測モデルの学習について記す。
このような汚水流量予測モデルでは、汚水流入量は先に述べたがポンプの清掃のような時以外は、時間別にほぼ同じ値を取るため、特異データを統計的判定方法にて除外した上で過去の時間に対応する実績値である汚水流入量データのデータベース化し、時間データ(たとえば4月10日・水曜日・春・午後10時)を入力すると対応する汚水流入量データを読み出して汚水流入量予測データを生成するモデルとしたため、月日・曜日・季節・時間を入力すると対応する汚水流入量予測データを選択して出力する。
The following describes learning of a predictive sewage inflow prediction model using case-based reasoning.
In such a sewage flow prediction model, the amount of sewage inflow is as described above, but it takes almost the same value by time except during cleaning of the pump. Create a database of sewage inflow data, which is the actual value corresponding to the past time, and input the time data (for example, April 10, Wednesday, Spring, 10:00 pm), read the corresponding sewage inflow data and read the sewage inflow Since it is a model that generates prediction data, when the date, day, season, and time are input, the corresponding sewage inflow prediction data is selected and output.

続いて上記した雨水流入量予測モデルについて説明する。
雨水流入量予測モデルは、具体的にはニューラルネットワークによる予測モデルである。図5は雨水流入量を予測するニューラルネットワークの構成図である。
ニューラルネットワークによるモデルは、例えば、入力層、中間層、出力層の3層構造のパーセプトロンモデルなどであり、入力に対して所望の出力をするように学習できる。
なお、図5のニューラルネットワークは、多入力1出力のニューラルネットワークであるが、もちろん多入力多出力のニューラルネットワークを構築しても良い。
Next, the rainwater inflow prediction model described above will be described.
The rainwater inflow prediction model is specifically a prediction model using a neural network. FIG. 5 is a block diagram of a neural network that predicts the amount of rainwater inflow.
The model based on the neural network is, for example, a perceptron model having a three-layer structure of an input layer, an intermediate layer, and an output layer, and can learn to produce a desired output with respect to the input.
The neural network of FIG. 5 is a multi-input single-output neural network, but a multi-input multi-output neural network may be constructed.

入力因子は気象データのうち特に雨量データを採用し、出力因子は雨水流入量データである。雨量データと雨水流入量データとの間には、通常、非線形な関係があり、線形的予測手法である回帰分析では必ずしも良好な予測結果を得ることができないが、このニューラルネットワークでは、データを学習するだけで適切な予測モデルを構築することができる。   The input factor is the rain data in particular, and the output factor is the rainwater inflow data. There is usually a non-linear relationship between rainfall data and stormwater inflow data, and regression analysis, which is a linear prediction method, does not always give good prediction results, but this neural network learns data. An appropriate prediction model can be constructed simply by doing.

学習段階で用いるニューラルネットワークの入力データは雨量データであり、通常は数時間分のデータを用いる。
なお、これらデータは通常は1時間単位の雨量データを用いるのが一般的であるが、10分単位の雨量データなど短い時間単位の雨量データを用いてもよい、また数時間の累積した雨量データにしてもよい。さらには各雨量データの差分値や絶対量などの変換を行なっても良い。
The input data of the neural network used in the learning stage is rainfall data, and usually several hours of data are used.
In general, these data usually use hourly rainfall data, but it is also possible to use short hourly rainfall data such as 10 minute rainfall data, or several hours of accumulated rainfall data. It may be. Furthermore, a difference value or an absolute amount of each rainfall data may be converted.

出力データとしては、分離手段80により求められてデータ蓄積手段60に登録されたデータであって、数分先から数時間先までの雨水流入量データである。
通常は1時間先までの雨水流入量データを出力するように設計するが、ニューラルネットワークを多出力化し数時間先までの雨水流入量データを一括で出力しても良い。また、雨水流入量データを加工し、現在からの差分値などにしても良い。
The output data is data obtained by the separating means 80 and registered in the data accumulating means 60, and is rainwater inflow data from several minutes to several hours ahead.
Normally, it is designed to output rainwater inflow data up to one hour ahead, but it is also possible to output the rainwater inflow data up to several hours in a lump by outputting multiple neural networks. Further, the rainwater inflow data may be processed to obtain a difference value from the present.

本形態では説明の具体化のため、図5で示すように、時間を異ならせた雨量データ(ある時点の雨量データ値、ある時点より1時間前の雨量データ、ある時点より1時間後の雨量データ)を入力因子とし、また、ある時点より1時間後の雨水流入量データを出力因子としている。   In the present embodiment, as shown in FIG. 5, for example, as shown in FIG. 5, rainfall data (rainfall data values at a certain point in time, rainfall data one hour before a certain point in time, rainfall amounts one hour after a certain point in time) Data) as an input factor, and rainwater inflow data one hour after a certain point in time as an output factor.

続いて上記した雨水流入量予測モデルの学習について説明する。ニューラルネットワークの学習では学習データを作成する必要がある。データ蓄積手段60により過去に蓄積した気象データに含まれる雨量データや、この雨量データに対応する雨水流入量データを時間別に収集し、上記ニューラルネットワークの入出力データセットを多数用意する。例えば、ある時点とある時点の一時間の前後の時刻の3点の雨量データを入力因子とし、また、ある時点から一時間先の雨水流入量データを出力因子とする入出力データセットである。このような入出力データセットを基準となる時点を異ならせて多数準備して、例えば周知のバックプロパゲーション法など周知の各種の学習アルゴリズムを用いてニューラルネットワークを構築する。   Next, learning of the above-described rainwater inflow prediction model will be described. In the learning of the neural network, it is necessary to create learning data. Rain data included in the weather data accumulated in the past by the data storage means 60 and rain water inflow data corresponding to the rain data are collected according to time, and a large number of input / output data sets of the neural network are prepared. For example, it is an input / output data set in which three points of rainfall data at a certain point in time and one hour before and after one hour are used as input factors, and rainwater inflow data one hour after a certain point in time is used as an output factor. A large number of such input / output data sets are prepared at different reference time points, and a neural network is constructed using various known learning algorithms such as a known back propagation method.

続いて上記した雨水流入量予測モデルの予測について説明する。予測段階では、上記にて作成したニューラルネットワークモデルに対し、学習時と同じ入力項目のデータを与えると所望の予測値が出力される。出力因子は学習時では雨水流入量データであったが、予測時は雨水流入量予測データであり、相違している。
例えば、図5で示す雨水流入量予測モデルでは、ある時刻を現在とし、現在の雨量データと1時間前の雨量データは、下水ポンプ場もしくは気象情報サービス会社から配信される実績の雨量データが入力される。1時間先の雨量データは、気象情報サービス会社から配信される予報の雨量データである。
Next, the prediction of the rainwater inflow prediction model described above will be described. In the prediction stage, if the data of the same input items as in the learning are given to the neural network model created above, a desired prediction value is output. The output factor was rainwater inflow data at the time of learning, but it was rainwater inflow prediction data at the time of prediction.
For example, in the rainwater inflow prediction model shown in FIG. 5, the current rainfall data and the rainfall data one hour ago are input as the actual rainfall data distributed from a sewage pump station or a weather information service company. Is done. The rainfall data for one hour ahead is forecast rainfall data distributed from a weather information service company.

なお、図5は1時間先の雨水流入量の説明図であるが、同じニューラルネットワークに対し、入力情報を変えることで1時間先よりさらに先の雨水流入量予測データを得ることができる。例えば、図5の1時間前・ある時点・1時間先という雨量データを、それぞれ現在・1時間後・2時間後の雨量データとすれば、2時間先の雨水流入量予測データを得ることができる。現在の雨量データは、下水ポンプ場もしくは気象情報サービス会社から配信される雨量実績データが入力される。1時間先および2時間先の雨量は、気象情報サービス会社から配信される予報の雨量データである。これにより2時間先の雨水流入量予測データを得ることができる。   FIG. 5 is an explanatory diagram of the amount of rainwater inflow one hour ahead. By changing input information for the same neural network, it is possible to obtain rainwater inflow prediction data further ahead of the hour. For example, if the rainfall data of 1 hour before, a certain time, and 1 hour ahead in FIG. 5 is the current, 1 hour, and 2 hours later, respectively, the rainwater inflow prediction data of 2 hours ahead can be obtained. it can. As current rainfall data, actual rainfall data distributed from a sewage pump station or a weather information service company is input. The 1-hour ahead and 2-hour ahead rainfall is forecast rain data distributed from a weather information service company. Thereby, rainwater inflow prediction data of 2 hours ahead can be obtained.

さらに雨水流入量予測手段102を改良してより短い間隔で予測値を得るようにしても良い。この改良について説明する。
雨水流入量予測手段102は、はn時間(本形態では1時間)後の雨水流量予測データを算出する手段であるが、現在の雨水流量データとn時間先(本形態では1時間先)の雨水流量予測データを、または、n時間(本形態では1時間)隔てて隣接する二つの雨水流量予測データを用いて補間し、時間間隔が短い雨水流量予測データに変換した上で生成する手段としても良い。
Further, the rainwater inflow prediction means 102 may be improved to obtain predicted values at shorter intervals. This improvement will be described.
The rainwater inflow prediction means 102 is a means for calculating rainwater flow rate prediction data after n hours (1 hour in this embodiment), but the current rainwater flow rate data and n hours ahead (1 hour ahead in this embodiment). As means for interpolating rainwater flow rate prediction data, or using two adjacent rainwater flow rate prediction data separated by n hours (in this embodiment, 1 hour), and converting them into rainwater flow rate prediction data having a short time interval. Also good.

電動式ポンプはすぐに起動するがエンジン式ポンプは起動するまでに約5分間要する。ポンプの制御を考慮した場合には5分間隔の下水流量予測に対するニーズがある。しかしながら、気象予報(降雨予報)による雨量データは、1時間単位のデータであるため、原理的に気象予報による雨量データから5分単位の予測はできない。そこで、本発明では、n時間単位の雨水流入量予測データを簡易的にそれよりも短い時間単位のデータに変換する。   The electric pump starts immediately, but the engine pump takes about 5 minutes to start. When pump control is considered, there is a need for sewage flow prediction every 5 minutes. However, since rainfall data based on weather forecasts (rainfall forecasts) are data in units of one hour, in principle, predictions in units of 5 minutes cannot be made from rainfall data based on weather forecasts. Therefore, in the present invention, the rainwater inflow prediction data in units of n hours is simply converted into data in units of time shorter than that.

以下図を参照しつつ具体的に説明する。図6は1時間単位の雨水流入量予測データを5分単位の雨水流入量予測データに補間する補間原理の説明図である。ここでは、ニューラルネットワークから出力される雨水流入量予測データが1時間単位であり、その値を5分単位に変換する例で説明する。
ニューラルネットワークの入出力データは、図6のように1時間単位の雨量データである。入力データを変えることにより1時間〜m時間先までの予測を行なう。例えば、2時間先の雨水流入量予測データを計算する場合には、現在雨量データ、1時間先予報雨量データ値、2時間先予報雨量データを入力する。3時間先の雨水流入量予測データを予測する場合には、1時間先予報雨量データ、2時間先予報雨量データ、3時間先予報雨量データを入力する。こうして1時間単位でm時間先までの雨水流入量予測データが得られる。次に、1時間単位の雨水流入量予測データを5分単位の雨水流入量予測データに線形補完する。図6で示すように、まず中央値である30分の雨水流入量予測データは雨水流入量予測モデルの1時間単位の雨水流入量予測データと一致させる。つぎに、その他の時間には、前後の雨水流入量予測データの値で線形補完する。例えば、次式である。
This will be specifically described below with reference to the drawings. FIG. 6 is an explanatory diagram of an interpolation principle for interpolating rainwater inflow prediction data in units of one hour into rainwater inflow prediction data in units of 5 minutes. Here, an example in which the rainwater inflow prediction data output from the neural network is in units of one hour and the value is converted into units of five minutes will be described.
The input / output data of the neural network is rainfall data in units of one hour as shown in FIG. By changing the input data, prediction is performed from 1 hour to m hours ahead. For example, when calculating rainwater inflow forecast data for two hours ahead, current rain data, one hour ahead forecast rainfall data value, and two hours ahead forecast rain data are input. When predicting rainwater inflow forecast data for 3 hours ahead, 1 hour ahead forecast rainfall data, 2 hours ahead forecast rain data, and 3 hours ahead forecast rain data are input. In this way, rainwater inflow prediction data up to m hours ahead is obtained in units of one hour. Next, the rainwater inflow prediction data in units of one hour is linearly complemented with the rainwater inflow prediction data in units of 5 minutes. As shown in FIG. 6, first, the 30-minute rainwater inflow prediction data, which is the median value, is made to coincide with the rainwater inflow prediction data in an hour unit of the rainwater inflow prediction model. Next, at other times, linear interpolation is performed with the values of preceding and following rainwater inflow prediction data. For example,

(数4)
5分先雨水流入量予測データ
=(5×現在雨水流入量データ+1時間先雨水流入量予測データ)/6
10分先雨水流入量予測データ
=(4×現在雨水流入量データ+2×1時間先雨水流入量予測データ)/6
15分先雨水流入量予測データ
=(3×現在雨水流入量データ+3×1時間先雨水流入量予測データ)/6
20分先雨水流入量予測データ
=(2×現在雨水流入量データ+4×1時間先雨水流入量予測データ)/6
25分先雨水流入量予測データ
=(1×現在雨水流入量データ+5×1時間先雨水流入量予測データ)/6
30分先雨水流入量予測データ
=1時間先雨水流入量予測データ
35分先雨水流入量予測データ
=(11×1時間先雨水流入量予測データ+2時間先雨水流入量予測データ)/12
40分先雨水流入量予測データ
=(10×1時間先雨水流入量予測データ+2×2時間先雨水流入量予測データ)/12
45分先雨水流入量予測データ
=(9×1時間先雨水流入量予測データ+3×2時間先雨水流入量予測データ)/12
50分先雨水流入量予測データ
=(8×1時間先雨水流入量予測データ+4×2時間先雨水流入量予測データ)/12
55分雨水流入量予測データ
=(7×1時間先雨水流入量予測データ+5×2時間先雨水流入量予測データ)/12
60分先雨水流入量予測データ
=(6×1時間先雨水流入量予測データ+6×2時間先雨水流入量予測データ)/12
65分先雨水流入量予測データ
=(7×1時間先雨水流入量予測データ+5×2時間先雨水流入量予測データ)/12
70分先雨水流入量予測データ
=(8×1時間先雨水流入量予測データ+4×2時間先雨水流入量予測データ)/12
75分先雨水流入量予測データ
=(9×1時間先雨水流入量予測データ+3×2時間先雨水流入量予測データ)/12
80分先雨水流入量予測データ
=(10×1時間先雨水流入量予測データ+2×2時間先雨水流入量予測データ)/12
85分先雨水流入量予測データ
=(11×1時間先雨水流入量予測データ+2時間先雨水流入量予測データ)/12
90分先雨水流入量予測データ
=(2時間先雨水流入量予測データ)
(Equation 4)
5 minutes ahead rainwater inflow prediction data = (5 x current rainwater inflow data + 1 hour ahead rainwater inflow prediction data) / 6
10-minute stormwater inflow prediction data = (4 x current stormwater inflow data + 2 x 1-hour stormwater inflow prediction data) / 6
15 minutes ahead stormwater inflow prediction data = (3 x current stormwater inflow data + 3 x 1 hour stormwater inflow prediction data) / 6
20-minute stormwater inflow prediction data = (2 x current stormwater inflow data + 4 x 1 hour stormwater inflow prediction data) / 6
25 minutes ahead stormwater inflow prediction data = (1 x current stormwater inflow data + 5 x 1 hour stormwater inflow prediction data) / 6
30-minute stormwater inflow prediction data = 1-hour stormwater inflow prediction data
35 minutes ahead stormwater inflow forecast data = (11 x 1 hour stormwater inflow forecast data + 2 hours stormwater inflow forecast data) / 12
40 minutes ahead stormwater inflow forecast data = (10 x 1 hour stormwater inflow forecast data + 2 x 2 hours stormwater inflow forecast data) / 12
45 minutes ahead rainwater inflow forecast data = (9 x 1 hour ahead rainwater inflow forecast data + 3 x 2 hours ahead rainwater inflow forecast data) / 12
50 minutes ahead rainwater inflow forecast data = (8 x 1 hour ahead rainwater inflow forecast data + 4 x 2 hours ahead rainwater inflow forecast data) / 12
55-minute rainwater inflow forecast data = (7 x 1 hour ahead rainwater inflow forecast data + 5 x 2 hours ahead rainwater inflow forecast data) / 12
60 minutes ahead stormwater inflow forecast data = (6 x 1 hour stormwater inflow forecast data + 6 x 2 hours stormwater inflow forecast data) / 12
65 minutes ahead rainwater inflow forecast data = (7 x 1 hour ahead rainwater inflow forecast data + 5 x 2 hours ahead rainwater inflow forecast data) / 12
70 minutes ahead rainwater inflow forecast data = (8 x 1 hour ahead rainwater inflow forecast data + 4 x 2 hours ahead rainwater inflow forecast data) / 12
75 minutes ahead rainwater inflow forecast data = (9 x 1 hour ahead rainwater inflow forecast data + 3 x 2 hours ahead rainwater inflow forecast data) / 12
80 minutes ahead rainwater inflow forecast data = (10 x 1 hour ahead rainwater inflow forecast data + 2 x 2 hours ahead rainwater inflow forecast data) / 12
85 minutes ahead stormwater inflow forecast data = (11 x 1 hour stormwater inflow forecast data + 2 hours stormwater inflow forecast data) / 12
90 minutes ahead rainwater inflow forecast data = (2 hours ahead rainwater inflow forecast data)

このように5分単位の予測値を得ることができる。
雨水流入量予測手段102はこのようなものである。
以上本形態の下水流入量予測装置1について説明した。本形態によれば、天候により変動しない汚水流入量と、天候により変動する雨水流入量と、を分けてそれぞれ予測モデルを構築して予想することで、正確な予測が可能となった。
Thus, a predicted value in units of 5 minutes can be obtained.
The rainwater inflow prediction means 102 is such.
The sewage inflow amount prediction device 1 of this embodiment has been described above. According to this embodiment, it is possible to accurately predict by constructing and predicting the prediction model by dividing the inflow amount of sewage that does not vary depending on the weather and the inflow amount of rainwater that varies depending on the weather.

続いて先に説明した予測手段をさらに改良した形態について図を参照しつつ説明する。図7は、他の予測手段100’の構成図である。図4の予測手段100に対し、さらに補正手段104を追加している点が相違している。
補正手段104は、汚水流入量予測データを過去の日時における下水流入量予測データと下水流入量データとから算出した誤差量より補正する手段である。
汚水流入量予測データと雨水流入量予測データとが合算された下水流入量予測データに対し、過去の予測誤差データを用いて補正し、補正下水流入量予測データを生成する。ここに予測誤差データとは過去にある時点の下水流入量を予測した下水流入量予測データと、同じ時点で実際の下水流入量データとの誤差を多数収集し、その平均値を誤差としたものである。補正式は次式のようになる。
Next, a mode in which the prediction means described above is further improved will be described with reference to the drawings. FIG. 7 is a configuration diagram of another prediction unit 100 ′. 4 is different from the prediction unit 100 of FIG. 4 in that a correction unit 104 is further added.
The correction means 104 is means for correcting the sewage inflow prediction data from an error amount calculated from the sewage inflow prediction data and the sewage inflow data at the past date and time.
The sewage inflow prediction data obtained by adding the sewage inflow prediction data and the rainwater inflow prediction data is corrected using the past prediction error data, and corrected sewage inflow prediction data is generated. Here, the prediction error data is a collection of many errors between the sewage inflow prediction data that predicted the sewage inflow at a certain time in the past and the actual sewage inflow data at the same time, and the average value was taken as the error. It is. The correction formula is as follows.

(数5)
補正下水流入量予測データ=下水流入量予測データ+予測誤差データ
(Equation 5)
Corrected sewage inflow prediction data = sewage inflow prediction data + prediction error data

これにより、下水流入量予測データの誤差を補正して下水流入量予測の精度を高めることができる。   Thereby, the error of the sewage inflow amount prediction data can be corrected to improve the accuracy of the sewage inflow amount prediction.

この予測手段100’をさらに改良し、過去の予測誤差データを用いて未来の下水流入量予測データを補正するときに、近い将来の下水流入量予測データの補正量を大きく、遠い将来の下水流入量予測データの補正量を小さくするような手段としても良い。この点について図を参照しつつ説明する。図8は、下水流入量の自己相関係数の特性図を示している。横軸は相関係数を取る二個の下水流入量の時間差である時間間隔を表し、縦軸は自己相関係数を表している。なお,自己相関係数とは、自分自身との相関係数であり、同じデータから時間をずらして2系列のデータを作成し、時間ごとの相関係数を算出した値である。例えば,1時間ごとのデータ{x1,x2,・・・,x48}とすると、1時間先の自己相関係数は、{x1,x2,・・・,x24}と{x2,x3,・・・x25}の2系列のデータを作成して算出する。2時間先の自己相関係数は、{x1,x2,・・・,x24}と{x3,x4,・・・x26}の2系列のデータを作成して算出する。
同様にして、24時間先までの自己相関係数が算出できる。自己相関係数が1に近ければ、自分自身の影響が強く、下水流入量予測の場合では、現在の下水流入量より、未来の流入量が簡単に予測できることを意味し、自己相関係数が0に近ければ、現在の下水流入量より、未来の流入量が予測できないことを意味する(雨量など他の要因が強くなる)。
図8に流量の自己相関係数は、近い将来では高いが、遠い将来になるに従い低下する。つまり、近い将来の下水流入量は現在の流入量に大きく影響されるが、遠い未来の流入量は現在の流入量にあまり影響されないことを示している。
そこで、本発明では、近い将来の予測値は大きく補正し、遠い将来の予測値は小さく補正することにより、無駄な誤差の補正を防止し、予測精度を高めるというものである。
When this prediction means 100 ′ is further improved and the future sewage inflow prediction data is corrected using the past prediction error data, the correction amount of the sewage inflow prediction data in the near future is increased, and the sewage inflow in the far future is increased. Means for reducing the correction amount of the amount prediction data may be used. This point will be described with reference to the drawings. FIG. 8 shows a characteristic diagram of the autocorrelation coefficient of the sewage inflow amount. The horizontal axis represents the time interval that is the time difference between the two sewage inflows taking the correlation coefficient, and the vertical axis represents the autocorrelation coefficient. The autocorrelation coefficient is a correlation coefficient with itself, and is a value obtained by creating two series of data by shifting time from the same data and calculating a correlation coefficient for each time. For example, assuming that the data for each hour {x1, x2,..., X48}, the autocorrelation coefficients one hour ahead are {x1, x2,..., X24} and {x2, x3,. -Create and calculate two series of data of x25}. The autocorrelation coefficient two hours ahead is calculated by creating two series of data {x1, x2,..., X24} and {x3, x4,.
Similarly, the autocorrelation coefficient up to 24 hours ahead can be calculated. If the autocorrelation coefficient is close to 1, it has a strong influence on itself, and in the case of sewage inflow prediction, it means that the future inflow can be predicted more easily than the current sewage inflow. If it is close to 0, it means that the future inflow cannot be predicted from the current inflow of sewage (other factors such as rainfall become stronger).
In FIG. 8, the autocorrelation coefficient of the flow rate is high in the near future, but decreases as the distant future is reached. In other words, the sewage inflow in the near future is greatly influenced by the current inflow, but the inflow in the far future is not significantly affected by the current inflow.
Therefore, in the present invention, the prediction value in the near future is corrected to a large value, and the prediction value in the far future is corrected to a small value, thereby preventing unnecessary error correction and improving the prediction accuracy.

続いて補正方法について説明する。まず過去n時間の平均予測誤差を次式により算出する。   Next, the correction method will be described. First, the average prediction error for the past n hours is calculated by the following equation.

(数6)
平均予測誤差データ=avg(下水流入量予測データ[T−t]−下水流入量データ[T−t])
(Equation 6)
Average prediction error data = avg (sewage inflow prediction data [T−t] −sewage inflow data [T−t])

もしくは次式を用いても良い。
(数7)
平均予測誤差データ=
avg(下水流入量予測データ[T−t]−下水流入量データ[T−t])/下水流入量データ[T−t]
Alternatively, the following equation may be used.
(Equation 7)
Average prediction error data =
avg (Sewage inflow prediction data [T−t] −Sewage inflow data [T−t]) / Sewage inflow data [T−t]

ここで、下水流入量予測データ[T−t]は、時刻t時間前における下水流入量予測データ
下水流入量データ[T−t]は、時刻t時間前における下水流入量データ
Here, the sewage inflow prediction data [Tt] is the sewage inflow prediction data before time t.
Sewage inflow data [T-t] is sewage inflow data before time t

次に平均予測誤差データを次式にて修正する。以下では未来のI時間先までの補正下水流入量予測データを修正する例を示している。数式6を用いた場合は以下のようになる。   Next, the average prediction error data is corrected by the following equation. In the following, an example of correcting the corrected sewage inflow prediction data up to the future I hour ahead is shown. When Equation 6 is used, it is as follows.

(数8)
補正下水流入量予測データi=
下水流入量予測データi+(I−i)×平均予測誤差データ/I
(Equation 8)
Corrected sewage inflow prediction data i =
Sewage inflow prediction data i + (I−i) × average prediction error data / I

iが大きくなる、つまり遠い将来になる程平均予測誤差データが小さくなっている。
また、数式7を使用した場合は次の数式のようになる。
The average prediction error data decreases as i increases, that is, in the far future.
When Expression 7 is used, the following expression is obtained.

(数9)
補正下水流入量予測データi=
下水流入量予測データi/(1+平均予測誤差データ×(I−i)/I)
(Equation 9)
Corrected sewage inflow prediction data i =
Sewage inflow prediction data i / (1 + average prediction error data × (I−i) / I)

ここで、0<i<I、下水流入量予測データiは、i時間先の予測値となる。   Here, 0 <i <I, and the sewage inflow prediction data i is a predicted value i hours ahead.

以上本発明の下水流入量予測装置1について説明した。本発明は上記各手段を含む一台のコンピュータとし、データ蓄積手段60をハードディスクとし、下水流入量データ算出手段70、分離手段80、モデル構築手段90、予測手段100をプログラムによりコンピュータに機能させるような手段としても良い。   The sewage inflow amount prediction device 1 of the present invention has been described above. The present invention is a computer including the above-described means, the data storage means 60 is a hard disk, and the sewage inflow data calculation means 70, the separation means 80, the model construction means 90, and the prediction means 100 are caused to function on a computer by a program. It may be a good means.

続いて本発明の下水流入量予測装置を用いて実際の予測を行った予測結果について図を参照しつつ詳しく説明する。図9は本実施例の下水流入量予測装置による予測結果の説明図である。
なお、今回の実施例1では、シミュレーションのみ実施したため、データ収集手段50によらず、データ蓄積手段60に蓄積されたデータを用いた。
Subsequently, a prediction result obtained by performing actual prediction using the sewage inflow amount prediction device of the present invention will be described in detail with reference to the drawings. FIG. 9 is an explanatory diagram of a prediction result by the sewage inflow amount prediction device of this embodiment.
In the first embodiment, since only simulation was performed, the data stored in the data storage unit 60 was used regardless of the data collection unit 50.

この予測条件は詳しくは以下のようになる
予測対象:合流式下水施設のある下水ポンプ場のポンプ井への下水流入量
予測時点:現時点より5分先〜12時間先(5分間隔)、
学習用データ:過去6ヶ月のデータ
(1)日付(1日単位)
(2)平日、休日(土曜含む)情報(1日単位)
(3)前日データとの相関係数(1日単位)
(4)雨量(5分単位)
(5)下水流入量(5分単位)
The prediction conditions are as follows in detail: Predicted object: Prediction of the amount of sewage flow into the pump well of the sewage pump station where the combined sewage facilities are located: 5 minutes to 12 hours ahead (5 minute intervals),
Data for learning: Data for the past 6 months (1) Date (1 day unit)
(2) Weekdays, holidays (including Saturday) information (1 day unit)
(3) Correlation coefficient with previous day data (1 day unit)
(4) Rainfall (in units of 5 minutes)
(5) Sewage inflow (5 minutes)

手順1.下水流入量データの算出(下水流入量データ算出手段70)
データ蓄積手段60に蓄積された水位データを用いて上記の数式1により下水流入量データの算出を行った。
Procedure 1. Calculation of sewage inflow data (sewage inflow data calculation means 70)
Using the water level data stored in the data storage means 60, the sewage inflow data was calculated according to the above Equation 1.

手順2.汚水・雨水の分離(分離手段80)
続いてこれら下水流入量データから、汚水流入量データと雨水流入量データとに分離生成する。分離のため、事例ベース推論を用い、月ごとに、平日、休日ごとに過去事例を抽出した。抽出条件は以下のとおりである。
(1)指定月を含まない過去4ヶ月
(2)指定曜日(平日or休日)
(3)相関係数0.9以上
(4)1日の雨量0mmの日
(5)下水流入量データはすべて正常にそろっており、欠測がないこと
Procedure 2. Separation of sewage and rainwater (separation means 80)
Subsequently, the sewage inflow data is separated into sewage inflow data and rainwater inflow data. For separation, case-based reasoning was used to extract past cases for each month, weekdays, and holidays. The extraction conditions are as follows.
(1) Past 4 months not including the specified month (2) Specified day of the week (weekdays or holidays)
(3) Correlation coefficient of 0.9 or more (4) Day with a daily rainfall of 0 mm (5) All sewage inflow data are normal and no missing measurements

上記条件では特に1日の雨量0mmの日、つまり雨が降っていないため、下水流入量データを汚水流入量データとみなせるデータを抽出している。このような下水流入量データを抽出すると、各月とも複数日が抽出される。そこで、平日、休日ごとに時間ごとの平均値をとり、それを汚水流入量データとした。今回抽出データは5分単位であるので、汚水流入量データは5分単位で24時間分のデータである。
そしてこのような汚水流入量データを用いて1日の雨量a[mm]の日の雨水流入量データは下水流入量データから上記算出した汚水流入量データを引いた値とした。
こうして、汚水流入量データは月ごとに平日・休日に分けて2種類、雨量別に雨水流入量データは1日ごとに多数作成される。
Under the above conditions, in particular, since the day has a daily rainfall of 0 mm, that is, it is not raining, data that can be regarded as sewage inflow data is extracted. When such sewage inflow data is extracted, a plurality of days are extracted for each month. Therefore, the average value for each hour on weekdays and holidays was taken and used as sewage inflow data. Since the extraction data this time is in units of 5 minutes, the sewage inflow data is data for 24 hours in units of 5 minutes.
Using such sewage inflow data, the daily stormwater inflow data for the daily rainfall a [mm] was obtained by subtracting the calculated sewage inflow data from the sewage inflow data.
In this way, two types of sewage inflow data are created for each month, weekdays and holidays, and a large amount of storm water inflow data is created for each day.

手順3.汚水流入量予測モデル構築(汚水流入量予測モデル構築手段91)
平日の1日24時間ごと汚水流入量データと、休日の24時間ごと汚水流入量データとを登録し、入力として、月、平日、休日、時刻(0時,1時,・・・,23時)を入力すると対応する汚水流入量データを出力する事例ベース推論によるモデルとした。
手順4.雨水流入量予測モデル構築(雨水流入量予測モデル構築手段92)
雨水流入量予測モデルはニューラルネットワークにて構築する。まず最初に学習用データを作成する。雨水流入量データは手順2で算出されているが、上記データの中には晴天時刻のデータ、下水流入量がマイナスのデータ(逆流した場合のデータ)など多数の特異データが含まれている。そこで、以下の条件を満たす雨水流入量データを抜き出した。
(1)雨水流入量0以上、
(2)降雨5mm/h以上
次に、抜き出したデータを以下のフォーマットに従って多数の学習セットを作成する。
1時間前雨量データmm/h、
ある時点の雨量データmm/h、
1時間先雨量データmm/h、
1時間先雨水流入量データmm/h
上記多数のデータセットを用いて、図5のニューラルネットワークをバックプロパゲーション法にて学習する。
Procedure 3. Sewage inflow prediction model construction (sewage inflow prediction model construction means 91)
Register sewage inflow data every 24 hours a day on weekdays and sewage inflow data every 24 hours on holidays, and input as month, weekday, holiday, time (0:00, 1:00, ..., 23:00) ) Is used as a model based on case-based reasoning that outputs corresponding sewage inflow data.
Procedure 4. Rainwater inflow prediction model construction (rainwater inflow prediction model construction means 92)
The rainwater inflow prediction model is constructed by a neural network. First, learning data is created. The rainwater inflow data is calculated in the procedure 2, but the above data includes a lot of specific data such as clear weather time data and negative sewage inflow data (data in the case of reverse flow). Therefore, rainwater inflow data that satisfies the following conditions was extracted.
(1) Rainwater inflow of 0 or more,
(2) Rain 5 mm / h or more Next, a large number of learning sets are created from the extracted data according to the following format.
1 hour before rain data mm / h,
Rain data at a certain point mm / h,
1 hour ahead rainfall data mm / h,
1 hour ahead rainwater inflow data mm 3 / h
The neural network shown in FIG. 5 is learned by the back propagation method using the large number of data sets.

続いて予測手段100’(図7参照)による予測を行う。
手順5.汚水流入量予測(汚水流入量予測手段101)
汚水流入量予測については、月、平日、休日、曜日、季節、時刻(0時,1時,・・・,23時)を入力すると時刻別の汚水流入量予測データを出力する。
手順6.雨水流入量予測(雨水流入量予測手段102)
先に構築したニューラルネットワークモデルの雨水流入量予測モデルを用いて、雨水流入量予測データを生成する。ニューラルネットワークに入力するデータは、1時間先雨水流入量を予測する場合には、以下のデータを入力する。
1時間前雨量データmm/h、
現在雨量データmm/h、
1時間先雨量データmm/h、
1時間先雨水流入量予測データmm/h
そして2時間先予測値を得る場合には、以下のデータを入力する。
現在雨量データmm/h、
1時間先雨量データmm/h、
2時間先雨量データmm/h
2時間先雨水流入量予測データmm/h
同様にして、12時間先雨水流入量予測データまで得る。ここで、未来の雨量については、気象予報値を用いるべきであるが、本実施例ではシミュレーションであるため、実際の雨量データを用いた。
ここで、ニューラルネットワークから出力される予測値は、1時間単位のデータである。そこで、数式4を用いて5分単位の雨水流入量予測データに変換して5分単位12時間先までの雨水流入量予測データを算出した。
Subsequently, prediction is performed by the prediction unit 100 ′ (see FIG. 7).
Procedure 5. Sewage inflow prediction (sewage inflow prediction means 101)
As for sewage inflow prediction, when a month, weekday, holiday, day of the week, season, and time (0:00, 1:00,..., 23:00) are input, sewage inflow prediction data for each time is output.
Procedure 6. Rainwater inflow prediction (rainwater inflow prediction means 102)
The rainwater inflow prediction data is generated using the rainwater inflow prediction model of the neural network model constructed previously. As data to be input to the neural network, the following data is input when the inflow amount of rainwater one hour ahead is predicted.
1 hour before rain data mm / h,
Current rainfall data mm / h,
1 hour ahead rainfall data mm / h,
1 hour ahead rainwater inflow prediction data mm 3 / h
And when obtaining a predicted value for 2 hours ahead, the following data are input.
Current rainfall data mm / h,
1 hour ahead rainfall data mm / h,
2-hour ahead rainfall data mm / h
2-hour ahead rainwater inflow prediction data mm 3 / h
Similarly, 12 hours ahead of stormwater inflow prediction data is obtained. Here, as for the future rainfall, the weather forecast value should be used. However, since this example is a simulation, actual rainfall data was used.
Here, the predicted value output from the neural network is data in units of one hour. Therefore, the rainwater inflow prediction data up to 12 hours ahead by 5 minutes is calculated by converting into rainwater inflow prediction data in 5 minutes using Equation 4.

手順7.汚水・雨水予測値合算(合算手段103)
汚水流入量予測データと雨水流入量予測データを合算した値を下水流入量予測データとした。また、この下水流入量予測データをデータ蓄積手段60に対し、補正前の下水流入量予測データとして保存した。
Step 7. Combined predicted sewage and rainwater (summing means 103)
The sum of sewage inflow prediction data and rainwater inflow prediction data was used as sewage inflow prediction data. Further, the sewage inflow prediction data was stored in the data storage means 60 as sewage inflow prediction data before correction.

手順8.予測の補正(補正手段104)
データ蓄積手段60から読み出した5分前〜1時間前の12点の補正前の下水流入量予測データと同時刻の実績の下水流入量データとの平均予測誤差データを数式6にて算出し、2時間先までの予測値を数式8にて補正した。ここで、2時間先までの予測値は24点あるため、数式8のIの値は24であり、
iの値は5分先で1、10分先で2、・・・2時間先で24となる。
Procedure 8. Prediction correction (correction means 104)
The average prediction error data between the sewage inflow prediction data before correction of 12 points from 5 minutes to 1 hour before read from the data storage means 60 and the actual sewage inflow data at the same time is calculated by Equation 6, The predicted value up to 2 hours ahead was corrected by Equation 8. Here, since there are 24 prediction values up to 2 hours ahead, the value of I in Equation 8 is 24,
The value of i is 1 in 5 minutes, 2 in 10 minutes, and 24 in 2 hours.

上記方法にて予測した値の一例を図9に示す。図9は、12時段階において、12時5分から24時まで予測した例である。合流式下水であるが、実績値と予測値とで際だった差がなく、適切に予測できている。特に降雨前の12時の段階において、未来の降雨情報を用いて流量が増加する傾向を適切に予測できている。   An example of values predicted by the above method is shown in FIG. FIG. 9 shows an example of prediction from 12: 5 to 24:00 at the 12:00 stage. Although it is combined sewage, there is no significant difference between the actual value and the predicted value, and it can be predicted appropriately. In particular, at the twelve o'clock stage before the rain, it is possible to appropriately predict the tendency of the flow rate to increase using the future rain information.

以上、本発明について説明した。本発明を用いることで、分流式下水管、合流式下水管、いずれの場合でも予測可能である。また人口・舗装率など経年変化にも自動的に予測モデルが変化し高精度な予測結果を得ることができる。また、本発明は特殊な方法を用いるものではなく、市販パッケージでも入手可能な事例ベース推論、ニューラルネットワークによる予測であり、容易に下水流入量予測装置を構築できる利点がある。   The present invention has been described above. By using the present invention, it is possible to predict in either case of a diverted sewer pipe or a combined sewer pipe. In addition, the prediction model automatically changes over time, such as population and pavement rate, and high-precision prediction results can be obtained. In addition, the present invention does not use a special method, but is based on case-based reasoning and neural network prediction that can be obtained in a commercial package, and has an advantage that a sewage inflow prediction device can be easily constructed.

本発明を実施するための最良の形態の下水流入量予測装置の構成図である。It is a block diagram of the sewage inflow amount prediction apparatus of the best form for implementing this invention. ポンプ井の説明図である。It is explanatory drawing of a pump well. モデル構築手段の構造図である。It is a structural diagram of a model construction means. 予測手段の構造図である。It is a structural diagram of a prediction means. 雨水流入量を予測するニューラルネットワークの構成図である。It is a block diagram of the neural network which estimates the amount of rainwater inflow. 1時間単位の雨水流入量予測データを5分単位の雨水流入量予測データに補間する補間原理の説明図である。It is explanatory drawing of the interpolation principle which interpolates the rainwater inflow prediction data of a 1-hour unit into the rainwater inflow prediction data of a 5-minute unit. 他の予測手段の構成図である。It is a block diagram of another prediction means. 下水流入量の自己相関係数の特性図である。It is a characteristic view of the autocorrelation coefficient of sewage inflow. 本実施例の下水流入量予測装置による予測結果の説明図である。It is explanatory drawing of the prediction result by the sewage inflow amount prediction apparatus of a present Example. 通常時の汚水流入量−時間特性図である。It is a sewage inflow amount-time characteristic figure at the normal time. 特異時の汚水流入量−時間特性図である。It is a sewage inflow amount-time characteristic figure at the time of peculiarity.

符号の説明Explanation of symbols

1:下水流入量予測装置
10:入力手段
20:出力手段
30:記録媒体読書手段
40:データ発信手段
50:データ収集手段
51:気象データ収集手段
52:暦データ収集手段
53:水位データ収集手段
54:ポンプ揚水量データ収集手段
60:データ蓄積手段
70:下水流入量データ算出手段
80:分離手段
90:モデル構築手段
91:雨水流入量予測モデル構築手段
92:汚水流入量予測モデル構築手段
100,100’:予測手段
101:汚水流入量予測手段
102:雨水流入量予測手段
103:合算手段
104:補正手段
2:ネットワーク
1: Sewage inflow prediction device 10: input means 20: output means 30: recording medium reading means 40: data transmission means 50: data collection means 51: meteorological data collection means 52: calendar data collection means 53: water level data collection means 54 : Pump pumping amount data collection means 60: data storage means 70: sewage inflow amount data calculation means 80: separation means 90: model construction means 91: rainwater inflow amount prediction model construction means 92: sewage inflow amount prediction model construction means 100, 100 ': Predictive means 101: Sewage inflow prediction means 102: Rainwater inflow prediction means 103: Summing means 104: Correction means 2: Network

Claims (6)

下水ポンプ場または終末処理場のポンプ施設へ流入する下水の流入量を予測する下水流入量予測装置において、
気象実績及び気象予報についての気象データを収集する気象データ収集手段と、
暦についての暦データを収集する暦データ収集手段と、
ポンプ井内の下水の水位についての水位データを収集する水位データ収集手段と、
ポンプ井の下水のポンプ揚水量(送水量)についての揚水量データを収集するポンプ揚水量データ収集手段と、
暦データ、気象データ、水位データおよびポンプ揚水量データを蓄積するデータ蓄積手段と、
暦データ、水位データおよび揚水量データに基づいてポンプ井への下水流入量である下水流入量データを算出する下水流入量データ算出手段と、
下水流入量データを分離して汚水流入量データと雨水流入量データとを算出する分離手段と、
暦データを入力因子とし、また、汚水流入量データを出力因子とするような汚水流入量予測モデルを構築する汚水流入量予測モデル構築手段と、
気象データを入力因子とし、また、雨水流入量データを出力因子とするような雨水流入量予測モデルを構築する雨水流入量予測モデル構築手段と、
予測に必要な暦データを汚水流入量予測モデルに入力して汚水流入量予測データを生成する汚水流入量予測手段と、
予測に必要な気象データを雨水流入量予測モデルに入力して雨水流入量予測データを生成する雨水流入量予測手段と、
生成した汚水流入量予測データと雨水流入量予測データとを合算して下水流入量予測データを生成する合算手段と、
を有することを特徴とする下水流入量予測装置。
In the sewage inflow prediction device that predicts the inflow of sewage flowing into the pump facility at the sewage pump station or the final treatment plant,
Meteorological data collection means for collecting meteorological data on weather results and weather forecasts;
Calendar data collection means for collecting calendar data about the calendar;
Water level data collection means for collecting water level data about the sewage water level in the pump well,
Pump pumping data collection means for collecting pumping volume data about pumping volume (pumped volume) of sewage of the pump well;
Data storage means for storing calendar data, meteorological data, water level data and pump pumping volume data;
Sewage inflow data calculation means for calculating sewage inflow data that is sewage inflow to the pump well based on the calendar data, the water level data, and the pumped water data;
Separation means for separating sewage inflow data and calculating sewage inflow data and rainwater inflow data;
A sewage inflow prediction model construction means for constructing a sewage inflow prediction model with calendar data as an input factor and sewage inflow data as an output factor;
A stormwater inflow prediction model construction means for constructing a stormwater inflow prediction model using weather data as an input factor and stormwater inflow data as an output factor;
A sewage inflow forecasting means for generating sewage inflow forecast data by inputting calendar data necessary for prediction into the sewage inflow forecast model,
A stormwater inflow forecasting means for generating stormwater inflow forecast data by inputting meteorological data necessary for the forecast into the stormwater inflow forecast model;
A summing means for generating the sewage inflow prediction data by adding the generated sewage inflow prediction data and the rainwater inflow prediction data;
A sewage inflow prediction apparatus characterized by comprising:
請求項1に記載の下水流入量予測装置において、
前記汚水流入量予測モデルは、
特異データが除去された上で、暦データと汚水流入量データとが連関して登録され、暦データが入力されると対応する汚水流入量データが抽出されることで予測する事例ベース推論による予測モデルであることを特徴とする下水流入量予測装置。
In the sewage inflow prediction device according to claim 1,
The sewage inflow prediction model is
Prediction based on case-based reasoning, in which calendar data and sewage inflow data are linked and registered after singular data is removed, and the corresponding sewage inflow data is extracted when calendar data is input A sewage inflow prediction device characterized by being a model.
請求項1または請求項2に記載の下水流入量予測装置において、
前記雨水流入量予測モデルはニューラルネットワークによる予測モデルであることを特徴とする下水流入量予測装置。
In the sewage inflow prediction device according to claim 1 or 2,
A sewage inflow prediction apparatus, wherein the rainwater inflow prediction model is a prediction model using a neural network.
請求項3に記載の下水流入量予測装置において、
前記雨水流入量予測モデルはn時間単位で雨水流量予測データを算出する予測モデルであり、
前記雨水流入量予測手段は、現在の雨水流量データと雨水流量予測データとを、または、隣接する二つの雨水流量予測データを補間することで時間間隔がn時間よりも短い雨水流量予測データを生成する手段であることを特徴とする下水流入量予測装置。
In the sewage inflow prediction device according to claim 3,
The rainwater inflow prediction model is a prediction model for calculating rainwater flow rate prediction data in units of n hours,
The rainwater inflow prediction means generates rainwater flow prediction data whose time interval is shorter than n hours by interpolating the current rainwater flow data and rainwater flow prediction data or two adjacent rainwater flow prediction data. A sewage inflow prediction device, characterized in that:
請求項1〜請求項4の何れか一項に記載の下水流入量予測装置において、
下水流入量予測データを、過去の日時における下水流入量予測データと下水流入量データとから算出した誤差量より補正する補正手段を備えることを特徴とする下水流入量予測装置。
In the sewage inflow prediction device according to any one of claims 1 to 4,
A sewage inflow prediction apparatus comprising correction means for correcting sewage inflow prediction data from an error amount calculated from sewage inflow prediction data and sewage inflow data at a past date and time.
請求項5に記載の下水流入量予測装置において、
前記補正手段は、過去の予測誤差を用いて未来の予測値を補正するときに、近い将来の予測値の補正量を大きく、遠い将来の予測値の補正量を小さくする手段であることを特徴とする下水流入量予測装置。
In the sewage inflow prediction device according to claim 5,
The correction unit is a unit that increases a correction amount of a near future prediction value and decreases a correction amount of a far future prediction value when correcting a future prediction value using a past prediction error. A sewage inflow prediction device.
JP2004178618A 2004-06-16 2004-06-16 Sewage inflow prediction device Active JP4326408B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2004178618A JP4326408B2 (en) 2004-06-16 2004-06-16 Sewage inflow prediction device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2004178618A JP4326408B2 (en) 2004-06-16 2004-06-16 Sewage inflow prediction device

Publications (2)

Publication Number Publication Date
JP2006002401A true JP2006002401A (en) 2006-01-05
JP4326408B2 JP4326408B2 (en) 2009-09-09

Family

ID=35771024

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2004178618A Active JP4326408B2 (en) 2004-06-16 2004-06-16 Sewage inflow prediction device

Country Status (1)

Country Link
JP (1) JP4326408B2 (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012172320A (en) * 2011-02-18 2012-09-10 Hitachi Ltd Influent discharge water control device, influent discharge water control method and influent discharge water control program
JP2016218599A (en) * 2015-05-18 2016-12-22 三菱電機株式会社 Used power amount management device
JP2019166431A (en) * 2018-03-22 2019-10-03 住友重機械エンバイロメント株式会社 Operation control device and operation control method for wastewater treatment facility
CN112020479A (en) * 2018-04-24 2020-12-01 懿华水处理技术有限责任公司 Sequencing batch reactor system and method

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012172320A (en) * 2011-02-18 2012-09-10 Hitachi Ltd Influent discharge water control device, influent discharge water control method and influent discharge water control program
JP2016218599A (en) * 2015-05-18 2016-12-22 三菱電機株式会社 Used power amount management device
JP2019166431A (en) * 2018-03-22 2019-10-03 住友重機械エンバイロメント株式会社 Operation control device and operation control method for wastewater treatment facility
JP7036490B2 (en) 2018-03-22 2022-03-15 住友重機械エンバイロメント株式会社 Operation control device and operation control method for wastewater treatment equipment
JP2022068358A (en) * 2018-03-22 2022-05-09 住友重機械エンバイロメント株式会社 Operation controller and operation control method for wastewater treatment apparatus
JP7344328B2 (en) 2018-03-22 2023-09-13 住友重機械エンバイロメント株式会社 Operation control device and operation control method for wastewater treatment equipment
CN112020479A (en) * 2018-04-24 2020-12-01 懿华水处理技术有限责任公司 Sequencing batch reactor system and method
US11472726B2 (en) 2018-04-24 2022-10-18 Evoqua Water Technologies Llc Sequencing batch reactor systems and methods
CN112020479B (en) * 2018-04-24 2023-07-25 懿华水处理技术有限责任公司 Sequencing batch reactor system and method

Also Published As

Publication number Publication date
JP4326408B2 (en) 2009-09-09

Similar Documents

Publication Publication Date Title
JP4807565B2 (en) Flow prediction device
KR102073768B1 (en) Drought information supply system based on portal
JP3352153B2 (en) Water distribution flow prediction device
CA3031517C (en) Method and apparatus for model-based control of a water distribution system
US8594936B1 (en) Providing current estimates of precipitation accumulations
Kjeldsen The revitalised FSR/FEH rainfall-runoff method
JP3795775B2 (en) Sewage inflow amount prediction apparatus and method, server apparatus
JP2004019583A (en) Method, device and system for estimating power generation output in wind power generation
CN113221439B (en) BP neural network-based drainage system real-time calibration and dynamic prediction method
EP2312506A1 (en) Computer-based method and device for automatically providing a prediction on a future energy demand to an energy source
JP2006092058A (en) Flow rate estimation device
KR101556463B1 (en) Methods, system, and program of estimating human and property losses by using gust and downpour indices of typhoons
JP4399318B2 (en) Pump control device
KR20160074325A (en) Electricity Demand Index (EDI) Forecasting System with respect to Weather Condition Change
JP4326408B2 (en) Sewage inflow prediction device
CN116740946B (en) Highway traffic flow big data analysis processing method
CN111815043B (en) Flood flow prediction method and system based on storm characteristics
JP4782656B2 (en) Rainwater runoff estimation device and program
KR100430219B1 (en) Prediction method for city water request of water service system
WO2022264422A1 (en) Rainy weather water infiltration rate estimation device, rainy weather water infiltration rate estimation method, and program
JP3268520B2 (en) How to forecast gas demand
Maheswari et al. Ramifications due to urbanization and climate change on stormwater runoff and conduit efficiency: an integrated water management approach
Bonilla et al. Digitalization in small water companies: The Case of Pamplona, Colombia
KR101500795B1 (en) Rainwater inflow facility management supporting method, rainwater inflow facility management supporting device, rainwater inflow facility management supporting system
Bjerregård et al. Probabilistic forecasting of rainfall response in a Danish stormwater tunnel

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20070215

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20080314

RD03 Notification of appointment of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7423

Effective date: 20080314

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20080317

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20080618

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20080618

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20090108

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090121

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090318

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

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

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120619

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 4326408

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

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130619

Year of fee payment: 4

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

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

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250