CN112611687A - Method and system for accurately positioning metal particles in GIL - Google Patents
Method and system for accurately positioning metal particles in GIL Download PDFInfo
- Publication number
- CN112611687A CN112611687A CN202011360057.4A CN202011360057A CN112611687A CN 112611687 A CN112611687 A CN 112611687A CN 202011360057 A CN202011360057 A CN 202011360057A CN 112611687 A CN112611687 A CN 112611687A
- Authority
- CN
- China
- Prior art keywords
- partial discharge
- gil
- matrix
- sensor
- sensors
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000002923 metal particle Substances 0.000 title claims abstract description 53
- 238000000034 method Methods 0.000 title claims abstract description 24
- 239000011159 matrix material Substances 0.000 claims abstract description 74
- 238000005259 measurement Methods 0.000 claims abstract description 13
- 238000001914 filtration Methods 0.000 claims description 5
- 230000001902 propagating effect Effects 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 230000004807 localization Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 238000011900 installation process Methods 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N5/00—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
- Y04S10/52—Outage or fault management, e.g. fault detection or location
Landscapes
- Chemical & Material Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Dispersion Chemistry (AREA)
- Gas-Insulated Switchgears (AREA)
- Testing Relating To Insulation (AREA)
- Locating Faults (AREA)
Abstract
The invention discloses a method for accurately positioning metal particles in GIL, which comprises the steps of obtaining initial partial discharge positions of the metal particles by adopting an extended time difference method according to actual field conditions, and calculating the partial discharge positions of the metal particles by combining multiple sensors to realize the accurate positioning of the metal particles. Firstly, after metal particles are subjected to partial discharge, the arrival time and the corresponding nodes detected by each sensor in the system are obtained, and a partial discharge positioning combination is formed; then, improving an European measurement matrix, and solving a primary position of partial discharge by combining a shortest path algorithm; and finally, calculating the occurrence time of the partial discharge, correcting the occurrence time by combining a GIL topological structure, and fusing multi-sensor information to finish the accurate positioning of the partial discharge. The invention provides a theoretical basis for the accurate positioning of the metal particles in the GIL and has good application value.
Description
Technical Field
The invention belongs to the technical field of metal particle monitoring in GIL, and particularly relates to a method for accurately positioning metal particles in GIL.
Background
GIL (Gas Insulated transmission Line) has been widely used in domestic and foreign power systems because it has the advantages of high voltage, large current, compact structure, flexible arrangement, stable operation, long service life, superior technical index, no external influence, etc.
The modular structural design of the GIL system makes the installation process complex, and once the GIL system fails, the emergency repair workload is large, and the power failure time is long. According to the statistics of the fault cases, insulation faults caused by the metal particles account for a large proportion of GIL faults, so that the method has very important practical significance for accurately monitoring the metal particles in the GIL.
The existing monitoring method is difficult to accurately position the positions of metal particles in the GIL, and the problems that the existing positioning technology has the defects of weak anti-interference capability, low operation stability, low measurement precision and the like when being applied on site, and the positioning is inaccurate or even cannot be performed are easily caused. Therefore, how to precisely locate the metal particles in the GIL is an urgent problem to be solved by those in the related art.
Disclosure of Invention
In order to solve the defects in the prior art, the invention aims to provide a method for accurately positioning metal particles in a GIL, and aims to improve a Euclidean distance measurement matrix, accurately position the accurate positions of the metal particles in the GIL, provide a theoretical basis for a GIL sensor arrangement scheme and accurate positioning, and have strong engineering practicability.
The invention adopts the following technical scheme. A method of accurately locating metal particles in GIL comprising the steps of:
S＝[S1 S2 … Sα … SN] (1)
in the formula:
s denotes a very high frequency sensor arrangement matrix,
Sαdenotes the alpha-th uhf sensor, alpha-1, 2, …, N,
n represents the number of ultrahigh frequency sensors;
and 4, calculating the occurrence time of the partial discharge according to the initial positioning result, obtaining a distance matrix of the partial discharge spread to each ultrahigh frequency sensor, and fusing information of multiple sensors to complete the precise positioning of the partial discharge.
Preferably, step 2 comprises:
step 2.1, after the metal particle partial discharge occurs, k ultrahigh frequency sensors in all N ultrahigh frequency sensors receive electromagnetic waves, k is more than or equal to 1 and less than or equal to N, and t isβRepresenting the propagation time of the electromagnetic wave detected by each uhf sensor matrix, β being 1,2, …, k, an electromagnetic wave propagation time matrix t is formed, which is expressed by the following equation (2),
t＝[t1 t2 … tβ … tk]T (2)
in the formula:
k represents the number of uhf sensors that receive the electromagnetic wave,
t represents a transposed matrix;
step 2.2, acquiring the minimum value of elements in the electromagnetic wave propagation time matrix t, and taking t as tminIndicating that the ultrahigh frequency sensor S corresponding to the minimum value of the arrival time of the electromagnetic wave is foundmAt a very high frequencymAll the nodes with ultrahigh frequency sensors arranged adjacently form an adjacent GIL node matrix SvExpressed by the following formula (3),
Sv＝[Sv1 Sv2 … Svr … Svw] (3)
in the formula:
w denotes an ultrahigh frequency sensor SmThe number of nodes adjacent to which the uhf sensor is disposed,
Svrrepresenting a matrix S of adjacent GIL nodesvAny one of the nodes in the network is provided with a plurality of nodes,
ultrahigh frequency sensor SmAnd adjacent GIL node matrix SvMay constitute w sets of moveout location combinations.
Preferably, step 3 comprises:
step 3.1, to deduct the nearest UHF sensor SmN-1 nodes and w adjacent nodes obtained after the nodes form w x (N-1) groups of lines which can be faulted, namely a time matrix D for the w x (N-1) groups of electromagnetic waves from a partial discharge position to k ultrahigh frequency sensors of the electromagnetic wavesfAs indicated by the general representation of the,
in the formula:
Dij＝[Dij1 Dij2 … Dijk]Trepresenting the theoretical time for the electromagnetic wave to reach each of the k uhf sensors when the ij line segment fails,
k represents the number of sensors receiving electromagnetic waves, and k is more than or equal to 1 and less than or equal to N;
step 3.2, constructing a Euclidean distance measure matrix E,
in the formula:
Eij＝‖Dij-t |, representing the similarity between theoretical and actual time, the smaller the number, the higher the similarity between the two sets of data, denoted as EminRepresenting the minimum value of the Euclidean distance measure matrix E;
step 3.3, obtaining a partial discharge preliminary positioning result from the minimum value in the measure matrix E, representing the corresponding partial discharge position by f, and obtaining the minimum value through the measure matrix, and then obtaining the secondary node SmPosition plus minimum EminF can be obtained.
Preferably, step 4 comprises:
step 4.1, according to the initial partial discharge position f and by combining the inherent structure of the GIL, the occurrence time t of partial discharge can be calculated0The expression is given by the following formula,
in the formula:
tmrepresents the minimum value in the matrix of equation (2),
Lbfindicating propagation of the preliminary partial discharge location f to the node SmCorresponding to the distance of the ultrahigh frequency sensor,
c represents the electromagnetic wave propagation speed;
step 4.2, using t0The distance matrix of the partial discharge signal propagating to each sensor is calculated as,
L＝[L1f L2f L3f … Lkf]T＝c(t-t0·m) (8)
in the formula:
the elements in the distance matrix L are the distances at which the partial discharge signal propagates to the sensors,
step 4.3, combining the topological structure of the GIL system to obtain a plurality of groups of partial discharge position information which is expressed by the following formula,
Sd＝[Sd1 Sd2 Sd3 … Sdk] (9)
step 4.4, fusing the information of each sensor, and finally determining the partial discharge distance,
the actual position of the metal particles can be finally determined according to the partial discharge distance.
Preferably, a filter is provided to perform 400MHz low pass filtering on the electromagnetic wave signals detected by each sensor, the filtered signals containing only TEM waves.
The present invention also provides a system for accurately locating metal particles in a GIL using said method, comprising: a plurality of ultrahigh frequency sensors, a time measuring module, a preliminary positioning module and an accurate positioning module,
a plurality of uhf sensors are arranged at a plurality of nodes in the GIL system, represented by the following formula (1),
S＝[S1 S2 … Sα … SN] (1)
in the formula:
s denotes a very high frequency sensor arrangement matrix,
Sαdenotes the alpha-th uhf sensor, alpha-1, 2, …, N,
n represents the number of uhf sensors and nodes where the uhf sensors are arranged;
the time measurement module is connected with the ultrahigh frequency sensors, and after the metal particle partial discharge occurs, the arrival time and the corresponding node detected by each ultrahigh frequency sensor in the GIL system are obtained;
a preliminary positioning module for ultrahigh frequency sensor S corresponding to minimum electromagnetic wave arrival timemAnd the adjacent sensors form a positioning combination to obtain a partial discharge initial position;
and the accurate positioning module is used for obtaining the initial position of the partial discharge to calculate the occurrence time of the partial discharge, correcting the partial discharge by combining the GIL topological structure and fusing multi-sensor information to finish the accurate positioning of the partial discharge.
Preferably, the time measurement module forms an electromagnetic wave propagation time matrix t after the partial discharge of the metal particles occurs, which is expressed by the following formula (2),
t＝[t1 t2 … tk]T (2)
in the formula:
k represents the number of uhf sensors that receive the electromagnetic wave,
t represents a transposed matrix;
obtaining the minimum value t of the elements in the electromagnetic wave propagation time matrix tminAnd a corresponding ultrahigh frequency sensor Sm。
Preferably, the preliminary localization module forms a Euclidean distance measure matrix E,
in the formula:
Eij＝‖Dij-t‖，Dij＝[Dij1 Dij2 … Dijk]Twhich represents the theoretical time for the electromagnetic wave to reach each of the k uhf sensors when a fault occurs in the ij line section.
Preferably, the precise positioning module averages partial discharge position information of all the uhf sensors receiving the electromagnetic wave signal to obtain precise positioning information.
Preferably, the system further comprises a filter for performing 400MHz low-pass filtering on the electromagnetic wave signals detected by the sensors, and the filtered signals only contain TEM waves.
Compared with the prior art, the method for accurately positioning the metal particles in the GIL has the advantages that the Euclidean distance measurement matrix is improved, the accurate positions of the metal particles in the GIL are accurately positioned, a theoretical basis is provided for the arrangement scheme and the accurate positioning of the GIL sensor, and the method has high engineering practicability.
Drawings
FIG. 1 is a flow chart of one embodiment of the present invention;
FIG. 2 is a schematic view of a proximity sensor according to the present invention;
fig. 3 is a diagram of an in-situ sensor mounting arrangement of the present invention.
Detailed Description
The present application is further described below with reference to the accompanying drawings. The following examples are only for illustrating the technical solutions of the present invention more clearly, and the protection scope of the present application is not limited thereby.
As shown in fig. 1, the present invention provides a method for accurately positioning metal particles in GIL, where accurate positioning in the present invention means a relative error of less than 2%, the method comprising the steps of:
S＝[S1 S2 … Sα … SN] (1)
in the formula:
s denotes a very high frequency sensor arrangement matrix,
Sαdenotes the alpha-th uhf sensor, alpha-1, 2, …, N,
n represents the number of uhf sensors.
It is understood that the number, spacing, and positions of the sensors may be arbitrarily set by those skilled in the art, and the GIL system may include nodes where the vhf sensors are not arranged, in addition to N nodes where the vhf sensors are arranged. For each node, its location information is known. As shown in fig. 3, an exemplary but non-limiting implementation is given by the embodiment of the present invention, where three sensors, a first sensor 1, a second sensor and a third sensor 3, are mounted on the GIL system housing.
And 2, acquiring the arrival time and the corresponding node detected by each ultrahigh frequency sensor in the GIL system after the metal particle partial discharge occurs, and forming a partial discharge positioning combination. The step 2 specifically comprises the following steps:
step 2.1, after the metal particle partial discharge occurs, k ultrahigh frequency sensors in all N ultrahigh frequency sensors receive electromagnetic waves, k is more than or equal to 1 and less than or equal to N, and t isβRepresenting the propagation time of the electromagnetic wave detected by each uhf sensor matrix, β being 1,2, …, k, an electromagnetic wave propagation time matrix t is formed, which is expressed by the following equation (2),
t＝[t1 t2 … tk]T (2)
in the formula:
k represents the number of uhf sensors that receive the electromagnetic wave,
t denotes a transposed matrix.
It is understood that a filter can be set to perform 400MHz low-pass filtering on the Electromagnetic Wave signal detected by each sensor, and the filtered signal only contains TEM waves (the electric vector and the magnetic vector are both perpendicular to the propagation direction).
Step 2.2, acquiring the minimum value of elements in the electromagnetic wave propagation time matrix t, and taking t as tminIndicating that the ultrahigh frequency sensor S corresponding to the minimum value of the arrival time of the electromagnetic wave is foundmAt a very high frequencymAll the nodes with ultrahigh frequency sensors arranged adjacently form an adjacent GIL node matrix SvExpressed by the following formula (3),
Sv＝[Sv1 Sv2 … Svr … Svw] (3)
in the formula:
w denotes an ultrahigh frequency sensor SmThe number of nodes adjacent to which the uhf sensor is disposed,
Svrrepresenting a matrix S of adjacent GIL nodesvAny one of the nodes, as shown in fig. 2.
Ultrahigh frequency sensor S when partial discharge of metal particles occurs in GIL systemmAnd adjacent GIL node matrix SvThe ultrahigh frequency sensors in the system can form w groups of time difference positioning combinations and are popularized to a GIL ultrahigh frequency sensor arrangement matrix S and the ultrahigh frequency sensors SmAnd N-1 time difference positioning combinations can be obtained with other N-1 ultrahigh frequency sensors. It is understood that in the following calculation, the matrix S is already usedvNode S and nodemThe formed positioning combination is unfolded.
And 3, constructing an European measurement matrix, and solving a primary position of the partial discharge by combining a shortest path algorithm. The step 3 specifically comprises the following steps:
and 3.1, in the GIL system, ij represents a line with a fault, and if the node is a node without a sensor, the signal can be analyzed after being transmitted to the nodes adjacent to the sensor. Therefore, in the analysis, the nearest ultrahigh frequency sensor S can be deductedmN-1 nodes and w adjacent nodes obtained after the nodes form w x (N-1) groups of lines which can be in fault, namely w x (N-1) groups of electromagnetic waves from a partial discharge position to k special nodes receiving the electromagnetic wavesTime matrix D of high-frequency sensorfAs indicated by the general representation of the,
in the formula:
Dij＝[Dij1 Dij2 … Dijk]Trepresenting the theoretical time for the electromagnetic wave to reach each of the k uhf sensors when the ij line segment fails,
k represents the number of sensors receiving electromagnetic waves, and k is more than or equal to 1 and less than or equal to N.
It is understood that, when calculating the theoretical time, the theoretical time is taken as a theoretical fault occurrence point by the midpoint of the line ij between the node i and the node j, and the theoretical time is equal to the distance from the fault occurrence point to the sensor divided by the propagation velocity of the electromagnetic wave.
Step 3.2, constructing a Euclidean distance measure matrix E,
in the formula:
Eij＝‖Dij-t |, representing the similarity between theoretical and actual time, the smaller the number, the higher the similarity between the two sets of data, denoted as EminRepresents the minimum value of the euclidean distance measure matrix E. Obtained EminAnd the subscript ij of (a) is the line where the fault point is located. EminIs a dimensionless value.
Step 3.3, from the minimum E in the measure matrix EminObtaining the initial positioning result of partial discharge, wherein the corresponding partial discharge position is represented by f, namely after the minimum value is obtained through the measure matrix, the node SmPlus the minimum value EminF can be obtained.
And 4, calculating the occurrence time of the partial discharge, correcting the occurrence time by combining a GIL topological structure, and fusing multi-sensor information to finish the accurate positioning of the partial discharge.
Step 4.1, according to the initial partial discharge position f and by combining the inherent structure of the GIL, the occurrence time t of partial discharge can be calculated0The expression is given by the following formula,
in the formula:
tminrepresents the minimum value in the matrix of equation (2),
Lbfindicating propagation of the preliminary partial discharge location f to the node SmCorresponding to the distance of the ultrahigh frequency sensor,
c represents the electromagnetic wave propagation velocity.
Step 4.2, using t0The distance matrix L of the partial discharge signal propagating to each sensor is calculated as,
L＝[L1f L2f L3f … Lkf]T＝c(t-t0·m) (8)
in the formula:
the elements in the distance matrix L are the distances at which the partial discharge signal propagates to the sensors,
step 4.3, the distance from the partial discharge signal to each sensor is calculated by the formula (8), the position of each sensor is known, the position information of the partial discharge signal is obtained by adding the distance to the position of the sensor, and a plurality of groups of partial discharge position information are obtained by combining the topological structure of the GIL system and are expressed by the following formula,
Sd＝[Sd1 Sd2 Sd3 … Sdk] (9)
step 4.4, fusing the information of each sensor, and finally determining the partial discharge distance,
the actual position of the metal particles can be finally determined according to the partial discharge distance.
The present invention also provides a system for accurately locating metal particles in a GIL using said method, comprising: a plurality of uhf sensors, a time measuring module, a preliminary location module and a fine location module, the plurality of uhf sensors being disposed at a plurality of nodes in the GIL system, expressed by the following formula (1),
S＝[S1 S2 … Sα … SN] (1)
in the formula:
s denotes a very high frequency sensor arrangement matrix,
Sαdenotes the alpha-th uhf sensor, alpha-1, 2, …, N,
n represents the number of uhf sensors and nodes where the uhf sensors are arranged;
the time measurement module is connected with the ultrahigh frequency sensors, and after the metal particle partial discharge occurs, the arrival time and the corresponding node detected by each ultrahigh frequency sensor in the GIL system are obtained;
a preliminary positioning module for ultrahigh frequency sensor S corresponding to minimum electromagnetic wave arrival timemAnd the adjacent sensors form a positioning combination to obtain a partial discharge initial position;
and the accurate positioning module is used for calculating the occurrence time of the partial discharge according to the obtained initial position of the partial discharge, correcting the partial discharge by combining the GIL topological structure and fusing multi-sensor information to finish accurate positioning of the partial discharge.
Application example:
as shown in fig. 3, an exemplary but non-limiting implementation is given by the embodiment of the present invention, with three sensors mounted on the GIL housing. The distance between the first sensor 1 and the second sensor 2 is 1.3 meters, and the distance between the second sensor 2 and the third sensor 3 is 3.2 meters.
For the example shown in fig. 3, the amplitude and waveform of the signal detected by each sensor are small, so as to obtain the time matrix T of the electromagnetic wave generated by partial discharge reaching the ultrahigh frequency sensor of the system, as shown in table 1.
TABLE 1 arrival time of electromagnetic waves at each sensor
| Sensor numbering | 1 | 2 | 3 |
| Time t/ns | 19.4 | 11.6 | 15.8 |
For the example shown in fig. 3, analyzing the sensor measurement data, it can be seen that the time for the electromagnetic wave to reach sensor No. 2 is the shortest. And sequentially carrying out time difference positioning miscalculation on the No. 2 sensor and the rest 2 ultrahigh frequency sensors, wherein the partial discharge positioning data is shown in a table 2.
TABLE 2 partial location discharge data
| Sensor combination | [2，1] | [2，3] |
| Euclidean distance measure matrix E | 0.0017 | 2.1517 |
| Partial discharge section | [1，2] | [2，3] |
| Partial discharge distance f/m | 2.770 | 1.280 |
It will be appreciated that table 2 is based on sensor segmentation, for example, where the electromagnetic wave arrives at sensor 1 and sensor 2 at different times for the [2,1] segment, 1 f is derived from the above steps, and where the arrival times at 2 and 3 are different for the [2,3] segment, then the other 1 f is derived. The corrected partial discharge distance is obtained from the equations (8) and (9), as shown in table 3.
TABLE 3 partial localization discharge data
| Sensor numbering | 1 | 2 | 3 |
| Partial discharge distance | 2.7323 | 2.7362 | 2.7358 |
| Absolute error | 0.0323 | 0.0362 | 0.0358 |
It is calculated according to equation (10) that the partial discharge occurs on segment 1-2, at a distance of sensor 2.7348m No. 1.
Through field exploration, the actual local discharge position is on the section 1-2, 2.7m away from the No. 1 sensor, the absolute error is 0.0348m, the relative error is 1.09%, the error requirement of field positioning is met, and the method has good engineering application value.
Compared with the prior art, the method for accurately positioning the metal particles in the GIL has the advantages that the Euclidean distance measurement matrix is improved, the accurate positions of the metal particles in the GIL are accurately positioned, a theoretical basis is provided for the arrangement scheme and the accurate positioning of the GIL sensor, and the method has high engineering practicability.
The present applicant has described and illustrated embodiments of the present invention in detail with reference to the accompanying drawings, but it should be understood by those skilled in the art that the above embodiments are merely preferred embodiments of the present invention, and the detailed description is only for the purpose of helping the reader to better understand the spirit of the present invention, and not for limiting the scope of the present invention, and on the contrary, any improvement or modification made based on the spirit of the present invention should fall within the scope of the present invention.
Claims (10)
1. A method of accurately locating metal particles in GIL comprising the steps of:
step 1, arranging N ultrahigh frequency sensors at N nodes in a GIL system, which is expressed by the following formula (1),
S＝[S1 S2 … Sα … SN] (1)
in the formula:
s denotes a very high frequency sensor arrangement matrix,
Sαdenotes the alpha-th uhf sensor, alpha-1, 2, …, N,
n represents the number of ultrahigh frequency sensors;
step 2, after the metal particle partial discharge occurs, acquiring the arrival time and the corresponding nodes detected by each ultrahigh frequency sensor in the GIL system, and forming a metal particle partial discharge positioning combination by the ultrahigh frequency sensor which firstly receives a partial discharge electromagnetic wave signal and each adjacent node thereof; with SmA sensor for indicating the first received partial discharge electromagnetic wave signal;
step 3, forming an Euclidean distance measure matrix according to the difference between the theoretical time and the actual time when the metal particles in each section of the line are partially discharged and reach each node, and using the minimum value in the Euclidean distance measure matrix and the sensor SmAdding the positions to obtain a partial discharge preliminary positioning result;
and 4, calculating the occurrence time of the partial discharge according to the initial positioning result, obtaining a distance matrix of the partial discharge spread to each ultrahigh frequency sensor, and fusing information of multiple sensors to complete the precise positioning of the partial discharge.
2. The method of accurately positioning metal particles in a GIL according to claim 1, characterized in that:
the step 2 comprises the following steps:
step 2.1, after the metal particle partial discharge occurs, k ultrahigh frequency sensors in all N ultrahigh frequency sensors receive electromagnetic waves, k is more than or equal to 1 and less than or equal to N, and t isβRepresenting the propagation time of the electromagnetic wave detected by each uhf sensor matrix, β being 1,2, …, k, an electromagnetic wave propagation time matrix t is formed, which is expressed by the following equation (2),
t＝[t1 t2 … tβ … tk]T (2)
in the formula:
k represents the number of uhf sensors that receive the electromagnetic wave,
t represents a transposed matrix;
step 2.2, acquiring the minimum value of elements in the electromagnetic wave propagation time matrix t, and taking t as tminIndicating that the ultrahigh frequency sensor S corresponding to the minimum value of the arrival time of the electromagnetic wave is foundmAt a very high frequencymAll the nodes with ultrahigh frequency sensors arranged adjacently form an adjacent GIL node matrix SvExpressed by the following formula (3),
Sv＝[Sv1 Sv2 … Svr … Svw] (3)
in the formula:
w denotes an ultrahigh frequency sensor SmThe number of nodes adjacent to which the uhf sensor is disposed,
Svrrepresenting a matrix S of adjacent GIL nodesvAny one of the nodes in the network is provided with a plurality of nodes,
ultrahigh frequency sensor SmAnd adjacent GIL node matrix SvMay constitute w sets of moveout location combinations.
3. The method of accurately positioning metal particles in a GIL according to claim 2, characterized in that:
the step 3 comprises the following steps:
step 3.1, to deduct the nearest UHF sensor SmN-1 nodes and w adjacent nodes obtained after the nodes form w x (N-1) groups of lines which can be faulted, namely a time matrix D for the w x (N-1) groups of electromagnetic waves from a partial discharge position to k ultrahigh frequency sensors of the electromagnetic wavesfAs indicated by the general representation of the,
in the formula:
Dij＝[Dij1 Dij2 … Dijk]Trepresenting the theoretical time for the electromagnetic wave to reach each of the k uhf sensors when the ij line segment fails,
k represents the number of sensors receiving electromagnetic waves, and k is more than or equal to 1 and less than or equal to N;
step 3.2, constructing a Euclidean distance measure matrix E,
in the formula:
Eij＝||Dij-t | |, which represents the similarity between theoretical time and actual time, the smaller the numerical value, the higher the similarity of two groups of data, denoted as EminRepresenting the minimum value of the Euclidean distance measure matrix E;
step 3.3, obtaining a partial discharge preliminary positioning result from the minimum value in the measure matrix E, representing the corresponding partial discharge position by f, and obtaining the minimum value through the measure matrix, and then obtaining the secondary node SmPosition plus minimum EminF can be obtained.
4. The method of accurately positioning metal particles in a GIL according to claim 3, characterized in that:
step 4 comprises the following steps:
step 4.1, according to the initial partial discharge position f, combining the inherent structure of the GIL to calculate the occurrence time t of partial discharge0The expression is given by the following formula,
in the formula:
tmrepresents the minimum value in the matrix of equation (2),
Lbfindicating propagation of the preliminary partial discharge location f to the node SmCorresponding to the distance of the ultrahigh frequency sensor,
c represents the electromagnetic wave propagation speed;
step 4.2, using t0The distance matrix of the partial discharge signal propagating to each sensor is calculated as,
L＝[L1f L2f L3f … Lkf]T＝c(t-t0·m) (8)
in the formula:
the elements in the distance matrix L are the distances at which the partial discharge signal propagates to the sensors,
step 4.3, combining the topological structure of the GIL system to obtain a plurality of groups of partial discharge position information which is expressed by the following formula,
Sd＝[Sd1 Sd2 Sd3 … Sdk] (9)
step 4.4, fusing the information of each sensor, and finally determining the partial discharge distance,
and finally determining the actual position of the metal particles according to the partial discharge distance.
5. The method of accurately positioning metal particles in a GIL according to claim 3, characterized in that:
the filters are arranged to carry out 400MHz low-pass filtering on the electromagnetic wave signals detected by the sensors, and the filtered signals only contain TEM waves.
6. A system for accurately locating metal particles in a GIL using the method of any one of claims 1-5, comprising: a plurality of superfrequency sensors, time measurement module, preliminary positioning module and accurate positioning module, its characterized in that:
a plurality of uhf sensors are arranged at a plurality of nodes in the GIL system, represented by the following formula (1),
S＝[S1 S2 … Sα … SN] (1)
in the formula:
s denotes a very high frequency sensor arrangement matrix,
Sαdenotes the alpha-th uhf sensor, alpha-1, 2, …, N,
n represents the number of uhf sensors and nodes where the uhf sensors are arranged;
the time measurement module is connected with the ultrahigh frequency sensors, and after the metal particle partial discharge occurs, the arrival time and the corresponding node detected by each ultrahigh frequency sensor in the GIL system are obtained;
a preliminary positioning module for ultrahigh frequency sensor S corresponding to minimum electromagnetic wave arrival timemAnd the adjacent sensors form a positioning combination to obtain a partial discharge initial position;
and the accurate positioning module is used for calculating the occurrence time of the partial discharge according to the obtained initial position of the partial discharge, correcting the partial discharge by combining the GIL topological structure and fusing multi-sensor information to finish accurate positioning of the partial discharge.
7. The system for accurately positioning metal particles in a GIL according to claim 6, wherein:
the time measurement module forms an electromagnetic wave propagation time matrix t after the metal particle partial discharge occurs, which is expressed by the following formula (2),
t＝[t1 t2 … tk]T (2)
in the formula:
k represents the number of uhf sensors that receive the electromagnetic wave,
t represents a transposed matrix;
obtaining the minimum value t of the elements in the electromagnetic wave propagation time matrix tminAnd a corresponding ultrahigh frequency sensor Sm。
8. The system for accurately positioning metal particles in a GIL according to claim 7, wherein:
the preliminary positioning module forms a euclidean distance measure matrix E,
in the formula:
Eij＝||Dij-t||，Dij＝[Dij1 Dij2 … Dijk]Twhich represents the theoretical time for the electromagnetic wave to reach each of the k uhf sensors when a fault occurs in the ij line section.
9. The system for accurately positioning metal particles in a GIL according to claim 8, wherein:
the accurate positioning module averages partial discharge position information of all the ultrahigh frequency sensors receiving electromagnetic wave signals to obtain accurate positioning information.
10. The system for accurately positioning metal particles in a GIL according to claim 8, wherein:
the system also comprises a filter, which is used for performing 400MHz low-pass filtering on the electromagnetic wave signals detected by each sensor, and the filtered signals only contain TEM waves.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011360057.4A CN112611687B (en) | 2020-11-27 | 2020-11-27 | Method and system for accurately positioning metal particles in GIL |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011360057.4A CN112611687B (en) | 2020-11-27 | 2020-11-27 | Method and system for accurately positioning metal particles in GIL |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112611687A true CN112611687A (en) | 2021-04-06 |
| CN112611687B CN112611687B (en) | 2023-05-30 |
Family
ID=75228056
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011360057.4A Active CN112611687B (en) | 2020-11-27 | 2020-11-27 | Method and system for accurately positioning metal particles in GIL |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112611687B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114217189A (en) * | 2021-12-30 | 2022-03-22 | 国网江苏省电力有限公司南通供电分公司 | GIL equipment fault positioning method adopting ultrahigh frequency transient current measurement |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA1154860A (en) * | 1980-07-10 | 1983-10-04 | Lev B. Kublanovsky | Method of locating a fault in a pressure pipe system from the coincidence of at least two sets of possible points |
| CN107037339A (en) * | 2017-04-28 | 2017-08-11 | 中国矿业大学 | A kind of uhf sensor points distributing method for monitoring GIS partial discharge |
| CN108169634A (en) * | 2017-12-19 | 2018-06-15 | 上海交通大学 | A kind of accurate method for obtaining the partial-discharge ultrahigh-frequency signal time difference |
| CN109782139A (en) * | 2018-12-25 | 2019-05-21 | 西安交通大学 | A kind of GIS ultrahigh frequency partial discharge monitoring system and its monitoring method |
| CN109932627A (en) * | 2019-04-30 | 2019-06-25 | 国网山东省电力公司烟台供电公司 | A kind of localization method and system of GIS pressure test shelf depreciation |
-
2020
- 2020-11-27 CN CN202011360057.4A patent/CN112611687B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA1154860A (en) * | 1980-07-10 | 1983-10-04 | Lev B. Kublanovsky | Method of locating a fault in a pressure pipe system from the coincidence of at least two sets of possible points |
| CN107037339A (en) * | 2017-04-28 | 2017-08-11 | 中国矿业大学 | A kind of uhf sensor points distributing method for monitoring GIS partial discharge |
| CN108169634A (en) * | 2017-12-19 | 2018-06-15 | 上海交通大学 | A kind of accurate method for obtaining the partial-discharge ultrahigh-frequency signal time difference |
| CN109782139A (en) * | 2018-12-25 | 2019-05-21 | 西安交通大学 | A kind of GIS ultrahigh frequency partial discharge monitoring system and its monitoring method |
| CN109932627A (en) * | 2019-04-30 | 2019-06-25 | 国网山东省电力公司烟台供电公司 | A kind of localization method and system of GIS pressure test shelf depreciation |
Non-Patent Citations (3)
| Title |
|---|
| 丁然: "D1 GIS局部放电检测机理与精确定位研究", 《CNKI博士电子期刊》 * |
| 周会峰等: "输电线路故障行波网络定位新方法", 《电力系统自动化》 * |
| 戴锋等: "考虑短线路折反射的配电网故障定位研究", 《电力系统保护与控制》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114217189A (en) * | 2021-12-30 | 2022-03-22 | 国网江苏省电力有限公司南通供电分公司 | GIL equipment fault positioning method adopting ultrahigh frequency transient current measurement |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112611687B (en) | 2023-05-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108375713B (en) | Novel power grid fault traveling wave positioning method and system | |
| CN109490706B (en) | Multi-branch power transmission line fault positioning method | |
| WO2019232595A1 (en) | A method of estimating the location of a fault on an electrical distribution network and an associated system | |
| Salehi et al. | Fault location on branched networks using mathematical morphology | |
| CN113721110B (en) | Regional power grid fault single/double-end positioning method based on node fault domain time difference method | |
| CN104330708A (en) | Fault location method for wide area traveling wave signal below power grid | |
| CN106771861A (en) | Complex electric network Fault Locating Method based on wide area traveling wave energy and time difference | |
| US20240044965A1 (en) | Parameter independent traveling wave-based fault location using unsynchronized measurements | |
| CN117192292B (en) | Lightning grounding electrode line fault distance measurement method and system | |
| CN112611687B (en) | Method and system for accurately positioning metal particles in GIL | |
| CN116930685B (en) | Single-end ranging method suitable for single-phase earth fault of power distribution network | |
| CN115097253A (en) | MMC-HVDC direct current transmission line fault distance measurement method and system | |
| CN113253052A (en) | High-voltage direct-current transmission line fault distance measurement method based on improved SMMG | |
| CN102540015B (en) | Single-phase earthing of small current earthing system distance-finding method and distance measuring equipment | |
| CN110058127A (en) | Non-effective earthing one-phase earthing failure in electric distribution network point recognition methods | |
| Kulikov et al. | Methods for implementing of the single-ended traveling wave fault locating on the transmission lines based on traveling wave pattern recognition | |
| CN111812456B (en) | Power transmission line fault positioning method, device and system | |
| CN115308524A (en) | Distribution line traveling wave positioning method based on zero mode and line mode traveling wave fusion | |
| Tian et al. | Accurate fault location of hybrid lines in distribution networks | |
| CN114755530A (en) | Robust fault positioning method for power transmission line | |
| CN108599115B (en) | Wide-area protection method and device for power transmission line | |
| CN111596226A (en) | Power distribution network ground fault positioning method and system | |
| CN110726904B (en) | Method for accurately positioning fault point of distributed T-connection power transmission circuit | |
| Chen et al. | Travelling Wave Single End Fault Location Method based on Network Information | |
| CN116609618B (en) | Method for achieving fault location based on primary and secondary fusion of breaker to obtain traveling wave |
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 |