CN111812454B - Automatic alignment and correction method and system based on wave recording data - Google Patents
Automatic alignment and correction method and system based on wave recording data Download PDFInfo
- Publication number
- CN111812454B CN111812454B CN202010592547.0A CN202010592547A CN111812454B CN 111812454 B CN111812454 B CN 111812454B CN 202010592547 A CN202010592547 A CN 202010592547A CN 111812454 B CN111812454 B CN 111812454B
- Authority
- CN
- China
- Prior art keywords
- sampling
- point
- value
- fault
- correction method
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/088—Aspects of digital computing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R35/00—Testing or calibrating of apparatus covered by the other groups of this subclass
- G01R35/005—Calibrating; Standards or reference devices, e.g. voltage or resistance standards, "golden" references
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
- Y04S10/52—Outage or fault management, e.g. fault detection or location
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Locating Faults (AREA)
Abstract
The invention relates to an automatic alignment and correction method and system based on wave recording data, and belongs to the technical field of relay protection of power systems. The method comprises the steps of firstly, primarily selecting a fault starting time point T0 by judging the position of a sudden change point, then, reversely deducing and correcting the T0 of wave recording data on two sides by utilizing the principle that the peak time on the two sides is the same through the time difference between the T0 and the peak time of adjacent voltage and current, and accurately determining the double-end distance measurement synchronous alignment point. The method realizes synchronous alignment of the wave recording data at two sides, overcomes the influence of common factors such as zero crossing point, gradual change development, large harmonic wave and distortion and the like on the simple judgment of the fault starting point, eliminates the time deviation possibly generated by alignment by directly judging the fault starting point, and is beneficial to improving the accuracy and reliability of double-end distance measurement.
Description
Technical Field
The invention belongs to the technical field of relay protection of power systems, and particularly relates to an automatic alignment and correction method and system based on wave recording data.
Background
The key factor of double-end distance measurement based on the wave recording data is that accurate synchronous alignment points are found by the wave recording data on the two sides. The existing method for synchronously aligning the wave recording data on two sides mainly comprises a self-contained time mark alignment method and an alignment method for judging the starting moment of a fault. The self-contained time mark alignment method is limited by external factors such as clock synchronization precision and the like, so that the risk is extremely high and the reliability is poor; and the alignment method for judging the fault starting time has the advantages of strong autonomy and flexibility, and is more suitable for practical application.
However, there are also unstable factors in determining the fault starting time, such as zero crossing point, gradual change of development, large harmonic and distortion of the fault starting time, which all affect the determination of the fault starting time to different degrees. Taking the recording data with a sampling rate of 1000Hz as an example, the judgment error of only 1 sampling point will generate 18 degrees phase angle deviation of the data on both sides, which will have serious influence on the double-end ranging result. Therefore, how to overcome the defects of the prior art is a problem which needs to be solved in the technical field of power system relay protection at present.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides a method and a system for automatically judging a fault starting time point based on wave recording data and realizing alignment correction by using adjacent wave crests. The method selects the most obvious and unique peak characteristic signal in the recording sampling curve as the alignment reference points at two sides, reversely deduces and corrects the fault starting time point T0 by using the principle that the time difference from the peak value to the T0 is the same, realizes synchronous alignment of the recording data at two sides, overcomes the influence of common factors such as zero crossing point, gradual change development, harmonic wave, large distortion and the like which are possibly applied to the simple judgment of the fault starting point, eliminates the time deviation which is possibly generated by directly judging the fault starting time point for alignment, and is beneficial to improving the accuracy and the reliability of double-end distance measurement.
In order to achieve the purpose, the technical scheme adopted by the invention is as follows:
an automatic alignment and correction method based on recorded wave data comprises the following steps:
step (1), judging a mutation point T0 by a single-phase sampling instantaneous value curve:
(1.1) on a sampling current instantaneous value curve, sampling points are sequentially taken from the ith sampling point to the right along the T axis at intervals of T from the ith sampling point, and 3 sampling points are taken in total and are respectively marked as Pi and Pi + T, pi +2T; namely, the ith sampling point is taken as the first point; the corresponding sampling current instantaneous values are respectively marked as V1, V2 and V3;
(1.2) calculating the instantaneous value of the current at three sampled points by the following steps
i) The difference value of the V1 and the V2 is recorded as det1, then det1= | V1-V2|;
ii) the difference between V2 and V3 is denoted as det2, then det2= | V2-V3|;
iii) The difference value of det1 and det12 is marked as detV, and then detV = | det1-det2|;
(1.3) if detV <0.1, the following judgment is performed:
i) If i >15T, T0= ∞, abandoning T0 of the local side to wait for the correction of the contralateral side in the step (2);
ii) if not, taking the next point from i to the right along the t axis, and repeating the steps (1.1) - (1.3);
(1.4) otherwise, if detV >0.1, T0= i +2T;
step (2), three-phase current and zero-phase current catastrophe point comprehensive T0:
adopting the method in the step (1) to sequentially calculate the mutation points of the A, B, C three-phase current instantaneous value curves Ia, ib and Ic, simultaneously calculating the mutation points of the zero-phase self-generated current instantaneous value curve I0, and respectively marking 4 calculated results as T0_ Ia, T0_ Ib, T0_ Ic and T0_ I0;
wherein I0= (Ia + Ib + Ic);
the sum T0 takes the minimum of 4 results, i.e.
Integrated T0= min { T0_ Ia, T0_ Ib, T0_ Ic, T0_ I0};
step (3), setting the two sides of the fault line as an M side and an N side respectively, and then correcting the starting time point of the fault at the other side by taking the comprehensive T0 point at the M side obtained in the step (2) as a reference or taking the comprehensive T0 point at the N side as a reference; here, the M-side integrated T0 point is denoted as MT0, and the N-side integrated T0 point is denoted as NT0.
Further, preferably, in the step (3), the N-side fault start time point NT0 is corrected with reference to the M-side fault start time point MT0, and a corrected value NT0' of NT0 is set, where the fault phase is X, and the fault voltage and the fault current are Ux and Ix; x is A, B or C; the method comprises the following steps:
(3.1.1) curve correction method based on I0 sampling instantaneous value:
NT0’=NT0+(NdetT_I0-MdetT_I0)；
(3.1.2) curve correction method based on Ix sampling instantaneous value:
NT0’=NT0-(NdetT_Ix-MdetT_Ix)；
(3.1.3) instantaneous value curve correction method based on Ux sampling:
NT0’=NT0-(NdetT_Ux-MdetT_Ux)；
wherein MdetT _ Ux represents the time difference of the first peak value on the left side of T0 between M side T0 and Ux;
NdetT _ Ux represents the time difference between the first peak value on the left side of T0 of the N side T0 and Ux;
MdetT _ Ix represents the time difference of the first peak value on the left side of T0 between M side T0 and Ix;
NdetT _ Ix represents the time difference of the first peak value of the N side T0 and Ix on the left side of T0;
MdetT _ I0 represents the time difference of the first peak value of the M side T0 and I0 on the right side of T0;
NdetT _ I0 represents the time difference of the first peak value of the N side T0 and the right side I0 of T0;
the logic of the use conditions of the correction methods (3.1.1) - (3.1.3) is as follows:
a) At the time of T0, if the secondary fundamental wave effective value of I0 is greater than 1/2 of the Ix fundamental wave effective value, using a curve correction method based on I0 sampling instantaneous values;
b) Otherwise, if the effective value of the secondary fundamental wave of Ix is greater than 0.1A, using a curve correction method based on Ix sampling instantaneous values;
c) Otherwise, curve correction based on Ux sampling transients is used.
Further, preferably, the data windows used for I0 are T0 right side T length, and the data windows used for Ix are T0 left side T length, for the fundamental effective values a) and b).
Further, in the step (3), preferably, the N-side fault start time point NT0 is used as a reference, the M-side fault start time point MT0 is corrected, the corrected value of MT0 is MT0', the fault phase is X, and the fault voltage and the fault current are Ux and Ix; x is A, B or C; the method comprises the following steps:
(3.2.1) curve correction method based on I0 sampling instantaneous value:
MT0’=MT0+(MdetT_I0-NdetT_I0)；
(3.2.2) curve correction method based on Ix sampling instantaneous value:
MT0’=MT0-(MdetT_Ix-NdetT_Ix)；
(3.2.3) instantaneous value curve correction method based on Ux sampling:
MT0’=MT0-(MdetT_Ux-NdetT_Ux)；
wherein MdetT _ Ux represents the time difference of the first peak value on the left side of T0 between M side T0 and Ux;
NdetT _ Ux represents the time difference between the first peak value on the left side of T0 of the N side T0 and Ux;
MdetT _ Ix represents the time difference of the first peak value on the left side of T0 between M side T0 and Ix;
NdetT _ Ix represents the time difference of the first peak value of the N side T0 and Ix on the left side of T0;
MdetT _ I0 represents the time difference of the first peak value of the M side T0 and I0 on the right side of T0;
NdetT _ I0 represents the time difference of the first peak value of the N side T0 and the right side I0 of T0;
the logic of the use conditions of the correction methods (3.2.1) - (3.2.3) is as follows:
(a) At the time of T0, if the secondary fundamental wave effective value of I0 is greater than 1/2 of the Ix fundamental wave effective value, using a curve correction method based on I0 sampling instantaneous values;
(b) Otherwise, if the secondary fundamental wave effective value of Ix is greater than 0.1A, using a curve correction method based on Ix sampling instantaneous values;
(c) Otherwise, curve rectification based on Ux sampling is used.
Further, it is preferable that the (a) and (b) fundamental effective values use data windows with a T length on the right side of T0 for I0, and a T length on the left side of T0 for Ix.
The invention also provides an automatic aligning and correcting system based on the wave recording data, which comprises the following components:
the data acquisition module is used for sequentially taking sampling points from the ith sampling point to the right along the T axis by taking T as an interval from the ith sampling point on a sampling current instantaneous value curve, and taking 3 sampling points in total, wherein the corresponding sampling current instantaneous values are respectively marked as V1, V2 and V3;
the first processing module is used for calculating the catastrophe points of A, B, C three-phase current instantaneous value curves Ia, ib and Ic according to the current instantaneous values of the obtained sampling points, and calculating the catastrophe points of a zero-phase self-generated current instantaneous value curve I0 to obtain comprehensive T0;
and the recording data alignment correction module is used for correcting the fault starting time point on the other side according to the obtained comprehensive T0 point on any one side of the two sides of the fault line as a reference.
The invention also provides electronic equipment which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor executes the program to realize the steps of the automatic alignment and correction method based on the recording data.
The present invention further provides a non-transitory computer-readable storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of the above-described recording data-based automatic alignment and correction method.
Setting the two sides of a fault line as an M side and an N side respectively, namely stations at two ends of a power transmission line respectively; the station is a transformer substation or a power plant, and a wave recorder and a protection device are arranged in the station.
Compared with the prior art, the invention has the following beneficial effects:
(1) An uncertain area at the starting moment of the fault is avoided, and instead, the adjacent peak value is directly positioned, the peak value belongs to a determined and unique signal, and the reliability is high;
(2) The method is not influenced by the difference of sampling rates of wave recording data on two sides; the recording data at two sides can independently judge the starting time point of the fault by using the inherent sampling rate of the recording data, respectively calculate the time difference between the starting point of the fault and the adjacent peak value, and the two sides can finish the correction of the respective starting points of the fault only by comparing the time difference;
(3) The fault tolerance rate of the fault judgment on the fault starting time point is high; the method can bear the deviation of at most 1/4 cycle of the starting point of the fault by automatically correcting the starting point of the fault through adjacent peak values;
(4) The correction certainty is high; three levels of zero-phase current, phase current and phase voltage are adopted for correction, so that the accuracy is ensured to the maximum extent, and the reliability is also ensured to 100%;
(5) The method is suitable for realizing computer programming and realizes the double-end ranging target without human intervention.
Drawings
FIG. 1 is a schematic diagram of a fault recording sampling instantaneous value curve; wherein, ua: a phase voltage instantaneous value curve; ia: a phase current instantaneous value curve; i0: a zero-phase self-generated current transient curve; t: sampling period; t0: a fault starting time point, namely an alignment reference point; detT _ Ua: the time difference of the first peak value on the left side of T0 between T0 and Ua; detT _ Ia: the time difference between the first peak value on the left side of T0 and the first peak value of Ia; detT _ I0: the time difference between the first peak value of T0 and the first peak value of I0 on the right side of T0;
FIG. 2 is a schematic diagram illustrating the sampling point positions for determining the mutation point T0;
FIG. 3 is a schematic diagram of the A, B, C three phase current and self-generated zero phase current;
FIG. 4 is a schematic structural diagram of an automatic alignment and correction system based on recorded wave data according to the present invention;
fig. 5 is a schematic structural diagram of an electronic device according to the present invention.
Detailed Description
The present invention will be described in further detail with reference to examples.
It will be appreciated by those skilled in the art that the following examples are illustrative of the invention only and should not be taken as limiting the scope of the invention. Those skilled in the art will recognize that the specific techniques or conditions, not specified in the examples, are according to the techniques or conditions described in the literature of the art or according to the product specification. The materials or equipment used are not indicated by manufacturers, and all are conventional products available by purchase.
Example 1
An automatic alignment and correction method based on recorded wave data comprises the following steps:
step (1), judging a mutation point T0 by a single-phase sampling instantaneous value curve:
(1.1) on a sampling current instantaneous value curve, sampling points are sequentially taken from the ith sampling point to the right along the T axis at intervals of T from the ith sampling point, and 3 sampling points are taken in total and are respectively marked as Pi and Pi + T, pi +2T; namely, the ith sampling point is taken as the first point; the corresponding sampling current instantaneous values are respectively marked as V1, V2 and V3;
(1.2) calculating the instantaneous value of the current at three sampled points by the following steps
i) The difference value of V1 and V2 is recorded as det1, then det1= | V1-V2|;
ii) the difference between V2 and V3 is denoted as det2, then det2= | V2-V3|;
iii) The difference value of det1 and det12 is marked as detV, and then detV = | det1-det2|;
(1.3) if detV <0.1, the following judgment is performed:
i) If i is greater than 15T, T0= ∞, abandoning the T0 of the local side to wait for the correction of the side in the step (2);
ii) if not, taking the next point from i to the right along the t axis, and repeating the steps (1.1) - (1.3);
(1.4) otherwise, if detV >0.1, T0= i +2T;
step (2), three-phase current and zero-phase current catastrophe point comprehensive T0:
calculating the mutation points of A, B, C three-phase current instantaneous value curves Ia, ib and Ic in sequence by adopting the method in the step (1), calculating the mutation points of a zero-phase self-generated current instantaneous value curve I0 at the same time, and respectively marking 4 calculated results as T0_ Ia, T0_ Ib, T0_ Ic and T0_ I0;
wherein I0= (Ia + Ib + Ic);
the sum T0 takes the minimum of 4 results, i.e.
Integrated T0= min { T0_ Ia, T0_ Ib, T0_ Ic, T0_ I0};
step (3), setting the two sides of the fault line as an M side and an N side respectively, and then correcting the starting time point of the fault at the other side by taking the comprehensive T0 point at the M side obtained in the step (2) as a reference or taking the comprehensive T0 point at the N side as a reference; the M-side integrated T0 point is denoted as MT0, and the N-side integrated T0 point is denoted as NT0.
Example 2
An automatic alignment and correction method based on recording data comprises the following steps:
step (1), judging a mutation point T0 by a single-phase sampling instantaneous value curve:
(1.1) on a sampling current instantaneous value curve, sampling points are sequentially taken from the ith sampling point to the right along the T axis at intervals of T from the ith sampling point, and 3 sampling points are taken in total and are respectively marked as Pi and Pi + T, pi +2T; namely, the ith sampling point is taken as the first point; the corresponding sampling current instantaneous values are respectively marked as V1, V2 and V3;
(1.2) calculating the instantaneous value of the current at three sampled points by the following steps
i) The difference value of the V1 and the V2 is recorded as det1, then det1= | V1-V2|;
ii) the difference between V2 and V3 is denoted as det2, then det2= | V2-V3|;
iii) The difference value of det1 and det12 is marked as detV, and then detV = | det1-det2|;
(1.3) if detV <0.1, the following judgment is performed:
i) If i is greater than 15T, T0= ∞, abandoning the T0 of the local side to wait for the correction of the side in the step (2);
ii) if not, taking the next point from i to the right along the t axis, and repeating the steps (1.1) - (1.3);
(1.4) otherwise, if detV >0.1, T0= i +2T;
step (2), three-phase current and zero-phase current catastrophe point comprehensive T0:
adopting the method in the step (1) to sequentially calculate the mutation points of the A, B, C three-phase current instantaneous value curves Ia, ib and Ic, simultaneously calculating the mutation points of the zero-phase self-generated current instantaneous value curve I0, and respectively marking 4 calculated results as T0_ Ia, T0_ Ib, T0_ Ic and T0_ I0;
wherein I0= (Ia + Ib + Ic);
the sum T0 takes the minimum of 4 results, i.e.
Integrated T0= min { T0_ Ia, T0_ Ib, T0_ Ic, T0_ I0};
step (3), setting the two sides of the fault line as an M side and an N side respectively, and then correcting the starting time point of the fault at the other side by taking the comprehensive T0 point at the M side obtained in the step (2) as a reference or taking the comprehensive T0 point at the N side as a reference; here, the M-side integrated T0 point is denoted as MT0, and the N-side integrated T0 point is denoted as NT0.
In the step (3), correcting the N-side fault starting time point NT0 by taking the M-side fault starting time point MT0 as a reference, setting a value NT0' after NT0 correction, wherein the fault phase is X, and the fault voltage and the fault current are Ux and Ix; x is A, B or C; the method comprises the following steps:
(3.1.1) curve correction method based on I0 sampling instantaneous value:
NT0’=NT0+(NdetT_I0-MdetT_I0)；
(3.1.2) curve correction method based on Ix sampling instantaneous value:
NT0’=NT0-(NdetT_Ix-MdetT_Ix)；
(3.1.3) instantaneous value curve correction method based on Ux sampling:
NT0’=NT0-(NdetT_Ux-MdetT_Ux)；
wherein MdetT _ Ux represents the time difference of the first peak value on the left side of T0 between M side T0 and Ux;
NdetT _ Ux represents the time difference of the first peak value of the N side T0 and Ux on the left side of T0;
MdetT _ Ix represents the time difference of the first peak value on the left side of T0 between M side T0 and Ix;
NdetT _ Ix represents the time difference of the first peak value of the N side T0 and Ix on the left side of T0;
MdetT _ I0 represents the time difference of the first peak values of the M side T0 and I0 on the right side of T0;
NdetT _ I0 represents the time difference of the first peak value of the N side T0 and the right side I0 of T0;
the logic of the use conditions of the correction methods (3.1.1) - (3.1.3) is as follows:
a) At the time of T0, if the secondary fundamental wave effective value of I0 is greater than 1/2 of the Ix fundamental wave effective value, using a curve correction method based on I0 sampling instantaneous values;
b) Otherwise, if the secondary fundamental wave effective value of Ix is greater than 0.1A, using a curve correction method based on Ix sampling instantaneous values;
c) Otherwise, curve correction based on Ux sampling transients is used.
Correcting the M-side fault starting time point MT0 by taking the N-side fault starting time point NT0 as a reference, setting the corrected value of MT0 as MT0', setting the fault phase as X, and setting the fault voltage and the fault current as Ux and Ix; x is A, B or C; the method comprises the following steps:
(3.2.1) curve correction method based on I0 sampling instantaneous value:
MT0’=MT0+(MdetT_I0-NdetT_I0)；
(3.2.2) curve correction method based on Ix sampling instantaneous value:
MT0’=MT0-(MdetT_Ix-NdetT_Ix)；
(3.2.3) instantaneous value curve correction method based on Ux sampling:
MT0’=MT0-(MdetT_Ux-NdetT_Ux)；
wherein MdetT _ Ux represents the time difference of the first peak value on the left side of T0 between M side T0 and Ux;
NdetT _ Ux represents the time difference between the first peak value on the left side of T0 of the N side T0 and Ux;
MdetT _ Ix represents the time difference of the first peak value on the left side of T0 between M side T0 and Ix;
NdetT _ Ix represents the time difference of the first peak value of the N side T0 and Ix on the left side of T0;
MdetT _ I0 represents the time difference of the first peak value of the M side T0 and I0 on the right side of T0;
NdetT _ I0 represents the time difference of the first peak value of the N side T0 and the right side I0 of T0;
the logic of the use conditions of the correction methods (3.2.1) - (3.2.3) is as follows:
(a) At the time of T0, if the secondary fundamental wave effective value of I0 is greater than 1/2 of the Ix fundamental wave effective value, using a curve correction method based on I0 sampling instantaneous values;
(b) Otherwise, if the secondary fundamental wave effective value of Ix is greater than 0.1A, using a curve correction method based on Ix sampling instantaneous values;
(c) Otherwise, curve rectification based on Ux sampling is used.
Example 3
An automatic alignment and correction method based on recorded wave data comprises the following steps:
step (1), judging a mutation point T0 by a single-phase sampling instantaneous value curve:
(1.1) on a sampling current instantaneous value curve, taking sampling points from the ith sampling point to the right along the T axis at intervals of T from the ith sampling point, and taking 3 sampling points in total and respectively marking as Pi and Pi + T, pi +2T; namely, the ith sampling point is taken as the first point; the corresponding sampling current instantaneous values are respectively marked as V1, V2 and V3;
(1.2) calculating the instantaneous value of the current at the three sampled points by the following steps
i) The difference value of the V1 and the V2 is recorded as det1, then det1= | V1-V2|;
ii) the difference between V2 and V3 is denoted as det2, then det2= | V2-V3|;
iii) The difference value of det1 and det12 is marked as detV, and then detV = | det1-det2|;
(1.3) if detV <0.1, the following judgment is performed:
i) If i is greater than 15T, T0= ∞, abandoning the T0 of the local side to wait for the correction of the side in the step (2);
ii) if not, taking the next point from i to the right along the t axis, and repeating the steps (1.1) - (1.3);
(1.4) otherwise, if detV >0.1, T0= i +2T;
step (2), three-phase current and zero-phase current catastrophe point comprehensive T0:
calculating the mutation points of A, B, C three-phase current instantaneous value curves Ia, ib and Ic in sequence by adopting the method in the step (1), calculating the mutation points of a zero-phase self-generated current instantaneous value curve I0 at the same time, and respectively marking 4 calculated results as T0_ Ia, T0_ Ib, T0_ Ic and T0_ I0;
wherein I0= (Ia + Ib + Ic);
the sum T0 takes the minimum of 4 results, i.e.
Integrated T0= min { T0_ Ia, T0_ Ib, T0_ Ic, T0_ I0};
step (3), setting the two sides of the fault line as an M side and an N side respectively, and then correcting the starting time point of the fault at the other side by taking the comprehensive T0 point at the M side obtained in the step (2) as a reference or taking the comprehensive T0 point at the N side as a reference; the M-side integrated T0 point is denoted as MT0, and the N-side integrated T0 point is denoted as NT0.
In the step (3), correcting the N-side fault starting time point NT0 by taking the M-side fault starting time point MT0 as a reference, setting a corrected value NT0' of NT0, wherein the fault phase is X, and the fault voltage and the fault current are Ux and Ix; x is A, B or C; the method comprises the following steps:
(3.1.1) curve correction method based on I0 sampling instantaneous value:
NT0’=NT0+(NdetT_I0-MdetT_I0)；
(3.1.2) curve correction method based on Ix sampling instantaneous value:
NT0’=NT0-(NdetT_Ix-MdetT_Ix)；
(3.1.3) instantaneous value curve correction method based on Ux sampling:
NT0’=NT0-(NdetT_Ux-MdetT_Ux)；
wherein MdetT _ Ux represents the time difference of the first peak value on the left side of T0 between M side T0 and Ux;
NdetT _ Ux represents the time difference between the first peak value on the left side of T0 of the N side T0 and Ux;
MdetT _ Ix represents the time difference of the first peak values of the M side T0 and Ix on the left side of T0;
NdetT _ Ix represents the time difference of the first peak value of the N side T0 and Ix on the left side of T0;
MdetT _ I0 represents the time difference of the first peak value of the M side T0 and I0 on the right side of T0;
NdetT _ I0 represents the time difference of the first peak value of the N side T0 and the right side I0 of T0;
the logic of the use conditions of the correction methods (3.1.1) - (3.1.3) is as follows:
a) At the time of T0, if the secondary fundamental wave effective value of I0 is greater than 1/2 of the Ix fundamental wave effective value, using a curve correction method based on I0 sampling instantaneous values;
b) Otherwise, if the secondary fundamental wave effective value of Ix is greater than 0.1A, using a curve correction method based on Ix sampling instantaneous values;
c) Otherwise, curve correction based on Ux sampling transients is used.
Correcting the M-side fault starting time point MT0 by taking the N-side fault starting time point NT0 as a reference, setting the corrected value of MT0 as MT0', setting the fault phase as X, and setting the fault voltage and the fault current as Ux and Ix; x is A, B or C; the method comprises the following steps:
(3.2.1) curve correction method based on I0 sampling instantaneous value:
MT0’=MT0+(MdetT_I0-NdetT_I0)；
(3.2.2) curve correction method based on Ix sampling instantaneous value:
MT0’=MT0-(MdetT_Ix-NdetT_Ix)；
(3.2.3) instantaneous value curve correction method based on Ux sampling:
MT0’=MT0-(MdetT_Ux-NdetT_Ux)；
wherein MdetT _ Ux represents the time difference of the first peak value on the left side of T0 between M side T0 and Ux;
NdetT _ Ux represents the time difference of the first peak value of the N side T0 and Ux on the left side of T0;
MdetT _ Ix represents the time difference of the first peak value on the left side of T0 between M side T0 and Ix;
NdetT _ Ix represents the time difference of the first peak value of the N side T0 and Ix on the left side of T0;
MdetT _ I0 represents the time difference of the first peak value of the M side T0 and I0 on the right side of T0;
NdetT _ I0 represents the time difference of the first peak value of the N side T0 and the right side I0 of T0;
the logic of the use conditions of the correction methods (3.2.1) - (3.2.3) is as follows:
(a) At the time of T0, if the effective value of the secondary fundamental wave of I0 is greater than 1/2 of the effective value of the Ix fundamental wave, using a curve correction method based on I0 sampling instantaneous values;
(b) Otherwise, if the secondary fundamental wave effective value of Ix is greater than 0.1A, using a curve correction method based on Ix sampling instantaneous values;
(c) Otherwise, curve correction based on Ux sampling transients is used.
Wherein, (a) and (b) the fundamental wave effective value, the data window used by I0 is T length on right side of T0, and the data window used by Ix is T length on left side of T0.
As shown in fig. 4, an automatic aligning and correcting system based on recording data includes:
the data acquisition module 101 is configured to sequentially take 3 sampling points from the ith sampling point to the right along a T axis at intervals of T on a sampling current instantaneous value curve from the ith sampling point, where corresponding sampling current instantaneous values are respectively marked as V1, V2, and V3;
the first processing module 102 is configured to perform catastrophe point calculation on the A, B, C three-phase current instantaneous value curves Ia, ib and Ic according to the current instantaneous value of the obtained sampling point, and calculate a catastrophe point of the zero-phase self-generated current instantaneous value curve I0 at the same time to obtain a comprehensive T0;
and the recording data alignment correction module 103 is configured to correct a fault starting time point on the other side based on the obtained comprehensive T0 point on any one of the two sides of the fault line.
In the embodiment of the present invention, the data acquisition module 101 sequentially takes 3 sampling points from the ith sampling point to the right along the T axis at intervals from the ith sampling point on the sampling current instantaneous value curve, and the corresponding sampling current instantaneous values are respectively marked as V1, V2, and V3; then, the first processing module 102 calculates the catastrophe points of the A, B, C three-phase current instantaneous value curves Ia, ib and Ic according to the current instantaneous values of the obtained sampling points, and calculates the catastrophe points of the zero-phase self-generated current instantaneous value curve I0 to obtain a comprehensive T0; finally, the recording data alignment correction module 103 corrects the failure start time point on the other side based on the obtained comprehensive T0 point on any one of the two sides of the failure line.
According to the automatic alignment and correction system based on the wave recording data, provided by the embodiment of the invention, an uncertain area of the starting time of the fault is avoided, instead, the adjacent peak value is directly positioned, the peak value belongs to a determined and unique signal, the reliability is high, the influence of the sampling rate difference of the wave recording data on two sides is avoided, the fault tolerance rate of the fault judgment at the starting time point of the fault is high, and the system is easy to popularize and apply.
The system provided by the embodiment of the present invention is used for executing the above method embodiments, and for details of the process and the details, reference is made to the above embodiments, which are not described herein again.
Fig. 5 is a schematic structural diagram of an electronic device according to an embodiment of the present invention, and referring to fig. 5, the electronic device may include: a processor (processor) 201, a communication Interface (communication Interface) 202, a memory (memory) 203 and a communication bus 204, wherein the processor 201, the communication Interface 202 and the memory 203 complete communication with each other through the communication bus 204. The processor 201 may call logic instructions in the memory 203 to perform the following method:
step (1), judging a mutation point T0 by a single-phase sampling instantaneous value curve:
(1.1) on a sampling current instantaneous value curve, sampling points are sequentially taken from the ith sampling point to the right along the T axis at intervals of T from the ith sampling point, and 3 sampling points are taken in total and are respectively marked as Pi and Pi + T, pi +2T; namely, the ith sampling point is taken as the first point; the corresponding sampling current instantaneous values are respectively marked as V1, V2 and V3;
(1.2) calculating the instantaneous value of the current at three sampled points by the following steps
i) The difference value of V1 and V2 is recorded as det1, then det1= | V1-V2|;
ii) the difference between V2 and V3 is denoted as det2, then det2= | V2-V3|;
iii) The difference value of det1 and det12 is marked as detV, and then detV = | det1-det2|;
(1.3) if detV <0.1, the following judgment is performed:
i) If i >15T, T0= ∞, abandoning T0 of the local side to wait for the correction of the contralateral side in the step (2);
ii) if not, taking the next point from i to the right along the t axis, and repeating the steps (1.1) - (1.3);
(1.4) otherwise, if detV >0.1, T0= i +2T;
step (2), three-phase current and zero-phase current catastrophe point comprehensive T0:
adopting the method in the step (1) to sequentially calculate the mutation points of the A, B, C three-phase current instantaneous value curves Ia, ib and Ic, simultaneously calculating the mutation points of the zero-phase self-generated current instantaneous value curve I0, and respectively marking 4 calculated results as T0_ Ia, T0_ Ib, T0_ Ic and T0_ I0;
wherein I0= (Ia + Ib + Ic);
the sum T0 takes the minimum of 4 results, i.e.
Integrated T0= min { T0_ Ia, T0_ Ib, T0_ Ic, T0_ I0};
step (3), setting the two sides of the fault line as an M side and an N side respectively, and then correcting the starting time point of the fault at the other side by taking the comprehensive T0 point at the M side obtained in the step (2) as a reference or taking the comprehensive T0 point at the N side as a reference; here, the M-side integrated T0 point is denoted as MT0, and the N-side integrated T0 point is denoted as NT0.
In addition, the logic instructions in the memory 203 may be implemented in the form of software functional units and stored in a computer readable storage medium when the logic instructions are sold or used as independent products. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
In another aspect, an embodiment of the present invention further provides a non-transitory computer-readable storage medium, on which a computer program is stored, where the computer program is implemented to perform the automatic alignment and correction method based on recorded wave data provided in the foregoing embodiments, for example, the method includes:
step (1), judging a mutation point T0 by a single-phase sampling instantaneous value curve:
(1.1) on a sampling current instantaneous value curve, sampling points are sequentially taken from the ith sampling point to the right along the T axis at intervals of T from the ith sampling point, and 3 sampling points are taken in total and are respectively marked as Pi and Pi + T, pi +2T; namely, the ith sampling point is taken as the first point; the corresponding sampling current instantaneous values are respectively marked as V1, V2 and V3;
(1.2) calculating the instantaneous value of the current at three sampled points by the following steps
i) The difference value of V1 and V2 is recorded as det1, then det1= | V1-V2|;
ii) the difference between V2 and V3 is recorded as det2, then det2= | V2-V3|;
iii) The difference value of det1 and det12 is marked as detV, and then detV = | det1-det2|;
(1.3) if detV <0.1, the following judgment is performed:
i) If i is greater than 15T, T0= ∞, abandoning the T0 of the local side to wait for the correction of the side in the step (2);
ii) if not, taking the next point from i to the right along the t axis, and repeating the steps (1.1) - (1.3);
(1.4) otherwise, if detV >0.1, T0= i +2T;
step (2), three-phase current and zero-phase current catastrophe point comprehensive T0:
calculating the mutation points of A, B, C three-phase current instantaneous value curves Ia, ib and Ic in sequence by adopting the method in the step (1), calculating the mutation points of a zero-phase self-generated current instantaneous value curve I0 at the same time, and respectively marking 4 calculated results as T0_ Ia, T0_ Ib, T0_ Ic and T0_ I0;
wherein I0= (Ia + Ib + Ic);
the sum T0 takes the minimum of 4 results, i.e.
Comprehensive T0= min { T0_ Ia, T0_ Ib, T0_ Ic, T0_ I0};
step (3), setting the two sides of the fault line as an M side and an N side respectively, and then correcting the starting time point of the fault at the other side by taking the comprehensive T0 point at the M side obtained in the step (2) as a reference or taking the comprehensive T0 point at the N side as a reference; here, the M-side integrated T0 point is denoted as MT0, and the N-side integrated T0 point is denoted as NT0.
The above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware. With this understanding in mind, the above-described technical solutions may be embodied in the form of a software product, which can be stored in a computer-readable storage medium such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods described in the embodiments or some parts of the embodiments.
Examples of the applications
As shown in fig. 1, the substations on both sides of the transmission line collect voltage and current signals in real time, convert the signals into secondary current and secondary voltage sampling instantaneous values through CT and PT, and store the sampled values by the in-station wave recorder and the protection device to form static wave recording data. Taking an a-phase fault as an example, a typical fault voltage, fault current and zero-phase current are shown in fig. 1.
(1) Single-phase sampling instantaneous value curve judgment mutation point T0
The method comprises the following steps:
a) As shown in fig. 2, on the sampling current instantaneous value curve, starting from the ith sampling point (the initial value of i is 1, i.e. the first sampling point), 3 sampling points are taken along the T axis in sequence at intervals of T, and are respectively marked as Pi, pi + T, pi +2T, and the corresponding sampling current instantaneous values are respectively marked as V1, V2 and V3;
b) The current instantaneous value of the sampled three points is calculated by the following steps
i) The difference between V1 and V2 is denoted as det1, then det1= | V1-V2 =
ii) the difference between V2 and V3 is denoted as det2, then det2= | V2-V3 =
iii) When the difference between det1 and det12 is detV, detV = | det1-det2
c) If detV <0.1, the following decision is made:
i) If i >15T, T0= ∞ (infinity), abandoning T0 of the local side to wait for correction of the side in the step (2);
ii) otherwise, i takes the next point to the right along the t axis, and repeats the steps a), b) and c);
d) Otherwise, if detV >0.1, T0= i +2T.
(2) Three-phase current and zero-phase current catastrophe point comprehensive T0
The method comprises the following steps:
a) Self-generating zero phase current. The method is formed by adding A, B, C three-phase current instantaneous values point by point, and the formula is as follows:
I0=(Ia+Ib+Ic)；
the resulting I0 curve is shown in fig. 3.
b) And (2) according to the method for judging the mutation point T0 by the single-phase sampling instantaneous value curve in the step (1), calculating the mutation points of Ia, ib, ic and I0 in sequence, and respectively marking 4 calculated results as T0_ Ia, T0_ Ib, T0_ Ic and T0_ I0.
c) Taking the minimum value of the 4 results in step b) as the comprehensive T0, i.e.
T0=min{T0_Ia，T0_Ib，T0_Ic，T0_I0}
(3) The recording data at two sides respectively judge the starting time point T0 of the fault
Let both sides be M side and N side respectively, the term definition of relevant trouble record wave curve is as table 1.
TABLE 1 two-sided Fault recording information expression symbol definition
The method comprises the following steps:
with MT0 as a reference, NT0 is corrected, and the corrected value of NT0 is: NT0'.
The method comprises the following steps:
1) Curve correcting method based on I0 sampling instantaneous value
NT0’=NT0+(NdetT_I0-MdetT_I0)
2) Ia sampling based instantaneous value curve correction method
NT0’=NT0-(NdetT_Ia-MdetT_Ia)
3) Ua sampling-based instantaneous value curve correction method
NT0’=NT0-(NdetT_Ua-MdetT_Ua)
The logic of the using condition of the correction method 1) 2) 3) is as follows:
a) At the time of T0, if the secondary fundamental wave effective value of I0 is greater than 1/2 of the Ia fundamental wave effective value, using a curve correction method based on I0 sampling instantaneous values;
b) Otherwise, if the secondary fundamental wave effective value of Ia is greater than 0.1A, using a curve correction method based on Ia sampling instantaneous values;
c) Otherwise, using a curve correction method based on the Ua sampling instantaneous value;
d) The fundamental wave effective values in the steps a) and b), the data window used by I0 is the length of T on the right side of T0, and the data window used by Ia is the length of T on the left side of T0.
The second method comprises the following steps:
correcting MT0 by taking NT0 as a reference, and setting the corrected value of MT0 as follows: MT0'.
The method comprises the following steps:
1) Curve correcting method based on I0 sampling instantaneous value
MT0’=MT0+(MdetT_I0-NdetT_I0)
2) Ia sampling based instantaneous value curve correction method
MT0’=MT0-(MdetT_Ia-NdetT_Ia)
3) Ua sampling-based instantaneous value curve correction method
MT0’=MT0-(MdetT_Ua-NdetT_Ua)
The logic of the using condition of the correction method 1) 2) 3) is as follows:
a) At the time of T0, if the secondary fundamental wave effective value of I0 is greater than 1/2 of the Ia fundamental wave effective value, using a curve correction method based on I0 sampling instantaneous values;
b) Otherwise, if the secondary fundamental wave effective value of Ia is greater than 0.1A, using a curve correction method based on Ia sampling instantaneous values;
c) Otherwise, using a curve correction method based on Ua sampling instantaneous value;
d) The fundamental wave effective values in the steps a) and b), the data window used by I0 is the length of T on the right side of T0, and the data window used by Ia is the length of T on the left side of T0.
Specifically, the following description is provided: the method is that the fault phase is assumed to be A, and the fault voltage and the fault current are taken to be Ua and Ia, and similarly, if the fault phase is B, the fault voltage and the fault current are taken to be Ub and Ib; and if the fault phase is C, the fault voltage and the fault current are Uc and Ic. Are all within the scope of protection of the patented process.
The foregoing shows and describes the general principles, essential features, and advantages of the invention. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, which are given by way of illustration of the principles of the present invention, but that various changes and modifications may be made without departing from the spirit and scope of the invention, and such changes and modifications are within the scope of the invention as claimed. The scope of the invention is defined by the appended claims and equivalents thereof.
Claims (6)
1. An automatic alignment and correction method based on recorded wave data is characterized by comprising the following steps:
step (1), judging a mutation point T0 by a single-phase sampling instantaneous value curve:
(1.1) on a sampling current instantaneous value curve, sampling points are sequentially taken from the ith sampling point to the right along the T axis at intervals of T from the ith sampling point, and 3 sampling points are taken in total and are respectively marked as Pi and Pi + T, pi +2T; namely, the ith sampling point is taken as the first point; the corresponding sampling current instantaneous values are respectively marked as V1, V2 and V3;
(1.2) calculating the instantaneous value of the current at three sampled points by the following steps
i) The difference value of V1 and V2 is recorded as det1, then det1= | V1-V2|;
ii) the difference between V2 and V3 is denoted as det2, then det2= | V2-V3|;
iii) The difference value of det1 and det2 is recorded as detV, then detV = | det1-det2|;
(1.3) if detV <0.1, the following judgment is performed:
i) If i >15T, T0= ∞, abandoning T0 of the local side to wait for the correction of the contralateral side in the step (2);
ii) if not, taking the next point from i to the right along the t axis, and repeating the steps (1.1) - (1.3);
(1.4) otherwise, if detV >0.1, T0= i +2T;
step (2), three-phase current and zero-phase current catastrophe point comprehensive T0:
calculating the mutation points of A, B, C three-phase current instantaneous value curves Ia, ib and Ic in sequence by adopting the method in the step (1), calculating the mutation points of a zero-phase self-generated current instantaneous value curve I0 at the same time, and respectively marking 4 calculated results as T0_ Ia, T0_ Ib, T0_ Ic and T0_ I0;
wherein I0= (Ia + Ib + Ic);
the sum of T0 takes the minimum of 4 results, i.e.
Comprehensive T0= min { T0_ Ia, T0_ Ib, T0_ Ic, T0_ I0};
step (3), setting the two sides of the fault line as an M side and an N side respectively, and then correcting the starting time point of the fault at the other side by taking the comprehensive T0 point at the M side obtained in the step (2) as a reference or taking the comprehensive T0 point at the N side as a reference; wherein, the comprehensive T0 point at the M side is marked as MT0, and the comprehensive T0 point at the N side is marked as NT0;
in the step (3), correcting the N-side fault starting time point NT0 by taking the M-side fault starting time point MT0 as a reference, setting a corrected value NT0' of NT0, wherein the fault phase is X, and the fault voltage and the fault current are Ux and Ix; x is A, B or C; the method comprises the following steps:
(3.1.1) curve correction method based on I0 sampling instantaneous value:
NT0’=NT0+(NdetT_I0-MdetT_I0)；
(3.1.2) curve correction method based on Ix sampling instantaneous value:
NT0’=NT0-(NdetT_Ix-MdetT_Ix)；
(3.1.3) instantaneous value curve correction method based on Ux sampling:
NT0’=NT0-(NdetT_Ux-MdetT_Ux)；
wherein MdetT _ Ux represents the time difference of the first peak value on the left side of T0 between M side T0 and Ux;
NdetT _ Ux represents the time difference between the first peak value on the left side of T0 of the N side T0 and Ux;
MdetT _ Ix represents the time difference of the first peak value on the left side of T0 between M side T0 and Ix;
NdetT _ Ix represents the time difference of the first peak value of the N side T0 and Ix on the left side of T0;
MdetT _ I0 represents the time difference of the first peak values of the M side T0 and I0 on the right side of T0;
NdetT _ I0 represents the time difference of the first peak value of the N side T0 and the right side I0 of T0;
the logic of the use conditions of the correction methods from (3.1.1) to (3.1.3) is as follows:
a) At the time of T0, if the effective value of the secondary fundamental wave of I0 is greater than 1/2 of the effective value of the Ix fundamental wave, using a curve correction method based on I0 sampling instantaneous values;
b) Otherwise, if the secondary fundamental wave effective value of Ix is greater than 0.1A, using a curve correction method based on Ix sampling instantaneous values;
c) Otherwise, using a curve correction method based on Ux sampling instantaneous values;
a) And b), when the secondary fundamental wave effective value and the Ix fundamental wave effective value of the I0 are acquired, the data window used by the I0 is the length of the T on the right side of the T0, and the data window used by the Ix is the length of the T on the left side of the T0.
2. The recording data-based automatic alignment and correction method of claim 1, wherein: in the step (3), correcting the M-side fault starting time point MT0 by taking the N-side fault starting time point NT0 as a reference, setting the corrected value of MT0 as MT0', the fault phase as X, and the fault voltage and the fault current as Ux and Ix; x is A, B or C; the method comprises the following steps:
(3.2.1) curve correction method based on I0 sampling instantaneous value:
MT0’=MT0+(MdetT_I0-NdetT_I0)；
(3.2.2) curve correction method based on Ix sampling instantaneous value:
MT0’=MT0-(MdetT_Ix-NdetT_Ix)；
(3.2.3) instantaneous value curve correction method based on Ux sampling:
MT0’=MT0-(MdetT_Ux-NdetT_Ux)；
wherein MdetT _ Ux represents the time difference of the first peak value on the left side of T0 between M side T0 and Ux;
NdetT _ Ux represents the time difference between the first peak value on the left side of T0 of the N side T0 and Ux;
MdetT _ Ix represents the time difference of the first peak values of the M side T0 and Ix on the left side of T0;
NdetT _ Ix represents the time difference of the first peak value of the N side T0 and Ix on the left side of T0;
MdetT _ I0 represents the time difference of the first peak value of the M side T0 and I0 on the right side of T0;
NdetT _ I0 represents the time difference of the first peak value of the N side T0 and the right side I0 of T0;
the logic of the use conditions of the correction methods from (3.2.1) to (3.2.3) is as follows:
(a) At the time of T0, if the secondary fundamental wave effective value of I0 is greater than 1/2 of the Ix fundamental wave effective value, using a curve correction method based on I0 sampling instantaneous values;
(b) Otherwise, if the effective value of the secondary fundamental wave of Ix is greater than 0.1A, using a curve correction method based on Ix sampling instantaneous values;
(c) Otherwise, curve correction based on Ux sampling transients is used.
3. The recording data-based automatic alignment and correction method of claim 2, wherein: (a) And (b) when the secondary fundamental wave effective value and the Ix fundamental wave effective value of the I0 are acquired, the data window used by the I0 is the length of the T on the right side of the T0, and the data window used by the Ix is the length of the T on the left side of the T0.
4. An automatic alignment and correction system based on recording data, which adopts the automatic alignment and correction method based on recording data of claim 1, and is characterized by comprising:
the data acquisition module is used for sequentially taking sampling points from the ith sampling point to the right along the T axis by taking T as an interval from the ith sampling point on a sampling current instantaneous value curve, and taking 3 sampling points in total, wherein the corresponding sampling current instantaneous values are respectively marked as V1, V2 and V3;
the first processing module is used for calculating the catastrophe points of A, B, C three-phase current instantaneous value curves Ia, ib and Ic according to the current instantaneous values of the obtained sampling points, and calculating the catastrophe points of a zero-phase self-generated current instantaneous value curve I0 to obtain comprehensive T0;
and the recording data alignment correction module is used for correcting the fault starting time point on the other side according to the obtained comprehensive T0 point on any one side of the two sides of the fault line as a reference.
5. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor executes the program to implement the steps of the method for automatic alignment and correction based on recorded wave data according to any one of claims 1 to 3.
6. A non-transitory computer-readable storage medium, on which a computer program is stored, wherein the computer program, when executed by a processor, implements the steps of the method for automatic alignment and rectification based on oscillometric data according to any of claims 1 to 3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010592547.0A CN111812454B (en) | 2020-06-24 | 2020-06-24 | Automatic alignment and correction method and system based on wave recording data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010592547.0A CN111812454B (en) | 2020-06-24 | 2020-06-24 | Automatic alignment and correction method and system based on wave recording data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111812454A CN111812454A (en) | 2020-10-23 |
| CN111812454B true CN111812454B (en) | 2023-01-24 |
Family
ID=72856461
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010592547.0A Active CN111812454B (en) | 2020-06-24 | 2020-06-24 | Automatic alignment and correction method and system based on wave recording data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111812454B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111679151B (en) * | 2020-05-27 | 2022-06-21 | 武汉华电顺承科技有限公司 | Synchronous alignment point vector calculation method for double-end distance measurement of wave recording data |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6255570A (en) * | 1985-09-04 | 1987-03-11 | Showa Electric Wire & Cable Co Ltd | Location of fault point for cable line |
| CN105116246A (en) * | 2015-07-27 | 2015-12-02 | 云南电力调度控制中心 | Power grid operation protection real-time online analysis method |
| CN106597215A (en) * | 2016-11-17 | 2017-04-26 | 云南电网有限责任公司电力科学研究院 | Fault location method and system based on transient travelling waves within short time window |
| CN107222280A (en) * | 2017-05-24 | 2017-09-29 | 云南电力试验研究院(集团)有限公司 | A kind of method that time synchronized is realized to different time scales fault recording system |
| CN108107321A (en) * | 2017-12-14 | 2018-06-01 | 科大智能电气技术有限公司 | A kind of electric power system fault waveform comparison method |
| CN109738762A (en) * | 2019-02-28 | 2019-05-10 | 国网湖北省电力有限公司黄石供电公司 | A kind of middle voltage distribution networks cable short trouble localization method based on both-end traveling wave method |
| WO2019097312A1 (en) * | 2017-11-17 | 2019-05-23 | Abb Schweiz Ag | Parameter free traveling wave based fault location for power transmission lines |
| CN109861187A (en) * | 2019-03-18 | 2019-06-07 | 湖南大学 | A kind of DC line protection method in DC distribution net |
| CN110780251A (en) * | 2019-10-31 | 2020-02-11 | 南京能云电力科技有限公司 | Power waveform comparison system and method based on point-by-point calculation of DELTA |
-
2020
- 2020-06-24 CN CN202010592547.0A patent/CN111812454B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6255570A (en) * | 1985-09-04 | 1987-03-11 | Showa Electric Wire & Cable Co Ltd | Location of fault point for cable line |
| CN105116246A (en) * | 2015-07-27 | 2015-12-02 | 云南电力调度控制中心 | Power grid operation protection real-time online analysis method |
| CN106597215A (en) * | 2016-11-17 | 2017-04-26 | 云南电网有限责任公司电力科学研究院 | Fault location method and system based on transient travelling waves within short time window |
| CN107222280A (en) * | 2017-05-24 | 2017-09-29 | 云南电力试验研究院(集团)有限公司 | A kind of method that time synchronized is realized to different time scales fault recording system |
| WO2019097312A1 (en) * | 2017-11-17 | 2019-05-23 | Abb Schweiz Ag | Parameter free traveling wave based fault location for power transmission lines |
| CN108107321A (en) * | 2017-12-14 | 2018-06-01 | 科大智能电气技术有限公司 | A kind of electric power system fault waveform comparison method |
| CN109738762A (en) * | 2019-02-28 | 2019-05-10 | 国网湖北省电力有限公司黄石供电公司 | A kind of middle voltage distribution networks cable short trouble localization method based on both-end traveling wave method |
| CN109861187A (en) * | 2019-03-18 | 2019-06-07 | 湖南大学 | A kind of DC line protection method in DC distribution net |
| CN110780251A (en) * | 2019-10-31 | 2020-02-11 | 南京能云电力科技有限公司 | Power waveform comparison system and method based on point-by-point calculation of DELTA |
Non-Patent Citations (3)
| Title |
|---|
| 《Traveling Wave-Based Solutions for Transmission Line Two-Terminal Data Time Synchronization》;Felipe V. Lopes等;《IEEE Transactions on Power Delivery》;20181231;全文 * |
| 《双端测距以及同步算法的实现》;吴瑞春等;《继电器》;20080301;第36卷(第5期);7-9 * |
| 《多变电站故障录波数据同步自动化方法与实现》;王海港等;《电力***保护与控制》;20150101;第43卷(第1期);102-107 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111812454A (en) | 2020-10-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108535543B (en) | Transformer area phase identification method based on synchronous sampling of concentrator and ammeter | |
| JP6636608B2 (en) | Method and apparatus for detecting voltage in a supply network | |
| Liang et al. | Fault location method in power network by applying accurate information of arrival time differences of modal traveling waves | |
| CN107222280A (en) | A kind of method that time synchronized is realized to different time scales fault recording system | |
| CN110927521A (en) | Single-ended traveling wave fault positioning method and device | |
| US20180191241A1 (en) | A Method for Correcting Effect of Saturation in Current Transformer and an Intelligent Electronic Device Therefor | |
| WO2003044547A1 (en) | Fault location using measurements from two ends of a line | |
| CN111812454B (en) | Automatic alignment and correction method and system based on wave recording data | |
| JPWO2008126240A1 (en) | Synchronous phasor measuring device and phase angle difference measuring device between buses using the same | |
| CN109767133B (en) | Electric quantity compensation method and device based on load recorded data | |
| EP3016339A1 (en) | Cycle slip detection method and device, and receiver | |
| Kuwalek | AM modulation signal estimation allowing further research on sources of voltage fluctuations | |
| CN111812455B (en) | Power transmission line double-end fault distance measurement method and system | |
| CN107179476B (en) | Distribution network fault distance measurement method | |
| CN109375048B (en) | Power transmission line parameter identification method and device based on fault recording data | |
| da Costa et al. | Real-time evaluation of impedance-based fault location algorithms | |
| Džafić et al. | Enhancing Prony’s method for fault location | |
| Torkaman et al. | Rearrangement method of reducing fault location error in tied uncompleted parallel lines | |
| JP2011072163A (en) | Ground fault direction finding apparatus | |
| CN113640615B (en) | Small-current ground fault line selection method based on evidence combination stable transient state information | |
| Redondo et al. | A strategy for improving the accuracy of flicker emission measurement from wind turbines | |
| CN111929630B (en) | Method and device for detecting saturation of current transformer | |
| Jovicic et al. | Linear state estimation considering refresh rates of RTU and PMU measurements | |
| JPH0345344B2 (en) | ||
| JPS58219463A (en) | Fault point location system for 4 terminal transmission line |
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 |