Analysis of Small Oscillations in Complex Electric Power Systems

In this article the mathematical model of complex regulated electric system in matrix form is developed. This mathematical model makes it possible to study the steady-state stability of a complex electrical system by determining the eigenvalues of the dynamics matrix. The model of an electrical system that reflects transient processes for small deviations is convenient, both algorithmically and computationally, in particular, in cases of their joint solution with steady-state equations—the equations of nodal voltages. The obtained results in the form of the eigenvalues of the matrix spectrum are qualitatively the same as the results of classical studies, which is a consequence of the adequacy of the proposed model and the correct reflection of the dynamic processes occurring in a real electrical system. In addition, the equations obtained are of independent importance for the analysis of various modes, including transient, electrical systems of any complexity.


Introduction
The present stage in the development of the power industry is characterized by the presence of large concentrated energy systems connected by relatively weak connections, in which the powers of distributed generation are actively included.
The change in the composition of generation and the structure of power consumption leads to a decrease in the permanent inertia of the elements of the power systems, increasing the sensitivity of the parameters of the regime of the power system as a whole to small perturbations.
As is known [1] [2] [3] [4], the dynamic properties of complex electrical systems can differ significantly from those of simple electric power systems (EPS), which is confirmed by numerous field and model experiments and computational and experimental studies.In a multi-machine electrical system, the choice of the parameters of the control devices is much more complicated than in the simplest EPS.Therefore, as a rule, in the case of a multi-machine EPS, one generator or one station is considered to be adjustable and the parameters of their automatic excitation controllers (AECs) are determined based on the task at hand-ensuring equal damping, the required stability factor, etc., and the parameters of AECs of other stations are assumed to be given, with constant emf.for a certain inductive resistance [2].
In this paper we study the dynamic properties of electrical systems for small deviations (steady-state stability), described by linearized differential equations with constant coefficients.
Verification of the stability of power systems consists in determining the possibility of the existence of a stable regime with small perturbations of the parameters of the regime with given values of the parameters of the power system, the mode of generating sources, the load of node points, and the tuning of automatic mode control devices [3].
The complication of modern electrical systems, the introduction of digital and logical control devices into their structure requires refined and in-depth studies of the modes of electrical systems.Such a problem can be successfully solved by matrix methods.The article suggests a matrix model of the electric system, resolved based on the absolute angles of the generators, which emphasizes the relevance of the task and the method for solving it [4] [5] [6].
In the article the mathematical model and equations of multi-machine electric system, resolved concerning absolute angles of load of generators are received.
On the basis of the obtained model, the results of steady-state stability analysis will be obtained using the example of a three-generator electric system.Matrix equations of the elements of the EPS and the whole system were compiled on the basis of the most widely obtained equations of state variables [4] [5] [6], which are small deviations of the mode parameters-the angles of the rotor load of the synchronous generator, busbar voltages, power and other operating parameters of the EPS.The considered matrix equations are used for analysis of transient processes and steady-state stability of EPS and for the synthesis of optimal parameters of regulators of synchronous machines operating in an electrical system.

Mathematical Model of Transients in a Complex Electrical System
This model describes the transient process in the electric system, taking into account the balance of the moments (powers) on the shaft of the i-th aggregate of the electric power system and has the form [5]: where ω 0 is the synchronous angular frequency; , , , ji i Ti Gi T P P δ -is the inertia constant of the i-th aggregate, the load angle of the i-th generator, the mechanical power of the i-th turbine, the electromagnetic power of the i-th synchronous generator, respectively.
The equation of electromagnetic power of the i-th synchronous generator in the positional idealization has the form [6]: ( ) where E i , E j -emf.i-th and j-th synchronous generators; y ii , y ij -intrinsic and mutual conductivity of the network; α ii , α ij are complementary angles.
) ( ) where It should be noted that the derivation of formula ( 4) uses the obvious relationships: ( ) ( ) ( ) sin and cos 1, valid for small deviations in the load angles of generators.
After transformations (2), taking into account (3), (4), Equation (1) takes the form: ( ) and taking into account the parameters of the initial regime and the relation finally leads to a differential equation in the deviations: ( ) where In the case of the damper contours of the rotor of the i-th synchronous generator, Equation ( 6) takes the form: A. K. Raximovich, M. T. Farkhadovich DOI: 10.4236/eng.2018.105017256 Engineering where P di is the coefficient of the generalized damper moment of the i-th generator.
If the deviation of the emf is taken into account.i-th synchronous generator, Equation ( 7) takes the form [7]: The peculiarity of Equation ( 8) is that it is allowed with respect to the absolute angles of the system generators and, for example, for the three-generator electric system has the form [6] [7] [8]: The equations of electromagnetic transient processes in the excitation circuit of the i-th synchronous machine in the deviations have the form [1] [2] [4]: ( ) ( ) where , , di ei pi T T T -the transition time constant of the excitation winding, the exciter time constant, the automatic excitation controller, respectively; , , -deviations of the synchronous, forced emf.and the voltage at the output of the automatic excitation controller, respectively.The formation of signals via the AEC ∆е i channels in an idealized form (provided that the constant times of the differentiating elements of the AEC are considered to be zero) can be represented in the form [9] [10]: ( ) ( ) where k 0Pk , k 1Pk , k 2Pk are the gain factors of the AEC on the deflection channels, the first and second derivatives of the regime parameters ΔP k , respectively, k is the number of adjustable mode parameters.
The advantage of Equations ( 7) and ( 8) is their dependence on the deviations of the absolute load angles of the generators (∆δ i ) rather than the relative angles (∆δ ij ), which provides computational convenience, since these equations can be joined to the equations of node voltages, whose solutions give absolute angles [11] [12] [13] [14] [15].
For small perturbations of the parameters of the regime, after the corresponding transformations (1) -( 13), it is possible to obtain a generalized block matrix А Σ of size (4n × 4n) for the dynamics of an electrical system with n generators having automatic strong excitation regulators in the form: The components of the matrix A Σ are defined in [6].
In this case, the vector-column of the state parameters containing the parameters of the electric system mode has the form: .
For example, for a three-generator EPS (Figure 1), assuming that the automatic excitation regulators react to voltage and load angle deviation of the generators (∆δ i , ∆U Gi ), as well as their first derivatives ( , ).The equation of the output of the automatic excitation controller for the i-th generator is: In this case, the matrix AΣ takes the form [6]: Vector-column of the space of states of parameters of the EPS regime: .
As can be seen, the generalized matrix A3 of the dynamics of the electrical system, consisting of 3 generators, is formed from the parameters of the system Engineering mode and the automatic regulation of the excitation of machines, and therefore fully characterizes the transient processes in this EPS.The matrix A3 is rather sparse, which is typical for a complex system containing n generators, so this fact determines the computational advantages of the proposed mathematical model in the calculation and experimental studies of EPS.

Example
As an example, consider the matrix (17) of the intrinsic dynamics of the three-generator electric system A3 (Figure 1).Table 1 shows the technical parameters of the electrical system under consideration.
Where X d -synchronous inductive resistance of the generator along the longitudinal axis; X c -resistance of branches (inductive resistance of power lines).
The result of the calculation is shown below.
The spectrum of the matrix of the intrinsic dynamics of the three-generator EPS A3 at the selected regime and system parameters is equal to: −0.0012 ±  Similar studies of the stability of complex electrical systems were carried out in [2] [5] [8], but they used mathematical models derived from the relative angles of the load of synchronous generators.In this study, the results were obtained with respect to the absolute angles of loads, which allows us to further determine and investigate a particular generator, the first approaching the stability limit.
Figure 2 shows a 3D visualization of the matrix a pseudo-spectrum obtained with the MATLAB software module EigTool.Note that the EigTool module is a software product developed by Oxford University [15].Its methodological basis is the method of Arnoldi computation of subspaces of A.N. Krylov [16].The horizontal axes in Figure 2 correspond to the axes of the complex plane.Logarithm of the norm of the resolvent function was postponed along the vertical axis.Peaks localize the eigenvalues of the matrix.

Conclusion
The dynamic properties of complex electrical systems can differ significantly from the properties of simple EPS, which is confirmed by numerous full-scale and model experiments and computational and experimental studies [12] [13]

Figure 1 .
Figure 1.Diagram of a three-generator electrical system.

Figure 2 .
Figure 2. 3D-visualization of the spectrum of a three-generator electrical system with a Hurwitz matrix A 3 .

Table 1 .
Parameters of elements of the three-generator power system.
has positive poles.The inclusion of AECs of other generators by deviation and the first derivative of the angle and voltage of the generators will ensure the stability of the system under study.A. K. Raximovich, M. T. Farkhadovich DOI: 10.4236/eng.2018.105017259 Engineering