JPH06121455A - Protective relay - Google Patents

Protective relay

Info

Publication number
JPH06121455A
JPH06121455A JP26331092A JP26331092A JPH06121455A JP H06121455 A JPH06121455 A JP H06121455A JP 26331092 A JP26331092 A JP 26331092A JP 26331092 A JP26331092 A JP 26331092A JP H06121455 A JPH06121455 A JP H06121455A
Authority
JP
Japan
Prior art keywords
vector
relay
current
calculated
impedance
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.)
Granted
Application number
JP26331092A
Other languages
Japanese (ja)
Other versions
JP3003421B2 (en
Inventor
Toshinobu Ebizaka
敏信 海老坂
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Priority to JP4263310A priority Critical patent/JP3003421B2/en
Publication of JPH06121455A publication Critical patent/JPH06121455A/en
Application granted granted Critical
Publication of JP3003421B2 publication Critical patent/JP3003421B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Emergency Protection Circuit Devices (AREA)

Abstract

PURPOSE:To perform interruption operation more accurately by calculating high frequency components generated in a power system by the occurrence of a fault. CONSTITUTION:An impedance vector is calculated by dividing a voltage vector by a current vector, or an operational vector is calculated by dividing the current vector by a fundamental wave reference vector in an arithmetic processing part 8, on the basis of various kinds of electric data sampled from a power system 12. And, the maximum and minimum values of this and their time are extracted and stored in a RAM 10. On the basis of these, the magnitude of the vector locus and the speed of rotation are calculated at the processing part 8. Besides, the ratios of high-frequency components contained in the vectors and their orders are calculated. Consequently, it becomes possible to perform interruption operation more accurately.

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【産業上の利用分野】本発明は、電力系統における事故
を検出して、該当する事故区間の除去指令を出力する保
護継電器に関し、特に電力系統からサンプリングした各
種電気量データをディジタル化して演算し、この結果に
応じて動作判定を行うディジタル演算形の保護継電器に
関するものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a protective relay which detects an accident in a power system and outputs a command to eliminate the corresponding accident section, and in particular, digitizes various electric quantity data sampled from the power system for calculation. The present invention relates to a digital operation type protective relay that makes an operation determination according to the result.

【0002】[0002]

【従来の技術】電力系統に短絡または接地などの事故が
発生した場合、あるいは他の系統に対して悪影響を与え
るような異常運転が行われた場合に、このような異常状
態を検出してその部分を系統から速やかに切り離すため
に、電力用保護継電器が電力系統の構成に応じて各所に
配置されている。この保護継電器には、その動作原理に
応じて種々の継電器が存在するが、その1つにディジタ
ル演算形の保護継電器(以下、ディジタルリレーとい
う)がある。これは、各種の系統電気量をディジタルデ
ータとしてサンプリングした後、所定の演算手順により
演算処理し、その結果に基づきリレーの動作判定を行う
ものであり、その技術的原理は、「ディジタルリレー実
務読本」(オーム社/H3.8.25初版発行)等で解
説されている。
2. Description of the Related Art When an accident such as a short circuit or grounding occurs in a power system, or when an abnormal operation that adversely affects other systems is performed, such an abnormal state is detected and detected. In order to quickly disconnect the parts from the grid, protective relays for electric power are arranged in various places according to the configuration of the power grid. There are various relays in this protective relay according to the operating principle thereof, and one of them is a digital operation type protective relay (hereinafter referred to as a digital relay). This is to sample various system electric quantities as digital data, to perform arithmetic processing according to a predetermined arithmetic procedure, and to judge the operation of the relay based on the result, the technical principle of which is "Digital Relay Practical Manual". "(Ohm Co./H3.8.25 first edition issued) etc.

【0003】さて、このようなディジタルリレーは、実
際の電力系統において特定区間の系統を保護する場合、
図6に示すように配置されている。同図に用いられてい
るリレーは、系統から得られる電圧と電流の比、すなわ
ちインピーダンスが所定値内となったとき動作するよう
に構成されており、系統におけるインピーダンスが送電
線の距離を表すことから、距離リレーとも呼ばれてい
る。従って、このリレーにより、系統に発生した事故の
発生方向と事故地点までの距離を判定することができ
る。
Now, such a digital relay protects a system of a specific section in an actual power system,
It is arranged as shown in FIG. The relay used in the figure is configured to operate when the ratio of the voltage and current obtained from the grid, that is, the impedance, is within a predetermined value, and the impedance in the grid represents the distance of the transmission line. Therefore, it is also called a distance relay. Therefore, with this relay, it is possible to determine the direction of the accident that has occurred in the system and the distance to the accident point.

【0004】図6において、61は電源を示しており、
62はリレー64に系統の電圧を導入するコンデンサ分
圧形電圧変成器、また63はリレー64に系統の電流を
導入する変流器である。リレー64はディジタル演算形
の系統保護リレー(DZ)である。65は線路側系統に
おけるリアクタンス分L、66は同じくコンデンサ分C
を示している。以上の構成により、変成器62からの電
圧および変流器63からの電流に応じて、これら電圧と
電流の比、すなわちリレー64から見た系統のインピー
ダンスが算出され、この算出結果に基づく動作判定によ
り、変流器63以降の線路側系統(保護区間)が保護さ
れている。
In FIG. 6, reference numeral 61 denotes a power source,
Reference numeral 62 is a capacitor voltage dividing type voltage transformer for introducing a system voltage to the relay 64, and 63 is a current transformer for introducing a system current to the relay 64. The relay 64 is a digital operation type system protection relay (DZ). Reference numeral 65 is a reactance component L in the line system, and 66 is a capacitor component C as well.
Is shown. With the above configuration, the ratio of the voltage and the current, that is, the impedance of the system viewed from the relay 64 is calculated according to the voltage from the transformer 62 and the current from the current transformer 63, and the operation determination based on the calculation result is performed. Thus, the line side system (protection section) after the current transformer 63 is protected.

【0005】このような一般的な電力系統において、近
年、電圧階級が上昇し、またケーブル系統が増大する傾
向にあり、この影響で線路側系統のリアクタンス分Lや
コンデンサ分Cが増大し、事故発生時に電源の基本周波
数fに対する1.5倍程度の周波数の低次高調波を発生
する系統が増加している。図6において、このリレー6
4により保護される線路側系統より背後のFa地点で事
故が発生し、線路側系統のリアクタンス分65およびコ
ンデンサ分66により、線路側からFa地点に向かって
nf分の電流Inf、また電源61からFa地点に向かっ
て電流I1fがそれぞれ発生したとする。なお、nf分と
は電源61の基本周波数f(=1f)に対するn次の高
調波分を示している。
In such a general electric power system, in recent years, the voltage class has risen and the cable system has tended to increase. Due to this influence, the reactance component L and the capacitor component C of the line side system increase and an accident occurs. An increasing number of systems generate low-order harmonics having a frequency of about 1.5 times the fundamental frequency f of the power source when they occur. In FIG. 6, this relay 6
An accident occurs at a point Fa behind the line side system protected by 4, and the reactance 65 and the capacitor 66 of the line side system cause a current Inf of nf toward the point Fa from the line side, and a power source 61. It is assumed that the current I 1f is generated from each point toward the point Fa. The nf component indicates the nth harmonic component with respect to the fundamental frequency f (= 1f) of the power supply 61.

【0006】ここで、前述のような系統状況から、電流
nfには1.5f分が多く含まれるため、これが電流と
して変流器63を介してリレー64へ導入され、一方、
この1.5f分に比べて電源61の影響を受け1f分を
多く含む電圧が変成器62を介してリレー64へ導入さ
れる。これに応じてリレー64では、前述のとおりイン
ピーダンスが算出され、動作判定がなされる。通常であ
れば、電流Infの方向が、保護されるべき線路側系統で
の事故発生時とは逆の方向、すなわち電源61へ向かっ
て流れることから、算出されたインピーダンスは所定値
以外となり、事故発生地点が保護区間外であると判断さ
れてリレー64は動作しない。
[0006] Here, from the system status as described above, the current I nf because it contains a lot 1.5f component, which is introduced to the relay 64 via the current transformer 63 as a current, whereas,
A voltage including a large amount of 1f, which is affected by the power supply 61 as compared with the 1.5f, is introduced into the relay 64 via the transformer 62. In response to this, in the relay 64, the impedance is calculated and the operation is determined as described above. Normally, the direction of the current Inf flows in the direction opposite to that at the time of the occurrence of an accident in the line side system to be protected, that is, toward the power supply 61, so the calculated impedance becomes a value other than the predetermined value, The relay 64 does not operate because it is determined that the accident occurrence point is outside the protection section.

【0007】しかし、電圧と電流の周波数に前述のよう
なズレがあるため、このインピーダンスベクトルが時間
とともに移動して、図9に示すようにリレーの本来動作
すべきベクトル領域へ進入し、リレー64を動作させ
る。すなわち、保護区間外である背後の系統で発生した
事故に対しても、リレー64を動作させてしまうことに
なる。従来、このような誤動作を防止するため、ディジ
タルリレーにおいて、図8に示されたようなフィルター
特性を持つディジタルフィルターを設けて、系統からサ
ンプリングした各種データから1.5f相当分の低次高
調波分を抽出し、別途抽出した1f分と比較して一定値
以上検出された場合には、リレー64の動作を強制的に
抑止する手段を設けていた。
However, since the frequency of the voltage and the current have the above-mentioned deviation, this impedance vector moves with time and enters the vector region where the relay should normally operate, as shown in FIG. To operate. That is, the relay 64 is operated even in the case of an accident that occurs in the back system outside the protection section. Conventionally, in order to prevent such malfunction, a digital filter having a filter characteristic as shown in FIG. 8 is provided in a digital relay, and various low-order harmonics equivalent to 1.5f are obtained from various data sampled from the system. A means for forcibly suppressing the operation of the relay 64 is provided when a minute is extracted and compared with a separately extracted 1f minute and a predetermined value or more is detected.

【0008】また、ディジタルリレーはこのような系統
保護リレーとしての配置される他に、図7に示すように
変圧器を保護する場合にも用いられる。同図に用いられ
るリレーは、保護する変圧器の両端の電流の差が一定値
以上になったとき動作するように構成されており、差動
リレーとも呼ばれている。図7において、71は電源を
示しており、72および73はリレー74に変成器77
の両端の電流を導入する変流器である。リレー74はデ
ィジタル演算形の変圧器保護リレー(TP)であり、変
圧器77は系統に直列に配置され、このリレー74によ
って保護される。75は線路側系統におけるリアクタン
ス分L、76は同じくコンデンサ分Cを示す。以上の構
成により、変流器72および73からの電流に応じて、
変圧器77への流入電流の和、すなわち差動電流が算出
され、この算出結果に基づく動作判定により、変圧器7
7が保護されている。
Further, the digital relay is used not only as such a system protection relay but also for protecting a transformer as shown in FIG. The relay used in the figure is configured to operate when the difference between the currents across the transformer to be protected exceeds a certain value, and is also called a differential relay. In FIG. 7, 71 indicates a power source, and 72 and 73 are relays 74 and transformers 77.
It is a current transformer that introduces current across both ends of. The relay 74 is a digital operation type transformer protection relay (TP), and the transformer 77 is arranged in series in the system and is protected by this relay 74. Reference numeral 75 represents a reactance component L in the line side system, and 76 represents a capacitor component C. With the above configuration, according to the currents from the current transformers 72 and 73,
The sum of the currents flowing into the transformer 77, that is, the differential current is calculated, and the transformer 7 is determined by the operation determination based on the calculation result.
7 are protected.

【0009】図7において、この変圧器保護リレー74
による保護区間内、すなわち変圧器77またはこの近傍
のFb地点で事故が発生し、線路側系統のリアクタンス
分75およびコンデンサ分76により、線路側系統から
Fb地点に向かってnf分の電流Infが、また電源71
からFb地点に向かって電流I1fがそれぞれ発生したと
する。ここで、前述のような系統状況から、電流Inf
1.5f分を多く含み、これが変流器73を介してリレ
ー74へ導入され、一方、電源71からの1f分の電流
1fが変流器72を介してリレー74へ導入される。こ
れに応じてリレー74では、前述のとおり差動電流が算
出され、動作判定がなされる。通常であれば電流Inf
増加することにより、算出された差動電流は所定値以上
となり、リレー74は動作する。
In FIG. 7, this transformer protection relay 74
An accident occurs in the protection section due to, that is, the transformer 77 or the Fb point in the vicinity thereof, and the reactance 75 and the capacitor 76 of the line side system cause a current Inf of nf amount from the line side system toward the Fb point. , Power supply 71 again
It is assumed that the currents I 1f are respectively generated from F to the point Fb. Here, from the system status as described above, the current I nf is rich in 1.5f component, which is introduced to the relay 74 via the current transformer 73, whereas, the 1f component of the current I 1f from the power source 71 It is introduced into the relay 74 via the current transformer 72. In response to this, the relay 74 calculates the differential current as described above and determines the operation. Normally, the current Inf increases, so that the calculated differential current becomes a predetermined value or more, and the relay 74 operates.

【0010】しかし、変圧器は電圧印加時の鉄心飽和に
より2f分が多く発生するという特性を持っており、従
来のディジタルリレーにおいては、このような変圧器に
特有な励磁突流現象の影響を除去するために、図8に示
すようなフィルター特性を持つディジタルフィルターを
設けて、系統からサンプリングした電流データから2f
相当分の低次高調波を抽出し、別途抽出した1f分との
比I2f/I1fの値が一定値以上である場合には、リレー
74の動作を強制的に抑止する手段を設けていた。
However, the transformer has a characteristic that a large amount of 2f is generated due to iron core saturation when a voltage is applied, and in the conventional digital relay, the influence of such an exciting surge current characteristic of the transformer is eliminated. In order to do so, a digital filter having the filter characteristics as shown in FIG. 8 is provided and 2f is obtained from the current data sampled from the system.
A means for forcibly suppressing the operation of the relay 74 is provided when a considerable amount of low-order harmonics is extracted and the ratio I 2f / I 1f with the separately extracted 1f is equal to or greater than a certain value. It was

【0011】[0011]

【発明が解決しようとする課題】従って、前者のような
従来の系統保護リレーにおいては、1.5f分を多く含
む低次高調波の検出に応じて、本来の保護動作を強制的
に抑止しているため、前述のような状況にある系統では
低次高調波に対する演算精度に欠け、事故地点の特定、
すなわち正確な電力系統の保護が困難となるという問題
があった。また、後者のような従来の変圧器保護リレー
においては、前述のような系統に変圧器が配置されてい
る場合、事故時に発生する1.5f分の低次高調波が2
f分として検出されて、リレーの動作が強制的に抑止さ
れるため、事故発生にもかかわらず保護動作が行われな
いという問題があった。本発明はこのような課題を解決
するためのものであり、電力系統からリレーに取り込ま
れた各種データに含まれる高調波の次数とその大きさを
算出することにより、遮断動作の的確な判定を可能とす
る各種保護継電器を提供することを目的としている。
Therefore, in the conventional system protection relay such as the former, the original protection operation is forcibly suppressed in response to the detection of low-order harmonics containing a large amount of 1.5f. Therefore, in the system in the above situation, the calculation accuracy for low-order harmonics is lacking, and the location of the accident point can be identified.
That is, there is a problem that it is difficult to accurately protect the power system. In the latter type of conventional transformer protection relay, when the transformer is arranged in the system as described above, 1.5 f of low-order harmonics generated at the time of accident occurs
Since it is detected as f minutes and the operation of the relay is forcibly suppressed, there is a problem that the protection operation is not performed despite the occurrence of the accident. The present invention is to solve such a problem, by calculating the order of harmonics and its magnitude contained in various data taken into the relay from the power system, the accurate determination of the interrupting operation. The purpose is to provide various protective relays that can be used.

【0012】[0012]

【課題を解決するための手段】このような目的を達成す
るために、本発明による保護継電器は、電力系統からサ
ンプリングされた電圧および電流データに基づき、電圧
ベクトルを電流ベクトルで除算したインピーダンスベク
トルを算出するベクトル演算手段と、前記インピーダン
スベクトルの最大値および最小値ならびにその時刻をそ
れぞれ抽出する抽出手段と、この抽出結果に基づき前記
インピーダンスベクトルに含まれる高調波分の割合とそ
の次数を算出する高調波分算出手段とを備えるものであ
る。また、電力系統からサンプリングされた電流データ
に基づき、電流ベクトルを電力系統の基本波基準ベクト
ルで除算した演算ベクトルを算出するベクトル演算手段
と、前記演算ベクトルの最大値および最小値ならびにそ
の時刻をそれぞれ抽出する抽出手段と、この抽出結果に
基づき前記演算ベクトルに含まれる高調波分の割合とそ
の次数を算出する高調波分算出手段とを備えるものであ
る。
In order to achieve such an object, a protective relay according to the present invention provides an impedance vector obtained by dividing a voltage vector by a current vector based on voltage and current data sampled from a power system. Vector calculating means for calculating, extraction means for extracting the maximum value and the minimum value of the impedance vector and the time thereof, respectively, and a harmonic for calculating the proportion of harmonics contained in the impedance vector and its order based on the extraction result. And a wavelet calculation means. Further, based on the current data sampled from the electric power system, vector operation means for calculating an operation vector by dividing the current vector by the fundamental wave reference vector of the power system, the maximum value and the minimum value of the operation vector and its time respectively. It is provided with an extracting means for extracting, and a harmonic component calculating means for calculating the ratio and the order of the harmonic components contained in the operation vector based on the extraction result.

【0013】[0013]

【作用】従って、電力系統状態を表す各種サンプリング
データに基づき、インピーダンスベクトルの軌跡がベク
トル演算手段により算出され、このベクトルの大きさと
回転速度からベクトルに含まれる高調波分の次数とその
大きさが演算されて、リレー動作の判定がなされる。ま
た、電力系統状態を表す各種サンプリングデータに基づ
き、抽出された差動電流ベクトルを電源基本波ベクトル
で除した演算ベクトルの軌跡がベクトル演算手段により
算出され、このベクトルの大きさと回転速度からベクト
ルに含まれる高調波分の次数とその大きさが演算され
て、リレー動作の判定がなされる。
Therefore, the locus of the impedance vector is calculated by the vector calculating means on the basis of various sampling data representing the state of the electric power system, and the order of the harmonic component contained in the vector and its magnitude are calculated from the magnitude and rotational speed of the vector. It is calculated and the relay operation is determined. Further, based on various sampling data representing the state of the power system, the vector operation means calculates the locus of the operation vector obtained by dividing the extracted differential current vector by the power source fundamental wave vector, and the vector is calculated from the size and rotation speed of this vector. The order of the contained harmonics and its magnitude are calculated to determine the relay operation.

【0014】[0014]

【実施例】次に、本発明について図面を参照して説明す
る。図1は本発明の一実施例であるディジタルリレーの
ブロック図である。図1において、1はディジタルリレ
ーであり、アナログ/ディジタル変換部2とディジタル
演算制御部7から構成されている。アナログ/ディジタ
ル変換部2は、電力系統から取り込む各種電気量の入力
数に応じて複数設けられる。12はこのディジタルリレ
ー1が保護する電力系統であり、また13は電力系統1
2から電流信号を取り込む変流器CTである。なお、こ
のディジタルリレーを、電力系統を保護するリレーとし
て使用する場合には、電流信号の他に電圧信号が必要と
なるため、変流器13とは別に変成器が設けられて、電
圧信号が取り込まれる。またディジタルリレーをこのよ
うに系統保護リレーとして使用する場合と変圧器保護リ
レーとして使用する場合とでは、後述の演算処理部8で
の演算および判定方法に違いがあるものの、各部の基本
的構成は同様である。
DESCRIPTION OF THE PREFERRED EMBODIMENTS Next, the present invention will be described with reference to the drawings. FIG. 1 is a block diagram of a digital relay which is an embodiment of the present invention. In FIG. 1, reference numeral 1 is a digital relay, which is composed of an analog / digital conversion unit 2 and a digital operation control unit 7. A plurality of analog / digital conversion units 2 are provided according to the number of inputs of various electric quantities taken in from the power system. 12 is a power system protected by this digital relay 1, and 13 is a power system 1
2 is a current transformer CT that takes in a current signal from 2. In addition, when this digital relay is used as a relay for protecting the power system, a voltage signal is required in addition to the current signal. Therefore, a transformer is provided separately from the current transformer 13, and the voltage signal is transmitted. It is captured. Although the digital relay is used as a system protection relay and the transformer protection relay is used as a transformer protection relay, there is a difference in the calculation and determination method in the calculation processing unit 8 described later, but the basic configuration of each unit is It is the same.

【0015】変流器13からの電流信号は入力変換器3
で後段の回路に適切な振幅の電圧信号に変換され、続く
アナログフィルタ4でこの信号に含まれる高次高調波分
(高調波ノイズ)および直流分が除去される。さらに、
サンプルホールド回路5で所定の時間間隔、例えば電力
基本周期(1/50Hz)の1/12の時間間隔(電気
角30゜)でサンプリングされ、次のサンプリングまで
出力保持される。続いてこのサンプリング信号はA/D
変換器6でアナログ信号からディジタル信号に変換さ
れ、ディジタル演算制御部7の演算処理部8に入力され
る。さらに、マイクロコンピュータから構成されている
演算処理部8では、ROM9に格納されたプログラムに
従って、RAM10に格納したデータに基づき後述する
演算が行なわれ、その結果に応じて対応する系統を遮断
すべきか否かが判断される。遮断すべきであると判断さ
れた場合には、11の外部インタフェースからこの事故
系統を遮断するための信号が対応する遮断器へ出力され
る。
The current signal from the current transformer 13 is supplied to the input converter 3
Is converted into a voltage signal having an appropriate amplitude by the circuit in the subsequent stage, and the high-order harmonic component (harmonic noise) and DC component contained in this signal are removed by the subsequent analog filter 4. further,
The sample and hold circuit 5 performs sampling at a predetermined time interval, for example, a time interval (electrical angle 30 °) of 1/12 of a power basic cycle (1/50 Hz), and the output is held until the next sampling. Then this sampling signal is A / D
The converter 6 converts the analog signal into a digital signal and inputs it to the arithmetic processing unit 8 of the digital arithmetic control unit 7. Further, in the arithmetic processing unit 8 composed of a microcomputer, the later-described arithmetic operation is performed based on the data stored in the RAM 10 according to the program stored in the ROM 9, and whether or not the corresponding system should be shut off according to the result. Is determined. When it is determined that the circuit should be cut off, a signal for cutting off the fault system is output from the external interface 11 to the corresponding breaker.

【0016】次に、本発明のディジタルリレーが図6の
ように電力系統保護リレーとして使用される場合、図1
の演算処理部8において、系統のインピーダンスに含ま
れる基本波分に対する高調波分の割合が算出され、その
結果によりリレーの動作判定がなされる。前述のよう
に、インダクタンス分Lあるいはコンデンサ分Cが増大
した系統においては、インピーダンスを算出するために
電力系統から取り込まれた電圧と電流には周波数のズレ
が生じるため、このインピーダンスベクトルが移動する
現象が発生する。この移動過程を複素数を用いて計算す
ると式(1)のようになる。式(1)において、Zはイ
ンピーダンスベクトル、V1 は事故時の電圧入力信号に
含まれる電源周波数の基本波分、同じくV2 は事故時の
電圧入力信号に含まれる高調波分をそれぞれ示してお
り、さらにIは事故時の電流入力を示している。また、
ωは電源の基本波角速度、tは時間、さらにnは高調波
の次数をそれぞれ示している。
Next, when the digital relay of the present invention is used as a power system protection relay as shown in FIG.
In the arithmetic processing unit 8, the ratio of the harmonic component to the fundamental component contained in the system impedance is calculated, and the relay operation is determined based on the result. As described above, in the system in which the inductance component L or the capacitor component C is increased, the voltage and the current taken in from the power system to calculate the impedance have a frequency shift, and the impedance vector moves. Occurs. When this moving process is calculated using a complex number, it becomes as shown in Expression (1). In equation (1), Z is the impedance vector, V 1 is the fundamental wave component of the power supply frequency included in the voltage input signal at the time of the accident, and V 2 is the harmonic component included in the voltage input signal at the time of the accident. Further, I indicates the current input at the time of the accident. Also,
ω is the angular velocity of the fundamental wave of the power source, t is time, and n is the order of harmonics.

【0017】[0017]

【数1】 [Equation 1]

【0018】ここで、この式(1)の形を見ると、イン
ピーダンスベクトルZは、固定ベクトルと回転ベクトル
とによって合成されたものであることがわかる。すなわ
ち、v2 /iは固定ベクトルを、またv1 /iは時計方
向に(n−1)時の高調波次数で回転するベクトルをそ
れぞれ意味しており、これを横軸を実軸,縦軸を虚軸と
したR−X座標に表現すると図3のようになる。同図に
おいて、インピーダンスベクトルZは、3a点で最大と
なり、3b点で最小となる軌跡を描くものであることが
わかる。また、これらの各点でのインピーダンスをZa
およびZb とすると、 最大値:Za =v2 /i+v1 /i 最小値:Zb =v2 /i−v1 /i となる。
Looking at the form of this equation (1), it can be seen that the impedance vector Z is a composite of a fixed vector and a rotation vector. That is, v 2 / i means a fixed vector, and v 1 / i means a vector that rotates in the harmonic order at (n-1) time in the clockwise direction. The horizontal axis is the real axis and the vertical axis is the vertical axis. When expressed in RX coordinates with the axis as an imaginary axis, it becomes as shown in FIG. In the figure, it can be seen that the impedance vector Z draws a locus having a maximum at the point 3a and a minimum at the point 3b. In addition, the impedance at each of these points is Z a
And Z b , the maximum value: Z a = v 2 / i + v 1 / i and the minimum value: Z b = v 2 / i−v 1 / i.

【0019】ここで、これらインピーダンスベクトルZ
a ,Zb 、およびそれぞれのベクトルが発生した時間T
za,Tzbは、電力系統から取り込んだ電圧信号および電
流信号のサンプリングデータから求めることが可能であ
ることから、v2 /i,v1/iがそれぞれ逆算可能と
なる。すなわち、v2 /iの大きさをZaba 、v1 /i
の大きさをZabs とすると、 Zaba =(Za +Zb )/2 Zabs =(Za −Zb )/2 と表すことができる。
Here, these impedance vectors Z
a , Z b , and the time T at which each vector occurred
Since za and T zb can be obtained from the sampling data of the voltage signal and the current signal fetched from the power system, v 2 / i and v 1 / i can be respectively back-calculated. That is, the magnitude of v 2 / i is Z aba , v 1 / i
Of the magnitude When Z abs, can be expressed as Z aba = (Z a + Z b) / 2 Z abs = (Z a -Z b) / 2.

【0020】さらに、インピーダンスベクトルZに含ま
れる高調波分の基本波分に対する割合Zf とその高調波
次数nは、f=50Hz基本波ベースでのZa ,Zb
の移動時間をTZab とすると、 Zf =Zaba /Zabszab =Tza−Tzb=1/{2・f(n−1)}(s) =1000/{2・50(n−1)}(ms) となり、これより高調波分の割合Zf とその次数nが算
出可能となる。
Further, the ratio Z f of the harmonic component contained in the impedance vector Z to the fundamental component and its harmonic order n are the moving time between Z a and Z b on the basis of f = 50 Hz fundamental component, T Zab Then, Z f = Z aba / Z abs T zab = T za −T zb = 1 / {2 · f (n−1)} (s) = 1000 / {2 · 50 (n−1)} (ms ), The ratio of harmonic components Z f and its order n can be calculated.

【0021】次に、本発明における演算処理の手順につ
いて、図2を参考に説明する。ステップ21において、
図1のアナログ/ディジタル変換部2でサンプリングさ
れた電圧および電流データをRAM10に格納する。次
にステップ22でこのデータおよびこれより前に取り込
まれ格納してあったデータに基づき、インピーダンスベ
クトルが算出され、RAM10に格納される。なお、サ
ンプリングデータからベクトルを算出する方法は、従来
からのベクトル計算法に基づくものである。
Next, the procedure of the arithmetic processing in the present invention will be described with reference to FIG. In step 21,
The voltage and current data sampled by the analog / digital converter 2 in FIG. 1 is stored in the RAM 10. Next, in step 22, an impedance vector is calculated based on this data and the data previously fetched and stored, and stored in the RAM 10. The method of calculating the vector from the sampling data is based on the conventional vector calculation method.

【0022】すなわち、ある電圧サンプル値v(m) より
n回前のサンプル値をv(m-n) とすると、サンプリン
グ間隔が電気角30゜の場合、 V2 =v(m)2+v(m-3)22 =i(m)2+i(m-3)2 VIcos θ=v(m)・i(m) +v(m-3)・i(m-3) VIsin θ=v(m-3)・i(m) −v(m)・i(m-3) となり、インピーダンスZ、抵抗R、リアクタンスX
は、 Z2 =V2 /I22 =Zcos θ=(VIcos θ)/I22 =Zsin θ=(VIsin θ)/I により算出される。
That is, assuming that a sample value n times before a certain voltage sample value v (m) is v (mn), V 2 = v (m) 2 + v (m- 3) 2 I 2 = i (m) 2 + i (m-3) 2 VIcos θ = v (m) ・ i (m) + v (m-3) ・ i (m-3) VIsin θ = v (m- 3) ・ i (m) -v (m) ・ i (m-3), impedance Z, resistance R, reactance X
Is calculated by: Z 2 = V 2 / I 2 R 2 = Zcos θ = (VIcos θ) / I 2 X 2 = Zsin θ = (VIsin θ) / I 2 .

【0023】続くステップ23,24で、これら算出し
た値のうち過去一定時間内に、すなわち事故発生からイ
ンピーダンスベクトルの軌跡がリレー動作領域を通過す
るのに充分な時間内における最大値および最小値が抽出
され、その値と発生した時刻とがステップ25,26で
それぞれ記憶される。そしてステップ27でこれらの値
から前述の演算方法に基づき、高調波分の割合とその次
数が算出される。以上の演算処理により、電力系統から
の各種電気量から、インピーダンスに含まれる高調波分
の割合とその次数を算出することが可能となる。従っ
て、これに基づき特定された事故の発生地点が保護区間
内か否かが判断され、電力系統保護の要否が判定され
る。
In the following steps 23 and 24, the maximum value and the minimum value of these calculated values within the past fixed time, that is, within the time sufficient for the locus of the impedance vector to pass through the relay operation area from the occurrence of the accident are determined. The extracted value and the time of occurrence are stored in steps 25 and 26, respectively. Then, in step 27, the ratio of the harmonic component and its order are calculated from these values based on the above-described calculation method. With the above calculation processing, it is possible to calculate the proportion of harmonics contained in impedance and its order from various amounts of electricity from the power system. Therefore, it is determined whether or not the accident occurrence point identified based on this is within the protection section, and whether or not the power system protection is necessary is determined.

【0024】また、本発明のディジタルリレーが、図7
のような変圧器保護リレーとして使用される場合、図1
の演算処理部8において、差動電流中の基本波分に対す
る高調波分の割合が算出され、その結果に応じてリレー
の動作判定がなされる。今、電源側からI1f,線路側
からInfの高調波分がリレーに取り込まれたとすると、
差動電流Id は、Id =I1f+Infとなる。このId
所定の電源基本波ベクトル(大きさ1,反時計方向に5
0回転/秒)で除算した演算ベクトルIddを複素数を用
いて表現すると式(2)のようになる。式(2)におい
て、ωは電源の基本波角速度、tは時間、さらにnは高
調波の次数をそれぞれ示している。
The digital relay of the present invention is shown in FIG.
When used as a transformer protection relay such as
In the arithmetic processing unit 8, the ratio of the harmonic component to the fundamental component in the differential current is calculated, and the operation of the relay is determined according to the result. Now, assuming that the harmonic components of I 1f from the power source side and In f from the line side are taken in by the relay,
The differential current I d is I d = I 1f + I nf . This I d is set to a predetermined power source fundamental wave vector (magnitude 1, 5 in the counterclockwise direction).
When the operation vector I dd divided by 0 revolutions / second) is expressed using a complex number, it becomes as shown in Expression (2). In Expression (2), ω indicates the fundamental wave angular velocity of the power source, t indicates time, and n indicates the order of harmonics.

【0025】[0025]

【数2】 [Equation 2]

【0026】ここで、この式(2)の形を見ると、演算
ベクトルIddは、固定ベクトルと回転ベクトルとによっ
て合成されたものであることがわかる。すなわち、I1f
は固定ベクトルを、またInfは反時計方向に回転するベ
クトルをそれぞれ意味しており、これをR−X座標に表
現すると図5のようになる。図5において、演算ベクト
ルIddは、5a点で最大となり、5b点で最小となる軌
跡を描くものであることがわかる。また、これらの各点
での演算ベクトルをIddmax およびIddmin とすると、 最大値:Iddmax =I1f+Inf 最小値:Iddmin =I1f−Inf となる。
Looking at the form of this equation (2), it can be seen that the operation vector I dd is a combination of a fixed vector and a rotation vector. That is, I 1f
Represents a fixed vector, and Inf represents a vector that rotates counterclockwise, which is expressed in RX coordinates as shown in FIG. In FIG. 5, it can be seen that the operation vector I dd draws a locus having a maximum at the point 5a and a minimum at the point 5b. Further, when the operation vectors at these points are I ddmax and I ddmin , the maximum value is I ddmax = I 1f + I nf The minimum value is I ddmin = I 1f −I nf .

【0027】ここで、このIddmax ,Iddmin 、および
それぞれのベクトルが発生した時間Tia ,Tib は、電
力系統から取り込んだ電流信号のサンプリングデータか
ら求めることが可能であることから、I1f,Infがそれ
ぞれ逆算可能となる。すなわち、I1fの大きさをI
dda 、Infの大きさをIdds とすると、 Idda =(Iddmax +Iddmin )/2 Idds =(Iddmax −Iddmin )/2 と表すことができる。
Since I ddmax and I ddmin and the times T ia and T ib at which the respective vectors occur can be obtained from the sampling data of the current signal taken from the power system, I 1f , I nf can be respectively back-calculated. That is, the size of I 1f is I
If the magnitudes of dda and I nf are I dds , then I dda = (I ddmax + I ddmin ) / 2 I dds = (I ddmax −I ddmin ) / 2 can be expressed.

【0028】さらに、演算ベクトルIddに含まれる高調
波分の基本波分に対する割合If とその次数nは、f=
50Hz基本波ベースでの5a,5b間の移動時間をT
iabとすると、 If =Idds /Iddaiab =Tia−Tib=1/{2・f(n−1)}(s) =1000/{2・50(n−1)}(ms) となり、これより高調波分の割合If とその次数nが算
出可能となる。
Further, the ratio If of the harmonics contained in the operation vector Idd to the fundamental wave and its order n are f =
The moving time between 5a and 5b based on the 50 Hz fundamental wave is T
If iab , then I f = I dds / I dda T iab = T ia −T ib = 1 / {2 · f (n−1)} (s) = 1000 / {2 · 50 (n−1)} ( ms), from which the ratio If of harmonics and its order n can be calculated.

【0029】次に、本発明における演算処理の手順につ
いて、図4を参考に説明する。ステップ41において、
図1のアナログ/ディジタル変換部2でサンプリングさ
れた電流データをRAM10に格納する。次にステップ
42でこのデータおよびこれより前に取り込まれ格納し
てあったデータに基づき、演算ベクトルが前述のベクト
ル計算法により算出され、RAM10に格納される。続
くステップ43,44で、これら算出した値のうち過去
一定時間内に、すなわち事故発生から演算ベクトルの軌
跡がリレー動作領域を通過するのに充分な時間内におけ
る最大値および最小値が抽出され、その値と発生した時
刻とがステップ45,46でそれぞれ記憶される。そし
てステップ47でこれらの値から前述の演算方法に基づ
き、高調波分の割合とその次数が算出される。
Next, the procedure of the arithmetic processing in the present invention will be described with reference to FIG. In step 41,
The current data sampled by the analog / digital converter 2 in FIG. 1 is stored in the RAM 10. Next, in step 42, an operation vector is calculated by the above-described vector calculation method based on this data and the data previously fetched and stored, and stored in the RAM 10. In the following steps 43 and 44, the maximum value and the minimum value are extracted from the calculated values within the past fixed time, that is, within the time sufficient for the trajectory of the operation vector to pass through the relay operation area from the accident occurrence, The value and the time of occurrence are stored in steps 45 and 46, respectively. Then, in step 47, the ratio of harmonics and the order thereof are calculated from these values based on the above-described calculation method.

【0030】以上の演算処理により、電力系統からの各
種電気量から、差動電流に含まれる高調波分の割合とそ
の次数を算出することが可能となる。従って、これに基
づき、この差動電流中の高調波分が、変圧器特有の単な
る電圧印可時の鉄心飽和による固有のものか、前述のよ
うなインダクタンス分Lあるいはコンデンサ分Cが増大
した系統における事故時のものかを判断することによ
り、変圧器事故発生の判定、すなわち電力系統保護の要
否が判定可能となる。
By the above arithmetic processing, it becomes possible to calculate the ratio and the order of the harmonic components contained in the differential current from various electric quantities from the electric power system. Therefore, based on this, whether the harmonic component in the differential current is peculiar to the transformer due to the saturation of the iron core when voltage is applied or in the system in which the inductance component L or the capacitor component C is increased as described above. By determining whether or not it is at the time of an accident, it is possible to determine whether a transformer accident has occurred, that is, whether or not power system protection is necessary.

【0031】[0031]

【発明の効果】以上説明したように、本発明は、電力系
統からサンプリングされた電圧および電流データに基づ
き、電圧ベクトルを電流ベクトルで除算したインピーダ
ンスベクトルを算出するベクトル演算手段を設けて、こ
のインピーダンスベクトルの最大値および最小値ならび
にその時刻を抽出し、このベクトルに含まれる高調波分
の割合とその次数を算出するようにしたので、インピー
ダンスベクトルの算出精度が改善される。従って、線路
側系統のインダクタンス分やコンデンサ分が増大してい
る系統に使用されている系統保護継電器において、事故
発生地点を特定する精度が向上するとともに、誤動作防
止のための強制抑止手段を必要とすることなく、系統保
護継電器の正確な動作判定が可能となるという顕著な効
果を奏するものである。
As described above, according to the present invention, the vector calculating means for calculating the impedance vector by dividing the voltage vector by the current vector based on the voltage and current data sampled from the power system is provided, and the impedance is calculated. Since the maximum value and the minimum value of the vector and the time thereof are extracted and the ratio of the harmonic components contained in this vector and its order are calculated, the calculation accuracy of the impedance vector is improved. Therefore, in the system protection relay used in the system where the inductance and the capacitor of the line side system are increasing, the accuracy of identifying the accident occurrence point is improved and the forced deterrent means for preventing malfunction is required. Without doing so, it is possible to accurately determine the operation of the system protection relay, which is a remarkable effect.

【0032】また、電力系統からサンプリングされた電
流データに基づき、電流ベクトルを電力系統の基本波基
準ベクトルで除算した演算ベクトルを算出するベクトル
演算手段を設けて、この演算ベクトルの最大値および最
小値ならびにその時刻を抽出し、このベクトルに含まれ
る高調波分の割合とその次数を算出するようにしたの
で、高調波分の判別精度が改善される。従って、変圧器
保護継電器において、変圧器特有の励磁突流現象に起因
する高調波分と、線路側系統のインダクタンス分やコン
デンサ分が増大している系統における事故時に発生する
高周波分とを精度よく判別することが可能となり、誤動
作防止のための強制抑止手段を必要とすることなく、変
圧器保護継電器の正確な動作判定が可能となるという顕
著な効果を奏するものである。
Further, based on the current data sampled from the electric power system, vector operation means for calculating an operation vector by dividing the current vector by the fundamental wave reference vector of the electric power system is provided, and the maximum value and the minimum value of this operation vector are provided. Also, since the time is extracted and the ratio of the harmonic components contained in this vector and its order are calculated, the accuracy of discrimination of the harmonic components is improved. Therefore, in the transformer protection relay, it is possible to accurately distinguish between the harmonic component caused by the exciting surge current phenomenon peculiar to the transformer and the high frequency component generated at the time of an accident in the system where the line side system inductance and capacitor are increasing. Therefore, there is a remarkable effect that the operation of the transformer protection relay can be accurately determined without the need of the force suppression means for preventing malfunction.

【図面の簡単な説明】[Brief description of drawings]

【図1】本発明の一実施例のブロック図である。FIG. 1 is a block diagram of an embodiment of the present invention.

【図2】本発明が系統保護リレーとして使用される場合
の、図1の演算処理部における演算処理手順を示す図で
ある。
FIG. 2 is a diagram showing an arithmetic processing procedure in the arithmetic processing unit of FIG. 1 when the present invention is used as a system protection relay.

【図3】図2の手順により算出されるインピーダンスベ
クトルの軌跡を示す図である。
FIG. 3 is a diagram showing a locus of an impedance vector calculated by the procedure of FIG.

【図4】本発明が変圧器保護リレーとして使用される場
合の、図1の演算処理部における演算処理手順を示す図
である。
FIG. 4 is a diagram showing an arithmetic processing procedure in the arithmetic processing unit of FIG. 1 when the present invention is used as a transformer protection relay.

【図5】図4の手順により算出される演算ベクトルの軌
跡を示す図である。
5 is a diagram showing a locus of a calculation vector calculated by the procedure of FIG.

【図6】系統保護リレーの接続状態を示す図である。FIG. 6 is a diagram showing a connection state of a system protection relay.

【図7】変圧器保護リレーの接続状態を示す図である。FIG. 7 is a diagram showing a connection state of a transformer protection relay.

【図8】従来のディジタルリレーにおける低次高調波分
を抽出するディジタルフィルタの利得/周波数特性を示
す図である。
FIG. 8 is a diagram showing gain / frequency characteristics of a digital filter for extracting low-order harmonic components in a conventional digital relay.

【図9】従来の距離リレーにおける動作特性とインピー
ダンス軌跡を示す図である。
FIG. 9 is a diagram showing operating characteristics and impedance loci in a conventional distance relay.

【符号の説明】[Explanation of symbols]

1 ディジタルリレー 2 アナログ/ディジタル変換部 7 ディジタル演算制御部 8 演算処理部 9 ROM 10 RAM 11 外部インターフェース 12 電力系統 13 変流器 1 Digital Relay 2 Analog / Digital Converter 7 Digital Operation Control 8 Operation Processing 9 ROM 10 RAM 11 External Interface 12 Power System 13 Current Transformer

Claims (2)

【特許請求の範囲】[Claims] 【請求項1】 電力系統からサンプリングした各種電気
量に対して演算を行うとともに、その演算結果に応じて
動作の要否を判定する保護継電器において、サンプリン
グされた電圧および電流データに基づき、電圧ベクトル
を電流ベクトルで除算したインピーダンスベクトルを算
出するベクトル演算手段と、前記インピーダンスベクト
ルの最大値および最小値ならびにその時刻をそれぞれ抽
出する抽出手段と、この抽出結果に基づき前記インピー
ダンスベクトルに含まれる高調波分の割合とその次数を
算出する高調波分算出手段とを備えることを特徴とする
保護継電器。
1. A protection relay that performs calculation on various amounts of electricity sampled from a power system and determines whether or not operation is required based on the calculation result, based on sampled voltage and current data, and voltage vector. Is calculated by dividing the current vector by an impedance vector, a vector arithmetic means for extracting the maximum value and the minimum value of the impedance vector and its time, and a harmonic component contained in the impedance vector based on the extraction result. And a harmonic component calculating means for calculating the order thereof and a protection relay.
【請求項2】 電力系統からサンプリングした各種電気
量に対して演算を行うとともに、その演算結果に応じて
動作の要否を判定する保護継電器において、サンプリン
グされた電流データに基づき、電流ベクトルを電力系統
の基本波基準ベクトルで除算した演算ベクトルを算出す
るベクトル演算手段と、前記演算ベクトルの最大値およ
び最小値ならびにその時刻をそれぞれ抽出する抽出手段
と、この抽出結果に基づき前記演算ベクトルに含まれる
高調波分の割合とその次数を算出する高調波分算出手段
とを備えることを特徴とする保護継電器。
2. A protection relay that performs an operation on various electric quantities sampled from a power system and determines whether or not an operation is necessary according to the operation result, and outputs a current vector based on the sampled current data. Vector operation means for calculating the operation vector divided by the fundamental wave reference vector of the system, extraction means for extracting the maximum value and the minimum value of the operation vector and the time thereof, respectively, and included in the operation vector based on the extraction result A protective relay, comprising: a harmonic component ratio and a harmonic component calculating means for calculating the order thereof.
JP4263310A 1992-10-01 1992-10-01 Protective relay Expired - Fee Related JP3003421B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP4263310A JP3003421B2 (en) 1992-10-01 1992-10-01 Protective relay

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP4263310A JP3003421B2 (en) 1992-10-01 1992-10-01 Protective relay

Publications (2)

Publication Number Publication Date
JPH06121455A true JPH06121455A (en) 1994-04-28
JP3003421B2 JP3003421B2 (en) 2000-01-31

Family

ID=17387707

Family Applications (1)

Application Number Title Priority Date Filing Date
JP4263310A Expired - Fee Related JP3003421B2 (en) 1992-10-01 1992-10-01 Protective relay

Country Status (1)

Country Link
JP (1) JP3003421B2 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008521374A (en) * 2004-11-18 2008-06-19 ドン・エナジー・セールス・アンド・ディストリビューション・アクティーゼルスカブ Compensation for a simple fiber optic Faraday effect sensor
US8447370B2 (en) 2007-01-04 2013-05-21 Bose Corporation Microphone techniques
CN105226607A (en) * 2015-10-22 2016-01-06 许继电气股份有限公司 A kind of extra-high voltage direct-current layer-specific access system 100Hz harmonic protection method

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS61161922A (en) * 1984-12-29 1986-07-22 株式会社東芝 Protective relay

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS61161922A (en) * 1984-12-29 1986-07-22 株式会社東芝 Protective relay

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008521374A (en) * 2004-11-18 2008-06-19 ドン・エナジー・セールス・アンド・ディストリビューション・アクティーゼルスカブ Compensation for a simple fiber optic Faraday effect sensor
JP4869245B2 (en) * 2004-11-18 2012-02-08 パワーセンス・アクティーゼルスカブ Compensation for a simple fiber optic Faraday effect sensor
US8447370B2 (en) 2007-01-04 2013-05-21 Bose Corporation Microphone techniques
CN105226607A (en) * 2015-10-22 2016-01-06 许继电气股份有限公司 A kind of extra-high voltage direct-current layer-specific access system 100Hz harmonic protection method
CN105226607B (en) * 2015-10-22 2017-09-12 许继电气股份有限公司 A kind of extra-high voltage direct-current layer-specific access system 100Hz harmonic protection methods

Also Published As

Publication number Publication date
JP3003421B2 (en) 2000-01-31

Similar Documents

Publication Publication Date Title
US11885843B2 (en) Method, system and apparatus for fault detection
CN107069734A (en) A kind of active filter and method with lightning protection isolation with arc-detection function
CN107632234A (en) A kind of deformation of transformer winding appraisal procedure based on recorder data
JP3808624B2 (en) System protection relay device
JP2004080839A (en) Ground direction relay and ground direction relay device
US4450497A (en) Ultra-high-speed relay
JP2004048855A (en) Distance relaying apparatus
JP3003421B2 (en) Protective relay
Mahmoud Detection and assessment scheme of voltage and current unbalance for three phase synchronous generators using dual numerical techniques based on correlation and deviation percentage concepts
EP1311866B1 (en) Current compensation method and device for power system protection
US11355919B2 (en) Current differential relay apparatus
EP0952654B1 (en) Bus-bar protective relay apparatus
Kennedy et al. Variable digital filter response time in a digital distance relay
JP3792893B2 (en) Busbar protection relay device
JP3160612B2 (en) Power system insulation deterioration detection method and apparatus
EP1078431B1 (en) Mimic high pass filter in a protective relay
JP4561366B2 (en) Digital ground fault current protection relay
US6154687A (en) Modified cosine filters
JP3132011B2 (en) Harmonic filter equipment protection relay
JP2863952B2 (en) Ground fault fault location method and device, ground fault distance relay
Moravej et al. Distance protection closed-loop testing using RTDS
JP2588414B2 (en) Power system accident judgment method
JP3558309B2 (en) Ratio differential relay
Samonto et al. Fuzzy Based Algorithm for Stage-I of Cascaded Intelligent Relaying
EP1295374B1 (en) Method and device for power system protection

Legal Events

Date Code Title Description
LAPS Cancellation because of no payment of annual fees