KR100840101B1

KR100840101B1 - Protection method of fault propagation using real time transient stability criterion in power systems

Info

Publication number: KR100840101B1
Application number: KR1020060092635A
Authority: KR
Inventors: 박은희; 심관식; 남해곤
Original assignee: 박은희; 심관식
Priority date: 2006-09-25
Filing date: 2006-09-25
Publication date: 2008-06-19
Also published as: KR20080027517A

Abstract

본 발명은 전력시스템에서 실시간으로 과도안정도를 판별하고 고장이 파급되는 것을 방지하는 방법에 관한 것으로, 실시간으로 취득한 이산데이터에서 전력을 계산하고, 이로부터 고장을 판단한다. 고장이 발생하였을 때, 실측한 기계적 입력과 전력의 편차로부터 가속력을 계산하고, 가속력으로부터 발전기에 누적되는 가속에너지를 계산한다. 가속에너지가 미리 설정된 설정값의 범위를 벗어나면 과도불안정으로 판단하여 발전기를 탈락함으로써 다른 발전기들은 정상운전을 가능케 한다. 지금까지 과도안정도 판별법에서는 전력위상각을 필요로 하므로 실시간 과도안정도 해석이 불가능하고, 시스템을 축약한 후 등면적법을 적용하여 과도안정도를 해석하였다. 그러나 본 발명에서는 위상각을 사용하지 않고 실시간으로 실측한 전력으로부터 데이터가 입력된 사이클 마다 과도안정도를 판별할 수 있어, 다른 발전기로 심각한 고장이 파급되는 것을 방지할 수 있다. 따라서 보다 안정적으로 전력을 생산, 공급할 수 있고, 불필요한 발전기 정지를 예방할 수 있어 막대한 비용절감도 기대할 수 있다.The present invention relates to a method for determining the transient stability in real time in the power system and preventing the failure from spreading. The present invention calculates power from discrete data acquired in real time and determines a failure therefrom. When a failure occurs, the acceleration force is calculated from the deviation of the measured mechanical input and power, and the acceleration energy accumulated in the generator is calculated from the acceleration force. If the acceleration energy is out of the preset setting range, it is judged as transient instability and the generator is dropped so that other generators can operate normally. Until now, the transient stability method requires power phase angle, so it is impossible to analyze the transient stability in real time. After reducing the system, the transient stability method is applied to analyze the transient stability. However, in the present invention, the transient stability can be determined for each cycle in which data is input from power measured in real time without using a phase angle, thereby preventing the occurrence of a serious failure with another generator. Therefore, it is possible to produce and supply electricity more stably, and to prevent unnecessary generator stoppages, thus enormous cost reduction can be expected.

Description

전력시스템에서 실시간 과도안정도 판별 방법을 이용한 고장파급방지법{ Protection method of fault propagation using real time transient stability criterion in power systems} Protection method of fault propagation using real time transient stability criterion in power systems

도 1은 본 발명을 수행하기 위한 하드웨어의 개략적인 블록 구성도,1 is a schematic block diagram of hardware for carrying out the present invention;

도 2는 등면적법을 설명하기 위한 전력-위상각 곡선,2 is a power-phase angle curve for explaining an isometric area method,

도 3은 전기적 출력과 기계적 입력 및 가속력을 나타내는 계통도,3 is a schematic diagram showing electrical output, mechanical input and acceleration force;

도 4는 과도고장시 가속에너지를 나타내는 도면,4 is a diagram showing acceleration energy during transient failure,

도 5는 전기적 에너지와 기계적 에너지 및 가속에너지를 나타내는 계통도,5 is a system diagram showing electrical energy, mechanical energy and acceleration energy;

도 6은 과도고장시 가속에너지와 전기적 에너지를 나타내는 도면,6 is a view showing acceleration energy and electrical energy during transient failure,

도 7은 과도상태에서 가속에너지를 나타내는 도면,7 is a view showing acceleration energy in a transient state,

도 8은 시간에 대한 가속에너지 도면,8 is an acceleration energy plot with respect to time,

도 9는 실시예를 위한 간단한 계통의 단선도,9 is a schematic diagram of a simple system for an embodiment,

도 10은 실시예를 위한 도9의 발전기들의 가속에너지,10 is the acceleration energy of the generators of FIG. 9 for an embodiment,

도 11은 도9의 발전기 1의 전력에 대한 시간응답 파형,FIG. 11 is a time response waveform for power of generator 1 of FIG. 9;

도 12는 도9의 발전기 5의 전력에 대한 시간응답 파형,12 is a time response waveform for the power of generator 5 of FIG.

도 13은 실시예를 위한 실측 전력 파형.13 is a measured power waveform for an embodiment.

<도면의 주요 부분에 대한 부호의 설명><Explanation of symbols for the main parts of the drawings>

10: 데이터 취득부 20: 입력 처리부10: data acquisition unit 20: input processing unit

30: 설정부 40: 연산부30: setting unit 40: calculating unit

60: 제어부 70: 출력부60: control unit 70: output unit

전력시스템은 발전기를 비롯한 다양한 기기와 제어기들이 네트워크를 통하여 서로 조합되어 운용되고 있다. 만일 발전기 모선이나 선로에 심각한 고장이 발생하였을 때, 특정 지역의 발전기들은 과도불안정 상태에 놓이게 되고 지역 전체 발전기가 정지될 수 있다. 전체 발전기의 탈락을 예방하기위하여 일부 발전기를 탈락시켜 나머지 발전기가 정지되는 것을 예방하고 있는데, 이러한 시스템을 고장파급방지장치라 한다. 고장파급방지장치로 발전기들을 안정시켜 계통영향을 최소화하도록 하고 있으며, 이러한 대책에도 불구하고 대단위 발전단지가 동시 정지할 경우 발전기 정지로 인한 시스템 주파수 저하정도에 따라 부하를 동시에 자동차단하여, 수급 균형을 유지함으로써 계통을 조기에 안정화시켜 광역정전을 방지하도록 하고 있다.In the power system, various devices and controllers including generators are operated in combination with each other through a network. If a serious failure of the generator bus or track occurs, the generators in a given area may be in an unstable state and the entire generator may be stopped. In order to prevent the fall of the entire generator, some generators are dropped to prevent the rest of the generators from being stopped. Such a system is called a fault spreading device. In order to minimize the system impact by stabilizing the generators with the fault spread prevention device, and in spite of such measures, if the large power generation complex stops simultaneously, the load is automatically cut off according to the system frequency drop caused by the generator stop, By maintaining it, the system is stabilized early to prevent widespread outages.

전력시스템에서 과도안정도는 단락사고나 발전력 또는 연계선의 갑작스런 상실 등 가혹한 외란이나 스위칭 조작에 의해서 발생되는 충격이 계통에 갑자기 가해질 때, 모든 발전기가 동기 상태를 유지할 수 있는지 여부가 관심이다. 외란의 크기와 돌발성 때문에 과도안정도 해석은 고장 발생으로부터 수초 간에 초점을 두며 짧은 해석시간으로 초기 동요(first swing) 해석이라 한다. 보통 과도안정도는 초기에 동요의 진폭이 감쇄되어 시스템이 안정하게 운전되지만 자동주파수 및 자동전압제어루프와 같은 제어 계통의 영향으로 수 분간 동요가 지속될 수도 있다.Transient stability in power systems is concerned about whether all generators can remain synchronized when the system is suddenly impacted by severe disturbances or switching operations, such as short-circuit accidents or sudden loss of generating power or linkage. Because of the magnitude and sudden nature of the disturbance, transient stability analysis focuses on the number of seconds from failure and is called the first swing analysis with a short analysis time. Normally, transient stability is initially reduced by the amplitude of the fluctuations, so that the system can operate stably, but fluctuations may last for several minutes under the influence of control systems such as automatic frequency and automatic voltage control loops.

지금까지 이러한 과도안정도는 일기무한모선 계통의 동요방정식과 등면적법(equal-area criterion)으로 해석되고, 대규모 계통에서도 계통축약을 통하여 일기무한모선 계통에 적용한 과도안정도 해석법을 적용하여 왔다. 등면적법을 이용한 과도안정도 해석은 전력위상각 곡선으로부터 임계고장제거 시간(critical clearing time)이나 임계전력위상각을 계산하는 방법이다. 즉 과도안정도 판별을 위해서는 임계전력위상각을 계산해야하는데, 고장 후 발전기의 가속면적와 감속면적이 같은 위상각을 임계전력위상각으로 선정할 수 있다. Until now, such transient stability has been interpreted as a fluctuation equation and equal-area criterion of the infinite busline system, and the transient stability analysis method applied to the weatherless busline system through system contraction has been applied even in large scale systems. The transient stability analysis using the isometry method is a method of calculating the critical clearing time or the critical power phase angle from the power phase angle curve. In other words, the critical power phase angle must be calculated to determine the transient stability. The phase power equal to the acceleration and deceleration areas of the generator after failure can be selected as the critical power phase angle.

등면적법은 일기무한모선에 적용할 수 있는 방법으로 기준모선의 위상각에 대한 상대적인 위상각을 사용하므로 많은 수의 발전기로 구성되어 있는 전력시스템에 직접 적용할 수 없다. 따라서 다기계통에서는 다수의 발전기들을 일기무한모선으로 축약해서 등면적법을 적용하고 있다. 등면적법은 전력위상각 곡선에 기초해서 과도안정도를 판별하므로 반드시 기준발전기에 대한 상대적 위상각에 대한 정보를 알아야 한다. 지금까지 전력시스템에서 과도안정도 판별은 컴퓨터 프로그램의 시뮬레이션에 의존해서 수행하였으며, 최근 동적안전성확보(dynamic security assessment) 측면에서 온라인 과도안정도 평가를 활발하게 수행하고 있으나 상기에 기술한 바와 같이 상대적 위상각을 측정해야 하므로 매우 짧은 시간에 과도안정도 판별을 수행해야 하는 과도안정도 해석에서는 부적당하며, 현재까지 실시간으로 과 도안정도를 판별할 수 기법은 없다.The iso-area method is a method that can be applied to the infinite bus, so it can not be directly applied to the power system composed of a large number of generators because it uses the relative phase angle to the phase angle of the reference bus. Therefore, the multi-machine barrel adopts the uniform area method by abbreviating many generators to infinite bus. Since the isometry method determines the transient stability based on the power phase angle curve, it is necessary to know information about the relative phase angle with respect to the reference generator. Until now, the transient stability of power system has been determined based on the simulation of computer program. Recently, online transient stability assessment has been actively performed in terms of dynamic security assessment. It is inadequate in transient stability analysis, which requires the measurement of transient stability in a very short time, and there is no technique to determine the transient stability in real time.

현재 발전기 인근에 심각한 고장이 발생한 경우, 룩업테이블(lookup table)에 의존해서 경험적으로 발전기들을 탈락시켜 발전기 그룹 전체의 탈조를 방지하고 있다. 그러나 경험적 기능만 부가되어 있고, 현재 전력시스템 상태가 전혀 반영되지 않으므로 불필요하게 많은 발전기를 탈락시키거나 부적절한 탈락으로 그룹 전체 발전기들을 탈조시킬 가능성이 있으며, 이로 인하여 전력시스템 전체를 불안정하게 할 수 있다.In the event of a serious failure in the vicinity of the current generator, it relies on a lookup table to empirically drop generators to prevent outages of the entire generator group. However, since only empirical functions are added and current power system status is not reflected at all, it is possible to drop out many generators unnecessarily or to disassemble group generators by improper dropout, which may cause the entire power system to become unstable.

본 발명의 내용과 관련된 종래의 기술로는 신경회로망을 이용해서 온라인으로 전력시스템의 과도안정 여유를 평가하는 방법(국내특허 제10-0199554호)과 전력시스템에 발생할 수 있는 중대한 고장에 대한 상정사고를 선택해서 발전기들의 상태와 안정도를 예측하는 과도안전도 평가에 대한 방법(국내특허 제10-0411817호), 그리고 복소전력의 시간당 변화율을 이용해서 동기탈조를 검출하는 방법(국내특허 제0599818호)등이 있다.Conventional techniques related to the contents of the present invention include a method for evaluating transient stability margin of a power system online using a neural network (Korean Patent No. 10-0199554) and an assumption about a serious failure that may occur in a power system. Method for evaluating transient safety to predict the state and stability of generators by selecting the (Korean Patent No. 10-0411817), and a method of detecting synchronous outage using the hourly rate of change of complex power (Korean Patent No. 0599818) Etc.

상기한 종래의 과도안정도 판별 방법에서는 등면적법과 상대적 위상각의 필요성으로 인하여 실시간 과도안정도 판별이 불가능하고, 과도상태에서 경험적으로 발전기를 탈락시킴으로써 과도불안정에 대한 대응이 부적절할 수 있는 문제점이 있다.In the conventional transient stability determination method, there is a problem in that it is impossible to determine the transient stability in real time due to the isometry method and the necessity of the relative phase angle, and the response to transient instability may be inappropriate by dropping the generator empirically in the transient state.

이에 본 발명은 상기와 종래기술의 문제점을 해소하기 위해 안출한 것으로, 위상각이나 상대적 위상각을 사용하지 않고 실측 전압과 전류로부터 계산한 전력으 로부터 얻은 가속력과 가속에너지를 사용하여 과도안정도를 판별하고 대응한다. 모든 과정을 과도안정도 해석에 허용할 수 있는 작은 시간 단위의 실시간으로 처리할 수 있도록 매우 단순한 과도안정도 판별 방법을 제시하고, 과도불안정시 탈락시킬 발전기를 선택하여 실시간으로 과도불안정을 제거하는 방법을 제공하는 데에 그 목적이 있다.Accordingly, the present invention has been made to solve the problems of the above and the prior art, and the transient stability is determined using the acceleration force and the acceleration energy obtained from the power calculated from the measured voltage and current without using the phase angle or the relative phase angle. And respond. It provides a very simple transient stability discrimination method to handle all processes in real time in small time units that can be accepted for transient stability analysis, and provides a method of eliminating transient instability in real time by selecting a generator to drop out in case of transient instability. Its purpose is to.

상기 목적을 달성하기 위한 본 발명의 실시간 전력시스템 과도안정도 판별과 고장파급방지 방법은, 실시간으로 발전기의 전압, 전류와 기계적 입력을 취득하는 스텝과, 상기 취득된 이산치 전압과 전류로부터 전력을 실시간으로 계산하여 고장을 판단하는 하는 스텝과, 고장이라 판단되면 기계적 입력(P_m)과 상기에서 계산한 전력(P_e)의 편차로부터 가속력(P_a)을 계산하는 스텝과, 상기 계산된 가속력를 일정시간동안 적분하여 가속에너지를 계산하는 스텝과, 상기 산출된 가속에너지와 설정값을 비교하는 스텝과, 상기 비교하는 스텝에서 가속에너지가 설정값의 지정된 범위를 벗어나면 과도불안정으로 판단하여 발전기를 탈락하는 스텝과, 상기 결과를 출력 및 저장하는 스텝을 포함하는 것을 특징으로 한다.The real-time power system transient stability determination and fault spreading prevention method of the present invention for achieving the above object comprises the steps of acquiring the voltage, current and mechanical input of the generator in real time, and real-time power from the obtained discrete voltage and current; The step of determining the failure by calculating the step, and the step of calculating the acceleration force (P _a ) from the deviation between the mechanical input (P _m ) and the calculated power (P _e ) if the failure is determined, and the calculated acceleration force A step of calculating acceleration energy by integrating for a time, comparing the calculated acceleration energy with a set value, and dropping the generator by judging transient instability if the acceleration energy is out of a specified range of the set value in the comparing step. And a step of outputting and storing the result.

이하, 본 발명의 원리를 첨부된 도면을 참조하여 상세히 설명하기로 한다.Hereinafter, the principle of the present invention will be described in detail with reference to the accompanying drawings.

[본 발명을 수행하기 위한 하드웨어][Hardware for Carrying Out the Invention]

도 1은 본 발명을 수행하기 위한 하드웨어의 개략적인 블록 구성도로서, 동도면을 참조하면 알 수 있듯이, 본 발명을 수행하기 위한 하드웨어는, 실시간 데이 터 취득부(10)와 입력 처리부(20)와 연산부(40)와 제어부(60)와 데이터 출력 및 저장부(70)를 포함하여 구성될 수 있다.1 is a schematic block diagram of hardware for carrying out the present invention. As can be seen from the drawings, the hardware for carrying out the present invention includes a real-time data acquisition unit 10 and an input processing unit 20. And a calculator 40, a controller 60, and a data output and storage unit 70.

상기 입력 처리부(20)는 입력되는 아날로그형태의 실시간 시계열 전압, 전류 데이터를 아날로그 투 디지털 변환(analog to digital convert)하여 연산부(40)로 입력한다.The input processor 20 analog-to-digital converts the real-time time series voltage and current data in an analog form and inputs the same to the calculator 40.

상기 설정부(30)는 과도불안정을 판별하기 위한 비교값으로 사용하기 위한 가속에너지나 누적가속력을 설정한다.The setting unit 30 sets acceleration energy or cumulative acceleration force for use as a comparison value for determining transient instability.

상기 연산부(40)는 입력 처리부(20)로부터 입력되는 이산치 전압과 전류로부터 전력과 가속력 및 가속에너지를 계산하고 과도안정도를 실시간으로 판별한다.The calculation unit 40 calculates power, acceleration force and acceleration energy from the discrete voltage and current input from the input processing unit 20 and determines the transient stability in real time.

상기 제어부(60)는 연산부의 결과에 따라 과도불안정시 최대가속에너지를 갖는 발전기를 정지한다.The control unit 60 stops the generator having the maximum acceleration energy in case of transient instability according to the result of the calculation unit.

상기 데이터 출력 및 저장부(70)는 상기 연산부(40)에 의해 계산된 결과와 제어동작의 결과 및 실시간으로 취득한 이산데이터를 출력하고 저장 장치에 저장한다.The data output and storage unit 70 outputs the result calculated by the operation unit 40, the result of the control operation, and the discrete data acquired in real time, and stores them in the storage device.

이하에서는, 상기한 바와 같이 구성된 하드웨어를 통해 수행되는 본 발명의 구체적인 내용에 대해 설명하기로 한다.Hereinafter, specific details of the present invention performed through the hardware configured as described above will be described.

본 발명은, 전력시스템에서 실시간으로 과도안정도를 판별하고 고장이 파급되는 것을 방지하는 방법에 관한 것으로, 실시간으로 취득한 이산데이터에서 전력을 계산하고, 이로부터 고장을 판단한다. 고장이 발생하였을 때, 실측한 기계적 입력과 전력의 편차로부터 가속력을 계산하고, 가속력으로부터 발전기에 누적되는 가 속에너지를 계산한다. 가속에너지가 미리 설정된 설정값의 범위를 벗어나면 과도불안정으로 판단하여 발전기를 탈락함으로써 다른 발전기들을 정상운전을 가능케 한다. The present invention relates to a method for determining transient stability in a real time in a power system and preventing a failure from spreading. The present invention calculates power from discrete data obtained in real time and determines a failure therefrom. When a failure occurs, the acceleration force is calculated from the measured mechanical input and power deviation, and the acceleration energy accumulated in the generator is calculated from the acceleration force. If the acceleration energy is out of the range of the preset setting value, it is judged as transient instability and the generator is dropped so that other generators can operate normally.

우선, 본 발명의 실시예에 적용되는 원리에 대해 설명하기로 한다.First, the principle applied to the embodiment of the present invention will be described.

지금까지 이러한 과도안정도는 일기무한모선 계통의 동요방정식과 등면적법(equal-area criterion)으로 해석되고, 대규모 계통에서도 계통축약을 통하여 일기무한모선 계통에 적용한 과도안정도 해석법을 적용하여 왔다. 기계적 입력과 전기적 출력을 각각 P_m, P_e라 하고, 관정정수와 회전자 위상각을 각각 H, δ라 하면, 동요방정식(swing equation)은 수학식1과 같다.Until now, such transient stability has been interpreted as a fluctuation equation and equal-area criterion of the infinite busline system, and the transient stability analysis method applied to the weatherless busline system through system contraction has been applied even in large scale systems. If the mechanical input and the electrical output are P _m and P _e , and the tube constant and the rotor phase angle are H and δ, respectively, the swing equation is as shown in Equation 1.

그리고 발전기 내부전압과 단자전압을 각각 E_f, V_t라 하고, 내부리액턴스를 X_T라 하면, 전기적 출력 P_e는 위상각의 비선형 함수인 수학식2와 같다.When the generator internal voltage and the terminal voltage are E _f and V _t , respectively, and the internal reactance is X _T , the electrical output P _e is expressed as Equation 2, which is a nonlinear function of the phase angle.

발전기의 기계적 입력은 빠르게 변화할 수 없으므로 외란 후, 위상각의 최대 변화는 전기적 출력 P_e≒0일 때 발생한다. 즉 외란 발생 후 초기동요는 전기적 출력 P_e에 의존하게 되고 따라서 초기과도상태의 동요는 정현파함수가 된다.Since the mechanical input of the generator cannot change quickly, the maximum change in phase angle after disturbance occurs at the electrical output P _e ≒ 0. In other words, after the disturbance occurs, the initial fluctuation depends on the electrical output P _e , and thus the initial transient fluctuation becomes the sinusoidal function.

동요방정식에서 등면적법은 임계고장제거 시간(critical clearing time) 또는 임계전력위상각을 계산하는 방법이다. 도2는 등면적법을 해석하기 위해서 나타낸 고장전후의 전력위상각 곡선이다. 도2에서 a점에서 정상운전 되고 있는 발전기가 고장이 발생하면 b점에서 c점으로 운전된다. 그리고 고장을 제거한 후에는 e점에서 f점을 따라 운전된다. 등면적법에 의한 과도안정도 판별을 위해서는 임계전력위상각을 계산해야하는데, 도2에서 가속면적 A_acc와 감속면적 A_dcc가 같을 때, 위상각 δ_c를 계산함으로써 임계전력위상각을 계산할 수 있다.In the fluctuation equation, the isometry method is a method of calculating the critical clearing time or the critical power phase angle. Fig. 2 is a power phase angle curve before and after failure shown for analyzing the isoarea method. In FIG. 2, when a generator in normal operation at point a fails, it is operated from point b to point c. After the fault is cleared, it operates from point e to point f. To determine the transient stability by the isometry method, the critical power phase angle must be calculated. When the acceleration area A _acc and the deceleration area A _dcc are the same in FIG. 2, the critical power phase angle can be calculated by calculating the phase angle δ _c . .

등면적법을 일기무한모선에 적용할 수 있는 방법으로 많은 수의 발전기로 구성되어 있는 전력시스템에 직접 적용할 수 없다. 따라서 다기계통에서는 다수의 발전기들을 일기무한모선으로 축약해서 등면적법을 적용하고 있다. It is not possible to apply the iso-area method directly to power systems consisting of a large number of generators. Therefore, the multi-machine barrel adopts the uniform area method by abbreviating many generators to infinite bus.

등면적법은 전력위상각 곡선(도2)에 기초해서 과도안정도를 판별하고 있다. 따라서 반드시 기준발전기에 대한 상대적 위상각에 대한 정보를 알아야 한다. 지금까지 전력시스템에서 과도안정도 판별은 컴퓨터 프로그램의 시뮬레이션에 의존해서 수행하였으며, 실시간으로 과도안정도를 판별하는 기법은 없다.The iso-area method discriminates transient stability based on the power phase angle curve (Fig. 2). Therefore, you must know information about the relative phase angle for the reference generator. So far, transient stability determination in power system is based on simulation of computer program, and there is no technique to determine transient stability in real time.

본 발명에서는 전력시스템에서 실시간으로 과도안정도를 판별하기 위한 새로운 개념으로써 가속에너지를 도입하였는데, 먼저 수학식1에 나타나있는 동요방정식(swing equation)과 가속에너지의 관계에 대해서 기술한다. In the present invention, the acceleration energy is introduced as a new concept for determining the transient stability in real time in the power system. First, the relationship between the swing equation and the acceleration energy shown in Equation 1 will be described.

동요방정식에서 외란에 의해서 전기적출력 P_e가 변화하면 발전기 내부 위상각이 변화한다. 그러므로 기계적 입력 P_m과 전기적 출력 P_e의 편차를 가속력(또는 감속력) P_a라 하면, 가속력 P_a=P_m-P_e가 성립하고 도3에서와 같이 나타낼 수 있다. 도3에서 알 수 있는 것과 같이 기계적 입력 P_m이 전기적 출력 P_e보다 크면(P_m>P_e), 발전기는 가속되어 발전기 내부에 가속에너지가 축적되고, 반대로 기계적 입력 P_m이 전기적 출력 P_e보다 작으면(P_m<P_e), 발전기는 감속되어 발전기 내부에 감속에너지가 축적된다. 이와 같은 가속력 및 감속력은 전기적 출력의 변화에 대한 발전기 내부에 저장되는 힘을 의미한다. In the fluctuation equation, when the electric output P _e is changed by disturbance, the phase angle inside the generator is changed. Therefore, if the deviation between the mechanical input P _m and the electrical output P _e is acceleration force (or deceleration force) P _a , acceleration force P _a = P _m -P _e is established and can be expressed as shown in FIG. 3. As can be seen in FIG. 3, when the mechanical input P _m is greater than the electrical output P _e (P _m > P _e ), the generator is accelerated to accumulate the acceleration energy inside the generator, on the contrary, the mechanical input P _m is the electrical output P _e. If smaller (P _m <P _e ), the generator is decelerated and deceleration energy is accumulated inside the generator. Such acceleration and deceleration forces mean the forces stored inside the generator for changes in electrical output.

따라서 일정시간 동안 가속력을 적분하면 발전기 내부에 축적되는 가속(감속)에너지에 대한 식을 얻을 수 있다. 가속력(감속력)을 P_a라 하고, 가속(감속)에너지를 E_a라 할 때, 시간 t₁초부터 t₂초까지 발전기 내부에 축적되는 가속에너지는 수학식3과 같이 나타낼 수 있다.Therefore, by integrating the acceleration force for a certain time, the equation for the acceleration (deceleration) energy accumulated in the generator can be obtained. When the acceleration force (deceleration force) is P _a and the acceleration energy (deceleration) energy is E _a , the acceleration energy accumulated in the generator from the time t ₁ second to t ₂ seconds can be expressed as Equation (3).

발전기 모선 인근에 고장이 발생하였을 때, 전기적 출력과 기계적 출력의 파형은 도4와 같다. 도4에서 외란이 발생한 후, 전기적 출력(410, 420)은 급변하는 반면, 기계적 출력(430)은 미세하게 변화함을 알 수 있다. 도4에서 전기적 출력(410)과 기계적 출력(430) 사이의 면적(400)이 발전기 내부에 축적된 가속에너지를 나타내고 있다. When a failure occurs near the generator busbar, the waveforms of the electrical output and the mechanical output are shown in FIG. In FIG. 4, after the disturbance occurs, the electrical outputs 410 and 420 change rapidly while the mechanical outputs 430 change slightly. In FIG. 4, the area 400 between the electrical output 410 and the mechanical output 430 represents the acceleration energy accumulated in the generator.

한편 발전기의 전기적 출력은 고장위치에 따라서 출력되는 파형이 다르다. 발전기 모선(300)에서 고장이 발생하면, 수학식2에서 모선 전압 V_t≒0이므로 전기적 출력은 거의 영에 가까게 감소한다(420). 반면에 발전기모선에서 비교적 원거리에서 고장이 발생하면, 모선 전압 V_t≠0이고, 고장위치와 발전기 사이의 부하에 전력을 공급하므로 전기적 출력은 감소하지만 영이 되지는 않는다(410).On the other hand, the electrical output of the generator is different from the waveform outputted according to the fault location. When a failure occurs in the generator bus 300, the bus voltage V _t ≒ 0 in Equation 2 decreases the electrical output to near zero (420). On the other hand, if a fault occurs in the generator bus at a relatively long distance, the bus voltage V _t ≠ 0, and the electric power is reduced but not zero since the power is supplied to the load between the fault position and the generator (410).

전력시스템에서 운전되고 있는 각 발전기가 수용할 수 있는 수학식3으로 표현되는 가속에너지의 양은 한정되어 있는데, 이것은 계통조건에 따라서 다르다. 고장 발생 시 수용할 수 있는 에너지의 최소량은 발전기 모선(300)에 고장이 발생하였을 때인데, 이것은 도3에서 전기적출력 P_e가 거의 0에 근접하여 있으므로 모든 가속에너지가 발전기내부에 축적되기 때문이다. 즉 짧은 시간에 발전기에 축적되는 에너지가 큰 값을 가지므로 수학식1로 표현되는 발전기 동요방정식에서 회전자속도(전력위상각 변화)가 순간적으로 빠르게 변화하고 발전기 동요에 민감하게 영향을 준다. 그러므로 발전기 모선 고장에서 발전기가 과도안정도를 유지하기 위한 가속에너지와 그 설정치를 얻을 수 있다.The amount of acceleration energy represented by Equation 3 that each generator operating in the power system can accommodate is limited, which depends on the system conditions. The minimum amount of energy that can be accommodated in the event of a failure is when a failure occurs in the generator bus 300. This is because in FIG. 3, all acceleration energy is accumulated in the generator since the electrical output _Pe is close to zero. . That is, since the energy accumulated in the generator in a short time has a large value, the rotor speed (change of the power phase angle) changes rapidly instantaneously in the generator shaking equation represented by Equation 1 and sensitively affects the generator shaking. Therefore, in case of generator bus failure, it is possible to obtain the acceleration energy and its set value for maintaining the transient stability of the generator.

본 발명은 상기 발전기 내부에 축적되는 가속에너지에 의해서 실시간 과도안 정도를 판별하는 방법에 관한 것이다.The present invention relates to a method for determining the real-time transient degree by the acceleration energy accumulated in the generator.

전력시스템에서 운전되고 있는 발전기들은 다수의 발전기들이 병렬연결 되어 그룹을 이루어 운전되고 있다. 따라서 각 발전기들은 상호 영향을 주고받으면서 안정한 동작점에서 운전되고 있기 때문에 네트워크 구조 및 고장의 위치에 따라서 각 발전기가 수용할 수 있는 가속에너지는 다르다. 임계 가속에너지(critical acceleration energy)는 임계상태에서의 발전기에 축적된 가속에너지로 발전기 인근 모선 고장 시에 최소가 되고, 발전기에서 원거리 고장 시에 발전기 인근 모선 고장 시 가속에너지보다 큰 값을 가지나 큰 차이는 없다. Generators in the power system are operated in groups with multiple generators connected in parallel. Therefore, since each generator operates at a stable operating point while exchanging influences with each other, the acceleration energy that each generator can accommodate depends on the network structure and the location of the failure. Critical acceleration energy is the acceleration energy accumulated in the generator in the critical state, which is the minimum at the failure of the bus near the generator, and has a larger value than the acceleration energy at the failure of the bus near the generator. There is no.

에너지보존법칙에 의해서 기계적으로 입력된 에너지(E_m)는 가속에너지(E_a)와 전기적출력 에너지(E_e)의 합과 같으므로 수학식4가 성립한다(도5). 도5에서 기계적 입력 에너지(E_m)가 일정하고, 네트워크에 공급된 에너지 E_e=0일 때, 가속에너지 E_a는 최대가 됨을 알 수 있다.Since the energy E _m mechanically input by the energy conservation law is equal to the sum of the acceleration energy E _a and the electrical output energy E _e , Equation 4 holds (FIG. 5). In FIG. 5, it can be seen that when the mechanical input energy E _m is constant and the energy E _e = 0 supplied to the network, the acceleration energy E _a becomes maximum.

따라서 과도상태에서 기계적 입력 에너지(E_m)가 일정하다고 할 때, 가속에너지(E_a)는 전기적출력 에너지(E_e)에 따라서 변화하므로 과도상태에서 가속에너지는 고장의 위치에 의해서 달라진다. 발전기 모선근처 고장이 발생하면, 소비되는 부하전력이 거의 없기 때문에 전기적 출력 에너지가 0에 근접한다(E_e≒0). 반면 발전기 모선에서 비교적 원거리에 고장이 발생하면 고장지점과 발전기 사이에 부하전력이 있기 때문에 전기적 출력 에너지는 0이 되지 않고 부하량에 비례한 전기적출력 에너지가 발생한다.Therefore, when the mechanical input energy E _m is constant in the transient state, the acceleration energy E _a varies according to the electrical output energy E _e , so the acceleration energy in the transient state varies depending on the location of the failure. If a failure occurs near the generator busbar, the electrical output energy is close to zero because little load power is consumed (E _e ≒ 0). On the other hand, if a fault occurs relatively far in the generator bus, there is a load power between the point of failure and the generator, so that the electrical output energy is not zero and the electrical output energy is proportional to the load.

도6에서 E_e(610)는 고장이 발생하였을 때, 과도상태의 전기적 출력 에너지, 즉 고장상태에서 부하에 공급한 에너지를 나타내고 있고, E_a(600)는 발전기 내부에 축적되는 가속에너지를 나타내고 있다. 가속에너지(600)와 전기적 출력 에너지(610)가 비슷한 값을 가지고 있어 발전기 모선 고장이 아님을 알 수 있다.In FIG. 6, E _e 610 represents the electrical output energy in the transient state, that is, the energy supplied to the load in the fault state when the failure occurs, and E _a 600 represents the acceleration energy accumulated in the generator. have. It can be seen that the acceleration energy 600 and the electrical output energy 610 have similar values, so that the generator bus is not a failure.

따라서 전기적 출력에너지(E_e)와 가속에너지(E_a) 사이의 관계에서 과도안정도를 판별할 수 있다. 과도안정도를 판별하기 위해서 각 발전기가 수용할 수 있는 가속에너지의 임계값을 결정해야한다. 상기한 바와 같이 발전기 모선 고장에서 E_e≒0이므로 이때 발전기에는 순수한 가속에너지만 축적된다. 따라서 특정 발전기가 수용할 수 있는 임계가속에너지(critical acceleration energy)는 발전기 모선 고장에서 얻을 수 있다.Therefore, it is possible to determine the transient stability in the relationship between the electrical output energy (E _e ) and the acceleration energy (E _a ). To determine the transient stability, it is necessary to determine the threshold of acceleration energy that each generator can accommodate. As described above, since E _e ≒ 0 in the generator bus failure, only pure acceleration energy is accumulated in the generator. Therefore, the critical acceleration energy that can be accommodated by a specific generator can be obtained from the generator bus failure.

가속에너지를 이용하여 실시간으로 전력시스템의 과도안정도를 판별하는데 있어, 가속에너지에 대한 기준값(설정값)이 필요하다. 설정값은 전력시스템 시뮬레이션 프로그램에서 발전기 모선에 대한 고장을 모의하고, 각 발전기가 수용할 수 있는 임계가속에너지를 계산한 후, 이를 실시간 과도안정도 판별을 위한 가속에너지 설정값으로 사용할 수 있다. 또는 발전기 모선에 고장이 발생하였을 때, 실측한 이산데이터에서 계산한 임계가속에너지를 설정값으로 사용할 수 있다. 발전기 모선 과 선로를 비롯해서 다양한 위치에 대한 고장을 모의해서 얻은 임계가속에너지들 중에서 최소값을 설정값으로 사용할 수 있는데, 다양한 상정사고에 대한 가속에너지를 E_a라 하면, 설정값 ε은 다음과 같이 정의한다.In determining the transient stability of the power system in real time using the acceleration energy, a reference value (set value) for the acceleration energy is required. The set value simulates the failure of the generator bus in the power system simulation program, calculates the critical acceleration energy that each generator can accept, and then uses it as the acceleration energy set value for real-time transient stability determination. Alternatively, when a failure occurs in the generator bus, the critical acceleration energy calculated from the measured discrete data can be used as the set value. The minimum value among the critical acceleration energies obtained by simulating faults at various locations including generator buses and lines can be used as the set value.If the acceleration energy for various assumed accidents is E _a , the set value ε is defined as do.

설정값 ε은 전력시스템의 외란에 대해서 특정발전기가 안정할 수 있는 임계값을 나타낸 것이다. 이 값은 발전기 모선 고장에 대한 값이므로 다른 발전기의 영향을 무시한 경우이므로 동일한 시간에 특정 발전기에 누적될 수 있는 가속에너지는 최대치가 된다. 그리고 임계가속에너지는 특정 발전기가 동기탈조를 일으키는 최소의 축적된 가속에너지이다.The set value ε represents the threshold at which a particular generator can be stabilized against disturbances in the power system. Since this value is for generator bus faults, neglecting the influence of other generators, the acceleration energy that can accumulate in a specific generator at the same time is the maximum. And the critical acceleration energy is the minimum accumulated acceleration energy that causes a specific generator to cause synchronous stoppage.

실시간 전력시스템에서 순시치 전압과 전류를 비롯하여 다양한 형태의 데이터 취득이 가능하다. 실시간 과도안정도는 보통 수 사이클(cycle, 1cycle=1/60초) 이내에 판단 및 해석이 가능해야 그 의미가 있는데, 따라서 기준발전기에 대한 위상각을 필요로 하는 등면적법은 실시간 적용이 어렵다. 반면에 가속에너지를 사용하면 실측되는 실시간의 이산데이터를 이용하므로 실시간 과도안정도 판별이 가능하다.Various types of data can be acquired, including instantaneous voltage and current, in real-time power systems. Real-time transient stability usually means that it can be judged and interpreted within a few cycles (cycle, 1cycle = 1/60 seconds). Therefore, real-time application of the isometry method, which requires a phase angle to the reference generator, is difficult. On the other hand, if acceleration energy is used, real-time transient stability can be determined because the measured real-time discrete data is used.

만일 전력시스템에서 실시간 데이터를 1 사이클 이하로 샘플링 하였을 때, 고장이 판단되면 즉시 발전기 회전자에 축적되는 에너지를 실시간으로 계산할 수 있다. 발전기에 축적되는 에너지는 가속력의 적분으로 계산할 수 있는데, 시간 t_s~t_e초 사이에 가속에너지는 수학식6과 같다. If real-time data is sampled in the power system less than 1 cycle, if the failure is determined, the energy accumulated in the generator rotor can be calculated in real time. The energy accumulated in the generator can be calculated as the integral of the acceleration force, and the acceleration energy is expressed by Equation 6 between the time t _s ~ t _e seconds.

실시간으로 취득한 데이터는 이산치이므로 P_a[0]를 과도상태가 처음 시작된 데이터에서 계산한 가속력이라 하면, 이산데이터에서 수학식6의 적분은 수학식7로 계산할 수 있다.Since the data acquired in real time is a discrete value, if P _a [0] is the acceleration force calculated from the data at which the transient state first started, the integral of Equation 6 in the discrete data can be calculated by Equation 7.

수학식7에서 N은 과도상태 후 취득한 데이터 수이고 Δt는 샘플링 구간이다. 그러므로 실시간 과도안정도 판별은 과도상태 임계가속에너지의 설정치 ε과 실측한 이산데이터에 의해서 계산한 가속에너지 E_a[n]을 비교해서 판별할 수 있다. 즉 수학식8과 같이 실측 계산한 가속에너지 E_a[n]이 설정된 임계가속에너지보다 크면 발전기는 과도불안정으로 판별할 수 있다. 그리고 실측 계산한 가속에너지 E_a[n]이 설정치 ε보다 작을 때, 안정으로 판별할 수 있다. 만일 고장 후 몇 사이클 이내에 고장이 제거되고, 감속에너지가 유입되면 발전기는 과도안정 상태로 운전을 지속할 수 있다. In Equation 7, N is the number of data acquired after the transient state and Δt is the sampling interval. Therefore, real-time transient stability determination can be determined by comparing the setpoint ε of the transient critical acceleration energy with the acceleration energy E _a [n] calculated from the measured discrete data. That is, as shown in Equation 8, if the measured acceleration energy E _a [n] is greater than the set threshold acceleration energy, the generator may determine transient instability. When the calculated acceleration energy E _a [n] is smaller than the set value ε, it can be determined as stable. If the fault is cleared and the deceleration energy is introduced within a few cycles after the fault, the generator can continue to operate in transient stability.

여기에서 E_a[n]은 고장발생 직후부터 n-번째 데이터까지 발전기에 축적된 가속에너지이다. 보통 고장제거는 고장 후 몇 사이클 이내에 이루어지는데, 샘플링 시간이 등간격으로 주어지므로 고장제거시간은 쉽게 계산할 수 있다. 예를 들면 데이터 샘플링을 0.5사이클로 할 때, n-번째 데이터에서 고장을 제거하면, 고장제거시간은 n/2 사이클이 된다. Where E _a [n] is the accelerating energy accumulated in the generator from just after the failure to the n-th data. Normally, fault elimination occurs within a few cycles of a fault. Sampling times are given at equal intervals, so fault elimination time can be easily calculated. For example, when the data sampling is 0.5 cycles and the fault is removed from the n-th data, the fault removal time is n / 2 cycles.

전력시스템에서는 동일 모선이나 인근 모선에 몇 기의 발전기들이 병렬로 연결되어 운전하고 있다. 동일 모선이나 인근 모선에 다기발전기가 병렬 운전될 때, 특정발전기가 과도불안정상태에 이르면 특정발전기를 제거해야 한다. 병렬 운전의 특성상 과도불안정에 이르면 다른 발전기들과 동시에 동기탈조를 유발하기 때문이다. 이때 탈락시켜야 할 발전기는 발전기 내부에 축적된 가속에너지가 가장 큰 발전기를 선택해야 다른 발전기 그룹의 탈조를 방지하는데 효과가 크다. 이것은 발전기 그룹의 전체 가속에너지에서 가장 큰 가속에너지를 제거하면 전체가속에너지는 상대적으로 크게 감소하고 또한 임계가속에너지에 대한 여유가 증가하기 때문이다. In the power system, several generators operate in parallel with the same bus or with a nearby bus. When multiple generators are operated in parallel on the same bus or on a nearby bus, the generator must be removed when the generator is transiently unstable. This is because when the transient instability is reached due to the characteristics of parallel operation, it causes synchronous breakdown at the same time as other generators. At this time, the generator to be dropped should be selected the generator with the largest acceleration energy accumulated in the generator, it is effective to prevent the stepping out of other generator groups. This is because if the largest acceleration energy is removed from the total acceleration energy of the generator group, the total acceleration energy decreases relatively and the margin for the critical acceleration energy increases.

발전기 그룹에서 i-번째 발전기의 가속에너지를 E_ai[n]라 할 때, 탈락시킬 발전기는 수학식9와 같이 최대 임계가속에너지를 가진 발전기를 선택한다. When the acceleration energy of the i-th generator in the generator group is E _ai [n], the generator to be dropped out selects the generator having the maximum critical acceleration energy as shown in Equation (9).

예를 들면 m-기의 발전기가 그룹 되어 있고, 이 그룹에 축적된 가속 에너지를 E_s[n]이라 하면 수학식10이 성립한다.For example, when m-group generators are grouped, and the acceleration energy accumulated in this group is E _s [n], Equation 10 is established.

만일 발전기 그룹에서 가장 큰 가속에너지 E_sel이 그룹 전체 가속에너지의 합 E_s[n]의 30%라 가정하였을 때, 이를 제거하면 발전기 그룹 측면에서는 0.3E_s[n] 만큼 가속에너지의 여유가 증가하므로 과도안정도의 여유가 증가한다. 도7에는 고장 중 시간에 대한 전력을 나타내고 있다. 도7에서 시간 t₀~t₁초 사이에 그룹발전기에 축적에너지는 A₁이고, 시간 t₁~t₂초 및 t₂~t₃초 사이에 발전기에 축적에너지는 각각 A₂와 A₃이다. 고장 중에 전기적 출력 변화가 작으면 축적된 가속에너지는 적분의 합에 대한 연산이 선형이기 때문에 도8과 같이 선형적으로 증가한다. 도8은 시간에 대한 발전기에 축적된 가속에너지를 나타내고 있는데, E₁, E₂, E₃은 가속에너지를 나타내고 있고, E_s는 각각 발전기그룹에 전체 가속에너지를 나타내고 있다. 그리고 E_sel은 발전기그룹에서 임의의 특정발전기에 축적된 가속에너지를 의미한다. 도8에 나타난 것과 같이 발전기 그룹 전체의 가속에너지에서 특정발전기에 축적된 가속에너지를 빼면 발전기 그룹에 남아있는 가속에너지 E_res가 된다. 따라서 과도상태에서 발전기 그룹의 과도안정도를 향상하기 위하여 특정발전기를 제거할 때에는 가속에너지가 가장 큰 발전기를 탈락시키는 것이 그룹의 과도안정도를 가장 크게 향상시킬 수 있다. 도8에서 E_sel이 크면 클수록 발전기 그룹에 남아있는 가속에너지 E_res가 작게 되므로 발전기 그룹의 과도안정도 여유가 커지게 된다.If the largest acceleration energy E _sel in the generator group is assumed to be 30% of the sum E _s [n] of the group's total acceleration energy, removing it increases the margin of acceleration energy by 0.3E _s [n] on the generator group side. Therefore, the margin of transient stability increases. Fig. 7 shows power versus time during failure. In FIG. 7, the accumulated energy in the group generator is A ₁ between the times t ₀ to t ₁ seconds, and the accumulated energy in the generator is between A ₂ and A ₃ between the times t ₁ to t ₂ seconds and t ₂ to t ₃ seconds, respectively. . If the change in electrical output during failure is small, the accumulated acceleration energy increases linearly as shown in Fig. 8 because the calculation for the sum of the integrals is linear. 8 shows the acceleration energy accumulated in the generator over time, E ₁ , E ₂ and E ₃ represent the acceleration energy, and E _s represents the total acceleration energy in the generator group, respectively. And E _sel refers to the acceleration energy accumulated in any particular generator in the generator group. As shown in FIG. 8, the acceleration energy accumulated in a specific generator is subtracted from the acceleration energy of the entire generator group, resulting in the acceleration energy E _res remaining in the generator group. Therefore, when removing a specific generator in order to improve the transient stability of the generator group in the transient state, dropping the generator with the largest acceleration energy can greatly improve the transient stability of the group. In Fig. 8, the larger the E _sel is, the smaller the acceleration energy E _res remaining in the generator group becomes, so that the transient stability margin of the generator group becomes larger.

발전기 그룹에서 과도안정도는 동일모선에 연결된 발전기 또는 하나의 발전기의 과도안정도에 좌우된다. 그러므로 그룹내부에 포함된 각 발전기에 대한 임계시간을 계산할 필요가 있다. 도8에서는 발전기 그룹에 대한 가속에너지를 나타내고 있다. 만일 도8에서 E₁, E₂, E₃가 특정발전기에 대한 축적된 가속에너지라 하면, 선형특성을 이용해서 임계시간을 추정할 수 있다. 설정된 가속에너지를 ε이라 하면, 임계시간 t_c는 수학식11로부터 계산할 수 있다. The transient stability in a generator group depends on the transient stability of the generator or one generator connected to the same bus. Therefore, it is necessary to calculate the threshold time for each generator included in the group. 8 shows the acceleration energy for the generator group. If E ₁ , E ₂ , and E ₃ in FIG. 8 are accumulated acceleration energy for a specific generator, the critical time may be estimated using the linear characteristic. If the set acceleration energy is ε, the threshold time t _c can be calculated from Equation (11).

지금까지는 실측한 전압과 전류로부터 계산한 전기적 출력 P_e로 모든 방정식을 기술하였다. 또한 전기적 출력에 동기속도를 곱하면 전기적 토크 T_e를 계산할 수 있다. 특히 전기적 출력(P_e)과 토크(T_e)를 단위법(per-unit)으로 나타내면 동기속도 ω_s=1.0이므로 전기적출력과 토크는 같다. 따라서 지금까지 기술한 P_e대신에 토크 T_e를 사용해도 된다.So far all equations have been described by the electrical output P _e calculated from the measured voltage and current. The electrical torque T _e can be calculated by multiplying the electrical output by the synchronous speed. In particular, if the electrical output (P _e ) and torque (T _e ) are expressed in per-unit, the synchronous speed ω _s = 1.0, so the electrical output and torque are the same. Therefore, you may use torque T _e instead of P _e described so far.

한편 과도상태에서 이산치 전기적 출력은 데이터가 입력되는 순간마다 누적할 수 있다. 누적된 가속력을 P_s[n]이라 하면, N개의 데이터에 대한 누적된 가속력은 수학식13으로부터 계산할 수 있다.On the other hand, in the transient state, discrete electrical output can accumulate every time data is input. If the accumulated acceleration force is P _s [n], the accumulated acceleration force for N pieces of data can be calculated from Equation 13.

그리고 이산치 데이터에서 샘플링 간격이 등간격이므로 가속에너지와 누적된 전기적출력이나 토크는 비례한다. In addition, since the sampling intervals are equally spaced in the discrete data, the acceleration energy and the accumulated electrical output or torque are proportional to each other.

따라서 상기 설정치를 가속력 P_s[n]을 사용해도 된다. 즉 설정치를 가속력 P_s[n]을 기준으로 사용하면 실시간으로 입력된 데이터에서 수학식13으로 표현되는 누적된 가속력과 비교하여 과도안정도를 판별할 수 있다. 즉 설정치 ε을 누적된 가속력 P_s[n]으로 설정하였을 때, 실측된 데이터로부터 누적한 가속력 P_s[n]과 설정 치 ε을 비교하여 과도안정도를 판별할 수 있다. 실측 계산한 가속력 P_s[n]이 설정치 ε보다 작을 때, 시스템은 안정으로 판별할 수 있고 반대로 수학식15와 같이 설정치 ε보다 클 때, 시스템은 불안정으로 판별할 수 있다. Therefore, the set value may be used for accelerating force P _s [n]. That is, when the set value is used as the acceleration force P _s [n] as a reference, the transient stability can be determined by comparing the accumulated acceleration force represented by Equation 13 in the data input in real time. That is, when the set value ε is set to the accumulated acceleration force P _s [n], the transient stability can be determined by comparing the accumulated acceleration force P _s [n] and the set value ε from the measured data. When the measured calculated acceleration force P _s [n] is smaller than the set value ε, the system can be determined as stable and, conversely, when the calculated acceleration force P _s [n] is larger than the set value ε as in Equation 15, the system can be determined as unstable.

여기에서 P_s[n]은 고장발생 직후부터 n-번째 데이터까지 누적한 가속력이고, 설정치는 과도안정도 시뮬레이션 프로그램 결과나 과거의 외란에서 실측한 이산데이터에서 계산한 누적가속력이다. 보통 고장제거는 고장 후 몇 사이클 이내에 이루어지는데, 샘플링 시간이 등간격으로 주어지므로 고장제거시간은 쉽게 계산할 수 있다. 예를 들면 데이터 샘플링을 0.5사이클로 할 때, n-번째 데이터에서 고장을 제거하면, 고장제거시간은 n/2 사이클이 된다. 즉 샘플링율을 T 사이클이라 하면 n-번째 데이터의 고장제거시간은 수학식16과 같다.Here, P _s [n] is the acceleration force accumulated from the immediately after the failure to the n-th data, and the set value is the accumulated acceleration calculated from the results of the transient stability simulation program or the discrete data measured in the past disturbance. Normally, fault elimination occurs within a few cycles of a fault. Sampling times are given at equal intervals, so fault elimination time can be easily calculated. For example, when the data sampling is 0.5 cycles and the fault is removed from the n-th data, the fault removal time is n / 2 cycles. That is, if the sampling rate is T cycles, the fault elimination time of the n-th data is expressed by Equation 16 below.

[실시예1]Example 1

이하는 상기한 본 발명의 내용을 실시예를 들어 설명한다. 도9는 6기의 발전기가 2기씩 동일 모선에 연결되어 있고, 모선과 모선이 2회선으로 연결되어 있으며 각 모선은 각각 부하모선에 연결되어 있는 계통의 단선도 이다.The following describes the contents of the present invention described above by way of examples. 9 is a single circuit diagram of a system in which six generators are connected to the same bus by two, buses and buses are connected in two lines, and each bus is connected to a load bus.

발전기 모선(900, 901)과 발전기 모선에서 비교적 원거리(902)에 고장을 모의하여 각 발전기들의 응동과 고장에 따른 발전기들의 누적된 임계에너지를 계산하여 도10에 나타내었다. 도10에서 발전기 G₁~G₄와 G₅~G₆는 서로 다른 응동 특성을 나타내고 있음을 알 수 있다.The faults are simulated at a relatively long distance 902 in the generator buses 900 and 901 and the generator buses, and the accumulated critical energy of the generators according to the response and failure of each generator is calculated and shown in FIG. 10. In FIG. 10, it can be seen that the generators G ₁ to G ₄ and G ₅ to G ₆ exhibit different response characteristics.

발전기들이 동일 모선에 연결(G₁-G₂, G₃-G₄, G₅-G₆)된 경우, 두 발전기는 거의 동일한 과도특성을 나타낸다. 도10에서 모든 고장에 대해서 동일 모선에 연결된 에너지는 같음을 알 수 있다.When the generators are connected to the same bus (G ₁ -G ₂ , G ₃ -G ₄ , G ₅ -G ₆ ), the two generators exhibit almost identical transients. It can be seen from FIG. 10 that all the faults have the same energy connected to the same bus.

가장 심각한 경우는 선로1-2 및 선로3-9에 고장이 발생한 경우로 발전기 모선 근처의 고장이므로 발전기 단자전압이 영에 가깝다. 이때 임계가속에너지는 1.287~1.362 사이에 형성되고 이로부터 과도안정도는 한 발전기 또는 동일모선에 연결된 발전기들의 응동이 좌우함을 알 수 있다. 따라서 ε_G1을 G₁ 발전기의 설정값이라 할 때, 각 발전기 임계가속에너지의 설정 값은 다음과 같이 설정할 수 있다.The most serious cases are the failures of lines 1-2 and 3-9, which are near the generator bus, so the generator terminal voltage is close to zero. At this time, the critical acceleration energy is formed between 1.287 and 1.362, and it can be seen from this that the transient stability depends on the response of generators connected to one generator or the same bus. Therefore, when ε _G1 is the set value of the G ₁ generator, the set value of the critical acceleration energy of each generator can be set as follows.

ε_G1=1.287, ε_G2=1.300, ε_G3=1.350, ε_G4=1.360, ε_G5=1.362, ε_G6=1.362ε _G1 = 1.287, ε _G2 = 1.300, ε _G3 = 1.350, ε _G4 = 1.360, ε _G5 = 1.362, ε _G6 = 1.362

또한 설정값의 총합은 발전기 그룹의 임계가속에너지 설정값 ε_gs이므로 다음과 같다.In addition, the sum of the set values is the critical acceleration energy set value ε _gs of the generator group as follows.

ε_gs=8.021ε _gs = 8.021

실측데이터에서 임계가속에너지는 도10에 나타난 설정 값을 계산할 때와 동일한 방법으로 계산하여 얻을 수 있는데, 과도안정도는 이들 설정 값과 실측데이터 에서 계산한 임계가속에너지를 비교함으로써 판단할 수 있다. 만일 G₁ 발전기에서 실측한 데이터에서 계산한 임계가속에너지 E_s가 1.287보다 크면, G₁ 발전기는 과도불안정 하다고 판단할 수 있다. 다른 발전기들에 대해서도 동일한 관계가 성립한다.In the measured data, the critical acceleration energy can be obtained by calculating the same method as when calculating the set values shown in FIG. 10. The transient stability can be determined by comparing these set values with the critical acceleration energy calculated from the measured data. If one threshold acceleration energy E _s calculated from the data measured at the generator G ₁ is larger than 1.287, the generator G ₁ may determine that the transient instability. The same relationship holds for other generators.

전체 발전기 그룹 임계가속에너지는 선로1-2 고장일 때는 6.341이고 선로3-9 고장일 때는 3.685이다. 따라서 그룹 임계가속에너지로부터 그룹의 과도안정도를 판별하는 것은 어렵고 또한 한 발전기의 과도불안정특성이 전체 그룹에 영향을 주기 때문에 그룹의 과도안정도를 판별하는 것은 불필요하다. 그러나 그룹 임계가속에너지로부터 발전기 탈락에 대한 여유를 계산할 수 있어 그룹의 안정도 여유를 판단할 수 있는 기준이 된다. The critical generator energy for the entire generator group is 6.341 for line 1-2 failure and 3.685 for line 3-9 failure. Therefore, it is difficult to determine the transient stability of the group from the group critical acceleration energy, and it is unnecessary to determine the transient stability of the group because the transient instability of one generator affects the entire group. However, the margin for generator dropout can be calculated from the group critical acceleration energy, which is a criterion for judging the stability margin of the group.

만일 선로 1-2 고장에서 발전기 G₁을 탈락하면, 1.287/6.341=0.203으로 그룹의 가속에너지는 20.3%의 여유가 생기고 G₁-G₂를 동시에 탈락하면, 2.587/6.341=0.408로 40.8%의 그룹 가속에너지의 여유가 발생한다. 따라서 선로 1-2 고장에서는 발전기 G₁-G₂ 또는 G₃-G₄를 동시에 탈락하면 다른 발전기의 과도안정도를 더욱 크게 향상시킬 수 있다. If the generator G ₁ is dropped in line 1-2 failure, the acceleration energy of the group will be 20.3% at 1.287 / 6.341 = 0.203, and if G ₁ -G ₂ is dropped at the same time, it will be 40.8% at 2.587 / 6.341 = 0.408. There is a margin of group acceleration energy. Therefore, in line 1-2 failure, dropping generators G ₁ -G ₂ or G ₃ -G ₄ simultaneously can greatly improve the transient stability of other generators.

또한 선로 3-9 고장에서 발전기 G₅를 탈락하면, 1.362/3.685=0.370으로 37.0%의 그룹 가속에너지 여유가 발생하고, G₅-G₆를 동시에 탈락하면, 2.724/3.685=0.740으로 74.0%의 그룹 가속에너지의 여유가 발생한다. 그러므로 선 로 3-9 고장으로 발전기 G₅-G₆가 과도불안정도 상태가 되면, 동시에 탈락하는 것이 나머지 발전기 G₁-G₄의 과도안정도에 향상에 중요함을 알 수 있다.Also, if the generator G ₅ is dropped in line 3-9 failure, a group acceleration energy margin of 37.0% occurs at 1.362 / 3.685 = 0.370, and if G ₅ -G ₆ is dropped at the same time, 74.0% at 2.724 / 3.685 = 0.740. There is a margin of group acceleration energy. Therefore, if the generator G ₅ -G ₆ becomes transient instability due to the failure of line 3-9, it can be seen that dropping out at the same time is important for improving the transient stability of the remaining generators G ₁ -G ₄ .

발전기 모선이 아닌 선로에서는 발전기 모선 고장에서의 가속에너지 보다 큰 임계가속에너지를 갖고 있다. 그러나 임계가속에너지는 큰 차이가 없음을 알 수 있다. 예를 들면 선로 1-2의 고장에서 G₁=1.287, G₂=1.300이지만 선로 5-7의 고장에서 G₁=1.398, G₂=1.406이다. 발전기모선의 고장이 더 작은 임계가속에너지를 가지고 있음을 알 수 있고 따라서 발전기 모선 고장에 대한 임계가속에너지를 설정치로 사용하는 것이 적절함을 알 수 있다. In non-generator bus lines, the critical acceleration energy is larger than the acceleration energy in the fault of the generator bus. However, it can be seen that the critical acceleration energy does not have a big difference. For example, G ₁ = 1.287, G ₂ = 1.300 at line 1-2, but G ₁ = 1.398, G ₂ = 1.406 at line 5-7. It can be seen that the failure of the generator bus has a smaller critical acceleration energy, and therefore it is appropriate to use the critical acceleration energy for the generator bus failure as a set point.

도11은 선로 3-9고장에서 발전기 G₁의 전력과 기계적입력의 시간응답을 나타내고 있고, 도12는 발전기 G₅의 시간응답 나타내고 있다.Fig. 11 shows the time response of the power and mechanical input of the generator G ₁ in the line 3-9 failure, and Fig. 12 shows the time response of the generator G ₅ .

각 발전기에 대한 임계가속에너지와 그룹에 대한 임계가속에너지의 조합이 다양하고 네트워크 구조가 모든 발전기 그룹에서 각각 다르므로 발전기 탈락에 대한 다양한 시나리오가 가능하다. 본 발명에서는 각 발전기의 임계가속에너지로 탈락여부를 결정하고, 그룹에서 탈락시켜야할 발전기는 임계가속에너지가 가장 큰 발전기를 선택한다.Various combinations of critical acceleration energy for each generator and critical acceleration energy for the group and the network structure are different in every generator group are possible. In the present invention, it is determined whether to drop out by the critical acceleration energy of each generator, the generator to be dropped in the group select the generator with the largest critical acceleration energy.

[실시예2]Example 2

이산치 신호가 연속적으로 입력될 때, 입력된 신호들의 최소의 조합으로 시스템의 응답이 정상상태인지의 여부를 판단해야한다. 본 발명에서는 고장 여부를 판단하는데 있어 주어진 시계열데이터에서 최소의 데이터수와 계산량으로 고장을 판단하는 방법을 위하여 표준편차 s_td를 도입하였다.When discrete signals are continuously input, it is necessary to determine whether the response of the system is normal with the minimum combination of the input signals. In the present invention, a standard deviation s _td is introduced for a method of determining a failure with a minimum number of data and a calculation amount in a given time series data.

수학식 17에서 N은 데이터 수이고, m은 N개로 샘플링한 전력의 평균이며, P_e[n]는 n-번째 데이터의 전력을 의미한다. 데이터 수 N은 샘플링 구간에 따라 임의로 설정할 수 있으나 0.01초의 샘플링 구간에서 N=10개 이상으로 선택하면 고장 구별이 가능하다. 상기 수학식17에서는 n-번째 전력데이터에 대한 표준편차를 정의하였으나 이를 가속력 P_a[n]과 그 평균을 이용하여 가속력에 대한 표준편차를 이용해도 동일한 결과를 얻을 수 있다.In Equation 17, N is the number of data, m is the average of the power sampled by N, P _e [n] means the power of the n-th data. The number of data N can be set arbitrarily according to the sampling interval, but failure can be distinguished by selecting N = 10 or more in the sampling interval of 0.01 second. In Equation 17, although the standard deviation of the n-th power data is defined, the same result can be obtained by using the standard deviation of the acceleration force using the acceleration force P _a [n] and the average thereof.

도13에는 실측데이터에서 과도상태 전후의 데이터를 나타낸 것으로 std1은 n=k~k+9까지 10개의 데이터에 대한 표준편차이고, std2는 n=k+1~k+10까지 10개의 데이터에 대한 표준편차이다. 만일 k+10번째 데이터가 과도상태가 직후 데이터라면 std1은 과도상태 직전의 표준편차를 의미하고, std2는 과도상태 직후 1개의 데이터가 포함된 표준편차를 의미한다. 과도상태 직전에 계산한 5개의 표준편차는 다음과 같다. Figure 13 shows the data before and after the transient state in the measured data, std1 is the standard deviation for 10 data from n = k ~ k + 9, std2 is for 10 data from n = k + 1 ~ k + 10 Standard deviation. If the k + 10th data is immediately after the transient state, std1 means the standard deviation immediately before the transient state, and std2 means the standard deviation including one data immediately after the transient state. The five standard deviations calculated immediately before the transient are:

no1, std1= 0.0260 std2=0.0213 no1, std1 = 0.0260 std2 = 0.0213

no2, std1= 0.0213 std2=0.0211 no2, std1 = 0.0213 std2 = 0.0211

no3, std1= 0.0211 std2=0.0209 no3, std1 = 0.0211 std2 = 0.0209

no4, std1= 0.0209 std2=0.0235 no4, std1 = 0.0209 std2 = 0.0235

no5, std1= 0.0255 std2=1.3969 no5, std1 = 0.0255 std2 = 1.3969

과도상태 직후 1개의 데이터가 포함된 데이터(no5)의 표준편차가 std2=1.3969로 다른 표준편차와 비교하면 매우 큰 값을 가지고 있음을 알 수 있다. 이 경우는 발전기 모선의 심각한 고장이 아닌 경우인데도 그 선별성이 매우 뛰어남을 알 수 있다. Immediately after the transient state, the standard deviation of the data (no5) containing one data is std2 = 1.3969, which is very large compared to other standard deviations. In this case, even if it is not a serious failure of the generator busbar, it can be seen that the selectivity is excellent.

이상에서와 같이 본 발명은 특정 실시예를 예시하여 설명하지만 본 발명이 상기 실시예에 한정되는 것은 아니며, 본 발명의 범주를 벗어나지 않는 범위내에서 여러 가지로 변형된 실시예도 가능하다.As described above, the present invention is illustrated by exemplifying specific embodiments, but the present invention is not limited to the above embodiments, and various modified embodiments may be made without departing from the scope of the present invention.

상술한 바와 같이 본 발명은, 전력시스템에서 실시간으로 과도안정도를 판별하고 고장이 파급되는 것을 방지하는 방법에 관한 것으로, 본 발명에서는 실시간으로 실측한 전력으로부터 데이터가 입력된 사이클 마다 과도안정도를 판별할 수 있어, 다른 발전기로 심각한 고장이 파급되는 것을 방지할 수 있어 보다 안정적으로 전력을 생산, 공급할 수 있다. 또한, 불필요한 발전기 정지 피할 수 있고, 대규모 정전 등을 예방할 수 있어 막대한 전력시스템 운용 비용을 절감할 수 있다.As described above, the present invention relates to a method for determining a transient stability in real time in a power system and preventing a failure from spreading. In the present invention, the transient stability can be determined for each cycle in which data is input from a measured power in real time. It is possible to prevent the spread of serious failures to other generators, thereby producing and supplying electricity more stably. In addition, unnecessary generator stoppages can be avoided, and large-scale power outages can be prevented, thereby reducing enormous power system operation costs.

Claims

실시간으로 전력시스템의 과도안정도를 판별하여 전체 전력시스템으로 고장이 파급되는 것을 방지하는데 있어서,In determining the transient stability of the power system in real time to prevent the failure from spreading to the entire power system,

실시간으로 발전기에서 출력되는 전압, 전류와 터빈에서 출력되어 발전기로 입력되는 기계적 입력을 취득하는 스텝과,A step of acquiring a voltage, a current output from the generator in real time, and a mechanical input output from the turbine and input to the generator;

상기 취득된 전압과 전류로부터 전력(P_e)을 실시간으로 계산하여 고장을 판단하는 하는 스텝과,Calculating a power _Pe from the acquired voltage and current in real time to determine a failure;

고장이라 판단되면 기계적 입력(P_m)과 상기에서 계산한 전력(P_e)의 편차로부터 가속력(P_a=P_m-P_e)을 계산하는 스텝과,If it is determined that the failure is a step of calculating the acceleration force (P _a = P _m -P _e ) from the deviation between the mechanical input (P _m ) and the electric power (P _e ) calculated above,

상기 계산된 가속력을 고장이 발생한 시간부터 적분하여 가속에너지(E_a)를 계산하는 스텝과,Calculating the acceleration energy E _a by integrating the calculated acceleration force from the time when the failure occurs;

상기 산출된 가속에너지(E_a)와 설정값을 비교하여 가속에너지가 설정값(ε) 보다 크면 과도불안정으로 판단하여 발전기를 탈락(trip)하는 스텝과,Comparing the calculated acceleration energy (E _a ) with the set value, and if the acceleration energy is greater than the set value (ε), determining that it is transient instability and tripping the generator;

상기 가속에너지와 설정값 및 탈락발전기 번호를 출력 및 저장하는 스텝을 포함하는 것을 특징으로 하는 전력시스템에서 실시간 과도안정도 판별 방법을 이용한 고장파급방지법.And a step of outputting and storing the acceleration energy, the set value, and the dropping generator number.

제 1항에 있어서,The method of claim 1,

설정값(ε)은 전력시스템의 과도안정도 해석 프로그램을 사용하여 동일한 전력시스템을 모의하는 스텝과, The set value (ε) is a step of simulating the same power system using the transient stability analysis program of the power system,

상기에서 모의 된 전력시스템의 각 모선에 고장을 발생시켜 발전기가 탈조 되기 직전시간인 임계고장제거시간(t_c)에서 각 모선 고장에 대응하는 가속에너지(E_a)를 계산하는 스텝과 Calculating the acceleration energy E _a corresponding to each bus failure at the critical failure elimination time t _c , which is the time immediately before the generator is disassembled by generating a failure in each bus of the simulated power system;

상기 계산된 가속에너지들 중 최소값을 설정값(ε)으로 치환하는 스텝을 포함하는 것을 특징으로 하는 전력시스템에서 실시간 과도안정도 판별 방법을 이용한 고장파급방지법.And a step of replacing the minimum value of the calculated acceleration energies with a set value (ε).

전력시스템의 과도불안정시 고장이 파급되는 것을 방지하기 위하여 탈락시킬 발전기를 선택하는데 있어서,In selecting a generator to drop out in order to prevent the failure from spreading in case of excessive instability of the power system,

모선에 연결된 각 발전기들의 가속에너지(E_ai)를 합산하여 동일 지역 발전기 그룹의 총 가속에너지(E_s)를 계산하는 스텝과,Calculating the total acceleration energy (E _s ) of the same area generator group by summing the acceleration energy (E _ai ) of each generator connected to the bus line;

상기 계산된 총 가속에너지에서 각 발전기의 가속에너지를 감산하여 각 발전기에 대한 과도안정도 여유(E_res=E_s-E_ai)를 계산하는 스텝과,Calculating a transient stability margin (E _res = E _s -E _ai ) for each generator by subtracting the acceleration energy of each generator from the calculated total acceleration energy;

상기 계산된 과도안정도 여유(E_res)가 가장 작은 발전기를 탈락(trip) 할 발전기로 선택하는 스텝을 포함하는 것을 특징으로 하는 전력시스템에서 실시간 과도안정도 판별 방법을 이용한 고장파급방지법.And a step of selecting the generator having the smallest calculated transient stability E _res as the generator to trip the fault, using the real-time transient stability determination method in the power system.

실시간으로 임계고장제거시간(t_c)을 계산하여 과도불안정을 판단하는데 있어서,In determining the transient instability (t _c ) in real time to determine the transient instability,

현재 시간(t₂)에서 계산한 가속에너지(E₂)와 이전 시간(t₁)에서 계산한 가속에너지(E₁)의 편차(E₂-E₁)로 시간편차(t₂-t₁)를 나눈 후 설정값(ε)을 곱하여 임계고장제거시간(t_c)을 계산하는 스텝과The current time (t ₂₎ an acceleration energy (E ₂₎ and the previous time (t ₁₎ an acceleration energy deviation (E ₂ -E ₁₎ time difference (t ₂ -t ₁₎ as the (E ₁₎ calculated in the calculation in Divide by and multiply the set value (ε) to calculate the critical fault removal time (t _c )

상기에서 계산된 임계고장제거시간(t_c)이 설정된 임계고장제거시간보다 크면 과도불안정으로 판단하는 스텝을 포함하는 것을 특징으로 하는 전력시스템에서 실시간 과도안정도 판별 방법을 이용한 고장파급방지법.And determining the transient instability if the threshold failure removal time t _c calculated above is greater than the set threshold failure removal time.

제 1항에서, In claim 1,

가속에너지(E_a)를 가속력의 누적 합인 누적가속력(P_s)으로 치환하여 과도안정도를 판별하는 스텝을 포함하는 것을 특징으로 하는 전력시스템에서 실시간 과도안정도 판별 방법을 이용한 고장파급방지법.And a step of determining the transient stability by substituting the acceleration energy (E _a ) with the cumulative acceleration force (P _s ), which is the cumulative sum of the acceleration forces, using the real-time transient stability determination method.

삭제delete

제 1항에 있어서, The method of claim 1,

실시간으로 전력을 계산하여 고장을 판단하는 스텝은, 연속적으로 계산된 10개 전력 데이터의 표준편차를 계산하는 스텝과,The step of determining the failure by calculating the power in real time comprises the steps of calculating the standard deviation of the ten consecutively calculated power data,

상기에서 계산한 표준편차가 이전 시간에서 계산한 표준편차보다 기준값 이상으로 크면 고장이라 판단하는 스텝을 포함하는 것을 특징으로 하는 전력시스템에서 실시간 과도안정도 판별 방법을 이용한 고장파급방지법.If the standard deviation calculated above is greater than the standard deviation calculated from the previous time is greater than the reference value failure failure prevention method using a real-time transient stability determination method comprising the step of determining that the failure.