JP7809603B2 - 電力制御システム、電力制御方法、および、プログラム - Google Patents

Description

本発明の実施形態は、電力制御システム、電力制御方法、および、プログラムに関する。

近年、太陽光や風力などの再生可能エネルギーを利用する電源設備の導入が拡大している。再生可能エネルギーを利用する電源設備は出力変動を伴うので、電力系統の需給バランスを維持するための調整力の重要性が増している。

例えば、分散リソース（電力に関する分散型のリソース）を利用したＶＰＰ（Virtual Power Plant：仮想発電所）の運用で調整力を創出することが求められる。分散リソースとしては、分散電源や分散蓄電池のほか、需要応動を提供する需要家リソース（需要家が所有する電力消費リソース）が考えられる。そのようなＶＰＰでは、応動量の確実性や精度が異なる複数の分散リソースの応動量を高精度に束ねる電力制御手法（アグリゲーション制御手法）が必要となる。

例えば、従来技術で、負荷情報とレスポンス率に基づいて需要家リソースに負荷調整を要請する電力調整システムや、カレンダー情報や気象予測情報を利用して電力調整を分散リソースに配分するアグリゲータシステムが知られている。

特許第５６０３３１８号公報

電力需給調整市場での市場取引を介した調整力提供においては、需給調整のための調整力指令に所定時間内に所定誤差で追従する調整力応動の創出が必要である。また、小規模系統の需給バランスを確保するための調整力提供において、逐次に変化する電力需要に追従する調整力の創出が必要である。

しかしながら、上述の従来技術は、主に分散リソースの事前情報をもとにリソース応動量を予想して運用するものであるが、リソースの状態に応じた応動量の調整や応動量指令変化に対する応動量の追従を高速に行うことには適していない。具体的には、例えば、分散リソースごとの動特性や伝達むだ時間（情報伝達に伴う遅延時間）を考慮して応動量を高速かつ安定に制御することはできていない。

そこで、本発明の実施形態は、複数の分散リソースに関して高速かつ安定な応動制御を実現可能な電力制御システム、電力制御方法、および、プログラムを提供することを課題とする。

実施形態の電力制御システムは、電力に関する複数の分散リソースの全体に対する応動指令である集約応動指令に、複数の前記分散リソースそれぞれによる応動量を合算した集約応動量を追従させるための電力制御システムであって、外部装置から取得した前記集約応動指令における応動量を、少なくとも前記分散リソースごとの応動可能量に基づいて、前記分散リソースごとに配分することで、複数のリソース配分指令を出力する配分演算部と、前記分散リソースごとに、前記リソース配分指令と、応動の観測量であるリソース観測量と、前記分散リソースの動特性および伝達むだ時間に関するモデルと、に基づいて、応動に関する前記分散リソースの状態量であるリソース状態量の推定値を算出して出力する状態推定部と、前記分散リソースごとに、前記リソース配分指令と、前記リソース状態量の推定値と、に基づいて、前記分散リソースに対する応動の制御指令であるリソース制御指令を出力する制御演算部と、を備える。

図１は、実施形態の電力制御システムの全体構成図である。 図２は、実施形態におけるリソース配分例を示す図である。 図３は、実施形態における各モデルを用いた演算の説明図である。 図４は、実施形態におけるオブザーバの演算の説明図である。 図５は、実施形態における応動可能量の説明図である。 図６は、実施形態における制御演算部の演算の説明図である。 図７は、実施形態の各分散リソースの性能等の例を示す図である。 図８は、図７を前提にしたシミュレーション結果の例を示す図である。 図９は、実施形態の電力制御システムによる処理の概要を示すフローチャートである。

以下、本発明の電力制御システム、電力制御方法、および、プログラムの実施形態について、図面を参照して説明する。なお、以下では、分散リソースを単にリソースと称する場合がある。

図１は、実施形態の電力制御システムＳの全体構成図である。電力制御システムＳは、電力に関する複数の分散リソース５０の全体に対する応動指令である集約応動指令Ｒ_ｓに、複数の分散リソース５０それぞれによる応動量を合算した集約応動量Ｙ_ｓを追従させるためのシステムである。電力制御システムＳは、配分演算部１０と、制御演算部２０と、状態推定部３０と、複数の分散リソース５０と、を備える。

なお、配分演算部１０、制御演算部２０、状態推定部３０は、情報処理を実行可能な装置（例えばコンピュータ装置）である。そして、それらのうち２つ以上が同一のハードウェアで構成されていてもよいし、あるいは、それらが別々のハードウェアで構成されていてもよい。ハードウェアとしては、少なくとも処理部と記憶部を備える。処理部は、例えば、ＣＰＵ（Central Processing Unit）によって構成され、各種演算処理を実行する。記憶部は、例えば、ＲＡＭ（Random Access Memory）、ＲＯＭ（Read Only Memory）、ＨＤＤ（Hard Disk Drive）、ＳＳＤ（Solid State Drive）などによって構成され、各種情報（電力制御に関する情報等）を記憶する。また、ハードウェアとして、入力部（例えば、マウス、キーボード、タッチパネル等）や表示部（ＬＣＤ（Liquid Crystal Display）等）などを備えていてもよい。

分散リソース５０として、図１の例では、分散リソース１～３を備えている。なお、分散リソース１～３を区別しないときは、分散リソース５０と称する場合もある。分散リソース１～３は、需要家が所有する需要家リソースであって、受電量を減少させることで、需要の予測値である基準値（過去データ等から予測）に対する応動量を創出する。

分散リソース１は、電力消費設備（不図示）のほかに、自家発ＰＧ１（自家発電設備）を備えている。分散リソース１は、制御演算部２０から受電量減少指示（リソース制御指令ｕ_１）を受けた場合に、自家発ＰＧ１を稼動させることで応動（応動量を創出）できる。分散リソース２は、分散リソース１と同様であるので、説明を省略する。

分散リソース３は、電力消費設備（不図示）のほかに、蓄電池ＳＢを備えている。分散リソース３は、制御演算部２０から受電量減少指示（リソース制御指令ｕ_３）を受けた場合に、蓄電池ＳＢから放電させることで応動（応動量を創出）できる。

なお、分散リソース５０によって応動量を創出する方法は、これらに限定されず、ほかに、例えば、電力消費設備の稼働調整（停止など）により応動量を創出するようにしてもよい。

配分演算部１０は、外部装置から取得した集約応動指令Ｒ_ｓにおける応動量を、少なくとも分散リソース５０ごとの応動可能量に基づいて、分散リソース５０ごとに配分することで、複数のリソース配分指令ｒ_ｉを出力（算出および出力。以下同様）する。その際に、分散リソース５０の使用順位（メリットオーダー）の情報も用いてもよい。使用順位は、例えば、応動量を創出する際の経済性や負担（例えば設備稼働停止に関する負担）の大きさ等に基づいて決定することができる。

ここで、図２は、実施形態におけるリソース配分例を示す図である。ここでは、図２（ａ）に示すように、使用順位の高い（運用コストの低い）ほうから順番に、分散リソース１、２、３となっている。また、分散リソース１の応動可能量の予測値は符号８１に示す幅である。その場合、集約応動指令Ｒ_ｓにおける応動量を、分散リソース１、２に配分することができる。

ただし、図２（ｂ）に示すように、分散リソース１の応動可能量の実際の値が、符号８１に示す幅よりも狭い符号８２に示す幅になる場合も考えられる。その場合、集約応動指令Ｒ_ｓにおける応動量を、分散リソース１、２、３に配分することで対応することができる。

なお、図２では、使用順位も考慮する場合について説明したが、これに限定されない。使用順位を考慮しない場合は、例えば、それぞれの分散リソース５０が均等の分担率になるように配分するようにしてもよい。

図１に戻って、制御演算部２０は、リソース制御器２１ａ～２１ｃを備える。なお、制御対象として分散リソース１～３を特に区別しないときには、リソース制御器２１ａ～２１ｃの代わりに制御演算部２０と称する場合もある。

制御演算部２０は、分散リソース５０ごとに、リソース配分指令ｒ_ｉ（i＝１，２，３）と、リソース状態量ｘ_ｉ（応動に関する分散リソース５０の状態量）の推定値と、に基づいて、分散リソース５０に対する応動の制御指令であるリソース制御指令ｕ_ｉを出力する。その際に、制御演算部２０は、分散リソース５０からフィードバックしたリソース観測量ｙ_Ｌｉ（応動の観測量）も併せて用いて、リソース配分指令ｒ_ｉとリソース応動量ｙ_ｉが一致するように（近づくように）リソース制御指令ｕ_ｉを算出する。

また、基準値と実際の需要量の偏差を含むリソース外乱も考慮する場合、制御演算部２０は、分散リソース５０ごとに、リソース配分指令ｒ_ｉと、リソース外乱ｄ_ｉを含むリソース状態量ｘ_ｉの推定値と、に基づいて、リソース制御指令ｕ_ｉを出力する。

状態推定部３０は、オブザーバ３１ａ～３１ｃを備える。なお、制御対象として分散リソース１～３を特に区別しないときには、オブザーバ３１ａ～３１ｃの代わりに状態推定部３０と称する場合もある。

状態推定部３０は、分散リソース５０ごとに、リソース配分指令ｒ_ｉと、リソース観測量ｙ_Ｌｉと、分散リソース５０の動特性および伝達むだ時間に関するモデルと、に基づいて、リソース外乱ｄ_ｉを含むリソース状態量ｘ_ｉの推定値を出力する。

なお、伝達むだ時間とは、例えば、図１において、符号ＩＮ１に示す入力むだ時間や、符号ＯＵＴ１に示す出力むだ時間である。符号ＩＮ１や符号ＯＵＴ１の箇所では、情報の伝達に遅延が生じる。

ここで、図３は、実施形態における各モデルを用いた演算の説明図である。図４は、実施形態におけるオブザーバ３１の演算の説明図である。

状態推定部３０は、３つのモデルと、外乱di(k)を状態量xi(k)に含めた式（１）の状態方程式と、を用いて、分散リソース５０ごとの状態量x_i(k)の推定値x^_i(k)を制御周期（例えば１０～６０秒程度）ごとに算出する（図４）。kは制御ステップである。また、x^_i(k)がx_i(k)の推定値を示すように、各変数について、「^」は推定値であることを示す記号である。ただし、推定値について「^」の記号を省略する場合もある。
ここで、各文字の定義は以下の式（2）～（5）の通りである。

例えば、入力むだ時間モデルM_uiについて、入力するリソース制御指令u_i(k)と出力するu_Li(k)との関係は、以下の式（６）、（７）の通りである。

また、リソースモデルM_riと出力むだ時間モデルM_yiについても、同様である。

また、式（１）におけるHiとしては、以下の式（８）に示すカルマンゲイン（最適オブザーバゲイン）を用いることができる。

ここで、Siは、以下の式（９）に示すリッカチ方程式の解である。Qoi、Roiはそれぞれ想定するシステムノイズ、観測ノイズの共分散行列であり、制御系の設計定数である。

オブザーバ３１では、計測できない情報を、他の計測できる情報からモデルを用いて推定するが、制御対象物の観測値にシステムノイズや観測ノイズが含まれており、カルマンゲインによってそれらのノイズがあっても精度よく推定できる。

例えば、システムノイズは、以下の式（１０）におけるｖ_i(k)である。

また、例えば、観測ノイズは、以下の式（１１）におけるｗ_i(k)である。

また、実時間（むだ時間がない場合の時間）のリソース応動量の推定値ｙ^_i(k)は、以下の式（１２）、（１３）で算出することができる。

ここで、図６は、実施形態における制御演算部２０の演算の説明図である。符号Ｒ１は、積分を示す。符号Ｒ２、Ｒ３は、乗算を示す。制御演算部２０は、図３の各モデルを対象に、以下の式（１４）～（１７）に示すように、ｒ_iとｙ_iを早く近づけるために積分値「∫（ｒ_i―ｙ_i）」を含む演算を行う。

そして、制御演算部２０は、以下の式（１８）～（２１）を用いて、最適サーボ系のリソース制御指令としてui(k)を算出する。Xi(k)はオブザーバの推定値x^_i(k)を用いて算出する。

なお、F_1iはx^_i(k)に乗算する項をまとめたものである。また、F_2iは∫(r_i(k)－y^i(k))に乗算する項をまとめたものである。

また、Piは、以下の式（２２）に示すリッカチ方程式の解である。なお、Qi、Riは重み行列で制御系の設計定数である。

上述のようにして、制御演算部２０は、分散リソース５０ごとに、リソース状態量ｘ_ｉの推定値に基づいてリソース応動量ｙ_ｉの推定値を算出し、実時間のリソース応動量の推定値ｙ_ｉがリソース配分指令ｒ_ｉに追従する（近づく）ようにリソース制御指令ｕ_ｉを算出する。

また、配分演算部１０は、リソース状態量ｘ_ｉの推定値におけるリソース外乱ｄ_ｉに基づいて、リソース配分指令ｒ_ｉを補正する。例えば、配分演算部１０は、分散リソース５０の需要変動や、基準値と実際の需要量の偏差の定常成分（例えば、１０～２０分程度の平均値）に合わせて、基準値を基準にした実質のリソース応動可能量に合うようにリソース配分指令ｒ_ｉを補正する。

ここで、図５は、実施形態における応動可能量の説明図である。例えば、分散リソース１について、需要の予測値である基準値ＳＶに対して、計画した応動可能量は符号Ｑ１に示す通りである。しかし、実際の需要量ＡＤは、基準値ＳＶより大きい値となっている。その場合、実際の応動可能量は、符号Ｑ２に示す通りになる。ただし、基準値ＳＶを基準にした実質の（見かけの）応動量は、符号Ｑ３に示す通りになる。その場合、配分演算部１０は、基準値ＳＶを基準にした実質のリソース応動可能量に合うようにリソース配分指令ｒ_ｉを補正する（図２参照）。

図１に戻って、また、状態推定部３０は、カルマンゲイン（式（８））を用いたオブザーバ３１に基づいて、リソース状態量ｘ_ｉの推定値を出力する。

また、状態推定部３０は、蓄電リソース（蓄電池ＳＢを含む分散リソース３）について、リソース状態量ｘ_ｉの推定値として、蓄電量および／または蓄電量変化を含む蓄電推定値を出力する。その場合、配分演算部１０は、蓄電リソースについて、蓄電推定値の低下に応じて応動可能量を小さくすることで、応動量が低下するようにリソース配分指令ｒ_ｉを変更する。

つまり、蓄電池ＳＢの放電で応動量を創出する蓄電池リソースの実時間の蓄電量を併せて推定し、蓄電量の低下時には蓄電リソースの配分を小さくし、配分に余裕がある他の分散リソースにその分の配分を大きくするように電力制御（アグリゲーション制御）を行うことができる。蓄電量の推定値は、特定時点の蓄電量の初期値とリソース制御指令ｕ_ｉ（放電出力の指令）の情報をもとに、以下の式（２３）に示すモデルで算出する。

ここで、z^i(k)は蓄電量zi(k)の推定値である。η_iは蓄電設備の放電効率である。Δtは制御周期である。

次に、図７、図８を参照して、シミュレーション例について説明する。図７は、実施形態の各分散リソースの性能等の例を示す図である。図８は、図７を前提にしたシミュレーション結果の例を示す図である。

図７に示すように、分散リソース１～３について、応動可能量、リソース動特性、リソース使用順位、需要変動が設定されている。そして、図８（ｂ）では、分散リソース１に関して、符号１１が基準値を示し、符号１２が基準値からのリソース制御指令を示し、符号１３が基準値からのリソース応動量を示す。また、図８（ｃ）では、分散リソース２に関して、符号２１が基準値を示し、符号２２がリソース制御指令を示し、符号２３がリソース応動量を示す。

また、図８（ｄ）では、分散リソース３に関して、符号３１が基準値を示し、符号３２がリソース制御指令を示し、符号３３がリソース応動量を示す。また、図８（ａ）において、符号Ｖ１が集約応動指令を示し、符号Ｖ２が集約応動量を示す。

図８（ａ）～（ｄ）からわかるように、複数の分散リソース１～３に関して高速かつ安定な応動制御を実現している。つまり、図８（ａ）に示すように、集約応動指令に対して集約応動量が高速に追従している。

次に、図９は、実施形態の電力制御システムＳによる処理の概要を示すフローチャートである。まず、ステップＳ１において、配分演算部１０は、外部装置から集約応動指令Ｒ_ｓを取得する。

次に、ステップＳ２において、配分演算部１０は、集約応動指令Ｒ_ｓにおける応動量を、分散リソース５０ごとの応動可能量や使用順位に基づいて、分散リソース５０ごとに配分することで、複数のリソース配分指令ｒ_ｉを出力する。

次に、ステップＳ３において、状態推定部３０は、分散リソース５０ごとに、リソース配分指令ｒ_ｉと、リソース観測量ｙ_Ｌｉと、分散リソース５０の動特性および伝達むだ時間に関するモデル（図３、図４）と、に基づいて、リソース外乱ｄ_ｉを含むリソース状態量ｘ_ｉの推定値を算出して出力する。

次に、ステップＳ４において、制御演算部２０は、分散リソース５０ごとに、リソース配分指令ｒ_ｉと、リソース外乱ｄ_ｉを含むリソース状態量ｘ_ｉの推定値と、に基づいて、リソース制御指令ｕ_ｉを出力する。

このように、本実施形態の電力制御システムＳによれば、分散リソース５０の動特性および伝達むだ時間に関するモデルを用いて上述の処理を行うことによって、複数の分散リソース５０に関して高速かつ安定な応動制御を実現することができる。

具体的には、実時間のリソース応動量の推定値に基づいて分散リソース５０の制御を行うので、各分散リソース５０の応動量情報の収集に時間遅延があるシステムにおいても、速応性に優れるリソース応動を実現することができる。

また、伝達むだ時間を含む分散リソース状態量の推定値を利用して制御するので、フィードバック安定性に優れる制御を実現することができる。

また、分散リソース応動量に影響を与える外乱成分の推定値に基づいて、応動量指令のリソース配分を補正するとともに、リソース配分に一致するようにリソース制御を行うので、需要家リソースの需要変動や基準値偏差の変化影響を小さく抑えることで高精度の応動を実現することができる。

さらに、例えば、需要設備とともに太陽光発電設備等を併設した需要家リソースでは、気象変化に伴う太陽光発電量の変化が外乱成分となるが、同様の構成でこの影響を小さく抑えて応動量の精度を向上させることができる。

また、蓄電設備の放電で応動量を創出する蓄電池リソースの実時間の蓄電量を推定し、蓄電量の低下に応じてリソース配分指令を小さくする、あるいは停止するので、蓄電量低下により蓄電リソースが放電を停止する際の集約応動量への影響を抑えて応動量の精度を向上させることができる。

本実施形態の電力制御システムＳで実行されるプログラムは、インストール可能な形式又は実行可能な形式のファイルでＣＤ－ＲＯＭ、フレキシブルディスク（FD）、ＣＤ－Ｒ、ＤＶＤ（Digital Versatile Disk）等のコンピュータで読み取り可能な記録媒体に記録されて提供されるようにしてもよい。

また、当該プログラムを、インターネット等のネットワークに接続されたコンピュータ上に格納し、ネットワーク経由でダウンロードさせることにより提供するように構成してもよい。また、当該プログラムをインターネット等のネットワーク経由で提供または配布するように構成してもよい。また、当該プログラムを、ＲＯＭ等に予め組み込んで提供するようにしてもよい。

本発明の実施形態を説明したが、この実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。この新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。この実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。

例えば、上述の実施形態では、需要家リソースを電力制御対象として取り上げたが、これに限定されず、例えば、応動量情報の収集に時間遅延がある発電リソース等のアグリゲーション制御に対して本発明を同様に適用することができる。

また、リソース応動量に影響を与える外乱が小さく、制御入力への追従性もよい分散リソース５０に対しては、制御演算部２０に一般のＰＩＤ制御を適用し、さらに状態推定部３０は省略してよい。すなわち、全ての分散リソース５０に状態推定部３０を備える構成としなくてもよい。

１０…配分演算部、２０…制御演算部、３０…状態推定部、３１…オブザーバ、５０…分散リソース、Ｓ…電力制御システム

Claims (9)

1. 電力に関する複数の分散リソースの全体に対する応動指令である集約応動指令に、複数の前記分散リソースそれぞれによる応動量を合算した集約応動量を追従させるための電力制御システムであって、
   外部装置から取得した前記集約応動指令における応動量を、少なくとも前記分散リソースごとの応動可能量に基づいて、前記分散リソースごとに配分することで、複数のリソース配分指令を出力する配分演算部と、
   前記分散リソースごとに、前記リソース配分指令と、応動の観測量であるリソース観測量と、前記分散リソースの動特性および伝達むだ時間に関するモデルと、に基づいて、応動に関する前記分散リソースの状態量であるリソース状態量の推定値を算出して出力する状態推定部と、
   前記分散リソースごとに、前記リソース配分指令と、前記リソース状態量の推定値と、に基づいて、前記分散リソースに対する応動の制御指令であるリソース制御指令を出力する制御演算部と、を備える電力制御システム。
2. 前記状態推定部は、前記分散リソースごとに、前記リソース状態量の推定値を算出する場合に、応動量に影響を与える外乱であるリソース外乱を含む前記リソース状態量の推定値を出力し、
   前記制御演算部は、前記分散リソースごとに、前記リソース配分指令と、前記リソース外乱を含む前記リソース状態量の推定値と、に基づいて、前記リソース制御指令を出力する、請求項１に記載の電力制御システム。
3. 前記分散リソースは、需要家が所有する前記分散リソースである需要家リソースであって、受電量を低減することで、需要の予測値である基準値に対する応動量を創出する需要家リソースであり、
   前記リソース外乱は、前記基準値と実際の需要量の偏差を含む、請求項２に記載の電力制御システム。
4. 前記制御演算部は、前記分散リソースごとに、前記リソース状態量の推定値に基づいて実時間の応動量であるリソース応動量の推定値を算出し、実時間の前記リソース応動量の推定値が前記リソース配分指令に追従するように前記リソース制御指令を出力する、請求項１に記載の電力制御システム。
5. 前記配分演算部は、前記リソース状態量の推定値における前記リソース外乱に基づいて、前記リソース配分指令を補正する、請求項２に記載の電力制御システム。
6. 前記状態推定部は、さらに、カルマンゲインを用いたオブザーバに基づいて、前記リソース状態量の推定値を算出する、請求項１に記載の電力制御システム。
7. 前記分散リソースは、蓄電設備の放電により応動量を創出する蓄電リソースを含み、
   前記状態推定部は、前記蓄電リソースについて、前記リソース状態量の推定値として、蓄電量および／または蓄電量変化を含む蓄電推定値を出力し、
   前記配分演算部は、前記蓄電リソースについて、前記蓄電推定値の低下に応じて、応動量が低下するように前記リソース配分指令を変更する、請求項１に記載の電力制御システム。
8. 電力に関する複数の分散リソースの全体に対する応動指令である集約応動指令に、複数の前記分散リソースそれぞれによる応動量を合算した集約応動量を追従させるための電力制御方法であって、
   外部装置から取得した前記集約応動指令における応動量を、少なくとも前記分散リソースごとの応動可能量に基づいて、前記分散リソースごとに配分することで、複数のリソース配分指令を出力する配分演算ステップと、
   前記分散リソースごとに、前記リソース配分指令と、応動の観測量であるリソース観測量と、前記分散リソースの動特性および伝達むだ時間に関するモデルと、に基づいて、応動に関する前記分散リソースの状態量であるリソース状態量の推定値を算出して出力する状態推定ステップと、
   前記分散リソースごとに、前記リソース配分指令と、前記リソース状態量の推定値と、に基づいて、前記分散リソースに対する応動の制御指令であるリソース制御指令を出力する制御演算ステップと、を備える電力制御方法。
9. 電力に関する複数の分散リソースの全体に対する応動指令である集約応動指令に、複数の前記分散リソースそれぞれによる応動量を合算した集約応動量を追従させるためのプログラムであって、
   コンピュータを、
   外部装置から取得した前記集約応動指令における応動量を、少なくとも前記分散リソースごとの応動可能量に基づいて、前記分散リソースごとに配分することで、複数のリソース配分指令を出力する配分演算部と、
   前記分散リソースごとに、前記リソース配分指令と、応動の観測量であるリソース観測量と、前記分散リソースの動特性および伝達むだ時間に関するモデルと、に基づいて、応動に関する前記分散リソースの状態量であるリソース状態量の推定値を算出して出力する状態推定部と、
   前記分散リソースごとに、前記リソース配分指令と、前記リソース状態量の推定値と、に基づいて、前記分散リソースに対する応動の制御指令であるリソース制御指令を出力する制御演算部と、して機能させるためのプログラム。