JP3719477B2

JP3719477B2 - Tidal current calculation method in radial distribution system.

Info

Publication number: JP3719477B2
Application number: JP05967998A
Authority: JP
Inventors: 良和福山; 洋大井
Original assignee: Fuji Electric Systems Co Ltd
Current assignee: Fuji Electric Co Ltd
Priority date: 1998-03-11
Filing date: 1998-03-11
Publication date: 2005-11-24
Anticipated expiration: 2018-03-11
Also published as: JPH11262180A

Description

【０００１】
【発明の属する技術分野】
本発明は、配電営業所の制御用コンピュータ等に付加され、放射状配電系統における電源容量、負荷容量、線路データ、変圧器データ、自動電圧調整器データ等が与えられた時に、三相潮流計算を高速で実行するための潮流計算方法に関する。
【０００２】
【従来の技術及び発明が解決しようとする課題】
従来、配電系統では、回路計算を用いた解析が多く行われてきたが、分散電源の導入、営業諸機能の高度化等により潮流計算による解析が必要になりつつある。特に、放射状配電系統は本質的に単相負荷が多く、また、末端部では不平衡負荷が多くなってきており、インバータエアコン等に代表される定電力負荷の増加により、従来多く利用されている簡略回路計算による損失計算や電圧制約条件では実態と合わなくなってきている。
【０００３】
このため、三相不平衡を考慮するとともに、定電力負荷等の各種の負荷特性を考慮した三相不平衡潮流計算を高速で行うことが要請されている。
一方、自動電圧調整器（ＳＶＲ）は二次側電圧の絶対値が既定値になるようにタップを自動調整する機能を有しているが、従来の潮流計算では系統のＳＶＲの存在を考慮しておらず、そのため、潮流計算の精度が正確性を欠き、実用性の点で問題があった。
【０００４】
そこで本発明は、系統に存在するＳＶＲが自動的にタップを調整する機能を定式化し、ＳＶＲの二次側電圧規定値から逆に求めたタップを用いて線路のＰ,Ｑ値を求めることにより、高精度な三相不平衡潮流計算を可能にした潮流計算方法を提供しようとするものである。
【０００５】
【課題を解決するための手段】
上記課題を解決するため、請求項１記載の発明は、自動電圧調整器が存在する放射状系統を、ブランチと負荷とをノードとするネットワークにより表現し、各ノードにおける有効電力、無効電力、電圧を状態方程式により求める三相潮流計算方法において、配電系統の電源電圧、負荷容量、線路データ、変圧器データ、自動電圧調整器データを入力するステップと、各分岐線の末端ノードから先頭ノードに向けて負荷量の総和を求めていき、先頭ノードへの流入電力を状態変数の初期値として計算するステップと、前記ブランチが送電線または変圧器である場合には、着目ノードの状態量としての有効電力、無効電力、電圧を、上流ノードにおける流出電力、電流、電圧及び着目ノードにおける負荷電力から求め、前記ブランチが自動電圧調整器である場合には、着目ノードの状態量としての有効電力、無効電力、電圧を、上流ノードの有効電力、無効電力、電圧の関数である状態方程式と、自動電圧調整器の二次側電圧規定値から求めたタップの値とを用いて各相ごとに求める前進計算ステップと、各分岐線の末端ノードにおける流出電力を先頭ノードにおける流入電力から差し引いて前記状態変数を更新する後進計算ステップと、各分岐線の末端ノードにおける流出電力がしきい値以下になったことを判定して計算を収束させる収束判定ステップと、を有するものである。
【０００７】
【発明の実施の形態】
以下、本発明の実施形態を説明する。
１．系統モデル
ブランチと負荷とをノードとする図１のネットワーク表現を用いる。ブランチモデルとして図示するようなπ形モデルを用いることにより、以下に示すように送電線や変圧器を取り扱うことが可能になる。
なお、図１において、Ｐ＋ｊＱは上流から流入する電力、ｒ＋ｊｘは線路のインピーダンス、Ｐ_L+ｊＱ_Lは定電力負荷、ａ_L+ｊｂ_Lは定電流負荷、Ｒ_L+ｊＸ_Lは定インピーダンス負荷、Ｙ₁,Ｙ₂は線路のアドミタンスを示す。
【０００８】
（１）送電線
送電線には、図２に示すような一般的なπ形モデルを用いる。図２において、Ｙは送電線のアドミタンスである。
【０００９】
（２）変圧器
変圧器には、図３に示すπ形モデルを用いる。図３における各値は、数式１〜数式３に示すとおりである。
【００１０】
【数１】
Ｚ₁=Ｚ_L/ｎ
【００１１】
【数２】
Ｙ₁=ｎ（ｎ−１）／Ｚ_L
【００１２】
【数３】
Ｙ₂=（１−ｎ）／Ｚ_L
【００１３】
なお、数式１〜数式３において、Ｚ_Lは変圧器の漏れインダクタンス、ｎは変成比（一次側を１．０とした時の二次側タップ比）である。
【００１４】
（３）負荷
負荷としては、定電力負荷（Ｓ＝Ｐ_L+ｊＱ_L:固定値）、定電流負荷（Ｉ＝ａ_L+ｊｂ_L:固定値）、定インピーダンス負荷（ｒ_L+ｊｘ_L:固定値）の３種類を取り扱うこととする。
定電力負荷の負荷電力Ｓ_kは数式４により、定電流負荷の負荷電力Ｓ_kは数式５により、定インピーダンス負荷の負荷電力Ｓ_kは数式６により、それぞれ表される。
【００１５】
【数４】

【００１６】
【数５】

【００１７】
【数６】

【００１８】
なお、数式５におけるｅ_k,ｆ_k,ａ_k,ｂ_kは電圧、電流ベクトルの実数部、虚数部を示し、また、数式５、数式６における「＊」は共役を示す（以下、同じ）。
【００１９】
（４）分散電源
系統に存在する分散電源は、一般に定電力を発生する機器として定義できることから、等価的に負の定電力負荷として取り扱うこととする。
【００２０】
２．状態方程式
ここでは、各ブランチに流入する電力及びノード電圧を状態量とする状態方程式を考える。
図４の２ノードによる系統において、ブランチが送電線である場合には、ノード２における有効電力Ｐ₂、無効電力Ｑ₂、電圧Ｖ₂,｜Ｖ₂|²に関して、以下の数式７〜数式１０が成り立つ。なお、Ｙ₁,Ｙ₂は線路充電容量のため、容量性と仮定する。
【００２１】
【数７】

【００２２】
この数式７において、Ｐ₁はノード１に流入する有効電力、Ｑ₁は同じく無効電力、Ｐ_loss,1はノード１における有効電力損失、Ｐ_L2はノード２における負荷の有効電力である。ここで、Ｐ₂はＰ₁,Ｑ₁,｜Ｖ₁|²の関数ｆ_Pによって表される。
【００２３】
【数８】

【００２４】
この数式８において、Ｑ_loss,1はノード１における無効電力損失、Ｑ_L2はノード２における負荷の無効電力である。
ここで、Ｑ₂もＰ₁,Ｑ₁,｜Ｖ₁|²の関数ｆ_Qによって表される。
【００２５】
【数９】

【００２６】
なお、数式９において、ｅ₁,ｅ₂はＶ₁のベクトルの実数部及び虚数部、ｆ₁,ｆ₂はＶ₂のベクトルの実数部及び虚数部である。
【００２７】
【数１０】

【００２８】
この数式１０において、｜Ｖ₂|²もＰ₁,Ｑ₁,｜Ｖ₁|²の関数ｆ_Vによって表される。
【００２９】
つまり、ノード１における有効電力Ｐ₁、無効電力Ｑ₁、電圧Ｖ₁を用いて、ノード２における状態量としての有効電力Ｐ₂、無効電力Ｑ₂、及び｜Ｖ₂|²を上記のように表すことができる。分岐がある場合には、上述した状態方程式に分岐線に流れる潮流分を考慮すればよい。
【００３０】
ブランチが変圧器である場合には、図５及び前述の数式１〜数式３を用いることにより、有効電力Ｐ₂、無効電力Ｑ₂、電圧Ｖ₂及び｜Ｖ₂|²は数式１１〜数式１４によって表される。ここで、Ｙ₁,Ｙ₂は容量性であると仮定したため、実際のＹ₁,Ｙ₂を数式２、数式３から代入する場合にはマイナスを付けることに注意する必要がある。
【００３１】
【数１１】

【００３２】
【数１２】

【００３３】
【数１３】

【００３４】
【数１４】

【００３５】
この場合にも、Ｐ₂,Ｑ₂,｜Ｖ₂|²は、Ｐ₁,Ｑ₁,｜Ｖ₁|²の関数ｆ_P,ｆ_Q,ｆ_Vによって表される。
以上のことから、電源電圧は一定であると仮定すると、放射状系統においては、上流と下流ノードとの関係を考慮すると、電源から流れる電力及び各分岐線への流出（流入）電力を状態変数とした潮流計算を行えばよく、他の変数については状態方程式を用いて電源端から逐次計算すれば良い。通常の高圧用潮流計算ではすべてのノードの状態量を状態変数とするの対し、本発明の方法では状態変数を配電線から各分岐線への流出（流入）電力とすることにより、状態変数を大幅に減少させて高速な潮流計算を実現することができる。
【００３６】
次に、系統に存在するＳＶＲの考慮方法につき説明する。
ＳＶＲでは、二次側電圧の絶対値が規定値になるようにタップを自動調整する。まず、上記タップが連続的に調整されると仮定する。
前記数式１、数式２、数式３及び数式１４から、数式１５が得られる。この数式１５において、｜Ｖ₂|はＳＶＲの二次側電圧規定値（絶対値）、｜Ｖ₁|は一次側電圧計算値（絶対値）である。
【００３７】
【数１５】

【００３８】
数式１５において、変成比ｎはｎ＞０であることから、数式１６が得られる。
【００３９】
【数１６】

【００４０】
ここで、潮流計算においては、数式１６の計算は単位法によって表現され、通常の運用状態を考えると、｜Ｖ₁|≒１.０、｜Ｖ₂|≒１.０、Ｚ_L≪１.０となる。従って、数式１６では、ｎ≒１.０、ｎ≒０.０という二つの解が求められる。このため、ｎには１に近い値を採用する。
実際に採り得るタップ値のうち、上記のｎに最も近いタップ値を用いて計算した二次側電圧が電圧規定値を超えるようなタップ値をｎ₁として採用する。そして、このタップ値ｎ₁を数式１〜数式３を経て数式１１、数式１２に代入することにより、Ｐ₂,Ｑ₂を計算する。
【００４１】
次いで、ＳＶＲを考慮した単相潮流計算方法について述べる。
上述したような方法により、上流側の状態量Ｐ₁,Ｑ₁,Ｖ₁に基づいて下流側の状態量Ｐ₂,Ｑ₂,Ｖ₂を計算すること（前進計算）が可能になる。
また、系統の末端ノードでは流出電力はない（Ｐ_n=Ｑ_n=０）という仮定から、状態変数を変更すること（後退計算）が可能である。
初期値として与える先頭ノードの初期状態量は、以下のように設定する。
ａ．電圧：設定する電源電圧とする。
ｂ．潮流：末端ノードからすべての負荷量の総和を求めて初期値とする。これは損失を考慮していないため、実際の潮流値よりも多少、少な目になる。
従って、末端ノードの電力値（Ｐ_n,Ｑ_n)分だけ先頭ノードの状態量を変化させ、末端電力がしきい値以下（例えば零）になるまで収束計算を行う。一般の分岐線のある放射状系統では、配電線の状態量により各分岐線の状態量も変更されるため、後退計算は、末端の分岐線から行う。
【００４２】
上述の前進計算、後退計算の概念を図６に示す。図６（ａ）の前進計算において、▲１▼の経路で上流側状態量から下流側状態量を求め、この経路から▲２▼で分岐した経路▲３▼に付き、同様に下流側状態量を求める。以下、同様にして▲４▼，▲５▼と計算を進める。
図６（ｂ）の後退計算では、▲１▼の経路の末端ノードの電力値に基づいて先頭ノードの状態変数を変化させながら収束計算を行い、経路▲２▼を経て分岐元の経路▲３▼について、同様に末端ノードの電力値に基づいて先頭ノードの状態変数を更新しながら収束計算を行なう。以下、同様にして▲４▼，▲５▼と計算を進める。
この方法は、状態変数が大幅に少なくなることから、すべてのノードの状態変数を用いる従来の方法に比べて高速計算が可能になる。
【００４３】
従って、単相潮流計算のアルゴリズムは以下のようになる。
（１）ステップ１：系統データの入力
電源電圧、負荷容量、線路データ、変圧器データ（定格電圧やタップ数、タップ幅等）、ＳＶＲデータ（二次側電圧規定値やタップ数、タップ幅等）を入力する。
（２）ステップ２：状態変数の初期値計算
各末端ノードから負荷量の総和を求めていき、各分岐線への流入電力（状態変数）の初期値をその分岐線以下の負荷量の総和とする。
（３）ステップ３：前進計算
ブランチが送電線であれば、前記数式７〜数式１０を用いて状態量を計算する。
ブランチが変圧器であれば、前記数式１１〜数式１４を用いて状態量を計算する。
ブランチがＳＶＲであれば、前記数式１１〜数式１４、数式１６を用いて状態量を計算する。
（４）ステップ４：後退計算
各末端ノードの流出電力分だけ各分岐線の先頭ノードへの流入電力から差し引くことにより、状態変数を更新する。
（５）ステップ５：収束判定
各末端ノードの電力値がしきい値以下になったら収束したと判定し、そうでなければ前記ステップ３以後の処理を繰り返す。
【００４４】
次に、三相潮流計算方法を以下に述べる。
三相潮流計算では、前記ステップ３における前進計算が数式１７のベクトル式となる。なお、図７は三相潮流計算を行う場合の２ノード系統を示している。この図において、Ｖ_kは上流ノードの電圧、Ｓ_kは上流ノードの流出電力、Ｉ_kは同じく電流、Ｚ_kは線路のインピーダンス、Ｖ_k+1は下流ノードの電圧、Ｓ_k+1は同じく流出電力、Ｓ_Lk+1は同じく負荷電力である。
【００４５】
【数１７】

【００４６】
Ｚ_kが送電線、変圧器である場合、それぞれのインピーダンス行列を代入する。ＳＶＲの場合には、各相ごとに前記数式１１〜数式１４、数式１６を用いて状態量を計算する。
【００４７】
最後に、本発明の参考例として、図８に示す単相回路に対する潮流計算例を以下に示す。図中、すべてのブランチ（ＳＶＲを除く）及びノードにおいて線路インピーダンスは０．００８＋ｊ０．０６６、定電力負荷量は０．０５４＋ｊ０．０２６１６となっていて同一であるとする。また、ノード１及びノード５における電圧は１．０（単位法）とする。
各ノードにおいて下流ノードに向かう有効電力Ｐ、無効電力Ｑは数式１８のようになる。なお、数式１８において、V＾２はノード電圧絶対値の二乗、｜V｜は線路電圧絶対値を示す。
【００４８】
【数１８】

【００４９】
この参考例の系統に対しては、８回の反復により収束している。また、ノード４，５間のＳＶＲについては、タップの解として１．０付近の解が１．１１８８８４、０．０付近の解が０．０００６４０となっており、両者のうち１．１１８８８４が解となる。
【００５０】
【発明の効果】
以上のように本発明では、各分岐線への流入電力を状態変数とし、他の状態変数については先頭ノードから逐次計算して求めると共に、ＳＶＲが自動的にタップを調整する機能を定式化し、二次側電圧規定値から求めたタップを用いて着目ノードの有効電力、無効電力を求めているので、状態変数が大幅に減少することとなって高速な潮流計算が可能になり、また、実系統に存在するＳＶＲを考慮した実用的な潮流計算を行うことができる。
【図面の簡単な説明】
【図１】系統モデルの説明図である。
【図２】送電線モデルの説明図である。
【図３】変圧器モデルの説明図である。
【図４】２ノード系統の送電線ブランチの説明図である。
【図５】２ノード系統の変圧器ブランチの説明図である。
【図６】前進計算、後退計算の概念の説明図である。
【図７】三相潮流計算に用いる２ノード系統の説明図である。
【図８】本発明の参考例における単相潮流計算用系統の説明図である。[0001]
BACKGROUND OF THE INVENTION
The present invention is added to a control computer or the like of a distribution office, and when a power capacity, load capacity, line data, transformer data, automatic voltage regulator data, etc. in a radial distribution system are given , a three- phase power flow calculation is performed. The present invention relates to a power flow calculation method for high-speed execution.
[0002]
[Prior art and problems to be solved by the invention]
Conventionally, in power distribution systems, many analyzes using circuit calculations have been performed. However, analysis based on power flow calculations is becoming necessary due to the introduction of distributed power supplies and the advancement of various sales functions. In particular, the radial power distribution system has essentially a single-phase load, and an unbalanced load has increased at the end, and has been widely used due to the increase in constant power load represented by inverter air conditioners and the like. Loss calculation by simplified circuit calculation and voltage constraint conditions are not consistent with the actual situation.
[0003]
For this reason, it is required to perform three-phase unbalanced power flow calculation at high speed in consideration of three-phase unbalance and various load characteristics such as constant power load.
On the other hand, the automatic voltage regulator (SVR) has a function to automatically adjust the tap so that the absolute value of the secondary side voltage becomes a predetermined value, but the conventional power flow calculation considers the existence of the SVR of the system. Therefore, the accuracy of the tidal current calculation is inaccurate, and there is a problem in terms of practicality.
[0004]
Therefore, the present invention formulates the function of automatically adjusting the tap by the SVR existing in the system, and calculates the P and Q values of the line by using the tap obtained reversely from the secondary voltage regulation value of the SVR. It is intended to provide a tidal current calculation method that enables highly accurate three- phase unbalanced tidal current calculations.
[0005]
[Means for Solving the Problems]
In order to solve the above-mentioned problem, the invention according to claim 1 expresses a radial system in which an automatic voltage regulator exists by a network having a branch and a load as nodes, and the active power, reactive power, and voltage at each node are expressed. In the three- phase power flow calculation method obtained from the state equation, the step of inputting the power supply voltage, load capacity, line data, transformer data, and automatic voltage regulator data of the distribution system, and from the end node of each branch line to the top node The step of calculating the total load amount and calculating the inflow power to the head node as the initial value of the state variable, and the active power as the state amount of the node of interest when the branch is a transmission line or a transformer , reactive power, voltage, outlet power in the upstream node, current, voltage, and determined from the load power at the target node, the branch automatic voltage regulation The active power, reactive power, and voltage as state quantities of the node of interest, the state equation that is a function of the active power, reactive power, and voltage of the upstream node and the secondary voltage of the automatic voltage regulator. A forward calculation step for each phase using a tap value obtained from a specified value; a backward calculation step for updating the state variable by subtracting the outflow power at the end node of each branch line from the inflow power at the top node; And a convergence determination step for determining that the outflow power at the end node of each branch line is equal to or less than a threshold value and for converging calculation.
[0007]
DETAILED DESCRIPTION OF THE INVENTION
Embodiments of the present invention will be described below.
1. The network representation of FIG. 1 having the system model branch and the load as nodes is used. By using a π-shaped model as illustrated as a branch model, it becomes possible to handle a transmission line and a transformer as shown below.
In FIG. 1, P + jQ is power flowing in from upstream, r + jx is line impedance, P _L + jQ _L is constant power load, a _L + jb _L is constant current load, R _L + jX _L is constant impedance load, Y ₁ and Y ₂ indicate line admittances.
[0008]
(1) Transmission line A general π-shaped model as shown in FIG. 2 is used for the transmission line. In FIG. 2, Y is the admittance of the power transmission line.
[0009]
(2) Transformer The π-type model shown in FIG. 3 is used for the transformer. Each value in FIG. 3 is as shown in Equations 1 to 3.
[0010]
[Expression 1]
Z ₁ = Z _L / n
[0011]
[Expression 2]
Y ₁ = n (n−1) / Z _L
[0012]
[Equation 3]
Y ₂ = (1-n) / Z _L
[0013]
In Equations 1 to 3, Z _L is a transformer leakage inductance, and n is a transformation ratio (secondary side tap ratio when the primary side is 1.0).
[0014]
(3) As load loads, constant power load (S = P _L + jQ _L : fixed value), constant current load (I = a _L + jb _L : fixed value), constant impedance load (r _L + jx _L : (Fixed value) will be handled.
The load power S _k The formula 4 of a constant power load, the load power S _k of the constant current load by Equation 5, the load power S _k of the constant impedance load by Equation 6, as represented respectively.
[0015]
[Expression 4]

[0016]
[Equation 5]

[0017]
[Formula 6]

[0018]
In Equation 5, e _k , f _k , a _k , and b _k represent the real part and imaginary part of the voltage and current vectors, and “*” in Equation 5 and Equation 6 represents the conjugate (hereinafter the same). .
[0019]
(4) Since a distributed power source existing in a distributed power system can be generally defined as a device that generates constant power, it is treated as a negative constant power load equivalently.
[0020]
2. State Equation Here, a state equation is considered in which the power and node voltage flowing into each branch are state quantities.
In systems with two nodes in FIG 4, when the branch is the transmission line, the effective power P ₂ at the node 2, the reactive power Q _2, the voltage V _2, | V ₂ | respect ^2, Equation 7 Equation 10 below Holds. Y ₁ and Y ₂ are assumed to be capacitive because of the line charging capacity.
[0021]
[Expression 7]

[0022]
In Equation 7, P ₁ is the active power flowing into the node 1, Q ₁ is the reactive power, P _{loss, 1} is the active power loss at the node 1, and P _L2 is the active power of the load at the node 2. Here, P ₂ is represented by a function f _P of P ₁ , Q ₁ , | V ₁ | ² .
[0023]
[Equation 8]

[0024]
In Equation 8, Q _{loss, 1} is a reactive power loss at node 1 and Q _L2 is a reactive power of a load at node 2.
Here, Q _{2 is} also expressed by a function f _{Q of} P ₁ , Q ₁ , | V ₁ | ² .
[0025]
[Equation 9]

[0026]
Note that in Equation 9, e _1, e ₂ are the real and imaginary parts of the vector of V _{_1,} f _1, f ₂ is the real part and the imaginary part of the vector of V _2.
[0027]
[Expression 10]

[0028]
In this equation 10, | V ₂ | ² also P _{_1,} Q _1, | is represented by ² function f _V | _V _1.
[0029]
That is, using the active power P ₁ , reactive power Q ₁ , and voltage V ₁ at node ₁ , active power P ₂ , reactive power Q ₂ , and | V ₂ | ² as state quantities at node 2 are as described above. Can be represented. In the case where there is a branch, the current flowing in the branch line may be taken into consideration in the above-described state equation.
[0030]
When the branch is a transformer, the active power P ₂ , the reactive power Q ₂ , the voltage V ₂ and | V ₂ | ² are expressed by Expressions 11 to 14 by using FIG. 5 and Expressions 1 to 3 described above. Represented by Here, since it is assumed that Y ₁ and Y ₂ are capacitive, it is necessary to pay attention to the fact that the actual Y ₁ and Y ₂ are assigned with a minus when substituting from the

expressions

2 and 3.
[0031]
[Expression 11]

[0032]
[Expression 12]

[0033]
[Formula 13]

[0034]
[Expression 14]

[0035]
In this _{_{case, P 2, Q 2, |}} V 2 | 2 _{_{is, P 1, Q 1, |}} V 1 | 2 function f _P, f _Q, represented by f _V.
From the above, assuming that the power supply voltage is constant, in the radial system, considering the relationship between the upstream and downstream nodes, the power flowing from the power supply and the outflow (inflow) power to each branch line are used as state variables. The other tidal current calculations may be performed, and other variables may be calculated sequentially from the power supply end using the state equation. In the normal high-pressure power flow calculation, the state variables of all nodes are used as state variables, whereas in the method of the present invention, the state variables are set to the outflow (inflow) power from the distribution line to each branch line. It can be greatly reduced and high-speed tidal current calculation can be realized.
[0036]
Next, a method for considering SVR existing in the system will be described.
In SVR, the tap is automatically adjusted so that the absolute value of the secondary side voltage becomes a specified value. First, assume that the tap is continuously adjusted.
From Equation 1, Equation 2, Equation 3, and Equation 14, Equation 15 is obtained. In Equation 15, | V ₂ | is a secondary voltage regulation value (absolute value) of SVR, and | V ₁ | is a primary voltage calculation value (absolute value).
[0037]
[Expression 15]

[0038]
In Formula 15, since the transformation ratio n is n> 0, Formula 16 is obtained.
[0039]
[Expression 16]

[0040]
Here, in the tidal current calculation, the calculation of Expression 16 is expressed by the unit method. Considering the normal operation state, | V ₁ | ≈1.0, | V ₂ | ≈1.0, Z _L << 1. 0. Therefore, in Equation 16, two solutions of n≈1.0 and n≈0.0 are obtained. For this reason, a value close to 1 is adopted as n.
Of the tap values that can be actually taken, a tap value such that the secondary side voltage calculated using the tap value closest to n described above exceeds the specified voltage value is adopted as n ₁ . Then, P ₂ and Q ₂ are calculated by substituting this tap value n ₁ into Equations 11 and 12 via Equations 1 to 3.
[0041]
Next, a single-phase power flow calculation method considering SVR will be described.
By the method as described above, it is possible to calculate the downstream state quantities P ₂ , Q ₂ , V ₂ based on the upstream state quantities P ₁ , Q ₁ , V ₁ (forward calculation).
Further, it is possible to change the state variable (reverse calculation) from the assumption that there is no outflow power at the end node of the system (P _n = Q _n = 0).
The initial state quantity of the head node given as the initial value is set as follows.
a. Voltage: The power supply voltage to be set.
b. Tidal current: Calculate the sum of all loads from the end node and use it as the initial value. This is a little less than the actual tidal current value because loss is not taken into account.
Therefore, the state quantity of the leading node is changed by the power value (P _n , Q _n ) of the terminal node, and convergence calculation is performed until the terminal power is equal to or less than a threshold value (for example, zero). In a radial system with a general branch line, the state quantity of each branch line is also changed depending on the state quantity of the distribution line, so the backward calculation is performed from the terminal branch line.
[0042]
The concept of the above-described forward calculation and backward calculation is shown in FIG. In the forward calculation in FIG. 6 (a), the downstream state quantity is obtained from the upstream state quantity in the path (1), and the downstream state quantity is similarly attached to the path (3) branched from this path in (2). Ask for. In the same manner, the calculation proceeds as (4) and (5).
In the backward calculation of FIG. 6B, convergence calculation is performed while changing the state variable of the leading node based on the power value of the terminal node of the route {circle around (1)}, and the branch source route {circle around (3)} via the route {circle around (2)}. Concerning ▼, convergence calculation is similarly performed while updating the state variable of the leading node based on the power value of the terminal node. In the same manner, the calculation proceeds as (4) and (5).
Since this method significantly reduces the number of state variables, high-speed calculation is possible as compared with the conventional method using the state variables of all nodes.
[0043]
Therefore, the algorithm for single-phase power flow calculation is as follows.
(1) Step 1: System data input power supply voltage, load capacity, line data, transformer data (rated voltage, number of taps, tap width, etc.), SVR data (secondary voltage specified value, number of taps, tap width, etc.) ).
(2) Step 2: Calculation of initial values of state variables The sum of load amounts is obtained from each end node, and the initial value of the inflow power (state variable) to each branch line is calculated as the sum of load amounts below that branch line. To do.
(3) Step 3: If the forward calculation branch is a power transmission line, the state quantity is calculated using Equation 7 to Equation 10.
If the branch is a transformer, the state quantity is calculated using Equations 11 to 14.
If the branch is an SVR, the state quantity is calculated using Equations 11 to 14 and 16.
(4) Step 4: Backward calculation The state variable is updated by subtracting from the inflow power to the head node of each branch line by the outflow power of each end node.
(5) Step 5: Convergence determination If the power value of each terminal node is less than or equal to the threshold value, it is determined that the power has converged, and if not, the processing after step 3 is repeated.
[0044]
Next, the three-phase power flow calculation method is described below.
In the three-phase power flow calculation, the forward calculation in Step 3 is a vector expression of Expression 17. FIG. 7 shows a two-node system when performing three-phase power flow calculation. In this figure, V _k is the voltage of the upstream node, S _k is the outflow power of the upstream node, I _k is the same current, Z _k is the impedance of the line, V _{k + 1} is the voltage of the downstream node, and S _{k + 1} is the same. The outflow power, S _{Lk + 1,} is also load power.
[0045]
[Expression 17]

[0046]
When Z _k is a transmission line or a transformer, the respective impedance matrices are substituted. In the case of SVR, the state quantity is calculated for each phase by using the equations 11 to 14 and 16.
[0047]
Finally, as a reference example of the present invention, a power flow calculation example for the single-phase circuit shown in FIG. 8 is shown below. In the figure, it is assumed that the line impedance is 0.008 + j0.066 and the constant power load amount is 0.054 + j0.02616 in all branches (except SVR) and nodes. The voltage at node 1 and node 5 is 1.0 (unit method).
The active power P and the reactive power Q that are directed to the downstream node in each node are expressed by Equation 18. In Equation 18, V ^ 2 is the square of the node voltage absolute value, and | V | is the line voltage absolute value.
[0048]
[Expression 18]

[0049]
For the system of this reference example, convergence is achieved by 8 iterations. As for the SVR between the nodes 4 and 5, the solution near 1.0 is 1.118884, the solution near 0.0 is 0.000640, and 1.118884 is the solution. It becomes.
[0050]
【The invention's effect】
As described above, in the present invention, the inflow power to each branch line is used as a state variable, and other state variables are sequentially calculated from the top node, and the function of SVR automatically adjusting the tap is formulated. Since the active power and reactive power of the node of interest are obtained using taps obtained from the secondary side voltage specification value, the state variables are greatly reduced, enabling high-speed power flow calculations. Practical power flow calculation considering SVR existing in the system can be performed.
[Brief description of the drawings]
FIG. 1 is an explanatory diagram of a system model.
FIG. 2 is an explanatory diagram of a power transmission line model.
FIG. 3 is an explanatory diagram of a transformer model.
FIG. 4 is an explanatory diagram of a transmission line branch of a two-node system.
FIG. 5 is an explanatory diagram of a transformer branch of a two-node system.
FIG. 6 is an explanatory diagram of the concept of forward calculation and reverse calculation.
FIG. 7 is an explanatory diagram of a two-node system used for three-phase power flow calculation.
FIG. 8 is an explanatory diagram of a single-phase power flow calculation system in a reference example of the present invention.

Claims

In the three- phase power flow calculation method that represents the radial system in which the automatic voltage regulator exists by a network having branches and loads as nodes, and obtains active power, reactive power, and voltage at each node by the state equation,
A step of inputting power supply voltage, load capacity, line data, transformer data, automatic voltage regulator data of the distribution system;
Calculating the sum of loads from the end node of each branch line toward the top node, and calculating the inflow power to the top node as an initial value of the state variable;
When the branch is a transmission line or a transformer, active power, reactive power, and voltage as state quantities of the target node are obtained from outflow power, current, voltage at the upstream node , and load power at the target node, and the branch Is an automatic voltage regulator, the active power, reactive power, and voltage as the state quantities of the node of interest are expressed as a state equation that is a function of the active power, reactive power, and voltage of the upstream node, and the automatic voltage regulator Advance calculation step for each phase using the tap value obtained from the secondary side voltage regulation value,
A backward calculation step of subtracting the outflow power at the end node of each branch line from the inflow power at the leading node to update the state variable;
A convergence determination step for determining that the outflow electric power at the end node of each branch line is equal to or lower than a threshold value and for converging the calculation;
A tidal current calculation method in a radial distribution system characterized by comprising: