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 PDF

Info

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
Application number
CN201810844520.9A
Other languages
Chinese (zh)
Other versions
CN108957130B (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.)
China Southern Power Grid Co Ltd
Original Assignee
China Southern Power Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Southern Power Grid Co Ltd filed Critical China Southern Power Grid Co Ltd
Priority to CN201810844520.9A priority Critical patent/CN108957130B/en
Publication of CN108957130A publication Critical patent/CN108957130A/en
Application granted granted Critical
Publication of CN108957130B publication Critical patent/CN108957130B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16Spectrum 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

A kind of trend of harmonic detection method of power based on instantaneous reactive power theory
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.
CN201810844520.9A 2018-07-27 2018-07-27 Electric power system harmonic detection method based on instantaneous reactive power theory Expired - Fee Related CN108957130B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
张杰成 等: "滑窗迭代DFT检测谐波和无功电流的新算法", 《电力系统及其自动化学报》 *
王文新 等: "一种用于有源滤波器的改进谐波检测法", 《科学技术与工程》 *
赵磊: "基于Prony的谐波间谐波检测方法和检测系统的研究与应用", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *
陈明亮 等: "一种完全谐波抑制正负序同步信号检测方法", 《电力电子技术》 *
陈雷: "力系统谐波检测与补偿方法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (1)

* Cited by examiner, † Cited by third party
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