(57)【特許請求の範囲】
【請求項１】
状態を測定可能な第１ノードと、状態を測定できない第２ノードと、を含む電力系統におけるノードの状態を推定する電力系統の状態推定方法であって、
前記第１ノードの状態を示す測定値を取得し、
第１時刻より前の第２時刻における前記第２ノードの状態を示す計算値を取得し、
前記第１時刻における前記測定値と、前記第２時刻における前記測定値と、の差分を、前記第２ノードに属する各ノードに按分し、前記各ノードに按分された値と、該当する前記各ノードについての前記計算値と、の和を、前記第１時刻における前記第２ノードの状態を示す値を算出するための初期値として設定し、
前記初期値と、前記第１時刻における前記測定値と、を用いて最適化計算を実行し、前記第１時刻における前記第２ノードの状態を示す値を算出する
ことを特徴とする電力系統の状態推定方法。
【請求項２】
前記差分は、前記各ノードに接続された負荷の容量に応じて按分される
ことを特徴とする請求項１に記載の電力系統の状態推定方法。
【請求項３】
前記第２時刻は、前記第１時刻における前記測定値との偏差が最小となる測定値が測定された時刻である
ことを特徴とする請求項１または請求項２に記載の電力系統の状態推定方法。
【請求項４】
状態を測定可能な第１ノードと、状態を測定できない第２ノードと、を含む電力系統におけるノードの状態を推定する電力系統の状態推定装置であって、
前記第１ノードの状態を示す測定値を取得する第１取得部と、
第１時刻より前の第２時刻における前記第２ノードの状態を示す計算値を取得する第２取得部と、
前記第１時刻における前記測定値と、前記第２時刻における前記測定値と、の差分を、前記第２ノードに属する各ノードに按分し、前記各ノードに按分された値と、該当する前記各ノードについての前記計算値と、の和を、前記第１時刻における前記第２ノードの状態を示す値を算出するための初期値として設定する初期値設定部と、
前記初期値と、前記第１時刻における前記測定値と、を用いて最適化計算を実行し、前記第１時刻における前記第２ノードの状態を示す値を算出する最適化計算部と、
を備えることを特徴とする電力系統の状態推定装置。
【請求項５】
状態を測定可能な第１ノードと、状態を測定できない第２ノードと、を含む電力系統におけるノードの状態を推定する電力系統の状態推定装置に対して、
前記第１ノードの状態を示す測定値を取得する第１機能と、
第１時刻より前の第２時刻における前記第２ノードの状態を示す計算値を取得する第２機能と、
前記第１時刻における前記測定値と、前記第２時刻における前記測定値と、の差分を、前記第２ノードに属する各ノードに按分し、前記各ノードに按分された値と、該当する前記各ノードについての前記計算値と、の和を、前記第１時刻における前記第２ノードの状態を示す値を算出するための初期値として設定する第３機能と、
前記初期値と、前記第１時刻における前記測定値と、を用いて最適化計算を実行し、前記第１時刻における前記第２ノードの状態を示す値を算出する第４機能と、
を実行させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、電力系統の状態推定方法、電力系統の状態推定装置、及びプログラムに関する。
【背景技術】
【０００２】
状態を計器類により測定可能なノードと、状態を測定できないノードと、を含む電力系統におけるノードの状態を推定するために、与えられた、系統設備データ（例えば電源容量、負荷容量、線路データ、変圧器データ）、及び、測定データ（例えば電圧値、電流値）の下で、当該電力系統内の各地点の未知の電圧・電流等の状態を示す値（状態値）を計算により算出することが知られている。
【０００３】
このような電力系統の状態を推定するための手法として、最適化計算により状態値の修正を繰り返し行い、測定データと計算結果との偏差が最小となるように電力系統の状態値を求める手法が知られている（例えば特許文献１、２）。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００２−５１４６４号公報
【特許文献２】特開２０１０−２６３７５４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ところで、最適化計算では、初期の状態値（初期値）を与える必要がある。特許文献１においては、負荷ノードの状態値を算出するための初期値は、電源からの送り出し電圧と電流の測定値を負荷ノードの契約電力に応じて按分することで割り当てられる。特許文献２は、負荷の現在の利用率を基準として初期値を計算している。
【０００６】
ここで、最適化計算の目的関数は複雑な非線形関数であり、かつ変数である状態値の自由度が大きい。しかし、特許文献１，２では、初期値が現時点における情報のみに基づいて与えられるため、過去の傾向や蓄積されたデータとの整合性が充分にとれていない。初期値が不適切に設定されると、得られる解が実際の系統運用とは異なるおそれがある。
【課題を解決するための手段】
【０００７】
上記課題を解決するための手段の一つは、状態を測定可能な第１ノードと、状態を測定できない第２ノードと、を含む電力系統におけるノードの状態を推定する電力系統の状態推定方法であって、前記第１ノードの状態を示す測定値を取得し、第１時刻より前の第２時刻における前記第２ノードの状態を示す計算値を取得し、前記第１時刻における前記測定値と、前記第２時刻における前記測定値と、の差分を、前記第２ノードに属する各ノードに按分し、前記各ノードに按分された値と、該当する前記各ノードについての前記計算値と、の和を、前記第１時刻における前記第２ノードの状態を示す値を算出するための初期値として設定し、前記初期値と、前記第１時刻における前記測定値と、を用いて最適化計算を実行し、前記第１時刻における前記第２ノードの状態を示す値を算出する。
【０００８】
その他、本願が開示する課題、及びその解決方法は、発明を実施するための形態の欄の記載、及び図面の記載等により明らかにされる。
【発明の効果】
【０００９】
本発明によれば、電力系統の状態を精度よく推定することが可能となる。
【図面の簡単な説明】
【００１０】
【図１】本実施形態が適用される電力系統の一例を示す図である。
【図２】本実施形態に係る電力系統の状態推定装置の構成を示すブロック図である。
【図３】本実施形態に係る電力系統の状態推定装置の動作手順を示すフローチャートである。
【発明を実施するための形態】
【００１１】
本明細書および添付図面の記載により、少なくとも以下の事項が明らかとなる。
【００１２】
添付図面を参照して、本実施形態に係る電力系統の状態推定装置及び推定方法を説明する。ここで、ノードの「状態」とは、例えば、該当するノードに流れる電流、該当するノードにおける電圧、有効電力、無効電力、及び、該当するノードにおける電流と電圧との間の位相角を含むものとする。また、状態値とは、上述した状態を示す値をいう。なお、測定値とは、該当するノードに設けられた測定器によって測定された、当該ノードの状態を示す値を指すものとする。
【００１３】
［本実施形態が適用される電力系統の一例］
図１を参照して、本実施形態が適用される電力系統の一例を説明する。電力系統９０は、例えば配電系統であって、図１に示されるように、ノードｎ０〜ｎ７、負荷Ｌ１〜Ｌ７、及び測定器９５を含む。
【００１４】
ノードｎ０は、下流側のノードに電力を供給する送り出しノードである。ノードｎ０には、ノードｎ０の状態を測定するための測定器９５が設けられている。測定器９５は、予め設定された時間間隔で、例えば、ノードｎ０から下流側に送り出される電流の絶対値｜Ｉ｜、ノードｎ０における電圧の絶対値｜Ｖ｜、ノードｎ０における電流と電圧との間の位相角θ、ノードｎ０における有効電力Ｐ及び無効電力Ｑを測定し、測定結果を状態推定装置１に出力する。
【００１５】
ノードｎ０の下流側には、ノードｎ１〜ｎ７が設けられている。具体的には、ノードｎ０にノードｎ１が接続され、ノードｎ１にはノード２が接続され、ノードｎ２にはノードｎ３が接続され、ノードｎ３にはノードｎ４及びノードｎ６が接続され、ノードｎ４にはノードｎ５が接続され、ノードｎ６にはノードｎ７が接続されている。
【００１６】
ノードｎ１〜ｎ７には、ノードｎ３を除き、負荷Ｌ１〜Ｌ７がそれぞれ接続されている。負荷Ｌ１〜Ｌ７は、ノードｎ１〜ｎ７のそれぞれに接続された複数の電力負荷を包括的に示しているものとする。ここで、ノードｎ１〜ｎ７から負荷Ｌ１〜Ｌ７に流れる電流の絶対値を｜Ｉ１｜〜｜Ｉ７｜とする。ただし、本実施形態では、ノードｎ１〜ｎ７の状態を測定するための測定器は設けられていないので、電流｜Ｉ１｜〜｜Ｉ７｜を含むノードｎ１〜ｎ７の状態を直接的に知ることはできない。
【００１７】
このように、本実施形態が適用される電力系統９０は、状態を測定可能なノードｎ０（第１ノード）と、状態を測定できないノードｎ１〜ｎ７（第２ノード）と、を含む。状態推定装置１は、このような電力系統９０の各ノードにおける状態値を算出する装置である。なお、図１に示された電力系統９０は、１カ所の送り出しノードｎ０と、そのノードｎ０に接続された負荷側のノードｎ１〜ｎ７と、から構成されているが、電力系統９０は、複数の送り出しノードと、それら送り出しノードのそれぞれに接続された負荷ノードと、から構成されてもよい。
【００１８】
［最適化計算及び初期値の設定］
状態推定装置１は、各ノードの状態値の算出のために最適化計算を実行するところ、かかる最適化計算のためには初期値の設定が必要である。本実施形態では、以下のように初期値が設定されたうえで、最適化計算が実行される。
【００１９】
＜初期値の設定＞
本実施形態では、初期値は、現在（第１時刻）の測定値、過去の所定の時刻（第２時刻）における測定値、及び、当該過去の所定の時刻における状態値に基づいて設定される。かかる「過去の所定の時刻」は、例えば次の時刻から適宜選択されてよい。
−  前回の状態値計算が実行された時刻。例えば、状態値の算出が１時間毎に実施される場合には、１時間前に測定された測定値、及び、１時間前に算出された状態値が採用される。なお、状態値の計算は、例えば１５分間隔のように、任意の時間間隔で実行されてよい。
−  前日の同時刻。例えば、最新の状態値の算出が午前１１時に実行される場合、前日の午前１１時に測定された測定値、及び、その時刻に算出された状態値が採用される。
−  最新の測定値と最も近似した傾向を示す過去の測定値が測定された時刻。例えば、最新の測定値との偏差が最小となるような過去の測定値が測定された時刻における測定値、及びその測定値に対応して算出された状態値が採用される。なお、過去の測定値とその測定値に対応する状態値との間には、整合が取れているものとする。
【００２０】
次いで、選択された過去の測定値とその測定値に対応する状態値との差分を計算する。そして、かかる差分を各ノードにおける過去の状態値に配分した値を、各ノードにおける状態値を算出するための初期値として設定する。かかる差分の配分の一例としては、負荷全体の容量（例えば契約容量や設備容量）に対する、負荷側の各ノードの容量の割合に応じて、上記差分を按分することが挙げられる。
【００２１】
上述した初期値の設定手法を、図１に例示される電力系統９０に則して具体的に説明する。まず、ノードｎｉ（１≦ｉ≦７）を流れる電流の初期値｜Ｉｉ｜（１≦ｉ≦７）を、以下の（式１）で計算する。
【数１】
0006642116.tif 000002
ただし、｜Ｉｉｂａｓｅ｜は、選択された過去の時刻においてノードｎｉを流れる電流の絶対値（状態値）、｜ΔＩｉ｜は、次の（式２）を用いて算出される、送り出しノードｎ０における測定値と状態値と間の差分の、ノードｎｉへの配分である。
【数２】
0006642116.tif 000003
ここで、｜Ｉｂａｓｅ｜は、選択された過去の時刻における送り出しノードｎ０での電流の絶対値（状態値）、Ｌｉはノードｎｉに接続された負荷の契約容量、ｎは状態推定の対象となるノード数である。図１の電力系統９０ではｎ＝７である。
【００２２】
このように、本実施形態では、選択された過去の時刻における送り出しノードｎ０での測定値の、対応する状態値からの変化量（｜Ｉ｜−｜Ｉｂａｓｅ｜）を、負荷の契約容量に応じて各ノードの電流の絶対値（状態値）｜Ｉｉｂａｓｅ｜に配分することで、初期値｜Ｉｉ｜を計算する。
【００２３】
上述した（式１）及び（式２）に基づいて計算された初期値を用いると、負荷側のノードｎｉに接続された負荷Ｌｉ（１≦ｉ≦７）における有効電力Ｐｉおよび無効電力Ｑｉは、以下の（式３）及び（式４）で計算される。
【数３】
0006642116.tif 000004
【数４】
0006642116.tif 000005
ここで、｜Ｖｉ｜はノードｎｉにおける電圧絶対値の初期値であり、θｉはノードｎｉにおける電圧Ｖｉと電流Ｉｉとの間の位相角の初期値である。初期値｜Ｖｉ｜としては、例えば、前回の状態値計算が実施された時刻における電圧の絶対値や、送り出しノードｎ０の電圧｜Ｖ０｜が使用されてよい。また、初期値θｉは任意に設定されてよい。
【００２４】
このようにしてノードｎｉに接続された負荷Ｌｉにおける有効電力Ｐｉおよび無効電力Ｑｉが設定された後、次の（式５）及び（式６）に示される潮流方程式を用いて潮流計算を行う。これにより、ノードｎｉにおける電圧絶対値｜Ｖｉ｜、及び位相角θｉが得られる。
【数５】
0006642116.tif 000006
【数６】
0006642116.tif 000007
ここで、｜Ｙｉ，ｊ｜は、ノードｎｉとノードｎｊ（１≦ｊ≦７）との間のアドミタンスの絶対値、δｉ，ｊはノードｎｉとノードｎｊとの間のインピーダンスの極座標形式における偏角である。また、＊は複素共役を、Ｒｅ｛  ｝は括弧内の複素数の実部を、Ｉｍ｛
｝は括弧内の複素数の虚部を、それぞれ示している。
【００２５】
上記（式５）及び（式６）を用いてノードｎｉにおける電圧の絶対値｜Ｖｉ｜を算出することで、上記（式３）及び（式４）からノードｎｉの電流の絶対値｜Ｉｉ｜を得ることができる。
【００２６】
そして、得られた電流の絶対値｜Ｉｉ｜及び電圧の絶対値｜Ｖｉ｜に基づいて、送り出しノードｎ０における状態値を算出することができる。なお、送り出しノードが複数存在する場合には、送り出しノード及びその下流の負荷ノードの組合せ毎に、上述した状態値の計算を行えばよい。
【００２７】
＜最適化計算＞
上記のようにして算出した初期値を用いて最適化計算を行い、送り出しノードｎ０における測定値と計算値（状態値）との差分が最小となるように、負荷側のノードｎｉにおける状態値を決定する。最適化計算の例を以下に示す。ただし、以下の例では、送り出しノードが複数存在するものとしている。
【００２８】
送り出しノードにおける状態値と測定値の差を評価関数とし、次の（式７）で示す。
【数７】
0006642116.tif 000008
ここで、ｆは評価関数、ｘは操作変数、ｋは送り出しノードの数（図１の電力系統９０ではｋ＝１）、Ｖｍ，ｉは送り出しノードｉにおける電圧の測定値、Ｖｉ（ｘ）は送り出しノードｉにおける電圧の計算値、Ｉｍ，ｉは送り出しノードｉを流れる電流の測定値、Ｉｉ（ｘ）は送り出しノードｉを流れる電流の計算値、θｍ，ｉは送り出しノードｉにおけ
る位相角、θｉ（ｘ）は送り出しノードｉにおける位相角の計算値、ｗ１は電圧絶対値の偏差に関する重み係数、ｗ２は電流絶対値の偏差に関する重み係数、ｗ３は位相角に関する重み係数である。本実施形態では、操作変数ｘの一例として、負荷側のノードｎｉに接続された負荷Ｌｉにおける有効電力Ｐｉおよび無効電力Ｑｉが用いられている。
【００２９】
そして、上記（式７）で示される評価関数の値が最小となるような操作変数ｘを計算する。計算の際、上記（式７）に、電圧、電流、および状態値に対して制約条件を与えてもよい。
【００３０】
この計算に最適化手法を適用する。最適化手法としては、例えば逐次二次計画法や発見的手法であるParticle Swarm Optimization(PSO)等を適用することができる。
【００３１】
［状態推定装置の機能構成］
上述した最適化計算を実行するための状態推定装置１の機能構成を、図２を参照しつつ説明する。状態推定装置１は、図２に示されるように、入力部１０、記憶部２０、初期値設定部３０、最適化計算部４０、及び出力部５０を含む。
【００３２】
入力部１０は、測定器９５との間のインターフェイスであるとともに、ユーザとの間のインターフェイスでもある。入力部１０は、測定器９５から出力された測定データを受信し、記憶部２０に記憶する。入力部１０はまた、ユーザ又は外部サーバから電力系統９０に関するデータを受け付けて、記憶部２０に記憶する。入力部１０は更に、計算を実行する時間間隔や出力方式などに関するユーザの指令を受け付けてもよい。
【００３３】
記憶部２０は、各種データを記憶するとともに、状態推定装置１の各機能を実行するためのプログラムを記憶している。かかる記憶部２０は、系統データ記憶部２１、測定データ記憶部２２、基準データ記憶部２３、初期値記憶部２４、及び状態値記憶部２５を含む。
【００３４】
系統データ記憶部２１は、電力系統９０を構成するノード及びブランチに関するデータを記憶する。ノードに関するデータは、例えば各ノードの種別である。ノードは、基準ノード（無限大母線）、ＰＱ指定ノード、ＰＶ指定ノードに分類されるので、基準ノードでは電圧の絶対値および位相を指定し、ＰＱ指定ノードではノードの有効電力および無効電力を指定し、ＰＶ指定ノードでは有効電力および電圧絶対値を指定することになる。また、各ノードの契約容量や設備容量も系統データ記憶部２１に記憶される。ブランチに関するデータは、例えばノード接続情報、各ノード間のインピーダンス、変圧器の有無である。
【００３５】
測定データ記憶部２２は、入力部１０が測定器９５から受信した、所定の時間間隔毎の測定データを記憶する。基準データ記憶部２３は、後述する状態値記憶部２５に記憶された所定の時間間隔毎の状態値のうち、選択された過去の時刻における状態値を記憶する。
【００３６】
初期値記憶部２４は、上述した最適化計算に利用される初期値を記憶する。状態値記憶部２５は、最適化計算部４０が所定の時間間隔毎に算出した状態値を記憶する。状態値記憶部２５は、測定データ記憶部２２が記憶している過去分の測定データに対応する状態値を記憶していることが望ましい。
【００３７】
なお、測定データ記憶部２２は、送り出しノードｎ０（第１ノード）の状態を示す測定値を取得する第１取得部に相当すると見ることができる。また、基準データ記憶部２３は、過去の時刻における負荷側のノードｎ１〜ｎ７の状態を示す計算値（状態値）を取得する第２取得部に相当すると見ることができる。あるいは、後述する初期値設定部３０を、
第１取得部及び第２取得部に相当するものと見てもよい。また、記憶部２０の全部又は一部は、状態推定装置１の外部のサーバに設けられてもよい。
【００３８】
初期値設定部３０は、先に手順を詳しく説明したように、現在（第１時刻）における送り出しノードｎ０での測定値と、選択された過去の時刻（第２時刻）における送り出しノードｎ０での測定値と、選択された過去の時刻における負荷側のノードｎ１〜ｎ７での状態値（計算値）と、に基づいて、現在におけるノードｎ１〜ｎ７の状態値を算出するための初期値を設定する。
【００３９】
最適化計算部４０は、先に手順を詳しく説明したように、初期値設定部３０において算出された初期値と、現在における送り出しノードｎ０での測定値と、を用いて最適化計算を実行し、現在におけるノードｎ１〜ｎ７の状態値を算出する。最適化計算部４０は、計算結果を出力部５０及び記憶部２０に出力する。
【００４０】
出力部５０は、最適化計算部４０における計算結果を、例えば表示装置、プリンタ、外部装置などに出力する。
【００４１】
このような機能構成を有する状態推定装置１は、ＣＰＵ、ＲＡＭ、及びＲＯＭを備えたコンピュータによって実現することができる。上述した初期値設定部３０及び最適化計算部４０の各機能はＣＰＵによって実行される。
【００４２】
［状態推定装置の動作］
上述した状態推定装置１の動作、つまり状態推定方法を、図３を参照しつつ説明する。
【００４３】
まずステップＳ１において、ノード、ブランチに関するデータが系統データ記憶部２１から読み込まれる。
【００４４】
次いで、ステップＳ２において、測定データが測定データ記憶部２２から読み込まれる。測定データは、例えば送り出しノードｎ０において測定された電圧絶対値、電流絶対値、及び力率である。ここ言う測定データは、上述したとおり、最新の測定データと、選択された過去の時刻における測定データと、を含む。
【００４５】
そして、ステップＳ３において、選択された過去の時刻における状態値が基準データ記憶部２３から読み出される。
【００４６】
そして、ステップＳ４において、ステップＳ１〜Ｓ３で読みだされた系統データ、測定データ及び状態値を用いて初期値を設定する。かかるステップＳ４は、初期値設定部３０において実行される。
【００４７】
ステップＳ５において、ステップＳ２で読みだされた最新の測定データと、ステップＳ４で設定された初期値と、を用いて最適化計算を実行し、状態値を更新する。かかるステップＳ５は、最適化計算部４０において実行される。
【００４８】
そして、ステップＳ６において、最適化計算で得られた状態値を解として出力し、一連の処理が終了する。
【００４９】
［本実施形態の効果］
上記のとおり、本実施形態では、過去の測定値及び状態値に基づいて状態推定計算の初期値を設定するため、初期値を電力系統９０の実際の運転状態に近づけることができる。したがって、電力系統９０の状態を精度よく算出することができる。また、状態推定計算
の収束性が向上するとともに、計算時間を短縮することが可能となる。
【００５０】
［まとめ］
以上説明したように、状態を測定可能な送り出しノードｎ０と、状態を測定できない負荷側のノードｎ１〜ｎ７と、を含む電力系統９０におけるノードの状態を推定する電力系統の状態推定方法は、送り出しノードｎ０の状態を示す測定値を取得し、過去の時刻におけるノードｎ１〜ｎ７の状態を示す計算値（状態値）を取得し、現在における測定値と、当該過去の時刻における測定値と、前記計算値（状態値）と、に基づいて、現在におけるノードｎ１〜ｎ７の状態を示す値（状態値）を算出するための初期値を設定し、その初期値と、現在における測定値と、を用いて最適化計算を実行し、現在におけるノードｎ１〜ｎ７の状態を示す値（状態値）を算出する。かかる実施形態によれば、最適化計算のための初期値を設定するに際して、現在における測定値のほかに、過去の測定値及び計算値（状態値）を考慮しているので、過去の系統運用を踏まえた初期値の設定が可能である。したがって、電力系統９０の状態を精度よく推定することができる状態推定方法を提供することができる。また、かかる状態推定方法では、最適化計算における解の収束性が向上するとともに、計算時間の短縮化を図ることができる。
【００５１】
また、現在における測定値と過去の時刻における測定値との差分を、ノードｎ１〜ｎ７の各ノードに按分し、各ノードに按分された値と、該当する各ノードについての計算値と、の和を、初期値として設定してもよい。例えば、上記差分は、各ノードに接続された負荷の容量に応じて按分されてもよい。かかる実施形態によれば、最新の測定値と過去の測定値との差分を各ノードに配分することで、過去の運用実績に則した適切な初期値を設定することが可能となる。その結果、計算の精度、解の収束性、計算時間の短縮化の更なる向上が期待される。
【００５２】
また、計算値（状態値）として、現在における測定値との間の誤差が最小となる計算値が選択されてもよい。かかる実施形態によれば、現在の系統の状態と近似する状態を示す過去のデータが用いられるので、初期値を適切に設定することが可能となる。このことは、電力系統９０の状態の更に精度よく推定することを可能とする。また、解の収束性の更なる向上及び計算時間の更なる短縮が可能となる。
【００５３】
あるいは、状態を測定可能な送り出しノードｎ０と、状態を測定できない負荷側のノードｎ１〜ｎ７と、を含む電力系統９０におけるノードの状態を推定する電力系統の状態推定装置１は、送り出しノードｎ０の状態を示す測定値を取得する測定データ記憶部２２と、過去の時刻におけるノードｎ１〜ｎ７の状態を示す計算値（状態値）を取得する基準データ記憶部２３と、現在における測定値と、過去の時刻における測定値と、上記計算値（状態値）と、に基づいて、現在におけるノードｎ１〜ｎ７の状態を示す値（状態値）を算出するための初期値を設定する初期値設定部３０と、前記初期値と、現在における測定値と、を用いて最適化計算を実行し、現在におけるノードｎ１〜ｎ７の状態を示す値（状態値）を算出する最適化計算部４０と、を備える。かかる実施形態によれば、最適化計算のための初期値を設定するに際して、現在における測定値のほかに、過去の測定値及び計算値（状態値）を考慮しているので、過去の系統運用を踏まえた適切な初期値の設定が可能である。したがって、電力系統９０の状態を精度よく推定することができる状態推定装置を提供することができる。また、かかる状態推定装置では、最適化計算における解の収束性が向上するとともに、計算時間の短縮化を図ることができる。
【００５４】
プログラムは、状態を測定可能な送り出しノードｎ０と、状態を測定できない負荷側のノードｎ１〜ｎ７と、を含む電力系統９０におけるノードの状態を推定する電力系統の状態推定装置１に対して、送り出しノードｎ０の状態を示す測定値を取得する第１機能と、過去の時刻におけるノードｎ１〜ｎ７の状態を示す計算値（状態値）を取得する第２機能
と、現在における測定値と、過去の時刻における測定値と、前記計算値（状態値）と、に基づいて、現在におけるノードｎ１〜ｎ７の状態を示す値（状態値）を算出するための初期値を設定する第３機能と、前記初期値と、現在における測定値と、を用いて最適化計算を実行し、現在におけるノードｎ１〜ｎ７の状態を示す値（状態値）を算出する第４機能と、を実行させる。かかる実施形態によれば、最適化計算のための初期値を設定するに際して、現在における測定値のほかに、過去の測定値及び計算値（状態値）を考慮しているので、過去の系統運用を踏まえた適切な初期値の設定が可能である。したがって、状態推定装置に対して電力系統９０の状態を精度よく推定させることが可能なプログラムを提供することができる。また、かかるプログラムの実行により、最適化計算における解の収束性が向上するとともに、計算時間の短縮化を図ることができる。
【００５５】
なお、上述した実施の形態は本発明の理解を容易にするためのものであり、本発明を限定して解釈するためのものではない。本発明はその趣旨を逸脱することなく変更、改良され得るとともに、本発明にはその等価物も含まれる。
【符号の説明】
【００５６】
１  状態推定装置
１０  入力部
２０  記憶部
３０  初期値設定部
４０  最適化計算部
５０  出力部
９０  電力系統
９５  測定器
ｎ０〜ｎ７  ノード
Ｌ１〜Ｌ７  負荷
【図１】
0006642116.tif 000009
【図２】
0006642116.tif 000010
【図３】
0006642116.tif 000011