CN112506867B - Recording comparison analysis method based on time deviation management - Google Patents
Recording comparison analysis method based on time deviation management Download PDFInfo
- Publication number
- CN112506867B CN112506867B CN202011296858.9A CN202011296858A CN112506867B CN 112506867 B CN112506867 B CN 112506867B CN 202011296858 A CN202011296858 A CN 202011296858A CN 112506867 B CN112506867 B CN 112506867B
- Authority
- CN
- China
- Prior art keywords
- comparison
- time
- wave recording
- ied
- sampling point
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 15
- 238000011214 deviation management Methods 0.000 title claims abstract description 14
- 238000007726 management method Methods 0.000 claims abstract description 55
- 230000035772 mutation Effects 0.000 claims abstract description 46
- 238000005070 sampling Methods 0.000 claims description 197
- 238000001914 filtration Methods 0.000 claims description 14
- 238000000034 method Methods 0.000 claims description 14
- 238000012937 correction Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 6
- 230000000630 rising effect Effects 0.000 claims description 6
- 230000005856 abnormality Effects 0.000 claims description 4
- 238000011084 recovery Methods 0.000 claims description 4
- 230000000052 comparative effect Effects 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 2
- 230000001052 transient effect Effects 0.000 claims 1
- 230000002159 abnormal effect Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/10—File systems; File servers
- G06F16/16—File or folder operations, e.g. details of user interfaces specifically adapted to file systems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/10—File systems; File servers
- G06F16/11—File system administration, e.g. details of archiving or snapshots
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F40/00—Handling natural language data
- G06F40/10—Text processing
- G06F40/194—Calculation of difference between files
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/20—Administration of product repair or maintenance
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Quality & Reliability (AREA)
- Computational Linguistics (AREA)
- Health & Medical Sciences (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- General Health & Medical Sciences (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Artificial Intelligence (AREA)
- Human Computer Interaction (AREA)
- Recording Measured Values (AREA)
Abstract
The invention discloses a recording comparison analysis method based on time deviation management, which comprises the steps of calculating the time deviation between a protection device and an acquisition unit and a management unit, matching a starting signal, searching a corresponding recording file, calculating a starting mutation point and a resetting mutation point, carrying out waveform alignment by using a zero crossing point, and carrying out analog quantity and switching quantity comparison alarm.
Description
Technical Field
The invention relates to the technical field of operation and maintenance of intelligent substations, in particular to a recording comparison analysis method based on time deviation management.
Technical Field
The sampling loop plays a crucial role for the proper operation of the protection device. The sampling loop is difficult to find due to the abnormality such as single-point grounding. Although some comparison methods based on the wave recording file exist at present, the problems exist, and the problems of abnormal time setting of the protection device, faults outside a double-circuit line region, high-order harmonic errors and the like are not solved. Therefore, in order to improve the abnormal positioning capability of the sampling loop of the transformer substation, a recording comparison analysis method based on time deviation management is researched, and the method has important significance for safe and reliable operation of the transformer substation.
Disclosure of Invention
Aiming at the defects in the prior art, the invention provides a recording comparison analysis method based on time deviation management, which can realize fault diagnosis and abnormal positioning of an intelligent substation sampling loop and improve the safety and reliability of the substation.
In order to achieve the purpose, the invention adopts the following technical scheme:
a recording comparison analysis method based on time deviation management comprises the following steps:
step 2, calculating the time deviation between the IED and the management unit;
step 3, the management unit acquires a comparison IED associated with each IED according to the channel configuration, matches the IED with a starting signal of the associated comparison IED, and acquires a corresponding recording file as a recording file to be compared;
step 4, correcting the recording time of the name of the wave recording file to be compared according to the time deviation of each IED, and grouping the IED and the wave recording files with overlapped wave recording file recording time sections of the compared IED to obtain wave recording files of a group;
step 5, respectively finding a starting mutation point and a resetting mutation point from wave recording files of the team in the associated analog quantity channel;
step 6, the current channel with the maximum effective current value in the wave recording files of the team is a basic channel for calculating a zero crossing point, and the accurate time deviation t _ diff between the IED and the comparison IED is calculated according to the starting break points of the basic channels of the wave recording files of the team;
7, finding out that all grouped start mutation points and all grouped recovery mutation points in the wave recording files of the group team have overlapped sections in time;
step 8, carrying out digital filtering on sampling point data of the wave recording file of the acquisition unit according to the filtering processing control word of the management unit;
step 9, selecting a wave recording file with a low sampling rate from the wave recording files of the team in the management unit as a reference wave recording file; the other wave recording file in the wave recording files of the team is a non-reference wave recording file, each sampling point in the reference wave recording file and each sampling point in the non-reference wave recording file are corrected according to the accurate time deviation t _ diff, the sampling point of an analog channel in the reference wave recording file is a reference sampling point, a corresponding comparison sampling point is found in the non-reference wave recording file, the reference sampling point is compared with the comparison sampling point, and the time of a switching value wave recording channel participating in comparison is corrected and aligned according to the accurate time deviation t _ diff; and comparing the time of the variable sites of the associated switching value recording channels.
The time offset in step 2 as described above is obtained by the following steps;
the management unit acquires and stores the current time of the protection device and the current time of the acquisition unit in an NTP mode, the time for the management unit to send an NTP request is defined as T1, the time for the protection device or the acquisition unit to receive the NTP request is defined as T2, the reply time of the protection device or the acquisition unit is defined as T3, and the reply time received by the management unit is defined as T4, so that the network delay d and the time deviation delta T between the protection device or the acquisition unit and the management unit are calculated according to the following formula:
d＝(T2-T1)+(T4-T3)；
△t＝[(T2-T1)+(T3-T4)]/2。
step 3 as described above comprises the steps of:
step 3.1, a starting signal buffer queue is built for each IED in the management unit, and starting signals of the IEDs are buffered in the starting signal buffer queue;
step 3.2, the management unit acquires a comparison IED associated with each IED according to the channel configuration;
step 3.3, checking whether the IED start signal buffer queue is empty every time the time length defined by the start signal time difference threshold, and if not, matching the start signal in the IED start signal buffer queue to the comparison IED start signal buffer queue, specifically:
if the difference value of the starting signal time of the starting signals in the starting signal cache queue of the IED and the starting signal time of the starting signal in the starting signal cache queue of the comparison IED is smaller than the configured starting signal time difference threshold value, matching the starting signals in the starting signal cache queue of the IED and the starting signal time in the starting signal cache queue of the comparison IED, wherein the starting signal time is the starting signal time after the time deviation correction obtained in the step 2;
step 3.4, the unmatched starting signals in the starting signal buffer queue are cleared if the unmatched starting signals are not matched after three times of service automatic starting cycles;
step 3.5, finding a corresponding recording file from the successfully called recording file according to the matched IED starting signal and the compared IED starting signal, and if the starting signal time and the recording file name recording time are less than 50ms, judging that the recording file corresponding to the starting signal is found;
step 3.6, when the IED and the comparison IED both find the wave recording file corresponding to the matched starting signal, taking the two wave recording files as wave recording files to be compared;
and 3.7, if the wave recording file corresponding to the starting signal is not found in the period of the current service automatic starting, the wave recording file is found again in the period of the next service automatic starting, and the wave recording file corresponding to the matched starting signal of the IED and the comparison IED is discarded if the wave recording file corresponding to the starting signal is not found in the period of the three continuous service automatic starting.
Step 4 as described above comprises the steps of:
step 4.1, the recording file name recording time of the recording files to be compared of the local IED and the comparison IED obtained in the step 3 is corrected according to the time deviation of the respective IED in the management unit, and the recording files to be compared of the local IED and the comparison IED are respectively sequenced according to the corrected recording file name recording time;
4.2, analyzing respective wave recording files to be compared of the IED and the comparison IED and acquiring a corresponding recording time section;
4.3, grouping the IED and the wave recording files with the overlapped wave recording file recording time sections of the comparison IED to obtain wave recording files of the group;
and 4.4, discarding the recording file of which the overlapping area is not found.
Step 5 as described above comprises the steps of:
step 5.1, selecting analog quantity channels with the association relationship established in the channel configuration of the wave recording files grouped in the step 4 in the management unit;
step 5.2, in the wave recording file of the team, in the associated analog quantity channel, removing the wave recording section with the sampling rate lower than the sampling rate invalid threshold value in any analog quantity channel;
step 5.3, in the wave recording file of the team, removing wave recording sections with discontinuous occurrence time in any analog quantity channel from the related analog quantity channel;
step 5.4, in the wave recording file of the team, removing sampling points with instantaneous values smaller than 0.001 times of rated values in any analog quantity channel in the related analog quantity channel;
and 5.5, in the wave recording file of the team, in each remaining wave recording section in the associated analog quantity channel, starting from the first sampling point of the third cycle, when the sudden change instantaneous values of the continuous three sampling points are greater than a sudden change starting threshold value, the first sampling point of the continuous three sampling points is a starting sudden change point, and when the sudden change instantaneous values of the sampling points of the continuous cycle after the starting of the sudden change point are less than the sudden change starting threshold value, the first sampling point of the continuous cycle is a repeated sudden change point.
The instantaneous value of the mutation in step 5.5 as described above is calculated as follows:
note ikFor the instantaneous sample value corresponding to the kth sample point, ik-NFor the instantaneous value, i, of the sampling point one sampling period before the kth sampling pointk-2NThe instantaneous value delta i of the sudden change of the kth sampling point is the instantaneous sampling value of the sampling point two sampling periods before the kth sampling pointkComprises the following steps:
Δik＝||ik-ik-N|-|ik-N-ik-2N||。
step 6 as described above comprises the steps of:
step 6.1, selecting current channels with established association relations in the channel configuration in the wave recording files of the team in the management unit, and finding the current channel with the maximum current effective value in the wave recording files of the team as a basic channel for calculating a zero crossing point;
step 6.2, starting mutation points of respective basic channels of the wave recording files of the team are taken, the time difference of the starting mutation points of the respective basic channels of the wave recording files of the team after time deviation correction firstly meets a pair of starting mutation points which are smaller than 5ms, the starting mutation points are respectively defined as a first starting mutation point A and a first starting mutation point B, and the time of the first starting mutation point A and the time of the first starting mutation point B after time deviation correction are respectively t1 and t 2;
step 6.3, finding a first rising trend zero-crossing point after the basic channels of the wave recording files of the team are added by 60ms by t 1;
step 6.4, if the zero-crossing point is not at the sampling point, calculating the time of the first rising trend zero-crossing point of the basic channels of the local IED and the comparative IED found in the step 6.3 by using an equal ratio difference method to be t3 and t4 respectively;
if (t3-t1) - (t4-t1) > < 10000us, t3 ═ t3-20000us, if (t4-t1) - (t3-t1) > < 10000us, t4 ═ t4-20000 us;
the precise time offset t _ diff between the present IED and the comparison IED is calculated as t3-t 4.
Step 9 as described above comprises the steps of:
step 9.1, selecting a wave recording file with a low sampling rate from the wave recording files of the team as a reference wave recording file in the management unit; the other wave recording file in the wave recording files of the team is a non-reference wave recording file, and each sampling point in the reference wave recording file and the non-reference wave recording file is corrected according to the accurate time deviation t _ diff obtained in the step 6.7;
step 9.2, taking the sampling point of the analog channel in the reference wave recording file as a reference sampling point, finding the sampling point by using the time of the reference sampling point of the analog channel of the reference wave recording file to the analog channel in the non-reference wave recording file, if the sampling point corresponding to the time is found, marking as a comparison sampling point, directly comparing the reference sampling point with the comparison sampling point, if the sampling point corresponding to the time cannot be found, calculating a comparison sampling point by using an interpolation algorithm, and comparing the reference sampling point with the comparison sampling point;
step 9.3, comparing the reference sampling point with the comparison sampling point comprises: comparing the reference sampling point with the instantaneous value of the comparison sampling point, comparing the effective amplitude of the next cycle from the reference sampling point with the effective amplitude of the next cycle from the comparison sampling point, and comparing the phase angle corresponding to the reference sampling point and the comparison sampling point;
9.4, comparing the difference value between the reference sampling point and the comparison sampling point instantaneous value obtained in the step 9.3 with an analog quantity instantaneous value comparison difference threshold value, and generating comparison abnormity when the amplitude difference between the effective amplitude of the next cycle from the reference sampling point and the effective amplitude of the next cycle from the comparison sampling point is greater than the analog quantity effective value amplitude comparison difference threshold value;
or the difference value between the reference sampling point and the comparison sampling point instantaneous value obtained in the step 9.3 is greater than the analog quantity instantaneous value comparison difference threshold, and the comparison abnormity is generated when the phase angle difference between the reference sampling point and the comparison sampling point is greater than the analog quantity effective value angle comparison difference threshold;
if the comparison abnormality occurs between the three continuous pairs of reference sampling points and the comparison sampling points, judging the channel comparison alarm;
step 9.6, acquiring the switching value wave recording channel for establishing the association relationship in the step 1 from the wave recording files participating in comparison;
9.7, correcting and aligning the time of the switching value recording channel participating in comparison according to the accurate time deviation t _ diff calculated in the step 6;
9.8, searching the time of all variable sites of the switching value recording channels participating in comparison;
and 9.9, comparing the time of the variable sites of the associated switching value recording channels, and giving an alarm when the time difference is greater than the switching value comparison difference threshold set in the step 1.
Compared with the prior art, the invention has the following beneficial effects:
1. the invention provides a recording comparison analysis method based on time deviation management. The accurate comparison of the homologous data of the wave recording file of the intelligent substation can be realized;
2. the invention can realize fault diagnosis and abnormal positioning of the intelligent substation sampling loop and improve the safety and reliability of the substation.
Drawings
FIG. 1 is a schematic diagram of the process steps of the present invention;
FIG. 2 is a schematic diagram of a mutation algorithm of the present invention;
FIG. 3 is a schematic diagram of a zero crossing point time calculation method of the present invention.
Detailed Description
The present invention will be described in further detail with reference to examples for the purpose of facilitating understanding and practice of the invention by those of ordinary skill in the art, and it is to be understood that the present invention has been described in the illustrative embodiments and is not to be construed as limited thereto.
As shown in fig. 1, a recording comparison analysis method based on time deviation management includes the following steps:
the management unit channel configuration needs to acquire a CFG file of the IED, the IED comprises an acquisition unit, a protection device and the like, and basic information of each channel is obtained through analysis. The channel configuration associates analog quantity channels of the same sampling source in the wave recording files of the protection device and the acquisition unit, and associates switching value channels reflecting the same switching value deflection. The protection device establishing a mapping relation with the acquisition unit at least comprises a line protection device, a main transformer protection device and a bus protection device;
when the bus connection form is a double-bus connection, a double-bus double-section connection or a double-bus single-section connection, the management unit needs to associate a voltage channel in a wave recording file of the line protection device and a voltage channel in a wave recording file of the main transformer protection device with position signals of the disconnecting switches at intervals where the line protection device and the main transformer protection device are located, and select the correspondingly connected bus voltage through the connection state of the disconnecting switches and the buses.
Step 2, calculating a time deviation delta t between the IED (including the protection device and the acquisition unit) and the management unit, and specifically comprising the following steps:
the management unit acquires and stores the current time of the protection device and the current time of the acquisition unit in an NTP mode. Defining the time of sending the NTP request by the management unit as T1, the time of receiving the NTP request by the protection device or the acquisition unit as T2, the reply time of the protection device or the acquisition unit as T3, and the time of receiving the reply by the management unit as T4, then the network delay d and the time deviation Deltat between the protection device or the acquisition unit and the management unit are calculated as follows:
d＝(T2-T1)+(T4-T3)；
△t＝[(T2-T1)+(T3-T4)]/2；
when the protection device does not support NTP, MMS information can be received through the management unit, the MMS message time is analyzed to obtain the current time of the protection device or the current time of the acquisition unit, the difference is made between the current time of the protection device and the current time of the acquisition unit and the local time of the management unit, the time deviation between the protection device and the management unit is calculated, and the time deviation between the acquisition device and the management unit is calculated. And meanwhile, the management unit corrects the received event sequence recording (SOE) time of the acquisition unit or the protection device into event sequence recording time under the time of the management unit through time deviation, wherein the event sequence recording time comprises starting signal time.
Step 3, the management unit acquires a comparison IED associated with each IED according to the channel configuration, matches the IED with a starting signal of the associated comparison IED, and acquires a corresponding recording file as a recording file to be compared, and the method comprises the following steps:
step 3.1, a starting signal buffer queue is built for each IED in the management unit, and starting signals of the IEDs are buffered in the starting signal buffer queue;
step 3.2, the management unit acquires a comparison IED associated with each IED according to the channel configuration;
step 3.3, checking whether the IED start signal buffer queue is empty every time length defined by the start signal time difference threshold, and if not, matching the start signal in the IED start signal buffer queue to the comparison IED start signal buffer queue, specifically:
and if the difference value of the starting signal time of the starting signals in the starting signal cache queue of the IED and the starting signal time of the comparison IED is smaller than the configured starting signal time difference threshold value, matching the starting signals in the starting signal cache queue of the IED and the starting signal time of the comparison IED. The starting signal time is the starting signal time after the time deviation correction obtained in the step 2;
step 3.4, the unmatched starting signals in the starting signal buffer queue are cleared if the unmatched starting signals are not matched after three times of service automatic starting cycles;
step 3.5, finding a corresponding recording file from the successfully called recording file according to the matched IED starting signal and the compared IED starting signal, and if the starting signal time and the recording file name recording time are less than 50ms, judging that the recording file corresponding to the starting signal is found;
step 3.6, only if the IED and the comparison IED find the wave recording file corresponding to the matched starting signal, the two wave recording files are used as the wave recording files to be compared;
step 3.7, if the wave recording file corresponding to the starting signal is not found in the period of the current service automatic starting, the wave recording file is found again in the period of the next service automatic starting, and the wave recording file corresponding to the matched starting signal of the IED and the comparison IED is discarded if the wave recording file is not found in the period of the three continuous service automatic starting;
step 4, correcting the recording time of the name of the wave recording file to be compared according to the time deviation of each IED, and grouping the IED and the wave recording files with the recording time sections of the wave recording files of the compared IED overlapped to obtain the wave recording files of the group, wherein the method specifically comprises the following steps:
step 4.1, the recording file name recording time of the recording files to be compared of the local IED and the comparison IED obtained in the step 3 is corrected according to the time deviation of the respective IED in the management unit, and the recording files to be compared of the local IED and the comparison IED are respectively sequenced according to the corrected recording file name recording time;
step 4.2, analyzing respective wave recording files to be compared of the IED and the comparison IED and acquiring corresponding recording time sections;
step 4.3, grouping the IED and the wave recording files with overlapped wave recording file recording time sections of the comparison IED to obtain wave recording files of the group;
and 4.4, discarding the recording file of which the overlapping area is not found.
Step 5, in the associated analog quantity channel, respectively finding a starting mutation point and a resetting mutation point from the wave recording files of the team, specifically comprising the following steps:
step 5.1, selecting analog quantity channels with the association relationship established in the channel configuration of the wave recording files grouped in the step 4 in the management unit;
step 5.2, in the wave recording file of the team, in the associated analog quantity channel, removing the wave recording section with the sampling rate lower than the sampling rate invalid threshold value in any analog quantity channel;
step 5.3, in the wave recording file of the team, removing wave recording sections with discontinuous occurrence time in any analog quantity channel from the related analog quantity channel;
step 5.4, in the wave recording file of the team, removing sampling points with instantaneous values smaller than 0.001 times of rated values in any analog quantity channel in the related analog quantity channel;
and 5.5, in the wave recording files of the team, in each remaining wave recording section in the associated analog quantity channel, starting a sudden change point and a resetting sudden change point from the first sampling point of the third cycle according to the following formulas. As shown in FIG. 2, where Δ ikCalculating the instantaneous value of the sudden change i for the fault component at the k-th sampling pointkFor the instantaneous sample value corresponding to the kth sample point, ik-NFor the instantaneous value of a sample point one sampling period before the kth sample point (one sampling period includes N sample points), ik-2NThe instantaneous sample value of the sample point two sample periods before the k-th sample point. When Δ ik、△ik+1、△ik+2When the sampling points are all larger than the start threshold of the mutation quantity, the k-th sampling point is a start mutation point, and a cycle wave delta i is continuously generated after the start mutation pointk<When the abrupt change quantity starts the threshold value, the k-th sampling point is a resetting abrupt change point corresponding to the starting abrupt change point; particularly, if two continuous sampling points which are less than 0.001-time rated value exist in the eight sampling points behind the current sampling point, the current sampling point is judged to be a reset catastrophe point;
Δik＝||ik-ik-N|-|ik-N-ik-2N||
step 5.6, the starting catastrophe point is not judged repeatedly between the first starting catastrophe point and the first resetting catastrophe point, the starting catastrophe point can be judged continuously after the first resetting catastrophe point, and the nth starting catastrophe point and the resetting catastrophe point of the channel are searched until the data of the last sampling point of the analog quantity channel;
step 6, the current channel with the maximum effective current value in the wave recording files of the team is the basic channel for calculating the zero crossing point, and the accurate time deviation t _ diff between the IED and the comparison IED is calculated according to the starting mutation point of the basic channel of the wave recording files of the team, and the method comprises the following steps:
step 6.1, selecting current channels with established association relations in the channel configuration in the wave recording files of the team in the management unit, and finding the current channel with the maximum current effective value in the wave recording files of the team as a basic channel for calculating a zero crossing point;
step 6.2, starting mutation points of respective basic channels of the wave recording files of the team are taken, the time difference of the starting mutation points of the respective basic channels of the wave recording files of the team after time deviation correction firstly meets a pair of starting mutation points which are smaller than 5ms, and the starting mutation points are respectively defined as a first starting mutation point A (startA) and a first starting mutation point B (startB), and the time of the first starting mutation point A (startA) and the time of the first starting mutation point B (startB) after time deviation correction are respectively t1 and t 2;
step 6.3, finding a first rising trend zero-crossing point after the basic channels of the wave recording files of the team are added by 60ms by t 1;
step 6.4, if the zero-crossing point is not at the sampling point, calculating the time of the first rising trend zero-crossing point of the basic channels of the local IED and the comparative IED found in the step 6.3 by using an equal ratio difference method to be t3 and t4 respectively;
step 6.5, the specific zero crossing point time calculation method is shown in fig. 3: searching two continuous points, when the instantaneous value of the Xth sampling point is less than 0 and the instantaneous value of the X +1 sampling point is greater than 0, the zero crossing point is between the two sampling points, and the time t0 of the zero crossing point can be obtained by the following formula, wherein txTime of the Xth sample point, tx+1Time of the X +1 th sampling point, fxInstantaneous value of the Xth sample point, fx+1Instantaneous value of the X +1 sampling point;
t0＝(fx+1*tx-fx*tx+1)/(fx+1-fx)
step 6.6, in order to prevent the error of the whole cycle wave, the following method can be used for correcting: if (t3-t1) - (t4-t1) > < 10000us, t3 ═ t3-20000us, if (t4-t1) - (t3-t1) > < 10000us, t4 ═ t4-20000 us;
step 6.7, calculating the accurate time deviation t _ diff between the local IED and the comparison IED t3-t 4, then shifting the time of the IED corresponding to the first start-up discontinuity point b (startb) by the accurate time deviation t _ diff, thereby completing the accurate alignment with the time of the IED where the first start-up discontinuity point a (starta) is located, so as to be used for the subsequent calculation.
And 7, precisely aligning the wave recording files of the team in the management unit based on the precise time deviation t _ diff of the step 6, combining overlapping sections with time coincidence between all grouped startup mutation points and all grouped recovery mutation points in the wave recording files of the team found in the step 5, wherein the combination is the time when the initial time of the overlapping sections is the earliest mutation point, and the termination time of the overlapping sections is the time of the latest recovery point.
And 8, the waveform of the wave recording file of the protection device is different from the waveform of the wave recording file of the acquisition unit under the normal condition due to the fact that the protection device physically has pre-filtering, and when the management unit filters and processes the control word 1, the sampling point data of the wave recording file of the acquisition unit is digitally filtered according to the threshold value of the filtering frequency. And when the filtering processing control word is set to be 0, the wave recording data of the acquisition unit is not subjected to digital filtering.
Step 9, selecting a wave recording file with a low sampling rate from the wave recording files of the team in the management unit as a reference wave recording file; the other wave recording file in the wave recording files of the team is a non-reference wave recording file, each sampling point in the reference wave recording file and each sampling point in the non-reference wave recording file are corrected according to the accurate time deviation t _ diff, the sampling point of an analog channel in the reference wave recording file is a reference sampling point, a corresponding comparison sampling point is found in the non-reference wave recording file, the reference sampling point is compared with the comparison sampling point, and the time of a switching value wave recording channel participating in comparison is corrected and aligned according to the accurate time deviation t _ diff; comparing the time of the associated variable sites of the switching value recording channels, comprising the following steps:
step 9.1, selecting a wave recording file with a low sampling rate from the wave recording files of the team as a reference wave recording file in the management unit; the other wave recording file in the wave recording files of the team is a non-reference wave recording file, and each sampling point in the reference wave recording file and the non-reference wave recording file is corrected according to the accurate time deviation t _ diff obtained in the step 6.7;
and 9.2, taking the sampling point of the analog channel in the reference wave recording file as a reference sampling point, finding the sampling point by using the time of the reference sampling point of the analog channel of the reference wave recording file to the analog channel in the non-reference wave recording file, if the sampling point corresponding to the time is found, marking as a comparison sampling point, directly comparing the reference sampling point with the comparison sampling point, if the sampling point corresponding to the time cannot be found, calculating a comparison sampling point by using an interpolation algorithm, and comparing the reference sampling point with the comparison sampling point. In particular, a cubic spline interpolation algorithm can be used to output the calculated instantaneous value at this time;
step 9.3, comparing the reference sampling point with the comparison sampling point comprises: comparing the reference sampling point with the instantaneous value of the comparison sampling point, comparing the effective amplitude of the next cycle from the reference sampling point with the effective amplitude of the next cycle from the comparison sampling point, and comparing the phase angle corresponding to the reference sampling point and the comparison sampling point; the phase angle of the reference sampling point is obtained by calculating from the reference sampling point to the sampling point of the next cycle, and the phase angle of the comparison sampling point is obtained by calculating from the comparison sampling point to the sampling point of the next cycle.
9.4, comparing the difference value between the reference sampling point and the comparison sampling point instantaneous value obtained in the step 9.3 with an analog quantity instantaneous value comparison difference threshold value, and generating comparison abnormity when the amplitude difference between the effective amplitude of the next cycle from the reference sampling point and the effective amplitude of the next cycle from the comparison sampling point is greater than the analog quantity effective value amplitude comparison difference threshold value;
or the difference value between the reference sampling point and the comparison sampling point instantaneous value obtained in the step 9.3 is greater than the analog quantity instantaneous value comparison difference threshold, and the comparison abnormity is generated when the phase angle difference between the reference sampling point and the comparison sampling point is greater than the analog quantity effective value angle comparison difference threshold;
if the comparison abnormality occurs between the three continuous pairs of reference sampling points and the comparison sampling points, judging the channel comparison alarm;
step 9.6, acquiring the switching value wave recording channel for establishing the association relationship in the step 1 from the wave recording files participating in comparison;
9.7, correcting and aligning the time of the switching value recording channel participating in comparison according to the accurate time deviation t _ diff calculated in the step 6;
9.8, searching the time of all variable sites of the switching value recording channels participating in comparison;
and 9.9, comparing the time of the variable sites of the associated switching value recording channels, and giving an alarm when the time difference is greater than the switching value comparison difference threshold set in the step 1.
The specific embodiments described herein are merely illustrative of the invention. Various modifications, additions and substitutions may be made by those skilled in the art to the specific embodiments described without departing from the spirit of the invention or exceeding the scope of the claims defined below.
Claims (8)
1. A recording comparison analysis method based on time deviation management is characterized by comprising the following steps:
step 1, setting parameters in a management unit, wherein the parameters comprise a service automatic start period, a start signal time difference threshold, a sudden change start threshold, a sampling rate invalid threshold, an analog quantity instantaneous value comparison difference threshold, an analog quantity effective value amplitude comparison difference threshold, an analog quantity effective value angle comparison difference threshold, a switching quantity comparison difference threshold, a filtering processing control word and a filtering frequency, and the channel configuration associates analog quantity channels of the same sampling source in wave recording files of a protection device and an acquisition unit and associates switching quantity channels reflecting the displacement of the same switching quantity;
step 2, calculating the time deviation between the IED and the management unit;
step 3, the management unit acquires a comparison IED associated with each IED according to the channel configuration, matches the IED with a starting signal of the associated comparison IED, and acquires a corresponding recording file as a recording file to be compared;
step 4, correcting the recording time of the name of the wave recording file to be compared according to the time deviation of each IED, and grouping the IED and the wave recording files with overlapped wave recording file recording time sections of the compared IED to obtain wave recording files of a group;
step 5, respectively finding a starting mutation point and a resetting mutation point from wave recording files of the team in the associated analog quantity channel;
step 6, the current channel with the maximum effective current value in the wave recording files of the team is a basic channel for calculating a zero crossing point, and the accurate time deviation t _ diff between the IED and the comparison IED is calculated according to the starting break points of the basic channels of the wave recording files of the team;
7, finding out that all grouped start mutation points and all grouped recovery mutation points in the wave recording files of the group team have overlapped sections in time;
step 8, carrying out digital filtering on sampling point data of the wave recording file of the acquisition unit according to the filtering processing control word of the management unit;
step 9, selecting a wave recording file with a low sampling rate from the wave recording files of the team in the management unit as a reference wave recording file; the other wave recording file in the wave recording files of the team is a non-reference wave recording file, each sampling point in the reference wave recording file and each sampling point in the non-reference wave recording file are corrected according to the accurate time deviation t _ diff, the sampling point of an analog channel in the reference wave recording file is a reference sampling point, a corresponding comparison sampling point is found in the non-reference wave recording file, the reference sampling point is compared with the comparison sampling point, and the time of a switching value wave recording channel participating in comparison is corrected and aligned according to the accurate time deviation t _ diff; and comparing the time of the variable sites of the associated switching value recording channels.
2. The recording comparison analysis method based on time deviation management of claim 1, wherein the time deviation in step 2 is obtained by the following steps;
the management unit acquires and stores the current time of the protection device and the current time of the acquisition unit in an NTP mode, the time for the management unit to send an NTP request is defined as T1, the time for the protection device or the acquisition unit to receive the NTP request is defined as T2, the reply time of the protection device or the acquisition unit is defined as T3, and the reply time received by the management unit is defined as T4, so that the network delay d and the time deviation delta T between the protection device or the acquisition unit and the management unit are calculated according to the following formula:
d＝(T2-T1)+(T4-T3)；
△t＝[(T2-T1)+(T3-T4)]/2。
3. the recording comparison analysis method based on time deviation management as claimed in claim 2, wherein said step 3 comprises the steps of:
step 3.1, a starting signal buffer queue is built for each IED in the management unit, and starting signals of the IEDs are buffered in the starting signal buffer queue;
step 3.2, the management unit acquires a comparison IED associated with each IED according to the channel configuration;
step 3.3, checking whether the IED start signal buffer queue is empty every time the time length defined by the start signal time difference threshold, and if not, matching the start signal in the IED start signal buffer queue to the comparison IED start signal buffer queue, specifically:
if the difference value of the starting signal time of the starting signals in the starting signal cache queue of the IED and the starting signal time of the starting signal in the starting signal cache queue of the comparison IED is smaller than the configured starting signal time difference threshold value, matching the starting signals in the starting signal cache queue of the IED and the starting signal time in the starting signal cache queue of the comparison IED, wherein the starting signal time is the starting signal time after the time deviation correction obtained in the step 2;
step 3.4, the unmatched starting signals in the starting signal buffer queue are cleared if the unmatched starting signals are not matched after three times of service automatic starting cycles;
step 3.5, finding a corresponding recording file from the successfully called recording file according to the matched IED starting signal and the compared IED starting signal, and if the starting signal time and the recording file name recording time are less than 50ms, judging that the recording file corresponding to the starting signal is found;
step 3.6, when the IED and the comparison IED both find the wave recording file corresponding to the matched starting signal, taking the two wave recording files as wave recording files to be compared;
and 3.7, if the wave recording file corresponding to the starting signal is not found in the period of the current service automatic starting, the wave recording file is found again in the period of the next service automatic starting, and the wave recording file corresponding to the matched starting signal of the IED and the comparison IED is discarded if the wave recording file corresponding to the starting signal is not found in the period of the three continuous service automatic starting.
4. The method according to claim 3, wherein the step 4 comprises the steps of:
step 4.1, the recording file name recording time of the recording files to be compared of the local IED and the comparison IED obtained in the step 3 is corrected according to the time deviation of the respective IED in the management unit, and the recording files to be compared of the local IED and the comparison IED are respectively sequenced according to the corrected recording file name recording time;
4.2, analyzing respective wave recording files to be compared of the IED and the comparison IED and acquiring a corresponding recording time section;
4.3, grouping the IED and the wave recording files with the overlapped wave recording file recording time sections of the comparison IED to obtain wave recording files of the group;
and 4.4, discarding the recording file of which the overlapping area is not found.
5. The method according to claim 4, wherein the step 5 comprises the steps of:
step 5.1, selecting analog quantity channels with the association relationship established in the channel configuration of the wave recording files grouped in the step 4 in the management unit;
step 5.2, in the wave recording file of the team, in the associated analog quantity channel, removing the wave recording section with the sampling rate lower than the sampling rate invalid threshold value in any analog quantity channel;
step 5.3, in the wave recording file of the team, removing wave recording sections with discontinuous occurrence time in any analog quantity channel from the related analog quantity channel;
step 5.4, in the wave recording file of the team, removing sampling points with instantaneous values smaller than 0.001 times of rated values in any analog quantity channel in the related analog quantity channel;
and 5.5, in the wave recording file of the team, in each remaining wave recording section in the associated analog quantity channel, starting from the first sampling point of the third cycle, when the sudden change instantaneous values of the continuous three sampling points are greater than a sudden change starting threshold value, the first sampling point of the continuous three sampling points is a starting sudden change point, and when the sudden change instantaneous values of the sampling points of the continuous cycle after the starting of the sudden change point are less than the sudden change starting threshold value, the first sampling point of the continuous cycle is a repeated sudden change point.
6. The recording comparison analysis method based on time deviation management according to claim 5, wherein the transient value of the mutation in step 5.5 is calculated as follows:
note ikFor the instantaneous sample value corresponding to the kth sample point, ik-NFor the instantaneous value, i, of the sampling point one sampling period before the kth sampling pointk-2NThe instantaneous value delta i of the sudden change of the kth sampling point is the instantaneous sampling value of the sampling point two sampling periods before the kth sampling pointkComprises the following steps:
Δik＝||ik-ik-N|-|ik-N-ik-2N||。
7. the method of claim 5, wherein the step 6 comprises the steps of:
step 6.1, selecting current channels with established association relations in the channel configuration in the wave recording files of the team in the management unit, and finding the current channel with the maximum current effective value in the wave recording files of the team as a basic channel for calculating a zero crossing point;
step 6.2, starting mutation points of respective basic channels of the wave recording files of the team are taken, the time difference of the starting mutation points of the respective basic channels of the wave recording files of the team after time deviation correction firstly meets a pair of starting mutation points which are smaller than 5ms, the starting mutation points are respectively defined as a first starting mutation point A and a first starting mutation point B, and the time of the first starting mutation point A and the time of the first starting mutation point B after time deviation correction are respectively t1 and t 2;
step 6.3, finding a first rising trend zero-crossing point after the basic channels of the wave recording files of the team are added by 60ms by t 1;
step 6.4, if the zero-crossing point is not at the sampling point, calculating the time of the first rising trend zero-crossing point of the basic channels of the local IED and the comparative IED found in the step 6.3 by using an equal ratio difference method to be t3 and t4 respectively;
if (t3-t1) - (t4-t1) > < 10000us, t3 ═ t3-20000us, if (t4-t1) - (t3-t1) > < 10000us, t4 ═ t4-20000 us;
the precise time offset t _ diff between the present IED and the comparison IED is calculated as t3-t 4.
8. The recording comparison analysis method based on time deviation management as claimed in claim 7, wherein said step 9 comprises the steps of:
step 9.1, selecting a wave recording file with a low sampling rate from the wave recording files of the team as a reference wave recording file in the management unit; the other wave recording file in the wave recording files of the team is a non-reference wave recording file, and each sampling point in the reference wave recording file and the non-reference wave recording file is corrected according to the accurate time deviation t _ diff obtained in the step 6.4;
step 9.2, taking the sampling point of the analog channel in the reference wave recording file as a reference sampling point, finding the sampling point by using the time of the reference sampling point of the analog channel of the reference wave recording file to the analog channel in the non-reference wave recording file, if the sampling point corresponding to the time is found, marking as a comparison sampling point, directly comparing the reference sampling point with the comparison sampling point, if the sampling point corresponding to the time cannot be found, calculating a comparison sampling point by using an interpolation algorithm, and comparing the reference sampling point with the comparison sampling point;
step 9.3, comparing the reference sampling point with the comparison sampling point comprises: comparing the reference sampling point with the instantaneous value of the comparison sampling point, comparing the effective amplitude of the next cycle from the reference sampling point with the effective amplitude of the next cycle from the comparison sampling point, and comparing the phase angle corresponding to the reference sampling point and the comparison sampling point;
9.4, comparing the difference value between the reference sampling point and the comparison sampling point instantaneous value obtained in the step 9.3 with an analog quantity instantaneous value comparison difference threshold value, and generating comparison abnormity when the amplitude difference between the effective amplitude of the next cycle from the reference sampling point and the effective amplitude of the next cycle from the comparison sampling point is greater than the analog quantity effective value amplitude comparison difference threshold value;
or the difference value between the reference sampling point and the comparison sampling point instantaneous value obtained in the step 9.3 is greater than the analog quantity instantaneous value comparison difference threshold, and the comparison abnormity is generated when the phase angle difference between the reference sampling point and the comparison sampling point is greater than the analog quantity effective value angle comparison difference threshold;
if the comparison abnormality occurs between the three continuous pairs of reference sampling points and the comparison sampling points, judging the channel comparison alarm;
step 9.6, acquiring the switching value wave recording channel for establishing the association relationship in the step 1 from the wave recording files participating in comparison;
9.7, correcting and aligning the time of the switching value recording channel participating in comparison according to the accurate time deviation t _ diff calculated in the step 6;
9.8, searching the time of all variable sites of the switching value recording channels participating in comparison;
and 9.9, comparing the time of the variable sites of the associated switching value recording channels, and giving an alarm when the time difference is greater than the switching value comparison difference threshold set in the step 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011296858.9A CN112506867B (en) | 2020-11-18 | 2020-11-18 | Recording comparison analysis method based on time deviation management |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011296858.9A CN112506867B (en) | 2020-11-18 | 2020-11-18 | Recording comparison analysis method based on time deviation management |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112506867A CN112506867A (en) | 2021-03-16 |
| CN112506867B true CN112506867B (en) | 2022-02-18 |
Family
ID=74956979
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011296858.9A Active CN112506867B (en) | 2020-11-18 | 2020-11-18 | Recording comparison analysis method based on time deviation management |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112506867B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113447745B (en) * | 2021-06-28 | 2024-03-01 | 中国船舶重工集团公司第七0四研究所 | Station domain virtual fault wave recording processing method |
| CN113805050B (en) * | 2021-08-31 | 2024-04-02 | 中国南方电网有限责任公司超高压输电公司广州局 | Phase selection closing angle monitoring method and device, computer equipment and storage medium |
| CN113722145B (en) * | 2021-11-04 | 2022-03-11 | 国网江西省电力有限公司电力科学研究院 | Method and device for rapidly identifying homologous anomaly in vacuum semantic environment |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108469561A (en) * | 2018-03-21 | 2018-08-31 | 成都滕达科技有限公司 | A kind of intelligent substation Commissioning of Relay Protection method |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5669629B2 (en) * | 2010-05-18 | 2015-02-12 | キヤノン株式会社 | Terahertz wave measuring apparatus and measuring method |
| CN103439607B (en) * | 2013-08-28 | 2016-11-23 | 三川电力设备股份有限公司 | By method and system and the Fault Locating Method of failure wave-recording identification element population parameter |
-
2020
- 2020-11-18 CN CN202011296858.9A patent/CN112506867B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108469561A (en) * | 2018-03-21 | 2018-08-31 | 成都滕达科技有限公司 | A kind of intelligent substation Commissioning of Relay Protection method |
Non-Patent Citations (1)
| Title |
|---|
| 基于广域故障录波信息的调度端电网故障诊断系统;李乃永;《电力系统自动化》;20140825;第38卷(第16期);第100-104页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112506867A (en) | 2021-03-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112506867B (en) | Recording comparison analysis method based on time deviation management | |
| US9128140B2 (en) | Detection of a fault in an ungrounded electric power distribution system | |
| EP2156203B1 (en) | Method and device to predict a state of a power system in the time domain | |
| CN103201923A (en) | Synchronization method for current differential protection | |
| US7583771B2 (en) | Systems and methods for resampling unreliable data | |
| CN109742741B (en) | Relay protection data synchronization method for active power distribution network | |
| CN110865278B (en) | Ground fault positioning method based on transient mutation energy capturing method | |
| US11258249B2 (en) | Primary and system protection for an electric power delivery system | |
| CN1656359A (en) | Instantaneous overcurrent element for heavily saturated current in a power system | |
| CN109975653B (en) | 10kV distribution line fault location method | |
| CN107872049A (en) | Method and apparatus for controlling the breaker in power supply network | |
| CN111426908B (en) | Single-phase earth fault protection method, device and system for small current earthing system | |
| US7742884B2 (en) | Sampling frequency control method and protective relay | |
| CN113281618A (en) | Low-voltage distribution line fault positioning method and device | |
| JP2015023742A (en) | Protection monitoring control system | |
| CN109142977B (en) | Phase automatic detection method of transient recording type fault indicator | |
| JP2002040067A (en) | Detection method for inter-order harmonics | |
| CN111987699B (en) | Current differential protection data synchronization method of alternating current line and related device | |
| CN113671308A (en) | Abnormal data screening method and system for digital acquisition device of power distribution switch | |
| CN110024249B (en) | Method for detecting a fault in an electric power transmission line and protection system using the method | |
| US11114892B2 (en) | Electric power system transducer failure monitor and measurement recovery | |
| JP4316103B2 (en) | Fault location system | |
| CN111157917A (en) | Low-current grounding line selection system and device and grounding fault positioning method | |
| JPH0345344B2 (en) | ||
| EP0957559A2 (en) | Method of synchronizing line differential protection device, and line differential protection device |
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 |