CN108957130A - A kind of trend of harmonic detection method of power based on instantaneous reactive power theory - Google Patents
A kind of trend of harmonic detection method of power based on instantaneous reactive power theory Download PDFInfo
- Publication number
- CN108957130A CN108957130A CN201810844520.9A CN201810844520A CN108957130A CN 108957130 A CN108957130 A CN 108957130A CN 201810844520 A CN201810844520 A CN 201810844520A CN 108957130 A CN108957130 A CN 108957130A
- Authority
- CN
- China
- Prior art keywords
- harmonic
- phase
- sequence
- frequency
- negative
- 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)
- Measurement Of Current Or Voltage (AREA)
Abstract
The invention discloses a kind of trend of harmonic detection method of power based on instantaneous reactive power theory: quickly being analyzed with fft algorithm power network signal first, tentatively obtains each harmonic wave and m-Acetyl chlorophosphonazo frequency content；And then frequency is established respectively and is two cordic phase rotators of harmonic wave positive sequence and harmonic wave negative phase-sequence, and is filtered respectively again after power grid three-phase signal is transformed to the two rotating coordinate systems respectively；Three-phase positive sequence and Negative sequence harmonic component can be obtained finally by coordinate inversion.The present invention can detect positive sequence harmonic but also detect Negative sequence harmonic, while can also detect the harmonic content of each frequency band including m-Acetyl chlorophosphonazo, and structure is simple, is easy to Digital Realization, have very strong practical value.
Description
Technical field
The present invention relates to a kind of trend of harmonic detection method of power, are based on instantaneous reactive power theory more particularly, to one kind
Trend of harmonic detection method of power.
Background technique
In recent years, installed capacity of wind-driven power is continuously increased, and HVDC transmission line is also more and more, so that electric system
Power electronics degree is higher and higher.And as power electronic equipment is increasing in the power system, so that in electric system
Harmonic content become complex, while there is also a large amount of m-Acetyl chlorophosphonazo, these are all that harmonic detecting technique urgently solves
Certainly the problem of.
Currently used harmonic detecting method mainly include filter method, frequency analysis method, wavelet analysis analytic approach and
Using the harmonic detecting method of instantaneous reactive power theory.
Filter method algorithm structure is simple, however accuracy and reliability is poor and band logical and bandstop filter engineering all
In be not easy to realize.Frequency analysis method generally use discrete Fourier transform (DFT) or Fast Fourier Transform (FFT) (FFT) to harmonic wave into
Row analysis, although can detecte out the amplitude and phase of each harmonic component, due to the limitation of its window width, real-time compared with
It is low, it applies in general to carry out off-line calculation, often precision and reliability are lower when for real-time harmonic detection.Based on wavelet transformation
Harmonic detecting method can detecte out fundamental active, reactive current and total harmonic current, but each harmonic content cannot be detected,
And the harmonic detecting method precision and reliability of wavelet transformation are also all poor.Using the harmonic detecting side of instantaneous reactive power theory
Method accurately can reliably detect system based on instantaneous reactive power theory, through coordinate transform and corresponding calculating
Harmonic content.
But the existing conventional harmonic detection algorithm based on instantaneous reactive power theory can only also detect that total harmonic wave contains
The content and phase angle that measures and cannot detect each harmonic, will no longer in certain harmonic waves for needing to obtain specific frequency and phase sequence
It is applicable in.
Summary of the invention
The technical problems to be solved by the invention exactly propose a kind of electric system based on instantaneous reactive power theory
Harmonic detecting method can detect positive sequence harmonic but also detect Negative sequence harmonic, while can also detect each frequency including m-Acetyl chlorophosphonazo
The harmonic content of rate section.
Above-mentioned technical problem is solved, the technical solution adopted in the present invention is as follows.
A kind of trend of harmonic detection method of power based on instantaneous reactive power theory, it is characterised in that include following step
It is rapid:
(1) to power grid three-phase current or voltage signal Xs(s=a, b, c) is sampled, its discrete sample signals is obtained；
(2) discrete sample signals X described to step (1)sIt is quickly analyzed using fft algorithm, tentatively obtains each harmonic wave
And m-Acetyl chlorophosphonazo frequency content；
(3) by XsTransforming to frequency respectively is fpPositive sequence and negative phase-sequence two-phase rotating coordinate system in, obtain component Xd、Xq；
Wherein fpFor harmonic frequency of interest in harmonic wave acquired in step (2) and m-Acetyl chlorophosphonazo；
(4) to step (3) Xd、XqLow-pass filtering is carried out, is obtained
(5) rightProgress frequency is fpPositive sequence and negative phase-sequence coordinate inversion, obtain fpThree-phase positive sequence and negative phase-sequence
Harmonic component.
The present invention further comprises following preferred embodiment:
Step (3) specifically includes the following contents:
Three-phase current or voltage are carried out abc/ α β coordinate system transformation by Clarke transform matrix (1) by 3.1:
Transforming to frequency is fpWhen the two-phase rotating coordinate system of positive sequence, phase angle theta is taken as 2 π fpT, transforming to frequency is fpNegative phase-sequence
Two-phase rotating coordinate system when, phase angle theta is taken as -2 π fpt；θ
3.2 carry out the coordinate system transformation of α β/dq by Park Transformation matrix (2).
Step (4) low-pass filtering is using low-pass first order filter, transmission function Glp(s) expression formula such as formula
(3):
In formula, ωcFor the cutoff frequency of low-pass first order filter, value is [0Hz, 50Hz], and s is Laplace operator.
In step (5), coordinate inversion is carried out according to transformation matrix (1), (2), (4), obtains three-phase positive sequence and negative phase-sequence
Harmonic component.
Cdq/abc=(Cαβ/dqCabc/αβ)-1 (4)
Wherein transformation matrix (1), (2) are the same as described in 3.1 and 3.2.
This method first quickly analyzes power network signal with fft algorithm, tentatively obtains each harmonic wave and m-Acetyl chlorophosphonazo frequency
Ingredient；And then frequency is established respectively and is two cordic phase rotators of harmonic wave positive sequence and harmonic wave negative phase-sequence, and power grid three-phase signal is distinguished
It is filtered respectively again after transforming to the two rotating coordinate systems；Three-phase positive sequence can be obtained finally by coordinate inversion
With Negative sequence harmonic component.
The present invention can obtain compared with the existing technology it is following the utility model has the advantages that
The present invention is based on the basic principles of instantaneous reactive power theory to propose a kind of novel Harmonic Detecting Algorithm, can
Detection positive sequence harmonic can detect Negative sequence harmonic again, while can also detect the harmonic content of each frequency band including m-Acetyl chlorophosphonazo, and
Structure is simple, is easy to Digital Realization, has very strong practical value.
Detailed description of the invention
Fig. 1 show the trend of harmonic detection method of power schematic diagram of the invention based on instantaneous reactive power theory；
Fig. 2 show harmonic detecting method schematic diagram of the embodiment of the present invention under hypothesis testing signal；
Specific embodiment
Technical solution of the present invention is further stated in detail with reference to the accompanying drawings of the specification.
A kind of trend of harmonic detection method of power based on instantaneous reactive power theory disclosed by the invention such as 1 institute of attached drawing
Show, comprising the following steps:
(1) to power grid three-phase current or voltage signal Xs(s=a, b, c) is sampled, its discrete sample signals is obtained；
(2) discrete sample signals X described to step (1)sIt is quickly analyzed using fft algorithm, tentatively obtains each harmonic wave
And m-Acetyl chlorophosphonazo frequency content；
(3) by XsTransforming to frequency respectively is fpIn the two-phase rotating coordinate system of positive sequence and negative phase-sequence, component X is obtainedd、Xq；
Wherein fpFor harmonic frequency of interest in harmonic wave acquired in step (2) and m-Acetyl chlorophosphonazo；
(4) to step (3) Xd、XqLow-pass filtering is carried out, is obtained
(5) rightProgress frequency is fpThe coordinate inversion of positive sequence and negative phase-sequence, obtains fpThree-phase positive sequence and negative phase-sequence
Harmonic component.
Step (3) includes the following contents:
Three-phase current/voltage is carried out abc/ α β coordinate system transformation by Clarke transform matrix (1) by 3.1:
Transforming to frequency is fpWhen the two-phase rotating coordinate system of positive sequence, it is 2 π f that phase, which takes θ,pT, transforming to frequency is fpNegative phase-sequence
Two-phase rotating coordinate system when, phase take θ be -2 π fpt；
3.2 carry out α β/dq coordinate system transformation by Park Transformation matrix (2):
Step (4) described low-pass filtering, using low-pass first order filter, transmission function Glp(s) expression formula such as formula
(3):
In formula, ωcFor the cutoff frequency of low-pass first order filter, value is [0Hz, 50Hz], and the embodiment of the present application is preferred
10Hz；S is Laplace operator.
In step (5), coordinate inversion is carried out according to transformation matrix (1), (2), (4), is finally obtaining to obtain three-phase just
Sequence and Negative sequence harmonic component.
Cdq/abc=(Cαβ/dqCabc/αβ)-1 (4)
Wherein transformation matrix (1), (2) are the same as described in 3.1 and 3.2.
Embodiment below with reference to Fig. 2 and table 1 illustrates the present invention.
The embodiment of the present invention tests signal using three-phase current, which includes the fundamental frequency of frequency 50Hz amplitude 100A
Forward-order current, the frequency tripling positive-sequence component of frequency 150Hz amplitude 10A, the fifth harmonic negative sequence component of frequency 250Hz amplitude 5A and
The negative phase-sequence m-Acetyl chlorophosphonazo component of frequency 70Hz amplitude 4A；The initial phase of each harmonic is zero.
Using harmonic detecting method as shown in Figure 2, harmonic detecting result as shown in table 1 is obtained.By 1 result of table as it can be seen that
Harmonic wave algorithm of the present invention can accurately detect positive sequence and Negative sequence harmonic under each frequency band, and can not only detect integral frequency harmonizing wave can also
M-Acetyl chlorophosphonazo is detected, traditional defect based on instantaneous reactive power theory Harmonic Detecting Algorithm is compensated for.
Table 1: harmonic detecting result of the embodiment of the present invention under hypothesis testing signal
Claims (4)
1. a kind of trend of harmonic detection method of power based on instantaneous reactive power theory, it is characterised in that the following steps are included:
(1) to Three-phase Power Systems current or voltage signal Xs(s=a, b, c) is sampled, its discrete sample signals is obtained；
(2) discrete sample signals X described to step (1)sIt is quickly analyzed using fft algorithm, tentatively obtains each harmonic wave and humorous
Wave frequency rate ingredient；
(3) by XsTransforming to frequency respectively is fpIn the two-phase rotating coordinate system of positive sequence and negative phase-sequence, component X is obtainedd、Xq；
Wherein fpFor in harmonic wave acquired in step (2) and m-Acetyl chlorophosphonazo, the harmonic frequency of care；
(4) to step (3) Xd、XqLow-pass filtering is carried out, is obtained
(5) rightProgress frequency is fpThe coordinate inversion of positive sequence and negative phase-sequence, obtains fpThree-phase positive sequence and Negative sequence harmonic
Component.
2. the trend of harmonic detection method of power according to claim 1 based on instantaneous reactive power theory, feature exist
The following contents is specifically included in the step (3):
Three-phase current or voltage are carried out abc/ α β coordinate system transformation, transformation matrix by Clarke transform matrix (1) by 3.1 are as follows:
When transforming to the two-phase rotating coordinate system that frequency is fp positive sequence, taking phase angle theta is 2 π fpT transforms to two that frequency is fp negative phase-sequence
When cordic phase rotator system, taking phase angle theta is -2 π fpt；
3.2 carry out α β/dq coordinate system transformation by Park Transformation matrix (2):
3. the trend of harmonic detection method of power according to claim 1 based on instantaneous reactive power theory, feature exist
In: the low-pass filtering of the step (4) is using low-pass first order filter, transmission function expression formula such as formula (3):
In formula, ωcFor the cutoff frequency of low-pass first order filter, value is [0Hz, 50Hz].
4. a kind of trend of harmonic detection method of power based on instantaneous reactive power theory according to claim 2, special
Sign is: in the step (5), carrying out coordinate inversion according to transformation matrix (1), (2), (4), finally obtains three-phase positive sequence
With Negative sequence harmonic component.
Cdq/abc=(Cαβ/dqCabc/αβ)-1 (4)
Wherein, transformation matrix (1), (2) are the same as described in 3.1 and 3.2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810844520.9A CN108957130B (en) | 2018-07-27 | 2018-07-27 | Electric power system harmonic detection method based on instantaneous reactive power theory |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810844520.9A CN108957130B (en) | 2018-07-27 | 2018-07-27 | Electric power system harmonic detection method based on instantaneous reactive power theory |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108957130A true CN108957130A (en) | 2018-12-07 |
| CN108957130B CN108957130B (en) | 2021-01-08 |
Family
ID=64465829
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810844520.9A Expired - Fee Related CN108957130B (en) | 2018-07-27 | 2018-07-27 | Electric power system harmonic detection method based on instantaneous reactive power theory |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108957130B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110531138A (en) * | 2019-07-08 | 2019-12-03 | 江苏科技大学 | A kind of Active Power Filter Harmonic Currents detection method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090237964A1 (en) * | 2006-09-15 | 2009-09-24 | Abb Schweiz Ag | Method for operating a converter circuit and apparatus for implementing the method |
| CN105162137A (en) * | 2015-09-12 | 2015-12-16 | 南昌航空大学 | Reactive and harmonic current detection algorithm in complex power network |
| CN105182106A (en) * | 2015-07-27 | 2015-12-23 | 国家电网公司 | Three-phase phase-controlled rectifier harmonic characteristic analysis method based on Park transformation |
| CN106771597A (en) * | 2017-03-03 | 2017-05-31 | 燕山大学 | A kind of no phase-locked loop formula SAI harmonic current detecting methods of omission Park conversion |
-
2018
- 2018-07-27 CN CN201810844520.9A patent/CN108957130B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090237964A1 (en) * | 2006-09-15 | 2009-09-24 | Abb Schweiz Ag | Method for operating a converter circuit and apparatus for implementing the method |
| CN105182106A (en) * | 2015-07-27 | 2015-12-23 | 国家电网公司 | Three-phase phase-controlled rectifier harmonic characteristic analysis method based on Park transformation |
| CN105162137A (en) * | 2015-09-12 | 2015-12-16 | 南昌航空大学 | Reactive and harmonic current detection algorithm in complex power network |
| CN106771597A (en) * | 2017-03-03 | 2017-05-31 | 燕山大学 | A kind of no phase-locked loop formula SAI harmonic current detecting methods of omission Park conversion |
Non-Patent Citations (5)
| Title |
|---|
| 张杰成 等: "滑窗迭代DFT检测谐波和无功电流的新算法", 《电力系统及其自动化学报》 * |
| 王文新 等: "一种用于有源滤波器的改进谐波检测法", 《科学技术与工程》 * |
| 赵磊: "基于Prony的谐波间谐波检测方法和检测系统的研究与应用", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
| 陈明亮 等: "一种完全谐波抑制正负序同步信号检测方法", 《电力电子技术》 * |
| 陈雷: "力系统谐波检测与补偿方法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110531138A (en) * | 2019-07-08 | 2019-12-03 | 江苏科技大学 | A kind of Active Power Filter Harmonic Currents detection method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108957130B (en) | 2021-01-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101793918B (en) | Voltage sag detection method | |
| CN103869162B (en) | Dynamic signal phasor measurement method based on time domain quasi-synchronization | |
| CN102621382B (en) | Method for detecting frequency, phase and amplitude of electric signal in electric power system | |
| CN101893651B (en) | Method for detecting positive sequence, negative sequence, idle and harmonic currents of power supply system | |
| CA2929382C (en) | Method of removing decaying dc component from power system fault signal | |
| CN104181374B (en) | Method for detecting and separating positive and negative sequence components of grid voltage of three-phase neutral-line-free system | |
| CN102385003B (en) | Method for detecting frequency and phase of voltage signals or current signals in electric power system | |
| CN103983889B (en) | Diagnosis method for active power filter switching device open-circuit fault based on model reference analysis | |
| CN105162137B (en) | Idle and harmonic current detection under complex electric network | |
| CN101509945B (en) | Real-time detection method for positive and negative sequence electricity quantity | |
| CN107632200B (en) | Harmonic detection method based on discrete Fourier transform | |
| CN104155520A (en) | Harmonic and inter-harmonic separation detecting method and device suitable for vessel electric power system | |
| CN105823921A (en) | Compensating current detection method based on instant space voltage vector orientation | |
| CN102221639A (en) | Positive and negative sequence current real-time detection method | |
| CN102062803B (en) | Active current extraction method based on orthogonal sine wave integral modulation | |
| CN109061345B (en) | Effective value measuring method and system suitable for power system | |
| CN106597217A (en) | MMC-HVDC alternating current side asymmetric fault diagnosis method | |
| CN104833853A (en) | Frequency-adaptive slide window DFT harmonic wave detection method | |
| CN105021872A (en) | Active and passive component detection method for different current components in distorted and asymmetric voltage state of electrical network | |
| CN108414838A (en) | A kind of inverter parallel system line impedance measurement method | |
| CN103149436B (en) | The m-Acetyl chlorophosphonazo detection system of three-phase active power distribution network | |
| CN108957130A (en) | A kind of trend of harmonic detection method of power based on instantaneous reactive power theory | |
| CN103546149A (en) | Phase locking method for three-phase power system | |
| CN103364634B (en) | Three-phase alternating-current supply frequency measurement method | |
| CN103176030A (en) | Method for detecting inter-harmonics of power distribution system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210108 Termination date: 20210727 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |