CN106645944A - Extraction method for instantaneous value of each component of three-phase power grid based on rotation transformation - Google Patents
Extraction method for instantaneous value of each component of three-phase power grid based on rotation transformation Download PDFInfo
- Publication number
- CN106645944A CN106645944A CN201610860915.9A CN201610860915A CN106645944A CN 106645944 A CN106645944 A CN 106645944A CN 201610860915 A CN201610860915 A CN 201610860915A CN 106645944 A CN106645944 A CN 106645944A
- Authority
- CN
- China
- Prior art keywords
- line voltage
- component
- phase
- sequence
- voltage
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses an extraction method for the instantaneous value of each component of a three-phase power grid based on rotation transformation. The extraction method comprises the following steps that power grid voltage under a three-phase ABC coordinate system is transformed to an alphabeta coordinate system through static coordinate transformation; then a zero sequence component, a positive sequence fundamental component, a negative sequence fundamental component and various harmonic components under the alphabeta coordinate system are acquired through a rotation transformation matrix; and finally the instantaneous positive sequence angle and the instantaneous negative sequence angle of the power grid are calculated by the acquired positive and negative sequence components. According to the method, various components of the power grid voltage can be extracted by the voltage sampling data of multiple historical periods with no requirement for the data of a complete power grid period so that the problems that the existing method can acquire the three-phase power grid instantaneous phase while the calculation method is complex and the delay is long can be solved. The method is easy to implement, low in delay and high in accuracy.
Description
Technical field
The present invention relates to electric field, more particularly to a kind of extracting method of the various components of line voltage, it is adaptable to electrical network
Asymmetrical voltage and during containing harmonic component, while obtain including zero-sequence component instantaneous value, positive sequence fundametal compoment instantaneous value, negative phase-sequence
Fundametal compoment instantaneous value, harmonic component instantaneous value.
Background technology
In application scenarios such as generations of electricity by new energy, it is desirable to which control system can quickly and accurately detect grid voltage amplitude
It is that control system feeds back in time accurate electric network information with instantaneous phase.At present, the phaselocked loop being widely used in control system
Mainly hardware phase lock high precision (PLL) and software phlase locking (SPLL) two ways.When line voltage is asymmetric and containing harmonic component,
This requirement can not be met well.Therefore it is required to compatible line voltage asymmetric and containing electricity during harmonic component
Net voltage phase-locked loop method, makes that control system tackles in time overcurrent that three-phase voltage asymmetry causes and grid-connected converter is straight
The fluctuation of stream side voltage.
Electrical network it is asymmetric and containing harmonic wave under, a key technical problem of Design of PLL is the various components of electrical network
Extract.At present, the solution that existing pertinent literature is proposed, mainly includes：Based on the decoupling method of double synchronous coordinate systems, prolong
When target signal filter method, moving average filter method etc..Decoupling method based on double synchronous coordinate systems is sat using the double synchronous of positive-negative sequence
Mark system, aligns, negative phase-sequence fundametal compoment realizes decoupling, and this method complex structure, the response time is longer and computationally intensive, realization
Get up relatively difficult.It is in synchronous rotating frame, by the voltage signal U containing nth harmonic that time delayed signal filters methoddTime delay
T/2 (n-1), then with UdIt is added, if containing various harmonic waves, needs are computed repeatedly by filtering the multiple of different harmonic waves.But these
Method all cannot be obtained including zero-sequence component instantaneous value, positive sequence fundametal compoment instantaneous value, negative phase-sequence fundametal compoment instantaneous value, humorous simultaneously
Wave component instantaneous value.
The content of the invention
The technical problem to be solved is the defect for overcoming prior art to exist, there is provided one kind is based on rotation
The extracting method of each component of three phase network of conversion, it is asymmetric and during containing harmonic component to solve line voltage, while obtaining
Ask including zero-sequence component instantaneous value, positive sequence fundametal compoment instantaneous value, negative phase-sequence fundametal compoment instantaneous value, harmonic component instantaneous value
Topic.
A kind of extracting method of each component of three phase network based on rotation transformation proposed by the present invention, comprises the steps：
Step 1, converted according to CLARKE, the line voltage under three-phase ABC coordinate systems is transformed to into α β coordinate systems, obtain α β
Line voltage instantaneous value under coordinate system, its expression formula is various component of voltage sums,
Write as plural form, such as following formula：
In formula：Emk--- the amplitude of line voltage k component of degree n ns；--- the starting phase angle of line voltage k component of degree n ns；
ω --- line voltage fundamental wave angular frequency；T --- the line voltage under three-phase ABC coordinate systems transforms to the time of α β coordinate systems
Point；E --- the line voltage vector under three-phase ABC coordinate systems；eɑ--- the voltage after line voltage α β coordinate transforms on α axles
Value；eβ--- the magnitude of voltage after line voltage α β coordinate transforms on β axles；J --- the imaginary part of symbol；K=0, ± 1, ± 2 ..., -50≤
k≤50。
Step 2, by the line voltage instantaneous value discretization obtained by coordinate transform, such as following formula：
In formula：TS--- the line voltage cycle；Em0--- the zero-sequence component amplitude of line voltage；--- line voltage
The phase angle of zero-sequence component；E-mk--- the peak value of line voltage-k component of degree n ns；--- line voltage-k component of degree n ns it is initial
Phase angle；
N is constant；The line voltage vector of e (n) --- current period；E (n-1) --- the current period previous cycle
Line voltage vector；E (n+1) --- the line voltage vector of a cycle after current period.
Step 3, complex column vector Ε (n) that the line voltage in front l cycle is expressed as an exponential form,
The instantaneous value of the various components of current period line voltage, is expressed as the complex column vector ε of an exponential form
(n), such as following formula：
In formula：----k time, k component of degree n ns of current period line voltage, k=1,2 ..., k≤50；
e-m1(n),em1The negative phase-sequence of (n) --- current period line voltage, positive sequence fundametal compoment；em0(n) --- current period electrical network electricity
The zero-sequence component of pressure.
Step 4, according to Ε (n) and the expression formula of ε (n), obtain a spin matrix R, as l=2k+1, further meter
Calculate the various components of current electric grid voltage, such as following formula：
ε (n)=R-1Ε(n) (7)
Step 5：By by the positive and negative sequence fundametal compoment of voltage for calculating acquisition, further calculate and obtain the positive and negative sequence phase angle of voltage
Sine and cosine value, such as following formula：
In formula：θ+,θ---- the instantaneous phase angle of positive and negative sequence fundametal compoment；
Re(em1(n)),Im(em1(n)) --- the real part and imaginary part of positive sequence fundametal compoment；
Re(e-m1(n)),Im(e-m1(n)) --- the real part and imaginary part of negative phase-sequence fundametal compoment.
Further, the instantaneous positive sequence angle and instantaneous negative phase-sequence angle of three phase network are obtained.
The invention has the beneficial effects as follows：By rotational transformation matrix, the various components of three-phase power grid voltage are extracted, by obtaining
Positive and negative sequence fundametal compoment calculate the instantaneous positive sequence angle and instantaneous negative phase-sequence angle of electrical network.This method only needs to multiple going through
The voltage sample data in history cycle can extract the various components of line voltage, and without the need for the number of a complete grid cycle
According to.Existing method Computational Methods complexity and the longer problem of time delay are avoided, it is an advantage of the current invention that：Realize it is simple, prolong
When little, high precision
Description of the drawings
Fig. 1:Based on the extracting method that the present invention is provided, voltage negative sequence component polar plot is obtained to obtain；
Fig. 2:Based on the extracting method that the present invention is provided, voltage positive-sequence component polar plot is obtained to obtain；
Fig. 3:Based on the extracting method that the present invention is provided, the order harmonic components polar plot of voltage 5 is obtained to obtain；
Fig. 4:Based on the extracting method that the present invention is provided, the negative sequence component of acquisition, positive-sequence component, the synthesis of 5 order harmonic components
Current period line voltage polar plot；
Fig. 5:Based on the extracting method that the present invention is provided, to electrical network in the case of 30 degree of three-phase power grid voltage phase place impact
The experimental result that voltage positive and negative sequence component and 5 order harmonic components are extracted.
Fig. 6:It is right in the case of three-phase power grid voltage BC two-phase voltage falls 50% based on the extracting method that the present invention is provided
The experimental result that line voltage positive and negative sequence component and 5 order harmonic components are extracted.
Fig. 7:Based on the extracting method that the present invention is provided, in three-phase power grid voltage 5 subharmonic feelings of the amplitude for 100V are injected
Line voltage positive and negative sequence component and 5 order harmonic components are extracted experimental result under condition.
Fig. 8:Based on the extracting method that the present invention is provided, fall in 30 degree of three-phase power grid voltage phase place impact, BC two-phase voltages
Fall 50% and injection amplitude for 100V 5 subharmonic, three kinds of situations it is simultaneous in the case of to line voltage positive and negative sequence component
And 5 order harmonic components extract experimental result.
Specific embodiment
With reference to embodiment and accompanying drawing, the present invention is described in further detail.
Principle of the invention and feature for convenience of description, to extract positive and negative sequence fundametal compoment and 5 order harmonic components as a example by,
The inventive method is illustrated, example only for convenience of being used to explain the present invention, is not intended to limit the scope of the present invention.
Embodiment：A kind of extraction side of each component instantaneous value of three phase network based on rotation transformation is proposed using the present invention
Method, extracts the positive and negative sequence fundametal compoment of three phase network and 5 order harmonic components.Its specific implementation process is as follows：
Step 1：Converted by static coordinate, the line voltage under three-phase ABC coordinate systems is transformed to into α β coordinate systems, obtained
Line voltage instantaneous value under α β coordinate systems, its expression formula is positive and negative sequence fundametal compoment and 5 order harmonic components sums, and is write as
Plural form, such as following formula：
In formula：E-m1,Em1,Em5--- the peak value of the positive and negative sequence fundametal compoment of line voltage and 5 order harmonic components；
--- the starting phase angle of the positive and negative sequence fundametal compoment of line voltage and 5 order harmonic components；
ω --- line voltage fundamental wave angular frequency.
Step 2：By the line voltage instantaneous value discretization obtained by coordinate transform, such as following formula：
In formula：TS--- the line voltage cycle；
Step 3：The line voltage in front 3 cycles is expressed as into column vector Ε (n), current period line voltage
Positive and negative sequence fundametal compoment and 5 order harmonic components are expressed as column vector ε (n), such as following formula：
In formula：e-m1(n),em1(n),em5The negative phase-sequence of (n) --- current period line voltage, positive sequence fundametal compoment and 5 times
Harmonic component.
Step 4：According to Ε (n) and the expression formula of ε (n), a 3 rank spin matrix R can be obtained, further calculate and work as
The positive and negative sequence fundametal compoment of front line voltage and 5 order harmonic components, such as following formula：
In formula：λ-m1,λ′-m1,λ″-m1--- the n-th cycle, the (n-1)th cycle, the n-th -2 cycle electrical network electricity when extracting negative sequence component
Amplitude gain of the pressure after rotation；
λm1,λ′m1,λ″m1--- the n-th cycle, the (n-1)th cycle, the n-th -2 cycle line voltage are passed through when extracting positive-sequence component
Amplitude gain after rotation；
λm5,λ′m5,λ″m5--- the n-th cycle, the (n-1)th cycle, the n-th -2 cycle line voltage Jing when extracting 5 order harmonic components
The amplitude gain crossed after rotation；
θ-m1,θ′-m1,θ″-m1--- the n-th cycle, the (n-1)th cycle, the n-th -2 cycle line voltage Jing when extracting negative sequence component
The phase shift crossed after rotation；
θm1,θ′m1,θ″m1--- the n-th cycle, the (n-1)th cycle, the n-th -2 cycle line voltage are passed through when extracting positive-sequence component
Phase shift after rotation；
θm5,θ′m5,θ″m5--- the n-th cycle, the (n-1)th cycle, the n-th -2 cycle line voltage Jing when extracting 5 order harmonic components
The phase shift crossed after rotation.
Step 5：By by the positive and negative sequence fundametal compoment of voltage for calculating acquisition, further calculate and obtain the positive and negative sequence phase angle of voltage
Sine and cosine value, such as following formula：
In formula：θ+,θ---- the instantaneous phase angle of positive and negative sequence fundametal compoment；
Re(em1(n)),Im(em1(n)) --- the real part and imaginary part of positive sequence fundametal compoment；
Re(e-m1(n)),Im(e-m1(n)) --- the real part and imaginary part of negative phase-sequence fundametal compoment.
As shown in Figure 1 to 4 the extracting method of positive and negative sequence component and 5 order harmonic components, with this rotation transformation side is illustrated
Method principle, λ represents the voltage magnitude gain after line voltage rotation in Fig. 1～Fig. 4, and θ represents phase after line voltage rotation
The skew of position.Fig. 1～Fig. 3 represents respectively negative sequence component, positive-sequence component, 5 subharmonic point obtained by the rotation transformation methods
The polar plot of amount, Fig. 4 represents the current period line voltage synthesized by the negative sequence component, positive-sequence component, 5 order harmonic components that obtain
Polar plot.
Fig. 5～Fig. 8 gives the experimental result for extracting line voltage positive and negative sequence component and 5 order harmonic components, and Fig. 5～
Fig. 7 sets forth 30 degree of three-phase power grid voltage phase place impact, BC two-phase voltages and fall 50%, three-phase power grid voltage injection amplitude
For the experimental result of 5 subharmonic of 100V, Fig. 8 is these three situations while the experimental result being added on normal electrical network.From experiment
As a result from the point of view of, asymmetric and during containing harmonic wave in line voltage, grid phase information stable and accurate within 2ms can be obtained
, embody the rapidity and accuracy of method proposed by the present invention.
The preferred embodiment of the invention is the foregoing is only, it is all at this not to limit the invention
Any modification, equivalent and improvement made within the spirit and principle of innovation and creation etc., should be included in the invention
Protection domain within.
Claims (3)
1. the extracting method of each component instantaneous value of a kind of three phase network based on rotation transformation, it is characterised in that：By α β coordinate systems
Under line voltage by a rotational transformation matrix, obtain zero-sequence component, positive sequence fundametal compoment, negative phase-sequence under α β coordinate systems
Fundametal compoment and each harmonic component, rotational transformation matrix is as follows：
Wherein：J-the imaginary part of symbol；K=0, ± 1, ± 2 ..., -50≤k≤50；ω-line voltage fundamental wave angular frequency；TS- electrical network
Voltage cycle.
2. the extracting method of each component instantaneous value of a kind of three phase network based on rotation transformation according to claim 1, its
It is characterised by：
The rotational transformation matrix R is obtained according to the following steps：
First, converted according to CLARKE, the line voltage under three-phase ABC coordinate systems is transformed to into α β coordinate systems, obtain α β coordinates
Line voltage instantaneous value under system：
In formula：Emk--- the amplitude of line voltage k component of degree n ns；--- the starting phase angle of line voltage k component of degree n ns；
ω --- line voltage fundamental wave angular frequency；T --- the line voltage under three-phase ABC coordinate systems is transformed to into the time of α β coordinate systems
Point；E --- the line voltage vector under three-phase ABC coordinate systems；eɑ--- the voltage after line voltage α β coordinate transforms on α axles
Value；eβ--- the magnitude of voltage after line voltage α β coordinate transforms on β axles；J --- the imaginary part of symbol；K=0, ± 1, ± 2 ..., -50≤
k≤50；
Step 2, by the line voltage instantaneous value discretization obtained by coordinate transform, such as following formula：
In formula：TS--- the line voltage cycle；Em0--- the zero-sequence component amplitude of line voltage；--- line voltage zero sequence point
The phase angle of amount；E-mk--- the peak value of line voltage-k component of degree n ns；--- the starting phase angle of line voltage-k component of degree n ns；
N is constant；The line voltage vector of e (n) --- current period, e (n-1) --- the line voltage in current period previous cycle
Vector, e (n+1) --- the line voltage vector of a cycle after current period；
Step 3, complex column vector Ε (n) that the line voltage in front l cycle is expressed as an exponential form,
The instantaneous value of the various components of current period line voltage, is expressed as complex column vector ε (n) of an exponential form, such as
Following formula：
In formula：e-mk(n),emk- k time, k component of degree n ns of (n) --- current period line voltage, k=0,1,2 ...；e-m1(n),em1
The negative phase-sequence of (n) --- current period line voltage, positive sequence fundametal compoment；em0(n) --- the zero sequence of current period line voltage point
Amount；
Mark：The line voltage vector of e (n) --- current period, the electrical network electricity in e (n-1) --- current period previous cycle
Pressure vector, e (n+1) --- the line voltage vector of a cycle after current period.
Step 4, according to Ε (n) and the expression formula of ε (n), obtain a spin matrix R：
As l=2k+1, the various components of current electric grid voltage are further calculated：
ε (n)=R-1Ε(n)。
3. the extracting method of each component instantaneous value of a kind of three phase network based on rotation transformation according to claim 1 and 2,
It is characterized in that：The instantaneous positive sequence angle and instantaneous negative phase-sequence angle of electrical network are further calculated by the positive and negative sequence component for obtaining：
In formula：θ+,θ---- the instantaneous phase angle of positive and negative sequence fundametal compoment；
Re(em1(n)),Im(em1(n)) --- the real part and imaginary part of positive sequence fundametal compoment；
Re(e-m1(n)),Im(e-m1(n)) --- the real part and imaginary part of negative phase-sequence fundametal compoment.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610860915.9A CN106645944B (en) | 2016-09-28 | 2016-09-28 | A kind of extracting method of each component instantaneous value of three phase network based on rotation transformation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610860915.9A CN106645944B (en) | 2016-09-28 | 2016-09-28 | A kind of extracting method of each component instantaneous value of three phase network based on rotation transformation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106645944A true CN106645944A (en) | 2017-05-10 |
| CN106645944B CN106645944B (en) | 2019-01-22 |
Family
ID=58853910
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610860915.9A Active CN106645944B (en) | 2016-09-28 | 2016-09-28 | A kind of extracting method of each component instantaneous value of three phase network based on rotation transformation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106645944B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109256789A (en) * | 2018-10-23 | 2019-01-22 | 深圳供电局有限公司 | A kind of method of three-phase imbalance regulating device and its current limliting |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS62249285A (en) * | 1986-04-23 | 1987-10-30 | Hitachi Ltd | Rotation arithmetic unit |
| CN101509945A (en) * | 2009-02-27 | 2009-08-19 | 国网电力科学研究院 | Real-time detection method for positive and negative sequence electricity quantity |
| CN101534015A (en) * | 2009-04-10 | 2009-09-16 | 保定天威集团有限公司 | Three-phase photovoltaic grid-connected inverting control method and device thereof |
| CN103647287A (en) * | 2013-11-13 | 2014-03-19 | 国网上海市电力公司 | Dynamic voltage reactive compensation method |
-
2016
- 2016-09-28 CN CN201610860915.9A patent/CN106645944B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS62249285A (en) * | 1986-04-23 | 1987-10-30 | Hitachi Ltd | Rotation arithmetic unit |
| CN101509945A (en) * | 2009-02-27 | 2009-08-19 | 国网电力科学研究院 | Real-time detection method for positive and negative sequence electricity quantity |
| CN101534015A (en) * | 2009-04-10 | 2009-09-16 | 保定天威集团有限公司 | Three-phase photovoltaic grid-connected inverting control method and device thereof |
| CN103647287A (en) * | 2013-11-13 | 2014-03-19 | 国网上海市电力公司 | Dynamic voltage reactive compensation method |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109256789A (en) * | 2018-10-23 | 2019-01-22 | 深圳供电局有限公司 | A kind of method of three-phase imbalance regulating device and its current limliting |
| CN109256789B (en) * | 2018-10-23 | 2021-07-02 | 深圳供电局有限公司 | Three-phase unbalance adjusting device and current limiting method thereof |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106645944B (en) | 2019-01-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101640534B (en) | Full digital phase-locked loop applying rapid frequency capture method | |
| CN107706929A (en) | Adaptive phase locked loop method and system based on minimum variance filtering | |
| CN103487652B (en) | Frequency self-adaptive real-time fractional harmonic wave detection method | |
| CN104181374B (en) | Method for detecting and separating positive and negative sequence components of grid voltage of three-phase neutral-line-free system | |
| CN101509945B (en) | Real-time detection method for positive and negative sequence electricity quantity | |
| CN109358223A (en) | A kind of sliding window DFT harmonic current detecting method and device, storage medium | |
| CN102401858A (en) | Method for detecting fundamental component and harmonic component of voltage of power grid | |
| CN104360137A (en) | Voltage sag detecting method suitable for dynamic voltage restorer | |
| CN105158592A (en) | Real-time three-phase electrical network impedance detection method on the asymmetric electrical network condition | |
| CN114123346B (en) | General DFF-SAI phase locking method under positive sequence rotation coordinate system | |
| McEwen et al. | A novel sampling theorem on the rotation group | |
| CN105021872A (en) | Active and passive component detection method for different current components in distorted and asymmetric voltage state of electrical network | |
| Rossi et al. | Two-fluid numerical simulations of turbulence inside Kelvin-Helmholtz vortices: Intermittency and reconnecting current sheets | |
| Mokeev | Filter synthesis for PMU | |
| CN102095915B (en) | Voltage signal detecting device using multiple synchronous reference coordinate system transformation | |
| CN104833853A (en) | Frequency-adaptive slide window DFT harmonic wave detection method | |
| CN103078633A (en) | Single-phase phase-locked loop | |
| CN103543331B (en) | A kind of method calculating electric signal harmonic wave and m-Acetyl chlorophosphonazo | |
| CN101820281A (en) | Single-phase phase locked loop based on double-park transformation phase discriminator and implement method thereof | |
| Pisarski et al. | Banks of templates for all-sky narrow-band searches of gravitational waves from spinning neutron stars | |
| CN106645944A (en) | Extraction method for instantaneous value of each component of three-phase power grid based on rotation transformation | |
| CN106483375B (en) | A kind of multi-frequency fractional harmonic wave detection method | |
| CN104868909A (en) | Floating frequency and phase lock loop based on voltage quadrature resonator (QR) and measuring method thereof | |
| CN102520246B (en) | Constant frequency phasor extraction method | |
| CN105067882A (en) | Power amplitude measuring method |
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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20190408 Address after: 211100 No. 155 General Avenue, Jiangning Economic Development Zone, Nanjing City, Jiangsu Province Patentee after: Nanjing Estun Automation Co., Ltd. Address before: 211100 No. 155 General Avenue, Jiangning Economic Development Zone, Nanjing City, Jiangsu Province Patentee before: Nanjing Estun Automatic Control Technology Co., Ltd. |
|
| TR01 | Transfer of patent right |