(57)【要約】  本願の電力制御装置（１）は、配電系統（８００）をリーマン多様体とみなし、配電系統（８００）内での電力損失を評価する関数（Ｌ）と、蓄電池の目標充電量からのずれを示す関数（Ｓ）の和を示す評価関数（Ｅ）を最小化する電力配分計画（ｕ０：ｔｆ、ｚ０：ｔｆ）を生成する配分計画部（１５）、および電力配分計画（ｕ０：ｔｆ、ｚ０：ｔｆ）と配分状態（ハットｚｔ）との間にずれが生じた際、電力配分計画と配分状態との距離が最小となるよう、接ベクトルと点を相互に滑らかに変換する写像を用いて電力配分計画の調整値を生成する配分調整部（１２）、を備えるように構成した。
2021014483.tif 000002
【特許請求の範囲】
【請求項１】
蓄電池を含む配電系統をリーマン多様体とみなし、前記配電系統において予測される電力供給量と電力需要から、前記蓄電池の残量を並べたベクトルと前記配電系統での制御可能な量を並べたベクトルとを変数とする前記配電系統での電力損失を評価する関数と、前記蓄電池の目標充電量からのずれを示す関数の和を示す評価関数を最小化する電力配分計画を生成する配分計画部、
前記配電系統における実際の配分状態を取得する配分状態取得部、
前記電力配分計画に基づいて前記配電系統での電力配分を制御する電力配分制御部、および
前記電力配分計画と前記配分状態との間にずれが生じた際、前記電力配分計画と前記配分状態との距離が最小となるよう、接ベクトルと点を相互に滑らかに変換する写像を用いて前記電力配分計画の調整値を生成する配分調整部、
を備えたことを特徴とする電力制御装置。
【請求項２】
前記配分調整部には、
前記配分状態の対数写像の値である第１の接ベクトルを計算する対数変換部、
前記対数変換部によって計算された第１の接ベクトルに係数を乗じ、接空間におけるフィードバック入力である第２の接ベクトルを計算する制御入力計算部、および
前記制御入力計算部によって計算された第２の接ベクトルの指数写像の値である前記調整値を計算する指数変換部、
が設けられていることを特徴とする請求項１に記載の電力制御装置。
【請求項３】
蓄電池を含む配電系統をリーマン多様体とみなし、前記配電系統において、電力供給量と電力需要を予測する予測ステップ、
前記予測された電力供給量と電力需要から、前記蓄電池の残量を並べたベクトルと前記配電系統での制御可能な量を並べたベクトルとを変数とする前記配電系統での電力損失を評価する関数と、前記蓄電池の目標充電量からのずれを示す関数の和を示す評価関数を最小化する電力配分計画を生成する電力配分計画生成ステップ、
前記電力配分計画に基づいて前記配電系統での電力配分を制御する電力配分制御ステップ、
前記配電系統における実際の配分状態を取得する配分状態取得ステップ、および
前記電力配分計画と前記配分状態との間にずれが生じた際、前記電力配分計画と前記配分状態との距離が最小となるよう、接ベクトルと点を相互に滑らかに変換する写像を用いて前記電力配分計画の調整値を生成する調整値生成ステップ、
を含むことを特徴とする電力配分計画設定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本願は、電力制御装置および電力配分計画設定方法に関するものである。
【背景技術】
【０００２】
再生可能エネルギー等の普及に伴い、配電側にも電源を持つようになっている。そこで、再生可能エネルギーの発電量も含めて電力需要を予測し、コストが最小になるように、配電系統の能力に応じて需要予測値と供給計画値の差分を分配する電力系統の制御システムが開示されている（例えば、特許文献１参照。）。しかし、再生可能エネルギーの発電可能量、または電力需要の予測値は、必ずしも実際の発電可能量、または電力需要と一致せず、通常は誤差を持つ。したがって、事前に計算した配電量に基づき配電系統を制御したとしても、実際の配電系統に関して最適な制御とはならない。【０００３】
この問題に対し、あらかじめ発電可能量、または電力需要に関する変動幅を確率分布として表示し、電力供給コストをより現実的に見積もる手法が開示されている（例えば、特許文献２参照。）。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】国際公開第２０１４／１１５５５６号（段落００１９〜００６１、図１〜図７）
【特許文献２】特開２０１２−３４４４４号公報（段落００１５〜００３６、図１〜図６）
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、見積もりが現実的であったとしても、確率分布の状態では、実際の配電量の制御に活かすことは困難であり、実態に応じた見積もりに基づく制御を行うことは困難であった。【０００６】
本願は、上記のような課題を解決するための技術を開示するものであり、実態に応じた見積もりに基づく制御を可能とする電力制御装置および電力配分計画設定方法を得ることを目的とする。
【課題を解決するための手段】
【０００７】
本願に開示される電力制御装置は、蓄電池を含む配電系統をリーマン多様体とみなし、前記配電系統において予測される電力供給量と電力需要から、前記蓄電池の残量を並べたベクトルと前記配電系統での制御可能な量を並べたベクトルとを変数とする前記配電系統での電力損失を評価する関数と、前記蓄電池の目標充電量からのずれを示す関数の和を示す評価関数を最小化する電力配分計画を生成する配分計画部、前記配電系統における実際の配分状態を取得する配分状態取得部、前記電力配分計画に基づいて前記配電系統での電力配分を制御する電力配分制御部、および前記電力配分計画と前記配分状態との間にずれが生じた際、前記電力配分計画と前記配分状態との距離が最小となるよう、接ベクトルと点を相互に滑らかに変換する写像を用いて前記電力配分計画の調整値を生成する配分調整部、を備えたことを特徴とする。【０００８】
本願に開示される電力配分計画設定方法は、蓄電池を含む配電系統をリーマン多様体とみなし、前記配電系統において、電力供給量と電力需要を予測する予測ステップ、前記予測された電力供給量と電力需要から、前記蓄電池の残量を並べたベクトルと前記配電系統での制御可能な量を並べたベクトルとを変数とする前記配電系統での電力損失を評価する関数と、前記蓄電池の目標充電量からのずれを示す関数の和を示す評価関数を最小化する電力配分計画を生成する電力配分計画生成ステップ、前記電力配分計画に基づいて前記配電系統での電力配分を制御する電力配分制御ステップ、前記配電系統における実際の配分状態を取得する配分状態取得ステップ、および前記電力配分計画と前記配分状態との間にずれが生じた際、前記電力配分計画と前記配分状態との距離が最小となるよう、接ベクトルと点を相互に滑らかに変換する写像を用いて前記電力配分計画の調整値を生成する調整値生成ステップ、を含むことを特徴とする。
【発明の効果】
【０００９】
本願に開示される電力制御装置あるいは電力配分計画設定方法によれば、予測と実態が一致しなかった場合でも、フィードバックが可能で、実態に応じた見積もりに基づく制御を可能とする。
【図面の簡単な説明】
【００１０】
【図１】実施の形態１にかかる電力制御装置の構成を説明するためのブロック図である。
【図２】実施の形態１にかかる電力制御装置の対数変換部の動作、あるいは電力配分計画設定方法における第１の接ベクトルを算出する手順を示す模式図である。
【図３】実施の形態１にかかる電力制御装置の配分計画部の動作、あるいは電力配分計画設定方法における電力配分計画生成ステップを示すフローチャートである。
【図４】実施の形態１にかかる電力制御装置の対数変換部の動作、あるいは電力配分計画設定方法における対数変換ステップを示すフローチャートである。
【図５】実施の形態１にかかる電力制御装置の動作、あるいは電力配分計画設定方法を示すフローチャートである。
【図６】実施の形態１にかかる電力制御装置、あるいは電力配分計画設定方法の効果を説明するためのグラフ形式の図である。
【図７】実施の形態１にかかる電力制御装置の演算処理を実行する部分の構成例を示すブロック図である。
【発明を実施するための形態】
【００１１】
実施の形態１．
図１〜図７は、実施の形態１にかかる電力制御装置の構成と動作、および電力配分計画設定方法について説明するためのものであり、図１は電力制御装置の構成を説明するためのブロック図、図２は電力制御装置の対数変換部の動作、あるいは電力配分計画設定方法における第１の接ベクトルを算出する手順を示す模式図、図３は配分計画部の動作、あるいは電力配分計画設定方法における電力配分計画生成ステップを示すフローチャート、図４は配分調整部を構成する対数変換部の動作、あるいは電力配分計画設定方法における対数変換ステップを示すフローチャートである。
【００１２】
また、図５は電力制御装置の全体動作、あるいは電力配分計画設定方法の全体動作を示すフローチャート、図６は電力制御装置、あるいは電力配分計画設定方法の効果を説明するための、状況変化が発生したときの、計画上の充電量の経時変化と、本実施の形態によるフィードバックを行った際の充電量の経時変化を示すグラフ形式の図である。さらに、図７は電力制御装置の演算処理を実行する部分の構成例を示すブロック図である。
【００１３】
実施の形態１にかかる電力制御装置１は、図１に示すように、電力供給量を予測する供給予測部１３、電力需要を予測する需要予測部１４、供給予測部１３と需要予測部１４の予測に基づき、電力配分計画を生成する配分計画部１５を備えている。そして、配電系統８００の配分状態を取得する配分状態取得部１６と、取得した配分状態と計算された電力配分計画とのずれを評価して、配分を調整する配分調整部１２と、調整した配分により、配電系統８００の電力配分を制御する電力配分制御部１１とを備える。なお、配電系統８００には、図示しない一つ以上の電力源と、負荷、および充放電を行う蓄電池が含まれているものとする。
【００１４】
供給予測部１３は、配電系統８００に含まれる太陽光発電、風力発電、燃料電池等の電力源による電力供給量を予測し、予測供給量Ｑｐを配分計画部１５に出力する。また、需要予測部１４は、配電系統８００に含まれる負荷等の電力需要を予測し、予測需要量Ｄｐを配分計画部１５に出力する。また、配分状態取得部１６は、配電系統８００における実際の電力供給と電力需要の状態ハットｚｔを取得し、配分調整部１２に出力する。
【００１５】
＜配分計画部＞
配分計画部１５は、供給予測部１３から受信した予測供給量Ｑｐと需要予測部１４から受信した予測需要量Ｄｐに基づき、最適な電力配分計画ｕｔ、ｚｔとフィードバック係数Ｈｔ，νν、Ｈｔ，νｚを算出する。電力配分計画ｕｔ、ｚｔの算出には、蓄電池の残量を並べたベクトルｘｔ、配電系統８００内での再生可能エネルギーでの発電量、融通される電力量等を並べたベクトルｙｔ、および蓄電池の充放電電力を含む電力制御装置１で制御可能な量を並べたベクトルである入力ｕｔを用いる。そして、ベクトルｘｔ、ｙｔをまとめたベクトルｚｔを用い、目標充電量からのずれを表す関数と、電力ロス、コストロス等の経済的損失を表す関数の和として定義される、後述する評価関数Ｅを最小化する電力配分計画ｕｔ、ｚｔを求める。そして、フィードバック係数Ｈｔ，νν、Ｈｔ，νｚとともに、配分調整部１２に出力する。
【００１６】
＜配分調整部＞
配分調整部１２は、第１の接ベクトルｄｚを算出する対数変換部１２１と、第２の接ベクトルｄｕを算出する制御入力計算部１２２と、フィードバック入力ｕを算出して、電力配分制御部１１に出力する指数変換部１２３とを備えている。そして、配電系統８００がなす多様体において、電力配分計画から配電系統８００の状態に向けた接ベクトルｄｚを用いて接空間におけるフィードバック入力ｄｕを算出し、それを実際のフィードバック入力ｕに変換して、電力配分制御部１１に出力する。
【００１７】
＜対数変換部＞
対数変換部１２１は、配分計画部１５から出力された現在の電力配分計画ｕｔ、ｚｔと、配分状態取得部１６から出力された現在の状態ハットｚｔから、対数写像の値である第１の接ベクトルｄｚを算出する。より具体的には、配電系統８００がなす多様体において、現在の電力配分計画ｕｔ、ｚｔから、現在の配電系統８００の現在の状態であるハットｚｔに向けた、第１の接ベクトルｄｚを算出し、制御入力計算部１２２に出力する。
【００１８】
＜制御入力計算部＞
制御入力計算部１２２は、対数変換部１２１から出力された第１の接ベクトルｄｚと、配分計画部１５から出力されたフィードバック係数Ｈｔ，νν、Ｈｔ，νｚから、接空間におけるフィードバック入力である第２の接ベクトルｄｕを算出する。そして、算出した第２の接ベクトルｄｕを指数変換部１２３に出力する。
【００１９】
＜指数変換部＞
指数変換部１２３は、第２の接ベクトルｄｕを指数変換して、実際のフィードバック入力ｕを算出し、電力配分制御部１１に出力する。以下、詳細に説明する。
【００２０】
本実施の形態１にかかる電力制御装置１あるいは電力配分計画設定方法における電力配分問題は、評価関数Ｅ（式１）を最小化する入力ｕｔ：ｔｆを求める問題として定式化される。
【００２１】
【数１】
2021014483.tif 000003
ただし、ｕ０：ｔｆとは、ｕ０，ｕ１，・・・，ｕｔｆ−１をまとめたものである。またｘｔ，ｙｔ，ｕｔは、方程式ｘｔ＋１＝ｆ（ｘｔ，ｙｔ，ｕｔ）、等式制約ｈ（ｘｔ，ｙｔ）＝０、およびｌ（ｘｔ，ｙｔ，ｕｔ）＝０を満たす。
【００２２】
ここで、ｘｔは１つまたは複数の蓄電池の残量を並べたベクトルである。また、ｙｔは太陽光、風力等の発電量、配電系統内で融通される電力等を並べたベクトルであってよい。入力ｕｔは蓄電池の充放電電力、発電機、あるいは上位電力系統からの受電量、買電量等、電力制御装置によって制御可能な量を並べたベクトルであってよい。関数Ｌは蓄電池の充放電損失、変換器の変換損失、電力あたりの買電価格等を表す関数であってよい。関数Ｓは蓄電池の目標充電量からのずれを表す関数であってよい。つまり、配電系統８００内での電力損失あるいは買電価格等を表す関数Ｌの積分と、蓄電池の目標充電量からのずれを示す関数Ｓとの和を示す評価関数Ｅを最小化する入力ｕ０：ｔｆを求めることになる。
【００２３】
このとき、任意のｘに対し、ｈ（ｘ，ｙ）＝０を満たすｙは少なくとも局所的に一意に定まる。ここでｘｔ、ｙｔをまとめたベクトルをｚｔとして、ｚｔとｕｔからｚｔ＋１を与える写像をＦとおく。評価関数Ｅの、時刻ｔからｔｆまでの部分和をＪｔ（ｚｔ，ｕｔ：ｔｆ）とおくと、部分和は式（２）を満たす。
Jt(zt,ut:tf)＝L(zt,ut)+Jt+1(F(zt,ut),ut+1:tf)  ・・・（２）
ただし、Ｊｔｆ（ｚｔｆ）＝Ｓ（ｚｔｆ）である。そして、動的計画法を適用することで、式（３）が成り立つ。
【００２４】
【数２】
2021014483.tif 000004
ここで、Ｑｔ（ｚ，ｕ）を式（３ｂ）のように定義する。ラグランジュ（Lagrange）の未定定数法より、ｌ（ｚ，ｕ）＝０を満たす範囲でＱｔ（ｚ，ｕ）を最小にするｕは、Ｈｔ（ｚ，ｕ，λ）：＝Ｑｔ（ｚ，ｕ）＋λΤｌ（ｚ，ｕ）を最小にするｕによって求められる。入力ｕと変数λをまとめたベクトルをνとおき、Ｈｔ（ｚ，ν）の２次までのテイラー（Taylor）展開をとると、式（４）であるから、これを最小にするｄνは式（５）で表される。
【００２５】
【数３】
2021014483.tif 000005
ただし、ｉ、ｊはνに関する添字、ｉ´、ｊ´は、ｚに関する添字である。
関数Ｖｔ（ｚｔ）は、Ｈｔ（ｚ，ν）のνに関する最小値であるから、式（５）を式（４）に代入して、式（６）を得る。
【００２６】
【数４】
2021014483.tif 000006
ここから、以下の関係式（式（７））を得る。
∇i´V＝∇j´Ht-(∇iHt)(∇ijHt)(∇ij´Ht)
∇i´j´V＝∇i´j´Ht-(∇i´jHt)(∇ijHt)(∇ij´Ht)  ・・・（７）
一方、式（８）の定義から、式群（９）および式群（１０）を得る。
【００２７】
【数５】
2021014483.tif 000007
ただし、ｉ、ｊはｕに関する添字であり、ｉ´´、ｊ´´は、λに関する添字である。ここでＦｋ´のｙに関する微分は、関数ｈ（ｘ，ｙ）が、ｙに関するヤコビ行列が正則であることから計算できる。以上の式を用いると、時刻ｔｆ−１から再帰的にｔ＝０までのＨｔ、Ｖｔの微分を求めることができる。
【００２８】
さて、本願の実施の形態にかかる電力制御装置１あるいは電力配分計画設定方法における特徴は、ある入力ｕ０：ｔｆに対し、更新後の入力ハットｕ０：ｔｆを式群（１１）によって求める点にある。
【００２９】
【数６】
2021014483.tif 000008
ただし、Ｅｘｐｕｔとは、等式制約ｈ（ｚ）＝０、およびｌ（ｚ，ｕ）＝０のなすリーマン多様体の、点ｕｔにおける指数写像であり、Ｌｏｇｚｔは等式制約ｈ（ｚ）＝０のなすリーマン多様体の、点ｚｔにおける対数写像である。
【００３０】
まず、等式制約ｈ（ｚ）＝０のなす多様体において、ｚの測地線方程式は、式（１２）で与えられる。
【００３１】
【数７】
2021014483.tif 000009
ただし、Δｋｌｈとは、Δｋｌｈ＝（δｉｈｋ）δｉｊ（δｊｈｌ）の逆行列であり、δｉｊはユークリッド計量δｉｊの逆行列である。この方程式の、初期条件ｚ（０）＝ｚｔ、ドットｚ（０）＝ｄｚにおけるｚ（１）の値が、指数写像Ｅｘｐｚｔ［ｄｚ］の値である。この値は公知の微分方程式の数値解法によって容易に計算することができる。
【００３２】
ここで、ｈ（ｚ）＝０のなすリーマン多様体上の非負かつ非退化な関数Ｄ（ｚ，ｚ´）で、少なくともｚ´の近傍でｚに関するヘッシアン（Hessian）が正定値となるようなものを考える。本実施形態において、Ｄ（ｚ，ｚ´）はユークリッド距離の２乗としてもよい。このＤ（ｚ，ｚ´）を用いて対数写像Ｌｏｇｚｔハットｚｔを効率的に求めることができる。
【００３３】
今、Ｄ（ｚ，ハットｚｔ）をｚの関数と考えると、Ｄ（ｚ，ハットｚｔ）は、ハットｚｔにおいて最小値０をとる凸関数である。従って、点ｚｔにおけるテイラー（Taylor）展開(式（１３）)に基づくＤ（ｚ，ハットｚｔ）の最小化方向ｄｚｉ＝−（∇ｉｊＤ）∇ｊＤは、Ｅｘｐｚｔ［ｄｚ］がハットｚｔとなるようなｄｚを近似する。
【００３４】
【数８】
2021014483.tif 000010
ここで、ｒ＝Ｅｘｐｚｔ［ｄｚ］において、再びＤ（ｚ，ハットｚｔ）をテイラー展開し、その最小化方向ｄチェックｚ（ｒ）を等式制約ｈ（ｚ）＝０のなす多様体上の平行移動方程式（式（１４））によって点ｚｔに引き戻すことで、新たな最小化方向ｄｚ＋ｄチェックｚ（ｚｔ）を構成することができる。
【００３５】
【数９】
2021014483.tif 000011
この操作を繰り返すことで、ｒは急速にハットｚｔに近づき、ほぼ一致した時点での最小化方向が対数写像Ｌｏｇｚｔハットｚｔの値を与える。この手順を図式化した図２を用いて説明する。
【００３６】
点ｚｔから距離関数Ｄ（ｚ，ハットｚｔ）を最小化する方向へのベクトルがｄｚであり、ｄｚ方向への測地線の先がｒである。ｒからＤ（ｚ，ハットｚｔ）を最小化する方向へのベクトルがｄチェックｚ（ｒ）であり、これを測地線にそって点ｚｔに引き戻したベクトルがｄチェックｚ（ｚｔ）である。新たな最小化方向はｄｚ＋ｄチェックｚ（ｚｔ）であり、その方向への測地線はｒよりもハットｚｔに近づくことがわかる。
【００３７】
等式制約ｈ（ｚ）＝０、およびｌ（ｚ，ｕ）＝０のなすリーマン多様体の、点ｕｔにおける指数写像Ｅｘｐｕｔは、ｚとｕを新たにｗとおくことで、式（１２）と同様の測地線方程式が得られる。得られた測地線方程式のうち、ドットｚに関する部分を、対数写像Ｌｏｇｚｔを求めたときの測地線とすることで得られる。
【００３８】
つぎに、図３のフローチャートを用いて、配分計画部１５における電力配分計画ｕ０：ｔｆ、ｚ０：ｔｆの生成手順について説明する。はじめに、初期状態ｘ０、ｙ０、および初期入力ｕ０：ｔｆを定める（ステップＳ２０１）。つづいて、ｚｔ＋１＝Ｆ（ｚｔ，ｕｔ）に従い、ｚｔとｕｔからｚｔ＋１を計算する（ステップＳ２０２）。なお、ステップＳ２０２は、ｔ＝０からｔｆ−１まで繰り返される。
【００３９】
次に、式（７）、（９）、（１０）を用い、Ｖｔ＋１の微分からＨｔとＶｔの微分を求める（ステップＳ２０３）。ステップＳ２０３は、ｔ＝ｔｆ−１からｔ＝０まで繰り返される。そして、式（１５）に示す勾配のノルムを算出し、これが所定のしきい値αより小さいかどうかを判定する（ステップＳ２０４）。
【００４０】
【数１０】
2021014483.tif 000012
勾配のノルムが、しきい値αよりも小さいならば（「Ｙｅｓ」）、現在のｕ０：ｔｆとｚ０：ｔｆ、およびＨｔとＶｔの微分（Ｈｔ，νν、Ｈｔ，νｚ）を電力配分計画、およびフィードバック係数として出力する（ステップＳ２０６）。一方、しきい値α以上であれば（「Ｎｏ」）、式群（１１）に従い、指数写像と対数写像によって、新たな入力ハットｕ０：ｔｆ、および状態ハットｚ０：ｔｆを計算し（ステップＳ２０５）、ステップＳ２０３の繰り返しに戻る。
【００４１】
つぎに、図４のフローチャートを用いて、対数変換部１２１での第１の接ベクトルｄｚの計算手順（対数変換ステップ）を表す。配分計画部１５が出力した時刻ｔでの電力配分計画ｕｔ、ｚｔと配分状態取得部１６から取得した現在の状態ハットｚｔから、点（ｕｔ，ｚｔ）におけるＤ（ｚ，ハットｚｔ）の最小化方向ｄｚを計算する（ステップＳ３０１）。つづいて、点（ｕｔ，ｚｔ）から、ｄｚ方向への測地線を計算し、ｒ＝Ｅｘｐｚｔ［ｄｚ］を求める（ステップＳ３０２）。
【００４２】
そして、Ｄ（ｒ，ハットｚｔ）のノルムが、所定のしきい値βより小さいかどうかを判定する（ステップＳ３０３）。Ｄ（ｒ，ハットｚｔ）のノルムが、しきい値βよりも小さいならば（「Ｙｅｓ」）、ｄｚを出力する（ステップＳ３０７）。一方、しきい値β以上であれば（「Ｎｏ」）、点ｒにおける最小化方向ｄチェックｚを計算し（ステップＳ３０４）、平行移動方程式（式（１４））にしたがい、ｄチェックｚを点ｚｔまで平行移動する（ステップＳ３０５）。続いて、ステップＳ３０５によって得られた方向ｄチェックｚをｄｚに足し合わせ、新たなｄｚとする(ステップＳ３０６)。
【００４３】
制御入力計算部１２２は、配分計画部１５が計算したフィードバック係数Ｈｔ，νν、Ｈｔ，νｚから、式（１６）により、接空間におけるフィードバック入力であるｄｕを計算する。
du＝-[Ht,νν-1Ht,νz]udz  ・・・（１６）
ここで、ｄｚは、第１の接ベクトルに相当し、ｄｕは、第２の接ベクトルに相当する。  
【００４４】
指数変換部１２３は、ｕ＝Ｅｘｐｕｔｄｕにより、接空間におけるフィードバック入力ｄｕから、実際のフィードバック入力ｕを計算し、電力配分制御部１１は、計算されたフィードバック入力ｕに基づき、配電系統８００における電力配分を制御する。
【００４５】
上述した動作の全体、つまり電力配分計画設定方法について図５のフローチャートに基づき説明する。はじめに、供給予測部１３による予測供給量Ｑｐの生成、需要予測部１４による予測需要量Ｄｐの生成が行われ、配分計画部１５に出力される（ステップＳ１００）。
【００４６】
配分計画部１５は、予測供給量Ｑｐと予測需要量Ｄｐをもとに、上述したベクトルｘｔ、ｙｔ、ｕｔを用いて、蓄電池の目標充電量からのずれを示す関数と、配電系統８００内での損失を示す関数の和として定義した評価関数Ｅを最小化する電力配分計画を求める。そして、写像Ｆによりｕ０：ｔｆ、ｚ０：ｔｆを、式（７）〜（１０）により、Ｈｔ、Ｖｔの微分を計算し、条件（||∇Ｈ||＜α）を満たすことを前提として、それぞれ電力配分計画、フィードバック係数として出力する（ステップＳ１１０（ステップＳ２０１〜Ｓ２０６））。
【００４７】
このとき、電力配分計画は、配分調整部１２の対数変換部１２１に出力され、フィードバック係数は、制御入力計算部１２２に出力される。なお、条件（||∇Ｈ||＜α）を満たさない場合は、式群（１１）を用いて、新たな入力と状態を算出し、条件を満たすまで再計算を行う。
【００４８】
電力配分制御部１１は、少なくとも電力配分計画が生成された直後は、生成された電力配分計画ｕ０：ｔｆ、ｚ０：ｔｆに基づいて、配電系統８００内の電力配分制御を実行する（ステップＳ１２０）。
【００４９】
一方、配分状態取得部１６は、配電系統８００から現実の配分状態ハットｚｔを取得し（ステップＳ１３０）、配分調整部１２の対数変換部１２１に出力する。対数変換部１２１は、配分状態ハットｚｔと電力配分計画ｕ０：ｔｆ、ｚ０：ｔｆを受信すると、計画と実態（電力配分計画ｕ０：ｔｆ、ｚ０：ｔｆと配分状態ハットｚｔ）が一致するか否かを判定（ステップＳ１４０）し、一致している場合（Ｙｅｓ）は、ステップＳ１２０へ戻る。一方、一致しなかった場合（Ｎｏ）は、電力配分計画の修正値であるフィードバック入力ｕを演算するため、現在の状態ハットｚｔから、対数写像の値である第１の接ベクトルｄｚを算出し（ステップＳ１５０（ステップＳ３０１〜Ｓ３０７））、制御入力計算部１２２に出力する。
【００５０】
制御入力計算部１２２は、第１の接ベクトルｄｚと、配分計画部１５からのフィードバック係数Ｈｔ，νν、Ｈｔ，νｚを用い、接空間におけるフィードバック入力である第２の接ベクトルｄｕを算出し（ステップＳ１６０）、指数変換部１２３に出力する。指数変換部１２３は、第２の接ベクトルｄｕを指数変換して、実際のフィードバック入力ｕを算出し（ステップＳ１７０）、電力配分制御部１１に出力し、ステップＳ１２０に戻る。
【００５１】
このような電力配分制御における効果について、図６を用いて説明する。図において、横軸は時刻であり、縦軸は蓄電池の充電量である。実線は配電計画時点における予測供給量と予測需要量に基づく最適な充放電曲線Ｐを表す。
【００５２】
実際の電力供給量が計画を上回る等の理由により、早期に蓄電池を充電したほうが目的関数を減少させられる状況が発生したとする。例えば、これまでの手法によれば、そのような状況が発生した時点（時刻t）で、最終時刻までの蓄電池の充放電計画を再計算しなければならない。しかし、本願の実施の形態１および以降の実施の形態２にかかる電力制御装置１または電力配分計画設定方法では、オンラインでのフィードバック制御によって、充放電計画を再計算することなく、目的関数（評価関数Ｅ）を減少させるような充放電曲線Ｍ（破線）を実現できる。またその際、蓄電池に関する制約（充電量は0から1の間）が考慮される。
【００５３】
ここで、配電系統８００が蓄電池を含む場合、最適な配電量は蓄電池の将来にわたる充電状態を考慮して決定されるため、誤差の補正もまた蓄電池の将来にわたる充電状態を考慮したものでなければならない。それに対して、上述した配分調整部１２を備えたことで、蓄電池を含む配電系統８００において、発電量または電力需要の予測が実際の発電量または需要と一致しない場合においても、蓄電池の将来にわたる充電状態を考慮に入れたフィードバック制御が実現できる。つまり、本実施の形態にかかる電力制御装置１あるいは電力配分計画設定方法では、蓄電池の将来にわたる充電状態を考慮して、最適な配電量を決定するので、実態への適応性と制御性を両立させることができる。
【００５４】
なお、本実施の形態に示す電力制御装置１は、オンラインの制御則が対数変換部１２１、制御入力計算部１２２、指数変換部１２３で構成されており、これらは微分方程式の初期値問題、およびその少数回の繰り返しと、単純な代数演算である。従って、最適化問題をオンラインで解く必要がなく、オンライン制御に適したものとなっている。
【００５５】
そして、本実施の形態で例示する配分計画部１５と、配分調整部１２（対数変換部１２１、制御入力計算部１２２、指数変換部１２３）は、例えば計算機内のプログラムとして実現されてもよい。その際、各ステップが必要とする情報がそれ以前のステップにて計算される限りにおいて、計算順序の入れ替えは許容される。また、指数変換部１２３、および対数変換部１２１における計算は、かならずしも指数写像と対数写像に限らなくてもよく、接ベクトルと多様体上の点を相互に滑らかに変換する写像であれば許容される。
【００５６】
その場合、配分計画部１５と配分調整部１２は、図７に示すように、プロセッサ１０１と記憶装置１０２を備えた一つのハードウェア１０によって構成することが考えられる。記憶装置１０２は、図示していないが、ランダムアクセスメモリ等の揮発性記憶装置と、フラッシュメモリ等の不揮発性の補助記憶装置とを具備する。また、フラッシュメモリの代わりにハードディスクの補助記憶装置を具備してもよい。プロセッサ１０１は、記憶装置１０２から入力されたプログラムを実行する。この場合、補助記憶装置から揮発性記憶装置を介してプロセッサ１０１にプログラムが入力される。また、プロセッサ１０１は、演算結果等のデータを記憶装置１０２の揮発性記憶装置に出力してもよいし、揮発性記憶装置を介して補助記憶装置にデータを保存してもよい。
【００５７】
実施の形態２．
本実施の形態２においては、電力配分問題が入力に関する不等式制約ｋ（ｘｔ，ｙｔ，ｕｔ）＞０を含む場合の電力配分における計算方法について説明する。なお、本実施の形態２にかかる電力制御装置あるいは電力配分計画設定方法では、不等式制約ｋ（ｘｔ，ｙｔ，ｕｔ）＞０を含むことにより、配分計画部と配分調整部で用いる計算式が異なる以外は、実施の形態１と同様である。そのため、実施の形態１で用いた図、あるいは同様部分の式については、援用し、同様部分の説明は繰り返さない。
【００５８】
蓄電池の残量に関する不等式制約等は、蓄電池の充放電電力、または充放電電流が変数ｕｔに含まれる場合、容易に入力に関する不等式制約に変換できる。本実施の形態２において、不等式制約は等式制約ｋ（ｘｔ，ｙｔ，ｕｔ）−ｓｔ＝０、および不等式制約ｓｔ＞０に置き換えられる。ただしｓｔ＞０とは、ｓｔの成分がすべて正であることを意味する。また、等式制約ｋ（ｘｔ，ｙｔ，ｕｔ）−ｓｔ＝０は、ｌ（ｚ，ｕ，ｓｔ）＝０に含めて取り扱う。
【００５９】
さて、ｓ＞０なるｓのなす多様体は、ｋ≧０とすると、計量ｇｉｊ（ｓ）＝ｋ＋（ｓｉ）−２を導入することで、リーマン多様体と考えることができる。組（ｚｔ，νｔ，ｓｔ）のなす多様体は、ｓｔのなす多様体と、ユークリッド空間の直積であり、これもまたリーマン多様体である。このリーマン多様体における、関数Ｈｔの微分は、νｔ＝（ｕｔ，λｔ，ｓｔ）を入力だとみなして、式（１７）と変わることを除き、実施の形態１と同様である。
∇ijHt＝∂ijHt-Γijk∇kHt  ・・・（１７）
【００６０】
ただし、δｉｊＨｔは、Ｈｔのνｉ、νｊに関する2階の偏微分であり、Γｉｊｋは、第２種クリストッフェル（Christoffel）記号である。ここでは、νｉがｓｊである場合、Γｉｉｉ＝−１／ｓｊであり、それ以外の添字の組み合わせに関してΓｉｊｋ＝０である。本実施の形態２では、実施の形態１においては∇ｉｊＨｔ＝δｉｊＨｔであった点に対して異なる。また、ｗ＝（ｚ，ｕ，ｓ）に関する測地線方程式については、ユークリッド計量δｉｊが計量ｇｉｊに置き換わることを除き、実施の形態１と同様である。
【００６１】
なお、本願は、様々な例示的な実施の形態および実施例が記載されているが、１つ、または複数の実施の形態に記載された様々な特徴、態様、および機能は特定の実施の形態の適用に限られるのではなく、単独で、または様々な組み合わせで実施の形態に適用可能である。従って、例示されていない無数の変形例が、本願明細書に開示される技術の範囲内において想定される。例えば、少なくとも１つの構成要素を変形する場合、追加する場合または省略する場合、さらには、少なくとも１つの構成要素を抽出し、他の実施の形態の構成要素と組み合わせる場合が含まれるものとする。
【００６２】
以上のように、各実施の形態にかかる電力制御装置１によれば、蓄電池を含む配電系統８００をリーマン多様体とみなし、配電系統８００内において予測される電力供給量（予測供給量Ｑｐ）と電力需要（予測需要量Ｄｐ）から、蓄電池の残量を並べたベクトルｘｔ（配電系統８００内で融通される電力等を並べたベクトルｙｔとまとめてベクトルｚｔ）と配電系統８００内での制御可能な量を並べたベクトルｕｔとを変数とする配電系統８００内での電力損失を評価する関数Ｌと、蓄電池の目標充電量からのずれを示す関数Ｓの和を示す評価関数Ｅ（式（１））を最小化する電力配分計画ｕ０：ｔｆ、ｚ０：ｔｆを生成する配分計画部１５、配電系統８００内における実際の配分状態ハットｚｔを取得する配分状態取得部１６、電力配分計画ｕ０：ｔｆ、ｚ０：ｔｆに基づいて配電系統８００内での電力配分を制御する電力配分制御部１１、および電力配分計画ｕ０：ｔｆ、ｚ０：ｔｆと配分状態ハットｚｔとの間にずれが生じた際、電力配分計画と配分状態との距離（例えば、Ｄ（ｚ，ハットｚｔ））が最小となるよう、接ベクトルと点を相互に滑らかに変換する写像を用いて電力配分計画の調整値（フィードバック入力ｕ）を生成する配分調整部１２、を備えるように構成したので、実態に応じた見積もりに基づく制御を可能とする。
【００６３】
とくに、配分調整部１２には、配分状態ハットｚｔの対数写像の値である第１の接ベクトルｄｚを計算する対数変換部１２１、対数変換部１２１によって計算された第１の接ベクトルｄｚに係数（フィードバック係数Ｈｔ，νν、Ｈｔ，νｚ）を乗じ、接空間におけるフィードバック入力である第２の接ベクトルｄｕを計算する制御入力計算部１２２、および制御入力計算部１２２によって計算された第２の接ベクトルｄｕの指数写像の値である調整値（フィードバック入力ｕ）を計算する指数変換部１２３、が設けられているように構成すれば、微分方程式の初期値問題、の小数回の繰り返しと単純な代数演算で調整値（フィードバック入力ｕ）を算出し、容易にオンライン制御できる。
【００６４】
また、各実施の形態にかかる電力配分計画設定方法によれば、蓄電池を含む配電系統８００をリーマン多様体とみなし、配電系統８００内において、電力供給量（予測供給量Ｑｐ）と電力需要（予測需要量Ｄｐ）を予測する予測ステップ（ステップＳ１００）、予測供給量Ｑｐと予測需要量Ｄｐから、蓄電池の残量を並べたベクトルｘｔ（配電系統８００内で融通される電力等を並べたベクトルｙｔとまとめてベクトルｚｔ）と配電系統８００内での制御可能な量を並べたベクトルｕｔとを変数とする配電系統８００内での電力損失を評価する関数Ｌと、蓄電池の目標充電量からのずれを示す関数Ｓの和を示す評価関数Ｅ（式（１））を最小化する電力配分計画ｕ０：ｔｆ、ｚ０：ｔｆを生成する電力配分計画生成ステップ（ステップＳ１１０（ステップＳ２０１〜Ｓ２０６））、電力配分計画ｕ０：ｔｆ、ｚ０：ｔｆに基づいて配電系統８００内での電力配分を制御する電力配分制御ステップ（ステップＳ１２０）、配電系統８００内における実際の配分状態ハットｚｔを取得する配分状態取得ステップ（ステップＳ１３０）、および電力配分計画ｕ０：ｔｆ、ｚ０：ｔｆと配分状態ハットｚｔとの間にずれが生じた際、電力配分計画と配分状態との距離例えば、Ｄ（ｚ，ハットｚｔ）が最小となるよう、接ベクトルと点を相互に滑らかに変換する写像を用いて電力配分計画の調整値（フィードバック入力ｕ）を生成する調整値生成ステップ（ステップＳ１４０〜Ｓ１７０）、を含むように構成したので、実態に応じた見積もりに基づく制御を可能とする。
【符号の説明】
【００６５】
１：電力制御装置、  １１：電力配分制御部、  １２：配分調整部、  １３：供給予測部、  １４：需要予測部、  １５：配分計画部、  １６：配分状態取得部、  １２１：対数変換部、  １２２：制御入力計算部、  １２３：指数変換部、  ８００：配電系統、  Ｄ（ｚ，ハットｚｔ）：距離、  Ｄｐ：予測需要量、  ｄｕ：第２の接ベクトル、  ｄｚ：第１の接ベクトル、  Ｈｔ，νν，Ｈｔ，νｚ：フィードバック係数、  Ｑｐ：予測供給量、  ｕ：フィードバック入力（調整値）、  ｘｔ，ｙｔ，ｚｔ：ベクトル、  ｕ０：ｔｆ,ｚ０：ｔｆ：電力配分計画、  ハットｚｔ：配分状態。
【図１】
2021014483.tif 000013
【図２】
2021014483.tif 000014
【図３】
2021014483.tif 000015
【図４】
2021014483.tif 000016
【図５】
2021014483.tif 000017
【図６】
2021014483.tif 000018
【図７】
2021014483.tif 000019
【手続補正書】
【提出日】令和3年7月5日(2021.7.5)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
蓄電池を含む配電系統をリーマン多様体とみなし、前記配電系統において予測される電力供給量と電力需要から、前記蓄電池の残量を並べたベクトルと前記配電系統での制御可能な量を並べたベクトルとを変数とする前記配電系統での電力損失を評価する関数と、前記蓄電池の目標充電量からのずれを示す関数の和を示す評価関数を最小化する電力配分計画を生成する配分計画部、
前記配電系統における実際の配分状態を取得する配分状態取得部、
前記電力配分計画に基づいて前記配電系統での電力配分を制御する電力配分制御部、および
前記電力配分計画が、前記配分状態に対してずれを有する際、フィードバック制御によって、前記ずれを有する電力配分計画と前記配分状態との距離が最小となるよう、接ベクトルと点を相互に滑らかに変換する写像を用いて前記ずれを有する電力配分計画の調整値を生成する配分調整部、
を備えたことを特徴とする電力制御装置。
【請求項２】
前記配分調整部には、
前記配分状態の対数写像の値である第１の接ベクトルを計算する対数変換部、
前記対数変換部によって計算された第１の接ベクトルに係数を乗じ、接空間におけるフィードバック入力である第２の接ベクトルを計算する制御入力計算部、および
前記制御入力計算部によって計算された第２の接ベクトルの指数写像の値である前記調整値を計算する指数変換部、
が設けられていることを特徴とする請求項１に記載の電力制御装置。
【請求項３】
蓄電池を含む配電系統をリーマン多様体とみなし、前記配電系統において、電力供給量と電力需要を予測する予測ステップ、
前記予測された電力供給量と電力需要から、前記蓄電池の残量を並べたベクトルと前記配電系統での制御可能な量を並べたベクトルとを変数とする前記配電系統での電力損失を評価する関数と、前記蓄電池の目標充電量からのずれを示す関数の和を示す評価関数を最小化する電力配分計画を生成する電力配分計画生成ステップ、
前記電力配分計画に基づいて前記配電系統での電力配分を制御する電力配分制御ステップ、
前記配電系統における実際の配分状態を取得する配分状態取得ステップ、および
前記電力配分計画が、前記配分状態に対してずれを有する際、フィードバック制御によって、前記ずれを有する電力配分計画と前記配分状態との距離が最小となるよう、接ベクトルと点を相互に滑らかに変換する写像を用いて前記ずれを有する電力配分計画の調整値を生成する調整値生成ステップ、
を含むことを特徴とする電力配分計画設定方法。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】０００７
【補正方法】変更
【補正の内容】
【０００７】
本願に開示される電力制御装置は、蓄電池を含む配電系統をリーマン多様体とみなし、前記配電系統において予測される電力供給量と電力需要から、前記蓄電池の残量を並べたベクトルと前記配電系統での制御可能な量を並べたベクトルとを変数とする前記配電系統での電力損失を評価する関数と、前記蓄電池の目標充電量からのずれを示す関数の和を示す評価関数を最小化する電力配分計画を生成する配分計画部、前記配電系統における実際の配分状態を取得する配分状態取得部、前記電力配分計画に基づいて前記配電系統での電力配分を制御する電力配分制御部、および前記電力配分計画が、前記配分状態に対してずれを有する際、フィードバック制御によって、前記ずれを有する電力配分計画と前記配分状態との距離が最小となるよう、接ベクトルと点を相互に滑らかに変換する写像を用いて前記ずれを有する電力配分計画の調整値を生成する配分調整部、を備えたことを特徴とする。
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】０００８
【補正方法】変更
【補正の内容】
【０００８】
本願に開示される電力配分計画設定方法は、蓄電池を含む配電系統をリーマン多様体とみなし、前記配電系統において、電力供給量と電力需要を予測する予測ステップ、前記予測された電力供給量と電力需要から、前記蓄電池の残量を並べたベクトルと前記配電系統での制御可能な量を並べたベクトルとを変数とする前記配電系統での電力損失を評価する関数と、前記蓄電池の目標充電量からのずれを示す関数の和を示す評価関数を最小化する電力配分計画を生成する電力配分計画生成ステップ、前記電力配分計画に基づいて前記配電系統での電力配分を制御する電力配分制御ステップ、前記配電系統における実際の配分状態を取得する配分状態取得ステップ、および前記電力配分計画が、前記配分状態に対してずれを有する際、フィードバック制御によって、前記ずれを有する電力配分計画と前記配分状態との距離が最小となるよう、接ベクトルと点を相互に滑らかに変換する写像を用いて前記ずれを有する電力配分計画の調整値を生成する調整値生成ステップ、を含むことを特徴とする。
【国際調査報告】
2021014483.tif 000020
2021014483.tif 000021
2021014483.tif 000022
2021014483.tif 000023