CN115267668A - Automatic linear positioning system and method for GIS partial discharge - Google Patents

Automatic linear positioning system and method for GIS partial discharge Download PDF

Info

Publication number
CN115267668A
CN115267668A CN202210824993.9A CN202210824993A CN115267668A CN 115267668 A CN115267668 A CN 115267668A CN 202210824993 A CN202210824993 A CN 202210824993A CN 115267668 A CN115267668 A CN 115267668A
Authority
CN
China
Prior art keywords
pulse
signal
partial discharge
positioning
signals
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202210824993.9A
Other languages
Chinese (zh)
Inventor
徐洪海
张勇
刘文松
胡飞
潘尚举
李煜
许超
潘宏晨
汤水成
王伟伟
于立岩
王玉玮
夏清普
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Jiangsu Electric Power Co ltd Nanjing Jiangning District Power Supply Branch
Jiangsu Hongyuan Electric Co Ltd
Original Assignee
State Grid Jiangsu Electric Power Co ltd Nanjing Jiangning District Power Supply Branch
Jiangsu Hongyuan Electric Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Jiangsu Electric Power Co ltd Nanjing Jiangning District Power Supply Branch, Jiangsu Hongyuan Electric Co Ltd filed Critical State Grid Jiangsu Electric Power Co ltd Nanjing Jiangning District Power Supply Branch
Priority to CN202210824993.9A priority Critical patent/CN115267668A/en
Publication of CN115267668A publication Critical patent/CN115267668A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
    • G01S5/06—Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
    • 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/12—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing
    • G01R31/1227—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing of components, parts or materials
    • G01R31/1254—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing of components, parts or materials of gas-insulated power appliances or vacuum gaps
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
    • G01S5/0278—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves involving statistical or probabilistic considerations

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Probability & Statistics with Applications (AREA)
  • Testing Relating To Insulation (AREA)

Abstract

The invention discloses a GIS partial discharge automatic linear positioning system and a method, wherein the system comprises a preprocessing module, a pulse extraction module, a signal sorting module and a positioning module; the method comprises the following steps: filtering each pulse signal sample; carrying out numerical value normalization processing; extracting a pulse part in each signal; adjusting the displacement of multiple data waveforms; obtaining a pulse initial position, and calculating a partial discharge source positioning value; and (5) positioning position statistics. The system and the method can process the acquired signals to obtain the accurate position of the partial discharge source.

Description

Automatic linear positioning system and method for GIS partial discharge
Technical Field
The invention relates to electrified detection of power equipment, in particular to a GIS partial discharge automatic linear positioning system and a method.
Background
Inside because of insulating decline, reasons such as manufacturing quality, can produce partial discharge under high-voltage electric field of GIS combined electrical apparatus, if overhauls inside the GIS air chamber, the expense is very expensive, consequently accurately confirms type, the position of defect in advance, is favorable to accurate quick development maintenance work.
The partial discharge signal can be measured by an ultrahigh frequency sensor placed at the pouring hole. The ultrahigh frequency sensors are arranged at the pouring holes on the two sides of the partial discharge source, the receiving time is inconsistent due to different positions away from the partial discharge source, and the position of the partial discharge source can be calculated according to the time difference of signals received by the two sensors.
However, the partial discharge signal is affected by the interference signal, the propagation path, the difference of the transmission parameters of the sensor, and the like, as shown in fig. 3, the received partial discharge pulse signal has a certain difference, and therefore, the time difference can be accurately obtained only when the pulse start point position needs to be obtained. Since the propagation speed of the electromagnetic wave signal is the speed of light, the selection of the starting point position of the pulse is inaccurate, and a positioning error of 300mm will be generated if the error of 1ns is generated, that is, whether the selection of the starting point position is accurate or not relates to the final error of the positioning result.
Disclosure of Invention
The purpose of the invention is as follows: the invention aims to provide a GIS partial discharge automatic linear positioning system and a method, so that collected signals are processed to obtain an accurate partial discharge source position.
The technical scheme is as follows: the invention relates to a GIS partial discharge automatic linear positioning system, which comprises the following modules:
a preprocessing module: the device is used for filtering and normalizing the signals;
a pulse extraction module: the partial discharge pulse signal extraction device is used for extracting partial discharge pulse signals from the collected signal pair samples S1 and S2;
a signal arrangement module: the device is used for carrying out interpolation and pulse alignment operations on signals;
a positioning module: estimating the starting point of the partial discharge pulse based on an average energy accumulation method, and performing positioning calculation; and the method is used for carrying out probability statistics on the multiple positioning results to obtain accurate positioning results.
A GIS partial discharge automatic linear positioning method comprises the following steps:
(1) Filtering each pulse signal sample;
(2) Carrying out numerical value normalization processing;
(3) Extracting a pulse part in each signal;
(4) Adjusting the displacement of multiple data waveforms;
(5) Obtaining a pulse initial position, and calculating a partial discharge source positioning value;
(6) And (5) positioning position statistics.
The step (1) is specifically as follows: filtering the signals S1 and S2 by an IIR (infinite impulse response) band-pass filtering calculation method, wherein S1 represents a plurality of groups of pulse signals collected by the sensor 1, and S1= { S = (zero mean of zero) is adopted11，S12,...,S1NS2 represents a plurality of groups of pulse signals collected by the sensor 2, S2= { S = }21，S22,...,S2NS1 or S2, or any element of SijI = {1,2}, which means signal 1, signal 2, j = {1, 2., N }, respectively, representing pulse train sequence numbers, where S = {1, 2., } represents the pulse train sequence numbers1j、S2jA set of pulse signal pairs, S, acquired by the sensor 1 and the sensor 2 at the same momentijM discrete sequences { a) representing a pulse signal waveform1,a2,...,aM}。
The step (2) is specifically as follows: taking the maximum value A in all pulse sequences of S1 and S2maxNormalizing each element according to a formula 1;
Ai＝ai/Amax (1)
wherein, aiFor a sample value in any one of the discrete sequences S1 or S2, AiIs a normalized value.
The step (3) is specifically as follows: setting the window width as p, moving according to the step length p, respectively calculating the sheath degree K through a formula 2 to quickly determine the position where the partial discharge pulse exists, and when the sheath degree is greater than 3, indicating that the position has signal mutation;
Figure BDA0003746119270000021
wherein μ is the signal Sijσ is the signal SijP is the window width; when more impact components exist in the signal, the kurtosis value is obviously increased, and the kurtosis value is larger when the impact is larger;
respectively moving forward and backward point by taking the central point of the window as a starting point and taking the window with the width of P1, integrating the value in the window P1 according to a formula 3 to obtain Em, and calculating EmAfter being respectively smaller than E1 and E2, the left side of the pulse is obtainedStarting number XijRight starting sequence number Yij，Xij、YijThe sequence between them is the partial discharge pulse signal to be extracted, and records each { X }ij,YijWhere i = {1,2}, which represents signals of the sensor 1 and the sensor 2, respectively, and j = {1, 2.., N }, which represents a set of signals acquired by the sensor 1 or the sensor 2;
Figure BDA0003746119270000022
taking X in all pulsesijTaking Y as the starting position of the minimum valueijAnd the medium maximum value is used as an end position, and all signals in the S1 and the S2 which are equal in length are uniformly intercepted.
The step (4) is specifically as follows:
(4.1) calculating a cross-correlation function between each pulse signal in S1 according to formula 4, respectively;
Figure BDA0003746119270000031
in the formula, i and j are respectively the ith sample and the jth sample of the same channel; k is a sample serial number; s is the length of the signal; m is the number of points of the shift, corresponding to the time shift A in the continuous signalj；AiAnd AjRepresenting waveform data acquired by two different sensors;
(4.2) calculating the sum M of the displacements between each sample and the remaining samples of S1 according to equation 5i；
Figure BDA0003746119270000032
Wherein i = {1,2, ·, n }, i denotes the ith pulse in S1, j = {1,2,. J,..., n-1}, j denotes the jth pulse in S1;
(4.3) with MiThe pulse in S1 corresponding to the maximum value is used as a reference, and the other samples are shift-aligned by the displacement thereof.
The step (5) is specifically as follows:
(5.1) calculating the pulse S of each set by equation 61j、S2jEnergy accumulation curve of (d);
Figure BDA0003746119270000033
Figure BDA0003746119270000034
wherein E iskIntegrating the 1 st sample to the kth sample, wherein delta is a deviation coefficient and is used for deflecting the energy accumulation curve in a negative direction, and the value range of a coefficient a is 1-5;
(5.2) respectively calculating an energy accumulation mean value curve of each group of pulse signals 1 and 2 through a formula 8;
Figure BDA0003746119270000035
in the formula 8, the first step is,
Figure BDA0003746119270000036
representing M groups of energy accumulation curves, wherein k = {1,2}, and represents energy accumulation curves corresponding to S1 and S2 signals respectively;
(5.3) are respectively paired
Figure BDA0003746119270000037
After the first derivative is obtained, the position of the minimum value is the position of the starting point of the pulse signal, and the starting time t of the pulse signal is obtained1、t2And obtaining the time difference between the two signals, delta t = t2-t1；
(5.4) calculating the position of the partial discharge source relative to the signal 1 sensor through a formula 9;
Figure BDA0003746119270000041
wherein, x is the distance of the partial discharge source relative to the sensor 1, L is the distance between the sensor 1 and the sensor 2, and v is the propagation velocity of the signal.
The step (6) is specifically as follows: repeating the steps (1) to (4), calculating to obtain each x value, if the number of x is N, dividing the length L into intervals with the width of w, and respectively counting the number P of the x values in the intervalsi，PiThe middle point of the maximum interval is the actual position of the partial discharge source.
A computer storage medium having stored thereon a computer program which, when executed by a processor, implements a method of automatic linear positioning for GIS partial discharges as described above.
A computer device comprises a storage, a processor and a computer program stored on the storage and running on the processor, wherein the processor implements the above method for automatic linear positioning of GIS partial discharge when executing the computer program.
Has the beneficial effects that: compared with the prior art, the invention has the following advantages: 1. the accuracy of positioning based on time difference is improved by multi-signal comprehensive positioning; 2. the GIS partial discharge positioning automation is realized, and the technical requirements and labor intensity of field operation are reduced.
Drawings
FIG. 1 is a flow chart of automatic positioning;
FIG. 2 is a schematic view of a partial discharge positioning;
FIG. 3 is a diagram showing a comparison of dual-channel partial discharge pulse waveforms;
FIG. 4 is a graph of energy mean;
fig. 5 is a positioning statistical characteristic diagram, in which fig. 5 (a) is a diagram of a single-point-calculation positioning result, fig. 5 (b) is a diagram of a multi-pulse-calculation positioning result, fig. 5 (c) is a diagram of a single-point-calculation positioning statistic, and fig. 5 (d) is a diagram of a multi-pulse-calculation positioning statistic.
Detailed Description
The technical scheme of the invention is further explained by combining the attached drawings.
A GIS partial discharge automatic linear positioning system comprises the following modules:
a preprocessing module: the device is used for filtering and normalizing the signals;
a pulse extraction module: the partial discharge pulse signal extraction device is used for extracting partial discharge pulse signals from the collected signal pair samples S1 and S2;
a signal arrangement module: the device is used for carrying out difference and pulse alignment operation on the signals;
a positioning module: estimating the starting point of the partial discharge pulse based on an average energy accumulation method, and performing positioning calculation; and the method is used for carrying out probability statistics on a plurality of positioning results to obtain an accurate positioning result.
The sensor placement method is shown in figure 2, in the figure S1 represents the groups of pulse signals collected by the sensor 1, S1= { S =11，S12,...,S1NS2 represents sets of pulse signals collected by the sensor 2, S2= { S = }21，S22,...,S2N}. Any one element S of S1 or S2ijI = {1,2}, which means signal 1, signal 2, j = {1, 2., N }, respectively, representing pulse train sequence numbers, where S = {1, 2., } represents the pulse train sequence numbers1j、S2jA set of pulse signal pairs, S, acquired by the sensor 1 and the sensor 2 at the same momentijM discrete sequences { a) representing a pulse signal waveform1,a2,...,aM}. And filtering the signals S1 and S2, wherein an IIR band-pass filtering calculation method is adopted in the invention.
As shown in fig. 1, a method for automatic linear positioning of GIS partial discharge includes the following steps:
(1) Each pulse signal sample is filtered.
Each pulse signal sample of S1, S2 is filtered. Various effective filtering calculation methods can be adopted, and a common band-pass filtering calculation method is adopted in the invention, wherein the stop band boundary is Ws = [250MHz,1550MHz ], and the band-pass boundary Wp = [300MHz,1500MHz ].
(2) And (6) carrying out numerical value normalization processing.
Taking the maximum value A in all pulse sequences of S1 and S2maxNormalizing each element according to a formula 1;
Ai＝ai/Amax (1)
wherein, aiFor a sample value in any one of the discrete sequences S1 or S2，AiIs aiNormalized values.
(3) The pulse part in each signal is extracted.
Setting the window width as p, moving according to the step length p, respectively calculating the sheath degree K through a formula 2 to quickly determine the position where the partial discharge pulse exists, and when the sheath degree is greater than 3, indicating that the position has signal mutation;
Figure BDA0003746119270000051
wherein μ is the signal Sijσ is the signal SijP is the window width; when more impact components exist in the signal, the kurtosis value is obviously increased, and the larger the impact, the larger the kurtosis value is.
Respectively moving forward and backward point by taking the central point of the window as a starting point and taking the window with the width of P1, integrating the value in the window P1 according to a formula 3 to obtain Em, and calculating EmRespectively obtaining the starting sequence number X at the left side of the pulse after being smaller than E1 and E2ijRight starting sequence number Yij，Xij、YijThe sequence between them is the partial discharge pulse signal to be extracted, and records each { X }ij,YijWhere i = {1,2}, which represents signals of the sensor 1 and the sensor 2, respectively, and j = {1, 2.., N }, which represents a set of signals acquired by the sensor 1 or the sensor 2.
Figure BDA0003746119270000061
Taking X in all pulsesijTaking Y as the starting position of the minimum valueijAnd the medium maximum value is used as an end position, and all signals in the S1 and the S2 which are equal in length are uniformly intercepted.
(4) And adjusting the displacement of multiple data waveforms.
(4.1) respectively calculating a cross-correlation function between each pulse signal in S1 according to formula 4;
Figure BDA0003746119270000062
in the formula, i and j are respectively the ith sample and the jth sample of the same channel; k is a sample serial number; s is the length of the signal; m is the number of points of the shift, corresponding to the time shift A in the continuous signalj；AiAnd AjRepresenting waveform data acquired by two different sensors.
(4.2) calculating the sum M of the displacements between each sample and the remaining samples of S1 according to equation 5i；
Figure BDA0003746119270000063
Wherein i = {1,2, ·, n }, i denotes the ith pulse in S1, j = {1,2,. J,..., n-1}, and j denotes the jth pulse in S1.
(4.3) with MiThe pulse in S1 corresponding to the maximum value is used as a reference, and the other samples are shift-aligned by the displacement thereof.
(5) And obtaining the initial position of the pulse, and calculating the positioning value of the partial discharge source.
(5.1) calculating the pulse S of each set by equation 61j、S2jThe energy accumulation curve of (a) is shown in fig. 3.
Figure BDA0003746119270000064
Figure BDA0003746119270000065
Wherein the content of the first and second substances,
Figure BDA0003746119270000066
integrating the 1 st sample to the kth sample, wherein delta is a deviation coefficient and is used for deflecting the energy accumulation curve in a negative direction, and the value range of a coefficient a is 1-5; in this example a =3.
(5.2) separately calculating the energy accumulation mean value curve of each set of pulse signal 1 and signal 2 by equation 8 is shown in fig. 4.
Figure BDA0003746119270000067
In the formula 8, the first and second groups of the compound,
Figure BDA0003746119270000071
the curves representing the M sets of energy accumulation mean values, k = {1,2}, which represent the energy accumulation curves corresponding to the S1 and S2 signals, respectively.
(5.3) are respectively paired
Figure BDA0003746119270000072
After the first derivative is obtained, the position of the minimum value is the position of the starting point of the pulse signal, and the starting time t of the pulse signal is obtained1、t2And obtaining the time difference between the two signals, delta t = t2-t1。
(5.4) calculating the position of the partial discharge source relative to the signal 1 sensor by formula 9.
Figure BDA0003746119270000073
Wherein, x is the distance of the partial discharge source relative to the sensor 1, L is the distance between the sensor 1 and the sensor 2, and v is the propagation velocity of the signal.
(6) And (5) positioning position statistics.
And (5) repeating the steps (1) to (4), and calculating to obtain each value x, wherein the distribution of the values x is shown in the attached figure 5. If the number of x is N, dividing the length L into intervals with the width of w =0.05m, and respectively counting the number P of x values in the intervalsi，PiThe middle point of the maximum interval is the actual position of the partial discharge source.

Claims (10)

1. The automatic linear positioning system for the GIS partial discharge is characterized by comprising the following modules:
a pretreatment module: the signal processing device is used for filtering and normalizing signals;
a pulse extraction module: the partial discharge pulse signal extraction device is used for extracting partial discharge pulse signals from the collected signal pair samples S1 and S2;
a signal arrangement module: the device is used for carrying out interpolation and pulse alignment operations on signals;
a positioning module: estimating the starting point of the partial discharge pulse based on an average energy accumulation method, and performing positioning calculation; and the method is used for carrying out probability statistics on the multiple positioning results to obtain accurate positioning results.
2. A method for automatic linear positioning of GIS partial discharge, said method using the system of claim 1, comprising the steps of:
(1) Filtering each pulse signal sample;
(2) Carrying out numerical value normalization processing;
(3) Extracting a pulse part in each signal;
(4) Adjusting the displacement of multiple data waveforms;
(5) Obtaining a pulse initial position, and calculating a partial discharge source positioning value;
(6) And (5) positioning position statistics.
3. The method for automatically linearly positioning the GIS partial discharge according to claim 2, wherein the step (1) is specifically as follows: filtering the signals S1 and S2 by adopting an IIR (infinite impulse response) band-pass filtering calculation method, wherein S1 represents a plurality of groups of pulse signals collected by the sensor 1, and S1= { S =11，S12,...,S1NS2 represents sets of pulse signals collected by the sensor 2, S2= { S = }21，S22,...,S2NS is any one element of S1 or S2ijI = {1,2}, which means signal 1, signal 2, j = {1, 2., N }, respectively, representing pulse train sequence numbers, where S = {1, 2., } represents the pulse train sequence numbers1j、S2jA set of pulse signal pairs, S, acquired by the sensors 1 and 2 at the same momentijM discrete sequences { a) representing a pulse signal waveform1,a2,...,aM}。
4. According to the rightThe method for automatically linearly positioning the GIS partial discharge according to claim 2, wherein the step (2) specifically comprises: taking the maximum value A in all pulse sequences of S1 and S2maxNormalizing each element according to a formula 1;
Ai＝ai/Amax(1)
wherein, aiIs a sample value of any one of the discrete sequences of S1 or S2, AiIs a normalized value.
5. The method for automatically linearly positioning the GIS partial discharge according to claim 2, wherein the step (3) is specifically as follows: setting the window width as p, moving according to the step length p, respectively calculating the sheath degree K through a formula 2 to quickly determine the position where the partial discharge pulse exists, and when the sheath degree is greater than 3, indicating that the position has signal mutation;
Figure FDA0003746119260000021
wherein μ is the signal Sijσ is the signal SijP is the window width; when more impact components exist in the signal, the kurtosis value is obviously increased, and the kurtosis value is larger when the impact is larger;
respectively moving forward and backward point by taking the central point of the window as a starting point and taking the window with the width of P1, integrating the value in the window P1 according to a formula 3 to obtain Em, and calculating EmRespectively less than E1 and E2, obtaining the starting sequence number X of the left side of the pulseijRight start sequence number Yij，Xij、YijThe sequence between them is the partial discharge pulse signal to be extracted, and records each { X }ij,YijWherein i = {1,2}, representing signals of the sensor 1 and the sensor 2, respectively, and j = {1, 2.. Times, N }, representing a set of signals acquired by the sensor 1 or the sensor 2;
Figure FDA0003746119260000022
taking X in all pulsesijTaking Y as the starting position of the minimum valueijAnd the medium maximum value is used as an end position, and all signals in the S1 and the S2 which are equal in length are uniformly intercepted.
6. The method for automatically linearly positioning the GIS partial discharge according to claim 2, wherein the step (4) is specifically as follows:
(4.1) respectively calculating a cross-correlation function between each pulse signal in S1 according to formula 4;
Figure FDA0003746119260000023
in the formula, i and j are respectively the ith sample and the jth sample of the same channel; k is a sample serial number; s is the length of the signal; m is the number of points of the shift, corresponding to the time shift A in the continuous signalj；AiAnd AjRepresenting waveform data obtained by two different sensors;
(4.2) calculating the sum M of the displacements between each sample and the remaining samples of S1 according to equation 5i；
Figure FDA0003746119260000024
Wherein i = {1,2, ·, n }, i denotes the ith pulse in S1, j = {1,2,. J,..., n-1}, j denotes the jth pulse in S1;
(4.3) with MiThe pulse in S1 corresponding to the maximum value is used as a reference, and other samples are shift-aligned by the displacement thereof.
7. The method according to claim 2, wherein the step (5) specifically comprises:
(5.1) calculating the pulse S of each set by equation 61j、S2jEnergy accumulation curve of (d);
Figure FDA0003746119260000031
Figure FDA0003746119260000032
wherein, EkIntegrating the 1 st sample to the kth sample, wherein delta is a deviation coefficient and is used for deflecting the energy accumulation curve in a negative direction, and the value range of a coefficient a is 1-5;
(5.2) respectively calculating an energy accumulation mean value curve of each group of pulse signals 1 and 2 through a formula 8;
Figure FDA0003746119260000033
in the formula 8, the first and second groups of the compound,
Figure FDA0003746119260000034
representing M groups of energy accumulation mean curves, wherein k = {1,2}, and represents energy accumulation curves corresponding to S1 and S2 signals respectively;
(5.3) are respectively paired
Figure FDA0003746119260000035
After the first derivative is obtained, the position of the minimum value is the position of the starting point of the pulse signal, and the starting time t of the pulse signal is obtained1、t2And obtaining the time difference between the two signals, delta t = t2-t1；
(5.4) calculating the position of the partial discharge source relative to the signal 1 sensor through a formula 9;
Figure FDA0003746119260000036
wherein, x is the distance of the partial discharge source relative to the sensor 1, L is the distance between the sensor 1 and the sensor 2, and v is the propagation velocity of the signal.
8. The method for automatically linearly positioning the GIS partial discharge according to claim 2, wherein the step (6) is specifically as follows: repeating the steps (1) to (4), calculating to obtain each x value, if the number of x is N, dividing the length L into intervals with the width of w, and respectively counting the number P of the x values in the intervalsi，PiThe middle point of the maximum interval is the actual position of the partial discharge source.
9. A computer storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements a method for automatic linear positioning of GIS partial discharges according to any of claims 2-8.
10. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor when executing the computer program implements a method for automatic linear positioning of GIS partial discharges according to any of claims 2-8.
CN202210824993.9A 2022-07-14 2022-07-14 Automatic linear positioning system and method for GIS partial discharge Pending CN115267668A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210824993.9A CN115267668A (en) 2022-07-14 2022-07-14 Automatic linear positioning system and method for GIS partial discharge

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210824993.9A CN115267668A (en) 2022-07-14 2022-07-14 Automatic linear positioning system and method for GIS partial discharge

Publications (1)

Publication Number Publication Date
CN115267668A true CN115267668A (en) 2022-11-01

Family

ID=83766219

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210824993.9A Pending CN115267668A (en) 2022-07-14 2022-07-14 Automatic linear positioning system and method for GIS partial discharge

Country Status (1)

Country Link
CN (1) CN115267668A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116879699A (en) * 2023-09-08 2023-10-13 北京大学 Target object determining method, device, arc detecting system and storage medium

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116879699A (en) * 2023-09-08 2023-10-13 北京大学 Target object determining method, device, arc detecting system and storage medium
CN116879699B (en) * 2023-09-08 2023-12-15 北京大学 Target object determining method, device, arc detecting system and storage medium

Similar Documents

Publication Publication Date Title
CN109085477B (en) Signal identification and positioning method for power cable distributed partial discharge monitoring system
CN109522801B (en) Method and system for separating overlapped peaks of time-of-flight mass spectrometer
CN106443203A (en) Pulse signal detection system and method
CN115267668A (en) Automatic linear positioning system and method for GIS partial discharge
CN110682159A (en) Cutter wear state identification method and device
CN104942490A (en) Ultrasonic sensor based method and system for recognition and positioning of fillet weld initial point
CN113901379A (en) Dynamic online rapid processing method for real-time data of edge terminal
CN111610428B (en) Parameter optimization method based on response aliasing measurement wavelet packet decomposition algorithm
CN111599348B (en) Automatic segmentation method and system for machine tool machining process monitoring signals
CN117213569A (en) Ultrasonic flow detection method
CN107192840A (en) The space flight of feature based signal is after day device fifth wheel recognition methods
CN116028790A (en) Quick identification and screening method for ship mechanical vibration characteristic line spectrum
CN108008262B (en) Rapid partial discharge detection method
CN111833905B (en) System and method for detecting quality of marked character based on audio analysis
CN115293188A (en) Fault diagnosis method and device for reciprocating mechanical equipment
CN112444562B (en) Non-reference guided wave damage monitoring method adopting symmetrical structure response signal analysis
Deuschle et al. Robust sensor spike detection method based on dynamic time warping
CN110926784B (en) GIS circuit breaker fault on-line monitoring device based on sound
CN114091368A (en) Method for identifying cavitation state of axial flow turbine
CN111123356B (en) Abnormal track intelligent identification method based on first arrival information
CN110376436B (en) Multi-scale noise power spectral line spectrum detection method
CN102624367A (en) Multi-channel pulse synchronization identification apparatus and method thereof
Hejtmanek et al. Comparison of two denoising techniques to improve UHF partial discharge localization
CN115616082B (en) Keyboard defect analysis method based on noise detection
CN108090432A (en) A kind of automatic identifying method of vibration signals spectrograph spectral peak

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