JP2015139320A - 配電系統の事故復旧方法、並びに配電系統の実負荷とその幅を推定する装置、方法、プログラム - Google Patents

Abstract

【課題】区間ごとの分散電源の出力と負荷が消費する有効電力を推定装置が正確に把握し、推定装置がその状態を考慮して区分開閉器だけを制御して事故復旧を行えるようにする。
【解決手段】入力されたデータの入力処理を行うとともに入力処理したデータを蓄積する(S1)。次いで、事故区間判定部により事故区間の判定を行う(S2)。そして事故直前における分散電源の出力と、負荷が消費する電力を推定する(S3)。本ステップの推定は、散布図、回帰式、若しくは回帰式と独立成分分析のいずれかを用いた推定を行う。ついで上記で推定した分散電源の出力と負荷が消費する電力についての推定値の幅を算出する(S4)。そのうえで事故復旧処理を行う(S5)。事故復旧処理にあたっては、配電系統の最大許容容量を限度に、推定装置が区分開閉器だけを制御して負荷の消費電力推定値分のみ事故復旧する。
【選択図】図１０

Description

本発明は、分散電源が連系した電力系統において、系統における事故時に系統を復旧させるための事故復旧方法、並びに、該事故復旧を確実に行うための配電系統における実負荷とその幅を推定する装置、方法、プログラムに関する。

近年、配電系統に、多くの分散電源（例えば、太陽光発電装置(PV：Photo Voltaic generation)，風力発電装置，燃料電池など）が導入されている。配電系統の適切な運用、特に事故復旧時の運用のためには、区間ごとの分散電源（例．太陽光発電装置）の出力と負荷が消費する有効電力を正確に把握し、その状態を考慮して事故復旧する必要がある。本明細書では、分散電源としてその代表例である太陽光発電装置(PV)を用いる場合について専ら説明するが、太陽光発電装置(PV)のみに限定されるものではない。なお、“区間”は、配電系統に事故が起こった場合にその周辺から電力を融通して事故復旧を行いやすくするために配電系統の要所に設けられる区分開閉器で区切られている電力接続域を意味するものである。

従来の配電系統では、太陽光発電装置(PV)が含まれていなかったため、
計測した電力＝負荷の消費する電力
の関係がなりたっていた。

そのため、従来の配電系統における事故復旧技術は、事故直前に計測した電力をまかなうだけの電力を融通するもの、つまり計測した電力をその周囲から融通して事故復旧させるものであった。

しかし、分散電源（例．太陽光発電装置）が大量導入された配電系統では、このような技術が使えない。例えば、負荷の消費する電力が１００ｋW、分散電源の出力量が１００ｋWの場合、計測される電力は０Wになる。事故がおきると分散電源は系統から解列し出力が０ｋWになる。事故復旧に必要な電力は、１００ｋWになるが、事故直前に計測した電力は０Wであるので、従来技術で融通しようとすると融通電力が不足することになる。

分散電源を考慮して事故復旧を図る従来の技術を開示した文献の幾つかを以下に例示する。すなわち、
下記特許文献１には、区分開閉器と分散電源を監視制御する配電制御装置を設けておき、配電制御装置が区分開閉器を操作するだけでなく、分散電源のオン／オフを制御して、負荷の合計が配電線の許容容量を超えないように制御する事故復旧方法が開示されている。

また下記特許文献２には、供給支障区間が発生しないように、区間内の分散電源の発電出力量を増加させて、各分散電源の総発電出力量と各需要設備の総負荷量とを一致させる事故復旧方法が開示されている。

また下記特許文献３には、負荷と分散電源とが混在する電力系統の事故復旧操作を行う配電系統監視制御装置において、配電系統監視制御装置があらかじめ選択した負荷を切り離しておき、負荷を十分に軽くした状態で配電線遮断器を投入して事故復旧させる方法が開示されている。

その一方、事故復旧に当たってはその前提として、事故直前における、実際に負荷が消費する電力と分散電源(例．太陽光発電装置)の出力量を正確に推定する技術が必要になる。そのため従来では、分散電源(例．太陽光発電装置)の出力量と負荷の消費する電力を推定する技術として、以下のような技術が知られている。すなわち、
下記特許文献４には、太陽光発電システムが導入される前など太陽光発電出力がないデータを用いて負荷が消費する電力と太陽光発電装置の出力量を推定する推定式を構築する技術が開示されている。

また下記特許文献５には、計測した有効電力と無効電力から、太陽光発電装置の出力する有効電力と、負荷が消費する有効電力とを独立成分分析を用いて分離する技術が開示されている。

特開２００７−０２８７６９号公報 特開２０１１−０６１９３１号公報 特開２００６−０６０８８５号公報 特開２０１２−１９１７７７号公報 特開２０１２−０９５４７８号公報

上記特許文献１ないし３は、いずれも、中央の制御装置(例．配電制御装置、配電系統監視制御装置)から分散電源（例．太陽光発電装置）、もしくは、負荷が消費する電力を直接制御するようにしている。しかし、実際の配電系統で分散電源（例．太陽光発電装置）や負荷が消費する電力を直接制御できることはまれであり、ほとんどの場合は分散電源の発電電力も負荷が消費する電力も共に制御できない。中央の制御装置が直接制御できるのは区分開閉器だけである。

また上記特許文献４は、太陽光発電出力がないデータがないと推定に適用できないという問題がある。太陽光発電出力がないのは雨天時などがあるが、日本の場合は、雨天が他の天気の日と比較すると少ないので、精度の高い推定式が構築できないという課題がある。

また上記特許文献５は、太陽光発電装置の出力する有効電力と、負荷が消費する有効電力を推定することはできるが、推定値の幅、つまり、どのぐらいその推定値がはずれそうなのかを提示することができないため、事故復旧の際に必要な電力を見誤るという課題がある。

実際の負荷の推定値が誤差なく推定できれば問題ないが、現状の技術レベルでは誤差無し推定は不可能である。よって、推定した値にある程度の余裕分をみて事故復旧時の融通電力とする必要がある。例えば、復旧に必要な負荷の推定値を９０ｋWと推定しても、余裕分をみて、１００ｋWや１１０ｋWの電力を融通する必要がある。しかしながら、この方法では推定した負荷と実際の負荷との間にどの程度の誤差があるのかが事前に分からないため、事故復旧の際の余裕分をどのくらい見積もればよいか正確に分からないという問題がある。

そこで本発明の課題の一つは、区間ごとの分散電源の出力と負荷が消費する有効電力を推定装置が正確に把握し、推定装置がその状態を考慮して区分開閉器だけを制御して事故復旧を行えるようにすることである。

また本発明の課題のもう一つは、分散電源である太陽光発電出力がある日でもない日でも、区間ごとに計測したデータを、太陽光発電装置の出力と負荷が消費する有効電力を推定するための回帰式、もしくは独立成分分析と回帰式を用いて、配電系統の実負荷及びその幅を正確に推定して事故復旧時の電力融通に利用可能とすることである。

上記課題を解決するために請求項１に記載の発明は、複数の分散電源と負荷が配電線に接続された電力系統の事故復旧方法において、
複数の区分開閉器を設けて前記配電線を複数の区間に区分するステップ、
前記配電線上の前記区間ごとの負荷を計測して推定手段に設けた記憶装置に記憶すると共に配電系統の最大許容容量を予め該記憶装置に記憶するステップ、
前記推定手段に設けた記憶装置に記憶された前記計測負荷を基に前記区間ごとの事故直前の分散電源出力と負荷の消費電力を推定するステップ、
事故後、前記配電系統の事故区間を特定するステップ、および、
前記配電系統の事故区間を特定した後、事故区間の負荷の消費電力推定値を用い、前記配電系統の最大許容容量を超えない範囲で、前記推定手段が前記区分開閉器を操作する事故復旧するステップ、
を含むことを特徴とする。

前記請求項１記載の発明において、前記推定手段が推定する前記負荷の消費電力推定値には、前記事故復旧時に過負荷とならない推定値の幅（余裕率）を持たせ、前記事故復旧するステップは、前記推定値の幅（余裕率）を含んだ前記負荷の消費電力推定値のみを事故復旧することを特徴とする（請求項２記載の発明）。

また前記請求項１に記載の発明において、前記区間ごとの事故直前の分散電源出力と負荷の消費電力を推定するステップは、前記記憶装置に記憶した前記区間ごとの前記負荷装置の消費電力、および該当時の気温を座標軸にプロットして得た散布図を用いて、前記負荷装置の消費電力と前記分散電源の出力を推定することを特徴とする（請求項３記載の発明）。

また前記請求項１ないし３のいずれか一項に記載の発明において、前記事故区間内の各区間の推定消費電力から前記分散電源の推定出力を除外した値を事故復旧に用いることを特徴とする（請求項４記載の発明）。

また前記請求項１ないし３のいずれか一項に記載の発明において、高圧側に連系される高圧系分散電源の出力を計測する手段を該高圧系分散電源装置が備える場合には、
前記事故区間の各区間の推定消費電力から分散電源の推定出力を減算し、その値に前記高圧系分散型電源の出力の計測値を加算した値を事故復旧に用いることを特徴とする（請求項５記載の発明）。

前記請求項４または５に記載の発明において、前記推定手段が推定する前記負荷の消費電力推定値には、前記事故復旧時に過負荷とならない推定値の幅（余裕率）を持たせ、前記負荷の消費電力推定値分のみ事故復旧するステップは、前記推定値の幅（余裕率）を含んだ前記負荷の消費電力推定値のみを事故復旧することを特徴とする（請求項６記載の発明）。

上記課題を解決するために請求項７に記載の発明は、複数の分散電源及び複数の負荷装置が配電線に接続された電力系統において、
複数の区分開閉器を設けて前記配電線を複数の区間に区分する手段と、
前記配電線上で前記区間ごとの有効電力および無効電力を計測して推定手段に設けた記憶装置に記憶する手段と、
前記記憶装置に記憶された前記区間ごとの前記分散電源および前記負荷装置の無効電力から前記区間ごとの負荷装置の有効電力を回帰式で推定する手段と、
前記回帰式で推定した前記区間ごとの負荷装置の有効電力推定値に対してその推定値の幅を算出する手段と、
を備えることを特徴とする。

前記請求項７に記載の発明において、前記回帰式で求められた区間ごとの負荷装置の有効電力推定値を基に、前記分散電源の有効電力を推定する手段と、
前記回帰式で推定した前記区間ごとの分散電源の有効電力推定値に対してその推定値の幅を算出する手段と、をさらに備えることを特徴とする（請求項８記載の発明）。

上記課題を解決するために請求項９に記載の発明は、複数の分散電源及び複数の負荷装置が配電線に接続された電力系統において、
複数の区分開閉器を設けて前記配電線を複数の区間に区分するステップ、
前記配電線上で前記区間ごとの有効電力および無効電力を計測して推定手段に設けた記憶装置に記憶するステップ、
前記記憶装置に記憶された前記区間ごとの前記分散電源および前記負荷装置の無効電力から前記区間ごとの負荷装置の有効電力を回帰式で推定するステップ、および、
前記回帰式で推定した前記区間ごとの負荷装置の有効電力推定値に対してその推定値の幅を算出するステップ、
を含むことを特徴とする。

前記請求項９に記載の発明において、前記回帰式で求められた区間ごとの負荷装置の有効電力推定値を基に、前記分散電源の有効電力を推定するステップ、
前記回帰式で推定した前記区間ごとの分散電源の有効電力推定値に対してその推定値の幅を算出するステップ、を含むことを特徴とする（請求項１０記載の発明）。

また前記請求項９または１０に記載の発明において、前記計測して有効電力と無効電力から独立成分分析を用いて、前記負荷装置の有効電力、もしくは、前記分散電源の出力を分離するステップと、
該分離した前記負荷装置の有効電力、もしくは、前記分散電源の有効電力と、前記計測した無効電力とから、前記回帰式を用いて、前記負荷装置の有効電力、もしくは、前記分散電源の有効電力を推定するステップと、
該推定した前記負荷装置の有効電力、もしくは、前記分散電源の有効電力の推定値の幅を算出するステップと、さらに含むことを特徴とする（請求項１１記載の発明）。

上記課題を解決するために請求項１２に記載の発明は、複数の分散電源及び複数の負荷装置が配電線に接続された電力系統において、
該電力系統に配電系統の実負荷を推定する推定装置を備え、該推定装置のコンピュータを、
複数の区分開閉器を設けて前記配電線を複数の区間に区分する手段と、
前記配電線上で前記区間ごとの有効電力および無効電力を計測して記憶装置に記憶する手段と、
前記記憶装置に記憶された前記区間ごとの前記分散電源および前記負荷装置の無効電力から前記区間ごとの負荷装置の有効電力を回帰式で推定する手段と、
前記回帰式で推定した前記区間ごとの負荷装置の有効電力推定値に対してその推定値の幅を算出する手段、として機能させるためのプログラムとしたものである。

上記課題を解決するために請求項１３に記載の発明は、複数の分散電源及び複数の負荷装置が配電線に接続された電力系統において、
該電力系統に配電系統の実負荷を推定する推定装置を備え、該推定装置のコンピュータを、
複数の区分開閉器を設けて前記配電線を複数の区間に区分する手段と、
前記配電線上で前記区間ごとの有効電力および無効電力を計測して記憶装置に記憶する手段と、
前記記憶装置に記憶された前記区間ごとの前記分散電源および前記負荷装置の無効電力から前記区間ごとの負荷装置の有効電力を回帰式で推定する手段と、
前記回帰式で推定した前記区間ごとの負荷装置の有効電力推定値に対してその推定値の幅を算出する手段と、
前記回帰式で求められた区間ごとの負荷装置の有効電力推定値を基に、前記分散電源の有効電力を推定する手段と、
前記回帰式で推定した前記区間ごとの分散電源の有効電力推定値に対してその推定値の幅を算出する手段、として機能させるためのプログラムとしたものである。

本発明によれば、分散電源が連系された電力系統において、配電系統に接続された分散電源や負荷の全ての出力を正確に計測しなくても、事故時に融通すべき負荷が消費する電力を推定でき、過負荷を起こすことなく事故復旧することが可能となる。また、負荷や分散電源を個別に制御することなく区分開閉器を制御するだけで事故復旧を図ることが可能である。

また本発明によれば、太陽光発電装置に代表される分散電源の出力量を個別に計測する計測装置を設置することなく、配電系統全体の太陽光出力の有効電力と、負荷が消費する有効電力とを分離して推定することができ、また、その推定値の幅を求めることができる。一方、上記特許文献４及び５に示された技術は、いずれも推定値の幅を求めることができないため、推定値に基づいて事故復旧する場合に、その推定誤差をどの程度みつもればよいか判断できないが、本発明によれば、推定値の推定誤差を定量的に求められるので適切な事故復旧を行うことが可能になる。

本発明の実施形態に係る配電系統の構成例を示す図である。 本発明の実施形態に係る電力推定装置(その１)の構成を示す図である。 本発明の実施形態にかかる計測負荷と雨天時の気温との関係を示す散布図(その１)である。 本発明の実施形態にかかる、多くのデータから作成された計測負荷と雨天時の気温との関係を示す散布図(その２)である。 本発明の実施形態に係る配電系統の構成例における配電系統の事故前の状態を示す図である。 本発明の実施形態に係る配電系統の構成例における配電系統の事故発生直後の状態を示す図である。 本発明の実施形態に係る配電系統の構成例における配電系統の事故復旧で想定される事故復旧の様子を示す図である。 本発明の実施形態に係る配電系統の構成例における配電系統の事故復旧の様子を示す図である。 本発明の実施形態に係る配電系統の構成例における配電系統が高圧連系と低圧連系とで構成されていた場合の事故復旧の様子を示す図である。 図８に示した配電系統の構成の一例を示す図である。 本発明の実施形態に係る電力推定装置(その１)の動作を説明する処理フローを示す図である。 本発明の実施形態に係る電力推定装置(その２)の構成を示す図である。 本発明の実施形態に係る電力推定装置(その２)の回帰式構築処理部が回帰式（式２）を得るための有効電力及び無効電力の関係を示す図である。 本発明の実施形態に係る電力推定装置(その２)における実負荷の有効電力推定値（回帰）が採る推定値の幅を示す概念図である。 本発明の実施形態に係る電力推定装置(その３)の構成を示す図である。 本発明の実施形態に係る電力推定装置(その３)の回帰式構築処理部が回帰式（式３）を得るための有効電力及び無効電力の関係を示す図である。 本発明の実施形態に係る電力推定装置(その３)における実負荷の有効電力推定値が採る推定値の幅を示す概念図である。

以下、本発明の実施の形態について、図面を参照しながら説明する。
図１は、本発明の実施形態に係る配電系統の構成例を示す図である。図１を用いて、本発明の実施形態に係る配電系統について説明する。

図１に示すように、変電所１から下流に位置するＡ系統及びＢ系統に分岐された配電線に、負荷や分散電源の代表例である太陽光発電装置(PV：Photo Voltaic generation )がつながる。そして配電線には、１もしくは複数の電力計（図示例では、電力計P1,・・・,P4）が設けられ、電力計で計測した有効電力と無効電力は定期的に電力推定装置１０に伝送される。なお電力計P1,・・・,P4では、上記した有効電力、無効電力のほか、電流、力率、開閉器の開閉状態等についても計測され、電力推定装置１０に伝送される。

図１について更に詳しく説明すると、変電所１に接続される配電線のＡ系統（フィーダともいう）には、一例として２棟の一般家庭と１つ工場が接続され、配電線のＢ系統（フィーダともいう）には、一例として３棟の一般家庭が接続されている。そして一般家庭の屋根等には、太陽光発電パネルが設置され、該太陽光発電パネルにより発電された電力はインバータを介して交流に変換されるとともに、家庭内の負荷で電力を消費したり、太陽光発電電力が余れば系統に電力を供給することが出来るよう構成している。なお、変電所１には発電所（不図示）から送電線を介して送電が行われている。

図１に示す配電系統の構成例は、単なる例示であり、実際は、これより多くの一般家庭、工場、メガソーラ等が配電線に接続されることになる。また図示例の工場は、太陽光発電装置(発電パネル)を有さずに、工場内のモータ、コンデンサ（リアクタンス）、及び、一般的な電力負荷、等から成るエネルギーを消費する負荷が接続され、配電線から供給される電力を消費するものとして記述されている。なお、無効電力補償装置200は、無効電力の調整を行うもので、それ自体の機能は当業者に知られているのでその説明を省略する。

また図示例の配電線には、開閉器S1の変電所１側に電力計P1が設けられ、電力計P1より下流の配電線に接続される負荷全てが消費する電力値が、また開閉器S2の変電所１側に電力計P2が設けられ、電力計P2より下流の配電線に接続される負荷全てが消費する電力値が、また開閉器S3の変電所１側に電力計P3が設けられ、電力計P3より下流の配電線に接続される負荷全てが消費する電力値が、さらに開閉器S4の変電所１側に電力計P4が設けられ、電力計P4より下流の配電線に接続される負荷全てが消費する電力値が、それぞれ計測され、電力計P1，電力計P2、および、電力計P3，電力計P4でそれぞれ計測した有効電力と無効電力の値が定期的に電力推定装置(その１)１０に伝送される。なお、図１では、S1〜S4を単に“開閉器”と表記しているが、正式には“区分開閉器”である。しかし図示上で図が煩雑になるため敢えて“開閉器”と表記することにしたものである。また、S5，S6は、系統を跨いで電力を融通するときに使用されるものであるので、“開閉器”として表記するものである。

このようにすることで各(区分)開閉器によって区切られた配電線の各区間ごとに有効電力と無効電力の値を計測することが可能となる。また配電線の各区間ごとに有効電力と無効電力の値を計測することによって、後述するように、散布図、回帰式、若しくは、回帰式と独立成分分析を用いて、区間ごとの実負荷の有効電力値や太陽光発電装置の有効電力値を推定することが可能となる。

図２は、本発明の実施形態に係る電力推定装置(その１)の構成を示す図である。図２に示す本発明の電力推定装置(その１)１０は、コンピュータ等の情報を処理する汎用の装置で構成されており、当該装置には、特に図示しないが、ＣＰＵ（中央処理ユニット）、記憶装置（ハードディスク）、メモリ、通信機能部、入出力インタフェース、入出力装置など当業者によく知られたハードウェア構成を備えている。そして上記記憶装置には予め所定のアプリケーションプログラムが記憶されており、上記ＣＰＵがこのアプリケーションプログラムを読出して実行することにより下記に記述する各種処理部１１〜１７の機能を実現する。

上記電力推定装置(その１)１０が備えるハードウェア構成の下に、
入力処理部１１は、配電系統に備えられた電力計P1，電力計P2、および、電力計P3，電力計P4で計測された有効電力、無効電力、電流、力率、開閉器の開閉状態、などの各値（データ）を取り込む。図示していないが、外部(例えば、気象庁等)により提供される、気象情報・暦情報など計測器以外のデータも取り込めるようにしており、また、必要に応じて、担当者(図示せず)が入力処理部１１から上記情報を直接入力することもできる。

データ蓄積処理部１２は、入力処理部１１で取り込んだ情報を蓄積するとともに、後述する電力推定装置(その１)１０内の各機能部１３〜１７で処理した情報を適宜蓄積する。
出力処理部１３は、データ蓄積処理部１２に蓄積された情報を画面に表示したり、紙等に印刷して出力する。また出力処理部１３は、電力推定装置(その１)１０内の各機能部１１〜１７で処理し蓄積した情報をネットワーク等を介して外部の装置(図示せず)に出力することができる。

上記した入力処理部１１及び出力処理部１３は、上述した情報処理装置の入出力装置により実現される。またデータ蓄積処理部１２は上述した情報処理装置の記憶装置（ハードディスク）、メモリ等により実現される。バス１８は上記処理部及び後記する処理部を接続するものである。

一方、電力推定装置(その１)１０の事故区間判定部１４は、事故が起きた区間を判定する処理を行う。本処理は、従来から我が国で実施されている通常の方法を用いる。すなわち我が国では、上流側(変電所１側)の開閉器を順次操作することで判定する方法が一般的である。つまり、事故発生時には事故が発生した区間を含むフィーダのすべての開閉器がいったん開く。その後、上流側の開閉器を閉じて過電流が流れなければ、さらに１つ下流の開閉器を閉じる。この操作を自動で継続して過電流が発生するまで行い、もし過電流が発生したならば当該開閉器より下流の区間を事故区間と判定するものである。

推定処理部１５は、分散電源（例．太陽光発電装置）の出力と、負荷が消費する電力を区間ごとに推定する処理を行う。推定方法としては、様々な方法があるが、本処理部では、以下の方法を用いる。もちろん本発明は、本例のみに限定されるものではない。すなわち、
図３に示されるように、計測された負荷（計測負荷）と雨天時の気温との関係を示す散布図(その１)を作成する。図３の散布図(その１)に用いられるデータとしては、例えば過去１か月の雨天時のデータが使用され、これを計測負荷及び温度から成る座標軸上にプロットすることで得られる。散布図は、例えば一日のうちの１時〜２４時まで１時間ごと、且つ平日及び休日ごとに１部ずつ作成する。上記では１時間ごとの刻みで作成する例を示すが、３０分ごとや１５分ごとの刻みで作成してもよい。なお図３ではプロットの数を割愛して示している。

ここで図３の散布図(その１)を参照し、該当時の気温と雨天時の当該散布図が示すプロットで交わるA点を負荷の消費電力として推定する。図３のB点に示す値は、事故発生当日の事故直前に上記した電力計で計測された計測値である。したがって、事故発生直前の分散電源（例．太陽光発電装置）の出力推定値は、A-Bとなる。この場合、気温が一致するA点がなければ、当該気温に一番近い雨天時の散布図上のデータを代用すればよい。そして上記計測値が、図１の電力計P1で計測されていれば、電力計P1以下の負荷が消費する電力と分散電源（例．太陽光発電装置）の出力による推定値となる。電力計P2で計測されていれば同様に、電力計P2以下の負荷が消費する電力と分散電源（例．太陽光発電装置）の出力による推定値となる。電力計の設置数が少ない場合には、区分ごとの負荷が消費する電力は正確には分からない。例えば、電力計P2がなければ電力計P1と電力計P2の間の区間の負荷が消費する電力は分からない。その場合は、各区間の契約容量から按分して推定することになる。

次に、図４を用いて推定幅処理部１６の処理について説明する。図４は、図３とは異なる散布図(その２)を示すものである。図３に示すような雨天時の散布図を多くのデータから作成すると、図４に示すように同じ気温であっても計測負荷に幅（両矢印線参照）が生じる。そこで、事故直前の気温において、図４に示すA点の幅をもとめ、その半分の値と、事前にデータベース(図示はしていないが、図２のデータ蓄積処理部１２内に作成される)上に設定される経験則で得られた任意の余裕率との加算により、以下に示される推定幅（余裕率）が算出される。すなわち、
推定値の幅（余裕率）＝ 散布図上のA点が示す幅 / 2 ＋ 任意の余裕率・・・式１
次いで、図２に示した事故復旧処理部１７の処理について、図５Ａ〜図８を用いて説明する。なお、図５Ａ〜図８に示す例では、任意の数の区間（上記した開閉器（図示例では矩形表示）によって区切られた区間）を設け、さらに、説明を簡略化するためにどの区間も、負荷が500kW、分散電源として太陽光発電装置（PV）が200kW出力しているものとする。つまり、各区間の見かけ上の負荷は300kWとなるものとする。そして図５Ａ〜図８において、白抜き矩形は開閉器が“開”の状態を示し、また黒矩形は開閉器が“閉”の状態を示すものとしている。

図５Ａの例では、各フィーダ（図１に示すＡ系統，Ｂ系統）に６区間が設けられる例を示しているので、変電所１から各フィーダにそれぞれ1800kWずつ送り出している。なお、変電所１から各フィーダに供給できる最大許容容量は4000ｋWとする。そして図５Ｂに示すような事故発生直後には、PV出力が無くなるため、もしも図６に示すように、正常なフィーダ(図示例ではＢ系統)から事故区間(吹き出し線表示)以外すべてを復旧させようとすると、合計の負荷が4300kWとなり、フィーダの最大許容容量4000kWを超えてしまうことになる。

そこで本発明の実施形態に係る事故復旧法に基づく事故復旧では、4000kWを超えない範囲で復旧させるので、図７に示すとおり、Ａ系統の事故区間及びその１つ下流の区間を除いた、合計3800kW分について復旧させるようにしている。なお、各区間の負荷については、上記した式１に示した「任意の余裕率」を含んだ、「負荷の消費電力推定値＋推定値の幅（余裕率）」に基づいて推定した値を使用している。

また本発明の実施形態では、上記に示す事故復旧法に加えて、図１に示す配電系統が図８に示されるような高圧連系の分散電源と、低圧連系の太陽光発電装置（PV）とで構成されている場合には、以下に示す事故復旧法に基づいて事故復旧を行う。すなわち、
分散電源が低圧に連系している太陽光発電装置（PV）の場合は、復旧後短時間で再連系されることが知られている。つまり短期的な過負荷状態では配電線は何らの問題も生じないので、低圧系の太陽光発電装置（PV）分については解列しないとみなしてよい。もし、分散電源がすべて低圧系の太陽光発電装置（PV）であるならば、停電区間の各区間の200kW分を負荷から除外できるので、各区間が300kWとみなして事故復旧処理を行ってもよい。

一方、高圧系に連系される分散電源が有る場合には、高圧系に連系される分散電源は、自動的に再連系することはなく、事故復旧後必ず手動で連系させるようになっている。よって高圧系に連系している分散電源は、事故時に解列されたとみなして事故復旧の処理を行うようにする。なお、高圧系に連系している分散電源は、通常、負荷が消費する電力および分散電源の出力が不図示の電力計で計測されているので、その分散電源の出力分については容易に把握できる。ここで高圧系と低圧系とは、通常、変電所１より給電される電圧が柱上変圧器を介さずに受電できるようにされている場合が高圧系とされ、一方、柱上変圧器を介して受電できるようにされている場合が低圧系とされている。

さらに具体例について図８を用いて説明する。図示例では、事故区間を含むＡ系統に高圧連系の分散電源が３区間あり、低圧連系の太陽光発電装置（PV）が２区間あるとする。高圧連系の分散電源は自動で再連系しないため、該当区間それぞれの融通すべき負荷は500ｋWである。低圧連系のPVは比較的短期間に自動的に連系するため、融通すべき負荷はPVの出力分減って300ｋWとなる。よって、融通すべき負荷電力は、事故区間を除く全負荷の1800ｋＷ＋600ｋＷ＋1500ｋＷ＝3900ｋWとなり、この分を正常なフィーダ（図示例ではＢ系統）を通じて復旧させることができる。つまり、この場合には、図６や図７に示す事故復旧処理では不可能であった、事故区間(吹き出し線表示)以外のすべての区間を復旧させることができる。

図９は、図８に示した高圧連系の分散電源と、低圧連系の太陽光発電装置（PV）とで配電線が構成される場合の具体的な配電系統の構成例を示す図である。図９では、各フィーダ、すなわちＡ系統，Ｂ系統のそれぞれにおいて、低圧連系の太陽光発電装置（PV）と、高圧連系の分散電源(例．メガソーラ)とが配電線に接続されるように構成されている。

以上における事故復旧動作を図１０の処理フローを用いて説明する。図１０は、本発明の実施形態に係る電力推定装置(その１)の動作を説明する処理フロー図である。図１０においてはステップを“S”と略記する。また図１、図２及び図９を適宜参照するものとする。

図１０において、まず、入力されたデータの入力処理を図２に示す入力処理部１１を用いて行うとともに入力処理したデータをデータ蓄積処理部１２に蓄積する（ステップS1）。次いで、ステップS2では、事故区間判定部１４により事故区間の判定を行う。

ステップS3では、推定処理部１５を用いて、事故直前における分散電源（例．太陽光発電装置）の出力と、負荷が消費する電力を推定する。本ステップでは、上述した散布図を用いた推定のほか、後述する回帰式あるいは、回帰式と独立成分分析とを用いた推定を行う。

次にステップS4では、推定幅処理部１６を用いて上記で推定した分散電源（例．太陽光発電装置）の出力と負荷が消費する電力についての推定値の幅を算出する。
そのうえでステップS5において、事故復旧処理を行う。事故復旧処理にあたっては、図７および図８で説明した事故復旧処理法が用いられる。

最後にステップS6において、図２に示す出力処理部１３を用いてデータの出力処理を行う。このようにして本発明の実施形態における一連の事故復旧処理が実行される。
以上説明したように、各区間に融通すべき負荷は、推定した負荷の消費電力から、推定した分散型電源の出力（高圧系および低圧系を含む）を減算し、その値に高圧系分散型電源の出力（計測値）を加算したものとして計算することができる。分散型電源の出力は高圧系分散型電源を含んだ値として推定されるが、高圧系分散型電源はすぐには復旧しないため、その出力値（測定値）をないものとして、自動で事故復旧する区間を見積もることが必要である。

図１１は、本発明の実施形態に係る電力推定装置(その２)の構成を示す図である。本電力推定装置(その２)２０は、図２に示した本発明の実施形態に係る電力推定装置(その１)１０の構成と同様、コンピュータ等の情報を処理する汎用の装置で構成されており、当該装置には、特に図示しないが、ＣＰＵ（中央処理ユニット）、記憶装置（ハードディスク）、メモリ、通信機能部、入出力インタフェース、入出力装置など当業者によく知られたハードウェア構成を備えている。そして上記記憶装置には予め所定のアプリケーションプログラムが記憶されており、上記ＣＰＵがこのアプリケーションプログラムを読出して実行することにより下記に記述する各種処理部２１〜２７の機能を実現する。なお、本発明の実施形態に係る配電系統の構成は図１に示した電力推定装置(その１)１０が電力推定装置(その２)２０に置き換わるのを除き、同様であるので、機能ブロック名が重複するものはその説明を省略するものとする。また上記各種処理部２１〜２７を接続するバスが設けられている。

本電力推定装置(その２)２０では、説明を簡単にするため１か所の電力計、例えば電力計P1、によって計測された値に基づく推定について説明する。この手法では計測地点から下流の負荷が消費する有効電力を推定するので、開閉器ごとに計測すれば、区間ごとの負荷が消費する有効電力を推定することができる。

電力推定装置(その２)２０のデータ抽出処理部２４は、後述する回帰式構築処理部２５にて必要となる回帰式構築用のデータを抽出する処理を行う。例えば、以下のような複数の抽出条件の組合せで抽出処理を行う。
抽出条件１：推定対象日から過去ｘ日以内の日(例．ｘ=30)
抽出条件２：t1時〜t2時(例．t1=10，t2=16)
抽出条件３：日射量時間0（つまり、PV発電量がない日）
ここで、日射量時間は、通常は気象庁にて計測している１時間当たりの日照時間であり、0〜10の値である。0のときには規定値以上の日射量がなく、規定値以上の日射量が6分ごとに１加算され、該当する１時間枠の中で60分あれば10となる。つまり、0は夜か雨天・曇天を示し、10は晴天を示すことになる。

なお、抽出したデータは、具体的には、以下のｎ個のデータになる。
有効電力計測値i ，無効電力計測値i
ただしi : 1〜ｎ
次に、回帰式構築処理部２５は、上記で抽出されたデータの無効電力計測値と、有効電力計測値を用いて以下の回帰式を構築する。回帰式の構築に当たっては、図１２に示すように、有効電力をＹ座標に、無効電力をＸ座標に設定し、当該座標上に上記データ抽出処理部２４によって抽出された有効電力と無効電力の値をプロットすることでその傾きから次のような回帰式を得る。つまり、式２のa0とa1を求める。この場合、データ抽出処理部２４が抽出条件３を選択することにすれば、日照時間が0のデータ、例えば雨天時のデータを得ることができ、これにより、計測データは負荷が消費する有効電力と無効電力を表すことになるから、負荷が消費する電力における有効電力と無効電力は、ほぼ直線的な関係としてグラフ表示されることが分かる。

有効電力計測値＝a0+a1無効電力計測値・・・式２
本例では、計測された有効電力と無効電力の値を用いているので、誤差が少ない（誤差があるとすれば計測誤差だけで、後述する独立成分分析による推定誤差の影響がなく）信頼性が高いデータである。その一方で、抽出条件が多く、少ない日のデータしか得られない。

上記において精度の高い回帰式を構築する（上記の係数a0,a1を求める）ためには、誤差が少なく多くのデータが必要であるが、実際にはそのようなデータは得ることができない。そこで、データの重きをいずれに置くかによって後述する独立成分分析による推定誤差を含む、つまり、データの多さに重きを置いて回帰式を構築するか、上記のように、少ないデータでも誤差の少なさに重きを置いて回帰式を構築するかの二者択一になる。

推定処理部２６は、回帰式を示す上記式２を用いて、抽出したデータの推定値を求める。すなわち、
実負荷の有効電力推定値（回帰）i＝a0+a1無効電力計測値i
太陽光発電の有効電力推定値i=実負荷の有効電力推定値（回帰）i
-有効電力計測値i
推定幅計算処理部２７は、上記した推定値の幅を求める。幅を求める式は以下の通りである。

実負荷の有効電力推定値（回帰）の幅i=実負荷の有効電力推定値（回帰）i±幅i
太陽光発電の有効電力推定値の幅i=実負荷の有効電力推定値（回帰）の幅i
-有効電力計測値i
幅i=t(ｎ−２，信頼区間) × √（Ve）×√（１＋1/n＋（無効電力計測値i-avg(無効電力計測値i)）²/Sxx）
Ve=Σ（実負荷の有効電力推定値（回帰）i-有効電力計測値i）²／(n-2)
Sxx=Σ（無効電力計測値i-avg(無効電力計測値i)）²
ここで、t(ｎ−２，信頼区間)は、ｔ分布表から求められる値であり、ｎが無限大の場合には正規分布になり、ｎが小さい場合には正規分布よりも幅が広くなる分布である。例えば、ｎが12で信頼区間が95％の場合には、2.228になり、ｎが無限大で信頼区間が95％の場合には1.960になる。なお、Veは、残差による分散であり、上式のとおり、実負荷の有効電力推定値（回帰）と有効電力計測値の差を意味する。Sxxは、偏差平方和と呼ばれる統計指標である。

図１３は、本発明の実施形態に係る電力推定装置(その２)における実負荷の有効電力推定値（回帰）が採る推定値の幅を示す概念図である。図１３において、○印は、有効電力計測値iを示し、△印は、実負荷の有効電力推定値（回帰）iを示し、実負荷の有効電力推定値（回帰）iは△印を結んだ直線上の値となる。また図１３において、左上及び右下の曲線は上記幅iを求める計算式によって既定される推定値の確からしさを担保する境界を示し、当該二つの曲線の乖離が両矢印線によって推定値の幅iとして示されている。

図１４は、本発明の実施形態に係る電力推定装置(その３)の構成を示す図である。本電力推定装置(その３)３０は、図２に示した本発明の実施形態に係る電力推定装置(その１)１０の構成と同様、コンピュータ等の情報を処理する汎用の装置で構成されており、当該装置には、特に図示しないが、ＣＰＵ（中央処理ユニット）、記憶装置（ハードディスク）、メモリ、通信機能部、入出力インタフェース、入出力装置など当業者によく知られたハードウェア構成を備えている。そして上記記憶装置には予め所定のアプリケーションプログラムが記憶されており、上記ＣＰＵがこのアプリケーションプログラムを読出して実行することにより下記に記述する各種処理部３１〜３９の機能を実現する。なお、本発明の実施形態に係る配電系統の構成は図１に示した電力推定装置(その１)1０が電力推定装置(その３)３０に置き換わるのを除き、同様であるので、機能ブロック名が重複するものはその説明を省略するものとする。また上記各種処理部３１〜３９を接続するバスが設けられている。

本電力推定装置(その３)３０では、説明を簡単にするため１か所の電力計、例えば電力計P1、によって計測された値に基づく推定について説明する。この手法では計測地点から下流の負荷が消費する有効電力を推定するので、開閉器ごとに計測すれば、区間ごとの負荷が消費する有効電力を推定することができる。

電力推定装置(その３)３０のデータ抽出処理部３４は、後述する復元行列構築処理部３５にて必要となる復元行列構築用のデータを抽出する処理を行う。例えば、以下のような複数の抽出条件の組合せでデータの抽出処理を行う。この抽出条件は図１１に示した電力推定装置(その２)２０のデータ抽出処理部２４で用いる抽出条件よりも少ないものになっている。
抽出条件１：推定対象日から過去ｘ日以内の日(例．ｘ=30)
抽出条件２：t1時〜t2時(例．t1=10，t2=16)
因みに通常の独立成分分析（ICA：Independent Component Analysis）では、推定したい時間を含む数十分もしくは数時間の時間的に連続したデータ（有効電力と無効電力）を推定に活用している。しかし本例では、データが時間的に連続していなくてもよく、例えば、前日の13時〜14時、前々日の13時〜14時のように時間的に連続していないデータを抽出してもよい。なお抽出したデータは、具体的には、以下のｎ個のデータである。

有効電力計測値i ，無効電力計測値i
ただしi : 1〜ｎ
復元行列構築処理部３５は、上記のようにして抽出したデータ（有効電力計測値iと無効電力計測値i）を用いて復元行列を構築する処理を行う。本復元行列構築処理は、通常のICAの一般的な処理であるため、詳細な説明を割愛する。

推定処理部１（３６）は、復元行列構築処理部３５で構築した復元行列を用いて実負荷が消費する電力(有効電力)と太陽光発電出力を分離する処理を行う。本推定処理部１（３６）の処理は、通常の独立成分分析を用いた分離処理であるので、詳細な説明を割愛する(必要であれば上記した特許文献５参照)。なお、本推定処理部１（３６）の処理によって出力されるデータは、以下のものになる。すなわち、
実負荷の有効電力推定値（独立成分分析）i (ただしi : 1〜ｎ)
回帰式構築処理部３７は、上記データ抽出処理部３４で抽出した無効電力計測値と、上記推定処理部１（３６）の処理によって出力された実負荷の有効電力推定値（独立成分分析）を用いて回帰式を構築する。回帰式の構築に当たっては、図１５に示すように、有効電力をＹ座標に、無効電力をＸ座標に設定し、当該座標上に上記推定処理部１（３６）の処理によって出力された実負荷の有効電力推定値（独立成分分析）をプロットすることでその傾きから次のような回帰式を得る。つまり、次式のa0とa1を求める。

実負荷の有効電力推定値（独立成分分析）＝a0+a1無効電力計測値・・・式３
そして推定処理部２（３８）は、上記の回帰式(式３)を用いて、抽出したデータ（データ抽出処理部３４で抽出した無効電力計測値，有効電力計測値）に基づく推定値を求める。つまり、
実負荷の有効電力推定値（回帰）i＝a0+a1無効電力計測値i
太陽光発電の有効電力推定値i=実負荷の有効電力推定値（回帰）i-有効電力計測値i
推定幅計算処理部３９は、上記推定処理部２（３８）が上記回帰式(式３)を用いて得た上記推定値の幅を求める。幅を求める式は、以下に示すとおりである。

実負荷の有効電力推定値（回帰）の幅i=実負荷の有効電力推定値（回帰）i±幅i
太陽光発電の有効電力推定値の幅i=実負荷の有効電力推定値（回帰）の幅i -有効電力計測値i
幅i=t(ｎ−２，信頼区間) × √（Ve）×√（１＋1/n＋（無効電力計測値i-avg(無効電力計測値i)）²/Sxx）
Ve=Σ（実負荷の有効電力推定値（回帰）i-実負荷の有効電力推定値（独立成分分析）i）²／(n-2)
Sxx=Σ（無効電力計測値i-avg(無効電力計測値i)）²
ここで、t(ｎ−２，信頼区間)は、ｔ分布表から求められる値であり、ｎが無限大の場合には正規分布になり、ｎが小さい場合には正規分布よりも幅が広くなる分布である。例えば、ｎが12で信頼区間が95％の場合には、2.228になり、ｎが無限大で信頼区間が95％の場合には1.960になる。なお、Veは、残差による分散であり、上式のとおり、実負荷の有効電力推定値（回帰）と独立成分分析で求めた有効電力推定値の差を意味する。また、Sxxは、偏差平方和と呼ばれる統計指標である。

図１６は、本発明の実施形態に係る電力推定装置(その３)における実負荷の有効電力推定値が採る推定値の幅を示す概念図である。図１６において、○印は、実負荷の有効電力推定値（独立成分分析）iを示し、△印は、実負荷の有効電力推定値（回帰）iを示し、実負荷の有効電力推定値（回帰）iは△印を結んだ直線上の値となる。また図１６において、左上及び右下の曲線は上記幅iを求める計算式によって既定される推定値の確からしさを担保する境界を示し、当該二つの曲線の乖離が両矢印線によって推定値の幅iとして示されている。

１ 変電所
10〜30 電力推定装置
11、21、31 入力処理部
12、22、32 データ蓄積処理部
13、23、33 出力処理部
14 事故区間判定部
15 推定処理部
16 推定幅処理部
17 事故復旧処理部
18 バス
24、34 データ抽出処理部
25、37 回帰式構築処理部
26 推定処理部
27、39 推定幅計算処理部
35 復元行列構築処理部
36 推定処理部１
38 推定処理部２
200 無効電力補償装置
300 開閉器S5
400 開閉器S6
P1〜P4 電力計
S1〜S4 (区分)開閉器

Claims (13)

1. 複数の分散電源と負荷が配電線に接続された電力系統の事故復旧方法において、
   複数の区分開閉器を設けて前記配電線を複数の区間に区分するステップ、
   前記配電線上の前記区間ごとの負荷を計測して推定手段に設けた記憶装置に記憶すると共に配電系統の最大許容容量を予め該記憶装置に記憶するステップ、
   前記推定手段に設けた記憶装置に記憶された前記計測負荷を基に前記区間ごとの事故直前の分散電源出力と負荷の消費電力を推定するステップ、
   事故後、前記配電系統の事故区間を特定するステップ、および、
   前記配電系統の事故区間を特定した後、事故区間の負荷の消費電力推定値を用い、前記配電系統の最大許容容量を超えない範囲で、前記推定手段が前記区分開閉器を操作する事故復旧するステップ、
   を含むことを特徴とする配電系統の事故復旧方法。
2. 前記推定手段が推定する前記負荷の消費電力推定値には、前記事故復旧時に過負荷とならない推定値の幅（余裕率）を持たせ、前記事故復旧するステップは、前記推定値の幅（余裕率）を含んだ前記負荷の消費電力推定値のみを事故復旧することを特徴とする請求項１に記載の配電系統の事故復旧方法。
3. 前記区間ごとの事故直前の分散電源出力と負荷の消費電力を推定するステップは、前記記憶装置に記憶した前記区間ごとの前記負荷装置の消費電力、および該当時の気温を座標軸にプロットして得た散布図を用いて、前記負荷装置の消費電力と前記分散電源の出力を推定することを特徴とする請求項１に記載の配電系統の事故復旧方法。
4. 前記事故区間内の各区間の推定消費電力から前記分散電源の推定出力を除外した値を事故復旧に用いることを特徴とする請求項１ないし３のいずれか一項に記載の配電系統の事故復旧方法。
5. 高圧側に連系される高圧系分散電源の出力を計測する手段を該高圧系分散電源装置が備える場合には、
   前記事故区間の各区間の推定消費電力から分散電源の推定出力を減算し、その値に前記高圧系分散型電源の出力の計測値を加算した値を事故復旧に用いることを特徴とする請求項１ないし３のいずれか一項に記載の配電系統の事故復旧方法。
6. 前記推定手段が推定する前記負荷の消費電力推定値には、前記事故復旧時に過負荷とならない推定値の幅（余裕率）を持たせ、前記負荷の消費電力推定値分のみ事故復旧するステップは、前記推定値の幅（余裕率）を含んだ前記負荷の消費電力推定値のみを事故復旧することを特徴とする請求項４または５に記載の配電系統の事故復旧方法。
7. 複数の分散電源及び複数の負荷装置が配電線に接続された電力系統において、
   複数の区分開閉器を設けて前記配電線を複数の区間に区分する手段と、
   前記配電線上で前記区間ごとの有効電力および無効電力を計測して推定手段に設けた記憶装置に記憶する手段と、
   前記記憶装置に記憶された前記区間ごとの前記分散電源および前記負荷装置の無効電力から前記区間ごとの負荷装置の有効電力を回帰式で推定する手段と、
   前記回帰式で推定した前記区間ごとの負荷装置の有効電力推定値に対してその推定値の幅を算出する手段と、
   を備えることを特徴とする配電系統の実負荷推定装置。
8. 前記回帰式で求められた区間ごとの負荷装置の有効電力推定値を基に、前記分散電源の有効電力を推定する手段と、
   前記回帰式で推定した前記区間ごとの分散電源の有効電力推定値に対してその推定値の幅を算出する手段と、
   をさらに備えることを特徴とする請求項７に記載の配電系統の実負荷推定装置。
9. 複数の分散電源及び複数の負荷装置が配電線に接続された電力系統において、
   複数の区分開閉器を設けて前記配電線を複数の区間に区分するステップ、
   前記配電線上で前記区間ごとの有効電力および無効電力を計測して推定手段に設けた記憶装置に記憶するステップ、
   前記記憶装置に記憶された前記区間ごとの前記分散電源および前記負荷装置の無効電力から前記区間ごとの負荷装置の有効電力を回帰式で推定するステップ、および、
   前記回帰式で推定した前記区間ごとの負荷装置の有効電力推定値に対してその推定値の幅を算出するステップ、
   を含むことを特徴とする配電系統の実負荷推定方法。
10. 前記回帰式で求められた区間ごとの負荷装置の有効電力推定値を基に、前記分散電源の有効電力を推定するステップ、
    前記回帰式で推定した前記区間ごとの分散電源の有効電力推定値に対してその推定値の幅を算出するステップ、
    を含むことを特徴とする請求項９に記載の配電系統の実負荷推定方法。
11. 前記計測して有効電力と無効電力から独立成分分析を用いて、前記負荷装置の有効電力、もしくは、前記分散電源の出力を分離するステップと、
    該分離した前記負荷装置の有効電力、もしくは、前記分散電源の有効電力と、前記計測した無効電力とから、前記回帰式を用いて、前記負荷装置の有効電力、もしくは、前記分散電源の有効電力を推定するステップと、
    該推定した前記負荷装置の有効電力、もしくは、前記分散電源の有効電力の推定値の幅を算出するステップと、
    さらに含むことを特徴とする請求項９または１０に記載の配電系統の実負荷推定方法。
12. 複数の分散電源及び複数の負荷装置が配電線に接続された電力系統において、
    該電力系統に配電系統の実負荷を推定する推定装置を備え、該推定装置のコンピュータを、
    複数の区分開閉器を設けて前記配電線を複数の区間に区分する手段と、
    前記配電線上で前記区間ごとの有効電力および無効電力を計測して記憶装置に記憶する手段と、
    前記記憶装置に記憶された前記区間ごとの前記分散電源および前記負荷装置の無効電力から前記区間ごとの負荷装置の有効電力を回帰式で推定する手段と、
    前記回帰式で推定した前記区間ごとの負荷装置の有効電力推定値に対してその推定値の幅を算出する手段、
    として機能させるためのプログラム。
13. 複数の分散電源及び複数の負荷装置が配電線に接続された電力系統において、
    該電力系統に配電系統の実負荷を推定する推定装置を備え、該推定装置のコンピュータを、
    複数の区分開閉器を設けて前記配電線を複数の区間に区分する手段と、
    前記配電線上で前記区間ごとの有効電力および無効電力を計測して記憶装置に記憶する手段と、
    前記記憶装置に記憶された前記区間ごとの前記分散電源および前記負荷装置の無効電力から前記区間ごとの負荷装置の有効電力を回帰式で推定する手段と、
    前記回帰式で推定した前記区間ごとの負荷装置の有効電力推定値に対してその推定値の幅を算出する手段と、
    前記回帰式で求められた区間ごとの負荷装置の有効電力推定値を基に、前記分散電源の有効電力を推定する手段と、
    前記回帰式で推定した前記区間ごとの分散電源の有効電力推定値に対してその推定値の幅を算出する手段、
    として機能させるためのプログラム。