(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2023179323
(43)【公開日】2023-12-19
(54)【発明の名称】情報処理装置、情報処理方法及び情報処理プログラム
(51)【国際特許分類】
H02J 3/00 20060101AFI20231212BHJP
G06Q 50/06 20120101ALI20231212BHJP
H02J 3/38 20060101ALI20231212BHJP
H02J 13/00 20060101ALI20231212BHJP
【ＦＩ】
H02J3/00 170
G06Q50/06
H02J3/38 110
H02J13/00 301A
【審査請求】未請求
【請求項の数】8
【出願形態】ＯＬ
(21)【出願番号】P 2022092593
(22)【出願日】2022-06-07
(71)【出願人】
【識別番号】000005234
【氏名又は名称】富士電機株式会社
(74)【代理人】
【識別番号】110000176
【氏名又は名称】弁理士法人一色国際特許事務所
(72)【発明者】
【氏名】大井 章弘
【テーマコード（参考）】
5G064
5G066
5L049
【Ｆターム（参考）】
5G064AC09
5G064DA03
5G066AA03
5G066HB06
5L049CC06
(57)【要約】 （修正有）
【課題】需要家の負荷の有効電力等と、発電設備の有効電力とを精度良く推定する。
【解決手段】電力系統１の所定区間において、負荷を有する需要家Ａと、余剰買取契約対象の発電設備を有する需要家Ｂと、全量買取契約対象の発電設備を有する需要家Ｃと、のいずれかを複数含む電力系統の状態を推定する情報処理装置２であって、有効電力及び無効電力を計算する第１計算部と、各需要家が需要家Ａ～Ｃのいずれに該当するかを示す分類情報と、各需要家の電力量計の測定値とに基づいて、各需要家のうちいずれかの有効電力を計算する第２計算部と、需要家Ａ及び需要家Ｂの有効電力の第１比率と、需要家Ｂ及び需要家Ｃの発電設備の有効電力の第２比率と、発電設備の力率と、有効電力及び無効電力と、各需要家のうちのいずれかの有効電力に基づいて、所定区間の負荷の有効電力及び無効電力と、所定区間の発電設備の有効電力と、を含む状態を推定する第１推定部を備える。
【選択図】図１
【解決手段】電力系統１の所定区間において、負荷を有する需要家Ａと、余剰買取契約対象の発電設備を有する需要家Ｂと、全量買取契約対象の発電設備を有する需要家Ｃと、のいずれかを複数含む電力系統の状態を推定する情報処理装置２であって、有効電力及び無効電力を計算する第１計算部と、各需要家が需要家Ａ～Ｃのいずれに該当するかを示す分類情報と、各需要家の電力量計の測定値とに基づいて、各需要家のうちいずれかの有効電力を計算する第２計算部と、需要家Ａ及び需要家Ｂの有効電力の第１比率と、需要家Ｂ及び需要家Ｃの発電設備の有効電力の第２比率と、発電設備の力率と、有効電力及び無効電力と、各需要家のうちのいずれかの有効電力に基づいて、所定区間の負荷の有効電力及び無効電力と、所定区間の発電設備の有効電力と、を含む状態を推定する第１推定部を備える。
【選択図】図１
【特許請求の範囲】
【請求項1】
電力系統の所定区間において、負荷を有する第１需要家と、余剰買取契約対象の発電設備を有する第２需要家と、全量買取契約対象の前記発電設備を有する第３需要家と、のいずれかに分類される需要家を複数含む前記電力系統の状態を推定する情報処理装置であって、
前記電力系統の配電線のセンサからの測定値に基づいて、前記所定区間における有効電力及び無効電力を計算する第１計算部と、
前記複数の需要家のそれぞれが前記第１～第３需要家のいずれに該当するかを示す分類情報と、前記複数の需要家の電力量計の測定値と、に基づいて、前記第１～第３需要家のうち少なくともいずれかの有効電力を計算する第２計算部と、
前記第１及び第２需要家のそれぞれの負荷における有効電力の第１比率と、前記第２及び３需要家のそれぞれの前記発電設備における有効電力の第２比率と、前記発電設備の力率と、前記所定区間における有効電力及び無効電力と、前記第１～第３需要家のうち少なくともいずれかの有効電力と、に基づいて、前記所定区間に含まれる負荷の有効電力及び無効電力と、前記所定区間に含まれる前記発電設備の有効電力と、を含む前記状態を推定する第１推定部と、
を備える、情報処理装置。
前記電力系統の配電線のセンサからの測定値に基づいて、前記所定区間における有効電力及び無効電力を計算する第１計算部と、
前記複数の需要家のそれぞれが前記第１～第３需要家のいずれに該当するかを示す分類情報と、前記複数の需要家の電力量計の測定値と、に基づいて、前記第１～第３需要家のうち少なくともいずれかの有効電力を計算する第２計算部と、
前記第１及び第２需要家のそれぞれの負荷における有効電力の第１比率と、前記第２及び３需要家のそれぞれの前記発電設備における有効電力の第２比率と、前記発電設備の力率と、前記所定区間における有効電力及び無効電力と、前記第１～第３需要家のうち少なくともいずれかの有効電力と、に基づいて、前記所定区間に含まれる負荷の有効電力及び無効電力と、前記所定区間に含まれる前記発電設備の有効電力と、を含む前記状態を推定する第１推定部と、
を備える、情報処理装置。
【請求項2】
請求項１に記載の情報処理装置であって、
前記第２計算部は、前記第１～第３需要家のうち少なくともいずれか２つの有効電力を計算し、
前記第１推定部は、前記第１比率と、前記第２比率と、前記力率と、前記所定区間における有効電力及び無効電力と、前記第１～第３需要家のうち少なくともいずれか２つの有効電力に基づいて、前記状態を推定する、
情報処理装置。
前記第２計算部は、前記第１～第３需要家のうち少なくともいずれか２つの有効電力を計算し、
前記第１推定部は、前記第１比率と、前記第２比率と、前記力率と、前記所定区間における有効電力及び無効電力と、前記第１～第３需要家のうち少なくともいずれか２つの有効電力に基づいて、前記状態を推定する、
情報処理装置。
【請求項3】
請求項２に記載の情報処理装置であって、
前記第２計算部は、前記第１～第３需要家の有効電力を計算し、
前記第１推定部は、前記第１比率と、前記第２比率と、前記力率と、前記所定区間における有効電力及び無効電力と、前記第１～第３需要家の有効電力に基づいて、前記状態を推定する、
情報処理装置。
前記第２計算部は、前記第１～第３需要家の有効電力を計算し、
前記第１推定部は、前記第１比率と、前記第２比率と、前記力率と、前記所定区間における有効電力及び無効電力と、前記第１～第３需要家の有効電力に基づいて、前記状態を推定する、
情報処理装置。
【請求項4】
請求項１に記載の情報処理装置であって、
前記第１推定部は、
前記センサからの測定値に基づいて計算された前記所定区間における有効電力、及び前記状態に基づいて計算された前記所定区間における有効電力の第１の差と、
前記センサからの測定値に基づいて計算された前記所定区間における無効電力、及び前記力率と、前記状態とに基づいて計算された前記所定区間における無効電力の第２の差と、
前記電力量計の測定値に基づいて計算された前記第１～第３需要家のうち少なくともいずれかの需要家の有効電力、及び前記第１比率と、前記第２比率と、前記状態とに基づいて計算された前記少なくともいずれかの需要家の有効電力の第３の差とに基づく第１目的関数を設定し、
前記状態を第１決定変数とした際に前記第１目的関数を最小化することにより、前記状態を推定する、
情報処理装置。
前記第１推定部は、
前記センサからの測定値に基づいて計算された前記所定区間における有効電力、及び前記状態に基づいて計算された前記所定区間における有効電力の第１の差と、
前記センサからの測定値に基づいて計算された前記所定区間における無効電力、及び前記力率と、前記状態とに基づいて計算された前記所定区間における無効電力の第２の差と、
前記電力量計の測定値に基づいて計算された前記第１～第３需要家のうち少なくともいずれかの需要家の有効電力、及び前記第１比率と、前記第２比率と、前記状態とに基づいて計算された前記少なくともいずれかの需要家の有効電力の第３の差とに基づく第１目的関数を設定し、
前記状態を第１決定変数とした際に前記第１目的関数を最小化することにより、前記状態を推定する、
情報処理装置。
【請求項5】
請求項１～４のいずれか一に記載の情報処理装置であって、
前記第１比率は、前記負荷の有効電力のうち、前記第１需要家の負荷の有効電力の割合を示す第１の値と、前記第２需要家の負荷の有効電力の割合を示す第２の値であり、
前記第２比率は、前記発電設備の有効電力のうち、前記第２需要家の発電設備の有効電力の割合を示す第３の値と、前記第３需要家の発電設備の有効電力の割合を示す第４の値であり、
前記第１比率と、前記第２比率とを推定する第２推定部を更に備え、
前記第２推定部は、
前記所定区間に含まれる負荷の有効電力、及び前記第１需要家の負荷の有効電力の測定値を前記第１の値で除した値の第４の差と、
前記所定区間に含まれる発電設備の有効電力、及び前記第３需要家の発電設備の有効電力の測定値を前記第４の値で除した値の第５の差とに基づく第２目的関数を設定し、
前記第１比率と、前記第２比率と、前記所定区間に含まれる負荷の有効電力と、前記所定区間に含まれる発電設備の有効電力とを第２決定変数とした際に前記第２目的関数を最小化することにより、前記第１比率と、前記第２比率とを推定する、
情報処理装置。
前記第１比率は、前記負荷の有効電力のうち、前記第１需要家の負荷の有効電力の割合を示す第１の値と、前記第２需要家の負荷の有効電力の割合を示す第２の値であり、
前記第２比率は、前記発電設備の有効電力のうち、前記第２需要家の発電設備の有効電力の割合を示す第３の値と、前記第３需要家の発電設備の有効電力の割合を示す第４の値であり、
前記第１比率と、前記第２比率とを推定する第２推定部を更に備え、
前記第２推定部は、
前記所定区間に含まれる負荷の有効電力、及び前記第１需要家の負荷の有効電力の測定値を前記第１の値で除した値の第４の差と、
前記所定区間に含まれる発電設備の有効電力、及び前記第３需要家の発電設備の有効電力の測定値を前記第４の値で除した値の第５の差とに基づく第２目的関数を設定し、
前記第１比率と、前記第２比率と、前記所定区間に含まれる負荷の有効電力と、前記所定区間に含まれる発電設備の有効電力とを第２決定変数とした際に前記第２目的関数を最小化することにより、前記第１比率と、前記第２比率とを推定する、
情報処理装置。
【請求項6】
請求項５に記載の情報処理装置であって、
前記第２推定部は、
時間依存性を有する、前記第１比率と、前記第２比率とを、前記第２決定変数とする、
情報処理装置。
前記第２推定部は、
時間依存性を有する、前記第１比率と、前記第２比率とを、前記第２決定変数とする、
情報処理装置。
【請求項7】
情報処理装置が、
電力系統の所定区間において、負荷を有する第１需要家と、余剰買取契約対象の発電設備を有する第２需要家と、全量買取契約対象の前記発電設備を有する第３需要家と、のいずれかに分類される需要家を複数含む前記電力系統の状態を推定する情報処理方法であって、
前記電力系統の配電線のセンサからの測定値に基づいて、前記所定区間における有効電力及び無効電力を計算するステップと、
前記複数の需要家のそれぞれが前記第１～第３需要家のいずれに該当するかを示す分類情報と、前記複数の需要家の電力量計の測定値と、に基づいて、前記第１～第３需要家のうち少なくともいずれかの有効電力を計算するステップと、
前記第１及び第２需要家のそれぞれの負荷における有効電力の第１比率と、前記第２及び第３需要家のそれぞれの前記発電設備における有効電力の第２比率と、前記発電設備の力率と、前記所定区間における有効電力及び無効電力と、前記第１～第３需要家のうち少なくともいずれかの有効電力と、に基づいて、前記所定区間に含まれる負荷の有効電力及び無効電力と、前記所定区間に含まれる前記発電設備の有効電力と、を含む前記状態を推定するステップと、
を含む処理を実行する、情報処理方法。
電力系統の所定区間において、負荷を有する第１需要家と、余剰買取契約対象の発電設備を有する第２需要家と、全量買取契約対象の前記発電設備を有する第３需要家と、のいずれかに分類される需要家を複数含む前記電力系統の状態を推定する情報処理方法であって、
前記電力系統の配電線のセンサからの測定値に基づいて、前記所定区間における有効電力及び無効電力を計算するステップと、
前記複数の需要家のそれぞれが前記第１～第３需要家のいずれに該当するかを示す分類情報と、前記複数の需要家の電力量計の測定値と、に基づいて、前記第１～第３需要家のうち少なくともいずれかの有効電力を計算するステップと、
前記第１及び第２需要家のそれぞれの負荷における有効電力の第１比率と、前記第２及び第３需要家のそれぞれの前記発電設備における有効電力の第２比率と、前記発電設備の力率と、前記所定区間における有効電力及び無効電力と、前記第１～第３需要家のうち少なくともいずれかの有効電力と、に基づいて、前記所定区間に含まれる負荷の有効電力及び無効電力と、前記所定区間に含まれる前記発電設備の有効電力と、を含む前記状態を推定するステップと、
を含む処理を実行する、情報処理方法。
【請求項8】
電力系統の所定区間において、負荷を有する第１需要家と、余剰買取契約対象の発電設備を有する第２需要家と、全量買取契約対象の前記発電設備を有する第３需要家と、のいずれかに分類される需要家を複数含む前記電力系統の状態を推定する情報処理プログラムであって、
コンピュータに、
前記電力系統の配電線のセンサからの測定値に基づいて、前記所定区間における有効電力及び無効電力を計算する第１計算部と、
前記複数の需要家のそれぞれが前記第１～第３需要家のいずれに該当するかを示す分類情報と、前記複数の需要家の電力量計の測定値と、に基づいて、前記第１～第３需要家のうち少なくともいずれかの有効電力を計算する第２計算部と、
前記第１及び第２需要家のそれぞれの負荷における有効電力の第１比率と、前記第２及び第３需要家のそれぞれの前記発電設備における有効電力の第２比率と、前記発電設備の力率と、前記所定区間における有効電力及び無効電力と、前記第１～第３需要家のうち少なくともいずれかの有効電力と、に基づいて、前記所定区間に含まれる負荷の有効電力及び無効電力と、前記所定区間に含まれる前記発電設備の有効電力と、を含む前記状態を推定する第１推定部と、
を実現させる、情報処理プログラム。
コンピュータに、
前記電力系統の配電線のセンサからの測定値に基づいて、前記所定区間における有効電力及び無効電力を計算する第１計算部と、
前記複数の需要家のそれぞれが前記第１～第３需要家のいずれに該当するかを示す分類情報と、前記複数の需要家の電力量計の測定値と、に基づいて、前記第１～第３需要家のうち少なくともいずれかの有効電力を計算する第２計算部と、
前記第１及び第２需要家のそれぞれの負荷における有効電力の第１比率と、前記第２及び第３需要家のそれぞれの前記発電設備における有効電力の第２比率と、前記発電設備の力率と、前記所定区間における有効電力及び無効電力と、前記第１～第３需要家のうち少なくともいずれかの有効電力と、に基づいて、前記所定区間に含まれる負荷の有効電力及び無効電力と、前記所定区間に含まれる前記発電設備の有効電力と、を含む前記状態を推定する第１推定部と、
を実現させる、情報処理プログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、情報処理装置、情報処理方法及び情報処理プログラムに関する。
【背景技術】
【0002】
電力系統において、例えば太陽光発電設備が存在する場合、配電線上に設けられたセンサの測定値は、需要家の負荷の有効電力と、太陽光発電設備の有効電力とが合算された値となる。そのため、センサの測定値からは、需要家の負荷の見かけ上の有効電力のみが把握される。
【0003】
しかしながら、電力系統の適切な運用計画や設備計画の観点から、需要家の負荷の有効電力等と、太陽光発電設備の有効電力とを個別に推定することが必要である。そこで、これらを個別に推定するための技術が知られている。
【0004】
例えば、特許文献１には、センサ内蔵自動開閉器の測定値と、日射量計の測定値とを用いて、需要家の負荷の有効電力等と、太陽光発電設備の有効電力とを推定する技術が開示されている。
【0005】
また、特許文献２には、配電線の送り出し点付近に設置されたセンサと、需要家のスマートメータの測定値とに基づいて、需要家の負荷の有効電力等と、太陽光発電設備の有効電力とを推定する技術が開示されている。
【先行技術文献】
【特許文献】
【0006】
【特許文献1】特開２０１２―０４４７４０号
【特許文献2】特開２０１６―０３９６５２号
【発明の概要】
【発明が解決しようとする課題】
【0007】
しかしながら、特許文献１に開示された技術においては、日射量計の導入や保守のための多大なコストが生じる。
【0008】
また、特許文献２に開示された技術は、全量買取契約対象である太陽光発電設備が存在することを前提とするため、このような発電設備が存在しない場合には適用することができない。
【0009】
本発明はこれらのような課題を鑑みてなされたものであり、需要家の負荷の有効電力等と、太陽光発電設備の有効電力とを精度良く推定することが可能な情報処理装置を提供することを目的とする。
【課題を解決するための手段】
【0010】
上記目的を達成するための一の発明は、電力系統の所定区間において、負荷を有する第１需要家と、余剰買取契約対象の発電設備を有する第２需要家と、全量買取契約対象の前記発電設備を有する第３需要家と、のいずれかに分類される需要家を複数含む前記電力系統の状態を推定する情報処理装置であって、前記電力系統の配電線のセンサからの測定値に基づいて、前記所定区間における有効電力及び無効電力を計算する第１計算部と、前記複数の需要家のそれぞれが前記第１～第３需要家のいずれに該当するかを示す分類情報と、前記複数の需要家の電力量計の測定値と、に基づいて、前記第１～第３需要家のうち少なくともいずれかの有効電力を計算する第２計算部と、前記第１及び第２需要家のそれぞれの負荷における有効電力の第１比率と、前記第２及び３需要家のそれぞれの前記発電設備における有効電力の第２比率と、前記発電設備の力率と、前記所定区間における有効電力及び無効電力と、前記第１～第３需要家のうち少なくともいずれかの有効電力と、に基づいて、前記所定区間に含まれる負荷の有効電力及び無効電力と、前記所定区間に含まれる前記発電設備の有効電力と、を含む前記状態を推定する第１推定部と、を備える、情報処理装置である。本発明の他の特徴については、本明細書の記載により明らかにする。
【発明の効果】
【0011】
本発明によれば、需要家の負荷の有効電力等と、発電設備の有効電力とを精度良く推定することが可能となる。
【図面の簡単な説明】
【0012】
【図1】実施形態の電力系統の一例と情報処理装置２を示す図である。
【図2】実施形態の情報処理装置２の機能ブロックを説明する図である。
【図3】比率を出力するまでの処理の流れを説明するフローチャートである。
【図4】実施形態の情報処理装置２が扱うデータの流れを説明する図である。
【図5】推定された状態ｘを出力するまでの処理の流れを説明するフローチャートである。
【図6】実施形態の情報処理装置２が扱うデータの流れを説明する図である。
【図7】最初の時間ステップの状態ｘを推定する処理の流れを説明するフローチャートである。
【図8】時間ステップごとの状態ｘを推定する処理の流れを説明するフローチャートである。
【図9】需要家の種別ごとの電力量計の測定値の合算の一例を示す図である。
【図10】推定された需要家の負荷の有効電力と、発電設備の有効電力と、これらの合算値の一例を示す図である。
【発明を実施するための形態】
【0013】
＝＝実施形態＝＝
＜＜電力系統１＞＞
図１は、後述する情報処理装置２が状態推定を行う電力系統１の一例を示す図である。電力系統１は、配電変電所１０と、配電線１１と、配電線１１の所定の位置に設置されたセンサ付開閉器ＳＷ１及びＳＷ２と、センサ付開閉器ＳＷ１，ＳＷ２の間の配電線１１の所定のノードＮｊ（ｊ＝１～ｎ）における負荷に設置された電力量計ＳＭｊ（ｊ＝１～ｎ）とを含む。
＜＜電力系統１＞＞
図１は、後述する情報処理装置２が状態推定を行う電力系統１の一例を示す図である。電力系統１は、配電変電所１０と、配電線１１と、配電線１１の所定の位置に設置されたセンサ付開閉器ＳＷ１及びＳＷ２と、センサ付開閉器ＳＷ１，ＳＷ２の間の配電線１１の所定のノードＮｊ（ｊ＝１～ｎ）における負荷に設置された電力量計ＳＭｊ（ｊ＝１～ｎ）とを含む。
【0014】
［配電変電所１０］
配電変電所１０は、送電線（不図示）から供給される電圧を変圧し、６．６ｋＶの電圧を配電線１１へと出力する。
配電変電所１０は、送電線（不図示）から供給される電圧を変圧し、６．６ｋＶの電圧を配電線１１へと出力する。
【0015】
［配電線１１］
配電線１１は、配電変電所１０を起点（送り出しノード）とし、配電変電所１０に放射状に接続されている。図１では、一の配電線１１のみが示されている。
配電線１１は、配電変電所１０を起点（送り出しノード）とし、配電変電所１０に放射状に接続されている。図１では、一の配電線１１のみが示されている。
【0016】
［ノード］
ノードは、配電線１１上に設けられている。本実施形態では、ノードは柱上変圧器（不図示）単位で管理されている集約単位である。配電変電所１０が配電線１１へ出力した電力は、ノードを介して需要家に供給される。なお、図１では、ｎ個のノードＮ１～ノードＮｎが設けられている。
ノードは、配電線１１上に設けられている。本実施形態では、ノードは柱上変圧器（不図示）単位で管理されている集約単位である。配電変電所１０が配電線１１へ出力した電力は、ノードを介して需要家に供給される。なお、図１では、ｎ個のノードＮ１～ノードＮｎが設けられている。
【0017】
［センサ付開閉器ＳＷ］
センサ付開閉器ＳＷ１，２のそれぞれは、設置点の少なくとも電圧、電流及び力率といった測定値を、定周期で測定可能なセンサを有する開閉器である。なお、センサ付き開閉器ＳＷ１と、センサ付き開閉器ＳＷ２とは同様であるため、以下、センサ付き開閉器ＳＷ１について説明する。本実施形態では、センサ付開閉器ＳＷ１が測定値を測定する周期は１分であるとする。
センサ付開閉器ＳＷ１，２のそれぞれは、設置点の少なくとも電圧、電流及び力率といった測定値を、定周期で測定可能なセンサを有する開閉器である。なお、センサ付き開閉器ＳＷ１と、センサ付き開閉器ＳＷ２とは同様であるため、以下、センサ付き開閉器ＳＷ１について説明する。本実施形態では、センサ付開閉器ＳＷ１が測定値を測定する周期は１分であるとする。
【0018】
また、本実施形態では、センサ付開閉器ＳＷ１は、１暦日においては、０：００から２３：５９まで、１分の測定周期で測定値を測定する。
【0019】
なお、必ずしもセンサ付開閉器ＳＷ１の測定周期は１分でなくてもよいし、全てのセンサ付開閉器ＳＷの測定周期が同じではなくてもよい。
【0020】
センサ付開閉器ＳＷ１によって測定された測定値は、測定値ＤＢ２２０（後述）に蓄積される。
【0021】
［需要家の分類］
需要家はそれぞれ、電力会社との電力買取契約の有無又は契約の内容に応じて、以下に説明する需要家Ａ、需要家Ｂ及び需要家Ｃのいずれかに分類される。つまり、電力系統１は、需要家Ａ需要家Ｂ及び需要家Ｃのいずれかに分類される需要家を複数含むことになる
需要家はそれぞれ、電力会社との電力買取契約の有無又は契約の内容に応じて、以下に説明する需要家Ａ、需要家Ｂ及び需要家Ｃのいずれかに分類される。つまり、電力系統１は、需要家Ａ需要家Ｂ及び需要家Ｃのいずれかに分類される需要家を複数含むことになる
【0022】
需要家Ａ（「第１需要家」に相当）は、電力会社と電力の買取契約を締結していない需要家である。
【0023】
需要家Ａは、電力を消費する負荷を有しているが、発電設備を有していない。負荷は、配電線１１から供給される電力を消費する。つまり、需要家Ａは、配電線１１から供給される電力を買電する。
【0024】
例えば、ノード１には、電力を消費する負荷を有する需要家Ａの設備として、集合住宅が接続されている。配電変電所１０が配電線１１へ出力した電力は、ノード１を介して集合住宅に供給される。
【0025】
需要家Ｂ（「第２需要家」に相当）は、電力会社と所謂余剰買取契約を締結している需要家である。
【0026】
需要家Ｂは、電力を消費する負荷と、余剰買取契約対象の発電設備とを有している。負荷は、配電線１１から供給される電力と、自己の発電設備が発電した電力の少なくとも一部とを消費する。発電設備としては、本実施形態では太陽光発電設備である。ただし、発電設備は太陽光発電設備に限られず、例えば、風力発電設備等でであってもよい
【0027】
つまり、需要家Ｂは、配電線１１から供給される電力を買電し、発電設備が発電した電力の少なくとも一部を消費する。更に、需要家Ｂは、発電設備が発電した電力のうち消費した電力を除く電力（余剰電力）を売電する。
【0028】
例えば、ノード２には、電力を消費する負荷と、発電設備とを有する需要家Ｂの設備として、太陽光発電設備が設置された一般住宅が接続されている。配電変電所１０が配電線１１へ出力した電力は、ノード２を介して一般住宅の負荷に供給される。更に、太陽光発電設備が発電した電力の少なくとも一部は、一般住宅の負荷に供給される。
【0029】
需要家Ｃ（「第３需要家」に相当）は、電力会社と所謂全量買取契約を締結している需要家である。
【0030】
需要家Ｃは、全量買取契約対象の発電設備を有している。需要家Ｃは、自己の発電設備が発電した電力の全量を売電する。
【0031】
例えば、ノード３には、発電設備を有する需要家Ｃの設備として、工場に設置された太陽光発電設備が接続されている。この需要家Ｃは、工場に設置された太陽光発電設備が発電した電力の全量を売電する。
【0032】
なお、この例では、工場は生産設備等の負荷を有するが、このような負荷は上述の需要家Ａに分類される。つまり、この例では、配電変電所１０が配電線１１へ出力した電力は、ノード３を介して工場の負荷に供給される。
【0033】
以下の説明では、電力系統１に含まれる需要家のうち、需要家Ａに分類される全ての需要家をまとめて、単に「需要家Ａ」と称する。需要家Ｂ及び需要家Ｃについても同様とする。
【0034】
［電力量計ＳＭ］
電力量計ＳＭは、電力系統１に含まれる需要家のそれぞれに設置される。電力量計ＳＭとしては、所謂スマートメータを用いることができる。
電力量計ＳＭは、電力系統１に含まれる需要家のそれぞれに設置される。電力量計ＳＭとしては、所謂スマートメータを用いることができる。
【0035】
需要家Ａ（電力買取契約を締結しておらず、負荷のみを有する需要家）については、電力を消費する負荷に対して１台の電力量計ＳＭが設置される。需要家Ａに設置された電力量計ＳＭは、負荷の使用電力量を、定周期で測定する。
【0036】
また、需要家Ｂ（余剰買取契約を締結しており、負荷及び発電設備を有する需要家）については、電力を消費する負荷と、発電設備とに対して一台の電力量計ＳＭが設置される。需要家Ｂに設置された電力量計ＳＭは、負荷の使用電力量と、自己の発電設備が発電した電力量との総計である電力量を、定周期で測定する。つまり、需要家Ｂに設置された電力量計ＳＭの測定値は、負荷の使用電力量と、自己の発電設備が発電した電力量とが混在したものである。
【0037】
また、需要家Ｃ（全量買取契約を締結しており、発電設備のみを有する需要家）については、発電設備に対して１台の電力量計ＳＭが設置される。需要家Ｃに設置された電力量計ＳＭは、自己の発電設備が発電した電力量を、定周期で測定する。
【0038】
本実施形態では、全ての電力量計ＳＭの測定周期は３０分であるとする。なお、電力量計ＳＭが測定した測定値を、測定周期で除した量は、測定周期の期間における平均の有効電力である。
【0039】
電力量計ＳＭは、１暦日においては、０：００から２３：３０まで、３０分の測定周期で測定値を測定する。
【0040】
なお、必ずしも電力量計ＳＭの測定周期は３０分でなくてもよいいし、全ての電力量計ＳＭの測定周期が同じではなくてもよい。
【0041】
電力量計ＳＭによって測定された測定値は、測定値ＤＢ２２１（後述）に蓄積される。
【0042】
＜＜情報処理装置２＞＞
情報処理装置２は、電力系統１の対象区間（「所定区間」に相当）において、電力系統１の状態を推定する装置である。なお、本実施形態において、対象区間は、図１に示す電力系統１のセンサ付開閉器ＳＷ１及びセンサ付開閉器ＳＷ２の間の区間であるとして説明する。
情報処理装置２は、電力系統１の対象区間（「所定区間」に相当）において、電力系統１の状態を推定する装置である。なお、本実施形態において、対象区間は、図１に示す電力系統１のセンサ付開閉器ＳＷ１及びセンサ付開閉器ＳＷ２の間の区間であるとして説明する。
【0043】
また、センサ付開閉器ＳＷ１及びＳＷ２の測定値のうち、情報処理装置２が扱う測定値は、センサ付開閉器ＳＷ１（配電変電所１０側）の測定値のみである。
【0044】
以下、情報処理装置２のハードウェア構成及び情報処理装置２の機能ブロックの順に説明する。
【0045】
＜情報処理装置２のハードウェア構成＞
図２は、本発明の一実施形態である情報処理装置２のハードウェア構成を説明する図である。情報処理装置２は、ＣＰＵ（Central Processing Unit）２００と、メモリ２０１と、通信装置２０２と、記憶装置２０３と、入力装置２０４と、出力装置２０５と、記録媒体読取装置２０６とを有するコンピュータである。
図２は、本発明の一実施形態である情報処理装置２のハードウェア構成を説明する図である。情報処理装置２は、ＣＰＵ（Central Processing Unit）２００と、メモリ２０１と、通信装置２０２と、記憶装置２０３と、入力装置２０４と、出力装置２０５と、記録媒体読取装置２０６とを有するコンピュータである。
【0046】
［ＣＰＵ２００］
ＣＰＵ２００は、メモリ２０１や記憶装置２０３に記憶された情報処理プログラムを実行することにより、情報処理装置２が有する様々な機能を実現する。
ＣＰＵ２００は、メモリ２０１や記憶装置２０３に記憶された情報処理プログラムを実行することにより、情報処理装置２が有する様々な機能を実現する。
【0047】
［メモリ２０１］
メモリ２０１は、例えばＲＡＭ（Random-Access Memory）等であり、様々なプログラムやデータ等の一時的な記憶領域として用いられる。
メモリ２０１は、例えばＲＡＭ（Random-Access Memory）等であり、様々なプログラムやデータ等の一時的な記憶領域として用いられる。
【0048】
［通信装置２０２］
通信装置２０２は、ネットワーク４を介して、他のコンピュータと各種プログラムやデータの受け渡しを行う。
通信装置２０２は、ネットワーク４を介して、他のコンピュータと各種プログラムやデータの受け渡しを行う。
【0049】
［記憶装置２０３］
記憶装置２０３は、ＣＰＵ２００によって、実行または処理される各種データを格納する非一時的な（例えば不揮発性の）記憶装置である。記憶装置２０３には、測定値ＤＢ２２０、測定値ＤＢ２２１、計算値ＤＢ２２２、計算値ＤＢ２２３、推定値ＤＢ２２４及び推定値ＤＢ２２５といったデータベースが格納されている（これらのデータベースに格納されているデータの詳細は後述する）。
記憶装置２０３は、ＣＰＵ２００によって、実行または処理される各種データを格納する非一時的な（例えば不揮発性の）記憶装置である。記憶装置２０３には、測定値ＤＢ２２０、測定値ＤＢ２２１、計算値ＤＢ２２２、計算値ＤＢ２２３、推定値ＤＢ２２４及び推定値ＤＢ２２５といったデータベースが格納されている（これらのデータベースに格納されているデータの詳細は後述する）。
【0050】
［入力装置２０４］
入力装置２０４は、ユーザによるコマンドやデータの入力を受け付ける装置であり、キーボード、タッチパネルディスプレイ上でのタッチ位置を検出するタッチセンサなどの入力インタフェースを含む。
入力装置２０４は、ユーザによるコマンドやデータの入力を受け付ける装置であり、キーボード、タッチパネルディスプレイ上でのタッチ位置を検出するタッチセンサなどの入力インタフェースを含む。
【0051】
［出力装置２０５］
出力装置２０５は、例えばディスプレイやプリンタなどの装置である。
出力装置２０５は、例えばディスプレイやプリンタなどの装置である。
【0052】
［記録媒体読取装置２０６］
記録媒体読取装置２０６は、ＳＤカードやＤＶＤ、ＣＤＲＯＭ等の記録媒体３に記録された情報処理プログラム等の様々なデータを読み取り、記憶装置２０３に格納する。
記録媒体読取装置２０６は、ＳＤカードやＤＶＤ、ＣＤＲＯＭ等の記録媒体３に記録された情報処理プログラム等の様々なデータを読み取り、記憶装置２０３に格納する。
【0053】
＜情報処理装置２の機能ブロック＞
情報処理装置２の機能ブロックを説明する。なお、ここでは概要を説明し、情報処理装置２が実行する詳細な処理については、後述するフローチャートを用いて詳細に説明する。
情報処理装置２の機能ブロックを説明する。なお、ここでは概要を説明し、情報処理装置２が実行する詳細な処理については、後述するフローチャートを用いて詳細に説明する。
【0054】
図２は、情報処理装置２の機能ブロックを示す図である。情報処理装置２は、取得部２１０と、取得部２１１と、計算部２１２，２１３と、推定部２１４，２１５と、出力部２１６とを含む。
【0055】
［取得部２１０］
取得部２１０は、記憶装置２０３の測定値ＤＢ２２０を参照し、センサ付開閉器ＳＷ１の所定期間Ｔの測定値を取得する。なお、ここでの所定期間Ｔとは、負荷における有効電力の比率及び発電設備における有効電力の比率（詳細は後述）を推定するための過去の期間ＴＡと、電力系統１の状態推定の対象となる期間（以下、「対象期間ＴＢ」と称する）とを含む。これらの期間の詳細は後述する。
取得部２１０は、記憶装置２０３の測定値ＤＢ２２０を参照し、センサ付開閉器ＳＷ１の所定期間Ｔの測定値を取得する。なお、ここでの所定期間Ｔとは、負荷における有効電力の比率及び発電設備における有効電力の比率（詳細は後述）を推定するための過去の期間ＴＡと、電力系統１の状態推定の対象となる期間（以下、「対象期間ＴＢ」と称する）とを含む。これらの期間の詳細は後述する。
【0056】
［取得部２１１］
取得部２１１は、記憶装置２０３の測定値ＤＢ２２１を参照し、電力量計ＳＭの所定期間の測定値を取得する。なお、ここでの所定期間Ｔとは、上述の取得部２１０についての所定期間Ｔと同様であり、詳細は後述する。
取得部２１１は、記憶装置２０３の測定値ＤＢ２２１を参照し、電力量計ＳＭの所定期間の測定値を取得する。なお、ここでの所定期間Ｔとは、上述の取得部２１０についての所定期間Ｔと同様であり、詳細は後述する。
【0057】
［計算部２１２］
計算部２１２（「第１計算部」に相当）は、取得部２１０が取得した、電力系統１の配電線１１のセンサ付開閉器ＳＷ１からの測定値に基づいて、対象区間における有効電力Ｐ及び無効電力Ｑを計算する。なお、ここで計算された有効電力Ｐは、電力系統１に含まれる需要家の負荷の有効電力と、発電設備の有効電力とが混在した値となる。
計算部２１２（「第１計算部」に相当）は、取得部２１０が取得した、電力系統１の配電線１１のセンサ付開閉器ＳＷ１からの測定値に基づいて、対象区間における有効電力Ｐ及び無効電力Ｑを計算する。なお、ここで計算された有効電力Ｐは、電力系統１に含まれる需要家の負荷の有効電力と、発電設備の有効電力とが混在した値となる。
【0058】
【0059】
ここで、分類情報とは、電力系統に含まれる複数の需要家のそれぞれが需要家Ａ～Ｃのいずれに該当するかを示す情報である。
【0060】
【0061】
［推定部２１４］
推定部２１４（「第２推定部）に相当）は、負荷における有効電力の比率Ｒ１（「第１比率」に相当）を推定する。負荷における有効電力の比率Ｒ１とは、需要家Ａ（電力買取契約を締結しておらず、負荷のみを有する需要家）及び需要家Ｂ（余剰買取契約を締結しており、負荷及び発電設備を有する需要家）のそれぞれの負荷における有効電力の比率である。
推定部２１４（「第２推定部）に相当）は、負荷における有効電力の比率Ｒ１（「第１比率」に相当）を推定する。負荷における有効電力の比率Ｒ１とは、需要家Ａ（電力買取契約を締結しておらず、負荷のみを有する需要家）及び需要家Ｂ（余剰買取契約を締結しており、負荷及び発電設備を有する需要家）のそれぞれの負荷における有効電力の比率である。
【0062】
推定部２１４は、更に、発電設備における有効電力の比率Ｒ２（「第２比率」に相当）を推定する。発電設備における有効電力の比率Ｒ２とは、需要家Ｂ（余剰買取契約を締結しており、負荷及び発電設備を有する需要家）及び需要家Ｃ（全量買取契約を締結しており、発電設備のみを有する需要家）のそれぞれの発電設備における有効電力の比率である。
【0063】
なお、以下の説明では、負荷における有効電力の比率Ｒ１と、発電設備における有効電力の比率Ｒ２とを総称して、「比率Ｒ」と称する場合がある。
【0064】
［推定部２１５］
推定部２１５（「第１推定部に相当」）は、対象区間における電力系統１の状態ｘを推定する。ここで、状態ｘとは、対象区間に含まれる負荷の有効電力ＰＬ及び無効電力ＱＬと、対象区間に含まれる発電設備の有効電力ＰＰＶと、を含む。なお、ここでの有効電力ＰＬ及び無効電力ＱＬと、発電設備の有効電力ＰＰＶとは、対象区間に含まれる全ての需要家についての合算値である。
推定部２１５（「第１推定部に相当」）は、対象区間における電力系統１の状態ｘを推定する。ここで、状態ｘとは、対象区間に含まれる負荷の有効電力ＰＬ及び無効電力ＱＬと、対象区間に含まれる発電設備の有効電力ＰＰＶと、を含む。なお、ここでの有効電力ＰＬ及び無効電力ＱＬと、発電設備の有効電力ＰＰＶとは、対象区間に含まれる全ての需要家についての合算値である。
【0065】
推定部２１５は、計算部２１３によって計算された負荷における有効電力の比率Ｒ１と、発電設備における有効電力の比率Ｒ２と、発電設備の力率と、需要家Ａ～Ｃの有効電力と、に基づいて、対象区間における状態ｘを推定する。
【0066】
［出力部２１６］
出力部２１６は、推定部２１４によって推定された比率Ｒを出力する。出力部２１６は、更に、推定部２１５によって推定された電力系統１の対象区間における状態ｘを出力する。
出力部２１６は、推定部２１４によって推定された比率Ｒを出力する。出力部２１６は、更に、推定部２１５によって推定された電力系統１の対象区間における状態ｘを出力する。
【0067】
＜＜比率Ｒを出力するまでの処理＞＞
図３及び図４を用いて、情報処理装置２が比率Ｒを出力するまでの処理の流れについて説明する。図３は、情報処理装置２が比率Ｒを出力するまでの処理の流れを説明するフローチャートである。図４は、情報処理装置２が扱うデータの流れを説明する図である。図３の処理は、ステップＳ１１～ステップＳ１５を含む。
図３及び図４を用いて、情報処理装置２が比率Ｒを出力するまでの処理の流れについて説明する。図３は、情報処理装置２が比率Ｒを出力するまでの処理の流れを説明するフローチャートである。図４は、情報処理装置２が扱うデータの流れを説明する図である。図３の処理は、ステップＳ１１～ステップＳ１５を含む。
【0068】
先ず、図３のステップＳ１１において、取得部２１１は、図４に示すように、測定値ＤＢ２２１から、過去の所定期間における電力量計ＳＭの測定値を取得する。ここで、取得部２１１は、対象区間に含まれる全ての需要家の電力量計ＳＭの測定値を取得する。
【0069】
次いで、ステップＳ１２において、計算部２１３は、ステップＳ１１において取得部２１１が取得した電力量計ＳＭの測定値を用いて、需要家Ａ～Ｃのそれぞれについて、電力量計ＳＭの測定値の合算値を計算する（図３のＳ１２）。
【0070】
ここで、計算部２１３は、過去の所定期間の各時間ステップに測定された全ての電力量計ＳＭの測定値を用いて、需要家Ａ～Ｃのそれぞれについて、電力量計ＳＭの測定値の合算値を計算する。
【0071】
次いで、ステップＳ１３において、計算部２１３は、需要家Ａ～Ｃのそれぞれの有効電力を計算する（図３のＳ１３）。
【0072】
ここで、計算部２１３は、ステップＳ１２において計算された需要家Ａ～Ｃのそれぞれについての合算値を、電力量計ＳＭの測定周期（３０分）で除することにより、需要家Ａ～Ｃのそれぞれの有効電力を計算する。このようにして計算された有効電力は、測定周期の期間における平均の有効電力である。
【0073】
計算部２１３が計算した有効電力は、図４に示す計算値ＤＢ２２３に蓄積される。
【0074】
次いで、ステップＳ１４において、推定部２１４は、負荷における有効電力の比率Ｒ１と、発電設備における有効電力の比率Ｒ２を推定する（図３のＳ１４）。
【0075】
ここで、負荷における有効電力の比率Ｒ１とは、対象区間に含まれる負荷の有効電力のうち、需要家Ａ（電力買取契約を締結しておらず、負荷のみを有する需要家）の負荷の有効電力の割合を示す値であるαＬ（「第１の値」に相当）と、需要家Ｂ（余剰買取契約を締結しており、負荷及び発電設備を有する需要家）の負荷の有効電力の割合を示す値であるαＬ＿ＰＶ（「第２の値」に相当）である。
【0076】
また、発電設備における有効電力の比率Ｒ２とは、対象区間に含まれる発電設備の有効電力のうち、需要家Ｂ（余剰買取契約を締結しており、負荷及び発電設備を有する需要家）の発電設備の有効電力の割合を示す値であるβＬ＿ＰＶ（「第３の値」に相当）と、需要家Ｃ（全量買取契約を締結しており、発電設備のみを有する需要家）の発電設備の有効電力の割合を示す値であるβＰＶ（「第４の値」に相当）である。
【0077】
これらの値の間には、以下の数１～数４の関係が成り立つ。
【数1】
【数2】
【数3】
【数4】
【0078】
このステップＳ１５において、推定部２１４は、所定の目的関数（「第２目的関数」に相当）を設定し、目的関数を最小化することにより、比率Ｒを推定する。
【0079】
目的関数の詳細な説明の前に、以下の数５～数７を用いて、目的関数に用いる変数について説明する。
【数5】
【数6】
【数7】
【0080】
【0081】
数５の右辺において、αＬは、需要家Ａの負荷の有効電力の割合を示す値であり、未知量である。また、本実施形態では、αＬはスカラー量であると仮定する。
【0082】
【0083】
【0084】
【0085】
数６の右辺において、βＰＶは、需要家Ｃの発電設備の有効電力の割合を示す値であり、未知量である。また、本実施形態では、βＰＶはスカラー量であると仮定する。
【0086】
【0087】
【0088】
本実施形態で用いる目的関数と、目的関数の最小化条件を以下の数８に示す。
【数8】
【0089】
【0090】
【0091】
【0092】
【0093】
また、このとき、推定部２１４は、以下の数９及び数１０に示す拘束条件を課して、目的関数を最小化する。
【数9】
【数10】
【0094】
【0095】
次いで、ステップＳ１５において、出力部２１６は、ステップＳ１４において推定部２１４が推定した比率Ｒ（αＬ、αＬ＿ＰＶ、βＬ＿ＰＶ及びβＰＶ）を、図４に示す推定値ＤＢ２２４に出力する。
【0096】
以上で、比率Ｒを出力するまでの処理が終了する。以上の処理によって出力された比率Ｒは、次に説明する状態ｘを推定する処理において用いられる。
【0097】
＜＜推定した状態ｘを出力するまでの処理＞＞
情報処理装置２が、推定した状態ｘを出力するまでの処理について説明する。以下では、状態推定の対象となる期間（対象期間）を、状態推定の対象となる日の０：００から２３：５９までとして説明する。
情報処理装置２が、推定した状態ｘを出力するまでの処理について説明する。以下では、状態推定の対象となる期間（対象期間）を、状態推定の対象となる日の０：００から２３：５９までとして説明する。
【0098】
情報処理装置２は、対象期間を有限の時間ステップに区切り、時間ステップごとの状態ｘを推定することにより、対象期間の状態の時系列ｘ（ｔ）を推定する。
【0099】
前述のように、センサ付開閉器ＳＷ１は、１暦日の０：００から、１分の間隔で測定値を測定する。以下の説明における時間ステップは、センサ付開閉器ＳＷ１が測定値を測定する時間とする。
【0100】
つまり、最初の時間ステップは、推定対象日の０：００である。２回目の時間ステップは、その１分後の０：０１である。最後の時間ステップは、１４４０回目の時間ステップであり、同日の２３：５９である。
【0101】
図５及び図６を用いて、情報処理装置２が、推定した状態ｘを出力するまでの処理の流れについて説明する。図５は、情報処理装置２が、推定した状態ｘを出力するまでの処理の流れを説明するフローチャートである。図６は、情報処理装置２が扱うデータの流れを説明する図である。この処理は、ステップＳ２１～ステップＳ２５を含む。
【0102】
先ず、図５のステップＳ２１において、取得部２１０は、測定値ＤＢ２２０から、対象期間におけるセンサ付開閉器ＳＷ１の測定値を取得する。
【0103】
本実施形態において、センサ付開閉器ＳＷ１の測定値とは、センサ付開閉器ＳＷ１の設置点における電圧、電流及び力率である。
【0104】
つまり、取得部２１０は、推定対象日の、０：００から２３：５９まで１分間隔で測定された、電圧、電流及び力率を取得する。
【0105】
次いで、ステップＳ２２において、取得部２１１は、測定値ＤＢ２２１から、対象期間における電力量計ＳＭの測定値を取得する。ここで、取得部２１１は、対象区間に含まれる全ての需要家の電力量計ＳＭの測定値を取得する。
【0106】
つまり、取得部２１１は、推定対象日の、０：００から２３：３０まで３０分間隔で測定された、全ての需要家の電力量計ＳＭの測定値を取得する。
【0107】
次いで、ステップＳ２３は、対象期間における最初の時間ステップ（つまり、０：００）の状態ｘ（０）を推定する処理である。なお、最初の時間ステップにおいては、センサ付開閉器ＳＷ１及び電力量計ＳＭは共に、測定値を測定している。
【0108】
図７は、ステップＳ２３の詳細を説明するフローチャートである。ステップＳ２３は、ステップＳ２３１～ステップＳ２３４を含む。
【0109】
先ず、図７のステップＳ２３１において、計算部２１２は、ステップＳ２１において取得部２１０が取得したセンサ付開閉器ＳＷ１の測定値に基づいて、対象区間における有効電力及び無効電力を計算する（図６のＳ２３１）。
【0110】
ここで、計算部２１２は、最初の時間ステップ（つまり、０：００）に測定された電圧、電流及び力率に基づいて、最初の時間ステップにおける有効電力及び無効電力を計算する。
【0111】
計算部２１２が計算した有効電力Ｐ（ｔ）及び無効電力Ｑ（ｔ）は、図６に示す計算値ＤＢ２２２に蓄積される。
【0112】
次いで、ステップＳ２３２において、計算部２１３は、ステップＳ２２において取得部２１１が取得した電力量計ＳＭの測定値を用いて、需要家Ａ～Ｃのそれぞれについて、電力量計ＳＭの測定値の合算値を計算する（図６のＳ２３２）。
【0113】
ここで、計算部２１３は、最初の時間ステップ（つまり、０：００）に測定された全ての電力量計ＳＭの測定値を用いて、需要家Ａ～Ｃのそれぞれについて、電力量計ＳＭの測定値の合算値を計算する。
【0114】
次いで、ステップＳ２３３において、計算部２１３は、需要家Ａ～Ｃのそれぞれの有効電力を計算する（図６のＳ２３３）。
【0115】
ここで、計算部２１３は、ステップＳ２３２において計算された需要家Ａ～Ｃのそれぞれについての合算値を、電力量計ＳＭの測定周期（３０分）で除することにより、需要家Ａ～Ｃのそれぞれの有効電力を計算する。このようにして計算された有効電力は、測定周期の期間における平均の有効電力である。
【0116】
計算部２１３が計算した平均の有効電力は、図６に示す計算値ＤＢ２２３に蓄積される。
【0117】
次いで、ステップＳ２３４において、推定部２１５は、対象期間における最初の時間ステップの状態ｘ（０）を推定する（図６のＳ２３４）。
【0118】
このステップＳ２３４において、推定部２１５は、所定の目的関数（「第１目的関数」に相当）を設定し、目的関数を最小化することにより、最初の時間ステップの状態ｘ（０）を推定する。
【0119】
目的関数の詳細な説明の前に、下記の数１１～数１５を用いて、推定部２１５が目的関数を設定する際に用いる変数について説明する。
【0120】
先ず、数１１は、センサ付開閉器ＳＷ１の測定値に基づいて計算された有効電力及び無効電力と、電力量計ＳＭの測定値に基づいて計算された需要家の有効電力とを示している。
【数11】
【0121】
【0122】
【0123】
次いで、数１２は、対象区間の状態と、センサ付開閉器ＳＷ１及び電力量計ＳＭの測定値との関係を示す式である。
【数12】
【0124】
【数13】
【数14】
【数15】
【0125】
数１２において、右辺の行列Ａは、数１４に示すように、ステップＳ１５（図４）で計算された比率Ｒ（αＬ、αＬ＿ＰＶ、βＬ＿ＰＶ及びβＰＶ）と、発電設備の力率角θＰＶによって定められる行列である。力率角θＰＶは、発電設備が力率を一定として運転されていることから、既知量である。従って、行列Ａは既知である。
【0126】
比率Ｒ（αＬ、αＬ＿ＰＶ、βＬ＿ＰＶ及びβＰＶ）は、図６に示す推定値ＤＢ２２４から取得される。
【0127】
【0128】
【0129】
本実施形態で用いる目的関数（「第１目的関数」に相当）と、目的関数の最小化条件を以下の数１６に、最小化条件を数１７に示す。
【数16】
【数17】
【0130】
【0131】
【0132】
【0133】
【0134】
【0135】
本実施形態の目的関数は、上述の差（数１３の各成分）に基づく関数である。
【0136】
推定部２１５は、数１６に示す目的関数を設定し、状態ｘ（ｔ）（数１５）を決定変数（「第１決定変数」に相当）とした際に目的関数を最小化することにより、状態ｘ（ｔ）を推定する（数１７）。
【0137】
ステップＳ２３４においては、推定部２１５は、時間ステップｔが０の場合の数１６に示す目的関数を最小化することにより、最初の時間ステップにおける状態ｘ（０）を推定する。
【0138】
推定部２１５が推定した状態ｘ（０）は、図６に示す推定値ＤＢ２２５に蓄積される。
【0139】
次いで、図５に戻り、ステップＳ２４は、時間ステップごとの状態ｘ（ｔ）を推定する処理である。情報処理装置２は、ステップＳ２４の処理を２回目の時間ステップから最後の時間ステップまで繰り返すことにより、状態ｘ（ｔ）の時系列を推定する。
【0140】
図８は、ステップＳ２４の詳細を説明するフローチャートである。以下の説明では、ｔ－１回目の時間ステップまでの状態ｘの推定が完了しており、ｔ回目の時間ステップの状態ｘ（ｔ）を推定する場合を説明する。
【0141】
ステップＳ２４は、ステップＳ２４１～ステップＳ２４６を含む。
【0142】
先ず、ステップＳ２４１において、計算部２１２は、センサ付開閉器ＳＷ１の測定値を更新する。つまり、計算部２１２は、以下の計算で扱う測定値の対象を、ｔ回目の時間ステップにおけるセンサ付開閉器ＳＷ１の測定値に更新する。
【0143】
次いで、ステップＳ２４２において、計算部２１２は、ステップＳ２４１において計算部２１２が更新したセンサ付開閉器ＳＷ１の測定値に基づいて、対象区間における有効電力及び無効電力を計算する。
【0144】
ここで、計算部２１２は、ｔ回目の時間ステップに測定された電圧、電流及び力率に基づいて、ｔ回目の時間ステップにおける有効電力及び無効電力を計算する。
【0145】
次いで、ステップＳ２４３において、計算部２１３は、ｔ回目の時間ステップにおいて、電力量計ＳＭの測定値の更新が有るか否かを判定する。
【0146】
前述のように、電力量計ＳＭの測定周期は３０分である。そのため、例えば、０：００（最初の時間ステップ）の次に電力量計ＳＭの測定値が更新される時刻は、０：３０（３１回目の時間ステップ）である。
【0147】
ステップＳ２４３において、計算部２１３が電力量計ＳＭの測定値が更新されたと判断した場合（Ｓ２４３：Ｙ）、計算部２１３は、以下の計算で扱う測定値の対象を、ｔ回目の時間ステップにおける電力量計ＳＭの測定値に更新し、ステップＳ２４４に進む。
【0148】
次いで、ステップＳ２４４において、計算部２１３は、ステップＳ２４３において更新した電力量計ＳＭの測定値を用いて、需要家Ａ～Ｃのそれぞれについて、電力量計ＳＭの測定値の合算値を計算する。
【0149】
ここで、計算部２１３は、ｔ回目の時間ステップに測定された全ての電力量計ＳＭの測定値を用いて、需要家Ａ～Ｃのそれぞれについて、電力量計ＳＭの測定値の合算値を計算する。
【0150】
次いで、ステップＳ２４５において、計算部２１３は、需要家Ａ～Ｃのそれぞれの有効電力を計算する。
【0151】
ここで、計算部２１３は、ステップＳ２３５において計算された需要家Ａ～Ｃのそれぞれについての合算値を、電力量計ＳＭの測定周期（３０分）で除することにより、需要家Ａ～Ｃのそれぞれの有効電力を計算する。このようにして計算された有効電力は、測定周期の期間における平均の有効電力である。
【0152】
図９は、計算部２１３が計算した需要家Ａ～Ｃの各時間ステップの有効電力の一例である。
【0153】
ステップＳ２４３において、計算部２１３が電力量計ＳＭの測定値が更新されていないと判断した場合（Ｓ２４３：Ｎ）、計算部２１３は、以下の計算で扱う測定値の対象を、ｔ－１回目の時間ステップにおける電力量計ＳＭの測定値のまま維持し、ステップＳ２３６に進む。
【0154】
ステップＳ２４６において、推定部２１５は、時間ステップごとの状態ｘ（ｔ）を推定する。以下では、推定部２１５がｔ回目の時間ステップの状態ｘ（ｔ）を推定する処理について説明する。
【0155】
ｔ回目の時間ステップの状態ｘ（ｔ）は、以下の数１８～数２０により推定することができる。
【数18】
【数19】
【数20】
【0156】
数１８において、左辺は、対象区間の状態ｘ（ｔ）であり、各成分は数１４に示されている。
【0157】
また、数１８において、右辺第１項及び第２項の状態ｘ（ｔ－１）は、既に推定されたｔ－１回目の時間ステップの状態である。
【0158】
右辺第２項の行列Ｋは、ゲイン行列であり、数１９に具体的に示されている。数１９のＷは、状態ｘの共分散行列であり、数１９においてはｔ－１回目の時間ステップにおける共分散行列が用いられている。
【0159】
共分散行列Ｗは、数１９により、時間ステップごとに更新される。数１９により、ｔ－１回目の時間ステップにおける共分散行列Ｗ（ｔ－１）に基づいて、ｔ回目の時間ステップにおける共分散行列Ｗ（ｔ）を更新することができる。
【0160】
また、数１８において、右辺第２項の行列Ａは、数１４に示された行列であり、既知の行列である。
【0161】
推定部２１５は、ステップＳ２３６の処理を繰り返すことにより、対象区間の状態ｘ（ｔ）（需要家の負荷の有効電力ＰＬ及び無効電力ＱＬと、発電設備の有効電力ＰＰＶ）を推定する。
【0162】
図１０は、推定された需要家の負荷の有効電力ＰＬと、発電設備の有効電力ＰＰＶと、これらの合算値を示す図である。ここでの合算値とは、数９及び数１０により、図９に示した需要家Ａ～Ｃの有効電力の合算値に等しい。
【0163】
前述のように、本実施形態では、発電設備が太陽光発電設備である態様を示している。そのため、図１０からわかるように、発電設備の有効電力ＰＰＶ（つまり、太陽光発電設備が発電した電力）は、日中の時間帯が大きく、正午（１２：００）付近で極大となる。また、夜間の時間帯は小さく、５：００以前及び１８：００以降は殆ど発電しない。
【0164】
次いで、ステップＳ２５において、出力部２１６は、ステップＳ２４において推定部２１４が推定した状態ｘ（ｔ）を（需要家の負荷の有効電力ＰＬ及び発電設備の有効電力ＰＰＶ）、推定値ＤＢ２２５に出力する。以上で、推定した状態ｘ（ｔ）を出力するまでの処理が終了する。
【0165】
以上のような処理によれば、低コストで、需要家の電力買取契約の形態に依らずに、需要家の負荷の有効電力等と、発電設備の有効電力とを推定することが可能となる。
【0166】
＝＝変形例１＝＝
実施形態においては、電力系統１に需要家Ａ～Ｃのそれぞれの電力量計ＳＭの測定値が得られる態様を示した。しかし、このような場合に限らず、需要家Ａ～Ｃのうちいずれか一の需要家の電力量計ＳＭの測定値が得られない場合であっても、対象区間の状態ｘ（ｔ）を推定することができる。
実施形態においては、電力系統１に需要家Ａ～Ｃのそれぞれの電力量計ＳＭの測定値が得られる態様を示した。しかし、このような場合に限らず、需要家Ａ～Ｃのうちいずれか一の需要家の電力量計ＳＭの測定値が得られない場合であっても、対象区間の状態ｘ（ｔ）を推定することができる。
【0167】
本先ず、一例として、需要家Ａ及びＢの電力量計ＳＭの測定値が得られ、需要家Ｃの電力量計ＳＭの測定値が得られない場合を説明する。
【0168】
この場合、実施形態で用いた数１１、数１３及び数１４に代えて、それぞれ以下の数２１、数２２及び数２３を用いればよい。
【数21】
【数22】
【数23】
【0169】
数２１は、数１１の右辺の３行目を省いたものに相当する。また、数２２は、数１３の右辺の３行目を省いたものに相当する。また、数２３は、数１４の右辺の３行目を省いたものに相当する。
【0170】
以上、需要家Ｃの電力量計ＳＭの測定値が得られない場合について説明したが、需要家Ａ又は需要家Ｂの電力量計ＳＭの測定値が得られない場合についても拡張することができる。
【0171】
需要家Ａの電力量計ＳＭの測定値が得られない場合、実施形態で用いた数１１、数１３及び数１４に代えて、それぞれの１行目を省いたものを用いればよい。
【0172】
需要家Ｂの電力量計ＳＭの測定値が得られない場合、実施形態で用いた数１１、数１３及び数１４に代えて、それぞれの２行目を省いたものを用いればよい。
【0173】
つまり、本変形例においては、計算部２１３は、需要家Ａ～Ｃのうち少なくともいずれか２つの有効電力を計算する。
【0174】
そして、推定部２１５は、比率Ｒ１（αＬ及びαＬ＿ＰＶ）と、比率Ｒ２（βＬ＿ＰＶ及びβＰＶ）と、力率と、対象区間における有効電力及び無効電力と、需要家Ａ～Ｃのうち少なくともいずれか２つの有効電力に基づいて、状態ｘ（ｔ）を推定する。
【0175】
＝＝変形例２＝＝
更に、需要家Ａ～Ｃのうち、いずれか２つの需要家の電力量計ＳＭの測定値が得られない場合（つまり、いずれか１つの需要家の電力量計ＳＭの測定値のみが得られる場合）であっても、対象区間の状態ｘ（ｔ）を推定することができる。
更に、需要家Ａ～Ｃのうち、いずれか２つの需要家の電力量計ＳＭの測定値が得られない場合（つまり、いずれか１つの需要家の電力量計ＳＭの測定値のみが得られる場合）であっても、対象区間の状態ｘ（ｔ）を推定することができる。
【0176】
需要家Ａの電力量計ＳＭの測定値のみが得られる場合、実施形態で用いた数１１、数１３及び数１４に代えて、それぞれの２行目及び３行目を省いたものを用いればよい。
【0177】
需要家Ｂの電力量計ＳＭの測定値のみが得られる場合、実施形態で用いた数１１、数１３及び数１４に代えて、それぞれの１行目及び３行目を省いたものを用いればよい。
【0178】
需要家Ｃの電力量計ＳＭの測定値のみが得られる場合、実施形態で用いた数１１、数１３及び数１４に代えて、それぞれの１行目及び２行目を省いたものを用いればよい。
【0179】
つまり、本変形例において、計算部２１３は、需要家Ａ～Ｃのうち少なくともいずれかの有効電力を計算する。
【0180】
そして、推定部２１５は、比率Ｒ１（αＬ及びαＬ＿ＰＶ）と、比率Ｒ２（βＬ＿ＰＶ及びβＰＶ）と、力率と、対象区間における有効電力及び無効電力と、需要家Ａ～Ｃのうち少なくともいずれかの有効電力に基づいて、対象区間の状態ｘ（ｔ）を推定することができる。
【0181】
＝＝変形例３＝＝
実施形態においては、負荷における有効電力の比率Ｒ１（αＬ及びαＬ＿ＰＶ）は一定値であり、時々刻々変化しないと仮定する態様を示した。
実施形態においては、負荷における有効電力の比率Ｒ１（αＬ及びαＬ＿ＰＶ）は一定値であり、時々刻々変化しないと仮定する態様を示した。
【0182】
しかしながら、実際には、負荷における有効電力の比率Ｒ１（αＬ及びαＬ＿ＰＶ）が時々刻々変化し、αＬ及びαＬ＿ＰＶを一定値と仮定すると十分な精度で状態ｘ（ｔ）を推定することが困難となる場合がある。
【0183】
そのような場合、需要家Ａと、需要家Ｂとを区別せずに取り扱うこととしてもよい。具体的には、実施形態で用いた数１１、数１３及び数１４に代えて、それぞれ以下の数２４、数２５及び数２６を用いればよい。
【数24】
【数25】
【数26】
【0184】
数２４は、数１１の１行目に２行目を加算し、２行目を省いたものに相当する。数２５は、数１３の１行目に２行目を加算し、２行目を省いたものに相当する。数２６は、数１４の１行目に２行目を加算し、２行目を省いたものに相当する。また、数２６の２つ目の等式において数１の関係を用いている。
【0185】
これによって、数２６の２つ目の等式からわかるように、αＬ及びαＬ＿ＰＶが消去されるため、これらの時間依存性を考慮する必要がない。
【0186】
つまり、本変形例において、計算部２１３は、需要家Ａ及び需要家Ｂの有効電力を合算した有効電力と、需要家Ｃの有効電力とを計算する。
【0187】
そして、推定部２１５は、比率Ｒ２（βＬ＿ＰＶ及びβＰＶ）と、力率と、対象区間における有効電力及び無効電力と、需要家Ａ及び需要家Ｂの有効電力を合算した有効電力と、需要家Ｃの有効電力とに基づいて、対象区間の状態ｘ（ｔ）を推定することができる。
【0188】
これによって、αＬ及びαＬ＿ＰＶが時々刻々変化することが想定される場合であっても、精度良く状態ｘ（ｔ）を推定することができる。
【0189】
＝＝変形例４＝＝
実施形態においては、発電設備における有効電力の比率Ｒ２（βＬ＿ＰＶ及びβＰＶ）は一定値であり、時々刻々変化しないと仮定する態様を示した。
実施形態においては、発電設備における有効電力の比率Ｒ２（βＬ＿ＰＶ及びβＰＶ）は一定値であり、時々刻々変化しないと仮定する態様を示した。
【0190】
しかしながら、実際には、発電設備における有効電力の比率Ｒ２（βＬ＿ＰＶ及びβＰＶ）が時々刻々変化し、βＬ＿ＰＶ及びβＰＶを一定値と仮定すると十分な精度で状態ｘ（ｔ）を推定することが困難となる場合がある。
【0191】
そのような場合、需要家Ｂと、需要家Ｃとを区別せずに取り扱うこととしてもよい。具体的には、実施形態で用いた数１１、数１３及び数１４に代えて、それぞれ以下の数２７、数２８及び数２９を用いればよい。
【数27】
【数28】
【数29】
【0192】
数２７は、数１１の２行目に３行目を加算し、３行目を省いたものに相当する。数２５は、数１３の２行目に３行目を加算し、３行目を省いたものに相当する。数２６は、数１４の２行目に３行目を加算し、３行目を省いたものに相当する。また、数２６の２つ目の等式において数３の関係を用いている。
【0193】
ここで、数２６の２つ目の等式からわかるように、βＬ＿ＰＶ及びβＰＶが消去されるため、これらの時間依存性を考慮する必要がない。
【0194】
つまり、本変形例において、計算部２１３は、需要家Ａの有効電力と、需要家Ｂ及びＣの有効電力を合算した有効電力とを計算する。
【0195】
そして、推定部２１５は、比率Ｒ１（αＬ及びαＬ＿ＰＶ）と、力率と、対象区間における有効電力及び無効電力と、需要家Ａ及び需要家Ｂの有効電力を合算した有効電力と、需要家Ａの有効電力と、需要家Ｂ及びＣの有効電力を合算した有効電力とに基づいて、状態ｘ（ｔ）を推定することができる。
【0196】
これによって、βＬ＿ＰＶ及びβＰＶが時々刻々変化することが想定される場合であっても、精度良く対象区間の状態ｘ（ｔ）を推定することができる。
【0197】
＝＝変形例５＝＝
実施形態においては、負荷における有効電力の比率Ｒ１（αＬ及びαＬ＿ＰＶ）と、発電設備における有効電力の比率Ｒ２（βＬ＿ＰＶ及びβＰＶ）とは、それぞれ一定値であり、時々刻々変化しないと仮定する態様を示した。
実施形態においては、負荷における有効電力の比率Ｒ１（αＬ及びαＬ＿ＰＶ）と、発電設備における有効電力の比率Ｒ２（βＬ＿ＰＶ及びβＰＶ）とは、それぞれ一定値であり、時々刻々変化しないと仮定する態様を示した。
【0198】
しかしながら、実際には、αＬ、αＬ＿ＰＶ、βＬ＿ＰＶ及びβＰＶが時々刻々変化し、これらを一定値と仮定すると十分な精度で状態ｘ（ｔ）を推定することが困難となる場合がある。
【0199】
そのような場合、αＬ、αＬ＿ＰＶ、βＬ＿ＰＶ及びβＰＶに時間依存性を付与することとしてもよい。具体的には、推定部２１４は、実施形態で用いた数８の目的関数の最小化において、時間依存性を有する、αＬ、αＬ＿ＰＶ、βＬ＿ＰＶ及びβＰＶを、決定変数とすればよい。
【0200】
これによって、αＬ、αＬ＿ＰＶ、βＬ＿ＰＶ及びβＰＶが時々刻と変化することが想定される場合であっても、精度良く状態ｘ（ｔ）を推定することができる。
【0201】
＝＝まとめ＝＝
実施形態の情報処理装置２は、電力系統１の所定区間において、負荷を有する需要家Ａと、余剰買取契約対象の発電設備を有する需要家Ｂと、全量買取契約対象の発電設備を有する需要家Ｃと、のいずれかに分類される需要家を複数含む電力系統１の状態ｘを推定する情報処理装置２であって、電力系統１の配電線１１のセンサ付開閉器ＳＷ１からの測定値に基づいて、所定区間における有効電力及び無効電力を計算する計算部２１２と、複数の需要家のそれぞれが需要家Ａ～Ｃのいずれに該当するかを示す分類情報と、複数の需要家の電力量計ＳＭの測定値と、に基づいて、需要家Ａ～Ｃのうち少なくともいずれかの有効電力を計算する計算部２１３と、需要家Ａ～Ｃのそれぞれの負荷における有効電力の比率Ｒ１と、需要家Ｂ及びＣのそれぞれの発電設備における有効電力の比率Ｒ２と、発電設備の力率と、所定区間における有効電力及び無効電力と、需要家Ａ～Ｃのうち少なくともいずれかの有効電力と、に基づいて、所定区間に含まれる負荷の有効電力及び無効電力と、所定区間に含まれる発電設備の有効電力と、を含む状態ｘを推定する推定部２１５と、を備える。
実施形態の情報処理装置２は、電力系統１の所定区間において、負荷を有する需要家Ａと、余剰買取契約対象の発電設備を有する需要家Ｂと、全量買取契約対象の発電設備を有する需要家Ｃと、のいずれかに分類される需要家を複数含む電力系統１の状態ｘを推定する情報処理装置２であって、電力系統１の配電線１１のセンサ付開閉器ＳＷ１からの測定値に基づいて、所定区間における有効電力及び無効電力を計算する計算部２１２と、複数の需要家のそれぞれが需要家Ａ～Ｃのいずれに該当するかを示す分類情報と、複数の需要家の電力量計ＳＭの測定値と、に基づいて、需要家Ａ～Ｃのうち少なくともいずれかの有効電力を計算する計算部２１３と、需要家Ａ～Ｃのそれぞれの負荷における有効電力の比率Ｒ１と、需要家Ｂ及びＣのそれぞれの発電設備における有効電力の比率Ｒ２と、発電設備の力率と、所定区間における有効電力及び無効電力と、需要家Ａ～Ｃのうち少なくともいずれかの有効電力と、に基づいて、所定区間に含まれる負荷の有効電力及び無効電力と、所定区間に含まれる発電設備の有効電力と、を含む状態ｘを推定する推定部２１５と、を備える。
【0202】
このような構成によれば、日射量計等の設備を要せずに、センサ付開閉器ＳＷ１の測定値と、電力量計ＳＭの測定値とを用いて、状態ｘを推定することができる。従って、このような構成によれば、低コストで、需要家の電力買取契約の形態に依らずに、需要家の負荷の有効電力等と、太陽光発電設備の有効電力とを推定することが可能となる。
【0203】
また、情報処理装置２において、計算部２１３は、需要家Ａ～Ｃのうち少なくともいずれか２つの有効電力を計算し、推定部２１５は、比率Ｒ１と、比率Ｒ２と、力率と、所定区間における有効電力及び無効電力と、需要家Ａ～Ｃのうち少なくともいずれか２つの有効電力に基づいて、状態ｘを推定してもよい。このような構成によれば、状態ｘを推定するのに用いる情報の量が増えるため、精度良く状態ｘを推定することができる。
【0204】
また、情報処理装置２において、計算部２１３は、需要家Ａ～Ｃの有効電力を計算し、推定部２１５は、比率Ｒ１と、比率Ｒ２と、力率と、所定区間における有効電力及び無効電力と、需要家Ａ～Ｃの有効電力に基づいて、状態ｘを推定してもよい。このような構成によれば、状態ｘを推定するのに用いる情報の量が更に増えるため、更に精度良く状態ｘを推定することができる。
【0205】
また、情報処理装置２において、推定部２１５は、センサ付開閉器ＳＷ１からの測定値に基づいて計算された所定区間における有効電力、及び状態ｘに基づいて計算された所定区間における有効電力の第１の差と、センサ付開閉器ＳＷ１からの測定値に基づいて計算された所定区間における無効電力、及び力率と、状態ｘとに基づいて計算された所定区間における無効電力の第２の差と、電力量計ＳＭの測定値に基づいて計算された需要家Ａ～Ｃのうち少なくともいずれかの需要家の有効電力、及び比率Ｒ１と、比率Ｒ２と、状態ｘとに基づいて計算された少なくともいずれかの需要家の有効電力の差とに基づく目的関数を設定し、状態ｘを決定変数とした際に目的関数を最小化することにより、状態ｘを推定する。このような構成によれば、更に精度良く状態ｘを推定することができる。
【0206】
また、情報処理装置２において、比率Ｒ１は、負荷の有効電力のうち、需要家Ａの負荷の有効電力の割合を示す値αＬと、需要家Ｂの負荷の有効電力の割合を示す値αＬ＿ＰＶであり、比率Ｒ２は、発電設備の有効電力のうち、需要家Ｂの発電設備の有効電力の割合を示す値βＬ＿ＰＶと、需要家Ｃの発電設備の有効電力の割合を示す値βＰＶであり、比率Ｒ１と、比率Ｒ２とを推定する推定部２１４を更に備え、推定部２１４は、所定区間に含まれる負荷の有効電力、及び需要家Ａの負荷の有効電力の測定値を値αＬで除した値の差と、所定区間に含まれる発電設備の有効電力、及び需要家Ｃの発電設備の有効電力の測定値を値βＰＶで除した値の差とに基づく目的関数を設定し、比率Ｒ１と、比率Ｒ２と、所定区間に含まれる負荷の有効電力と、所定区間に含まれる発電設備の有効電力とを決定変数とした際に目的関数を最小化することにより、比率Ｒ１と、比率Ｒ２とを推定する。このような構成によれば、比率Ｒ１と比率Ｒ２との精度が向上するため、更に精度良く状態ｘを推定することができる。
【0207】
また、情報処理装置２において、推定部２１４は、時間依存性を有する、比率Ｒ１と、比率Ｒ２とを、第２決定変数とする。このような構成によれば、時刻に応じた比率Ｒ１と比率Ｒ２とを用いることができるため、更に精度良く状態ｘを推定することができる。
【0208】
また、本実施形態の方法によれば、日射量計等の設備を要せずに、センサ付開閉器ＳＷ１の測定値と、電力量計ＳＭの測定値とを用いて、状態ｘを推定することができる。従って、このような構成によれば、低コストで、需要家の電力買取契約の形態に依らずに、需要家の負荷の有効電力等と、発電設備の有効電力とを推定することが可能となる。
【0209】
また、本実施形態の情報処理装置２に用いられるプログラムによれば、日射量計等の設備を要せずに、センサ付開閉器ＳＷ１の測定値と、電力量計ＳＭの測定値とを用いて、状態ｘを推定することができる。従って、このような構成によれば、低コストで、需要家の電力買取契約の形態に依らずに、需要家の負荷の有効電力等と、発電設備の有効電力とを推定することが可能となる。
【符号の説明】
【0210】
１：電力系統
１０：配電変電所
１１：配電線
２：情報処理装置
２００：ＣＰＵ
２０１：メモリ
２０２：通信装置
２０３：記憶装置
２０４：入力装置
２０５：出力装置
２０６：記録媒体読取装置
２１０：取得部
２１１：取得部
２１２：計算部
２１３：計算部
２１４：推定部
２１５：推定部
２１６：出力部
１０：配電変電所
１１：配電線
２：情報処理装置
２００：ＣＰＵ
２０１：メモリ
２０２：通信装置
２０３：記憶装置
２０４：入力装置
２０５：出力装置
２０６：記録媒体読取装置
２１０：取得部
２１１：取得部
２１２：計算部
２１３：計算部
２１４：推定部
２１５：推定部
２１６：出力部
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
