(57)【要約】【課題】  配電系統の限られた電圧・電流観測値から配電系統の潮流状態を精度良く推定する。【解決手段】  配電系統状態推定装置を配電系統の潮流状態を推定する配電系統状態推定装置であって、配電系統の計測値を格納する計測データベースと、配電系統に接続される負荷または発電量の上下限値，日変化パターンまたは線路インピーダンス修正量上下限値を格納する設備データベースと、前記潮流計算データベースから読み出した配電系統の計測値及び前記設備データベースから読み出した配電系統に接続される負荷または発電量の上下限値，日変化パターンまたは線路インピーダンス修正量上下限値に基づいて前記配電系統状態の推定値を算出する配電系統状態推定部と、前記配電系統状態推定部が算出した配電系統状態の推定値を表示する表示部とで構成する。【選択図】図１
2006087177.tif 000002
【特許請求の範囲】
【請求項１】
配電系統の潮流状態を推定する配電系統状態推定装置であって、配電系統の計測値を格納する計測データベースと、配電系統に接続される負荷の上下限値，発電量の上下限値，日変化パターンまたは線路インピーダンス修正量上下限値を格納する設備データベースと、前記計測データベースから読み出した配電系統の計測値及び前記設備データベースから読み出した配電系統に接続される負荷または発電量の上下限値，日変化パターンまたは線路インピーダンス修正量上下限値に基づいて前記配電系統状態の推定値を算出する配電系統状態推定部と、前記配電系統状態推定部が算出した配電系統状態の推定値を表示する表示部とを備えることを特徴とする配電系統の状態推定装置。
【請求項２】
請求項１において、負荷，発電量またはインピーダンスの変化に対する潮流状態の変化量を表す感度係数を算出する感度係数算出部を備え、前記配電系統状態推定部は前記感度係数算出部が算出した感度係数に基づいて配電系統状態の推定値を算出することを特徴とする配電系統の状態推定装置。
【請求項３】
請求項１または２において、前記感度係数算出部が算出する感度係数が負荷，発電量，インピーダンス変化に対する線路有効電力，線路無効電力，線路電流，負荷有効電力，発電量有効電力，負荷無効電力，発電量無効電力または母線電圧のいずれかの変化量であることを特徴とする配電系統の状態推定装置。
【請求項４】
請求項１から３のいずれかにおいて、前記設備データベースは、状態推定計算に用いる測定データのリスト，修正量を求める負荷のリストまたは修正量を求める発電設備のリストを有することを特徴とする配電系統の状態推定装置。
【請求項５】
請求項１から４のいずれかにおいて、前記設備データベースは、状態推定計算に用いる測定データ，修正量を求める負荷または修正量を求める発電設備を入力するユーザインターフェイスを有することを特徴とする配電系統の状態推定装置。
【請求項６】
請求項１から５のいずれかにおいて、前記状態推定結果データを、配電系統の測定点，ノード，線路，負荷または発電設備と対応させて画面に表示する表示部を有することを特徴とする配電系統の状態推定装置。
【請求項７】
配電系統の潮流状態を推定する配電系統状態推定方法であって、配電系統の計測値を計測データベースに記憶する手順と、配電系統に接続される負荷の上下限値，発電量の上下限値，日変化パターンまたは線路インピーダンス修正量上下限値を設備データベースに記憶する手順と、前記計測データベースから読み出した配電系統の計測値及び前記設備データベースから読み出した配電系統に接続される負荷上下限値，発電量の上下限値，日変化パターンまたは線路インピーダンス修正量上下限値に基づいて前記配電系統状態の推定値を算出する手順と、前記配電系統状態推定部が算出した配電系統状態の推定値を表示部に表示させる手順とを備えることを特徴とする配電系統状態推定方法。
【請求項８】
コンピュータを配電系統の計測値を格納する計測データベース，配電系統に接続される負荷または発電量の上下限値，日変化パターンまたは線路インピーダンス修正量上下限値を格納する設備データベース，前記潮流計算データベースから読み出した配電系統の計測値及び前記設備データベースから読み出した配電系統に接続される負荷または発電量の上下限値，日変化パターンまたは線路インピーダンス修正量上下限値に基づいて前記配電系統状態の推定値を算出する配電系統状態推定部，前記配電系統状態推定部が算出した配電系統状態の推定値を表示する表示部として機能させる配電系統状態推定プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、観測情報と設備情報を元に配電系統の電圧，電流等の潮流状態を推定する潮流状態推定装置，方法及びプログラムに関するものである。
【背景技術】
【０００２】
配電系統への分散型電源導入時の電圧維持や電力会社の設備新設・保守の省力化の要求に対し、配電系統の実状を把握した電圧・電流管理の必要性が増している。電力系統の状態推定技術は、特に高圧送電系統において各変電所（ノード）の測定値が十分に(冗長に)えられる場合を中心に技術開発が進められてきた。例えば、特公平７−１６２８９号公報に示されるように、電力系統の複数の地点で観測される潮流の観測値に重みを持たせ、重み付き最小二乗法によって系統の状態を推定する方法が示されている。
【０００３】
一方、測定情報が十分得られない配電系統の潮流状態推定方法については、これまでに、特開２００３−７９０７１号公報に配電変電所の送り出し電流からフィーダ内各区間の電流を推定する技術が示されている。これは、配電系統の送り出し電流を測定し、各区間の負荷特性データベースに基づいて各区間電流として配分し、各区間の負荷電流や電力を推定する方法である。
【０００４】
【特許文献１】特公平７−１６２８９号公報
【特許文献２】特開２００３−７９０７１号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
前述の方法では、観測値が十分に得られない場合、潮流状態を推定することが困難であったり、また配電系統内の各地点の電圧情報が得られた場合にその情報を活用して推定精度を向上させることが困難である課題があった。本発明はこれらの課題を解決し、配電系統の限られた電圧・電流観測値から配電系統の潮流状態を精度良く推定する状態推定装置，方法及びプログラムを提供することにある。
【課題を解決するための手段】
【０００６】
本発明の特徴は、配電系統の潮流状態を推定する配電系統状態推定装置であって、配電系統の計測値を格納する計測データベースと、配電系統に接続される負荷または発電量の上下限値，日変化パターンまたは線路インピーダンス修正量上下限値を格納する設備データベースと、潮流計算データベースから読み出した配電系統の計測値及び設備データベースから読み出した配電系統に接続される負荷または発電量の上下限値，日変化パターンまたは線路インピーダンス修正量上下限値に基づいて配電系統状態の推定値を算出する配電系統状態推定部と、配電系統状態推定部が算出した配電系統状態の推定値を表示する表示部とを備えた点にある。
【発明の効果】
【０００７】
本発明によれば、配電系統に設置された限られたセンサからの電圧・電流等の観測値から配電系統の負荷・発電電力，線路潮流，母線電圧などの潮流状態や線路インピーダンスを、得られた測定情報に応じて精度良く推定することができる。
【発明を実施するための最良の形態】
【０００８】
以下、本発明の実施例を説明する。
【実施例１】
【０００９】
図１は、本発明の一実施例を示す図である。図１は配電系統１００と、状態推定装置
１０とそれらを結ぶ通信ネットワーク２００の構成を示した図である。配電系統１００は、配電変電所１１０とノード（母線）１３０およびそれらを接続する線路１４０，ノード１３０に接続される負荷１５０や発電設備１６０，配電線路に設置されるセンサ１７０で構成される。配電変電所１１０は配電変電所母線１２０や、センサ１７０を含む。センサは、線路の電流，電流力率，有効電力Ｐ，無効電力Ｑ，ノード電圧Ｖを測定し、通信ネットワーク２００を介して状態推定装置１０に情報を送る。
【００１０】
次に状態推定装置１０の構成について説明する。表示装置１１，キーボードやマウス等の入力手段１２，コンピュータ（ＣＰＵ）１３，通信手段１４，ＲＡＭ１５，各種データベース（潮流計算データベース２１，計測データベース２２，設備データベース２３，状態推定結果データベース２４，プログラムデータベース２５）がバス線３０に接続されている。コンピュータ（ＣＰＵ）１３は、計算プログラムを実行して表示すべき画像データの指示，各種データベース内のデータの検索等を行う。ＲＡＭ１５は表示用の画像データ，潮流計算結果，最適化計算結果，状態推定結果データベース２４等の計算結果データを一旦格納するメモリーであり、ＣＰＵ１３によって必要な画像データを生成して表示装置１１（例えば表示ディスプレイ画面）に表示する。
【００１１】
状態推定装置内のメモリには大きく分けて５つのデータベースが格納される。潮流計算データベース２１には、線路１４０のインピーダンスを示す線路定数Ｚ(＝Ｒ＋ｊＸ)，負荷・発電量，系統の線路やノードの接続状況を表す系統構成データが記憶されている。計測データベース２２には、配電系統１００内のセンサ１７０で計測された線路の電流，電流力率，有効電力Ｐ，無効電力Ｑ，負荷や発電量，ノード電圧Ｖなどの情報が格納される。これらのデータは、通信線路２１０や通信ネットワークや状態推定装置１０の通信手段１４を介して伝送される。設備データベース２３は、配電系統内の負荷１５０や発電設備１６０の有効・無効電力の上下限や日変化パターンや、配電系統内の線路１４０の線路定数Ｚ(＝Ｒ＋ｊＸ)の上下限のデータを格納する。状態推定結果データベース２４は、状態推定計算の結果である線路の有効電力，無効電力，電流，電流力率，負荷や発電量，ノード電圧，線路定数などの計算結果を格納する。プログラムデータベース２５は、計算プログラムである潮流計算プログラム，状態推定計算プログラム，最適化計算プログラムを格納する。これらのプログラムは必要に応じてＣＰＵ１３に読み出され、計算実行される。
【００１２】
次に図２，図３を用いて設備データベース２３に格納されるデータについて説明する。図２は負荷の上下限，日変化パターンを示す説明図である。横軸は時間(１日分の２４時間)を表し、縦軸は負荷の有効電力量を示す。ＰＬmax２１１は各時間帯の負荷最大値を、ＰＬmin２１２ は同様に負荷最小値を表す。これらのデータは、各時間帯の負荷の取りうる最大値，最小値と見なすことも可能であり、後述の状態推定計算アルゴリズムで上下限制約値として扱われる。ＰＬａ２１０は負荷の有効電力の日変化パターンを示すグラフであり、各負荷毎に設けられる。負荷の有効電力が測定されず不明である場合は、この日変化パターンの該当時刻ｔの有効電力値ＰＬａ(ｔ)を疑似測定値として用いることができる。また、状態推定計算の過程で、負荷の有効電力を修正する必要が生じた場合の上下限は、時刻ｔではそれぞれＰＬmax(ｔ)，ＰＬmin(ｔ)とすることで、推定される負荷量を上下限内に納めることができる。
【００１３】
図３は発電の上下限，日変化パターンを示す説明図である。図２と同様に横軸は時間
（１日分の２４時間）を表し、縦軸は発電の有効電力量を示す。ＰＧmax２２１ は各時間帯の発電最大値を、ＰＧmin２２２ は同様に発電最小値を表す。これらのデータは、各時間帯の発電の取りうる最大値，最小値と見なすことも可能であり、後述の状態推定計算アルゴリズムで上下限制約値として扱われる。ＰＧａ２２０は発電の有効電力の日変化パターンを示すグラフであり、各発電設備毎に設けられる。発電の有効電力が測定されず不明である場合は、この日変化パターンの該当時刻ｔの有効電力値ＰＧａ(ｔ)を疑似測定値として用いることができる。また、状態推定計算の過程で、発電の有効電力を修正する必要が生じた場合の上下限は、時刻ｔではそれぞれＰＧmax(ｔ)，ＰＧmin(ｔ)とすることで、推定される発電量を上下限内に納めることができる。ここでは、有効電力に関するデータについて示したが、無効電力のデータを用意することで、無効電力の状態推定に用いることができる。
【００１４】
次に、状態推定装置１０の計算処理内容について説明する。図４は状態推定処理アルゴリズムを示すフローチャートである。ここでは、ある時刻において有効・無効電力量が不明である負荷・発電量を推定し、そのデータを用いて潮流計算を行うことで、系統内各所の電圧，電流，有効無効潮流を計算する、状態推定アルゴリズムの例を示している。以下、処理の流れを説明する。
【００１５】
まず、ステップＳ１で系統条件の設定を行う。ここでは、潮流計算に必要となる線路インピーダンスＺ(＝Ｒ＋ｊＸ)，負荷・発電量Ｐ，Ｑ初期値を、潮流計算データベース２１や設備データベース２３や入力手段１２のユーザ入力からＲＡＭ１５に読み出す。ステップＳ２では負荷・発電量の上下限制約を設定する。上下限データは設備データベース２３から読み出すか、ユーザ入力から読み出す。ここで、図２，図３を用いた説明における上下限制約を用いてもよい。また、負荷・発電力を修正する対象であるか否かについても、予め設備データベース２３に設定しておくか、ユーザが入力手段１２を介して指定しておく。ステップＳ３で、ステップＳ１で設定したデータを用いて潮流計算を行い、各ノードの電圧，線路潮流(有効・無効電力潮流，線路電流，電流力率など)を計算し、計算結果をＲＡＭ１５に格納する。ステップＳ４で、計測データベース２２に格納されたデータか、通信手段１４と通信ネットワークを介してセンサ１７０で測定されたデータをＲＡＭ１５に読み出す。ここで、測定データのなかで状態推定計算に用いるデータであるか使用しないデータであるかは、計測データベース２２に予め設定されているか、ユーザが入力手段１２を介して指定しておく。ステップＳ５では潮流計算による感度係数α，βを計算する。この感度係数は、負荷・発電量のＰ，Ｑの変化に対するノード電圧，線路潮流等の変化を表す係数であり、計算方法の詳細は後述する。ステップＳ６では、ステップＳ３で計算された潮流計算結果と、ステップＳ４で読み出した計測データとの偏差を計算する。例えば、あるノードｍの電圧については、測定値Ｖｍ′と潮流計算結果Ｖｍの偏差ΔＶｍを
ΔＶｍ＝Ｖｍ′−Ｖｍのように求めておく。電流，有効・無効電力等についても、状態推定計算に用いる測定値がある場合は同様に計算しておく。ステップＳ７で、負荷・発電量の修正量を計算する。ここで、修正量は感度係数等を用いて問題を二次計画問題として定式化を行い、それを最適化計算プログラムで解くことで計算することができる。定式化方法の詳細は後述する。最適化計算プログラムとしては一般的な二次計画法を解くアルゴリズムを用いることで、一意に解が得られる。ステップＳ８では、ステップＳ７で計算した負荷・発電量の修正量に基づいて、潮流計算入力データの負荷・発電量を修正する。修正後の入力データを用いて潮流計算を再度実施し、各ノードの電圧，線路潮流（有効・無効電力潮流，線路電流，電流力率など）を計算し、計算結果を状態推定結果データベース
２４に格納する。この潮流計算結果と、ステップＳ４で読み出した計測データとの偏差を再度計算し、偏差を求める。この偏差が基準値以上であれば、ステップＳ９からステップＳ５に戻り、再度負荷・発電量の修正量を求める。偏差が基準値以内であれば、系統の状態推定結果が得られたものとして終了する。
【００１６】
次に、図４のステップＳ７で行う最適化計算のための定式化の詳細について説明する。負荷・発電量の修正量計算問題は、二次計画法として定式化される。まず目的関数について説明する。目的関数は、潮流・負荷・発電量などの測定値と推定値の誤差に重み係数をかけたものの総和とし、（式１）で定義される。
【００１７】
【数１】
2006087177.tif 000003
【００１８】
ここで、ｙは電流Ｉ，電圧Ｖ，有効電力Ｐ，無効電力Ｑの、測定値と計算値の偏差を定義する中間変数であり、（式２）で表される。
【００１９】
【数２】
2006087177.tif 000004
【００２０】
これらの式は、二次計画問題の中の等式制約として扱われる。ここで、ΔＬｐj はノードｊに接続される負荷または発電の有効電力の修正量を、ΔＬｑj は同じく無効電力の修正量を表す。ΔＬｐj，ΔＬｑjとｙが最適化問題の変数となる。
【００２１】
次に、不等式制約として、負荷・発電量の測定値，推定値の修正量上限制約を定式化する。これは、各負荷・発電量の修正量に上限を設ける制約である。ここで、例えば、測定された値は測定誤差程度に、測定値がない負荷などは大きい値とし、誤差のしわとりを行う負荷を限定させることも可能となる。負荷・発電量の上下限制約から、負荷・発電量の修正量ΔＬｐ，ΔＬｑの上下限も容易に計算できるため、これらの上下限制約として問題を（式３ａ）〜（式３ｄ）のように定義することができる。
【００２２】
【数３】
2006087177.tif 000005
【００２３】
次に(式２)で用いた感度係数α，βについて説明する。ここで、感度係数は、図８に示すように、例えばノードｊの負荷・発電量のΔＬj を変化させた場合に、ノードｉ１３０の電圧変化量ΔＶijや線路ｋ１４０に流れる電流変化量ΔＩkj，有効電力変化量ΔＰkj，無効電量変化量ΔＱkjを示す指標である。これらの指標は、初期系統条件から負荷（発電量）を微少量変化させて潮流計算を行い、初期条件での潮流計算結果との偏差を求めることで把握することができる。負荷・発電の有効電力変化に対する感度係数α，無効電力変化に対する感度係数βの計算式は（式４ａ）〜（式４ｄ）で表される。
【００２４】
【数４】
2006087177.tif 000006
【００２５】
以上のように定式化された二次計画問題を解くことで、負荷・発電修正量ΔＬp，ΔＬqを得ることができる。
【００２６】
次に、配電線路インピーダンスの推定方法について示す。配電線路インピーダンスは通常の潮流計算では定数として扱われるが、実際は配電線路の敷設状況によっておなじ配電線種、亘長でもインピーダンスが１割以上異なることもあり得る。そこで、配電線路の計測データと前述の状態推定計算アルゴリズムによって、配電線インピーダンスの実際の定数を推定するアルゴリズムが必要となる。以下、配電線インピーダンス推定アルゴリズムについて説明する。配電線インピーダンス推定は、複数時間断面の状態推定計算結果を用いて計算を行う。まず、図５で計算に用いる各種データの計算方法を説明する。図５(ａ)はあるノードｋの負荷・発電量の一例を示したグラフである。横軸は時間を表す。実線
２３１は有効電力Ｐｆ、点線２３２は無効電力Ｑｆの状態推定計算結果を表す。ここで、時刻ｔ１，ｔ２の二つの時刻の状態推定計算結果データに着目する。時刻ｔ１の有効電力Ｐｆ(ｔ１)と時刻ｔ２の有効電力Ｐｆ(ｔ２)の偏差をΔＰｆとする。同様に、無効電力についてもＱｆ(ｔ１)とＱｆ(ｔ２)の偏差をΔＱｆと定義する。一方、図５（ｂ）は同時刻のあるノード（母線）ｉの電圧の一例を示したグラフである。実線２３３は各時刻の状態推定計算結果の電圧推定値を示している。時刻ｔ１の電圧推定値Ｖｎ(ｔ１)と時刻ｔ２の電圧推定値ΔＶｎ(ｔ２)の偏差をΔＶｎとする。また、図中の黒丸は電圧測定値を表す。時刻ｔ１の電圧測定値Ｖｎ′(ｔ１)と時刻ｔ２の電圧測定値ΔＶｎ′(ｔ２)の偏差を
ΔＶｎ′と定義する。電圧推定値の偏差ΔＶｎと電圧測定値の偏差ΔＶｎ′の差分を
ΔＶｚと定義する。
【００２７】
次に、これらの値を用いて線路インピーダンス修正量を計算するアルゴリズムについて図６を用いて説明する。
【００２８】
まず、ステップＳ２１で系統条件の設定を行う。ここでは、潮流計算に必要となる線路インピーダンスＺ（＝Ｒ＋ｊＸ），負荷・発電量Ｐ，Ｑ初期値を、潮流計算データベース２１や設備データベース２３や入力手段１２のユーザ入力からＲＡＭ１５に読み出す。ステップＳ２２では線路インピーダンスを修正する線路の選択およびインピーダンスの上下限制約を設定する。上下限データは設備データベース２３から読み出すか、ユーザ入力から読み出す。線路インピーダンスを修正する対象であるか否かについても、予め設備データベース２３に設定しておくか、ユーザが入力手段１２を介して指定しておく。ステップＳ２３で、推定に用いる時間断面ｔ１，ｔ２を選択する。これは、予め設備データベース２３に設定しておくか、ユーザが指定する。ステップＳ２４で、時間断面ｔ１，ｔ２の状態推定計算を行い、各ノードの電圧，線路潮流（有効・無効電力潮流など）を計算しておく。ステップＳ２５で、図５で説明した線路潮流変化量ΔＰｆ，ΔＱｆを計算する。ステップＳ２６では潮流計算による感度係数γｒ，γｘを計算する。この感度係数は、線路インピーダンスのｒ，ｘの変化に対するノード電圧偏差ΔＶｚの変化を表す係数である。計算方法は前述の感度係数α，βの計算と同様に複数回の潮流計算によって計算することができる。具体的には、ｔ１の時間断面での状態推定計算結果である初期条件（条件１）から線路抵抗ｒ（またはリアクタンスｘ）を微少量変化させ、なおかつ各負荷（または発電量）をΔＰｆ，ΔＱｆほど変化させた条件（条件２）で潮流計算を行い、条件１，２の電圧差分ΔＶｎ″を計算し、初期条件でのｔ１，ｔ２での状態推定結果から求めたΔＶｎ′との偏差を求めることで把握することができる。すなわち、感度γｒ＝(ΔＶｎ′−
ΔＶｎ″)／Δｒ，感度γｘ＝（ΔＶｎ′−ΔＶｎ″）／Δｘ、の式で求められる。ステップＳ２７では、図５（ｂ）で説明した、計測データを潮流計算結果データの偏差ΔＶｚ（＝ΔＶｎ−ΔＶｎ′）の計算を行う。この偏差を０とするようなインピーダンスの修正量が求める答えとなる。ステップＳ２８で、線路インピーダンスの修正量を計算する。ここで、修正量は感度係数等を用いて問題を二次計画問題として定式化を行い、それを最適化計算プログラムで解くことで計算することができる。定式化方法は前述の（式１）〜
（式４）の定式化において、変数ΔＬｐ，ΔＬｑをインピーダンス修正量Δｒ，Δｘに、感度係数α，βをγｒ，γｘと置き換えることで同様に定式化できる。ステップＳ２９では、ステップＳ２８で計算した線路インピーダンスの修正量に基づいて、潮流計算入力データの線路インピーダンスを修正する。修正後の線路インピーダンスデータとした条件１のデータから、各負荷（または発電量）をΔＰｆ，ΔＱｆほど変化させた条件（条件２）で潮流計算を再度実施し、偏差ΔＶｚを再度計算する。このΔＶｚが基準値以上であれば、ステップＳ３０からステップＳ２６に戻り、再度負荷・発電量の修正量を求める。偏差が基準値以内であれば、系統の線路インピーダンス推定結果が得られたものとして終了する。
【００２９】
次に、図７を用いて状態推定計算結果表示の一例を説明する。図７は配電系統の状態推定計算結果の表示装置１１への表示例を示す説明図である。ここでは、ディスプレイ画面への表示を考える。状態推定を行う系統図が画面に表示され、各ノード１３０に符号Ｎ１〜Ｎ５が、また各センサ１７０に符号ａ〜ｉが割り付けて表示され、負荷１５０・発電設備１６０がノードの符号で、センサ測定情報がセンサの符号で特定できるように表示される。各センサの測定結果と、状態推定計算結果は、例えば一覧表で対比して表示される。
ここで、測定値が得られている値に対しても推定値を計算している理由は、測定値に測定誤差（精度的なもの，時間ずれ的なもの）などが含まれるために、その誤差を含めてもっともつじつまのあう値を計算するために必要となる。例えば、すぐとなりの地点で測定しても、異なる値が情報として収集されるような場合があり、その対策として推定値を計算することが必要となる。ここでは、表９１に各測定点での潮流を、表９２に各ノードの負荷を、表９３に各ノードの発電量を、表９４に各線路のインピーダンスを、それぞれ測定値(または初期値)と状態推定計算結果である推定値の対比で示している。このように表示することで、ユーザに状態推定計算結果をわかりやすく伝えることが可能となる。ここでは、画面への出力例を示したが、書類等に印刷可能なフォーマットのデータとしてユーザに提供してもよい。
【００３０】
上記本発明の実施例によれば、配電系統に設置された限られたセンサからの電圧・電流等の観測値から配電系統の負荷・発電電力，線路潮流，母線電圧などの潮流状態や線路インピーダンスを、得られた測定情報に応じて精度良く推定することができる。
【００３１】
また、配電系統の母線電圧範囲などが限られた条件範囲内に保つために必要となる、負荷・発電量の範囲を求めるための設備計画支援用途に利用できる利点がある。
【００３２】
また、各時刻における潮流計算を行うためのデータを用意できるために、将来時刻に発生する系統条件変化に対する系統電圧・電流変化量を予め潮流計算によって把握することができることによって、配電系統の信頼度評価システムや配電自動化システムの機能として用いることも可能となる利点がある。
【産業上の利用可能性】
【００３３】
配電系統に設置された限られたセンサからの電圧・電流等の観測値から配電系統の負荷・発電電力，線路潮流，母線電圧などの潮流状態や線路インピーダンスを、得られた測定情報に応じて精度良く推定するための状態推定装置に適用できる。また、配電系統の母線電圧範囲などが限られた条件範囲内に保つために必要となる、負荷・発電量の範囲を求めるための設備計画支援用途にも活用できる。また、各時刻における潮流計算を行うためのデータを用意できるために、将来時刻に発生する系統条件変化に対する系統電圧・電流変化量を予め潮流計算によって把握することができるため、配電系統の信頼度評価システムや配電自動化システムの機能として用いることも可能である。
【図面の簡単な説明】
【００３４】
【図１】配電系統と状態推定装置の構成例を示した説明図である。
【図２】負荷の上下限，日変化パターンを示す説明図である。
【図３】発電の上下限，日変化パターンを示す説明図である。
【図４】状態推定処理アルゴリズムを示すフローチャートである。
【図５】負荷・発電量およびノード電圧の一例を示したグラフである。
【図６】線路インピーダンス修正量を計算するアルゴリズムを示すフローチャートである。
【図７】配電系統の状態推定計算結果の表示例を示す説明図である。
【図８】感度係数の求め方を補足する説明図である。
【符号の説明】
【００３５】
１０…状態推定装置、１１…表示装置、１２…入力手段、１３…コンピュータ、１４…通信手段、１５…ＲＡＭ、２１…潮流計算データベース、２２…計測データベース、２３…設備データベース、２４…状態推定結果データベース、２５…プログラムデータベース、３０…バス線、１００…配電系統、１１０…配電変電所、１３０…ノード、１４０…線路、１５０…負荷、１６０…発電設備、１７０…センサ、２００…通信ネットワーク、
２１０…通信線路、２２０…ＰＧａ、２２１…ＰＧmax、２２２…ＰＧmin。
【図１】
2006087177.tif 000007
【図２】
2006087177.tif 000008
【図３】
2006087177.tif 000009
【図４】
2006087177.tif 000010
【図５】
2006087177.tif 000011
【図６】
2006087177.tif 000012
【図７】
2006087177.tif 000013
【図８】
2006087177.tif 000014