CN113690859A - Starting and acting method and system for power grid differential protection - Google Patents

Starting and acting method and system for power grid differential protection Download PDF

Info

Publication number
CN113690859A
CN113690859A CN202110892722.2A CN202110892722A CN113690859A CN 113690859 A CN113690859 A CN 113690859A CN 202110892722 A CN202110892722 A CN 202110892722A CN 113690859 A CN113690859 A CN 113690859A
Authority
CN
China
Prior art keywords
value
dtw
differential protection
starting
current
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110892722.2A
Other languages
Chinese (zh)
Other versions
CN113690859B (en
Inventor
吴通华
濮宏飞
姚刚
洪丰
郑小江
王威威
李新东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nari Technology Co Ltd
Original Assignee
Nari Technology Co Ltd
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 Nari Technology Co Ltd filed Critical Nari Technology Co Ltd
Priority to CN202110892722.2A priority Critical patent/CN113690859B/en
Publication of CN113690859A publication Critical patent/CN113690859A/en
Application granted granted Critical
Publication of CN113690859B publication Critical patent/CN113690859B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/26—Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured

Landscapes

  • Emergency Protection Circuit Devices (AREA)

Abstract

The invention discloses a starting and acting method and a system of power grid differential protection, wherein the method comprises the steps of carrying out DTW calculation on a normalized current sampling sequence and carrying out unitization processing on an obtained DTW value; aiming at the DTW value after the unitization processing, calculating a starting value of a starting method and a setting value of an action method by combining a short window algorithm; judging whether a fault occurs or not and whether the fault is in a line area or not based on the calculated starting value and setting value; if the judgment result is that the power grid line has a fault, starting differential protection; if the judgment result is that the circuit area is in fault, the differential protection is controlled to immediately output an action signal. The invention can reliably start differential protection when different types of faults occur at different positions of a power grid line, and the system can greatly reduce the oscillation amplitude of the DTW distance waveform of the action method in an alternating current system and enable the differential protection to rapidly output an action signal.

Description

Starting and acting method and system for power grid differential protection
Technical Field
The invention belongs to the technical field of differential protection, and particularly relates to a starting and acting method and system of power grid differential protection.
Background
With the massive access of distributed new energy, the traditional staged current protection of the power distribution network line is not suitable any more. In order to effectively identify faults inside and outside the area and improve the power supply reliability of a power grid, a differential protection principle needs to be introduced and a communication channel capable of meeting the requirement of protecting data transmission is supplemented.
At present, both optical fiber and 5G communication technologies can be used for differential protection in data transmission performance. Compared with optical fiber communication, the base station in the 5G network is low in laying difficulty and flexible in site selection, so that the 5G communication is selected as a protection channel, the differential protection is promoted, and the reliability of power grid protection under the condition of high-proportion new energy access is improved. However, to use 5G communication as a data channel for differential protection, the problems of data synchronization and packet loss still need to be solved.
Differential protection based on Dynamic Time Warping (DTW) algorithm requires calculating the DTW distance between two sets of time sequences. When the problem of sampling point synchronization error increase or data loss occurs, the influence on the DTW distance is small, and the DTW algorithm weakens the adverse effect brought by 5G communication in principle.
Journal article "a differential protection algorithm for resisting communication delay and packet loss" provides a differential protection action method based on a DTW algorithm by using a DTW distance between phase currents on two sides. Journal article "positioning method for small current grounding fault section based on dynamic time bending distance" plans a differential protection action method through the DTW distance between the zero sequence currents on both sides, and can effectively detect the small current grounding fault. The chinese patent application No. 202010369821.8 discloses a method for calculating a DTW value using a normalized current sequence, which can effectively identify a single-phase ground fault. In the method in the prior art, the length of the data window is greater than or equal to one power frequency period, and in an alternating current system, if the length of the data window is reduced, extreme values occur in the DTW distance near the zero crossing point and the peak valley value of the current, and the corresponding waveform greatly oscillates. In addition, the protection starting method based on the DTW algorithm is not researched in the methods.
Disclosure of Invention
In order to solve the problems, the invention provides a starting and acting method and a system of power grid differential protection, which can realize the function of reliably starting the differential protection when different types of faults occur at different positions of a power grid line, and the system can greatly reduce the oscillation amplitude of a DTW distance waveform of an acting method in an alternating current system and enable the differential protection to rapidly output an acting signal.
In order to achieve the technical purpose and achieve the technical effects, the invention is realized by the following technical scheme:
in a first aspect, the present invention provides a method for starting and operating a power grid differential protection, including:
performing DTW calculation on the normalized current sampling sequence, and performing unitization processing on the obtained DTW value;
aiming at the DTW value after the unitization processing, calculating a starting value of a starting method and a setting value of an action method by combining a short window algorithm;
judging whether a fault occurs or not and whether the fault is in a line area or not based on the calculated starting value and setting value;
if the judgment result is that the power grid line has a fault, starting differential protection;
if the judgment result is that the circuit area is in fault, the differential protection is controlled to immediately output an action signal.
Optionally, the normalized current sampling sequence is obtained by:
after receiving a phase current sampling sequence of an opposite end, processing the phase current sampling sequence of the opposite end by using a substitution processing method, wherein the substitution processing method is to use previous effective sampling point data to substitute current abnormal sampling point data;
and selecting two current sampling sequences from all the phase current sampling sequences and the zero sequence current sampling sequences at two ends of the line, and normalizing the two current sampling sequences.
Optionally, the DTW calculation specifically includes:
intercepting the last T of the two normalized current sampling sequences0Sampling point data;
and calculating the DTW value between the sampling points of the intercepted parts of the two current sampling sequences by using a DTW algorithm.
Optionally, the calculation formula adopted by the alternative processing method is as follows:
Ii'＝Ii-1,i＝2,3,4...
wherein, Ii-1Is the amplitude of the current at the I-1 th effective sampling pointi' is the magnitude of the current at the ith sample point after the substitution process.
Optionally, the calculation formula of the DTW value is:
Dij＝dij+min{Di,j-1,Di-1,j,Di-1,j-1}
wherein d isijRepresenting the Euclidean distance between the current amplitude of the ith sampling point of the first current sampling sequence and the current amplitude of the jth sampling point of the second current sampling sequence; dijRepresenting the ith-T of a sequence of samples taken by the first current0+1，i-T0+2, …, i sample point current amplitudes and j-T of the second sequence of current samples0+1，j-T0The DTW value between +2, …, j sample point current amplitudes.
Optionally, the formula used for unitizing the obtained DTW value is:
Figure BDA0003196583120000021
d is a DTW actual accumulated value between the two current sampling sequences;
Figure BDA0003196583120000022
is the DTW value after the unit.
Optionally, the calculation formula of the short window algorithm is:
Figure BDA0003196583120000023
Figure BDA0003196583120000031
wherein the content of the first and second substances,
Figure BDA0003196583120000032
for the first time after being optimized by a short window algorithm
Figure BDA0003196583120000033
The DTW difference of the phase is the DTW difference at two sampling points separated by 1/4 power frequency cycles
Figure BDA0003196583120000034
Average value of (d);
Figure BDA0003196583120000035
for the first time after being optimized by a short window algorithm
Figure BDA0003196583120000036
DTW braking amount of phase, which is DTW braking amount
Figure BDA0003196583120000037
Average value of (d); t is the total number of sampling points in a power frequency period, and k is a sampling serial number.
Optionally, the promoter method comprises promoter method 1, promoter method 2, promoter method 3;
the promoter value of the promoter method 1 was calculated by the following steps:
defining a current sampling sequence to be divided into three groups of AB, BC and CA;
for group AB, then:
Figure BDA0003196583120000038
Figure BDA0003196583120000039
Figure BDA00031965831200000310
in the formula (I), the compound is shown in the specification,
Figure BDA00031965831200000311
for the optimized DTW difference between the two phases AB,
Figure BDA00031965831200000312
for the optimized DTW braking amount of the two AB phases through the short window algorithm,
Figure BDA00031965831200000313
and
Figure BDA00031965831200000314
are respectively as
Figure BDA00031965831200000315
And
Figure BDA00031965831200000316
the amount of change in (c); t is the total number of sampling points in a power frequency period, and k is a sampling serial number;
if the starting value is
Figure BDA00031965831200000317
Greater than a threshold value Dop1If yes, differential protection is started;
for the BC group, then:
Figure BDA00031965831200000318
Figure BDA00031965831200000319
Figure BDA00031965831200000320
in the formula (I), the compound is shown in the specification,
Figure BDA00031965831200000321
for the optimized BC two-phase DTW difference momentum through the short window algorithm,
Figure BDA00031965831200000322
for the optimized BC two-phase DTW braking amount through the short window algorithm,
Figure BDA00031965831200000323
and
Figure BDA00031965831200000324
are respectively as
Figure BDA00031965831200000326
And
Figure BDA00031965831200000325
the amount of change in (c); t is the total number of sampling points in a power frequency period, and k is a sampling serial number;
if the starting value is
Figure BDA00031965831200000327
Greater than a threshold value Dop1If yes, differential protection is started;
for the CA group, then:
Figure BDA00031965831200000328
Figure BDA00031965831200000329
Figure BDA0003196583120000041
in the formula (I), the compound is shown in the specification,
Figure BDA0003196583120000042
for the DTW difference between two phases of CA after the optimization of the short window algorithm,
Figure BDA0003196583120000043
for the optimized DTW braking amount between two phases of CA by the short window algorithm,
Figure BDA0003196583120000044
and
Figure BDA0003196583120000045
are respectively as
Figure BDA0003196583120000047
And
Figure BDA0003196583120000046
the amount of change in (c); t is the total number of sampling points in a power frequency period, and k is a sampling serial number;
if the starting value is
Figure BDA0003196583120000048
Greater than a threshold value Dop1If yes, differential protection is started;
promoter method 2 promoter values were calculated by the following steps:
for phase A current, split it into odd sequences IA1And even sequence IA2；
IA1(k)＝IA(2k-1),k＝1,2...
-IA2(k)＝-IA(2k),k＝1,2...
Calculation of IA1and-IA2DTW value of
Figure BDA0003196583120000049
Calculation of IA1And IA2DTW value of
Figure BDA00031965831200000410
Calculated by short window algorithm
Figure BDA00031965831200000411
And
Figure BDA00031965831200000412
then there are:
Figure BDA00031965831200000413
if the starting value is
Figure BDA00031965831200000414
Greater than a threshold value Dop2If yes, differential protection is started;
for the phase B current, split it into odd sequences IB1And even sequence IB2；
IB1(k)＝IB(2k-1),k＝1,2...
-IB2(k)＝-IB(2k),k＝1,2...
Calculation of IB1and-IB2DTW value of
Figure BDA00031965831200000415
Calculation of IB1And IB2DTW value of
Figure BDA00031965831200000416
Calculated by short window algorithm
Figure BDA00031965831200000417
And
Figure BDA00031965831200000418
then there are:
Figure BDA00031965831200000419
if the starting value is
Figure BDA00031965831200000420
Greater than a threshold value Dop2If yes, differential protection is started;
for the C-phase current, split it into odd sequences IC1And even sequence IC2；
IC1(k)＝IC(2k-1),k＝1,2...
-IC2(k)＝-IC(2k),k＝1,2...
Calculation of IC1and-IC2DTW value of
Figure BDA00031965831200000421
Calculation of IC1And IC2DTW value of
Figure BDA00031965831200000422
Calculated by short window algorithm
Figure BDA00031965831200000423
And
Figure BDA00031965831200000424
then there are:
Figure BDA00031965831200000425
if the starting value is
Figure BDA00031965831200000426
Greater than a threshold value Dop2If yes, differential protection is started;
the promoter value for promoter method 3 was calculated by the following steps:
adding the three-phase current sampling sequences and dividing by 3 to obtain zero-sequence current Im0Is shown bym0Divided into odd sequences Im01And even sequence Im02Calculating Im01and-Im02DTW value of
Figure BDA0003196583120000051
Calculation of Im01And Im02DTW value of
Figure BDA0003196583120000052
Recalculation
Figure BDA0003196583120000053
And
Figure BDA0003196583120000054
then there are:
Figure BDA0003196583120000055
if the starting value is
Figure BDA0003196583120000056
Greater than a threshold value Dop3Then the differential protection is started, m represents one of the two ends of the line, and 0 represents zero sequence.
Optionally, if the determination result is that the power grid line fails, the differential protection is started, specifically:
and if the starting value of any one of the method 1, 2 and 3 is greater than the threshold value, judging that the power grid line has a fault, and starting differential protection.
Optionally, the action method comprises action submethod 1 and action submethod 2;
the setting value of the action sub-method 1 is calculated by the following steps:
obtaining the DTW differential quantity optimized by the short window algorithm by using two current sampling sequences of the same phase at two ends of the line
Figure BDA0003196583120000057
And amount of braking
Figure BDA0003196583120000058
The time interval of two sampling points in the short window algorithm is 1/8 power frequency periods, and then:
Figure BDA0003196583120000059
if the setting value
Figure BDA00031965831200000510
Greater than a threshold value Dop4If so, the differential protection action is carried out, K is a braking coefficient and belongs to [0,1]]；
The setting value of the action sub-method 2 is calculated by the following steps:
calculating zero sequence current at two ends of the line, dividing the zero sequence current into odd-even sequences, and obtaining the DTW differential quantity optimized by the short window algorithm
Figure BDA00031965831200000511
And amount of braking
Figure BDA00031965831200000512
The time interval of two sampling points in the short window algorithm is 1/8 power frequency periods, and then:
Figure BDA00031965831200000513
if the setting value
Figure BDA00031965831200000514
Greater than a threshold value Dop5Then, differential protection is performed.
Optionally, if the determination result is an internal fault of the line area, controlling the differential protection to immediately output the action signal, specifically:
if the setting value of any one of the sub-methods 1 and 2 is larger than the threshold value, the occurrence of the fault in the line area is judged, and the differential protection immediately acts.
In a third aspect, the present invention provides a system for starting and operating a power grid differential protection, including a storage medium and a processor;
the storage medium is used for storing instructions;
the processor is configured to operate in accordance with the instructions to perform the method according to any one of the first aspects.
Compared with the prior art, the invention has the beneficial effects that:
the invention provides a starting and acting method and a system of power grid differential protection, which can realize the function of reliably starting the differential protection when different types of faults occur at different positions of a power grid line, and the system can greatly reduce the oscillation amplitude of a DTW distance waveform of an acting method in an alternating current system and enable the differential protection to rapidly output an acting signal.
The short window algorithm can effectively solve the problem that the waveform of the DTW algorithm in the alternating current system greatly vibrates during short window calculation, and realizes differential protection setting based on the DTW algorithm, so that the required data storage capacity is small.
The substitution processing method can greatly reduce the influence of the DTW algorithm on the problem sampling points. When continuous multipoint data loss occurs or frames are lost due to check errors, differential protection can be prevented from malfunction by increasing the threshold value, and reliability is improved.
Drawings
In order that the present disclosure may be more readily and clearly understood, reference is now made to the following detailed description of the present disclosure taken in conjunction with the accompanying drawings, in which:
FIG. 1 is a flow diagram of a differential protection system comprising submethods 1 and 3 of a startup method and an action method;
FIG. 2 is a flow chart of the differential protection system consisting of submethods 2 and 3 of the startup method and the action method;
FIG. 3 is a schematic diagram of a power distribution network model in which an embodiment is implemented;
FIG. 4(a) shows a sub-method 1 of the action method without introducing the short window algorithm at f2Setting value waveform curve when ABCG fault occurs;
FIG. 4(b) shows a sub-method 1 of the action method without introducing the short window algorithm at f2A normalized current waveform curve when an ABCG fault occurs;
FIG. 5 is a diagram of a startup methodSub-method 1
Figure BDA0003196583120000061
And
Figure BDA0003196583120000062
at f4A waveform curve when an AG fault occurs;
FIG. 6 shows the startup value at f in submethod 2 of the startup method3A waveform curve when an AB fault occurs;
FIG. 7 is a submethod 3 of the startup method at f2A waveform curve when an AG1000 omega fault occurs;
FIG. 8 shows that sub-method 1 of the operation method in the normal 5G communication has a setting value f under the condition of introducing and not introducing the short window algorithm4Waveform when an AB fault occurs;
FIG. 9 shows that the setting value in FIG. 8 is f when continuous 3-point data packet loss occurs2A waveform at which an ABG fault occurs;
FIG. 10 shows that the setting value in FIG. 8 is f when continuous 12-point data packet loss occurs6Waveform when an AB fault occurs;
FIG. 11 shows a sub-method 2 of the method after the short window algorithm is introduced when the 5G communication is normal and continuous 3-point data packet loss occurs, at f1ABG failure.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail with reference to the following embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and do not limit the scope of the invention.
The following detailed description of the principles of the invention is provided in connection with the accompanying drawings.
Example 1
The embodiment of the invention provides a starting and acting method of power grid differential protection, which specifically comprises the following steps:
performing DTW calculation on the normalized current sampling sequence, and performing unitization processing on the obtained DTW value;
aiming at the DTW value after the unitization processing, calculating a starting value of a starting method and a setting value of an action method by combining a short window algorithm;
judging whether a fault occurs or not and whether the fault is in a line area or not based on the calculated starting value and setting value;
if the judgment result is that the power grid line has a fault, starting differential protection;
if the judgment result is that the circuit area is in fault, the differential protection is controlled to immediately output an action signal.
Referring to fig. 1-2, the detailed description of the startup and operation method in the embodiment of the present invention is provided in conjunction with the specific implementation process and the power distribution network model in fig. 3.
Step 1: sampling three-phase current values at two ends of a line in real time at a certain sampling frequency, and sending an obtained sampling sequence to a sampling selection device at one end by using a 5G network;
step 2: after receiving the current sampling sequence of the opposite end, carrying out normalization processing on each current sequence, wherein the effect is shown in fig. 4(b), the current waveform amplitude values are all in the interval of [0,1], the substitution processing method of the current abnormal sampling point is replaced by the previous effective sampling point data, the data packet loss problem is processed, and the lengths of 2 sampling sequences in the subsequent DTW calculation are ensured to be equal; the calculation formula of the alternative processing method is as follows:
I′i＝Ii-1,i＝2,3,4...
wherein, Ii-1Is the amplitude of the current at the I-1 th effective sampling pointi' is the magnitude of the current at the ith sample point after the substitution process.
And step 3: respectively intercepting the last T of 2 sampling sequences required by the subsequent DTW calculation according to the requirements of the differential protection starting and action method0Calculating corresponding DTW values by using short window algorithm, as shown in FIG. 6, FIG. 8, FIG. 9, FIG. 10 and FIG. 11, and then performing unitization processing on the obtained DTW values, and listing the starting values of the starting method and the setting values of the action method; the DTW calculation specifically comprises the steps of firstly calculating Euclidean distances d between sampling points in the intercepted parts of 2 current sampling sequencesijThen to dijRule makingThe DTW value D can be obtained by sexual accumulationij：
Dij＝dij+min{Di,j-1,Di-1,j,Di-1,j-1}
Wherein d isijRepresenting the Euclidean distance between the current amplitude of the ith sampling point of the first current sampling sequence and the current amplitude of the jth sampling point of the second current sampling sequence; dijRepresenting the ith-T of a sequence of samples taken by the first current0+1，i-T0+2, …, i sample point current amplitudes and j-T of the second sequence of current samples0+1，j-T0The DTW value between +2, …, j sample point current amplitudes.
The short window algorithm is characterized in that the interval between the peak-valley value and the zero crossing point of the fault current is 1/4 power frequency periods before and after the peak-valley value and the zero crossing point of the fault current occur, and the maximum value and the minimum value of the DTW waveform also occur at the time interval of 1/4 power frequency periods. Therefore, the waveform oscillation can be reduced by a method of first obtaining the DTW values after unitization at two sampling points near the maximum value and the minimum value, and then obtaining the average value of the DTW values:
Figure BDA0003196583120000081
Figure BDA0003196583120000082
wherein the content of the first and second substances,
Figure BDA0003196583120000083
for the first time after being optimized by a short window algorithm
Figure BDA0003196583120000084
The DTW difference of the phase is the DTW difference at two sampling points separated by 1/4 power frequency cycles
Figure BDA0003196583120000085
Average value of (d);
Figure BDA0003196583120000086
for the first time after being optimized by a short window algorithm
Figure BDA0003196583120000087
DTW braking amount of phase, which is DTW braking amount
Figure BDA0003196583120000088
Average value of (d); t is the total number of sampling points in a power frequency period, and k is a sampling serial number.
The formula used for unitizing the obtained DTW value is as follows:
Figure BDA0003196583120000089
d is a DTW actual accumulated value between the two current sampling sequences;
Figure BDA00031965831200000810
is the DTW value after the unit.
And 4, step 4: as for the starting method, as shown in fig. 5, 6 and 7, if the starting method includes a promoter sub method 1, a promoter sub method 2 and a promoter sub method 3, when the starting value of any one sub method is greater than a threshold value, a fault occurs in a power grid line, and differential protection is started (the starting method of the local terminal starts the differential protection of the local terminal, and the starting method of the opposite terminal starts the differential protection of the opposite terminal);
and 5: as for the operation method, as shown in fig. 8, 9, 10 and 11, if the setting value of any one of the submethods is greater than the threshold value in the operation submethods 1 and 2, an intra-line fault occurs, and the differential protection immediately sends an operation signal;
step 6: if packet loss occurs to the peer-to-peer data and the number of consecutive packet loss does not exceed 3, as shown in fig. 8, 9, and 11, the threshold values of 2 sub-methods of the action method may be increased to ensure that the action method does not malfunction when an out-of-area fault occurs; if the number of consecutive lost packets exceeds 3, as shown in fig. 10, the adjustment threshold value cannot be set reliably, so that the differential protection needs to be locked for a short time to prevent malfunction.
As shown in fig. 6, 8, 9, 10 and 11, the starting method sub-method 2 and the operation method sub-method 2 of the present invention calculate the DTW value and then calculate the starting value and the setting value by the short window algorithm.
The starting method comprises a starter method 1, as shown in fig. 5, dividing a three-phase current sequence at the m side of a home terminal into three groups of AB, BC and CA;
for group AB, then:
Figure BDA0003196583120000091
Figure BDA0003196583120000092
Figure BDA0003196583120000093
in the formula (I), the compound is shown in the specification,
Figure BDA0003196583120000094
for the optimized DTW difference between the two phases AB,
Figure BDA0003196583120000095
for the optimized DTW braking amount of the two AB phases through the short window algorithm,
Figure BDA0003196583120000096
and
Figure BDA0003196583120000097
are respectively as
Figure BDA0003196583120000098
And
Figure BDA0003196583120000099
the amount of change in (c);
t is the total number of sampling points in a power frequency period, and k is a sampling serial number;
if the starting value is
Figure BDA00031965831200000910
Greater than a threshold value Dop1If yes, differential protection is started;
for the BC group, then:
Figure BDA00031965831200000911
Figure BDA00031965831200000912
Figure BDA00031965831200000913
in the formula (I), the compound is shown in the specification,
Figure BDA00031965831200000914
for the optimized BC two-phase DTW difference momentum through the short window algorithm,
Figure BDA00031965831200000915
for the optimized BC two-phase DTW braking amount through the short window algorithm,
Figure BDA00031965831200000916
and
Figure BDA00031965831200000917
are respectively as
Figure BDA00031965831200000918
And
Figure BDA00031965831200000919
the amount of change in (c);
t is the total number of sampling points in a power frequency period, and k is a sampling serial number;
if the starting value is
Figure BDA00031965831200000920
Greater than a threshold value Dop1If yes, differential protection is started;
for the CA group, then:
Figure BDA00031965831200000921
Figure BDA00031965831200000922
Figure BDA00031965831200000923
in the formula (I), the compound is shown in the specification,
Figure BDA00031965831200000924
for the DTW difference between two phases of CA after the optimization of the short window algorithm,
Figure BDA00031965831200000925
for the optimized DTW braking amount between two phases of CA by the short window algorithm,
Figure BDA00031965831200000926
and
Figure BDA00031965831200000927
are respectively as
Figure BDA00031965831200000928
And
Figure BDA00031965831200000929
the amount of change in (c); t is the total number of sampling points in a power frequency period, and k is a sampling serial number;
if the starting value is
Figure BDA00031965831200000930
Greater than a threshold value Dop1Then the differential protection is enabled.
The starting method comprises a starter method 2, as shown in fig. 6, splitting a three-phase current sequence at one end into odd-even sequences respectively;
for phase A current, split it into odd sequences IA1And even sequence IA2；
IA1(k)＝IA(2k-1),k＝1,2...
-IA2(k)＝-IA(2k),k＝1,2...
Calculation of IA1and-IA2DTW value of
Figure BDA0003196583120000101
Calculation of IA1And IA2DTW value of
Figure BDA0003196583120000102
Calculated by short window algorithm
Figure BDA0003196583120000103
And
Figure BDA0003196583120000104
then there are:
Figure BDA0003196583120000105
if the starting value is
Figure BDA0003196583120000106
Greater than a threshold value Dop2If yes, differential protection is started;
for the phase B current, split it into odd sequences IB1And even sequence IB2；
IB1(k)＝IB(2k-1),k＝1,2...
-IB2(k)＝-IB(2k),k＝1,2...
Calculation of IB1and-IB2DTW value of
Figure BDA0003196583120000107
Calculation of IB1And IB2DTW value of
Figure BDA0003196583120000108
Calculated by short window algorithm
Figure BDA0003196583120000109
And
Figure BDA00031965831200001010
then there are:
Figure BDA00031965831200001011
if the starting value is
Figure BDA00031965831200001012
Greater than a threshold value Dop2If yes, differential protection is started;
for the C-phase current, split it into odd sequences IC1And even sequence IC2；
IC1(k)＝IC(2k-1),k＝1,2...
-IC2(k)＝-IC(2k),k＝1,2...
Calculation of IC1and-IC2DTW value of
Figure BDA00031965831200001013
Calculation of IC1And IC2DTW value of
Figure BDA00031965831200001014
Calculated by short window algorithm
Figure BDA00031965831200001015
And
Figure BDA00031965831200001016
then there are:
Figure BDA00031965831200001017
if the starting value is
Figure BDA00031965831200001018
Greater than a threshold value Dop2Then the differential protection is enabled.
The starting method comprises a starter method 3, as shown in fig. 7, taking the m side of the line as an example, calculating the zero sequence current I of the m sidem0Is shown bym0Divided into odd sequences Im01And even sequence Im02Calculating Im01and-Im02DTW value of
Figure BDA00031965831200001019
Calculation of Im01And Im02DTW value of
Figure BDA00031965831200001020
Recalculation
Figure BDA00031965831200001021
And
Figure BDA00031965831200001022
then there are:
Figure BDA00031965831200001023
if the left starting value is larger than the threshold value Dop3Then the differential protection is enabled.
The action method comprises an action submethod 1, as shown in fig. 8, 9 and 10, and is characterized in that two current sequences of the same phase at two ends of a line are used for obtaining a DTW difference quantity and a braking quantity, and then a short window algorithm is used for obtaining the DTW difference quantity and the braking quantity
Figure BDA0003196583120000111
And
Figure BDA0003196583120000112
the time interval between two sampling points in the short window algorithm here is 1/8 power frequency cycles. Then there are:
Figure BDA0003196583120000113
if the upper left setting value is larger than the threshold value Dop4Then the differential protection outlet action signal.
The action method comprises an action submethod 2, as shown in fig. 11, calculating zero sequence current at two ends of the line and dividing the zero sequence current into odd-even sequences, calculating DTW difference and braking quantity, and calculating by using a short window algorithm
Figure BDA0003196583120000114
And
Figure BDA0003196583120000115
the time interval between two sampling points in the short window algorithm here is 1/8 power frequency cycles. Then there are:
Figure BDA0003196583120000116
if the upper left setting value is larger than the threshold value Dop5Then the differential protection outlet action signal.
As shown in fig. 5, in steady state, the start value of the sub-method 1 of the start method is a very small positive number and is basically stable, when the line fails, at least one of the two DTW values in the maximum function corresponding to the start value will increase rapidly, and the threshold value D is takenop1Is 0.003 if the starting value is greater than Dop1Then the differential protection is enabled. Submethod 1 of the startup method can start differential protection within 2ms after the occurrence of a fault and can be used to detect various types of faults other than a high-resistance ground fault.
As shown in FIG. 6, in sub-method 2 of the start-up method, when a line fault occurs, the start-up value increases rapidly, and the threshold value D is takenop2Is 0.028, if the starting value is greater than Dop2Then the differential protection is enabled. Sub-method 2 of the start-up method can start differential protection within 8ms after the occurrence of a fault and can be used to detect various types of faults other than a single-phase earth fault.
As shown in fig. 7, sub-method 3 of the start-up method also increases the start-up value rapidly when a line fault occurs, taking threshold value Dop3Is 0.02, if the starting value is greater than Dop3Then the differential protection is enabled. Submethod 2 of the startup method can start differential protection within 2ms after the fault occurs, can be used to detect asymmetric ground faults, and is used as a complement to submethods 1 and 2 of the startup method.
As shown in fig. 8, 9, and 10, in the submethod 1 of the operation method, when an out-of-range fault occurs, the setting value is lower than the threshold value, and the differential protection does not operate. For the fault in the area, the threshold value D is taken except that the high-resistance earth fault can not be setop4And the packet loss is 0.01, and the differential protection can output the action signal within 6.75ms as long as the packet loss is not more than 3 continuously. If the number of consecutive lost packets exceeds 3, the operation speed gradually decreases.
As shown in fig. 11, if there are no more than 3 consecutive lost packets, sub-method 2 of the action method can adjust all asymmetric ground faults and can export the action signal within 6.75ms as well.
Example 2
Based on the same inventive concept as embodiment 1, the embodiment of the invention provides a starting and acting system for power grid differential protection, which comprises a storage medium and a processor;
the storage medium is used for storing instructions;
the processor is configured to operate in accordance with the instructions to perform the method of any of embodiment 1.
The foregoing shows and describes the general principles and broad features of the present invention and advantages thereof. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, which are described in the specification and illustrated only to illustrate the principle of the present invention, but that various changes and modifications may be made therein without departing from the spirit and scope of the present invention, which fall within the scope of the invention as claimed. The scope of the invention is defined by the appended claims and equivalents thereof.

Claims (12)

1. A starting and action method for power grid differential protection is characterized by comprising the following steps:
performing DTW calculation on the normalized current sampling sequence, and performing unitization processing on the obtained DTW value;
aiming at the DTW value after the unitization processing, calculating a starting value of a starting method and a setting value of an action method by combining a short window algorithm;
judging whether a fault occurs or not and whether the fault is in a line area or not based on the calculated starting value and setting value;
if the judgment result is that the power grid line has a fault, starting differential protection;
if the judgment result is that the circuit area is in fault, the differential protection is controlled to immediately output an action signal.
2. A method according to claim 1, wherein the normalized current sampling sequence is obtained by:
after receiving a phase current sampling sequence of an opposite end, processing the phase current sampling sequence of the opposite end by using a substitution processing method, wherein the substitution processing method is to use previous effective sampling point data to substitute current abnormal sampling point data;
and selecting two current sampling sequences from all the phase current sampling sequences and the zero sequence current sampling sequences at two ends of the line, and normalizing the two current sampling sequences.
3. A method according to claim 2, wherein the method further comprises: the DTW calculation specifically comprises the following steps:
intercepting the last T of the two normalized current sampling sequences0Sampling point data;
and calculating the DTW value between the sampling points of the intercepted parts of the two current sampling sequences by using a DTW algorithm.
4. A method according to claim 3, wherein the DTW value is processed in units of the formula:
Figure FDA0003196583110000011
d is a DTW actual accumulated value between the two current sampling sequences;
Figure FDA0003196583110000012
is the DTW value after the unit.
5. The method for starting and operating the power grid differential protection according to claim 2, wherein the substitution processing method adopts a calculation formula as follows:
I′i＝Ii-1,i＝2,3,4...
wherein, Ii-1Is the amplitude of the current of the ith-1 valid sampling point, I'iInstead of the magnitude of the current at the ith sample point after processing.
6. The method for starting and operating grid differential protection according to claim 1, wherein the DTW value is calculated by the formula:
Dij＝dij+min{Di,j-1,Di-1,j,Di-1,j-1}
wherein d isijRepresenting the Euclidean distance between the current amplitude of the ith sampling point of the first current sampling sequence and the current amplitude of the jth sampling point of the second current sampling sequence; dijRepresenting the ith-T of a sequence of samples taken by the first current0+1，i-T0+2, …, i sample point current amplitudes and j-T of the second sequence of current samples0+1，j-T0The DTW value between +2, …, j sample point current amplitudes.
7. The method for starting and operating the power grid differential protection according to claim 1, wherein the short window algorithm has a calculation formula as follows:
Figure FDA0003196583110000021
Figure FDA0003196583110000022
wherein the content of the first and second substances,
Figure FDA0003196583110000023
for the first time after being optimized by a short window algorithm
Figure FDA0003196583110000024
The DTW difference of the phase is the DTW difference at two sampling points separated by 1/4 power frequency cycles
Figure FDA0003196583110000025
Average value of (d);
Figure FDA0003196583110000026
for the first time after being optimized by a short window algorithm
Figure FDA0003196583110000027
DTW braking amount of phase, which is DTW braking amount
Figure FDA0003196583110000028
Average value of (d); t is the total number of sampling points in a power frequency period, and k is a sampling serial number.
8. A method according to claim 1, wherein the method further comprises the steps of: the starting method comprises a promoter method 1, a promoter method 2 and a promoter method 3;
the promoter value of the promoter method 1 was calculated by the following steps:
defining a current sampling sequence to be divided into three groups of AB, BC and CA;
for group AB, then:
Figure FDA0003196583110000029
Figure FDA00031965831100000210
Figure FDA00031965831100000211
in the formula (I), the compound is shown in the specification,
Figure FDA00031965831100000212
for the optimized DTW difference between the two phases AB,
Figure FDA00031965831100000213
for the optimized DTW braking amount of the two AB phases through the short window algorithm,
Figure FDA00031965831100000214
and
Figure FDA00031965831100000215
are respectively as
Figure FDA00031965831100000216
And
Figure FDA00031965831100000217
the amount of change in (c); t is the total number of sampling points in a power frequency period, and k is a sampling serial number;
if the starting value is
Figure FDA00031965831100000218
Greater than a threshold value Dop1If yes, differential protection is started;
for the BC group, then:
Figure FDA00031965831100000219
Figure FDA00031965831100000220
Figure FDA00031965831100000221
in the formula (I), the compound is shown in the specification,
Figure FDA0003196583110000031
for the optimized BC two-phase DTW difference momentum through the short window algorithm,
Figure FDA0003196583110000032
for the optimized BC two-phase DTW braking amount through the short window algorithm,
Figure FDA0003196583110000033
and
Figure FDA0003196583110000034
are respectively as
Figure FDA0003196583110000035
And
Figure FDA0003196583110000036
the amount of change in (c); t is the total number of sampling points in a power frequency period, and k is a sampling serial number;
if the starting value is
Figure FDA0003196583110000037
Greater than a threshold value Dop1If yes, differential protection is started;
for the CA group, then:
Figure FDA0003196583110000038
Figure FDA0003196583110000039
Figure FDA00031965831100000310
in the formula (I), the compound is shown in the specification,
Figure FDA00031965831100000311
for the DTW difference between two phases of CA after the optimization of the short window algorithm,
Figure FDA00031965831100000312
for the optimized DTW braking amount between two phases of CA by the short window algorithm,
Figure FDA00031965831100000313
and
Figure FDA00031965831100000314
are respectively as
Figure FDA00031965831100000315
And
Figure FDA00031965831100000316
the amount of change in (c); t is the total number of sampling points in a power frequency period, and k is a sampling serial number;
if the starting value is
Figure FDA00031965831100000317
Greater than a threshold value Dop1If yes, differential protection is started;
promoter method 2 promoter values were calculated by the following steps:
for phase A current, split it into odd sequences IA1And even sequence IA2；
IA1(k)＝IA(2k-1),k＝1,2...
-IA2(k)＝-IA(2k),k＝1,2...
Calculation of IA1and-IA2DTW value of
Figure FDA00031965831100000318
Calculation of IA1And IA2DTW value of
Figure FDA00031965831100000319
Calculated by short window algorithm
Figure FDA00031965831100000320
And
Figure FDA00031965831100000321
then there are:
Figure FDA00031965831100000322
if the starting value is
Figure FDA00031965831100000323
Greater than a threshold value Dop2If yes, differential protection is started;
for the phase B current, split it into odd sequences IB1And even sequence IB2；
IB1(k)＝IB(2k-1),k＝1,2...
-IB2(k)＝-IB(2k),k＝1,2...
Calculation of IB1and-IB2DTW value of
Figure FDA00031965831100000324
Calculation of IB1And IB2DTW value of
Figure FDA00031965831100000325
Calculated by short window algorithm
Figure FDA00031965831100000326
And
Figure FDA00031965831100000327
then there are:
Figure FDA00031965831100000328
if the starting value is
Figure FDA00031965831100000329
Greater than a threshold value Dop2If yes, differential protection is started;
for the C-phase current, split it into odd sequences IC1And even sequence IC2；
IC1(k)＝IC(2k-1),k＝1,2...
-IC2(k)＝-IC(2k),k＝1,2...
Calculation of IC1and-IC2DTW value of
Figure FDA0003196583110000041
Calculation of IC1And IC2DTW value of
Figure FDA0003196583110000042
Calculated by short window algorithm
Figure FDA0003196583110000043
And
Figure FDA0003196583110000044
then there are:
Figure FDA0003196583110000045
if the starting value is
Figure FDA0003196583110000046
Greater than a threshold value Dop2If yes, differential protection is started;
the promoter value for promoter method 3 was calculated by the following steps:
adding the three-phase current sampling sequences and dividing by 3 to obtain zero-sequence current Im0Is shown bym0Divided into odd sequences Im01And even sequence Im02Calculating Im01and-Im02DTW value of
Figure FDA0003196583110000047
Calculation of Im01And Im02DTW value of
Figure FDA0003196583110000048
Recalculation
Figure FDA0003196583110000049
And
Figure FDA00031965831100000410
then there are:
Figure FDA00031965831100000411
if the starting value is
Figure FDA00031965831100000412
Greater than a threshold value Dop3Then the differential protection is started, m represents one of the two ends of the line, and 0 represents zero sequence.
9. A method according to claim 8, wherein the method further comprises: if the judgment result is that the power grid line fails, starting differential protection, specifically:
and if the starting value of any one of the method 1, 2 and 3 is greater than the threshold value, judging that the power grid line has a fault, and starting differential protection.
10. A method according to claim 1, wherein the method further comprises the steps of: the action method comprises an action submethod 1 and an action submethod 2;
the setting value of the action sub-method 1 is calculated by the following steps:
obtaining the DTW differential quantity optimized by the short window algorithm by using two current sampling sequences of the same phase at two ends of the line
Figure FDA00031965831100000413
And amount of braking
Figure FDA00031965831100000414
The time interval of two sampling points in the short window algorithm is 1/8 power frequency periods, and then:
Figure FDA00031965831100000415
if the setting value
Figure FDA00031965831100000416
Greater than a threshold value Dop4If so, the differential protection action is carried out, K is a braking coefficient and belongs to [0,1]](ii) a The setting value of the action submethod 2 is measured by the following stepsCalculating:
calculating zero sequence current at two ends of the line, dividing the zero sequence current into odd-even sequences, and obtaining the DTW differential quantity optimized by the short window algorithm
Figure FDA00031965831100000417
And amount of braking
Figure FDA00031965831100000418
The time interval of two sampling points in the short window algorithm is 1/8 power frequency periods, and then:
Figure FDA00031965831100000419
if the setting value
Figure FDA0003196583110000051
Greater than a threshold value Dop5Then, differential protection is performed.
11. The method according to claim 1, wherein if the determination result is an internal fault of the line area, the method controls an immediate exit operation signal of the differential protection, specifically:
if the setting value of any one of the sub-methods 1 and 2 is larger than the threshold value, the occurrence of the fault in the line area is judged, and the differential protection immediately acts.
12. A starting and action system for power grid differential protection is characterized by comprising a storage medium and a processor;
the storage medium is used for storing instructions;
the processor is configured to operate in accordance with the instructions to perform the method of any of claims 1-11.
CN202110892722.2A 2021-08-04 2021-08-04 Starting and acting method and system for power grid differential protection Active CN113690859B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110892722.2A CN113690859B (en) 2021-08-04 2021-08-04 Starting and acting method and system for power grid differential protection

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110892722.2A CN113690859B (en) 2021-08-04 2021-08-04 Starting and acting method and system for power grid differential protection

Publications (2)

Publication Number Publication Date
CN113690859A true CN113690859A (en) 2021-11-23
CN113690859B CN113690859B (en) 2023-02-21

Family

ID=78578794

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110892722.2A Active CN113690859B (en) 2021-08-04 2021-08-04 Starting and acting method and system for power grid differential protection

Country Status (1)

Country Link
CN (1) CN113690859B (en)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU736254A1 (en) * 1977-03-09 1980-05-25 Предприятие П/Я В-2616 Device for short-circuiting protection of ac generator of two-channel power supply system of aircraft in distributing network
CN104078948A (en) * 2014-07-15 2014-10-01 国家电网公司 Transmission line positive sequence current split-phase differential protection method based on short data window
CN105353271A (en) * 2015-11-03 2016-02-24 云南电网有限责任公司电力科学研究院 Current differential protection saturation distinguishing method
CN105529688A (en) * 2016-02-25 2016-04-27 三峡大学 Transformer excitation inrush current and fault differential current recognition method based on Hausdorff distance algorithm
CN108562820A (en) * 2018-05-04 2018-09-21 国网甘肃省电力公司电力科学研究院 Single-phase ground fault line selecting method of small-electric current grounding system and device
WO2020144600A1 (en) * 2019-01-08 2020-07-16 Panoramic Power Ltd. A methods for determining an alternate current motor fault in a non-variable frequency device based on current analysis
CN111478297A (en) * 2020-04-30 2020-07-31 深圳供电局有限公司 Calculation method for power grid differential protection
CN111987687A (en) * 2020-07-14 2020-11-24 许继集团有限公司 Quick-break differential protection method for preventing surge interference unwanted operation
CN112054482A (en) * 2020-08-06 2020-12-08 三峡大学 Converter transformer zero sequence differential protection anti-misoperation method based on DTW algorithm
CN113013852A (en) * 2021-02-26 2021-06-22 许继集团有限公司 Line protection optical fiber channel transmission data processing method and device
CN113034305A (en) * 2021-02-10 2021-06-25 上海千居智科技有限公司 Non-invasive load monitoring event classification method and storage medium

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU736254A1 (en) * 1977-03-09 1980-05-25 Предприятие П/Я В-2616 Device for short-circuiting protection of ac generator of two-channel power supply system of aircraft in distributing network
CN104078948A (en) * 2014-07-15 2014-10-01 国家电网公司 Transmission line positive sequence current split-phase differential protection method based on short data window
CN105353271A (en) * 2015-11-03 2016-02-24 云南电网有限责任公司电力科学研究院 Current differential protection saturation distinguishing method
CN105529688A (en) * 2016-02-25 2016-04-27 三峡大学 Transformer excitation inrush current and fault differential current recognition method based on Hausdorff distance algorithm
CN108562820A (en) * 2018-05-04 2018-09-21 国网甘肃省电力公司电力科学研究院 Single-phase ground fault line selecting method of small-electric current grounding system and device
WO2020144600A1 (en) * 2019-01-08 2020-07-16 Panoramic Power Ltd. A methods for determining an alternate current motor fault in a non-variable frequency device based on current analysis
CN111478297A (en) * 2020-04-30 2020-07-31 深圳供电局有限公司 Calculation method for power grid differential protection
CN111987687A (en) * 2020-07-14 2020-11-24 许继集团有限公司 Quick-break differential protection method for preventing surge interference unwanted operation
CN112054482A (en) * 2020-08-06 2020-12-08 三峡大学 Converter transformer zero sequence differential protection anti-misoperation method based on DTW algorithm
CN113034305A (en) * 2021-02-10 2021-06-25 上海千居智科技有限公司 Non-invasive load monitoring event classification method and storage medium
CN113013852A (en) * 2021-02-26 2021-06-22 许继集团有限公司 Line protection optical fiber channel transmission data processing method and device

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李海林: "时间序列数据挖掘中的动态时间弯曲研究综述", 《控制与决策》 *

Also Published As

Publication number Publication date
CN113690859B (en) 2023-02-21

Similar Documents

Publication Publication Date Title
CN105656038B (en) A kind of meter and power network receive the probability load flow calculation method of ability
CN101562332A (en) Self-adaptive current fast tripping protection method for distribution network comprising inverse distributed power
CN109038664B (en) High-frequency turn setting method considering fan grid-involvement protection
WO2022242020A1 (en) Method and system for optimizing direct current transmission limit within planning period
CN104361068B (en) Parallel method of partition and system during a kind of data deduplication
CN113690859B (en) Starting and acting method and system for power grid differential protection
CN115877129B (en) Low-voltage distribution network neutral line disconnection comprehensive judgment method based on current and voltage
CN103107530A (en) Multi-direct-current urgent coordination control method based on sensitivity
CN111541280A (en) Power grid wind power maximum permeability evaluation method considering static voltage stability constraint
CN108736921A (en) A kind of power line carrier communication lead detecting method for resisting random impulsive noise
CN110149221A (en) A kind of construction method and device of the adaptive network based on alliance's chain
CN109861855B (en) Method and device for determining importance of nodes in power communication network
Peng et al. A new approximation algorithm for computing 2-restricted disjoint paths
EP2476211B1 (en) A method for estimating transmission properties of a telecommunications transmission line
CN112102611B (en) Urban road network key path identification method based on interruption probability influence model
Luo et al. Single-Ended Time Domain Fault Location Based on Transient Signal Measurements of Transmission Lines
CN114725903B (en) Self-adaptive full-line current protection method, device, medium and system for power distribution network
CN112561137B (en) Power grid high-frequency generator tripping optimization method for new energy unit
CN111541266B (en) Load flow calculation method of power distribution system with distributed power supply based on forward-backward substitution idea
Deshpande et al. A statistical approach to anomaly detection in interdomain routing
Sun et al. A new routing selection method considering network capacity and path cost
Ahuja Performance based reliability Optimization for computer networks
CN109578208A (en) Wind power generating set and its Yaw control method and device
CN117880208A (en) Speed limiting method for realizing lock-free speed limiting by multi-core dynamic speed adjustment
Chen et al. An Improved Backward/Forward Algorithm for Distribution Network Power Flow Calculation Considering Meteorological Factors

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