Detailed Description
The details of the invention will be more clearly understood in conjunction with the accompanying drawings and description of specific embodiments of the invention. The specific embodiments of the invention described herein are for purposes of illustration only and are not to be construed as limiting the invention in any way. Given the teachings of the present invention, one of ordinary skill in the related art will contemplate any possible modification based on the present invention, and such should be considered to be within the scope of the present invention.
As shown in fig. 1 to 5, the present invention provides a method for determining damping characteristics of a generator excitation system, the method comprising:
determining an oscillation frequency of an oscillating unit (fault unit) based on synchronous phasor measurement data (PMU data) of the oscillating unit;
Determining an excitation voltage phase of the excitation system based on the synchrophasor measurement data (PMU data), and determining a hysteresis characteristic of the excitation system based on the oscillation frequency;
And determining the damping moment phase of the excitation system (PSS) according to the excitation voltage phase and the hysteresis characteristic of the excitation system, and determining the damping characteristic of the excitation system of the generator according to the damping moment phase of the excitation system.
The judging method provided by the invention can automatically analyze and judge the damping characteristic of the excitation system in the low-frequency oscillation based on the data of the wide area synchronous measurement system (WAMS) so as to provide an important basis for fault positioning.
The judging method provided by the invention not only can reduce the influence of low data quality on damping characteristic analysis and greatly save manpower and material resources for accident analysis, but also can obviously shorten the accident analysis time, is beneficial to quickly determining the oscillation cause, and has important significance for quickly recovering power production and improving the safety and reliability of a power system.
In an alternative embodiment of the invention, the synchrophasor measurement data (PMU data) comprises at least the excitation voltage, the active power and the rotational speed of the oscillating group.
In an alternative embodiment of the present invention, after the synchrophasor measurement data (PMU data) is collected, the synchrophasor measurement data (PMU data) is cleaned, and then the cleaned data is used for subsequent calculation and analysis.
In an alternative example of this embodiment, the synchrophasor measurement data (PMU data) D (Δt) at time t 1~t2 (where t 1 is the initial time of the interception of the data and t 2 is the end time of the interception of the data) is taken as the analysis object (t 2-t1 =Δt), fourier analysis is performed on it, only the fundamental component data D 1 (Δt) is retained, and other harmonic components are removed.
In an alternative example of this embodiment, the specific calculation of the oscillation frequency of the oscillating group from synchrophasor measurement data (PMU data) is as follows:
The zero crossing times in the excitation voltage U f (Δt), the active power P e (Δt), and the rotation speed ω (Δt) are recorded as T ' Uf_1~t'Uf_n、t'Pe_1~t'Pe_n and T ' ω_1~t'ω_n (n=1.2.3.) respectively, and the oscillation period average value t=2 (T ' Pe_n-t'Pe_1)/(n-1) is obtained, thereby finally obtaining the oscillation frequency f 0 =1/T.
In an alternative embodiment of the invention, the synchrophasor measurement data comprises the active power of the oscillating unit, wherein the determining the excitation voltage phase of the excitation system based on the synchrophasor measurement data comprises determining the voltage phase of the excitation system according to the active power.
Specifically, determining the phase of the excitation system based on the active power includes determining the spatial phase of the PSS output signal U PSS in the Δδ - Δω coordinate system (the excitation voltage U f is approximately the same as the phase of the PSS output signal U PSS in the Δω - Δδ coordinate system) based on the power deviation ΔP e in the PMU data, i.e., ΔP e is the starting point of the Δδ - Δω coordinate system phase 0, and calculating the phase of the excitation voltage U f relative to ΔP e (Lead positive and lag negative) as shown in equation (1):
in an alternative embodiment of the invention, the synchrophasor measurement data comprises the rotational speed of the oscillating group, wherein the determining the excitation voltage phase of the excitation system based on the synchrophasor measurement data comprises determining the voltage phase of the excitation system according to the rotational speed.
Specifically, determining the phase of the excitation system based on the rotational speed includes determining the spatial phase of the PSS output signal U PSS in the ΔΔω - Δω coordinate system (the excitation voltage U f is approximately the same as the phase of the PSS output signal U PSS in the Δω - Δδ coordinate system) based on the rotational speed Δω in the PMU data, i.e., Δω is the Δδ - Δω coordinate system phase 0 start, and calculating the excitation voltage U f relative to Δω phase(Lead positive and lag negative) is as shown in formula (2):
In the present invention, the exciting voltage phase may also be obtained by other methods well known to those skilled in the art, and will not be described herein.
In an alternative embodiment of the invention, determining the hysteresis characteristic of the excitation system based on the oscillation frequency comprises obtaining a phase frequency function by data fitting according to the phase of the damping moment of the excitation system relative to the output signal of the excitation system, and determining the hysteresis characteristic of the excitation system by the phase frequency function at the oscillation frequency.
In an alternative example of this embodiment, the gratuitous compensation characteristic of the excitation system may be obtained by direct fitting based on field test results.
In another alternative example of this embodiment, the gratuitous compensation characteristic of the excitation system may be obtained from theoretical calculations.
In this embodiment, since the exciter system of the oscillating unit has no compensation characteristic, which is numerically equal to the phase (positive lead and negative lag) of the damping torque ΔT PSS of the exciter system (PSS) relative to the output signal U PSS of the exciter system (PSS), the phase-frequency function can be obtained by data fittingCalculating uncompensated characteristic phase of corresponding excitation system under current f 0
In an alternative embodiment of the invention, determining the damping torque phase of the exciter system based on the exciter voltage phase and the hysteresis characteristic of the exciter system comprises summing the exciter voltage phase and the uncompensated characteristic phase. Namely, the excitation voltage phase and the corresponding excitation system uncompensated characteristic phase under the current f 0 Adding the phases of the damping moment of the excitation system (PSS) in the rotating speed and the power angle (delta omega-delta) coordinate system of the generator
In an alternative example of this embodiment, the excitation voltage phase is calculated by the active power P e (Δt), thenI.e. ΔΓ PSS phase (lead positive and lag negative) relative to Δp e.
In another alternative example of this embodiment, the excitation voltage phase is calculated from the rotational speed ω (Δt), thenI.e. ΔΓ PSS phase with respect to Δω (lead positive and lag negative).
In an alternative embodiment of the present invention, if the damping torque phase of the excitation system (PSS) is in the first quadrant and the second quadrant in the Δω - Δδ coordinate system, it is considered that the PSS provides a positive damping torque, and the larger the projection of ΔΣ PSS in the Δω axis positive direction, the larger the PSS provides a positive damping torque;
If the PSS damping torque phase is in the three-quadrant and the four-quadrant in the delta omega-delta coordinate system, the PSS is considered to provide negative damping torque, and the larger the projection of delta T PSS in the delta omega axis negative direction is, the larger the PSS provides negative damping torque.
In the invention, the basic principle of the PSS (power system stabilizer) is that one or more signals of the rotation speed deviation delta omega, the frequency deviation delta f and the electric power deviation delta P e of a generator are used as the input of the control of an excitation system, and positive damping moment is generated after a phase correction link, so that the purpose of restraining the low-frequency oscillation of the system is realized. Taking PSS2B with wider application in a unit as an example, taking DeltaP e and Deltaω as input signals, generating an output signal U PSS to act on an AVR after a phase compensation link to generate an additional moment DeltaT pss, and after proper PSS gain and phase compensation are determined, enabling the projection of a resultant moment of DeltaT pss and DeltaT e2(ΔTe2 on a Deltaω axis to be positive, wherein the excitation system provides positive damping, and the basic principle of the action is shown in figure 2.
PSS adopts the existing Philips-Heffron (Phillips-Fulong) single machine infinite system model. As shown in FIG. 3, where ΔT M is mechanical torque, ΔT e1 and ΔT e2 are both electromagnetic torque, Δω and Δδ are generator rotational speed and power angle variation, respectively, ΔU PSS is PSS output variation, ΔU REF is terminal voltage reference variation, ΔE fd is excitation voltage variation, ΔE q 'is generator transient potential variation, T j is unit rotational inertia, ω 0 is generator rated rotational speed, T d0' is generator time constant, and D is mechanical damping.
When the disturbance on the system side or the prime motor side causes power oscillation, the disturbance of the generator end voltage induction system generates a variable quantity DeltaU t, the PSS acquires generator active and rotating speed signals to generate a PSS output signal DeltaU PSS, the PSS, deltaU t and DeltaU REF act through an excitation regulator (AVR) to generate an excitation voltage variable quantity DeltaE fd, then act with a transfer function K 3(1+K3Td0 'to generate a DeltaE q', and finally act through a gain K 2 to generate a phase characteristic of an electromagnetic torque DeltaT e2,ΔTe2, so that the excitation system has the effect of inhibiting or aggravating the power oscillation. Wherein:
In the Δδ - Δω coordinate system, as shown in fig. 4a and 4b, when Δt e2 projects Δt e2 'on the Δω axis as positive, the excitation system provides positive damping, i.e. suppressing the low frequency oscillation, and when Δt e2 projects Δt e2' on the Δω axis as negative, the excitation system provides negative damping, i.e. exacerbating the low frequency oscillation. The judgment of the damping characteristic of the excitation system can be realized by analyzing the space phase of the electromagnetic torque delta T e2 in the delta-delta omega coordinate system.
The invention also provides a judging device of the damping characteristic of the generator excitation system, which comprises:
the oscillation frequency calculation module is used for determining and calculating the oscillation frequency of the oscillating unit based on synchronous phasor measurement data of the oscillating unit of the generator;
The exciting voltage phase calculation module is used for determining the exciting voltage phase of the exciting system based on the synchronous phasor measurement data;
a hysteresis characteristic calculation module of the magnetic system, which determines hysteresis characteristics of the excitation system based on the vibration frequency;
And the damping moment phase analysis module is used for determining the damping moment phase of the excitation system according to the excitation voltage phase and the hysteresis characteristic of the excitation system and judging the damping characteristic of the excitation system of the generator according to the damping moment phase of the excitation system.
The invention also proposes a computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the decision method as described above when executing the computer program.
In an alternative embodiment of the invention, the computer device (analysis means) makes the determination according to the following procedure:
(1) And (3) data import, namely, exporting PMU data collected on exciting voltage, active power and rotating speed of the oscillating unit (fault unit), and reading the PMU data into computer equipment (analysis device).
(2) And (5) data cleaning. At time t 1~t2 (t 1 is the start of data interception, t 2 is the end time of data interception), PMU data D (Δt) is the analysis object (t 2-t1 =Δt), fourier analysis is performed on the PMU data D (Δt), only fundamental component data D 1 (Δt) is retained, and other harmonic components are removed.
(3) The oscillation frequency f 0 is calculated. The zero crossing times in the records U f(Δt)、Pe (Δt) and ω (Δt) are T ' Uf_1~t'Uf_n、t'Pe_1~t'Pe_n and T ' ω_1~t'ω_n respectively (n=1.2.3.), the oscillation period average value t=2 (T ' Pe_n-t'Pe_1)/(n-1) is calculated to obtain the oscillation frequency f 0 =1/T.
(4) The excitation voltage U f phase is determined.
Method 1 determining the spatial phase of U PSS in the Δδ - Δω coordinate system based on ΔP e in the PMU data, i.e., ΔP e is the starting point of Δδ - Δω coordinate system phase 0, calculating the phase of U f relative to ΔP e (Lead positive and lag negative) as shown in equation (1):
method 2 determining the spatial phase of U PSS in the Δδ - Δω coordinate system based on Δω in the PMU data, i.e., Δω is the Δδ - Δω coordinate system phase 0 start, calculating the phase of U f relative to Δω (Lead positive and lag negative) is as shown in formula (2):
wherein, the excitation voltage U f and the PSS output signal U PSS are approximately the same phase in the Δω - Δδ coordinate system.
(5) And calculating hysteresis characteristics of the excitation system. The uncompensated characteristic of the excitation system of the oscillating unit is obtained through on-site actual measurement or theoretical calculation, the value of the uncompensated characteristic is equal to the phase (positive lead and negative lag) of PSS damping moment DeltaT PSS relative to the PSS output signal U PSS, and then a phase-frequency function is obtained through data fittingCalculating uncompensated characteristic phase of corresponding excitation system under current f 0
(6) The PSS damping torque DeltaT PSS phase is determined. Adding the calculation results of (4) and (5) to obtain the phase of the PSS damping moment in the delta omega-delta coordinate systemIf the calculation of method 1 in (4) is adopted, thenI.e. ΔT PSS phase relative to ΔP e (positive lead and negative lag), if calculated by method 2 in (4)I.e. ΔΓ PSS phase with respect to Δω (lead positive and lag negative).
(7) PSS damping characteristics analysis. If vector ΔT PSS is in the first and second quadrants of the Δω - Δδ coordinate system, the PSS is considered to provide positive damping torque, and the larger the projection of ΔT PSS in the Δω axis positive direction, the larger the PSS provides positive damping torque, and if vector ΔT PSS is in the third and fourth quadrants of the Δω - Δδ coordinate system, the PSS is considered to provide negative damping torque, and the larger the projection of ΔT PSS in the Δω axis negative direction, the larger the PSS provides negative damping torque.
The invention also proposes a computer readable storage medium storing a computer program which, when executed by a processor, implements the method of any one of claims 1 to 8.
The specific implementation process of the method and the device for judging the damping characteristic of the generator excitation system provided by the invention is described in detail by combining an embodiment, so that the effectiveness of the method and the device for judging the damping characteristic of the generator excitation system is proved:
The low-frequency oscillation phenomenon occurs in the deep peak shaving and phase advancing operation of a certain power plant #4 unit, the PSS is withdrawn in the oscillation process, the power oscillation subsides, and a technician preliminarily judges that the low-frequency oscillation event is related to abnormal output of an excitation system (PSS). In the following, the damping characteristics provided by the PSS are analyzed by the rotational speed variation Δω and the excitation voltage U f in the PMU data, taking this low frequency oscillation event as an example.
Step 1, the #4 unit PMU data (as shown in fig. 6) is led into an analysis device, data cleaning is performed, and the active oscillation frequency f 0 =1.80 Hz is calculated.
Step 2, selecting method 2 to calculate the phase of U f. The phase of U f relative to Deltaomega is calculated based on the Deltaomega of the #4 unitSince the excitation voltage U f is approximately the same as the PSS output signal U PSS in phase in the Δω - Δδ coordinate system, it is known that the PSS output signal U PSS is in phase with the Δω axis.
Step 3, the uncompensated characteristic of the excitation system of the #4 unit is measured through field test, and as shown in table 1, a phase frequency function is obtained through data fittingCalculating the uncompensated characteristic phase of the excitation system corresponding to f 0 =1.80 Hz I.e., f 0 =1.80 Hz, the PSS damping torque Δt PSS phase lags the PSS output signal U PSS by 111.8 °.
No. 1#4 unit excitation system uncompensated characteristic
| f(Hz) |
0.1 |
0.2 |
0.3 |
0.4 |
0.5 |
0.6 |
0.7 |
0.8 |
0.9 |
1 |
| Φ(°) |
-70.6 |
-103.1 |
-100.6 |
-101.6 |
-97.7 |
-92.5 |
-91.7 |
-92.7 |
-88.7 |
-90.3 |
| f(Hz) |
1.1 |
1.2 |
1.3 |
1.4 |
1.5 |
1.6 |
1.7 |
1.8 |
1.9 |
2 |
| Φ(°) |
-93.9 |
-101.3 |
-106.5 |
-106.5 |
-108.0 |
-116.4 |
-117.1 |
-111.8 |
-108.2 |
-103.8 |
Step 4, adding the calculation results of the step 2 and the step 3 to obtain the phase of PSS damping moment of the #4 unit relative to delta omegaI.e. ΔΓ PSS lags Δω111.8 °, in the Δω - Δδ coordinate system in quadrant IV, as shown in fig. 7, from which it is determined that the #4 set PSS provides negative damping during active oscillation.
And 5, deducing that the PSS is abnormal in output in the process because the PSS provides negative damping in the active oscillation process, and carrying out fault investigation on an excitation system. After investigation, the parameter setting in the excitation system PSS is incorrect, so that the PSS output is abnormal, and further power oscillation is caused.
The method and the device for judging the damping characteristic of the generator excitation system are based on a Philips-Heffron single machine infinite system mathematical model, PMU data of a low-frequency oscillation accident unit are led into the excitation system damping characteristic judging device, after data cleaning, the damping characteristic of the excitation system is judged by analyzing the space phase of electromagnetic moment delta T PSS generated by PSS in a delta omega coordinate system, and therefore the analysis of the cause of the low-frequency oscillation accident is assisted.
The detailed explanation of the embodiments described above is only for the purpose of explaining the present invention so as to enable a better understanding of the present invention, but these descriptions should not be construed as limiting the present invention in any way, and in particular, the individual features described in the different embodiments may be arbitrarily combined with each other to constitute other embodiments, and these features should be understood as being applicable to any one embodiment, except as explicitly stated to the contrary, without being limited to only the described embodiment.