(57)【要約】【課題】アンサンブル予報を基に適切に需要予測を行う。【解決手段】ユーザからの指定に従って、配信期間内を配信日としており、需要予測対象日時と指定配信日の所定配信時刻との差である予報間隔を予報時刻とする、指定地点における複数シナリオ分の気象予報値を抽出し、複数シナリオ分の気象予報値の中央値又は平均値によって需要予測値分布の高低を調整する項と、気象予報値と当該気象予報値の中央値又は平均値との差によって需要予測値分布の幅を調整する項とを有する需要予測モデルの係数を、算出される需要予測値の分布がユーザ設定に従うように算出し、需要予測モデルの係数を設定した需要予測モデルに対して、指定配信日における予報間隔を予報時刻とする、指定地点における複数シナリオ分の気象予報値を適用して複数シナリオ分の需要予測値を算出する。【選択図】図１
2009225550.tif 000002
【特許請求の範囲】
【請求項１】
電力需要実績値を単位時間毎に格納する電力需要実績データ格納部と、１又は複数の地点における複数予報時刻の気象予報値を複数シナリオ分配信日毎に格納する気象予報データ格納部とにアクセス可能なコンピュータにより実行される需要予測処理方法であって、
ユーザから指定された処理対象配信期間と指定地点と需要予測対象日時と需要予測のための指定配信日とに従って、前記気象予報データ格納部から、前記処理対象配信期間内を配信日としており、前記需要予測対象日時と前記指定配信日の所定配信時刻との差である予報間隔を予報時刻とする、前記指定地点における前記複数シナリオ分の気象予報値を抽出し、抽出気象予報値格納部に格納するステップと、
前記電力需要実績データ格納部から、前記処理対象配信期間内の各日の前記所定配信時刻と前記予報間隔とを加算することによって得られる仮想需要予測日時における電力需要実績値を抽出し、抽出電力需要実績値格納部に格納するステップと、
前記複数シナリオ分の気象予報値の中央値又は平均値によって需要予測値分布の高低を調整する項と、前記気象予報値と前記気象予報値の中央値又は平均値との差によって前記需要予測値分布の幅を調整する項とを有する需要予測モデルの係数を、前記抽出気象予報値格納部に格納されている前記気象予報値と前記需要予測モデルとによって算出される需要予測値の分布が、前記抽出電力需要実績値格納部に格納されている前記電力需要実績値に基づくユーザ設定に従うように算出し、需要予測モデルデータ格納部に格納する需要予測モデル生成ステップと、
前記需要予測モデルデータ格納部に格納されている前記需要予測モデルの係数を設定した前記需要予測モデルに対して、前記気象予報データ格納部に格納されており且つ前記指定配信日における前記予報間隔を予報時刻とする、前記指定地点における前記複数シナリオ分の気象予報値を適用して、前記複数セットの需要予測値を算出し、需要予測シナリオデータ格納部に格納するステップと、
を含む需要予測処理方法。
【請求項２】
前記需要予測値分布の幅を調整する項が、
前記需要予測値分布において前記中央値又は平均値以上の部分の分布の幅を調整する項と、
前記需要予測値分布において前記中央値又は平均値より下の部分の分布の幅を調整する項と
を含む請求項１記載の需要予測処理方法。
【請求項３】
前記需要予測モデルが、
前記気象予報値と前記気象予報値の中央値又は平均値との差の二乗についての項
をさらに含む請求項１記載の需要予測処理方法。
【請求項４】
前記ユーザ設定が、
前記電力需要実績値の最大値と最小値をカバーする可能な限り幅の狭い分布
である請求項１記載の需要予測処理方法。
【請求項５】
前記需要予測モデル生成ステップにおいて、ＬＰ法を用いて前記需要予測モデルの係数を算出する
請求項４記載の需要予測処理方法。
【請求項６】
前記ユーザ設定が、
前記抽出気象予報値格納部に格納されている前記気象予報値と前記需要予測モデルとによって算出される需要予測値の昇順系列又は降順系列を順に区間閾値とする各区間又は特定区間における前記電力需要実績値の発生頻度、発生確率、又は累積発生確率に関する設定値
である請求項１記載の需要予測処理方法。
【請求項７】
前記需要予測モデル生成ステップにおいて、ランダム法を用いて前記需要予測モデルの係数を算出する
請求項６記載の需要予測処理方法。
【請求項８】
請求項１乃至７のいずれか１つに記載の需要予測処理方法をコンピュータに実行させるためのプログラム。
【請求項９】
電力需要実績値を単位時間毎に格納する電力需要実績データ格納部と、
１又は複数の地点における複数予報時刻の気象予報値を複数シナリオ分配信日毎に格納する気象予報データ格納部と、
ユーザから指定された処理対象配信期間と指定地点と需要予測対象日時と需要予測のための指定配信日とに従って、前記気象予報データ格納部から、前記処理対象配信期間内を配信日としており、前記需要予測対象日時と前記指定配信日の所定配信時刻との差である予報間隔を予報時刻とする、前記指定地点における前記複数シナリオ分の気象予報値を抽出し、抽出気象予報値格納部に格納する手段と、
前記電力需要実績データ格納部から、前記処理対象配信期間内の各日の前記所定配信時刻と前記予報間隔とを加算することによって得られる仮想需要予測日時における電力需要実績値を抽出し、抽出電力需要実績値格納部に格納する手段と、
前記複数シナリオ分の気象予報値の中央値又は平均値によって需要予測値分布の高低を調整する項と、前記気象予報値と前記気象予報値の中央値又は平均値との差によって前記需要予測値分布の幅を調整する項とを有する需要予測モデルの係数を、前記抽出気象予報値格納部に格納されている前記気象予報値と前記需要予測モデルとによって算出される需要予測値の分布が、前記抽出電力需要実績値格納部に格納されている前記電力需要実績値に基づくユーザ設定に従うように算出し、需要予測モデルデータ格納部に格納する需要予測モデル生成手段と、
前記需要予測モデルデータ格納部に格納されている前記需要予測モデルの係数を設定した前記需要予測モデルに対して、前記気象予報データ格納部に格納されており且つ前記指定配信日における前記予報間隔を予報時刻とする、前記指定地点における前記複数シナリオ分の気象予報値を適用して、前記複数セットの需要予測値を算出し、需要予測シナリオデータ格納部に格納する手段と、
を有する需要予測処理装置。
【請求項１０】
発電実績値を単位時間毎に格納する発電実績データ格納部と、１又は複数の地点における複数予報時刻の気象予報値を複数シナリオ分配信日毎に格納する気象予報データ格納部とにアクセス可能なコンピュータにより実行される発電予測処理方法であって、
ユーザから指定された処理対象配信期間と指定地点と発電予測対象日時と発電予測のための指定配信日とに従って、前記気象予報データ格納部から、前記処理対象配信期間内を配信日としており、前記発電予測対象日時と前記指定配信日の所定配信時刻との差である予報間隔を予報時刻とする、前記指定地点における前記複数シナリオ分の気象予報値を抽出し、抽出気象予報値格納部に格納するステップと、
前記発電実績データ格納部から、前記処理対象配信期間内の各日の前記所定配信時刻と前記予報間隔とを加算することによって得られる仮想発電予測日時における発電実績値を抽出し、抽出発電実績値格納部に格納するステップと、
前記複数シナリオ分の気象予報値の中央値又は平均値によって発電予測値分布の高低を調整する項と、前記気象予報値と前記気象予報値の中央値又は平均値との差によって前記発電予測値分布の幅を調整する項とを有する発電予測モデルの係数を、前記抽出気象予報値格納部に格納されている前記気象予報値と前記発電予測モデルとによって算出される発電予測値の分布が、前記抽出発電実績値格納部に格納されている前記発電実績値に基づくユーザ設定に従うように算出し、発電予測モデルデータ格納部に格納する発電予測モデル生成ステップと、
前記発電予測モデルデータ格納部に格納されている前記発電予測モデルの係数を設定した前記発電予測モデルに対して、前記気象予報データ格納部に格納されており且つ前記指定配信日における前記予報間隔を予報時刻とする、前記指定地点における前記複数シナリオ分の気象予報値を適用して、前記複数セットの発電予測値を算出し、発電予測シナリオデータ格納部に格納するステップと、
を含む発電予測処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、電力の需要予測技術及び電力の発電予測技術に関する。
【背景技術】
【０００２】
電力需要の予測は安定供給と経済性の追求のために不可欠である。このため電力会社では当日、翌日の需要や１週間先、１ヶ月先、１年先の需要など、さまざまな時間断面を対象とした予測を行っている。
【０００３】
短期的な時間断面を対象とした需要予測で一般的な方法として、電力需要と気象状況との関係に着目した予測が行われている。すなわち過去の両者の関係を回帰式などで需要モデルとして表し、これに予測対象日の気象に関する予報値を与え予測を行う。ところが、この気象予報には常に誤差（予報誤差）が伴う。このため気象予報に基づいた需要予測には不確実性が含まれることとなり、電力の発電機運用等で問題とされている。
【０００４】
例えば、気象予報が外れて、需要予測値が低めに大きく外れた場合には発電機の予備力が極端に減少し電力安定供給に向けてのマージン低下原因となり問題となる。また、需要予測値が高めに大きく外れた場合には、安定供給上は問題ないものの、需要に対して過剰となった発電機を発電効率の悪い低出力で運転することとなり経済性が悪化する原因となる。
【０００５】
従来の需要予測では気象庁等が配信する単一の気象予報をもとに需要予測を行っていた。このため需要値も一点のみ予測するに限られていた。危機管理と経済性の双方を満足した発電機等の電力設備の運用に向けては、電力需要を従来の１点のみの予測ではなく幅をもって予測することが望ましい。しかしながら、日々の気象予報自身の不確実性を定量的に評価することの困難性もあり、将来の電力需要を幅や確率分布などの形態で客観的に予測するシステムや方法は現在にいたるまで構築に至ってはいない。
【０００６】
近年、従来の単独予報ではなく複数の気象予報シナリオを作成、提供するアンサンブル予報の開発と運用が実施されている（古川、酒井、「アンサンブル予報」、東京堂出版、２００４年）。
【０００７】
なお、気象庁などが実施する気象予報を大別すると天気予報者が気象学的知識や経験に基づいて地上天気図や高層天気図など解析し予報する“主観的予報”、数値気象予報モデルによる“客観予報”がある。
【０００８】
後者の数値気象予報モデルは、基礎的な物理法則（数式で書かれた）に基づいて大気の力学要素である気圧、気温、風速、風向などの時間変化を定量的に算出し、これらの変数の将来の３次元分布を求めるものである。計算では大気を立体格子で区切り、その格子点毎に地上および高層観測などから得た大気の実測データなどから初期値を設定し、計算機上で模擬した物理法則に従い格子毎に各気象変数の将来予測値を算出する。
【０００９】
一方で大気の運動は非線形な特質をもっている。この非線形性のため将来の動きは初期値に非常に敏感な特質を有することとなる。また、予報期間が長くなればなるほど敏感性は大きくなり、わずかな初期値の違いが気象予報の大きな違いとなる。気象予報において予報の土台となる大気の完全な初期値を観測から得ることは不可能であり、初期値に関わる設定誤差は気象予報誤差の一要因となっている。
【００１０】
アンサンブル予報では、初期値に観測誤差と同程度の小さな誤差を与えた多数の初期値を用意する。その後、それぞれの初期値毎に独立に数値予報を行い、異なる大気の将来値を多数得るものである。
【００１１】
例えば特開平６−１０５４６５号公報記載の需要予測装置に、複数の気象予報を入力すれば複数個の需要予測値を得ることは可能である。特にアンサンブル予報など客観的な手法に基づいた予報値を入力情報とすることで日々一貫性のある幅による需要予測が可能となる。
【特許文献１】特開平６−１０５４６５号公報
【発明の開示】
【発明が解決しようとする課題】
【００１２】
しかしながら、アンサンブル予報等を需要予測に利用するためには課題がある。まず、アンサンブル予報が扱う気象と、電力需要と関係する地上の気象（特に気温など）とは異なる点である。アンサンブル予報は大気や地形位置を離散表現した数値気象予報モデルに従い気象予報を実施する。得られる予報値は地表や大気をメッシュ分割した格子点での気象予報値である。したがってアンサンブル予報が対象とする気象と実際の需要地域の気象値とには位置的な違い、さらには気象モデル自体の表現精度などにより必ず差（バイアス）が発生する。また、アンサンブル予報で求められる複数の気象予報シナリオのばらつきに関した課題もある。アンサンブル予報は初期値の与え方により気象予報のばらつきも大きく変わる。たとえばアンサンブル予報を分析したこれまでの報告（「気温と気圧に関する週間アンサンブル気象予報の精度の予報確率の検証」、電力中央研究所報告、Ｎ０５０４８、平成１８年）によれば、得られた予報のばらつきが実際の気象値の範囲より狭いという指摘もある。
【００１３】
このような気象予報値を需要予測に応用した際には結果として得られる需要予測値の幅が狭すぎたり、逆に広すぎたり、さらには高め側や低め側のどちらかに偏りすぎたり、といったバイアスが発生し、単にアンサンブル予報を需要予測の入力情報にするだけでは実務利用には供せない。
【００１４】
そこで、本発明の目的は、アンサンブル予報を基に適切に需要予測又は発電予測を行うための技術を提供することである。
【００１５】
また、本発明の他の目的は、気象予報と電力需要又は発電量との間の各種バイアスを自動補正した上で需要予測又は発電予測を行うための技術を提供することである。
【課題を解決するための手段】
【００１６】
本発明に係る需要予測処理方法は、電力需要実績値を単位時間毎に格納する電力需要実績データ格納部と、１又は複数の地点における複数予報時刻の気象予報値を複数シナリオ分配信日毎に格納する気象予報データ格納部とにアクセス可能なコンピュータにより実行される。そして、（Ａ）ユーザから指定された処理対象配信期間と指定地点と需要予測対象日時と需要予測のための指定配信日とに従って、気象予報データ格納部から、処理対象配信期間内を配信日としており、需要予測対象日時と指定配信日の所定配信時刻との差である予報間隔を予報時刻とする、指定地点における複数シナリオ分の気象予報値を抽出し、抽出気象予報値格納部に格納するステップと、（Ｂ）電力需要実績データ格納部から、処理対象配信期間内の各日の所定配信時刻と予報間隔とを加算することによって得られる仮想需要予測日時における電力需要実績値を抽出し、抽出電力需要実績値格納部に格納するステップと、（Ｃ）複数シナリオ分の気象予報値の中央値又は平均値によって需要予測値分布の高低を調整する項と、気象予報値と当該気象予報値の中央値又は平均値との差によって需要予測値分布の幅を調整する項とを有する需要予測モデルの係数を、抽出気象予報値格納部に格納されている気象予報値と需要予測モデルとによって算出される需要予測値の分布が、抽出電力需要実績値格納部に格納されている電力需要実績値に基づくユーザ設定に従うように算出し、需要予測モデルデータ格納部に格納する需要予測モデル生成ステップと、（Ｄ）需要予測モデルデータ格納部に格納されている需要予測モデルの係数を設定した需要予測モデルに対して、気象予報データ格納部に格納されており且つ指定配信日における予報間隔を予報時刻とする、指定地点における複数シナリオ分の気象予報値を適用して、複数シナリオ分の需要予測値を算出し、需要予測シナリオデータ格納部に格納するステップとを含む。
【００１７】
このようにすれば、使用すべき複数シナリオ分の気象予報値と対応する電力需要実績値とから、ユーザ設定に従って、需要予測値分布についてのバイアスを調整する項を含む需要予測モデルの係数を算出して、適切な複数シナリオ分の需要予測値を算出することができるようになる。
【００１８】
また、上で述べた需要予測値分布の幅を調整する項が、需要予測値分布において中央値又は平均値以上の部分の分布の幅を調整する項と、需要予測値分布において中央値又は平均値より下の部分の分布の幅を調整する項とを含むようにしてもよい。このようにすれば、より適切にバイアスを調整することができるようになる。
【００１９】
さらに、上で述べた需要予測モデルが、気象予報値と当該気象予報値の中央値又は平均値との差の二乗についての項をさらに含むようにしてもよい。このようにすれば、より精密にバイアスを調整することができるようになる。
【００２０】
また、上で述べたユーザ設定が、電力需要実績値の最大値と最小値をカバーする可能な限り幅の狭い分布である場合もある。このようにすれば、算出される複数シナリオ分の需要予測値を、上限値及び下限値を示すものとして解釈して取り扱うことができるようになる。なお、このような場合、上で述べた需要予測モデル生成ステップにおいて、ＬＰ法を用いて需要予測モデルの係数を算出するようにしてもよい。
【００２１】
また、上で述べたユーザ設定が、抽出気象予報値格納部に格納されている気象予報値と需要予測モデルとによって算出される需要予測値の昇順系列又は降順系列を順に区間閾値とする各区間又は特定区間における電力需要実績値の発生頻度、発生確率、又は累積発生確率に関する設定値である場合もある。このようにすれば、例えば、上位特定区間や下位特定区間の発生確率や累積発生確率を指定することによって、ある程度所望の分布形状を有する需要予測モデルを生成することができるようになる。なお、このような場合、上で述べた需要予測モデル生成ステップにおいて、ランダム法（ランダム探索法とも呼ぶ）を用いて需要予測モデルの係数を算出するようにしてもよい。
【００２２】
本発明に係る発電予測処理方法は、発電実績値を単位時間毎に格納する発電実績データ格納部と、１又は複数の地点における複数予報時刻の気象予報値を複数シナリオ分配信日毎に格納する気象予報データ格納部とにアクセス可能なコンピュータにより実行される発電予測処理方法である。そして、（Ａ）ユーザから指定された処理対象配信期間と指定地点と発電予測対象日時と発電予測のための指定配信日とに従って、気象予報データ格納部から、処理対象配信期間内を配信日としており、発電予測対象日時と指定配信日の所定配信時刻との差である予報間隔を予報時刻とする、指定地点における複数シナリオ分の気象予報値を抽出し、抽出気象予報値格納部に格納するステップと、（Ｂ）発電実績データ格納部から、処理対象配信期間内の各日の所定配信時刻と予報間隔とを加算することによって得られる仮想発電予測日時における発電実績値を抽出し、抽出発電実績値格納部に格納するステップと、（Ｃ）複数シナリオ分の気象予報値の中央値又は平均値によって発電予測値分布の高低を調整する項と、気象予報値と気象予報値の中央値又は平均値との差によって発電予測値分布の幅を調整する項とを有する発電予測モデルの係数を、抽出気象予報値格納部に格納されている気象予報値と発電予測モデルとによって算出される発電予測値の分布が、抽出発電実績値格納部に格納されている発電実績値に基づくユーザ設定に従うように算出し、発電予測モデルデータ格納部に格納する発電予測モデル生成ステップと、（Ｄ）発電予測モデルデータ格納部に格納されている発電予測モデルの係数を設定した発電予測モデルに対して、気象予報データ格納部に格納されており且つ指定配信日における予報間隔を予報時刻とする、指定地点における複数シナリオ分の気象予報値を適用して、複数セットの発電予測値を算出し、発電予測シナリオデータ格納部に格納するステップとを含む。
【００２３】
なお、本方法をコンピュータに実行させるためのプログラムを作成することができ、このプログラムは、例えばフレキシブルディスク、ＣＤ−ＲＯＭ、光磁気ディスク、半導体メモリ、ハードディスク等の記憶媒体又は記憶装置に格納される。また、ネットワークなどを介してデジタル信号として配信される場合もある。尚、中間的な処理結果はメインメモリ等の記憶装置に一時保管される。
【発明の効果】
【００２４】
本発明によれば、アンサンブル予報を基に適切に需要予測又は発電予測を行うことができるようになる。
【００２５】
また、本発明の他の側面によれば、気象予報と電力需要又は発電量との間の各種バイアスを自動補正した上で需要予測又は発電予測を行うことができるようになる。
【発明を実施するための最良の形態】
【００２６】
本発明の一実施の形態に係るシステム概要を図１に示す。例えば、インターネットなどのネットワーク１には、アンサンブル予報等のデータ生成及び配信を行う気象予報サーバ３と、本実施の形態に係る主要な処理を実施する電力需要予測装置５とが接続されている。また、電力需要に対する供給源である各発電機の出力（発電量。ここでは需要実績値とする）を所定の時間間隔（例えば３０分や１時間）で計測する装置である電力データ計測装置１１と、当該電力データ計測装置１１に接続されており、各発電機出力を所定時間間隔で合計して、所定時間毎の電力需要実績値を算出する電力計測データ収集装置９も存在している。そして、例えば電力会社のＬＡＮ（Local Area Network）であるネットワーク７を介して、電力計測データ収集装置９は、電力需要予測装置５と接続される。なお、電力計測データ収集装置９及び電力データ計測装置１１は、本実施の形態の主旨ではないので、これ以上述べない。また、ネットワーク１又は７には、電力需要予測装置５のユーザが操作するユーザ端末が接続されることもある。
【００２７】
気象予報サーバ３は、例えば気象予報機関が運用している１台のコンピュータであっても良いし複数のコンピュータで構成されても良いが、数値気象予測モデルに、気象観測値などから作成した初期値を数値気象モデルの格子点毎に与えて、将来の各格子点の大気変量を予測することによって、数値気象予報データを生成する数値気象予報生成部３１と、生成された数値気象予報データを例えば要求に応じてネットワーク１を介して配信する数値気象予報配信部３３とを有する。数値気象予報生成部３１は、模式的に示せば、図２に示すように、地球表面（緯度方向及び経度方向）を所定サイズの格子に分割すると共に、鉛直方向も同様に所定間隔で同様の格子に分割して、各格子点における大気変量（気温、湿度、気圧など）の状態を、将来の時系列データとして予測する。時系列であるから、図２右下のグラフ（縦軸は気温を表し、横軸は時間を表す）に示すように、１つの格子点において１つの大気変量（例えば気温）について、現時点より将来の複数（例えば５１通り）の気象シナリオが生成される。なお、気象シナリオの生成方法については、本実施の形態の主旨ではないので、これ以上述べない。
【００２８】
電力需要予測装置５は、電力計測データ収集装置９から電力需要実績データを取得する電力需要実績データ受信部５１１と気象予報サーバ３の数値気象予報配信部３３から数値気象予報データを取得する数値気象予報データ受信部５１３とを有するインターフェース部５１と、電力需要実績データ受信部５１１が取得した電力需要実績データを格納する電力需要実績データ格納部５３と、数値気象予報データ受信部５１３が取得した数値気象予報データを格納する数値気象予報データ格納部５４と、ユーザからの需要予測設定データの入力を受け付け、一時的に保持する需要予測設定入力部５７と、電力需要実績データ格納部５３と数値気象予報データ格納部５４に格納されているデータを用いて需要予測設定入力部５７によって入力された需要予測設定データに従って需要予測モデル生成処理を実施する需要予測モデル生成部５５と、需要予測モデル生成部５５によって生成された需要予測モデルデータ（需要予測モデルの係数値）を格納する需要予測モデルデータ格納部５９と、数値気象予報データ受信部５１３によって取得された最新の数値気象予報データ（需要予測に用いられる数値気象予報データ）を格納する最新数値気象予報データ格納部６１と、需要予測モデルデータ格納部５９と最新数値気象予報データ格納部６１に格納されているデータを用い、需要予測設定入力部５７によって入力された需要予測設定データに従って複数の需要予測値（需要予測シナリオ）を計算する需要予測シナリオ計算部６３と、需要予測シナリオ計算部６３によって計算された需要予測シナリオデータを格納する需要予測シナリオデータ格納部６５と、需要予測シナリオデータ格納部６５に格納されたデータを出力する処理を実施する需要予測シナリオ出力部６７とを有する。
【００２９】
なお、電力需要実績データ受信部５１１は、例えば定期的に電力計測データ収集装置９からデータを受信して、電力需要実績データ格納部５３に格納する。また、数値気象予報データ受信部５１３も、例えば定期的に気象予報サーバ３からデータを受信して、数値気象予報データ格納部５４に格納する。但し、需要予測設定入力部５７から需要予測に用いる配信日のデータを受け取ると、当該配信日のデータについては、最新数値気象予報データ格納部６１に格納する。
【００３０】
また、需要予測モデル生成部５５は、図３に示すように、需要予測設定入力部５７によって入力された需要予測設定データに従って必要となる過去の気象予報データを数値気象予報データ格納部５４から抽出する過去気象予報データ抽出部５５１と、過去気象予報データ抽出部５５１によって抽出されたデータを格納する抽出過去気象予報データ格納部５５２と、需要予測設定入力部５７によって入力された需要予測設定データに従って必要となる過去の電力需要実績データを抽出過去気象予報データ格納部５５２に格納されているデータに基づき電力需要実績データ格納部５３から抽出する過去電力需要実績データ抽出部５５３と、過去電力需要実績データ抽出部５５３によって抽出されたデータを格納する抽出電力需要実績データ格納部５５４と、抽出過去気象予報データ格納部５５２と抽出電力需要実績データ格納部５５４に格納されたデータを用い、設定入力部５７によって入力された需要予測設定データに従って需要予測モデルの係数値を算出する係数算出部５５５とを有する。係数算出部５５５は、以下で述べるような処理を実施するＬＰ法処理部５５５１とランダム法処理部５５５３とを有する。
【００３１】
以下、図４乃至図２９を用いて、図１及び図３に示した電力需要予測装置５の処理内容を説明する。まず、電力需要実績データ受信部５１１は、ネットワーク７を介して電力需要実績データを取得して、電力需要実績データ格納部５３に格納する（ステップＳ１）。電力需要実績データ格納部５３には、例えば図５のようなデータが格納される。図５の例では、需要日（年月日）、需要時刻（時。ここでは１時間毎）、需要（ＭＷ）が登録されるようになっている。
【００３２】
次に、数値気象予報データ受信部５１３は、気象予報サーバ３の数値気象予報配信部３３から、ネットワーク１を介して数値気象予報データを取得して、数値気象予報データ格納部５４に格納する（ステップＳ３）。例えば、数値気象予報データ格納部５４には、図６に示すような形式でデータが格納される。図６の例では、予報配信日時（配信日、配信時間、予報初期値の日時）毎に、予報シナリオ番号で特定される各予報シナリオのデータが登録される。また、各予報シナリオのデータは、予報対象地点（緯度及び経度）と予報対象高度（垂直高度又は気圧面）とで特定される各格子点のデータが登録される。さらに、各格子点のデータは、予報対象時刻（予報初期値からの時間間隔（予報間隔とも呼ぶ））毎に、気象要素（気温、相対湿度、降水量、風速など）が登録されるようになっている。このように階層化されて、以下で述べる抽出処理で、必要なデータが抽出される。
【００３３】
そして、需要予測設定入力部５７は、需要予測の結果を欲するユーザから計算条件等の需要予測設定データの入力を受け付け、例えばメインメモリなどの記憶装置に格納する（ステップＳ５）。需要予測設定データは、（１）需要予測モデルの種類、（２）需要予測モデルで利用する気象予報変数（気温や湿度など）、（３）需要予測モデルで利用する格子点、（４）需要予測モデルの係数算出に用いる過去データ範囲、（５）需要予測対象日時及び需要予測に利用する気象予報の配信日、（６）需要予測モデルの係数算出アルゴリズム及び当該係数算出アルゴリズムの制御パラメータなどを含む。
【００３４】
ここで、本実施の形態においてユーザが指定する（１）需要予測モデルの種類について説明しておく。ここでは、例えば以下の（１）乃至（６）式で示される需要予測モデルが選択される。
【００３５】
【数１】
2009225550.tif 000003
なお、ＬＯＡＤ（ｔ，ｉ）は、気象予報毎にｎ個得られる気象予報シナリオの第ｉ番目の予報から得られた需要予測値を表す。なお、１≦ｉ≦ｎである。また、ｔは、予測対象時刻を表す。予報配信日時から予報対象時刻までの時間間隔で指定される。ｍは、需要予測のために選んだ格子点の数を表す。Ｘk,iは、格子点ｋの気象要素（気温など）の第ｉ番目の気象予報シナリオでの値を表す。Ｘ＾k（＾はＸの上のバーを表す）は、格子点ｋの気象要素（気温など）の全気象予報シナリオ（ｎ個）の中央値（又は平均値。以下、同じ。）を表す。ΔＸk,iは、上記の中央値と、第ｉ番目の気象予報シナリオの値との差（Ｘk,i−Ｘ＾k）を表す。Ｓ+は、ΔＸk,iが正（０を含む）の時に１となり、負の時０となるスイッチング変数を表す。Ｓ-は、ΔＸk,iが正（０を含む）の時に０となり、負の時に１となるスイッチング変数を表す。ａ，ｂk，ｃk，ｄk，ｅk，ｆkは、以下の処理で算出する未知係数を表す。ＮＮ１tは、入力をＸ^k,i（ｋ＝１，．．ｍ）及びΔＸk,i（ｋ＝１，．．ｍ）、出力をＬＯＡＤ（ｔ，ｉ）とする階層型ニューラルネットワーク、ＮＮ２tは、入力をＸk,i（ｋ＝１，．．ｍ）及びその中央値、出力をＬＯＡＤ（ｔ，ｉ）とする階層型ニューラルネットワークを表す。
【００３６】
本発明の課題の欄でも説明したが、アンサンブル予報の予報値と電力需要との間のバイアスが問題となる。バイアスには大気モデル（格子）と需要予測対象である電力供給エリアの位置的なずれ、アンサンブル予報自体の変動幅の偏りなどがある。図７を用いて、電力需要予測から見たアンサンブル予報のバイアスの影響を説明する。図７では、アンサンブル予報に基づいて、予測計算実施時点である時刻ｔ0から複数の需要予測値のシーケンスを計算している。
【００３７】
ここで、時刻ｔ1、ｔ2、ｔ3で需要実績値が「×」印で表される値になったとする。この場合、時刻ｔ2では、需要実績値は需要予測シナリオの範囲内に位置している。すなわち、複数の需要予測シナリオが、需要実績値を、範囲内に捉えた結果となっている。これに対して、時刻ｔ1、ｔ3では需要実績値が需要予測シナリオの範囲から上側、下側に大きく外れている。言い換えるならば、時刻ｔ1については需要予測シナリオの幅は狭すぎ、又は下側にシフトしているといえる。時刻ｔ3については需要予測シナリオの幅が狭すぎ、又は上側にシフトしているということになる。
【００３８】
これらのずれは、数値気象予報データを生成する数値気象予報生成部３１側等に由来するものである。こうしたずれが系統的に発生すると想定して、本実施の形態ではこの気象予報の系統的なずれによる需要予測への影響を需要予測モデルによって補正し、需要予測シナリオとして予測幅による需要予測が適切に可能となるよう需要予測モデルを生成する。このため上で述べた（２）乃至（４）及び（５）式では、従来の需要予測モデルとは異なり気象予報のバイアスを補正するための項を特別に設けている。
【００３９】
ここで（３）式の需要予測モデルを例にしてバイアスの補正について具体的に説明する。図８では、縦軸は電力需要を表し、横軸は時間を表す。図８の例では、実線ａで表すように電力需要が推移しているが、バイアス補正を行う前の電力需要予測式で算出した将来の電力需要予測が分布Ａで表され、一方、実際の電力需要の大きさとその大きさの需要が発生する確率からなる実際の分布が分布Ｂで表されるものとする。ここで、（３）式の第１項及び第２項は、電力需要とアンサンブル予報が対象とする気象（気温など）との一般的な関係を捉えるための項であり、格子点の気象と電力需要との位置的関係の違いなどによるバイアスを調整するためのものである。そして、図８では、Ｐで表すように、分布の中央値部分同士のバイアスを補正するように作用する。また、（３）式の第３項は、アンサンブル予報による需要予測分布における中央値より上側のばらつきの大きさを補正するための項であり、図８では、Ｑに示すように、分布の中央値より上側の幅を調整するように作用する。例えば、アンサンブル予報による需要予測分布Ａにおいて中央値より上側に十分な広がりが存在しない場合には、本項のｃkを適切な値に設定することにより、需要予測シナリオにおいても十分な幅を有するようになる。さらに、（３）式の第４項は、アンサンブル予報による需要予測分布Ａにおける中央値より下側のばらつきの大きさを補正するための項であり、図８では、Ｒに示すように、分布の中央値より下側の幅を調整するように作用する。すなわち、係数ｄkにより需要予測シナリオにおいて中央値より下側の分布の幅を実際の電力需要分布Ｂにおける幅に合わせることができるようになる。
【００４０】
このようにして、（３）式によれば、図９に示すように、実際の電力需要分布Ｂに重なるように、需要予測シナリオのバイアスを調整することができる。
【００４１】
なお、（２）式については、（３）式の第２項及び第３項を統合して、全体の幅という観点でバイアスを調整する需要予測式である。
【００４２】
また、（４）式については、ΔＸ2k,iを用いる２つの項を（３）式に追加しており、上で述べたバイアスの調整をより精密に実施するようにする式である。
【００４３】
（１）式は、気象要素（気温など）と該当する係数の積和をベースとした基本的な需要予測モデルであり、バイアスの調整は含まれていない。以下では、（１）式については、上で述べた本実施の形態の特徴を有しないので、特に触れないものとする。
【００４４】
また、（５）式及び（６）式についても、本実施の形態とは別の従来技術を用いて階層型ニューラルネットワークを構築しておくものとして、以下では詳述しない。但し、（５）式、（６）式で示すように、気象予報シナリオの中央値又は平均値、同中央値又は平均値からの各予報シナリオの気象要素の偏差値を説明変数とするものについては、従来とは明らかに異なる。
【００４５】
上で述べた（３）式のような需要予測式は一例であって、本実施の形態の特徴は、需要予測モデルを、各気象予報シナリオの気象要素の値と、その中央値又は平均値を説明変数とすることである。また、幅の広がりの大きさに関するバイアスを調整するための各気象予報シナリオの気象要素の値の中央値又は平均値からの偏差値を説明変数とすることである。なお、偏差値についてであるが（５）（６）式のような内部演算処理機能をもつ予測式では入力情報としては与えずに、内部処理により偏差値を作成することも可能である。（６）式はこれに相当する。式で表すと以下のとおりである。
【００４６】
【数２】
2009225550.tif 000004
なお、ｆtは、需要予測関数を表す。Ｗt,iは、予想対象時刻ｔにおける、アンサンブル予報の第ｉ番目のメンバーのベクトルであり、指定される気象要素及び格子点の値を要素として含む。図１０の行列のｉ番目の行ベクトルに相当する。Ｗt＾（＾はＷの上のバーを表す）は、予測対象時刻ｔにおける、指定気象要素及び格子点毎の気象予報値の中央値又は平均値からなるベクトルである。具体的には、図１１に示すようなベクトルであり、図１１ではｍｅｄｉａｎは中央値を表す関数である。
【００４７】
なお、電力需要は、気象要素以外に曜日の違い（平日、土曜、日曜）などの影響も受ける。これらについては、以下で述べる過去データの抽出処理で予測対象日と同じ曜日の気象予報値及び需要実績値を抽出するようにしたり、需要予測式において曜日の違いを表すフラグ変数を追加するようにしても良い。これらについては、従来から行われており、ここでは説明を省略する。
【００４８】
需要予測設定入力部５７によって入力されるデータの説明に戻って、（２）需要予測モデルで利用する気象予報変数については、気温や湿度などを指定する。
【００４９】
（３）需要予測モデルで利用する格子点については、緯度、経度、高度等によって、１又は複数の格子点を指定する。
【００５０】
（４）需要予測モデルの係数算出に用いる過去データ範囲については、例えば開始年月日及び終了年月日を、配信日として指定する。その他の指定方法であっても良い。
【００５１】
（５）需要予測対象日時及び需要予測に利用する気象予報の配信日については、電力需要予測対象日及び時刻と、その際に用いる気象予報の配信日とを指定する。
【００５２】
（６）需要予測モデルの係数算出アルゴリズム及び当該係数算出アルゴリズムの制御パラメータなどについては、例えばＬＰ法とランダム法のいずれかと、以下で述べる詳細な計算条件（目的関数、制約条件、評価関数、計算終了条件など）を入力するものである。
【００５３】
図４の説明に戻って、次に、需要予測モデル生成部５５は、需要予測モデル生成処理を実施する（ステップＳ７）。この需要予測モデル生成処理の詳細については、以下で述べるが、本ステップによって、ユーザによって指定された需要予測モデルにおける係数値が算出され、需要予測モデルデータ格納部５９に格納される。
【００５４】
そして、数値気象予報データ受信部５１３は、需要予測設定入力部５７によって入力された需要予測に利用する気象予報の配信日のデータに従って、該当日の気象予報データを気象予報サーバ３の数値気象予報配信部３３から取得し、最新数値気象予報データとして最新数値気象予報データ格納部６１に格納する（ステップＳ９）。
【００５５】
その後、需要予測シナリオ計算部６３は、需要予測モデルデータ格納部５９に格納されている係数値によって特定され且つユーザによって指定された種類の需要予測モデルと、最新数値気象予報データ格納部６１に格納されている最新数値気象予報データとから、需要予測シナリオを算出し、需要予測シナリオデータ格納部６５に格納する（ステップＳ１１）。具体的には例えば（１）乃至（６）式のうちいずれか指定された式に従って、各気象予報シナリオに応じて需要予測値を算出する。なお、需要予測対象日時を複数にして、それぞれについて各気象予報シナリオに応じた需要予測値を算出するようにしても良い。
【００５６】
そして、需要予測シナリオ出力部６７は、需要予測シナリオデータ格納部６５に格納されているデータを、ユーザに対して出力する（ステップＳ１３）。例えば、電力需要予測装置５に表示装置が接続されている場合には、当該表示装置に表示したり、印刷装置などに出力するようにしても良い。さらに、電力需要予測装置５が接続されているネットワーク経由でユーザ端末に出力するようにしても良い。表示例については、以下で述べる需要予測モデル生成処理に関連する部分があるので、需要予測モデル生成処理を説明した後に説明することとする。
【００５７】
次に、図１２乃至図２７を用いて、需要予測モデル生成処理について詳細に説明する。この処理においては、まず、過去気象予報データ抽出部５５１は、過去気象予報データ抽出処理を実施する（図１２：ステップＳ２１）。
【００５８】
過去気象予報データ抽出処理については、図１３を用いて説明する。過去気象予報データ抽出部５５１は、指定された気象予報変数情報を検索条件１として設定する（図１３：ステップＳ３１）。上で述べた（２）で指定された気温などが検索条件１として設定される。
【００５９】
次に、過去気象予報データ抽出部５５１は、指定された格子点情報を検索条件２として設定する（ステップＳ３３）。上で述べた（３）で指定された格子点（緯度・経度、高度などで指定）が検索条件２として設定される。
【００６０】
さらに、過去気象予報データ抽出部５５１は、気象予報シナリオ過去データ範囲の取得開始日と取得終了日を検索条件３として設定する（ステップＳ３５）。上で述べた（４）で指定された取得開始日と取得終了日が検索条件３として設定される。
【００６１】
そして、過去気象予報データ抽出部５５１は、需要予測対象日時及び利用する気象予報シナリオ配信日から予報間隔を算出し、検索条件４として設定する（ステップＳ３７）。図１４を用いて予報間隔について説明する。図１４においては、左から右に時間が流れているものとする。現在Ｔ１のタイミングであるとして、Ｔ３は需要予測対象日時として指定される。また、Ｔ２が気象予報シナリオ配信日（具体的には一日における配信時刻）として指定される。そうすると、需要予測対象日時における需要予測を行うためには、気象予報シナリオ配信日Ｔ２において配信される一連の将来時刻の気象予報のうち（Ｔ３−Ｔ２）後の気象予報を取得する必要がある。もし該当する時刻の気象予報が存在しない場合には、最も近い時刻の気象予報を特定する必要がある。この（Ｔ３−Ｔ２）を、予報間隔とする。以下の需要予測モデル生成処理においても、気象予報シナリオ過去データ範囲（取得開始日から取得終了日までの各日）を配信日として、予報間隔後の時刻における気象予報シナリオのデータを用いることによって、適切な需要予測モデルの係数値を算出する。従って、この予報間隔を検索条件４とする。
【００６２】
その後、過去気象予報データ抽出部５５１は、検索条件１乃至４をＡＮＤ条件として、数値気象予報データ格納部５４に格納されているデータから、需要予測モデル生成に利用する気象予報データを抽出して、抽出過去気象予報データ格納部５５２に格納する（ステップＳ３９）。
【００６３】
このようにして、需要予測モデル生成に利用する気象予報データのみが抽出されることとなる。
【００６４】
図１２の説明に戻って、次に、過去電力需要実績データ抽出部５５３は、過去電力需要実績データ抽出処理を実施する（ステップＳ２３）。この過去電力需要実績データ抽出処理については図１５を用いて説明する。
【００６５】
まず、過去電力需要実績データ抽出部５５３は、抽出過去気象予報データ格納部５５２に格納されている気象予報データを用いて予測対象日時を算出する（図１５：ステップＳ４１）。この処理については、図１６を用いて説明する。図１６のテーブルにおいて、予報配信日の欄には、気象予報シナリオ過去データ範囲の開始日から終了日まで列挙されており、予報対象日時の欄には、予報間隔（例えば＋１２０時間）後の日時が登録されている。すなわち、予報配信日の固定の予報配信時刻に予報間隔を加算して得られる日時を算出する。この日時が、ユーザによって需要予測対象日時及び気象予報シナリオ配信日で指定された予測態様に従った過去の予測対象日時となる。このようにして得られた過去の予測対象日時に対応する電力需要実績データを抽出する。
【００６６】
よって、過去電力需要実績データ抽出部５５３は、電力需要実績データ格納部５３から、過去予測対象日時に対応する電力需要実績データを抽出し、抽出電力需要実績データ格納部５５４に格納する（ステップＳ４３）。
【００６７】
このようにして、電力需要実績データについても必要なデータのみが抽出されることとなる。
【００６８】
図１２の説明に戻って、次に、係数算出部５５５は、気象予報データ加工処理を実施する（ステップＳ２５）。この処理は、需要予測モデル（２）乃至（４）で頻繁に用いられる中央値Ｘk＾及び中央値からの偏差ΔＸk,iを予め算出する処理である。
【００６９】
この気象予報データ加工処理については、図１７乃至図１９を用いて説明する。図１７で示すように、係数算出部５５５は、中央値算出処理を実施し、処理結果を抽出過去気象予報データ格納部５５２に格納する（ステップＳ５１）。さらに、偏差算出処理を実施し、処理結果を抽出過去気象予報データ格納部５５２に格納する（ステップＳ５３）。
【００７０】
図１８を用いて中央値算出処理について説明する。係数算出部５５５は、抽出過去気象予報データ格納部５５２を参照して、需要予測モデル生成で利用する気象予報の格子点の数をｍに代入する（ステップＳ６１）。また、抽出過去気象予報データ格納部５５２を参照して、抽出されている気象予報データにおける総シナリオ数をｎに代入する（ステップＳ６３）。さらに、カウンタｋを１に初期化する（ステップＳ６５）。
【００７１】
そして、係数算出部５５５は、抽出過去気象予報データ格納部５５２に格納されている抽出された過去気象予報データから格子点ｋの全気象予報シナリオ値｛Ｘk,1，Ｘk,2，．．Ｘk,n｝を抽出する（ステップＳ６７）。また、抽出した気象予報シナリオ値を昇順にソートする（ステップＳ６９）。その結果として、ｎが奇数の場合には（ｎ＋１）／２番目の値を中央値として特定し、ｎが偶数の場合にはｎ／２＋１番目の値を格子点ｋの中央値として特定し、抽出過去気象予報データ格納部５５２に格納する（ステップＳ７１）。
【００７２】
その後、係数算出部５５５は、ｋを１インクリメントして（ステップＳ７３）、ｋ≦ｍであるか判断する（ステップＳ７５）。この条件を満たす場合にはステップＳ６７に戻る。一方、この条件を満たさなくなった場合には、元の処理に戻る。なお、これらの一連の処理は抽出した全ての日のデータに対して日毎に実施するものである。
【００７３】
さらに、図１９を用いて偏差算出処理について説明する。係数算出部５５５は、カウンタｋを１に初期化し（ステップＳ８１）、カウンタｉを１に初期化する（ステップＳ８３）。そして、抽出された過去気象予報データから格子点ｋのｉ番目の気象予報シナリオ値を抽出する（ステップＳ８５）。そして、抽出過去気象予報データ格納部５５２から格子点ｋの中央値を取得し（ステップＳ８７）、ｉ番目の気象予報シナリオ値と中央値との差を、偏差ΔＸk,iとして算出し、抽出過去気象予報データ格納部５５２に格納する（ステップＳ８９）。その後、ｉを１インクリメントする（ステップＳ９１）。そして、ｉ≦ｎであるか判断する（ステップＳ９３）。ｎは総シナリオ数である。ｉがこの条件を満たしている場合には、ステップＳ８５に戻る。一方、ｉがこの条件を満たさなくなった場合には、ｋを１をインクリメントする（ステップＳ９５）。そして、ｋ≦ｍであるか判断する（ステップＳ９７）。この条件を満たす場合にはステップＳ８３に戻る。一方、この条件を満たさなくなった場合には元の処理に戻る。なお、これらの一連の処理は抽出した全ての日のデータに対して日毎に実施するものである。
【００７４】
図１２の説明に戻って、係数算出部５５５は、引き続き需要予測モデル生成処理を実施する（ステップＳ２７）。この需要予測モデル生成処理については、図２０乃至図２７を用いて説明する。
【００７５】
まず、係数算出部５５５は、指定需要予測モデル種別（例えば（３）式など）を特定する（ステップＳ１０１）。また、指定処理アルゴリズム（例えばＬＰ法又はランダム法）を特定する（ステップＳ１０３）。さらに、指定制御パラメータなどを特定する（ステップＳ１０５）。目的関数、制約条件、終了条件等各種の制御パラメータを特定する。そして、指定処理アルゴリズムがＬＰ法であれば（ステップＳ１０６：Ｙｅｓルート）、抽出過去気象予報データ格納部５５２から気象予報データを読み出し、抽出電力需要実績データ格納部５５４から需要実績データを読み出し（ステップＳ１０７）、ＬＰ法による需要予測モデル係数推定処理を実施する（ステップＳ１０９）。
【００７６】
ＬＰ法による係数推定法では、以下で述べるように、制約条件によって需要実績値が需要予測範囲を超えることの無いように需要予測モデルの係数を推定するものである。すなわち、図７の時刻ｔ1及びｔ3のように需要予測範囲の幅が狭くなりすぎないように係数を推定している。また、以下で述べるような目的関数を採用して、目的関数を最小化するということは、需要予測の振れ幅を、制約条件の範囲内でできるだけ小さく狭めるように働く。従って、トータルとしてＬＰ法を採用することによって、需要予測幅の上下限が需要実績が発生する範囲の上下限幅と可能な限り一致するような需要予測モデルが求められる。このようにして生成された需要予測モデルを用いて予測計算を実施した場合には、将来の需要予測値の高め側最大値、低め側最小値の目安情報を得ることができるようになる。
【００７７】
一方、指定処理アルゴリズムがランダム法であれば（ステップＳ１０６：Ｎｏルート）、係数算出部５５５は、抽出過去気象予報データ格納部５５２から気象予報データを読み出し、抽出電力需要実績データ格納部５５４から需要実績データを読み出し（ステップＳ１１１）、ランダム法による需要予測モデル係数推定処理を実施する（ステップＳ１１３）。
【００７８】
ランダム法による係数推定法では、生成された需要予測モデルによって得られるそれぞれの需要予測シナリオに確率を与えることが可能となる。例えば、需要予測モデルを生成する際には、需要予測値の上側５番目、下側５番目が、それぞれ上側１０％、下側１０％の確率点となるように目的関数を設定すれば、生成された需要予測モデルから算出される需要予測値は、設定に従った分布を有するようになる。
【００７９】
上で述べた処理フローでは、ＬＰ法又はランダム法のいずれかを実施することになっているが、両方を実施するようにしても良い。
【００８０】
以下、それぞれの方法について詳述する。
【００８１】
（Ａ）ＬＰ法
ＬＰ法による需要予測モデル係数推定処理について図２１及び図２２を用いて説明する。まず、係数算出部５５５のＬＰ法処理部５５５１は、需要予測設定入力部５７によって入力されたデータに従って、指定需要予測モデルに対応する目的関数を設定する（ステップＳ１２１）。例えば、（３）式の需要予測モデルを採用する場合には、以下のような目的関数を設定し、Ｊを最小化する。
【００８２】
【数３】
2009225550.tif 000005
Ｊを最小化するということは、需要予測の幅をできる限り狭くするように働く。
【００８３】
（２）式や（４）式を採用する場合には、それぞれの式の気象シナリオの偏差に関わる係数であるｂk，ｃk，ｄk，ｅk，ｆkを変数としてこの目的関数の形を同じように変えればよい。
【００８４】
次に、ＬＰ法処理部５５５１は、需要予測設定入力部５７によって入力されたデータに従って、指定需要予測モデルに対する制約条件１及び２を設定する（ステップＳ１２３）。例えば、（３）式を採用する場合には、以下のような制約条件を設定する。
【００８５】
【数４】
2009225550.tif 000006
【００８６】
【数５】
2009225550.tif 000007
【００８７】
（８）式は、過去の気象予報と需要予測値の関係において、日（ｄ）毎に複数得られる需要予測値の最大値を実際の需要実績値Ｙ（ｄ）が超過することを許さないという制約に相当する。ΔＸk,imax_d(d)の添え字ｉmax_dは、ｄ日において需要予測値が最大となるであろう需要予測シナリオ（対応する気象予測シナリオ）の番号を表す。この番号の与え方については、以下で述べる。
【００８８】
（９）式は、過去の気象予報と需要予測値の関係において、日（ｄ）毎に複数得られる需要予測値の最小値を実際の需要予測値Ｙ（ｄ）が下回ることを許さないという制約に相当する。ΔＸk,imin_d(d)の添え字ｉmin_dは、ｄ日において需要予測値が最小となるであろう需要予測シナリオ（対応する気象予測シナリオ）の番号を表す。この番号の与え方については、以下で述べる。
【００８９】
（８）及び（９）式は、抽出された需要実績値の各日についての制約式となる。従って、例えば３０日分の過去データを抽出した場合には、（８）式及び（９）式それぞれについて３０本ずつとなる。
【００９０】
（２）式や（４）式を採用する場合には、この制約条件の値を同じように変えればよい。
【００９１】
さらに、ＬＰ法処理法５５５１は、制約条件３を、ａ≧０に設定する（ステップＳ１２５）。ＬＰ法では、係数値は、正又は負（処理上は負の場合には係数にマイナスをかけ、正の条件にして解く）のいずれか一方でなければならない。このため、本処理では、正及び負のいずれについても処理を実施し、目的関数の値が最も小さくなる方の係数値を需要予測式の係数とする。ここでは、正とする。
【００９２】
その後、ＬＰ法処理部５５５１は、制約条件１の２つのインデックスｉmax_d及びｉmin_dを仮設定する（ステップＳ１２７）。本来ｉmax_d及びｉmin_dは需要予測モデルの係数が得られた後に定まる値である。このため、本処理では、まず初期値として各日の最大需要予測値、最小需要予測値を算出する需要予測シナリオ（気象予報シナリオに対応する需要予測値）の番号を仮設定する。具体的には、夏季において気温が高いほど高い需要が、低いほど低い需要が発生するため、最も高気温の気象予報シナリオの番号をｉmax_dと設定し、最も低い気温の気象予報シナリオの番号をｉmin_dと設定する。なお、格子点として複数地点を選んだ場合には、以下で述べる式のように、選んだ地点の気象予報シナリオにおける気象要素（気温）を平均し、平均後に最も高気温の気象予報シナリオの番号をｉmax_dと設定し、最も低気温の気象予報シナリオの番号をｉmin_dと設定する。
【００９３】
【数６】
2009225550.tif 000008
【００９４】
仮設定に基づき、ＬＰ問題として係数ａ，ｂk，ｃk，ｄkを求める。その後、求められた係数値を基に各日の需要予測値の中から最大値及び最小値を特定する。この最大値及び最小値に対応した需要予測シナリオを新たにｉmax_d及びｉmin_dとおき、再度ＬＰ問題として係数ａ，ｂk，ｃk，ｄkを求める。これらを繰り返し制約条件を満たし且つ目的関数を最小化する係数を求めて行く。なお、ｉmax_d及びｉmin_dを設定した後の係数推定処理には単体法（関根、「数理計画法」、岩波書店、１９７６年を参照のこと。）などのＬＰ専用手法を利用する。
【００９５】
そして、ＬＰ法処理部５５５１は、目的関数及び制約条件を単体法ルーチンで処理して係数を算出し、メインメモリなどの記憶装置に格納する（ステップＳ１２９）。単体法ルーチンについては、周知であるからここでは説明を省略する。
【００９６】
その後、ＬＰ法処理部５５５１は、算出された係数を用いて指定需要予測モデル式で各日の需要予測値を気象予報シナリオ毎に算出し、メインメモリなどの記憶装置に格納する（ステップＳ１３１）。そして、各日について最大需要予測値及び最小需要予測値を特定する（ステップＳ１３３）。さらに、各日について制約条件１（すなわち、（８）及び（９）式）が成立しているか判断する（ステップＳ１３５）。処理は端子Ａを介して図２２の処理に移行する。
【００９７】
図２２の処理の説明に移行して、ＬＰ法処理部５５５１は、ステップＳ１３３で特定された結果が制約条件１を逸脱したか判断する（ステップＳ１３７）。逸脱している場合には、各日について、最大需要予測を生じさせた気象予報シナリオの番号でｉmax_dを、最小需要予測を生じさせた気象予報シナリオの番号でｉmin_dを設定する（ステップ１３９）。そして処理は端子Ｂを介してステップＳ１２９に戻る。
【００９８】
一方、制約条件１を逸脱していないと判断された場合には、ＬＰ法処理部５５５１は、この場合の目的関数値Ｊ及び係数の値を、需要予測モデルデータ格納部５９に仮に格納する（ステップＳ１４１）。そして、制約条件３が、ａ＜０となっているか判断する（ステップＳ１４３）。ａ＜０となっていない場合には、制約条件３をａ＜０に設定して（ステップＳ１４５）端子Ｃを介してステップＳ１２７に戻る。
【００９９】
一方、制約条件がａ＜０に設定されている場合には、ａについて正の場合及び負の場合について係数の値が算出されたことになるので、ＬＰ法処理部５５５１は、ａが正の場合と負の場合の目的関数値Ｊを比較して、小さい方の場合の係数の値を、採用する係数の値として需要予測モデルデータ格納部５９に格納する（ステップＳ１４７）。
【０１００】
このような処理を実施することによって、上で述べたように、需要予測幅が需要実績が発生する範囲と一致するような需要予測モデルが生成されたこととなる。これによって、需要予測を行った際に得られた需要予測値の最大値及び最小値は、需要予測対象日時における高め側最大値、低め側最小値といった情報に相当するものとして取り扱うことができる。
【０１０１】
（Ｂ）ランダム法
次に、図２３乃至図２９を用いて、ランダム法による需要予測モデル係数推定処理について説明する。なお、ランダム法は、ランダム探索（馬場、坂和、「数理計画法入門」、共立出版社、１９８９年などを参照のこと）による方法である。
【０１０２】
ランダム法を具体的に述べる前に、ランダム法で使用される発生確率、累積発生確率について説明しておく。
【０１０３】
本実施の形態で利用する需要予測モデルは電力需要が発生する範囲を予測するものである。これに対して将来需要の期待値を予測することを目的とした従来方式で主に用いられている最小二乗法などの手法をそのまま適用することはできない。なお、最小二乗法では予測モデルから得られる予測値と回帰データとして用いる実績値との差（残差）を特徴量として、この残差の２乗和が最小となるように予測モデルの係数を推定する。
【０１０４】
本実施の形態では従来の「残差」やその「２乗和」といった尺度ではなく、「発生度数」「発生確率」「累積発生確率」などの指標を需要予測モデルの係数推定のために新たに定義し、その上で過去データを用いて需要予測モデルの係数を求める。
【０１０５】
これらの指標については前処理として図２３で説明するように算出される。まず、ランダム法処理部５５５３は、抽出過去気象予報データ格納部５５２に格納されている気象予報データの中で最も過去の日の気象予報シナリオｎ個を抽出すると共に、該当日の需要実績値を抽出電力需要実績データ格納部５５４から抽出する（ステップＳ２０１）。この段階では、例えば図２４（ａ）のように、ｎ個の気象予報シナリオが抽出される。
【０１０６】
そして、ランダム法処理部５５５３は、気象予報シナリオｎ個を需要予測モデル（ここでは係数値が決定されているものとする）に対して適用して、ｎ個の需要予測値を算出し、例えばメインメモリなどの記憶装置に格納する（ステップＳ２０３）。この段階では、例えば図２４（ｂ）のように、ｎ個の需要予測値が算出される。そして、算出されたｎ個の需要予測値をソートする（ステップＳ２０５）。この段階では、図２４（ｃ）に示すように降順にソートされる。
【０１０７】
さらに、ランダム法処理部５５５３は、抽出した需要実績値と、ソート後の需要予測値とを比較する（ステップＳ２０７）。例えば図２４（ｄ）に示すように、需要実績値が７９５０であれば、図２４（ｃ）に示したソート後の需要予測値と比較すると、第２順位の８１５２と、第３順位の７８９５の間であることが分かる。
【０１０８】
そして、ランダム法処理部５５５３は、抽出された需要実績値の発生位置を特定し、発生頻度分布テーブル（図２４（ｅ））における度数を更新する（ステップＳ２０９）。ステップＳ２０７では、第２順位と第３順位の間に需要実績値が該当するので、発生位置は３と特定され、図２４（ｅ）に示すように、発生位置３の度数が１インクリメントされる。図２４（ｅ）の例では２０が２１になる。なお、発生位置は、第１順位より大きい場合に発生位置１、第１順位と第２順位の間の場合に発生位置２、．．．第ｎ−１順位と第ｎ順位の間の場合に発生位置ｎ、第ｎ順位より小さい場合に発生位置ｎ＋１となる。
【０１０９】
その後、ランダム法処理部５５５３は、抽出過去気象予報データ格納部５５２に格納されているデータにおいて全ての日について処理したか判断する（ステップＳ２１１）。未処理の日が存在する場合には、処理日を１日後に設定し（ステップＳ２１３）、処理日の気象予報シナリオｎ個を抽出過去気象予報データ格納部５５２から読み出し、処理日の需要実績値を抽出電力需要実績データ格納部５５４から抽出する（ステップＳ２１５）。そしてステップＳ２０３に戻る。
【０１１０】
例えば図２４（ｆ）のように、ｎ個の気象予報シナリオが新たに抽出された場合には、図２４（ｇ）のように、各気象予報シナリオに対して需要予測値が算出される。そして、図２４（ｈ）のように、需要予測値をソートして、図２４（ｉ）に示すような需要実績値と比較する。この場合、需要予測値３８００は、第１順位より大きいので、発生位置１と特定され、図２４（ｊ）に示すように、発生位置１の度数を１インクリメントする。
【０１１１】
このような処理を繰り返すと、発生位置毎に度数が特定される。
【０１１２】
そして、全ての日について処理した場合には、ランダム法処理部５５５３は、発生頻度分布テーブルにおいて、発生確率、上側累積発生確率、下側累積発生確率を算出し、登録する（ステップＳ２１７）。
【０１１３】
発生確率Ｐiは、以下の（１１）式で表される。
【０１１４】
【数７】
2009225550.tif 000009
Ｘiは発生頻度テーブルにおける発生位置ｉの度数（発生頻度）である。
【０１１５】
上側累積発生確率Ｕiは、以下の（１２）式で表される。
【０１１６】
【数８】
2009225550.tif 000010
【０１１７】
下側累積発生確率Ｌiは、以下の（１３）式で表される。
【０１１８】
【数９】
2009225550.tif 000011
【０１１９】
ステップＳ２１７を実行することによって、図２５に示されるような発生頻度分布テーブルが生成され、例えば需要予測モデルデータ格納部５９に格納される。図２５の例では、発生位置毎に、発生度数Ｘi、発生確率Ｐi、上側累積発生確率Ｕi、下側累積発生確率Ｌiが登録されるようになっている。
【０１２０】
このような前提の下、ランダム法による需要予測モデル係数推定処理を説明する。ランダム法処理部５５５３は、需要予測設定入力部５７によって入力されたデータに従って、目的関数を設定する（図２６：ステップＳ１５１）。
【０１２１】
発生確率についての目的関数を採用する場合には、例えば以下のような目的関数を設定する。
【０１２２】
【数１０】
2009225550.tif 000012
Ｐiは、需要実績値が、ソート後の需要予測シナリオにおける発生位置ｉに発生した確率を表す。Ｐtarget,iは、上記発生確率の目標値（例えば、Ｐtarget,i＝１／（ｎ＋１）×１００）を表す。ｗiは、重み係数を表す。なお、Ｐtarget,i＝１／（ｎ＋１）×１００は、需要予測モデルから得られる需要予測シナリオの発生確率が等しくなるように需要予測モデルの係数を求める場合に設定される。また、この際ｗiはｉによらず固定値（例えば１）とする。
【０１２３】
（１４）式を用いるのは、需要予測モデルで得られる需要予測値の確率分布を、目標の形状となるようにする場合である。
【０１２４】
一方、累積確率についての目的関数を採用する場合には、例えば以下のような目的関数を設定する。
【０１２５】
【数１１】
2009225550.tif 000013
Ｕiは、需要実績値がソート後の需要予測シナリオにおける上位ｉ番目の値を超えた確率を表す。Ｌiは、需要実績値がソート後の需要予測シナリオにおける下位ｉ−１番目の値以下となった確率を表す。Ｕtarget,iは、Ｕiの目標値を表す。さらに、Ｌtarget,iは、Ｌiの目標値を表す。ｗi及びｖiは重み係数を表す。
【０１２６】
例えばＵtarget,iは、以下のように設定される。
【０１２７】
【数１２】
2009225550.tif 000014
ｎはシナリオの数である。
【０１２８】
例えばＬtarget,iは、以下のように設定される。
【０１２９】
【数１３】
2009225550.tif 000015
ｎはシナリオの数である。
【０１３０】
本方法では、特に需要予測シナリオの高め側や低め側などのある特定順位のシナリオのみに着目し、これを例えば需要実績の上側１０％確率点、下側１０％確率点などと設定するような場合に適している。すなわち、危機管理などの観点から、予め定めた特定の確率点を遵守するように将来の需要予測を行う際に適している。例えば、重み係数は着目した順位の重みを１、それ以外は０とおく。なお、着目方法であるが上下一点ずつではなく複数点を指定してもよい。例えば上側５番目、１０番目の順位のシナリオをそれぞれ上側１０％、２０％点に、下側１０番目、５番目の順位のシナリオをそれぞれ、下側２０％、１０％点に指定するなどである。
【０１３１】
なお、上の例では目的関数は、発生確率や累積確率とその目標値との差の２乗値からなる関数であるが、２乗値に代えて絶対値を利用することも可能である。さらには制約条件を制約式として与える代わりに、以下に述べるような式のように、制約に対する逸脱量をペナルティー項として目的関数に加えるといった方法も可能である。
【０１３２】
【数１４】
2009225550.tif 000016
【０１３３】
ランダム法処理部５５５３は、需要予測設定入力部５７によって入力されたデータに従って、制約条件を設定する（ステップＳ１５３）。例えば、（１０）式のような制約条件を採用する。
【０１３４】
そして、ランダム法処理部５５５３は、需要予測設定入力部５７によって入力されたデータに従って、探索終了条件を設定する（ステップＳ１５５）。例えば、最小目的関数値Ｊmin及び探索上限回数Ｎmaxを設定する。また、カウンタＮを１に初期化する（ステップＳ１５７）。
【０１３５】
さらに、ランダム法処理部５５５３は、係数初期値を設定する（ステップＳ１５９）。係数初期値には、例えば予め定められた一定範囲内からランダムに与えるようにしても良い。また、一旦ＬＰ法などで係数値を算出して用いるようにしても良い。
【０１３６】
その後、ランダム法処理部５５５３は、図２３に示した前処理を実施する（ステップＳ１６１）。係数初期値を指定需要予測モデルに代入して計算を行う。さらに、前処理の結果である発生頻度分布テーブルを用いて、目的関数値Ｊ０を算出し、例えばメインメモリなどの記憶装置に格納する（ステップＳ１６３）。
【０１３７】
また、ランダム法処理部５５５３は、係数値の変更量をランダムに生成し、係数値に加算する（ステップＳ１６５）。例えば、正規分布の乱数を発生させて、元の係数値に加算する。例えばΔａ，Δｂk，Δｃk，Δｄkを生成し、ａ＋Δａ，ｂk＋Δｂk，ｃk＋Δｃk，ｄk＋Δｄkを算出する。そして、このように算出された係数値は、係数値に関する制約条件（例えば（１０）式）を満たしているか判断する（ステップＳ１６７）。制約条件を満たしていない場合には、係数値を元に戻し（ステップＳ１６９）、ステップＳ１６５に戻る。一方、制約条件を満たしている場合には、端子Ｄを介して図２７の処理に移行する。
【０１３８】
さらに、ランダム法処理部５５５３は、抽出過去気象予報データ格納部５５２に格納されているデータ並びに算出された係数値及び指定需要予測モデルを用いて需要予測値を算出し、例えばメインメモリなどの記憶装置に格納する（ステップＳ１７０）。そして、算出された需要予測値を用いて、前処理を実施する（ステップＳ１７１）。さらに、前処理によって生成された発生頻度分布テーブルを用いて、目的関数値Ｊ１を算出し、例えばメインメモリなどの記憶装置に格納する（ステップＳ１７３）。そして、Ｊ１＜Ｊ０であるか判断する（ステップＳ１７５）。Ｊ１＜Ｊ０は、目的関数の値が改善しているか判断するものである。Ｊ１＜Ｊ０でない場合には端子Ｅを介してステップＳ１６９を経由し係数値を１ステップ前に戻した上でステップＳ１６５に戻る。
【０１３９】
一方Ｊ１＜Ｊ０である場合には、ランダム法処理部５５５３は、Ｊ０＝Ｊ１として設定し、さらにＮを１インクリメントする（ステップＳ１７７）。その後、探索終了条件を満たしたか判断する（ステップＳ１７９）。すなわちＪがＪminに達したか、ＮがＮmaxに達したか判断し、両方達していない場合には、ステップＳ１６５で生成した変更量を、さらに係数値に加算して（ステップＳ１８１）、ステップＳ１７０に戻る。
【０１４０】
一方、ＪがＪminに達したか、又はＮがＮmaxに達した場合には、ランダム法処理部５５５３は、ステップＳ１６５で算出された係数値を、需要予測モデルデータ格納部５９に格納する（ステップＳ１８３）。
【０１４１】
以上のように処理すれば需要予測モデルが目的に応じた形で生成されるようになる。
【０１４２】
次に、図２８を用いて、ステップＳ１３で表示される画面の一例を説明する。図２８では、ＬＰ法もランダム法も両方とも実施した場合の画面例を示す。図２８の例では、需要予測対象日と、確率情報と、個別予測結果（ＬＰ法）と、個別予測結果（ランダム法）とが表示されるようになっている。確率情報には、ＬＰ法を使用した場合における需要上限値と需要下限値とが含む。また、ランダム法を使用した場合における上側１０％点の需要予測値と下側１０％の需要予測値とを含む。さらに、個別予測結果については、降順にソートされた需要予測値が列挙されている。
【０１４３】
なお、ランダム法だけ又はＬＰ法だけの場合には、該当するデータのみを表示するようにしても良い。さらに、個別予測結果については、全てのシナリオについて表示する必要はない。
【０１４４】
さらに、連続的に需要予測対象日時を設定して、連続して需要予測値を算出するようにしてもよい。その場合には、図２９に示すような需要予測値のカーブとして表示することも可能になる。
【０１４５】
以上本発明の実施の形態を説明したが、本発明はこれに限定されるものではない。例えば、ＬＰ法とランダム法を用いる例を示したが、他の方法を採用するようにしても良い。また、目的関数や制約条件については、その目的に応じて変更することも可能である。
【０１４６】
さらに、図１及び図３で示した機能ブロック図は一例であって、必ずしも実際のプログラムモジュール構成と一致しない場合がある。さらに図１では、電力需要予測装置５は、ネットワーク１及び７に接続されていることを前提としているが、ネットワークに接続せずに、ＣＤ−ＲＯＭなどの記憶媒体経由でデータを取得するようにしても良い。
【０１４７】
また、本発明は電力需要の予測だけではなく、風力発電量の予測や太陽光発電量の予測などの分散電源の発電出力予測などにも応用可能である。すなわち、電力需要に関するデータの代わりに、発電量のデータ（発電実績値や発電予測値）を用いて上で述べたような処理を実施すればよい。また電力取引などにおいて需要と電力価格が相関を示す際には価格予測への応用が可能である。
【０１４８】
なお、需要予測装置５、気象予報サーバ３などはコンピュータ装置であって、図３０に示すように当該コンピュータ装置においては、メモリ２５０１（記憶部）とＣＰＵ２５０３（処理部）とハードディスク・ドライブ（ＨＤＤ）２５０５と表示装置２５０９に接続される表示制御部２５０７とリムーバブル・ディスク２５１１用のドライブ装置２５１３と入力装置２５１５とネットワークに接続するための通信制御部２５１７とがバス２５１９で接続されている。オペレーティング・システム（ＯＳ）及びＷｅｂブラウザを含むアプリケーション・プログラムは、ＨＤＤ２５０５に格納されており、ＣＰＵ２５０３により実行される際にはＨＤＤ２５０５からメモリ２５０１に読み出される。必要に応じてＣＰＵ２５０３は、表示制御部２５０７、通信制御部２５１７、ドライブ装置２５１３を制御して、必要な動作を行わせる。また、処理途中のデータについては、メモリ２５０１に格納され、必要があればＨＤＤ２５０５に格納される。このようなコンピュータは、上で述べたＣＰＵ２５０３、メモリ２５０１などのハードウエアとＯＳ及び必要なアプリケーション・プログラムとが有機的に協働することにより、上で述べたような各種機能を実現する。なお、コンピュータ装置は複数台のコンピュータによって構成される場合もある。
【図面の簡単な説明】
【０１４９】
【図１】本発明の実施の形態に係るシステムの概要図である。
【図２】数値気象予報データを説明するための概念図である。
【図３】需要予測モデル生成部の機能ブロック図である。
【図４】本発明の実施の形態におけるメインフローを示す図である。
【図５】電力需要実績データ格納部に格納されるデータの一例を示す図である。
【図６】数値気象予報データ格納部に格納されるデータの一例を示す図である。
【図７】複数気象予報シナリオに基づいた予測の一例を示す図である。
【図８】需要予測式のバイアス調整について説明するための図である。
【図９】需要予測式のバイアス調整について説明するための図である。
【図１０】需要予測式を説明するための図である。
【図１１】需要予測式を説明するための図である。
【図１２】需要予測モデル生成処理の処理フローを示す図である。
【図１３】過去気象予報データ抽出処理の処理フローを示す図である。
【図１４】予報間隔の算出を説明するための図である。
【図１５】過去電力需要実績データ抽出処理の処理フローを示す図である。
【図１６】抽出すべき需要実績値の日時の算出を説明するための図である。
【図１７】気象予報データ加工処理の処理フローを示す図である。
【図１８】中央値算出処理の処理フローを示す図である。
【図１９】偏差算出処理の処理フローを示す図である。
【図２０】需要予測モデル生成処理の処理フローを示す図である。
【図２１】ＬＰ法による需要予測モデル係数推定処理の処理フローを示す図である。
【図２２】ＬＰ法による需要予測モデル係数推定処理の処理フローを示す図である。
【図２３】前処理の処理フローを示す図である。
【図２４】発生頻度分布テーブルの生成について説明するための図である。
【図２５】発生頻度分布テーブルの一例を示す図である。
【図２６】ランダム法による需要予測モデル係数指定処理の処理フローを示す図である。
【図２７】ランダム法による需要予測モデル係数指定処理の処理フローを示す図である。
【図２８】需要予測結果の表示例を示す図である。
【図２９】需要予測結果の他の表示例を示す図である。
【図３０】コンピュータの機能ブロック図である。
【符号の説明】
【０１５０】
１，７  ネットワーク    ３  気象予報サーバ
５  電力需要予測装置    ７  電力計測データ収集装置
９  電力計測データ収集装置    １１  電力データ計測装置
３１  数値気象予報生成部    ３３  数値気象予報配信部
５１  インターフェース部    ５３  電力需要実績データ格納部
５４  数値気象予報データ格納部
５５  需要予測モデル生成部    ５７  需要予測設定入力部
５９  需要予測モデルデータ格納部    ６１  最新数値気象予報データ格納部
６３  需要予測シナリオ計算部    ６５  需要予測シナリオデータ格納部
６７  需要予測シナリオ出力部
５１１  電力需要実績データ受信部    ５１３  数値気象予報データ受信部
【図１】
2009225550.tif 000017
【図２】
2009225550.tif 000018
【図３】
2009225550.tif 000019
【図４】
2009225550.tif 000020
【図５】
2009225550.tif 000021
【図６】
2009225550.tif 000022
【図７】
2009225550.tif 000023
【図８】
2009225550.tif 000024
【図９】
2009225550.tif 000025
【図１０】
2009225550.tif 000026
【図１１】
2009225550.tif 000027
【図１２】
2009225550.tif 000028
【図１３】
2009225550.tif 000029
【図１４】
2009225550.tif 000030
【図１５】
2009225550.tif 000031
【図１６】
2009225550.tif 000032
【図１７】
2009225550.tif 000033
【図１８】
2009225550.tif 000034
【図１９】
2009225550.tif 000035
【図２０】
2009225550.tif 000036
【図２１】
2009225550.tif 000037
【図２２】
2009225550.tif 000038
【図２３】
2009225550.tif 000039
【図２４】
2009225550.tif 000040
【図２５】
2009225550.tif 000041
【図２６】
2009225550.tif 000042
【図２７】
2009225550.tif 000043
【図２８】
2009225550.tif 000044
【図２９】
2009225550.tif 000045
【図３０】
2009225550.tif 000046