CN102608422B - Method for computing alternating current component fundamental amplitude value in microcomputer protection - Google Patents

Method for computing alternating current component fundamental amplitude value in microcomputer protection Download PDF

Info

Publication number
CN102608422B
CN102608422B CN201210080890.2A CN201210080890A CN102608422B CN 102608422 B CN102608422 B CN 102608422B CN 201210080890 A CN201210080890 A CN 201210080890A CN 102608422 B CN102608422 B CN 102608422B
Authority
CN
China
Prior art keywords
frequency
voltage amplitude
current component
fundamental
fundamental 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.)
Active
Application number
CN201210080890.2A
Other languages
Chinese (zh)
Other versions
CN102608422A (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.)
NARI Group Corp
Nari Technology Co Ltd
Original Assignee
Nari Technology Co Ltd
Nanjing NARI Group Corp
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 Nari Technology Co Ltd, Nanjing NARI Group Corp filed Critical Nari Technology Co Ltd
Priority to CN201210080890.2A priority Critical patent/CN102608422B/en
Publication of CN102608422A publication Critical patent/CN102608422A/en
Application granted granted Critical
Publication of CN102608422B publication Critical patent/CN102608422B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Emergency Protection Circuit Devices (AREA)
  • Measurement Of Current Or Voltage (AREA)

Abstract

The invention discloses a method for computing alternating current component fundamental amplitude value in microcomputer protection, belonging to the microcomputer protection field of a power system, comprising the steps of obtaining the average value of the fundamental amplitude computed by Fourier series in the alternating current component 1/2 period, and then multiplying a coefficient related to frequency to obtain the actual fundamental amplitude value of the alternating current component. The method solves the problems in the prior art that as the alternating-current component fundamental amplitude value and the actual fundamental amplitude value of the alternating current component have a certain error caused by frequency shift, the more the frequency shift is, the bigger the error is, and the fundamental amplitude value error caused by the frequency shift can influence the protection action behavior of the microcomputer protection part. The method for computing alternating current component fundamental amplitude value in microcomputer protection can accurately compute the actual fundamental amplitude value of the alternating current component in the condition that the alternating current component fundamental frequency is not uniform with the rated frequency of the power system.

Description

The of ac fundamental voltage amplitude computing method used in a kind of Microcomputer Protection
Technical field
The invention belongs to electric system Microcomputer Protection field; relate to the of ac fundamental voltage amplitude computing method used in a kind of Microcomputer Protection, particularly use the fundamental voltage amplitude that calculates of Fourier series to average within of ac 1/2 cycle and be multiplied by one and ask the computing method of the actual fundamental voltage amplitude of of ac with the coefficient of frequency dependence.
Background technology
The general fundamental voltage amplitude using Fourier series to calculate of ac in current Microcomputer Protection.First calculate the coefficient of Fourier series according to the rated frequency of electric system and the sampling rate of device in Microcomputer Protection code and exist in a data form.These coefficients are always constant in the operational process of whole microcomputer protecting device.If the fundamental frequency of of ac and the rated frequency of electric system of input microcomputer protecting device have certain skew, so use these coefficient calculations fundamental voltage amplitude out will have certain error with actual fundamental voltage amplitude.
The first-harmonic composition of of ac can be expressed as x (t)=Msin (2 π ft+2 π fa), and wherein M, f, t, 2 π fa are respectively the fundamental voltage amplitude of of ac, the frequency of of ac, time, the initial phase angle of of ac.The computing formula utilizing Fourier series to ask for its fundamental voltage amplitude is:
F = | 2 f 0 ∫ - 1 / 2 f 0 1 / 2 f 0 M sin ( 2 πft + 2 πfa ) e - j 2 π f 0 t dt | - - - ( 1 )
F in formula (1) 0for the fundamental frequency chosen when calculating Fourier series.Abbreviation formula (1) can obtain:
F = M | 2 f 0 sin ( πf / f 0 ) / π / ( f 0 2 - f 2 ) | ( f 0 2 + f 2 ) / 2 + ( f 0 2 - f 2 ) / 2 * cos ( 4 πfa )
(2)
For moment t, the phase angle of of ac x (t) is 2 π ft+2 π fa, and therefore the pass of its fundamental voltage amplitude and time is:
A ( f , t ) = M | 2 f 0 sin ( πf / f 0 ) / π / ( f 0 2 - f 2 ) | ( f 0 2 + f 2 ) / 2 + ( f 0 2 - f 2 ) / 2 * cos ( 4 πft + 4 πfa ) - - - ( 3 )
The fundamental voltage amplitude calculated with the ratio of actual fundamental voltage amplitude is:
A b ( f , t ) = A ( f , t ) / M
= | 2 f 0 sin ( πf / f 0 ) / π / ( f 0 2 - f 2 ) | ( f 0 2 + f 2 ) / 2 + ( f 0 2 - f 2 ) / 2 * cos ( 4 πft + 4 πfa ) - - - ( 4 )
As can be seen from formula (3); when the fundamental frequency of of ac and the rated frequency of electric system that input microcomputer protecting device have certain skew time; the fundamental voltage amplitude using Fourier series to calculate is not equal to the fundamental voltage amplitude of of ac reality; but one with the function of of ac frequency and time correlation, the ratio of the fundamental voltage amplitude calculated and actual fundamental voltage amplitude is such as formula shown in (4).
In Microcomputer Protection, some protection operates in very large frequency range, and the fundamental voltage amplitude error caused due to frequency shift (FS) will affect the action behavior of this part protection.
Summary of the invention
In order to the fundamental voltage amplitude solving of ac fundamental voltage amplitude and of ac reality in prior art has certain error; frequency shift (FS) is more; error is larger; and the fundamental voltage amplitude error caused due to frequency shift (FS) will affect the problem of this computer part protection act behavior; the invention provides the of ac fundamental voltage amplitude computing method used in a kind of Microcomputer Protection, can when of ac fundamental frequency and electric system rated frequency inconsistent correctly calculate the fundamental voltage amplitude of of ac reality.
In order to solve the problem, the technical solution used in the present invention is:
The of ac fundamental voltage amplitude computing method used in Microcomputer Protection, comprise the following steps, the fundamental frequency f chosen 0when use Fourier series to calculate the fundamental voltage amplitude A (f, t) of of ac, t is sampled point; Adopt Frequency Measurement Algorithm to measure the frequency f of of ac, it is characterized in that: further comprising the steps of:
1), coefficient Coe (f) relevant to of ac frequency f is asked for;
1.1) ratio of the fundamental voltage amplitude that calculates of Fourier series and actual fundamental voltage amplitude is in theory:
A b ( f , t ) = | 2 f 0 sin ( πf / f 0 ) / π / ( f 0 2 - f 2 ) | ( f 0 2 + f 2 ) / 2 + ( f 0 2 - f 2 ) / 2 * cos ( 4 πft + 4 πfa ) - - - ( 5 )
Wherein f, t are respectively the frequency of of ac, the time, and 4 π fa+4 π ft are the twice of of ac x (t) at t phase angle, f 0for the fundamental frequency chosen when calculating Fourier series,
1.2) A is asked for b(f, t) mean value A within the 1/2f time bmean(f),
A bmean ( f ) = 2 f ∫ t t + 1 / 2 f A b ( f , t ) - - - ( 6 )
1.3) relevant to of ac frequency f coefficient Coe (f) is A bmeanthe inverse of (f),
Coe ( f ) = 1 / A bmean ( f ) - - - ( 7 ) ;
2) the fundamental voltage amplitude A (f, t), calculated Fourier series averages A within the 1/2f time mean(f):
A mean ( f ) = 2 f ∫ t t + 1 / 2 f A ( f , t ) - - - ( 8 ) ;
3), step 2) in the A that obtains mean(f) and step 1.3) in the Coe (f) that obtains be multiplied and namely obtain the fundamental voltage amplitude M of of ac reality:
M=A mean(f)/A bmean(f)=A mean(f)*Coe(f)(9)。
The invention has the beneficial effects as follows: by computing method of the present invention, can when of ac fundamental frequency and electric system rated frequency inconsistent correctly calculate the fundamental voltage amplitude of of ac reality.
Embodiment
In microcomputer protective relay device, data acquisition and processing (DAP) is all discrete, and each above-mentioned computation process all uses discrete method to process in protective relaying device.Frequency f and the fundamental voltage amplitude A (f, n) (n is current sampling point) of of ac has been tried to achieve in Microcomputer Protection.Therefore as follows to the computing method of of ac fundamental voltage amplitude in Microcomputer Protection:
1. suppose that the scope of the of ac fundamental frequency processed in Microcomputer Protection is between 5Hz ~ 70Hz; be that interval carry out discretize to frequency with 0.1Hz to this frequency range of 5Hz ~ 70Hz, ask for the Coe (f) at various discrete frequency place in advance and exist inside a data form.I-th data inside data form are the Coe (f) that frequency (5+0.1i) Hz is corresponding, and the computing method for the Coe (f) at each discrete frequency f place are as follows:
1.1) ratio of the fundamental voltage amplitude that calculates of Fourier series and actual fundamental voltage amplitude is in theory:
A b ( f , t ) = | 2 f 0 sin ( πf / f 0 ) / π / ( f 0 2 - f 2 ) | ( f 0 2 + f 2 ) / 2 + ( f 0 2 - f 2 ) / 2 * cos ( 4 πft + 4 πfa ) - - - ( 10 )
Wherein f, t, 4 π fa are respectively the frequency of of ac, time, the initial phase angle of twice of ac, f 0for the fundamental frequency chosen when calculating Fourier series.
1.2) A is asked for b(f, t) mean value A within the 1/2f time bmean(f).
A bmean ( f ) = 2 f ∫ t t + 1 / 2 f A b ( f , t ) - - - ( 11 )
1.3) relevant to of ac frequency f coefficient Coe (f) is A bmeanthe inverse of (f).
Coe(f)=1/A bmean(f) (12)
2. ask for the mean value A of A (f, n) within the 1/2f time mean(f)
A mean ( f ) = 1 N 0 Σ k = n - N 0 + 1 n A ( f , k ) - - - ( 13 )
N in formula 0for the sampling number of Microcomputer Protection within the 1/2f time, n is current sampling point.
3. for the frequency f of of ac, solving equation (5+0.1i)=f obtains i=(f-5) * 10, carries out round obtain i to i 0, in the form calculated in the 1st step i-th 0individual data are Coe (f) corresponding to frequency f.
4. step 2) in the A that obtains mean(f) and step 3) in the Coe (f) that inquires be multiplied and namely obtain the fundamental voltage amplitude M of of ac reality.
More than show and describe ultimate principle of the present invention, principal character and advantage.The technician of the industry should understand; the present invention is not restricted to the described embodiments; what describe in above-described embodiment and instructions just illustrates principle of the present invention; without departing from the spirit and scope of the present invention; the present invention also has various changes and modifications, and these changes and improvements all fall in the claimed scope of the invention.Application claims protection domain is defined by appending claims and equivalent thereof.

Claims (1)

1. the of ac fundamental voltage amplitude computing method used in Microcomputer Protection, comprise the following steps, the fundamental frequency f chosen 0when use Fourier series to calculate the fundamental voltage amplitude A (f, t) of of ac, t is sampled point; Adopt Frequency Measurement Algorithm to measure the frequency f of of ac, it is characterized in that: further comprising the steps of:
1), coefficient Coe (f) relevant to of ac frequency f is asked for;
1.1) ratio of the fundamental voltage amplitude that calculates of Fourier series and actual fundamental voltage amplitude is in theory:
A b ( f , t ) = | 2 f 0 sin ( πf / f 0 ) / π / ( f 0 2 - f 2 ) | ( f 0 2 + f 2 ) / 2 + ( f 0 2 - f 2 ) / 2 * cos ( 4 πft + 4 πfa ) - - - ( 5 )
Wherein f, t are respectively the frequency of of ac, the time, and 4 π fa+4 π ft are the twice of of ac x (t) at t phase angle, f 0for the fundamental frequency chosen when calculating Fourier series,
1.2) A is asked for b(f, t) mean value A within the 1/2f time bmean(f),
A bmean ( f ) = 2 f ∫ t t + 1 / 2 f A b ( f , t ) - - - ( 6 )
1.3) relevant to of ac frequency f coefficient Coe (f) is A bmeanthe inverse of (f),
Coe(f)=1/A bmean(f) (7);
2) the fundamental voltage amplitude A (f, t), calculated Fourier series averages A within the 1/2f time mean(f):
A mean ( f ) = 2 f ∫ t t + 1 / 2 f A ( f , t ) - - - ( 8 )
3), step 2) in the A that obtains mean(f) and step 1.3) in the Coe (f) that obtains be multiplied and namely obtain the fundamental voltage amplitude M of of ac reality:
M=A mean(f)/A bmean(f)=A mean(f)*Coe(f) (9),
For the frequency f of of ac, solving equation (5+0.1i)=f obtains i=(f-5) * 10, carries out round obtain i to i 0, in the form that formula (7) calculates i-th 0individual data are Coe (f) corresponding to frequency f.
CN201210080890.2A 2012-03-23 2012-03-23 Method for computing alternating current component fundamental amplitude value in microcomputer protection Active CN102608422B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210080890.2A CN102608422B (en) 2012-03-23 2012-03-23 Method for computing alternating current component fundamental amplitude value in microcomputer protection

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210080890.2A CN102608422B (en) 2012-03-23 2012-03-23 Method for computing alternating current component fundamental amplitude value in microcomputer protection

Publications (2)

Publication Number Publication Date
CN102608422A CN102608422A (en) 2012-07-25
CN102608422B true CN102608422B (en) 2015-05-20

Family

ID=46525949

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210080890.2A Active CN102608422B (en) 2012-03-23 2012-03-23 Method for computing alternating current component fundamental amplitude value in microcomputer protection

Country Status (1)

Country Link
CN (1) CN102608422B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1314997A1 (en) * 2001-11-24 2003-05-28 EADS Deutschland Gmbh Method for HPRF-radar measurement
CN101561462A (en) * 2009-04-27 2009-10-21 国网电力科学研究院 Phasor amplitude measurement method of power system

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012029093A1 (en) * 2010-09-02 2012-03-08 トヨタ自動車株式会社 Frequency control apparatus and charging system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1314997A1 (en) * 2001-11-24 2003-05-28 EADS Deutschland Gmbh Method for HPRF-radar measurement
CN101561462A (en) * 2009-04-27 2009-10-21 国网电力科学研究院 Phasor amplitude measurement method of power system

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
一种用于频率偏移时有效值计算的修正方法;刘益青等;《电力系统自动化》;20080125;第32卷(第02期);正文第80-81页第1-2节 *
消除傅氏算法中衰减非周期分量的方法在信号频率偏离额定值时的特性研究;张艳霞等;《继电器》;20030415;第31卷(第04期);第15-17页 *
电网基波幅值检测方法改进;杨保海等;《九江学院学报(自然科学版)》;20051231(第03期);第7-10,110页 *

Also Published As

Publication number Publication date
CN102608422A (en) 2012-07-25

Similar Documents

Publication Publication Date Title
CN104049144B (en) A kind of implementation method of the synchronous phasor measurement for filtering attenuating dc component
CN102435844B (en) Sinusoidal signal phasor calculating method being independent of frequency
CN104020352B (en) A kind of synchronous phasor measuring method being applicable to M class PMU unit
CN102818921B (en) Method for calculating amplitude and phase of alternating current signal based on iterative Fourier transform
WO2015157989A1 (en) Synchronous phasor measurement method applicable to p-type phasor measurement unit (mpu)
CN107589299B (en) Electric power signal synchronous phasor measuring method based on multi-frequency the measures model
CN102288804B (en) Method for calculating resistive current of arrester based on orthogonal transformation
Köse et al. Interharmonics analysis of power signals with fundamental frequency deviation using Kalman filtering
CN106970264A (en) A kind of improvement phase difference correction method for considering mains frequency rate of change
CN102508022B (en) Method for detecting power grid frequency by using optimal multiplier Newton algorithm
CN105137180A (en) High precision harmonic wave analysis method based on six item cosine window four spectral line interpolation
CN104360156B (en) The measuring method of signal phase angle when a kind of power system frequency offsets
CN109490630A (en) A kind of dynamic phasor measurement method based on pencil of matrix
KR101547266B1 (en) Apparatus and method for correcting phase between voltage and current signal of power meter and the power meter using the same
Jiekang et al. High-accuracy, wide-range frequency estimation methods for power system signals under nonsinusoidal conditions
CN105486921A (en) Kaiser third-order mutual convolution window triple-spectrum-line interpolation harmonic wave and inter-harmonic wave detection method
CN104090143B (en) Substation bus bar voltage magnitude instantaneous value measuring method
CN103543331A (en) Method for calculating harmonics and inter-harmonics of electric signal
CN103605904B (en) Self compensation power system amplitude arithmetic based on error estimation
CN104101781B (en) Substation bus bar voltage phase angle instantaneous value measuring method
CN102608422B (en) Method for computing alternating current component fundamental amplitude value in microcomputer protection
CN111830348B (en) Dynamic synchronous phasor measurement method and device under low-frequency oscillation
CN102012456A (en) Reactive power measurement method
CN104407197A (en) Signal phasor measurement method based on trigonometric function iteration
EP3199956B1 (en) Apparatus for determination of the frequency of an electrical signal and associated method

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

Address after: High road high tech Development Zone Nanjing city Jiangsu province 210061 No. 20

Co-patentee after: NARI Group Corp.

Patentee after: NARI TECHNOLOGY Co.,Ltd.

Address before: High road high tech Development Zone Nanjing city Jiangsu province 210061 No. 20

Co-patentee before: NARI Group CORPORATION STATE GRID ELECTRIC POWER INSTITUTE

Patentee before: NARI TECHNOLOGY Co.,Ltd.

CP01 Change in the name or title of a patent holder