JP4453791B2 - 配水量予測装置 - Google Patents
配水量予測装置 Download PDFInfo
- Publication number
- JP4453791B2 JP4453791B2 JP2000239577A JP2000239577A JP4453791B2 JP 4453791 B2 JP4453791 B2 JP 4453791B2 JP 2000239577 A JP2000239577 A JP 2000239577A JP 2000239577 A JP2000239577 A JP 2000239577A JP 4453791 B2 JP4453791 B2 JP 4453791B2
- Authority
- JP
- Japan
- Prior art keywords
- prediction
- value
- water distribution
- distribution amount
- temperature
- 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.)
- Expired - Lifetime
Links
Landscapes
- Feedback Control In General (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Description
【発明の属する技術分野】
本発明は、取水施設、浄水施設、配水設備の運用計画を立てるために必要な配水量予測を、精度良く行うことが出来る配水量予測装置に関するものである。
【0002】
【従来の技術】
従来の配水量予測装置の一例を図2のブロック図に示す( 特開平8-239868) 。図において、1は河川やダムなどより浄水施設へ水を供給する取水施設、2は河川、ダムなどより取水した水から水質基準に適合した必要量の浄水を安定して供給する浄水施設、3は各家庭、工場などで使用される上水を自然流下、ポンプなどにより需要家に配水する配水施設である。4は需要家に配水された配水流量を測定する流量計、5は温度計、6はオペレータにより実績天候、予想天候、予想気温を入力する天候入力装置である。7は計測値蓄積装置、8は平均変動パターン作成装置、9は雨天日補正係数作成装置、10はシステム変数データ作成装置、11は自己回帰モデル作成装置、14は残差部分の予測値演算装置、15は配水量予測値演算装置である。つぎに、この配水量予測装置の動作について述べる。
(1)計測値を蓄積する。流量や配水量など計測器によって測定された計測値及びオペレータにより入力された天候データなどを計測値蓄積装置7に蓄積する。
(2)配水量の平均変動パターンを作成する。計測値蓄積装置7に蓄積された配水量を、曜日、時刻、天候によって層別して平均的な変動パターンを平均変動パターン作成装置8により作成する。すなわち、平均変動パターンは、各時刻単位の配水量を曜日毎の直近数日間の平均値から算出する。また、平均変動パターンの更新は、予測日により近い計測値を使用することを目的に毎日行う。しかし、配水量は、雨の影響を受けるため、雨天日の配水量は、平均変動パターンの更新に使用しない。
(3)雨天日の配水量の補正を行なう。過去の配水量から雨天日と雨天日以外の配水量の平均値の比を、各時刻単位に雨天日補正係数作成装置9により計算する。すなわち、雨天日の平均変動パターンは、雨天日以外の平均変動パターンに雨天日補正係数を各時刻毎に乗じて求める。
(4)天候、曜日の影響を除いたデータを作成する。計測値蓄積装置7に記憶している計測値、平均変動パターン作成装置8で作成した平均変動パターン、雨天日補正係数作成装置9で作成した雨天日補正係数をシステム変数データ作成装置10に入力して、天候、曜日の影響を除いたデータを作成する。配水量から平均変動パターンを引いた部分( 以下、これを残差部分呼ぶ) と最高気温( 1日の最高気温) を対象とするデータを作成する。残差部分のデータは、式(1) により求められる。
Qst(n) =Qt(n)−Qave(n)・ C(n) (n=1,2,3,・・・,24) ・・・・(1)
ただし、 n :時刻、
Qst(n) :時刻n における残差部分の計測値、
Qt(n) :時刻n における配水量、
Qave(n):時刻n における平均変動パターン、
C(n) :時刻n における雨天日補正係数である。
(5)自己回帰モデルを作成する。システム変数データ作成装置10で作成したデータを自己回帰モデル作成装置11に入力して自己回帰モデルを作成する。いま、時刻n におけるプロセスの状態をk次元の全変数ベクトルX(n) 、時刻n よりm 時点前の全変数ベクトルをX(n-m) 、白色ノイズベクトルをU(n) 、自己回帰モデルの回帰係数をA(m) 、自己回帰モデルの最適次数をM で表すと、その自己回帰表現は、式(2) により表される。
【0003】
【数1】
【0004】
従って自己回帰モデルの作成とは、自己回帰係数、白色ノイズベクトルの分散および自己回帰モデルの最適次数決定に帰結される。自己回帰係数A(m)は、要素をAij(m)とし、次の連立方程式をi=1,2,3,・・・・,kについて解くことにより求められる。ただし、Xi、Xjの相互分散をRij(l)、自己回帰係数の要素をAij(m)とすると式(3) のようになる。
【0005】
【数2】
【0006】
白色ノイズベクトルU(n) の要素をεi(n)とすると、その残差分散値σi2は、式(4) で表される。
【0007】
【数3】
【0008】
なお、モデルの最適次数Mは、予測誤差を表す式(5) のMFPE(M) を最小にする値である。
【0009】
【数4】
【0010】
ただし、N はデータ数、‖dM‖はU(n)の分散共分散行列推定値である。また、MFPEはMultiple Final Prediction Error の頭文字である。このようにして自己回帰係数、白色ノイズの分散および最適モデル次数が求められ、自己回帰モデルが作成される。従って、残差部分の予測を行うために必要な、残差部分と最高気温との関係式を自己回帰モデルから求めることができる。自己回帰モデルの更新は、直近の計測値を使用することを目的に直近数十日の計測値を使用して1日1回行う。
(6)残差部分の予測値を演算する。自己回帰モデル作成装置11で作成した自己回帰モデルとシステム変数データ作成装置10で作成したデータから統計的に類推可能な残差部分の24時間予測を残差部分の予測値演算装置14により演算する。自己回帰モデルを用いた時の24時間予測は、式(6) のように表される。
【0011】
【数5】
【0012】
ただし、Qst(n)p:時刻n における残差部分の予測値、
Qst(n) :時刻n における残差部分の計測値、
Temp(n) :時刻n の最高気温、
A11(m) :残差部分の予測値に対する残差部分の自己回帰係数、
A12(m) :残差部分の予測値に対する最高気温の自己回帰係数、n :時刻
である。
しかし、1点先以上の予測が必要なため、1点先以上の予測には、残差部分は予測値を使用し、最高気温は予想最高気温を使用する。このようにして得られた残差部分の予測値Qst(0)p, Qst(1)p・・・・Qst(23)pを予測値演算装置15に出力する。
(7)配水量の予測値を演算する。予測日の天候、曜日、時刻から予測日に該当する平均変動パターンを選択し、近の配水量、最高気温、予想最高気温から求めた残差部分の予測値との和を配水量予測値演算装置15により演算する。配水量予測式は、式(7) のように決定される。
Qt(n)p =Qst(n)p+Qave(n)・ C(n) (n=1,2,3,・・・,24) ・・・・(7)
ただし、Qt(n)p は予測配水量である。このように、この方法による配水量予測装置では、自己回帰モデルのモデル次数は1点先の予測を想定した次数を設定し、また、自己回帰モデルのシステム変数に適用する気温も日最高気温を使用しており、浄水場の送水運転で数年間問題なく稼働している。
【0013】
【発明が解決しようとする課題】
ところが、このように従来の配水量予測装置では、自己回帰モデルのモデル次数は1点先の予測を想定した次数を設定し、自己回帰モデルのシステム変数に適用する気温も日最高気温を使用していたため、数日先の長期予測には予測精度がなお不十分であった。
そこで、本発明は数日先の配水量の長期予測を高精度に行なうことのできる配水量予測装置を提供することを目的とする。
【0014】
【課題を解決するための手段】
上記課題を解決するため、浄水場の運用計画を立てるために必要な配水量を予測する配水量予測装置において、
各時刻の配水量の計測値、天候、気温、曜日を含む配水量に影響する条件の計測値を蓄積する計測値蓄積装置と、
配水量を曜日、時刻、天候によって層別し、各層別された雨天日以外の日の配水量の曜日毎の直近数日間について時刻毎に平均値を求める平均変動パターン作成装置と、
各時刻について雨天日と雨天日以外の日との配水量の比を求める補正係数作成装置と、
前記配水量と前記平均変動パターンの差である残差部分と一日の最高気温をシステム変数とするデータを作成するシステム変数データ作成装置と、
前記システム変数データ作成装置の変数を使用して自己回帰モデルを作成する自己回帰モデル作成装置と、
前記自己回帰モデルと、前記システム変数データ作成装置により作成されたデータとから残差部分の予測値を演算する残差部分の予測値演算装置と、
前記平均変動パターンと、前記残差部分の予測値との和を配水量予測値とする配水量予測値演算装置とを備え、
前記自己回帰モデル作成装置と前記残差部分の予測値演算装置との間に配置され、前記自己回帰モデルおよび前記システム変数データ作成装置の変数から残差部分の24時間予測を演算し、その予測における誤差が最小となる最適モデル次数を求めるモデル次数演算装置と、
前記計測値蓄積装置の気温データに基づいて各時刻の予想気温を演算し、前記残差部分の予測値演算装置にシステム変数として与える予想気温演算装置とを備えたものである。
【0015】
【発明の実施の形態】
以下、本発明の実施例を図を参照して詳細に説明する。図1は、本発明の配水量予測装置を示すブロック図である。図において、12は24点予測時の最適モデル次数を演算するモデル次数演算装置、13は各時刻の気温を気象庁から発表された予想最高気温、予想最低気温を使用して演算する予想気温演算装置である。その他の符号の説明は、従来技術と同様であるため省略する。つぎに、本配水量予測装置の動作について述べる。
[1]計測値を計測値蓄積装置7に蓄積する(従来技術と同じ)。
[2]配水量の平均変動パターンを平均変動パターン作成装置8により作成する(従来技術と同じ)。
[3]雨天日の配水量の補正を雨天日補正係数作成装置9により行なう(従来技術と同じ)。
[4]天候、曜日の影響を除いたデータをシステム変数データ作成装置10二より作成する(従来技術と同じ)。
[5]自己回帰モデルを自己回帰モデル作成装置11により作成する(従来技術と同じ)。
[6]自己回帰モデルの最適モデル次数を演算する。
24点予測時の最適モデル次数をモデル次数演算装置12により演算する。自己回帰モデル作成装置11で作成した自己回帰モデル係数とシステム変数データ作成装置10で作成したデータから統計的に類推可能な残差部分の24時間予測を演算する。自己回帰モデルを用いた時の24時間予測は、前述の(6) 式で表される。ただし、1点先以上の予測が必要なため、1点先以上の予測には、残差部分は予測値を使用し、気温は予測した気温を使用する。最適モデル次数を求めるために先ず残差部分の予測誤差を(8) 式で求める。
【0016】
【数6】
【0017】
ただし、M はモデル次数、Y(M)は残差部分の予測誤差である。
(8) 式においてモデル次数Mを1からデータ数の3分の1程度まで変動させ、最も予測誤差Y(M)が小さくなるモデル次数Mを求める。このように長期予測時における最適モデル次数を求めることができる。この最適モデル次数を残差部分の予測値演算装置14に入力する。
[7]各時刻の予想気温を演算する。予想気温は、各時刻の気温を気象庁から発表された予想最高気温、予想最低気温を用いて、予想気温演算装置13により演算する。最高気温の係数を1、最低気温の係数を0として過去の気温データから各時刻の気温係数を(9) 式にて算出する。
【0018】
【数7】
【0019】
ただし、TempConf(n) :時刻n の気温係数、
TempMAX :直近数ヶ月の日平均最高気温、
TempMIN :直近数ヶ月の日平均最低気温、
Temp(n) :直近数ヶ月の時刻n の平均気温である。
このようにして求められた各時刻の気温係数と予想最高気温、予想最低気温から各時刻の予想気温を(10)式によって求める。
【0020】
【数8】
【0021】
ただし、TempMAXp:予想最高気温、
TempMINp:予想最低気温、
Temp(n)p:時刻n における予想気温である。
このように各時刻の予想気温を演算することができる。この予想気温を自己回帰モデルのシステム変数として残差部分の予測値演算装置14に入力する。
[8]残差部分の予測値を残差部分の予測値演算装置14により演算する(従来技術と同じ)。ただし、気温は予想気温を使用する。
[9]配水量の予測値を配水量予測値演算装置15により演算する(従来技術と同じ)。
本実施例では、24時間の予測例について述べたが、本発明を用いることにより3日先程度の長期配水量予測についても精度良く行うことができる。
【0022】
【発明の効果】
以上述べたように、本発明によれば最適モデル次数を求めるモデル次数演算装置と各時刻の予想気温を自動的に演算し気温の影響を考慮する予想気温演算装置とをもうけたので、数日先の長期予測に対して精度の高い配水量予測装置を得る効果がある。さらに、本発明の配水量予測装置により、最適な送水計画を立案できる効果もある。
【図面の簡単な説明】
【図1】本発明の配水量予測装置を示すブロック図である。
【図2】従来の配水量予測装置を示すブロックである。
【符号の説明】
1 取水施設
2 浄水施設
3 配水施設
4 流量計
5 気温計
6 天候入力装置
7 計測値蓄積装置
8 平均変動パターン作成装置
9 雨天日補正係数作成装置
10 システム変数データ作成装置
11 自己回帰モデル作成装置
12 モデル次数演算装置
13 予想気温演算装置
14 残差部分の予測値演算装置
15 配水量予測値演算装置
Claims (1)
- 各時刻の配水量の計測値、天候、気温、曜日を含む配水量に影響する条件の計測値を蓄積する計測値蓄積装置と、
配水量を曜日、時刻、天候によって層別し、各層別された雨天日以外の日の配水量の曜日毎の直近数日間について時刻毎に平均値を求める平均変動パターン作成装置と、
各時刻について雨天日と雨天日以外の日との配水量の比を求める補正係数作成装置と、
前記配水量と前記平均変動パターンの差である残差部分と一日の最高気温をシステム変数とするデータを作成するシステム変数データ作成装置と、
前記システム変数データ作成装置の変数を使用して自己回帰モデルを作成する自己回帰モデル作成装置と、
前記自己回帰モデルと、前記システム変数データ作成装置により作成されたデータとから残差部分の予測値を演算する残差部分の予測値演算装置と、
前記平均変動パターンと、前記残差部分の予測値との和を配水量予測値とする配水量予測値演算装置とを備え、浄水場の運用計画を立てるために必要な配水量を予測する配水量予測装置において、
前記自己回帰モデル作成装置と前記残差部分の予測値演算装置との間に配置され、前記自己回帰モデルおよび前記システム変数データ作成装置の変数から残差部分の24時間予測を演算し、その予測における誤差が最小となる最適モデル次数を求めるモデル次数演算装置と、
前記計測値蓄積装置の気温データに基づいて各時刻の予想気温を演算し、前記残差部分の予測値演算装置にシステム変数として与える予想気温演算装置とを備えたことを特徴とする配水量予測装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000239577A JP4453791B2 (ja) | 2000-08-08 | 2000-08-08 | 配水量予測装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000239577A JP4453791B2 (ja) | 2000-08-08 | 2000-08-08 | 配水量予測装置 |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2002054182A JP2002054182A (ja) | 2002-02-20 |
JP2002054182A5 JP2002054182A5 (ja) | 2005-11-10 |
JP4453791B2 true JP4453791B2 (ja) | 2010-04-21 |
Family
ID=18731091
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2000239577A Expired - Lifetime JP4453791B2 (ja) | 2000-08-08 | 2000-08-08 | 配水量予測装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4453791B2 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103485386A (zh) * | 2013-09-10 | 2014-01-01 | 温州大学 | 一种基于灰色关联法的变频恒压供水系统控制方法 |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103487187B (zh) * | 2013-09-10 | 2015-05-20 | 温州大学 | 一种基于反求法的变频供水系统运行效率在线检测方法 |
CN103488082B (zh) * | 2013-09-10 | 2015-12-09 | 温州大学 | 一种基于反求法的高效变频恒压供水系统控制方法 |
CN103487186B (zh) * | 2013-09-10 | 2015-04-15 | 温州大学 | 基于灰色关联法的变频供水系统运行效率在线检测方法 |
CN103556677B (zh) * | 2013-09-10 | 2015-07-22 | 台州神能电器有限公司 | 一种高效变频恒压供水系统的控制方法 |
JP5886339B2 (ja) * | 2014-02-28 | 2016-03-16 | 中国電力株式会社 | 気温予測システム、気温予測方法 |
JP6467953B2 (ja) * | 2015-01-30 | 2019-02-13 | 中国電力株式会社 | 気温予測システム、気温予測方法およびプログラム |
-
2000
- 2000-08-08 JP JP2000239577A patent/JP4453791B2/ja not_active Expired - Lifetime
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103485386A (zh) * | 2013-09-10 | 2014-01-01 | 温州大学 | 一种基于灰色关联法的变频恒压供水系统控制方法 |
CN103485386B (zh) * | 2013-09-10 | 2014-12-24 | 温州大学 | 一种基于灰色关联法的变频恒压供水系统控制方法 |
Also Published As
Publication number | Publication date |
---|---|
JP2002054182A (ja) | 2002-02-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gustafsson | Control-theoretic techniques for stepsize selection in implicit Runge-Kutta methods | |
Trémolet | Incremental 4D-Var convergence study | |
JP4453791B2 (ja) | 配水量予測装置 | |
JP3336135B2 (ja) | プロセスに指令を与えるための方法および装置 | |
JPH059809B2 (ja) | ||
CN113256005A (zh) | 一种基于神经网络模型的电站水位过程预测方法及设备 | |
CN108830415A (zh) | 一种基于非线性关联模型的气象敏感负荷功率估算方法 | |
JP2008154430A (ja) | 供給計画策定システム、供給計画策定プログラム及び供給計画策定方法 | |
JPH0995917A (ja) | 河川水位予測装置 | |
Veersé et al. | Limited-memory BFGS diagonal preconditioners for a data assimilation problem in meteorology | |
JP4407420B2 (ja) | 配水量予測装置 | |
JP2006333687A (ja) | 電力負荷曲線生成方法及び熱負荷曲線生成方法 | |
JPH08239868A (ja) | 配水量予測方法と装置 | |
JP3626851B2 (ja) | 河川水位予測装置 | |
CN109059875A (zh) | 一种驱动完美模型开展月尺度径流预报的方法 | |
JP4745759B2 (ja) | 河川汚濁負荷源推定装置、方法、およびプログラム | |
JP5242084B2 (ja) | 受水流量算出方法および装置 | |
Alexe et al. | On the discrete adjoints of adaptive time stepping algorithms | |
KR101622068B1 (ko) | 설정 계산 시스템의 학습 장치 및 학습 방법 | |
CN117824772B (zh) | 一种天然气流量计量自适应补偿方法、系统、终端及介质 | |
JP2712688B2 (ja) | 配水量の需要予測装置 | |
JP2001140287A (ja) | 配水量予測システム及び配水量予測方法 | |
CN113158394B (zh) | 基于蒸发误差补偿的流域水文模型误差修正方法和装置 | |
JP2023174396A (ja) | ろ過池洗浄スケジューリング装置、ろ過池洗浄スケジューリングシステム及びろ過池洗浄スケジューリング方法 | |
CN114757036A (zh) | 一种基于逐次逼近法的河道断面自动插补方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050905 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050905 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050905 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20080117 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080515 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080703 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20081222 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090109 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20090706 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090925 |
|
A911 | Transfer to examiner for re-examination before appeal (zenchi) |
Free format text: JAPANESE INTERMEDIATE CODE: A911 Effective date: 20091013 |
|
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: 20100113 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130212 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4453791 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20100126 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140212 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20150212 Year of fee payment: 5 |
|
EXPY | Cancellation because of completion of term |