CN111781439A - Power cable partial discharge signal detection method and device - Google Patents

Power cable partial discharge signal detection method and device Download PDF

Info

Publication number
CN111781439A
CN111781439A CN202010466552.7A CN202010466552A CN111781439A CN 111781439 A CN111781439 A CN 111781439A CN 202010466552 A CN202010466552 A CN 202010466552A CN 111781439 A CN111781439 A CN 111781439A
Authority
CN
China
Prior art keywords
partial discharge
power cable
discharge signal
transformation
value
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
Application number
CN202010466552.7A
Other languages
Chinese (zh)
Other versions
CN111781439B (en
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.)
Wuzhou Power Supply Bureau of Guangxi Power Grid Co Ltd
Original Assignee
Wuzhou Power Supply Bureau of Guangxi 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 Wuzhou Power Supply Bureau of Guangxi Power Grid Co Ltd filed Critical Wuzhou Power Supply Bureau of Guangxi Power Grid Co Ltd
Priority to CN202010466552.7A priority Critical patent/CN111781439B/en
Publication of CN111781439A publication Critical patent/CN111781439A/en
Application granted granted Critical
Publication of CN111781439B publication Critical patent/CN111781439B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R19/00Arrangements for measuring currents or voltages or for indicating presence or sign thereof
    • G01R19/0046Arrangements for measuring currents or voltages or for indicating presence or sign thereof characterised by a specific application or detail not covered by any other subgroup of G01R19/00
    • G01R19/0053Noise discrimination; Analog sampling; Measuring transients
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16Spectrum analysis; Fourier analysis
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
    • Y04S10/52Outage or fault management, e.g. fault detection or location

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Testing Relating To Insulation (AREA)

Abstract

The invention discloses a method and a device for detecting a partial discharge signal of a power cable, wherein the method comprises the following steps: performing wavelet transformation on the collected partial discharge signals of the power cable based on the selected wavelet basis and the decomposition layer number to obtain corresponding wavelet coefficients; obtaining a threshold value through calculation based on an algorithm selected by the threshold value, and performing threshold value processing on the wavelet coefficient based on the threshold value to obtain an estimated wavelet coefficient; reconstructing the estimated wavelet coefficient to obtain an estimated value of the partial discharge signal of the power cable; carrying out H-order cross overlapping differential transformation on the estimated value of the power cable partial discharge signal, and calculating to obtain the transformation quantity of the power cable partial discharge signal; and calculating to obtain the average energy value and the maximum value of the transformation quantity of the partial discharge signal of the power cable based on the selected data window. In the implementation of the invention, the identification accuracy of the local discharge signal in a complex noise environment is improved, and the detection efficiency is improved.

Description

Power cable partial discharge signal detection method and device
Technical Field
The invention relates to the technical field of partial discharge, in particular to a method and a device for detecting a partial discharge signal of a power cable.
Background
Power cables are increasingly widely used due to the advantages of safe and reliable power supply, urban beautification and the like. The power cable has high safety operation reliability under the condition of not being influenced by external environment factors and artificial factors, but accessories such as an intermediate joint, a terminal and the like of the cable are easy to have problems in the cable laying and accessory installation and manufacturing processes. After the cable breaks down, the fault must be detected in time, otherwise serious economic loss and social influence can be caused. Therefore, the quick and effective detection of the insulation condition of the cable accessory has important significance for guaranteeing the safe and reliable operation of the cable.
Because partial discharge can well reflect the insulation state of the cable, the online monitoring method based on the partial discharge signal is widely applied and is the most effective method for judging the insulation state of the cable at present. However, when the partial discharge is monitored on line, the electrical equipment is in a normal working state, strong electromagnetic interference exists in the environment, and a partial discharge signal generated by a cable insulation defect is weak and is often submerged in the strong interference. Due to the existence of complex noise interferences such as white noise, periodic narrow-band signals and the like in the field environment, reliable detection and extraction of partial discharge signals are difficult to realize only based on a single time-frequency analysis method. Therefore, how to accurately identify and extract the partial discharge signal of the cable in various noise interferences needs to be solved.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides a method and a device for detecting a partial discharge signal of a power cable, so that the identification accuracy of the partial discharge signal in a complex noise environment is improved.
In order to solve the above technical problem, an embodiment of the present invention provides a method for detecting a partial discharge signal of a power cable, where the method includes:
performing wavelet transformation on the collected partial discharge signals of the power cable based on the selected wavelet basis and the decomposition layer number to obtain corresponding wavelet coefficients;
obtaining a threshold value through calculation based on an algorithm selected by the threshold value, and performing threshold value processing on the wavelet coefficient based on the threshold value to obtain an estimated wavelet coefficient;
reconstructing the estimated wavelet coefficient to obtain an estimated value of the partial discharge signal of the power cable;
carrying out H-order cross overlapping differential transformation on the estimated value of the power cable partial discharge signal, and calculating to obtain the transformation quantity of the power cable partial discharge signal;
and calculating to obtain the average energy value and the maximum value of the transformation quantity of the partial discharge signal of the power cable based on the selected data window.
Optionally, the threshold value is obtained by calculating the algorithm selected based on the threshold value, and the specific calculation formula is as follows:
Figure BDA0002512825300000021
wherein, TvTo representThreshold value, σnRepresenting the standard deviation of the noise and N the length of the signal.
Optionally, the performing H-order cross-lap differential transformation on the estimated value of the power cable partial discharge signal, and obtaining a transformation quantity of the power cable partial discharge signal through calculation includes:
obtaining a transformation coefficient array through calculation based on the selected order H of the cross overlapping differential transformation;
and performing H-order cross overlapping differential transformation on the estimated value of the partial discharge signal of the power cable based on the transformation coefficient array, and calculating to obtain the transformation quantity of the partial discharge signal of the power cable.
Optionally, the obtaining a transform coefficient array by calculation based on the selected order H of the cross-lapped differential transform includes:
the transformation coefficient array of the H order is translated to the left by one bit to obtain a new array;
and calculating to obtain the transform coefficient corresponding to the H order based on the new array.
Optionally, based on the transform coefficient array, performing H-order cross-lapped differential transform on the estimated value of the power cable partial discharge signal, and obtaining a transform quantity of the power cable partial discharge signal through calculation, where a specific calculation formula is as follows:
Figure BDA0002512825300000022
wherein S isH.f(n) represents the amount of conversion of the partial discharge signal of the power cable, and (c)j)HAn array of transform coefficients is represented,
Figure BDA0002512825300000031
represents an estimate of the power cable partial discharge signal, H represents the order, and j represents the number of decomposition layers.
Optionally, the obtaining, by calculation, an average energy value and a maximum value of a transformation quantity of the partial discharge signal of the power cable based on the selected data window further includes: and if the maximum value is larger than the fixed value, detecting that the local discharge signal of the power cable is detected.
Optionally, the setting value is set by using an adaptive threshold, including:
if the discharge signal is not detected, the fixed value E6.setting=Krel.1max(E6.normal) In which E6.normalIs the average energy value of the cable under normal conditions, Krel.1A reliability coefficient greater than 1;
if a discharge signal is detected, the value is fixed
Figure BDA0002512825300000032
Wherein EH(nmax) Is a maximum value, Krel.2A reliability factor greater than 1.
In addition, an embodiment of the present invention further provides a power cable partial discharge signal detection apparatus, where the apparatus includes:
a wavelet transformation module: the system is used for performing wavelet transformation on the collected partial discharge signals of the power cable based on the selected wavelet basis and the decomposition layer number to obtain corresponding wavelet coefficients;
a threshold processing module: the method comprises the steps that a threshold value is obtained through calculation based on an algorithm selected by the threshold value, and threshold value processing is carried out on the wavelet coefficient based on the threshold value to obtain an estimated wavelet coefficient;
a reconstruction module: the wavelet coefficient is used for reconstructing the estimated wavelet coefficient to obtain an estimated value of the partial discharge signal of the power cable;
a cross-lapped differential transform module: the device is used for carrying out H-order cross overlapping differential transformation on the estimated value of the power cable partial discharge signal and obtaining the transformation quantity of the power cable partial discharge signal through calculation;
a calculation module: and the average energy value and the maximum value of the transformation quantity of the partial discharge signal of the power cable are obtained through calculation based on the selected data window.
In the implementation of the invention, a threshold denoising method of wavelet transformation is utilized to suppress white noise interference, and on the basis, the transient characteristics of partial discharge signals are extracted by utilizing wavelet transformation and cross-lapped differential transformation, so that the identification accuracy of the partial discharge signals in a complex noise environment is improved; the cross-lapped differential transformation is suitable for extracting high-frequency transient signals, and detecting partial discharge signals by combining wavelet transformation and cross-lapped differential transformation, can effectively inhibit complex noise interference, and has higher detection sensitivity; the cross-over differential transformation algorithm is simple, and the efficiency of the detection method is greatly improved; in addition, the criterion adopts a self-adaptive threshold, so that the detection sensitivity is ensured, and the detection reliability is improved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
Fig. 1 is a schematic flow chart of a power cable partial discharge signal detection method in an embodiment of the present invention;
fig. 2 is a schematic structural component diagram of a power cable partial discharge signal detection device in an embodiment of the invention;
fig. 3 is a schematic flow chart of a power cable partial discharge signal detection method in an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. 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.
Example one
Referring to fig. 1, fig. 1 is a schematic flow chart of a power cable partial discharge signal detection method according to an embodiment of the present invention.
As shown in fig. 1, a method for detecting a partial discharge signal of a power cable includes:
s11: performing wavelet transformation on the collected partial discharge signals of the power cable based on the selected wavelet basis and the decomposition layer number to obtain corresponding wavelet coefficients;
in the specific implementation process of the invention, wavelet basis and decomposition layer number J are selected, and wavelet transformation is firstly carried out on the acquired signals f (n) to obtain corresponding wavelet coefficients Wj,k
S12: obtaining a threshold value through calculation based on an algorithm selected by the threshold value, and performing threshold value processing on the wavelet coefficient based on the threshold value to obtain an estimated wavelet coefficient;
in the specific implementation process of the invention, a general threshold value selection algorithm is adopted to calculate the threshold value TvAnd for the wavelet coefficient W obtained by decompositionj,kSoft threshold value processing is carried out to obtain estimated wavelet coefficient
Figure BDA0002512825300000051
Specifically, the algorithm selected based on the threshold is calculated to obtain the threshold, and the specific calculation formula is as follows:
Figure BDA0002512825300000052
wherein, TvRepresents the threshold value, σnRepresenting the standard deviation of the noise and N the length of the signal, i.e. the number of sample points.
S13: reconstructing the estimated wavelet coefficient to obtain an estimated value of the partial discharge signal of the power cable;
in the implementation of the invention, the estimated wavelet coefficients are processed
Figure BDA0002512825300000053
Reconstructing to obtain the estimated value of the partial discharge signal of the power cable
Figure BDA0002512825300000054
S14: carrying out H-order cross overlapping differential transformation on the estimated value of the power cable partial discharge signal, and calculating to obtain the transformation quantity of the power cable partial discharge signal;
in a specific implementation process of the present invention, the performing H-order cross-lap differential transformation on the estimated value of the power cable partial discharge signal, and obtaining a transformation quantity of the power cable partial discharge signal through calculation includes: obtaining a transformation coefficient array through calculation based on the selected order H of the cross overlapping differential transformation; and performing H-order cross overlapping differential transformation on the estimated value of the partial discharge signal of the power cable based on the transformation coefficient array, and calculating to obtain the transformation quantity of the partial discharge signal of the power cable.
Wherein, the step of obtaining a transform coefficient array by calculation based on the selected order H of the cross-lapped differential transform comprises: the transformation coefficient array of the H order is translated to the left by one bit to obtain a new array; and calculating to obtain the transform coefficient corresponding to the H order based on the new array. Specifically, order (c)j)1=[1 -1 0 … 0]1×(H+1)Represents an array (c)j)1The new array resulting after a bit shift to the left is then (c'j)1=[0 1 -1 0 … 0]1×(H+1)J is 1, …, and H +1 represents the serial number of the element in the array; then the transform coefficients corresponding to each order can be calculated as: (c)j)h=(cj)h-1-(c′j)h-1,h=2,…,H。
In addition, in the step of performing H-order cross-lapped differential transform on the estimated value of the power cable partial discharge signal based on the transform coefficient array, and obtaining the transform quantity of the power cable partial discharge signal through calculation, a specific calculation formula is as follows:
Figure BDA0002512825300000061
wherein S isH.f(n) represents the amount of conversion of the partial discharge signal of the power cable, and (c)j)HAn array of transform coefficients is represented,
Figure BDA0002512825300000062
represents an estimate of the power cable partial discharge signal, H represents the order, and j represents the number of decomposition layers.
S15: and calculating to obtain the average energy value and the maximum value of the transformation quantity of the partial discharge signal of the power cable based on the selected data window.
In a specific implementation process of the present invention, the calculating, based on the selected data window, to obtain an average energy value and a maximum value of a transformation quantity of the partial discharge signal of the power cable further includes: and if the maximum value is larger than the fixed value, detecting that the local discharge signal of the power cable is detected.
Specifically, a certain time window T is taken to calculate SH.fThe average energy value of (n) is:
Figure BDA0002512825300000063
p is the number of sampling points corresponding to the time window T; and then calculate EHMaximum value E of (n)H(nmax) Wherein n ismaxThe sampling point corresponding to the maximum value is satisfied with EH(nmax)>EH.settingThen detecting as a partial discharge signal, recording nmaxAnd a corresponding maximum value EH(nmax)。
It should be noted that, the fixed value is set by using an adaptive threshold, which includes:
if the discharge signal is not detected, i.e. q is 0, the fixed value E is set6.setting=Krel.1max(E6.normal) In which E6.normalIs the average energy value of the cable under normal conditions, Krel.1To a coefficient of reliability greater than 1, i.e. Krel.1=1.5-2.5;
If a discharge signal is detected, i.e. q > 0, the value is fixed
Figure BDA0002512825300000064
Wherein EH(nmax) Is a maximum value, Krel.2Reliability system of more than 1Number, i.e. Krel.2=1.5-2.5。
In the implementation of the invention, a threshold denoising method of wavelet transformation is utilized to suppress white noise interference, and on the basis, the transient characteristics of partial discharge signals are extracted by utilizing wavelet transformation and cross-lapped differential transformation, so that the identification accuracy of the partial discharge signals in a complex noise environment is improved; the cross-lapped differential transformation is suitable for extracting high-frequency transient signals, and detecting partial discharge signals by combining wavelet transformation and cross-lapped differential transformation, can effectively inhibit complex noise interference, and has higher detection sensitivity; the cross-over differential transformation algorithm is simple, and the efficiency of the detection method is greatly improved; in addition, the criterion adopts a self-adaptive threshold, so that the detection sensitivity is ensured, and the detection reliability is improved.
Example two
Referring to fig. 2, fig. 2 is a schematic structural composition diagram of a power cable partial discharge signal detection device in an embodiment of the present invention.
As shown in fig. 2, a power cable partial discharge signal detecting apparatus, the apparatus comprising:
the wavelet transform module 11: the system is used for performing wavelet transformation on the collected partial discharge signals of the power cable based on the selected wavelet basis and the decomposition layer number to obtain corresponding wavelet coefficients;
the threshold processing module 12: the method comprises the steps that a threshold value is obtained through calculation based on an algorithm selected by the threshold value, and threshold value processing is carried out on the wavelet coefficient based on the threshold value to obtain an estimated wavelet coefficient;
a reconstruction module 13: the wavelet coefficient is used for reconstructing the estimated wavelet coefficient to obtain an estimated value of the partial discharge signal of the power cable;
cross-lapped differential transform module 14: the device is used for carrying out H-order cross overlapping differential transformation on the estimated value of the power cable partial discharge signal and obtaining the transformation quantity of the power cable partial discharge signal through calculation;
the calculation module 15: and the average energy value and the maximum value of the transformation quantity of the partial discharge signal of the power cable are obtained through calculation based on the selected data window.
Specifically, the working principle of the device related function module according to the embodiment of the present invention may refer to the description related to the first method embodiment, and is not described herein again.
In the implementation of the invention, a threshold denoising method of wavelet transformation is utilized to suppress white noise interference, and on the basis, the transient characteristics of partial discharge signals are extracted by utilizing wavelet transformation and cross-lapped differential transformation, so that the identification accuracy of the partial discharge signals in a complex noise environment is improved; the cross-lapped differential transformation is suitable for extracting high-frequency transient signals, and detecting partial discharge signals by combining wavelet transformation and cross-lapped differential transformation, can effectively inhibit complex noise interference, and has higher detection sensitivity; the cross-over differential transformation algorithm is simple, and the efficiency of the detection method is greatly improved; in addition, the criterion adopts a self-adaptive threshold, so that the detection sensitivity is ensured, and the detection reliability is improved.
EXAMPLE III
In the specific implementation, the sampling frequency is set to be 100MHz, a white noise signal and sinusoidal signals with different frequencies are superposed in a partial discharge signal to simulate the mixing of white noise and periodic narrow-band interference in the discharge signal, and a single-exponential and double-exponential damped oscillation model is adopted in the discharge signal without noise pollution; in addition, the signal-to-noise ratio snr of the white noise signal is 5; the mathematical expression for the periodic narrowband interfering signal is:
Figure BDA0002512825300000081
wherein, n is 6, fj100kHz, 200kHz, 500kHz, 1MHz, 3M Hz and 4.5M Hz respectively;
as shown in fig. 3, the specific implementation process is as follows:
the first step, selecting wavelet base db6 and decomposing layer number J as 4, then making wavelet transform on the collected signal f (n) to obtain corresponding wavelet coefficient Wj,k
Secondly, calculating a threshold value T by adopting a general threshold value selection methodvThe method specifically comprises the following steps:
Figure BDA0002512825300000082
wherein σnIs the standard deviation of the noise, and N is the length of the signal, i.e. the number of sampling points;
thirdly, soft threshold value processing is carried out on the wavelet coefficient obtained by decomposition to obtain an estimated wavelet coefficient
Figure BDA0002512825300000083
Step four, for wavelet coefficient
Figure BDA0002512825300000084
Reconstructing to obtain the estimated value of the original signal
Figure BDA0002512825300000085
Fifthly, selecting the order H of the cross overlapping differential transformation to be 6, and calculating the transformation coefficient array (c)j)6The specific calculation method comprises the following steps:
a. order (c)j)1=[1 -1 0 0 0 0 0],(c′j)1Representative array (c)j)1The new array resulting after a bit shift to the left is then (c'j)1=[0 1 -1 0 0 0 0]J 1, …, 7 represents the sequence deficiency of the elements in the array;
b. then the transform coefficients corresponding to H ═ 6 can be calculated as: (c)j)6=[1 -6 15 -20 15 -6 1];
Sixth, using the obtained coefficient (c)j)6To the signal
Figure BDA0002512825300000086
Performing 6-order cross overlap difference conversion to obtain conversion quantity S6.f(n), the specific calculation formula is as follows:
Figure BDA0002512825300000087
seventhly, selecting a time window T equal to 1us corresponding to 100 sampling points P, and setting the discharge frequency q equal to 0;
calculating S6.fThe average energy value of (n) is:
Figure BDA0002512825300000088
calculation of E6Maximum value E of (n)6(nmax) Wherein n ismaxThe sampling point corresponding to the maximum value is satisfied with E6(nmax)>E6.settingThen, the partial discharge signal is detected, the number of discharges q is q +1, and n is recordedq.maxAnd a corresponding maximum value E6(nq.max);
Wherein the fixed value adopts a self-adaptive threshold:
a. before the discharge signal is not detected (q is 0), take E6.settins=Krel.1max(E6.normal),E6.normalTaking the reliability coefficient K for the average energy value of the cable under the normal conditionrel.11.5-2.5, in this example Krel.1=2;
b. If a discharge signal (q > 0) is detected, then fetch
Figure BDA0002512825300000091
Coefficient of reliability Krel.21.5-2.5, in this example Krel.2=2。
In the implementation of the invention, a threshold denoising method of wavelet transformation is utilized to suppress white noise interference, and on the basis, the transient characteristics of partial discharge signals are extracted by utilizing wavelet transformation and cross-lapped differential transformation, so that the identification accuracy of the partial discharge signals in a complex noise environment is improved; the cross-lapped differential transformation is suitable for extracting high-frequency transient signals, and detecting partial discharge signals by combining wavelet transformation and cross-lapped differential transformation, can effectively inhibit complex noise interference, and has higher detection sensitivity; the cross-over differential transformation algorithm is simple, and the efficiency of the detection method is greatly improved; in addition, the criterion adopts a self-adaptive threshold, so that the detection sensitivity is ensured, and the detection reliability is improved.
Those skilled in the art will appreciate that all or part of the steps in the methods of the above embodiments may be implemented by associated hardware instructed by a program, which may be stored in a computer-readable storage medium, and the storage medium may include: read-only Memory (ROM), Random Access Memory (RAM), magnetic or optical disks, and the like.
In addition, the above detailed description is provided for the method and apparatus for detecting a partial discharge signal of a power cable according to the embodiments of the present invention, and a specific example should be used herein to explain the principle and the implementation of the present invention, and the description of the above embodiment is only used to help understanding the method and the core idea of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, there may be variations in the specific embodiments and the application scope, and in summary, the content of the present specification should not be construed as a limitation to the present invention.

Claims (8)

1. A power cable partial discharge signal detection method, the method comprising:
performing wavelet transformation on the collected partial discharge signals of the power cable based on the selected wavelet basis and the decomposition layer number to obtain corresponding wavelet coefficients;
obtaining a threshold value through calculation based on an algorithm selected by the threshold value, and performing threshold value processing on the wavelet coefficient based on the threshold value to obtain an estimated wavelet coefficient;
reconstructing the estimated wavelet coefficient to obtain an estimated value of the partial discharge signal of the power cable;
carrying out H-order cross overlapping differential transformation on the estimated value of the power cable partial discharge signal, and calculating to obtain the transformation quantity of the power cable partial discharge signal;
and calculating to obtain the average energy value and the maximum value of the transformation quantity of the partial discharge signal of the power cable based on the selected data window.
2. The method for detecting a partial discharge signal of a power cable according to claim 1, wherein the algorithm selected based on the threshold is calculated to obtain the threshold, and the specific calculation formula is as follows:
Figure FDA0002512825290000011
wherein, TvRepresents the threshold value, σnRepresenting the standard deviation of the noise and N the length of the signal.
3. The method for detecting a partial discharge signal of a power cable according to claim 1, wherein said performing H-order cross-lapped differential transformation on the estimated value of the partial discharge signal of the power cable and calculating a transformation amount of the partial discharge signal of the power cable comprises:
obtaining a transformation coefficient array through calculation based on the selected order H of the cross overlapping differential transformation;
and performing H-order cross overlapping differential transformation on the estimated value of the partial discharge signal of the power cable based on the transformation coefficient array, and calculating to obtain the transformation quantity of the partial discharge signal of the power cable.
4. The method of claim 3, wherein the obtaining the transform coefficient array by calculation based on the selected order H of the cross-lapped differential transform comprises:
the transformation coefficient array of the H order is translated to the left by one bit to obtain a new array;
and calculating to obtain the transform coefficient corresponding to the H order based on the new array.
5. The method according to claim 3, wherein the conversion amount of the partial discharge signal of the power cable obtained by performing H-order cross-lapped differential conversion on the estimated value of the partial discharge signal of the power cable based on the conversion coefficient array is calculated by using a specific calculation formula as follows:
Figure FDA0002512825290000021
wherein S isH.f(n) represents the amount of conversion of the partial discharge signal of the power cable, and (c)j)HAn array of transform coefficients is represented,
Figure FDA0002512825290000022
represents an estimate of the power cable partial discharge signal, H represents the order, and j represents the number of decomposition layers.
6. The power cable partial discharge signal detection method of claim 1, wherein the calculating the mean energy value and the maximum value of the transformation amount of the power cable partial discharge signal based on the selected data window further comprises: and if the maximum value is larger than the fixed value, detecting that the local discharge signal of the power cable is detected.
7. The power cable partial discharge signal detection method of claim 6, wherein the constant value is set by an adaptive threshold, comprising:
if the discharge signal is not detected, the fixed value E6.setting=Krel.1max(E6.normal) In which E6.normalIs the average energy value of the cable under normal conditions, Krel.1A reliability coefficient greater than 1;
if a discharge signal is detected, the value is fixed
Figure FDA0002512825290000023
Wherein EH(nmax) Is a maximum value, Krel.2A reliability factor greater than 1.
8. A power cable partial discharge signal detection apparatus, the apparatus comprising:
a wavelet transformation module: the system is used for performing wavelet transformation on the collected partial discharge signals of the power cable based on the selected wavelet basis and the decomposition layer number to obtain corresponding wavelet coefficients;
a threshold processing module: the method comprises the steps that a threshold value is obtained through calculation based on an algorithm selected by the threshold value, and threshold value processing is carried out on the wavelet coefficient based on the threshold value to obtain an estimated wavelet coefficient;
a reconstruction module: the wavelet coefficient is used for reconstructing the estimated wavelet coefficient to obtain an estimated value of the partial discharge signal of the power cable;
a cross-lapped differential transform module: the device is used for carrying out H-order cross overlapping differential transformation on the estimated value of the power cable partial discharge signal and obtaining the transformation quantity of the power cable partial discharge signal through calculation;
a calculation module: and the average energy value and the maximum value of the transformation quantity of the partial discharge signal of the power cable are obtained through calculation based on the selected data window.
CN202010466552.7A 2020-05-28 2020-05-28 Power cable partial discharge signal detection method and device Active CN111781439B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010466552.7A CN111781439B (en) 2020-05-28 2020-05-28 Power cable partial discharge signal detection method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010466552.7A CN111781439B (en) 2020-05-28 2020-05-28 Power cable partial discharge signal detection method and device

Publications (2)

Publication Number Publication Date
CN111781439A true CN111781439A (en) 2020-10-16
CN111781439B CN111781439B (en) 2023-04-14

Family

ID=72754382

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010466552.7A Active CN111781439B (en) 2020-05-28 2020-05-28 Power cable partial discharge signal detection method and device

Country Status (1)

Country Link
CN (1) CN111781439B (en)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1853932A2 (en) * 2005-03-04 2007-11-14 University Of Strathclyde Detecting partial discharge in high voltage cables
CN103338023A (en) * 2013-03-19 2013-10-02 甘肃省电力公司检修公司 Filter eliminating external noises irrelevant to partial discharge during high-frequency partial discharge measurement
CN103336226A (en) * 2013-05-08 2013-10-02 清华大学 Identification method of various partial discharge source types in gas insulated substation (GIS)
CN104407270A (en) * 2014-12-04 2015-03-11 国家电网公司 Online fault monitoring device for cable connector in 10-35kV power distribution network and method for evaluating system state
CN104459498A (en) * 2014-12-22 2015-03-25 国网上海市电力公司 Method for suppressing interference on partial discharge test of ultra-high voltage AC transformer
CN204256093U (en) * 2014-12-04 2015-04-08 国家电网公司 The On-line Fault monitoring device of cable splice in a kind of 10 ~ 35kV power distribution network
CN106093724A (en) * 2016-06-20 2016-11-09 山东大学 A kind of shelf depreciation Wavelet noise-eliminating method based on hybrid particle swarm
CN107728018A (en) * 2017-09-20 2018-02-23 国网浙江省电力公司绍兴供电公司 A kind of noise-reduction method of power cable scene local discharge signal
CN108983051A (en) * 2018-07-25 2018-12-11 武汉科技大学 Based on the synchronous shelf depreciation kind identification method for squeezing wavelet transformation
CN109239554A (en) * 2018-09-28 2019-01-18 山东康威通信技术股份有限公司 A kind of denoising of power cable partial discharge signal and useful signal extracting method and system
CN110783896A (en) * 2019-09-25 2020-02-11 昆明理工大学 Distribution network single-phase grounding protection method based on weak fault active and passive joint detection
CN111046836A (en) * 2019-12-24 2020-04-21 杭州电力设备制造有限公司 Method, system, equipment and storage medium for filtering, denoising and analyzing partial discharge signal

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1853932A2 (en) * 2005-03-04 2007-11-14 University Of Strathclyde Detecting partial discharge in high voltage cables
US20090012727A1 (en) * 2005-03-04 2009-01-08 Wah Hoon Siew Detecting Partial Discharge in High Voltage Cables
CN103338023A (en) * 2013-03-19 2013-10-02 甘肃省电力公司检修公司 Filter eliminating external noises irrelevant to partial discharge during high-frequency partial discharge measurement
CN103336226A (en) * 2013-05-08 2013-10-02 清华大学 Identification method of various partial discharge source types in gas insulated substation (GIS)
CN204256093U (en) * 2014-12-04 2015-04-08 国家电网公司 The On-line Fault monitoring device of cable splice in a kind of 10 ~ 35kV power distribution network
CN104407270A (en) * 2014-12-04 2015-03-11 国家电网公司 Online fault monitoring device for cable connector in 10-35kV power distribution network and method for evaluating system state
CN104459498A (en) * 2014-12-22 2015-03-25 国网上海市电力公司 Method for suppressing interference on partial discharge test of ultra-high voltage AC transformer
CN106093724A (en) * 2016-06-20 2016-11-09 山东大学 A kind of shelf depreciation Wavelet noise-eliminating method based on hybrid particle swarm
CN107728018A (en) * 2017-09-20 2018-02-23 国网浙江省电力公司绍兴供电公司 A kind of noise-reduction method of power cable scene local discharge signal
CN108983051A (en) * 2018-07-25 2018-12-11 武汉科技大学 Based on the synchronous shelf depreciation kind identification method for squeezing wavelet transformation
CN109239554A (en) * 2018-09-28 2019-01-18 山东康威通信技术股份有限公司 A kind of denoising of power cable partial discharge signal and useful signal extracting method and system
CN110783896A (en) * 2019-09-25 2020-02-11 昆明理工大学 Distribution network single-phase grounding protection method based on weak fault active and passive joint detection
CN111046836A (en) * 2019-12-24 2020-04-21 杭州电力设备制造有限公司 Method, system, equipment and storage medium for filtering, denoising and analyzing partial discharge signal

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
JUNHYUCK SEO: "Signal extraction methods for online partial discharge (PD) measurement of power transformers", 《2015 IEEE 11TH INTERNATIONAL CONFERENCE ON THE PROPERTIES AND APPLICATIONS OF DIELECTRIC MATERIALS (ICPADM)》 *
唐炬等: "SF6局部放电分解组分光声检测信号交叉响应处理技术", 《高电压技术》 *
尚海昆: "电力变压器局部放电信号的特征提取与模式识别方法研究", 《中国博士学位论文全文数据库 工程科技II辑》 *
顾媛媛: "电力电缆局部放电信号理论及处理研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 *

Also Published As

Publication number Publication date
CN111781439B (en) 2023-04-14

Similar Documents

Publication Publication Date Title
Song et al. Second generation wavelet transform for data denoising in PD measurement
Cunha et al. A new wavelet selection method for partial discharge denoising
CN110501631B (en) Online intermittent fault detection and diagnosis method
CN103018629A (en) Method for analyzing power system fault recording data based on Marla algorithm
Chang et al. Denoising of partial discharge signals in wavelet packets domain
CN109581270B (en) Method for rapidly identifying CVT fault based on wavelet multi-scale transformation
CN111077386A (en) Early fault signal noise reduction method for electrical equipment
CN110991481A (en) High-voltage shunt reactor internal loosening fault diagnosis method based on cross wavelet transformation
CN111239565B (en) Oil-filled casing partial discharge pulse signal processing method and system based on layered denoising model
CN108089100A (en) The detection method of small current neutral grounding system arc light resistance earth fault
Yang et al. Partial discharge signal extraction method based on EDSSV and low rank RBF neural network
CN110161351B (en) Transformer winding fault test system and diagnosis method under oscillatory wave
CN108594156B (en) Improved current transformer saturation characteristic identification method
CN114038470A (en) Method, device, equipment and medium for extracting voiceprint features of transformer signals
CN111044865B (en) Partial discharge wavelet denoising method for direct-current gas insulation electrical equipment
CN113642417A (en) Improved wavelet algorithm-based denoising method for partial discharge signals of insulated overhead conductor
CN111781439B (en) Power cable partial discharge signal detection method and device
CN109324268B (en) Power distribution network early fault detection method and device based on Bayesian inference
CN114200232B (en) Method and system for detecting fault traveling wave head of power transmission line
CN117009870A (en) Pump cavitation state identification method for frequency domain improved SDP diagram
CN114779028B (en) Generator partial discharge online monitoring device and monitoring method
CN114462452B (en) Asynchronous motor rotor broken bar fault diagnosis method using successive variable mode decomposition algorithm
CN115774199A (en) Fault diagnosis method and system for new energy battery, electronic device and storage medium
CN107797025B (en) Fault positioning method and device for power system
Madhu et al. A novel algorithm for denoising of simulated partial discharge signals using adaptive wavelet thresholding methods

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