CN112834849A - Ultrasonic positioning method and device for partial discharge source of transformer - Google Patents
Ultrasonic positioning method and device for partial discharge source of transformer Download PDFInfo
- Publication number
- CN112834849A CN112834849A CN202110006160.7A CN202110006160A CN112834849A CN 112834849 A CN112834849 A CN 112834849A CN 202110006160 A CN202110006160 A CN 202110006160A CN 112834849 A CN112834849 A CN 112834849A
- Authority
- CN
- China
- Prior art keywords
- ultrasonic
- transformer
- node
- nodes
- estimated
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
-
- 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/1209—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 using acoustic measurements
-
- 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/18—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
- G01S5/22—Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Testing Relating To Insulation (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
The invention discloses an ultrasonic positioning method and device for a transformer local discharge source, which are used for solving the problem that when the ultrasonic positioning of the local discharge source is carried out, the positioning has larger errors due to the influence of metal parts in a transformer. According to the transformer structure, a transformer node numerical model is constructed; receiving actual arrival time differences measured by a plurality of ultrasonic sensors arranged on a transformer; traversing a plurality of nodes in the transformer node numerical model, and determining estimated arrival time differences of ultrasonic signals respectively arriving at a plurality of ultrasonic sensors from the plurality of nodes through an ultrasonic path search algorithm based on an A-x routing algorithm; and determining the position of the partial discharge source from the plurality of nodes through a mixed frog-leaping algorithm according to the actual time difference of arrival and the estimated time difference of arrival. The invention can improve the positioning precision of the local discharge source, avoid the interference of internal metal parts and improve the positioning accuracy of the local discharge source.
Description
Technical Field
The invention relates to the technical field of electricity, in particular to an ultrasonic positioning method and device for a partial discharge source of a transformer.
Background
The power transformer is an important device in the power system, and the insulation state of the power transformer directly influences the safe and stable operation of the power system. In the transformer, there may be a phenomenon of partial discharge. Such non-penetrating discharge occurring in the equipment insulation is a cause of deterioration of the insulation of the electric power equipment. Therefore, in the initial stage of discharge, the local discharge source is detected and positioned, and the method is very important for operation maintenance and insulation failure early warning of equipment.
The generation of partial discharge is accompanied by a series of physicochemical phenomena such as generation of ultrasonic waves, electromagnetic radiation, electric pulses, luminescence, heat generation, new products, and the like.
Currently, when a partial discharge source in a transformer is positioned, an insulation sensor is usually arranged, and an ultrasonic method is used for detection. However, in actual situations, metal components such as an iron core and a winding exist in the transformer, and when the ultrasonic wave propagates along a straight path and passes through the metal components, a large amount of energy is attenuated due to refraction and reflection, so that an ultrasonic wave signal cannot be received by the insulation sensor, and positioning errors of the local discharge source are caused.
Disclosure of Invention
The invention provides an ultrasonic positioning method and device for a transformer local discharge source, which are used for solving the problem that when the ultrasonic positioning of the local discharge source is carried out, the positioning has larger errors due to the influence of metal parts in a transformer.
Therefore, the invention adopts the following technical scheme: an ultrasonic positioning method of a transformer partial discharge source comprises the following steps:
constructing a transformer node numerical model according to the transformer structure;
receiving actual arrival time differences measured by a plurality of ultrasonic sensors arranged on the transformer, wherein the actual arrival time differences are actual arrival time differences of ultrasonic signals generated by a local discharge source and arriving at the plurality of ultrasonic sensors;
traversing a plurality of nodes in the transformer node numerical model, and determining estimated arrival time differences of the ultrasonic signals from the plurality of nodes to the plurality of ultrasonic sensors respectively through an ultrasonic path search algorithm based on an A-x routing algorithm;
and determining the position of the partial discharge source from the plurality of nodes by a mixed frog-leaping algorithm according to the actual time difference of arrival and the estimated time difference of arrival.
In an implementation manner of the present invention, constructing a transformer node numerical model according to a transformer structure specifically includes: determining nodes corresponding to different structures in the constructed corresponding transformer node numerical model and position parameters, price parameters and speed parameters corresponding to the nodes respectively according to different structures of different transformers; the node numerical model corresponds to the transformer, and each node in the node numerical model corresponds to a corresponding position in the transformer respectively; the position parameter represents the position of a node, the price parameter represents the difficulty of the ultrasonic signal propagating in the node, the speed parameter represents the propagation speed of the ultrasonic signal in the node, and the price parameter and the speed parameter are related to the structure of the transformer.
In one implementation of the present invention, the a-routing algorithm includes an estimated price distance H value of a node, where the estimated price distance H value represents an estimated distance from the node to an end point of a corresponding shortest propagation path, and the estimated price distance H value is related to the price parameter; the estimated price distance H value is determined by: determining a plurality of corresponding estimated price distances between the nodes and a terminal point according to the directions of a plurality of nodes adjacent to the nodes in the transformer node numerical model; and determining the minimum estimated price distance from the plurality of estimated price distances as the H value of the node.
In one implementation of the invention, the estimated price distance H is given by Hi＝min(Hi1,Hi2,Hi3,Hi4,Hi5,Hi6) It is determined that,
wherein the content of the first and second substances,
Hi1,Hi2,Hi3,Hi4,Hi5,Hi6representing the estimated price distance, x, corresponding to the six directions of the nodei、yi、ziRepresenting coordinates of nodes, m representing the number of nodes corresponding to the x-axis direction, n representing the number of nodes corresponding to the y-axis direction, l representing the number of nodes corresponding to the z-axis direction, and xend、yend、zendCoordinates representing the end point, poilAnd the price parameter of the node corresponding to the transformer oil in the transformer structure is represented.
In an implementation manner of the present invention, determining, by an ultrasonic path search algorithm based on the a-x routing algorithm, estimated arrival time differences of the ultrasonic signals respectively arriving at the ultrasonic sensors from the nodes specifically includes: determining the shortest propagation paths of the ultrasonic signals from the nodes to the corresponding ultrasonic sensors respectively through an ultrasonic path search algorithm based on an A-x routing algorithm; calculating to obtain the propagation time corresponding to each shortest propagation path according to the determined shortest propagation path; and determining estimated arrival time differences of the ultrasonic signals respectively arriving at the ultrasonic sensors from the nodes according to the calculated propagation time.
In an implementation manner of the present invention, calculating, according to the determined shortest propagation path, a propagation time corresponding to each shortest propagation path includes: according to
Calculating the propagation time; where t denotes the propagation time, x1、y1、z1Coordinates, x, representing the starting point of the shortest propagation pathj、yj、zjCoordinates, v, representing nodes in the shortest propagation pathjAnd expressing the speed parameters of the nodes, num expresses the number of all nodes in the shortest propagation path, and dl expresses the corresponding side length of the nodes in the transformer node data model.
In an implementation manner of the present invention, determining the position of the partial discharge source from the plurality of nodes by a mixed frog-leaping algorithm according to the actual time difference of arrival and the estimated time difference of arrival specifically includes: calculating the fitness value corresponding to each node respectively through a mixed frog-leaping algorithm and the difference between the actual arrival time difference and the estimated arrival time difference; iteration and swarm optimization are carried out through a mixed frog-leaping algorithm, the fitness value corresponding to each node is recalculated, and the node with the highest fitness value is determined from the plurality of nodes and serves as the position of the local discharge source.
In one implementation manner of the invention, the number of the ultrasonic sensors is several; the calculating the fitness value corresponding to each node respectively through a mixed frog-leaping algorithm and the difference between the actual arrival time difference and the estimated arrival time difference specifically comprises the following steps: by passing
Calculating the fitness value corresponding to each node; wherein, PiA frog representing the calculated fitness value;
representing the actual time difference of arrival of the ultrasonic signals received by the first and second sensors,
indicating that the source of partial discharge is located at PiWhen the position of the ultrasonic sensor is detected, the estimated wavelength time difference of the ultrasonic signals received by the first sensor and the second sensor is calculated;
representing the actual time difference of arrival of the ultrasonic signals received by the first and third sensors,
indicating that the source of partial discharge is located at PiWhen the position of the sensor is within the preset range, the estimated wavelength time difference of the ultrasonic signals received by the first sensor and the third sensor is calculated;
representing the actual time difference of arrival of the ultrasonic signals received by the first and fourth sensors,
indicating that the source of partial discharge is located at PiThe estimated wavelength time difference of the received ultrasonic signals of the first sensor and the fourth sensor is calculated.
In one implementation of the invention, the transformer comprises transformer oil, internal metal components, a metal tank wall; determining the shortest propagation path of the ultrasonic signal from the plurality of nodes to the corresponding ultrasonic sensor respectively through an ultrasonic path search algorithm based on an a-x routing algorithm, specifically comprising: and determining the shortest propagation path of the ultrasonic signals from the plurality of nodes to the corresponding ultrasonic sensor respectively by an ultrasonic path search algorithm based on the A-th routing algorithm, wherein the shortest propagation path does not include the node corresponding to the internal metal component.
The invention also provides an ultrasonic positioning device of the transformer partial discharge source, which comprises:
the building module is used for building a transformer node numerical model according to the transformer structure;
the receiving module is used for receiving actual arrival time differences measured by a plurality of ultrasonic sensors arranged on the transformer, wherein the arrival time differences are actual arrival time differences of ultrasonic signals generated by a local discharge source and arriving at the ultrasonic sensors;
the path determining module is used for traversing a plurality of nodes in the transformer node numerical model and determining estimated arrival time differences of the ultrasonic signals from the nodes to the ultrasonic sensors respectively through an ultrasonic path searching algorithm based on an A-x routing algorithm;
and the position determining module is used for determining the position of the partial discharge source from the plurality of nodes through a mixed frog-leaping algorithm according to the actual time difference of arrival and the estimated time difference of arrival.
The invention has the following beneficial effects:
according to the method, the propagation path of the ultrasonic signal is considered under the blocking effect of metal components in the transformer on the ultrasonic signal through the A-path searching algorithm, and the estimated arrival time differences corresponding to all positions of the transformer oil are optimized through the mixed frog-leaping algorithm, so that the position where the local discharge source is most likely to exist is determined, and the positioning of the local discharge source is realized. Each position of the transformer oil is accurately divided in a numerical model building mode, so that the accuracy of the positioning of the partial discharge source is improved, the interference of metal parts in the transformer is avoided, iteration optimization is combined, the accuracy of the positioning of the partial discharge source is improved, and the method is high in feasibility and easy to realize.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the invention and not to limit the invention. In the drawings:
fig. 1 is a flowchart of an ultrasonic positioning method for a partial discharge source of a transformer according to an embodiment of the present invention;
FIG. 2 is a flow chart of another method for ultrasonic localization of a transformer partial discharge source according to an embodiment of the present invention;
fig. 3 is a schematic structural diagram of an ultrasonic positioning device of a transformer partial discharge source according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the technical solutions of the present invention will be clearly and completely described below with reference to the specific embodiments of the present invention and the accompanying drawings. It is to be understood that the described embodiments are merely exemplary of the invention, and not restrictive of the full scope of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Fig. 1 is a flowchart of an ultrasonic positioning method for a transformer partial discharge source according to an embodiment of the present invention, which specifically includes the following steps:
s101: and constructing a transformer node numerical model according to the transformer structure.
In the embodiment of the invention, for convenience of calculation, the structure of the transformer to be detected can be simplified, and the corresponding transformer node numerical model is constructed according to the simplified transformer structure. The method is applicable to transformers with different structures.
Specifically, a transformer node numerical model can be established through gridding, and the solid structure of the transformer is converted into nodes in the transformer node numerical model, wherein each node corresponds to the solid structure at a corresponding position in the transformer.
In one possible implementation, each node may be represented as a cube with a certain side length, and the transformer node numerical model may be divided into a plurality of regularly arranged cube nodes. Wherein the side length dl of the node may be 0.01 meter.
Each node may include several attributes including a location parameter, a price parameter, a speed parameter, and the like. The position parameter indicates the position of the node, and may be represented by (x, y, z), where x is 0,1,2, … …, m, y is 0,1,2, … …, n, z is 0,1,2, … …, and l. It can be known that the transformer node numerical model includes m nodes at most in the x-axis direction, n nodes at most in the y-axis direction, and l nodes at most in the z-axis direction. The price parameter represents the difficulty of the ultrasonic signal in the node transmission, the price parameter of the corresponding node is different according to different materials of the transformer structure, and the larger the price parameter is, the greater the difficulty of the ultrasonic signal in the corresponding node transmission is represented. The speed parameter represents the propagation speed of the ultrasonic signal in the node, and the speed parameter of the corresponding node is different according to different materials of the transformer structure. The price parameter is related to the speed parameter, which is in a negative correlation relationship.
In one embodiment, the transformer structure may include transformer oil, internal metal components, and a metal tank wall. Since the ultrasonic signal is hard to propagate through the internal metal member, the price parameter of the node corresponding to the internal metal member is set to 1000, and the speed parameter is set to 0.1 m/s. According to the propagation condition of the ultrasonic signal in the metal box wall, the price parameter of the metal box wall can be set to be 1, and the speed parameter can be expressed as vsteelParticularly 6100m/s is preferable. The price parameter of the corresponding node of the transformer oil can be expressed as poil，
voil1390m/s can be specifically selected for the speed parameter corresponding to the transformer oil node, when v issteel and voilWhen fixed, canDetermination of poilThe numerical value of (c).
S102: actual arrival time differences measured by a number of ultrasonic sensors provided on the transformer are received.
A number of ultrasonic sensors may be provided on the transformer for detecting ultrasonic signals generated by a local discharge source within the transformer. Accordingly, the arrival time difference of the ultrasonic signal at each ultrasonic sensor can be calculated from the time at which each ultrasonic sensor detects the ultrasonic signal.
It can be determined that the time difference of arrival is actual data corresponding to the actually measured ultrasonic signal generated by the local discharge source.
It should be noted that the ultrasonic sensors are arranged outside the metal box wall of the transformer, the number and the positions of the ultrasonic sensors are not limited in the invention, but in practical application, more than 3 ultrasonic sensors are arranged, and the ultrasonic sensors are arranged in different directions of the transformer, which is beneficial to accurately positioning the local discharge source.
S103: and traversing a plurality of nodes in the transformer node numerical model, and determining estimated arrival time differences of ultrasonic signals respectively arriving at the ultrasonic sensors from the nodes through an ultrasonic path search algorithm based on an A-x routing algorithm.
In the embodiment of the present invention, in order to realize the positioning of the partial discharge source, the actual arrival time difference measured in S102 needs to be simulated to determine the most likely position of the partial discharge source corresponding to the arrival time difference. That is, when it is necessary to determine where the partial discharge source is located in the transformer, the actual arrival time difference data is generated.
Based on the method, a plurality of nodes in the transformer node numerical model can be traversed, and the shortest propagation path and the propagation time of the corresponding ultrasonic signal to each ultrasonic sensor and the corresponding arrival time difference are respectively determined when each node is used as a local discharge source. Since the partial discharge source is only generated in the transformer oil, the nodes represent the nodes corresponding to the transformer oil, and the arrival time difference is the estimation data.
Specifically, the process of determining the shortest propagation path of the ultrasonic signal from the node to the corresponding ultrasonic sensor through the ultrasonic path search algorithm based on the a-x routing algorithm mainly comprises the following steps;
firstly, for each node corresponding to the transformer oil, the shortest propagation path from the position of the ultrasonic wave to each ultrasonic sensor is determined. Thus, in calculating each of the shortest propagation paths separately, one of the plurality of nodes may be determined as a starting point of the shortest propagation path, and one of the plurality of ultrasonic sensors may be determined as an ending point of the shortest propagation path.
Secondly, determining the coordinates corresponding to the starting point, setting the corresponding F value to be 0, and adding the F value into an opening list of the A-path searching algorithm. The nodes in the opening list represent reachable nodes, the value F represents the sum of the price distance G value of the node from the starting point and the estimated price distance H value from the end point, and the price distance G value and the estimated price distance H value are related to the price parameters of the nodes.
And thirdly, determining the node with the lowest F value in the open list as the current node, deleting the current node from the open list, and adding the current node into the closed list. Wherein the nodes in the closed list represent the nodes to reach.
Fourthly, from the adjacent nodes of the current node, the node which is not added into the opening list and the closing list and can propagate the ultrasonic wave is determined, the node is added into the opening list, the F value of the node is calculated, and the current node is taken as the father node of the node. In the numerical model of the transformer node, the current node comprises six adjacent nodes corresponding to six surfaces corresponding to the cubic shape of the node. The nodes capable of transmitting ultrasonic waves are embodied by the price parameters of the nodes.
If the neighbor node is already in the closed list, it is ignored as indicating that the node has been processed. If it is difficult for an adjacent node to propagate ultrasonic waves, for example, a node corresponding to the internal metal component in the embodiment of the present invention, it may be determined that the price parameter of the node is 1000 (i.e., the node corresponding to the internal metal component), and the node is ignored.
And if the adjacent node is added into the opening list, taking the current node as the father node of the node, and calculating the G value of the node. If the calculated G value is smaller than the existing G value of the node, which indicates that when the current node is taken as the father node of the node, the price distance of the node from the starting point is shorter and the shortest propagation path requirement is better met, the G value and the F value of the node are recalculated, and the father node of the node is changed into the current node.
And repeating the searching process of the adjacent nodes corresponding to the third step and the fourth step until the end point is added into the closing list.
Fifthly, father nodes are searched forward from the end point in sequence, and the path formed by all the searched nodes is used as the shortest propagation path, so that the shortest propagation path of the ultrasonic signal generated when one transformer oil node is used as the starting point to reach one ultrasonic sensor can be determined.
Further, when the price distance G value of the node is determined, the price distance G is passed throughi＝Gi-1+Hi，G0And (5) determining the value as 0. Wherein G isi-1Representing the G value of the parent node of the node.
When the estimated price distance H value of the node is determined, a plurality of corresponding estimated price distances between the current node and the terminal can be determined according to the directions of all nodes adjacent to the current node in the transformer node numerical model. And then, determining the minimum estimated price distance from the plurality of estimated price distances, and taking the minimum estimated price distance which meets the requirement of the shortest propagation path as the H value of the current node.
When the node is a cube, the adjacent nodes are six nodes corresponding to the six faces of the cube respectively.
Further, the estimated price distance H value is passed through Hi＝min(Hi1,Hi2,Hi3,Hi4,Hi5,Hi6) It is determined that,
wherein the content of the first and second substances,
Hi1,Hi2,Hi3,Hi4,Hi5,Hi6representing six estimated price distances, x, corresponding to six directions of the nodei、yi、ziRepresenting coordinates of nodes, m representing the maximum number of nodes included in the x-axis direction in the model, n representing the maximum number of nodes included in the y-axis direction, l representing the maximum number of nodes included in the z-axis direction, and xend、yend、zendCoordinates representing the end point, poilRepresenting a price parameter for the transformer oil node.
Therefore, in the process of calculating the estimated price distance H value, the corresponding estimated price distance propagated in the transformer oil node is calculated according to the price parameter of the transformer oil node based on the path of the ultrasonic signal propagated from the transformer oil to the outside, the corresponding estimated price distance propagated in the metal tank wall node is calculated according to the price parameter of the metal tank wall, and the estimated price distances propagated to the ultrasonic sensor by the ultrasonic signal are obtained by adding. If the ultrasonic signal is hard to propagate through the internal metal member, the ultrasonic signal is not included in the calculation range, and p is based onoilThe same principle of multiplication, since the price parameter of the metal box wall node is 1, the expression of multiplication by 1 is omitted.
By the method, the shortest propagation path of the generated ultrasonic signals reaching each ultrasonic sensor arranged outside the transformer can be determined when the node corresponding to each transformer oil is used as the position of the partial discharge source.
Then, for a plurality of nodes where there may be a local discharge source, the propagation time of the ultrasonic signal from the node to each ultrasonic sensor can be calculated according to the shortest propagation path corresponding to each node. Further, the estimated arrival time difference of the node at each ultrasonic sensor can be calculated from the calculated different propagation times.
Specifically, the calculation of the propagation time corresponding to each shortest propagation path may be implemented by the following formula:
where t denotes the propagation time, x1、y1、z1Coordinates, x, representing the starting point of the shortest propagation pathj、yj、zjCoordinates, v, representing nodes in the shortest propagation pathjThe speed parameters representing the nodes are different according to the transformer structure, and comprise voil、vsteelNum represents the number of all nodes in the shortest propagation path, and dl represents the corresponding side length of the nodes in the transformer node data model.
S104: and determining the position of the local discharge source from a plurality of nodes by a mixed frog-leaping algorithm according to the actual time difference of arrival and the estimated time difference of arrival.
In the embodiment of the invention, the smaller the difference between the measured actual arrival time difference and the calculated estimated arrival time difference is, the more likely the node corresponding to the estimated arrival time difference is to be the position of the partial discharge source.
Thus, the fitness value corresponding to each node can be calculated by the leapfrog algorithm and the difference between the actual time difference of arrival and the estimated time difference of arrival. And iteration and swarm optimization are carried out through a mixed frog-leaping algorithm, and the fitness value corresponding to each node is recalculated, so that the node with the highest fitness value is determined from a plurality of nodes possibly having local discharge sources and serves as the position of the local discharge source. Wherein, the higher the fitness value, the higher the probability that the position of the corresponding node is the position of the partial discharge source.
Specifically, the embodiment of the present invention takes 4 ultrasonic sensors as an example, and describes in detail a process of determining the position of the partial discharge source by the leapfrog algorithm.
First, M frogs are randomly generated, and the ith frog represents the solution of the problem as Pi＝(xi,yi,zi)。
Calculating the fitness value corresponding to each node according to the difference between the actual arrival time difference and the estimated arrival time difference, which can be realized by the following formula:
wherein, PiA frog representing the calculated fitness value;
an actual value representing the difference in arrival time of the ultrasonic signals received by the sensor 1 and the sensor 2,
indicating that the source of partial discharge is located at PiThe estimated value of the wavelength time difference of the ultrasonic signals received by the sensor 1 and the sensor 2 is obtained by calculation;
an actual value representing the difference in arrival time of the ultrasonic signals received by the sensor 1 and the sensor 3,
indicating that the source of partial discharge is located at PiThe estimated value of the wavelength time difference of the ultrasonic signals received by the sensor 1 and the sensor 3 is obtained by calculation;
an actual value representing the difference in arrival time of the ultrasonic signals received by the sensor 1 and the sensor 4,
indicating that the source of partial discharge is located at PiAt the position of (3), the estimated value of the wavelength time difference between the ultrasonic signals received by the sensor 1 and the sensor 4 is calculated.
Similarly, when there are a plurality of ultrasonic sensors, a corresponding fitness value may be calculated by performing a series of operations such as subtraction, square, and addition of the corresponding actual arrival time difference and the estimated arrival time difference based on the same calculation principle as the formula.
And secondly, arranging the M frogs from good to bad according to the fitness value, and dividing the M frogs into N sub-groups. Wherein, the frogs ranked 1 are divided into the 1 st sub-group, the frogs ranked 2 are divided into the 2 nd sub-group, the frogs ranked N are divided into the N th sub-group, the frogs ranked N +1 are divided into the 1 st sub-group, and so on until all frogs are divided.
Then, the local number of iterations c within the subgroup is given1Number of iterations mixed with global c2Maximum leapfrog step length smaxWith the minimum leapfrog step length smin. Repeating the following step c for each sub-population1Secondly:
first, the best individual P in the current subgroup is updated according to fitness valuebWith the global best individual PgDetermining the position P of the worst individual in the current iteration sub-populationw。
Second, P is treated according to the following strategywUpdating:
(1) updating the leapfrog step length according to the following formula:
six＝int(rand(0,1)×(xb-xw))
siy＝int(rand(0,1)×(yb-yw))
siz＝int(rand(0,1)×(zb-zw))
(||smax||≤||si||≤||smin||)
wherein x isb、yb、zbRepresents the best individual P in the current subgroupbCoordinate of (a), xw、yw、zwRepresents the worst individual P in the current subgroupwThe coordinates of (a).
new Pw＝Pw+si
If new PwThe adaptability value of the method is superior to the original PwThen use the new PwBy substitution of the original Pw。
(2) Otherwise, updating the leapfrog step length according to the following formula:
six＝int(rand(0,1)×(xg-xw))
siy＝int(rand(0,1)×(yg-yw))
siz＝int(rand(0,1)×(zg-zw))
(||smax||≤||si||≤||smin||)
wherein x isg、yg、zgRepresenting the globally best individual PgCoordinate of (a), xw、yw、zwRepresents the worst individual P in the current subgroupwThe coordinates of (a).
new Pw＝Pw+si
If new PwThe adaptability value of the method is superior to the original PwThen use the new PwBy substitution of the original Pw。
(3) Otherwise, a new P is randomly generatedw。
Thirdly, after all sub-groups complete the local depth search, if the global mixed iteration times c are satisfied2And when the evolution process is finished, outputting a global optimal value, namely the position of the local discharge source. Otherwise, remixing, sorting and dividing all the frog individuals, and turning to the first step for recalculation.
In the embodiment of the invention, the propagation path of the ultrasonic signal is considered under the blocking effect of the metal components in the transformer on the ultrasonic signal through the A-path searching algorithm, and the estimated arrival time differences corresponding to all positions of the transformer oil are optimized through the mixed frog-leaping algorithm so as to determine the most possible position of the local discharge source and realize the positioning of the local discharge source. Each position of the transformer oil is accurately divided in a numerical model building mode, so that the accuracy of the positioning of the partial discharge source is improved, the interference of metal parts in the transformer is avoided, iteration optimization is combined, the accuracy of the positioning of the partial discharge source is improved, and the method is high in feasibility and easy to realize.
Fig. 2 is a flowchart of another ultrasonic positioning method for a partial discharge source of a transformer according to an embodiment of the present invention, which specifically includes the following steps:
and acquiring the structure and the size of the transformer to be detected, and constructing a corresponding transformer node numerical model.
The method comprises the steps of arranging a plurality of ultrasonic sensors on the outer wall of a transformer, detecting ultrasonic signals generated by a local discharge source, and determining the actual arrival time difference among the sensors.
And traversing each node of the transformer node numerical model, wherein the node may have a partial discharge source, and determining the shortest propagation path of the corresponding ultrasonic signal to each ultrasonic sensor, the corresponding propagation time and the estimated arrival time difference when each node is taken as the position of the partial discharge source respectively based on an A-path searching algorithm.
And performing iterative optimization based on a mixed frog-leaping algorithm according to the actual time difference of arrival and the estimated time difference of arrival to determine a global optimal solution, namely a node with the highest fitness value as the position of the local discharge source.
It should be noted that the principle of the ultrasonic positioning method for the transformer partial discharge source shown in fig. 1 and fig. 2 is substantially the same, and therefore, the parts not described in detail in fig. 2 may specifically refer to the related description of fig. 1, and the details of the present invention are not repeated herein.
Based on the same inventive concept, the above ultrasonic positioning method for the transformer partial discharge source provided in the embodiment of the present invention further provides a corresponding ultrasonic positioning device for the transformer partial discharge source, as shown in fig. 3.
Fig. 3 is a schematic structural diagram of an ultrasonic positioning apparatus for a transformer partial discharge source according to an embodiment of the present invention, which specifically includes:
a building module 301, which builds a transformer node numerical model according to the transformer structure;
a receiving module 302, configured to receive actual arrival time differences measured by a plurality of ultrasonic sensors disposed on the transformer, where the arrival time differences are actual arrival time differences of ultrasonic signals generated by a local discharge source reaching the plurality of ultrasonic sensors;
a path determining module 303, configured to traverse a plurality of nodes in the transformer node numerical model, and determine estimated arrival time differences of the ultrasonic signals respectively arriving at the ultrasonic sensors from the plurality of nodes through an ultrasonic path search algorithm based on an a-x routing algorithm;
and the position determining module 304 is used for determining the position of the partial discharge source from the plurality of nodes through a mixed frog-leaping algorithm according to the actual time difference of arrival and the estimated time difference of arrival.
The embodiments of the present invention are described in a progressive manner, and the same and similar parts among the embodiments can be referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, as for the apparatus embodiment, since it is substantially similar to the method embodiment, the description is relatively simple, and for the relevant points, reference may be made to the partial description of the method embodiment.
It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.
The above description is only an example of the present invention, and is not intended to limit the present invention. Various modifications and alterations to this invention will become apparent to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention should be included in the scope of the claims of the present invention.
Claims (10)
1. An ultrasonic positioning method for a partial discharge source of a transformer is characterized by comprising the following steps:
constructing a transformer node numerical model according to the transformer structure;
receiving actual arrival time differences measured by a plurality of ultrasonic sensors arranged on the transformer, wherein the actual arrival time differences are actual arrival time differences of ultrasonic signals generated by a local discharge source and arriving at the plurality of ultrasonic sensors;
traversing a plurality of nodes in the transformer node numerical model, and determining estimated arrival time differences of the ultrasonic signals from the plurality of nodes to the plurality of ultrasonic sensors respectively through an ultrasonic path search algorithm based on an A-x routing algorithm;
and determining the position of the partial discharge source from the plurality of nodes by a mixed frog-leaping algorithm according to the actual time difference of arrival and the estimated time difference of arrival.
2. The method according to claim 1, wherein constructing a transformer node numerical model based on the transformer structure specifically comprises:
determining nodes corresponding to different structures in the constructed corresponding transformer node numerical model and position parameters, price parameters and speed parameters corresponding to the nodes respectively according to different structures of different transformers; the node numerical model corresponds to the transformer, and each node in the node numerical model corresponds to a corresponding position in the transformer respectively;
the position parameter represents the position of a node, the price parameter represents the difficulty of the ultrasonic signal propagating in the node, the speed parameter represents the propagation speed of the ultrasonic signal in the node, and the price parameter and the speed parameter are related to the structure of the transformer.
3. The method of claim 2, wherein the a-routing algorithm includes an estimated price distance H value for a node, the estimated price distance H value representing an estimated distance of the node to reach an end of the corresponding shortest propagation path, the estimated price distance H value being related to the price parameter;
the estimated price distance H value is determined by:
determining a plurality of corresponding estimated price distances between the nodes and a terminal point according to the directions of a plurality of nodes adjacent to the nodes in the transformer node numerical model;
and determining the minimum estimated price distance from the plurality of estimated price distances as the H value of the node.
4. The method of claim 3, wherein the estimated price distance H is represented by Hi＝min(Hi1,Hi2,Hi3,Hi4,Hi5,Hi6) It is determined that,
5. The method according to claim 2, wherein determining the estimated arrival time differences of the ultrasonic signals respectively arriving at the ultrasonic sensors from the nodes by an ultrasonic path search algorithm based on an a-x-ray-seeking algorithm comprises:
determining the shortest propagation paths of the ultrasonic signals from the nodes to the corresponding ultrasonic sensors respectively through an ultrasonic path search algorithm based on an A-x routing algorithm;
calculating to obtain the propagation time corresponding to each shortest propagation path according to the determined shortest propagation path;
and determining estimated arrival time differences of the ultrasonic signals respectively arriving at the ultrasonic sensors from the nodes according to the calculated propagation time.
6. The method according to claim 5, wherein the step of calculating the propagation time corresponding to each shortest propagation path according to the determined shortest propagation path specifically comprises:
according to
Calculating the propagation time; where t denotes the propagation time, x1、y1、z1Coordinates, x, representing the starting point of the shortest propagation pathj、yj、zjCoordinates, v, representing nodes in the shortest propagation pathjAnd expressing the speed parameters of the nodes, num expresses the number of all nodes in the shortest propagation path, and dl expresses the corresponding side length of the nodes in the transformer node data model.
7. The method according to claim 1, wherein determining the location of the partial discharge source from the plurality of nodes by a mixed frog-jump algorithm based on the actual time difference of arrival and the estimated time difference of arrival comprises:
calculating the fitness value corresponding to each node respectively through a mixed frog-leaping algorithm and the difference between the actual arrival time difference and the estimated arrival time difference;
iteration and swarm optimization are carried out through a mixed frog-leaping algorithm, the fitness value corresponding to each node is recalculated, and the node with the highest fitness value is determined from the plurality of nodes and serves as the position of the local discharge source.
8. The method of claim 7, wherein the ultrasonic sensor is a plurality of;
the calculating the fitness value corresponding to each node respectively through a mixed frog-leaping algorithm and the difference between the actual arrival time difference and the estimated arrival time difference specifically comprises the following steps:
by passing
Calculating the fitness value corresponding to each node;
wherein, PiA frog representing the calculated fitness value;
representing the actual time difference of arrival of the ultrasonic signals received by the first and second sensors,
indicating that the source of partial discharge is located at PiWhen the position of the ultrasonic sensor is detected, the estimated wavelength time difference of the ultrasonic signals received by the first sensor and the second sensor is calculated;
representing the actual time difference of arrival of the ultrasonic signals received by the first and third sensors,
indicating that the source of partial discharge is located at PiWhen the position of the sensor is within the preset range, the estimated wavelength time difference of the ultrasonic signals received by the first sensor and the third sensor is calculated;
representing the actual time difference of arrival of the ultrasonic signals received by the first and fourth sensors,
indicating that the source of partial discharge is located at PiAt the position of (2), the first sensor obtained by calculationAnd an estimated wavelength time difference of the ultrasonic signal received by the fourth sensor.
9. The method of claim 5, wherein the transformer comprises transformer oil, internal metal components, and a metal tank wall;
determining the shortest propagation path of the ultrasonic signal from the plurality of nodes to the corresponding ultrasonic sensor respectively through an ultrasonic path search algorithm based on an a-x routing algorithm, specifically comprising:
and determining the shortest propagation path of the ultrasonic signals from the plurality of nodes to the corresponding ultrasonic sensor respectively by an ultrasonic path search algorithm based on the A-th routing algorithm, wherein the shortest propagation path does not include the node corresponding to the internal metal component.
10. An ultrasonic locating device of a transformer partial discharge source is characterized by comprising:
the building module is used for building a transformer node numerical model according to the transformer structure;
the receiving module is used for receiving actual arrival time differences measured by a plurality of ultrasonic sensors arranged on the transformer, wherein the actual arrival time differences are actual arrival time differences of ultrasonic signals generated by a local discharge source and arriving at the ultrasonic sensors;
the path determining module is used for traversing a plurality of nodes in the transformer node numerical model and determining estimated arrival time differences of the ultrasonic signals from the nodes to the ultrasonic sensors respectively through an ultrasonic path searching algorithm based on an A-x routing algorithm;
and the position determining module is used for determining the position of the partial discharge source from the plurality of nodes through a mixed frog-leaping algorithm according to the actual time difference of arrival and the estimated time difference of arrival.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110006160.7A CN112834849B (en) | 2021-01-05 | 2021-01-05 | Ultrasonic positioning method and device for partial discharge source of transformer |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110006160.7A CN112834849B (en) | 2021-01-05 | 2021-01-05 | Ultrasonic positioning method and device for partial discharge source of transformer |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112834849A true CN112834849A (en) | 2021-05-25 |
| CN112834849B CN112834849B (en) | 2022-07-19 |
Family
ID=75927601
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110006160.7A Active CN112834849B (en) | 2021-01-05 | 2021-01-05 | Ultrasonic positioning method and device for partial discharge source of transformer |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112834849B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001033510A (en) * | 1999-07-23 | 2001-02-09 | Kansai Electric Power Co Inc:The | Method and equipment for searching generation position of distribution line partial discharge |
| US20090012727A1 (en) * | 2005-03-04 | 2009-01-08 | Wah Hoon Siew | Detecting Partial Discharge in High Voltage Cables |
| KR20170052205A (en) * | 2015-11-04 | 2017-05-12 | 한국전기연구원 | Performance testing system and method for motor drive controller for machine tools |
| CN107132459A (en) * | 2017-03-31 | 2017-09-05 | 国网浙江省电力公司电力科学研究院 | A kind of partial discharge of transformer ultrasound locating method |
| CN108896893A (en) * | 2018-09-13 | 2018-11-27 | 国网安徽省电力有限公司电力科学研究院 | The positioning system and localization method of a kind of Partial Discharge Sources in electrical equipment |
| CN109142996A (en) * | 2018-08-24 | 2019-01-04 | 西南交通大学 | A kind of test method of oil-immersed and high-voltage casing moisture point of entry |
| CN110288075A (en) * | 2019-06-26 | 2019-09-27 | 电子科技大学 | A kind of feature selection approach based on improvement shuffled frog leaping algorithm |
-
2021
- 2021-01-05 CN CN202110006160.7A patent/CN112834849B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001033510A (en) * | 1999-07-23 | 2001-02-09 | Kansai Electric Power Co Inc:The | Method and equipment for searching generation position of distribution line partial discharge |
| US20090012727A1 (en) * | 2005-03-04 | 2009-01-08 | Wah Hoon Siew | Detecting Partial Discharge in High Voltage Cables |
| KR20170052205A (en) * | 2015-11-04 | 2017-05-12 | 한국전기연구원 | Performance testing system and method for motor drive controller for machine tools |
| CN107132459A (en) * | 2017-03-31 | 2017-09-05 | 国网浙江省电力公司电力科学研究院 | A kind of partial discharge of transformer ultrasound locating method |
| CN109142996A (en) * | 2018-08-24 | 2019-01-04 | 西南交通大学 | A kind of test method of oil-immersed and high-voltage casing moisture point of entry |
| CN108896893A (en) * | 2018-09-13 | 2018-11-27 | 国网安徽省电力有限公司电力科学研究院 | The positioning system and localization method of a kind of Partial Discharge Sources in electrical equipment |
| CN110288075A (en) * | 2019-06-26 | 2019-09-27 | 电子科技大学 | A kind of feature selection approach based on improvement shuffled frog leaping algorithm |
Non-Patent Citations (1)
| Title |
|---|
| 邵先军 等: "GIS特高频局部放电检测与诊断技术的研究进展", 《西安交通大学学报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112834849B (en) | 2022-07-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kim et al. | Landmark-guided subgoal generation in hierarchical reinforcement learning | |
| Moore et al. | Room geometry estimation from a single channel acoustic impulse response | |
| CN103888205B (en) | A kind of method for predicting electromagnetic wave propagation and device | |
| CN103529364B (en) | Many Partial Discharge Sources localization method in oil based on ultrasonic thinned array sensor | |
| CN110375753A (en) | Map-matching method, device, server and storage medium | |
| CN101960329A (en) | Method and device for trilateration using los link prediction and pre-measurement los path filtering | |
| CN108376116A (en) | Based on the method for generating test case for improving particle cluster algorithm | |
| CN111368881A (en) | Low-frequency GPS track road network matching method based on multidimensional data fusion analysis | |
| CN110823226B (en) | Unmanned aerial vehicle intelligent route planning test method based on metamorphic test technology | |
| Dudani et al. | Partial discharge detection in transformer using adaptive grey wolf optimizer based acoustic emission technique | |
| Wang et al. | Partial discharge localisation methodology for power transformers based on improved acoustic propagation route search algorithm<? show [AQ ID= Q1]?> | |
| Feng et al. | Target tracking based on improved square root cubature particle filter via underwater wireless sensor networks | |
| Kelly et al. | Detecting arrivals in room impulse responses with dynamic time warping | |
| CN114742307A (en) | Wave element prediction method and system | |
| CN104849573A (en) | Simulation prediction method for strong electromagnetic pulse environment of ship platform in marine environment | |
| James et al. | Pekeris waveguide comparisons of methods for predicting acoustic field amplitude uncertainty caused by a spatially uniform environmental uncertainty (L) | |
| CN112834849B (en) | Ultrasonic positioning method and device for partial discharge source of transformer | |
| Chen et al. | Local path searching based map matching algorithm for floating car data | |
| Han et al. | Research on node location algorithm of Zigbee based on optimized neural network | |
| CN108594241B (en) | A kind of AUV Sound stealth method for situation assessment | |
| CN105388460A (en) | Indoor underwater target positioning method based on genetic algorithm | |
| KR101436858B1 (en) | Path generating method for autonomous mobile robot using uncertainty | |
| CN112763868B (en) | Local discharge source positioning method and system based on immune particle swarm | |
| CN112084655A (en) | Ground penetrating radar parameter inversion method based on non-monotonic line search | |
| CN110602635B (en) | Indoor map matching enhanced positioning method, device and storage 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 |