JP4112161B2 - Transmission limit search method and recording medium storing transmission limit search program - Google Patents
Transmission limit search method and recording medium storing transmission limit search program Download PDFInfo
- Publication number
- JP4112161B2 JP4112161B2 JP2000227755A JP2000227755A JP4112161B2 JP 4112161 B2 JP4112161 B2 JP 4112161B2 JP 2000227755 A JP2000227755 A JP 2000227755A JP 2000227755 A JP2000227755 A JP 2000227755A JP 4112161 B2 JP4112161 B2 JP 4112161B2
- Authority
- JP
- Japan
- Prior art keywords
- stability
- calculation
- transmission limit
- determination
- limit value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Description
【0001】
【発明の属する技術分野】
この発明は、発電機や負荷などから構成される電力系統の送電限界値を算出するに際して、計算量の軽減や精度の向上を図る送電限界探索方法及び送電限界探索プログラムが記録された記録媒体に関するものである。
【0002】
【従来の技術】
図7は従来の送電限界探索方法を示すフローチャートである。図において、ST1は詳細安定度計算に必要なデータを作成するステップ、ST2はステップST1で作成されたデータを用いて詳細安定度計算を実行するステップ、ST3は詳細安定度計算の計算結果を解析するステップ、ST4は電力系統が安定であるか否かを判定する判定ステップ、ST5は電力系統が安定である場合、潮流量を増やすステップ、ST6は電力系統が不安定である場合、電力系統が不安定になる直前の潮流量を送電限界値と認定するステップである。
【0003】
次に動作について説明する。
送電限界値を算出する場合、まず、詳細安定度計算に必要なデータを作成し(ステップST1)、そのデータを用いて詳細安定度計算を実行する(ステップST2)。
【0004】
そして、その詳細安定度計算の計算結果を解析して(ステップST3)、電力系統が安定であるか否かを判定する(ステップST4)。
電力系統が安定である場合、特定の発電機や負荷の値を増やすなどして潮流量を増やし(ステップST5)、上記の詳細安定度計算を再実行させる。
一方、電力系統が不安定になると、電力系統が不安定になる直前の潮流量を送電限界値と認定する(ステップST6)。
【0005】
【発明が解決しようとする課題】
従来の送電限界探索方法は以上のように構成されているので、詳細安定度計算を繰り返し実行すれば、送電限界値を求めることができるが、詳細安定度計算の計算が完了するまでに長時間を要するため、速やかに送電限界値を求めることができない課題があった。
【0006】
この発明は上記のような課題を解決するためになされたもので、正確な送電限界値を速やかに求めることができる送電限界探索方法及び送電限界探索プログラムが記録された記録媒体を得ることを目的とする。
【0007】
【課題を解決するための手段】
この発明に係る送電限界探索方法は、安定度判別の判別結果が安定であることを示す場合には、送電量を増やして安定度計算と安定度判別の並列実行を継続する一方、その判別結果が不安定であることを示す場合には、安定度計算と安定度判別の実行を中止し、不安定になる直前の送電量を送電限界値と認定するようにしたものである。
【0008】
この発明に係る送電限界探索方法は、安定度計算における事故除去直後の計算結果から大まかな安定度判別を実行して不安定なケースをスクリーニングし、その不安定なケースについて詳細な安定度判別を実行するようにしたものである。
【0009】
この発明に係る送電限界探索方法は、電力系統の運動エネルギーと位置エネルギーを算出し、その運動エネルギーと位置エネルギーを比較して安定度を判別するようにしたものである。
【0010】
この発明に係る送電限界探索方法は、送電限界値が認定されると、安定度計算を実行して、その送電限界値の正当性を確認するようにしたものである。
【0011】
この発明に係る送電限界探索プログラムが記録された記録媒体は、安定度判別の判別結果が安定であることを示す場合には、送電量を増やして安定度計算と安定度判別の実行を継続させる一方、その判別結果が不安定であることを示す場合には、安定度計算と安定度判別の実行を中止させて、不安定になる直前の送電量を送電限界値と認定する送電限界値認定処理手順を設けたものである。
【0012】
この発明に係る送電限界探索プログラムが記録された記録媒体は、安定度判別処理手順が安定度計算における事故除去直後の計算結果から大まかな安定度判別を実行して不安定なケースをスクリーニングし、その不安定なケースについて詳細な安定度判別を実行するようにしたものである。
【0013】
この発明に係る送電限界探索プログラムが記録された記録媒体は、安定度判別処理手順が電力系統の運動エネルギーと位置エネルギーを算出し、その運動エネルギーと位置エネルギーを比較して安定度を判別するようにしたものである。
【0014】
この発明に係る送電限界探索プログラムが記録された記録媒体は、送電限界値が認定されると、安定度計算を実行して、その送電限界値の正当性を確認する確認処理手順を設けたものである。
【0015】
【発明の実施の形態】
以下、この発明の実施の一形態を説明する。
実施の形態1.
図1はこの発明の実施の形態1による送電限界探索方法を示すフローチャートである。図において、ST11は詳細安定度計算に必要なデータを作成するステップ、ST12はステップST11で作成されたデータを用いて詳細安定度計算を実行するステップ、ST13は詳細安定度計算と並列に実行し、その詳細安定度計算の計算途中結果に基づいて短絡電流を算定するステップ、ST14は詳細安定度計算と並列に実行し、その詳細安定度計算の計算途中結果に基づいて過負荷判別を行うステップ、ST15は詳細安定度計算と並列に実行し、その詳細安定度計算の計算途中結果に基づいて電圧安定度判別を行うステップ、ST16は詳細安定度計算と並列に実行し、その詳細安定度計算の計算途中結果に基づいて過渡安定度判別を行うステップである。
【0016】
ST17はステップST13〜ST16における判別結果から電力系統が安定であるか否かを判定するステップ、ST18は電力系統が安定である場合、送電量を増やすステップ、ST19は電力系統が不安定である場合、電力系統が不安定になる直前の送電量を送電限界値と認定するステップである。
【0017】
図2は送電限界探索方法を実行するプログラムの構築例を示す説明図であり、図において、1はステップST13の内容を実行する短絡電流算定部、2はステップST14の内容を実行する過負荷判別部、3はステップST15の内容を実行する電圧判別部、4はステップST16の内容を実行する過渡安定度判別部である。
なお、送電限界探索方法を実行するプログラム(安定度計算処理手順、安定度判別処理手順及び送電限界値認定処理手順から構成される送電限界探索プログラム)は、計算機が読み取ることが可能なハードディスクやCDROMなどの記録媒体に記録されるものとする。
【0018】
次に動作について説明する。
送電限界値を算出する場合、まず、詳細安定度計算に必要なデータを作成し(ステップST11)、そのデータを用いて詳細安定度計算を実行する(ステップST12)。
【0019】
そして、詳細安定度計算の実行を開始すると、その詳細安定度計算と並列に安定度判別を実行する(ステップST13〜ST16)。
即ち、その詳細安定度計算の計算途中結果に基づく短絡電流の算定と、その詳細安定度計算の計算途中結果に基づく過負荷の判定と、その詳細安定度計算の計算途中結果に基づく電圧安定度の判別と、その詳細安定度計算の計算途中結果に基づく過渡安定度の判別とを行う。
【0020】
ステップST13〜ST16における判別結果を参照し、電力系統が安定であるか否かを判定する(ステップST17)。
電力系統が安定である場合、特定の発電機や負荷の値を増やすなどして送電量を増やし(ステップST18)、上記の安定度判別を再実行させる。
一方、電力系統が不安定になると、詳細安定度計算と安定度判別の実行を中止し、電力系統が不安定になる直前の送電量を送電限界値と認定する(ステップST19)。
【0021】
以上で明らかなように、この実施の形態1によれば、安定度判別の判別結果が安定であることを示す場合には、送電量を増やして安定度計算と安定度判別の並列実行を継続する一方、その判別結果が不安定であることを示す場合には、安定度計算と安定度判別の実行を中止し、不安定になる直前の送電量を送電限界値と認定するように構成したので、詳細安定度計算の計算が完了する前に、安定度の判別結果を得て送電限界値を求めることができる。したがって、正確な送電限界値を速やかに求めることができる効果を奏する。
【0022】
実施の形態2.
図3は想定される事故ケースや潮流断面から明らかに安定と判別されるケースを効率的に絞り込むスクリーニング処理を示すフローチャートであり、図4はスクリーニングの結果、不安定であると判定されたケースに対して行う詳細な安定度判別処理を示すフローチャートである。
【0023】
上記実施の形態1では、詳細安定度計算と安定度判別処理を並列に実行するものについて示したが、安定度判別処理を実行するに際して、最初にスクリーニング処理を実行してから詳細な安定度判別処理を実行するようにしてもよい。
【0024】
最初にスクリーニング処理の内容から説明する。
スクリーニング処理とは、図5に示すような対象系統内の発電機を加速側と減速側の2グループに分割して、各グループを1台の等価発電機に縮約し、等面積法を適用して過渡安定度判別を行うものである。
【0025】
まず、事故除去後の詳細安定度結果から、系統の慣性中心δcenterを求める。
【数1】
次に、各発電機の位相角δi と慣性中心δcenterとの大小比較を実行し、慣性中心δcenterより位相角δi が大きい発電機群を加速側グループ、慣性中心δcenterより位相角δi が小さい発電機群を減速側グループとして、これらの発電機群を(2)式〜(6)式に示すように、等価発電機に集約する。
【数2】
【0026】
次に、上記の2機系を安定度的に等価な1機系に縮約する。
【数3】
(9)式の機械的出力PE は次のように表すことができる。
PE =P0 +P1 ・sinδ+P2 ・cosδ (12)
ここで、P0 =0と仮定すると、(12)式は次のように変形することができる。
PE =P1 ・sinδ+P2 ・cosδ=PPEAKsin(δ+δ0 )(13)
ただし、PPEAK=(P1 2+P2 2)1/2
【0027】
次に、等価1機発電機の機械的出力PE と位相角δから事故後の電力相差角曲線:PE =PPEAKsin(δ+δ0 )を求める。ただし、初期位相角δ0 はわからないため、δ0 =0と仮定する。
そして、事故後の電力相差角曲線より、運動エネルギーVk を下記に示すように算出する。
Vk =Mω2 (t= Δ T) /2ω0 (14)
【0028】
δ0 =0と仮定しているため、電力相差角曲線はPE =PPEAKsinδとなり、PPEAK及び不安定平衡点δu は下記のようにして算出することができる。
【数4】
そして、(15)式及び(16)式を下式に代入して、位置エネルギーである臨界エネルギーVc を算出する。
【数5】
【0029】
上記のようにして、運動エネルギーVk と臨界エネルギーVc を算出すると、これらを大小比較することにより、安定であるか不安定であるかを判別し、スクリーニング処理を終了する。
Vk <Vc :安定
Vk ≧Vc :不安定
【0030】
次に詳細な安定度判別処理の内容について説明する。
上記のスクリーニング処理において、不安定であると判定されたケースを対象にして、詳細な安定度判別処理を一定間隔毎に行う。
【0031】
まず、スクリーニング処理と同様に、図5に示すような対象系統内の発電機を加速側と減速側の2グループに分割して、各グループを1台の等価発電機に縮約し、等面積法を適用して過渡安定度判別を行う。
詳細過渡安定度判別を行う時点(時刻t)の詳細安定度計算結果から系統の慣性中心δcenterを求める。
【数6】
【0032】
次に、各発電機の位相角δi (t)と慣性中心δcenterとの大小比較を実行し、慣性中心δcenterより位相角δi (t)が大きい発電機群を加速側グループ、慣性中心δcenterより位相角δi (t)が小さい発電機群を減速側グループとして、これらの発電機群を(19)式〜(23)式に示すように、等価発電機に集約する。
【数7】
【0033】
次に、上記の2機系を安定度的に等価な1機系に縮約する。
【数8】
等価1機発電機の電力相差角曲線はPE =PPEAKsin(δ+δ0 )で表されるが、これを変換すると次のようになる。
【数9】
【0034】
次に、詳細安定度計算のサンプリングデータから最小2乗法を用いて、(29)式の未知の項P1 ,P2 を算出する。
なお、サンプリングデータは、事故除去の過渡的な値を避けるため、事故除去から一定時間後のデータを最低5点、最大30点用いる。
【数10】
(30)式より、未知のパラメータX(P1 ,P2 )は、最小2乗法を用いて、以下の式より推定することができる。
X=(AT A)-1AT b (31)
【0035】
(31)式より算出した未知の項P1 ,P2 より、PPEAK及び不安定平衡点δu を算出する。
【数11】
運動エネルギーVk と臨界エネルギーVc は下記のようにして算出する。
【数12】
【0036】
上記のようにして、運動エネルギーVk と臨界エネルギーVc を算出すると、これらを大小比較することにより、安定であるか不安定であるかを判別し、詳細な安定度判別処理を終了する。
Vk <Vc :安定
Vk ≧Vc :不安定
【0037】
以上で明らかなように、この実施の形態2によれば、安定度計算における事故除去直後の計算結果から大まかな安定度判別を実行して不安定なケースをスクリーニングし、その不安定なケースについて詳細な安定度判別を実行するように構成したので、すべてのケースについて詳細な安定度判別を実行する必要がなくなる。したがって、上記実施の形態1よりも更に速やかに送電限界値を求めることができる効果を奏する。
【0038】
実施の形態3.
上記実施の形態1,2では、安定度判別を実行して送電限界値を求めるものについて示したが、図6に示すように、さらに、安定度計算を実行して、その送電限界値の正当性を確認するようにしてもよい。
【0039】
具体的には、上記実施の形態1,2と同様にして、暫定的に送電限界値を算出すると(ステップST19)、その送電限界値で送電する条件の下で、詳細な安定度計算を実行し(ステップST20)、確認計算を実行する(ステップST21)。
【0040】
そして、当該条件下での計算結果が安定であることを示し、かつ、当該送電限界値を1ステップ増加して、確認計算をしたとき不安定であることを示す場合には、その送電限界値が正しいものと認定する(ステップST22〜ST24)。
一方、当該条件下での計算結果が不安定であることを示す場合、あるいは、当該条件下での計算結果が安定であることを示し、かつ、当該送電限界値を1ステップ増加して、確認計算をしたとき安定であることを示す場合には、当該送電限界値が正しくないものと認定し、適宜送電量を変化させて正しい送電限界値を算出する(ステップST22〜ST24)。
これにより、この実施の形態3によれば、精度の高い送電限界値を得ることができる効果を奏する。
【0041】
【発明の効果】
以上のように、この発明によれば、安定度判別の判別結果が安定であることを示す場合には、送電量を増やして安定度計算と安定度判別の並列実行を継続する一方、その判別結果が不安定であることを示す場合には、安定度計算と安定度判別の実行を中止し、不安定になる直前の送電量を送電限界値と認定するように構成したので、正確な送電限界値を速やかに求めることができる効果がある。
【0042】
この発明によれば、安定度計算における事故除去直後の計算結果から大まかな安定度判別を実行して不安定なケースをスクリーニングし、その不安定なケースについて詳細な安定度判別を実行するように構成したので、更に速やかに送電限界値を求めることができる効果がある。
【0043】
この発明によれば、電力系統の運動エネルギーと位置エネルギーを算出し、その運動エネルギーと位置エネルギーを比較して安定度を判別するように構成したので、処理の複雑化を招くことなく、正確な安定度判別を行うことができる効果がある。
【0044】
この発明によれば、送電限界値が認定されると、安定度計算を実行して、その送電限界値の正当性を確認するように構成したので、精度の高い送電限界値を得ることができる効果がある。
【0045】
この発明によれば、安定度判別の判別結果が安定であることを示す場合には、送電量を増やして安定度計算と安定度判別の実行を継続させる一方、その判別結果が不安定であることを示す場合には、安定度計算と安定度判別の実行を中止させて、不安定になる直前の送電量を送電限界値と認定する送電限界値認定処理手順を設けるように構成したので、正確な送電限界値を速やかに求めることができる効果がある。
【0046】
この発明によれば、安定度判別処理手順が安定度計算における事故除去直後の計算結果から大まかな安定度判別を実行して不安定なケースをスクリーニングし、その不安定なケースについて詳細な安定度判別を実行するように構成したので、更に速やかに送電限界値を求めることができる効果がある。
【0047】
この発明によれば、安定度判別処理手順が電力系統の運動エネルギーと位置エネルギーを算出し、その運動エネルギーと位置エネルギーを比較して安定度を判別するように構成したので、処理の複雑化を招くことなく、正確な安定度判別を行うことができる効果がある。
【0048】
この発明によれば、送電限界値が認定されると、安定度計算を実行して、その送電限界値の正当性を確認する確認処理手順を設けるように構成したので、精度の高い送電限界値を得ることができる効果がある。
【図面の簡単な説明】
【図1】 この発明の実施の形態1による送電限界探索方法を示すフローチャートである。
【図2】 送電限界探索方法を実行するプログラムの構築例を示す説明図である。
【図3】 想定される事故ケースや潮流断面から明らかに安定と判別されるケースを効率的に絞り込むスクリーニング処理を示すフローチャートである。
【図4】 スクリーニングの結果、不安定であると判定されたケースに対して行う詳細な安定度判別処理を示すフローチャートである。
【図5】 等価発電機の縮約を説明する説明図である。
【図6】 この発明の実施の形態3による送電限界探索方法を示すフローチャートである。
【図7】 従来の送電限界探索方法を示すフローチャートである。
【符号の説明】
1 短絡電流算定部、2 過負荷判別部、3 電圧判別部、4 過渡安定度判別部。[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a transmission limit search method for reducing a calculation amount and improving accuracy when calculating a transmission limit value of a power system including a generator, a load, and the like, and a recording medium on which a transmission limit search program is recorded. Is.
[0002]
[Prior art]
FIG. 7 is a flowchart showing a conventional power transmission limit search method. In the figure, ST1 is a step of creating data necessary for detailed stability calculation, ST2 is a step of executing detailed stability calculation using the data created in step ST1, and ST3 is an analysis of a calculation result of the detailed stability calculation. ST4 is a determination step for determining whether or not the power system is stable, ST5 is a step for increasing the tidal flow when the power system is stable, and ST6 is a step when the power system is unstable when the power system is unstable. In this step, the tidal flow just before the instability is recognized as the transmission limit value.
[0003]
Next, the operation will be described.
When calculating the power transmission limit value, first, data necessary for the detailed stability calculation is created (step ST1), and the detailed stability calculation is executed using the data (step ST2).
[0004]
Then, the calculation result of the detailed stability calculation is analyzed (step ST3), and it is determined whether or not the power system is stable (step ST4).
When the power system is stable, the flow rate is increased by increasing the value of a specific generator or load (step ST5), and the above detailed stability calculation is re-executed.
On the other hand, when the power system becomes unstable, the tidal flow immediately before the power system becomes unstable is recognized as the transmission limit value (step ST6).
[0005]
[Problems to be solved by the invention]
Since the conventional transmission limit search method is configured as described above, the transmission limit value can be obtained by repeatedly executing the detailed stability calculation, but it takes a long time to complete the calculation of the detailed stability calculation. Therefore, there is a problem that the transmission limit value cannot be obtained promptly.
[0006]
The present invention has been made to solve the above-described problems, and an object of the present invention is to obtain a transmission limit search method and a recording medium on which a transmission limit search program is recorded, in which an accurate transmission limit value can be quickly obtained. And
[0007]
[Means for Solving the Problems]
In the transmission limit search method according to the present invention, when the determination result of the stability determination indicates that the determination is stable, the power transmission amount is increased and the parallel calculation of the stability calculation and the stability determination is continued while the determination result When it indicates that is unstable, the calculation of stability and the determination of stability are stopped, and the transmission amount immediately before becoming unstable is recognized as the transmission limit value.
[0008]
The transmission limit search method according to the present invention performs a rough stability determination from the calculation result immediately after the accident removal in the stability calculation, screens an unstable case, and performs a detailed stability determination for the unstable case. It is something to be executed.
[0009]
The power transmission limit search method according to the present invention calculates kinetic energy and potential energy of an electric power system, and compares the kinetic energy and potential energy to determine stability.
[0010]
In the transmission limit search method according to the present invention, when the transmission limit value is recognized, the stability calculation is executed and the validity of the transmission limit value is confirmed.
[0011]
In the recording medium on which the transmission limit search program according to the present invention is recorded, when the determination result of the stability determination indicates that it is stable, the power transmission amount is increased and the execution of the stability calculation and the stability determination is continued. On the other hand, if the judgment result shows that it is unstable, the execution of stability calculation and stability judgment is stopped, and the power transmission limit value authorization that authorizes the power transmission amount immediately before becoming unstable as the power transmission limit value. A processing procedure is provided.
[0012]
The recording medium in which the transmission limit search program according to the present invention is recorded, the stability determination processing procedure performs rough stability determination from the calculation result immediately after the accident removal in the stability calculation, and screens an unstable case, Detailed stability determination is performed for the unstable case.
[0013]
In the recording medium recorded with the transmission limit search program according to the present invention, the stability determination processing procedure calculates the kinetic energy and potential energy of the power system, and compares the kinetic energy and potential energy to determine the stability. It is a thing.
[0014]
The recording medium in which the transmission limit search program according to the present invention is recorded is provided with a confirmation processing procedure for confirming the validity of the transmission limit value by executing the stability calculation when the transmission limit value is recognized. It is.
[0015]
DETAILED DESCRIPTION OF THE INVENTION
An embodiment of the present invention will be described below.
FIG. 1 is a flowchart showing a transmission limit search method according to
[0016]
ST17 is a step for determining whether or not the power system is stable from the determination results in steps ST13 to ST16, ST18 is a step for increasing the transmission amount when the power system is stable, and ST19 is a case where the power system is unstable. In this step, the transmission amount immediately before the power system becomes unstable is recognized as the transmission limit value.
[0017]
FIG. 2 is an explanatory diagram showing a construction example of a program for executing the transmission limit search method. In the figure, 1 is a short-circuit current calculation unit that executes the contents of step ST13, and 2 is an overload determination that executes the contents of step ST14. 3 and 3 are a voltage determination unit that executes the content of step ST15, and 4 is a transient stability determination unit that executes the content of step ST16.
Note that a program for executing a transmission limit search method (a transmission limit search program including a stability calculation processing procedure, a stability determination processing procedure, and a transmission limit value recognition processing procedure) is a hard disk or CDROM that can be read by a computer. It is recorded on a recording medium such as
[0018]
Next, the operation will be described.
When calculating the power transmission limit value, first, data necessary for the detailed stability calculation is created (step ST11), and the detailed stability calculation is executed using the data (step ST12).
[0019]
When the execution of the detailed stability calculation is started, the stability determination is executed in parallel with the detailed stability calculation (steps ST13 to ST16).
That is, calculation of short-circuit current based on the calculation result of the detailed stability calculation, determination of overload based on the calculation result of the detailed stability calculation, and voltage stability based on the calculation result of the detailed stability calculation And transient stability based on the result of the calculation of the detailed stability.
[0020]
With reference to the determination results in steps ST13 to ST16, it is determined whether or not the power system is stable (step ST17).
If the power system is stable, the power transmission amount is increased by increasing the value of a specific generator or load (step ST18), and the above stability determination is re-executed.
On the other hand, when the power system becomes unstable, the detailed stability calculation and the stability determination are stopped, and the transmission amount immediately before the power system becomes unstable is recognized as the transmission limit value (step ST19).
[0021]
As is apparent from the above, according to the first embodiment, when the determination result of the stability determination indicates that it is stable, the power transmission amount is increased and the parallel calculation of the stability calculation and the stability determination is continued. On the other hand, when the determination result indicates that it is unstable, the stability calculation and stability determination are stopped, and the transmission amount immediately before becoming unstable is recognized as the transmission limit value. Therefore, before the calculation of the detailed stability calculation is completed, it is possible to obtain the transmission limit value by obtaining the stability determination result. Therefore, there is an effect that an accurate power transmission limit value can be quickly obtained.
[0022]
FIG. 3 is a flowchart showing a screening process for efficiently narrowing down an assumed accident case or a case that is clearly determined to be stable from a tidal current section, and FIG. 4 shows a case where the result of screening is determined to be unstable. It is a flowchart which shows the detailed stability determination process performed with respect to it.
[0023]
In the first embodiment, the detailed stability calculation and the stability determination process are performed in parallel. However, when the stability determination process is performed, the screening process is first performed and then the detailed stability determination is performed. Processing may be executed.
[0024]
First, the contents of the screening process will be described.
The screening process is to divide the generator in the target system as shown in Fig. 5 into two groups, the acceleration side and the deceleration side, and reduce each group to one equivalent generator and apply the equal area method. Thus, the transient stability determination is performed.
[0025]
First, the inertia center δ center of the system is obtained from the detailed stability result after the accident removal.
[Expression 1]
Next, the phase angle δ i of each generator And the inertia center δ center are compared, and the phase angle δ i from the inertia center δ center is calculated. Generator group with a large value is the acceleration side group, the phase angle δ i from the center of inertia δ center A generator group having a small value is set as a deceleration side group, and these generator groups are aggregated into equivalent generators as shown in equations (2) to (6).
[Expression 2]
[0026]
Next, the above-mentioned two-machine system is reduced to a one-machine system that is equivalent in terms of stability.
[Equation 3]
The mechanical output P E in equation (9) can be expressed as follows.
P E = P 0 + P 1 · sin δ + P 2 · cos δ (12)
Here, assuming that P 0 = 0, equation (12) can be modified as follows.
P E = P 1 · sin δ + P 2 · cos δ = P PEAK sin (δ + δ 0 ) (13)
However, P PEAK = (P 1 2 + P 2 2 ) 1/2
[0027]
Next, a post-accident power phase difference curve: P E = P PEAK sin (δ + δ 0 ) is obtained from the mechanical output P E of the equivalent single generator and the phase angle δ. However, since the initial phase angle δ 0 is unknown, it is assumed that δ 0 = 0.
Then, the kinetic energy V k is calculated as shown below from the power phase difference angle curve after the accident.
V k = Mω 2 (t = ΔT ) / 2ω 0 (14)
[0028]
Since it is assumed that δ 0 = 0, the power phase difference angle curve is P E = P PEAK sin δ, and P PEAK and the unstable equilibrium point δ u can be calculated as follows.
[Expression 4]
Then, the critical energy V c that is potential energy is calculated by substituting the equations (15) and (16) into the following equation.
[Equation 5]
[0029]
When the kinetic energy V k and the critical energy V c are calculated as described above, they are compared to determine whether they are stable or unstable, and the screening process is terminated.
V k <V c : Stable V k ≧ V c : Unstable
Next, the details of the stability determination process will be described.
In the screening process described above, a detailed stability determination process is performed at regular intervals for cases determined to be unstable.
[0031]
First, as in the screening process, the generator in the target system as shown in FIG. 5 is divided into two groups, the acceleration side and the deceleration side, and each group is reduced to one equivalent generator, and the same area is obtained. Apply the method to determine the transient stability.
The inertia center δ center of the system is obtained from the detailed stability calculation result at the time (time t) when the detailed transient stability determination is performed.
[Formula 6]
[0032]
Next, the phase angle δ i of each generator (T) and the inertia center δ center are compared, and the phase angle δ i is calculated from the inertia center δ center. The generator group with a large (t) is the acceleration side group, the phase angle δ i from the inertia center δ center A generator group having a small (t) is set as a deceleration side group, and these generator groups are aggregated into equivalent generators as shown in equations (19) to (23).
[Expression 7]
[0033]
Next, the above-mentioned two-machine system is reduced to a one-machine system that is equivalent in terms of stability.
[Equation 8]
The power phase difference angle curve of the equivalent single-machine generator is expressed as P E = P PEAK sin (δ + δ 0 ), and this is converted as follows.
[Equation 9]
[0034]
Next, the unknown terms P 1 and P 2 of the equation (29) are calculated from the sampling data of the detailed stability calculation using the least square method.
The sampling data uses a minimum of 5 points and a maximum of 30 points after a certain period of time since the accident is removed in order to avoid a transient value for the accident removal.
[Expression 10]
From the equation (30), the unknown parameter X (P 1 , P 2 ) can be estimated from the following equation using the least square method.
X = (A T A) −1 A T b (31)
[0035]
P PEAK and unstable equilibrium point δ u are calculated from unknown terms P 1 and P 2 calculated from equation (31).
[Expression 11]
The kinetic energy V k and the critical energy V c are calculated as follows.
[Expression 12]
[0036]
When the kinetic energy V k and the critical energy V c are calculated as described above, the magnitude is compared to determine whether the energy is stable or unstable, and the detailed stability determination processing is terminated.
V k <V c : Stable V k ≧ V c : Unstable
As apparent from the above, according to the second embodiment, an unstable case is screened by performing rough stability determination from the calculation result immediately after the accident removal in the stability calculation. Since detailed stability determination is performed, it is not necessary to perform detailed stability determination for all cases. Therefore, the power transmission limit value can be obtained more quickly than in the first embodiment.
[0038]
In the first and second embodiments, the stability determination is performed to determine the transmission limit value. However, as shown in FIG. 6, the stability calculation is further performed to verify the transmission limit value. You may make it confirm sex.
[0039]
Specifically, in the same manner as in the first and second embodiments, when a transmission limit value is tentatively calculated (step ST19), detailed stability calculation is performed under the condition of transmitting power at the transmission limit value. (Step ST20) and confirmation calculation is executed (Step ST21).
[0040]
If the calculation result under the condition is stable and the transmission limit value is increased by one step to indicate that the calculation result is unstable when the confirmation calculation is performed, the transmission limit value Is recognized as correct (steps ST22 to ST24).
On the other hand, if the calculation result under the condition is unstable, or the calculation result under the condition is stable and the transmission limit value is increased by one step When it is shown that the calculation is stable, it is recognized that the transmission limit value is not correct, and a correct transmission limit value is calculated by appropriately changing the transmission amount (steps ST22 to ST24).
Thereby, according to this
[0041]
【The invention's effect】
As described above, according to the present invention, when the determination result of the stability determination indicates that the determination is stable, the power transmission amount is increased and the parallel calculation of the stability calculation and the stability determination is continued, while the determination If the result shows that the result is unstable, the stability calculation and stability determination are stopped, and the transmission amount immediately before becoming unstable is recognized as the transmission limit value. There is an effect that the limit value can be obtained quickly.
[0042]
According to the present invention, an unstable case is screened by executing a rough stability determination from the calculation result immediately after the accident removal in the stability calculation, and a detailed stability determination is executed for the unstable case. Since it comprised, there exists an effect which can obtain | require a power transmission limit value more rapidly.
[0043]
According to the present invention, the kinetic energy and potential energy of the power system are calculated, and the kinetic energy and potential energy are compared to determine the stability. There is an effect that the stability determination can be performed.
[0044]
According to the present invention, when the power transmission limit value is certified, the stability calculation is executed and the validity of the power transmission limit value is confirmed, so that a highly accurate power transmission limit value can be obtained. effective.
[0045]
According to the present invention, when the determination result of the stability determination indicates that it is stable, the power transmission amount is increased to continue the calculation of the stability and the determination of the stability, while the determination result is unstable. In order to show that, because it was configured to stop the calculation of stability and stability determination, and to provide a transmission limit value certification processing procedure to authorize the transmission amount immediately before becoming unstable as the transmission limit value, There is an effect that an accurate transmission limit value can be quickly obtained.
[0046]
According to the present invention, the stability determination processing procedure performs rough stability determination from the calculation result immediately after the accident removal in the stability calculation to screen the unstable case, and the detailed stability for the unstable case is screened. Since the determination is performed, the power transmission limit value can be obtained more quickly.
[0047]
According to the present invention, the stability determination processing procedure calculates the kinetic energy and potential energy of the power system, and compares the kinetic energy and potential energy to determine the stability. There is an effect that accurate stability determination can be performed without incurring.
[0048]
According to the present invention, when the power transmission limit value is certified, the stability calculation is executed and the confirmation processing procedure for confirming the validity of the power transmission limit value is provided. There is an effect that can be obtained.
[Brief description of the drawings]
FIG. 1 is a flowchart showing a power transmission limit search method according to
FIG. 2 is an explanatory diagram showing a construction example of a program for executing a power transmission limit search method.
FIG. 3 is a flowchart showing a screening process for efficiently narrowing down an assumed accident case or a case that is clearly determined to be stable from a tidal current cross section.
FIG. 4 is a flowchart showing detailed stability determination processing performed for a case determined to be unstable as a result of screening;
FIG. 5 is an explanatory diagram for explaining contraction of an equivalent generator.
FIG. 6 is a flowchart showing a power transmission limit search method according to
FIG. 7 is a flowchart showing a conventional power transmission limit search method.
[Explanation of symbols]
1 short-circuit current calculation unit, 2 overload determination unit, 3 voltage determination unit, 4 transient stability determination unit.
Claims (8)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000227755A JP4112161B2 (en) | 2000-07-27 | 2000-07-27 | Transmission limit search method and recording medium storing transmission limit search program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000227755A JP4112161B2 (en) | 2000-07-27 | 2000-07-27 | Transmission limit search method and recording medium storing transmission limit search program |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2002051469A JP2002051469A (en) | 2002-02-15 |
JP4112161B2 true JP4112161B2 (en) | 2008-07-02 |
Family
ID=18721160
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2000227755A Expired - Lifetime JP4112161B2 (en) | 2000-07-27 | 2000-07-27 | Transmission limit search method and recording medium storing transmission limit search program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4112161B2 (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007189840A (en) * | 2006-01-13 | 2007-07-26 | Toshiba Corp | Power system stabilizing apparatus |
JP5986827B2 (en) * | 2012-07-03 | 2016-09-06 | 一般財団法人電力中央研究所 | Power system stabilization analysis device, power system stabilization analysis method, and power system stabilization analysis program |
WO2020095401A1 (en) * | 2018-11-08 | 2020-05-14 | 株式会社日立製作所 | Available power transmission capacity analysis device, available power transmission capacity analysis method, and computer program |
-
2000
- 2000-07-27 JP JP2000227755A patent/JP4112161B2/en not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
JP2002051469A (en) | 2002-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Castanier et al. | Characteristic constraint modes for component mode synthesis | |
Feiner et al. | Mistuning identification of bladed disks using a fundamental mistuning model—part II: application | |
Salvat et al. | Modeling of abradable coating removal in aircraft engines through delay differential equations | |
Hackenberg et al. | An approach for estimating the effect of transient sweep through a resonance | |
Kurstak et al. | Parametric reduced order models for bladed disks with mistuning and varying operational speed | |
Saito et al. | Effects of a cracked blade on mistuned turbine engine rotor vibration | |
Holland et al. | Testing and calibration procedures for mistuning identification and traveling wave excitation of blisks | |
JP2019512774A (en) | Data processing to characterize battery thermal behavior | |
US6636817B2 (en) | Method and apparatus for enhanced mechanical signature analysis | |
Kurstak et al. | Multistage blisk and large mistuning modeling using fourier constraint modes and PRIME | |
Pesaresi et al. | Numerical and experimental investigation of an underplatform damper test rig | |
Zhu et al. | Velocity estimation of motion systems based on low-resolution encoders | |
Kurstak et al. | Parametric reduced order models for bladed disks with mistuning and varying operational speed | |
JP4112161B2 (en) | Transmission limit search method and recording medium storing transmission limit search program | |
Buchwald et al. | On the influence of blade aspect ratio on aerodynamic damping | |
JP4505392B2 (en) | Critical failure removal time calculation method, critical failure removal time calculation program, and recording medium | |
CN112414697A (en) | Wind power blade degree of freedom fatigue loading control method, device and system | |
Schwerdt et al. | A model reduction method for bladed disks with large geometric mistuning using a partially reduced intermediate system model | |
Jha et al. | Spatial discretization of axially moving media vibration problems | |
Rivas-Guerra et al. | Identification of mistuning characteristics of bladed disks from free response data—part ii | |
CN106940735B (en) | Method for predicting elastic vibration stability of induction motor | |
Kielb et al. | Flutter of low pressure turbine blades with cyclic symmetric modes: A preliminary design method | |
Mei et al. | Active control of regenerative chatter during metal cutting process | |
Chong et al. | Practical and robust control for precision motion: AR‐CM NCTF control of a linear motion mechanism with friction characteristics | |
JP3865608B2 (en) | Power transmission limit calculation method for power system and computer-readable recording medium recorded with transmission limit calculation program |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20050803 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20050803 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20061002 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20080222 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20080311 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080409 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110418 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120418 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120418 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130418 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130418 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140418 Year of fee payment: 6 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |