CN102830333A - Transformer substation local discharge positioning method based on electromagnetic antenna array signal processing - Google Patents

Transformer substation local discharge positioning method based on electromagnetic antenna array signal processing Download PDF

Info

Publication number
CN102830333A
CN102830333A CN2012103328979A CN201210332897A CN102830333A CN 102830333 A CN102830333 A CN 102830333A CN 2012103328979 A CN2012103328979 A CN 2012103328979A CN 201210332897 A CN201210332897 A CN 201210332897A CN 102830333 A CN102830333 A CN 102830333A
Authority
CN
China
Prior art keywords
msub
mrow
antenna
msubsup
mtr
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
CN2012103328979A
Other languages
Chinese (zh)
Other versions
CN102830333B (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.)
Shanghai Michael Energy Technology Co Ltd
Original Assignee
Shanghai Jiaotong 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 Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN201210332897.9A priority Critical patent/CN102830333B/en
Publication of CN102830333A publication Critical patent/CN102830333A/en
Application granted granted Critical
Publication of CN102830333B publication Critical patent/CN102830333B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

The invention discloses a transformer substation local discharge positioning method based on electromagnetic antenna array signal processing. The transformer substation local discharge positioning method is finished by a local discharge detection device, and the local discharge detection device comprises an omnidirectional antenna array, a built-in amplifier, a super-high-speed data sampling unit and a data processing unit which are mounted in a transformer substation. The positioning method comprises the steps as follows: constructing a matrix according to a concept of working out a signal wave arrival direction by a rotation invariant technology based on L-shaped antenna array signal processing; obtaining an azimuth angle of a local discharge source arriving at an antenna; and finishing planar location of the local discharge source of the transformer substation. A time delay sequence of a signal does not need to be calculated, so that the requirement on a sampling rate of an acquisition system can be lowered; and a plane coordinate of the local discharge source is obtained by working out a linear cross point in two wave arrival directions, namely working out a linear equation set in two unknowns, so that calculation of a nonlinear equation set is avoided.

Description

Transformer substation partial discharge positioning method based on electromagnetic wave antenna array signal processing
Technical Field
The invention belongs to the technical field of insulation monitoring of high-voltage power equipment, and relates to a transformer substation partial discharge positioning method based on electromagnetic wave antenna array signal processing.
Background
Partial discharge is both a sign and an expression form of insulation degradation and a reason for further insulation degradation, so partial discharge monitoring becomes an effective means for timely finding insulation defects of electric power equipment and avoiding insulation breakdown faults.
The positioning of the partial discharge position is beneficial to formulating a targeted maintenance treatment scheme, the power failure time is reduced, and the maintenance efficiency is improved. The ultrahigh frequency electromagnetic wave has the advantages of strong anti-interference performance, high sensitivity, stable propagation speed and the like, and is a current hotspot problem when being applied to partial discharge positioning and fault diagnosis. A group of ultra wide band omnidirectional antenna arrays are arranged in the space of a whole substation of a transformer substation to receive electromagnetic wave signals radiated by partial discharge, the position of a partial discharge source is obtained by analyzing and processing the received electromagnetic wave signals, and preliminary research and verification are carried out on the feasibility of partial discharge positioning of the space of the transformer substation by the university of Strathclyde in the UK.
The propagation speed of electromagnetic waves in the air is stable and approximate to the light speed, so that signals received by a sensor array can be analyzed and processed to obtain a signal time delay sequence, and the accurate position of a partial discharge source is obtained by solving a positioning equation system based on the time delay sequence. In practical application, under the condition of low signal-to-noise ratio, the signal delay is difficult to calculate, and a large positioning error is brought, even the positioning cannot be carried out. The equation set to be solved by the positioning method based on the signal time delay sequence is nonlinear, the computer is complex to solve and large in calculation amount, and the measurement error and the time delay calculation error can cause no solution or uncertain solution of the equation set.
The objective of the antenna array side-to-side problem is to locate the signal source by determining the direction of arrival (DOA) of the signals of the two antenna arrays and using the intersection of straight lines in the two directions. The method proposed by Roy et al for obtaining the direction of arrival of a signal by using an invariant rotation technique (ESPRIT) estimates the signal parameters by using the invariant rotation technique, and obtains a generalized eigenvalue on a unit circle by decomposing the matrix beam to determine the direction of arrival of each signal. The method has the advantage of small calculation amount, and can effectively separate independent signal sources and noise sources. The array structure restricts the performance of the estimation of the direction of arrival, and Sarkar T K and the like research the lower bound of DOA estimation variance and prove that the performance of the L-shaped array is optimal.
Disclosure of Invention
The invention aims to overcome the defects of the prior art, and provides a transformer substation partial discharge positioning method based on the L-shaped antenna array signal processing and the idea of solving the signal arrival direction by using the rotation invariant technology.
The technical solution of the invention is as follows:
the invention is composed of an omnidirectional antenna, a built-in amplifier, a high-speed data sampling unit and a data processing unit (computer) which are installed in a transformer substation. The omnidirectional antenna is used for receiving electromagnetic waves excited by partial discharge, the bandwidth of the omnidirectional antenna used for the system is 300-2000M, the received signals are amplified by the preamplifier with the same bandwidth, and the signals are transmitted to the ultra-high speed data acquisition unit through the RF coaxial shielded cable; the ultra-high speed data acquisition unit synchronously acquires 4 paths of output signals, and the sampling frequency of each channel is more than 2 GS/s; the data processing unit processes the acquired 4 paths of digital signals by using the array signals to obtain the direction of arrival angle, and the accurate position of the partial discharge source is obtained by calculation.
2M-1 antennas are respectively distributed on the x axis and the y axis of a space three-dimensional coordinate system at equal intervals to form an L-shaped array, and K (K < M) statistically independent signals are assumed in the space and are respectively directed from two-dimensional angles
Figure BDA00002117217500021
Reaches the antenna array SkIs a signal source of S'kIs SkProjection onto the xoy plane, as shown in fig. 2.
Taking the first array element of the x-axis array (i.e. the array element at the origin position) as a reference point, the sampling length of the signal is N, and the output signals of the x-axis array and the y-axis array are respectively expressed as:
X=AxS+Nx (1)
Y=AyS+Ny (2)
wherein X, Y ∈ CM×NObserved signals of two arrays respectively, x is signal, Nx,Ny∈CM×NIs the observed noise of the array. A. thex,Ay∈CM×KThe matrix is a Vandermonde matrix which is composed of array manifold vectors of K signals, and the matrix is expressed by a formula (3) and a formula (4).
Figure BDA00002117217500031
Wherein d isx,dyThe array element spacing of the x-axis array and the y-axis array is respectively, and lambda is the wavelength of the signal. The observed signals X and Y are blocked as follows:
<math> <mrow> <mi>X</mi> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>X</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mover> <mi>x</mi> <mo>&RightArrow;</mo> </mover> <mi>M</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mover> <mi>x</mi> <mo>&RightArrow;</mo> </mover> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>X</mi> <mn>2</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow> </math>
<math> <mrow> <mi>Y</mi> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>Y</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mover> <mi>y</mi> <mo>&RightArrow;</mo> </mover> <mi>M</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mover> <mi>y</mi> <mo>&RightArrow;</mo> </mover> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Y</mi> <mn>2</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein
Figure BDA00002117217500034
And
Figure BDA00002117217500035
respectively row 1 and row M signals of signal X,
Figure BDA00002117217500036
and
Figure BDA00002117217500037
respectively, the 1 st and mth row signals of signal Y. Matrix AxAnd AyThe corresponding blocks are as follows:
<math> <mrow> <msub> <mi>A</mi> <mi>x</mi> </msub> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>A</mi> <mrow> <mi>x</mi> <mn>1</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mover> <mi>a</mi> <mo>&RightArrow;</mo> </mover> <mi>xM</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mover> <mi>a</mi> <mo>&RightArrow;</mo> </mover> <mrow> <mi>x</mi> <mn>1</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>A</mi> <mrow> <mi>x</mi> <mn>2</mn> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow> </math>
<math> <mrow> <msub> <mi>A</mi> <mi>y</mi> </msub> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>A</mi> <mrow> <mi>y</mi> <mn>1</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mover> <mi>a</mi> <mo>&RightArrow;</mo> </mover> <mi>yM</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mover> <mi>a</mi> <mo>&RightArrow;</mo> </mover> <mrow> <mi>y</mi> <mn>1</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>A</mi> <mrow> <mi>y</mi> <mn>2</mn> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow> </math>
the following cross-covariance matrix is calculated:
C 1 = E [ X 1 Y 1 H ] = A x 1 R s A y 1 H + N 1 - - - ( 9 )
<math> <mrow> <msub> <mi>C</mi> <mn>2</mn> </msub> <mo>=</mo> <mi>E</mi> <mo>[</mo> <msub> <mi>X</mi> <mn>2</mn> </msub> <msubsup> <mi>Y</mi> <mn>1</mn> <mi>H</mi> </msubsup> <mo>]</mo> <mo>=</mo> <msub> <mi>A</mi> <mrow> <mi>x</mi> <mn>1</mn> </mrow> </msub> <msub> <mi>&Phi;</mi> <mi>x</mi> </msub> <msub> <mi>R</mi> <mi>s</mi> </msub> <msubsup> <mi>A</mi> <mrow> <mi>y</mi> <mn>1</mn> </mrow> <mi>H</mi> </msubsup> <mo>+</mo> <msub> <mi>N</mi> <mn>2</mn> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow> </math>
<math> <mrow> <msub> <mi>C</mi> <mn>3</mn> </msub> <mo>=</mo> <mi>E</mi> <mo>[</mo> <msub> <mi>X</mi> <mn>1</mn> </msub> <msubsup> <mi>Y</mi> <mn>2</mn> <mi>H</mi> </msubsup> <mo>]</mo> <mo>=</mo> <msub> <mi>A</mi> <mrow> <mi>x</mi> <mn>1</mn> </mrow> </msub> <msub> <mi>R</mi> <mi>s</mi> </msub> <msubsup> <mrow> <msub> <mi>&Phi;</mi> <mi>y</mi> </msub> <mi>A</mi> </mrow> <mrow> <mi>y</mi> <mn>1</mn> </mrow> <mi>H</mi> </msubsup> <mo>+</mo> <msub> <mi>N</mi> <mn>3</mn> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow> </math>
<math> <mrow> <msub> <mi>C</mi> <mn>4</mn> </msub> <mo>=</mo> <mi>E</mi> <mo>[</mo> <msub> <mi>X</mi> <mn>2</mn> </msub> <msubsup> <mi>Y</mi> <mn>2</mn> <mi>H</mi> </msubsup> <mo>]</mo> <mo>=</mo> <msub> <mi>A</mi> <mrow> <mi>x</mi> <mn>1</mn> </mrow> </msub> <msub> <mi>&Phi;</mi> <mi>x</mi> </msub> <msub> <mi>R</mi> <mi>s</mi> </msub> <mrow> <msub> <mi>&Phi;</mi> <mi>y</mi> </msub> <msubsup> <mi>A</mi> <mrow> <mi>y</mi> <mn>1</mn> </mrow> <mi>H</mi> </msubsup> </mrow> <mo>+</mo> <msub> <mi>N</mi> <mn>4</mn> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein E [. C]Representing a mathematical expectation, NiI is 1,2,3,4 is the cross covariance matrix of the noise, Rs=E[SSH]∈CK×KIs a cross-covariance matrix of signals, R, formed by assuming that K signals are independent of each othersAs a diagonal matrix, phixAnd phiyFor the rotation matrix:
Figure BDA000021172175000314
Figure BDA000021172175000315
the construction matrix C ═ C1,C2,C3,C4)TAnd calculating RC=CCH. For matrix RCPerforming eigenvalue decomposition to obtain a signal subspace expressed by a matrix, wherein the maximum K eigenvalues are { f }kI K |, 1,2, …, K } to correspond to the feature value fkThe feature vector of (a) is a matrix composed of column vectors, and the block representation is:
<math> <mrow> <mi>E</mi> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>E</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>E</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>E</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>E</mi> <mn>4</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>A</mi> <mrow> <mi>x</mi> <mn>1</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>A</mi> <mrow> <mi>x</mi> <mn>1</mn> </mrow> </msub> <msub> <mi>&Phi;</mi> <mi>x</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>A</mi> <mrow> <mi>x</mi> <mn>1</mn> </mrow> </msub> <msubsup> <mi>&Phi;</mi> <mi>y</mi> <mi>H</mi> </msubsup> </mtd> </mtr> <mtr> <mtd> <msub> <mi>A</mi> <mrow> <mi>x</mi> <mn>1</mn> </mrow> </msub> <msub> <mi>&Phi;</mi> <mi>x</mi> </msub> <msubsup> <mi>&Phi;</mi> <mi>y</mi> <mi>H</mi> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow> </math>
since the matrix T is reversible, ΩyAnd
Figure BDA00002117217500042
having the same characteristic value, ΩxAnd phixHas the same characteristic value and is omega obtained by the formula (15)yAnd ΩxThe least squares solution of (c) is:
Figure BDA00002117217500043
upper label
Figure BDA00002117217500044
Represents the Moore-Penrose Pseudo inverse of the matrix.
Figure BDA00002117217500045
The K eigenvalues of (a) are respectively:
Figure BDA00002117217500046
Figure BDA00002117217500047
the K eigenvalues of (a) are respectively:k=1,2,…,K。
note the book
Figure BDA00002117217500049
And
Figure BDA000021172175000410
from the above analysis, ukAnd vkCan be respectively composed of matrix
Figure BDA000021172175000411
And
Figure BDA000021172175000412
the eigenvalues of (2) are solved. So that a signal S can be obtainedkThe two-dimensional direction of arrival angle of (d) is:
Figure BDA000021172175000413
<math> <mrow> <msub> <mi>&theta;</mi> <mi>k</mi> </msub> <mo>=</mo> <mi>arcsin</mi> <msqrt> <msubsup> <mi>u</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>v</mi> <mi>k</mi> <mn>2</mn> </msubsup> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>18</mn> <mo>)</mo> </mrow> </mrow> </math>
compared with the prior art, the method does not need to calculate the time delay sequence of the signal, can reduce the requirement on the sampling rate of the acquisition system, can obtain the plane coordinate of the local discharge source by solving the intersection point of the straight lines in the two directions of arrival, namely solving the equation set of the first order of two elements, and avoids solving the nonlinear equation set.
Drawings
Fig. 1 is a schematic structural diagram of a radio frequency antenna array partial discharge detection apparatus according to the present invention.
Fig. 2 is a schematic diagram of an L-shaped antenna array.
Fig. 3 is a flowchart of the partial discharge detection apparatus according to the present invention.
Fig. 4 is a schematic diagram of a planar antenna array and a signal azimuth in an embodiment of the invention.
Fig. 5 shows a waveform and a frequency spectrum of a pulse current for simulation.
Fig. 6 is a schematic diagram (unit: m) of a simulated antenna array and a partially discharged source plane.
Fig. 7 is a simulated waveform of an electromagnetic wave signal received by an antenna array.
Fig. 8 is a field test antenna array layout, wherein (a) is a field test picture and (b) is an antenna coordinate system diagram.
Fig. 9 is a partial discharge signal waveform of a substation site.
Detailed Description
An embodiment of the present invention is given below with reference to the accompanying drawings. The present embodiment is implemented on the premise of the technical solution of the present invention, and a detailed implementation manner and a process are given, but the scope of the present invention is not limited to the following embodiments.
Examples
A planar antenna array composed of 4 antennas is placed in a substation, assuming that the projection of a local discharge source in the plane of the antenna array is P', firstly, an L-shaped array composed of antennas 1,2 and 4 is considered, a three-dimensional rectangular coordinate system is established by taking the antenna 1 as the origin of coordinates, and the xoy plane of the three-dimensional rectangular coordinate system is shown in FIG. 4.
This is the case when M is 2 and K is 1 in fig. 1, and the observed signal is a real signal. The following gives the local discharge source azimuth angle by calculation
Figure BDA00002117217500051
Further obtaining the position of the partial discharge source:
step 1: c is calculated from the signals observed by the antennas 1,2, 4 according to the equations (9) to (12)1、C2、C3And C4And constructing a matrix: c ═ C1,C2,C3,C4)T
Step 2: since K is 1, the matrix R is calculatedC=CCHThe eigenvector corresponding to the largest eigenvalue of (a) is the matrix E shown in equation (15).
And step 3: according to equation (16)
Figure BDA00002117217500052
And
Figure BDA00002117217500053
calculated to obtain
Figure BDA00002117217500054
And
Figure BDA00002117217500055
is a real number, i.e.:
Figure BDA00002117217500056
Figure BDA00002117217500057
u and v are obtained from formulas (19) and (20):
Figure BDA00002117217500058
Figure BDA00002117217500059
from equations (21, 22), and the two-dimensional direction angle of arrival of the signal
Figure BDA000021172175000510
The calculation formula (17, 18) of theta is known as the azimuth angle
Figure BDA000021172175000511
Independent of the wavelength of the signal, and theta is dependent on wavelength, discussed below by calculation
Figure BDA000021172175000512
The method for positioning the partial discharge plane of the transformer substation is realized.
And 4, step 4: calculation of azimuth from equation (17)
Figure BDA000021172175000513
By using the same method, considering the L-shaped array composed of the antennas 3, 2, 4, and taking the antenna 3 as the origin of coordinates, a rectangular coordinate system is established according to the dashed line in fig. 2, so as to obtain the azimuth angle of the local discharge source in the coordinate system x ' o ' y ', and the complementary angle is the azimuth angle in fig. 4
Figure BDA00002117217500061
Let the coordinates of the antennas 1,2,3,4 in the coordinate system xoy be (0,0), (d) respectivelyx,0)、(dx,dy)、(0,dy) The projection P' of the partial discharge source in the xoy plane has coordinates (x, y).
By azimuth angle
Figure BDA00002117217500062
And
Figure BDA00002117217500063
and the coordinates of antenna 1 and antenna 3 may yield two linear equations:
Figure BDA00002117217500064
and solving a linear equation system of two variables consisting of the formula (24) and the formula (25) to obtain the coordinates (x, y) of the projection P' of the local discharge source in the xoy plane, namely realizing the plane positioning of the local discharge source.
Simulation and test verification of method
In order to verify the accuracy and feasibility of the ultrahigh frequency array signal processing-based substation partial discharge positioning method, the partial discharge signals obtained by electromagnetic wave simulation software and the ultrahigh frequency electromagnetic wave signals actually measured by the substation are analyzed and processed respectively, and the position coordinates of the partial discharge source are calculated by using the method.
Simulation verification
Ansoft HFSS 13.0version software is adopted to perform simulation of partial discharge radiation electromagnetic waves. The simulation partial discharge source adopts Gaussian pulse with the amplitude of pulse current of 1A and the pulse width of 4ns, the current direction is the positive direction of the z axis, and the waveform diagram and the frequency spectrum analysis of the pulse current are shown in figure 5.
Establishing a three-dimensional coordinate system in a simulation substation space, wherein the position coordinate of a simulation partial discharge source is P (1.5,6,2), the coordinates of a radio frequency antenna array are respectively (0,0,0), (3,4,0) and (0,4,0), and the unit is as follows: m, assuming that the electromagnetic wave signal radiated by the partial discharge can reach the receiving antenna without a shielding straight line. A schematic plan view of the simulated antenna array and the partial discharge source is shown in fig. 6.
The results of the electromagnetic wave simulation software obtaining the waveforms of the electromagnetic wave signals received by the 4 antennas are shown in fig. 7. By using the method, the azimuth angle of the partial discharge source in fig. 4 can be obtained by analyzing and processing the simulated signal
Figure BDA00002117217500066
And
Figure BDA00002117217500067
75.3 deg. and 126.9 deg., respectively. Thereby obtaining a system of linear equations shown in formulas (26) and (27), and solving to obtain the plane position of the partial discharge source as (1.55, 5.93), unit: m, compared with the real plane position (1.5,6), the error is 8.6cm, and the requirement of positioning precision is met.
Verification of field test result of transformer substation
In order to verify the effectiveness of the algorithm under the condition of on-site interference of a transformer substation, an ultrahigh frequency signal actually measured by an antenna array partial discharge detection positioning system at a certain 500kV transformer substation is explained. The four ultrahigh frequency omnidirectional antennas of the system form a planar rectangular array, the length and the width of the array are respectively 1.58m and 1.16m, and the layout of the antenna array is shown in fig. 8. A coordinate system is established by taking the antenna 1 as a coordinate origin, and coordinates of the four antennas in the xoy plane are respectively (0,0), (1.58,1.16) and (0,1.16), and the unit is: and m is selected.
A simulated partial discharge source is arranged in the transformer substation, and the coordinate position is (-4,5, 5) m. The high-speed oscilloscope with the storage function is used as a signal acquisition system to synchronously acquire and store ultrahigh frequency electromagnetic waves radiated by the same discharge source and received by an antenna array, fig. 9 shows signal waveforms received by a group of 4 sensors synchronously acquired by the oscilloscope, the sampling rate is 5G, and the signal length is 2 mus (10000 sampling points).
By using the array signal processing method, the local discharge ultrahigh frequency signal acquired on site is analyzed and processed, and the azimuth angle of the partial discharge source in a coordinate system established by taking the antenna 1 as the origin and the antenna 3 as the origin can be obtained
Figure BDA00002117217500071
And
Figure BDA00002117217500072
126.2 deg. and 143.8 deg., respectively. Thereby obtaining a system of linear equations shown in formulas (26) and (27), and solving to obtain the plane position of the partial discharge source as (-3.83, 5.23), unit: and m, compared with the real plane position (-4,5), the error is 28.6cm, and the requirement of the substation partial discharge positioning precision is met. An antenna coordinate system (unit: m) is applied.

Claims (1)

1. A transformer substation partial discharge positioning method based on electromagnetic wave antenna array signal processing is completed by a partial discharge detection device, wherein the device is composed of an omnidirectional antenna array, a built-in amplifier, an ultra-high speed data sampling unit and a data processing unit which are arranged in a transformer substation; the omnidirectional antenna is used for receiving electromagnetic waves excited by partial discharge, the bandwidth of the omnidirectional antenna array is 300-2000M, and the received signals are amplified by the preamplifier with the same bandwidth and then transmitted to the ultra-high speed data acquisition unit; the method is characterized in that the positioning method comprises the following specific steps:
establishing a three-dimensional direct coordinate system xoy by taking a first antenna as an origin, a first antenna and a second antenna as an x axis and a first antenna and a third antenna as a y axis, wherein the first antenna, the second antenna and the third antenna form a first L-shaped antenna array;
calculating a cross covariance matrix according to signals received by the L-type antenna array:
C 1 = E [ X 1 Y 1 H ] = A x 1 R s A y 1 H + N 1 - - - ( 9 )
<math> <mrow> <msub> <mi>C</mi> <mn>2</mn> </msub> <mo>=</mo> <mi>E</mi> <mo>[</mo> <msub> <mi>X</mi> <mn>2</mn> </msub> <msubsup> <mi>Y</mi> <mn>1</mn> <mi>H</mi> </msubsup> <mo>]</mo> <mo>=</mo> <msub> <mi>A</mi> <mrow> <mi>x</mi> <mn>1</mn> </mrow> </msub> <msub> <mi>&Phi;</mi> <mi>x</mi> </msub> <msub> <mi>R</mi> <mi>s</mi> </msub> <msubsup> <mi>A</mi> <mrow> <mi>y</mi> <mn>1</mn> </mrow> <mi>H</mi> </msubsup> <mo>+</mo> <msub> <mi>N</mi> <mn>2</mn> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow> </math>
<math> <mrow> <msub> <mi>C</mi> <mn>3</mn> </msub> <mo>=</mo> <mi>E</mi> <mo>[</mo> <msub> <mi>X</mi> <mn>1</mn> </msub> <msubsup> <mi>Y</mi> <mn>2</mn> <mi>H</mi> </msubsup> <mo>]</mo> <mo>=</mo> <msub> <mi>A</mi> <mrow> <mi>x</mi> <mn>1</mn> </mrow> </msub> <msub> <mi>R</mi> <mi>s</mi> </msub> <msubsup> <mrow> <msub> <mi>&Phi;</mi> <mi>y</mi> </msub> <mi>A</mi> </mrow> <mrow> <mi>y</mi> <mn>1</mn> </mrow> <mi>H</mi> </msubsup> <mo>+</mo> <msub> <mi>N</mi> <mn>3</mn> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow> </math>
<math> <mrow> <msub> <mi>C</mi> <mn>4</mn> </msub> <mo>=</mo> <mi>E</mi> <mo>[</mo> <msub> <mi>X</mi> <mn>2</mn> </msub> <msubsup> <mi>Y</mi> <mn>2</mn> <mi>H</mi> </msubsup> <mo>]</mo> <mo>=</mo> <msub> <mi>A</mi> <mrow> <mi>x</mi> <mn>1</mn> </mrow> </msub> <msub> <mi>&Phi;</mi> <mi>x</mi> </msub> <msub> <mi>R</mi> <mi>s</mi> </msub> <mrow> <msub> <mi>&Phi;</mi> <mi>y</mi> </msub> <msubsup> <mi>A</mi> <mrow> <mi>y</mi> <mn>1</mn> </mrow> <mi>H</mi> </msubsup> </mrow> <mo>+</mo> <msub> <mi>N</mi> <mn>4</mn> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein, E [. C]Representing a mathematical expectation, NiI is 1,2,3,4 is the cross covariance matrix of the noise, Rs=E[SSH]∈CK×KIs a cross-covariance matrix of signals, R, formed by assuming that K signals are independent of each othersAs a diagonal matrix, phixAnd phiyFor the rotation matrix:
Figure FDA00002117217400015
Figure FDA00002117217400016
the construction matrix C ═ C1,C2,C3,C4)T
Calculating RC=CCH
For matrix RCPerforming eigenvalue decomposition to obtain a signal subspace expressed by a matrix, wherein the maximum K eigenvalues are { f }kI K ═ 1,2,. K }, to correspond to the characteristic value fkThe feature vector of (a) is a matrix composed of column vectors, and the block representation is:
<math> <mrow> <mi>E</mi> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>E</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>E</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>E</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>E</mi> <mn>4</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>A</mi> <mrow> <mi>x</mi> <mn>1</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>A</mi> <mrow> <mi>x</mi> <mn>1</mn> </mrow> </msub> <msub> <mi>&Phi;</mi> <mi>x</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>A</mi> <mrow> <mi>x</mi> <mn>1</mn> </mrow> </msub> <msubsup> <mi>&Phi;</mi> <mi>y</mi> <mi>H</mi> </msubsup> </mtd> </mtr> <mtr> <mtd> <msub> <mi>A</mi> <mrow> <mi>x</mi> <mn>1</mn> </mrow> </msub> <msub> <mi>&Phi;</mi> <mi>x</mi> </msub> <msubsup> <mi>&Phi;</mi> <mi>y</mi> <mi>H</mi> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein T is a reversible matrix;
' Ji
Figure FDA00002117217400022
Figure FDA00002117217400023
Calculate ΩyAnd ΩxThe least squares solution of (a) is as follows:
Figure FDA00002117217400024
upper label
Figure FDA00002117217400025
Represents the Moore-Penrose Pseudo inverse of the matrix,
Figure FDA00002117217400026
the K eigenvalues of (a) are respectively:
Figure FDA00002117217400027
Figure FDA00002117217400028
the K eigenvalues of (a) are respectively:
Figure FDA00002117217400029
k=1,2,…,K;
calculating signal SkThe two-dimensional azimuth angle to the first antenna is:
Figure FDA000021172174000210
<math> <mrow> <msub> <mi>&theta;</mi> <mi>k</mi> </msub> <mo>=</mo> <mi>arcsin</mi> <msqrt> <msubsup> <mi>u</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>v</mi> <mi>k</mi> <mn>2</mn> </msubsup> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>18</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein,
Figure FDA000021172174000212
and
Figure FDA000021172174000213
establishing a second three-dimensional direct coordinate system x ' o ' y ' by taking the fourth antenna as an original point, the fourth antenna and the third antenna as an x axis and the fourth antenna and the second antenna as a y axis, wherein the fourth antenna, the second antenna and the third antenna form a second L-shaped antenna array;
and seventhly, calculating a two-dimensional azimuth angle of the signal Sk to the fourth antenna:
Figure FDA000021172174000214
calculating the projection coordinate of the local discharge source in the xoy coordinate plane to realize the plane location of the local discharge source, wherein the formula is as follows:
Figure FDA000021172174000215
wherein, (x, y) is the plane coordinate of the partial discharge source in the transformer substation, and dx,dyThe spacing of the antenna array x-axis and y-axis antennas respectively.
CN201210332897.9A 2012-09-10 2012-09-10 Transformer substation local discharge positioning method based on electromagnetic antenna array signal processing Active CN102830333B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210332897.9A CN102830333B (en) 2012-09-10 2012-09-10 Transformer substation local discharge positioning method based on electromagnetic antenna array signal processing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210332897.9A CN102830333B (en) 2012-09-10 2012-09-10 Transformer substation local discharge positioning method based on electromagnetic antenna array signal processing

Publications (2)

Publication Number Publication Date
CN102830333A true CN102830333A (en) 2012-12-19
CN102830333B CN102830333B (en) 2015-06-10

Family

ID=47333536

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210332897.9A Active CN102830333B (en) 2012-09-10 2012-09-10 Transformer substation local discharge positioning method based on electromagnetic antenna array signal processing

Country Status (1)

Country Link
CN (1) CN102830333B (en)

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103592584A (en) * 2013-11-19 2014-02-19 国家电网公司 Substation discharging source locating method based on electromagnetic wave detection
CN103713242A (en) * 2013-12-27 2014-04-09 上海交通大学 Novel ultrahigh frequency sensor for positioning local discharge source space and array thereof
CN103713243A (en) * 2013-12-30 2014-04-09 国家电网公司 10-V switch cabinet internal discharge detection device based on electromagnetic wave antenna receiving array
CN103913680A (en) * 2014-03-17 2014-07-09 上海交通大学 Partial discharging positioning method based on time delay error stochastic simulation and statistic analysis
CN103969559A (en) * 2014-05-15 2014-08-06 国家电网公司 Method and system for positioning local discharge space of power equipment of transformer substation
CN104133164A (en) * 2014-07-07 2014-11-05 国家电网公司 Visible projection method in space partial-discharge positioning
CN104574536A (en) * 2014-11-18 2015-04-29 深圳市金溢科技股份有限公司 RSU (Road Side Unit) equipment of ETC (Electronic Toll Collection) system
CN104614611A (en) * 2015-01-30 2015-05-13 电子科技大学 Method for detecting damaged element of receiving antenna array online
CN104614653A (en) * 2015-02-09 2015-05-13 国家电网公司 Array antenna based multiple local discharge point positioning and distinguishing method for local discharge detection device
CN105223480A (en) * 2015-10-22 2016-01-06 西安交通大学 The Positioning Error Simulation method of aerial array time difference method location transformer station Partial Discharge Sources
CN105866636A (en) * 2016-03-23 2016-08-17 国家电网公司 Transformer substation positioning method based on time difference positioning
CN108445447A (en) * 2018-02-27 2018-08-24 国家电网有限公司 A kind of station domain space Mutual coupling system of substation's discharge source
CN109375074A (en) * 2018-11-16 2019-02-22 广州供电局有限公司 Detection method for local discharge, device, equipment and system based on ultrahigh-frequency signal
CN110146842A (en) * 2019-06-14 2019-08-20 哈尔滨工业大学 Signal carrier frequency and two dimension DOA method for parameter estimation based on lack sampling
CN110297214A (en) * 2019-07-17 2019-10-01 南京航空航天大学 Mostly relatively prime array cooperates with indoor radiation source positioning device and method
CN112986768A (en) * 2021-03-22 2021-06-18 保定天威保变电气股份有限公司 Partial discharge positioning device and method
CN116819262A (en) * 2023-08-31 2023-09-29 陕西公众电气股份有限公司 Insulation defect identification and positioning method for high-voltage electrical equipment

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101571568A (en) * 2009-05-15 2009-11-04 河北省电力研究院 Ultrasonic planar array sensor suitable for partial discharge detection of transformer
CN101702000A (en) * 2009-11-26 2010-05-05 华北电力大学(保定) Phased ultrasound locating method and system of partial discharge of transformer
JP2011242244A (en) * 2010-05-18 2011-12-01 Mitsui Eng & Shipbuild Co Ltd Partial discharge locator and partial discharge location method
CN202141781U (en) * 2011-06-21 2012-02-08 山东电力研究院 Online monitoring and positioning device for transformer station partial discharge signals

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101571568A (en) * 2009-05-15 2009-11-04 河北省电力研究院 Ultrasonic planar array sensor suitable for partial discharge detection of transformer
CN101702000A (en) * 2009-11-26 2010-05-05 华北电力大学(保定) Phased ultrasound locating method and system of partial discharge of transformer
JP2011242244A (en) * 2010-05-18 2011-12-01 Mitsui Eng & Shipbuild Co Ltd Partial discharge locator and partial discharge location method
CN202141781U (en) * 2011-06-21 2012-02-08 山东电力研究院 Online monitoring and positioning device for transformer station partial discharge signals

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
侯慧娟等: "基于超高频电磁波的变电站局部放电空间定位", 《高电压技术》, vol. 38, no. 6, 30 June 2012 (2012-06-30), pages 1334 - 1340 *
俞燕乐等: "电力变压器局部放电超宽带射频定位研究", 《机电信息》, no. 9, 25 March 2012 (2012-03-25), pages 149 - 150 *
刁鸣等: "基于修正ESPRIT算法的二维DOA估计", 《哈尔滨工程大学学报》, vol. 29, no. 4, 15 April 2008 (2008-04-15), pages 407 - 410 *

Cited By (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103592584B (en) * 2013-11-19 2016-02-10 国家电网公司 A kind of transformer station's discharge source localization method based on Electromagnetic Wave Detection
CN103592584A (en) * 2013-11-19 2014-02-19 国家电网公司 Substation discharging source locating method based on electromagnetic wave detection
CN103713242A (en) * 2013-12-27 2014-04-09 上海交通大学 Novel ultrahigh frequency sensor for positioning local discharge source space and array thereof
CN103713243A (en) * 2013-12-30 2014-04-09 国家电网公司 10-V switch cabinet internal discharge detection device based on electromagnetic wave antenna receiving array
CN103713243B (en) * 2013-12-30 2016-08-24 国家电网公司 10V inside switch cabinet discharging detection device based on electromagnetic antenna receiving array
CN103913680A (en) * 2014-03-17 2014-07-09 上海交通大学 Partial discharging positioning method based on time delay error stochastic simulation and statistic analysis
CN103913680B (en) * 2014-03-17 2016-09-07 上海交通大学 Partial discharge positioning method based on time delay error stochastic simulation Yu statistical analysis
CN103969559A (en) * 2014-05-15 2014-08-06 国家电网公司 Method and system for positioning local discharge space of power equipment of transformer substation
CN104133164A (en) * 2014-07-07 2014-11-05 国家电网公司 Visible projection method in space partial-discharge positioning
CN104574536A (en) * 2014-11-18 2015-04-29 深圳市金溢科技股份有限公司 RSU (Road Side Unit) equipment of ETC (Electronic Toll Collection) system
CN104574536B (en) * 2014-11-18 2017-02-22 深圳市金溢科技股份有限公司 RSU (Road Side Unit) equipment of ETC (Electronic Toll Collection) system
CN104614611A (en) * 2015-01-30 2015-05-13 电子科技大学 Method for detecting damaged element of receiving antenna array online
CN104614611B (en) * 2015-01-30 2017-06-16 电子科技大学 A kind of online test method of the impaired array element of receiving antenna array
CN104614653A (en) * 2015-02-09 2015-05-13 国家电网公司 Array antenna based multiple local discharge point positioning and distinguishing method for local discharge detection device
CN105223480B (en) * 2015-10-22 2018-01-19 西安交通大学 The Positioning Error Simulation method of aerial array time difference method positioning transformer station Partial Discharge Sources
CN105223480A (en) * 2015-10-22 2016-01-06 西安交通大学 The Positioning Error Simulation method of aerial array time difference method location transformer station Partial Discharge Sources
CN105866636B (en) * 2016-03-23 2018-09-11 国家电网公司 Transformer substation positioning method based on time difference positioning
CN105866636A (en) * 2016-03-23 2016-08-17 国家电网公司 Transformer substation positioning method based on time difference positioning
CN108445447A (en) * 2018-02-27 2018-08-24 国家电网有限公司 A kind of station domain space Mutual coupling system of substation's discharge source
CN108445447B (en) * 2018-02-27 2021-09-28 国家电网有限公司 Substation area space direction of arrival estimation system of substation discharge source
CN109375074A (en) * 2018-11-16 2019-02-22 广州供电局有限公司 Detection method for local discharge, device, equipment and system based on ultrahigh-frequency signal
CN110146842A (en) * 2019-06-14 2019-08-20 哈尔滨工业大学 Signal carrier frequency and two dimension DOA method for parameter estimation based on lack sampling
CN110297214A (en) * 2019-07-17 2019-10-01 南京航空航天大学 Mostly relatively prime array cooperates with indoor radiation source positioning device and method
CN112986768A (en) * 2021-03-22 2021-06-18 保定天威保变电气股份有限公司 Partial discharge positioning device and method
CN112986768B (en) * 2021-03-22 2022-03-15 保定天威保变电气股份有限公司 Partial discharge positioning device and method
CN116819262A (en) * 2023-08-31 2023-09-29 陕西公众电气股份有限公司 Insulation defect identification and positioning method for high-voltage electrical equipment
CN116819262B (en) * 2023-08-31 2023-11-10 陕西公众电气股份有限公司 Insulation defect identification and positioning method for high-voltage electrical equipment

Also Published As

Publication number Publication date
CN102830333B (en) 2015-06-10

Similar Documents

Publication Publication Date Title
CN102830333B (en) Transformer substation local discharge positioning method based on electromagnetic antenna array signal processing
US10663563B2 (en) On-site calibration of array antenna systems
CN102662132B (en) A kind of local discharge of electrical equipment localization method and system
CN103399291B (en) Fast sparse recovery-based super-resolution direction-of-arrival estimation method
CN101702000A (en) Phased ultrasound locating method and system of partial discharge of transformer
CN106556783A (en) A kind of shelf depreciation direction-finding method in transformer station based on superfrequency phased array principle
CN109581078B (en) Directional diagram measuring system and method suitable for antenna in half-space environment
CN110297213A (en) Radiation source positioning device and method based on the unmanned aerial vehicle platform for loading relatively prime linear array
CN107783086B (en) Method for diagnosing distorted position of antenna array aperture amplitude phase field
CN104614653A (en) Array antenna based multiple local discharge point positioning and distinguishing method for local discharge detection device
CN107085171B (en) A kind of optimization placement method of the superfrequency phased array for substation&#39;s shelf depreciation direction finding
CN106443587A (en) High-resolution rapid deconvolution sound source imaging algorithm
WO2014106414A1 (en) Locating method of partial discharge source and locating system of partial discharge source
CN103983903A (en) Whole substation partial discharge locating method for utilizing oscilloscope to detect radio-frequency signals
CN108181557B (en) Method for determining ultrahigh frequency partial discharge signal direction
CN103913680A (en) Partial discharging positioning method based on time delay error stochastic simulation and statistic analysis
CN103235286B (en) High-precision locating method for electric noise sources
CN107515956B (en) Large finite plane array analysis method based on HFSS unit method
CN105866636B (en) Transformer substation positioning method based on time difference positioning
CN110515038A (en) It is a kind of based on the adaptive passive location device of unmanned plane-array and implementation method
CN103777186A (en) Method for calculating near-field radar echo characteristics of moving object based on beam decomposition and local irradiation
Nannuru et al. Sparse planar arrays for azimuth and elevation using experimental data
CN103048641B (en) Vector array steady focusing processing method with combined constraint optimization form
CN114563756A (en) Linear interpolation method for direction finding sample of dual-channel correlation interferometer
CN103592584A (en) Substation discharging source locating method based on electromagnetic wave detection

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20160316

Address after: 200240 Shanghai city Minhang District Jianchuan Road No. 951 Building 5 floor A1065 room 1

Patentee after: Shanghai Judian Electric Equipment Co., Ltd.

Address before: 200240 Dongchuan Road, Shanghai, No. 800, No.

Patentee before: Shanghai Jiao Tong University

TR01 Transfer of patent right

Effective date of registration: 20190118

Address after: Room 3070, Building B, 555 Dongchuan Road, Minhang District, Shanghai, 200241

Patentee after: Shanghai Michael Energy Technology Co., Ltd.

Address before: Room A1065, 1st floor, No. 951 Jianchuan Road, Minhang District, Shanghai, 200240

Patentee before: Shanghai Judian Electric Equipment Co., Ltd.

TR01 Transfer of patent right