Next Article in Journal
Hysteresis Modeling of Piezoelectric Actuators Based on a T-S Fuzzy Model
Previous Article in Journal
A Self-Tuned Method for Impedance-Matching of Planar-Loop Resonators in Conformable Wearables
Previous Article in Special Issue
Task Scheduling in Cloud Computing Environment Using Advanced Phasmatodea Population Evolution Algorithms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of Improved Quasi-Affine Transformation Evolutionary Algorithm in Power System Stabilizer Optimization

School of Electronic, Electrical Engineering and Physics, Fujian University of Technology, Fuzhou 350118, China
*
Author to whom correspondence should be addressed.
Electronics 2022, 11(17), 2785; https://doi.org/10.3390/electronics11172785
Submission received: 27 May 2022 / Revised: 15 August 2022 / Accepted: 24 August 2022 / Published: 4 September 2022
(This article belongs to the Special Issue Machine Learning in the Industrial Internet of Things)

Abstract

:
This paper proposes a parameter coordination optimization design of a power system stabilizer (PSS) based on an improved quasi-affine transformation evolutionary (QUATRE) algorithm to suppress low-frequency oscillation and improve the dynamic stability of power systems. To begin, the simulated annealing (SA) algorithm randomly updates the globally optimal solution of each QUATRE iteration and matches the inferior solution with a certain probability to escape the local extreme point. This new algorithm is first applied to the power system. Since the damping ratio is one of the criteria with which to measure the dynamic stability of the power system, this paper sets the objective function according to the principle of maximization of the damping coefficient of the electromechanical mode, and uses SA-QUATRE to search a group of global optimal PSS parameter combinations to improve the safety factor of the system as much as possible. Finally, the method’s rationality and validity were validated by applying it to the simulation examples of the IEEE 4-machine 2-area system with different operation states. The comparison with the traditional optimization algorithm shows that the proposed method has more advantages for multi-machine PSS parameter coordination optimization, can restrain the low-frequency oscillation of the power system more effectively and can enhance the system’s stability.

1. Introduction

In recent years, the continuous operation of large generator sets and the large-scale use of high-gain and fast excitation systems have improved the transient level of the system, but also caused the problem of low frequency oscillation of the power system due to insufficient damping, which has increasingly become an important factor restricting the overall interconnection of the power grid [1,2,3]. A power system stabilizer (PSS) can not only offset the negative damping torque generated by the regulator, but also provide additional positive damping for the system, which is the most efficient and economical measure for suppressing the low-frequency oscillation of the power system [4,5]. Especially in the context of a large power grid, PSS has been widely used. The configuration of PSS, that is, the research theory of the installation site, has been perfected, and another key issue is the parameter coordination, which is also one of the hot topics discussed by scholars. It has been shown that the damping of local and interval oscillation modes can be optimized to a certain extent if the PSS can obtain appropriate parameters [6].
In control theory, the phase compensation approach can be used to design a power system stabilizer, but for multiple power system stabilizers in the system, the traditional phase compensation method usually carries out parameter setting sequentially, which will inevitably cause the interaction effect between each stabilizer [7]. To overcome this phenomenon, some studies have pointed out that the parameter tuning of PSS is equivalent to an optimal problem based on the eigenvalues of electromechanical modes, that is, to move the eigenvalues to the left half plane of the complex plane as much as possible. In recent years, using the global optimization ability of artificial intelligence algorithms to solve this phenomenon has provided a new idea for the coordination design of a controller and has become a hot research topic. Particle swarm optimization (PSO), a basic and powerful approach that many researchers have studied and applied to power systems, is a commonly used optimization algorithm [8]. In Ref. [9], PSS parameters were optimized based on a PSO algorithm. Although the optimization principle is simple, the optimized particles in the optimization iteration process are prone to precocity and easily fall into local optimum, which is also the common fault of most algorithms. For this reason, Ref. [10] proposed an improved PSO algorithm based on particle swarm theory. Although the performance of the algorithm has been improved to a certain extent, the particles do not completely tend to the global optimal value, but tend to a relatively good value. Ref. [11] proposed the gray wolf optimization algorithm to optimize PSS parameters, which has fewer algorithm parameters and faster convergence speed. However, in the later optimization period, with the decrease of population number, the optimization rate decreases, and at the same time, it is possible to fall into local optimal. Therefore, Ref. [12] proposed an improved gray wolf optimization algorithm, which overcomes the problem that it is easy to fall into local optimum and strengthens the robustness of the algorithm. However, the optimization algorithm search process tends to be complicated, so the reliability of the optimization result needs to be further improved.
In this paper, a method combining the quasi-affine transformation evolutionary algorithm [13] and the simulated annealing algorithm is proposed to coordinate the optimization design of stabilizer parameters of multi-machine power system. The main contribution of the paper lies in the following:
The proposed SA-QUATRE algorithm not only overcomes the defect of slow convergence of particle cluster optimization algorithm, it also has stronger cooperation and reduces time complexity, and accepts inferior solutions with a certain probability to avoid falling into local minimum in the process of searching.
To overcome the interaction between power system stabilizers, in this paper, the parameter coordination tuning problem of the power system stabilizer is transformed into the optimal problem based on the characteristic value of oscillation mode and damping coefficient, and the dynamic stability of the power system is improved.
The validity and feasibility of the proposed method are verified by simulation with a test function and the IEEE 4-machine 11-node system. It can be found that, compared with the existing multi-machine power system with PSS parameters, the multi-machine and multi-node power systems with the SA-QUATRE algorithm designed and optimized under stabilizer suppression has better system dynamic performance.
The rest of the paper is organized as follows: Section 2 introduces PSS structure and each link, as well as its transfer function, objective function and constraints of the parameter optimization of PSS for the multi-machine system. Section 3 introduces the proposed SA-QUATRE algorithm in detail and PSS parameter optimization based on it. Section 4 discusses the simulation example and the experiment results. Section 5 presents the final conclusion.

2. Parameter Optimization of PSS for Multi-Machine System

Voltage stabilizers suppress oscillations by adding a signal to the excitation system that produces a positive damping torque to offset the negative damping torque produced by the Voltage regulator [14]. The principle and excitation systems are shown in Figure 1, where, taking the rotor angle deviation Δ w of the generator as the input signal, the PSS and excitation system transfer function can be expressed as [15]:
U P S S , i = K s , i s T 5 , i 1 + s T 5 , i ( 1 + s T 1 , i 1 + s T 2 , i ) ( 1 + s T 3 , i 1 + s T 4 , i ) Δ w ,
E f d , i = K A , i 1 + s T A , i ( U R E F , i U C + U P S S , i ) ,
where the subscript i represents the generator i; K s is the amplification gain, T 5 is the high-pass filter time constant, usually set as 5 s or 10 s; in this paper, T 5 = 10   s . Two lead-lag compensators are used to eliminate the delay between excitation and electromagnetic torque. In practical applications, the two lead-lag compensators can compensate the low frequency and high frequency phases respectively. After each input signal passes through PSS, its output signal can provide the corresponding reference voltage for the excitation system, and its reference voltage serves as the reference modulation signal provided to the excitation system, so that the negative damping or weak damping caused by the fast excitation system can be compensated accordingly [16,17].
PSS plays an effective role mainly by adjusting K s and the inertia time constant of the lead-lag compensators, and the other parameters for a generator are usually unchanged. Therefore, it is necessary to tune parameters reasonably, otherwise the effect of suppressing low frequency oscillation will be counterproductive.
In general, a linear differential equation can be used to describe the dynamic behavior of a power system:
x ˙ = f ( x , u ) ,
where the column vector x is the state vector, u is the external input vector, and the derivative of the state variable to time is represented by x ˙ . When the variables of the system remain constant, the system is in equilibrium. The nonlinear system equation can then be linearized at the equilibrium point [18], where the equation is as follows:
{ Δ x ˙ = A Δ x + B Δ u , Δ y = C Δ x + D Δ u ,
where Δ x is the state vector, Δ y is the output vector, Δ u is the input vector, A is the state matrix, B is the control matrix, C is the output matrix, D is the feedforward matrix. For the power system, y is not a direct function of u , that is D = 0 . The stability of the system can be determined by judging the position of eigenvalues of the state matrix A on the complex plane [19].
Many studies show that the generator controller is effective in changing the real part of the oscillation mode; in the meantime, it has little influence on the imaginary part. As shown in Figure 2, the dotted line is equal to the damping ratio line. After a lot of engineering practice, the Ontario Electric Power Bureau of Canada proposed that ξ = 0.03 is the critical state in the normal operation of the power system. Thus,  ξ min is set as 0.03 in this paper [20]. In fact, the optimization of PSS is to move the eigenvalues of the state matrix of the system to the left part of the complex plane through continuous optimization, as far away from the virtual axis as possible. According to this rule, we can define the optimization objective function of the stabilizer and measure the dynamic characteristics of the system by the damping ratio of the system.
In the complex plane, each such conjugate complex root λ j = σ j + j w d j has the following relationship with the system damping ratio ξ j :
ξ j = σ j σ j 2 + w d j 2 .
Define the objective function f:
f = j = 1 n i = 1 k ( ξ min ξ i , j ) ,
where n represents the number of operating modes, k is the number of oscillation modes, ξ min is the preset minimum damping coefficient, ξ i , j represents the ith electromechanical oscillation mode damping coefficient in the jth operating mode. Combined with other constraints, the parameter coordination of PSS can be expressed as the following eigenvalue optimization form:
{ min f , s . t .   K i min K i K i m a x , T 1 i min T 1 i T 1 i m a x , T 3 i min T 3 i T 3 i m a x ,
where K i max = 50 , K i min = 0.1 , T 1 i m a x = T 3 i m a x = 1, K i , T 1 i and T 3 i are three variables to be optimized and the remaining parameters are given in advance.

3. PSS Parameter Optimization Based on SA-QUATRE Algorithm

The evolution formula of the quasi-affine transformation is analogous to the affine transformation in geometry. The affine transformation function ( f : X Y ) in geometry is as follows:
X M X + B .
The evolutionary structure of the QUATRE algorithm is to use X M X + M ¯ B . Supposing a population of N particles is searching in a D-dimensional space, and X in Formula (8) is used to represent the particle’s position. If the position of particle i is X i = [ x 1 , x 2 , , x D ] , and the population size is ps, the population position can be expressed as X = [ X 1 , X 2 , , X p s ] T . Then the particle position update formula is as follows:
{ X M X + M ¯ B , B = X g b e s t + c ( X r 1 X r 2 ) ,
where B represents the evolutionary guidance matrix, represents the bitwise multiplication of matrix elements. X r 1 and X r 2 are generated by randomly arranging the row vectors of the matrix X. Their difference is the difference matrix, which is used to represent the particle search radius. This search method helps to adapt to different search dimensions, c is the coefficient factor or step size of the differential matrix. If the ith particle obtains the best fitness value, it is recorded as X b e s t i . Then the global best coordinate matrix for each particle is shown below, and its size is ps * 1.
X b e s t = [ X b e s t i , X b e s t i , , X b e s t i ] T .
M is the cooperative search matrix. M ¯ represents the incidence matrix of M which is the core of the algorithm. This algorithm transforms all the ordinary individuals of the population into the optimal individuals and extends the global optimal individual search method to the whole population. Although this operation achieves individual equivalence, the complexity of the algorithm will become n times. Therefore, the collaborative architecture is introduced to reduce the search complexity of the algorithm, which is implemented by the M matrix.
The M matrix is transformed from its initialization matrix M i n i t . M i n i t is a D-dimensional Boolean matrix whose elements of the lower triangular matrix are all 1 and the rest elements are 0, as shown in Equation (11). Its stacking method is determined according to the population size ps, and the dimension D. There are three stacking methods: if ps = D, if ps = n * D, M i n i t is n times the previous case of vertical stack, if ps = n * D+ k and ps%D = k (% means remainder), the n * ps row of M i n i t is consistent with the second method, and the last k rows are the first k rows of the Boolean matrix in the first method. After the stacking is completed, the conversion is achieved through two consecutive operations. The first step is to randomly arrange the row elements of the matrix M i n i t and perform independent operations on each row element, The second step is to randomly arrange the row vectors of the matrix without changing the row elements [21]. So, M is compared with M i n i t , the only thing that changes is the position of 0 and 1 in each row and the position of each row vector.
Μ i n i t = [ 1 1 1 1 1 1 ] ~ [ 1 1 1 1 1 ] = M .
This cooperative structure can effectively solve the defects that individuals in the population cannot achieve information sharing and high coordination in the evolution process due to the existence of two states of global search and local search.
Although QUATRE can gain performance advantages and reduce time complexity by multiplying the population, it will also fall into a local optimum in the later stage of the algorithm due to the reduction of population diversity. The simulated annealing algorithm is a global search algorithm, which accepts the difference with a certain probability during the search process and jumps out of the local optimum. Based on selecting the appropriate temperature parameter T, this paper uses a simulated annealing algorithm to improve QUATRE, which can improve the deficiencies of QUATRE’s later local optimization, speed up the algorithm’s process, and better design the parameters of the power system stabilizer. Figure 3 shows the flow of the improved quasi-affine transformation evolution algorithm. The steps are summarized as follows:
Step1: Initialize the population and randomly initialize the particle position X = [ X 1 , X 2 , , X N ] T , limit the search range, set the difference matrix coefficient factor c, and set the initial position as the individual historical optimal solution X p b e s t of each particle as and the optimal global solution X g b e s t ;
Step2: Calculate the objective function f i t n e s s ( X i ) of each particle;
Step3: Find the optimal value f i t n e s s min ( X ) in the value of the initial particle objective function, and calculate the initial temperature T of the simulated annealing, T = f i t n e s s min ( X ) / log ( 0.2 ) ;
Step4: Calculate the sudden change probability at each temperature. The method is first to calculate the sudden change probability of each population corresponding to the best individual and then divide by the sum of the jump probability of all populations;
Step5: Randomly set and generate a probability. If the jump probability of a specific population is greater than or equal to this random value, update the X g b e s t ;
Step6: Evolve the parent population according to formula (10), update the position and generate the offspring population;
Step7: Calculate the optimal fitness value after each update, update the temperature of the simulated annealing algorithm. Output the optimal particles after meeting the conditions of the number of iterations or convergence accuracy.
This article applies the SA-QUATRE algorithm introduced above to the coordinated design of power system stabilizers. For a multi-machine system, all the PSS parameter tuning and optimization processes are carried out simultaneously. The parameters to be tuned for each PSS are T 1 and T 3 . Figure 4 shows the flow chart of optimizing PSS parameters by SA-QUATRE algorithm. Under the objective function and constraint conditions of Equation (6), the optimization steps are summarized as follows:
(1)
First, import the basic power flow and dynamic data of the grid and configure the number and position of power system stabilizers according to the electromechanical oscillation characteristics of the power system to be analyzed. In the power system, installing power system stabilizers for each generator is unrealistic. So, the residue method is used to choose the installation location of PSS in this article [22];
(2)
Set the parameters and operating mode of the power system stabilizer, linearize the system model and calculate the initial damping ratio in this mode through eigenvalues. If all are greater than ξmin, then end, otherwise continue;
(3)
Initialize the parameters of the SA-QUATRE algorithm and simulated annealing, and at the same time, generate a set of initial solutions X = { X 1 , X 2 , , X D } T ;
(4)
Assign each solution Xi to the variable in the PSS and calculate the damping coefficient of each electromechanical oscillation mode under this operating state;
(5)
Evaluate the group according to the objective function based on characteristic value;
(6)
Use the improved quasi-affine transformation evolutionary algorithm to continuously update and search, generate offspring populations, and update next-generation candidate solutions;
(7)
The operation ends when reaching the maximum number of iterations, otherwise, it returns to the (4) step and enters the next cycle;
(8)
Finally, obtain the global optimal parameter combination of PSS and the fitness value of the objective function. The specific process is as follows:

4. Simulation Example

4.1. Algorithm Robustness Analysis

In order to preliminarily verify that the improved QUATRE algorithm in this paper has better convergence accuracy and global search ability, eight internationally-used benchmark functions were used for testing. The function expressions are shown in Table 1. The number of populations of different algorithms is 25, the maximum number of iterations is 1000, and the dimension of the function is 10. Figure 5 shows the convergence curve of the function. It can be seen that the improved QUATRE algorithm can approach the theoretical optimal value at the beginning of the iteration, which shows its fast convergence.
In order to avoid the chance of experiments, each function ran 50 times independently. The average value and error result of each algorithm are shown in Table 2. It is clear from the data in the table that the optimal fitness value of this method is the closest to the target value and its standard deviation is the smallest compared with the other three algorithms. In general, the algorithm in this paper has the ability of global optimization, and the overall performance of the algorithm has a certain superiority.

4.2. IEEE 4-Machine 2-Area System

In order to verify the feasibility of the method proposed in this article in suppressing oscillations, the IEEE 4-machine 2-area system [23] is now used as an example, as shown in Figure 6. The generators G 1 and G 2 are installed in area 1. The generators G 3 and G 4 are installed in area 2, and a tie line connects the two areas. Because the two regions are connected by a contact line, the stability of the system itself is very weak. As long as there is a little disturbance outside the system, static instability easily occurs, and even leads to dynamic instability. Therefore, in practice, the system is always monitoring and protecting; it is sensitive and complex. After calculation and analysis, the PSS is configured in G 2 and G 4 , and the dimension of the parameters to be optimized is six, which include the gain ( K s [ 0.1 , 50 ] ) and phase compensation ( T 1 [ 0.01 , 1 ] , T 3 [ 0.01 , 1 ] ) of each PSS respectively.
The effectiveness of the proposed method is analyzed in the following two operation states of the system.

4.2.1. Small Disturbance

At the excitation reference voltages on the generators G 2 and G 4 , a square wave pulse disturbance signal with an amplitude of 0.05 appears at 1 s, and the duration is 0.1 s. In this disturbing operation mode, the electromechanical modes without power system stabilizers are shown in Table 3. It can be seen that the system has three oscillation modes, of which, the two modes with a frequency of about 1.1 Hz are regional modes, and the one with a frequency of 0.6409 Hz is the interval mode. The eigenvalues of mode 1 are in the right half of the complex plane, showing negative damping, which is the most crucial control target in the coordinated design of multi-machine PSS.
The simulation time of the model is 15 s, and the objective function is iterated 1000 times with different algorithms. Finally, the global optimal PSS parameter combination is searched for as shown in Table 4. In this operating mode, the system tie-line oscillation power after each method is adjusted for PSS as shown in Figure 7. The G1 power angle signal (take G4 as a reference) is shown in Figure 8.
After a small disturbance is applied, the power on the tie-line between the two regions of the system will increase and oscillate. If it is not suppressed, the system will be seriously destabilized. Under the PSS parameters under the system default parameter configuration, that is, when the unoptimized PSS is applied, the power oscillation is suppressed to a greater extent than when the PSS is not installed, but the time required is about 13 s. The PSS, under the coordinated design of SA-QUATRE, can effectively calm down in the early stage of oscillation and reflect the rapidity of this method. From Figure 7 and Figure 8, we can also see that, in comparison with other methods, before the oscillation is basically subsided, the overall amplitude of the PSS optimized by the method in this paper is more minor. In Figure 7a, you can see that under the traditional lead-lag compensators phase compensation setting, the power oscillation amplitude of the PSS appears as a local downdraft at 3~4 s, indicating that the overshoot is too large and the system damping is relatively tiny. This phenomenon is effectively alleviated after optimization. In summary, the PSS parameters and suppression effect designed by the algorithm in this paper are more superior.
Use the above methods to coordinate the design of PSS parameters and compare them with the case where PSS is not installed. Since the suppression effect and the speed changes of generators G1~G4 are similar, only G1 speed changes are attached in this article, as shown in Figure 9. In the absence of PSS, the oscillation of the system cannot be attenuated to equilibrium, but it can be effectively suppressed with the stabilizer. The proposed method is better for the synthesis of the inhibition time and the overall amplitude. Table 5 shows the characteristic roots and damping of the system after the PSS function under the coordinated design of various methods. From the table, we can see that the installation of PSS can generally enhance the damping of each oscillation mode of the system, especially the dominant mode, that is, the negative damping mode mentioned in Table 2. It is further verified that the system after the PSS optimized by the method in this paper is installed under the disturbance of this method can have higher stability.
In order to see the system changes more intuitively, Figure 10 shows the eigenvalue distribution of the electromechanical mode when the system is configured with the PSS optimized by the method in this paper and when the PSS is not installed. It can be clearly seen from the figure that in the state of no PSS, the electromechanical modes are close to the imaginary axis and even reach the right half of the complex plane, indicating that the system has been unstable or has reached a critical stable state. After configuring the stabilizer set according to the method in this paper, the eigenvalues are all moved to the desired area, which enhances the damping ratio of each mode.

4.2.2. Large Disturbance

This paper designs and verifies the transient performance of the power system stabilizer to further prove the applicability of the method proposed in this study. For the system shown in Figure 6, a large disturbance is set, in which occurs a three-phase short-circuit at 1 s that stops at 1/6 s. To coordinate the design of the PSS, the method described in this paper and the traditional phase compensation method were used, respectively, the comparison was made with the controller not the installed scenario. Figure 11 and Figure 12 depict the change in active power on the tie-line and the rotation speed of generator G1 in this mode of operation. Figure 13 depicts the G2 power angle signal change using the G4 power angle as a reference.
In the case of this large disturbance, the system has an amplitude-enhancing oscillation, and the PSS adjusted by the traditional phase compensation method is adopted. Although it can attenuate the oscillation to balance within a certain period of time, it is slightly inferior in speed and amplitude compared to the coordinated design method proposed in this article. The power system stabilizer optimized in this paper can effectively suppress it quickly, and the overall amplitude of the oscillation is smaller in the early stage of oscillation. Table 6 shows the various electromechanical modes of the system with the power system stabilizer designed by this method and without PSS. The table shows that in this disturbance operation state, a negative damping mode will be generated when the stabilizer is not configured, which threatens the system’s stable operation. After inputting the PSS tuned by the method proposed in this paper, not only is the negative damping of the dominant oscillation mode increased to the positive damping, but the damping of the other two modes is also enhanced at the same time.
The main objective of this paper is to maximize the damping ratio by moving the eigenvalues of the system to the left of the imaginary axis. From the characteristic distribution of electromechanical mode in Figure 14, it can be seen clearly and intuitively that the coordinated design of PSS using the method in this paper can also achieve the goal and effectively improve the system’s stability in a large disturbance state.

5. Conclusions

This paper proposes a multi-machine power system stabilizer parameter coordinated optimization design based on the SA-QUATRE algorithm to suppress power system low-frequency oscillation and enhance system stability. For the first time, the new algorithm is applied to power system problems, and the following conclusions are reached through experiments:
(1)
In the case of system oscillation instability, the installation of the power system stabilizer on the generator can effectively calm the oscillation, so that the system can be restored to normal operation in a specified time. However, the effect of PSS presented under different design methods is also different. Simulations show that using the method in this paper to design the parameters of PSS can convert the negative damping electromechanical mode that threatens the stability of the system to positive damping, and at the same time enhance the damping ratio of other oscillation modes, effectively improving the stability and robustness of the system, regardless of the operation mode of small or large disturbance;
(2)
For multiple power system stabilizers in the system, the traditional phase compensation method usually carries out parameter settings sequentially, which will inevitably cause the interaction effect between each stabilizer. This paper uses the performance target based on eigenvalues instead of the traditional phase compensation method to optimize the PSS parameters to achieve more effective suppression of the low-frequency oscillation of the power system;
(3)
The SA-QUATRE algorithm proposed in this paper can jump out of the local optimum and has better convergence and accuracy. This algorithm not only overcomes the defect of slow convergence of the particle cluster optimization algorithm, it also has stronger cooperation, reduces time complexity, and accepts inferior solutions with a certain probability to avoid falling into local minimum in the process of searching;
(4)
The SA-QUATRE algorithm has good applicability. The applicability of this method to the coordinated design of the damping controller of the flexible AC transmission system is still under study, and it is also worthy of further popularization and application to other fields;
(5)
In the future, PSS can be used to suppress the low-frequency oscillation caused by the addition of renewable energy to the traditional power grid and analyze the adaptability of PSS in new power systems.

Author Contributions

Formal analysis, C.Z.; Funding acquisition, J.H.; Investigation, J.L. and Y.K.; Methodology, J.H. and S.W.; Software, J.L.; Supervision, J.H. and C.Z.; Writing—original draft, Y.K.; Writing—review & editing, S.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China, grant number 51977039, grant number 61872095, grant number 61571139, grant number 61872128, in part by Fujian Institute of Technology Research Startup Fund Project, grant number GY-Z18060, in part by Fujian Science Fund for Distinguished Young Scholars, grant number 2020J06043.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xu, G.; Fang, Y.; Li, Z.; Wu, X.; Liu, F.; Li, W. Phenomenon and mechanism of interval tie-line power oscillation in ultra-low frequency oscillation. Power Syst. Autom. 2020, 44, 69–82. [Google Scholar]
  2. Hu, H.; Zhou, Y.; Li, X.; Lei, K. Low-frequency oscillation in electric railway depot: A comprehensive review. IEEE Trans. Power Electron. 2020, 36, 295–314. [Google Scholar] [CrossRef]
  3. Bai, F.; Wang, B. Influence of photovoltaic grid connection on low-frequency oscillation. J. Sol. Energy 2020, 41, 255–261. [Google Scholar]
  4. Obaid, Z.A.; Mejeed, R.A.; Al-Mashhadani, A. Investigating the Impact of using Modern Power System Stabilizers on Frequency Stability in Large Dynamic Multi-Machine Power System. In Proceedings of the 2020 55th International Universities Power Engineering Conference (UPEC), Turin, Italy, 1–4 September 2020; pp. 1–6. [Google Scholar]
  5. Mehrzad, A.; Darmiani, M.; Mousavi, Y.; Shafie-Khah, M.; Aghamohammadi, M. An Efficient Rapid Method for Generators Coherency Identification in Large Power Systems. IEEE Open Access J. Power Energy 2022, 9, 151–160. [Google Scholar] [CrossRef]
  6. Gude, M.K.; Salma, U. A novel approach of PSS optimal parameter tuning in a multi-area power system using hybrid butterfly optimization algorithm-particle swarm optimization. Int. J. Syst. Assur. Eng. Manag. 2022, 1–10. [Google Scholar] [CrossRef]
  7. Pan, X.; Zhang, L.; Zhang, W.; Xu, Y.; Bian, H. parameter coordination optimization method of multi-operation mode power system stabilizer based on optimization algorithm of moth-extinguishing. Power Grid Technol. 2020, 44, 3038–3046. [Google Scholar]
  8. Alzubi, I.; Al-Masri, H.M.K.; Abuelrub, A. Modified Particle Swarm Optimization Algorithms for Solving Economic Load Dispatch. In Proceedings of the 2022 3rd International Conference on Smart Grid and Renewable Energy (SGRE), Doha, Qatar, 20–22 March 2022; pp. 1–5. [Google Scholar]
  9. H, S.J.; Laly, M.J. Optimal Design of Power System Stabilizer for Damping Low Frequency Oscillations in a Multi-Machine Power System. In Proceedings of the 2020 International Conference on Power Electronics and Renewable Energy Applications (PEREA), Kannur, India, 27–28 November 2020; pp. 1–6. [Google Scholar]
  10. Rodrigues, F.; Molina, Y.; Silva, C.; Ñaupari, Z. Simultaneous tuning of the AVR and PSS parameters using particle swarm optimization with oscillating exponential decay. Int. J. Electr. Power Energy Syst. 2021, 133, 107215. [Google Scholar] [CrossRef]
  11. Devarapalli, R.; Bhattacharyya, B. A hybrid modified grey wolf optimization-sine cosine algorithm-based power system stabilizer parameter tuning in a multimachine power system. Optim. Control. Appl. Methods 2020, 41, 1143–1159. [Google Scholar] [CrossRef]
  12. Gurung, S.; Jurado, F.; Naetiladdanon, S.; Sangswang, A. Comparative analysis of probabilistic and deterministic approach to tune the power system stabilizers using the directional bat algorithm to improve system small-signal stability. Electr. Power Syst. Res. 2020, 181, 106176. [Google Scholar] [CrossRef]
  13. Meng, Z.; Pan, J.-S.; Xu, H. Quasi-Affine TRansformation Evolutionary (QUATRE) algorithm: A cooperative swarm based algorithm for global optimization. Knowl.-Based Syst. 2016, 97, 104–121. [Google Scholar] [CrossRef]
  14. Zhao, F.; Guo, C.; Si, J. Research on PSS parameter optimization based on bat algorithm. Control. Eng. 2019, 1, 21–25. [Google Scholar]
  15. Zhang, C.; Qiu, B. Parameters coordinated optimization design of power system stabilizer by social learning based particle swarm optimization algorithm. Adv. Technol. Electr. Eng. 2022, 41, 24–33. [Google Scholar]
  16. Li, Q.; Wu, C.; Chen, L.; Duan, R.; Huang, W. Parameter optimization of power system stabilizer for suppressing frequency oscillation. Power Syst. Autom. 2020, 44, 93–99. [Google Scholar]
  17. Saadatmand, M.; Mozafari, B.; Gharehpetian, G.B.; Soleymani, S. Optimal coordinated tuning of power system stabilizers and wide-area measurement-based fractional-order PID controller of large-scale PV farms for LFO damping in smart grids. Int. Trans. Electr. Energy Syst. 2021, 31, e12612. [Google Scholar] [CrossRef]
  18. Chen, J.; Jin, T.; Zhu, X.; Li, Z.; Zhang, K. Parameter optimization of power system stabilizer based on SOGWO. Power Syst. Prot. Control. 2020, 48, 159–164. [Google Scholar]
  19. Sreedivya, K.M.; Jeyanthy, P.A.; Devaraj, D. Improved design of interval type-2 fuzzy based wide area power system stabilizer for inter-area oscillation damping. Microprocess. Microsyst. 2021, 83, 103957. [Google Scholar] [CrossRef]
  20. Liu, Q. Power System Stability and Generator Excitation Control; China Electric Power Press: Beijing, China, 2007. [Google Scholar]
  21. Zhenyu, M. Research on QUasi-Affine Transformation Evolutionary Algorithm with Cooperative Structure; Harbin Institute of Technology: Harbin, China, 2018. [Google Scholar]
  22. Du, W.; Wang, H. The Theory and Method of Power System Low-Frequency Power Oscillation Mode Analysis; Science Press: Beijing, China, 2017. [Google Scholar]
  23. Chen, J.; Jin, T.; Mohamed, M.A.; Wang, M. An adaptive TLS-ESPRIT algorithm based on an S-G filter for analysis of low frequency oscillation in wide area measurement systems. IEEE Access 2019, 7, 47644–47654. [Google Scholar] [CrossRef]
Figure 1. Excitation system with traditional lead-lag compensation type PSS.
Figure 1. Excitation system with traditional lead-lag compensation type PSS.
Electronics 11 02785 g001
Figure 2. The left region of s-plane.
Figure 2. The left region of s-plane.
Electronics 11 02785 g002
Figure 3. Flow chart of SA-QUATRE algorithm.
Figure 3. Flow chart of SA-QUATRE algorithm.
Electronics 11 02785 g003
Figure 4. Flow chart of optimizing PSS parameters by SA-QUATRE algorithm.
Figure 4. Flow chart of optimizing PSS parameters by SA-QUATRE algorithm.
Electronics 11 02785 g004
Figure 5. Iterative optimization of f1 test function.
Figure 5. Iterative optimization of f1 test function.
Electronics 11 02785 g005
Figure 6. IEEE 4−machine 2−area system, structural diagram and Simulink model.
Figure 6. IEEE 4−machine 2−area system, structural diagram and Simulink model.
Electronics 11 02785 g006
Figure 7. Active power output of tie-line: (a) Comparison between the method in this paper and the unoptimized PSS; (b) Comparison of multiple algorithms to optimize PSS.
Figure 7. Active power output of tie-line: (a) Comparison between the method in this paper and the unoptimized PSS; (b) Comparison of multiple algorithms to optimize PSS.
Electronics 11 02785 g007
Figure 8. Power angle signal: (a) Comparison between the method in this paper and the unoptimized PSS; (b) Comparison of multiple algorithms to optimize PSS.
Figure 8. Power angle signal: (a) Comparison between the method in this paper and the unoptimized PSS; (b) Comparison of multiple algorithms to optimize PSS.
Electronics 11 02785 g008
Figure 9. IEEE 4-machine 2-area system.
Figure 9. IEEE 4-machine 2-area system.
Electronics 11 02785 g009
Figure 10. Distribution of eigenvalues.
Figure 10. Distribution of eigenvalues.
Electronics 11 02785 g010
Figure 11. Active power output of tie-line.
Figure 11. Active power output of tie-line.
Electronics 11 02785 g011
Figure 12. Speed change.
Figure 12. Speed change.
Electronics 11 02785 g012
Figure 13. Power angle signal.
Figure 13. Power angle signal.
Electronics 11 02785 g013
Figure 14. Distribution of eigenvalues.
Figure 14. Distribution of eigenvalues.
Electronics 11 02785 g014
Table 1. Eight general benchmark functions.
Table 1. Eight general benchmark functions.
Functional ExpressionSearch AreaOptimum Value
f 1 = i = 1 D x i 2 [−100,100]0
f 2 = x i 2 + 10 6 i = 2 D x i 2 [−100,100]01
f 3 = 10 6 x i 2 + i = 2 D x i 2 [−100,100]0
f 4 = i = 1 D ( x i 2 10 cos ( 2 π x i ) + 10 ) [−100,100]0
f 5 = i = 1 D 1 ( 100 ( x i 2 x i + 1 ) 2 + ( x i 1 ) 2 ) [−100,100]0
f 6 ( x ) = 1 4000 i = 1 D x i 2 Π i = 1 D cos ( x i i ) + 1 [−100,100]0
f 7 = 20 exp ( 0.2 1 D i = 1 D x i 2 ) exp ( 1 D i = 1 D cos ( 2 π x i ) ) + 20 + e [−100,100]0
f 8 = ( 1 D 1 i = 1 D 1 ( x i + x i sin 2 ( 50 x i 0.2 ) ) ) 2 [−100,100]0
Table 2. Results of eight general benchmark functions.
Table 2. Results of eight general benchmark functions.
FunctionSA-QUATREPSOQUATRESA
Mean ValueStandard DeviationMean ValueStandard DeviationMean ValueStandard DeviationMean ValueStandard Deviation
f1000.00120.00321.08 × 10−263.10 × 10−251.69 × 104262.974
f20.29611.20473.82 × 1088.52 × 1080.34301.33702.90 × 10191.53 × 1019
f33.09 × 10−73.05 × 10−7854.7821483.983.403 × 10−75.441 × 10−76.51 × 1063.80 × 106
f40.04840.01951.00 × 1064.01 × 1050.63400.54771.17 × 1071.02 × 106
f52.74243.763816.07322.9833.74414.76381516.1651.090
f60.15550.08862.88083.19240.20100.17732.881 × 10374.165
f720.0610.077720.3230.087520.4540.092020.5560.0880
f80.77532.672548.42622.2620.88903.55211.35 × 1073.41 × 106
Table 3. Electromechanical modes without PSS.
Table 3. Electromechanical modes without PSS.
Mode 1Mode 2Mode 3
characteristic root0.1029 ± 4.0251i−0.6743 ± 7.2483i−0.6804 ± 7.0347i
f/Hz0.64091.15421.1202
ξ −0.02560.09260.0927
Table 4. Electromechanical modes without PSS.
Table 4. Electromechanical modes without PSS.
MethodG2G4
KST1T3KST1T3
PSO49.940.570.9949.740.390.99
SA46.050.810.9642.261.000.58
QUATRE42.670.650.6638.600.750.39
This paper49.660.960.69349.1460.5160.98
Table 5. Characteristic root and damping ratio of the system after designing PSS by different methods.
Table 5. Characteristic root and damping ratio of the system after designing PSS by different methods.
TypeCharacteristic RootDamping RatioMinimum Damping Ratio
this paper−1.6173 ± 5.2151i0.29620.1933
−1.3001 ± 4.9332i0.2548
−0.7476 ± 3.7957i0.1933
QUATRE−1.6761 ± 5.3423i0.29940.1786
−1.4534 ± 4.5378i0.3050
−0.6497 ± 3.5789i0.1786
SA−1.8203 ± 6.3656i0.27490.0671
−0.2659 ± 3.9544i0.0671
−0.5773 ± 0.8160i0.1932
PSO−1.9048 ± 5.6322i0.32040.1350
−1.8905 ± 4.7963i0.3667
−0.5381 ± 3.9493i0.1350
Table 6. System characteristic root and damping without and with PSS installed.
Table 6. System characteristic root and damping without and with PSS installed.
ModeNot Installed PSSDamping RatioMethod of This ArticleDamping Ratio
Characteristic RootCharacteristic Root
10.1051 ± 4.0277i−0.0261−0.6497 ± 3.5789i0.1786
2−0.6683 ± 7.2631i0.0916−1.6173 ± 5.2152i0.2962
3−0.6756 ± 7.0499i0.0954−1.3001 ± 4.9333i0.2548
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Huang, J.; Liu, J.; Zhang, C.; Kuang, Y.; Weng, S. Application of Improved Quasi-Affine Transformation Evolutionary Algorithm in Power System Stabilizer Optimization. Electronics 2022, 11, 2785. https://doi.org/10.3390/electronics11172785

AMA Style

Huang J, Liu J, Zhang C, Kuang Y, Weng S. Application of Improved Quasi-Affine Transformation Evolutionary Algorithm in Power System Stabilizer Optimization. Electronics. 2022; 11(17):2785. https://doi.org/10.3390/electronics11172785

Chicago/Turabian Style

Huang, Jing, Jiajing Liu, Cheng Zhang, Yu Kuang, and Shaowei Weng. 2022. "Application of Improved Quasi-Affine Transformation Evolutionary Algorithm in Power System Stabilizer Optimization" Electronics 11, no. 17: 2785. https://doi.org/10.3390/electronics11172785

APA Style

Huang, J., Liu, J., Zhang, C., Kuang, Y., & Weng, S. (2022). Application of Improved Quasi-Affine Transformation Evolutionary Algorithm in Power System Stabilizer Optimization. Electronics, 11(17), 2785. https://doi.org/10.3390/electronics11172785

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop