CN107271842A - A kind of parallel erected on same tower double-circuit line Fault Locating Method based on positive-sequence component in the same direction - Google Patents
A kind of parallel erected on same tower double-circuit line Fault Locating Method based on positive-sequence component in the same direction Download PDFInfo
- Publication number
- CN107271842A CN107271842A CN201710413795.2A CN201710413795A CN107271842A CN 107271842 A CN107271842 A CN 107271842A CN 201710413795 A CN201710413795 A CN 201710413795A CN 107271842 A CN107271842 A CN 107271842A
- Authority
- CN
- China
- Prior art keywords
- same
- line
- fault
- sequence component
- sequence
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000009434 installation Methods 0.000 claims abstract description 5
- 238000005259 measurement Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 claims description 2
- 238000012546 transfer Methods 0.000 claims description 2
- 230000001681 protective effect Effects 0.000 claims 1
- 230000004807 localization Effects 0.000 abstract 1
- 230000005540 biological transmission Effects 0.000 description 9
- 238000000691 measurement method Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000014509 gene expression Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
- G01R31/085—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution lines, e.g. overhead
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Locating Faults (AREA)
Abstract
The present invention relates to a kind of parallel erected on same tower double-circuit line Fault Locating Method based on positive-sequence component in the same direction; it is characterized in that; the six-phase voltage amount and the magnitude of current at parallel erected on same tower double-circuit line system two ends after occurring first with the electric power mutual-inductor collection failure of protection installation place; then six sequence decouplings are carried out to the fundametal compoment at double-circuit line two ends using six-sequence component; fault localization function finally is constructed using positive-sequence component in the same direction, fault location is carried out according to the phase characteristic of range function.
Description
Technical Field
The invention relates to the technical field of power system relay protection, in particular to a same-tower parallel-frame double-circuit line fault positioning method based on a same-direction positive sequence component.
Background
With the development of power technology and the continuous expansion of power grid capacity, long-distance large-capacity transregional power transmission becomes the development direction of electric energy transmission. The double-circuit line parallel to the same tower is greatly applied to actual operation and planning construction and becomes the development trend of ultra-high voltage transmission due to the advantages of large transmission capacity, low engineering cost, narrow outgoing line corridor, small occupied area, short construction period and the like. Due to the wide application of the double circuit lines on the same tower and the requirement of system safety and stability, the research on corresponding fault location has also become a hot problem of people's attention. After a fault occurs, quick and accurate fault location has very important significance for timely repairing a line, reducing economic loss, recovering power supply as soon as possible and improving the reliability of a power system.
Because the distance between the two circuits of the same tower and the double circuit lines is very close, the two circuits not only have mutual inductance between phases, but also have mutual inductance between the lines. Besides single-circuit faults of any circuit, overline faults caused by factors such as lightning and pole falling can also occur. The problems of serious coupling between lines, various fault types, complexity and the like bring certain difficulty to fault location of the double-circuit line which is erected on the same tower, and the fault location method of the single-circuit line cannot be simply applied to the double-circuit line.
The existing double-circuit line fault distance measurement algorithm mainly comprises the following steps according to the principle: an intelligent method, a traveling wave distance measurement method and a fault analysis method. The artificial intelligence method utilizes the advanced results of related disciplines to carry out fault location, but is not perfect and mature in principle and is still in an exploration stage at present. The traveling wave distance measurement method utilizes the transmission property of fault transient traveling waves to measure the distance without being influenced by line types, grounding impedance, system oscillation and system parameters at two sides, but needs to be provided with special high-speed sampling equipment, and has large hardware investment and cost and complex technology. The fault analysis method utilizes system related parameters and voltage and current measured during line fault to obtain fault distance through analysis and calculation, and can be divided into a single-end measurement method and a double-end measurement method according to different information sources. The single-end measurement method utilizes the single-end electric quantity of the line to measure the distance, does not need communication equipment to transmit opposite-end information, but is difficult to eliminate the influence of system impedance and transition resistance on an algorithm, and the distance measurement precision cannot meet the requirement. The double-end method overcomes the defects in principle, but the calculation process is usually complex, iterative search and pseudo-root identification are required, and the calculation amount is large.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides a fault positioning method for a double-circuit line with a same tower and a same frame based on a same-direction positive sequence component.
A same-tower parallel-frame double-circuit line fault positioning method based on a same-direction positive sequence component is characterized in that a power transformer at a protection installation position is used for collecting six-phase voltage quantity and current quantity at two ends of a same-tower parallel-frame double-circuit line system after a fault occurs, then a six-sequence component method is adopted for conducting six-sequence decoupling on fundamental wave components at two ends of a double-circuit line, finally a fault distance measurement function is constructed by using the same-direction positive sequence component, and fault positioning is conducted according to phase characteristics of the distance measurement function. The method comprises the following specific steps:
(1) firstly, a power transformer at a protection installation position is used for collecting six-phase voltage quantities at two ends of a double-circuit line system which is erected on the same tower in parallel after a fault occursAnd the amount of six-phase currentWherein I represents a first loop, II represents a second loop, and phi is phase AOr phase B or phase C;
(2) converting the six-phase electrical quantities at two ends of the double-circuit line into six independent sequence components by using an M matrix given by the following formula: a T1 homodromous positive sequence component, a T2 homodromous negative sequence component, a T0 homodromous zero sequence component, an F1 reverse positive sequence component, an F2 reverse negative sequence component, and an F0 reverse zero sequence component.
In the formula a=ej120°°
(3) Using the forward voltage magnitude of the N-terminal of the lineMagnitude of equidirectional positive sequence currentCalculating the same-direction positive sequence electric quantity of the M end of the lineThe calculation formula is as follows:
in the formula,for the line forward sequence transfer constant in the same direction,is the line co-directional positive sequence wave impedance, Z1Is a forward sequence impedance of unit length of the line, Y1For the same-direction positive-sequence admittance of unit length of the line,/MNIs the length of the double-loop line MN.
(4) Calculating a ranging functionWherein lMfThe distance from the fault point to the M terminal.
(5) The distance from the fault point to the M end is calculated by the method (4)Wherein β is the constant gamma of line syntropy positive sequence transmissionT1The imaginary part of (1), angle (·), is the phase function.
The invention has the beneficial effects that:
(1) the homodromous positive sequence network decomposed by the six-sequence component method is used for ranging, the influence of mutual inductance between parallel lines is avoided, and faults can be located without mutual inductance parameters between the lines; because the same-direction positive sequence component exists in the whole fault period of all fault types of the same tower parallel double-circuit line, the distance measurement algorithm can accurately position any fault.
(2) The linear phase characteristic of the ranging function is utilized to carry out ranging, a pseudo root does not appear in principle, and the ranging result is not influenced by factors such as transition resistance, system impedance, load current and the like.
(3) The principle is simple, the fault position is directly obtained according to the phase characteristics of the ranging function, iterative search and pseudo root identification are not needed, the required operand is small, and the ranging precision is high.
Drawings
FIG. 1 is a schematic diagram of a double circuit line on the same tower
FIG. 2 is a schematic diagram of a double-circuit forward network
The meaning of each reference number in the drawings and the text:
is the potential of the m-terminal power supply,is n terminal power supply potential;
is the same-direction positive sequence voltage quantity of the line M side;
the amount of forward sequence current in the same direction of the line M side;
the voltage quantity of the N side of the line is the same-direction positive sequence voltage quantity;
the quantity of the same-direction positive sequence current at the N side of the line;
for line M side flowThe magnitude of the equidirectional positive sequence current of the fault point;
the amount of equidirectional positive sequence current flowing to a fault point on the N side of the line;
the positive sequence voltage quantity of the fault point in the same direction;
is the same-direction positive-sequence component of the fault current.
Detailed Description
The invention will be further described in detail with reference to the drawings.
Fig. 1 is a schematic diagram of a double-circuit line with the same tower and the same frame applying the invention. Because the coupling between the double-circuit lines of the same tower parallel frame is serious, the double-circuit lines with the mutual inductance between phases and the mutual inductance between the lines are decoupled into six mutually independent sequence nets by adopting a six-sequence component method. The same-direction positive sequence component exists in all fault types of the double-loop, so the invention adopts the same-direction positive sequence network to carry out fault location, and fig. 2 shows the double-loop same-direction positive sequence network after the fault occurs.
When a fault occurs at the point f of the double-circuit line system which is erected on the same tower in parallel, according to a long-line equation, the forward sequence voltage in the same direction at the fault point f and the forward sequence current in the same direction flowing into the fault point at the N side are respectively as follows:
wherein,andrespectively equal positive sequence voltage quantity and equal positive sequence current quantity of the N side of the circuit;the forward sequence voltage quantity of the fault point is the same;the amount of equidirectional positive sequence current flowing to a fault point on the N side of the line; gamma rayT1And ZcT1Respectively homodromous positive sequence transmission constant and wave impedance.
As can be seen from fig. 1:
whereinThe amount of forward sequence current flowing to the fault point on the line M side;is the same-direction positive-sequence component of the fault current.
The actual value formula of the homodromous positive sequence voltage at the end M is as follows:
substituting the formula (1) and the formula (2) into the formula (3) to obtain:
wherein,the actual equidirectional positive sequence voltage quantity of the end M of the line is obtained;the voltage quantity of the M end homodromous positive sequence is derived from the N end homodromous positive sequence component according to a long line equation; lMfAnd lNfThe distances from the M end and the N end of the line to a fault point, lMf+lNf=lMN。
From equation (3) and equation (4):
the actual value of the homodromous positive sequence current flowing from the line side to the M end is shown as (6):
substituting the formula (1) and the formula (2) into the formula (6) to obtain:
wherein,the actual equidirectional positive sequence current quantity flowing to the M end from the line side;the quantity of the M-end homodromous positive sequence current is derived from the N-end homodromous positive sequence component according to a long-line equation.
At the M end, there are:
substituting the formula (8) into the formula (7) and arranging to obtain:
using equation (5) and equation (9), the following functional expression representing the fault current can be written:
and constructing a fault distance measurement function according to the formula (5) and the formula (9), namely the formula (11).
Wherein lMfIs a positive number, and the number of the positive number,andrespectively, the same-direction positive sequence voltage quantity and the same-direction positive sequence current quantity on the line M side.Andthe electrical expressions of (a) are respectively:
the transmission constant of a high-voltage overhead line can be written asWherein r is0、L0、C0The resistance, inductance and capacitance values of the circuit per unit length are respectively, α is an attenuation constant, and β is a phase coefficient.
Using abs [ f (l)Mf)]And angle [ f (l)Mf)]Representing the magnitude and phase angle of the ranging function, respectively, the magnitude and phase angle of the ranging function can be further reduced to real number expressions of equation (12) and equation (13).
Since the line resistance of the high-voltage transmission line is very small, that is, the attenuation constant α is very small and close to zero, a certain error is caused by using the amplitude characteristic of the ranging function, that is, the formula (12) to perform ranging. In order to avoid measurement errors caused by transmission constants, the invention adopts the phase characteristic of the ranging function, namely formula (13), to carry out fault ranging, namely the phase angle value of the ranging function is in direct proportion to the fault distance.
The distance from the fault point to the M end is obtained by the method (13)
Claims (1)
1. A same-tower parallel-frame double-circuit line fault positioning method based on a same-direction positive sequence component is characterized in that a power transformer at a protection installation position is used for collecting six-phase voltage quantity and current quantity at two ends of a same-tower parallel-frame double-circuit line system after a fault occurs, then six-sequence decoupling is carried out on fundamental wave components at two ends of a double-circuit line by adopting a six-sequence component method, finally a fault distance measurement function is constructed by using the same-direction positive sequence component, and fault positioning is carried out according to the phase characteristics of the distance measurement function. The method comprises the following steps:
(1) firstly, the power transformer at the protective installation position is utilized to collect fault occurrenceSix-phase voltage quantity at two ends of rear same tower parallel double-circuit line systemAnd the amount of six-phase currentWherein I represents a first loop, II represents a second loop, and phi is phase A, phase B or phase C;
(2) converting the six-phase electrical quantities at two ends of the double-circuit line into six independent sequence components by using an M matrix given by the following formula: a T1 homodromous positive sequence component, a T2 homodromous negative sequence component, a T0 homodromous zero sequence component, an F1 reverse positive sequence component, an F2 reverse negative sequence component, an F0 reverse zero sequence component;
in the formula a=ej120°;
(3) Using the forward voltage magnitude of the N-terminal of the lineMagnitude of equidirectional positive sequence currentCalculating the equidirectional positive sequence electrical quantity of the M end of the lineThe calculation formula is as follows:
in the formula,for the line forward sequence transfer constant in the same direction,is the line co-directional positive sequence wave impedance, Z1Is a forward sequence impedance, Y, of the same direction per unit length of the line1For the same-direction positive-sequence admittance of unit length of the line,/MNThe length of the double-circuit line MN;
(4) calculating a ranging functionWherein lMfThe distance from the fault point to the M end;
(5) the distance from the fault point to the M end is calculated by the method (4)Wherein β is the constant gamma of line syntropy positive sequence transmissionT1The imaginary part of (1), angle (·), is the phase function.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710413795.2A CN107271842B (en) | 2017-06-05 | 2017-06-05 | same-tower parallel-frame double-circuit line fault positioning method based on same-direction positive sequence component |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710413795.2A CN107271842B (en) | 2017-06-05 | 2017-06-05 | same-tower parallel-frame double-circuit line fault positioning method based on same-direction positive sequence component |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107271842A true CN107271842A (en) | 2017-10-20 |
CN107271842B CN107271842B (en) | 2019-12-10 |
Family
ID=60065541
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710413795.2A Expired - Fee Related CN107271842B (en) | 2017-06-05 | 2017-06-05 | same-tower parallel-frame double-circuit line fault positioning method based on same-direction positive sequence component |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107271842B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108092245A (en) * | 2018-01-09 | 2018-05-29 | 国网江苏省电力有限公司无锡供电分公司 | A kind of high voltage transmission line line protection device and method |
CN108362978A (en) * | 2018-03-08 | 2018-08-03 | 国网山东省电力公司滨州供电公司 | A kind of T-type power transmission line fault locating method based on range function phase characteristic |
CN109655710A (en) * | 2018-12-03 | 2019-04-19 | 国网江苏省电力有限公司检修分公司 | Double circuits on same tower transmission line of electricity both ends measure the bearing calibration of phasor, apparatus and system |
CN109839569A (en) * | 2019-02-27 | 2019-06-04 | 燕山大学 | A kind of double circuit lines distance measuring method and system based on six-sequence component |
CN110361632A (en) * | 2019-08-05 | 2019-10-22 | 大唐水电科学技术研究院有限公司 | A kind of part coupling asynchronous fault distance-finding method of double-circuit line on same pole |
CN112363009A (en) * | 2020-10-14 | 2021-02-12 | 中国电力科学研究院有限公司 | Single-ended fault location method and system for same-tower line grounding fault |
CN112684374A (en) * | 2020-11-11 | 2021-04-20 | 大唐水电科学技术研究院有限公司 | Fault analysis method for same-phase overline disconnection of four circuit lines on same tower |
CN115166411A (en) * | 2022-06-20 | 2022-10-11 | 天津大学 | Single-circuit T-shaped power transmission line fault location method based on fault branch judgment |
CN115267641A (en) * | 2022-09-22 | 2022-11-01 | 武汉格蓝若智能技术有限公司 | Method and system for identifying error abnormity of current transformer in same-tower double-circuit power transmission line |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2752363A1 (en) * | 2010-09-16 | 2012-03-16 | Schweitzer Engineering Laboratories, Inc. | Fault location in a non-homogeneous electric power line |
CN102914726A (en) * | 2012-11-07 | 2013-02-06 | 华北电力大学(保定) | Fault positioning method for common-tower double-circuit line |
CN105929302A (en) * | 2016-04-13 | 2016-09-07 | 上海交通大学 | Sequence component relation based range finding method for single-end fault of power transmission line |
-
2017
- 2017-06-05 CN CN201710413795.2A patent/CN107271842B/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2752363A1 (en) * | 2010-09-16 | 2012-03-16 | Schweitzer Engineering Laboratories, Inc. | Fault location in a non-homogeneous electric power line |
CN102914726A (en) * | 2012-11-07 | 2013-02-06 | 华北电力大学(保定) | Fault positioning method for common-tower double-circuit line |
CN105929302A (en) * | 2016-04-13 | 2016-09-07 | 上海交通大学 | Sequence component relation based range finding method for single-end fault of power transmission line |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108092245B (en) * | 2018-01-09 | 2023-07-21 | 国网江苏省电力有限公司无锡供电分公司 | High-voltage transmission line protection device and method |
CN108092245A (en) * | 2018-01-09 | 2018-05-29 | 国网江苏省电力有限公司无锡供电分公司 | A kind of high voltage transmission line line protection device and method |
CN108362978A (en) * | 2018-03-08 | 2018-08-03 | 国网山东省电力公司滨州供电公司 | A kind of T-type power transmission line fault locating method based on range function phase characteristic |
CN109655710A (en) * | 2018-12-03 | 2019-04-19 | 国网江苏省电力有限公司检修分公司 | Double circuits on same tower transmission line of electricity both ends measure the bearing calibration of phasor, apparatus and system |
CN109655710B (en) * | 2018-12-03 | 2020-12-01 | 国网江苏省电力有限公司检修分公司 | Correction method, device and system for phasor measurement at two ends of same-pole double-circuit transmission line |
CN109839569A (en) * | 2019-02-27 | 2019-06-04 | 燕山大学 | A kind of double circuit lines distance measuring method and system based on six-sequence component |
CN110361632A (en) * | 2019-08-05 | 2019-10-22 | 大唐水电科学技术研究院有限公司 | A kind of part coupling asynchronous fault distance-finding method of double-circuit line on same pole |
CN112363009A (en) * | 2020-10-14 | 2021-02-12 | 中国电力科学研究院有限公司 | Single-ended fault location method and system for same-tower line grounding fault |
CN112363009B (en) * | 2020-10-14 | 2023-12-01 | 中国电力科学研究院有限公司 | Single-ended fault location method and system for same-tower line ground fault |
CN112684374A (en) * | 2020-11-11 | 2021-04-20 | 大唐水电科学技术研究院有限公司 | Fault analysis method for same-phase overline disconnection of four circuit lines on same tower |
CN115166411A (en) * | 2022-06-20 | 2022-10-11 | 天津大学 | Single-circuit T-shaped power transmission line fault location method based on fault branch judgment |
CN115267641A (en) * | 2022-09-22 | 2022-11-01 | 武汉格蓝若智能技术有限公司 | Method and system for identifying error abnormity of current transformer in same-tower double-circuit power transmission line |
CN115267641B (en) * | 2022-09-22 | 2023-01-17 | 武汉格蓝若智能技术有限公司 | Method and system for identifying error abnormity of current transformer in same-tower double-circuit power transmission line |
Also Published As
Publication number | Publication date |
---|---|
CN107271842B (en) | 2019-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107271842B (en) | same-tower parallel-frame double-circuit line fault positioning method based on same-direction positive sequence component | |
CN105067950B (en) | Two Terminal Fault Location method based on longitudinal impedance | |
CN102200563B (en) | Line single-phase earth fault single-terminal ranging method based on positioning function amplitude characteristics | |
Suonan et al. | A novel fault-location method for HVDC transmission lines | |
CN102135571B (en) | Anti-interference measurement method for zero sequence impedance of super-high-voltage/ultrahigh-voltage multi-loop power transmission line | |
CN105738769B (en) | Series compensation double line down localization method based on distributed parameter model | |
CN106124927B (en) | The asynchronous fault distance-finding method of double-circuit lines on the same pole containing series compensation | |
CN103792465A (en) | Power distribution network one-phase grounding fault location method based on zero sequence voltage | |
CN110850236A (en) | Power distribution network fault positioning method based on parameter estimation | |
CN102175954A (en) | Circuit inter-phase fault single-end ranging method | |
CN109541392A (en) | A kind of one-end fault ranging method suitable for flexible HVDC transmission system | |
CN107632238B (en) | Multi-end transmission line fault location method based on WAMS system | |
CN113671314B (en) | Method for positioning and ranging single-phase earth fault section of ring network of power distribution network | |
CN107817420B (en) | Non-synchronous data fault location method for non-whole-course same-tower double-circuit power transmission line | |
Wang et al. | Novel location algorithm for single‐line‐to‐ground faults in transmission line with distributed parameters | |
CN107271841A (en) | The T-shaped fault positioning method for transmission line of common-tower double-return based on forward-order current difference | |
CN108362978A (en) | A kind of T-type power transmission line fault locating method based on range function phase characteristic | |
Chen et al. | A novel traveling wave based fault location scheme for power distribution grids with distributed generations | |
CN103217623B (en) | The line double-end fault distance-finding method had nothing to do with transition resistance and load current | |
Nasab et al. | A hybrid scheme for fault locating for transmission lines with TCSC | |
Lin et al. | Fault location for three-ended ring-topology power system using minimum GPS-based measurements and CVT/CT sensing | |
Li et al. | Augmented state estimation method for fault location based on on-line parameter identification of PMU measurement data | |
Chen et al. | A new scheme for fault location of three‐terminal parallel transmission lines without transposer | |
CN110456229A (en) | A kind of Single Terminal Traveling Wave Fault Location device of distribution line | |
CN106199192A (en) | The positive sequence on-line testing method of parallel many back transmission lines Zero sequence parameter |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191210 Termination date: 20210605 |
|
CF01 | Termination of patent right due to non-payment of annual fee |