CN103269211B - The digital filtering method of three-phase alternating current system measurement data - Google Patents

The digital filtering method of three-phase alternating current system measurement data Download PDF

Info

Publication number
CN103269211B
CN103269211B CN201310157826.4A CN201310157826A CN103269211B CN 103269211 B CN103269211 B CN 103269211B CN 201310157826 A CN201310157826 A CN 201310157826A CN 103269211 B CN103269211 B CN 103269211B
Authority
CN
China
Prior art keywords
phase
voltage signal
under
base
alternating current
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201310157826.4A
Other languages
Chinese (zh)
Other versions
CN103269211A (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.)
Electric Power Research Institute of Guangdong Power Grid Co Ltd
Original Assignee
Wuhan University WHU
Electric Power Research Institute of Guangdong 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 Wuhan University WHU, Electric Power Research Institute of Guangdong Power Grid Co Ltd filed Critical Wuhan University WHU
Priority to CN201310157826.4A priority Critical patent/CN103269211B/en
Publication of CN103269211A publication Critical patent/CN103269211A/en
Application granted granted Critical
Publication of CN103269211B publication Critical patent/CN103269211B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The digital filtering method of a kind of three-phase alternating current system measurement data, including step: by three-phase voltage signal VA、VB、VCBe converted to the voltage signal V under biphase α β coordinate systemα、Vβ;According to the transmission function set respectively to voltage signal Vα、VβCarry out low-pass filtering and obtain voltage signal V'α、V'β;Calculate described voltage signal V'α、V'βBase Ld phase δ;The base Ld phase delay angle δ of low-pass filtering is calculated according to transferometer0;According to base Ld phase δ by voltage signal V'α、V'βBe converted to rotate the DC voltage V under dq two phase coordinate systemd、Vq;With δ δ0For angle of transformation by DC voltage Vd、VqBe converted to the voltage signal V' under three phase coordinatesA、V'B、V'C.The digital filtering method of above-mentioned three-phase alternating current system measurement data, can efficiently solve the base Ld phase delay issue of tradition filtering mode, improves the phase controlling precision of power system.

Description

The digital filtering method of three-phase alternating current system measurement data
Technical field
The present invention relates to digital filtering technique field, particularly relate to the numeral filter of a kind of three-phase alternating current system measurement data Wave method.
Background technology
Along with high-power switch device is in the application in the fields such as wind-power electricity generation, photovoltaic generation, Survey of Flexible AC Transmission System, power train System presents more and more significant non-linear, especially applies voltage source converter (VSC) based on pulsewidth modulation (PWM) technology Occasion, high order characteristic harmonics amplitude is close, possibly even beyond fundamental frequency amplitude, in order to solve the problems referred to above, can take be System primary side installs the scheme of wave filter additional, but due to the popularity of higher hamonic wave frequency spectrum, uncertainty, and filter apparatus cost Higher, the obtained effect of this technology is the most undesirable.
The wave distortion caused due to higher hamonic wave makes system normally to control and accurately measures, therefore, and existing electricity Force system typically can carry out the digital filtering of single order or high-order lowpass in the measurement links of secondary system side to sampled data, currently Main research is all the parameter designing being transfer function H (s).But, actual central any physically realizable system all can be deposited Postponing, no matter how the coefficient in transfer function H (s) designs, three-phase alternating current system measurement data after filtering after, be all difficult to Avoiding filtering forward and backward base Ld phase and postpone situation, three-phase alternating current system measurement data percent harmonic distortion (THD) is the highest, There is severe deviations in the phase controlling making power system.
Summary of the invention
Based on this, it is necessary to for the problems referred to above, it is provided that the digital filtering method of a kind of three-phase alternating current system measurement data.
The digital filtering method of a kind of three-phase alternating current system measurement data, comprises the steps:
By three-phase voltage signal VA、VB、VCBe converted to the voltage signal V under biphase α β coordinate systemα、Vβ
According to the transmission function set respectively to described voltage signal Vα、VβCarry out low-pass filtering and obtain voltage signal V'α、 V'β
Calculate described voltage signal V'α、V'βBase Ld phase δ;
The base Ld phase delay angle δ of described low-pass filtering is calculated according to described transferometer0
According to described base Ld phase δ by described voltage signal V'α、V'βBe converted to rotate the unidirectional current under dq two phase coordinate system Pressure Vd、Vq
With δ-δ0For angle of transformation by described DC voltage Vd、VqBe converted to the voltage signal V' under three phase coordinatesA、V'B、V'C
The digital filtering method of above-mentioned three-phase alternating current system measurement data, first sits three-phase voltage signal from ABC three-phase Mark is transformed to α β two phase coordinates, then carries out low-pass filtering under α β coordinate system, and filtered fundamental frequency value is transformed to by α β coordinate DC quantity under dq rotating coordinate system, and base Ld phase produced by low-pass filtering postpones δ0Sit being returned ABC by dq coordinate transform Timestamp compensates, thus efficiently solves the base Ld phase delay issue of tradition filtering mode, improves the phase place of power system Control accuracy.
Accompanying drawing explanation
Fig. 1 is the digital filtering method flow chart of the three-phase alternating current system measurement data of an embodiment;
Fig. 2 is simulation example laboratory model construction schematic diagram;
Fig. 3 is the simulation experiment result schematic diagram.
Detailed description of the invention
Specific embodiment party to the digital filtering method of the three-phase alternating current system measurement data of the present invention below in conjunction with the accompanying drawings Formula is described in detail.
Fig. 1 shows the digital filtering method flow chart of the three-phase alternating current system measurement data of an embodiment, including such as Lower step:
Step S101: by three-phase voltage signal VA、VB、VCBe converted to the voltage signal V under biphase α β coordinate systemα、Vβ
In one embodiment, also include before step S101: the measurement data of three-phase alternating current system is carried out sampling and obtains Obtain the voltage signal V under three phase coordinatesA、VB、VC
Specifically, the voltage data V that sampling obtains will be measured in step S101A、VB、VC, from static ABC tri-phase coordinates Being transformed to static α β two phase coordinates, the formula of conversion is as follows:
Wherein, ABC → α β coordinate is transformed to orthogonal transformation, arranges α phase and A same axis.
Step S102: according to the transmission function set respectively to described voltage signal Vα、VβCarry out low-pass filtering and obtain voltage Signal V'α、V'β
Specifically, calculating process includes below equation:
In formula, H (s) is the transmission function set,
Step S103: calculate described voltage signal V'α、V'βBase Ld phase δ.
Specifically, calculating process includes below equation:
δ=arctg(V'β/V'α);
In formula, δ is base Ld phase, V'α、V'βIt is respectively the voltage signal of α β coordinate after filtering.
Step S104: calculate the base Ld phase delay angle δ of described low-pass filtering according to described transferometer0
Specifically, phase delay angle δ that angle is low-pass filtering link of transfer function H (j ω)0, i.e. through low pass filtered After ripple, under biphase α β coordinate system, the base Ld phase δ of voltage signal postpones relative to the voltage signal of A, B, C three-phase under three phase coordinates δ0, δ0Calculating process includes below equation:
δ 0 = arctg ( ImH ( jω ) ReH ( jω ) ) ;
The transmission function of described low-pass filtering is H (j ω), and its formula can be expressed as:
H ( jω ) = a 0 + a 1 jω b 0 + b 1 jω - b 2 ω 2 ;
Wherein, a0、a1、b0、b1、b2For default constant, equal to 1, ω is fundamental frequency angular frequency to the modulus value of H (j ω) | H (j ω) | Rate.
Step S105: according to described base Ld phase δ by described voltage signal V'α、V'βBe converted to rotate dq two phase coordinate system Under DC voltage Vd、Vq
Specifically, calculating process includes below equation:
Vd+jVq=(Vα+jVβ)e-jδ
In formula, V'α、V'βBeing respectively the voltage signal of α β coordinate after filtering, δ is base Ld phase, Vd、VqIt is respectively dq to rotate The voltage signal of coordinate system.
Step S106: with δ-δ0For angle of transformation by described DC voltage Vd、VqBe converted to the alternating voltage letter under three phase coordinates Number V'A、V'B、V'C
Specifically, variable is returned static ABC tri-phase coordinates from rotating the biphase coordinate transform of dq, q axle advanced d axle is set, become Changing angle is δ-δ0, its calculating process includes the following:
In formula, θ=δ-δ0, by above-mentioned processing procedure, during low-pass filtering, produced base Ld phase postpones δ0At dq → ABC coordinate transform link is compensated.
The digital filtering method of the three-phase alternating current system measurement data of the present invention, for the measurement data of three-phase alternating current system Process, first the three-phase voltage signal (containing higher hamonic wave) of sampling is transformed to α β two phase coordinates from ABC three-phase coordinate, then Carrying out low-pass filtering under α β coordinate system, filtered fundamental frequency value is transformed to the DC quantity under dq rotating coordinate system by α β coordinate, And base Ld phase produced by low-pass filtering postpones δ0Compensate when being returned ABC coordinate by dq coordinate transform, thus efficient solution The certainly base Ld phase delay issue of tradition filtering mode.
In order to become apparent from technical scheme, illustrate a simulation example below in conjunction with the accompanying drawings.
Shown in Figure 2, Fig. 2 is simulation example laboratory model construction schematic diagram, and wherein the transmission function of low-pass filtering is adopted With simplest first-order filtering H (j ω)=G/ (1+Ts), parameter G=1.606, T=0.004 second.
The ABC three-phase voltage signal then inputted is VA、VB、VC, fundamental frequency (50Hz) phase voltage virtual value is 1.0,7 subharmonic Content be 10%, 19 subharmonic content be 180%, total harmonic distortion factor (THD) is 180.278%, after filtering output voltage signal For V'A、V'B、V'C
As it is shown on figure 3, Fig. 3 is the simulation experiment result schematic diagram, V in figureAFor input A phase voltage signal waveform, wherein, V'AFor the A phase voltage signal waveform of output, V' after the filtering method of the present invention is filteredAFundamental frequency virtual value be 1.0, And phase place and VAIdentical, 7 subharmonic content be 1.812%, 19 subharmonic content be 12.007%, overall percent harmonic distortion (THD) It is 12.143%.
Even if by the experimental result of above-mentioned simulation example it can be seen that transfer function H (s) uses the most optimized and the simplest Single first-order low-pass wave function, the filtering method of the present invention also can make total harmonic distortion factor (THD) reduce by more than 90%, and base The amplitude of frequency and phase place keep constant.
Embodiment described above only have expressed the several embodiments of the present invention, and it describes more concrete and detailed, but also Therefore the restriction to the scope of the claims of the present invention can not be interpreted as.It should be pointed out that, for those of ordinary skill in the art For, without departing from the inventive concept of the premise, it is also possible to make some deformation and improvement, these broadly fall into the guarantor of the present invention Protect scope.Therefore, the protection domain of patent of the present invention should be as the criterion with claims.

Claims (5)

1. the digital filtering method of a three-phase alternating current system measurement data, it is characterised in that comprise the steps:
By three-phase voltage signal VA、VB、VCBe converted to the voltage signal V under biphase α β coordinate systemα、Vβ
According to the transmission function set respectively to described voltage signal Vα、VβCarry out low-pass filtering and obtain voltage signal V'α、V'β
Calculate described voltage signal V'α、V'βBase Ld phase δ;
The base Ld phase delay angle δ of described low-pass filtering is calculated according to described transferometer0;Including:
δ 0 = a r c t g ( Im H ( j ω ) Re H ( j ω ) ) ;
Wherein, H (j ω) is transmission function, and its formula can be expressed as:
H ( j ω ) = a 0 + a 1 j ω b 0 + b 1 j ω - b 2 ω 2 ;
Wherein, a0、a1、b0、b1、b2For default constant, equal to 1, ω is fundamental frequency angular frequency to the modulus value of H (j ω) | H (j ω) |, j Imaginary part for plural number;
According to described base Ld phase δ by described voltage signal V'α、V'βBe converted to rotate the DC voltage under dq two phase coordinate system Vd、Vq
With δ-δ0For angle of transformation by described DC voltage Vd、VqBe converted to the voltage signal V' under three phase coordinatesA、V'B、V'C
The digital filtering method of three-phase alternating current system measurement data the most according to claim 1, it is characterised in that described By three-phase voltage signal VA、VB、VCBe converted to the voltage signal V under biphase α β coordinate systemα、VβStep also include:
The measurement data of three-phase alternating current system is carried out the voltage signal V under sampling acquisition three phase coordinatesA、VB、VC
The digital filtering method of three-phase alternating current system measurement data the most according to claim 1, it is characterised in that described meter Calculate described voltage signal V'α、V'βThe step of base Ld phase δ include:
δ=arctg (V'β/V'α);
In formula, δ is base Ld phase, V'α、V'βIt is respectively the voltage signal under α β coordinate after filtering.
The digital filtering method of three-phase alternating current system measurement data the most according to claim 1, it is characterised in that described According to described base Ld phase δ by described voltage signal V'α、V'βBe converted to rotate the DC voltage V under dq two phase coordinate systemd、Vq's Step includes:
Vd+jVq=(Vα+jVβ)e-jδ
In formula, V'α、V'βBeing respectively the voltage signal under α β coordinate after filtering, δ is base Ld phase, Vd、VqIt is respectively the biphase seat of dq Voltage signal under Biao, j is the imaginary part of plural number.
The digital filtering method of three-phase alternating current system measurement data the most according to claim 1, it is characterised in that with δ-δ0 For angle of transformation by described DC voltage Vd、VqBe converted to the voltage signal V' under three phase coordinatesA、V'B、V'CStep include:
Wherein, θ=δ-δ0
CN201310157826.4A 2013-04-28 2013-04-28 The digital filtering method of three-phase alternating current system measurement data Active CN103269211B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310157826.4A CN103269211B (en) 2013-04-28 2013-04-28 The digital filtering method of three-phase alternating current system measurement data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310157826.4A CN103269211B (en) 2013-04-28 2013-04-28 The digital filtering method of three-phase alternating current system measurement data

Publications (2)

Publication Number Publication Date
CN103269211A CN103269211A (en) 2013-08-28
CN103269211B true CN103269211B (en) 2016-08-10

Family

ID=49012823

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310157826.4A Active CN103269211B (en) 2013-04-28 2013-04-28 The digital filtering method of three-phase alternating current system measurement data

Country Status (1)

Country Link
CN (1) CN103269211B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111308173B (en) * 2020-02-24 2022-05-03 广东海悟科技有限公司 Method, electronic device, and computer-readable storage medium for solving phase lag caused by three-phase current sampling filtering
CN112583319B (en) * 2020-12-02 2022-03-18 美的威灵电机技术(上海)有限公司 Phase voltage detection method and device of motor, electric appliance and readable storage medium

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101673952A (en) * 2009-08-14 2010-03-17 燕山大学 Precise phase locking method based on cross decoupling self-adaptive complex filter
CN102223100A (en) * 2011-06-17 2011-10-19 北京中能清源科技有限公司 Control method of three-phase grid-connected inverter based on modified proportional resonant regulator

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4297953B2 (en) * 2007-06-22 2009-07-15 三洋電機株式会社 Motor control device and compressor

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101673952A (en) * 2009-08-14 2010-03-17 燕山大学 Precise phase locking method based on cross decoupling self-adaptive complex filter
CN102223100A (en) * 2011-06-17 2011-10-19 北京中能清源科技有限公司 Control method of three-phase grid-connected inverter based on modified proportional resonant regulator

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
三相电压不平衡条件下锁相环的设计与实现;林百娟;《中国优秀硕士学位论文全文数据库》;20091215;正文第9页第2.1节及图2.2 *

Also Published As

Publication number Publication date
CN103269211A (en) 2013-08-28

Similar Documents

Publication Publication Date Title
CN103227581B (en) Inverter parallel harmonic wave ring current restraining method for controlling harmonic wave droop
CN104953606B (en) Networked layered compensation method for voltage unbalance of PCC (Point of Common Coupling) of islanded microgrid
CN102857133B (en) Current control method and current control system of single-phase single-stage photovoltaic inverter
CN102570476B (en) Repetitive-control-based method for controlling compensation current of DSTATCOM (Distribution Static Synchronous Compensator)
CN102097824B (en) Control method of LCL (Lower Control Limit) type grid-connected inverter
CN102904568B (en) Self-adaptive grid-tied converter single phase soft phase-locked loop
CN110739678B (en) Control method for series virtual impedance of grid-connected converter
CN104953801B (en) Voltage-source type combining inverter current harmonics elimination device and method
US9184652B2 (en) Method and apparatus for inverter output current harmonic reduction
CN104600748B (en) A kind of isolated island micro-capacitance sensor multi-inverter control system possessing active power filtering function and method of work thereof
CN106787910A (en) It is applied to the improvement repetitive controller method for designing of combining inverter current control
CN102801346B (en) Three-phase inverter with no-signal interconnecting lines connected in parallel and control method of three-phase inverter
CN105763094A (en) Inverter control method based on voltage feedforward and recombination current control
CN106712100A (en) Perturbation observation-based control method for virtual grid-connected synchronous inverter and control system
CN102628894B (en) Selective harmonic wave extracting method
CN102136738A (en) Control method of grid-connected inverter of large-scale grid-connected photovoltaic power station
CN103199539A (en) Voltage dip compensation method of unified power quality conditioner with zero input active power
CN106410858A (en) Software digital phase-locking method based on dual dq coordination conversion
CN104821601A (en) Three-phase photovoltaic grid-connected inverter control device
CN108599257B (en) Current control method suitable for high phase-locked loop bandwidth
CN105896539A (en) Phase sequence identification adaptive method applied to electric energy quality management device
CN103605038A (en) Protective error detecting and verifying system of electronic current transformer
Lee et al. Performance improvement of grid-connected inverter systems under unbalanced and distorted grid voltage by using a PR controller
CN102437753A (en) Three-ring control method of three-phase PWM (Pulse-Width Modulation) rectifier for LCL (Lower Control Limit) filtering
CN103269211B (en) The digital filtering method of three-phase alternating current system measurement data

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
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 510080 Dongfeng East Road, Dongfeng, Guangdong, Guangzhou, Zhejiang Province, No. 8

Co-patentee after: Wuhan University

Patentee after: ELECTRIC POWER RESEARCH INSTITUTE, GUANGDONG POWER GRID CO., LTD.

Address before: 510080 Dongfeng East Road, Dongfeng, Guangdong, Guangzhou, Zhejiang Province, No. 8

Co-patentee before: Wuhan University

Patentee before: Electrical Power Research Institute of Guangdong Power Grid Corporation

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210122

Address after: 510080 water Donggang 8, Dongfeng East Road, Yuexiu District, Guangzhou, Guangdong.

Patentee after: Electric Power Research Institute of Guangdong Power Grid Co.,Ltd.

Address before: 510080 water Donggang 8, Dongfeng East Road, Yuexiu District, Guangzhou, Guangdong.

Patentee before: Electric Power Research Institute of Guangdong Power Grid Co.,Ltd.

Patentee before: WUHAN University