CN116169646A - Bus protection method and system based on current sampling value - Google Patents

Bus protection method and system based on current sampling value Download PDF

Info

Publication number
CN116169646A
CN116169646A CN202211137833.3A CN202211137833A CN116169646A CN 116169646 A CN116169646 A CN 116169646A CN 202211137833 A CN202211137833 A CN 202211137833A CN 116169646 A CN116169646 A CN 116169646A
Authority
CN
China
Prior art keywords
current
value
sampling
sudden change
preset
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.)
Pending
Application number
CN202211137833.3A
Other languages
Chinese (zh)
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.)
China Electric Power Research Institute Co Ltd CEPRI
Original Assignee
China Electric Power Research Institute Co Ltd CEPRI
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 China Electric Power Research Institute Co Ltd CEPRI filed Critical China Electric Power Research Institute Co Ltd CEPRI
Priority to CN202211137833.3A priority Critical patent/CN116169646A/en
Publication of CN116169646A publication Critical patent/CN116169646A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02H—EMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
    • H02H7/00—Emergency protective circuit arrangements specially adapted for specific types of electric machines or apparatus or for sectionalised protection of cable or line systems, and effecting automatic switching in the event of an undesired change from normal working conditions
    • H02H7/22—Emergency protective circuit arrangements specially adapted for specific types of electric machines or apparatus or for sectionalised protection of cable or line systems, and effecting automatic switching in the event of an undesired change from normal working conditions for distribution gear, e.g. bus-bar systems; for switching devices
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
    • G01R19/0038—Circuits for comparing several input signals and for indicating the result of this comparison, e.g. equal, different, greater, smaller (comparing pulses or pulse trains according to amplitude)
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
    • G01R19/10—Measuring sum, difference or ratio
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
    • G01R19/14—Indicating direction of current; Indicating polarity of voltage
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
    • G01R19/30—Measuring the maximum or the minimum value of current or voltage reached in a time interval
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02H—EMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
    • H02H1/00—Details of emergency protective circuit arrangements
    • H02H1/0007—Details of emergency protective circuit arrangements concerning the detecting means

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Emergency Protection Circuit Devices (AREA)

Abstract

The invention discloses a bus protection method and a bus protection system based on a current sampling value, wherein the bus protection method comprises the following steps: when the starting element meets a preset starting criterion and is started, current grouping is carried out on the basis of current abrupt change value sampling values of each side of the electric element at each sampling moment, and a first current abrupt change value and a second current abrupt change value corresponding to each sampling moment are determined; performing polarity conversion on the first current abrupt change amount and the second current abrupt change amount corresponding to each sampling time to obtain a first current abrupt change amount conversion value and a second current abrupt change amount conversion value corresponding to each sampling time; performing point-by-point integration based on the first current abrupt change conversion value and the second current abrupt change conversion value corresponding to each sampling moment to obtain a first integral value and a second integral value; and determining the fault type according to a preset bus protection criterion based on the first integral value and the second integral value, so as to perform bus protection based on the determined fault type.

Description

Bus protection method and system based on current sampling value
Technical Field
The invention relates to the technical field of relay protection of power systems, in particular to a bus protection method based on a current sampling value.
Background
And a bus in the transformer substation is connected with a circuit and a transformer, so that the influence range of bus faults is large, and the fault bus needs to be rapidly cut off by bus protection. The existing bus protection mostly adopts the principle of industrial frequency differential motion protection, and the action speed of protection is limited by a fixed sampling data window. In recent years, a large number of new energy sources are connected into an alternating current power grid, after an alternating current system fails, short circuit current provided by the new energy sources is influenced by a control strategy, the amplitude of the failure current is rapidly reduced, and the rapidity and the sensitivity of bus protection are seriously influenced.
In addition, the protection principle based on transient state quantity is mainly used for fault identification by adopting the wave head characteristics of traveling waves or transient state energy of different frequency bands on a frequency domain, but the sampling rate is high, the calculation amount of a mathematical algorithm for extracting the transient state component is large, the protection action performance is limited by various factors such as fault time, transition resistance and the like, and the protection reliability is required to be further improved.
Disclosure of Invention
The invention provides a bus protection method and system based on a current sampling value, which are used for solving the problem of how to quickly and accurately identify bus faults.
In order to solve the above problems, according to an aspect of the present invention, there is provided a bus bar protection method based on a current sampling value, the method comprising:
When the starting element meets a preset starting criterion and is started, current grouping is carried out on the basis of current abrupt change value sampling values of each side of the electric element at each sampling moment, and a first current abrupt change value and a second current abrupt change value corresponding to each sampling moment are determined;
performing polarity conversion on the first current abrupt change amount and the second current abrupt change amount corresponding to each sampling time to obtain a first current abrupt change amount conversion value and a second current abrupt change amount conversion value corresponding to each sampling time;
performing point-by-point integration based on the first current abrupt change conversion value and the second current abrupt change conversion value corresponding to each sampling moment to obtain a first integral value and a second integral value;
and determining the fault type according to a preset bus protection criterion based on the first integral value and the second integral value, so as to perform bus protection based on the determined fault type.
Preferably, the determining the first current mutation amount and the second current mutation amount corresponding to each sampling time based on the current mutation amount sampling values of each side of the electrical element at each sampling time includes:
for any sampling time, selecting the current abrupt change sampling value with the largest amplitude in the current abrupt change sampling values of all sides of the electric element at the any sampling time as the first current abrupt change delta i corresponding to the any sampling time m (t)＝max(|Δi 1 (t)|,|Δi 2 (t)|,…,|Δi h (t) |); determining a second current abrupt change delta i corresponding to any sampling moment according to the difference value between the sum of the current abrupt change sampling values of each side of the electric element at any sampling moment and the first current abrupt change n (t)＝Δi Σ (t)-Δi m (t)；
Wherein, |Δi h (t) is the amplitude of the current abrupt change sampling value at the sampling time t at the h side; Δi Σ And (t) is the sum of current mutation sampling values of each side of the electric element at the sampling time t.
Preferably, the performing polarity conversion on the first current abrupt change amount and the second current abrupt change amount corresponding to each sampling time to obtain a first current abrupt change amount conversion value and a second current abrupt change amount conversion value corresponding to each sampling time includes:
for any sampling time, if the first current abrupt change delta i corresponding to any sampling time t m If the polarity of (t) is negative, determining Δi M (t)＝-Δi m (t)，Δi N (t)＝-Δi n (t); otherwise, determining Δi M (t)＝Δi m (t),Δi N (t)＝Δi n (t); wherein Δi n (t) is the sampling timet corresponds to a second current abrupt change; Δi M (t) and Δi N (t) the first current abrupt change value and the second current abrupt change value, respectively.
Preferably, the determining the fault type according to a preset bus protection criterion based on the first integral value and the second integral value includes:
If the first current abrupt change conversion value and the second current abrupt change conversion value meet a preset bus protection criterion, determining that the fault type is a fault in a bus area; otherwise, determining the fault type as the fault outside the bus area;
the bus protection criterion is preset, and the bus protection criterion comprises:
Figure BDA0003852112740000031
wherein ΣΔi M (t) and ΣΔi N (t) a first current-abrupt change-over value and a second current-abrupt change-over value, respectively; i e Is the rated current of the bus bar protection device.
Preferably, wherein the method further comprises:
when the starting element meets the first preset promoter criterion or the second preset promoter criterion, determining that the starting element meets the preset starting criterion;
wherein the first preset promoter criterion comprises:
Δf 1 (t)＝f 1 (t)-f 1 (t-T/4)>f 1.set ，
the second preset promoter criterion comprises:
Δf 2 (t)＝f 2 (t)-f 2 (t-T/4)>f 2.set ，
wherein Δf 1 (t) is a first variation; Δf 2 (t) is a second variation; f (f) 1 (t) is first data at time t; f (f) 2 (t) is second data at time t; f (f) 1 (T-T/4) is f before quarter cycle 1 (t)；f 2 (T-T/4) is f before quarter cycle 2 (t)；f 1.set A first preset threshold value; f (f) 2.set Is a second predetermined threshold.
Preferably, wherein the method further comprises:
calculating abrupt change delta i of each branch current sampling value of bus based on each branch current sampling value of bus φ.k (t) comprising:
Δi φ.k (t)＝i φ.k (t)-i φ.k (t-T s )，
calculating differential current i of each phase based on mutation amount of current sampling value of each branch φ.Σ (t) comprising:
Figure BDA0003852112740000032
calculating first data based on each phase difference stream, comprising:
f 1 (t)＝(i a.Σ (t)-i b.Σ (t)) 2 +(i b.Σ (t)-i c.Σ (t)) 2 +(i c.Σ (t)-i a.Σ (t) 2 ，
wherein i is φ.k (t) is the current sampling value of phi phase k branch at t time, phi is A, B, C, i φ.k (t-T s ) T-T being the kth branch of phi phase s Current sample value at moment, T s The time corresponding to the one-cycle wave; n is the number of branches; f (f) 1 (t) is first data; i.e a.Σ (t)、i b.Σ (t) and i c.Σ (t) is the differential flow of three phases A, B, C at time t respectively;
dividing the current of each branch of the bus into two groups, comparing the mutation amounts of the current sampling values of all branches, and determining the branch current delta i with the maximum mutation amount of the current sampling value φ.max (t)，Δi φ.max (t)＝max(Δi φ.1 (t),Δi φ.2 (t)，…，Δi φ.k (t))；
Calculate the division Δi φ.max (t) the sum Δi of the abrupt amounts of the current sampling values of the other branches outside the corresponding branch φ.s (t)，Δi φ.s (t)＝i φ.Σ (t)-Δi φ.max (t)；
Based on i φ.s (t) calculating second data comprises:
f 2 (t)＝(i a.s (t)-i b.s (t)) 2 +(i b.s (t)-i c.s (t)) 2 +(i c.s (t)-i a.s (t)) 2 ，
wherein Δi φ.k (t) is a current sampling value of the kth branch of phi phase at the moment t, and phi is A, B, C; i.e φ.Σ (t) is the sum of the current sample values of all branches of phi phase; f (f) 2 (t) is second data.
According to another aspect of the present invention, there is provided a bus bar protection system based on current sampling values, the system comprising:
the current abrupt change determining unit is used for grouping currents based on current abrupt change sampling values of each side of the electric element at each sampling moment when the starting element meets a preset starting criterion and determining a first current abrupt change and a second current abrupt change corresponding to each sampling moment;
The polarity conversion unit is used for respectively carrying out polarity conversion on the first current abrupt change amount and the second current abrupt change amount corresponding to each sampling time so as to obtain a first current abrupt change amount conversion value and a second current abrupt change amount conversion value corresponding to each sampling time;
the integration unit is used for carrying out point-to-point integration based on the first current mutation amount conversion value and the second current mutation amount conversion value corresponding to each sampling time to obtain a first integration value and a second integration value;
and the bus protection unit is used for determining the fault type according to a preset bus protection criterion based on the first integral value and the second integral value so as to conduct bus protection based on the determined fault type.
Preferably, the current abrupt change amount determining unit performs current grouping based on current abrupt change amount sampling values of each side of the electric element at each sampling time, and determines a first current abrupt change amount and a second current abrupt change amount corresponding to each sampling time, including:
for any sampling time, each side of the electrical element is selected at the sampling timeThe current abrupt change sampling value with the largest amplitude in the current abrupt change sampling values is the first current abrupt change delta i corresponding to any sampling time m (t)＝max(|Δi 1 (t)|,|Δi 2 (t)|,…,|Δi h (t) |); determining a second current abrupt change delta i corresponding to any sampling moment according to the difference value between the sum of the current abrupt change sampling values of each side of the electric element at any sampling moment and the first current abrupt change n (t)＝Δi Σ (t)-Δi m (t)；
Wherein, |Δi h (t) is the amplitude of the current abrupt change sampling value at the sampling time t at the h side; Δi Σ And (t) is the sum of current mutation sampling values of each side of the electric element at the sampling time t.
Preferably, the polarity conversion unit performs polarity conversion on the first current abrupt change amount and the second current abrupt change amount corresponding to each sampling time, so as to obtain a first current abrupt change amount conversion value and a second current abrupt change amount conversion value corresponding to each sampling time, and the polarity conversion unit includes:
for any sampling time, if the first current abrupt change delta i corresponding to any sampling time t m If the polarity of (t) is negative, determining Δi M (t)＝-Δi m (t)，Δi N (t)＝-Δi n (t); otherwise, determining Δi M (t)＝Δi m (t),Δi N (t)＝Δi n (t); wherein Δi n (t) is a second current abrupt change amount corresponding to the sampling time t; Δi M (t) and Δi N (t) the first current abrupt change value and the second current abrupt change value, respectively.
Preferably, the bus protection unit determines the fault type according to a preset bus protection criterion based on the first integrated value and the second integrated value, and includes:
If the first current abrupt change conversion value and the second current abrupt change conversion value meet a preset bus protection criterion, determining that the fault type is a fault in a bus area; otherwise, determining the fault type as the fault outside the bus area;
the bus protection criterion is preset, and the bus protection criterion comprises:
Figure BDA0003852112740000051
wherein ΣΔi M (t) and ΣΔi N (t) a first current-abrupt change-over value and a second current-abrupt change-over value, respectively; i e Is the rated current of the bus bar protection device.
Preferably, wherein the system further comprises:
the starting judgment unit is used for determining that the starting element meets the preset starting criterion when the starting element meets the first preset promoter criterion or the second preset promoter criterion;
wherein the first preset promoter criterion comprises:
Δf 1 (t)＝f 1 (t)-f 1 (t-T/4)>f 1.set ，
the second preset promoter criterion comprises:
Δf 2 (t)＝f 2 (t)-f 2 (t-T/4)>f 2.set ，
wherein Δf 1 (t) is a first variation; Δf 2 (t) is a second variation; f (f) 1 (t) is first data at time t; f (f) 2 (t) is second data at time t; f (f) 1 (T-T/4) is f before quarter cycle 1 (t)；f 2 (T-T/4) is f before quarter cycle 2 (t)；f 1.set A first preset threshold value; f (f) 2.set Is a second predetermined threshold.
Preferably, the start-up judging unit further includes:
Calculating abrupt change delta i of each branch current sampling value of bus based on each branch current sampling value of bus φ.k (t) comprising:
Δi φ.k (t)＝i φ.k (t)-i φ.k (t-T s )，
calculating differential current i of each phase based on mutation amount of current sampling value of each branch φ.Σ (t) comprising:
Figure BDA0003852112740000061
calculating first data based on each phase difference stream, comprising:
f 1 (t)＝(i a.Σ (t)-i b.Σ (t)) 2 +(i b.Σ (t)-i c.Σ (t)) 2 +(i c.Σ (t)-i a.Σ (t) 2 ，
wherein i is φ.k (t) is the current sampling value of phi phase k branch at t time, phi is A, B, C, i φ.k (t-T s ) T-T being the kth branch of phi phase s Current sample value at moment, T s The time corresponding to the one-cycle wave; n is the number of branches; f (f) 1 (t) is first data; i.e a.Σ (t)、i b.Σ (t) and i c . Σ (t) is the differential flow of three phases A, B, C at time t respectively;
dividing the current of each branch of the bus into two groups, comparing the mutation amounts of the current sampling values of all branches, and determining the branch current delta i with the maximum mutation amount of the current sampling value φ.max (t)，Δi φ.max (t)＝max(Δi φ . 1 (t),Δi φ.2 (t)，…，Δi φ.k (t))；
Calculate the division Δi φ.max (t) the sum Δi of the abrupt amounts of the current sampling values of the other branches outside the corresponding branch φ . s (t)，Δi φ.s (t)＝i φ.Σ (t)-Δi φ.max (t)；
Based on i φ.s (t) calculating second data comprises:
f 2 (t)＝(i a.s (t)-i b.s (t)) 2 +(i b.s (t)-i c.s (t)) 2 +(i c.s (t)-i a.s (t)) 2 ，
wherein Δi φ.k (t) is a current sampling value of the kth branch of phi phase at the moment t, and phi is A, B, C; i.e φ.Σ (t) is the sum of the current sample values of all branches of phi phase; f (f) 2 (t) is second data.
The invention provides a bus protection method and system based on a current sampling value, comprising the following steps: when the starting element meets a preset starting criterion and is started, current grouping is carried out on the basis of current abrupt change value sampling values of each side of the electric element at each sampling moment, and a first current abrupt change value and a second current abrupt change value corresponding to each sampling moment are determined; performing polarity conversion on the first current abrupt change amount and the second current abrupt change amount corresponding to each sampling time to obtain a first current abrupt change amount conversion value and a second current abrupt change amount conversion value corresponding to each sampling time; performing point-by-point integration based on the first current abrupt change conversion value and the second current abrupt change conversion value corresponding to each sampling moment to obtain a first integral value and a second integral value; and determining the fault type according to a preset bus protection criterion based on the first integral value and the second integral value, so as to perform bus protection based on the determined fault type. The starting element is responsible for detecting the occurrence time of faults, the protection element is responsible for identifying internal and external faults of a bus area, the sudden change characteristics of the current after the faults can be extracted rapidly by utilizing the current sampling values of three phases at the same time, the sudden change amount of the current sampling values of multiple branches of the bus and the starting element and the sudden change amount of the grouping current sampling values are constructed by utilizing the current sampling values of multiple branches of the bus, the occurrence time of the internal and external faults of the bus area is detected rapidly according to the action behaviors and time sequences of the two starting elements, the current of the multiple branches is divided into two groups according to the characteristics of the current sampling values of the multiple branches in the internal and external faults of the bus area, the motion trail of the current sampling values after the faults is visually reflected on the plane through the polarity conversion, the sampling value integration and other data processing, and the bus faults are identified rapidly and accurately through dividing the internal and external fault identification areas of the bus.
Drawings
Exemplary embodiments of the present invention may be more completely understood in consideration of the following drawings:
FIG. 1 is a flow chart of a current sample value based bus bar protection method 100 according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a bus bar in-zone fault and out-of-zone fault identification zone according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of a bus fault simulation model according to an embodiment of the present invention;
FIG. 4 is a schematic diagram of current waveforms for each branch of a fault bus in a zone according to an embodiment of the present invention;
FIG. 5 is a schematic diagram of the result of an action of a fault bus initiation element within a zone according to an embodiment of the present invention;
FIG. 6 is a schematic diagram of an in-zone fault bus protection action scenario in accordance with an embodiment of the present invention;
FIG. 7 is a schematic diagram of current waveforms for each branch of an out-of-zone fault bus according to an embodiment of the present invention;
FIG. 8 is a schematic diagram of the operation of an out-of-zone fault bus initiation element according to an embodiment of the present invention;
FIG. 9 is a schematic diagram of an out-of-zone fault bus protection action scenario in accordance with an embodiment of the present invention;
FIG. 10 is a schematic diagram of current waveforms of each branch of a bus at the time of an out-of-zone fault to an in-zone fault according to an embodiment of the present invention;
FIG. 11 is a schematic diagram of the operation of the out-of-zone fault to in-zone fault bus initiation element according to an embodiment of the present invention;
FIG. 12 is a schematic diagram of an out-of-zone fault to in-zone fault bus protection action scenario in accordance with an embodiment of the present invention;
fig. 13 is a schematic structural diagram of a bus bar protection system 1300 based on current sampling values according to an embodiment of the present invention.
Detailed Description
The exemplary embodiments of the present invention will now be described with reference to the accompanying drawings, however, the present invention may be embodied in many different forms and is not limited to the examples described herein, which are provided to fully and completely disclose the present invention and fully convey the scope of the invention to those skilled in the art. The terminology used in the exemplary embodiments illustrated in the accompanying drawings is not intended to be limiting of the invention. In the drawings, like elements/components are referred to by like reference numerals.
Unless otherwise indicated, terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art. In addition, it will be understood that terms defined in commonly used dictionaries should be interpreted as having a meaning that is consistent with their meaning in the context of the relevant art and will not be interpreted in an idealized or overly formal sense.
Fig. 1 is a flow chart of a bus bar protection method 100 based on current sample values according to an embodiment of the present invention. As shown in fig. 1, in the bus protection method based on current sampling values provided in the embodiment of the invention, a starting element is responsible for detecting the occurrence time of faults, a protection element is responsible for identifying internal faults and external faults of a bus area, the sudden change characteristics of current after faults can be rapidly extracted by using three-phase current sampling values at the same time, multiple-branch current sampling value sudden change amounts of the multiple-branch current sampling values and starting elements and grouping current sampling value sudden change starting elements are constructed by using bus multiple-branch current sampling values, the occurrence time of internal faults and external faults of the bus area are rapidly detected according to the action behaviors and time sequences of the two starting elements, then the multiple-branch current is divided into two groups according to the characteristics of the multiple-branch current sampling values in the internal faults and the external faults of the bus area, the motion trail of the current sampling values after faults is intuitively reflected on the plane through data processing such as polarity conversion and sampling value integration, and the bus fault identification area is rapidly and accurately identified through dividing the internal faults and external faults. Starting from step 101, when a starting element meets a preset starting criterion and starts in step 101, current grouping is performed based on current mutation value sampling values of each sampling time on each side of an electric element, and a first current mutation value and a second current mutation value corresponding to each sampling time are determined.
Preferably, wherein the method further comprises:
when the starting element meets the first preset promoter criterion or the second preset promoter criterion, determining that the starting element meets the preset starting criterion;
wherein the first preset promoter criterion comprises:
Δf 1 (t)＝f 1 (t)-f 1 (t-T/4)>f 1.set ，
the second preset promoter criterion comprises:
Δf 2 (t)＝f 2 (t)-f 2 (t-T/4)>f 2.set ，
wherein Δf 1 (t) is a first variation; Δf 2 (t) is a second variation; f (f) 1 (t) is first data at time t; f (f) 2 (t) is second data at time t; f (f) 1 (T-T/4) is f before quarter cycle 1 (t)；f 2 (T-T/4) is f before quarter cycle 2 (t)；f 1.set A first preset threshold value; f (f) 2.set Is a second predetermined threshold.
Preferably, wherein the method further comprises:
calculating abrupt change delta i of each branch current sampling value of bus based on each branch current sampling value of bus φ.k (t) comprising:
Δi φ.k (t)＝i φ.k (t)-i φ.k (t-T s )，
calculating differential current i of each phase based on mutation amount of current sampling value of each branch φ.Σ (t) comprising:
Figure BDA0003852112740000091
calculating first data based on each phase difference stream, comprising:
f 1 (t)＝(i a.Σ (t)-i b.Σ (t)) 2 +(i b.Σ (t)-i c.Σ (t)) 2 +(i c.Σ (t)-i a.Σ (t) 2 ，
wherein i is φ.k (t) is the current sampling value of phi phase k branch at t time, phi is A, B, C, i φ.k (t-T s ) T-T being the kth branch of phi phase s Current sample value at moment, T s The time corresponding to the one-cycle wave; n is the number of branches; f (f) 1 (t) is first data; i.e a.Σ (t)、i b.Σ (t) and i c . Σ (t) is the differential flow of three phases A, B, C at time t respectively;
dividing the current of each branch of the bus into two groups, comparing the mutation amounts of the current sampling values of all branches, and determining the branch current delta i with the maximum mutation amount of the current sampling value φ.max (t)，Δi φ.max (t)＝max(Δi φ.1 (t),Δi φ.2 (t)，…，Δi φ.k (t))；
Calculate the division Δi φ.max (t) the sum Δi of the abrupt amounts of the current sampling values of the other branches outside the corresponding branch φ.s (t)，Δi φ.s (t)＝i φ.Σ (t)-Δi φ.max (t)；
Based on i φ.s (t) calculating second data comprises:
f 2 (t)＝(i a.s (t)-i b.s (t)) 2 +(i b.s (t)-i c.s (t)) 2 +(i c.s (t)-i a.s (t)) 2 ，
wherein Δi φ.k (t) is a current sampling value of the kth branch of phi phase at the moment t, and phi is A, B, C; i.e φ.Σ (t) is the sum of the current sample values of all branches of phi phase; f (f) 2 (t) is second data.
Preferably, the determining the first current mutation amount and the second current mutation amount corresponding to each sampling time based on the current mutation amount sampling values of each side of the electrical element at each sampling time includes:
for any sampling time, selecting the current abrupt change sampling value with the largest amplitude in the current abrupt change sampling values of all sides of the electric element at the any sampling time as the first current abrupt change delta i corresponding to the any sampling time m (t)＝max(|Δi 1 (t)|,|Δi 2 (t)|,…,|Δi h (t) |); determining a second current abrupt change delta i corresponding to any sampling moment according to the difference value between the sum of the current abrupt change sampling values of each side of the electric element at any sampling moment and the first current abrupt change n (t)＝Δi Σ (t)-Δi m (t)；
Wherein, |Δi h (t) is the amplitude of the current abrupt change sampling value at the sampling time t at the h side; Δi Σ (t) sampling the electrical component on each sideThe sum of the current burst value samples at time t.
The invention utilizes the busbar protection of the current abrupt change sampling value track characteristic to realize based on a starting element and a protection element. The starting element comprises a multi-branch current sampling value abrupt change quantity and a starting element, and a grouping current sampling value abrupt change quantity starting element.
(1) Principle of abrupt change of sampling value and starting element of multi-branch current
And constructing a fault starting element in the area by utilizing the abrupt change of the current sampling value of each branch of the bus.
Firstly, calculating the abrupt change delta i of the current sampling value of each branch of the bus φ.k (t)，
Δi φ.k (t)＝i φ.k (t)-i φ.k (t-T s ) (1)
Wherein i is φ.k (t) is the current sampling value of phi phase k branch at t time, phi is A, B, C, i φ.k (t-T s ) T-T being the kth branch of phi phase s Current sample value at moment, T s Is the corresponding time of a cycle wave (20 ms).
Calculating each phase difference stream i by using mutation quantity of current sampling value of each branch φ.Σ (t)，
Figure BDA0003852112740000111
Will i φ.Σ (t) substituting f (t) to obtain
f 1 (t)＝(i a.Σ (t)-i b.Σ (t)) 2 +(i b.Σ (t)-i c.Σ (t)) 2 +(i c.Σ (t)-i a.Σ (t) 2 (3)
The abrupt change amount of the multi-branch current sampling value and the starting criterion I (namely a first preset promoter criterion) are as follows:
Δf 1 (t)＝f 1 (t)-f 1 (t-T/4)>f 1.set (4)
wherein f 1 (T-T/4) is f before one quarter cycle (5 ms) 1 (t)，f 1.set Is fixed for action, and f 1.set ＝f 1. Floating device +f 1. Fixing ，f 1. Floating device For floating threshold, the value is T E [ T-T, T-T/2 ]]Δf in between 1 (t) maximum value. f (f) 1. Fixing Is a fixed threshold.
(2) Principle of grouping current sampling value abrupt change starting element
Firstly, dividing the current of each branch of a bus into two groups, comparing the mutation amounts of current sampling values of all branches, and extracting the branch current delta i with the largest mutation amount of the current sampling values φ.max (t)。
Δi φ.max (t)＝max(Δi φ.1 (t),Δi φ.2 (t)…Δi φ.k (t)) (5)
Let Δi φ.max (t) corresponding branch currents are one group, and the mutation amounts of the sampling values of the rest branch currents are summed to form another group delta i φ.s (t) obtainable according to formula (2),
Δi φ.s (t)＝i φ.Σ (t)-Δi φ.max (t) (6)
will i φ.s (t) substituting f (t) to obtain
f 2 (t)＝(i a.s (t)-i b.s (t)) 2 +(i b.s (t)-i c.s (t)) 2 +(i c.s (t)-i a.s (t)) 2 (7)
The grouping current sampling value mutation quantity starting criterion II (namely a second preset promoter criterion) is as follows:
Δf 2 (t)＝f 2 (t)-f 2 (t-T/4)>f 2.set (8)
wherein f 2 (T-T/4) is f before one quarter cycle (5 ms) 2 (t)，f 2.set Is fixed for action, and f 2.set ＝f 2. Floating device +f 2. Fixing ，f 2. Floating device For floating threshold, the value is Deltaf 2 (T) at t.epsilon.t-T, T-T/2]Maximum value between. f (f) 2. Fixing Is a fixed threshold.
In the invention, when the abrupt change of the multi-branch current sampling value and the starting element meet a first preset starting criterion or the abrupt change of the grouping current sampling value meets a second preset starting criterion, the starting element is determined to meet the preset starting criterion, at the moment, the bus can be determined to have faults, but the fault type cannot be determined, at the moment, the starting element is started, and the protection element is triggered to execute bus protection, so that the fault type is determined.
Specifically, the current grouping is performed first. Specifically, comparing the sampled value amplitude of the current abrupt change at each sampling time of each side of the electric element, and selecting a maximum amplitude branch max (|Δi) 1 (t)|,|Δi 2 (t)|,…,|Δi h (t) |), let Δi m (t) sampling the current mutation of the branch with the maximum amplitude, and adding delta i to the current mutation m (t) as a first current argument, the remaining branch current sample value arguments are summed to Δi n (t)＝Δi Σ (t)-Δi m (t) adding Δi to the reaction mixture n (t) as a second current abrupt amount; wherein, |Δi h (t) is the amplitude of the current abrupt change sampling value at the sampling time t at the h side; Δi Σ And (t) is the sum of current mutation sampling values of each side of the electric element at the sampling time t.
In step 102, the polarity of the first current abrupt change amount and the second current abrupt change amount corresponding to each sampling time are respectively converted to obtain a first current abrupt change amount conversion value and a second current abrupt change amount conversion value corresponding to each sampling time.
Preferably, the performing polarity conversion on the first current abrupt change amount and the second current abrupt change amount corresponding to each sampling time to obtain a first current abrupt change amount conversion value and a second current abrupt change amount conversion value corresponding to each sampling time includes:
for any sampling time, if the first current abrupt change delta i corresponding to any sampling time t m If the polarity of (t) is negative, determining Δi M (t)＝-Δi m (t)，Δi N (t)＝-Δi n (t); otherwise, determining Δi M (t)＝Δi m (t),Δi N (t)＝Δi n (t); wherein Δi n (t) is a second current abrupt change amount corresponding to the sampling time t; Δi M (t) and Δi N (t) the first current abrupt change value and the second current abrupt change value, respectively.
In step 103, point-by-point integration is performed based on the first current abrupt change conversion value and the second current abrupt change conversion value corresponding to each sampling time, and a first integrated value and a second integrated value are obtained.
In step 104, a fault type is determined according to a preset bus protection criterion based on the first integral value and the second integral value, so as to perform bus protection based on the determined fault type.
Preferably, the determining the fault type according to a preset bus protection criterion based on the first integral value and the second integral value includes:
if the first current abrupt change conversion value and the second current abrupt change conversion value meet a preset bus protection criterion, determining that the fault type is a fault in a bus area; otherwise, determining the fault type as the fault outside the bus area;
the bus protection criterion is preset, and the bus protection criterion comprises:
Figure BDA0003852112740000131
wherein ΣΔi M (t) and ΣΔi N (t) a first current-abrupt change-over value and a second current-abrupt change-over value, respectively; i e Is the rated current of the bus bar protection device.
In the present invention, after completing the current grouping, polarity conversion is performed, including: analysis of the first processed (. DELTA.i) m (t)，Δi n (t)) if Δi m (t) polarity "-", let Δi M (t)＝-Δi m (t),Δi N (t)＝-Δi n (t) if Δi m (t) polarity "+", let Δi M (t)＝Δi m (t),Δi N (t)＝Δi n (t). After polarity conversion, sample value integration is performed, including: the data (Δi) obtained after the polarity conversion is converted M (t)，Δi N (t)) is subjected to point-by-point integration to obtain (ΣΔi) M (t)，ΣΔi N (t)). Finally, determining ΣΔi M (t) and ΣΔi N And (t) whether the bus protection criterion is met, if so, determining that the internal faults of the bus area occur, otherwise, determining that the external faults occur. The principle of judging the faults inside and outside the zone is shown in fig. 2.
Wherein, bus protection criterion is:
Figure BDA0003852112740000132
wherein ΣΔi M (t) and ΣΔi N (t) a first current-abrupt change-over value and a second current-abrupt change-over value, respectively; i e Is the rated current of the bus bar protection device.
In the invention, f when the bus bar fails in the zone 1 (t-T/4)＝0，f 1 (t)>0，Δf 1 (t)>f 1.set The criterion I is satisfied. f (f) 2 (T-T/4) is a fixed value, f 2 (t) superimposing the 100Hz component, Δf, to a fixed value 2 (t)>f 2.set Criterion II is satisfied.
F when the bus fails outside the area 1 (t-T/4)＝0，f 1 (t)＝0，Δf 1 (t)<f 1.set The criterion I is not satisfied. f (f) 2 (T-T/4) is a fixed value, f 2 (t) superimposing the 100Hz component, Δf, to a fixed value 2 (t)>f 2.set Criterion II is satisfied.
When the bus generates a fault in the external fault transfer area, the time of the external fault of the bus area is not satisfied, and the criterion I is satisfied; the fault time in the bus area is satisfied by both criteria I and II.
And for the protection element, through data preprocessing, firstly, the internal and external faults are concentrated in the 90-degree interval range of the first quadrant and the second quadrant, so that the internal and external faults can be conveniently identified. Secondly, the monotonic change characteristics of the currents in the internal and external fault processes can be intuitively displayed through the track of the integration value of the abrupt variable current sampling values in the fault period.
And when the bus area fails, the operation track of the integrated value of the current mutation sampling value rapidly enters an out-of-area failure recognition area of the second quadrant from the original point, so that the out-of-area failure is accurately recognized, and the correct and non-action is protected.
When the faults occur in the bus area, the operation track of the integrated value of the current mutation quantity sampling value rapidly enters the fault identification area in the area of the first quadrant from the original point, so that the faults in the area are accurately identified, and correct actions are protected.
When the bus generates an out-of-zone fault to an in-zone fault, after the out-of-zone fault, the operation track of the integrated value of the current mutation value rapidly enters an out-of-zone fault identification zone of the second quadrant from the original point, so that the out-of-zone fault is accurately identified, and the correct and non-action is protected. When the fault develops into an in-area fault, the operation track of the integrated value of the current mutation sampling value rapidly enters the in-area fault identification area from the out-area fault identification area, accurately identifies the in-area fault and cuts off a fault bus.
The following specifically exemplifies embodiments of the present invention
The simulation model is shown in fig. 3, wherein S1 is a wind field, the capacity is 250MW, S2-S4 are conventional alternating current systems, B1 is a bus, F1 is a fault point in a bus area, and F2 is a fault point outside the bus area.
(1) Failure in a zone
When the bus is in an A-phase grounding fault at the point F1, the primary value of fault phase current of each branch of the bus B1 is shown as in fig. 4, and it can be seen that the short-circuit current provided by the wind farm is the smallest (branch 1) and the current phases of the branches 2-4 are the same.
When the bus area fails, the action condition of the starting element criterion provided by the patent is shown in fig. 5, a red curve in the diagram is a starting element fixed value, and therefore, when the bus area fails, the criterion I, II acts 0.86ms after the failure (the first sampling value after the failure).
When the fault occurs in the busbar zone, the action track of the busbar protection criterion provided by the patent is shown in fig. 6, and after the fault, the operation track of the integrated value of the current mutation value rapidly enters the fault identification zone in the first quadrant from the origin, so that the fault in the zone is accurately identified.
(2) Out-of-zone failure
When the bus at the point F2 has an A-phase ground fault, the primary value of the fault phase current of each branch of the bus B1 is shown as in fig. 7, and it can be seen that the short-circuit current provided by the wind farm is the smallest (branch 1) during the out-of-zone fault, and the current phases of the branches 2 and 3 are opposite to the current phase of the branch 4.
When the bus area fails, the action condition of the starting element criterion provided by the patent is shown in fig. 8, a red curve in the diagram is a starting element fixed value, and it can be seen that when the bus fails outside the area, the criterion I does not act, and the criterion II acts 0.86ms after the failure (the first sampling value after the failure).
When the bus zone fails, the action track of the bus protection criterion provided by the patent is shown as a figure 9, and after the fault, the operation track of the integrated value of the current mutation value rapidly enters the out-of-zone fault identification zone of the second quadrant from the origin, so that the out-of-zone fault is accurately identified, and the correct and non-action is protected.
(3) Out-of-zone fault to in-zone fault
When the bus line F2 has an A-phase ground fault, the bus line F1 has an A-phase ground fault after the fault lasts for 80ms, the primary value of the fault phase current of each branch of the bus line B1 is shown in fig. 10, and the short circuit current provided by the wind farm is minimum (branch 1) during the fault period. When the fault occurs outside the area, the current phase of the branch circuits 2 and 3 is opposite to the current phase of the branch circuit 4; in the event of a fault in the zone, the current phases of the branches 2, 3 are identical.
When the bus fails outside the area and changes into the area, the action condition of the starting element criterion provided by the patent is shown in fig. 11, a red curve in the diagram is the fixed value of the starting element, and as can be seen, when the bus fails outside the area, the criterion I does not act, the criterion II acts after 0.86ms after the failure (the first sampling value after the failure), when the area fails, the criterion I acts after 0.86ms after the area, and meanwhile, the criterion II acts again.
When the bus generates the out-of-zone fault and changes into the in-zone fault, the action track of the bus protection criterion provided by the patent is shown as a figure 12, and the action track of the integrated value of the current mutation value rapidly enters the out-of-zone fault identification zone of the second quadrant from the origin after the out-of-zone fault, so that the out-of-zone fault is accurately identified, and the protection is correct and not actuated. When the fault develops into an in-area fault, the operation track of the integrated value of the current mutation sampling value rapidly enters the in-area fault identification area from the out-area fault identification area, accurately identifies the in-area fault and cuts off a fault bus.
Fig. 13 is a schematic structural diagram of a bus bar protection system 1300 based on current sampling values according to an embodiment of the present invention. As shown in fig. 13, a bus bar protection system 1300 based on current sampling values according to an embodiment of the present invention includes: a current abrupt amount determination unit 1301, a polarity conversion unit 1302, an integration unit 1303, and a bus bar protection unit 1304.
Preferably, the current abrupt change amount determining unit 1301 is configured to, when the starting element meets a preset starting criterion, perform current grouping based on current abrupt change amount sampling values of each side of the electrical element at each sampling time, and determine a first current abrupt change amount and a second current abrupt change amount corresponding to each sampling time.
Preferably, the current abrupt change amount determining unit 1301 performs current grouping based on current abrupt change amount sampling values of each side of the electrical element at each sampling time, and determines a first current abrupt change amount and a second current abrupt change amount corresponding to each sampling time, including:
for any sampling time, selecting the current abrupt change sampling value with the largest amplitude in the current abrupt change sampling values of all sides of the electric element at the any sampling time as the first current abrupt change delta i corresponding to the any sampling time m (t)＝max(|Δi 1 (t)|,|Δi 2 (t)|,…,|Δi h (t) |); determining a second current abrupt change delta i corresponding to any sampling moment according to the difference value between the sum of the current abrupt change sampling values of each side of the electric element at any sampling moment and the first current abrupt change n (t)＝Δi Σ (t)-Δi m (t)；
Wherein, |Δi h (t) is the amplitude of the current abrupt change sampling value at the sampling time t at the h side; Δi Σ And (t) is the sum of current mutation sampling values of each side of the electric element at the sampling time t.
Preferably, the polarity conversion unit 1302 is configured to perform polarity conversion on the first current abrupt change amount and the second current abrupt change amount corresponding to each sampling time, so as to obtain a first current abrupt change amount conversion value and a second current abrupt change amount conversion value corresponding to each sampling time.
Preferably, the polarity conversion unit 1302 performs polarity conversion on the first current abrupt change amount and the second current abrupt change amount corresponding to each sampling time to obtain a first current abrupt change amount conversion value and a second current abrupt change amount conversion value corresponding to each sampling time, which includes:
for any sampling time, if the first current abrupt change delta i corresponding to any sampling time t m If the polarity of (t) is negative, determining Δi M (t)＝-Δi m (t)，Δi N (t)＝-Δi n (t); otherwise, determining Δi M (t)＝Δi m (t),Δi N (t)＝Δi n (t); wherein Δi n (t) is a second current abrupt change amount corresponding to the sampling time t; Δi M (t) and Δi N (t) the first current abrupt change value and the second current abrupt change value, respectively.
Preferably, the integrating unit 1303 is configured to perform point-by-point integration based on the first current abrupt change conversion value and the second current abrupt change conversion value corresponding to each sampling time, and obtain a first integrated value and a second integrated value.
Preferably, the bus protection unit 1304 is configured to determine a fault type according to a preset bus protection criterion based on the first integrated value and the second integrated value, so as to perform bus protection based on the determined fault type.
Preferably, the bus protection unit 1304, based on the first integral value and the second integral value, determines a fault type according to a preset bus protection criterion, including:
If the first current abrupt change conversion value and the second current abrupt change conversion value meet a preset bus protection criterion, determining that the fault type is a fault in a bus area; otherwise, determining the fault type as the fault outside the bus area;
the bus protection criterion is preset, and the bus protection criterion comprises:
Figure BDA0003852112740000171
wherein ΣΔi M (t) and ΣΔi N (t) a first current-abrupt change-over value and a second current-abrupt change-over value, respectively; i e Is the rated current of the bus bar protection device.
Preferably, wherein the system further comprises:
the starting judgment unit is used for determining that the starting element meets the preset starting criterion when the starting element meets the first preset promoter criterion or the second preset promoter criterion;
wherein the first preset promoter criterion comprises:
Δf 1 (t)＝f 1 (t)-f 1 (t-T/4)>f 1.set ，
the second preset promoter criterion comprises:
Δf 2 (t)＝f 2 (t)-f 2 (t-T/4)>f 2.set ，
wherein Δf 1 (t) is a first variation; Δf 2 (t) is a second variation; f (f) 1 (t) is first data at time t; f (f) 2 (t) is second data at time t; f (f) 1 (T-T/4) is f before quarter cycle 1 (t)；f 2 (T-T/4) is f before quarter cycle 2 (t)；f 1.set A first preset threshold value; f (f) 2.set Is a second predetermined threshold.
Preferably, the start-up judging unit further includes:
Calculating abrupt change delta i of each branch current sampling value of bus based on each branch current sampling value of bus φ.k (t) comprising:
Δi φ.k (t)＝i φ.k (t)-i φ.k (t-T s )，
calculating differential current i of each phase based on mutation amount of current sampling value of each branch φ.Σ (t) comprising:
Figure BDA0003852112740000181
calculating first data based on each phase difference stream, comprising:
f 1 (t)＝(i a.Σ (t)-i b.Σ (t)) 2 +(i b.Σ (t)-i c.Σ (t)) 2 +(i c.Σ (t)-i a.Σ (t) 2 ，
wherein i is φ.k (t) is the current sampling value of phi phase k branch at t time, phi is A, B, C, i φ.k (t-T s ) T-T being the kth branch of phi phase s Current sample value at moment, T s The time corresponding to the one-cycle wave; n is the number of branches; f (f) 1 (t) is first data; i.e a.Σ (t)、i b.Σ (t) and i c . Σ (t) is the differential flow of three phases A, B, C at time t respectively;
dividing the current of each branch of the bus into two groups, comparing the mutation amounts of the current sampling values of all branches, and determining the branch current delta i with the maximum mutation amount of the current sampling value φ.max (t)，Δi φ.max (t)＝max(Δi φ.1 (t),Δi φ.2 (t)，…，Δi φ.k (t))；
Calculate the division Δi φ.max (t) the sum Δi of the abrupt amounts of the current sampling values of the other branches outside the corresponding branch φ.s (t)，Δi φ.s (t)＝i φ.Σ (t)-Δi φ.max (t)；
Based on i φ.s (t) calculating second data comprises:
f 2 (t)＝(i a.s (t)-i b.s (t)) 2 +(i b.s (t)-i c.s (t)) 2 +(i c.s (t)-i a.s (t)) 2 ，
wherein Δi φ.k (t) is a current sampling value of the kth branch of phi phase at the moment t, and phi is A, B, C; i.e φ.Σ (t) is the sum of the current sample values of all branches of phi phase; f (f) 2 (t) is second data.
The bus protection system 1300 based on current sampling values according to the embodiment of the present invention corresponds to the bus protection method 100 based on current sampling values according to another embodiment of the present invention, and will not be described herein.
The invention has been described with reference to a few embodiments. However, as is well known to those skilled in the art, other embodiments than the above disclosed invention are equally possible within the scope of the invention, as defined by the appended patent claims.
Generally, all terms used in the claims are to be interpreted according to their ordinary meaning in the technical field, unless explicitly defined otherwise therein. All references to "a/an/the [ means, component, etc. ]" are to be interpreted openly as referring to at least one instance of said means, component, etc., unless explicitly stated otherwise. The steps of any method disclosed herein do not have to be performed in the exact order disclosed, unless explicitly stated.
It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Finally, it should be noted that: the above embodiments are only for illustrating the technical aspects of the present invention and not for limiting the same, and although the present invention has been described in detail with reference to the above embodiments, it should be understood by those of ordinary skill in the art that: modifications and equivalents may be made to the specific embodiments of the invention without departing from the spirit and scope of the invention, which is intended to be covered by the claims.

Claims (12)

1.一种基于电流采样值的母线保护方法，其特征在于，所述方法包括：1. A busbar protection method based on current sampling value, is characterized in that, described method comprises: 当启动元件满足预设启动判据启动时，基于电气元件各侧在每一采样时刻的电流突变量采样值进行电流分组，确定每一采样时刻对应的第一电流突变量和第二电流突变量；When the starting element satisfies the preset starting criterion to start, the current is grouped based on the sampling value of the current mutation amount on each side of the electrical component at each sampling moment, and the first current mutation amount and the second current mutation amount corresponding to each sampling moment are determined. ; 分别对每一采样时刻对应的第一电流突变量和第二电流突变量进行极性转换，以获取每一采样时刻对应的第一电流突变量转换值和第二电流突变量转换值；Respectively perform polarity conversion on the first sudden change in current and the second sudden change in current corresponding to each sampling moment, so as to obtain the converted value of the first sudden change in current and the converted value of the second sudden change in current corresponding to each sampling moment; 基于每一采样时刻对应的第一电流突变量转换值和第二电流突变量转换值进行逐点积分，获取第一积分值和第二积分值；performing point-by-point integration based on the converted value of the first sudden change in current and the converted value of the second sudden change in current corresponding to each sampling moment, to obtain the first integral value and the second integral value; 基于所述第一积分值和第二积分值，按照预设母线保护判据确定故障类型，以基于确定的故障类型进行母线保护。Based on the first integral value and the second integral value, the fault type is determined according to the preset busbar protection criterion, so as to perform busbar protection based on the determined fault type. 2.根据权利要求1所述的方法，其特征在于，所述基于电气元件各侧在每一采样时刻的电流突变量采样值进行电流分组，确定每一采样时刻对应的第一电流突变量和第二电流突变量，包括：2. The method according to claim 1, characterized in that, performing current grouping based on the sampling values of current mutations at each sampling moment on each side of the electrical component, and determining the first current mutation corresponding to each sampling moment and The second sudden change in current, including: 对于任一采样时刻，选取电气元件各侧在该任一采样时刻的电流突变量采样值中幅值最大的电流突变量采样值为该任一采样时刻对应的第一电流突变量Δim(t)＝max(|Δi1(t)|,|Δi2(t)|,…,|Δih(t)|)；根据电气元件各侧在该任一采样时刻的电流突变量采样值之和与所述第一电流突变量的差值，确定该任一采样时刻对应的第二电流突变量Δin(t)＝ΔiΣ(t)-Δim(t)；For any sampling moment, select the first current mutation amount Δi m (t )＝max(|Δi 1 (t)|,|Δi 2 (t)|,...,|Δi h (t)|); According to the sum of the sampling values of the sudden change of current on each side of the electrical component at any sampling moment The difference between the first sudden change in current and the second sudden change in current Δi n (t)=Δi Σ (t)-Δi m (t) corresponding to any sampling moment is determined; 其中，|Δih(t)|为h侧在采样时刻t时的电流突变量采样值的幅值；ΔiΣ(t)为电气元件各侧在采样时刻t时的电流突变量采样值之和。Among them, |Δi h (t)| is the amplitude of the sampling value of the current mutation at the sampling time t on the h side; Δi Σ (t) is the sum of the sampling values of the current mutation at each side of the electrical component at the sampling time t . 3.根据权利要求1所述的方法，其特征在于，所述分别对每一采样时刻对应的第一电流突变量和第二电流突变量进行极性转换，以获取每一采样时刻对应的第一电流突变量转换值和第二电流突变量转换值，包括:3. The method according to claim 1, characterized in that the polarity conversion is performed on the first current sudden change and the second current sudden change corresponding to each sampling moment, so as to obtain the first current sudden change corresponding to each sampling moment A current sudden change conversion value and a second current sudden change conversion value, including: 对于任一采样时刻，若该任一采样时刻t对应的第一电流突变量Δim(t)的极性为负，则确定ΔiM(t)＝-Δim(t)，ΔiN(t)＝-Δin(t)；反之，则确定ΔiM(t)＝Δim(t),ΔiN(t)＝Δin(t)；其中，Δin(t)为采样时刻t对应的第二电流突变量；ΔiM(t)和ΔiN(t)分别为第一电流突变量转换值和第二电流突变量转换值。For any sampling moment, if the polarity of the first current mutation Δi m (t) corresponding to the arbitrary sampling moment t is negative, then determine Δi M (t)=-Δi m (t), Δi N (t )=-Δi n (t); otherwise, determine Δi M (t)=Δi m (t), Δi N (t)=Δi n (t); where Δi n (t) is the corresponding The second sudden change in current; Δi M (t) and Δi N (t) are the converted value of the first sudden change in current and the converted value of the second sudden change in current, respectively. 4.根据权利要求1所述的方法，其特征在于,所述基于所述第一积分值和第二积分值，按照预设母线保护判据确定故障类型，包括：4. The method according to claim 1, characterized in that, based on the first integral value and the second integral value, determining the fault type according to the preset busbar protection criterion comprises: 若第一电流突变量转换值和第二电流突变量转换值满足预设母线保护判据，则确定故障类型为母线区内故障；反之，则确定故障类型为母线区外故障；If the converted value of the first sudden change in current and the converted value of the second sudden change in current meet the preset busbar protection criterion, the fault type is determined to be a busbar internal fault; otherwise, the fault type is determined to be a busbar external fault; 其中，预设母线保护判据，包括：Among them, the preset busbar protection criteria include:
Figure FDA0003852112730000021
Figure FDA0003852112730000021
 其中，ΣΔiM(t)和ΣΔiN(t)分别为第一电流突变量转换值和第二电流突变量转换值；Ie为母线保护装置的额定电流。Wherein, ΣΔi M (t) and ΣΔi N (t) are the converted value of the first sudden change of current and the converted value of the second sudden change of current respectively; I e is the rated current of the bus protection device. 5.根据权利要求1所述的方法，其特征在于，所述方法还包括：5. The method according to claim 1, wherein the method further comprises: 当启动元件满足第一预设启动子判据或第二预设启动子判据时，确定所述启动元件满足预设启动判据；When the activation element satisfies the first preset promoter criterion or the second preset promoter criterion, it is determined that the activation element satisfies the preset activation criterion; 其中，所述第一预设启动子判据，包括：Wherein, the first preset promoter criterion includes: Δf1(t)＝f1(t)-f1(t-T/4)>f1.set，Δf 1 (t)=f 1 (t)-f 1 (tT/4)>f 1.set , 所述第二预设启动子判据，包括：The second preset promoter criterion includes: Δf2(t)＝f2(t)-f2(t-T/4)>f2.set，Δf 2 (t)=f 2 (t)-f 2 (tT/4)>f 2.set , 其中，Δf1(t)为第一变化量；Δf2(t)为第二变化量；f1(t)为t时刻的第一数据；f2(t)为t时刻的第二数据；f1(t-T/4)为四分之一周波之前的f1(t)；f2(t-T/4)为四分之一周波之前的f2(t)；f1.set为第一预设门槛值；f2.set为第二预设门槛值。Wherein, Δf 1 (t) is the first variation; Δf 2 (t) is the second variation; f 1 (t) is the first data at time t; f 2 (t) is the second data at time t; f 1 (tT/4) is f 1 (t) before a quarter cycle; f 2 (tT/4) is f 2 (t) before a quarter cycle; f 1.set is the first preset Set the threshold value; f 2.set is the second preset threshold value. 6.根据权利要求5所述的方法，其特征在于，所述方法还包括：6. The method according to claim 5, further comprising: 基于母线各支路电流采样值计算母线各支路电流采样值的突变量Δiφ.k(t)，包括：Calculate the sudden change Δi φ.k (t) of the current sampling value of each branch of the bus based on the current sampling value of each branch of the bus, including: Δiφ.k(t)＝iφ.k(t)-iφ.k(t-Ts)，Δi φ.k (t)=i φ.k (t)-i φ.k (tT s ), 基于各支路电流采样值的突变量计算各相的差流iφ.Σ(t)，包括：Calculate the differential current i φ.Σ (t) of each phase based on the sudden change of the current sampling value of each branch, including:
Figure FDA0003852112730000031
Figure FDA0003852112730000031
 基于各相差流计算第一数据，包括：calculating the first data based on each phase difference flow, including: f1(t)＝(ia.Σ(t)-ib.Σ(t))2+(ib.Σ(t)-ic.Σ(t))2+(ic.Σ(t)-ia.Σ(t)2，f 1 (t)＝(i a.Σ (t)-i b.Σ (t)) 2 +(i b.Σ (t)-i c.Σ (t)) 2 +(i c.Σ ( t)-i a.Σ (t) 2 , 其中，iφ.k(t)为φ相第k条支路t时刻的电流采样值，φ取值为A、B、C，iφ.k(t-Ts)为φ相第k条支路t-Ts时刻的电流采样值，Ts为一周波对应的时间；N为支路的数量；f1(t)为第一数据；ia.Σ(t)、ib.Σ(t)和ic.Σ(t)分别为t时刻A、B、C三相的差流；Among them, i φ.k (t) is the current sampling value of the kth branch of φ phase at time t, and the values of φ are A, B, C, and i φ.k (tT s ) is the kth branch of φ phase Current sampling value at time tT s , T s is the time corresponding to one cycle; N is the number of branches; f 1 (t) is the first data; i a.Σ (t), i b.Σ (t) and i c.Σ (t) are the differential currents of A, B, and C phases at time t; 将母线各支路电流分为两组，比较所有支路的电流采样值的突变量，确定电流采样值的突变量最大的支路电流Δiφ.max(t)，Δiφ.max(t)＝max(Δiφ.1(t),Δiφ.2(t)，…，Δiφ.k(t))；Divide the current of each branch of the bus into two groups, compare the sudden changes of the current sampling values of all branches, and determine the branch current with the largest sudden change of the current sampling values Δi φ.max (t), Δi φ.max (t) = max(Δi φ.1 (t), Δi φ.2 (t), ..., Δi φ.k (t)); 计算除Δiφ.max(t)对应支路外其余支路的电流采样值的突变量之和Δiφ.s(t)，Δiφ.s(t)＝iφ.Σ(t)-Δiφ.max(t)；Calculate the sum of the sudden changes of the current sampling values of the other branches except the branch corresponding to Δi φ.max (t) , Δi φ.s ( t), Δi φ.s (t)=i φ.Σ (t)-Δi φ.max (t); 基于iφ.s(t)计算第二数据包括：Calculating the second data based on i φ.s (t) includes: f2(t)＝(ia.s(t)-ib.s(t))2+(ib.s(t)-ic.s(t))2+(ic.s(t)-ia.s(t))2，f 2 (t)＝(i as (t)-i bs (t)) 2 +(i bs (t)-i cs (t)) 2 +(i cs (t)-i as (t)) 2 , 其中，Δiφ.k(t)为φ相第k条支路t时刻的电流采样值，φ取值为A、B、C；iφ.Σ(t)为φ相所有支路的电流采样值之和；f2(t)为第二数据。Among them, Δi φ.k (t) is the current sampling value of the kth branch of φ phase at time t, and the values of φ are A, B, C; i φ.Σ (t) is the current sampling value of all branches of φ phase sum of values; f 2 (t) is the second data. 7.一种基于电流采样值的母线保护系统，其特征在于，所述系统包括：7. A bus protection system based on current sampling values, characterized in that the system includes: 电流突变量确定单元，用于当启动元件满足预设启动判据启动时，基于电气元件各侧在每一采样时刻的电流突变量采样值进行电流分组，确定每一采样时刻对应的第一电流突变量和第二电流突变量；The current sudden change determination unit is used to perform current grouping based on the sampling values of the current sudden change at each sampling moment on each side of the electrical component when the starting element meets the preset starting criterion to start, and determine the first current corresponding to each sampling moment mutation amount and the second current mutation amount; 极性转换单元，用于分别对每一采样时刻对应的第一电流突变量和第二电流突变量进行极性转换，以获取每一采样时刻对应的第一电流突变量转换值和第二电流突变量转换值；A polarity conversion unit, configured to perform polarity conversion on the first current sudden change and the second current sudden change corresponding to each sampling moment, so as to obtain the first current sudden change conversion value and the second current corresponding to each sampling moment mutation conversion value; 积分单元，用于基于每一采样时刻对应的第一电流突变量转换值和第二电流突变量转换值进行逐点积分，获取第一积分值和第二积分值；The integration unit is used to perform point-by-point integration based on the first current mutation value conversion value and the second current mutation value conversion value corresponding to each sampling moment, and obtain the first integration value and the second integration value; 母线保护单元，用于基于所述第一积分值和第二积分值，按照预设母线保护判据确定故障类型，以基于确定的故障类型进行母线保护。The busbar protection unit is configured to determine a fault type according to a preset busbar protection criterion based on the first integral value and the second integral value, so as to perform busbar protection based on the determined fault type. 8.根据权利要求7所述的系统，其特征在于，所述电流突变量确定单元，基于电气元件各侧在每一采样时刻的电流突变量采样值进行电流分组，确定每一采样时刻对应的第一电流突变量和第二电流突变量，包括：8. The system according to claim 7, characterized in that the current mutation amount determining unit performs current grouping based on the current mutation amount sampling values of each side of the electrical component at each sampling moment, and determines the current mutation amount corresponding to each sampling moment The first sudden change in current and the second sudden change in current include: 对于任一采样时刻，选取电气元件各侧在该任一采样时刻的电流突变量采样值中幅值最大的电流突变量采样值为该任一采样时刻对应的第一电流突变量Δim(t)＝max(|Δi1(t)|,|Δi2(t)|,…,|Δih(t)|)；根据电气元件各侧在该任一采样时刻的电流突变量采样值之和与所述第一电流突变量的差值，确定该任一采样时刻对应的第二电流突变量Δin(t)＝ΔiΣ(t)-Δim(t)；For any sampling moment, select the first current mutation amount Δi m (t )＝max(|Δi 1 (t)|,|Δi 2 (t)|,...,|Δi h (t)|); According to the sum of the sampling values of the sudden change of current on each side of the electrical component at any sampling moment The difference between the first sudden change in current and the second sudden change in current Δi n (t)=Δi Σ (t)-Δi m (t) corresponding to any sampling moment is determined; 其中，|Δih(t)|为h侧在采样时刻t时的电流突变量采样值的幅值；ΔiΣ(t)为电气元件各侧在采样时刻t时的电流突变量采样值之和。Among them, |Δi h (t)| is the amplitude of the sampling value of the current mutation at the sampling time t on the h side; Δi Σ (t) is the sum of the sampling values of the current mutation at each side of the electrical component at the sampling time t . 9.根据权利要求7所述的系统，其特征在于，所述极性转换单元，分别对每一采样时刻对应的第一电流突变量和第二电流突变量进行极性转换，以获取每一采样时刻对应的第一电流突变量转换值和第二电流突变量转换值，包括:9. The system according to claim 7, wherein the polarity conversion unit performs polarity conversion on the first current sudden change and the second current sudden change corresponding to each sampling moment, so as to obtain each The first current sudden change value conversion value and the second current sudden change value conversion value corresponding to the sampling moment include: 对于任一采样时刻，若该任一采样时刻t对应的第一电流突变量Δim(t)的极性为负，则确定ΔiM(t)＝-Δim(t)，ΔiN(t)＝-Δin(t)；反之，则确定ΔiM(t)＝Δim(t),ΔiN(t)＝Δin(t)；其中，Δin(t)为采样时刻t对应的第二电流突变量；ΔiM(t)和ΔiN(t)分别为第一电流突变量转换值和第二电流突变量转换值。For any sampling moment, if the polarity of the first current mutation Δi m (t) corresponding to the arbitrary sampling moment t is negative, then determine Δi M (t)=-Δi m (t), Δi N (t )=-Δi n (t); otherwise, determine Δi M (t)=Δi m (t), Δi N (t)=Δi n (t); where Δi n (t) is the corresponding The second sudden change in current; Δi M (t) and Δi N (t) are the converted value of the first sudden change in current and the converted value of the second sudden change in current, respectively. 10.根据权利要求7所述的系统，其特征在于,所述母线保护单元，基于所述第一积分值和第二积分值，按照预设母线保护判据确定故障类型，包括：10. The system according to claim 7, wherein the busbar protection unit, based on the first integral value and the second integral value, determines the fault type according to a preset busbar protection criterion, including: 若第一电流突变量转换值和第二电流突变量转换值满足预设母线保护判据，则确定故障类型为母线区内故障；反之，则确定故障类型为母线区外故障；If the converted value of the first sudden change in current and the converted value of the second sudden change in current meet the preset busbar protection criterion, the fault type is determined to be a busbar internal fault; otherwise, the fault type is determined to be a busbar external fault; 其中，预设母线保护判据，包括：Among them, the preset busbar protection criteria include:
Figure FDA0003852112730000051
Figure FDA0003852112730000051
 其中，ΣΔiM(t)和ΣΔiN(t)分别为第一电流突变量转换值和第二电流突变量转换值；Ie为母线保护装置的额定电流。Wherein, ΣΔi M (t) and ΣΔi N (t) are the converted value of the first sudden change of current and the converted value of the second sudden change of current respectively; I e is the rated current of the bus protection device. 11.根据权利要求7所述的系统，其特征在于，所述系统还包括：11. The system according to claim 7, further comprising: 启动判断单元，用于当启动元件满足第一预设启动子判据或第二预设启动子判据时，确定所述启动元件满足预设启动判据；An activation judging unit, configured to determine that the activation element satisfies a preset activation criterion when the activation element satisfies a first preset promoter criterion or a second preset promoter criterion; 其中，所述第一预设启动子判据，包括：Wherein, the first preset promoter criterion includes: Δf1(t)＝f1(t)-f1(t-T/4)>f1.set，Δf 1 (t)=f 1 (t)-f 1 (tT/4)>f 1.set , 所述第二预设启动子判据，包括：The second preset promoter criterion includes: Δf2(t)＝f2(t)-f2(t-T/4)>f2.set，Δf 2 (t)=f 2 (t)-f 2 (tT/4)>f 2.set , 其中，Δf1(t)为第一变化量；Δf2(t)为第二变化量；f1(t)为t时刻的第一数据；f2(t)为t时刻的第二数据；f1(t-T/4)为四分之一周波之前的f1(t)；f2(t-T/4)为四分之一周波之前的f2(t)；f1.set为第一预设门槛值；f2.set为第二预设门槛值。Wherein, Δf 1 (t) is the first variation; Δf 2 (t) is the second variation; f 1 (t) is the first data at time t; f 2 (t) is the second data at time t; f 1 (tT/4) is f 1 (t) before a quarter cycle; f 2 (tT/4) is f 2 (t) before a quarter cycle; f 1.set is the first preset Set the threshold value; f 2.set is the second preset threshold value. 12.根据权利要求11所述的系统，其特征在于，所述启动判断单元，还包括：12. The system according to claim 11, wherein the startup judging unit further comprises: 基于母线各支路电流采样值计算母线各支路电流采样值的突变量Δiφ.k(t)，包括：Calculate the sudden change Δi φ . k (t) of the current sampling value of each branch of the bus based on the current sampling value of each branch of the bus, including: Δiφ.k(t)＝iφ.k(t)-iφ.k(t-Ts)，Δi φ.k (t)=i φ.k (t)-i φ.k (tT s ), 基于各支路电流采样值的突变量计算各相的差流iφ.Σ(t)，包括：Calculate the differential current i φ.Σ (t) of each phase based on the sudden change of the current sampling value of each branch, including:
Figure FDA0003852112730000052
Figure FDA0003852112730000052
 基于各相差流计算第一数据，包括：calculating the first data based on each phase difference flow, including: f1(t)＝(ia.Σ(t)-ib.Σ(t))2+(ib.Σ(t)-ic.Σ(t))2+(ic.Σ(t)-ia.Σ(t)2，f 1 (t)＝(i a.Σ (t)-i b.Σ (t)) 2 +(i b.Σ (t)-i c.Σ (t)) 2 +(i c.Σ ( t)-i a.Σ (t) 2 , 其中，iφ.k(t)为φ相第k条支路t时刻的电流采样值，φ取值为A、B、C，iφ.k(t-Ts)为φ相第k条支路t-Ts时刻的电流采样值，Ts为一周波对应的时间；N为支路的数量；f1(t)为第一数据；ia.Σ(t)、ib.Σ(t)和ic.Σ(t)分别为t时刻A、B、C三相的差流；Among them, i φ.k (t) is the current sampling value of the kth branch of φ phase at time t, and the values of φ are A, B, C, and i φ.k (tT s ) is the kth branch of φ phase Current sampling value at time tT s , T s is the time corresponding to one cycle; N is the number of branches; f 1 (t) is the first data; i a.Σ (t), i b.Σ (t) and i c . Σ (t) are the differential currents of A, B, and C phases at time t; 将母线各支路电流分为两组，比较所有支路的电流采样值的突变量，确定电流采样值的突变量最大的支路电流Δiφ.max(t)，Δiφ.max(t)＝max(Δiφ.1(t),Δiφ.2(t)，…，Δiφ.k(t))；Divide the current of each branch of the bus into two groups, compare the sudden changes of the current sampling values of all branches, and determine the branch current with the largest sudden change of the current sampling values Δi φ.max (t), Δi φ.max (t) = max(Δi φ . 1 (t), Δi φ.2 (t), ..., Δi φ.k (t)); 计算除Δiφ.max(t)对应支路外其余支路的电流采样值的突变量之和Δiφ.s(t)，Δiφ.s(t)＝iφ.Σ(t)-Δiφ.max(t)；Calculate the sum of the sudden changes of the current sampling values of the other branches except the branch corresponding to Δi φ.max (t) , Δi φ.s ( t), Δi φ.s (t)=i φ.Σ (t)-Δi φ.max (t); 基于iφ.s(t)计算第二数据包括：Calculating the second data based on i φ.s (t) includes: f2(t)＝(ia.s(t)-ib.s(t))2+(ib.s(t)-ic.s(t))2+(ic.s(t)-ia.s(t))2，f 2 (t)＝(i as (t)-i bs (t)) 2 +(i bs (t)-i cs (t)) 2 +(i cs (t)-i as (t)) 2 , 其中，Δiφ.k(t)为φ相第k条支路t时刻的电流采样值，φ取值为A、B、C；iφ.Σ(t)为φ相所有支路的电流采样值之和；f2(t)为第二数据。Among them, Δi φ.k (t) is the current sampling value of the kth branch of φ phase at time t, and the values of φ are A, B, C; i φ.Σ (t) is the current sampling value of all branches of φ phase sum of values; f 2 (t) is the second data.
CN202211137833.3A 2022-09-19 2022-09-19 Bus protection method and system based on current sampling value Pending CN116169646A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211137833.3A CN116169646A (en) 2022-09-19 2022-09-19 Bus protection method and system based on current sampling value

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211137833.3A CN116169646A (en) 2022-09-19 2022-09-19 Bus protection method and system based on current sampling value

Publications (1)

Publication Number Publication Date
CN116169646A true CN116169646A (en) 2023-05-26

Family

ID=86417039

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211137833.3A Pending CN116169646A (en) 2022-09-19 2022-09-19 Bus protection method and system based on current sampling value

Country Status (1)

Country Link
CN (1) CN116169646A (en)

Similar Documents

Publication Publication Date Title
CN109444640B (en) Power distribution network single-phase high-resistance earth fault detection method, system and storage medium
Hu et al. High-frequency voltage injection based stator interturn fault detection in permanent magnet machines
Elbouchikhi et al. Generalized likelihood ratio test based approach for stator-fault detection in a PWM inverter-fed induction motor drive
CN109507516A (en) Earth-fault detecting method, system and storage medium based on steady state fault amount
CN110907761B (en) Continuous line selection method and system for single-phase earth fault
CN110994553B (en) A fixed value adjustment method for transformer differential protection based on waveform decay characteristics
US10551425B2 (en) Method for quickly identifying disconnection of CT in protection of 3/2 connection mode based bus
Pietrzak et al. Stator Winding Fault Detection of Permanent Magnet Synchronous Motors Based on the Short-Time Fourier Transform
CN106054016A (en) Method for judging short circuit fault of inverter
Raj et al. A novel open-transistor fault detection method in symmetric cascaded H-bridge multilevel inverter
CN116169646A (en) Bus protection method and system based on current sampling value
CN110568300A (en) A single-phase-to-ground fault identification method for distribution network based on multi-source information
CN112736845B (en) CT trailing current identification method and device based on current phase angle difference calculation and failure protection method
CN107658852A (en) A kind of guard method of AC excitation synchronous machinery excitation winding phase fault and system
CN109245049B (en) Over-current protection method and device for rotor winding of alternating-current excitation synchronous motor
CN112180206A (en) Power distribution network single-phase earth fault positioning method and device
CN113437732B (en) Pilot protection method and system for photovoltaic power generation grid-connected tie line
CN113708343B (en) Stator grounding protection method and device for salient pole generator and electronic equipment
CN115144700A (en) New energy grid-connected system fault protection method based on injection signal amplitude
Liu et al. Primary admittance based fault detection for inter-turn short circuit in 3-phase power transformers
Gomathy et al. IPSO based fault analysis in power transformer
Li et al. Fault Detection in Ring-bus Multi-terminal DC System Based on Modified S Transform
CN111796203A (en) Identification method and device for single-phase earth fault at network bridge side of static frequency converter
CN117554750B (en) A method and system for determining the direction of phase-to-phase fault in a distribution network
CN112670951B (en) Method and device for judging action of apparent power proportional differential element of phase-splitting transformer

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