【手続補正書】
【提出日】平成17年9月5日(2005.9.5)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】０００２
【補正方法】変更
【補正の内容】
【０００２】          【従来の技術】従来の配水量予測装置の一例を図２のブロック図に示す( 特開平8-239868) 。図において、１は河川やダムなどより浄水施設へ水を供給する取水施設、２は河川、ダムなどより取水した水から水質基準に適合した必要量の浄水を安定して供給する浄水施設、３は各家庭、工場などで使用される上水を自然流下、ポンプなどにより需要家に配水する配水施設である。４は需要家に配水された配水流量を測定する流量計、５は温度計、６はオペレータにより実績天候、予想天候、予想気温を入力する天候入力装置である。７は計測値蓄積装置、８は平均変動パターン作成装置、９は雨天日補正係数作成装置、１０はシステム変数データ作成装置、１１は自己回帰モデル作成装置、１４は残差部分の予測値演算装置、１５は配水量予測値演算装置である。つぎに、この配水量予測装置の動作について述べる。 （１）計測値を蓄積する。流量や配水量など計測器によって測定された計測値及びオペレータにより入力された天候データなどを計測値蓄積装置７に蓄積する。 （２）配水量の平均変動パターンを作成する。計測値蓄積装置７に蓄積された配水量を、曜日、時刻、天候によって層別して平均的な変動パターンを平均変動パターン作成装置８により作成する。すなわち、平均変動パターンは、各時刻単位の配水量を曜日毎の直近数日間の平均値から算出する。また、平均変動パターンの更新は、予測日により近い計測値を使用することを目的に毎日行う。しかし、配水量は、雨の影響を受けるため、雨天日の配水量は、平均変動パターンの更新に使用しない。 （３）雨天日の配水量の補正を行なう。過去の配水量から雨天日と雨天日以外の配水量の平均値の比を、各時刻単位に雨天日補正係数作成装置９により計算する。すなわち、雨天日の平均変動パターンは、雨天日以外の平均変動パターンに雨天日補正係数を各時刻毎に乗じて求める。 （４）天候、曜日の影響を除いたデータを作成する。計測値蓄積装置７に記憶している計測値、平均変動パターン作成装置８で作成した平均変動パターン、雨天日補正係数作成装置９で作成した雨天日補正係数をシステム変数データ作成装置１０に入力して、天候、曜日の影響を除いたデータを作成する。配水量から平均変動パターンを引いた部分( 以下、これを残差部分呼ぶ) と最高気温( １日の最高気温) を対象とするデータを作成する。残差部分のデータは、式(1) により求められる。      Ｑst(n) ＝Ｑt(n)−Ｑave(n)・ Ｃ(n)    (n=1,2,3,・・・,24)     ・・・・(1)         ただし、 n    ：時刻、        Ｑst(n) ：時刻n における残差部分の計測値、        Ｑt(n)  ：時刻n における配水量、        Ｑave(n)：時刻n における平均変動パターン、        Ｃ(n)   ：時刻n における雨天日補正係数である。  （５）自己回帰モデルを作成する。システム変数データ作成装置１０で作成したデータを自己回帰モデル作成装置１１に入力して自己回帰モデルを作成する。いま、時刻n におけるプロセスの状態をｋ次元の全変数ベクトルＸ(n) 、時刻n よりm 時点前の全変数ベクトルをＸ(n-m) 、白色ノイズベクトルをＵ(n) 、自己回帰モデルの回帰係数をＡ(m) 、自己回帰モデルの最適次数をM で表すと、その自己回帰表現は、式(2) により表される。  【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００１０
【補正方法】変更
【補正の内容】
【００１０】ただし、N はデータ数、‖dM‖はU(n)の分散共分散行列推定値である。また、MFPEはMultiple Final Prediction Error の頭文字である。このようにして自己回帰係数、白色ノイズの分散および最適モデル次数が求められ、自己回帰モデルが作成される。従って、残差部分の予測を行うために必要な、残差部分と最高気温との関係式を自己回帰モデルから求めることができる。自己回帰モデルの更新は、直近の計測値を使用することを目的に直近数十日の計測値を使用して１日１回行う。 （６）残差部分の予測値を演算する。自己回帰モデル作成装置１１で作成した自己回帰モデルとシステム変数データ作成装置１０で作成したデータから統計的に類推可能な残差部分の24時間予測を残差部分の予測値演算装置１４により演算する。自己回帰モデルを用いた時の24時間予測は、式(6) のように表される。
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】００１２
【補正方法】変更
【補正の内容】
【００１２】ただし、Ｑst(n)p：時刻n における残差部分の予測値、        Ｑst(n) ：時刻n における残差部分の計測値、        Temp(n) ：時刻n の最高気温、        Ａ11(m) ：残差部分の予測値に対する残差部分の自己回帰係数、        Ａ12(m) ：残差部分の予測値に対する最高気温の自己回帰係数、n ：時刻        である。          しかし、１点先以上の予測が必要なため、１点先以上の予測には、残差部分は予測値を使用し、最高気温は予想最高気温を使用する。このようにして得られた残差部分の予測値Ｑst(0)p, Ｑst(1)p・・・・Ｑst(23)pを予測値演算装置１５に出力する。 （７）配水量の予測値を演算する。予測日の天候、曜日、時刻から予測日に該当する平均変動パターンを選択し、近の配水量、最高気温、予想最高気温から求めた残差部分の予測値との和を配水量予測値演算装置１５により演算する。配水量予測式は、式(7) のように決定される。      Ｑt(n)p ＝Ｑst(n)p＋Ｑave(n)・ Ｃ(n)      (n=1,2,3,・・・,24) ・・・・(7)           ただし、Ｑt(n)p は予測配水量である。このように、この方法による配水量予測装置では、自己回帰モデルのモデル次数は１点先の予測を想定した次数を設定し、また、自己回帰モデルのシステム変数に適用する気温も日最高気温を使用しており、浄水場の送水運転で数年間問題なく稼働している。
【手続補正４】
【補正対象書類名】明細書
【補正対象項目名】００１５
【補正方法】変更
【補正の内容】
【００１５】        【発明の実施の形態】以下、本発明の実施例を図を参照して詳細に説明する。図１は、本発明の配水量予測装置を示すブロック図である。図において、１２は24点予測時の最適モデル次数を演算するモデル次数演算装置、１３は各時刻の気温を気象庁から発表された予想最高気温、予想最低気温を使用して演算する予想気温演算装置である。その他の符号の説明は、従来技術と同様であるため省略する。つぎに、本配水量予測装置の動作について述べる。 ［１］計測値を計測値蓄積装置７に蓄積する（従来技術と同じ）。 ［２］配水量の平均変動パターンを平均変動パターン作成装置８により作成する（従来技術と同じ）。 ［３］雨天日の配水量の補正を雨天日補正係数作成装置９により行なう（従来技術と同じ）。 ［４］天候、曜日の影響を除いたデータをシステム変数データ作成装置10二より作成する（従来技術と同じ）。 ［５］自己回帰モデルを自己回帰モデル作成装置１１により作成する（従来技術と同じ）。 ［６］自己回帰モデルの最適モデル次数を演算する。        24点予測時の最適モデル次数をモデル次数演算装置１２により演算する。自己回帰モデル作成装置１１で作成した自己回帰モデル係数とシステム変数データ作成装置１０で作成したデータから統計的に類推可能な残差部分の24時間予測を演算する。自己回帰モデルを用いた時の24時間予測は、前述の(6) 式で表される。ただし、１点先以上の予測が必要なため、１点先以上の予測には、残差部分は予測値を使用し、気温は予測した気温を使用する。最適モデル次数を求めるために先ず残差部分の予測誤差を(8) 式で求める。【手続補正５】
【補正対象書類名】明細書
【補正対象項目名】００１７
【補正方法】変更
【補正の内容】
【００１７】ただし、M はモデル次数、Y(M)は残差部分の予測誤差である。        (8) 式においてモデル次数Ｍを１からデータ数の３分の１程度まで変動させ、最も予測誤差Y(M)が小さくなるモデル次数Ｍを求める。このように長期予測時における最適モデル次数を求めることができる。この最適モデル次数を残差部分の予測値演算装置１４に入力する。 ［７］各時刻の予想気温を演算する。予想気温は、各時刻の気温を気象庁から発表された予想最高気温、予想最低気温を用いて、予想気温演算装置１３により演算する。最高気温の係数を１、最低気温の係数を０として過去の気温データから各時刻の気温係数を(9) 式にて算出する。
【手続補正６】
【補正対象書類名】明細書
【補正対象項目名】００２１
【補正方法】変更
【補正の内容】
【００２１】ただし、TempMAXp：予想最高気温、        TempMINp：予想最低気温、        Temp(n)p：時刻n における予想気温である。        このように各時刻の予想気温を演算することができる。この予想気温を自己回帰モデルのシステム変数として残差部分の予測値演算装置１４に入力する。 ［８］残差部分の予測値を残差部分の予測値演算装置１４により演算する（従来技術と同じ）。ただし、気温は予想気温を使用する。 ［９］配水量の予測値を配水量予測値演算装置１５により演算する（従来技術と同じ）。        本実施例では、２４時間の予測例について述べたが、本発明を用いることにより３日先程度の長期配水量予測についても精度良く行うことができる。
【手続補正７】
【補正対象書類名】明細書
【補正対象項目名】図面の簡単な説明
【補正方法】変更
【補正の内容】
【図面の簡単な説明】                【図１】本発明の配水量予測装置を示すブロック図である。        【図２】従来の配水量予測装置を示すブロックである。        【符号の説明】        １  取水施設        ２  浄水施設        ３  配水施設        ４  流量計        ５  気温計        ６  天候入力装置        ７  計測値蓄積装置        ８  平均変動パターン作成装置        ９  雨天日補正係数作成装置        １０  システム変数データ作成装置        １１  自己回帰モデル作成装置        １２  モデル次数演算装置        １３  予想気温演算装置        １４  残差部分の予測値演算装置        １５  配水量予測値演算装置