CN113406453B - PRPD/PRPS map data processing method and detection device based on MCU - Google Patents

PRPD/PRPS map data processing method and detection device based on MCU Download PDF

Info

Publication number
CN113406453B
CN113406453B CN202110735971.0A CN202110735971A CN113406453B CN 113406453 B CN113406453 B CN 113406453B CN 202110735971 A CN202110735971 A CN 202110735971A CN 113406453 B CN113406453 B CN 113406453B
Authority
CN
China
Prior art keywords
loc
signal
array
seg
mcu
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.)
Active
Application number
CN202110735971.0A
Other languages
Chinese (zh)
Other versions
CN113406453A (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.)
Pingdingshan University
Original Assignee
Pingdingshan University
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 Pingdingshan University filed Critical Pingdingshan University
Priority to CN202110735971.0A priority Critical patent/CN113406453B/en
Publication of CN113406453A publication Critical patent/CN113406453A/en
Application granted granted Critical
Publication of CN113406453B publication Critical patent/CN113406453B/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
    • G01R31/12Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16Spectrum analysis; Fourier analysis
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/04Programme control other than numerical control, i.e. in sequence controllers or logic controllers
    • G05B19/042Programme control other than numerical control, i.e. in sequence controllers or logic controllers using digital processors
    • G05B19/0423Input/output

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Mathematical Physics (AREA)
  • Testing Relating To Insulation (AREA)

Abstract

The invention discloses a PRPD/PRPS map data processing method and a detection device based on MCU, wherein the operation of the processing method is only an addition and subtraction algorithm, the operation amount is small, the method can be realized on MCU with limited resources, and the detection device is lighter and more convenient to carry and operate compared with the traditional partial discharge detection equipment; the MCU processes the digital power frequency signal and the digital discharge signal, and the amplitude and the phase of the obtained pulse signal are uploaded to the upper computer, and the upper computer does not need to carry out a large amount of data operation, so that the upper computer can adopt terminals with lower operation performance such as a PDA or a smart phone, and the like, and the complexity of the used device is further reduced; the detection device is used for realizing the processing method, the device comprises an MCU, the MCU comprises a RAM and a serial port, the RAM is connected with an AD module through a data bus, and the MCU is connected with an upper computer through the serial port.

Description

PRPD/PRPS map data processing method and detection device based on MCU
Technical Field
The invention relates to the technical field of high-voltage electrical equipment insulation detection, in particular to a PRPD/PRPS map data processing method and a detection device based on MCU.
Background
When the high-voltage electrical equipment has partial discharge faults, various signals such as high frequency, ultrasonic, transient ground wave, ultrahigh frequency and the like are generated, and whether the partial discharge occurs or not can be effectively judged by detecting the signals. The common partial discharge detection equipment mainly comprises a high-frequency partial discharge detector, an ultrasonic partial discharge detector, a transient-state ground wave partial discharge detector and an ultrahigh-frequency partial discharge detector. The basic structure of the instrument has the same structure and mainly comprises a partial discharge sensor module, a power frequency signal sensor module, a signal filtering and amplifying conditioning module, a signal acquisition module, a signal analysis module and a result presentation module. The current mainstream partial discharge analysis method is to process the detected signal, draw a PRPD (phase resolved partial discharge)/PRPS (phase resolved pulse sequence) map related to the power frequency period, and perform the basis of partial discharge diagnosis through the characteristics of the map. The drawing of the partial discharge PRPD/PRPS map requires the following two key parameters, namely the amplitude of the partial discharge pulse signal and the phase relation between the pulse signal generation time and the power frequency signal. The signal analysis module in the partial discharge detector can extract the two key parameters from the detected partial discharge signal and power frequency signal. Because the frequency of partial discharge is high, the frequency of the signal frequency band is wide and the frequency is high, the data volume of the detected original signal is large, meanwhile, the partial discharge is closely related to the power frequency period, a large amount of calculation is often required to obtain the two types of parameters, and the signal analysis module of the traditional detection instrument is usually finished in a computer to obtain a more accurate partial discharge PRPD/PRPS map, so that the detection instrument is usually heavy and the field operation is relatively complex. The detection instruments play a good role in traditional power failure detection, but when live inspection of substation total station equipment is carried out, a plurality of carrying and operation are inconvenient, and the main reason is that hardware resources are excessively occupied by a data analysis module.
Disclosure of Invention
In order to solve the defects in the prior art, the invention provides the PRPD/PRPS map data processing method and the detection device based on the MCU, wherein the operation of the processing method is only an addition and subtraction algorithm, the operation amount is small, the operation difficulty is low, the processing process is carried out on a single chip microcomputer, and compared with the traditional detection instrument, the detection device is lighter and convenient to carry and operate.
In order to achieve the above purpose, the present invention provides the following technical solutions: a PRPD/PRPS map data processing method based on MCU, MCU electric connection has AD module and host computer, MCU includes RAM, its characterized in that: the processing method comprises the following steps:
s1, collecting partial discharge signals and power frequency signals;
s2, converting the partial discharge signal and the power frequency signal into a digital discharge signal and a digital power frequency signal, and storing the digital discharge signal and the digital power frequency signal in a RAM;
s3, extracting a digital power frequency signal, and calculating the 0-degree position of the digital power frequency signal;
s4, extracting a digital discharge signal, and acquiring the position and the amplitude of a pulse signal in the digital discharge signal;
s5, adjusting the position of the pulse signal according to the 0-degree position of the digital power frequency signal;
and S6, uploading the amplitude of the pulse signal and the position of the adjusted pulse signal to an upper computer for map display.
As a further optimization of the above-described MCU-based PRPD/PRPS profile data processing method: the specific steps of S2 are as follows:
s201, setting the sampling rate of AD module as f S/S, the sampling bit number as D bit, setting the empty array MER_A and the empty array MER_B in RAM;
s202, an AD module converts a partial discharge signal and a power frequency signal into a digital discharge signal x [ n ] and a digital power frequency signal g [ n ];
s203, storing the digital discharging signal x [ n ] and the digital power frequency signal g [ n ] in the empty array MER_A [ ], storing the digital discharging signal x [ n ] and the digital power frequency signal g [ n ] in the empty array MER_B [ ] after the empty array MER_A [ ] is full, and executing S3 after any one of the empty array MER_A [ ] and the empty array MER_B [ ] is full.
As a further optimization of the above-described MCU-based PRPD/PRPS profile data processing method: the specific steps of S3 are as follows:
s301, reading a digital power frequency signal g [ n ], and setting the signal time to be 40ms;
s302, digital power frequency signal g [ n ]]Is equally divided into signals g a [n]Sum signal g b [n]Signal g a [n]Divided into signal segments g a_1 [n]Signal section g a_2 [n]Signal section g a_3 [n]Signal section g a_4 [n]The sampling point number of each signal segment is
Figure BDA0003141703830000021
S303, orthogonal components p and q are calculated:
Figure BDA0003141703830000022
Figure BDA0003141703830000023
s304, judging signal g a [n]Middle 0 deg. position loc a0 Is defined by the quadrants of:
p and q are both positive, and the quadrant at the 0 position is 1; p is positive, q is negative, and the quadrant at the 0 DEG position is 2; p is negative, q is positive, and the quadrant at the 0 position is 3; when p and q are both negative, the quadrant at the 0 position is 4.
S305, calculating 0-degree position loc a0
When quadrant is 1, signal g a [n]The position of the minimum value is subtracted
Figure BDA0003141703830000024
The sampling points are 0 degree position loc a0 The method comprises the steps of carrying out a first treatment on the surface of the When quadrant is 2, signal g a [n]The position of the minimum minus +.>
Figure BDA0003141703830000025
The sampling points are 0 degree position loc a0 The method comprises the steps of carrying out a first treatment on the surface of the When quadrant is 3, signal g a [n]The position of the minimum value is added with +.>
Figure BDA0003141703830000026
The sampling points are 0 degree position loc a0 The method comprises the steps of carrying out a first treatment on the surface of the Quadrant 4, signal g a [n]The position of the minimum value is added with +.>
Figure BDA0003141703830000031
The sampling points are 0 degree position loc a0
S306, calculating signal g b [n]In 0 deg. position loc b0
Figure BDA0003141703830000032
As a further optimization of the above-described MCU-based PRPD/PRPS profile data processing method: the specific steps of S4 are as follows:
s401, reading digital discharge signal x [ n ]]Setting a signal threshold value x Thr And signal time, signal time is 40ms, and an array amp is built in RAM]And array loc_pha [];
S402, setting the width of the rectangular data window as L, L as even number, and setting the digital discharge signal x [ n ]]Starting at x t]From x [ t ]]Starting to continuously read L sampling points, and recording the signal section where the L sampling points are positioned as x L [n]Wherein t=0 to (n-L);
s403, statistics x L [n]X of the amplitude maximum of (x) L_max And its position loc:
if x L_max ≤x Thr Then, t is added with 1, S402 is executed, and t is the number of repeated execution times; if x L_max >x Thr And x is L_max The amplitude is larger than the amplitude of L/2 sampling points before loc and L/2 sampling points after loc, then x is L [n]The amplitude x of the pulse signal is recorded as a pulse signal L_max And position loc;
s404, search array loc_pha []If the position loc has been stored, the amplitude x obtained in S403 is then calculated L_max And position loc discard, otherwise amplitude x L_max Store array amp]The position loc is stored in an array loc_pha []I represents the obtained pulse ordinal index; s405, adding 1 to t, executing S402 until t=n-L, array amp []And array loc_pha []Respectively store x [ n ]]The amplitude and the position of the pulse signals contained in the pulse signal processing circuit are set to be I.
As a further optimization of the above-described MCU-based PRPD/PRPS profile data processing method: the specific steps of S5 are as follows:
s501, establishing an array A_seg [ ], an array B_seg [ ], an array C_seg [ ] and an array C_seg_plus [ ];
s502, extracting the position loc_phai [ I ] of the pulse signals stored in the array loc_phaj, wherein i=0 to (I-1);
s503, loc_phai]And loc a0 And loc b0 Comparison was performed:
if loc_phai]Less than loc a0 Then calculate
Figure BDA0003141703830000033
And (loc) a0 -loc_pha[i]) And assign the difference to loc_phai]Amp [ i ]]And assigned loc_phai]Sequentially storing array A_seg [ []The method comprises the steps of carrying out a first treatment on the surface of the If loc_phai]Greater than loc a0 And is less than loc b0 Then calculate loc_phai]And loc a0 A difference value, which is assigned to loc_phai]Amp [ i ]]And assigned loc_phai]Sequentially storing an array B_seg [ into ] a memory]The method comprises the steps of carrying out a first treatment on the surface of the If loc_phai]Greater than loc b0 Then calculate loc_phai]And loc b0 Is assigned to loc_phai [ i ]]Amp [ i ]]And assigned loc_phai]Sequentially storing the array C_seg [ into ] a memory];
S504, combining the values in the array A_seg [ ] and the array C_seg_plus [ ] and then storing the values in the array A_seg [ ] and simultaneously clearing the array C_seg_plus [ ];
s505, storing the array C_seg [ ] into the array C_seg_plus [ ].
As a further optimization of the above-described MCU-based PRPD/PRPS profile data processing method: the specific method of S6 is as follows: uploading the positions and the amplitudes of the pulse signals stored in the arrays A_seg [ ] and B_seg [ ] to an upper computer, clearing the arrays A_seg [ ] and B_seg [ ] after uploading is finished, and returning to S2 for repeated execution.
The PRPD/PRPS map data detection device based on the MCU is used for the PRPD/PRPS map data processing method based on the MCU, and comprises the MCU, wherein the MCU comprises a RAM and a serial port, the RAM is connected with an AD module through a data bus, and the MCU is connected with an upper computer through the serial port.
The beneficial effects are that: the PRPD/PRPS map data processing method and the detection device based on the MCU provided by the invention have the advantages that the processing method is operated by an addition and subtraction algorithm, the operation amount is small, the detection device can be realized on the MCU with limited resources, and compared with the traditional partial discharge detection equipment, the detection device is lighter and convenient to carry and operate; MCU processes the digital power frequency signal and the digital discharge signal, and the amplitude and the phase of the obtained pulse signal are uploaded to the upper computer, and the upper computer does not need to carry out a large amount of data operation, so that the upper computer can adopt terminals with lower operation performance such as PDA or smart phone, and the like, and the complexity of the detection device is further reduced.
Drawings
FIG. 1 is a flow chart of a data processing method of the present invention;
FIG. 2 is a block diagram of the detection device of the present invention;
fig. 3 is a graph of sample points versus amplitude.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Referring to fig. 1 to 3, an MCU-based PRPD/PRPS map data processing method, the MCU is electrically connected with an AD module and an upper computer, and the MCU includes a RAM, and is characterized in that: the processing method comprises the following steps:
s1, collecting partial discharge signals and power frequency signals.
In the step S1, the local discharge signal and the power frequency signal can be collected as a local discharge equipment sensor, an ultrahigh frequency antenna, an ultrasonic sensor and the like, the collected local discharge signal and the power frequency signal are conditioned, the conditioned local discharge signal is a signal after envelope detection and frequency reduction, the conditioned power frequency signal is a power frequency signal after band-pass filtering, and the AD module receives the conditioned local discharge signal and the power frequency signal and executes the step S2.
In order to meet the requirements of the processing method, the sampling rate of the AD module is up to 100kS/s and above, the sampling bit number is up to 12bit and above, the MCU working frequency is up to 80MHz and above, and the MCU has a DMA (direct memory access) function.
S2, converting the partial discharge signals and the power frequency signals into digital discharge signals and digital power frequency signals, and storing the digital discharge signals and the digital power frequency signals in a RAM (memory).
The specific steps of S2 are as follows:
s201, setting the sampling rate of AD module as f S/S, sampling bit number as D bit, setting up empty array MER_A and empty array MER_B in RAM.
S202, the AD module converts the partial discharge signal and the power frequency signal into a digital discharge signal x [ n ] and a digital power frequency signal g [ n ].
S203, storing the digital discharging signal x [ n ] and the digital power frequency signal g [ n ] in the empty array MER_A [ ], storing the digital discharging signal x [ n ] and the digital power frequency signal g [ n ] in the empty array MER_B [ ] after the empty array MER_A [ ] is full, and executing S3 after any one of the empty array MER_A [ ] and the empty array MER_B [ ] is full.
S2, empty array MER_A [ []And empty array MER_B []Is of the size of
Figure BDA0003141703830000051
Byte (+)>
Figure BDA0003141703830000052
To round up the operator). The AD module converts the partial discharge signal and the power frequency signal into digital discharge signals x [ n ]]And digital power frequency signal g [ n ]]Then, the DMA function of MCU sequentially outputs to the empty array MER_A [ MER_A ]]And empty array MER_B []And storing, namely storing the method as S203.
S3, extracting a digital power frequency signal, and calculating the 0-degree position of the digital power frequency signal.
The specific steps of S3 are as follows:
s301, reading a digital power frequency signal g [ n ], and setting the signal time to be 40ms.
S302, digital power frequency signal g [ n ]]Is equally divided into signals g a [n]Sum signal g b [n]Signal g a [n]Divided into signal segments g a_1 [n]Signal section g a_2 [n]Signal section g a_3 [n]Signal section g a_4 [n]The sampling point number of each signal segment is
Figure BDA0003141703830000053
Figure BDA0003141703830000054
S303, orthogonal components p and q are calculated:
Figure BDA0003141703830000055
Figure BDA0003141703830000056
s304, judging signal g a [n]Middle 0 deg. position loc a0 Is defined by the quadrants of:
p and q are both positive, and the quadrant at the 0 position is 1; p is positive, q is negative, and the quadrant at the 0 DEG position is 2; p is negative, q is positive, and the quadrant at the 0 position is 3; when p and q are both negative, the quadrant at the 0 position is 4.
S305, calculating 0-degree position loc a0
When quadrant is 1, signal g a [n]The position of the minimum value is subtracted
Figure BDA0003141703830000057
The sampling points are 0 degree position loc a0 The method comprises the steps of carrying out a first treatment on the surface of the When quadrant is 2, signal g a [n]The position of the minimum minus +.>
Figure BDA0003141703830000058
The sampling points are 0 degree position loc a0 The method comprises the steps of carrying out a first treatment on the surface of the When quadrant is 3, signal g a [n]The position of the minimum value is added with +.>
Figure BDA0003141703830000061
The sampling points are 0 degree position loc a0 The method comprises the steps of carrying out a first treatment on the surface of the Quadrant 4, signal g a [n]The position of the minimum value is added with +.>
Figure BDA0003141703830000062
The sampling points are 0 degree position loc a0
S306, calculating signal g b [n]In 0 deg. position loc b0
Figure BDA0003141703830000063
And S3, processing the digital power frequency signal g [ n ]:
empty array MER_A []And empty array MER_B []After any one of the storage is full, the digital power frequency signal g [ n ] is read from the empty array of the storage is full]The signal time is set to 40ms, and the digital power frequency signal g [ n ] is based on the digital orthogonal algorithm]Is estimated at 0 deg. position. Digital power frequency signal g [ n ]]Is divided into two sections, respectively signal g a [n]Sum signal g b [n]Signal g a [n]Sum signal g b [n]Each containing 20ms of signal. Digital power frequency signal g [ n ]]Is divided into two parts:
for signal g a [n]0 of (2)The position of the degree, signal g a [n]Dividing into four segments, and marking the four signal segments as signal segment g a_1 [n]Signal section g a_2 [n]Signal section g a_3 [n]Signal section g a_4 [n]Calculating the signal g from the four signal segments a [n]After calculation, quadrant judgment is performed on the 0 ° position according to the method shown in table 1, and the 0 ° position is calculated according to the quadrant judgment as shown in S305.
For signal g b [n]0 DEG position of (2), signal g b [n]Without processing, the 0 ° position can be directly calculated as S306.
Table 1: digital power frequency signal 0 degree position quadrant judgment table
p + + - -
q + - + -
Quadrant with a plurality of quadrants 1 4 2 3
S4, extracting a digital discharge signal, and acquiring the position and the amplitude of a pulse signal in the digital discharge signal.
The specific steps of S4 are as follows:
s401, reading digital discharge signal x [ n ]]Setting a signal threshold value x Thr And signal time, signal time is 40ms, and an array amp is built in RAM]And array loc_pha []。
S402, setting the width of the rectangular data window as L, L as even number, and setting the digital discharge signal x [ n ]]Starting at x t]From x [ t ]]Starting to continuously read L sampling points, and recording the signal section where the L sampling points are positioned as x L [n]Wherein t=0 to (n-L).
S403, statistics x L [n]X of the amplitude maximum of (x) L_max And its position loc:
if x L_max ≤x Thr Then, t is added with 1, S402 is executed, and t is the number of repeated execution times; if x L_max >x Thr And x is L_max The amplitude is larger than the amplitude of L/2 sampling points before loc and L/2 sampling points after loc, then x is L [n]The amplitude x of the pulse signal is recorded as a pulse signal L_max And a location loc.
S404, search array loc_pha []If the position loc has been stored, the amplitude x obtained in S403 is then calculated L_max And position loc discard, otherwise amplitude x L_max Store array amp]The position loc is stored in an array loc_pha []I represents the obtained pulse ordinal index.
S405, adding 1 to t, executing S402 until t=n-L, storing the amplitude and position of the pulse signal contained in x [ t ] in the array amp [ and the array loc_phaj ] respectively, and setting the number of the pulse signals as I.
In S4, the digital discharge signal x [ n ] is processed:
reading digital discharge signals x [ n ] from a full-memory empty array]Setting a signal threshold value x Thr And signal time, which is also 40ms, while establishing an array amp in RAM]And array loc_pha []. Setting a rectangular data window with width L and digital discharge signal x [ n ]]Is marked as x t]In x [ t ]]Starting from a digital discharge signal x [ n ] with a data window]Performing frame selection, wherein the signal segment selected by the frame is marked as x L [n]Statistics of x L [n]Amplitude maximum x in (a) L_max And the position loc, if x L_max Both greater than the signal threshold x Thr And the amplitude value is larger than the amplitude value of L/2 sampling points before loc and L/2 sampling points after loc, then x is L [n]Is a pulse signal. Search array loc_pha []If the position loc is not stored, the amplitude x L_max Store array amp]The position loc is stored in an array loc_pha []. According to this procedure, from x [ t ]]Initially, digital discharge signals x [ n ] are sequentially discharged using a data window]And performing box selection, and moving one sampling point at a time until t=n-L.
S5, adjusting the position of the pulse signal according to the 0-degree position of the digital power frequency signal.
The specific steps of S5 are as follows:
s501, establishing an array A_seg [ ], an array B_seg [ ], an array C_seg [ ] and an array C_seg_plus [ ];
s502, extracting the position loc_phai [ I ] of the pulse signal stored in the array loc_phaj [ I ], wherein i=0 to (I-1).
S503, loc_phai]And loc a0 And loc b0 Comparison was performed:
if loc_phai]Less than loc a0 Then calculate
Figure BDA0003141703830000071
And (loc) a0 -loc_pha[i]) And assign the difference to loc_phai]Amp [ i ]]And assigned loc_phai]Sequentially storing array A_seg [ []The method comprises the steps of carrying out a first treatment on the surface of the If loc_phai]Greater than loc a0 And is less than loc b0 Then calculate loc_phai]And loc a0 A difference value, which is assigned to loc_phai]Amp [ i ]]And assigned loc_phai]Sequentially storing an array B_seg [ into ] a memory]The method comprises the steps of carrying out a first treatment on the surface of the If loc_phai]Greater than loc b0 Then calculate loc_phai]And loc b0 Is assigned to loc_phai [ i ]]Amp [ i ]]And assigned loc_phai]Sequentially storing the array C_seg [ into ] a memory];
S504, combining the values in the array A_seg [ ] and the array C_seg_plus [ ] and then storing the values in the array A_seg [ ] and simultaneously clearing the array C_seg_plus [ ];
s504, combining the values in the array A_seg [ ] and the array C_seg_plus [ ] and then storing the values in the array A_seg [ ] and simultaneously clearing the array C_seg_plus [ ];
s505, storing the array C_seg [ ] into the array C_seg_plus [ ].
S5, according to the digital power frequency signal g [ n ] obtained in S3]0 deg. position loc of (2) a0 And loc b0 For the position loc_phai [ i ] of the pulse signal obtained in S4]And (3) adjusting:
four arrays are built, respectively, array A_seg [ []Array B_seg []Array C_seg []And array C_seg_plus []Four arrays are used to segment the pulse signal in full period. Extracting loc_phai [ i ] obtained in S4]By loc_phai]And loc a0 And loc b0 Comparison was performed:
if loc_phai]Less than loc a0 Then calculate
Figure BDA0003141703830000081
And (loc) a0 -loc_pha[i]) And assign the difference to loc_phai]Amp [ i ]]And assigned loc_phai]Sequentially storing array A_seg [ []The method comprises the steps of carrying out a first treatment on the surface of the If loc_phai]Greater than loc a0 And is less than loc b0 Then calculate loc_phai]And loc a0 A difference value, which is assigned to loc_phai]Amp [ i ]]And assigned loc_phai]Sequentially storing an array B_seg [ into ] a memory]The method comprises the steps of carrying out a first treatment on the surface of the If loc_phai]Greater than loc b0 Then calculate loc_phai]And loc b0 Is assigned to loc_phai [ i ]]Amp [ i ]]And assigned loc_phai]Sequentially storing the array C_seg [ into ] a memory]。
After the comparing and assigning steps are completed, the values in C_seg_plus are combined to A_seg and C_seg_plus is emptied, and the values in C_seg are stored in C_seg_plus.
And S6, uploading the amplitude of the pulse signal and the position of the adjusted pulse signal to an upper computer for map display.
The specific method of S6 is as follows: uploading the positions and the amplitudes of the pulse signals stored in the arrays A_seg [ ] and B_seg [ ] to an upper computer, clearing the arrays A_seg [ ] and B_seg [ ] after uploading is finished, and returning to S2 for repeated execution.
Because the processing processes are all carried out on the MCU, and the operation is only addition and subtraction, the upper computer does not need to carry out the operation, so that terminals such as a PDA (palm computer) or a smart phone with lower operation performance can be adopted for the upper computer, and the complexity of the required device is reduced.
The utility model provides a PRPD/PRPS map data detection device based on MCU, is used for above-mentioned PRPD/PRPS map data processing method based on MCU, its characterized in that: the device comprises an MCU, wherein the MCU comprises a RAM and a serial port, the RAM is connected with an AD module through a data bus, and the MCU is connected with an upper computer through the serial port.
Working principle: the AD module adopts an AD analog-to-digital converter with the model of AD7606, the sampling rate is 100kS/s, the MCU adopts a singlechip with the model of STM32F407ZGT6, and the MCU has a DMA function. The AD module receives the conditioned partial discharge signals and the conditioned power frequency signals, converts the partial discharge signals and the power frequency signals into digital discharge signals and digital power frequency signals, the MCU reads the digital discharge signals and the digital power frequency signals in the AD module by using the DAM function, processes the digital discharge signals and the digital power frequency signals, and the amplitude and the position of the pulse signals obtained after the processing are transmitted to the upper computer through the serial port, and the upper computer displays a map according to the amplitude and the position of the pulse signals.
The present invention will be described in detail with reference to specific examples.
In order to meet the requirements of the processing method, an AD module adopts an AD analog-to-digital converter with the model of AD7606, and an MCU adopts a singlechip with the model of STM32F407ZGT 6.
FIG. 3 shows the data stored in the array MER_A [ ] for a partial discharge detection, the sampling rate is 100kS/s, and the array MER_A [ ] stores the power frequency signals of 4000 sampling points and the partial discharge signals of 4000 sampling points.
Processing the digital power frequency signal according to S3:
signal g a [n]Divided into signal segments g a_1 [n]Signal section g a_2 [n]Signal section g a_3 [n]Signal section g a_4 [n]According to S303, the calculated quadrature component p is 1633.15 and the quadrature component q is-940.17, the signal g a [n]Is in quadrant 4, signal g a [n]The minimum value of (2) is at 1334 th sample point, signal g a [n]0 deg. position loc of (2) a0 The method comprises the following steps: loc a0 =1348+500=1848,
According to S306, a calculated loc is obtained b0 =3848。
Processing the digital discharge signal according to S4:
setting the width l=10 of the rectangular data window, the signal threshold x Thr =0.015V, first extraction x [ n ]]Wherein the maximum amplitude value is 0.009V, less than the signal threshold x Thr . Extracting the digital discharge signal until 623 th extraction, at which time signal segment x L [n]X of the maximum value of (x) L_max First exceeding the signal threshold x Thr . After all the digital discharge signals were processed, 4 pulse signals were obtained in total, and the amplitude and position thereof were recorded as shown in table 2.
Table 2: amplitude and position of pulse signal obtained by searching
Numbering device 1 2 3 4
loc_pha[] 632 2175 2942 3708
amp[] 0.407 0.425 0.819 0.587
Adjusting the position of the pulse signal according to S5:
due to loc a0 =1848、loc b0 The position and amplitude of the first pulse signal are stored in the array A_seg =3848]The second, third and fourth pulses are stored in array B_seg [ []In array C_seg []Is a null array. Due to array C_seg_plus [ sic ]]The middle is also empty. Final array A_seg []And B_seg []The data stored in the memory are shown in table 3.
Table 3: data stored in A_seg [ ] and B_seg [ ]
A_seg[] 632 0.407
B_seg[] 2175 0.425 2942 0.819 3708 0.587
According to S6, the MCU uploads the data stored in the arrays A_seg [ ] and B_seg [ ] to the upper computer through the serial port, and the upper computer displays the PRPD/PRPS map.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims (5)

1. A PRPD/PRPS map data processing method based on MCU, MCU electric connection has AD module and host computer, MCU includes RAM, its characterized in that: the processing method comprises the following steps:
s1, collecting partial discharge signals and power frequency signals;
s2, converting the partial discharge signal and the power frequency signal into a digital discharge signal and a digital power frequency signal, and storing the digital discharge signal and the digital power frequency signal in a RAM;
the specific steps of S2 are as follows:
s201, setting the sampling rate of the AD module as fS/S, the sampling bit number as D bit, and setting a null array MER_A and a null array MER_B in the RAM;
s202, an AD module converts a partial discharge signal and a power frequency signal into a digital discharge signal x [ n ] and a digital power frequency signal g [ n ];
s203, storing the digital discharging signal x [ n ] and the digital power frequency signal g [ n ] in a null array MER_A [ ], storing the digital discharging signal x [ n ] and the digital power frequency signal g [ n ] in a null array MER_B [ ] after the null array MER_A [ ] is full, and executing S3 after any one of the null array MER_A [ ] and the null array MER_B [ ] is full; s3, extracting a digital power frequency signal, and calculating the 0-degree position of the digital power frequency signal;
the specific steps of S3 are as follows:
s301, reading a digital power frequency signal g [ n ], and setting the signal time to be 40ms;
s302, digital power frequency signal g [ n ]]Is equally divided into signals g a [n]Sum signal g b [n]Signal g a [n]Divided into signal segments g a_1 [n]Signal section g a_2 [n]Signal section g a_3 [n]Signal section g a_4 [n]The sampling point number of each signal segment is
Figure QLYQS_1
S303, orthogonal components p and q are calculated:
Figure QLYQS_2
s304, judging signal g a [n]Middle 0 deg. position loc a0 Is defined by the quadrants of:
p and q are both positive, and the quadrant at the 0 position is 1; p is positive, q is negative, and the quadrant at the 0 DEG position is 2; p is negative, q is positive, and the quadrant at the 0 position is 3; when p and q are both negative, the quadrant at the 0 DEG position is 4;
s305, calculating 0-degree position loc a0
When quadrant is 1, signal g a [n]The position of the minimum value is subtracted
Figure QLYQS_3
The sampling points are 0 degree position loc a0 The method comprises the steps of carrying out a first treatment on the surface of the When quadrant is 2, signal g a [n]The position of the minimum minus +.>
Figure QLYQS_4
The sampling points are 0 degree position loc a0 The method comprises the steps of carrying out a first treatment on the surface of the When quadrant is 3, signal g a [n]The position of the minimum value is added with +.>
Figure QLYQS_5
Figure QLYQS_6
The sampling points are 0 degree position loc a0 The method comprises the steps of carrying out a first treatment on the surface of the Quadrant 4, signal g a [n]The position of the minimum value is added with +.>
Figure QLYQS_7
Figure QLYQS_8
The sampling points are 0 degree position loc a0;
S306, calculating signal g b [n]In 0 deg. position loc b0
Figure QLYQS_9
S4, extracting a digital discharge signal, and acquiring the position and the amplitude of a pulse signal in the digital discharge signal;
s5, adjusting the position of the pulse signal according to the 0-degree position of the digital power frequency signal;
and S6, uploading the amplitude of the pulse signal and the position of the adjusted pulse signal to an upper computer for map display.
2. The MCU-based PRPD/PRPS pattern data processing method of claim 1, wherein: the specific steps of S4 are as follows:
s401, reading digital discharge signal x [ n ]]Setting a signal threshold value x Thr And signal time, signal time is 40ms, and an array amp is built in RAM]And array loc_pha [];
S402, setting the width of the rectangular data window as L, L as even number, and setting the digital discharge signal x [ n ]]Starting at x t]From x [ t ]]Initial connectionReading L sampling points continuously, and marking the signal section where the L sampling points are positioned as x L [n]Wherein t=0 to (n-L);
s403, statistics x L [n]X of the amplitude maximum of (x) L_max And its position loc:
if x L_max ≤x Thr Then, t is added with 1, S402 is executed, and t is the number of repeated execution times; if x L_max >x Thr And x is L_max The amplitude is larger than the amplitude of L/2 sampling points before loc and L/2 sampling points after loc, then x is L [n]The amplitude x of the pulse signal is recorded as a pulse signal L_max And position loc;
s404, search array loc_pha []If the position loc has been stored, the amplitude x obtained in S403 is then calculated L_max And position loc discard, otherwise amplitude x L_max Store array amp]The position loc is stored in an array loc_pha []I represents the obtained pulse ordinal index;
s405, adding 1 to t, executing S402 until t=n-L, storing the amplitude and position of the pulse signal contained in x [ n ] in the array amp [ ] and the array loc_phaj respectively, and setting the number of the pulse signals as I.
3. The MCU-based PRPD/PRPS pattern data processing method of claim 2, wherein: the specific steps of S5 are as follows:
s501, establishing an array A_seg [ ], an array B_seg [ ], an array C_seg [ ] and an array C_seg_plus [ ];
s502, extracting the position loc_phai [ I ] of the pulse signals stored in the array loc_phaj, wherein i=0 to (I-1);
s503, loc_phai]And loc a0 And loc b0 Comparison was performed:
if loc_phai]Less than loc a0 Then calculate
Figure QLYQS_10
And (loc) a0 -loc_pha[i]) And assign the difference to loc_phai]Amp [ i ]]And assigned loc_phai]Sequentially storing array A_seg [ []The method comprises the steps of carrying out a first treatment on the surface of the If loc_phai]Greater than loc a0 And is less than loc b0 Then calculateloc_pha[i]And loc a0 A difference value, which is assigned to loc_phai]Amp [ i ]]And assigned loc_phai]Sequentially storing an array B_seg [ into ] a memory]The method comprises the steps of carrying out a first treatment on the surface of the If loc_phai]Greater than loc b0 Then calculate loc_phai]And loc b0 Is assigned to loc_phai [ i ]]Amp [ i ]]And assigned loc_phai]Sequentially storing the array C_seg [ into ] a memory];
S504, combining the values in the array A_seg [ ] and the array C_seg_plus [ ] and then storing the values in the array A_seg [ ] and simultaneously clearing the array C_seg_plus [ ];
s505, storing the array C_seg [ ] into the array C_seg_plus [ ].
4. A method for processing PRPD/PRPS pattern data based on an MCU as claimed in claim 3, wherein: the specific method of S6 is as follows: uploading the positions and the amplitudes of the pulse signals stored in the arrays A_seg [ ] and B_seg [ ] to an upper computer, clearing the arrays A_seg [ ] and B_seg [ ] after uploading is finished, and returning to S2 for repeated execution.
5. An MCU-based PRPD/PRPS pattern data detection apparatus for use in an MCU-based PRPD/PRPS pattern data processing method of claim 1, characterized in that: the device comprises an MCU, wherein the MCU comprises a RAM and a serial port, the RAM is connected with an AD module through a data bus, and the MCU is connected with an upper computer through the serial port.
CN202110735971.0A 2021-06-30 2021-06-30 PRPD/PRPS map data processing method and detection device based on MCU Active CN113406453B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110735971.0A CN113406453B (en) 2021-06-30 2021-06-30 PRPD/PRPS map data processing method and detection device based on MCU

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110735971.0A CN113406453B (en) 2021-06-30 2021-06-30 PRPD/PRPS map data processing method and detection device based on MCU

Publications (2)

Publication Number Publication Date
CN113406453A CN113406453A (en) 2021-09-17
CN113406453B true CN113406453B (en) 2023-06-23

Family

ID=77680545

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110735971.0A Active CN113406453B (en) 2021-06-30 2021-06-30 PRPD/PRPS map data processing method and detection device based on MCU

Country Status (1)

Country Link
CN (1) CN113406453B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114325250B (en) * 2021-11-16 2024-02-27 国网天津市电力公司电力科学研究院 Partial discharge detection device and method integrating positioning detection and map detection functions

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4741023A (en) * 1986-12-23 1988-04-26 General Electric Company On-line test and diagnostic system for power system stabilizer
JPH02296166A (en) * 1989-05-10 1990-12-06 Nippon Telegr & Teleph Corp <Ntt> Electric-signal-waveform observing apparatus utilizing electrooptical effect
US5519300A (en) * 1993-06-29 1996-05-21 Liberty Technologies, Inc. Method and apparatus for analysis of polyphase electrical motor systems
CN1346985A (en) * 2001-09-28 2002-05-01 清华大学 Digital statistic information method for measuring ultra-long local discharge
CN202939265U (en) * 2012-10-19 2013-05-15 上海市电力公司 Vehicle-mounted transformer substation partial discharge positioning system
CN106124939A (en) * 2016-06-12 2016-11-16 西安博源电气有限公司 Distributed high tension cable partial discharge monitoring and alignment system
CN108375391A (en) * 2018-05-04 2018-08-07 苏州汇川技术有限公司 Cosine and sine signal processing method and system
CN209460370U (en) * 2018-12-26 2019-10-01 河北工业大学 A kind of squirrel cage asynchronous motor mechanical back to back test platform based on double PWM converter four quadrant runnings

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7061276B2 (en) * 2004-04-02 2006-06-13 Teradyne, Inc. Digital phase detector
CN104375067B (en) * 2014-11-18 2017-03-01 深圳供电局有限公司 A kind of loop-network switch cabinet local discharge detection device and method
CN105606975B (en) * 2016-03-09 2018-08-03 武汉华威众科电力有限公司 A kind of orientable superfrequency cable local discharge detection method and device
CN107422232A (en) * 2017-06-20 2017-12-01 国网山东省电力公司诸城市供电公司 A kind of digital live detection instrument of power distribution network terminal equipment
EP3588110B1 (en) * 2018-06-29 2022-08-17 ABB Schweiz AG Partial discharge measurement system
CN110389284B (en) * 2019-08-12 2024-04-26 平顶山学院 Detection system and detection method for equipment fault discharge on distribution network column

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4741023A (en) * 1986-12-23 1988-04-26 General Electric Company On-line test and diagnostic system for power system stabilizer
JPH02296166A (en) * 1989-05-10 1990-12-06 Nippon Telegr & Teleph Corp <Ntt> Electric-signal-waveform observing apparatus utilizing electrooptical effect
US5519300A (en) * 1993-06-29 1996-05-21 Liberty Technologies, Inc. Method and apparatus for analysis of polyphase electrical motor systems
CN1346985A (en) * 2001-09-28 2002-05-01 清华大学 Digital statistic information method for measuring ultra-long local discharge
CN202939265U (en) * 2012-10-19 2013-05-15 上海市电力公司 Vehicle-mounted transformer substation partial discharge positioning system
CN106124939A (en) * 2016-06-12 2016-11-16 西安博源电气有限公司 Distributed high tension cable partial discharge monitoring and alignment system
CN108375391A (en) * 2018-05-04 2018-08-07 苏州汇川技术有限公司 Cosine and sine signal processing method and system
CN209460370U (en) * 2018-12-26 2019-10-01 河北工业大学 A kind of squirrel cage asynchronous motor mechanical back to back test platform based on double PWM converter four quadrant runnings

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于正交分解器的介质损耗角测试仪的设计;郭松林;朱金玲;;电测与仪表(第05期);全文 *
电力电缆局部放电的高频与特高频联合检测;陈庆国等;电机与控制学报(第04期);全文 *
语音信号动态特征分析及其可视化的关键技术研究;薛丽芳;信息科技(第08期);全文 *

Also Published As

Publication number Publication date
CN113406453A (en) 2021-09-17

Similar Documents

Publication Publication Date Title
CN109387813B (en) High-precision frequency measurement method based on broadband radar signal reception
CN108256524A (en) A kind of automatic reading method of multilist index formula instrument
CN102539864B (en) Digital oscilloscope and signal measurement method
CN106646165B (en) GIS internal insulation defect classification and positioning method and system
CN113406453B (en) PRPD/PRPS map data processing method and detection device based on MCU
CN105259410B (en) The frequency estimating methods and its device of a kind of lack sampling waveform under very noisy interference
CN113608023A (en) Sampling frequency self-adaptive harmonic detection method and system
CN116520095B (en) Fault location method, system and computer readable storage medium
CN102175918B (en) Universal digital geophone implementation method and device
CN104132884B (en) The immediate processing method of a kind of signal base line in signal processing system and device
CN106777505A (en) The frequency estimating methods and device of the robust of the undersampled signal based on frequency deviation identification
CN111487462A (en) Ultra-fast frequency measurement method
CN111611686A (en) Detection method for communication signal time-frequency domain
CN112116917B (en) Phase jump degree-based method for separating acoustic signals of reactor body and fan
CN111596228B (en) Electric leakage detection method and device
CN115905782A (en) Signal processing system and method of real-time fast Fourier transform spectrometer based on FPGA
CN110736879B (en) Short-wave radio frequency direct-acquisition bridge type vector impedance detection method
CN106405341A (en) Method used for reducing partial discharge ultrahigh frequency signal
Li et al. A wideband reconnaissance receiver design based on real-time spectrum analysis technology
CN206349984U (en) A kind of broadband signal sampling playback system that compensation is divided based on frequency range
CN104320094B (en) Pulse signal standing-wave protecting circuit
CN203012130U (en) Superspeed sampling rate acquisition apparatus
CN109212576B (en) Multi-mode rapid nuclide detection system
CN206116013U (en) A pronunciation filter circuit and pronunciation drawing system for feature extraction
CN105974283B (en) Cable partial discharge feature extraction method based on wavelet packet survival index singular entropy

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