CN109142966A - Fault distance-finding method, device, equipment and medium based on route measured data - Google Patents

Fault distance-finding method, device, equipment and medium based on route measured data Download PDF

Info

Publication number
CN109142966A
CN109142966A CN201810774896.7A CN201810774896A CN109142966A CN 109142966 A CN109142966 A CN 109142966A CN 201810774896 A CN201810774896 A CN 201810774896A CN 109142966 A CN109142966 A CN 109142966A
Authority
CN
China
Prior art keywords
waveform
moment
fault
wave
zero mould
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201810774896.7A
Other languages
Chinese (zh)
Inventor
冯瑞发
张巍
廖永力
张贵峰
田治仁
龚博
黄增浩
张志强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CSG Electric Power Research Institute
Research Institute of Southern Power Grid Co Ltd
Original Assignee
Research Institute of Southern Power Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Research Institute of Southern Power Grid Co Ltd filed Critical Research Institute of Southern Power Grid Co Ltd
Priority to CN201810774896.7A priority Critical patent/CN109142966A/en
Publication of CN109142966A publication Critical patent/CN109142966A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/088Aspects of digital computing

Abstract

The present invention discloses the fault distance-finding method based on route measured data, is the zero mould current component for extract after triumphant human relations boolean transformation to actual measurement fault current traveling wave data in monitoring point including obtaining;And it is filtered;Filtered zero mould current component is carried out two points of recursion SVD to decompose, obtains the first layer details coefficients signal of decomposition;At the time of obtaining the Mintrop wave shape of the maximum absolute value of first waveforms amplitude in details coefficients signal waveform, and obtain corresponding according to Mintrop wave shape, as the first moment;In preset data in window, in details coefficients signal waveform, identical waveform carries out correlation processing with the polarity of Mintrop wave shape, obtains correlation coefficient value;At the time of the polarity of waveform includes corresponding to the maximum value of wave crest and trough acquisition correlation coefficient value, as the second moment;The distance between fault point and monitoring point is calculated according to the first, second moment and zero mould electric current velocity of wave.Solve the problems, such as that existing positioning is not accurate, ranging reliability is high.

Description

Fault distance-finding method, device, equipment and medium based on route measured data
Technical field
The present invention relates to fault localization technology examination field more particularly to a kind of fault localization sides based on route measured data Method, device, equipment and storage medium.
Background technique
Transmission line of electricity has a very wide distribution, and span is big, and locating orographic condition is different, often due to external world's lightning stroke, external force The reasons such as destruction, bird pest, pollution flashover, screen of trees break down, after failure occurs, if route operation maintenance personnel accurately can not quickly be found out Abort situation will directly affect the fast power restoration of faulty line, not only influences the safe operation of power grid, can also give the people Production and living bring strong influence.
Existing circuit fault distance measurement mainly have fault analytical method based on power frequency quality with based on high frequency transient The Travelling Wave Fault Location method of traveling wave.Fault analytical method is based on route distributed parameter model, by establishing false voltage electric current and event Hinder the functional relation of distance to be positioned, this method is easy by system operation mode, fault point transition resistance and line parameter circuit value The influence of accuracy, range accuracy is not high, is generally used in failure wave-recording at present.Travelling wave ranging method is since there is no failure The ranging defect of analytic approach, range measurement principle is simple, and range accuracy is higher, is widely used in standing on site at present in ranging.But From the point of view of practical ranging situation, field failure range unit ranging accuracy rate is only 70% or so, at present Single Ended Fault Location Cause wave head arrival time accurate since attenuation distortion occurs for fault waveform persistent oscillation, the transmission of fault point back wave wave head Fault point back wave arrival time caused by extracting is difficult to demarcate, and often occurs that range accuracy is poor or even the situation of ranging failure.
Summary of the invention
In view of this, the embodiment of the present invention provides a kind of fault distance-finding method based on route measured data, device, equipment And storage medium, to solve present in existing travelling wave ranging technology since fault traveling wave persistent oscillation, waveform exist centainly The reasons such as attenuation distortion cause fault point back wave arrival time to be difficult to the problem of demarcating, and realize and are accurately positioned.
In a first aspect, the embodiment of the invention provides a kind of fault distance-finding methods based on route measured data, including obtain Take zero mould current component of monitoring point;Wherein, the zero mould current component is to carry out triumphant human relations to actual measurement fault current traveling wave data It extracts and obtains after boolean's transformation;
The zero mould current component is filtered;
The filtered zero mould current component is carried out two points of recursion SVD to decompose, obtains the first layer details point of decomposition Measure signal;
The Mintrop wave shape of the maximum absolute value of first waveforms amplitude in the details coefficients signal waveform is obtained, and according to institute At the time of stating the maximum absolute value acquisition correspondence of waveforms amplitude, wherein as the first moment at the time of described corresponding;
According to the polarity of the Mintrop wave shape, in preset data in window, in the details coefficients signal waveform with institute The identical waveform of polarity for stating Mintrop wave shape carries out correlation processing, to obtain corresponding correlation coefficient value;Wherein, the waveform Polarity includes wave crest and trough;
At the time of acquisition corresponding to the maximum value of the correlation coefficient value, as the second moment;
Fault point and the prison is calculated according to first moment, second moment and the zero mould electric current velocity of wave The distance between measuring point.
In the first possible implementation of first aspect, the polarity according to the Mintrop wave shape, in preset number According to when window in, in the details coefficients signal waveform, identical waveform carries out correlation processing with the polarity of the Mintrop wave shape, To obtain corresponding correlation coefficient value;Wherein, the polarity of the waveform includes wave crest and trough includes:
In the details coefficients signal waveform, extraction and the consistent waveform of the first polarity of wave, as wave to be analyzed Shape;
According to window when the preset data, the Mintrop wave shape is subjected to correlation with each waveform to be analyzed respectively Processing obtains the corresponding correlation coefficient value of each waveform to be analyzed.
The possible implementation of with reference to first aspect the first, in second of possible implementation of first aspect, institute It states according to window when the preset data, the Mintrop wave shape is subjected to correlation processing with each waveform to be analyzed respectively, Obtaining the corresponding correlation coefficient value of each waveform to be analyzed includes:
Discretization is carried out to the Mintrop wave shape and each waveform to be analyzed, respectively obtains the corresponding number of the Mintrop wave shape Value signal sequence and the corresponding digital signal sequences of each waveform to be analyzed;
In preset data in window, according to related coefficient calculation formula to the corresponding numerical signal sequence of the Mintrop wave shape Digital signal sequences corresponding with each waveform to be analyzed carry out correlation analysis, obtain each waveform pair to be analyzed The correlation coefficient value answered.
The possible implementation of second with reference to first aspect, in the third possible implementation of first aspect, institute Stating related coefficient calculation formula includes:
Wherein, N is the preset data time window length, x (n) For the corresponding numerical signal sequence of the Mintrop wave shape, y (n) is the corresponding digital signal sequences of the waveform to be analyzed, ρxyFor phase Coefficient values.
In the 4th kind of possible implementation of first aspect, the zero mould current component for obtaining monitoring point;Wherein, institute State zero mould current component be extracted after actual measurement fault current traveling wave data are carried out with triumphant human relations boolean transformation obtain include:
A cycle current traveling wave signal of A, B and C three-phase after the failure of acquisition monitoring point acquisition, wherein the electric current of A phase Travelling wave signal is ia, the current traveling wave signal of B phase is ib, the current traveling wave signal of C phase is ic
Obtain zero mould i of the electric current0
To the ia, the ib, the ic, the i0, the i1With the i2Triumphant human relations boolean transformation is carried out, zero mould point is obtained Measure signal
It is described that the filtered zero mould current component is carried out two in the 5th kind of possible implementation of first aspect Recursion SVD is divided to decompose, the first layer details coefficients signal for obtaining decomposition includes:
Obtain filtered zero mould current component i0=[x1,x2,x3,...,xN], wherein N is sampling number;
Hankel matrix H is constructed to the filtered zero mould current component, wherein
SVD decomposition is carried out to the H, obtains two singular value σaAnd σj, wherein the σaCorresponding approximation component signal Aj, The σjCorresponding details coefficients signal Dj
Obtain the first layer details coefficients signal D decomposedj
In the 6th kind of possible implementation of first aspect, it is described according to first moment, second moment and The zero mould electric current velocity of wave is calculated the distance between fault point and the monitoring point and includes:
Failure is calculated according to first moment, second moment and the zero mould electric current velocity of wave and range formula The distance between point and the monitoring point;Wherein, the range formula isWherein, t1When being described first It carves, t2For first moment, v0For the zero mould electric current velocity of wave, l is the distance between the fault point and the monitoring point.
Second aspect, the embodiment of the invention provides a kind of fault location devices based on route measured data, comprising:
Zero mould current component obtains module, for obtaining zero mould current component of monitoring point;Wherein, the zero mould electric current point Amount obtains extract after triumphant human relations boolean converts to actual measurement fault current traveling wave data;
Filter module, for being filtered to the zero mould current component;
Details coefficients signal acquisition module, for the filtered zero mould current component to be carried out two points of recursion SVD points Solution, obtains the first layer details coefficients signal of decomposition;
First moment obtained module, for obtaining the absolute value of first waveforms amplitude in the details coefficients signal waveform Maximum Mintrop wave shape, and at the time of obtain corresponding according to the maximum absolute value of the waveforms amplitude, wherein at the time of described corresponding As the first moment;
Correlation coefficient value obtains module, for the polarity according to the Mintrop wave shape, in preset data in window, to described Identical waveform carries out correlation processing with the polarity of the Mintrop wave shape in details coefficients signal waveform, to obtain corresponding correlation Coefficient value;Wherein, the polarity of the waveform includes wave crest and trough;
Second moment obtained module, at the time of for obtaining corresponding to the maximum value of the correlation coefficient value, as second Moment;
Distance calculation module, for being calculated according to first moment, second moment and the zero mould electric current velocity of wave Obtain the distance between fault point and the monitoring point.
The third aspect, the embodiment of the invention also provides a kind of fault localization equipment based on route measured data, special Sign is, including processor, memory and stores in the memory and is configured as the meter executed by the processor Calculation machine program, the processor realize that the failure based on route measured data is surveyed as described above when executing the computer program Away from method.
Fourth aspect, the embodiment of the invention also provides a kind of computer readable storage mediums, which is characterized in that the meter Calculation machine readable storage medium storing program for executing includes the computer program of storage, wherein controls the calculating in computer program operation Equipment executes the fault distance-finding method as described above based on route measured data where machine readable storage medium storing program for executing.
Compared with prior art, the fault distance-finding method disclosed by the invention based on route measured data obtains prison first Zero mould current component of measuring point;Wherein, the zero mould current component is to carry out triumphant human relations boolean to actual measurement fault current traveling wave data It extracts and obtains after transformation;The zero mould current component is filtered;The filtered zero mould current component is carried out two points Recursion SVD is decomposed, and obtains the first layer details coefficients signal of decomposition;Obtain first waveform in the details coefficients signal waveform The Mintrop wave shape of the maximum absolute value of amplitude, and at the time of obtain corresponding according to the maximum absolute value of the waveforms amplitude, wherein institute As the first moment at the time of stating corresponding;According to the polarity of the Mintrop wave shape, in preset data in window, to the details point The identical waveform progress correlation processing with the polarity of the Mintrop wave shape in signal waveform is measured, to obtain corresponding related coefficient Value;Wherein, the polarity of the waveform includes wave crest and trough;At the time of acquisition corresponding to the maximum value of the correlation coefficient value, As the second moment;According to first moment, second moment and the zero mould electric current velocity of wave be calculated fault point with The distance between described monitoring point.It is filtered based on field measurement fault current data using the zero _exit of fault current Afterwards, two points of recursion svd algorithms and binding signal correlation analysis method realize the accurate positioning of fault point.This method can be solved effectively Due to can not the subsequent oscillation wave of the initial traveling wave of accurate recognition failure and fault point back wave in certainly traditional Single Terminal Traveling Wave Fault Location method Wave head arrival time and cause positioning accuracy not high, or even the problem of positioning failure, range accuracy are high, and Range finding reliability is high, Meet Practical Project requirement, increases cost without additional when application, it is good in economic efficiency.
Detailed description of the invention
By reading a detailed description of non-restrictive embodiments in the light of the attached drawings below, of the invention other Feature, objects and advantages will become more apparent upon:
Fig. 1 is the schematic diagram for the fault localization equipment based on route measured data that the embodiment of the present invention one provides;
Fig. 2 is the flow diagram of the fault distance-finding method provided by Embodiment 2 of the present invention based on route measured data;
Fig. 3 is the structural schematic diagram for the fault location device based on route measured data that the embodiment of the present invention three provides.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other Embodiment shall fall within the protection scope of the present invention.
Those skilled in the art will appreciate that the present invention can be implemented as equipment, method or computer program product. Therefore, the present disclosure may be embodied in the following forms, it may be assumed that can be complete hardware, be also possible to complete software (including Firmware, resident software, microcode etc.), it can also be the form that hardware and software combines, referred to generally herein as " circuit ", " mould Block " or " system ".In addition, in some embodiments, the present invention is also implemented as in one or more computer-readable mediums In computer program product form, include computer-readable program code in the computer-readable medium.
It can be using any combination of one or more computer-readable media.Computer-readable medium can be calculating Machine readable signal medium or computer readable storage medium.Computer readable storage medium for example can be --- but it is unlimited In system, device or the device of --- electricity, magnetic, optical, electromagnetic, infrared ray or semiconductor, or any above combination.It calculates The more specific example (non exhaustive list) of machine readable storage medium storing program for executing includes: electrical connection with one or more conducting wires, just Taking formula computer disk, hard disk, random access memory (RAM), read-only memory (ROM), erasable type may be programmed read-only storage Device (EPROM or flash memory), optical fiber, portable compact disc read-only memory (CD-ROM), light storage device, magnetic memory device, Or above-mentioned any appropriate combination.In this document, computer readable storage medium can be it is any include or storage journey The tangible medium of sequence, the program can be commanded execution system, device or device use or in connection.
Computer-readable signal media may include in a base band or as carrier wave a part propagate data-signal, Wherein carry computer-readable program code.The data-signal of this propagation can take various forms, including --- but It is not limited to --- electromagnetic signal, optical signal or above-mentioned any appropriate combination.Computer-readable signal media can also be Any computer-readable medium other than computer readable storage medium, which can send, propagate or Transmission is for by the use of instruction execution system, device or device or program in connection.
The program code for including on computer-readable medium can transmit with any suitable medium, including --- but it is unlimited In --- wireless, electric wire, optical cable, RF etc. or above-mentioned any appropriate combination.
The computer for executing operation of the present invention can be write with one or more programming languages or combinations thereof Program code, described program design language include object oriented program language-such as Java, Smalltalk, C++, Further include conventional procedural programming language-such as " C " language or similar programming language.Program code can be with It fully executes, partly execute on the user computer on the user computer, being executed as an independent software package, portion Divide and partially executes or executed on a remote computer or server completely on the remote computer on the user computer.? Be related in the situation of remote computer, remote computer can pass through the network of any kind --- including local area network (LAN) or Wide area network (WAN)-be connected to subscriber computer, or, it may be connected to outer computer (such as mentioned using Internet service It is connected for quotient by internet).
Below with reference to the method for the embodiment of the present invention, the flow chart of equipment (system) and computer program product and/or The block diagram description present invention.It should be appreciated that each box in each box and flowchart and or block diagram of flowchart and or block diagram Combination, can be realized by computer program instructions.These computer program instructions can be supplied to general purpose computer, dedicated The processor of computer or other programmable data processing units, to produce a kind of virtual machine, these computer programs refer to It enables and being executed by computer or other programmable data processing units, produced in the box in implementation flow chart and/or block diagram The device of defined function/operation.
These computer program instructions can also be stored in can make computer or other programmable data processing units In computer-readable medium operate in a specific manner, in this way, the instruction of storage in computer-readable medium just produces one Function/operation command device specified in a box including in implementation flow chart and/or block diagram (instructionmeans) manufacture (manufacture).
Computer program instructions can also be loaded into computer, other programmable data processing units or other equipment On, so that series of operation steps are executed in computer, other programmable data processing units or other equipment, in terms of generating The process that calculation machine is realized, so that the instruction executed on a computer or other programmable device is capable of providing implementation flow chart And/or function/operation process specified in the box in block diagram.
Embodiment one
Referring to Figure 1, Fig. 1 is showing for the fault localization equipment based on route measured data that the embodiment of the present invention one provides It is intended to;In executing the fault distance-finding method provided in an embodiment of the present invention based on route measured data, as shown in Figure 1, this is based on The fault localization equipment of route measured data includes: at least one processor 11, such as CPU, at least one network interface 14 or Person's other users interface 13, memory 15, at least one communication bus 12, communication bus 12 is for realizing between these components Connection communication.Wherein, user interface 13 optionally may include USB interface and other standards interface, wireline interface.Network connects Mouth 14 may include optionally Wi-Fi interface and other wireless interfaces.Memory 15 may include high speed RAM memory, It may further include non-labile memory (non-volatilememory), a for example, at least magnetic disk storage.Memory 15 optionally may include at least one storage device for being located remotely from aforementioned processor 11.
In some embodiments, memory 15 stores following element, executable modules or data structures, or Their subset or their superset:
Operating system 151 includes various system programs, for realizing various basic businesses and hardware based of processing Business;
Program 152.
Specifically, processor 11 executes described in the embodiment of the present invention for calling the program 152 stored in memory 15 Fault distance-finding method based on route measured data.
Alleged processor can be central processing unit (Central Processing Unit, CPU), can also be it His general processor, digital signal processor (Digital Signal Processor, DSP), specific integrated circuit (Application Specific Integrated Circuit, ASIC), ready-made programmable gate array (Field- Programmable Gate Array, FPGA) either other programmable logic device, discrete gate or transistor logic, Discrete hardware components etc..General processor can be microprocessor or the processor is also possible to any conventional processor Deng the processor is the control centre of the fault distance-finding method based on route measured data, utilizes various interfaces and line The various pieces of the entire fault distance-finding method based on route measured data of road connection.
The memory can be used for storing the computer program and/or module, and the processor is by operation or executes Computer program in the memory and/or module are stored, and calls the data being stored in memory, realization is based on The various functions of the electronic device of the fault localization of route measured data.The memory can mainly include storing program area and deposit Store up data field, wherein storing program area can application program needed for storage program area, at least one function (for example sound is broadcast Playing function, text conversion function etc.) etc.;Storage data area, which can be stored, uses created data (such as audio according to mobile phone Data, text message data etc.) etc..In addition, memory may include high-speed random access memory, it can also include non-volatile Property memory, such as hard disk, memory, plug-in type hard disk, intelligent memory card (Smart Media Card, SMC), secure digital (Secure Digital, SD) card, flash card (Flash Card), at least one disk memory, flush memory device or other Volatile solid-state part.
Wherein, if the integrated module of the fault localization based on route measured data is in the form of SFU software functional unit Realize and when sold or used as an independent product, can store in a computer readable storage medium.Based on this The understanding of sample, the present invention realize all or part of the process in above-described embodiment method, can also be referred to by computer program Relevant hardware is enabled to complete, the computer program can be stored in a computer readable storage medium, the computer Program is when being executed by processor, it can be achieved that the step of each embodiment of the method for the invention.Wherein, the computer program includes Computer program code, the computer program code can for source code form, object identification code form, executable file or certain A little intermediate forms etc..The computer-readable medium may include: any entity that can carry the computer program code Or device, recording medium, USB flash disk, mobile hard disk, magnetic disk, CD, computer storage, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), electric carrier signal, telecommunication signal and software Distribution medium etc..
The side of the fault localization according to an embodiment of the present invention based on route measured data described below with reference to accompanying drawings Method.
Embodiment two
Referring to fig. 2, Fig. 2 is a kind of method of fault localization based on route measured data provided by Embodiment 2 of the present invention Flow diagram, comprising:
S11, the zero mould current component for obtaining monitoring point;Wherein, the zero mould current component is to actual measurement fault current row Wave number is extracted after converting according to the triumphant human relations boolean of progress and is obtained;
S12, the zero mould current component is filtered;
S13, the filtered zero mould current component is carried out to two points of recursion SVD decomposition, the first layer for obtaining decomposition is thin Save component signal;
S14, the Mintrop wave shape and root for obtaining the maximum absolute value of first waveforms amplitude in the details coefficients signal waveform At the time of obtaining corresponding according to the maximum absolute value of the waveforms amplitude, wherein as the first moment at the time of described corresponding;
S15, according to the polarity of the Mintrop wave shape, in preset data in window, in the details coefficients signal waveform Identical waveform carries out correlation processing with the polarity of the Mintrop wave shape, to obtain corresponding correlation coefficient value;Wherein, the wave The polarity of shape includes wave crest and trough;
S16, at the time of obtain corresponding to the maximum value of the correlation coefficient value, as the second moment;
S17, fault point and institute are calculated according to first moment, second moment and the zero mould electric current velocity of wave State the distance between monitoring point.
Specifically, step S11, the zero mould current component for obtaining monitoring point;Wherein, the zero mould current component is pair Actual measurement fault current traveling wave data extract after triumphant human relations boolean transformation obtain include:
A cycle current traveling wave signal of A, B and C three-phase after the failure of acquisition monitoring point acquisition, wherein the electric current of A phase Travelling wave signal is ia, the current traveling wave signal of B phase is ib, the current traveling wave signal of C phase is ic
Obtain zero mould i of the electric current0
To the ia, the ib, the ic, the i0, the i1With the i2Triumphant human relations boolean transformation is carried out, zero mould point is obtained Measure signal
It should be noted that carrying out triumphant human relations boolean to it after live traveling wave ranging device acquisition fault current traveling wave data Decoupling transformation respectively obtains zero mould, 1 mould and 2 mold component signals of transient current, due to single-phase when route is run in practice Ground fault is the most common fault type, and when 1 mould of use or the progress data analysis of 2 mold components, 1 mold component can not identify C phase fault out, 2 mold components can not identify B phase fault, therefore the present invention is calculated using zero _exit, for A, B, C three-phase Ground fault, zero _exit can identify well.
Specifically, step S13, it is described the filtered zero mould current component be subjected to two points of recursion SVD decompose, obtain The first layer details coefficients signal of decomposition includes:
Obtain filtered zero mould current component i0=[x1,x2,x3,...,xN], wherein N is sampling number;
Hankel matrix H is constructed to the filtered zero mould current component, wherein
SVD decomposition is carried out to the H, obtains two singular value σaAnd σj, wherein the σaCorresponding approximation component signal Aj, The σjCorresponding details coefficients signal Dj
Obtain the first layer details coefficients signal D decomposedj
Specifically, step S15, the polarity according to the Mintrop wave shape, in preset data in window, to the details Identical waveform carries out correlation processing with the polarity of the Mintrop wave shape in component signal waveform, to obtain corresponding related coefficient Value;Wherein, the polarity of the waveform includes wave crest and trough includes:
In the details coefficients signal waveform, extraction and the consistent waveform of the first polarity of wave, as wave to be analyzed Shape;
According to window when the preset data, the Mintrop wave shape is subjected to correlation with each waveform to be analyzed respectively Processing obtains the corresponding correlation coefficient value of each waveform to be analyzed.
Preferably, described according to window when the preset data, by the Mintrop wave shape respectively with each wave to be analyzed Shape carries out correlation processing, obtains the corresponding correlation coefficient value of each waveform to be analyzed and includes:
Discretization is carried out to the Mintrop wave shape and each waveform to be analyzed, respectively obtains the corresponding number of the Mintrop wave shape Value signal sequence and the corresponding digital signal sequences of each waveform to be analyzed;
In preset data in window, according to related coefficient calculation formula to the corresponding numerical signal sequence of the Mintrop wave shape Digital signal sequences corresponding with each waveform to be analyzed carry out correlation analysis, obtain each waveform pair to be analyzed The correlation coefficient value answered.
Preferably, the related coefficient calculation formula includes:
Wherein, N is the preset data time window length, x (n) For the corresponding numerical signal sequence of the Mintrop wave shape, y (n) is the corresponding digital signal sequences of the waveform to be analyzed, ρxyFor phase Coefficient values.
It should be noted that in embodiments of the present invention, when the Mintrop wave shape is wave crest, then believing the details coefficients Subsequent wave crest extracts analysis in number waveform, and when the Mintrop wave shape is trough, then to the details coefficients signal wave Subsequent trough extracts analysis in shape.
Preferably, described that event is calculated according to first moment, second moment and the zero mould electric current velocity of wave The distance between barrier point and the monitoring point include:
Failure is calculated according to first moment, second moment and the zero mould electric current velocity of wave and range formula The distance between point and the monitoring point;Wherein, the range formula isWherein, t1When being described first It carves, t2For first moment, v0For the zero mould electric current velocity of wave, l is the distance between the fault point and the monitoring point.
The whole flow process of the embodiment of the present invention is specifically described below:
Collection in worksite fault current traveling wave data carry out triumphant human relations to it after acquiring actual measurement fault current traveling wave data Boolean's transformation, extracts transient current zero _exit signal, then using wavelet threshold denoising method to the zero _exit of extraction into Row filtering, it is assumed that collected A, B, C three-phase fault current traveling wave of traveling wave ranging device is respectively ia、ib、ic, i0、i1、i2Respectively For zero mould, 1 mould and 2 mold components of electric current, it is based on triumphant human relations boolean transformation matrix, available i0、i1、i2With ia、ib、icBetween Relationship, it may be assumed thatIt obtainsZero _exit signal is thus obtained.It will be upper The i stated0Component carries out wavelet threshold denoising, utilizes filtered i0Component construction Hankel matrix carries out two points of recursion SVD points Solution, and retain the details coefficients signal after decomposing for subsequent analysis, using wavelet threshold denoising method to zero mould after decoupling Component signal is filtered, and obtains a series of low frequencies at 3 layers, after decomposition to signal decomposition using db3 wavelet basis when filtering Then component signal and high frequency signal carry out hard -threshold processing to high fdrequency component, i.e., all enabling high frequency coefficient is zero, only protect Low frequency coefficient is stayed, then low frequency coefficient is reconstructed, obtains filtered zero _exit signal.Assuming that filtered current zero Mold component i0=[x1,x2,x3,...,xN], N is sampling number, it can then construct following Hankel matrix:SVD decomposition is carried out to H, obtains two singular value σa、σj, respectively correspond the main body general picture letter of signal Cease AjWith detailed information Dj, then proceed to utilize σaCorresponding approximation component signal construction H-matrix carries out SVD decomposition.Assuming that carrying out The approximation component obtained after j-1 decomposition is Aj-1, enable Aj-1=[aj-1,1,aj-1,2,...,aj-1,N], it can then construct following Hankel matrix:To HjSVD decomposition is carried out, is obtained:Wherein Uj= (uj1,uj2), Uj∈R2×2, Vj=(vj1,vj2,...,vj(N-1)), Vj∈R(N-1)×(N-1), Sj=(diag (σajdj), o), Sj∈ R2×(N-1), σaj、σdjSingular value respectively after jth time SVD decomposition, i.e.,U in formulaji∈R2 ×1, vji=R(N-1)×1, i=1,2,For the approximate part of signal,It, then can benefit for the detail section of signal The approximation component A after signal jth time is decomposed is found out with themjWith details coefficients Dj, calculated approximation component AjContinue to construct Hankel matrix simultaneously repeats the above process the decomposition carried out next time, until i0Resolve into a series of details coefficients and approximation component Until, i.e. i0=Aj+D1+D2+...+Dj.The first layer details coefficients waveform D that decomposition is obtained1It is analyzed, extracts described the One layer of details coefficients waveform D1In first waveform values maximum absolute value Mintrop wave shape, and the absolute value according to the waveform values At the time of maximum obtains corresponding, wherein as the first moment at the time of described corresponding, for example, the details coefficients signal waveform D1Mintrop wave shape be a wave crest, then in preset data window, extract with the first consistent waveform of polarity of wave, that is, extract with The consistent wave crest of Mintrop wave peak polarity, the wave crest extracted makees signal correlation processing with the Mintrop wave peak respectively, described thin Save component signal waveform D1Mintrop wave shape be a trough, then in preset data window, extract with it is described head polarity of wave it is consistent Waveform, that is, extract with the consistent trough of Mintrop wave paddy polarity, the trough extracted makees signal phase with the Mintrop wave paddy respectively The processing of closing property, and calculate each autocorrelation coefficient values.Discretization is carried out to the wave crest or trough extracted, wherein discretization can adopt It is carried out with discrete integration formula, for the two discrete digital signal sequences x (n) and y (n) obtained after discrete, wherein N is The preset data time window length, x (n) are the corresponding numerical signal sequence of the Mintrop wave shape, and y (n) is the wave to be analyzed The corresponding digital signal sequences of shape, ρxyIt is the correlation coefficient ρ in the data window of N in length for correlation coefficient valuexyCalculation formula Are as follows:It is ρ that all correlation coefficient value, which are calculated,xy=(ρ12,..., ρk), k is relevant calculation number, enables ρ0=max (ρ12,...,ρk), it is ρ by correlation coefficient value0Wave crest (paddy) it is corresponding when It carves and was used as the second moment, i.e., be the first moment t at the time of the Mintrop wave peak (paddy) is corresponding1, most with Mintrop wave peak (paddy) related coefficient It is the second moment t at the time of big wave crest (paddy) is corresponding2, zero mould velocity of wave of traveling wave is v0, then fault point and path monitoring can be calculated The distance between point l are as follows:
Compared with prior art, the fault distance-finding method disclosed by the invention based on route measured data obtains prison first Zero mould current component of measuring point;Wherein, the zero mould current component is to carry out triumphant human relations boolean to actual measurement fault current traveling wave data It extracts and obtains after transformation;The zero mould current component is filtered;The filtered zero mould current component is carried out two points Recursion SVD is decomposed, and obtains the first layer details coefficients signal of decomposition;Obtain first waveform in the details coefficients signal waveform The Mintrop wave shape of the maximum absolute value of amplitude, and at the time of obtain corresponding according to the maximum absolute value of the waveforms amplitude, wherein institute As the first moment at the time of stating corresponding;According to the polarity of the Mintrop wave shape, in preset data in window, to the details point The identical waveform progress correlation processing with the polarity of the Mintrop wave shape in signal waveform is measured, to obtain corresponding related coefficient Value;Wherein, the polarity of the waveform includes wave crest and trough;At the time of acquisition corresponding to the maximum value of the correlation coefficient value, As the second moment;According to first moment, second moment and the zero mould electric current velocity of wave be calculated fault point with The distance between described monitoring point.It is filtered based on field measurement fault current data using the zero _exit of fault current Afterwards, two points of recursion svd algorithms and binding signal correlation analysis method realize the accurate positioning of fault point.This method can be solved effectively Due to can not the subsequent oscillation wave of the initial traveling wave of accurate recognition failure and fault point back wave in certainly traditional Single Terminal Traveling Wave Fault Location method Wave head arrival time and cause positioning accuracy not high, or even the problem of positioning failure, range accuracy are high, and Range finding reliability is high.
Embodiment three
It is a kind of fault location device based on route measured data that the embodiment of the present invention three provides referring to Fig. 3, Fig. 3 Structural schematic diagram;Include:
Zero mould current component obtains module 31, for obtaining zero mould current component of monitoring point;Wherein, the zero mould electric current Component is to extract to obtain after carrying out actual measurement fault current traveling wave data triumphant human relations boolean transformation;
Filter module 32, for being filtered to the zero mould current component;
Details coefficients signal acquisition module 33, for the filtered zero mould current component to be carried out two points of recursion SVD It decomposes, obtains the first layer details coefficients signal of decomposition;
First moment obtained module 34, for obtaining in the details coefficients signal waveform the absolute of first waveforms amplitude At the time of being worth maximum Mintrop wave shape, and obtain corresponding according to the maximum absolute value of the waveforms amplitude, wherein when described corresponding It carves and was used as the first moment;
Correlation coefficient value obtains module 35, for the polarity according to the Mintrop wave shape, in preset data in window, to institute The identical waveform progress correlation processing with the polarity of the Mintrop wave shape in details coefficients signal waveform is stated, to obtain corresponding phase Coefficient values;Wherein, the polarity of the waveform includes wave crest and trough;
Second moment obtained module 36, at the time of for obtaining corresponding to the maximum value of the correlation coefficient value, as the Two moment;
Distance calculation module 37, for according to first moment, second moment and the zero mould electric current wave-velocity meter Calculation obtains the distance between fault point and the monitoring point.
Preferably, the correlation coefficient value acquisition module 35 includes:
Waveform acquiring unit to be analyzed, in the details coefficients signal waveform, extracting and the first polarity of wave Consistent waveform, as waveform to be analyzed;
Correlation coefficient value acquiring unit, for according to window when the preset data, by the Mintrop wave shape respectively with it is each The waveform to be analyzed carries out correlation processing, obtains the corresponding correlation coefficient value of each waveform to be analyzed.
Preferably, the correlation coefficient value acquiring unit includes:
Discretization is carried out to the Mintrop wave shape and each waveform to be analyzed, respectively obtains the corresponding number of the Mintrop wave shape Value signal sequence and the corresponding digital signal sequences of each waveform to be analyzed;
In preset data in window, according to related coefficient calculation formula to the corresponding numerical signal sequence of the Mintrop wave shape Digital signal sequences corresponding with each waveform to be analyzed carry out correlation analysis, obtain each waveform pair to be analyzed The correlation coefficient value answered.
Preferably, the related coefficient calculation formula includes:
Wherein, N is the preset data time window length, x (n) For the corresponding numerical signal sequence of the Mintrop wave shape, y (n) is the corresponding digital signal sequences of the waveform to be analyzed, ρxyFor phase Coefficient values.
Preferably, the zero mould current component acquisition module 31 includes:
A cycle current traveling wave signal of A, B and C three-phase after the failure of acquisition monitoring point acquisition, wherein the electric current of A phase Travelling wave signal is ia, the current traveling wave signal of B phase is ib, the current traveling wave signal of C phase is ic
Obtain zero mould i of the electric current0
To the ia, the ib, the ic, the i0, the i1With the i2Triumphant human relations boolean transformation is carried out, zero mould point is obtained Measure signal
Preferably, the details coefficients signal acquisition module 33 includes:
Obtain filtered zero mould current component i0=[x1,x2,x3,...,xN], wherein N is sampling number;
Hankel matrix H is constructed to the filtered zero mould current component, wherein
SVD decomposition is carried out to the H, obtains two singular value σaAnd σj, wherein the σaCorresponding approximation component signal Aj, The σjCorresponding details coefficients signal Dj
Obtain the first layer details coefficients signal D decomposedj
Preferably, the distance calculation module 37 includes:
It is described that fault point and institute are calculated according to first moment, second moment and the zero mould electric current velocity of wave Stating the distance between monitoring point includes:
Failure is calculated according to first moment, second moment and the zero mould electric current velocity of wave and range formula The distance between point and the monitoring point;Wherein, the range formula isWherein, t1When being described first It carves, t2For first moment, v0For the zero mould electric current velocity of wave, l is the distance between the fault point and the monitoring point.
Compared with prior art, the embodiment of the present invention has the following beneficial effects: the zero mould electric current for obtaining monitoring point first Component;Wherein, the zero mould current component is to extract to obtain after carrying out actual measurement fault current traveling wave data triumphant human relations boolean transformation; The zero mould current component is filtered;The filtered zero mould current component is carried out two points of recursion SVD to decompose, is obtained The first layer details coefficients signal of decomposition;Obtain the maximum absolute value of first waveforms amplitude in the details coefficients signal waveform Mintrop wave shape, and at the time of obtain corresponding according to the maximum absolute value of the waveforms amplitude, wherein conduct at the time of described corresponding First moment;According to the polarity of the Mintrop wave shape, in preset data in window, in the details coefficients signal waveform with institute The identical waveform of polarity for stating Mintrop wave shape carries out correlation processing, to obtain corresponding correlation coefficient value;Wherein, the waveform Polarity includes wave crest and trough;At the time of acquisition corresponding to the maximum value of the correlation coefficient value, as the second moment;According to institute State the first moment, second moment and the zero mould electric current velocity of wave be calculated between fault point and the monitoring point away from From.Based on field measurement fault current data, using the zero _exit of fault current, it is filtered rear, two points of recursion svd algorithms And binding signal correlation analysis method realizes the accurate positioning of fault point.This method can effectively solve traditional Single Terminal Traveling Wave Fault Location In method due to can not the subsequent oscillation wave of the initial traveling wave of accurate recognition failure and fault point back wave wave head arrival time and make It is not high at positioning accuracy, or even the problem of positioning failure, range accuracy is high, and Range finding reliability is high.
It should be noted that the apparatus embodiments described above are merely exemplary, wherein described be used as separation unit The unit of explanation may or may not be physically separated, and component shown as a unit can be or can also be with It is not physical unit, it can it is in one place, or may be distributed over multiple network units.It can be according to actual It needs that some or all of the modules therein is selected to achieve the purpose of the solution of this embodiment.In addition, device provided by the invention In embodiment attached drawing, the connection relationship between module indicate between them have communication connection, specifically can be implemented as one or A plurality of communication bus or signal wire.Those of ordinary skill in the art are without creative efforts, it can understand And implement.
The above is a preferred embodiment of the present invention, it is noted that for those skilled in the art For, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also considered as Protection scope of the present invention.
It should be noted that in the above-described embodiments, all emphasizing particularly on different fields to the description of each embodiment, in some embodiment In the part that is not described in, reference can be made to the related descriptions of other embodiments.Secondly, those skilled in the art should also know It knows, the embodiments described in the specification are all preferred embodiments, and related movement and simulation must be that the present invention must Must.

Claims (10)

1. a kind of fault distance-finding method based on route measured data characterized by comprising
Obtain zero mould current component of monitoring point;Wherein, the zero mould current component be to actual measurement fault current traveling wave data into It extracts and obtains after the triumphant human relations boolean transformation of row;
The zero mould current component is filtered;
The filtered zero mould current component is carried out two points of recursion SVD to decompose, obtains the first layer details coefficients letter of decomposition Number;
The Mintrop wave shape of the maximum absolute value of first waveforms amplitude in the details coefficients signal waveform is obtained, and according to the wave At the time of the maximum absolute value of shape amplitude obtains corresponding, wherein as the first moment at the time of described corresponding;
According to the polarity of the Mintrop wave shape, in preset data in window, in the details coefficients signal waveform with the head The identical waveform of the polarity of waveform carries out correlation processing, to obtain corresponding correlation coefficient value;Wherein, the polarity of the waveform Including wave crest and trough;
At the time of acquisition corresponding to the maximum value of the correlation coefficient value, as the second moment;
Fault point and the monitoring point is calculated according to first moment, second moment and the zero mould electric current velocity of wave The distance between.
2. the fault distance-finding method according to claim 1 based on route measured data, which is characterized in that described according to institute The polarity for stating Mintrop wave shape, in preset data in window, to the polarity in the details coefficients signal waveform with the Mintrop wave shape Identical waveform carries out correlation processing, to obtain corresponding correlation coefficient value;Wherein, the polarity of the waveform include wave crest and Trough includes:
In the details coefficients signal waveform, extraction and the consistent waveform of the first polarity of wave, as waveform to be analyzed;
According to window when the preset data, the Mintrop wave shape is carried out at correlation with each waveform to be analyzed respectively Reason obtains the corresponding correlation coefficient value of each waveform to be analyzed.
3. the fault distance-finding method according to claim 2 based on route measured data, which is characterized in that described according to institute The Mintrop wave shape is carried out correlation processing with each waveform to be analyzed respectively, obtained each by window when stating preset data The corresponding correlation coefficient value of the waveform to be analyzed includes:
Discretization is carried out to the Mintrop wave shape and each waveform to be analyzed, respectively obtains the corresponding numerical value letter of the Mintrop wave shape Number sequence and the corresponding digital signal sequences of each waveform to be analyzed;
In preset data in window, according to related coefficient calculation formula to the corresponding numerical signal sequence of the Mintrop wave shape and respectively The corresponding digital signal sequences of a waveform to be analyzed carry out correlation analysis, and it is corresponding to obtain each waveform to be analyzed Correlation coefficient value.
4. the fault distance-finding method according to claim 3 based on route measured data, which is characterized in that the phase relation Counting calculation formula includes:
Wherein, N is the preset data time window length, and x (n) is described The corresponding numerical signal sequence of Mintrop wave shape, y (n) are the corresponding digital signal sequences of the waveform to be analyzed, ρxyFor related coefficient Value.
5. the fault distance-finding method according to claim 1 based on route measured data, which is characterized in that the acquisition prison Zero mould current component of measuring point;Wherein, the zero mould current component is to carry out triumphant human relations boolean to actual measurement fault current traveling wave data After transformation extract obtain include:
A cycle current traveling wave signal of A, B and C three-phase after the failure of acquisition monitoring point acquisition, wherein the current traveling wave of A phase Signal is ia, the current traveling wave signal of B phase is ib, the current traveling wave signal of C phase is ic
Obtain zero mould i of the electric current0
To the ia, the ib, the ic, the i0, the i1With the i2Triumphant human relations boolean transformation is carried out, zero _exit letter is obtained Number
6. the fault distance-finding method according to claim 1 based on route measured data, which is characterized in that described to filter The zero mould current component afterwards carries out two points of recursion SVD and decomposes, and the first layer details coefficients signal for obtaining decomposition includes:
Obtain filtered zero mould current component i0=[x1,x2,x3,...,xN], wherein N is sampling number;
Hankel matrix H is constructed to the filtered zero mould current component, wherein
SVD decomposition is carried out to the H, obtains two singular value σaAnd σj, wherein the σaCorresponding approximation component signal Aj, described σjCorresponding details coefficients signal Dj
Obtain the first layer details coefficients signal D decomposedj
7. the fault distance-finding method according to claim 1 based on route measured data, which is characterized in that described according to institute It states the first moment, second moment and the zero mould electric current velocity of wave and the distance between fault point and the monitoring point is calculated Include:
According to first moment, second moment and the zero mould electric current velocity of wave and range formula be calculated fault point with The distance between described monitoring point;Wherein, the range formula isWherein, t1For first moment, t2 For first moment, v0For the zero mould electric current velocity of wave, l is the distance between the fault point and the monitoring point.
8. a kind of fault location device based on route measured data characterized by comprising
Zero mould current component obtains module, for obtaining zero mould current component of monitoring point;Wherein, the zero mould current component is Actual measurement fault current traveling wave data extract after triumphant human relations boolean converts and are obtained;
Filter module, for being filtered to the zero mould current component;
Details coefficients signal acquisition module is decomposed for the filtered zero mould current component to be carried out two points of recursion SVD, is obtained Take the first layer details coefficients signal of decomposition;
First moment obtained module, for obtaining the maximum absolute value of first waveforms amplitude in the details coefficients signal waveform Mintrop wave shape, and at the time of obtain corresponding according to the maximum absolute value of the waveforms amplitude, wherein conduct at the time of described corresponding First moment;
Correlation coefficient value obtains module, for the polarity according to the Mintrop wave shape, in preset data in window, to the details Identical waveform carries out correlation processing with the polarity of the Mintrop wave shape in component signal waveform, to obtain corresponding related coefficient Value;Wherein, the polarity of the waveform includes wave crest and trough;
Second moment obtained module, at the time of for obtaining corresponding to the maximum value of the correlation coefficient value, as the second moment;
Distance calculation module, for being calculated according to first moment, second moment and the zero mould electric current velocity of wave The distance between fault point and the monitoring point.
9. a kind of fault localization equipment based on route measured data, which is characterized in that including processor, memory and storage In the memory and it is configured as the computer program executed by the processor, the processor executes the computer The fault distance-finding method based on route measured data as described in claim 1 to 8 any one is realized when program.
10. a kind of computer readable storage medium, which is characterized in that the computer readable storage medium includes the calculating of storage Machine program, wherein equipment where controlling the computer readable storage medium in computer program operation is executed as weighed Benefit require any one of 1 to 8 described in the fault distance-finding method based on route measured data.
CN201810774896.7A 2018-07-13 2018-07-13 Fault distance-finding method, device, equipment and medium based on route measured data Pending CN109142966A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810774896.7A CN109142966A (en) 2018-07-13 2018-07-13 Fault distance-finding method, device, equipment and medium based on route measured data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810774896.7A CN109142966A (en) 2018-07-13 2018-07-13 Fault distance-finding method, device, equipment and medium based on route measured data

Publications (1)

Publication Number Publication Date
CN109142966A true CN109142966A (en) 2019-01-04

Family

ID=64800444

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810774896.7A Pending CN109142966A (en) 2018-07-13 2018-07-13 Fault distance-finding method, device, equipment and medium based on route measured data

Country Status (1)

Country Link
CN (1) CN109142966A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109901021A (en) * 2019-04-08 2019-06-18 珠海妙微科技有限公司 Based on the multiple dimensioned SVD converter technique of two-dimentional Hankel matrix
CN113010141A (en) * 2021-03-04 2021-06-22 湖北三赫智能科技有限公司 Method, system and device for extracting signal head wave information and readable storage medium
CN113131453A (en) * 2021-05-25 2021-07-16 华北电力大学 Single-ended traveling wave protection method for flexible direct current transmission line
WO2024004114A1 (en) * 2022-06-30 2024-01-04 三菱電機株式会社 Conductor abnormality sensing device and conductor abnormality sensing method

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105353268A (en) * 2015-10-10 2016-02-24 电子科技大学 Method for judging and positioning distributed traveling wave fault of power transmission line
CN108020754A (en) * 2017-11-24 2018-05-11 山东理工大学 Single-ended traveling wave fault location method based on Waveform Reconstructing

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105353268A (en) * 2015-10-10 2016-02-24 电子科技大学 Method for judging and positioning distributed traveling wave fault of power transmission line
CN108020754A (en) * 2017-11-24 2018-05-11 山东理工大学 Single-ended traveling wave fault location method based on Waveform Reconstructing

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
姚李孝: "基于小波相关性的配电网单相接地故障测距", 《电力自动化设备》 *
赵峰 等: "基于二分递推 SVD 的牵引网故障测距研究", 《控制工程》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109901021A (en) * 2019-04-08 2019-06-18 珠海妙微科技有限公司 Based on the multiple dimensioned SVD converter technique of two-dimentional Hankel matrix
CN109901021B (en) * 2019-04-08 2021-07-16 珠海妙微科技有限公司 Multi-scale SVD (singular value decomposition) transformation method based on two-dimensional Hankel matrix
CN113010141A (en) * 2021-03-04 2021-06-22 湖北三赫智能科技有限公司 Method, system and device for extracting signal head wave information and readable storage medium
CN113131453A (en) * 2021-05-25 2021-07-16 华北电力大学 Single-ended traveling wave protection method for flexible direct current transmission line
CN113131453B (en) * 2021-05-25 2022-08-09 华北电力大学 Single-ended traveling wave protection method for flexible direct current transmission line
WO2024004114A1 (en) * 2022-06-30 2024-01-04 三菱電機株式会社 Conductor abnormality sensing device and conductor abnormality sensing method

Similar Documents

Publication Publication Date Title
CN109142966A (en) Fault distance-finding method, device, equipment and medium based on route measured data
CN101907437B (en) Wavelet difference algorithm-based cable fault localization method
Taniguchi et al. Observation results of lightning shielding for large-scale transmission lines
CN106405339A (en) Power transmission line fault reason identification method based on high and low frequency wavelet feature association
CN106646034B (en) A kind of arrester structural health on-line monitoring analysis system
CN108508320A (en) Arc grounding fault identification method based on harmonic energy and wave distortion feature
CN109712179B (en) Method and device for determining jumper wire electrical clearance of power transmission line
CN105184674B (en) A kind of electric energy metering error prediction technique and equipment
CN110543921A (en) cable early fault identification method based on waveform learning
CN101778333B (en) Detection method and device of microphone state
CN108632832B (en) Network coverage analysis method and system
CN107884739A (en) A kind of powered analysis system of electronic mutual inductor
CN108347698A (en) A kind of on-line off-line event trace analysis method, apparatus and system
CN109581126A (en) A kind of electric power electrification detection system and method
CN104931793A (en) Substation grounding grid grounding resistance acquisition method
CN109116196A (en) A kind of power cable fault discharging sound intelligent identification Method
CN108845177A (en) A kind of lightning stroke monitoring for power distribution network and recognition methods, device and equipment
CN104407278A (en) Extraction method and device of ultrasonic wave characteristic parameters
Stracqualursi et al. Tower models for power systems transients: A review
Chai et al. Tampering detection of digital recordings using electric network frequency and phase angle
CN112595876A (en) Rocket lightning triggering thunder and lightning data acquisition device and method and lightning triggering line length determination method
Diendorfer et al. Ground flash density and lightning exposure of power transmission lines
CN103618287A (en) Lightning protection method and system
CN106599865A (en) Disconnecting link state recognition device and method
CN108646145A (en) A kind of transmission line of electricity flashover tower localization method and system

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20190104