JP3003421B2 - Protective relay - Google Patents

Protective relay

Info

Publication number
JP3003421B2
JP3003421B2 JP4263310A JP26331092A JP3003421B2 JP 3003421 B2 JP3003421 B2 JP 3003421B2 JP 4263310 A JP4263310 A JP 4263310A JP 26331092 A JP26331092 A JP 26331092A JP 3003421 B2 JP3003421 B2 JP 3003421B2
Authority
JP
Japan
Prior art keywords
vector
relay
current
calculating
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.)
Expired - Fee Related
Application number
JP4263310A
Other languages
Japanese (ja)
Other versions
JPH06121455A (en
Inventor
敏信 海老坂
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)

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 for detecting an accident in an electric power system and outputting a removal command for a corresponding accident section, and in particular, digitizes and calculates various electric quantity data sampled from the electric power system. The present invention relates to a digital operation type protection relay that determines an operation 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 an electric 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 separate the part from the grid, power protection relays are arranged at various places according to the configuration of the grid. Various types of protection relays exist according to the operation principle. One of them is a digital operation type protection relay (hereinafter referred to as a digital relay). This involves sampling various system electrical quantities as digital data, performing arithmetic processing according to a predetermined arithmetic procedure, and determining the operation of the relay based on the results. The technical principle is as follows: (Ohmsha / H3.8.25 first edition).

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

【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 supply.
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). 65 is a reactance L in the line side system, and 66 is a capacitor C similarly.
Is shown. According to the above configuration, a ratio between 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. Thereby, 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次の高
調波分を示している。
[0005] In such a general power system, in recent years, the voltage class has increased and the cable system has tended to increase, and as a result, the reactance L and the capacitor C of the line side system have increased, resulting in an accident. There is an increasing number of systems that generate low-order harmonics at a frequency of about 1.5 times the fundamental frequency f of the power supply at the time of generation. In FIG. 6, this relay 6
An accident occurs at a point Fa behind the line-side system protected by the line 4, and a current Inf of nf from the line side to the point Fa and a power supply 61 due to the reactance 65 and the capacitor 66 of the line-side system. It is assumed that a current I 1f is generated from to the point Fa. The term “nf” indicates an 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,
The voltage which is affected by the power supply 61 and includes more 1 f than the 1.5 f is introduced to the relay 64 via the transformer 62. In response to this, in the relay 64, the impedance is calculated as described above, and the operation is determined. Normally, the direction of the current Inf flows in the opposite direction to that at the time of occurrence of an accident in the line-side system to be protected, that is, flows toward the power supply 61, so that the calculated impedance becomes a value other than the predetermined value, It is determined that the accident location is outside the protection zone, and relay 64 does not operate.

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

【0008】また、ディジタルリレーはこのような系統
保護リレーとしての配置される他に、図7に示すように
変圧器を保護する場合にも用いられる。同図に用いられ
るリレーは、保護する変圧器の両端の電流の差が一定値
以上になったとき動作するように構成されており、差動
リレーとも呼ばれている。図7において、71は電源を
示しており、72および73はリレー74に変成器77
の両端の電流を導入する変流器である。リレー74はデ
ィジタル演算形の変圧器保護リレー(TP)であり、変
圧器77は系統に直列に配置され、このリレー74によ
って保護される。75は線路側系統におけるリアクタン
ス分L、76は同じくコンデンサ分Cを示す。以上の構
成により、変流器72および73からの電流に応じて、
変圧器77への流入電流の和、すなわち差動電流が算出
され、この算出結果に基づく動作判定により、変圧器7
7が保護されている。
[0008] In addition to being arranged as such a system protection relay, a digital relay is also used for protecting a transformer as shown in FIG. The relay used in the figure is configured to operate when the difference between the currents at both ends of the transformer to be protected exceeds a certain value, and is also called a differential relay. In FIG. 7, reference numeral 71 denotes a power supply, and reference numerals 72 and 73 denote a transformer 77 and a relay 74.
Is a current transformer that introduces a current at both ends of the current transformer. The relay 74 is a digital operation type transformer protection relay (TP). The transformer 77 is arranged in series in a system, and is protected by the relay 74. Reference numeral 75 denotes a reactance L in the line side system, and reference numeral 76 denotes a capacitor C in the same manner. With the above configuration, according to the currents from current transformers 72 and 73,
The sum of the inflow currents into the transformer 77, that is, the differential current is calculated, and the operation is determined 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, the transformer protection relay 74
An accident occurs at the Fb point in or near the transformer 77, that is, the transformer 77, and the current I nf for nf from the line-side system toward the Fb point due to the reactance 75 and the capacitor 76 of the line-side system. , And power supply 71
Suppose that a current I 1f is generated from 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 The current is introduced to the relay 74 via the current transformer 72. In response to this, in the relay 74, the differential current is calculated as described above, and the operation is determined. Normally, when the current Inf increases, 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 a conventional digital relay, the influence of an exciting sudden flow phenomenon peculiar to such a transformer is eliminated. For this purpose, a digital filter having a filter characteristic as shown in FIG.
Means are provided for extracting a considerable amount of low-order harmonics and forcibly inhibiting the operation of the relay 74 when the value of the ratio I 2f / I 1f with respect to the separately extracted 1f is equal to or more than a certain value. Was.

【0011】[0011]

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

【0012】[0012]

【課題を解決するための手段】このような目的を達成す
るために、本発明による保護継電器は、電力系統からサ
ンプリングされた電圧および電流データに基づき、電圧
ベクトルを電流ベクトルで除算したインピーダンスベク
トルを算出するベクトル演算手段と、インピーダンスベ
クトルの最大値および最小値ならびにこれら発生時刻
それぞれ抽出する抽出手段と、この抽出手段により抽出
されたインピーダンスベクトルの最大値および最小値か
らインピーダンスベクトルに含まれる高調波分の割合を
算出するとともに、インピーダンスベクトルの最大値お
よび最小値の発生時刻から高調波分の次数を算出する
調波分算出手段とを備えるものである。また、電力系統
からサンプリングされた電流データに基づき、電流ベク
トルを電力系統の基本波基準ベクトルで除算した演算ベ
クトルを算出するベクトル演算手段と、演算ベクトルの
最大値および最小値ならびにこれら発生時刻をそれぞれ
抽出する抽出手段と、この抽出手段により抽出された演
算ベクトルの最大値および最小値から演算ベクトルに含
まれる高調波分の割合を算出するとともに、演算ベクト
ルの最大値および最小値の発生時刻から高調波分の次数
を算出する高調波分算出手段とを備えるものである。
In order to achieve the above object, a protection 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. a vector calculating means for calculating for, extracting means for extracting respective maximum and minimum values and their time of occurrence of the impedance vector extracted by the extraction means
Maximum and minimum values of the calculated impedance vector
The ratio of the harmonic components contained in the impedance vector
Calculation and the maximum value of the impedance vector
And a harmonic component calculating means for calculating the order of harmonic components from the occurrence time of the minimum value . Further, based on the current data sampled from the power system, a vector operation means for calculating an operation vector obtained 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 the time at which these are generated , respectively. Extraction means for extracting, and the performances extracted by the extraction means.
From the maximum and minimum values of the
Calculate the percentage of harmonics
The order of the harmonic component from the time when the maximum and minimum values of the
And a harmonic component calculating means for calculating

【0013】[0013]

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

【0014】[0014]

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

【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
Then, the signal is converted into a voltage signal having an amplitude suitable for the circuit at the subsequent stage, and the high-order harmonic component (harmonic noise) and the DC component contained in the signal are removed by the analog filter 4. further,
The sample and hold circuit 5 performs sampling at a predetermined time interval, for example, at a time interval of 1/12 (electrical angle 30 °) of the basic power cycle (1/50 Hz), and the output is held until the next sampling. Subsequently, this sampling signal is A / D
The signal is converted from an analog signal to a digital signal by the converter 6 and is input to the arithmetic processing unit 8 of the digital arithmetic control unit 7. Further, in an arithmetic processing unit 8 composed of a microcomputer, an operation described later is performed based on the data stored in the RAM 10 in accordance with a program stored in the ROM 9 and, depending on the result, whether or not the corresponding system should be shut off. Is determined. If it is determined that the fault should be interrupted, a signal for shutting down the accident system is output from the external interface 11 to the corresponding circuit 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 included in the impedance of the system is calculated, and the operation of the relay is determined based on the calculation result. As described above, in a system in which the inductance L or the capacitor C is increased, the voltage and the current taken from the power system for calculating the impedance have a frequency difference, so that the impedance vector shifts. Occurs. When this movement process is calculated using a complex number, it becomes as shown in equation (1). In the equation (1), Z is an impedance vector, V 1 is a fundamental component of a power supply frequency included in the voltage input signal at the time of an accident, and V 2 is a harmonic component included in the voltage input signal at the time of an accident. And I indicates the current input at the time of the accident. Also,
ω indicates the fundamental wave angular velocity of the power supply, t indicates time, and n indicates 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 となる。
Here, looking at the form of the equation (1), it can be seen that the impedance vector Z is a combination of the fixed vector and the rotation vector. That is, v 2 / i means a fixed vector, and v 1 / i means a vector rotating clockwise at the harmonic order at (n−1). FIG. 3 shows the R-X coordinates with the axis being an imaginary axis. In the figure, it can be seen that the impedance vector Z draws a locus which becomes maximum at the point 3a and becomes minimum at the point 3b. The impedance at each of these points is represented by Z a
And Z b , the maximum value: Z a = v 2 / i + v 1 / i 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 taken in from the power system, v 2 / i and v 1 / i can be back calculated. That is, the magnitude of v 2 / i is Z aba , v 1 / i
Letting the magnitude of Z abs be 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が算
出可能となる。
Furthermore, the ratio Z f and its harmonic order n of the fundamental wave component of the harmonic component contained in the vector impedance Z is, f = 50 Hz fundamental basis of Z a, the travel time between Z b 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 ), From which the ratio of harmonic components Zf 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 of FIG. Next, at step 22, an impedance vector is calculated based on this data and the data previously captured and stored, and stored in the RAM 10. The method of calculating a vector from sampling data is based on a 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), when the sampling interval is an electrical angle of 30 °, V 2 = v (m) 2 + v (m− 3) 2 I 2 = i (m) 2 + i (m−3) 2 VI cos θ = v (m) · i (m) + v (m−3) · i (m−3) VI sin θ = v (m− 3) · i (m) −v (m) · i (m−3), impedance Z, resistance R, reactance X
Is calculated as follows: Z 2 = V 2 / I 2 R 2 = Z cos θ = (VI cos θ) / I 2 X 2 = Z sin θ = (VI sin θ) / 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 a certain period of time in the past, that is, within a time sufficient for the trajectory 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. Through the above arithmetic processing, it is possible to calculate the ratio of harmonic components included in the impedance and the order thereof from various amounts of electricity from the power system. Therefore, based on this, it is determined whether or not the point of occurrence of the specified accident is within the protection section, and it is determined whether or not power system protection is necessary.

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

【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 となる。
Here, looking at the form of this equation (2), it can be seen that the operation vector Idd is a combination of the fixed vector and the rotation vector. That is, I 1f
Represents a fixed vector, and Inf represents a vector that rotates in a counterclockwise direction, and this is expressed as R-X coordinates as shown in FIG. In FIG. 5, it can be seen that the operation vector Idd draws a locus which becomes maximum at the point 5a and becomes minimum at the point 5b. If 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 と表すことができる。
[0027] Here, since the I ddmax, I ddmin, and the time T ia that each vector has occurred, T ib may be determined from the sampling data of the captured current signal from the power system, I 1f , Inf can be calculated back. That is, the magnitude of I 1f is
dda, when the size of the I nf and I dds, can be expressed as I dda = (I ddmax + I ddmin) / 2 I dds = (I ddmax -I ddmin) / 2.

【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 harmonic contained in the operation vector Idd to the fundamental wave and its order n are represented by f =
The travel time between 5a and 5b based on the 50 Hz fundamental wave is represented by T
When iab, 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 the harmonic component 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. Next, in step 42, an operation vector is calculated by the above-described vector calculation method based on this data and the data which has been previously captured 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 a fixed time in the past, that is, within a time sufficient for the trajectory of the operation vector to pass through the relay operation area from the occurrence of the accident, The value and the time of occurrence are stored in steps 45 and 46, respectively. In step 47, the ratio of the harmonic component and its order are calculated from these values based on the above-described calculation method.

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

【0031】[0031]

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

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

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

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

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

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

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

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

【図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 a lower harmonic component in a conventional digital relay.

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

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

1 ディジタルリレー 2 アナログ/ディジタル変換部 7 ディジタル演算制御部 8 演算処理部 9 ROM 10 RAM 11 外部インターフェース 12 電力系統 13 変流器 DESCRIPTION OF SYMBOLS 1 Digital relay 2 Analog / digital conversion part 7 Digital operation control part 8 Operation processing part 9 ROM 10 RAM 11 External interface 12 Power system 13 Current transformer

───────────────────────────────────────────────────── フロントページの続き (58)調査した分野(Int.Cl.7,DB名) H02H 3/38 - 3/52 H02H 7/04 - 7/045 ──────────────────────────────────────────────────続 き Continued on the front page (58) Field surveyed (Int.Cl. 7 , DB name) H02H 3/38-3/52 H02H 7 /04-7/045

Claims (2)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】 電力系統からサンプリングした各種電気
量に対して演算を行うとともに、その演算結果に応じて
動作の要否を判定する保護継電器において、 サンプリングされた電圧および電流データに基づき、電
圧ベクトルを電流ベクトルで除算したインピーダンスベ
クトルを算出するベクトル演算手段と、 前記インピーダンスベクトルの最大値および最小値なら
びにこれら発生時刻をそれぞれ抽出する抽出手段と、この抽出手段により抽出されたインピーダンスベクトル
の最大値および最小値からインピーダンスベクトルに含
まれる高調波分の割合を算出するとともに、インピーダ
ンスベクトルの最大値および最小値の発生時刻から前記
高調波分の次数を算出する 高調波分算出手段とを備える
ことを特徴とする保護継電器。
1. A protection relay that performs an operation on various amounts of electricity sampled from a power system and determines whether or not operation is necessary according to a result of the operation. , A vector calculating means for calculating an impedance vector obtained by dividing a current vector by a current vector; an extracting means for extracting a maximum value and a minimum value of the impedance vector and their occurrence times ; and an impedance vector extracted by the extracting means.
From the maximum and minimum values of
Calculate the proportion of harmonics
From the occurrence times of the maximum and minimum values of the
And a harmonic component calculating means for calculating the order of the harmonic component.
【請求項2】 電力系統からサンプリングした各種電気
量に対して演算を行うとともに、その演算結果に応じて
動作の要否を判定する保護継電器において、 サンプリングされた電流データに基づき、電流ベクトル
を電力系統の基本波基準ベクトルで除算した演算ベクト
ルを算出するベクトル演算手段と、 前記演算ベクトルの最大値および最小値ならびにこれら
発生時刻をそれぞれ抽出する抽出手段と、この抽出手段により抽出された演算ベクトルの最大値お
よび最小値から演算ベクトルに含まれる高調波分の割合
を算出するとともに、演算ベクトルの最大値および最小
値の発生時刻から前記高調波分の次数を算出する 高調波
分算出手段とを備えることを特徴とする保護継電器。
2. A protection relay that performs an operation on various amounts of electricity sampled from an electric power system and determines whether or not operation is necessary according to the operation result, wherein a current vector is calculated based on the sampled current data. a vector calculating means for calculating an operation vector obtained by dividing the fundamental wave reference vector of the system, the maximum and minimum values of the operational vector as well as their
Extraction means for extracting each occurrence time, and the maximum value and the like of the operation vector extracted by the extraction means;
And the ratio of harmonic components included in the operation vector from the minimum value
And the maximum and minimum values of the operation vector
And a harmonic component calculating means for calculating an order of the harmonic from a value generation time .
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 JPH06121455A (en) 1994-04-28
JP3003421B2 true 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)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
PL1820034T3 (en) * 2004-11-18 2010-03-31 Powersense As Compensation of simple fiberoptic faraday effect sensors
US7920903B2 (en) 2007-01-04 2011-04-05 Bose Corporation Microphone techniques
CN105226607B (en) * 2015-10-22 2017-09-12 许继电气股份有限公司 A kind of extra-high voltage direct-current layer-specific access system 100Hz harmonic protection methods

Family Cites Families (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

Also Published As

Publication number Publication date
JPH06121455A (en) 1994-04-28

Similar Documents

Publication Publication Date Title
US3984737A (en) Protective relaying system
Bo et al. A new technique for transformer protection based on transient detection
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
JP3003421B2 (en) Protective relay
TWI222257B (en) Distance relay apparatus
US4450497A (en) Ultra-high-speed 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
CN104062555A (en) Identification method for distributing line high-resistance ground fault feature harmonic waves
CN112986817B (en) Brushless exciter exciting winding turn-to-turn short circuit fault on-line detection device and method
Kennedy et al. Variable digital filter response time in a digital distance relay
Panova et al. Development of mathematical models of microprocessor-based relay protection devices for 220/110 kV nodal distribution substation in Matlab/Simulink
EP1311866A1 (en) Current compensation method and device for power system protection
JPH1141793A (en) Exciting inrush current discriminator
JP3160612B2 (en) Power system insulation deterioration detection method and apparatus
Moravej et al. Distance protection closed-loop testing using RTDS
Fakruddin et al. Microprocessor-based algorithm for symmetrical-components calculation
JP3255691B2 (en) Negative-phase overcurrent relay
CA1144595A (en) Rotor thermal protection relay
JPH0297223A (en) Trouble zone sectional detector
Samonto et al. Fuzzy Based Algorithm for Stage-I of Cascaded Intelligent Relaying
JPH06292322A (en) Overcurrent relay
JPH02234071A (en) Ground fault detector
Padam et al. Designing A Protective Relay for Differential Protection of Power Transformers Using Clarks Transform and S Transform

Legal Events

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