JP2023010451A - 状態推定装置、状態推定装置の制御方法及びプログラム - Google Patents

Abstract

【課題】配電系統において、推定を行いたい箇所（推定点）の数に比べ測定器の数が少ない場合であっても配電系統の状態を推定する。【解決手段】配電系統において複数のノードを有して構成される区間の状態を求める状態推定装置であって、前記区間の端部における潮流測定値、及び前記複数のノードにおける有効電力測定値を取得する測定値取得部と、前記複数のノードにおける無効電力の比率を表す分布係数を取得する分布係数取得部と、前記区間における無効電力の合計値を前記分布係数を用いて前記複数のノードにおける無効電力に換算することで前記区間の状態を推定可能に構成された第１方程式に、前記潮流測定値、及び前記有効電力測定値を入力することで、前記区間における無効電力の合計値と、前記複数のノードにおける有効電力と、を前記区間の状態として求める状態推定部と、を備える。【選択図】図３

Description

本発明は、状態推定装置、状態推定装置の制御方法及びプログラムに関する。

例えば、配電系統の状態を、計測された系統情報に基づいて推定する技術が知られている（例えば、特許文献１参照）。

特開２０１８－０８５８７９号公報

しかしながら、配電系統において推定を行いたい箇所（推定点）の数に比べ測定器の数が少ない場合のように、測定データに基づき推定点の状態値を一意に決定することが困難な場合がある。

本発明はこのような課題を鑑みてなされたものであり、配電系統において推定を行いたい箇所（推定点）の数に比べ測定器の数が少ない場合であっても、配電系統の状態を推定可能とする状態推定装置、状態推定装置の制御方法及びプログラムを提供することを目的とする。

前述した課題を解決する一態様は、配電系統において複数のノードを有して構成される区間の状態を求める状態推定装置であって、前記区間の端部における潮流測定値、及び前記複数のノードにおける有効電力測定値を取得する測定値取得部と、前記複数のノードにおける無効電力の比率を表す分布係数を取得する分布係数取得部と、前記区間における無効電力の合計値を前記分布係数を用いて前記複数のノードにおける無効電力に換算することで前記区間の状態を推定可能に構成された第１方程式に、前記潮流測定値、及び前記有効電力測定値を入力することで、前記区間における無効電力の合計値と、前記複数のノードにおける有効電力と、を前記区間の状態として求める状態推定部と、を備える。

本発明は、配電系統において推定を行いたい箇所（推定点）の数に比べ測定器の数が少ない場合であっても、配電系統の状態を推定することができる。

状態推定装置の構成を示す図である。 記憶装置の構成を示す図である。 配電系統及び状態推定装置を説明するための図である。 配電系統を説明するための図である。 配電系統を説明するための図である。 状態推定装置の機能構成を示す図である。 状態推定装置で実行される処理の一例を示すフローチャートである。

本明細書及び添付図面の記載により、少なくとも以下の事項が明らかとなる。
＜＜＜状態推定装置２００の構成＞＞＞

図１は、本発明の一実施形態である状態推定装置２００の構成を示す図である。状態推定装置２００は、図３～図５等に示す配電系統１０００の状態を推定する装置であり、ＣＰＵ（Central Processing Unit）２１０、メモリ２２０、通信装置２３０、記憶装置２４０、入力装置２５０、出力装置２６０、及び記録媒体読取装置２７０を有するコンピュータである。

ＣＰＵ２１０は、メモリ２２０や記憶装置２４０に記憶された状態推定装置制御プログラム７００を実行することにより、状態推定装置２００が有する様々な機能を実現する。

メモリ２２０は、例えばＲＡＭ（Random-Access Memory）等であり、様々なプログラムやデータ等の一時的な記憶領域として用いられる。

記憶装置２４０は、ＣＰＵ２１０によって、実行または処理される各種データを格納する非一時的な（例えば不揮発性の）記憶装置である。

記憶装置２４０に状態推定装置制御プログラム７００及び系統情報テーブル６００が記憶されている様子を図２に示す。

記憶装置２４０に記憶されている状態推定装置制御プログラム７００や系統情報テーブル６００等の各種のデータがメモリ２２０に読み出されてＣＰＵ２１０によって実行あるいは処理されることにより、状態推定装置２００の各種機能が実現される。

また状態推定装置制御プログラム７００は、本実施形態に係る状態推定装置２００が有する機能を実現するためのプログラムを総称しており、例えば、状態推定装置２００上で動作するアプリケーションプログラムやＯＳ（Operating System）、種々のライブラリ等を含む。

系統情報テーブル６００は、配電系統１０００の構成や、配電系統１０００を構成する機器の電気的特性などを記録したテーブルである。

系統情報テーブル６００には、例えば、配電系統１０００を状態方程式や潮流方程式などの方程式を用いて模擬する際に必要なデータが記録されている。詳細は後述するが、状態推定装置２００は、センサ付き開閉器ＳＷの潮流の測定値と、スマートメータＳＭの有効電力測定値を取得した状態で、系統情報テーブル６００を用いて、以下に詳述する状態推定処理を実行する。状態推定処理では、区間Ｋをはじめとする配電系統１０００の状態、つまり、配電系統１０００における電圧、電流、電力等の状態値が得られる。

図１に戻って、入力装置２５０は、ユーザによるコマンドやデータの入力を受け付ける装置であり、キーボード、タッチパネルディスプレイ上でのタッチ位置を検出するタッチセンサなどの入力インタフェースを含む。

出力装置２６０は、例えばディスプレイやプリンタなどの装置である。

通信装置２３０は、ネットワーク５００を介して、他のコンピュータと各種プログラムやデータの受け渡しを行う。

記録媒体読取装置２７０は、ＳＤカードやＤＶＤ、ＣＤＲＯＭ等の記録媒体８００に記録された状態推定装置制御プログラム７００や系統情報テーブル６００等の様々なデータを読み取り、記憶装置２４０に格納する。
＜＜＜配電系統１０００の一例＞＞＞
図３～図５は、状態推定装置２００が状態推定を行う配電系統１０００の一例を示す図である。配電系統１０００は、例えば、６．６ｋＶ系の高圧系統であり、配電変電所１１００、配電線１２００、センサ付き開閉器ＳＷ、スマートメータＳＭを含む。

センサ付き開閉器ＳＷは、配電系統１０００の一つの区間Ｋの端部に配置される。つまり、センサ付き開閉器ＳＷは、配電系統１０００を複数の区間に区切る。センサ付き開閉器ＳＷは、測定点（区間Ｋの端部）の電圧、電流及び位相、または有効電力及び無効電力を計測可能な計測器である。以下では、センサ付き開閉器ＳＷによって測定された電圧電流及び位相、または有効電力および無効電力を、包括的に「潮流測定値」または「潮流の測定値」と称する。

スマートメータＳＭは、配電系統１０００に接続されている需要家が消費した有効電力を計測可能な計測器である。

需要家は、それぞれ電力の負荷となる。需要家には、配電線１２００から供給される電力を消費する設備（例えば、工場）が含まれ得るとともに、配電線１２００に対し、電力を供給するインバータ等の発電設備（不図示）も含まれる。このため、配電線１２００には、配電線１２００からの電力を消費する設備と、配電線１２００に電力を供給する発電設備と、が接続されていることになる。ここでは、消費、供給の別を問わず、包括して負荷と称する。

また、配電系統１０００の区間Ｋ内には複数のノードＮが設けられているが、ノードＮは、区間Ｋをさらに細かく区分した各範囲を表し、区間Ｋに設置されている変圧器や配電線１２００の分岐点などのような、適宜定められた配電系統１０００上の位置により特定される。そして本実施形態では、区間Ｋ内の各ノードＮには、１以上の需要家のスマートメータＳＭが集約されているものとする。

以下では、スマートメータＳＭによって測定され、各ノードＮに集約された有効電力を、包括的に「有効電力測定値」と称する。

配電変電所１１００は、送電線（不図示）から供給される電圧を変圧し、６．６ｋＶの電圧を配電線１２００へと出力する。電力は、ノードＮを介して需要家（不図示）に供給される。

なお配電系統１０００には、その他にも様々な設備やセンサ等が含まれているが、便宜上、ここでは簡素化した配電系統１０００を一例として図示している。

例えば、図４の例を参照しながら配電系統１０００について説明する。

配電系統１０００には、配電変電所１１００を起点（送り出しノード）とし、配電線１２００が放射状に接続されている。配電線１２００にはセンサ付き開閉器ＳＷ１～ＳＷ５が設置され、これらのセンサ付き開閉器ＳＷ１～ＳＷ５で分割された区間が区間Ｋ１～区間Ｋ４として定められている。

ただし、区間Ｋ４のように配電系統１０００の末端側にセンサ付き開閉器ＳＷが設置されていない場合には、送り出し側のセンサ付き開閉器ＳＷから末端ノードまでを区間として定める。

また、配電線１２００上にはノードＮを定義する。ノードＮは柱上変圧器単位や指定された需要家単位で管理されている集約単位である。複数の需要家のスマートメータＳＭからの計測値は各ノードＮごとに集約され、状態推定で用いられる。

センサ付き開閉器（ＳＷ１～ＳＷ５）では、定周期Ｔｓ１（例えば１分）で設置点の電圧、電流、位相が計測される。これらの計測値は、データベース（例えば記憶装置２４０に構築される）に格納される。

この時、電圧、電流、位相から有効電力潮流、無効電力潮流には一意に変換することができるものとし、以降、センサ付き開閉器ＳＷからの測定値は、電圧、有効電力潮流、無効電力潮流として取り扱う。

スマートメータＳＭでは、各ノード（Ｎ１～Ｎ１２）単位で、定周期Ｔｓ２（例えば３０分）で需要家の電力使用量（有効電力）が測定される。これらの測定値は、データベース（例えば記憶装置２４０に構築される）に格納される。

周期毎の電力使用量を測定周期で割ることで、有効電力量の平均値を計算することができるため、スマートメータＳＭからの測定値は有効電力平均値として取り扱う。
＜＜＜情報処理装置の機能ブロック＞＞＞

図６は、状態推定装置２００の機能ブロックを示す図である。状態推定装置２００は、測定値取得部２０１、分布係数取得部２０２、状態推定部２０３及び無効電力算出部２０４の各機能を有する。

これらの各機能は、状態推定装置２００のハードウェアによって本実施形態に係る状態推定装置制御プログラム７００が実行されることにより実現される。

測定値取得部２０１は、配電系統１０００の区間Ｋの端部における潮流測定値、及び区間Ｋ内の複数のノードＮにおける有効電力測定値を取得する。上述したように、本実施形態では測定値取得部２０１は、センサ付き開閉器ＳＷから１分毎に潮流測定値を取得し、スマートメータＳＭから３０分毎に有効電力測定値を取得する。

分布係数取得部２０２は、区間Ｋ内の複数のノードＮにおける無効電力の比率を表す分布係数を取得する。例えば分布係数取得部２０２は、これらの複数のノードＮに接続される設備に関するデータ、及び複数のノードＮにおける過去の無効電力の測定データの両方又は一方を用いて所定のデータ解析を行うことにより、分布係数を取得する。詳しくは後述する。

状態推定部２０３は、区間Ｋにおける無効電力の合計値を分布係数を用いて複数のノードＮにおける無効電力に換算することで区間Ｋの状態を推定可能に構成された第１方程式に、潮流測定値及び有効電力測定値を入力することで、区間Ｋにおける無効電力の合計値と、前記複数のノードにおける有効電力と、を区間の状態として求める。

このような態様により、状態推定装置２００は、配電系統１０００において推定を行いたい箇所（推定点）の数に比べ測定器の数が少ない場合であっても配電系統１０００の状態を推定することが可能となる。

例えば状態推定部２０３は、区間Ｋ内の各ノードＮの無効電力が分からなくても、区間Ｋの無効電力の合計値が分かっていれば、これを第１方程式に入力して解くことで、区間Ｋの状態を推定することが可能となる。

そのため、状態推定装置２００は、配電系統１０００の区間Ｋ内の各ノードＮの無効電力が得られない場合であっても、配電系統１０００の状態を推定することが可能となる。

たとえば、状態推定装置２００は、配電系統１０００の状態推定を行う際に、状態値の要素を各ノードＮの有効電力及び無効電力とはせずに、各ノードＮの有効電力と、指定された区間Ｋ内のノードの無効電力合計値と定義しているため、測定値の要素数に対する自由度が制約され、状態推定を行う際に状態値を一意に決定することができる。

また電力の線路ロスに関しては、各ノードＮの有効電力と、指定した区間Ｋの全ノードＮの無効電力合計値を状態値として取り扱い、第１方程式による測定値の計算を行うことで線路ロスが考慮された状態推定を行うことが可能となっている。

無効電力算出部２０４は、配電系統１０００の区間Ｋにおける無効電力の合計値と、分布係数と、を用いて、区間Ｋ内の複数のノードＮにおける無効電力を算出する。

この場合、状態推定部２０３は、さらに、区間Ｋ内の複数のノードＮにおける無効電力と、区間Ｋ内の複数のノードにおける有効電力と、を用いて所定の第２方程式を計算し、複数のノードＮにおける電圧値を算出することにより、区間Ｋの状態を推定するようにしてもよい。

このような態様によれば、状態推定装置２００は、配電系統１０００において推定を行いたい箇所（推定点）の状態を推定することが可能となる。

図３を参照しながらもう少し詳しく説明する。

図３は、状態推定装置２００が、センサ付き開閉器ＳＷから電圧、電流、位相に関する定周期（例えば１分）の測定値（潮流測定値）を取得し、スマートメータＳＭから電力使用量に関する定周期（例えば３０分）の測定値（有効電力測定値）を取得し、記憶装置２４０内のデータベースに格納済みの系統情報テーブル６００（系統設備データ）を取得して、第１方程式を用いて潮流計算を行うことで、事前に定義したノードＮ単位の有効電力値と、事前に定義した区間Ｋのノードの無効電力合計値を、配電系統１０００の状態として推定する様子を示している（状態推定１）。

また、状態推定装置２００は、無効電力の区間Ｋ単位での分布モデル（分布係数）を用い、区間Ｋでの無効電力合計値を各ノードＮに展開することで、各ノードＮの有効電力と無効電力を得て、第２方程式を用いて潮流計算へこれらを入力することで、各ノードＮの電圧分布を推定する（状態推定２）。

このように、本実施形態に係る状態推定装置２００は、配電系統１０００上に設置されたセンサ付き開閉器ＳＷの電圧、電流、位相と、需要家のスマートメータＳＭからの電力量を測定値とし、配電系統１０００の負荷容量や太陽光発電を含む電源容量、線路データ、変圧器データ等で構成される系統構成データ（系統情報テーブル６００）を用いて、配電系統１０００内の未測定地点の電圧、電流（潮流）を定周期で推定する。

なお状態推定部２０３は、潮流測定値、及び有効電力測定値に対して、測定周期が長いほど重みが小さくなるように定められる重み係数を乗じた上で第１方程式を計算するようにするとよい。

つまり、本実施形態では、潮流測定値は１分周期、有効電力測定値は３０分周期で取得可能な測定値であるから、状態推定部２０３は、潮流測定値よりも有効電力測定値の方が小さな重み係数を乗じてから第１方程式を計算するとよい。

このような態様により、測定周期が遅い測定値に関しては状態推定結果への影響度合いを小さくし、より高精度に状態推定を行うことが可能となる。

また状態推定部２０３は、所定期間毎（例えば１分ごと）に、最新の潮流測定値、及び有効電力測定値を用いて第１方程式を計算することにより、配電系統１０００の区間Ｋにおける無効電力の合計値と、区間Ｋ内の複数のノードＮにおける有効電力と、を繰り返し算出するようにするとよい。

このような態様により、常に配電系統１０００の最新の状態を推定することが可能となる。
＜＜＜処理の詳細＞＞＞

状態推定装置２００が実行する状態推定処理の第１の具体例を、図７の処理フローを用いて説明する。

なお状態推定装置２００は、配電系統１０００全体の区間Ｋを一括してモデル化して計算を実施するが、説明の簡易化のため、図５に示す配電系統１０００の区間Ｋにおける状態推定処理について説明する。

図５に示す配電系統１０００は、２台のセンサ付き開閉器ＳＷ、２つのノードＮ、１つの区間Ｋを有して構成され、各ノードＮおよびセンサ付き開閉器ＳＷは、配電線１２００で接続されている。

また、ノードＮ毎にスマートメータＳＭの測定値（有効電力測定値）は集約され、記憶装置２４０のデータベースに格納され、状態推定に用いられる。

以下に示す物理量は、規格化されたｐｕ単位を用いる。時刻ｔは離散化された値を用いる。

この時、時刻ｔにおける全体の測定値ｙ_ｔは以下の式（３）となる。

ただし、測定周期は、センサ付き開閉器ＳＷを１分、スマートメータＳＭは３０分とし、状態推定の演算周期を１分とする。

そのため、センサ付き開閉器ＳＷからの測定値は演算周期毎に更新されるが、スマートメータＳＭからの測定値は３０回の演算中は更新されず、一定値を用いることになる。

例えば、Ｌを以下の（６）のように作成する。

したがって、ノードＮ１とノードＮ２の無効電力は、（８）のように表される。

状態推定装置２００は、例えば配電系統１０００内の負荷設備の容量、太陽光発電設備の容量、電力契約形態等の設備に関連するデータや、過去の無効電力の測定データを用いたパラメータ調整等のデータ解析により分布係数を決定する（Ｓ１０２０）。

ただし、Ｈはｈの線形化行列である。この時、Ｆ＝HL とした場合、状態方程式Ｆ（第１方程式）は以下の（１２）となる。このようにして状態方程式Ｆは、分布係数を内在させ、そしてこれにより、区間Ｋにおける無効電力の合計値を複数のノードＮにおける無効電力に換算可能に構成される。

そのため、センサ付き開閉器ＳＷでの測定周期を状態推定の演算周期とした場合、スマートメータＳＭからの測定値は更新されない期間があり、その場合には、実際に発生し得る変動をとらえることができない。

重みについては、センサ付き開閉器ＳＷからの測定値の重みと比較し、スマートメータＳＭからの測定値の重みを小さくすることで、状態推定による状態値の推定結果に対し、スマートメータからの測定値の影響を小さくすることができる。

これにより、スマートメータＳＭからの測定値が更新されない間においても、センサ付き開閉器ＳＷからの測定値を用いて変動を詳細に状態値に反映することができる。

状態推定装置２００は、過去の測定値データを用いた状態推定のシミュレーションによる解析等を行うことで、重みの具体的な値を決定する。

例えば、状態推定装置２００は、状態推定手法として最小二乗法を適用した場合、以下の式（１５）で状態値を得る。

なお状態推定手法としては、遂次最小二乗法やカルマンフィルタなど、各種状態手法を適用することが可能である。

とし、状態値を各ノードＮ１、Ｎ２の有効電力と無効電力に展開し、各ノードＮ１、Ｎ２の有効電力と無効電力の値から、各ノードＮ１、Ｎ２の電圧（振幅、位相）を出力する潮流方程式（第２方程式）によって潮流計算を行い、各ノードＮ１、Ｎ２の電圧を計算することができる（Ｓ１０４０）。

なお本実施形態では、単一の区間Ｋの状態推定を行う場合の例を示したが、複数の区間Ｋの状態推定を行う場合でも、同様の計算により状態値を推定可能である。

次に、状態推定装置２００が行う状態推定処理の第２の具体例について説明する。

本実施形態では、センサ付き開閉器ＳＷからの測定値（潮流測定値）と、スマートメータＳＭからの測定値（有効電力測定値）と、で構成される測定値ベクトルを定義し、事前に設定した複数のノードＮの有効電力と、事前に設定した区間Ｋ内の複数のノードＮの無効電力合計値と、で構成される状態値ベクトルを定義し、状態値ベクトルを用いて測定点の値を計算する状態方程式（第１方程式）を定義する。

ただし、センサ付き開閉器ＳＷからの測定値は、実際の測定値（例えば電圧、電流、位相）を想定し、その値から電圧、有効電力潮流、無効電力潮流に変換している。また、スマートメータＳＭからの測定値は使用電力量を想定し、その値から有効電力平均値に変換している。また、スマートメータＳＭからの測定値は事前に設定したノード単位に集約され、測定値ベクトルに与えられる。

重み行列Ｗは、以下の（２２）に示すように、重みベクトルｗを対角要素にもつ行列である。

仮に重みをすべて１とした場合には、すべての測定値が均等に考慮された場合の状態値が計算される。しかし、センサ付き開閉器ＳＷとスマートメータＳＭでは測定周期が異なり、一般的にはスマートメータＳＭからの測定値の方が測定周期が遅い。測定周期が短い測定値に合わせて状態推定を計算する際、測定周期が長い測定値は、その値が更新されるまでは更新前の値を測定値として与えることになる。

また、スマートメータＳＭからの測定値は、電力量を有効電力の平均値に変換された値となっており、さらに時間遅れの要素が含まれる。そのため、センサ付き開閉器ＳＷでの測定周期に合わせて状態推定を行う場合には、スマートメータＳＭでの測定値はデータとしての信頼性は低いと考えられる。

これらを考慮した場合、重みベクトルのスマートメータＳＭに関する値を、センサ付き開閉器ＳＷに関する値と比較し小さい値を与えることで、スマートメータＳＭの測定値の状態推定への影響度合いを弱めることができ、主にノード有効電力に関する変動要素を状態値に反映することが可能となる。

状態推定装置２００は、過去データを用いた解析を行うことで重みベクトルを事前に設定し、状態推定を行う。

なお、行列Ｌには無効電力の分布を表すモデルが必要であるが、例えば過去データを用いた解析や負荷設備、太陽光発電設備に関する値を用いて作成することができる。

行列Ｌを用いるようにした場合、線形化された潮流方程式は以下の（２６）のように表される。

以上、本実施形態に係る状態推定装置２００、状態推定装置２００の制御方法、及びプログラム７００について詳細に説明したが、本実施形態によれば、配電系統１０００において推定を行いたい箇所（推定点）の数に比べ測定器の数が少ない場合であっても、配電系統１０００の状態を推定することが可能となる。

なお上述した実施の形態は本発明の理解を容易にするためのものであり、本発明を限定して解釈するためのものではない。本発明はその趣旨を逸脱することなく変更、改良され得るとともに、本発明にはその等価物も含まれる。

例えば、最小二乗法、重み付き最小二乗法、負荷修正型状態推定法といった各時間断面で独立した状態推定手法の他に、カルマンフィルタのような時系列で遂次処理を行う状態推定手法においても本発明の適用が可能である。

２００ 状態推定装置
２０１ 測定値取得部
２０２ 分布係数取得部
２０３ 状態推定部
２０４ 無効電力算出部
２１０ ＣＰＵ
２２０ メモリ
２３０ 通信装置
２４０ 記憶装置
２５０ 入力装置
２６０ 出力装置
２７０ 記録媒体読取装置
５００ ネットワーク
６００ 系統情報テーブル
７００ 状態推定装置制御プログラム
８００ 記録媒体
１０００ 配電系統
１１００ 配電変電所
１２００ 配電線
ＳＷ センサ付き開閉器
ＳＭ スマートメータ
Ｎ ノード
Ｋ 区間

Claims (7)

1. 配電系統において複数のノードを有して構成される区間の状態を求める状態推定装置であって、
   前記区間の端部における潮流測定値、及び前記複数のノードにおける有効電力測定値を取得する測定値取得部と、
   前記複数のノードにおける無効電力の比率を表す分布係数を取得する分布係数取得部と、
   前記区間における無効電力の合計値を前記分布係数を用いて前記複数のノードにおける無効電力に換算することで前記区間の状態を推定可能に構成された第１方程式に、前記潮流測定値、及び前記有効電力測定値を入力することで、前記区間における無効電力の合計値と、前記複数のノードにおける有効電力と、を前記区間の状態として求める状態推定部と、
   を備える状態推定装置。
2. 請求項１に記載の状態推定装置であって、
   前記区間における無効電力の合計値と、前記分布係数と、を用いて、前記複数のノードにおける無効電力を算出する無効電力算出部と、
   を備え、
   前記状態推定部は、さらに、前記複数のノードにおける無効電力と、前記複数のノードにおける有効電力と、を用いて所定の第２方程式を計算することで、前記複数のノードにおける電圧値を、前記区間の状態として求める、
   状態推定装置。
3. 請求項１又は２に記載の状態推定装置であって、
   前記状態推定部は、前記潮流測定値、及び前記有効電力測定値に対して、測定周期が長いほど重みが小さくなるように定められる重み係数を乗じた上で前記第１方程式を計算する、
   状態推定装置。
4. 請求項１～３のいずれかに記載の状態推定装置であって、
   前記分布係数取得部は、前記複数のノードに接続される設備に関するデータ、及び前記複数のノードにおける過去の無効電力の測定データの少なくともいずれかを用いて所定のデータ解析を行うことにより前記分布係数を取得する、
   状態推定装置。
5. 請求項１～４のいずれかに記載の状態推定装置であって、
   前記状態推定部は、所定期間毎に、最新の前記潮流測定値、及び前記有効電力測定値を用いて前記第１方程式を計算することにより、前記区間における無効電力の合計値と、前記複数のノードにおける有効電力と、を繰り返し算出する、
   状態推定装置。
6. 配電系統において複数のノードを有して構成される区間の状態推定を行う状態推定装置の制御方法であって、
   前記状態推定装置が、
   前記区間の端部における潮流測定値、及び前記複数のノードにおける有効電力測定値を取得し、
   前記複数のノードにおける無効電力の比率を表す分布係数を取得し、
   前記区間における無効電力の合計値を前記分布係数を用いて前記複数のノードにおける無効電力に換算することで前記区間の状態を推定可能に構成された第１方程式に、前記潮流測定値、及び前記有効電力測定値を入力することで、前記区間における無効電力の合計値と、前記複数のノードにおける有効電力と、を前記区間の状態として求める、
   状態推定装置の制御方法。
7. 配電系統において複数のノードを有して構成される区間の状態推定を行うコンピュータに、
   前記区間の端部における潮流測定値、及び前記複数のノードにおける有効電力測定値を取得する手順と、
   前記複数のノードにおける無効電力の比率を表す分布係数を取得する手順と、
   前記区間における無効電力の合計値を前記分布係数を用いて前記複数のノードにおける無効電力に換算することで前記区間の状態を推定可能に構成された第１方程式に、前記潮流測定値、及び前記有効電力測定値を入力することで、前記区間における無効電力の合計値と、前記複数のノードにおける有効電力と、を前記区間の状態として求める手順と、
   を実行させるためのプログラム。
   

Images