JP3620930B2

JP3620930B2 - Power system characteristic estimation apparatus and characteristic estimation method

Info

Publication number: JP3620930B2
Application number: JP18048496A
Authority: JP
Inventors: 孝浩野畑; 勇一中尾; 孝典角田; 吉村　　隆志; 哲也福永
Original assignee: Tokyo Electric Power Co Inc; Nissin Electric Co Ltd
Current assignee: Tokyo Electric Power Co Inc; Nissin Electric Co Ltd
Priority date: 1996-07-10
Filing date: 1996-07-10
Publication date: 2005-02-16
Anticipated expiration: 2016-07-10
Also published as: JPH1028327A

Description

【０００１】
【発明の属する技術分野】
この発明は、電力系統の特性を推定する装置および方法に関し、特に調相設備容量や無効電力の系統負荷特性定数の推定を容易にした電力系統の特性推定装置および特性推定方法に関する。
【０００２】
【従来の技術】
従来より、力率改善および調相のために負荷系統に調相設備が設けられているが、その容量がどれだけであるか、および系統の負荷特性の把握が、電力系統の運用上または電力系統の計画上重要なポイントの一つとなる。
【０００３】
一般に、実際の負荷系統は一次変電所の二次側母線以下が縮約して模擬され、有効電力と無効電力との関係は或る一定の方程式で関係付けられる。実際の負荷系統の典型的な等価回路を図２１のように表せば、有効電力Ｐと無効電力Ｑの関係は、Ｑ＝ＡＰ^２＋ＢＰ−Ｃで表現される。ここでＡは電力系統の線路のリアクタンスに関係する値であり、ＡＰ^２は線路のリアクタンスによる無効電力、Ｂは負荷の力率に関係する値であり、ＢＰは負荷による無効電力、Ｃは電力系統の充電容量と調相設備の容量による無効電力に相当する。従って有効電力と無効電力の直角座標平面（以下「Ｐ−Ｑ平面」という。）に、一定周期で繰り返し測定した有効電力と無効電力の対をプロットすれば、例えば図２２に示すように、理想的には上記２次方程式で示される曲線（以下「Ｐ−Ｑカーブ」という。）上に並ぶことになる。
【０００４】
そこで従来は、３相電力系統の有効電力と無効電力を分析する負荷特性記録装置を用いて得た有効電力と無効電力の対から、最小二乗法によって上記２次方程式の係数Ａ，Ｂ，Ｃをそれぞれ求めている。そして、特に係数Ｃ（以下「Ｃ項」という。）に基づいて調相設備の容量を推定し、これを例えば電圧一定化のための情報として利用して、電力系統の適正な運用を計っている。
【０００５】
【発明が解決しようとする課題】
ところが実際の負荷系統は図２３の等価回路に示すように負荷力率が変動するため、上記ＢＰ項の値が変動し、重要なＣ項の推定が困難となる場合があった。特に工場負荷系統では負荷力率の変動が激しい。ここで、工場負荷系統における有効電力と無効電力の推移を具体的に表せば図２４に示すようになり、時間経過に伴って大きく変動する。そのため、有効電力と無効電力の対をＰ−Ｑ平面上にプロットしても、図２５に示すように大きくばらつく。（この図２５に示した例では、近似曲線を２本引くことができるので、２４時間の間に調相設備が投入／遮断されているものと予想される。）従ってこのようなデータを基に最小二乗法によって上記２次方程式の係数を求めても、信頼性の低い値しか求められないことになる。
【０００６】
また、従来は一定周期で求めた有効電力と無効電力のデータをそのまま用いているが、一般に有効電力と無効電力はいずれも電圧変動に比例して変化するとは限らず、負荷の特性に応じて電圧変動に従って有効電力と無効電力は変動する。そのため電圧変動が生じた場合、上記２次方程式の係数を高い精度で求められないという問題があった。
【０００７】
また、Ｃ項とともに系統負荷の特性そのものを知ることが、系統解析を行う上で重要であるが、従来の方法ではＣ項のみしか得られないという問題があった。この発明の目的は、電力系統の特性を表す上記２次方程式の係数を正確に求められるようにした電力系統の特性推定装置および特性推定方法を提供することにある。
【０００８】
また、この発明の他の目的は、特に無効電力の系統負荷の特性定数を正確に求めることのできる電力系統の特性推定方法を提供することにある。
【０００９】
【課題を解決するための手段】
この発明は、負荷力率の変動による有効電力と無効電力の時間経過に伴う変動が生じても、Ｐ−Ｑ平面上のプロット位置のばらつきが生じないようにするため、請求項１および請求項３に記載の通り、電力系統の有効電力と無効電力を一定周期で繰り返し分析して、これらの時系列データを求めるとともに、これらの時系列データを有効電力と無効電力についてそれぞれ平滑化し、この平滑化した有効電力と無効電力の複数の対から、これらの関係を示す方程式の係数を電力系統における特性値として求める。
【００１０】
このように電力系統の有効電力と無効電力を一定周期で繰り返し求めるとともに、これをそれぞれ平滑化することによって、Ｐ−Ｑ平面上のプロット位置のばらつきが抑えられる。すなわちこれらの平滑化した有効電力と無効電力の複数の対から、それらの関係を示す方程式の係数を例えば最小二乗法により求めた場合に、それらの係数は電力系統の現実の特性値をより正確に表すものとして扱えるようになる。ここで、移動平均による有効電力と無効電力の推移の例、およびこれをＰ−Ｑ平面上へプロットした例を図４および図５に示す。サンプルデータは図２４および図２５を作成したものと同一であり、図２４および図２５に示した例と比較すれば明らかなように、実際の急激な負荷力率の変動による影響を受けずに、有効電力と無効電力の移動平均値は時間経過にともない緩やかに推移し、Ｐ−Ｑ平面上の各点のばらつきは少なくなって、各々の点が線上に並ぶ。（図５に示したとおり、この例では２４時間の間に調相設備が投入／遮断されているため、近似曲線が２本引ける。）なお、上記方程式の典型例は前述のとおり、Ｑ＝ＡＰ^２＋ＢＰ−Ｃとして２次方程式で表現されるが、一般に係数Ａは小さいので、Ｑ＝ＢＰ−Ｃとして１次方程式で表現することもできる。
【００１１】
また、この発明は電圧変動が生じても、時系列データとして求めた有効電力と無効電力の複数の対からそれらの関係を示す方程式の係数を高精度に求めるため、請求項２および請求項４に記載の通り、電力系統の電圧、有効電力および無効電力を一定周期で繰り返し求めるとともに、測定による有効電力，無効電力，電圧をそれぞれＰ，Ｑ，Ｖとし、定格電圧をＶ′とし、有効電力の負荷特性定数、無効電力の負荷特性定数をそれぞれｎｐ，ｎｑとしたとき、Ｐ′＝Ｐ（Ｖ′／Ｖ）^ｎｐで示される補正後の有効電力Ｐ′、およびＱ′＝Ｑ（Ｖ′／Ｖ）^ｎｑで示される補正後の無効電力Ｑ′をそれぞれ求め、この補正後の有効電力と無効電力の複数の対から、これらの関係を示す方程式の係数を電力系統における特性値として求める。
【００１２】
図７は上記補正によるＰ−Ｑカーブの変化の例を示す図であり、×印は補正前、□印は補正後のプロット位置である。このように電圧変動が生じても、定格電圧の下で求めた場合と等価な有効電力と無効電力を求め、この補正後の有効電力と無効電力の複数の対からそれらの関係を示す方程式の係数を求めることによって、電圧変動による影響を受けずに電力系統の特性を推定できるようになる。
【００１３】
また、この発明は上述した負荷力率の変動による問題と電圧変動による問題の双方を解消するために、請求項５に記載の通り、電力系統の電圧、有効電力および無効電力を一定周期で繰り返し求めるとともに、測定による有効電力，無効電力，電圧をそれぞれＰ，Ｑ，Ｖとし、定格電圧をＶ′とし、有効電力の負荷特性定数、無効電力の負荷特性定数をそれぞれｎｐ，ｎｑとしたとき、Ｐ′＝Ｐ（Ｖ′／Ｖ）^ｎｐで示される補正後の有効電力Ｐ′およびＱ′＝Ｑ（Ｖ′／Ｖ）^ｎｑで示される補正後の無効電力Ｑ′をそれぞれ求め、この補正後の有効電力と無効電力の時系列データをそれぞれ平滑化し、この平滑化した有効電力と無効電力の複数の対から、これらの関係を示す方程式の係数を電力系統における特性値として求める。
【００１４】
この方法によって負荷力率の変動および電圧変動が比較的激しい場合であっても、電力系統の特性をより正確に推定できるようになる。
【００１５】
また、この発明は上記有効電力および無効電力の負荷特性定数を自動的に求めるようにするため、請求項６に記載の通り、電力系統の有効電力と無効電力を一定周期で分析してこれらの時系列データを求め、電力系統の電圧の擾乱を検出するとともに、その擾乱前の有効電力、無効電力および電圧をＰｏ，Ｑｏ，Ｖｏとし、擾乱後の有効電力、無効電力および電圧をＰ１，Ｑ１，Ｖ１としたとき、Ｐ１＝Ｐｏ（Ｖ１／Ｖｏ）^ｎｐで示されるｎｐ、およびＱ１＝Ｑｏ（Ｖ１／Ｖｏ）^ｎｑで示されるｎｑをそれぞれ求め、請求項４または５に記載の有効電力の負荷特性定数ｎｐおよび無効電力の負荷特性定数ｎｑにそれぞれ適用する。
【００１６】
このように電圧のステップ変化などの擾乱があったときに、その擾乱前後の有効電力、無効電力および電圧から有効電力および無効電力の負荷特性定数ｎｐ，ｎｑが自動的に求められて、これを用いて電圧変動に応じた有効電力と無効電力の補正が行われる。
【００１７】
さらに、この発明は、系統負荷の無効電力負荷特性定数を正確に求めるようにするために、請求項７に記載の通り、電力系統の有効電力、無効電力の時系列データと、電圧ステップ変化データとから下記のステップ（１）〜（３）により、無効電力の負荷特性定数ｎｑを推定することを特徴とする。
【００１８】
（１）Ｐ′＝Ｐ（Ｖ′／Ｖ）^ｎｐ
Ｑ′＝｛Ｑ＋Ｃ′（Ｖ／Ｖ′）^２｝（Ｖ′／Ｖ）^ｎｑ−Ｃ′
より（Ｐ，Ｑ，Ｖ）を定格電圧（Ｖ′）時に補正してＣ′を推定する。但し、最初はＣ′，ｎｑを任意の値とし、ｎｐは、
Ｐ_１＝Ｐ_０（Ｖ_１／Ｖ_０）^ｎｐより求める。
【００１９】
（２）上記（１）で推定したＣ′と電圧ステップ変化データとにより、次式からｎｑを推定する。
【００２０】
Ｑ_１＝（Ｑ_０＋Ｃ′）（Ｖ_１／Ｖ_０）^ｎｑ−Ｃ′（Ｖ_１／Ｖ_０）^２
（３）上記（１）で推定したＣ′と（２）で推定したｎｑを用いて（１）以下を所定回数繰り返してｎｑを求める。
【００２１】
但し、
Ｐ，Ｑ，Ｖ：測定（実）データ（有効電力、無効電力、電圧）
Ｐ′，Ｑ′，Ｖ′，Ｃ′：定格電圧時補正データ（有効電力、無効電力、電圧、負荷系統容量値（Ｃ項））
Ｐ_０，Ｑ_０，Ｖ_０：電圧ステップ変化前データ（有効電力、無効電力、電圧）
Ｐ_１，Ｑ_１，Ｖ_１：電圧ステップ変化後データ（有効電力、無効電力、電圧）
ｎｐ：有効電力の負荷特性定数
ｎｑ：無効電力の負荷特性定数
いま、バンク二次側の６６ｋＶ系以下を縮約した系統のモデルを示すと図１０に示すようになる。Ｐ，Ｑ，Ｖはそれぞれ有効電力、無効電力、電圧であり、Ｃ項は、このモデルにおいて−ｊＱｃとして示している。
【００２２】
図１０に示すモデルにおいて、系統の電圧−有効電力特性、電圧−無効電力特性は、電圧ＶがＶ_０からＶ_１まで擾乱によってステップ状に変化した場合、次式で表される。
【００２３】
【数１】

【００２４】
上記式において、負荷の電圧特性は、それぞれ、定電力特性、定電流特性、定インピーダンス特性に分けることができる。
【００２５】
（ａ）定電力特性・・・電圧が変化しても、電力が一定な負荷特性。（ｎｐ、ｎｑ＝０）
（ｂ）定電流特性・・・電圧が変化したとき、電力が電圧に比例する負荷特性。（ｎｐ、ｎｑ＝１）
（ｃ）定インピーダンス特性・・・電圧が変化したとき、電力が電圧の２乗に比例する負荷特性。（ｎｐ、ｎｑ＝２）
ここで、−ｊＱ_ｃは定インピーダンス特性であるから、Ｃ項についてはｎｑ＝２となる。従って、Ｖ_０からＶ_１へのステップ電圧変化時には、Ｃ（Ｖ_１／Ｖ_０）^２のＱ（無効電力）が−ｊＱ_ｃにおいて発生する。
【００２６】
よって、ステップ電圧変化時におけるＣ項の変動分を考慮した場合の負荷Ｌにおける上記（２）式は次のようになる。ただし、ここではバンク分については無視する。
【００２７】
【数２】

【００２８】
同式において、右項の（Ｑ_０＋Ｃ）はステップ電圧変化前のＱ、左項はステップ電圧変化後のＱを表す。上記（３）式より、無効電力の負荷特性定数ｎｑは、次の（４）式で求めることができる。
【００２９】
【数３】

【００３０】
なお、系統の電圧−有効電力特性を表す上記（１）式についてはＣ項を考慮する必要がない。従って、有効電力の負荷特性定数ｎｐは、上記（１）式においてそれぞれ電圧ステップ変化データを代入することによりｎｐを求めることができる。
【００３１】
ところで、上記（３）式は、母線電圧（６６ｋＶ）の変動を考慮していないために、ｎｑの値は正確ではない。そこで、請求項４、請求項５に示すように、母線電圧を定格電圧に補正するのが望ましい。この場合、Ｃ項とｎｑは相互に関連し且つ未知数である。今、測定データ（実データ）を、Ｐ，Ｑ，Ｖ、定格電圧補正後のデータをＰ′，Ｑ′，Ｖ′とし、さらに定格電圧補正後のＣ項をＣ′とすると、この時のＱ′は上記（３）式に基づいて（５）式のように表される。
【００３２】
【数４】

【００３３】
上記（５）式の右項のＣ′（Ｖ／Ｖ′）^２は、Ｃ′が定格電圧補正後の値であるために、これを実データに対応する値にするためである。すなわち、Ｃ′に（Ｖ／Ｖ′）^２を乗ずることにより実データに対応するＣ項を求めるようにしている。なお、有効電力分についてはＣ項は無関係であるために、上記電圧ステップ変化データを用いて（６）式より有効電力の負荷特性定数ｎｐを求め、このｎｐを用いて（７）式よりＰ′を求める。
【００３４】
ところで、上記（５）式においてはＣ′およびｎｑが未知数かつ相互に関連する値である。そこで、最初にこのＣ′およびｎｑにそれぞれ適当な値を代入して定格電圧補正後のＰ−Ｑ相関データＰ′，Ｑ′，Ｖ′を、実データＰ，Ｑ，Ｖより求める。Ｃ′，ｎｑの適当な値としては、一例として、Ｃ′＝０，ｎｑ＝０，１０とする。
【００３５】
上記のようにして定格電圧補正後のＰ−Ｑ相関データＰ′，Ｑ′，Ｖ′により図１１に示すようにＰ′−Ｑ′相関カーブからＣ項（Ｃ′）を推定する。この段階で、Ｃ項の値が最初に設定した任意の値よりもより正確な値に近づく。続いて、推定したＣ項（Ｃ′）と電圧ステップ変化データとによりｎｑを推定する。このｎｑ推定は、上記（３）式を用いる。すなわち、（３）式のＣを上記のようにして推定したＣ′に置き換えることにより、上記（４）式によってｎｑを推定する。推定されたｎｑは、最初に設定した任意の値よりもより真の値に近づいている。
【００３６】
続いて、このｎｑとＣ′を用いて上記（５）および（７）式により定格電圧補正後のＰ′−Ｑ′相関データを求め、再びＣ項（Ｃ′）を推定する。また、上記（３）式（４）式により再びｎｑを推定する。この動作を繰り返していくことにより、Ｃ項（Ｃ′）とｎｑがより真の値に近づき収束していく。ｎｑが適当な範囲まで集束した段階でｎｑおよびＣ項を確定する。なお、以上の繰り返し演算において、未知数のＣ′，ｎｑを最初に適当な数に設定するが、この値が真の値から大きく離れていたとしても、Ｃ′やｎｑが発散しないことがシミュレーションにより確かめられている。
【００３７】
以上の方法により、Ｃ項とｎｑを高精度に求めることができる。
【００３８】
【発明の実施の形態】
この発明の実施形態である負荷特性記録装置の構成を図１〜図９を基に以下説明する。
【００３９】
図１は装置全体の構成を示すブロック図である。図１において１００は現場に設けられた端末装置、１０１は制御室に設けられた監視装置である。絶縁・変換回路１は３相電圧・電流信号を絶縁入力し、所定の変成比で電圧信号に変換する。フィルタ２はサンプリングによる折り返し誤差による影響を受けないために、所定の周波数以上の成分を除去する。サンプルホールド回路３は入力信号を所定タイミングでサンプリングしホールドする。マルチプレクサ４は複数のサンプルホールド回路３のうち１つを選択してＡＤコンバータ５へ与える。ＡＤコンバータ５はこれをデジタルデータに変換する。マイクロコンピュータ６はマルチプレクサ４を選択するとともに、ＡＤコンバータ５により変換されたデータを読み取る。通信インタフェース７は端末装置１００と監視装置１０１との間でデータ伝送を行うために用い、マイクロコンピュータ６は通信インタフェース７を介して監視装置１０１へ各種計測データなどを伝送する。一方、監視装置１０１において通信インタフェース８は端末装置１００との間でデータ伝送制御を行い、マイクロコンピュータ９は通信インタフェース８を介して端末装置１００から各種計測データなどを受信する。表示器１１およびプロッタ１３は計測値などを出力するために用い、マイクロコンピュータ９は表示インタフェース１０およびプロッタインタフェース１２を介してデータの出力を行う。
【００４０】
図２はサンプリングおよび特性値算出の処理手順を示すフローチャートである。まず、あらかじめ定めたサンプリングタイミングとなれば、（例えば５分毎に、）３相電圧・電流信号のサンプリングを行い、ＡＤ変換値を読み取る。そして有効電力Ｐおよび無効電力Ｑを算出し、更に次式で有効電力と無効電力の負荷特性定数ｎｐ，ｎｑによる補正（以下「定格電圧補正」という。）後の有効電力Ｐ′と無効電力Ｑ′をそれぞれ求める。
【００４１】
Ｐ′＝Ｐ（Ｖ′／Ｖ）^ｎｐ
Ｑ′＝Ｑ（Ｖ′／Ｖ）^ｎｑ
ここで、Ｖは測定電圧、Ｖ′は定格電圧である。
【００４２】
その後、有効電力Ｐ′および無効電力Ｑ′の移動平均値を求める。すなわちｋ回目に求めた有効電力をＰ（ｋ）’、無効電力をＱ（ｋ）’とすれば、時刻ｔにおける移動平均値Ｐｔ’，Ｑｔ’ はそれぞれ次式で示される。
【００４３】

【００４４】
上式でｍは平均化の幅であり、図４に示した例ではｍ＝９としている。
【００４５】
このようにして求めた有効電力と無効電力の移動平均値をＰ−Ｑ平面上にプロットする。これにより、図５に示したように、各点のばらつきは少なくなって、各々の点が線上に並ぶ。
【００４６】
以上の処理をサンプリングタイミングになる毎に繰り返す。そして、特性値算出タイミングとなれば、例えば２４時間が経過する毎（５分毎にサンプリングを行った場合、２８８点分のサンプリングを行う毎）に、有効電力と無効電力の移動平均値の対を基に、最小二乗法によりＱ＝ＡＰ^２＋ＢＰ−Ｃの曲線のあてはめを行い、各係数Ａ，Ｂ，Ｃをそれぞれ算出する。そして、上記曲線をＰ−Ｑカーブとして描画するとともに、係数Ａ，Ｂ，Ｃの値をそれぞれ出力する。
【００４７】
図３は有効電力および無効電力の負荷特性定数算出手順を示すフローチャートである。同図に示すように、もし測定電圧にステップ変化（擾乱）が生じたことを判定すれば、その擾乱前後の有効電力、無効電力および電圧の値から有効電力の負荷特性定数ｎｐおよび無効電力の負荷特性定数ｎｑをそれぞれ算出する。図６に示すように、例えば５分毎に求めた電圧値が前回の電圧値に比較して一定値以上の差が発生すれば、擾乱があったものと見なして、各擾乱が生じる毎にｎｐ，ｎｑを算出する。すなわちステップ変化前の有効電力、無効電力および電圧をＰｏ，Ｑｏ，Ｖｏとし、ステップ変化後の有効電力、無効電力および電圧をＰ１，Ｑ１，Ｖ１としたとき、次式を満足するｎｐ，ｎｑを求める。
【００４８】
Ｐ１＝Ｐｏ（Ｖ１／Ｖｏ）^ｎｐ
Ｑ１＝Ｑｏ（Ｖ１／Ｖｏ）^ｎｑ
例えば
ｎｐ＝ｌｏｇ（Ｐ１／Ｐｏ）／ｌｏｇ（Ｖ１／Ｖｏ）
ｎｑ＝ｌｏｇ（Ｑ１／Ｑｏ）／ｌｏｇ（Ｖ１／Ｖｏ）
を算出することによってｎｐ，ｎｑを求める。その後、一定時間（例えば２４時間）経過するまで上記処理を繰り返し、一定時間が経過すれば昼間と夜間とでそれぞれｎｐ，ｎｑの平均値を算出する。図２に示した有効電力と無効電力の補正ステップでは、昼間と夜間とでそれぞれのｎｐ，ｎｑを用いて補正を行う。
【００４９】
尚、上述の例では、時系列での移動平均処理と定格電圧補正の双方を行う例を示したが、これらの処理は互いに独立であるため、いずれか一方のみの処理を行っても、電力系統の特性を表す上記２次方程式の係数の推定精度を高めることができる。例えば、住宅負荷系統のように負荷力率変動の少ない系統について、特性推定を行う場合は、移動平均処理を行わずに、定格電圧補正のみを行っても、その効果は大きい。その例を図８および図９に示す。すなわち、図８および図９は住宅地に対する定格電圧６６ｋＶのフィーダについて計測した際の、Ｐ−Ｑ平面上のプロット例およびＰ−Ｑカーブであり、図８は定格電圧補正を行わない場合、図９は行った場合である。ここでは負荷を有効電力に対しては定電流特性、すなわちｎｐ＝１とし、無効電力については定インピーダンス特性、すなわちｎｑ＝２として補正を行った。両図に示すように、この例では定格電圧補正を行わないで求めたＣ項は＋７．４９２Ｍｖａｒであるのに対し、補正を行った求めたＣ項は＋９．４５４Ｍｖａｒとなって、大きく変化することがわかる。換言すれば、それだけ特性推定精度が高まることがわかる。
【００５０】
次に、実Ｐ−Ｑ相関データＰ，Ｑ，Ｖを定格電圧補正し、精度の高いＣ項と無効電力負荷特性定数ｎｑを求める実施形態について説明する。
【００５１】
装置の構成については図１に示すものと同一である。図１２は、Ｃ項およびｎｑを推定するための処理手順を示すフローチャートである。
【００５２】
有効電力Ｐおよび無効電力Ｑを算出するまでは図２に示す動作と同じであり、何回か繰り返すことによって、時系列の、実Ｐ−Ｑ相関データＰ，Ｑ，Ｖが算出される。そして、ｎｑ，Ｃ項を推定値出力する適当なタイミングになれば、以下の手順を繰り返す。
【００５３】
（１）実Ｐ−Ｑ相関データより、定格電圧補正後のＰ′−Ｑ′相関データを求める（（５）および（７）式参照）。但し、最初は、同式におけるＣ項（Ｃ′）およびｎｑが未知数であるために、適当な値（ここではＣ′＝０、ｎｑ＝０，１０）に設定する。
【００５４】
（２）Ｐ′−Ｑ′相関カーブを描画する。
【００５５】
（３）Ｃ項（Ｃ′）を求め、推定値とする。
【００５６】
（４）電圧ステップ変化データ（Ｖ_０→Ｖ_１）を用いて、ｎｑを求め推定値とする（（３）式参照）。なお電圧ステップ変化データは、測定（実）電圧にステップ変化（擾乱）が生じたときの変化データＶ_０，Ｖ_１を用いる。
【００５７】
（５）ｎｑが所定のレベルに収束するまで（ここでは、０．９５＜（ｎｑ_ｉ／ｎｑ_ｉ＋１）＜１．０５になるまで）、上記（１）〜（４）を繰り返す。
【００５８】
（６）一定レベルまで収束した段階で、そのときのｎｑおよびＣ項（Ｃ′）を出力する。
【００５９】
なお、上記Ｐ′−Ｑ′相関カーブを描画してＣ項（Ｃ′）を推定するときに、図２で示したような移動平均処理を行うことも可能である。
【００６０】
図１３〜図２０にシミュレーション結果を示す。このシミュレーションでは真のｎｐ，ｎｑをそれぞれ２、４としてバンク２次側電圧の変動がない場合とある場合、およびバンクＺ_ｓのない場合とある場合、並びに線路インピーダンスＺ_Ｌがない場合とある場合のそれぞれのシミュレーション結果を示すものである。なお初回のｎｑを０，１０としている。これらのシミュレーション結果に示すように、ｎｑの推定演算を少なくとも５回行うまでにｎｑが収束することがわかる。すなわち、電圧にばらつきのあるＰ−Ｑ相関データにおいて、ｎｑを繰り返し演算することにより、電圧補正しない場合に比べてＣ項とｎｑを精度よく求めることができる。
【００６１】
【発明の効果】
請求項１および請求項３に記載の発明によれば、電力系統の有効電力と無効電力を一定周期で繰り返し求めるとともに、これをそれぞれ平滑化したため、負荷力率の変動が激しくても、これらの平滑化した有効電力と無効電力の複数の対から、それらの関係を示す方程式の係数を推定した場合に、それらの係数は電力系統の現実の特性値をより正確に表すものとして扱えるようになる。
【００６２】
請求項２および請求項４に記載の発明によれば、常に定格電圧の下で求めた場合と等価な有効電力と無効電力を求め、この補正後の有効電力と無効電力の複数の対からそれらの関係を示す方程式の係数を求めるようにしたため、電圧変動が生じても電力系統の特性推定の精度が高まる。
【００６３】
請求項５に記載の発明によれば、負荷力率の変動および電圧変動の双方が激しい場合であっても、電力系統の特性をより正確に推定できるようになる。
【００６４】
請求項６に記載の発明によれば、電圧のステップ変化などの擾乱があったときに、その擾乱前後の有効電力、無効電力および電圧から有効電力および無効電力の負荷特性定数ｎｐ，ｎｑが自動的に求められるため、有効電力および無効電力の負荷特性定数を求めるための特別な計測が不要となる。
【００６５】
請求項７に記載の発明によれば、電力系統に電圧変動があった場合でも、ｎｑを繰り返し演算することにより高精度のｎｑおよびＣ項を求めることができ、系統負荷の特性を正確に把握することができる。
【図面の簡単な説明】
【図１】負荷特性記録装置の構成を示すブロック図である。
【図２】サンプリングおよび特性値算出の処理手順を示すフローチャートである。
【図３】有効電力・無効電力の負荷特性定数算出手順を示すフローチャートである。
【図４】有効電力と無効電力の時間経過にともなう推移の例を示す図である。
【図５】Ｐ−Ｑ平面上のプロットの例を示す図である。
【図６】負荷特性定数の算出タイミングの例を示す図である。
【図７】負荷特性定数による定格電圧補正を行う前後のＰ−Ｑカーブの変化の例を示す図である。
【図８】負荷特性定数による定格電圧補正前のＰ−Ｑカーブの例を示す図である。
【図９】負荷特性定数による定格電圧補正後のＰ−Ｑカーブの例を示す図である。
【図１０】バンク２次側以下を６６ｋＶに縮約した場合の電力系統のモデル図である。
【図１１】Ｐ−Ｑ相関データを定格電圧補正し、Ｃ′を求める方法を説明する図である。
【図１２】ｎｑおよびＣ′を繰り返し演算により収束させて推定する手順を示すフローチャートである。
【図１３】〜
【図２０】シミュレーション結果を示す図である。
【図２１】負荷系統の等価回路図である。
【図２２】Ｐ−Ｑカーブの例を示す図である。
【図２３】変動負荷系統の等価回路図である。
【図２４】有効電力と無効電力の時間経過にともなう推移の例を示す図である。
【図２５】Ｐ−Ｑ平面上のプロットの例を示す図である。
【符号の説明】
１００−端末装置
１０１−監視装置[0001]
BACKGROUND OF THE INVENTION
The present invention relates to an apparatus and method for estimating power system characteristics, and more particularly, to a power system characteristic estimation apparatus and characteristic estimation method that facilitate estimation of phase-adjusted equipment capacity and system load characteristic constants of reactive power.
[0002]
[Prior art]
Conventionally, a phase adjustment facility has been provided in the load system for power factor improvement and phase adjustment. However, it is important to know how much capacity is available and how the load characteristics of the system are. This is one of the important points in system planning.
[0003]
In general, an actual load system is simulated by contracting below the secondary bus of the primary substation, and the relationship between active power and reactive power is related by a certain equation. If a typical equivalent circuit of an actual load system is represented as shown in FIG. 21, the relationship between the active power P and the reactive power Q is expressed by Q = AP ² + BP-C. Here, A is a value related to the reactance of the line of the power system, AP ² is a reactive power due to the reactance of the line, B is a value related to the power factor of the load, BP is a reactive power due to the load, and C is a power. It corresponds to the reactive power due to the charging capacity of the grid and the capacity of the phase adjusting equipment. Therefore, if a pair of active power and reactive power measured repeatedly at a constant period is plotted on a rectangular coordinate plane (hereinafter referred to as “PQ plane”) of active power and reactive power, an ideal power is obtained as shown in FIG. Specifically, they are arranged on a curve represented by the above quadratic equation (hereinafter referred to as “PQ curve”).
[0004]
Therefore, conventionally, coefficients A, B, and C of the above quadratic equation are obtained by a least square method from a pair of active power and reactive power obtained by using a load characteristic recording device that analyzes active power and reactive power of a three-phase power system. Are seeking each. In particular, the capacity of the phase adjusting equipment is estimated on the basis of a coefficient C (hereinafter referred to as “C term”), and this is used as information for stabilizing the voltage, for example, for proper operation of the power system. Yes.
[0005]
[Problems to be solved by the invention]
However, in the actual load system, the load power factor fluctuates as shown in the equivalent circuit of FIG. 23, so that the value of the BP term fluctuates and it may be difficult to estimate an important C term. Especially, the load power factor fluctuates greatly in the factory load system. Here, if the transition of the active power and the reactive power in the factory load system is specifically expressed, it becomes as shown in FIG. 24, and it fluctuates greatly with the passage of time. Therefore, even if a pair of active power and reactive power is plotted on the PQ plane, it varies greatly as shown in FIG. (In the example shown in FIG. 25, since two approximate curves can be drawn, it is expected that the phase adjusting equipment is turned on / off during 24 hours.) Therefore, based on such data Even if the coefficient of the quadratic equation is obtained by the least square method, only a value with low reliability can be obtained.
[0006]
Conventionally, active power and reactive power data obtained at regular intervals are used as they are, but in general, both active power and reactive power do not always change in proportion to voltage fluctuations, depending on the characteristics of the load. The active power and reactive power fluctuate according to the voltage fluctuation. Therefore, when voltage fluctuation occurs, there is a problem that the coefficient of the quadratic equation cannot be obtained with high accuracy.
[0007]
Further, knowing the characteristics of the system load itself together with the C term is important for system analysis, but there is a problem that only the C term can be obtained by the conventional method. An object of the present invention is to provide a power system characteristic estimation device and a characteristic estimation method in which the coefficient of the quadratic equation representing the characteristics of the power system can be accurately obtained.
[0008]
Another object of the present invention is to provide a power system characteristic estimation method capable of accurately obtaining a characteristic constant of a reactive power system load.
[0009]
[Means for Solving the Problems]
In order to prevent variations in plot positions on the PQ plane from occurring even when fluctuations occur in the active power and reactive power with the passage of time due to fluctuations in the load power factor, the present invention does not cause variations in plot positions. As described in FIG. 3, the active power and reactive power of the power system are repeatedly analyzed at regular intervals to obtain these time series data, and these time series data are smoothed for the active power and reactive power respectively. The coefficient of the equation showing these relations is obtained as a characteristic value in the power system from the plurality of pairs of activated active power and reactive power.
[0010]
In this way, the active power and reactive power of the power system are repeatedly obtained at a constant period, and by smoothing each of them, variations in plot positions on the PQ plane can be suppressed. In other words, when the coefficients of the equations indicating the relationship between these smoothed active power and reactive power are obtained by, for example, the least square method, these coefficients more accurately represent the actual characteristic values of the power system. Can be handled as Here, FIG. 4 and FIG. 5 show an example of transition of active power and reactive power by moving average, and an example of plotting this on the PQ plane. The sample data is the same as that created in FIG. 24 and FIG. 25, and as is clear from comparison with the examples shown in FIG. 24 and FIG. 25, it is not affected by the actual rapid change of the load power factor. The moving average values of the active power and the reactive power gradually change with time, the variation of the points on the PQ plane is reduced, and the points are arranged on a line. (As shown in FIG. 5, in this example, since the phase adjusting equipment is turned on / off for 24 hours, two approximate curves can be drawn.) Note that a typical example of the above equation is Q = Although it is expressed by a quadratic equation as AP ² + BP-C, since the coefficient A is generally small, it can also be expressed by a linear equation as Q = BP-C.
[0011]
Further, according to the present invention, in order to obtain the coefficient of the equation indicating the relationship from a plurality of pairs of active power and reactive power obtained as time series data with high accuracy even when voltage fluctuation occurs,

claims

2 and 4 As described in, the power system voltage, active power, and reactive power are repeatedly obtained at regular intervals, and the active power, reactive power, and voltage are P, Q, and V, and the rated voltage is V '. When the load characteristic constants and reactive power load characteristic constants are np and nq, respectively, P ′ = P (V ′ / V) corrected active power P ′ indicated by ^np and Q ′ = Q (V ′ / V) Reactive reactive power Q ′ indicated by ^nq is obtained, and coefficients of equations indicating these relationships are obtained as characteristic values in the power system from the corrected active power and reactive power pairs.
[0012]
FIG. 7 is a diagram showing an example of a change in the PQ curve due to the above correction, where x marks are before correction, and □ are plot positions after correction. Even if voltage fluctuations occur in this way, the active power and reactive power equivalent to those obtained under the rated voltage are obtained, and the equation indicating the relationship between these corrected active power and reactive power pairs is obtained. By obtaining the coefficient, it becomes possible to estimate the characteristics of the power system without being affected by voltage fluctuations.
[0013]
Further, in order to solve both the problem caused by the fluctuation of the load power factor and the problem caused by the voltage fluctuation, the present invention repeats the voltage, active power and reactive power of the power system at a constant cycle as described in claim 5. When the active power, reactive power, and voltage by measurement are P, Q, and V, the rated voltage is V ′, the load characteristic constant of active power and the load characteristic constant of reactive power are np and nq, respectively, P ′ = P (V ′ / V) The corrected effective power P ′ indicated by ^np and the corrected reactive power Q ′ indicated by Q ′ = Q (V ′ / V) ^nq are respectively obtained, and this corrected The active power and reactive power time-series data are each smoothed, and coefficients of equations indicating these relationships are obtained as characteristic values in the power system from the smoothed active power and reactive power pairs.
[0014]
This method makes it possible to more accurately estimate the characteristics of the power system even when the load power factor fluctuation and voltage fluctuation are relatively severe.
[0015]
Further, in order to automatically obtain the load characteristic constants of the active power and reactive power according to the present invention, the active power and reactive power of the power system are analyzed at regular intervals as described in claim 6 to obtain these Time series data is obtained, voltage disturbance of the power system is detected, active power, reactive power and voltage before the disturbance are set to Po, Qo and Vo, and active power, reactive power and voltage after the disturbance are P1, Q1. , V1, and P1 = Po (V1 / Vo) ^np indicated by np and Q1 = Qo (V1 / Vo) ^nq indicated by nq, respectively, and the load of active power according to

claim

4 or 5 This is applied to the characteristic constant np and the reactive power load characteristic constant nq, respectively.
[0016]
When there is a disturbance such as a step change in voltage in this way, the load characteristic constants np and nq of the active power and reactive power are automatically obtained from the active power, reactive power and voltage before and after the disturbance. The active power and the reactive power are corrected according to the voltage fluctuation.
[0017]
Furthermore, in order to accurately obtain the reactive power load characteristic constant of the system load, the present invention provides the power system active power, reactive power time series data, and voltage step change data as described in claim 7. The reactive power load characteristic constant nq is estimated by the following steps (1) to (3).
[0018]
(1) P ′ = P (V ′ / V) ^np
Q ′ = {Q + C ′ (V / V ′) ² } (V ′ / V) ^nq −C ′
Then, (P, Q, V) is corrected at the rated voltage (V ') to estimate C'. However, initially C ′ and nq are arbitrary values, and np is
P ₁ = P ₀ (V ₁ / V ₀ ) ^Obtained from ^np .
[0019]
(2) nq is estimated from the following equation using C ′ and voltage step change data estimated in (1) above.
[0020]
Q ₁ = (Q ₀ + C ′) (V ₁ / V ₀ ) ^nq −C ′ (V ₁ / V ₀ ) ²
(3) Using the C ′ estimated in the above (1) and the nq estimated in (2), (1) is repeated a predetermined number of times to obtain nq.
[0021]
However,
P, Q, V: Measurement (actual) data (active power, reactive power, voltage)
P ′, Q ′, V ′, C ′: Correction data at rated voltage (active power, reactive power, voltage, load system capacity value (C term))
P ₀ , Q ₀ , V ₀ : Data before voltage step change (active power, reactive power, voltage)
P ₁ , Q ₁ , V ₁ : Data after voltage step change (active power, reactive power, voltage)
np: Load characteristic constant of active power nq: Load characteristic constant of reactive power Now, FIG. 10 shows a model of a system in which the 66 kV system or less on the bank secondary side is reduced. P, Q, and V are active power, reactive power, and voltage, respectively, and C term is shown as -jQc in this model.
[0022]
In the model shown in FIG. 10, the voltage-active power characteristic and the voltage-reactive power characteristic of the system are expressed by the following equations when the voltage V changes stepwise due to disturbance from V ₀ to V ₁ .
[0023]
[Expression 1]

[0024]
In the above equation, the voltage characteristics of the load can be divided into constant power characteristics, constant current characteristics, and constant impedance characteristics, respectively.
[0025]
(A) Constant power characteristics: Load characteristics in which power is constant even when the voltage changes. (Np, nq = 0)
(B) Constant current characteristics: Load characteristics in which power is proportional to voltage when the voltage changes. (Np, nq = 1)
(C) Constant impedance characteristics: Load characteristics in which the power is proportional to the square of the voltage when the voltage changes. (Np, nq = 2)
Here, since -jQ _c is a constant impedance characteristic, nq = 2 for the C term. Therefore, when the step voltage changes from V ₀ to V ₁ , Q (reactive power) of C (V ₁ / V ₀ ) ² is generated at −jQ _c .
[0026]
Therefore, the above equation (2) for the load L when considering the variation of the C term when the step voltage changes is as follows. However, the bank is ignored here.
[0027]
[Expression 2]

[0028]
In the equation, (Q ₀ + C) in the right term represents Q before the step voltage change, and the left term represents Q after the step voltage change. From the above equation (3), the load characteristic constant nq of reactive power can be obtained by the following equation (4).
[0029]
[Equation 3]

[0030]
In addition, it is not necessary to consider C term about the said (1) Formula showing the voltage-active power characteristic of a system | strain. Therefore, the load characteristic constant np of the active power can be obtained by substituting the voltage step change data in the above equation (1).
[0031]
By the way, since the above equation (3) does not consider the fluctuation of the bus voltage (66 kV), the value of nq is not accurate. Therefore, as shown in

claims

4 and 5, it is desirable to correct the bus voltage to the rated voltage. In this case, the C term and nq are interrelated and unknown. Assuming that the measured data (actual data) is P, Q, V, the rated voltage corrected data is P ′, Q ′, V ′, and the C voltage after the rated voltage correction is C ′, Q ′ is expressed as in equation (5) based on equation (3) above.
[0032]
[Expression 4]

[0033]
C ′ (V / V ′) ^{2 in} the right term of the above expression (5) is for making C ′ a value corresponding to actual data because C ′ is a value after the rated voltage correction. That is, the C term corresponding to the actual data is obtained by multiplying C ′ by (V / V ′) ² . Since the C term is not related to the active power, the load characteristic constant np of the active power is obtained from the equation (6) using the voltage step change data, and P is obtained from the equation (7) using this np. Find ′.
[0034]
By the way, in the above equation (5), C ′ and nq are unknown values and mutually related values. Therefore, first, appropriate values are substituted into C ′ and nq, respectively, and the PQ correlation data P ′, Q ′, V ′ after the rated voltage correction is obtained from the actual data P, Q, V. Appropriate values for C ′ and nq are, for example, C ′ = 0 and nq = 0,10.
[0035]
The C term (C ′) is estimated from the P′-Q ′ correlation curve as shown in FIG. 11 from the PQ correlation data P ′, Q ′, V ′ after the rated voltage correction as described above. At this stage, the value of the C term approaches a more accurate value than the initially set arbitrary value. Subsequently, nq is estimated from the estimated C term (C ′) and voltage step change data. This nq estimation uses the above equation (3). That is, nq is estimated by the above equation (4) by replacing C in the equation (3) with C ′ estimated as described above. The estimated nq is closer to a true value than the initially set arbitrary value.
[0036]
Subsequently, using this nq and C ′, P′-Q ′ correlation data after rated voltage correction is obtained by the above equations (5) and (7), and the C term (C ′) is estimated again. Further, nq is estimated again by the above formula (3) and formula (4). By repeating this operation, the C term (C ′) and nq become closer to the true value and converge. When nq converges to an appropriate range, nq and C terms are determined. In the above iterative calculation, the unknown numbers C ′ and nq are set to appropriate numbers at first. However, even if this value is far from the true value, C ′ and nq do not diverge by simulation. It has been confirmed.
[0037]
By the above method, the C term and nq can be obtained with high accuracy.
[0038]
DETAILED DESCRIPTION OF THE INVENTION
The configuration of the load characteristic recording apparatus according to the embodiment of the present invention will be described below with reference to FIGS.
[0039]
FIG. 1 is a block diagram showing the overall configuration of the apparatus. In FIG. 1, 100 is a terminal device provided in the field, and 101 is a monitoring device provided in the control room. The insulation / conversion circuit 1 receives a three-phase voltage / current signal by insulation and converts it into a voltage signal at a predetermined transformation ratio. Since the filter 2 is not affected by the aliasing error due to sampling, a component having a predetermined frequency or higher is removed. The sample hold circuit 3 samples and holds the input signal at a predetermined timing. The multiplexer 4 selects one of the plurality of sample and hold circuits 3 and supplies it to the AD converter 5. The AD converter 5 converts this into digital data. The microcomputer 6 selects the multiplexer 4 and reads the data converted by the AD converter 5. The communication interface 7 is used for data transmission between the terminal device 100 and the monitoring device 101, and the microcomputer 6 transmits various measurement data to the monitoring device 101 via the communication interface 7. On the other hand, in the monitoring device 101, the communication interface 8 performs data transmission control with the terminal device 100, and the microcomputer 9 receives various measurement data and the like from the terminal device 100 via the communication interface 8. The display 11 and the plotter 13 are used to output measurement values and the like, and the microcomputer 9 outputs data via the display interface 10 and the plotter interface 12.
[0040]
FIG. 2 is a flowchart showing processing procedures for sampling and characteristic value calculation. First, when a predetermined sampling timing is reached, sampling of a three-phase voltage / current signal is performed (for example, every 5 minutes), and an AD conversion value is read. Then, the active power P and the reactive power Q are calculated, and the active power P ′ and the reactive power Q after correction (hereinafter referred to as “rated voltage correction”) of the active power and the reactive power by the load characteristic constants np and nq are calculated by the following equations. ′ Is obtained.
[0041]
P ′ = P (V ′ / V) ^np
Q ′ = Q (V ′ / V) ^nq
Here, V is a measurement voltage and V ′ is a rated voltage.
[0042]
Thereafter, moving average values of the active power P ′ and the reactive power Q ′ are obtained. That is, if the active power obtained at the k-th time is P (k) ′ and the reactive power is Q (k) ′, the moving average values Pt ′ and Qt ′ at time t are represented by the following equations, respectively.
[0043]

[0044]
In the above equation, m is the width of averaging, and m = 9 in the example shown in FIG.
[0045]
The moving average values of the active power and the reactive power thus obtained are plotted on the PQ plane. Thereby, as shown in FIG. 5, the variation of each point decreases, and each point is arranged on a line.
[0046]
The above processing is repeated every time the sampling timing is reached. Then, at the characteristic value calculation timing, for example, every 24 hours (when sampling is performed every 5 minutes, every time 288 points are sampled), a pair of active power and reactive power moving average values Based on the above, a curve of Q = AP ² + BP-C is applied by the least square method, and each coefficient A, B, C is calculated. And while drawing the said curve as a PQ curve, the value of coefficient A, B, C is each output.
[0047]
FIG. 3 is a flowchart showing a procedure for calculating load characteristics constants of active power and reactive power. As shown in the figure, if it is determined that a step change (disturbance) has occurred in the measured voltage, the load characteristic constant np and reactive power of the active power are calculated from the values of active power, reactive power and voltage before and after the disturbance. Load characteristic constants nq are respectively calculated. As shown in FIG. 6, for example, if a voltage value obtained every 5 minutes differs from a previous voltage value by a certain value or more, it is considered that there is a disturbance, and each time a disturbance occurs. np and nq are calculated. That is, when active power, reactive power, and voltage before step change are Po, Qo, and Vo, and active power, reactive power, and voltage after step change are P1, Q1, and V1, np and nq satisfying the following equations are obtained. Ask.
[0048]
P1 = Po (V1 / Vo) ^np
Q1 = Qo (V1 / Vo) ^nq
For example, np = log (P1 / Po) / log (V1 / Vo)
nq = log (Q1 / Qo) / log (V1 / Vo)
Np and nq are obtained by calculating. Thereafter, the above process is repeated until a certain time (for example, 24 hours) elapses, and when the certain time elapses, average values of np and nq are calculated for daytime and nighttime, respectively. In the correction step of active power and reactive power shown in FIG. 2, correction is performed using np and nq at daytime and nighttime.
[0049]
In the above example, an example of performing both the time-series moving average process and the rated voltage correction is shown. However, since these processes are independent of each other, even if only one of the processes is performed, the power The estimation accuracy of the coefficient of the above quadratic equation representing the characteristics of the system can be increased. For example, when performing characteristic estimation for a system with a small load power factor variation such as a residential load system, the effect is great even if only the rated voltage correction is performed without performing the moving average process. Examples thereof are shown in FIGS. 8 and 9 are a plot example and a PQ curve on the PQ plane when measuring a feeder with a rated voltage of 66 kV for a residential area. FIG. 9 is a case where it went. Here, the load is corrected with a constant current characteristic for active power, that is, np = 1, and with respect to reactive power, with a constant impedance characteristic, that is, nq = 2. As shown in both figures, in this example, the C term obtained without the rated voltage correction is +7.492 Mvar, whereas the obtained C term after the correction is +9.454 Mvar, which varies greatly. I understand that. In other words, it can be seen that the characteristic estimation accuracy increases accordingly.
[0050]
Next, an embodiment will be described in which actual P-Q correlation data P, Q, and V are rated voltage corrected to obtain a highly accurate C term and reactive power load characteristic constant nq.
[0051]
The configuration of the apparatus is the same as that shown in FIG. FIG. 12 is a flowchart showing a processing procedure for estimating the C term and nq.
[0052]
The operation is the same as that shown in FIG. 2 until the active power P and the reactive power Q are calculated, and time series real PQ correlation data P, Q, and V are calculated by repeating several times. Then, the following procedure is repeated at an appropriate timing for outputting the estimated values of the nq and C terms.
[0053]
(1) P′-Q ′ correlation data after rated voltage correction is obtained from actual PQ correlation data (see equations (5) and (7)). However, initially, since the C term (C ′) and nq in the equation are unknown numbers, appropriate values (here, C ′ = 0, nq = 0, 10) are set.
[0054]
(2) A P′-Q ′ correlation curve is drawn.
[0055]
(3) Obtain the C term (C ′) and use it as an estimated value.
[0056]
(4) Using the voltage step change data (V ₀ → V ₁ ), nq is obtained and used as an estimated value (see equation (3)). The voltage step change data uses change data V ₀ and V ₁ when a step change (disturbance) occurs in the measured (actual) voltage.
[0057]
(5) The above (1) to (4) are repeated until nq converges to a predetermined level (here, 0.95 <(nq _i / nq _{i + 1} ) <1.05).
[0058]
(6) At the stage of convergence to a certain level, the nq and C term (C ′) at that time are output.
[0059]
When the P′-Q ′ correlation curve is drawn to estimate the C term (C ′), the moving average process as shown in FIG. 2 can be performed.
[0060]
Simulation results are shown in FIGS. True np In this simulation, when there and without fluctuation of the bank secondary voltage as 2,4 respectively nq, and if there a case where there is no bank Z _s, and if there and if there is no line impedance Z _L Each simulation result is shown. The initial nq is set to 0 and 10. As shown in these simulation results, it can be seen that nq converges before nq is estimated at least five times. That is, by repeatedly calculating nq in PQ correlation data with variations in voltage, the C term and nq can be obtained with higher accuracy than when voltage correction is not performed.
[0061]
【The invention's effect】
According to the first and third aspects of the present invention, the active power and reactive power of the power system are repeatedly obtained at a constant period, and each of them is smoothed. When the coefficients of the equations representing their relationship are estimated from multiple smoothed active and reactive power pairs, these coefficients can be treated as more accurately representing the actual characteristic values of the power system. .
[0062]
According to the invention described in claim 2 and claim 4, the active power and the reactive power equivalent to those obtained under the rated voltage are always obtained, and from these corrected active power and reactive power pairs, Since the coefficient of the equation indicating the relationship is obtained, the accuracy of the power system characteristic estimation is improved even if the voltage fluctuates.
[0063]
According to the fifth aspect of the present invention, it is possible to more accurately estimate the characteristics of the power system even when both the load power factor fluctuation and the voltage fluctuation are severe.
[0064]
According to the sixth aspect of the present invention, when there is a disturbance such as a step change in voltage, the load characteristic constants np and nq of the active power and reactive power are automatically calculated from the active power, reactive power and voltage before and after the disturbance. Therefore, special measurement for obtaining the load characteristic constants of active power and reactive power is not required.
[0065]
According to the seventh aspect of the present invention, even when there is a voltage fluctuation in the power system, it is possible to obtain highly accurate nq and C terms by repeatedly calculating nq, and accurately grasp the characteristics of the system load. can do.
[Brief description of the drawings]
FIG. 1 is a block diagram illustrating a configuration of a load characteristic recording apparatus.
FIG. 2 is a flowchart showing a processing procedure for sampling and characteristic value calculation;
FIG. 3 is a flowchart showing a procedure for calculating a load characteristic constant of active power / reactive power.
FIG. 4 is a diagram illustrating an example of transition of active power and reactive power over time.
FIG. 5 is a diagram illustrating an example of a plot on a PQ plane.
FIG. 6 is a diagram illustrating an example of load characteristic constant calculation timing.
FIG. 7 is a diagram illustrating an example of a change in a PQ curve before and after performing rated voltage correction using a load characteristic constant.
FIG. 8 is a diagram illustrating an example of a PQ curve before a rated voltage correction by a load characteristic constant.
FIG. 9 is a diagram illustrating an example of a PQ curve after the rated voltage is corrected by a load characteristic constant.
FIG. 10 is a model diagram of the power system when the bank secondary side and below are reduced to 66 kV.
FIG. 11 is a diagram for explaining a method of obtaining C ′ by correcting rated voltage of PQ correlation data.
FIG. 12 is a flowchart showing a procedure for estimating nq and C ′ by convergence by iterative calculation.
FIG.
FIG. 20 is a diagram illustrating simulation results.
FIG. 21 is an equivalent circuit diagram of a load system.
FIG. 22 is a diagram illustrating an example of a PQ curve.
FIG. 23 is an equivalent circuit diagram of the variable load system.
FIG. 24 is a diagram illustrating an example of transition of active power and reactive power over time.
FIG. 25 is a diagram illustrating an example of a plot on a PQ plane.
[Explanation of symbols]
100-terminal device 101-monitoring device

Claims

While obtaining the time series data of the active power and reactive power of the power system, and means for smoothing these time series data for the active power and reactive power respectively, and a plurality of pairs of the smoothed active power and reactive power, A power system characteristic estimation device comprising means for obtaining coefficients of equations indicating these relationships as characteristic values in the power system.

Obtain time series data of power system voltage, active power and reactive power, and set active power, reactive power and voltage by measurement as P, Q and V, rated voltage as V ', load characteristic constant of active power, Assuming that the reactive power load characteristic constants are np and nq, respectively, P ′ = P (V ′ / V) ^np after correction, and Q ′ = Q (V ′ / V) ^nq And a means for obtaining each of the corrected reactive power Q ′ and a means for obtaining the coefficient of an equation representing the relationship as a characteristic value in the power system from a plurality of pairs of the corrected active power and reactive power. Power system characteristic estimation device.

Obtain time series data of the active power and reactive power of the power system, smooth the time series data for active power and reactive power, respectively, and the relationship between these smoothed active power and reactive power pairs. The characteristic estimation method of the electric power system which calculates | requires the coefficient of the equation which shows as a characteristic value in an electric power system.

Obtain time series data of power system voltage, active power and reactive power, and set active power, reactive power and voltage by measurement as P, Q and V, rated voltage as V ', load characteristic constant of active power, Assuming that the reactive power load characteristic constants are np and nq, respectively, P ′ = P (V ′ / V) ^np after correction, and Q ′ = Q (V ′ / V) ^nq Power system characteristic estimation method for obtaining each of the corrected reactive power Q 'shown and obtaining coefficients of equations representing these relationships as characteristic values in the power system from a plurality of corrected active power and reactive power pairs .

Obtain time series data of power system voltage, active power and reactive power, and set active power, reactive power and voltage by measurement as P, Q and V, rated voltage as V ', load characteristic constant of active power, When the reactive power load characteristic constants are np and nq, respectively, the corrected active power P ′ and Q ′ = Q (V ′ / V) ^nq indicated by P ′ = P (V ′ / V) ^np Each corrected reactive power Q ′ is obtained, time series data of the corrected active power and reactive power is smoothed, and the relationship between these smoothed active power and reactive power pairs is shown. A method for estimating characteristics of a power system in which the coefficient of the equation is obtained as a characteristic value in the power system.

Obtain time series data of the active power and reactive power of the power system, detect the disturbance of the voltage of the power system, and set the active power, reactive power and voltage before the disturbance as Po, Qo, Vo, and the effective power after the disturbance. And when the reactive power and voltage are P1, Q1, and V1, respectively, P1 = Po (V1 / Vo) ^np indicated by np and Q1 = Qo (V1 / Vo) ^nq indicated by nq are obtained, respectively. Or a power system characteristic estimation method applied to each of the active power load characteristic constant np and the reactive power load characteristic constant nq described in 5.

The reactive power load characteristic constant nq is estimated from the time series data of the active power and reactive power of the power system and the voltage step change data by the following steps (1) to (3). Characteristic estimation method.
(1) P ′ = P (V ′ / V) ^np
Q ′ = {Q + C ′ (V / V ′) ² } (V ′ / V) ^nq −C ′
Then, (P, Q, V) is corrected at the rated voltage (V ') to estimate C'. However, initially C ′ and nq are arbitrary values, and np is
P ₁ = P ₀ (V ₁ / V ₀ ) ^Obtained from ^np .
(2) nq is estimated from the following equation using C ′ and voltage step change data estimated in (1) above.
Q ₁ = (Q ₀ + C ′) (V ₁ / V ₀ ) ^nq −C ′ (V ₁ / V ₀ ) ²
(3) Using the C ′ estimated in the above (1) and the nq estimated in (2), (1) is repeated a predetermined number of times to obtain nq.
However,
P, Q, V: Measurement (actual) data (active power, reactive power, voltage)
P ′, Q ′, V ′, C ′: Correction data at rated voltage (active power, reactive power, voltage, load system capacity value (C term))
P ₀ , Q ₀ , V ₀ : Data before voltage step change (active power, reactive power, voltage)
P ₁ , Q ₁ , V ₁ : Data after voltage step change (active power, reactive power, voltage)
np: Load characteristic constant of active power nq: Load characteristic constant of reactive power