CN1402013A - Method for pure electrical measuring no-load potential phasor of synchronous electric generator - Google Patents
Method for pure electrical measuring no-load potential phasor of synchronous electric generator Download PDFInfo
- Publication number
- CN1402013A CN1402013A CN02130955A CN02130955A CN1402013A CN 1402013 A CN1402013 A CN 1402013A CN 02130955 A CN02130955 A CN 02130955A CN 02130955 A CN02130955 A CN 02130955A CN 1402013 A CN1402013 A CN 1402013A
- Authority
- CN
- China
- Prior art keywords
- value
- sin
- cos
- represent
- ref
- 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
- 238000000034 method Methods 0.000 title claims abstract description 69
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 50
- 238000005070 sampling Methods 0.000 claims abstract description 35
- 230000001052 transient effect Effects 0.000 claims description 51
- XOFYZVNMUHMLCC-ZPOLXVRWSA-N prednisone Chemical compound O=C1C=C[C@]2(C)[C@H]3C(=O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 XOFYZVNMUHMLCC-ZPOLXVRWSA-N 0.000 claims description 48
- 230000008569 process Effects 0.000 claims description 38
- 238000005259 measurement Methods 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 16
- 238000000691 measurement method Methods 0.000 claims description 15
- 238000004804 winding Methods 0.000 claims description 15
- 238000013016 damping Methods 0.000 claims description 12
- 239000013598 vector Substances 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 239000013604 expression vector Substances 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 claims description 3
- ODKSFYDXXFIFQN-UHFFFAOYSA-M argininate Chemical compound [O-]C(=O)C(N)CCCNC(N)=N ODKSFYDXXFIFQN-UHFFFAOYSA-M 0.000 claims description 3
- 229910052802 copper Inorganic materials 0.000 claims description 3
- 239000010949 copper Substances 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 230000009897 systematic effect Effects 0.000 claims description 3
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims description 2
- 230000003068 static effect Effects 0.000 claims description 2
- 230000008901 benefit Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 8
- 230000000875 corresponding effect Effects 0.000 description 7
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000005622 photoelectricity Effects 0.000 description 2
- 238000013341 scale-up Methods 0.000 description 2
- 230000003442 weekly effect Effects 0.000 description 2
- 241000143229 Idaea consanguinaria Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000005674 electromagnetic induction Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Landscapes
- Control Of Eletrric Generators (AREA)
Abstract
A electric method for measuring the idling potential of synchronous electric generator includes sampling output phase voltage and phase current at intervals, finding the real and imaginary parts of the positive-order components for voltage and current, calculating active voltage value and output electromagnetic power, prescribing dynamic procedure of synchronous electric generator, creating dynamic model, and calculating idling potential. Its advantage is high precision.
Description
Technical field
A kind of pure electric measurement method of synchronous generator no-load emf phasor belongs to the Power System and its Automation technical field.
Background technology
No-load emf E
qAnd phasing degree (being also referred to as rotor electrical angle, merit angle) δ (below be referred to as no-load emf phasor E
q∠ δ) is the important parameter that characterizes the synchronous generator running status and differentiate stability of power system.For many years, numerous research workers have explored measuring of the method that can accurately measure it, particularly merit angle δ always and have paid attention to widely and deep research.In recent years, along with the theoretical research of phasor measurement unit (PMU, Phasor Measurement Unit) and carrying out in a deep going way of application, the measurement of no-load emf phasor is had more importance and urgency.
Existing measuring method is divided into two big classes from principle: class methods realize measuring by non electrical quantity sensor (comprising photoelectricity or magnetoelectricity conversion).The common practice is, mechanical measuring point or measure speed gears are set on armature spindle, photoelectricity or calutron are installed on peritrochanteric, the latter receives the pulse signal that is produced by the former or utilizes electromagnetic induction principle to obtain the amount of certain and rotor-position or velocity correlation, and then realizes the measurement at merit angle by certain conversion.These class methods often need the generator body is carried out in various degree transformation, complex process, and, need signal Processing and Error Compensation Technology by means of more complicated owing to adopt the non electrical quantity sensor, to remove various structural errors, from measurement data, to obtain required characteristic quantity; And, be difficult to be applicable to other generator at the case proposition, cause realizing that cost is bigger.Another kind is pure electric measurement method, and the output voltage, electric current of promptly gathering synchronous generator be or/and other electric parameters, and then by theoretical analysis with calculate and obtain the amplitude of no-load emf or/and the merit angle.The advantage of this method is not contain non-electric sensor, universalization degree height, realize that cost is low.The realization technology that a kind of typical case of these class methods also is simultaneously simplification is based on the analytical Calculation method of motor steady-state operation equation or phasor graph, when it moves at systematic steady state, can calculate the no-load emf phasor more exactly, but in system's transient state process, because the analytic formula that method relied on is untenable, causes the bigger error of calculation.As, according to Simulation results, in the transient state process after generator end generation three phase short circuit fault, its measuring error can reach 50%.
Summary of the invention
The object of the present invention is to provide a kind of pure electric measurement method of synchronous generator no-load emf phasor, this method has been taken into account the dynamic behaviour characteristic of synchronous generator, thereby all has higher measuring accuracy in the stable state of motor and transient state (comprising time transient state) operational process.
The present invention is a kind of pure electric measurement method of synchronous generator no-load emf phasor, comprises the Park equation model of describing the synchronous generator dynamic process, it is characterized in that the method includes the steps of:
At first carry out initialization of variable:
N: AC sampling is counted, the integer greater than 12;
Δ T: in the AC sampling cycle, initial value is the 20/N millisecond, can adopt the sampling of fixed sampling interval technique or changing distance;
The initial zero clearing of each sampling buffer;
With Δ T is the value of gathering output phase voltage and output current phase at interval, converts them to per unit value, when sampling number is less than N, then continues to wait for new sampled point; When sampling number is more than or equal to N, whenever obtain a new sampled point, it and preceding N-1 sampled point form a new sampling period;
To the N point data in each new sampling period, order is carried out following steps;
1) adopt the positive sequence conversion to ask for the positive-sequence component real part and the imaginary part of voltage, electric current; And the sampled data of utilizing the new sampling period is used integral mean method calculating voltage effective value U
tWith output electromagnetic power P
e
2) adopt the Park equation model to describe the dynamic process of synchronous generator, establish on its rotor d-axis (also claiming the d axle) and comprise field copper and damping winding, hand over axle (also claiming the q axle) to comprise two damping winding, the q axle is leading d axle 90 degree on electric phase place; Ignore the influence of stator winding transient state process, rotation speed change to the influence of speed electromotive force and damping winding electric current to no-load emf amplitude E
qWith d axle transient potential increment Delta E '
dInfluence, set up the practical dynamic model of measuring, this model comprises:
A. Δ E '
dDynamic process adopt following transport function to describe, wherein input is I
q, output is Δ E '
d:
X
q: the reactance of q axle, adopt per unit value to represent;
X '
q: q axle subtranient reactance, adopt per unit value to represent;
T
Q0: q axle open circuit time constant, unit are second;
I
q: rotor q shaft current, adopt per unit value to represent;
B.E
q, E '
qRepresent by following simultaneous equations with the relation of θ:
E
qsinθ=E′
qsinθ+(X
d-X′
d)I
dsinθ,
E
qcosθ=E′
qcosθ+(X
d-X′
d)I
dcosθ,
E′
qsinθ=ΔE′
dcosθ-A+(X′
d-X
q)sinθ[cosθI
r+sinθI
i],
E′
qcosθ=-ΔE′
dsinθ+B+(X′
d-X
q)cosθ[cosθI
r+sinθI
i],
Parameter A wherein, B is defined as follows:
A=U
r+RI
r-X
qI
i,
B=U
i+X
qI
r+RI
i;
R: stator winding resistance, adopt per unit value to represent;
U
r: the real part of terminal voltage positive-sequence component, adopt per unit value to represent;
U
i: the imaginary part of terminal voltage positive-sequence component, adopt per unit value to represent;
I
r: the real part of output current positive-sequence component, adopt per unit value to represent;
I
i: the imaginary part of output current positive-sequence component, adopt per unit value to represent;
E '
q: q axle transient potential, adopt per unit value to represent;
ω: angular frequency, adopt per unit value to represent;
ω
0: specified angular frequency, for the 50Hz system, constant is 100 π radian per seconds;
θ: the absolute value of the electric anglec of rotation of no-load emf vector, unit is a radian, that is:
θ wherein
0Be corresponding t
0Initial value constantly;
I
d: rotor d shaft current, adopt per unit value to represent;
X
d: the reactance of d axle, adopt per unit value to represent;
X '
d: d axle subtranient reactance, adopt per unit value to represent.
3) measure the no-load emf phasor E that dynamic model calculates synchronous generator according to the practicality of synchronous generator
q∠ δ;
4) utilize and practically to measure no-load emf phasor that dynamic model calculates and in the transient state of synchronous generator and steady-state process, have excellent precision, but after system suffers big disturbance, last tens of to hundreds of milliseconds inferior transient state process, the measurement at merit angle has certain error, needs further to revise measured value.
The practicality according to synchronous generator in the described step 3) is measured dynamic model and is calculated the no-load emf phasor, comprises the steps:
A. calculating parameter A, B:
A=U
r+RI
r-X
qI
i,
B=U
i+X
qI
r+RI
i;
B. making the d axle transient potential increment of current time temporarily is the value of last sampled point, i.e. Δ E '
d(t)=Δ E '
d(t-Δ T), cycle counter Counter=0;
C. cancellation E ' from following simultaneous equations
q:
E′
qsinθ=ΔE′
dcosθ-A+(X′
d-X
q)sinθ[cosθI
r+sinθI
i],
E′
qcosθ=-ΔE′
dsinθ+B+(X′
d-X
q)cosθ[cosθI
r+sinθI
i],
Obtain equation Acos θ+Bsin θ=Δ E '
d
If its discriminant
Then make Δ E '
d=0, Counter is that maximum cycle adds 1, from equation Acos θ+Bsin θ=Δ E '
dIn solve sin θ and cos θ;
If its discriminant
Then can be from equation Acos θ+Bsin θ=Δ E '
dIn solve sin θ and cos θ; After obtaining sin θ and cos θ, calculate E ' according to aforementioned simultaneous equations
qSin θ and E '
qCos θ;
D. calculate q shaft current I according to the Parker transform
q(t);
E. make Δ E ' temporarily
d(t)
0=Δ E '
d(t), according to Δ E '
dThe dynamic process transport function, utilize numerical method to upgrade Δ E '
d(t), wherein input is I
q, output is Δ E '
d:
If Δ E '
d(t) with Δ E '
d(t) error between 0 exceeds maximum cycle less than limits of error η or Counter value, then continues next step f, otherwise, make Counter increase 1, go to step c;
F. according to Parker transformation calculations d shaft current I
d(t), and then according to following formula calculate E
qSin θ and E
qCos θ:
E
qsinθ=E′
qsinθ+(X
d-X′
d)I
dsinθ,E
qcosθ=E′
qcosθ+(X
d-X′
d)I
d?cosθ;
G. with E
qSin θ and E
qCos θ sequence multiply by modulates complex carrier signal
, and ask the mean value of its one-period, and obtain corresponding two vectors, be made as C
1+ jD
1, C
2+ jD
2Wherein the CALCULATION OF PARAMETERS formula is as follows:
T=2 π/ω wherein
Ref
In assumed condition | Δ ω |<<ω
RefThink E in [t-T, the t] period down,
qConstant with δ, and make t
0Constantly, θ
Ref0=0, then following relational expression is set up:
E
q=| C
1+ jD
1|, δ=arg (C
1+ jD
1), E
q=| C
2+ jD
2|, δ=arg (C
2+ jD
2)-90 °; Wherein || the mould of expression vector, the angle of arg () expression vector;
ω
Ref: system angle frequency reference value, adopt per unit value to represent, can get and make systematic steady state/specified angular frequency or center of inertia angular frequency;
Δ ω=ω-ω
Ret: angular frequency deviation, adopt per unit value to represent;
θ
Ref: the reference value of the electric anglec of rotation of no-load emf vector, unit is a radian, and
θ wherein
Ref0Be corresponding t
0Initial value constantly;
δ=θ-θ
Ref: the relative value of the electric anglec of rotation of no-load emf vector, i.e. merit angle, unit is a radian; Get E
q=(| C
1+ jD
1|+| C
2+ jD
2|)/2, δ=[arg (C
1+ jD
1)+arg (C
2+ jD
2)-90 °]/2, thereby the no-load emf phasor E of synchronous generator obtained
q∠ δ.
Utilizing practical measurement dynamic model to calculate in the process of no-load emf phasor, the limits of error η of employing is 10
-4, maximum cycle is 10.
In the pure electric measurement method of synchronous generator no-load emf phasor, the method for revising the measuring power angle value may further comprise the steps:
At first, order is δ (t) according to the merit angle that the practical measurement of synchronous generator dynamic model calculates
0
Secondly, wave equation according to generator, promptly following transport function equation with numerical calculations angular frequency deviation Δ ω, wherein is input as P
e, be output as Δ ω:
T wherein
wBe that a value is 2~10 seconds a time constant;
H: inertia constant, unit is second;
D: ratio of damping, no unit;
Once more, calculate a coefficient lambda that shows the current running status of system away from inferior transient state degree, the method of calculating λ is: λ → 0 expression system is in time transient state process, inferior transient state process after λ → 1 expression disturbance finishes substantially, the reflection running status that λ can be level and smooth between [0,1] is away from the degree of inferior transient state;
At last, revise the merit angle:
If λ (t)=1 then make merit angle nonce δ (t)=δ (t)
0, otherwise make δ ' (t)=δ ' (t-Δ T)+ω
0Δ ω Δ T;
And then calculate the merit angle according to following formula:
δ(t)=λδ(t)
0+(1-λ)δ′(t),
Obtain the no-load emf phasor E of synchronous generator time transient state process
q∠ δ.
The coefficient lambda computing method are as follows:
λ(t)=min{1,λ′(t)},
Wherein integer M ≈ 2max (T "
D0, T "
Q0)/Δ T;
U
T, ref: the reference value of terminal voltage, adopt per unit value to represent;
T "
Q0: q axle open circuit time transient state time constant, unit is second;
T "
D0: d axle open circuit time transient state time constant, unit is second;
Determine coefficient k
1, k
2Method be: when the variation value is in the static error scope that allows, promptly | Δ U
t(t) |<1% o'clock, λ ' was (t) in 0.95~1 scope; When the variation value exceeds the dynamic error scope of permission, promptly | Δ U
t(t) |>20% o'clock, λ ' was (t) in 0~0.01 scope.
Coefficient k
1Preferred value be 1.5, k
2Preferred value be 40.
Operate in software in the Industrial PC Computer and realized the pure electric measurement method of this synchronous generator no-load emf phasor, in steady-state process, at the parameter of electric machine accurately under the prerequisite, when sampling rate during for 48 in ripple weekly, the steady-state error that algorithm itself causes is no more than 2%; In the inferior transient state process in 300 milliseconds under the generator output rating after the serious disturbance of generation, when sampling rate during for 48 in ripple weekly, the Algorithm Error of measuring power angle is no more than 10%, and in 300 milliseconds of transient state processes in addition after disturbance, the Measurement Algorithm error is no more than 5%.As seen, the pure electric measurement method of this synchronous generator no-load emf phasor all has higher measuring accuracy in the stable state of motor and transient state (comprising time transient state) operational process
Description of drawings
Fig. 1 is for realizing a cover measurement mechanism synoptic diagram of the present invention.
Fig. 2 is the main program block diagram of the embodiment of the invention.
Fig. 3 is for calculating E '
qSin θ, E '
qThe subroutine block diagram of cos θ.
Fig. 4 is for upgrading Δ E '
d(t) subroutine block diagram.
Fig. 5 is for calculating E
qSubroutine block diagram with δ.
The subroutine block diagram of Fig. 6 for merit angle δ is revised.
Embodiment
The pure electric measurement method of the synchronous generator no-load emf phasor that the present invention puts forward can adopt the multiple hardwares scheme to realize, this example is introduced a kind of measuring system that realizes based on Industrial PC Computer more intuitively, comprises relevant hardware configuration and software flow.
The hardware configuration of measuring system comprises Industrial PC Computer as shown in Figure 1, ABC threephase potential transformer (PT) summation current transformer (CT), and from the A/D capture card on the Industrial PC Computer of being installed in of the secondary side of little PT and little CT output image data.
Synchronous generator outlet bus place, PT and CT be measuring machine end ABC three-phase voltage and electric current respectively, obtains ± the interior ac voltage signal of 120V scope, delivers to little PT and little CT, is converted to ± the interior voltage signal of 5V scope by it; The A/D capture card carries out low-pass filtering, AC sampling and digital-to-analog conversion with little PT and little CT secondary side voltage signal, and the digital signal that obtains is transferred to Industrial PC Computer CPU by pci bus, handle by corresponding software programs, finish the measurement of no-load emf phasor, and then output to user's monitoring interface.
Measuring method proposed by the invention mainly is embodied on the software that operates in the Industrial PC Computer, the master routine of software as shown in Figure 2,
Measuring method may further comprise the steps:
1. initialization
A. the parameter of pointing out the user must want from the computer interface input comprises the synchronous generator parameter: stator winding resistance R, q axle reactance X
q, d axle reactance X
d, q axle subtranient reactance X '
q, d axle subtranient reactance X '
d, q axle open circuit time constant T
Q0, q axle open circuit time transient state time constant T "
Q0, d axle open circuit time transient state time constant T "
D0, inertia constant H, ratio of damping D, the specified angular frequency of system
0, the system reference angular frequency
Ref(default value is ω
0); If input is famous value, then the user also must provide voltage base value U
BaseWith power base value S
Base, the famous value that Automatic Program will be correlated with changes into per unit value, and conversion formula is: per unit value=famous value/base value.The selection of base value is according to voltage base value U
BaseBe the specified phase voltage of generator, power base value S
BaseBe the generator rated capacity, electric current base value I
Base=P
Base/ (3U
Base), the impedance base value
The base value of angular frequency is ω
0The R=0.0181 Ω of system in this example,
S
Base=30/0.95MVA, the then base value of impedance
The per unit value of R is 0.0181/6.0306=0.003.
B. the sampling number N that points out the user to set every power frequency period (50Hz system), this routine N=48 represents every power frequency period (20 milliseconds) sampling 48 points, thus sampling interval Δ T=20/48 millisecond.
C. point out the relevant configured parameter of user's input measurement device, comprising: (this example is the PT no-load voltage ratio
, CT no-load voltage ratio (this example is 600A/120V), little PT no-load voltage ratio (this example is 120V/5V), little CT no-load voltage ratio (this example is 120V/5V), filtering scale-up factor (this example is 0.95) and A/D conversion coefficient (this example is 5V/3FFFh) etc.
D. parameter initialization in the algorithm: d axle transient potential increment initial value Δ E '
d(0)=0, calculates rotor q shaft current initial value I
q(0); Each sampling buffer zero clearing: u
a(t-k Δ T)=0, u
b(t-k Δ T)=0, u
c(t-k Δ T)=0, i
a(t-k Δ T)=0, i
b(t-k Δ T)=0, i
c(t-k Δ T)=0, k=0 ..., N-1,
They are used for circulating and deposit the synchronous generator output phase voltage u that collects
a, u
b, u
c, current i
a, i
b, i
cSignal forms the sampled point sequence, and wherein t is current sampling instant, and k is the numbering of sample sequence, and that k=0 represents is nearest, be t sampled point constantly.
2. under the driving of AC sampling data, carry out following link successively:
A. obtain new sampling number certificate from A/D sampling card, at first they are marked change, that is: consider the factors such as no-load voltage ratio, filtering scale-up factor and A/D conversion coefficient of PT/CT and little PT/CT, the digital quantity that obtains is converted to the actual physics value of synchronous generator outlet bus place correspondence, again according to formula: per unit value=famous value/base value obtains corresponding per unit value.
After obtaining the per unit value of voltage, electric current, with one of the value of voltage/current sample sequence pusher successively, nearest sampled point leaves on the data point of k=0 correspondence.
B. when sampling number is less than N, then continue to wait for new sampled point, otherwise, whenever obtaining a new sampled point, it and preceding N-1 sampled point form a new sampling period; To the N point data in each new sampling period, order is carried out following steps and is measured calculating.
1) calculates following parameter.
A. utilize the positive sequence conversion to ask for the real part and the imaginary part of the positive-sequence component of voltage, electric current, be respectively U
r, U
iAnd I
r, I
i, computing formula is:
B. utilize integral mean method calculating voltage effective value U by the sampled data of nearest one-period
tWith output electromagnetic power P
e, computing formula is:
k=0,1,...,N-1;
k=0,1,...,N-1;
2) establish on the rotor d axle of synchronous generator and comprise field copper and damping winding, the q axle comprises two damping winding, ignores the influence of stator winding transient state process, rotation speed change to the influence of speed electromotive force and damping winding electric current to no-load emf amplitude E
qWith d axle transient potential increment Delta E '
dInfluence, set up the practical dynamic model of measuring:
A. Δ E '
dDynamic process adopt following transport function to describe, wherein input is I
q, output is Δ E '
d:
B.E
q, E '
qRepresent by following simultaneous equations with the relation of θ:
E
qsinθ=E′
qsinθ+(X
d-X′
d)I
dsinθ,
E
qcosθ=E′
qcosθ+(X
d-X′
d)I
dcosθ, E′
qsinθ=ΔE′
dcosθ-A+(X′
d-X
q)sinθ[cosθI
r+sinθI
i],
E′
qcosθ=-ΔE′
dsinθ+B+(X′
d-X
q)cosθ[cosθI
r+sinθI
i],
Parameter A wherein, B is defined as follows:
A=U
r+RI
r-X
qI
i,
B=U
i+X
qI
r+RI
i;
E '
q: q axle transient potential, adopt per unit value to represent;
ω: angular frequency, adopt per unit value to represent;
θ: the absolute value of the electric anglec of rotation of no-load emf vector, unit is a radian, that is:
θ wherein
0Be corresponding t
0Initial value constantly;
3) measure the no-load emf phasor E that dynamic model calculates synchronous generator according to the practicality of synchronous generator
q∠ δ:
3.1) calculating parameter A, B
A=U
r+RI
r-X
qI
i,
B=U
i+X
qI
r+RI
i;
3.2) to make the d axle transient potential increment of current time temporarily be the value of last sampled point, i.e. Δ E '
d(t)=Δ E '
d(t-Δ T), counter Counter=0; Maximum cycle is 10 in this example.
3.3) cancellation E ' from following simultaneous equations
q:
E′
qsinθ=ΔE′
dcosθ-A+(X′
d-X
q)sinθ[cosθI
r+sinθI
i],
E′
qcosθ=-ΔE′
dsinθ+B+(X′
d-X
q)cosθ[cosθI
r+sinθI
i],
Obtain equation Acos θ+Bsin θ=Δ E '
d, therefrom solve sin θ and cos θ.
As shown in Figure 3,
Calculate sin θ and cos θ according to x:
sinθ=2x/(1+x
2),cosθ=(1-x
2)/(1+x
2)。
Note, when x=∞, sin θ=0, cos θ=-1.
And then according to above-mentioned simultaneous equations formula calculating E '
qSin θ, E '
qCos θ.
3.4) according to Parker transformation calculations q shaft current I
q(t), computing formula is:
3.5) make Δ E ' temporarily
d(t)
0=Δ E '
d(t), according to Δ E '
dThe dynamic process transport function, utilize numerical method to upgrade Δ E '
d(t), as shown in Figure 4, wherein input is I
q, output is Δ E '
d:
Promptly calculate according to following formula successively, obtain new Δ E '
d(t), z wherein, pz is a temporary variable.
z(t-ΔT)=ΔE′
d(t-ΔT)/(X
q-X′
q)+I
q(t-ΔT),
pz(t-ΔT)=1/T
q0[z(t-ΔT)-I
q(t-ΔT)],
z(t)=ΔE′
d(t)/(X
q-X′
q)+I
q(t),
pz(t)=1/T
q0[z(t)-I
q(t)],
z(t)=z(t-ΔT)+[pz(t-ΔT)+pz(t)]/2·ΔT,
Δ E '
d(t)=-(X
q-X '
q) [I
q(t)-z (t)]; If | Δ E '
d(t)-Δ E '
d(t) (η is the limits of error to 0|<η, and this example gets 10
-4) or Counter>10 item continue next step 3.6), otherwise, make Counter=Counter+1, go to 3.3);
3.6) according to Parker transformation calculations d shaft current I
d(t), computing formula is:
And then according to following formula calculating E
qSin θ and E
qCos θ:
E
qsinθ=E′
qsinθ+(X
d-X′
d)I
dsinθ,
E
qcosθ=E′
qcosθ+(X
d-X′
d)I
dcosθ。
3.7) calculating E
qAnd δ, its subroutine block diagram as shown in Figure 5:
A. with E
qSin θ and E
qCos θ sequence multiply by modulates complex carrier signal
, and ask the mean value of its one-period, and obtain corresponding two vectors, be made as C
1+ jD
1, C
2+ jD
2, wherein the CALCULATION OF PARAMETERS formula is as follows:
T=2 π/ω wherein
Ref
B. in multi-machine power system, because the relative value at merit angle is just meaningful, at its definition δ=θ-θ
RefIn, θ
RefTo choose mainly be to decide on convenient, it can be taken as certain with reference to the rotor electrical angle of unit or the equivalent phasing degree of system inertia center (COI), or even the integration of stable state (specified) angular frequency.Because
Thereby
。Make t
0Constantly, θ
Ref0=0, then
。In assumed condition | Δ ω |<<ω
RefThink E in [t-T, the t] period down,
qConstant with δ, thus there is following relational expression: C
1=E
qCos δ, D
1=E
qSin δ; C
2=-E
qSin δ, D
2=E
qCos δ; By C
1(C
2) be real part, D
1(D
2) for imaginary part constitutes vector, then have: E
q=| C
1+ jD
1|, δ=arg (C
1+ jD
1), E
q=| C
2+ jD
2|, δ=arg (C
2+ jD
2)-90 °; Wherein || the mould of expression vector, the angle of arg () expression vector; Get: E
q=(| C
1+ jD
1|+| C
2+ jD
2|)/2, δ=[arg (C
1+ jD
1)+arg (C
2+ jD
2)-90 °]/2;
Thereby obtain the no-load emf phasor E of synchronous generator
q∠ δ.
4) the no-load emf phasor E that obtains according to above-mentioned steps
q∠ δ has excellent precision in the transient state of synchronous generator and steady-state process; But after system suffers big disturbance, last tens ofly to hundreds of milliseconds inferior transient state processes, the measurement at merit angle has certain error, needs further correction measured value, the subroutine block diagram that merit angle δ is revised as shown in Figure 6:
A. making the merit angle that calculates according to the practical measurement of synchronous generator dynamic model is δ (t)
0
B. wave equation according to generator, promptly following transport function equation with numerical calculations angular frequency deviation Δ ω, wherein is input as P
e, be output as Δ ω:
T wherein
wBe that a value is 2~10 seconds a time constant;
C. calculate one and show that the current running status of system is made as λ away from the coefficient of inferior transient state degree.The method of calculating λ is: λ → 0 expression system is in time transient state process, inferior transient state process after λ → 1 expression disturbance finishes substantially, [0,1] λ can smoothly hang down the reflection running status away from inferior transient state degree between, can adopt parameters such as voltage, inferior transient state time constant in the λ computing method.The λ computing method that this example adopts are as follows:
k
1=1.5,k
2=40,
λ(t)=min{1,λ′(t)};
Wherein integer M ≈ 2max (T "
D0, T "
Q0)/Δ T;
D. revise the merit angle:
If λ (t)=1 would make merit angle nonce δ ' (t)=δ (t)
0, otherwise make δ ' (t)=δ ' (t-Δ T)+ω
0Δ ω Δ T; And then according to following formula correction merit angle:
δ(t)=λδ(t)
0+(1-λ)δ′(t),
Obtain the no-load emf phasor E of synchronous generator time transient state process
q∠ δ.
Above-mentioned measuring method all can obtain very high precision at inferior transient state, transient state and the steady-state process of synchronous generator, and can be on the monitor of Industrial PC Computer dynamic refresh measurement result E
q∠ δ.
Claims (6)
1. the pure electric measurement method of a synchronous generator no-load emf phasor comprises the Park equation model of describing the synchronous generator dynamic process, it is characterized in that the method includes the steps of:
At first carry out initialization of variable:
N: AC sampling is counted, the integer greater than 12;
Δ T: in the AC sampling cycle, initial value is the 20/N millisecond, can adopt the sampling of fixed sampling interval technique or changing distance;
The initial zero clearing of each sampling buffer;
With Δ T is the value of gathering output phase voltage and output current phase at interval, converts them to per unit value, when sampling number is less than N, then continues to wait for new sampled point; When sampling number is more than or equal to N, whenever obtain a new sampled point, it and preceding N-1 sampled point form a new sampling period;
To the N point data in each new sampling period, order is carried out following steps;
1) adopt the positive sequence conversion to ask for the positive-sequence component real part and the imaginary part of voltage, electric current; And the sampled data of utilizing the new sampling period is used integral mean method calculating voltage effective value U
tWith output electromagnetic power P
e
2) adopt the Park equation model to describe the dynamic process of synchronous generator, establish on its rotor d axle and comprise field copper and damping winding, the q axle comprises two damping winding, and the q axle is leading d axle 90 degree on electric phase place; Ignore the influence of stator winding transient state process, rotation speed change to the influence of speed electromotive force and damping winding electric current to no-load emf amplitude E
qWith d axle transient potential increment Delta E '
dInfluence, set up the practical dynamic model of measuring, this model comprises:
A. Δ E '
dDynamic process adopt following transport function to describe, wherein input is I
q, output is Δ E '
d:
X '
q: the reactance of q axle, adopt per unit value to represent; X '
q: q axle subtranient reactance, adopt per unit value to represent; T
Q0: q axle open circuit time constant, unit are second; I
q: rotor q shaft current, adopt per unit value to represent; B.E
q, E '
qRepresent by following simultaneous equations with the relation of θ: E
qSin θ=E '
qSin θ+(X
d-X
d) I
dSin θ, E
qCos θ=E '
qCos θ+(X
d-X '
d) I
dCos θ, E '
qSin θ=Δ E '
dCos θ-A+ (X '
d-X
q) sin θ [cos θ I
r+ sin θ I
i], E '
qCos θ=-Δ E '
dSin θ+B+ (X '
d-X
q) cos θ [cos θ I
r+ sin θ I
i], parameter A wherein, B is defined as follows: A=U
r+ RI
r-X
qI
i, B=U
i+ X
qI
r+ RI
iR: stator winding resistance, adopt per unit value to represent; U
r: the real part of terminal voltage positive-sequence component, adopt per unit value to represent;
U
i: the imaginary part of terminal voltage positive-sequence component, adopt per unit value to represent;
I
r: the real part of output current positive-sequence component, adopt per unit value to represent;
I
i: the imaginary part of output current positive-sequence component, adopt per unit value to represent;
E '
q: q axle transient potential, adopt per unit value to represent;
ω: angular frequency, adopt per unit value to represent;
ω
0: specified angular frequency, for the 50Hz system, constant is 100 π radian per seconds;
θ: the absolute value of the electric anglec of rotation of no-load emf vector, unit is a radian, that is:
θ wherein
0Be corresponding t
0Initial value constantly;
I
d: rotor d shaft current, adopt per unit value to represent;
X
d: the reactance of d axle, adopt per unit value to represent;
X '
d: d axle subtranient reactance, adopt per unit value to represent;
3) measure the no-load emf phasor E that dynamic model calculates synchronous generator according to the practicality of synchronous generator
q∠ δ;
4) utilize and practically to measure no-load emf phasor that dynamic model calculates and in the transient state of synchronous generator and steady-state process, have excellent precision, but after system suffers big disturbance, last tens of to hundreds of milliseconds inferior transient state process, the measurement at merit angle has certain error, needs further to revise measured value.
2. the pure electric measurement method of synchronous generator no-load emf phasor according to claim 1 is characterized in that dynamic model calculating no-load emf phasor is measured in the practicality according to synchronous generator in the described step 3), comprises the steps:
A. calculating parameter A, B:
A=U
r+RI
r-X
qI
i,
B=U
i+X
qI
r+RI
i;
B. making the d axle transient potential increment of current time temporarily is the value of last sampled point, i.e. Δ E '
d(t)=Δ E '
d(t-Δ T), cycle counter Counter=0;
C. cancellation Δ E ' from following simultaneous equations
q:
E′
qsinθ=ΔE′
dcosθ-A+(X′
d-X
q)sinθ[cosθI
r+sinθI
i],
E′
qcosθ=-ΔE′
dsinθ+B+(X′
d-X
q)cosθ[cosθI
r+sinθI
i],
Obtain equation Acos θ+Bsin θ=Δ E '
d
If its discriminant
Then make Δ E '
d=0, Counter is that maximum cycle adds 1, from equation Acos θ+Bsin θ=Δ E '
dIn solve sin θ and cos θ;
If its discriminant
Then can be from equation Acos θ+Bsin θ=Δ E '
dIn solve sin θ and cos θ;
After obtaining sin θ and cos θ, calculate E ' according to aforementioned simultaneous equations
qSin θ and E '
qCos θ;
D. calculate q shaft current I according to the Parker transform
q(t);
E. make Δ E ' temporarily
d(t)
0=Δ E '
d(t), according to Δ E '
dThe dynamic process transport function, utilize numerical method to upgrade Δ E '
d(t), wherein input is I
q, output is Δ E '
d:
If Δ E '
d(t) with Δ E '
d(t)
0Between error exceed maximum cycle less than limits of error η or Counter value, then continue next step f, otherwise, make Counter increase 1, go to step c;
F. according to Parker transformation calculations d shaft current I
d(t), and then according to following formula calculate E
qSin θ and E
qCos θ:
E
qsinθ=E′
qsinθ+(X
d-X′
d)I
dsinθ,E
qcosθ=E′
qcosθ+(X
d-X′
d)I
dcosθ;
G. with E
qSin θ and E
qCos θ sequence multiply by modulates complex carrier signal
, and ask the mean value of its one-period, and obtain corresponding two vectors, be made as C
1+ jD
1, C
2+ jD
2Wherein the CALCULATION OF PARAMETERS formula is as follows:
T=2 π/ω wherein
Ref
In assumed condition | Δ ω |<<ω
RefThink E in [t-T, the t] period down,
qConstant with δ, and make t
0Constantly, θ
Ref0=0, then following relational expression is set up:
E
q=|C
1+jD
1|,δ=arg(C
1+jD
1),E
q=|C
2+jD
2|,δ=arg(C
2+jD
2)-90°;
Wherein || the mould of expression vector, the angle of arg () expression vector;
ω
Ref: system angle frequency reference value, adopt per unit value to represent, can get and make systematic steady state/specified angular frequency or center of inertia angular frequency;
Δ ω=ω-ω
Ref: angular frequency deviation, adopt per unit value to represent;
θ
Ref: the reference value of the electric anglec of rotation of no-load emf vector, unit is a radian, and
, θ wherein
Ref0Be corresponding t
0Initial value constantly;
δ=θ-θ
Ref: the relative value of the electric anglec of rotation of no-load emf vector, i.e. merit angle, unit is a radian; Get E
q=(| C
1+ jD
1|+| C
2+ jD
2|)/2, δ=[arg (C
1+ jD
1)+arg (C
2+ jD
2)-90 °]/2, thereby the no-load emf phasor E of synchronous generator obtained
q∠ δ.
3. the pure electric measurement method of synchronous generator no-load emf phasor according to claim 2 is characterized in that described limits of error η is 10
-4, described maximum cycle is 10.
4. the pure electric measurement method of synchronous generator no-load emf phasor according to claim 1 is characterized in that may further comprise the steps the method for the described correction measuring power angle value in the step 4):
At first, order is δ (t) according to the merit angle that the practical measurement of synchronous generator dynamic model calculates
0
Secondly, wave equation according to generator, promptly following transport function equation with numerical calculations angular frequency deviation Δ ω, wherein is input as P
e, be output as Δ ω:
T wherein
wBe that a value is 2~10 seconds a time constant;
H: inertia constant, unit is second;
D: ratio of damping, no unit;
Once more, calculate a coefficient lambda that shows the current running status of system away from inferior transient state degree, the method of calculating λ is: λ → 0 expression system is in time transient state process, inferior transient state process after λ → 1 expression disturbance finishes substantially, the reflection running status that λ can be level and smooth between [0,1] is away from the degree of inferior transient state;
At last, revise the merit angle:
If λ (t)=1 would make merit angle nonce δ ' (t)=δ (t)
0, otherwise make δ ' (t)=δ ' (t-Δ T)+ω
0Δ ω Δ T;
And then calculate the merit angle according to following formula:
δ(t)=λδ(t)
0+(1-λ)δ′(t),
Obtain the no-load emf phasor E of synchronous generator time transient state process
q∠ δ.
5. the pure electric measurement method of synchronous generator no-load emf phasor according to claim 4 is characterized in that the coefficient lambda computing method are as follows:
λ(t)=min{1,λ′(t)},
Wherein integer M ≈ 2max (T "
D0, T "
Q0)/Δ T;
U
T, ref: the reference value of terminal voltage, adopt per unit value to represent;
T "
Q0: q axle open circuit time transient state time constant, unit is second;
T "
D0: d axle open circuit time transient state time constant, unit is second;
Determine coefficient k
1, k
2Method be: when the variation value is in the static error scope that allows, promptly | Δ U
t(t) |<1% o'clock, λ ' was (t) in 0.95~1 scope; When the variation value exceeds the dynamic error scope of permission, promptly | Δ U
t(t) |>20% o'clock, λ ' was (t) in 0~0.01 scope.
6. the pure electric measurement method of synchronous generator no-load emf phasor according to claim 5 is characterized in that coefficient k
1Preferred value be 1.5, k
2Preferred value be 40.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB021309558A CN1185496C (en) | 2002-09-23 | 2002-09-23 | Method for pure electrical measuring no-load potential phasor of synchronous electric generator |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB021309558A CN1185496C (en) | 2002-09-23 | 2002-09-23 | Method for pure electrical measuring no-load potential phasor of synchronous electric generator |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1402013A true CN1402013A (en) | 2003-03-12 |
CN1185496C CN1185496C (en) | 2005-01-19 |
Family
ID=4746525
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB021309558A Expired - Fee Related CN1185496C (en) | 2002-09-23 | 2002-09-23 | Method for pure electrical measuring no-load potential phasor of synchronous electric generator |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1185496C (en) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101614554B (en) * | 2009-07-08 | 2011-07-20 | 保定市三川电气有限责任公司 | Measuring device and measuring method for continuous physical quantity |
CN103412139A (en) * | 2013-07-22 | 2013-11-27 | 南京南瑞继保电气有限公司 | Synchronous generator rotating speed measuring method |
CN107677960A (en) * | 2017-09-21 | 2018-02-09 | 国网福建省电力有限公司 | The Parameter Identification Method of Synchronous Generator measured based on PMU |
CN108092581A (en) * | 2018-01-12 | 2018-05-29 | 广东美的暖通设备有限公司 | Phase-voltage measurement method, apparatus, air conditioner and the medium of three-phase inverting circuit |
CN109490622A (en) * | 2018-11-10 | 2019-03-19 | 王浩 | A kind of electricity data calibration technology based on tri-state information |
CN109782013A (en) * | 2018-12-27 | 2019-05-21 | 深圳柴火创客教育服务有限公司 | Motor speed measuring method, apparatus, readable storage medium storing program for executing and electrical equipment |
CN110988676A (en) * | 2019-11-25 | 2020-04-10 | 北京昊鹏智能技术有限公司 | Mechanical equipment fault diagnosis method and device and storage medium |
CN111537785A (en) * | 2020-05-25 | 2020-08-14 | 国网山西省电力公司营销服务中心 | Three-phase secondary circuit no-load on-line inspection tester |
-
2002
- 2002-09-23 CN CNB021309558A patent/CN1185496C/en not_active Expired - Fee Related
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101614554B (en) * | 2009-07-08 | 2011-07-20 | 保定市三川电气有限责任公司 | Measuring device and measuring method for continuous physical quantity |
CN103412139A (en) * | 2013-07-22 | 2013-11-27 | 南京南瑞继保电气有限公司 | Synchronous generator rotating speed measuring method |
CN107677960A (en) * | 2017-09-21 | 2018-02-09 | 国网福建省电力有限公司 | The Parameter Identification Method of Synchronous Generator measured based on PMU |
CN107677960B (en) * | 2017-09-21 | 2020-01-10 | 国网福建省电力有限公司 | Synchronous generator parameter identification method based on PMU measurement |
CN108092581A (en) * | 2018-01-12 | 2018-05-29 | 广东美的暖通设备有限公司 | Phase-voltage measurement method, apparatus, air conditioner and the medium of three-phase inverting circuit |
CN108092581B (en) * | 2018-01-12 | 2020-06-12 | 广东美的暖通设备有限公司 | Phase voltage detection method and device of three-phase inverter circuit, air conditioner and medium |
CN109490622A (en) * | 2018-11-10 | 2019-03-19 | 王浩 | A kind of electricity data calibration technology based on tri-state information |
CN109782013A (en) * | 2018-12-27 | 2019-05-21 | 深圳柴火创客教育服务有限公司 | Motor speed measuring method, apparatus, readable storage medium storing program for executing and electrical equipment |
CN109782013B (en) * | 2018-12-27 | 2022-07-05 | 深圳柴火创客教育服务有限公司 | Motor speed measuring method and device, readable storage medium and electric equipment |
CN110988676A (en) * | 2019-11-25 | 2020-04-10 | 北京昊鹏智能技术有限公司 | Mechanical equipment fault diagnosis method and device and storage medium |
CN110988676B (en) * | 2019-11-25 | 2022-02-22 | 北京昊鹏智能技术有限公司 | Mechanical equipment fault diagnosis method and device and storage medium |
CN111537785A (en) * | 2020-05-25 | 2020-08-14 | 国网山西省电力公司营销服务中心 | Three-phase secondary circuit no-load on-line inspection tester |
Also Published As
Publication number | Publication date |
---|---|
CN1185496C (en) | 2005-01-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1207679C (en) | High-accuracy synchronous phasor measuring method | |
CN1109390C (en) | Device and method for stabilizing electric power system | |
CN1211907C (en) | Induction motor driver and its parameter estimation method | |
CN1463484A (en) | Inverter control method and its device | |
CN1042081C (en) | system for controlling a brushless DC motor | |
CN1558253A (en) | Method and apparatus for measuring complex impedance of cells and batteries | |
CN1838518A (en) | Power inverter system and method of correcting supply voltage of the same | |
CN1180275C (en) | Method for measuring motor constant of induction motor | |
CN1171378C (en) | Device for detecting rotor position of brushless d.c. motor | |
CN1747306A (en) | Motor drive control apparatus, motor drive control method and program of the same | |
CN1375124A (en) | Initial magnetic pole estimating device for AC synchronous motor | |
CN1050467C (en) | Motor control apparatus equipped with control means for controlling rotational position of motor | |
CN1956317A (en) | Field weakening vector controller for permanent magnet synchronous motor and module | |
CN1667942A (en) | Rotor position presuming method and apparatus, motor control method, compressor and program | |
CN101043194A (en) | Vector control method and apparatus for induction motor | |
CN1874140A (en) | Speed control apparatus of vector controlled alternating current motor | |
CN1402013A (en) | Method for pure electrical measuring no-load potential phasor of synchronous electric generator | |
CN101039093A (en) | Vektorsteuerungsvorrichtung fur dauermagnetmotor | |
CN1230973C (en) | Safety speed monitoring system for three-phase AC transmission without sensor | |
CN1825752A (en) | Sliding moding structure direct torque servo-driver | |
CN1588793A (en) | Method or controlling permanent magnet synchronous motor-air conditioner compressor system without speed sensor | |
CN87101640A (en) | The device of stabilizing torsional-vibration of turbogenerator shaft | |
CN101042426A (en) | Methods and apparatus for frequency rectification | |
CN1578111A (en) | Controller for AC motor and ac motor system | |
CN1517666A (en) | Measuring device and method for electric machine |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20050119 Termination date: 20180923 |
|
CF01 | Termination of patent right due to non-payment of annual fee |